17.15. Rasterlagen clippen en samenvoegen¶

Voor deze les gaan we een hellingslaag berekenen voor een gebied dat een stedelijk gebied omsluit, wat is opgegeven in een vectorlaag met één enkele polygoon. De basis DEM is opgedeeld in twee rasterlagen die, tezamen, een veel groter gebied bedekken dan dat rondom de stad waarmee we willen werken. Als u het project opent dat overeenkomt met deze les, zult u iets zoals dit zien.

Deze lagen hebben twee problemen:

• Zij bedekken ene gebied dat veel te groot is voor wat we willen (we zijn geïnteresseerd in een kleinere regio rondom het centrum van de stad)

• Zij staan in twee verschillende bestanden. (De stadsgrenzen vallen precies binnen één enkele rasterlaag, maar, zoals gezegd, we willen er wat extra gebied rondom bij).

Beide zijn eenvoudig op te lossen met behulp van de toepasselijke geo-algoritmen.

Eerst willen we een rechthoek maken die het gebied definieert dat we willen. We maken, om dit te doen, een laag die het begrenzingsvak van de laag met de grenzen van het stedelijk gebied bevat, en dan bufferen we die, om een rasterlaag te verkrijgen die iets meer bedekt dan strikt noodzakelijk.

We kunnen het algoritme Polygon from layer extent gebruiken om het begrenzingsvak te berekenen.

We gebruiken het algoritme Fixed distance buffer, met de volgende waarden voor de parameters om het te bufferen.

Hier is het resulterende begrenzingsvak dat is verkregen met behulp van de hierboven weergegeven parameters

Het is een afgerond vak, maar we kunnen eenvoudig het equivalente vak met vierkante hoeken verkrijgen door het algoritme Polygon from layer extent er op uit te voeren. We zouden eerst de stadsgrenzen hebben kunnen bufferen, en dan de rechthoek van het bereik berekenen, wat ons één stap bespaart.

Het zal u opvallen dat de rasters een verschillende projectie dan de vector hebben. We zouden ze daarom opnieuw moeten projecteren, met behulp van het gereedschap Warp (Herprojecteren), voordat we verder gaan.

met deze laag die het begrenzingsvak bevat van de rasterlaag die we willen verkrijgen, kunnen we beiden rasterlagen bijsnijden, met behulp van het algoritme Clip Grid with Polygons algorithm.

Als de lagen eenmaal zijn bijgesneden, kunnen zij worden samengevoegd met behulp van het algoritme Merge raster layers.

Een celgrootte is nodig voor de samengevoegde laag. We zullen dezelfde gebruiken als voor de invoer. U hoeft dat niet vooraf te weten voordat u het algoritme aanroept. Klik eenvoudigweg op de knoop aan de rechterkant van het tekstveld en u zult een dialoogvenster krijgen om eenvoudige rekenkundige formules in te voeren, en een lijst van frequent gebruikte waarden, met daarnaast de celgrootten en coördinaten voor de begrenzing van alle beschikbare lagen.

Opmerking: U kunt tijd besparen door eerst samen te voegen en dan bij te snijden, en u zult voorkomen dat u het algoritme voor clippen tweemaal moet aanroepen. Als er echter meerdere lagen zijn die moeten worden samengevoegd en zij zijn van een nogal aanzienlijke grootte, dan zult u eindigen met een grote laag die later misschien moeilijk te verwerken is. In dat geval dient u misschien het algoritme voor clippen meerder malen aan te roepen, wat tijd zou kunnen vergen, maar maak u geen zorgen, we zullen snel zien dat er enkele aanvullende gereedschappen zijn om die bewerking te automatiseren. In dit voorbeeld hebben we slechts twee lagen, dus u zou zich daar nu geen zorgen over moeten maken.

Daarmee krijgen we de uiteindelijke DEM die we willen.

Nu is het tijd om de hellingslaag te berekenen.

Een hellingslaag kan worden berekend met het algoritme Slope, Aspect, Curvature, maar de in de laatste stap verkregen DEM is niet geschikt als invoer, omdat de waarden voor de hoogte in meters zijn, maar de celgrootte is niet uitgedrukt in meters (de laag gebruikt een CRS met geografische coördinaten). Opnieuw projecteren is nodig. Het algoritme Warp (herprojecteren) kan opnieuw worden gebruikt om een rasterlaag opnieuw te projecteren. We projecteren opnieuw naar een CRS met meters als eenheden (bijv. 3857), zodat we daarna de helling correct kunen berekenen, met ofwel SAGA of GDAL.

Hier is de opnieuw geprojecteerde DEM.

Met de nieuwe DEM kan de helling nu worden berekend.

En hier is de uiteindelijke hellingslaag.

De door het algoritme Slope, Aspect, Curvature geproduceerde helling is uitgedrukt in radialen, maar graden zijn een meer praktische en algemenere eenheid. Het algoritme Metric conversions zal ons helpen de conversie uit te voeren (maar voor het geval u niet wist dat dat algoritme bestond, kunt u de Raster calculator gebruiken die we al eerder hebben gebruikt).

We krijgen de uiteindelijke laag die we willen door de geconverteerde hellingslaag opnieuw terug te projecteren met Reproject raster layer.

De processen voor opnieuw projecteren hebben er voor gezorgd dat de uiteindelijke laag gegevens bevat die buiten het begrenzingsvak liggen dat we in een van de eerste stappen hebben berekend. Dit kan worden opgelost door het opnieuw te clippen, zoals we hebben gedaan om de basis DEM te krijgen.