In order to enhance the spatial resolution of multispectral images and avoid severe spectral distortion,the amplified multispectral images and the original panchromatic image were decomposed by non-subsampled contourlet transform with the same layers.A 5×5sliding window around the fusion pixel of the matrix in the low frequency sub-band of the decomposed images was selected to compute the local correlation coefficient and fourth-order correlation coefficient.If the calculated local correlation coefficient w...