ossimBonneProjection.h

Go to the documentation of this file.

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

#include <ossim/projection/ossimMapProjection.h>

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

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

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

   
   mutable double Bonn_a;                 /* Semi-major axis of ellipsoid in meters */
   mutable double Bonn_f;         /* Flattening of ellipsoid */
   mutable double es2;        /* Eccentricity (0.08181919084262188000) squared  */
   mutable double es4;         /* es2 * es2 */
   mutable double es6;         /* es4 * es2 */
   mutable double M1;               /* Bonn_M(Bonna,c0lat,c1s2lat,c2s4lat,c3s6lat) */
   mutable double m1;               /* Bonn_m(coslat,sinlat,es2) */
   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 */
   mutable double a0;             /* 3.0 * e1 / 2.0 - 27.0 * e3 / 32.0 */
   mutable double a1;          /* 21.0 * e2 / 16.0 - 55.0 * e4 / 32.0 */
   mutable double a2;          /* 151.0 * e3 / 96.0 */
   mutable double a3;          /* 1097.0 * e4 / 512.0 */
   
   
/* Bonne projection Parameters */
   mutable double  Bonn_Origin_Lat;           /* Latitude of origin in radians     */
   mutable double  Bonn_Origin_Long;                          /* Longitude of origin in radians    */
   mutable double  Bonn_False_Northing;                       /* False northing in meters          */
   mutable double  Bonn_False_Easting;                        /* False easting in meters           */
   mutable double  Sin_Bonn_Origin_Lat;           /* sin(Bonn_Origin_Lat)              */
   mutable double  Bonn_am1sin;                   /* Bonn_a * m1 / Sin_Bonn_Origin_Lat */
   
/* Maximum variance for easting and northing values for WGS 84.
 */
   mutable double  Bonn_Max_Easting;
   mutable double  Bonn_Min_Easting;
   mutable double  Bonn_Delta_Northing;
   
  long Set_Bonne_Parameters(double a,
                            double f,
                            double Origin_Latitude,
                            double Central_Meridian,
                            double False_Easting,
                            double False_Northing);


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


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


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

TYPE_DATA
};

#endif