OSSIM - Open Source Software Image Map  Version 1.9.0 (20180803)
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ossimTransMercatorProjection Class Reference

#include <ossimTransMercatorProjection.h>

Inheritance diagram for ossimTransMercatorProjection:
ossimMapProjection ossimProjection ossimObject ossimErrorStatusInterface ossimReferenced

Public Member Functions

 ossimTransMercatorProjection (const ossimEllipsoid &ellipsoid=ossimEllipsoid(6378137, 6356752.3142), const ossimGpt &origin=ossimGpt())
 
 ossimTransMercatorProjection (const ossimEllipsoid &ellipsoid, const ossimGpt &origin, double falseEasting, double falseNorthing, double scaleFactor)
 
 ~ossimTransMercatorProjection ()
 
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 setScaleFactor (double scaleFactor)
 
void setFalseEastingNorthing (double falseEasting, double falseNorthing)
 
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 std::ostream & print (std::ostream &out) const
 
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 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
 

Protected Member Functions

long Set_Transverse_Mercator_Parameters (double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double Scale_Factor)
 
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
 
long Convert_Geodetic_To_Transverse_Mercator (double Latitude, double Longitude, double *Easting, double *Northing) const
 
long Convert_Transverse_Mercator_To_Geodetic (double Easting, double Northing, double *Latitude, double *Longitude) const
 
- Protected Member Functions inherited from ossimMapProjection
virtual ~ossimMapProjection ()
 
void updateFromTransform ()
 
- Protected Member Functions inherited from ossimReferenced
virtual ~ossimReferenced ()
 

Protected Attributes

double TranMerc_a
 
double TranMerc_f
 
double TranMerc_es
 
double TranMerc_ebs
 
double TranMerc_Origin_Lat
 
double TranMerc_Origin_Long
 
double TranMerc_False_Northing
 
double TranMerc_False_Easting
 
double TranMerc_Scale_Factor
 
double TranMerc_ap
 
double TranMerc_bp
 
double TranMerc_cp
 
double TranMerc_dp
 
double TranMerc_ep
 
double TranMerc_Delta_Easting
 
double TranMerc_Delta_Northing
 
- 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 ossimTransMercatorProjection.h.

Constructor & Destructor Documentation

◆ ossimTransMercatorProjection() [1/2]

ossimTransMercatorProjection::ossimTransMercatorProjection ( const ossimEllipsoid &  ellipsoid = ossimEllipsoid(6378137, 6356752.3142),
const ossimGpt &  origin = ossimGpt() 
)

Definition at line 53 of file ossimTransMercatorProjection.cpp.

References update().

55  :
56  ossimMapProjection(ellipsoid, origin),
57  TranMerc_a(6378137.0),
58  TranMerc_f(1.0/298.257223563),
59  TranMerc_es(0.0066943799901413800),
60  TranMerc_ebs(0.0067394967565869),
61  TranMerc_Origin_Lat(origin.latr()),
62  TranMerc_Origin_Long(origin.lonr()),
63  TranMerc_False_Northing(0.0),
64  TranMerc_False_Easting(0.0),
65  TranMerc_Scale_Factor(1.0),
66  TranMerc_ap(6367449.1458008),
67  TranMerc_bp(16038.508696861),
68  TranMerc_cp(16.832613334334),
69  TranMerc_dp(0.021984404273757),
70  TranMerc_Delta_Easting(40000000.0),
71  TranMerc_Delta_Northing(40000000.0)
72 {
73  update();
74 }
ossimMapProjection::origin
virtual ossimGpt origin() const
Definition: ossimMapProjection.cpp:84
ossimMapProjection::ossimMapProjection
ossimMapProjection(const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())
Definition: ossimMapProjection.cpp:40
ossimGpt::lonr
double lonr() const
Returns the longitude in radian measure.
Definition: ossimGpt.h:76
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66

◆ ossimTransMercatorProjection() [2/2]

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

Definition at line 76 of file ossimTransMercatorProjection.cpp.

References update().

81  :
82  ossimMapProjection(ellipsoid, origin),
83  TranMerc_a(6378137.0),
84  TranMerc_f(1.0/298.257223563),
85  TranMerc_es(0.0066943799901413800),
86  TranMerc_ebs(0.0067394967565869),
87  TranMerc_Origin_Lat(origin.latr()),
88  TranMerc_Origin_Long(origin.lonr()),
89  TranMerc_False_Northing(falseNorthing),
90  TranMerc_False_Easting(falseEasting),
91  TranMerc_Scale_Factor(scaleFactor),
92  TranMerc_ap(6367449.1458008),
93  TranMerc_bp(16038.508696861),
94  TranMerc_cp(16.832613334334),
95  TranMerc_dp(0.021984404273757),
96  TranMerc_Delta_Easting(40000000.0),
97  TranMerc_Delta_Northing(40000000.0)
98 {
99  update();
100 }
ossimMapProjection::origin
virtual ossimGpt origin() const
Definition: ossimMapProjection.cpp:84
ossimMapProjection::ossimMapProjection
ossimMapProjection(const ossimEllipsoid &ellipsoid=ossimEllipsoid(), const ossimGpt &origin=ossimGpt())
Definition: ossimMapProjection.cpp:40
ossimGpt::lonr
double lonr() const
Returns the longitude in radian measure.
Definition: ossimGpt.h:76
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66

◆ ~ossimTransMercatorProjection()

ossimTransMercatorProjection::~ossimTransMercatorProjection ( )
inline

Definition at line 28 of file ossimTransMercatorProjection.h.

28 {}

Member Function Documentation

◆ Convert_Geodetic_To_Transverse_Mercator()

long ossimTransMercatorProjection::Convert_Geodetic_To_Transverse_Mercator ( double  Latitude,
double  Longitude,
double *  Easting,
double *  Northing 
) const
protected

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 392 of file ossimTransMercatorProjection.cpp.

References M_PI, SPHSN, SPHTMD, TranMerc_ebs, TranMerc_False_Easting, TranMerc_False_Northing, TRANMERC_NO_ERROR, TranMerc_Origin_Lat, TranMerc_Origin_Long, TranMerc_Scale_Factor, and TWO_PI.

Referenced by forward(), and Set_Transverse_Mercator_Parameters().

397 { /* BEGIN Convert_Geodetic_To_Transverse_Mercator */
398 
399  /*
400  * The function Convert_Geodetic_To_Transverse_Mercator converts geodetic
401  * (latitude and longitude) coordinates to Transverse Mercator projection
402  * (easting and northing) coordinates, according to the current ellipsoid
403  * and Transverse Mercator projection coordinates. If any errors occur, the
404  * error code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is
405  * returned.
406  *
407  * Latitude : Latitude in radians (input)
408  * Longitude : Longitude in radians (input)
409  * Easting : Easting/X in meters (output)
410  * Northing : Northing/Y in meters (output)
411  */
412 
413  double c; /* Cosine of latitude */
414  double c2;
415  double c3;
416  double c5;
417  double c7;
418  double dlam; /* Delta longitude - Difference in Longitude */
419  double eta; /* constant - TranMerc_ebs *c *c */
420  double eta2;
421  double eta3;
422  double eta4;
423  double s; /* Sine of latitude */
424  double sn; /* Radius of curvature in the prime vertical */
425  double t; /* Tangent of latitude */
426  double tan2;
427  double tan3;
428  double tan4;
429  double tan5;
430  double tan6;
431  double t1; /* Term in coordinate conversion formula - GP to Y */
432  double t2; /* Term in coordinate conversion formula - GP to Y */
433  double t3; /* Term in coordinate conversion formula - GP to Y */
434  double t4; /* Term in coordinate conversion formula - GP to Y */
435  double t5; /* Term in coordinate conversion formula - GP to Y */
436  double t6; /* Term in coordinate conversion formula - GP to Y */
437  double t7; /* Term in coordinate conversion formula - GP to Y */
438  double t8; /* Term in coordinate conversion formula - GP to Y */
439  double t9; /* Term in coordinate conversion formula - GP to Y */
440  double tmd; /* True Meridional distance */
441  double tmdo; /* True Meridional distance for latitude of origin */
442  long Error_Code = TRANMERC_NO_ERROR;
443 // double temp_Origin;
444 // double temp_Long;
445 
446 // if ((Latitude < -MAX_LAT) || (Latitude > MAX_LAT))
447 // { /* Latitude out of range */
448 // Error_Code|= TRANMERC_LAT_ERROR;
449 // }
450  if (Longitude > M_PI)
451  Longitude -= TWO_PI;
452 // if ((Longitude < (TranMerc_Origin_Long - MAX_DELTA_LONG))
453 // || (Longitude > (TranMerc_Origin_Long + MAX_DELTA_LONG)))
454 // {
455 // if (Longitude < 0)
456 // temp_Long = Longitude + TWO_PI;
457 // else
458 // temp_Long = Longitude;
459 // if (TranMerc_Origin_Long < 0)
460 // temp_Origin = TranMerc_Origin_Long + TWO_PI;
461 // else
462 // temp_Origin = TranMerc_Origin_Long;
463 // if ((temp_Long < (temp_Origin - MAX_DELTA_LONG))
464 // || (temp_Long > (temp_Origin + MAX_DELTA_LONG)))
465 // Error_Code|= TRANMERC_LON_ERROR;
466 // }
467  if (!Error_Code)
468  { /* no errors */
469 
470  /*
471  * Delta Longitude
472  */
473  dlam = Longitude - TranMerc_Origin_Long;
474 
475 // if (fabs(dlam) > (9.0 * M_PI / 180))
476 // { /* Distortion will result if Longitude is more than 9 degrees from the Central Meridian */
477 // Error_Code |= TRANMERC_LON_WARNING;
478 // }
479 
480  if (dlam > M_PI)
481  dlam -= TWO_PI;
482  if (dlam < -M_PI)
483  dlam += TWO_PI;
484  if (fabs(dlam) < 2.e-10)
485  dlam = 0.0;
486 
487  s = sin(Latitude);
488  c = cos(Latitude);
489  c2 = c * c;
490  c3 = c2 * c;
491  c5 = c3 * c2;
492  c7 = c5 * c2;
493  t = tan (Latitude);
494  tan2 = t * t;
495  tan3 = tan2 * t;
496  tan4 = tan3 * t;
497  tan5 = tan4 * t;
498  tan6 = tan5 * t;
499  eta = TranMerc_ebs * c2;
500  eta2 = eta * eta;
501  eta3 = eta2 * eta;
502  eta4 = eta3 * eta;
503 
504  /* radius of curvature in prime vertical */
505  sn = SPHSN(Latitude);
506 
507  /* True Meridianal Distances */
508  tmd = SPHTMD(Latitude);
509 
510  /* Origin */
511  tmdo = SPHTMD (TranMerc_Origin_Lat);
512 
513  /* northing */
514  t1 = (tmd - tmdo) * TranMerc_Scale_Factor;
515  t2 = sn * s * c * TranMerc_Scale_Factor/ 2.e0;
516  t3 = sn * s * c3 * TranMerc_Scale_Factor * (5.e0 - tan2 + 9.e0 * eta
517  + 4.e0 * eta2) /24.e0;
518 
519  t4 = sn * s * c5 * TranMerc_Scale_Factor * (61.e0 - 58.e0 * tan2
520  + tan4 + 270.e0 * eta - 330.e0 * tan2 * eta + 445.e0 * eta2
521  + 324.e0 * eta3 -680.e0 * tan2 * eta2 + 88.e0 * eta4
522  -600.e0 * tan2 * eta3 - 192.e0 * tan2 * eta4) / 720.e0;
523 
524  t5 = sn * s * c7 * TranMerc_Scale_Factor * (1385.e0 - 3111.e0 *
525  tan2 + 543.e0 * tan4 - tan6) / 40320.e0;
526 
527  *Northing = TranMerc_False_Northing + t1 + pow(dlam,2.e0) * t2
528  + pow(dlam,4.e0) * t3 + pow(dlam,6.e0) * t4
529  + pow(dlam,8.e0) * t5;
530 
531  /* Easting */
532  t6 = sn * c * TranMerc_Scale_Factor;
533  t7 = sn * c3 * TranMerc_Scale_Factor * (1.e0 - tan2 + eta ) /6.e0;
534  t8 = sn * c5 * TranMerc_Scale_Factor * (5.e0 - 18.e0 * tan2 + tan4
535  + 14.e0 * eta - 58.e0 * tan2 * eta + 13.e0 * eta2 + 4.e0 * eta3
536  - 64.e0 * tan2 * eta2 - 24.e0 * tan2 * eta3 )/ 120.e0;
537  t9 = sn * c7 * TranMerc_Scale_Factor * ( 61.e0 - 479.e0 * tan2
538  + 179.e0 * tan4 - tan6 ) /5040.e0;
539 
540  *Easting = TranMerc_False_Easting + dlam * t6 + pow(dlam,3.e0) * t7
541  + pow(dlam,5.e0) * t8 + pow(dlam,7.e0) * t9;
542  }
543  return (Error_Code);
544 } /* END OF Convert_Geodetic_To_Transverse_Mercator */
SPHSN
#define SPHSN(Latitude)
Definition: ossimTransMercatorProjection.cpp:45
SPHTMD
#define SPHTMD(Latitude)
Definition: ossimTransMercatorProjection.cpp:41
M_PI
#define M_PI
Definition: ossimConstants.h:446
TRANMERC_NO_ERROR
#define TRANMERC_NO_ERROR
Definition: ossimTransMercatorProjection.cpp:23
TWO_PI
#define TWO_PI
Definition: ossimConstants.h:449

◆ Convert_Transverse_Mercator_To_Geodetic()

long ossimTransMercatorProjection::Convert_Transverse_Mercator_To_Geodetic ( double  Easting,
double  Northing,
double *  Latitude,
double *  Longitude 
) const
protected

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 547 of file ossimTransMercatorProjection.cpp.

References M_PI, RAD_PER_DEG, SPHSN, SPHSR, SPHTMD, TranMerc_ebs, TranMerc_False_Easting, TranMerc_False_Northing, TRANMERC_NO_ERROR, TranMerc_Origin_Lat, TranMerc_Origin_Long, TranMerc_Scale_Factor, and TWO_PI.

Referenced by inverse().

551 { /* BEGIN Convert_Transverse_Mercator_To_Geodetic */
552 
553  /*
554  * The function Convert_Transverse_Mercator_To_Geodetic converts Transverse
555  * Mercator projection (easting and northing) coordinates to geodetic
556  * (latitude and longitude) coordinates, according to the current ellipsoid
557  * and Transverse Mercator projection parameters. If any errors occur, the
558  * error code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is
559  * returned.
560  *
561  * Easting : Easting/X in meters (input)
562  * Northing : Northing/Y in meters (input)
563  * Latitude : Latitude in radians (output)
564  * Longitude : Longitude in radians (output)
565  */
566 
567  double c; /* Cosine of latitude */
568  double de; /* Delta easting - Difference in Easting (Easting-Fe) */
569  double dlam; /* Delta longitude - Difference in Longitude */
570  double eta; /* constant - TranMerc_ebs *c *c */
571  double eta2;
572  double eta3;
573  double eta4;
574  double ftphi; /* Footpoint latitude */
575  int i; /* Loop iterator */
576  // double s; /* Sine of latitude */
577  double sn; /* Radius of curvature in the prime vertical */
578  double sr; /* Radius of curvature in the meridian */
579  double t; /* Tangent of latitude */
580  double tan2;
581  double tan4;
582  double t10; /* Term in coordinate conversion formula - GP to Y */
583  double t11; /* Term in coordinate conversion formula - GP to Y */
584  double t12; /* Term in coordinate conversion formula - GP to Y */
585  double t13; /* Term in coordinate conversion formula - GP to Y */
586  double t14; /* Term in coordinate conversion formula - GP to Y */
587  double t15; /* Term in coordinate conversion formula - GP to Y */
588  double t16; /* Term in coordinate conversion formula - GP to Y */
589  double t17; /* Term in coordinate conversion formula - GP to Y */
590  double tmd; /* True Meridional distance */
591  double tmdo; /* True Meridional distance for latitude of origin */
592  long Error_Code = TRANMERC_NO_ERROR;
593 
594 // if ((Easting < (TranMerc_False_Easting - TranMerc_Delta_Easting))
595 // ||(Easting > (TranMerc_False_Easting + TranMerc_Delta_Easting)))
596 // { /* Easting out of range */
597 // Error_Code |= TRANMERC_EASTING_ERROR;
598 // }
599 // if ((Northing < (TranMerc_False_Northing - TranMerc_Delta_Northing))
600 // || (Northing > (TranMerc_False_Northing + TranMerc_Delta_Northing)))
601 // { /* Northing out of range */
602 // Error_Code |= TRANMERC_NORTHING_ERROR;
603 // }
604 
605  if (!Error_Code)
606  {
607  /* True Meridional Distances for latitude of origin */
608  tmdo = SPHTMD(TranMerc_Origin_Lat);
609 
610  /* Origin */
611  tmd = tmdo + (Northing - TranMerc_False_Northing) / TranMerc_Scale_Factor;
612 
613  /* First Estimate */
614  sr = SPHSR(0.e0);
615  ftphi = tmd/sr;
616 
617  for (i = 0; i < 5 ; i++)
618  {
619  t10 = SPHTMD (ftphi);
620  sr = SPHSR(ftphi);
621  ftphi = ftphi + (tmd - t10) / sr;
622  }
623 
624  /* Radius of Curvature in the meridian */
625  sr = SPHSR(ftphi);
626 
627  /* Radius of Curvature in the meridian */
628  sn = SPHSN(ftphi);
629 
630  /* Sine Cosine terms */
631  // s = sin(ftphi);
632  c = cos(ftphi);
633 
634  /* Tangent Value */
635  t = tan(ftphi);
636  tan2 = t * t;
637  tan4 = tan2 * tan2;
638  eta = TranMerc_ebs * pow(c,2);
639  eta2 = eta * eta;
640  eta3 = eta2 * eta;
641  eta4 = eta3 * eta;
642  de = Easting - TranMerc_False_Easting;
643  if (fabs(de) < 0.0001)
644  de = 0.0;
645 
646  /* Latitude */
647  t10 = t / (2.e0 * sr * sn * pow(TranMerc_Scale_Factor, 2));
648  t11 = t * (5.e0 + 3.e0 * tan2 + eta - 4.e0 * pow(eta,2)
649  - 9.e0 * tan2 * eta) / (24.e0 * sr * pow(sn,3)
650  * pow(TranMerc_Scale_Factor,4));
651  t12 = t * (61.e0 + 90.e0 * tan2 + 46.e0 * eta + 45.E0 * tan4
652  - 252.e0 * tan2 * eta - 3.e0 * eta2 + 100.e0
653  * eta3 - 66.e0 * tan2 * eta2 - 90.e0 * tan4
654  * eta + 88.e0 * eta4 + 225.e0 * tan4 * eta2
655  + 84.e0 * tan2* eta3 - 192.e0 * tan2 * eta4)
656  / ( 720.e0 * sr * pow(sn,5) * pow(TranMerc_Scale_Factor, 6) );
657  t13 = t * ( 1385.e0 + 3633.e0 * tan2 + 4095.e0 * tan4 + 1575.e0
658  * pow(t,6))/ (40320.e0 * sr * pow(sn,7) * pow(TranMerc_Scale_Factor,8));
659  *Latitude = ftphi - pow(de,2) * t10 + pow(de,4) * t11 - pow(de,6) * t12
660  + pow(de,8) * t13;
661 
662  t14 = 1.e0 / (sn * c * TranMerc_Scale_Factor);
663 
664  t15 = (1.e0 + 2.e0 * tan2 + eta) / (6.e0 * pow(sn,3) * c *
665  pow(TranMerc_Scale_Factor,3));
666 
667  t16 = (5.e0 + 6.e0 * eta + 28.e0 * tan2 - 3.e0 * eta2
668  + 8.e0 * tan2 * eta + 24.e0 * tan4 - 4.e0
669  * eta3 + 4.e0 * tan2 * eta2 + 24.e0
670  * tan2 * eta3) / (120.e0 * pow(sn,5) * c
671  * pow(TranMerc_Scale_Factor,5));
672 
673  t17 = (61.e0 + 662.e0 * tan2 + 1320.e0 * tan4 + 720.e0
674  * pow(t,6)) / (5040.e0 * pow(sn,7) * c
675  * pow(TranMerc_Scale_Factor,7));
676 
677  /* Difference in Longitude */
678  dlam = de * t14 - pow(de,3) * t15 + pow(de,5) * t16 - pow(de,7) * t17;
679 
680  /* Longitude */
681  (*Longitude) = TranMerc_Origin_Long + dlam;
682  while (*Latitude > (90.0 * RAD_PER_DEG))
683  {
684  *Latitude = M_PI - *Latitude;
685  *Longitude += M_PI;
686  if (*Longitude > M_PI)
687  *Longitude -= TWO_PI;
688  }
689 
690  while (*Latitude < (-90.0 * RAD_PER_DEG))
691  {
692  *Latitude = - (*Latitude + M_PI);
693  *Longitude += M_PI;
694  if (*Longitude > M_PI)
695  *Longitude -= TWO_PI;
696  }
697  if (*Longitude > TWO_PI)
698  *Longitude -= TWO_PI;
699  if (*Longitude < -M_PI)
700  *Longitude += TWO_PI;
701 
702 // if (fabs(dlam) > (9.0 * M_PI / 180))
703 // { /* Distortion will result if Longitude is more than 9 degrees from the Central Meridian */
704 // Error_Code |= TRANMERC_LON_WARNING;
705 // }
706  }
707  return (Error_Code);
708 } /* END OF Convert_Transverse_Mercator_To_Geodetic */
SPHSN
#define SPHSN(Latitude)
Definition: ossimTransMercatorProjection.cpp:45
SPHTMD
#define SPHTMD(Latitude)
Definition: ossimTransMercatorProjection.cpp:41
M_PI
#define M_PI
Definition: ossimConstants.h:446
SPHSR
#define SPHSR(Latitude)
Definition: ossimTransMercatorProjection.cpp:48
TRANMERC_NO_ERROR
#define TRANMERC_NO_ERROR
Definition: ossimTransMercatorProjection.cpp:23
TWO_PI
#define TWO_PI
Definition: ossimConstants.h:449
RAD_PER_DEG
#define RAD_PER_DEG
Definition: ossimConstants.h:458

◆ dup()

virtual ossimObject* ossimTransMercatorProjection::dup ( ) const
inlinevirtual

Implements ossimProjection.

Definition at line 30 of file ossimTransMercatorProjection.h.

30 {return new ossimTransMercatorProjection(*this);}
ossimTransMercatorProjection::ossimTransMercatorProjection
ossimTransMercatorProjection(const ossimEllipsoid &ellipsoid=ossimEllipsoid(6378137, 6356752.3142), const ossimGpt &origin=ossimGpt())
Definition: ossimTransMercatorProjection.cpp:53

◆ forward()

ossimDpt ossimTransMercatorProjection::forward ( const ossimGpt &  worldPoint) const
virtual

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

Implements ossimMapProjection.

Definition at line 182 of file ossimTransMercatorProjection.cpp.

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

183 {
184  double easting = 0.0;
185  double northing = 0.0;
186  ossimGpt gpt = latLon;
187 
188  if (theDatum)
189  {
190  if (theDatum->code() != latLon.datum()->code())
191  {
192  gpt.changeDatum(theDatum); // Shift to our datum.
193  }
194  }
195 
196  Convert_Geodetic_To_Transverse_Mercator(gpt.latr(),
197  gpt.lonr(),
198  &easting,
199  &northing);
200 
201  return ossimDpt(easting, northing);
202 }
ossimDatum::code
virtual const ossimString & code() const
Definition: ossimDatum.h:57
ossimGpt::changeDatum
void changeDatum(const ossimDatum *datum)
This will actually perform a shift.
Definition: ossimGpt.cpp:316
ossimGpt::lonr
double lonr() const
Returns the longitude in radian measure.
Definition: ossimGpt.h:76
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66
ossimTransMercatorProjection::Convert_Geodetic_To_Transverse_Mercator
long Convert_Geodetic_To_Transverse_Mercator(double Latitude, double Longitude, double *Easting, double *Northing) const
Definition: ossimTransMercatorProjection.cpp:392
ossimMapProjection::theDatum
const ossimDatum * theDatum
This is only set if we want to have built in datum shifting.
Definition: ossimMapProjection.h:318

◆ Get_Transverse_Mercator_Parameters()

void ossimTransMercatorProjection::Get_Transverse_Mercator_Parameters ( double *  a,
double *  f,
double *  Origin_Latitude,
double *  Central_Meridian,
double *  False_Easting,
double *  False_Northing,
double *  Scale_Factor 
) const
protected

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 355 of file ossimTransMercatorProjection.cpp.

References TranMerc_a, TranMerc_f, TranMerc_False_Easting, TranMerc_False_Northing, TranMerc_Origin_Lat, TranMerc_Origin_Long, and TranMerc_Scale_Factor.

363 { /* BEGIN Get_Tranverse_Mercator_Parameters */
364  /*
365  * The function Get_Transverse_Mercator_Parameters returns the current
366  * ellipsoid and Transverse Mercator projection parameters.
367  *
368  * a : Semi-major axis of ellipsoid, in meters (output)
369  * f : Flattening of ellipsoid (output)
370  * Origin_Latitude : Latitude in radians at the origin of the (output)
371  * projection
372  * Central_Meridian : Longitude in radians at the center of the (output)
373  * projection
374  * False_Easting : Easting/X at the center of the projection (output)
375  * False_Northing : Northing/Y at the center of the projection (output)
376  * Scale_Factor : Projection scale factor (output)
377  */
378 
379  *a = TranMerc_a;
380  *f = TranMerc_f;
381  *Origin_Latitude = TranMerc_Origin_Lat;
382  *Central_Meridian = TranMerc_Origin_Long;
383  *False_Easting = TranMerc_False_Easting;
384  *False_Northing = TranMerc_False_Northing;
385  *Scale_Factor = TranMerc_Scale_Factor;
386 
387  return;
388 } /* END OF Get_Tranverse_Mercator_Parameters */

◆ getFalseEasting()

double ossimTransMercatorProjection::getFalseEasting ( ) const
inlinevirtual
Returns
The false easting.

Reimplemented from ossimMapProjection.

Definition at line 64 of file ossimTransMercatorProjection.h.

64 {return TranMerc_False_Easting;}

◆ getFalseNorthing()

double ossimTransMercatorProjection::getFalseNorthing ( ) const
inlinevirtual
Returns
The false northing.

Reimplemented from ossimMapProjection.

Definition at line 65 of file ossimTransMercatorProjection.h.

65 {return TranMerc_False_Northing;}

◆ getScaleFactor()

double ossimTransMercatorProjection::getScaleFactor ( ) const
inline

Definition at line 66 of file ossimTransMercatorProjection.h.

Referenced by ossimStatePlaneProjectionInfo::matchesProjection(), and ossimGeoTiff::parseProjection().

66 {return TranMerc_Scale_Factor;}

◆ inverse()

ossimGpt ossimTransMercatorProjection::inverse ( const ossimDpt &  projectedPoint) const
virtual

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

Implements ossimMapProjection.

Definition at line 169 of file ossimTransMercatorProjection.cpp.

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

170 {
171  double lat = 0.0;
172  double lon = 0.0;
173 
174  Convert_Transverse_Mercator_To_Geodetic(eastingNorthing.x,
175  eastingNorthing.y,
176  &lat,
177  &lon);
178 
179  return ossimGpt(lat*DEG_PER_RAD, lon*DEG_PER_RAD, 0.0, theDatum);
180 }
DEG_PER_RAD
#define DEG_PER_RAD
Definition: ossimConstants.h:452
ossimTransMercatorProjection::Convert_Transverse_Mercator_To_Geodetic
long Convert_Transverse_Mercator_To_Geodetic(double Easting, double Northing, double *Latitude, double *Longitude) const
Definition: ossimTransMercatorProjection.cpp:547
ossimMapProjection::theDatum
const ossimDatum * theDatum
This is only set if we want to have built in datum shifting.
Definition: ossimMapProjection.h:318

◆ loadState()

bool ossimTransMercatorProjection::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 214 of file ossimTransMercatorProjection.cpp.

References ossimKeywordlist::find(), ossimMapProjection::loadState(), ossimKeywordNames::SCALE_FACTOR_KW, setDefaults(), STATIC_TYPE_NAME, ossimMapProjection::theFalseEastingNorthing, ossimString::toDouble(), TranMerc_False_Easting, TranMerc_False_Northing, TranMerc_Scale_Factor, ossimKeywordNames::TYPE_KW, update(), ossimDpt::x, and ossimDpt::y.

216 {
217  bool flag = ossimMapProjection::loadState(kwl, prefix);
218  const char* type = kwl.find(prefix, ossimKeywordNames::TYPE_KW);
219  const char* scaleFactor = kwl.find(prefix, ossimKeywordNames::SCALE_FACTOR_KW);
220 
221  setDefaults();
222  if(ossimString(type) == STATIC_TYPE_NAME(ossimTransMercatorProjection))
223  {
224  TranMerc_False_Easting = theFalseEastingNorthing.x;
225  TranMerc_False_Northing = theFalseEastingNorthing.y;
226 
227  if(scaleFactor)
228  {
229  double d = ossimString(scaleFactor).toDouble();
230  if (d > 0.0) // Check to avoid divide by zero.
231  {
232  TranMerc_Scale_Factor = d;
233  }
234  }
235  }
236 
237  update();
238 
239  return flag;
240 }
ossimMapProjection::loadState
virtual bool loadState(const ossimKeywordlist &kwl, const char *prefix=0)
Method to the load (recreate) the state of an object from a keyword list.
Definition: ossimMapProjection.cpp:923
ossimKeywordlist::find
const char * find(const char *key) const
Definition: ossimKeywordlist.cpp:599
ossimDpt::y
double y
Definition: ossimDpt.h:165
ossimKeywordNames::TYPE_KW
static const char * TYPE_KW
Definition: ossimKeywordNames.h:190
STATIC_TYPE_NAME
#define STATIC_TYPE_NAME(T)
Definition: ossimRtti.h:325
ossimString::toDouble
double toDouble() const
Definition: ossimString.cpp:870
ossimDpt::x
double x
Definition: ossimDpt.h:164
ossimKeywordNames::SCALE_FACTOR_KW
static const char * SCALE_FACTOR_KW
Definition: ossimKeywordNames.h:165
ossimMapProjection::theFalseEastingNorthing
ossimDpt theFalseEastingNorthing
Hold the false easting northing.
Definition: ossimMapProjection.h:349

◆ operator==()

bool ossimTransMercatorProjection::operator== ( const ossimProjection &  projection) const
virtual

Returns TRUE if principal parameters are within epsilon tolerance.

Reimplemented from ossimMapProjection.

Definition at line 722 of file ossimTransMercatorProjection.cpp.

References ossim::almostEqual(), and TranMerc_Scale_Factor.

723 {
724  if (!ossimMapProjection::operator==(proj))
725  return false;
726 
727  const ossimTransMercatorProjection* p =
728  dynamic_cast<const ossimTransMercatorProjection*>(&proj);
729  if (!p) return false;
730 
731  if (!ossim::almostEqual(TranMerc_Scale_Factor,p->TranMerc_Scale_Factor)) return false;
732 
733  return true;
734 }
ossim::almostEqual
bool almostEqual(T x, T y, T tolerance=FLT_EPSILON)
Definition: ossimCommon.h:53

◆ print()

std::ostream & ossimTransMercatorProjection::print ( std::ostream &  out) const
virtual

Prints data members to stream. Returns stream&.

Reimplemented from ossimMapProjection.

Definition at line 710 of file ossimTransMercatorProjection.cpp.

References ossimMapProjection::print(), ossimKeywordNames::SCALE_FACTOR_KW, and TranMerc_Scale_Factor.

711 {
712  out << setiosflags(ios::fixed) << setprecision(15)
713  << "// ossimTransMercatorProjection::print\n"
714  << ossimKeywordNames::SCALE_FACTOR_KW << ": " << TranMerc_Scale_Factor
715  << endl;
716  return ossimMapProjection::print(out);
717 }
ossimKeywordNames::SCALE_FACTOR_KW
static const char * SCALE_FACTOR_KW
Definition: ossimKeywordNames.h:165
ossimMapProjection::print
virtual std::ostream & print(std::ostream &out) const
Prints data members to stream.
Definition: ossimMapProjection.cpp:1366

◆ saveState()

bool ossimTransMercatorProjection::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 204 of file ossimTransMercatorProjection.cpp.

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

205 {
206  kwl.add(prefix,
207  ossimKeywordNames::SCALE_FACTOR_KW,
208  TranMerc_Scale_Factor,
209  true);
210 
211  return ossimMapProjection::saveState(kwl, prefix);
212 }
ossimKeywordlist::add
void add(const char *prefix, const ossimKeywordlist &kwl, bool overwrite=true)
Definition: ossimKeywordlist.cpp:177
ossimMapProjection::saveState
virtual bool saveState(ossimKeywordlist &kwl, const char *prefix=0) const
Method to save the state of an object to a keyword list.
Definition: ossimMapProjection.cpp:811
ossimKeywordNames::SCALE_FACTOR_KW
static const char * SCALE_FACTOR_KW
Definition: ossimKeywordNames.h:165

◆ Set_Transverse_Mercator_Parameters()

long ossimTransMercatorProjection::Set_Transverse_Mercator_Parameters ( double  a,
double  f,
double  Origin_Latitude,
double  Central_Meridian,
double  False_Easting,
double  False_Northing,
double  Scale_Factor 
)
protected

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 248 of file ossimTransMercatorProjection.cpp.

References Convert_Geodetic_To_Transverse_Mercator(), M_PI, MAX_DELTA_LONG, MAX_LAT, TranMerc_a, TranMerc_ap, TranMerc_bp, TranMerc_cp, TranMerc_Delta_Easting, TranMerc_Delta_Northing, TranMerc_dp, TranMerc_ebs, TranMerc_ep, TranMerc_es, TranMerc_f, TranMerc_False_Easting, TranMerc_False_Northing, TRANMERC_NO_ERROR, TranMerc_Origin_Lat, TranMerc_Origin_Long, TranMerc_Scale_Factor, and TWO_PI.

Referenced by update().

256 { /* BEGIN Set_Tranverse_Mercator_Parameters */
257  /*
258  * The function Set_Tranverse_Mercator_Parameters receives the ellipsoid
259  * parameters and Tranverse Mercator projection parameters as inputs, and
260  * sets the corresponding state variables. If any errors occur, the error
261  * code(s) are returned by the function, otherwise TRANMERC_NO_ERROR is
262  * returned.
263  *
264  * a : Semi-major axis of ellipsoid, in meters (input)
265  * f : Flattening of ellipsoid (input)
266  * Origin_Latitude : Latitude in radians at the origin of the (input)
267  * projection
268  * Central_Meridian : Longitude in radians at the center of the (input)
269  * projection
270  * False_Easting : Easting/X at the center of the projection (input)
271  * False_Northing : Northing/Y at the center of the projection (input)
272  * Scale_Factor : Projection scale factor (input)
273  */
274 
275  double tn; /* True Meridianal distance constant */
276  double tn2;
277  double tn3;
278  double tn4;
279  double tn5;
280  double dummy_northing;
281  double TranMerc_b; /* Semi-minor axis of ellipsoid, in meters */
282 // double inv_f = 1 / f;
283  long Error_Code = TRANMERC_NO_ERROR;
284 
285 // if (a <= 0.0)
286 // { /* Semi-major axis must be greater than zero */
287 // Error_Code |= TRANMERC_A_ERROR;
288 // }
289 // if ((inv_f < 250) || (inv_f > 350))
290 // { /* Inverse flattening must be between 250 and 350 */
291 // Error_Code |= TRANMERC_INV_F_ERROR;
292 // }
293 // if ((Origin_Latitude < -MAX_LAT) || (Origin_Latitude > MAX_LAT))
294 // { /* origin latitude out of range */
295 // Error_Code |= TRANMERC_ORIGIN_LAT_ERROR;
296 // }
297 // if ((Central_Meridian < -M_PI) || (Central_Meridian > TWO_PI))
298 // { /* origin longitude out of range */
299 // Error_Code |= TRANMERC_CENT_MER_ERROR;
300 // }
301 // if ((Scale_Factor < MIN_SCALE_FACTOR) || (Scale_Factor > MAX_SCALE_FACTOR))
302 // {
303 // Error_Code |= TRANMERC_SCALE_FACTOR_ERROR;
304 // }
305  if (!Error_Code)
306  { /* no errors */
307  TranMerc_a = a;
308  TranMerc_f = f;
309  TranMerc_Origin_Lat = 0;
310  TranMerc_Origin_Long = 0;
311  TranMerc_False_Northing = 0;
312  TranMerc_False_Easting = 0;
313  TranMerc_Scale_Factor = 1;
314 
315  /* Eccentricity Squared */
316  TranMerc_es = 2 * TranMerc_f - TranMerc_f * TranMerc_f;
317  /* Second Eccentricity Squared */
318  TranMerc_ebs = (1 / (1 - TranMerc_es)) - 1;
319 
320  TranMerc_b = TranMerc_a * (1 - TranMerc_f);
321  /*True meridianal constants */
322  tn = (TranMerc_a - TranMerc_b) / (TranMerc_a + TranMerc_b);
323  tn2 = tn * tn;
324  tn3 = tn2 * tn;
325  tn4 = tn3 * tn;
326  tn5 = tn4 * tn;
327 
328  TranMerc_ap = TranMerc_a * (1.e0 - tn + 5.e0 * (tn2 - tn3)/4.e0
329  + 81.e0 * (tn4 - tn5)/64.e0 );
330  TranMerc_bp = 3.e0 * TranMerc_a * (tn - tn2 + 7.e0 * (tn3 - tn4)
331  /8.e0 + 55.e0 * tn5/64.e0 )/2.e0;
332  TranMerc_cp = 15.e0 * TranMerc_a * (tn2 - tn3 + 3.e0 * (tn4 - tn5 )/4.e0) /16.0;
333  TranMerc_dp = 35.e0 * TranMerc_a * (tn3 - tn4 + 11.e0 * tn5 / 16.e0) / 48.e0;
334  TranMerc_ep = 315.e0 * TranMerc_a * (tn4 - tn5) / 512.e0;
335  Convert_Geodetic_To_Transverse_Mercator(MAX_LAT,
336  MAX_DELTA_LONG,
337  &TranMerc_Delta_Easting,
338  &TranMerc_Delta_Northing);
339  Convert_Geodetic_To_Transverse_Mercator(0,
340  MAX_DELTA_LONG,
341  &TranMerc_Delta_Easting,
342  &dummy_northing);
343  TranMerc_Origin_Lat = Origin_Latitude;
344  if (Central_Meridian > M_PI)
345  Central_Meridian -= TWO_PI;
346  TranMerc_Origin_Long = Central_Meridian;
347  TranMerc_False_Northing = False_Northing;
348  TranMerc_False_Easting = False_Easting;
349  TranMerc_Scale_Factor = Scale_Factor;
350  } /* END OF if(!Error_Code) */
351  return (Error_Code);
352 } /* END of Set_Transverse_Mercator_Parameters */
M_PI
#define M_PI
Definition: ossimConstants.h:446
MAX_DELTA_LONG
#define MAX_DELTA_LONG
Definition: ossimTransMercatorProjection.cpp:37
TRANMERC_NO_ERROR
#define TRANMERC_NO_ERROR
Definition: ossimTransMercatorProjection.cpp:23
TWO_PI
#define TWO_PI
Definition: ossimConstants.h:449
ossimTransMercatorProjection::Convert_Geodetic_To_Transverse_Mercator
long Convert_Geodetic_To_Transverse_Mercator(double Latitude, double Longitude, double *Easting, double *Northing) const
Definition: ossimTransMercatorProjection.cpp:392

◆ setDefaults()

void ossimTransMercatorProjection::setDefaults ( )

Definition at line 160 of file ossimTransMercatorProjection.cpp.

References TranMerc_Delta_Easting, TranMerc_Delta_Northing, TranMerc_False_Easting, TranMerc_False_Northing, and TranMerc_Scale_Factor.

Referenced by loadState().

161 {
162  TranMerc_False_Easting = 0.0;
163  TranMerc_False_Northing = 0.0;
164  TranMerc_Scale_Factor = 1.0;
165  TranMerc_Delta_Easting = 40000000.0;
166  TranMerc_Delta_Northing = 40000000.0;
167 }

◆ setFalseEasting()

void ossimTransMercatorProjection::setFalseEasting ( double  falseEasting)

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

Definition at line 118 of file ossimTransMercatorProjection.cpp.

References TranMerc_False_Easting, and update().

119 {
120  TranMerc_False_Easting = falseEasting;
121 
122  update();
123 }

◆ setFalseEastingNorthing()

void ossimTransMercatorProjection::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 133 of file ossimTransMercatorProjection.cpp.

References TranMerc_False_Easting, TranMerc_False_Northing, and update().

135 {
136  TranMerc_False_Easting = falseEasting;
137  TranMerc_False_Northing = falseNorthing;
138 
139  update();
140 }

◆ setFalseNorthing()

void ossimTransMercatorProjection::setFalseNorthing ( double  falseNorthing)

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

Definition at line 126 of file ossimTransMercatorProjection.cpp.

References TranMerc_False_Northing, and update().

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode().

127 {
128  TranMerc_False_Northing = falseNorthing;
129 
130  update();
131 }

◆ setParameters()

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

Definition at line 149 of file ossimTransMercatorProjection.cpp.

References TranMerc_False_Easting, TranMerc_False_Northing, TranMerc_Scale_Factor, and update().

152 {
153  TranMerc_False_Easting = falseEasting;
154  TranMerc_False_Northing = falseNorthing;
155  TranMerc_Scale_Factor = scaleFactor;
156 
157  update();
158 }

◆ setScaleFactor()

void ossimTransMercatorProjection::setScaleFactor ( double  scaleFactor)

Definition at line 142 of file ossimTransMercatorProjection.cpp.

References TranMerc_Scale_Factor, and update().

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode(), and ossimWktProjectionFactory::doTransverseMercator().

143 {
144  TranMerc_Scale_Factor = scaleFactor;
145 
146  update();
147 }

◆ update()

void ossimTransMercatorProjection::update ( )
virtual

Reimplemented from ossimMapProjection.

Definition at line 102 of file ossimTransMercatorProjection.cpp.

References ossimEllipsoid::getA(), ossimEllipsoid::getFlattening(), ossimGpt::latr(), ossimGpt::lonr(), Set_Transverse_Mercator_Parameters(), ossimMapProjection::theEllipsoid, ossimMapProjection::theFalseEastingNorthing, ossimMapProjection::theOrigin, TranMerc_False_Easting, TranMerc_False_Northing, TranMerc_Scale_Factor, ossimMapProjection::update(), ossimDpt::x, and ossimDpt::y.

Referenced by ossimEpsgProjectionFactory::createProjFromAutoCode(), loadState(), ossimTransMercatorProjection(), setFalseEasting(), setFalseEastingNorthing(), setFalseNorthing(), setParameters(), and setScaleFactor().

103 {
104  Set_Transverse_Mercator_Parameters(theEllipsoid.getA(),
105  theEllipsoid.getFlattening(),
106  theOrigin.latr(),
107  theOrigin.lonr(),
108  TranMerc_False_Easting,
109  TranMerc_False_Northing,
110  TranMerc_Scale_Factor);
111 
112  theFalseEastingNorthing.x = TranMerc_False_Easting;
113  theFalseEastingNorthing.y = TranMerc_False_Northing;
114 
115  ossimMapProjection::update();
116 }
ossimDpt::y
double y
Definition: ossimDpt.h:165
ossimEllipsoid::getA
const double & getA() const
Definition: ossimEllipsoid.h:66
ossimGpt::lonr
double lonr() const
Returns the longitude in radian measure.
Definition: ossimGpt.h:76
ossimDpt::x
double x
Definition: ossimDpt.h:164
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66
ossimEllipsoid::getFlattening
const double & getFlattening() const
Definition: ossimEllipsoid.h:68
ossimMapProjection::theEllipsoid
ossimEllipsoid theEllipsoid
This method verifies that the projection parameters match the current pcs code.
Definition: ossimMapProjection.h:312
ossimMapProjection::theFalseEastingNorthing
ossimDpt theFalseEastingNorthing
Hold the false easting northing.
Definition: ossimMapProjection.h:349
ossimTransMercatorProjection::Set_Transverse_Mercator_Parameters
long Set_Transverse_Mercator_Parameters(double a, double f, double Origin_Latitude, double Central_Meridian, double False_Easting, double False_Northing, double Scale_Factor)
Definition: ossimTransMercatorProjection.cpp:248

Member Data Documentation

◆ TranMerc_a

double ossimTransMercatorProjection::TranMerc_a
protected

Definition at line 94 of file ossimTransMercatorProjection.h.

Referenced by Get_Transverse_Mercator_Parameters(), and Set_Transverse_Mercator_Parameters().

◆ TranMerc_ap

double ossimTransMercatorProjection::TranMerc_ap
protected

Definition at line 107 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_bp

double ossimTransMercatorProjection::TranMerc_bp
protected

Definition at line 108 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_cp

double ossimTransMercatorProjection::TranMerc_cp
protected

Definition at line 109 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_Delta_Easting

double ossimTransMercatorProjection::TranMerc_Delta_Easting
protected

Definition at line 114 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters(), and setDefaults().

◆ TranMerc_Delta_Northing

double ossimTransMercatorProjection::TranMerc_Delta_Northing
protected

Definition at line 115 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters(), and setDefaults().

◆ TranMerc_dp

double ossimTransMercatorProjection::TranMerc_dp
protected

Definition at line 110 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_ebs

double ossimTransMercatorProjection::TranMerc_ebs
protected

Definition at line 97 of file ossimTransMercatorProjection.h.

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

◆ TranMerc_ep

double ossimTransMercatorProjection::TranMerc_ep
protected

Definition at line 111 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_es

double ossimTransMercatorProjection::TranMerc_es
protected

Definition at line 96 of file ossimTransMercatorProjection.h.

Referenced by Set_Transverse_Mercator_Parameters().

◆ TranMerc_f

double ossimTransMercatorProjection::TranMerc_f
protected

Definition at line 95 of file ossimTransMercatorProjection.h.

Referenced by Get_Transverse_Mercator_Parameters(), and Set_Transverse_Mercator_Parameters().

◆ TranMerc_False_Easting

double ossimTransMercatorProjection::TranMerc_False_Easting
protected

Definition at line 103 of file ossimTransMercatorProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), loadState(), Set_Transverse_Mercator_Parameters(), setDefaults(), setFalseEasting(), setFalseEastingNorthing(), setParameters(), and update().

◆ TranMerc_False_Northing

double ossimTransMercatorProjection::TranMerc_False_Northing
protected

Definition at line 102 of file ossimTransMercatorProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), loadState(), Set_Transverse_Mercator_Parameters(), setDefaults(), setFalseEastingNorthing(), setFalseNorthing(), setParameters(), and update().

◆ TranMerc_Origin_Lat

double ossimTransMercatorProjection::TranMerc_Origin_Lat
protected

Definition at line 100 of file ossimTransMercatorProjection.h.

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

◆ TranMerc_Origin_Long

double ossimTransMercatorProjection::TranMerc_Origin_Long
protected

Definition at line 101 of file ossimTransMercatorProjection.h.

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

◆ TranMerc_Scale_Factor

double ossimTransMercatorProjection::TranMerc_Scale_Factor
protected

Definition at line 104 of file ossimTransMercatorProjection.h.

Referenced by Convert_Geodetic_To_Transverse_Mercator(), Convert_Transverse_Mercator_To_Geodetic(), Get_Transverse_Mercator_Parameters(), loadState(), operator==(), print(), saveState(), Set_Transverse_Mercator_Parameters(), setDefaults(), setParameters(), setScaleFactor(), and update().

The documentation for this class was generated from the following files:
Generated on Fri Aug 3 2018 08:47:16 for OSSIM - Open Source Software Image Map by   doxygen 1.8.14