利用matlab编写的区域生长的程序，用于图像分割_20130520213934.rar

clear all clc I = imread('C:\Documents and Settings\Administrator\桌面\新的极实验\1_256.tif'); % I = imread('lena_B.jpg'); [m,n,l]=size(I); if l==3 I=rgb2gray(I); end figure(1), imshow(I,[]) title('初始图像'); figure(2), imhist(I) title('初始图像的直方图'); T=30;%可以修改的阈值 f=double(I); [m,n]=size(I); shed1=zeros(3,round(m*n/2)); %存储区域生长方向上的点和该区域的均值的绝对差值和该点的坐标 sp1=0; % sp1 相当于指针，指向 shed1 中的最后放入的值和坐标 shed2=zeros(2,m); %存储符合生长条件的点的坐标 sp2=0; % sp2 指针，指向 shed2 中的最后放入点的坐标 Cut=ones(size(f))*260; % Cut 为区域生长后的新图像 %%Cut=ones(m,n)*260; hmin=min(min(f)); hmax=max(max(f)); % if hmax>253 % hmax=253; % end h1 = waitbar(0,'Please wait...'); for q=hmin:hmax for i=1:m for j=1:n if (f(i,j)==q&Cut(i,j)==260) %确定该点满足作为种子的条件，且未并入已% 有生长区域 % Cut(i,j)=259; % 259 时，标记该点在原图像的对应点已并入生长区域 ave=f(i,j); %确定新区域的均值的起始值 k=1; %设置生成的区域的象素个数 [Cut,shed1,sp1]=ruzhan(f,Cut,shed1,sp1,ave,i,j,m,n); %把周围的 8 个点入%栈 [shed1,sp1]=arrange(shed1,sp1); %对栈 shed1 的数据进行由大到小的排序 [shed1,sp1,shed2,sp2]=listed(shed1,sp1,T,shed2,sp2); %% 确定符合条件的生%%长点，将它从 shed1 中取出，并放入 shed2 中 end % 根据生长点开始用 8 连通方式进行生长 while (sp2~=0) %当 sp2=0 时表示找不到符合的点， % if (sp2~=0) %当有新的值加入区域时，求新的平均值 sum=ave*k; for t=1:sp2 x=shed2(1,t); y=shed2(2,t); sum=sum+f(x,y); k=k+1; end ave=sum/k;%%%%合并后的区域的均值 % end for t=1:sp2 %合并栈 shed2 中的点，生成新的区域 x=shed2(1,t); y=shed2(2,t); Cut(x,y)=q; [Cut,shed1,sp1]=ruzhan(f,Cut,shed1,sp1,ave,x,y,m,n); end sp2=0; [shed1,sp1]=arrange(shed1,sp1); [shed1,sp1,shed2,sp2]=listed(shed1,sp1,T,shed2,sp2); end % 在一片区域生成之后，对栈 shed1(1,:)中未并入区域的值进行处理 if (sp1~=0) for t=1:sp1 x=shed1(2,t); y=shed1(3,t); Cut(x,y)=260; end sp1=0; shed1=zeros(3,round(m*n/2)); shed2=zeros(2,m); end end end waitbar(q/(hmax-1)) end max(max(Cut)) % I=find(Cut>=254); % Cut(I)=0; % Cut=im2uint8(Cut); Cut=uint8(Cut); imshow(Cut) BW2 = edge(Cut,'sobel',0.000000000000001); figure(22), imshow(BW2) %%%%Cut为区域生长之后的图像 %%%%%%将Cut中的区域的个数放在Cut1中，查看区域生长之后区域的个数和分布 clear k Cut1=zeros(m,n); Cut2=zeros(m,n); Cut2=im2bw(Cut2); k=1; a=min(min(Cut)); b=max(max(Cut)); a=double(a); b=double(b); for i=a:b II=find(Cut==i); if length(II)~=0 Cut2(II)=1; [L,num]=bwlabel(Cut2,8); for j=1:num III=find(L==j); Cut1(III)=k; k=k+1; end Cut2(II)=0; end end a1=max(max(Cut1)) a2=min(min(Cut1)) %%%%%%%将Cut1中的区域进行合并，将合并之后的区域放在Z4中 a1=double(a1); f=I; f=double(f); Z1=Cut1; Z2=Z1; Z3=Z1; [m,n]=size(Z1); for i=1:a1 [L,H]=find(Z1==i); if length(L)<=200%若区域太小，则合并 j=1; for i1=1:length(L) if (L(i1)-1>0)&(H(i1)-1>0) if Z1(L(i1)-1,H(i1)-1)~=i zhan(1,j)=Z1(L(i1)-1,H(i1)-1); j=j+1; end end if L(i1)-1>0 if Z1(L(i1)-1,H(i1))~=i zhan(1,j)=Z1(L(i1)-1,H(i1)); j=j+1; end end if (L(i1)-1>0)&(H(i1)+1<=n) if Z1(L(i1)-1,H(i1)+1)~=i zhan(1,j)=Z1(L(i1)-1,H(i1)+1); j=j+1; end end if (H(i1)-1>0) if Z1(L(i1),H(i1)-1)~=i zhan(1,j)=Z1(L(i1),H(i1)-1); j=j+1; end end if (H(i1)+1<=n) if Z1(L(i1),H(i1)+1)~=i zhan(1,j)=Z1(L(i1),H(i1)+1); j=j+1; end end if (H(i1)+1<=n&L(i1)+1<=m) if Z1(L(i1)+1,H(i1)+1)~=i zhan(1,j)=Z1(L(i1)+1,H(i1)+1); j=j+1; end end if (L(i1)+1<=m) if Z1(L(i1)+1,H(i1))~=i zhan(1,j)=Z1(L(i1)+1,H(i1)); j=j+1; end end if (L(i1)+1<=m&H(i1)-1>0) if Z1(L(i1)+1,H(i1)-1)~=i zhan(1,j)=Z1(L(i1)+1,H(i1)-1); j=j+1; end end end %%%%%%%%%%%%%%将和这些点相连的那些点放入矩阵zhan中 if ~isempty(zhan) y(1,1)=zhan(1,1); j1=2; for i3=2:length(zhan) for i4=1:length(y) if zhan(1,i3)~=y(1,i4) k=1; break else k=0; end end if k==1 y(1,j1)=zhan(1,i3); j1=j1+1; end end end %%%%%%以上将zhan里面的元素不同的元素放在y矩阵中 %%%以下将当前区域合并到与之相邻的区域中，规则为平均灰度值最接近,合并后的矩阵为Z2； if length(y)==1 for i5=1:length(L) Z2(L(i5),H(i5))=y(1,1); end else s1=0; for i6=1:length(L) s1=s1+f(L(i6),H(i6));%%%%为要处理的像素的总灰度值 end s=0; II=find(Z1==y(1,1)); for i7=1:length(II) s=s+f(II(i7)); %%%%%为第一个总灰度值 end cha=abs(s-s1); log=y(1,1); for i8=2:length(y) II=find(Z1==y(1,i8)); s2=0; for i9=1:length(II) s2=s2+f(II(i9)); end cha1=abs(s1-s2); if cha1<cha cha=cha1; log=y(1,i8); end end for i10=1:length(L) Z2(L(i10),H(i10))=log;%%%%合成之后的为Z2矩阵 end end end Z1=Z2; zhan=[]; y=[]; k=0; end Z4=zeros(m,n); k=1; m1=min(min(Z1)); m2=max(max(Z1)); for i11=m1:m2 II=find(Z1==i11); if ~isempty(II) Z4(II)=k; k=k+1; end end max(max(Z4)) %%%%%分成的区域的个数 min(min(Z4)) %%%%a1为Cut中的区域个数， BW2 = edge(Z4,'sobel',0.000000000000001); figure(30),imshow(BW2); Lrgb = label2rgb(Z4); figure(33), imshow(Lrgb), title('Lrgb') imwrite(Lrgb,'C:\Documents and Settings\Administrator\桌面\lena30_200.tif') III=find(BW2==1); I(III)=255; figure(32),imshow(I) imwrite(I,'C:\Documents and Settings\Administrator\桌面\lena_I.tif') % im1=imread('C:\Documents and Settings\Administrator\桌面\2.tif'); im2=xiaos_ipolartrans(Z4, 256, 256, 128, 128,'log','valid'); figure(44),imshow(im2,[]) X44=mat2gray(im2) imwrite(X44,'C:\Documents and Settings\Administrator\桌面\lena_huanyuan.tif')