QGIS, OrfeoToolbox i indeksy spektralne

     Niedawno pisałem o możliwości obliczenia indeksów wegetacyjnych (spektralnych) w Quantum GIS przy użyciu Kalkulatora Rastrów. W programie istnieje również inny sposób określenia indeksów. W aplikacji, mamy dostęp do potężnego narzędzia – Sextante, które jest zintegrowane z QGIS-em pod nazwą Geoprocesing (dawniej była to wtyczka). Wśród jego modułów znajduje się Orfeo Toolbox, dzięki któremu możemy korzystać z dużej liczby algorytmów w celu przetwarzania obrazów teledetekcyjnych (opis jak skonfigurować moduł Orfeo Toolbox znajduje się w poście Teledetekcja w QGIS). W OT obliczanie indeksów wegetacyjnych dla wczytanego zobrazowania zaimplementowano w Feature Extraction > Radiometric Indices.

  

Widoczne indeksy zostały podzielone na 3 grupy:

1. Wegetacja:

• NDVI - Normalized Difference Vegetation Index
• TNDVI – Transformed Normalized Difference Vegetation Index
• RVI - Ratio Vegetation Index
• SAVI - Soil Adjusted Vegetation Index
• TSAVI - Transformed Soil Adjusted Vegetation Index
• MSAVI - Modified Soil Adjusted Vegetation Index
• MSAVI2 - Modified Soil Adjusted Vegetation Index 2
• GEMI - Global Environment Monitoring Index
• IPVI - Infrared Percentage Vegetation Index
• LAIFromNDVILog - Leaf Area Index from NDVI
• LAIFromRefLinear - Leaf Area Index from reflectance
• LAIFromNDVIFormo - Leaf Area Index from reflectance. Dotyczy tylko Formosata 2.

2. Woda:

• NDWI - Normalized Difference Water Index
• NDWI2 - Normalized Difference Water Index
• MNDWI - Modified Normalized Difference Water Index
• NDPI - Normalized Difference Pond Index
• NDTI - Normalized Difference Turbidity Index

3. Gleba:

• RI - Redness Index
• CI - Color Index
• BI - Brilliance Index
• BI2 - Brilliance Index

     W celu obliczenia interesującego nas indeksu w oknie wczytujemy obraz wielospektralny (jak wykonać taki obraz w QuantumGIS z poszczególnych kanałów spektralnych opisałem w poście Kompozycje barwne w QGIS) i zaznaczamy odpowiednią kolejność zakresów naszego zobrazowania np. dla Landsata 7 kanałem czerwonym będzie – 3, a kanałem podczerwonym – 4. Następnie z rozwijalnej listy wybieramy indeks wegetacyjny. Po określeniu miejsca zapisu pliku wynikowego i wygenerowaniu finalnego obrazu mamy do dyspozycji scenę z obliczonym indeksem.