mayidePCNN.rar_PCNN图像_pcnn_pcnn马德_马义德

function[Edge,Numberofaera]=PCNN(X) %X:输入的灰度图像,Edge:检测到的一些边界点,Numberofaera则表明了在各次迭代时激活的块区域 clear; clc; I=imread('lena1.bmp'); [Xa,Ya]=size(I); subplot(1,2,1); imshow(I); %imshow(I); X=double(I); %X=double(imread('lena.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本来是激发的像素,现在边界被置零, %也不能说白激发了,B矩阵有纪录!当然下次就休想再被激发 [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(1,2,2); imshow(S); end %while