CN106290285B - A kind of non-intrusion type laser scanning imaging method based on stochastical sampling - Google Patents
A kind of non-intrusion type laser scanning imaging method based on stochastical sampling Download PDFInfo
- Publication number
- CN106290285B CN106290285B CN201610834865.7A CN201610834865A CN106290285B CN 106290285 B CN106290285 B CN 106290285B CN 201610834865 A CN201610834865 A CN 201610834865A CN 106290285 B CN106290285 B CN 106290285B
- Authority
- CN
- China
- Prior art keywords
- matrix
- fluorescence intensity
- laser scanning
- imaging method
- intrusion type
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/645—Specially adapted constructive features of fluorimeters
- G01N21/6456—Spatial resolved fluorescence measurements; Imaging
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The non-intrusion type laser scanning imaging method based on stochastical sampling that the invention discloses a kind of, include: S1: generating the random matrix Φ that element only includes 0 and 1, wherein the size of random matrix Φ is identical as the size of fluorescence intensity matrix I and corresponding scanning angle matrix Θ；S2: the angle on the corresponding scanning angle matrix Θ of element for being 1 by laser scanning random matrix Φ intermediate value collects the fluorescence intensity on the corresponding fluorescence intensity matrix I of element that random matrix Φ intermediate value is 1, obtains incomplete fluorescence intensity matrixWherein ο symbol indicates that corresponding element is multiplied；S3: by solving the reconstruction model under low-rank constraint, to find out complete fluorescence intensity matrix I=U；S4: the picture of object is recovered from the complete fluorescence intensity matrix reconstructed.Non-intrusion type laser scanning imaging method proposed by the present invention based on stochastical sampling and low-rank reconstruct, the time required to capable of substantially reducing scanning process, while keeping higher image quality.
Description
Technical field
The present invention relates to calculating imaging, biomedical imaging, the technical fields such as image reconstruction, more particularly to it is a kind of be based on
The non-intrusion type laser scanning imaging method of machine sampling.
Background technique
In fields such as medical imaging, industrial detections, generally require to small knots such as biological tissue cell, industrial chips
Structure is imaged using the foundation as analyzing and diagnosing and detection.But since these image forming mediums are often translucent scattering
Layer, traditional imaging method based on geometric optics are no longer applicable in, except non-demolition scattering layer or toward injection in scattering layer assist at
The substance of picture, but these means are easy to damage observed object.In recent years, a kind of based on the non-of laser speckle scanning
Intrusive imaging method is suggested, it can obtain the object for being hidden in scattering layer behind under the premise of not destroying scattering layer
Clearly as.In this approach, beam of laser is mapped to the fixation position of scattering layer, forms laser speckle simultaneously under scattering process
It is radiated in the plane where fluorescent object.The fluorescence that the speckle fallen on fluorescent object inspires be reflected back toward scattering layer and by
Acquisition, total fluorescence volume of the sum that collected fluorescent intensity adds up as laser under the incident angle.By laser according to sweeping
The angle retouched in angle matrix is scanned one by one, available corresponding fluorescence intensity matrix, i.e., every in fluorescence intensity matrix
One element value is laser according to total fluorescence volume under the angle incidence of corresponding element in scanning angle matrix.It finally can use phase
Bit recovery algorithm recovers the picture of object from fluorescence intensity matrix.In this traditional method, in order to obtain preferably at
As effect, the size of scan matrix needs sufficiently large, it is therefore desirable to total number of angles of laser scanning is huge, and due to
Fluorescence signal is weaker in practice, needs to increase the time for exposure of fluorescent collecting device, finally makes the acquisition of total imaging data
Time is very very long.Due to prolonged laser irradiation, it be easy to cause the damage of tested article, it is therefore necessary to it is multiple to reduce scanning
Miscellaneous degree reduces sweep time.
Fig. 1 show traditional non-intrusion type imaging device schematic diagram based on laser speckle scanning, including translucent scattered
Penetrate layer 1, fluorescent object 2, translucent scattering layer 3 and optical filter 4.When beam of laser vertically injects scattering layer, in fluorescent object
Place plane (being denoted as u-v plane) forms speckle pattern, is denoted as S (u, v).When laser with normal (θ=(θ into θ anglex,θy)) enter
When penetrating, " memory effect " existing for speckle translates the speckle in u-v plane at this time only, and pattern does not occur substantially
Change, i.e., S ' at this time=S (u-d1θx,v-d1θy).To in total fluorescence volume at fluorescent collecting end are as follows:
I (θ)=∫ ∫ O (u, v) S (u-d1θx,v-d1θy) dudv=[O*S] (θ).
Incident laser is scanned one by one according to the angle on scanning angle matrix Θ, is recorded under each incident angle
Total fluorescence volume, ultimately form fluorescence intensity matrix I, then restored from the fluorescence intensity matrix using Phase Retrieve Algorithm
To the image of the fluorescent object.
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 solve the above technical problems, the present invention proposes a kind of non-intrusion type laser based on stochastical sampling and low-rank reconstruct
Scan imaging method the time required to capable of substantially reducing scanning process, while keeping higher image quality.
In order to achieve the above objectives, the invention adopts the following technical scheme:
The non-intrusion type laser scanning imaging method based on stochastical sampling that the invention discloses a kind of, comprising the following steps:
S1: the random matrix Φ that element only includes 0 and 1 is generated, wherein the size of random matrix Φ and fluorescence intensity matrix I
And the size of corresponding scanning angle matrix Θ is identical；
S2: the angle on the corresponding scanning angle matrix Θ of element for being 1 by laser scanning random matrix Φ intermediate value is received
Collect the fluorescence intensity on the corresponding fluorescence intensity matrix I of element that random matrix Φ intermediate value is 1, obtains incomplete fluorescence intensity
MatrixWhereinSymbol indicates that corresponding element is multiplied；
S3: by solving the reconstruction model under following low-rank constraint, to find out complete fluorescence intensity matrix I=U；
Wherein, TV is total variance function, PkFor the matrix-block of p × p, w (Pk) it is matrix-block PkCorresponding weighting function, |
|·||*For nuclear norm, η is regular parameter, and Δ is all matrix-block P of matrix U to be askedkSet, Ω is the element acquired
Index set；
S4: the picture of object is recovered from the complete fluorescence intensity matrix reconstructed in step S3.
Preferably, 1 proportion of element in step S1 in random matrix Φ is 20% or more.
Preferably, laser scanning is carried out using galvanometer scanning system in step S2.
Preferably, the matrix-block P in step S3kCorresponding weighting function w (Pk) value is carried out according to the following formula:
Wherein, ξkFor matrix-block PkIn the number of element that acquires, λ is threshold value.
Preferably, wherein threshold value λ=0.2p2。
Preferably, it is specifically included in step S4: the complete fluorescence reconstructed from step S3 using Phase Retrieve Algorithm
The picture of object is recovered in intensity matrix.
Compared with prior art, the beneficial effects of the present invention are: the non-intrusion types of the invention based on stochastical sampling to swash
Optical scanning imaging method obtains incomplete fluorescence intensity by stochastical sampling using the method for stochastical sampling and low-rank reconstruct
Matrix then by rebuilding fluorescence intensity matrix with low-rank constraint and the restructing algorithm of total variance constraint, and eventually passes through phase
Bit recovery obtains the picture of object；By means of the present invention, the acquisition time of imaging data is greatly shortened, while can also be kept
Higher image quality.
Detailed description of the invention
Fig. 1 is traditional non-intrusion type imaging device schematic diagram based on laser speckle scanning；
Fig. 2 is that the process of the non-intrusion type laser scanning imaging method based on stochastical sampling of the preferred embodiment of the present invention is shown
It is intended to.
Fig. 3 is that the non-intrusion type laser scanning imaging method based on stochastical sampling of preferred embodiment through the invention carries out
The effect diagram of imaging.
Specific embodiment
Below against attached drawing and in conjunction with preferred embodiment, the invention will be further described.
As shown in Fig. 2, the preferred embodiment of the present invention discloses a kind of non-intrusion type laser scanning imaging based on stochastical sampling
Method, comprising the following steps:
S1: the random matrix Φ that element only includes 0 and 1 is generated, wherein the size of random matrix Φ and fluorescence intensity matrix I
And the size of corresponding scanning angle matrix Θ is identical；
Wherein, in random matrix Φ 1 proportion of element be ρ, then this method spent time be only conventional method ρ ×
100%；Ratio shared by the probability distribution of random matrix Φ, element 1 can be depending on actual imaging effect, more preferably
In embodiment, ρ >=0.2.
S2: the angle on the corresponding scanning angle matrix Θ of element for being 1 by laser scanning random matrix Φ intermediate value is received
Collect the fluorescence intensity on the corresponding fluorescence intensity matrix I of element that random matrix Φ intermediate value is 1, obtains incomplete fluorescence intensity
MatrixWhereinSymbol indicates that corresponding element is multiplied；
Wherein, complete fluorescence intensity matrix is used at random under the guidance of the step S1 random matrix Φ generated,
That is the element Φ that only acquisition random matrix Φ intermediate value is 1i,jCorresponding fluorescence intensity Ii,j；Laser scanning device is scanned using speckle
The device of imaging, wherein traditional galvanometer scanning system based on process control can be used, then control system only scans and refers to
Fixed angle and skip the angle for not needing acquisition；It can also be rotated using laser light source is directly controlled.
S3: by solving the reconstruction model under following low-rank constraint, to find out complete fluorescence intensity matrix I=U；
Wherein, TV (Total Variation) is total variance function, PkFor the matrix-block of p × p, w (Pk) it is matrix-block Pk
Corresponding weighting function, | | | |*For nuclear norm, η is regular parameter, and Δ is all matrix-block P of matrix U to be askedkSet,
Ω is the set of the index of the element acquired；Indicate optimized-typeU when being minimized, s.t. are the abbreviations of subject to, mean constraint condition.
Wherein, the constraint condition of the reconstruction model belongs to fidelity term, it is intended that the fluorescence intensity matrix of reconstruct connects as far as possible
Nearly imperfect fluorescence intensity matrix collected；In objective function, the nuclear norm of total variance function and the matrix-block with weighting
With for regular terms, it is intended that each local matrix-block of the matrix reconstructed is smooth and local low-rank, regular parameter η use
Regulate and control the specific gravity of two regular terms, is optimal the structure of final reconstruct.In addition, reconstruction model is because of different matrix-blocks
In known fluorescence intensity element proportion it is different and assign different matrix-blocks different weights, when specific implementation, can basis
Quality reconstruction optimum selecting weighting function.
When specific implementation, matrix-block PkCorresponding weighting function w (Pk) value is carried out according to the following formula:
Wherein, ξkFor matrix-block PkIn the number of element that acquires, λ is threshold value, can be set as 0.2p2, i.e. matrix
The selection of/5th of the total element number of block, actual weighting function and threshold value should be depending on practical quality reconstructions.
S4: the picture of object is recovered from the complete fluorescence intensity matrix reconstructed in step S3.
Specifically, object is recovered from the complete fluorescence intensity matrix reconstructed in step S3 using Phase Retrieve Algorithm
The picture of body may further use following steps:
1) auto-correlation I ★ I=(O*S) ★ (O*)=(O ★ O) * (S ★ S) for calculating fluorescence intensity matrix, due to speckle
Auto-correlation is peaking function, therefore has I ★ I ≈ O ★ O；Wherein, I is IntegratedIntensity (total light intensity), O Object
(object to be imaged), S are Speckle (laser speckle), and * is convolution symbol (Convolution), ★ cross-correlation symbol
(Correlation)。
(2) doing Fourier transformation to auto-correlation can obtain
It (3) can be from using Hybrid-Input-Output algorithm (Mixed design output algorithm)And O it is non-
O is recovered in the priori conditions such as negativity.
As shown in figure 3, be non-intrusion type laser scanning imaging method based on stochastical sampling through the invention carry out at
The effect diagram of picture, the leftmost side one is classified as the picture that traditional method all acquired obtains in Fig. 3, remaining is respectively classified as use
The obtained picture of method of the invention, wherein sample rate ρ from right to left is respectively 0.1,0.15,0.2,0.25, namely imaging
Acquisition time needed for data is only traditional 0.1 times, 0.15 times, 0.2 times and 0.25 times, that is, greatly shortens adopting for imaging data
Collect the time；And it can be seen from the figure that resulting image keeps higher image quality by means of the present invention, especially exist
ρ=0.2 and when ρ=0.25, has little effect image quality.I.e. in the application present invention, when ρ >=0.2, greatly shorten
The time required to scanning process, while not influencing image quality also.
Non-intrusion type laser scanning imaging method based on stochastical sampling of the invention is reconstructed using stochastical sampling and low-rank
Method, incomplete fluorescence intensity matrix is obtained by stochastical sampling, then by with low-rank constraint with total variance constrain
Restructing algorithm rebuild fluorescence intensity matrix, and eventually pass through phase recovery and obtain the picture of object；By means of the present invention, greatly
The big acquisition time for shortening imaging data, while higher image quality can also be kept.
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 (6)
1. a kind of non-intrusion type laser scanning imaging method based on stochastical sampling, which comprises the following steps:
S1: generating the random matrix Φ that element only includes 0 and 1, wherein the size of random matrix Φ and fluorescence intensity matrix I and
The size of corresponding scanning angle matrix Θ is identical；
S2: by laser scanning random matrix Φ intermediate value be 1 the corresponding scanning angle matrix Θ of element on angle, collect with
The fluorescence intensity on the corresponding fluorescence intensity matrix I of element that machine matrix Φ intermediate value is 1, obtains incomplete fluorescence intensity matrixWhereinSymbol indicates that corresponding element is multiplied；
S3: by solving the reconstruction model under following low-rank constraint, to find out complete fluorescence intensity matrix I=U；
Wherein, TV is total variance function, PkFor the matrix-block of p × p, w (Pk) it is matrix-block PkCorresponding weighting function, | | | |*
For nuclear norm, η is regular parameter, and Δ is all matrix-block P of matrix U to be askedkSet, Ω is the index of the element acquired
Set；
S4: the picture of object is recovered from the complete fluorescence intensity matrix reconstructed in step S3.
2. non-intrusion type laser scanning imaging method according to claim 1, which is characterized in that random matrix in step S1
1 proportion of element in Φ is 20% or more.
3. non-intrusion type laser scanning imaging method according to claim 1, which is characterized in that use galvanometer in step S2
Scanning system carries out laser scanning.
4. non-intrusion type laser scanning imaging method according to claim 1, which is characterized in that the matrix-block in step S3
PkCorresponding weighting function w (Pk) value is carried out according to the following formula:
Wherein, ξkFor matrix-block PkIn the number of element that acquires, λ is threshold value.
5. non-intrusion type laser scanning imaging method according to claim 4, which is characterized in that wherein threshold value λ=0.2
p2。
6. non-intrusion type laser scanning imaging method according to any one of claims 1 to 5, which is characterized in that step S4
In specifically include: object is recovered from the complete fluorescence intensity matrix reconstructed in step S3 using Phase Retrieve Algorithm
Picture:
(1) auto-correlation I ★ I=(O*S) ★ (O*S)=(O ★ O) * (S ★ S) of fluorescence intensity matrix is calculated, wherein I is total light
By force, O is object to be imaged, and S is laser speckle, and * is convolution symbol, ★ cross-correlation symbol；
(2) doing Fourier transformation to auto-correlation can obtain
(3) using Mixed design output algorithm fromAnd the priori conditions such as O nonnegativity in recover O.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610834865.7A CN106290285B (en) | 2016-09-20 | 2016-09-20 | A kind of non-intrusion type laser scanning imaging method based on stochastical sampling |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610834865.7A CN106290285B (en) | 2016-09-20 | 2016-09-20 | A kind of non-intrusion type laser scanning imaging method based on stochastical sampling |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106290285A CN106290285A (en) | 2017-01-04 |
| CN106290285B true CN106290285B (en) | 2019-04-16 |
Family
ID=57711414
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610834865.7A Active CN106290285B (en) | 2016-09-20 | 2016-09-20 | A kind of non-intrusion type laser scanning imaging method based on stochastical sampling |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106290285B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104865234A (en) * | 2015-06-03 | 2015-08-26 | 清华大学深圳研究生院 | Imaging method for noninvasive semitransparent imaging device |
| CN105182333A (en) * | 2015-08-24 | 2015-12-23 | 西安电子科技大学 | Sparse scene down-sampling SAR imaging method based on matrix filling |
| CN105259155A (en) * | 2015-11-16 | 2016-01-20 | 清华大学深圳研究生院 | Rapid non-invasive type semitransparent imaging method and device |
| CN105388135A (en) * | 2015-10-28 | 2016-03-09 | 清华大学深圳研究生院 | Non-invasive laser scanning imaging method |
-
2016
- 2016-09-20 CN CN201610834865.7A patent/CN106290285B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104865234A (en) * | 2015-06-03 | 2015-08-26 | 清华大学深圳研究生院 | Imaging method for noninvasive semitransparent imaging device |
| CN105182333A (en) * | 2015-08-24 | 2015-12-23 | 西安电子科技大学 | Sparse scene down-sampling SAR imaging method based on matrix filling |
| CN105388135A (en) * | 2015-10-28 | 2016-03-09 | 清华大学深圳研究生院 | Non-invasive laser scanning imaging method |
| CN105259155A (en) * | 2015-11-16 | 2016-01-20 | 清华大学深圳研究生院 | Rapid non-invasive type semitransparent imaging method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106290285A (en) | 2017-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Buchholz et al. | Cryo-CARE: Content-aware image restoration for cryo-transmission electron microscopy data | |
| US11368608B2 (en) | Compressed sensing based object imaging system and imaging method therefor | |
| WO2018213723A1 (en) | Imaging signal extraction apparatus and methods of using same | |
| CN108169173A (en) | A kind of big visual field high-resolution three dimensional diffraction chromatography micro imaging method | |
| CN101909525A (en) | Image processing method, device and laminagraph device | |
| CN107084789A (en) | Single pixel detector spectrum reflectivity reconstructing method based on sparse prior | |
| CN109752377B (en) | Spectroscopic bimodal projection tomography tissue blood vessel imaging device and method | |
| CN105455780A (en) | Fluorescent molecular tomography imaging reconstruction method based on limited projection of double grids | |
| US10830701B2 (en) | Scanned line angular projection microscopy | |
| CN105388135B (en) | A kind of non-intrusion type laser scanning imaging method | |
| JP2015522172A (en) | Method and system for restoring optical properties of a diffractive object immersed in a liquid medium | |
| Ma et al. | Cascade neural approximating for few-shot super-resolution photoacoustic angiography | |
| Dammers et al. | Automatic identification of gray and white matter components in polarized light imaging | |
| Wu et al. | Deep learning enhanced NIR-II volumetric imaging of whole mice vasculature | |
| Song et al. | Light-field microscopy for the optical imaging of neuronal activity: When model-based methods meet data-driven approaches | |
| CN106290285B (en) | A kind of non-intrusion type laser scanning imaging method based on stochastical sampling | |
| TWI554740B (en) | Optical system for fast three-dimensional imaging | |
| CN107121419B (en) | Non-invasive imaging method and device | |
| CN116503258A (en) | Super-resolution computing imaging method, device, electronic equipment and storage medium | |
| US20170000438A1 (en) | Energy modulated luminescence tomography | |
| Sathiya Priya et al. | Enhanced retina blood vessel segmentation by super resolution generative adversarial networks based U-net | |
| CN105259155B (en) | A kind of translucent imaging method of fast noninvasive and device | |
| US12112452B2 (en) | Holographic ultra resolution imaging | |
| Anandasabapathy et al. | An optical, endoscopic brush for high-yield diagnostics in esophageal cancer | |
| CN110044474A (en) | A kind of relevance imaging optimization method based on interpolation algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |