In remote sensing community, IHS (intensity, hue, and saturation) transform is one of the most commonly used fusion algorithm. A study on IHS fusion indicates that the color distortion cannot be avoided. Meanwhile, wavelet decomposition has a property of frequency division in transform domain. And the statistical property of wavelet coefficient reflects those significant features. So, a united optimal fusion method, which using the statistical property of wavelet decomposition and IHS transform on pixel and feature levels, is proposed. That is, the high frequency of intensity component I is fused on feature level with multi-resolution wavelet in IHS space, and the low frequency of intensity component I is fused on pixel level with optimal weight coefficients. Spectral information and spatial resolution are two performance indexes of optimal weight coefficients. Experiment results show that it is a practical and effective method.