geodeticCoord.h

Go to the documentation of this file.

/*******************************************************************************
** Copyright (c) 1992-2010 VT MAK
** All rights reserved.
*******************************************************************************/


#ifndef DtGeodeticCoord_H_
#define DtGeodeticCoord_H_

#include "vlVector.h"
#include "vlMath.h"
#include "referenceEllipsoid.h"
#include "geoidDatum.h"
#include <iosfwd>



#ifdef DtUSE_UTILITIES_NAMESPACE
namespace DtUSE_UTILITIES_NAMESPACE
{
#endif

class DT_DLL_MATRIX DtGeodeticCoord : public DtVector
{
public:

   DtGeodeticCoord( const DtReferenceEllipsoid* ellipse = DtGlobalEllipsoid() );

   DtGeodeticCoord( double lat, double lon, double alt,
      const DtReferenceEllipsoid* ellipse = DtGlobalEllipsoid() );

   virtual ~DtGeodeticCoord();

   DtGeodeticCoord( const DtGeodeticCoord& src );

   DtGeodeticCoord& operator=( const DtGeodeticCoord& src );

   virtual DtGeodeticCoord& setGeocentric( const DtVector& geocentric );

   virtual void getGeocentric( DtVector& geocentric ) const;

   virtual DtVector geocentric() const;

   virtual void printDataToStream(std::ostream &str) const;

   virtual void setLat( DtFloat64 lat );
   virtual DtFloat64 lat() const;

   virtual void setLon( DtFloat64 lon );
   virtual DtFloat64 lon() const;

   virtual void setAlt( DtFloat64 alt);
   virtual DtFloat64 alt() const;

   virtual const DtMapDatum& mapDatum() const;

   virtual void setGeoidDatumType( DtGeoidType geoid );

protected:
   const DtReferenceEllipsoid* myEllipsoid;
   DtGeoidDatum* myGeoidDatum;
}; 

inline DtFloat64 DtGeodeticCoord::lat() const 
{ 
   return (*this)[0];
}

inline DtFloat64 DtGeodeticCoord::lon() const 
{ 
   return (*this)[1];
}

inline DtFloat64 DtGeodeticCoord::alt() const 
{ 
   if ( !myGeoidDatum )
   {
      return (*this)[2];
   }
   else
   {
      //convert long and lat to degrees
      double latDegrees = DtRad2Deg((*this)[0]);
      double lonDegrees = DtRad2Deg((*this)[1]);
      //Find undulation and subtract it from the altitude.
      double und = myGeoidDatum->undulation(latDegrees, lonDegrees);
      return (*this)[2] - und;
   }
}

inline void DtGeodeticCoord::setLat( DtFloat64 lat ) 
{ 
   (*this)[0] = lat; 
}

inline void DtGeodeticCoord::setLon( DtFloat64 lon ) 
{ 
   (*this)[1] = lon; 
}

inline void DtGeodeticCoord::setAlt( DtFloat64 alt ) 
{ 
   if ( !myGeoidDatum )
   {
      (*this)[2] = alt;
   }
   else
   {
      //convert long and lat to degrees
      double latDegrees = DtRad2Deg((*this)[0]);
      double lonDegrees = DtRad2Deg((*this)[1]);
      //Find undulation and add it to the altitude value.
      double und = myGeoidDatum->undulation(latDegrees, lonDegrees);
      (*this)[2] = alt + und;
   }
}

DT_DLL_MATRIX std::ostream & operator << (std::ostream &str, const  DtGeodeticCoord &sr);

DT_DLL_MATRIX void DtGeocToGeod( const DtVector& in, DtGeodeticCoord& out );

DT_DLL_MATRIX void DtGeodToGeoc( const DtGeodeticCoord& in, DtVector& out );

#ifdef DtUSE_UTILITIES_NAMESPACE
}
#endif

#endif