图像分割基于matlab纹理度量图像分割【含Matlab源码 2329期】.zip

function exp3a global blocksize; global halfsize; global nblockx; global nblocky; blocksize=16; % window size halfsize=blocksize/2; img=imread('test1.bmp'); [height,width]=size(img); nblockx=width*2/blocksize-1; nblocky=height*2/blocksize-1; [MENT,EENT,MASM,EASM,MCON,ECON]=texture_feature(img); % segmentation by kmean delta1=0.5; delta2=1.5; A=mean(MENT); a=mean(A); B=mean(EENT); b=mean(B); C=mean(MASM); c=mean(C); D=mean(EASM); d=mean(D); E=mean(MCON); e=mean(E); F=mean(ECON); f=mean(F); z01=[a*delta1,b*delta1,c*delta1,d*delta1,e*delta1,f*delta1]; z02=[a*delta2,b*delta2,c*delta2,d*delta2,e*delta2,f*delta2]; for i=1:nblocky for j=1:nblockx c(i,j)=0; Z=[MENT(i,j),EENT(i,j),MASM(i,j),EASM(i,j),MCON(i,j),ECON(i,j)]; Zz01=0; Zz02=0; for k=1:6 Zz01=Zz01+(Z(k)-z01(k))*(Z(k)-z01(k)); Zz02=Zz02+(Z(k)-z02(k))*(Z(k)-z02(k)); end if (Zz01>Zz02) c(i,j)=1; end end end m=0; n=0; z11=zeros(1,6); z12=zeros(1,6); for i=1:nblocky for j=1:nblockx Z=[MENT(i,j),EENT(i,j),MASM(i,j),EASM(i,j),MCON(i,j),ECON(i,j)]; if (c(i,j)==0) z11=z11+Z; m=m+1; else z12=z12+Z; n=n+1; end end end z11=z11/m; z12=z12/n; while (z01~=z11 | z02~=z12) z01=z11; z02=z12; for i=1:nblocky for j=1:nblockx c(i,j)=0; Z=[MENT(i,j),EENT(i,j),MASM(i,j),EASM(i,j),MCON(i,j),ECON(i,j)]; Zz01=0; Zz02=0; for k=1:6 Zz01=Zz01+(Z(k)-z01(k))*(Z(k)-z01(k)); Zz02=Zz02+(Z(k)-z02(k))*(Z(k)-z02(k)); end if (Zz01>Zz02) c(i,j)=1; end end end m=0; n=0; z11=zeros(1,6); z12=zeros(1,6); for i=1:nblocky for j=1:nblockx Z=[MENT(i,j),EENT(i,j),MASM(i,j),EASM(i,j),MCON(i,j),ECON(i,j)]; if (c(i,j)==0) z11=z11+Z; m=m+1; else z12=z12+Z; n=n+1; end end end z11=z11/m; z12=z12/n; end % generate the segmentation resulting image S=zeros(height,width); for i=1:nblocky for j=1:nblockx if (c(i,j)==0) for m=1:blocksize for n=1:blocksize S(halfsize*i-halfsize+m,halfsize*j-halfsize+n)=255; end end else for m=1:blocksize for n=1:blocksize S(halfsize*i-halfsize+m,halfsize*j-halfsize+n)=0; end end end end end H=histeq(img); figure,subplot(1,3,1),imshow(img); title('原始图片'); subplot(1,3,2),imshow(H); title('图片增强'); subplot(1,3,3),imshow(S); title('分割结果');