[GRASS-user] raster map calculations give faulty output

Annekatrien Debien annekatrien.debien at gmail.com
Thu Feb 26 04:04:54 EST 2009


Dear all,

when working with Grass 6.3 as a native Windows version, every raster map
calculation I make goes wrong. The whole calculation seems to go well, I
don't get any errors, but when I take a look at r.info for the newly
calculated map, the range is always "min = -1.#QNAN0  max = -1.#QNAN0 ",
regardless whether the original map is CELL or FCELL. In the original map
however, the range is correct. This miscalculation happens with different
commands:

r.mapcalc "red_rad=((235+4.5)/(255-1))*(mar2001_B3-1)-4.5"
r.resamp.stats ndvi_mar2001_idr out=ndvi_landsat method=average
r.mapcalc "red=(mar2001_B3-10)"
r.mapcalc "ndvi_mar2001=1.0*(mar2001_B4-red)/(mar2001_B4+red)"
i.tasscap -7 band1=mar2001_B1 band2=mar2001_B2 band3=mar2001_B3
band4=mar2001_B4 band5=mar2001_B5 band7=mar2001_B7 outprefix=tasscap

Did anyone ever encounter this problem and found a solution to it? I've also
tried it in Grass 6.4 on Windows, and the same happens.

Thanks!
Annekatrien
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20090226/146066e3/attachment-0001.html


More information about the grass-user mailing list