#include <ossimRsmModel.h>

Inheritance diagram for ossimRsmModel:

Public Member Functions
	ossimRsmModel ()
	default constructor More...

	ossimRsmModel (const ossimRsmModel &obj)
	copy constructor More...

const ossimRsmModel &	operator= (const ossimRsmModel &rhs)
	assignment operator More...

virtual void	worldToLineSample (const ossimGpt &world_point, ossimDpt &image_point) const
	worldToLineSample() Overrides base class implementation. More...

virtual void	lineSampleToWorld (const ossimDpt &image_point, ossimGpt &world_point) const
	lineSampleToWorld() Overrides base class pure virtual. More...

virtual void	lineSampleHeightToWorld (const ossimDpt &image_point, const double &heightEllipsoid, ossimGpt &worldPoint) const
	lineSampleHeightToWorld() Overrides base class pure virtual. More...

virtual void	imagingRay (const ossimDpt &image_point, ossimEcefRay &image_ray) const
	imagingRay() Overrides base class pure virtual. More...

virtual void	updateModel ()

virtual void	initAdjustableParameters ()

virtual ossimObject *	dup () const
	dup() Returns pointer to a new instance, copy of this. More...

virtual bool	useForward () const

virtual std::ostream &	print (std::ostream &out) const
	print() Extends base-class implementation. More...

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
	saveState Fulfills ossimObject base-class pure virtuals. More...

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
	loadState Fulfills ossimObject base-class pure virtuals. More...

Public Member Functions inherited from ossimSensorModel
	ossimSensorModel ()

	ossimSensorModel (const ossimSensorModel &copy_this)

	ossimSensorModel (const ossimKeywordlist &geom_kwl)

const ossimSensorModel &	operator= (const ossimSensorModel &rhs)
	assignment operator More...

virtual ossimObject *	getBaseObject ()

virtual const ossimObject *	getBaseObject () const

virtual ossimGpt	origin () const

virtual ossimDpt	getMetersPerPixel () const

virtual const double &	getNominalPosError () const
	Returns the estimated Absolute horizontal position error (CE90) of the sensor model. More...

virtual const double &	getRelativePosError () const
	Returns the estimated RELATIVE horizontal position error (CE90) of the sensor model. More...

virtual void	setNominalPosError (const double &ce90)
	Assigns the absolute image position error uncertainty (abs CE90) More...

virtual void	setRelativePosError (const double &ce90)
	Assigns the relative image position error uncertainty (rel CE90) More...

void	setRefImgPt (const ossimDpt &pt)
	Sets the center line sampe of the image. More...

void	setRefGndPt (const ossimGpt &pt)
	Sets the center latitude, longitude, height of the image. More...

void	setImageRect (const ossimDrect &imageRect)

void	setGroundRect (const ossimGpt &ul, const ossimGpt &ur, const ossimGpt &lr, const ossimGpt &ll)

ossimDpt	imageSize () const

void	setImageSize (const ossimDpt &size)

virtual void	adjustableParametersChanged ()

virtual bool	insideImage (const ossimDpt &p) const

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

const ossimString &	getImageID () const
	Access methods: More...

const ossimDrect &	getImageClipRect () const

virtual ossim_uint32	degreesOfFreedom () const

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

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

virtual ossimDpt	getForwardDeriv (int parmIdx, const ossimGpt &gpos, double hdelta=1e-11)

virtual ossimGpt	getInverseDeriv (int parmIdx, const ossimDpt &ipos, double hdelta=1e-11)

virtual ossimSensorModel::CovMatStatus	getObsCovMat (const ossimDpt &ipos, NEWMAT::SymmetricMatrix &Cov, const ossim_float64 defPointingSigma=0.5) const
	Gives 2X2 covariance matrix of observations. More...

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

void	computeGsd ()
	This method computes the ground sample distance(gsd) and sets class attributes theGSD and theMeanGSD by doing a lineSampleHeightToWorld on four points and calculating the distance from them. More...

virtual bool	getImageGeometry (const ossimString &, const ossimString &, ossimKeywordlist &) const
	Extracts geometry info from a non-ossim key,value pair to an ossim keyword list. 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	setupOptimizer (const ossimString &setup)
	setupFromString() Derived classes should implement as needed. Initialize parameters needed for optimizeFit and degreesOfFreedom More...

Public Member Functions inherited from ossimAdjustableParameterInterface
	ossimAdjustableParameterInterface ()

	ossimAdjustableParameterInterface (const ossimAdjustableParameterInterface &rhs)

virtual	~ossimAdjustableParameterInterface ()

void	newAdjustment (ossim_uint32 numberOfParameters=0)

void	setAdjustmentDescription (const ossimString &description)

ossimString	getAdjustmentDescription () const

ossimString	getAdjustmentDescription (ossim_uint32 adjustmentIdx) const
	Returns adjustmentDescription of specific adjustmentInfo. More...

void	setCurrentAdjustment (ossim_uint32 adjustmentIndex, bool notify=false)

bool	setCurrentAdjustment (const ossimString &description, bool notify=false)
	Sets the current adjustment to the adjustment with a matching description. More...

void	eraseAdjustment (bool notify)

void	eraseAdjustment (ossim_uint32 idx, bool notify)

void	resetAdjustableParameters (bool notify=false)

void	copyAdjustment (ossim_uint32 idx, bool notify)

void	copyAdjustment (bool notify=false)

void	keepAdjustment (ossim_uint32 idx, bool createCopy)

virtual void	keepAdjustment (bool createCopy=true)

const ossimAdjustableParameterInterface &	operator= (const ossimAdjustableParameterInterface &rhs)

void	removeAllAdjustments ()

ossim_uint32	getNumberOfAdjustableParameters () const

double	getAdjustableParameter (ossim_uint32 idx) const

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, bool notify=false)

virtual void	setAdjustableParameter (ossim_uint32 idx, double value, double sigma, bool notify=false)

double	getParameterSigma (ossim_uint32 idx) const

void	setParameterSigma (ossim_uint32 idx, double value, bool notify=false)

ossimUnitType	getParameterUnit (ossim_uint32 idx) const

void	setParameterUnit (ossim_uint32 idx, ossimUnitType unit)

void	setParameterUnit (ossim_uint32 idx, const ossimString &unit)

void	setParameterCenter (ossim_uint32 idx, double center, bool notify=false)

double	getParameterCenter (ossim_uint32 idx) const

double	computeParameterOffset (ossim_uint32 idx) const

void	setParameterOffset (ossim_uint32 idx, ossim_float64 value, bool notify=false)

ossimString	getParameterDescription (ossim_uint32 idx) const

void	setParameterDescription (ossim_uint32 idx, const ossimString &descrption)

ossim_int32	findParameterIdxGivenDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

ossim_int32	findParameterIdxContainingDescription (ossim_uint32 adjustmentIdx, const ossimString &name) const

bool	isParameterLocked (ossim_uint32 idx) const

void	setParameterLockFlag (ossim_uint32 idxParam, bool flag)

bool	getParameterLockFlag (ossim_uint32 idx) const

void	lockAllParametersCurrentAdjustment ()

void	unlockAllParametersCurrentAdjustment ()

void	lockAllParameters (ossim_uint32 idxAdjustment)

void	unlockAllParameters (ossim_uint32 idxAdjustment)

void	resizeAdjustableParameterArray (ossim_uint32 numberOfParameters)

void	setAdjustment (const ossimAdjustmentInfo &adj, bool notify=false)

void	setAdjustment (ossim_uint32 idx, const ossimAdjustmentInfo &adj, bool notify=false)

void	addAdjustment (const ossimAdjustmentInfo &adj, bool notify)

void	getAdjustment (ossimAdjustmentInfo &adj) const

void	getAdjustment (ossim_uint32 idx, ossimAdjustmentInfo &adj) const

ossim_uint32	getNumberOfAdjustments () const

ossim_uint32	getCurrentAdjustmentIdx () const

void	setDirtyFlag (bool flag=true)

void	setAllDirtyFlag (bool flag=true)

bool	hasDirtyAdjustments () const

virtual void	saveCurrentAdjustmentOnly (ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))
	Saves the current active adjustment to the KWL file. More...

bool	saveAdjustments (ossimKeywordlist &kwl, const ossimString &prefix=ossimString("")) const
	Save all adjustments to the KWL file. More...

bool	loadAdjustments (const ossimKeywordlist &kwl, const ossimString &prefix=ossimString(""))

std::ostream &	print (std::ostream &out) const
	Dumps the currently active adjustment to ostream. More...

Protected Member Functions
ossim_uint32	getPcaIndex (const double &x, const double &y, const double &z) const
	Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given ground point. More...

ossim_uint32	getPcaIndex (const ossimDpt &ipt, bool shiftPoint) const
	Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given image point. More...

void	lowOrderPolynomial (const double &x, const double &y, const double &z, ossimDpt &ipt) const
	Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given ground point. More...

double	polynomial (const double &x, const double &y, const double &z, const ossim_uint32 &maxx, const ossim_uint32 &maxy, const ossim_uint32 &maxz, std::vector< double > pcf) const

double	dPoly_dLat (const double &x, const double &y, const double &z, const ossim_uint32 &maxx, const ossim_uint32 &maxy, const ossim_uint32 &maxz, std::vector< double > pcf) const

double	dPoly_dLon (const double &x, const double &y, const double &z, const ossim_uint32 &maxx, const ossim_uint32 &maxy, const ossim_uint32 &maxz, std::vector< double > pcf) const

double	dPoly_dHgt (const double &x, const double &y, const double &z, const ossim_uint32 &maxx, const ossim_uint32 &maxy, const ossim_uint32 &maxz, std::vector< double > pcf) const

bool	validate () const
	Performs sanity check on key/required rsm data. More...

virtual	~ossimRsmModel ()
	virtual destructor More...

Protected Member Functions inherited from ossimSensorModel
virtual	~ossimSensorModel ()

virtual ossimDpt	extrapolate (const ossimGpt &gp) const

virtual ossimGpt	extrapolate (const ossimDpt &ip, const double &height=ossim::nan()) const

void	buildNormalEquation (const ossimTieGptSet &tieSet, NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &residue, NEWMAT::ColumnVector &projResidue, double pstep_scale)

NEWMAT::ColumnVector	getResidue (const ossimTieGptSet &tieSet)

NEWMAT::ColumnVector	solveLeastSquares (NEWMAT::SymmetricMatrix &A, NEWMAT::ColumnVector &r) const

NEWMAT::Matrix	invert (const NEWMAT::Matrix &m) const
	stable invert stolen from ossimRpcSolver More...

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Member Functions inherited from ossimAdjustableParameterInterface
bool	paramChanged (ossim_uint32 param_idx) const
	Returns true if specified parameter has been modified since last setAllChangeFlag(false) call. More...

void	setAllChangeFlags (bool areChanged)
	Sets all the change flags to the boolean indicated to indicate parameters are changed (TRUE) or not (FALSE). More...

void	initChangeFlags ()
	Initializes the change flags to TRUE. More...

Protected Attributes
ossimRsmida	m_ida

ossimRsmpia	m_pia

std::vector< ossimRsmpca >	m_pca

Protected Attributes inherited from ossimSensorModel
ossimIpt	theImageSize

ossimDpt	theSubImageOffset

ossimString	theImageID

ossimString	theSensorID

ossimDpt	theGSD

ossim_float64	theMeanGSD

ossimGpt	theRefGndPt

ossimDpt	theRefImgPt

ossimPolygon	theBoundGndPolygon

ossimDrect	theImageClipRect

ossim_float64	theRelPosError

ossim_float64	theNominalPosError

ossimDpt	theParWRTx
	Partials for current point. More...

ossimDpt	theParWRTy

ossimDpt	theParWRTz

ossimDpt	theObs
	Observations & residuals for current point. More...

ossimDpt	theResid

ossimRefPtr< ossimProjection >	theSeedFunction
	Used as an initial guess for iterative solutions and a guess for points outside the support bounds. More...

bool	theExtrapolateImageFlag

bool	theExtrapolateGroundFlag

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Additional Inherited Members
Public Types inherited from ossimSensorModel
enum	CovMatStatus { COV_INVALID = 0, COV_PARTIAL = 1, COV_FULL = 2 }

enum	DeriveMode { OBS_INIT =-99, EVALUATE =-98, P_WRT_X = -1, P_WRT_Y = -2, P_WRT_Z = -3 }

Static Public Member Functions inherited from ossimSensorModel
static void	writeGeomTemplate (ostream &os)

Detailed Description

Definition at line 17 of file ossimRsmModel.h.

Constructor & Destructor Documentation

◆ ossimRsmModel() [1/2]

ossimRsmModel::ossimRsmModel ( )

default constructor

Definition at line 41 of file ossimRsmModel.cpp.

References initAdjustableParameters().

Referenced by dup().

    :
    ossimSensorModel(),
    m_ida(),
    m_pia(),   
    m_pca()
 {
    initAdjustableParameters();
    
 }

◆ ossimRsmModel() [2/2]

ossimRsmModel::ossimRsmModel ( const ossimRsmModel & obj )

copy constructor

Definition at line 52 of file ossimRsmModel.cpp.

    :
    ossimSensorModel( obj ),
    m_ida( obj.m_ida ),
    m_pia( obj.m_pia ),
    m_pca( obj.m_pca )
 {
    
 }

◆ ~ossimRsmModel()

ossimRsmModel::~ossimRsmModel ( )

protectedvirtual

virtual destructor

Definition at line 74 of file ossimRsmModel.cpp.

75 {

76 }

Member Function Documentation

◆ dPoly_dHgt()

double ossimRsmModel::dPoly_dHgt	(	const double &	x,
		const double &	y,
		const double &	z,
		const ossim_uint32 &	maxx,
		const ossim_uint32 &	maxy,
		const ossim_uint32 &	maxz,
		std::vector< double >	pcf
	)		const

protected

Definition at line 643 of file ossimRsmModel.cpp.

References x, and y.

 {
    double dr = 0.0;
    ossim_uint32 index = 0;
    for (ossim_uint32 k = 0; k <= maxz; ++k)
    {
       for (ossim_uint32 j = 0; j <= maxy; ++j)
       {
          for (ossim_uint32 i = 0; i <= maxx; ++i)
          {
             if (k>0)
             {
                dr += k*pcf[index]*std::pow(x,(double)i) *
                   std::pow(y,(double)j)*std::pow(z,(double)(k-1));
             }
             ++index;
          } 
       }     
    }
    return dr;
 }

◆ dPoly_dLat()

double ossimRsmModel::dPoly_dLat	(	const double &	x,
		const double &	y,
		const double &	z,
		const ossim_uint32 &	maxx,
		const ossim_uint32 &	maxy,
		const ossim_uint32 &	maxz,
		std::vector< double >	pcf
	)		const

protected

Definition at line 595 of file ossimRsmModel.cpp.

References x, and y.

Referenced by lineSampleHeightToWorld().

 {
    double dr = 0.0;
    ossim_uint32 index = 0;
    for (ossim_uint32 k = 0; k <= maxz; ++k)
    {
       for (ossim_uint32 j = 0; j <= maxy; ++j)
       {
          for (ossim_uint32 i = 0; i <= maxx; ++i)
          {
             if (j>0)
             {
                dr+=j*pcf[index]*std::pow(x,(double)i)*std::pow(y,(double)(j-1))*std::pow(z,(double)k);
             }
             ++index;
          }
       }
    }
    return dr;
 }

◆ dPoly_dLon()

double ossimRsmModel::dPoly_dLon	(	const double &	x,
		const double &	y,
		const double &	z,
		const ossim_uint32 &	maxx,
		const ossim_uint32 &	maxy,
		const ossim_uint32 &	maxz,
		std::vector< double >	pcf
	)		const

protected

Definition at line 619 of file ossimRsmModel.cpp.

References x, and y.

Referenced by lineSampleHeightToWorld().

 {
    double dr = 0.0;
    ossim_uint32 index = 0;
    for (ossim_uint32 k = 0; k <= maxz; ++k)
    {
       for (ossim_uint32 j = 0; j <= maxy; ++j)
       {
          for (ossim_uint32 i = 0; i <= maxx; ++i)
          {
             if (i>0)
             {
                dr += i*pcf[index]*std::pow(x,(double)(i-1)) *
                   std::pow(y,(double)j)*std::pow(z,(double)k);
             }
             ++index;
          } 
       }     
    }
    return dr;
 }

◆ dup()

ossimObject * ossimRsmModel::dup ( ) const

virtual

dup() Returns pointer to a new instance, copy of this.

Implements ossimProjection.

Reimplemented in ossimNitfRsmModel.

Definition at line 397 of file ossimRsmModel.cpp.

References ossimRsmModel().

 {
    return new ossimRsmModel(*this);
 }

◆ getPcaIndex() [1/2]

ossim_uint32 ossimRsmModel::getPcaIndex	(	const double &	x,
		const double &	y,
		const double &	z
	)		const

protected

Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given ground point.

Parameters

x	Longitude in radians.
y	Latitude, in radians.
z	Height in meters.

Returns: Index into m_pca vector.

Definition at line 521 of file ossimRsmModel.cpp.

References lowOrderPolynomial(), x, and y.

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

 {
    ossimDpt ipt;
    lowOrderPolynomial( x, y, z, ipt );
    return getPcaIndex( ipt, false );
 }

◆ getPcaIndex() [2/2]

ossim_uint32 ossimRsmModel::getPcaIndex	(	const ossimDpt &	ipt,
		bool	shiftPoint
	)		const

protected

Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given image point.

Note: (0,0) in ossim space is (0.5, 0.5) in rsm space.

Parameters

ipt	Image line, sample.
shiftPoint	If true ipt is shifted +0.5 to go from ossim space(0 is center of pixel) to RSM space(0 is upper left of pixel). If false point is coming from rsm space and does not need to be shifted.

Returns: Index into m_pca vector.

Definition at line 529 of file ossimRsmModel.cpp.

References ossimRsmpia::m_cnis, ossimRsmpia::m_cssiz, m_ida, ossimRsmida::m_minc, ossimRsmida::m_minr, m_pia, ossimRsmpia::m_rnis, ossimRsmpia::m_rssiz, ossimDpt::x, and ossimDpt::y.

 {
    //---
    //  RSM (0,0) is upper left corner of pixel(0,0). OSSIM (0,0) is
    //  center of the pixel; hence, the shift 0.5 if coming from ossim.
    //---
    double shift = shiftPoint ? 0.5 : 0.0;
    
    // Row section number:
    double rsn = std::floor( ( ipt.y + shift - (double)(m_ida.m_minr) ) /
                             (double)(m_pia.m_rssiz) );
    if ( rsn < 0.0 )
    {
       rsn = 0.0;
    }
    else if ( rsn > (m_pia.m_rnis-1) )
    {
       rsn = m_pia.m_rnis-1;
    }
    // Column section number:
    double csn = std::floor( ( ipt.x + shift - (double)(m_ida.m_minc) ) / (double)(m_pia.m_cssiz) );
    if ( csn < 0.0 )
    {
       csn = 0.0;
    }
    else if ( csn > (m_pia.m_cnis-1) )
    {
       csn = m_pia.m_cnis-1;
    }
 
    //return static_cast<ossim_uint32>(rsn) * m_pia.m_rnis + static_cast<ossim_uint32>(csn);
    return static_cast<ossim_uint32>(rsn) * m_pia.m_cnis + static_cast<ossim_uint32>(csn);
 }

◆ imagingRay()

void ossimRsmModel::imagingRay	(	const ossimDpt &	image_point,
		ossimEcefRay &	image_ray
	)		const

virtual

imagingRay() Overrides base class pure virtual.

Reimplemented from ossimSensorModel.

Definition at line 356 of file ossimRsmModel.cpp.

References getPcaIndex(), lineSampleHeightToWorld(), and m_pca.

Referenced by lineSampleToWorld().

 {
    ossim_uint32 pcaIndex = getPcaIndex( imagePoint, true );
    
    //---
    // For "from point", "to point" we want the image ray to be from above the
    // ellipsoid down to Earth.
    // 
    // It appears the ray "from point" must be above the ellipsiod for the
    // ossimElevSource::intersectRay method; ultimately, the
    // ossimEllipsoid::nearestIntersection method, else it goes off in the
    // weeds...
    //---
    double vectorLength = m_pca[pcaIndex].m_znrmsf * 2.0;
 
    ossimGpt gpt;
 
    // "from" point
    double intHgt = m_pca[pcaIndex].m_znrmo + vectorLength;
    lineSampleHeightToWorld(imagePoint, intHgt, gpt);
    ossimEcefPoint intECFfrom(gpt);
 
    // "to" point
    lineSampleHeightToWorld(imagePoint, m_pca[pcaIndex].m_znrmo, gpt);
    ossimEcefPoint intECFto(gpt);
 
    // Construct ray
    ossimEcefRay ray(intECFfrom, intECFto);
 
    imageRay = ray;
 }

◆ initAdjustableParameters()

void ossimRsmModel::initAdjustableParameters ( )

virtual

Reimplemented from ossimAdjustableParameterInterface.

Definition at line 393 of file ossimRsmModel.cpp.

Referenced by ossimRsmModel().

394 {

395 }

◆ lineSampleHeightToWorld()

void ossimRsmModel::lineSampleHeightToWorld	(	const ossimDpt &	image_point,
		const double &	heightEllipsoid,
		ossimGpt &	worldPoint
	)		const

virtual

lineSampleHeightToWorld() Overrides base class pure virtual.

Height understood to be relative to standard ellipsoid.

Implements ossimSensorModel.

Definition at line 198 of file ossimRsmModel.cpp.

References dPoly_dLat(), dPoly_dLon(), getPcaIndex(), ossimGpt::hgt, ossim::isnan(), ossimGpt::lat, ossimGpt::lon, m_pca, ossimNotify(), ossimNotifyLevel_WARN, polynomial(), ossim::radiansToDegrees(), ossimGpt::wrap(), ossimDpt::x, and ossimDpt::y.

Referenced by imagingRay().

 {
    // Borrowed from ossimRpcModel algorithm to converge on polynomial roots
 
    //---
    // Constants for convergence tests:
    //---
    // SPEC says 1/20 of a pixel for polynomial fit, so converge to at least that point
    static const int    MAX_NUM_ITERATIONS  = 100;
    static const double CONVERGENCE_EPSILON = 0.05;  // pixels
 
    ossim_uint32 pcaIndex = getPcaIndex( image_point, true );
    
    // Image point of 0 to ossim is 0.5 to RSM.
    // double U    = (image_point.y-m_rnrmo) / (m_rnrmsf);
    // double V    = (image_point.x-m_cnrmo) / (m_cnrmsf);
 
    double U = (image_point.y+0.5-m_pca[pcaIndex].m_rnrmo) / (m_pca[pcaIndex].m_rnrmsf);
    double V = (image_point.x+0.5-m_pca[pcaIndex].m_cnrmo) / (m_pca[pcaIndex].m_cnrmsf);
 
    //---
    // Initialize quantities to be used in the iteration for ground point:
    //---
    double nlat = 0.0;  // normalized latitude
    double nlon = 0.0;  // normalized longitude
    double nhgt;
 
    if(ossim::isnan(ellHeight))
    {
      nhgt = (- m_pca[pcaIndex].m_znrmo) / m_pca[pcaIndex].m_znrmsf;  // norm height
    }
    else
    {
       nhgt = (ellHeight - m_pca[pcaIndex].m_znrmo) / m_pca[pcaIndex].m_znrmsf;  // norm height
    }
 
    double epsilonU = CONVERGENCE_EPSILON/m_pca[pcaIndex].m_rnrmsf;
    double epsilonV = CONVERGENCE_EPSILON/m_pca[pcaIndex].m_cnrmsf;
    int    iteration = 0;
    //---
    // Declare variables only once outside the loop. These include:
    // * polynomials (numerators Pu, Pv, and denominators Qu, Qv),
    // * partial derivatives of polynomials wrt X, Y,
    // * computed normalized image point: Uc, Vc,
    // * residuals of normalized image point: deltaU, deltaV,
    // * partial derivatives of Uc and Vc wrt X, Y,
    // * corrections to normalized lat, lon: deltaLat, deltaLon.
    //---
    double Pu, Qu, Pv, Qv;
    double dPu_dLat, dQu_dLat, dPv_dLat, dQv_dLat;
    double dPu_dLon, dQu_dLon, dPv_dLon, dQv_dLon;
    double Uc, Vc;
    double deltaU, deltaV;
    double dU_dLat, dU_dLon, dV_dLat, dV_dLon, W;
    double deltaLat, deltaLon;
 
    //---
    // Now iterate until the computed Uc, Vc is within epsilon of the desired
    // image point U, V:
    //---
    do
    {
       //---
       // Calculate the normalized line and sample Uc, Vc as ratio of
       // polynomials Pu, Qu and Pv, Qv:
       //---
       Pu = polynomial(nlon, nlat, nhgt, m_pca[pcaIndex].m_rnpwrx, m_pca[pcaIndex].m_rnpwry, m_pca[pcaIndex].m_rnpwrz, m_pca[pcaIndex].m_rnpcf);
       Qu = polynomial(nlon, nlat, nhgt, m_pca[pcaIndex].m_rdpwrx, m_pca[pcaIndex].m_rdpwry, m_pca[pcaIndex].m_rdpwrz, m_pca[pcaIndex].m_rdpcf);
       Pv = polynomial(nlon, nlat, nhgt, m_pca[pcaIndex].m_cnpwrx, m_pca[pcaIndex].m_cnpwry, m_pca[pcaIndex].m_cnpwrz, m_pca[pcaIndex].m_cnpcf);
       Qv = polynomial(nlon, nlat, nhgt, m_pca[pcaIndex].m_cdpwrx, m_pca[pcaIndex].m_cdpwry, m_pca[pcaIndex].m_cdpwrz, m_pca[pcaIndex].m_cdpcf);
       Uc = Pu/Qu;
       Vc = Pv/Qv;
 
       //---
       // Compute residuals between desired and computed line, sample:
       //---
       deltaU = U - Uc;
       deltaV = V - Vc;
 
       //---
       // Check for convergence and skip re-linearization if converged:
       //---
       if ((fabs(deltaU) > epsilonU) || (fabs(deltaV) > epsilonV))
       {
          //---
          // Analytically compute the partials of each polynomial wrt lat, lon:
          //---
          dPu_dLat = dPoly_dLat(nlon, nlat, nhgt, m_pca[pcaIndex].m_rnpwrx, m_pca[pcaIndex].m_rnpwry, m_pca[pcaIndex].m_rnpwrz, m_pca[pcaIndex].m_rnpcf);
          dQu_dLat = dPoly_dLat(nlon, nlat, nhgt, m_pca[pcaIndex].m_rdpwrx, m_pca[pcaIndex].m_rdpwry, m_pca[pcaIndex].m_rdpwrz, m_pca[pcaIndex].m_rdpcf);
          dPv_dLat = dPoly_dLat(nlon, nlat, nhgt, m_pca[pcaIndex].m_cnpwrx, m_pca[pcaIndex].m_cnpwry, m_pca[pcaIndex].m_cnpwrz, m_pca[pcaIndex].m_cnpcf);
          dQv_dLat = dPoly_dLat(nlon, nlat, nhgt, m_pca[pcaIndex].m_cdpwrx, m_pca[pcaIndex].m_cdpwry, m_pca[pcaIndex].m_cdpwrz, m_pca[pcaIndex].m_cdpcf);
          dPu_dLon = dPoly_dLon(nlon, nlat, nhgt, m_pca[pcaIndex].m_rnpwrx, m_pca[pcaIndex].m_rnpwry, m_pca[pcaIndex].m_rnpwrz, m_pca[pcaIndex].m_rnpcf);
          dQu_dLon = dPoly_dLon(nlon, nlat, nhgt, m_pca[pcaIndex].m_rdpwrx, m_pca[pcaIndex].m_rdpwry, m_pca[pcaIndex].m_rdpwrz, m_pca[pcaIndex].m_rdpcf);
          dPv_dLon = dPoly_dLon(nlon, nlat, nhgt, m_pca[pcaIndex].m_cnpwrx, m_pca[pcaIndex].m_cnpwry, m_pca[pcaIndex].m_cnpwrz, m_pca[pcaIndex].m_cnpcf);
          dQv_dLon = dPoly_dLon(nlon, nlat, nhgt, m_pca[pcaIndex].m_cdpwrx, m_pca[pcaIndex].m_cdpwry, m_pca[pcaIndex].m_cdpwrz, m_pca[pcaIndex].m_cdpcf);
 
          //---
          // Analytically compute partials of quotients U and V wrt lat, lon:
          //---
          dU_dLat = (Qu*dPu_dLat - Pu*dQu_dLat)/(Qu*Qu);
          dU_dLon = (Qu*dPu_dLon - Pu*dQu_dLon)/(Qu*Qu);
          dV_dLat = (Qv*dPv_dLat - Pv*dQv_dLat)/(Qv*Qv);
          dV_dLon = (Qv*dPv_dLon - Pv*dQv_dLon)/(Qv*Qv);
 
          W = dU_dLon*dV_dLat - dU_dLat*dV_dLon;
 
          //---
          // Now compute the corrections to normalized lat, lon:
          //---
          deltaLat = (dU_dLon*deltaV - dV_dLon*deltaU) / W;
          deltaLon = (dV_dLat*deltaU - dU_dLat*deltaV) / W;
          nlat += deltaLat;
          nlon += deltaLon;
       }
 
       ++iteration;
 
    } while (((fabs(deltaU)>epsilonU) || (fabs(deltaV)>epsilonV))
             && (iteration < MAX_NUM_ITERATIONS));
    //---
    // Test for exceeding allowed number of iterations. Flag error if so:
    //---
    if (iteration == MAX_NUM_ITERATIONS)
    {
       ossimNotify(ossimNotifyLevel_WARN)
          << "WARNING ossimRsmModel::lineSampleHeightToWorld:\n"
          << "Max number of iterations reached in ground point "
          << "solution. Results are inaccurate." << endl;
    }
 
    //---
    // Now un-normalize the ground point lat, lon and establish return quantity:
    //
    // lon will 0 to 2PI when image is near PI radians as specified in RSMIDA
    // GRNDD field when value is "H" versus "G".   OSSIMGPT wrap handles this
    // automatically, so no need to worry about it.
    //---
    gpt.lat = ossim::radiansToDegrees(nlat*m_pca[pcaIndex].m_ynrmsf + m_pca[pcaIndex].m_ynrmo);
    gpt.lon = ossim::radiansToDegrees(nlon*m_pca[pcaIndex].m_xnrmsf + m_pca[pcaIndex].m_xnrmo);
    gpt.hgt = (nhgt * m_pca[pcaIndex].m_znrmsf) + m_pca[pcaIndex].m_znrmo; //ellHeight;
 
    //---
    // Note: See above note. Added in wrap call. Longitude was coming out 242
    // when should have been -118. (drb - 22 May 2015)
    //---
    gpt.wrap();
    
 } // End: ossimRsmModel::lineSampleHeightToWorld( ... )

◆ lineSampleToWorld()

void ossimRsmModel::lineSampleToWorld	(	const ossimDpt &	image_point,
		ossimGpt &	world_point
	)		const

virtual

lineSampleToWorld() Overrides base class pure virtual.

Intersects DEM.

Reimplemented from ossimSensorModel.

Definition at line 167 of file ossimRsmModel.cpp.

References ossimDpt::hasNans(), imagingRay(), ossimElevManager::instance(), ossimElevSource::intersectRay(), ossimGpt::makeNan(), ossimDpt::x, and ossimDpt::y.

 {
    if(!imagePoint.hasNans())
    {
       ossimEcefRay ray;
       
       //---
       // Note:
       // RSM (0,0) is upper left corner of pixel(0,0). OSSIM (0,0) is
       //  center of the pixel; hence, the + 0.5. (drb 22 May 2015)
       //---
       // imagingRay(imagePoint, ray);
       imagingRay(ossimDpt(imagePoint.x+0.5, imagePoint.y+0.5), ray);
       ossimElevManager::instance()->intersectRay(ray, worldPoint);
    }
    else
    {
       worldPoint.makeNan();
    }
 }

◆ loadState()

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

virtual

loadState Fulfills ossimObject base-class pure virtuals.

Loads and saves geometry KWL files. Returns true if successful.

Reimplemented from ossimSensorModel.

Reimplemented in ossimNitfRsmModel.

Definition at line 454 of file ossimRsmModel.cpp.

 {
    static const char MODULE[] = "ossimRsmModel::loadState";
    if (traceExec())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << MODULE << " entered...\n";
    }
 
    bool status = false;
    
    // Check for type match before preceeding:
    std::string pfx = ( prefix ? prefix : "" );
    std::string type = kwl.findKey( pfx, std::string(ossimKeywordNames::TYPE_KW) );
    if ( (type == "ossimNitfRsmModel" ) || ( type == MODEL_TYPE_KW ) )
    {
       // Pass on to the base-class for parsing first:
       if ( ossimSensorModel::loadState(kwl, prefix) )
       {
          if ( m_ida.loadState( kwl, pfx ) )
          {
             if ( m_pia.loadState( kwl, pfx ) )
             {
                m_pca.clear();
                
                for ( ossim_uint32 tagIndex = 0; tagIndex < m_pia.m_tnis; ++tagIndex )
                {
                   ossimRsmpca pca;
                   if ( pca.loadState( kwl, pfx, tagIndex ) )
                   {
                      m_pca.push_back( pca );
                   }
                   else
                   {
                      ossimNotify(ossimNotifyLevel_WARN)
                         << "WARNING! RSMPCA[" << tagIndex << "] intitialization failed!"
                         << std::endl;
                      break; // Get out...
                   }
                }
 
                // Should now have a rsmpca record for each segment.
                if ( m_pia.m_tnis == (ossim_uint32)m_pca.size() )
                {
                   // Set the status for downstream code.
                   status = true;
                   
                   updateModel();
                }
                
             } // Matches: if ( m_pia.loadState( kwl, pfx ) )
             
          } // Matches:if ( m_ida.loadState( kwl, pfx ) ) 
          
       } // Matches: if ( ossimSensorModel::loadState(kwl, prefix) )
       
    } // Matches: if ( (type == "ossimNitfRsmModel" ) || ...
    
    if (traceExec())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << MODULE << " exit status = " << (status?"true":"false") << "\n";
    }
    
    return status;
 }

◆ lowOrderPolynomial()

void ossimRsmModel::lowOrderPolynomial	(	const double &	x,
		const double &	y,
		const double &	z,
		ossimDpt &	ipt
	)		const

protected

Gets index into RSM Polynomial Coefficients(rsmpca) container array for a given ground point.

Parameters

x	Longitude in radians.
y	Latitude, in radians.
z	Height in meters. ipt Initialized by this.

Definition at line 563 of file ossimRsmModel.cpp.

References ossimRsmpia::m_c0, ossimRsmpia::m_cx, ossimRsmpia::m_cxx, ossimRsmpia::m_cxy, ossimRsmpia::m_cxz, ossimRsmpia::m_cy, ossimRsmpia::m_cyy, ossimRsmpia::m_cyz, ossimRsmpia::m_cz, ossimRsmpia::m_czz, m_pia, ossimRsmpia::m_r0, ossimRsmpia::m_rx, ossimRsmpia::m_rxx, ossimRsmpia::m_rxy, ossimRsmpia::m_rxz, ossimRsmpia::m_ry, ossimRsmpia::m_ryy, ossimRsmpia::m_ryz, ossimRsmpia::m_rz, ossimRsmpia::m_rzz, ossimDpt::x, x, ossimDpt::y, and y.

Referenced by getPcaIndex().

 {
    ipt.y = m_pia.m_r0 + m_pia.m_rx * x + m_pia.m_ry * y + m_pia.m_rz * z +
       m_pia.m_rxx * x * x + m_pia.m_rxy * x * y + m_pia.m_rxz * x * z +
       m_pia.m_ryy * y * y + m_pia.m_ryz * y * z + m_pia.m_rzz * z * z;
 
    ipt.x = m_pia.m_c0 + m_pia.m_cx * x + m_pia.m_cy * y + m_pia.m_cz * z +
       m_pia.m_cxx * x * x + m_pia.m_cxy * x * y + m_pia.m_cxz * x * z +
       m_pia.m_cyy * y * y + m_pia.m_cyz * y * z + m_pia.m_czz * z * z;
 }

◆ operator=()

const ossimRsmModel & ossimRsmModel::operator= ( const ossimRsmModel & rhs )

assignment operator

Definition at line 62 of file ossimRsmModel.cpp.

References m_ida, m_pca, m_pia, and ossimSensorModel::operator=().

Referenced by ossimNitfRsmModel::operator=().

 {
    if (this != &rhs)
    {
       ossimSensorModel::operator=(rhs);
       m_ida = rhs.m_ida;
       m_pia = rhs.m_pia;      
       m_pca = rhs.m_pca;
    }
    return *this;  
 }

◆ polynomial()

double ossimRsmModel::polynomial	(	const double &	x,
		const double &	y,
		const double &	z,
		const ossim_uint32 &	maxx,
		const ossim_uint32 &	maxy,
		const ossim_uint32 &	maxz,
		std::vector< double >	pcf
	)		const

protected

Definition at line 575 of file ossimRsmModel.cpp.

References x, and y.

Referenced by lineSampleHeightToWorld(), and worldToLineSample().

 {
    double r = 0.0;
    ossim_uint32 index = 0;
    for (ossim_uint32 k = 0; k <= maxz; ++k)
    {
       for (ossim_uint32 j = 0; j <= maxy; ++j)
       {
          for (ossim_uint32 i = 0; i <= maxx; ++i)
          {
             r+=pcf[index]*std::pow(x,(double)i)*std::pow(y,(double)j)*std::pow(z,(double)k);
             ++index;
          }
       }
    }
    return r;
 }

◆ print()

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

virtual

print() Extends base-class implementation.

Dumps contents of object to ostream.

Reimplemented from ossimSensorModel.

Definition at line 407 of file ossimRsmModel.cpp.

References saveState().

 {
    std::string prefix = "";
    ossimKeywordlist kwl;
    saveState( kwl, prefix.c_str() );
    out << kwl;
    return out;
 }

◆ saveState()

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

virtual

saveState Fulfills ossimObject base-class pure virtuals.

Loads and saves geometry KWL files. Returns true if successful.

Reimplemented from ossimSensorModel.

Reimplemented in ossimNitfRsmModel.

Definition at line 416 of file ossimRsmModel.cpp.

Referenced by print().

 {
    static const char MODULE[] = "ossimRsmModel::saveState";
    if (traceExec())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << MODULE << " entered...\n";
    }
 
    kwl.add(prefix, ossimKeywordNames::TYPE_KW, MODEL_TYPE_KW.c_str());
 
    //---
    // Hand off to base class for common stuff:
    //---
    ossimSensorModel::saveState(kwl, prefix);
 
    std::string pfx = (prefix ? prefix : "" );
 
    // IDA:
    m_ida.saveState( kwl, prefix );
 
    // PIA:
    m_pia.saveState( kwl, prefix );  
 
    // PCA:
    for ( ossim_uint32 i = 0; i < m_pca.size(); ++i )
    {
       m_pca[i].saveState( kwl, pfx, i );
    }
    
    if (traceExec())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << MODULE << " exited...\n";
    }
 
    return true;
 }

◆ updateModel()

void ossimRsmModel::updateModel ( )

virtual

VIRTUAL METHOD: updateModel() Following a change to the adjustable parameter set, this virtual is called to permit instances to compute derived quantities after parameter change.

Reimplemented from ossimSensorModel.

Definition at line 389 of file ossimRsmModel.cpp.

390 {

391 }

◆ useForward()

virtual bool ossimRsmModel::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 70 of file ossimRsmModel.h.

70 {return false;}

◆ validate()

bool ossimRsmModel::validate ( ) const

protected

Performs sanity check on key/required rsm data.

Returns: true on success, false on error.

Definition at line 667 of file ossimRsmModel.cpp.

References ossimRsmpia::m_cssiz, ossimRsmida::m_grndd, m_ida, m_pca, m_pia, ossimRsmpia::m_rnis, ossimRsmpia::m_rssiz, ossimRsmpia::m_tnis, ossimNotify(), ossimNotifyLevel_WARN, and status.

 {
    static const char MODULE[] = "ossimRsmModel::validate";
 
    bool status = true;
 
    if (  (m_pia.m_rnis == 0) ||  (m_pia.m_rnis == 0) || (m_pia.m_tnis == 0) )
    {
       status = false;
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmpia must have at least one section!" << std::endl;
    }
 
    if ( m_pca.size() != m_pia.m_tnis )
    {
       status = false;
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmpca array not equal to section count!" << std::endl;
    }
    if ( ( m_ida.m_grndd != 'G' ) && ( m_ida.m_grndd != 'H' ) && ( m_ida.m_grndd != 'R' ) )
    {
       status = false;
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmida grndd Ground Domain Form not set!" << std::endl;
    }
    if ( m_ida.m_grndd == 'R' )
    {
       status = false;
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmida grndd Rectangular Ground Domain not supported!" << std::endl;
    }
    if ( m_pia.m_rssiz == 0 )
    {
       status = false; // divide by zero.
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmpia rrsiz Section row size cannot be zero!" << std::endl;
    }
    if ( m_pia.m_cssiz == 0 )
    {
       status = false; // divide by zero.
       ossimNotify(ossimNotifyLevel_WARN)
          << MODULE
          << " ERROR: rsmpia cssiz Section column size cannot be zero!" << std::endl;
    }
 
    return status;
 }

◆ worldToLineSample()

void ossimRsmModel::worldToLineSample	(	const ossimGpt &	world_point,
		ossimDpt &	image_point
	)		const

virtual

worldToLineSample() Overrides base class implementation.

Directly computes line-sample from the polynomials.

Reimplemented from ossimSensorModel.

Definition at line 84 of file ossimRsmModel.cpp.

References ossim::degreesToRadians(), getPcaIndex(), ossimGpt::hgt, ossimGpt::isHgtNan(), ossimGpt::isLatNan(), ossimGpt::isLonNan(), ossim::isnan(), ossimGpt::lat, ossimDpt::line, ossimGpt::lon, ossimRsmida::m_grndd, m_ida, m_pca, ossimDpt::makeNan(), polynomial(), ossimDpt::samp, x, and y.

 {
    if(ground_point.isLatNan() || ground_point.isLonNan() )
    {
       img_pt.makeNan();
       return;
    }
 
    //---
    // RSMIDA GRNDD Ground Domain Form:
    // G: Geodetic: range x is -pi to pi, y is -pi/2 to pi/2
    // H: Geodetic: range x is 0 to 2pi, y is -pi/2 to pi/2 (where image is close to pi.
    // R: Rectangular (not supported).
    //---
 
    // Initial xyz for computing the pca index.
    double x;
    if ( m_ida.m_grndd == 'H' )
    {
       x = ossim::degreesToRadians((ground_point.lon >= 0.0) ?
                                   ground_point.lon : ground_point.lon + 360.0);
    }
    else
    {
       x = ossim::degreesToRadians( ground_point.lon );
    }
       
    double y = ossim::degreesToRadians(ground_point.lat);
    double z = ground_point.hgt;
    if ( ossim::isnan( z ) )
    {
       z = 0.0; // ??? drb
    }
 
    ossim_uint32 pcaIndex = getPcaIndex( x, y, z );
 
    //---
    // Normalize the lat, lon, hgt:
    // a_norm = (a-offset)/scalefactor
    //---
    y = (y - m_pca[pcaIndex].m_ynrmo) / m_pca[pcaIndex].m_ynrmsf;
    x = (x - m_pca[pcaIndex].m_xnrmo) / m_pca[pcaIndex].m_xnrmsf;   
    if( ground_point.isHgtNan() )
    {
      z = ( - m_pca[pcaIndex].m_znrmo) / m_pca[pcaIndex].m_znrmsf;
    }
    else
    {
      z = (ground_point.hgt - m_pca[pcaIndex].m_znrmo) / m_pca[pcaIndex].m_znrmsf;
    }
 
    double rnNrm =  polynomial(x, y, z, m_pca[pcaIndex].m_rnpwrx, m_pca[pcaIndex].m_rnpwry, m_pca[pcaIndex].m_rnpwrz, m_pca[pcaIndex].m_rnpcf); 
    double rdNrm =  polynomial(x, y, z, m_pca[pcaIndex].m_rdpwrx, m_pca[pcaIndex].m_rdpwry, m_pca[pcaIndex].m_rdpwrz, m_pca[pcaIndex].m_rdpcf);
    double cnNrm =  polynomial(x, y, z, m_pca[pcaIndex].m_cnpwrx, m_pca[pcaIndex].m_cnpwry, m_pca[pcaIndex].m_cnpwrz, m_pca[pcaIndex].m_cnpcf);
    double cdNrm =  polynomial(x, y, z, m_pca[pcaIndex].m_cdpwrx, m_pca[pcaIndex].m_cdpwry, m_pca[pcaIndex].m_cdpwrz, m_pca[pcaIndex].m_cdpcf);
    
    double rNrm  = rnNrm / rdNrm;
    double cNrm  = cnNrm / cdNrm;
 
    //---
    //  Unnormalize the computed value
    //  a = (a_norm * scalefactor) + offset
    //
    //  Note:
    //
    //  RSM (0,0) is upper left corner of pixel(0,0). OSSIM (0,0) is
    //  center of the pixel; hence, the - 0.5. (drb 22 May 2015)
    //---
    
    // img_pt.line = (rNrm * m_rnrmsf) + m_rnrmo; 
    // img_pt.samp = (cNrm * m_cnrmsf) + m_cnrmo; 
    img_pt.line = (rNrm * m_pca[pcaIndex].m_rnrmsf) + m_pca[pcaIndex].m_rnrmo - 0.5; 
    img_pt.samp = (cNrm * m_pca[pcaIndex].m_cnrmsf) + m_pca[pcaIndex].m_cnrmo - 0.5; 
 
 
 } // End: ossimRsmModel::worldToLineSample( ... )

Member Data Documentation

◆ m_ida

ossimRsmida ossimRsmModel::m_ida

protected

Definition at line 165 of file ossimRsmModel.h.

Referenced by getPcaIndex(), operator=(), validate(), and worldToLineSample().

◆ m_pca

std::vector<ossimRsmpca> ossimRsmModel::m_pca

protected

Definition at line 167 of file ossimRsmModel.h.

Referenced by imagingRay(), lineSampleHeightToWorld(), operator=(), validate(), and worldToLineSample().

◆ m_pia

ossimRsmpia ossimRsmModel::m_pia

protected

Definition at line 166 of file ossimRsmModel.h.

Referenced by getPcaIndex(), lowOrderPolynomial(), operator=(), and validate().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimRsmModel() [1/2]

◆ ossimRsmModel() [2/2]

◆ ~ossimRsmModel()

Member Function Documentation

◆ dPoly_dHgt()

◆ dPoly_dLat()

◆ dPoly_dLon()

◆ dup()

◆ getPcaIndex() [1/2]

◆ getPcaIndex() [2/2]

◆ imagingRay()

◆ initAdjustableParameters()

◆ lineSampleHeightToWorld()

◆ lineSampleToWorld()

◆ loadState()

◆ lowOrderPolynomial()

◆ operator=()

◆ polynomial()

◆ print()

◆ saveState()

◆ updateModel()

◆ useForward()

◆ validate()

◆ worldToLineSample()

Member Data Documentation

◆ m_ida

◆ m_pca

◆ m_pia