arcpy批量以表格显示分区统计（以DEM数据为例）

对tif文件进行区域统计的时候，通常情况下要先利用shp文件进行裁剪，然后再尽区域统计，当然区域统计输出的是tif文件，尽管每个区域有每个区域的值，但是对于数据量大的数据进行统计的时候并不方便，因此可以直接利用以表格显示分区统计。
这里主要记录一下多个shp文件对于一个tif文件的区域统计。

# -*- coding: utf-8 -*-
import arcpy
# 这里是遍历多个shp文件进行区域统计
# shp文件所在文件夹
z_shp = "G:/beijing/buffer_clip/buf_3000/"
arcpy.env.workspace = z_shp
# 搜索该路径中所有的shp文件
shps = arcpy.ListFiles("*.shp")
# dem文件所在路径
dem_file = "G:/beijing/DEM/tif/dem"
# 输出Raster要在的文件夹
out_dem_file = "G:/beijing/DEM/DEM/buf_3000/"
out_dem_statistic = "G:/beijing/DEM/DEM_sta/buf_3000/"
for shp in shps:
    catchment_name = shp[:-4]
    out_raster = out_dem_file + catchment_name
    out_table_sta = out_dem_statistic + catchment_name + ".dbf"
    # 注：这里要加.dbf是因为本来以表格显示分区统计这个功能输出的结果是info表，并不是dbf表，加个后缀可以直接转化，个人认为dbf表比info表好用
    # 裁剪tif文件，这一步也可以去掉，如果你不需要裁剪结果的话
    clip_raster = arcpy.gp.ExtractByMask_sa(dem_file, shp, out_raster)
    # 进行以表格形式显示分区统计
    outZSaT = arcpy.gp.ZonalStatisticsAsTable_sa(shp, "Id", out_raster, out_table_sta, "DATA", "MEAN")
    print "finished"

熟悉arcgis的朋友应该知道在info表的情况下，尽管你是输出多个表，但是在info的格式下就会出现多个表只有一个info文件夹，info文件夹下也没有明显的不同表格，因此如果你要在进行一些其他操作的话其实是有点不太好用的。