Python ArcGIS geometric_median用法及代码示例

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

用法:

arcgis.raster.functions.geometric_median(rasters, epsilon=0.001, max_iteration=10, extent_type='FirstOf', cellsize_type='FirstOf')

返回：

应用了函数的输出栅格。

计算多波段图像时间序列中像素的几何中值。

该函数的参数如下：

Parameter	Description
rasters	Raster / ImageryLayer 对象所需的列表。
epsilon	可选浮点数。指定两次连续迭代之间的收敛值。当epsilon小于或等于指定值时，迭代将停止，并使用上次迭代的结果。
max_iteration	可选整数。指定要完成的最大迭代次数。一旦达到该值，计算就会结束，无论 epsilon 设置如何。
extent_type	可选字符串。指定要用于函数的范围。 “FirstOf” - 使用第一个输入栅格的范围来确定处理范围。这是默认设置。 “IntersectionOf” - 使用重叠像素的范围来确定处理范围。 “UnionOf” - 使用所有栅格的范围来确定处理范围。 “LastOf” - 使用最后一个输入栅格的范围来确定处理范围。
cellsize_type	可选字符串。指定要用于函数的像元大小。 “FirstOf” - 使用输入栅格的第一个像元大小。这是默认设置。 “MinOf” - 使用所有输入栅格中最小的像元大小。 “MaxOf” - 使用所有输入栅格中最大的像元大小。 “MeanOf” - 使用所有输入栅格的平均像元大小。 “LastOf” - 使用输入栅格的最后一个像元大小。

例子：

# Usage Example 1: Calculates the geometric_median function on a list of input rasters.

geometric_median_raster = geometric_median([raster1, raster2, raster3], epsilon=0.001, max_iteration=10)