Laserdata

1(12)
D atum: 2015-10-01
D ok umentversion: 2.2
Produktbeskrivning:
Laserdata
LANTMÄTERIET
2015-10-01
2 (12)
Innehållsförteckning
1
Allmän beskrivning ..................................................................... 3
1.1
Innehåll ................................................................................................ 3
1.2
Geografisk täckning .............................................................................. 3
1.3
Geografiskt utsnitt ................................................................................ 4
1.4
Koordinatsystem ................................................................................... 4
2
Kvalitetsbeskrivning ................................................................... 4
3
Leveransens innehåll .................................................................. 5
3.1
Katalogstruktur i leverans .................................................................... 5
3.2
Leveransformat .................................................................................... 6
3.3
Filuppsättning och innehåll ................................................................... 6
3.3.1
Beskrivning av bildfiler i metadata ............................................................. 7
3.3.2
Beskrivning av innehåll i XML-filen ............................................................ 10
4
Förändringsförteckning ............................................................ 12
4.1
Senaste förändring ............................................................................. 12
4.2
Tidigare förändringar .......................................................................... 12
Bilaga A – Produktions- och skanningsområden
Bilaga B – Kvalitetsbeskrivning nationell höjdmodell
LANTMÄTERIET
2015-10-01
3 (12)
1 Allmän beskrivning
Lantmäteriet har sedan 2009 på regeringens uppdrag arbetat med att framställa en
nationell höjdmodell med hög noggrannhet. Laserskanning valdes som metod för insamlingen. Vid den efterföljande bearbetningen har laserpunkterna klassificerats
som punkter på mark, vatten eller broar. Övriga punkter har lämnats oklassificerade.
1.1 Innehåll
Produkten Laserdata utgörs av ett punktmoln där varje punkt är klassificerad till
mark, vatten, bro eller oklassificerad (vilket även inkluderar eventuella felaktiga
punkter). Till produkten levereras metadata som redovisar ursprung och status för
genomförd bearbetning, samt fullständiga metadata från skanningen, inklusive flygstråk.
För den användare som behöver kan även filer med överlappande data, ”kantfiler”,
beställas. Kantfiler innehåller data i en överlappszon på ca 200 meter in på angränsande skanningsområde. Om dessa filer önskas måste detta särskilt anges vid beställning.
Vad som ingår i datasetet framgår under rubriken 3.3 Filuppsättning och innehåll.
1.2 Geografisk täckning
Med laserdata som grund skapas en rikstäckande terrängmodell (en avbildning av
markens form). Framväxten av laserskanningen redovisas på Lantmäteriets hemsida,
se www.geolex.lm.se.
·
Leverantörens veckorapport - I en fil redovisas mer detaljerad information
om skanningen, med tidpunkter för olika delmoment, utrustning och antal
skannade stråk.
·
Produktionsstatus och klart i lager – Det förstnämnda redovisas för att ge
användare möjlighet till fältinventering i nära anslutning till att skanning genomförs. Översikten redovisar områden i tre olika statusnivåer; när stråkplanering är godkänd, skanning påbörjad och skanning preliminärt avslutad
(omskanning kan komma ifråga om kvaliteten underkänns i den efterföljande
kontrollen). Efter att skanningen är avslutad kan det dröja ytterligare upp till
ca 6 månader innan laserdata finns klart i lager.
Klart i lager redovisar områden, med klassificeringsnivå på respektive skanningsområde, som är lagrade i grunddatalagret och klara för leverans till användare.
Sverige är indelat i produktions- och skanningsområden. Se bilaga A.
LANTMÄTERIET
2015-10-01
4 (12)
1.3 Geografiskt utsnitt
Minsta enhet för bearbetning och leverans motsvarar en ruta om 2,5 x 2,5 km, en
s.k. bearbetningsruta, anpassad geografiskt till bladindelningen i indexsystemet i
SWEREF 99 TM. Rutan har en egen bladbeteckning med nedre vänstra hörnets koordinat i hundratal meter åtföljd av sidans längd i hundratal meter, t.ex. 67275_5975_25.
1.4 Koordinatsystem
Plan: SWEREF 99 TM
Höjd: RH 2000
2 Kvalitetsbeskrivning
Se bilaga B Kvalitetsbeskrivning nationell höjdmodell.
LANTMÄTERIET
2015-10-01
5 (12)
3 Leveransens innehåll
3.1 Katalogstruktur i leverans
I en leverans sorteras lasfilerna inklusive metadata in under en katalog med samma
namn som skanningsområdet de tillhör, för att hålla samman metadata som är gemensamma.
Nedanstående exempel visar katalogstrukturen. Metadata för hela skanningsområdet ligger i en underkatalog till skanningsområdet.
LANTMÄTERIET
2015-10-01
6 (12)
3.2 Leveransformat
Laserdata levereras i LAS-format version 1.2, i punktdataformat 1 med intensitet och
veckotid. Lasfilerna, och eventuellt beställda kantfiler, levereras packade till RARformat (programvara för att packa upp filer finns att ladda ner gratis från internet),
vilket innebär att filstorleken minskas till ca 1/5 av ursprunglig storlek.
3.3 Filuppsättning och innehåll
Filnamn (exempel)
Beskrivning
09P001_67475_5875_25.las
(levereras packat i .rar)
I filnamnet ingår identiteten för skanningsområdet, koordinaterna för rutans nedre vänstra hörn, storleken på rutan i 100tal meter samt filformat.
09P001_67475_5875_25_density.tif
Rasterfil med punkttäthet för de laserpunkter som klassificerats som mark, en för varje 2,5 km-ruta. Se beskrivning 3.3.1.2.
09P001_67475_5875_25_density.tfw
Georefereringsfil till bildfilen ovan.
09P001_las_21850.xml
(skanningsområde_las_ID-nr.xml)
Metadata som redovisar ursprung och grad av bearbetning (se
avsnitt 3.3.2).
09P001_coverage.tif
En rasterfil som redovisar punkttäthet i laserpunktmolnet,
sista och enda eko, för hela skanningsområdet. Se beskrivning
3.3.1.1.
09P001_coverage.tfw
Georefereringsfil till bildfilen ovan.
09P001_density.tif
En rasterfil för hela skanningsområdet som visar punkttäthet
för de laserpunkter som klassificerats som mark.
09P001_density.tfw
Georefereringsfil till bildfilen ovan.
09P001_tasq_bild.tif
Redovisar höjdavvikelse i överlappet mellan stråken inom
skanningsområdet. Se beskrivning 3.3.1.3.
09P001_tasq_bild.tfw
Georefereringsfil till bildfilen ovan.
561444_561794.trj
Bandata (trajectory data), en fil för varje stråk ingående i aktuellt skanningsområde. Filnamnet kan variera beroende på vilken utrustning som använts, men slutar alltid på .trj. I filen
med stråkmetadata, finns namnet för respektive trajectoryfil
kopplad till stråk.
Skanningenheter.txt
Alla skanningsenheter som har använts vid skanning av landet. Inte bara de som berör aktuell leverans.
21850.kml (ID-nr.kml)
Fil som visar områdets utbredning i Google Earth, i SWEREF
99 TM.
09P001_67475_5875_25.las
(levereras packat i .rar)
Kantfiler. Innehåller laserdata i en överlappszon på ca 200 meter in på angränsande skanningsområde. Filerna har samma
namn som ordinarie lasfiler, men ligger i en egen katalog med
namn enligt ”09P001_kant”.
OBS! Filerna måste beställas särskilt!
LANTMÄTERIET
2015-10-01
7 (12)
3.3.1 Beskrivning av bildfiler i metadata
3.3.1.1 Punkttäthet - sista och enda eko
För varje skanningsområde redovisas punkttätheten i laserpunktmolnet som en bildfil, med 10 m upplösning, benämnd t.ex. 09P001_coverage.tif. Endast punkter från
sista av flera eller enda eko redovisas, det vill säga de punkter som kan representera
markytan.
Punkttätheten i laserpunktmolnet åskådliggörs med färger enligt tabellen nedan:
Färg
Punkttäthet
Blått
> 1 punkt/m2
Grönt
> 0.5 punkter/m2
Gult
> 0.25 punkter/m2
Rött
< 0.25 punkter/m2
Svart
0 punkter/m2
Kommentar
Svart färg i bilden beror dels på att vattenytorna har maskats bort, dels
på hål i laserpunktmolnet p.g.a. dålig reflektion. Även del av skanningsområde som inte skannats för att det är hav eller utanför riksgränsen, redovisas som svart.
Bilden redovisar punkttäthet i laserpunktmolnet, sista samt enda eko, för hela skanningsområdet.
LANTMÄTERIET
2015-10-01
8 (12)
3.3.1.2 Punkttäthet - markklassificerade laserpunkter
En bild med 10 m upplösning redovisar den genomsnittliga punkttätheten av laserpunkter som klassificerats som mark. Varje fil täcker en bearbetningsruta om
2,5 x 2,5 km. Punkttätheten åskådliggörs med färger enligt tabellen nedan. Filen heter
t.ex. 09P001_67475_5875_25_density.tif. Motsvarande rasterbild levereras även för
hela skanningsområdet.
Färg
Punkttäthet
Kommentar
Blått
> 0.5 pkt/m2
På öppna ytor och i överlappen mellan stråk kan det bli fler markträffar
än minimikravet på 0.5 pkt/m2.
Grönt
0.25-0.5 pkt/m2
I genomsnitt finns det minst en markträff inom ett område om 2 × 2 m.
Gult
0.0625-0.25 pkt/m2
I genomsnitt finns minst en markträff inom ett område om 4 × 4 m. Terrängmodellen kan ha försämrad detaljeringsgrad.
Rött
< 0.0625 pkt/m2
I genomsnitt finns det mindre än en markträff inom ett område om
4 × 4 m. Orsaken kan exempelvis vara tät skog eller vatten. Terrängmodellen kan ha kraftigt försämrad detaljeringsgrad.
Svart
0 pkt/m2
Svart färg i bilden beror dels på att vattenytorna har maskats bort, dels
på hål i laserpunktmolnet. Hål i laserpunktmolnet beror på dålig reflektion eller tät vegetation, vilket kan orsaka fullständigt bortfall av markträffar. Dålig reflektion kan förekomma på t.ex. vattenytor och nylagd
asfalt.
Exempel: Färgerna i densitetsbilden representerar olika punkttäthet på mark i laserdata.
LANTMÄTERIET
2015-10-01
9 (12)
3.3.1.3 Höjdnoggrannhet i stråköverlapp
Statistik som visar avvikelser i höjd mellan flygstråken inom varje skanningsområde
redovisas i en bildfil med 10 m upplösning. Filens namn är t.ex. 09P001_tasq_bild.tif.
Avvikelser i höjd mellan flygstråken. Blå färg representerar <0.10 m, grön 0.10-0.25 m,
gul 0.25-0.50 m och röd >0.50 m avvikelse mellan stråken.
LANTMÄTERIET
2015-10-01
10 (12)
3.3.2 Beskrivning av innehåll i XML-filen
Följande är ett exempel på innehåll i metadatafilen skanningsområde_las_ID-nr.xml.
XML-schema finns att hämta här: http://namespace.lantmateriet.se/gdshojd/.
I första delen finns information som är gemensam för hela skanningsområdet.
Fält
Exempel
Förklaring
Områdesnamn
09P001
Identitet på skanningsområde.
Ursprung
1
1 = Lantmäteriets laserskanning 2009-2016.
Höjdnoggrannhet
0.05 m
Utfall från kontroll mot kända punkter.
Punkterna ligger på öppna plana hårdgjorda
ytor.
Kontrollytor i höjd
9
Antal kontrollytor i höjd.
Kontrollytor i plan
7
Antal kontrollytor i plan.
Kvalitetsanmärkning
Stor mängd felaktigt
markklassificerad låg
vegetation.
Anger att man under bearbetningen upptäckt avvikelser eller problem som användaren bör vara uppmärksam på.
Klassificeringsprogramvara
TerraScan 009.006
Programvara som använts för att klassificera
alla rutor inom skanningsområdet till klassificeringsnivå 1. Annan programvara eller
version kan förekomma för enstaka rutor
med en högre klassificeringsnivå.
Ursprungligt klassificeringsdatum
2009-11-20
Datum för den ursprungliga klassificeringen.
Därefter följer information om respektive levererad ruta som ingår i skanningsområdet:
Fält
Exempel
Förklaring
Ruta
67475_5875_25
Koordinatangivelse för rutans nedre vänstra
hörn och utbredning på marken i 100-tal meter.
Skanningsdatum
2009-05-29,
2009-05-30
Datum när rutan skannades. En ruta innehåller normalt punkter från mer än ett stråk.
Stråken kan ibland ha olika datum.
Klassificeringsprogramvara
TerraScan 009.006
Programvara och version som senast använts
för att klassificera laserpunkterna. Kan avse
automatisk eller manuell klassificering.
Klassificeringsnivå
1
Omfattning av klassificering. Se kvalitetsbeskrivningen, bilaga B.
Senaste klassificeringsdatum
2009-12-01
Datum när senaste klassificering är genomförd.
Klassnummer
1
Klasser som kan ingå. Se bilaga B.
LANTMÄTERIET
2015-10-01
11 (12)
Antal punkter
5143301
Antal punkter i varje klass.
Totalt antal punkter
12 023 602
Antal laserpunkter totalt.
Lägsta höjd
61,96
Lägsta höjdvärde i punktmolnet.
Medelhöjd
96,1
Medelhöjd i punktmolnet.
Högsta höjd
162,33
Högsta höjdvärde i punktmolnet.
Antal enda ekon
7 874 089
Antal laserpulser som bara gett ett enda eko.
Antal första ekon
1 978 341
Antal ekon som varit det första av flera.
Antal mellanekon
195 842
Kan vara både andra och tredje ekot, men
inte första eller sista.
Antal sista ekon
1 975 330
Antal ekon som varit det sista av flera.
Lägsta intensitet
0
Lägsta registrerade intensitet.
Medelintensitet
75
Genomsnittlig intensitet.
Högsta intensitet
255
Högsta registrerade intensitet.
På Lantmäteriets hemsida, Hojddata/Fakta-om-laserskanning/Planer-utfall-ochstodmaterial, under Läge och skanningsdatum för stråk, finns en fil för att koppla
ihop laserpunkt med datum för laserskanning.
LANTMÄTERIET
2015-10-01
12 (12)
4 Förändringsförteckning
Senaste förändring av produktbeskrivningen har en detaljerad beskrivning. Denna
tas bort när ny ändring tillkommer. I tabellen anges i vilken version av produktbeskrivning för Laserdata, ändringen införts. Datumet anger från vilken dag ändringen
gäller.
4.1 Senaste förändring
Version
2.1 - 2.2
Datum
2015-10-01
Orsak samt ändring mot tidigare version
Justering av texten om geografisk täckning då förändringar har skett i
Geolex, samt enstaka korrektur, även i bilaga B.
4.2 Tidigare förändringar
Version
Datum
Orsak samt ändring mot tidigare version
2.0
2015-02-25
Klassificeringsnivå 3 tillagd, kvalitetsbeskrivning som bilaga samt i samband med det större omstuvningar och revidering.
1.7
2014-01-01
”Ny” borttaget ur Ny nationell höjdmodell.
Lagt till ordet ”öppna” i avsnitt 2.4 Lägesnoggrannhet.
1.6
2013-06-01
Kompletterat med information om kantfiler, fjällskanning, intensitet, klassificeringsnivå och leverans i tabellformat.
Redaktionella ändringar samt framsida.
1.5
2012-12-04
Uppdaterat med länkar till nya lantmateriet.se.
1.4
2012-04-02
Kompletterat med information omklassificering för broar och dammar.
Språkliga och redaktionella korrigeringar.
1.3
2011-12-01
Ändringar rörande nytt leveransformat, XML, för metadata, kapitel 3, och
ändrad mappstruktur vid leverans.
Nytt fält i metadatafilen: kvalitetsanmärkning och några nya värden i kapitel 2.6.
1.2
2010-04-26
Uppgifter har lagts till under punkt 2.4 och 2.5.
Stavfel och andra mindre rättningar har gjorts.
1.1
2010-03-24
Metadata har kompletterats.
Bilaga A
Produktbeskrivning för Laserdata och Grid2+
Produktions- och skanningsområden
Produktionsområde
Sverige är indelat i ovanstående angivna produktionsområden (kartan till vänster).
Område A, B och C skannades med några få undantag under icke vegetationsperiod.
Övriga områden skannas oberoende av årstid.
Skanningsområde
Varje produktionsområde delas in i ett antal skanningsområden (blått i kartan till
höger) som i normalfallet täcker ett område om 25 x 50 km. I fjällen kan skanningsområdena ha en annan form, då skanningen har anpassats efter terrängen.
Varje skanningsområde får en unik beteckning, t.ex. 09P001, där de första siffrorna är
initialt planerat produktionsår, sedan följer produktionsområdets bokstavsbeteckning efterföljt av ett löpnummer.
Bilaga B
Produktbeskrivning för Laserdata och Grid2+
2015-10-01
Kvalitetsbeskrivning nationell höjdmodell
LANTMÄTERIET
Innehållsförteckning
1.
Inledning
3
2.
Generellt om kvaliteten
3
3.
Insamlingsmetod
4
4.
Beskrivning av uppnådd kvalitet
5
4.1.
Fullständighet (punkttäthet/täckning) - brist
6
4.2.
Tematisk noggrannhet – klassificeringsnoggrannhet
8
4.2.1.
Klassificeringsnivå 1
8
4.2.2.
Klassificeringsnivå 2
9
4.2.3.
Klassificeringsnivå 3
9
4.2.4.
Kända brister i klassificeringen
4.3.
6.
17
4.3.1.
Absolut noggrannhet
17
4.3.2.
Relativ noggrannhet
17
4.3.3.
Lägesnoggrannhet – för grid
17
4.4.
5.
Lägesnoggrannhet
10
Aktualitet
Övrigt som påverkar kvaliteten
19
20
5.1.
Skannat vid flera tillfällen/vid olika årstider
20
5.2.
Olika typer av skannersystem
22
Användbarhet
23
2(23)
LANTMÄTERIET
1. Inledning
I den här bilagan beskrivs den kvalitet som nationella höjdmodellen
generellt har, samt ges en översikt av brister som kan förekomma.
Kvaliteten redovisas med de kvalitetsparametrar som beskrivs i
standard SS-EN ISO 19157:2013 Geografisk information –
Datakvalitet.
Bilagan är gemensam för båda produkterna som utgör den nationella
höjdmodellen, dvs. Laserdata och Grid 2+ (terrängmodell). Terrängmodellens kvalitet är beroende av kvaliteten på laserdata och är
därför svår att beskriva separat.
2. Generellt om kvaliteten
Den uppmätta eller uppskattade kvaliteten gäller hela den hittills
skannade ytan. Det är därför viktigt att för aktuellt området ta del av
metadata så som skanningstidpunkt och punkttäthet. Det finns också
ett fritextfält i metadata i XML som kan innehålla kommentarer om
eventuella lokala brister som iakttagits vid bearbetningen av
skannerdata. God kännedom om lokala terrängförhållanden är också
värdefullt vid tolkning av data.
En bra terrängmodell kräver bra redovisning av markpunkter. För
att få bästa markrepresentation behöver skanning ske när så många
laserpulser som möjligt har en chans att nå ner till marken. Att
skanna på sommaren innebär att löv och undervegetation hindrar
laserpulser från att nå marken. På hösten kan höga växter hindra
laserpulserna från att nå marken, även om löven har fallit från
träden. Vintertid ger ett snötäcke tjockare än några centimeter en
felaktig bild av markytans form. På våren kan smältvatten i stora
mängder vara ett hinder för skanning, då vatten döljer markytan. Det
innebär att tidsfönstret för en bra skanning inte är så stort. Bästa
förutsättningar för skanning ges utanför vegetationssäsong, dvs. på
våren och sent på hösten. De delar av Sverige som i huvudsak täcks
av lövträd, dvs. södra Sverige (område A, B och C i bilaga A i
produktbeskrivningen) har med några få undantag skannats under
icke vegetationssäsong. Resten av landet har skannats oavsett årstid
och kan därmed ha områden med relativt tät vegetation.
Viktigt att förstå är att kvaliteten på terrängmodellen varierar med
terrängen och skanningstidpunkten.
3(23)
LANTMÄTERIET
3. Insamlingsmetod
Höjddata samlas in genom att marken skannas med laser från
flygplan.
Några fakta om skanningen (ungefärliga värden):
·
·
·
·
·
Punkttäthet: 0,5-1 punkt per kvadratmeter (ner till 0,25
punkter per kvadratmeter i kalfjällsområden).
Flyghöjd över marken: 1 700-2 300 meter (upp till 3 500 meter i
fjällområden).
Skanningsvinkel: ± 20º.
Stråkövertäckning: 20 %.
Träffyta på mark (footprint): 0,4–0,8 meter, beroende på
flyghöjd.
Laserskanningen görs med instrument som kan ge minst fyra ekon
från samma laserpuls. Alla punkter bevaras genom hela produktionskedjan, även felaktiga punkter på hög eller låg höjd.
Intensitet (styrkan hos den reflekterade laserpulsen) registreras för
varje eko. Det gäller dock ej Leica-instrumentet, där bara de tre första
ekonas intensitet registreras.
4(23)
LANTMÄTERIET
4. Beskrivning av uppnådd kvalitet
Kvalitetstema
Kvalitetsparameter
Uppnådd kvalitet
Datamängd
Beskrivs
i avsnitt
Fullständighet
Brist
98 % av den totalt skannade ytan uppnår:
Laserdata
4.1
(Punkttäthet/täckning)
Minst 0,5 punkter/m2. Vatten undantaget.
Minst 0,25 punkter/m2 i kalfjällsområden utan infrastuktur.
Endast sista ekot räknas.
Tematisk noggrannhet
Klassificeringsnoggrannhet
Riktigheten i klassificeringen är överlag bra, men viss felklassning förekommer.
Laserdata
4.2
Lägesnoggrannhet
Absolut noggrannhet
På öppna plana hårdgjorda ytor är medelfelet:
Laserdata
4.3.1
0,05 m i höjd
0,25 m i plan
På kalfjäll kan noggrannheten vara något lägre på grund av
högre flyghöjd.
Aktualitet
Relativ noggrannhet
0,1 m i medelfel för höjdskillnad mellan angränsande stråk.
Laserdata
4.3.2
Lägesnoggrannhet - för
grid
Vid interpoleringen till terrängmodell i gridformat sker en
försämring av noggrannhet i höjd, men terrängmodellens
lägesnoggrannhet ligger väldigt nära laserdatas absoluta
noggrannhet på 0,05 m i höjd. Försämringen är störst i kraftigt
kuperad terräng.
Terrängmodell
4.3.3
Se beskrivning.
Se beskrivning.
Laserdata
4.4
Terrängmodell
Användbarhet
Se beskrivning.
HMK standardnivå 1
Laserdata
6
Terrängmodell
5(23)
LANTMÄTERIET
4.1. Fullständighet (punkttäthet/täckning) - brist
Punkttätheten på mark varierar med terrängtyp, typ av vegetation,
under vilken årstid laserskanning skedde och en rad andra faktorer.
Denna variation innebär att det i vissa områden finns brister i punkttätheten, medan det på andra platser är hög punkttäthet. Av den
skannade ytan uppnår 98 % önskad punkttäthet eller tätare. Siffran
är ett medelvärde på punkttätheten på rutor om 10 x 10 m i varje
skanningsområde. Metadata redovisar hur punkttätheten varierar i
området, från mer än 1 punkt ner till 0 punkter per kvadratmeter (se
kap 3.3.1 respektive 4.4.1 i produktbeskrivning för Laserdata
respektive Grid2+). Det är viktigt att ta del av denna information,
samt skanningstidpunkt, för att kunna förstå kvaliteten och
användbarheten i områdets höjddata.
Tätast med punkter på mark är det på öppna ytor och där skanningsstråken överlappar varandra.
Brister i punkttätheten på marken finns framför allt i områden med
tät vegetation som har skannats under vegetationsperiod. Den täta
vegetationen har inte släppt igenom laserpulserna som istället har
stannat på vegetationen. I de fallen kan punkttätheten vara god i
laserpunktmolnet, men dålig i det klassificerade markpunktskiktet:
träffar på mark kan saknas helt inom en viss yta (se Figur 1, och
under avsnitt 4.2.4 Vegetation nedan).
I starkt kuperad terräng kan träffar på branta terrängformer (t.ex.
vertikala bergväggar) saknas helt. Detta beror på skanningsvinkeln
när flygplanet passerade. Om det dessutom är tät skog eller annan
tät vegetation i den branta terrängen försämras givetvis möjligheterna till markträffar (se Figur 2).
Ekon från laserpulser som träffat vatten, vissa hustak (Figur 21),
nylagd asfalt eller andra objekt som ger dålig reflektion kan saknas
helt. Då uppstår det hål i laserpunktmolnet.
Vatten absorberar de flesta pulserna, utom i rät vinkel mot ytan då
pulserna ofta reflekteras, som i en spegel. Det ger en oregelbunden
punkttäthet över vatten.
6(23)
LANTMÄTERIET
Figur 1: På grund av den täta vegetationen har inga laserpunkter nått ner till marken i
svackan. Klassificeringen blivit rätt, men det saknas markpunkter. När terrängmodellen
skapas kommer denna yta att interpolers. Vita punkter är oklassificerade, orange är mark.
Figur 2: I riktigt ojämn terräng
blir markytan generellt svår
att definiera, speciellt med låg
punkttäthet. Eftersom terrängmodellen skapas utifrån markklassificerade punkter kan det
i den här typen av bergssluttningar finnas stora lokala felaktigheter i terrängmodellen.
Foto: Andreas Rönnberg
7(23)
LANTMÄTERIET
4.2. Tematisk noggrannhet – klassificeringsnoggrannhet
Klassificering av Lantmäteriets laserdata sker i olika steg, med
succesiv kvalitetsförbättring, vilket resulterar i olika klassificeringsnivåer.
De klassificeringsnivåer som används är:
Klassificeringsnivå
Förklaring
1
Automatiserad klassificering av mark, vatten och övrigt.
2
Klassificering av broar och säkrad markklassificering av
dammar.
3
Förbättrad klassificering av vatten samt utjämning av
vattenytor i terrängmodellen.
Klassificeringsnivå redovisas i metadata per skanningsområde.
Resultatet från laserskanningen är en mängd laserpunkter med känt
läge i plan och höjd. Alla typer av objekt på och ovan markytan finns
representerade i detta punktmoln.
För att kunna framställa en terrängmodell som representerar markytan måste punkter på marken särskiljas från övriga laserpunkter. I
lasfilerna ingår dock samtliga laserpunkter, även eventuella felaktiga
punkter (hög- och lågpunkter). Inga punkter tas alltså bort.
De klasser som används är:
01
Oklassificerade punkter
02
Punkter på mark
09
Punkter på vatten
11
Punkter på broar
4.2.1.
KLASSIFICERINGSNIVÅ 1
Klassificering av laserpunkterna till mark, vatten eller övrigt görs
med automatiska metoder. Som stöd används bland annat vattenpolygoner och byggnadsytor från GGD (Lantmäteriets Grundläggande Geografiska Data). Tre olika parameteruppsättningar
används för olika terrängtyper utifrån en grov uppdelning av landet.
Den automatiska klassificeringen blir dock inte alltid rätt. En heltäckande kontroll av klassificeringen går inte att göra, men en begränsad manuell översyn, utifrån felsignaler och med stickprov, görs
efter den automatiska klassificeringen. Det är svårt att upptäcka om
punkter har hamnat i fel klass, så en liten mängd felaktigt klassificerade punkter kan kvarstå.
8(23)
LANTMÄTERIET
4.2.2.
KLASSIFICERINGSNIVÅ 2
Vid den automatiska klassificeringen hamnar broar i olika klasser.
Större broar kan bli klassificerade som mark, broar över vatten kan få
klassen vatten och i andra fall kan broar bli helt oklassificerade.
Därför görs en, till största delen automatisk, omklassificering av
broar.
För att identifiera broar används Trafikverkets databaser BaTMan
(Bridge and Tunnel Management), NVDB (Nationell vägdatabas) och
GGD, samt automatiska algoritmer. Klassificeringen av broar
innebär att laserpunkter på den del av vägbanan som har luft under
sig klassificeras som bro. Punkter på broar med mindre spännvidd
än 3 meter, liksom tunnlar, klassificeras inte till bro.
I samband med klassificeringen av broar görs även en säkrad markklassificering av dammar. Det innebär manuell klassificering så att
markklassificerade punkter finns längs hela fördämningens överkant. För att identifiera dammar används SMHI:s databas SVAR
(Svenskt vattenarkiv), samt GGD.
4.2.3.
KLASSIFICERINGSNIVÅ 3
För att förbättra klassificeringen av vad som verkligen var mark vid
tidpunkten för skanning och för att åstadkomma släta och väldefinierade vattenytor i terrängmodellen, behövs noggrant karterade
strandlinjer som överensstämmer med vattenståndet vid skanningstillfället. Karteringen baseras helt på laserdata, eftersom flygbilder
inte insamlas samtidigt med laserskanningen. Sjöar och vattendrag
nykarteras och höjdsätts utifrån punkter på vattenytan. Vattenyta för
hav, sjöar mindre än 0,25 km² och vattendrag smalare än 6 meter
ingår inte.
De nya polygonerna används sedan som brytlinjer vid framställningen av terrängmodellen. Det innebär att sjöar och vattendrag får
släta och väldefinierade ytor med en enda höjd i terrängmodellen.
Undantaget gäller för vatten som skannats vid olika vattenstånd, där
kan det bli nivåskillnader (se Figur 19). Strömmande vatten får luta,
men rinnriktningen är konstant och ytan slät. Forsar jämnas inte ut,
men får en väldefinierad strand.
Den ursprungliga klassificeringen till vatten i laserdata görs om för
hela området. Klassningen görs med hjälp av de nykarterade strandlinjerna och där dessa inte finns används GGD.
Arbetet påbörjades under hösten 2014 och ger klassificeringsnivå 3.
9(23)
LANTMÄTERIET
4.2.4.
KÄNDA BRISTER I KLASSIFICERINGEN
Vegetation
Det finns två orsaker till att vegetation ibland felaktigt kan klassificeras som mark. Dels gör begränsningar i tekniken för laserskanning att ekon på kort avstånd (cirka 2 meter) efter ett föregående
eko inte kan registreras, dels är vegetationen ibland så tät, t.ex.
spannmål eller buskage, att inga laserpulser kan tränga igenom till
marken (Figur 3-Figur 7). När något av dessa fenomen orsakar fullständigt bortfall av punkter på markytan, kan vegetationen ovanför
bli klassificerad som mark (Figur 7-Figur 8).
Figur 3: På detta grustag som är överväxt
med tätt lövsly har terrängmodellen hamnat
2–3 m över den verkliga markytan.
Foto: Andreas Rönnberg
Figur 4: Täta snår av lövsly växer snabbt
upp på många hyggen. På det här hygget
har höjdmodellen delvis hamnat cirka 2 m
över den verkliga markytan.
Foto: Andreas Rönnberg
Figur 5: Även på öppen mark med tätt
växande grödor kan höjdmodellen hamna
ovanför den verkliga markytan.
Foto: Andreas Rönnberg
Figur 6: I denna glesa men fuktiga tallskog
täcks marken helt av en relativt hög
blandad örtvegetation. Skanningen tränger
inte ner till marken när denna vegetation är
kraftig medan skanning under vår och sen
höst ger ett bättre resultat.
Foto: Andreas Rönnberg
10(23)
LANTMÄTERIET
10 m
Figur 7: Täta låga buskage, som enbuskarna på Ölands alvar, kan vara redovisade som
mark, dvs. små kullar i terrängmodellen. Eftersom inga laserpulser har trängt igenom
buskaget till marken, och buskarna har en ”mjuk övergång” till marken, uppfattas buskarna
som mark av algoritmen vid den automatiska klassificeringen.
11(23)
LANTMÄTERIET
Figur 8: Punkter på mark eller vatten saknas ofta i tät vass, och vassen har vid den
automatiska klassificeringen fått klassen mark. Vatten är här klassificerat utifrån
vattenpolygon i GGD (klassificeringsnivå 2). Röda strecket är profilen i bilden under.
Höjdskillnaden på ett par meter mellan mark och vass syns tydligt.
12(23)
LANTMÄTERIET
Brant terräng
I kraftigt kuperad terräng med skarpa höjdförändringar kan
markklassificeringen bli ofullständig. Orsaken är att algoritmen för
markklassificering kan uppfatta kraftigt avvikande punkter som
byggnader eller vegetation (Figur 9-Figur 13).
Grus- och bergtäkter och liknande föränderliga miljöer justeras inte i
klassificeringen.
Figur 9: Den här typen av terräng med branter och blockmark kan vara bristfälligt redovisad. Laserpunkter på branterna kan vara oklassificerade, och branterna blir därmed något
utslätade i terrängmodellen. Foto: Stig Lövborg.
10 m
Figur 10: Bergvägg vid väg som inte har blivit klassificerad till mark, sannolikt p.g.a. den
branta vinkeln. Profilen är dragen enligt det röda strecket i bilden nedan.
Figur 11:
Bergvägg intill
väg i markpunktskiktet visat med
terrängskuggning. Före och
efter rättning av
klassificeringen.
Profilen i bilden
ovan är dragen
vid det röda
strecket.
10 m
13(23)
LANTMÄTERIET
50m
Figur 12: Ytterligare ett exempel på en kulle som inte klassificerats till mark. Profilen är
dragen vid det röda strecket i bilden nedan.
100 m
Figur 13: Terrängskuggad bild över kullen i föregående bild. Till vänster är hur terrängmodellen ser ut efter den automatiska klassificeringen. Bilden till höger visar hur kullen
egentligen ser ut, dvs. efter omklassificering.
14(23)
LANTMÄTERIET
Byggnader
Byggnader kan bli felaktigt markklassificerade om de täcker en stor
yta eller om de har en mjuk anslutning till markytan (Figur 14-Figur
15). För att undvika detta används byggnadsytor från GGD för att
välja bort punkter på stora byggnader innan marklassificeringen.
Övriga byggnader hanteras i de flesta fall korrekt, och de fel som
ändå uppstår åtgärdas efter markklassificeringen med halvautomatiska metoder.
Inom tätbebyggda områden kan t.ex. upphöjda innergårdar och
parkeringshus med låga platta tak bli klassificerade som mark. Det
kan vara svårt, både med automatiska och manuella metoder, att
avgöra vad som faktiskt är mark, och var marknivån egentligen ska
gå (se Figur 16).
Figur 14: Göransson Arena, Sandviken,
var nybyggd och fanns inte med i GGD.
Byggnaden var inte heller tillräckligt
kantig för att uppfattas som byggnad
vid den automatiska klassificeringen
och blev därmed klassificerad som
mark. Den här typen av fel dyker oftast
upp i kvalitetskontrollerna och rättas av
operatörerna. Bilden visar laserdata
med terrängskuggning.
Figur 15: Annat exempel på en byggnad som inte ingått i byggnadspolygonerna från GGD,
och som vid den automatiska klassificeringen inte har uppfattats som byggnad.
15(23)
LANTMÄTERIET
Figur 16: Taket på ett parkeringshus har klassificerats som mark. Frågan är dock var marken
ska gå.
Lågpunkter
Lågpunkter uppstår vanligtvis på grund av reflexer i speglande ytor,
vilket gör att de är vanligast intill byggnader. Men lågpunkter kan
också uppstå på grund av tillfälliga fel i laserskannerns längdmätning. Lågpunkter elimineras normalt under klassificeringen, men
kan i sällsynta fall påverka slutprodukterna.
Vatten
Redovisningen av vatten har brister. En anledning till detta är att
vattnets reflekterande och absorberande egenskaper ger oregelbunden och bristfällig punkttäthet. Andelen felaktiga hög- och lågpunkter verkar också vara större över vatten.
På klassificeringsnivå 1 och 2 är vatten klassificerat med vattenpolygoner från GGD som mask. I några områden där GGD inte fanns
användes vägkartans vattenpolygoner. Dessa polygoner stämmer
inte alltid överens med den strandlinje som rådde vid skanningen.
Det kan bero på naturliga förändringar i strandzonen eller skillnader
i vattenstånd. Punkter i strandzonen kan därmed vara felaktigt
klassificerade som mark eller vatten.
Från klassificeringsnivå 3 är vattendrag bredare än 6 meter och sjöar
större än 0,25 km² omklassificerade och vattenytorna utjämnade i
terrängmodellen med hjälp av utifrån laserdata nykarterade strandlinjer. Övrigt vatten på nivå 3 har klassificerats om med aktuella data
från GGD.
Observera att om laserdata har samlats in vid olika tillfällen kan
vattenståndet ha förändrats och trappstegseffekter uppstå (se Figur
19).
16(23)
LANTMÄTERIET
4.3. Lägesnoggrannhet
4.3.1.
Absolut noggrannhet
För att säkerställa sömlösa höjddata (jämn övergång mellan skanningsområden) med hög noggrannhet över hela landet görs en inpassning av laserpunktmolnet. Som stöd för inpassningen används
punkter inmätta på marken (stödpunkter) i både plan och höjd i
SWEREF 99 TM. Därefter verifieras noggrannheten i både plan och
höjd mot inmätta kontrollpunkter.
Mätningar visar att medelfelet i höjd, på öppna plana hårdgjorda
ytor, ligger på 0,05 meter i genomsnitt. Medelfelet i plan ligger i
genomsnitt på 0,25 meter. Noggrannheten i plan hos enskilda laserpunkter är normalt många gånger sämre än i höjd. I någorlunda plan
terräng är detta inget problem, men i starkt sluttande terräng inverkar detta på noggrannheten i höjd, som försämras i takt med att
lutningen ökar.
I de kalfjällsområden som flygs på högre höjd kan noggrannheten bli
något sämre på grund av den högre flyghöjden.
4.3.2.
Relativ noggrannhet
Kontroller visar att skillnaden i höjd mellan angränsande skanningsstråk ligger på 0,1 meter i medelfel.
4.3.3.
Lägesnoggrannhet – för grid
Vid interpoleringen till en terrängmodell i gridformat sker en viss
försämring av noggrannheten i höjd. Störst är försämringen i kraftigt
kuperad terräng, medan terrängmodellens lägesnoggrannhet i de
flesta terrängtyper ligger väldigt nära laserdatas absoluta noggrannhet på 0,05 meter i höjd.
Försämringar vid interpolering står i relation till gridet storlek – ju
glesare grid, desto sämre kvalitet. Modellens kvalitet beror också på
och varierar med de här beskrivna kvaliteterna och brister i laserdata
- sämre noggrannhet i laserdata ger sämre noggrannhet i terrängmodellen.
Terrängmodellen skapas genom linjär interpolering i TIN
(Triangulated Irregular Network), med mark- och vattenklassificerade punkter som noder. I områden som saknar markpunkter kan
interpolering ske över stora avstånd. Se exempelvis Figur 1, där
markpunkter saknas och interpolering måste ske. Metadata i form av
punkttäthet på mark är därför mycket viktig eftersom den visar vilka
17(23)
LANTMÄTERIET
områden som har lägre detaljeringsgrad (se kap 4.4.1 i produktbeskrivning för Grid2+).
Om det saknas punkter på mer eller mindre vertikala bergväggar
däremot, behöver det inte vara ett problem, då interpoleringen där
kan bli korrekt, om omgivande punkter ligger bra (Figur 17).
På klassificeringsnivå 1-2 kan terrängmodellens vattenytor ha
störande artefakter från interpoleringen, trots att noggrannheten i
höjd hos laserpunkterna är hög. Eftersom terrängmodellen framställs
med linjär interpolering i TIN, kan mycket långa trianglar bildas
mellan stranden och närmaste punkt på vattenytan (Figur 18). Inom
sådana trianglar saknas laserpunkter helt, och om markpunkter
saknas närmast vattenytan kan den interpolerade höjden avvika
kraftigt från vattenytans verkliga höjd. På klassificeringsnivå 3
kvarstår inte detta problem för vattendrag bredare än 6 m och sjöar
större än 0,25 km², då dessa vattenytor har utjämnats. Däremot kan
en trappstegseffekt istället uppstå, om en vattenyta har skannats vid
olika tillfällen och vattenståndet har varit olika (se Figur 19).
Figur 17: I den här bergsbranten har klassificeringen inte blivit helt rätt. Men när terrängmodellen skapas och denna yta interpoleras blir terrängmodellen ändå ganska rätt.
18(23)
LANTMÄTERIET
Figur 18: Efter interpolering till terrängmodell i gridformat ansluter trianglar till punkter på
vattenytan eller i vissa fall till punkter vid den motsatta stranden. Det ger ojämna vattenytor. Gäller klassificeringsnivå 1-2, samt små sjöar, smalare vattendrag och hav på klassificeringsnivå 3.
4.4. Aktualitet
Laserdata
Laserpunktmolnet är en ögonblicksbild som Lantmäteriet inte
kommer att ajourhålla. Punktmolnet kan få nytt klassificeringsdatum
vid rättningar av felaktig klassificering, men punktmolnet i sig får
inte ny aktualitet. I den mån laserskanning används för ajourhållning
av terrängmodellen kan vissa områden komma att skannas på nytt
och laserdata får då ett nytt aktualitetsdatum i berörda områden.
Redovisningen sker i metadata.
Terrängmodell
Det finns ännu ingen beslutad ambitionsnivå för ajourhållning av
terrängmodellen.
19(23)
LANTMÄTERIET
5. Övrigt som påverkar kvaliteten
5.1. Skannat vid flera tillfällen/vid olika årstider
Ett helt skanningsområde kan sällan slutföras vid ett enda tillfälle.
Om tiden mellan de olika skanningstillfällena är lång kan t.ex.
vegetation och vattenstånd ha förändrats mellan olika skanningsstråk (se Figur 19). Det innebär att intilliggande områden kan se olika
ut i data, trots att det i verkligheten är ett homogent landskap. Det
kan också vara mänsklig påverkan som gör att landskapet skiljer sig
från en tid till en annan, exempelvis vid byggandet av en väg. Det är
alltså viktigt att kontrollera i metadata för det aktuella området om
det har skannats i flera sessioner och under långa tidsintervall. Om
så är fallet gäller det att vara särskilt uppmärksam på att skillnader i
bland annat vegetationen och vattenstånd mellan skanningstillfällena kan ha påverkat resultatet. Se Figur 20.
100 m
Figur 19: Horisontellt genom översta bilden går en gräns mellan två skanningsområden som
har skannats vid olika tillfällen. Vid den dragna profilen syns nivåskillnaden på vattnet i
laserpunktmolnet, vilket resulterar i trappsteg i terrängmodellen.
20(23)
LANTMÄTERIET
Figur 20: Överst detalj ur differensbild mellan två skanningar utförda vid olika årstider.
Nederst ett ortofoto över samma terrängavsnitt. Uppe till höger ett område med vass, i
mitten åkermark, och nere till vänster vass och annan låg vegetation. Differensbilden är
grön vid 0 m skillnad och blir helt röd respektive blå vid ±0.4 m skillnad i höjd.
21(23)
LANTMÄTERIET
5.2. Olika typer av skannersystem
Skanningen genomförs med flera typer av skannersystem. Terrängmodellen blir helt jämförbar men laserdata har vissa karakteristiska
skillnader som man bör vara uppmärksam på. Det som skiljer är bl.a.
att de olika skannrarna använder olika skalor för registrering av
intensitet, samt att de har olika skanningsmönster. I de flesta tillämpningar saknar dess skillnader praktisk betydelse.
De flesta system lagrar diskreta punkter i realtid, medan vissa lagrar
hela vågformen som sedan efterbearbetas till diskreta punkter. Den
sistnämnda tekniken ger högre punkttäthet i vegetation, samt bättre
redovisning av mark under låg vegetation, vilket minskar risken för
felaktigt markklassificerad vegetation.
För laserdata från Leica ALS70 är intensiteten korrigerad för effekten
av AGC (Automatic Gain Control).
De olika typer av skannersystem som har använts hittills är:
·
·
·
·
·
Leica ALS50-II
Leica ALS60
Leica ALS70-HA
Optech ALTM Gemini
Riegl LMS-Q780 (full vågform)
I metadata framgår vilket system som har använts för aktuell
skanning.
22(23)
LANTMÄTERIET
6. Användbarhet
Lantmäteriets laserdata och terrängmodell räknas som standardnivå 1 i HMK:s indelningssystem. Det betyder att trots att
Lantmäteriets höjddata håller hög kvalitet, innebär den förhållandevis låga punkttätheten att punktmolnet och terrängmodellen bäst
lämpar sig för nationell och regional översiktlig planering.
Figur 21: Ett exempel på hur glest laserträffarna (samtliga ekon) hamnar och den verkliga
träffytan. På bilden syns också ett exempel på att punkter saknas över ett tak.
Nationella höjdmodellen har sin grund i behovet av data för att
bättre kunna bedöma följderna av ett förändrat klimat och genomföra klimatanpassningsåtgärder. Men nationella höjdmodellen har
också använts i en rad andra tillämpningsområden såsom jord- och
skogsbruk, krisberedskap, projektering, exploatering, geologi,
arkeologi och orienteringskartor.
Läs gärna användares rapporter på Lantmäteriets hemsida, där
redovisas olika erfarenheter kring användandet av Lantmäteriets
höjddata.
23(23)