I am attempting to calculate zonal statistics on a single raster for multiple shapefiles. When I load all the shapefiles into the contents pane the code runs (tested on 5 files), but I have thousands of shapefiles I need to conduct this calculation and cannot load them all into the contents pane. When I dont load the files to the contents pane only the zonal stat for the first file is calculated, after that an error is returned:
ExecuteError Traceback (most recent call last)
In [10]:
Line 23: outzone = ZonalStatistics(shape[i], "presence", Raster("C:pathway to file/.tif"), "MAXIMUM", "NODATA")
File c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Functions.py, in ZonalStatistics:
Line 6530: process_as_multidimensional)
File c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Utils.py, in swapper:
Line 53: result = wrapper(*args, **kwargs)
File c:\program files\arcgis\pro\Resources\arcpy\arcpy\sa\Functions.py, in Wrapper:
Line 6522: process_as_multidimensional)
File c:\program files\arcgis\pro\Resources\arcpy\arcpy\geoprocessing\_base.py, in <lambda>:
Line 511: return lambda *args: val(*gp_fixargs(args, True))
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000865: Input raster or feature zone data: Abrawayaomys_chebezi.shp does not exist.
WARNING 001000: Zone field: Field presence does not exist
Failed to execute (ZonalStatistics).
The full code is listed below:
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("Spatial")
indir = "C:/pathway to folder"
outdir = "C:/pathway to output/"
env.workspace = indir
shape = arcpy.ListFiles("*.shp")
i = 0
while i < len(shape):
outzone = ZonalStatistics(shape[i], "presence", Raster("pathway to raster"), "MAXIMUM", "NODATA")
outzone.save(outdir+shape[i]+"_Jan_Max.tif")
del outzone
i += 1
print ("Sampling iteration #" + str(i))
print ("Complete")
the pathways were changed just for this question.
is there any suggestions for how to avoid this error?
