<HTML dir=ltr><HEAD><TITLE>Re: [Proj] large difference while projecting utm to utm back and forth</TITLE>
<META content="text/html; charset=unicode" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.18999"></HEAD>
<BODY>
<DIV dir=ltr id=idOWAReplyText99387>
<DIV dir=ltr><FONT color=#000000 size=2 face="Times New Roman">Christoph,</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT> </DIV>
<DIV dir=ltr><FONT size=2>For that much of a zone extension, you need to have additional terms in the classical expansion, say, to the fourteenth derivative for <SPAN style="FONT-FAMILY: 'Times New Roman','serif'; FONT-SIZE: 12pt; mso-ansi-language: EN-US; mso-fareast-font-family: 'Times New Roman'; mso-fareast-language: EN-US; mso-bidi-language: AR-SA">±24°</SPAN>of the Central Meridian; - or - use the closed form with nested iterative loops as now implemented in ESRI software. Most people don't use that because of the extreme convergence angle between Grid North and True North is counter-productive for large-scale mapping. If it's not large-scale (1:50,000 or larger), why use UTM?</FONT></DIV>
<DIV dir=ltr><FONT size=2></FONT><FONT color=#000000 size=2 face="Times New Roman"></FONT> </DIV></DIV>
<DIV dir=ltr id=idSignature6423>
<DIV><FONT color=#000000 size=2 face="Times New Roman"><SPAN style="FONT-SIZE: 10pt">
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt; mso-ansi-language: DE" lang=DE>Clifford J. Mugnier, C.P., C.M.S.<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /><o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Chief of Geodesy,<o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><B><SPAN style="FONT-VARIANT: small-caps; FONT-SIZE: 10pt">Center for GeoInformatics<o:p></o:p></SPAN></B></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Department of Civil Engineering <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Patrick F. Taylor Hall 3223A<o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><B><SPAN style="FONT-SIZE: 10pt">LOUISIANA STATE UNIVERSITY <o:p></o:p></SPAN></B></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Baton Rouge, LA<SPAN style="mso-spacerun: yes"> </SPAN>70803<o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Voice and Facsimile:<SPAN style="mso-spacerun: yes"> </SPAN>(225) 578-8536 [Academic] <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Voice and Facsimile:<SPAN style="mso-spacerun: yes"> </SPAN>(225) 578-4578 [Research] <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Cell: (225) 328-8975 [Academic & Research]<o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Honorary Life Member of the <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Louisiana Society of Professional Surveyors <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Fellow Emeritus of the ASPRS <o:p></o:p></SPAN></DIV>
<DIV style="MARGIN: 0in 0in 0pt" class=MsoNormal><SPAN style="FONT-SIZE: 10pt">Member of the Americas Petroleum Survey Group<o:p></o:p></SPAN></DIV><BR></SPAN></FONT></DIV></DIV>
<DIV dir=ltr><BR>
<HR tabIndex=-1>
<FONT size=2 face=Tahoma><B>From:</B> proj-bounces@lists.maptools.org on behalf of Mikhail Tchernychev<BR><B>Sent:</B> Wed 16-Feb-11 11:13<BR><B>To:</B> PROJ.4 and general Projections Discussions<BR><B>Subject:</B> Re: [Proj] large difference while projecting utm to utm back and forth<BR></FONT><BR></DIV>
<DIV><BR>
<P><FONT size=2> Hello Christoph,<BR><BR>I might be wrong but I'm afraid it is not simply possible. In my<BR>experience with Transverse Mercator<BR>formulas (not with cs2cs) they are ill-posed when you get far from the<BR>central meridian.<BR>In your case you have 5 zones difference, which is 30 degrees. In other<BR>words transformation<BR>is just not defined when you get so far.<BR><BR>It would work if you have points at the border of the zones, say 30 and<BR>31, but<BR>35 is way far.<BR><BR>Best Regards,<BR>Mikhail<BR><BR>On 2/16/2011 5:19 AM, Christoph Dohmen wrote:<BR>> Dear list,<BR>><BR>> I know, I should not do it but I have to :-(<BR>> Given is a point 600000 east and 8260000 north in utm zone 35.<BR>> I have to project this location to utm zone 30.<BR>> Doing some calculations I have to project this value back from<BR>> zone 30 to zone 35.<BR>> But I get differences of about 3km east and more than 200m north.<BR>><BR>> Those are the steps I made using cs2cs:<BR>> cs2cs +proj=utm +zone=35 +datum=WGS84 +to +proj=utm +zone=30 +datum=WGS84<BR>> 600000 8260000<BR>> 1451785.04 8532758.68 0.00<BR>><BR>> and the way back:<BR>> cs2cs +proj=utm +zone=30 +datum=WGS84 +to +proj=utm +zone=35 +datum=WGS84<BR>> 1451785.04 8532758.68<BR>> 597050.49 8260246.21 0.00<BR>><BR>> Can anyone describe a way to achive the startvalues?<BR>><BR>> Best regards<BR>><BR>> Christoph<BR>><BR>> P.S.<BR>> I use proj 4.7.1 and the latest from svn.<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>