ZonalStatisticsAsTable

Question

Hello everyone! I'm starting with Python and I hope someone will help me see through the mistakes I might have made.

I'm making a zona lstatistics as table loop with python and i would like to make ane average of all TIF files pixels according to a grid (.shp).  When I run my script I've the error message bellow:  ERROR 000876: Output raster: E:\climat_data	b_51N_Precip_2014-01-15.dbf's extension is invalid for the output raster format.
Failed to execute (ZonalStatistics).

I'm leaving my script here... maybe someone can help me figure it out:

arcpy.CheckOutExtension("Spatial")
#Set local variables
in_zone_data = 'grid51N.shp'
zone_field = 'id'

raster = arcpy.ListRasters("51N*","ALL")
for stat in raster:
    print (stat)

outDir = r"E:\climat_data"

for stat in raster:
        raster_name = os.path.basename(stat).rstrip(os.path.splitext(stat)[1])
        outTable = os.path.join(outDir,"tb_"+raster_name+".dbf")
        arcpy.gp.ZonalStatistics_sa(in_zone_data, zone_field, stat, outTable, "MEAN", "NODATA")
        print 'Done'+raster_name

Thank you for any help.

Luke_Pinner · Answer

arcpy.sa.ZonalStatistics outputs a raster not a table. You need arcpy.sa.ZonalStatisticsAsTable

IrèneDaubard · Answer

I'm using arcgis pro... I made a model builder and export the script just to know the different tools I've to use. # -*- coding: utf-8 -*- """Generated by ArcGIS ModelBuilder on: 2020-07-24 11:41:30 All ModelBuilder functionality may not be exported. Edits may be required for equivalency with the original model. """ import arcpy # To allow overwriting the outputs change the overwrite option to true. arcpy.env.overwriteOutput = False # Local variables: grid51N_shp__4_ = r"H:\CG_mosq\grid\grid51N.shp" grid51N_shp__10_ = grid51N_shp__4_ grid51N_shp = r"H:\CG_mosq\grid\grid51N.shp" hourlyPrecipRate_2014_01_01_2014_01_14_tif__3_ = r"H:\CG_mosq\climat_data\precipitation\hourlyPrecipRate_2014-01-01_2014-01-14\hourlyPrecipRate_2014-01-01_2014-01-14.tif" hourlyPrecipRate_20140101_20 = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\hourlyPrecipRate_20140101_202" hourlyPrecipRate_20140101_201 = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\hourlyPrecipRate_20140101_203" precip_01_01_reproj_resampled_grid_dbf = r"H:\CG_mosq\climat_data\precipitation\precip-01-01_reproj_resampled_grid.dbf" grid51N_shp__5_ = r"H:\CG_mosq\grid\grid51N.shp" grid51N_shp__7_ = grid51N_shp__5_ grid51N_shp__2_ = r"H:\CG_mosq\grid\grid51N.shp" Tair_f_inst_2014_01_01_2014_01_14_tif = r"H:\CG_mosq\climat_data\air_temperature\Tair_f_inst_2014-01-01_2014-01-14\Tair_f_inst_2014-01-01_2014-01-14.tif" Tair_f_inst_20140101_2014011 = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\Tair_f_inst_20140101_2014011" Tair_f_inst_20140101_20140111 = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\Tair_f_inst_20140101_20140111" tair_01_01_reproj_resampled_grid_dbf = r"H:\CG_mosq\climat_data\air_temperature air-01-01_reproj_resampled_grid.dbf" grid51N_shp__6_ = r"H:\CG_mosq\grid\grid51N.shp" grid51N_shp__9_ = grid51N_shp__6_ grid51N_shp__3_ = r"H:\CG_mosq\grid\grid51N.shp" v3f0acee6d99dbb08f5a84d3bdddf49c4_EVI_median_tif__2_ = r"H:\CG_mosq\climat_data\vegetation_index\EVI_2014-01-01_2014-01-14\3f0acee6d99dbb08f5a84d3bdddf49c4.EVI_median.tif" c3f0acee6d99dbb08f5a84d3bddd = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\c3f0acee6d99dbb08f5a84d3bddd" c3f0acee6d99dbb08f5a84d3bddd1 = r"H:\CG_mosq\MyProject_arcgispro_modelbuilder\MyProject_arcgispro_modelbuilder.gdb\c3f0acee6d99dbb08f5a84d3bddd1" evi_01_01_reproj_resampled_grid_dbf = r"H:\CG_mosq\climat_data\vegetation_index\evi-01-01_reproj_resampled_grid.dbf" sr = arcpy.SpatialReference(32651)#"WGS 84 UTM zone 51N" # Process: Project Raster (2) arcpy.ProjectRaster_management(in_raster=hourlyPrecipRate_2014_01_01_2014_01_14_tif__3_, out_raster=hourlyPrecipRate_20140101_20, out_coor_system="PROJCS['WGS_1984_UTM_Zone_51N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',123.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", resampling_type="NEAREST", cell_size="11171.1925206363 11171.1925206363", geographic_transform="", Registration_Point="", in_coor_system="GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]", vertical="NO_VERTICAL") # Process: Resample (2) arcpy.Resample_management(in_raster=hourlyPrecipRate_20140101_20, out_raster=hourlyPrecipRate_20140101_201, cell_size="4000 4000", resampling_type="NEAREST") # Process: Zonal Statistics as Table (2) arcpy.gp.ZonalStatisticsAsTable_sa(in_zone_data=grid51N_shp, zone_field="id", in_value_raster=hourlyPrecipRate_20140101_201, out_table=precip_01_01_reproj_resampled_grid_dbf, ignore_nodata="DATA", statistics_type="MEAN") # Process: Join Field arcpy.JoinField_management(in_data=grid51N_shp__4_, in_field="id", join_table=precip_01_01_reproj_resampled_grid_dbf, join_field="MEAN", fields="") # Process: Project Raster arcpy.ProjectRaster_management(in_raster=Tair_f_inst_2014_01_01_2014_01_14_tif, out_raster=Tair_f_inst_20140101_2014011, out_coor_system=sr) # Process: Resample arcpy.Resample_management(in_raster=Tair_f_inst_20140101_2014011, out_raster=Tair_f_inst_20140101_20140111, cell_size="4000 4000", resampling_type="NEAREST") # Process: Zonal Statistics as Table arcpy.gp.ZonalStatisticsAsTable_sa(in_zone_data=grid51N_shp__2_, zone_field="id", in_value_raster=Tair_f_inst_20140101_20140111, out_table=tair_01_01_reproj_resampled_grid_dbf, ignore_nodata="DATA", statistics_type="MEAN") Process: Join Field (2) arcpy.JoinField_management(in_data=grid51N_shp__5_, in_field="id", join_table=tair_01_01_reproj_resampled_grid_dbf, join_field="id", fields="") # Process: Project Raster (3) arcpy.ProjectRaster_management(in_raster=v3f0acee6d99dbb08f5a84d3bdddf49c4_EVI_median_tif__2_, out_raster=c3f0acee6d99dbb08f5a84d3bddd, out_coor_system="PROJCS['WGS_1984_UTM_Zone_51N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',123.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", resampling_type="NEAREST", cell_size="464.671253451976 464.671253451976", geographic_transform="", Registration_Point="", in_coor_system="GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]", vertical="NO_VERTICAL") # Process: Resample (3) arcpy.Resample_management(in_raster=c3f0acee6d99dbb08f5a84d3bddd, out_raster=c3f0acee6d99dbb08f5a84d3bddd1, cell_size="4000 4000", resampling_type="NEAREST") # Process: Zonal Statistics as Table (3) arcpy.gp.ZonalStatisticsAsTable_sa(in_zone_data=grid51N_shp__3_, zone_field="id", in_value_raster=c3f0acee6d99dbb08f5a84d3bddd1, out_table=evi_01_01_reproj_resampled_grid_dbf, ignore_nodata="DATA", statistics_type="MEAN") # Process: Join Field (3) arcpy.JoinField_management(in_data=grid51N_shp__6_, in_field="id", join_table=evi_01_01_reproj_resampled_grid_dbf, join_field="id", fields="") ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ I'll switch the "51N_"+data and get back to you just in case... I currently no longer have access to my office. Thank you TanuHoque_0-1636058088983.png

DanPatterson · Answer

Are you using arcmap or pro? your print statements aren't consistent with python 3. Justd as an FYI..., esri grids can't begin with a number "51N_"+data so that would have been an issue without the *.tif file extension.