Zobrazení zdroje stránky 153ZODH / 1. cvičení

{{Cvičení|153ZODH|1|Úvod do GRASS GIS}}
== Osnova ==

Cvičení je zaměřeno na seznámení se systémem [[GRASS GIS]]. Více informací v rámci předmětu [[153YFSG|Free software GIS]].

Vysvětleny jsou základní pojmy jako je ''databanka'', ''location'' (lokace), ''mapset''. Spuštění systému GRASS, vizualizace jednotlivých mapových vrstev, demonstrace základních příkazů  s důrazem na modulární strukturu systému. Studenti se seznámí s [http://en.wikipedia.org/wiki/Command-line_interface rozhraním příkazové řádky] (CLI) a [http://en.wikipedia.org/wiki/Graphical_user_interface grafickým uživatelským rozhraním] (GUI) systému GRASS.

=== Seznam použitých příkazů ===

* {{GrassPrikaz|g.list}}
* {{GrassPrikaz|g.manual}}
* {{GrassPrikaz|d.mon}}
* {{GrassPrikaz|g.region}}
* {{GrassPrikaz|d.rast}}
* {{GrassPrikaz|d.vect}}
* {{GrassPrikaz|g.gui}}

== Spuštění systému GRASS, základní pojmy ==

''Poznámka: V dalším textu budeme předpokládat verzi [http://grass.osgeo.org/download/index.php GRASS 6.5]. Pokud používáte jinou verzi, očekávejte drobné rozdíly...''

GRASS spustíme příkazem

 grass65

GRASS nastartuje v ''grafickém'' či ''textovém'' módu (podle nastavení v souboru <tt>~/.grassrc6</tt>, toto chování lze ovlivnit parametry <code>-gui</code> a <code>-text</code>).

 grass65 -text

[[Soubor:ZOD-cv1-text-64.png|frame|center|Textová uvítací obrazovka GRASS 6.4]]

 grass65 -gui

[[Soubor:ZOD-cv1-wxgui.png|thumb|center|400px|Grafická uvítací obrazovka GRASS 6.4 (wxGUI)]]

=== GRASS Database ===

;'''DATABASE''' (databanka): je adresář umístěný na disku, např. <code>$HOME/grassdata</code>. V tomto adresáři jsou v podstatě uložena ''veškerá'' data, se kterými GRASS pracuje (tedy rastrové a vektorové mapové vrstvy, atributové tabulky, popisky apod.). Jednu z mála výjimek představují atributová (popisná) data skladovaná v některém z externích databázových systémů (např. [http://www.postgresql.org PostgreSQL], [http://www.mysql.org MySQL]) apod.

;'''LOCATION''' (lokace): je adresář umístěný v databance. Obsahuje data, která souvisejí s daným projektem. Je definována souřadnicovým systémem (referenční elipsoid, kartografické zobrazení, mapové jednotky) a velikostí zájmového území.

;'''MAPSET''' (mapset): je souborem map, které tvoří jakýsi logický celek v rámci lokace (daného projektu). Může např. odpovídat jednotlivým uživatelům (to pokud na projektu pracuje více lidí) nebo uceleným analýzám (studium vegetace, záplavová území, ...). Každá lokace musí obsahovat alespoň jeden mapset s unikátním názvem PERMANENT. Ten většinou obsahuje základní datové vrstvy, ostatní mapsety jsou pracovní (zpracování vstupních dat, jejich analýza...).

Příklad GRASS databanky:

 $HOME
  |  
  +--- grassdata
  |    |  
  |    +--- zod
  |    |    |
  |    |    +--- PERMANENT
  |    |    |
  |    |    +--- landa
  |    |
  |    +--- nc_spm_08
  |    |    |
  |    |    +--- PERMANENT

=== Uživatelská rozhraní ===

==== Příkazová řádka ====

 grass65 -text

[[Image:grass-cli.png|thumb|center|640px|GRASS příkazová řádka]]

Grafické uživatelské rozhraní lze spustit z příkazové řádky pomocí příkazu {{GrassPrikaz|g.gui}}, např.

<source lang="bash">
 g.gui gui=wxpython
</source>

==== GIS Manager (Tcl/Tk) ====

Tcl/Tk GUI je starší GUI systému GRASS, které bylo v současnosti plně nahrazeno novým wxGUI (viz níže).

 grass65 -tcltk

[[Image:grass-tcltk.png|thumb|640px|center|GRASS TCL/TK GUI]]

==== [http://grass.osgeo.org/wiki/WxGUI wxGUI] (wxPython) ====

 grass65 -wxpython

nebo (výchozím GUI pro GRASS je wxGUI)

 grass65 -gui

[[Image:grass-wxgui.png|thumb|640px|center|GRASS wxGUI]]

== Příkazy systému GRASS ==

GRASS je ''modulární'' systém, disponuje poměrně rozsáhlou množinou malých, ale výkonných programů. To odpovídá koncepci UNIXu jako takového, daný program má za úkol vyřešit dílčí problém, měl by být co nejmenší a poměrně jednoduchý.

Jednotlivé příkazy (tzv. moduly) systému GRASS mají až na výjimky podobnou syntax a jsou jako celek poměrně konzistentní. Jejich jména se skládají z předpony označující skupinu příkazů a krátkého názvu napovídající účel modulu. Tak například modul {{GrassPrikaz|d.mon}} patří do skupiny "Display" a je určen pro manipulaci s tzv. GRASS monitorem (otevření, uzavření).

<center>
{|class="border"
|+Skupiny GRASS modulů
|-
! prefix !! skupina !! krátký popis
|-
|<tt>db.</tt>
|database
|práce s&nbsp;externími databázovými systémy
|-
|<tt>d.</tt>
|display
|grafické výstupy a&nbsp;vizuální dotazy
|-
|<tt>g.</tt>
|general
|obecné příkazy pro manipulaci s daty
|-
|<tt>r3.</tt>
|grid 3D
|zpracování 3D rastrových dat (voxels)
|-
|<tt>i.</tt>
|imagery
|zpracování obrazových dat
|-
|<tt>ps.</tt>
|postscript
|příkazy pro tvorbu mapových výstupů ve formátu PostScript
|-
|<tt>r.</tt>
|raster
|zpracování rastrových dat
|-
|<tt>v.</tt>
|vector
|zpracování vektorových dat
|}
</center>

Pokud tedy napíšete do konzole "g." a stisknete klávesu TAB (automatické doplnění příkazu), vypíší se všechny příkazy ze skupiny
modulů "General".

První modul, se kterým se seznámíme, bude {{GrassPrikaz|g.list}}. Na něm si ukážeme standardní chování GRASS modulů -- pokud spustíme modul bez parametrů, objeví se GUI dialog modulu.

[[soubor:Grass-g_list-wxgui.png|thumb|400px|center|GUI modulu g.list]]

Velmi užitečným parametrem je <code>help</code>, který zobrazí krátký popis modulu, jeho syntaxi (názvy parametrů a přepínačů), tak např.:

<source lang="bash">
 g.list --help
</source>

<pre>
Description:
 Lists available GRASS data base files of the user-specified data type to standard output.

Keywords:
 general, map management

Usage:
 g.list [-f] type=datatype[,datatype,...] [mapset=string] [--verbose]
   [--quiet]

Flags:
  -f   Verbose listing (also list map titles)
 --v   Verbose module output
 --q   Quiet module output

Parameters:
    type   Data type
           options: rast,rast3d,vect,oldvect,asciivect,icon,labels,sites,
                    region,region3d,group,3dview
  mapset   Mapset to list (default: current search path)
</pre>

Rozsáhlejší nápovědu nabízí {{GrassPrikaz|g.manual}} s přepínačem <code>-m</code> zobrazí nápovědu ve formátu manuálové stránky.

Vraťme se k modulu {{GrassPrikaz|g.list}}, který vypíše seznam dostupných datových vrstev lokace. Seznam rastových a vektorových vrstev vypíšeme příkazem:

<source lang="bash">
 #seznam dostupných rastrových mapových vrstev
 #
 g.list rast
 
 ----------------------------------------------
 raster files available in mapset PERMANENT:
 rad1  rad2  spot1  spot2  spot3  tm1  tm2  tm3  tm4  tm5  tm6  tm7
 
 ----------------------------------------------

 #seznam dostupných vektorových mapových vrstev
 #
 g.list vect
 
 ----------------------------------------------
 vector files available in mapset PERMANENT:
 silnice  voda  zelez 
 
 ----------------------------------------------
</source>

V našem případě data obsahuje pouze mapset PERMANENT, pracovní mapset je zatím prázdný. Mapové vrstvy z mapsetu PERMANENT můžeme z pracovního mapsetu pouze číst, nelze je modifikovat. Mapset v tomto ohledu (do jisté míry) zohledňuje uživatelská přístupová práva (právo čtení, zápisu).

== Vizualizace mapových vrstev ==

=== Příkazová řádka ===

Data lze vizualizovat v grafickém okně tzv. "GRASS monitoru", těchto grafických oken může být celkově spuštěno 6 (označují se x0 až x6).
Pro manipulaci s GRASS monitory slouží modul {{GrassPrikaz|d.mon}}, pro otevření okna slouží parametr <tt>start</tt>, pro výběr <tt>select</tt> a pro zavření monitoru <tt>stop</tt>. ''Poznámka:'' Tyto příkazy nefungují pod MS Windows.

<source lang="bash">
 # otevřít první GRASS monitor
 #
 d.mon start=x0
 #
 # otevřít druhý GRASS monitor
 #
 d.mon start=x1
 #
 # vybrat první GRASS monitor
 #
 d.mon select=x0
 #
 # zavřít první GRASS monitor
 #
 d.mon stop=x1
</source>

Názvy parametrů či přepínačů lze libovolně zkracovat či dokonce vypouštět (to v případě, že se jedná o první parametr). Výše uvedený příklad lze tedy alternativně zapsat:

<source lang="bash">
 d.mon x0
 d.mon x1
 d.mon se=x0
 d.mon sto=x1
</source>

Důležitým pojem je tzv. '''region'''. Region je část zájmového území (tvar pravidelného čtyřúhelníka), které je definováno hraničními souřadnicemi SEVER JIH ZÁPAD VÝCHOD. Kromě rozsahu území nese informace o prostorovém rozlišení (ve směru sever-jih, východ-západ) pro rastrová data. ''Výchozím regionem'' je region, který pokrývá (měl by) zájmové území lokace. Pro manipulaci s aktuálním regionem slouží modul {{GrassPrikaz|g.region}}. Mezi podstatné volby modulu patří přepínač <tt>-d</tt>, který nastaví výchozí region a parametry <tt>rast</tt> a <tt>vect</tt>, které nastaví region podle dané mapové vrstvy.

Chceme-li tedy zobrazit družicový snímek '''tm1''', je vhodné nejprve nastavit aktuální region podle této rastrové vrstvy (tedy území, které
bude zobrazeno v GRASS monitoru) - rozsah a rozlišení.

<source lang="bash">
 # nastavit aktivní region podle mapové vrstvy
 #
 g.region rast=tm1
</source>

Pro zobrazení rastrových dat se používá modul {{GrassPrikaz|d.rast}}, pro vektorová data existuje podobně {{GrassPrikaz|d.vect}}. Oba moduly  disponují velkou množinou přepínačů a parametrů, které ovlivňují zobrazený výsledek na monitoru. 

Demonstrační příklad:

<source lang="bash">
 # zobrazit v GRASS monitoru vybranou rastrovou a vektorovou mapovou vrstvu
 #
 d.rast tm1
 d.vect zelez color=yellow
</source>

[[Soubor:ZOD-cv1-d.rast_vect.png|frame|center|Vizualizace rastrových a vektorových dat v GRASS monitoru]]

=== wxGUI ===

<!--[[Soubor:grass-d_rgb-wxgui.png|thumb|center|640px|GRASS wxGUI - přidání RGB vrstvy]]-->
[[Soubor:wxgui-tm453.png|thumb|center|640px|GRASS wxGUI - přidání RGB vrstvy]]

----
< [[153YZOD Zpracování obrazových dat|Stránky předmětu]] {{bullet}} [[153YZOD Zpracování obrazových dat - cvičení 2|Další cvičení]]

{{ZOD}}
{{GRASS}}