# Bodová pole do CSV

## Úvod
Mějme prázdný projekt `ceska_republika.qgz`, projekt osahuje pouze vrstvy z WMS ČUZK (Základní mapa a Letecké snímky) a podkladovou mapu OpenStreetMap. Pozor, OpenStreetMap dlaždice jsou původně v jiném souřadnicovém systému a do S-JTSK jsou transformovány. Proto jsou zobrazen "kostrbatě".

Další vrstva, kterou projekt zobrazuje, je vrstva digitálního modelu České republiky ze zdroje https://land.copernicus.eu/imagery-in-situ/eu-dem/eu-dem-v1.1

Připravený projekt stáhnete z adresy https://gismentors.cz/tmp/ceska_republika.zip . Projekt má cca 1GB - je to způsobeno digitálním modelem terénu.

## Tvorba vrstevnic

Sice lze vypočítat vrstevnice na celém území ČR, ale trvalo by to nějakou dobu (v řádu desítek minut). Proto si nejprve uděláme výřez zájmového území.

1. zazoomuje na zájmovou oblast
2. v menu `Processing` zaškrtněte `Toolbox`
3. na pravé straně QGIS okna se objeví tzv. processing tool box, v něm vyhledejte nástroj `clip raster by extent` - ořez rastru podle zadaných souřadnic.


4. v dialogovém okně nástroje zvolíme vstupní vrstvu, ze které chceme vyříznout zájmovou oblast (`eudem`), souřadnice zájmové oblasti (v menu zvolte `map extent`) a v tzv. pokročilých parametrech (`advanced parameters`) určitě zvolte vysokou kompresi (high compression)


Na první pohled se v mapě nic nezměnilo - ale v seznamu vrstev přibyla vrstva `Clipped` - je to právě ořízlé zájmové území.

5. V processing toolboxu (napravo) najdeme nástroj `Contour`

6. V dialogovém okně zvolíme: vstupní vrsvu náš ořez `clipped`, zbytek můžeme ponechat (interval vrstevnic 10m, atribut výšky nového vektorového souboru `ELEV` a další.).

7. Po výpočtu se mi vrstva přidala do mapy, ale nikde se mi neukázala. Bylo to proto, že vrstva neměla nastavený souřadnicový systém. Musel jsem jít přes `pravé tlačítko myši` do menu vrstvy a v záložce `zdroj (source)` nastavit souřadnicový systém na S-JTSK

8. V symbologii vrstvy jsem nastavil nejprve `Graduated` symbol, podle atributu `ELEV`, zvolil jsem barevný přechod a dole dialogu kliknul na tlačítko `Classify (klasifikovat)` - tím se mi vytvořila legenda vrstevnic v pěti skupinách.

Po uložení se mi objevily vrstevnice barevně naškálované podle výšky terénu.


## Náhodná bodová vrstva
Nejprve vytvoříme náhodný vzorek bodů. Počítám, že byste raději používala body z nějakého datového zdroje (např. Shapefile) nebo vlastní editací.

Jak editovat data ručně máme popsáno v https://training.gismentors.eu/qgis-zacatecnik/vektorova_data/editace.html

Jak otevřít datovou vrstvu CSV, máme popsáno v https://training.gismentors.eu/qgis-zacatecnik/vektorova_data/import_delim.html?highlight=csv

1. V `toolboxu` jsem vyhledal nástroj `Random points in extent` (náhodné body na území)

2. V dialogu jsem zvolil použití aktuálního mapového výřezu (`Use Map Canvas Extent`) a nastavil jsem počet bodů na větší než 1 (v mém případě 200)

3. Po potvrzení vznilka nová datová vrstva `Random points` - čistě náhodně rozmístěné body.

## Uložení výšky terénu do bodové vrstvy a export do CSV
1. V toolboxu najděte nástroj `Sample raster values` (vzorkovat rastrové hodnoty)

2. V nástroji pak zvolte vstupní bodovu vrstvu - v mém případě `Random points` a rastrovou vrstvu - `eudem`.


3. Po potvrzení vznikne nová vrstva `Sampled`, jejíž atributová tabulka obsahuje jak identifikátor prvku `id` tak hodnotu z rastrové mapy `eudem` ve sloupečku `SAMPLE_1`.

4. Pravým tlačítkem na vrsvě, najděte v menu `Export` a uložte prvky jako soubor CSV.

5. V dialaogu pro uložení do souboru zvolte formát CSV, název výstupního souboru a ve volbách vrstvy `Layer options` nastavte, že chcete uložit také geometrii `GEOMETRY` jako `XY` (tedy `AS_XY`).

6. Vznilký soubor lze pak načíst v tabulkovém kalkulátoru


## Tvorba pravidelného gridu

1. V nástroji `toolbox` najděte nástroj `Create grid`

2. V dialogu nástroje zvolte, že chcete vytvořit bodový grid (`Point`), nastavte rozsah zájmového území (podle vrstvy, podle mapového okna, vlatsní) a že horizontální a vertikální vzdálenost mezi body bude 100 (metrů).

3. Opět se vytvořila nová vrstva. Pokud jste byla hodně "odzoomovaná", vidíte v mapě hodně bodů, které se navzájem překrývají. Pokud chcete něco vidět, zvětšete si zájmovou oblast.


Doplnit výšku k bodům gridu už jsme popisovali výše.

# Isoplochy

Po nějakém zkoušení jsme došli k názoru, že nejlepší bude opravdu převod na rastrová data a zpětná interpolace izoploch. Celý proces bude nejlepší "zabalit" do hotového nástroje. Nejprve jednotlivé kroky:

1. Otevření vrstvy BaP.shp
2. převod čar na body pomocí nástroje `Extrahovat lomové body`
3. Vznikne nová vrstva `Uzly` - musíte jí nastavit správný souř. systém (EPSG:5514)
4. Pomocí nástroje ze skupiny `GDAL` se jménem `Rastrová analýza/Grid (Lineární)` nastavte rastrovou interpolaci
5. Vstup `Uzly`, Z hodnota `VALUE`, Profil `Vysoká komprese`
Výsledkem je rastrová mapa
6. Pomocí nástroje `GDAL / Extrakce rastru / Izoliniové polygony` nastavte
Vstupní vrstva: `Interpolovaný (Linerární)`, Interval mezi vrstevnicemi `0.0001` (nebo jiný), Atributy pro minimální a maximální výšku jsem změnil na `bapmin` a `bapmax`. Zaškrtnul jsem ještě 3D vektor - praktický význam to nemá žádný, ale třeba by se Vám to hodilo někdy do budoucna
Výsledkem jsem isoplochy
7. Isoplochy nastylujeme (`pravé tlačítko na vrstvu -> Vlastnosti -> Symbologie`) - `Kategorizovaný symbol`, hodnota `bapmin`, a smazat hodnotu 0 a prázdnou hodnotu - maximální v legendě

# Tvorba vlastního nástroje na isoplochy

Výše popsaný postup jsem v nástroji `Návrhář modelů` namodeloval - akorát tomu chybí nastavení souř. systému výsledné vrstvy na EPSG:5514
Nástroj je dostupný v `Toolboxu`
Po nastavení vstupního souboru a hodnoty atributu, který nás zajímá, provede všechny výpočty automaticky a přidá novou vrstvu `plochy` do mapy

Tímhle postupem už nebudete muset dělat všechny kroky ručně, stačí pouze spouštět vlastní nástroj. Připravám Vám projekt, kde je nástroj dostupný (ale můžete si ho uložit pro všechny budoucí projekty - a samozřejmě upravit).

# Export digitálnímo modelu do TXT souboru

Je to trochu hack, výstupní formát se jmenuje AAIGrid – Arc/Info ASCII Grid a QGIS ho moc neumí, jde ho ale dostat ven následujícím způsobem:

1. Najděte nástroj "Oříznout rastr podle rozsahu" 
2. Nastavte rozsah ořezu
3. Doplňkové parametry příkazové řádky zadejte `-of AAIGrid` 

Vytvoří se dočasný soubor, sice se přidá do QGIS, ale nejspíš se nic nezobrazí.

4. Ve vlastnostech nové vrstvy se podívejte na záložku `Informace` a tam je uvedena cesta k souboru, který bude končit jako `OUTPUT.tif` 
5. Soubor obsahuje data, která potřebujete, stačí ho zkopírovat a přejmenovat.

## Poznámky
* Všechny vytvořené vrstvy jsou dočasné a po zavření QGIS budou ztraceny - pokud je neuložíte (`pravé tlačítko -> uložit vrstvu`). V našich příkladech jsme uložili pouze jednu vrstvu jako CSV. I v tomto malém příkladu vzniklo hodně odčasných vrstev.
* Projekt `ceska_republika.qgz` je prázdný - obsahuje pouze pokladové mapy a vstupní digitální model reliéfu.
* Ostatní vrstvy jsou přidány jako služby serverů třetích stran (ČUZK, OpenStreetMap) - lze je použít pro orientaci, pro mapový výstup, ale nelze je využít pro další analýzu
* ČUZK i OpenStreetMap datasety lze získat v jejich "datové" podobě a lze nad nimi tvořit další analýzy. To ale není obsaženo v tomto příkladu