<div dir="ltr">Dear all, i was able to compute zonal statistics for multiple areas using this script. I used <div>> import grass.script as gscript</div><div>> v_list = gscript.read_command("g.list" , type = "vect").strip().split("\n")</div><div># i have 5 vectors in my list</div><div>> method = ("sum" , "average")</div><div>> for v in v_list:</div><div> >    gscript.run_command("v.rast.stats", map_=v, raster = "LC1", column_prefix = "LC1", method = methods)</div><div><br></div><div>Now i have like 4 rasters that i would like to use as input to compute my statistics. Is there a way to loop this process for rasters and vector at the same time? i tried creating a list for rasters and using the name list as raster input but it seems not to work! cheers, Alessandro</div></div><br><div class="gmail_quote"><div dir="ltr">Il giorno mar 30 ott 2018 alle ore 11:52 Alessandro Sebastiani <<a href="mailto:alessandro.sebastiani@uniroma1.it">alessandro.sebastiani@uniroma1.it</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thanks a lot for your help!</div><br><div class="gmail_quote"><div dir="ltr">Il giorno mar 30 ott 2018 alle ore 11:49 Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net" target="_blank">nik@nikosalexandris.net</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">* Alessandro Sebastiani <<a href="mailto:alessandro.sebastiani@uniroma1.it" target="_blank">alessandro.sebastiani@uniroma1.it</a>> [2018-10-29 21:32:51 +0100]:<br>
<br>
>thank you all for your suggestions.<br>
>Dear Nikos, I have some 0-1 rasters that represents presence-absence of<br>
>different land covers. Each raster's resolution is 10x10m, thus i want to<br>
>compute the sum in order to obtain the surface covered by each land cover<br>
>within my buffers. As regards to the buffer, i want to use the 0-100,<br>
>0-200, 0-300, 0-400 and 0-500 for computing stats. I can do it separately<br>
>with the module v.rast.stats. I just don't know how to loop it in a python<br>
>script<br>
<br>
Dear Alessandro,<br>
<br>
I don't have time to write back as I would like. Instead, I will send<br>
you off-list some unpublished script.  With a bit of reading, there, you<br>
might identify some ideas that will eventually help in what you are<br>
doing.<br>
<br>
Kind regards, Nikos<br>
</blockquote></div>
</blockquote></div>