Tutorial: Introductie GDAL

GDAL kan een rasterbestand herprojecteren met de volgende syntaxis:

gdal raster reproject -d EPSG:... <input> <output>

Het argument -d geeft het doelcoördinatensysteem aan. Als het broncoördinatensysteem onbekend is, moet dit worden opgegeven met -s. EPSG:… staat voor de EPSG‑code van de projectie. <input> en <output> zijn respectievelijk de invoer‑ en uitvoerbestanden.

We gaan nu een Digital Elevation Model (DEM) herprojecteren dat is verkregen via de Shuttle Radar Topography Mission (SRTM). Je kunt DEM’s voor je eigen studiegebied downloaden via USGS Earth Explorer of de OpenTopography DEM Downloader plugin in QGIS. In deze oefening gebruiken we de dataset die bij de cursus wordt geleverd.

Om het DEM te herprojecteren van WGS‑84 lat/lon naar Amersfoort / RD New, gebruiken we het volgende commando:

gdal raster reproject -d EPSG:XXXXX srtm_37_02.tif dem_rd.tif

1. Vervang XXXXX door de juiste EPSG‑code voor Amersfoort/RD New (zie één van je eerdere antwoorden met gdal info).

2. Voer het volgende commando uit:

gdal raster reproject -d EPSG:28992 srtm_37_02.tif dem_rd.tif <ENTER>

3. Visualiseer het resultaat in QGIS.

Meer informatie over het gdal raster reproject‑commando vind je hier.

Merk op dat standaard nearest neighbour resampling wordt gebruikt. Bekijk de argumenten van het commando in de documentatie als je een andere resampling‑methode wilt gebruiken of de uitvoercelgrootte wilt instellen.