<div dir="ltr"><div dir="ltr"><div dir="ltr"><br><br>On Tue, Sep 25, 2018 at 7:38 PM Helmut Kudrnovsky <<a href="mailto:hellik@web.de">hellik@web.de</a>> wrote:<br>><br>> fyi see<br>> <a href="https://lists.osgeo.org/pipermail/qgis-developer/2018-September/054644.html">https://lists.osgeo.org/pipermail/qgis-developer/2018-September/054644.html</a><br>><br>> with GRASS mentioned<br>> ------------------<br>> Kristian Evers:<br>><br>> Right, here are the calculated areas as returned by a number of different<br>> GIS applications and the planimeter app of GeographicLib for reference:<br>><br>> Caris LOTS: 14.737 km^2<br>> ArcMap:     14.727,446 km^2<br>> MapInfo:    14.727,352 km^2<br>> GeoMedia:   14.726,443 km^2<br>> Planimeter: 14.722,522 km^2<br>> QGIS 3.2:   14.652,181 km^2<br>> QGIS 2.8:   14.652,181 km^2<br><div><br></div><div><div>adding to the confusion:</div><div><br></div><div>I used the geographiclib API as included in PROJ 5.2.0 following the example for geod_polygonarea() in geodesic.h and get</div><div>geographiclib: 14,737.935 km^2<br></div><div>quite different from</div><div>Planimeter: 14,722.522 km^2</div><div><br></div><div>GRASS native gives 14,718.097679</div><div>as also reported by Helmut and Stefan</div><div><br></div><div>Markus M<br></div></div><div>></div>> The polygon that I have used to get the numbers above can be found here:<br>> <a href="https://gist.github.com/kbevers/207b5bcb9be20e7554abe5f56742ec2c">https://gist.github.com/kbevers/207b5bcb9be20e7554abe5f56742ec2c</a><br>><br>> I am quite confident that GeographicLib delivers the most accurate result<br>> (if you have doubts, this reference [0] should convince you). As can be seen<br>> from the table above all but QGIS come fairly close. I expect some variation<br>> in the results as these are numerical approximations, although I think QGIS<br>> is too far of the mark. My suspicion is that the geodesic algorithm used by<br>> QGIS (and apparently GRASS) is to blame here.<br>><br>> /Kristian<br>><br>> [0] <a href="https://arxiv.org/pdf/1102.1215.pdf">https://arxiv.org/pdf/1102.1215.pdf</a><br>> -----------------<br>><br>><br>><br>><br>><br>> -----<br>> best regards<br>> Helmut<br>> --<br>> Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Dev-f3991897.html">http://osgeo-org.1560.x6.nabble.com/Grass-Dev-f3991897.html</a><br>> _______________________________________________<br>> grass-dev mailing list<br>> <a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-dev">https://lists.osgeo.org/mailman/listinfo/grass-dev</a></div></div></div>