Fizyczne Podstawy Teledetekcji Wykład 7

1. Fizyczne Podstawy TeledetekcjiWykład 7 Krzysztof Markowicz Instytut Geofizyki Uniwersytet Warszawski kmark@igf.fuw.edu.pl

2. Teledetekcja powierzchni Ziemi Rozpoznawanie typów podłoża na podstawie współczynników odbicia promieniowania w zakresie widzialnym i bliskiej podczerwieni. Współczynnik odbicia podłoża zależy nie tylko od rodzaju powierzchni ale również od położenia Słońca oraz detektora satelitarnego lub lotniczego. W ogólnym przypadku są to dwa kąty zenitalne i dwa kąty azymutalne. Dlatego wprowadza się pojęcie BRDFu - dwukierunkowego współczynnika odbicia.

3. Współczynnik 2-kierunkowego odbicia w przeciwieństwie od albeda zdefiniowanego dla strumieni promieniowania nie zależy od własności optycznych atmosfery a jedynie od własności samej powierzchni odbijającej. (,)  d o o(,) Definicja BRDF-u na górnej granicy atmosfery Definicja BRDF-u powierzchni ziemi 3

4. Większość typów powierzchni ziemi wykazuje własności optyczne pomiędzy dwoma skrajnymi typami: idealnie gładka (odbicie zwierciadlane) oraz powierzchnia szorstką (Lambertsowska) W pierwszym przypadku współczynnik odbicia jest niezerowy jedynie dla kąta padania równego kątowi odbicia (odbicie Fresnela) W drugim przypadku radiancja promieniowania odbitego jest izotropowa. Określenie BRDF-u powierzchni ziemi jest kluczowe dla większości metod teledetekcyjnych. Czasami funkcję BRDF definiuje się dla całego zakresu falowego promieniowania. 4

5. Odbicie promieniowania 5

6. 6

7. Zmienność spektralna współczynnika odbicia powierzchni ziemi

8. ASTER baza danych współczynników odbicia podłoża http://speclib.jpl.nasa.gov/ Typy podłoża i ilość dostępnych ich radzajów Minerals (1748) Lunar (17) Vegetation(4) Rocks(473) Meteorites (60) Water/Snow/Ice(9) Soils (69)Man-Made (84)

11. Brown loamy fine sand

14. Zmiany współczynnika odbicia w bliskiej podczerwieni (NIR) opisywane są przez indeks NDVI (znormalizowany, różnicowy indeks wegetacyjny) I(NIR) oraz I(VIS) są radiancja promieniowania odbitego w bliskiej podczerwieni i w obszarze widzialnym. NDVI mówi nam o zawartości składnika wegetacyjnego w danym typie podłoża. Im większy indeks tym zawartość biomasy większa. Przykład W przyrządzie AVHRR(Advance Very Heigh Resolution Radiometer) w celu wyznaczenia NDVI wykorzystuje się kanał 1 (0.54-0.68 m) oraz kanał 2 (0.73-1.10 m).

15. Typowe wartości indeksu NDVI

16. Korekcja atmosferyczna • Wartości I(NIR) oraz I(VIS) są mierzone na górnej granicy atmosfery a więc zawierają również przyczynek od atmosfery • Imeas=Isurf+I* • I*poprawka atmosferyczna często ogranicza się tylko do rozpraszania Rayleigha

17. The enhanced vegetation index (EVI) (NIR), (RED), (BLUE) – skorygowane na rozpraszanie Rayleigha i absorpcję przez ozon współczynniki odbicia w kanałach: bliskiej podczerwieni, czerwonym oraz niebieskim. Współczynniki C1, C2,, L uwzględniają wpływ aerozolu i dla przyrządu MODIS-EVI wynoszą odpowiednio 6.0, 7.5 oraz 1.0. Parametr G nosi nazwę współczynnika wzmocnienia i wynosi 2.5.

18. Zoptymalizowany indeks EVI został opracowany w celu zwiększenia czułości detekcji powierzchni wegetacyjnych. Indeks NDVI jest czuły na zawartość chlorofilu w podłożu podczas gdy EVI na różne rodzaje podłoża w tym indeksu LAI. Ponadto EVI minimalizuje wpływ atmosfery w porównaniu do klasycznego indeksu NDVI.

19. Indeks LAI (Leaf Area Index) • Jest zdefiniowany jako stosunek górnej powierzchni liści do całkowitej powierzchni. Współczynnik przyjmuje wartości od 0 dla gleby do 6 dla gęstego lasu. • Parametr te zależy od powierzchni oraz kształtu liści i jest związany z całkowitą zawartością biomasy.

20. EVI

25. CERES/SARB własności podłożahttp://www-surf.larc.nasa.gov/surf/pages/explan.html rozdzielczość spektralna bazy danych to 1/6o (około 20 km) Typy powierzchni 1. Evergreen Needleleaf Forest2. Evergreen Broadleaf Forest3. Deciduous Needleleaf Forest4. Deciduous Broadleaf Forest5. Mixed Deciduous Forest6. Closed Shrubland7. Open Shrubland8. Woody Savanna9. Savanna10. Grassland11. Permenant Wetland12. Cropland13. Urban14. Crop/Natural Veg. Mosaic15. Permanent Snow/Ice16. Barren/Desert17. Water Bodies18. Tundra19. Fresh Snow20. Sea Ice

26. International Geosphere/Biosphere Programme (IGBP) scene types

27. Albedo powierzchni

28. Zdolność emisyjna

29. Konieczność walidacji • Pomiary zdalne obarczone są błędami, które wynikają z wpływu atmosfery. Zmienność atmosfery wpływa na niepewności w wyznaczaniu współczynników odbicia. • Dlatego pomiary zdalne powierzchni ziemi wymagają walidacji przy użyciu pomiarów in-situ. • Prowadzi się ją: - przy pomocy urządzeń (spektrometrów) do pomiarów współczynników odbicia tuż nad powierzchnia ziemi. - przez analizę składu wycinka powierzchni ziemi. • - innymi metodami.

30. Korekcja własności podłoża ze względu na geometrię skanowania. • Pomiary pokazują, że dla danego położenia Słońca współczynniki odbicia mogą zmieniać się o czynnik 2 w obszarze bliskiej podczerwieni i więcej w obszarze widzialnym. • Zmienności współczynników odbicia powierzchni Ziemi wraz ze zmianą kątów padania promieniowania słonecznego oraz katów określających położenie detektorów satelitarnych jest źródłem problemów w standaryzacji pomiarów np. indeksu NDVI, EVI itd. • Wynika z tego, że NDVI czy EVI zależą nie tylko od własności podłoża ale również od położenia Słońca i detektora skanującego na orbicie (samolocie). • Wymaga to więc stosowania korekcji sygnału.

31. BRDF jako funkcja kąta zenitalnego i azymutalnego dla lasu iglastego[Schaaf et al.]

32. Parametryzacja BRDFu model Ross-Li-Maignan • s – kąt zenitalny Słońca  v – kąt zenitalny detektora  - względny kąt azymutalny Słońca i detektora, k0 , k1, k2 współczynniki, F1 określa cześć objętościową jądra rozpraszania z korekcją na tzw. hot spot, F2 opisuje cześć geometryczną. gdzie  jest kątem rozpraszania, zaś o opisuje szerokość kątową hot spot. Zakłada się, że wynosi ona 1.5 stopnia.

33. kąt rozpraszania:

34. Uwaga do geometrii… SAT v SUN • Ponieważ kąty w układzie sferycznymi licznymi są od zenitu ważny jest zwrot wektora opisującego położenie detektora satelitarnego oraz propagacje promieniowania słonecznego. • W ostatnim przypadku wektor skierowany jest ku ziemi co oznacza, że: s 180-s • Gdy satelita i Słońce znajdują po przeciwnych stronach (tak jak na rysunku) to wszystkich wzorach zamiast s powinno być 180-s , zaś =0. • Gdy satelita i Słońce znajdują się po tej samej stronie to wówczas zamiast s powinno być -s , zaś =180. • W przypadku innego ustawienia kąt  przybiera wartości z przedziału 0-180.

35. Kształt funkcji F1 i F2 Maksima na wykresie F2 opisują efekt hot spot

36. Bacour et al., 2005

37. Rozkłady prawdopodobieństwa dla parametrów k

38. Model Rahmana

39. W dalszej części będziemy posługiwać się tylko pierwszym modelem. Definiujemy współczynniki: Wyniki pomiarów pokazują, że kształt BRDF zmienia się znacznie wolniej niż jego amplituda. Pozwala to na oszacowanie parametrów V i R a następnie używanie ich do skalowania zmian BRDFu. Korekcja współczynników odbicia będzie wykonywana dla kąta zenitalnego Słońca 45 stopni oraz skanowania z kierunku nadiru.

40. Seria czasowa współ. odbicia (2000 to 2004) wykonanych przy użyciu MODISaw kanałach REDoraz NIR nad Afrykańska sawanną. Szum (noise) odpowiada fluktuacjom związanym z różną geometrią w czasie pomiaru.

41. Wyznaczanie BRDFu W klasycznym przybliżeniu zakłada się, że współ. odbicia nie zmienia się w czasie. minimalizacja: Nowe przybliżenie zaproponowane przez Vermote et al., 2009 uwzględnia powolną zmianę współ. odbicia w czasie. minimalizacja:

42. Minimalizacja sprowadza się do rozwiązania problemu liniowego w postaci: gdzie:

43. Dane nieskorygowane

44. Skorygowany współczynnik odbicia (sawanna tropikalna) przy użyciu klasycznej metody.

45. Skorygowane dane przy użyciu nowej metody

46. Parametry V i R jako funkcję indeksu NDVI Red: band 2 Blue: band 1

47. Wyniki z nową korekcją przy założeniu, że R i V zależą liniowo od NDVI.

48. Oryginalne dane NDVI

49. NDVI na podstawie klasycznej korekcji

50. NDVI na podstawie nowej korekcji ze stałymi wielkościami V i R.