Jak stáhnout výškopis a udělat z něj vrstevnice nebo stínovaný reliéf
České výškové modely z LiDARu jsou od roku 2023 zdarma, ale přijdou jako mračno bodů (LAZ), ne jako hotová mapa. Ukážeme, který model vybrat, jak z něj nástroji GDAL vyrobit vrstevnice, stínovaný reliéf i sklonitost, a jak na to bez příkazové řádky v QGIS.
Vrstevnice, stínovaný reliéf nebo sklonitost svahu nikde ke stažení hotové nenajdete, protože jde o odvozeniny, které si každý počítá podle svých potřeb (jiný interval vrstevnic, jiný úhel osvětlení). Vstupem je vždy výškopis, tedy digitální model terénu. Ten český je od 1. 7. 2023 otevřená data zdarma, ale přijde jako mračno bodů z laserového skenování, ne jako hotová mapa. Tenhle návod ukazuje cestu od staženého mračna k použitelné vrstvě.
1. Který model: DMR, nebo DMP?
- DMR 5G, model reliéfu (holá zem), nejpřesnější (0,18 m mimo les). Pro vrstevnice a sklonitost terénu je to první volba, ale je jen jako mračno bodů LAZ.
- DMR 4G je méně přesný (5m mřížka), zato ho ČÚZK nabízí i jako hotový raster TIFF, takže odpadá krok s interpolací mračna bodů.
- DMP 1G je model povrchu, tedy včetně budov a vegetace. Hodí se na zastínění nebo výšku zástavby, ne na tvar terénu. Pozor, jeho distribuce v této generaci končí 1. 1. 2027.
2. Stažení: mapové listy přes ATOM
Data se nestahují pro celou republiku najednou, ale po mapových listech SM5 (2,5 × 2 km). Na geoportálu ČÚZK najdete stahovací službu ATOM, kde vyberete dlaždice podle zájmového území. Pro menší lokalitu (pozemek, obec) stačí pár listů, pro kraj už desítky, počítejte s tím při výběru rozsahu.
3. Z mračna bodů (LAZ) na rastr
LAZ je komprimované mračno bodů, každý bod má souřadnici X, Y, Z. Pro vrstevnice nebo hillshade potřebujete pravidelný rastr. Převod zvládne PDAL, který LAZ čte nativně a umí ho rovnou zapsat jako GeoTIFF s interpolací mezi body:
# LAZ (mracno bodu) -> GeoTIFF rastr, rozliseni 1 m, IDW interpolace
pdal translate SM5_list.laz dmr.tif \
--writers.gdal.resolution=1 --writers.gdal.output_type=idw4. Vrstevnice, stínovaný reliéf a sklonitost (GDAL)
Z hotového rastru už je to jen jeden příkaz na výstup. Nástroje gdal_contour a gdaldem jsou součástí GDAL, které máte, pokud jste podle návodu o QGIS nainstalovali QGIS (GDAL je jeho součástí).
# Vrstevnice po 1 m do sloupce 'vyska'
gdal_contour -a vyska -i 1 dmr.tif vrstevnice.shp
# Stinovany relief (svetlo od severozapadu, 45 stupnu nad obzorem)
gdaldem hillshade dmr.tif stin.tif -z 1.5 -az 315 -alt 45
# Sklonitost svahu v procentech
gdaldem slope dmr.tif sklon.tif -pInterval vrstevnic (-i) volte podle měřítka: 1 m pro detailní pozemek, 5 nebo 10 m pro celou obec, jinak výstup zahltí mapu čarami. Parametr -z u hillshade je vertikální zveličení, hodnota kolem 1 až 2 dá realistický, ne přehnaně dramatický reliéf.
5. Bez příkazové řádky: totéž v QGIS
Stejné GDAL nástroje jsou dostupné i s klikáním, bez terminálu. V QGIS je najdete v Processing Toolboxu (menu Zpracování, Panel nástrojů) pod GDAL, ve skupině Rastrová analýza: Vrstevnice, Kopcovitý stín (hillshade) a Sklon. Mračno bodů LAZ navíc QGIS od verze 3.18 umí otevřít přímo jako vrstvu bodového mračna a prohlížet ho ve 3D, i bez převodu na rastr.
Mimo hranice ČR
DMR 5G pokrývá jen území Česka. Pro příhraniční nebo mezinárodní analýzy použijte globální model, třeba Copernicus DEM (30 m, novější a přesnější) nebo starší SRTM. Postup s gdal_contour a gdaldem je na nich stejný, jen bez kroku s LAZ, protože oba se distribuují rovnou jako rastr.
