|
|
| Line 1: |
Line 1: |
| [[Image:Usgs map traverse mercator.PNG|right|thumb|346px|A transverse Mercator projection]]
| | As they call me Gabrielle. I am a very cashier and I'm ordering pretty good financially. As a definite girl what I simply like is going within order to karaoke but I haven't made a dime you'll. My husband and I chose to reside in Guam but I will own to move in a new year or two. See all that is new on my website here: http://prometeu.net<br><br>Feel free to visit my weblog - [http://prometeu.net clash of clans cheats iphone] |
| The '''transverse Mercator''' [[map projection]] is an adaptation of the standard [[Mercator projection]]. The transverse version is widely used in national and international mapping systems around the world, including the [[Universal Transverse Mercator coordinate system|UTM]]. When paired with a suitable [[Datum (geodesy)|geodetic datum]], the transverse Mercator delivers high accuracy in zones less than a few degrees in east-west extent.
| |
| | |
| ==Standard and transverse aspects==
| |
| The standard (or ''Normal'') Mercator and the transverse Mercator are two different ''[[map projection#Aspects of the projection|aspects]]'' of the same mathematical construction. Because of the common foundation, the transverse Mercator inherits many traits from the normal Mercator:
| |
| | |
| * Both [[map projection|projections]] are [[map projection#Cylindrical|cylindrical]]: for the Normal Mercator, the axis of the cylinder coincides with the polar axis and the line of tangency with the equator. For the transverse Mercator, the axis of the cylinder lies in the equatorial plane, and the line of tangency is any chosen meridian, thereby designated the ''[[Meridian (geography)|central meridian]]''.
| |
| * Both projections may be modified to secant forms, which means the scale has been reduced so that the cylinder slices through the model globe.
| |
| * Both exist in spherical and [[Figure of the earth|ellipsoidal]] versions.
| |
| * Both projections are [[conformal map|conformal]], so that the [[Scale (map)|point scale]] is independent of direction and ''local'' shapes are well preserved;
| |
| * Both projections have constant scale the line of tangency (the equator for the normal Mercator and the central meridian for the transverse).
| |
| | |
| Since the central meridian of the transverse Mercator can be chosen at will, it may be used to construct highly accurate maps (of narrow width) anywhere on the globe. The secant, ellipsoidal form of the transverse Mercator is the most widely applied of all projections for accurate large scale maps.
| |
| | |
| ==Spherical transverse Mercator==
| |
| In constructing a map on any projection, a [[Spherical Earth|sphere]] is normally chosen to model the earth when the extent of the mapped region exceeds a few hundred kilometers in length in both dimensions. For maps of smaller regions, an [[Figure of the Earth|ellipsoidal model]] must be chosen if greater accuracy is required; see next section. The spherical form of the transverse Mercator projection was one of the seven 'new' projections presented, in 1772, by [[Johann Heinrich Lambert]]<ref name=lambert >Lambert, Johann Heinrich. 1772. ''Ammerkungen und Zusatze zurder Land und Himmelscharten Entwerfung''. In [http://www.kuttaka.org/~JHL/L1772a.pdf Beyträge zum Gebrauche der Mathematik und deren Anwendung, part 3], section 6)</ref><ref name=wangerin>Albert Wangerin
| |
| (Editor), 1894. ''Ostwald's Klassiker der exacten Wissenschaften''
| |
| ('''54'''). Published by Wilhelm Engelmann. This is Lambert's paper with additional comments by the editor. Available at the [http://name.umdl.umich.edu/ABR2581.0001.001 University of Michigan Historical Math Library].</ref> (also available in a modern English translation<ref name=tobler>Tobler, Waldo R, ''Notes and Comments on the Composition of Terrestrial and Celestial Maps'', 1972 (University of Michigan Press). Reprinted (2010) by Esri: [http://store.esri.com/esri/showdetl.cfm?SID=2&Product_ID=1284&Category_ID=38]</ref>) Lambert did not name his projections; the name ''transverse Mercator'' dates from the second half of the nineteenth century.<ref name=flattening>{{cite book | author=Snyder, John P. | title=Flattening the Earth: Two Thousand Years of Map Projections | publisher =University of Chicago Press| year=1993|isbn=0-226-76747-7 | page=82}} This is an excellent survey of virtually all known projections from antiquity to 1993.</ref> The principal properties of the transverse projection are here presented in comparison with the properties of the normal projection.
| |
| | |
| ===Normal and transverse spherical projections===
| |
| <!--=============================================================-->
| |
| {| style="text-align:left" style="margin: 1em auto 1em auto"
| |
| |-valign=top
| |
| ! width="1%" |
| |
| ! width="36%"|Normal Mercator
| |
| ! width="3%"|
| |
| ! width="1%" |
| |
| ! width="36%" |Transverse Mercator
| |
| |-valign=top
| |
| |
| |
| | align="center" width="300px" | [[Image:MercNormSph.png|center|thumb|upright|300px|Spherical Normal (equatorial) Mercator (truncated at ''y'' = ±π, corresponding to approximately 85 degrees).]]
| |
| |
| |
| |
| |
| | align="center" width="300px" | [[Image:MercTranSph.png|center|thumb|upright|300px|Spherical transverse Mercator (truncated at ''x'' = ±''π'' in units of Earth radius.)]]
| |
| |-valign=top
| |
| | •
| |
| |The central meridian projects to the straight line ''x'' = 0. Other meridians project to straight lines with ''x'' constant.
| |
| |
| |
| | •
| |
| |The central meridian projects to the straight line ''x'' = 0. Meridians 90 degrees east and west of the central meridian project to lines of constant ''y'' through the projected poles. All other meridians project to complicated curves.
| |
| |-valign=top
| |
| | •
| |
| | The equator projects to the straight line ''y'' = 0 and parallel circles project to straight lines of constant ''y''.
| |
| |
| |
| | •
| |
| |The equator projects to the straight line ''y'' = 0 but all other parallels are complicated closed curves.
| |
| |-valign=top
| |
| | •
| |
| |Projected meridians and parallels intersect at right angles.
| |
| |
| |
| | •
| |
| |Projected meridians and parallels intersect at right angles.
| |
| |-valign=top
| |
| | •
| |
| |The projection is unbounded in the ''y'' direction. The poles lie at infinity.
| |
| |
| |
| | •
| |
| |The projection is unbounded in the ''x'' direction. The points on the equator at ninety degrees from the central meridian are projected to infinity.
| |
| |-valign=top
| |
| | •
| |
| |The projection is conformal. The shapes of small elements are well preserved.
| |
| |
| |
| | •
| |
| |The projection is conformal. The shapes of small elements are well preserved.
| |
| |-valign=top
| |
| | •
| |
| |Distortion increases with ''y''. The projection is not suited for world maps. Distortion is small near the equator and the projection (particularly in its ellipsoidal form) is suitable for accurate mapping of equatorial regions.
| |
| |
| |
| | •
| |
| |Distortion increases with ''x''. The projection is not suited for world maps. Distortion is small near the central meridian and the projection (particularly in its ellipsoidal form) is suitable for accurate mapping of narrow regions.
| |
| |-valign=top
| |
| | •
| |
| |Greenland is almost as large as Africa; the actual area is about one thirteenth that of Africa.
| |
| |
| |
| | •
| |
| |Greenland and Africa are both near to the central meridian; their shapes are good and the ratio of the areas is a good approximation to actual values.
| |
| |-valign=top
| |
| | •
| |
| |The [[Scale (map)|point scale factor]] is independent of direction. It is a function of ''y'' on the projection. (On the sphere it depends on latitude only.) The scale is true on the equator.
| |
| |
| |
| | •
| |
| |The point scale factor is independent of direction. It is a function of ''x'' on the projection. (On the sphere it depends on both latitude and longitude.) The scale is true on the central meridian.
| |
| |-valign=top
| |
| | •
| |
| |The projection is reasonably accurate near the equator. Scale at an angular distance of 5° (in latitude) away from the equator is less than 0.4% greater than scale at the equator, and is about 1.54% greater at an angular distance of 10°.
| |
| |
| |
| | •
| |
| |The projection is reasonably accurate near the central meridian. Scale at an angular distance of 5° (in longitude) away from the central meridian is less than 0.4% greater than scale at the central meridian, and is about 1.54% at an angular distance of 10°.
| |
| |-valign=top
| |
| | •
| |
| |In the secant version the scale is reduced on the equator and it is true on two lines parallel to the projected equator (and corresponding to two parallel circles on the sphere).
| |
| |
| |
| | •
| |
| |In the secant version the scale is reduced on the central meridian and it is true on two lines parallel to the projected central meridian.
| |
| |-valign=top
| |
| | •
| |
| |Convergence (the angle between projected meridians and grid lines with ''x'' constant) is identically zero. Grid north and true north coincide.
| |
| |
| |
| | •
| |
| |Convergence is zero on the equator and non-zero everywhere else. It increases as the poles are approached. Grid north and true north do not coincide.
| |
| |-valign=top
| |
| | •
| |
| |[[Rhumb line]]s (of constant azimuth on the sphere) project to straight lines.
| |
| |
| |
| | •
| |
| |Rhumb lines project to complex curves.
| |
| <!--|-valign=top
| |
| | •
| |
| |left
| |
| |
| |
| | •
| |
| |right
| |
| -->
| |
| |}
| |
| | |
| ==Ellipsoidal transverse Mercator==
| |
| The ellipsoidal form of the transverse Mercator projection was developed by [[Carl Friedrich Gauss]] in 1825<ref name=gauss1825>Gauss, Karl Friedrich, 1825. "Allgemeine Auflösung der Aufgabe: die Theile einer gegebnen Fläche auf einer andern gegebnen Fläche so abzubilden, daß die Abbildung dem Abgebildeten in den kleinsten Theilen ähnlich wird" Preisarbeit der Kopenhagener Akademie 1822. [http://books.google.co.jp/books?id=q2I_AAAAcAAJ&printsec=frontcover&source=gbs_atb#v=onepage&q&f=false Schumacher Astronomische Abhandlungen, Altona, no. 3], p. 5–30. [Reprinted, 1894, Ostwald’s Klassiker der Exakten Wissenschaften, no. 55: Leipzig, Wilhelm Engelmann, p. 57–81, with editing by Albert Wangerin, pp. 97–101. Also in Herausgegeben von der Gesellschaft der Wissenschaften zu Göttingen in Kommission bei Julius Springer in Berlin, 1929, v. 12, pp. 1–9.]</ref> and further analysed by [[Johann Heinrich 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> The projection is known by several names: Gauss Conformal or Gauss-Krüger in Europe; the transverse Mercator in the US; or Gauss-Krüger transverse Mercator generally. The projection is conformal with a constant scale on the central meridian. (There are other conformal generalisations of the transverse Mercator from the sphere to the ellipsoid but only Gauss-Krüger has a constant scale on the central meridian.) Throughout the twentieth century the Gauss-Krüger transverse Mercator was adopted, in one form or another, by many nations (and international bodies);<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> in addition it provides the basis for the [[Universal Transverse Mercator]] series of projections. The Gauss-Krüger projection is now the most widely used projection in accurate large scale mapping.
| |
| | |
| The projection, as developed by Gauss and Krüger, was expressed in terms of low order power series which were assumed to diverge in the east-west direction, exactly as in the spherical version. This was proved to be untrue by British cartographer E.H. Thompson, whose unpublished exact (closed form) version of the projection, reported by L.P. Lee in 1976,<ref name=lee_exact>Lee, L.P. (1976). ''Conformal Projections Based on Elliptic Functions''. Supplement No. 1 to Canadian Cartographer, Vol 13. (Designated as Monograph 16). Toronto: Department of Geography, York University. A report of unpublished analytic formulae involving incomplete elliptic integrals obtained by E.H. Thompson in 1945. The article may be purchased from University of Toronto [http://www.utpjournals.com/]. At the present time (2010) it is necessary to purchase several units in order to obtain the relevant pages: pp 1–14, 92–101 and 107–114.</ref> showed that the ellipsoidal projection is finite (below). This is the most striking difference between the spherical and ellipsoidal versions of the transverse Mercator projection: Gauss-Krüger gives a reasonable projection of the ''whole'' ellipsoid to the plane, although its principal application is to accurate large scale mapping "close" to the central meridian.
| |
| | |
| [[Image:MercTranEll.png|center|thumb|350px|Ellipsoidal Transverse Mercator: a finite projection.]]
| |
| | |
| ===Features===
| |
| :*Near the central meridian (Greenwich in the above example) the projection has low distortion and the shapes of Africa, western Europe, Britain, Greenland, Antarctica compare favourably with a globe.
| |
| :* The central regions of the transverse projections on sphere and ellipsoid are indistinguishable on the small scale projections shown here.
| |
| :*The meridians at 90° east and west of the chosen central meridian project to horizontal lines through the poles. The more distant hemisphere is projected above the north pole and below the south pole.
| |
| :*The equator bisects Africa, crosses South America and then continues onto the complete outer boundary of the projection; the top and bottom edges and the right and left edges must be identified (i.e. they represent the same lines on the globe). (Indonesia is bisected).
| |
| :*Distortion increases towards the right and left boundaries of the projection but it does not increase to infinity. Note the Galapagos Islands where the 90° west meridian meets the equator at bottom left.
| |
| :*The map is conformal. Lines intersecting at any specified angle on the ellipsoid project into lines intersecting at the same angle on the projection. In particular parallels and meridians intersect at 90°.
| |
| :*The point scale factor is independent of direction at any point so that the shape of a '' small'' region is reasonably well preserved. The necessary condition is that the magnitude of scale factor must not vary too much over the region concerned. Note that while South America is distorted greatly the island of Ceylon is small enough to be reasonably shaped although it is far from the central meridian.
| |
| :*The choice of central meridian greatly affects the appearance of the projection. If 90°W is chosen then the whole of the Americas is reasonable. If 145°E is chosen the Far East is good and Australia is oriented with north up.
| |
| | |
| In most applications the [[Gauss–Krüger coordinate system|Gauss–Krüger]] is applied to a narrow strip near the central meridians where the differences between the spherical and ellipsoidal versions are small, but nevertheless important in accurate mapping. Direct series for scale, convergence and distortion are functions of eccentricity and both latitude and longitude on the ellipsoid: inverse series are functions of [[Eccentricity (mathematics)#Ellipses|eccentricity]] and both x and y on the projection. In the secant version the lines of true scale on the projection are no longer parallel to central meridian; they curve slightly. The convergence angle between projected meridians and the x constant grid lines is no longer zero (except on the equator) so that a grid bearing must be corrected to obtain an azimuth from true north. The difference is small, but not negligible, particularly at high latitudes.
| |
| | |
| ===Implementations of the Gauss–Krüger projection===
| |
| In his 1912<ref name="kruger"/> paper, Krüger presented two distinct solutions, distinguished here by the expansion parameter:
| |
| * '''Krüger–n''' (paragraphs 5 to 8). Formulae for the direct projection, giving the coordinates ''x'' and ''y'', are fourth order expansions in terms of the third flattening, ''n'' (the ratio of the difference and sum of the major and minor axes of the ellipsoid). The coefficients are expressed in terms of latitude (φ), longitude (λ), major axis (''a'') and eccentricity (''e''). The inverse formulae for φ and λ are also fourth order expansions in ''n'' but with coefficients expressed in terms of ''x'', ''y'', ''a'' and ''e''. (See [[Transverse Mercator: flattening series]])
| |
| * '''Krüger–λ''' (paragraphs 13 and 14). Formulae giving the projection coordinates ''x'' and ''y'' are expansions (of orders 5 and 4 respectively) in terms of the longitude λ, expressed in radians: the coefficients are expressed in terms of φ, ''a'' and ''e''. The inverse projection for φ and λ are sixth order expansions in terms of the ratio ''x/a'', with coefficients expressed in terms of ''y'', ''a'' and ''e''. (See [[Transverse Mercator: Redfearn series]])
| |
| | |
| The Krüger–λ series were the first to be implemented, possibly because they were much easier to evaluate on the hand calculators of the mid twentieth century.
| |
| *'''Lee–Redfearn–OSGB'''. In 1945 L.P.Lee<ref name=lee_series>Lee L P, (1945). 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> confirmed the λ expansions of Krüger and proposed their adoption by the OSGB<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> but Redfearn (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> pointed out that they were not accurate because of (a) the relatively high latitudes of Great Britain and (b) the great width of the area mapped, over 10 degrees of longitude. Redfearn extended the series to eighth order and examined which terms were necessary to attain an accuracy of 1mm (ground measurement). The [[Transverse Mercator: Krüger–λ series|Redfearn]] series are still the basis of the OSGB map projections.<ref name=osgb />
| |
| *'''Thomas–UTM''' The λ expansions of Krüger were also confirmed by Paul 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> they are readily available in Snyder.<ref name=snyder /> His projection formulae, completely equivalent to those presented by Redfearn, were adopted by the United States Defence Mapping Agency as the basis for the [[Universal Transverse Mercator coordinate system|UTM]].<ref name=utm>{{cite paper |first=J. W. |last=Hager |first2=J. F. |last2=Behensky |first3=B. W. |last3=Drew |year=1989 |title=The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS) |work=Technical Report TM 8358.2, Defense Mapping Agency |url=http://earth-info.nga.mil/GandG/publications/tm8358.2/TM8358_2.pdf }}</ref> They are also incorporated into the Geotrans<ref>Geotrans, 2010, Geographic translator, version 3.0, URL http://earth-info.nga.mil/GandG/geotrans/</ref> coordinate converter made available by the United States National Geospatial-Intelligence Agency [http://earth-info.nga.mil/GandG/geotrans/#zzzb1|NGA].
| |
| *'''Other countries'''. The Redfearn series are the basis for geodetic mapping in many countries: Australia, Germany, Canada, South Africa to name but a few. (A list is given in Appendix A.1 of Stuifbergen 2009.)<ref>N. Stuifbergen, 2009, Wide zone transverse Mercator projection, Technical Report 262, Canadian Hydrographic Service, URL http://www.dfo-mpo.gc.ca/Library/337182.pdf.</ref>
| |
| *Many variants of the Redfearn series have been proposed but only those adopted by national cartographic agencies are of importance. For an example of modifications which do not have this status see [[Transverse Mercator: Bowring series]]). All such modifications have been eclipsed by the power of modern computers and the development of high order ''n''-series outlined below. The precise Redfearn series, although of low order, cannot be disregarded as they are still enshrined in the quasi-legal definitions of OSGB and UTM etc.
| |
| | |
| The Krüger–''n'' series are described on the page [[Transverse Mercator: series in ''n'' (third flattening)]]. They have been implemented (to fourth order in ''n'') by the following nations.
| |
| *'''France'''<ref>http://geodesie.ign.fr/contenu/fichiers/documentation/algorithmes/notice/NTG_76.pdf</ref>
| |
| *'''Finland'''<ref>R. Kuittinen, T. Sarjakoski, M. Ollikainen, M. Poutanen, R. Nuuros, P. Tätilä, J. Peltola, R. Ruotsalainen, and M. Ollikainen, 2006, ETRS89—järjestelmään liittyvät karttaprojektiot, tasokoordinaatistot ja karttalehtijako, Technical Report JHS 154, Finnish Geodetic Institute, Appendix 1, Projektiokaavart, URL
| |
| http://docs.jhs-suositukset.fi/jhs-suositukset/JHS154/JHS154_liite1.pdf.</ref>
| |
| *'''Sweden'''<ref>
| |
| http://www.lantmateriet.se/Global/Kartor%20och%20geografisk%20information/GPS%20och%20m%C3%A4tning/Geodesi/Formelsamling/Gauss_Conformal_Projection.pdf</ref>
| |
| *'''Japan'''<ref>
| |
| http://psgsv2.gsi.go.jp/koukyou/jyunsoku/pdf/H25_junsoku_furoku6.pdf#page=22</ref>
| |
| Higher order versions of the Krüger–''n'' series have been implemented to seventh order by Ensager and Poder<ref name=poder>K. E. Engsager and K. Poder, 2007, [http://icaci.org/documents/ICC_proceedings/ICC2007/documents/doc/THEME%202/oral%201/2.1.2%20A%20HIGHLY%20ACCURATE%20WORLD%20WIDE%20ALGORITHM%20FOR%20THE%20TRANSVE.doc A highly accurate world wide algorithm for the transverse Mercator mapping (almost)], in Proc. XXIII Intl. Cartographic Conf. (ICC2007), Moscow, p. 2.1.2.</ref> and to tenth order by Kawase.<ref name=kawase>Kawase, K. (2011): [http://www.gsi.go.jp/common/000062452.pdf A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection], Bulletin of the [[Geospatial Information Authority of Japan]], '''59''', pp 1–13</ref> Apart from a series expansion for the transformation between latitude and conformal latitude, Karney has implemented the series to thirtieth order.<ref name=karney />
| |
| | |
| ===Exact Gauss-Krüger and accuracy of the truncated series===
| |
| The exact solution of E. H. Thompson, described by L.P. Lee,<ref name=lee_exact /> is summarized on the page [[Transverse Mercator: exact solution]]. It is constructed in terms of elliptic functions (defined in chapters 19 and 22 of the NIST<ref>F. W.J. Olver, D.W. Lozier, R.F. Boisvert, and C.W. Clark, editors,2010, NIST Handbook of Mathematical Functions (Cambridge University Press), available online at URL http://dlmf.nist.gov.</ref> handbook) which can be calculated to arbitrary accuracy using algebraic computing systems such as Maxima.<ref>Maxima, 2009, A computer algebra system, version 5.20.1, URL http://maxima.sf.net.</ref> Such an implementation of the exact solution is described 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>
| |
| | |
| The exact solution is a valuable tool in assessing the accuracy of the truncated ''n'' and λ series. For example, the original 1912 Krüger–''n'' series compares very favourably with the exact values: they differ by less than 0.31 μm within 1000 km of the central meridian and by less than 1 mm out to 6000 km. On the other hand the difference of the Redfearn series used by Geotrans and the exact solution 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. Thus the Krüger–''n'' series are very much better than the Redfearn λ series.
| |
| | |
| 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 while the span in longitude reaches almost 50 degrees. Krüger–''n'' is accurate to within 1mm but the Redfearn version of the Krüger–λ series has a maximum error of 1 '''kilometre'''.
| |
| | |
| Karney's own 8th order (in ''n'') series is accurate to 5 nm within 3900 km of the central meridian.
| |
| | |
| ==Formulae for the spherical transverse Mercator==
| |
| | |
| ===Spherical normal Mercator revisited===
| |
| [[File:Cylindrical Projection basics.svg|thumb|400px|center|The normal aspect of a tangent cylindrical projection of the sphere]]
| |
| The normal cylindrical projections are described in relation to a cylinder tangential at the equator with axis along the polar axis of the sphere. The cylindrical projections are constructed so that all points on a meridian are projected to points with <math>x=a\lambda</math> and <math>y</math> a prescribed function of <math>\phi</math>. For a tangent Normal Mercator projection the (unique) formulae which guarantee conformality are:<ref name=merc/>
| |
| :<math>x = a\lambda\,,\qquad
| |
| y = a\ln \bigg[\tan \bigg(\frac{\pi}{4} + \frac{\phi}{2} \bigg)\bigg]
| |
| = \frac{a}{2}\ln\left[\frac{1+\sin\phi}{1-\sin\phi}\right].
| |
| </math>
| |
| Conformality implies that the [[Scale (map)|point scale]], <math>k</math>, is independent of direction: it is a function of latitude only:
| |
| :<math>k(\phi)=\sec\phi.\,</math>
| |
| For the secant version of the projection there is a factor of <math>k_0</math> on the right hand side of all these equations: this ensures that the scale is equal to <math>k_0</math> on the equator.
| |
| | |
| ===Normal and transverse graticules===
| |
| [[Image:Transverse mercator graticules.svg|thumb|400px|center|Transverse mercator graticules]]
| |
| <!--=============================================================-->
| |
| The figure on the left shows how a transverse cylinder is related to the conventional graticule on the sphere. It is tangential to some arbitrarily chosen meridian and its axis is perpendicular to that of the sphere. The <math>x</math> and <math>y</math> axes defined on the figure are related to the equator and central meridian exactly as they are for the normal projection. In the figure on the right a rotated graticule is related to the transverse cylinder in the same way that the normal cylinder is related to the standard graticule. The 'equator', 'poles' (E and W) and 'meridians' of the rotated graticule are identified with the chosen central meridian, points on the equator 90 degrees east and west of the central meridian, and great circles through those points.
| |
| | |
| [[File:Transverse mercator geometry.svg|thumb|170px|right|Transverse mercator geometry]]
| |
| {{space|5}}The position of an arbitrary point <math>(\phi,\lambda)</math> on the standard graticule can also be identified in terms of angles on the rotated graticule: <math>\phi'</math> (angle M'CP) is an effective latitude and <math>-\lambda'</math> (angle M'CO) becomes an effective longitude. (The minus sign is necessary so that <math>(\phi',\lambda')</math> are related to the rotated graticule in the same way that <math>(\phi,\lambda)</math> are related to the standard graticule). The Cartesian <math>(x',y')</math> axes are related to the rotated graticule in the same way that the axes <math>(x,y)</math> axes are related to the standard graticule.
| |
| | |
| {{space|5}}The tangent transverse Mercator projection defines the coordinates <math>(x',y')</math> in terms of <math>-\lambda'</math> and <math>\phi'</math> by the transformation formulae of the tangent Normal Mercator projection:
| |
| ::<math>x' = -a\lambda'\,\qquad
| |
| y' = \frac{a}{2}
| |
| \ln\left[\frac{1+\sin\phi'}{1-\sin\phi'}\right].
| |
| </math>
| |
| This transformation projects the central meridian to a straight line of finite length and at the same time projects the great circles through E and W (which include the equator) to infinite straight lines perpendicular to the central meridian. The true parallels and meridians (other than equator and central meridian) have no simple relation to the rotated graticule and they project to complicated curves.
| |
| | |
| ===The relation between the graticules===
| |
| The angles of the two graticules are related by using [[spherical trigonometry]] on the spherical triangle NM'P defined by the true meridian through the origin, OM'N, the true meridian through an arbitrary point, MPN, and the great circle WM'PE. The results are:<ref name=merc/>
| |
| ::<math>
| |
| \begin{align}
| |
| \sin\phi'&=\sin\lambda\cos\phi,\\
| |
| \tan\lambda'&=\sec\lambda\tan\phi.
| |
| \end{align}
| |
| </math>
| |
| | |
| ===Direct transformation formulae===
| |
| The direct formulae giving the Cartesian coordinates <math>(x,y)</math> follow immediately from the above. Setting <math>x=y'</math> and <math>y=-x'</math> (and restoring factors of <math>k_0</math> to accommodate secant versions)
| |
| ::<math>
| |
| \begin{align}
| |
| x(\lambda,\phi)&= \frac{1}{2}k_0a
| |
| \ln\left[
| |
| \frac{1+\sin\lambda\cos\phi}
| |
| {1-\sin\lambda\cos\phi}\right],\\
| |
| y(\lambda,\phi)&= k_0 a\arctan\left[\sec\lambda\tan\phi\right],
| |
| \end{align}
| |
| </math>
| |
| The above expressions are given in Lambert<ref name=lambert/> and also (without derivations) 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> Maling<ref name=maling>{{cite book | author=Maling, Derek Hylton | title=Coordinate Systems and Map Projections | publisher =Pergamon Press| year=1992|isbn=0-08-037233-3 |edition=second}}.</ref> and Osborne<ref name=merc>[http://www.mercator99.webspace.virginmedia.com The Mercator Projections] Detailed derivations of all formulae quoted in this article</ref> (with full details).
| |
| | |
| ===Inverse transformation formulae===
| |
| Inverting the above equations gives
| |
| ::<math>
| |
| \begin{align}
| |
| \lambda(x,y)&
| |
| = \arctan\bigg[ \sinh\frac{x}{k_0a}
| |
| \sec\frac{y}{k_0a} \bigg],
| |
| \\[1ex]
| |
| \phi(x,y)&= \arcsin\bigg[ \mbox{sech}\;\frac{x}{k_0a}
| |
| \sin\frac{y}{k_0a} \bigg].
| |
| \end{align}
| |
| </math>
| |
| | |
| ===Point scale===
| |
| In terms of the coordinates with respect to the rotated graticule the [[Scale (map)|point scale]] factor is given by <math>k=\sec\phi'</math>: this may be expressed either in terms of the geographical coordinates or in terms of the projection coordinates:
| |
| ::<math>
| |
| \begin{align}
| |
| k(\lambda,\phi)&=\frac{k_0}{(1-\sin^2\lambda\cos^2\phi)^{1/2}},\\
| |
| k(x,y)&=k_0\cosh\bigg(\frac{x}{k_0a}\bigg).
| |
| \end{align}
| |
| </math>
| |
| The second expression shows that the scale factor is simply a function of the distance from the central meridian of the projection. A typical value of the scale factor is <math>k_0=0.9996</math> so that <math>k=1</math> when <math>x</math> is approximately 180 km. When <math>x</math> is approximately 255 km and <math>k_0=1.0004</math>: the scale factor is within 0.04% of unity over a strip of about 510 km wide.
| |
| | |
| {{anchor|convergence}}
| |
| | |
| ===Convergence===
| |
| [[Image:Transverse mercator convergence.svg|right|thumb|upright|150px|The angle of convergence]]
| |
| | |
| The convergence angle <math>\gamma</math> at a point on the projection is defined by the angle measured ''from'' the projected meridian, which defines true north, ''to'' a grid line of constant x, defining grid north. Therefore <math>\gamma</math> is positive in the quadrant north of the equator and east of the central meridian and also in the quadrant south of the equator and west of the central meridian. The convergence must be added to a grid bearing to obtain a bearing from true north. For the secant transverse Mercator the convergence may be expressed<ref name=merc/> either in terms of the geographical coordinates or in terms of the projection coordinates:
| |
| ::<math>
| |
| \begin{align}
| |
| \gamma(\lambda,\phi)&=\arctan(\tan\lambda\sin\phi),\\
| |
| \gamma(x,y)&=\arctan\bigg(\tanh\frac{x}{k_0a}\tan\frac{y}{k_0a}\bigg).
| |
| \end{align}
| |
| </math>
| |
| | |
| ==Formulae for the ellipsoidal transverse Mercator==
| |
| Details of actual implementations
| |
| | |
| * Gauss-Kruger series in longitude: [[Transverse Mercator: Redfearn series]]
| |
| * Gauss-Kruger series in ''n'' (third flattening): [[Transverse Mercator: flattening series]] {{See also|Universal Transverse Mercator coordinate system#Simplified formulas}}
| |
| * Exact (closed form) transverse Mercator projection: [[Transverse Mercator: exact solution]]
| |
| * Fourth order Redfearn series by concise formulae (example): [[Transverse Mercator: Bowring series]]
| |
| | |
| ==See also==
| |
| {{Portal|Atlas}}
| |
| * [[List of map projections]]
| |
| * [[Mercator projection]]
| |
| * [[Map projection]]
| |
| * [[Scale (map)]]
| |
| * [[Universal Transverse Mercator coordinate system]]
| |
| | |
| ==References==
| |
| {{Reflist|30em}}
| |
| <!--{{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]-->
| |
| | |
| ==External links==
| |
| * The projections used to illustrate this article were prepared using Geocart which is available from http://www.mapthematics.com
| |
| | |
| {{Map Projections}}
| |
| | |
| {{DEFAULTSORT:Transverse Mercator Projection}}
| |
| [[Category:Cartographic projections]]
| |
| [[Category:Conformal mapping]]
| |
| [[Category:Geocodes]]
| |
| | |
| [[de:UTM-Koordinatensystem]]
| |
| [[es:Coordenadas UTM]]
| |
| [[fr:Transverse Universelle de Mercator]]
| |
| [[no:UTM-koordinater]]
| |
| [[nn:UTM-koordinatar]]
| |
| [[pt:Universal Transversa de Mercator]]
| |
| [[tr:UTM Projeksiyonu]]
| |