为了尽可能地保留全色图像的空间信息和多光谱图像的光谱信息，提出了一种基于改进梯度投影非负矩阵分解和复Contourlet变换的遥感图像融合方法.首先，以多光谱图像的强度分量图像为标准，对全色图像做直方图匹配，得到新的全色图像；然后，利用复Contourlet变换分别分解多光谱图像的强度分量图像和新的全色图像，得到各自对应的低频分量和高频分量；接着，对两幅低频分量图像采用改进梯度投影的非负矩阵分解作为融合规则获取新的低频分量，并对两幅高频分量图像使用系数绝对值较大法获取新的高频分量；最后，通过逆复Contourlet变换和逆色调-饱和度-强度变换获得融合后的图像.大量实验结果表明，与HSI方法、NMF与无下采样Contourlet变换结合的方法以及提升小波变换与HSI结合的方法相比，本文方法获得的融合图像具有更高的空间分辨率和更多的光谱信息.

In order to preserve both spectral and spatial information simultaneously in fused image, a method of remote sensing image fusion based on improved gradient projection Non-negative Matrix Factorization (NMF) and Complex Contourlet Transform (CCT) was proposed. Firstly, the panchromatic image was histogram matched to the intensity component of multispectral image to obtain a new panchromatic image. Then, the intensity component of multispectral image and the new panchromatic image were decomposed by CCT, respectively. Corresponding low frequency images and high frequency images were obtained. Next, the new low frequency image was obtained by fusion of two low frequency images using improved gradient projection NMF, and the new high frequency image was obtained by fusion of two high frequency images with absolute maximum criterion. Finally, the fused image was reconstructed by inverse CCT and inverse Hue-Saturation-Intensity (HSI) transform. A large number of experimental results show that, compared with other existing relative methods such as HSI method, the method combining NMF with nonsubsampled contourlet transform and the method combining lifting wavelet with HSI, the images fused by proposed method have higher spatial resolution and more spectral information.