【Arcpy】-【属性表运算】按条件分类并计数-双重条件选择

现需要对属性表中两个字段分别进行阈值划分，并组合判断不同类别的数量，如果在Arcgis里用按属性选择将进行多次繁琐操作，所以考虑用代码完成。

实现思路是进行双重条件判断，通过循环计数。

代码实例：

1. 设置工作环境和阈值范围的字典

# encoding: utf-8
import arcpy

arcpy.env.workspace = r"E:\02-Nuclear_Pop_Risk\03-2001_2020World\02-PlantYearSHP_FB_FB\Buchong\Tongji"
shapefile_path = r"BC2015-2020.shp"
col_name = '20'#这里是我的字段名称需要

#设置阈值
count_dict = {wn_range: {'< -0.02': 0, '-0.02 to 0': 0, '0 to 0.04': 0, '0.04 to 0.08': 0, '> 0.08': 0}
              for wn_range in ['0-100', '100-500', '500-1000', '1000-2000', '2000-']}

我这里的先决条件是wn_range 在['0-100', '100-500', '500-1000', '1000-2000', '2000-']这几个范围内，再计算第二个条件属于'< -0.02'，'-0.02 to 0'……来进行计数。

2. 计算

# 使用SearchCursor遍历特征
with arcpy.da.SearchCursor(shapefile_path, ["WN_MulHaz{}".format(col_name), "VPD{}_50".format(col_name)]) as cursor:
    for row in cursor:
        WN_MulHaz, vpd50_value = row
        # 确定WN_MulHaz的范围
        if 0 <= WN_MulHaz< 100:
            wn_range = '0-100'
        elif 100 <= WN_MulHaz< 500:
            wn_range = '100-500'
        elif 500 <= WN_MulHaz< 1000:
            wn_range = '500-1000'
        elif 1000 <= WN_MulHaz< 2000:
            wn_range = '1000-2000'
        elif 2000 <= WN_MulHaz:
            wn_range = '2000-'
        else:
            continue  # 如果不在指定范围内，则跳过

        # 根据VPD50的值分类并计数
        if vpd50_value < -0.02:
            count_dict[wn_range]['< -0.02'] += 1
        elif -0.02 <= vpd50_value < 0:
            count_dict[wn_range]['-0.02 to 0'] += 1
        elif 0 < vpd50_value <= 0.04:
            count_dict[wn_range]['0 to 0.04'] += 1
        elif 0.04 < vpd50_value <= 0.08:
            count_dict[wn_range]['0.04 to 0.08'] += 1
        elif vpd50_value > 0.08:
            count_dict[wn_range]['> 0.08'] += 1

3. 结果输出

for wn_range, vpd50_counts in count_dict.items():
    print("Haz Range: {}".format(wn_range))
    for vpd50_range, count in vpd50_counts.items():
        print("VPD50 {}: {}".format(vpd50_range,count))
    print("")