Viktigt
Översättning är en gemenskapsinsats du kan gå med i. Den här sidan är för närvarande översatt till 100.00%.
17.15. Klippning och sammanslagning av rasterlager
Observera
I den här lektionen kommer vi att se ett annat exempel på förberedelse av rumsliga data för att fortsätta använda geoalgoritmer i verkliga scenarier.
I den här lektionen ska vi beräkna ett lutningslager för ett område som omger en stad, vilket ges i ett vektorlager med en enda polygon. DEM-basen är uppdelad i två rasterlager som tillsammans täcker ett område som är mycket större än området runt den stad som vi vill arbeta med. Om du öppnar projektet som motsvarar den här lektionen kommer du att se ungefär så här.
Dessa lager har två problem:
De täcker ett område som är för stort för det vi vill ha (vi är intresserade av en mindre region runt stadens centrum)
De ligger i två olika filer (stadsgränsen ligger i ett enda rasterlager, men som sagt, vi vill ha lite extra yta runtomkring).
Båda dessa problem kan enkelt lösas med lämpliga geoalgoritmer.
Först skapar vi en rektangel som definierar det område vi vill ha. För att göra det skapar vi ett lager som innehåller begränsningsrutan för lagret med gränserna för stadsområdet, och sedan buffrar vi det så att vi får ett rasterlager som täcker lite mer än vad som är absolut nödvändigt.
För att beräkna begränsningsrutan kan vi använda algoritmen Polygon från lagerutbredning
För att buffra den använder vi algoritmen Fixed distance buffer med följande parametervärden.
Varning
Syntaxen har ändrats i de senaste versionerna; ställ in både Distance och Arc vertex till .25
Här är den resulterande begränsningsrutan som erhålls med hjälp av parametrarna som visas ovan
Det är en rundad låda, men vi kan enkelt få en motsvarande låda med fyrkantiga vinklar genom att köra algoritmen Polygon från lagrets utsträckning på den. Vi kunde ha buffrat stadsgränserna först och sedan beräknat omfattningsrektangeln, vilket sparar ett steg.
Du kommer att märka att rasterna har en annan projektion än vektorn. Vi bör därför omprojicera dem innan vi går vidare, med hjälp av verktyget Warp (reproject).
Observera
De senaste versionerna har ett mer komplext gränssnitt. Se till att minst en komprimeringsmetod är vald.
Med detta lager som innehåller gränsboxen för det rasterlager som vi vill få fram kan vi beskära båda rasterlagren med hjälp av algoritmen Clip raster with polygon.
När lagren har beskurits kan de slås samman med hjälp av SAGA Mosaic raster layers-algoritmen.
Observera
Du kan spara tid genom att först sammanfoga och sedan beskära, och du undviker att anropa beskärningsalgoritmen två gånger. Men om det finns flera lager att slå samman och de har en ganska stor storlek, kommer du att få ett stort lager som det senare kan vara svårt att bearbeta. I det fallet kan du behöva anropa klippalgoritmen flera gånger, vilket kan vara tidskrävande, men oroa dig inte, vi kommer snart att se att det finns några ytterligare verktyg för att automatisera den operationen. I det här exemplet har vi bara två lager, så du behöver inte oroa dig för det nu.
Därmed får vi den slutliga DEM som vi vill ha.
Nu är det dags att beräkna lutningslageret.
Ett lutningslager kan beräknas med algoritmen Slope, Aspect, Curvature, men den DEM som erhölls i det sista steget är inte lämplig som indata, eftersom höjdvärdena är i meter men cellstorleken inte uttrycks i meter (lageret använder ett CRS med geografiska koordinater). En omprojektering behövs. För att återprojicera ett rasterlager kan algoritmen Warp (reproject) användas igen. Vi omprojicerar till ett CRS med meter som enhet (t.ex. 3857), så att vi sedan kan beräkna lutningen korrekt, med antingen SAGA eller GDAL.
Med den nya DEM kan lutningen nu beräknas.
Och här är det resulterande sluttningslageret.
Lutningen som produceras av algoritmen Slope, Aspect, Curvature kan uttryckas i grader eller radianer; grader är en mer praktisk och vanlig enhet. Om du beräknade den i radianer hjälper algoritmen Metric conversions oss att göra omvandlingen (men om du inte visste att den algoritmen fanns kan du använda rasterkalkylatorn som vi redan har använt).
Genom att återprojicera det konverterade sluttningslageret tillbaka med * Återprojicera rasterlager * får vi det slutliga lageret vi ville ha.
Varning
att göra: Lägg till bild
Reprojektionsprocesserna kan ha lett till att det slutliga lagret innehåller data utanför den avgränsande box som vi beräknade i ett av de första stegen. Detta kan lösas genom att klippa det igen, precis som vi gjorde för att få fram bas-DEM.