#include <ossimMercatorProjection.h>

Inheritance diagram for ossimMercatorProjection:

Public Member Functions
	ossimMercatorProjection (const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())

	ossimMercatorProjection (const ossimEllipsoid &ellipsoid, const ossimGpt &origin, double falseEasting, double falseNorthing, double scaleFactor)

	~ossimMercatorProjection ()

virtual ossimObject *	dup () const

virtual ossimGpt	inverse (const ossimDpt &eastingNorthing) const
	Will take a point in meters and convert it to ground. More...

virtual ossimDpt	forward (const ossimGpt &latLon) const
	All map projections will convert the world coordinate to an easting northing (Meters). More...

virtual void	update ()

void	setFalseEasting (double falseEasting)

void	setFalseNorthing (double falseNorthing)

void	setFalseEastingNorthing (double falseEasting, double falseNorthing)

void	setScaleFactor (double scaleFactor)

void	setParameters (double falseEasting, double falseNorthing, double scaleFactor)

void	setDefaults ()

double	getFalseEasting () const

double	getFalseNorthing () const

double	getScaleFactor () const

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
	Returns TRUE if principal parameters are within epsilon tolerance. More...

Public Member Functions inherited from ossimMapProjection
	ossimMapProjection (const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())

	ossimMapProjection (const ossimMapProjection &src)

virtual ossimGpt	origin () const

virtual ossimDpt	worldToLineSample (const ossimGpt &worldPoint) const

virtual void	worldToLineSample (const ossimGpt &worldPoint, ossimDpt &lineSample) const

virtual ossimGpt	lineSampleToWorld (const ossimDpt &projectedPoint) const

virtual void	lineSampleToWorld (const ossimDpt &projectedPoint, ossimGpt &gpt) const

virtual void	lineSampleHeightToWorld (const ossimDpt &lineSampPt, const double &heightAboveEllipsoid, ossimGpt &worldPt) const
	This is the pure virtual that projects the image point to the given elevation above ellipsoid, thereby bypassing reference to a DEM. More...

virtual void	lineSampleToEastingNorthing (const ossimDpt &liineSample, ossimDpt &eastingNorthing) const

virtual void	eastingNorthingToLineSample (const ossimDpt &eastingNorthing, ossimDpt &lineSample) const

virtual void	eastingNorthingToWorld (const ossimDpt &eastingNorthing, ossimGpt &worldPt) const

virtual double	getStandardParallel1 () const
	Derived classes should implement as needed. More...

virtual double	getStandardParallel2 () const
	Derived classes should implement as needed. More...

virtual void	setPcsCode (ossim_uint32 pcsCode)

virtual ossim_uint32	getPcsCode () const
	Returns the EPSG PCS code or 32767 if the projection is a custom (non-EPSG) projection. More...

virtual ossimString	getProjectionName () const
	Returns the projection name. More...

virtual double	getA () const
	ACCESS METHODS: More...

virtual double	getB () const

virtual double	getF () const

virtual ossimDpt	getMetersPerPixel () const

virtual const ossimDpt &	getDecimalDegreesPerPixel () const
	Returns decimal degrees per pixel as an ossimDpt with "x" representing longitude and "y" representing latitude. More...

virtual const ossimDpt &	getUlEastingNorthing () const

virtual const ossimGpt &	getUlGpt () const

virtual const ossimDatum *	getDatum () const

const ossimEllipsoid &	getEllipsoid () const

const ossimGpt &	getOrigin () const

virtual bool	isGeographic () const

virtual void	applyScale (const ossimDpt &scale, bool recenterTiePoint)
	Applies scale to theDeltaLonPerPixel, theDeltaLatPerPixel and theMetersPerPixel data members (eg: theDeltaLonPerPixel *= scale.x). More...

virtual void	setEllipsoid (const ossimEllipsoid &ellipsoid)
	SET METHODS: More...

virtual void	setAB (double a, double b)

virtual void	setDatum (const ossimDatum *datum)
	Sets theDatum to datum. More...

virtual void	setOrigin (const ossimGpt &origin)
	Sets theOrigin to origin. More...

virtual void	setMetersPerPixel (const ossimDpt &gsd)

virtual void	setDecimalDegreesPerPixel (const ossimDpt &gsd)

virtual void	setUlTiePoints (const ossimGpt &gpt)

virtual void	setUlTiePoints (const ossimDpt &eastingNorthing)

virtual void	setUlEastingNorthing (const ossimDpt &ulEastingNorthing)

virtual void	setUlGpt (const ossimGpt &ulGpt)

virtual void	assign (const ossimProjection &aProjection)

virtual std::ostream &	print (std::ostream &out) const
	Prints data members to stream. More...

virtual void	computeDegreesPerPixel ()
	Computes the approximate resolution in degrees/pixel. More...

virtual void	computeMetersPerPixel ()
	This will go from the ground point and give you an approximate meters per pixel. More...

void	setMatrix (double rotation, const ossimDpt &scale, const ossimDpt &translation)

void	setMatrixScale (const ossimDpt &scale)

void	setMatrixRotation (double rotation)

void	setMatrixTranslation (const ossimDpt &translation)

void	snapTiePointTo (ossim_float64 multiple, ossimUnitType unitType)
	Utility method to snap the tie point to some multiple. More...

void	snapTiePointToOrigin ()

void	setElevationLookupFlag (bool flag)

bool	getElevationLookupFlag () const

ossimUnitType	getModelTransformUnitType () const

void	setModelTransformUnitType (ossimUnitType unit)

bool	hasModelTransform () const

virtual bool	isAffectedByElevation () const
	Implementation of pure virtual ossimProjection::isAffectedByElevation method. More...

void	setProjectionUnits (ossimUnitType units)

ossimUnitType	getProjectionUnits () const
	OSSIM considers all map projection coordinates (including false eastings and northings) to be in meters. More...

virtual bool	isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const

Public Member Functions inherited from ossimProjection
	ossimProjection ()

virtual	~ossimProjection ()

virtual void	getRoundTripError (const ossimDpt &imagePoint, ossimDpt &errorResult) const

virtual void	getRoundTripError (const ossimGpt &groundPoint, ossimDpt &errorResult) const

virtual void	getGroundClipPoints (ossimGeoPolygon &gpts) const

virtual bool	operator!= (const ossimProjection &projection) const

Public Member Functions inherited from ossimObject
	ossimObject ()

virtual	~ossimObject ()

virtual ossimString	getShortName () const

virtual ossimString	getLongName () const

virtual ossimString	getDescription () const

virtual ossimString	getClassName () const

virtual RTTItypeid	getType () const

virtual bool	canCastTo (ossimObject *obj) const

virtual bool	canCastTo (const RTTItypeid &id) const

virtual bool	canCastTo (const ossimString &parentClassName) const

virtual void	accept (ossimVisitor &visitor)

Public Member Functions inherited from ossimReferenced
	ossimReferenced ()

	ossimReferenced (const ossimReferenced &)

ossimReferenced &	operator= (const ossimReferenced &)

void	ref () const
	increment the reference count by one, indicating that this object has another pointer which is referencing it. More...

void	unref () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

void	unref_nodelete () const
	decrement the reference count by one, indicating that a pointer to this object is referencing it. More...

int	referenceCount () const

Public Member Functions inherited from ossimErrorStatusInterface
	ossimErrorStatusInterface ()

virtual	~ossimErrorStatusInterface ()

virtual ossimErrorCode	getErrorStatus () const

virtual ossimString	getErrorStatusString () const

virtual void	setErrorStatus (ossimErrorCode error_status) const

virtual void	setErrorStatus () const

virtual void	clearErrorStatus () const

bool	hasError () const

Private Member Functions
long	Set_Mercator_Parameters (double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double *Scale_Factor)

void	Get_Mercator_Parameters (double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double *Scale_Factor) const

long	Convert_Geodetic_To_Mercator (double Latitude, double Longitude, double Easting, double Northing) const

long	Convert_Mercator_To_Geodetic (double Easting, double Northing, double Latitude, double Longitude) const

Private Attributes
double	Merc_a

double	Merc_f

double	Merc_e

double	Merc_es

double	Merc_Origin_Lat

double	Merc_Origin_Long

double	Merc_False_Northing

double	Merc_False_Easting

double	Merc_Scale_Factor

double	Merc_ab

double	Merc_bb

double	Merc_cb

double	Merc_db

double	Merc_Delta_Easting

double	Merc_Delta_Northing

bool	theSphericalFlag

Additional Inherited Members
Protected Member Functions inherited from ossimMapProjection
virtual	~ossimMapProjection ()

void	updateFromTransform ()

Protected Member Functions inherited from ossimReferenced
virtual	~ossimReferenced ()

Protected Attributes inherited from ossimMapProjection
ossimEllipsoid	theEllipsoid
	This method verifies that the projection parameters match the current pcs code. More...

ossimGpt	theOrigin

const ossimDatum *	theDatum
	This is only set if we want to have built in datum shifting. More...

ossimDpt	theMetersPerPixel
	Holds the number of meters per pixel. More...

ossimDpt	theDegreesPerPixel
	Hold the decimal degrees per pixel. More...

ossimGpt	theUlGpt
	Hold tie point in decimal degrees. More...

ossimDpt	theUlEastingNorthing
	Hold tie point as easting northing. More...

ossimDpt	theFalseEastingNorthing
	Hold the false easting northing. More...

ossim_uint32	thePcsCode
	Projection Coordinate System(PCS) code. More...

bool	theElevationLookupFlag

ossimMatrix4x4	theModelTransform

ossimMatrix4x4	theInverseModelTransform

ossimUnitType	theModelTransformUnitType

ossimUnitType	theProjectionUnits
	Linear units of the projection as indicated in the projection's specification: More...

Protected Attributes inherited from ossimErrorStatusInterface
ossimErrorCode	theErrorStatus

Detailed Description

Definition at line 18 of file ossimMercatorProjection.h.

Constructor & Destructor Documentation

◆ ossimMercatorProjection() [1/2]

ossimMercatorProjection::ossimMercatorProjection	(	const ossimEllipsoid &	ellipsoid = `ossimEllipsoid()`,
		const ossimGpt &	origin = `ossimGpt()`
	)

Definition at line 40 of file ossimMercatorProjection.cpp.

References setDefaults(), and update().

    :ossimMapProjection(ellipsoid, origin)
 {
    setDefaults();
    update();
 }

◆ ossimMercatorProjection() [2/2]

ossimMercatorProjection::ossimMercatorProjection	(	const ossimEllipsoid &	ellipsoid,
		const ossimGpt &	origin,
		double	falseEasting,
		double	falseNorthing,
		double	scaleFactor
	)

Definition at line 48 of file ossimMercatorProjection.cpp.

References Merc_False_Easting, Merc_False_Northing, Merc_Scale_Factor, setDefaults(), and update().

    :ossimMapProjection(ellipsoid, origin)
 {
    setDefaults();
    Merc_False_Easting  = falseEasting;
    Merc_False_Northing = falseNorthing;
    Merc_Scale_Factor   = scaleFactor;
 
    update();
 }

◆ ~ossimMercatorProjection()

ossimMercatorProjection::~ossimMercatorProjection ( )

inline

Definition at line 30 of file ossimMercatorProjection.h.

30 {}

Member Function Documentation

◆ Convert_Geodetic_To_Mercator()

long ossimMercatorProjection::Convert_Geodetic_To_Mercator	(	double	Latitude,
		double	Longitude,
		double *	Easting,
		double *	Northing
	)		const

private

The function Convert_Geodetic_To_Mercator converts geodetic (latitude and longitude) coordinates to Mercator projection easting, and northing coordinates, according to the current ellipsoid and Mercator projection parameters. If any errors occur, the error code(s) are returned by the function, otherwise MERC_NO_ERROR is returned.

Latitude : Latitude (phi) in radians (input) Longitude : Longitude (lambda) in radians (input) Easting : Easting (X) in meters (output) Northing : Northing (Y) in meters (output)

Definition at line 373 of file ossimMercatorProjection.cpp.

References MAX_LAT, Merc_a, Merc_e, Merc_False_Easting, Merc_False_Northing, MERC_LAT_ERROR, MERC_NO_ERROR, Merc_Origin_Long, Merc_Scale_Factor, and PI.

Referenced by forward(), and Set_Mercator_Parameters().

 { /* BEGIN Convert_Geodetic_To_Mercator */
    long Error_Code = MERC_NO_ERROR;
 /*
  * The function Convert_Geodetic_To_Mercator converts geodetic (latitude and
  * longitude) coordinates to Mercator projection (easting and northing)
  * coordinates, according to the current ellipsoid and Mercator projection
  * parameters.  If any errors occur, the error code(s) are returned by the
  * function, otherwise MERC_NO_ERROR is returned.
  *
  *    Latitude          : Latitude (phi) in radians           (input)
  *    Longitude         : Longitude (lambda) in radians       (input)
  *    Easting           : Easting (X) in meters               (output)
  *    Northing          : Northing (Y) in meters              (output)
  */
 
   double ctanz2;        /* Cotangent of z/2 - z - Isometric colatitude     */
   double e_x_sinlat;    /* e * sin(Latitude)                               */
   double Delta_Long;    /* Difference in origin longitude and longitude    */
   double tan_temp;
   double pow_temp;
 
 
    if ((Latitude < -MAX_LAT) || (Latitude > MAX_LAT))
    { /* Latitude out of range */
      Error_Code |= MERC_LAT_ERROR;
    }
 //   if ((Longitude < -PI) || (Longitude > (2*PI)))
 //   { /* Longitude out of range */
 //     Error_Code |= MERC_LON_ERROR;
 //   }
   if (!Error_Code)
   { /* no errors */
     if (Longitude > PI)
       Longitude -= (2*PI);
     e_x_sinlat = Merc_e * sin(Latitude);
     tan_temp = tan(PI / 4.e0 + Latitude / 2.e0);
     pow_temp = pow( ((1.e0 - e_x_sinlat) / (1.e0 + e_x_sinlat)),
                     (Merc_e / 2.e0) );
     ctanz2 = tan_temp * pow_temp;
     *Northing = Merc_Scale_Factor * Merc_a * log(ctanz2) + Merc_False_Northing;
     Delta_Long = Longitude - Merc_Origin_Long;
 //     if (Delta_Long > PI)
 //       Delta_Long -= (2 * PI);
 //     if (Delta_Long < -PI)
 //       Delta_Long += (2 * PI);
     *Easting = Merc_Scale_Factor * Merc_a * Delta_Long
                + Merc_False_Easting;
   }
 
   return (Error_Code);
 } /* END OF Convert_Geodetic_To_Mercator */

◆ Convert_Mercator_To_Geodetic()

long ossimMercatorProjection::Convert_Mercator_To_Geodetic	(	double	Easting,
		double	Northing,
		double *	Latitude,
		double *	Longitude
	)		const

private

The function Convert_Mercator_To_Geodetic converts Mercator projection easting and northing coordinates to geodetic (latitude and longitude) coordinates, according to the current ellipsoid and Mercator projection coordinates. If any errors occur, the error code(s) are returned by the function, otherwise MERC_NO_ERROR is returned.

Easting : Easting (X) in meters (input) Northing : Northing (Y) in meters (input) Latitude : Latitude (phi) in radians (output) Longitude : Longitude (lambda) in radians (output)

Definition at line 430 of file ossimMercatorProjection.cpp.

References M_PI, Merc_a, Merc_ab, Merc_bb, Merc_cb, Merc_db, Merc_False_Easting, Merc_False_Northing, MERC_NO_ERROR, Merc_Origin_Long, Merc_Scale_Factor, PI, and theSphericalFlag.

Referenced by inverse().

 { /* BEGIN Convert_Mercator_To_Geodetic */
 /*
  * The function Convert_Mercator_To_Geodetic converts Mercator projection
  * (easting and northing) coordinates to geodetic (latitude and longitude)
  * coordinates, according to the current ellipsoid and Mercator projection
  * coordinates.  If any errors occur, the error code(s) are returned by the
  * function, otherwise MERC_NO_ERROR is returned.
  *
  *    Easting           : Easting (X) in meters                  (input)
  *    Northing          : Northing (Y) in meters                 (input)
  *    Latitude          : Latitude (phi) in radians              (output)
  *    Longitude         : Longitude (lambda) in radians          (output)
  */
 
   double dx;     /* Delta easting - Difference in easting (easting-FE)      */
   double dy;     /* Delta northing - Difference in northing (northing-FN)   */
   double xphi;   /* Isometric latitude                                      */
    long Error_Code = MERC_NO_ERROR;
 
 #if 0
    if(theSphericalFlag)
    {
       *Latitude = M_PI*.5 - 2.0 * atan(exp(-Easting / Merc_a));
       *Longitude = Easting/Merc_a;
       
       return Error_Code;
    }
 #endif
 //   if ((Easting < (Merc_False_Easting - Merc_Delta_Easting))
 //       || (Easting > (Merc_False_Easting + Merc_Delta_Easting)))
 //   { /* Easting out of range */
 //     Error_Code |= MERC_EASTING_ERROR;
 //   }
 //   if ((Northing < (Merc_False_Northing - Merc_Delta_Northing))
 //       || (Northing > (Merc_False_Northing + Merc_Delta_Northing)))
 //   { /* Northing out of range */
 //     Error_Code |= MERC_NORTHING_ERROR;
 //   }
   if (!Error_Code)
   { /* no errors */
     dy = Northing - Merc_False_Northing;
     dx = Easting - Merc_False_Easting;
     *Longitude = Merc_Origin_Long + dx / (Merc_Scale_Factor * Merc_a);
     xphi = PI / 2.e0 
            - 2.e0 * atan(1.e0 / exp(dy / (Merc_Scale_Factor * Merc_a)));
     *Latitude = xphi + Merc_ab * sin(2.e0 * xphi) + Merc_bb * sin(4.e0 * xphi)
                 + Merc_cb * sin(6.e0 * xphi) + Merc_db * sin(8.e0 * xphi);
 //     if (*Longitude > PI)
 //       *Longitude -= (2 * PI);
 //     if (*Longitude < -PI)
 //       *Longitude += (2 * PI);
   }
   return (Error_Code);
 } /* END OF Convert_Mercator_To_Geodetic */

◆ dup()

virtual ossimObject* ossimMercatorProjection::dup ( ) const

inlinevirtual

Implements ossimProjection.

Definition at line 31 of file ossimMercatorProjection.h.

31 {return new ossimMercatorProjection(*this);}

ossimMercatorProjection::ossimMercatorProjection

ossimMercatorProjection(const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())

Definition: ossimMercatorProjection.cpp:40

◆ forward()

ossimDpt ossimMercatorProjection::forward ( const ossimGpt & worldPoint ) const

virtual

All map projections will convert the world coordinate to an easting northing (Meters).

Implements ossimMapProjection.

Definition at line 154 of file ossimMercatorProjection.cpp.

References ossimGpt::changeDatum(), ossimDatum::code(), Convert_Geodetic_To_Mercator(), ossimGpt::datum(), ossimGpt::latd(), ossimGpt::latr(), ossimGpt::lond(), ossimGpt::lonr(), M_PI, ossimDpt::makeNan(), Merc_a, MERC_NO_ERROR, ossimMapProjection::theDatum, and theSphericalFlag.

 {
    double easting  = 0.0;
    double northing = 0.0;
    ossimGpt gpt = latLon;
    ossimDpt result;
    if (theDatum)
    {
       if (theDatum->code() != latLon.datum()->code())
       {
          gpt.changeDatum(theDatum); // Shift to our datum.
       }
    }
    if(theSphericalFlag)
    {
       double lat = latLon.latd();
       double lon = latLon.lond();
       double shift = M_PI * Merc_a;
       easting = lon * shift / 180.0;
       northing = log( tan((90 + lat) * M_PI / 360.0 )) / (M_PI / 180.0);
       
       northing = northing * shift / 180.0;
       result = ossimDpt(easting, northing);
    }
    else 
    {
      long errorCode = Convert_Geodetic_To_Mercator(gpt.latr(),
                                    gpt.lonr(),
                                    &easting,
                                    &northing);
      if(errorCode!=MERC_NO_ERROR)
      {
       result.makeNan();
      }
      else
      {
       result = ossimDpt(easting, northing);
      }
    }      
 
    
    return result;
 }

◆ Get_Mercator_Parameters()

void ossimMercatorProjection::Get_Mercator_Parameters	(	double *	a,
		double *	f,
		double *	Origin_Latitude,
		double *	Central_Meridian,
		double *	False_Easting,
		double *	False_Northing,
		double *	Scale_Factor
	)		const

private

The function Get_Mercator_Parameters returns the current ellipsoid parameters, Mercator projection parameters, and scale factor.

a : Semi-major axis of ellipsoid, in meters (output) f : Flattening of ellipsoid (output) Origin_Latitude : Latitude in radians at which the (output) point scale factor is 1.0 Central_Meridian : Longitude in radians at the center of (output) the projection False_Easting : A coordinate value in meters assigned to the central meridian of the projection. (output) False_Northing : A coordinate value in meters assigned to the origin latitude of the projection (output) Scale_Factor : Multiplier which reduces distances in the projection to the actual distance on the ellipsoid (output)

Definition at line 334 of file ossimMercatorProjection.cpp.

References Merc_a, Merc_f, Merc_False_Easting, Merc_False_Northing, Merc_Origin_Lat, Merc_Origin_Long, and Merc_Scale_Factor.

 { /* BEGIN Get_Mercator_Parameters */
 /*
  * The function Get_Mercator_Parameters returns the current ellipsoid
  * parameters, Mercator projection parameters, and scale factor.
  *
  *    a                 : Semi-major axis of ellipsoid, in meters   (output)
  *    f                 : Flattening of ellipsoid                               (output)
  *    Origin_Latitude   : Latitude in radians at which the          (output)
  *                          point scale factor is 1.0
  *    Central_Meridian  : Longitude in radians at the center of     (output)
  *                          the projection
  *    False_Easting     : A coordinate value in meters assigned to the
  *                          central meridian of the projection.     (output)
  *    False_Northing    : A coordinate value in meters assigned to the
  *                          origin latitude of the projection       (output)
  *    Scale_Factor      : Multiplier which reduces distances in the 
  *                          projection to the actual distance on the
  *                          ellipsoid                               (output)
  */
 
   *a = Merc_a;
   *f = Merc_f;
   *Origin_Latitude = Merc_Origin_Lat;
   *Central_Meridian = Merc_Origin_Long;
   *False_Easting = Merc_False_Easting;
   *False_Northing = Merc_False_Northing;
   *Scale_Factor = Merc_Scale_Factor;
   
   return;
 } /* END OF Get_Mercator_Parameters */

◆ getFalseEasting()

double ossimMercatorProjection::getFalseEasting ( ) const

inlinevirtual

Returns: The false easting.

Reimplemented from ossimMapProjection.

Definition at line 70 of file ossimMercatorProjection.h.

70 {return Merc_False_Easting;}

ossimMercatorProjection::Merc_False_Easting

double Merc_False_Easting

Definition: ossimMercatorProjection.h:103

◆ getFalseNorthing()

double ossimMercatorProjection::getFalseNorthing ( ) const

inlinevirtual

Returns: The false northing.

Reimplemented from ossimMapProjection.

Definition at line 71 of file ossimMercatorProjection.h.

71 {return Merc_False_Northing;}

ossimMercatorProjection::Merc_False_Northing

double Merc_False_Northing

Definition: ossimMercatorProjection.h:102

◆ getScaleFactor()

double ossimMercatorProjection::getScaleFactor ( ) const

inline

Definition at line 72 of file ossimMercatorProjection.h.

72 {return Merc_Scale_Factor;}

ossimMercatorProjection::Merc_Scale_Factor

double Merc_Scale_Factor

Definition: ossimMercatorProjection.h:104

◆ inverse()

ossimGpt ossimMercatorProjection::inverse ( const ossimDpt & projectedPoint ) const

virtual

Will take a point in meters and convert it to ground.

Implements ossimMapProjection.

Definition at line 126 of file ossimMercatorProjection.cpp.

References Convert_Mercator_To_Geodetic(), M_PI, ossim::radiansToDegrees(), ossimMapProjection::theDatum, theSphericalFlag, ossimDpt::x, and ossimDpt::y.

 {
    double lat = 0.0;
    double lon = 0.0;
 
    if(theSphericalFlag)
    {
       double shift = M_PI * 6378137.0;
       lon = (eastingNorthing.x / shift) * 180.0;
       lat = (eastingNorthing.y / shift) * 180.0;
       
       lat = 180 / M_PI * (2 * atan( exp( lat * M_PI / 180.0)) - M_PI / 2.0);   
    }
    else 
    {
       Convert_Mercator_To_Geodetic(eastingNorthing.x,
                                    eastingNorthing.y,
                                    &lat,
                                    &lon);
       lat = ossim::radiansToDegrees(lat);
       lon = ossim::radiansToDegrees(lon);
    }
 
    
    
    return ossimGpt(lat, lon, 0.0, theDatum);  
 }

◆ loadState()

bool ossimMercatorProjection::loadState	(	const ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)

virtual

Method to the load (recreate) the state of an object from a keyword list. Return true if ok or false on error.

Reimplemented from ossimMapProjection.

Definition at line 209 of file ossimMercatorProjection.cpp.

References ossimKeywordlist::find(), ossimMapProjection::loadState(), Merc_False_Easting, Merc_False_Northing, Merc_Scale_Factor, ossimKeywordNames::SCALE_FACTOR_KW, setDefaults(), STATIC_TYPE_NAME, ossimMapProjection::theFalseEastingNorthing, ossimString::toDouble(), ossimKeywordNames::TYPE_KW, update(), ossimDpt::x, and ossimDpt::y.

 {
    bool flag = ossimMapProjection::loadState(kwl, prefix);
 
    const char* type          = kwl.find(prefix, ossimKeywordNames::TYPE_KW);
    const char* scaleFactor   = kwl.find(prefix, ossimKeywordNames::SCALE_FACTOR_KW);
    
    setDefaults();
 
    if(ossimString(type) == STATIC_TYPE_NAME(ossimMercatorProjection))
    {
       Merc_False_Easting  = theFalseEastingNorthing.x;
       Merc_False_Northing = theFalseEastingNorthing.y;
 
       if(scaleFactor)
       {
          Merc_Scale_Factor = ossimString(scaleFactor).toDouble();
       }
    }
    update();
    
    return flag;
 }

◆ operator==()

bool ossimMercatorProjection::operator== ( const ossimProjection & projection ) const

virtual

Returns TRUE if principal parameters are within epsilon tolerance.

Reimplemented from ossimMapProjection.

Definition at line 492 of file ossimMercatorProjection.cpp.

References ossim::almostEqual(), and Merc_Scale_Factor.

 {
    if (!ossimMapProjection::operator==(proj))
       return false;
 
    const ossimMercatorProjection* p = dynamic_cast<const ossimMercatorProjection*>(&proj);
    if (!p) return false;
    
    if (!ossim::almostEqual(Merc_Scale_Factor,p->Merc_Scale_Factor)) return false;
 
    return true;
 }

◆ saveState()

bool ossimMercatorProjection::saveState	(	ossimKeywordlist &	kwl,
		const char *	prefix = `0`
	)		const

virtual

Method to save the state of an object to a keyword list. Return true if ok or false on error.

Reimplemented from ossimMapProjection.

Definition at line 199 of file ossimMercatorProjection.cpp.

References ossimKeywordlist::add(), Merc_Scale_Factor, ossimMapProjection::saveState(), and ossimKeywordNames::SCALE_FACTOR_KW.

 {
    kwl.add(prefix,
            ossimKeywordNames::SCALE_FACTOR_KW,
            Merc_Scale_Factor,
            true);
 
    return ossimMapProjection::saveState(kwl, prefix);
 }

◆ Set_Mercator_Parameters()

long ossimMercatorProjection::Set_Mercator_Parameters	(	double	a,
		double	f,
		double	Origin_Latitude,
		double	Central_Meridian,
		double	False_Easting,
		double	False_Northing,
		double *	Scale_Factor
	)

private

The function Set_Mercator_Parameters receives the ellipsoid parameters and Mercator projcetion parameters as inputs, and sets the corresponding state variables. It calculates and returns the scale factor. If any errors occur, the error code(s) are returned by the function, otherwise MERC_NO_ERROR is returned.

a : Semi-major axis of ellipsoid, in meters (input) f : Flattening of ellipsoid (input) Origin_Latitude : Latitude in radians at which the (input) point scale factor is 1.0 Central_Meridian : Longitude in radians at the center of (input) the projection False_Easting : A coordinate value in meters assigned to the central meridian of the projection. (input) False_Northing : A coordinate value in meters assigned to the origin latitude of the projection (input) Scale_Factor : Multiplier which reduces distances in the projection to the actual distance on the ellipsoid (output)

Definition at line 241 of file ossimMercatorProjection.cpp.

References ossim::almostEqual(), Convert_Geodetic_To_Mercator(), MAX_LAT, Merc_a, Merc_ab, Merc_bb, Merc_cb, Merc_db, Merc_Delta_Easting, Merc_Delta_Northing, Merc_e, Merc_es, Merc_f, Merc_False_Easting, Merc_False_Northing, MERC_NO_ERROR, Merc_Origin_Lat, Merc_Origin_Long, Merc_Scale_Factor, PI, and theSphericalFlag.

Referenced by update().

 { /* BEGIN Set_Mercator_Parameters */
 /*
  * The function Set_Mercator_Parameters receives the ellipsoid parameters and
  * Mercator projection parameters as inputs, and sets the corresponding state 
  * variables.  It calculates and returns the scale factor.  If any errors
  * occur, the error code(s) are returned by the function, otherwise Merc_NO_ERROR
  * is returned.
  *
  *    a                 : Semi-major axis of ellipsoid, in meters   (input)
  *    f                 : Flattening of ellipsoid                               (input)
  *    Origin_Latitude   : Latitude in radians at which the          (input)
  *                          point scale factor is 1.0
  *    Central_Meridian  : Longitude in radians at the center of     (input)
  *                          the projection
  *    False_Easting     : A coordinate value in meters assigned to the
  *                          central meridian of the projection.     (input)
  *    False_Northing    : A coordinate value in meters assigned to the
  *                          origin latitude of the projection       (input)
  *    Scale_Factor      : Multiplier which reduces distances in the 
  *                          projection to the actual distance on the
  *                          ellipsoid                               (output)
  */
 
   double es2;   /* Eccentricity squared of ellipsoid to the second power    */
   double es3;   /* Eccentricity squared of ellipsoid to the third power     */
   double es4;   /* Eccentricity squared of ellipsoid to the fourth power    */
   double sin_olat; /* sin(Origin_Latitude), temp variable */
 //  double inv_f = 1 / f;
   long Error_Code = MERC_NO_ERROR;
 
    theSphericalFlag = ossim::almostEqual(f, 0.0);
 
 //   if (a <= 0.0)
 //   { /* Semi-major axis must be greater than zero */
 //     Error_Code |= MERC_A_ERROR;
 //   }
 //   if ((inv_f < 250) || (inv_f > 350))
 //   { /* Inverse flattening must be between 250 and 350 */
 //     Error_Code |= MERC_INV_F_ERROR;
 //   }
 //   if ((Origin_Latitude < -MAX_LAT) || (Origin_Latitude > MAX_LAT))
 //   { /* origin latitude out of range */
 //     Error_Code |= MERC_ORIGIN_LAT_ERROR;
 //   }
 //   if ((Central_Meridian < -PI) || (Central_Meridian > (2*PI)))
 //   { /* origin longitude out of range */
 //     Error_Code |= MERC_CENT_MER_ERROR;
 //   }
   if (!Error_Code)
   { /* no errors */
     Merc_a = a;
     Merc_f = f;
     Merc_Origin_Lat = Origin_Latitude;
 //     if (Central_Meridian > PI)
 //       Central_Meridian -= (2*PI);
     Merc_Origin_Long = Central_Meridian;
     Merc_False_Northing = False_Northing;
     Merc_False_Easting = False_Easting;
     Merc_es = 2 * Merc_f - Merc_f * Merc_f;
     Merc_e = sqrt(Merc_es);
     sin_olat = sin(Origin_Latitude);
     Merc_Scale_Factor = 1.0 / ( sqrt(1.e0 - Merc_es * sin_olat * sin_olat) 
                                 / cos(Origin_Latitude) );
     es2 = Merc_es * Merc_es;
     es3 = es2 * Merc_es;
     es4 = es3 * Merc_es;
     Merc_ab = Merc_es / 2.e0 + 5.e0 * es2 / 24.e0 + es3 / 12.e0
               + 13.e0 * es4 / 360.e0;
     Merc_bb = 7.e0 * es2 / 48.e0 + 29.e0 * es3 / 240.e0 
               + 811.e0 * es4 / 11520.e0;
     Merc_cb = 7.e0 * es3 / 120.e0 + 81.e0 * es4 / 1120.e0;
     Merc_db = 4279.e0 * es4 / 161280.e0;
     *Scale_Factor = Merc_Scale_Factor;
     Convert_Geodetic_To_Mercator( MAX_LAT, (Merc_Origin_Long + PI),                                           
                                   &Merc_Delta_Easting, &Merc_Delta_Northing);
     if (Merc_Delta_Easting < 0)
       Merc_Delta_Easting = -Merc_Delta_Easting;
     Merc_Delta_Easting *= 1.01;
     Merc_Delta_Easting -= Merc_False_Easting;
     Merc_Delta_Northing *= 1.01;
     Merc_Delta_Northing -= Merc_False_Northing;
   } /* END OF if(!Error_Code) */
   return (Error_Code);
 } /* END OF Set_Mercator_Parameters */

◆ setDefaults()

void ossimMercatorProjection::setDefaults ( )

Definition at line 117 of file ossimMercatorProjection.cpp.

References Merc_Delta_Easting, Merc_Delta_Northing, Merc_False_Easting, Merc_False_Northing, and Merc_Scale_Factor.

Referenced by loadState(), and ossimMercatorProjection().

 {
    Merc_False_Easting  = 0.0;
    Merc_False_Northing = 0.0;
    Merc_Delta_Easting  = 20237883.0;
    Merc_Delta_Northing = 23421740.0;
    Merc_Scale_Factor   = 1.0;
 }

◆ setFalseEasting()

void ossimMercatorProjection::setFalseEasting ( double falseEasting )

SetFalseEasting. The value is in meters. Update is then called so we can pre-compute paramters

Definition at line 80 of file ossimMercatorProjection.cpp.

References Merc_False_Easting, and update().

 {
    Merc_False_Easting = falseEasting;
    update();
 }

◆ setFalseEastingNorthing()

void ossimMercatorProjection::setFalseEastingNorthing	(	double	falseEasting,
		double	falseNorthing
	)

Sets both false easting and northing values. The values are expected to be in meters. Update is then called so we can pre-compute paramters

Definition at line 98 of file ossimMercatorProjection.cpp.

References Merc_False_Easting, Merc_False_Northing, and update().

 {
    Merc_False_Easting = falseEasting;
    Merc_False_Northing = falseNorthing;
    update();
 }

◆ setFalseNorthing()

void ossimMercatorProjection::setFalseNorthing ( double falseNorthing )

SetFalseNorthing. The value is in meters. Update is then called so we can pre-compute paramters

Definition at line 86 of file ossimMercatorProjection.cpp.

References Merc_False_Northing, and update().

 {
    Merc_False_Northing = falseNorthing;
    update();
 }

◆ setParameters()

void ossimMercatorProjection::setParameters	(	double	falseEasting,
		double	falseNorthing,
		double	scaleFactor
	)

Allows one to set all parameters for this projections. parallels are in degrees and easting northings are in meters. Update is then called so we can pre-compute paramters.

Definition at line 106 of file ossimMercatorProjection.cpp.

References Merc_False_Easting, Merc_False_Northing, Merc_Scale_Factor, and update().

 {
    Merc_False_Easting  = falseEasting;
    Merc_False_Northing = falseNorthing;
    Merc_Scale_Factor   = scaleFactor;
    
    update(); 
 }

◆ setScaleFactor()

void ossimMercatorProjection::setScaleFactor ( double scaleFactor )

sets the scale.

Definition at line 92 of file ossimMercatorProjection.cpp.

References Merc_Scale_Factor, and update().

Referenced by ossimWktProjectionFactory::doMercator().

 {
    Merc_Scale_Factor = scaleFactor;
    update();
 }

◆ update()

void ossimMercatorProjection::update ( )

virtual

Reimplemented from ossimMapProjection.

Definition at line 63 of file ossimMercatorProjection.cpp.

References ossimEllipsoid::getA(), ossimEllipsoid::getFlattening(), ossimGpt::latr(), ossimGpt::lonr(), Merc_False_Easting, Merc_False_Northing, Merc_Scale_Factor, Set_Mercator_Parameters(), ossimMapProjection::theEllipsoid, ossimMapProjection::theFalseEastingNorthing, ossimMapProjection::theOrigin, ossimMapProjection::update(), ossimDpt::x, and ossimDpt::y.

Referenced by loadState(), ossimMercatorProjection(), setFalseEasting(), setFalseEastingNorthing(), setFalseNorthing(), setParameters(), and setScaleFactor().

 {
    
    Set_Mercator_Parameters(theEllipsoid.getA(),
                            theEllipsoid.getFlattening(),
                            theOrigin.latr(),
                            theOrigin.lonr(),
                            Merc_False_Easting,
                            Merc_False_Northing,
                            &Merc_Scale_Factor);
 
    theFalseEastingNorthing.x = Merc_False_Easting;
    theFalseEastingNorthing.y = Merc_False_Northing;
 
    ossimMapProjection::update();
 }

Member Data Documentation

◆ Merc_a

double ossimMercatorProjection::Merc_a

mutableprivate

Definition at line 94 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), Convert_Mercator_To_Geodetic(), forward(), Get_Mercator_Parameters(), and Set_Mercator_Parameters().

◆ Merc_ab

double ossimMercatorProjection::Merc_ab

mutableprivate

Definition at line 107 of file ossimMercatorProjection.h.

Referenced by Convert_Mercator_To_Geodetic(), and Set_Mercator_Parameters().

◆ Merc_bb

double ossimMercatorProjection::Merc_bb

mutableprivate

Definition at line 108 of file ossimMercatorProjection.h.

Referenced by Convert_Mercator_To_Geodetic(), and Set_Mercator_Parameters().

◆ Merc_cb

double ossimMercatorProjection::Merc_cb

mutableprivate

Definition at line 109 of file ossimMercatorProjection.h.

Referenced by Convert_Mercator_To_Geodetic(), and Set_Mercator_Parameters().

◆ Merc_db

double ossimMercatorProjection::Merc_db

mutableprivate

Definition at line 110 of file ossimMercatorProjection.h.

Referenced by Convert_Mercator_To_Geodetic(), and Set_Mercator_Parameters().

◆ Merc_Delta_Easting

double ossimMercatorProjection::Merc_Delta_Easting

mutableprivate

Definition at line 114 of file ossimMercatorProjection.h.

Referenced by Set_Mercator_Parameters(), and setDefaults().

◆ Merc_Delta_Northing

double ossimMercatorProjection::Merc_Delta_Northing

mutableprivate

Definition at line 115 of file ossimMercatorProjection.h.

Referenced by Set_Mercator_Parameters(), and setDefaults().

◆ Merc_e

double ossimMercatorProjection::Merc_e

mutableprivate

Definition at line 96 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), and Set_Mercator_Parameters().

◆ Merc_es

double ossimMercatorProjection::Merc_es

mutableprivate

Definition at line 97 of file ossimMercatorProjection.h.

Referenced by Set_Mercator_Parameters().

◆ Merc_f

double ossimMercatorProjection::Merc_f

mutableprivate

Definition at line 95 of file ossimMercatorProjection.h.

Referenced by Get_Mercator_Parameters(), and Set_Mercator_Parameters().

◆ Merc_False_Easting

double ossimMercatorProjection::Merc_False_Easting

mutableprivate

Definition at line 103 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), Convert_Mercator_To_Geodetic(), Get_Mercator_Parameters(), loadState(), ossimMercatorProjection(), Set_Mercator_Parameters(), setDefaults(), setFalseEasting(), setFalseEastingNorthing(), setParameters(), and update().

◆ Merc_False_Northing

double ossimMercatorProjection::Merc_False_Northing

mutableprivate

Definition at line 102 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), Convert_Mercator_To_Geodetic(), Get_Mercator_Parameters(), loadState(), ossimMercatorProjection(), Set_Mercator_Parameters(), setDefaults(), setFalseEastingNorthing(), setFalseNorthing(), setParameters(), and update().

◆ Merc_Origin_Lat

double ossimMercatorProjection::Merc_Origin_Lat

mutableprivate

Definition at line 100 of file ossimMercatorProjection.h.

Referenced by Get_Mercator_Parameters(), and Set_Mercator_Parameters().

◆ Merc_Origin_Long

double ossimMercatorProjection::Merc_Origin_Long

mutableprivate

Definition at line 101 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), Convert_Mercator_To_Geodetic(), Get_Mercator_Parameters(), and Set_Mercator_Parameters().

◆ Merc_Scale_Factor

double ossimMercatorProjection::Merc_Scale_Factor

mutableprivate

Definition at line 104 of file ossimMercatorProjection.h.

Referenced by Convert_Geodetic_To_Mercator(), Convert_Mercator_To_Geodetic(), Get_Mercator_Parameters(), loadState(), operator==(), ossimMercatorProjection(), saveState(), Set_Mercator_Parameters(), setDefaults(), setParameters(), setScaleFactor(), and update().

◆ theSphericalFlag

bool ossimMercatorProjection::theSphericalFlag

private

Definition at line 118 of file ossimMercatorProjection.h.

Referenced by Convert_Mercator_To_Geodetic(), forward(), inverse(), and Set_Mercator_Parameters().

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimMercatorProjection() [1/2]

◆ ossimMercatorProjection() [2/2]

◆ ~ossimMercatorProjection()

Member Function Documentation

◆ Convert_Geodetic_To_Mercator()

◆ Convert_Mercator_To_Geodetic()

◆ dup()

◆ forward()

◆ Get_Mercator_Parameters()

◆ getFalseEasting()

◆ getFalseNorthing()

◆ getScaleFactor()

◆ inverse()

◆ loadState()

◆ operator==()

◆ saveState()

◆ Set_Mercator_Parameters()

◆ setDefaults()

◆ setFalseEasting()

◆ setFalseEastingNorthing()

◆ setFalseNorthing()

◆ setParameters()

◆ setScaleFactor()

◆ update()

Member Data Documentation

◆ Merc_a

◆ Merc_ab

◆ Merc_bb

◆ Merc_cb

◆ Merc_db

◆ Merc_Delta_Easting

◆ Merc_Delta_Northing

◆ Merc_e

◆ Merc_es

◆ Merc_f

◆ Merc_False_Easting

◆ Merc_False_Northing

◆ Merc_Origin_Lat

◆ Merc_Origin_Long

◆ Merc_Scale_Factor

◆ theSphericalFlag