Tutoriel: Introduction à GDAL
Conditions d’achèvement
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
).