Optimal transforms for multispectral and multilayer image coding

Multispectral images are composed of a series of images at differing optical wavelengths. Since these images can be quite large, they invite efficient source coding schemes for reducing storage and transmission requirements. Because multispectral images include a third (spectral) dimension with nonstationary behavior, these multilayer data sets require specialized coding techniques. The authors develop both a theory and specific methods for performing optimal transform coding of multispectral images. The theory is based on the assumption that a multispectral image may be modeled as a set of jointly stationary Gaussian random processes. Therefore, the methods may be applied to any multilayer data set which meets this assumption. Although the authors do not assume the autocorrelation has a separable form, they show that the optimal transform for coding has a partially separable structure. In particular, they prove that a coding scheme consisting of a frequency transform within each layer followed by a separate KL transform across the layers at each spatial frequency is asymptotically optimal as the block size becomes large. Two simplifications of this method are also shown to be asymptotically optimal if the data can be assumed to satisfy additional constraints. The proposed coding techniques are then implemented using subband filtering methods, and the various algorithms are tested on multispectral images to determine their relative performance characteristics.