<div dir="ltr"><br><br>On Wed, Aug 1, 2018 at 11:59 AM, Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br>><br>> `r.statistics` which was/is meant for categorical raster maps, does not<br>> provide for a 'count' method.<br>><br>> `r.stats.zonal`, on the other hand, which is meant for floating point<br>> raster maps, does so.<br>><br>> This is rather unexpected. 'Count'ing pixels should rather derive 'CELL'<br>> maps and not 'FCELL'.<br>><br>> At the moment, I require an extra mapcalc using the int() function.<br>><br><div>> Any comments?</div><div><br></div><div>Counting pixels can produce a result that is too large for CELL (32 bit integer), therefore storing the count as DCELL is a bit safer because it can exactly store integers larger than the 32 bit maximum.</div><div><br></div><div>In short, use r.stats.zonal rather than r.statistics.</div><div><br></div><div>Markus M</div><div><br></div>><br>> Nikos<br>><br>> _______________________________________________<br>> grass-dev mailing list<br>> <a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-dev">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br><br></div>