Difference between revisions of "Prostorové databáze - cvičení"
m (→Cvičení 4 - Vytvoření tabulky s prostorovými daty) |
m (→Cvičení 4 - Vytvoření tabulky s prostorovými daty) |
||
Line 6: | Line 6: | ||
==Cvičení 4 - Vytvoření tabulky s prostorovými daty== | ==Cvičení 4 - Vytvoření tabulky s prostorovými daty== | ||
− | Příklad vytvoření tabulky s prostorovými daty | + | Princip práce s prostorovými daty v prostředí SŘBD Oracle. |
+ | |||
+ | Příklad vytvoření tabulky s prostorovými daty, včetně tvorby prostorového indexu: příklad [[Parcely]] | ||
==Příkládek== | ==Příkládek== |
Revision as of 09:53, 31 March 2014
Contents
- 1 Obsah a cíle cvičení z předmětu Objektově orientované technologie
- 2 Cvičení 1 - Instalace prostorové databáze a vizualizačního nástroje
- 3 Cvičení 2 - Prostředí klienta
- 4 Cvičení 3 - Opakování SQL
- 5 Cvičení 4 - Vytvoření tabulky s prostorovými daty
- 6 Příkládek
- 7 Cvičení 5 - Prostorové dotazy
- 8 Cvičení 6 - Import dat s prostorovým kontextem
- 9 Cvičení 7 - Vizualizace dat s prostorovým kontextem
- 10 Cvičení 8 -
- 11 Cvičení 9 -
- 12 Cvičení 10 - Zpracování vlastního projektu
Obsah a cíle cvičení z předmětu Objektově orientované technologie
Cvičení 1 - Instalace prostorové databáze a vizualizačního nástroje
Cvičení 2 - Prostředí klienta
Cvičení 3 - Opakování SQL
Cvičení 4 - Vytvoření tabulky s prostorovými daty
Princip práce s prostorovými daty v prostředí SŘBD Oracle.
Příklad vytvoření tabulky s prostorovými daty, včetně tvorby prostorového indexu: příklad Parcely
Příkládek
Vytvořte jednoduchou strukturu dat, která založí evidenci Vlastníků a jejich Parcel. Promyslete tak, aby bylo možné vlastnit více parcel (současně i postupně v čase).
Vytvoření tabulek a vložení dat
create table Vlastnik( id_vla NUMBER PRIMARY KEY, jmenovla varchar2(25) not null, prijmenivla varchar2(25) not null, adresavla varchar2(30) ); insert into vlastnik values (1, 'Jan', 'Trčka', 'Frýdek-Místek');
create table parcela (idParcela number primary key, ku varchar2(8), tvarParcely sdo_geometry );
INSERT INTO parcela VALUES( 0, 'ova01', SDO_GEOMETRY(2003, NULL, NULL, SDO_ELEM_INFO_ARRAY(1,1003,3), --obdelnik SDO_ORDINATE_ARRAY(1,1,5,10) ) );
INSERT INTO parcela VALUES( 1, 'bu02', SDO_GEOMETRY(2003, NULL, NULL, SDO_ELEM_INFO_ARRAY(1,1003,3), --obdelnik SDO_ORDINATE_ARRAY(2,6, 9,9) ) );
INSERT INTO parcela VALUES( 20, 'ova005', SDO_GEOMETRY( 2003, NULL, NULL, SDO_ELEM_INFO_ARRAY(1,1003,1), -- polygon SDO_ORDINATE_ARRAY(5,10 , 4,10 , 2,10 , 9,15 , 9,16 , 20,10, 5,10 ) ) );
create table vlastnictvi( idVlastnictvi number primary key, idVlastnik number(6) not null, idParcela number(6) not null, odVlastni date not null, -- datum nabyti parcely doVlastni date, -- datum pozbyti parcely constraints vla_vla foreign key (idVlastnik) references vlastnik, constraints vla_pa foreign key (idParcela) references parcela );
Vložení záznamů do tabulky Vlastnictví - ukázka změn vlastnictví v čase. Tatáž parcela je postupně vlastněna jedním druhým a opět prvním vlastníkem. Pro vložení hodnot datumu je využita konverzní funkce to_date, která umožňuje zadávat datum ve zvoleném formátu (zde DD-MM-YY).
insert into vlastnictvi values(1, 2, 0, to_date('01-01-07', 'DD-MM-YY'), to_date('01-01-08', 'DD-MM-YY')); insert into vlastnictvi values((2, 4, 0, to_date('01-01-09', 'DD-MM-YY'),to_date('01-01-11', 'DD-MM-YY')); insert into vlastnictvi values(3, 2, 0, to_date('01-01-12', 'DD-MM-YY'),null);
Vyhledání aktuálního vlastníka parcely s číslem id=0.
select y.jmenovla, y.prijmenivla from vlastnictvi x, vlastnik y where x.id_vla = y.id_vla and x.doVlastni is null and x.id_par = 0;
Vyhledání aktuálního vlastníka parcely s číselným označením 253.
select y.jmenovla, y.prijmenivla from vlastnictvi x, vlastnik y, parcela z where x.id_par = z.id_par and x.id_vla = y.id_vla and x.doVlastni is null and z.cisloparcely = 253;
Tvorba prostorového indexu
Nejprve je nutno "obeznámit" SŘBD s tím, že existuje tabulka s geometrickou složkou. Tyto informace se ukládají v systémové tabulce user_sdo_geom_metadata, kam je nutno zapsat název tabulky, jméno sloupce obsahujícího geosložku, rámec prostoru, nad kterým bude vytvořen prostorový index a souřadnicový systém, ke kterému jsou data vztažena.
insert into user_sdo_geom_metadata values ('parcela', 'tvarparcela', sdo_dim_array(sdo_dim_element('X',0,30,0.005), sdo_dim_element('Y',0,30,0.005)), null );
Ověření správně uložených informací:
select * from user_sdo_geom_metadata;
Požadavek na prostorové zaindexování geometrie parcel:
create index sp_par_idx on parcela(tvarparcela) indextype is mdsys.spatial_index;
Teprve nyní je možné klást prostorové dotazy a zobrazit geometrickou složku v GeoRaptoru.
Dotazy nad daty o vlastnících a parcelách
Jak bude vypadat průnik dvou parcel? Například, budu chtít postavit na dvou pozemcích (zakreslím si nový objekt přes hranice) a o kolik např. zemědělské půdy přijdu?
select p1.cisloparcely, p2.cisloparcely, SDO_GEOM.SDO_INTERSECTION(p1.tvarparcela, p2.tvarparcela,0.005) from parcela p1, parcela p2 where p1.cisloparcely = 1 AND p2.cisloparcely = xx; --xx je nový geoobjekt - tvar chystané stavby
možná i sum
vzdálenost mezi parcelami stejného typu - použít alter table a doplnit do Parcela nový sloupec
ALTER TABLE table_name ADD column_name column-definition; ALTER TABLE Parcela ADD typParcela varchar2(15);
další
nakrelsit si v Georaptoru nějakou parcelu - napr trojuhelnik, kruh zkopirovat si hodnoty geometrie a vlozit do tabulky parcela, zjisti si definici vloženého prvku - zapsat do tabulky parcel = pomocí insert = ukázat propojení definice objektu a metod, které s nimi pracují spočítat si její obsah - create shape = dole je výsledek update obsahu prvně vložených parcel - přes výpočet GeoRaptor update parcela set vymera = 34 where cisloparcely=253; dotázky - kolik je vech parcle, jak velké území zabírají, kolik pozemků ma jeden vlasník kolik pozemku vlastni jeden vlasnik? select count (*) from vlastnik group by id_vla;
Protina se zamyslena stavba s moji parcelou?
select SDO_GEOM.SDO_INTERSECTION(p1.tvarparcela, p2.tvarparcela, 0.005) from parcela p1, parcela p2 where p1.cisloparcely = 258 AND p2.cisloparcely = 111; -- p1 = 253 je moje parcela, p2 = 111 stavba
U kterých parcel dojde k záboru půdy při stavbě?
select sdo_geom.sdo_intersection(p.tvarparcely,p2.tvarparcely,0.005) from parcela p, parcela p2 where p2.cisloparcela=111;