Polya urn model: Difference between revisions
en>Felsenst |
|||
Line 1: | Line 1: | ||
The article [[Transverse Mercator projection]] restricts itself to general features of the projection. This article describes in detail one of the (two) implementations developed by Louis Krüger in 1912;<ref name=kruger>Krüger, L. (1912). ''[http://dx.doi.org/10.2312/GFZ.b103-krueger28 Konforme Abbildung des Erdellipsoids in der Ebene]''. Royal Prussian Geodetic Institute, New Series 52.</ref> that expressed as a power series in the longitude difference from the central meridian. These series were recalculated by Lee in 1946,<ref name=lee_series>Lee L P, (1946). Survey Review, Volume '''8''' (Part 58), pp 142–152. [http://www.ingentaconnect.com/content/maney/sre/1945/00000008/00000058/art00004 ''The transverse Mercator projection of the spheroid'']. (Errata and comments in Volume '''8''' (Part 61), pp 277–278.</ref> by Redfearn in 1948,<ref name=Redfearn>Redfearn, J C B (1948). Survey Review, Volume '''9''' (Part 69), pp 318–322, [http://www.ingentaconnect.com/content/maney/sre/1948/00000009/00000069/art00005 ''Transverse Mercator formulae''].</ref> and by Thomas in 1952.<ref>Thomas, Paul D (1952). ''Conformal Projections in Geodesy and Cartography''. Washington: U.S. Coast and Geodetic Survey Special Publication 251.</ref><ref name=snyder /> They are often referred to as the Redfearn series, or the Thomas series. This implementation is of great importance since it is widely used in the U.S. State Plane Coordinate System,<ref name=snyder /> in national (Britain,<ref name=osgb/> Ireland<ref>See [[Irish grid reference system]] and [[Irish Transverse Mercator]]</ref> <!--Germany,<ref>australian</ref> Australia<ref>australian</ref>--> and many others) and also international<ref>{{cite web | url = http://eusoils.jrc.ec.europa.eu/Projects/Alpsis/Docs/ref_grid_sh_proc_draft.pdf | title = Short Proceedings of the 1st European Workshop on Reference Grids, Ispra, 27–29 October 2003 | publisher = [[European Environment Agency]] | date = 2004-06-14 | accessdate = 2009-08-27 | page = 6}}The EEA recommends the Transverse Mercator for conformal pan-European mapping at scales larger than 1:500,000</ref> mapping systems, including the [[Universal Transverse Mercator coordinate system]] (UTM).<ref name=utm1>Defense Mapping Agency Technical Report TM 8358.1: Datums, Ellipsoids, Grids and Grid Reference Systems. URL: http://earth-info.nga.mil/GandG/publications/</ref><ref name=utm2>J. W. Hager, J.F. Behensky, and B.W. Drew, 1989. Defense Mapping Agency Technical Report TM 8358.2. The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS). URL http://earth-info.nga.mil/GandG/publications/</ref> They are also incorporated into the Geotrans coordinate converter made available by the United States National Geospatial-Intelligence Agency.<ref>Geotrans, 2010, Geographic translator, version 3.0, URL http://earth-info.nga.mil/GandG/publications/</ref> When paired with a suitable [[Datum (geodesy)|geodetic datum]], the series deliver high accuracy in zones less than a few degrees in east-west extent. | |||
{{TOC limit|2}} | |||
<!--=============================================================--> | |||
==Preliminaries I: datum and ellipsoid parameters== | |||
The series must be used with a [[Datum (geodesy)|geodetic datum]] which specifies the position, orientation and shape of a [[Figure of the Earth|Reference ellipsoid]]. Although the projection formulae depend only on the shape parameters of the reference ellipsoid the full set of datum parameters is necessary to link the projection coordinates to true positions in three dimensional space. The datums and reference ellipsoids associated with particular implementations of the Redfearn formulae are listed [[#Implementations|below]]. A comprehensive list of important ellipsoids is given in the article on the [[Figure of the Earth#table|Figure of the Earth]]. | |||
{{space|5}}In specifying ellipsoids it is normal to give the [[semi-major axis]] (equatorial axis), <math>a</math>, along with either the [[Flattening|inverse flattening]], <math>1/f</math>, or the [[semi-minor axis]] (polar axis), <math>b</math>, or sometimes both. The series presented below use the eccentricity, <math>e</math>, in preference to the flattening, <math>f</math>. In addition they use the parameters <math>n</math>, called the [[Flattening#First and second flattening|third flattening]], and <math>e'</math>, the [[Eccentricity (mathematics)|second eccentricity]]. There are only two independent shape parameters and there are many relations between them: in particular | |||
::<math> | |||
\begin{align} | |||
f&=\frac{a-b}{a}, \qquad e^2=2f-f^2, \qquad e'^2=\frac{e^2}{1-e^2}\\ | |||
b&=a(1-f)=a(1-e^2)^{1/2},\qquad n=\frac{a-b}{a+b}. | |||
\end{align} | |||
</math> | |||
The projection formulae also involve <math>\rho(\phi)</math>, the [[Radius of curvature (applications)|radius of curvature]] of the meridian (at latitude <math>\phi</math>), and <math>\nu(\phi) </math>, the radius of curvature in the [[prime vertical]]. | |||
(The prime vertical is the vertical plane orthogonal to the meridian plane at a point on the ellipsoid). The radii of curvature are defined as follows: | |||
:<math> | |||
\nu(\phi)=\frac{a}{\sqrt{1-e^2\sin^2\phi}}, | |||
\qquad | |||
\rho(\phi)=\frac{\nu^3(1-e^2)}{a^2}. | |||
</math> | |||
In addition the functions <math>\beta(\phi)</math> and <math>\eta(\phi)</math> are defined as: | |||
:<math> | |||
\beta(\phi)=\frac{\nu(\phi)}{\rho(\phi)},\qquad | |||
\eta^2=\beta-1={e'^2\cos^2\!\phi}. | |||
</math> | |||
For compactness it is normal to introduce the following abbreviations: | |||
::<math> | |||
s=\sin\phi,\qquad | |||
c=\cos\phi,\qquad | |||
t=\tan\phi. | |||
</math> | |||
<!--=============================================================--> | |||
==Preliminaries II: meridian distance== | |||
===Meridian distance=== | |||
The article on [[Meridian arc]] describes several methods of computing <math>m(\phi)</math>, the meridian distance from the equator to a point at latitude <math>\phi</math> : the expressions given below are those used in the '''actual'' implementation of the Transverse Mercator projection by the OSGB.<ref name=osgb /> The truncation error is less than 0.1mm so the series is certainly accurate to within 1mm, the design tolerance of the OSGB implementation. | |||
::<math> | |||
\begin{align} | |||
m(\phi)&=B_0\phi+B_2\sin 2\phi+B_4\sin4\phi+B_6\sin6\phi+\cdots, | |||
\end{align} | |||
</math> | |||
where the coefficients are given to order <math>n^3</math> (order <math>e^6</math>) by | |||
::<math> | |||
\begin{align} | |||
B_0 &= b\bigg(1+n+\frac{5}{4}n^2+\frac{5}{4}n^3 \bigg), | |||
\qquad | |||
B_4 = b\bigg(\frac{15}{16}n^2+\frac{15}{16}n^3 \bigg),\\ | |||
B_2 &= - b\bigg(\frac{3}{2}n+\frac{3}{2}n^2+\frac{21}{16}n^3 \bigg), | |||
\qquad | |||
B_6 = - b\bigg(\frac{35}{48}n^3\bigg). | |||
\end{align} | |||
</math> | |||
The meridian distance from equator to pole is | |||
::<math>m_p=m(\pi/2)=\pi B_0/2\,.</math> | |||
The form of the series specified for UTM is a variant of the above exhibiting higher order terms with a truncation error of 0.03mm. | |||
<!--=============================================================--> | |||
====Inverse meridian distance==== | |||
Neither the OSGB nor the UTM implementations define an inverse series for the meridian distance; instead they use an iterative scheme. For a given meridian distance <math>M</math> first set <math>\phi_0=M/B_0</math> and then iterate using | |||
::<math> | |||
\begin{align} | |||
\phi_n=\phi_{n-1}+ \frac{M-m(\phi_{n-1})}{B_0}, \qquad n=1,2,3,\ldots | |||
\end{align} | |||
</math> | |||
until <math>|M-m(\phi_{n-1})|<0.01</math>mm. | |||
{{space|5}}The inversion ''can'' be effected by a series, presented here for later reference. For a given meridian distance, <math>M</math>, define the [[latitude|rectifying latitude]] by | |||
::<math> | |||
\mu=\frac{\pi M}{2 m_p}. | |||
</math> | |||
The geodetic latitude corresponding to <math>M</math> is (Snyder<ref name=snyder /> page 17)<!--(Krüger,<ref name=kruger /> page 12–13)-->: | |||
::<math> | |||
\begin{align} | |||
\phi&=\mu+D_2\sin 2\mu+D_4\sin4\mu+D_6\sin6\mu+D_8\sin8\mu+\cdots,\\ | |||
\end{align} | |||
</math> | |||
where, to <math>O(n^4)</math>, | |||
::<math> | |||
\begin{align} | |||
D_2 & = \frac{3}{2}n-\frac{27}{32}n^3, & | |||
D_4 & = \frac{21}{16}n^2-\frac{55}{32}n^4,\\[8pt] | |||
D_6 & = \frac{151}{96}n^3, & | |||
D_8 & = \frac{1097}{512}n^4. | |||
\end{align} | |||
</math> | |||
<!--=============================================================--> | |||
==An outline of the method== | |||
<!--=============================================================--> | |||
The normal aspect of the Mercator projection of a sphere of radius <math>R</math> is described by the equations | |||
:<math> | |||
x = R\lambda, \qquad\qquad y = R\psi, | |||
</math> | |||
where <math>\psi</math>, the [[Latitude#Isometric latitude|isometric latitude]], is given by | |||
::<math> | |||
\begin{align} | |||
\psi &= \ln \left[\tan \left(\frac{\pi}{4} + \frac{\phi}{2} \right) \right]. | |||
\end{align} | |||
</math> | |||
<!-- & = \frac {1}{2} \ln \left[ \frac {1 + \sin\phi}{1 - \sin\phi} \right] \\[2ex] | |||
& = \sinh^{-1} \left( \tan\phi\right) | |||
& = \tanh^{-1} \left( \sin\phi\right)--> | |||
On the ellipsoid the isometric latitude becomes | |||
::<math> | |||
\begin{align} | |||
\psi &= \ln \left[\tan \left(\frac{\pi}{4} + \frac{\phi}{2} \right) \right] | |||
-\frac{e}{2} \ln \left[ \frac {1 + e\sin\phi}{1 - e\sin\phi} \right]. | |||
\end{align} | |||
</math> | |||
By construction, the projection from the geodetic coordinates (<math>\phi</math>,<math>\lambda</math>) to the coordinates (<math>\psi</math>,<math>\lambda</math>) is conformal. If the coordinates (<math>\psi</math>,<math>\lambda</math>) are used to define a point <math>\zeta=\psi+i\lambda</math> in the complex plane, then any analytic function <math>f(\zeta)</math> will define another conformal projection. Kruger's method involves seeking the specific <math>f(\zeta)</math> which generates a uniform scale along the central meridian, <math>\lambda=0</math>. He achieved this by investigating a Taylor series approximation with the projection coordinates given by: | |||
::<math> | |||
\begin{align} | |||
y+ix&=f(\zeta)=f(\psi+i\lambda)\\ | |||
&= f(\psi+i.0) + A_1\lambda+ A_2\lambda^2+ A_3\lambda^3+ \ldots, | |||
\end{align} | |||
</math> | |||
where the real part of <math>f(\psi+i.0)</math> must be proportional to the meridian distance function <math>m(\phi)</math>. The (complex) coefficients <math>A_n</math> depend on derivatives of <math>f(\zeta)</math> which can be reduced to derivatives of <math>m(\phi)</math> with respect to <math>\psi</math>, (not <math>\phi</math>). The derivatives are straightforward to evaluate in principle but the expressions become very involved at high orders because of the complicated relation between <math>\psi</math> and <math>\phi</math>. Separation of real and imaginary parts gives the series for <math>x</math> and <math>y</math> and further derivatives give the scale and convergence factors. | |||
<!--=============================================================--> | |||
==The series in detail== | |||
This section presents the eighth order series as published by Redfearn<ref name=Redfearn/> (but with <math>x</math> and <math>y</math> interchanged and the longitude difference from the central meridian denoted by <math>\lambda</math> instead of <math>\omega</math>). Equivalent eighth order series, with different notations, can be found in Snyder<ref name=snyder>{{Cite book| author=Snyder, John P. | title=Map Projections – A Working Manual. U.S. Geological Survey Professional Paper 1395 | publisher =United States Government Printing Office, Washington, D.C. | year=1987}}This paper can be downloaded from [http://pubs.er.usgs.gov/pubs/pp/pp1395 USGS pages.] It gives full details of most projections, together with interesting introductory sections, but it does not derive any of the projections from first principles.</ref> (pages 60–64) and at many web sites such as that for the Ordnance Survey of Great Britain.<ref name=osgb>A guide to coordinate systems in Great Britain. This is available as a pdf document at http://www.ordnancesurvey.co.uk/oswebsite/gps/information/coordinatesystemsinfo/guidecontents/</ref> | |||
{{space|5}}The direct series are developed in terms of the longitude difference from the central meridian, expressed in radians: the inverse series are developed in terms of the ratio <math>x/a</math>. The projection is normally restricted to narrow zones (in longitude) so that both of the expansion parameters are typically less than about 0.1, guaranteeing rapid [[Convergence (mathematics)|convergence]]. For example in each [[Universal Transverse Mercator coordinate system|UTM]] zone these expansion parameters are less than 0.053 and for the British national grid ([[British national grid reference system|NGGB]]) they are less than 0.09. All of the direct series giving <math>x</math>, <math>y</math>, scale <math>k</math>, convergence <math>\gamma</math> are functions of both latitude and longitude and the parameters of the ellipsoid: all inverse series giving <math>\phi</math>, <math>\lambda</math>, <math>k</math>, <math>\gamma</math> are functions of both <math>x</math> and <math>y</math> and the parameters of the ellipsoid. | |||
===Direct series=== | |||
In the following series <math>\lambda</math> is the ''difference'' of the longitude of an arbitrary point and the longitude of the chosen central meridian: <math>\lambda</math> is in radians and is positive east of the central meridian. The W coefficients are functions of <math>\phi</math> listed [[#The coefficients for all series|below]]. The series for <math>y</math> reduces to the scaled meridian distance when <math>\lambda=0</math>. | |||
::<math> | |||
\begin{align} | |||
x(\lambda,\phi)&=k_0\nu\left[\lambda c | |||
+\frac{\lambda^3c^3 W_3}{3!} | |||
+\frac{\lambda^5c^5 W_5}{5!} | |||
+\frac{\lambda^7c^7 W_7}{7!}\right] , \\[1ex] | |||
y(\lambda,\phi)&=k_0\left[m(\phi) +\frac{\lambda^2 \nu c^2t}{2} | |||
+\frac{\lambda^4 \nu c^4tW_4}{4!} | |||
+\frac{\lambda^6 \nu c^6tW_6}{6!} | |||
+\frac{\lambda^8 \nu c^8tW_8}{8!} \right], | |||
\end{align} | |||
</math> | |||
<!--=============================================================--> | |||
===Inverse series=== | |||
The inverse series involve a further construct: the '''footpoint latitude'''. Given a point <math>(x,y)</math> on the projection the '''footpoint''' is defined as the point on the central meridian with coordinates <math>(0,y)</math>. Since the scale on the central meridian is <math>k_0</math> the meridian distance from the equator to the footpoint is equal to <math>m=y/k_0</math>. The corresponding footpoint latitude, <math>\phi_1</math>, is calculated by iteration or the inverse meridian distance series as described above. | |||
::<math> | |||
\begin{align} | |||
\mu&=\frac{\pi y}{2 m_p k_0},\\ | |||
\phi_1&=\mu+D_2\sin 2\mu+D_4\sin4\mu+D_6\sin6\mu+D_8\sin8\mu+\cdots,\\ | |||
\end{align} | |||
</math> | |||
Denoting functions evaluated at <math>\phi_1</math> by a subscript '1', the inverse series are: | |||
::<math> | |||
\begin{align} | |||
\lambda(x,y) | |||
&= | |||
\frac{x}{c_1(k_0\nu_1)} | |||
-\frac{x^{3}V_3}{3!c_1(k_0\nu_1)^3} | |||
-\frac{x^{5}V_5}{5!c_1(k_0\nu_1)^5} | |||
-\frac{x^{7}V_7}{7!c_1(k_0\nu_1)^7}, | |||
\\ | |||
\phi(x,y)&=\phi_1 | |||
-\frac{x^2 \beta_1t_1 }{2(k_0\nu_1)^2} | |||
-\frac{x^4 \beta_1t_1 U_4}{4!(k_0\nu_1)^4} | |||
-\frac{x^6 \beta_1t_1 U_6}{6!(k_0\nu_1)^6} | |||
-\frac{x^8 \beta_1t_1 U_8}{8!(k_0\nu_1)^8}. | |||
\end{align} | |||
</math> | |||
<!--=============================================================--> | |||
===Point scale and convergence=== | |||
The point scale <math>k</math> is independent of direction for a conformal transformation. It may be calculated in terms of geographic or projection coordinates. Note that the series for <math>k</math> reduce to <math>k_0</math> when either <math>\lambda=0</math> or <math>x=0</math> . The convergence <math>\gamma</math> may also be calculated (in radians) in terms of geographic or projection coordinates: | |||
::<math> | |||
\begin{align} | |||
k(\lambda,\phi) | |||
&=k_0\left[1 | |||
+\frac{\lambda^2c^2H_2}{2} | |||
+\frac{\lambda^4c^4H_4}{24} | |||
+\frac{\lambda^6c^6H_6}{720}\right],\\ | |||
\gamma(\lambda,\phi) | |||
&=\lambda s | |||
+\frac{\lambda^3c^3t H_3}{3} | |||
+\frac{\lambda^5c^5t H_5}{15} | |||
+\frac{\lambda^7c^7tH_7}{315},\\ | |||
k(x,y)&=k_0\left[1 | |||
+\frac{x^2K_2}{2(k_0\nu_1)^2} | |||
+\frac{x^4K_4}{24(k_0\nu_1)^4} | |||
+\frac{x^6K_6}{720(k_0\nu_1)^6}\right]\\ | |||
\gamma(x,y) &= | |||
\frac{xt_1}{k_0\nu_1} | |||
+ \frac{x^3t_1 K_3}{3(k_0\nu_1)^3} | |||
+ \frac{x^5t_1 K_5}{15(k_0\nu_1)^5} | |||
+ \frac{x^7t_1 K_7}{315(k_0\nu_1)^7}. | |||
\end{align} | |||
</math> | |||
===The coefficients for all series=== | |||
::<math> | |||
\begin{align} | |||
W_3 &= \beta-t^2 \\ | |||
W_5 &= 4\beta^3(1-6t^2)+\beta^2(1+8t^2) -2\beta t^2 +t^4 \\ | |||
W_7 &= 61-479t^2+179t^4-t^6+O(e^2) \\ | |||
W_4 &= 4\beta^2+\beta-t^2 \\ | |||
W_6 &= 8\beta^4(11{-}24t^2)-28\beta^3(1{-}6t^2) +\beta^2(1{-}32t^2) | |||
-2\beta t^2 +t^4 \\ | |||
W_8 &= 1385-3111t^2+543t^4-t^6 +O(e^2)\\ | |||
V_3&= \beta_1+2t_1^2 \\ | |||
V_5&= 4\beta_1^3(1-6t_1^2)-\beta_1^2(9-68t_1^2) | |||
-72\beta_1 t_1^2 -24t_1^4\\ | |||
V_7&= 61+662t_1^2+1320t_1^4+720t_1^6 \\ | |||
U_4&= 4\beta_1^2-9\beta_1(1-t_1^2)-12t_1^2\\ | |||
U_6&= 8\beta_1^4(11-24t_1^2)-12\beta_1^3(21-71t_1^2) | |||
+15\beta_1^2(15-98t_1^2+15t_1^4) \\ | |||
&\qquad\qquad +180\beta_1(5t_1^2-3t_1^4)+360t_1^4\\ | |||
U_8&=-1385-3633t_1^2-4095t_1^4-1575t_1^6\\ | |||
H_2&= \beta\\ | |||
H_4&= 4\beta^3(1-6t^2)+\beta^2(1+24t^2)-4\beta t^2\\ | |||
H_6&=61-148t^2+16t^4\\ | |||
H_3&=2\beta^2-\beta\\ | |||
H_5&=\beta^4(11-24t^2)-\beta^3(11-36t^2) | |||
+\beta^2(2-14t^2)+\beta t^2\\ | |||
H_7&=17-26t^2+2t^4\\ | |||
K_2&=\beta_1 \\ | |||
K_4&=4\beta_1^3(1-6t_1^2)-3\beta_1^2(1-16t_1^2) | |||
-24\beta_1 t_1^2\\ | |||
K_6&=1\\ | |||
K_3&=2\beta_1^2-3\beta_1 -t_1^2\\ | |||
K_5&=\beta_1^4(11-24t_1^2)-3\beta_1^3(8-23t_1^2) | |||
+5\beta_1^2(3-14t_1^2)+30\beta_1 t_1^2+3t_1^4\\ | |||
K_7&=-17-77t_1^2-105t_1^4-45t_1^6 | |||
\end{align} | |||
</math> | |||
<!-- V_4&= 4\beta_1^2-9\beta_1-6t_1^2\\ | |||
V_6&= 8\beta_1^4(11{-}24t_1^2){-}84\beta_1^3(3{-}8t_1^2) | |||
{+}225\beta_1^2(1{-}4t_1^2) | |||
{+}600\beta_1 t_1^2 {+}120t_1^4 \\ | |||
V_8&=-1385-7266t_1^2-10920t_1^4-5040t_1^6\\ | |||
--> | |||
===Accuracy of the series=== | |||
The exact solution of Lee-Thompson,<ref name=lee_exact>Lee, L.P. (1976). ''Conformal Projections Based on Elliptic Functions'' (Supplement No. 1 to Canadian Cartographer, Vol 13.) pp. 1–14, 92–101 and 107–114. Toronto: Department of Geography, York University. A report of unpublished analytic formulae involving incomplete elliptic integrals obtained by E. H. Thompson in 1945. Available from the [http://utpjournals.metapress.com/content/120327/?sortorder=asc&p_o=113|University of Toronto Press].</ref> implemented by Karney (2011),<ref name=karney>C. F. F. Karney (2011), ''[http://dx.doi.org/10.1007/s00190-011-0445-3 Transverse Mercator with an accuracy of a few nanometers],'' J. Geodesy 85(8), 475-485 (2011); preprint of paper and C++ implementation of algorithms are available at [http://geographiclib.sourceforge.net/tm.html tm.html].</ref> is of great value in assessing the accuracy of the truncated Redfearn series. It confirms that the truncation error of the (eighth order) Redfearn series is less than 1 mm out to a longitude difference of 3 degrees, corresponding to a distance of 334 km from the central meridian at the equator but a mere 35 km at the northern limit of an UTM zone. | |||
{{space|5}}The Redfearn series become much worse as the zone widens. Karney discusses Greenland as an instructive example. The long thin landmass is centred on 42W and, at its broadest point, is no more than 750 km from that meridian whilst the span in longitude reaches almost 50 degrees. The Redfearn series attain a maximum error of 1 '''kilometre'''. | |||
==Coordinates, grids, eastings and northings== | |||
The '''projection coordinates''' as defined above are Cartesian coordinates such that the central meridian corresponds to <math>x=0</math> and the equator corresponds <math>y=0</math>. Both ''x'' and ''y'' are defined for all values of <math>\lambda</math> and <math>\phi</math> (although <math>x</math> tends to (±) infinity as <math>\lambda</math> tends to ±90° or ±π/2 radians). The projection does not define a grid: the grid is an independent construct which could be defined arbitrarily. In practice the national implementations, and UTM, use grids aligned with the Cartesian axes of the projection but they are of finite extent, with origins which need not coincide with the origin of the projection coordinates where the central meridian intersects the equator. (Note that in the overlap regions of the UTM zones at neighbouring grids are inclined with respect to each other). | |||
{{space|5}}The '''true grid origin''' is always taken on the central meridian so that grid coordinates will be negative west of the central meridian. To avoid such negative grid coordinates standard practice defines a '''false origin''' to the west (and possibly north or south) of the grid origin: the coordinates relative to the false origin define '''Eastings''' and '''Northings''' which will always be positive. The '''false Easting''', <math>E0</math>, is the distance of the true grid origin east of the false origin. The '''false Northing''', <math>N0</math>, is the distance of the true grid origin north of the false origin. If the true origin of the grid is at latitude <math>\phi_0</math> on the central meridian and the scale factor the central meridian is <math>k_0</math> then these definitions give Eastings and Northings by: | |||
:: <math> | |||
\begin{align} | |||
E&=E0+x(\lambda,\phi),\\[1ex] | |||
N&=N0+y(\lambda,\phi)-k_0 m(\phi_0). | |||
\end{align} | |||
</math> | |||
The inverse series are calculated from the Redfearn series with <math>x</math> replaced by <math>E-E0</math> and the footpoint latitude calculated by inversion of | |||
::<math> | |||
m(\phi_1)=\frac{y}{k_0}=\frac{N-N0}{k_0}+m(\phi_0) | |||
</math> | |||
The terms Eastings and Northings are ''misnomers''. Grid lines of the transverse projection, other than the axes <math>x=0</math> and <math>y=0,</math> do '''not''' run north-south or east-west as defined by parallels and meridians. This is evident from the global projections in [[Transverse Mercator projection]]. Near the central meridian the differences are small but measurable. The difference between the north-south grid lines and the true meridians is the [[Transverse Mercator projection#Convergence|angle of convergence]]. | |||
==Implementations== | |||
The implementations give below are examples of the use of the Redfearn series. The defining documents in various countries differ slightly in notation and, more importantly, in the neglect of some of the small terms. The analysis of small terms depends on the latitude and longitude ranges in the various grids. There are also slight differences in the formulae utilised for meridian distance: one extra term is sometimes added to the formula specified above but such a term is less than 0.1mm. | |||
===OSGB=== | |||
The implementation of the transverse Mercator projection in Great Britain is fully described in the [[OSGB]] document [http://badc.nerc.ac.uk/help/coordinates/OSGB.pdf A guide to coordinate systems in Great Britain], Appendices A.1, A.2 and C.<ref name=osgb/> | |||
::datum: OSGB36 | |||
::ellipsoid: Airy 1830 | |||
:: major axis: 6 377 563.396 | |||
:: minor axis: 6 356 256.909 | |||
::central meridian longitude: 2°W | |||
::central meridian scale factor : 0.9996012717 | |||
:: projection origin: 2°W and 0°N | |||
::true grid origin: 2°W and 49°N | |||
::false easting of true grid origin, E0 (metres): 400,000 | |||
::false northing of true grid origin, N0 (metres): -100,000 | |||
::E = E0 + x = 400000 + x | |||
::N = N0 + y -k0*m(49°)= y - 5527063 | |||
The extent of the grid is 300 km to the east and 400 km to the west of the central meridian and 1300 km north from the ''false'' origin, (OSGB<ref name=osgb/> Section 7.1), but with the exclusion of parts of Northern Ireland, Eire and France. A [[grid reference]] is denoted by the pair (E,N) where E ranges from slightly over zero to 800000m and N ranges from zero to 1300000m. To reduce the number of figures needed to give a grid reference, the grid is divided into 100 km squares, which each have a two-letter code. National Grid positions can be given with this code followed by an easting and a northing both in the range 0 and 99999m. | |||
{{space|5}}The projection formulae differ slightly from the Redfearn formulae presented here. They have been simplified by neglect of most terms of seventh and eighth order in <math>\lambda</math> or <math>x/a</math>: the only exception is seventh order term in the series for <math>\lambda</math> in terms of <math>x/a</math>. This simplification is based on the examination of the Redfearn terms over the ''actual'' extent of the grid. The only other differences are (a) the absorption of the central scale factor into the [[Radius of curvature (mathematics)|radii of curvature]] and [[Meridian arc|meridian distance]], (b) the replacement of the parameter <math>\beta</math> by the parameter <math>\eta</math> (defined [[#Preliminaries I: datum and ellipsoid|above]]). | |||
The OSGB manual<ref name=osgb/> includes a discussion of the [[Helmert transformation]]s which are required to link [[geodetic]] coordinates on Airy 1830 [[Reference ellipsoid|ellipsoid]] and on WGS84. | |||
===UTM=== | |||
The article on the [[Universal Transverse Mercator coordinate system|Universal Transverse Mercator projection]] gives a general survey, but the full specification is defined in U.S. Defence Mapping Agency Technical Manuals TM8358.1<ref name=utm1 /> and TM8358.2.<ref name=utm2 /> This section provides details for '''zone 30''' as another example of the Redfearn formulae (usually termed Thomas formulae in the United States.) | |||
::ellipsoid: International 1924 (a.k.a. Hayford 1909) | |||
:: major axis: 6 378 388.000 | |||
:: minor axis: 6 356 911.946 | |||
::central meridian longitude: 3°W | |||
:: projection origin: 3°W and 0°N | |||
::central meridian scale factor: 0.9996 | |||
::true grid origin: 3°W and 0°N | |||
::false easting of true grid origin, E0: 500,000 | |||
::E = E0 + x = 500000 + x | |||
::northern hemisphere false northing of true grid origin N0: 0 | |||
::northern hemisphere: N = N0 + y = y | |||
::southern hemisphere false northing of true grid origin N0: 10,000,000 | |||
::southern hemisphere: N = N0 + y = 10,000,000 + y | |||
The series adopted for the meridian distance incorporates terms of fifth order in <math>n</math> but the manual states that these are less than 0.03 mm (TM8358.2<ref name=utm2 /> Chapter 2). The projection formulae use, <math>e'</math>, the second eccentrity (defined [[#Preliminaries I: datum and ellipsoid|above]]) instead of <math>n</math>. The grid reference schemes are defined in the article [[Universal Transverse Mercator coordinate system]]. The accuracy claimed for the UTM projections is 10 cm in grid coordinates and 0.001 arc seconds for geodetic coordinates. | |||
===Ireland=== | |||
The transverse Mercator projection in Eire and Northern Ireland (an international implementation spanning one country and part of another) is currently implemented in two ways: | |||
[[Irish grid reference system]] | |||
::datum: Ireland 1965 | |||
::ellipsoid: Airy 1830 modified | |||
:: major axis: 6 377 340.189 | |||
:: minor axis: 6 356 034.447 | |||
::central meridian scale factor: 1.000035 | |||
::true origin: 8°W and 53.5°N | |||
::false easting of true grid origin, E0: 200,000 | |||
::false northing of true grid origin, N0: 250,000 | |||
The Irish grid uses the OSGB projection formulae. | |||
[[Irish Transverse Mercator]] | |||
::datum: Ireland 1965 | |||
::ellipsoid: [[GRS80]] | |||
:: major axis: ? | |||
:: minor axis: ? | |||
::central meridian scale factor: 0.999820 | |||
::true origin: 8°W and 53.5°N | |||
::false easting of true grid origin, E0: 600,000 | |||
::false northing of true grid origin, N0: 750,000 | |||
This is an interesting example of the transition between use of a traditional ellipsoid and a modern global ellipsoid. The adoption of radically different false origins helps to prevent confusion between the two systems. | |||
==See also== | |||
* [[Mercator projection]] | |||
* [[Map projection]] | |||
* [[Scale (map)]] | |||
* [[Universal Transverse Mercator coordinate system]] | |||
==References== | |||
{{Reflist}} | |||
<!--{{cite book | author=Snyder, John P. | title=Map Projections – A Working Manual. U.S. Geological Survey Professional Paper 1395 | publisher =United States Government Printing Office, Washington, D.C | year=1987 | id = }} This paper can be downloaded from [http://pubs.er.usgs.gov/pubs/pp/pp1395 USGS pages]--> | |||
{{DEFAULTSORT:Transverse Mercator: Redfearn Series}} | |||
[[Category:Cartographic projections]] | |||
[[Category:Conformal mapping]] | |||
[[Category:Geocodes]] |
Latest revision as of 14:42, 29 January 2014
The article Transverse Mercator projection restricts itself to general features of the projection. This article describes in detail one of the (two) implementations developed by Louis Krüger in 1912;[1] that expressed as a power series in the longitude difference from the central meridian. These series were recalculated by Lee in 1946,[2] by Redfearn in 1948,[3] and by Thomas in 1952.[4][5] They are often referred to as the Redfearn series, or the Thomas series. This implementation is of great importance since it is widely used in the U.S. State Plane Coordinate System,[5] in national (Britain,[6] Ireland[7] and many others) and also international[8] mapping systems, including the Universal Transverse Mercator coordinate system (UTM).[9][10] They are also incorporated into the Geotrans coordinate converter made available by the United States National Geospatial-Intelligence Agency.[11] When paired with a suitable geodetic datum, the series deliver high accuracy in zones less than a few degrees in east-west extent.
To safe the long run, there are quite a few funding options, one out of which is investing in new launch apartment property. In case you are planning to make investment in Singapore condominiums to safeguard your future, then now we have provide you with sure pros & cons of making funding in the apartment property for your help. Read these pros and cons before you resolve to make the funding at Duo Residences , Thomson View , or different condos.
In the end, New Hyde Park residents obtained their want to preserve the historical spot, and McDonald's had no other option but to revive the property to its former glory. Clara Kirk, who runs two ladies's shelters in Englewood, instructed DNAinfo the price of growing a property into something greater than a backyard or expanded yard might be a problem. The Mayor's officer informed HuffPost that under Chicago's Large Heaps program, candidates would wish to own a property for a minimum of five years before selling. Jade Residences is a new rare freehold residential property launching close to Serangoon MRT and situated at Lew Lian Vale. Benefit from the 50m Lap pool, Kids's Pool in this beautiful improvement. You are not obliged to proceed to buy. Woodlands EC Forestville @ Woodlands
Property investors may be typically categorized into 2 broad classes. These investing with a give attention to rental yield and those investing with a concentrate on capital beneficial properties. Seasoned investors tend to favor accomplished properties as it will probably instantly generate cash move when it comes to rental collections. But when your focus is on capital good points, new launches are inclined to have an important appreciation in worth by the time it HIGH. When a neighborhood developer launched a property in the east in mid-2011, as a substitute of using the standard trade apply minidvd.nl of balloting for the sequence to select a flat, their marketing agent requested interested buyers to line up in front of the gross sales gallery. ROI = (month-to-month rental – mortgage reimbursement – upkeep charge – property tax) x 12 ÷ preliminary funding
Reductions are usually given in the course of the preview, like 5~10% of the list value. Regardless of the xx% discounts (fluctuate from developers to developers), the bottom line is the enticing PSF you may enjoy to buy on the very Preview day. Another benefit is that you've the precedence to decide on your selection unit should no one else choose the same. In a hot property market, many units are snapped up at VIP Preview day, while some initiatives are absolutely sold even earlier than the public come to know about it.
After acquiring a Sale License (subject to government circumstances meant to protect folks shopping for property in Singapore), he might proceed to promote models in his development. Property Launches Listings Map (Singapore & Iskandar Malaysia) - All Properties Listed In SGDevelopersale.com Buyers do NOT, and will NOT, need to pay any agent any price, when buying property in Singapore. PropertyLaunch.sg , is a web site with the aims to supply quality, correct and nicely presented info to all our clients or anybody who are keen on new property launches in Singapore. Marina One Residences By Malaysia Khazanah & Singapore Temasek. SINGAPORE – Singaporeans are streets ahead of every other group of foreigners snapping up property developed by UEM Dawn at Iskandar Malaysia.
RC Suites is an elegant residential and industrial development for the discerning individual who values trendy dwelling in stunning environment. With a daring up to date facade housing 45 cosy flats, providing the perfect residing spaces for younger upwardly cellular professionals.28 RC Suites is near NE8 Farrer Park MRT Station and never removed from several Faculties such as Farrer Park Primary College,Stamford Primary College and Service Hospital Connexion.
After close to 1 / 4-century of doing enterprise at the Denton House in New Hyde Park, although, McDonald's appears to be getting alongside simply high quality with the area people, and residents who originally had an issue with the franchise being there have made the correct changes to get by. As a part of the Inexperienced and Wholesome Chicago Neighborhoods initiative approved by the Chicago Plan Fee Thursday, the specific Giant Tons pilot program will allow qualifying residents and nonprofits to buy metropolis-owned vacant lots for $1 in the Englewood neighborhood on the South Facet. For years, $1 lot programs have cropped up in different cities around the nation. The principles and necessities range, but what they all have in frequent is the next-to-nothing worth.
Preliminaries I: datum and ellipsoid parameters
The series must be used with a geodetic datum which specifies the position, orientation and shape of a Reference ellipsoid. Although the projection formulae depend only on the shape parameters of the reference ellipsoid the full set of datum parameters is necessary to link the projection coordinates to true positions in three dimensional space. The datums and reference ellipsoids associated with particular implementations of the Redfearn formulae are listed below. A comprehensive list of important ellipsoids is given in the article on the Figure of the Earth.
Template:SpaceIn specifying ellipsoids it is normal to give the semi-major axis (equatorial axis), , along with either the inverse flattening, , or the semi-minor axis (polar axis), , or sometimes both. The series presented below use the eccentricity, , in preference to the flattening, . In addition they use the parameters , called the third flattening, and , the second eccentricity. There are only two independent shape parameters and there are many relations between them: in particular
The projection formulae also involve , the radius of curvature of the meridian (at latitude ), and , the radius of curvature in the prime vertical. (The prime vertical is the vertical plane orthogonal to the meridian plane at a point on the ellipsoid). The radii of curvature are defined as follows:
In addition the functions and are defined as:
For compactness it is normal to introduce the following abbreviations:
Preliminaries II: meridian distance
Meridian distance
The article on Meridian arc describes several methods of computing , the meridian distance from the equator to a point at latitude : the expressions given below are those used in the 'actual implementation of the Transverse Mercator projection by the OSGB.[6] The truncation error is less than 0.1mm so the series is certainly accurate to within 1mm, the design tolerance of the OSGB implementation.
where the coefficients are given to order (order ) by
The meridian distance from equator to pole is
The form of the series specified for UTM is a variant of the above exhibiting higher order terms with a truncation error of 0.03mm.
Inverse meridian distance
Neither the OSGB nor the UTM implementations define an inverse series for the meridian distance; instead they use an iterative scheme. For a given meridian distance first set and then iterate using
Template:SpaceThe inversion can be effected by a series, presented here for later reference. For a given meridian distance, , define the rectifying latitude by
The geodetic latitude corresponding to is (Snyder[5] page 17):
An outline of the method
The normal aspect of the Mercator projection of a sphere of radius is described by the equations
where , the isometric latitude, is given by
On the ellipsoid the isometric latitude becomes
By construction, the projection from the geodetic coordinates (,) to the coordinates (,) is conformal. If the coordinates (,) are used to define a point in the complex plane, then any analytic function will define another conformal projection. Kruger's method involves seeking the specific which generates a uniform scale along the central meridian, . He achieved this by investigating a Taylor series approximation with the projection coordinates given by:
where the real part of must be proportional to the meridian distance function . The (complex) coefficients depend on derivatives of which can be reduced to derivatives of with respect to , (not ). The derivatives are straightforward to evaluate in principle but the expressions become very involved at high orders because of the complicated relation between and . Separation of real and imaginary parts gives the series for and and further derivatives give the scale and convergence factors.
The series in detail
This section presents the eighth order series as published by Redfearn[3] (but with and interchanged and the longitude difference from the central meridian denoted by instead of ). Equivalent eighth order series, with different notations, can be found in Snyder[5] (pages 60–64) and at many web sites such as that for the Ordnance Survey of Great Britain.[6]
Template:SpaceThe direct series are developed in terms of the longitude difference from the central meridian, expressed in radians: the inverse series are developed in terms of the ratio . The projection is normally restricted to narrow zones (in longitude) so that both of the expansion parameters are typically less than about 0.1, guaranteeing rapid convergence. For example in each UTM zone these expansion parameters are less than 0.053 and for the British national grid (NGGB) they are less than 0.09. All of the direct series giving , , scale , convergence are functions of both latitude and longitude and the parameters of the ellipsoid: all inverse series giving , , , are functions of both and and the parameters of the ellipsoid.
Direct series
In the following series is the difference of the longitude of an arbitrary point and the longitude of the chosen central meridian: is in radians and is positive east of the central meridian. The W coefficients are functions of listed below. The series for reduces to the scaled meridian distance when .
Inverse series
The inverse series involve a further construct: the footpoint latitude. Given a point on the projection the footpoint is defined as the point on the central meridian with coordinates . Since the scale on the central meridian is the meridian distance from the equator to the footpoint is equal to . The corresponding footpoint latitude, , is calculated by iteration or the inverse meridian distance series as described above.
Denoting functions evaluated at by a subscript '1', the inverse series are:
Point scale and convergence
The point scale is independent of direction for a conformal transformation. It may be calculated in terms of geographic or projection coordinates. Note that the series for reduce to when either or . The convergence may also be calculated (in radians) in terms of geographic or projection coordinates:
The coefficients for all series
Accuracy of the series
The exact solution of Lee-Thompson,[12] implemented by Karney (2011),[13] is of great value in assessing the accuracy of the truncated Redfearn series. It confirms that the truncation error of the (eighth order) Redfearn series is less than 1 mm out to a longitude difference of 3 degrees, corresponding to a distance of 334 km from the central meridian at the equator but a mere 35 km at the northern limit of an UTM zone.
Template:SpaceThe Redfearn series become much worse as the zone widens. Karney discusses Greenland as an instructive example. The long thin landmass is centred on 42W and, at its broadest point, is no more than 750 km from that meridian whilst the span in longitude reaches almost 50 degrees. The Redfearn series attain a maximum error of 1 kilometre.
Coordinates, grids, eastings and northings
The projection coordinates as defined above are Cartesian coordinates such that the central meridian corresponds to and the equator corresponds . Both x and y are defined for all values of and (although tends to (±) infinity as tends to ±90° or ±π/2 radians). The projection does not define a grid: the grid is an independent construct which could be defined arbitrarily. In practice the national implementations, and UTM, use grids aligned with the Cartesian axes of the projection but they are of finite extent, with origins which need not coincide with the origin of the projection coordinates where the central meridian intersects the equator. (Note that in the overlap regions of the UTM zones at neighbouring grids are inclined with respect to each other).
Template:SpaceThe true grid origin is always taken on the central meridian so that grid coordinates will be negative west of the central meridian. To avoid such negative grid coordinates standard practice defines a false origin to the west (and possibly north or south) of the grid origin: the coordinates relative to the false origin define Eastings and Northings which will always be positive. The false Easting, , is the distance of the true grid origin east of the false origin. The false Northing, , is the distance of the true grid origin north of the false origin. If the true origin of the grid is at latitude on the central meridian and the scale factor the central meridian is then these definitions give Eastings and Northings by:
The inverse series are calculated from the Redfearn series with replaced by and the footpoint latitude calculated by inversion of
The terms Eastings and Northings are misnomers. Grid lines of the transverse projection, other than the axes and do not run north-south or east-west as defined by parallels and meridians. This is evident from the global projections in Transverse Mercator projection. Near the central meridian the differences are small but measurable. The difference between the north-south grid lines and the true meridians is the angle of convergence.
Implementations
The implementations give below are examples of the use of the Redfearn series. The defining documents in various countries differ slightly in notation and, more importantly, in the neglect of some of the small terms. The analysis of small terms depends on the latitude and longitude ranges in the various grids. There are also slight differences in the formulae utilised for meridian distance: one extra term is sometimes added to the formula specified above but such a term is less than 0.1mm.
OSGB
The implementation of the transverse Mercator projection in Great Britain is fully described in the OSGB document A guide to coordinate systems in Great Britain, Appendices A.1, A.2 and C.[6]
- datum: OSGB36
- ellipsoid: Airy 1830
- major axis: 6 377 563.396
- minor axis: 6 356 256.909
- central meridian longitude: 2°W
- central meridian scale factor : 0.9996012717
- projection origin: 2°W and 0°N
- true grid origin: 2°W and 49°N
- false easting of true grid origin, E0 (metres): 400,000
- false northing of true grid origin, N0 (metres): -100,000
- E = E0 + x = 400000 + x
- N = N0 + y -k0*m(49°)= y - 5527063
The extent of the grid is 300 km to the east and 400 km to the west of the central meridian and 1300 km north from the false origin, (OSGB[6] Section 7.1), but with the exclusion of parts of Northern Ireland, Eire and France. A grid reference is denoted by the pair (E,N) where E ranges from slightly over zero to 800000m and N ranges from zero to 1300000m. To reduce the number of figures needed to give a grid reference, the grid is divided into 100 km squares, which each have a two-letter code. National Grid positions can be given with this code followed by an easting and a northing both in the range 0 and 99999m.
Template:SpaceThe projection formulae differ slightly from the Redfearn formulae presented here. They have been simplified by neglect of most terms of seventh and eighth order in or : the only exception is seventh order term in the series for in terms of . This simplification is based on the examination of the Redfearn terms over the actual extent of the grid. The only other differences are (a) the absorption of the central scale factor into the radii of curvature and meridian distance, (b) the replacement of the parameter by the parameter (defined above).
The OSGB manual[6] includes a discussion of the Helmert transformations which are required to link geodetic coordinates on Airy 1830 ellipsoid and on WGS84.
UTM
The article on the Universal Transverse Mercator projection gives a general survey, but the full specification is defined in U.S. Defence Mapping Agency Technical Manuals TM8358.1[9] and TM8358.2.[10] This section provides details for zone 30 as another example of the Redfearn formulae (usually termed Thomas formulae in the United States.)
- ellipsoid: International 1924 (a.k.a. Hayford 1909)
- major axis: 6 378 388.000
- minor axis: 6 356 911.946
- central meridian longitude: 3°W
- projection origin: 3°W and 0°N
- central meridian scale factor: 0.9996
- true grid origin: 3°W and 0°N
- false easting of true grid origin, E0: 500,000
- E = E0 + x = 500000 + x
- northern hemisphere false northing of true grid origin N0: 0
- northern hemisphere: N = N0 + y = y
- southern hemisphere false northing of true grid origin N0: 10,000,000
- southern hemisphere: N = N0 + y = 10,000,000 + y
The series adopted for the meridian distance incorporates terms of fifth order in but the manual states that these are less than 0.03 mm (TM8358.2[10] Chapter 2). The projection formulae use, , the second eccentrity (defined above) instead of . The grid reference schemes are defined in the article Universal Transverse Mercator coordinate system. The accuracy claimed for the UTM projections is 10 cm in grid coordinates and 0.001 arc seconds for geodetic coordinates.
Ireland
The transverse Mercator projection in Eire and Northern Ireland (an international implementation spanning one country and part of another) is currently implemented in two ways:
- datum: Ireland 1965
- ellipsoid: Airy 1830 modified
- major axis: 6 377 340.189
- minor axis: 6 356 034.447
- central meridian scale factor: 1.000035
- true origin: 8°W and 53.5°N
- false easting of true grid origin, E0: 200,000
- false northing of true grid origin, N0: 250,000
The Irish grid uses the OSGB projection formulae.
- datum: Ireland 1965
- ellipsoid: GRS80
- major axis: ?
- minor axis: ?
- central meridian scale factor: 0.999820
- true origin: 8°W and 53.5°N
- false easting of true grid origin, E0: 600,000
- false northing of true grid origin, N0: 750,000
This is an interesting example of the transition between use of a traditional ellipsoid and a modern global ellipsoid. The adoption of radically different false origins helps to prevent confusion between the two systems.
See also
References
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
- ↑ Krüger, L. (1912). Konforme Abbildung des Erdellipsoids in der Ebene. Royal Prussian Geodetic Institute, New Series 52.
- ↑ Lee L P, (1946). Survey Review, Volume 8 (Part 58), pp 142–152. The transverse Mercator projection of the spheroid. (Errata and comments in Volume 8 (Part 61), pp 277–278.
- ↑ 3.0 3.1 Redfearn, J C B (1948). Survey Review, Volume 9 (Part 69), pp 318–322, Transverse Mercator formulae.
- ↑ Thomas, Paul D (1952). Conformal Projections in Geodesy and Cartography. Washington: U.S. Coast and Geodetic Survey Special Publication 251.
- ↑ 5.0 5.1 5.2 5.3 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534This paper can be downloaded from USGS pages. It gives full details of most projections, together with interesting introductory sections, but it does not derive any of the projections from first principles. - ↑ 6.0 6.1 6.2 6.3 6.4 6.5 A guide to coordinate systems in Great Britain. This is available as a pdf document at http://www.ordnancesurvey.co.uk/oswebsite/gps/information/coordinatesystemsinfo/guidecontents/
- ↑ See Irish grid reference system and Irish Transverse Mercator
- ↑ Template:Cite webThe EEA recommends the Transverse Mercator for conformal pan-European mapping at scales larger than 1:500,000
- ↑ 9.0 9.1 Defense Mapping Agency Technical Report TM 8358.1: Datums, Ellipsoids, Grids and Grid Reference Systems. URL: http://earth-info.nga.mil/GandG/publications/
- ↑ 10.0 10.1 10.2 J. W. Hager, J.F. Behensky, and B.W. Drew, 1989. Defense Mapping Agency Technical Report TM 8358.2. The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS). URL http://earth-info.nga.mil/GandG/publications/
- ↑ Geotrans, 2010, Geographic translator, version 3.0, URL http://earth-info.nga.mil/GandG/publications/
- ↑ Lee, L.P. (1976). Conformal Projections Based on Elliptic Functions (Supplement No. 1 to Canadian Cartographer, Vol 13.) pp. 1–14, 92–101 and 107–114. Toronto: Department of Geography, York University. A report of unpublished analytic formulae involving incomplete elliptic integrals obtained by E. H. Thompson in 1945. Available from the of Toronto Press.
- ↑ C. F. F. Karney (2011), Transverse Mercator with an accuracy of a few nanometers, J. Geodesy 85(8), 475-485 (2011); preprint of paper and C++ implementation of algorithms are available at tm.html.