PCNN去噪函数

function I=denoise_gj(X,c)%PCNN检测噪声并去除（子函数）if nargout == 1 c=2;end[w,h]=size(X);Weight =[0.707 1 0.707;1 0 1;0.707 1 0.707];beta=0.1;alpha=0.1;threshold0=280;threshold=zeros(w,h);Y = zeros(w+1,h+1);for c1=c:-1:1 for i=2:w-1 for j=2:h-1 F=X(i,j); V=[Y(i-1,j-1) Y(i-1,j) Y(i-1,j+1); Y(i,j-1) Y(i,j) Y(i,j+1); Y(i+1,j-1) Y(i+1,j) Y(i+1,j+1)]; L=sum(sum(V.*Weight)); U=double(F)*(1+beta*L); if U>threshold(i,j) %点火产生脉冲 Y(i,j)=1; if c1~=2 M=medfilt2(X(i-1:i+1,j-1:j+1)); X(i,j)=M(5); end else Y(i,j)=0; end threshold(i,j)=exp(-alpha)*threshold(i,j)+threshold0*Y(i,j); end endendI=X;end