[gdal-dev] Re: GDAL interpolation of Gtiff / ESRI polygon.shp

[RolandB]

Dear Kyle,

thank you for your fast response. I had been trying adapt gdalwarp() for
some time now and it seems it does not want to interpolate.

This what I did
$user:gdal_rasterize -a_nodata -99 -a kennummer -te 5.93 50.39 6.74722
51.22886 -ts 59 94 -ot Int32 -l Bk50_rur_wgs ./Input/bk50_rur_wgs.shp
./tmp/output1.tif
0...10...20...30...40...50...60...70...80...90...100 - done.

$user: gdal_translate -of netCDF ./tmp/output1.tif ./tmp/output1.nc
Input file size is 59, 94

$user: gdalwarp -ot Int32 -ts 59 94 -r near ./tmp/output1.tif
./tmp/output.tif
Creating output file that is 59P x 94L.
Processing input file ./tmp/output1.tif.
Using internal nodata values (eg. -99) for image ./tmp/output1.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.

$user: gdal_translate -of netCDF ./tmp/output.tif ./tmp/interpolated.nc
Input file size is 59, 94

The nodata which is previously marked by -99 as in file ./tmp/output1.tif &
./tmp/output1.nc are set to 0 by the gdalwarp and show up like that in the
interpolated.nc. Not any more as nodata, as in ./tmp/output1.nc

Is there anything you perceive which I should change to get an interpolated
grid data set over the values of -99 with the second last and the last
command line?

Best regards,
Roland

--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/GDAL-interpolation-of-Gtiff-ESRI-polygon-shp-tp7114807p7150275.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.