[GIS] Clipping raster with vector layer using GDAL

gdalpython

I have installed GDAL using the Osgeo installer. How can I clip a raster layer with a vector layer programmatically? Is there any GDAL API available which can help me with this? I am using Python.

Best Answer

I'm not sure about the gdal api, there is void* GDALWarpOptions::hCutline in the Warp Options referenced from the Warp API tutorial, but no explicit examples. Are you sure you need a programmatic answer? The command line utilities can do it out of the box:

create a shapefile containing just the area of interest clipping polygon
use ogrinfo to determine the extent of the clipping shapefile
use gdal_translate to clip to the shape extents
use gdalwarp with -cutline parameter

Steps 2 & 3 are for optimization, you could get by with just gdalwarp -cutline ....

See Clipping rasters with GDAL using polygons from Linfinity for linux based solution all wrapped up in one script. Another cutline example can be seen in Michael Corey's tutorial creating hillshades for Mapnik.

Related Solutions

[GIS] GDAL and Python: how to create new files for (e.g.) clipping with gdalwarp

I finally came up with a quite nice result calling GDAL expressions from Python using "Call ()" method and specifying the path with "os.chdir ()" method. Here is a script example:

import os
from subprocess import call
call(["ls", "-l"])

path= 'C:/users/unimi/Desktop/'

os.chdir(path)

cmd= 'gdalwarp -dstnodata 0 -q -cutline \
    study_area_foscagno.shp -crop_to_cutline -of GTiff sept2000_B5_refl.tif clip_sept2000_B5_refl.tif'
call (cmd)

The "Call ()" method is here used to do all the GDAL stuff without even importing GDAL in Python ("Call ()" method is even better than "os.system ()" method, see "https://stackoverflow.com/questions/89228/calling-an-external-command-in-python"). "os.chdir ()" method than have let me set the path where the program should look for my data. Hope this could help.

GDAL Raster Clipping – How to Clip a Raster with Another Raster Using GDAL

I don't know if it's possible to clip a raster with an other raster but you could use gdaltindex to build the shapefile with the extent of your raster.

http://www.gdal.org/gdaltindex.html

Best Answer

Related Solutions

[GIS] GDAL and Python: how to create new files for (e.g.) clipping with gdalwarp

GDAL Raster Clipping – How to Clip a Raster with Another Raster Using GDAL

Related Question