06. Komplexní úlohy zpracování geoprostorových dat v Python

Komplexní úloha kombinující knihovny Fiona, Shapely, Rasterio, PDAL a GDAL.

Datové zdroje

• RUIAN

  • VDP

• DMR5G

Zadání

První část

1. Stáhněte data ve výměnném formátu RÚIAN (VFR) nejprve manuálně (VDP) a prozkoumejte jej pomocí ogrinfo -so -al (OSGeo4W Shell)
2. Vyřešte automatizované stahování VFR z URL adresy url = 'https://vdp.cuzk.cz/vymenny_format/soucasna/{datum}_OB_{obec}_UKSH.xml.zip'
3. Vyřešte problém s více geometriemi (Fiona nepodporuje zdroje s vice geometriemi) - pomocí GDAL převeďte vybranou geometrii (OriginalniHranice) do nové datové vrstvy (tip: gdal.VectorTranslate)
4. Data zpracujte pomocí Fiony: vytvořte novou datovou vrstvu, která bude obsahovat pouze parcely, pro které platí, že rozdíl výměry SGI a SPI je větší než vámi zadaná tolerance. Atributová tabulka bude obsahovat parcelní číšlo (kmenové číslo a pododdělení čísla), výměru SPI a rozdíl výměr.
5. Zjistěte kolik parcel nad tolerancí mělo plochu měřěnu až po roce 2020?

Druhá část

1. Stáhněte DMR5G Geoportál pro vámi zvolené území odpovídající první části zadání
2. Vytvořte DMR v rastrové reprezentaci pomocí PDAL
3. Pro každou parcelu (z první části) vypočítejte rozdíl vodorovné a skloněné plochy. Vypište celkový počet parcel a počet s lišící se plochou.