In the past years, one common way of enhancing the spatial resolution of a hyperspectral (HS) image has been to fuse it with complementary information coming from multispectral (MS) or panchromatic images. This paper proposes a new method for reconstructing a high-spatial, high-spectral image from measurements acquired after compressed sensing by multiple sensors of different spectral ranges and spatial resolutions, with specific attention to HS and MS compressed images. To solve this problem, we introduce a fusion model based on the linear spectral unmixing model classically used for HS images and investigate an optimization algorithm based on a block coordinate descent strategy. The nonnegative and sum-to-one constraints resulting from the intrinsic physical properties of abundances as well as a total variation penalization are used to regularize this ill-posed inverse problem. Simulation results conducted on realistic compressed HS and MS images show that the proposed algorithm can provide fusion results that are very close to those obtained with uncompressed images, with the advantage of using a significantly reduced number of measurements.