**************
Georeferencing
**************
**Kroky nutné při georeferencingu**

* Založení cílové Location (cr)
* (Založení zdrojové  Location (xy))
* Import rastrových dat do GRASSu (command:`r.in.*`)
* Vytvoření skupiny mapových vrstev určených k transformaci (:command:`i.group, i.target`)
* Přiřazení bodů o známých souřadnicích (:command:`i.points`)
* Vlastní proces transformace (:command:`i.rectify`)

===============================
Založení zdrojové location (xy)
===============================
::

    GRASS> exit
    GRASS> grass63 -text

Založte location v souřadném systému XY::

    ...
    LOCATION: xy
    MAPSET:   jachym
    ...

    projection: 0 (x,y)
    zone: 0
    north:       1000
    south:       0
    east:        1000
    west:        0
    
    e-w res:     1
    n-s res:     1
    
    total rows: 1000
    total cols: 1000
    total cells: 1,000,000

    Do you accept this region? (y/n) [y] >

Import rastrového souboru::

    GRASS> r.in.gdal in=data/snih.png out=snehovamapa

Import se nezdařil, protože projekce Location neodpovídá projekci rastru,
přidejte parametr `-o`::

    GRASS> r.in.gdal in=data/snih.png out=snehovamapa -o

    GRASS> g.list rast

    ----------------------------------------------
    raster soubory dostupné v mapsetu jachym:
    snehovamapa.blue  snehovamapa.green snehovamapa.red   
    
    ----------------------------------------------

Složte tři barevené kanály do jednoho souboru::

	GRASS> g.region -p rast=snehovamapa.red
        
        ...
        north:      4962
        south:      0
        west:       0
        east:       7016
        ...

	GRASS> r.composite --help
	GRASS> r.composite r=snehovamapa.red g=snehovamapa.green b=snehovamapa.blue out=snehovamapa
	GRASS> d.rast snehovamapa


.. image:: ../../images/snih.png

Nyní jsou data naimportována do *zdrojové* LOCATION, vyskočte z GRASSu a
přepněte se do *cílové* LOCATION::

    GRASS> exit && grass63 grassdata/cr/jachym

Přepnout se můžete taky pomocí modulu :command:`g.mapset`::

    GRASS> g.mapset jachym location=cr

    Erasing monitors ...
    Cleaning up temporary files ...
    UPOZORNĚNÍ: Your shell continues to use the history for the old mapset.
    You can switch the history by commands:
    history -w; history -r /home/jachym/grassdata/cr/jachym/.bashhistory; HISTFILE=/home/jachym/grassdata/cr/jachym/.bashhistory

    GRASS> g.gisenv

Naimportojte vektorovou mapu hranic České republiky, poslouží nám jako
podklad::

        GRASS> v.in.ogr --help
	GRASS> v.in.ogr dsn=data/cr.shp out=cr
	GRASS> v.in.ogr dsn=data/cr.shp out=cr -o
	GRASS> g.region vect=cr res=50 -p
	GRASS> d.vect cr # zobrazit pomocí GIS Manageru

.. image:: ../../images/cr-vektor.png


Zkontrolovat nastavení regionu cílové location::

    GRASS> g.region vect=cr res=50



==============
Georeferencing
==============

Princip georeferencingu

.. image:: ../../images/georeferencing.png


Spuštění modulu: :menuselection:`Soubor --> Georektifikace`

* Spustit georeferencer
* Zvolit zdrojový mapset `/home/jachym/grassdata/xy/jachym`

.. image:: ../../images/georeferencing1.png

Založit novou skupinu rastrových map :file:`snehovemapy` a přidat do ní
mapy :file:`snehovamapa,snehovamapa.red, snehovamapa.blue, snehovamapa.green`

.. image:: ../../images/georeferencing2.png

Doplnit název mapové skupiny :file:`snehovemapy` a název mapy, která
bude zobrazena `snehovamapa`

.. image:: ../../images/georeferencing3.png

Nastavení Ground Control Points (Interaktivně nebo ručně)

.. image:: ../../images/georeferencing4.png

**Polynom *n*-tého řádu:**

#. Natažení, posunutí 
#. ,,Ohnutí`` kolem jedné osy
#. ,,Ohnutí`` kolem více os

Počet nutných bodů pro transformaci = ((n + 1) x (n + 2))/2

**Střední kvadratická chyba** by neměla přesáhnout polovinu rozlišení hodnoty cílové *Location*

Výsledek

.. image:: ../../images/georeferencing5.png

========================================
Co dělat, když GRASS Georeferencer selže
========================================

* Přepnout se do zdrojové location (xy)
* Zkontrolovat, je-li vytvořena skupina rastrových map ::

    GRASS> g.list group

* Nastavit cílovou location pro skupinu map ::
  
   GRASS> i.target group=snehovemapy location=cr mapset=jachym

* Vytvořit soubor s GCP (:command:`gis.m` nebo :command:`i.points`)