#include <ossimLandSatModel.h>

Inheritance diagram for ossimLandSatModel:

Public Types
enum	ProjectionType { UNKNOWN_PROJECTION = 0, UTM_MAP, UTM_ORBIT, SOM_MAP, SOM_ORBIT }

enum	AdjustParamIndex { INTRACK_OFFSET = 0, CRTRACK_OFFSET, LINE_GSD_CORR, SAMP_GSD_CORR, ROLL_OFFSET, YAW_OFFSET, YAW_RATE, MAP_ROTATION, NUM_ADJUSTABLE_PARAMS }

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 }

Public Member Functions
	ossimLandSatModel ()

	ossimLandSatModel (const ossimFfL7 &head)

	ossimLandSatModel (const ossimFilename &init_file)

	ossimLandSatModel (const ossimKeywordlist &geom_kwl)

	ossimLandSatModel (const ossimLandSatModel &rhs)

virtual ossimObject *	dup () const

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

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

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

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

virtual void	imagingRay (const ossimDpt &image_point, ossimEcefRay &image_ray) const

virtual void	updateModel ()

virtual bool	useForward () const

virtual bool	setupOptimizer (const ossimString &init_file)
	image to ground faster 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...

virtual void	lineSampleToWorld (const ossimDpt &image_point, ossimGpt &world_point) const

virtual void	worldToLineSample (const ossimGpt &world_point, ossimDpt &image_point) const

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)

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...

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

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

Protected Member Functions
virtual	~ossimLandSatModel ()
	uses file path to init model More...

void	initFromHeader (const ossimFfL7 &head)

virtual void	initAdjustableParameters ()

void	initMapProjection ()

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
double	theIllumAzimuth

double	theIllumElevation

double	theOrbitAltitude

double	theOrbitInclination

int	theMapZone

ossimDpt	theMapOffset

int	theWrsPathNumber

int	theWrsRowNumber

double	theMeridianalAngle

double	thePositionError

ProjectionType	theProjectionType

ossimRefPtr< ossimMapProjection >	theMapProjection

double	theMapAzimAngle

double	theMapAzimCos

double	theMapAzimSin

double	theMap2IcRotAngle

double	theMap2IcRotCos

double	theMap2IcRotSin

double	theIntrackOffset

double	theCrtrackOffset

double	theLineGsdCorr

double	theSampGsdCorr

double	theRollOffset

double	theYawOffset

double	theYawRate

double	theMapRotation

NEWMAT::Matrix	theRollRotMat

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

Detailed Description

Definition at line 39 of file ossimLandSatModel.h.

Member Enumeration Documentation

◆ AdjustParamIndex

enum ossimLandSatModel::AdjustParamIndex

Enumerator
INTRACK_OFFSET
CRTRACK_OFFSET
LINE_GSD_CORR
SAMP_GSD_CORR
ROLL_OFFSET
YAW_OFFSET
YAW_RATE
MAP_ROTATION
NUM_ADJUSTABLE_PARAMS

Definition at line 61 of file ossimLandSatModel.h.

    {
       INTRACK_OFFSET = 0,
       CRTRACK_OFFSET,
       LINE_GSD_CORR,
       SAMP_GSD_CORR,
       ROLL_OFFSET,
       YAW_OFFSET,
       YAW_RATE,
       MAP_ROTATION,
       NUM_ADJUSTABLE_PARAMS // not an index
    };

◆ ProjectionType

enum ossimLandSatModel::ProjectionType

Enumerator
UNKNOWN_PROJECTION
UTM_MAP
UTM_ORBIT
SOM_MAP
SOM_ORBIT

Definition at line 52 of file ossimLandSatModel.h.

    {
       UNKNOWN_PROJECTION = 0,
       UTM_MAP,
       UTM_ORBIT,
       SOM_MAP,
       SOM_ORBIT
    };

Constructor & Destructor Documentation

◆ ossimLandSatModel() [1/5]

ossimLandSatModel::ossimLandSatModel ( )

CONSTRUCTORS:

Definition at line 104 of file ossimLandSatModel.cpp.

Referenced by dup().

    :
    ossimSensorModel(),
    theIntrackOffset     (0.0),
    theCrtrackOffset     (0.0),
    theLineGsdCorr       (0.0),   
    theSampGsdCorr       (0.0),
    theRollOffset        (0.0),
    theYawOffset         (0.0),
    theYawRate           (0.0),
    theMapRotation       (0.0)
 
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel: entering..." << std::endl;
 
    initAdjustableParameters();
    
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel: returning..." << std::endl;
 }

◆ ossimLandSatModel() [2/5]

ossimLandSatModel::ossimLandSatModel ( const ossimFfL7 & head )

Definition at line 128 of file ossimLandSatModel.cpp.

    :
    ossimSensorModel(),
    theIntrackOffset     (0.0),
    theCrtrackOffset     (0.0),
    theLineGsdCorr       (0.0),   
    theSampGsdCorr       (0.0),
    theRollOffset        (0.0),
    theYawOffset         (0.0),
    theYawRate           (0.0),
    theMapRotation       (0.0)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(head): entering..." << std::endl;   
 
    initFromHeader(head);
 
    if (traceDebug()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(head): Exited..." << std::endl;
 }

◆ ossimLandSatModel() [3/5]

ossimLandSatModel::ossimLandSatModel ( const ossimFilename & init_file )

Definition at line 155 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(init_file): entering..." << std::endl;
 
    setupOptimizer(init_file);
 
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(init_file): Exited..." << std::endl;
 }

◆ ossimLandSatModel() [4/5]

ossimLandSatModel::ossimLandSatModel ( const ossimKeywordlist & geom_kwl )

Definition at line 173 of file ossimLandSatModel.cpp.

    :
    ossimSensorModel(),
    theIntrackOffset     (0.0),
    theCrtrackOffset     (0.0),
    theLineGsdCorr       (0.0),   
    theSampGsdCorr       (0.0),
    theRollOffset        (0.0),
    theYawOffset         (0.0),
    theYawRate           (0.0),
    theMapRotation       (0.0)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(geom_kwl): entering..." << std::endl;
 
    initAdjustableParameters();
 
    //***
    // Parse keywordlist for geometry:
    //***
    loadState(geom_kwl);
 
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(geom_kwl): Exited..." << std::endl;
 }

◆ ossimLandSatModel() [5/5]

ossimLandSatModel::ossimLandSatModel ( const ossimLandSatModel & rhs )

Definition at line 200 of file ossimLandSatModel.cpp.

    :
    ossimSensorModel   (rhs),
    theIllumAzimuth    (rhs.theIllumAzimuth),
    theIllumElevation  (rhs.theIllumElevation),
    theOrbitAltitude   (rhs.theOrbitAltitude),
    theOrbitInclination(rhs.theOrbitInclination),
    theMapZone         (rhs.theMapZone),
    theMapOffset       (rhs.theMapOffset),
    theWrsPathNumber   (rhs.theWrsPathNumber),
    theWrsRowNumber    (rhs.theWrsRowNumber),
    theMeridianalAngle (rhs.theMeridianalAngle),
    thePositionError   (rhs.thePositionError),
    theProjectionType  (rhs.theProjectionType),
    theMapProjection   (rhs.theMapProjection.valid()?
                        (ossimMapProjection*)rhs.theMapProjection->dup():
                        (ossimMapProjection*)0),
    theMapAzimAngle    (rhs.theMapAzimAngle),
    theMapAzimCos      (rhs.theMapAzimCos),
    theMapAzimSin      (rhs.theMapAzimSin),
    theMap2IcRotAngle  (rhs.theMap2IcRotAngle),
    theMap2IcRotCos    (rhs.theMap2IcRotCos),
    theMap2IcRotSin    (rhs.theMap2IcRotSin),
    theIntrackOffset   (rhs.theIntrackOffset),
    theCrtrackOffset   (rhs.theCrtrackOffset),
    theLineGsdCorr     (rhs.theLineGsdCorr),
    theSampGsdCorr     (rhs.theSampGsdCorr),
    theRollOffset      (rhs.theRollOffset),
    theYawOffset       (rhs.theYawOffset),
    theYawRate         (rhs.theYawRate),
    theMapRotation     (rhs.theMapRotation),
    theRollRotMat      (rhs.theRollRotMat)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(rhs): entering..." << std::endl;
    
    initAdjustableParameters();
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::ossimLandSatModel(rhs): returning..." << std::endl;
 }

◆ ~ossimLandSatModel()

ossimLandSatModel::~ossimLandSatModel ( )

protectedvirtual

uses file path to init model

Definition at line 247 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::~ossimLandSatModel: entering..." << std::endl;
 
    theMapProjection = 0;
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::~ossimLandSatModel: returning..." << std::endl;
 }

Member Function Documentation

◆ dup()

ossimObject * ossimLandSatModel::dup ( ) const

virtual

Implements ossimProjection.

Definition at line 239 of file ossimLandSatModel.cpp.

References ossimLandSatModel().

 {
    return new ossimLandSatModel(*this);
 }

◆ imagingRay()

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

virtual

Given an image point, returns a ray originating at some arbitrarily high point (ideally at the sensor position) and pointing towards the target.

Reimplemented from ossimSensorModel.

Definition at line 459 of file ossimLandSatModel.cpp.

 {
 #if 0
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "ossimLandSatModel::imagingRay: entering..." << std::endl;
 
    bool debug_flag = false; // setable by interactive debugger
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "inImgPt = " << inImgPt << std::endl;
    }
 #endif
    
    //***
    // Get ground point for given map image point:
    //***
    ossimDpt rot_img_pt(-inImgPt.line*theMapAzimSin+inImgPt.samp*theMapAzimCos,
                         inImgPt.line*theMapAzimCos+inImgPt.samp*theMapAzimSin);
    ossimDpt map_point 
       (theMapOffset.x + rot_img_pt.samp*(theGSD.samp+theSampGsdCorr), 
        theMapOffset.y - rot_img_pt.line*(theGSD.line+theLineGsdCorr));
 
    ossimGpt inGndPt (theMapProjection->inverse(map_point));
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t theMapOffset="<<theMapOffset<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t rot_img_pt="<<rot_img_pt<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t image point map_point="<<map_point<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t inGndPt="<<inGndPt<<endl;
    }
 #endif
    //***
    // Rotate the image map coordinates by the map-to-IC rotation
    // to arrive at intrack/crosstrack coordinate:
    //***
    ossimDpt offInMapPt (inImgPt - theRefImgPt);
    ossimDpt icInPt
       (offInMapPt.line*theMap2IcRotSin + offInMapPt.samp*theMap2IcRotCos,
        offInMapPt.line*theMap2IcRotCos - offInMapPt.samp*theMap2IcRotSin);
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t offInMapPt="<<offInMapPt<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t icInPt="<<icInPt<<endl;
    }
 #endif
    
    //***
    // Establish an image map point at vehicle NADIR corresponding to imaging
    // line:
    //***
    ossimDpt icNdrPt (0.0, icInPt.line);
    ossimDpt offNdrMapPt
       (-icNdrPt.line*theMap2IcRotSin + icNdrPt.samp*theMap2IcRotCos,
         icNdrPt.line*theMap2IcRotCos + icNdrPt.samp*theMap2IcRotSin);
    ossimDpt ndrMapPt(offNdrMapPt + theRefImgPt);
    ossimDpt rotNdrMapPt
       (-ndrMapPt.line*theMapAzimSin + ndrMapPt.samp*theMapAzimCos,
        ndrMapPt.line*theMapAzimCos + ndrMapPt.samp*theMapAzimSin);
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t icNdrPt="<<icNdrPt<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t offNdrMapPt="<<offNdrMapPt<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t ndrMapPt="<<ndrMapPt<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t rotNdrMapPt="<<rotNdrMapPt<<endl;
    }
 #endif
 
    //***
    // Establish the vehicle position at time of line imaging:
    //***
    map_point.y =theMapOffset.y+rotNdrMapPt.x*(theGSD.x+theSampGsdCorr);
    if ((theProjectionType == SOM_ORBIT) || (theProjectionType == SOM_MAP))
       map_point.x = theMapOffset.x+rotNdrMapPt.y*(theGSD.y+theLineGsdCorr);
    else
       map_point.x = theMapOffset.x-rotNdrMapPt.y*(theGSD.y+theLineGsdCorr);
       
    ossimGpt vehiclePlhPos(theMapProjection->inverse(map_point));
    vehiclePlhPos.hgt = theOrbitAltitude;
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t map_point="<<map_point<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t vehiclePlhPos="<<vehiclePlhPos<<endl;
    }
 #endif
    
    //***
    // Establish an LSR space at the vehicle with X along the intrack direction:
    //***
    ossimLsrSpace icrSpace (vehiclePlhPos, theMeridianalAngle-90.0);
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t icrSpace="<<icrSpace<<endl;
    }
 #endif
 
    //***
    // Now establish a ray pointing to the imaged ground point from the vehicle
    // in ICR space:
    //***
    ossimLsrPoint lsrInPt (inGndPt, icrSpace);
    ossimLsrPoint vehicleLsrPos (0.0, 0.0, 0.0, icrSpace);
    ossimLsrRay   initLsrImgRay (vehicleLsrPos, lsrInPt);
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t initLsrImgRay="<<initLsrImgRay<<endl;
    }
 #endif
 
    //***
    // Establish the attitude rotation matrix considering the attitude biases
    // and rates:
    //***
    double cos, sin;
    double norm_line = inImgPt.line/theImageSize.line;
    double yaw = theYawOffset + theYawRate*norm_line;
    cos = ossim::cosd(yaw);
    sin = ossim::sind(yaw);
    NEWMAT::Matrix T_yaw = ossimMatrix3x3::create( cos,-sin, 0.0,
                                                   sin, cos, 0.0,
                                                   0.0, 0.0, 1.0);
    NEWMAT::Matrix attRotMat = T_yaw * theRollRotMat;
 
    //***
    // Now apply the perturbation to the ray due to the adjustable parameters,
    // and arrive at the ECF adjusted imaging ray:
    //***
    ossimLsrVector adjLsrImgRayDir (attRotMat*initLsrImgRay.direction().data(),
                                    icrSpace);
    ossimLsrPoint  adjLsrImgRayOrg (theIntrackOffset,
                                    theCrtrackOffset,
                                    0.0,  // no radial adjustment of position
                                    icrSpace);
    ossimLsrRay adjLsrImgRay (adjLsrImgRayOrg, adjLsrImgRayDir);
    image_ray = ossimEcefRay(adjLsrImgRay);
 
 #if 0
    if (traceDebug() || debug_flag)
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t adjLsrImgRay="<<adjLsrImgRay<<endl;
       ossimNotify(ossimNotifyLevel_DEBUG) << "\t image_ray="<<image_ray<<endl;
    }
 
    if (traceExec())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::imagingRay: Returning..." << std::endl;
    }
 #endif
 }

◆ initAdjustableParameters()

void ossimLandSatModel::initAdjustableParameters ( )

protectedvirtual

Reimplemented from ossimAdjustableParameterInterface.

Definition at line 1107 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::initAdjustableParameters: entering..." << std::endl;
 
    //
    // Allocate storage for adjustables and assign their names and units strings
    //
    resizeAdjustableParameterArray(NUM_ADJUSTABLE_PARAMS);
    int numParams = getNumberOfAdjustableParameters();
    //***
    // Initialize base-class adjustable parameter array:
    //***
    for (int i=0; i<numParams; i++)
    {
       setAdjustableParameter(i, 0.0);
       setParameterDescription(i, PARAM_NAMES[i]);
       setParameterUnit(i,PARAM_UNITS[i]);
    }
    //***
    // Initialize base-class parameter sigma array:
    //***
    setParameterSigma(INTRACK_OFFSET, 500.0); //change for Landsat 5
    setParameterSigma(CRTRACK_OFFSET, 500.0); //change for Landsat 5
    setParameterSigma(LINE_GSD_CORR, 0.005);  
    setParameterSigma(SAMP_GSD_CORR, 0.005);  
    setParameterSigma(ROLL_OFFSET, 0.01);  
    setParameterSigma(YAW_OFFSET, 0.01);  
    setParameterSigma(YAW_RATE, 0.05);  
    setParameterSigma(MAP_ROTATION, 0.1);
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::initAdjustableParameters: returning..." << std::endl;
 }

◆ initFromHeader()

void ossimLandSatModel::initFromHeader ( const ossimFfL7 & head )

protected

Initializes the model given a fast format header.

Definition at line 260 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::initFromHeader: entering..." << std::endl;
    
    //***
    // Initialize base-class data members:
    //***
    theRefGndPt         = head.theCenterGP;
    theRefImgPt         = head.theCenterImagePoint;
    theImageSize.x      = head.thePixelsPerLine;
    theImageSize.y      = head.theLinesPerBand;
    theImageID          = head.theRequestNumber;
    theImageClipRect    = ossimDrect(0, 0, theImageSize.x-1, theImageSize.y-1);
    theGSD.samp         = head.theGsd;
    theGSD.line         = head.theGsd;
    theIllumAzimuth     = head.theSunAzimuth;
    theIllumElevation   = head.theSunElevation;
    theMapZone          = head.theUsgsMapZone;
    theWrsPathNumber    = head.thePathNumber;
    theWrsRowNumber     = head.theRowNumber;
    theRollOffset       = head.theOffNadirAngle;
    theMeanGSD          = head.theGsd;
    
    //satellite orbit and accuracy
    ossimString satname(head.theSatName);
    if (satname.contains("7"))
    {
       theOrbitAltitude    = L7_ORBIT_ALTITUDE;
       theOrbitInclination = L7_ORBIT_INCLINATION;
       theNominalPosError  = L7_NOMINAL_POS_ERROR;
    } else if (satname.contains("5"))
    {
       theOrbitAltitude    = L5_ORBIT_ALTITUDE;
       theOrbitInclination = L5_ORBIT_INCLINATION;
       theNominalPosError  = L5_NOMINAL_POS_ERROR;
    } else {
       theErrorStatus = 1; //MODEL_ERROR
       if(traceDebug())
       {
          ossimNotify(ossimNotifyLevel_WARN) << "WARNING ossimLandSatModel::initFromHeader: " << "Unknown satellite name : " << satname << std::endl;
       }
    }   
    
    
    //***
    // Assign the meridianal angle. The Meridianal angle is
    // the angle between true north and the vehicles ground track. The
    // orientation angle is the rotation of the raster image from north (this
    // will be 0 for map-oriented imagery).
    //
    // Compute meridianal angle between the along-track direction
    // and the meridian going through the scene center excluding earth rotation
    // effects. The spherical triangle equation is Napier's rule for right-angle
    // spherical triangles.  This angle is a negative number from -8.2 to -90
    // depending on scene latitude:
    //***
    double phi_c = ossim::atand(ossim::tand(theRefGndPt.lat)/
                                GEODETIC_2_GEOCENTRIC_FACTOR);
    double cos_phi_c   = ossim::cosd(phi_c);
    theMeridianalAngle = -ossim::asind(ossim::cosd(theOrbitInclination) / cos_phi_c);
    theMapAzimAngle = head.theOrientationAngle;
 
    //***
    // Fetch the corner points from the header:
    //***
    ossimDpt v[4]; // temporarily holds vertices for ground polygon
    v[0] = head.theUL_Corner;
    v[1] = head.theUR_Corner;
    v[2] = head.theLR_Corner;
    v[3] = head.theLL_Corner;
    theBoundGndPolygon = ossimPolygon(4, v);
 
    //***
    // Initialize map projection parameters:
    // Determine the projection/orientation type:
    //***
    ossimString orient_type = head.theProductType;
    ossimString proj_type   = head.theMapProjectionName;
    if (proj_type.contains("SOM"))
    {
       //***
       // Inverse the line direction GSD for SOM, since projection "intrack"
       // axis (x) is descending:
       //***
       if (orient_type.contains("ORBIT"))
          theProjectionType  = SOM_ORBIT;
       else
          theProjectionType = SOM_MAP;
 
       // this is a hack.  After testing two it is consistently off
       // by 90 degrees.
       //
       theMapAzimAngle -= 90.0;
    }
    else if (proj_type.contains("UTM"))
    {
       if (orient_type.contains("ORBIT"))
          theProjectionType  = UTM_ORBIT;
       else
          theProjectionType = UTM_MAP;
    }
    else
    {
       theErrorStatus = 1; //MODEL_ERROR
       if(traceDebug())
       {
          ossimNotify(ossimNotifyLevel_WARN) << "WARNING ossimLandSatModel::initFromHeader: "
                                              << "Unknown projection/orientation type." << std::endl;
       }
       return;
    }
 
    //***
    // Establish the map projection:
    //***
    initMapProjection();
 
       
    theMapOffset = theMapProjection->forward(head.theUL_Corner);
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG)
          << "DEBUG ossimLandSatModel::initFromHeader:"
          << "\ntheMapProjection:\n";
 
       theMapProjection->print(ossimNotify(ossimNotifyLevel_DEBUG));
 
       ossimNotify(ossimNotifyLevel_DEBUG)
          << "\nHeader upper left ground point:  " << head.theUL_Corner
          << std::endl;
    }
 
    //*** 
    // initialize remaining data members:
    //***
    initAdjustableParameters();
    updateModel();
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::initFromHeader: returning..." << std::endl;
 }

◆ initMapProjection()

void ossimLandSatModel::initMapProjection ( )

protected

Definition at line 1012 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) <<  "DEBUG ossimLandSatModel::writeGeomTemplate: entering... " << std::endl;
 
    theMapProjection = 0;   
    //*** 
    // Instantiate the proper map projection:
    //***
    if ((theProjectionType == SOM_ORBIT) || (theProjectionType == SOM_MAP))
    {
       theMapProjection = new ossimSpaceObliqueMercatorProjection(
          ossimSpaceObliqueMercatorProjection::SOM_TYPE_LANDSAT_7,
          (double)theWrsPathNumber);
    }
    else 
    {
       ossimUtmProjection* utm   = new ossimUtmProjection(theMapZone);
       theMapProjection = utm;
       if(theRefGndPt.latd() < 0.0)
       {
          utm->setHemisphere('S');
       }
       else
       {
          utm->setHemisphere('N');
       }
    }
 
    //***
    // Initialize angle sines/cosines used in phiLambda projection:
    //***
    theMap2IcRotAngle = theMeridianalAngle + theMapAzimAngle;
    theMap2IcRotCos   = ossim::cosd(theMap2IcRotAngle);
    theMap2IcRotSin   = ossim::sind(theMap2IcRotAngle);
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) <<  "DEBUG ossimLandSatModel::writeGeomTemplate: returning... " << std::endl;
    return;
 }

◆ lineSampleHeightToWorld()

void ossimLandSatModel::lineSampleHeightToWorld	(	const ossimDpt &	lineSampPt,
		const double &	heightEllipsoid,
		ossimGpt &	worldPt
	)		const

virtual

METHOD: lineSampleHeightToWorld This is the pure virtual that performs the actual work of projecting the image point to the given elevation above Ellipsoid.

Implements ossimSensorModel.

Definition at line 409 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::lineSampleHeightToWorld: entering..." << std::endl;
 
 #if 0
    //***
    // Extrapolate if point is outside image:
    //***
    if (!insideImage(image_point))
    {
       gpt = extrapolate(image_point, height);
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::lineSampleHeightToWorld: returning..." << std::endl;
       return;
    }
 #endif
    
    //***
    // First establish imaging ray from image point:
    //***
    ossimEcefRay imaging_ray;
    imagingRay(image_point, imaging_ray);
    ossimEcefPoint Pecf (imaging_ray.intersectAboveEarthEllipsoid(height));
    gpt = ossimGpt(Pecf);
 }

◆ loadState()

bool ossimLandSatModel::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 ossimSensorModel.

Definition at line 718 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::loadState: entering..." << std::endl;
 
    if (traceDebug())
    {
       ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::loadState:"
                                           << "\nInput kwl:  " << kwl
                                           << std::endl;
    }
 
    const char* value = NULL;
    const char* keyword =NULL;
    bool success;
 
    //***
    // Assure this keywordlist contains correct type info:
    //***
    value = kwl.find(prefix, ossimKeywordNames::TYPE_KW);
    if (!value || (strcmp(value, TYPE_NAME(this)))) 
      {
        theErrorStatus = 1;
        return false;
      }
 
    //
    // Clear out any existing adjustable params:
    //
    if(getNumberOfAdjustableParameters() != NUM_ADJUSTABLE_PARAMS)
      {
        initAdjustableParameters();
      }
    
    //***
    // Pass on to the base-class for parsing first:
    //***
    success = ossimSensorModel::loadState(kwl, prefix);
    if (!success) 
      {
        theErrorStatus++;
        return false;
      }
  
    keyword = PROJECTION_TYPE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
 
    theProjectionType = (ProjectionType) atoi(value);
  
    keyword = MAP_ZONE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMapZone = atoi(value);
  
    keyword = MAP_OFFSET_X_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMapOffset.x = atof(value);
  
    keyword = MAP_OFFSET_Y_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMapOffset.y = atof(value);
  
    keyword = WRS_PATH_NUMBER_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theWrsPathNumber = atoi(value);
 
    keyword = ROW_NUMBER_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theWrsRowNumber = atoi(value);
  
    keyword = ILLUM_AZIMUTH_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theIllumAzimuth = atof(value);
  
    keyword = ILLUM_ELEVATION_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theIllumElevation = atof(value);
  
    keyword = MERIDIANAL_ANGLE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMeridianalAngle = atof(value);
  
    keyword = ORBIT_ALTITUDE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theOrbitAltitude = atof(value);
    
    keyword = ORBIT_INCLINATION_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theOrbitInclination = atof(value);
    
    keyword = MAP_AZIM_ANGLE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMapAzimAngle = atof(value);
  
    keyword = MAP_2Ic_ROT_ANGLE_KW;
    value = kwl.find(prefix, keyword);
    if (!value)
      {
        theErrorStatus++;
        return false;
      }
    theMap2IcRotAngle = atof(value);
 
    theIntrackOffset = 0.0;
    theCrtrackOffset = 0.0;
    theLineGsdCorr   = 0.0;
    theSampGsdCorr   = 0.0;
    theRollOffset    = 0.0;
    theYawOffset     = 0.0;
    theYawRate       = 0.0;
    theMapRotation   = 0.0;
 
    //***
    // Adjustable parameters are optional keywords:
    //***
 //    keyword = INTRACK_OFFSET_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theIntrackOffset = ossimString(value).toDouble();
 //    else
 //     theIntrackOffset = 0.0;
    
 //    keyword = CRTRACK_OFFSET_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theCrtrackOffset = ossimString(value).toDouble();
 //    else
 //      theCrtrackOffset = 0.0;
 
 //    keyword = LINE_GSD_CORR_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theLineGsdCorr = ossimString(value).toDouble();
 //    else
 //       theLineGsdCorr = 0.0;
 
 //    keyword = SAMP_GSD_CORR_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theSampGsdCorr = ossimString(value).toDouble();
 //    else
 //       theSampGsdCorr = 0.0;
    
 //    keyword = ROLL_OFFSET_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theRollOffset = ossimString(value).toDouble();
 //    else
 //       theRollOffset = 0.0;
    
 //    keyword = YAW_OFFSET_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theYawOffset = ossimString(value).toDouble();
 //    else
 //      theYawOffset = 0.0;
    
 //    keyword = YAW_RATE_KW;
 //    value = kwl.find(prefix, keyword);
 //    if(value)
 //       theYawRate = ossimString(value).toDouble();
 //    else
 //       theYawRate = 0.0;
 
 //    keyword = MAP_ROTATION_KW;
 //    value = kwl.find(prefix,  keyword);
 //    if(value)
 //       theMapRotation = ossimString(value).toDouble();
 //    else
 //       theMapRotation = 0.0;
 
    //***
    // Initialize given parameters read:
    //***
    initMapProjection();
    updateModel();
    
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::loadState: returning..." << std::endl;
    return true;
 }

◆ print()

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

virtual

Extends base-class implementation. Dumps contents of object to ostream.

Reimplemented from ossimSensorModel.

Definition at line 626 of file ossimLandSatModel.cpp.

References ossimString::chars(), ossimDpt::line, ossimDpt::samp, ossimSensorModel::theGSD, ossimSensorModel::theImageID, ossimSensorModel::theImageSize, ossimSensorModel::theRefGndPt, and ossimSensorModel::theRefImgPt.

 {
    os << "\nDump of ossimLandSatModel object at "
       << hex << this << ":\n"
       << "\nLandSatModel -- Dump of all data members: "
       << "\n         theImageID: " << theImageID.chars()
       << "\n       theImageSize: " << theImageSize
       << "\n        theRefImgPt: " << theRefImgPt
       << "\n        theRefGndPt: " << theRefGndPt
       << "\n        theGSD.line: " << theGSD.line
       << "\n        theGSD.samp: " << theGSD.samp
       << "\n  theProjectionType: " << PROJ_TYPE[theProjectionType]
       << "\n         theMapZone: " << theMapZone
       << "\n       theMapOffset: " << theMapOffset
       << "\n   theWrsPathNumber: " << theWrsPathNumber
       << "\n    theWrsRowNumber: " << theWrsRowNumber
       << "\n    theIllumAzimuth: " << theIllumAzimuth
       << "\n  theIllumElevation: " << theIllumElevation
       << "\n   thePositionError: " << thePositionError
       << "\n theMeridianalAngle: " << theMeridianalAngle
       << "\n   theOrbitAltitude: " << theOrbitAltitude
       << "\ntheOrbitInclination: " << theOrbitInclination
       << "\n    theMapAzimAngle: " << theMapAzimAngle
       << "\n  theMap2IcRotAngle: " << theMap2IcRotAngle
       << "\n   theIntrackOffset: " << theIntrackOffset
       << "\n   theCrtrackOffset: " << theCrtrackOffset
       << "\n     theLineGsdCorr: " << theLineGsdCorr
       << "\n     theSampGsdCorr: " << theSampGsdCorr
       << "\n      theRollOffset: " << theRollOffset
       << "\n       theYawOffset: " << theYawOffset
       << "\n         theYawRate: " << theYawRate
       << "\n     theMapRotation: " << theMapRotation
       << endl;
 
    return ossimSensorModel::print(os);
 }

◆ saveState()

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

virtual

Fulfills ossimObject base-class pure virtuals. Loads and saves geometry KWL files. Returns true if successful.

Reimplemented from ossimSensorModel.

Definition at line 669 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::saveState: entering..." << std::endl;
 
    kwl.add(prefix, ossimKeywordNames::TYPE_KW, TYPE_NAME(this));
 
    //***
    // Hand off to base class for common stuff:
    //***
    ossimSensorModel::saveState(kwl, prefix);
 
    //***
    // Save off data members:
    //***
    kwl.add(prefix, PROJECTION_TYPE_KW,   theProjectionType, true);
    kwl.add(prefix, MAP_ZONE_KW,          theMapZone, true);
    kwl.add(prefix, MAP_OFFSET_X_KW,      theMapOffset.x, true);
    kwl.add(prefix, MAP_OFFSET_Y_KW,      theMapOffset.y, true);
    kwl.add(prefix, WRS_PATH_NUMBER_KW,   theWrsPathNumber, true);
    kwl.add(prefix, ROW_NUMBER_KW,        theWrsRowNumber, true);
    kwl.add(prefix, ILLUM_AZIMUTH_KW,     theIllumAzimuth, true);
    kwl.add(prefix, ILLUM_ELEVATION_KW,   theIllumElevation, true);
    kwl.add(prefix, MERIDIANAL_ANGLE_KW,  theMeridianalAngle, true);
    kwl.add(prefix, ORBIT_ALTITUDE_KW,    theOrbitAltitude, true);
    kwl.add(prefix, ORBIT_INCLINATION_KW, theOrbitInclination, true);
    kwl.add(prefix, MAP_AZIM_ANGLE_KW,    theMapAzimAngle, true);
    kwl.add(prefix, MAP_2Ic_ROT_ANGLE_KW, theMap2IcRotAngle, true);
 
 //    kwl.add(prefix, INTRACK_OFFSET_KW ,   theIntrackOffset, true);
 //    kwl.add(prefix, CRTRACK_OFFSET_KW,    theCrtrackOffset, true);
 //    kwl.add(prefix, LINE_GSD_CORR_KW,     theLineGsdCorr, true);
 //    kwl.add(prefix, SAMP_GSD_CORR_KW,     theSampGsdCorr, true);
 //    kwl.add(prefix, ROLL_OFFSET_KW,       theRollOffset, true);
 //    kwl.add(prefix, YAW_OFFSET_KW,        theYawOffset, true);
 //    kwl.add(prefix, YAW_RATE_KW,          theYawRate, true);
 //    kwl.add(prefix, MAP_ROTATION_KW,      theMapRotation, true);
    
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::saveState: returning..." << std::endl;
    return true;
 }

◆ setupOptimizer()

bool ossimLandSatModel::setupOptimizer ( const ossimString & init_file )

virtual

image to ground faster

Reimplemented from ossimOptimizableProjection.

Definition at line 1140 of file ossimLandSatModel.cpp.

 {
    //init model using file path
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::setupOptimizer(init_file): entering..." << std::endl;
 
    theMapProjection = 0;
    theIntrackOffset = 0.0;
    theCrtrackOffset = 0.0;
    theLineGsdCorr   = 0.0;   
    theSampGsdCorr   = 0.0;
    theRollOffset    = 0.0;
    theYawOffset     = 0.0;
    theYawRate       = 0.0;
    theMapRotation   = 0.0;
 
    ossimRefPtr<ossimFfL7> ff_headerp;
    if (ossimString::downcase(init_file).contains("header.dat"))
    {
       ossimRefPtr<ossimFfL5> h = new ossimFfL5(init_file); 
       ff_headerp = h.get();
       
       if (!ff_headerp->getErrorStatus())
       {
          double d = fabs(h->revb()->theUlEasting - h->revb()->theCenterEasting)/h->theGsd;
          h->theCenterImagePoint.x = static_cast<ossim_int32>(d); // d + 0.5 ???
          
          d = fabs(h->revb()->theUlNorthing - h->revb()->theCenterNorthing)/h->theGsd;
          h->theCenterImagePoint.y = static_cast<ossim_int32>(d); // d + 0.5 ???
          initFromHeader(*ff_headerp);
          
          theMapOffset.x = h->revb()->theUlEasting;
          theMapOffset.y = h->revb()->theUlNorthing;
          
       }
       else
       {
          ff_headerp = 0;
       }
    }
    else
    {
       ff_headerp=new ossimFfL7(init_file);
       if (!ff_headerp->getErrorStatus())
       {
          initFromHeader(*ff_headerp);
       }
       else
       {
          ff_headerp = 0;
       }
    }
    if(!ff_headerp.valid())
    {
       //
       // If not header, then check for possible KWL file. The loadState sets the
       // error status:
       //
       ossimFilename init_filename(init_file);
       ossimKeywordlist kwl(init_filename);
       loadState(kwl);
    }
 
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::setupOptimizer(init_file): Exited..." << std::endl;
    return true;
 }

◆ updateModel()

void ossimLandSatModel::updateModel ( )

virtual

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 1054 of file ossimLandSatModel.cpp.

References ossimAdjustableParameterInterface::computeParameterOffset(), ossim::cosd(), ossimMatrix3x3::create(), CRTRACK_OFFSET, INTRACK_OFFSET, LINE_GSD_CORR, MAP_ROTATION, ROLL_OFFSET, SAMP_GSD_CORR, ossim::sind(), theCrtrackOffset, theIntrackOffset, theLineGsdCorr, theMapAzimAngle, theMapAzimCos, theMapAzimSin, theMapRotation, theProjectionType, theRollOffset, theRollRotMat, theSampGsdCorr, theYawOffset, theYawRate, UTM_ORBIT, YAW_OFFSET, and YAW_RATE.

 {
    //***
    // Update the geometry...
    // Adjusted = Initial + adj_parm[i]*adj_sigma[i]
    // Some parameters have a zero initial value and that term is not included.
    //***
   theIntrackOffset = computeParameterOffset(INTRACK_OFFSET);
   theCrtrackOffset = computeParameterOffset(CRTRACK_OFFSET);
   theLineGsdCorr   = computeParameterOffset(LINE_GSD_CORR);
   theSampGsdCorr   = computeParameterOffset(SAMP_GSD_CORR);
   theRollOffset    = computeParameterOffset(ROLL_OFFSET);
   theYawOffset     = computeParameterOffset(YAW_OFFSET);
   theYawRate       = computeParameterOffset(YAW_RATE);
   theMapRotation   = computeParameterOffset(MAP_ROTATION);
 
 //   int numParams = getNumberOfAdjustableParameters();
    //***
    // Now the initial values have been updated, zero out the adjustment:
    //***
 //    for (int i=0; i<numParams; i++)
 //       setAdjustableParameter(i, 0.0);
 
    //***
    // Compute image-common values dependent on adjustable parameters:
    //***
    if (theProjectionType == UTM_ORBIT)
    {
       theMapAzimCos = ossim::cosd(-theMapAzimAngle + theMapRotation);
       theMapAzimSin = ossim::sind(-theMapAzimAngle + theMapRotation);
    }
    else
    {
       theMapAzimCos = ossim::cosd(theMapAzimAngle + theMapRotation);
       theMapAzimSin = ossim::sind(theMapAzimAngle + theMapRotation);
    }
 
    double cos = ossim::cosd(theRollOffset);
    double sin = ossim::sind(theRollOffset);
    theRollRotMat = ossimMatrix3x3::create( 1.0, 0.0, 0.0,
                                            0.0, cos,-sin,
                                            0.0, sin, cos);
    
    
 }

◆ useForward()

virtual bool ossimLandSatModel::useForward ( ) const

inlinevirtual

ossimOptimizableProjection

Implements ossimOptimizableProjection.

Definition at line 120 of file ossimLandSatModel.h.

120 {return false;}

◆ writeGeomTemplate()

void ossimLandSatModel::writeGeomTemplate ( ostream & os )

static

Writes a template of geom keywords processed by loadState and saveState to output stream.

Definition at line 964 of file ossimLandSatModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::writeGeomTemplate: entering..." << std::endl;
 
    os <<
       "//**************************************************************\n"
       "// Template for LandSat model keywordlist\n"
       "//**************************************************************\n"
       << ossimKeywordNames::TYPE_KW << ": " << "ossimLandSatModel" << endl;
 
    ossimSensorModel::writeGeomTemplate(os);
    
    os << "//\n"
       << "// Derived-class ossimLandSatModel Keywords:\n"
       << "//\n"
       << PROJECTION_TYPE_KW     << ": <float>\n"
       << MAP_ZONE_KW            << ": <float>\n"
       << MAP_OFFSET_X_KW        << ": <float>\n"
       << MAP_OFFSET_Y_KW        << ": <float>\n"
       << WRS_PATH_NUMBER_KW     << ": <float>\n"
       << ROW_NUMBER_KW          << ": <float>\n"
       << ILLUM_AZIMUTH_KW       << ": <float>\n"
       << ILLUM_ELEVATION_KW     << ": <float>\n"
       << MERIDIANAL_ANGLE_KW    << ": <float>\n"
       << ORBIT_ALTITUDE_KW      << ": <float>\n"
       << ORBIT_INCLINATION_KW   << ": <float>\n"
       << MAP_AZIM_ANGLE_KW      << ": <float>\n"
       << MAP_2Ic_ROT_ANGLE_KW   << ": <float>\n"
       << INTRACK_OFFSET_KW      << ": <float> [optional]\n"
       << CRTRACK_OFFSET_KW      << ": <float> [optional]\n"
       << LINE_GSD_CORR_KW       << ": <float> [optional]\n"
       << SAMP_GSD_CORR_KW       << ": <float> [optional]\n"
       << ROLL_OFFSET_KW         << ": <float> [optional]\n"
       << YAW_OFFSET_KW          << ": <float> [optional]\n"
       << YAW_RATE_KW            << ": <float> [optional]\n"
       << MAP_ROTATION_KW        << ": <float> [optional]\n"
       << endl;
    os << "\n" <<endl;
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimLandSatModel::writeGeomTemplate: returning..." << std::endl;
    return;
 }

Member Data Documentation

◆ theCrtrackOffset

double ossimLandSatModel::theCrtrackOffset

protected

Definition at line 162 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theIllumAzimuth

double ossimLandSatModel::theIllumAzimuth

protected

Definition at line 137 of file ossimLandSatModel.h.

◆ theIllumElevation

double ossimLandSatModel::theIllumElevation

protected

Definition at line 138 of file ossimLandSatModel.h.

◆ theIntrackOffset

double ossimLandSatModel::theIntrackOffset

protected

Definition at line 161 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theLineGsdCorr

double ossimLandSatModel::theLineGsdCorr

protected

Definition at line 163 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theMap2IcRotAngle

double ossimLandSatModel::theMap2IcRotAngle

protected

Definition at line 154 of file ossimLandSatModel.h.

◆ theMap2IcRotCos

double ossimLandSatModel::theMap2IcRotCos

protected

Definition at line 155 of file ossimLandSatModel.h.

◆ theMap2IcRotSin

double ossimLandSatModel::theMap2IcRotSin

protected

Definition at line 156 of file ossimLandSatModel.h.

◆ theMapAzimAngle

double ossimLandSatModel::theMapAzimAngle

protected

Definition at line 151 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theMapAzimCos

double ossimLandSatModel::theMapAzimCos

protected

Definition at line 152 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theMapAzimSin

double ossimLandSatModel::theMapAzimSin

protected

Definition at line 153 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theMapOffset

ossimDpt ossimLandSatModel::theMapOffset

protected

Definition at line 142 of file ossimLandSatModel.h.

◆ theMapProjection

ossimRefPtr<ossimMapProjection> ossimLandSatModel::theMapProjection

protected

Definition at line 149 of file ossimLandSatModel.h.

◆ theMapRotation

double ossimLandSatModel::theMapRotation

protected

Definition at line 168 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theMapZone

int ossimLandSatModel::theMapZone

protected

Definition at line 141 of file ossimLandSatModel.h.

◆ theMeridianalAngle

double ossimLandSatModel::theMeridianalAngle

protected

Definition at line 145 of file ossimLandSatModel.h.

◆ theOrbitAltitude

double ossimLandSatModel::theOrbitAltitude

protected

Definition at line 139 of file ossimLandSatModel.h.

◆ theOrbitInclination

double ossimLandSatModel::theOrbitInclination

protected

Definition at line 140 of file ossimLandSatModel.h.

◆ thePositionError

double ossimLandSatModel::thePositionError

protected

Definition at line 146 of file ossimLandSatModel.h.

◆ theProjectionType

ProjectionType ossimLandSatModel::theProjectionType

protected

Definition at line 148 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theRollOffset

double ossimLandSatModel::theRollOffset

protected

Definition at line 165 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theRollRotMat

NEWMAT::Matrix ossimLandSatModel::theRollRotMat

protected

Definition at line 173 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theSampGsdCorr

double ossimLandSatModel::theSampGsdCorr

protected

Definition at line 164 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theWrsPathNumber

int ossimLandSatModel::theWrsPathNumber

protected

Definition at line 143 of file ossimLandSatModel.h.

◆ theWrsRowNumber

int ossimLandSatModel::theWrsRowNumber

protected

Definition at line 144 of file ossimLandSatModel.h.

◆ theYawOffset

double ossimLandSatModel::theYawOffset

protected

Definition at line 166 of file ossimLandSatModel.h.

Referenced by updateModel().

◆ theYawRate

double ossimLandSatModel::theYawRate

protected

Definition at line 167 of file ossimLandSatModel.h.

Referenced by updateModel().

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ AdjustParamIndex

◆ ProjectionType

Constructor & Destructor Documentation

◆ ossimLandSatModel() [1/5]

◆ ossimLandSatModel() [2/5]

◆ ossimLandSatModel() [3/5]

◆ ossimLandSatModel() [4/5]

◆ ossimLandSatModel() [5/5]

◆ ~ossimLandSatModel()

Member Function Documentation

◆ dup()

◆ imagingRay()

◆ initAdjustableParameters()

◆ initFromHeader()

◆ initMapProjection()

◆ lineSampleHeightToWorld()

◆ loadState()

◆ print()

◆ saveState()

◆ setupOptimizer()

◆ updateModel()

◆ useForward()

◆ writeGeomTemplate()

Member Data Documentation

◆ theCrtrackOffset

◆ theIllumAzimuth

◆ theIllumElevation

◆ theIntrackOffset

◆ theLineGsdCorr

◆ theMap2IcRotAngle

◆ theMap2IcRotCos

◆ theMap2IcRotSin

◆ theMapAzimAngle

◆ theMapAzimCos

◆ theMapAzimSin

◆ theMapOffset

◆ theMapProjection

◆ theMapRotation

◆ theMapZone

◆ theMeridianalAngle

◆ theOrbitAltitude

◆ theOrbitInclination

◆ thePositionError

◆ theProjectionType

◆ theRollOffset

◆ theRollRotMat

◆ theSampGsdCorr

◆ theWrsPathNumber

◆ theWrsRowNumber

◆ theYawOffset

◆ theYawRate