Tutoriel: Introduction à GDAL: Requêtes spatiales de données vectorielles

4. Convertir des formats SIG

4.3. Requêtes spatiales de données vectorielles

Pour notre carte de Delft, nous voulons faire l'analyse SIG suivante:

Sélectionnez la commune de Delft sur la carte des communes et l´enregistrez dans un nouveau fichier de formes;
Intersectez les limites de la commune de Delft avec la carte routière des Pays-Bas.

Nous pouvons utiliser une requête spatiale pour sélectionner une entité à partir d'une carte vectorielle.

1. Quel est l'attribut de la carte des communes contenant les noms des communes? Vous pouvez utiliser ogrinfo ou QGIS pour répondre à cette question.

2. Exécutez la commande suivante:

ogr2ogr
-f "ESRI Shapefile" -where GM_NAAM='Delft' 
-a_srs EPSG:28992 delft.shp gem_2011_gn1.shp

Cela enregistrera la fonctionnalité avec GM_NAAM Delft dans un nouveau fichier appelé delft.shp. L'argument -a_srs EPSG:28992 est utilisé pour affecter la projection Amersfoort / RD New au fichier de sortie. L'argument -f définit le format de sortie.

3. Ouvrez maintenant dans QGIS le DEM reprojeté (dem_rd.tif), la carte de route reprojetée (roadsreprojected.shp) et la commune de Delft (delft.shp).