[gdal-dev] how to find gdal_fillnodata

[Bryan Keith]

> Bryan Keith wrote:
>> Hello,
>>
>> I combined two DEM datasets with gdal_merge, and there's a broken strip
>> of
>> nodata between the two.  What gdal utility can I use to fix this?  I
>> found
>> references to gdal_fillnodata on the web, but I don't have it.  I
>> installed FWTools 2.2.8 on windows, but that doesn't seem to include
>> gdal_fillnodata either.  On Linux I have GDAL 1.4.4.0, released
>> 2007/11/23.  I looked at the source code (from the daily svn snapshot),
>> and it looks I need gdal.FillNoData but that's not part of the bindings
>> that I have on either Windows or Linux.  How can I get this utility or
>> is
>> there another utility that I can use?
>
> Bryan,
>
> This utility only exists in the most recent builds from trunk.  I see it
> was not in FWTools 2.3.0, but it is in the gdal-dev package of OSGeo4W.
> On linux you could build from svn with python support and you should have
> it.
>

Frank,

I installed OSGeo4W this afternoon.  I had to do a little futzing to get
python to work with gdal-dev (change the windows %path% variable and
change the Python sys.path list).  However, when I run gdal_fillnodata, I
end up with a tif that has all 0s.  What am I doing wrong?  Commands are
pasted below.

Bryan

C:\Projects>gdalinfo dem.tif -stats
Driver: GTiff/GeoTIFF
Files: dem.tif
Size is 2200, 1411
Coordinate System is:
Origin = (1749100.677571530000000,304762.295903306220000)
Pixel Size = (32.820441700754138,-32.820441700754138)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 1749100.678,  304762.296) (122d52'36.09"W, 38d30'0.63"N)
Lower Left  ( 1749100.678,  258452.653) (122d52'30.48"W, 38d22'22.86"N)
Upper Right ( 1821305.649,  304762.296) (122d37'27.85"W, 38d30'6.53"N)
Lower Right ( 1821305.649,  258452.653) (122d37'23.86"W, 38d22'28.75"N)
Center      ( 1785203.163,  281607.474) (122d44'59.57"W, 38d26'14.94"N)
Band 1 Block=2200x1 Type=Float32, ColorInterp=Gray
  Minimum=10.000, Maximum=563.300, Mean=78.248, StdDev=78.791
  NoData Value=-32766
  Metadata:
    STATISTICS_MINIMUM=10
    STATISTICS_MAXIMUM=563.29998779297
    STATISTICS_MEAN=78.248205065279
    STATISTICS_STDDEV=78.791386753552

C:\Projects>\OSGeo4W\apps\gdal-dev\bin\gdal_fillnodata.bat dem.tif
 dem_nnd.tif
0...10...20...30...40...50...60...70...80...90...100 - done.

C:\Projects>gdalinfo dem_nnd.tif -stats
Driver: GTiff/GeoTIFF
Files: dem_nnd.tif
Size is 2200, 1411
Coordinate System is:
Origin = (1749100.677571530000000,304762.295903306220000)
Pixel Size = (32.820441700754138,-32.820441700754138)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 1749100.678,  304762.296) (122d52'36.09"W, 38d30'0.63"N)
Lower Left  ( 1749100.678,  258452.653) (122d52'30.48"W, 38d22'22.86"N)
Upper Right ( 1821305.649,  304762.296) (122d37'27.85"W, 38d30'6.53"N)
Lower Right ( 1821305.649,  258452.653) (122d37'23.86"W, 38d22'28.75"N)
Center      ( 1785203.163,  281607.474) (122d44'59.57"W, 38d26'14.94"N)
Band 1 Block=2200x1 Type=Float32, ColorInterp=Gray
  Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
  Metadata:
    STATISTICS_MINIMUM=0
    STATISTICS_MAXIMUM=0
    STATISTICS_MEAN=0
    STATISTICS_STDDEV=0

C:\Projects>