Спутник ДЗЗ

Код, который все это рассчитывает в Google Earth Engine, будет сегодня несколько позже.

174 views05:38

GEE-26. ForWarn: мониторинг состояния леса

Продолжение разговора, начатого здесь и здесь.

Код: https://code.earthengine.google.com/9f473f5bde9b236eda92cbe28bbfd243

Скрипт рассчитывает два вида оценок изменения состояния леса:

1. 1/3/5/10yr Baseline — для выявления нарушений.
2. All-Year Baseline — для сравнения текущего состояния с климатической нормой.

Маска леса строится по данным ESA WorldCover, по состоянию на январь 2020 года.

Для определения причин изменения состояния леса создаются слои данных:

* карта очагов возгораний (данные FIRMS) — слой Fires
* месячные медианные RGB-композиты Sentinel-2 (Natural Colors) — слой Month Median Composite.

В скрипте использованы идеи, которые уже применялись нами для расчета Vegetation Condition Index. Из новинок в коде: оператор выбора switch — такой же, как в других языках программирования, и условие where, примененное при создании маски леса.

Настройка скрипта выполняется в разделе “Пользовательские параметры”:


// AOI.
var point = ee.Geometry.Point(-80.2265, 34.3676);
var width = 20000; // ширина буфера вокруг точки, м
var geometry = point.buffer({'distance': width}).bounds();

// Дата наблюдения.
var date = ee.Date('2019-09-01');
// Базовый период.
var base = '1'; // Значения: 1/3/5/10 - для max; ALL - для mean.

geometry — геометрия AOI (один полигон). Нужно загрузить в Assets шейпфайл AOI (CRS: WGS 84), нарисовать AOI прямо в GEE или задать точку, как в примере. Базовый период base выбирается в зависимости от задачи: для поиска нарушений или для сравнения с климатической нормой. base = '1' — это 1yr Baseline, аналогично задаются другие интервалы. base = 'ALL' — это сравнение с климатической нормой (за весь период наблюдений MODIS).

Маскирование леса можно отменить, закомментировав строку percent_change = ... в разделе "Расчет":


// Маскирование леса.
var forest = getForestMask();
percent_change = percent_change.multiply(forest);

Учебник по Google Earth Engine

#лес #GEE

209 viewsedited 11:54