Python ArcGIS band_arithmetic用法及代码示例

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

用法:

arcgis.raster.functions.band_arithmetic(raster, band_indexes=None, astype=None, method=0)

返回：

应用了函数的输出栅格。

band_arithmetic 函数对栅格的波段执行算术运算。欲了解更多信息，请参阅Band Arithmetic function

Parameter	Description
raster	需要输入 Raster / ImageryLayer 对象。
band_indexes	可选字符串或列表。带索引可以以空格分隔的字符串或整数或浮点值列表的形式给出。例如，“4 3”或 [4,3]。对于用户定义的方法，波段索引可以以表达式形式给出，例如“(B3 - B1)/(B3 + B1)”。
astype	可选字符串。指定输出像素类型。可用选项是 - “C128” | “C64” | “F32” | “F64” | “S16” | “S32” | “S8” | “U1” | “U16” | “U2” | “U32” | “U4” | “U8”。默认为无。
method	可选的整数。您要部署的波段算术算法的类型。您可以定义自定义算法，或选择预定义的索引。 0 = UserDefined 1 = NDVI 2 = SAVI 3 = TSAVI 4 = MSAVI 5 = GEMI 6 = PVI 7 = GVITM 8 = Sultan 9 = VARI 10 = GNDVI 11 = SR 12 = NDVIre 13 = SRre 14 = MTVI2 15 = RTVICore 16 = CIre 17 = CIg 18 = NDWI 19 = EVI 20 = IronOxide 21 = FerrousMinerals 22 = ClayMinerals 23 = WNDWI 24 = BAI 25 = NBR 26 = NDBI 27 = NDMI 28 = NDSI 29 = MNDW

例子：

# Usage Example 1: Apply the 'NDVI' band arithmetic algorithm on the input raster.

band_op = band_arithmetic(raster, method=1)

# Usage Example 2: Apply user defined band arithmetic algorithm on the input raster.

band_operation = "(B3 - B1)/(B3 + B1)"
band_op = band_arithmetic(raster, band_indexes=band_operation, method=0)