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
                                                          )