#include <ossimHdf5GridModel.h>

Inheritance diagram for ossimHdf5GridModel:

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

virtual	~ossimHdf5GridModel ()
	virtual destructor More...

bool	initialize (ossimHdf5 *hdf5, const ossimString &projDataPath="")
	Initializes from an open HDF5 file. More...

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix) const
	Makes sure that the "type" keyword entry reflects the base class, not this one. More...

Public Member Functions inherited from ossimCoarseGridModel
	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	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...

Protected Member Functions
bool	initCoarseGrid (const char *datasetName, ossimDblGrid &coarseGrid)

bool	crossesDateline ()

Protected Member Functions inherited from ossimCoarseGridModel
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
ossimPolygon	m_boundGndPolygon

ossimRefPtr< ossimHdf5 >	m_hdf5

ossimIpt	m_imageSize

ossimString	m_projDataPath

Protected Attributes inherited from ossimCoarseGridModel
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

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

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

Static Public Member Functions inherited from ossimCoarseGridModel
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)

Static Protected Attributes inherited from ossimCoarseGridModel
static double	theInterpolationError = .1

static ossim_int32	theMinGridSpacing = 100

Detailed Description

Definition at line 35 of file ossimHdf5GridModel.h.

Constructor & Destructor Documentation

◆ ossimHdf5GridModel()

ossimHdf5GridModel::ossimHdf5GridModel ( )

default constructor.

Definition at line 44 of file ossimHdf5GridModel.cpp.

References ossimDblGrid::SAWTOOTH_90, ossimDblGrid::setDomainType(), and ossimCoarseGridModel::theLatGrid.

 : ossimCoarseGridModel()
 {
    theLatGrid.setDomainType(ossimDblGrid::SAWTOOTH_90);
 }

◆ ~ossimHdf5GridModel()

ossimHdf5GridModel::~ossimHdf5GridModel ( )

virtual

virtual destructor

Definition at line 50 of file ossimHdf5GridModel.cpp.

51 {

52 }

Member Function Documentation

◆ crossesDateline()

bool ossimHdf5GridModel::crossesDateline ( )

protected

Definition at line 362 of file ossimHdf5GridModel.cpp.

References ossimDblGrid::getNode(), ossimDblGrid::setDomainType(), ossimDblGrid::size(), size, ossimCoarseGridModel::theLonGrid, ossimDblGrid::WRAP_180, ossimDblGrid::WRAP_360, x, and y.

Referenced by initialize().

 {
    bool crossesDateline = false;
 
    ossim_int32 longitude = 0;
    bool found179 = false;
    bool found181 = false;
 
    ossimIpt size (theLonGrid.size());
    //double left, right;
    //int xr = size.x-1;
 
    for (ossim_uint32 y=0; (y<(ossim_uint32)size.y) && !crossesDateline; ++y )
    {
 #if 0
       left  = theLonGrid.getNode(0,  y);
       right = theLonGrid.getNode(xr, y);
       if (left > right)
          crossesDateline = true;
 #endif
 
       for ( ossim_uint32 x = 0; x<(ossim_uint32)size.x; ++x)
       {
          longitude = (ossim_int32) theLonGrid.getNode(x,  y); // Cast to integer.
 
          // look for 179 -> -179...
          if ( !found179 )
          {
             if ( longitude == 179 )
             {
                found179 = true;
                continue;
             }
          }
          else // found179 == true
          {
             if ( longitude == 178 )
             {
                break; // Going West, 179 -> 178
             }
             else if ( longitude == -179 )
             {
                crossesDateline = true;
                break;
             }
          }
 
          // look for -179 -> 179...
          if ( !found181 )
          {
             if ( longitude == -179 )
             {
                found181 = true;
                continue;
             }
          }
          else // found181 == true
          {
             if ( longitude == -178 )
             {
                break; // Going East -179 -> -178
             }
             else if ( longitude == 179 )
             {
                crossesDateline = true;
                break;
             }
          }
       }
    }
 
    if ( crossesDateline )
       theLonGrid.setDomainType(ossimDblGrid::WRAP_360);
    else
       theLonGrid.setDomainType(ossimDblGrid::WRAP_180);
 
    return crossesDateline;
 }

◆ initCoarseGrid()

bool ossimHdf5GridModel::initCoarseGrid	(	const char *	datasetName,
		ossimDblGrid &	coarseGrid
	)

protected

Definition at line 98 of file ossimHdf5GridModel.cpp.

References ossimString::chars(), ossimHdf5::findDatasetByName(), ossimHdf5::findGroupByName(), m_hdf5, m_imageSize, m_projDataPath, ossimIpt::x, and ossimIpt::y.

Referenced by initialize().

 {
    ostringstream xmsg;
 
    // Convention used: (u,v) is file space, (x,y) is CG space
    Group* group = m_hdf5->findGroupByName(m_projDataPath.chars(), 0, true);
    DataSet* dataset  = m_hdf5->findDatasetByName(datasetName, group, true);
    if (dataset == NULL)
    {
       xmsg  << "ossimHdf5GridModel:"<<__LINE__
             <<" ERROR: Could not find dataset \""<<datasetName<<"\" in file.";
       throw ossimException(xmsg.str());
    }
 
    // Get dataspace of the dataset.
    DataSpace dataSpace = dataset->getSpace();
    const ossim_int32 DIM_COUNT = dataSpace.getSimpleExtentNdims();
    if ( DIM_COUNT != 2  )
       return false;
 
    // Get the extents. dimsOut[0] is height, dimsOut[1] is width:
    std::vector<hsize_t> dimsOut(DIM_COUNT);
    dataSpace.getSimpleExtentDims( &dimsOut.front(), 0 );
    m_imageSize.y = dimsOut[0];
    m_imageSize.x = dimsOut[1];
 
    // Initialize the ossimDblGrid. Round up if size doesn't fall on end pixel.
    ossimDpt dspacing (GRID_SAMPLING_INTERVAL, GRID_SAMPLING_INTERVAL);
    ossim_uint32 gridRows = m_imageSize.y / GRID_SAMPLING_INTERVAL + 1;
    ossim_uint32 gridCols = m_imageSize.x / GRID_SAMPLING_INTERVAL + 1;
    if ( m_imageSize.y % GRID_SAMPLING_INTERVAL)
       ++gridRows;
    if ( m_imageSize.x % GRID_SAMPLING_INTERVAL)
       ++gridCols;
    ossimIpt gridSize (gridCols, gridRows);
 
    // The grid as used in base class, has UV-space always at 0,0 origin
    ossimDpt gridOrigin(0.0,0.0);
    coarseGrid.setNullValue(ossim::nan());
    coarseGrid.initialize(gridSize, gridOrigin, dspacing);
 
    std::vector<hsize_t> inputCount(DIM_COUNT);
    std::vector<hsize_t> inputOffset(DIM_COUNT);
 
    inputOffset[0] = 0; // y_img is set below.
    inputOffset[1] = 0;
    inputCount[0] = 1; // y_img
    inputCount[1] = (hsize_t)m_imageSize.x; // x_img
 
    // Output dataspace dimensions. Reading a line at a time.
    const ossim_int32 OUT_DIM_COUNT = 3;
    std::vector<hsize_t> outputCount(OUT_DIM_COUNT);
    outputCount[0] = 1;    // band
    outputCount[1] = 1;    // y_img
    outputCount[2] = m_imageSize.x; // x_img
 
    // Output dataspace offset.
    std::vector<hsize_t> outputOffset(OUT_DIM_COUNT);
    outputOffset[0] = 0;
    outputOffset[1] = 0;
    outputOffset[2] = 0;
 
    ossimScalarType scalar = m_hdf5->getScalarType( *dataset);
    if ( scalar != OSSIM_FLOAT32 )
       return false;
 
    // See if we need to swap bytes:
    ossimEndian endian;
    bool needSwap = false;
    if (m_hdf5->getByteOrder(dataset) != ossim::byteOrder())
       needSwap = true;
    DataType dataType = dataset->getDataType();
 
    // Output dataspace always the same, width of one line.
    DataSpace bufferDataSpace( OUT_DIM_COUNT, &outputCount.front());
    bufferDataSpace.selectHyperslab( H5S_SELECT_SET,
                                     &outputCount.front(),
                                     &outputOffset.front() );
 
    //  Arrays to hold a single line of latitude longitude values.
    vector<ossim_float32> values(m_imageSize.x);
    ossim_float32 val = 0;
    hsize_t y_img = 0;
 
    // Line loop:
    for ( ossim_uint32 y = 0; y < gridRows; ++y )
    {
       // y_img = line in image space
       y_img = y*GRID_SAMPLING_INTERVAL;
       if ( y_img < (ossim_uint32) m_imageSize.y )
       {
          inputOffset[0] = y_img;
          dataSpace.selectHyperslab( H5S_SELECT_SET, &inputCount.front(), &inputOffset.front() );
 
          // Read data from file into the buffer.
          dataset->read( &(values.front()), dataType, bufferDataSpace, dataSpace );
          if ( needSwap )
             endian.swap( &(values.front()), m_imageSize.x );
 
          // Sample loop:
          hsize_t x_img = 0;
          for ( ossim_uint32 x = 0; x < gridCols; ++x )
          {
             // x_img = sample in image space
             x_img = x*GRID_SAMPLING_INTERVAL;
             if ( x_img < (ossim_uint32) m_imageSize.x )
             {
                val = values[x_img];
                if (ossim::isnan(val)) // Nulls in grid!
                {
                   xmsg  << "ossimHdf5GridModel:"<<__LINE__<<" encountered nans!";
                   throw ossimException(xmsg.str());
                }
             }
             else // Last column is outside of image bounds.
             {
                // Delta between last two latitude grid values.
                ossim_float32 val1 = coarseGrid.getNode( x-2, y );
                ossim_float32 val2 = coarseGrid.getNode( x-1, y );
                ossim_float32 spacing = val2 - val1;
                val = val2 + spacing;
 
 #if 0 /* Please leave for debug. (drb) */
                cout << "val1: " << val1 << " val2 " << val2<<endl;;
 #endif
             }
 
             coarseGrid.setNode( x, y, val );
 
 #if 0 /* Please leave for debug. (drb) */
             cout << "x,y,x_img,y_img,val:" << x << "," << y << ","<< x_img << "," << y_img << ","
                   << coarseGrid.getNode(x, y) << endl;
 #endif
          } // End sample loop.
       }
       else // Row is outside of image bounds:
       {
          // Sample loop:
          for ( ossim_uint32 x = 0; x < gridCols; ++x )
          {
             ossim_float32 val = ossim::nan();
             ossim_float32 val1 = coarseGrid.getNode( x, y-2 );
             ossim_float32 val2 = coarseGrid.getNode( x, y-1 );
             ossim_float32 spacing = val2 - val1;
             val = val2 + spacing;
            coarseGrid.setNode( x, y, val );
 
 #if 0 /* Please leave for debug. (drb) */
             hsize_t x_img = x*GRID_SPACING; // Sample in image space
             cout << "val1: " << val1 << " val2 " << val2
                   << "\nx,y,x_img,y_img,val:" << x << "," << y << ","
                   << x_img << "," << y_img << "," << val << endl;
 #endif
          } // End sample loop.
       } // Matches if ( y_img < m_imageSize.y ){...}else{
    } // End line loop.
 
    dataSpace.close();
 
    return true;
 
 #if 0
    // Original refactor code ###############################################
 
    // Verify dimensions:
    DataSpace dataSpace = dataset->getSpace();
    if (dataSpace.getSimpleExtentNdims() != 2)
    {
       xmsg << "ossimHdf5GridModel:"<<__LINE__<<" ERROR: lat/lon grid dataspace rank != 2.";
       throw ossimException(xmsg.str());
    }
 
    // Fetch size of grid in file. Extents are assumed to be the same for both lat and lon grids:
    hsize_t datExtents[2];
    dataSpace.getSimpleExtentDims(datExtents);
    if ((datExtents[0] < 2) || (datExtents[1] < 2))
    {
       xmsg << "ossimHdf5GridModel:"<<__LINE__<<" ERROR: lat/lon grid size is < 2.";
       throw ossimException(xmsg.str());
    }
 
    // Initialize the base class coarse grids:
    theImageSize  = ossimDpt(datExtents[0], datExtents[1]);
    ossimDrect uvRect(0, 0, theImageSize.u-1, theImageSize.v-1);
    ossimDpt cgSpacing(GRID_SAMPLING_INTERVAL, GRID_SAMPLING_INTERVAL);
    coarseGrid.initialize(uvRect, cgSpacing, ossim::nan());
 
    // Declare data of interest in file (the whole thing):
    hsize_t offset[2] = { 0, 0 };
    dataSpace.selectHyperslab( H5S_SELECT_SET, datExtents, offset );
 
    // Initialize dataspace for memory buffer needed by dataset read operation:
    hsize_t bufExtents[2] = { (hsize_t) theImageSize.u, 1 };
    DataSpace bufSpace(2, bufExtents);
    DataType dataType = dataset->getDataType();
    string cname = dataType.fromClass();
    cout << cname<<endl;
    if (dataType.getClass() != H5T_FLOAT)
    {
       xmsg << "ossimHdf5GridModel:"<<__LINE__<<" ERROR: lat/lon grid datatype must be float.";
       throw ossimException(xmsg.str());
    }
    ossim_float32* buffer = new ossim_float32 [bufExtents[0]]; // assumes float datatype
 
    // See if we need to swap bytes:
    ossimEndian* endian = 0;
    AtomType* atomType = dynamic_cast<AtomType*>(&dataType);
    if(atomType)
    {
       ossimByteOrder ossimByteOrder = ossim::byteOrder();
       H5T_order_t h5order = atomType->getOrder();
       if( ((h5order == H5T_ORDER_LE) && (ossimByteOrder != OSSIM_LITTLE_ENDIAN)) ||
             ((h5order == H5T_ORDER_BE) && (ossimByteOrder != OSSIM_BIG_ENDIAN)))
          endian = new ossimEndian();
    }
 
    // Loop over input grid rows, sampling according to desired interval to fill output
    // coarse grid, for latitude:
    ossimIpt cgGridSize (coarseGrid.size());
    ossim_uint32 x=0, y=0, u=0, v=0;
    for ( ; (y<cgGridSize.y) && (v<theImageSize.v); ++y, v+=GRID_SAMPLING_INTERVAL )
    {
       offset[1] = v; // offset[0] always = 0
       dataSpace.selectHyperslab( H5S_SELECT_SET, bufExtents, offset);
       dataset->read( buffer, dataType,  bufSpace, dataSpace );
       if ( endian )
          endian->swap( buffer, bufExtents[0] );
 
       // Need to subsample the input row and save into coarse grid:
       for ( x=0, u=0; (x<cgGridSize.x)&&(u<theImageSize.u); ++x, u+=GRID_SAMPLING_INTERVAL)
       {
          { // TODO REMOVE DEBUG BLOCK
             cout<<datasetName<<" ("<<x<<", "<<y<<"): "<<buffer[u]<<endl;
          }
 
          if ( ossim::isnan(buffer[u]))
          {
             xmsg << "ossimHdf5GridModel:"<<__LINE__<<" ERROR: encountered nans in lat/lon grid.";
             throw ossimException(xmsg.str());
          }
          coarseGrid.setNode( x, y, buffer[u] );
       }
 
       // Check if last column is outside of image bounds.
       if (x < cgGridSize.x)
          coarseGrid.setNode( x, y, coarseGrid.getNode( x-1, y ) );
    }
 
    // Check if last row is outside of image bounds.
    if (y < cgGridSize.y)
    {
       for ( x=0; x<cgGridSize.x; ++x)
          coarseGrid.setNode( x, y, coarseGrid.getNode( x, y-1 ) );
    }
 
    delete dataset;
    delete buffer;
    delete endian;
 
    return true;
 #endif
 }

◆ initialize()

bool ossimHdf5GridModel::initialize	(	ossimHdf5 *	hdf5,
		const ossimString &	projDataPath = `""`
	)

Initializes from an open HDF5 file.

Definition at line 54 of file ossimHdf5GridModel.cpp.

References crossesDateline(), initCoarseGrid(), ossimDblGrid::initialize(), ossimCoarseGridModel::initializeModelParams(), m_hdf5, m_imageSize, m_projDataPath, ossimDblGrid::maxValue(), ossimDblGrid::minValue(), ossimDblGrid::origin(), ossimNotify(), ossimNotifyLevel_FATAL, ossimDblGrid::size(), ossimDblGrid::spacing(), ossimCoarseGridModel::theDlatDhGrid, ossimCoarseGridModel::theDlonDhGrid, ossimCoarseGridModel::theHeightEnabledFlag, ossimSensorModel::theImageSize, ossimCoarseGridModel::theLatGrid, ossimCoarseGridModel::theLonGrid, ossimSensorModel::theSeedFunction, ossimIpt::u, ossimIpt::v, ossimIpt::x, x, and ossimIpt::y.

Referenced by ossimHdf5ProjectionFactory::createProjection(), and ossimViirsHandler::getImageGeometry().

 {
    if (!hdf5)
       return false;
 
    m_hdf5 = hdf5;
    m_projDataPath = projDataPath;
    theHeightEnabledFlag = false;
 
    try
    {
       initCoarseGrid("Latitude",  theLatGrid);
       initCoarseGrid("Longitude", theLonGrid);
    }
    catch (ossimException& x)
    {
       ossimNotify(ossimNotifyLevel_FATAL)<<x.what();
       return false;
    }
 
    theDlatDhGrid.initialize(theLatGrid.size(), theLatGrid.origin(), theLatGrid.spacing(), 0.0);
    theDlonDhGrid.initialize(theLonGrid.size(), theLonGrid.origin(), theLonGrid.spacing(), 0.0);
 
    // Check for dateline crossing among the longitude grid:
    crossesDateline();
 
    ossimGpt ulg (theLatGrid.maxValue(), theLonGrid.minValue());
    ossimGpt urg (theLatGrid.maxValue(), theLonGrid.maxValue());
    ossimGpt lrg (theLatGrid.minValue(), theLonGrid.maxValue());
    ossimGpt llg (theLatGrid.minValue(), theLonGrid.minValue());
    ossimDrect imageRect(0, 0, m_imageSize.x-1, m_imageSize.y-1);
    theSeedFunction = new ossimBilinearProjection(imageRect.ul(), imageRect.ur(),
                                                  imageRect.lr(), imageRect.ll(),
                                                  ulg, urg, lrg, llg);
 
    // Bileaner projection to handle
    initializeModelParams(imageRect);
 
    ossimIrect bounds (0, 0, theImageSize.u-1, theImageSize.v-1);
    initializeModelParams(bounds);
 
    return true;
 }

◆ saveState()

bool ossimHdf5GridModel::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix
	)		const

virtual

Makes sure that the "type" keyword entry reflects the base class, not this one.

Reimplemented from ossimCoarseGridModel.

Definition at line 442 of file ossimHdf5GridModel.cpp.

References ossimKeywordlist::add(), ossimCoarseGridModel::saveState(), and ossimKeywordNames::TYPE_KW.

 {
    bool stat = ossimCoarseGridModel::saveState(kwl, prefix);
    kwl.add(prefix, ossimKeywordNames::TYPE_KW, "ossimCoarseGridModel", true);
 
    return stat;
 }

Member Data Documentation

◆ m_boundGndPolygon

ossimPolygon ossimHdf5GridModel::m_boundGndPolygon

protected

Definition at line 60 of file ossimHdf5GridModel.h.

◆ m_hdf5

ossimRefPtr<ossimHdf5> ossimHdf5GridModel::m_hdf5

protected

Definition at line 62 of file ossimHdf5GridModel.h.

Referenced by initCoarseGrid(), and initialize().

◆ m_imageSize

ossimIpt ossimHdf5GridModel::m_imageSize

protected

Definition at line 63 of file ossimHdf5GridModel.h.

Referenced by initCoarseGrid(), and initialize().

◆ m_projDataPath

ossimString ossimHdf5GridModel::m_projDataPath

protected

Definition at line 64 of file ossimHdf5GridModel.h.

Referenced by initCoarseGrid(), and initialize().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimHdf5GridModel()

◆ ~ossimHdf5GridModel()

Member Function Documentation

◆ crossesDateline()

◆ initCoarseGrid()

◆ initialize()

◆ saveState()

Member Data Documentation

◆ m_boundGndPolygon

◆ m_hdf5

◆ m_imageSize

◆ m_projDataPath