检测像素值变化

使用连续数据的变化检测基于逐个像素地计算光学影像或建模数据(例如温度或风速)之间的差异。

比较连续建模数据的目的通常是识别在一段时间内幅度或特定方向发生变化的区域。 ArcGIS AllSource 可提供多种用于比较多个连续栅格数据集的方法。

两个植被指数栅格之间的像素值变化检测

计算变化栅格函数可即时计算两个栅格图层之间的差异。 仅处理可见范围内的像素,如果缩小到超出输入数据的源分辨率,则像素将重新采样为更大尺寸以加快处理速度。 这是一种计算变化的方法,用于快速可视化和评估变化。 对于更复杂的工作流,可以将此函数与其他栅格函数结合使用。

计算变化栅格地理处理工具计算两个栅格数据集之间的差异,并生成一个包含变化信息的新栅格数据集。 对于更复杂的工作流,可以将此工具与其他工具结合使用。

变化检测向导将工具与函数相结合,引导您完成类别变化检测、像素值变化检测、光谱变化检测或时间序列变化检测的执行过程。 您可以比较两个栅格数据集或多维栅格数据集中的两个时间片。

以下部分将详细介绍执行连续变化检测时变化检测向导中的每个窗格。

变化检测向导

可从分析选项卡栅格组中的变化检测下拉按钮启动“变化检测”向导。 如果您不在 2D 地图场景中使用,该按钮将不可用。

配置窗格

变化检测向导中的第一个窗格是配置窗格,您可以在其中设置变化检测方法值。 要比较连续数据集,请将变化检测方法设置为像素值变化

在启动向导后,如果在内容窗格中选择连续栅格,变化检测方法选项将默认设置为像素值变化

参数描述

由栅格转出

计算中要使用的第一个栅格。 要评估从时间 1(较早)到时间 2(较晚)的变化,请在此处输入时间 1 栅格。

转为栅格

计算中要使用的第二个栅格。 要评估从时间 1(较早)到时间 2(较晚)的变化,请在此处输入时间 2 栅格。

要比较多维数据集中的两个时间片,请输入多维栅格作为起始栅格参数。 配置窗格将使用下表中的参数进行更新。

参数描述

变量

要分析的变量。

尺寸

用于比较时间片的维度字段。

起始时间片

计算中要使用的第一个时间片。 要评估从时间 1(较早)到时间 2(较晚)的变化,请在此处输入时间 1 时间片。

终止时间片

计算中要使用的第二个时间片。 要评估从时间 1(较早)到时间 2(较晚)的变化,请在此处输入时间 2 时间片。

处理范围

输出的处理范围。

X 跳跃因子

x 跳跃因子,用于计算统计数据和变化栅格的直方图。

Y 跳跃因子

y 跳跃因子,用于计算统计数据和变化栅格的直方图。

波段差异窗格

波段差异窗格可用于指定比较输入连续栅格数据集的方法。 如果数据为多波段影像,可以比较两个影像的单个波段,或首先计算波段指数,然后再执行比较操作。

参数描述

差异类型

计算中使用的差异方法。

  • 绝对 - 起始栅格参数中的像素值和目标栅格参数中的像素值之间的数学差或两者做减法。
    输出 = (目标栅格) - (起始栅格)
    这是默认设置。
  • 相对 - 像素值的差异(考虑要比较的值的数量)。
    输出 = (目标栅格 - 起始栅格) / max (目标栅格, 起始栅格)

波段差异方法

选择要比较连续栅格的方法。

  • 单个波段差异 - 计算每个栅格中单个波段的波段值之间的差异。
  • 波段指数差异 - 计算每个栅格的波段指数,并计算指数值之间的差异。

起始栅格波段目标栅格波段

起始栅格参数和目标栅格参数中分别选择要在计算中比较的波段。 当波段差异方法设置为单个波段差异时,此参数为必需项。

通常,您将比较两个影像中同一波段的像素值。 例如,比较两个日期的近红外波段值。 有时,两个影像的波段指数不同,例如,如果比较 Landsat 5 影像和 Landsat 8 影像,因此指定每个影像的波段至关重要。

波段索引

选择要为每个栅格计算的波段指数,然后计算其间差异。 当波段差异方法设置为单个波段差异时,此参数为必需项。

选项列表与指数库中的指数列表相同。

您必须指定为每个影像生成指数所需的波段。

像元大小类型

选择输出栅格中使用的像元大小。 如果所有输入像元大小均相同,则所有选项都会产生相同的结果。

  • 第一个 - 使用输入栅格的第一个像元大小。
  • 最小值 - 使用所有输入栅格中的最小像元大小。
  • 最大值 - 使用所有输入栅格中的最大像元大小。 这是默认设置。
  • 平均值 - 使用所有输入栅格的平均像元大小。
  • 最后一个 - 使用输入栅格的最后一个像元大小。

范围类型

选择输出栅格中应使用的范围。

  • 第一个 - 使用第一个输入栅格的范围来确定处理范围。

  • 交集 - 使用叠置像素的范围来确定处理范围。 这是默认设置。
  • 并集 - 使用所有栅格的范围来确定处理范围。
  • 最后一个 - 使用最后一个输入栅格的范围来确定处理范围。

单击预览,使用您指定的选项将预览图层添加到地图。 您可以修改选项,然后单击更新预览以使用新选项刷新预览图层。

注:

预览图层使用栅格函数生成。 缩小到超出源分辨率时,栅格函数将使用重采样像素大小来处理数据。 为确保预览在最终结果以您希望的方式显示,请缩放到数据的源分辨率。 右键单击预览图层,然后单击缩放至源分辨率

当您在波段差异窗格中单击下一步时,向导将计算波段指数(如果适用)、计算栅格之间的差异并计算差异值的直方图。 (此操作可能需要一些时间,具体取决于数据集大小。)

“差异分类”窗格

借助差异分类窗格,您可以实时探索地图中的像素值差异,并使用该信息对结果进行分类。 例如,您可以将 NDVI 中的变化分类为亏损或盈利类,也可以对 NIR 波段值差异进行分类以提取火灾后迹地,并将其分类为燃烧面积类。

探索差异部分占据差异分类窗格的上半部分。

探索差异部分和交互式直方图
显示差异值直方图以及用于可视化变化的交互式控点。 打开此窗格后,默认情况下,系统会将 Preview_Mask 图层添加到地图中。 当您将最小值和最大值直方图控点拖动到直方图中的不同值处时,最小值和最大值之间的像素值会显示在地图上的 Preview_Mask 图层中。 所有其他像素值均会被掩膜且不可见。 如果地图中的所有其他图层均已关闭,可能更容易查看这些值。 您还可以修改 Preview_Mask 图层的符号系统。

要跳过分类并保留原始像素值差异,请取消选中对值差异进行分类参数,然后单击下一步转到下一个窗格。

使用交互式直方图生成类

要使用交互式直方图生成类,请完成以下步骤:

  1. 拖动直方图中的最小值和最大值控点,直到确定分析所需的重要变化信息(例如,强烈负变化)。
  2. 单击添加新类按钮 插入下拉箭头

    系统会将在直方图中选择的最小值和最大值作为最小和最大类值添加到窗格底部的输出分类表中。

  3. 分类输出表中,将类值(整数)添加到 Output 字段中,并将类名称添加到 Class Name 字段中。

    这些项是必需项。

  4. 指定用于符号化类的颜色。

手动生成类

要手动输入类信息,而无需使用交互式直方图,请完成以下步骤:

  1. 单击生成按钮 配色方案
  2. 生成窗格中,将最大值参数设置为要生成的类数量。

    可以为与每个唯一值相关联的类指定基本名称。 您还可以为类符号系统指定配色方案和透明度。

  3. 完成后,单击确定
  4. 分类输出表中,为每个类手动输入 MinimumMaximum 差异字段值。

单击预览,使用您指定的选项将预览图层添加到地图。 您可以修改选项,然后单击更新预览以使用新选项刷新预览图层。

注:

预览图层使用栅格函数生成。 缩小到超出源分辨率时,栅格函数将使用重采样像素大小来处理数据。 为确保预览在最终结果以您希望的方式显示,请缩放到数据的源分辨率。 右键单击预览图层,然后单击缩放至源分辨率

后处理窗格

可通过后处理窗格平滑结果并生成输出。

参数描述

平滑邻域

用于平滑结果的焦点邻域的大小(以像素行和列为单位)。

  • - 不执行平滑。 这是默认设置。
  • 3 × 3 - 基于 3 像素行乘 3 像素列的邻域平滑结果。
  • 5 × 5 - 基于 5 像素行乘 5 像素列的邻域平滑结果。
  • 7 × 7 - 基于 7 像素行乘 7 像素列的邻域平滑结果。

统计数据填充方法

用于重新计算像素值以获得更平滑结果的统计方法。

  • 最小值 - 计算邻域内像素的最小值。
  • 最大值 - 计算邻域内像素的最大值。
  • 均值 - 计算邻域内像素的平均值。 这是默认设置。
  • 标准差 - 计算邻域内像素的标准差值。
  • 中值 - 计算邻域内像素的中值。
  • 众数 - 计算邻域内像素的众数值或出现次数最多的值。
  • 少数 - 计算邻域内像素的少数值或出现次数最少的值。

结果另存为

指定要生成的输出类型。

  • 栅格数据集 - 将结果另存为栅格数据集。
  • 栅格函数模板 - 将结果另存为栅格函数模板,以用于其他处理。
  • 要素类 - 将结果另存为面要素类。

输出数据集

输出数据集的名称。 此参数仅在结果另存为被设置为栅格数据集要素类时可用。 如果输出为栅格数据集,请指定文件扩展名以生成受支持格式的栅格。

对于其他平滑选项,您可以对结果使用统计栅格函数。 在向导中选择作为平滑邻域,并对输出栅格数据集使用“统计”栅格函数。

单击预览,使用您指定的选项将预览图层添加到地图。 您可以修改选项,然后单击更新预览以使用新选项刷新预览图层。

注:

预览图层使用栅格函数生成。 缩小到超出源分辨率时,栅格函数将使用重采样像素大小来处理数据。 为确保预览在最终结果以您希望的方式显示,请缩放到数据的源分辨率。 右键单击预览图层,然后单击缩放至源分辨率

提取土壤调整植被指数 (SAVI) 值的损失

以下示例将提取两个 Landsat 8 影像的土壤调整植被指数 (SAVI) 值的损失。

  1. 向地图添加两个 Landsat 8 图像。
  2. 内容窗格中选择较早的影像图层后,从影像选项卡的分析组中启动变化检测向导。

    配置窗格中,将起始栅格参数设置为较早的栅格图层。

  3. 目标栅格参数设置为较新的栅格图层。
  4. 单击下一步
  5. 波段差异窗格中,设置用于计算 SAVI 值变化的分析。
    1. 差异类型参数设置为绝对
    2. 波段差异方法参数设置为波段指数差异
    3. 波段指数参数 - 起始栅格部分中,将近红外波段指数设置为波段 5,并将红色波段指数设置为波段 4。 保留土壤亮度校正系数的默认值 0.5
    4. 波段指数参数 - 目标栅格部分中,将近红外波段指数设置为波段 5,并将红色波段指数设置为波段 4。 保留土壤亮度校正系数的默认值 0.5
    5. 像元大小类型参数设置为最大,并将范围类型参数设置为交点
  6. 单击下一步
  7. 差异分类窗格中,将最大值直方图控点拖动到 0,并将最小值直方图控点保持在最小负值位置。 单击添加新类
  8. 分类输出表中,将 Output 字段值设置为 1,并将 Class Name 字段设置为 Loss,然后使用颜色选取器选择红色。
  9. 单击下一步
  10. 后处理窗格中,将平滑邻域参数设置为 3 x 3,并将统计数据填充方法参数设置为众数
  11. 对于输出数据集参数,输入 SAVI_Loss.tif。 单击运行
  12. 将数据集添加到地图后,单击完成关闭变化检测向导。

在本主题中
  1. 变化检测向导