Python ArcGIS analyze_changes_using_ccdc用法及代码示例

本文简要介绍 python 语言中 arcgis.raster.analytics.analyze_changes_using_ccdc 的用法。

用法:

arcgis.raster.analytics.analyze_changes_using_ccdc(input_multidimensional_raster=None, bands_for_detecting_change=[], bands_for_temporal_masking=[], chi_squared_threshold=0.99, min_anomaly_observations=6, update_frequency=1, output_name=None, context=None, *, gis=None, future=False, **kwargs)

返回：

影像图层项目

函数使用 CCDC 算法评估像素值随时间的变化，并生成包含模型结果的多维栅格。 ArcGIS Image Server 10.8.1 及更高版本中可用的函数。

Parameter	Description
input_multidimensional_raster	必需的 ImageryLayer 对象。输入多维栅格。传送门物品可以通过。输入多维栅格必须至少有 12 个切片，跨度至少为 1 年。
bands_for_detecting_change	可选列表。用于更改检测的波段 ID。如果未提供波段 ID，则将使用输入栅格数据集中的所有波段。列表中的每个元素都应在 1 到 n 的范围内，其中 n 是输入栅格的波段数。 例子： [1,2,3,4,6]
bands_for_temporal_masking	可选列表。绿色波段和 SWIR 波段的波段 ID，用于遮蔽云、云影和雪。如果未提供波段 ID，则不会发生屏蔽。列表中的每个元素都应在 1 到 n 的范围内，其中 n 是输入栅格的波段数。 例子： [1,2]
chi_squared_threshold	可选浮点数。卡方变化概率阈值。如果观测值的计算变化概率高于此阈值，则将其标记为异常，这是潜在的变化事件。默认值为 0.99。 例子： 0.99
min_anomaly_observations	可选整数。在将事件视为更改之前必须发生的最小连续异常观察次数。在将像素视为真正的变化之前，必须将像素标记为指定数量的连续时间片的异常。默认值为 6。
update_frequency	可选浮点数。表示更新频率的值。默认值为 1。
output_name	可选的。如果未提供，则由该方法创建影像服务并将其用作输出栅格。您可以从 GIS 中传入现有的影像服务项目来使用它。或者，您可以传入应通过此方法创建的输出图像服务的名称，以用作工具的输出。如果该名称的服务已经存在，则会引发 RuntimeError
context	上下文包含影响任务执行的其他设置。 context 参数覆盖通过arcgis.env 参数设置的值 此函数具有以下设置： Extent(范围)：定义分析区域的边界框。 例子： {“extent”: {“xmin”: -122.68, “ymin”: 45.53, “xmax”: -122.45, “ymax”: 45.6, “spatialReference”: {“wkid”: 4326}}} 输出空间参考 (outSR)：输出栅格将投影到输出空间参考中。 例子： {“outSR”: {spatial reference}} 捕捉栅格 (snapRaster)：输出栅格的像元将与指定的捕捉栅格对齐。 例子： {‘snapRaster’: {‘url’: ‘<image_service_url>’}} 像元大小(cellSize)：输出栅格将具有由像元大小指定的分辨率。 例子： {‘cellSize’: 11} or {‘cellSize’: {‘url’: <image_service_url>}} or {‘cellSize’: ‘MaxOfIn’} 并行处理因子(parallelProcessingFactor)：控制光栅处理(CPU)服务实例。 例子： 具有指定数量的处理实例的语法示例： {“parallelProcessingFactor”: “2”} 具有指定百分比的总处理实例的语法示例： {“parallelProcessingFactor”: “60%”}
gis	可选的 GIS 对象。如果未指定，则使用当前活动的连接。
future	仅关键字参数。可选的布尔值。如果为 True，则结果将是一个 GPJob 对象，并且结果将异步返回。
folder	仅关键字参数。可选的 str 或 dict。使用给定的文件夹名称在门户中创建一个文件夹(如果不存在)，并将输出保留在此文件夹中。 create_folder() 返回的字典也可以作为输入传入。 例子： {‘username’: ‘user1’, ‘id’: ‘6a3b77c187514ef7873ba73338cf1af8’, ‘title’: ‘trial’}
tiles_only	仅关键字参数。可选的布尔值。在ArcGIS Online 中，此函数的默认输出影像服务将是平铺影像图层。 要在 ArcGIS Online 上创建动态影像图层作为输出，请将 tiles_only 参数设置为 False。 函数将不支持ArcGIS Enterprise 中的tiles_only 参数，默认情况下会生成动态影像图层。

示例 1：

# Usage Example 1: This example performs continuous change detection where only one band is used in the change detection
# and the chi-squared probability threshold is 0.90.
analyze_changes_using_ccdc_op = analyze_changes_using_ccdc(input_multidimensional_raster=input_multidimensional_raster,
                                                           bands_for_detecting_change=[0],
                                                           bands_for_temporal_masking=[],
                                                           chi_squared_threshold=0.99,
                                                           min_anomaly_observations=6,
                                                           update_frequency=1,
                                                           future=False,
                                                           gis=gis
                                                          )

示例 2：

# Usage Example 2: This example performs continuous change detection where bands 3 and 7 (indexed at 2 and 6)
# are used as snow, cloud, and cloud shadow mask.
analyze_changes_using_ccdc_op = analyze_changes_using_ccdc(input_multidimensional_raster=input_multidimensional_raster,
                                                           bands_for_detecting_change=[0,1,2,3,4,5,6],
                                                           bands_for_temporal_masking=[2,6],
                                                           chi_squared_threshold=0.99,
                                                           min_anomaly_observations=3,
                                                           update_frequency=1,
                                                           future=False,
                                                           gis=gis
                                                          )