Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

LED-based compressive spectral-temporal imaging

Open Access Open Access

Abstract

A compressive spectral-temporal imaging system is reported. A multi-spectral light-emitting diode array is used for target illumination and spectral modulation, while a digital micro-mirror device (DMD) encodes the spatial and temporal frames. Several encoded video frames are captured in a snapshot of an integrating focal plane array (FPA). A high-frame-rate spectral video is reconstructed from the sequence of compressed measurements captured by the grayscale low-frame-rate camera. The imaging system is optimized through the design of the DMD patterns based on the forward model. Laboratory implementation is conducted to validate the performance of the proposed imaging system. We experimentally demonstrate the video acquisition with eight spectral bands and six temporal frames per FPA snapshot, and thus a 256 × 256 × 8 × 6 4D cube is reconstructed from a single 2D measurement.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

Full Article  |  PDF Article
More Like This
Snapshot compressive imaging based digital image correlation: temporally super-resolved full-resolution deformation measurement

Wenwu Chen, Bo Zhang, Liuning Gu, Haibo Liu, Jinli Suo, and Xinxing Shao
Opt. Express 30(19) 33554-33573 (2022)

Compressive temporal imaging using a rolling shutter camera array

Felipe Guzmán, Pablo Meza, and Esteban Vera
Opt. Express 29(9) 12787-12800 (2021)

Coded aperture compression temporal imaging based on a dual-mask and deep denoiser

YouRan Ge, GangRong Qu, YuHao Huang, and Duo Liu
J. Opt. Soc. Am. A 40(7) 1468-1477 (2023)

Supplementary Material (17)

NameDescription
Visualization 1       Video of the simulation input 4D data-cube with size 256 by 256 by 8 by 32.
Visualization 2       Video of the simulation reconstruction result with 12.5% transmittance random coding.
Visualization 3       Video of the simulation reconstruction result with 50% transmittance random coding.
Visualization 4       Video of the simulation reconstruction result with 12.5% transmittance blue noise coding.
Visualization 5       Video of the simulation reconstruction result with 50% transmittance blue noise coding.
Visualization 6       Video of the simulation reconstruction result with 6 temporal frames compression.
Visualization 7       Video of the simulation reconstruction result with 10 temporal frames compression.
Visualization 8       Video of the simulation reconstruction result with GPSR4D algorithm.
Visualization 9       Video of the simulation reconstruction result with GAPTV4D algorithm.
Visualization 10       Video of the simulation reconstruction result with DeSCI4D algorithm.
Visualization 11       Video of the expeirmental reconstruction result of natural plants with gaptv4D algorithm.
Visualization 12       Video of the expeirmental reconstruction result of color checker with 4 temporal frames compression and DeSCI4D algorithm.
Visualization 13       Video of the expeirmental reconstruction result of color checker with 6 temporal frames compression.
Visualization 14       Video of the expeirmental reconstruction result of Super Mario with 4 temporal frames compression.
Visualization 15       Video of the expeirmental reconstruction result of Super Mario with 6 temporal frames compression.
Visualization 16       Video of the expeirmental reconstruction result of natural plants with 4 temporal frames compression.
Visualization 17       Video of the expeirmental reconstruction result of natural plants with 6 temporal frames compression.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (11)

Fig. 1.
Fig. 1. LeSTI testbed. Spectral LED illumination is used for spectral scene modulation. An imaging lens focuses the spectrally modulated scene onto the DMD imaging plane, where spatial modulations are imposed. A second objective lens focuses the spatial modulated scene onto the imaging plane, where a grayscale detector captures multiple temporally modulated video frames.
Fig. 2.
Fig. 2. (a) LeSTI system architecture. Multi-spectral LEDs response for illumination and spectral modulations; a DMD is used for imposing spatial binary patterns. A grayscale sensor is utilized for image capture. (b) Eight types of LED’s spectra shown in the 400-700nm wavelength range. (c) The designed LED illuminator with multi-spectral LEDs used in our system. Eight different types of LEDs are used and evenly placed on a circuit board with three elements of each type.
Fig. 3.
Fig. 3. (a) Timing sequence in the LeSTI sensing. (b) LeSTI system sensing process (bottom row). A pixel on the scene shows the spectral modulation in the system (middle row).
Fig. 4.
Fig. 4. PSNR and SSIM plots with different system configurations. Notice that the GAP-TV reconstruction algorithm is applied for reconstructions. Here ‘white’ refers to random binary random patterns used in the DMD and ‘Blue’ means the blue noise patterns are used. The following percentage numbers denote the pattern transmittance.
Fig. 5.
Fig. 5. Original image cube (row 1, Visualization 1) and simulation reconstruction results of $N_t = 2$ with random (white noise) coding transmittance 12.5% (row 2, Visualization 2), random coding 50% (row 3, Visualization 3), blue noise coding 12.5% (row 4, Visualization 4), and blue noise coding 50% (row 5, Visualization 5). In the first row, original data is presented with spectral images, reproduced RGB images and zoomed information. Reconstruction results with different configurations (rows 2-5) are shown with the same layout. The compressed measurements are shown at the bottom left (last row). Spectral plots of selected regions are shown in the bottom right.
Fig. 6.
Fig. 6. Simulation reconstruction results of 25% transmittance blue noise with different temporal compression, i.e. $N_t = 6$ (row 2, Visualization 6), $10$ (row 3, Visualization 7). Notice that the P3 interest point of reconstruction $N_t = 6$ is on frame 4, and for reconstruction $N_t = 10$ it is on frame 6. The coordinates of P3 on these two frames move with LEGO.
Fig. 7.
Fig. 7. Simulation reconstruction with moving LEGO and static background shown in reproduced RGB images. Simulation measurement and the reconstruction of 4 frames from a snapshot when the objects have (a) small motion, (b) fast motion, and (c) very fast motion with respect to the LED scanning. Notice that the GAP-TV algorithm is applied to generate the reconstruction results.
Fig. 8.
Fig. 8. Simulation reconstruction results with different reconstruction algorithms, i.e. GPSR4D (row 2, Visualization 8), GAP-TV4D (row 3, Visualization 9), DeSCI4D (row 4, Visualization 10). The results are based on 25% transmittance blue noise with $N_t = 8$ frames. Four reproduced RGB images are selected to be shown.
Fig. 9.
Fig. 9. Three imaging scenes: (a) color checker, (b) Super Mario and background, and (c) natural plants.
Fig. 10.
Fig. 10. Experimental reconstruction spectral video excerpts with $N_t = 4$ using GAPTV4D (top-right, Visualization 11) and DeSCI4D (bottom-right, Visualization 12) algorithm.
Fig. 11.
Fig. 11. Real data reconstruction spectral video excerpts with different temporal compression ratios, i.e. Color checker with $N_t = 4$ (Visualization 12) and $N_t = 6$ (Visualization 13), Super Mario with $N_t = 4$ (Visualization 14) and $N_t = 6$ (Visualization 15), Natural plants with $N_t = 4$ (Visualization 16) and $N_t = 6$ (Visualization 17).

Equations (15)

Equations on this page are rendered with MathJax. Learn more.

g(x,y)=f(x,y,λ,t)m(x,y,λ,t)dλdt,
fnx,ny,nt=Sfnx,ny,nt.
G=nl=1Nlnt=1NtCnl,ntFnl,nt+E,
Mnλ,nt=nl=1Nlsnl,nλCnl,nt.
G=nλ=1Nλnt=1NtMnλ,ntFnλ,nt+E.
H=[diag(M1,1),diag(M2,1),,diag(MNλ,1),diag(M1,2),,diag(MNλ,Nt)],
g=Hf+e.
f^=argmaxfexp{12σe2gHf22+logp(f)}=argminf12gHf22σe2logp(f).
f^=argminf12gHf22+λR(f).
fk+1=argminf12gHf22+ρ2f(zkuk)22,
zk+1=argminfλR(z)+ρ2z(fk+1+uk)22,
uk+1=uk+(fk+1zk+1),
fk+1=(zkuk)+H(gH(zkuk))(Diag(HH)+ρ),
zk+1=Dσ^k(fk+1+uk),
uk+1=uk+(fk+1zk+1),
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.