demo_multi_dy.m

% ========================================================================
% Detail-Preserving Multi-Exposure Fusion with Edge-Preserving Structural
% Patch Decomposition, IEEE TCSVT,2021
% algorithm Version 1.0
% Copyright(c) 2021, Hui Li, Tsz Nam Chan, Xianbiao Qi and Wuyuan Xie
% All Rights Reserved.
% ----------------------------------------------------------------------
% Permission to use, copy, or modify this software and its documentation
% for educational and research purposes only and without fee is hereby
% granted, provided that this copyright notice and the original authors'
% names appear on all copies and supporting documentation. This program
% shall not be used, rewritten, or adapted as the basis of a commercial
% software or hardware product without first obtaining permission of the
% authors. The authors make no representations about the suitability of
% this software for any purpose. It is provided "as is" without express
% or implied warranty.
%----------------------------------------------------------------------
% This is a dynamic scene implementation of "Detail-Preserving Multi-Exposure Fusion with Edge-Preserving Structural Patch Decomposition"
% Please refer to the following paper:
% H. Li et al., "Detail-Preserving Multi-Exposure Fusion with Edge-Preserving Structural Patch Decomposition, 2021" In press
% IEEE Transactions on Circuits and Systems for Video Technology
% Please kindly report any suggestions or corrections to xiaohui102788@126.com
%----------------------------------------------------------------------
clear all;
close all;
addpath(genpath(pwd));

folderNames = {'ArchSequence'};

for i = 1:size(folderNames,2)
% for i=9
    folderName = folderNames{i};
    
 Dir =sprintf('%s/',folderName);
    
    %       Dir =sprintf(' E:/20150816????????????/TIP_MSPD/TIP_MSPDcode/image_sequence(downsampled)/%s/',folderName);  %static scenes
    imgSeqColor= loadImg(Dir); % [0,1]
    imgSeqColor = downSample(imgSeqColor, 1024);
         %% motion detection      
r=10;
  
tic
  imgSeqColor = detection_color(imgSeqColor,r);
    r1=4;
    

    T=0.5;
    lambda=0.25;
    %% single scale
    % C_out= SPD_fast_single3(imgSeqColor,r,t,eps);
    %% multi-scale scale
    [ D1,i_mean1,aa1,N1] = scale_fine(imgSeqColor,r1,lambda);
    

    [w,h,~,~]=size(imgSeqColor);
nlev = floor(log(min(w,h)) / log(2))-5;

D2 = cell(nlev,1);
aa2= cell(nlev,1);
N2= cell(nlev,1);
    r2=4;


     lambda=lambda*T;

for ii=1:nlev
    [ D2{ii},i_mean2,aa2{ii},N2{ii}] = scale_interm(i_mean1,r2,lambda);
    i_mean1=i_mean2;
    lambda=lambda*T;

end

   
%% the coarsest  scale
    r3=3;
    t=1;

[fI3,i_mean3,aa3,N3] = scale_coarse(i_mean2,r3,lambda);

%% reconstruct
%% Intermediate layers
for ii=nlev:-1:1
    temp=aa2{ii};
    fI=zeros(size(temp));
    fI(1:2:size(temp,1),1:2:size(temp,2))=fI3;
    B2=boxfilter(fI, r2)./ N2{ii}+D2{ii};
    
    fI3=B2;
end
%% finest layers
fI=zeros(size(aa1));
fI(1:2:size(aa1,1),1:2:size(aa1,2))=B2;
B1=boxfilter(fI, r1)./ N1;
% C_out=repmat(B1,[1 1 3])+2/pi*atan(2.3*D1);
C_out=repmat(B1,[1 1 3])+D1;
toc
    
%     % scale=1;
%     % for ii=1:scale
%     % [ fI2_detail(:,:,ii),i_mean2] = scale2(i_mean1,r,t);
%     % i_mean1=i_mean2;
%     % end
%     % fI2_detail_sum=sum(fI2_detail,3);
%     
%     
%     r3=5;
%     t=1;
%     scale=0.08;
%     [fI3,i_mean3,aa3,N3] = scale3(i_mean2,r3,t,scale);
%     
%     %% reconstruct
%     
%     fI=zeros(size(aa2));
%     fI(1:2:size(aa2,1),1:2:size(aa2,2))=fI3;
%     B2=boxfilter(fI, r2)./ N2+D2;
%     
%     fI=zeros(size(aa1));
%     fI(1:2:size(aa1,1),1:2:size(aa1,2))=B2;
%     B1=boxfilter(fI, r1)./ N1;
%     C_out=repmat(B1,[1 1 3])+D1;
    
    % C_out=repmat(fI2,[1 1 3])+repmat(fI2_detail_sum,[1 1 3])+fI1_detail;
    
    % C_out=repmat(fI2,[1 1 3])+repmat(fI2_detail_sum,[1 1 3])+fI1_detail;
    
    %  toc
    figure,imshow(C_out,[])
    
    
    
    %% estimating images
%     filename = sprintf('%s',folderName);
%     %  imwrite(C_out,strcat('./exr_result',filename,'.png'));
%     filename = ['/Users/huili/Desktop/Multi_mean_ifv/comparison methods/ADF/res_adf_test/',filepath1(i).name(1:end-4),'_adf',filepath1(i).name(end-5:end)];
%     imwrite(C_out,strcat(filename,'.png'));
    
    
     filename = [folderName,'_espd','.png'];
        %  imwrite(C_out,strcat('./exr_result',filename,'.png'));
     imwrite(C_out,filename);

   
    
end