mgFlatEartTranslation.h

Go to the documentation of this file.

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

#include "geodCoord.h"

#ifdef DtUSE_UTILITIES_NAMESPACE
namespace DtUSE_UTILITIES_NAMESPACE
{
#endif

inline void DtGocentricToMgFlatEarth(const DtVector& in, DtVector& out, const DtFloat64 &refLat, DtFloat64 &refLon ) const 
{ 
   const double theMultigenFlatEarthRadius = 6366707.02; 
   DtGeodeticCoord geodeticPosition; 

   geodeticPosition.setGeocentric(in); 
   double lat = DtRad2Deg(geodeticPosition.lat()); 
   double lon = DtRad2Deg(geodeticPosition.lon()); 

   out.setY(((lat - DtRad2Deg(refLat)) * 
      (theMultigenFlatEarthRadius * M_PI)) / 180.0); 

   out.setX(((lon - DtRad2Deg(refLon)) * 
      (theMultigenFlatEarthRadius * M_PI * cos(refLat))) / 180.0); 

   out.setZ(geodeticPosition.alt()); 
} 

#ifdef DtUSE_UTILITIES_NAMESPACE
}
#endif