[GRASS-user] Filling Null values in a MODIS image

Hanlie Pretorius hanlie.pretorius at gmail.com
Wed Aug 4 04:29:16 EDT 2010


Hi,

I have a MODIS land cover image that contains NULL values according to r.univar:
-----
r.univar map=modis_land_cover_type1_2001_primary at PERMANENT
total null and non-null cells: 41884
total null cells: 17145
-----

So, I tried to fill it with r.fillnulls and I get the following:
-----
r.fillnulls input=MCD12Q1.A2001001.005_land_cover_type1 at PERMANENT
output=MCD12Q1.A2001001.005_land_cover_type1_filled
Locating and isolating NULL areas...
Reading input raster map <r_fillnulls_1911 at PERMANENT>...
Writing output raster map <r_fillnulls_1911.buf>...
Creating interpolation points...
Extracting points...
Building topology for vector map <vecttmp_fillnulls_1911>...
Registering primitives...
0 primitives registered
0 vertices registered
Building areas...
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
Number of nodes: 0
Number of primitives: 0
Number of points: 0
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
r.to.vect complete.
Interpolating 0 points
Not sufficient points to interpolate. Maybe no hole(s) to fill in the
current map region?
Removing raster <MASK>
Removing raster <r_fillnulls_1911>
Removing raster <r_fillnulls_1911.buf>
Removing raster <r_fillnulls_1911_filled>
Raster map <r_fillnulls_1911_filled> not found
<r_fillnulls_1911_filled> nothing removed
Removing vector <vecttmp_fillnulls_1911>
(Wed Aug  4 10:18:22 2010) Command finished (0 sec)
-----

Output of g.region -p:
-----
g.region -p rast=modis_land_cover_type1_2001_primary at PERMANENT res=0.00423094
projection: 3 (Latitude-Longitude)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      27:19:08.833466S
south:      28:30:53.475952S
west:       28:03:57.840271E
east:       28:41:37.727966E
nsres:      0:00:15.210751
ewres:      0:00:15.269511
rows:       283
cols:       148
cells:      41884
-----

Output of r.info:
-----
r.info map=modis_land_cover_type1_2001_primary at PERMANENT
 +----------------------------------------------------------------------------+
 | Layer:    modis_land_cover_type1_2001_p  Date: Thu Jul 29 15:17:55 2010
 | Mapset:   PERMANENT                      Login of Creator: hanlie
 | Location: world_wgs84
 | DataBase: /media/0847147784/grassdata
 | Title:     ( MCD12Q1.A2001001.005_land_cover_type1 )
 | Timestamp: none
 |----------------------------------------------------------------------------
 |
 |   Type of Map:  raster               Number of Categories: 14
 |   Data Type:    CELL
 |   Rows:         11957
 |   Columns:      6277
 |   Total Cells:  75054089
 |        Projection: Latitude-Longitude
 |            N: 27:19:08.833466S    S: 28:30:53.475952S   Res: 0:00:00.36001
 |            E: 28:41:37.727966E    W: 28:03:57.840271E   Res: 0:00:00.36002
 |   Range of data:    min = 1  max = 14
 |
 |   Data Description:
 |    generated by r.in.gdal
 |
 |   Comments:
 |    r.in.gdal input="F:\Hanlie\UCT\M.Sc\Data\modis\MOD12Q1\MCD12Q1.A2001\
 |    001.h20v11.005.2009342153812_land_cover_type1.tif" output="MCD12Q1.A\
 |    2001001.005_land_cover_type1"
 |
 +----------------------------------------------------------------------------+


I have tried to see the NULLS by displaying everything in white,
except for NULLS displayed in red, but I didn't see anything.

Does anyone know how I can find out if there are actually NULL value
and where they are?

Thanks
Hanlie


More information about the grass-user mailing list