Druhé cvičení - vstup údajů

Cviko v ArcGISu

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

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

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

• kontrola topologie

• 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