Arcpy 实现NDVI最大值合成/平均值合成

当研究NDVI产品数据时，需要日产品月最大值合成/月平均值合成，或者月产品数据季度最大值合成/平均值合成，ArcGIS很方便的实现这些功能，Arcpy辅助更加便捷的处理批量数据。

首先介绍一下Arcpy函数CellStatistics（像元统计）：

根据多个栅格数据计算每个像元的统计数据。

可用的统计数据有：众数、最大值、均值、中位数、最小值、少数、范围、标准差、总和及变异度。

用法

• 使用此工具时输入栅格的顺序不重要。

• 将多波段栅格指定为输入栅格或常量值（Python 中的 in_rasters_or_constants）之一时，将使用所有波段。

  要处理一系列来自多波段栅格的波段，可先用波段合成工具创建由这些特定波段组成的新栅格数据集，并在输入栅格或常量值（Python 中的 in_rasters_or_constants）列表中使用该结果。

• 对于最大值、最小值、均值、中位数、众数、少数及总和这些统计数据类型，如果将单个栅格数据用作输入，则输出像元值将与输入像元值相同。对于范围和 STD这些统计数据类型，输出像元值全部为 0；对于变异度这种统计数据，输出像元值均为 1。

python源代码：

outCellStatistics = CellStatistics([inRaster01, inRaster02, inRaster03], "RANGE", "NODATA")
outCellStatistics.save("C:/sapyexamples/output/cellstats")

至关重要(加在关键代码之前)：

arcpy.CheckOutExtension("Spatial")

效果展示：

 

平均值合成：