图像处理中的模糊C均值聚类matlab算法及源代码

function [IX2]=fcm(IM); %IM是输入的源图像 %IX2是分类结果 IM=imread('wwwlt163com200806191143.jpg'); IM=IM(:,:,1); figure(1) imshow(uint8(IM)) [maxX,maxY]=size(IM); IM=double(IM); IMM=cat(4,IM,IM,IM,IM); %初始化聚类中心（4类） cc1=7; cc2=70; cc3=180; cc4=240; ttFcm=0; while(ttFcm<15) ttFcm=ttFcm+1; c1=repmat(cc1,maxX,maxY); c2=repmat(cc2,maxX,maxY); c3=repmat(cc3,maxX,maxY); c4=repmat(cc4,maxX,maxY); c=cat(4,c1,c2,c3,c4); ree=repmat(0.000001,maxX,maxY); ree1=cat(4,ree,ree,ree,ree); distance=IMM-c; distance=distance.*distance+ree1; daoshu=1./distance; daoshu2=daoshu(:,:,1)+daoshu(:,:,2)+daoshu(:,:,3)+daoshu(:,:,4); %计算隶属度u distance1=distance(:,:,1).*daoshu2; u1=1./distance1; distance2=distance(:,:,2).*daoshu2; u2=1./distance2; distance3=distance(:,:,3).*daoshu2; u3=1./distance3; distance4=distance(:,:,4).*daoshu2; u4=1./distance4; %计算聚类中心z ccc1=sum(sum(u1.*u1.*IM))/sum(sum(u1.*u1)); ccc2=sum(sum(u2.*u2.*IM))/sum(sum(u2.*u2)); ccc3=sum(sum(u3.*u3.*IM))/sum(sum(u3.*u3)); ccc4=sum(sum(u4.*u4.*IM))/sum(sum(u4.*u4)); tmpMatrix=[abs(cc1-ccc1)/cc1,abs(cc2-ccc2)/cc2,abs(cc3-ccc3)/cc3,abs(cc4-ccc4)/cc4]; pp=cat(3,u1,u2,u3,u4); for i=1:maxX for j=1:maxY if max(pp(i,j,:))==u1(i,j) IX(i,j)=1; elseif max(pp(i,j,:))==u2(i,j) IX2(i,j)=2; elseif max(pp(i,j,:))==u3(i,j) IX2(i,j)=3; else IX2(i,j)=4; end end end %判结束条件 if max(tmpMatrix)<0.0001 break; else cc1=ccc1; cc2=ccc2; cc3=ccc3; cc4=ccc4; end for i=1:maxX for j=1:maxY if IX2(i,j)==1 IMMM(i,j)=240; elseif IX2(i,j)==3 IMMM(i,j)=180; elseif IX2(i,j)==2 IMMM(i,j)=70; else IMMM(i,j)=7; end end end %显示每次聚类分割结果 figure(2); imshow(uint8(IMMM)); end for i=1:maxX for j=1:maxY if IX2(i,j)==4 IMMM(i,j)=240; elseif IX2(i,j)==3 IMMM(i,j)=180; elseif IX2(i,j)==2 IMMM(i,j)=70; else IMMM(i,j)=7; end end end %显示最终分类结果 IMMM=uint8(IMMM); figure(3); imshow(IMMM); end