Paikkatietojärjestelmät/Esimerkit/Vinovalovarjosteen tuottaminen Maanmittauslaitoksen laserkeilausaineistosta
Maanmittauslaitoksen laserkeilausaineistojen käsittely
muokkaaTä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
muokkaaMaanmittauslaitos 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- 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- 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- 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- Laserkeilaus Maanmittauslaitos.
- Laserkeilausaineisto Maanmittauslaitos.
- Jakob Ventin: Ohjelmia Maanmittauslaitoksen tuottaman laserkeilausaineistojen prosessointiin ja visualisointiin
Viitteet
muokkaa- ↑ Laserkeilausaineisto Maanmittauslaitos.
- ↑ LIDAR GRASS-Wiki.