#include <ossimUtmProjection.h>

Inheritance diagram for ossimUtmProjection:

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

	ossimUtmProjection (const ossimEllipsoid &ellipsoid, const ossimGpt &origin, ossim_int32 zone, char hemisphere)

	ossimUtmProjection (ossim_int32 zone)

	ossimUtmProjection (const ossimUtmProjection &src)

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 ()

virtual void	setOrigin (const ossimGpt &origin)
	This will set the utm zone and utm origin base on origin passed in. More...

void	setZone (const ossimGpt &ground)

void	setZone (ossim_int32 zone)

void	setHemisphere (const ossimGpt &ground)

void	setHemisphere (char hemisphere)

ossim_int32	getZone () const

char	getHemisphere () const

virtual bool	saveState (ossimKeywordlist &kwl, const char *prefix=0) const
	Method to save the state of an object to a keyword list. More...

virtual bool	loadState (const ossimKeywordlist &kwl, const char *prefix=0)
	Method to the load (recreate) the state of an object from a keyword list. More...

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

virtual double	getFalseEasting () const

virtual double	getFalseNorthing () const

virtual bool	operator== (const ossimProjection &projection) const
	Returns TRUE if principal parameters are within epsilon tolerance. More...

virtual ossim_uint32	getPcsCode () const
	Returns the EPSG PCS code based on datum, zone, and hemisphere. 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 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	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 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

Static Public Member Functions
static ossim_int32	computeZone (const ossimGpt &gpt)

static double	computeZoneMeridian (ossim_int32 zone)
	Return in decimal degrees the zone meridian. More...

Private Member Functions
ossim_int32	Set_Transverse_Mercator_Parameters (double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double Scale_Factor)
	The function Set_Tranverse_Mercator_Parameters receives the ellipsoid parameters and Tranverse Mercator projection parameters as inputs, and sets the corresponding state variables. More...

void	Get_Transverse_Mercator_Parameters (double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double *Scale_Factor) const
	The function Get_Transverse_Mercator_Parameters returns the current ellipsoid and Transverse Mercator projection parameters. More...

ossim_int32	Convert_Geodetic_To_Transverse_Mercator (double Latitude, double Longitude, double Easting, double Northing) const
	The function Convert_Geodetic_To_Transverse_Mercator converts geodetic (latitude and longitude) coordinates to Transverse Mercator projection (easting and northing) coordinates, according to the current ellipsoid and Transverse Mercator projection coordinates. More...

ossim_int32	Convert_Transverse_Mercator_To_Geodetic (double Easting, double Northing, double Latitude, double Longitude) const
	The function Convert_Transverse_Mercator_To_Geodetic converts Transverse Mercator projection (easting and northing) coordinates to geodetic (latitude and longitude) coordinates, according to the current ellipsoid and Transverse Mercator projection parameters. More...

Private Attributes
double	theTranMerc_a

double	theTranMerc_f

double	theTranMerc_es

double	theTranMerc_ebs

double	theTranMerc_Origin_Lat

double	theTranMerc_Origin_Long

double	theTranMerc_False_Northing

double	theTranMerc_False_Easting

double	theTranMerc_Scale_Factor

double	theTranMerc_ap

double	theTranMerc_bp

double	theTranMerc_cp

double	theTranMerc_dp

double	theTranMerc_ep

double	theTranMerc_Delta_Easting

double	theTranMerc_Delta_Northing

ossim_int32	theZone
	zone can be from 1 through 60 (0 == NOT SET) More...

char	theHemisphere
	can be N or S. More...

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 16 of file ossimUtmProjection.h.

Constructor & Destructor Documentation

◆ ossimUtmProjection() [1/4]

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

Definition at line 61 of file ossimUtmProjection.cpp.

References ossimMapProjection::origin(), setHemisphere(), setZone(), and update().

Referenced by dup().

    :
       ossimMapProjection(ellipsoid, origin),
       theTranMerc_a(6378137.0),
       theTranMerc_f(1.0/298.257223563),
       theTranMerc_es(0.0066943799901413800),
       theTranMerc_ebs(0.0067394967565869),
       theTranMerc_Origin_Lat(0.0),
       theTranMerc_Origin_Long(0.0),
       theTranMerc_False_Northing(0.0),
       theTranMerc_False_Easting(500000.0),
       theTranMerc_Scale_Factor(0.9996),
       theTranMerc_ap(6367449.1458008),
       theTranMerc_bp(16038.508696861),
       theTranMerc_cp(16.832613334334),
       theTranMerc_dp(0.021984404273757),
       theTranMerc_Delta_Easting(40000000.0),
       theTranMerc_Delta_Northing(40000000.0),
       theZone(1),
       theHemisphere('N')
       
 {
    setZone(origin);
    setHemisphere(origin);
    update();
 }

◆ ossimUtmProjection() [2/4]

ossimUtmProjection::ossimUtmProjection	(	const ossimEllipsoid &	ellipsoid,
		const ossimGpt &	origin,
		ossim_int32	zone,
		char	hemisphere
	)

Definition at line 89 of file ossimUtmProjection.cpp.

References setHemisphere(), setZone(), and update().

    :
       ossimMapProjection(ellipsoid, origin),
       theTranMerc_a(6378137.0),
       theTranMerc_f(1.0/298.257223563),
       theTranMerc_es(0.0066943799901413800),
       theTranMerc_ebs(0.0067394967565869),
       theTranMerc_Origin_Lat(0.0),
       theTranMerc_Origin_Long(0.0),
       theTranMerc_False_Northing(0.0),
       theTranMerc_False_Easting(500000.0),
       theTranMerc_Scale_Factor(0.9996),
       theTranMerc_ap(6367449.1458008),
       theTranMerc_bp(16038.508696861),
       theTranMerc_cp(16.832613334334),
       theTranMerc_dp(0.021984404273757),
       theTranMerc_Delta_Easting(40000000.0),
       theTranMerc_Delta_Northing(40000000.0),
       theZone(zone),
       theHemisphere(hemisphere)      
 {
    setZone(zone);
    setHemisphere(hemisphere);
    update();
 }

◆ ossimUtmProjection() [3/4]

ossimUtmProjection::ossimUtmProjection ( ossim_int32 zone )

Definition at line 118 of file ossimUtmProjection.cpp.

References setZone(), and update().

    :
       ossimMapProjection(),
       theTranMerc_a(6378137.0),
       theTranMerc_f(1.0/298.257223563),
       theTranMerc_es(0.0066943799901413800),
       theTranMerc_ebs(0.0067394967565869),
       theTranMerc_Origin_Lat(0.0),
       theTranMerc_Origin_Long(0.0),
       theTranMerc_False_Northing(0.0),
       theTranMerc_False_Easting(500000.0),
       theTranMerc_Scale_Factor(0.9996),
       theTranMerc_ap(6367449.1458008),
       theTranMerc_bp(16038.508696861),
       theTranMerc_cp(16.832613334334),
       theTranMerc_dp(0.021984404273757),
       theTranMerc_Delta_Easting(40000000.0),
       theTranMerc_Delta_Northing(40000000.0),
       theZone(zone),
       theHemisphere('N')      
 {
    setZone(zone);
    update();
 }

◆ ossimUtmProjection() [4/4]

ossimUtmProjection::ossimUtmProjection ( const ossimUtmProjection & src )

Definition at line 143 of file ossimUtmProjection.cpp.

References setHemisphere(), setZone(), theHemisphere, theZone, and update().

       :ossimMapProjection(src),
       theTranMerc_a(src.theTranMerc_a),
       theTranMerc_f(src.theTranMerc_f),
       theTranMerc_es(src.theTranMerc_es),
       theTranMerc_ebs(src.theTranMerc_ebs),
       theTranMerc_Origin_Lat(src.theTranMerc_Origin_Lat),
       theTranMerc_Origin_Long(src.theTranMerc_Origin_Long),
       theTranMerc_False_Northing(src.theTranMerc_False_Northing),
       theTranMerc_False_Easting(src.theTranMerc_False_Easting),
       theTranMerc_Scale_Factor(src.theTranMerc_Scale_Factor),
       theTranMerc_ap(src.theTranMerc_ap),
       theTranMerc_bp(src.theTranMerc_bp),
       theTranMerc_cp(src.theTranMerc_cp),
       theTranMerc_dp(src.theTranMerc_dp),
       theTranMerc_Delta_Easting(src.theTranMerc_Delta_Easting),
       theTranMerc_Delta_Northing(src.theTranMerc_Delta_Northing),
       theZone(src.theZone),
       theHemisphere(src.theHemisphere)      
 {
    setZone(theZone);
    setHemisphere(theHemisphere);
    update();
 }

Member Function Documentation

◆ computeZone()

ossim_int32 ossimUtmProjection::computeZone ( const ossimGpt & gpt )

static

Definition at line 304 of file ossimUtmProjection.cpp.

References ossimGpt::latd(), ossimGpt::lond(), ossimGpt::lonr(), and M_PI.

Referenced by loadState(), setZone(), and ossimSonomaSensor::updateModel().

 {
    ossim_int32 result = 0;
 
    double longitude = ground.lonr();
    double lat_Degrees  = (ossim_int32)( (ground.latd()) + 0.00000005);
    double long_Degrees = (ossim_int32)( (ground.lond()) + 0.00000005);
    
    if (longitude < M_PI)
       result = (ossim_int32)( (31 + ((180 * longitude) / (6 * M_PI)) ) + 0.00000005);
    else
       result = (ossim_int32)( (((180 * longitude) / (6 * M_PI)) - 29) + 0.00000005);
    if (result > 60)
       result = 1;
     /* UTM special cases */
    if ((lat_Degrees > 55) && (lat_Degrees < 64) && (long_Degrees > -1)
        && (long_Degrees < 3))
       result = 31;
    if ((lat_Degrees > 55) && (lat_Degrees < 64) && (long_Degrees > 2)
        && (long_Degrees < 12))
       result = 32;
    if ((lat_Degrees > 71) && (long_Degrees > -1) && (long_Degrees < 9))
       result = 31;
    if ((lat_Degrees > 71) && (long_Degrees > 8) && (long_Degrees < 21))
       result = 33;
    if ((lat_Degrees > 71) && (long_Degrees > 20) && (long_Degrees < 33))
       result = 35;
    if ((lat_Degrees > 71) && (long_Degrees > 32) && (long_Degrees < 42))
       result = 37;
 
    return result;
 }

◆ computeZoneMeridian()

double ossimUtmProjection::computeZoneMeridian ( ossim_int32 zone )

static

Return in decimal degrees the zone meridian.

Definition at line 337 of file ossimUtmProjection.cpp.

Referenced by loadState(), setZone(), and update().

 {
    return (6.0 * zone - 183.0);;
 }

◆ Convert_Geodetic_To_Transverse_Mercator()

ossim_int32 ossimUtmProjection::Convert_Geodetic_To_Transverse_Mercator	(	double	Latitude,
		double	Longitude,
		double *	Easting,
		double *	Northing
	)		const

private

The function Convert_Geodetic_To_Transverse_Mercator converts geodetic (latitude and longitude) coordinates to Transverse Mercator projection (easting and northing) coordinates, according to the current ellipsoid and Transverse Mercator projection coordinates.

If any errors occur, the error code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is returned.

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

Definition at line 558 of file ossimUtmProjection.cpp.

References M_PI, SPHSN, SPHTMD, theTranMerc_ebs, theTranMerc_False_Easting, theTranMerc_False_Northing, theTranMerc_Origin_Lat, theTranMerc_Origin_Long, theTranMerc_Scale_Factor, TWO_PI, and UTM_NO_ERROR.

Referenced by forward(), and Set_Transverse_Mercator_Parameters().

 {      /* BEGIN Convert_Geodetic_To_Transverse_Mercator */
 
   /*
    * The function Convert_Geodetic_To_Transverse_Mercator converts geodetic
    * (latitude and longitude) coordinates to Transverse Mercator projection
    * (easting and northing) coordinates, according to the current ellipsoid
    * and Transverse Mercator projection coordinates.  If any errors occur, the
    * error code(s) are returned by the function, otherwise UTM_NO_ERROR is
    * returned.
    *
    *    Latitude      : Latitude in radians                         (input)
    *    Longitude     : Longitude in radians                        (input)
    *    Easting       : Easting/X in meters                         (output)
    *    Northing      : Northing/Y in meters                        (output)
    */
 
   double c;       /* Cosine of latitude                          */
   double c2;
   double c3;
   double c5;
   double c7;
   double dlam;    /* Delta longitude - Difference in Longitude       */
   double eta;     /* constant - theTranMerc_ebs *c *c                   */
   double eta2;
   double eta3;
   double eta4;
   double s;       /* Sine of latitude                        */
   double sn;      /* Radius of curvature in the prime vertical       */
   double t;       /* Tangent of latitude                             */
   double tan2;
   double tan3;
   double tan4;
   double tan5;
   double tan6;
   double t1;      /* Term in coordinate conversion formula - GP to Y */
   double t2;      /* Term in coordinate conversion formula - GP to Y */
   double t3;      /* Term in coordinate conversion formula - GP to Y */
   double t4;      /* Term in coordinate conversion formula - GP to Y */
   double t5;      /* Term in coordinate conversion formula - GP to Y */
   double t6;      /* Term in coordinate conversion formula - GP to Y */
   double t7;      /* Term in coordinate conversion formula - GP to Y */
   double t8;      /* Term in coordinate conversion formula - GP to Y */
   double t9;      /* Term in coordinate conversion formula - GP to Y */
   double tmd;     /* True Meridional distance                        */
   double tmdo;    /* True Meridional distance for latitude of origin */
   ossim_int32    Error_Code = UTM_NO_ERROR;
 //  double temp_Origin;
 //  double temp_Long;
 
 //   if ((Latitude < -MAX_LAT) || (Latitude > MAX_LAT))
 //   {  /* Latitude out of range */
 //     Error_Code|= UTM_LAT_ERROR;
 //   }
    if (Longitude > M_PI)
      Longitude -= TWO_PI;
 //   if ((Longitude < (theTranMerc_Origin_Long - MAX_DELTA_LONG))
 //       || (Longitude > (theTranMerc_Origin_Long + MAX_DELTA_LONG)))
 //   {
 //     if (Longitude < 0)
 //       temp_Long = Longitude + TWO_PI;
 //     else
 //       temp_Long = Longitude;
 //     if (theTranMerc_Origin_Long < 0)
 //       temp_Origin = theTranMerc_Origin_Long + TWO_PI;
 //     else
 //       temp_Origin = theTranMerc_Origin_Long;
 //     if ((temp_Long < (temp_Origin - MAX_DELTA_LONG))
 //         || (temp_Long > (temp_Origin + MAX_DELTA_LONG)))
 //       Error_Code|= UTM_LON_ERROR;
 //   }
   if (!Error_Code)
   { /* no errors */
 
     /* 
      *  Delta Longitude
      */
     dlam = Longitude - theTranMerc_Origin_Long;
 
 //     if (fabs(dlam) > (9.0 * M_PI / 180))
 //     { /* Distortion will result if Longitude is more than 9 degrees from the Central Meridian */
 //       Error_Code |= UTM_LON_WARNING;
 //     }
 
      if (dlam > M_PI)
        dlam -= TWO_PI;
      if (dlam < -M_PI)
        dlam += TWO_PI;
      if (fabs(dlam) < 2.e-10)
        dlam = 0.0;
 
     s = sin(Latitude);
     c = cos(Latitude);
     c2 = c * c;
     c3 = c2 * c;
     c5 = c3 * c2;
     c7 = c5 * c2;
     t = tan (Latitude);
     tan2 = t * t;
     tan3 = tan2 * t;
     tan4 = tan3 * t;
     tan5 = tan4 * t;
     tan6 = tan5 * t;
     eta = theTranMerc_ebs * c2;
     eta2 = eta * eta;
     eta3 = eta2 * eta;
     eta4 = eta3 * eta;
 
     /* radius of curvature in prime vertical */
     sn = SPHSN(Latitude);
 
     /* True Meridianal Distances */
     tmd = SPHTMD(Latitude);
 
     /*  Origin  */
     tmdo = SPHTMD (theTranMerc_Origin_Lat);
 
     /* northing */
     t1 = (tmd - tmdo) * theTranMerc_Scale_Factor;
     t2 = sn * s * c * theTranMerc_Scale_Factor/ 2.e0;
     t3 = sn * s * c3 * theTranMerc_Scale_Factor * (5.e0 - tan2 + 9.e0 * eta 
                                                 + 4.e0 * eta2) /24.e0; 
 
     t4 = sn * s * c5 * theTranMerc_Scale_Factor * (61.e0 - 58.e0 * tan2
                                                 + tan4 + 270.e0 * eta - 330.e0 * tan2 * eta + 445.e0 * eta2
                                                 + 324.e0 * eta3 -680.e0 * tan2 * eta2 + 88.e0 * eta4 
                                                 -600.e0 * tan2 * eta3 - 192.e0 * tan2 * eta4) / 720.e0;
 
     t5 = sn * s * c7 * theTranMerc_Scale_Factor * (1385.e0 - 3111.e0 * 
                                                 tan2 + 543.e0 * tan4 - tan6) / 40320.e0;
 
     *Northing = theTranMerc_False_Northing + t1 + pow(dlam,2.e0) * t2
                 + pow(dlam,4.e0) * t3 + pow(dlam,6.e0) * t4
                 + pow(dlam,8.e0) * t5; 
 
     /* Easting */
     t6 = sn * c * theTranMerc_Scale_Factor;
     t7 = sn * c3 * theTranMerc_Scale_Factor * (1.e0 - tan2 + eta ) /6.e0;
     t8 = sn * c5 * theTranMerc_Scale_Factor * (5.e0 - 18.e0 * tan2 + tan4
                                             + 14.e0 * eta - 58.e0 * tan2 * eta + 13.e0 * eta2 + 4.e0 * eta3 
                                             - 64.e0 * tan2 * eta2 - 24.e0 * tan2 * eta3 )/ 120.e0;
     t9 = sn * c7 * theTranMerc_Scale_Factor * ( 61.e0 - 479.e0 * tan2
                                              + 179.e0 * tan4 - tan6 ) /5040.e0;
 
     *Easting = theTranMerc_False_Easting + dlam * t6 + pow(dlam,3.e0) * t7 
                + pow(dlam,5.e0) * t8 + pow(dlam,7.e0) * t9;
   }
   return (Error_Code);
 } /* END OF Convert_Geodetic_To_Transverse_Mercator */

◆ Convert_Transverse_Mercator_To_Geodetic()

ossim_int32 ossimUtmProjection::Convert_Transverse_Mercator_To_Geodetic	(	double	Easting,
		double	Northing,
		double *	Latitude,
		double *	Longitude
	)		const

private

The function Convert_Transverse_Mercator_To_Geodetic converts Transverse Mercator projection (easting and northing) coordinates to geodetic (latitude and longitude) coordinates, according to the current ellipsoid and Transverse Mercator projection parameters.

If any errors occur, the error code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is returned.

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

Definition at line 713 of file ossimUtmProjection.cpp.

References M_PI, RAD_PER_DEG, SPHSN, SPHSR, SPHTMD, theTranMerc_ebs, theTranMerc_False_Easting, theTranMerc_False_Northing, theTranMerc_Origin_Lat, theTranMerc_Origin_Long, theTranMerc_Scale_Factor, TWO_PI, and UTM_NO_ERROR.

Referenced by inverse().

 {      /* BEGIN Convert_Transverse_Mercator_To_Geodetic */
 
   /*
    * The function Convert_Transverse_Mercator_To_Geodetic converts Transverse
    * Mercator projection (easting and northing) coordinates to geodetic
    * (latitude and longitude) coordinates, according to the current ellipsoid
    * and Transverse Mercator projection parameters.  If any errors occur, the
    * error code(s) are returned by the function, otherwise UTM_NO_ERROR is
    * returned.
    *
    *    Easting       : Easting/X in meters                         (input)
    *    Northing      : Northing/Y in meters                        (input)
    *    Latitude      : Latitude in radians                         (output)
    *    Longitude     : Longitude in radians                        (output)
    */
 
   double c;       /* Cosine of latitude                          */
   double de;      /* Delta easting - Difference in Easting (Easting-Fe)    */
   double dlam;    /* Delta longitude - Difference in Longitude       */
   double eta;     /* constant - theTranMerc_ebs *c *c                   */
   double eta2;
   double eta3;
   double eta4;
   double ftphi;   /* Footpoint latitude                              */
   int    i;       /* Loop iterator                   */
   /*double s;        Sine of latitude                        */
   double sn;      /* Radius of curvature in the prime vertical       */
   double sr;      /* Radius of curvature in the meridian             */
   double t;       /* Tangent of latitude                             */
   double tan2;
   double tan4;
   double t10;     /* Term in coordinate conversion formula - GP to Y */
   double t11;     /* Term in coordinate conversion formula - GP to Y */
   double t12;     /* Term in coordinate conversion formula - GP to Y */
   double t13;     /* Term in coordinate conversion formula - GP to Y */
   double t14;     /* Term in coordinate conversion formula - GP to Y */
   double t15;     /* Term in coordinate conversion formula - GP to Y */
   double t16;     /* Term in coordinate conversion formula - GP to Y */
   double t17;     /* Term in coordinate conversion formula - GP to Y */
   double tmd;     /* True Meridional distance                        */
   double tmdo;    /* True Meridional distance for latitude of origin */
   ossim_int32 Error_Code = UTM_NO_ERROR;
 
 //   if ((Easting < (theTranMerc_False_Easting - theTranMerc_Delta_Easting))
 //       ||(Easting > (theTranMerc_False_Easting + theTranMerc_Delta_Easting)))
 //   { /* Easting out of range  */
 //     Error_Code |= UTM_EASTING_ERROR;
 //   }
 //   if ((Northing < (theTranMerc_False_Northing - theTranMerc_Delta_Northing))
 //       || (Northing > (theTranMerc_False_Northing + theTranMerc_Delta_Northing)))
 //   { /* Northing out of range */
 //     Error_Code |= UTM_NORTHING_ERROR;
 //   }
 
   if (!Error_Code)
   {
     /* True Meridional Distances for latitude of origin */
     tmdo = SPHTMD(theTranMerc_Origin_Lat);
 
     /*  Origin  */
     tmd = tmdo +  (Northing - theTranMerc_False_Northing) / theTranMerc_Scale_Factor; 
 
     /* First Estimate */
     sr = SPHSR(0.e0);
     ftphi = tmd/sr;
 
     for (i = 0; i < 5 ; i++)
     {
       t10 = SPHTMD (ftphi);
       sr = SPHSR(ftphi);
       ftphi = ftphi + (tmd - t10) / sr;
     }
 
     /* Radius of Curvature in the meridian */
     sr = SPHSR(ftphi);
 
     /* Radius of Curvature in the meridian */
     sn = SPHSN(ftphi);
 
     /* Sine Cosine terms */
     // s = sin(ftphi);
     c = cos(ftphi);
 
     /* Tangent Value  */
     t = tan(ftphi);
     tan2 = t * t;
     tan4 = tan2 * tan2;
     eta = theTranMerc_ebs * pow(c,2);
     eta2 = eta * eta;
     eta3 = eta2 * eta;
     eta4 = eta3 * eta;
     de = Easting - theTranMerc_False_Easting;
     if (fabs(de) < 0.0001)
       de = 0.0;
 
     /* Latitude */
     t10 = t / (2.e0 * sr * sn * pow(theTranMerc_Scale_Factor, 2));
     t11 = t * (5.e0  + 3.e0 * tan2 + eta - 4.e0 * pow(eta,2)
                - 9.e0 * tan2 * eta) / (24.e0 * sr * pow(sn,3) 
                                        * pow(theTranMerc_Scale_Factor,4));
     t12 = t * (61.e0 + 90.e0 * tan2 + 46.e0 * eta + 45.E0 * tan4
                - 252.e0 * tan2 * eta  - 3.e0 * eta2 + 100.e0 
                * eta3 - 66.e0 * tan2 * eta2 - 90.e0 * tan4
                * eta + 88.e0 * eta4 + 225.e0 * tan4 * eta2
                + 84.e0 * tan2* eta3 - 192.e0 * tan2 * eta4)
           / ( 720.e0 * sr * pow(sn,5) * pow(theTranMerc_Scale_Factor, 6) );
     t13 = t * ( 1385.e0 + 3633.e0 * tan2 + 4095.e0 * tan4 + 1575.e0 
                 * pow(t,6))/ (40320.e0 * sr * pow(sn,7) * pow(theTranMerc_Scale_Factor,8));
     *Latitude = ftphi - pow(de,2) * t10 + pow(de,4) * t11 - pow(de,6) * t12 
                 + pow(de,8) * t13;
 
     t14 = 1.e0 / (sn * c * theTranMerc_Scale_Factor);
 
     t15 = (1.e0 + 2.e0 * tan2 + eta) / (6.e0 * pow(sn,3) * c * 
                                         pow(theTranMerc_Scale_Factor,3));
 
     t16 = (5.e0 + 6.e0 * eta + 28.e0 * tan2 - 3.e0 * eta2
            + 8.e0 * tan2 * eta + 24.e0 * tan4 - 4.e0 
            * eta3 + 4.e0 * tan2 * eta2 + 24.e0 
            * tan2 * eta3) / (120.e0 * pow(sn,5) * c  
                              * pow(theTranMerc_Scale_Factor,5));
 
     t17 = (61.e0 +  662.e0 * tan2 + 1320.e0 * tan4 + 720.e0 
            * pow(t,6)) / (5040.e0 * pow(sn,7) * c 
                           * pow(theTranMerc_Scale_Factor,7));
 
     /* Difference in Longitude */
     dlam = de * t14 - pow(de,3) * t15 + pow(de,5) * t16 - pow(de,7) * t17;
 
     /* Longitude */
     (*Longitude) = theTranMerc_Origin_Long + dlam;
     while (*Latitude > (90.0 * RAD_PER_DEG))
     {
       *Latitude = M_PI - *Latitude;
       *Longitude += M_PI;
       if (*Longitude > M_PI)
         *Longitude -= TWO_PI;
     }
 
      while (*Latitude < (-90.0 * RAD_PER_DEG))
      {
        *Latitude = - (*Latitude + M_PI);
        *Longitude += M_PI;
        if (*Longitude > M_PI)
          *Longitude -= TWO_PI;
      }
      if (*Longitude > TWO_PI)
        *Longitude -= TWO_PI;
      if (*Longitude < -M_PI)
        *Longitude += TWO_PI;
 
 //     if (fabs(dlam) > (9.0 * M_PI / 180))
 //     { /* Distortion will result if Longitude is more than 9 degrees from the Central Meridian */
 //       Error_Code |= UTM_LON_WARNING;
 //     }
   }
   return (Error_Code);
 } /* END OF Convert_Transverse_Mercator_To_Geodetic */

◆ dup()

ossimObject * ossimUtmProjection::dup ( ) const

virtual

Implements ossimProjection.

Definition at line 227 of file ossimUtmProjection.cpp.

References ossimUtmProjection().

 {
    ossimUtmProjection* proj = new ossimUtmProjection(*this);
    return proj;
 }

◆ forward()

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

virtual

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

Implements ossimMapProjection.

Definition at line 205 of file ossimUtmProjection.cpp.

References ossimGpt::changeDatum(), ossimDatum::code(), Convert_Geodetic_To_Transverse_Mercator(), ossimGpt::datum(), ossimGpt::latr(), ossimGpt::lonr(), and ossimMapProjection::theDatum.

Referenced by ossimApplanixUtmModel::updateModel(), and ossimSonomaSensor::updateModel().

 {
    double easting  = 0.0;
    double northing = 0.0;
    ossimGpt gpt = latLon;
    
    if (theDatum)
    {
       if (theDatum->code() != latLon.datum()->code())
       {
          gpt.changeDatum(theDatum); // Shift to our datum.
       }
    }
 
    Convert_Geodetic_To_Transverse_Mercator(gpt.latr(),
                                            gpt.lonr(),
                                            &easting,
                                            &northing);
    
    return ossimDpt(easting, northing);
 }

◆ Get_Transverse_Mercator_Parameters()

void ossimUtmProjection::Get_Transverse_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_Transverse_Mercator_Parameters returns the current ellipsoid and Transverse Mercator projection parameters.

a : Semi-major axis of ellipsoid, in meters (output) f : Flattening of ellipsoid (output) Origin_Latitude : Latitude in radians at the origin of the (output) projection Central_Meridian : Longitude in radians at the center of the (output) projection False_Easting : Easting/X at the center of the projection (output) False_Northing : Northing/Y at the center of the projection (output) Scale_Factor : Projection scale factor (output)

Definition at line 521 of file ossimUtmProjection.cpp.

References theTranMerc_a, theTranMerc_f, theTranMerc_False_Easting, theTranMerc_False_Northing, theTranMerc_Origin_Lat, theTranMerc_Origin_Long, and theTranMerc_Scale_Factor.

 { /* BEGIN Get_Tranverse_Mercator_Parameters  */
   /*
    * The function Get_Transverse_Mercator_Parameters returns the current
    * ellipsoid and Transverse Mercator projection parameters.
    *
    *    a                 : Semi-major axis of ellipsoid, in meters    (output)
    *    f                 : Flattening of ellipsoid                              (output)
    *    Origin_Latitude   : Latitude in radians at the origin of the   (output)
    *                         projection
    *    Central_Meridian  : Longitude in radians at the center of the  (output)
    *                         projection
    *    False_Easting     : Easting/X at the center of the projection  (output)
    *    False_Northing    : Northing/Y at the center of the projection (output)
    *    Scale_Factor      : Projection scale factor                    (output) 
    */
 
   *a = theTranMerc_a;
   *f = theTranMerc_f;
   *Origin_Latitude = theTranMerc_Origin_Lat;
   *Central_Meridian = theTranMerc_Origin_Long;
   *False_Easting = theTranMerc_False_Easting;
   *False_Northing = theTranMerc_False_Northing;
   *Scale_Factor = theTranMerc_Scale_Factor;
   
   return;
 } /* END OF Get_Tranverse_Mercator_Parameters */

◆ getFalseEasting()

double ossimUtmProjection::getFalseEasting ( ) const

virtual

Returns: The false easting.

Reimplemented from ossimMapProjection.

Definition at line 886 of file ossimUtmProjection.cpp.

References theTranMerc_False_Easting.

 {
    return theTranMerc_False_Easting;
 }

◆ getFalseNorthing()

double ossimUtmProjection::getFalseNorthing ( ) const

virtual

Returns: The false northing.

Reimplemented from ossimMapProjection.

Definition at line 891 of file ossimUtmProjection.cpp.

References theTranMerc_False_Northing.

 {
    return theTranMerc_False_Northing;
 }

◆ getHemisphere()

char ossimUtmProjection::getHemisphere ( ) const

Definition at line 238 of file ossimUtmProjection.cpp.

References theHemisphere.

Referenced by ossimEpsgProjectionDatabase::getCodeFromUtmProj(), ossimPdfWriter::getLgiDictProjectionType(), ossimGpkgSpatialRefSysRecord::init(), ossimChipProcTool::newIdentityProjection(), ossimGeoTiff::parseProjection(), and ossimNitfWriterBase::writeGeometry().

 {
    return theHemisphere;
 }

◆ getPcsCode()

ossim_uint32 ossimUtmProjection::getPcsCode ( ) const

virtual

Returns the EPSG PCS code based on datum, zone, and hemisphere.

This overrides ossimMapProjection::getPcsCode(). Sets ossimMapProjection::thePcsCode on success.

Returns: code on success, 0 if datum is not set.

Reimplemented from ossimMapProjection.

Definition at line 927 of file ossimUtmProjection.cpp.

References ossimEpsgProjectionDatabase::getCodeFromUtmProj(), and ossimEpsgProjectionDatabase::instance().

 {
    // Always recompute in case origin, zone, or hemisphere were changed.
    return ossimEpsgProjectionDatabase::instance()->getCodeFromUtmProj( this );
 }

◆ getZone()

ossim_int32 ossimUtmProjection::getZone ( ) const

Definition at line 233 of file ossimUtmProjection.cpp.

References theZone.

Referenced by ossimEpsgProjectionDatabase::getCodeFromUtmProj(), ossimPdfWriter::getLgiDictProjectionType(), ossimGpkgSpatialRefSysRecord::init(), ossimChipProcTool::newIdentityProjection(), ossimGeoTiff::parseProjection(), and ossimNitfWriterBase::writeGeometry().

 {
    return theZone;
 }

◆ inverse()

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

virtual

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

Implements ossimMapProjection.

Definition at line 192 of file ossimUtmProjection.cpp.

References Convert_Transverse_Mercator_To_Geodetic(), DEG_PER_RAD, ossimMapProjection::theDatum, ossimDpt::x, and ossimDpt::y.

Referenced by ossimSonomaSensor::lineSampleHeightToWorld(), ossimSonomaSensor::lineSampleToWorld(), and ossimApplanixUtmModel::updateModel().

 {
    double lat = 0.0;
    double lon = 0.0;
    
    Convert_Transverse_Mercator_To_Geodetic(eastingNorthing.x,
                                            eastingNorthing.y,
                                            &lat,
                                            &lon);
    
    return ossimGpt(lat*DEG_PER_RAD, lon*DEG_PER_RAD, 0.0, theDatum);
 }

◆ loadState()

bool ossimUtmProjection::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 342 of file ossimUtmProjection.cpp.

References ossimString::c_str(), computeZone(), computeZoneMeridian(), ossimKeywordlist::find(), ossimKeywordNames::HEMISPHERE_KW, ossimGpt::latd(), ossimMapProjection::loadState(), ossimGpt::lond(), setHemisphere(), ossimMapProjection::theOrigin, theZone, ossimString::trim(), ossimString::upcase(), update(), and ossimKeywordNames::ZONE_KW.

 {
    const char* zone       = kwl.find(prefix, ossimKeywordNames::ZONE_KW);
    const char* hemisphere = kwl.find(prefix, ossimKeywordNames::HEMISPHERE_KW);
 
    ossimMapProjection::loadState(kwl, prefix);
 
    // initialize zone to a dummy value.
    //
    theZone = 0;
 //   if(ossimString(type) == STATIC_TYPE_NAME(ossimUtmProjection))
    {
       if(!zone)
       {
          theZone = computeZone(theOrigin);
       }
       else if (zone) 
       {
          theZone = atoi(zone);
 
          if( (theZone  < 1) || (theZone > 60) )
          {
             theZone = computeZone(theOrigin);
          }
          else
          {
 //            if(!kwl.find(prefix, ossimKeywordNames::CENTRAL_MERIDIAN_KW))
 //            {
                theOrigin.lond(computeZoneMeridian(theZone));
 //            }
 //            if(!kwl.find(prefix, ossimKeywordNames::ORIGIN_LATITUDE_KW))
 //            {
                theOrigin.latd(0);
 //            }
          }
       }
       if (hemisphere)
       {
          ossimString s = hemisphere;
          s = s.trim();
          s = s.upcase();
          
          setHemisphere(*s.c_str());
       }
       else
       {
          char hemisphere = theOrigin.latd()<0?'S':'N';
          setHemisphere(hemisphere);
       }
    }
 
    update();
    
    return true;
 }

◆ operator==()

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

virtual

Returns TRUE if principal parameters are within epsilon tolerance.

Reimplemented from ossimMapProjection.

Definition at line 899 of file ossimUtmProjection.cpp.

References ossimMapProjection::operator==(), theHemisphere, and theZone.

 {
    bool result = false;
    if ( this == &proj )
    {
       result = true; // Pointer addresses the same.
    }
    else
    {
       //---
       // Check our stuff first.  No sense going onto ossimMapProjection::operator==
       // if we are not a utm projection.
       //---
       const ossimUtmProjection* p = dynamic_cast<const ossimUtmProjection*>(&proj);
       if ( p )
       {
          if ( theZone == p->theZone )
          {
             if ( theHemisphere == p->theHemisphere )
             {
                result = ossimMapProjection::operator==(proj);
             }
          }
       }
    }
    return result;   
 }

◆ print()

std::ostream & ossimUtmProjection::print ( std::ostream & out ) const

virtual

Prints data members to stream.

Returns stream&.

Reimplemented from ossimMapProjection.

Definition at line 876 of file ossimUtmProjection.cpp.

References ossimMapProjection::print(), theHemisphere, and theZone.

 {
    out << setiosflags(ios::fixed) << setprecision(15)
        << "// ossimUtmProjection::print"
        << "\ntheZone:  " << theZone
        << "\ntheHemisphere:  " << theHemisphere
        << endl;
    return ossimMapProjection::print(out);
 }

◆ saveState()

bool ossimUtmProjection::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 399 of file ossimUtmProjection.cpp.

References ossimKeywordlist::add(), ossimKeywordNames::HEMISPHERE_KW, ossimMapProjection::saveState(), theHemisphere, theZone, and ossimKeywordNames::ZONE_KW.

 {
    kwl.add(prefix,
            ossimKeywordNames::ZONE_KW,
            theZone,
            true);
    
    kwl.add(prefix,
            ossimKeywordNames::HEMISPHERE_KW,
            theHemisphere,
            true);
    
    return ossimMapProjection::saveState(kwl, prefix);
 }

◆ Set_Transverse_Mercator_Parameters()

ossim_int32 ossimUtmProjection::Set_Transverse_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_Tranverse_Mercator_Parameters receives the ellipsoid parameters and Tranverse Mercator projection parameters as inputs, and sets the corresponding state variables.

If any errors occur, the error code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is returned.

a : Semi-major axis of ellipsoid, in meters (input) f : Flattening of ellipsoid (input) Origin_Latitude : Latitude in radians at the origin of the (input) projection Central_Meridian : Longitude in radians at the center of the (input) projection False_Easting : Easting/X at the center of the projection (input) False_Northing : Northing/Y at the center of the projection (input) Scale_Factor : Projection scale factor (input)

Definition at line 414 of file ossimUtmProjection.cpp.

References Convert_Geodetic_To_Transverse_Mercator(), M_PI, MAX_DELTA_LONG, MAX_LAT, theTranMerc_a, theTranMerc_ap, theTranMerc_bp, theTranMerc_cp, theTranMerc_Delta_Easting, theTranMerc_Delta_Northing, theTranMerc_dp, theTranMerc_ebs, theTranMerc_ep, theTranMerc_es, theTranMerc_f, theTranMerc_False_Easting, theTranMerc_False_Northing, theTranMerc_Origin_Lat, theTranMerc_Origin_Long, TWO_PI, and UTM_NO_ERROR.

Referenced by update().

 { /* BEGIN Set_Tranverse_Mercator_Parameters */
   /*
    * The function Set_Tranverse_Mercator_Parameters receives the ellipsoid
    * parameters and Tranverse Mercator projection parameters as inputs, and
    * sets the corresponding state variables. If any errors occur, the error
    * code(s) are returned by the function, otherwise UTM_NO_ERROR is
    * returned.
    *
    *    a                 : Semi-major axis of ellipsoid, in meters    (input)
    *    f                 : Flattening of ellipsoid                              (input)
    *    Origin_Latitude   : Latitude in radians at the origin of the   (input)
    *                         projection
    *    Central_Meridian  : Longitude in radians at the center of the  (input)
    *                         projection
    *    False_Easting     : Easting/X at the center of the projection  (input)
    *    False_Northing    : Northing/Y at the center of the projection (input)
    *    Scale_Factor      : Projection scale factor                    (input) 
    */
 
   double tn;        /* True Meridianal distance constant  */
   double tn2;
   double tn3;
   double tn4;
   double tn5;
   double dummy_northing;
   double TranMerc_b; /* Semi-minor axis of ellipsoid, in meters */
 //  double inv_f = 1 / f;
   ossim_int32 Error_Code = UTM_NO_ERROR;
 
 //   if (a <= 0.0)
 //   { /* Semi-major axis must be greater than zero */
 //     Error_Code |= UTM_A_ERROR;
 //   }
 //   if ((inv_f < 250) || (inv_f > 350))
 //   { /* Inverse flattening must be between 250 and 350 */
 //     Error_Code |= UTM_INV_F_ERROR;
 //   }
 //   if ((Origin_Latitude < -MAX_LAT) || (Origin_Latitude > MAX_LAT))
 //   { /* origin latitude out of range */
 //     Error_Code |= UTM_ORIGIN_LAT_ERROR;
 //   }
 //   if ((Central_Meridian < -M_PI) || (Central_Meridian > TWO_PI))
 //   { /* origin longitude out of range */
 //     Error_Code |= UTM_CENT_MER_ERROR;
 //   }
 //   if ((Scale_Factor < MIN_SCALE_FACTOR) || (Scale_Factor > MAX_SCALE_FACTOR))
 //   {
 //     Error_Code |= UTM_SCALE_FACTOR_ERROR;
 //   }
   if (!Error_Code)
   { /* no errors */
     theTranMerc_a = a;
     theTranMerc_f = f;
     theTranMerc_Origin_Lat = 0;
     theTranMerc_Origin_Long = 0;
     theTranMerc_False_Northing = 0;
     theTranMerc_False_Easting = 0; 
     // theTranMerc_Scale_Factor = 1;
 
     /* Eccentricity Squared */
     theTranMerc_es = 2 * theTranMerc_f - theTranMerc_f * theTranMerc_f;
     /* Second Eccentricity Squared */
     theTranMerc_ebs = (1 / (1 - theTranMerc_es)) - 1;
 
     TranMerc_b = theTranMerc_a * (1 - theTranMerc_f);    
     /*True meridianal constants  */
     tn = (theTranMerc_a - TranMerc_b) / (theTranMerc_a + TranMerc_b);
     tn2 = tn * tn;
     tn3 = tn2 * tn;
     tn4 = tn3 * tn;
     tn5 = tn4 * tn;
 
     theTranMerc_ap = theTranMerc_a * (1.e0 - tn + 5.e0 * (tn2 - tn3)/4.e0
                                 + 81.e0 * (tn4 - tn5)/64.e0 );
     theTranMerc_bp = 3.e0 * theTranMerc_a * (tn - tn2 + 7.e0 * (tn3 - tn4)
                                        /8.e0 + 55.e0 * tn5/64.e0 )/2.e0;
     theTranMerc_cp = 15.e0 * theTranMerc_a * (tn2 - tn3 + 3.e0 * (tn4 - tn5 )/4.e0) /16.0;
     theTranMerc_dp = 35.e0 * theTranMerc_a * (tn3 - tn4 + 11.e0 * tn5 / 16.e0) / 48.e0;
     theTranMerc_ep = 315.e0 * theTranMerc_a * (tn4 - tn5) / 512.e0;
     Convert_Geodetic_To_Transverse_Mercator(MAX_LAT,
                                             MAX_DELTA_LONG,
                                             &theTranMerc_Delta_Easting,
                                             &theTranMerc_Delta_Northing);
     Convert_Geodetic_To_Transverse_Mercator(0,
                                             MAX_DELTA_LONG,
                                             &theTranMerc_Delta_Easting,
                                             &dummy_northing);
     theTranMerc_Origin_Lat = Origin_Latitude;
     if (Central_Meridian > M_PI)
       Central_Meridian -= TWO_PI;
     theTranMerc_Origin_Long = Central_Meridian;
     theTranMerc_False_Northing = False_Northing;
     theTranMerc_False_Easting = False_Easting; 
     // theTranMerc_Scale_Factor = Scale_Factor;
   } /* END OF if(!Error_Code) */
   return (Error_Code);
 }  /* END of Set_Transverse_Mercator_Parameters  */

◆ setHemisphere() [1/2]

void ossimUtmProjection::setHemisphere ( const ossimGpt & ground )

Definition at line 275 of file ossimUtmProjection.cpp.

References ossimGpt::latd().

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode(), ossimEpsgProjectionDatabase::createProjFromUtmCode(), ossimChipperUtil::getNewUtmProjection(), loadState(), ossimChipProcTool::newIdentityProjection(), ossimChipProcTool::newUtmProjection(), ossimUtmProjection(), setOrigin(), and ossimSonomaSensor::updateModel().

 {
    char hemisphere = ground.latd()<0.0?'S':'N';
    setHemisphere(hemisphere);
 }

◆ setHemisphere() [2/2]

void ossimUtmProjection::setHemisphere ( char hemisphere )

Definition at line 281 of file ossimUtmProjection.cpp.

References ossimGpt::latd(), theHemisphere, ossimMapProjection::theOrigin, and theTranMerc_False_Northing.

 {
    hemisphere = toupper(hemisphere);
    if((hemisphere != 'N') &&
       (hemisphere != 'S'))
    {
       theHemisphere = theOrigin.latd() < 0?'S':'N';
    }
    else
    {
       theHemisphere = hemisphere;
    }
 
    if (theHemisphere == 'N')
    {
       theTranMerc_False_Northing = 0.0;
    }
    else
    {
       theTranMerc_False_Northing  = 10000000.0;
    }
 }

◆ setOrigin()

void ossimUtmProjection::setOrigin ( const ossimGpt & origin )

virtual

This will set the utm zone and utm origin base on origin passed in.

Note: If the origin latitude is not 0.0 it will also set the hemisphere. If 0.0 hemisphere is unchanged from previous stored value.

Reimplemented from ossimMapProjection.

Definition at line 243 of file ossimUtmProjection.cpp.

References ossimGpt::latd(), ossimMapProjection::origin(), setHemisphere(), ossimMapProjection::setOrigin(), and setZone().

Referenced by ATP::AtpGenerator::constructChain(), and ossimEpsgProjectionFactory::createProjFromAutoCode().

 {
    setZone(origin);
    // NOTE: We will not set the hemisphere if the origin latitude is 0.0.
    if (origin.latd() != 0.0)
    {
       setHemisphere(origin);
    }
    ossimMapProjection::setOrigin(origin);
 }

◆ setZone() [1/2]

void ossimUtmProjection::setZone ( const ossimGpt & ground )

Definition at line 254 of file ossimUtmProjection.cpp.

References computeZone().

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode(), ossimEpsgProjectionDatabase::createProjFromUtmCode(), ossimChipperUtil::getNewUtmProjection(), ossimChipProcTool::newIdentityProjection(), ossimChipProcTool::newUtmProjection(), newUtmView(), ossimUtmProjection(), setOrigin(), ossimApplanixUtmModel::updateModel(), and ossimSonomaSensor::updateModel().

 {
    setZone(computeZone(ground));
 }

◆ setZone() [2/2]

void ossimUtmProjection::setZone ( ossim_int32 zone )

Definition at line 259 of file ossimUtmProjection.cpp.

References computeZone(), computeZoneMeridian(), ossimGpt::datum(), ossimGpt::latd(), ossimGpt::lond(), ossimGpt::lonr(), ossimMapProjection::theDatum, ossimMapProjection::theOrigin, theTranMerc_Origin_Long, and theZone.

 {
    if( (zone < 1) || (zone > 60) )
    {
       theZone = computeZone(theOrigin);
    }
    else
    {
       theZone = zone;
    }
    theOrigin.lond(computeZoneMeridian(theZone));
    theOrigin.latd(0);
    theOrigin.datum(theDatum);
    theTranMerc_Origin_Long = theOrigin.lonr();
 }

◆ update()

void ossimUtmProjection::update ( )

virtual

Reimplemented from ossimMapProjection.

Definition at line 168 of file ossimUtmProjection.cpp.

References computeZoneMeridian(), ossimEllipsoid::getA(), ossimEllipsoid::getFlattening(), ossimGpt::latd(), ossimGpt::latr(), ossimGpt::lond(), ossimGpt::lonr(), ossimMapProjection::origin(), Set_Transverse_Mercator_Parameters(), ossimMapProjection::theEllipsoid, ossimMapProjection::theFalseEastingNorthing, theHemisphere, ossimMapProjection::theOrigin, theTranMerc_False_Easting, theTranMerc_False_Northing, theTranMerc_Scale_Factor, theZone, ossimMapProjection::update(), ossimDpt::x, and ossimDpt::y.

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode(), loadState(), and ossimUtmProjection().

 {
    ossimGpt origin = theOrigin;
    origin.lond(computeZoneMeridian(theZone));
    origin.latd(0.0);
    double falseNorthing = 10000000.0;
    if (theHemisphere == 'N')
    {
       falseNorthing = 0.0;
    }
    Set_Transverse_Mercator_Parameters(theEllipsoid.getA(),
                                       theEllipsoid.getFlattening(),
                                       origin.latr(),
                                       origin.lonr(),
                                       theTranMerc_False_Easting,
                                       falseNorthing,
                                       theTranMerc_Scale_Factor);
 
    theFalseEastingNorthing.x = theTranMerc_False_Easting;
    theFalseEastingNorthing.y = theTranMerc_False_Northing;
 
    ossimMapProjection::update();
 }

Member Data Documentation

◆ theHemisphere

char ossimUtmProjection::theHemisphere

private

can be N or S.

Definition at line 223 of file ossimUtmProjection.h.

Referenced by getHemisphere(), operator==(), ossimUtmProjection(), print(), saveState(), setHemisphere(), and update().

◆ theTranMerc_a

double ossimUtmProjection::theTranMerc_a

private

Definition at line 188 of file ossimUtmProjection.h.

Referenced by Get_Transverse_Mercator_Parameters(), and Set_Transverse_Mercator_Parameters().

◆ theTranMerc_ap

double ossimUtmProjection::theTranMerc_ap

private

Definition at line 205 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_bp

double ossimUtmProjection::theTranMerc_bp

private

Definition at line 206 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_cp

double ossimUtmProjection::theTranMerc_cp

private

Definition at line 207 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_Delta_Easting

double ossimUtmProjection::theTranMerc_Delta_Easting

private

Definition at line 212 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_Delta_Northing

double ossimUtmProjection::theTranMerc_Delta_Northing

private

Definition at line 213 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_dp

double ossimUtmProjection::theTranMerc_dp

private

Definition at line 208 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_ebs

double ossimUtmProjection::theTranMerc_ebs

private

Definition at line 194 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), and Set_Transverse_Mercator_Parameters().

◆ theTranMerc_ep

double ossimUtmProjection::theTranMerc_ep

private

Definition at line 209 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_es

double ossimUtmProjection::theTranMerc_es

private

Definition at line 192 of file ossimUtmProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ theTranMerc_f

double ossimUtmProjection::theTranMerc_f

private

Definition at line 189 of file ossimUtmProjection.h.

Referenced by Get_Transverse_Mercator_Parameters(), and Set_Transverse_Mercator_Parameters().

◆ theTranMerc_False_Easting

double ossimUtmProjection::theTranMerc_False_Easting

private

Definition at line 201 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), getFalseEasting(), Set_Transverse_Mercator_Parameters(), and update().

◆ theTranMerc_False_Northing

double ossimUtmProjection::theTranMerc_False_Northing

private

Definition at line 200 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), getFalseNorthing(), Set_Transverse_Mercator_Parameters(), setHemisphere(), and update().

◆ theTranMerc_Origin_Lat

double ossimUtmProjection::theTranMerc_Origin_Lat

private

Definition at line 198 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), and Set_Transverse_Mercator_Parameters().

◆ theTranMerc_Origin_Long

double ossimUtmProjection::theTranMerc_Origin_Long

private

Definition at line 199 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), Set_Transverse_Mercator_Parameters(), and setZone().

◆ theTranMerc_Scale_Factor

double ossimUtmProjection::theTranMerc_Scale_Factor

private

Definition at line 202 of file ossimUtmProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), and update().

◆ theZone

ossim_int32 ossimUtmProjection::theZone

private

zone can be from 1 through 60 (0 == NOT SET)

Definition at line 218 of file ossimUtmProjection.h.

Referenced by getZone(), loadState(), operator==(), ossimUtmProjection(), print(), saveState(), setZone(), and update().

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ossimUtmProjection() [1/4]

◆ ossimUtmProjection() [2/4]

◆ ossimUtmProjection() [3/4]

◆ ossimUtmProjection() [4/4]

Member Function Documentation

◆ computeZone()

◆ computeZoneMeridian()

◆ Convert_Geodetic_To_Transverse_Mercator()

◆ Convert_Transverse_Mercator_To_Geodetic()

◆ dup()

◆ forward()

◆ Get_Transverse_Mercator_Parameters()

◆ getFalseEasting()

◆ getFalseNorthing()

◆ getHemisphere()

◆ getPcsCode()

◆ getZone()

◆ inverse()

◆ loadState()

◆ operator==()

◆ print()

◆ saveState()

◆ Set_Transverse_Mercator_Parameters()

◆ setHemisphere() [1/2]

◆ setHemisphere() [2/2]

◆ setOrigin()

◆ setZone() [1/2]

◆ setZone() [2/2]

◆ update()

Member Data Documentation

◆ theHemisphere

◆ theTranMerc_a

◆ theTranMerc_ap

◆ theTranMerc_bp

◆ theTranMerc_cp

◆ theTranMerc_Delta_Easting

◆ theTranMerc_Delta_Northing

◆ theTranMerc_dp

◆ theTranMerc_ebs

◆ theTranMerc_ep

◆ theTranMerc_es

◆ theTranMerc_f

◆ theTranMerc_False_Easting

◆ theTranMerc_False_Northing

◆ theTranMerc_Origin_Lat

◆ theTranMerc_Origin_Long

◆ theTranMerc_Scale_Factor

◆ theZone