CN107315176B - 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
- CN107315176B CN107315176B CN201710517013.XA CN201710517013A CN107315176B CN 107315176 B CN107315176 B CN 107315176B CN 201710517013 A CN201710517013 A CN 201710517013A CN 107315176 B CN107315176 B CN 107315176B
- Authority
- CN
- China
- Prior art keywords
- image
- frequency
- function
- spatial light
- scattering
- 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.)
- Active
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 devices under the conditions of a kind of powerful gas scattering, the imaging device is used to that object to be imaged 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 projects, the laser optical path is penetrated on the object after the reflection of the micro mirror of the spatial light modulator using after first lens transmission and through the scattering medium, then the laser optical path is penetrated in described image sensor after transmiting after the object reflects and penetrates the scattering medium using second lens the object to be imaged.The invention also discloses the imaging methods 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 scattering of powerful gas.
Description
Technical field
Imaging device and side under the conditions of being scattered the present invention relates to picture imaging techniques field more particularly to a kind of powerful gas
Method.
Background technique
The industries such as aviation, navigation and highway communication have widely the imaging in the strong scattering mediums such as haze, misty rain
Demand.Imaging method in existing common strong scattering medium is the imaging of 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, under the conditions of powerful gas scatters, image quality substantially reduces；
The image quality under the conditions of 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 auxiliary and understands design and technical solution of the invention, not necessarily
The prior art for belonging to present patent application, no tangible proof show above content present patent application the applying date
In disclosed situation, above-mentioned background technique should not be taken to the novelty and creativeness of evaluation the application.
Summary of the invention
In order to improve the image quality under the conditions of scattering of powerful gas, under the conditions of the present invention proposes that a kind of powerful gas scatters
Imaging device and method.
In order to achieve the above object, the invention adopts the 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,
In, the spatial light modulator includes multiple turnover micro mirrors, and the spatial light modulator setting is projected in the laser
Laser optical path on, the laser optical path through the spatial light modulator the micro mirror reflection after using first lens
It is penetrated on the object after transmission and through the scattering medium, then the laser optical path reflects and saturating through the object
Cross after the scattering medium using penetrated after second lens transmission in described image sensor with to the object into
Row imaging.
Preferably, the laser uses wavelength for the laser light source of 720~904nm.
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, using above-mentioned imaging device carry out at
Picture, comprising the following steps:
S1: the calculation matrix for M × N that one is all 1 is input to the spatial light modulator, in described image sensor
The first image of upper generation, wherein 1 in the calculation matrix indicates to turn over the micro mirror corresponding in the spatial light modulator
It goes to and the laser optical path that the laser projects is reflected 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
Calculation matrix and the received intensity signal of corresponding described image sensor, reduction generates the second image, wherein the measurement square
0 in battle array indicates to overturn the micro mirror corresponding in the spatial light modulator to the laser light not projected the laser
Road is reflected on the object；
The first image: being weighted with second image by the way of frequency domain weighting and be added by S3, is generated most
The synthetic image of the whole object.
Preferably, step S1 further includes being filtered to the first image, generates filtered first image,
The first image in step S3 is filtered the first image.
Preferably, pass through the group calculation matrix and the received light intensity letter of corresponding described image sensor in step S2
Breath, reduction generate the second image and specifically include: use 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 receives is total
Intensity, Φ are calculation matrix collection, and m is the matrix quantity of calculation matrix described in one group, n=M × N；It can root by above-mentioned formula
X is generated according to Φ and y reconstruct, i.e. reduction generates second image.
Preferably, wherein the algorithm reconstructed uses OMP algorithm.
Preferably, step S3 is specifically included:
S31: the first image and second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
S32: using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as the 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 the first image and
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 into
Row inversefouriertransform generates the synthetic image of the final object；
Wherein in step s 32:
Described first two-dimentional piecewise function w1It is as follows with the relationship of frequency f:
Described second two-dimentional piecewise function w2It is as follows with the relationship of frequency f:
Wherein, F is highest frequency.
Preferably, step S3 is specifically included:
S31: the first image and second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
S32: using the first two-dimensional Gaussian function and the second two-dimensional Gaussian function respectively as the first image and described
The weighting function of second image, wherein first two-dimensional Gaussian function and second two-dimensional Gaussian function pass through normalizing respectively
Change, and the sum of first two-dimensional Gaussian function and second two-dimensional Gaussian function are 1；
S33: by first two-dimensional Gaussian function and second two-dimensional Gaussian function respectively with the first image and
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 into
Row inversefouriertransform generates 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 greater than second
When 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 greater than second
When 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 at
As device can be realized simultaneously the imaging mode of two kinds of different principles, including active illumination imaging method and the imaging of compressed sensing ghost
Method is allowed and obtains the first image obtained using active illumination imaging method and use simultaneously by the imaging device
The second image that compressed sensing ghost imaging method obtains, so as to which the first image and the second image are further carried out General Office
Reason, to obtain preferably synthetic image, to substantially increase the image quality under the conditions of scattering of powerful gas.
In further embodiment, laser uses wavelength for the laser light source of 720nm~904, so that laser issues
Laser optical path to the scattering medium in atmosphere have better penetrability, and keep diffraction effect will not occur.Spatial light tune
Device processed includes the turnover micro mirror of M × N number of matrix arrangement, so as to which the calculation matrix of M × N is input to space light modulation
Device plays the role of modulated light source, by the way that one group of calculation matrix being randomly generated is carried out corresponding the second image of Self -adaptive,
Reduce the sampling number for generating the second image.
In further scheme, the imaging in conjunction with active illumination imaging method and compressed sensing ghost imaging method is special
Property, it can weight letter by Gaussian function or specific piecewise function as the first image and the second image in the present invention
Number, the method for line frequency domain weighting summation of going forward side by side obtain final synthetic image and are superior to the first image and the second image.
Detailed description of the invention
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 is the spectrogram of the original image 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 is the schematic diagram of the first two-dimensional Gaussian function of the embodiment of the present invention one；
Fig. 5 b is the schematic diagram of the second two-dimensional Gaussian function of the embodiment of the present invention one；
Fig. 6 a is the result that the first two-dimensional Gaussian function of the embodiment of the present invention one is multiplied with the frequency-domain function of the first image
Schematic diagram；
Fig. 6 b is that the result that the second two-dimensional Gaussian function of inventive embodiments one is multiplied with the frequency-domain function of the second image is shown
It is intended to；
Fig. 6 c is the result of Fig. 6 a with Fig. 6 b being added；
Fig. 7 a is the schematic diagram for the first image that the embodiment of the present invention one obtains；
Fig. 7 b is the schematic diagram for the second image that the embodiment of the present invention one obtains；
Fig. 7 c is the schematic diagram for the synthetic image that the embodiment of the present invention one obtains；
Fig. 8 a is the schematic diagram of the original image image of two object of the embodiment of the present invention；
Fig. 8 b is the schematic diagram for the first image that the embodiment of the present invention two obtains；
Fig. 8 c is the schematic diagram for filtered first image that Fig. 8 b is obtained by gaussian filtering；
Fig. 8 d is the schematic diagram for the second image that the embodiment of the present invention two obtains；
Fig. 8 e is the schematic diagram for the synthetic image that the embodiment of the present invention two obtains；
Fig. 9 a is the schematic diagram of the original image image of three object of the embodiment of the present invention；
Fig. 9 b is the schematic diagram for the first image that the embodiment of the present invention three obtains；
Fig. 9 c is the schematic diagram for filtered first image that Fig. 9 b is obtained by gaussian filtering；
Fig. 9 d is the schematic diagram for the second image that the embodiment of the present invention three obtains；
Fig. 9 e is the schematic diagram for the synthetic image that the embodiment of the present invention three obtains.
Specific embodiment
Below against attached drawing and in conjunction with 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, by the imaging device to object 60 carry out
Imaging, wherein there are scattering mediums 70 between the imaging device and object 60.The wherein primary structure of the imaging device are as follows:
Spatial light modulator 20 includes multiple turnover micro mirrors, and the laser optical path of the injection of laser 10 is arranged in spatial light modulator 20
On, laser optical path is after the reflection of the micro mirror of spatial light modulator 20 using after the transmission of the first lens 30 and through scattering medium 70
It penetrates on object 60, then laser optical path is reflected using object 60 and after transmission scattering medium 70 by the second lens 40
It is penetrated after transmission on imaging sensor 50 so that object 60 to be imaged.Wherein, laser 10 uses in some embodiments
Wavelength is the laser light source of 720~904nm, 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 the near-infrared laser of 808nm
Light source has relatively good penetrability to the misty rain in atmosphere, and spatial light modulator 20 includes that M × N number of matrix arrangement turns over
The micro mirror turned.Object is imaged by the imaging device, comprising the following steps:
S1: being input to spatial light modulator 20 for the calculation matrix for M × N that one is all 1, raw on imaging sensor 50
At the first image, wherein 1 in calculation matrix indicates to overturn micro mirror corresponding in spatial light modulator 20 to by laser 10
The laser optical path of injection is reflected on object 60；
At this point, spatial light modulator 20 reflects all light, entire optical 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 some embodiments, also the first image is filtered, filtered first image is generated, wherein filtering
Processing can use gaussian filtering method.
S2: one group of calculation matrix (can be and be randomly generated) comprising 0 and 1 M × N is input to spatial light modulator
20, by this group of calculation matrix and the received intensity signal of corresponding imaging sensor 50, reduction generates the second image, wherein surveying
0 expression in moment matrix overturns micro mirror corresponding 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 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, realized to light
The modulation in source；It determines whether the light in the region reflects to object 60 whether the overturning of micro mirror array, and then determines target
Whether the corresponding region on 60 surface of object is illuminated namely the practical area for having corresponded to body surface and being illuminated of each calculation matrix
Domain；By multiple calculation matrix in one group and the received intensity signal of corresponding imaging sensor, reduction generates the second image, 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 receives is always strong
Degree, Φ are calculation matrix collection, and m is the matrix quantity of one group of calculation matrix, every a line that n=M × N namely calculation matrix are concentrated
One group of coding (a corresponding calculation matrix) of i.e. corresponding primary sampling spatial light modulator；It can basis by above-mentioned formula
Φ and y reconstruct generates x, i.e. reduction generates the second image；The algorithm wherein reconstructed can use OMP algorithm (orthogonal matching pursuit
Algorithm).
S3: being weighted addition for the first image and the second image by the way of frequency domain weighting, generates final target
The synthetic image of object.
Pass through the spectrogram to active illumination imaging method under the conditions of different scattering coefficients and compressed sensing ghost imaging method
It makes comparisons respectively with the spectrogram of the original image of object, (abscissa is frequency to the spectrogram of original image, 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
The spectrogram of image is as shown in Figure 2 b, and the spectrogram for the second image that compressed sensing ghost imaging method obtains is as shown in Figure 2 c, leads to
Cross compare it can be seen that entire wave band be nearly all active illumination imaging method advantage section, highest frequency 1/10th with
Upper this section (5-45) may be considered absolute section；Under high scattering coefficient (6.5), active illumination imaging method is obtained
The spectrogram of first image is as shown in Figure 2 d, the spectrogram for the second image that compressed sensing ghost imaging method obtains such as Fig. 2 e institute
Show, by comparing it can be seen that being two on the basis of high band still should be by active illumination imaging results, but in highest frequency
/ 10th or less (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 compressed sensing ghost imaging method, in some embodiments,
The first two-dimentional piecewise function and the second two-dimentional piecewise function can be used respectively as the weight letter of the first image and the second image
Number directly uses the frequency spectrum of advantage method in the absolute predominance section of the two to carry out frequency domain weighting, smoothed in intermediate region
It crosses, specific as follows:
First two-dimentional piecewise function w1It is as follows with the relationship of frequency f:
Second two-dimentional piecewise function w2It is as follows with the relationship of frequency f, as shown in Figure 3:
Wherein, F is highest frequency, which 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: the first image and the second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
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: the 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 obtains is multiplied, and is then added, the frequency-domain function of obtained synthesis, then carry out anti-Fourier's change
It changes, that is, generates the synthetic image of final object.
It can be seen that the function similar to Gaussian function, therefore, in this hair from the schematic diagram of piecewise function shown in Fig. 3
In other bright embodiments, Gaussian function also can be used as weighting function to carry out frequency domain weighting, as shown in figures 4 a and 4b,
It is compared using piecewise function and Gaussian function as weighting function processing, Fig. 4 a is to pass through piecewise function integrated treatment
Synthetic image, SSIM 0.76053, Fig. 4 b are the synthetic image by Gaussian function integrated treatment, and SSIM 0.78426 can
To find out that the comprehensive effect of piecewise function integrated treatment and Gaussian function almost without difference, namely passes 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 carrys out value, i.e. step S3 can also specifically be wrapped
It includes:
S31: the first image and the second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
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 be respectively by normalization, and the first two dimension is high
The sum of this function and the second two-dimensional Gaussian function are 1；
S33: the 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 obtains is multiplied, and is then added, the frequency-domain function of obtained synthesis, then carry out anti-Fourier's change
It changes, that is, generates 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 greater than the
When two predetermined values, 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 it is predetermined that frequency is greater than second
When 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 is as shown in Figure 5 a, the schematic diagram of the second two-dimensional Gaussian function such as Fig. 5 b institute
Show, that the two is added and be 1, the two passes through the frequency-domain function phase that Fourier transformation obtains with the first image and the second image respectively
Multiply (frequency-domain function multiplied result such as Fig. 6 a institute that the first two-dimensional Gaussian function and the first image are obtained by passing through Fourier transformation
Show, the second two-dimensional Gaussian function and the second image pass through the frequency-domain function multiplied result such as Fig. 6 b institute obtained by Fourier transformation
Show), it is then added, obtains comprehensive frequency-domain function as fig. 6 c, using Fourier transformation is returned, that is, generate final target
The synthetic image of object is as shown in Figure 7 c, wherein by step S1 and S2 the first image respectively obtained and the second image respectively as schemed
Shown in 7a and Fig. 7 b, Fig. 7 c is 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 original image image of object is as shown in Figure 8 a, under conditions of the forward scattering coefficient bd of scattering medium is 4.5, with
Evaluation criterion of the structural similarity (SSIM) as picture quality.The first image obtained according to step S1 is as shown in Figure 8 b,
SSIM is 0.7564, and the first image after gaussian filtering is as shown in Figure 8 c, SSIM 0.89081；It is obtained according to step S2
The second image as shown in figure 8d, SSIM 0.60799；According to step S3, using the first dimensional Gaussian of such as Fig. 5 a and Fig. 5 b
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 0.88449.
Embodiment three:
The original image image of object as illustrated in fig. 9, the forward scattering coefficient bd of scattering medium be 5.5 under conditions of, with
Evaluation criterion of the structural similarity (SSIM) as picture quality.The first image obtained according to step S1 as shown in figure 9b,
SSIM is 0.20654, the first image after gaussian filtering as is shown in fig. 9 c, SSIM 0.43112；It is obtained according to step S2
The second image as shown in figure 9d, SSIM 0.37954；According to step S3, using the first dimensional Gaussian of such as Fig. 5 a and Fig. 5 b
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 0.54104.
Embodiment two and embodiment three are the synthesis quality reconstruction under different scattering coefficients respectively, it can be seen that are passed through
It, can be very close to preferable imaging effect, both bad in two methods gap great disparity in conjunction with two kinds of imaging methods
When, available effect is better than the image of the two simultaneously.
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 obtains, wherein active illumination imaging method and compressed sensing ghost imaging method are in principle
There is very big difference, while also variant on imaging characteristic.By applicant's the study found that compressed sensing ghost imaging method pair
Random noise is insensitive, but more sensitive to whole 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 range
It is interior to retain medium-high frequency information well.In order to inhibit noise, also active illumination is imaged in the preferred embodiment of the present invention
The image that method obtains is filtered, such as gaussian filtering.
In current picture imaging techniques field, research mainly concentrates on the restructing algorithm and reality fortune of compressed sensing
With being analyzed almost without such as frequency domain characteristic of the people to its imaging results；In addition, 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, often it is limited to replace another kind with a kind of imaging method,
However 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
In conjunction with the frequency domain characteristic of two methods, to obtain the result for being better than two methods simultaneously.And the present invention overcomes the prior arts
Two kinds of imaging modes are integrated into same set of imaging device by the prejudice of middle research, so that the imaging device may be implemented two kinds
The imaging mode of different principle, and the image respectively obtained better than the two is obtained by above-mentioned specific algorithm, it substantially increases
Image quality under the conditions of powerful gas scattering.
The above content is a further detailed description of the present invention in conjunction with specific preferred embodiments, and it cannot be said that
Specific implementation of the invention is only limited to these instructions.For those skilled in the art to which the present invention belongs, it is not taking off
Under the premise of from present inventive concept, several equivalent substitute or obvious modifications can also be made, and performance or use is identical, all answered
When being considered as belonging to protection scope of the present invention.
Claims (8)
1. the imaging method under the conditions of a kind of powerful gas scattering, which is characterized in that using the imaging dress under the conditions of the scattering of powerful gas
It sets and is imaged, wherein the imaging device under the conditions of the scattering of powerful gas is used to that object to be imaged through scattering medium, wrapped
Include laser, spatial light modulator, the first lens, the second lens and imaging sensor, wherein the spatial light modulator includes
Multiple turnover micro mirrors, the spatial light modulator are arranged on the laser optical path that the laser projects, the laser light
Road is situated between after the reflection of the micro mirror of the spatial light modulator using after first lens transmission and through the scattering
Matter is penetrated on the object, then the laser optical path after the object reflects and penetrates the scattering medium using
It penetrates after the second lens transmission the object to be imaged in described image sensor, the spatial light modulator
Including the turnover micro mirror of M × N number of matrix arrangement；
Imaging method under the conditions of powerful gas scattering the following steps are included:
S1: being input to the spatial light modulator for the calculation matrix for M × N that one is all 1, raw in described image sensor
At the first image, wherein in the calculation matrix 1 indicate by the micro mirror corresponding in the spatial light modulator overturn to
The laser optical path that the laser projects is reflected 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
Matrix and the received intensity signal of corresponding described image sensor, reduction generates the second image, wherein in the calculation matrix
0 indicate the micro mirror corresponding 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: the first image is weighted with second image by the way of frequency domain weighting and is added, is generated final
The synthetic image of the object.
2. imaging method according to claim 1, which is characterized in that the laser uses wavelength for 720~904nm's
Laser light source.
3. imaging method according to claim 1, which is characterized in that step S1 further includes carrying out to the first image
Filtering processing generates filtered first image, and the first image in step S3 is filtered the first image.
4. imaging method according to claim 1, which is characterized in that pass through the group calculation matrix and phase in step S2
The received intensity signal of the described image sensor answered, reduction generate the second image and specifically include: use 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 receives is always strong
Degree, Φ is calculation matrix collection, and m is the matrix quantity of calculation matrix described in one group, n=M × N；It can basis by above-mentioned formula
Φ and y reconstruct generates x, i.e. reduction generates second image.
5. imaging method according to claim 4, which is characterized in that the algorithm wherein reconstructed uses OMP algorithm.
6. imaging method according to any one of claims 1 to 5, which is characterized in that step S3 is specifically included:
S31: the first image and second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
S32: using the first two-dimentional piecewise function and the second two-dimentional piecewise function respectively as the 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 the first image and described
Second image is multiplied by the frequency-domain function that Fourier transformation obtains, and is then added, obtains comprehensive frequency-domain function, then carries out anti-
Fourier transformation generates the synthetic image of the final object；
Wherein in step s 32:
Described first two-dimentional piecewise function w1It is as follows with the relationship of frequency f:
Described second two-dimentional piecewise function w2It is as follows with the relationship of frequency f:
Wherein, F is highest frequency.
7. imaging method according to any one of claims 1 to 5, which is characterized in that step S3 is specifically included:
S31: the first image and second image are passed through into Fourier transformation respectively and obtain frequency-domain function；
S32: using the first two-dimensional Gaussian function and the second two-dimensional Gaussian function respectively as the first image and described second
The weighting function of image, wherein first two-dimensional Gaussian function and second two-dimensional Gaussian function pass through normalization respectively,
And the sum of first two-dimensional Gaussian function and second two-dimensional Gaussian function are 1；
S33: by first two-dimensional Gaussian function and second two-dimensional Gaussian function respectively with the first image and described
Second image is multiplied by the frequency-domain function that Fourier transformation obtains, and is then added, obtains comprehensive frequency-domain function, then carries out anti-
Fourier transformation generates the synthetic image of the final object.
8. imaging method according to claim 7, which is characterized in that in step S32:
When frequency is less than first predetermined value in first two-dimensional Gaussian function, corresponding weight is 0, and it is predetermined that frequency is greater than second
When 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 it is predetermined that frequency is greater than second
When 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 CN107315176A (en) | 2017-11-03 |
| CN107315176B true 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) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108594429B (en) * | 2018-04-13 | 2021-06-08 | 中国科学院光电研究院 | Cloud-penetrating imaging method based on wavefront correction |
| CN109520969B (en) * | 2018-10-26 | 2021-03-09 | 中国科学院国家空间科学中心 | Distributed scattering imaging method based on atmospheric medium self-modulation |
| CN110132901B (en) * | 2019-05-21 | 2020-07-31 | 北京理工大学 | System and method for synthetic aperture through scattering media imaging |
| CN110865391B (en) * | 2019-11-14 | 2021-09-21 | 清华大学 | Millimeter wave terahertz multi-polarization imaging method and imaging system for target enhancement |
| CN111352126B (en) * | 2020-03-11 | 2022-03-08 | 中国科学院国家空间科学中心 | Single-pixel imaging method based on atmospheric scattering medium modulation |
| CN111652925B (en) * | 2020-06-29 | 2023-04-07 | 合肥中科迪宏自动化有限公司 | Method for extracting target global feature Hu invariant moment by using single-pixel imaging |
| CN113992840B (en) * | 2021-09-15 | 2023-06-23 | 中国航天科工集团第二研究院 | Large-view-field high-resolution imaging method and device based on compressed sensing |
| 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 (4)
| 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 |
| 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 (3)
| 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 |
| US9599805B2 (en) * | 2011-10-19 | 2017-03-21 | National Synchrotron Radiation Research Center | Optical imaging system using structured illumination |
-
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 (4)
| 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 |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2019000659A1 (en) | 2019-01-03 |
| CN107315176A (en) | 2017-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107315176B (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 | |
| Odegard et al. | Wavelet-based SAR speckle reduction and image compression | |
| KR20150008877A (en) | Image processing device, image processing method, program, print medium, and recording medium | |
| CN109636766A (en) | Polarization differential and intensity image Multiscale Fusion method based on marginal information enhancing | |
| CN104978724A (en) | Infrared polarization fusion method based on multi-scale transformation and pulse coupled neural network | |
| CN101359399A (en) | Cloud-removing method for optical image | |
| Zhang et al. | Fusion of infrared and visible light images based on nonsubsampled shearlet transform | |
| Adu et al. | Image fusion based on visual salient features and the cross-contrast | |
| CN110349101B (en) | Multi-scale combined wavelet image denoising method | |
| Sharumathi et al. | A survey on various image enhancement techniques for underwater acoustic images | |
| CN106504208A (en) | Based on orderly minima and the high-spectrum image width destriping method of wavelet filtering | |
| Gao et al. | Infrared and visible image fusion using dual-tree complex wavelet transform and convolutional sparse representation | |
| Su-xia et al. | Image fusion based on regional energy and standard deviation | |
| Liu et al. | Infrared and visible image fusion for shipborne electro-optical pod in maritime environment | |
| Huang et al. | An NSST-Based Fusion Method for Airborne Dual-Frequency, High-Spatial-Resolution SAR Images | |
| Chen et al. | Modification of standard image compression methods for correlation-based pattern recognition | |
| Widjaja | Noisy face recognition using compression-based joint wavelet-transform correlator | |
| Benedetto et al. | Wavelet packets for multi-and hyper-spectral imagery | |
| El-Tawel et al. | Fusion of multispectral and full polarimetric SAR images in NSST domain | |
| Schmalz et al. | Data compression techniques for underwater imagery | |
| Bansal et al. | Image Enhancement In HSI Space Using Wavelet Transform | |
| CN111199530A (en) | Fusion method of SAR image and visible light image |
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 |