[GRASS-user] MASK does not do it's work (?) with r.mapcalc

nikos.alexandris at felis.uni-freiburg.de nikos.alexandris at felis.uni-freiburg.de
Tue Dec 18 20:13:24 EST 2007


OK... I am confused!

(An example) With r.mapcalc I create:

(a) NDVIDIFF0 based on "if(NDVIDIFF > 0, StudyArea_Raster,
null())"

{ Data type and r.describe of NDVIDIFF0: CELL, * 1-23 25-30
32-42 }

(b) and NDVIDIFF0_fromNDVI based on "if(NDVIDIFF > 0,
NDVIDIFF, null())"

{ Data type and r.descrive of NDVIDIFF0_fromNDVI: DCELL, *
0.000000-0.636731 0.642682-0.702189 0.714091-0.743845
0.749795-0.755746 
[...] 
1.309167-1.315117 1.315117-1.321068 1.451985-1.457936
1.511492-1.517443 }

Creating the MASK with r.mask or with r.mapcalc=MASK based
on (a) is always successful. Creating the MASK with (b) is
not what I expect and, in addition, it gives different
results when using r.mapcalc MASK=... instead of r.mask.

Do the DCELL values of (b) affect the process... ?



Thank you,

Nikos.


More information about the grass-user mailing list