Sitk图像分割

import cv2
import SimpleITK as sitk
import matplotlib.pylab as pylab
def show_image(img,title=None):
    nda=sitk.GetArrayViewFromImage(img)
    pylab.imshow(nda,cmap='gray')
    pylab.axis('off')
    if(title):
        pylab.title(title,size=20)
img=cv2.imread('C:/Users/xpp/Desktop/Lena.png')
gray=cv2.cvtColor(img,cv2.COLOR_BGR2GRAY)#将彩色图片转换为灰度图片
img=255*gray
img_T1=sitk.GetImageFromArray(img)
img_T1_255=sitk.Cast(sitk.RescaleIntensity(img_T1),sitk.sitkUInt8)#sitk分割
seed=(100,120)#种子
for upper in [80,85,90]:
    pylab.figure(figsize=(18,20))
    pylab.subplot(221)
    show_image(img_T1,"Original Image")
    pylab.scatter(seed[0],seed[1],color='red',s=50)
    pylab.subplot(222)
    seg=sitk.ConnectedThreshold(img_T1,seedList=[seed],lower=40,upper=upper)
    show_image(seg,"Region Growing")
    pylab.subplot(223)
    show_image(sitk.LabelOverlay(img_T1_255,seg),"Connected Threshold")
    pylab.axis('off')
    pylab.tight_layout()
    pylab.show()