![View this Page](/schemes/simple/viewie.gif)
![Edit this Page](/schemes/simple/edit.gif)
![Uploads to this Page](/schemes/simple/upload.gif)
![History of this Page](/schemes/simple/history.gif)
![Top of the Swiki](/schemes/simple/top.gif)
![Recent Changes](/schemes/simple/recent.gif)
![Search the Swiki](/schemes/simple/search.gif)
Druhé cvičení - vstup údajů
Načtení obrazu do grassu
- vytvoření nové lokace X-Y
- nazvěte ji brno-xy
- velikost ponechte implicitní, ta se nastaví později podle velikosti importovaného obrazu
- stáhněte obrázky: Mapa1 Mapa2
- r.in.gdal -e input=map1.tiff output=map1 (import obrazu, -e rozšíří hranice lokace)
- r.in.gdal -e input=map2.1.tiff output=map2
- g.region -p
- g.region rast=map1.red
- g.list rast (mrkněte, co r.in.gdal vyprodukoval)
- r.composite red=map1.red blue=map1.blue green=map1.green output=map1
- r.composite red=map2.red blue=map2.blue green=map2.green output=map2
- mapy map1 a map2 už jsou v barvách
- g.region rast=map1.red
- g.region -p
Příprava geokoordinace
- smyslem je projektovat obrázky do lokace, která má prostorový charakter (na rozdíl od X-Y)
- založte novou lokaci v Lat-Lon (jméno lokace brno-latlon, datum=wgs84)
- okraje dáme N=50, S=49, W=16:20, E=16:50
- implicitní rozlišení 1/10 vteřiny v obou směrech (00:00:0.1)
- vytvořte lokaci a vraťte se k brno-xy
- zadání geo-datumu
Postup geokoordinace
- založíme skupinu map pro projekci ven z lokace: i.group group=brno input=map1,map2
- zadáme skupině její cíl projekce: i.target group=brno location=brno-latlon mapset=PERMANENT
- zadávání "vlícovacích" (referenčních) bodů. Jejich počet závisí na zvoleném stupni transformačního polynomu (v našem případě je stupeň 1 a tomu odpovídá potřeba čtyřech bodů)
- i.points
- v příkazové řádce zvolte skupinu "brno"
- v monitoru vyberte jako podklad pro zadávání bodů mapu map2
- body jsou 4, jejich poloha je uprostřed písmena "O" v nápisu "GEODIS" (jsou zvýrazněné červeně)
![Uploaded Image: crtek.jpg](/vyuka-gis/uploads/26/crtek.jpg)
- souřadnice bodů jsou
- 1: 16:28:43.34E 49:14:0.55N
- 2: 16:28:41.35E 49:12:54N
- 3: 16:34:32.68E 49:12:46.24N
- 4: 16:32:20E 49:14:33.21N
- kontrola korektnosti zadaných bodů: svítí zeleně, funkce "ANALYZE" dole v okně (ukazuje i odhad chyby)
- kontrola souborů: ~/grass/brno-xy/$vas_mapset/group/brno
- rektifikace: i.rectify -a group=brno extension=obr order=1
- přesuneme se do cílové lokace, tam by měl být výsledek projekce
V lokaci brno-latlon
- máme tu hotové mapy map1obr, map2obr
- vykreslete je, upravte patřičně region
- prověřte geografickou správnost mapové vrstvy (d.where, d.measure)
- půjdeme na vektorizaci bodů, linií a ploch
Vstup geoTiffu (je to DEM pro danou lokalitu)
- stáhněte: demb.tiff
- r.in.gdal -e input=demb.tiff output=dem
- zobrazte v monitoru v kontextu geokoordinovaných map
- nviz
Vektorizace
- k vektorům budeme pořizovat tabulky s atributy
Bodová vrstva "obcebod"
- v.digit -n map=obcebod
- nastavení digitalizačního prkna (druhé tlačítko zprava):
- background - napsat příkaz pro vykreslení podkladu (d.rast map1)
- table - zadat atributy tabulky - nazev, obyv (int). Pak "Create table"
- panel, tlačítko "redraw"
- panel, tlačítko "+"
- db.select sql="select from obcebod"
cat|nazev|obyv
1|Kninicky|1000
2||0
3|Komin|20000
4|Zabovresky|20000
5|Medlanky|1000
- d.vect obcebod display=attr attrcol=nazev
Liniová vrstva "silnice"
- v.digit -n map=silnice bgcmd="d.rast map1obr"
- v nastavení vytvořte tabulku (bez dalších atributů, ty lze dodat později)
(dodělání sloupce)
$cat >l
alter table silnice add ajdy integer
$db.execute input=l
$
![Uploaded Image: topo1.jpg](/vyuka-gis/uploads/26/topo1.jpg)
- vložení nového bodu (funkce "split line")
Vrstva polygonů "parcely" (atribut číslo, majitel, ...)
- vytvořte tabulku
- budou se zadávat boundaries v módu "no category" (kategorie se zadává až pro centroid)
- pořiďte několik ohraničení, pak vytvořte centroidy (v módu s kategor