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

Moritz Lennert mlennert at club.worldonline.be
Tue Jan 22 11:58:32 EST 2008


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


More information about the grass-user mailing list