# -*- coding: UTF-8 -*-
import os
import arcpy
from arcpy.sa import *
import csv
arcpy.CheckOutExtension("ImageAnalyst") # 检查许可
arcpy.CheckOutExtension("spatial")
# 输入路径
inws = r"E:\Data\grassGpp\allgrassGpp\year"
# 文本txt
# 更改要存储结果的路径,自动生成表格
# OutputFile = open(r"E:\Data_Shuju\mean\allgrass/GppMean.csv".decode('utf-8'), 'w')
# 输出路径 把文件名和参数写入csv py 2用的wb 3用w newline='' 创建一个csv文件
mycsv = open(r"E:\Data_Shuju\mean\allgrass\GppMean2.csv", 'wb') # 新建并打开csv
csv_write = csv.writer(mycsv, dialect='excel')
# 获取路径下的tif
arcpy.env.workspace = inws
rasters = arcpy.ListRasters("*", "tif")
print(rasters)
'''
# 利用glob包,将inws下的所有tif文件读存放到rasters中
rasters = glob.glob(os.path.join(inws, '*.tif'))
path = os.path.join(inws, '*.tif')
print(path)
'''
# whereClause = "VALUE = -9999999" # 去除异常值。如果无异常值可删去.
# 循环rasters中的所有影像,进行“求平均值”操作
for ras in rasters:
# outSetNull = SetNull(ras, ras, whereClause) # 去除异常值;如果无异常值可删去
meanValueInfo = arcpy.GetRasterProper
获取文件夹下所有TIF的均值导出为csv文件
最新推荐文章于 2025-08-27 22:56:26 发布

最低0.47元/天 解锁文章
2287

被折叠的 条评论
为什么被折叠?



