Доступно с лицензией Image Analyst.
Доступно с лицензией Spatial Analyst.
Сводка
Создает растровый объект, содержащий аномальные значения пикселов входного многомерного растра на основе временного интервала и метода вычисления аномалий.
Описание
Воспользуйтесь функцией Anomaly для вычисления аномальных значений пикселов для временных переменных многомерного растрового объекта. Например, если у вас есть ежедневные данные об объеме выпавших осадков, вы сможете найти пикселы, значения которых сильно отличаются от значений среднемесячного объема выпавших осадков по всему вашему набору данных.
Существует три математических метода, использующихся для вычисления аномальных значений:
- Отличие от среднего = x - µ
- x = значение пиксела в срезе
- µ = среднее значение этого пиксела за заданный временной интервал
- Процентное отличие от среднего = |x - µ| /[(x + µ)/2]
- x = значение пиксела в срезе
- µ = среднее значение этого пиксела за заданный временной интервал
- |x - µ| = абсолютное значение разницы между значением и средним значением
- Процент от усредненного = x / µ
- x = значение пиксела в срезе
- µ = среднее значение этого пиксела за заданный временной интервал
Указанный набор растровых данных является временным для растрового объекта. Чтобы сделать его постоянным, можно вызвать метод save растровго объекта.
Синтаксис
Anomaly (in_raster, {dimension_name}, {anomaly_type}, {temporal_interval}, {ignore_nodata})
Parameter | Объяснение | Тип данных |
in_raster | Входной многомерный набор растровых данных. | Raster |
dimension_name | Имя измерения, вдоль которого будет выполняться поиск аномалии. (Значение по умолчанию — StdTime) | String |
anomaly_type | Определяет метод, используемый для расчета аномалий.
(Значение по умолчанию — DIFFERENCE_FROM_MEAN) | String |
temporal_interval | Задает временной интервал, использующийся для вычисления среднего значения.
(Значение по умолчанию — None) | String |
ignore_nodata | Определяет, будут ли игнорироваться значения NoData в анализе.
(Значение по умолчанию — True) | Boolean |
Тип данных | Объяснение |
Raster | Выходной многомерный растр аномалий. |
Пример кода
Вычисляет аномалии объема выпавших осадков по всем срезам для среднемесячных значений.
import arcpy
from arcpy.ia import *
arcpy.CheckOutExtension("ImageAnalyst")
in_raster = Raster('ClimateData_Daily.nc', True)
# Choose the precipitation variable
prcp_raster = Subset(in_raster, variables = 'prcp')
# Calculate percent of mean anomalies,
# where the mean_value is the mean of all slices
prcp_anomaly = Anomaly(prcp_raster, anomaly_type = 'PERCENT_OF_MEAN')
# Calculate the difference from mean,
# where the mean_value is the mean of each year
prcp_monthly_anomaly = Anomaly(
prcp_raster, temporal_interval = 'YEARLY')