Tutorial: Introductie GDAL: Ruimtelijke queries met vectordata

4. Converteren van GIS formaten

4.3. Ruimtelijke queries met vectordata

Voor onze kaart van Delft willen we de volgende GIS analyses doen:

Selecteer de gemeente Delft van de gemeentenkaart en sla het op in een nieuwe shapefile;
Gebruik de gemeentegrens van Delft om de wegen uit te knippen, zodat we alleen de wegen van Delft hebben.

Met een ruimtelijke querie kunnen we een object uit een vectorlaag selecteren.

1. Welk attribuut wordt in de gemeentenkaart gebruikt voor de gemeentenaam? Gebruik ogrinfo om het antwoord te vinden.

2. Voer het volgende commando uit:

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

Dit commando slaat het object met GM_NAAM Delft op in een bestand met de naam delft.shp. Het argument -a_srs EPSG:28992 gebruiken we om

de Amersfoort/RD Nieuw projectie toe te wijzen aan het uitvoerbestand. Met het argument -f definiëren we het uitvoerformaat.

3. We kunnen de roadsreprojected.shp vector nu knippen, zodat we alleen de wegen in de gemeente Delft over houden. We gebruiken hiervoor het -clipsrc argument:

ogr2ogr -clipsrc delft.shp roadsdelft.shp

roadsreprojected.shp

4. Open nu in QGIS het geherprojecteerde hoogtemodel (dem_rd.tif), de geknipte wegenkaart (roadsdelft.shp) en de gemeentegrens van Delft ( delft.shp).