function Y = fuse_fsd(M1， M2， zt， ap， mp) %Y = fuse_fsd(M1， M2， zt， ap， mp) image fusion with fsd pyramid % % M1 - input image A % M2 - input image B % zt - maximum decomposition level % ap - coefficient selection highpass (see selc.m) % mp - coefficient selection base image (see sel