fun_pcnn.rar_Fun_ Fun_ Fun_pcnn

% function [Fea]=fun_pcnn(I); function B=fun_pcnn(I,Bate,V_E,Alpha_E,W); % pcnn processing %I = imread('D:\毕业设计\艾永好课题相关材料\艾永好课题相关材料\图像样本\aluminium（铝板带）\凹坑\13.bmp') [r,c] = size(I); B=zeros(r,c); %Y=zeros(r,c); S=zeros(r+2,c+2); M=zeros(r,c); E=zeros(r,c); % T=zeros(r+2,c+2); % T(2:r+1,2:c+1)=I; % J=I; Bate=0.1; V_L=1; % V_E=1.1; % V_E=20; V_E=256; %Alpha_E=0.05; Alpha_E=0.3; % W=[0.5 1 0.5;1 0 1;0.5 1 0.5]; W=[0.7071 1 0.7071;1 0 1;0.7071 1 0.7071]; n=1; while sum(sum(B))==0 for i=2:r+1 for j=2:c+1 V=[S(i-1,j-1) S(i-1,j) S(i-1,j+1); S(i,j-1) S(i,j) S(i,j+1); S(i+1,j-1) S(i+1,j) S(i+1,j+1)]; L=V_L*sum(sum(V.*W)); F=I(i-1,j-1); U=double(F)*(1+Bate*double(L)); if U>=E(i-1,j-1) %|E(i-1,j-1)<2 M(i-1,j-1)=1; E(i-1,j-1)=V_E; %Y(i-1,j-1)=1; if n==1 B(i-1,j-1)=0; else B(i-1,j-1)=1; E(i-1,j-1)=100000; %限制已经点火的神经元不再点火 end else M(i-1,j-1)=0; % Y(i-1,j-1)=0; end end end E(find(B~=1))=exp(-Alpha_E)*E(find(B~=1));%+VE*sum(sum(V)); S=zeros(r+2,c+2); if n>1 S(2:r+1,2:c+1)=M; end figure,imshow(mat2gray(M)) figure,imshow(B) n=n+1; end % [hy,ly]=find(B~=0); % len=length(hy); % for c=1:len % k=medfilt2(T(hy(c):hy(c)+2,ly(c):ly(c)+2)); % J(hy(c),ly(c))=k(5); % end % figure,imshow(J) %B=fun_pcnn(1-J(:,:,1),0.2,20,0.3)+fun_pcnn(J(:,:,1),0.2,20,0.3); %B=fun_pcnn(1-J(:,:,2),0.2,20,0.3)+fun_pcnn(J(:,:,2),0.2,20,0.3); % B=fun_pcnn(1-J(:,:,3),0.2,20,0.3)+fun_pcnn(J(:,:,2),0.2,20,0.3);