Suite

Requête du fichier de formes pour crop_to_cutline à l'aide de gdalwarp


J'ai un fichier de formes avec plusieurs polygones. Je voudrais découper une série de rasters mais n'utiliser qu'une seule entité de la classe d'entités shapefile. J'utilise gdal_warp. Il existe une option pour interroger la classe d'entités à l'aide de sql en utilisant l'option -csql dans l'outil.

Quelqu'un peut-il fournir un exemple de ce à quoi cela ressemble dans la ligne de commande?

Peut-être n'est-il utilisé que lors de la connexion à une base de données spatialement activée en tant que source de données ?


Voici ma réponse :

Bien que je ne comprenne pas complètement comment fonctionnent les calques dans gdal, j'ai trouvé qu'ils sont presque toujours le nom de fichier sans chemin ni extension. Par exemple : "C:/Folder/AllRiverMileSegments.shp" serait chargé en tant que couche "AllRiverMileSegments".

L'écriture d'une instruction SQL pour gdalwarp utilise cette couche comme nom de table et appelle les colonnes de la table attributaire en tant que colonnes de cette table. Dans mon fichier de formes ci-dessus, je souhaite interroger des segments de miles fluviaux. Dans ma table attributaire, j'ai une colonne appelée "RiverMile" que j'utiliserai pour interroger les segments de miles fluviaux de la classe d'entités. La section de gdalwarp qui appelle cette instruction ressemblera à ceci :

-csql "SELECT * FROM AllRiverSegments O RiverMile = $RIVERMILE"

Où $ RIVERMILE est le nombre de milles fluviaux que je veux interroger.

Merci à MerseyViking de m'avoir fait démarrer ici.


Voir la vidéo: Read and write vector files with GDALOGR in Python (Octobre 2021).