arcgis python 栅格计算 历年最大值

ArcGIS Python提供了强大的栅格计算功能，可以轻松实现历年最大值的计算。下面是使用ArcGIS Python进行历年最大值计算的步骤：

1. 导入所需模块：在Python脚本中，首先需要导入ArcGIS相关的模块，如arcpy和arcpy.sa。

2. 定义输入路径和文件名：指定包含历年栅格数据的文件夹路径和栅格数据的文件名。

3. 创建栅格对象：利用arcpy.sa模块中的Raster函数，将栅格数据加载到内存中，并创建一个栅格对象。

4. 创建一个空的结果栅格：使用arcpy.sa模块中的Con函数，创建一个与输入栅格相同大小的空白栅格，用于存储历年最大值。

5. 循环计算最大值：通过for循环，分别计算每个像元在历年中的最大值。可以使用arcpy.sa模块中的CellStatistics函数，在每次循环中计算历年最大值，并将结果存储到结果栅格中。

6. 保存结果：使用arcpy.sa模块中的Raster函数，将结果栅格保存为新的栅格数据。

以下是一个示例代码，用于计算历年降水数据的最大值：

import arcpy
from arcpy.sa import *

# 定义输入路径和文件名
input_folder = r"C:\data"
input_files = ['rainfall2000.tif', 'rainfall2001.tif', 'rainfall2002.tif']

# 创建空白结果栅格
result_raster = arcpy.sa.Con(Raster(input_files[0]) > 0, Raster(input_files[0]) * 0, -9999)

# 循环计算最大值
for input_file in input_files:
raster = Raster(input_folder + "\\" + input_file)
result_raster = arcpy.sa.CellStatistics([result_raster, raster], "MAXIMUM")

# 保存结果
result_raster.save(input_folder + "\\max_rainfall.tif")

以上代码会将指定文件夹下的降水数据进行历年最大值计算，并将结果保存为max_rainfall.tif文件。你可以根据需要修改输入路径和文件名，以适应自己的数据。