Алгоритм обнаружения изменений с помощью экспоненциально взвешенного скользящего среднего. Этот алгоритм вычисляет гармоническую модель для «обучающей» части входных данных и вычитает её из исходных результатов. Затем остатки подвергаются построению диаграмм Шеухарта X-bar и построению экспоненциально взвешенного скользящего среднего. Пиксели с отклонениями указываются, когда диаграммы сигнализируют об отклонении от заданных контрольных пределов.
В результате получается пятиканальное изображение, содержащее следующие каналы:
- ewma: одномерный массив значений EWMA для каждого входного изображения. Отрицательные значения указывают на нарушение, а положительные — на восстановление.
- harmonicCoefficients: Одномерный массив пар вычисленных гармонических коэффициентов. Коэффициенты упорядочены в следующем порядке: [constant, sin0, cos0, sin1, cos1...]
- rmse: среднеквадратичная ошибка (RMSE) гармонической регрессии.
- rSquared: значение коэффициента детерминации (r-квадрат), полученное в результате гармонической регрессии.
- остатки: одномерный массив остатков гармонической регрессии.
См.: Brooks, EB, Wynne, RH, Thomas, VA, Blinn, CE и Coulston, JW, 2014. Обнаружение изменений в режиме реального времени с использованием статистических контрольных диаграмм качества и данных Landsat. IEEE Transactions on Geoscience and Remote Sensing, 52(6), pp.3316-3332.
| Использование | Возвраты | ee.Algorithms.TemporalSegmentation.Ewmacd(timeSeries, vegetationThreshold, trainingStartYear, trainingEndYear, harmonicCount , xBarLimit1 , xBarLimit2 , lambda , lambdasigs , rounding , persistence ) | Изображение |
| Аргумент | Тип | Подробности | timeSeries | Коллекция изображений | Коллекция, из которой будут извлечены изображения EWMA. Предполагается, что эта коллекция будет содержать по одному изображению за каждый год и будет отсортирована по времени. |
vegetationThreshold | Плавать | Пороговое значение для растительности. Значения ниже этого порога считаются отсутствием растительности. |
trainingStartYear | Целое число | Начало учебного года включительно. |
trainingEndYear | Целое число | Завершение года обучения (эксклюзивно). |
harmonicCount | Целое число, по умолчанию: 2 | Количество использованных пар гармонических функций (синус и косинус). |
xBarLimit1 | Число с плавающей запятой, значение по умолчанию: 1.5 | Пороговое значение для ограничения xBar при начальной тренировке. |
xBarLimit2 | Целое число, значение по умолчанию: 20 | Пороговое значение для достижения лимита xBar. |
lambda | Число с плавающей запятой, значение по умолчанию: 0.3 | Параметр настройки «лямбда» определяет вес нового года по сравнению со скользящим средним. |
lambdasigs | Число с плавающей запятой, значение по умолчанию: 3 | Границы контроля EWMA, выраженные в единицах стандартных отклонений. |
rounding | Логическое значение, по умолчанию: true | Следует ли округлять значения EWMA? |
persistence | Целое число, по умолчанию: 3 | Минимальное количество наблюдений, необходимое для рассмотрения возможности изменения. |