17.15. Rasterlagen clippen en samenvoegen¶
Notitie
In deze les zullen we een ander voorbeeld van ruimtelijke voorbereiding zien, om het gebruiken van geo-algoritmen in scenario’s van de echte wereld te vervolgen.
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.
Notitie
Recent versions have a more complex interface. Select at least one compression method.
With this layer that contains the bounding box of the raster layer that we want to obtain, we can crop both of the raster layers, using the Clip grid with polygon algorithm.
Als de lagen eenmaal zijn bijgesneden, kunnen zij worden samengevoegd met behulp van het algoritme van GDAL Merge.
Waarschuwing
The following paragraph is obsolete.
A cellsize is needed for the merged layer. We will use the same one of the input ones. You do not need to know it in advance before calling the algorithm. Just click on the button in the right–hand size of the text field and you will have a dialog to enter small mathematical formulas, and a list of frequently used values, among them the cellsizes and bounding coordinates of all available layers.
Notitie
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.
Here is the reprojected DEM.
Met de nieuwe DEM kan de helling nu worden berekend.
En hier is de uiteindelijke hellingslaag.
The slope produced by the Slope, Aspect, Curvature algorithm is expressed in radians, but degrees are a more practical and common unit. The Metric conversions algorithm will help us to do the conversion (but in case you didn’t know that algorithm existed, you could use the raster calculator that we have already used).
We krijgen de uiteindelijke laag die we willen door de geconverteerde hellingslaag opnieuw terug te projecteren met Reproject raster layer.
Waarschuwing
todo: Afbeelding toevoegen
The reprojection processes have caused the final layer to contain data outside the bounding box that we calculated in one of the first steps. This can be solved by clipping it again, as we did to obtain the base DEM.