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

#include <ossimWgs72Datum.h>

Inheritance diagram for ossimWgs72Datum:
ossimThreeParamDatum ossimDatum ossimObject ossimReferenced

Public Member Functions

 ossimWgs72Datum ()
 
virtual ~ossimWgs72Datum ()
 
virtual ossimGpt shift (const ossimGpt &aPt) const
 
virtual ossimGpt shiftToWgs84 (const ossimGpt &aPt) const
 
virtual ossimGpt shiftFromWgs84 (const ossimGpt &aPt) const
 
- Public Member Functions inherited from ossimThreeParamDatum
 ossimThreeParamDatum (const ossimString &code, const ossimString &name, const ossimEllipsoid *anEllipsoid, ossim_float64 sigmaX, ossim_float64 sigmaY, ossim_float64 sigmaZ, ossim_float64 westLongitude, ossim_float64 eastLongitude, ossim_float64 southLatitude, ossim_float64 northLatitude, ossim_float64 param1, ossim_float64 param2, ossim_float64 param3)
 
virtual ~ossimThreeParamDatum ()
 
virtual ossim_float64 param1 () const
 
virtual ossim_float64 param2 () const
 
virtual ossim_float64 param3 () const
 
virtual ossim_float64 param4 () const
 
virtual ossim_float64 param5 () const
 
virtual ossim_float64 param6 () const
 
virtual ossim_float64 param7 () const
 
virtual bool isEqualTo (const ossimObject &obj, ossimCompareType compareType=OSSIM_COMPARE_FULL) const
 
- Public Member Functions inherited from ossimDatum
 ossimDatum (const ossimString &alpha_code, const ossimString &name, const ossimEllipsoid *anEllipsoid, ossim_float64 sigmaX, ossim_float64 sigmaY, ossim_float64 sigmaZ, ossim_float64 westLongitude, ossim_float64 eastLongitude, ossim_float64 southLatitude, ossim_float64 northLatitude)
 
virtual bool isTheSameAs (const ossimDatum *aDatum) const
 
virtual const ossimString & code () const
 
virtual const ossimString & name () const
 
virtual ossim_uint32 epsgCode () const
 
virtual const ossimEllipsoid * ellipsoid () const
 
virtual ossim_float64 sigmaX () const
 
virtual ossim_float64 sigmaY () const
 
virtual ossim_float64 sigmaZ () const
 
virtual ossim_float64 westLongitude () const
 
virtual ossim_float64 eastLongitude () const
 
virtual ossim_float64 southLatitude () const
 
virtual ossim_float64 northLatitude () const
 
virtual bool isInside (ossim_float64 latitude, ossim_float64 longitude) const
 
bool operator== (const ossimDatum &rhs) const
 
bool operator!= (const ossimDatum &rhs) const
 
- Public Member Functions inherited from ossimObject
 ossimObject ()
 
virtual ~ossimObject ()
 
virtual ossimObject * dup () const
 
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 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
 Generic print method. More...
 
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 Attributes

 TYPE_DATA
 

Additional Inherited Members

- Protected Member Functions inherited from ossimDatum
virtual ~ossimDatum ()
 Only friend ossimDatumFactory is permitted to delete. More...
 
virtual void molodenskyShift (double a, double da, double f, double df, double dx, double dy, double dz, double Lat_in, double Lon_in, double Hgt_in, double &Lat_out, double &Lon_out, double &Hgt_out) const
 
bool withinMolodenskyRange (const ossimGpt &pt) const
 
- Protected Member Functions inherited from ossimReferenced
virtual ~ossimReferenced ()
 
- Protected Attributes inherited from ossimThreeParamDatum
ossim_float64 theParam1
 
ossim_float64 theParam2
 
ossim_float64 theParam3
 
 TYPE_DATA
 
- Protected Attributes inherited from ossimDatum
ossimString theCode
 
ossim_uint32 theEpsgCode
 
ossimString theName
 
const ossimEllipsoid * theEllipsoid
 
ossim_float64 theSigmaX
 
ossim_float64 theSigmaY
 
ossim_float64 theSigmaZ
 
ossim_float64 theWestLongitude
 
ossim_float64 theEastLongitude
 
ossim_float64 theSouthLatitude
 
ossim_float64 theNorthLatitude
 
 TYPE_DATA
 

Detailed Description

Definition at line 18 of file ossimWgs72Datum.h.

Constructor & Destructor Documentation

◆ ossimWgs72Datum()

ossimWgs72Datum::ossimWgs72Datum ( )

Definition at line 32 of file ossimWgs72Datum.cpp.

References ossimDatum::ellipsoid().

33  :ossimThreeParamDatum("WGD",
34  "World Geodetic System 1972",
35  ossimEllipsoidFactory::instance()->wgs72(),
36  0.0,
37  0.0,
38  0.0,
39  -180,
40  180,
41  -90,
42  90,
43  0.0,
44  0.0,
45  0.0)
46 {
47  if(!ellipsoid())
48  {
49  //ERROR
50  }
51 }
ossimEllipsoidFactory::instance
static ossimEllipsoidFactory * instance()
Definition: ossimEllipsoidFactory.cpp:17
ossimDatum::ellipsoid
virtual const ossimEllipsoid * ellipsoid() const
Definition: ossimDatum.h:60
ossimThreeParamDatum::ossimThreeParamDatum
ossimThreeParamDatum(const ossimString &code, const ossimString &name, const ossimEllipsoid *anEllipsoid, ossim_float64 sigmaX, ossim_float64 sigmaY, ossim_float64 sigmaZ, ossim_float64 westLongitude, ossim_float64 eastLongitude, ossim_float64 southLatitude, ossim_float64 northLatitude, ossim_float64 param1, ossim_float64 param2, ossim_float64 param3)
Definition: ossimThreeParamDatum.h:35

◆ ~ossimWgs72Datum()

virtual ossimWgs72Datum::~ossimWgs72Datum ( )
inlinevirtual

Definition at line 23 of file ossimWgs72Datum.h.

23 {}

Member Function Documentation

◆ shift()

ossimGpt ossimWgs72Datum::shift ( const ossimGpt &  aPt) const
virtual

Will shift a ground point to this datum. If they are the same datum then nothing happens.

Reimplemented from ossimThreeParamDatum.

Definition at line 53 of file ossimWgs72Datum.cpp.

References ossimGpt::datum(), ossimDatum::ellipsoid(), ossimEllipsoid::getA(), ossimEllipsoid::getB(), ossimGpt::height(), ossimGpt::latd(), ossimGpt::lond(), shiftFromWgs84(), and ossimDatum::shiftToWgs84().

54 {
55  const ossimDatum *aDatum = aPt.datum();
56 
57  if( (ellipsoid()->getA()== aPt.datum()->ellipsoid()->getA())&&
58  (ellipsoid()->getB()== aPt.datum()->ellipsoid()->getB()))
59  {
60  return ossimGpt(aPt.latd(), aPt.lond(), aPt.height(), this);
61  }
62 
63  if(aDatum)
64  {
65  return shiftFromWgs84(aDatum->shiftToWgs84(aPt));
66  }
67 
68  return aPt;
69 }
ossimDatum::shiftToWgs84
virtual ossimGpt shiftToWgs84(const ossimGpt &aPt) const =0
ossimGpt::lond
double lond() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:97
ossimGpt::latd
double latd() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:87
ossimGpt::datum
const ossimDatum * datum() const
datum().
Definition: ossimGpt.h:196
ossimEllipsoid::getA
const double & getA() const
Definition: ossimEllipsoid.h:66
ossimDatum::ellipsoid
virtual const ossimEllipsoid * ellipsoid() const
Definition: ossimDatum.h:60
ossimGpt::height
double height() const
Definition: ossimGpt.h:107
ossimEllipsoid::getB
const double & getB() const
Definition: ossimEllipsoid.h:67
ossimWgs72Datum::shiftFromWgs84
virtual ossimGpt shiftFromWgs84(const ossimGpt &aPt) const
Definition: ossimWgs72Datum.cpp:139

◆ shiftFromWgs84()

ossimGpt ossimWgs72Datum::shiftFromWgs84 ( const ossimGpt &  aPt) const
virtual

Reimplemented from ossimThreeParamDatum.

Definition at line 139 of file ossimWgs72Datum.cpp.

References ossimEllipsoid::a(), DEG_PER_RAD, ossimEllipsoid::flattening(), ossimGpt::height(), ossimEllipsoidFactory::instance(), ossimDatumFactory::instance(), ossimGpt::isHgtNan(), ossimGpt::latd(), ossimGpt::latr(), ossimGpt::lond(), M_PI, ossimNotify(), ossimNotifyLevel_WARN, SEC_PER_RAD, ossimEllipsoidFactory::wgs72(), ossimDatumFactory::wgs72(), ossimEllipsoidFactory::wgs84(), and ossimDatumFactory::wgs84().

Referenced by shift().

140 {
141  /* Begin Geodetic_Shift_WGS84_To_WGS72 */
142  /*
143  * This function shifts a geodetic coordinate (latitude, longitude in radians
144  * and height in meters) relative to WGS84 to a geodetic coordinate
145  * (latitude, longitude in radians and height in meters) relative to WGS72.
146  *
147  * WGS84_Lat : Latitude in radians relative to WGS84 (input)
148  * WGS84_Lon : Longitude in radians relative to WGS84 (input)
149  * WGS84_Hgt : Height in meters relative to WGS84 (input)
150  * WGS72_Lat : Latitude in radians relative to WGS72 (output)
151  * WGS72_Lon : Longitude in radians relative to WGS72 (output)
152  * WGS72_Hgt : Height in meters relative to WGS72 (output)
153  */
154  double Delta_Lat;
155  double Delta_Lon;
156  double Delta_Hgt;
157  double WGS84_a; /* Semi-major axis of WGS84 ellipsoid */
158  double WGS84_f; /* Flattening of WGS84 ellipsoid */
159  double WGS72_a; /* Semi-major axis of WGS72 ellipsoid */
160  double WGS72_f; /* Flattening of WGS72 ellipsoid */
161  double da; /* WGS72_a - WGS84_a */
162  double df; /* WGS72_f - WGS84_f */
163  double Q;
164  double sin_Lat;
165  double sin2_Lat;
166  const ossimDatum *wgs84 = ossimDatumFactory::instance()->wgs84();
167  const ossimDatum *wgs72 = ossimDatumFactory::instance()->wgs72();
168  const ossimEllipsoid *wgs84Ellipsoid = ossimEllipsoidFactory::instance()->wgs84();
169  const ossimEllipsoid *wgs72Ellipsoid = ossimEllipsoidFactory::instance()->wgs72();
170 
171  if(!wgs84 || !wgs72 || !wgs72Ellipsoid || !wgs84Ellipsoid)
172  {
173  ossimNotify(ossimNotifyLevel_WARN) << "ossimWgs72Datum::shiftFromWgs84, NULL pointer found and no shift will be performed\n";
174 
175  return (aPt);
176  }
177 
178  WGS84_a = wgs84Ellipsoid->a();
179  WGS84_f = wgs84Ellipsoid->flattening();
180  WGS72_a = wgs72Ellipsoid->a();
181  WGS72_f = wgs72Ellipsoid->flattening();
182 
183  da = WGS72_a - WGS84_a;
184  df = WGS72_f - WGS84_f;
185  Q = M_PI / 648000;
186  sin_Lat = sin(aPt.latr());
187  sin2_Lat = sin_Lat * sin_Lat;
188 
189  Delta_Lat = (-4.5 * cos(aPt.latr())) / (WGS84_a*Q)
190  + (df * sin(2*aPt.latr())) / Q;
191  Delta_Lat /= SEC_PER_RAD;
192  Delta_Lon = -0.554 / SEC_PER_RAD;
193  Delta_Hgt = -4.5 * sin_Lat + WGS84_a * df * sin2_Lat - da - 1.4;
194 
195  if(aPt.isHgtNan())
196  {
197 
198  return ossimGpt(aPt.latd() + Delta_Lat*DEG_PER_RAD,
199  aPt.lond() + Delta_Lon*DEG_PER_RAD,
200  Delta_Hgt,
201  this);
202  }
203  return ossimGpt(aPt.latd() + Delta_Lat*DEG_PER_RAD,
204  aPt.lond() + Delta_Lon*DEG_PER_RAD,
205  aPt.height() + Delta_Hgt,
206  this);
207 }
ossimEllipsoidFactory::wgs72
const ossimEllipsoid * wgs72() const
Definition: ossimEllipsoidFactory.h:31
ossimEllipsoidFactory::instance
static ossimEllipsoidFactory * instance()
Definition: ossimEllipsoidFactory.cpp:17
ossimGpt::lond
double lond() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:97
DEG_PER_RAD
#define DEG_PER_RAD
Definition: ossimConstants.h:452
ossimEllipsoid::flattening
double flattening() const
Definition: ossimEllipsoid.h:76
SEC_PER_RAD
#define SEC_PER_RAD
Definition: ossimConstants.h:455
ossimGpt::latd
double latd() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:87
ossimDatumFactory::wgs72
const ossimDatum * wgs72() const
Definition: ossimDatumFactory.h:50
ossimGpt::isHgtNan
bool isHgtNan() const
Definition: ossimGpt.h:143
ossimEllipsoidFactory::wgs84
const ossimEllipsoid * wgs84() const
Definition: ossimEllipsoidFactory.h:30
M_PI
#define M_PI
Definition: ossimConstants.h:446
ossimEllipsoid::a
const double & a() const
Definition: ossimEllipsoid.h:63
ossimGpt::height
double height() const
Definition: ossimGpt.h:107
ossimDatumFactory::instance
static ossimDatumFactory * instance()
Definition: ossimDatumFactory.cpp:55
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66
ossimDatumFactory::wgs84
const ossimDatum * wgs84() const
Definition: ossimDatumFactory.h:49
ossimNotify
OSSIMDLLEXPORT std::ostream & ossimNotify(ossimNotifyLevel level=ossimNotifyLevel_WARN)
Definition: ossimNotify.cpp:241

◆ shiftToWgs84()

ossimGpt ossimWgs72Datum::shiftToWgs84 ( const ossimGpt &  aPt) const
virtual

Reimplemented from ossimThreeParamDatum.

Definition at line 71 of file ossimWgs72Datum.cpp.

References ossimEllipsoid::a(), DEG_PER_RAD, ossimEllipsoid::flattening(), ossimGpt::height(), ossimEllipsoidFactory::instance(), ossimDatumFactory::instance(), ossimGpt::isHgtNan(), ossimGpt::latd(), ossimGpt::latr(), ossimGpt::lond(), M_PI, ossimNotify(), ossimNotifyLevel_WARN, SEC_PER_RAD, ossimEllipsoidFactory::wgs72(), ossimDatumFactory::wgs72(), ossimEllipsoidFactory::wgs84(), and ossimDatumFactory::wgs84().

72 {
73 /* Begin Geodetic_Shift_WGS72_To_WGS84 */
74  /*
75  * This function shifts a geodetic coordinate (latitude, longitude in radians
76  * and height in meters) relative to WGS72 to a geodetic coordinate
77  * (latitude, longitude in radians and height in meters) relative to WGS84.
78  *
79  * WGS72_Lat : Latitude in radians relative to WGS72 (input)
80  * WGS72_Lon : Longitude in radians relative to WGS72 (input)
81  * WGS72_Hgt : Height in meters relative to WGS72 (input)
82  * WGS84_Lat : Latitude in radians relative to WGS84 (output)
83  * WGS84_Lon : Longitude in radians relative to WGS84 (output)
84  * WGS84_Hgt : Height in meters relative to WGS84 (output)
85  */
86  double Delta_Lat;
87  double Delta_Lon;
88  double Delta_Hgt;
89  double WGS84_a; /* Semi-major axis of WGS84 ellipsoid */
90  double WGS84_f; /* Flattening of WGS84 ellipsoid */
91  double WGS72_a; /* Semi-major axis of WGS72 ellipsoid */
92  double WGS72_f; /* Flattening of WGS72 ellipsoid */
93  double da; /* WGS84_a - WGS72_a */
94  double df; /* WGS84_f - WGS72_f */
95  double Q;
96  double sin_Lat;
97  double sin2_Lat;
98 
99  const ossimDatum *wgs84 = ossimDatumFactory::instance()->wgs84();
100  const ossimDatum *wgs72 = ossimDatumFactory::instance()->wgs72();
101  const ossimEllipsoid *wgs84Ellipsoid = ossimEllipsoidFactory::instance()->wgs84();
102  const ossimEllipsoid *wgs72Ellipsoid = ossimEllipsoidFactory::instance()->wgs72();
103 
104  if(!wgs84 || !wgs72 || !wgs72Ellipsoid || !wgs84Ellipsoid)
105  {
106  ossimNotify(ossimNotifyLevel_WARN) << "ossimWgs72Datum::shiftToWgs84, NULL pointer found and no shift will be performed\n";
107  return (aPt);
108  }
109 
110  WGS84_a = wgs84Ellipsoid->a();
111  WGS84_f = wgs84Ellipsoid->flattening();
112  WGS72_a = wgs72Ellipsoid->a();
113  WGS72_f = wgs72Ellipsoid->flattening();
114  da = WGS84_a - WGS72_a;
115  df = WGS84_f - WGS72_f;
116  Q = M_PI / 648000;
117  sin_Lat = sin(aPt.latr());
118  sin2_Lat = sin_Lat * sin_Lat;
119 
120  Delta_Lat = (4.5 * cos(aPt.latr())) / (WGS72_a*Q) + (df * sin(2*aPt.latr())) / Q;
121  Delta_Lat /= SEC_PER_RAD;
122  Delta_Lon = 0.554 / SEC_PER_RAD;
123  Delta_Hgt = 4.5 * sin_Lat + WGS72_a * df * sin2_Lat - da + 1.4;
124 
125  if(aPt.isHgtNan())
126  {
127  return ossimGpt(aPt.latd() + Delta_Lat*DEG_PER_RAD,
128  aPt.lond() + Delta_Lon*DEG_PER_RAD,
129  Delta_Hgt,
130  wgs84);
131  }
132  return ossimGpt(aPt.latd() + Delta_Lat*DEG_PER_RAD,
133  aPt.lond() + Delta_Lon*DEG_PER_RAD,
134  aPt.height() + Delta_Hgt,
135  wgs84);
136  /* End Geodetic_Shift_WGS72_To_WGS84 */
137 }
ossimEllipsoidFactory::wgs72
const ossimEllipsoid * wgs72() const
Definition: ossimEllipsoidFactory.h:31
ossimEllipsoidFactory::instance
static ossimEllipsoidFactory * instance()
Definition: ossimEllipsoidFactory.cpp:17
ossimGpt::lond
double lond() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:97
DEG_PER_RAD
#define DEG_PER_RAD
Definition: ossimConstants.h:452
ossimEllipsoid::flattening
double flattening() const
Definition: ossimEllipsoid.h:76
SEC_PER_RAD
#define SEC_PER_RAD
Definition: ossimConstants.h:455
ossimGpt::latd
double latd() const
Will convert the radian measure to degrees.
Definition: ossimGpt.h:87
ossimDatumFactory::wgs72
const ossimDatum * wgs72() const
Definition: ossimDatumFactory.h:50
ossimGpt::isHgtNan
bool isHgtNan() const
Definition: ossimGpt.h:143
ossimEllipsoidFactory::wgs84
const ossimEllipsoid * wgs84() const
Definition: ossimEllipsoidFactory.h:30
M_PI
#define M_PI
Definition: ossimConstants.h:446
ossimEllipsoid::a
const double & a() const
Definition: ossimEllipsoid.h:63
ossimGpt::height
double height() const
Definition: ossimGpt.h:107
ossimDatumFactory::instance
static ossimDatumFactory * instance()
Definition: ossimDatumFactory.cpp:55
ossimGpt::latr
double latr() const
latr().
Definition: ossimGpt.h:66
ossimDatumFactory::wgs84
const ossimDatum * wgs84() const
Definition: ossimDatumFactory.h:49
ossimNotify
OSSIMDLLEXPORT std::ostream & ossimNotify(ossimNotifyLevel level=ossimNotifyLevel_WARN)
Definition: ossimNotify.cpp:241

Member Data Documentation

◆ TYPE_DATA

ossimWgs72Datum::TYPE_DATA

Definition at line 28 of file ossimWgs72Datum.h.

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