<?xml version="1.0" encoding="UTF-8"?>
<rss xmlns:content="http://purl.org/rss/1.0/modules/content/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:taxo="http://purl.org/rss/1.0/modules/taxonomy/" version="2.0">
  <channel>
    <title>topic Can't save Zonal Statistics As Table in Geoprocessing Questions</title>
    <link>https://community.esri.com/t5/geoprocessing-questions/can-t-save-zonal-statistics-as-table/m-p/681108#M22618</link>
    <description>&lt;HTML&gt;&lt;HEAD&gt;&lt;/HEAD&gt;&lt;BODY&gt;&lt;SPAN&gt;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?&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;Here is the code:&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;import arcpy&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;import sys, os, string&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;from arcpy.sa import *&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;from arcpy import env&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;#Get an specific raster image file from the folder&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;LstRst = arcpy.ListRasters("*", "img")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;fileprefix = "aug07_ndvi"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;for Rst in LstRst:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; Img1 = string.find(Rst, fileprefix + "_stanislaus.img")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; if Img1 &amp;gt; -1:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; print Rst&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt; #Get an specific field from a feature class file&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;LstFC = arcpy.ListFeatureClasses("*")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;for fc in LstFC:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print fc&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;inputFC = "StanislausAgFields.shp"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;field = "Field_ID"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;lstFlds = arcpy.ListFields(inputFC, field)&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;if len(lstFlds) == 0:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddField_management(inputFC, field, "Text")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;else:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print "field already exists"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.CheckOutExtension("Spatial")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.overwriteOutput = 1&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;try:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; arcpy.AddMessage("Processing zonal means...")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; OutZsaT1 = ZonalStatisticsAsTable(inputFC, field, Rst, "outTable1", "DATA", "MEAN")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;except:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddMessage("Problem in attempting to calculate NDVI means. Check for existence of rasters")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print arcpy.GetMessages()&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddMessage(arcpy.GetMessages())&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print arcpy.GetMessages()[/SIZE]&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;I tried both ways to save zonal statistics using the followings:&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;OutZsaT1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;However, I don't see "aug07_zonalsat.dbf" table in that directory. &lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;Would appreciate for your help. thanks&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;RG&lt;/SPAN&gt;&lt;/BODY&gt;&lt;/HTML&gt;</description>
    <pubDate>Wed, 18 Jan 2012 16:19:12 GMT</pubDate>
    <dc:creator>RameshGautam</dc:creator>
    <dc:date>2012-01-18T16:19:12Z</dc:date>
    <item>
      <title>Can't save Zonal Statistics As Table</title>
      <link>https://community.esri.com/t5/geoprocessing-questions/can-t-save-zonal-statistics-as-table/m-p/681108#M22618</link>
      <description>&lt;HTML&gt;&lt;HEAD&gt;&lt;/HEAD&gt;&lt;BODY&gt;&lt;SPAN&gt;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?&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;Here is the code:&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;import arcpy&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;import sys, os, string&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;from arcpy.sa import *&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;from arcpy import env&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;#Get an specific raster image file from the folder&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;LstRst = arcpy.ListRasters("*", "img")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;fileprefix = "aug07_ndvi"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;for Rst in LstRst:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; Img1 = string.find(Rst, fileprefix + "_stanislaus.img")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; if Img1 &amp;gt; -1:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; print Rst&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt; #Get an specific field from a feature class file&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;LstFC = arcpy.ListFeatureClasses("*")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;for fc in LstFC:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print fc&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;inputFC = "StanislausAgFields.shp"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;field = "Field_ID"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;lstFlds = arcpy.ListFields(inputFC, field)&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;if len(lstFlds) == 0:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddField_management(inputFC, field, "Text")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;else:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print "field already exists"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; &lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.CheckOutExtension("Spatial")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;arcpy.env.overwriteOutput = 1&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;try:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; arcpy.AddMessage("Processing zonal means...")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; OutZsaT1 = ZonalStatisticsAsTable(inputFC, field, Rst, "outTable1", "DATA", "MEAN")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp; outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;except:&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddMessage("Problem in attempting to calculate NDVI means. Check for existence of rasters")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print arcpy.GetMessages()&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; arcpy.AddMessage(arcpy.GetMessages())&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;&amp;nbsp;&amp;nbsp;&amp;nbsp; print arcpy.GetMessages()[/SIZE]&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;I tried both ways to save zonal statistics using the followings:&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;outTable1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;OutZsaT1.save("D:\Project\Stanislaus\Timeseries\aug07_zonalstat.dbf")&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;However, I don't see "aug07_zonalsat.dbf" table in that directory. &lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;Would appreciate for your help. thanks&lt;/SPAN&gt;&lt;BR /&gt;&lt;BR /&gt;&lt;SPAN&gt;RG&lt;/SPAN&gt;&lt;/BODY&gt;&lt;/HTML&gt;</description>
      <pubDate>Wed, 18 Jan 2012 16:19:12 GMT</pubDate>
      <guid>https://community.esri.com/t5/geoprocessing-questions/can-t-save-zonal-statistics-as-table/m-p/681108#M22618</guid>
      <dc:creator>RameshGautam</dc:creator>
      <dc:date>2012-01-18T16:19:12Z</dc:date>
    </item>
    <item>
      <title>Re: Can't save Zonal Statistics As Table</title>
      <link>https://community.esri.com/t5/geoprocessing-questions/can-t-save-zonal-statistics-as-table/m-p/681109#M22619</link>
      <description>&lt;HTML&gt;&lt;HEAD&gt;&lt;/HEAD&gt;&lt;BODY&gt;&lt;SPAN&gt;For every file path listed in the script, eg.&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;"D:\Project\Stanislaus\Timeseries\NDVI"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;you need to replace the backslashes with forward slashes or double backslashes ( \\ ) or place an r before the string ie.&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;r"D:\Project\Stanislaus\Timeseries\NDVI"&lt;/SPAN&gt;&lt;BR /&gt;&lt;SPAN&gt;check the python documentation on string requirements and escape codes.&lt;/SPAN&gt;&lt;/BODY&gt;&lt;/HTML&gt;</description>
      <pubDate>Wed, 18 Jan 2012 19:35:09 GMT</pubDate>
      <guid>https://community.esri.com/t5/geoprocessing-questions/can-t-save-zonal-statistics-as-table/m-p/681109#M22619</guid>
      <dc:creator>DanPatterson_Retired</dc:creator>
      <dc:date>2012-01-18T19:35:09Z</dc:date>
    </item>
  </channel>
</rss>

