图像模糊C均值聚类分割matlab代码

亲测可用！图像模糊C均值聚类分割matlab代码，聚类分割后显示图像。 仅需要自己修改读图路径。

IM=imread（‘D:\文件及下载相关\桌面\代码 - 副本\代码\Image\ALLgray\monarchNEW_heibaiedge.bmp‘）; 
 
IM=IM（::1）; 
figure（1） 
imshow（uint8（IM）） 
[maxXmaxY]=size（IM）; 
IM=double（IM）; 
 
IMM=cat（4IMIMIMIM）; 
%??????????????4?? 
cc1=7; 
cc2=70; 
cc3=180; 
cc4=240; 
ttFcm=0; 
while（ttFcm<15） 
    ttFcm=ttFcm+1; 
    c1=repmat（cc1maxXmaxY）; 
    c2=repmat（cc2maxXmaxY）; 
    c3=repmat（cc3maxXmaxY）; 
    c4=repmat（cc4maxXmaxY）; 
    c=cat（4c1c2c3c4）; 
    ree=repmat（0.000001maxXmaxY）; 
    ree1=cat（4reereereeree）; 
    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）/cc1abs（cc2-ccc2）/cc2abs（cc3-ccc3）/cc3abs（cc4-ccc4）/cc4]; 
     
        pp=cat（3u1u2u3u4）; 
        for i=1:maxX 
            for j=1:maxY 
                if max（pp（ij:））==u1（ij） 
                    IX（ij）=1; 
                elseif max（pp（ij:））==u2（ij） 
                    IX2（ij）=2; 
                elseif max（pp（ij:））==u3（ij） 
                    IX2（ij）=3; 
                else 
                    IX2（ij）=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（ij）==1 
                    IMMM（ij）=240; 
                elseif IX2（ij）==3 
                    IMMM（ij）=180; 
                elseif IX2（ij）==2  
                    IMMM（ij）=70; 
                else 
                    IMMM（ij）=7; 
                end 
            end 
        end 
         
        %?????ξ??????? 
        figure（2）; 
        imshow（uint8（IMMM））; 
end 
    for i=1:maxX 
            for j=1:maxY 
                if IX2（ij）==4 
                    IMMM（ij）=240; 
                elseif IX2（ij）==3 
                    IMMM（ij）=180; 
                elseif IX2（ij）==2  
                    IMMM（ij）=70; 
                else 
                    IMMM（ij）=7; 
                end 
            end 
        end 
%???????????? 
IMMM=uint8（IMMM）; 
figure（3）; 
imshow（IMMM）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2938  2019-06-10 19:54  fcm2.m