阈值分割及形态学处理.doc

阈值分割及形态学处理：
close all;
clear all;
I=imread('');
figure
imshow(I),title('原图')
level = graythresh(I);
bw = 阈值分割及形态学处理：
close all;
clear all;
I=imread('');
figure
imshow(I),title('原图')
level = graythresh(I);
bw = im2bw(I,level);
se = strel('disk',5);
bw2 = imopen(bw,se);
figure
imshow(bw),title('Thresholded image')
figure
imshow(bw2),title('After opening')
[L,NUM] = bwlabel(bw2,4);
figure
hist(L);title('L的直方图')
[wid,hig]=size(L);
num=1;
mm1=0;mm2=0;
for m=1:NUM
for i=1:wid
for j=1:hig
if L(i,j)==m,mm2=mm2+1;
end
end
end
if mm2>mm1,num=m;
end
mm1=max(mm1,mm2);
mm2=0;
end %寻找最大num循环
for i=1:wid
for j=1:hig
if L(i,j)==num,bw2(i,j)=1;
else bw2(i,j)=0;
end
end
end
figure
imshow(bw2),title('最大连通成分提取')
I1=rgb2gray(I);
K=immultiply(I1,bw2);
figure
imshow(K),title('最后结果')
图谱分析与所得结果: