数字图像处理实验（四，图像分割）（python实现）

import matplotlib.pyplot as plt
import numpy as np


plt.rc("font", family='Microsoft YaHei')
img = plt.imread('../img/4.1.jpg')
img=img[:,:,0]*255 # 像素归一画
plt.subplot(131), plt.title('原始图像')
plt.imshow(img, 'gray')
best_th=1 # 选出来的最好的阈值
best_s=0
bar=[]
for th in range(1,255):#阈值设置为[1,254]
    u1=0
    u1_cnt=0
    u2=0 #计算俩类的灰度值
    u2_cnt=0
    for i in range(img.shape[0]):
        for j in range(img.shape[1]):
            if img[i][j]best_s:
        best_s=s
        best_th=th
img_new=np.zeros(img.shape)
for i in range(img.shape[0]):
    for j in range(img.shape[1]):
        if img[i][j]