Maximum a Posteriori Fusion Method Based on Gradient Consistency Constraint for Multispectral/Panchromatic Remote Sensing Images
MENG Xiang-chao1, SHEN Huan-feng1*, ZHANG Hong-yan2, ZHAGN Liang-pei2, LI Hui-fang1
1. School of Resource and Environmental Sciences, Wuhan University, Wuhan 430079, China 2. State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
Abstract：Multispectral (MS) images with high spatial resolution (HR) can be obtained by fusing MS images and panchromatic (PAN) image, the HR MS images have an important significance in image interpretation and classification, etc. In the present paper, a new image fusion method based on gradient consistency constraint for MS/PAN images is developed. The method is based on maximum a posteriori (MAP) framework. The relationship of desired HR MS images and PAN image is formulated by gradient consistency constraint. Observation model of MS images and the Huber-Markov priori are combined to solve the fused image by gradient descent algorithm. In the proposed method, gradient consistency constraint is introduced, and defect of band number restriction is overcomed in conventional model-based fusion methods. Iterative step for every band is solved adaptively, and spectral characteristics of each band are fully taken into account, so it not only ensures the spectral information fidelity, but also improves the integration degree of spatial information of fused image. The proposed method has been tested using IKONOS and WorldView-2 images. It is compared with GS, AIHS and AMBF fusion methods from both qualitative and quantitative aspects. Experimental results show that the proposed method can better preserve spectral information while enhance spatial resolution, and it has broader applicability and better fusion result than other methods.
Key words：Image fusion;Multispectral;Panchromatic;Maximum a posteriori (MAP);Gradient consistency
[1] Shah V P, Younan N H, King R L. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(5): 1323. [2] Massip P, Blanc P, Wald L. IEEE Transactions on Geoscience and Remote Sensing, 2012, 50(3): 800. [3] Wang Wencheng, Chang Faliang. Journal of Computers, 2011, 6(12): 2559. [4] Otazu X, Gonzalez-Audicana M, Fors O, et al. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(10): 2376. [5] Da Cunha A L, Zhou J, Do M N. IEEE Transactions on Image Processing, 2006, 15(10): 3089. [6] Jiang Cheng, Zhang Hongyan, Shen Huanfeng, et al. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2013, DOI: 10.1109/JSTARS.2013.2283236. [7] Jiang Cheng, Zhang Hongyan, Shen Huanfeng, et al. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 629. [8] Joshi M V, Bruzzone L, Chaudhuri S. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(9): 2549. [9] Joshi M, Jalobeanu A. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(3): 1245. [10] Molina R, Vega M, Mateos J, et al. Applied and Computational Harmonic Analysis, 2008, 24(2): 251. [11] Li Z, Leung H. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(5): 1482. [12] Zhang Liangpei, Shen Huanfeng, Gong Wei, et al. IEEE Transactions on Systems Man, and Cybernetics, Part B: Cybernetics, 2012, 42(6): 1693. [13] Rahmani S, Strait M, Merkurjev D, et al. IEEE Geoscience and Remote Sensing Letters, 2010, 7(4): 746. [14] Zhang Yun. U. S. Patent, 7340099B2, 2008. [15] Thomas C, Ranchin T, Wald L, et al. IEEE Transactions on Geoscience and Remote Sensing, 2008, 46(5): 1301.
