ossimPolyconicProjection.h

Go to the documentation of this file.

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

#ifndef ossimPolyconicProjection_HEADER
#define ossimPolyconicProjection_HEADER

#include <ossim/projection/ossimMapProjection.h>

class OSSIMDLLEXPORT ossimPolyconicProjection : public ossimMapProjection
{
public:   
   ossimPolyconicProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
                            const ossimGpt& origin = ossimGpt());
   ossimPolyconicProjection(const ossimEllipsoid& ellipsoid,
                            const ossimGpt& origin,
                            const double falseEasting,
                            const double falseNorthing);
   ~ossimPolyconicProjection(){}
   virtual ossimObject *dup()const
      {
         return new ossimPolyconicProjection(*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();

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

   virtual bool loadState(const ossimKeywordlist& kwl,
                          const char* prefix=0);
private:
   
   mutable double Poly_a;                      /* Semi-major axis of ellipsoid in meters */
   mutable double Poly_f;              /* Flattening of ellipsoid */
   mutable double es2;             /* Eccentricity (0.08181919084262188000) squared         */
   mutable double es4;              /* es2 * es2 */
   mutable double es6;              /* es4 * es2 */
   mutable double M0;
   mutable double c0;      /* 1 - es2 / 4.0 - 3.0 * es4 / 64.0 - 5.0 * es6 / 256.0 */
   mutable double c1; /* 3.0 * es2 / 8.0 + 3.0 * es4 / 32.0 + 45.0 * es6 / 1024.0 */
   mutable double c2;             /* 15.0 * es4 / 256.0 + 45.0 * es6 / 1024.0 */
   mutable double c3;               /* 35.0 * es6 / 3072.0 */
   
/* Polyconic projection Parameters */
   mutable double Poly_Origin_Lat;                   /* Latitude of origin in radians     */
   mutable double Poly_Origin_Long;                  /* Longitude of origin in radians    */
   mutable double Poly_False_Northing;               /* False northing in meters          */
   mutable double Poly_False_Easting;                /* False easting in meters           */
   
/* Maximum variance for easting and northing values for WGS 84.
 */
   mutable double Poly_Max_Easting;
   mutable double Poly_Max_Northing;
   mutable double Poly_Min_Easting;
   mutable double Poly_Min_Northing;

  long Set_Polyconic_Parameters(double a,
                                double f,
                                double Origin_Latitude,
                                double Central_Meridian,
                                double False_Easting,
                                double False_Northing);



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



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


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

TYPE_DATA
};

#endif