Multiscale reconstruction for computational spectral imaging

In this work we develop a spectral imaging system and associated reconstruction methods that have been designed to exploit the theory of compressive sensing. Recent work in this emerging field indicates that when the signal of interest is very sparse (i.e. zero-valued at most locations) or highly compressible in some basis, relatively few incoherent observations are necessary to reconstruct the most significant non-zero signal components. Conventionally, spectral imaging systems measure complete data cubes and are subject to performance limiting tradeoffs between spectral and spatial resolution. We achieve single-shot full 3D data cube estimates by using compressed sensing reconstruction methods to process observations collected using an innovative, real-time, dual-disperser spectral imager. The physical system contains a transmissive coding element located between a pair of matched dispersers, so that each pixel measurement is the coded projection of the spectrum in the corresponding spatial location in the spectral data cube. Using a novel multiscale representation of the spectral image data cube, we are able to accurately reconstruct 256×256×15 spectral image cubes using just 256×256 measurements.

[1]  Mark A Neifeld,et al.  Feature-specific imaging. , 2003, Applied optics.

[2]  Alex Zehnder,et al.  The Reuven Ramaty high-energy solar spectroscopic imager (RHESSI) mission , 2003, SPIE Optics + Photonics.

[3]  R. Nowak,et al.  Multiscale likelihood analysis and complexity penalized estimation , 2004, math/0406424.

[4]  N. Sloane,et al.  Hadamard transform optics , 1979 .

[5]  David J. Brady,et al.  High-throughput hyperspectral microscopy , 2006, SPIE BiOS.

[6]  W. Denk,et al.  Two-photon laser scanning fluorescence microscopy. , 1990, Science.

[7]  Xiaobai Sun,et al.  Compressive sampling strategies for integrated microspectrometers , 2006, SPIE Defense + Commercial Sensing.

[8]  Robert D. Nowak,et al.  Signal Reconstruction From Noisy Random Projections , 2006, IEEE Transactions on Information Theory.

[9]  Jonathan Martin Mooney,et al.  High-throughput hyperspectral infrared camera , 1997 .

[10]  Mark Neifeld,et al.  Multispectral principal component imaging. , 2003, Optics express.

[11]  Robert D. Nowak,et al.  Platelets: a multiscale approach for recovering edges and surfaces in photon-limited medical imaging , 2003, IEEE Transactions on Medical Imaging.

[12]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[13]  C. Burrus,et al.  Noise reduction using an undecimated discrete wavelet transform , 1996, IEEE Signal Processing Letters.

[14]  Robert D. Nowak,et al.  Fast multiresolution photon-limited image reconstruction , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[15]  Jane R. Foster,et al.  Application of imaging spectroscopy to mapping canopy nitrogen in the forests of the central Appalachian Mountains using Hyperion and AVIRIS , 2003, IEEE Trans. Geosci. Remote. Sens..

[16]  Xiaobai Sun,et al.  Sensor-layer image compression based on the quantized cosine transform , 2005, SPIE Defense + Commercial Sensing.

[17]  David J. Brady,et al.  Compressive imaging spectrometers using coded apertures , 2006, SPIE Defense + Commercial Sensing.

[18]  N. P. Pitsianis,et al.  Static 2 D aperture coding for multimodal multiplex spectroscopy , 2005 .

[19]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[20]  M. Descour,et al.  Computed-tomography imaging spectrometer: experimental calibration and reconstruction results. , 1995, Applied optics.

[21]  Rainer Riesenberg,et al.  Sensitive Hadamard transform imaging spectrometer with a simple MEMS , 2003, SPIE Remote Sensing.

[22]  Yaakov Tsaig,et al.  Extensions of compressed sensing , 2006, Signal Process..

[23]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[24]  L. Lucy An iterative technique for the rectification of observed distributions , 1974 .

[25]  William L. Smith,et al.  Hyperspectral remote sensing of atmospheric profiles from satellites and aircraft , 2001, SPIE Asia-Pacific Remote Sensing.

[26]  Douglas L. Jones,et al.  Wavelet-based hyperspectral image estimation , 2003, IGARSS 2003. 2003 IEEE International Geoscience and Remote Sensing Symposium. Proceedings (IEEE Cat. No.03CH37477).

[27]  G Oskarsdottir,et al.  Fourier-transform infrared imaging using a rapid-scan spectrometer. , 1999, Optics letters.

[28]  Paul A. Bernhardt,et al.  Direct reconstruction methods for hyperspectral imaging with rotational spectrotomography , 1995 .

[29]  Robert D. Nowak,et al.  Multiscale Poisson Intensity and Density Estimation , 2007, IEEE Transactions on Information Theory.

[30]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[31]  Thomas D. Nielsen,et al.  Hyperspectral imaging: a novel approach for microscopic analysis. , 2001, Cytometry.

[32]  Jonathan Martin Mooney,et al.  Convex projections algorithm for restoration of limited-angle chromotomographic images , 1999 .