CN107315176A - Imaging device and method under the conditions of a kind of powerful gas scattering - Google Patents
Imaging device and method under the conditions of a kind of powerful gas scattering Download PDFInfo
- Publication number
- CN107315176A CN107315176A CN201710517013.XA CN201710517013A CN107315176A CN 107315176 A CN107315176 A CN 107315176A CN 201710517013 A CN201710517013 A CN 201710517013A CN 107315176 A CN107315176 A CN 107315176A
- Authority
- CN
- China
- Prior art keywords
- image
- mrow
- mtd
- frequency
- function
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/003—Bistatic lidar systems; Multistatic lidar systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/481—Constructional features, e.g. arrangements of optical elements
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses the imaging device under the conditions of a kind of powerful gas scattering, the imaging device is used to be imaged object through scattering medium, including laser, spatial light modulator, first lens, second lens and imaging sensor, wherein, the spatial light modulator includes multiple turnover micro mirrors, the spatial light modulator is arranged on the laser optical path that the laser is projected, after first lens transmission and penetrated again on the object through the scattering medium after micro mirror reflection of the laser optical path through the spatial light modulator, then the laser optical path reflects through the object and penetrated again after second lens transmission in described image sensor to be imaged to the object after passing through the scattering medium.The invention also discloses the imaging method under the conditions of a kind of powerful gas scattering.Imaging device and method under the conditions of powerful gas scattering proposed by the present invention, substantially increase the image quality under the conditions of the scattering of powerful gas.
Description
Technical field
The present invention relates to the imaging device under the conditions of picture imaging techniques field, more particularly to a kind of powerful gas scattering and side
Method.
Background technology
The industries such as aviation, navigation and highway communication have widely to the imaging in the strong scattering mediums such as haze, misty rain
Demand.Imaging method in existing common strong scattering medium is imaged for near-infrared laser active illumination, utilizes specific wavelength
Light realizes preferable imaging effect to the penetrability of atmospheric scattering, and this method is substantially exactly to use near-infrared laser as light source
Traditional optical imaging concept, but active illumination imaging method powerful gas scattering under the conditions of, image quality is substantially reduced；
The image quality under the conditions of the scattering of powerful gas is improved, is the direction of those skilled in the art's effort.
The disclosure of background above technology contents is only used for design and the technical scheme that auxiliary understands the present invention, and it is not necessarily
Belong to the prior art of present patent application, without tangible proof show the above present patent application the applying date
In the case of disclosed, above-mentioned background technology should not be taken to evaluate the novelty and creativeness of the application.
The content of the invention
In order to improve the image quality under the conditions of the scattering of powerful gas, the present invention is proposed under the conditions of a kind of powerful gas scattering
Imaging device and method.
In order to achieve the above object, the present invention uses following technical scheme：
The invention discloses the imaging device under the conditions of a kind of powerful gas scattering, the imaging device is used to be situated between through scattering
Confrontation object is imaged, including laser, spatial light modulator, the first lens, the second lens and imaging sensor, its
In, the spatial light modulator includes multiple turnover micro mirrors, and the spatial light modulator is arranged on the laser and projected
Laser optical path on, the laser optical path through the spatial light modulator the micro mirror reflection after again pass through first lens
Penetrated after transmission and through the scattering medium on the object, then the laser optical path reflects and saturating through the object
Cross after the scattering medium and penetrated again after second lens transmission in described image sensor to enter to the object
Row imaging.
Preferably, the laser uses wavelength for 720~904nm LASER Light Source.
Preferably, the spatial light modulator includes the turnover micro mirror of M × N number of matrix arrangement.
The invention also discloses the imaging method under the conditions of a kind of powerful gas scattering, carried out into using above-mentioned imaging device
Picture, comprises the following steps：
S1：The calculation matrix for the M × N for being all 1 by one is input to the spatial light modulator, in described image sensor
The first image of upper generation, wherein 1 in the calculation matrix represents to turn over the corresponding micro mirror in the spatial light modulator
The laser optical path for projecting the laser is gone to reflex on the object；
S2：One group of calculation matrix comprising 0 and 1 M × N is input to the spatial light modulator, by described in the group
The intensity signal that calculation matrix and corresponding described image sensor are received, reduction the second image of generation, wherein the measurement square
0 in battle array represents to overturn the corresponding micro mirror in the spatial light modulator to the laser light not projected the laser
Road is reflexed on the object；
S3：Described first image is weighted with second image by the way of frequency domain weighting and is added, generation is most
The synthetic image of the whole object.
Preferably, step S1 also includes, and processing is filtered to described first image, generates filtered first image,
Described first image in step S3 is filtered described first image.
Preferably, the light intensity received in step S2 by the group calculation matrix and corresponding described image sensor is believed
Breath, reduction the second image of generation is specifically included：Using following calculation formula：
Y=Φ x
Wherein, x is the image raw information of one-dimensional, and y is that the reflected light that m sampling described image sensor is received is total
Intensity, Φ is calculation matrix collection, and m is the matrix quantity of calculation matrix described in one group, n=M × N；Can root by above-mentioned formula
According to Φ and y reconstruct generation x, i.e. reduction generation second image.
Preferably, wherein the algorithm reconstructed uses OMP algorithms.
Preferably, step S3 is specifically included：
S31：Described first image and second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as described first image and described
The weighting function of second image；
S33：By the described first two-dimentional piecewise function and the second two-dimentional piecewise function respectively with described first image and
The frequency-domain function that second image is obtained by Fourier transformation is multiplied, and is then added, obtains comprehensive frequency-domain function, then enter
Row inversefouriertransform, that is, generate the synthetic image of the final object；
Wherein in step s 32：
Described first two-dimentional piecewise function w1Relation with frequency f is as follows：
Described second two-dimentional piecewise function w2Relation with frequency f is as follows：
Wherein, F is highest frequency.
Preferably, step S3 is specifically included：
S31：Described first image and second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimensional Gaussian function and the second two-dimensional Gaussian function respectively as described first image and described
The weighting function of second image, wherein first two-dimensional Gaussian function and second two-dimensional Gaussian function are respectively through normalizing
Change, and first two-dimensional Gaussian function and the second two-dimensional Gaussian function sum are 1；
S33：By first two-dimensional Gaussian function and second two-dimensional Gaussian function respectively with described first image and
The frequency-domain function that second image is obtained by Fourier transformation is multiplied, and is then added, obtains comprehensive frequency-domain function, then enter
Row inversefouriertransform, that is, generate the synthetic image of the final object.
Preferably, in step S32：
When frequency is less than first predetermined value in first two-dimensional Gaussian function, corresponding weight is 0, and frequency is more than second
During predetermined value, corresponding weight is 1, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, corresponding
Weight is bigger；
When frequency is less than first predetermined value in second two-dimensional Gaussian function, corresponding weight is 1, and frequency is more than second
During predetermined value, corresponding weight is 0, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, corresponding
Weight is smaller.
Compared with prior art, the beneficial effects of the present invention are：Under the conditions of powerful gas scattering proposed by the present invention into
As device can realize the imaging mode of two kinds of different principles, including active illumination imaging method and compressed sensing ghost imaging simultaneously
Method so that the first image and the use that can obtain obtaining using active illumination imaging method simultaneously by the imaging device
The second image that compressed sensing ghost imaging method is obtained, so as to which the first image and the second image further are carried out into General Office
Reason, to obtain preferably synthetic image, so as to substantially increase the image quality under the conditions of the scattering of powerful gas.
In further scheme, laser uses wavelength for the LASER Light Source of 720nm~904 so that laser is sent
Laser optical path there is more preferable penetrability to the scattering medium in air, and keep will not occurring diffraction effect.Spatial light is adjusted
Device processed includes the turnover micro mirror of M × N number of matrix arrangement, so as to which M × N calculation matrix is input into space light modulation
Device, plays a part of modulated light source, by the way that one group of calculation matrix randomly generated is carried out into the corresponding image of Self -adaptive second,
Reduce the sampling number for producing the second image.
In further scheme, the imaging with reference to active illumination imaging method and the terrible imaging method of compressed sensing is special
Property, the present invention in the first image and the weight letter of the second image can be used as by Gaussian function or specific piecewise function
Number, the method for line frequency domain weighting summation of going forward side by side obtains final synthetic image and is superior to the first image and the second image.
Brief description of the drawings
Fig. 1 is the schematic diagram of the imaging device under the conditions of the powerful gas scattering of the preferred embodiment of the present invention；
Fig. 2 a are the spectrograms of the artwork of object；
Fig. 2 b and Fig. 2 c are the spectrograms of the first image and the second image under low scattering coefficient；
Fig. 2 d and Fig. 2 e are the spectrograms of the first image and the second image under high scattering coefficient；
Fig. 3 is the schematic diagram of the second two-dimentional piecewise function in some embodiments of the invention；
Fig. 4 a and Fig. 4 b are the signal for the synthetic image for handling Gaussian function and piecewise function as weighting function respectively
Figure；
Fig. 5 a are the schematic diagrames of the first two-dimensional Gaussian function of the embodiment of the present invention one；
Fig. 5 b are the schematic diagrames of the second two-dimensional Gaussian function of the embodiment of the present invention one；
Fig. 6 a are the results that are multiplied with the frequency-domain function of the first image of the first two-dimensional Gaussian function of the embodiment of the present invention one
Schematic diagram；
Fig. 6 b are that the second two-dimensional Gaussian function of inventive embodiments one shows with the result that the frequency-domain function of the second image is multiplied
It is intended to；
Fig. 6 c are Fig. 6 a and Fig. 6 b result being added；
Fig. 7 a are the schematic diagrames for the first image that the embodiment of the present invention one is obtained；
Fig. 7 b are the schematic diagrames for the second image that the embodiment of the present invention one is obtained；
Fig. 7 c are the schematic diagrames for the synthetic image that the embodiment of the present invention one is obtained；
Fig. 8 a are the schematic diagrames of the artwork image of the object of the embodiment of the present invention two；
Fig. 8 b are the schematic diagrames for the first image that the embodiment of the present invention two is obtained；
Fig. 8 c are the schematic diagrames for filtered first image that Fig. 8 b are obtained by gaussian filtering；
Fig. 8 d are the schematic diagrames for the second image that the embodiment of the present invention two is obtained；
Fig. 8 e are the schematic diagrames for the synthetic image that the embodiment of the present invention two is obtained；
Fig. 9 a are the schematic diagrames of the artwork image of the object of the embodiment of the present invention three；
Fig. 9 b are the schematic diagrames for the first image that the embodiment of the present invention three is obtained；
Fig. 9 c are the schematic diagrames for filtered first image that Fig. 9 b are obtained by gaussian filtering；
Fig. 9 d are the schematic diagrames for the second image that the embodiment of the present invention three is obtained；
Fig. 9 e are the schematic diagrames for the synthetic image that the embodiment of the present invention three is obtained.
Embodiment
Below against accompanying drawing and with reference to preferred embodiment the invention will be further described.
As shown in figure 1, the imaging device under the conditions of the powerful gas scattering of the preferred embodiment of the present invention includes laser 10, sky
Between optical modulator 20, the first lens 30, the second lens 40 and imaging sensor 50, object 60 is carried out by the imaging device
Imaging, wherein there is scattering medium 70 between the imaging device and object 60.The primary structure of the wherein imaging device is：
Spatial light modulator 20 includes multiple turnover micro mirrors, and spatial light modulator 20 is arranged on the laser optical path of the injection of laser 10
On, after micro mirror reflection of the laser optical path through spatial light modulator 20 again after the transmission of the first lens 30 and through scattering medium 70
Penetrate on object 60, then laser optical path is again after object 60 reflects and passes through scattering medium 70 by the second lens 40
Penetrate on imaging sensor 50 to be imaged object 60 after transmission.Wherein, laser 10 is used in some embodiments
Wavelength is 720~904nm LASER Light Source, and spatial light modulator 20 includes the turnover micro mirror of M × N number of matrix arrangement.
In the specific embodiment of the invention, the laser 10 of the imaging device uses wavelength for 808nm near-infrared laser
Light source, has relatively good penetrability to the misty rain in air, and spatial light modulator 20 includes turning over for M × N number of matrix arrangement
The micro mirror turned.Object is imaged by the imaging device, comprised the following steps：
S1：The calculation matrix for the M × N for being all 1 by one is input to spatial light modulator 20, raw on imaging sensor 50
Represent to overturn corresponding micro mirror in spatial light modulator 20 to by laser 10 into the first image, 1 wherein in calculation matrix
The laser optical path of injection is reflexed on object 60；
Now, spatial light modulator 20 reflects all light, and whole light path is exactly a laser active illumination imaging optical path,
What is received in the image planes of imaging sensor 50 is exactly the two dimensional image of object；
In certain embodiments, processing is also filtered to the first image, filtered first image is generated, wherein filtering
Processing can use gaussian filtering method.
S2：One group of calculation matrix (can randomly generate) comprising 0 and 1 M × N is input to spatial light modulator
20, the intensity signal received by this group of calculation matrix and corresponding imaging sensor 50, reduction the second image of generation, wherein surveying
0 expression in moment matrix overturns corresponding micro mirror in spatial light modulator 20 anti-to the laser optical path not projected laser 1
It is mapped on object 60；
Wherein, in the present embodiment, spatial light modulator 20 is made up of M × N number of turnover micro mirror, specific by inputting
Calculation matrix, some micromirrors on its surface can be allowed to overturn so that the light of particular spatial location could be reflected, realize to light
The modulation in source；Whether overturning for micro mirror array determines that the just no of the region is reflected to object 60, and then determines target
Whether the corresponding region on the surface of thing 60 is illuminated, namely each calculation matrix is actual has corresponded to the area that body surface is illuminated
Domain；The intensity signal received by multiple calculation matrix in one group and corresponding imaging sensor, reduction the second image of generation, tool
Body uses following calculation formula：
Y=Φ x
Wherein, x is the image raw information of one-dimensional, and y is that the reflected light that m sampled images sensor 50 is received is always strong
Degree, Φ is calculation matrix collection, and m is every a line that the matrix quantity of one group of calculation matrix, n=M × N, namely calculation matrix are concentrated
Correspond to one group of coding (one calculation matrix of correspondence) of once sampling spatial light modulator；Can basis by above-mentioned formula
Φ and y reconstruct generation x, i.e. the second image of reduction generation；The algorithm wherein reconstructed can use OMP algorithm (orthogonal matching pursuits
Algorithm).
S3：The first image and the second image are weighted addition by the way of frequency domain weighting, final target is generated
The synthetic image of thing.
Pass through the spectrogram to active illumination imaging method under the conditions of different scattering coefficients and the terrible imaging method of compressed sensing
The spectrogram with the artwork of object is made comparisons respectively, and (abscissa is frequency to the spectrogram of artwork, and ordinate is as shown in Figure 2 a
The amplitude of frequency-domain function after Fourier transformation), under low scattering coefficient (3.5), active illumination imaging method obtain first
As shown in Figure 2 b, the spectrogram for the second image that compressed sensing ghost imaging method is obtained leads to the spectrogram of image as shown in Figure 2 c
Cross compare it can be seen that whole wave band be nearly all active illumination imaging method advantage it is interval, highest frequency 1/10th with
Upper (5-45) this interval may be considered absolute interval；Under high scattering coefficient (6.5), active illumination imaging method is obtained
The spectrogram of first image as shown in Figure 2 d, spectrogram such as Fig. 2 e institutes of the second image that compressed sensing ghost imaging method is obtained
Show, by comparing it can be seen that still should be two on the basis of active illumination imaging results, but in highest frequency in high band
Less than/10th (in a figures 2.5 within), compressed sensing ghost imaging method has shown advantage.
Therefore according to the characteristic of active illumination imaging method and the terrible imaging method of compressed sensing, in certain embodiments,
Can be using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as the first image and the weight letter of the second image
Count to carry out frequency domain weighting, it is smoothed in intermediate region in both absolute predominance intervals directly using the frequency spectrum of advantage method
Cross, it is specific as follows：
First two-dimentional piecewise function w1Relation with frequency f is as follows：
Second two-dimentional piecewise function w2Relation with frequency f is as follows, as shown in Figure 3：
Wherein, F is highest frequency, and the value depends on the size of the image of object, and actual value is the image of object
The half of catercorner length.
Further, step S3 can specifically include：
S31：First image and the second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as the first image and the second image
Weighting function；
S33：First two-dimentional piecewise function and the second two-dimentional piecewise function are passed through with the first image and the second image respectively
The frequency-domain function that Fourier transformation is obtained is multiplied, and is then added, the frequency-domain function of obtained synthesis, then carry out anti-Fourier's change
Change, that is, generate the synthetic image of final object.
The function is can be seen that similar to Gaussian function from the schematic diagram of the piecewise function shown in Fig. 3, therefore, in this hair
In other bright embodiments, also frequency domain weighting as weighting function can be carried out using Gaussian function, as shown in figures 4 a and 4b,
Piecewise function and Gaussian function are compared as weighting function processing, Fig. 4 a are to pass through piecewise function integrated treatment
Synthetic image, SSIM is 0.76053, Fig. 4 b for by the synthetic image of Gaussian function integrated treatment, SSIM is 0.78426, can
To find out that the effect that piecewise function integrated treatment and Gaussian function are integrated almost is not different, namely pass through piecewise function and Gauss
Function has obtained preferably effect as weighting function, or even the effect of Gaussian function can be more preferably.
Therefore, in another embodiment, the first two-dimentional piecewise function and the second two dimension segmentation letter in above-mentioned steps S32
The first two-dimensional Gaussian function can also be respectively adopted in number and the second two-dimensional Gaussian function comes value, i.e. step S3 and can also specifically wrapped
Include：
S31：First image and the second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimensional Gaussian function and the second two-dimensional Gaussian function respectively as the first image and the second image
Weighting function, wherein the first two-dimensional Gaussian function and the second two-dimensional Gaussian function are respectively through normalization, and the first two dimension is high
This function and the second two-dimensional Gaussian function sum are 1；
S33：First two-dimensional Gaussian function and the second two-dimensional Gaussian function are passed through with the first image and the second image respectively
The frequency-domain function that Fourier transformation is obtained is multiplied, and is then added, the frequency-domain function of obtained synthesis, then carry out anti-Fourier's change
Change, that is, generate the synthetic image of final object.
Wherein, when frequency is less than first predetermined value in the first two-dimensional Gaussian function, corresponding weight is 0, and frequency is more than the
During two predetermined values, corresponding weight is 1, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, correspondence
Weight it is bigger；
When frequency is less than first predetermined value in second two-dimensional Gaussian function, corresponding weight is 1, and frequency is more than second and made a reservation for
During value, corresponding weight is 0, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, corresponding weight
It is smaller.
Embodiment one：
The schematic diagram of first two-dimensional Gaussian function as shown in Figure 5 a, schematic diagram such as Fig. 5 b institutes of the second two-dimensional Gaussian function
Show, that both are added and for 1, both pass through the frequency-domain function phase that Fourier transformation is obtained with the first image and the second image respectively
Multiply (frequency-domain function multiplied result such as Fig. 6 a institutes that the first two-dimensional Gaussian function and the first image are obtained by Fourier transformation
Show, frequency-domain function multiplied result such as Fig. 6 b institutes that the second two-dimensional Gaussian function and the second image are obtained by Fourier transformation
Show), then it is added, obtains comprehensive frequency-domain function as fig. 6 c, then by returning Fourier transformation, that is, generates final target
The synthetic image of thing as shown in Figure 7 c, wherein by step S1 and S2 the first image respectively obtained and the second image respectively as scheme
Shown in 7a and Fig. 7 b, Fig. 7 c are made comparisons with Fig. 7 a and Fig. 7 b respectively, it can be seen that the effect of synthetic image is better than the first image
With the second image.
Embodiment two：
The artwork image of object as shown in Figure 8 a, scattering medium forward scattering coefficient bd be 4.5 under conditions of, with
Structural similarity (SSIM) as picture quality evaluation criterion.The first image obtained according to step S1 as shown in Figure 8 b,
SSIM is 0.7564, and as shown in Figure 8 c, SSIM is 0.89081 to the first image after gaussian filtering；Obtained according to step S2
The second image as shown in figure 8d, SSIM is 0.60799；According to step S3, using such as Fig. 5 a and Fig. 5 b the first dimensional Gaussian
Function and the second two-dimensional Gaussian function carry out frequency domain weighting with the first image and the second image respectively and are added, and obtain synthetic image such as
Shown in Fig. 8 e, SSIM is 0.88449.
Embodiment three：
The artwork image of object as illustrated in fig. 9, scattering medium forward scattering coefficient bd be 5.5 under conditions of, with
Structural similarity (SSIM) as picture quality evaluation criterion.The first image obtained according to step S1 as shown in figure 9b,
SSIM is 0.20654, and as is shown in fig. 9 c, SSIM is 0.43112 to the first image after gaussian filtering；Obtained according to step S2
The second image as shown in figure 9d, SSIM is 0.37954；According to step S3, using such as Fig. 5 a and Fig. 5 b the first dimensional Gaussian
Function and the second two-dimensional Gaussian function carry out frequency domain weighting with the first image and the second image respectively and are added, and obtain synthetic image such as
Shown in Fig. 9 e, SSIM is 0.54104.
Embodiment two and embodiment three are the synthesis quality reconstruction under different scattering coefficients respectively, it can be seen that passed through
With reference to two kinds of imaging methods, when two methods gap is greatly different can closely preferable imaging effect, both bad
When, effect can be obtained while better than both images.
In the present invention, the first image is the image obtained according to active illumination imaging method, and the second image is according to pressure
Contracting perceives the image that terrible imaging method is obtained, wherein, active illumination imaging method and the terrible imaging method of compressed sensing are in principle
There is very big difference, while also variant on imaging characteristic.Found by the research of applicant, compressed sensing ghost imaging method pair
Random noise is insensitive, but more sensitive to overall fluctuation of optical field intensity, and active illumination imaging method is on the contrary；Compressed sensing
Terrible imaging method can preferably retain the low-frequency information of image, and active illumination imaging side rule is in certain scattering coefficient scope
It is interior to retain medium-high frequency information well.In order to suppress that also active illumination is imaged in noise, the preferred embodiments of the present invention
The image that method is obtained is filtered processing, such as gaussian filtering.
In current picture imaging techniques field, research mainly concentrates on the restructing algorithm and actual fortune of compressed sensing
With almost nobody is analyzed such as frequency domain characteristic of its imaging results；Further, since ghost imaging is all often using saturating
Penetrate the structure of formula, conventional active illumination imaging is then reflective structure, in existing technology, also without both respectively and
The system that conventional imaging mode is integrated.And in research before this, be often limited to a kind of imaging method replace it is another,
But those skilled in the art do not notice both of the above for the information of image have it is different stress, the present invention in initiate ground
With reference to the frequency domain characteristic of two methods, so as to obtain while better than the result of two methods.And instant invention overcomes prior art
The prejudice of middle research, two kinds of imaging modes is integrated into same set of imaging device so that the imaging device can realize two kinds
The imaging mode of different principle, and the image for obtaining respectively obtaining better than both by above-mentioned specific algorithm, are substantially increased
Image quality under the conditions of powerful gas scattering.
Above content is to combine specific preferred embodiment further description made for the present invention, it is impossible to assert
The specific implementation of the present invention is confined to these explanations.For those skilled in the art, do not taking off
On the premise of from present inventive concept, some equivalent substitutes or obvious modification can also be made, and performance or purposes are identical, all should
When being considered as belonging to protection scope of the present invention.
Claims (10)
1. the imaging device under the conditions of a kind of powerful gas scattering, it is characterised in that the imaging device is used to pass through scattering medium
Object is imaged, including laser, spatial light modulator, the first lens, the second lens and imaging sensor, wherein,
The spatial light modulator includes multiple turnover micro mirrors, and the spatial light modulator is arranged on swashing for the laser injection
In light light path, the laser optical path is transmitted by first lens again after the micro mirror reflection through the spatial light modulator
Penetrate afterwards and through the scattering medium on the object, then the laser optical path reflects through the object and passes through institute
State and penetrated again after second lens transmission in described image sensor to be carried out into the object after scattering medium
Picture.
2. imaging device according to claim 1, it is characterised in that the laser uses wavelength for 720~904nm's
LASER Light Source.
3. imaging device according to claim 1 or 2, it is characterised in that the spatial light modulator includes M × N number of square
The turnover micro mirror of battle array arrangement.
4. the imaging method under the conditions of a kind of powerful gas scattering, it is characterised in that entered using the imaging device described in claim 3
Row imaging, comprises the following steps：
S1：The calculation matrix for the M × N for being all 1 by one is input to the spatial light modulator, raw in described image sensor
Into the first image, wherein in the calculation matrix 1 represent by the corresponding micro mirror in the spatial light modulator overturn to
The laser optical path that the laser is projected is reflexed on the object；
S2：One group of calculation matrix comprising 0 and 1 M × N is input to the spatial light modulator, passes through the group measurement
The intensity signal that matrix and corresponding described image sensor are received, reduction the second image of generation, wherein in the calculation matrix
0 represent the corresponding micro mirror in the spatial light modulator is overturn it is anti-to the laser optical path not projected the laser
It is mapped on the object；
S3：Described first image is weighted with second image by the way of frequency domain weighting and is added, is generated finally
The synthetic image of the object.
5. imaging method according to claim 4, it is characterised in that step S1 also includes, is carried out to described first image
Filtering process, it is filtered described first image to generate the described first image in filtered first image, step S3.
6. imaging method according to claim 4, it is characterised in that pass through the group calculation matrix and phase in step S2
The intensity signal that the described image sensor answered is received, reduction the second image of generation is specifically included：Using following calculation formula：
Y=Φ x
Wherein, x is the image raw information of one-dimensional, and y is that the reflected light that m sampling described image sensor is received is always strong
Degree, Φ is calculation matrix collection, and m is the matrix quantity of calculation matrix described in one group, n=M × N；Can basis by above-mentioned formula
Φ and y reconstruct generation x, i.e. reduction generation second image.
7. imaging method according to claim 6, it is characterised in that the algorithm wherein reconstructed uses OMP algorithms.
8. the imaging method according to any one of claim 4 to 7, it is characterised in that step S3 is specifically included：
S31：Described first image and second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as described first image and described second
The weighting function of image；
S33：By the described first two-dimentional piecewise function and the second two-dimentional piecewise function respectively with described first image and described
The frequency-domain function that second image is obtained by Fourier transformation is multiplied, and is then added, and obtains comprehensive frequency-domain function, then carry out anti-
Fourier transformation, that is, generate the synthetic image of the final object；
Wherein in step s 32：
Described first two-dimentional piecewise function w1Relation with frequency f is as follows：
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>f</mi>
<mo><</mo>
<mn>0.1</mn>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>f</mi>
<mo>-</mo>
<mn>0.1</mn>
<mi>F</mi>
</mrow>
<mrow>
<mn>0.067</mn>
<mi>F</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0.1</mn>
<mi>F</mi>
<mo>&le;</mo>
<mi>f</mi>
<mo><</mo>
<mn>0.167</mn>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0.167</mn>
<mi>F</mi>
<mo>&le;</mo>
<mi>f</mi>
<mo>&le;</mo>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Described second two-dimentional piecewise function w2Relation with frequency f is as follows：
<mrow>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>f</mi>
<mo><</mo>
<mn>0.1</mn>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mn>0.167</mn>
<mi>F</mi>
<mo>-</mo>
<mi>f</mi>
</mrow>
<mrow>
<mn>0.067</mn>
<mi>F</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0.1</mn>
<mi>F</mi>
<mo>&le;</mo>
<mi>f</mi>
<mo><</mo>
<mn>0.167</mn>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0.167</mn>
<mi>F</mi>
<mo>&le;</mo>
<mi>f</mi>
<mo>&le;</mo>
<mi>F</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Wherein, F is highest frequency.
9. the imaging method according to any one of claim 4 to 7, it is characterised in that step S3 is specifically included：
S31：Described first image and second image are obtained into frequency-domain function respectively through Fourier transformation；
S32：Using the first two-dimensional Gaussian function and the second two-dimensional Gaussian function respectively as described first image and described second
The weighting function of image, wherein first two-dimensional Gaussian function and second two-dimensional Gaussian function are respectively through normalization,
And first two-dimensional Gaussian function and the second two-dimensional Gaussian function sum are 1；
S33：By first two-dimensional Gaussian function and second two-dimensional Gaussian function respectively with described first image and described
The frequency-domain function that second image is obtained by Fourier transformation is multiplied, and is then added, and obtains comprehensive frequency-domain function, then carry out anti-
Fourier transformation, that is, generate the synthetic image of the final object.
10. imaging method according to claim 9, it is characterised in that in step S32：
When frequency is less than first predetermined value in first two-dimensional Gaussian function, corresponding weight is 0, and frequency is more than second and made a reservation for
During value, corresponding weight is 1, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, corresponding weight
It is bigger；
When frequency is less than first predetermined value in second two-dimensional Gaussian function, corresponding weight is 1, and frequency is more than second and made a reservation for
During value, corresponding weight is 0, and when frequency is between first predetermined value and the second predetermined value, frequency is bigger, corresponding weight
It is smaller.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710517013.XA CN107315176B (en) | 2017-06-29 | 2017-06-29 | Imaging device and method under the conditions of a kind of powerful gas scattering |
| PCT/CN2017/102512 WO2019000659A1 (en) | 2017-06-29 | 2017-09-20 | Imaging device and method under strong atmospheric scattering condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710517013.XA CN107315176B (en) | 2017-06-29 | 2017-06-29 | Imaging device and method under the conditions of a kind of powerful gas scattering |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107315176A true CN107315176A (en) | 2017-11-03 |
| CN107315176B CN107315176B (en) | 2019-04-26 |
Family
ID=60181279
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710517013.XA Active CN107315176B (en) | 2017-06-29 | 2017-06-29 | Imaging device and method under the conditions of a kind of powerful gas scattering |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN107315176B (en) |
| WO (1) | WO2019000659A1 (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108594429A (en) * | 2018-04-13 | 2018-09-28 | 中国科学院光电研究院 | Saturating cloud and mist imaging method based on wavefront correction |
| CN109520969A (en) * | 2018-10-26 | 2019-03-26 | 中国科学院国家空间科学中心 | A kind of distributed diffusion imaging method based on atmospheric medium automodulation |
| CN110132901A (en) * | 2019-05-21 | 2019-08-16 | 北京理工大学 | The system and method that synthetic aperture wears scattering medium imaging |
| CN110865391A (en) * | 2019-11-14 | 2020-03-06 | 清华大学 | Millimeter wave terahertz multi-polarization imaging method and imaging system for target enhancement |
| CN111352126A (en) * | 2020-03-11 | 2020-06-30 | 中国科学院国家空间科学中心 | Single-pixel imaging method based on atmospheric scattering medium modulation |
| CN113992840A (en) * | 2021-09-15 | 2022-01-28 | 中国航天科工集团第二研究院 | Large-view-field high-resolution imaging method and device based on compressed sensing |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111652925B (en) * | 2020-06-29 | 2023-04-07 | 合肥中科迪宏自动化有限公司 | Method for extracting target global feature Hu invariant moment by using single-pixel imaging |
| CN115220061B (en) * | 2022-07-15 | 2024-05-10 | 哈工大机器人（合肥）国际创新研究院 | Orthogonal normalization-based deep learning polarization ghost imaging method and system |
| CN116609794B (en) * | 2023-07-21 | 2023-09-26 | 中国人民解放军国防科技大学 | Single-pixel imaging method, device and equipment based on radial Chebyshev light field |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1474162A (en) * | 2003-08-07 | 2004-02-11 | 中国科学技术大学 | Optical imaging method and device for invisible image |
| EP1480441B1 (en) * | 2003-05-21 | 2008-11-12 | Esko-Graphics Imaging GmbH | Method and apparatus for multi-track imaging using single-mode beams and diffraction-limited optics |
| US20130100525A1 (en) * | 2011-10-19 | 2013-04-25 | Su Yu CHIANG | Optical imaging system using structured illumination |
| CN103363924A (en) * | 2013-07-15 | 2013-10-23 | 中国科学院空间科学与应用研究中心 | Compressing three-dimension calculation ghost imaging system and method |
| CN105223582A (en) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | A kind of laser infrared radar imaging device based on compressed sensing and formation method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102192787A (en) * | 2010-03-04 | 2011-09-21 | 范冰清 | Infrared imaging detection system |
| CN102486410A (en) * | 2010-12-06 | 2012-06-06 | 中国科学院微电子研究所 | Optical imaging device |
-
2017
- 2017-06-29 CN CN201710517013.XA patent/CN107315176B/en active Active
- 2017-09-20 WO PCT/CN2017/102512 patent/WO2019000659A1/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1480441B1 (en) * | 2003-05-21 | 2008-11-12 | Esko-Graphics Imaging GmbH | Method and apparatus for multi-track imaging using single-mode beams and diffraction-limited optics |
| CN1474162A (en) * | 2003-08-07 | 2004-02-11 | 中国科学技术大学 | Optical imaging method and device for invisible image |
| US20130100525A1 (en) * | 2011-10-19 | 2013-04-25 | Su Yu CHIANG | Optical imaging system using structured illumination |
| CN103363924A (en) * | 2013-07-15 | 2013-10-23 | 中国科学院空间科学与应用研究中心 | Compressing three-dimension calculation ghost imaging system and method |
| CN105223582A (en) * | 2015-09-01 | 2016-01-06 | 西安交通大学 | A kind of laser infrared radar imaging device based on compressed sensing and formation method |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108594429A (en) * | 2018-04-13 | 2018-09-28 | 中国科学院光电研究院 | Saturating cloud and mist imaging method based on wavefront correction |
| CN109520969A (en) * | 2018-10-26 | 2019-03-26 | 中国科学院国家空间科学中心 | A kind of distributed diffusion imaging method based on atmospheric medium automodulation |
| CN109520969B (en) * | 2018-10-26 | 2021-03-09 | 中国科学院国家空间科学中心 | Distributed scattering imaging method based on atmospheric medium self-modulation |
| CN110132901A (en) * | 2019-05-21 | 2019-08-16 | 北京理工大学 | The system and method that synthetic aperture wears scattering medium imaging |
| CN110132901B (en) * | 2019-05-21 | 2020-07-31 | 北京理工大学 | System and method for synthetic aperture through scattering media imaging |
| CN110865391A (en) * | 2019-11-14 | 2020-03-06 | 清华大学 | Millimeter wave terahertz multi-polarization imaging method and imaging system for target enhancement |
| CN111352126A (en) * | 2020-03-11 | 2020-06-30 | 中国科学院国家空间科学中心 | Single-pixel imaging method based on atmospheric scattering medium modulation |
| CN111352126B (en) * | 2020-03-11 | 2022-03-08 | 中国科学院国家空间科学中心 | Single-pixel imaging method based on atmospheric scattering medium modulation |
| CN113992840A (en) * | 2021-09-15 | 2022-01-28 | 中国航天科工集团第二研究院 | Large-view-field high-resolution imaging method and device based on compressed sensing |
| CN113992840B (en) * | 2021-09-15 | 2023-06-23 | 中国航天科工集团第二研究院 | Large-view-field high-resolution imaging method and device based on compressed sensing |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107315176B (en) | 2019-04-26 |
| WO2019000659A1 (en) | 2019-01-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107315176A (en) | Imaging device and method under the conditions of a kind of powerful gas scattering | |
| Hunt et al. | Karhunen-Loeve multispectral image restoration, part I: Theory | |
| Lu et al. | Underwater image descattering and quality assessment | |
| CN105139367A (en) | Visible-light polarization image fusion method based on non-subsampled shearlets | |
| US10237519B2 (en) | Imaging apparatus, imaging system, image generation apparatus, and color filter | |
| Odegard et al. | Wavelet-based SAR speckle reduction and image compression | |
| CN106033599A (en) | Visible light enhancement method based on polarized imaging | |
| CN102044070A (en) | Retinex based nonlinear color image enhancement method | |
| CN101950412A (en) | Method for enhancing details and compressing dynamic range of infrared image | |
| CN101359399A (en) | Cloud-removing method for optical image | |
| CN106504208A (en) | Based on orderly minima and the high-spectrum image width destriping method of wavelet filtering | |
| US10446600B2 (en) | Imaging system and imaging device having a random optical filter array | |
| Petrovic et al. | Cross-band pixel selection in multiresolution image fusion | |
| CN106991670A (en) | One kind is without reference noise image quality evaluating method and system | |
| CN116739958B (en) | Dual-spectrum polarization super-resolution fusion detection method and system | |
| Roy et al. | A survey on visibility enhancement techniques in degraded atmospheric outdoor scenes | |
| Jain et al. | Edge-based prediction for lossless compression of hyperspectral images | |
| Wu et al. | Remote sensing image data fusion based on ihs and local deviation of wavelet transformation | |
| Widjaja | Noisy face recognition using compression-based joint wavelet-transform correlator | |
| CN112950507B (en) | Method for improving single-pixel color imaging performance under scattering environment based on deep learning | |
| Usha et al. | Atmospheric correction of remotely sensed multispectral satellite images in transform domain | |
| Hu et al. | Stripe Noise Removal for Infrared Image by Regularized Spectral Separation | |
| Bansal et al. | Image Enhancement In HSI Space Using Wavelet Transform | |
| Talal et al. | Fusion-based Resolution Enhancement of Satellite Images: Comparative Study and Performance Evaluation | |
| Regan et al. | Enhancing Dehazing Performance of Single Optical Satellite Images using Gamma Correction and Improved DCP |
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 |