referenceEllipsoid.h

Go to the documentation of this file.

/*******************************************************************************
** Copyright (c) 1992-2025 MAK Technologies, Inc
** All rights reserved.
*******************************************************************************/


#ifndef GeodeticEllipsoid_H_
#define GeodeticEllipsoid_H_

#include <vlutil/vlNetTypes.h>
#include "vlVector.h"
#include <vlutil/vlMachineTypes.h>
#include "mapDatum.h"


#ifdef DtUSE_UTILITIES_NAMESPACE
namespace DtUSE_UTILITIES_NAMESPACE
{
#endif


class DT_DLL_MATRIX DtReferenceEllipsoid
{
public:
    DtReferenceEllipsoid( const DtMapDatum& datum );
    virtual ~DtReferenceEllipsoid();

    void useMapDatum( const DtMapDatum& datum );
    const DtMapDatum& mapDatum() const;

    void geocentricToGeodetic( const DtVector& gc, DtVector& gd ) const;
    void geodeticToGeocentric( const DtVector& gd, DtVector& gc ) const;

protected:

    void initialize_ellipsoid( double a, double f );
    void geo_to_rad( double deg, double min, double sec, double* radian ) const;
    void rad_to_geo( double radian, double* deg, double* min, double* sec ) const;



    void geodetic_to_geocentric( double lat, double lon, double height,
       double* x, double* y, double* z ) const;
    
    void geocentric_to_geodetic( double xp, double yp, double zp,
       double* lat, double* lon, double* height ) const;

public:

   static void setMaxIterationsForConvergence(unsigned int maxItr);

   static unsigned int maxIterationsForConvergence();

   static void setThrowExceptionOnConversionFrom0Geocentric(bool enable);
   
   static void setNotifyConversionFrom0Geocentric(bool enable);

   static void setNotifyOnceConversionFrom0Geocentric(bool enableNotifyOnce);
   
protected:

   static unsigned int theMaxIterationsForConvergence;

   static bool theThrowExceptionOnConversionFrom0Geocentric;
   static bool theNotifyOnConversionFrom0Geocentric;
   static bool theNotifyOnceOnConversionFrom0Geocentric;

   DtMapDatum myMapDatum;
   int myInited;

   double ellipse_a;           
   double ellipse_a_sq;        
   double ellipse_f;           
   double ellipse_b;           
   double ellipse_b_sq;        
   double ellipse_a_sq_over_b; 
   double ellipse_e;           
   double ellipse_c1;          
};

DT_DLL_MATRIX DtReferenceEllipsoid* DtGlobalEllipsoid();

#ifdef DtUSE_UTILITIES_NAMESPACE
}
#endif

#endif