gdal/ogr: How to really crop a shapefile?

gis gdal shapefile ogr
Hugolpz picture Hugolpz · Aug 19, 2013 · Viewed 7.8k times · Source

Given a SHP file corresponding to European countries, and...

Given defined area corresponting to France such :

West : 005° 48' W
East : 010° E
North : 051° 30' N
South : 041° N

How to get only the dots/geometries which intersects my defined area with gdal ? so my crop indeed is a real crop, which just keep the necessary geometries. Strong preference to gdal, ogr or console solutions.

Answer

Mike T picture Mike T · Aug 21, 2013

Use GDAL's ogr2ogr command-line utility. If you have a file Europe.shp which has a spatial reference with units of degrees, then use the -clipsrc option with decimal degrees to make a new shapefile:

ogr2ogr -clipsrc -5.8 41 10 51.5 France.shp Europe.shp

There is also an example at the ogr2ogr page with France used as an example.

Related questions

Rasterizing a GDAL layer

Edit Here is the proper way to do it, and the documentation: import random from osgeo import gdal, ogr RASTERIZE_COLOR_FIELD = "__color__" def rasterize(pixel_size=25): # Open the data source orig_data_source = ogr.Open("test.shp") # Make a …

Read More
How to use GDAL to create geotiff from tiff and 4 corners latitude and longitude

I have an image (map) without geo data in TIFF format. I need to get GeoTIFF file from my image. I have latitude and longitude for each corner of my map. How can I add my geo data to my …

Read More
Reading hdf files into R and converting them to geoTIFF rasters

I'm trying to read MODIS 17 data files into R, manipulate them (cropping etc.) and then save them as geoTIFF's. The data files come in .hdf format and there doesn't seem to be an easy way to read them into R. …

Read More