Python实现迭代阈值法

import matplotlib.pylab as plt
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

def Read_raw(filename,w,h,m):
    #初始化image
    image=[]
    for i in range(m*h):
        image.append([])
        for j in range(w):
            image[i].append(255)
    #读取raw
    with open(filename,'rb') as f:
        data=f.read()
    #写入raw
    for i in range(m*h):
        for j in range(w):
            image[i][j]=data[i*w+j]
    print('读取文件{0}成功'.format(filename))
    return image

def Write_raw(filename,image):
    #将数组image转为bytes
    bt=bytearray()
    for i in range(len(image)):
        for j in range(len(image[i])):
            bt.append(image[i][j])
    #写入raw
    with open(filename,'wb') as f:
        f.write(bt)
    f.close()
    print('写入文件{0}成功'.format(filename))

def iteration(image):
    def itera(image,T):
        T1=[]
        T2=[]
        for i in range(len(image)):
                for j in range(len(image[i])):
                    if image[i][j]=T:
                        T2.append(image[i][j])
        u1=sum(T1)/len(T1)
        u2=sum(T2)/len(T2)
        return u1,u2
                
    T1=255/2
    T0=0.01
    while True:
        u1,u2=itera(image,T1)
        T2=(u1+u2)/2
        if (abs(T1-T2)/(max(T1,T2))