Lekce 3 - Data v GIS Obsah 1. Cíle lekce...........................................................................................................................................1 2. Geografická data...............................................................................................................................1 3. Rastrový a vektorový model prostorových dat..................................................................................2 4. Data v rastrovém modelu..................................................................................................................4 5. Data ve vektorovém modelu..............................................................................................................5 5.1 Špagetový (spaghetti)) model prostorových dat........................................................................5 5.2 Topologický model prostorových dat.........................................................................................5 6. Základní entity prostorových dat.......................................................................................................5 6.1 Geometrické entity.....................................................................................................................5 6.2 Topologické entity......................................................................................................................6 7. Vztahy mezi prostorovými daty.........................................................................................................8 8. Příklady datových modelů prostorových dat.....................................................................................9 8.1 Obecný model ploch..................................................................................................................9 8.2 Plochy parcel platné katastrální mapy.......................................................................................9 8.3 Plochy využití území (land-use)..............................................................................................10 8.4 Digitální mapa města..............................................................................................................10 1. Cíle lekce • vysvětlit pojem geograficky vztažená data, prostorová data • popsat vektorový a rastrový model prostorových dat • vymezit základní entity prostorových dat a vztahy mezi nimi 2. Geografická data Geografická proměnlivost reálného světa je nekonečně komplexní, k jejímu přesnému zachycení bychom potřebovali nekonečně velkou databázi. Data tedy musí být redukována na konečnou a spravovatelnou velikost. Pravidla, která jsou užita pro konverzi reálného světa do diskrétních objektů, vytvářejí datový model. Datové modelování v GIS se neliší od modelování běžných IS. Datový model každého informačního systému obsahuje objekty (entity) a relace mezi objekty (entitami)1. Mezi relacemi je významná relace hierarchie (vyjadřuje nadtyp - podtyp objektu), tato relace je acyklická2, nemusí to být strom (les). V datovém modelu GIS jsou vymezeny speciální - prostorové - objekty a speciální - prostorové - relace. Prostorové objekty a relace jsou pokud možno zobecňovány a standardizovány: architekt systému vybírá z předem stanovených tříd prostorových objektů a z předem stanovených tříd prostorových relací. Prostorové objekty se dále dělí dále na dva základní podtypy • geometrický objekt, který popisuje tvar objektu a lokalizaci objektu v prostoru - například souřadnicemi • topologický objekt, který vchází do prostorové relace s jinými (topologickými) objekty 1 Omezíme se na binární relace na množině všech objektů obsažených v datovém modelu. Relace mezi objekty (entitami) mohou mít různé názvy v závislosti na použité metodice. Například v návrzích objektově orientovaných systémů se často nazývají asociace. 2 Binární relace r c: MxM je acyklická, pokud platí: je-li n > 2, m^, n^,..., n^ e M, (m, ,1^),..., (nijj.i, mn) e MxM, pak mi * m^. 1 Příklady topologických relací: • objekt typu hrana začíná v objektu typu uzel a končí v objektu typu uzel • objekt typu prostorový objekt leží uvnitř objektu typu areál • objekt typu hrana tvoří hranici mezi dvěma objekty typu areál • objekt typu prostorový objekt leží {nalevo, napravo} od objektu typu hrana Příklad hierarchie jednoduchého datového modelu: třídy topologických entit třídy geometrických entit hierarchie (nadtyp-podtyp) topologický vztah vztah topologických a geometrických entit Příklad struktury jednoduchého datového modelu GIS 3. Rastrový a vektorový model prostorových dat Základními modely prostorové složky geografických dat jsou rastrové a vektorové modely. Rastrový model dělí modelovanou oblast (WOI) do pravoúhlé sítě buněk, každá buňka obsahuje jednotlivou hodnotu sledovaného údaje. Buňce odpovídá plocha sledované oblasti (území). Jedna množina buněk spolu s asociovanými hodnotami vytváří vrstvu. V databázi může být uloženo mnoho vrstev. Vektorový model využívá pro vymezení lokalizace diskrétní bodové nebo liniové elementy - úsečky. Diskrétní objekty vektorového datového modelu vznikají spojením úseček. Na rozdíl o rastrového datového modelu nemusí vektorový datový model pokrývat celé území. (Rastrový model říká co se vyskytuje kdekoliv - v každém místě oblasti, vektorový model říká, kde se cokoliv vyskytuje, zobrazuje umístění každého objektu.) V obou modelech je prostorová informace reprezentována základními (homogenními) jednotkami. V rastrových modelech jsou to pixely (buňky), ve vektorových modelech jsou homogenními jednotkami body, lomené čáry a polygony. Pokud rastrový a vektorový model popisují shodné území ve srovnatelné přesnosti, je v rastrovém modelu řádově více základních jednotek než ve vektorovém modelu. Rozměr základních jednotek v rastrovém modelu je konstantní, rozměr základních jednotek ve vektorových modelech je velice proměnlivý - lomená čára se může skládat ze dvou bodů nebo také z několika desítek tisíc bodů. 2 rastrový model vektorový model Příklad rastrové a vektorové reprezentace Vektorový model Rastrový model Výhody - jednoduchá datová struktura - snadné provádění některých operací- například mezi překrývajícími se vrstvami (overlay) - viz popis analytických funkcí - zachycuje data s proměnlivou prostorovou složkou (barevný digitální obraz) - nemusí být žádná primární znalost jevu pro získání dat, vzorkování lze provést automaticky Nevýhody - málo kompaktní (řídké matice) - obtížná reprezentace topologických vlastností - neestetický výstup reprezentovaných dat (pouze za cenu zvýšení hustoty dat) - nemožnost libovolného zvětšování obrazu Výhody - kompaktní datová struktura výhodné uchování topologických vlastností blíží se ručně kresleným mapám, možnost estetického výstupu dat Nevýhody - složitější datová struktura - obtížnější overlay operace - nevýhodná reprezentace dat s proměnlivou prostorovou složkou (digitální obraz) Volba datového modelu prostorových dat a datového modelu popisných dat je dána následujícími požadavky na vlastnosti implementace GIS: • stupněm integrace geografického informačního systému do "vnějšího" informačního systému (datový model GIS musí být přizpůsoben okolí, pokud je GIS komponentou širšího informačního systému) • základní architekturou informačního systému - centralizovaný x distribuovaný systém (komunikace v distribuovaném systému je složitější a tedy nevýhodná pro komplikované datové modely) • cenou pořízení dat a cenou údržby dat ve zvoleném modelu (pořízení i údržba dat ve složitějších modelech je dražší) • softwarovým vybavením, které je k dispozici • zkušenostmi a znalostmi klíčových pracovníků • historickými precedenty v dané oblasti Volba mezi rastrovým nebo vektorovým GIS není řešením otázky, který je lepší, ale otázky, za jakých podmínek je kterých z nich výhodnější a jakou máme možnost využít nejvhodnější řešení v konkrétních případech. Lze diskutovat o čtyřech oblastech: • přesnost souřadnic • rychlost analytických funkcí požadavky na velikost paměti 3 • charakteristiky jevu, který má GIS modelovat Přesnost souřadnic • u rastru závisí na přesnosti pořízení dat (družicového snímku, skenovacího zařízení), obtížně se reprezentují liniové objekty s malým průměrem • u rastru není často známo, jaký vztah má poloha hraničních bodů nebo středu buňky k reálnému kartografickému souřadnému systému • u vektorových dat může být přesnost prakticky libovolná, závisí pouze na metodě reprezentace souřadnic (počet bitů na jednu souřadnici) a na velikosti zobrazovaného území, velká přesnost však často nemá smysl, je nutné zohlednit vlastnosti vstupujících dat • přesnost vektorových dat se uplatní při zpracování dat z přesných geodetických měření, při vytváření map z těchto měření • některé jevy (hranice parcel, administrativních jednotek, osy vedení inženýrských sítí) jsou vhodné pro reprezentaci vektory s velkou přesností, u některých jevů (hranice druhů lesních porostů, rozšíření živočichů a pod.) jsou přesné hranic zavádějící Výpočetní rychlost • rastrová data mohou být velmi rychle zpracovány při odpovědích na většinu analytických dotazů včetně překryvů (vrstev), blízkosti, často stačí porovnání buněk v různých vrstvách • ve vektorových systémech je pro stejné odpovědi vyžadována větší výpočetní kapacita - musí být řešeny komplexní geometrické problémy například výpočet polohy a průsečíku dvou úseček • při výpočtech průniků polygonů vznikají singulární polygony Požadavky na paměť • v nejjednodušším režimu vyžadují rastrová data na jednu buňku jeden nebo dva byty, v těchto případech je podstatně omezen počet použitých řádků a sloupců • rastrová data je možné komprimovat několika metodami, stupeň komprese závisí na prostorové proměnlivosti dat • u vektorových dat jsou nároky na paměť malé, závisí na složitosti objektů a přesnosti souřadnic Charakteristika modelovaného jevu • rastrovou reprezentaci je možné použít tam, kde není apriori známa prostorová variabilita jevu • typická data v rastrovém formátu pocházejí ze satelitních nebo leteckých snímků a ze skenování • vektorová reprezentace poskytuje lepší možnosti, pokud se prostorová variabilita liší v různých oblastech území • typická data ve vektorovém formátu pocházejí z geodetických měření, zobrazují socio-ekonomické jevy, objekty utilit, vlastnické hranice, apod. 4. Data v rastrovém modelu Hodnoty uchované v buňce (pixelu) • hodnota uchovaná v buňce (pixelu) závisí jednak na kódované realitě, jednak na vlastnostech SW • možnosti: celá čísla, čísla v pohyblivé čárce, alfanumerické hodnoty • každá buňka - pixel - má pouze jednu hodnotu, což nepřesně interpretuje realitu (hranice jevu prochází buňkou), některé systémy umožňují uchovat více hodnot v jedné buňce a jejich podíly • rastrová data mohou být vizualizována jako množina mapových vrstev (mapová vrstva je množina dat popisující jednoduchou charakteristiku pro každou buňku uvnitř sledované oblasti) 4 • typické rastrové databáze obsahují více než sto vrstev, důležitou charakteristikou vrstvy je rozlišení (rozměr buňky), orientace (azimut, úhel mezi severem a směrem sloupců rastru) 5. Data ve vektorovém modelu 5.1 Špagetový (spaghetti)) model prostorových dat Pojem špagetový model vektorových dat\e odvozen od způsobu digitalizace papírové mapy a uchování digitalizovaných dat. Papírová mapa je převedena po čarách do posloupností bodů definovaných souřadnicemi X a Y - struktur typu lomená čára. Bodové prvky jsou převedeny do dvojic souřadnic -struktur typu bod. Soubor lomených čar a bodů nemá žádnou vnitřní strukturu - odtud termín špagetový model. V širším pojetí lze pojem špagetový model použít při charakteristice datového modelu, ve kterém jsou objekty prostorových dat uchovány bez vazeb na jiné objekty prostorových dat stejné nebo jiných tříd (například v sekvenčním souboru typu DXF nebo DGN). Ve špagetovém datovém modelu tedy nejsou topologické objekty. Špagetový model nepostihuje relace mezi prostorovými daty, vztahy musí být odvozeny výpočtem. Špagetový model je jednoduchá struktura vhodná pro počítačovou reprodukci map, ale i pro uchovávání prostorových dat, se kterými nebudou prováděny prostorové operace a které mají charakter pozadí - podkladu pro práci s jinými prostorovými daty. Výhodou špagetového modelu je levná údržba prostorových dat. 5.2 Topologický model prostorových dat Topologický model je nejrozšířenějším modelem, který uchovává prostorová data a vztahy mezi nimi. Tímto modelem lze reprezentovat areály i sítě. Základními objekty v tomto modelu jsou: • rovinný graf (arc - node model), • areálový rovinný graf • síťový graf s identifikovanými hranami • síťový graf s identifikovanými uzly Topologický model obsahuje topologické objekty, nemusí však obsahovat geometrické objekty (v tomto málo častém případě probíhá údržba prostorových dat bez jejich grafického vyjádření). Mezi hlavní výhody topologického modelu je možnost provádět analýzy prostorových dat bez užití souřadnic - bez práce s geometrickými objekty. Nevýhodou je nutnost údržby topologické struktury, která se často - zejména pro grafy s velkým počtem prvků - provádí dávkově. 6. Základní entity prostorových dat 6.1 Geometrické entity Bod typedef struct { long x; // souřadnice x long y; // souřadnice y } pointT; příklady: střed sídla, definiční bod budovy, šachta kanalizační sítě Lomená čára • posloupnost bodů • efektivní reprezentace lomených čar: přírůstky (offset) ve vektorech, Freemanův řetězový kód v rastrech • příklady: osa komunikační sítě, osa koleje, osa kanalizační sítě, vodovodní sítě, plynovodní 5 sítě, elektrické sítě, vodního toku Polygon, areál • polygon - uzavřená lomená čára • areál - posloupnost orientovaných polygonů • příklady: hranice fyzickogeografických oblastí (typy půdy, lesa, užití Země), hranice administrativních jednotek (okresy, obce, katastrální území) Sledované území může být pokryto plochami (areály) dvěma způsoby: 1. entity jsou izolované plochy, mohou se překrývat (příklad: plochy zasažené různými lesními požáry), prostor nemusí být plně pokryt entitami 2. každé místo patří právě jedné entitě (až na body hraničních čar), plochy plně pokrývají prostor První způsob pokrytí území může být konvertován do druhého. Z toho vyplývají dvě možnosti uchování ploch: • polygony - každý polygon je uložen jako posloupnost bodů, většina bodů je uložena dvakrát • diagraf rovinného grafu (arc-node model) -jsou uloženy hranice mezi areály, každý jen jednou s ukazateli na areály Pixel • primitiv - základní jednotka rastrových dat, důležitá je jeho velikost Řetěz pixelů • řetěz pixelů (např. ve Freemanově kódu) vymezuje liniový prvek v rastrovém modelu • příklady: osa liniového prvku (inženýrské sítě, komunikace), skelet liniového prvku Oblast pixelů • množina pixelů - vymezuje areálový prvek v rastrovém modelu • příklady: fyzickogeografické oblasti, oblasti land-use 6.2 Topologické entity Definice rovinného grafu Orientováný graf (U,H) je binární relace H na konečné množině U (H je tedy množina dvojic prvků z množiny U). Prvky množiny U nazýváme uzly (uzlové body), prvky množiny H hrany grafu. Grafy je běžné zaznamenávat pomocí kresby, pro tuto kresbu budeme používat termín diagraf. diagrafem grafu (U,H) budeme rozumět množinu bodů euklidovské roviny, které vzájemně jednoznačně korespondují množině U a množinu linií (oblouků) vzájemně jednoznačně korespondující množině H. Platí přitom, že žádná z linií jako svůj vnitřní bod neobsahuje bod korespondující s některým z uzlů a že linie odpovídající hraně tvořené dvěma uzly má za koncové body ty body, které s uzly korespondují. Stupeň uzlového bodu\e definován jako počet hran, se kterými uzlový bod inciduje (počet dvojic - prvků množiny H, ve kterých se uzlový bod vyskytuje). Graf\e rovinný, pokud existuje jeho diagraf takový, že linie diagrafu korespondující hranám grafu se protínají právě jen v uzlových bodech. Rovinný graf nazveme rovinným areálovým grafem, pokud jeho hrany tvoří hranice areálů v rovině. Každá hrana grafu je identifikována dvěma údaji: jmény (identifikacemi) areálů, jejichž hranici tvoří. Datové struktury areálového grafu Hrany typedef struct { long numVert; // poěet vrcholu 6 pointT *vert; graphNodT *ubl; graphNodT *ub2; idT idl; idT id2; } graphEdgeT; // vrcholy // poěáteěni uzlový bod // koncový uzlový bod // identifikace levého areálu // identifikace pravého areálu Uzly typedef struct { pointT int graphEdgeT } graphNodT; Hranice areálu typedef { long graphEdgeT } boundaryT; node; degree; edge[]; // souřadnice uzlového bodu // stupeň uzlového bodu // hrany incidující s uzlovým bodem numEdge; // počet hran tvořících hranici edge [ ] ; // seznam hran tvořících hranici Areál typedef struct { idT long boundaryT } areaT; id; numBoundary; boundary[ ] ; // identifikace areálu // počet hranic tvořících are // seznam hranic tvořících are 7 uzlové body O stupně 2 (8) Q stupně 3 (a,ß,s,()),ri,Q £ stupně 4 (y) hrany identifikace hrany identifikace a = (a, TJ) 11,0 h= (y,<|>) 11,15/2 b = (a,y) 11,12 i = (iO 21,15/2 c= (a,ß) 12,0 j = ( \ N \ \ \ K lomená čára, která tvoří hranici jedné nebo dvou ploch definiční bod plochy 8.2 Plochy parcel platné katastrální mapy cd Parcela / polohopis technické mapy plocha parcely +je tvořena +tvořI digitální katastrální mapa referenční katastrální mapa hranice parcely r lež f 1 ^1 v / \ le živ / 1 \l Pokud popisné číslo parcely neleží uvnitř parcely, musí být asociováno sliniovým popisným prvkem typu "vynášecí čára", jejíž jeden konec musí ležet uvnitř parcely. digitální katastrálni mapa referenční katastrálni mapa referenční číslo parcely digitální katastrální mapa referenční katastrální mapa popisné číslo parcely —- I 9 8.3 Plochy využití území (land-use) cd Příklad 2 - využití území (landuse) / plocha využití území definiční bod bloku území definiční bod bloku území::bodový popisný prvek definiční bod bloku území:: bodový popisný prvek:: definiční bod budovy 1..* 1.2 {ordered} hranice využiti územi hranice využiti územi::hranice parcely hranice využiti územi::liniový prvek polohopisu hranice využiti územi::liniový prvek polohopisu::obvod budovy 1 8.4 Digitální mapa města Z cdDMP digitálni katastrálni mapa referenční katastrálni mapa polohopis technické mapy výškopis inženýrská síť technické mapy autorizovaná hranice ochranného pásma hranice technického ochranného pásma autorizovaná inženýrská síť 10 cd vysvětlivky 1 dědičnost třída prv ní třída potomků dědičnost <- dědičnost dědičnost druhá třída potomků <- dědičnost třída potomků prv ní třídy agregace |o---- další třída třída potomků první I druhé třídy Z cdDMP-DKM referenční katastrální mapa referenční číslo parcely referenční katastrální mapa popisné číslo parcely polohopis technické mapy referenční katastrální mapa výškopis textový popisný prvek referenční katastrální mapa hraníce parcely DMP digitální katastrální mapa polohopis technické mapy referenční katastrální mapa výškopis liniový popisný prvek polohopis technické mapy referenční katastrální mapa obvod budovy spojení bodů polohopisu polohopis technické mapy referenční katastrální mapa výškopis liniový prvek polohopisu polohopis technické mapy | referenční katastrální mapa výškopis bodový popisný prvek polohopis technické mapy referenční katastrální mapa výškopis bod bodového pole d polohopis technické mapy referenční katastrální mapa výškopis bodový prvek polohopisu polohopis technické mapy referenční katastrální mapa výškopis bod polohopisu 11 cd spojení bodů polohopisu y digitální katastrální mapa polohopis technické mapy referenční katastrální mapa výškopis bodový prvek polohopisu 1 0..1 digitální katastrální mapa polohopis technické mapy 1 spojení bodů polohopisu digitální katastrální mapa polohopis technické mapy referenční katastrální mapa výškopis bod polohopisu referenční katastrální mapa výškopis liniový prvek polohopisu 1..* 1..* li..* 0..1 | cd DMP -polohopis technické mapy / digitální katastrální mapa referenční katastrální mapa výškopis textový popisný prvek plocha parcely digitální katastrální mapa referenční katastrální mapa výškopis bodový popisný prvek plocha budovy SL DMP polohopis technické mapy digitální katastrální mapa referenční katastrální mapa obvod budovy A A gitálnl katastrální mapa referenční katastrální mapa výškopis bodový prvek polohopisu digitálni katastrální mapa referenční katastrální mapa výškopis liniový prvek polohopisu ] spojení bodů polohopisu digitálni katastrální mapa referenční katastrální mapa výškopis liniový popisný prvek plocha využiti území plocha administrativně správní jednotky hranice městské části digitálni katastrální mapa referenční katastrální mapa výškopis bod bodového pole digitálni katastrální mapa referenční katastrální mapa výškopis bod polohopisu 12 cd DMP - inženýrské sítě technické mapy / autorizovaná inženýrská síť vodovod autorizovaná inženýrská síťm autorizovaná inženýrská síť silnoproud autorizovaná inženýrská síť slaboproud DMP inženýrská síť technické mapy neurčená síť autorizovaná inženýrská síťm teplovod autorizovaná inženýrská síť p rod u kto v od autorizovaná inženýrská síť potrubní pošta autorizovaná inženýrská síťí plynovod autorizovaná inženýrská síť kolektor cd DMP - hranice technického ochranného pásma / autorizovaná hranice ochranného pásma hranice pásma plynovodu DMP -£> hranice technického ochranného pásma autorizovaná hranice ochranného pásma hranice pásma dálnice <- autorizovaná hranice ochranného pásma hranice pásma celostátní dráhy autorizovaná hranice ochranného pásma \ hranice pásma radiokomunikaci autorizovaná hranice ochranného pásma hranice pásma letiště autorizovaná hranice ochranného pásma \ hranice pásma metra 13 cd DMP - kanalizace / neurčená síť neurčená síť neurčená síť plynovod plynovod plynovod potrubní pošta potrubní pošta potrubní pošta produktovod produktovod produktovod silnoproud silnoproud silnoproud slaboproud slaboproud slaboproud teplovod teplovod teplovod vodovod vodovod vodovod bodový popisný textový popisný liniový popisný prvek IS prvek IS prvek IS neurčená síť\ autorizovaná inženýrská síť neurčená síť plynovod inženýrská síť technické rrapy plynovod potrubní pošta -o kanalizace potrubní pošta produktovod silnoproud k1- plynovod potrubní pošta produktovod slaboproud slaboproud teplovod teplovod vodovod liniový prvek vodovod bodový prvek A polohopisu IS polohopisu IS neurčená si plynovod potrubní pošta produktovod slaboproud teplovod vodovod bod polohopisu IS spojení bodů polohopisu IS Připomínky a dotazy k obsahu lekce posílej, prosím, na adresu: Rudolf Richter, richter@fi.muni.cz 14