<div dir="ltr"><div dir="ltr">Ken<br><br>The scale is the quotient between the length in the map and the length on the ellipsoid, so essentially you’re right. But note that the figures you obtain are valid for infinitesimal items only.<br><br>For large lengths/areas you should do the computations with geodesics on the ellipsoid, using e.g. Charles Karney’s Planimeter tool over at <a href="https://geographiclib.sourceforge.io/cgi-bin/Planimeter">https://geographiclib.sourceforge.io/cgi-bin/Planimeter</a> (online service, includes link to source etc.)<br><br>Also, a day where you have too much time on your hand, see the discussion over at <a href="http://osgeo-org.1560.x6.nabble.com/Re-Qgis-user-New-Features-in-Shape-Tools-3-2-0-td5378898.html">http://osgeo-org.1560.x6.nabble.com/Re-Qgis-user-New-Features-in-Shape-Tools-3-2-0-td5378898.html</a><br><br>/Thomas<br></div><div dir="ltr"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den tir. 23. apr. 2019 kl. 06.14 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">Dear Proj list,<br>
<br>
Please excuse duplicates - Even sent this with a different subject "Fwd: proj list issue" and I have not gotten any replies. I'm resending with an improved subject in case people skipped the email with that subject.<br>
<br>
I'm following up on a helpful thread from this list from 2017: <a href="https://lists.osgeo.org/pipermail/proj/2017-November/007496.html" rel="noreferrer" target="_blank">https://lists.osgeo.org/pipermail/proj/2017-November/007496.html</a><br>
<br>
I'm calculating map projection errors, and want to check my interpretation of the results with someone who knows proj better. I'm using EPSG:3413 projection and the following proj string:<br>
<br>
+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1<br>
<br>
To estimate the error in Southern Greenland, I do:<br>
<br>
echo -40 60 | proj -VS ${PROJSTR}<br>
<br>
and see<br>
<br>
#Stereographic<br>
# Azi, Sph&Ell<br>
# lat_ts=<br>
# +proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +no_defs<br>
# +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1<br>
#Final Earth figure: ellipsoid<br>
# Major axis (a): 6378137.000<br>
# 1/flattening: 298.257224<br>
# squared eccentricity: 0.006694379990<br>
Longitude: 40dW [ -40 ]<br>
Latitude: 60dN [ 60 ]<br>
Easting (x): 289632.50<br>
Northing (y): -3310514.64<br>
Meridian scale (h) : 1.03942808 ( 3.943 % error )<br>
Parallel scale (k) : 1.03942808 ( 3.943 % error )<br>
Areal scale (s): 1.08041073 ( 8.041 % error )<br>
Angular distortion (w): 0.000<br>
Meridian/Parallel angle: 90.00000<br>
Convergence : 5d [ 5.00000000 ]<br>
Max-min (Tissot axis a-b) scale error: 1.03943 1.03943<br>
<br>
The key values for Northern Greenland are:<br>
<br>
echo -40 82 | proj -VS ${PROJSTR}<br>
...<br>
...<br>
Meridian scale (h) : 0.97460025 ( -2.54 % error )<br>
Parallel scale (k) : 0.97460025 ( -2.54 % error )<br>
Areal scale (s): 0.94984565 ( -5.015 % error )<br>
<br>
My question is how to interpret the Areal scale error<br>
<br>
I think that in this projection N. Greenland is over-represented by 8% area and any areal feature there should be divided by 1.0804 or linear features divided by 1.0394 to bring them to true scale. Similarly, in S. Greenland I should adjust a areal feature by dividing by 0.94984565. Is this correct?<br>
<br>
Essentially, does "-2.54 % error" mean that the EPSG:3413 projection is smaller than it should be, or that the "true" spheroid is smaller than EPSG:3413 here?<br>
<br>
Thanks,<br>
<br>
-k.<br>
<br>
_______________________________________________<br>
PROJ mailing list<br>
<a href="mailto:PROJ@lists.osgeo.org" target="_blank">PROJ@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/proj</a><br>
</blockquote></div>