[gdal-dev] how to find gdal_fillnodata

[Bryan Keith]

Hello,

See below.  I'm sending this again since I still don't have the issue
resolved.

Bryan

>> 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>
>