ossimDemHeader Class Reference

#include <ossimDemHeader.h>

Public Member Functions
	ossimDemHeader ()
bool	open (const ossimFilename &file)
	open method that takes a file. More...
bool	open (std::shared_ptr< ossim::istream > &str, const std::string &connectionString)
ossim::istream &	open (ossim::istream &is)
	open method that takes a stream. More...
std::ostream &	print (std::ostream &out) const
	Print method. More...
std::string const &	getQuadName () const
std::string const &	getProcessInfo () const
double	getSEGeoCornerX () const
double	getSEGeoCornerY () const
ossim_int32	getProcessCode () const
std::string const &	getSectionIndicator () const
std::string const &	getMappingCenterCode () const
ossim_int32	getLevelCode () const
ossim_int32	getElevPattern () const
ossim_int32	getGroundRefSysCode () const
ossimString	getProjectionType () const
ossim_int32	getGroundRefSysZone () const
ossim_int32	getGroundRefSysUnits () const
ossim_int32	getElevationUnits () const
ossim_int32	getNumPolySides () const
const ossimDemPointVector &	getDEMCorners () const
double	getCounterclockAngle () const
ossim_int32	getElevAccuracyCode () const
double	getMinimumElev () const
double	getMaximumElev () const
double	getSpatialResX () const
double	getSpatialResY () const
double	getSpatialResZ () const
ossim_int32	getProfileRows () const
ossim_int32	getProfileColumns () const
ossim_int32	getLargeContInt () const
ossim_int32	getMaxSourceUnits () const
ossim_int32	getSmallContInt () const
ossim_int32	getMinSourceUnits () const
ossim_int32	getSourceDate () const
ossim_int32	getInspRevDate () const
std::string const &	getInspFlag () const
ossim_int32	getDataValFlag () const
ossim_int32	getSuspectVoidFlag () const
ossim_int32	getVertDatum () const
ossim_int32	getHorizDatum () const
ossimString	getHorizDatumCode () const
ossim_int32	getDataEdition () const
ossim_int32	getPerctVoid () const
ossim_int32	getWestEdgeFlag () const
ossim_int32	getNorthEdgeFlag () const
ossim_int32	getEastEdgeFlag () const
ossim_int32	getSouthEdgeFlag () const
double	getVertDatumShift () const
bool	getImageGeometry (ossimKeywordlist &kwl, const char *prefix=NULL) const

Private Attributes
std::string	_quadName
std::string	_processInfo
double	_seGeoCornerX
double	_seGeoCornerY
ossim_int32	_processCode
std::string	_sectionIndicator
std::string	_mapCenterCode
ossim_int32	_levelCode
ossim_int32	_elevPattern
ossim_int32	_groundRefSysCode
ossim_int32	_groundRefSysZone
ossim_int32	_groundRefSysUnits
ossim_int32	_elevUnits
ossim_int32	_numPolySides
ossimDemPointVector	_demCorners
double	_counterclockAngle
ossim_int32	_elevAccuracyCode
double	_minElevation
double	_maxElevation
double	_spatialResX
double	_spatialResY
double	_spatialResZ
ossim_int32	_profileRows
ossim_int32	_profileColumns
ossim_int32	_largeContInt
ossim_int32	_maxSourceUnits
ossim_int32	_smallContInt
ossim_int32	_minSourceUnits
ossim_int32	_sourceDate
ossim_int32	_inspRevDate
std::string	_inspFlag
ossim_int32	_valFlag
ossim_int32	_suspectVoidFlg
ossim_int32	_vertDatum
ossim_int32	_horizDatum
ossim_int32	_dataEdition
ossim_int32	_perctVoid
ossim_int32	_westEdgeFlag
ossim_int32	_northEdgeFlag
ossim_int32	_eastEdgeFlag
ossim_int32	_southEdgeFlag
double	_vertDatumShift

Friends
OSSIM_DLL std::ostream &	operator<< (ossim::ostream &s, const ossimDemHeader &header)
OSSIM_DLL ossim::istream &	operator>> (ossim::istream &s, ossimDemHeader &header)

Detailed Description

Definition at line 30 of file ossimDemHeader.h.

Constructor & Destructor Documentation

ossimDemHeader()

ossimDemHeader::ossimDemHeader

(

)

Definition at line 78 of file ossimDemHeader.cpp.

   : _seGeoCornerX(0.0),
     _seGeoCornerY(0.0),
     _processCode(0),
     _levelCode(0),
     _elevPattern(0),
     _groundRefSysCode(0),
     _groundRefSysZone(0),
     _groundRefSysUnits(0),
     _elevUnits(0),
     _numPolySides(0),
     _counterclockAngle(0.0),
     _elevAccuracyCode(0),
     _minElevation(0),
     _maxElevation(0),
     _spatialResX(0.0),
     _spatialResY(0.0),
     _spatialResZ(0.0),
     _profileRows(0),
     _profileColumns(0),
     _largeContInt(0),
     _maxSourceUnits(0),
     _smallContInt(0),
     _minSourceUnits(0),
     _sourceDate(0),
     _inspRevDate(0),
     _valFlag(0),
     _suspectVoidFlg(0),
     _vertDatum(0),
     _horizDatum(0),
     _dataEdition(0),
     _perctVoid(0),
     _westEdgeFlag(0),
     _northEdgeFlag(0),
     _eastEdgeFlag(0),
     _southEdgeFlag(0),
     _vertDatumShift(0.0)
{
}

Member Function Documentation

getCounterclockAngle()

double ossimDemHeader::getCounterclockAngle

(

)

const

In Radians the counterclockwise angle from the primary axis of ground planimetric reference to the primary axis of the DEM local reference system. Set to 0 to align with the coordinate system

Definition at line 209 of file ossimDemHeader.cpp.

References _counterclockAngle.

{
   return _counterclockAngle;
}

getDataEdition()

ossim_int32 ossimDemHeader::getDataEdition

(

)

const

01-99 Primarily a DMA specfic field. (For USGS use set to 01)

Definition at line 329 of file ossimDemHeader.cpp.

References _dataEdition.

{
   return _dataEdition;
}

getDataValFlag()

ossim_int32 ossimDemHeader::getDataValFlag

(

)

const

0 == No Validation 1 == RSME computed from test points, no quantative test or interactive editing/review 2 == Batch process water body edit and RSME computed from test points 3 == Review and edit, including water edit. No RSME computed from test points 4 == Level 1 DEM's reviewed & edited. RSME computed from test points 5 == Level 2 & 3 DEM's reviewed & edited. RSME computed from test points

Definition at line 305 of file ossimDemHeader.cpp.

References _valFlag.

{
   return _valFlag;
}

getDEMCorners()

ossimDemPointVector const & ossimDemHeader::getDEMCorners

(

)

const

The coordinates of the quadrangle corners are ordered in a clockwise direction beginning with the southwest corner.

Definition at line 203 of file ossimDemHeader.cpp.

References _demCorners.

{
   return _demCorners;
}

getEastEdgeFlag()

ossim_int32 ossimDemHeader::getEastEdgeFlag

(

)

const

East edge match staus flag

Definition at line 353 of file ossimDemHeader.cpp.

References _eastEdgeFlag.

{
   return _eastEdgeFlag;
}

getElevAccuracyCode()

ossim_int32 ossimDemHeader::getElevAccuracyCode

(

)

const

0 == unknown accuracy 1 == accuracy info is given in logical record type C

Definition at line 215 of file ossimDemHeader.cpp.

References _elevAccuracyCode.

{
   return _elevAccuracyCode;
}

getElevationUnits()

ossim_int32 ossimDemHeader::getElevationUnits

(

)

const

Code defining unit of measure for elevation coordinates used throughout the file. 1 == Feet 2 == Meters

Definition at line 191 of file ossimDemHeader.cpp.

References _elevUnits.

{
   return _elevUnits;
}

getElevPattern()

ossim_int32 ossimDemHeader::getElevPattern

(

)

const

1 == Regular 2 == Random (reserved for future use).

Definition at line 167 of file ossimDemHeader.cpp.

References _elevPattern.

{
   return _elevPattern;
}

getGroundRefSysCode()

ossim_int32 ossimDemHeader::getGroundRefSysCode

(

)

const

Code defining ground planimetric reference system. 0 == Geographic 1 == UTM 2 == State Plane

Definition at line 173 of file ossimDemHeader.cpp.

References _groundRefSysCode.

Referenced by getImageGeometry(), getProjectionType(), and ossimDemGrid::read().

{
   return _groundRefSysCode;
}

getGroundRefSysUnits()

ossim_int32 ossimDemHeader::getGroundRefSysUnits

(

)

const

Code defining unit of measure for ground planimetric coordinates used throughout the file. 0 == Radians 1 == Feet 2 == Meters 3 == Arc Seconds

Definition at line 185 of file ossimDemHeader.cpp.

References _groundRefSysUnits.

Referenced by getImageGeometry().

{
   return _groundRefSysUnits;
}

getGroundRefSysZone()

ossim_int32 ossimDemHeader::getGroundRefSysZone

(

)

const

Code defining zone in ground planimetric reference system.

Definition at line 179 of file ossimDemHeader.cpp.

References _groundRefSysZone.

Referenced by getImageGeometry().

{
   return _groundRefSysZone;
}

getHorizDatum()

ossim_int32 ossimDemHeader::getHorizDatum

(

)

const

1 == North Amercian Datum 1927 (NAD 27) 2 == World Geoditic System 1972 (WGS 72) 3 == WGS 84 4 == NAD 83 5 == Old Hawaii Datum 6 == Puerto Rico Datum

Definition at line 323 of file ossimDemHeader.cpp.

References _horizDatum.

Referenced by getHorizDatumCode().

{
   return _horizDatum;
}

getHorizDatumCode()

ossimString ossimDemHeader::getHorizDatumCode

(

)

const

Returns

The ossim datum code as a string.

Definition at line 842 of file ossimDemHeader.cpp.

References getHorizDatum().

Referenced by getImageGeometry().

{
   ossimString code;
   switch(getHorizDatum())
   {
      case 1:
         code = "NAS-C";
         break;
      case 2:
         code = "WGD";
         break;
      case 3:
         code = "WGE";
        break;         
      case 4:
         code = "NAR-C";
        break;         
      case 5:
         code = "OHA-M";
        break;         
      case 6:
         code = "PUR";
         break;
      default:
         code = "unknown";
        break;         
   }
   return code;
}

getImageGeometry()

bool ossimDemHeader::getImageGeometry	(	ossimKeywordlist &	kwl,
		const char *	prefix = NULL
	)		const

Populates keyword list with geometry info.

Returns

true on success, false on error.

Definition at line 739 of file ossimDemHeader.cpp.

References _demCorners, abs, ossimKeywordlist::add(), ossimKeywordNames::DATUM_KW, ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LAT, ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LON, getGroundRefSysCode(), getGroundRefSysUnits(), getGroundRefSysZone(), getHorizDatumCode(), getProjectionType(), getSpatialResX(), ossimKeywordNames::HEMISPHERE_KW, M_PI, ossimKeywordNames::METERS_PER_PIXEL_X_KW, ossimKeywordNames::METERS_PER_PIXEL_Y_KW, ossimKeywordNames::TIE_POINT_EASTING_KW, ossimKeywordNames::TIE_POINT_LAT_KW, ossimKeywordNames::TIE_POINT_LON_KW, ossimKeywordNames::TIE_POINT_NORTHING_KW, ossimKeywordNames::TYPE_KW, ossim::usft2mtrs(), and ossimKeywordNames::ZONE_KW.

{
   ossimString type = getProjectionType();
   ossimString datum = getHorizDatumCode();
   if ( (type == "unknown") || (datum == "unknown") )
   {
      return false;
   }

   if (_demCorners.size() != 4)
   {
      return false;
   }

   double tieX   = _demCorners[1].getX();
   double tieY   = _demCorners[1].getY();
   double scaleX = getSpatialResX();
   double scaleY = getSpatialResX();
   
   // Add these for all projections.
   kwl.add(prefix,
           ossimKeywordNames::TYPE_KW,
           type);
   
   kwl.add(prefix,
           ossimKeywordNames::DATUM_KW,
           datum);

   if (getGroundRefSysCode() == 0) // Geographic.
   {
      // ESH 10/2008 -- The Dem ground units can be either radians or 
      // arc seconds, so we have to convert parameters in these units 
      // to degrees which is what OSSIM is assuming.
      bool bIsArcSecs = (getGroundRefSysUnits() == 3) ? true : false;
      bool bIsRadians = (getGroundRefSysUnits() == 0) ? true : false;

      double convertFactor = 1.0;
      if ( bIsArcSecs == true )
      {
         convertFactor = 1.0 / 3600;
      }
      else if ( bIsRadians == true )
      {
         convertFactor = 180.0 / M_PI;
      }

      kwl.add(prefix,
              ossimKeywordNames::TIE_POINT_LON_KW,
              (tieX * convertFactor) );
      kwl.add(prefix,
              ossimKeywordNames::TIE_POINT_LAT_KW,
              (tieY * convertFactor) );
      kwl.add(prefix,
              ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LON,
              (scaleX * convertFactor) );
      kwl.add(prefix,
              ossimKeywordNames::DECIMAL_DEGREES_PER_PIXEL_LAT,
              (scaleY * convertFactor) );
   }
   else if (getGroundRefSysCode() == 1) // UTM
   {
      // Get the zone.
      kwl.add(prefix,
              ossimKeywordNames::ZONE_KW,
              abs(getGroundRefSysZone()));

      // Set the hemisphere.
      kwl.add(prefix,
              ossimKeywordNames::HEMISPHERE_KW,
              (getGroundRefSysZone()<0?"S":"N"));

      // Set the tie.
      if (getGroundRefSysUnits() == 1) // Feet...
      {
         tieX = ossim::usft2mtrs(tieX);
         tieY = ossim::usft2mtrs(tieY);
         scaleX = ossim::usft2mtrs(scaleX);
         scaleY = ossim::usft2mtrs(scaleY);
         
      }
      kwl.add(prefix,
              ossimKeywordNames::TIE_POINT_EASTING_KW,
              tieX);
      kwl.add(prefix,
              ossimKeywordNames::TIE_POINT_NORTHING_KW,
              tieY);
      kwl.add(prefix,
              ossimKeywordNames::METERS_PER_PIXEL_X_KW,
              scaleX);
      kwl.add(prefix,
              ossimKeywordNames::METERS_PER_PIXEL_Y_KW,
              scaleY);
   }
   else if (getGroundRefSysCode() == 2) // State Plane
   {
      return false; // Need to code...
   }

   
   return true;
}

getInspFlag()

std::string const & ossimDemHeader::getInspFlag

(

)

const

"I" indicates all processes of part3, Quality Control have been performed.

Definition at line 299 of file ossimDemHeader.cpp.

References _inspFlag.

{
   return _inspFlag;
}

getInspRevDate()

ossim_int32 ossimDemHeader::getInspRevDate

(

)

const

YYY

Definition at line 293 of file ossimDemHeader.cpp.

References _inspRevDate.

{
   return _inspRevDate;
}

getLargeContInt()

ossim_int32 ossimDemHeader::getLargeContInt

(

)

const

Present only if two or more primary intervals exist (level 2 DEM's only)

Definition at line 263 of file ossimDemHeader.cpp.

References _largeContInt.

{
   return _largeContInt;
}

getLevelCode()

ossim_int32 ossimDemHeader::getLevelCode

(

)

const

1 == DEM-1 2 == DEM-2 3 == DEM-3

Definition at line 161 of file ossimDemHeader.cpp.

References _levelCode.

{
   return _levelCode;
}

getMappingCenterCode()

std::string const & ossimDemHeader::getMappingCenterCode

(

)

const

Mapping Center origin code. "EMC","WMC","MCMC","RMMC","FS" and "GPM2"

Definition at line 155 of file ossimDemHeader.cpp.

References _mapCenterCode.

{
   return _mapCenterCode;
}

getMaximumElev()

double ossimDemHeader::getMaximumElev

(

)

const

Maximum elevation in the DEM.

Definition at line 227 of file ossimDemHeader.cpp.

References _maxElevation.

{
   return _maxElevation;
}

getMaxSourceUnits()

ossim_int32 ossimDemHeader::getMaxSourceUnits

(

)

const

0 == NA 1 == feet 2 == meters (level two DEM's only)

Definition at line 269 of file ossimDemHeader.cpp.

References _maxSourceUnits.

{
   return _maxSourceUnits;
}

getMinimumElev()

double ossimDemHeader::getMinimumElev

(

)

const

Minimum elevation in the DEM.

Definition at line 221 of file ossimDemHeader.cpp.

References _minElevation.

{
   return _minElevation;
}

getMinSourceUnits()

ossim_int32 ossimDemHeader::getMinSourceUnits

(

)

const

1 == feet 2 == meters

Definition at line 281 of file ossimDemHeader.cpp.

References _minSourceUnits.

{
   return _minSourceUnits;
}

getNorthEdgeFlag()

ossim_int32 ossimDemHeader::getNorthEdgeFlag

(

)

const

North edge match status flag

Definition at line 347 of file ossimDemHeader.cpp.

References _northEdgeFlag.

{
   return _northEdgeFlag;
}

getNumPolySides()

ossim_int32 ossimDemHeader::getNumPolySides

(

)

const

Set n == 4

Definition at line 197 of file ossimDemHeader.cpp.

References _numPolySides.

{
   return _numPolySides;
}

getPerctVoid()

ossim_int32 ossimDemHeader::getPerctVoid

(

)

const

Contains the percentage of nodes in the file set to void

Definition at line 335 of file ossimDemHeader.cpp.

References _perctVoid.

{
   return _perctVoid;
}

getProcessCode()

ossim_int32 ossimDemHeader::getProcessCode

(

)

const

1 == Autocorrelation resample simple bilinear 2 == Manual profile GRIDEM simple bilinear 3 == DLG/hypsography CTOG 8-direction bilinear 4 == Interpolation from photogrammetic system contours DCASS 4-direction bilinear 5 == DLG/hypsography LINETRACE, LT4X complex linear 6 == DLG/hypsography CPS-3, ANUDEM, GRASS complex polynomial 7 == Electronic imaging (non-photogrammetric), active or passive sensor systems

Definition at line 143 of file ossimDemHeader.cpp.

References _processCode.

Referenced by print().

{
   return _processCode;
}

getProcessInfo()

std::string const & ossimDemHeader::getProcessInfo

(

)

const

Free format desciptor field containing process information.

Definition at line 125 of file ossimDemHeader.cpp.

References _processInfo.

Referenced by print().

{
   return _processInfo;
}

getProfileColumns()

ossim_int32 ossimDemHeader::getProfileColumns

(

)

const

Definition at line 257 of file ossimDemHeader.cpp.

References _profileColumns.

Referenced by ossimDemGrid::fillGeographic(), and ossimDemGrid::fillUTM().

{
   return _profileColumns;
}

getProfileRows()

ossim_int32 ossimDemHeader::getProfileRows

(

)

const

The number of rows and columns in profiles in the DEM. NOTE: When profile rows == 1, profile columns is the number of profiles in the DEM.

Definition at line 251 of file ossimDemHeader.cpp.

References _profileRows.

{
   return _profileRows;
}

getProjectionType()

ossimString ossimDemHeader::getProjectionType

(

)

const

Returns

The projection type as a string.

Definition at line 872 of file ossimDemHeader.cpp.

References getGroundRefSysCode().

Referenced by getImageGeometry().

{
   ossimString type;
   switch(getGroundRefSysCode())
   {
      case 0:
         type = "ossimEquDistCylProjection";
         break;
      case 1:
         type = "ossimUtmProjection";
         break;
      case 2:
      default:
         type = "unknown";
         break;
   }
   return type;
}

getQuadName()

std::string const & ossimDemHeader::getQuadName

(

)

const

Quadrangle name.

Definition at line 119 of file ossimDemHeader.cpp.

References _quadName.

Referenced by print().

{
   return _quadName;
}

getSectionIndicator()

std::string const & ossimDemHeader::getSectionIndicator

(

)

const

Used for 30 minute DEMs - Identifies 1:100,000 scale sections Formatted as XNN, wher X is "S" == 7.5 min. & "F" == 15 min. NN is a two digit sequence number

Definition at line 149 of file ossimDemHeader.cpp.

References _sectionIndicator.

{
   return _sectionIndicator;
}

getSEGeoCornerX()

double ossimDemHeader::getSEGeoCornerX

(

)

const

SE geographic corner orderd as: x == Longitude == SDDDMMSS.SSSS (neg sign (s) right justified, no leading zeros, plus sign (s) implied

Definition at line 131 of file ossimDemHeader.cpp.

References _seGeoCornerX.

Referenced by print().

{
   return _seGeoCornerX;
}

getSEGeoCornerY()

double ossimDemHeader::getSEGeoCornerY

(

)

const

SE geographic corner orderd as: y == Latitude == SDDMMSS.SSSS (neg sign (s) right justified, no leading zeros, plus sign (s) implied

Definition at line 137 of file ossimDemHeader.cpp.

References _seGeoCornerY.

Referenced by print().

{
   return _seGeoCornerY;
}

getSmallContInt()

ossim_int32 ossimDemHeader::getSmallContInt

(

)

const

Smallest or only primary contour interval (level 2 DEM's only)

Definition at line 275 of file ossimDemHeader.cpp.

References _smallContInt.

{
   return _smallContInt;
}

getSourceDate()

ossim_int32 ossimDemHeader::getSourceDate

(

)

const

YYY

Definition at line 287 of file ossimDemHeader.cpp.

References _sourceDate.

{
   return _sourceDate;
}

getSouthEdgeFlag()

ossim_int32 ossimDemHeader::getSouthEdgeFlag

(

)

const

South edge match status flag

Definition at line 359 of file ossimDemHeader.cpp.

References _southEdgeFlag.

{
   return _southEdgeFlag;
}

getSpatialResX()

double ossimDemHeader::getSpatialResX

(

)

const

Definition at line 233 of file ossimDemHeader.cpp.

References _spatialResX.

Referenced by ossimDemGrid::getGroundCoords(), and getImageGeometry().

{
   return _spatialResX;
}

getSpatialResY()

double ossimDemHeader::getSpatialResY

(

)

const

Definition at line 239 of file ossimDemHeader.cpp.

References _spatialResY.

Referenced by ossimDemGrid::fillGeographic(), ossimDemGrid::fillUTM(), and ossimDemGrid::getGroundCoords().

{
   return _spatialResY;
}

getSpatialResZ()

double ossimDemHeader::getSpatialResZ

(

)

const

Definition at line 245 of file ossimDemHeader.cpp.

References _spatialResZ.

Referenced by ossimUsgsDemTileSource::fillBuffer().

{
   return _spatialResZ;
}

getSuspectVoidFlag()

ossim_int32 ossimDemHeader::getSuspectVoidFlag

(

)

const

0 == none 1 == suspect areas 2 == void areas 3 == suspect & void areas

Definition at line 311 of file ossimDemHeader.cpp.

References _suspectVoidFlg.

{
   return _suspectVoidFlg;
}

getVertDatum()

ossim_int32 ossimDemHeader::getVertDatum

(

)

const

1 == local mean sea level 2 == National Geodetic Vertical Datum 1929 3 == North American Vertical Datum 1988

Definition at line 317 of file ossimDemHeader.cpp.

References _vertDatum.

{
   return _vertDatum;
}

getVertDatumShift()

double ossimDemHeader::getVertDatumShift

(

)

const

Vertical datum shift - value is in the form of SFFF.DD Value is the average shift value for the four quadrangle corners obtained from program VERTCOM.

Definition at line 365 of file ossimDemHeader.cpp.

References _vertDatumShift.

{
   return _vertDatumShift;
}

getWestEdgeFlag()

ossim_int32 ossimDemHeader::getWestEdgeFlag

(

)

const

West edge match status flag

Definition at line 341 of file ossimDemHeader.cpp.

References _westEdgeFlag.

{
   return _westEdgeFlag;
}

open() [1/3]

bool ossimDemHeader::open

(

const ossimFilename &

file

)

open method that takes a file.

Parameters

file	File name to open.

Returns

true on success false on error.

Definition at line 376 of file ossimDemHeader.cpp.

References ossimString::c_str(), ossim::StreamFactoryRegistry::instance(), and ossimDemUtil::isUsgsDem().

Referenced by operator>>(), and ossimDemInfo::print().

{
   bool result = ossimDemUtil::isUsgsDem(file);
   std::string connectionString = file.c_str();
   std::shared_ptr<ossim::istream> is = ossim::StreamFactoryRegistry::instance()->
      createIstream(connectionString, std::ios_base::in | std::ios_base::binary);
   if(result)
   {
      if ( is && is->good() )
      {
         open( is, connectionString );
      }
      else
      {
         result = false;
      }    
   }

   return result;
}

open() [2/3]

bool ossimDemHeader::open	(	std::shared_ptr< ossim::istream > &	str,
		const std::string &	connectionString
	)

Definition at line 399 of file ossimDemHeader.cpp.

References _counterclockAngle, _dataEdition, _demCorners, _eastEdgeFlag, _elevAccuracyCode, _elevPattern, _elevUnits, _groundRefSysCode, _groundRefSysUnits, _groundRefSysZone, _horizDatum, _inspFlag, _inspRevDate, _largeContInt, _levelCode, _mapCenterCode, _maxElevation, _maxSourceUnits, _minElevation, _minSourceUnits, _northEdgeFlag, _numPolySides, _perctVoid, _processCode, _processInfo, _profileColumns, _profileRows, _quadName, _sectionIndicator, _seGeoCornerX, _seGeoCornerY, _smallContInt, _sourceDate, _southEdgeFlag, _spatialResX, _spatialResY, _spatialResZ, _suspectVoidFlg, _valFlag, _vertDatum, _vertDatumShift, _westEdgeFlag, ossimDemUtil::getDouble(), ossimDemUtil::getLong(), ossimDemUtil::getRecord(), ossimDemUtil::isUsgsDem(), x, and y.

{
   bool result = ossimDemUtil::isUsgsDem(ossimFilename(connectionString));
   if ( result&&str&& str->good() )
   {

      char* bufstr = new char[1024];
      char* temp   = new char[1024];
      ossim_int32 i;
      
      ossimDemUtil::getRecord(*str, bufstr);
      
      strncpy(temp, bufstr, 40);
      temp[40] = '\0';
      _quadName = temp;
      
      strncpy(temp,bufstr+40,40);
      temp[40] = '\0';
      _processInfo = temp;
      
      ossimDemUtil::getDouble(bufstr, 109, 13, _seGeoCornerX);
      ossimDemUtil::getDouble(bufstr, 122, 13, _seGeoCornerY);
      _processCode = ossimDemUtil::getLong(bufstr, 135, 1);
      
      strncpy(temp,bufstr+137,3);
      temp[3] = '\0';
      _sectionIndicator = temp;
      
      strncpy(temp,bufstr+140,4);
      temp[4] = '\0';
      _mapCenterCode = temp;
      
      _levelCode = ossimDemUtil::getLong(bufstr, 144, 6);
      _elevPattern = ossimDemUtil::getLong(bufstr, 150, 6);
      _groundRefSysCode = ossimDemUtil::getLong(bufstr, 156, 6);
      _groundRefSysZone = ossimDemUtil::getLong(bufstr, 162, 6);
      _groundRefSysUnits = ossimDemUtil::getLong(bufstr, 528, 6);
      _elevUnits = ossimDemUtil::getLong(bufstr, 534, 6);
      _numPolySides = ossimDemUtil::getLong(bufstr, 540, 6);
      
      for (i = 0; i < 4; i++)
      {
         double x,y;
         ossim_int32 pos = 546 + (i * 48);
         ossimDemUtil::getDouble(bufstr, pos, 24, x);
         ossimDemUtil::getDouble(bufstr, pos + 24, 24, y);
         _demCorners.push_back(ossimDemPoint(x,y));
      }
      
      ossimDemUtil::getDouble(bufstr, 738, 24, _minElevation);
      ossimDemUtil::getDouble(bufstr, 762, 24, _maxElevation);
      ossimDemUtil::getDouble(bufstr, 786, 24, _counterclockAngle );
      _elevAccuracyCode = ossimDemUtil::getLong(bufstr, 810, 6);
      ossimDemUtil::getDouble(bufstr, 816, 12, _spatialResX);
      ossimDemUtil::getDouble(bufstr, 828, 12, _spatialResY);
      ossimDemUtil::getDouble(bufstr, 840, 12, _spatialResZ);
      _profileRows = ossimDemUtil::getLong(bufstr, 852, 6);
      _profileColumns = ossimDemUtil::getLong(bufstr, 858, 6);
      _largeContInt = ossimDemUtil::getLong(bufstr, 864, 5);
      _maxSourceUnits = ossimDemUtil::getLong(bufstr, 869, 1);
      _smallContInt = ossimDemUtil::getLong(bufstr, 870, 5);
      _minSourceUnits = ossimDemUtil::getLong(bufstr, 875, 1);
      _sourceDate = ossimDemUtil::getLong(bufstr, 876, 4);
      _inspRevDate = ossimDemUtil::getLong(bufstr, 880, 4);
      
      strncpy(temp, bufstr+884,1);
      temp[1]='\0';
      _inspFlag = temp;
      
      _valFlag = ossimDemUtil::getLong(bufstr, 885, 1);
      _suspectVoidFlg = ossimDemUtil::getLong(bufstr, 886, 2);
      _vertDatum = ossimDemUtil::getLong(bufstr, 888, 2);
      _horizDatum = ossimDemUtil::getLong(bufstr, 890, 2);
      if (_horizDatum == 0)
         _horizDatum = 1;   // Default to NAD27
      
      _dataEdition = ossimDemUtil::getLong(bufstr, 892, 4);
      _perctVoid = ossimDemUtil::getLong(bufstr, 896, 4);
      _westEdgeFlag = ossimDemUtil::getLong(bufstr, 900, 2);
      _northEdgeFlag = ossimDemUtil::getLong(bufstr, 902, 2);
      _eastEdgeFlag = ossimDemUtil::getLong(bufstr, 904, 2);
      _southEdgeFlag = ossimDemUtil::getLong(bufstr, 906, 2);
      ossimDemUtil::getDouble(bufstr, 908, 7, _vertDatumShift);

      delete [] bufstr;
      delete [] temp;
      bufstr = 0;
      temp = 0;
   }
   return result;
}

open() [3/3]

std::istream & ossimDemHeader::open

(

ossim::istream &

is

)

open method that takes a stream.

Parameters

is	The input stream to read from.

Returns

stream

Definition at line 493 of file ossimDemHeader.cpp.

References _counterclockAngle, _dataEdition, _demCorners, _eastEdgeFlag, _elevAccuracyCode, _elevPattern, _elevUnits, _groundRefSysCode, _groundRefSysUnits, _groundRefSysZone, _horizDatum, _inspFlag, _inspRevDate, _largeContInt, _levelCode, _mapCenterCode, _maxElevation, _maxSourceUnits, _minElevation, _minSourceUnits, _northEdgeFlag, _numPolySides, _perctVoid, _processCode, _processInfo, _profileColumns, _profileRows, _quadName, _sectionIndicator, _seGeoCornerX, _seGeoCornerY, _smallContInt, _sourceDate, _southEdgeFlag, _spatialResX, _spatialResY, _spatialResZ, _suspectVoidFlg, _valFlag, _vertDatum, _vertDatumShift, _westEdgeFlag, ossimDemUtil::getDouble(), ossimDemUtil::getLong(), ossimDemUtil::getRecord(), x, and y.

{
   if(in.good())
   {
      char* bufstr = new char[1024];
      char* temp   = new char[1024];
      ossim_int32 i;
      
      ossimDemUtil::getRecord(in, bufstr);
      
      strncpy(temp, bufstr, 40);
      temp[40] = '\0';
      _quadName = temp;
      
      strncpy(temp,bufstr+40,40);
      temp[40] = '\0';
      _processInfo = temp;
      
      ossimDemUtil::getDouble(bufstr, 109, 13, _seGeoCornerX);
      ossimDemUtil::getDouble(bufstr, 122, 13, _seGeoCornerY);
      _processCode = ossimDemUtil::getLong(bufstr, 135, 1);
      
      strncpy(temp,bufstr+137,3);
      temp[3] = '\0';
      _sectionIndicator = temp;
      
      strncpy(temp,bufstr+140,4);
      temp[4] = '\0';
      _mapCenterCode = temp;
      
      _levelCode = ossimDemUtil::getLong(bufstr, 144, 6);
      _elevPattern = ossimDemUtil::getLong(bufstr, 150, 6);
      _groundRefSysCode = ossimDemUtil::getLong(bufstr, 156, 6);
      _groundRefSysZone = ossimDemUtil::getLong(bufstr, 162, 6);
      _groundRefSysUnits = ossimDemUtil::getLong(bufstr, 528, 6);
      _elevUnits = ossimDemUtil::getLong(bufstr, 534, 6);
      _numPolySides = ossimDemUtil::getLong(bufstr, 540, 6);
      
      for (i = 0; i < 4; i++)
      {
         double x,y;
         ossim_int32 pos = 546 + (i * 48);
         ossimDemUtil::getDouble(bufstr, pos, 24, x);
         ossimDemUtil::getDouble(bufstr, pos + 24, 24, y);
         _demCorners.push_back(ossimDemPoint(x,y));
      }
      
      ossimDemUtil::getDouble(bufstr, 738, 24, _minElevation);
      ossimDemUtil::getDouble(bufstr, 762, 24, _maxElevation);
      ossimDemUtil::getDouble(bufstr, 786, 24, _counterclockAngle );
      _elevAccuracyCode = ossimDemUtil::getLong(bufstr, 810, 6);
      ossimDemUtil::getDouble(bufstr, 816, 12, _spatialResX);
      ossimDemUtil::getDouble(bufstr, 828, 12, _spatialResY);
      ossimDemUtil::getDouble(bufstr, 840, 12, _spatialResZ);
      _profileRows = ossimDemUtil::getLong(bufstr, 852, 6);
      _profileColumns = ossimDemUtil::getLong(bufstr, 858, 6);
      _largeContInt = ossimDemUtil::getLong(bufstr, 864, 5);
      _maxSourceUnits = ossimDemUtil::getLong(bufstr, 869, 1);
      _smallContInt = ossimDemUtil::getLong(bufstr, 870, 5);
      _minSourceUnits = ossimDemUtil::getLong(bufstr, 875, 1);
      _sourceDate = ossimDemUtil::getLong(bufstr, 876, 4);
      _inspRevDate = ossimDemUtil::getLong(bufstr, 880, 4);
      
      strncpy(temp, bufstr+884,1);
      temp[1]='\0';
      _inspFlag = temp;
      
      _valFlag = ossimDemUtil::getLong(bufstr, 885, 1);
      _suspectVoidFlg = ossimDemUtil::getLong(bufstr, 886, 2);
      _vertDatum = ossimDemUtil::getLong(bufstr, 888, 2);
      _horizDatum = ossimDemUtil::getLong(bufstr, 890, 2);
      if (_horizDatum == 0)
         _horizDatum = 1;   // Default to NAD27
      
      _dataEdition = ossimDemUtil::getLong(bufstr, 892, 4);
      _perctVoid = ossimDemUtil::getLong(bufstr, 896, 4);
      _westEdgeFlag = ossimDemUtil::getLong(bufstr, 900, 2);
      _northEdgeFlag = ossimDemUtil::getLong(bufstr, 902, 2);
      _eastEdgeFlag = ossimDemUtil::getLong(bufstr, 904, 2);
      _southEdgeFlag = ossimDemUtil::getLong(bufstr, 906, 2);
      ossimDemUtil::getDouble(bufstr, 908, 7, _vertDatumShift);

      delete [] bufstr;
      delete [] temp;
      bufstr = 0;
      temp = 0;
   }
   return in;
}

print()

std::ostream & ossimDemHeader::print

(

std::ostream &

out

)

const

Print method.

Parameters

out	Stream to print to.

Returns

std::ostream&

Definition at line 584 of file ossimDemHeader.cpp.

References getProcessCode(), getProcessInfo(), getQuadName(), getSEGeoCornerX(), and getSEGeoCornerY().

Referenced by operator<<(), and ossimDemInfo::print().

{
   const int W = 24; // format width
   const int CW = W-8; // format corner width
   
   // Capture the original flags then set float output to full precision.
   std::ios_base::fmtflags f = out.flags();

   // Note:  This is only a partial print for now...
   std::string prefix = "usgs_dem.";
   
   out << std::setiosflags(std::ios_base::fixed|std::ios_base::left)
       << std::setprecision(10)
      
       << prefix << std::setw(W)
       << "quadrangle_name:" << getQuadName().c_str() << "\n"
       << prefix << std::setw(W)
       << "process_info:" << getProcessInfo().c_str() << "\n"
       << prefix << std::setw(W)
       << "se_geo_corner_x:" << getSEGeoCornerX() << "\n"
       << prefix << std::setw(W)
       << "se_geo_corner_y:" << getSEGeoCornerY() << "\n"
       << prefix << std::setw(W)
       << "process_code:";
   
   ossim_int32 tmpl = getProcessCode() - 1;
   if ( tmpl >= 0 && tmpl <= MAX_PROCESS_CODE_INDEX)
   {
      out << PROCESS_CODE[tmpl] << "\n";
   }
   else
   {
      out << "unknown" << "\n";
   }

   out << prefix << std::setw(W)
       << "section_indicator: " << getSectionIndicator().c_str() << "\n"
       << prefix << std::setw(W)
       << "mapping_center_code: " << getMappingCenterCode().c_str() << "\n"
       << prefix << std::setw(W)
       << "level_code: " << getLevelCode() << "\n"
       << prefix << std::setw(W)
       << "elev_pattern: " << getElevPattern() << "\n"
       << prefix << std::setw(W) << "ground_ref_sys: ";
   
   tmpl = getGroundRefSysCode();
   if ( tmpl >= 0 && tmpl <= MAX_GROUND_REF_SYSTEM_INDEX)
   {
      out << GROUND_REF_SYSTEM[tmpl] << "\n";
   }
   else
   {
      out << "unknown\n";
   }
   
   out << prefix << std::setw(W)
       << "ground_ref_sys_zone: " << getGroundRefSysZone() << "\n"
       << prefix << std::setw(W)
       << "ground_ref_sys_units:";
   tmpl = getGroundRefSysUnits();
   if ( tmpl >= 0 && tmpl <= MAX_GROUND_REF_SYSTEM_UNITS_INDEX)
   {
      out << GROUND_REF_SYSTEM_UNITS[tmpl] << "\n";
   }
   else
   {
      out << "unknown\n";
   }
         
   out << prefix << std::setw(W)
       << "elevation_units: ";
   tmpl = getElevationUnits();
   if ( tmpl >= 0 && tmpl <= MAX_GROUND_REF_SYSTEM_UNITS_INDEX)
   {
      out << GROUND_REF_SYSTEM_UNITS[tmpl] << "\n";
   }
   else
   {
      out << "unknown\n";
   }

   out << prefix << std::setw(W)
       << "number_poly_sides: " << getNumPolySides() << "\n"
       << prefix << std::setw(W)
       << "counterclock_angle: " << getCounterclockAngle()<< "\n"
       << prefix << std::setw(W)
       << "elev_accuracy_code: " << getElevAccuracyCode()<< "\n"
       << prefix << std::setw(W)
       << "minimum_elevation: " << getMinimumElev() << "\n"
       << prefix << std::setw(W)
       << "maximum_elevation: " << getMaximumElev() << "\n"
       << prefix << std::setw(W)
       << "spatial_res_x: " << getSpatialResX() << "\n"
       << prefix << std::setw(W)
       << "spatial_res_y: " << getSpatialResY() << "\n"
       << prefix << std::setw(W)
       << "spatial_res_z:" << getSpatialResZ() << "\n"
       << prefix << std::setw(W)
       << "profile_rows: " << getProfileRows() << "\n"
       << prefix << std::setw(W)
       << "profile_columns:" << getProfileColumns() << "\n"
       << prefix << std::setw(W)
       << "source_date:" << getSourceDate() << "\n"
       << prefix << std::setw(W)
       << "revision_date:" << getInspRevDate() << "\n"
       << prefix  << std::setw(W)
       << "vertical_datum:";
   
   tmpl = getVertDatum() - 1;
   if ( tmpl >= 0 && tmpl <= MAX_VERTICAL_DATUM_INDEX)
   {
      out << VERTICAL_DATUM[tmpl] << "\n";
   }
   else
   {
      out << "unknown\n";
   }

   out << prefix << std::setw(W)
       << "vertical_datum_shift:" << getVertDatumShift() << "\n"
       << prefix << std::setw(W)
       << "horizontal_datum:";
   tmpl = getHorizDatum() - 1;
   if ( tmpl >= 0 && tmpl <= MAX_HORIZONTAL_DATUM_INDEX)
   {
      out << HORIZONTAL_DATUM[tmpl] << "\n";
   }
   else
   {
      out << "unknown\n";
   }

   const ossimDemPointVector CORNERS = getDEMCorners();
   for (unsigned int i=0; i < CORNERS.size(); ++i)
   {
      out << prefix << "corner[" << i
          << std::setw(CW)<< "].x:  " << CORNERS[i].getX() << "\n"
          << prefix << "corner[" << i
          << std::setw(CW)<< "].y:  " << CORNERS[i].getY() << "\n"
          << std::endl;
   }

   out << std::endl;

   // Reset flags.
   out.setf(f);   

   return out;
}

Friends And Related Function Documentation

operator<<

OSSIM_DLL std::ostream& operator<< ( ossim::ostream &  s, const ossimDemHeader &  header  )

friend

operator>>

OSSIM_DLL ossim::istream& operator>> ( ossim::istream &  s, ossimDemHeader &  header  )

friend

Definition at line 734 of file ossimDemHeader.cpp.

{
    return header.open(s);
}

Member Data Documentation

_counterclockAngle

double ossimDemHeader::_counterclockAngle

private

Definition at line 358 of file ossimDemHeader.h.

Referenced by getCounterclockAngle(), and open().

_dataEdition

ossim_int32 ossimDemHeader::_dataEdition

private

Definition at line 378 of file ossimDemHeader.h.

Referenced by getDataEdition(), and open().

_demCorners

ossimDemPointVector ossimDemHeader::_demCorners

private

Definition at line 357 of file ossimDemHeader.h.

Referenced by getDEMCorners(), getImageGeometry(), and open().

_eastEdgeFlag

ossim_int32 ossimDemHeader::_eastEdgeFlag

private

Definition at line 382 of file ossimDemHeader.h.

Referenced by getEastEdgeFlag(), and open().

_elevAccuracyCode

ossim_int32 ossimDemHeader::_elevAccuracyCode

private

Definition at line 359 of file ossimDemHeader.h.

Referenced by getElevAccuracyCode(), and open().

_elevPattern

ossim_int32 ossimDemHeader::_elevPattern

private

Definition at line 351 of file ossimDemHeader.h.

Referenced by getElevPattern(), and open().

_elevUnits

ossim_int32 ossimDemHeader::_elevUnits

private

Definition at line 355 of file ossimDemHeader.h.

Referenced by getElevationUnits(), and open().

_groundRefSysCode

ossim_int32 ossimDemHeader::_groundRefSysCode

private

Definition at line 352 of file ossimDemHeader.h.

Referenced by getGroundRefSysCode(), and open().

_groundRefSysUnits

ossim_int32 ossimDemHeader::_groundRefSysUnits

private

Definition at line 354 of file ossimDemHeader.h.

Referenced by getGroundRefSysUnits(), and open().

_groundRefSysZone

ossim_int32 ossimDemHeader::_groundRefSysZone

private

Definition at line 353 of file ossimDemHeader.h.

Referenced by getGroundRefSysZone(), and open().

_horizDatum

ossim_int32 ossimDemHeader::_horizDatum

private

Definition at line 377 of file ossimDemHeader.h.

Referenced by getHorizDatum(), and open().

_inspFlag

std::string ossimDemHeader::_inspFlag

private

Definition at line 373 of file ossimDemHeader.h.

Referenced by getInspFlag(), and open().

_inspRevDate

ossim_int32 ossimDemHeader::_inspRevDate

private

Definition at line 372 of file ossimDemHeader.h.

Referenced by getInspRevDate(), and open().

_largeContInt

ossim_int32 ossimDemHeader::_largeContInt

private

Definition at line 367 of file ossimDemHeader.h.

Referenced by getLargeContInt(), and open().

_levelCode

ossim_int32 ossimDemHeader::_levelCode

private

Definition at line 350 of file ossimDemHeader.h.

Referenced by getLevelCode(), and open().

_mapCenterCode

std::string ossimDemHeader::_mapCenterCode

private

Definition at line 349 of file ossimDemHeader.h.

Referenced by getMappingCenterCode(), and open().

_maxElevation

double ossimDemHeader::_maxElevation

private

Definition at line 361 of file ossimDemHeader.h.

Referenced by getMaximumElev(), and open().

_maxSourceUnits

ossim_int32 ossimDemHeader::_maxSourceUnits

private

Definition at line 368 of file ossimDemHeader.h.

Referenced by getMaxSourceUnits(), and open().

_minElevation

double ossimDemHeader::_minElevation

private

Definition at line 360 of file ossimDemHeader.h.

Referenced by getMinimumElev(), and open().

_minSourceUnits

ossim_int32 ossimDemHeader::_minSourceUnits

private

Definition at line 370 of file ossimDemHeader.h.

Referenced by getMinSourceUnits(), and open().

_northEdgeFlag

ossim_int32 ossimDemHeader::_northEdgeFlag

private

Definition at line 381 of file ossimDemHeader.h.

Referenced by getNorthEdgeFlag(), and open().

_numPolySides

ossim_int32 ossimDemHeader::_numPolySides

private

Definition at line 356 of file ossimDemHeader.h.

Referenced by getNumPolySides(), and open().

_perctVoid

ossim_int32 ossimDemHeader::_perctVoid

private

Definition at line 379 of file ossimDemHeader.h.

Referenced by getPerctVoid(), and open().

_processCode

ossim_int32 ossimDemHeader::_processCode

private

Definition at line 347 of file ossimDemHeader.h.

Referenced by getProcessCode(), and open().

_processInfo

std::string ossimDemHeader::_processInfo

private

Definition at line 344 of file ossimDemHeader.h.

Referenced by getProcessInfo(), and open().

_profileColumns

ossim_int32 ossimDemHeader::_profileColumns

private

Definition at line 366 of file ossimDemHeader.h.

Referenced by getProfileColumns(), and open().

_profileRows

ossim_int32 ossimDemHeader::_profileRows

private

Definition at line 365 of file ossimDemHeader.h.

Referenced by getProfileRows(), and open().

_quadName

std::string ossimDemHeader::_quadName

private

Definition at line 343 of file ossimDemHeader.h.

Referenced by getQuadName(), and open().

_sectionIndicator

std::string ossimDemHeader::_sectionIndicator

private

Definition at line 348 of file ossimDemHeader.h.

Referenced by getSectionIndicator(), and open().

_seGeoCornerX

double ossimDemHeader::_seGeoCornerX

private

Definition at line 345 of file ossimDemHeader.h.

Referenced by getSEGeoCornerX(), and open().

_seGeoCornerY

double ossimDemHeader::_seGeoCornerY

private

Definition at line 346 of file ossimDemHeader.h.

Referenced by getSEGeoCornerY(), and open().

_smallContInt

ossim_int32 ossimDemHeader::_smallContInt

private

Definition at line 369 of file ossimDemHeader.h.

Referenced by getSmallContInt(), and open().

_sourceDate

ossim_int32 ossimDemHeader::_sourceDate

private

Definition at line 371 of file ossimDemHeader.h.

Referenced by getSourceDate(), and open().

_southEdgeFlag

ossim_int32 ossimDemHeader::_southEdgeFlag

private

Definition at line 383 of file ossimDemHeader.h.

Referenced by getSouthEdgeFlag(), and open().

_spatialResX

double ossimDemHeader::_spatialResX

private

Definition at line 362 of file ossimDemHeader.h.

Referenced by getSpatialResX(), and open().

_spatialResY

double ossimDemHeader::_spatialResY

private

Definition at line 363 of file ossimDemHeader.h.

Referenced by getSpatialResY(), and open().

_spatialResZ

double ossimDemHeader::_spatialResZ

private

Definition at line 364 of file ossimDemHeader.h.

Referenced by getSpatialResZ(), and open().

_suspectVoidFlg

ossim_int32 ossimDemHeader::_suspectVoidFlg

private

Definition at line 375 of file ossimDemHeader.h.

Referenced by getSuspectVoidFlag(), and open().

_valFlag

ossim_int32 ossimDemHeader::_valFlag

private

Definition at line 374 of file ossimDemHeader.h.

Referenced by getDataValFlag(), and open().

_vertDatum

ossim_int32 ossimDemHeader::_vertDatum

private

Definition at line 376 of file ossimDemHeader.h.

Referenced by getVertDatum(), and open().

_vertDatumShift

double ossimDemHeader::_vertDatumShift

private

Definition at line 384 of file ossimDemHeader.h.

Referenced by getVertDatumShift(), and open().

_westEdgeFlag

ossim_int32 ossimDemHeader::_westEdgeFlag

private

Definition at line 380 of file ossimDemHeader.h.

Referenced by getWestEdgeFlag(), and open().

---

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

• ossimDemHeader.h
• ossimDemHeader.cpp