CN101803908A - Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method - Google Patents

Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method Download PDF

Info

Publication number
CN101803908A
CN101803908A CN 201010115512 CN201010115512A CN101803908A CN 101803908 A CN101803908 A CN 101803908A CN 201010115512 CN201010115512 CN 201010115512 CN 201010115512 A CN201010115512 A CN 201010115512A CN 101803908 A CN101803908 A CN 101803908A
Authority
CN
China
Prior art keywords
sample
reference arm
arm
optical fiber
fiber coupler
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.)
Pending
Application number
CN 201010115512
Other languages
Chinese (zh)
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN 201010115512 priority Critical patent/CN101803908A/en
Publication of CN101803908A publication Critical patent/CN101803908A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention discloses a dispersive modulation-based non-mirror image optimal frequency domain imaging system and a dispersive modulation-based non-mirror image optimal frequency domain imaging method. In the system, a transmission-type optimal scanning delay line is arranged in a reference arm for modulating dispersion so as to ensure the identity of the axial positions of a sample under different dispersion states; by changing a rotation angle of an oscillating mirror in the transmission optimal scanning delay line, the rapid modulation of dispersion is realized, and two group of interference spectrum signals of the same sample are acquired under two different dispersion states; the two group of acquired interference spectrum signals are respectively multiplied by corresponding dispersion compensation factors so as to accurately compensate the dispersion of imaginary part subjected to corresponding Fourier transform; subtraction operation is performed on the two group of interference spectrum signals subjected to dispersion compensation, then, the imaginary part of the corresponding repeated reflection signal disappears, and the real part still contains the dispersion factor; and finally, the interference spectrum signals subjected to subtraction are subjected to dispersive compensation again and then subjected to Fourier inverse transformation so as to obtain the actual reflection signal of the sample for non-mirror image reconstruction of the sample image.

Description

Based on synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion and method
Technical field
The present invention relates to optical coherent chromatographic imaging (OCT) technology, relate in particular to a kind of based on synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion and method.
Background technology
Optical coherence tomography (Optical Coherence Tomography, abbreviation OCT) imaging technique is a kind of novel optical image technology, can carry out non-intruding, noncontact, high-resolution imaging in vivo to the organizational structure and the physiological function of tested live body sample interior, at the early diagnosis of disease with in body biopsy field extensive application prospect.
Optimal frequency domain imaging system is a kind of pattern of optical coherence tomography system, by adopting high speed frequency-sweeping laser source and the axial interference spectrum signal of point probe collected specimens, reference arm is made of fixed plane mirror, do not need to carry out axial scan, by axial interference spectrum signal is carried out the axial depth information that inverse fourier transform can obtain sample, have at a high speed and highly sensitive characteristics.But also there is its inherent defect in the optimal frequency domain imaging, when the interference signal of sample depth information is carried in collection, also collect mutual interference signal between each layer of sample, each layer of sample itself from coherent interference signal, reference light itself from coherent noises such as coherent interference signals.And owing to what collect is the real part of interference spectrum signal, rather than plural interference spectrum signal, it is Hermitian conjugate that the real part of this interference spectrum signal is carried out the result that inverse fourier transform obtains, and has caused having produced in image being superimposed upon on the sample real number picture about the zero complete symmetric complex conjugate picture in light path position.
In order to differentiate the real number picture of sample, usually by regulating the reference arm light path 1 light path point is moved on to outside the sample surfaces, can make real number picture and complex conjugate picture not overlapping on image like this, but because near the fringe visibility zero light path is the highest, be that image sensitivity is the highest, the way that zero light path is removed in employing causes highly sensitive image-region can't obtain utilizing, and because zero light path is positioned at outside the sample, has caused the imaging depth of OFDI system only to utilize half.Eliminate the complex conjugate picture of Fourier domain OCT, near the high sensitivity zone can better utilization zero light path, and make imaging depth expand one times, external a lot of scientific research institutions have all carried out the research of this respect.People such as M.Wojtkowski propose to utilize piezoelectric ceramic actuator to move the method for the reflecting mirror of reference arm, the Zhongping Chen group in Irving branch school, University of California proposes to adopt the method for electro-optic phase modulator, Duke University Izatt group proposes to adopt the method for 3 * 3 fiber couplers, by between adjacent axial interference spectrum signal, introducing fixed additive phase, adopt plural interference spectrum recovery algorithms to reconstruct the plural form of interference spectrum signal, carry out inverse fourier transform again, thereby eliminate the complex conjugate picture.Medical college G.J.Tearney group of Harvard University proposes to adopt the acousto-optic phase-modulator that the interference spectrum signal is carried out the method for carrier frequency and the method for polarization encoder is removed the complex conjugate picture.People such as Hofer propose to adopt chromatic dispersion material that chromatic dispersion is provided and eliminate the method for complex conjugate picture with complicated iterative algorithm, and people such as S.Witte remove the complex conjugate picture to the elimination spike algorithm that also adopts chromatic dispersion material and carry out the chromatic dispersion coding and propose to simplify.
Above-mentioned these methods, all there is its inherent defect, as utilizes piezoelectric ceramic actuator to carry out the multistep phase-moving method need to take multiple measurements same position, reduce image taking speed, and, be subjected to of the influence of various environmental perturbations easily to phase place to the stability requirement height of phase place; Utilize the method for electro-optic phase modulator harmony optical phase modulator need introduce more complicated and expensive instrument and equipment, and the system data picking rate has been proposed harsh requirement; Utilizing the method for 3 * 3 fiber couplers to be subjected to temperature easily causes and recovers the inaccurate of plural interference spectrum, the complex conjugate suppression ratio that influence is whole the influence of the coefficient of coup; The existing method of utilizing chromatic dispersion to remove the complex conjugate picture that proposes need rely on complicated iterative algorithm, and the simplification of algorithm has also been reduced the complex conjugate suppression ratio.Therefore be necessary to study and be easy to realize and the mirror method that disappears that the complex conjugate suppression ratio is high.
Summary of the invention
The objective of the invention is at the deficiencies in the prior art, provide a kind of based on synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion and method.Reference arm at optimal frequency domain imaging system is provided with the transmission-type optical scan delay-line, realizes that the synthetic light path of keeping reference arm simultaneously of chromatic dispersion is constant, guarantees the homogeneity of sample axial location under different chromatic dispersion states.By the anglec of rotation of galvanometer in the quick change transmission-type optical scan delay-line, obtain two groups of interference spectrum signals of same sample under two kinds of chromatic dispersion states.By the multiplied by corresponding dispersion compensation factor chromatic dispersion of corresponding imaginary part is compensated, implement dispersion compensation once more after two groups of interference spectrums behind the dispersion compensation are subtracted each other, carry out inverse fourier transform then, just can obtain the real reflected signal of sample, the no mirror image that is used for sample image is rebuild, and makes one times of system imaging degree of depth expansion.
The objective of the invention is to be achieved through the following technical solutions:
A kind of based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion, it comprises: swept light source, first broadband optical fiber coupler, second broadband optical fiber coupler, sample ami light circulator, reference arm light circulator, sample arm Polarization Controller, sample arm collimating mirror, reference arm Polarization Controller, reference arm collimating mirror, reference arm condenser lens, reference arm plane mirror, sample arm scanning galvanometer, sample arm condenser lens, the 3rd broadband optical fiber coupler, transmission-type optical scan delay-line, balance detection device, Mach-Zehnder interferometers, data collecting card, computer.Wherein, swept light source links to each other with first broadband optical fiber coupler; First broadband optical fiber coupler connects second broadband optical fiber coupler and Mach-Zehnder interferometers respectively; Second broadband optical fiber coupler connects sample ami light circulator and reference arm light circulator respectively; The sample ami light circulator links to each other successively with sample arm Polarization Controller, sample arm collimating mirror; The collimated light beam degree of being of sample arm scanning galvanometer and the outgoing of sample arm collimating mirror is placed, and the light beam of turning back is radiated on the sample by the sample arm condenser lens; The reference arm light circulator links to each other successively with reference arm Polarization Controller, reference arm collimating mirror, reference arm condenser lens, reference arm plane mirror; The output port of reference arm light circulator is connected with the transmission-type optical scan delay-line again, the output port of transmission-type optical scan delay-line and sample ami light circulator links to each other with the 3rd broadband optical fiber coupler respectively, and two output ports of the 3rd broadband optical fiber coupler link to each other with two input ports of balance detection device; The output port of balance detection device, Mach-Zehnder interferometers links to each other with data collecting card respectively, and data collecting card links to each other with computer.
A kind of formation method of using above-mentioned based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion may further comprise the steps:
(1) the transmission-type optical scan delay-line is set in the reference arm of optimal frequency domain imaging system, looses in order to color to be provided;
(2) regulate the transmission-type optical scan delay-line, produce the chromatic dispersion that changes with galvanometer 15 anglecs of rotation, and it is constant to keep the reference arm light path;
(3) gather interference spectrum signal under two kinds of chromatic dispersion states, eliminate mirror image with dispersion compensation and the algorithm that subtracts each other, the acquisition gamut does not have the sample image of mirror image.
Compare with background technology, the present invention has following technique effect:
1, quick image cancellation algorithm.In the transmission-type optical scan delay-line, by changing the anglec of rotation of galvanometer, realize fast modulation to chromatic dispersion, obtain two groups of interference spectrum signals of same sample under two kinds of chromatic dispersion states thus.By these two groups of interference spectrum signal dispersion are compensated, reduce the imaginary part of common interflection signal, then the dispersion compensation second time is carried out in the residual dispersion of the interference spectrum signal of corresponding real part, can obtain plural interference spectrum signal, at last this plural interference spectrum signal is carried out inverse fourier transform, promptly obtain the optical coherence tomography image that gamut does not have mirror image.
2, can guarantee to realize under the constant condition of reference arm light path the fast modulation of chromatic dispersion.By the transversal displacement y between galvanometer rotating shaft and the fourier transform lens optical axis in the adjusting transmission-type optical scan delay-line 0, keep the reference arm light path constant in the time of can being implemented in the galvanometer anglec of rotation change in the transmission-type scanning delay line; Regulate the angle theta between grating normal and the optical axis g, can realize the fast modulation of chromatic dispersion by changing the anglec of rotation of galvanometer in the transmission-type optical scan delay-line, the light path of keeping reference arm is constant, guarantees the homogeneity of sample axial location under different chromatic dispersion states.Realize zero group delay and chromatic dispersion modulation simultaneously by introducing the transmission-type optical scan delay-line, be easy to realize, and guaranteed the compactedness and the reliability of optimal frequency domain imaging system.
Description of drawings
Fig. 1 is the system schematic of the specific embodiment based on the synthetic non-mirror image optimal frequency domain formation method of chromatic dispersion of the present invention;
Fig. 2 is the structural representation based on the transmission-type optical scan delay-line in the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion of the present invention;
Fig. 3 is the sequencing contro figure based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion of the present invention;
Fig. 4 is the algorithm flow chart based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion of the present invention;
Among the figure: 1, swept light source, 2, broadband optical fiber coupler, 3, optical circulator, 4, Polarization Controller, 5, the sample arm collimating mirror, 6, the sample arm scanning galvanometer, 7, condenser lens, 8, sample, 9, the reference arm collimating mirror, 10, condenser lens, 11, plane mirror, 12, collimating mirror, 13, balzed grating,, 14, fourier transform lens, 15, galvanometer, 16, corner cube prism, 17, receive reflecting mirror, 18, receive collimating mirror, 19, the balance detection device, 20, Mach-Zehnder interferometers, 21, data collecting card, 22, computer, 23, the transmission-type optical scan delay-line.
The specific embodiment
The present invention is based on the synthetic non-mirror image optimal frequency domain formation method of chromatic dispersion, reference arm at optimal frequency domain imaging system is provided with the transmission-type optical scan delay-line, realize that the synthetic light path of keeping reference arm simultaneously of chromatic dispersion is constant, guarantee the homogeneity of sample axial location under different chromatic dispersion states.By the anglec of rotation of galvanometer in the quick change transmission-type optical scan delay-line, obtain two groups of interference spectrum signals of same sample under two kinds of chromatic dispersion states.By the multiplied by corresponding dispersion compensation factor chromatic dispersion of corresponding imaginary part is compensated, implement dispersion compensation once more after two groups of interference spectrums behind the dispersion compensation are subtracted each other, carry out inverse fourier transform then, just can obtain the real reflected signal of sample, rebuild with the no mirror image that is used for sample image, make one times of system imaging degree of depth expansion.Its concrete steps are as follows:
1, by the galvanometer in the transmission-type optical scan delay-line in the synchronizing sequential circuit control reference arm, the scanheads of sample arm and the collection of data collecting card are synchronous, gather the interference spectrum signal of determining sequential.
2, by regulating the angle of grating normal and optical axis in the transmission-type optical scan delay-line, generation can be with the chromatic dispersion of galvanometer anglec of rotation change, by regulating the transversal displacement of galvanometer rotating shaft and optical axis, the light path that guarantees reference arm when the galvanometer anglec of rotation changes is constant, guarantees the homogeneity of sample axial location under different chromatic dispersion states.
3, in the transmission-type optical scan delay-line in the one-period of vibration mirror scanning, gather the interference spectrum signal of two groups of corresponding different galvanometer anglecs of rotation by synchronizing sequential circuit control data capture card, obtain two groups of interference spectrum signals of same sample under two kinds of chromatic dispersion states thus, by implementing the interference spectrum behind the dispersion compensation is subtracted each other and carries out once more the algorithm of dispersion compensation, reconstruct the complex expression of interference spectrum signal, by inverse fourier transform, can reconstruct the gamut optical coherence tomography image of no mirror image at last.
The present invention is based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion, comprise swept light source, broadband optical fiber coupler, optical circulator, Polarization Controller, the sample arm that constitutes by collimating mirror, scanning galvanometer, two gummed achromat, sample, the reference arm that constitutes by collimating mirror, two gummed achromat, plane mirror, transmission-type optical scan delay-line, by the demarcation light path that Mach-Zehnder interferometers constitutes, balance detection device, data collecting card, computer.The low-coherent light that comes out from swept light source, enter demarcation light path and main interference instrument light path respectively through broadband optical fiber coupler, the light that enters the demarcation light path produces a road sign through Mach-Zehnder interferometers and decides signal, the light that enters main interference instrument light path enters reference arm and sample arm respectively after the broadband optical fiber coupler beam split, the coherent light of two arms interferes back formation interference signal to be surveyed through the balance detection device at the broadband optical fiber coupler place, interference signal and demarcation signal are through the data collecting card collection, handle at computer at last, produce the interference signal of plural form, by the inverse fourier transform reconstructed image.
The transmission-type optical scan delay-line is made up of collimating mirror, balzed grating,, fourier transform lens, galvanometer, corner cube prism, reception reflecting mirror, reception collimating mirror, and wherein there is an adjustable angle theta in the optical axis of the normal of balzed grating, and fourier transform lens g, balzed grating, is positioned on the front focal plane of fourier transform lens, and galvanometer is positioned on the back focal plane of fourier transform lens, and the rotating shaft of galvanometer and the spacing of optical axis are an adjustable side-play amount y 0
The present invention is further illustrated below in conjunction with drawings and Examples, and it is more obvious that purpose of the present invention and effect will become.
Figure 1 shows that based on the structure chart of a specific embodiment of the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion, comprise swept light source 1, first broadband optical fiber coupler 2, second broadband optical fiber coupler 3, sample ami light circulator 4, reference arm light circulator 5, sample arm Polarization Controller 6, sample arm collimating mirror 7, reference arm Polarization Controller 8, reference arm collimating mirror 9, reference arm condenser lens 10, reference arm plane mirror 11, sample arm scanning galvanometer 19, sample arm condenser lens 20, sample 21, the 3rd broadband optical fiber coupler 22, transmission-type optical scan delay-line 23, balance detection device 24, Mach-Zehnder interferometers 25, data collecting card 26, computer 27.
Wherein, swept light source 1 links to each other with first broadband optical fiber coupler 2; First broadband optical fiber coupler 2 connects second broadband optical fiber coupler 3 and Mach-Zehnder interferometers 25 respectively; Second broadband optical fiber coupler 3 connects sample ami light circulator 4 and reference arm light circulator 5 respectively; Sample ami light circulator 4 links to each other successively with sample arm Polarization Controller 6, sample arm collimating mirror 7; Sample arm scanning galvanometer 19 is 45 degree placements with the collimated light beam of sample arm collimating mirror 7 outgoing, and the light beam of turning back is radiated on the sample 21 by sample arm condenser lens 20; Reference arm light circulator 5 links to each other successively with reference arm Polarization Controller 8, reference arm collimating mirror 9, reference arm condenser lens 10, reference arm plane mirror 11; The output port of reference arm light circulator 5 is connected with transmission-type optical scan delay-line 23 again, the output port of transmission-type optical scan delay-line 23 and sample ami light circulator 4 links to each other with the 3rd broadband optical fiber coupler 22 respectively, and two output ports of the 3rd broadband optical fiber coupler 22 link to each other with two input ports of balance detection device 24; The output port of balance detection device 24, Mach-Zehnder interferometers 25 links to each other with data collecting card 26 respectively, and data collecting card 26 links to each other with computer 27.
As shown in Figure 1, the low-coherent light that sends from swept light source 1, enter demarcation light path and main interference instrument light path respectively through first broadband optical fiber coupler 2, enter the light of demarcating light path and decide signal through Mach-Zehnder interferometers 25 generations one road sign, the light that enters main interference instrument light path enters reference arm and sample arm respectively after 3 beam split of second broadband optical fiber coupler, light in the reference arm enters transmission-type optical scan delay-line 23 after collimating through collimating mirror 12, reaching sample light that the 3rd broadband optical fiber coupler 22 and sample arm return converges the back and interferes, enter balance detection device 24, the demarcation signal that OCT interference spectrum signal that forms and Mach-Zehnder interferometers 25 produce is gathered by data collecting card 26 simultaneously, and these interference spectrum signals import into and carry out date processing and image reconstruction in the computer 27 at last.
Of the present invention based on the synthetic non-mirror image optimal frequency domain formation method of chromatic dispersion, may further comprise the steps:
1, the transmission-type optical scan delay-line is set in the reference arm of optimal frequency domain imaging system, in order to chromatic dispersion to be provided.
The light path that light among Fig. 1 in the sample arm is walked is a constant, uses L in the following derivation of equation SamExpression, the light path that the light in the reference arm is walked can be divided into two parts: the light path L except that transmission-type optical scan delay-line 23 RefExpression, the light path L in the transmission-type optical scan delay-line 23 RsodExpression, then the phase contrast of interference spectrum signal can be expressed as:
Φ=kL sam-kL ref-kL rsod (1)
Wherein k is light wave number and k=2 π/λ.Light is by transmission-type optical scan delay-line 23, and its phase change amount can be expressed as:
k·L rsod=φ R+k·l (2)
Wherein, l is galvanometer anglec of rotation γ in the transmission-type optical scan delay-line 23 and grating defocusing amount (L-f) light path when being zero.φ RAdditional phase place, then this additive phase φ under the non-vanishing state of galvanometer anglec of rotation γ in the transmission-type optical scan delay-line 23 and grating defocusing amount (L-f) RExpression is:
φ R ( ω ) = 4 y 0 ωγ c + 4 Lωγ sin α c cos α + 4 Lω cos α c + 4 fω c cos α + 8 πm ( L - f ) sin α p cos β - 8 fω c - - - ( 3 )
Wherein ω is light frequency and ω=kc, y 0Be the transversal displacements of galvanometer 15 rotating shafts apart from optical axis, γ is the anglec of rotation of galvanometer 15, and L is the distance between balzed grating, 13 and the fourier transform lens 14, and f is the focal length of fourier transform lens 14, and β is by p sin β=m (λ-λ 0) decision, p is a grating constant, m is that the order of diffraction is inferior, λ 0Be centre wavelength, α is by α=β+θ gDecision, θ gBe the angle that is become between the normal of balzed grating, 13 and the optical axis.φ R(ω) Taylor expansion is:
φ R ( ω ) = φ R ( ω 0 ) + φ R ′ ( ω 0 ) · ( ω - ω 0 ) + 1 2 ! φ R ′ ′ ( ω 0 ) · ( ω - ω 0 ) 2 + · · · - - - ( 4 )
Wherein, φ R0) be phase retardation, φ R' (ω 0) be group delay, φ R "0) be GVD.Its expression formula is respectively:
φ R ( ω 0 ) = 4 ω 0 y 0 γ c + 4 ω 0 ( L - f ) c - - - ( 5 )
φ R ′ ( ω 0 ) = 4 y 0 γ c - 8 πmLγ p ω 0 cos θ g + 4 ( L - f ) c - - - ( 6 )
φ R ′ ′ ( ω 0 ) = 16 π 2 m 2 c sin θ g Lγ p 2 ω 0 3 cos 3 θ g + 16 π 2 m 2 c ( L - f ) p 2 ω 0 3 cos 2 θ g - - - ( 7 )
Only consider φ R(ω) second order term of Taylor expansion can get interference spectrum phase contrast expression formula to (2), (4) formula substitution (1) formula and is:
Φ = ω c · [ L sam - L ref - l - φ R ′ ( ω 0 ) · c ] - [ φ R ( ω 0 ) - φ R ′ ( ω 0 ) · ω 0 + 1 2 φ R ′ ′ ( ω 0 ) · ( ω - ω 0 ) 2 ] - - - ( 8 )
2, regulate the transmission-type optical scan delay-line, produce the chromatic dispersion that changes with galvanometer 15 anglecs of rotation, and it is constant to keep the reference arm light path.
The transversal displacement of regulating between galvanometer 15 rotating shafts and optical axis is a following expression:
y 0 = m λ 0 L P · cos θ g - - - ( 9 ) Group delay φ then R' (ω 0) not the anglec of rotation γ with galvanometer 15 change.Again because when the light path l of the defocusing amount (L-f) of the anglec of rotation γ of the galvanometer in the transmission-type optical scan delay-line 23 15 and balzed grating, 13 when all being zero be constant, and the light path L except that transmission-type optical scan delay-line 23 in the reference arm RefBe constant,, be made as so first in (8) formula is only relevant with the sample depth z:
z=L sam-L ref-l-φ R′(ω 0)·c (10)
(8) second GVD φ in the formula R" (ω 0) change with the anglec of rotation γ of scanning galvanometer 15, promptly form the quantitative variable phase of the chromatic dispersion generation that causes by transmission-type optical scan delay-line 23 d, its expression formula is:
φ d = - [ φ R ( ω 0 ) - φ R ′ ( ω 0 ) · ω 0 + 1 2 φ R ′ ′ ( ω 0 ) · ( ω - ω 0 ) 2 ] - - - ( 11 )
The OCT interference signal that the reference light that sample light that sample arm is returned and reference arm return interferes the back to form, the interference spectrum signal expression that is collected by data collecting card 21 is:
Figure GSA00000045388800075
3, gather interference spectrum signal under two kinds of chromatic dispersion states, eliminate mirror image with dispersion compensation and the algorithm that subtracts each other, the acquisition gamut does not have the sample image of mirror image.
By the galvanometer 15 in the synchronous sequence signal control transmission-type optical scan delay-line 23, corresponding to galvanometer 15 different rotary angles, transmission-type optical scan delay-line 23 provides two kinds of different chromatic dispersions respectively, computer 21 produces two trigger collection signal triggering data collecting cards 20 respectively and carries out data acquisition, and two interference spectrum signals that collect are:
Figure GSA00000045388800076
Because chromatic dispersion phase place
Figure GSA00000045388800078
Can separate the parcel algorithm by phase place and try to achieve by place plane mirror in sample arm, each layer of sample relevant certainly
Figure GSA00000045388800081
Since intensity very I to ignore reference arm intensity | E r| 2Can remove DC terms by the balance detection device.To the interference spectrum signal times with dispersion compensation factors
Figure GSA00000045388800082
Can be compensated the interference spectrum signal I of chromatic dispersion c(ω):
Figure GSA00000045388800083
Figure GSA00000045388800084
The processing that subtracts of (6) formula that dispersion compensation is obtained and (7) formula, the interference spectrum signal of the corresponding interflection signal imaginary part that can be eliminated:
Figure GSA00000045388800085
To (8) formula poor divided by twice dispersion compensation factors square
Figure GSA00000045388800086
The interference spectrum signal that can obtain plural form is:
Figure GSA00000045388800087
At last resulting plural interference spectrum signal is carried out just can be eliminated an axial scan depth information of mirror image of Fourier transform, can realize the more good utilisation of zero light path place high sensitivity zone and make investigation depth enlarge one times.
Figure 2 shows that the sketch map of transmission-type optical scan delay-line.Wherein the optical axis of the normal of balzed grating, 13 and fourier transform lens 14 forms a θ gAngle, the light of reference arm enters the transmission-type optical scan delay-line after collimating through collimating mirror 12, collimated beam is via 13 beam split of reflection-type balzed grating,, each spectral components after the beam split focuses on the galvanometer 15 by fourier transform lens 14, after galvanometer 15 reflections, turn back to balzed grating, 13 by fourier transform lens 14, and merged into a branch of collimated light behind the diffraction once more by balzed grating, 13 and project on the corner cube prism 16, the light of process corner cube prism 16 is after staggering with the vertical direction in plane, optical table place, turn back to balzed grating, 13 along former direction, pass through fourier transform lens 14 and galvanometer 15 again, passed through fourier transform lens 14 and balzed grating, 13 once more after galvanometer 15 reflections, again merge into emergent light behind 13 4 diffraction of balzed grating,, the outgoing beam that is received reflecting mirror 17 reflections at last enters and receives collimating mirror 18.
Figure 3 shows that sequencing contro figure based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion.Computer 27 produces synchronous sequence signals and controls galvanometer 15 in the transmission-type optical scan delay-line 23 and the scanning galvanometer 19 in the sample arm.Computer 27 produces the scanning galvanometer 19 in one road triangular signal driving sample arm, and computer 27 produces the galvanometer 15 in one road square-wave signal driving transmission-type optical scan delay-line 23.Low level in the transmission-type optical scan delay-line in the square wave driving signal of 23 galvanometer 15 and the respectively corresponding two kinds of different chromatic dispersion states of high level.Utilize computer 27 to produce synchronous sequence, the ascent stage correspondence that the triangular wave of the scanning galvanometer 19 in sample arm drives signal in the transmission-type optical scan delay-line low level section in the square wave driving signal of 23 galvanometer 15, and the descending branch correspondence that the triangular wave of the scanning galvanometer 19 in sample arm drives signal in the transmission-type optical scan delay-line high level section in the square wave driving signal of 23 galvanometer 15.Computer 27 produces one tunnel trigger collection signal, is used for trigger data acquisition card 26 and carries out data acquisition, collects under two kinds of chromatic dispersion states corresponding to two groups of interference spectrum signals of same sample, imports at last and carries out date processing in the computer 27.
Figure 4 shows that algorithm flow chart based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion.Interference spectrum signal 1 and the interference spectrum signal 2 that collects be multiply by corresponding dispersion compensation factors respectively, make the imaginary part chromatic dispersion fine compensation behind the corresponding Fourier transform, obtain interference spectrum signal 1 behind the dispersion compensation and the interference spectrum signal 2 behind the dispersion compensation, with two groups of interference spectrum signal subtractions behind the above-mentioned dispersion compensation, the imaginary part of then corresponding interflection signal disappears, and real part still includes the chromatic dispersion factor.Interference spectrum signal after subtracting each other is implemented dispersion compensation once more, promptly obtain plural interference spectrum signal, carry out the real reflected signal that Fourier transform just can obtain sample at last, can be used for rebuilding the gamut optical coherence tomography image of no mirror image.

Claims (2)

1. one kind based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion, it is characterized in that it comprises: swept light source (1), first broadband optical fiber coupler (2), second broadband optical fiber coupler (3), sample ami light circulator (4), reference arm light circulator (5), sample arm Polarization Controller (6), sample arm collimating mirror (7), reference arm Polarization Controller (8), reference arm collimating mirror (9), reference arm condenser lens (10), reference arm plane mirror (11), sample arm scanning galvanometer (19), sample arm condenser lens (20), the 3rd broadband optical fiber coupler (22), transmission-type optical scan delay-line (23), balance detection device (24), Mach-Zehnder interferometers (25), data collecting card (26), computer (27).Wherein, swept light source (1) links to each other with first broadband optical fiber coupler (2); First broadband optical fiber coupler (2) connects second broadband optical fiber coupler (3) and Mach-Zehnder interferometers (25) respectively; Second broadband optical fiber coupler (3) connects sample ami light circulator (4) and reference arm light circulator (5) respectively; Sample ami light circulator (4) links to each other successively with sample arm Polarization Controller (6), sample arm collimating mirror (7); Sample arm scanning galvanometer (19) is the placement of (45) degree with the collimated light beam of sample arm collimating mirror (7) outgoing, and the light beam of turning back is radiated on the sample (21) by sample arm condenser lens (20); Reference arm light circulator (5) links to each other successively with reference arm Polarization Controller (8), reference arm collimating mirror (9), reference arm condenser lens (10), reference arm plane mirror (11); The output port of reference arm light circulator (5) is connected with transmission-type optical scan delay-line (23) again, the output port of transmission-type optical scan delay-line (23) and sample ami light circulator (4) links to each other with the 3rd broadband optical fiber coupler (22) respectively, and two output ports of the 3rd broadband optical fiber coupler (22) link to each other with two input ports of balance detection device (24); The output port of balance detection device (24), Mach-Zehnder interferometers (25) links to each other with data collecting card (26) respectively, and data collecting card (26) links to each other with computer (27).
2. an application rights requires 1 described formation method based on the synthetic non-mirror image optimal frequency domain imaging system of chromatic dispersion, may further comprise the steps:
(1) the transmission-type optical scan delay-line is set in the reference arm of optimal frequency domain imaging system, in order to chromatic dispersion to be provided.
(2) regulate the transmission-type optical scan delay-line, produce the chromatic dispersion that changes with galvanometer 15 anglecs of rotation, and it is constant to keep the reference arm light path.
(3) gather interference spectrum signal under two kinds of chromatic dispersion states, eliminate mirror image with dispersion compensation and the algorithm that subtracts each other, the acquisition gamut does not have the sample image of mirror image.
CN 201010115512 2010-03-01 2010-03-01 Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method Pending CN101803908A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010115512 CN101803908A (en) 2010-03-01 2010-03-01 Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010115512 CN101803908A (en) 2010-03-01 2010-03-01 Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method

Publications (1)

Publication Number Publication Date
CN101803908A true CN101803908A (en) 2010-08-18

Family

ID=42605988

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010115512 Pending CN101803908A (en) 2010-03-01 2010-03-01 Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method

Country Status (1)

Country Link
CN (1) CN101803908A (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062902A (en) * 2010-12-24 2011-05-18 华南师范大学 Mach-Zehnder interferometer based tunable flat-top multi-channel optical fiber filter
CN102151121A (en) * 2011-01-28 2011-08-17 浙江大学 Method and system for calibrating spectrum based on interference spectrum phase information
CN102188237A (en) * 2011-05-26 2011-09-21 浙江大学 Phase-multiplexing-based full-range sweep frequency OCT (Optical Coherence Tomography) imaging method and system
CN103070669A (en) * 2013-01-18 2013-05-01 杭州电子科技大学 Light-spectrum phase calibration system and method based on cascade Mach-Zehnder interferometer
CN103845039A (en) * 2014-02-07 2014-06-11 中国科学院上海光学精密机械研究所 Spectrograph for frequency domain OCT (Optical Coherence Tomography) system
CN104523239A (en) * 2015-01-12 2015-04-22 南京理工大学 Full-depth spectral domain optical coherent tomography device and method
CN104634282A (en) * 2015-02-28 2015-05-20 清华大学 Phase error compensation device for phase interferometry by optical probe
CN104688188A (en) * 2015-03-12 2015-06-10 清华大学 Spectral optical coherence imaging system based on optic computation
CN104990497A (en) * 2015-06-15 2015-10-21 上海交通大学 Rear-focal bias eccentric light beam method initiative defocusing amount detection apparatus without knife edge and method
CN105030201A (en) * 2015-07-27 2015-11-11 浙江大学 Correction method and system of sweep frequency OCT digital phases based on information on static regions
CN105559756A (en) * 2016-02-05 2016-05-11 浙江大学 Microangiography method and system based on total space modulation spectrum segmentation angle combining
CN105852816A (en) * 2016-05-16 2016-08-17 浙江大学 Angle composite speckle denoising method and system for full-channel modulation encoding
CN106510620A (en) * 2016-10-19 2017-03-22 苏州安莱光电科技有限公司 Full-optical non-contact vital sign monitoring device
CN106821323A (en) * 2017-02-21 2017-06-13 苏州大学 The adjustable optical coherence tomography method and system of image planes form
CN106841117A (en) * 2017-01-05 2017-06-13 上海理工大学 Ultrafast imaging method based on difference detecting method and Run- Length Coding
CN108287132A (en) * 2017-12-18 2018-07-17 首都师范大学 A kind of Terahertz asynchronous high-speed scanning system trigger signal generation device and method
CN109758119A (en) * 2019-01-28 2019-05-17 广东唯仁医疗科技有限公司 A kind of tumor-microvessel image processing apparatus
CN110006356A (en) * 2019-04-29 2019-07-12 中国科学院苏州生物医学工程技术研究所 Based on the real-time calibration apparatus and method in SS-OCT distance measurement system
CN110811287A (en) * 2018-08-10 2020-02-21 佛山市爱米吧科技有限公司 Anti-rotation-blockage control system of grain quantifying mechanism, electric cooker and anti-rotation-blockage control method
CN112244759A (en) * 2020-09-24 2021-01-22 温州医科大学 OCT imaging system suitable for whole-cycle dynamic determination of Schlemm's tube of room corner outflow channel
CN113040722A (en) * 2021-04-30 2021-06-29 电子科技大学 Method for improving frequency domain coherence tomography imaging depth
CN113503901A (en) * 2021-09-08 2021-10-15 武汉昊衡科技有限公司 Device and method for eliminating measurement signal jitter of white light interferometer
CN113984715A (en) * 2021-10-28 2022-01-28 上海盛晃光学技术有限公司 Coherence tomography device and method
CN114279568A (en) * 2022-03-04 2022-04-05 清华大学 Multispectral imaging method, device and equipment for encoding compression based on chromatic dispersion

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201328803Y (en) * 2008-12-22 2009-10-21 浙江大学 Ultra-wideband spectrum detecting system for OCT (optical coherence tomography) of spectral domain with ultrahigh resolution
CN101617935A (en) * 2009-08-06 2010-01-06 浙江大学 Method and system for wide-spectrum and high-resolution detection based on space-time light splitting in OCT

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201328803Y (en) * 2008-12-22 2009-10-21 浙江大学 Ultra-wideband spectrum detecting system for OCT (optical coherence tomography) of spectral domain with ultrahigh resolution
CN101617935A (en) * 2009-08-06 2010-01-06 浙江大学 Method and system for wide-spectrum and high-resolution detection based on space-time light splitting in OCT

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《浙江大学硕士学位论文》 20100101 徐磊 透射式快速扫描光学延迟线的研制和应用研究 23,25,26,31,41-47,52-53 1-2 , 2 *

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062902A (en) * 2010-12-24 2011-05-18 华南师范大学 Mach-Zehnder interferometer based tunable flat-top multi-channel optical fiber filter
CN102151121A (en) * 2011-01-28 2011-08-17 浙江大学 Method and system for calibrating spectrum based on interference spectrum phase information
CN102151121B (en) * 2011-01-28 2012-08-15 浙江大学 Method and system for calibrating spectrum based on interference spectrum phase information
CN102188237A (en) * 2011-05-26 2011-09-21 浙江大学 Phase-multiplexing-based full-range sweep frequency OCT (Optical Coherence Tomography) imaging method and system
CN102188237B (en) * 2011-05-26 2012-11-14 浙江大学 Phase-multiplexing-based full-range sweep frequency OCT (Optical Coherence Tomography) imaging method and system
CN103070669A (en) * 2013-01-18 2013-05-01 杭州电子科技大学 Light-spectrum phase calibration system and method based on cascade Mach-Zehnder interferometer
CN103070669B (en) * 2013-01-18 2015-05-13 杭州电子科技大学 Light-spectrum phase calibration system and method based on cascade Mach-Zehnder interferometer
CN103845039B (en) * 2014-02-07 2015-09-16 中国科学院上海光学精密机械研究所 For the spectrogrph of frequency domain OCT system
CN103845039A (en) * 2014-02-07 2014-06-11 中国科学院上海光学精密机械研究所 Spectrograph for frequency domain OCT (Optical Coherence Tomography) system
CN104523239A (en) * 2015-01-12 2015-04-22 南京理工大学 Full-depth spectral domain optical coherent tomography device and method
CN104634282A (en) * 2015-02-28 2015-05-20 清华大学 Phase error compensation device for phase interferometry by optical probe
CN104688188A (en) * 2015-03-12 2015-06-10 清华大学 Spectral optical coherence imaging system based on optic computation
CN104990497A (en) * 2015-06-15 2015-10-21 上海交通大学 Rear-focal bias eccentric light beam method initiative defocusing amount detection apparatus without knife edge and method
CN104990497B (en) * 2015-06-15 2017-11-14 上海交通大学 Eccentric pencil method active defocus amount detector and method without burnt biasing after the edge of a knife
CN105030201A (en) * 2015-07-27 2015-11-11 浙江大学 Correction method and system of sweep frequency OCT digital phases based on information on static regions
CN105030201B (en) * 2015-07-27 2017-10-13 浙江大学 Frequency sweep OCT digit phases antidote and system based on static region information
CN105559756A (en) * 2016-02-05 2016-05-11 浙江大学 Microangiography method and system based on total space modulation spectrum segmentation angle combining
CN105559756B (en) * 2016-02-05 2019-11-15 浙江大学 Based on the compound microangiography method and system of total space modulation spectrum segmentation angle
CN105852816A (en) * 2016-05-16 2016-08-17 浙江大学 Angle composite speckle denoising method and system for full-channel modulation encoding
CN106510620A (en) * 2016-10-19 2017-03-22 苏州安莱光电科技有限公司 Full-optical non-contact vital sign monitoring device
CN106841117A (en) * 2017-01-05 2017-06-13 上海理工大学 Ultrafast imaging method based on difference detecting method and Run- Length Coding
CN106841117B (en) * 2017-01-05 2019-04-19 上海理工大学 Ultrafast imaging method based on difference detecting method and Run- Length Coding
CN106821323B (en) * 2017-02-21 2023-10-13 苏州大学 Optical coherence tomography method and system with adjustable image plane morphology
CN106821323A (en) * 2017-02-21 2017-06-13 苏州大学 The adjustable optical coherence tomography method and system of image planes form
CN108287132A (en) * 2017-12-18 2018-07-17 首都师范大学 A kind of Terahertz asynchronous high-speed scanning system trigger signal generation device and method
CN110811287B (en) * 2018-08-10 2021-02-09 佛山市爱米吧科技有限公司 Anti-rotation-blockage control system of grain quantifying mechanism, electric cooker and anti-rotation-blockage control method
CN110811287A (en) * 2018-08-10 2020-02-21 佛山市爱米吧科技有限公司 Anti-rotation-blockage control system of grain quantifying mechanism, electric cooker and anti-rotation-blockage control method
CN109758119A (en) * 2019-01-28 2019-05-17 广东唯仁医疗科技有限公司 A kind of tumor-microvessel image processing apparatus
CN110006356B (en) * 2019-04-29 2021-04-02 中国科学院苏州生物医学工程技术研究所 Real-time calibration device and method based on SS-OCT (SS-OCT) distance measurement system
CN110006356A (en) * 2019-04-29 2019-07-12 中国科学院苏州生物医学工程技术研究所 Based on the real-time calibration apparatus and method in SS-OCT distance measurement system
CN112244759A (en) * 2020-09-24 2021-01-22 温州医科大学 OCT imaging system suitable for whole-cycle dynamic determination of Schlemm's tube of room corner outflow channel
WO2022062050A1 (en) * 2020-09-24 2022-03-31 温州医科大学 Oct imaging system suitable for all-round dynamic determination of schlemm's canal of anterior chamber angle outflow channel
CN113040722A (en) * 2021-04-30 2021-06-29 电子科技大学 Method for improving frequency domain coherence tomography imaging depth
CN113040722B (en) * 2021-04-30 2022-03-15 电子科技大学 Method for improving frequency domain coherence tomography imaging depth
CN113503901A (en) * 2021-09-08 2021-10-15 武汉昊衡科技有限公司 Device and method for eliminating measurement signal jitter of white light interferometer
CN113984715A (en) * 2021-10-28 2022-01-28 上海盛晃光学技术有限公司 Coherence tomography device and method
CN114279568A (en) * 2022-03-04 2022-04-05 清华大学 Multispectral imaging method, device and equipment for encoding compression based on chromatic dispersion

Similar Documents

Publication Publication Date Title
CN101803908A (en) Dispersive modulation-based non-mirror image optimal frequency domain imaging system and method
CN102028454B (en) Mirror-image separation method and system based on group delay wavenumber carrier frequency
CN103344569B (en) Polarization complex frequency domain optical coherence tomography imaging method and system
US7742174B2 (en) Methods, systems and computer program products for removing undesired artifacts in fourier domain optical coherence tomography (FDOCT) systems using continuous phase modulation and related phase modulators
US6775007B2 (en) Frequency-encoded parallel OCT and associated systems and methods
CN102151121B (en) Method and system for calibrating spectrum based on interference spectrum phase information
CN102188237B (en) Phase-multiplexing-based full-range sweep frequency OCT (Optical Coherence Tomography) imaging method and system
CN100464696C (en) Spectral coverage OCT imaging method based on optical scanning delay line and the system
CN202027563U (en) Spectrum calibration system based on interference spectrum phase information
CN101832817B (en) Parallel complex frequency domain optical coherence tomography imaging method and system
US7697145B2 (en) System for fourier domain optical coherence tomography
CN101297750B (en) Complex spectral domain optical coherence tomography method and system
CN102657518B (en) Method of complex frequency-domain optical coherence tomography using differential sinusoidal phase modulation
CN103070669B (en) Light-spectrum phase calibration system and method based on cascade Mach-Zehnder interferometer
US7394546B2 (en) Method and apparatus for full phase interferometry
CN102818786B (en) Sinusoidal phase modulation parallel complex frequency domain optical coherence tomography imaging system and method
CN104706322B (en) A kind of frequency sweep optical coherence imaging systems based on optical oomputing
CN101793823A (en) Apparatus and method for rangings and noise reduction of low coherence interferometry (LCI) and optical coherence tomography (OCT) signals
CN2916623Y (en) Full depth detecting frequency domain optical interference chromatographic imaging device
CN103439295A (en) Full-range Fourier-domain Doppler optical coherence tomography method
CN105147241A (en) Method and system based on double-space carrier frequency technology for increasing OCT (optical coherence tomography) imaging depth
CN103792192B (en) Based on the polarization domain optical coherence tomography system of simple detector
CN201814557U (en) Mirror image-free optical frequency domain imaging system based on chromatic dispersion modulation
CN108931478A (en) Single acquisition non-dispersive phase shift whole-field optically coherent chromatographic imaging device and method
CN206497027U (en) Monochromatic light spectrometer polarization domain optical coherence tomography system based on photoswitch

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20100818