Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It only takes a minute to sign up.
Sign up to join this community
I am writing a shapefile after some processing; however, in the final step I need to do some selection by attributes from the the shapefile.
I am using the command in a shell with the aim of using it in a Python script when it works.
ogr2ogr -f "ESRI Shapefile" -select * where ID="1" outfile.shp infile.shp
I am getting the error message:
FAILURE:
Unable to open datasource `Downloads' with the following drivers.
What could I be doing wrong?
You miss a minus sign before where and the select is not necessary, so it should be:
ogr2ogr -where ID="1" outfile.shp infile.shp
or if you have to do more complex query on your input data:
ogr2ogr -sql "SELECT * FROM infile WHERE ID='1'" outfile.shp infile.shp
If ID is a field of Integer type, substitute ID='1' with ID=1.
Notes:
-f "ESRI Shapefile" is not necessary because "ESRI Shapefile" is the ogr2ogr default output format; If using GDAL >= 2.3, the format is guessed from the extension.-select and use directly the -where clause when you want to select all the fields.
os.system('''ogr2ogr ... ''') . be sure to import os at the top of your python script
ogr2ogr -sql "SELECT * FROM infile WHERE ID='1'" outfile.shp infile.shp instead of just ID=1 is there a way to provide a vector or a set of numbers instead, e.g. ogr2ogr -sql "SELECT * FROM infile WHERE ID IN ['1','5','29']" outfile.shp infile.shp?
ogr2ogr -sql "SELECT * FROM infile WHERE ID IN [1,5,29]" outfile.shp infile.shp Second question: to get the coordinates of the selected geometry without creating another shp: ogrinfo -sql "SELECT OGR_GEOMETRY FROM infile WHERE ID IN [1,5,29]" infile.shp
May 15, 2017 at 8:43