Hi Alessandro, You could try v.rast.bufferstats [1], if you do not have hundrets of thousands of points….
If the buffers don`t overlap, you could use v.rast.stats [2] from the upcoming GRASS 7.6 release, that will allow multiple raster input… Otherwise have a look at the general introduction to GRASS and Python [3]. Cheers Stefan 1: https://grass.osgeo.org/grass74/manuals/addons/v.rast.bufferstats.html 2: https://grass.osgeo.org/grass76/manuals/v.rast.stats.html 3: https://grasswiki.osgeo.org/wiki/GRASS_and_Python From: grass-user <grass-user-boun...@lists.osgeo.org> On Behalf Of Alessandro Sebastiani Sent: mandag 29. oktober 2018 16:06 To: grass-user@lists.osgeo.org Subject: [GRASS-user] zonal statistics for multiple areas Hello to everybody, I hope my question is appropriate for this mail list. I have created 5 different buffer layers (d=100,200,300,400,500 m) from a point vector. Now i want to compute some zonal statistics using different rasters as input. I know how to do that separately, but i was wondering how could i automate this procedure using a python script. I know python's basics, but i was not able to do so following guidelines that i found on the internet, e.g GRASS tutorial. Thank you in advance, A
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
