CN108663118A - Infrared broadband EO-1 hyperion calculates imaging devices and methods therefor - Google Patents
Infrared broadband EO-1 hyperion calculates imaging devices and methods therefor Download PDFInfo
- Publication number
- CN108663118A CN108663118A CN201810050205.9A CN201810050205A CN108663118A CN 108663118 A CN108663118 A CN 108663118A CN 201810050205 A CN201810050205 A CN 201810050205A CN 108663118 A CN108663118 A CN 108663118A
- Authority
- CN
- China
- Prior art keywords
- infrared
- light modulator
- haplotype
- spatial light
- matrix
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims description 15
- 230000003595 spectral effect Effects 0.000 claims abstract description 47
- 102000054766 genetic haplotypes Human genes 0.000 claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 32
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 36
- 230000005855 radiation Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000003287 optical effect Effects 0.000 claims description 13
- 238000004611 spectroscopical analysis Methods 0.000 claims description 11
- 239000013598 vector Substances 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 9
- 238000000701 chemical imaging Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000005538 encapsulation Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 4
- 241000287196 Asthenes Species 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- SBIBMFFZSBJNJF-UHFFFAOYSA-N selenium;zinc Chemical compound [Se]=[Zn] SBIBMFFZSBJNJF-UHFFFAOYSA-N 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000004044 response Effects 0.000 abstract description 3
- 238000013480 data collection Methods 0.000 abstract description 2
- 230000007812 deficiency Effects 0.000 abstract description 2
- 238000003331 infrared imaging Methods 0.000 abstract description 2
- 230000035945 sensitivity Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
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/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
-
- 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/021—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using plane or convex mirrors, parallel phase plates, or particular reflectors
-
- 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
- G01J2003/283—Investigating the spectrum computer-interfaced
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of infrared broadband EO-1 hyperions to calculate imaging device comprising preposition infrared lens, coding templet generation load, reflective spatial light modulator, infrared relay lens group, haplotype spectral radiometer, infrared spectrum are calculated as module, hardware synchronization controller；It also discloses infrared broadband EO-1 hyperion and calculates imaging method.The present invention has the advantage of wider spectral response range and lower system noise using haplotype spectral radiometer, mechanical scanning drawback is not only avoided, reduces data collection capacity, shorten the imaging required time, but also improves spectral resolution and spatial resolution；The present invention not only has the features such as the real-time, quick, accurate of Fourier spectrum radiometer, high sensitivity, also effectively overcomes the deficiencies of mechanical scanning speed existing for conventional elements non-imaged type infrared imaging spectrometer is slow, and data redundancy amount is big.
Description
Technical field
The invention belongs to spectral imaging technology fields, and in particular to a kind of infrared broadband EO-1 hyperion calculate imaging device and
Its method.
Background technology
Infrared spectroradio meter is the instrument and equipment for measuring object in infrared band emission spectrum.By obtaining the red of object
External spectrum radiates, and comes the radiation characteristic and physical attribute of object analysis, such as determines peak strength, the position of object infra-red emission
It sets and shape, analyzes Infrared spectra adsorption characteristic of gas etc.；Therefore, it be widely used in agricultural, environment, geology,
The fields such as oil, security protection, military affairs.
Infrared spectroradio meter can be divided into two class of unit non-imaged type and imaging-type.Unit non-imaged type infrared spectral radiant
Meter acquires the infrared spectral radiant distributed data in a visual field every time as unit of system visual field, and data do not include target
Spatial information；Its outstanding feature is that have very wide spectral response range and lower system noise；But due to being that machinery is swept
Formula is retouched, system stability is required high；And gathered data amount is huge, is unfavorable for transimission and storage.Imaging-type infrared spectrum
The photoelectric sensing module of radiometer is made of a large amount of cell sensor, and what each cell sensor completed that it is observed instantaneous regards
The infrared spectral radiant of field measures, and the infrared spectral radiant for being achieved in spatial decomposition measures, and obtains infrared in system visual field
The space distribution information of radiation；However, imaging-type infrared spectroradio meter needs repeatedly to expose target in each dimension
Light and Multiple-Scan could obtain the complete three-dimensional stereo data of target, cause observation time to extend, and involve great expense.
Recently as the appearance of the spatial light modulation devices such as digital micromirror array (DMD), liquid crystal light valve, it is based on notch
The research of diameter, the staring imaging spectral technique for calculating imaging has become a research hotspot.
Invention content
In view of this, the main purpose of the present invention is to provide a kind of infrared broadband EO-1 hyperion calculate imaging device and its
Method.
In order to achieve the above objectives, the technical proposal of the invention is realized in this way：
The embodiment of the present invention provides a kind of infrared broadband EO-1 hyperion calculating imaging device comprising preposition infrared lens,
Coding templet generates load, reflective spatial light modulator, infrared relay lens group, haplotype spectral radiometer, infrared spectrum
Calculate image-forming module, hardware synchronization controller；
The preposition infrared lens module, for target scene infra-red radiation to be converged to reflective spatial light modulator
On target surface；
The coding templet maker module, it is random for the Bernoulli Jacob for repeatedly measuring and loading to be generated in advance by computer
Coding templet matrix, the value of matrix element are " -1 " and " 1 ", and obedience Bernoulli Jacob's random distribution, and by random coded template
Pass to reflective spatial light modulator；
The reflective spatial light modulator is located on the focal plane of preposition infrared lens, and optical window carries out ZnSe
1-15 μm of infrared band after material secondary encapsulation in coverage goal scene infra-red radiation；Further according to real-time random coded template pair
The random reflected of target scene infra-red radiation completes the coded modulation to target scene；
The infrared relay lens group module, for converging to the modulated light path of the reflective spatial light modulator
Haplotype spectral radiometer；
The haplotype spectral radiometer, spatial distribution for being measured according to the light path and is sent to infrared calculating
Image-forming module；
The infrared calculating image-forming module, for reconstructing target field according to the spatial distribution and calculation code image-forming principle
The profile information of scape different-waveband obtains infrared broadband high spectrum image；
The hardware synchronization controller makes reflective spatial light modulator and haplotype spectrum spoke for passing through timing control
It penetrates meter and realizes synchronous coding modulation and data acquisition.
In said program, the infrared relay lens group module includes speculum and Infrared Lens group, and the speculum is used
Pass through the modulated light path of light field in adjustment；The Infrared Lens group is made of multiple lens that focus and optical axis overlap,
The focal length of lens and clear aperture of the Infrared Lens group are adjusted according to haplotype spectral radiometer Receiver aperture and position, so as to
It is received by haplotype spectral radiometer.
In said program, the reflective spatial light modulator, be specifically used for by circuit control micro-reflector according to add
It carries the random coded template generated and changes rollover states, wherein " -1 " indicates that micro-reflector overturns -12 °, and " 1 " indicates micro- reflection
Mirror overturns+12 °.
The embodiment of the present invention also provides a kind of infrared broadband EO-1 hyperion calculating imaging method, and this method is：It is exerted according to uncle
Sharp random coded template carries out compression sampling to infrared target scene random coded, to echo signal after coding, finally according to weight
The compressed encoding sampled signal that the reconstruct of structure algorithm obtains obtains original infrared high spectrum image.
In said program, this method specifically includes following steps：
Step 1：Assuming that original infrared spectroscopic imaging data cube VP×Q×L, resolution sizes are P × Q, and spectrum is logical
Road number is L, is now expressed as (i, j) a image pixel intensities of first of spectrum channel in cubeThe then infrared spoke of target scene
It penetrates and is represented by vectorTarget scene infra-red radiation converged at by preposition infrared lens reflective
On the target surface of spatial light modulator；
Step 2：The reflective spatial light modulator loads the random matrix g that a resolution ratio is P × Q, passes through change
Its pixel unit on off state, makes infra-red radiation be reflected at random, completes the coded modulation to target scene；
Step 3：Modulated infrared target scene adjusts beam path alignment by infrared relay lens group；
Step 4：Haplotype spectral radiometer acquires the modulation light field after collimation, and often modulation is primary, exports one group of pressure
Contracting sampling spectroscopic data
Step 5：And so on, the reflective spatial light modulator is by D modulation, the haplotype spectral radiance
Collected compression sampling spectroscopic data is counted to be followed successively by：
W=[w1,w2,…wD] (2)
Step 6：Finally, the spectroscopic data that infrared calculating image-forming module obtains haplotype spectral radiometer is reconstructed
Obtain the infrared broadband high spectrum image of target and scene.
In said program, the step 6 is realized especially by following steps：
Step 101：Entire data acquisition is expressed as matrix form：
W=GV (3)
In formula, W=[w1,w2,…wD]TFor the compressed spectrum data obtained after the measurement of D sublinears, G=[g1,
g2,…,gD]TUsed modulation function is measured for D sublinears, V is initial data cube；
Step 7：Sparse transformation is carried out to initial data V, it is made to meet the precondition of compressive sensing theory, its is dilute
It dredges and is expressed as,
V=Ψ θ (4)
In formula, θ is sparse transformation vectors of the original signal V at base Ψ, and θ is that Q × L ties up matrix, and Ψ is that Q × Q dimensions are sparse
Transformation matrix, if having k nonzero element, wherein k in θ<<Q, it is the k rarefaction representations of original signal to claim θ；
Step 102：Then, formula W=GV is expressed as again：
W=GV=G Ψ θ=Θ θ (5)
In formula, Θ be P × Q tie up sensing matrix, by the cubical refutation process of entire spectral image data be considered as one it is right
The convex Optimization Solution problem of original sparse coefficient：
In formula, the l of first item difference between model and measurement data2Norm minimum, Section 2 are the l of reconstruction coefficients1
Norm indicates its sparsity, parameter γ>0 is regulatory factor；
Step 103：The optimum solution of above formula is sought by compressed sensing restructing algorithmThus it is red to obtain target scene for inverting
External spectrum image can similarly acquire the high spectrum image of arbitrary wave band.
In said program, the step 2 is specially：The reflective spatial light modulator load a resolution ratio be P ×
The coding templet of Q, the coding templet are Bernoulli Jacob's random matrix, meet and are limited equidistant property (RIP) i.e.In formula, calculation matrix column vector is all unit length, and signal v is that degree of rarefication is K
Sparse signal, δKBe a value range be (0,1) constant, calculation matrix g tieed up for P × Q, the value of matrix element be -1 or
Person 1, and obeys Bernoulli Jacob's random distribution：
In formula, f isProbability.
Compared with prior art, the present invention has wider spectral response range and relatively low using haplotype spectral radiometer
System noise advantage, not only avoid mechanical scanning drawback, reduce data collection capacities, shorten and be imaged the required time, and
And improve spectral resolution and spatial resolution；The present invention is using calculating imaging theory and unit non-imaged type infrared spectral radiant
The infrared broadband EO-1 hyperion for counting the method development combined calculates the reality that imaging device not only has Fourier spectrum radiometer
When, quickly, accurate, high sensitivity the features such as, also effectively overcome the presence of conventional elements non-imaged type infrared imaging spectrometer
Mechanical scanning speed it is slow, the deficiencies of data redundancy amount is big.
Description of the drawings
Fig. 1 provides a kind of connection block diagram of infrared broadband EO-1 hyperion calculating imaging device for the embodiment of the present invention；
Fig. 2, which for the embodiment of the present invention provides a kind of infrared broadband EO-1 hyperion and calculates the corresponding optical texture of imaging device, to be shown
It is intended to；
Fig. 3 is Bernoulli Jacob's random coded template that reflective spatial light modulator loads in the present invention；
Fig. 4 provides a kind of flow chart of infrared broadband EO-1 hyperion calculating imaging method for the embodiment of the present invention；
Fig. 5 is the infrared high spectrum image obtained using OMP restructing algorithms in the present invention.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
The embodiment of the present invention provides a kind of infrared broadband EO-1 hyperion calculating imaging device, as shown in Figure 1 comprising preposition
Infrared lens, coding templet generate load, reflective spatial light modulator, infrared relay lens group, haplotype spectral radiometer,
Infrared spectrum is calculated as module, hardware synchronization controller；
The preposition infrared lens module, for target scene infra-red radiation to be converged to reflective spatial light modulator
On target surface；
The coding templet maker module, it is random for the Bernoulli Jacob for repeatedly measuring and loading to be generated in advance by computer
Coding templet matrix, the value of matrix element are " -1 " and " 1 ", and obedience Bernoulli Jacob's random distribution, and by random coded template
Pass to reflective spatial light modulator；
The reflective spatial light modulator is located on the focal plane of preposition infrared lens, optical window material secondary
1-15 μm of infrared band after encapsulation in coverage goal scene infra-red radiation；Further according to real-time random coded template to target scene
The random reflected of infra-red radiation completes the coded modulation to target scene；
Specifically, it can be covered after the optical window progress ZnSe material secondary encapsulation of the reflective spatial light modulator red
1-15 μm of wave section.It is generated according to load by the micro-reflector of the circuit controls reflective spatial light modulator such as FPGA and DDC
Random coded template change rollover states, wherein " -1 " indicate micro-reflector overturn -12 °, " 1 " indicate micro-reflector overturning+
12°。
The reflective spatial light modulator can cover infrared broadband after carrying out secondary encapsulation to its optical window, break through
Existing space optical modulator is only used for the limitation of visible light near infrared light spectral imaging technology.
The infrared relay lens group module, for converging to the modulated light path of the reflective spatial light modulator
Haplotype spectral radiometer；
The haplotype spectral radiometer, spatial distribution for being measured according to the light path and is sent to infrared calculating
Image-forming module；
Specifically, the haplotype spectral radiometer is passed through using non-imaged type infrared spectroradio meter acquisition target scene
Modulated spectral information.
The infrared calculating image-forming module, for reconstructing target field according to the spatial distribution and calculation code image-forming principle
The profile information of scape different-waveband obtains infrared broadband high spectrum image；
Specifically, the infrared relay lens group module includes speculum and Infrared Lens group, and the speculum is for adjusting
Warping crosses the modulated light path of light field；The Infrared Lens group is made of multiple lens that focus and optical axis overlap, described
The focal length of lens and clear aperture of Infrared Lens group are adjusted according to haplotype spectral radiometer Receiver aperture and position, so as to coverlet
First type spectral radiometer receives.
The hardware synchronization controller makes reflective spatial light modulator and haplotype spectrum spoke for passing through timing control
It penetrates meter and realizes synchronous coding modulation and data acquisition.
Imaging device, which is calculated, based on the infrared broadband EO-1 hyperion builds Optical devices as shown in Figure 2, it is described preposition
The target scene infra-red radiation of acquisition is emitted on reflective spatial light modulator (DMD) by infrared lens module, the reflection
The light path of formula spatial light modulator reflection is converged to by infrared relay lens group on haplotype spectral radiometer.
The embodiment of the present invention also provides a kind of infrared broadband EO-1 hyperion calculating imaging method, and this method is：It is exerted according to uncle
Sharp random coded template carries out compression sampling to infrared target scene random coded, to echo signal after coding, finally according to weight
The compressed encoding sampled signal that the reconstruct of structure algorithm obtains obtains original infrared high spectrum image.
Especially by following implemented：
(1) original infrared spectroscopic imaging data cube V is setP×Q×L, image size is P × Q, and spectrum channel number is L；
Now (i, j) a image pixel intensities of first of spectrum channel in cube are expressed asThen target scene infra-red radiation is writeable
At the form of vectorIt is converged on reflective spatial light modulator target surface by preposition infrared lens.
The infrared band that can be received on the target surface of the reflective spatial light modulator target surface is rung by infrared lens spectrum
Range is answered to determine.
(2) reflective spatial light modulator meets limited equidistant property (RIP) by loading one, and resolution ratio is P × Q
Coding templet as calculation matrix, i.e.,In formula, calculation matrix column vector is all unit
Length, signal v are the sparse signal that degree of rarefication is K, δKBe a value range be (0,1) constant.Calculation matrix g is tieed up for P × Q,
The value of matrix element is -1 or 1, and obeys Bernoulli Jacob's random distribution：I.e.
In formula, f isProbability makes infra-red radiation be reflected at random, completes by changing its pixel unit on off state
To the coded modulation of target scene；
(3) modulated infrared target scene adjusts beam path alignment by infrared relay lens group；
(4) haplotype spectral radiometer acquires the modulation light field after collimation, and often modulation is primary, and one compression of output is adopted
Sample spectroscopic data
(5) and so on, reflective spatial light modulator is by D modulation, the collected pressure of haplotype spectral radiometer
Contracting sampling spectroscopic data is followed successively by：W=[w1,w2,…wD]；
(6) finally, acquisition is reconstructed in the spectroscopic data that infrared calculating image-forming module obtains haplotype spectral radiometer
The infrared broadband high spectrum image of target and scene.
As shown in figure 3, generating Bernoulli Jacob's random coded mould that 4000 groups of resolution ratio are 128 × 96 in advance using MATLAB
Plate, value are -1 or 1, and " -1 " and " 1 " obeys Bernoulli Jacob's random distribution.Due to spatial light modulator resolution ratio be 1024 ×
768, therefore make a pixel of 8 × 8 pixel compositions " super pixel " as coding templet on spatial light modulator minute surface, lead to
It crosses " -1 " and indicates that micro-reflector overturns -12 °, " 1 " indicates that micro-reflector overturns+12 °, changes spatial light modulator micro-reflector
Rollover states make the light of different location enter subsequent optical system, realize and are modulated to target scene.
As shown in figure 4, a kind of infrared broadband EO-1 hyperion of the embodiment of the present invention calculates imaging method, this method specifically includes
Following steps：
(1) original infrared spectroscopic imaging data cube V is setP×Q×L, image size is P × Q, and spectrum channel number is L.
Now (i, j) a image pixel intensities of first of spectrum channel in cube are expressed asThen target scene infra-red radiation is writeable
At the form of vectorIt is converged on reflective spatial light modulator target surface by preposition infrared lens.
(2) coding templet that one resolution ratio of reflective spatial light modulator load is P × Q is measured as calculation matrix
Matrix, which meets, is limited equidistant property (RIP) i.e.In formula, calculation matrix column vector is all that unit is long
Degree, signal v are the sparse signal that degree of rarefication is K, δKBe a value range be (0,1) constant.Calculation matrix g, square are tieed up for P × Q
The value of array element element is -1 or 1, and obeys independent Bernoulli Jacob's random distribution：I.e.
In formula, f isProbability makes infra-red radiation be reflected at random, completes by changing its pixel unit on off state
To the coded modulation of target scene；
(3) modulated infrared target scene adjusts beam path alignment by infrared relay lens group；
(4) haplotype spectral radiometer acquires the modulation light field after collimation, and often modulation is primary, and one compression of output is adopted
Sample spectroscopic data
(5) and so on, reflective spatial light modulator is by D modulation, the collected pressure of haplotype spectral radiometer
Contracting sampling spectroscopic data is followed successively by：W=[w1,w2,…wD]；
(6) entire data acquisition is expressed as matrix form：
W=GV
In formula, W=[w1,w2,…wD]TFor the compressed spectrum data obtained after the measurement of D sublinears, G=[g1,
g2,…,gD]TUsed modulation function is measured for D sublinears, V is initial data cube；
(7) it is ill signal V theoretically to be recovered from measured value W, but the sparsity that signal has provides thus
May, compressive sensing theory frame gives in the presence of the necessary and sufficient condition for determining solution, that is, limits isometry (RIP).It is sparse or
Compressible signal can be by less than the linear measurement of traditional sampling number by Accurate Reconstruction.To ensure data cube
Exact inversion needs to carry out sparse transformation to initial data V, it is made to meet the precondition of compressive sensing theory, its is sparse
It is expressed as such as formula：
V=Ψ θ
In formula, θ is sparse transformation vectors of the original signal V at base Ψ, and θ is that Q × L ties up matrix, and Ψ is that Q × Q dimensions are sparse
Transformation matrix.If having k nonzero element, wherein k in θ<<N, it is the k rarefaction representations of original signal to claim θ；
(8) then, formula W=GV can be expressed as again：
W=GV=G Ψ θ=Θ θ
In formula, Θ is that P × Q ties up sensing matrix, is based on above-mentioned model, can be by the cubical inverting of entire spectral image data
Process is considered as a convex Optimization Solution problem to original sparse coefficient：
In formula, the l of first item difference between model and measurement data2Norm minimum, Section 2 are the l of reconstruction coefficients1
Norm indicates its sparsity, parameter γ>0 is regulatory factor；
(9) optimum solution of above formula is sought by compressed sensing restructing algorithm such as OMP algorithmsThus inverting obtains target field
Scape infrared spectroscopic imaging can similarly acquire the high spectrum image of arbitrary wave band.
As shown in figure 5, the present invention is 32.5% completion to target by continuously loading n group random coded templates, sample rate
The compression sampling of scene, the broadband infrared high spectrum image obtained using image reconstruction algorithm.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the scope of the present invention.
Claims (7)
1. a kind of infrared broadband EO-1 hyperion calculates imaging device, which is characterized in that it includes preposition infrared lens, coding templet
Generation load, reflective spatial light modulator, infrared relay lens group, haplotype spectral radiometer, infrared spectrum are calculated as picture
Module, hardware synchronization controller；
The preposition infrared lens module, the target surface for target scene infra-red radiation to be converged to reflective spatial light modulator
On；
The coding templet maker module, for the Bernoulli Jacob's random coded for repeatedly measuring load to be generated in advance by computer
The value of pattern matrix, matrix element is " -1 " and " 1 ", and obeys Bernoulli Jacob's random distribution, and random coded template is transmitted
To reflective spatial light modulator；
The reflective spatial light modulator is located on the focal plane of preposition infrared lens, and optical window carries out ZnSe materials
1-15 μm of infrared band after secondary encapsulation in coverage goal scene infra-red radiation；Further according to real-time random coded template to target
The random reflected of scene infra-red radiation completes the coded modulation to target scene；
The infrared relay lens group module, for the modulated light path of the reflective spatial light modulator to be converged to unit
Type spectral radiometer；
The haplotype spectral radiometer, spatial distribution for being measured according to the light path and is sent to infrared calculating imaging
Module；
The infrared calculating image-forming module, for reconstructing target scene not according to the spatial distribution and calculation code image-forming principle
Profile information with wave band obtains infrared broadband high spectrum image；
The hardware synchronization controller makes reflective spatial light modulator and haplotype spectral radiometer for passing through timing control
Realize synchronous coding modulation and data acquisition.
2. infrared broadband EO-1 hyperion according to claim 1 calculates imaging device, which is characterized in that the infrared relaying
Lens group module includes speculum and Infrared Lens group, and the speculum passes through the modulated light path of light field for adjusting；It is described
Infrared Lens group is made of multiple lens that focus and optical axis overlap, the focal length of lens and light admission port of the Infrared Lens group
Diameter is adjusted according to haplotype spectral radiometer Receiver aperture and position, to be received by haplotype spectral radiometer.
3. infrared broadband EO-1 hyperion according to claim 1 or 2 calculates imaging device, which is characterized in that the reflection
Formula spatial light modulator is specifically used for changing overturning according to the random coded template that load generates by circuit control micro-reflector
State, wherein " -1 " indicates that micro-reflector overturns -12 °, and " 1 " indicates that micro-reflector overturns+12 °.
4. a kind of infrared broadband EO-1 hyperion calculates imaging method, which is characterized in that this method is：According to Bernoulli Jacob's random coded
Template carries out compression sampling to infrared target scene random coded, to echo signal after coding, is finally reconstructed according to restructing algorithm
The compressed encoding sampled signal of acquisition obtains original infrared high spectrum image.
5. infrared broadband EO-1 hyperion according to claim 4 calculates imaging method, which is characterized in that this method is specifically wrapped
Include following steps：
Step 1：Assuming that original infrared spectroscopic imaging data cube VP×Q×L, resolution sizes are P × Q, spectrum channel number
For L, now (i, j) a image pixel intensities of first of spectrum channel in cube are expressed asThen target scene infra-red radiation
It is represented by vectorTarget scene infra-red radiation is converged at into Reflective spatial by preposition infrared lens
On the target surface of optical modulator；
Step 2：The reflective spatial light modulator loads the random matrix g that a resolution ratio is P × Q, by changing its picture
Plain unit switch state, makes infra-red radiation be reflected at random, completes the coded modulation to target scene；
Step 3：Modulated infrared target scene adjusts beam path alignment by infrared relay lens group；
Step 4：Haplotype spectral radiometer acquires the modulation light field after collimation, and often modulation is primary, and one group of compression of output is adopted
Sample spectroscopic data
Step 5：And so on, the reflective spatial light modulator is adopted by D modulation, the haplotype spectral radiometer
The compression sampling spectroscopic data collected is followed successively by：
W=[w1,w2,…wD] (2)
Step 6：Finally, acquisition is reconstructed in the spectroscopic data that infrared calculating image-forming module obtains haplotype spectral radiometer
The infrared broadband high spectrum image of target and scene.
6. infrared broadband EO-1 hyperion according to claim 5 calculates imaging method, which is characterized in that the step 6 tool
Body is realized by following steps：
Step 101：Entire data acquisition is expressed as matrix form：
W=GV (3)
In formula, W=[w1,w2,…wD]TFor the compressed spectrum data obtained after the measurement of D sublinears, G=[g1,g2,…,gD]T
Used modulation function is measured for D sublinears, V is initial data cube；
Step 7：Sparse transformation is carried out to initial data V, so that it is met the precondition of compressive sensing theory, by its sparse table
It is shown as,
V=Ψ θ (4)
In formula, θ is sparse transformation vectors of the original signal V at base Ψ, and θ is that Q × L ties up matrix, and Ψ is that Q × Q ties up sparse transformation
Matrix, if having k nonzero element, wherein k in θ<<Q, it is the k rarefaction representations of original signal to claim θ；
Step 102：Then, formula W=GV is expressed as again：
W=GV=G Ψ θ=Θ θ (5)
In formula, Θ is that P × Q ties up sensing matrix, and the cubical refutation process of entire spectral image data is considered as one to original
The convex Optimization Solution problem of sparse coefficient：
In formula, the l of first item difference between model and measurement data2Norm minimum, Section 2 are the l of reconstruction coefficients1Model
Number, indicates its sparsity, parameter γ>0 is regulatory factor；
Step 103：The optimum solution of above formula is sought by compressed sensing restructing algorithmThus inverting obtains target scene infrared light
Spectrogram picture can similarly acquire the high spectrum image of arbitrary wave band.
7. infrared broadband EO-1 hyperion according to claim 5 or 6 calculates imaging method, which is characterized in that the step
Two are specially：The reflective spatial light modulator loads the coding templet that a resolution ratio is P × Q, which is primary
Sharp random matrix is exerted, meets and is limited equidistant property (RIP) i.e.In formula, calculation matrix row
Vector is all unit length, and signal v is the sparse signal that degree of rarefication is K, δKBe a value range be (0,1) constant, for P
× Q ties up calculation matrix g, and the value of matrix element is -1 or 1, and obeys Bernoulli Jacob's random distribution：
In formula, f isProbability.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810050205.9A CN108663118B (en) | 2018-01-18 | 2018-01-18 | Infrared broadband hyperspectral calculation imaging device and method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810050205.9A CN108663118B (en) | 2018-01-18 | 2018-01-18 | Infrared broadband hyperspectral calculation imaging device and method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108663118A true CN108663118A (en) | 2018-10-16 |
| CN108663118B CN108663118B (en) | 2020-10-16 |
Family
ID=63784850
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810050205.9A Active CN108663118B (en) | 2018-01-18 | 2018-01-18 | Infrared broadband hyperspectral calculation imaging device and method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108663118B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109828285A (en) * | 2019-01-07 | 2019-05-31 | 北京理工大学 | A kind of two waveband time domain data compression perception high speed imaging method and device |
| CN113008370A (en) * | 2021-03-02 | 2021-06-22 | 北京理工大学 | Three-dimensional self-adaptive compression reconstruction method based on liquid crystal hyperspectral calculation imaging system |
| CN114442195A (en) * | 2022-01-28 | 2022-05-06 | 苏州浪潮智能科技有限公司 | Short wave infrared image generation device, method, computer equipment and storage medium |
| CN114882235A (en) * | 2022-05-12 | 2022-08-09 | 北京中星时代科技有限公司 | Infrared image target automatic identification system and method based on image data characteristics |
-
2018
- 2018-01-18 CN CN201810050205.9A patent/CN108663118B/en active Active
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109828285A (en) * | 2019-01-07 | 2019-05-31 | 北京理工大学 | A kind of two waveband time domain data compression perception high speed imaging method and device |
| CN113008370A (en) * | 2021-03-02 | 2021-06-22 | 北京理工大学 | Three-dimensional self-adaptive compression reconstruction method based on liquid crystal hyperspectral calculation imaging system |
| CN114442195A (en) * | 2022-01-28 | 2022-05-06 | 苏州浪潮智能科技有限公司 | Short wave infrared image generation device, method, computer equipment and storage medium |
| CN114442195B (en) * | 2022-01-28 | 2023-11-17 | 苏州浪潮智能科技有限公司 | Short wave infrared image generation device, method, computer equipment and storage medium |
| CN114882235A (en) * | 2022-05-12 | 2022-08-09 | 北京中星时代科技有限公司 | Infrared image target automatic identification system and method based on image data characteristics |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108663118B (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108663118A (en) | Infrared broadband EO-1 hyperion calculates imaging devices and methods therefor | |
| CN103472457B (en) | Three-dimensional imaging system and method for calculating correlation flight time by means of sparse aperture compression | |
| CN103472455B (en) | Four-dimensional spectral imaging system and method for calculating correlation flight time by means of sparse aperture compression | |
| CN106772430B (en) | The single pixel photon counting 3-D imaging system and method approached based on multiresolution wavelet | |
| CN103472456B (en) | Active imaging system and method based on sparse aperture compressing calculation correlation | |
| CN104121990B (en) | Compressed sensing broadband Hyperspectral imager based on random grating | |
| CN110471082A (en) | Single pixel laser calculates imaging device and method | |
| Jin et al. | Polarimetric scattering and SAR information retrieval | |
| CN103913229B (en) | Coding template multi-target active imaging spectrum system and method | |
| CN109520969B (en) | Distributed scattering imaging method based on atmospheric medium self-modulation | |
| CN105675146B (en) | Dual wavelength three-dimensional temperature field imaging device, system and method based on compressed sensing | |
| CN106993121B (en) | A kind of compressed sensing based high spectrum image acquisition imaging system and control method | |
| CN106053433B (en) | Laser-induced spectral analysis method and apparatus are perceived based on light modulation compression dimensionality reduction | |
| CN104019898B (en) | Ultrasensitive spectral imaging astronomical telescope and astronomical spectral imaging method | |
| CN109297925A (en) | A kind of Terahertz high-resolution fast imaging device based on splits' positions perception | |
| US20230125131A1 (en) | Ultrafast light field tomography | |
| EP4357741A1 (en) | Method and system for measuring spectral reflectivity | |
| CN108007574A (en) | The fast illuminated image spectrum linear polarization detection device of resolution ratio adjustable type and method | |
| CN103954357B (en) | The acquisition methods of compressed spectrum imaging system calculation matrix | |
| CN105737992B (en) | Dual wavelength temperature field imaging device, system and method based on compressed sensing | |
| CN114659634A (en) | Miniature snapshot type compressed spectrum imaging detection device and detection method | |
| CN112987021B (en) | Structured light three-dimensional imaging system and method integrating time-of-flight method and structured light method | |
| CN106872977A (en) | A kind of chromatography SAR three-D imaging methods based on the weak orthogonal matching pursuit of segmentation | |
| CN117848503A (en) | Multi-spectral-band polarization intelligent detection device and method for three-dimensional restoration of high-precision target | |
| CN114264370A (en) | Compressed sensing computed tomography spectrometer system and imaging method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |