使用Arcgis Pro编写python代码进行分区统计,如果任务比较多的话很方便。Arcgis Pro永远的神,饱受Arcmap折磨的伙伴可以考虑试试(本代码没有在Arcmap中测试,不知道能不能用)
import numpy as np
shengs = ['JL_','LN_']#省份
datas = ['ENN','LSI','I'] #数据后缀名
crops = ["corn","rice","soybean"]#作物类型
for sheng in shengs:
for data in datas:
for crop in crops:
data_path = sheng+data+".dat"# 赋值栅格
outpath = r"输入你的文件输出路径"
outpathtable1 = r"表格1路径"
outpathtable2 = r"表格2路径(最后加上.csv)"
classdata = "条件栅格"
print(data_path)
# con函数,满足classdata>0的部分赋予data_path的值(注意如果是等于不要用两个等号,一个等号就可以了)
output_raster = arcpy.sa.Con(classdata,data_path,where_clause="value>0")
output_raster.save(outpath)
#按表格形式显示分区统计,如果实在不知道参数怎么设置就自己先运行一个,在历史消息中右键复制python命令就有python代码了
arcpy.sa.ZonalStatisticsAsTable("你的统计矢量", "统计矢量字段", outpath, outpathtable1, "DATA", "MEAN(统计均值可以修改)", "CURRENT_SLICE", [90], "AUTO_DETECT", "ARITHMETIC", 360)
#导出为csv表格
arcpy.ExportTable_conversion(outpathtable1, outpathtable2)