Belangrijk

Vertalen is een inspanning van de gemeenschap waaraan u deel kunt nemen. Deze pagina is momenteel voor 66.67% vertaald.

17.15. Rasterlagen clippen en samenvoegen

Notitie

In deze les zullen we een ander voorbeeld van ruimtelijke voorbereiding zien, om het gebruiken van geo-algoritmes 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 een 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-algoritmes.

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.

To calculate the bounding box , we can use the Extract layer extent algorithm

To buffer it, we use the Buffer algorithm, with the following parameter values.

Waarschuwing

Syntaxis is in recente versies gewijzigd; stel zowel Afstand als Boogpunt in op .25

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

It is a rounded box, but we can easily get the equivalent box with square angles, by running the Extract layer extent algorithm on it. We could have buffered the city limits first, and then calculate the extent rectangle, saving one step.

You will notice that the rasters has a different projection from the vector. We should therefore reproject them before proceeding further, using the Warp (reproject) tool.

Notitie

Recente versies hebben een meer complexe interface. Zorg er voor dat tenminste één methode voor compressie is geselecteerd.

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 raster by mask layer algorithm.

Once the layers have been cropped, they can be merged using the SAGA Mosaic raster layers algorithm.

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.

A slope layer can be computed with the Slope algorithm, but the DEM obtained in the last step is not suitable as input, since elevation values are in meters but cellsize is not expressed in meters (the layer uses a CRS with geographic coordinates). A reprojection is needed. To reproject a raster layer, the Warp (reproject) algorithm can be used again. We reproject into a CRS with meters as units (e.g. 3857), so we can then correctly calculate the slope, with either SAGA or GDAL.

Met de nieuwe DEM kan de helling nu worden berekend.

En hier is de uiteindelijke hellingslaag.

The slope produced by the Slope algorithm can be expressed in degrees or radians; degrees are a more practical and common unit. In case you calculated it in radians, 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).

Reprojecting the converted slope layer back with the Reproject raster layer, we get the final layer we wanted.

De processen voor opnieuw projecteren zouden er voor hebben kunnen zorgen 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.