【谷速代碼】matlab源碼 nsct與pcnn結合的圖像融合算法

1、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% P

2、ara.alpha_Theta=0.2;% Para.beta=3;% 0.2 or 3% Para.vL=1.0;% Para.vTheta=20;%=============================================================disp('PCNN is processing...')%=============================================

3、================% 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;%==============================

4、===============================[p,q]=size(matrix);F_NA=abs(matrix);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Initialize the parameters. % You'd better change them according to your applications % alpha

5、_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

6、_arrange/2);center_y=round(link_arrange/2);High_Coeffs_Rule='SF-PCNN'; %'max','PCNN','SF-PCNN'%%ori_A=imread('clockA.tif');ori_B=imread('clockB.tif');A=double(ori_A)/255;B=

7、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

8、=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 %%%%%%%%

9、%%%% % save yA yA% save yB yB%%%%%%%%%%%% load the NSCT coefficients %%%%%%%%%%%% load yAload yB%%n = length(yA);%% Initialized the coefficients of fused imageFused=yA;%%%=============================================% Lo