[GRASS-dev] strange r.stats output -- bug?

Kirk Wythers kirk.wythers at gmail.com
Wed Oct 19 01:56:11 EDT 2011


Can anyone explain the behavior of r.stats when applied to a raster that has been aggregated using r.resamp.stats with the "mode" option? It seems to me that the output should look like the example below from r.stats on the original resolution map (250 meter resolution) in the fist example. However, this is not the case. When applied to the aggregated map (1000 meter resolution), r.stats reports binned cell values and counts. 

GRASS 6.4.1 (northcentralus_albersequalarea):~ > r.resamp.stats -w input=mnwimifnfftgk7wt175v3_250m output=mnwimifnfftgk7wt175v3_1000m method=mode --verbose 100%

This all appears to work properly. If I export the raster to MATLAB and simply view the cell values they are all whole number like 500, 600, 700, etc....

However, when I run r.stats on each map, I get:

GRASS 6.4.1 (northcentralus_albersequalarea):~ > r.stats -cl input=mnwimifnfftgk7wt175v3_250m 100%
0  893025
100  18419
120  34024
170  125
260  18
380  527
400  4269
500  36610
600  163
700  12479
800  51090
900  53455
990  5
999  527
* no data 874

GRASS 6.4.1 (northcentralus_albersequalarea):~ > r.stats -cl input=mnwimifnfftgk7wt175v3_1000m
100%
0-3.917647 from  to  929109
97.941176-101.858824 from  to  11467
117.529412-121.447059 from  to  32643
168.458824-172.376471 from  to  1
258.564706-262.482353 from  to  1
376.094118-380.011765 from  to  30
399.6-403.517647 from  to  383
497.541176-501.458824 from  to  24231
599.4-603.317647 from  to  12
697.341176-701.258824 from  to  2184
799.2-803.117647 from  to  50355
897.141176-901.058824 from  to  55190
995.082353-999 from  to  4

The output to r.stats from the 250 meter resolution map look fine... but doesn't make sense to me that the 1000 meter resolution map would have it's count and cell values binned into groups like 0-3.9, and 97.9-101.9, when I was using the "method=mode" flag in r.resamp.stats???

What am I missing here?

Thanks in advance



More information about the grass-dev mailing list