CN102920438A - High-resolution optical scanning holographic slice imaging method based on variable pupils - Google Patents

High-resolution optical scanning holographic slice imaging method based on variable pupils Download PDF

Info

Publication number
CN102920438A
CN102920438A CN2012104226436A CN201210422643A CN102920438A CN 102920438 A CN102920438 A CN 102920438A CN 2012104226436 A CN2012104226436 A CN 2012104226436A CN 201210422643 A CN201210422643 A CN 201210422643A CN 102920438 A CN102920438 A CN 102920438A
Authority
CN
China
Prior art keywords
psi
matrix
pupil
slice
optical scanning
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
Application number
CN2012104226436A
Other languages
Chinese (zh)
Other versions
CN102920438B (en
Inventor
欧海燕
王秉中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201210422643.6A priority Critical patent/CN102920438B/en
Publication of CN102920438A publication Critical patent/CN102920438A/en
Application granted granted Critical
Publication of CN102920438B publication Critical patent/CN102920438B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Holo Graphy (AREA)

Abstract

The invention discloses a high-resolution optical scanning holographic slice imaging method based on variable pupils, belongs to the field of optical scanning and mainly solves problems of large out-of-focus noises existing when images of any two-dimensional slice are reconstructed in the prior art. According to the method, a two-dimensional scanning mirror is used for controlling deflection of a first Fresnel zone plate so as to achieve primary two-dimensional scanning of an object to be detected to obtain a first matrix equation; after a second pupil is deflected, the two-dimensional scanning mirror controls deflection of a second Fresnel zone plate so as to achieve secondary two-dimensional scanning of the object to be detected to obtain a second matrix equation; and the two matrix equations are combined, and a conjugate gradient algorithm is introduced to achieve slice imaging. According to the high-resolution optical scanning holographic slice imaging method based on the variable pupils, by means of the technical scheme, high-accuracy slice imaging is achieved, the out-of-focus noises are greatly reduced, and the method is applicable to each optical field for acquiring slice imaging.

Description

A kind of high-resolution optical scanning holography slice imaging method based on variable pupil
Technical field
The invention belongs to the optical scanning field, specifically, relate to a kind of high-resolution optical scanning holography slice imaging method based on variable pupil.
Background technology
The optical scanning holographic technique, be called for short OSH, it is a kind of non-traditional imaging technique based on the scanning of Fresnel plate, namely by the high resolution three-dimensional imaging of 2-D optical scanning realization to target, it all has wide practical use in fields such as biomedical imaging, fluorescent object imaging, 3D hologram television system and optical remote sensings.
And the two-dimensional hologram that obtains by the optical scanning holographic technique, the three-dimensional information that has comprised complete object, therefore in the optical scanning holographic technique to a slice imaging that the important analysis treatment step is exactly object of object hologram, i.e. the image reconstruction of object random two-dimensional tangent plane.And the difficult point in the reconstruct of object random two-dimensional tangent plane picture is how to eliminate the noise from other aspects of object, i.e. the out of focus noise.Slice imaging is the inverse problem during a typical image is processed, and also is an ill-posed problem simultaneously.
' Optical Scanning Holography with MATLAB ' has proposed a kind of traditional slice imaging method to document, namely the hologram with object carries out convolution algorithm with treating the Fresnel plate conjugation at reformatted slices place, thereby realization slice imaging, but owing to can't eliminate noise isolation, so its application is subject to great restriction.
' Three-dimensional microscopy and sectional image reconstruction using optical scanning holography ' has introduced a kind of contrary imaging algorithm to document, this iterative algorithm can realize that axial resolution is the slice imaging of 1 millimeter, and energy establishment out of focus noise, but it can't realize good imaging under less axial dimension.
' Depth resolution enhancement in optical scanning holography with a dual-wavelength laser source ' has proposed a kind of method of utilizing dual laser to improve the slice imaging axial resolution to document, it utilizes output wavelength to be respectively the laser instrument of 632nm and 543nm, obtain two groups of object holograies, and then axial resolution is increased to about 2.5 microns, but because two different wave lengths of working simultaneously in optical system have been introduced larger noise, cause its practicality to be extremely restricted.
Summary of the invention
The object of the invention is to reduce the out of focus noise in the slice imaging, a kind of high-resolution optical scanning holography slice imaging method based on variable pupil is proposed, utilize two different Fresnel plates that same object under test is scanned, obtain two groups of holograies, introduced more system of linear equations thereby decide inverse problem for this discomfort of slice imaging, realized and high-resolution slice imaging.
The technical solution used in the present invention is as follows:
A kind of high-resolution optical scanning holography slice imaging method based on variable pupil may further comprise the steps:
The light that (1) first polarization beam apparatus sends same light source is divided into two bundles, and light beam forms plane wave by the first pupil, and the second bundle light produces through acousto-optic modulator
Figure DEST_PATH_585119DEST_PATH_IMAGE001
Translation after form spherical wave by the second pupil again, plane wave and the first spherical wave by the second polarization beam apparatus optically focused, are produced to interfere at object under test after the polymerization and form the first Fresnel plate;
(2) utilize the two-dimensional scan mirror to control the deflection of the first Fresnel plate, thereby realize the two-dimensional scan first time to object under test, obtain comprising slice information The first matrix equation;
(3) voltage of adjusting spatial light modulator makes the second pupil produce skew;
The light that (4) first polarization beam apparatus send same light source is divided into two bundles, and light beam forms spherical wave by the first pupil, and the second bundle light produces by acousto-optic modulator
Figure DEST_PATH_95046DEST_PATH_IMAGE001
Translation after form the second spherical wave by the second pupil after the skew again, plane wave and the second spherical wave by the second polarization beam apparatus optically focused, are produced to interfere on the object under test again after the polymerization and form the second Fresnel plate;
(5) utilize the two-dimensional scan mirror to control the deflection of the second Fresnel plate, thereby realize the two-dimensional scan second time to object under test, obtain comprising slice information The second matrix equation;
(6) the first matrix equation and the second matrix equation are integrated, made the slice imaging process be converted into a minimal linear equation, and according to conjugate gradient algorithms, solve slice information
Figure DEST_PATH_366945DEST_PATH_IMAGE002
Wherein, the first Fresnel plate is described in the described step (1):
Figure DEST_PATH_369667DEST_PATH_IMAGE003
（1）
Wherein x, y, zRepresent space coordinates, kBe the wave number of light, zBe the distance of object under test to the two-dimensional scan mirror.
Further, object under test is the set of two discrete slices in the described step (2), and the axial location of two sections is respectively
Figure DEST_PATH_269490DEST_PATH_IMAGE004
With
Figure DEST_PATH_897917DEST_PATH_IMAGE005
, the specific implementation that therefore obtains the first matrix equation is as follows:
(2a) object under test is carried out the two-dimensional scan first time, obtains the first two-dimensional hologram:
Figure DEST_PATH_199586DEST_PATH_IMAGE006
（2）
Complex function wherein
Figure DEST_PATH_373209DEST_PATH_IMAGE007
Be the amplitude information of this object under test, * represents two-dimensional convolution simultaneously;
(2b) the Fresnel plate is existed With
Figure DEST_PATH_192446DEST_PATH_IMAGE005
The value at place is converted to respectively matrix
Figure DEST_PATH_83042DEST_PATH_IMAGE008
With
Figure DEST_PATH_681427DEST_PATH_IMAGE009
(2c) with the first two-dimensional hologram and matrix
Figure DEST_PATH_290263DEST_PATH_IMAGE008
With
Figure DEST_PATH_526073DEST_PATH_IMAGE009
Combine and obtain the first matrix equation:
Figure DEST_PATH_536754DEST_PATH_IMAGE010
（3）
Wherein
Figure DEST_PATH_317759DEST_PATH_IMAGE011
Be white Gaussian noise, this white Gaussian noise is that length is
Figure DEST_PATH_476208DEST_PATH_IMAGE012
The one dimension vector matrix.
In order accurately to try to achieve slice information, behind the regulation voltage, the offset distance of pupil is
Figure DEST_PATH_187812DEST_PATH_IMAGE013
, therefore the second Fresnel plate is:
Figure DEST_PATH_866049DEST_PATH_IMAGE014
（4）。
Further, the specific implementation that obtains the second matrix in the described step (5) is identical with the specific implementation that step (2) obtains the first matrix, and therefore the second matrix equation is following formula:
Figure DEST_PATH_4907DEST_PATH_IMAGE015
（5）。
Further, solve slice information in the described step (6)
Figure DEST_PATH_650652DEST_PATH_IMAGE002
Method as follows:
(6a) the first matrix equation and the second matrix equation are integrated, are obtained:
Figure DEST_PATH_900367DEST_PATH_IMAGE016
（6）
(6b) following formula be converted into minimize linear equation:
（7）
Wherein || || represent second order ranks norm,
Figure DEST_PATH_8449DEST_PATH_IMAGE018
For penalty factor and
Figure DEST_PATH_141490DEST_PATH_IMAGE018
＞0, CBe Laplace operator, this solution that minimizes linear equation is expressed as:
Figure DEST_PATH_194896DEST_PATH_IMAGE019
（8）
Wherein
Figure DEST_PATH_847726DEST_PATH_IMAGE020
For
Figure DEST_PATH_328386DEST_PATH_IMAGE021
Conjugate transpose;
(4c) can minimize linear equation to this by the introducing conjugate gradient algorithms and find the solution, obtain slice information
Figure DEST_PATH_683144DEST_PATH_IMAGE002
Value.
Compared with prior art, the present invention has following beneficial effect:
(1) the present invention has used the optical scanning holographic technique, namely pass through the deflection realization of two Fresnel plates to the twice sweep of object under test, obtain its high-resolution three-dimensional imaging, Fresnel plate range of application comprises biology, machinery, optics, electricity etc., therefore the present invention is applicable to every field, and range of application is very wide;
(2) the present invention only by regulating the voltage of spatial light modulator, can make the position of the second pupil exist
Figure DEST_PATH_540241DEST_PATH_IMAGE022
Upper skew, and the pupil after the second bundle light is by skew can obtain new spherical wave, obtain the second Fresnel plate, compare with a run-down, obtained a two-dimensional hologram more, therefore decide inverse problem for this discomfort of slice imaging and introduce more system of linear equations, so that finding the solution slice information
Figure DEST_PATH_234528DEST_PATH_IMAGE002
The time its result who obtains further approach slice information
Figure DEST_PATH_699138DEST_PATH_IMAGE002
Actual value, this is the reason that fundamentally improves slice imaging resolution;
(3) the present invention does not need mobile object under test, only needs to regulate the voltage of spatial light modulator, is easy to realize, has very strong operability;
(4) the first matrix equation and the second matrix equation that will try to achieve of the present invention combines, and is translated into and minimizes linear equation, and can solve the high slice information of degree of accuracy in conjunction with conjugate gradient algorithms
Figure DEST_PATH_744455DEST_PATH_IMAGE002
Value, whole computational process is simple, practical operation is easy, has greatly simplified slice information
Figure DEST_PATH_201981DEST_PATH_IMAGE002
Solution procedure;
(5) the present invention accurately tries to achieve slice information
Figure DEST_PATH_750774DEST_PATH_IMAGE002
Process in simple to operate, implementation is easy, has very strong practicality and operability, be fit to promote the use of.
Description of drawings
Fig. 1 is structural representation of the present invention.
Fig. 2 is testing sample section schematic diagram among the present invention-embodiment.
Fig. 3 is pupil different among the present invention-embodiment
Figure DEST_PATH_651865DEST_PATH_IMAGE023
Corresponding Fresnel figure.
Fig. 4 is the cosine hologram that for the first time scanning obtains among the present invention-embodiment.
Fig. 5 is the sinusoidal hologram that for the first time scanning obtains among the present invention-embodiment.
Fig. 6 is the cosine hologram that for the second time scanning obtains among the present invention-embodiment.
Fig. 7 is the sinusoidal hologram that for the second time scanning obtains among the present invention-embodiment.
Fig. 8 adopts diverse ways to exist among the present invention-embodiment
Figure DEST_PATH_918898DEST_PATH_IMAGE004
With
Figure DEST_PATH_117798DEST_PATH_IMAGE005
The slice imaging figure that the place obtains.
The specific embodiment
The invention will be further described below in conjunction with accompanying drawing and embodiment, and embodiments of the present invention include but not limited to the following example.
Embodiment
As shown in Figure 2, for the simplification problem, Fig. 2 is object under test schematic diagram among the present invention, and this object under test only comprises two slice information.In the optical scanning first time, the axial location of two sections of object is respectively
Figure DEST_PATH_583415DEST_PATH_IMAGE024
,
Figure DEST_PATH_842358DEST_PATH_IMAGE025
,
Figure DEST_PATH_409736DEST_PATH_IMAGE004
With
Figure DEST_PATH_146748DEST_PATH_IMAGE005
Axial distance be 10nm, and each section is of a size of
Figure DEST_PATH_732450DEST_PATH_IMAGE026
, matrix size is , wherein can realize scanning process according to step as follows:
Step 1 variable pupil is set to
Figure DEST_PATH_951390DEST_PATH_IMAGE028
, obtain first Fresnel plate FZP 1, object is carried out the scanning first time
(1) as shown in Figure 1, the angular frequency that is sent by same light source is that the light of ω is divided into two bundles by the first polarization beam apparatus BS1, and wherein light beam is by the first pupil
Figure DEST_PATH_492093DEST_PATH_IMAGE029
Form plane wave Pass through again the first pupil after the frequency displacement of light beam process acousto-optic modulator generation Ω
Figure DEST_PATH_533047DEST_PATH_IMAGE031
, the second pupil
Figure DEST_PATH_261969DEST_PATH_IMAGE031
Be variable pupil, realized by transmissive spatial light modulator, for the first time during scanning samples, will Be set to Thereby, form the first spherical wave.Use therein single wavelength light source centre wavelength is 632nm.
(2) the processed light of this two bundle is combined through the second polarization beam apparatus BS2, produces to interfere at testee to form the first Fresnel plate FZP 1, then utilize the deflection of this two-dimensional scan mirror 2D Scanning control TD-FZP, thereby realize the two-dimensional scan first time to three-dimensional body.
When object under test was scanned, the photodetector reception propagated into the light wave of test surface, and produced outer spill current output, and photogenerated current is processed through electricity such as mixing, amplifications, and the generation demodulating information also is stored in the computer.
The space impulse response of this optical scanning holophotal system, namely the first Fresnel plate can be expressed as:
Figure DEST_PATH_908610DEST_PATH_IMAGE032
（1）
Wherein
Figure DEST_PATH_187145DEST_PATH_IMAGE033
Represent space coordinates, Wave number for light.Find out from (1) formula, for a certain axial location
Figure DEST_PATH_422134DEST_PATH_IMAGE035
, the first Fresnel plate be one about
Figure DEST_PATH_177732DEST_PATH_IMAGE036
Two-dimensional function.
Suppose complex function
Figure DEST_PATH_615666DEST_PATH_IMAGE037
Represent the amplitude information of object, this object can be expressed as through the first two-dimensional hologram that obtains after the optical system scans:
Figure DEST_PATH_98600DEST_PATH_IMAGE038
（9）
Wherein * represents two-dimensional convolution, if object is regarded as the set of series of discrete section, and can be to axial coordinate zCarry out discretization and process, be expressed as
Figure DEST_PATH_571170DEST_PATH_IMAGE039
, represent respectively the axial location at different sections places.The first two-dimensional hologram of (9) formula sign can be expressed as so:
Figure DEST_PATH_232089DEST_PATH_IMAGE040
（10）
Owing to only having two sections, following formula can be reduced to following form so:
（2）
Analyze for (2) formula is reduced to one group of linear equation, we will
Figure DEST_PATH_443945DEST_PATH_IMAGE042
With Be converted to respectively the one dimension vector matrix
Figure DEST_PATH_868421DEST_PATH_IMAGE044
With
Figure DEST_PATH_280948DEST_PATH_IMAGE045
Because the section of object under test is one
Figure DEST_PATH_105685DEST_PATH_IMAGE046
Matrix,
Figure DEST_PATH_287267DEST_PATH_IMAGE044
With
Figure DEST_PATH_555569DEST_PATH_IMAGE045
Then for length be
Figure DEST_PATH_455391DEST_PATH_IMAGE047
The one dimension vector matrix.Equally, the two-dimensional hologram of object
Figure DEST_PATH_21502DEST_PATH_IMAGE048
Also can be converted into a length is
Figure DEST_PATH_119908DEST_PATH_IMAGE049
The one dimension vector matrix
Figure DEST_PATH_746061DEST_PATH_IMAGE050
(2) convolution algorithm of formula can be expressed as matrix operations, and the first Fresnel plate exists for this reason With The value at place is respectively
Figure DEST_PATH_280662DEST_PATH_IMAGE053
Matrix
Figure DEST_PATH_77717DEST_PATH_IMAGE054
With
Figure DEST_PATH_499602DEST_PATH_IMAGE053
Matrix
Figure DEST_PATH_673095DEST_PATH_IMAGE055
, two shown in being constructed as follows respectively
Figure DEST_PATH_683776DEST_PATH_IMAGE056
Two-dimensional matrix
Figure DEST_PATH_714049DEST_PATH_IMAGE057
With
Figure DEST_PATH_810181DEST_PATH_IMAGE058
:
Figure DEST_PATH_69255DEST_PATH_IMAGE059
（11）
Figure DEST_PATH_262339DEST_PATH_IMAGE060
（12）
Then can be expressed as the first matrix equation by (11) formula and (12) formula:
Figure DEST_PATH_214246DEST_PATH_IMAGE061
（3）
Wherein
Figure DEST_PATH_797674DEST_PATH_IMAGE062
With
Figure DEST_PATH_109706DEST_PATH_IMAGE063
The white Gaussian noise of representative system is that length is N 2 The one dimension vector matrix.
So-called slice imaging, namely will from G 1 In recover slice information
Figure DEST_PATH_94980DEST_PATH_IMAGE002
Information, this is that a discomfort is decided inverse problem, so we recycle the second Fresnel plate same object under test is scanned, and improves the slice imaging axial resolution.
Step 2 utilizes the second Fresnel plate that same object under test is scanned
Regulate the voltage of spatial light modulator, will
Figure DEST_PATH_217788DEST_PATH_IMAGE064
Be set to:
Figure DEST_PATH_288512DEST_PATH_IMAGE065
（13）
Thereby obtain the spherical wave of center skew, wherein
Figure DEST_PATH_404236DEST_PATH_IMAGE066
With
Figure DEST_PATH_244016DEST_PATH_IMAGE067
Represent the second pupil
Figure DEST_PATH_537725DEST_PATH_IMAGE064
Spatial offset on the axle.The second spherical wave and plane wave are combined through the second polarization beam apparatus, produce at testee and interfere formation the second Fresnel plate FZP 2
Therefore the space impulse of this optical scanning holophotal system response, namely the second Fresnel plate can be expressed as:
Figure DEST_PATH_749580DEST_PATH_IMAGE068
（4）
(4) formula and the contrast of (1) formula can be found out, by changing the second pupil
Figure DEST_PATH_443867DEST_PATH_IMAGE069
, can utilize same set of optical system, design realizes different Fresnel plates.
By same object under test is scanned, can obtain the second object under test hologram, same process can be characterized by matrix equation:
Figure DEST_PATH_908477DEST_PATH_IMAGE070
（5）
Can find out from (6) formula, obtain again by scanning us for the second time
Figure DEST_PATH_953794DEST_PATH_IMAGE047
Individual system of linear equations.In the twice sweep, the Fresnel plate of record sample holographic information is different, namely , as seen utilizing second Fresnel plate is different to the system of linear equations that same object under test scans acquisition.We have just added more effectively system of linear equations for finding the solution of this ill-posed problem like this.
Step 3 utilizes the hologram that scans rear acquisition by twice Fresnel plate to carry out slice imaging.
So-called slice imaging is immediately known
Figure DEST_PATH_960113DEST_PATH_IMAGE072
Situation under, find the solution slice information
Figure DEST_PATH_861204DEST_PATH_IMAGE073
The the first rectangle equation and the second matrix equation that at first twice two-dimensional scan are obtained combine, and are expressed as:
Figure DEST_PATH_128237DEST_PATH_IMAGE074
（6）
Finding the solution of this problem can be converted into the following linear equation that minimizes, that is:
Figure DEST_PATH_123875DEST_PATH_IMAGE075
（7）
Wherein in (8) || || expression second order norm, λ〉0Be penalty factor, CIt is Laplace operator.The solution of this minimization problem can be expressed as:
Figure DEST_PATH_792754DEST_PATH_IMAGE076
（8）
Wherein
Figure DEST_PATH_864746DEST_PATH_IMAGE020
Be matrix
Figure DEST_PATH_353496DEST_PATH_IMAGE021
Conjugate transpose.By introducing conjugate gradient algorithms, can find the solution for (9) formula, and this convergence depends on matrix For the object that comprises two sections, this matrix can be expressed as:
Figure DEST_PATH_941790DEST_PATH_IMAGE078
（14）
By (14) formula as can be known, this matrix is positive definite symmetric matrices, then finds the solution by conjugate gradient algorithms, thereby obtains slice information
Figure DEST_PATH_919104DEST_PATH_IMAGE002
, wherein said conjugate gradient algorithms are prior art.
In the twice sweep, the Fresnel figure that different pupils obtain as shown in Figure 3.Wherein Fig. 4 is the cosine hologram that scans for the first time, and Fig. 5 is the sinusoidal hologram that scans for the first time, and Fig. 6 is the cosine hologram of the scanning second time, and Fig. 7 is the sinusoidal hologram of the scanning second time.By Fig. 4～Fig. 7 as seen, we can't obtain the slice information of any testing sample from hologram.
Compare three kinds of slice imaging methods: (1) traditional slice imaging method, namely carry out convolution with Fresnel zone plate conjugate function and the sample holograms of section place; (2) based on the fixing slice imaging method of pupil single sweep operation; (3) based on the slice imaging method of variable pupil twice sweep.Fig. 8 (a)-(f) has showed respectively the result who carries out slice imaging with three kinds of distinct methods.Can be found out by Fig. 8 (a)-(b), traditional slice imaging method not only can not be distinguished two sections, and has introduced larger out of focus noise; And Fig. 8 (c) although-(d) show can suppress to a certain extent the out of focus noise based on the fixing slice imaging method of pupil single sweep operation, can't reach higher 10nm axial resolution; And adopt method of the present invention, and then can distinguish axial distance fully and be two degree of depth sections of 10nm, its result this means shown in Fig. 8 (e)-(f), the inventive method has been brought up to 10nm with the axial resolution of optical scanning holography.
According to above-described embodiment, just can well realize the present invention.

Claims (6)

1. the high-resolution optical scanning holography slice imaging method based on variable pupil is characterized in that, may further comprise the steps:
The light that (1) first polarization beam apparatus sends same light source is divided into two bundles, light beam forms plane wave by the first pupil, form the first spherical wave by the second pupil again after the second translation of bundle light through acousto-optic modulator generation Ω, plane wave and the first spherical wave are passed through the second polarization beam apparatus optically focused, produce interference at object under test after the polymerization and form the first Fresnel plate;
(2) utilize the two-dimensional scan mirror to control the deflection of the first Fresnel plate, thereby realize the two-dimensional scan first time to object under test, obtain comprising the first matrix equation of slice information ψ;
(3) voltage of adjusting spatial light modulator makes the second pupil produce skew;
The light that (4) first polarization beam apparatus send same light source is divided into two bundles, light beam forms plane wave by the first pupil, form the second spherical wave by the second pupil after the skew again after the second translation of bundle light by acousto-optic modulator generation Ω, plane wave and the second spherical wave are passed through the second polarization beam apparatus optically focused, produce interference at object under test after the polymerization and form the second Fresnel plate;
(5) utilize the two-dimensional scan mirror to control the deflection of the second Fresnel plate, thereby realize the two-dimensional scan second time to object under test, obtain comprising the second matrix equation of slice information ψ;
(6) the first matrix equation and the second matrix equation are integrated, made the slice imaging process be converted into a minimal linear equation, and according to conjugate gradient algorithms, solve slice information ψ.
2. a kind of high-resolution optical scanning holography slice imaging method based on variable pupil according to claim 1 is characterized in that, the first Fresnel plate is described in the described step (1):
h 1 ( x , y , z ) = - j k 2 πz exp ( j k 2 z ( x 2 + y 2 ) ) - - - ( 1 )
X wherein, y, z represents space coordinates, and k is the wave number of light, and z is the distance that object under test arrives the two-dimensional scan mirror.
3. a kind of high-resolution optical scanning holography slice imaging method based on variable pupil according to claim 2 is characterized in that, object under test is the set of two discrete slices in the described step (2), and the axial location of two sections is respectively z 1And z 2, the specific implementation that therefore obtains the first matrix equation is as follows:
(2a) object under test is carried out the two-dimensional scan first time, obtains the first two-dimensional hologram:
Figure FDA00002329076400021
Complex function wherein
Figure FDA00002329076400022
Be the amplitude information of this object under test, * represents two-dimensional convolution simultaneously;
(2b) with the first Fresnel plate at z 1And z 2The value at place is converted to respectively matrix H 1(z 1) and H 1(z 2);
(2c) with the first two-dimensional hologram and matrix H 1(z 1) and H 1(z 2) combine and obtain the first matrix equation:
G 1 = H 1 ( z 1 ) ψ 1 + H 1 ( z 2 ) ψ 2 + n 1 = H 1 ( z 1 ) H 1 ( z 2 ) · ψ 1 ψ 2 + n 1 = H 1 ψ + n 1 - - - ( 3 )
N wherein 1Be white Gaussian noise, this white Gaussian noise is that length is N 2The one dimension vector matrix.
4. a kind of high-resolution optical scanning holography slice imaging method based on variable pupil according to claim 3 is characterized in that, behind the regulation voltage, the offset distance of pupil is x 0, y 0, therefore the second Fresnel plate is:
h 2 ( x , y , z ) = - j k 2 πz exp ( j k 2 z ( ( x - x 0 ) 2 + ( y - y 0 ) 2 ) ) - - - ( 4 ) .
5. a kind of high-resolution optical scanning holography slice imaging method based on variable pupil according to claim 4, it is characterized in that, the specific implementation that obtains the second matrix in the described step (5) is identical with the specific implementation that step (2) obtains the first matrix, and therefore the second matrix equation is following formula:
G 2 = H 2 ( z 1 ) ψ 1 + H 2 ( z 2 ) ψ 2 + n 2 = H 2 ( z 1 ) H 2 ( z 2 ) · ψ 1 ψ 2 + n 2 = H 2 ψ + n 2 - - - ( 5 ) .
6. a kind of high-resolution optical scanning holography slice imaging method based on variable pupil according to claim 5 is characterized in that, the method that solves slice information ψ in the described step (6) is as follows:
(6a) the first matrix equation and the second matrix equation are integrated, are obtained:
G = G 1 G 2 = H 1 ( z 1 ) H 1 ( z 2 ) H 2 ( z 1 ) H 2 ( z 2 ) ψ 1 ψ 2 + n 1 n 2 = Hψ + n ,
H = H 1 ( z 1 ) H 1 ( z 2 ) H 2 ( z 1 ) H 2 ( z 2 ) , ψ = ψ 1 ψ 2 , n = n 1 n 2 - - - ( 6 )
(6b) following formula be converted into minimize linear equation:
f ( ψ ) = | | Hψ - G | | 2 2 + λ | | Cψ | | 2 2 - - - ( 7 )
Wherein ‖ ‖ represents second order ranks norm, and λ is penalty factor and λ＞0, and C is Laplace operator, and this solution that minimizes linear equation is expressed as:
(H +H+λC +C)f(ψ)＝H +G （8）
H wherein +Conjugate transpose for H;
(6c) can minimize linear equation to this by the introducing conjugate gradient algorithms and find the solution, obtain the value of slice information ψ.
CN201210422643.6A 2012-10-30 2012-10-30 High-resolution optical scanning holographic slice imaging method based on variable pupils Expired - Fee Related CN102920438B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210422643.6A CN102920438B (en) 2012-10-30 2012-10-30 High-resolution optical scanning holographic slice imaging method based on variable pupils

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210422643.6A CN102920438B (en) 2012-10-30 2012-10-30 High-resolution optical scanning holographic slice imaging method based on variable pupils

Publications (2)

Publication Number Publication Date
CN102920438A true CN102920438A (en) 2013-02-13
CN102920438B CN102920438B (en) 2014-07-16

Family

ID=47635479

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210422643.6A Expired - Fee Related CN102920438B (en) 2012-10-30 2012-10-30 High-resolution optical scanning holographic slice imaging method based on variable pupils

Country Status (1)

Country Link
CN (1) CN102920438B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104159094A (en) * 2014-07-09 2014-11-19 四川大学 Method for improving optical scanning holographic tomography effect
CN104614970A (en) * 2015-02-15 2015-05-13 电子科技大学 Optical scanning holographic image edge extracting method based on double-hole pupil
CN105204311A (en) * 2015-07-06 2015-12-30 电子科技大学 Gaussian apodization based optical scanning holographic edge detection method
CN105404128A (en) * 2016-01-05 2016-03-16 中国科学院光电研究院 Multiframe phase shift digital holography method and device
CN105549370A (en) * 2016-02-23 2016-05-04 中国科学院光电研究院 Synthetic aperture digital holographic method and device based on multichannel low-frequency heterodyne
CN107015466A (en) * 2017-04-19 2017-08-04 电子科技大学 The holographic one-point positioning method of optical scanner based on TR MUSIC algorithms
CN107240074A (en) * 2017-05-15 2017-10-10 电子科技大学 Based on the hot-tempered sound removing method of the two-dimentional optimal defocus of Entropic method and genetic algorithm
CN108613110A (en) * 2016-12-22 2018-10-02 深圳市绎立锐光科技开发有限公司 A kind of light-source system, automobile lighting system and control method
CN109884869A (en) * 2019-04-22 2019-06-14 重庆邮电大学 One kind calculating holographic imaging systems and method based on dicoria heterodyne
CN110907415A (en) * 2019-11-01 2020-03-24 浙江大学 Three-dimensional sub-ten nanometer positioning method and device based on parallel detection

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6038041A (en) * 1997-09-17 2000-03-14 Virginia Tech Intellectual Properties, Inc. Three-dimensional holographic fluorescence microscope
US20060132799A1 (en) * 2001-06-29 2006-06-22 Franck Dubois Digital holographic microscope
CN201101518Y (en) * 2007-07-03 2008-08-20 浙江大学 Common channel type endoscopic optical coherent chromatography imaging system
CN101430428A (en) * 2008-11-25 2009-05-13 中国科学院微电子研究所 Super-resolution Fresnel zone plate
US20090181396A1 (en) * 2008-01-10 2009-07-16 Pacific Biosciences Of California, Inc. Methods and systems for analysis of fluorescent reactions with modulated excitation
CN101694532A (en) * 2009-10-22 2010-04-14 中国科学院光电技术研究所 Phase type photon sieve
CN102393569A (en) * 2011-11-25 2012-03-28 中国科学院上海光学精密机械研究所 Polarized optical scanner

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6038041A (en) * 1997-09-17 2000-03-14 Virginia Tech Intellectual Properties, Inc. Three-dimensional holographic fluorescence microscope
US20060132799A1 (en) * 2001-06-29 2006-06-22 Franck Dubois Digital holographic microscope
CN201101518Y (en) * 2007-07-03 2008-08-20 浙江大学 Common channel type endoscopic optical coherent chromatography imaging system
US20090181396A1 (en) * 2008-01-10 2009-07-16 Pacific Biosciences Of California, Inc. Methods and systems for analysis of fluorescent reactions with modulated excitation
CN101430428A (en) * 2008-11-25 2009-05-13 中国科学院微电子研究所 Super-resolution Fresnel zone plate
CN101694532A (en) * 2009-10-22 2010-04-14 中国科学院光电技术研究所 Phase type photon sieve
CN102393569A (en) * 2011-11-25 2012-03-28 中国科学院上海光学精密机械研究所 Polarized optical scanner

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
EDMUND Y.LAM ET AL.: "Three-dimensional microscopy and sectional image reconstruction using optical scanning holography", 《APPLIED OPTICS》, vol. 48, no. 34, 1 December 2009 (2009-12-01), pages 113 - 119, XP001550264, DOI: doi:10.1364/AO.48.00H113 *
孙萍等: "基于涅耳波带板扫描全息术的光学层析成像", 《光学学报》, vol. 24, no. 1, 31 January 2004 (2004-01-31), pages 110 - 114 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104159094A (en) * 2014-07-09 2014-11-19 四川大学 Method for improving optical scanning holographic tomography effect
CN104614970A (en) * 2015-02-15 2015-05-13 电子科技大学 Optical scanning holographic image edge extracting method based on double-hole pupil
CN105204311B (en) * 2015-07-06 2018-05-18 电子科技大学 A kind of optical scanner holography edge detection method based on Gauss apodization
CN105204311A (en) * 2015-07-06 2015-12-30 电子科技大学 Gaussian apodization based optical scanning holographic edge detection method
CN105404128A (en) * 2016-01-05 2016-03-16 中国科学院光电研究院 Multiframe phase shift digital holography method and device
CN105404128B (en) * 2016-01-05 2018-04-13 中国科学院光电研究院 Multiframe phase-shifted digital holographic method and device
CN105549370A (en) * 2016-02-23 2016-05-04 中国科学院光电研究院 Synthetic aperture digital holographic method and device based on multichannel low-frequency heterodyne
CN105549370B (en) * 2016-02-23 2018-11-09 中国科学院光电研究院 Synthetic aperture digital hologram method and device based on multichannel low frequency heterodyne
CN108613110A (en) * 2016-12-22 2018-10-02 深圳市绎立锐光科技开发有限公司 A kind of light-source system, automobile lighting system and control method
CN107015466A (en) * 2017-04-19 2017-08-04 电子科技大学 The holographic one-point positioning method of optical scanner based on TR MUSIC algorithms
CN107015466B (en) * 2017-04-19 2019-07-19 电子科技大学 Optical scanner holography one-point positioning method based on TR-MUSIC algorithm
CN107240074A (en) * 2017-05-15 2017-10-10 电子科技大学 Based on the hot-tempered sound removing method of the two-dimentional optimal defocus of Entropic method and genetic algorithm
CN109884869A (en) * 2019-04-22 2019-06-14 重庆邮电大学 One kind calculating holographic imaging systems and method based on dicoria heterodyne
CN110907415A (en) * 2019-11-01 2020-03-24 浙江大学 Three-dimensional sub-ten nanometer positioning method and device based on parallel detection

Also Published As

Publication number Publication date
CN102920438B (en) 2014-07-16

Similar Documents

Publication Publication Date Title
CN102920438B (en) High-resolution optical scanning holographic slice imaging method based on variable pupils
CN102805613B (en) High-resolution optical scanning holographic slice imaging method based on two-time scanning
CN105204311B (en) A kind of optical scanner holography edge detection method based on Gauss apodization
CN103822577B (en) Single pixel Terahertz holographic imaging apparatus and method
CN108089425B (en) A method of the elimination optical scanner holography defocus noise based on deep learning
CN104459971B (en) A kind of Structured Illumination micro imaging system based on integrated photon chip
CN108761752A (en) Fourier overlapping associations imaging system, device and method
CN110675451B (en) Digital self-adaptive correction method and system based on phase space optics
CN105184295B (en) A kind of holoscan space length extracting method based on wavelet transformation and connected domain
CN103743707A (en) Microscopy method and device for controlled high-speed chromatography phase position
Yang et al. Synthetic Bessel light needle for extended depth-of-field microscopy
CN110360924A (en) A kind of laminated micro imaging method of double wave and system
Li et al. Interleaved segment correction achieves higher improvement factors in using genetic algorithm to optimize light focusing through scattering media
CN109188881A (en) A kind of THz wave digital hologram imaging method and system
CN103760689B (en) Far field beam focal spot position control method is organized in expection based on optical phased array more
CN109238131A (en) A kind of optical coherence tomography method and system of transverse direction super-resolution
CN114488513B (en) Full-vector modulation single-fiber high-signal-to-noise-ratio three-dimensional imaging method and device
CN107835074B (en) A method of eliminating accidental enciphering optical scanner holography defocus noise
CN108153132B (en) A kind of optical scanner holography self-focusing method based on mean value gradient function
CN110133849A (en) A kind of point spread function reconstruction method
CN113009801B (en) High-speed multi-directional line confocal digital holographic three-dimensional microscopic imaging method and device
CN107145053B (en) Optical scanner holography axially position method based on TR-MUSIC algorithm
Leng et al. Speckle suppression in digital holographic imaging with random phases and different wavelengths
Xia et al. Exploring angular-steering illumination-based eyebox expansion for holographic displays
Chen et al. Resolution enhancement with improved range Doppler algorithm in high numerical aperture OCT

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140716

Termination date: 20141030

EXPY Termination of patent right or utility model