Tag Archive | SRTM

Sverige mot oanade höjder

Detta inlägg kommer att bli en riktig ”höjdare”! Jag har en ”hög” ambitionsnivå! Inlägget kommer att nå oanade ”höjder”…

Nej nu räcker det med TV4 humor…

Nu skall det handla om höjddata och vad som framträder när man skalar bort allt annat.

När man tittar på kartor så är det oftast sjöar, skogar, städer, vägar och järnvägar som framträder. Det är sällan man sätter sig ner och tittar på en karta med enbart höjdinformation.

Screenshot from 2014-11-30 10:02:59Eftersom det numera går att hämta hem SRTM (radarinhämtade höjddata från rymdfärjan) med 1 sekunds upplösning, vilket är ca 30 meter stora pixlar, så finns det nu ganska högupplösta höjddata fritt tillgängliga över halva världen och resten är på väg (tidigare inlägg). I bilden ovan har jag kombinerat höjdinformation med ett lager med höjdskuggning, vilket jag tycker framhäver topografin på ett snyggt sätt.

Alla bilderna nedan kommer att vara ”nära” optimal pixelmatchning för vad som är möjligt med SRTM, medan bilden ovan mest skall ses som en översikt. Om du klickar på en bild så kommer den att visas i full upplösning på skärmen. En begränsning med SRTM för Sveriges del är att norr om Uppsala, eller 60° nordlig latitud, så tar det slut. Det genomfördes inga mätningar norr om detta.

Screenshot from 2014-11-30 10:16:16Var hittar vi ovanstående bild?

Nu kan du ju göra det här till en personlig utmaning så därför kommer här en förvarning. Om du skrollar ned lite så avslöjas platsen…

Lidköping! Det som stack ut var Högkullen till höger i bild, men även den väldigt ”vassa” udden till vänster får mig att fundera över de krafter som formade landet.

Screenshot from 2014-11-30 10:23:00Ny bild, ny utmaning.

Någonstans i denna bild så går det en järnvägstunnel rakt genom åsen. Det är nämligen Hallandsåsen och Båstad ligger i bukten precis norr om åsen.

Screenshot from 2014-11-30 10:28:33Den här platsen är kanske mest intressant av de tre hittills. Här ser man linjer i nord-sydriktning, men även skärningar i öst-väst. Hur detta gått till kan man också fundera över.

Var det är? Uppe till höger ligger Norrköping och strax under bilden i nedre vänstra kanten Linköping.

Screenshot from 2014-11-30 10:37:03Färgskalan i denna bild är lite blekare, vilket bör skvallra om att vi nu befinner oss på lite högre höjd. Man kan även se tydliga spår av människan med de linjer som framför allt är tydliga i nord-sydriktning något till höger om bildens mitt.

Det är E4 och Värnamo ligger något söderut. Till vänster om motorvägen ligger Store Mosse. Nationalparken är Sveriges största orörda myrmark söder om Lappland.

Jag nöjer mig så och kan bara uppmana till eget utforskande av höjddata. Om man inte är så intresserad av Sverige, så varför inte ta en tur till Grand Canyon, eller Rio De Janeiro?

Äh, vad sjutton vi avslutar med en bild över Grand Canyon i samma stil och skala som bilderna tidigare.

Screenshot from 2014-11-30 10:55:35

Du kan hämta data från:

http://earthexplorer.usgs.gov

Mer SRTM 1 släppt

Screenshot from 2014-11-24 12:38:44Nu finns ännu mera SRTM med 1 sekunds upplösning (ca 30 meter) att ladda hem fritt från exempelvis USGS. (Tidigare inlägg)

Tidigare var endast större delen av Afrika tillgängligt, men nu är även Nord- och Sydamerika samt Europa tillgänglia. Resten av Afrika och Asien kommer så fort det går.

”Mer data åt folket”

[EDIT]

Screenshot from 2014-11-24 13:24:43I bilden ovan är Jönköping uppe till vänster och Eksjö nere till höger. Kan du hitta Nässjö?

Höjddata är så bra att det tydligt framgår var vägarna går. Detta är främst möjligt då detta är ”reflekterad” höjd, vilket innebär att vegetationen är inräknad i höjdvärdena. Marken kan därför bli väldigt ”knottrig”, men detta är då förklaringen.

Klicka på bilden för att se den i bästa skärmupplösning.

Nya höjddata

För inte så länge sedan gav president Obama beskedet att SRTM data med 1 bågsekunds (ca 30 meter) upplösning skulle släppas fritt.

Inledningsvis så är det Afrika som släpps och sedan kommer resten av världen att följa.

Screenshot from 2014-10-05 15:03:20Täckningen just nu beskrivs av bilden ovan.

Då jag har ett professionellt intresse av en del regioner i Afrika så tänkte jag testa vad dessa data kan vara användbara till.

Jag loggar in med mitt USGS konto och väljer ut ett litet område som innebär att jag skall ladda hem fyra filer med höjddata. I det här fallet väljer jag GeoTIFF, men det går att välja BIL eller DTED också.

Då jag vill ha en enhetlig symbolsättning så använder jag virtuella raster (menyn Raster – Övrigt) för att läsa in de fyra filerna:

Screenshot from 2014-10-05 15:18:31Virtuella raster är mycket användbara i fall som detta och fungerar precis som ett vanligt rasterlager i allt väsentligt.

Screenshot from 2014-10-05 15:21:38Väl inläst i QGIS så ser det helt OK ut (bild ovan). Det framgår exempelvis väldigt tydligt var vattendragen finns.

Lite bearbetning då…

Screenshot from 2014-10-05 15:25:46Jag börjar med en hillshade. Observera z-faktorn i dialogen ovan, vilken behövs för att kompensera för att koordinater och höjddata är i olika format (jag använder 113 200 x cos (latitud) ).

Screenshot from 2014-10-05 15:29:50Utsnittet ovan är ca 50 km brett och visar tydligt på ett antal bergsformationer.

Genom att kombinera skuggningen med en satellitbild så blir det en ganska häftig effekt.

Screenshot from 2014-10-05 15:36:15Då skall vi se om det går att göra höjdkurvor också.

I bilden ovan så är höjdskillnaden lite drygt 200 meter, så jag provar med höjdkurvor med ekvidistans 20 meter till att börja med. Det kommer att ta någon minut för hela området, men jag har tid att vänta lite.

Screenshot from 2014-10-05 15:41:37Inte helt perfekt men för att vara obearbetat så ser det mycket bra ut (bild ovan).

Om jag zoomar in på en av formationerna och sätter ut etiketter på höjdkurvorna så ser det faktiskt väldigt bra ut.

Screenshot from 2014-10-05 15:45:42

Berget ovan är ca 10 kilometer från väst till öst och med tanke på storlek på skärm med mera så drar jag slutsatsen att denna typ av höjdkurvor kommer att bli utmärkta på kartor runt skala 1:50 000.

Screenshot from 2014-10-05 16:28:02Letar man lite, och man behöver inte anstränga sig särskilt mycket, så kan man hitta väldigt intressanta formationer i princip överallt.

Jag tror att dessa höjddata med skaplig upplösning kommer att bli väldigt användbara, inte minst inom olika hjälporganisationer runt om i världen, där det spelar roll om data kostar pengar.  Hela Afrika finns som sagt tillgängligt nu och resten av världen kommer så småningom. Planen är att allt skall vara tillgängligt senast under 2015.

Egna höjdkurvor.

Jag tänkte laborera lite med höjdkurvor som jag gör själv i detta inlägg.

Jag börjar med att ladda hem lite höjddata från SRTM i *.hgt format (länk). I den komprimerade filen man tar hem så finns det flera hgt-filer som när man lägger till dem symbolsätts individuellt. Ett enkelt sätt att skapa ett enda hanterbart lager med samma symbolsättning är att skapa ett virtuellt rasterlager.

Gå till ”Raster”, ”Övrigt” och ”Bygg Virtuellt Raster (Katalog)”.

Skärmbild från 2013-11-22 15:21:30

När jag försöker välja sökvägar med knapparna ”Välj…” i verktyget så totalkraschar QGIS, men det går bra att skriva eller klistra in sökvägarna till katalogen med hgt-filer och till den vrt-fil man vill skapa. Jag lägger vrt-filen på samma ställe som datat, jag tycker det blir enklast så just nu. Observera att jag även valt att öppna en katalog med filer och inte enskilda filer.

När verktyget kört så läggs den virtuella katalogen till i lagerlistan, därefter är det bara att klicka på ”Stäng” för att ta bort verktyget (varför kan det inte stängas automatiskt?).

SRTM

Bilden ovan visar ”tiles” till vänster och den virtuella katalogen till höger.

Sedan så kör man ”bara” verktyget för att skapa höjdkurvor.

Skärmbild från 2013-11-22 15:42:31

Du hittar verktyget under ”Raster”/”Extrahera”.

Jag kryssade även i rutan att skapa attributnamn för höjdkurvorna, så att dessa även fick ett höjdvärde. Kan vara bra att ha senare.

Höjdkurvorna gick ganska snabbt att få fram, med hänsyn till att de beräknades över hela Götaland och del av Svealand (i princip), och resultatet efter lite stilval och etiketter ser du här nedanför.

HöjdkurvorKlicka på bilden för att se den i full skärmupplösning.

Om jag får säga det själv så är jag ganska nöjd med resultatet. Jag använder en regelbaserad stil där de tunna linjerna är en ofiltrerad regel och de tjockare en regel med ”Regular Expression” som filtrerar fram rader på jämna 50-tal meter:

 regexp_match(  tostring(   toint( ”ELEV”) ),’.[05].’)

”ELEV” är mitt höjdattribut och för att testa attributet som text så tar jag bort decimaler med ”toint” och gör om talen till text med ”tostring”. Uttrycket matchar sedan strängar med ett valfritt tecken  på position ett och tre, och en femma eller nolla på position två. Uttrycket fungerar således bara på höjder under 1000 meter, men det är ganska lätt att modifiera.

Jag använder ett villkor med samma urval för att välja ut vilka linjer som skall få etiketter. Jag lägger även till ett ”m” efter höjdangivelsen.

CASE WHEN regexp_match( tostring( toint( ”ELEV”) ),’.[05].’) THEN tostring(ELEV) || ‘ m’ END

I bakgrunden så har jag den virtuella höjddatakatalogen som låg till grund för beräkningen. Och om du undrar så är det Vättern och Jönköping som syns på bilden.

Höjddata i SRTM

Många som håller på med GIS har hört talas om SRTM höjddata.

Detta är NASA ”Shuttle Radar Topography Mission” där man använt radar från rymdfärjan för att mäta upp hela jordytan.

Dessa finns att ladda hem från olika platser på nätet (exempelvis här eller här). Ett problem är att dessa ibland har ”hål” där data saknas. Men även att data norr om 60° ofta saknas.

Det finns även sidor där man kan ladda ner ”lagade” data. En av dessa sidor är View Finder Panoramas.

Där kan man ladda ner filer med upplösningen 1 och 3 bågsekunder över valda delar av världen (även norr om 60°). En bågsekund är ungefär 30 meter, och således är 3 bågsekunder omkring 100 meter. Det är ganska skaplig upplösning för att vara fritt tillgängligt på nätet.

Nedladdade filer i formatet hgt går att öppna och använda i exempelvis QGIS. Från andra sidor kan man även ladda hem data i exempelvis geotiff.