Apart from the region resolution, it's land cover classes, r.fillnulls
will produce nonsense data: floating point. Rather use r.neighbors
with a modal filter, patch the original with the filtered map, that
would replace NULLs with the most common surrounding land cover type.
Markus M
Hanlie Pretorius wrote:
Hi,
I have a MODIS land cover image that contains NULL values according to
r.univar:
-
r.univar map=modis_land_cover_type1_2001_prim...@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_ty...@permanent
output=MCD12Q1.A2001001.005_land_cover_type1_filled
Locating and isolating NULL areas...
Reading input raster map r_fillnulls_1...@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_prim...@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_prim...@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
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user