<HTML dir=ltr><HEAD><TITLE>[Proj] Transverse Mercator algorithm with good accuracy/speedtrade-off?</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.6000.16640" name=GENERATOR></HEAD>
<BODY>
<DIV id=idOWAReplyText94165 dir=ltr>
<DIV dir=ltr><FONT face="Times New Roman" color=#000000 size=2>Dozier came up with an iterative forward and inverse Transverse Mercator that's good to almost 90° from the Central Meridian with sub-millimeter round-trip precision. It was written up in an obscure NOAA pub some 20+ years ago in pre-ANSI "C" language. ESRI has implemented it with substantial code and algebraic corrections, and I understand that it works very well.</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Clifford J. Mugnier, C.P., C.M.S.</FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>LSU Student ASCE Chapter Faculty Advisor</FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>and<BR>Past National Director (2006-2008),<BR>Photogrammetric Applications Division<BR>American Society for Photogrammetry and Remote Sensing </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>and </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Chief of Geodesy,<BR>CENTER FOR GEOINFORMATICS </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Department of Civil Engineering </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Patrick F. Taylor Hall 3223A<BR>LOUISIANA STATE UNIVERSITY </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Baton Rouge, LA 70803<BR>Voice and Facsimile: (225) 578-8536 [Academic] </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Voice and Facsimile: (225) 578-4474 [Research] </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Honorary Life Member of the </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Louisiana Society of Professional Surveyors </FONT></DIV>
<DIV dir=ltr><FONT face="Times New Roman" size=2>Member Emeritus of the ASPRS </FONT></DIV>
<DIV dir=ltr><FONT size=2><FONT face="Times New Roman">Member of the Americas Petroleum Survey Group<BR>======================================================<BR></FONT><A href="https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.lsu.edu/exchweb/bin/redir.asp?URL=https://email.ls" target=_blank><FONT face="Times New Roman">http://www.asprs.org/resources/GRIDS/</FONT></A><BR><FONT face="Times New Roman">======================================================</FONT></FONT></DIV></DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT face=Tahoma size=2><B>From:</B> proj-bounces@lists.maptools.org on behalf of Mikael Rittri<BR><B>Sent:</B> Thu 22-May-08 04:04<BR><B>To:</B> proj@lists.maptools.org<BR><B>Subject:</B> [Proj] Transverse Mercator algorithm with good accuracy/speedtrade-off? <BR></FONT><BR></DIV>
<DIV>
<P><FONT size=2>Hello,<BR>I am thinking about extending the accuracy of our Transverse Mercator. But is it worth the trouble? <BR><BR>Summary:<BR>-------<BR>tmerc is fast, but inaccurate far from the central line.<BR>I have implemented a quite exact algorithm, but it was much slower than tmerc. Is there a happy compromise?<BR><BR>Long-winded motivation:<BR>----------------------<BR>The tmerc of Proj4 is efficient. But it has millimeter accuracy only out to 6° of longitude from the central meridian [Oscar van Vlijmen, <http://lists.maptools.org/pipermail/proj/2006-June/002316.html]. <BR><BR>Sometimes, a Transverse Mercator is used farther away: Sweden uses UTM zone 33 (lon_0 = 15°E) for the whole country, and the eastmost point Kataja is at 65°42'N, 24°10'E. I think Norway does the same, with the same projection, and its eastmost point Hornøya is at 70°22'N, 31°10'E. I have found the following errors for tmerc:<BR><BR> errors forward inverse<BR> Kataja 0.25 mm 16.5 mm<BR> Hornøya 7.51 mm 2715 mm<BR><BR>The local scale factor is only 1.0018 at Kataja and 1.0040 at Hornøya, so one can argue that ellipsoid formulas should indeed be used here, rather than spherical formulas. Snyder writes that ellipsoid formulas can be motivated to some 10° to 15° of arc distance from the central meridian [Map projections: A working manual, page 48].<BR>Although Hornøya is 16°10' of longitude away from lon_0, it is only 599 km or 5°24' of arc distance away, because of the high latitude: cos(70°22')=0.336.<BR> So, I am slightly worried about these errors (even though there are other Swedish and Norwegian projections intended for local use at detailed scales). <BR> A more trivial reason for higher accuracy is that I would like to make our geographic toolkit more foolproof.<BR>The mismatch between forward and inverse can cause problems.<BR><BR>Oscar van Vlijmen wrote<BR>[<A href="http://lists.maptools.org/pipermail/proj/2006-June/002296.html">http://lists.maptools.org/pipermail/proj/2006-June/002296.html</A>]:<BR>> At least three geodetic services use routines approximating the exact<BR>> TM better than tmerc, DMA/NIMA/NGA and the like:<BR>> 1) French IGN<BR>> <<A href="http://www.ign.fr/telechargement/MPro/geodesie/CIRCE/NTG_76.pdf">http://www.ign.fr/telechargement/MPro/geodesie/CIRCE/NTG_76.pdf</A>><BR>> 2) Swedish Lantmäteriet<BR>> <<A href="http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/geodesi/Formelsamling/Gauss_Conformal_Projection.pdf">http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/geodesi/Formelsamling/Gauss_Conformal_Projection.pdf</A>> <BR>> 3) Finnish JHS <BR>> <<A href="http://www.jhs-suositukset.fi/intermin/hankkeet/jhs/home.nsf/files/JHS154/$file/JHS154.pdf">http://www.jhs-suositukset.fi/intermin/hankkeet/jhs/home.nsf/files/JHS154/$file/JHS154.pdf</A>><BR>> Each follow a slightly different route, but the<BR>> differences in the results are small.<BR><BR>(The Finnish report has been moved, and is now on<BR> <A href="http://docs.jhs-suositukset.fi/jhs-suositukset/JHS154/JHS154.pdf">http://docs.jhs-suositukset.fi/jhs-suositukset/JHS154/JHS154.pdf</A>)<BR><BR>I have implemented the Swedish proposal[*], which was not too difficult, and I got it to agree with Oscar van Vlijmen's test point. But my implementation is about 4.5 times slower than tmerc in the forward direction, and about 15 times slower in the inverse direction. I doubt that I can improve the speed much.<BR><BR>So, finally, I wonder if there is some happy compromise:<BR>an algorithm that is somewhat more accurate than tmerc, but not too much slower?<BR><BR>Best regards,<BR> Mikael<BR><BR>PS: I noted that Clifford Mugnier has implemented a more accurate version of the tmerc/Snyder formulas, giving millimeter accuracy 24 degrees away, based on US Army Tech. Manual TM 5-241-10, "Universal Transverse Mercator Grid: Extension of zone to zone transformation tables", 1962, 642 pages.<BR>[<A href="http://lists.maptools.org/pipermail/proj/2003-January/thread.html">http://lists.maptools.org/pipermail/proj/2003-January/thread.html</A>]<BR>But as I understood it, this version is cumbersome to implement (and the report is quite hard to find).<BR><BR>[*] footnote: the Swedish note gives no literature references, but I have learned from Lars E. Engberg that the formulas come from<BR> L. Krüger's book:<BR> Konforme Abbildung des Erdellipsoids in der Ebene,<BR> Veröffentlichung des Könglich Preuzischen Geodätischen Institutes,<BR> Neue Folge 52, Leipzig 1912.<BR>The formulas also appear in<BR> Kurt Bretterbauer:<BR> Gebrauchsformeln für die UTM-Projektion nach Krüger,<BR> Österr. Zeitschr. für Vermessung und Geoinformation, 3/2003;<BR> abstract in German and English at <A href="http://www.ovg.at/index.php?id=222">http://www.ovg.at/index.php?id=222</A>.<BR><BR>--<BR>Mikael Rittri<BR>Carmenta AB<BR>Box 11354<BR>SE-404 28 Göteborg<BR>Visitors: Sankt Eriksgatan 5<BR>SWEDEN<BR>Tel: +46-31-775 57 37<BR>Mob: +46-703-60 34 07<BR>mikael.rittri@carmenta.com<BR>www.carmenta.com<BR><BR>_______________________________________________<BR>Proj mailing list<BR>Proj@lists.maptools.org<BR><A href="http://lists.maptools.org/mailman/listinfo/proj">http://lists.maptools.org/mailman/listinfo/proj</A><BR></FONT></P></DIV></BODY></HTML>