geodaten mit gdal

Eine vollständige Dokumentation aller möglichen Befehle ist unter https://gdal.org zu finden

Raster

Infos über Datensätze gdalinfo

gdalinfo raster.tif   # Gibt Angaben zum raster.tif aus

Virtual Raster

Mit virtuellen Raster können schnell viele Einzeldateien in einer Datei zusammengefasst werden. Es handelt sich dabei um eine XML Datei, in welcher die Verknüpfungen zu den Einzeldateien (zB. GeoTiffs) definiert sind. Virtuelle Raster können mit QGIS oder GDAL als "normale" Raster geladen und bearbeiten werden. Der Befehl dazu ist gdalbuildvrt.

gdalbuildvrt OutputRaster.vrt data/*.tif

(dieser Befehl erstellt ein aus allen .tif Dateien im Unterodner data virtuelles Raster mit dem Namen OutputRaster.vrt )

Konvertieren von Rastern mit gdal_translate

Mit gdal_translate können Raster zwischen vielen Formaten konvertiert werden.

Die generelle Form des Befehl ist

gdal_translate -of <format> inputraster outputraster  # -of -> output format

Bei der Konvertierung können in gleichen Schritte folgende zusätzliche Parameter gesetzt werden:

Ausschnitt Raster:

gdal_translate -of GTIFF -projwin 2644807 1249396  2646299 1248478 -projwin_srs EPSG:2056 --oo TILED=YES swissimage_exar.vrt ausschnitt.tif

Höhenmodell ausschneiden, resize, reprojekt:

gdal_translate -of GTIFF -projwin 2598466 1233972  2612226 1226772 -projwin_srs EPSG:2056 -tr 2 2 ch.so.agi.lidar_2019.dtm.tif out.tif

Erstellen von xyz aus höhenmodell mit resize

gdal_translate -of XYZ -projwin 2598466 1233972  2612226 1226772 -projwin_srs EPSG:2056 -tr 5 5 ch.so.agi.lidar_2019.dtm.tif out.xyz

Höhenlinien mit gdal_contour

Beispiel (Aequidistanz = 5m)

gdal_contour -a elev raster.tif contour.shp -i 5

Reproject Raster

Mit gdalwarp können Raster zwischen verschiedenen Koordinatensystemen (SRS = spatial reference system) konvertiert werden. Gängige SRS sind:

  • EPSG:21781 -> LV03
  • EPSG:2056 -> LV95 (default in CH)
  • EPSG:4326 -> WSG84 (zB. GPS Koordinaten)
gdalwarp -t_srs EPSG:4326 input.tif output.tif

Vektordaten

Infos über Datensätze ogrinfo

ogrinfo vectordaten.gpkg  # oder *.shp, dxf ... 

Gängige Parameter:

  • -so -> summary only
  • -al -> all layer (zB. bei Geopackages relevant, welche mehrere Layer enthalten

Konvertieren von Vektordaten mit ogr2ogr

ToDo ...

space