现需要对属性表中两个字段分别进行阈值划分,并组合判断不同类别的数量,如果在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("")