Tutoriel: Introduction à GDAL

4. Convertir des formats SIG

4.3. Requêtes spatiales sur des données vectorielles

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

Sélectionner la commune de Delft sur la carte des communes et la sauvegarder dans un nouveau fichier de forme (shapefile);
Intersecter 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 caractéristique dans une carte vectorielle.

1. Quel est l'attribut de la carte des communes contenant les noms des communes ? Vous pouvez utiliser soit 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

Ceci sauvegardera l'entité 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. Nous pouvons découper le vecteur roadsreprojected.shp afin qu'il ne couvre que la municipalité de Delft. Nous utilisons pour cela l'argument -clipsrc.

ogr2ogr -clipsrc delft.shp roadsdelft.shp roadsreprojected.shp

4. Ouvrez maintenant dans QGIS le DEM (MNE) reprojeté (dem_rd.tif), la carte routière découpée (roadsdelft.shp) et la commune de Delft (delft.shp).