<div dir="ltr"><div><div><div><div>Hello,<br><br></div>Thank you for your investigations. <br><br></div>Is there anything else I can do about this? I trust my QGIS measurement because I know the size of the object I have drawn around. However I would really like to use GRASS to automate my workflow... <br><br></div>Thanks<br><br></div>James<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 19 May 2017 at 13:41, Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 19/05/17 11:25, James Duffy wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Just to clear a few things up.<br>
<br>
The shapefile contains only a single polygon (I created it as a simple<br>
test to work out this discrepancy).<br>
<br>
OTF is turned off in QGIS, I am working in EPSG 4326.<br>
<br>
The 0.127m2 comes from using the Measure tool (OTF off) in QGIS. And it<br>
is set to display m2.<br>
<br>
Using the $area calculator in the field calculator in QGIS does indeed<br>
show something that is more likely degrees squared.<br>
</blockquote>
<br></span>
Below some results from quick tests.<br>
<br>
Summary: personally I would trust GRASS GIS a bit more, but I'm biaised...<br>
<br>
***************************<br>
Test 1:<br>
<br>
GRASS trunk in NC demo location (EPSG 3358)<br>
<br>
g.region n=200100 s=200000 w=500000 e=500100 -ap<br>
v.in.region out=nc_testpolygon<br>
v.to.db -p nc_testpolygon op=area --q<br>
1|10000<br>
<br>
GRASS trunk in EPSG 4326 location<br>
v.proj location=nc_spm_08 mapset=user1 input=nc_testpolygon<br>
v.to.db -p nc_testpolygon op=area --q<br>
1|10002.503642603<br>
v.out.ogr nc_testpolygon out=nc_testpolygon.shp (i.e. shapefile is in degress, EPSG 4326)<br>
<br>
QGIS (using the exported shapefile)<br>
<br>
No OTF<br>
measuring tool: 12.324,737 m²<br>
$area in calculator: 12320.2346026<br>
<br>
OTF to EPSG 4326, with datum WGS84<br>
measuring tool: 10.072,355 m²<br>
$area in calculator: 10046.3679510<br>
<br>
<br>
******************************<wbr>**<br>
Test 2:<br>
<br>
GRASS trunk, in ESPG 4326 location<br>
<br>
g.region n=74 s=73 w=0 e=1 res=1 -p<br>
v.in.region region_1degree_73_74N<br>
v.to.db -p region_1degree_73_74N op=area --q<br>
1|3539348689.40454<br>
<br>
v.out.ogr region_1degree_73_74N out=region_1degree_73_74N.shp<br>
<br>
<br>
GRASS trunk in EPSG 3575 location<br>
v.proj location=LL_WGS84 mapset=mlennert input=region_1degree_73_74N<br>
v.to.db -p region_1degree_73_74N op=area --q<br>
1|3539337483.97153<br>
<br>
<br>
QGIS (using the exported shapefile)<br>
no OTF projection<br>
<br>
measuring tool 1,002 m² (sic)<br>
export geometry column: 1<br>
$area in calculator <a href="tel:12392029030" value="+12392029030" target="_blank">12392029030</a>.5000000<br>
<br>
with OTF projection, datum for measurements WGS84<br>
<br>
measuring tool 3.503.322.774,337 m²<br>
export geometry column: 1<br>
$area in calculator: 3499702798.3400998<br>
<br>
reproject to EPSG 3575<br>
with OTF, datum WGS84<br>
<br>
export geometry column 3539169026.18628<br>
export geoemtry column (column named otf) 3499702798.00000<br>
$area in calculator: 3499702798.6579294<br>
<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span style="color:rgb(39,78,19)"><span style="font-family:tahoma,sans-serif"><font size="4"><b>James Duffy</b></font><br>PhD Researcher<br>Environment and Sustainability Institute<br>Penryn Campus<br>University of Exeter<br>Penryn<br>Cornwall<br>TR10 9FE</span></span><br></div></div></div></div></div></div>
</div>