<HTML dir=ltr><HEAD><TITLE>[Proj] Comments of tmerc, etmerc and ftmerc errors</TITLE>
<META http-equiv=Content-Type content="text/html; charset=unicode">
<META content="MSHTML 6.00.6000.16674" name=GENERATOR></HEAD>
<BODY>
<DIV id=idOWAReplyText43637 dir=ltr>
<DIV dir=ltr><FONT face="Times New Roman" color=#000000 size=2>
<DIV id=idOWAReplyText56060 dir=ltr>
<DIV dir=ltr><FONT face="Times New Roman" color=#000000 size=2>I recommend that error contour maps ( showing the graticule only), be posted to the PROJ4 documentation site - perhaps as "pdf" files.  This would be a superb graphic of "why this and why not that."  Similar things were immortalized in some of the German language geometric geodesy texts of the 1940s - 1950s, but never translated or duplicated into British or American English texts.  </FONT></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr>--------------------------------------------</DIV></DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr>Furthermore, I would like to suggest that a special routine be developed in PROJ4 for "SAFETY FAIRWAYS" that can accommodate the United States Code of Federal Regulations, Volume 33, Chapter 1, Part 166 (33 CFR Ch.1 Pt. 166, Subparts A & B) per the authority of 33 U.S.C. 1223; 49 CFR 1.46.  </DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr>This covers the concept of what the United States government has established as international guarantees that certain shipping lanes within our Exclusive Economic Zone are free of hazards to navigation.  Established by the U.S. Army Corps of Engineers and enforced by the U.S. Coast Guard, this affects the Oil Exploration and Production Industry within U.S. Waters, primarily in the Gulf of Mexico.  The same concept applies internationally to "TRAFFIC SEPARATION SCHEMES" as published by the British Admiralty Office, UK. for the rest of the world.</DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr>When prospecting for hydrocarbons in offshore areas, national governments award exploration and production concessions or ("Lease Blocks in U.S. waters), to Oil Companies according to various national or provincial Grid systems such as UTM, Lambert Conformal Conic (both fully conformal conic as well as the French Army Truncated Cubic), Cassini-Soldner, Ellipsoidal Bonne (really!), Oblique Stereographic, Transverse Mercator, <EM>etc.</EM>  In the Gulf of Mexico, the State Plane Coordinate Systems of Louisiana, Texas, and the UTM are used by the U.S. Bureau of Land Management through the administration of the U.S. Geological Survey, Minerals Management Service for putting out for public bid leases for Offshore Lease Blocks.  These concession leases are defined by a particular grid system, but a senior right exists for any and all Safety Fairways and/or Anchorages as published in 33 CFR, Ch.1, Pt.1.  Interestingly, these Safety Fairways (as well as the international Traffic Separation Schemes), are defined by end points given in latitude and longitude and represent ellipsoidal loxodromes (rhumb lines) in between the end points !!!  The geodetic/cartographic conundrum facing the GIS Specialist is, " where does a particular ellipsoidal loxodrome intersect a specific Grid line defining a Lease Block?"</DIV>
<DIV dir=ltr> </DIV>
<DIV dir=ltr>The conundrum exists because a Lease Block Concession Grid line is only a straight line in the projection space, (such as UTM, Lambert Conformal Conic (both fully conformal conic as well as the French Army Truncated Cubic), Cassini-Soldner, Ellipsoidal Bonne (really!), Oblique Stereographic, Transverse Mercator, <EM>etc.</EM>), but in that <U>same</U> projection space the LOXODROME is <STRONG>curved!</STRONG>  Likewise when the loxodrome is a straight line in isometric normal Mercator space, ... , <STRONG>all other Grid lines are curved!</STRONG></DIV>
<DIV dir=ltr><BR>The solution is iterative.  The loxodrome is easily defined in isometric normal Mercator space by using the two end points as arguments to simply form the equation of the <EM>first </EM>line, (y = mx + b)<SUB><SPAN style="FONT-SIZE: 10pt; COLOR: black">1</SPAN></SUB>, find its approximate intersection with a particular Grid line by taking two points that define an enclosing section of that Grid line, (enclosing the loxodrome), take the inverse of those two Grid points to obtain the latitude and longitude, re-compute them into isometric normal Mercator space, form the equation of that second line, <SPAN style="FONT-SIZE: 10pt; COLOR: black">(y = mx + b)<SUB>2</SUB>, and then solve for the intersection of the two lines.  Then continue to iterate on smaller and smaller sections of the Grid line until one finds the exact intersection point.</SPAN>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">How important is that?  Some years ago, one Oil company to remain unnamed - but still in business nowadays, won a Lease Block bid that was due South of Houston in the Gulf of Mexico.  They filed an application for permit to drill from the Corps of Engineers and were denied a permit to drill in their own Lease Block because of the presence of a protected Anchorage (defined by loxodromes)!!!  They sued the Federal Government, the Minerals Management Service, the Corps of Engineers, the Coast Guard, <EM>etc</EM>., and LOST.  The Federal Judge told them that it was their responsibility to check to see if any protected waters according to 33 CFR, Ch.1, Pt.1. , encroached on any Lease Blocks they were interested in.  It cost that Oil company many, many millions of dollars for that little "geodetic/cartographic conundrum." </SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">Currently, all offshore exploration and production surveying done is with GPS.  The National Datum of the United States is the North American Datum of 1983, identical to the WGS84 Datum at the one-meter level.  However, all Lease Blocks in the Gulf of Mexico <EM>are still defined on</EM> the North American Datum of 1927!  Perplexing, isn't it?</SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">When one looks at current Oil Industry-oriented maps of the entire Gulf of Mexico, there's four different map projections exhibited.  The UTM lease blocks in Federal waters are straight only on the Transverse Mercator, the two State lease block areas are straight only  on Lambert (Texas and Louisiana) zones , and the Safety Fairways are straight only on the normal Mercator.  Funny thing ... everything is forced to be a straight line ... a mathematical impossibility.</SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">Such are some of the geodetic/cartographic tribulations of today's GIS Specialists in the Oil Patch.  Similar conundrums are found elsewhere in the world in association with Traffic Separation Schemes.  It seems to me that this would be a nifty subroutine package for PROJ4 with respect to practical applications.</SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">There's my two cents' worth for the day.</SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black"></SPAN> </P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">Cliff Mugnier</SPAN></P>
<P class=MsoNormal style="MARGIN: 0in 0in 0pt"><SPAN style="FONT-SIZE: 10pt; COLOR: black">LOUISIANA STATE UNIVERSITY</SPAN></P></DIV></FONT></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 Gerald I. Evenden<BR><B>Sent:</B> Thu 12-Jun-08 11:17<BR><B>To:</B> PROJ.4 and general Projections Discussions<BR><B>Subject:</B> [Proj] Comments of tmerc, etmerc and ftmerc errors<BR></FONT><BR></DIV>
<DIV>
<P><FONT size=2>Two testing grids were made on tmerc, etmerc and ftmerc: in a geographic range<BR>of 0 to 6 degrees longitude by 0 to 84 degrees latitude and 0 to 500km<BR>Easting to 9,700km Northing.  The most apparent problem is that for<BR>reasonable demonstration the results should be presented as contour maps of<BR>the errors but inclusion of same is not possible in this email list.<BR><BR>A simple summary can be made: in the lat-lon tests tmerc maintained an<BR>accuracy of <0.7mm (as would be expected occurred at the equator) and<BR>comparison of etmerc and ftmerc was always less than 1micron.  In the<BR>Easting/Northing test, tmerc failed dismally at 21 meters at about 58degrees<BR>longitude/84.78 degree latitude.  Again, ftmerc and etmerc maintained<BR>agreement to within 1micron.<BR><BR>I believe it is quite safe to say that tmerc remains the projection of choice<BR>for UTM applications: adequate accuracy performance and best computational<BR>speed. Note that the lat/lon test was for 12 degree zones and thus tmerc is<BR>would appear suitable for what I understand of the 10 degree zones used for<BR>some Canadian province grids.<BR><BR>As for extended range, I will probably check the agreement of etmerc and<BR>ftmerc over extended range but this is of lesser priority.  Extended<BR>geographic range usage of any conformal projection is a contentious issue as<BR>any resultant grid system has sufficiently large scale errors as to make the<BR>Cartesian usage of the grid very questionable.  Usage of conformal maps for<BR>global mapping is very questionable because of the distortion.  Normal<BR>Mercator persists only because of its usefulness in rhumbline navigation<BR>which is probably fading due to GPS.  Lastly, equal-area mapping is far<BR>preferable for thematic mapping as it at least give a better perspective of<BR>the area extent of regions that is, in most cases, more important that an<BR>often poor preservation of regional shape.<BR>--<BR>The whole religious complexion of the modern world is due<BR>to the absence from Jerusalem of a lunatic asylum.<BR>-- Havelock Ellis (1859-1939)  British psychologist<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>