CN111829985A - Bimodal full-field optical coherence tomography device and imaging method - Google Patents

Bimodal full-field optical coherence tomography device and imaging method Download PDF

Info

Publication number
CN111829985A
CN111829985A CN202010649446.2A CN202010649446A CN111829985A CN 111829985 A CN111829985 A CN 111829985A CN 202010649446 A CN202010649446 A CN 202010649446A CN 111829985 A CN111829985 A CN 111829985A
Authority
CN
China
Prior art keywords
sample
coherent
optical coherence
full
coherence tomography
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
CN202010649446.2A
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.)
Hangzhou Jinyu Information Technology Co ltd
Original Assignee
Hangzhou Jinyu Information Technology Co ltd
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 Hangzhou Jinyu Information Technology Co ltd filed Critical Hangzhou Jinyu Information Technology Co ltd
Priority to CN202010649446.2A priority Critical patent/CN111829985A/en
Publication of CN111829985A publication Critical patent/CN111829985A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The utility model provides a bimodulus full field optical coherence tomography device, includes low coherent light source, the output of low coherent light source be connected with beam splitter BS, the one end of beam splitter BS be connected with area array CCD detector through imaging Lens, the other end of beam splitter BS be connected with piezoceramics PZT through objective MO, objective MO is the reference mirror, objective MO set up on the reference arm, the output of reference arm is provided with the sample arm, the arrangement has the sample on the sample arm. The phase difference between the reference arm and the sample arm is kept unchanged, coherent images are continuously acquired and are subjected to standard deviation calculation to obtain coherent images of the weak scattering characteristic structure, and then the coherent images are fused into a full-field optical coherence tomography image, so that submicron-order coherent imaging can be performed on the strong scattering characteristic and weak scattering characteristic structure in the sample in real time, the sample is not required to be subjected to slicing, dyeing, fluorescent labeling and other processing, and the method has a wide clinical application prospect.

Description

Bimodal full-field optical coherence tomography device and imaging method
Technical Field
The invention relates to the technical field of optical coherence tomography, in particular to a bimodal full-field optical coherence tomography device and an imaging method.
Background
Full Field Optical Coherence Tomography (FFOCT) was first proposed in 1998 by Dubois A, Vabre L, Bocara AC et al as a non-invasive and non-destructive micron-scale Optical Tomography technique that allows observation of the microstructure inside biological tissues, and has the characteristics of higher resolution and faster imaging compared to conventional OCT (Optical Coherence Tomography), and has great application prospects in the biomedical Field.
The full-field optical coherence tomography has good optical coherence tomography effect on collagen and fiber with strong scattering property in biological tissues, however, the structure with weak scattering property, such as dynamic change in cells, which provides the same important information in pathological diagnosis is difficult to observe.
Therefore, in order to make up for the defect that the existing full-field optical coherence tomography cannot image the structure with weak scattering characteristics, a bimodal full-field optical coherence tomography device needs to be designed.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a bimodal full-field optical coherence tomography device capable of performing submicron-order coherent imaging on structures with strong scattering characteristics and weak scattering characteristics in a sample in real time and an imaging method using the device.
The technical problem to be solved by the invention is realized by adopting the following technical scheme:
the utility model provides a bimodulus full field optical coherence tomography device, includes low coherent light source, the output of low coherent light source be connected with beam splitter BS, one end of beam splitter BS be connected with area array CCD detector through imaging Lens, the other end of beam splitter BS be connected with piezoceramics PZT through objective MO, objective MO be the reference mirror, objective MO set up on the reference arm, the output of reference arm be provided with the sample arm, the sample arm on settle have the sample.
A bimodal full-field optical coherence tomography method comprises the following steps:
firstly, driving piezoelectric ceramic PZT, shifting the phase of a reference mirror, synchronously acquiring four coherent images by an area array CCD detector, and resolving by utilizing a four-step phase shifting algorithm to obtain a static full-field optical coherence tomography image I at a sample axial depth position DS
Secondly, setting the driving voltage of the piezoelectric ceramic PZT to be zero, keeping the phase difference between the reference arm and the sample arm unchanged, and continuously acquiring N coherent images I by the area array CCD detector within the time Tt1,It2……ItN
Thirdly, calculating and obtaining dynamic change imaging I of the weak scattering structure at the axial depth position D of the sample by using a standard deviation algorithmD
The fourth step, imaging the dynamic change IDMarked green and superimposed on the static full-field optical coherence tomography image ISAnd obtaining the bimodal full-field coherent tomography image at the axial depth position D of the sample.
Preferably, the calculation process of the four-step algorithm in the first step is as follows:
firstly, taking the coherent phase difference phi value as 0,
Figure BDA0002574337470000021
Pi and
Figure BDA0002574337470000022
collecting according to the selected coherent phase difference to obtain four coherent light intensity images respectively I (x, y; 0),
Figure BDA0002574337470000023
I (x, y; pi) and
Figure BDA0002574337470000024
wherein: x is the abscissa of the pixel point of the coherent light intensity image, y is the ordinate of the pixel point of the coherent light intensity image, and then each pixel position of the coherent light intensity image is scanned in sequence according to the positionFormula (II)
Figure BDA0002574337470000025
Calculating to obtain a static full-field optical coherence tomography image IS
Preferably, the calculation formula of the standard deviation algorithm in the third step is
Figure BDA0002574337470000031
Wherein D (x, y) is the dynamic signal strength at the pixel point (x, y),<S(x,y)>for the mean value of the interference signal intensity, S (x, y, t), calculated from a small sample size of the M imagesi) Is the signal strength of the pixel point (x, y) at the N coherent images.
The invention has the advantages and positive effects that:
the phase difference between the reference arm and the sample arm is kept unchanged, coherent images are continuously acquired and are subjected to standard deviation calculation to obtain coherent images of a weak scattering characteristic structure, and then the coherent images are fused into a full-field optical coherence tomography image, so that submicron-order coherent imaging can be performed on the strong scattering characteristic and weak scattering characteristic structures in the sample in real time, the sample is not required to be subjected to slicing, dyeing, fluorescent labeling and the like, frozen sections and paraffin sections with the defects of complex section making, low resolution and the like can be replaced, and the method has a wide clinical application prospect.
Drawings
FIG. 1 is a schematic structural diagram of a bimodal full-field optical coherence tomography apparatus of the present invention;
FIG. 2 is a schematic step diagram of the bimodal full field optical coherence tomography method of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It will be understood that when an element is referred to as being "secured to" another element, it can be directly on the other element or intervening elements may also be present. When a component is referred to as being "connected" to another component, it can be directly connected to the other component or intervening components may also be present. When a component is referred to as being "disposed on" another component, it can be directly on the other component or intervening components may also be present. The terms "vertical," "horizontal," "left," "right," and the like as used herein are for illustrative purposes only.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
The embodiments of the invention will be described in further detail below with reference to the accompanying drawings:
as shown in FIG. 1, the bimodal full-field optical coherence tomography apparatus of the present invention includes a low coherence light source, an output end of the low coherence light source is connected with a beam splitter BS, one end of the beam splitter BS is connected with an area array CCD detector through an imaging Lens, the other end of the beam splitter BS is connected with a piezoelectric ceramic PZT through an objective MO, the objective MO is a reference mirror, the objective MO is disposed on a reference arm, an output end of the reference arm is disposed with a sample arm, and a sample is disposed on the sample arm.
As shown in FIG. 2, the bimodal full-field optical coherence tomography method of the present invention includes the following steps:
firstly, driving piezoelectric ceramic PZT to shift the phase of the reference mirror, and detecting by planar array CCDThe device synchronously acquires four coherent images, and the static full-field optical coherence tomography image I at the axial depth position D of the sample is obtained by resolving through a four-step phase-shifting algorithmSThe calculation process of the four-step moving term algorithm is as follows:
firstly, taking the coherent phase difference phi value as 0,
Figure BDA0002574337470000041
Pi and
Figure BDA0002574337470000042
collecting according to the selected coherent phase difference to obtain four coherent light intensity images respectively I (x, y; 0),
Figure BDA0002574337470000043
I (x, y; pi) and
Figure BDA0002574337470000044
wherein: x is the abscissa of the pixel point of the coherent light intensity image, y is the ordinate of the pixel point of the coherent light intensity image, then each pixel position of the coherent light intensity image is scanned in sequence and according to the formula
Figure BDA0002574337470000045
Calculating to obtain a static full-field optical coherence tomography image IS
Secondly, setting the driving voltage of the piezoelectric ceramic PZT to be zero, keeping the phase difference between the reference arm and the sample arm unchanged, and continuously acquiring N coherent images I by the area array CCD detector within the time Tt1,It2……ItN
Thirdly, calculating and obtaining dynamic change imaging I of the weak scattering structure at the axial depth position D of the sample by using a standard deviation algorithmDThe standard deviation algorithm formula is as follows:
Figure BDA0002574337470000051
wherein D (x, y) is at the pixelThe dynamic signal strength at point (x, y),<S(x,y)>for the mean value of the interference signal intensity, S (x, y, t), calculated from a small sample size of the M imagesi) Signal strength of pixel point (x, y) at N coherent images;
the fourth step, imaging the dynamic change IDMarked green and superimposed on the static full-field optical coherence tomography image ISAnd obtaining the bimodal full-field coherent tomography image at the axial depth position D of the sample.
Furthermore, the objective lens MO adopts the umplan FLN model of Olympus, the I/O driving control card used in the bimodal full-field optical coherence tomography device system adopts an NI6233 board card of NI, and the piezoelectric ceramic PZT adopts the PST150/7/20VS12 model of Piezomechanik.
Furthermore, the area array CCD detector is a silicon-based multi-channel array detector and can detect ultraviolet, visible and near infrared light. The sensor is a high-sensitivity semiconductor device and is suitable for analyzing weak Raman signals, and a CCD detector allows multi-channel operation (the whole spectrum can be detected in one collection), so the sensor is very suitable for detecting the Raman signals. CCD detectors are widely used, with CCD detectors being used as sensors in digital cameras, and higher-level CCD detectors being used in scientific spectrometers for better sensitivity, uniformity and noise characteristics. CCD detectors are typically one-dimensional (linear) or two-dimensional (planar) arrays consisting of thousands of individual detector elements (also called pixels). The more charge each element is exposed to light, the more intense the light is and the longer the time is, the more charge is generated. The final read-out electronics draws charge from the picture elements so that each charge is read out for measurement. In a common raman spectrometer, raman scattering is first dispersed by a diffraction grating and then projected onto a long axis of a CCD array, a first pixel detects a spectral low-wave-number start signal, a second pixel detects a next spectral position signal, and so on, and the last pixel will detect a spectral high-wave-number end signal. The CCD detector needs to be cooled to a lower temperature to acquire a high-quality spectrum, and two cooling modes are generally adopted, wherein one cooling mode is semiconductor refrigeration and can achieve the lowest temperature of-90 ℃, the other cooling mode is liquid nitrogen low-temperature refrigeration, and the lowest temperature can achieve-196 ℃.
Further, piezoelectric ceramics are information functional ceramic materials capable of converting mechanical energy and electrical energy to each other, and have dielectric properties, elasticity, and the like in addition to piezoelectric properties, and have been widely used in medical imaging, acoustic sensors, acoustic transducers, ultrasonic motors, and the like. The piezoelectric ceramic is manufactured by utilizing the piezoelectric effect that the material causes the relative displacement of the centers of positive and negative charges in the material to generate polarization under the action of mechanical stress, so that bound charges with opposite signs appear on the surfaces of two ends of the material, and has the sensitive characteristic, and the piezoelectric ceramic has the piezoelectric characteristic and is mainly different from a typical piezoelectric quartz crystal which does not contain ferroelectric components: in order to make the piezoelectric ceramic show macroscopic piezoelectric property, the piezoelectric ceramic must be fired and then placed under a strong direct current electric field after being repolarized on the end face, so that the respective polarization vectors of the original disordered orientations are preferentially oriented along the direction of the electric field, and the piezoelectric ceramic after being polarized can retain certain macroscopic residual polarization intensity after the electric field is cancelled, thereby leading the ceramic to have certain piezoelectric property.
In specific implementation, a reference mirror is subjected to phase shifting by adopting a piezoelectric ceramic PZT (piezoelectric ceramic) according to a white light low coherence interference principle, an optical path between the reference mirror and a sample mirror is changed, an interference signal is received by using an area array CCD (charge coupled device) detector, a strong scattering characteristic static two-dimensional chromatographic image at a certain depth position of a sample is obtained by computer processing, a phase is kept unchanged by controlling the piezoelectric ceramic PZT, the area array CCD detector continuously collects the interference signal, and a weak scattering characteristic dynamic two-dimensional chromatographic image at the same depth position of the sample is obtained by computer processing.
The phase difference between the reference arm and the sample arm is kept unchanged, coherent images are continuously acquired and are subjected to standard deviation calculation to obtain coherent images of a weak scattering characteristic structure, and then the coherent images are fused into a full-field optical coherence tomography image, so that submicron-order coherent imaging can be performed on the strong scattering characteristic and weak scattering characteristic structures in the sample in real time, the sample is not required to be subjected to slicing, dyeing, fluorescent labeling and the like, frozen sections and paraffin sections with the defects of complex section making, low resolution and the like can be replaced, and the method has a wide clinical application prospect.
It should be emphasized that the embodiments described herein are illustrative rather than restrictive, and thus the present invention is not limited to the embodiments described in the detailed description, but other embodiments derived from the technical solutions of the present invention by those skilled in the art are also within the scope of the present invention.

Claims (4)

1. A kind of bimodal whole field optical coherence tomography imaging device, characterized by that: the device comprises a low-coherence light source, wherein the output end of the low-coherence light source is connected with a beam splitter BS, one end of the beam splitter BS is connected with an area array CCD detector through an imaging Lens, the other end of the beam splitter BS is connected with a piezoelectric ceramic PZT through an objective MO, the objective MO is a reference mirror and is arranged on a reference arm, the output end of the reference arm is provided with a sample arm, and a sample is arranged on the sample arm.
2. A bimodal full-field optical coherence tomography method is characterized in that: the method comprises the following steps:
firstly, driving piezoelectric ceramic PZT, shifting the phase of a reference mirror, synchronously acquiring four coherent images by an area array CCD detector, and resolving by utilizing a four-step phase shifting algorithm to obtain a static full-field optical coherence tomography image I at a sample axial depth position DS
Secondly, setting the driving voltage of the piezoelectric ceramic PZT to be zero, keeping the phase difference between the reference arm and the sample arm unchanged, and continuously acquiring N coherent images I by the area array CCD detector within the time Tt1,It2……ItN
Thirdly, calculating and obtaining the axial depth of the sample by using a standard deviation algorithmDynamic change imaging of weakly scattering structures at degree position D ID
The fourth step, imaging the dynamic change IDMarked green and superimposed on the static full-field optical coherence tomography image ISAnd obtaining the bimodal full-field coherent tomography image at the axial depth position D of the sample.
3. The bi-modal full-field optical coherence tomography method of claim 2, wherein: the calculation process of the four-step algorithm in the first step is as follows:
firstly, taking the coherent phase difference phi value as 0,
Figure FDA0002574337460000011
Pi and
Figure FDA0002574337460000012
collecting according to the selected coherent phase difference to obtain four coherent light intensity images respectively I (x, y; 0),
Figure FDA0002574337460000013
I (x, y; pi) and
Figure FDA0002574337460000014
wherein: x is the abscissa of the pixel point of the coherent light intensity image, y is the ordinate of the pixel point of the coherent light intensity image, then each pixel position of the coherent light intensity image is scanned in sequence and according to the formula
Figure FDA0002574337460000021
Calculating to obtain a static full-field optical coherence tomography image IS
4. The bi-modal full-field optical coherence tomography method of claim 2, wherein: the calculation formula of the standard deviation algorithm in the third step is
Figure FDA0002574337460000022
Wherein D (x, y) is the dynamic signal strength at the pixel point (x, y),<S(x,y)>for the mean value of the interference signal intensity, S (x, y, t), calculated from a small sample size of the M imagesi) Is the signal strength of the pixel point (x, y) at the N coherent images.
CN202010649446.2A 2020-07-08 2020-07-08 Bimodal full-field optical coherence tomography device and imaging method Pending CN111829985A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010649446.2A CN111829985A (en) 2020-07-08 2020-07-08 Bimodal full-field optical coherence tomography device and imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010649446.2A CN111829985A (en) 2020-07-08 2020-07-08 Bimodal full-field optical coherence tomography device and imaging method

Publications (1)

Publication Number Publication Date
CN111829985A true CN111829985A (en) 2020-10-27

Family

ID=72899702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010649446.2A Pending CN111829985A (en) 2020-07-08 2020-07-08 Bimodal full-field optical coherence tomography device and imaging method

Country Status (1)

Country Link
CN (1) CN111829985A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114813695A (en) * 2022-03-08 2022-07-29 华南师范大学 Bimodal imaging system and method based on Raman scattering spectrum and quantitative phase

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150109623A1 (en) * 2012-04-23 2015-04-23 Ben-Gurion University Of The Negev Research And Development Authority True-spectroscopic dual mode high resolution full-field optical coherence tomography using liquid crystal devices
CN106175689A (en) * 2016-07-29 2016-12-07 东北大学秦皇岛分校 Whole-field optically coherent chromatography method and device
CN107452029A (en) * 2017-07-31 2017-12-08 中国医学科学院生物医学工程研究所 A kind of optics microvascular blood flow imaging method
CN107743582A (en) * 2015-04-10 2018-02-27 Ll技术管理公司 Method and system for whole audience interference micro-imaging
CN107788950A (en) * 2017-09-30 2018-03-13 浙江大学 Blood flow imaging method and system based on adaptive threshold fuzziness
CN108245130A (en) * 2016-12-28 2018-07-06 南京理工大学 A kind of optical coherence tomography angiographic apparatus and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150109623A1 (en) * 2012-04-23 2015-04-23 Ben-Gurion University Of The Negev Research And Development Authority True-spectroscopic dual mode high resolution full-field optical coherence tomography using liquid crystal devices
CN107743582A (en) * 2015-04-10 2018-02-27 Ll技术管理公司 Method and system for whole audience interference micro-imaging
CN106175689A (en) * 2016-07-29 2016-12-07 东北大学秦皇岛分校 Whole-field optically coherent chromatography method and device
CN108245130A (en) * 2016-12-28 2018-07-06 南京理工大学 A kind of optical coherence tomography angiographic apparatus and method
CN107452029A (en) * 2017-07-31 2017-12-08 中国医学科学院生物医学工程研究所 A kind of optics microvascular blood flow imaging method
CN107788950A (en) * 2017-09-30 2018-03-13 浙江大学 Blood flow imaging method and system based on adaptive threshold fuzziness

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡斌;江兴方;石晓灏;刘磊;邓林红;: "全场光学相干层析在细胞成像中的应用研究", 量子电子学报, no. 04, pages 437 - 443 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114813695A (en) * 2022-03-08 2022-07-29 华南师范大学 Bimodal imaging system and method based on Raman scattering spectrum and quantitative phase

Similar Documents

Publication Publication Date Title
CN109124615B (en) Selective area high dynamic laser speckle blood flow imaging device and method
Ovcharenko et al. A novel test rig for in situ and real time optical measurement of the contact area evolution during pre-sliding of a spherical contact
US7787112B2 (en) Depth of field extension for optical tomography
US7088455B1 (en) Methods and apparatus for material evaluation using laser speckle
US7330574B2 (en) Best-focus estimation by lateral scanning
JP4317743B2 (en) Mechanical property measurement method based on in-situ determination of optical indenter contact surface and its test equipment
CN106705911A (en) Thermal wave imaging film layer thickness detection system and method
CN111829985A (en) Bimodal full-field optical coherence tomography device and imaging method
CN113390919B (en) Method for observing material phase boundary by phase-locked infrared imaging
CN109085119A (en) A kind of the copolymerization coke 3-D imaging system and implementation method of the detection of Raman tomographic spectroscopy
CN113794839A (en) Focal plane automatic focusing method of imaging ellipsometer
JP3365474B2 (en) Polarizing imaging device
CN111982825B (en) Full-polarization fast dynamic Stokes imaging method
KR102483017B1 (en) Measuring path delay through a liquid-crystal variable retarder at non-uniform retardance intervals
Li et al. High-throughput measurement of coefficient of thermal expansion using a high-resolution digital single-lens reflex camera and digital image correlation
KR101762384B1 (en) Observation system and method for 2-dimensional materials
JP2545209B2 (en) Crystal defect inspection method and inspection apparatus therefor
CN111487190B (en) Single-ion imaging detection method and device
CN108692676A (en) The three-dimensional shape measuring method of sweep type white interference microscope is used
JPS63222206A (en) Image processor for measuring structure of optical fiber
CN117470125B (en) Phase contrast optical coherence elastography strain calculation method and system
Sanderson et al. Quantification of ciliary beat frequency and metachrony by high-speed digital video
RU2649247C1 (en) Method of the active thermal non-destructive control results analysis of articles from polymer composite materials
Zhang et al. Incidence-modularized 3× 4 Mueller matrix polarimetry: A targeted tool with high temporal resolution for in-vivo tissue imaging
CN114216562A (en) Intensity-time integration type rapid polarization imaging method and device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20201027

RJ01 Rejection of invention patent application after publication