
<HTML><HEAD>

<META content="text/html; charset=utf-8" http-equiv=Content-Type>

<META name=GENERATOR content="MSHTML 8.00.6001.23532"></HEAD>

<BODY style="MARGIN: 4px 4px 1px; FONT: 10pt Tahoma">

<DIV>In my proposed comment, I said "<FONT face="Courier New">This is not a straight line on the ellipsoid.</FONT>" which begs the question </DIV>

<DIV>what else did I mean by "straight" on an ellipsoid since I'd already said the sides weren't geodesics?</DIV>

<DIV> </DIV>

<DIV>Since the side of a polygon is defined by a linear relationship between latitude and longitude, it is a rhumb </DIV>

<DIV>line/loxodrome, isn't it?</DIV>

<DIV> </DIV>

<DIV>On reflection, I'd change the second sentence of my proposed comment to:</DIV>

<DIV> </DIV>

<DIV><FONT face="Courier New">The sides of the polygon are rhumb lines and, in general, not geodesics.</FONT></DIV>

<DIV><FONT face="Courier New"></FONT> </DIV>

<DIV>and delete the sentence about the "<FONT face="Courier New">straight line on the ellipsoid</FONT>".</DIV>

<DIV> </DIV>

<DIV>Richard</DIV>

<DIV> </DIV>

<DIV>

<DIV><FONT size=2 face=Arial>

<DIV><FONT color=#ff0000><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>

<P style="MARGIN: 0cm 0cm 0pt" class=MsoNormal><FONT face=Arial><FONT size=2><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>Dr. R. E. Roger</SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US> | Senior Spatial Analyst</SPAN></FONT></FONT></P>

<P style="MARGIN: 0cm 0cm 0pt" class=MsoNormal><FONT face=Arial><FONT size=2><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>NSW Department of Primary Industries </SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US>|</SPAN><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US> NSW Office of Water</SPAN></SPAN></FONT></FONT><SPAN lang=EN-US><FONT color=#003e7e><FONT face=Arial><FONT size=2> <BR><?xml:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" /><st1:Street w:st="on"><st1:address w:st="on">181 <SPAN style="mso-bidi-font-family: Arial">-</SPAN> 183 Anson Street</st1:address></st1:Street></FONT></FONT></FONT></SPAN><FONT face=Arial><FONT size=2><SPAN style="COLOR: #002664" lang=EN-US> | Orange NSW 2800 | <st1:address w:st="on"><st1:Street w:st="on">PO Box</st1:Street> 53</st1:address> | Orange NSW 2800<BR></SPAN><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>T</SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US>: 02 6363 8729 |</SPAN><SPAN lang=EN-US><FONT color=#003e7e> </FONT></SPAN><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>F</SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US>: 02 6361 3839 |</SPAN><SPAN lang=EN-US><FONT color=#003e7e> </FONT></SPAN><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>E</SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US>: <A href="mailto:Richard.Roger@water.nsw.gov.au">Richard.Roger@water.nsw.gov.au</A></SPAN></FONT></FONT><FONT face=Arial><FONT size=2><SPAN lang=EN-US><FONT color=#003e7e> <BR></FONT></SPAN><SPAN class=Redtext><SPAN style="COLOR: #c60c30; mso-bidi-font-family: 'Times New Roman'" lang=EN-US>W</SPAN></SPAN><SPAN style="COLOR: #002664" lang=EN-US>: <A href="http://www.trade.nsw.gov.au/"><SPAN style="COLOR: #002664; TEXT-DECORATION: none; text-underline: none">www.water.nsw.gov.au</SPAN></A></SPAN></FONT></FONT></P>

<P style="MARGIN: 0cm 0cm 0pt" class=MsoNormal><FONT face=Arial></FONT></SPAN></SPAN></P></FONT></DIV></FONT><BR></DIV></DIV>

<DIV><BR>>>> Glynn Clements <glynn@gclements.plus.com> 4/11/2013 12:57 pm >>><BR><BR>Markus Neteler wrote:<BR><BR>> > I hope this is the correct place to post a request for an update to two<BR>> > comments in /lib/gis/area_poly1.c, and replace it with something more<BR>> > informative.<BR>> ><BR>> > My remarks are based on reading /lib/gis/area_poly1.c from<BR>> > grass-7.0.svn_src_snapshot_2013_07_13.<BR>> <BR>> ... I think that his suggestions deserve a feedback, anyone of the<BR>> grass devs knowledgeable with G_begin_ellipsoid_polygon_area()?<BR><BR>It's been a while since I looked at it, but I'm fairly sure that he's<BR>correct. ISTR that the calculation is essentially the area integral of<BR>the scale factor (which is a function of latitude) over a polygon in<BR>the Euclidean plane.<BR><BR>-- <BR>Glynn Clements <glynn@gclements.plus.com><BR></DIV></BODY></HTML>