With the rapid advent of satellite imaging, remote sensing images are now widely used for object classification and recognition [1]}, tracking [2]}, environmental monitoring [3]}, etc. Multiband imaging sensors acquire information as a band of two-dimensional images, where each image captures a narrow wavelength of the electromagnetic spectrum. The precision of these sensors is determined by their spectral and spatial resolutions, i.e., how densely the spectral and spatial domains are sampled. However, existing sensors have a fundamental tradeoff between spatial and spectral resolutions [1]}, [5]}. For example, hyperspectral (HS) images have hundreds of bands but low spatial resolution, while multispectral (MS) images have fewer bands but higher spatial resolution. An extreme case is a panchromatic (PAN) image which has a single band, but its spatial resolution is much higher than HS and MS images. Reconstruction of high-resolution multiband images by fusing images with complementary spatial and spectral resolutions is thus of great interest [1]}.
