Temat 4
(określenie wybranych wskaźników ilościowych z obrazów
satelitarnych)
Damian Nowak
Gr 5
Nr 4
Sprawozdanie techniczne
I.
Dane formalno-prawne
1.
ZLECENIODAWCA: Wydział Geodezji Górniczej i Inżynierii Środowiska Akademii Górniczo-
Hutniczej w Krakowie
2.
WYKONAWCA: Damian Nowak
3.
PRZEDMIOT ZLECENIA: Określenie wybranych wskaźników ilościowych z obrazów
satelitarnych
4.
TERMIN WYKONANIA ZLECENIA: 10.03.2011.
II.
Wykorzystane materiały
Obrazy satelitarne z satelity Landsat TM przedstawiające okolice Krakowa – odrębne obrazy
dla 8 kanałów spektralnych: niebieskiego, zielonego, czerwonego, bliskiej podczerwieni, średniej
podczerwieni, kanał termalny, kanał panchromatyczny.
III.
Wykonanie zlecenia
Przebieg ćwiczenia
o Stworzenie kompozycji barwnej dla kanałów RED(krak3), GREEN(krak2),
BLUE(krak1) i zapisanie jej pod nazwą krak123
o Za pomocą narzędzia Map Calculation dokonano konwersji jasności
zarejestrowanej przez detektor w satelicie na luminancję
o Ponownie za pomocą funkcji Map Calculation dokonano konwersji luminancji na
temperaturę radiacyjną korzystając z „odwróconego wzoru Plancka”.
Temperatura uzyskana została w stopniach Kelwina.
o Za pomocą narzędzia Map Calculation przeliczono temperaturę na stopnie
Celsjusza.
o Analiza mapy radiacyjnej
Pomiar temperatury w miejscach o różnym pokryciu terenu
Rodzaj użytkowania
terenu
Zakres
temperatury [°C]
Średnia
temperatura [°C]
Zabudowa zwarta
25-27
26
Zabudowa rozproszona
23-26.5
25
Użytki rolne
20-23
22
Lasy
19-21
20
Rzeki
17-20
18
Zbiorniki wodne
16-18
17
Odpowiedź na pytania:
Jakie tereny charakteryzują się wyższymi temperaturami? –
Wyższe temperatury notowane są w większych osiedlach
ludzkich, na terenach przemysłowych
Jakie tereny charakteryzują się niższymi temperaturami? –
Temperaturami niższymi charakteryzują się duże skupiska drzew,
lasy, woda stojąca i płynąca.
Obliczenie powierzchni obszaru, dla którego temperatura radiacyjna jest
wyższa niż 20.8 °C
Za pomocą Map Calculation dokonano odpowiedniej konwersji
zdjęcia satelitarnego
Przy pomocy histogramu odczytano potrzebne informacje
Powierzchnia = 45716679 m
2
Procent całości = 88,57%
o Opracowanie wskaźnika indeksu wegetacyjnego NDVI
Za pomocą Map Calculation przeprowadzono konwersję w oparciu o wzór:
o Analiza mapy wskaźnika wegetacyjnego NDVI
Pomiar wskaźnika NDVI w miejscach o różnym pokryciu terenu
Rodzaj użytkowania
terenu
Zakres wskaźnika
NDVI
Średnia wartosć
wskaźnika NDVI
Zbiorniki wodne
od -0.22 do -0.40
-0.31
Rzeki
od -0.30 do -0.40
-0.35
Starorzecza
od -0.30 do 0
-0.15
Użytki zielone
od 0.11 do 0.25
0.18
Lasy
od 0.36 do 0.57
0.46
Gleba odkryta
od 0 do 0.20
0.1
Zabudowa jednorodzinna od -0.10 do 0.20
0.05
Zabudowa przemysłowo-
handlowa
od -0.15 do -0.34
-0.25
Obliczenie powierzchni obszaru, dla którego wskaźnik NDVI jest większy
od 0.15
Za pomocą Map Calculation dokonano odpowiedniej konwersji
zdjęcia satelitarnego
Przy pomocy histogramu odczytano potrzebne informacje
Powierzchnia = 28821734 m
2
Procent całości = 56,94%
Do sprawozdania zostały dołączone następujące załączniki:
1.
Mapa luminancji
2.
Mapa temperatury radiacyjnej w stopniach Kelwina
3.
Mapa temperatury radiacyjnej w stopniach Celsjusza
4.
Mapa obszaru o temperaturze większej od 20.8 stopni Celsjusza
5.
Mapa wskaźnika wegetacyjnego NDVI
6.
Mapa obszaru o wskaźniku NDVI>0.15
Załącznik 1 – mapa luminancji
Załącznik 2 - Mapa temperatury radiacyjnej w stopniach Kelwina
Załącznik 3 - Mapa temperatury radiacyjnej w stopniach Celsjusza
Załącznik 4 – Mapa obszaru o temperaturze większej od 20.8 stopni Celsjusza
Załącznik 5 - Mapa wskaźnika wegetacyjnego NDVI
Załacznik 6 - Mapa obszaru o wskaźniku NDVI>0.15