14. Komplexní úlohy zpracování dat katastru nemovitostí
Komplexní úlohy zpracovávající data katastru ve formátu VFK v databázi PostGIS s rozšířením na automatizaci analýz pomocí Pythonu.
Datové zdroje
Použijeme anonymizované VFK soubory poskytované ČÚZK: https://services.cuzk.cz/vfk/anonym/
Potřebujeme dva VFK soubory obsahující data SGI a SPI pro vybrané katastrální území, např.:
Načtení dat v prostředí QGIS
Načtení VFK dat pomocí VFK pluginu
Do prostředí QGISu nainstalujeme VFK plugin (Plugins -> Manage and
Install Plugins
). V záložce All
najdeme VFK Plugin a
nainstalujeme jej.
Plugin otevřeme (Plugins -> VFK -> Otevřít prohlížeč VFK
) a zadáme
jako vstup dva výše uvedené VFK soubory (první soubor mapa a poté
spi_s_jpv).
Import vybraných tabulek do databáze PostGIS
Do databáze PostGIS naimportujeme tři tabulky: PAR
, DRUPOZ
a ZPVYPO
.
Vrstva PAR
byla načtena do mapové okna automaticky pomocí VFK
pluginu. Naimportovat tuto vrstvu můžeme do PostGISu standardním
postupem ve správci databází:
Note
Nezapomene rovnou vytvořit prostorový index pomocí volby Create spatial index
.
Important
Databázový systém PostgreSQL nerozlišuje ve výchozích chování malá a velká písmena. Proto je v tomto případě nutné zadávat názvy tabulek a atributů v uvozovkách, např. select count(*) from "PAR" where "PAR_TYPE" = 'PKN'
Užitečné může být cílovou tabulku pojmenovat malými písmeny a zaškrtnout volbu Convert field names to lowercase
. V tomto případě nemusíme v SQL příkazech používat uvozovky v názvech tabulek a sloupců.
Abychom mohli do PostGISu naimportovat tabulku DRUPOZ
, tak ji
nejprve musíme načíst do QGIS projektu. V okně Browser
najdeme na
disku soubor s příponou *.db
, který vytvořil VFK plugin v
kořenovém adresáři vstupních VFK souborů. Vybranou tabulku přidáme do
projektu a naimportujeme pomocí správce databází podobně jako vrstvu
PAR
. Stejným způsobem naimportujeme i tabulku ZPVYPO
.
Správce databází poté zobrazuje v cílovém schématu importované tabulky.
Tip
Import dat můžeme jednoduše provést z příkazové řádky (pod MS Windows použijte OSGeo4W Shell):
ogr2ogr -f PostgreSQL "PG:dbname=uzpr host=gislab.fsv.cvut.cz user=xxx password=xxx" 602515_stav.db PAR DRUPOZ ZPVYPO
Případně můžete do databáze PostGIS importovat celý obsah datasetu VFK:
Zadání
- První část (SQL) - pro řešení úloh využijte dokumentace formátu VFK
- Druhá část (Python+SQL)