Paikkatietojärjestelmät/Esimerkit/Vinovalovarjosteen tuottaminen Maanmittauslaitoksen laserkeilausaineistosta

Maanmittauslaitoksen laserkeilausaineistojen käsittely

Tässä esimerkissä kuvataan eri tapoja hyödyntää Maanmittauslaitoksen laserkeilausaineistoja. Ensisijaisesti kuvataa kuinka voi tuottaa maanpintamallin ja vinovalovarjostekuvia.

Maanmittauslaitoksen laserkeilausaineistoja voi ladata Maanmittauslaitoksen latatuspalvelusta. Aineisto on saatavilla .laz-tiedostomuodossa.

.laz-tiedostojen käsittely

Maanmittauslaitos tarjoaa laserkeilausaineistot 3×3 km ruuduissa. Tiedosto on pakattu LAS-tiedosto, jonka pääte on .laz. Jotkut ohjelmat pystyvät suoraan lukemaan .laz-tiedostoja, toiset vain purettuja .las-tiedostoja ja joitain ohjelmia varten .las-tiedosto pitää muuttaa tekstitiedostoksi. .laz-tiedostojen käsittelyyn tarvitsee joko LAStools-ohjelmakokoelman tai LASzip ja libLAS -kirjastot.

GRASS GIS-ohjelmiston kehitysversion (7) voi kääntää LASzip ja libLAS -kirjastojen kanssa, jolloin se tukee suoraan .laz-tiedostoja. Nykyinen version (6.4.2) edellyttää libLASin tai LAStoolsin käyttöä LAS-tietojen muuttamiseksi tekstimuotoon.

Pistepilviaineiston tuominen paikkatieto-ohjelmistoon

GRASS GIS 6.4    

Ohjelmistoversio: 6.4.2

---

GRASS GIS 7.2    

Ohjelmistoversio: 7.2 (svn: )

---

(Rajaa tarkasteltava alue g.region-komennolla ensin.) v.edit map=vtemp tool=create --overwrite ; v.to.3d input=vtemp output=voutput type=point height=0 --overwrite ; v.db.addtable map=voutput --verbose columns="x double precision, y double precision, z double precision, class integer, time double precision, angle integer, intensity integer, return integer, returns integer, source integer" --overwrite ; for i in *.las ; do las2txt --parse xyzctairnp --stdout $i | v.in.ascii -znbri x=1 y=2 z=3 input=- output=vtemp separator=comma columns="x double precision, y double precision, z double precision, class integer, time double precision, angle integer, intensity integer, return integer, returns integer, source integer" --overwrite ; v.patch -naebz input=vtemp output=voutput --overwrite ; done ; g.remove type=vector name=vtemp

ArcMap 10    

ArcMap ohjelmistoversio: 10.0

ArcGIS Laajennokset: 3D-Analyst

---

1. Lisää LAStools-työkalulaatikko ArcToolboxiin.

Tarkempi ohje: Using the ArcToolbox window.

2. Pura .laz-tiedosto .las-tiedostoksi käyttämällä LAStools:iin kuuluvaa laszip-työkalua. Löydät sen ArcToolboxista.

3. Tuo LAS-aineisto pistepilveksi ArcMapiin LAS to Multipoint-työkalulla. Jos haluat tuoda vain maanpintapisteet valitse Input Class Code 2.

Tarkempia ohjeita:

- How to use the LAS To Multipoint tool

- LAS To Multipoint (3D Analyst)

Maanpintamallin luominen

GRASS    

ArcMap 10    

ArcMap ohjelmistoversio: 10.0

ArcGIS Laajennokset: 3D-Analyst

---

a. Interpoloi rasteripohjainen korkeusmalli tuoduista maanpintapisteistä 3D-Analystin Raster Interpolation työkaluilla.

Tarkemmat ohjeet: An overview of the Raster Interpolation toolset

Esimerkiksi Inverse Distance Weighted menetelmällä (IDW):

1) Valitse Z value fieldiin arvoksi Shape.Z

2) Valitse Output cellsize halutun tarkkuuden mukaan. Ota alkuperäisaineiston pistetiheys huomioon!

Tarkempi ohje: IDW (3D Analyst)

b. Interpoloi vektoripohjainen korkeusmalli Create TIN työkalulla

Tarkempi ohje: Create TIN (3D Analyst)

Vinovalovarjostekuvan luominen

ArcMap 10    

ArcMap ohjelmistoversio: 10.0

ArcGIS Laajennokset: 3D-Analyst

---

1. Luo vinovalovarjoste rasteripohjaisesta pintamallista Hillshade työkalulla. (3D Analyst → Raster Surface)

Tarkemmat ohjeet: Hillshade (3D Analyst)

Muualla verkossa

• Laserkeilaus Maanmittauslaitos.
• Laserkeilausaineisto Maanmittauslaitos.
• Jakob Ventin: Ohjelmia Maanmittauslaitoksen tuottaman laserkeilausaineistojen prosessointiin ja visualisointiin