ossimUtmProjection.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
//
// Author:  Garrett Potts
//
// Description:
//
// Calls Geotrans Utm projection code.  
//*******************************************************************
//  $Id: ossimUtmProjection.h 17815 2010-08-03 13:23:14Z dburken $
#ifndef ossimUtmProjection_HEADER
#define ossimUtmProjection_HEADER
#include <ossim/projection/ossimMapProjection.h>

class OSSIMDLLEXPORT ossimUtmProjection : public ossimMapProjection
{
public:
   ossimUtmProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
                      const ossimGpt& origin = ossimGpt());
   
   ossimUtmProjection(const ossimEllipsoid& ellipsoid,
                      const ossimGpt& origin,
                      ossim_int32 zone,
                      char hemisphere);
   ossimUtmProjection(ossim_int32 zone);
   ossimUtmProjection(const ossimUtmProjection& src);
   
   virtual ossimObject* dup()const;

   virtual ossimGpt inverse(const ossimDpt &eastingNorthing)const;
   virtual ossimDpt forward(const ossimGpt &latLon)const;
   virtual void update();

   virtual void setOrigin(const ossimGpt& origin);
   
   void setZone(const ossimGpt& ground);
   void setZone(ossim_int32 zone);
   void setHemisphere(const ossimGpt& ground);
   void setHemisphere(char hemisphere);
   static ossim_int32 computeZone(const ossimGpt& gpt);

   static double computeZoneMeridian(ossim_int32 zone);

   ossim_int32 getZone()const;

   char getHemisphere()const;

   virtual bool saveState(ossimKeywordlist& kwl,
                          const char* prefix=0)const;
   
   virtual bool loadState(const ossimKeywordlist& kwl,
                          const char* prefix=0);

   virtual std::ostream& print(std::ostream& out) const;   

   virtual double getFalseEasting() const;

   virtual double getFalseNorthing() const;

   virtual bool operator==(const ossimProjection& projection) const;

   virtual ossim_uint32 getPcsCode() const;
   
private:
   
   /*_____________GEOTRANS_______________*/
   
   ossim_int32 Set_Transverse_Mercator_Parameters(double a,      
                                                  double f,
                                                  double Origin_Latitude,
                                                  double Central_Meridian,
                                                  double False_Easting,
                                                  double False_Northing,
                                                  double Scale_Factor);
   void Get_Transverse_Mercator_Parameters(double *a,
                                           double *f,
                                           double *Origin_Latitude,
                                           double *Central_Meridian,
                                           double *False_Easting,
                                           double *False_Northing,
                                           double *Scale_Factor)const;
   
   
   ossim_int32 Convert_Geodetic_To_Transverse_Mercator (double Latitude,
                                                        double Longitude,
                                                        double *Easting,
                                                        double *Northing)const;
   
   
   ossim_int32 Convert_Transverse_Mercator_To_Geodetic(double Easting,
                                                       double Northing,
                                                       double *Latitude,
                                                       double *Longitude)const;
   
   double theTranMerc_a;  /* Semi-major axis of ellipsoid i meters */
   double theTranMerc_f;  /* Flattening of ellipsoid  */
   
   /* Eccentricity (0.08181919084262188000) squared */
   double theTranMerc_es;
   
   double theTranMerc_ebs;   /* Second Eccentricity squared */
   
   /* Transverse_Mercator projection Parameters */
   
   double theTranMerc_Origin_Lat;     /* Latitude of origin in radians */
   double theTranMerc_Origin_Long;    /* Longitude of origin in radians */
   double theTranMerc_False_Northing; /* False northing in meters */
   double theTranMerc_False_Easting;  /* False easting in meters */
   double theTranMerc_Scale_Factor;   /* Scale factor  */
   
   /* Isometeric to geodetic latitude parameters, default to WGS 84 */
   double theTranMerc_ap;
   double theTranMerc_bp;
   double theTranMerc_cp;
   double theTranMerc_dp;
   double theTranMerc_ep;
   
   /* Maximum variance for easting and northing values for WGS 84. */
   double theTranMerc_Delta_Easting;
   double theTranMerc_Delta_Northing;  
   
   ossim_int32 theZone;
   
   char   theHemisphere;
   
TYPE_DATA
};

#endif