Tutorial: Introductie GDAL
3. Herprojecteren van GIS data
3.1. Herprojecteer rasterdata
GDAL heeft de mogelijkheid om een rastercoördinatensysteem te wijzigen met behulp van de volgende syntaxis:
gdalwarp
-t_srs EPSG:... <input> <output>
Het -t_srs argument specificeert het doelcoördinatensysteem. Als het broncoördinatensysteem onbekend is, moet het worden opgegeven met het -s_srs
argument. EPSG:... specificeert de EPSG-code van de projectie. <input> and <output> zijn respectievelijk de invoer- en uitvoergegevens.
We gaan nu een digitaal hoogtemodel van de Shuttle Radar Topography Mission (SRTM) herprojecteren. Je kan hoogtemodellen voor je eigen gebied downloaden van USGS Earth Explorer. Hier gebruiken we de data behorende bij deze cursus.
Om het hoogtemodel te herprojecteren van WGS-84 lat/lon naar Amersfoort/RD Nieuw gebruiken we dit commando:
gdalwarp -t_srs EPSG:XXXXX srtm_37_02.tif dem_rd.tif
XXXXX met de EPSG code voor Amersfoort/RD Nieuw
(zie een van je eerdere antwoorden met ogrinfo).2. Voer dit commando uit:
gdalwarp -t_srs EPSG:28992 srtm_37_02.tif dem_rd.tif <ENTER>
3. Visualiseer het resultaat in QGIS.
Meer informatie over het gdalwarp-commando is hier te vinden.