基于Python实现自然间断点批量处理

对excel数据进行批量自然间断点处理

为了解决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')
  • 作者:seibert(联系作者)
  • 发表时间:2020-05
  • 版权声明:未经站长允许,不可转载
  • 公众号转载:请在文末添加作者公众号二维码
  • 评论

    可能感兴趣内容