为了解决Arcgis进行自然间断点无法批量处理,本人结合Python实现对原始的excel数据进行批量自然间断点处理,希望可以帮助一些人,从重复性劳力中解放出来!
import xlrd
import xlwt
import jenkspy
Table= xlrd.open_workbook(excel数据路径)
#打印Table相关信息
print(Table.sheet_names())
sheet1 = Table.sheets()[0]
print(sheet1.name,sheet1.nrows, sheet1.ncols)
print(sheet1.col_values(1))
f = xlwt.Workbook()
#遍历Table的每个sheet
for k in range(0,len(Table.sheets())):
sheetSave = f.add_sheet(Table.sheets()[k].name,cell_overwrite_ok = True)
sheet = Table.sheets()[k]
#遍历sheet的列数据
for j in range(1,sheet.ncols):
col = sheet.col_values(j)
List = col[1:len(col)]
ListID = col[1:len(col)]#对类别进行标识
#计算间断点
breaks = jenkspy.jenks_breaks(List,nb_class = 8)
#根据计算的间断点,判断每个原始数据属于哪个类别
for i in range(len(List)):
if List[i] < breaks[1]:
ListID[i] = 1
elif List[i] < breaks[2] and List[i] >= breaks[1]:
ListID[i] = 2
elif List[i] < breaks[3] and List[i] >= breaks[2]:
ListID[i] = 3
elif List[i] < breaks[4] and List[i] >= breaks[3]:
ListID[i] = 4
elif List[i] < breaks[5] and List[i] >= breaks[4]:
ListID[i] = 5
elif List[i] < breaks[6] and List[i] >= breaks[5]:
ListID[i] = 6
elif List[i] < breaks[7] and List[i] >= breaks[6]:
ListID[i] = 7
else:
ListID[i] = 8
print(breaks)
print(ListID)
for u in range(0,len(ListID)):
sheetSave.write(u,j,ListID[u])#写数据
print(j,'------')
print(k,'======')
#保存数据成excel表
f.save(r'D:\demo.xls')
评论