On 22/01/08 14:22, Nikos Alexandris wrote:
I still can't explain why I can't create a mask with the following raster:
+----------------------------------------------------------------------------+
| Layer: NDVI20060.35 Date: Tue Jan 22 13:53:02 2008
|
| Mapset: nik Login of Creator: nik
|
| Location: Peloponnese_HGRS87
|
| DataBase: /home/nik/grassWorks/Peloponnese_DamageAssessment2007
|
| Title: ( NDVI20060.35 )
|
| Timestamp: none
|
|----------------------------------------------------------------------------|
|
|
| Type of Map: raster Number of Categories: 255
|
| Data Type: DCELL
|
| Rows: 1048
|
| Columns: 972
|
| Total Cells: 1018656
|
| Projection: Transverse Mercator
|
| N: 4255830 S: 3993750 Res: 250.07633588
|
| E: 468810 W: 225720 Res: 250.09259259
|
| Range of data: min = 0.001304 max = 0.349991
|
|
|
| Data Description:
|
| generated by r.mapcalc
|
|
|
| Comments:
|
| if((NDVI2006 > 0 && NDVI2006 < 0.350000), NDVI2006, null())
|
|
|
+----------------------------------------------------------------------------+
Whenever I "r.mask input=NDVI20060.35" or with "rmapcalc MASK=NDVI20060.35",
I get nothing when I try to view or to use r.mapcalc on others rasters.
Any ideas about this?
Maybe the value of the MASK map has to be 1 ? I see that your MASK
varies from 0.001304 to 0.349991.
Moritz
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
