arcgispro编写python代码使用con函数并且进行按表格形式分区统计

使用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)