生成趋势栅格 (影像分析工具)

汇总

用于面向多维栅格中一个或多个变量估计每个像素沿维度的趋势。

用法

该工具的一个应用示例是：假设您拥有 40 年的月度海洋温度数据，并希望为每个像素拟合一条趋势线，以观察温度随时间的变化趋势及其具体变化情况。
受支持的多维栅格数据集包括云栅格格式 (CRF)、多维镶嵌数据集或者由 netCDF、GRIB 或 HDF 格式文件生成的多维栅格图层。
该工具将生成云栅格格式 (CRF) 的多维栅格数据集。目前，尚不支持其他输出格式。
默认情况下，将使用 LZ77 压缩类型对多维栅格输出进行压缩。但是，建议将压缩类型更改为 LERC 并根据数据来调整最大误差值。例如，如果您希望分析结果精确到三位小数，对于最大误差值，请使用 0.001。避免不必要的精度要求，因为它们将增加处理时间和存储大小。
该工具可用于沿线性、谐波或多项式趋势线拟合数据，也可用于使用 Mann-Kendall 或 Seasonal-Kendall 测试执行趋势检测。
Mann-Kendall 和 Seasonal-Kendall 测试可用于确定数据中是否存在单调性趋势。上述测试为非参数测试，这意味着它们不会假设具体的数据分布。Mann-Kendall 测试不考虑序列相关性或季节性影响。如果数据为季节性数据，则 Seasonal-Kendall 测试更为适合。
如果使用该工具执行 Mann-Kendall 或 Seasonal-Kendall 测试，则输出将为五波段栅格，如下所示：
- 波段 1 = Sen 斜率
- 波段 2 = p 值
- 波段 3 = Mann-Kendall 得分 (S)
- 波段 4 = S 方差
- 波段 5 = Z 得分
Mann-Kendall 或 Seasonal-Kendall 测试的输出可用于确定多维时间序列中的哪些像素具有统计显著性趋势。您可以将此信息与线性、谐波或多项式趋势分析结合使用，以提取时间序列中的重要趋势。您可以生成包含像素（像素带有显著 p 值）的掩膜，将该掩膜应用于多维栅格，以及将此掩膜多维栅格用作工具的输入，以执行线性、谐波或多项式趋势分析。
系统提供三个趋势线选项用于沿维度拟合变量值的趋势：线性、谐波和多项式。三个趋势拟合选项如下所述。
- Linear—线性趋势线是最佳拟合直线，可用于估计简单的线性关系。线性趋势可突出显示以稳定速率增加或减少的变化率。线性趋势线的公式如下：
  - y = 像素变量值
  - x = 维度值
  - ß0 = y 截距
  - ß1 = 线性斜率或变化率
    
    ß1 > 0 表示增加趋势
    
    ß1 < 0 表示减少趋势
- Harmonic—谐波趋势线是周期性重复的曲线，最适用于描述遵循周期性模式的数据，例如季节性温度变化。谐波趋势线的公式如下：
  - y = 像素变量值
  - t = 儒略日
  - ß0 = y 截距
  - ß1 = 变化率
  - α, γ = 年际或年内变化的系数
  - ω = i
  - f = 谐波频率
- Polynomial—多项式趋势线是一条曲线，可用于波动的数据。在这种情况下，多项式的阶数值用于指示发生的最大波动数。多项式趋势线的公式如下：
  - y = 像素变量值
  - x = 维度值
  - ß0, ß1, ß2, ß3, ..., ßn = 常系数
如果该工具用于线性、谐波或多项式趋势分析，输出趋势栅格可以作为使用趋势栅格预测工具的输入。趋势栅格属于多维栅格，其中每个剖切均是一个包含趋势线相关信息的多波段栅格。如果要分析包含单个维度（例如，时间）的数据集的单个变量趋势，则输出数据集中将具有单个剖切。如果要分析包含多个维度（例如，时间和深度）的数据集的单个变量，则每个剖切将包含沿分析中未包含维度的每个维度值的趋势信息。

对于线性趋势分析，输出包含三波段栅格，其中：
- 波段 1 = 斜率
- 波段 2 = 截距
- 波段 3 = 均方根误差 (RMSE) 或最佳拟合线周围的误差
对于谐波趋势分析，输出中的波段数取决于谐波频率。当频率设置为 1 时，输出将为五波段栅格，如下所示：
- 波段 1 = 斜率
- 波段 2 = 截距
- 波段 3 = Harmonic_sin1
- 波段 4 = Harmonic_cos1
- 波段 5 = RMSE
当频率设置为 2 时，输出将为七波段栅格，如下所示
- 波段 1 = 斜率
- 波段 2 = 截距
- 波段 3 = Harmonic_sin1
- 波段 4 = Harmonic_cos1
- 波段 5 = Harmonic_sin2
- 波段 6 = Harmonic_cos2
- 波段 7 = RMSE
对于多项式趋势分析，输出中的波段数取决于多项式阶数。二阶多项式拟合将生成四波段栅格，如下所示：
- 波段 1 = Polynomial_2
- 波段 2 = Polynomial_1
- 波段 3 = Polynomial_0
- 波段 4 = RMSE
三阶多项式拟合将生成五波段栅格，如下所示：
- 波段 1 = Polynomial_3
- 波段 2 = Polynomial_2
- 波段 3 = Polynomial_1
- 波段 4 = Polynomial_0
- 波段 5 = RMSE
谐波趋势分析的周期长度参数用于表示您希望在一天或全年的数据中看到的周期数量和周期长度。例如，如果您希望数据在一年中经历两个变化周期，则周期长度将为 182.5 天或 0.5 年。如果您每三小时采集一次温度数据，并且每天存在一个变化周期，则周期长度将为 1 天。
谐波趋势分析的频率参数用于描述谐波模型以拟合至数据。如果频率设置为 1，则将使用线性和一阶谐波曲线的组合来拟合模型。如果频率为 2，则将使用线性、一阶谐波和二阶谐波曲线的组合来拟合数据。如果频率为 3，则将使用附加的三阶谐波曲线对数据进行建模，依此类推。
可以生成模型拟合优度统计数据作为线性、谐波或多项式趋势拟合的可选输出。可以计算均方根误差 (RMSE)、R 平方和趋势斜率 p 值并显示在输出栅格统计数据部分的属性窗口中。也可以使用 RGB 符号表示输出趋势栅格，并将统计数据指定为红色、绿色、蓝色波段来显示统计数据。

标注	说明	数据类型
输入多维栅格	输入多维栅格数据集。	Raster Dataset; Raster Layer; Mosaic Dataset; Mosaic Layer; Image Service; File
维度	将沿此维度为分析中所选的一个或多个变量提取趋势。	String
变量 [维度信息]（描述）（可选）	将计算趋势的一个或多个变量。如果未指定变量，则将分析多维栅格中的第一个变量。	String
趋势类型（可选）	指定对沿某一维度的像素值执行的趋势分析类型。线性—将沿线性趋势线拟合变量像素值。这是默认设置。多项式—将沿二阶多项式趋势线拟合变量像素值。谐波—将沿谐波趋势线拟合变量像素值。 Mann-Kendall—变量像素值将使用 Mann-Kendall 趋势测试进行评估。 Seasonal-Kendall—变量像素值将使用 Seasonal-Kendall 趋势测试进行评估。	String
频率 / 多项式阶数（可选）	趋势拟合中要使用的频率或多项式阶数。如果趋势类型为多项式，此参数指定多项式阶数。如果趋势类型为谐波，此参数指定用于拟合趋势的模型数量。仅当所分析的维度为时间时，该参数才会包含在趋势分析中。如果趋势类型参数设置为谐波，默认值为 1，意味着使用一阶谐波曲线来拟合模型。如果趋势类型参数设置为多项式，默认值为 2，即二阶多项式。	Long
忽略 NoData （可选）	指定分析中是否忽略 NoData 值。选中—分析将包括沿给定维度的所有有效像素，并忽略 NoData 像素。这是默认设置。未选中—如果沿给定维度的像素包含 NoData 值，则分析结果将变为 NoData。	Boolean
周期长度（可选）	要进行建模的周期性变化的长度。当趋势类型设置为谐波时，需要使用此参数。例如，叶绿度通常在一年中具有一个较强的变化周期，因此周期长度为 1 年。每小时温度数据在一天中具有一个较强的变化周期，因此周期长度为 1 天。对于每年周期变化的数据，默认长度为 1 年。	Double
周期单位（可选）	指定用于谐波周期长度的时间单位。天—谐波周期的长度单位为天。年—谐波周期的长度单位为年。这是默认设置。	String
RMSE （可选）	指定是否将计算趋势拟合线的均方根误差 (RMSE)。选中—将计算 RMSE。这是默认设置。未选中—不计算 RMSE。	Boolean
R 平方（可选）	指定是否将为趋势拟合线计算 R 平方拟合优度统计数据。选中—将计算 R 平方值。未选中—不计算 R 平方值。这是默认设置。	Boolean
斜率系数的 P 值（可选）	指定是否将计算趋势线斜率系数的 p 值统计数据。选中—将计算 p 值。未选中—不计算 p 值。这是默认设置。	Boolean
季节性期间（可选）	指定在执行 Seasonal-Kendall 测试时用于季节性期间长度的时间单位。天—季节性期间的长度单位为天。这是默认设置。月—季节性期间的长度单位为月。	String

返回值

标注	说明	数据类型
输出多维栅格	输出云栅格格式 (CRF) 多维栅格数据集。	Raster

GenerateTrendRaster(in_multidimensional_raster, dimension, {variables}, {line_type}, {frequency}, {ignore_nodata}, {cycle_length}, {cycle_unit}, {rmse}, {r2}, {slope_p_value}, {seasonal_period})

名称	说明	数据类型
in_multidimensional_raster	输入多维栅格数据集。	Raster Dataset; Raster Layer; Mosaic Dataset; Mosaic Layer; Image Service; File
dimension	将沿此维度为分析中所选的一个或多个变量提取趋势。	String
variables [variables,...] （可选）	将计算趋势的一个或多个变量。如果未指定变量，则将分析多维栅格中的第一个变量。	String
line_type （可选）	指定对沿某一维度的像素值执行的趋势分析类型。 `LINEAR`—将沿线性趋势线拟合变量像素值。这是默认设置。 `POLYNOMIAL`—将沿二阶多项式趋势线拟合变量像素值。 `HARMONIC`—将沿谐波趋势线拟合变量像素值。 `MANN-KENDALL`—变量像素值将使用 Mann-Kendall 趋势测试进行评估。 `SEASONAL-KENDALL`—变量像素值将使用 Seasonal-Kendall 趋势测试进行评估。	String
frequency （可选）	趋势拟合中要使用的频率或多项式阶数。如果趋势类型为多项式，此参数指定多项式阶数。如果趋势类型为谐波，此参数指定用于拟合趋势的模型数量。仅当所分析的维度为时间时，该参数才会包含在趋势分析中。如果 `line_type` 参数为 `HARMONIC`，默认值为 1，意味着使用一阶谐波曲线来拟合模型。如果 `line_type` 参数为 `POLYNOMIAL`，默认值为 2，即二阶多项式。	Long
ignore_nodata （可选）	指定分析中是否忽略 NoData 值。 `DATA`—分析将包括沿给定维度的所有有效像素，并忽略 NoData 像素。这是默认设置。 `NODATA`—如果沿给定维度的像素包含 NoData 值，则分析结果将变为 NoData。	Boolean
cycle_length （可选）	要进行建模的周期性变化的长度。当 `line_type` 设置为 `HARMONIC`时，需要使用此参数。例如，叶绿度通常在一年中具有一个较强的变化周期，因此周期长度为 1 年。每小时温度数据在一天中具有一个较强的变化周期，因此周期长度为 1 天。对于每年周期变化的数据，默认长度为 1 年。	Double
cycle_unit （可选）	指定用于谐波周期长度的时间单位。 `DAYS`—谐波周期的长度单位为天。 `YEARS`—谐波周期的长度单位为年。这是默认设置。	String
rmse （可选）	指定是否将计算趋势拟合线的均方根误差 (RMSE)。 `RMSE`—将计算 RMSE。这是默认设置。 `NO_RMSE`—不计算 RMSE。	Boolean
r2 （可选）	指定是否将为趋势拟合线计算 R 平方拟合优度统计数据。 `R2`—将计算 R 平方值。 `NO_R2`—不计算 R 平方值。这是默认设置。	Boolean
slope_p_value （可选）	指定是否将计算趋势线斜率系数的 p 值统计数据。 `SLOPEPVALUE`—将计算 p 值。 `NO_SLOPEPVALUE`—不计算 p 值。这是默认设置。	Boolean
seasonal_period （可选）	指定在执行 Seasonal-Kendall 测试时用于季节性期间长度的时间单位。 `DAYS`—季节性期间的长度单位为天。这是默认设置。 `MONTHS`—季节性期间的长度单位为月。	String

返回值

名称	说明	数据类型
out_multidimensional_raster	输出云栅格格式 (CRF) 多维栅格数据集。	Raster

代码示例

GenerateTrendRaster 示例 1（Python窗口）

本示例计算降水量数据随时间变化的多项式趋势拟合。

# Import system modules
import arcpy
from arcpy.ia import *

# Run
trendCoeffMultidim = arcpy.ia.GenerateTrendRaster("model_time_series.crf",
    "StdTime", "precip", "POLYNOMIAL", 2, "DATA", None, '', "RMSE","R2", None)

# Save output
trendCoeffMultidim.save("C:/data/polynomial_trend_coefficients.crf")

GenerateTrendRaster 示例 2（独立脚本）

本示例计算温度和盐度数据随深度变化的线性趋势拟合。

# Import system modules
import arcpy
from arcpy.ia import *

# Define input parameters
inFile = "C:/Data/ClimateVariables.crf"
dimensionName = "Depth"
variables = "salinity;temp"
trend_model = "LINEAR"
ignore_nodata = "DATA"

# Run - fit a linear model for salinity and temperature
trendCoeffMultidim = GenerateTrendRaster(inFile, dimensionName,
    variables, trend_model, None, ignore_nodata)

# Save output
trendCoeffMultidim.save("C:/data/linear_trend_coefficients.crf")

GenerateTrendRaster 示例 3（独立脚本）

本示例测试单调性趋势，然后计算沿 NDVI 时间序列的谐波趋势拟合。

# Import system modules
import arcpy
from arcpy.ia import *

# Define input parameters

inFile = "C:/data/ndvi_time_series.crf"
dimensionName = "StdTime"
variables = "ndvi"
trend_test = "SEASONAL_KENDALL"
trend_model = "HARMONIC"
harmonic_frequency = 1
ignore_nodata = "DATA"
length_of_cycle = 365.25
rmse = RMSE
r_square = NO_r2
p_value = NO_SLOPEPVALUE

# Test for monotonic trend

trendtest = GenerateTrendRaster(inFile, dimensionName, variables,
    trend_test, '', ignore_Nodata, '', '', '', '', '', 'MONTHS')

# Run
trendCoeffMultidim = GenerateTrendRaster(inFile, dimensionName,
    variables, trend_model, harmonic_frequency, ignore_nodata,
    length_of_cycle,rmse,r_square,p_value)

# Save output
trendCoeffMultidim.save("C:/data/harmonic_trend_coefficients.crf")

环境

像元大小, 压缩, 当前工作空间, 地理变换, NoData, 输出配置关键字, 输出坐标系, 范围, 并行处理因子, 金字塔, 栅格统计, 重采样方法, 临时工作空间, 捕捉栅格, 切片大小

许可信息

基本: 需要 Image Analyst
标准: 需要 Image Analyst
高级: 需要 Image Analyst

生成趋势栅格 (影像分析工具)

汇总

用法

参数

返回值

环境

许可信息

相关主题