PRACA DYPLOMOWA - AGH
PRACA DYPLOMOWA - AGH
PRACA DYPLOMOWA - AGH
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Akademia Górniczo-Hutnicza<br />
im. Stanisława Staszica w Krakowie<br />
_________________________<br />
Wydział Geodezji Górniczej<br />
i Inżynierii Środowiska<br />
Katedra Geoinformacji, Fotogrametrii i Teledetekcji Środowiska<br />
<strong>PRACA</strong> <strong>DYPLOMOWA</strong><br />
„Wielokryterialna analiza GIS<br />
lokalizacji ośrodków narciarskich w okolicach<br />
Krakowa<br />
z uwzględnieniem aspektów środowiskowych”<br />
„Multicriteria GIS analysis<br />
of the ski resorts location nearby the Cracow City<br />
according to the environmental aspects”<br />
Imię i nazwisko: Dorota Żwawa<br />
Kierunek studiów: Inżynieria Środowiska<br />
Specjalność: Monitoring Środowiska<br />
Ocena:<br />
Recenzent<br />
dr inż. Wojciech Drzewiecki<br />
Opiekun pracy<br />
dr inż. Tomasz Pirowski<br />
Oświadczam, świadoma odpowiedzialności karnej za poświadczenie nieprawdy, że<br />
niniejszą pracę dyplomową wykonałam osobiście i samodzielnie i że nie<br />
korzystałam ze źródeł innych niż wymienione w pracy.<br />
.......................................................................<br />
Czytelny podpis autora pracy<br />
1<br />
Kraków 2011
Podziękowania dla<br />
Pana dr inż. Tomasza Pirowskiego<br />
za poświęcony czas i cenne wskazówki<br />
przy realizacji pracy.<br />
2
Streszczenie<br />
W pracy zaprezentowano możliwość wykorzystania Systemów Informacji<br />
Geograficznej (GIS) do wspomagania podejmowania decyzji, a mianowicie do wyboru<br />
najlepszej lokalizacji obiektu narciarskiego. Jako obszar analizy przyjęto teren<br />
stanowiący północną część powiatu krakowskiego oraz niewielkie fragmenty<br />
przylegających do niego powiatów olkuskiego i miechowskiego. Wykorzystano dane<br />
pochodzące m.in. z mapy geośrodowiskowej, glebowo-rolniczej, hydrograficznej oraz<br />
z Bazy Danych Topograficznych (TBD).<br />
W części teoretycznej przybliżono obowiązujące, prawne regulacje dotyczące<br />
pozyskiwania i systematyzowania danych. Przedstawiono strukturę tworzenia informacji<br />
przestrzennej oraz zasoby danych wybranych instytucji. W tej części opisano także<br />
metody decyzyjne m.in. wielokryterialną analizę MCE (Multi-Criteria Evaluation), której<br />
działanie wykorzystano w dalszej części pracy.<br />
W części praktycznej poddano analizie wybrany teren pod kątem wielu czynników<br />
mających wpływ na lokalizację stoku narciarskiego. Opisano etapy przygotowania map<br />
kryteriów, najistotniejsze z nich to mapy zalegania pokrywy śnieżnej, czasu dojazdu,<br />
ekspozycji. Określono ważności poszczególnych czynników na tle pozostałych. Na<br />
podstawie powstałych dwóch zestawów wag wykonano dwie analizy MCE - pierwszą<br />
odzwierciedlającą aspekt ekonomiczny oraz drugą - aspekt turystyczny. W ich wyniku<br />
wyszukano obszary potencjalnie odpowiednie na stoki zjazdowe. Następnie wykonano<br />
trzecią analizę uzupełniającą, dzięki której uzyskano informacje o możliwości<br />
zlokalizowania w pobliżu wyselekcjonowanych terenów obiektów infrastruktury<br />
turystycznej jak parking, restauracja, hotel.<br />
W wyniku przeprowadzonych analiz wyznaczono wstępne lokalizacje spełniające<br />
postawione kryteria. Manualna weryfikacja wyselekcjonowanych obszarów, bazująca na<br />
podkładach map topograficznych, pozwoliła wyeliminować błędne wyniki a także<br />
pozwoliła wyciągnąć wnioski co do usprawnienia sposobu selekcji terenów.<br />
3
Summary<br />
In the dissertation the possibility of using The Systems of Geographical<br />
Information (GIS) is presented. It was used for supporting of the decisions, mainly for<br />
choosing the best location of the ski objects. The area of the analysis was located in the<br />
northern part of Cracow’s district and small parts of Olkusz’s and Miechów’s districts as<br />
well Facts from geo – environmental, soil, hydrographic maps were used, as well as<br />
Datas from Topograpical Base (TBD).<br />
In the theoretical part the legal regulations connected with gaining and<br />
systematizing of datas were described. The structure of making some spatial information<br />
and resources of chosen institutions was described. The decision’s methods such as the<br />
complex analysis MCE (Multi – Criteria Evaluation), which actions were used in<br />
a different part of the dissertation.<br />
In the practical part the described area was analysed according to some factors<br />
connected with the location of the ski object. The stages of preparing maps of different<br />
factors were described. The most important maps were connected with snow coverage,<br />
a time of travelling to that place, the exposure of the place. The factors were described<br />
from the least important to the most important ones. Two analysis (MCE) were done: the<br />
first one connected with an economical aspect and the other one connected with<br />
a tourism. According to these factors the most appropriate areas were found. Another<br />
analysis was done for finding the some information connected with an infrastructure of<br />
tourism such as: hotels, restaurants and car parks.<br />
According to these analysis the locations of the ski objects were found. The<br />
manual selection of the areas, based on some maps, led to the right conclusions and<br />
getting rid of wrong ideas.<br />
4
Spis treści<br />
1. Wstęp ........................................................................................................................... 7<br />
2. GIS w zarządzaniu środowiskiem ............................................................................. 9<br />
2.1. INSPIRE (Infrastructure for Spatial Information in Europe) ............................... 9<br />
2.2. Ustawa o Infrastrukturze Informacji Przestrzennej (IIP) .................................... 12<br />
2.3. Krajowy System Informacji Geograficznej......................................................... 14<br />
2.3.1. Wkład samorządu terytorialnego w tworzenie Krajowego Systemu<br />
Informacji Geograficznej na przykładzie gminy ........................................................ 17<br />
2.3.2. Regionalne systemy informacji przestrzennej ............................................. 18<br />
2.4. GIS w monitoringu i zarządzaniu środowiskiem ................................................ 20<br />
2.4.1. Koncepcja Ekoinfonetu ................................................................................ 22<br />
2.4.2. Ekoinfonet Osady......................................................................................... 23<br />
2.5. Wykorzystanie informacji o środowisku (przykład) – wybór miejsca lokalizacji<br />
biogazowni rolniczej w Gminie Zagórz” ....................................................................... 24<br />
3. GIS jako system wspomagania decyzji ................................................................... 29<br />
3.1. Deterministyczne metody podejmowania decyzji............................................... 29<br />
3.1.1. Analiza wieloparametryczna MCE .............................................................. 30<br />
3.1.2. Ocena wielopodmiotowa MOE.................................................................... 34<br />
3.2. Niedeterministyczne metody podejmowania decyzji .......................................... 34<br />
3.3. AHP ..................................................................................................................... 36<br />
4. Opracowanie projektu GIS...................................................................................... 39<br />
4.1. Przepisy prawne dotyczące lokalizacji wyciągów .............................................. 39<br />
4.2. Charakterystyka terenu badań ............................................................................. 41<br />
4.3. Oprogramowanie ................................................................................................. 42<br />
4.4. Dane .................................................................................................................... 44<br />
4.5. Opracowanie przestrzennej bazy danych ............................................................ 46<br />
5. Wielokryterialna analiza MCE ............................................................................... 49<br />
5.1. Podział kryteriów ................................................................................................ 49<br />
5.2. Opracowanie map czynników ............................................................................. 51<br />
5.2.1. Mapa czasu podróży .................................................................................... 51<br />
5
5.2.2. Mapa zalegania pokrywy śnieżnej ............................................................... 54<br />
5.2.3. Ekspozycja ................................................................................................... 59<br />
5.2.4. Mapa spadków ............................................................................................. 60<br />
5.2.5. Mapa przydatności rolniczej gleb ................................................................ 61<br />
5.2.6. Mapa odległości od dróg .............................................................................. 64<br />
5.2.7. Mapa terenów podmokłych .......................................................................... 67<br />
5.2.8. Mapa form ochrony przyrody ...................................................................... 67<br />
5.2.9. Mapy odległości od formy ochrony przyrody i pomników przyrody. ......... 69<br />
5.2.10. Mapa odległości od zabudowy zwartej ........................................................ 70<br />
5.2.11. Mapa odległości od wód .............................................................................. 71<br />
5.2.12. Mapa odległości od cmentarzy .................................................................... 72<br />
5.3. Opracowanie mapy ograniczeń ........................................................................... 73<br />
5.4. Lokalizacja ośrodka narciarskiego ...................................................................... 77<br />
5.4.1. Wagowanie .................................................................................................. 77<br />
5.4.2. Analiza nr 1 – aspekt ekonomiczny ............................................................. 81<br />
5.4.3. Analiza nr 2 – aspekt atrakcyjności turystycznej ......................................... 86<br />
5.4.4. Napotkane nieprawidłowości, próby ich uniknięcia .................................... 88<br />
5.5. Analiza nr 3 – tereny infrastruktury turystycznej towarzyszącej stokowi<br />
narciarskiemu. ................................................................................................................ 90<br />
5.6. Kompleksowa waloryzacja wyselekcjonowanych obszarów .............................. 92<br />
6. Podsumowanie i wnioski ........................................................................................ 102<br />
7. Literatura ................................................................................................................ 105<br />
Spis tabel i rysunków ..................................................................................................... 109<br />
Załączniki........................................................................................................................ 112<br />
6
1. Wstęp<br />
Obecnie wiele ludzi decyduje się na aktywną formę wypoczynku. Narciarstwo<br />
stało się bardzo popularnym sportem, rośnie więc zapotrzebowanie na obiekty typu<br />
ośrodki narciarskie. Dobrze utrzymany stok o odpowiedniej grubości pokrywy śnieżnej<br />
kusi amatorów tego typu rozrywki, zwłaszcza jeśli jest niedaleko miejsca zamieszkania.<br />
Wyjazd na narty nie musi być wtedy jedynie weekendową rozrywką.<br />
Odpowiednia lokalizacja to uwzględnienie wielu czynników stanowiących<br />
o atrakcyjności bądź opłacalności lokalizowanej inwestycji. Dodatkowo, w przypadku<br />
inwestycji mogących w znaczący sposób oddziaływać na środowisko bądź<br />
społeczeństwo, powstaje problem sprzeczności interesów często kilku grup. Rozwiązanie<br />
takich sporów to duży problem dla jednostek decyzyjnych.<br />
Proces podejmowania decyzji można wspierać wykorzystując narzędzia<br />
zaimplementowane w Systemach Informacji Geograficznych (pol. SIG, ang. GIS).<br />
Istnieją oprogramowania posiadające wbudowane procedury wielokryterialnego<br />
podejmowania decyzji. Jedynym z nich jest IDRISI stworzone przez Clark University.<br />
Program ten wykorzystuje metodę analitycznego procesu hierarchicznego (AHP).<br />
Korzystanie z tego typu oprogramowań usprawnia pracę, wymaga jednak równie<br />
dużego zaangażowania i wiedzy eksperckiej. Jest to bowiem praca o charakterze<br />
interdyscyplinarnym, gdyż uwzględnia wiele różnorodnych czynników (przykładowo:<br />
dotyczące zagospodarowanie terenu, czynniki środowiskowe, kulturowe, ekonomiczne,<br />
techniczne). Korzystanie z systemów GIS niesie z sobą wiele ułatwień, są one jednak<br />
wciąż mało rozpowszechnione. Mimo, iż gminy stanowią podstawę struktury tworzenia<br />
informacji przestrzennej, wciąż w pełni nie korzystają z możliwości jej zastosowania.<br />
Początkowe rozdziały niniejszej pracy opisują zależność poszczególnych<br />
jednostek administracyjnych w budowaniu infrastruktury przestrzennej, przydatność<br />
tworzonych przez nie danych które dzięki obowiązującej ustawie o infrastrukturze<br />
przestrzennej mają charakter interoperacyjny. Opisana została również różnorodność<br />
magazynowanych danych o środowisku pochodzących z podsystemów Państwowego<br />
Monitoringu Środowiska, planowanie ich zestandaryzowania i utworzenia geoportali<br />
7
tematycznych. Są to niezwykle istotne informacje, gdyż odpowiedni zasób danych to<br />
pierwszy krok do podjęcia decyzji.<br />
Część teoretyczna zawiera również opis metod decyzyjnych. W zależności od<br />
warunków podejmowania decyzji rozróżnia się m.in. deterministyczne, do których<br />
należy wielokryterialna analiza MCE. Metoda ta polega na wyselekcjonowaniu<br />
obszarów, które spełniają wszystkie postawione twarde kryteria (ograniczenia).<br />
Dodatkowo w wyniku łączenie miękkich kryteriów otrzymuje się ocenę terenu w skali<br />
przydatności. W analizie tej kryteria mogą być jednakowo ważne, ale jest możliwe<br />
również ich wagowanie.<br />
Praktycznym celem pracy było wyszukanie najlepszej lokalizacji dla ośrodka<br />
narciarskiego przy wykorzystaniu metody wielokryterialnej analizy podejmowania<br />
decyzji (MCE). Analizowanym terenem jest północna część powiatu krakowskiego oraz<br />
mały południowo-wschodni fragment powiatu olkuskiego i południowo-zachodni<br />
powiatu miechowskiego.<br />
Wykonano dwie analizy MCE uwzględniające dwa różne aspekty. Pierwszy<br />
czysto ekonomiczny oraz drugi w większym stopniu uwzględniający atrakcyjność<br />
turystyczną terenu. W celu wygenerowania wag kryteriów zastosowano metodę<br />
porównywania parami (AHP). Obie analizy uwzględniają kryterium ukształtowania<br />
rzeźby terenu jako najważniejszy z perspektywy możliwości zrealizowania inwestycji.<br />
Dodatkowo uwzględniają szereg czynników środowiskowych, jak choćby przydatność<br />
rolnicza gleb, obecność terenów chronionych. W sumie zestaw kryteriów dla obu<br />
wariantów składa się z tych samych 14 czynników oraz 10 barier.<br />
Wykonano jeszcze jedną analizę MCE w celu znalezienia terenów odpowiednich<br />
dla infrastruktury turystycznej towarzyszącej stokowi narciarskiemu. Uwzględnia ona<br />
jedynie kilka kluczowych dla niej kryteriów ze zmianami ich wag w stosunku do<br />
poprzednich analiz.<br />
W końcowej części pracy zestawiono najważniejsze informacje związane<br />
z trudnościami napotkanymi na etapie selekcji obszarów (uznanych za przydatne<br />
i posiadających wymaganą powierzchnię). Zwieńczeniem pracy jest kompleksowy opis<br />
wybranych terenów,<br />
8
2. GIS w zarządzaniu środowiskiem<br />
Coraz większe zapotrzebowanie w zakresie dostępu do informacji przestrzennych,<br />
wykorzystania ich w procesach decyzyjnych oraz rosnąca liczba danych odniesionych do<br />
różnych zjawisk przyrodniczych, antropogenicznych wymusiły budowę infrastruktury<br />
informacji przestrzennej.<br />
2.1. INSPIRE (Infrastructure for Spatial Information in Europe)<br />
Zmiany klimatyczne sprawiły, że zmienił się sposób prowadzenia polityki przez<br />
kraje wysoko rozwinięte, podstawowym paradygmatem stał się „zrównoważony rozwój”.<br />
Taki model gospodarki świadomie kształtuje relacje między rozwojem gospodarczym,<br />
dbałością o środowisko przyrodnicze a zdrowiem człowieka. Zrównoważony rozwój<br />
umożliwia zaspokojenie potrzeb obecnego pokolenia bez umniejszania szans przyszłym.<br />
Warunkiem stało się ograniczenie szkodliwego wpływu przemysłu i konsumpcji na stan<br />
środowiska i zasoby przyrodnicze. Wymaga to stałego monitoringu środowiska. Rozwój<br />
infrastruktury informacji przestrzennej związany jest ze środowiskiem przyrodniczym<br />
i koordynacją prac w zakresie jego monitorowania (Gotlib i in., 2007).<br />
W czerwcu 2004 roku Komisja Europejska wystąpiła z propozycją dyrektywy<br />
ustanawiającej ramy prawne dla ustanowienia i działania Infrastruktury Informacji<br />
Przestrzennej w Europie. Dyrektywa INSPIRE została przyjęta przez Parlament<br />
Europejski i Radę Unii 14 marca 2007 roku. Jej celem jest wspieranie wspólnotowej<br />
polityki ochrony środowiska oraz innych działań mogących oddziaływać na środowisko.<br />
Dyrektywa definiuje wymagania jakie powinna spełniać budowana w krajach<br />
członkowskich infrastruktura informacji przestrzennej np. poprzez zapewnienie<br />
mechanizmów interoperacyjności (rys. 2.1). Zastosowanie tych mechanizmów pozwala<br />
na współdzielenie zasobów danych pochodzących z różnych źródeł, wykorzystania ich<br />
przez różnorodnych użytkowników, zarówno prywatnych jak i państwowych. Nie wnika<br />
przy tym w szczegółowe rozwiązania organizacyjne i technologiczne w jej tworzeniu<br />
(http://www.akademiainspire.pl/ [2]).<br />
9
Rozproszone,<br />
nieujednolicone<br />
bazy danych<br />
Zestandaryzowane<br />
bazy danych<br />
Zestandaryzowana<br />
publikacja danych<br />
Ogólnoeuropejski<br />
dostęp do danych<br />
Rys. 2.1 Schemat działania Dyrektywy (Jarząbek, 2011).<br />
Nowy akt prawny obliguje państwa członkowskie do stworzenia geoportali.<br />
Celem ich tworzenia jest udostępnienie zasobów danych przestrzennych w formie usług<br />
sieciowych, metadanych oraz serwisów katalogowych. Wprowadzenie przepisów<br />
Dyrektywy ma na celu ułatwienie i przyspieszenie dostępu do danych, zwiększenie<br />
możliwości ich porównywania wśród krajów Unii, spowoduje wzrost zainteresowania<br />
danymi przestrzennymi jako narzędziem do zarządzania szeroko pojętym środowiskiem.<br />
Kraje członkowskie powinny utworzyć i obsługiwać (http://www.geoportal.gov.pl/ [10]):<br />
• usługi wyszukiwania, umożliwiające wyszukiwanie zbiorów oraz usług danych<br />
przestrzennych na podstawie zawartości odpowiadających im metadanych oraz<br />
umożliwiające wyświetlanie zawartości metadanych;<br />
• usługi przeglądania, umożliwiające co najmniej: wyświetlanie, nawigowanie,<br />
powiększanie i pomniejszanie, przesuwanie lub nakładanie na siebie zbiorów<br />
danych przestrzennych oraz wyświetlanie informacji z legendy i wszelkiej istotnej<br />
zawartości metadanych;<br />
• usługi pobierania, umożliwiające pobieranie kopii całych zbiorów danych<br />
przestrzennych lub części takich zbiorów oraz, gdy jest to wykonalne, dostęp<br />
bezpośredni;<br />
• usługi przekształcania, umożliwiające przekształcenie zbiorów danych<br />
przestrzennych w celu osiągnięcia interoperacyjności;<br />
• usługi umożliwiające uruchamianie usług danych przestrzennych.<br />
Główne założenia INSPIRE (http://www.akademiainspire.pl/ [2]):<br />
• dane powinny być pozyskiwane tylko jeden raz oraz przechowywane i zarządzane<br />
w sposób najbardziej poprawny i efektywny przez odpowiednie instytucje<br />
i służby;<br />
• powinna być zapewniona ciągłość przestrzenna danych tak, aby było możliwe<br />
pozyskanie różnych zasobów, z różnych źródeł oraz aby możliwe było ich<br />
udostępnianie wielu użytkownikom i do różnorodnych zastosowań;<br />
10
• dane przestrzenne powinny być przechowywane na odpowiednim (jednym)<br />
poziomie administracji publicznej i udostępniane podmiotom na wszystkich<br />
pozostałych poziomach;<br />
• dane przestrzenne niezbędne do odpowiedniego zarządzania przestrzenią na<br />
wszystkich poziomach administracji publicznej powinny być powszechnie<br />
dostępne (tj. bez warunków ograniczających i/lub utrudniających ich swobodne<br />
wykorzystanie);<br />
• powinien być zapewniony dostęp do informacji o tym, jakie dane przestrzenne są<br />
dostępne i na jakich warunkach, a także informacja umożliwiająca użytkownikowi<br />
ocenę przydatności tych danych do swoich celów.<br />
Początkowo INSPIRE miała być ściśle ukierunkowane na środowisko<br />
przyrodnicze (Gotlib i in., 2007), rozszerzono jednak zakres tematyczny, co w istotny<br />
sposób zwiększyło funkcjonalność systemu. Zbiory danych infrastruktury reprezentują<br />
bogatą treść informacyjną podzieloną na 34 tematy zawarte w trzech aneksach ww.<br />
Dyrektywy. Obejmują one następujące zagadnienia:<br />
Aneks I:<br />
1. Systemy odniesienia za pomocą współrzędnych<br />
2. Systemy siatek geograficznych<br />
3. Nazwy geograficzne<br />
4. Jednostki administracyjne<br />
5. Adresy<br />
6. Działki katastralne<br />
7. Sieci transportowe<br />
8. Hydrografia<br />
9. Obszary chronione<br />
Aneks II:<br />
1. Ukształtowanie terenu<br />
2. Użytkowanie terenu<br />
3. Ortofotomapa<br />
4. Geologia<br />
Aneks III:<br />
1. Jednostki statystyczne<br />
2. Budynki<br />
3. Gleba<br />
4. Zagospodarowanie przestrzenne<br />
11
5. Zdrowie i bezpieczeństwo ludzi<br />
6. Usługi użyteczności publicznej i służby państwowe<br />
7. Urządzenia do monitorowania środowiska<br />
8. Obiekty produkcyjne i przemysłowe<br />
9. Obiekty rolnicze i akwakultury<br />
10.Rozmieszczenie ludności – demografia<br />
11. Gospodarowanie obszarem/strefy ograniczone/regulacyjne<br />
i jednostki sprawozdawcze<br />
12. Strefy zagrożenia naturalnego<br />
13. Warunki atmosferyczne<br />
14. Warunki meteorologiczno-geograficzne<br />
15. Warunki oceanograficzno-geograficzne<br />
16. Regiony morskie<br />
17. Regiony biogoegraficzne<br />
18. Siedliska i obszary przyrodniczo jednorodne<br />
19. Rozmieszczenie gatunków<br />
20. Zasoby energetyczne<br />
21. Zasoby mineralne<br />
Tak szeroki zakres informacji wskazuje na możliwość wszechstronnego<br />
zastosowania infrastruktury do celów gospodarczych, społecznych i politycznych.<br />
2.2. Ustawa o Infrastrukturze Informacji Przestrzennej (IIP)<br />
Dyrektywa INSPIRE zobligowała do ustanowienia krajowej infrastruktury<br />
informacji przestrzennej wchodzącej w skład infrastruktury europejskiej wszystkich 27<br />
krajów członkowskich.<br />
Ustawa o Infrastrukturze Informacji Przestrzenne z dnia 4 marca 2010r. jest<br />
osobnym aktem prawnym, zmieniającym jednocześnie dotychczas obowiązujące ustawy.<br />
Zmiany dotyczą ustawy z dnia 17 maja 1989 r. — Prawo geodezyjne i kartograficzne,<br />
ustawy z dnia 4 lutego 1994 r. — Prawo geologiczne i górnicze, ustawy z dnia 29<br />
czerwca 1995 r. o statystyce publicznej, ustawy z dnia 27 kwietnia 2001 r. — Prawo<br />
ochrony środowiska, ustawy z dnia 16 kwietnia 2004 r. o ochronie przyrody (zmiany<br />
zawarte w art. 23 – art. 28 IIP). Nowelizacja Prawa Geodezyjnego i Kartograficznego<br />
w szczególności miała na celu dostosowanie przepisów tej ustawy do potrzeb związanych<br />
z tworzeniem w Polsce infrastruktury informacji przestrzennej<br />
(http://www.bip.mswia.gov.pl/ [6]).<br />
„Celem ustawy IIP jest dostosowana do potrzeb i łatwo dostępna informacja<br />
przestrzenna dla każdego” (Gaździcki, 2009). Jej zadaniem jest również optymalizacja<br />
pozyskiwanych danych przestrzennych oraz utrzymanie tych zasobów przez<br />
12
administrację publiczną. Koordynatorem działań związanych z tworzeniem,<br />
utrzymywaniem i rozwijaniem infrastruktury jest minister właściwy ds. administracji<br />
publicznej (rys. 2.2). Jego zadaniem jest również informowanie Komisji Europejskiej<br />
o postępach tworzenia i funkcjonowania infrastruktury. Główny Geodeta Kraju<br />
wykonuje określone zadania koordynacyjne zawarte w art. 19 IIP. Organem<br />
opiniodawczym i doradczym jest Rada Infrastruktury Przestrzennej (art. 22). Dla<br />
każdego tematu danych przestrzennych (zawartych w załączniku do ustawy, wyznaczone<br />
przez INSPIRE) przewidziany jest organ wiodący. Może nim być minister lub centralny<br />
organ administracji, którego zadaniem jest koordynowanie prac w zakresie określonego<br />
przez ustawę tematu.<br />
„Infrastruktura jest tworzona, utrzymywana i rozwijana, a także funkcjonuje<br />
w wyniku współdziałania współtworzących ją organów administracji oraz osób trzecich”<br />
(art. 17 ust. 1 IIP).<br />
Rys. 2.2 Struktura koordynacyjna (http://www.akademiainspire.pl/[1]).<br />
13
Podsumowując, ustawa definiuje zasady tworzenia infrastruktury, wyznacza<br />
organy właściwe w tej sprawie, definiuje zasady współpracy i koordynacji<br />
funkcjonowania infrastruktury. Określa również zasady korzystania z niej (dostęp do<br />
zbiorów danych przestrzennych, form udostępniania danych).<br />
2.3. Krajowy System Informacji Geograficznej<br />
Adresatem Dyrektywy INSPIRE jest administracja publiczna wszystkich szczebli<br />
(rys. 2.3), toteż informacje na przestrzeni kraju zbierane są na różnych poziomach<br />
i o różnej szczegółowości tworząc Krajowy System Informacji Geograficznej (KSIG).<br />
Europa<br />
Polska<br />
Województwo<br />
Powiat<br />
Gmina<br />
Rys. 2.3 Struktura tworzenia geoinformacji (Jarząbek, 2011)<br />
INSPIRE nie jest nowym programem zbierania danych. Zakłada on optymalizację<br />
zakresu wykorzystywania już dostępnych, poprzez określenie wymagań odnośnie<br />
dokumentacji istniejących danych, wdrożenia usług ukierunkowanych na łatwiejszy do<br />
nich dostęp (Preuss, 2004a). Krajowy system tworzony na trzech poziomach powstaje<br />
w oparciu o dane z różnych źródeł (Preuss, 2004b):<br />
• krajowy – nadzorowany przez Głównego Geodetę Kraju w oparciu o Bazę<br />
Danych Ogólnogeograficznych w skali 1: 250 000;<br />
• wojewódzki - przez marszałków województw w oparciu o TBD w skali 1:10 000<br />
i VMap 2 w skali 1:50 000;<br />
• powiatowy - przez starostów w oparciu o ewidencję gruntów.<br />
Baza Danych Ogólnogeograficznych BDO zawiera podstawowe informacje<br />
zapisane z dokładnością właściwą skali 1:250 000. Są one uporządkowane w ośmiu<br />
tematach podzielonych na dwadzieścia warstw informacyjnych<br />
14
(http://www.codgik.gov.pl/ [8]). Baza ta powstała w oparciu o wiele różnych materiałów<br />
źródłowych, wśród których najistotniejszym był zasób bazy VMap poziomu 1,<br />
udostępniony przez Zarząd Geografii Wojskowej Sztabu Generalnego Wojska Polskiego.<br />
Celem tej bazy jest umożliwienie prezentacji obiektów i zjawisk w skali kraju przez różne<br />
instytucje i urzędy (Osada, 2009).<br />
Baza Danych Topograficznych TBD zawiera informacje odpowiadające<br />
szczegółowości mapy topograficznej w skali 1:10 000. Jej treść obejmuje 10 obszarów<br />
tematycznych, z których każdy zapisany jest w kilku warstwach (jednostki podziału<br />
administracyjnego, sieci dróg i kolei, budowle i urządzenia, kompleksy pokrycia terenu,<br />
kompleksy użytkowania terenu, sieci cieków, tereny chronione, osnowa, sieci uzbrojenia<br />
terenu, punkty adresowe). Podzielona jest na następujące komponenty: wektorowa baza<br />
danych topograficznych, mapa cyfrowa, baza Ortofotomapa cyfrowych, numeryczny<br />
model rzeźby terenu (Osada, 2009). Producentami składowych Bazy Danych<br />
Topograficznych są firmy komercyjne, zwykle geodezyjne a dysponentem danych -<br />
poprzez CODGiK oraz WODGiK, jest GUGIK (Osada, 2009).<br />
Idąc dalej w kierunku integracji baz danych tworzonych w skali kraju,<br />
wykorzystano wojskowe mapy wektorowe. Vmap L2, czyli wojskowa mapa poziomu<br />
drugiego (odpowiadająca skali 1:50 000) obejmuje obszar całego kraju, opracowana<br />
została w formacie VPF, który jest standardem NATO. Treść tego produktu obejmuje 110<br />
klas obiektów zgrupowanych w dziewięciu użytkowych warstwach tematycznych:<br />
granice administracyjne, rzeźba terenu, fizjografia, transport, zabudowa, hydrografia,<br />
roślinność, infrastruktura obsługi ruchu powietrznego i przemysł (Bielecka, 2006).<br />
Krajowy system uzupełniają również mapy topograficzne w skalach 1:50 000<br />
i 1:10 000, mapa zasadnicza, ortofotomapa oraz dane z Ewidencji Gruntów i Budynków<br />
(EGiB) (rys. 2.4).<br />
15
Rys. 2.4 Składniki Krajowego Systemu Informacji Geograficznej (Brożek, 2011).<br />
Rysunek 2.5 przedstawia ideę interoperacyjności pomiędzy różnymi zasobami<br />
danych przestrzennych prowadzonymi na poziomie krajowym, wojewódzkim<br />
i powiatowym w ramach Krajowego Systemu Informacji o Terenie. Rysunek obrazuje jak<br />
np. budynek z bazy EGiB przechodzi do bazy TBD a budynek z bazy TBD może zasilić<br />
bazę danych ogólnogeograficznych. Przedstawiony jest również przepływ informacji<br />
w drugą stronę. Zasoby różnych baz uzupełniają się dając spójny wynik (Osada, 2009).<br />
Rys. 2.5 Idea interoperacyjności pomiędzy różnymi danym prowadzonych na różnych<br />
poziomach (Osada, 2009).<br />
16
2.3.1. Wkład samorządu terytorialnego w tworzenie Krajowego Systemu<br />
Informacji Geograficznej na przykładzie gminy<br />
Nowelizowane przepisy art. 47a i art. 47b wypełniają lukę prawną dotyczącą<br />
treści i zasad prowadzenia przez gminy ewidencji numeracji porządkowej nieruchomości,<br />
która jest jednym z referencyjnych rejestrów krajowej infrastruktury informacji<br />
przestrzennej (http://www.bip.mswia.gov.pl/ [6]).<br />
Rolę katastru w Polsce pełni Ewidencja Gruntów i Budynków (EGiB)<br />
w powiązaniu z systemem ksiąg wieczystych. W EGiB zapisana jest identyfikacja<br />
nieruchomości w przestrzeni, w księgach wieczystych prawa rzeczowe do nieruchomości.<br />
Powiązaniem między tymi systemami jest numer księgi wieczystej, zapisany jako atrybut<br />
działek ewidencyjnych w EGiB, i identyfikator geodezyjny działki zapisany w księdze<br />
wieczystej (Gotlib i in., 2007).<br />
Prowadzenie EGiB jest zadaniem rządowym powierzonym samorządom (Prawo<br />
Geodezyjne i Kartograficzne, art. 22). Z reguły EGiB prowadzą władze samorządowe<br />
powiatu, czasem gminy. Kontrolę nad EGiB sprawuje Główny Geodeta Kraju przez<br />
wojewódzkich inspektorów Nadzoru Geodezyjnego i Kartograficznego<br />
(http://www.zcpwz.pl/ [18]).<br />
Najcenniejszym elementem każdego systemu informacji o terenie są dane. Dane<br />
EGiB pokrywające całą powierzchnię kraju są znakomitym punktem wyjścia do<br />
tworzenia takich systemów. Szczególną rolę odgrywają w systemach informacji o terenie<br />
na szczeblu powiatowym i miejskim, dla którego są naturalną warstwą referencyjną.<br />
Dane z tej ewidencji są wykorzystywane przez urzędy powiatowe i gminne, stanowią<br />
podstawę do planowania przestrzennego, wydawania decyzji administracyjnych<br />
i zarządzania infrastrukturą.<br />
W urzędach powiatowych z danych EGiB korzystają m.in. (Gotlib i in., 2007):<br />
• wydziały urbanistyki, architektury i budownictwa do ustalania stron postępowania<br />
przy wydawaniu pozwoleń na budowę, tzn. właścicieli działek graniczących<br />
z działką na której ma być prowadzona budowa;<br />
• wydziały gospodarki nieruchomościami i mieniem powiatu do stworzenia bazy<br />
danych o nieruchomościach Skarbu Państwa oraz bazy danych<br />
o nieruchomościach stanowiących własność powiatu;<br />
17
• wydziały ochrony środowiska, rolnictwa i leśnictwa do wydawania decyzji<br />
o wyłączeniu gruntu w produkcji rolnej i naliczaniu opłat za to wyłączenie,<br />
pozwoleń wodnoprawnych oraz zezwoleń na wycięcie drzew;<br />
• dane z EGiB wykorzystywane są przy kontrolach terenowych prowadzonych<br />
z urzędu lub w odpowiedzi na skargi ludności, do planowania lokalizacji<br />
inwestycji ważnych dla środowiska.<br />
W urzędach miast i gmin wykorzystywane są dane EGiB przy (Gotlib i in., 2007):<br />
• naliczaniu podatków od nieruchomości, podatku rolnego i podatku leśnego;<br />
• naliczaniu opłat za wieczyste użytkowanie gruntów komunalnych;<br />
• tworzeniu bazy danych o nieruchomościach tworzących zasób gminny (rejestr<br />
mienia komunalnego);<br />
• regulacji stanu prawnego;<br />
• opracowaniu danych na potrzeby statystyki;<br />
• wydawaniu decyzji o warunkach zabudowy;<br />
• śledzeniu zmian w terenie;<br />
• określaniu nowych terenów pod inwestycje;<br />
• tworzeniu map tematycznych (drogi gminne, uzbrojenie terenu).<br />
2.3.2. Regionalne systemy informacji przestrzennej<br />
Regionalny System Informacji Przestrzennej jest częścią składową Krajowego<br />
Systemu Informacji o Terenie. Tworzenie ich ma na celu głównie optymalizację<br />
procesów decyzyjnych administracji publicznej wszystkich szczebli, zwiększenie<br />
skuteczności ochrony i monitorowania środowiska oraz usprawnienie procesów<br />
planistycznych. Ważnym aspektem jest również usprawnienie działania służb zarządzania<br />
kryzysowego. Dodatkowo możliwym staje się wykorzystanie baz danych w dziedzinach<br />
(Gotlib i in., 2007):<br />
• ratownictwo medyczne;<br />
• statystyka publiczna;<br />
• projektowanie i utrzymywanie sieci uzbrojenia terenu i urządzeń technicznych;<br />
• nadzorowanie stanu płatności za dostarczone usługi użytkownikom;<br />
• ujawnianie prawidłowości zjawisk przyrodniczych i społeczno-gospodarczych;<br />
18
• rozwój geobiznesu i geomarketingu;<br />
• przygotowanie opracowań fizjograficznych w obrębie planowanych przebiegów<br />
nowych i modernizowanych tras komunikacyjnych;<br />
• przygotowanie wariantów koncepcyjnych lokalizacji inwestycji.<br />
Systemy informacji przestrzennej tworzone w skali lokalnej ułatwiają władzom<br />
realizację wielu zadań ustawowych oraz ułatwiają życie mieszkańcom, dlatego wiele<br />
województw, powiatów i miast opracowało portale geoinformacyjne. Zazwyczaj są one<br />
tworzone przez specjalistyczne firmy GIS we współpracy z ośrodkami dokumentacji<br />
geodezyjno-kartograficznej (Gotlib i in., 2007).<br />
Przykłady regionalnych systemów informacji przestrzennej:<br />
• Regionalny System Informacji Przestrzennej Województwa Łódzkiego;<br />
• Mazowiecki System Informacji Przestrzennej;<br />
• Regionalny System Informacji Przestrzennej Województwa Śląskiego;<br />
• Podlaski System Informacji Przestrzennej;<br />
• Kujawsko-Pomorski System Informacji Przestrzennej;<br />
• Dolnośląski System Informacji Przestrzennej;<br />
• Opolski System Informacji Przestrzennej;<br />
• Wrota Małopolski.<br />
Wrota Małopolski<br />
Jednym z bardziej znanych wojewódzkich serwisów internetowych są „Wrota<br />
Małopolski”. Celem jego realizacji było zgromadzenie różnorodnych zasobów<br />
i udostępnienie ich w jednolity zintegrowany sposób. Głównymi modułami tematycznymi<br />
portalu są: „aktualności”, „wizytówka Małopolski”, „edukacja”, „kultura”,<br />
„niepełnosprawni”, „praca”, „przedsiębiorczość”, „środowisko i obszary wiejskie”,<br />
„turystyka”, „zdrowie”, „mapy GIS-GPS”. Zasoby kartograficzne udostępniane<br />
w ostatnim z wymienionych modułów pochodzą ze zbiorów Wojewódzkiego Ośrodka<br />
Dokumentacji Geodezyjnej i Kartograficznej. System udostępniono aby upowszechnić<br />
informacje geograficzne o województwie małopolskim (http://mapy.wrotamalopolski.pl/<br />
[12]).<br />
Moduł „mapy GIS-GPS” składa się z kilku elementów. Najważniejszymi są<br />
aplikacje „mapa interaktywna” i „Małopolska z lotu ptaka” umożliwiające szybkie<br />
19
przeglądanie bardzo dużych zasobów danych rastrowych. Istnieje możliwość<br />
dynamicznego przesuwania obrazu, zmiany skali, wyszukiwanie zdefiniowanych<br />
punktów (kościoły, skrzyżowania, urzędy gmin, szkoły).<br />
Serwis udostępnia trzy rodzaje map: klasyczne mapy topograficzne opracowane<br />
w skali 1:10000, mapy fotograficzne wykonane na podstawie zdjęć lotniczych oraz mapa<br />
historyczna całego województwa. Pomijając mapę historyczną aktualność map jest<br />
zróżnicowana, najstarszymi są dwukolorowe mapy topograficzne (1985-95), mapy<br />
topograficzne wielokolorowe i mapy fotograficzne wykonane w ostatnich dziesięciu<br />
latach (Gotlib i in., 2007).<br />
2.4. GIS w monitoringu i zarządzaniu środowiskiem<br />
„Komisja Europejska do tej pory wdrożyła już wiele dyrektyw mających na celu<br />
poprawę stanu środowiska naturalnego. Na przykład ESE Natura 2000 działająca<br />
w sposób bezpośredni, bądź dyrektywa 90/313/EWG w sposób pośredni poprzez<br />
zagwarantowanie dostępu do informacji o środowisku. Wejście w życie dyrektywy<br />
INSPIRE miało na celu stworzenie jednolitej infrastruktury danych przestrzennych,<br />
stanowiącej podstawę monitoringu i oceny stanu środowiska naturalnego, a w<br />
konsekwencji źródła informacji dla potrzeb kreowania polityki UE w tym zakresie”<br />
(Wężyk i in., 2009).<br />
Organem utworzonym w celu zapewnienia wiarygodnych informacji o stanie<br />
środowiska jest Państwowy Monitoring Środowiska (PMŚ). Ustawa Prawo ochrony<br />
środowiska (P.o.ś.) definiuje PMŚ jako obejmujący nie tylko diagnozę stanu środowiska,<br />
ale także jego prognozę, nałożyła również obowiązek systematycznego gromadzenia,<br />
przetwarzania i rozpowszechniania informacji o środowisku (art. 25 ust. 1 i 2).<br />
System informatyczny funkcjonujący w Inspektoracie Ochrony Środowiska na<br />
potrzeby Państwowego Monitoringu Środowiska w Województwie Małopolskim składa<br />
się z kilku podsystemów związanych z poszczególnymi komponentami środowiska.<br />
Systemy te rejestrują, przetwarzają i udostępniają dane o środowisku przy użyciu technik<br />
informatycznych.<br />
W województwie małopolskim w ramach PMŚ funkcjonują bazy danych<br />
(http://www.krakow.pios.gov.pl/ [11]):<br />
• w podsystemie monitoringu powietrza:<br />
20
system zbierania danych z monitoringu automatycznego powietrza oparty na<br />
oprogramowaniu XR francuskiej firmy ISEO,<br />
JPOAT – baza danych jakości powietrza, służy do przekazywania zgodnie<br />
z wymaganiami ustawowymi informacji o stacjach i stanowiskach<br />
pomiarowych oraz wynikach z pomiarów,<br />
• w podsystemie monitoringu wód<br />
JAWO – baza danych o jakości wód powierzchniowych płynących,<br />
• w podsystemie monitoringu hałasu<br />
rejestr wyników badań klimatu akustycznego (arkusze w formacie Excel),<br />
• w podsystemie monitoringu pól elektromagnetycznych<br />
rejestr wyników badań pól elektromagnetycznych (arkusze w formacie Excel).<br />
Dla potrzeb PMŚ gromadzone są informacje również w innych bazach<br />
funkcjonujących w Inspektoracie (http://www.krakow.pios.gov.pl/ [11]):<br />
• w Wojewódzkiej Bazie Informacji o Środowisku (WBIoŚ) funkcjonującej<br />
w systemie EKOINFONET, realizowanej zgodnie z obowiązkiem określonym<br />
w art. 286a ustawy P.o.ś. Baza danych prowadzona jest w oparciu o dane<br />
o emisjach zanieczyszczeń do powietrza, poborze wód oraz odprowadzaniu<br />
ścieków, ujętych w formularzach przesyłanych przez podmioty korzystające ze<br />
środowiska,<br />
• w bazie składowisk i spalarni – rejestrowane są wszystkie zmiany stanu<br />
ilościowego oraz prawnego składowisk przemysłowych i komunalnych oraz dane<br />
o instalacjach do termicznego przekształcania odpadów.<br />
Obecnie funkcjonujące bazy danych nie spełniają wymogów sprawozdawczości<br />
międzynarodowych, wobec czego Główny Inspektorat Ochrony Środowiska (GIOŚ)<br />
planuje budowę nowych baz danych PMŚ w oparciu o aktualizację koncepcji<br />
SI EKOINFONET, opracowanie koncepcji Systemu Informacji Przestrzennej (GIS) dla<br />
SI EKOINFONET oraz wykonanie projektów wstępnych czterech baz danych PMŚ:<br />
• monitoringu powietrza;<br />
• monitoringu hałasu;<br />
• monitoringu wód powierzchniowych (rzeki);<br />
• monitoringu promieniowania niejonizującego.<br />
21
Program Państwowego Monitoringu Środowiska województwa małopolskiego na<br />
lata 2010-2012 przewiduje zestandaryzowanie informacji przestrzennej uwzględniającej<br />
interoperacyjność i harmonizację danych przestrzennych zgodnie z prawodawstwem<br />
krajowym i aktami wykonawczymi KE do dyrektywy 2007/2/WE ustanawiającej<br />
infrastrukturę informacji przestrzennej we Wspólnocie (INSPIRE). Dotychczas<br />
tematyczny geoportal funkcjonuje jedynie w obszarze osadów.<br />
2.4.1. Koncepcja Ekoinfonetu<br />
Koncepcja działania systemu Ekoinfonet polega na doskonaleniu form i sposobów<br />
gromadzenia, przetwarzania i wizualizacji danych o stanie środowiska (rys. 2.6)<br />
Rys. 2.6 Schemat działania Ekoinfonetu (Kraszewski, 2002).<br />
„Koncepcja rysuje zamysł systemu baz danych tworzących Krajowe<br />
Repozytorium Danych Inspekcji Ochrony Środowiska (KRDIOŚ), gdzie byłyby<br />
archiwizowane wszystkie wyniki pomiarów i obserwacji, jakie są zbierane<br />
w podsystemach PMŚ, jak również informacje dotyczące pozostałych dziedzin<br />
działalności Inspekcji” (Kraszewski, 2003).<br />
Rysunek 2.7 przedstawia schemat zbierania danych pochodzących ze stacji<br />
automatycznych pomiarów, z pomiarów manualnych, a także już istniejących<br />
22
magazynowanych w archiwach różnych instytucji i uczelni. Zebrane dane zostają<br />
ujednolicone w jednej bazie danych, przez co stają się szerzej i częściej wykorzystywane.<br />
Rys. 2.7 Podsystem zbierania danych, wprowadzenie danych do KRDIOŚ<br />
(Kraszewski, 2002).<br />
2.4.2. Ekoinfonet Osady<br />
W ramach podsystemu Państwowego Monitoringu Środowiska – Monitoring<br />
jakości śródlądowych wód powierzchniowych serwis internetowy Osady prezentuje dane<br />
dotyczące badań osadów wodnych rzek i jezior. Projekt został zrealizowany przez<br />
Państwowy Instytut Geologiczny (obecnie Państwowy Instytut Badawczy) dla Głównego<br />
Inspektoratu Ochrony Środowiska. Bezpośredni nadzór nad realizacją programu badań<br />
sprawuje Departament Monitoringu GIOŚ. Serwis składa się z dwóch sekcji<br />
(http://ekoinfonet.gios.gov.pl/):<br />
• „Wyniki badań” w formie tabeli dla wybranego roku prezentuje informacje<br />
o zanieczyszczeniach oraz wynikach ocen;<br />
• „Mapa”- aktywna przeglądarka danych przestrzennych pozwala na wyszukanie<br />
i dostęp do atrybutów monitorowanego punktu oraz wizualizację na tle<br />
wybranych warstw referencyjnych.<br />
23
Zgodnie z dyrektywą INSPIRE dane prezentowane w serwisie dostępne są<br />
również w postaci usług WMS i WFS.<br />
Poniższy rysunek 2.8 przedstawia zakres udostępnionych informacji<br />
o przykładowym monitorowanym punkcie na rzece Dunajec. Punkt zaprezentowano na<br />
warstwie bazowej - Mapa Wektorowa Poziomu 0. Pozostałymi dostępnymi warstwami<br />
referencyjnymi są: Państwowy Rejestr Granic oraz Granice Regionalnych Zarządów<br />
Gospodarki Wodnej.<br />
Rys. 2.8 Ekoinfonet Osady – przykład geoportalu tematycznego zawierającego dane<br />
z jednego z podsystemów PMŚ (zakres udostępnianych danych).<br />
2.5. Wykorzystanie informacji o środowisku (przykład) – wybór<br />
miejsca lokalizacji biogazowni rolniczej w Gminie Zagórz”<br />
Większość decyzji w gminach podejmowanych jest na podstawie informacji<br />
przestrzennych. Podjęcie decyzji wymaga zestawienia ze sobą informacji z różnych map.<br />
Najczęściej wykorzystywanymi informacjami przestrzennymi są: mapa ewidencyjna,<br />
mapa zasadnicza, miejscowe plany zagospodarowania przestrzennego, studium<br />
24
uwarunkowań i kierunków zagospodarowania przestrzennego, mapa topograficzna, mapa<br />
obszarów chronionych, mapa glebowo – rolnicza, ortofotomapa. Przydatnymi są też<br />
informacje o rozmieszczaniu pozwoleń na budowę, rozmieszczeniu decyzji o warunkach<br />
zabudowy oraz decyzji o lokalizacji inwestycji celu publicznego, rozmieszczenie<br />
ludności.<br />
Projekt Gminy Zagórz miał na celu wskazanie odpowiedniej lokalizacji pod<br />
inwestycję biogazowi (http://www.akademiainspire.pl [3]) . Rozważono kilka<br />
wariantów inwestycji spełniających kryteria:<br />
• wybrane miejsce lokalizacji nie będzie ingerowało w obszary chronione,<br />
• dostępność do niezbędnej infrastruktury technicznej,<br />
• dostępność substratów,<br />
• teren przeznaczony w miejscowym planie zagospodarowania przestrzennego pod<br />
prowadzenie inwestycji o charakterze przemysłowym oraz zgodność z kierunkami<br />
polityki przestrzennej Gminy Zagórz,<br />
• w przypadku braku miejscowego planu, teren przewidziany pod inwestycje<br />
powinien spełniać przesłanki niezbędne do wydania decyzji lokalizacyjnej,<br />
• teren lokalizacji inwestycji położony z dala od zabudowy mieszkaniowej,<br />
• planowana inwestycja zlokalizowana będzie na działkach stanowiących własność<br />
gminy.<br />
Wykorzystane dane przestrzenne<br />
Dane nr 1<br />
Nazwa: Ortofotomapa<br />
Źródło: serwis WMS Geoportal.gov.pl<br />
Rodzaj danych: raster<br />
Przeznaczenie: podkład do wizualizacji projektu studium przypadku<br />
Dane nr 2<br />
Nazwa: warstwa działek ewidencyjnych<br />
Źródło: serwis WMS Geoportal.gov.pl<br />
Rodzaj danych: wektor<br />
Przeznaczenie: podkład do analizy granic i pow. działek ewidencyjnych na<br />
których możliwa będzie realizacja inwestycji, usytuowanie przestrzenne<br />
działki oraz podkład wizualizacji projektu studium przypadku<br />
Dane nr 3<br />
Nazwa: dane katastralne<br />
25
Źródło: PODGiK<br />
Rodzaj danych: wektor<br />
Przeznaczenie: ustalenie właścicieli działek niezbędnych do wykupienia<br />
celem realizacji inwestycji oraz dokładnej powierzchni i sposobu<br />
użytkowania gruntów<br />
Dane nr 4<br />
Nazwa: dane adresowe<br />
Źródło: zasób własny<br />
Rodzaj danych: tekstowe<br />
Przeznaczenie: ustalenie adresów działek niezbędnych do wykupienia celem<br />
realizacji inwestycji<br />
Dane nr 5<br />
Nazwa: studium ochrony przeciwpowodziowej<br />
Źródło: Regionalny Zarząd Gospodarki Wodnej w Krakowie<br />
Rodzaj danych: wektor<br />
Przeznaczenie: wizualizacja wyznaczonych obszarów bezpośredniego<br />
zagrożenia powodzią oraz zasięgów stref zalewowych dla wód<br />
powodziowych o prawdopodobieństwie przewyższenia p=0,5%, p=1%<br />
i p=2% tle planowanej inwestycji<br />
Dane nr 6<br />
Nazwa: Studium Uwarunkowań i Kierunki Zagospodarowania Przestrzennego<br />
Gminy Zagórz<br />
Źródło: zasób własny<br />
Rodzaj danych: raster<br />
Przeznaczenie: wizualizacja uwarunkowań i kierunki polityki przestrzennej<br />
dzielnicy Zasław na tle planowanej inwestycji – sprawdzenie, czy rozważane<br />
warianty lokalizacji biogazowi są zgodne z kierunkami – polityki<br />
przestrzennej Gminy Zagórz<br />
Dane nr 7<br />
Nazwa: Miejscowy Plan Zagospodarowania Przestrzennego dla miasta<br />
Zagórz o nazwie „Zastaw-I”<br />
Źródło: zasób własny<br />
Rodzaj danych: raster<br />
Przeznaczenie: ustalenie przeznaczenie terenów objętych planem i ich<br />
wizualizacja w stosunku do rozważanych wariantów lokalizacji biogazowi.<br />
Inwestycja wymagała współpracy z wieloma instytucjami m.in. z Urzędem Miasta<br />
i Gminy Zagórz (Referat Architektury, Gospodarki Nieruchomościami i Ochrony<br />
Środowiska) który opiniował lokalizację pod względem urbanistyczno-architektonicznym<br />
i estetycznym; szkodliwości dla środowiska i zdrowia ludzi. Referat Organizacyjny<br />
i Spraw Obywatelskich udostępniał informacje o danych adresowych właścicieli działek.<br />
26
Do zadań Referatu Rozwoju i Gospodarki Komunalnej należało opiniowanie lokalizacji<br />
w zakresie istniejącej i projektowanej infrastruktury technicznej oraz w zakresie<br />
dostępności terenu inwestycji do dróg publicznych. Regionalny Dyrektor Ochrony<br />
Środowiska w Rzeszowie – Wydział Ocen Oddziaływania na Środowisko to instytucja<br />
odpowiedzialna za obszar Natura 2000, na który inwestycja może potencjalnie znacząco<br />
oddziaływać. Instytucją do której należała odpowiedzialność za utrzymanie dobrego stanu<br />
wód oraz konsultacje w zakresie terenów bezpośredniego zagrożenia powodziowego był<br />
Regionalny Zarząd Gospodarki Wodnej w Krakowie, Zarząd Zlewni Sanu z/s<br />
w Przemyślu.<br />
Spośród czterech wariantów (rys. 2.9) wybrano rozwiązanie najkorzystniejsze dla<br />
środowiska. Wybrany teren cechował się dobrą infrastrukturą komunikacyjną. Wybrane<br />
miejsce lokalizacji nie ingeruje w obszar Natura 2000 i inne obszary chronione, nie miał<br />
przeszkód prawnych zarówno do wydania decyzji lokalizacyjnej jak i sporządzenia<br />
miejscowego planu zagospodarowania przestrzennego, nie narusza własności osób<br />
prywatnych (rys. 2.10).<br />
Miejsce lokalizacji położone jest poza terenami bezpośredniego zagrożenia<br />
powodziowego oraz poza granicami zalewu wodą (rys. 2.11). Inwestycję zlokalizowano<br />
na gruncie w pełni uzbrojonym w niezbędną infrastrukturę techniczną, nie sąsiaduje<br />
również bezpośrednio z terenami zabudowy jednorodzinnej.<br />
Rys. 2.9 Lokalizacja scenariuszy na tle obszarów chronionych<br />
(http://www.akademiainspire.pl/[3])<br />
27
Rys. 2.10 Najkorzystniejsza lokalizacja biogazowni na tle katastru<br />
(http://www.akademiainspire.pl/[3])<br />
Rys. 2.11 Najkorzystniejsza lokalizacja biogazowni na tle projektu Studium Zagrożenia<br />
Powodziowego (http://www.akademiainspire.pl/[3]).<br />
28
3. GIS jako system wspomagania decyzji<br />
Podjęcie decyzji ma swój początek w gromadzeniu informacji. Tworzenie baz danych<br />
przestrzennych to początek na drodze do rozwiązania problemów lokalizacyjnych. Baza<br />
w formie elektronicznej umożliwia sortowanie danych, a przez to szybki dostęp do tych które<br />
nas interesują. To właśnie selektywny dostęp do danych GIS pozwala na tworzenie map<br />
tematycznych. Wielu użytkowników może korzystać z tej samej bazy i wykorzystywać ją do<br />
tworzenia map na własne potrzeby. System GIS umożliwia również przeprowadzanie analiz<br />
przestrzennych których wyniki przedstawia się w postaci map (Hejmanowska, 2005). Taka<br />
analiza wymaga doboru odpowiednich kryteriów, np. teren pod wyciąg narciarski powinien<br />
uwzględniać m.in. następujące:<br />
• nachylenie terenu – przykładowo minimum 10%;<br />
• zaleganie pokrywy śnieżnej – przykładowo powyżej 70 dni/rok;<br />
• minimalna powierzchnia terenu – przykładowo nie mniej niż10 ha.<br />
Właściwy dobór kryteriów jest bardzo istotny w trakcie podejmowania decyzji.<br />
W zależności od rodzaju kryteriów przyjętych do analiz rozróżnia się dwie metody<br />
wykorzystywane do rozwiązania danego problemu: twarde, miękkie.<br />
Zastosowanie kryterium np. powierzchnia terenu minimum 10 ha, charakteryzuje<br />
metodę twardą. Barierą nazywa się wartość progową przyjmowaną w analizach. Korzystając<br />
z twardych kryteriów analizy uzyskuje się mapę przedstawiającą tereny spełniające oraz<br />
niespełniające postawionych warunków. Stosując kryteria miękkie otrzymuje się stopień<br />
przydatności określonych obszarów dla określonego w trakcie analizy celu. Kryterium to nie<br />
jest definiowane jako jedna liczba, ale jako przestrzenny rozkład parametru związanego na<br />
przykład z nachyleniem terenu (Hejmanowska, 2004).<br />
3.1. Deterministyczne metody podejmowania decyzji<br />
Uwzględniając rodzaj posiadanych informacji, problemy decyzyjne dzieli się na<br />
podejmowane w warunkach (Hejmanowska, 2004):<br />
• pewności – decyzja pociąga za sobą określone, znane konsekwencje;<br />
29
• ryzyka – decyzja pociąga za sobą więcej niż jedną konsekwencję, znany jest zbiór<br />
możliwych konsekwencji i prawdopodobieństwa ich wystąpienia;<br />
• niepewności – gdy nie znamy prawdopodobieństw wystąpienia konsekwencji danej<br />
decyzji.<br />
Deterministyczne metody podejmowania decyzji pojawiają się w warunkach<br />
pewności. W procesie decyzyjnym wybierany jest wariant optymalny. Następnie definiuje się<br />
zestaw kryteriów lub warianty zestawu kryteriów, które są analizowane w trakcie procesu<br />
podejmowania decyzji z wykorzystaniem GIS. W procesie oceny kryteriów można stosować<br />
(Hejmanowska, 2005):<br />
• analizę wieloparametryczną MCE (Multi-Criteria Evaluation),<br />
• analizę z uwzględnieniem kilku celów MOE (Multi-Objective Evaluation).<br />
3.1.1. Analiza wieloparametryczna MCE<br />
Analiza MCE polega na wyselekcjonowaniu obszarów które jednocześnie spełniają<br />
wszystkie postawione kryteria. W tej analizie kryteria mogą być jednakowo ważne, ale jest<br />
możliwe również ich wagowanie. W przypadku kryteriów twardych w wyniku analizy<br />
otrzymuje się obszary spełniające i niespełniające zadane warunki. Analiza z miękkimi<br />
kryteriami wyznacza stopień przydatności danego obszaru.<br />
Analiza kryteriów może odbywać się trzema metodami (Hejmanowska, 2005):<br />
• Nakładkowania,<br />
jest to metoda twarda, polega na pomnożeniu przez siebie lub dodaniu do siebie<br />
wszystkich wyników analiz spełniających pojedyncze kryteria. Wynikiem takiej analizy<br />
jest albo część wspólna (operacja logiczna AND) albo suma logiczna (OR) składowych<br />
wyników analiz. Wynikiem operacji AND jest obszar dla którego zostały spełnione<br />
wszystkie kryteria, natomiast dla operacji OR wystarczy, że zostało spełnione<br />
przynajmniej jedno kryterium.<br />
• Liniowego wagowania WLC (Weighted Linear Combination),<br />
metoda miękka, polega na pomnożeniu każdego kryterium przez wagę tego kryterium<br />
i zsumowaniu wyników tego mnożenia<br />
30
(3.1)<br />
Gdzie:<br />
S – przydatność,<br />
w – waga kryterium,<br />
x – wartość parametru,<br />
i – kryterium,<br />
n – liczba kryteriów.<br />
Uzyskany wynik jest następnie mnożony przez bariery zgodnie ze wzorem, jeśli<br />
są dodatkowo zdefiniowane kryteria w postaci barier.<br />
(3.2)<br />
Gdzie:<br />
c j – j-ta bariera<br />
W analizie typu WLC czynniki (miękkie kryteria) standaryzuje się do ciągłej skali<br />
przydatności mieszczącej się w wartościach od 0 (najmniej przydatne) do 255<br />
(najbardziej przydatne). Takie przeskalowanie pozwala na połączenie i porównanie<br />
kryteriów ze sobą. Miękkie kryteria definiowane są jako obszary przydatności o nieostrej<br />
granicy pomiędzy terenami przydatnymi i nieprzydatnymi. Natomiast bariery mają<br />
charakter kryteriów twardych, mapa przydatności jest zero – jedynkowa, obszar jest<br />
przydatny (1) lub nieprzydatny (0).<br />
Standaryzacja<br />
Do standaryzacji czynników zazwyczaj wykorzystuje się teorię zbiorów<br />
rozmytych (fuzzy sets), gdzie liniowa funkcja wzrostu/spadku przydatności jest jedną<br />
z możliwości. Najważniejszym problemem standaryzacji jest kwestia wyboru punktów<br />
końcowych, których przynależność do zbioru osiąga 0 lub 1 (0 lub 255).<br />
31
Istnieją cztery typy standaryzacji (Eastman, 2006):<br />
Funkcja sigmoidalna („s-kształtna”),<br />
jest najbardziej popularna. Charakteryzują ją cztery punkty przegięcia (a, b, c,<br />
d) wzdłuż osi x regulujące kształt krzywej (a – wzrost powyżej 0, b – osiągnięcie 1,<br />
c – spadek poniżej 1, d – osiągnięcie 0). Cztery warianty funkcji: monotonicznie<br />
rosnąca (a, b=c=d), monotonicznie malejąca (a=b=c, d), rosnąca – malejąca (a, b=c,<br />
d), rosnąca-stała-malejąca (a, b, c, d).<br />
Rys. 3.1 Rodzaje funkcji sigmoidalnej (Eastman, 2006).<br />
Funkcja „J-kształtna”(wykładnicza),<br />
również ma szerokie zastosowanie, choć w większości przypadków<br />
trafniejszym wydaje się zastosowanie funkcji sigmoidalnej. Wymaga podania<br />
położenia wartości czterech punktów: a – wzrost powyżej 0,5, b – osiągnięcie 1, c –<br />
spadek poniżej 1, d – spadek poniżej 0,5.<br />
Wykresy przedstawiają różne warianty funkcji: monotonicznie rosnąca (a, b=c=d),<br />
monotonicznie malejąca (a=b=c, d), rosnąca – malejąca (a, b=c, d), rosnąca – stała –<br />
malejąca (a, b, c, d).<br />
Rys. 3.2 Rodzaje funkcji „j-shaped” (Eastman, 2006).<br />
Funkcja liniowa,<br />
stosowana bardzo często w urządzeniach elektronicznych wykorzystujących<br />
logikę zbiorów rozmytych ze względu na prostotę i wykorzystywanie sygnałów<br />
czujnika – najczęściej liniowego. Warianty funkcji podobne jak w dwóch poprzednich<br />
typach funkcji.<br />
32
Rys. 3.3 Rodzaje funkcji liniowej (Eastman, 2006).<br />
Funkcja zdefiniowana przez użytkownika<br />
W przypadku gdy żadna ze standardowych funkcji nie może być zastosowana<br />
użytkownik podaje funkcję własną. Może zawierać dowolną ilość punktów<br />
kontrolnych (przegięć funkcji). Pomiędzy punktami wartość przynależności jest<br />
interpolowana liniowo.<br />
Rys. 3.4 Funkcja zdefiniowana przez użytkownika (Eastman, 2006).<br />
Wagi kryteriów<br />
Wagi kryteriów mogą być opracowywane na wiele sposobów. W prostych<br />
przypadkach można na podstawie swojej wiedzy i doświadczeń przypisać wartości<br />
wag (podział 1,0 na poszczególne kryteria). Jeśli pojawia się większa ilość kryteriów<br />
i należy rozważyć wszystkie ich wzajemne relacje stosuje się metodę porównania<br />
parami. Zastosowanie takiego sposobu daje możliwość otrzymywania pewniejszych<br />
ocen względnej ważności poszczególnych kryteriów. Technika porównywania parami<br />
jest częścią procedury decyzyjnej AHP której szerszy opis znajduje się w rozdziale<br />
3.3.<br />
• Uporządkowanej średniej ważonej OWA (Ordered Weighted Average),<br />
metoda miękka, stanowi rozwinięcie metody liniowego wagowania WLC. W metodzie tej<br />
dodatkowo wprowadza się dodatkowy zestaw wag porządkujących. Działanie wag<br />
porządkujących ilustruje rysunek 3.5.<br />
33
Rys. 3.5 Kompromis w zależności od ryzyka (Hejmanowska, 2004).<br />
Metoda OWA jest najogólniejszym przypadkiem analizy, natomiast metoda<br />
nakładkowania i metoda liniowego wagowania są szczególnymi przypadkami metody<br />
OWA.<br />
3.1.2. Ocena wielopodmiotowa MOE<br />
Ocena uwzględniająca kilka celów może dotyczyć wspólnych celów lub celów<br />
przeciwstawnych. Proces decyzyjny w pierwszym przypadku może być przeprowadzony przy<br />
wykorzystaniu jednej z metod MCE. W modelach decyzyjnych w przypadku sytuacji<br />
konfliktowej obszar może zostać przydzielony zgodnie z jednym z celów. Rozwiązanie<br />
sytuacji konfliktowej może przebiegać etapowo:<br />
• uszeregowanie celów zgodnie z priorytetem ważności;<br />
• wybranie celu o najwyższym priorytecie i wykonanie dla niego analizy przydatności.<br />
Wynik takiej analizy można potraktować jako kryterium typu bariera i kontynuować<br />
analizę biorąc pod uwagę następne cele (Hejmanowska, 2005).<br />
3.2. Niedeterministyczne metody podejmowania decyzji<br />
Najczęściej jednak proces decyzyjny przebiega w warunkach niepewności.<br />
Wykorzystywane dane mają różną dokładność, co od samego początku obciąża proces<br />
decyzyjny. Innym źródłem błędów jest niemożność idealnego wymodelowania relacji między<br />
danymi. Wybrana reguła decyzyjna również stanowi źródło błędów. Istnieje różna<br />
kombinacja uwzględniania błędności baz danych GIS w procach decyzyjnych (rys 3.6).<br />
Uwzględnić można tylko niedokładność danych, dalszą analizę przeprowadza się<br />
uwzględniając jedynie twarde kryteria. Można również przeprowadzić pełną analizę<br />
34
dokładności i przeprowadzić „miękką” regułę decyzyjną, dającą prawdopodobieństwo<br />
wystąpienia przydatności jakiegoś zjawiska (Hejmanowska, 2004).<br />
Rys. 3.6 Problem niedokładności w procesie podejmowania decyzji (Hejmanowska, 2004).<br />
Dokładność danych źródłowych<br />
Proces podejmowania decyzji należy rozpocząć od rozpoznania dokładności danych<br />
źródłowych. W przypadku danych pomiarowych dokładność może być określana przez błąd<br />
średni. Może on dotyczyć niedokładności związanej z lokalizacją przestrzenną lub<br />
z niedokładnością atrybutów. Dokładność niemierzalnych danych wyraża się w sposób<br />
proporcjonalny, np. jako proporcjonalny błąd nadmiaru lub/i niedomiaru wynikowej klasy<br />
uzyskanej w rezultacie automatycznej klasyfikacji obrazu (Hejmanowska, 2005)<br />
Błędy modelu<br />
Proces wspomagania decyzji systemami GIS wymaga zgromadzenia danych<br />
wejściowych (np. DTM, mapa pokrycia terenu) które są obciążone błędami. Dane źródłowe<br />
w trakcie wykonywanych analiz są poddawane transformacjom. Błędy danych wejściowych<br />
przenoszą się na wyniki wszystkich wykonywanych analiz, stając się błędem modelu. Każda<br />
analiza jest modelowaniem jakiejś innej wielkości np. na podstawie współrzędnych oblicza<br />
się powierzchnię obiektu, z DTM modelowane są nachylenia. System GIS wykorzystuje<br />
w tych przypadkach wybrany algorytm obliczeniowy, który mniej lub bardziej adekwatnie<br />
modeluje wartości powierzchni terenu czy nachylenia w danym punkcie. Nieadekwatność<br />
modelu w stosunku do rzeczywistości oraz błąd danych źródłowych sumują się i określają<br />
błąd modelu (Hejmanowska, 2004).<br />
35
Wiarygodność reguł decyzyjnych<br />
Analiza przestrzenna składa się z poszczególnych analiz cząstkowych. Końcowym<br />
etapem procesu decyzyjnego jest zastosowanie wybranej reguły decyzyjnej. Oprócz<br />
„twardych” metod decyzyjnych, w których nie uwzględnia się niedokładności reguły<br />
decyzyjnej, można wykorzystać metody na definiowanych poziomach wiarygodności<br />
(Hejmanowska, 2005)<br />
• prawdopodobieństwie (metoda „rozmyta” – fuzzy);<br />
• prawdopodobieństwie warunkowym (metoda Bayesa);<br />
• ufności i wiarygodności (metoda Dempstera-Shafera).<br />
„Miękkie” reguły decyzyjne wykorzystuje się gdy brana jest pod uwagę<br />
niedokładność danych lub relacji. Ostatecznie wybrana metoda decyzyjna ma wpływ na<br />
wiarygodność wyniku analizy.<br />
3.3. AHP<br />
Analityczny proces hierarchiczny (AHP) (rys. 3.7) jest jedną z metod służących do<br />
wspomagania procesów decyzyjnych. Metoda opracowana przez Saaty’ego, ułatwia<br />
dokonywanie optymalnych wyborów w przypadku wielokryterialnych problemów<br />
decyzyjnych poprzez ich redukcje do serii porównań parami. W efekcie pozwala to na<br />
liczbowe przedstawienie ważności analizowanych elementów. Dzięki niej możliwe jest<br />
porównanie czynników zarówno jakościowych, jak i ilościowych. Nadanie hierarchii<br />
wykonuje się z wielu różnych punktów widzenia (w różnych kategoriach) uwzględniając<br />
priorytety oceniającego, pozwala to na wybór najlepszego wariantu<br />
(https://www.gddkia.gov.pl/ [9]).<br />
36
Rys. 3.7 Schemat hierarchiczny w metodzie Saaty’ego (https://www.gddkia.gov.pl/ [9]).<br />
Model AHP skład się z następujących etapów (Bohatkieiwicz i in., 2008):<br />
• określenie wariantów przedsięwzięcia.<br />
• Określenie celu analizy<br />
W etapie tym określa się cel analizy, jaki chce się uzyskać oraz spełniających<br />
do kryteriów, według których będą oceniane warianty.<br />
• Ustalenie skali ocen.<br />
Zwyczajowo w metodzie AHP używa się poniższej skali ocen wariantów:<br />
9 – pierwszy element jest zdecydowanie korzystniejszy do drugiego,<br />
7 – pierwszy element jest dużo korzystniejszy od drugiego,<br />
5 – pierwszy element jest wyraźnie korzystniejszy od drugiego,<br />
3 – pierwszy element jest nieznacznie korzystniejszy od drugiego,<br />
1 – oba elementy są jednakowo korzystne,<br />
1/3 – pierwszy element jest nieznacznie mniej korzystny od drugiego,<br />
1/5 – pierwszy element jest wyraźnie mniej korzystny od drugiego,<br />
1/7 – pierwszy element jest dużo mniej korzystny od drugiego,<br />
1/9 – pierwszy element jest zdecydowanie mniej korzystny od drugiego.<br />
Liczbami 2, 4, 6, 8, 1/2, 1/4, 1/6, 1/8 oznacza się oceny pośrednie.<br />
37
• Określenie istotności kryteriów.<br />
Do późniejszych obliczeń koniecznym jest określenie, które z analizowanych<br />
kryteriów jest najistotniejsze, a które mniej istotne. Każdemu z kryteriów przypisuje<br />
się odpowiednią wartość – wagę.<br />
• Porównanie wariantów.<br />
W następnym kroku porównuje się ze sobą wszystkie warianty przedsięwzięcia<br />
(na zasadzie każdy z każdym) względem każdego z kryteriów. Przy analizie,<br />
porównując każdą parę elementów danego poziomu, rozważa się następujące aspekty:<br />
Który z porównywanych wariantów jest ważniejszy/lepszy<br />
Jak bardzo wariant jest ważny<br />
• Analiza hierarchii rozwiązań.<br />
W kolejnym etapie wykonując obliczenia, mnoży się wynik oceny każdego<br />
wariantu dla danego kryterium oraz wagę danego kryterium. Wyniki mnożenia dla<br />
każdego wariantu dodaje się do siebie otrzymując jego globalna ocenę w aspekcie<br />
analizowanych kryteriów. Wykonując tę czynność dla każdego z wariantów otrzymuje<br />
się wynik w postaci hierarchii rozwiązań (macierzy).<br />
• Sprawdzenie poprawności otrzymanych wyników.<br />
Stosując metodę AHP należy sprawdzić, czy otrzymane wyniki nie naruszają<br />
zasady stałości preferencji – prawidłowego porównania wariantów oraz kryteriów<br />
między sobą. W tym celu należy obliczyć współczynnik niespójności. Przy<br />
prawidłowo skonstruowanej macierzy wartość współczynnika niespójności nie<br />
przekracza 0,1.<br />
Metoda ta mimo wielu zalet nie jest pozbawiona wad, o których napisali w swej<br />
publikacji Bohatkiewicz i in. (2008):<br />
• ograniczoną liczbę porównywanych elementów na tym samym poziomie hierarchii;<br />
• założenie pełnej porównywalności kryteriów i wariantów występujących w modelu, co<br />
nie jest zawsze zgodne z rzeczywistością (np. elementy ilościowe i jakościowe);<br />
• utrudnione uwzględnienie zależności pomiędzy cząstkowymi kryteriami celu.<br />
38
4. Opracowanie projektu GIS<br />
Kompleksowe opracowanie projektu GIS wymaga wiedzy eksperckiej z zakresu<br />
wymagań prawnych co do właściwej lokalizacji danego obiektu. Nie jest łatwym<br />
określenie wszystkich kryteriów mających wpływ na lokalizację, gdyż jest ich dużo,<br />
dodatkowo opisane są w wielu aktach prawnych. Ponadto należy uwzględnić również<br />
kryteria odnajdujące uzasadnienie jedynie w regułach wynikających ze sztuki<br />
projektowania. W rozdziale 4.1 zostały wymienione akty prawne, których zapisy<br />
warunkowały powstanie niektórych kryteriów.<br />
Istotnym jest znajomość analizowanego terenu, jego lokalizacja np. względem<br />
większych skupisk ludzi (ważne w przypadku analiz pod kątem odległości i czasu<br />
dojazdu z miasta do obiektu narciarskiego). Dokładne poznanie terenu może przyczynić<br />
się do zauważenia istotności nowych kryteriów, jak choćby lokalizacja obiektu<br />
rekreacyjnego w pobliżu atrakcyjnych turystycznie miejsc (występowanie parków<br />
narodowych, rezerwatów). Analizy terenu pod tym kątem dokonano w rozdziale 4.2.<br />
Kolejny krok w tworzeniu projektu GIS to pozyskanie danych o ukształtowaniu<br />
terenu, jego dotychczasowym zagospodarowaniu, planowanym zagospodarowaniu oraz<br />
szeregu informacji o środowiskowych uwarunkowaniach (opis danych źródłowych<br />
zawiera rozdział 4.4). Przy wykorzystaniu ww. danych źródłowych powstała przestrzenna<br />
baza danych potrzebnych do właściwych analiz (etapy powstania bazy zawiera rozdział<br />
4.5).<br />
Wykorzystanie danych, wraz z odpowiednim oprogramowaniem (którego opis<br />
i sposób wykorzystania opisuje rozdział 4.3), wyposaża użytkownika w komplet narzędzi<br />
potrzebnych do powstania projektu.<br />
4.1. Przepisy prawne dotyczące lokalizacji wyciągów<br />
Aby podjąć decyzję o lokalizacji należy uwzględnić wiele różnorodnych<br />
kryteriów, które opisane są w poniższych ustawach i rozporządzeniach.<br />
39
Ustawy:<br />
• Ustawa z dnia 3 lutego 1995 r. o ochronie gruntów rolnych i leśnych (Dz.U.<br />
z 2004 Nr 121, poz. 1266 z późn. zm),<br />
• Ustawa z dnia 18 lipca 2001 r. Prawo wodne (Dz.U. 2001 Nr 115 poz. 1229),<br />
• Ustawa z dnia 31 stycznia 1959 r. o cmentarzach i chowaniu zmarłych (Dz.U.<br />
z 2000 Nr 23, poz. 295 z późn. zm),<br />
• Ustawa z dnia 16 kwietnia 2004 r. o ochronie przyrody (Dz.U. Nr 92, poz. 880<br />
z późn. zm),<br />
• Ustawa z dnia 23 lipca 2003 r. o ochronie zabytków i opiece nad zabytkami<br />
(Dz.U. Nr 162, poz. 1568 z późn. zm),<br />
• Ustawa z dnia 21 marca 1985 r. o drogach publicznych (Dz.U. z 2007 Nr 19, poz.<br />
115 z późn. zm),<br />
• Ustawa z dnia 27 kwietnia 2001 r. Prawo ochrony środowiska (Dz.U. z 2008 Nr<br />
25, poz. 150 z późn. zm).<br />
Rozporządzenia:<br />
• Rozporządzenie Ministra Infrastruktury z dnia 7 marca 2008 r. w sprawie<br />
wymagań, jakie muszą spełniać cmentarze, groby i inne miejsca pochówku zwłok<br />
i szczątków (Dz. U. z dnia 21 marca 2008 r.),<br />
• Rozporządzenie Ministra Infrastruktury z dnia 11 grudnia 2003 r. w sprawie<br />
zasadniczych wymagań dla kolei linowych przeznaczonych do przewozu osób<br />
(Dziennik Ustaw z dnia 3 lutego 2004 r. Nr 15 poz. 130),<br />
• Rozporządzenie Rady Ministrów z dnia 6 maja 1997 r. w sprawie określenia<br />
warunków bezpieczeństwa osób przebywających w górach, pływających,<br />
kąpiących się i uprawiających sporty wodne. (Dz. U. z 7 czerwca 1997 r.),<br />
• Rozporządzenie Ministra Gospodarki Komunalnej z dnia 25 sierpnia 1959 r.<br />
w sprawie określenia, jakie tereny pod względem sanitarnym są odpowiednie na<br />
cmentarze (Dz.U. 1959 Nr 52, poz. 315),<br />
• Rozporządzenie Ministra Transportu i Gospodarki Morskiej z dnia 2 marca 1999<br />
r. w sprawie warunków technicznych, jakim powinny odpowiadać drogi<br />
publiczne i ich usytuowanie (Dz.U. 1999 Nr 43, poz. 430).<br />
40
Nie wszystkie kryteria mają swoje prawne uwarunkowania, w analizie uwzględniono<br />
również takie które związane są z dobrą sztuką projektowania. Przykładem takiego<br />
kryterium jest ekspozycja, czy zaleganie pokrywy śnieżnej.<br />
4.2. Charakterystyka terenu badań<br />
Analizowany w niniejszej pracy obszar (rys. 4.1) obejmuje gminy: Miasto<br />
Kraków, Michałowice, Zielonki, Wielka Wieś, Kocmyrzów-Luborzyca, Słomniki –<br />
obszar wiejski, Iwanowice, Zabierzów, Jerzmanowice-Przeginia, Skała-miasto, Skała –<br />
obszar wiejski, Trzyciąż, Gołcza, Wolbrom – obszar wiejski, Olkusz – obszar wiejski.<br />
powierzchnia całego obszaru wynosi 66 157 ha. W obrębie terenu znajduje się północna<br />
część Krakowa (dzielnice: Stare Miasto, Grzegórzki, Prądnik Czerwony, Prądnik Biały,<br />
Mistrzejowice, Bieńczyce; przylegające do nich Krowodrza, Bronowice, Wzgórza<br />
Krzesławickie, Nowa Huta we fragmentach) (http://www.bip.krakow.pl/ [5]) oraz miasto<br />
Skała. W bliskim sąsiedztwie granic znajdują się miasta: Słomniki, Krzeszowice, Olkusz,<br />
Wolbrom, Miechów.<br />
Najludniejszym miastem jest Kraków (750 tys. osób), 37 tys. to liczba ludności<br />
w Olkuszu, w Krzeszowicach, Wolbromiu oraz Miechowie wynosi ok. 10 tys.<br />
W najmniejszych miastach Skała i Słomniki mieszka kolejno 3,5 tys., 4,3 tys. ludności<br />
(http://www.stat.gov.pl/ [14]).<br />
Teren usytuowany jest w pasie Wyżyny Olkuskiej (będącej częścią Wyżyny<br />
Krakowsko-Częstochowskiej) i Wyżyny Miechowskiej (http://www.widoczek.nets.pl/<br />
[16]). Wyżyny te to obszary typowo rolnicze, pokryte utworami lessowymi sprzyjającymi<br />
uprawom. Duży procent analizowanego terenu stanowią gleby o wysokiej przydatności<br />
rolniczej kompleksu 2 (pszennego dobrego), prawie 14% gleb to gleby kompleksu 1<br />
pszennego bardzo dobrego (patrz rozdział 5.2.5 rys. 5.14).<br />
Wyżyna Olkuska zbudowana jest z wapieni jurajskich. Ma bardzo urozmaiconą<br />
rzeźbę, oprócz ostańców krasowych występują tu wąwozy krasowe. Do interesujących<br />
form należą znajdujące się w dnach dolin iglice skalne m.in. Maczuga Herkulesa i Igła<br />
Deotymy w Dolinie Prądnika (Kondracki, 2000).<br />
W pobliżu Krakowa istnieją dwa kompleksy narciarskie: w Sieprawiu oddalonym<br />
o 22 km oraz w Podstolicach o 16 km. Oba położone są na południe od Krakowa.<br />
41
4.3. Oprogramowanie<br />
Rys. 4.1 Podział administracyjny analizowanego terenu.<br />
W niniejszej pracy wykorzystano oprogramowanie IDRISI Andes oraz GeoMedia<br />
Professional. Poniżej znajduje się ich krótka charakterystyka.<br />
IDRISI Andes<br />
IDRISI należy do programów z rodziny Systemów Informacji Geograficznej<br />
(GIS). Został opracowany na Uniwersytecie Clark w Stanach Zjednoczonych w ramach<br />
utworzonego tam Projektu IDRISI. Obecnie rozwój oraz dystrybucja oprogramowania<br />
odbywa się poprzez specjalną jednostkę Uniwersytetu Clark, zajmującą się zagadnieniami<br />
GIS – Clark Labs. Ponadto jednostka ta prowadzi prace badawcze z zakresu geografii<br />
i GIS, jej działalność na charakter niekomercyjny.<br />
IDRISI jest programem rastrowym, o wartości dydaktycznej, mimo to posiada<br />
zdolności analityczne na profesjonalnym poziomie. Dokonuje się w nim m.in. analiz<br />
danych przestrzennych w formacie rastrowym, interpolacji przestrzennych. Możliwa jest<br />
42
analiza cyfrowych modeli wysokości. Posiada elementy relacyjnej bazy danych. Cechuje<br />
się łatwym zarządzaniem danymi oraz łatwością wizualizacji wyników. Umożliwia<br />
konwersję i eksport do wielu plików. Ma zastosowanie również w teledetekcji<br />
(http://beskid.geo.uj.edu.pl/ [4]).<br />
GeoMedia Professional<br />
GeoMedia Professional jest oprogramowaniem GIS tworzonym przez firmę<br />
Intergraph. Przede wszystkim wykorzystywany jest jako narzędzie do analiz danych<br />
przestrzennych. Wyposażony jest w zestaw narzędzi do prowadzenia analiz takich jak:<br />
kwerendy dotyczące atrybutów i położenia obiektów, wyznaczanie stref buforowych,<br />
przestrzenne nakładanie warstw tematycznych. Dane mogą być przechowywane<br />
w plikach formatu GeoMedia, jak również w systemach zarządzania bazami danych jak<br />
Oracle, Microsoft Access czy SQL. Program obsługuje wszystkie ważniejsze formaty<br />
plików z danymi przestrzennymi (CAD oraz standardowe formaty relacyjnych baz<br />
danych). Program umożliwia prowadzenie analiz obrazów rastrowych, tworzenie modeli<br />
DEM (Longley, 2006).<br />
Przygotowanie danych oraz ich obróbkę wykonano w programie GeoMedia.<br />
Przygotowanie polegało na zaimportowaniu danych źródłowych (w formacie<br />
wektorowym) do programu, następnie wykonaniu analiz przestrzennych polegających na<br />
zapytaniach do baz danych (tzw. zapytania atrybutowe). W ten sposób wyodrębniono np.<br />
kategorie dróg (każda droga posiadała swój atrybut przypisujący ją do określonej<br />
kategorii). Pogram ten posłużył również do eksportu przygotowanych danych do formatu<br />
SHAPE (format reprezentacji danych wektorowych, opracowany przez firmę ESRI).<br />
W końcowym etapie zaprezentowano w nim wyniki analiz.<br />
W programie IDRISI wykonano wszystkie zasadnicze etapy pracy.<br />
Przeprowadzono szereg operacji matematycznych na obrazach rastrowych. Jedną<br />
z najczęściej wykonywanych operacji były analizy odległości (etapem pośrednim<br />
standaryzacji czynników często było wyznaczenie odległości od danego obiektu;<br />
wyznaczanie buforów; generowanie mapy kosztów). Korzystano również z analiz<br />
sąsiedztwa (filtracja, mapa spadków, ekspozycji, analiza kosztów). W efekcie utworzono<br />
mapy wszystkich kryteriów będących składowymi analizy MCE wykonanej również<br />
w tym programie.<br />
43
4.4. Dane<br />
Dane wykorzystane w pracy pochodzą z następujących źródeł:<br />
• Numeryczna mapa hydrograficzna Polski w skali 1:50 000 z Wojewódzkiego<br />
Zasobu Geodezyjnego i Kartograficznego w Krakowie<br />
• Mapa Geośrodowiskowa Polski w skali 1:50 000 pochodząca z Państwowego<br />
Instytutu Geologicznego w Warszawie<br />
• Mapa glebowo-rolnicza w skali 1:25000 z Instytutu Uprawy Nawożenia<br />
i Gleboznawstwa w Puławach.<br />
• Baza Danych Topograficznych - drogi oraz koleje w formacie SHAPE (SKJZ_L,<br />
SKKL_L),<br />
• VMap w formacie SHAPE<br />
Dane z powyżej wymienionych źródeł są danymi wektorowymi. Dodatkowo<br />
w pracy korzystano z mapy użytkowania terenu w 19 kategoriach oraz Numerycznego<br />
Model Terenu o rozdzielczości geometrycznej 30 m. Są to mapy rastrowe. Szczegółowy<br />
wykaz użytych danych i ich źródeł zawiera tabela 4.1. Dodatkowo z VMAP-y<br />
wykorzystano dane dotyczące granic administracyjnych. Posłużyły one do prezentacji<br />
wyników z poszczególnych etapów pracy.<br />
Wszystkie wymienione wyżej dane otrzymano z Katedry Geoinformacji,<br />
Fotogrametrii i Teledetekcji Środowiska Wydziału Geodezji Górniczej i Inżynierii<br />
Środowiska <strong>AGH</strong> w Krakowie.<br />
44
Tab. 4.1 Wykaz użytych danych i ich źródeł do opracowania kryteriów.<br />
L.p. Kryteria Wykorzystane dane<br />
1. Mapa czasu podróży<br />
2.<br />
Mapa zalegania pokrywy<br />
śnieżnej<br />
3. ekspozycja<br />
4. Mapa spadków<br />
5.<br />
Mapa przydatności<br />
rolniczej gleb<br />
6. Mapa odległości od dróg<br />
Mapa pokrycia<br />
Mapa dróg (SKJZ_L.shp)<br />
NMT_30<br />
Gleby prądnik, baza_gleby.mdb<br />
Mapa dróg (SKJZ_L.shp)<br />
Pochodzenie<br />
danych<br />
TBD<br />
Mapa glebowo<br />
rolnicza<br />
TBD<br />
7.<br />
8.<br />
9.<br />
10.<br />
11.<br />
Mapa terenów<br />
podmokłych<br />
Mapa form ochrony<br />
przyrody<br />
Mapa odległości od<br />
terenów chronionych<br />
Mapa odległości od<br />
pomników przyrody<br />
i obiektów dziedzictwa<br />
kulturowego<br />
Mapa odległości od<br />
zwartej zabudowy<br />
12. Mapa odległości od wód<br />
Tereny podmokłe - tereny_p<br />
Bagna, trzęsawiska -<br />
bagno_trzesawisko_p<br />
Parkia narodowe - f32c_park_n<br />
Otuliny parków narodowych -<br />
f32c_stref_pn<br />
Rezerwaty - f35c_rezerwat<br />
Obszar chronionego krajobrazu -<br />
f34c_ochk<br />
Parkia narodowe - f32c_park_n<br />
Otuliny parków narodowych -<br />
f32c_stref_pn<br />
Rezerwaty - f35c_rezerwat<br />
Pomniki przyrody ożywionej -<br />
f37l_przyroda, f37p_przyroda<br />
Pomniki przyrody nieożywionej -<br />
f38p_ob_prz_n<br />
Obiekty dziedzictwa<br />
kulturowego - f39c_dz_kult<br />
Zabudowa zwarta<br />
Rzeki - rzeka_strumien_a,<br />
rzeka_strumien_l;<br />
Jeziora i stawy –<br />
jezioro_staw_a<br />
45<br />
Hydrologiczna<br />
mapa Polski<br />
VMap<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
Mapa pokrycia<br />
VMap<br />
hydrologiczna<br />
mapa Polski
13.<br />
Mapa odległości od<br />
cmentarzy<br />
14. Mapa ograniczeń<br />
Zbiorniki retencyjne -<br />
f19c_zbiornik<br />
cmentarz_a, cmentarz_p<br />
Sieć dróg (SKJZ_L.shp)<br />
Sieć torów kolejowych<br />
(SKKL_L.shp)<br />
Rzeki - rzeka_strumien_a,<br />
rzeka_strumien_l;<br />
Jeziora i stawy - jezioro_staw_a<br />
Zbiorniki retencyjne -<br />
f19c_zbiornik<br />
Parkia narodowe - f32c_park_n<br />
Otuliny parków narodowych -<br />
f32c_stref_pn<br />
Rezerwaty - f35c_rezerwat<br />
Tereny zalewowe -<br />
f20c_powodz<br />
Pomniki przyrody ożywionej -<br />
f37l_przyroda, f37p_przyroda<br />
Pomniki przyrody nieożywionej -<br />
f38p_ob_prz_n<br />
Obiekty dziedzictwa<br />
kulturowego - f39c_dz_kult<br />
Strefy ochronne ujęć wody -<br />
f22p_ujecie<br />
Mapa<br />
geośrodowiskowa<br />
VMap<br />
TBD<br />
TBD<br />
VMap<br />
hydrologiczna<br />
mapa Polski<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
Mapa<br />
geośrodowiskowa<br />
4.5. Opracowanie przestrzennej bazy danych<br />
Przygotowanie danych polegało na wczytaniu istniejących danych z różnych<br />
źródeł do programu Geomedia. Następnie przetworzono je (tj. wykonano analizy<br />
przestrzenne mające na celu uzyskanie nowych informacji np. z mapy gleb<br />
wyodrębniono kompleksy) i wyeksportowano do formatu SHAPE.<br />
Pierwszym etapem pracy było utworzenie nowej geoprzestrzeni, w której<br />
wyświetlono dostępne dane.<br />
Create a new GeoWorkspace using → Blank GeoWorkspace<br />
46
Następnie zdefiniowano system współrzędnych w układzie „1992”<br />
View → GeoWorkspace Coordinate System<br />
W poszczególnych zakładkach wybrano następujące parametry:<br />
• Ogólne informacje (General):<br />
Projection - układ współrzędnych prostokątnych płaskich,<br />
• Przestrzeń odwzorowania (Projection Space):<br />
odwzorowanie Transverse Mercator o parametrach:<br />
południk centralny: 19°,<br />
przesunięcie współrzędnej poziomej: 500 000m,<br />
przesunięcie współrzędnej pionowej: - 5 300 000m,<br />
skala na południku centralnym: 0,9993,<br />
• Przestrzeń geograficzna (Geographic Space):<br />
Geodezyjny układ odniesienia: WGS84<br />
Elipsoida referencji: WGS84<br />
Kolejny krok to utworzenie geohurtowni, w której umieszczono dane potrzebne<br />
do dalszej pracy.<br />
Warehouse → New Warehouse (wykorzystano istniejący szablon normal.mdt).<br />
Program wyświetla dane zawarte w geohurtowni wykorzystując system<br />
współrzędnych zdefiniowany dla geoprzestrzeni. Zmiana systemu współrzędnych nie<br />
wpływa na dane w geohurtowni a jedynie na dane w oknie mapy.<br />
W następnym etapie załadowano dane do projektu<br />
Warehouse → New Connection.<br />
Za pomocą tego narzędzia dodano do geohurtowni kolejno:<br />
• mapę geośrodowiskową – typ połączenia Access,<br />
• mapę gleby – typ połączenia Access,<br />
• mapę hydrograficzną – typ połączenia MapInfo,<br />
• TBD – typ połączenia ArcView<br />
• VMap – typ połączenia ArcView<br />
47
Połączenie typu Access oznacza, że treść mapy jest przechowywana w postaci<br />
bazy danych o rozszerzeniu MDB. Połączenie MapInfo zastosowano w przypadku mapy<br />
postaci bazy danych o rozszerzeniu DBF, natomiast typu ArcView – gdy plik był<br />
formatu SHAPE.<br />
Aby wyświetlić obiekty wyżej wymienionych map należało dodać je do legendy<br />
Legend → Add Legend Entries.<br />
Przykładową, wielowarstwową kompozycję map przestawiono na rys 4.2.<br />
Rys. 4.2 Wielowarstwowa kompozycja map w programie GeoMedia.<br />
48
5. Wielokryterialna analiza MCE<br />
5.1. Podział kryteriów<br />
Kryteria w analizie MCE dzieli się na kryteria typu bariera i kryteria typu<br />
parametr (czynnik). Podziału na bariery i parametry dokonano na podstawie ważności<br />
kryteriów lokalizacyjnych odnalezionych w przepisach prawnych oraz wynikających<br />
z zasad dobrego projektowania. W tabeli 5.1 zostały wymienione czynniki wraz ze<br />
skróconą informacją dotyczącą standaryzacji. W tabeli 5.2 wymieniono bariery.<br />
czynnik<br />
Mapa czasu dojazdu<br />
(patrz 5.2.1)<br />
Mapa zalegania pokrywy śnieżnej<br />
(patrz 5.2.2)<br />
Mapa ekspozycji<br />
(patrz 5.2.3)<br />
Mapa spadków<br />
(patrz 5.2.4)<br />
Mapa przydatności rolniczej gleb<br />
(patrz 5.2.5)<br />
Tab. 5.1 Podział czynników<br />
standaryzacja<br />
Mapa przydatności (fuzzy)<br />
Funkcja liniowa monotonicznie malejąca, c = 1, d = 123<br />
Mapa przydatności (fuzzy)<br />
funkcja liniowa monotonicznie rosnącą, a = 75, b = 100<br />
Mapa przydatności (reclass)<br />
Azymuty N, NW – 255, NE – 100, E, SE, S, SW, W – 50<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna monotonicznie rosnącą, a = 6, b = 20<br />
Mapa przydatności (reclass)<br />
Kompleksy 1,2,1z, lasy – 10; kompleksy 3, 4, 5, 8, 2z –<br />
100, kompleksy 6, 7, 9, 14, 3z – 230; gleby rolniczo<br />
nieprzydatne (RN), nieużytki – 255.<br />
Mapa odległości od dróg: A, S, GP<br />
(patrz 5.2.6)<br />
Mapa odległości od dróg: G, L, Z<br />
(patrz 5.2.6)<br />
Mapa terenów podmokłych<br />
(patrz 5.2.7)<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna monotonicznie rosnąca, a=1000,<br />
b=2000<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna symetryczna,<br />
a=30, b=100, c=1000, d=1500<br />
Mapa przydatności (reclass)<br />
bagna, trzęsawiska – 10, tereny podmokłe – 50,<br />
reszta – 255<br />
49
Mapa form ochrony przyrody<br />
[park narodowy, otulina parku<br />
narodowego, rezerwaty, obszar<br />
chronionego krajobrazu]<br />
(patrz 5.2.8)<br />
Mapa odległości od form ochrony<br />
przyrody<br />
[park narodowy, strefa parku<br />
narodowego, rezerwaty]<br />
(patrz 5.2.9)<br />
Mapa odległości od<br />
pomników przyrody<br />
[ożywionej, nieożywionej, obiekty<br />
dziedzictwa kulturowego]<br />
(patrz 5.2.9)<br />
Mapa odległości od zabudowy<br />
zwartej<br />
(patrz 5.2.10)<br />
Mapa odległości od wód<br />
[rzeki, zbiorniki wodne (jeziora,<br />
stawy, zbiorniki retencyjne)]<br />
(patrz 5.2.11)<br />
Mapa odległości od cmentarzy<br />
(patrz 5.2.12)<br />
Mapa przydatności (reclass)<br />
park narodowy, jego otulina oraz rezerwaty – 0, obszar<br />
chronionego krajobrazu – 150, pozostałe tereny - 255<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna symetryczna, a=5, b=10, c=3000,<br />
d=25000<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna symetryczna,<br />
punkty kontrolne: a=100, b=150, c=1000, d=4000<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna monotonicznie rosnąca, a=3000,<br />
b=4000<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna monotonicznie rosnąca a=100, b=200<br />
Mapa przydatności (fuzzy)<br />
Funkcja sigmoidalna monotonicznie rosnąca, a=150,<br />
b=1000<br />
Tab. 5.2 Podział barier<br />
bariera<br />
standaryzacja<br />
Drogi drogi – 0, pozostałe tereny – 1<br />
Rzeki rzeki – 0, pozostały teren – 1<br />
Zbiorniki wodne zbiorniki wodne – 0, pozostałe tereny – 1<br />
Tereny zalewowe tereny zalewowe – 0, pozostałe tereny – 1<br />
Park narodowy park narodowy – 0, pozostałe tereny – 1<br />
Otulina parku narodowego otulina parku narodowego – 0, pozostałe tereny – 1<br />
Rezerwaty rezerwaty – 0, pozostałe tereny – 1<br />
Zabudowa zabudowa – 0, pozostałe tereny – 1<br />
Ujęcia wody (bufor 100 m) ujęcia wody – 0, pozostałe tereny – 1<br />
Pomniki przyrody (bufor 100 m) pomniki przyrody – 0, pozostałe tereny – 1<br />
50
5.2. Opracowanie map czynników<br />
5.2.1. Mapa czasu podróży<br />
Pierwszym z czynników mogących wpływać na opłacalność inwestycji jest jej<br />
dostępność względem aglomeracji. W związku z tym postanowiono zamodelować czas<br />
dotarcia do danego miejsca z Krakowa zakładając, że podróżujący korzysta z samochodu,<br />
a opuszczając pojazd porusza się pieszo poza drogami.<br />
Prędkość przemieszczania w dużym stopniu zależy od terenu, po którym odbywa<br />
się ruch. Poruszanie się po drogach jest szybsze i wygodniejsze niż po bezdrożach.<br />
Szczególnie ważnym aspektem podróżowania są koszty, bardziej ekonomiczna jest jazda<br />
po drogach głównych, niż po drogach lokalnych.<br />
W systemach GIS obliczanie mapy kosztów rozpoczyna się od określenia kosztów<br />
jednostkowych dla różnych elementów pokrycia terenu (mapa tarcia). Jako koszt bazowy<br />
przyjęto 60 km/h (tarcie = 1j), jest to średnia prędkość poruszania się po drogach<br />
kategorii A oraz S, ruch pieszy jest wielokrotnie niższy np. 4 km/h (tarcie = 15j). Ujemna<br />
wartość kosztów jednostkowych oznacza, że na danym terenie nie można się<br />
przemieszczać (np. rzeka), takie tereny określa się mianem barier i przypisuje wartość -1.<br />
Wartości tarcia poszczególnych elementów pokrycia terenu zawiera tabela nr 5.3.<br />
drogi<br />
pokrycie<br />
Tab. 5.3 Koszty jednostkowe elementów pokrycia terenu.<br />
Element pokrycia terenu<br />
Średnia<br />
prędkość<br />
[km/h]<br />
Wartość<br />
tarcia<br />
1 klasa A 60 1<br />
2 klasa S 60 1<br />
3 klasa GP 60 1<br />
4 klasa G 30 2<br />
5 klasa Z 30 2<br />
6 klasa L 30 2<br />
1 zabudowa zwarta miejska 5 12<br />
2 zabudowa luźna wielorodzinna 5 12<br />
3 zabudowa luźna jednorodzinna 5 12<br />
4 zabudowa letniskowa 3 20<br />
5 zabudowa przemysłowa 3 20<br />
6 zabudowa handlowa 3 20<br />
51
7<br />
zabudowa użyteczności<br />
publicznej<br />
3 20<br />
8 Tereny komunikacyjne 10 6<br />
9 kopalnie i zwałowiska - -1<br />
10 budowy - -1<br />
11 zieleń urządzona 3 20<br />
12 tereny rolnicze 2 30<br />
13 użytki zielone 3 20<br />
14 nieużytki 3 20<br />
15 lasy 4 15<br />
16 zadrzewienia i zakrzewienia 4 15<br />
17<br />
tereny otwarte ze sporadyczną<br />
roślinnością<br />
4 15<br />
18 podmokłości 1 60<br />
19 wody - -1<br />
Mapy dróg oraz pokrycia przeskalowano stosując narzędzie RECLASS. Powstały<br />
w ten sposób dwa obrazy przedstawiające rozmieszczenie wartości tarcia (zał. 1 a. –<br />
pokrycie trenu, zał. 1 b. – drogi). Dodatkowo utworzono mapę barier (zał. 1 c.), na którą<br />
składają się wody i kolej, wokół których zastosowano 30-metrowy bufor. Obrazy<br />
pokrycia terenu i dróg z wartościami tarcia oraz bariery połączono stosując narzędzie<br />
Image Calculate → COVER.<br />
Wynik wykorzystano jako powierzchnię tarcia przy tworzeniu mapy kosztów.<br />
GIS Analysis → Distance Operators → COST<br />
Z dostępnych algorytmów wybrano Cost grow (rys. 5.1).<br />
Rys. 5.1 Zastosowanie algorytmu COST.<br />
52
Algorytm ten wymaga zdefiniowania punktu startowego, od którego będzie<br />
liczony koszt podróży. Jest nim punkt zlokalizowany w centrum Krakowa. Wynikiem<br />
analizy jest mapa kosztów (rys. 5.2)<br />
Rys. 5.2 Mapa kosztów.<br />
Kolejny etap to transformacja mapy kosztów na mapę czasu dojazdu.<br />
Wykorzystano do tego wzór:<br />
1,8<br />
mapa _ dojazdu = mapa _ kosztów∗<br />
60<br />
(5.1)<br />
Licznik współczynnika 1,8/60 to rozdzielczość piksela (30 m) podzielona przez<br />
przyjętą w analizie średnią prędkość poruszania się po drogach (17 m/s), mianownik<br />
przeskalował wynik na minuty.<br />
Powstały obraz przedstawia rozmieszczenie obszarów do których dojazd<br />
z centrum Krakowa może zająć od kilku minut do dwóch godzin (rys. 5.3).<br />
53
Rys. 5.3 Czas dojazdu [min].<br />
Rys. 5.4 Czas dojazdu – przydatność.<br />
Mapę czasu dojazdu przeskalowano narzędziem FUZZY stosując funkcję liniową<br />
monotonicznie malejącą (Linear Monotonically decreasing) o punktach kontrolnych<br />
c = 1, d = 123 (rys. 5.4). Tak przygotowana mapa jest jednym z czynników analizy MCE.<br />
5.2.2. Mapa zalegania pokrywy śnieżnej<br />
Zaleganie pokrywy śnieżnej warunkuje uprawianie narciarstwa. Duże znaczenie<br />
dla utrzymania odpowiedniej pokrywy śnieżnej mają wysokość terenu n.p.m., stosunki<br />
termiczne, anemometryczne, radiacyjne i wilgotnościowe (Nowosad, 1992).<br />
Odpowiednio niska temperatura pozwala na dłuższe utrzymanie pokrywy śnieżnej<br />
również w przypadku sztucznego naśnieżania.<br />
Nowosad (1992) w swym opracowaniu podaje wzór na liczbę dni z pokrywą<br />
śnieżną. Skonstruowane wzory empiryczne są uzależnione od charakterystyk<br />
meteorologicznych, od szerokości i długości geograficznej oraz wysokości terenu.<br />
Wykorzystanie modelu na zaleganie pokrywy daje pełny obraz rozmieszczenia<br />
najzimniejszych terenów. Na potrzebę przeprowadzonych w dalszej części pracy analiz,<br />
powstałe dane zastępują dane z nieposiadanej mapy temperatur.<br />
Opis metody zastosowanej w opracowaniu:<br />
Rozważania przeprowadzone zostały na przykładzie średniej rocznej liczby dni<br />
z pokrywą śnieżną w 56 stacjach i posterunkach zlokalizowanych w części polskiej<br />
Karpat. Ponownej analizie została poddana średnia z 42 stacji i posterunków położonych<br />
54
w Beskidach i na Pogórzach. Zastosowana metoda najmniejszych kwadratów pozwoliła<br />
na skonstruowanie wielomianów 1-ego, 2-ego i 3-ego stopnia. Stosowano tu tradycyjną<br />
metodę, wykorzystując średnie ze wszystkich stacji.<br />
Średnie pogrupowano w zależności od wysokości stacji. Dla przedziałów 100 m<br />
obliczono średnią wysokość n. p. m. stacji i posterunków z tego przedziału oraz średnią<br />
ze średnich liczb dni z pokrywą śnieżną. Dopiero tak utworzone pary liczb posłużyły do<br />
skonstruowania poniższych wzorów:<br />
• dla Beskidów i Pogórzy:<br />
wklęsłe formy terenu (209 – 710 m n.p.m.):<br />
l. dni = 0.1031 * h + 48.7 (5.2)<br />
<br />
<br />
wypukłe formy terenu (237 – 1022 m n.p.m.):<br />
l. dni = 0.0846 * h + 56.6 (5.3)<br />
42 stacje łącznie (209 – 1022 m n.p.m.):<br />
l. dni = 0.0872 * h + 56.3 (5.4)<br />
• dla polskiej części Karpat:<br />
wklęsłe formy terenu (209 – 1520 m n.p.m.):<br />
l. dni = 0.0889 * h + 53.5 (5.5)<br />
wypukłe formy terenu (237 – 1991 m n.p.m.):<br />
l. dni = 0.0821 * h + 58.2 (5.6)<br />
<br />
56 stacji łącznie (209 – 1991 m n.p.m.):<br />
l. dni = 0.0847 * h + 56.5 (5.7)<br />
Przedmiotem rozważań niniejszej pracy dyplomowej są tereny usytuowane<br />
w pasie Wyżyny Krakowsko-Częstochowskiej i Wyżyny Miechowskiej. W dostępnej<br />
literaturze nie znaleziono modelu opisującego zaleganie pokrywy śnieżnej dla tego<br />
obszaru. Zastosowano więc model dla Beskidów i Pogórzy obliczony przez Nowosada<br />
(1992). Okres zalegania pokrywy śnieżnej głównie zależy od wysokości terenu. Dla<br />
analizowanego obszaru wahają się one w przedziale od 200 do 500 m n.p.m. (wzór<br />
obejmuje wysokości 200 do 1200 m n.p.m.). Kolejnym czynnikiem mającym wpływ na<br />
zaleganie jest klimat. Klimat rejonu wyżyn jest nieco różny od klimatu Pogórza. Cechą<br />
55
charakterystyczną warunków klimatycznych Wyżyny Krakowsko-Częstochowskiej<br />
i Wyżyny Miechowskiej jest oddziaływanie krakowskiej miejskiej wyspy ciepła. Jej<br />
wpływ przejawia się podwyższeniem średniej rocznej temperatury o 1,2 st. C oraz<br />
przedłużeniem termicznych pór roku w stosunku do Pogórza Karpackiego<br />
(http://www.wrotamalopolski.pl [17]). Wynikiem zastosowania wzoru dla Pogórzy mogą<br />
być pewne przekłamania w ilości dni zalegania pokrywy śnieżnej na terenach<br />
wyżynnych.<br />
Zastosowanie wzorów wymaga rozgraniczenia terenów na wypukłe i wklęsłe.<br />
W tym celu użyto filtru uśredniającego (Mean Filter) 7x7 (rys. 5.5).<br />
GIS Analysis → Context Operators → FILTER<br />
Filtracja sama w sobie polega na przekształceniu obrazu źródłowego w celu<br />
wydobycia z niego informacji które są potrzebne do dalszej obróbki. Filtr uśredniający<br />
jest filtrem wygładzającym, jego klasycznym zadaniem jest przestrzenne uśrednienie<br />
wartości obrazu prowadzące do redukcji szumów. W tym przypadku jest to uproszczona<br />
metoda wyselekcjonowania obszarów leżących powyżej lub poniżej lokalnej przeciętnej<br />
wysokości terenu.<br />
Rys. 5.5 NMT po filtracji uśredniającej (obszar<br />
gm. Michałowice).<br />
Rys. 5.6 Formy wklęsłe (odcienie zieleni)<br />
i wypukłe (odcienie żółci i czerwieni) (obszar<br />
gm. Michałowice).<br />
Dlatego kolejny etap to odjęcie nmt po filtracji od nmt przed filtracją. Wynikiem<br />
jest obraz przedstawiający formy wklęsłe i wypukłe (Rys. 5.6).<br />
56
Modułem RECLASS rozdzielono formy wklęsłe i wypukłe na dwa niezależne<br />
obrazy (o wartościach binarnych), pomnożone przez NMT przedstawiają rozkład<br />
wysokości (Rys. 5.7, Rys. 5.8).<br />
Rys. 5.7 Wysokości form wklęsłych (obszar<br />
gm. Michałowice).<br />
Rys. 5.8 Wysokości form wypukłych (obszar<br />
gm. Michałowice).<br />
Powstałe obrazy wysokości form wklęsłych i wypukłych podstawiono kolejno do<br />
wzorów na liczbę dni ze stałą pokrywą śnieżną dla form wklęsłych (wzór 5.2) oraz dla<br />
form wypukłych (wzór 5.3).<br />
Wynik zastosowania modułu dla obu form pomnożono przez obraz tych form<br />
o wartościach 0 i 1 (np. dla form wklęsłych: formy wklęsłe 1, formy wypukłe 0),<br />
zsumowane przedstawiają przestrzenne rozmieszczenie ilości dni ze stałą pokrywą<br />
śnieżną (rys. 5.10 w obszarze gminy Michałowice, rys. 5.9 w obszarze całego<br />
analizowanego terenu). Najmniejszą liczbą dni charakteryzują się południowe tereny<br />
form wklęsłych. Model obliczył dla Krakowa 70 dni zalegania.<br />
57
Rys. 5.9 Ilość dni ze stałą pokrywą śnieżną.<br />
Rys. 5.10 Ilość dni ze stałą pokrywą śnieżną<br />
(obszar gm. Michałowice).<br />
Rys. 5.11 Zaleganie pokrywy śnieżnej – mapa przydatności.<br />
Przy pomocy FUZZY przeskalowano powyższą mapę funkcją liniową<br />
monotonicznie rosnącą (Linear Monotonically increasing) o punktach kontrolnych<br />
a = 75, b = 100 (rys. 5.11).<br />
58
5.2.3. Ekspozycja<br />
Duży wpływ na warunki śnieżne i długość zalegania pokrywy śnieżnej ma<br />
ekspozycja stoku (kierunek jego nachylenia). Najlepsze są stoki północno-zachodnie<br />
i północne, ponieważ są najmniej nasłonecznione co przyczynia się do dłuższego<br />
zalegania na nich śniegu (http://www.tatrzanska.com/ [15]).<br />
Ekspozycję wyznaczono za pomocą funkcji ASPECT. Wynikiem działania<br />
modułu jest obraz przedstawiony w wartościach azymutu.<br />
GIS Analysis → Surface Analysis → Topographic Variables → ASPECT<br />
Wynik zreklasyfikowano (tab. 5.4) tworząc obraz przedstawiający teren<br />
z ekspozycją na osiem kierunków świata (rys. 5.12).<br />
L.p. Kierunek<br />
Tab. 5.4 Ekspozycje i współczynniki ich opłacalności.<br />
Azymut<br />
Współczynnik<br />
opłacalności (0-255)<br />
1. N 337.5 – 0 oraz 0 – 22.5 255<br />
2. NE 22.5 – 67.5 100<br />
3. E 67.5 – 112.5 50<br />
4. SE 112.5 – 157.5 50<br />
5. S 157.5 – 202.5 50<br />
6. SW 202.5 – 247.5 50<br />
7. W 247.5 – 292.5 50<br />
8. NW 292.5 – 337.5 255<br />
59
Rys. 5.12 Mapa ekspozycji stoków względem stron świata.<br />
Kolejny etap to przygotowanie mapy ekspozycji do analizy MCE w tym celu<br />
dokonano następnej reklasyfikacji przypisując wartości zamieszczone w tabeli 5.4, wynik<br />
reklasyfikacji znajduje się w załączniku 3 a.<br />
5.2.4. Mapa spadków<br />
Jednym z najistotniejszych kryteriów dotyczącym lokalizacji wyciągu jest<br />
nachylenie terenu. Mapa spadków została wygenerowana w stopniach na podstawie<br />
analizy powierzchni numerycznego modelu terenu za pomocą funkcji SLOPE<br />
GIS Analysis → Surface Analysis → Topographic Variables → SLOPE<br />
Pojedynczy piksel powstałego obrazu zawiera informację na temat wartości<br />
nachylenia w danym miejscu. Największe spadki znajdują się na terenie Ojcowskiego<br />
60
Parku Narodowego, jest to jednak obszar wykluczony z możliwości inwestowania.<br />
Zróżnicowanie pozostałego analizowanego terenu jest niewielkie. Tereny o nachyleniu<br />
powyżej 12 stopni, stanowią jedynie 3% powierzchni całego obszaru (powierzchnia<br />
terenów powyżej 12 st. wynosi 1237,5 ha, całkowita powierzchnia obszaru 41698,7 ha,<br />
obliczenia nie uwzględniają terenów będących barierą).<br />
Współczynnik opłacalności został przyporządkowany rosnąco w stosunku do<br />
rosnącej wartości spadków. Użyto modułu FUZZY, obraz przeskalowano funkcją<br />
sigmoidalną monotonicznie rosnącą o punktach przegięcia a = 6, b = 20 (zał. 3 b),<br />
oznacza to, że spadki powyżej 20 stopni mają najwyższy współczynnik przydatności 255.<br />
5.2.5. Mapa przydatności rolniczej gleb<br />
W ustawie o ochronie gruntów rolnych i leśnych można znaleźć kryterium<br />
odnośnie wykorzystania gruntów.<br />
Na cele nierolnicze i nieleśne można przeznaczać przede wszystkim grunty oznaczone<br />
w ewidencji gruntów jako nieużytki, a w razie ich braku – inne grunty o najniższej<br />
przydatności produkcyjnej. (art. 6 ust 1.)<br />
Przeznaczenia gruntów rolnych i leśnych na cele nierolnicze i nieleśne zawiera<br />
miejscowy plan zagospodarowania przestrzennego uwzględniający wymagane zgody<br />
określone w ust. 2 art. 7.<br />
Przeznaczenie na cele nierolnicze i nieleśne:<br />
1) gruntów rolnych stanowiących użytki rolne klas I–III, jeżeli ich zwarty obszar<br />
projektowany do takiego przeznaczenia przekracza 0,5 ha – wymaga uzyskania zgody<br />
Ministra Rolnictwa i Gospodarki Żywnościowej,<br />
2) gruntów leśnych stanowiących własność Skarbu Państwa – wymaga uzyskania zgody<br />
Ministra Ochrony Środowiska, Zasobów Naturalnych i Leśnictwa lub upoważnionej przez<br />
niego osoby,<br />
5) pozostałych gruntów leśnych<br />
– wymaga uzyskania zgody marszałka województwa wyrażanej po uzyskaniu opinii izby<br />
rolniczej.(art. 7 ust.2)<br />
61
Źródłem danych do pozyskania mapy przydatności rolniczej gleb jest Mapa<br />
glebowo-rolnicza z IUNG w Puławach w skali 1:25000. Poszczególne kompleksy gleb<br />
wyodrębniono w programie GeoMedia Professional wykorzystując zapytania trybutowe<br />
Analysis → Attribute Query.<br />
Obraz filtrowano narzędziem Filter gdzie jako atrybut wybrano „KOMPLEKS”<br />
do którego kolejno przypisywano wartości oznaczające kompleksy gleb.<br />
Rys. 5.13 Filtracja mapy - etap tworzenia mapy kompleksów gleb GeoMedia Professional<br />
Kolejnym krokiem było przejście z danymi z programu GeoMedia Professional do<br />
IDRISI. W tym celu wykonano eksport poszczególnych kompleksów w GeoMediach do<br />
formatu SHAPE.<br />
Warhouse →Export to → Export to Shapefile<br />
Następnie import plików SHAPE w programie IDRISI.<br />
File → Import → Software → Specific Formats → ESRI Formats → SHAPEIDR<br />
Powstałe obrazy zamieniono z formatu wektorowego na rastrowy.<br />
Reformat → RASTERVECTOR<br />
W przypadku kompleksów gleb jako opcję konwersji zaznaczono polygon to<br />
raster. Nazwę obrazu wyjściowego skopiowano z obrazu wejściowego. Modułem<br />
INITIAL skopiowano parametry przestrzenne z mapy pokrycia .<br />
Wymienione powyżej narzędzia programu IDRISI przedstawiono na zdjęciu obrazu<br />
z okna programu (zał. 2).<br />
62
Każdy kolejny kompleks reklasyfikowano przypisując im wartości od 1 do 16.<br />
Numery kompleksów oraz ich pełne nazwy wraz z współczynnikami opłacalności<br />
zestawiono w tabeli 5.5.<br />
L.p.<br />
Tab. 5.5 Kompleksy przydatności rolniczej<br />
Kompleksy przydatności rolniczej na<br />
gruntach ornych<br />
Współczynnik<br />
opłacalności (0-<br />
255)<br />
1. Kompleks 1 Pszenny bardzo dobry 10<br />
2. Kompleks 2 Pszenny dobry 10<br />
3. Kompleks 3 Pszenny wadliwy 100<br />
4. Kompleks 4 Żytni bardzo dobry 100<br />
5. Kompleks 5 Żytni dobry 100<br />
6. Kompleks 6 Żytni słaby 230<br />
7. Kompleks 7 Żytni bardzo słaby 230<br />
8. Kompleks 8<br />
Zbożowo-pastewny<br />
mocny<br />
100<br />
9. Kompleks 9 Zbożowo-pastewny słaby 230<br />
10. Kompleks 14<br />
Gleby orne przeznaczone<br />
pod użytki zielone<br />
230<br />
Kompleksy przydatności rolniczej na<br />
użytkach zielonych<br />
11. Kompleks 1z Bardzo dobry i dobry 10<br />
12. Kompleks 2z Średni 100<br />
13. Kompleks 3z Słaby i bardzo słaby 230<br />
14. Gleby rolniczo nieprzydatne (RN) 255<br />
15. Nieużytki 255<br />
16. Lasy 10<br />
Zsumowano wszystkie klasy w jedną mapę (rys. 5.14) wykorzystując kalkulator<br />
obrazów (Image Calculator). Tak przygotowaną mapę przeskalowano w sposób, aby<br />
tereny na których najbardziej opłaca się inwestować (gleby o najniższych klasach)<br />
otrzymały jak największy współczynnik opłacalności (wartości współczynników<br />
opłacalności zawiera tabela 5.5, obraz po przeskalowaniu zawiera załącznik 3 j).<br />
63
Rys. 5.14 Mapa przydatności rolniczej gleb.<br />
5.2.6. Mapa odległości od dróg<br />
Artykuł 43 ustawy z dnia 21 marca 1985 r. o drogach publicznych określa<br />
minimalne odległości obiektów budowlanych od zewnętrznej krawędzi jezdni (tab. 5.6).<br />
Tab. 5.6 Minimalne odległości obiektów budowlanych od zewnętrznej krawędzi jezdni (Ustawa<br />
o drogach publicznych).<br />
W rozporządzeniu w sprawie warunków technicznych, jakim powinny<br />
odpowiadać drogi publiczne i ich usytuowanie znajduje się paragraf 4, który<br />
przyporządkowuje klasy dróg do odpowiednich kategorii<br />
64
„§ 4. 1. W celu określenia wymagań technicznych i użytkowych wprowadza się<br />
następujące klasy dróg:<br />
1) autostrady, oznaczone dalej symbolem "A",<br />
2) ekspresowe, oznaczone dalej symbolem "S",<br />
3) główne ruchu przyspieszonego, oznaczone dalej symbolem "GP",<br />
4) główne, oznaczone dalej symbolem "G",<br />
5) zbiorcze, oznaczone dalej symbolem "Z",<br />
6) lokalne, oznaczone dalej symbolem "L",<br />
7) dojazdowe, oznaczone dalej symbolem "D".<br />
2. Drogi zaliczone do jednej z kategorii, w rozumieniu przepisów o drogach<br />
publicznych, powinny mieć parametry techniczne i uŜytkowe odpowiadające<br />
następującym klasom dróg:<br />
- drogi krajowe - klasy A, S, GP i wyjątkowo klasy G,<br />
- drogi wojewódzkie - klasy G, Z i wyjątkowo klasy GP,<br />
- drogi powiatowe - klasy G, Z i wyjątkowo klasy L,<br />
- drogi gminne - klasy L, D i wyjątkowo klasy Z.”<br />
Klasy dróg podobnie jak kompleksy gleb wyodrębniono w programie GeoMedia<br />
stosując zapytania atrybutowe. Przejście z danymi z GeoMediów do IDRISI przebiegało<br />
analogicznie jak w przypadku mapy przydatności rolniczej gleb. Mapę sieci dróg<br />
przedstawiono w 6 kategoriach A – autostrady, S – drogi ekspresowe, GP – drogi główne<br />
ruchu przyspieszonego, G – drogi główne, Z – drogi zbiorcze, L – drogi lokalne (rys.<br />
5.15)<br />
Rys. 5.15 Sieć dróg, centrum Krakowa i obszar na północny zachód od Krakowa.<br />
65
Klasy dróg podzielono na dwie kategorie. Pierwsza z nich, czyli drogi A, S, GP są<br />
to drogi krajowe łączące duże miasta, umożliwiają szybką komunikację między nimi.<br />
Jako najważniejszy powód odsunięcia inwestycji od tych dróg uznano hałas. Optymalna<br />
odległość od tych dróg stanowi czynnik zwiększający atrakcyjność terenu. Optymalna -<br />
rozumiana jako jednocześnie ułatwiająca dotarcie turystom do tras zjazdowych i nie<br />
wpływająca niekorzystnie na walory turystyczne.<br />
Druga kategoria do której zaliczono drogi klas G, L oraz Z to drogi o mniejszym<br />
natężeniu ruchu i hałasu. W projekcie inwestycji przewidziane jako drogi bezpośredniego<br />
dojazdu pod wyciąg narciarski.<br />
Oba czynniki poddano standaryzacji. Pierwszym krokiem do standaryzacji jest obliczenie<br />
odległości od obiektów, wykorzystano w tym celu funkcję DISTANCE (Przykładowy<br />
obraz odległości od dróg A, S, GP przedstawia rysunek 5.16).<br />
GIS Analysis → Distance Operators → Distance<br />
Proces standaryzacji przebiegł następująco:<br />
Czynnik drogi A, S, GP standaryzowano przy użyciu funkcji sigmoidalnej<br />
monotonicznie rosnącej (Sigmoidal Monotonically increasing) o punktach przegięcia<br />
a=1000, b=2000.<br />
GIS Analysis → Desicion Support → FUZZY.<br />
Efektem standaryzacji jest rozkład przydatności o wartościach w zakresie od 0 do<br />
255. W odległości do 1 km teren jest nieprzydatny, jego przydatność wzrasta wraz ze<br />
wzrostem odległości do 2 km, dalej teren jest tak samo przydatny (rys 5.17).<br />
Czynnik drogi G, Z, L standaryzowano funkcją sigmoidalną (Sigmoidal<br />
Symmetric) o punktach przegięcia a=30, b=100, c=1000, d=1500. Efektem takiego<br />
doboru parametrów standaryzacji jest obraz przedstawiony na rysunku 5.18. Teren jest<br />
nieprzydatny w odległości 30 m od drogi, jego przydatność wzrasta do odległości 100 m<br />
i utrzymuje się do 1 km, po przekroczeniu 1 km maleje do osiągnięcia odległości 1.5 km.<br />
66
Rys. 5.16 Odległość od dróg A,<br />
S, GP.<br />
Rys. 5.17 Odległość od dróg<br />
A, S, GP – przydatność.<br />
Rys. 5.18 Odległość od dróg G,<br />
Z, L – przydatność.<br />
5.2.7. Mapa terenów podmokłych<br />
Inwestowanie na terenach podmokłych generuje dodatkowe koszty osuszania<br />
dlatego przypisano im małe wartości przydatności.<br />
Mapa terenów podmokłych powstała dwóch elementów:<br />
• tereny podmokłe - pochodzące z numerycznej mapy hydrograficznej Polski,<br />
• bagna, trzęsawiska – obiekt VMapy.<br />
Kolejny etap to reklasyfikacja, poszczególnym elementom przypisano nowe<br />
wartości, charakteryzujące opłacalność inwestowania na danym terenie: bagna,<br />
trzęsawiska-10, tereny podmokłe-50, reszta-255 (zał. 3 c).<br />
5.2.8. Mapa form ochrony przyrody<br />
Elementami mapy form ochrony przyrody są (rys. 5.19):<br />
• park narodowy,<br />
• otulina parku narodowego,<br />
• rezerwaty,<br />
• obszary chronionego krajobrazu.<br />
Mapa form ochrony przyrody została zreklasyfikowana, każdemu z elementów<br />
przypisano współczynnik opłacalności: Ojcowski Park Narodowy i jego otulina, a także<br />
rezerwaty – 0, obszary chronionego krajobrazu – 150, pozostałe tereny – 255 (zał. 3 d).<br />
67
W obszarach parku narodowego, rezerwatów, otuliny parku narodowego zabrania<br />
się budowania nowych obiektów budowlanych, dlatego te tereny otrzymały wartość 0.<br />
Obszary chronionego krajobrazu są terenami wartościowymi ze względu na możliwość<br />
zaspokajania potrzeb związanych z turystyką i wypoczynkiem, podlegają ochronie na<br />
podstawie ustawy z 16 kwietnia 2004 r. o ochronie przyrody. Na takim obszarze może<br />
być wprowadzony zakaz realizacji przedsięwzięć mogących znacząco oddziaływać na<br />
środowisko w rozumieniu przepisów ustawy z 3 października 2008 r. o udostępnianiu<br />
informacji o środowisku i jego ochronie, udziale społeczeństwa w ochronie środowiska<br />
oraz o ocenach oddziaływania na środowisko. Zakaz nie dotyczy realizacji przedsięwzięć,<br />
dla których przeprowadzona ocena oddziaływania na środowisko wykazała brak znacząco<br />
negatywnego wpływu na ochronę przyrody obszaru chronionego krajobrazu. Ta forma<br />
ochrony obejmuje znaczną część analizowanego terenu.<br />
Rys. 5.19 Mapa form ochrony przyrody.<br />
68
5.2.9. Mapy odległości od formy ochrony przyrody i pomników przyrody.<br />
Sąsiedztwo form ochrony przyrody i pomników przyrody stanowi o atrakcyjności<br />
terenów mogących stać się terenami tras narciarskich. Ładne widoki, cisza to wielkie<br />
zalety, dlatego w projekcie uwzględniono bliskie sąsiedztwo parków i rezerwatów<br />
przypisując takim terenom wysokie wartości przydatności.<br />
Mapa odległości od terenów chronionych.<br />
Na mapę form ochrony przyrody w tym przypadku składają się elementy:<br />
• park narodowy,<br />
• strefa ochronna parku narodowego,<br />
• rezerwaty.<br />
Od sumy tych elementów policzono odległość. Standaryzację przeprowadzono<br />
wykorzystując funkcję sigmoidalną (Sigmoidal Symmetric) o punktach przegięcia a=5,<br />
b=10, c=3000, d=20000 (zał. 3 e).<br />
Mapa odległości od pomników przyrody i obiektów dziedzictwa kulturowego<br />
Pomniki przyrody jako formę ochrony środowiska definiuje ustawa o ochronie<br />
środowiska:<br />
Pomnikami przyrody są pojedyncze twory przyrody żywej i nieożywionej lub ich skupiska<br />
o szczególnej wartości przyrodniczej, naukowej, kulturowej, historycznej lub<br />
krajobrazowej oraz odznaczające się indywidualnymi cechami, wyróżniającymi je wśród<br />
innych tworów, okazałych rozmiarów drzewa, krzewy gatunków rodzimych lub obcych,<br />
źródła, wodospady, wywierzyska, skałki, jary, głazy narzutowe oraz jaskinie (art. 40 ust.<br />
1).<br />
Na terenach niezabudowanych, jeżeli nie stanowi to zagrożenia dla ludzi lub mienia,<br />
drzewa stanowiące pomniki przyrody podlegają ochronie aż do ich samoistnego,<br />
całkowitego rozpadu (art. 40 ust. 2).<br />
O obiektach dziedzictwa kulturowego mówi ustawa o ochronie zabytków i opiece<br />
nad zabytkami. Ustawa definiuje zabytek jako nieruchomość lub rzecz ruchomą, ich<br />
części lub zespoły, będące dziełem człowieka lub związane z jego działalnością<br />
i stanowiące świadectwo minionej epoki bądź zdarzenia, których zachowanie leży<br />
69
w interesie społecznym ze względu na posiadaną wartość historyczną, artystyczną lub<br />
naukową(art. 3 ust. 1)<br />
Ochronie i opiece podlegają, bez względu na stan zachowania:<br />
1) zabytki nieruchome będące, w szczególności:<br />
g) parkami, ogrodami i innymi formami zaprojektowanej zieleni,<br />
h) miejscami upamiętniającymi wydarzenia historyczne bądź działalność wybitnych<br />
osobistości lub instytucji (art. 6 ust. 1)<br />
Stanowiska archeologiczne zawierające się w obiektach dziedzictwa kulturowego mapy<br />
można zaliczyć do zabytków archeologicznych na podstawie art.6 ust.1 pkt.3.<br />
Dla powyżej zdefiniowanych obiektów ustanowiono barierę 100 m, która zapewni<br />
ich nienaruszalność jak również stref bezpośrednio do nich przylegających.<br />
Elementy mapy:<br />
• pomniki przyrody ożywionej,<br />
• pomniki przyrody nieożywionej,<br />
• obiekty dziedzictwa kulturowego.<br />
Źródłem danych jest mapa geośrodowiskowa.<br />
Zastosowano standaryzację typu sigmoidalna (Sigmoidal Symmetric) o punktach<br />
przegięcia a=100, b=150, c=1000, d=4000 (zał. 3 f). Do 100 m to tereny nieprzydatne, dla<br />
terenów od 100 do 150 m przydatność rośnie, a powyżej 150 m aż do 1000 m<br />
przydatność jest najwyższa, po przekroczeniu 1000 m spada aż do osiągnięcia 4000 m,<br />
tereny powyżej 4 km to tereny nieprzydatne.<br />
5.2.10. Mapa odległości od zabudowy zwartej<br />
Kompleks narciarski aby zaspokoić wymogi turystyki rekreacyjnej powinien<br />
znajdować się w odpowiedniej odległości od dużych skupisk ludzi (rys. 5.20). Możliwość<br />
aktywnego wypoczynku z dala od hałasu i wielkomiejskiej zabudowy zapewni terenom<br />
większe zainteresowanie wśród narciarzy. Jako minimalną odpowiednią odległość<br />
zaproponowano 3 km.<br />
70
Rys. 5.20 Rozmieszczenie zabudowy zwartej i luźnej.<br />
Zabudowa zwarta została wyodrębniona z mapy pokrycia. Czynnik ten<br />
zestandaryzowano przy użyciu funkcji sigmoidalnej monotonicznie rosnącej (Sigmoidal<br />
Monotonically increasing) o punktach przegięcia a=3000, b=4000 (zał. 3 f). Przydatność<br />
terenów rośnie od 3 km do 4 km, powyżej 4 km tereny mają taką samą najwyższą<br />
przydatność.<br />
5.2.11. Mapa odległości od wód<br />
Wody analizowanego obszaru to w większości rzeki. Tereny wzdłuż rzek to<br />
zawsze potencjalne tereny zalewowe, dlatego zdecydowano się na odsunięcie inwestycji<br />
od wód. Oddalenie stało się istotnym również ze względu na częste występowanie wzdłuż<br />
rzek stoków nieatrakcyjnych pod kątem planowanej inwestycji tj. krótkich, bardzo<br />
stromych oraz nachylonych w kierunku rzeki. Z powyższych przyczyn uznano teren do<br />
100 m od rzek za nieprzydatne.<br />
71
Na mapę wód składają się:<br />
• rzeki;<br />
• zbiorniki wodne : jeziora i stawy; zbiorniki retencyjne.<br />
Standaryzację wykonano przy użyciu funkcji sigmoidalnej monotonicznie<br />
rosnącej (Sigmoidal Monotonically increasing) o punktach przegięcia a=100, b=200 (zał.<br />
3 h). Przydatność terenów rośnie od 100m do 200 m, powyżej 200 m tereny mają taką<br />
samą najwyższą przydatność.<br />
Poniżej na rysunku 5.21 przedstawiono południowo-wschodnią część<br />
analizowanej powierzchni z fragmentem terenu zalewowego na rzece Dłubni oraz<br />
zbiornikiem Zesławice; obrazuje również zbiorniki retencyjne będące w fazie projektu na<br />
rzece Baranówka.<br />
Rys. 5.21 Rozmieszczenie elementów mapy wód w południowo-wschodniej części<br />
analizowanego obszaru.<br />
5.2.12. Mapa odległości od cmentarzy<br />
Pewne obostrzenia, co do lokalizacji w sąsiedztwie cmentarzy znajdujemy<br />
w prawie.<br />
72
Na podstawie rozporządzenia ministra gospodarki komunalnej z dn. 25 sierpnia<br />
1959 r. w sprawie określenia, jakie tereny pod względem sanitarnym są odpowiednie na<br />
cmentarze wyznacza się strefy sanitarne od cmentarzy.<br />
(§ 3 ust.1): „Odległość cmentarza od zabudowań mieszkalnych, od zakładów<br />
produkujących artykuły żywności, zakładów żywienia zbiorowego bądź zakładów<br />
przechowujących artykuły żywności oraz studzien, źródeł i strumieni, służących do<br />
czerpania wody do picia i potrzeb gospodarczych, powinna wynosić co najmniej 150 m;<br />
odległość ta może być zmniejszona do 50 m pod warunkiem, że teren w granicach od 50<br />
do 150 m odległości od cmentarza posiada sieć wodociągową i wszystkie budynki<br />
korzystające z wody są do tej sieci podłączone.”<br />
Mapa cmentarzy zawiera dane typu powierzchnia oraz typu punkt (cmentarz_a –<br />
nazwa obiektu uzupełniona o literkę „a” oznacza , że jego powierzchnia jest odpowiednio<br />
duża by przedstawić go w skali mapy, cmentarz_p – literka „p” informuje o tym, że jest<br />
to mały obiekt którego lokalizację jedynie zaznaczono w formie punktu na mapie – nie<br />
w skali). Dane o cmentarzach pochodzą z VMapy.<br />
Do standaryzacji zastosowano funkcję sigmoidalną (Sigmoidal Monotonically<br />
increasing) o punktach przegięcia a=150, b=1000 (zał. 3, i). Zastosowanie tej funkcji dało<br />
obraz przydatności terenu, w którym obszar w odległości do 150 m od cmentarzy jest<br />
nieprzydatny, przydatność wzrasta do osiągnięcia 1km, dalej osiąga stałą, najwyższą<br />
wartość.<br />
5.3. Opracowanie mapy ograniczeń<br />
Ograniczenia lokalizacyjne dotyczą terenów, na których z przyczyn fizycznych<br />
bądź zakazów prawnych nie można budować nowych obiektów. Najoczywistszymi dla<br />
niniejszej inwestycji są wody, tereny zabudowane, drogi, tory kolejowe. Dodatkowo<br />
prawo chroni takie tereny jak parki narodowe i ich otuliny, rezerwaty, strefy ujęć wody,<br />
pomniki przyrody czy tereny zalewowe.<br />
Kolej<br />
Artykuł 53 ust. 2 ustawy o transporcie kolejowym określa odległość usytuowania<br />
obiektów budowlanych w odległości nie mniejszej niż 20 m od osi skrajnego toru.<br />
Dodatkowo: Odległości, o których mowa w ust. 2, dla budynków mieszkalnych, szpitali,<br />
73
domów opieki społecznej, obiektów rekreacyjno-sportowych, budynków związanych<br />
z wielogodzinnym pobytem dzieci i młodzieży powinny być zwiększone, w zależności od<br />
przeznaczenia budynku, w celu zachowania norm dopuszczalnego hałasu w środowisku,<br />
określonych w odrębnych przepisach(art. 53 ust. 3).<br />
Park narodowy i rezerwat<br />
Park narodowy oraz rezerwat są formami ochrony przyrody określonymi<br />
w ustawie o ochronie przyrody (art. 6 ust. 1).<br />
Park narodowy obejmuje obszar wyróżniający się szczególnymi wartościami<br />
przyrodniczymi, naukowymi, społecznymi, kulturowymi i edukacyjnymi, o powierzchni nie<br />
mniejszej niż 1 000 ha, na którym ochronie podlega cała przyroda oraz walory<br />
krajobrazowe(art. 8 ust. 1).<br />
Rezerwat przyrody obejmuje obszary zachowane w stanie naturalnym lub mało<br />
zmienionym, ekosystemy, ostoje i siedliska przyrodnicze, a także siedliska roślin, siedliska<br />
zwierząt i siedliska grzybów oraz twory i składniki przyrody nieożywionej, wyróżniające<br />
się szczególnymi wartościami przyrodniczymi, naukowymi, kulturowymi lub walorami<br />
krajobrazowym (art. 13 ust. 1).<br />
Na obszarach graniczących z parkiem narodowym wyznacza się otulinę parku<br />
Narodowego (art. 11 ust. 1). Ustawa definiuje otulinę jako:<br />
strefę ochronną graniczącą z formą ochrony przyrody i wyznaczoną indywidualnie dla<br />
formy ochrony przyrody w celu zabezpieczenia przed zagrożeniami zewnętrznymi<br />
wynikającymi z działalności człowieka (art. 5).<br />
W parkach narodowych oraz w rezerwatach przyrody zabrania się:<br />
1) budowy lub rozbudowy obiektów budowlanych i urządzeń technicznych, z wyjątkiem<br />
obiektów i urządzeń służących celom parku narodowego albo rezerwatu przyrody;<br />
9) niszczenia gleby lub zmiany przeznaczenia i użytkowania gruntów;<br />
22) wykonywania prac ziemnych trwale zniekształcających rzeźbę terenu (art. 15 ust. 1).<br />
Ujęcia wody<br />
Ograniczenia lokalizacyjne można odnaleźć w Prawie wodnym.<br />
W celu zapewnienia odpowiedniej jakości wody w ujęciu wody dla ludności wymagane<br />
jest ustanowienie stref ochronnych ujęć wody (podziemnej i powierzchniowej) oraz<br />
obszarów ochronnych zbiorników śródlądowych.<br />
74
Strefę ochronną ujęć stanowi obszar, na którym obowiązują zakazy i ograniczenia<br />
użytkowania gruntów i korzystania z wody. Strefa ochronna została podzielona na teren<br />
ochrony bezpośredniej i pośredniej.<br />
Na terenach ochrony pośredniej może być zabronione lub ograniczone wykonywanie<br />
robót oraz innych czynności powodujących zmniejszenie przydatności ujmowanej wody<br />
lub wydajności ujęcia, a w szczególności:<br />
wykonywanie robót melioracyjnych oraz wykopów ziemnych (art. 54 ust. 1 pkt. 6).<br />
(…) oprócz zakazów lub ograniczeń, o których mowa w ust. 1, może być zabronione lub<br />
ograniczone:<br />
lokalizowanie budownictwa mieszkalnego oraz turystycznego (art. 54 ust.3 pkt. 1).<br />
Strefa ochronna ujęcia wody z potoku górskiego lub z górnego biegu rzeki może<br />
obejmować całą zlewnię cieku powyżej ujęcia wody (art. 56 ust. 3).<br />
Tereny zalewowe<br />
Na obszarach szczególnego zagrożenia powodzią zabrania się wykonywania robót<br />
oraz czynności utrudniających ochronę przed powodzią lub zwiększających<br />
zagrożenie powodziowe, w tym:<br />
1) wykonywania urządzeń wodnych oraz budowy innych obiektów budowlanych;<br />
2) sadzenia drzew lub krzewów, z wyjątkiem plantacji wiklinowych na potrzeby regulacji<br />
wód oraz roślinności stanowiącej element zabudowy biologicznej dolin rzecznych lub<br />
służącej do wzmacniania brzegów, obwałowań lub odsypisk;<br />
3) zmiany ukształtowania terenu, składowania materiałów oraz wykonywania innych<br />
robót, z wyjątkiem robót związanych z regulacją lub utrzymywaniem wód oraz brzegu<br />
morskiego, a także utrzymywaniem, odbudową, rozbudową lub przebudową wałów<br />
przeciwpowodziowych wraz z obiektami związanymi z nimi funkcjonalnie (art. 88l ust. 1<br />
Prawo wodne).<br />
Podsumowując, na mapę ograniczeń składają się elementy (ich obrazy kolejno<br />
zamieszczone są w załączniku 4):<br />
• sieć dróg (zał. 4 a);<br />
• sieć torów kolejowych wraz z buforem 30 m (zał. 4 b);<br />
• wody tj. rzeki oraz zbiorniki wodne (zał. 4 c);<br />
75
• park narodowy (zał. 4 d);<br />
• otulina parku narodowego (zał. 4 e);<br />
• rezerwat (zał. 4 f);<br />
• tereny zalewowe z mapy stref zalewowych (zał. 4 g);<br />
• tereny zabudowane (zał. 4 h);<br />
• pomniki przyrody wraz z buforem 100 m (zał. 4 i);<br />
• strefy ochronne ujęć wody wraz z buforem 100 m (zał. 4 j).<br />
Mapa ograniczeń jest mapą o wartościach binarnych. Wartość 0 przypisano<br />
terenom stanowiącym barierę lokalizacyjną, natomiast wartość 1 takim, na których<br />
istnieje możliwość lokalizacji inwestycji (rys. 5.22).<br />
Rys. 5.22 Mapa ograniczeń.<br />
76
5.4. Lokalizacja ośrodka narciarskiego<br />
Poniżej zaprezentowano sposób postępowania, jaki zastosowano do wyszukania<br />
lokalizacji ośrodka narciarskiego. Najlepsze położenie takiego obiektu to połączenie<br />
dwóch sąsiadujących z sobą terenów o różnych uwarunkowaniach morfologicznych.<br />
Jeden o dużych spadkach, odpowiedni na stok zjazdowy oraz drugi płaski, na którym<br />
można umieścić całą towarzyszącą infrastrukturę jak hotel, parking.<br />
W związku z powyższym, etapy poszukiwania odpowiedniej lokalizacji<br />
podzielono na dwie części. Celem pierwszej było znalezienie stoków (patrz rozdz. 5.4.1,<br />
5.4.2 oraz 5.4.3). Wykorzystano w niej czynniki i bariery przygotowane w sposób<br />
opisany w powyższych podrozdziałach rozdziału 5. Celem drugiego było znalezienie<br />
terenów pod pozostałą zabudowę turystyczną. Analiza drugiej części wymagała nieco<br />
innej konfiguracji czynników (patrz rozdz. 5.4.4).<br />
5.4.1. Wagowanie<br />
Część pierwsza obejmowała dwie niezależne analizy MCE zróżnicowane pod<br />
względem znaczenia poszczególnych kryteriów. Aby wygenerować wagi kryteriów<br />
zastosowano metodę porównywania parami (AHP). Idea oraz sposób przypisywania wag<br />
w tej metodzie zostały opisane w rozdziale 3.3. Tabele prezentujące etap porównania<br />
kryteriów (na zasadzie „każdy z każdym”) zawiera załącznik 5 (dla analizy ze<br />
zwiększonym znaczeniem aspektu ekonomicznego) oraz załącznik 6 (dla analizy ze<br />
zwiększonym znaczeniem aspektu atrakcyjności turystycznej).<br />
Pierwszą analizę nazwano „ekonomiczną” z racji kierowania się w trakcie<br />
przypisywania wag surowymi zasadami opłacalności inwestowania. Skupiono się na tym,<br />
aby inwestor poniósł jak najmniejsze koszty w trakcie budowy, starano się również tak<br />
wagować by koszty eksploatacji były niskie.<br />
Drugi nazwany „turystycznym” to w rzeczywistości aspekt również w głównej<br />
mierze nastawiony na zysk, stawia jednak na atrakcyjność nie tylko samego terenu<br />
inwestycyjnego, ale również jego okolic. Uwzględnia potrzeby i wymagania<br />
potencjalnego klienta.<br />
77
O ile pierwszy aspekt ukierunkowano na zyski – oszczędności na kosztach, o tyle<br />
drugi przewiduje ewentualne poniesienie większych nakładów finansowych na budowę<br />
i eksploatację, upatruje natomiast zysków z wysokiej frekwencji użytkowników.<br />
W obu przypadkach istotną rolę odegrały czynniki środowiskowe. Przykładowo,<br />
czynnik przydatności rolniczej gleb wagowano, często przyrównując jego ważność do<br />
ważności innych kierunkowych czynników poszczególnych aspektów (załączniki 5 i 6<br />
przedstawiają etap porównań). Takie działanie ma na celu ochronę gleb wysokich klas<br />
oraz terenów leśnych przed degradacją.<br />
Wagowanie w analizie 1 – aspekt ekonomiczny<br />
W pierwszej analizie kryteria wagowano uwzględniając głównie aspekt<br />
ekonomiczny (opłacalność inwestowania), dlatego najistotniejszymi są:<br />
• odpowiednio długo zalegająca pokrywa śnieżna – w efekcie wydłużenie sezonu,<br />
mniejsze koszty związane z naśnieżaniem;<br />
• duże nachylenie – warunkuje możliwość powstania tras zjazdowych;<br />
• bliska odległość do dróg dojazdowych – łatwość dojazdu zapewnia większą<br />
frekwencję.<br />
Tańszym będzie również wykupienie działek mających status nierolniczy bądź<br />
nieleśny, zatem dużą wagę przypisano kryterium przydatności rolniczej gleb.<br />
Wzdłuż rzek występuje dużo stromych lecz krótkich, a zatem nieprzydatnych<br />
stoków, dlatego uznano również za ważny czynnik odległość od wód (zarówno dla<br />
aspektu ekonomicznego jak i turystycznego). Wysoką wagę uzasadniono wynikami<br />
przeprowadzonego eksperymentu w którym, czynnikowi odległość od wód przypisano<br />
małą wagę 0.02. Rysunki 5.23 a. i b. przedstawiają przykłady wyników takiej analizy<br />
o przydatności minimum 170 i powierzchni minimum 6 ha (te same dolne wartości<br />
zarówno przydatności jak i powierzchni występują we wszystkich analizach tej pracy).<br />
Widać na nich stoki w bezpośrednim sąsiedztwie rzek. Dla porównania przedstawiono<br />
przykładowe wyniki analiz: aspektu turystycznego w którym waga czynnika wynosiła<br />
0.0740 (rysunek 5.23 c – lokalizacja 1) oraz aspektu ekonomicznego, waga czynnika to<br />
0.0796 (rysunek 5.23 d – lokalizacja 3), w odniesieniu do lokalizacji względem rzek.<br />
Efektem zastosowania wysokich wag jest oddalenie wyników od rzek.<br />
78
a. b.<br />
c.<br />
d.<br />
Rys. 5.23 Przykładowe wyniki: a) i b) analizy eksperymentalnej w której czynnikowi<br />
odległość od wód przypisano niską wagę; c) analizy aspektu turystycznego oraz d) aspektu<br />
ekonomicznego w których czynnikowi przypisano wysoką wagę.<br />
Względnie dużą wagę nadano również terenom podmokłym, których osuszenie<br />
generowałoby dodatkowe koszty.<br />
Wagowanie w analizie 2 – aspekt turystyczno-społeczny<br />
Druga analiza jako ważniejsze uwzględniła czynniki określające atrakcyjność<br />
turystyczną terenu. Do takich zalicza się bliskość parków narodowych, rezerwatów<br />
i pomników przyrody. W założeniu idei bliskiej lokalizacji w stosunku do parku, klient<br />
(narciarz) nie tylko będzie mógł spędzić aktywnie czas na stoku, ale także pozwiedzać<br />
okoliczne tereny. Może to stanowić swego rodzaju rekompensatę w przypadku braku<br />
pogody narciarskiej.<br />
Większą wagę ma również czas dojazdu do mającego powstać obiektu. Bliskie<br />
sąsiedztwo może zachęcić narciarzy, dla których cenniejszy jest czas spędzony na stoku<br />
niż w podróży.<br />
79
Mimo znaczącego oddziaływania tych czynników najważniejszymi wciąż<br />
pozostały nachylenie stoku, długość zalegania pokrywy śnieżnej oraz bliska odległość od<br />
dróg dojazdowych. Ich duże wagi stanowią o opłacalności inwestycji.<br />
Poniższa tabela 5.7 zawiera zestawienie wag kryteriów obu analiz.<br />
Tab. 5.7 Wagi kryteriów generowane w oparciu o metodą porównywania parami (AHP).<br />
Czynnik<br />
Wagi analizy<br />
nr 1 (aspekt<br />
ekonomiczny)<br />
Wagi analizy<br />
nr 2 (aspekt<br />
turystyczny)<br />
Pokrywa śnieżna 0.2103 0.1638<br />
Nachylenie terenu 0.2091 0.2299<br />
Odległość od dróg kategorii G, Z, L 0.1644 0.1377<br />
Odległość od dróg kategorii A, S, GP 0.0148 0.0387<br />
Przydatność rolnicza gleb 0.1221 0.1022<br />
Formy ochrony przyrody 0.0336 0.0158<br />
Odległość od terenów podmokłych 0.0422 0.0116<br />
Odległość od wód 0.0796 0.0740<br />
Czas podróży 0.0282 0.0448<br />
Odległość od zwartej zabudowy 0.0166 0.0252<br />
Odległość od cmentarzy 0.0144 0.0409<br />
Ekspozycja 0.0433 0.0184<br />
Odl. od terenów chronionych (parki narodowe,<br />
rezerwaty)<br />
Odl. od pomników przyrody i obiektów<br />
dziedzictwa kulturowego<br />
0.0108 0.0526<br />
0.0106 0.0445<br />
Współczynnik spójności 0.10 0.07<br />
Największe zmiany w wagach pomiędzy dwoma analizami nastąpiły dla czynnika<br />
„pokrywa śnieżna” (różnica wynosi 0.0465), stosunkowo duża była również dla czynnika<br />
odległości od terenów podmokłych (0.0306). W obu tych przypadkach wagi większe<br />
cechowały aspekt ekonomiczny, wynika to z tego iż aspekt ekonomiczny to lokalizacja,<br />
której budowa i eksploatacja przewiduje jak najniższe koszty. Drugi z aspektów, zysków<br />
dopatruje się w innych możliwościach.<br />
Przewagę w ważności czynnika „czas podróży” zauważyć można w analizie<br />
drugiej. Ideą doboru wag analizy nr 2 „turystycznej” jest zachęcenie do korzystania<br />
z usług ośrodka, „wyposażając” go w dodatkowe walory, jakim jest np. krótszy czas<br />
dojazdu. Warto zauważyć również różnicę w wagach czynnika „odległość od cmentarzy”,<br />
jest to zmiana o 0.0265, także z przewagą dla analizy nr 2. Waga w pierwszej analizie<br />
80
odzwierciedla jedynie obowiązek prawny odsunięcia inwestycji od cmentarzy, wartość<br />
w drugiej to zapewnienie wyeliminowania widoku na cmentarz.<br />
Śledząc tabelę widać dużo wyższe wagi kryteriów „odległość od terenów<br />
chronionych” oraz „odległość od pomników przyrody i obiektów dziedzictwa<br />
kulturowego” w analizie aspektu turystycznego. Bliska odległość do parku narodowego<br />
czy pomników przyrody w zamierzeniach powinna stanowić dodatkową atrakcję,<br />
zachęcić do dłuższego pobytu.<br />
Bardzo podobne wagi można zauważyć w przypadku czynnika „odległość od<br />
wód” (różnica wag to jedynie 0.0056). Jest to czynnik o stałej ważności, nie zależy od<br />
założeń żadnego z aspektów.<br />
5.4.2. Analiza nr 1 – aspekt ekonomiczny<br />
Wynik analizy MCE nr 1 uwzględniający aspekt ekonomiczny przedstawia<br />
rysunek 5.24. Intensywnie czerwone obszary są terenami o najwyższej przydatności, dla<br />
których wartość wynosi 225 w skali od 0 do 255. Przydatność 200 charakteryzuje tylko<br />
kilka terenów o relatywnie małej powierzchni, nieprzydatnej pod inwestycję wyciągu<br />
narciarskiego. Dopiero dla przydatności o minimalnym progu 170 następuje<br />
wyselekcjonowanie odpowiednio dużych obszarów.<br />
Aby obliczyć powierzchnię wybranych obszarów zastosowano następujące<br />
funkcje:<br />
GROUP w celu pogrupowania obiektów,<br />
GIS Analysis → Context Operators → GROUP<br />
AREA to funkcja obliczająca powierzchnię pojedynczych obiektów.<br />
GIS Analysis → Database Query → AREA<br />
Stosując kalkulator obrazów (Image Calculator) wybrano obiekty powyżej 6 ha.<br />
Uznano to jako wystarczającą (minimalną) powierzchnię dla tego typu inwestycji.<br />
Proponowana wielkość nie odbiega bowiem od wielkości powierzchni stoków innych<br />
ośrodków zlokalizowanych blisko Krakowa np. ośrodek w Podstolicach, którego<br />
powierzchnia stoków zjazdowych niewiele przekracza 6 ha (zał. 7).<br />
81
Rys. 5.24 Wynik analizy MCE, mapa przydatności – aspekt ekonomiczny.<br />
Następnie zastosowano funkcję PERIM w celu wykluczenia spośród wyników<br />
terenów nieprzydatnych, o kształcie nieregularnym będących często stokami szerokimi,<br />
ale krótkimi.<br />
GIS Analysis → Database Query → PERIM<br />
Funkcja ta obliczyła dla wcześniej pogrupowanych wyników proporcję pomiędzy<br />
obwodem a powierzchnią obszarów, co jest równe wskaźnikowi regularności. Im wyższy<br />
wskaźnik tym obszar bardziej regularny. W tej analizie uznano jako odpowiedni<br />
o wartości 0.09 (tj. na tyle dużej by wyeliminować stoki najbardziej nieregularne –<br />
najniższa wartość 0.059).<br />
Rysunek 5.25 przedstawia efekt selekcji z użyciem PERIM. Dla łatwiejszej<br />
orientacji na mapie jako punkt odniesienia przedstawiono lokalizację Ojcowskiego Parku<br />
Narodowego. Jak widać, uzyskano 7 potencjalnych terenów pod inwestycję. Poniżej<br />
każdy z nich opatrzono szerszym komentarzem, charakteryzując je na podstawie ich<br />
82
manualnej analizy na tle map każdego z czynników, bądź weryfikacji na podkładzie<br />
mapy topograficznej.<br />
Wyniki końcowe w postaci map i najistotniejszych parametrów stoków<br />
przedstawiono w rozdziale podsumowującym (patrz rozdz. 5.6, tab. 5.8), gdzie znalazły<br />
się również wnioski z dodatkowej analizy związanej z możliwością lokalizacji w pobliżu<br />
infrastruktury obiektu.<br />
Rys. 5.25 Lokalizacja obszarów wynikowych analizy nr 1 (aspekt ekonomiczny).<br />
Lokalizacja nr 1 – aspekt ekonomiczny<br />
Teren wysunięty najdalej na zachód od Ojcowskiego Parku Narodowego, znajduje<br />
się w miejscowości Jerzmanowice. Największa różnica wzniesień wynosi 40 m,<br />
najdłuższy stok mierzy ok. 350 m. Spadki niewielkie w granicach 10 – 11 stopni, ostatnie<br />
100 m stoku charakteryzuje niewielki wzrost nachylenia do 14 stopni. Bliska odległość<br />
od drogi dojazdowej(ok. 100 m), długość zalegania pokrywy śnieżnej to od 94 do 98 dni.<br />
83
Czas dojazdu z centrum Krakowa to około pół godziny. Większa część obszaru zawiera<br />
się w glebach kompleksu 3 pszennego wadliwego, pozostała mniejsza część to kompleks<br />
2 pszenny dobry.<br />
Lokalizacja nr 2 – aspekt ekonomiczny<br />
Lokalizacja w wsi Wierzchowice, bliskie południowe sąsiedztwo z parkiem<br />
narodowym. Spadki wahają się między 12 a 18 stopniami. Różnica wzniesień wynosi 45<br />
stopni, najdłuższa możliwa trasa zjazdowa to około 300 m. Bardzo bliska odległość od<br />
drogi lokalnej, do drogi zbiorczej ok. 600 m. Zaleganie pokrywy śnieżnej to ok. 90 dni.<br />
Czas dojazdu – 35 minut. Cała powierzchnia znajduje na glebach kompleksu 3 oraz 14 są<br />
to gleby orne przeznaczone pod użytki zielone.<br />
Lokalizacja nr 3 – aspekt ekonomiczny<br />
Obszar w granicach wsi Iwanowice, wysunięty na wschód od parku narodowego.<br />
Około 70 % powierzchni terenu to gleby kompleksu 3 pszennego wadliwego, pozostała<br />
część - gleby kompleksu 10 orne przeznaczone pod użytki zielone. Długość zalegania<br />
pokrywy śnieżnej wynosi około 81 - 84 dni, stoki krótkie (maksymalna długość ok. 200<br />
m), spadki 11 – 16 stopni, maksymalna różnica wzniesień 35 m. Bliskość zarówno do<br />
dróg lokalnych jak i zbiorczych, 1,5 km do drogi głównej, czas dojazdu około pół<br />
godziny.<br />
Lokalizacja nr 4 – aspekt ekonomiczny<br />
Największy obszar, jego powierzchnia to aż 28 ha. Spadki są tutaj bardzo małe<br />
1 - 5 stopni, a więc teren całkowicie nieprzydatny pod tego rodzaju inwestycję. Jego duża<br />
przydatność (z analizy MCE) o wartości 170 wynika z dużych wag przypisanych dwóm<br />
najistotniejszym czynnikom z punktu opłacalności inwestycji: spadkom i czasowi<br />
zalegania pokrywy śnieżnej. Jest to obszar wysunięty daleko na północ, gdzie pokrywa<br />
śnieżna zalega nawet przez 95 dni. Duże znaczenie w przypisaniu wysokiej przydatności<br />
w analizie miała również wysoka waga czynnika przydatności rolniczej gleb, rozważany<br />
teren w 100 % zajmuje kompleksy o niskiej użyteczności rolniczej, stąd ich duża<br />
przydatność inwestycyjna (żytni słaby oraz żytni bardzo słaby).<br />
84
Lokalizacja nr 5 – aspekt ekonomiczny<br />
Kolejny teren wysunięty na północ, charakteryzuje się niewielkimi spadkami 3 do<br />
11 stopni, jego wysoka przydatność spowodowana jest podobnie jak w przypadku<br />
poprzednio opisywanego terenu długością zaleganiem pokrywy śnieżnej (95 dni). Stok<br />
o nachyleniu 9 do 11 stopni to 200 m potencjalnej trasy zjazdowej przy różnicy wzniesień<br />
około 20 m. Znajduje się blisko drogi zbiorczej, czas dojazdu z centrum Krakowa<br />
zajmuje godzinę.<br />
Lokalizacja nr 6 i nr 7 – aspekt ekonomiczny<br />
Obszary nr 6 i nr 7 sąsiadują ze sobą, są najbardziej wysuniętymi na północ.<br />
Długość zalegania pokrywy śnieżnej to nawet 98 dni. Oba charakteryzują się podobnymi<br />
rozmieszczeniami spadków 6 do 12 stopni, fragmenty na wyższej wysokości mają<br />
mniejsze spadki, im niżej nachylenie rośnie. Są to duże stoki 9 i 10,5 ha. Trasa zjazdowa<br />
mogłaby mierzyć 250 m ze spadkami 11-12 stopni przy różnicy wzniesień 35 m<br />
w przypadku obszaru nr 7, nr 6 to stok o długości 300 m ze spadkami 9-11 stopni, różnica<br />
wzniesień 40 m. Najbliższa droga dojazdowa jest oddalona o około kilometr. Czas<br />
dojazdu 1h 15 min. Oba tereny zalegają na glebach kompleksów 3, 4 i 6 o dużej<br />
przydatności inwestycyjnej.<br />
Spośród powyższych wyników najatrakcyjniejszym okazała się lokalizacja nr 2.<br />
Odznacza się najróżnorodniejszymi i największymi spadkami (12-18 stopni). Na tle<br />
innych charakteryzuje się dużą różnicą wzniesień (45 m) i długim stokiem (300 m).<br />
Kolejną atrakcyjną jest lokalizacja nr 1 - posiada mniejsze spadki (10-14), nieco mniejszą<br />
różnicę wzniesień (40 m), jest jednak dłuższym stokiem od poprzedniego o 50 m.<br />
85
5.4.3. Analiza nr 2 – aspekt atrakcyjności turystycznej<br />
Analiza nr 2 (rys. 5.26) wykazała najwyższą przydatność 234 w skali od 0 do<br />
255. Z tych samych względów, co w analizie nr 1 wyselekcjonowano jednak obszary<br />
o mniejszej przydatności równej bądź wyższej od 170 i o minimalnej powierzchni<br />
6 ha. Zastosowano również funkcję PERIM, w celu eliminacji stoków nieregularnych,<br />
szerokich ale krótkich.<br />
Rys. 5.26 Wynik analizy MCE, mapa przydatności – aspekt turystyczny.<br />
Analiza ta, przypisująca większą wagę czynnikom mającym wpływ na walory<br />
turystyczne, powieliła generalnie selekcję i wskazanie obszarów, jakie wystąpiło<br />
w analizie nr 1. Istotną zmianą są dwa dodatkowe tereny położone bliżej Krakowa.<br />
Rys. 5.27, przedstawia rozmieszczenie wyników, podobnie jak w poprzedniej<br />
analizie, dla łatwiejszej orientacji na mapie jako punkt odniesienia przedstawiono<br />
lokalizację Ojcowskiego Parku Narodowego. Poniżej podano charakterystykę nowych<br />
terenów.<br />
86
Rys. 5.27 Lokalizacja obszarów wynikowych analizy nr 2 (aspekt atrakcyjności<br />
turystycznej).<br />
Lokalizacja nr 1 – aspekt atrakcyjności turystycznej<br />
Obszar znajduje się we wsi Kozierów w gminie Michałowice. Jest szerokim<br />
stokiem o powierzchni 12,6 ha, spadkach 9 do 15 stopni. Przykładowa trasa zjazdowa to<br />
trasa o spadkach 10 przez 12 do 15 stopni o długości 300 m przy różnicy wzniesień 40 m.<br />
Pokrywa śnieżna zalega tu ok. 85 dni a gleby są kompleksów 2 i 3. Droga lokalna<br />
znajduje się kilka metrów od stóp zbocza, czas dojazdu z centrum Krakowa 20 minut.<br />
W analizie nr 1 teren nie został wybrany z uwagi na wyższą wagę czynnika<br />
zalegania pokrywy śnieżnej względem analizy nr 2 (różnica wag to prawie 22%).<br />
Lokalizacja ta charakteryzuje się krótko zalegającą pokrywą w porównaniu z wynikami<br />
poprzedniej analizy (w której prawie wszystkie wyniki cechują się ponad 90 dniowym<br />
zaleganiem).<br />
87
Lokalizacja nr 2 – aspekt atrakcyjności turystycznej<br />
Drugi obszar znajduje się we wsi Masłomiąca, w tej samej gminie. Jest to stok<br />
szeroki, ale krótki (ok. 150 m) o powierzchni 10 ha i spadkach 10-14 stopni. Droga<br />
lokalna znajduje się w odległości 500 m, czas dojazdu – 20 minut. Śnieg zalega tu 82 dni,<br />
wyznaczony obszar obejmuje gleby kompleksu 3.<br />
W analizie nr 1 teren nie został wybrany z uwagi na sporą różnicę wag czynnika<br />
długości zalegania pokrywy śnieżnej w porównaniu do analizy nr 2. Sytuacja podobna do<br />
opisanej w lokalizacji nr 1.<br />
Spośród opisanych wyników analizy aspektu turystycznego lepszą jest lokalizacja<br />
nr 1 gdyż charakteryzuje się znacznie dłuższym stokiem (długości 350 m podczas gdy<br />
drugi stok to jedynie 150 m).<br />
5.4.4. Napotkane nieprawidłowości, próby ich uniknięcia<br />
Przy okazji weryfikacji obszaru nr 2 analizy nr 2 (aspekt turystyczny) (rys. 5.28,<br />
c) na mapie topograficznej zauważono sąsiedztwo obszaru o atrakcyjniejszym nachyleniu<br />
i dłuższym stoku (rys. 5.28, b). Analiza MCE przypisała temu terenowi przydatność<br />
o wartości 200, jest więc interesującym z punktu widzenia inwestora. Jego kształt jest<br />
jednak bardzo nieregularny. Pojawiły się wątpliwości o nie do końca słusznym<br />
zastosowaniu funkcji PERIM gdyż ta przypisując mu wartość 0.090 wyeliminowała go<br />
uznając za stok nieprzydatny. Podobnie wykluczono wiele innych terenów, których<br />
eliminacja z kolei okazała się jednak słuszna, jak choćby przedstawiony na rysunku<br />
oznaczony literką „a” którego wartość PERIM wynosi 0.097. Jest to stok bardzo szeroki<br />
i krótki. Z kolei obszar „b” jest stokiem o długości 400 m, różnicy wzniesień 45 m<br />
i spadkach od 10 do 18 stopni. Pominięcie tego wyniku to utrata interesującego obszaru<br />
spełniającego wszelkie wymagania morfologiczne i lokalizacyjne.<br />
88
a.<br />
b.<br />
c.<br />
Rys. 5.28 Lokalizacja stoków analizy aspektu atrakcyjności turystycznej – efekty zastosowania<br />
funkcji PERIM.<br />
Aby uniknąć podobnych nieprawidłowości podjęto próbę zastosowania filtru<br />
modalnego. Jego działanie polega na obliczeniu w obszarze maski filtru, wartości<br />
modalnej – wartości najbardziej prawdopodobnej w otoczeniu piksela bieżącego. Efektem<br />
zastosowania filtru oraz funkcji MAX (w celu połączenia obrazu po filtracji i obrazu<br />
wejściowego) było połączenie dwóch blisko siebie leżących obszarów, bądź<br />
uregulowanie ich kształtów. Jednak częstym niepożądanym efektem było łączenie<br />
krótkich, szerokich stoków o przeciwnych ekspozycjach (przykład – rys. 5.29).<br />
Zastosowanie filtru nie przyniosło oczekiwanych rezultatów, dlatego z niego<br />
zrezygnowano.<br />
89
a.<br />
c.<br />
b.<br />
Rys. 5.29 Przykład niepożądanego efektu filtracji modalnej. a) obraz przed filtracją; b) obraz po<br />
filtracji; c) obraz po filtracji na podkładzie mapy topograficznej.<br />
Powyższe dwa przykłady przedstawiają problem braku uwzględnienia czynnika<br />
długości stoków. Niestety program IDRISI wykorzystywany do przeprowadzania analiz<br />
nie dysponuje funkcją umożliwiającą ich obliczanie. Konieczne byłoby zastosowanie<br />
innych narzędzi, jak program GRASS, jednak przekroczyło to zakres niniejszej pracy.<br />
Spośród sposobów rozwiązania problemu długości stoków najmniej zawodnym<br />
okazało się zastosowanie funkcji PERIM. Należy jednak traktować jej działanie jako<br />
wymagające kontroli ze strony użytkownika. Zawężona liczba wyników przedstawionych<br />
w pracy dla analiz obu aspektów jest rezultatem zastosowania właśnie tej funkcji.<br />
5.5. Analiza nr 3 – tereny infrastruktury turystycznej towarzyszącej<br />
stokowi narciarskiemu.<br />
Wykonanie dodatkowej analizy MCE nr 3 ma na celu znalezienie terenów<br />
odpowiednich na zagospodarowanie pod wypożyczalnie sprzętu narciarskiego, hotelu,<br />
restauracji i parkingu. Uwzględnia ona jedynie kilka istotnych kryteriów (ze zmianami<br />
ich wag w stosunku do poprzednich analiz a w niektórych przypadkach, jak spadków,<br />
również zmiany samej waloryzacji terenu) oraz wszystkie ograniczenia zastosowane<br />
w poprzednich analizach.<br />
90
Kryteria<br />
Mapy spadków, odległości od dróg klas G, Z, L oraz przydatności rolniczej gleb<br />
zostały przygotowane w taki sam sposób jak do poprzednich analiz MCE. Zmiany<br />
nastąpiły w sposobie standaryzacji w przypadku czynników: spadki oraz odległość od<br />
dróg klas G, Z, L.<br />
Mapa spadków<br />
Współczynnik opłacalności zostały przyporządkowany malejąco w stosunku do<br />
rosnącej wartości spadków. Użyto modułu FUZZY, obraz przeskalowano funkcją<br />
sigmoidalną monotonicznie malejącą o punktach przegięcia a = 6, b = 20 (zał. 8 a),<br />
oznacza to, że spadki poniżej 6 stopni mają najwyższy współczynnik przydatności 255.<br />
Standaryzacja tej analizy przypisuje najwyższą przydatność terenom płaskim,<br />
w poprzednich rzecz działa się przeciwnie. Przydatność wzrastała wraz z rosnącą<br />
wartością spadków.<br />
Mapa odległości od dróg klas G, Z, L<br />
Czynnik drogi G, Z, L standaryzowano funkcją sigmoidalną monotonicznie<br />
malejącą (Sigmoidal Monotonically Decreasing) o punktach przegięcia c=1, d=1000 (zał.<br />
8 b). Przydatność terenu maleje wraz z odległością do 1000 m, po przekroczeniu tej<br />
odległości przyjmuje wartość 0.<br />
W przypadku poprzednich analiz czynnik ten zestandaryzowano przypisując<br />
przydatność 0 terenom do 30 m, od 100 m do 1 km teren cechował się najwyższą<br />
przydatnością, malała do osiągnięcia odległości 1.5 km.<br />
Mapa przydatności rolniczej gleb<br />
analizach.<br />
Wartości przydatności kompleksów gleb pozostały takie same jak w poprzednich<br />
Z powodu małej ilości czynników nie zastosowano metody AHP do<br />
wygenerowania wag. Przypisano je na podstawie własnej wiedzy i doświadczeń.<br />
Najwyższą wagę nadano odległościom od dróg klas G, Z, L 0.4, pozostałym czynnikom<br />
po 0.3. Analiza wykazała wysoką przydatność dużej ilości obszarów. Istotnym było<br />
jednak znalezienie obszarów o dużej przydatności i o odpowiednio dużej powierzchni,<br />
a przy tym leżących w bliskim sąsiedztwie wybranych wcześniej stoków. Jako<br />
91
wystarczającą powierzchnię pod zabudowę tego typu obiektów uznano 0.8 ha. Z wyniku<br />
analizy wyselekcjonowano obszary o przydatności 170 lub wyższej.<br />
5.6. Kompleksowa waloryzacja wyselekcjonowanych obszarów<br />
W niniejszym rozdziale przedstawiono efekt poszukiwań najlepszej lokalizacji<br />
stoku oraz towarzyszącej mu infrastruktury. Każdy z dziewięciu wyników (siedem<br />
uzyskanych przy analizie aspektu ekonomicznego oraz dodatkowe dwa przy analizie<br />
aspektu turystycznego) przedstawiono na podkładzie mapy topograficznej w układzie<br />
1992 i 1965 w skali 1:10 000. Oprócz obszarów przydatnych jako stoki narciarskie<br />
i wyciągi, zaprezentowano powiązane z nimi obszary przeznaczone na niezbędną<br />
infrastrukturę dla tego typu przedsięwzięć (odpowiadające wymogom analizy nr 3)<br />
Lokalizacja nr 1 – aspekt ekonomiczny<br />
Rys. 5.30 Końcowy wynik analizy; lokalizacja nr 1 – aspekt ekonomiczny.<br />
92
Przedstawiony stok jest stokiem długim o spadkach 10-14 stopni (dokładniejszy<br />
opis stoku w rozdziale 5.4.2). Ponadto przylega do niego obszar odpowiedni pod<br />
zabudowę towarzyszącą (o powierzchni 0.8 ha) w bliskiej odległości od drogi dojazdowej<br />
co stanowi dodatkowy atut tej lokalizacji.<br />
Lokalizacja nr 2 – aspekt ekonomiczny<br />
Rys. 5.31 Końcowy wynik analizy; lokalizacja nr 2 – aspekt ekonomiczny.<br />
Jest to stok bardzo atrakcyjny pod względem nachylenia oraz różnicy wzniesień.<br />
Jego długość to 300 m, zlokalizowany jest bardzo blisko drogi dojazdowej (dokładniejszy<br />
opis w rozdziale 5.4.2). Do wyznaczonego stoku bezpośrednio przylega również teren<br />
odpowiedni pod infrastrukturę, wątpliwym pozostaje jednak jego lokalizacja na szczycie<br />
stoku.<br />
93
Lokalizacja nr 3 – aspekt ekonomiczny<br />
Rys. 5.32 Końcowy wynik analizy; lokalizacja nr 3 – aspekt ekonomiczny.<br />
Stok tej lokalizacji mierzy jedynie 200 m długości. Posiada natomiast inne walory<br />
jak np. spore spadki (nawet 16 stopni) czy bezpośrednio przylegające tereny odpowiednie<br />
pod zabudowę typu parkig, restauracja. Warto zwrócić uwagę, że spośród wyników<br />
analizy aspektu ekonomicznego charakteryzuje się najkrócej zalegającą pokrywą śnieżną.<br />
94
Lokalizacja nr 4 – aspekt ekonomiczny<br />
Rys. 5.33 Końcowy wynik analizy; lokalizacja nr 4 – aspekt ekonomiczny.<br />
Przedstawiony teren jest terenem całkowicie nieprzydatnym ze względu na bardzo<br />
małe nachylenie (1-5 stopni). Charakteryzuje się natomiast długo zalegającą pokrywą<br />
śnieżną (do 95 dni), jest to czynnik o bardzo dużej wadze. Pojawienie się obszaru<br />
niemalże płaskiego jest wynikiem przypisania prawie tak samo dużej wagi zarówno<br />
czynnikowi długości zalegania pokrywy śnieżnej jak i wysokim spadkom. Uznanie<br />
przydatności analizy MCE na poziomie 170 za wystarczającą nie eliminuje tego typu<br />
problematycznych wyników.<br />
Z uwagi na nieprzydatność terenu jako stok, na rysunku nie przedstawiono<br />
wyników analizy wyznaczającej tereny pod zabudowę pod pozostałą infrastrukturę<br />
turystyczną.<br />
95
Lokalizacja nr 5 – aspekt ekonomiczny<br />
Rys. 5.34 Końcowy wynik analizy; lokalizacja nr 5 – aspekt ekonomiczny.<br />
Obszar charakteryzuje się krótkim stokiem o niewielkich spadkach od 3 do 11<br />
stopni. Jego zaletą jest jedynie długo zalegająca pokrywa śnieżna oraz sąsiedztwo<br />
obszarów pod zabudowę towarzyszącą blisko drogi dojazdowej.<br />
96
Lokalizacja nr 6 i nr 7 – aspekt ekonomiczny<br />
Rys. 5.35 Końcowy wynik analizy; lokalizacja nr 6 i nr 7 – aspekt ekonomiczny.<br />
Tereny nr 6 i nr 7 przedstawiono na jednym rysunku dlatego, iż sąsiadują ze sobą,<br />
poza tym cechują się bardzo podobnymi parametrami. Podobnie jak na obszarach<br />
lokalizacji nr 4 i nr 5 śnieg zalega bardzo długo, na obydwu nawet 97 dni. Występują na<br />
nich spadki o małych wartościach głównie w przedziale 5-8 stopni (lokalnie 12 stopni),<br />
a czas dojazdu jest jednym z najdłuższych i wynosi 75 minut.<br />
Obszary pozostałej infrastruktury turystycznej w sporej części pokrywają się ze<br />
stokami. Warto zauważyć, że dokładnie te części stoków które nie pokrywają się<br />
z wynikami analizy dla infrastruktury są najbardziej interesujące pod względem spadków.<br />
Analiza nr 3 zadziałała w tym wypadku jako dodatkowy czynnik selekcjonujący obszary<br />
pod kątem atrakcyjności dla stoków.<br />
97
Lokalizacja nr 1 – aspekt atrakcyjności turystycznej<br />
Rys. 5.36 Końcowy wynik analizy; lokalizacja nr 1 – aspekt atrakcyjności turystycznej.<br />
Pierwszy z wyników analizy aspektu turystycznego to stok o spadkach nawet do<br />
15 stopni. Jest to atrakcyjny teren dodatkowo ze względu na sporą różnicę wysokości (40<br />
m) oraz długość stoku 300 m (dokładniejszy opis stoku w rozdziale 5.4.3). Posiada<br />
również atrakcyjne, przyległe tereny pod zabudowę typu parking, restaurację w bliskiej<br />
odległości od dróg dojazdowych. Warto zauważyć, że obszar ten znajduje się bardzo<br />
blisko Krakowa, czas dojazdu zajmuje jedynie 20 minut.<br />
98
Lokalizacja nr 2 – aspekt atrakcyjności turystycznej<br />
Rys. 5.37 Końcowy wynik analizy; lokalizacja nr 2 – aspekt atrakcyjności turystycznej.<br />
Przedstawiony stok charakteryzuje się spadkami w przedziale 10 – 14 stopni.<br />
Zlokalizowany jest blisko Krakowa, dodatkową zaletą jest bliskie sąsiedztwo terenów<br />
pod pozostałą zabudowę turystyczną. Zdecydowanie najgorszym parametrem,<br />
decydującym o nieatrakcyjności tej lokalizacji jest bardzo krótki stok mierzący zaledwie<br />
150 m.<br />
Tak krótki stok nie został wyeliminowany na etapie analizy MCE (w której<br />
przewidziano występowanie podobnych stoków wzdłuż rzek, dlatego czynnikowi<br />
odległości od wód przypisano dużą wagę), gdyż stok ten nie znajduje się w bliskiej<br />
odległości od rzeki. Nie wyeliminowano go również poprzez zastosowanie funkcji<br />
PERIM, gdyż posiadał zbyt duży współczynnik. Przyczyną pojawienia się takiego stoku<br />
wśród wyników jest brak modułu w programie IDRISI wyliczającego długość stoków.<br />
Ostatnim sposobem na eliminację tego typu niechcianych wyników jest weryfikacja na<br />
mapie topograficznej i odrzucenie ich przez analityka.<br />
99
Zbiorcza informacja na temat podstawowych parametrów wyszukanych<br />
lokalizacji stoku obu aspektów zestawiono w tabeli 5.8.<br />
Tab. 5.8 Podstawowe parametry wybranych stoków.<br />
Lokalizacja<br />
Nachylenie<br />
[°]<br />
Różnica<br />
wzniesień<br />
[m]<br />
Długość<br />
stoku [m]<br />
Długość<br />
zalegania<br />
pokrywy<br />
śnieżnej<br />
[dni]<br />
Powierzchnia<br />
[ha]<br />
Czas<br />
dojazdu<br />
[min]<br />
Analiza nr 1 – aspekt ekonomiczny<br />
lok. nr 1 10-14 40 350 94-98 7.5 30<br />
lok. nr 2 12-18 45 300 90 6.5 35<br />
lok. nr 3 11-16 35 200 81-84 6.5 30<br />
lok. nr 4 1-5 10 600 95 28 70<br />
lok. nr 5 3-11 20 200 95 8 60<br />
lok. nr 6 6-12 65 500 91-97 9 75<br />
lok. nr 7 6-12 45 400 94-98 10.5 75<br />
Analiza nr 2 – aspekt turystyczny<br />
lok. nr 1 9-15 40 300 85 12.6 20<br />
lok. nr 2 10-14 30 150 82 10 20<br />
Spośród wyżej zaproponowanych lokalizacji za najlepszą uznano lokalizację nr 1<br />
aspektu ekonomicznego. Za tym wyborem przemawiają dobre parametry stoku,<br />
dodatkowo u jego podstawy w bardzo dobrej lokalizacji tj. blisko drogi dojazdowej,<br />
znajduje się teren pod pozostałą infrastrukturę turystyczną. Atutem tego obszaru jest<br />
także bardzo długo zalegająca pokrywa śnieżna.<br />
Lepszym stokiem pod względem parametrów jest stok lokalizacji nr 2 aspektu<br />
ekonomicznego. Bezpośrednio do niego przylega również teren pod pozostałą<br />
infrastrukturę, wątpliwym jest jednak jego położenie względem stoku, gdyż znajduje się<br />
na jego szczycie a nie jak w przypadku lokalizacji nr 1 u stóp zbocza.<br />
Interesującym wynikiem jest również lokalizacja nr 1 aspektu turystycznego. Stok<br />
posiada dobrą możliwość lokalizacji infrastruktury towarzyszącej (w pobliżu drogi<br />
dojazdowej oraz u dołu stoku). Wynik ten nieco gorzej wypada względem poprzednich<br />
100
pod względem długości zalegania pokrywy śnieżnej, najlepiej natomiast wypada, co do<br />
czasu dojazdu z Krakowa (zajmuje jedynie 20 minut).<br />
Zdecydowanym faworytem pozostaje więc lokalizacja nr 1 aspektu<br />
ekonomicznego, gdyż zapewnia możliwość lokalizacji obu niezbędnych elementów (tj.<br />
stoku oraz infrastruktury towarzyszącej), wspólnie stanowiących w pełni kompleksowy<br />
obiekt narciarski. Przy tym, jest kompromisem pomiędzy nieco mniej stromym stokiem,<br />
ale dłuższym i dłużej utrzymującym pokrywę śnieżną.<br />
Powyższe porównania najciekawszych terenów utwierdzają w przekonaniu, że<br />
pomimo nieuwzględnienia istotnego kryterium, jakim jest długość stoków,<br />
przeprowadzone analizy doprowadziły do znalezienia wartych zainteresowania<br />
lokalizacji.<br />
101
6. Podsumowanie i wnioski<br />
Po przeglądzie dostępnych informacji przestrzennych, do analizy wykorzystano dane<br />
charakteryzujące analizowany teren pod kątem przydatności rolniczej gleb, pokrycia<br />
użytkowego (zabudowa, drogi), pokrycia różnymi formami ochrony przyrody, rozmieszczenia<br />
obiektów kulturowych, pokrycia siecią rzeczną, wysokości topograficznej (NMT). Na<br />
podstawie mapy NMT oraz dysponując odpowiednimi narzędziami programu IDRISI<br />
uzyskano dodatkowe informacje o terenie: jego nachylenia, ekspozycja stoków oraz długość<br />
zalegania pokrywy śnieżnej.<br />
Przydatnymi danymi, których zabrakło, mogłyby okazać się dane z mapy temperatur.<br />
Uzupełniłyby one analizę o informacje na temat najzimniejszych obszarów, na których śnieg<br />
zalegałby najdłużej. Istotną informację stanowić by mogła również mapa potencjalnych<br />
użytkowników obiektu (rozmieszczenie ludności lub – jeszcze lepiej – rozmieszczenie<br />
ludności aktywnie uprawiających sporty zimowe). Z braku takich danych przyjęto jako<br />
jedynych odbiorców mieszkańców Krakowa.<br />
Od strony narzędziowej brakowało możliwości obliczenia długości stoków,<br />
z problemem tym próbowano radzić sobie na kilka sposobów m.in. zastosowano funkcję<br />
PERIM. Funkcja ta przypisując wartości obiektom na podstawie obliczeń ich obwodów,<br />
umożliwia porównanie ich między sobą pod względem regularności kształtów. To działanie<br />
w zamierzeniu miało doprowadzić do wyeliminowania stoków szerokich, lecz krótkich<br />
(których w analizowanym terenie występuje dużo). Zastosowanie tego sposobu w dużym<br />
stopniu poradziło sobie z problemem. Metoda ta nie jest jednak w pełni niezawodna.<br />
Potwierdziła to weryfikacja wyników, gdy napotkano na atrakcyjny obszar wyeliminowany<br />
przez jego zastosowanie. Zastosowanie funkcji PERIM to dobry kierunek w stronę<br />
automatyzacji selekcji, wymaga jednak dopracowania.<br />
Innym zastosowanym sposobem na wykluczenie nieprzydatnych stoków było<br />
zwiększenie (na etapie przygotowywania map czynników) buforu terenów nieprzydatnych<br />
w odległości od rzek, gdyż spora część tego typu terenów znajduje się wzdłuż nich. Pozostałe<br />
krótkie stoki eliminowano w trakcie weryfikacji na mapie topograficznej.<br />
102
Pomimo zestawienia wielu danych o terenie w kryteria warunkujące lokalizację stoku,<br />
to wśród wyników analizy MCE znalazło się kilka terenów, które są całkowicie nieprzydatne<br />
dla rozpatrywanego rodzaju inwestycji. Są to tereny o bardzo małych spadkach. Ich<br />
pojawienie się wśród wyników należy tłumaczyć tym iż przyznano porównywalne, duże wagi<br />
dwóm kryteriom: długości zalegania pokrywy śnieżnej oraz nachyleniu terenu (są to kryteria<br />
warunkujące opłacalność inwestycji). Oznacza to, że teren o niekoniecznie dużym spadku, ale<br />
z długo zalegającą pokrywą śnieżną stanowi ten sam stopień użyteczności, w tym przypadku<br />
na poziomie przydatności 170 analiz aspektu ekonomicznego i turystycznego. Teoretycznie<br />
problem takich wyników można rozwiązać dodając jako ograniczenie tereny o niższych<br />
spadkach niż 10 stopni. Metoda taka wyklucza jednak również takie stoki, w których<br />
miejscowo pojawia się mniejszy spadek rozbijający ciągłość powierzchni. Taki teren na<br />
dalszym etapie selekcji mógłby zostać wyeliminowany z powodu nieregularnych kształtów<br />
lub zbyt małego pola powierzchni.<br />
Lepszym podejściem wydaje się zastosowanie mniej restrykcyjnych warunków.<br />
Analiza wyłania wówczas większą ilość wyników, wśród których najwłaściwszą selekcję<br />
zapewni weryfikacja analityka. Ten oceni, w jakim stopniu każdy wynik spełnia wymagania<br />
poszczególnych kryteriów. Należy podkreślić, że takie działanie ma sens jedynie, gdy ilość<br />
wyników nie jest nadmiernie duża, czyli w przypadku niewielkiego terenu, jak ten<br />
analizowany w niniejszej pracy. W przypadku większych obszarów należy postawić na<br />
bardziej rygorystyczne warunki i w ten sposób ograniczyć liczbę wyników.<br />
Wyniki przeprowadzonych dwóch analiz (pierwsza – aspekt ekonomiczny, druga –<br />
uwarunkowania turystyczne) nie różnią się od siebie w znaczący sposób. W obu<br />
uwzględniono wszystkie opracowane czynniki, zmieniano jedynie ich wagi. Analizowany<br />
obszar analiz jest zbyt mały by zauważyć znaczące różnice. Dodatkowo charakter inwestycji<br />
z góry przesądza o dużych wagach czynników morfologicznych, toteż spora część terenu<br />
zostaje odrzucona z powodu małych nachyleń w każdym z rozpatrywanych wariantów analizy<br />
MCA. Najwięcej atrakcyjnych spadków znajduje się w obrębie Ojcowskiego Parku<br />
Narodowego i jego otulinie, które zostały wyłączone z możliwości inwestowania (w ich<br />
obszarze znajduje się aż 50 % stoków o nachyleniu powyżej 10 stopni dla przyjętego obszaru<br />
badań). Na pozostałe nieliczne stoki w mniejszym stopniu oddziaływała reszta czynników,<br />
stąd kilka pokrywających się wyników.<br />
103
Przeprowadzona została również trzecia analiza MCE, której celem było znalezienie<br />
terenów będących odpowiednimi pod zabudowę typu parking, hotel, restauracja, czy<br />
wypożyczalnia sprzętu narciarskiego. Było to dodatkowe kryterium dla lokalizacji ośrodka<br />
narciarskiego. Ważnym i istotnie różnym od poprzednich analiz czynnikiem było<br />
preferowanie niewielkiego spadku, którego to występowanie stwarza dogodniejsze warunki<br />
budowy. Brak takich obszarów w pobliżu wybranych stoków nie jest jednak warunkiem<br />
eliminującym, gdyż nawet na terenach o dużym spadku można otrzymać pozwolenie na<br />
budowę obiektów gastronomicznych itp.<br />
Podsumowując system GIS to możliwość analizowania pod kątem określonych<br />
inwestycji dużych obszarów, możliwość uwzględnienia wielu czynników mogących<br />
w mniejszym bądź większym stopniu wpływać na ich lokalizację. Pomimo zastosowania<br />
zaawansowanych technik systemu GIS ważnym okazuje się również weryfikacja<br />
otrzymanego wyniku przez analityka, odniesienie go do rzeczywistości.<br />
System GIS może więc stanowić doskonałe narzędzie selekcji informacji na etapie<br />
wstępnym, natomiast ostateczne decyzje lokalizacyjne wymagają wnikliwej analizy wielu<br />
aspektów trudnych na dzień dzisiejszy do pełnego uwzględnienia metodami automatycznymi.<br />
104
7. Literatura<br />
Bibliografia<br />
<br />
<br />
<br />
<br />
<br />
Bielecka E., 2006 „Systemy informacji geograficznej. Teoria i zastosowanie”;<br />
Polsko-Japońska Wyższa Szkoła Technik Komputerowych<br />
Bohatkiewicz J., 2008 „Podręcznik dobrych praktyk wykonywania opracowań<br />
środowiskowych dla dróg krajowych” wersja 1.1;<br />
http://edroga.pl/images/stories/os_procedury/podrecznik_dobrych_praktyk/podrec<br />
znik_dobrych_praktyk.pdf dnia 28.06.2011.<br />
Brożek M., 2011 „Krajowy System Informacji Geograficznej”, materiały dla<br />
studentów; http://matrix.ar.krakow.pl/~mbrozek/materialy/sip/KSIG.pdf<br />
Dygaszewski J., 2006 „GEOPRTAL.GOV.PL jako system udostępniania danych<br />
z Państwowego Zasobu Geodezyjnego i Kartograficznego”;<br />
http://www.e-administracja.org.pl/dwumiesiecznik/artykul.phpart=36<br />
Eastman J. R., 2006 „IDRISI Andes Guide to GIS and Image Processing” Clark<br />
Labs Clark University, USA<br />
Gaździcki J., 2009 „INSPIRE – od idei do realizacji”<br />
http://poklinspire.gugik.gov.pl/files/prezentacje/02_gazdzickij_pokl_inspire_17gr<br />
ud2009.pdf<br />
<br />
<br />
<br />
Gotlib D., Iwaniak A., Olszewski R., 2007 „GIS Obszary zastosowań”<br />
Wydawnictw Naukowe PWN, Warszawa<br />
Hejmanowska B., 2004 „Ryzyko procesów decyzyjnych w aspekcie dokładności<br />
baz danych GIS” Akademia Górniczo – Hutnicza im. Stanisława Staszica,<br />
Kraków<br />
Hejmanowska B., 2005 „Wpływ jakości danych na ryzyko procesów decyzyjnych<br />
wspieranych analizami GIS” Rozprawy, monografie 141, <strong>AGH</strong> Uczelniane<br />
Wydawnictwa Naukowo-Dydaktyczne, Kraków<br />
105
Jarząbek J., 2011 „Infrastruktura informacji przestrzennej w zadaniach<br />
samorządów”, prezentacja z konferencji "Zainspiruj się! Infrastruktura informacji<br />
przestrzennej w jednostkach samorządowych", Warszawa<br />
http://www.akademiainspire.pl/konferencja-2011/prezentacje-z-konferencji/1-<br />
sesja/113-infrastruktura-informacji-przestrzennej-w-zadaniach-samorzdow<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
Kondracki J. A., 2000 „Geografia regionalna Polski” Wydawnictwo Naukowe<br />
PWN, Warszawa<br />
Kraszewski A., 2002 „Koncepcja systemu informacyjnego Inspekcji Ochrony<br />
Środowiska Ekoinfonet” Instytut Systemów Inżynierii Środowiska, Politechnika<br />
Warszawska; http://www.is.pw.edu.pl/plik/202/Ekoinfonet.pdf<br />
Kraszewski A., 2003 „System informacyjny Inspekcji Ochrony Środowiska - baza<br />
danych monitoringu powietrza JPOAT.” Konferencja Monitoring środowiska<br />
w samorządowych sieciach pomiarowych, Warszawa<br />
Longley P.A., Goodchild M.F., Maguire D.J, Rhind D.W., 2006 „GIS Teoria<br />
i praktyka” Wydawnictwo Naukowe PWN, Warszawa<br />
Nowosad M., 1992 „Uwagi do metody obliczania liczby dni z pokrywą śnieżną”;<br />
http://serwisy.umcs.lublin.pl/marek.nowosad/publications%20and%20other%20pa<br />
pers/1992%20Uwagi%20do%20metody%20obliczania%20liczby%20dni%20z%2<br />
0pokrywa%20sniezna.pdf<br />
Osada E., 2009 „Krajowy system informacji o terenie” Wydawnictwo Naukowe<br />
Dolnośląska Szkoła Wyższa,<br />
Preuss R., 2004a „Potrzeba budowy portalu obrazowego w Polsce” Archiwum<br />
Fotogrametrii, Kartografii i Teledetekcji vol. 14, Warszawa;<br />
http://www.sgp.geodezja.org.pl/ptfit/wydawnictwa/bialobrzegi/Bialobrzegi2004/5<br />
6-preuss.doc<br />
Preuss R., 2004b „Zakres zastosowań cyfrowej ortofotomapy w Polsce”<br />
Archiwum Fotogrametrii, Kartografii i Teledetekcji vol. 14, Warszawa;<br />
http://www.sgp.geodezja.org.pl/ptfit/wydawnictwa/bialobrzegi/Bialobrzegi2004/4<br />
9-preuss.doc<br />
106
Wężyk P., Mrugała M., Wńczyk R., Szwałko P., 2009 „Portal mapowy Zielony<br />
Kraków jako element realizacji INSPIRE” Polskie Towarzystwo Informacji<br />
Przestrzennej, Roczniki Geomatyki 2007, Tom VII, Zeszyt 6, Warszawa<br />
Netografia<br />
[1] http://www.akademiainspire.pl/dyrektywa-inspire/ustawa-o-iip - Ustawa<br />
o Infrastrukturze Informacji Przestrzennej<br />
[2] http://www.akademiainspire.pl/dyrektywa-inspire - Dyrektywa INSPIRE<br />
[3] http://www.akademiainspire.pl/najlepsze-gminy-akademii-inspire/zagorz -<br />
prezentacja studium przypadku „Wybór miejsca lokalizacji Biogazowni rolniczej<br />
w Gminie Zagórz” w ramach projektu „Geoinformacja w praktyce – Akademia<br />
INSPIRE”.<br />
[4] http://www.beskid.geo.uj.edu.pl/idrisi/ - Centrum IDRISI - Polska<br />
[5] http://www.bip.krakow.pl/sub_dok_id=15464&metka=1 - mapa podziału<br />
Krakowa na dzielnice<br />
[6] http://www.bip.mswia.gov.pl/download.phps=4&id=4514 - „Uzasadnienie<br />
projektu ustawy o infrastrukturze informacji przestrzennej”<br />
[7] http://www.biznes.gazetaprawna.pl/artykuly/325774,czy_na_obszarze_chronioneg<br />
o_krajobrazu_mozliwe_sa_komercyjne_inwestycje.html – Artykuł „Czy na<br />
obszarze chronionego krajobrazu możliwe są komercyjne inwestycje”,<br />
[8] http://www.codgik.gov.pl/baza-danych-ogolnogeograficznych.html - Baza<br />
Danych Ogólnogeograficznych<br />
[9] http://www.gddkia.gov.pl/userfiles/articles/p/posiedzenie-zespolu-ocenyprzeds_3124//documents/porwnanie-wariantw.pdf<br />
- informacje na temat metody<br />
porównania parami (AHP).<br />
[10] http://www.geoportal.gov.pl/index.phpoption=com_content&view=article&id=<br />
7&Itemid=20 – „Ramy prawne – Dyrektywa INSPIRE”<br />
[11] http://www.krakow.pios.gov.pl/publikacje/2010/ppms2010_2012.pdf - „Program<br />
państwowego monitoringu środowiska województwa małopolskiego na lata 2010-<br />
2012”<br />
107
[12] http://www.mapy.wrotamalopolski.pl/index.htm - informacje na temat portalu<br />
Wrota Małopolski<br />
[13] http://www.ski.podstolice.com.pl/pl/o_podstolice_ski,19,91 - strona ośrodka<br />
narciarskiego Podstolice-Ski<br />
[14] http://www.stat.gov.pl/cps/rde/xbcr/gus/PUBL_L_powierzchnia_ludnosc_teryt_<br />
2011.pdf; Główny Urząd statystyczny „Powierzchnia i ludność w przekroju<br />
terytorialnym w 2011 r.”, Warszawa<br />
[15] http://www.tatrzanska.com/archiwum/kotelnica/spolka.html - strona ośrodka<br />
narciarskiego Kotelnica Białczańska; kryteria lokalizacji stacji narciarskiej<br />
[16] http://widoczek.nets.pl/bieszczady/mezoregion/mapy/duze/mezoregiony_kondra<br />
ckiego.png - Mapa mezoregionów fizycznogeograficznych Polski na tle<br />
szczegółowego podziału administracyjnego.<br />
[17] http://www.wrotamalopolski.pl/root_SIOW/Region/stan+srodowiska/warunki+k<br />
limatyczne.htm – Informacje na temat warunków klimatycznych województwa<br />
małopolskiego.<br />
[18] http://www.zcpwz.pl/projekt/2010/n1125.html - „Modernizacja i aktualizacja<br />
ewidencji gruntów i budynków - aspekty formalne”, Departament Informacji<br />
o Nieruchomościach Głównego Urzędu Geodezji i Kartografii.<br />
108
Spis tabel i rysunków<br />
Tabele<br />
Tab. 4.1 Wykaz użytych danych i ich źródeł do opracowania kryteriów. ......................... 45<br />
Tab. 5.1 Podział czynników ............................................................................................... 49<br />
Tab. 5.2 Podział barier ....................................................................................................... 50<br />
Tab. 5.3 Koszty jednostkowe elementów pokrycia terenu. ............................................... 51<br />
Tab. 5.4 Ekspozycje i współczynniki ich opłacalności. ..................................................... 59<br />
Tab. 5.5 Kompleksy przydatności rolniczej ...................................................................... 63<br />
Tab. 5.6 Minimalne odległości obiektów budowlanych od zewnętrznej krawędzi jezdni<br />
(Ustawa o drogach publicznych). ...................................................................................... 64<br />
Tab. 5.7 Wagi kryteriów generowane w oparciu o metodą porównywania parami (AHP).<br />
............................................................................................................................................ 80<br />
Tab. 5.8 Podstawowe parametry wybranych stoków....................................................... 100<br />
Rysunki<br />
Rys. 2.1 Schemat działania Dyrektywy (Jarząbek, 2011). ................................................ 10<br />
Rys. 2.2 Struktura koordynacyjna (http://www.akademiainspire.pl/[1])........................... 13<br />
Rys. 2.3 Struktura tworzenia geoinformacji (Jarząbek, 2011)........................................... 14<br />
Rys. 2.4 Składniki Krajowego Systemu Informacji Geograficznej (Brożek, 2011). ......... 16<br />
Rys. 2.5 Idea interoperacyjności pomiędzy różnymi danym prowadzonych na różnych<br />
poziomach (Osada, 2009). ................................................................................................. 16<br />
Rys. 2.6 Schemat działania Ekoinfonetu (Kraszewski, 2002). .......................................... 22<br />
Rys. 2.7 Podsystem zbierania danych, wprowadzenie danych do KRDIOŚ (Kraszewski,<br />
2002). ................................................................................................................................. 23<br />
Rys. 2.8 Ekoinfonet Osady – przykład geoportalu tematycznego zawierającego dane<br />
z jednego z podsystemów PMŚ (zakres udostępnianych danych). .................................... 24<br />
Rys. 2.9 Lokalizacja scenariuszy na tle obszarów chronionych<br />
(http://www.akademiainspire.pl/[3]).................................................................................. 27<br />
109
Rys. 2.10 Najkorzystniejsza lokalizacja biogazowni na tle katastru<br />
(http://www.akademiainspire.pl/[3]).................................................................................. 28<br />
Rys. 2.11 Najkorzystniejsza lokalizacja biogazowni na tle projektu Studium Zagrożenia<br />
Powodziowego (http://www.akademiainspire.pl/[3])........................................................ 28<br />
Rys. 3.1 Rodzaje funkcji sigmoidalnej (Eastman, 2006). .................................................. 32<br />
Rys. 3.2 Rodzaje funkcji „j-shaped” (Eastman, 2006). ..................................................... 32<br />
Rys. 3.3 Rodzaje funkcji liniowej (Eastman, 2006). ......................................................... 33<br />
Rys. 3.4 Funkcja zdefiniowana przez użytkownika (Eastman, 2006). .............................. 33<br />
Rys. 3.5 Kompromis w zależności od ryzyka (Hejmanowska, 2004). .............................. 34<br />
Rys. 3.6 Problem niedokładności w procesie podejmowania decyzji (Hejmanowska,<br />
2004). ................................................................................................................................. 35<br />
Rys. 3.7 Schemat hierarchiczny w metodzie Saaty’ego (https://www.gddkia.gov.pl/ [9]).<br />
............................................................................................................................................ 37<br />
Rys. 4.1 Podział administracyjny analizowanego terenu. .................................................. 42<br />
Rys. 4.2 Wielowarstwowa kompozycja map w programie GeoMedia. ............................. 48<br />
Rys. 5.1 Zastosowanie algorytmu COST. .......................................................................... 52<br />
Rys. 5.2 Mapa kosztów. ..................................................................................................... 53<br />
Rys. 5.3 Czas dojazdu [min].............................................................................................. 54<br />
Rys. 5.4 Czas dojazdu – przydatność................................................................................. 54<br />
Rys. 5.5 NMT po filtracji uśredniającej (obszar gm. Michałowice). ................................ 56<br />
Rys. 5.6 Formy wklęsłe (odcienie zieleni) i wypukłe (odcienie żółci i czerwieni) (obszar<br />
gm. Michałowice). ............................................................................................................. 56<br />
Rys. 5.7 Wysokości form wklęsłych (obszar gm. Michałowice). ..................................... 57<br />
Rys. 5.8 Wysokości form wypukłych (obszar gm. Michałowice). .................................... 57<br />
Rys. 5.9 Ilość dni ze stałą pokrywą śnieżną. ...................................................................... 58<br />
Rys. 5.10 Ilość dni ze stałą pokrywą śnieżną (obszar gm. Michałowice). ........................ 58<br />
Rys. 5.11 Zaleganie pokrywy śnieżnej – mapa przydatności. ........................................... 58<br />
Rys. 5.12 Mapa ekspozycji stoków względem stron świata. ............................................. 60<br />
Rys. 5.13 Filtracja mapy - etap tworzenia mapy kompleksów gleb GeoMedia Professional<br />
............................................................................................................................................ 62<br />
Rys. 5.14 Mapa przydatności rolniczej gleb. ..................................................................... 64<br />
Rys. 5.15 Sieć dróg, centrum Krakowa i obszar na północny zachód od Krakowa. ......... 65<br />
Rys. 5.16 Odległość od dróg A, S, GP............................................................................... 67<br />
110
Rys. 5.17 Odległość od dróg A, S, GP – przydatność. ...................................................... 67<br />
Rys. 5.18 Odległość od dróg G, Z, L – przydatność. ......................................................... 67<br />
Rys. 5.19 Mapa form ochrony przyrody. ........................................................................... 68<br />
Rys. 5.20 Rozmieszczenie zabudowy zwartej i luźnej. ..................................................... 71<br />
Rys. 5.21 Rozmieszczenie elementów mapy wód w południowo-wschodniej części<br />
analizowanego obszaru. ..................................................................................................... 72<br />
Rys. 5.22 Mapa ograniczeń. ............................................................................................... 76<br />
Rys. 5.23 Przykładowe wyniki: a) i b) analizy eksperymentalnej w której czynnikowi<br />
odległość od wód przypisano niską wagę; c) analizy aspektu turystycznego oraz d)<br />
aspektu ekonomicznego w których czynnikowi przypisano wysoką wagę. ...................... 79<br />
Rys. 5.24 Wynik analizy MCE, mapa przydatności – aspekt ekonomiczny. .................... 82<br />
Rys. 5.25 Lokalizacja obszarów wynikowych analizy nr 1 (aspekt ekonomiczny). ......... 83<br />
Rys. 5.26 Wynik analizy MCE, mapa przydatności – aspekt turystyczny. ....................... 86<br />
Rys. 5.27 Lokalizacja obszarów wynikowych analizy nr 2 (aspekt atrakcyjności<br />
turystycznej). ...................................................................................................................... 87<br />
Rys. 5.28 Lokalizacja stoków analizy aspektu atrakcyjności turystycznej – efekty<br />
zastosowania funkcji PERIM. ............................................................................................ 89<br />
Rys. 5.29 Przykład niepożądanego efektu filtracji modalnej. a) obraz przed filtracją; b)<br />
obraz po filtracji; c) obraz po filtracji na podkładzie mapy topograficznej. ...................... 90<br />
Rys. 5.30 Końcowy wynik analizy; lokalizacja nr 1 – aspekt ekonomiczny. ................... 92<br />
Rys. 5.31 Końcowy wynik analizy; lokalizacja nr 2 – aspekt ekonomiczny. .................... 93<br />
Rys. 5.32 Końcowy wynik analizy; lokalizacja nr 3 – aspekt ekonomiczny. .................... 94<br />
Rys. 5.33 Końcowy wynik analizy; lokalizacja nr 4 – aspekt ekonomiczny. .................... 95<br />
Rys. 5.34 Końcowy wynik analizy; lokalizacja nr 5 – aspekt ekonomiczny. .................... 96<br />
Rys. 5.35 Końcowy wynik analizy; lokalizacja nr 6 i nr 7 – aspekt ekonomiczny. .......... 97<br />
Rys. 5.36 Końcowy wynik analizy; lokalizacja nr 1 – aspekt atrakcyjności turystycznej. 98<br />
Rys. 5.37 Końcowy wynik analizy; lokalizacja nr 2 – aspekt atrakcyjności turystycznej. 99<br />
111
Załączniki<br />
Zał. 1 Rozmieszczenie wartości tarcia w obszarze gminy Michałowic.<br />
a. Mapa pokrycia. b. Drogi. c. Mapa barier (rzeki, kolej).<br />
Zał. 2 Import danych formatu SHAPE do formatu programu.<br />
112
Zał. 3 Mapy przydatności czynników.<br />
a. Ekspozycja. b. Spadki.<br />
c. Tereny podmokłe. d. Formy ochrony przyrody.<br />
e. Odległość od terenów chronionych. f. Odległość od pomników przyrody.<br />
113
g. Odległość od zabudowy zwartej h. Odległość od wód.<br />
i. Odległość od cmentarzy. j. Przydatność rolnicza gleb.<br />
114
Zał. 4 Mapy ograniczeń.<br />
a. Ograniczenie – sieć dróg. b. Ograniczenie – sieć kolei.<br />
c. ograniczenie – wody. d. Ograniczenie – park narodowy.<br />
e. Ograniczenie – otulina parku narodowego. f. Ograniczenie – rezerwaty przyrody.<br />
115
g. Ograniczenie – teren zalewowy. h. Ograniczenie – zabudowa.<br />
i. Ograniczenie – pomniki przyrody. j. Ograniczenie - strefy ujęć wody.<br />
Zał. 5 Porównywanie czynników parami (AHP) – aspekt ekonomiczny.<br />
pok_sn<br />
iez<br />
spadki_<br />
st_trud<br />
g_z_l<br />
pok_sniez 1<br />
spadki_st_tru<br />
d<br />
1 1<br />
g_z_l 1/3 1/3 1<br />
a_s_gp<br />
przydat<br />
_rol_gle<br />
b<br />
ter_chro<br />
nione<br />
a_s_gp 1/9 1/9 1/9 1<br />
przydat_rol_g<br />
leb<br />
1/3 1/3 1/3 9 1<br />
ter_chronione 1/7 1/7 1/7 5 1/7 1<br />
ter_pod<br />
m<br />
ter_podm 1/7 1/6 1/6 3 1/5 3 1<br />
wody 1/5 1/6 1/5 7 1/2 3 5 1<br />
m_dojazdu 1/9 1/7 1/7 5 1/7 1/3 1/3 1/5 1<br />
wody<br />
m_doja<br />
zdu<br />
zab_zwar 1/9 1/9 1/9 1/2 1/9 1/5 1/3 1/7 1/3 1<br />
zab_zw<br />
ar<br />
cmentarz 1/9 1/9 1/7 1/3 1/7 1/3 1/5 1/5 1/3 1/3 1<br />
cmentar<br />
z<br />
ekspozycja 1/7 1/7 1/7 5 1/6 3 1 1/3 3 3 5 1<br />
ekspozy<br />
cja<br />
parki_nar 1/9 1/9 1/9 1 1/9 1/3 1/5 1/7 1/5 1/3 1/3 1/4 1<br />
pom_przyr 1/9 1/9 1/9 1 1/9 1/3 1/5 1/7 1/5 1/3 1/3 1/5 1 1<br />
parki_n<br />
ar<br />
pom_pr<br />
zyr<br />
116
Zał. 6 Porównywania czynników parami (AHP) – aspekt atrakcyjności turystycznej.<br />
pok_sni<br />
ez<br />
spadki_<br />
st_trud<br />
g_z_l<br />
pok_sniez 1<br />
spadki_st_tru<br />
d<br />
3 1<br />
g_z_l 1 1/3 1<br />
a_s_gp<br />
przydat<br />
_rol_gle<br />
b<br />
ter_chro<br />
nione<br />
a_s_gp 1/6 1/5 1/5 1<br />
przydat_rol_g<br />
leb<br />
1/4 1/3 1/4 4 1<br />
ter_chronione 1/7 1/7 1/7 1/4 1/6 1<br />
ter_pod<br />
m<br />
ter_podm 1/7 1/7 1/7 1/6 1/6 2 1<br />
wody 1/3 1/5 1/3 2 1/3 4 6 1<br />
m_dojazdu 1/4 1/5 1/3 1 1/3 5 6 1/3 1<br />
wody<br />
m_dojaz<br />
du<br />
zab_zwar 1/7 1/7 1/3 1/4 1/5 3 4 1/3 1/3 1<br />
zab_zw<br />
ar<br />
cmentarz 1/6 1/6 1/3 1 1/3 4 4 1/3 1 3 1<br />
cmentar<br />
z<br />
ekspozycja 1/6 1/7 1/6 1/3 1/6 3 4 1/3 1/3 2 1/3 1<br />
ekspozy<br />
cja<br />
parki_nar 1/4 1/6 1/3 1 1/3 5 5 1/3 1 5 1 3 1<br />
pom_przyr 1/4 1/6 1/3 1 1/3 5 5 1/3 1 5 1 3 1/3 1<br />
parki_n<br />
ar<br />
pom_pr<br />
zyr<br />
Zał. 7 Uzasadnienie uznania 6 ha powierzchni stoku jako wystarczające.<br />
a.<br />
b.<br />
Krawędzie mierzonej powierzchni (a.) wyznaczono na podstawie rysunku stoków (b.)<br />
zamieszczonych na stronie ośrodka narciarskiego w Podstolicach. Jest to jedynie szacowana wielkość<br />
powierzchni (oficjalna strona obiektu nie podaje takiej informacji). Pomiaru dokonano na mapie<br />
topograficznej pochodzącej z geoportalu.<br />
117
Zał. 8 Mapy przydatności czynników dla analizy nr 3.<br />
a. Spadki<br />
b. Odległość od dróg kategorii G, Z, L.<br />
118