#include <ossimBilinearProjection.h>

Inheritance diagram for ossimBilinearProjection:

Public Member Functions
	ossimBilinearProjection ()

	ossimBilinearProjection (const ossimBilinearProjection &rhs)

	ossimBilinearProjection (const ossimDpt &ul, const ossimDpt &ur, const ossimDpt &lr, const ossimDpt &ll, const ossimGpt &ulg, const ossimGpt &urg, const ossimGpt &lrg, const ossimGpt &llg)

virtual ossimObject *	dup () const

virtual ossimGpt	origin () const

virtual void	worldToLineSample (const ossimGpt &worldPoint, ossimDpt &lineSampPt) const

virtual void	lineSampleToWorld (const ossimDpt &lineSampPt, ossimGpt &worldPt) const

virtual void	lineSampleHeightToWorld (const ossimDpt &lineSampPt, const double &heightAboveEllipsoid, ossimGpt &worldPt) const

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)

virtual bool	operator== (const ossimProjection &projection) const

virtual ossimDpt	getMetersPerPixel () const

virtual std::ostream &	print (std::ostream &out) const
	Outputs theErrorStatus as an ossimErrorCode and an ossimString. More...

virtual ossim_float64	setTiePoints (const std::vector< ossimDpt > &lsPt, const std::vector< ossimGpt > &geoPt)

virtual bool	setupOptimizer (const ossimString &setup)
	setupFromString() Derived classes should implement as needed. Initialize parameters needed for optimizeFit and degreesOfFreedom More...

virtual bool	useForward () const

virtual ossim_uint32	degreesOfFreedom () const
	better go from image to ground, also means that errors variance are in squared meters More...

virtual double	optimizeFit (const ossimTieGptSet &tieSet, double *targetVariance=0)

virtual bool	isAffectedByElevation () const
	Implementation of pure virtual ossimProjection::isAffectedByElevation method. More...

void	getTiePoints (std::vector< ossimDpt > &lsPt, std::vector< ossimGpt > &geoPt) const
	Access method for tie point information. More...

Public Member Functions inherited from ossimProjection
	ossimProjection ()

virtual	~ossimProjection ()

virtual ossimDpt	forward (const ossimGpt &wp) const

virtual ossimGpt	inverse (const ossimDpt &pp) const

virtual void	getRoundTripError (const ossimDpt &imagePoint, ossimDpt &errorResult) const

virtual void	getRoundTripError (const ossimGpt &groundPoint, ossimDpt &errorResult) const

virtual void	getGroundClipPoints (ossimGeoPolygon &gpts) const

virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const

virtual bool	operator!= (const ossimProjection &projection) const

Public Member Functions inherited from ossimObject
	ossimObject ()

virtual	~ossimObject ()

virtual ossimString	getShortName () const

virtual ossimString	getLongName () const

virtual ossimString	getDescription () const

virtual ossimString	getClassName () const

virtual RTTItypeid	getType () const

virtual bool	canCastTo (ossimObject *obj) const

virtual bool	canCastTo (const RTTItypeid &id) const

virtual bool	canCastTo (const ossimString &parentClassName) const

virtual void	accept (ossimVisitor &visitor)

Public Member Functions inherited from ossimReferenced
	ossimReferenced ()

	ossimReferenced (const ossimReferenced &)

ossimReferenced &	operator= (const ossimReferenced &)

void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...

void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

int	referenceCount () const

Public Member Functions inherited from ossimErrorStatusInterface
	ossimErrorStatusInterface ()

virtual	~ossimErrorStatusInterface ()

virtual ossimErrorCode	getErrorStatus () const

virtual ossimString	getErrorStatusString () const

virtual void	setErrorStatus (ossimErrorCode error_status) const

virtual void	setErrorStatus () const

virtual void	clearErrorStatus () const

bool	hasError () const

Public Member Functions inherited from ossimOptimizableProjection
	ossimOptimizableProjection ()

	ossimOptimizableProjection (const ossimOptimizableProjection &source)

virtual	~ossimOptimizableProjection ()

virtual ossimOptimizableProjection &	operator= (const ossimOptimizableProjection &source)

virtual bool	needsInitialState () const
	needsInitialState() More...

Protected Member Functions
virtual	~ossimBilinearProjection ()

void	initializeBilinear ()

bool	dPtsHaveNan () const
	Checks theLineSamplePt for nans. More...

bool	gPtsHaveNan () const
	Checks theGeographicPt for nans. More...

ossimDpt	midLineSamplePt () const

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes
std::vector< ossimDpt >	theLineSamplePt

std::vector< ossimGpt >	theGeographicPt

ossimLeastSquaresBilin	theLatFit

ossimLeastSquaresBilin	theLonFit

ossimLeastSquaresBilin	theXFit

ossimLeastSquaresBilin	theYFit

bool	theInterpolationPointsHaveNanFlag

bool	theInverseSupportedFlag

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Detailed Description

Definition at line 19 of file ossimBilinearProjection.h.

Constructor & Destructor Documentation

◆ ossimBilinearProjection() [1/3]

ossimBilinearProjection::ossimBilinearProjection ( )

Definition at line 37 of file ossimBilinearProjection.cpp.

Referenced by dup().

    :
       ossimOptimizableProjection(),
       theLineSamplePt(0),
       theGeographicPt(0),
       theLatFit(),
       theLonFit()
 {
 #ifdef OSSIM_ID_ENABLED
    if (traceDebug()) ossimNotify(ossimNotifyLevel_DEBUG) << OSSIM_ID << endl;
 #endif    
 }

◆ ossimBilinearProjection() [2/3]

ossimBilinearProjection::ossimBilinearProjection ( const ossimBilinearProjection & rhs )

Definition at line 50 of file ossimBilinearProjection.cpp.

    :
       ossimOptimizableProjection(rhs),
       theLineSamplePt(rhs.theLineSamplePt),
       theGeographicPt(rhs.theGeographicPt),
       theLonFit(rhs.theLonFit)
 {
 }

◆ ossimBilinearProjection() [3/3]

ossimBilinearProjection::ossimBilinearProjection	(	const ossimDpt &	ul,
		const ossimDpt &	ur,
		const ossimDpt &	lr,
		const ossimDpt &	ll,
		const ossimGpt &	ulg,
		const ossimGpt &	urg,
		const ossimGpt &	lrg,
		const ossimGpt &	llg
	)

Definition at line 59 of file ossimBilinearProjection.cpp.

References initializeBilinear(), theGeographicPt, and theLineSamplePt.

    :
       ossimOptimizableProjection(),
       theLineSamplePt(4),
       theGeographicPt(4),
       theLatFit(),
       theLonFit()
 {
    theLineSamplePt[0] = ul;
    theLineSamplePt[1] = ur;
    theLineSamplePt[2] = lr;
    theLineSamplePt[3] = ll;
 
    theGeographicPt[0] = ulg;
    theGeographicPt[1] = urg;
    theGeographicPt[2] = lrg;
    theGeographicPt[3] = llg;
    
    initializeBilinear();
 }

◆ ~ossimBilinearProjection()

ossimBilinearProjection::~ossimBilinearProjection ( )

protectedvirtual

Definition at line 87 of file ossimBilinearProjection.cpp.

88 {

89 }

Member Function Documentation

◆ degreesOfFreedom()

ossim_uint32 ossimBilinearProjection::degreesOfFreedom ( ) const

virtual

better go from image to ground, also means that errors variance are in squared meters

Implements ossimOptimizableProjection.

Definition at line 534 of file ossimBilinearProjection.cpp.

 {
    return 2*4; //height not used
 }

◆ dPtsHaveNan()

bool ossimBilinearProjection::dPtsHaveNan ( ) const

protected

Checks theLineSamplePt for nans.

Returns: true if any point has a nan.

Note: If theLineSamplePt size is 0 this returns false.

Definition at line 391 of file ossimBilinearProjection.cpp.

References theLineSamplePt.

Referenced by initializeBilinear(), and midLineSamplePt().

 {
    if (theLineSamplePt.size() == 0)
    {
       return false;
    }
 
    vector<ossimDpt>::const_iterator i = theLineSamplePt.begin();
    while (i != theLineSamplePt.end())
    {
       if ( (*i).hasNans() )
       {
          return true;
       }
       ++i;
    }
    return false;
 }

◆ dup()

ossimObject * ossimBilinearProjection::dup ( ) const

virtual

Implements ossimProjection.

Definition at line 91 of file ossimBilinearProjection.cpp.

References ossimBilinearProjection().

 {
    return new ossimBilinearProjection(*this);
 }

◆ getMetersPerPixel()

ossimDpt ossimBilinearProjection::getMetersPerPixel ( ) const

virtual

ACCESS METHODS:

Implements ossimProjection.

Definition at line 313 of file ossimBilinearProjection.cpp.

References lineSampleToWorld(), ossimEcefVector::magnitude(), midLineSamplePt(), and ossimDpt::x.

 {
    ossimGpt centerG;
    ossimGpt rightG;
    ossimGpt topG;
 
    ossimDpt midPoint = midLineSamplePt();
                      
    lineSampleToWorld(midPoint, centerG);
    lineSampleToWorld(midPoint+ossimDpt(1,0), rightG);
    lineSampleToWorld(midPoint+ossimDpt(0,-1), topG);
    
    ossimEcefPoint centerP = centerG;
    ossimEcefPoint rightP  = rightG;
    ossimEcefPoint topP    = topG;
 
    ossimEcefVector horizontal = rightP-centerP;
    ossimEcefVector vertical   = topP-centerP;
 
    ossimDpt result(horizontal.magnitude(),
                    vertical.magnitude());
 
    result.x = (result.x + result.y)/2.0;
    result.y = result.x;
    
    return result;
 }

◆ getTiePoints()

void ossimBilinearProjection::getTiePoints	(	std::vector< ossimDpt > &	lsPt,
		std::vector< ossimGpt > &	geoPt
	)		const

Access method for tie point information.

Definition at line 552 of file ossimBilinearProjection.cpp.

References theGeographicPt, and theLineSamplePt.

 {
    lsPt  = theLineSamplePt;
    geoPt = theGeographicPt;
 }

◆ gPtsHaveNan()

bool ossimBilinearProjection::gPtsHaveNan ( ) const

protected

Checks theGeographicPt for nans.

Returns: true if any point has a nan.

Note: If theLineSamplePt size is 0 this returns false.

Definition at line 410 of file ossimBilinearProjection.cpp.

References theGeographicPt.

Referenced by initializeBilinear(), and origin().

 {
    //---
    // NOTE:  This method ignores nans in the height field.
    //---
    
    if (theGeographicPt.size() == 0)
    {
       return false;
    }
 
    vector<ossimGpt>::const_iterator i = theGeographicPt.begin();
    while (i != theGeographicPt.end())
    {
       if ( (*i).isLatNan() || (*i).isLonNan() )
       {
          return true;
       }
       ++i;
    }
    return false;
 }

◆ initializeBilinear()

void ossimBilinearProjection::initializeBilinear ( )

protected

Definition at line 341 of file ossimBilinearProjection.cpp.

References ossimLeastSquaresBilin::addSample(), ossimLeastSquaresBilin::clear(), dPtsHaveNan(), ossimProjection::getRoundTripError(), gPtsHaveNan(), ossimDpt::length(), ossimLeastSquaresBilin::solveLS(), theGeographicPt, theInterpolationPointsHaveNanFlag, theInverseSupportedFlag, theLatFit, theLineSamplePt, theLonFit, theXFit, theYFit, x, and y.

Referenced by ossimBilinearProjection(), and setTiePoints().

 {
    theInterpolationPointsHaveNanFlag = dPtsHaveNan()||gPtsHaveNan();
    theInverseSupportedFlag = true;
    if(!theInterpolationPointsHaveNanFlag)
    {
       theLatFit.clear();
       theLonFit.clear();
       theXFit.clear();
       theYFit.clear();
       
       const ossim_uint32 SIZE = (ossim_uint32)theLineSamplePt.size();
       if (SIZE != theGeographicPt.size())
       {
          return;
       }
       
       for (ossim_uint32 i = 0; i < SIZE; ++i)
       {
          theLatFit.addSample(theLineSamplePt[i].x,
                              theLineSamplePt[i].y,
                              theGeographicPt[i].latd());
          
          theLonFit.addSample(theLineSamplePt[i].x,
                              theLineSamplePt[i].y,
                              theGeographicPt[i].lond());
          
          theXFit.addSample(theGeographicPt[i].lond(),
                            theGeographicPt[i].latd(),
                            theLineSamplePt[i].x);
          theYFit.addSample(theGeographicPt[i].lond(),
                            theGeographicPt[i].latd(),
                            theLineSamplePt[i].y);
          
       }
       
       theLatFit.solveLS();
       theLonFit.solveLS();
       theXFit.solveLS();
       theYFit.solveLS();
       ossimDpt errorResult;
       getRoundTripError(theLineSamplePt[0],
                         errorResult);
       if(errorResult.length() > 1)
       {
          theInverseSupportedFlag = false;
       }
    }
 }

◆ isAffectedByElevation()

virtual bool ossimBilinearProjection::isAffectedByElevation ( ) const

inlinevirtual

Implementation of pure virtual ossimProjection::isAffectedByElevation method.

Returns: false.

Implements ossimProjection.

Definition at line 90 of file ossimBilinearProjection.h.

90 { return false; }

◆ lineSampleHeightToWorld()

void ossimBilinearProjection::lineSampleHeightToWorld	(	const ossimDpt &	lineSampPt,
		const double &	heightAboveEllipsoid,
		ossimGpt &	worldPt
	)		const

virtual

METHOD: lineSampleHeightToWorld This is the pure virtual that projects the image point to the given elevation above ellipsoid, thereby bypassing reference to a DEM. Useful for projections that are sensitive to elevation (such as sensor models).

Implements ossimProjection.

Definition at line 155 of file ossimBilinearProjection.cpp.

References ossimGpt::datum(), ossimGpt::hgt, ossim::isnan(), ossimGpt::lat, ossimGpt::lon, ossimLeastSquaresBilin::lsFitValue(), ossimGpt::makeNan(), theGeographicPt, theInterpolationPointsHaveNanFlag, theLatFit, theLonFit, ossimDpt::x, and ossimDpt::y.

Referenced by lineSampleToWorld().

 {
    worldPt.makeNan();
 
    if (theInterpolationPointsHaveNanFlag)
    {
       return;
    }
 
    
    worldPt.lat = theLatFit.lsFitValue(lineSampPt.x, lineSampPt.y);
    worldPt.lon = theLonFit.lsFitValue(lineSampPt.x, lineSampPt.y);
    if (ossim::isnan(heightAboveEllipsoid) == false)
    {
       worldPt.hgt = heightAboveEllipsoid;
    }
    if (theGeographicPt.size())
    {
       worldPt.datum(theGeographicPt[0].datum());
    }
 }   

◆ lineSampleToWorld()

void ossimBilinearProjection::lineSampleToWorld	(	const ossimDpt &	lineSampPt,
		ossimGpt &	worldPt
	)		const

virtual

METHOD: lineSampleToWorld() Performs the inverse projection from line, sample to ground (world):

Implements ossimProjection.

Definition at line 146 of file ossimBilinearProjection.cpp.

References lineSampleHeightToWorld(), and ossim::nan().

Referenced by getMetersPerPixel().

 {
    lineSampleHeightToWorld(lineSampPt,
                            ossim::nan(),
                            worldPt);
    
 }

◆ loadState()

bool ossimBilinearProjection::loadState	(	const ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)

virtual

Method to the load (recreate) the state of the object from a keyword list. Return true if ok or false on error.

Reimplemented from ossimProjection.

Definition at line 228 of file ossimBilinearProjection.cpp.

References ossimString::empty(), ossimKeywordlist::find(), ossimGpt::height(), ossimGpt::isHgtNan(), ossimProjection::loadState(), ossimKeywordlist::numberOf(), theGeographicPt, theLineSamplePt, ossimDpt::toPoint(), ossimGpt::toPoint(), ossimString::toString(), and ossim::toVector().

 {
    // Load the base class.
    ossimProjection::loadState(kwl, prefix);
 
    // Start with clear lists.
    theLineSamplePt.clear();
    theGeographicPt.clear();
 
    ossimString imagePoints  = kwl.find(prefix, "image_points");
    ossimString groundPoints = kwl.find(prefix, "ground_points");
    
    if(!imagePoints.empty()&&!groundPoints.empty())
    {
       ossim::toVector(theLineSamplePt, imagePoints);
       ossim::toVector(theGeographicPt, groundPoints);
    }
    else 
    {
       //---
       // Get the number of points.
       // If 0 or gpt size not equal to dpt size get out.
       //---
       const ossim_uint32 SIZE = kwl.numberOf(prefix, "gpt");
       if ( (SIZE == 0) || (SIZE != kwl.numberOf(prefix, "dpt")) )
       {
          return false;
       }
       
       for (ossim_uint32 i = 0; i < SIZE; ++i)
       {
          const char* lookup;
          ossimString index_string = ossimString::toString(i);
          
          // Get the geographic point.
          ossimString kw = "gpt";
          kw += index_string;
          lookup = kwl.find(prefix, kw);
          if (lookup)
          {
             ossimGpt gp;
             gp.toPoint(std::string(lookup));
             
             //---
             // Allow for "nan" height values by substituting with 0.0 so the
             // hasNans() will work.  "nan"s will get placed in the point if the
             // user doesn't have the elevation manager preferences set up
             // correctly.
             //---
             if (gp.isHgtNan())
             {
                gp.height(0.0);
             }
             theGeographicPt.push_back(gp);
          }
          
          // Get the line sample point.
          kw = "dpt";
          kw += index_string;
          lookup = kwl.find(prefix, kw);
          if (lookup)
          {
             ossimDpt dp;
             dp.toPoint(std::string(lookup));
             theLineSamplePt.push_back(dp);
          }
       }
    }
 
    if (traceDebug())
    {
       print(ossimNotify(ossimNotifyLevel_DEBUG));
    }
    
    initializeBilinear();
 
    return true;
 }

◆ midLineSamplePt()

ossimDpt ossimBilinearProjection::midLineSamplePt ( ) const

protected

Returns: The mid point of theLineSamplePts.

Definition at line 433 of file ossimBilinearProjection.cpp.

References dPtsHaveNan(), ossimDpt::makeNan(), theLineSamplePt, ossimDpt::x, x, ossimDpt::y, and y.

Referenced by getMetersPerPixel().

 {
    ossimDpt result;
    
    if ( (theLineSamplePt.size() == 0) || dPtsHaveNan())
    {
       result.makeNan();
       return result;
    }
 
    double x = 0.0;
    double y = 0.0;
    vector<ossimDpt>::const_iterator i = theLineSamplePt.begin();
    while (i != theLineSamplePt.end())
    {
       x += (*i).x;
       y += (*i).y;
       ++i;
    }
    const double SIZE = theLineSamplePt.size();
    result.x = x / SIZE;
    result.y = y / SIZE;
 
    return result;
 }

◆ operator==()

bool ossimBilinearProjection::operator== ( const ossimProjection & projection ) const

virtual

OPERATOR: == Compares this instance with arg projection.

Implements ossimProjection.

Definition at line 308 of file ossimBilinearProjection.cpp.

 {
    return false;
 }

◆ optimizeFit()

double ossimBilinearProjection::optimizeFit	(	const ossimTieGptSet &	tieSet,
		double *	targetVariance = `0`
	)

virtual

optimizeFit() adjusts projection to GLOBALLY fit best to ALL tie points (no outlier rejection) can select an optional threshold for accuracy

RETURNS : resulting image error variance of fitted model, or negative value for error

IMPORTANT : the unit for *targetVariance and the returned value is : meters^2 when useForward() is true OR pixel^2 if useForward() is false

Implements ossimOptimizableProjection.

Definition at line 540 of file ossimBilinearProjection.cpp.

References ossimTieGptSet::getSlaveMasterPoints(), and setTiePoints().

 {
    //NOTE : IGNORE targetVariance
    std::vector<ossimDpt> imagePoints;
    std::vector<ossimGpt> groundPoints;
    tieSet.getSlaveMasterPoints(imagePoints, groundPoints);
    return setTiePoints(imagePoints, groundPoints); //variance in meters
 }

◆ origin()

ossimGpt ossimBilinearProjection::origin ( ) const

virtual

METHOD: origin() Returns projection's ground point origin. That is the GP corresponding to line=0, sample=0.

Implements ossimProjection.

Definition at line 96 of file ossimBilinearProjection.cpp.

References ossimGpt::datum(), gPtsHaveNan(), ossimGpt::height(), ossimGpt::latd(), ossimGpt::lond(), ossimGpt::makeNan(), and theGeographicPt.

 {
    ossimGpt result;
    result.makeNan();
    if ( (theGeographicPt.size() == 0) || gPtsHaveNan() )
    {
       return result;
    }
 
    double lat        = 0.0;
    double lon        = 0.0;
    const double SIZE = theGeographicPt.size();
    
    vector<ossimGpt>::const_iterator i = theGeographicPt.begin();
    while (i != theGeographicPt.end())
    {
       lat += (*i).latd();
       lon += (*i).lond();
       ++i;
    }
 
    result.latd(lat/SIZE);
    result.lond(lon/SIZE);
    result.height(0.0);
    result.datum(theGeographicPt[0].datum());
    
    return result;
 }

◆ print()

std::ostream & ossimBilinearProjection::print ( std::ostream & out ) const

virtual

Outputs theErrorStatus as an ossimErrorCode and an ossimString.

Returns: std::ostream&

Note: Derived classes should only have to implement a virtual print, not an operator<< funtion as it's implemented here.

Reimplemented from ossimProjection.

Definition at line 459 of file ossimBilinearProjection.cpp.

References ossimNotify(), ossimNotifyLevel_INFO, ossimProjection::print(), theGeographicPt, and theLineSamplePt.

 {
    ossimNotify(ossimNotifyLevel_INFO)
       << "ossimBilinearProjection::print\n";
 
    ossim_uint32 index = 0;
    vector<ossimDpt>::const_iterator di = theLineSamplePt.begin();
    while (di != theLineSamplePt.end())
    {
       ossimNotify(ossimNotifyLevel_INFO)
          << "theLineSamplePt[" << index << "]:  "
          << (*di) << endl;
       ++di;
       ++index;
    }
 
    index = 0;
    vector<ossimGpt>::const_iterator gi = theGeographicPt.begin();
    while (gi != theGeographicPt.end())
    {
       ossimNotify(ossimNotifyLevel_INFO)
          << "theGeographicPt[" << index << "]:  "
          << (*gi) << endl;
       ++gi;
       ++index;
    }
 
    return ossimProjection::print(out);
 }

◆ saveState()

bool ossimBilinearProjection::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)		const

virtual

METHODS: saveState, loadState Fulfills ossimObject base-class pure virtuals.

Reimplemented from ossimProjection.

Definition at line 180 of file ossimBilinearProjection.cpp.

References ossimKeywordlist::add(), ossimString::c_str(), os2, ossimProjection::saveState(), theGeographicPt, theLineSamplePt, ossimString::toString(), and ossim::toStringList().

 {
    if (theLineSamplePt.size() != theGeographicPt.size())
    {
       // Should never happen.
       return false;
    }
 
    ossimProjection::saveState(kwl, prefix);
 
    ossimString imagePoints;
    ossimString groundPoints;
    ossim::toStringList(imagePoints, theLineSamplePt);
    ossim::toStringList(groundPoints, theGeographicPt);
    kwl.add(prefix, 
            "image_points",
            imagePoints,
            true);
    kwl.add(prefix, 
            "ground_points",
            groundPoints,
            true);
 #if 0
    const ossim_uint32 SIZE = (ossim_uint32)theLineSamplePt.size();
 
    for (ossim_uint32 i = 0; i < SIZE; ++i)
    {
       ossimString index_string = ossimString::toString(i);
 
       // Add the geographic point.
       ossimString kw = "gpt";
       kw += index_string;
       ostringstream os1;
       os1 << theGeographicPt[i];
       kwl.add(prefix, kw, os1.str().c_str());
 
       // Add the sample line..
       kw = "dpt";
       kw += index_string;
       ostringstream os2;
       os2 << theLineSamplePt[i];
       kwl.add(prefix, kw, os2.str().c_str());
    }
 #endif
    return true;
 }

◆ setTiePoints()

ossim_float64 ossimBilinearProjection::setTiePoints	(	const std::vector< ossimDpt > &	lsPt,
		const std::vector< ossimGpt > &	geoPt
	)

virtual

METHOD: setTiePoints add any number of tie points and calculate bilinear fit returns the ground error variance (=RMS^2), in meters^2

Definition at line 489 of file ossimBilinearProjection.cpp.

References ossimGpt::hgt, initializeBilinear(), ossimGpt::lat, ossimGpt::lon, ossimLeastSquaresBilin::lsFitValue(), ossimNotify(), ossimNotifyLevel_INFO, theGeographicPt, theLatFit, theLineSamplePt, and theLonFit.

Referenced by ossim_hdf5::getBilinearProjection(), and optimizeFit().

 {
    if (lsPt.size() != geoPt.size())
    {
       ossimNotify(ossimNotifyLevel_INFO)
          << "mismatch in image and ground point number" << endl;
       return -1.0;
    }
    if (lsPt.size() < 4)
    {
       ossimNotify(ossimNotifyLevel_INFO)
          << "not enough tie points - need at least 4" << endl;
       return -1.0;
    }
 
    theLineSamplePt = lsPt;
    theGeographicPt = geoPt;
 
    //compute fit
    initializeBilinear();
 
    //compute variance
    ossim_float64 sumerr2=0.0;
    vector<ossimDpt>::const_iterator i;
    vector<ossimGpt>::const_iterator j;
    ossimGpt gres;
    for(i=theLineSamplePt.begin() , j=theGeographicPt.begin() ; i != theLineSamplePt.end() ; ++i, ++j )
    {
       gres.lat  = theLatFit.lsFitValue(i->x, i->y);
       gres.lon  = theLonFit.lsFitValue(i->x, i->y);
       gres.hgt  = j->hgt; //same height as ground point
 
       sumerr2 += ( ossimEcefPoint(gres) - ossimEcefPoint(*j) ).norm2(); //add squared error in meters
    }
 
    return sumerr2 / theLineSamplePt.size(); //variance in meter^2
 }

◆ setupOptimizer()

bool ossimBilinearProjection::setupOptimizer ( const ossimString & setup )

virtual

setupFromString() Derived classes should implement as needed. Initialize parameters needed for optimizeFit and degreesOfFreedom

Returns: This implementation does nothing and returns false.

Reimplemented from ossimOptimizableProjection.

Definition at line 528 of file ossimBilinearProjection.cpp.

 {
    return false;
 }

◆ useForward()

virtual bool ossimBilinearProjection::useForward ( ) const

inlinevirtual

useForward() return true when it's better (more accurate / fast) to use forward (from ground to image) than inverse(from image to ground)

Implements ossimOptimizableProjection.

Definition at line 79 of file ossimBilinearProjection.h.

79 {return false;}

◆ worldToLineSample()

void ossimBilinearProjection::worldToLineSample	(	const ossimGpt &	worldPoint,
		ossimDpt &	lineSampPt
	)		const

virtual

METHOD: worldToLineSample() Performs the forward projection from ground point to line, sample.

Implements ossimProjection.

Definition at line 125 of file ossimBilinearProjection.cpp.

References ossimGpt::latd(), ossimGpt::lond(), ossimLeastSquaresBilin::lsFitValue(), ossimDpt::makeNan(), theInterpolationPointsHaveNanFlag, theInverseSupportedFlag, theXFit, theYFit, ossimProjection::worldToLineSample(), ossimDpt::x, and ossimDpt::y.

 {
    if(!theInverseSupportedFlag)
    {
       ossimProjection::worldToLineSample(worldPoint, lineSampPt);
    }
    else 
    {
       lineSampPt.makeNan();
       
       if (!theInterpolationPointsHaveNanFlag)
       {
          lineSampPt.x = theXFit.lsFitValue(worldPoint.lond(),
                                            worldPoint.latd());
          lineSampPt.y = theYFit.lsFitValue(worldPoint.lond(),
                                            worldPoint.latd());
       }
    }
 }

Member Data Documentation

◆ theGeographicPt

std::vector<ossimGpt> ossimBilinearProjection::theGeographicPt

protected

Definition at line 119 of file ossimBilinearProjection.h.

Referenced by getTiePoints(), gPtsHaveNan(), initializeBilinear(), lineSampleHeightToWorld(), loadState(), origin(), ossimBilinearProjection(), print(), saveState(), and setTiePoints().

◆ theInterpolationPointsHaveNanFlag

bool ossimBilinearProjection::theInterpolationPointsHaveNanFlag

protected

Definition at line 125 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), lineSampleHeightToWorld(), and worldToLineSample().

◆ theInverseSupportedFlag

bool ossimBilinearProjection::theInverseSupportedFlag

protected

Definition at line 126 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), and worldToLineSample().

◆ theLatFit

ossimLeastSquaresBilin ossimBilinearProjection::theLatFit

protected

Definition at line 120 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), lineSampleHeightToWorld(), and setTiePoints().

◆ theLineSamplePt

std::vector<ossimDpt> ossimBilinearProjection::theLineSamplePt

protected

Definition at line 118 of file ossimBilinearProjection.h.

Referenced by dPtsHaveNan(), getTiePoints(), initializeBilinear(), loadState(), midLineSamplePt(), ossimBilinearProjection(), print(), saveState(), and setTiePoints().

◆ theLonFit

ossimLeastSquaresBilin ossimBilinearProjection::theLonFit

protected

Definition at line 121 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), lineSampleHeightToWorld(), and setTiePoints().

◆ theXFit

ossimLeastSquaresBilin ossimBilinearProjection::theXFit

protected

Definition at line 122 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), and worldToLineSample().

◆ theYFit

ossimLeastSquaresBilin ossimBilinearProjection::theYFit

protected

Definition at line 123 of file ossimBilinearProjection.h.

Referenced by initializeBilinear(), and worldToLineSample().

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ossimBilinearProjection() [1/3]

◆ ossimBilinearProjection() [2/3]

◆ ossimBilinearProjection() [3/3]

◆ ~ossimBilinearProjection()

Member Function Documentation

◆ degreesOfFreedom()

◆ dPtsHaveNan()

◆ dup()

◆ getMetersPerPixel()

◆ getTiePoints()

◆ gPtsHaveNan()

◆ initializeBilinear()

◆ isAffectedByElevation()

◆ lineSampleHeightToWorld()

◆ lineSampleToWorld()

◆ loadState()

◆ midLineSamplePt()

◆ operator==()

◆ optimizeFit()

◆ origin()

◆ print()

◆ saveState()

◆ setTiePoints()

◆ setupOptimizer()

◆ useForward()

◆ worldToLineSample()

Member Data Documentation

◆ theGeographicPt

◆ theInterpolationPointsHaveNanFlag

◆ theInverseSupportedFlag

◆ theLatFit

◆ theLineSamplePt

◆ theLonFit

◆ theXFit

◆ theYFit