chiduxuanze.m

用于选择结构元素的大小

%%% 尺度选择 program %%%
clear 
close all
%I=imread('f:\X\骨折图片\17.jpg');
I=imread('22.jpg');
I=imresize(I,[512 512]);
I=rgb2gray(I);
imshow(I)
figure,imhist(I)
se0=strel('disk',1);
t0=imtophat(I,se0);
se1=strel('disk',2);
t1=imtophat(I,se1);
se2=strel('disk',3);
t2=imtophat(I,se2);
se3=strel('disk',4);
t3=imtophat(I,se3);
se4=strel('disk',5);
t4=imtophat(I,se4);
se5=strel('disk',6);
t5=imtophat(I,se5);
se6=strel('disk',7);
t6=imtophat(I,se6);
se7=strel('disk',8);
t7=imtophat(I,se7);
se8=strel('disk',9);
t8=imtophat(I,se8);
se9=strel('disk',10);
t9=imtophat(I,se9);
se10=strel('disk',11);
t10=imtophat(I,se10);
se11=strel('disk',12);
t11=imtophat(I,se11);
se12=strel('disk',13);
t12=imtophat(I,se12);
se13=strel('disk',14);
t13=imtophat(I,se13);
se14=strel('disk',16);
t14=imtophat(I,se14);

[L0,num0]=bwlabel(t0);
[L1,num1]=bwlabel(t1);
%figure,imshow(label2rgb(L1,@jet,'k'))
[L2,num2]=bwlabel(t2);
%figure,imshow(label2rgb(L2,@jet,'k'))
[L3,num3]=bwlabel(t3);
[L4,num4]=bwlabel(t4);
[L5,num5]=bwlabel(t5);
[L6,num6]=bwlabel(t6);
[L7,num7]=bwlabel(t7);
[L8,num8]=bwlabel(t8);
[L9,num9]=bwlabel(t9);
[L10,num10]=bwlabel(t10);
[L11,num11]=bwlabel(t11);
[L12,num12]=bwlabel(t12);
[L13,num13]=bwlabel(t13);
[L14,num14]=bwlabel(t14);
%figure,imshow(label2rgb(L3,@jet,'k'))
figure,plot([0 2 3 4 5 6 7 8 9 10 11 12 13 14 16],[num0 num1 num2 num3 num4 num5 num6 num7 num8 num9 num10 num11 num12 num13 num14])

se0=strel('disk',1);
t0=imbothat(I,se0);
se1=strel('disk',2);
t1=imbothat(I,se1);
se2=strel('disk',3);
t2=imbothat(I,se2);
se3=strel('disk',4);
t3=imbothat(I,se3);
se4=strel('disk',5);
t4=imbothat(I,se4);
se5=strel('disk',6);
t5=imbothat(I,se5);
se6=strel('disk',7);
t6=imbothat(I,se6);
se7=strel('disk',8);
t7=imbothat(I,se7);
se8=strel('disk',9);
t8=imbothat(I,se8);
se9=strel('disk',10);
t9=imbothat(I,se9);
se10=strel('disk',11);
t10=imbothat(I,se10);
se11=strel('disk',12);
t11=imbothat(I,se11);
se12=strel('disk',13);
t12=imbothat(I,se12);
se13=strel('disk',14);
t13=imbothat(I,se13);
se14=strel('disk',16);
t14=imbothat(I,se14);

[L0,num0]=bwlabel(t0);
[L1,num1]=bwlabel(t1);
%figure,imshow(label2rgb(L1,@jet,'k'))
[L2,num2]=bwlabel(t2);
%figure,imshow(label2rgb(L2,@jet,'k'))
[L3,num3]=bwlabel(t3);
[L4,num4]=bwlabel(t4);
[L5,num5]=bwlabel(t5);
[L6,num6]=bwlabel(t6);
[L7,num7]=bwlabel(t7);
[L8,num8]=bwlabel(t8);
[L9,num9]=bwlabel(t9);
[L10,num10]=bwlabel(t10);
[L11,num11]=bwlabel(t11);
[L12,num12]=bwlabel(t12);
[L13,num13]=bwlabel(t13);
[L14,num14]=bwlabel(t14);
hold on
plot([0 2 3 4 5 6 7 8 9 10 11 12 13 14 16],[num0 num1 num2 num3 num4 num5 num6 num7 num8 num9 num10 num11 num12 num13 num14],'r')