utmReferencePoint.h

Go to the documentation of this file.

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

#ifdef DtUSE_UTILITIES_NAMESPACE
namespace DtUSE_UTILITIES_NAMESPACE
{
#endif

enum DtPolarDirection { DtEast, DtWest, DtNorth, DtSouth };


class DtDegMinSec
{
public:
   double deg, min, sec;
   DtPolarDirection direction;
};

class DT_DLL_MATRIX DtUtmReferencePoint
{
public:

   DtUtmReferencePoint();

   DtUtmReferencePoint( const DtDegMinSec& latRef,
                        const DtDegMinSec& lonRef,
                        int standardOffsets=1, 
                        long zone = 0,
                        const DtMapDatum& mapDatum = DtCurrentMapDatum(),
                        double scaleFactor = 0.9996);

   virtual ~DtUtmReferencePoint();

   virtual void geod2utm( double alat, double alon, double* x, double* y );

   virtual void utm2geod( double x, double y, double* lat, double* lon );

   virtual const DtMapDatum& mapDatum() const;

   virtual void useMapDatum( const DtMapDatum& mapDatum );

   virtual const DtReferenceEllipsoid* ellipsoid() const;

   void init(const DtDegMinSec& latRef,
             const DtDegMinSec& lonRef,
             int standardOffsets=1, 
             long zone = 0,
             const DtMapDatum& mapDatum = DtCurrentMapDatum() );

   virtual void uninitialize();
   
   virtual bool isInitialized() const;

   virtual double scale() const;
   virtual void setScale(double scale);

   virtual double xOrigin() const;
   virtual void setXOrgin(double x_origin);

   virtual double yOrigin() const;
   virtual void setYOrgin(double y_origin);

   virtual int standardOffsets() const;
   virtual void setStandardOffsets(int so);

   virtual DtDegMinSec latRef() const;
   virtual void setLatRef(const DtDegMinSec &lat_ref);

   virtual DtDegMinSec longRef() const;
   virtual void setLongRef(const DtDegMinSec &lat_ref);

   virtual long zoneOrigin() const;
   virtual void setZoneOrigin(long zo);
   
   virtual double standardFalseEasting() const;
   virtual void setStandardFalseEasting(double sfe);

   virtual double standardFalseNorthing() const;
   virtual void setStandardFalseNorthing(double sfn);

   virtual double standardFalseNorthingCrossover() const;
   virtual void setStandardFalseNorthingCrossover(double sfnc);

protected:

   void initDatumConstants();
   void initUtm( const double lat[4], const double lon[4], long zone = 0 );
   void degree( const double a[3], double* b ) const ;





   virtual void geod2utm( double lat, double lon,
      long* zone, double* northing, double* easting );
   virtual void utm2geod( double northing, double easting,
      long zone, double* lat, double* lon );

protected:

   int    myInited;
   int    myStandardOffsets;
   double myXOrigin;
   double myYOrigin; 

   long        myZoneOrigin;
   DtDegMinSec myLatRef;
   DtDegMinSec myLonRef;

   DtReferenceEllipsoid myEllipsoid;

   DtFloat64 mySemiMajor;
   DtFloat64 mySemiMinor;
   DtFloat64 myE_sq;
   DtFloat64 myEp_sq;
   DtFloat64 myEf;

   DtFloat64 myScale;
   double myStandardFalseEasting;
   double myStandardFalseNorthing;
   double myStandardFalseNorthingCrossover;

protected:


   static void setLongitudeFudgeFactor(double val);
   
   static double theLongitudeFudgeFactor;
};

DT_DLL_MATRIX DtUtmReferencePoint* DtGlobalUtmReferencePoint();

int DtUtmGridZoneFromLatLon(double latitude, double longitude);

#ifdef DtUSE_UTILITIES_NAMESPACE
}
#endif