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. _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
