图像拉普拉斯金字塔融合算法源代码

function Y = fuse_lap（M1M2ztflag）
%UNtitleD Summary of this function goes here
%   Detailed explanation goes here

[z1s1]=size（M1）;
TM1 = double（M1）;

[z2s2]=size（M2）;
TM2=double（M2）;

if（z1~=z2）||（s1~=s2）
error（‘Input image are not of the same size!‘）;
return;
end

w=1/256*[1 4 6 4 1;
        4 16 24 16 4;
        6 24 36 24 6;
        4 16 24 16 4;
        1 4 6 4 1];
 w=double（w）;

E1=cell（1zt）;
E2=cell（1zt）;

for i1=1:zt-1
    [zs]=size（TM1）;
    zl（i1）=z;
    sl（i1）=s;

    G1=conv2（TM1w‘same‘）;
    G2=conv2（TM2w‘same‘）;
    
    M1T = conv2（（undec2（dec2（G1）））（4*w）‘same‘）;
	M2T = conv2（（undec2（dec2（G2）））（4*w）‘same‘）;

    [m_t1n_t1]=size（TM1）;
    [m_t2n_t2]=size（M1T）;
    m_t1= max（m_t1m_t2）;
    n_t1 = max（n_t1n_t2）;   
    TM1（m_t1n_t1）=0;
    M1T（m_t1n_t1）=0;
    
    [m_t3n_t3]=size（M2T）;
    [m_t4n_t4]=size（TM2）;
    m_t3= max（m_t3m_t4）;
    n_t3 = max（n_t3n_t4）;   
    M2T（m_t3n_t3）=0;
    TM2（m_t3n_t3）=0;

   
    E1（i1） = {TM1 - M1T};
    E2（i1） = {TM2 - M2T};
    
    TM1 = dec2（G1）;
    TM2 = dec2（G2）;   

end
%循环重构
%    [pq]=s