By spotter

2015-09-24 19:40:13 8 Comments

I want to calculate zonal statistics as table for multiple rasters, but can't figure out how. What I am trying is:

import arcpy, os
arcpy.env.workspace=r'F:\Sheyenne\Atmospherically Corrected Landsat\Indices\Main\NDWI\main_master'
#zones to calculate stats on
#where to save output files
destination_path=r'F:\Sheyenne\Atmospherically Corrected Landsat\Indices\Main\NDWI\zonal_stats'
 #list of rasters to calculate statistics on
rasters = arcpy.ListRasters("*water_.tif*")
#start loop
for raster in rasters:
    destination_path = os.path.join(destination_path, raster + ".dbf"), "OBJECTID",raster, destination_path,"DATA","ALL")

The error returned is :

Traceback (most recent call last):

  File "<ipython-input-10-9522776a92dd>", line 1, in <module>
    runfile('F:/python codes/', wdir='F:/python codes')

  File "C:\Users\spotter\AppData\Local\Continuum\Anaconda\lib\site-packages\spyderlib\widgets\externalshell\", line 685, in runfile
    execfile(filename, namespace)

  File "C:\Users\spotter\AppData\Local\Continuum\Anaconda\lib\site-packages\spyderlib\widgets\externalshell\", line 71, in execfile
    exec(compile(scripttext, filename, 'exec'), glob, loc)

  File "F:/python codes/", line 21, in <module>, "OBJECTID", raster, destination_raster,"DATA","ALL")

  File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\geoprocessing\", line 498, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))

ExecuteError: ERROR 999999: Error executing function.
Create output table failed
Failed to execute (ZonalStatisticsAsTable).


@radouxju 2015-09-24 20:13:02

for a single tool, you could also right-click on "zonal stat as table" and run it in batch.

Concerning your code, if it runs without errors nor results, you should first check the layer being processed ("print raster"in the loop), then make sure that there is no projection issue.

@nicksan 2015-09-24 19:50:13

The fourth parameter to ZonalStatisticsAsTable should be the output name - since you've put "DATA" there (which should be parameter 5), you're likely getting no output. You'll need to put your destination_path variable in that location.

But I'll also caution you that, as written, destination_path is going to keep growing since you're appending to it in each round of your loop. Try assigning to a new variable. Something like

for raster in rasters:
    destination_raster = os.path.join(destination_path, raster + ".dbf"), "OBJECTID", raster,destination_raster,"DATA","ALL")

That style should preserve destination_path for each iteration and let you reuse a variable to provide to the Zonal Stats function call.

@spotter 2015-09-24 19:51:08

I actually just noticed that, but some thing, code runs no output

@nicksan 2015-09-24 20:38:52

I'd definitely recommend checking your inputs and your outputs, as mentioned in a few other comments then - add a "print raster" to the loop as well as a "print destination_raster" (or whatever you decided to name it) and make sure that both make sense. You could also try just running it for a single iteration and seeing what happens (effectively removing the loop)

@spotter 2015-09-24 21:02:44

the rasters populate as expected, but you are right, when I run it using only 1 raster without loop it doesn't work either

Related Questions

Sponsored Content

2 Answered Questions

[SOLVED] Error converting Json file to a feature class

0 Answered Questions

RemapRange for Reclassify Error in ArcPy

1 Answered Questions

Arguments error when running LineDensity scripts in arcpy

1 Answered Questions

[SOLVED] Make raster layer arcpy error 000865

2 Answered Questions

[SOLVED] Assertion Error, python

  • 2016-09-29 14:46:44
  • JMap
  • 240 View
  • 1 Score
  • 2 Answer
  • Tags:   arcpy

2 Answered Questions

[SOLVED] Python, Describe.spatialReference errors

1 Answered Questions

1 Answered Questions

Sponsored Content