<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Just project your cells using ellipsoidal cylindric equal area projection and compute the area of the resulting rectangles. This gives you perfect accuracy at the minimum possible computational effort.<br><br>Cheers,<br><div id="AppleMailSignature" dir="ltr">— daan</div><div dir="ltr"><br>On Apr 23, 2019, at 10:00, Thomas Knudsen <<a href="mailto:knudsen.thomas@gmail.com">knudsen.thomas@gmail.com</a>> wrote:<br><br></div><blockquote type="cite"><div dir="ltr"><div dir="ltr"><div dir="ltr">Ken,<br><br>I think I understand your case better now: Your base data are gridded (either model output or EO), and you want to know the “real” area of a cluster of grid cells each nominally being 200 m x 200 m.<br><br>Hence, nominally, your cells are 40000 m2, but really, they differ. So I presume your expected accuracy is +/- 1 cell, i.e. +/-40000 m2.<br><br>Doing the computation using the geodesic planimeter, you would need to trace the outline of each cell cluster (corner coordinate of each outer cell), convert to geographical WGS84, then computing the area of the resulting tiny-sided polygon using the planimeter tool.<br><br>To get a feeling of the precision you may obtain, try comparing the geodesic area of a cell expressed as a polygon defined by its corner coordinates with the properly scaled value for your original approach: Charles Karney states that the accuracy of his algorithm is on the order of 0.1 m2/vertice, i.e. 0.4 m2 for a polygon tracing the outline of a grid cell of nominally 40000 m2 - i.e. a relative accuracy of approximately 1e-5.<br><br>Since the original mapping to the projected grid may be dubious, it may even be a (line-of-code wise) fine shortcut to compute individual geodesic areas for each raster cell, then summing up everything (I suppose the grid parameters are static, and you will only have to do this once).<br><br>/Thomas<br><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den tir. 23. apr. 2019 kl. 11.47 skrev Ken Mankoff <<a href="mailto:mankoff@3m411.com">mankoff@3m411.com</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Thomas,<br>
<br>
Thanks for the reply.<br>
<br>
> But note that the figures you obtain are valid for infinitesimal items<br>
> only. For large lengths/areas you should do the computations with<br>
> geodesics on the ellipsoid<br>
<br>
Can you clarify what you mean by this? Or more specifically, am I doing it incorrectly? I'm exporting a very large grid but at 200 x 200 m resolution to a CSV file, then using proj to calculate the error every 200 m. Is that close enough to infinitesimal? Or define "large lengths" please. I'm then using the errors to correct linear features at the ~200 m cell resolution.<br>
<br>
> using e.g. Charles Karney’s Planimeter tool over at<br>
> <a href="https://geographiclib.sourceforge.io/cgi-bin/Planimeter" rel="noreferrer" target="_blank">https://geographiclib.sourceforge.io/cgi-bin/Planimeter</a> (online<br>
> service, includes link to source etc.)<br>
<br>
I've downloaded and compiled the code. Planimeter takes input in lat/lon, UTM/UPS, or MGRS. I'm working in EPSG:3413. I can provide lon,lat at each cell rather than x,y but I'm not sure if this introduces new errors.<br>
<br>
> Also, a day where you have too much time on your hand, see the<br>
> discussion over at<br>
> <a href="http://osgeo-org.1560.x6.nabble.com/Re-Qgis-user-New-Features-in-Shape-Tools-3-2-0-td5378898.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/Re-Qgis-user-New-Features-in-Shape-Tools-3-2-0-td5378898.html</a><br>
<br>
A nice read - not that I understood all of it - and one reason I use and contribute to OSS/FS.<br>
<br>
-k.<br>
<br>
</blockquote></div>
</div></blockquote><blockquote type="cite"><div dir="ltr"><span>_______________________________________________</span><br><span>PROJ mailing list</span><br><span><a href="mailto:PROJ@lists.osgeo.org">PROJ@lists.osgeo.org</a></span><br><span><a href="https://lists.osgeo.org/mailman/listinfo/proj">https://lists.osgeo.org/mailman/listinfo/proj</a></span><br></div></blockquote></body></html>