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:

ogr2ogr -f PostgreSQL "PG:dbname=uzpr host=gislab.fsv.cvut.cz user=xxx password=xxx" 602515_stav.db

Zadání¶

První část (SQL) - pro řešení úloh využijte dokumentace formátu VFK
Druhá část (Python+SQL)