PCNN.zip_PCNN 图像分割_PCNN分割程序_PCNN图像分割_pcnn 分割_pcnn分割

function[Edge,Numberofaera]=PCNN(X) %X:输入的是彩色图像并转换为灰度图像,Edge:检测到的一些边界点,Numberofaera则表明了在各次迭代时激活的块区域 clear; clc; A=imread('flower.jpg'); subplot(2,2,1) imshow(A) title('原始图像') I=rgb2gray(A); [Xa,Ya]=size(I); subplot(2,2,2); imshow(I); title('灰度图像') %imshow(I); X=double(I); %X=double(imread('flower.bmp')); Weight=[0.5 1 0.5;1 0 1;0.5 1 0.5]; Beta=0.32; Yuzhi=200; Decay=0.31; [a,b]=size(X); Threshold=zeros(a,b); S=zeros(a+2,b+2); B=zeros(a,b); %标记样板,表明该pixel是否被激活过; Y=zeros(a,b); Edge=zeros(a,b);Numberofaera=zeros(a,b);Numberofaera_1=zeros(a,b); Num_1=0;Num=0; n=1; while(sum(sum(B))~=Xa*Ya) %若采用128*128的图像，须注意。 for i0=2:a+1 for i1=2:b+1 V=[S(i0-1,i1-1) S(i0-1,i1) S(i0-1,i1+1); S(i0,i1-1) S(i0,i1) S(i0,i1+1); S(i0+1,i1-1) S(i0+1,i1) S(i0+1,i1+1)]; L=sum(sum(V.*Weight)); F=X(i0-1,i1-1); U=double(F)*(1+Beta*double(L)); if U>=Threshold(i0-1,i1-1)|Threshold(i0-1,i1-1)<95 T(i0-1,i1-1)=1; Threshold(i0-1,i1-1)=Yuzhi; Y(i0-1,i1-1)=1; if n==1 B(i0-1,i1-1)=0; %避免第一次全部激发造成的影响 else B(i0-1,i1-1)=1; %已发射过的标记 Threshold(i0-1,i1-1)=1000000;%相当于不会被第二次激活 end else T(i0-1,i1-1)=0; %no use? Y(i0-1,i1-1)=0; end end end Threshold(find(B~=1))=exp(-Decay)*Threshold(find(B~=1)); %被激活过的像不再参与迭代过程 if n~=1 Edge=Edge+judge_edge(Y,n); Y(find(Edge<0))=0; %边界点被置零,Y本来是激发的像素,现在边界被置零, [Numberofaera_1,Num_1]=bwlabel(Y,4); for i=1:a for j=1:b if Numberofaera_1(i,j)~=0 Numberofaera_1(i,j)=Numberofaera_1(i,j)+Num; end end end Numberofaera=Numberofaera+Numberofaera_1; Num=Num_1; end if n==1 S=zeros(a+2,b+2); else S=Bianhuan(T); end n=n+1; Numberofaera_1=zeros(a,b); subplot(2,2,3); imshow(S); title('分割后图像') end %while