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 PDFInfo
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims description 23
- 239000006185 dispersion Substances 0.000 claims abstract description 85
- 238000001228 spectrum Methods 0.000 claims abstract description 50
- 230000003287 optical effect Effects 0.000 claims description 64
- 239000013307 optical fiber Substances 0.000 claims description 35
- 230000010287 polarization Effects 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 abstract 1
- 230000009466 transformation Effects 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 75
- 238000012014 optical coherence tomography Methods 0.000 description 6
- 230000001427 coherent effect Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000001105 regulatory effect Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 208000007479 Orofaciodigital syndrome type 1 Diseases 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 201000003455 orofaciodigital syndrome I Diseases 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
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
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:
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:
Wherein, φ
R(ω
0) be phase retardation, φ
R' (ω
0) be group delay, φ
R "(ω
0) be GVD.Its expression formula is respectively:
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:
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:
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:
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:
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:
Because chromatic dispersion phase place
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
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
Can be compensated the interference spectrum signal I of chromatic dispersion
c(ω):
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:
To (8) formula poor divided by twice dispersion compensation factors square
The interference spectrum signal that can obtain plural form is:
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.
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 (26)
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 |
CN110907402A (en) * | 2019-12-24 | 2020-03-24 | 上海雄博精密仪器股份有限公司 | Time-domain optical coherence tomography 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 |
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 |
CN114931358A (en) * | 2022-04-06 | 2022-08-23 | 上海健康医学院 | Large depth of field micrometer resolution optical coherence tomography endoscopic probe |
Citations (2)
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 |
-
2010
- 2010-03-01 CN CN 201010115512 patent/CN101803908A/en active Pending
Patent Citations (2)
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)
Title |
---|
《浙江大学硕士学位论文》 20100101 徐磊 透射式快速扫描光学延迟线的研制和应用研究 23,25,26,31,41-47,52-53 1-2 , 2 * |
Cited By (41)
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 |
CN103845039A (en) * | 2014-02-07 | 2014-06-11 | 中国科学院上海光学精密机械研究所 | Spectrograph for frequency domain OCT (Optical Coherence Tomography) system |
CN103845039B (en) * | 2014-02-07 | 2015-09-16 | 中国科学院上海光学精密机械研究所 | For the spectrogrph of frequency domain OCT 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 |
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 |
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 |
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 |
CN106841117B (en) * | 2017-01-05 | 2019-04-19 | 上海理工大学 | Ultrafast imaging method based on difference detecting method and Run- Length Coding |
CN106841117A (en) * | 2017-01-05 | 2017-06-13 | 上海理工大学 | Ultrafast imaging method based on difference detecting method and Run- Length Coding |
CN106821323A (en) * | 2017-02-21 | 2017-06-13 | 苏州大学 | The adjustable optical coherence tomography method and system of image planes form |
CN106821323B (en) * | 2017-02-21 | 2023-10-13 | 苏州大学 | Optical coherence tomography method and system with adjustable image plane morphology |
CN108287132A (en) * | 2017-12-18 | 2018-07-17 | 首都师范大学 | A kind of Terahertz asynchronous high-speed scanning system trigger signal generation device and 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 |
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 |
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 |
CN110006356B (en) * | 2019-04-29 | 2021-04-02 | 中国科学院苏州生物医学工程技术研究所 | Real-time calibration device and method based on SS-OCT (SS-OCT) distance measurement system |
CN110907402A (en) * | 2019-12-24 | 2020-03-24 | 上海雄博精密仪器股份有限公司 | Time-domain optical coherence tomography system |
CN110907402B (en) * | 2019-12-24 | 2024-07-30 | 上海雄博精密仪器股份有限公司 | Time domain optical coherence tomography system |
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 |
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 |
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 |
CN114931358A (en) * | 2022-04-06 | 2022-08-23 | 上海健康医学院 | Large depth of field micrometer resolution optical coherence tomography endoscopic probe |
CN114931358B (en) * | 2022-04-06 | 2024-07-26 | 上海健康医学院 | Large depth of field micron resolution optical coherence tomography endoscopic probe |
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 | |
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 | |
CN100464696C (en) | Spectral coverage OCT imaging method based on optical scanning delay line and the system | |
CN102188237B (en) | Phase-multiplexing-based full-range sweep frequency OCT (Optical Coherence Tomography) imaging method and system | |
CN101832817B (en) | Parallel complex frequency domain optical coherence tomography method and system | |
CN101297750B (en) | Complex spectral domain optical coherence tomography method and system | |
US9448056B2 (en) | System for fourier domain optical coherence tomography | |
CN102657518B (en) | Method of complex frequency-domain optical coherence tomography using differential sinusoidal phase modulation | |
JP4344829B2 (en) | Polarized light receiving image measuring device | |
CN100520361C (en) | Frequency domain optical coherence tomography method and system for full-depth detection | |
CN102818786B (en) | Sinusoidal phase modulation parallel complex frequency domain optical coherence tomography imaging system and method | |
US7394546B2 (en) | Method and apparatus for full phase interferometry | |
CN104706322B (en) | A kind of frequency sweep optical coherence imaging systems based on optical oomputing | |
CN103070669B (en) | Light-spectrum phase calibration system and method based on cascade Mach-Zehnder interferometer | |
CN101793823A (en) | Apparatus and method for rangings and noise reduction of low coherence interferometry (LCI) and optical coherence tomography (OCT) signals | |
CN2916623Y (en) | Frequency domain optical coherence tomography device for full-depth detection | |
CN105147241A (en) | Method and system based on double-space carrier frequency technology for increasing OCT (optical coherence tomography) imaging depth | |
JP2013064748A (en) | Process, arrangement and system for providing polarization sensitive optical frequency domain imaging of sample | |
CN203016921U (en) | Spectral phase calibration system based on cascaded mach-zehnder interferometer | |
CN104568842B (en) | The device and method of biological tissue's group delay dispersion is measured with reference to FDOCT | |
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 | |
CN103792192B (en) | Based on the polarization domain optical coherence tomography system of simple detector |
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 |