<div dir="ltr"><div><div><br><br>On Tue, Nov 21, 2017 at 9:26 AM, Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote:<br>><br>> On 20/11/17 21:20, Ken Mankoff wrote:<br>>><br>>> Dear Helmut,<br>>><br>>> On 2017-11-02 at 07:17, Helmut Kudrnovsky <<a href="mailto:hellik@web.de">hellik@web.de</a>> wrote:<br>>>><br>>>> Ken Mankoff wrote<br>>>>><br>>>>> Does GRASS (or GDAL?, proj4?) have a function (or addon) to provide the<br>>>>> scale error  for a given cell?<br>>>><br>>>><br>>>> AFAIR there is now an area function in r.mapcalc in GRASS trunk to calculate<br>>>> the pixel area.<br>>><br>>><br>>> There is. I've just built GRASS 7.5. But the area() function just returns nsres*ewres for each cell, it does not consider projection errors.<br>><br>><br>> You could run area() both in your projected location and in a lat-long location (preferably using a similar datum), although you would have to be careful in the definition of the resolution in each location. <br><br></div>You can avoid problems with raster resolution by creating a vector grid with areas. Calculate the area sizes for the vector grid, then reproject the vector grid and calculate the area sizes again. Maybe the differences in area sizes for the same vector area give you what you are looking for.<br><br></div>Markus M<br><div><div><br>> Comparing the output should give you a notion of projection error.<br>><br>> For visualization you might want to try the Indicatrix mapper Plugin in QGIS...<br>><br>> Moritz<br>><br>><br>> _______________________________________________<br>> grass-user mailing list<br>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br><br></div></div></div>