lung_test.m

tic
I=imread（‘C:\Users\zdh\Desktop\bishe\image\9.jpg‘）;
% figureimshow（I）title（‘原图像‘）;
I=rgb2gray（I）;
I1=double（I）;
T=（min（I1（:））+max（I1（:）））/2;
d=false;
i=0;
while~d
    u1=find（I1<=T）;
    u2=find（I1>T）;
    Tn=（mean（I1（u1））+mean（I1（u2）））/2;
    d=abs（Tn-T）<1;
    T=Tn;
    i=i+1;
end
I1（u1）=0;
I1（u2）=1;
% figureimshow（I1）title（‘二值图像‘）;
se=strel（‘disk‘18）;
A=imerode（I1se）