ossimGeoidNgsHeader Class Reference

#include <ossimGeoidNgsHeader.h>

Public Member Functions
	ossimGeoidNgsHeader ()
	ossimGeoidNgsHeader (const ossimFilename &fileName, ossimByteOrder byteOrder=OSSIM_LITTLE_ENDIAN)
bool	initialize (const ossimFilename &fileName, ossimByteOrder byteOrder=OSSIM_LITTLE_ENDIAN)
double	southernMostLat () const
double	westernMostLon () const
double	latDelta () const
double	lonDelta () const
int	rows () const
int	cols () const
int	dataType () const
bool	pointWithin (double lat, double lon) const
const ossimFilename &	filename () const
int	headerSize () const
int	dataTypeSize () const
double	getHeightDelta (double lat, double lon) const

Private Attributes
ossimFilename	theFilename
ossimByteOrder	theByteOrder
double	theSouthernMostLatitude
double	theWesternMostLongitude
double	theLatDelta
double	theLonDelta
int	theRows
int	theCols
int	theDataType

Friends
OSSIM_DLL std::ostream &	operator<< (std::ostream &out, const ossimGeoidNgsHeader &data)

Detailed Description

Definition at line 16 of file ossimGeoidNgsHeader.h.

Constructor & Destructor Documentation

ossimGeoidNgsHeader() [1/2]

ossimGeoidNgsHeader::ossimGeoidNgsHeader

(

)

Definition at line 25 of file ossimGeoidNgsHeader.cpp.

   :theFilename(""),
    theByteOrder(OSSIM_LITTLE_ENDIAN),
    theSouthernMostLatitude(0.0),
    theWesternMostLongitude(0.0),
    theLatDelta(0.0),
    theLonDelta(0.0),
    theRows(0),
    theCols(0),
    theDataType(0)
{   
}

ossimGeoidNgsHeader() [2/2]

ossimGeoidNgsHeader::ossimGeoidNgsHeader	(	const ossimFilename &	fileName,
		ossimByteOrder	byteOrder = OSSIM_LITTLE_ENDIAN
	)

Definition at line 38 of file ossimGeoidNgsHeader.cpp.

References ossim::byteOrder(), and initialize().

{
   initialize(fileName, byteOrder);
}

Member Function Documentation

cols()

int ossimGeoidNgsHeader::cols ( ) const

inline

Definition at line 33 of file ossimGeoidNgsHeader.h.

Referenced by initialize().

{return theCols;}

dataType()

int ossimGeoidNgsHeader::dataType ( ) const

inline

Definition at line 34 of file ossimGeoidNgsHeader.h.

{return theDataType;}

dataTypeSize()

int ossimGeoidNgsHeader::dataTypeSize ( ) const

inline

Definition at line 39 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta().

{return theDataType==1?4:0;}

filename()

const ossimFilename& ossimGeoidNgsHeader::filename ( ) const

inline

Definition at line 36 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta().

{return theFilename;}

getHeightDelta()

double ossimGeoidNgsHeader::getHeightDelta	(	double	lat,
		double	lon
	)		const

Returns

Height delta or ossim::nan() if not found.

Definition at line 108 of file ossimGeoidNgsHeader.cpp.

References dataTypeSize(), filename(), FLT_EPSILON, ossimEndian::getSystemEndianType(), headerSize(), ossim::nan(), ossimNotify(), ossimNotifyLevel_FATAL, ossimEndian::swap(), theByteOrder, theCols, theFilename, theLatDelta, theLonDelta, theSouthernMostLatitude, and theWesternMostLongitude.

{
   float result=0.0;

   // note the headers go from 0 to 360 degrees starting at the prime meridian.
   // ours goes from -180 to 180 degrees.  We will need to shift
   //
   if(lon < 0.0) lon += 360;
   double latSpaces = (lat - theSouthernMostLatitude)/theLatDelta;
   double lonSpaces = (lon - theWesternMostLongitude)/theLonDelta;

   if(latSpaces < -FLT_EPSILON) return ossim::nan();
   if(lonSpaces < -FLT_EPSILON) return ossim::nan();

   long latSpace0 = (long)std::floor(latSpaces);
   long latSpace1 = latSpace0 + 1;
   long lonSpace0 = (long)std::floor(lonSpaces);
   long lonSpace1 = lonSpace0 + 1;

   ifstream input(filename().c_str(), ios::in|ios::binary);
   if(!input)
   {
      ossimNotify(ossimNotifyLevel_FATAL) << "FATAL: " << "ossimGeoidNgsHeader::getGeoidHeight, " << "unable to open file: " << theFilename << "\n";
      return 0;
   }

   long offset00 = (long)(headerSize() + (latSpace0*theCols + lonSpace0)*dataTypeSize());
   long offset01 = (long)(headerSize() + (latSpace0*theCols + lonSpace1)*dataTypeSize());
   long offset11 = (long)(headerSize() + (latSpace1*theCols + lonSpace1)*dataTypeSize());
   long offset10 = (long)(headerSize() + (latSpace1*theCols + lonSpace0)*dataTypeSize());

   double tLat = latSpaces - floor(latSpaces);
   double tLon = lonSpaces - floor(lonSpaces);


   float v00;
   float v01;
   float v11;
   float v10;
   
   input.seekg(offset00,ios::beg);
   input.read((char*)&v00, 4);
   
   input.seekg(offset01,ios::beg);
   input.read((char*)&v01, 4);
   
   input.seekg(offset11,ios::beg);
   input.read((char*)&v11, 4);

   input.seekg(offset10,ios::beg);
   input.read((char*)&v10, 4);
   ossimEndian endian;
   if(endian.getSystemEndianType() != theByteOrder)
   {
      endian.swap(v00);
      endian.swap(v01);
      endian.swap(v11);
      endian.swap(v10);
   }
   double bottom = (v00 + (v01 - v00)*tLat);
   double top    = (v10 + (v11 - v10)*tLat);
   result        = bottom + (top - bottom)*tLon;

   return result;
   
}

headerSize()

int ossimGeoidNgsHeader::headerSize ( ) const

inline

Definition at line 38 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta().

{return 44;}

initialize()

bool ossimGeoidNgsHeader::initialize	(	const ossimFilename &	fileName,
		ossimByteOrder	byteOrder = OSSIM_LITTLE_ENDIAN
	)

Definition at line 44 of file ossimGeoidNgsHeader.cpp.

References ossim::byteOrder(), ossimString::c_str(), cols(), FLT_EPSILON, ossimEndian::getSystemEndianType(), rows(), ossimEndian::swap(), theByteOrder, theCols, theDataType, theFilename, theLatDelta, theLonDelta, theRows, theSouthernMostLatitude, and theWesternMostLongitude.

Referenced by ossimGeoidNgs::addFile(), and ossimGeoidNgsHeader().

{
   ifstream input(fileName.c_str(), ios::in|ios::binary);
   if(!input)
   {
//      cerr << MODULE << " Cannot open geoid file in " << fileName.c_str()
//           << endl;
       return false;
   }

   theByteOrder = byteOrder;
   double latOrigin  = 0.0;
   double lonOrigin  = 0.0;
   double latSpacing = 0.0;
   double lonSpacing = 0.0;
   int    rows       = 0;
   int    cols       = 0;
   int    type       = 0;
   
   input.read((char*)&latOrigin, 8);
   input.read((char*)&lonOrigin, 8);
   input.read((char*)&latSpacing, 8);
   input.read((char*)&lonSpacing, 8);
   input.read((char*)&rows, 4);
   input.read((char*)&cols, 4);
   input.read((char*)&type, 4);

   ossimEndian endian;

   if(endian.getSystemEndianType() != byteOrder)
   {
      endian.swap(latOrigin);
      endian.swap(lonOrigin);
      endian.swap(latSpacing);
      endian.swap(lonSpacing);
      endian.swap(rows);
      endian.swap(cols);
      endian.swap(type);
   }
   
   theFilename = fileName;
   theSouthernMostLatitude = latOrigin;
   theWesternMostLongitude = lonOrigin;
   theLatDelta             = latSpacing;
   theLonDelta             = lonSpacing;
   theRows                 = rows;
   theCols                 = cols;
   theDataType             = type;

   return ((theSouthernMostLatitude >= -FLT_EPSILON)&&
           (theSouthernMostLatitude <= 90.0)&&
           (theWesternMostLongitude >= -FLT_EPSILON)&&
           (theWesternMostLongitude <= 360.0));
}

latDelta()

double ossimGeoidNgsHeader::latDelta ( ) const

inline

Definition at line 30 of file ossimGeoidNgsHeader.h.

{return theLatDelta;}

lonDelta()

double ossimGeoidNgsHeader::lonDelta ( ) const

inline

Definition at line 31 of file ossimGeoidNgsHeader.h.

{return theLonDelta;}

pointWithin()

bool ossimGeoidNgsHeader::pointWithin	(	double	lat,
		double	lon
	)		const

Definition at line 100 of file ossimGeoidNgsHeader.cpp.

References theCols, theLatDelta, theLonDelta, theRows, theSouthernMostLatitude, and theWesternMostLongitude.

{
   return( (lat >= theSouthernMostLatitude) &&
           (lat <= (theSouthernMostLatitude + theRows*theLatDelta)) &&
           (lon >= theWesternMostLongitude) &&
           (lon <= theWesternMostLongitude + theCols*theLonDelta));
}

rows()

int ossimGeoidNgsHeader::rows ( ) const

inline

Definition at line 32 of file ossimGeoidNgsHeader.h.

Referenced by initialize().

{return theRows;}

southernMostLat()

double ossimGeoidNgsHeader::southernMostLat ( ) const

inline

Definition at line 28 of file ossimGeoidNgsHeader.h.

{return theSouthernMostLatitude;}

westernMostLon()

double ossimGeoidNgsHeader::westernMostLon ( ) const

inline

Definition at line 29 of file ossimGeoidNgsHeader.h.

{return theWesternMostLongitude;}

Friends And Related Function Documentation

operator<<

OSSIM_DLL std::ostream& operator<< ( std::ostream &  out, const ossimGeoidNgsHeader &  data  )

friend

Member Data Documentation

theByteOrder

ossimByteOrder ossimGeoidNgsHeader::theByteOrder

private

Definition at line 46 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), and initialize().

theCols

int ossimGeoidNgsHeader::theCols

private

theCols specifies the number of lon increments

Definition at line 74 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), operator<<(), and pointWithin().

theDataType

int ossimGeoidNgsHeader::theDataType

private

Specifies the datatype. if this is 1 then the data that follows is 4 byte floats.

Definition at line 80 of file ossimGeoidNgsHeader.h.

Referenced by initialize(), and operator<<().

theFilename

ossimFilename ossimGeoidNgsHeader::theFilename

private

Definition at line 45 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), and operator<<().

theLatDelta

double ossimGeoidNgsHeader::theLatDelta

private

Specifies the spacing of the latitude direction.

Definition at line 59 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), operator<<(), and pointWithin().

theLonDelta

double ossimGeoidNgsHeader::theLonDelta

private

longitude spacing.

Definition at line 64 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), operator<<(), and pointWithin().

theRows

int ossimGeoidNgsHeader::theRows

private

theRows specifies the number of lat increments

Definition at line 69 of file ossimGeoidNgsHeader.h.

Referenced by initialize(), operator<<(), and pointWithin().

theSouthernMostLatitude

double ossimGeoidNgsHeader::theSouthernMostLatitude

private

Definition at line 47 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), operator<<(), and pointWithin().

theWesternMostLongitude

double ossimGeoidNgsHeader::theWesternMostLongitude

private

Important: The geoid.bin files expresses this as a possitive value.

example: -90.0 is 270 degrees.

Definition at line 54 of file ossimGeoidNgsHeader.h.

Referenced by getHeightDelta(), initialize(), operator<<(), and pointWithin().

---

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

• ossimGeoidNgsHeader.h
• ossimGeoidNgsHeader.cpp