ossimCylEquAreaProjection.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
// 
// Author:  Garrett Potts
//
// Description:
//
// Calls Geotrans Cylinder Equal Area  projection code.  
//*******************************************************************
//  $Id: ossimCylEquAreaProjection.h 17815 2010-08-03 13:23:14Z dburken $

#ifndef ossimCylEquAreaProjection_HEADER
#define ossimCylEquAreaProjection_HEADER

#include <ossim/projection/ossimMapProjection.h>

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

   virtual ossimDpt forward(const ossimGpt &worldPoint)    const;
   virtual ossimGpt inverse(const ossimDpt &projectedPoint)const;
   virtual void update();
   
   void setFalseEasting(double falseEasting);

   void setFalseNorthing(double falseNorthing);

   void setFalseEastingNorthing(double falseEasting, double falseNorthing);

   double getFalseEasting()const{return  Cyeq_False_Easting;}
   double getFalseNorthing()const{return Cyeq_False_Northing;}

   void setDefaults();

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

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


/* Ellipsoid Parameters, default to WGS 84 */
   mutable double Cyeq_a;                     /* Semi-major axis of ellipsoid in meters */
   mutable double Cyeq_f;             /* Flattening of ellipsoid */
   mutable double es2;            /* Eccentricity (0.08181919084262188000) squared  */
   mutable double es;                  /* Sqrt(es2) */
   mutable double es4;             /* es2 * es2 */
   mutable double es6;             /* es4 * es2 */
   mutable double k0;
   mutable double Cyeq_a_k0;                  /* Cyeq_a * k0 */
   mutable double two_k0;                           /* 2.0 * k0 */
   mutable double c0;                 /* es2 / 3.0 + 31.0 * es4 / 180.0 + 517.0 * es6 / 5040.0 */
   mutable double c1;              /* 23.0 es4 / 360.0 + 251.0 * es6 / 3780.0 */
   mutable double c2;              /* 761.0 * es6 / 45360.0 */
   
/* Cylindrical Equal Area projection Parameters */
   mutable double Cyeq_Origin_Lat;                  /* Latitude of origin in radians     */
   mutable double Cyeq_Origin_Long;                 /* Longitude of origin in radians    */
   mutable double Cyeq_False_Northing;              /* False northing in meters          */
   mutable double Cyeq_False_Easting;               /* False easting in meters           */
   
/* Maximum variance for easting and northing values for WGS 84.
 */
   mutable double Cyeq_Max_Easting;
   mutable double Cyeq_Min_Easting;
   mutable double Cyeq_Delta_Northing;

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

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



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


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

TYPE_DATA
};

#endif