Tutorial: Introductie GDAL

De laatste oefening uit de vorige sectie zou efficiënter zijn geweest als je alles in één keer had kunnen uitvoeren, zonder het tussentijds wegschrijven van de datatype‑conversie. Voor dat soort workflows bestaat het commando gdal raster pipeline.

Laten we het toepassen op de volledige workflow, die de volgende stappen uitvoert zonder tussentijdse bestanden op schijf te schrijven:

Eerst moeten we de bestaande tussentijdse resultaten verwijderen om zeker te zijn dat we nieuwe data aanmaken.

1. Voer het volgende commando uit om alle bestandsnamen met dem erin te verwijderen:
   del dem*.*

2. Nu kunnen we de pipeline typen:

gdal raster pipeline read srtm_37_02.tif ! reproject -d EPSG:28992 ! set-type --ot Float32 ! write dem_rd.map -f PCRaster --co VS=SCALAR

Uitleg:

• Het commando om een raster‑pipeline te maken is gdal raster pipeline.
• We moeten beginnen met het inlezen van het invoerraster: read srtm_37_02.tif.
• Elke volgende bewerkingsstap moet worden voorafgegaan door !.
• Eerst herprojecteren we: reproject -d EPSG:28992.
• Daarna zetten we het datatype van het resultaat van die stap naar Float32: ! set-type --ot Float32.
• Vervolgens schrijven we het resultaat weg naar PCRaster‑formaat met het datatype scalar: ! write dem_rd.map -f PCRaster --co VS=SCALAR.

Concluderend: pipelines zijn een eenvoudige manier om workflows te maken die je in één keer kunt uitvoeren. In combinatie met batchbestanden kan dit heel gebruiksvriendelijk en krachtig zijn, zonder dat je scripts in Python hoeft te schrijven.