Inconsistency d.what.rast & r.stats ?

Agustin Lobo alobo at
Fri Jul 16 12:32:42 EDT 1999

I'm computing the average values of a raster within  each category of 
raster and I think that I've found some inconsistent results.

First, I compute the per-class averages:

alobo/grscripts> r.average base=medvgeugoode1km.rec cover=test

Second, I check the result for, i.e., class 23:

alobo/grscripts> d.what.rast map=medvgeugoode1km.rec,
 Left:  what's here
 Right: quit

2676500(E) 4409500(N)
medvgeugoode1km.rec in user1  (23)S_SUPRAMEDITERRANEAN_ZONE
  in user1, quant   (111)
  in user1, actual  (110.794075)

Third, I compute the statistics for each class, but the results are not
the same as given by d.what.rast. For example, for class 23:

alobo/grscripts> r.stats -c in=medvgeugoode1km.rec,
r.stats:  100%
23 110.299171-110.744906 6082

While d.what.rast was giving 110.794075. Is this a bug?


Dr. Agustin Lobo
Instituto de Ciencias de la Tierra (CSIC)
Lluis Sole Sabaris s/n
08028 Barcelona SPAIN
tel 34 93409 5410
fax 34 93411 0012
alobo at

More information about the grass-user mailing list