Igår påbörjade jag ett arbete med ”fjärranalys” av Landsat 8 data. Detta inlägg är en direkt fortsättning på det inlägget och nu är det dags att se vad resultatet blev.
Har du inte läst gårdagens inlägg så skall du nog börja med det, annars så kör vi på…
Det sista jag gjorde i går var att köra igång en vektorisering av det rasterlager som var resultatet av SCP. Detta tog riktigt lång tid men det var många ganska små områden och en riktigt stor yta att beräkna.
Men då QGIS i mitt tycke är mycket bättre på att göra många saker samtidigt än ArcGIS så kunde jag fortsätta att jobba med andra lager medan beräkningen höll på. Dialogrutan i bilden nedan flyter ovanpå QGIS men hindrar mig inte att fortsätta jobba i övrigt.
Den pan-skärpta Landsat bilden är en jättebra referens när man kontrollerar resultatet och lite transparens på olika lager från exempelvis Open Street Map gör att det ser ganska bra ut redan.
Notera att jag flyttat från det område jag använda för att ”träna” min fjärranalys i bilden ovan. Detta gör jag för att kontrollera att inställningarna fungerar över hela ytan. Det är nog ändå klokt att begränsa sig till mindre områden när man gör fjärranalyser, inte minst om man vill behöva rätta så lite som möjligt i de data som genereras.
Av någon anledning så blev mitt tänkta polygonlager ett linjelager…
Inte hela världen då det enkelt går att göra om linjelager till polygoner med geometriverktygen.
Jag vill även förenkla dessa vektordata då de är lite väl ”pixeliserade”.
En förenkling med ”20 meter” reducerar antalet brytpunkter rejält, men ger också upphov till så kallade ”slivers” eller glipor mellan ytorna.
Detta går att fixa men just nu så duger det för mina experiment.
Om jag enbart ställer in färger för mitt nya lager, lägger till vägar och järnvägar från Open Street Map och etiketter från Lantmäteriets Översiktskarta så blir det som i bilden ovan. Egentligen fantastiskt bra redan.
Lite justeringar på exempelvis stilarna för vägar och en mycket lätt transparens för vissa markytor så att Landsat bilden syns igenom något, samt ytterligare några utvalda objekt från Open Street Map så är det klart.
Jämför jag med OSM i original, och då framför allt terrängen i området så är kanske inte Halmstad det bästa exemplet att jämföra med eftersom OSM redan är ganska bra där. För Eksjö så är dock skillnaden stor.
Bilden ovan är inte helt rättvisande då jag även lagt på de höjdlager jag jobbade med förra veckan, baserade på Lantmäteriets öppna höjddata, men jämför man enbart terrängen utanför vägar så kan det här vara ett sätt att gå vidare för att förbättra Open Street Map över stora ytor, och nej då menar jag inte att publicera data baserat på Lantmäteriets produkter utan enbart sådant som är bearbetat med fjärranalys i QGIS från Landsat 8 data, vilka så långt jag kunnat läsa mig till är ”public domain”.