ння NDVI для них.
Тип об'ектаОтраженіе в червоній області спектраОтраженіе в інфрачервоній області спектраЗначеніе NDVIГустая растітельность0, 10,50,7 Розріджена растітельность0, 10,30,5 Відкрита почва0, 250,30,025 Облака0, 250,250 Сніг і лед0 ,3750,35-0 , 05Вода0 ,020,01-0, 25Іскусственние матеріали (бетон, асфальт) 0,30,1-0,5
У наявних наборах даних вже були надані розраховані NDVI зображення. Тим не менш, було прийнято рішення отримати їх самостійно - причиною тому послужило те, що наявні NDVI зображення були нормалізовані, і механізм цієї нормалізації не був достатньо докладно описаний в специфікації знімків.
Для розрахунку NDVI використовувалася наведена вище формула і відповідні спектральні канали знімків. Після розрахунку для кожного з отриманих зображень було проведено розбиття на класи у відповідності зі специфікою досліджуваного об'єкта. Всього вийшло 11 класів із значеннями від 0 до 10 в цілих числах. Клас із значенням 0 призначений для об'єктів з NDV від - 1 до 0 включно, тобто для хмар, снігу, льоду, води і техногенних об'єктів. Далі класи йдуть з кроком 0,1 значення NDVI. Таким чином, в клас 1 потрапляє відкрита грунт, і далі класи від 2 і до 10 представляють рослинність в порядку підвищення її щільності. Після класифікації до цих зображень може бути застосована єдина палітра розмальовки [1], що дозволяє візуально порівняти отримані зображення. Слід також згадати про те, що на зображення були накладені маски хмарності і водних об'єктів, що додаються до вихідних даних. Значення NDVI в областях, що закриваються цими масками, було прийнято за 0. Також для знімків 2009 року було застосовано маски для областей, що не мають даних. Подібний підхід дозволяє точніше відсікти не цікавлять нас області.
На практиці описані вище процедури були виконані засобами декількох програмних пакетів. Спочатку за допомогою програми Global Mapper 12.0 з HDR - файлів знімків були витягнуті шари 3 і 4 каналів і збережені потім в окремі TIF - файли. Аналогічним чином з вихідних файлів були витягнуті маски хмарності і водних об'єктів. Паралельно для кожного з цих шарів були прописані вірні параметри проекції. Маски перепусток витягувати не виникло необхідності, так як вони вже були накладені на всі необхідні шари.
Потім в програмному пакеті ERDAS Imagine 9.0 були складені дві моделі обробки растрових даних. Перша модель на підставі 3 і 4 каналів знімка проводила розрахунок NDVI, маючи в якості граничного умови тільки нерівність 0 знаменника у формулі для його розрахунку. Друга модель класифікувала розрахований індекс за описаною вище схемою, паралельно накладаючи маски хмарності і водних об'єктів. Для результуючих растрів засобами цього ж програмного пакета відповідних класів присвоювалися відповідні кольори.
Нарешті, після обробки результуючі зображення поміщалися в ArcMap 10.1, де відбувалося їх остаточно оформлення, а саме вибір масштабу і відображеному на карті області, нанесення шкали і координатної сітки, а також складання і нанесення на карту легенди.
Як приклад нижче наведена пара таких зображень, отриманих для знімка 2002 року.
Малюнок 1 - Фрагмент NDVI зображення для знімка 2002 з застосованою атмосферної корекцією
Вже на підставі малюнків 1 і 2 можна впевнено сказати, що знімк...
