GeographicLib
1.37
|
Transverse Mercator projection. More...
#include <GeographicLib/TransverseMercator.hpp>
Public Member Functions | |
TransverseMercator (real a, real f, real k0) | |
void | Forward (real lon0, real lat, real lon, real &x, real &y, real &gamma, real &k) const |
void | Reverse (real lon0, real x, real y, real &lat, real &lon, real &gamma, real &k) const |
void | Forward (real lon0, real lat, real lon, real &x, real &y) const |
void | Reverse (real lon0, real x, real y, real &lat, real &lon) const |
Inspector functions | |
Math::real | MajorRadius () const |
Math::real | Flattening () const |
Math::real | CentralScale () const |
Static Public Member Functions | |
static const TransverseMercator & | UTM () |
Friends | |
class | Ellipsoid |
Transverse Mercator projection.
This uses Krüger's method which evaluates the projection and its inverse in terms of a series. See
Krüger's method has been extended from 4th to 6th order. The maximum error is 5 nm (5 nanometers), ground distance, for all positions within 35 degrees of the central meridian. The error in the convergence is 2 × 10−15" and the relative error in the scale is 6 − 10−12%%. See Sec. 4 of arXiv:1002.1417 for details. The speed penalty in going to 6th order is only about 1%. TransverseMercatorExact is an alternative implementation of the projection using exact formulas which yield accurate (to 8 nm) results over the entire ellipsoid.
The ellipsoid parameters and the central scale are set in the constructor. The central meridian (which is a trivial shift of the longitude) is specified as the lon0 argument of the TransverseMercator::Forward and TransverseMercator::Reverse functions. The latitude of origin is taken to be the equator. There is no provision in this class for specifying a false easting or false northing or a different latitude of origin. However these are can be simply included by the calling function. For example, the UTMUPS class applies the false easting and false northing for the UTM projections. A more complicated example is the British National Grid (EPSG:7405) which requires the use of a latitude of origin. This is implemented by the GeographicLib::OSGB class.
See TransverseMercator.cpp for more information on the implementation.
See Transverse Mercator projection for a discussion of this projection.
Example of use:
TransverseMercatorProj is a command-line utility providing access to the functionality of TransverseMercator and TransverseMercatorExact.
Definition at line 80 of file TransverseMercator.hpp.
GeographicLib::TransverseMercator::TransverseMercator | ( | real | a, |
real | f, | ||
real | k0 | ||
) |
Constructor for a ellipsoid with
[in] | a | equatorial radius (meters). |
[in] | f | flattening of ellipsoid. Setting f = 0 gives a sphere. Negative f gives a prolate ellipsoid. If f > 1, set flattening to 1/f. |
[in] | k0 | central scale factor. |
GeographicErr | if a, (1 − f) a, or k0 is not positive. |
Definition at line 53 of file TransverseMercator.cpp.
References GeographicLib::Math::isfinite(), and GeographicLib::Math::sq().
void GeographicLib::TransverseMercator::Forward | ( | real | lon0, |
real | lat, | ||
real | lon, | ||
real & | x, | ||
real & | y, | ||
real & | gamma, | ||
real & | k | ||
) | const |
Forward projection, from geographic to transverse Mercator.
[in] | lon0 | central meridian of the projection (degrees). |
[in] | lat | latitude of point (degrees). |
[in] | lon | longitude of point (degrees). |
[out] | x | easting of point (meters). |
[out] | y | northing of point (meters). |
[out] | gamma | meridian convergence at point (degrees). |
[out] | k | scale of projection at point. |
No false easting or northing is added. lat should be in the range [−90°, 90°]; lon and lon0 should be in the range [−540°, 540°).
Definition at line 314 of file TransverseMercator.cpp.
References GeographicLib::Math::AngDiff(), GeographicLib::Math::AngNormalize(), GeographicLib::Math::asinh(), GeographicLib::Math::degree(), GeographicLib::Math::hypot(), GeographicLib::Math::pi(), and GeographicLib::Math::sq().
Referenced by GeographicLib::OSGB::Forward(), GeographicLib::UTMUPS::Forward(), and main().
void GeographicLib::TransverseMercator::Reverse | ( | real | lon0, |
real | x, | ||
real | y, | ||
real & | lat, | ||
real & | lon, | ||
real & | gamma, | ||
real & | k | ||
) | const |
Reverse projection, from transverse Mercator to geographic.
[in] | lon0 | central meridian of the projection (degrees). |
[in] | x | easting of point (meters). |
[in] | y | northing of point (meters). |
[out] | lat | latitude of point (degrees). |
[out] | lon | longitude of point (degrees). |
[out] | gamma | meridian convergence at point (degrees). |
[out] | k | scale of projection at point. |
No false easting or northing is added. lon0 should be in the range [−540°, 540°). The value of lon returned is in the range [−180°, 180°).
Definition at line 476 of file TransverseMercator.cpp.
References GeographicLib::Math::AngNormalize(), GeographicLib::Math::degree(), GeographicLib::Math::hypot(), GeographicLib::Math::pi(), and GeographicLib::Math::sq().
Referenced by main(), GeographicLib::OSGB::Reverse(), and GeographicLib::UTMUPS::Reverse().
|
inline |
TransverseMercator::Forward without returning the convergence and scale.
Definition at line 168 of file TransverseMercator.hpp.
|
inline |
TransverseMercator::Reverse without returning the convergence and scale.
Definition at line 177 of file TransverseMercator.hpp.
|
inline |
Definition at line 190 of file TransverseMercator.hpp.
|
inline |
Definition at line 196 of file TransverseMercator.hpp.
|
inline |
Definition at line 210 of file TransverseMercator.hpp.
|
static |
A global instantiation of TransverseMercator with the WGS84 ellipsoid and the UTM scale factor. However, unlike UTM, no false easting or northing is added.
Definition at line 207 of file TransverseMercator.cpp.
References GeographicLib::Constants::UTM_k0(), GeographicLib::Constants::WGS84_a(), and GeographicLib::Constants::WGS84_f().
Referenced by GeographicLib::UTMUPS::Forward(), and GeographicLib::UTMUPS::Reverse().
|
friend |
Definition at line 113 of file TransverseMercator.hpp.