Sélection d'entités par attributs à l'aide d'ogr2ogr?

J'écris un fichier de formes après un traitement; cependant, dans la dernière étape, je dois faire une sélection par attributs à partir du fichier de formes.

J'utilise la commande dans un shell dans le but de l'utiliser dans un script Python quand cela fonctionne.

ogr2ogr -f "ESRI Shapefile" -select * where ID="1" outfile.shp infile.shp

Je reçois le message d'erreur:

FAILURE: 
Unable to open datasource `Downloads' with the following drivers.

Que pourrais-je faire de mal?

— user2757128
source

Est-ce vraiment le message exact pour cette commande?

— BradHards

Vous manquez un signe moins avant whereet ce selectn'est pas nécessaire, il devrait donc être:

ogr2ogr -where ID="1" outfile.shp infile.shp

ou si vous devez effectuer une requête plus complexe sur vos données d'entrée:

ogr2ogr -sql "SELECT * FROM infile WHERE ID='1'" outfile.shp infile.shp

Si IDest un champ de type Entier, remplacez ID='1'par ID=1.

Remarques:

-f "ESRI Shapefile"n'est pas nécessaire car il "ESRI Shapefile"s'agit du ogr2ogrformat de sortie par défaut;
il est pratique d'ignorer -selectet d'utiliser directement la -whereclause lorsque vous souhaitez sélectionner tous les champs.

— Antonio Falciano
source

Merci beaucoup ça marche. J'utilise: ogr2ogr -where "ID = '1'" output.shp input.shp

— user2757128

Comment puis-je utiliser ce code: ogr2ogr -sql "SELECT * FROM infile WHERE ID = '1'" outfile.shp infile.shp "en python?

— Shiuli Pervin

@ShiuliPervin vous pouvez utiliser l'instruction en python en utilisant os.system('''ogr2ogr ... ''') . assurez-vous d'être import osen haut de votre script python

— geoeye

@afalciano, pour l'instruction ogr2ogr -sql "SELECT * FROM infile WHERE ID='1'" outfile.shp infile.shpau lieu de simplement ID = 1, existe-t-il un moyen de fournir un vecteur ou un ensemble de nombres, par exemple ogr2ogr -sql "SELECT * FROM infile WHERE ID IN ['1','5','29']" outfile.shp infile.shp?

— hlm

@afalciano une autre chose ... au lieu d'écrire un autre fichier shp, existe-t-il un moyen d'obtenir simplement les coordonnées des polygones ID sélectionnés?

— hlm