<div dir="ltr"><div><div>Thanks all to have pointed me to the right pages. I solved with the r.univar command:<br><br>r.univar -t map=forestArea zones=regionId separator=comma output=forestArea_by_region.csv<br><br></div>Cheers,<br></div>  /Antonello<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On 19 September 2016 at 21:59, Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</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 Mon, Sep 19, 2016 at 6:41 PM, Antonello Lobianco<br>
<<a href="mailto:antonello@lobianco.org">antonello@lobianco.org</a>> wrote:<br>
> Let's say in grass I have a raster layer "regionId" and a raster layer<br>
> "forestArea".<br>
><br>
> In the regionId I have stored the ids of the regions, in the forestArea I<br>
> stored the square km of forest within the pixel.<br>
><br>
> How can I print a table summing up the forest area by region, that is the<br>
> sum of the values of a layer aggregated ("group by") by a second layer ?<br>
<br>
</span>Please take a look at<br>
<a href="https://grasswiki.osgeo.org/wiki/Zonal_statistics" rel="noreferrer" target="_blank">https://grasswiki.osgeo.org/<wbr>wiki/Zonal_statistics</a><br>
<br>
There are several options.<br>
<br>
Best<br>
<span class="HOEnZb"><font color="#888888">Markus<br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Antonello Lobianco<br>AgroParisTech, Laboratoire d'Economie Forestière<br>14 Rue Girardet - 54000 Nancy, France<br>Tel: +33.383396865<br>Email: <a href="mailto:antonello.lobianco@nancy.inra.fr" target="_blank">antonello.lobianco@nancy.inra.fr</a><br><a href="http://antonello.lobianco.org/" target="_blank">http://antonello.lobianco.org</a></div></div>
</div>