WO2020117245A1 - Spectromètre d'imagerie reconfigurable par ordinateur - Google Patents
Spectromètre d'imagerie reconfigurable par ordinateur Download PDFInfo
- Publication number
- WO2020117245A1 WO2020117245A1 PCT/US2018/064265 US2018064265W WO2020117245A1 WO 2020117245 A1 WO2020117245 A1 WO 2020117245A1 US 2018064265 W US2018064265 W US 2018064265W WO 2020117245 A1 WO2020117245 A1 WO 2020117245A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- mask
- encoded
- spectral components
- moving image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0229—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using masks, aperture plates, spatial light modulators or spatial filters, e.g. reflective filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/12—Generating the spectrum; Monochromators
- G01J3/14—Generating the spectrum; Monochromators using refracting elements, e.g. prisms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/12—Generating the spectrum; Monochromators
- G01J3/18—Generating the spectrum; Monochromators using diffraction elements, e.g. grating
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
Definitions
- Hyperspectral imaging spectrometers have been utilized for a wide range of applications, including, but not limited to environmental monitoring, biomedical imaging, surveillance, food safety inspection and control, agricultural monitoring, and mineralogy. Both spacebome and airborne hyperspectral imagers can provide critical data for a wide range of scientific, economic and national security needs. Despite the high expected value of Earth-scale spectral information, high cost and complexity has limited the number of fielded spaceborne hyperspectral imagers, making global, persistent coverage with rapid revisits infeasible.
- a computational imaging approach to hyperspectral imaging spectroscopy yields a spectrometer that can be miniaturized while maintaining high sensitivity.
- This type of spectrometer called a Computational Reconfigurable Imaging SPectrometer (CRISP)
- CRISP Computational Reconfigurable Imaging SPectrometer
- It uses a dual-disperser reimaging design, with a static coding mask in the image plane, and relies on platform motion or sensor scanning to achieve coding diversity.
- an inventive spectrometer can make spatial-spectral coded measurements of a three-dimensional (3D) data cube on a two-dimensional (2D) detector array, which has fewer detector elements than the mask has resolvable features.
- 3D three-dimensional
- 2D two-dimensional
- This measurement approach provides significant radiometry advantages over traditional pushbroom spectrometers. These radiometric advantages grow as the format size of the focal plane array increases and are especially pronounced for detector-noise limited measurement systems, such as uncooled microbolometer arrays.
- inventive imaging spectrometry extends these radiometric advantages further by providing an advantage called“coding gain”: As the along-track format size of the focal plane array grows beyond the minimum number of platform shifts needed to reconstruct the spatial/spectral data cube, the measurement becomes more overdetermined, and the SNR of the reconstructed spatial/spectral data cube continues to increase. The result is that the SNR can grow ever-larger in proportion the square root of the along-track focal plane array format size, for detector- and shot-noise-limited systems.
- an inventive imaging spectrometer makes many measurements in parallel, so it realizes a multiplex advantage (Fellgetf s advantage) that increases its SNR. Operating in the longwave infrared regime with noisy uncooled microbolometers further increases the impact of the multiplex advantage.
- the combined impact of the coding gain and the multiplex advantage results hyperspectral radiometric performance with an uncooled detector array that is competitive with conventional pushbroom imaging spectrometers with cooled detectors.
- inventive computational reconfigurable imaging spectrometers have other advantages as well. For example, they can make fully determined or overdetermined
- an inventive spectrometer can be used on a moving platform, such as a satellite or airplane. As a result, an inventive spectrometer doesn’t require moving parts, whereas an imaging spectrometer for a microscope may need the spectrograph hardware to move relative to a stationary microscope objective.
- any application which uses multi- or hyperspectral imaging can benefit from the potential for miniaturization offered by the inventive technology, particularly with the advent of small drones and satellites.
- a weather-sensing or terrain-mapping satellite, drone, balloon, airplane, or helicopter could use a computational reconfigurable imaging spectrometer.
- An inventive spectrometer could be used for remote sensing of land from space or airborne vehicles are also domains where hyperspectral imaging has been applied.
- Chemical hazard or trace gas detection is another application of hyperspectral sensing that could benefit from CRISP.
- Biomedical imaging applications include tissue spectroscopy, microscopy, and frequency- domain optical coherence tomography.
- Wide-area surveillance instruments which scan to cover a wide field of regard can also use CRISP, and be used to secure buildings, airports, or other territories.
- Urban imaging ground-based or airborne
- for monitoring of power usage, pollutants, and trace gases is another application for the inventive technology.
- An inventive spectrometer also called a hyperspectral imaging system, may include a first dispersive element, a coding mask in optical communication with the first dispersive element, a second dispersive element in optical communication with the coding mask, and a detector array, such as an uncooled microbolometer array configured to detect long-wave infrared radiation, in optical communication with the second dispersive element and having M3 N detector elements, where M and N are positive integers (e.g., MIN > 10 or 100).
- the first dispersive element disperses an image moving laterally with respect to an optical axis of the hyperspectral imaging system into N spectral components.
- the coding mask encodes the N spectral components with a predetermined code, such as a binary Walsh-Hadamard S-matrix code or a random binary code.
- the second dispersive element recombines the N spectral components into an encoded image.
- the detector array detects the encoded image.
- the coding mask can be reconfigured, e.g., to improve spatial and/or spectral resolution when reconstructing the laterally moving image from the laterally moving encoded image.
- Such a system may also include a processor, operably coupled to the detector array, to form a reconstructed version of the image from the encoded image with a coding gain ofjM /2N . If the image moves fast enough to cause motion blur in the reconstructed version of the image, the processor may estimate the motion blur from the reconstructed version of the image based on the predetermined code.
- the detector array can also acquire a sequence of encoded images as the image moves laterally with respect to the optical axis of the hyperspectral imaging system, in which case the processor can estimate a trajectory of the image with respect to the hyperspectral imaging system based on the sequence of encoded images.
- the hyperspectral imaging system may also include a first lens in optical communication with the coding mask and at least one second lens in optical communication with the coding mask and the detector array.
- the first lens focuses the N spectral components onto the coding mask.
- the second lens relays the encoded image to the detector array.
- An inventive spectrometer can also be used to carry out inventive methods of
- hyperspectral imaging An example of these inventive methods includes separating a moving image into N spectral components propagating orthogonal to a direction of motion of the moving image.
- a lens or other optical element focuses the N spectral components onto at least M
- a detector array such as an uncooled bolometer array, samples at least spatially resolvable spots of the encoded moving image at a frame rate / to form a hyperspectral data cube.
- a processor forms a reconstruction of the moving image with a coding gain of M/2/V based on the hyperspectral data cube.
- the binary static mask may apply a predetermined code to data in the hyperspectral data cube represented by a measurement matrix f, in which cases the processor forms the
- the processor may minimize a mean squared error (MSE) equal to (s 2 /M) Trace( p + p) -1 , where s 2 represents a variance of measurement noise associated with sampling the at least spatially resolvable spots of the encoded moving image.
- MSE mean squared error
- the processor can also estimate a trajectory of the moving image based on at least a portion of the hyperspectral data cube.
- FIG. 1 A shows a computational reconfigurable imaging spectrometer dispersing and encoding true images x with a static mask to form coded measurements.
- FIG. IB shows decoding of the coded measurements acquired with a computational reconfigurable imaging spectrometer to produce a final image.
- FIG. 1C illustrates how a computational reconfigurable imaging spectrometer measures the binary code of the mask shifted as a function of wavelength.
- FIG. ID illustrates how a single column of the binary mask in a computational reconfigurable imaging spectrometer contributes to the final image acquired by a computational reconfigurable imaging spectrometer.
- FIG. 2A shows different“fine” binary masks whose feature sizes remain constant with spectral bin width.
- FIG. 2B shows different“coarse” binary masks whose feature sizes scale with spectral bin width.
- FIG. 2C is a plot of noise equivalent temperature difference (NETD) versus spectral bin width for the fine binary masks of FIG. 2 A (upper trace) and the coarse binary masks of FIG. 2B (lower trace).
- FIG. 3 A shows a first composite image of the ground acquired by a computational reconfigurable imaging spectrometer aboard an airplane.
- FIG. 3B is a second composite image of the ground acquired by a computational reconfigurable imaging spectrometer aboard an airplane.
- FIG. 3C is a cross-correlation of the first composite image of FIG. 3A and the second composite image of FIG. 3B.
- FIG. 3D is a plot of pixel shift versus frame number for the derived from cross correlations like those in FIG. 3C for x motion (upper trace) and y motion (lower trace).
- FIG. 4A shows a schematic of a computational reconfigurable imaging spectrometer built on an optical breadboard.
- FIG. 4B shows a photograph of the computational reconfigurable imaging spectrometer shown in FIG. 4A.
- FIG. 5 A shows a randomly coded binary mask suitable for use in a computational reconfigurable imaging spectrometer.
- FIG. 5B shows a binary mask encoded with a Hadamard S-matrix suitable for use in a computational reconfigurable imaging spectrometer.
- FIG. 5C is a plot of simulated mean squared error (MSE) versus level of determination ( M/Nl ) of a computational reconfigurable imaging spectrometer measurement using the randomly patterned coding mask of FIG. 5 A versus the more efficient S-matrix (Hadamard) coding mask of FIG. 5B.
- MSE mean squared error
- M/Nl level of determination
- FIG. 6 A shows an image of a 100 °C blackbody (Cl-Sy stems SR-800) acquired with a computational reconfigurable imaging spectrometer (the arrow indicates direction of motion of the scene).
- FIG. 6B is a plot of raw blackbody spectra for blackbodies at 40 °C, 100 °C, and 138 °C (100-pixel by 100-pixel averages) acquired with a computational reconfigurable imaging spectrometer.
- FIG. 6C is shows the blackbody spectra of FIG. 6B adjusted for focal plane array wavelength sensitivity (the solid lines represent theoretically calculated blackbody curves (Planck’s law) for each temperature).
- FIG. 6D shows blackbody spectra as seen through 7.90 pm and 10.27 pm Spectrogon narrow bandpass filters with a computational reconfigurable imaging spectrometer (solid black lines represent filter transmission data provided by the vendor).
- FIG. 7 is a surface plot of the ratio of the signal-to-noise ratio (SNR) for a computational reconfigurable imaging spectrometer to the SNR of a pushbroom spectrometer with the same lenses and focal plane array versus number of image frames for blackbodies at 40 °C (upper trace) and 138 °C (lower trace).
- SNR signal-to-noise ratio
- FIG. 8A is a transmission spectrum of about 2000 ppm of R-134a gas at 21 °C (room temperature) acquired with a computational reconfigurable imaging spectrometer.
- FIG. 8B is an emission spectrum (left axis) and (1 - Transmission) (right axis) of about 2000 ppm of R- 134a gas at 50 °C acquired with a computational reconfigurable imaging spectrometer.
- FIG. 8C shows an image acquired with a computational reconfigurable imaging spectrometer of two polyethylene bags, one filled with air and one that contains -2000 ppm of R-134a gas mixed with air, placed side-by-side in front of a 100 °C blackbody
- FIG. 9 is a plot of the noise-equivalent temperature difference (NETD) of reconstructed spectra versus array/field-of-view size for predicted and simulated computational reconfigurable imaging spectrometers, for a spectrometer with an uncooled microbolometer detector array, and for NASA specifications.
- NETD noise-equivalent temperature difference
- FIGS. 1A and IB show a computational reconfigurable imaging spectrometer 100 that exploits platform motion, dispersive elements 110 and 150, and a reconfigurable coded focal - plane mask 130 to make a time series of encoded measurements of the optical spectrum at each pixel in a focal plane array 170.
- This encoding is inverted using specialized processing to recover the spectrum.
- this approach has sensitivity and other advantages over existing imaging spectrometer designs, enabling miniaturization and improved area coverage.
- the first dispersive element 110 is a prism, grating, or other optical element that spectrally disperses light from a scene x.
- the minimum resolvable spectral feature is limited by the product of the spatial point spread function extent in distance on the mask 130 and the change in wavelength per unit distance due to dispersion.
- the number of spectral unknowns should not be bigger than the number of imager rows, and the amount of spectral spread for the spectral bins should not be larger than the focal plane array 170.
- the spectral dispersion can be performed by the first dispersive element 110 in the along- track dimension or cross-track dimension (i.e., in the same direction as the scene motion or orthogonal to the direction of scene motion and the spectrometer’s optical axis z), with a suitable choice of coding mask 130.
- a first lens 120 creates a dispersed image of the scene on the mask 130, which encodes the dispersed image with a known pattern.
- the second dispersive element 140 e.g., another prism or grating
- a telescope formed of lenses 140 and 160 images the encoded, dispersed image F onto the focal plane array 170, which has fewer pixels (detector elements) than the mask.
- the focal plane array 170 samples the resulting 2D radiance field, which contains light summed from the whole spectrum, but with a position-dependent weighting as a function of wavelength to create coded measurements y.
- Platform motion or scanning in a direction orthogonal to the spectrometer’s optical axis z causes the scene to move relative to the mask 130 and the focal plane array 170.
- the frame rate is equal to the time of flight over the instantaneous field of view of a single row of pixels in the focal plane array 170, but other scenarios are possible. This produces a time modulated intensity pattern which is sampled at the frame rate of the focal plane array 170.
- the spectrometer’s spatial and spectral characteristics are reconfigurable, simply by swapping or changing the mask 130, for example, with a wheel or translation stage.
- the mask 130 itself can be static, so the spectrometer 100 does not require any active components.
- the spectrometer 100 can operate over a broad range of wavelengths, limited primarily by the responsivity of the focal plane array 170. If the focal plane array 170 is an uncooled (micro)bolometer array, for instance, the spectrometer 100 may operate in the long-wave infrared (LWIR) region of the electromagnetic spectrum (i.e., a wavelength range of about 8-14 microns). And if the focal plane array 170 is a CCD or CMOS array, the spectrometer 100 may operate in the visible or near-infrared (NIR) regions of the visible spectrum.
- LWIR long-wave infrared
- NIR near-infrared
- FIGS. 1C and ID illustrate how the spectrometer 100 of FIGS. 1 A and IB modulates the spectrum with a binary code from the mask 130.
- FIG. 1C shows that code is shifted as a function of wavelength.
- coded images 131, 132, and 133 there are three coded images 131, 132, and 133, respectively, each shifted vertically with respect to the other coded images.
- a superposition 134 of these coded images 131-134 is recombined by the second dispersive element 150 and sensed by the focal plane array 170.
- FIG. ID shows that in an along-track measurement, wavelength- shifted modulation combined with along-track platform motion leads to a system of equations constructed from the mask’ features.
- one column of the binary mask yields vertically shifted measurements as in FIG. 1C.
- Along-track motion converts this code into a“known measurement matrix cp,” which, when multiplied by the spatial/spectral data cube x yields the CRISP coded measurement y.
- the measurement matrix f is well-determined or overdetermined.
- the Walsh-Hadamard S-matrix code is an optimal encoding in principle but in practice can lack robustness if, for example, the number of measurements is not a power of 2 (minus one) or if there is deviation from perfect binary modulation, grid alignment, and rectilinear motion.
- a processor 180 coupled to focal plane array 170 uses the pseudoinverse of the measurement matrix, cp + , to compute a reconstructed data cube x, also shown in FIG. IB.
- the reconstructed data cube, or spectra, x minimize the mean squared errors between the
- the forward matrix model is derived from the code used in the binary mask.
- the processor 180 can perform the reconstruction efficiently, in part by exploiting repetitive and cyclic patterns within the forward matrix model f and/or S-matrix orthogonality to reduce computation time. More specifically, the S-matrix is derived from an orthonormal matrix (with scale and offset), so its inverse is straightforward to compute. In general, the spectra in a given column of the spectral data cube are modulated by the same column of the mask, so each column of the spectral data cube has its own x N system matrix f. Thus, for each column of the spatial data cube, the processor computes a single pseudoinverse f + and uses it to solve every spectrum in that column.
- the spectrometer 100 measures an image of the mask under monochromatic light before (or during or after) the scene measurement. This measured image of the mask is then used to determine how the spectrometer 100 actually encodes the image, and this assumption is used in processing. Use of the measured image of the mask in the processing, rather than the assumed“as designed” binary mask, is more robust because it accounts for the point spread, aberrations, distortions, and actual mask (mis)alignment relative to the pixel grid.
- the real“as-measured” mask may not be perfectly binary but may include some gray-scale features due to the real effects.
- the number of features or pixels in the mask can be larger or smaller than the number of pixels in the focal plane array, with corresponding tradeoffs.
- Masks with more features (“pixels”) than detector elements in the focal plane array can enable increased spatial and spectral resolution.
- Masks with fewer features reduce resolution but enable higher signal-to-noise ratio measurements over a more coarsely sampled grid.
- the number of frames F times the number of pixels P should be greater than or equal to NP.
- the compressive case can be accomplished by minimizing the residuals on the measurements while minimizing a data regularization term such as the Frobenius norm of the recovered spectral matrix x (a low rank prior).
- a minimization can be carried out in many ways, including via the alternating direction method of multipliers (ADMM) method.
- a good nominal operating point (though not the only operating point) is to match the number of mask features in the mask 130 to the number of detectors in the focal plane array 170 plus the number of wavelength bins minus 1. For example, if the focal plane array 170 in FIGS.
- the mask can be reconfigured between measurements. That is, the mask can be replaced, translated laterally, rotated, or changed between measurements. Shifting or rotating the mask or replacing the mask with a different mask causes the computational reconfigurable imaging spectrometer to modulate spectral data cube with a different code. A different code may be advantageous for a number of reasons. By swapping the mask with a different one, the spectrometer can be reconfigured for optimal performance at different spectral or spatial resolutions.
- the imager (focal plane array) has M rows of detector elements (pixels), the spectrum is divided into N spectral bins, and the mask contains M + N- 1 rows of features. If the system doesn’t need to operate with its full spectral resolution, the spectrum can be divided into coarser bins (e.g., into N/ 2, N/ 4, or L75 spectral bins), with each bin treated as an unknown to be reconstructed. In that case, models suggest that a mask which has coarser spatial features (e.g., 2, 4, or 5 times the width of the detector element) can encode the coarser spectrum such that it can be recovered with lower mean square error.
- coarser spatial features e.g., 2, 4, or 5 times the width of the detector element
- FIGS. 2A and 2B show different mask options.
- the mask has fine features whose sizes remain constant while the wavelength bins grow wider.
- the mask has coarse features whose sizes scale with the wavelength bin width.
- spatial masks with finer spatial features may be better for reconstructing data at subpixel spatial or fmer-than-a-bin spectral resolutions, while spatial masks with coarser spectral features may be better optimal for reconstructing data at coarser-than-a-bin spectral resolutions.
- the mask’s predetermined code may be selected to improve
- the mask’s predetermined code may enable improved reconstructions at high spectral resolution (e.g., finer spectral bins than the usual N spectral bins) or low spectral resolution (e.g., coarser spectral bins than the usual N spectral bins).
- a mask may also have spatial features of different sizes (e.g., smaller than an imager pixel or larger than an imager pixel) in the direction parallel to dispersion.
- FIG. 2C is a plot of the noise equivalent temperature difference (NETD) versus spectral bin width for binary masks with fine features (upper trace) and coarse features (lower trace) with the same optics.
- the upper axis shows spectral bin width normalized to pixel size and the It shows that, for the same optics, varying the mask feature size improves performance for recovering a coarser spectrum.
- the upper trace was derived by using a single static mask (e.g., as in FIG. 2A), solving for the spectrum with N bins, and spectrally averaging and downsampling to get the spectrum at successively coarser spectral bin sizes.
- the lower trace was derived by using a reconfigurable mask with coarser spatial features (e.g., as in FIG. 2B) and solving for the successively coarser spectral bins directly.
- the performance, shown as NETD, is better with the reconfigurable mask.
- the mask features to be sub-pixel in size enables digital super resolution for systems that are not diffraction limited.
- a mask with features smaller than a detector pixel will encode the spectral cube with fine temporal resolution, and thus enable reconstruction of the spatial data cube at subpixel spatial resolutions.
- the spatio-spectral modulation encodes information that enables recovery of spectrum at a subpixel resolution.
- the point spread function can be larger than a focal plane array pixel, even if mask feature is same size as focal plane array pixel. In some circumstances, the PSF may limit spatial and spectral resolution.
- a computational reconfigurable imaging spectrometer can be more sensitive than other imaging spectrometers, leading to improvements in area coverage rates and allowing for lower- size/weight/power instruments to achieve the same SNR as other imaging spectrometers.
- the CRISP approach for hyperspectral imaging enables two sensitivity gains, the Fellgett or multiplex advantage found in other multiplexed measurement (e.g., Hadamard spectrometer) designs and a“coding gain.” These gains can be obtained simultaneously in practical designs and applications.
- the sensitivity enhancement can be engineered by increasing the number of along-track measurements (hence imager rows), for background- and dark-noise limited spectrometers. To see how, consider a spectrometer with N l wavelength bins that makes M along-track
- CRISP gains a multiplex sensitivity advantage (Fellgett’ s advantage) in SNR in the case where signal-independent detector noise is dominant.
- the SNR gain is L /L ⁇ /2 relative to a dispersive slit spectrometer that measures the signal in each spectral bin separately.
- SNR is defined as average signal intensity over the signal reconstruction error standard deviation.
- MSE The expected mean squared error for such a system
- the number of wavelengths of interest N l can be exceeded by the number of measurements M in part due to the format and frame rates achieved by modern focal plane arrays.
- M > N l an additional coding gain of M/N l is achievable.
- M/N l is an integer and the measurement matrix consists of multiple vertically concatenated S-matrices, i.e., [0067] As shown below, this gain can also be obtained for M sufficiently larger than N l using random binary matrices.
- the multiplex and coding gains can be obtained simultaneously, resulting in an SNR enhancement for CRISP versus dispersive slit designs of M /2.
- This combined sensitivity gain enables sensitive measurements using lower performance detectors such as uncooled microbolometers.
- SNR enhancement relative to existing FTIR and dispersive slit spectrometers is possible as well via the coding gain.
- FTIR Fourier- transform infrared
- dispersive spectrometers are possible as well. Specifically, if M > N l then the reconstruction problem is overdetermined, and SNR is enhanced by a factor called the coding gain, J M/2N l .
- FTIR and dispersive spectrometer types can produce coding gain only with prohibitive difficulty, such as replicating large portions of the instrument to get an overdetermined system like CRISP.
- FTIR systems possess a Fellgett advantage of L [Nc/ (4v/2) over dispersive systems when dark-noise limited.
- this advantage is less than that of CRISP by a factor of 2.8 ⁇ M/N due to losses and inefficiencies resulting from the beam splitter, cosine modulation, and symmetric interferograms.
- imaging FTIR spectrometers typically require complex moving parts. They must also sweep through collection of the entire interferogram before the scene shifts out of view, which becomes impractical on a moving platform without mounting the spectrometer or imaging optics on a gimbal.
- variable-gap etalon spectrometer has about twice the loss of an FTIR spectrometer, giving it an SNR that is at least about 5.6 times smaller than the SNR of a computational reconfigurable imaging spectrometer. It also has limitations on spectral resolution and spectral coverage versus CRISP.
- a computational reconfigurable imaging spectrometer has several advantages over existing technologies, including conventional dispersive spectrometers and FTIR spectrometers.
- a computational reconfigurable imaging spectrometer can achieve both a multiplex (Fellgett) gain and a coding gain, improving the system SNR beyond what is possible via alternative architectures assuming a sufficiently large focal plane array.
- the coding gain enables SNR improvements even under shot-noise limited imaging conditions, under which traditional multiplexed spectral imaging architectures can at best offer no loss relative to non-multiplexed designs.
- a computational reconfigurable imaging spectrometer therefore provides a sensitivity gain in both shot-noise limited and detector-noise limited cases compared to existing
- a computational reconfigurable imaging spectrometer does not require moving parts. These advantages translate directly into reduced size, weight, and power (SWaP) requirements, as expensive cooling components can be either be reduced or eliminated, and compact optical designs may be used.
- SWaP reduced size, weight, and power
- a computational reconfigurable imaging spectrometer can be reconfigured by swapping or changing the coding mask mid-mission.
- a computational reconfigurable imaging spectrometer sensitivity gains can yield improvements in area coverage rates and allow for smaller, lighter, more efficient instruments with the same SNR as existing systems.
- the computational reconfigurable imaging spectrometer 100 produces a raw product (the coded measurement y) that includes (a) scene imagery x and (b) fixed pattern imprinting from the mask features.
- This raw image product contains direct 2D scene imagery enables the use of computer-vision motion tracking to remove jitter and nonlinear motion. This is in direct contrast to other hyperspectral imaging modalities.
- Pushbroom system do not produce 2D images, making any form of computer- vision motion tracking very difficult or impossible.
- Michelson systems produce interferograms again making motion tracking very difficult.
- each CRISP frame constitutes a complete x-y image (in contrast to pushbroom spectrometry, where each frame is a single spatial column with the orthogonal axis being spectral wavelength)
- the cross-correlation between successive frames can be used to provide registration and remove motion artifacts due to camera motion caused by movement of its platform.
- Such platform motion can be especially detrimental to any technique that creates a final image product from consecutive frames (e.g. FTIR-based spectrometers or CRISP spectrometers).
- FIGS. 3 A-3D illustrates actual motion compensation used in a flight test of the CRISP system.
- FIGS. 3A and 3B show complete x-y images of the ground taken with an airborne CRISP system. The image in FIG. 3B was acquired 50 frames after the image in FIG. 3 A.
- FIG. 3C shows the cross-correlation between the images in FIGS. 3A and 3B.
- FIG. 3D is a plot of the location of the cross-correlation peak for cross-correlation functions repeated for many successive frames. It represents the relative motion track of the scene. This track of the scene is used to register the acquired frames to a common reference, thereby removing motion artifacts.
- a digital focal plane could be used to enhance deblurring.
- the method for doing this is to change the gain in the detector during the blurring and to vary the gains chosen across blur measurements. This ability to modulate the blur functions enriches the types of encoding that can be achieved and can yield superior reconstruction quality.
- Motion deblurring works by design - the modulation of the blur function encodes the high frequency information that would otherwise be lost in the low pass filter of the blur function. Essentially, this is a spreading code that is applied to the signal prior to the blur, thus ensuring that high frequency information gets mapped into the low frequency region admitted by the motion blur.
- the relationship between the unknown 3D spatial/spectral data cube x and the measurements y can be modeled via a linear function.
- the size of the vector y is the number of pixels times the number of frames.
- the unknown 3D spatial/spectral data cube has a size equal to the number of detector elements by number of wavelength bins.
- FIGS. 4A and 4B show a schematic and a photograph, respectively, of a computational reconfigurable imaging spectrometer 400 built using commercial off-the-shelf (COTS) components.
- the spectrometers 400 was built for concept demonstration and thus has not been optimized for optical throughput. It operates in the longwave infrared (LWIR) regime using an uncooled microbolometer camera 470 (significantly noisier than state-of-the-art HgCdTe detectors) because a detector-noise limited system best shows its advantages.
- LWIR longwave infrared
- a first ZnSe dispersive prism 410 (without antireflective coating) separates light into different wavelength bins.
- a first 52 mm, f/1 lens 420 (aspheric triplets; FLIR model A655sc) focuses the dispersed light onto an antireflective coated ZnSe coding mask 430 with a random checkerboard pattern applied using contact lithography.
- the feature sizes of the coding mask 430 are matched to the pixel size of the imager 470.
- a second 52 mm, f/1 lens 440 and a second ZnSe dispersive prism 450 form a spectrally combined, encoded image that is sampled by the microbolometer camera 470 at 50 frames per second (fps) with an array of 640 x 480 pixels, each of with is 17 pm square. At 17 pm, the pixels of the imager 470 and the coding mask 430 are smaller than the diffraction-limited point spread function (PSF) in the LWIR, which is about 20-34 pm.
- PSF diffraction-limited point spread function
- FIG. 5A shows a random checkerboard pattern like the one applied to the mask 430.
- FIG. 5B shows a better-designed mask - in this case, a binary Hadamard S-matrix.
- the mask“truth” i.e., a picture of the mask used to decode the hyperspectral data cube
- MSE expected mean squared error
- FIGS. 6A-6D show spectral measurements of an extended source blackbody (Cl-Systems SR-800) acquired with the computation reconfigurable imaging spectrometer 400 of FIGS. 4 A and 4B.
- FIG. 6 A shows a single frame of the calibrated blackbody source at 100 °C imaged through the spectrometer 400.
- the blackbody source was positioned about 3 m from the spectrometer 400 and translated via a motorized translation stage to simulate platform motion (the white arrow in FIG. 6A indicates the direction of motion).
- the spectrometer 400 captured a total of 335 frames of the blackbody source moving across the image plane, making the dataset overdetermined by a factor of five.
- FIG. 6B shows the resultant raw spectra (after decoding the dataset) of the calibrated blackbody source at 40 °C, 100 °C, and 138 °C. After adjusting for the wavelength specific sensitivity of the microbolometer 470, the resulting spectral shape matches perfectly to blackbody radiation theory (solid black lines), as shown in FIG. 6C.
- FIG. 6C shows spectra of a pair of narrow bandpass filters (Spectrogon Inc.) between the spectrometer 400 and the blackbody source.
- the filters had center wavelengths of 7.90 pm and 10.27 mih.
- the resulting spectral shapes agree well with the filter transmission data (black lines) provided by the vendor.
- the spectra shown in FIGS. 6B-6D were obtained by subtracting a background measurement (in this case, a measurement with the spectrometer input blocked) from the raw spectra. This background subtraction helps to reduce or eliminate trends in the baseline or any clutter in the spectrum due to narcissus or other effects.
- a background measurement in this case, a measurement with the spectrometer input blocked
- the coding gains enable enhanced sensitivity for both shot-noise limited and detector-noise limited systems. These gains become apparent when comparing the performance of a computational reconfigurable imaging spectrometer to the performance of a pushbroom spectrometer.
- the pushbroom spectrometer is the spectrometer 400 in FIGS. 4A and 4B reconfigured for pushbroom operation, i.e., with the first prism 410 replaced with an uncoated ZnSe window of similar size (to keep the optical throughput consistent) and the coding mask 430 replaced with a slit whose width is matched to the pixel size of the imager 470. Both configurations were used to image the same 40 °C and 138 °C blackbody sources.
- FIG. 7 shows the ratio of SNRs of the spectrometer 400 (SNRCRISP) compared to same spectrometer reconfigured to operate as a pushbroom spectrometer (SNRpushbroom) versus number of measurements (frames) M for the CRISP measurements (the pushbroom measurements are single-shot measurements).
- the CRISP SNR gains grow as M (shown by black lines). After 335 frames, CRISP outperforms the pushbroom by a factor of five for the 40 °C blackbody source and by a factor of three for the 138 °C blackbody source.
- the measured performance enhancements at 335 frames for both the 40 °C and 138 °C blackbody sources are within the theoretical upper-bound expected for a detector-noise limited measurement ( M / 2 ⁇ 9), and the theoretical lower-bound expected for a shot-noise limited measurement ⁇
- the spectrometer 400 has an optical throughput of 11%, which suggests 40-70 mK effective shot noise for the measurements of the 40 °C and 138 °C blackbody sources.
- the measured detector noise is 100-200 mK, which puts the measurements in a modestly detector- noise limited regime.
- the SNR should improve as M, which is clearly shown in FIG. 7. Larger format imagers should increase the CRISP SNR. Detector-noise limited performance can be approached via a combination of higher optical throughput and a more carefully measured mask function used to decode the hyperspectral cube.
- FIGS. 8A-8D illustrate the spectrometer’s ability to perform spectroscopic analysis of trace gases.
- the spectrometer 400 made hyperspectral measurements of a thin polyethylene bag (about 80% transmissive in the LWIR) with a trace amount (about 2000 ppm) of 1,1,1,2-tetrafluoroethane (R-134a) gas, buffered with nitrogen to 1 atm.
- the bag was positioned between the spectrometer and a blackbody source set to 100 °C.
- FIG. 8A shows measured (squares) and computed (line) transmission spectra of the bag at room temperature.
- the measured transmission spectrum is referenced to a measurement of a polyethylene bag filled with only nitrogen and closely resembles the computed transmission spectrum, which was computed from data in the Pacific Northwest National Laboratory (PNNL) database available at https://nwir.pnl.gov/.
- PNNL Pacific Northwest National Laboratory
- FIG. 8B shows measured and computed emission spectra for the same bag of R- 134a heated to a temperature of about 50 °C.
- the measured emission spectrum in FIG. 8B is the difference between the raw measurement and a room-temperature blackbody spectrum baseline. Again, the measured data closely matches the emission spectrum computed from the PNNL database.
- the spectra displayed in FIGS. 8A and 8B are the result of 100-pixel-by- 100-pixel averaging.
- FIGS. 8C and 8D show the computational reconfigurable imaging spectrometer’s ability to detect and locate a trace gas within the scene.
- the bag of trace R-134a gas was placed next to a bag full of room air between a 100 °C blackbody radiator, shown at left in FIG. 8C, and the spectrometer 400.
- the bag of air in this case simply exists as the null experiment.
- FIG. 8D is a heat map produced from the data of FIG. 8C that shows the location of the R-134a in the scene.
- the heat map is produced by simply plotting the height of a known peak in the transmission spectrum of the R-134a gas.
- a white dotted border depicts the perimeter of the extended blackbody source.
- FIG. 9 illustrates the performance potential of a computational reconfigurable imaging spectrometer for space-based infrared sounding. It shows the predicted and simulated noise equivalent temperature difference (NETD) for a computational reconfigurable imaging spectrometer and for a standard pushbroom spectrometer versus detector array size for the NASA Atmospheric Infrared Sounder (AIRS) specification.
- the NASA AIRS specification is for an LWIR temperature sounding sensor with roughly 12 nm spectral resolution in a spectral band from 13.5 pm to 15 pm, and 8 km spatial resolution.
- FIG. 9 shows that a computational reconfigurable imaging spectrometer with an uncooled microbolometer array as the focal plane array over a pushbroom spectrometer with an uncooled microbolometer.
- This analysis assumes that the both spectrometers operate in the presence of the same spatial and temporal noise.
- the NETD of a CRISP system is more than an order of magnitude lower for modest focal plane sizes and can match NASA sounder sensitivity requirements.
- Analyses of other hyperspectral and multispectral imaging systems show similar gains, covering a wide gamut of spatial and spectral ranges and resolution capabilities.
- the computational reconfigurable imaging spectrometer would also have significant SWaP advantages over a spectrometer with a cooled HgCdTe focal plane array.
- inventive embodiments are presented by way of example only and that, within the scope of the appended claims and equivalents thereto, inventive embodiments may be practiced otherwise than as specifically described and claimed.
- inventive embodiments of the present disclosure are directed to each individual feature, system, article, material, kit, and/or method described herein.
- a computer may be embodied in any of a number of forms, such as a rack-mounted computer, a desktop computer, a laptop computer, or a tablet computer. Additionally, a computer may be embedded in a device not generally regarded as a computer but with suitable processing capabilities, including a Personal Digital Assistant (PDA), a smart phone or any other suitable portable or fixed electronic device.
- PDA Personal Digital Assistant
- a computer may have one or more input and output devices. These devices can be used, among other things, to present a user interface. Examples of output devices that can be used to provide a user interface include printers or display screens for visual presentation of output and speakers or other sound generating devices for audible presentation of output.
- Examples of input devices that can be used for a user interface include keyboards, and pointing devices, such as mice, touch pads, and digitizing tablets.
- a computer may receive input information through speech recognition or in other audible format.
- Such computers may be interconnected by one or more networks in any suitable form, including a local area network or a wide area network, such as an enterprise network, and intelligent network (IN) or the Internet.
- networks may be based on any suitable technology and may operate according to any suitable protocol and may include wireless networks, wired networks or fiber optic networks.
- the various methods or processes may be coded as software that is executable on one or more processors that employ any one of a variety of operating systems or platforms. Additionally, such software may be written using any of a number of suitable programming languages and/or programming or scripting tools, and also may be compiled as executable machine language code or intermediate code that is executed on a framework or virtual machine.
- inventive concepts may be embodied as a computer readable storage medium (or multiple computer readable storage media) (e.g., a computer memory, one or more floppy discs, compact discs, optical discs, magnetic tapes, flash memories, circuit configurations in Field Programmable Gate Arrays or other semiconductor devices, or other non-transitory medium or tangible computer storage medium) encoded with one or more programs that, when executed on one or more computers or other processors, perform methods that implement the various embodiments of the invention discussed above.
- the computer readable medium or media can be transportable, such that the program or programs stored thereon can be loaded onto one or more different computers or other processors to implement various aspects of the present invention as discussed above.
- program or“software” are used herein in a generic sense to refer to any type of computer code or set of computer-executable instructions that can be employed to program a computer or other processor to implement various aspects of embodiments as discussed above. Additionally, it should be appreciated that according to one aspect, one or more computer programs that when executed perform methods of the present invention need not reside on a single computer or processor but may be distributed in a modular fashion amongst a number of different computers or processors to implement various aspects of the present invention.
- Computer-executable instructions may be in many forms, such as program modules, executed by one or more computers or other devices.
- program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types.
- functionality of the program modules may be combined or distributed as desired in various embodiments.
- data structures may be stored in computer-readable media in any suitable form.
- data structures may be shown to have fields that are related through location in the data structure. Such relationships may likewise be achieved by assigning storage for the fields with locations in a computer-readable medium that convey relationship between the fields.
- any suitable mechanism may be used to establish a relationship between information in fields of a data structure, including through the use of pointers, tags or other mechanisms that establish relationship between data elements.
- inventive concepts may be embodied as one or more methods, of which an example has been provided.
- the acts performed as part of the method may be ordered in any suitable way. Accordingly, embodiments may be constructed in which acts are performed in an order different than illustrated, which may include performing some acts simultaneously, even though shown as sequential acts in illustrative embodiments.
- a reference to“A and/or B”, when used in conjunction with open-ended language such as“comprising” can refer, in one embodiment, to A only (optionally including elements other than B); in another embodiment, to B only (optionally including elements other than A); in yet another embodiment, to both A and B (optionally including other elements); etc.
- “or” should be understood to have the same meaning as“and/or” as defined above.
- “or” or“and/or” shall be interpreted as being inclusive, i.e., the inclusion of at least one, but also including more than one, of a number or list of elements, and, optionally, additional unlisted items. Only terms clearly indicated to the contrary, such as“only one of’ or“exactly one of,” or, when used in the claims,“consisting of,” will refer to the inclusion of exactly one element of a number or list of elements.
- the phrase“at least one,” in reference to a list of one or more elements, should be understood to mean at least one element selected from any one or more of the elements in the list of elements, but not necessarily including at least one of each and every element specifically listed within the list of elements and not excluding any combinations of elements in the list of elements.
- This definition also allows that elements may optionally be present other than the elements specifically identified within the list of elements to which the phrase“at least one” refers, whether related or unrelated to those elements specifically identified.
- “at least one of A and B” can refer, in one embodiment, to at least one, optionally including more than one, A, with no B present (and optionally including elements other than B); in another embodiment, to at least one, optionally including more than one, B, with no A present (and optionally including elements other than A); in yet another embodiment, to at least one, optionally including more than one, A, and at least one, optionally including more than one, B (and optionally including other elements); etc.
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
La présente invention concerne un système d'imagerie hyperspectrale comprenant : un premier élément de dispersion pour disperser une image se déplaçant latéralement par rapport à un axe optique du système d'imagerie hyperspectrale en composants spectraux N, un masque de codage, en communication optique avec le premier élément de dispersion, pour coder les composants spectraux N avec un code prédéterminé, un second élément de dispersion, en communication optique avec le masque de codage, pour recombiner les composants spectraux N en une image codée, une rangée de détecteurs, en communication optique avec le second élément de dispersion et ayant des éléments détecteurs M >/= N, pour détecter l'image codée, M et N étant des nombres entiers positifs et un processeur, fonctionnellement accouplé à la rangée de détecteurs, pour former une version reconstruite de l'image à partir de l'image codée, le processeur configuré pour estimer le flou de mouvement à partir de la version reconstruite de l'image sur la base du code prédéterminé.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/US2018/064265 WO2020117245A1 (fr) | 2018-12-06 | 2018-12-06 | Spectromètre d'imagerie reconfigurable par ordinateur |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/US2018/064265 WO2020117245A1 (fr) | 2018-12-06 | 2018-12-06 | Spectromètre d'imagerie reconfigurable par ordinateur |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2020117245A1 true WO2020117245A1 (fr) | 2020-06-11 |
Family
ID=70975408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/US2018/064265 Ceased WO2020117245A1 (fr) | 2018-12-06 | 2018-12-06 | Spectromètre d'imagerie reconfigurable par ordinateur |
Country Status (1)
| Country | Link |
|---|---|
| WO (1) | WO2020117245A1 (fr) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113418873A (zh) * | 2021-06-23 | 2021-09-21 | 中国科学院长春光学精密机械与物理研究所 | 一种高光谱成像系统及重构光谱成像方法 |
| CN113791037A (zh) * | 2021-08-19 | 2021-12-14 | 南京航空航天大学 | 一种基于生成对抗网络的硅基傅里叶变换光谱测量方法 |
| CN115330649A (zh) * | 2022-07-29 | 2022-11-11 | 融智生物科技(青岛)有限公司 | 质谱图的合成方法、系统、装置及介质 |
| CN115880199A (zh) * | 2023-03-03 | 2023-03-31 | 湖南大学 | 一种长波红外高光谱与多光谱图像融合方法、系统及介质 |
| CN118411441A (zh) * | 2024-05-08 | 2024-07-30 | 西湖智能视觉科技(杭州)有限公司 | 基于特种光纤单曝光光谱成像的内窥检测系统及其方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020057431A1 (en) * | 1999-04-09 | 2002-05-16 | Fateley William G. | System and method for encoded spatio-spectral information processing |
| US20070097363A1 (en) * | 2005-10-17 | 2007-05-03 | Brady David J | Coding and modulation for hyperspectral imaging |
| US20160165155A1 (en) * | 2013-07-18 | 2016-06-09 | Southwest Jiaotong University | Method for reducing image fuzzy degree of tdi-ccd camera |
| WO2016191157A1 (fr) * | 2015-05-22 | 2016-12-01 | Microsoft Technology Licensing, Llc | Système d'imagerie hyperspectrale de mouvement à double capteur |
-
2018
- 2018-12-06 WO PCT/US2018/064265 patent/WO2020117245A1/fr not_active Ceased
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020057431A1 (en) * | 1999-04-09 | 2002-05-16 | Fateley William G. | System and method for encoded spatio-spectral information processing |
| US20070097363A1 (en) * | 2005-10-17 | 2007-05-03 | Brady David J | Coding and modulation for hyperspectral imaging |
| US20160165155A1 (en) * | 2013-07-18 | 2016-06-09 | Southwest Jiaotong University | Method for reducing image fuzzy degree of tdi-ccd camera |
| WO2016191157A1 (fr) * | 2015-05-22 | 2016-12-01 | Microsoft Technology Licensing, Llc | Système d'imagerie hyperspectrale de mouvement à double capteur |
Non-Patent Citations (2)
| Title |
|---|
| MAHGOUB, A. ET AL.: "Aligning the frames of a non stationary imaging Fourier transform spectrometer for spectrum retrieva", 2009 16TH IEEE INTERNATIONAL CONFERENCE ON IMAGE PROCESSING, November 2009 (2009-11-01), pages 573, XP031628545 * |
| SULLENBERGER, R. ET AL.: "Computational reconfigurable imaging spectrometer", OPT. EXPRESS, vol. 25, no. 25, 11 December 2017 (2017-12-11), pages 31960 - 31969 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113418873A (zh) * | 2021-06-23 | 2021-09-21 | 中国科学院长春光学精密机械与物理研究所 | 一种高光谱成像系统及重构光谱成像方法 |
| CN113418873B (zh) * | 2021-06-23 | 2022-05-17 | 中国科学院长春光学精密机械与物理研究所 | 一种高光谱成像系统及重构光谱成像方法 |
| CN113791037A (zh) * | 2021-08-19 | 2021-12-14 | 南京航空航天大学 | 一种基于生成对抗网络的硅基傅里叶变换光谱测量方法 |
| CN115330649A (zh) * | 2022-07-29 | 2022-11-11 | 融智生物科技(青岛)有限公司 | 质谱图的合成方法、系统、装置及介质 |
| CN115880199A (zh) * | 2023-03-03 | 2023-03-31 | 湖南大学 | 一种长波红外高光谱与多光谱图像融合方法、系统及介质 |
| CN118411441A (zh) * | 2024-05-08 | 2024-07-30 | 西湖智能视觉科技(杭州)有限公司 | 基于特种光纤单曝光光谱成像的内窥检测系统及其方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10909670B2 (en) | Computational reconfigurable imaging spectrometer | |
| US11778289B2 (en) | Multi-camera imaging systems | |
| WO2020117245A1 (fr) | Spectromètre d'imagerie reconfigurable par ordinateur | |
| Scott et al. | The NN-explore exoplanet stellar speckle imager: Instrument description and preliminary results | |
| Yang et al. | Compressive hyperspectral imaging via sparse tensor and nonlinear compressed sensing | |
| Hu et al. | Practical snapshot hyperspectral imaging with DOE | |
| Yang et al. | Deep‐learning based on‐chip rapid spectral imaging with high spatial resolution | |
| Mu et al. | Snapshot hyperspectral imaging polarimetry with full spectropolarimetric resolution | |
| Sullenberger et al. | Computational reconfigurable imaging spectrometer | |
| Oktem et al. | High-resolution multi-spectral imaging with diffractive lenses and learned reconstruction | |
| US11895405B2 (en) | Method and system for performing high speed optical image detection | |
| Hagen et al. | Real-time quantitative hydrocarbon gas imaging with the gas cloud imager (GCI) | |
| Coluccia et al. | Optical compressive imaging technologies for space big data | |
| Barducci et al. | Developing a new hyperspectral imaging interferometer for earth observation | |
| Bennett et al. | Hyperspectral imaging in the infrared using LIFTIRS | |
| Murguia et al. | Compact visible/near-infrared hyperspectral imager | |
| Barducci et al. | ALISEO: a new stationary imaging interferometer | |
| Kumar et al. | Advances in spaceborne hyperspectral imaging systems. | |
| Vujkovic-Cvijin et al. | Adaptive spectral imager for space-based sensing | |
| Bézy et al. | The Copernicus Sentinel-5 mission for operational atmospheric monitoring: status and developments | |
| Barducci et al. | Simulation of the performance of a stationary imaging interferometer for high-resolution monitoring of the Earth | |
| Hinojosa et al. | Compressive spectral imaging using multiple snapshot colored-mosaic detector measurements | |
| Pagnutti et al. | Methods, procedures, and example results for evaluating commercial off-the-shelf hyperspectral instruments for use as satellite cross-calibration radiometers | |
| Szajewski et al. | Review on infrared single-pixel imaging | |
| Silny et al. | Large format imaging spectrometers for future hyperspectral Landsat mission |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 18942382 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 18942382 Country of ref document: EP Kind code of ref document: A1 |