#include <ossimCoarseGridModel.h>

Inheritance diagram for ossimCoarseGridModel:

Public Member Functions
	ossimCoarseGridModel ()

	ossimCoarseGridModel (const ossimCoarseGridModel &copy_this)

	ossimCoarseGridModel (const ossimFilename &geom_file)
	Accepts name of geometry file. More...

	ossimCoarseGridModel (const ossimKeywordlist &geom_kwl)
	Accepts OSSIM keywordlist geometry file. More...

	~ossimCoarseGridModel ()

virtual void	buildGrid (const ossimDrect &imageBounds, ossimProjection *proj, double heightDelta=500.0, bool enableHeightFlag=false, bool makeAdjustableFlag=true)
	This method will build a grid from any projector. More...

virtual void	buildGrid (const ossimDrect &imageBounds, ossimImageGeometry *geom, double heightDelta=500.0, bool enableHeightFlag=false, bool makeAdjustableFlag=true)
	Assigns the grid data given a geometry. More...

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

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

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

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

bool	saveCoarseGrid (const ossimFilename &cgFileName) const
	Saves the coarse grid to the specified file. More...

bool	loadCoarseGrid (const ossimFilename &cgFileName)
	Loads the coarse grid from the specified file. More...

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

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

virtual void	lineSampleHeightToWorld (const ossimDpt &image_point, const double &heightEllipsoid, ossimGpt &world_pt) const
	This is the virtual that performs the actual work of projecting the image point to the earth at some specified elevation. More...

virtual void	initAdjustableParameters ()

virtual bool	useForward () const

virtual bool	isAffectedByElevation () const
	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	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 void	updateModel ()

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

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

Static Public Member Functions
static void	setInterpolationError (double error=.1)
	This is used when building a grid from a projector. More...

static void	setMinGridSpacing (ossim_int32 minSpacing=100)

static void	writeGeomTemplate (ostream &os)
	Writes a template of geometry keywords processed by loadState and saveState to output stream. More...

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

Protected Member Functions
void	reallocateGrid (const ossimIpt &size)
	Deletes existing allocated memory and reallocates new space. More...

void	initializeModelParams (ossimIrect irect)
	Initializes base class data members after grids have been assigned. More...

virtual ossimGpt	extrapolate (const ossimDpt &imgPt, const double &height=ossim::nan()) const
	Implements its own extrapolation since this can be handled by ossimDblGrid. More...

Protected Member Functions inherited from ossimSensorModel
virtual	~ossimSensorModel ()

virtual ossimDpt	extrapolate (const ossimGpt &gp) 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
ossimFilename	theGridFilename

ossimDblGrid	theLatGrid

ossimDblGrid	theLonGrid

ossimDblGrid	theDlatDhGrid

ossimDblGrid	theDlonDhGrid

ossimDblGrid *	theDlatDparamGrid

ossimDblGrid *	theDlonDparamGrid

ossimAdjustmentInfo	theInitialAdjustment

bool	theHeightEnabledFlag

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

Static Protected Attributes
static double	theInterpolationError = .1

static ossim_int32	theMinGridSpacing = 100

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 }

Detailed Description

Definition at line 34 of file ossimCoarseGridModel.h.

Constructor & Destructor Documentation

◆ ossimCoarseGridModel() [1/4]

ossimCoarseGridModel::ossimCoarseGridModel ( )

Definition at line 58 of file ossimCoarseGridModel.cpp.

    :
       ossimSensorModel(),
       theDlatDparamGrid (0),
       theDlonDparamGrid (0),
       theHeightEnabledFlag(true)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel: entering..."
       << std::endl;
 
    theLatGrid.setDomainType(ossimDblGrid::SAWTOOTH_90);
    theLonGrid.setDomainType(ossimDblGrid::WRAP_180);
    theLatGrid.enableExtrapolation();
    theLonGrid.enableExtrapolation();
 
    setErrorStatus();
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel: returning..."
       << std::endl;
 }

◆ ossimCoarseGridModel() [2/4]

ossimCoarseGridModel::ossimCoarseGridModel ( const ossimCoarseGridModel & copy_this )

Definition at line 85 of file ossimCoarseGridModel.cpp.

    :
       ossimSensorModel  (model),
       theGridFilename   (model.theGridFilename),
       theLatGrid        (model.theLatGrid),
       theLonGrid        (model.theLonGrid),
       theDlatDhGrid     (model.theDlatDhGrid),
       theDlonDhGrid     (model.theDlonDhGrid),
       theDlatDparamGrid (0),
       theDlonDparamGrid (0),
       theHeightEnabledFlag(true)
 {
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel(model): entering..."
       << std::endl;
 
    int numberOfParams = getNumberOfAdjustableParameters();
    if(numberOfParams)
    {
       //***
       // Allocate adjustable parameter partials grids then assign:
       //***
       theDlatDparamGrid = new ossimDblGrid [numberOfParams];
       theDlonDparamGrid = new ossimDblGrid [numberOfParams];
       
       for (int i=0; i<numberOfParams; i++)
       {
          theDlatDparamGrid[i] = model.theDlatDparamGrid[i];
          theDlonDparamGrid[i] = model.theDlonDparamGrid[i];
       }
    }
    
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel: returning..."
       << std::endl;
 }

◆ ossimCoarseGridModel() [3/4]

ossimCoarseGridModel::ossimCoarseGridModel ( const ossimFilename & geom_file )

Accepts name of geometry file.

This can be either MET ECG geom file, or OSSIM keywordlist geometry file.

Definition at line 128 of file ossimCoarseGridModel.cpp.

    :
       ossimSensorModel(),
       theDlatDparamGrid (0),
       theDlonDparamGrid (0),
       theHeightEnabledFlag(true)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel(geom_file): entering..." << std::endl;
    
    theLatGrid.setDomainType(ossimDblGrid::SAWTOOTH_90);
    theLonGrid.setDomainType(ossimDblGrid::WRAP_180);
    theLatGrid.enableExtrapolation();
    theLonGrid.enableExtrapolation();
 
    ossimKeywordlist kwl;
    if(geom_file.exists()&&kwl.addFile(geom_file))
    {
       loadState(kwl);
       theGridFilename = geom_file.path();
    }
    else
    {
       ++theErrorStatus;
    }
 
    if (traceExec())   ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel(geom_file): returning..." << std::endl;
    return;
 }

◆ ossimCoarseGridModel() [4/4]

ossimCoarseGridModel::ossimCoarseGridModel ( const ossimKeywordlist & geom_kwl )

Accepts OSSIM keywordlist geometry file.

Definition at line 163 of file ossimCoarseGridModel.cpp.

    :
       ossimSensorModel(),
       theDlatDparamGrid (0),
       theDlonDparamGrid (0),
       theHeightEnabledFlag(true)
 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::ossimCoarseGridModel(geom_kwl): entering..." << std::endl;
 
    theLatGrid.setDomainType(ossimDblGrid::SAWTOOTH_90);
    theLonGrid.setDomainType(ossimDblGrid::WRAP_180);
    theLatGrid.enableExtrapolation();
    theLonGrid.enableExtrapolation();
 
   // Parse keywordlist for geometry:
    loadState(geom_kwl);
 }

◆ ~ossimCoarseGridModel()

ossimCoarseGridModel::~ossimCoarseGridModel ( )

Definition at line 459 of file ossimCoarseGridModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::~ossimCoarseGridModel: entering..."
       << std::endl;
 
    if(theDlatDparamGrid&&theDlonDparamGrid)
    {
       //***
       // Deallocate memory:
       //***
       delete [] theDlatDparamGrid;
       delete [] theDlonDparamGrid;
       theDlatDparamGrid = NULL;
       theDlonDparamGrid = NULL;
    }
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG)
       << "DEBUG ossimCoarseGridModel::~ossimCoarseGridModel: returning..."
       << std::endl;
 }

Member Function Documentation

◆ buildGrid() [1/2]

void ossimCoarseGridModel::buildGrid	(	const ossimDrect &	imageBounds,
		ossimProjection *	proj,
		double	heightDelta = `500.0`,
		bool	enableHeightFlag = `false`,
		bool	makeAdjustableFlag = `true`
	)

virtual

This method will build a grid from any projector.

Assigns the grid data given a projection (typically a rigorous sensor model)

The accuracy of the grid can be controlled by the static method setInterpolationError().

Parameters

imageBounds	Must be the image space bounds for the projection.
proj	The projector that will be used to approximate a bilinear grid over.

Definition at line 184 of file ossimCoarseGridModel.cpp.

References ossimRefPtr< T >::get(), and ossimImageGeometry::setProjection().

 {
    ossimRefPtr<ossimImageGeometry> geom = new ossimImageGeometry();
    geom->setProjection(proj);
    buildGrid(imageBounds, geom.get(), heightDelta, enableHeightFlag, makeAdjustableFlag);
 }  

◆ buildGrid() [2/2]

void ossimCoarseGridModel::buildGrid	(	const ossimDrect &	imageBounds,
		ossimImageGeometry *	geom,
		double	heightDelta = `500.0`,
		bool	enableHeightFlag = `false`,
		bool	makeAdjustableFlag = `true`
	)

virtual

Assigns the grid data given a geometry.

Definition at line 199 of file ossimCoarseGridModel.cpp.

References ossimImageGeometry::getProjection(), ossimDrect::hasNans(), ossimDrect::height(), ossimSensorModel::imageSize(), ossimDatumFactory::instance(), ossimImageGeometry::localToWorld(), ossimDrect::midPoint(), ossimAdjustableParameterInterface::resizeAdjustableParameterArray(), theDlatDparamGrid, theDlonDparamGrid, theHeightEnabledFlag, ossimDrect::ul(), ossimDatumFactory::wgs84(), ossimDrect::width(), ossimIpt::x, and ossimIpt::y.

 {
    theHeightEnabledFlag =  enableHeightFlag;
    
    if (!geom->getProjection() || imageBounds.hasNans())
       return;
 
    // don't let it get any smaller than 100, 100 pixels
    // on the input projector
    //
    // may want this to be adjusted by outside
    //
    const ossimDatum* targetDatum = ossimDatumFactory::instance()->wgs84();
    ossimIpt gridSize(2,2);
    ossimDpt gridOrigin(0,0);
    ossimGpt gpt;
    ossimGpt gpt2;
    ossimGpt bilinearGpt;
    resizeAdjustableParameterArray(0);
    double normSplit = 1.0;
    ossimIpt imageSize = ossimIpt(imageBounds.width(), imageBounds.height());
    double error = 0.0;
 
    ossimIpt imageOrigin = imageBounds.ul();
    ossimDpt spacing ((double)(imageBounds.width()-1)/(gridSize.x-1),
       (double)(imageBounds.height()-1)/(gridSize.y-1));
 
    if(theDlatDparamGrid)
    {
       delete [] theDlatDparamGrid;
       theDlatDparamGrid = NULL;
    }
    if(theDlonDparamGrid)
    {
       delete [] theDlonDparamGrid;
       theDlonDparamGrid = NULL;
    }
 
    geom->localToWorld(imageBounds.midPoint(), gpt);
 
    do
    {
       if(traceDebug())
       {
          ossimNotify(ossimNotifyLevel_DEBUG) << "Checking grid size " << gridSize << std::endl;
       }
 
       spacing = ossimDpt((double)(imageBounds.width()-1)/(gridSize.x-1),
          (double)(imageBounds.height()-1)/(gridSize.y-1));
 
       theLatGrid.setNullValue(ossim::nan());
       theLonGrid.setNullValue(ossim::nan());
       theDlatDhGrid.setNullValue(0.0);
       theDlonDhGrid.setNullValue(0.0);
       theLatGrid.setDomainType(ossimDblGrid::SAWTOOTH_90);
       theLonGrid.setDomainType(ossimDblGrid::WRAP_180);
       theLatGrid.initialize(gridSize, gridOrigin, spacing);
       theLonGrid.initialize(gridSize, gridOrigin, spacing);
       theDlatDhGrid.initialize(gridSize, gridOrigin, spacing);
       theDlonDhGrid.initialize(gridSize, gridOrigin, spacing);
       ossim_int32 x, y;
 
       for(y = 0; y < gridSize.y; ++y)
       {
          for(x = 0; x < gridSize.x; ++x)
          {
             ossimDpt norm((double)x/(double)(gridSize.x-1),
                (double)y/(double)(gridSize.y-1));
 
             ossimDpt pt(imageOrigin.x + norm.x*(imageSize.x-1),
                imageOrigin.y + norm.y*(imageSize.y-1));
 
             geom->localToWorld(pt, gpt);
             double h = gpt.height();
             if(ossim::isnan(h))
             {
                h += heightDelta;
             }
             ossimDpt fullPt;
             geom->rnToFull(pt, 0, fullPt);
             geom->getProjection()->lineSampleHeightToWorld(fullPt, h, gpt2);
             gpt.changeDatum(targetDatum);
             gpt2.changeDatum(targetDatum);
 
             theLatGrid.setNode(x, y, gpt.latd());
             theLonGrid.setNode(x, y, gpt.lond());
 
             theDlatDhGrid.setNode(x, y, (gpt2.latd() - gpt.latd())/heightDelta);
             theDlonDhGrid.setNode(x, y, (gpt2.lond() - gpt.lond())/heightDelta);
          }
       }
       ossim_int32 upperY = 2*gridSize.y;
       ossim_int32 upperX = 2*gridSize.x;
       error = 0.0;
 
       // Set all base-class data members needed for subsequent calls to projection code:
       initializeModelParams(imageBounds);
 
       for(y = 0; ((y < upperY)&&(error < theInterpolationError)); ++y)
       {
          for(x = 0; ((x < upperX)&&(error<theInterpolationError)); ++x)
          {
             ossimDpt norm((double)x/(double)(upperX-1),
                (double)y/(double)(upperY-1));
 
             ossimDpt imagePoint(imageOrigin.x + norm.x*(imageSize.x-1),
                imageOrigin.y + norm.y*(imageSize.y-1));
             ossimDpt testIpt;
 
             geom->localToWorld(imagePoint, gpt);
             worldToLineSample(gpt, testIpt);
             error = (testIpt-imagePoint).length();
          }
       }
 
       gridSize.x *= 2;
       gridSize.y *= 2;
       normSplit *= .5;
    } while((error > theInterpolationError) &&
            ((imageSize.x*normSplit) > theMinGridSpacing) &&
            ((imageSize.y*normSplit) > theMinGridSpacing));
 
    gridSize = theLatGrid.size();
 
    ossimAdjustableParameterInterface* adjustableParameters = 
       PTR_CAST(ossimAdjustableParameterInterface, geom->getProjection());
    removeAllAdjustments();
    if(adjustableParameters&&makeAdjustableFlag)
    {
       if(adjustableParameters->getNumberOfAdjustableParameters() > 0)
       {
          newAdjustment(adjustableParameters->getNumberOfAdjustableParameters());
 
          int numberOfParams = getNumberOfAdjustableParameters();
          if(numberOfParams)
          {
             //***
             // Allocate adjustable parameter partials grids then assign:
             //***
             theDlatDparamGrid = new ossimDblGrid [numberOfParams];
             theDlonDparamGrid = new ossimDblGrid [numberOfParams];
             for(int paramIdx = 0; paramIdx < numberOfParams; ++ paramIdx)
             {
                theDlonDparamGrid[paramIdx].setNullValue(0.0);
                theDlatDparamGrid[paramIdx].setNullValue(0.0);
                theDlatDparamGrid[paramIdx].initialize(gridSize, gridOrigin, spacing);
                theDlonDparamGrid[paramIdx].initialize(gridSize, gridOrigin, spacing);
                setAdjustableParameter(paramIdx, 0.0);
                setParameterSigma(paramIdx, adjustableParameters->getParameterSigma(paramIdx));
                setParameterUnit(paramIdx, adjustableParameters->getParameterUnit(paramIdx));
                setParameterCenter(paramIdx, 0.0);
                setParameterDescription(paramIdx,
                   adjustableParameters->getParameterDescription(paramIdx));
 
                double oldParameter = adjustableParameters->getAdjustableParameter(paramIdx);
                adjustableParameters->setAdjustableParameter(paramIdx, 1.0, true);
                double adjust = adjustableParameters->computeParameterOffset(paramIdx);
                double deltaLat = 0;
                double deltaLon = 0;
                if(adjust != 0.0)
                {
                   for(int y = 0; y < gridSize.y; ++y)
                   {
                      for(int x = 0; x < gridSize.x; ++x)
                      {
                         ossimDpt norm((double)x/(double)(gridSize.x-1),
                            (double)y/(double)(gridSize.y-1));
 
                         ossimDpt pt(imageOrigin.x + norm.x*(imageSize.x-1),
                            imageOrigin.y + norm.y*(imageSize.y-1));
                         geom->localToWorld(pt, gpt);
 
                         gpt.changeDatum(targetDatum);
                         gpt2.latd(theLatGrid(pt));
                         gpt2.lond(theLonGrid(pt));
                         deltaLat = gpt.latd()-gpt2.latd();
                         deltaLon = gpt.lond()-gpt2.lond();
 
                         theDlatDparamGrid[paramIdx].setNode(x, y, deltaLat/adjust);
                         theDlonDparamGrid[paramIdx].setNode(x, y, deltaLon/adjust);
                      }
                   }
 
                   // The partials grids for this parameter are initialized, now initialize the
                   // grid's extrapolator:
                   theDlatDparamGrid[paramIdx].enableExtrapolation();
                   theDlonDparamGrid[paramIdx].enableExtrapolation();
                }
                adjustableParameters->setAdjustableParameter(paramIdx, oldParameter, true);
             }
          }
       }
    }
    getAdjustment(theInitialAdjustment);
 }

◆ dup()

virtual ossimObject* ossimCoarseGridModel::dup ( ) const

inlinevirtual

Returns pointer to a new instance, copy of this.

Implements ossimProjection.

Definition at line 85 of file ossimCoarseGridModel.h.

85 { return new ossimCoarseGridModel(*this); }

ossimCoarseGridModel::ossimCoarseGridModel

ossimCoarseGridModel()

Definition: ossimCoarseGridModel.cpp:58

◆ extrapolate()

ossimGpt ossimCoarseGridModel::extrapolate	(	const ossimDpt &	imgPt,
		const double &	height = `ossim::nan()`
	)		const

protectedvirtual

Implements its own extrapolation since this can be handled by ossimDblGrid.

Reimplemented from ossimSensorModel.

Definition at line 1048 of file ossimCoarseGridModel.cpp.

References lineSampleHeightToWorld().

Referenced by lineSampleToWorld(), and ossimH5GridModel::worldToLineSample().

 {
    ossimGpt gpt;
    lineSampleHeightToWorld(local_ip, height, gpt);
    return gpt;
 }

◆ imagingRay()

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

virtual

METHOD: imagingRay(image_point, &ossimEcefRay) Given an image point, returns a ray originating at some arbitrarily high point (ideally at the sensor position) and pointing towards the target. Implemented here but should be overriden for more efficient solution.

Reimplemented from ossimSensorModel.

Definition at line 578 of file ossimCoarseGridModel.cpp.

References ossimSensorModel::imagingRay().

 {
    ossimSensorModel::imagingRay(image_point, image_ray);
 }

◆ initAdjustableParameters()

void ossimCoarseGridModel::initAdjustableParameters ( )

virtual

Reimplemented from ossimAdjustableParameterInterface.

Definition at line 563 of file ossimCoarseGridModel.cpp.

References ossimAdjustableParameterInterface::addAdjustment(), ossimAdjustableParameterInterface::getNumberOfAdjustableParameters(), ossimAdjustableParameterInterface::setAdjustment(), and theInitialAdjustment.

 {
    if(getNumberOfAdjustableParameters() < 1)
    {
       addAdjustment(theInitialAdjustment, true);
    }
    else
    {
       setAdjustment(theInitialAdjustment, true);
    }
 }

◆ initializeModelParams()

void ossimCoarseGridModel::initializeModelParams ( ossimIrect imageBounds )

protected

Initializes base class data members after grids have been assigned.

It is assumed that theImageSize and the origin image point were already set.

Definition at line 413 of file ossimCoarseGridModel.cpp.

Referenced by ossimHdf5GridModel::initialize(), and ossimEnviCgModel::loadEnviGeocFile().

 {
    // NOTE: it is assumed that the grid size and spacing is the same for ALL grids:
    ossimIpt gridSize (theLatGrid.size());
    ossimDpt spacing  (theLatGrid.spacing());
    ossimDpt v[4];
    v[0].lat = theLatGrid.getNode(0,0);
    v[0].lon = theLonGrid.getNode(0,0);
    v[1].lat = theLatGrid.getNode(gridSize.x-1, 0);
    v[1].lon = theLonGrid.getNode(gridSize.x-1, 0);
    v[2].lat = theLatGrid.getNode(gridSize.x-1, gridSize.y-2);
    v[2].lon = theLonGrid.getNode(gridSize.x-1, gridSize.y-2);
    v[3].lat = theLatGrid.getNode(0, gridSize.y-2);
    v[3].lon = theLonGrid.getNode(0, gridSize.y-2);
 
    // Guaranty longitude values are -180 to 180
    for (int i=0; i<4; i++)
    {
       if (v[i].lon > 180.0)
          v[i].lon -= 360.0;
    }
 
    theBoundGndPolygon = ossimPolygon(4, v);
    
    theImageSize  = ossimDpt(imageBounds.width(), imageBounds.height());
    theRefImgPt   = imageBounds.midPoint();
    theRefGndPt.lat = theLatGrid(theRefImgPt);
    theRefGndPt.lon = theLonGrid(theRefImgPt);
    
    ossimDpt ref_ip_dx (theRefImgPt.x+1.0, theRefImgPt.y    );
    ossimDpt ref_ip_dy (theRefImgPt.x    , theRefImgPt.y+1.0);
    ossimGpt ref_gp_dx (theLatGrid(ref_ip_dx), theLonGrid(ref_ip_dx));
    ossimGpt ref_gp_dy (theLatGrid(ref_ip_dy), theLonGrid(ref_ip_dy));
 
    theGSD.x   = theRefGndPt.distanceTo(ref_gp_dx);
    theGSD.y   = theRefGndPt.distanceTo(ref_gp_dy);
 
    theMeanGSD = (theGSD.line + theGSD.samp)/2.0;
    theImageClipRect  = imageBounds;
    theSubImageOffset = imageBounds.ul();
 }

◆ isAffectedByElevation()

bool ossimCoarseGridModel::isAffectedByElevation ( ) const

virtual

image to ground faster

Overrides base ossimSensorModel::isAffectedByElevation method.

Returns: true if height enabled, false if not.

Reimplemented from ossimSensorModel.

Definition at line 1055 of file ossimCoarseGridModel.cpp.

References theHeightEnabledFlag.

 {
    return theHeightEnabledFlag;
 }

◆ lineSampleHeightToWorld()

void ossimCoarseGridModel::lineSampleHeightToWorld	(	const ossimDpt &	image_point,
		const double &	heightEllipsoid,
		ossimGpt &	world_pt
	)		const

virtual

This is the virtual that performs the actual work of projecting the image point to the earth at some specified elevation.

Implements ossimSensorModel.

Definition at line 513 of file ossimCoarseGridModel.cpp.

Referenced by extrapolate(), lineSampleToWorld(), and ossimH5GridModel::worldToLineSample().

 {
    if (traceExec()) ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::lineSampleHeightToWorld: entering..." << std::endl;
    
    if(theLatGrid.size().x < 1 ||
       theLatGrid.size().y < 1)
    {
       worldPt.makeNan();
       return;
    }
 
    double height = (ossim::isnan(arg_hgt_above_ellipsoid)) ? 0.0 : arg_hgt_above_ellipsoid;
 
    // Note that there is no check for image point outside of valid image rect because this model
    // uses the extrapolation inherent to the ossimDblGrid.
 
    // The image point may correspond to an offset sub-image. Need to convert
    // to full image space before anything:
    ossimDpt ip = lineSampPt + theSubImageOffset;
    
    // Establish the interpolated values from the grids:
    worldPt.lat = theLatGrid(ip);
    worldPt.lon = theLonGrid(ip);
    worldPt.hgt = height;
 
    if(theHeightEnabledFlag)
    {
       // Adjust horizontally due to elevation:
       worldPt.lat += theDlatDhGrid(ip)*height;
       worldPt.lon += theDlonDhGrid(ip)*height;
    }
    int numberOfParams = getNumberOfAdjustableParameters();
  
    // Now add increments due to adjustable parameter deltas:
    for (int p=0; p<numberOfParams; p++)
    {
        worldPt.lat += (theDlatDparamGrid[p](ip) * computeParameterOffset(p));
        worldPt.lon += (theDlonDparamGrid[p](ip) * computeParameterOffset(p));
    }
 
    worldPt.limitLonTo180();
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::lineSampleHeightToWorld: returning..." << std::endl;
 }

◆ lineSampleToWorld()

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

virtual

Implementation of base-class pure virtual projection methods. These methods may be overriden by derived classes if those have more efficient schemes. The implementations here are iterative (relatively slow). Both depend on calls to the pure virtual lineSampleHeightToWorld() method.

Reimplemented from ossimSensorModel.

Definition at line 483 of file ossimCoarseGridModel.cpp.

References extrapolate(), ossimSensorModel::insideImage(), lineSampleHeightToWorld(), ossimSensorModel::lineSampleToWorld(), and theHeightEnabledFlag.

 {
    if(!theHeightEnabledFlag)
    {
       //
       // Extrapolate if image point is outside image:
       //
       if (!insideImage(image_point))
       {
          gpt = extrapolate(image_point);
          return;
       }
       
       lineSampleHeightToWorld(image_point, 0.0, gpt);
    }
    else
    {
       ossimSensorModel::lineSampleToWorld(image_point, gpt);
    }
 }

◆ loadCoarseGrid()

bool ossimCoarseGridModel::loadCoarseGrid ( const ossimFilename & cgFileName )

Loads the coarse grid from the specified file.

Returns: Returns true if successful.

Definition at line 855 of file ossimCoarseGridModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadCoarseGrid: entering..." << std::endl;
 
    ossimDpt v[4];
    ossimIpt grid_size;
 
    //***
    // Open existing grid file:
    //***
    ifstream instream (cgFileName.chars());
    if (!instream.is_open())
    {
       theErrorStatus++;
       if(traceDebug())
       {
          ossimNotify(ossimNotifyLevel_FATAL) << "FATAL ossimCoarseGridModel::loadCoarseGrid: Error encountered opening coarse grid file <" << cgFileName
          << ">. Check that the file exists and is readable." << std::endl;
       }
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "CEBUG ossimCoarseGridModel::loadCoarseGrid: returning with error..." << std::endl;
       return false;
    }
    theGridFilename = cgFileName;
    if(theDlatDparamGrid)
      {
        delete [] theDlatDparamGrid;
        theDlatDparamGrid = NULL;
      }
    if(theDlonDparamGrid)
      {
        delete [] theDlonDparamGrid;
        theDlonDparamGrid = NULL;
      }
    //
    // Reallocate memory:
    //
    int numberOfParams = getNumberOfAdjustableParameters();
    if(numberOfParams)
    {
       theDlatDparamGrid = new ossimDblGrid [numberOfParams];
       theDlonDparamGrid = new ossimDblGrid [numberOfParams];
    }
    //***
    // Let each grid object read itself from the input file:
    //***
    if (!theLatGrid.load(instream))
    {
       ++theErrorStatus;
       return false;
    }
    if (!theLonGrid.load(instream))
    {
       ++theErrorStatus;
       return false;
    }
    if (!theDlatDhGrid.load(instream))
    {
       ++theErrorStatus;
       return false;
    }
    if (!theDlonDhGrid.load(instream))
    {
       ++theErrorStatus;
       return false;
    }
 
    for (int p=0; p<numberOfParams; p++)
    {
       if (!theDlatDparamGrid[p].load(instream))
       {
          ++theErrorStatus;
          return false;
       }
       if (!theDlonDparamGrid[p].load(instream))
       {
          ++theErrorStatus;
          return false;
       }
    }
 
    //***
    // Initialize the bounding ground rectangle:
    //***
    grid_size = theLatGrid.size();
 
    v[0].lat = theLatGrid(0,0);
    v[0].lon = theLonGrid(0,0);
    v[1].lat = theLatGrid(theImageSize.x-1, 0);
    v[1].lon = theLonGrid(theImageSize.x-1, 0);
    v[2].lat = theLatGrid(theImageSize.x-1, theImageSize.y-1);
    v[2].lon = theLonGrid(theImageSize.x-1, theImageSize.y-1);
    v[3].lat = theLatGrid(0, theImageSize.y-1);
    v[3].lon = theLonGrid(0, theImageSize.y-1);
 
    theBoundGndPolygon = ossimPolygon(4, v);
    
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadCoarseGrid: returning..." << std::endl;
    return true;
 }

◆ loadState()

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

virtual

Fulfills ossimObject base-class pure virtuals.

Loads geometry KWL files.

Returns: Returns true if successful.

Reimplemented from ossimSensorModel.

Reimplemented in ossimH5GridModel.

Definition at line 680 of file ossimCoarseGridModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadState: entering..." << std::endl;
 
    clearErrorStatus();
    const char* value;
    bool success;
    
    //***
    // Assure this keywordlist contains correct type info:
    //***
    value = kwl.find(prefix, ossimKeywordNames::TYPE_KW);
    if (!value || (strcmp(value, "ossimCoarseGridModel")))
    {
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadState:  returning..." << std::endl;
       theErrorStatus++;
       return false;
    }
    value = kwl.find(prefix, "height_enabled_flag");
    if(value)
    {
       theHeightEnabledFlag = ossimString(value).toBool();
    }
 
    //***
    // Pass on to the base-class for parsing first:
    //***
    success = ossimSensorModel::loadState(kwl, prefix);
    if (!success)
    {
       theErrorStatus++;
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadState: returning with error..." << std::endl;
       return false;
    }
 
    //***
    // Look for geom filename or explicit grid filename to establish path to grid:
    //***
    theGridFilename = kwl.find(prefix, GRID_FILE_NAME_KW);
 
    if (!theGridFilename.isReadable())
    {
       //---
       // This is set in ossimProjectionFactoryBase::createProjectionFromGeometryFile
       // so we can derive the dot.ocg from it.
       //---
       ossimFilename alt_path_to_grid = kwl.find("kwl_source");
       theGridFilename = alt_path_to_grid.setExtension(DEFAULT_GRID_FILE_EXT);
       
       if (!theGridFilename.isReadable())
       {
          ossimFilename alt_path_to_grid = kwl.find(prefix, ossimKeywordNames::GEOM_FILE_KW);
          theGridFilename = alt_path_to_grid.setExtension(DEFAULT_GRID_FILE_EXT);
       }
    }
    
    if (!theGridFilename.isReadable())
    {
       ossimNotify(ossimNotifyLevel_FATAL) << "ossimCoarseGridModel::loadState() -- Error "
          "encountered opening coarse grid file at "<< "<" <<theGridFilename << ">." << std::endl;
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadState: returning with error..." << std::endl;
       theErrorStatus++;
       return false;
    }
 
    // Load the coarse grid file specified in KWL. This method resets the
    // theErrorStatus to OSSIM_OK if successful:
    if (!loadCoarseGrid(theGridFilename))
    {
       theErrorStatus++;
       return false;
    }
 
    // crossesDateline legacy. No longer saved.
    bool crossesDateline = false;
    kwl.getBoolKeywordValue(crossesDateline, CROSSES_DATELINE_KW, prefix);
    if (crossesDateline)
       theLonGrid.setDomainType(ossimDblGrid::WRAP_360);
 
    // Add the coarse grid filename to list of support files being referenced for logging purposes:
    ossimSupportFilesList::instance()->add(theGridFilename.expand());
 
    ossimString initAdjPrefix = ossimString(prefix) + "init_adjustment.";
    theInitialAdjustment.loadState(kwl, initAdjPrefix.c_str());
 
    if((ossim::isnan(theRefGndPt.hgt)) ||
       (theRefGndPt.hgt == 0))
    {
       theRefGndPt.hgt = ossimElevManager::instance()->getHeightAboveEllipsoid(theRefGndPt);
       if(theRefGndPt.hgt < 0)
       {
          theRefGndPt.hgt = fabs(theRefGndPt.hgt);
       }
    }
       
    if(theInitialAdjustment.getNumberOfAdjustableParameters() < 1)
    {
       getAdjustment(theInitialAdjustment);
    }
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::loadState: returning..." << std::endl;
    if (theErrorStatus)
       return false;
    
    return true;
 }

◆ print()

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

virtual

Extends base-class implementation.

Dumps contents of object to ostream.

Reimplemented from ossimSensorModel.

Definition at line 591 of file ossimCoarseGridModel.cpp.

References ossimDblGrid::getNode(), ossimDblGrid::size(), size, ossimDblGrid::spacing(), ossimSensorModel::theBoundGndPolygon, theDlatDhGrid, theDlonDhGrid, theGridFilename, ossimSensorModel::theGSD, ossimSensorModel::theImageID, ossimSensorModel::theImageSize, theLatGrid, theLonGrid, ossimSensorModel::theRefGndPt, ossimSensorModel::theRefImgPt, ossimSensorModel::theSensorID, ossimIpt::x, and ossimIpt::y.

 {
    out << "\nDump of ossimCoarseGridModel object at: " << this << "\n"
        << "\n           Grid File Name: " << theGridFilename 
        << "\n                 Image ID: " << theImageID 
        << "\n                   Sensor: " << theSensorID
        << "\n  Image Size (rows, cols): " << theImageSize
        << "\n      Ref Pt (samp, line): " << theRefImgPt 
        << "\n   Ref Pt (lat, lon, hgt): " << theRefGndPt 
        << "\n           GSD (row, col): " << theGSD 
        << "\n  Bounding Ground Polygon: " << theBoundGndPolygon << endl;
 //      << "\n         Number of Params: " << theNumParams << "\n"<<endl;
    
    char buf[256];
    ossimIpt size (theLatGrid.size());
    ossimDpt spacing (theLatGrid.spacing());
    int line, samp;
    ossimIpt node;
    
    out << "[ line,  samp]        lat        lon         dLat/dH      dLon/dH\n"
        << "-------------------------------------------------------------------"
        <<endl;
 
    for (node.y=0; node.y<size.y; node.y++)
    {
       line = (int) (node.y*spacing.y);
       
       for (node.x=0; node.x<size.x; node.x++)
       {
          samp = (int) (node.x*spacing.x);
          
          sprintf(buf, "[%5d, %5d]    %+9.5f  %+10.5f    %+11.4e  %+11.4e",
                   line, samp,
                   theLatGrid.getNode(node),
                   theLonGrid.getNode(node),
                   theDlatDhGrid.getNode(node),
                   theDlonDhGrid.getNode(node));
          out << buf << endl;
       }
       out <<"-----------------------------------------------------------------"
           <<endl;
    }
 
    out << "\n\nDump of lat/lon Partials w.r.t. Adjustable Parameters:"<<endl;
    out << "\nEnd Dump of ossimCoarseGridModel.\n" <<  endl;
    return out;
 }

◆ reallocateGrid()

void ossimCoarseGridModel::reallocateGrid ( const ossimIpt & size )

protected

Deletes existing allocated memory and reallocates new space.

This may happen if a new grid is loaded over an existing one.

Definition at line 962 of file ossimCoarseGridModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::reallocateGrid:entering..." << endl;
 
 
    //***
    // Deallocate memory:
    //***
    if(theDlatDparamGrid)
      {
        delete [] theDlatDparamGrid;
        theDlatDparamGrid = NULL;
      }
    if(theDlonDparamGrid)
      {
        delete [] theDlonDparamGrid;
        theDlonDparamGrid = NULL;
      }
    //***
    // determine grid spacing given new info:
    //***
    ossimDpt spacing ((double)(theImageSize.x-1)/(double)(grid_size.x-1),
                      (double)(theImageSize.y-1)/(double)(grid_size.y-1));
 
    //***
    // Allocate all:
    //***
    ossimDpt grid_origin(0.0, 0.0);
    theLatGrid.setNullValue(ossim::nan());
    theLonGrid.setNullValue(ossim::nan());
    theDlatDhGrid.setNullValue(0.0);
    theDlonDhGrid.setNullValue(0.0);
    theLatGrid.initialize(grid_size, grid_origin, spacing);
    theLonGrid.initialize(grid_size, grid_origin, spacing);
    theDlatDhGrid.initialize(grid_size, grid_origin, spacing);
    theDlonDhGrid.initialize(grid_size, grid_origin, spacing);
    
    int numberOfParams = getNumberOfAdjustableParameters();
    if(numberOfParams)
    {
       
       theDlatDparamGrid = new ossimDblGrid [numberOfParams];
       theDlonDparamGrid = new ossimDblGrid [numberOfParams];
    }
    for (int p=0; p<numberOfParams; p++)
    {
       theDlonDparamGrid[p].setNullValue(0.0);
       theDlatDparamGrid[p].setNullValue(0.0);
       theDlatDparamGrid[p].initialize(grid_size, grid_origin, spacing);
       theDlonDparamGrid[p].initialize(grid_size, grid_origin, spacing);
    }
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::reallocateGrid: returning..." << std::endl;
    return;
 }

◆ saveCoarseGrid()

bool ossimCoarseGridModel::saveCoarseGrid ( const ossimFilename & cgFileName ) const

Saves the coarse grid to the specified file.

Returns: Returns true if successful.

Definition at line 795 of file ossimCoarseGridModel.cpp.

Referenced by ossimHdf5Tool::execute().

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveCoarseGrid: entering..." << std::endl;
 
    // Create and open grid file as stream:
    theGridFilename = fileName.expand();
    theGridFilename.setExtension(DEFAULT_GRID_FILE_EXT);
    ofstream outstream (theGridFilename.chars());
    if (!outstream.is_open())
    {
       theErrorStatus++;
       ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveCoarseGrid: Error "
          "encountered creating coarse grid file <" << theGridFilename<< ">. Check that directory "
          "exists and is writable." << std::endl;
       if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveCoarseGrid: returning with error..." << std::endl;
       return false;
    }
 
    // Let each grid object write itself to the output file:
    theLatGrid.save(outstream, "Latitude Grid");
    theLonGrid.save(outstream, "Longitude Grid");
    theDlatDhGrid.save(outstream, "dLat/dH Grid");
    theDlonDhGrid.save(outstream, "dLon_dH Grid");
 
    ossimString descr;
    int numberOfParams = getNumberOfAdjustableParameters();
    for (int p=0; p<numberOfParams; p++)
    {
       descr = getParameterDescription(p) + " dLat_dParam Grid";
       theDlatDparamGrid[p].save(outstream, descr.chars());
       descr = getParameterDescription(p) + " dLon_dParam Grid";
       theDlonDparamGrid[p].save(outstream, descr.chars());
    }
    
    // Since the geom file is needed in the same path as the grid file, take this opportunity to 
    // write the geom file out as well:
    ossimFilename geom_file (theGridFilename);
    geom_file.setExtension(DEFAULT_GEOM_FILE_EXT);
    ossimKeywordlist kwl;
    saveState(kwl);
    kwl.write(geom_file);
 
    // Add to the list of support files referenced (though technically it has not yet been 
    // referenced, but will be next time this image is opened):
    ossimSupportFilesList::instance()->add(geom_file);
    ossimSupportFilesList::instance()->add(theGridFilename);
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveCoarseGrid: returning..." << std::endl;
    return true;
 }

◆ saveState()

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

virtual

Fulfills ossimObject base-class pure virtuals.

Saves geometry KWL files.

Returns: Returns true if successful.

Reimplemented from ossimSensorModel.

Reimplemented in ossimH5GridModel, ossimHdf5GridModel, and ossimEnviCgModel.

Definition at line 647 of file ossimCoarseGridModel.cpp.

Referenced by ossimEnviCgModel::saveState(), ossimHdf5GridModel::saveState(), and ossimH5GridModel::saveState().

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveState: entering..." << std::endl;
 
    kwl.add(prefix, ossimKeywordNames::TYPE_KW, "ossimCoarseGridModel");
 
    //---
    // Save just the file part so that the geometry stuff is relocatable.
    // Code was added to ossimProjectionFactoryBase::createProjectionFromGeometryFile
    // to handle this.
    //---
    kwl.add( prefix, GRID_FILE_NAME_KW, theGridFilename.file() );
    
    kwl.add(prefix, "height_enabled_flag", theHeightEnabledFlag, true);
    ossimSensorModel::saveState(kwl, prefix);
    ossimString initAdjPrefix = ossimString(prefix) + "init_adjustment.";
    theInitialAdjustment.saveState(kwl, initAdjPrefix);
    
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::saveState: returning..." << std::endl;
 
    return true;
 }

◆ setInterpolationError()

void ossimCoarseGridModel::setInterpolationError ( double error = .1 )

static

This is used when building a grid from a projector.

You can set the interpolation error. The default is subpixel accuracy (within .1 of a pixel).

Definition at line 399 of file ossimCoarseGridModel.cpp.

References theInterpolationError.

Referenced by ossimplugins::ossimGeometricSarSensorModel::createReplacementOCG().

 {
    theInterpolationError = error;
 }

◆ setMinGridSpacing()

void ossimCoarseGridModel::setMinGridSpacing ( ossim_int32 minSpacing = 100 )

static

Definition at line 404 of file ossimCoarseGridModel.cpp.

References theMinGridSpacing.

 {
    theMinGridSpacing = minSpacing;
 }

◆ useForward()

virtual bool ossimCoarseGridModel::useForward ( ) const

inlinevirtual

ossimOptimizableProjection

Implements ossimOptimizableProjection.

Definition at line 111 of file ossimCoarseGridModel.h.

111 {return false;}

◆ writeGeomTemplate()

void ossimCoarseGridModel::writeGeomTemplate ( ostream & os )

static

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

Definition at line 1024 of file ossimCoarseGridModel.cpp.

 {
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::writeGeomTemplate: entering..." << std::endl;
 
    os <<
       "//**************************************************************\n"
       "// Template for OCG model kewordlist\n"
       "//**************************************************************\n"
       << ossimKeywordNames::TYPE_KW << ": " << MODEL_TYPE << endl;
 
    ossimSensorModel::writeGeomTemplate(os);
    
    os << "//\n"
       << "// Derived-class ossimCoarseGridModel Keywords:\n"
       << "//\n"
       << GRID_FILE_NAME_KW << ": <string>\n" << endl;
 
    if (traceExec())  ossimNotify(ossimNotifyLevel_DEBUG) << "DEBUG ossimCoarseGridModel::writeGeomTemplate: returning..." << std::endl;
    return;
 }

Member Data Documentation

◆ theDlatDhGrid

ossimDblGrid ossimCoarseGridModel::theDlatDhGrid

protected

Definition at line 132 of file ossimCoarseGridModel.h.

Referenced by ossimHdf5GridModel::initialize(), ossimEnviCgModel::loadEnviGeocFile(), and print().

◆ theDlatDparamGrid

ossimDblGrid* ossimCoarseGridModel::theDlatDparamGrid

protected

Definition at line 134 of file ossimCoarseGridModel.h.

Referenced by buildGrid().

◆ theDlonDhGrid

ossimDblGrid ossimCoarseGridModel::theDlonDhGrid

protected

Definition at line 133 of file ossimCoarseGridModel.h.

Referenced by ossimHdf5GridModel::initialize(), ossimEnviCgModel::loadEnviGeocFile(), and print().

◆ theDlonDparamGrid

ossimDblGrid* ossimCoarseGridModel::theDlonDparamGrid

protected

Definition at line 135 of file ossimCoarseGridModel.h.

Referenced by buildGrid().

◆ theGridFilename

ossimFilename ossimCoarseGridModel::theGridFilename

mutableprotected

Definition at line 129 of file ossimCoarseGridModel.h.

Referenced by print().

◆ theHeightEnabledFlag

bool ossimCoarseGridModel::theHeightEnabledFlag

protected

Definition at line 140 of file ossimCoarseGridModel.h.

Referenced by buildGrid(), ossimHdf5GridModel::initialize(), ossimH5GridModel::initializeModelParams(), isAffectedByElevation(), lineSampleToWorld(), and ossimEnviCgModel::loadEnviGeocFile().

◆ theInitialAdjustment

ossimAdjustmentInfo ossimCoarseGridModel::theInitialAdjustment

protected

Definition at line 139 of file ossimCoarseGridModel.h.

Referenced by initAdjustableParameters().

◆ theInterpolationError

double ossimCoarseGridModel::theInterpolationError = .1

staticprotected

Definition at line 137 of file ossimCoarseGridModel.h.

Referenced by setInterpolationError().

◆ theLatGrid

ossimDblGrid ossimCoarseGridModel::theLatGrid

protected

Definition at line 130 of file ossimCoarseGridModel.h.

Referenced by ossimH5GridModel::debugDump(), ossimH5GridModel::getWktFootprint(), ossimHdf5GridModel::initialize(), ossimH5GridModel::initializeModelParams(), initializeModelParams(), ossimEnviCgModel::loadEnviGeocFile(), ossimH5GridModel::ossimH5GridModel(), ossimHdf5GridModel::ossimHdf5GridModel(), print(), and ossimH5GridModel::setGridNodes().

◆ theLonGrid

ossimDblGrid ossimCoarseGridModel::theLonGrid

protected

Definition at line 131 of file ossimCoarseGridModel.h.

Referenced by ossimHdf5GridModel::crossesDateline(), ossimH5GridModel::debugDump(), ossimH5GridModel::getWktFootprint(), ossimHdf5GridModel::initialize(), ossimH5GridModel::initializeModelParams(), initializeModelParams(), ossimEnviCgModel::loadEnviGeocFile(), print(), and ossimH5GridModel::setGridNodes().

◆ theMinGridSpacing

ossim_int32 ossimCoarseGridModel::theMinGridSpacing = 100

staticprotected

Definition at line 138 of file ossimCoarseGridModel.h.

Referenced by setMinGridSpacing().

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimCoarseGridModel() [1/4]

◆ ossimCoarseGridModel() [2/4]

◆ ossimCoarseGridModel() [3/4]

◆ ossimCoarseGridModel() [4/4]

◆ ~ossimCoarseGridModel()

Member Function Documentation

◆ buildGrid() [1/2]

◆ buildGrid() [2/2]

◆ dup()

◆ extrapolate()

◆ imagingRay()

◆ initAdjustableParameters()

◆ initializeModelParams()

◆ isAffectedByElevation()

◆ lineSampleHeightToWorld()

◆ lineSampleToWorld()

◆ loadCoarseGrid()

◆ loadState()

◆ print()

◆ reallocateGrid()

◆ saveCoarseGrid()

◆ saveState()

◆ setInterpolationError()

◆ setMinGridSpacing()

◆ useForward()

◆ writeGeomTemplate()

Member Data Documentation

◆ theDlatDhGrid

◆ theDlatDparamGrid

◆ theDlonDhGrid

◆ theDlonDparamGrid

◆ theGridFilename

◆ theHeightEnabledFlag

◆ theInitialAdjustment

◆ theInterpolationError

◆ theLatGrid

◆ theLonGrid

◆ theMinGridSpacing