%% by Vincent  path(path,'PCNN_toolbox/') path(path,'FusionEvaluation/') path(path,'fusetool/') %% %Low_Coeffs_Rule='PCNN'; %'ave','PCNN','SF-PCNN' %High_Coeffs_Rule='PCNN'; %'max','PCNN','SF-PCNN' % %% Parameters for Curvelet im1=double(imread('lena_1.png'));  im2=double(imread('lena_2.png'));  im3=double(imread('lena_3.png'));  im4=double(imread('lena_4.png'));  im1=im1(:,:,1);  im2=im2(:,:,1);  Im3=im3(:,:,1);  Im4=im4(:,:,1);  %% Parameters for PCNN Para.iterTimes=200; Para.link_arrange=3; Para.alpha_L=0.06931;% 0.06931 Or 1 Para.alpha_Theta=0.2; Para.beta=3;% 0.2 or 3 Para.vL=1.0; Para.vTheta=20; %% disp('Decompose the image via curvelet ...') yA = fdct_wrapping(im1,1,2);  yB = fdct_wrapping(im2,1,2); yC = fdct_wrapping(im1,1,2);  yD = fdct_wrapping(im2,1,2);  n = length(yA);  %% Initialized the coefficients of fused image Fused=yA;  for l = 1:n  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%     for d = 1:length(yA{l})         %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%         Ahigh = yA{l}{d};         Bhigh = yB{l}{d};         Chigh = yC{l}{d};         Dhigh = yD{l}{d};        [Abeta,Bbeta] = weighting(yA{l}{d},yB{l}{d});        [Cbeta,Dbeta] = weighting(yC{l}{d},yD{l}{d});        Fused{l}{d} = Four_PCNN(Ahigh, Bhigh,Chigh,Dhigh,Abeta,Bbeta,Cbeta,Dbeta);               %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%     end end %============================================= f_im = real(ifdct_wrapping(Fused,1));  disp('Reconstruct is ended...') %% Fusion=Fusion*255; Fusion(Fusion<0)=0; %disp('F>255') Fusion(Fusion>255)=255; Fusion=round(Fusion); %% Eval=Evaluation(double(ori_A),double(ori_B),Fusion,256); disp(Eval) %%