[GRASS-user] odd r.stats output, after r.resample.stats
Kirk Wythers
kirk.wythers at gmail.com
Tue Oct 18 15:31:31 EDT 2011
I am trying to aggregate resample a 250 meter resolution raster to 1000m resolution as the most frequently occurring value in the 250m map. I am using:
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???
What am I missing here?
Thanks in advance
More information about the grass-user
mailing list