ossimLambertConformalConicProjection.h

Go to the documentation of this file.

//*******************************************************************
//
// License:  See top level LICENSE.txt file.
//
// Author:  Garrett Potts
//
// Description:
//
// Calls Lamberts Conformal Conic projection code.  
//*******************************************************************
//  $Id: ossimLambertConformalConicProjection.h 17815 2010-08-03 13:23:14Z dburken $

#ifndef ossimLambertConformalConicProjection_HEADER
#define ossimLambertConformalConicProjection_HEADER

#include <iosfwd>
#include <ossim/projection/ossimMapProjection.h>

class OSSIMDLLEXPORT ossimLambertConformalConicProjection : public ossimMapProjection
{
public:
   ossimLambertConformalConicProjection(const ossimEllipsoid& ellipsoid = ossimEllipsoid(),
                                        const ossimGpt& origin = ossimGpt());

   ossimLambertConformalConicProjection(const ossimEllipsoid& ellipsoid,
                                        const ossimGpt& origin,
                                        double stdParallel1,
                                        double stdParallel2,
                                        double falseEasting,
                                        double falseNorthing);
   
   virtual ~ossimLambertConformalConicProjection();
   
   virtual ossimObject *dup()const;
   
   virtual ossimGpt inverse(const ossimDpt &eastingNorthing)const;
   virtual ossimDpt forward(const ossimGpt &latLon)const;
   virtual void update();
   
   void setStandardParallel1(double degree);

   void setStandardParallel2(double degree);

   void setStandardParallels(double parallel1Degree,
                             double prallel2Degree);

   void setFalseEasting(double falseEasting);

   void setFalseNorthing(double falseNorthing);

   void setFalseEastingNorthing(double falseEasting, double falseNorthing);

   void setParameters(double parallel1,
                      double parallel2,
                      double falseEasting,
                      double falseNorthing);

   virtual double getStandardParallel1() const;

   virtual double getStandardParallel2() const;
   
   void setDefaults();

  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 bool operator==(const ossimProjection& projection) const;

private:

   //_______________GEOTRANS____________________
   mutable double Lambert_a;               /* Semi-major axis of ellipsoid, in meters */
   mutable double Lambert_f;       /* Flattening of ellipsoid */
   mutable double es;         /* Eccentricity of ellipsoid */
   mutable double es2;         /* Eccentricity squared */
   mutable double es_OVER_2;        /* Eccentricity / 2.0 */
   mutable double n;                 /* Ratio of angle between meridians */
   mutable double F;                 /* Flattening of ellipsoid */
   mutable double rho0;              /* Height above ellipsoid */
   mutable double Lambert_aF;        /* Lambert_a * F */
   
/* Lambert_Conformal_Conic projection Parameters */
   mutable double Lambert_Std_Parallel_1; /* Lower std. parallel, in radians */
   mutable double Lambert_Std_Parallel_2; /* Upper std. parallel, in radians */
   mutable double Lambert_Origin_Lat;     /* Latitude of origin, in radians */
   mutable double Lambert_Origin_Long;                /* Longitude of origin, in radians */
   mutable double Lambert_False_Northing;             /* False northing, in meters */
   mutable double Lambert_False_Easting;              /* False easting, in meters */

  long Set_Lambert_Parameters(double a,
                              double f,
                              double Origin_Latitude,
                              double Central_Meridian,
                              double Std_Parallel_1,
                              double Std_Parallel_2,
                              double False_Easting,
                              double False_Northing);




  void Get_Lambert_Parameters(double *a,
                              double *f,
                              double *Origin_Latitude,
                              double *Central_Meridian,
                              double *Std_Parallel_1,
                              double *Std_Parallel_2,
                              double *False_Easting,
                              double *False_Northing)const;

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

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


TYPE_DATA
};

#endif