ArcGIS Pro python 获取一百多幅栅格的平均值

目标

我需要计算流域内的平均等效水柱高，等效水柱高的格网已经创建好，如下：

每个月一个这样的栅格格网，现在有近两百个格网需要求平均等效水柱高
要求：按照时间顺序求取每一个格网的平均等效水柱高，并生成时间序列表单，如下：

看到arcpro上显示栅格资料时，有平均值这一属性，于是打算用arcpy实现它

步骤

打开pycharm，创建项目，修改python环境为arcpro的python环境

问题1：import arcpy 报错
E:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\lib\site-packages\numpy_init_.py:140: UserWarning: mkl-service package failed to import, therefore Intel® MKL initialization ensuring its correct out-of-the box operation under condition when Gnu OpenMP had already been loaded by Python process is not assured. Please install mkl-service package, see http://github.com/IntelPython/mkl-service from . import _distributor_init
解决：将py文件放入include中

代码：

# -*- coding:UTF-8 -*- # 要对字符编码格式进行修改，否则无法让中文合法
# *************************************************************************
# 目的：读取一百多个栅格的平均值属性
# 时间：2021年3月6日
# 作者：vivianooooo
# 思路：通过获取栅格属性函数获取栅格的cell mean值
# 创建一个数组存储时间和mean值，组成时间序列，并输出到txt文件
# *************************************************************************

import arcpy
from arcpy import env

def Model3():
    arcpy.env.workspace = "D:\\data"
    # 获取工作空间中的栅格列表
    grid_list = arcpy.ListRasters("*","GRID")
    file = open('D:\\data\\tws.txt', 'w')
    file.write('time' + '\t' + 'tws_cm'+'\n')
    # 迭代栅格
    for raster in grid_list :
    	# 获取栅格数据属性
        Property = arcpy.management.GetRasterProperties(in_raster=raster, property_type="MEAN", band_index=raster)[0]
        # 取原文件名的标识时间部分输出时间序列值
        file.write(raster[5:12]+'\t'+Property+'\n')

if __name__ == '__main__':
    Model3()

结果

arcpy.management.GetRasterProperties函数说明