13.02.2015 Views

PRACA DYPLOMOWA - AGH

PRACA DYPLOMOWA - AGH

PRACA DYPLOMOWA - AGH

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!