ossimMillerProjection.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
// 
// Author:  Garrett Potts (gpotts@imagelinks.com)
//
// Description:
//
// Calls Geotrans Miller projection code.  
//*******************************************************************
//  $Id: ossimMillerProjection.h 17815 2010-08-03 13:23:14Z dburken $

#ifndef ossimMillerProjection_HEADER
#define ossimMillerProjection_HEADER

#include <ossim/projection/ossimMapProjection.h>

class OSSIMDLLEXPORT ossimMillerProjection : public ossimMapProjection
{
public:
   ossimMillerProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
                         const ossimGpt& origin = ossimGpt());
   ossimMillerProjection(const ossimEllipsoid& ellipsoid,
                         const ossimGpt& origin,
                         const double falseEasting,
                         const double falseNorthing);
   
   ~ossimMillerProjection(){}
   
   virtual ossimObject *dup()const{return new ossimMillerProjection(*this);}

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

   void setFalseEasting(double falseEasting);

   void setFalseNorthing(double falseNorthing);

   void setFalseEastingNorthing(double falseEasting, double falseNorthing);

   void setDefaults();

   void setCentralMeridian(double centralMeridian);

   double getFalseEasting()const{return Mill_False_Easting;}
   double getFalseNorthing()const{return Mill_False_Northing;}
   virtual bool saveState(ossimKeywordlist& kwl,
                          const char* prefix=0)const;

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

   //________________GEOTRANS STUFF___________________

   mutable double Mill_a;   /* Semi-major axis of ellipsoid in meters */
   mutable double Mill_f;   /* Flattening of ellipsoid */
   mutable double es2;      /* Eccentricity (0.08181919084262188000) squared         */
   mutable double es4;      /* es2 * es2    */
   mutable double es6;      /* es4 * es2   */
   mutable double Ra;       /* Spherical Radius */
   
/* Miller projection Parameters */
   mutable double Mill_Origin_Long;                  /* Longitude of origin in radians    */
   mutable double Mill_False_Easting;
   mutable double Mill_False_Northing;
   mutable double Mill_Delta_Northing;
   mutable double Mill_Max_Easting;
   mutable double Mill_Min_Easting;
   
  long Set_Miller_Parameters(double a,
                             double f,
                             double Central_Meridian,
                             double False_Easting,
                             double False_Northing);


  void Get_Miller_Parameters(double *a,
                             double *f,
                             double *Central_Meridian,
                             double *False_Easting,
                             double *False_Northing)const;


  long Convert_Geodetic_To_Miller (double Latitude,
                                   double Longitude,
                                   double *Easting,
                                   double *Northing)const;


  long Convert_Miller_To_Geodetic(double Easting,
                                  double Northing,
                                  double *Latitude,
                                  double *Longitude)const;

TYPE_DATA
};

#endif