blood?=?imread（‘blood1.BMP‘）;
[xy]=size（blood）;????????????????%?求出圖象大小
b=double（blood）;??????????????????
N?=sqrt（100）?*?randn（xy）;????????%??生成方差為10的白噪聲
I=b+N;????????????????????????????%??噪聲干擾圖象

for?i=1:x?????????????????????????%?實際圖象的灰度為0～255
????for?j=1:y
????????if?（I（ij）>255）
????????????I（ij）=255;
????????end
????????if?（I（ij）<0）
????????????I（ij）=0;
????????end
????end
end????


z0=max（max（I））;???????????????????%?求出圖象中最大的灰度
z1=min（min（I））;???????????????????%?最小的灰度?

T=（z0+z1）/2;??????????????????????
TT=0;
S0=0;?n0=0;
S1=0;?n1=0;
allow=0.5;???????????????????????%?新舊閾值的允許接近程度
d=abs（T-TT）;
count=0;?????????????????????????%?記錄幾次循環

while（d>=allow）?????????????????%?迭代最佳閾值分割算法
????count=count+1;
????for?i=1:x
????????for?j=1:y
????????????if?（I（ij）>=T）
????????????????S0=S0+I（ij）;
????????????????n0=n0+1;
????????????end
????????????if?（I（ij）????????????????S1=S1+I（ij）;
????????????????n1=n1+1;
????????????end
????????end
????end?
????T0=S0/n0;
????T1=S1/n1;
????TT=（T0+T1）/2;
????d=abs（T-TT）;
????T=TT;
end

Seg=zeros（xy）;
for?i=1:x
????for?j=1:y
????????if（I（ij）>=T）
????????????Seg（ij）=1;???????????????%?閾值分割的圖象
????????end
????end
end

SI=1-Seg;????????????????????????????%?閾值分割后的圖象求反，便于用腐蝕算法求邊緣
se1=strel（‘square‘3）;???????????????%?定義腐蝕算法的結構
SI1=imerode（SIse1）;?????????????????%?腐蝕算法
BW=SI-SI1;???????????????????????????%?邊緣檢測
%=====傳統的邊緣檢測方法======%
I=uint8（I）;
BW1=edge（I‘sobel‘?0.11）;
BW2=edge（I‘log‘?0.015）;
%===========圖象顯示==========%
figure（1）;
imshow（I）;title（‘Original‘）?????????????????????????%?顯示閾值分割的圖象
figure（2）;
imshow（BW2）;title（‘Soble‘）?????????????????????????%?顯示新算法的邊緣圖象
figure（3）
imshow（BW1）;title（‘Gauss-Laplace‘）?
figure（4）
imshow（BW）;title（‘New?algorithm‘）?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????73160??2003-11-08?18:39??細胞圖像分割源代碼\blood1.BMP

?????文件???????1962??2004-03-28?16:19??細胞圖像分割源代碼\cancer.m

????..A.SH.?????11264??2011-10-14?14:51??細胞圖像分割源代碼\Thumbs.db

?????目錄??????????0??2011-10-14?14:51??細胞圖像分割源代碼

-----------?---------??----------?-----??----

????????????????86386????????????????????4