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 (i.group, i.target)
• Přiřazení bodů o známých souřadnicích (i.points)
• Vlastní proces transformace (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

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 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

Zkontrolovat nastavení regionu cílové location:

GRASS> g.region vect=cr res=50

Georeferencing

Princip georeferencingu

Spuštění modulu: Soubor ‣ Georektifikace

• Spustit georeferencer
• Zvolit zdrojový mapset /home/jachym/grassdata/xy/jachym

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

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

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

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

1. Natažení, posunutí
2. ,,Ohnutí`` kolem jedné osy
3. ,,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

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 (gis.m nebo i.points)