17.2. Rasteranalyse

17.2.1. Rasterberekeningen

Rasterberekeningen, in het menu Raster, stelt u in staat berekeningen uit te voeren op basis van bestaande raster pixelwaarden (zie Fig. 17.21). De resultaten worden weggeschreven naar een nieuw raster, in een indeling die wordt ondersteund door GDAL.

../../../_images/raster_calculator1.png

Fig. 17.21 Rasterberekeningen

De lijst Rasterbanden bevat alle geladen rasterlagen die kunnen worden gebruikt. Dubbelklik, om een raster aan het expressieveld van Rasterberekeningen toe te voegen, op de naam in de lijst Velden. U kunt de operatoren gebruiken om expressies voor berekeningen te construeren, of u kunt ze eenvoudigweg in het vak typen.

In het gedeelte Resultaatlaag moet u een uitvoerlaag definiëren. U kunt:

  • checkbox Direct raster maken in plaats van de laag naar schijf te schrijven:

    • Indien niet geselecteerd wordt de uitvoer opgeslagen op schijf als een nieuw gewoon bestand. Een pad voor Uitvoerlaag en een Indeling uitvoer zijn vereist.

    • Indien geselecteerd wordt een virtuele rasterlaag, d.i. een rasterlaag die wordt gedefinieerd door zijn URI en waarvan de pixels direct worden berekend, gemaakt. Het is geen nieuw bestand op schijf; de virtuele laag is nog steeds verbonden met de rasters die zijn gebruikt voor de berekeningen, wat betekent dat het verwijderen of verplaatsen van deze rasters het zouden verbreken. Een Laagnaam kan worden opgegeven, anders wordt de expressie voor de berekening als zodanig gebruikt. Verwijderen van de virtuele laag uit het project verwijdert het, en het kan persistent als bestand worden gemaakt door het contextmenu van de laag Exporteren ► Opslaan als… te gebruiken.

  • Definieer het Ruimtelijk bereik van de berekening, gebaseerd op een bereik van een invoer rasterlaag, of op aangepaste coördinaten X,Y

  • Stel de Resolutie voor de laag in met kolom- en rijnummer. Als de invoerlaag een andere resolutie heeft, zullen de waarden worden geresampled met het algoritme Nearest neighbor.

  • Met het keuzevak checkbox Voeg resultaat toe aan project,zal de resultaatlaag automatisch worden toegevoegd aan het gebied van de Legenda en kan worden gevisualiseerd. Standaard geselecteerd voor virtuele rasters.

Het gedeelte Operatoren bevat alle beschikbare operatoren. Klik, om een operator toe te voegen aan het expressievak van Rasterberekeningen, op de toepasselijk knop. Rekenkundige berekeningen (+, -, *, … ) en trigonometrische functies (sin, cos, tan, … ) zijn beschikbaar. Voorwaardelijke expressies (=, !=, <, >=, … ) geven ofwel 0 voor false of 1 voor true terug, en kunnen daarom worden gebruikt met andere operatoren en functies.

Hint

Bekijk ook het algoritme Raster calculator.

17.2.1.1. Voorbeelden

Hoogtewaarden van meter naar voet omzetten

Voor het omzetten van een hoogteraster van meters naar voet kunt u de omrekeningsfactor 3.28 gebruiken. De expressie is:

"elevation@1" * 3.28

Een uitknipmasker gebruiken

Wanneer u van een hoogterasterkaart dat deel wilt uitknippen dat boven 0 meter hoogte ligt – kunt u de volgende expressie gebruiken om in één keer een uitknipmasker te maken en het resultaat weg te schrijven naar een nieuwe rasterkaart.

("elevation@1" >= 0) * "elevation@1"

Met andere woorden: voor elke cel met een waarde groter dan of gelijk aan 0 evalueert de voorwaardelijke expressie tot 1, wat de originele waarde behoudt door die te vermenigvuldigen met 1. Anders evalueert de voorwaardelijke expressie tot 0, wat de rasterwaarde instelt op 0. Dit maakt direct het uitknipmasker.

Een raster classificeren

Als u een raster wilt classificeren – bijvoorbeeld in twee hoogteklassen, kunt u de volgende expressie gebruiken om, in één stap, een raster te maken met twee waarden 1 en 2.

("elevation@1" < 50) * 1 + ("elevation@1" >= 50) * 2

Met andere woorden: voor elke cel kleiner dan 50 wordt zijn waarde ingesteld op 1. Voor elke cel groter dan of gelijk aan 50 wordt de waarde ingesteld op 2.

Of u kunt de operator IF gebruiken.

if ( elevation@1 < 50 , 1 , 2 )

17.2.2. Raster uitlijnen

Dit gereedschap is in staat om meerdere rasters als invoer te gebruiken en ze perfect uit te lijnen, wat betekent:

  • opnieuw projecteren naar hetzelfde CRS,

  • opnieuw samplen naar dezelfde celgrootte en verschuiving in het raster,

  • een interessegebied te clippen,

  • waarden, indien nodig, opnieuw op schaal brengen.

Alle rasters zullen worden opgeslagen in andere bestanden.

Open eerst de gereedschappen via Raster ► Rasters uitlijnen… en klik op de knop symbologyAdd Nieuw raster toevoegen om een bestaand raster in QGIS te selecteren. Selecteer een uitvoerbestand om het raster na het uitlijnen op te slaan, de methode voor opnieuw samplen en of de gereedschappen Waarden opnieuw schalen naar celgrootte moeten uitvoeren. De methode voor resamplen kan zijn (zie Fig. 17.22):

  • Nearest Neighbor

  • Bilineair (2x2 kernel)

  • Kubisch (4x4 kernel): Kubische convolutie benadering

  • Kubisch B-Spline (4x4 kernel): Kubische B-Spline benadering

  • Lanczos (6x6 kernel): Lanczos windowed sinc interpolatie

  • Gemiddelde: berekent het gemiddelde van alle bijdragende niet-NODATA pixels

  • Modus: selecteert de waarde die het meest voorkomt van alle monsterpunten

  • Maximum, Minimum, Mediaan, Eerste kwartiel (Q1) of Derde kwartiel (Q3) van alle bijdragende niet-NODATA pixels

../../../_images/raster_align_edit.png

Fig. 17.22 Opties voor resamplen van raster selecteren

In het hoofddialoogvenster Raster uitlijnen kunt u nog steeds symbologyEdit Bestandsinstellingen bewerken of symbologyRemove Bestaand bestand verwijderen uit de lijst met rasterlagen. U kunt ook een of meer andere opties kiezen (zie Fig. 17.23):

  • De Referentielaag selecteren,

  • Naar een nieuw CRS transformeren,

  • Een andere Celgrootte instellen,

  • Een andere Ruitennetverspringing instellen,

  • Tot bereik clippen: kan gebruikergedefinieerd zijn of gebonden aan een laag of aan het kaartvenster

  • Uitvoergrootte

  •  Uitgelijnde rasters toevoegen aan kaartvenster.

../../../_images/raster_align.png

Fig. 17.23 Raster uitlijnen