【谷速代码】matlab源码 NSCT与PCNN结合的图像融合算法.docx

function R=PCNN_withParameters(matrix,Para)% Calulate the firing times of neurons in PCNN% Please set the Para as follows:% Parameters for PCNN% Para.iterTimes=200;% Para.link_arrange=3;% Para.alpha_L=1;% 0.06931 Or 1% Para.alpha_Theta=0.2;% Para.beta=3;% 0.2 or 3% Para.vL=1.0;% Para.vTheta=20;%=disp(PCNN is processing.)%=% get the parameters of PCNNlink_arrange=Para.link_arrange;np=Para.iterTimes;alpha_L=Para.alpha_L;alpha_Theta=Para.alpha_Theta ;beta=Para.beta;vL=Para.vL;vTheta=Para.vTheta;%=p,q=size(matrix);F_NA=abs(matrix);% Initialize the parameters. % Youd better change them according to your applications % alpha_L=1;% alpha_Theta=0.2;% beta=3;% vL=1.00;% vTheta=20;% Generate the null matrix that could be usedL=zeros(p,q);U=zeros(p,q);Y=zeros(p,q);Y0=zeros(p,q);Theta=zeros(p,q);% Compute the linking strength.center_x=round(link_arrange/2);center_y=round(link_arrange/2);W=zeros(link_arrange,link_arrange);for i=1:link_arrange for j=1:link_arrange if (i=center_x)&(j=center_y) W(i,j)=0; else W(i,j)=1./sqrt(i-center_x).2+(j-center_y).2); end endend%F=F_NA;for n=1:np K=conv2(Y,W,same); L=exp(-alpha_L)*L+vL*K; Theta=exp(-alpha_Theta)*Theta+vTheta*Y; U=F.*(1+beta*L); Y=im2double(UTheta); Y0=Y0+Y; endR=Y0;%=% If you use the code, please cite the paper as follows:% References:%1 Qu Xiao-Bo, YAN Jing-Wen, XIAO Hong-Zhi, ZHU Zi-Qian. Image Fusion Algorithm Based on Spatial % Frequency-Motivated Pulse Coupled Neural Networks in Nonsubsampled Contourlet Transform Domain. % Acta Automatica Sinica, Vol.34, No.12, pp: 1508-1514.Dec.2008.%=%clc;close all;clear all;%path(path,PCNN_toolbox/)path(path,nsct_toolbox)path(path,FusionEvaluation/)%Low_Coeffs_Rule=SF-PCNN; %ave,PCNN,SF-PCNNHigh_Coeffs_Rule=SF-PCNN; %max,PCNN,SF-PCNN%ori_A=imread(clockA.tif);ori_B=imread(clockB.tif);A=double(ori_A)/255;B=double(ori_B)/255;% Parameters for NSCTpfilt = 9-7;dfilt = pkva;nlevs = 0,1,3,4,4;% default% Parameters for PCNNPara.iterTimes=200;Para.link_arrange=3;Para.alpha_L=1;% 0.06931 Or 1Para.alpha_Theta=0.2;Para.beta=3;% 0.2 or 3Para.vL=1.0;Para.vTheta=20;%disp(Decompose the image via nsct .)% yA=nsctdec(A,nlevs,dfilt,pfilt);% yB=nsctdec(B,nlevs,dfilt,pfilt);% save the NSCT coefficients % % save yA yA% save yB yB% load the NSCT coefficients % load yAload yB%n = length(yA);% Initialized the coefficients of fused imageFused=yA;%=% Lowpass subbanddisp(Process in Lowpass subband.)ALow1= yA1;BLow1 =yB1;ALow2= yA2;BLow2 =yB2;%switch Low_Coeffs_Rule case average Fused1=(ALow1+BLow1)/2; Fused2=(ALow2+BLow2)/2; case PCNN Fused1=fusion_NSCT_PCNN(ALow1,BLow1,Para); Fused2=fusion_NSCT_PCNN(ALow2,BLow2,Para); case SF-PCNN Fused1=fusion_NSCT_SF_PCNN(ALow1,BLow1,Para); Fused2=fusion_NSCT_SF_PCNN(ALow2,BLow2,Para);end%=% Bandpass subbandsdisp(Process in Bandpass subbands.)for l = 3:n % for d = 1:length(yAl) % Ahigh = yAld; Bhigh = yBld; switch High_Coeffs_Rule case max decision_map=(abs(Ahigh)=abs(Bhigh); Fusedld=decision_map.*Ahigh + (decision_map).*Bhigh; case PCNN Fusedld=fusion_NSCT_PCNN(Ahigh,Bhigh,Para); case SF-PCNN Fusedld=fusion_NSCT_SF_PCNN(Ahigh,Bhigh,Para); end % endenddisp(High frequecy field process is ended)%=disp(Reconstruct the image via nsct .)F=nsctrec(Fused, dfilt, pfilt);disp(Reconstruct is ended.)%F=F*255;F(F255)F(F255)=255;F=round(F);%