Алгоритм поиска временных интервалов геохимических аномалий на основе продуктов TERRA/MODIS для Нижне-Нойбинский рудного узла

Глобальной задачей всего исследования является поиск корреляционных связей между геохимическими аномалиями и спектральными характеристиками объектов, которые располагаются на участках с повышенной концентрацией химических элементов в почве. В рамках данной работы была поставлена цель, выделить временные интервалы геохимических аномалий. В качестве пространственных данных для исследования использовались продукты MODIS, а именно EVI, NDVI, GPP и LAI.
В основе исследования лежит анализ данных и построение временных рядов по данным полученным за последние 13 лет. Для исследования был выбран период с мая по сентябрь включительно для анализа всего вегетационного цикла. Были определен перспективный участок Нойбинского месторождения - Нижне-Нойбинский рудный узел.
Было принято решение анализировать данные растительного покрова по 5 месяцам (май, июнь, июль, август, сентябрь), для того чтобы охватить весь вегетационный цикл растительности.
Для визуализации и дальнейшей обработки полученных значений использовался пакет прикладных программ для решения задач технических вычислений MatLab. При поиске аномальных значений данных по каждому из продуктов MODIS, за отличительный признак были приняты скачкообразные изменения статистических данных. «Скорость» этих изломов значительно превосходит относительно плавные сезонные изменения показателей растительности. Для выделения резко изламывающихся сигналов было применено двойное дифференцирование. В результате этой операции был получен график отображающий крутизну “изломов”.
Вторым этапом введение порогового детектора, т.е. значения “крутизна” которых не превышала естественного сезонного изменения приравнивались к нулю. Все значимые отклонения от нормы фиксировались на графике.
Завершающим этапом является выделение тех временных интервалов которые содержат одновременно аномальные значения по 4 показателям характеризующих растительность. Таким образом, было выделено 8 временных интервалов.
В дальнейшем, планируется применять данные со спутников среднего пространственного разрешения, съемка которых проводилась в указанные промежутки времени.