For循环中区域统计中的多个栅格？答案

【问题标题】：Multiple Rasters in Zonal Statistics in a For-loop?For循环中区域统计中的多个栅格？
【发布时间】：2016-08-27 10:34:18
【问题描述】：

我正在为我的初学者 Python 课程做一个项目，但遇到了一点困难。我有三个要进行区域统计的 .tif 文件，但出现错误。这是我的脚本：

import arcpy
import os
from arcpy import env
from arcpy.sa import *
env.workspace = r'C:\Users\alvaremi\Documents\Final Project_Python'
path = r'C:\Users\alvaremi\Documents\Final Project_Pythonn'
env.overwriteOutput = 1
arcpy.CheckOutExtension('Spatial')

in_zone_data = 'counties_in_cog.shp'
zone_field = 'NAME'

impervious = os.listdir(env.workspace + '\ImpvClipped')
print impervious

for year in impervious:
  if year.endswith(".tif"):
    outZonalStatistics = ZonalStatistics(in_zone_data, zone_field, year, "MEAN", "NODATA")
    outZonalStatistics.save(year[:8] + 'zonalstats')

print 'Done'

当我运行它时，我得到这个错误： ExecuteError：执行失败。参数无效。错误 000865：输入值栅格：2001impvclipped.tif 不存在。执行失败（区域统计）。

我也不确定如何保存新文件以便它们保留日期。我要运行区域统计的文件是“2001impclipped”、“2006impclipped”和“2011impclipped”。

谢谢！

【问题讨论】：

标签： python loops for-loop arcpy

【解决方案1】：

您需要将完整的目录路径添加到文件名中，以便 Python 找到它。

fileName = env.workspace + '\ImpvClipped\' + year
ZonalStatistics(in_zone_data, zone_field, fileName, "MEAN", "NODATA")

【讨论】：