0

I am trying to change the extent of a rasterized shapefile to match that of another raster so it can act as a mask for the latter. Other posts suggest using gdalwarp. However, when using gdalwarp to create a new output it throws me a TypeError: argument of type 'float' is not iterable. Is gdalwarp unable to use float values? Although I can't see this written anywhere as a limitation of the function. Below is my code:

slope_src = r'C:\data\extraction\slope.tif'
mask_src = r'C:\data\buildings\blds.tif'
mask_new = r'C:\data\buildings\blds_mask.tif'
slope_raster = gdal.Open(slope_src)
ulx, xres, xskew, uly, yskew, yres = slope_raster.GetGeoTransform()
xres = 0.4
yres = 0.4
lrx = ulx + (slope_raster.RasterXSize * xres)
lry = uly + (slope_raster.RasterYSize * yres)
subprocess.call(["C:\\PRogram Files\\GRASS GIS 7.3\\extrabin\\gdalwarp.exe", "-te", lrx, lry, ulx, uly, mask_src, mask_new])
Improve this question
4
  • 1
    It seems to me you have wrong order in te key (-te xmin ymin xmax ymax). ulx == xmin, lrx == xmax, lry == ymin, uly == ymax. Try to swap x in your command line like this: subprocess.call(["C:\\PRogram Files\\GRASS GIS 7.3\\extrabin\\gdalwarp.exe", "-te", ulx, lry, lrx, uly, mask_src, mask_new])
    – Dmitry Baryshnikov
    Aug 7, 2017 at 8:02
  • This removed the error, but now the issue is that 'mask_new' is not created by gdalwarp.
    – Tins
    Aug 7, 2017 at 8:42
  • 1
    From your code I don't understand the sense of using gdalwarp. To extract part of the image the gdal_translate with -projwin ulx uly lrx lry key are applicable. gdalwarp is for reprojecting from one SRS to another or for reference based on GCP.
    – Dmitry Baryshnikov
    Aug 7, 2017 at 9:27
  • With the use of '.SetGeoTransform()' I was able to match the extents of the rasterized shapefile by passing it with the results of '.GetGeoTransform() applied to the raster to be masked.
    – Tins
    Aug 15, 2017 at 7:06

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.