Abstract: To deal with the problem of panchromatic and multispectral remote sensing images fusion,a combinatorial optimization image fusion method(COFM) is proposed in this paper.First,the HSI transform is used to obtain the intensity component of a multispectral image;Second,Contourlet transform is adopted to decompose the panchromatic image and the intensity component of the multispectral image,respectively.Consequently the high-and low-frequency of sub-images can be obtained.For high-frequency sub-images,fractal features are extracted and the maximum fusion rule is utilized to fuse these features.For low-frequency sub-images,energy features are extracted,and the non-dominated sorting genetic algorithm-II(NSGA-II) is applied to select the weight coefficients before sub-images are fused with weighted model.The experimental results demonstrate that the performance of COFM is better than that of traditional approaches.COFM can increase the spatial resolution of the fused image while preserving original image's spectral information.