1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
|
#pragma once
/**
* \file NETGeographicLib/TransverseMercator.h
* \brief Header for NETGeographicLib::TransverseMercator class
*
* NETGeographicLib is copyright (c) Scott Heiman (2013)
* GeographicLib is Copyright (c) Charles Karney (2010-2012)
* <charles@karney.com> and licensed under the MIT/X11 License.
* For more information, see
* http://geographiclib.sourceforge.net/
**********************************************************************/
namespace NETGeographicLib
{
/**
* \brief .NET wrapper for GeographicLib::TransverseMercator.
*
* This class allows .NET applications to access GeographicLib::TransverseMercator.
*
* This uses Krüger's method which evaluates the projection and its
* inverse in terms of a series. See
* - L. Krüger,
* <a href="http://dx.doi.org/10.2312/GFZ.b103-krueger28"> Konforme
* Abbildung des Erdellipsoids in der Ebene</a> (Conformal mapping of the
* ellipsoidal earth to the plane), Royal Prussian Geodetic Institute, New
* Series 52, 172 pp. (1912).
* - C. F. F. Karney,
* <a href="http://dx.doi.org/10.1007/s00190-011-0445-3">
* Transverse Mercator with an accuracy of a few nanometers,</a>
* J. Geodesy 85(8), 475--485 (Aug. 2011);
* preprint
* <a href="http://arxiv.org/abs/1002.1417">arXiv:1002.1417</a>.
*
* 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<sup>−15</sup>" and the relative error in the scale
* is 6 − 10<sup>−12</sup>%%. See Sec. 4 of
* <a href="http://arxiv.org/abs/1002.1417">arXiv:1002.1417</a> 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 \e 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 (<a href="http://www.spatialreference.org/ref/epsg/7405/">
* EPSG:7405</a>) which requires the use of a latitude of origin. This is
* implemented by the GeographicLib::OSGB class.
*
* See GeographicLib::TransverseMercator.cpp for more information on the
* implementation.
*
* See \ref transversemercator for a discussion of this projection.
*
* C# Example:
* \include example-TransverseMercator.cs
* Managed C++ Example:
* \include example-TransverseMercator.cpp
* Visual Basic Example:
* \include example-TransverseMercator.vb
*
* <B>INTERFACE DIFFERENCES:</B><BR>
* A default constructor is provided that assumes WGS84 parameters and
* a UTM scale factor.
*
* The MajorRadius, Flattening, and CentralScale functions are
* implemented as properties.
**********************************************************************/
public ref class TransverseMercator
{
private:
// pointer to the unmanaged GeographicLib::TransverseMercator.
const GeographicLib::TransverseMercator* m_pTransverseMercator;
// the finalizer frees the unmanaged memory when the object is destroyed.
!TransverseMercator(void);
public:
/**
* Constructor for a ellipsoid with
*
* @param[in] a equatorial radius (meters).
* @param[in] f flattening of ellipsoid. Setting \e f = 0 gives a sphere.
* Negative \e f gives a prolate ellipsoid. If \e f > 1, set flattening
* to 1/\e f.
* @param[in] k0 central scale factor.
* @exception GeographicErr if \e a, (1 − \e f ) \e a, or \e k0 is
* not positive.
**********************************************************************/
TransverseMercator(double a, double f, double k0);
/**
* The default constructor assumes a WGS84 ellipsoid and a UTM scale
* factor.
**********************************************************************/
TransverseMercator();
/**
* The destructor calls the finalizer.
**********************************************************************/
~TransverseMercator()
{ this->!TransverseMercator(); }
/**
* Forward projection, from geographic to transverse Mercator.
*
* @param[in] lon0 central meridian of the projection (degrees).
* @param[in] lat latitude of point (degrees).
* @param[in] lon longitude of point (degrees).
* @param[out] x easting of point (meters).
* @param[out] y northing of point (meters).
* @param[out] gamma meridian convergence at point (degrees).
* @param[out] k scale of projection at point.
*
* No false easting or northing is added. \e lat should be in the range
* [−90°, 90°]; \e lon and \e lon0 should be in the
* range [−540°, 540°).
**********************************************************************/
void Forward(double lon0, double lat, double lon,
[System::Runtime::InteropServices::Out] double% x,
[System::Runtime::InteropServices::Out] double% y,
[System::Runtime::InteropServices::Out] double% gamma,
[System::Runtime::InteropServices::Out] double% k);
/**
* Reverse projection, from transverse Mercator to geographic.
*
* @param[in] lon0 central meridian of the projection (degrees).
* @param[in] x easting of point (meters).
* @param[in] y northing of point (meters).
* @param[out] lat latitude of point (degrees).
* @param[out] lon longitude of point (degrees).
* @param[out] gamma meridian convergence at point (degrees).
* @param[out] k scale of projection at point.
*
* No false easting or northing is added. \e lon0 should be in the range
* [−540°, 540°). The value of \e lon returned is in
* the range [−180°, 180°).
**********************************************************************/
void Reverse(double lon0, double x, double y,
[System::Runtime::InteropServices::Out] double% lat,
[System::Runtime::InteropServices::Out] double% lon,
[System::Runtime::InteropServices::Out] double% gamma,
[System::Runtime::InteropServices::Out] double% k);
/**
* TransverseMercator::Forward without returning the convergence and scale.
**********************************************************************/
void Forward(double lon0, double lat, double lon,
[System::Runtime::InteropServices::Out] double% x,
[System::Runtime::InteropServices::Out] double% y);
/**
* TransverseMercator::Reverse without returning the convergence and scale.
**********************************************************************/
void Reverse(double lon0, double x, double y,
[System::Runtime::InteropServices::Out] double% lat,
[System::Runtime::InteropServices::Out] double% lon);
/** \name Inspector functions
**********************************************************************/
///@{
/**
* @return \e a the equatorial radius of the ellipsoid (meters). This is
* the value used in the constructor.
**********************************************************************/
property double MajorRadius { double get(); }
/**
* @return \e f the flattening of the ellipsoid. This is the value used in
* the constructor.
**********************************************************************/
property double Flattening { double get(); }
/**
* @return \e k0 central scale for the projection. This is the value of \e
* k0 used in the constructor and is the scale on the central meridian.
**********************************************************************/
property double CentralScale { double get(); }
///@}
};
} // namespace NETGeographicLib
|
