Paikkatietojärjestelmät/Esimerkit/Vinovalovarjosteen tuottaminen Maanmittauslaitoksen laserkeilausaineistosta
Maanmittauslaitoksen laserkeilausaineistojen käsittely
[muokkaa | muokkaa wikitekstiä]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
[muokkaa | muokkaa wikitekstiä]Maanmittauslaitos tarjoaa laserkeilausaineistot 3×3 km ruuduissa.[1] 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.[2]
Pistepilviaineiston tuominen paikkatieto-ohjelmistoon
[muokkaa | muokkaa wikitekstiä]- Ohjelmistoversio: 6.4.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 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:
Maanpintamallin luominen
[muokkaa | muokkaa wikitekstiä]- 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
[muokkaa | muokkaa wikitekstiä]- 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
[muokkaa | muokkaa wikitekstiä]- Laserkeilaus Maanmittauslaitos.
- Laserkeilausaineisto Maanmittauslaitos.
- Jakob Ventin: Ohjelmia Maanmittauslaitoksen tuottaman laserkeilausaineistojen prosessointiin ja visualisointiin
Viitteet
[muokkaa | muokkaa wikitekstiä]- ↑ Laserkeilausaineisto Maanmittauslaitos.
- ↑ LIDAR GRASS-Wiki.