This article describes a subspace clustering strategy for the spectral compression of multispectral images. Unlike standard principal component analysis, this approach finds clusters in several different subspaces of different dimension. Consequently, instead of representing all spectra in a single low‐dimensional subspace of a fixed dimension, spectral data are assigned to multiple subspaces having a range of dimensions from one to eight. In other words, this strategy allows us to distribute spectra into different subspaces thereby obtaining the best fit for each. As a result, more resources can be allocated to those spectra that need many dimensions for accurate representation and fewer resources to those that can be modeled using fewer dimensions. For a given compression ratio, this trade off reduces the overall reconstruction error. In the case of compressing multispectral images, this initial compression method is followed by JPEG2000 compression in order to remove the spatial redundancy in the data as well. © 2015 Wiley Periodicals, Inc. Col Res Appl, 41, 7–15, 2016