Is it possible to clip a shapefile from a raster layer using gdalwarp?

Question

I have a shapefile with a single geometry as the following:

I also have a GeoTIFF raster that I can put behind this vector layer... I'd like to find a way of programmatically clipping this shapefile feature from my GeoTIFF file, in a way that I get a png image as a result.

The closest solution I found to this problem is this answer. This solution solves the problem using gdalwarp with the following command line:

gdalwarp -of GTiff -cutline area_of_interest.shp -cl area_of_interest  -crop_to_cutline PCE_in_gw.asc  data_masked7.tiff

In my case, I only have the Shapefile and the GeoTIFF file. I didn't get what I'm supposed to put on the .asc file and the area_of_interest parameter in this command line sample. What extra procedures do I have to do to clip my GeoTIFF with a Shapefile while using gdalwarp?

score 2 · Answer 1 · edited Dec 17 '21 at 22:25

2

Look this example in Pyhton and Jupyter Notebook:

gdal.Warp(destNameOrDestDS = outputpath, 
          srcDSOrSrcDSTab  = intputpath, 
          cutlineDSName    = shapefile or geojson, 
          cropToCutline    = True,
          copyMetadata     = True,
          dstNodata        = 0)

edited Dec 17 '21 at 22:25

Kadir Şahbaz

76,800
56
247
389

answered Dec 17 '21 at 22:24

Helios

427
2
8

score 2 · Accepted Answer · answered Dec 17 '21 at 22:42

2

From the command line you can do the following:

gdalwarp -of GTiff -cutline input.shp -crop_to_cutline input.tif input_clip.tif

Another alternative is with pkcrop from pktools

pkcrop -i input.tif -e input.shp -cut -o input_clip.tif

answered Dec 17 '21 at 22:42

dmci

4,882
2
19
31

Is it possible to clip a shapefile from a raster layer using gdalwarp?

2 Answers2