ＭＡＴＬＡＢ实现ＦＣＭ算法，简单代码实现

先来看看结果好了，这是一幅原图，经过ＦＣＭ算法之后，假设选择３个ｃｅｎｔｅｒ，将这幅彩色图像转化后得到一个只有三种灰度值的图像，还能将每一种灰度单独提取出来，得到对应的只含有一种灰度的三个图。

它的数学原理为　　　　　　　　

这里Ｖｋ取ｖ１，ｖ２，ｖ３三个任意初始值，经过反复迭代得到稳定的ｖ１，ｖ２，ｖ３即可。

原理方面不再赘述。

详细实现代码如下。

一开始是一个调入图像命令，读者可以根据自己需要修改参数。

BB=imread('C:\Users\User\Desktop\san.jpeg');
CC=rgb2gray(BB);


A=rgb2gray(BB);


A=CC;
A=double(A);
length=412;
width=268;
p=zeros(length,width,3);
u=zeros(length,width,3);
v1=A(8,78);
v2=A(11,82);
v3=A(2,86);
i=1;
j=1;
while(i<=length)
while(j<=width)
p(i,j,1)=abs(double(A(i,j))-double(v1));
p(i,j,2)=abs(double(A(i,j))-double(v2));
p(i,j,3)=abs(double(A(i,j))-double(v3));
j=j+1;
end
j=1;
i=i+1;
end
V1=-100;
V2=-100;
V3=-100;
x=0;
while(abs(v1-V1)>1e-5||abs(v2-V2)>1e-5||abs(v3-V3)>1e-5)
V1=v1;
V2=v2;
V3=v3;
i=1;
j=1;
while(i<=length)
while(j<=width)
if(p(i,j,1)-0<1e-5)
 u(i,j,1)=1;
elseif(p(i,j,2)-0<1e-5) 
u(i,j,2)=1;
elseif(p(i,j,3)-0<1e-5)
 u(i,j,3)=1;
else
 u(i,j,1)=1/[(p(i,j,1)/p(i,j,1))^2+(p(i,j,1)/p(i,j,2))^2+(p(i,j,1)/p(i,j,3))^2];
 u(i,j,2)=1/[(p(i,j,2)/p(i,j,1))^2+(p(i,j,2)/p(i,j,2))^2+(p(i,j,2)/p(i,j,3))^2];
 u(i,j,3)=1/[(p(i,j,3)/p(i,j,1))^2+(p(i,j,3)/p(i,j,2))^2+(p(i,j,3)/p(i,j,3))^2];
end
j=j+1;
end
j=1;
i=i+1;
end
i=1;
j=1;
a1=0;
a2=0;
b1=0;
b2=0;
c1=0;
c2=0;
while(i<=length)
while(j<=width)
a1=a1+u(i,j,1)^2*A(i,j);
a2=a2+u(i,j,1)^2;
b1=b1+u(i,j,2)^2*A(i,j);
b2=b2+u(i,j,2)^2;
c1=c1+u(i,j,3)^2*A(i,j);
c2=c2+u(i,j,3)^2;
j=j+1;
end
j=1;
i=i+1;
end
v1=a1/a2;
v2=b1/b2;
v3=c1/c2;
i=1;
j=1;
while(i<=length)
while(j<=width)
p(i,j,1)=abs(double(A(i,j))-double(v1));
p(i,j,2)=abs(double(A(i,j))-double(v2));
p(i,j,3)=abs(double(A(i,j))-double(v3));
j=j+1;
end
j=1;
i=i+1;
end
x=x+1;
end
Uzz=zeros(length,width);
i=1;
j=1;
while(i<=length)
while(j<=width)
max=u(i,j,1);
k=1;
l=k;
while(k<=3)
if(u(i,j,k)>max)
max=u(i,j,k);
l=k;
end
k=k+1;
end
Uzz(i,j)=l;
j=j+1;
end
j=1;
i=i+1;
end
i=1;
j=1;
while(i<=length)
while(j<=width)
if(Uzz(i,j)==1)
A(i,j)=v1;
elseif(Uzz(i,j)==2)
A(i,j)=v2;
else
A(i,j)=v3;
end
j=j+1;
end
j=1;
i=i+1;
end
A=uint8(A);
figure(1);

imshow(A);／／到这里得到ＦＣＭ模糊结果图像

cluster1=zeros(length,width);
cluster2=zeros(length,width);
cluster3=zeros(length,width);
i=1;
j=1;
while(i<=length)
while(j<=width)
if(double(A(i,j))==6)
cluster1(i,j)=v1;
elseif(double(A(i,j))==66)
cluster2(i,j)=v2;
else
cluster3(i,j)=v3;
end
j=j+1;
end
j=1;
i=i+1;
end


figure(2);
imshow(cluster1)
figure(3);
imshow(cluster2)
figure(4);
imshow(cluster3);／／到这里得到对应只含一个模糊中心的三个不同图像。