I am trying to save Zonal Statistics Table, which is a mean of a raster images based on fields identified using Field_ID. Although, the code runs well, but it does not produce output result. Can you please find what???s wrong with this code?
Here is the code:
import arcpy
import sys, os, string
from arcpy.sa import *
from arcpy import env
#Get an specific raster image file from the folder
arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
LstRst = arcpy.ListRasters("*", "img")
fileprefix = "aug07_ndvi"
for Rst in LstRst:
Img1 = string.find(Rst, fileprefix + "_stanislaus.img")
if Img1 > -1:
print Rst
#Get an specific field from a feature class file
arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries"
LstFC = arcpy.ListFeatureClasses("*")
for fc in LstFC:
print fc
inputFC = "StanislausAgFields.shp"
field = "Field_ID"
lstFlds = arcpy.ListFields(inputFC, field)
if len(lstFlds) == 0:
arcpy.AddField_management(inputFC, field, "Text")
else:
print "field already exists"
arcpy.CheckOutExtension("Spatial")
arcpy.env.overwriteOutput = 1
try:
arcpy.AddMessage("Processing zonal means...")
OutZsaT1 = ZonalStatisticsAsTable(inputFC, field, Rst, "outTable1", "DATA", "MEAN")
outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")
except:
arcpy.AddMessage("Problem in attempting to calculate NDVI means. Check for existence of rasters")
print arcpy.GetMessages()
arcpy.AddMessage(arcpy.GetMessages())
print arcpy.GetMessages()[/SIZE]
I tried both ways to save zonal statistics using the followings:
outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")
OutZsaT1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")
However, I don't see "aug07_zonalsat.dbf" table in that directory.
Would appreciate for your help. thanks
RG
