[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN113406008B - Photoacoustic imaging device and method based on white light interference - Google Patents

Photoacoustic imaging device and method based on white light interference Download PDF

Info

Publication number
CN113406008B
CN113406008B CN202110683906.8A CN202110683906A CN113406008B CN 113406008 B CN113406008 B CN 113406008B CN 202110683906 A CN202110683906 A CN 202110683906A CN 113406008 B CN113406008 B CN 113406008B
Authority
CN
China
Prior art keywords
sample
photoacoustic
optical fiber
light source
excitation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110683906.8A
Other languages
Chinese (zh)
Other versions
CN113406008A (en
Inventor
王毅
郭振宇
马振鹤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN202110683906.8A priority Critical patent/CN113406008B/en
Publication of CN113406008A publication Critical patent/CN113406008A/en
Application granted granted Critical
Publication of CN113406008B publication Critical patent/CN113406008B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/1702Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
    • 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
    • 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
    • G01N2021/451Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods for determining the optical absorption

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

一种基于白光干涉的光声成像装置及方法,装置包括光学干涉检测系统、光声激发系统、扫描系统及计算机;光学干涉检测系统一端与计算机进行电连接,光学干涉检测系统另一端与扫描系统进行光连接;光声激发系统一端与计算机进行电连接,光声激发系统另一端与扫描系统进行光连接;扫描系统与样品进行光连接。本发明的基于白光干涉的光声成像方法,利用白光干涉检测样品内光声信号激发点的反射光强变化,可以得到同一时刻来自样品不同深度的光声信号,而且具有一定的深度分辨能力,并且采用傅里叶变化进行信号解调,相对于目前所采用的时域解调方式,可有效提高探测灵敏度。

Figure 202110683906

A photoacoustic imaging device and method based on white light interference, the device comprises an optical interference detection system, a photoacoustic excitation system, a scanning system and a computer; one end of the optical interference detection system is electrically connected with the computer, and the other end of the optical interference detection system is connected with the scanning system Optical connection is performed; one end of the photoacoustic excitation system is electrically connected to the computer, and the other end of the photoacoustic excitation system is optically connected to the scanning system; the scanning system is optically connected to the sample. The white light interference-based photoacoustic imaging method of the present invention uses the white light interference to detect the reflected light intensity change of the excitation point of the photoacoustic signal in the sample, and can obtain the photoacoustic signals from different depths of the sample at the same time, and has a certain depth resolution capability, And using Fourier transform to demodulate the signal, compared with the time domain demodulation method currently used, can effectively improve the detection sensitivity.

Figure 202110683906

Description

一种基于白光干涉的光声成像装置及方法A photoacoustic imaging device and method based on white light interference

技术领域technical field

本发明属于光声成像技术领域,特别是涉及一种基于白光干涉的光声成像装置及方法。The invention belongs to the technical field of photoacoustic imaging, and in particular relates to a photoacoustic imaging device and method based on white light interference.

背景技术Background technique

光声成像(Photoacoustic imaging,PAI)是近几年迅速发展的医学成像技术,PAI结合了纯光学成像的高对比度和纯超声成像的深穿透性,是一种非侵入性成像模式,可用于生物组织的结构、功能成像。PAI的原理是光声效应,当脉冲激光照射到生物组织时,组织吸收光能并产生热弹性膨胀,并由此产生相应的超声波,检测超声波可得到组织的吸收分布图像。Photoacoustic imaging (PAI) is a rapidly developing medical imaging technology in recent years. PAI combines the high contrast of pure optical imaging and the deep penetration of pure ultrasound imaging. It is a non-invasive imaging mode that can be used for Structural and functional imaging of biological tissues. The principle of PAI is the photoacoustic effect. When the pulsed laser is irradiated to the biological tissue, the tissue absorbs the light energy and generates thermoelastic expansion, and thus generates the corresponding ultrasonic wave. The absorption distribution image of the tissue can be obtained by detecting the ultrasonic wave.

目前,基于压电换能器的接触式光声成像技术作为一种比较成熟的技术而广泛使用,该技术使用压电换能器直接检测光声信号,但是由于超声波在不同介质中的声阻抗不同,这也使得超声在两种介质面产生强烈的反射,所以为了提高灵敏度及减小损耗,声耦合介质是必要的,这也从原理上限制了该技术的使用范围。At present, contact photoacoustic imaging technology based on piezoelectric transducers is widely used as a relatively mature technology. This technology uses piezoelectric transducers to directly detect photoacoustic signals, but due to the acoustic impedance of ultrasonic waves in different media Differently, this also makes the ultrasound produce strong reflection on the two medium surfaces, so in order to improve the sensitivity and reduce the loss, an acoustic coupling medium is necessary, which also limits the scope of application of this technology in principle.

对于烧伤、脑科等特殊情况的检查,就需要一种非接触的检测方式,因此基于光学干涉的非接触检测方法作为压电换能器的一种优化方法被提出。该方法可以实现非接触并获得样品信息,与压电换能器相比,光学干涉检测具有非接触、小型化、高灵敏度等特点。For the inspection of special cases such as burns and brain diseases, a non-contact detection method is required. Therefore, a non-contact detection method based on optical interference is proposed as an optimization method of piezoelectric transducers. This method can realize non-contact and obtain sample information. Compared with piezoelectric transducers, optical interference detection has the characteristics of non-contact, miniaturization and high sensitivity.

例如,申请号为201510881786.7的中国专利申请,其公开了一种基于光学干涉法的非接触光声探测方法及装置,并提出了一种基于光学干涉的检测方法,但是该光学干涉的光声成像方法仍存在缺陷。为了解决由于组织样品表面粗糙而导致的反射的探测光的强度及相位较弱并出现随机变化,该方法在样品表面涂抹水层,水层产生一个均匀的反射面,但是这种在样品表面添加水层的方式,实质并没有做到完全非接触,在应用上同样存在不便之处。For example, the Chinese patent application with the application number of 201510881786.7 discloses a non-contact photoacoustic detection method and device based on optical interferometry, and proposes a detection method based on optical interference, but the photoacoustic imaging of the optical interference The method still has flaws. In order to solve the weak and random variation in the intensity and phase of the reflected probe light caused by the rough surface of the tissue sample, this method applies a water layer on the surface of the sample, and the water layer produces a uniform reflection surface, but this method adds a layer of water on the surface of the sample. The method of the water layer is not completely non-contact in essence, and there are also inconveniences in application.

此外,申请号为201910587193.8中国专利申请,其公开了一种非接触光声成像装置及方法,并提出了一种基于3×3光纤耦合器解调的非接触光声成像系统和方法,利用光学干涉方法直接检测样品内光声信号激发点的反射光强变化,并使激发光和样品光的焦点重合于样品内部,焦点处的吸收体吸收激光能量导致该位置的光学折射率变化,进而引起背向散射的样品光光强增大,再利用3×3光纤耦合器解调的干涉方法测量这种光强变化进行成像。虽然该方法解决了水层问题,也提升了方法的灵活性,但是该方法依然存在缺陷,即在同一时刻对于样品内部的深度信息是模糊的,也就无法得到样品中光声信号激发位置的深度信息,并且解调结果过度依赖3×3光纤耦合器的制造精度,否则会产生较大误差。In addition, the application number is 201910587193.8 Chinese patent application, which discloses a non-contact photoacoustic imaging device and method, and proposes a non-contact photoacoustic imaging system and method based on 3×3 fiber coupler demodulation, using optical The interference method directly detects the change of the reflected light intensity of the excitation point of the photoacoustic signal in the sample, and makes the focus of the excitation light and the sample light coincide with the inside of the sample. The light intensity of the backscattered sample light increases, and the interferometric method of demodulation by a 3×3 fiber coupler is used to measure this light intensity change for imaging. Although this method solves the problem of water layer and improves the flexibility of the method, it still has defects, that is, the depth information inside the sample is ambiguous at the same time, so it is impossible to obtain the excitation position of the photoacoustic signal in the sample. depth information, and the demodulation result is overly dependent on the manufacturing accuracy of the 3×3 fiber coupler, otherwise there will be large errors.

发明内容SUMMARY OF THE INVENTION

针对现有技术存在的问题,本发明提供一种基于白光干涉的光声成像装置及方法,利用白光干涉检测样品内光声信号激发点的反射光强变化,可以得到同一时刻来自样品不同深度的光声信号,而且具有一定的深度分辨能力,并且采用傅里叶变化进行信号解调,相对于目前所采用的时域解调方式,可有效提高探测灵敏度。In view of the problems existing in the prior art, the present invention provides a photoacoustic imaging device and method based on white light interference, which utilizes white light interference to detect changes in the reflected light intensity of the excitation point of the photoacoustic signal in the sample, and can obtain images from different depths of the sample at the same time. The photoacoustic signal has a certain depth resolution capability, and the Fourier transform is used for signal demodulation, which can effectively improve the detection sensitivity compared with the currently used time domain demodulation method.

为了实现上述目的,本发明采用如下技术方案:一种基于白光干涉的光声成像装置,包括光学干涉检测系统、光声激发系统、扫描系统及计算机;所述光学干涉检测系统一端与计算机进行电连接,光学干涉检测系统另一端与扫描系统进行光连接;所述光声激发系统一端与计算机进行电连接,光声激发系统另一端与扫描系统进行光连接;所述扫描系统与样品进行光连接。In order to achieve the above purpose, the present invention adopts the following technical scheme: a photoacoustic imaging device based on white light interference, including an optical interference detection system, a photoacoustic excitation system, a scanning system and a computer; one end of the optical interference detection system is electrically connected to the computer. The other end of the optical interference detection system is optically connected to the scanning system; one end of the photoacoustic excitation system is electrically connected to the computer, and the other end of the photoacoustic excitation system is optically connected to the scanning system; the scanning system is optically connected to the sample .

所述光学干涉检测系统包括低相干光源、光纤隔离器、光开关、2×2光纤耦合器、第一准直器、第一透镜、第一反射镜及光谱仪;所述低相干光源通过光纤隔离器与光开关进行光连接,光开关一路与计算机进行电连接,光开关另一路与2×2光纤耦合器通过光纤进行连接;所述第一准直器一端与2×2光纤耦合器通过光纤进行连接,第一准直器另一端发出的激光通过第一透镜射向第一反射镜;所述光谱仪一端与计算机进行电连接,光谱仪另一端与2×2光纤耦合器通过光纤进行连接。The optical interference detection system includes a low-coherence light source, an optical fiber isolator, an optical switch, a 2×2 optical fiber coupler, a first collimator, a first lens, a first mirror and a spectrometer; the low-coherence light source is isolated by an optical fiber The first collimator is optically connected to the optical switch, one way of the optical switch is electrically connected to the computer, and the other way of the optical switch is connected to the 2×2 fiber optic coupler through an optical fiber; one end of the first collimator is connected to the 2×2 fiber optic coupler through an optical fiber For connection, the laser light emitted by the other end of the first collimator is directed to the first mirror through the first lens; one end of the spectrometer is electrically connected to the computer, and the other end of the spectrometer is connected to the 2×2 fiber coupler through optical fibers.

所述光声激发系统包括激发光源及第二反射镜;所述激发光源与计算机进行电连接,激发光源发出的激发光直接射向第二反射镜。The photoacoustic excitation system includes an excitation light source and a second reflection mirror; the excitation light source is electrically connected to a computer, and the excitation light emitted by the excitation light source is directed toward the second reflection mirror.

所述扫描系统包括第二准直器、二向色镜、Y向扫描振镜、X向扫描振镜及第二透镜;所述第二准直器一端与2×2光纤耦合器通过光纤进行连接,第二准直器另一端发出的激光依次通过二向色镜、Y向扫描振镜、X向扫描振镜及第二透镜射向样品;由所述第二反射镜反射的激光直接射向二向色镜,并依次通过Y向扫描振镜、X向扫描振镜及第二透镜射向样品。The scanning system includes a second collimator, a dichroic mirror, a Y-direction scanning galvanometer, an X-direction scanning galvanometer and a second lens; one end of the second collimator is connected to a 2×2 fiber coupler through an optical fiber. connected, the laser light emitted from the other end of the second collimator is directed to the sample through the dichroic mirror, the Y-direction scanning galvanometer, the X-direction scanning galvanometer and the second lens in turn; the laser reflected by the second mirror directly irradiates the sample The dichroic mirror is directed to the sample through the Y-direction scanning galvanometer, the X-direction scanning galvanometer and the second lens in turn.

一种基于白光干涉的光声成像方法,采用了所述的基于白光干涉的光声成像装置,包括如下步骤:A photoacoustic imaging method based on white light interference adopts the photoacoustic imaging device based on white light interference, comprising the following steps:

步骤一:启动低相干光源,由低相干光源发出的探测激光依次经光纤隔离器及光开关进入2×2光纤耦合器,再通过2×2光纤耦合器分两路输出,一路作为参考光依次通过第一准直器、第一透镜及第一反射镜后并原路返回2×2光纤耦合器,另一路作为样品光依次经过第二准直器、二向色镜、Y向扫描振镜、X向扫描振镜、第二透镜及样品后并原路返回2×2光纤耦合器;Step 1: Start the low-coherence light source, and the detection laser emitted by the low-coherence light source enters the 2×2 fiber coupler through the fiber isolator and the optical switch in turn, and then outputs two paths through the 2×2 fiber coupler, and one path is used as the reference light in turn. After passing through the first collimator, the first lens and the first reflector, and returning to the 2×2 fiber coupler in the same way, the other way as the sample light passes through the second collimator, the dichroic mirror and the Y-direction scanning galvanometer in turn. , X-direction scanning galvanometer, second lens and sample and return to 2×2 fiber coupler in the same way;

步骤二:当参考光和样品光原路返回2×2光纤耦合器后,会通过2×2光纤耦合器直接进入光谱仪,之后由计算机进行光谱分析;Step 2: When the reference light and the sample light return to the 2×2 fiber coupler, they will directly enter the spectrometer through the 2×2 fiber coupler, and then the computer will perform spectral analysis;

步骤三:由计算机向激发光源发出触发信号,促使激发光源启动,当激发光触发信号为高电平时,由激发光源输出激发光,激发光依次经过第二反射镜、二向色镜、Y向扫描振镜、X向扫描振镜及第二透镜射向样品,并与样品光在样品内部汇聚在一点上,样品吸收激光能量后会产生光声压,光声压又会促使样品内部激发点的光学折射率发生变大,进而使反射光强增大,并产生光声信号;反之,当激发光触发信号为低电平时,激发光源暂停输出激发光;Step 3: The computer sends a trigger signal to the excitation light source to prompt the excitation light source to start. When the excitation light trigger signal is at a high level, the excitation light is output from the excitation light source, and the excitation light passes through the second reflecting mirror, the dichroic mirror, and the Y direction in turn. The scanning galvanometer, X-direction scanning galvanometer and the second lens are directed to the sample, and converge with the sample light at a point inside the sample. After the sample absorbs the laser energy, a photoacoustic pressure will be generated, and the photoacoustic pressure will stimulate the excitation point inside the sample. The optical index of refraction increases, which in turn increases the reflected light intensity and generates a photoacoustic signal; on the contrary, when the excitation light trigger signal is at a low level, the excitation light source stops outputting the excitation light;

步骤四:由计算机向光开关发出触发信号,当光开关触发信号为高电平时,由低相干光源发出的探测激光可通过光开关直达样品;反之,当光开关触发信号为低电平时,由低相干光源发出的探测激光不可通过光开关;Step 4: Send a trigger signal to the optical switch by the computer. When the trigger signal of the optical switch is at a high level, the detection laser emitted by the low-coherence light source can directly reach the sample through the optical switch; on the contrary, when the trigger signal of the optical switch is at a low level The detection laser emitted by the low-coherence light source cannot pass through the optical switch;

步骤五:当激发光触发信号处于低电平周期且光开关触发信号处于高电平周期时,将检测到的光声信号记为S0(K);此外,当激发光触发信号和光开关触发信号均处于高电平周期时,将检测到的光声信号记为S1(K);其中,K表示光谱仪的波数坐标;Step 5: When the excitation light trigger signal is in the low level period and the optical switch trigger signal is in the high level period, the detected photoacoustic signal is recorded as S0(K); in addition, when the excitation light trigger signal and the optical switch trigger signal are When both are in the high-level period, the detected photoacoustic signal is recorded as S1(K); wherein, K represents the wavenumber coordinate of the spectrometer;

步骤六:对光声信号S0(K)和光声信号S1(K)进行预处理,消除直流分量及强度归一化,然后进行快速傅里叶变换,可得到光声信号S0(K)的幅度谱F0(u)以及光声信号S1(K)的幅度谱F1(u),其中u表示频率,频率u和深度z成正比,且z=au,式中,a为比例系数,且比例系数a为已知量,可通过测量一个已知深度的样品进行确定;最后,通过公式z=au可以得到不同深度的样品反射光强分布F0(z)和F1(z),其中,F0(z)为激发光源暂停输出激发光时样品的反射光强,F1(z)为激发光源输出激发光时样品的反射光强;Step 6: Preprocess the photoacoustic signal S0(K) and the photoacoustic signal S1(K), eliminate the DC component and normalize the intensity, and then perform fast Fourier transform to obtain the amplitude of the photoacoustic signal S0(K) The spectrum F0(u) and the amplitude spectrum F1(u) of the photoacoustic signal S1(K), where u represents the frequency, the frequency u is proportional to the depth z, and z=au, where a is the proportionality coefficient, and the proportionality coefficient a is a known quantity, which can be determined by measuring a sample with a known depth; finally, the reflected light intensity distributions F0(z) and F1(z) of samples at different depths can be obtained by the formula z=au, where F0(z ) is the reflected light intensity of the sample when the excitation light source stops outputting the excitation light, and F1(z) is the reflected light intensity of the sample when the excitation light source outputs the excitation light;

步骤七:计算样品不同深度的光声信号P(z),计算公式为P(z)=F1(z)-F0(z);Step 7: Calculate the photoacoustic signal P(z) at different depths of the sample, the calculation formula is P(z)=F1(z)-F0(z);

步骤八:通过X向扫描振镜和Y向扫描振镜实现二维扫描,在计算机中实现二维成像。Step 8: realize two-dimensional scanning through the X-direction scanning galvanometer and the Y-direction scanning galvanometer, and realize two-dimensional imaging in the computer.

本发明的有益效果:Beneficial effects of the present invention:

本发明的基于白光干涉的光声成像装置及方法,利用白光干涉检测样品内光声信号激发点的反射光强变化,可以得到同一时刻来自样品不同深度的光声信号,而且具有一定的深度分辨能力,并且采用傅里叶变化进行信号解调,相对于目前所采用的时域解调方式,可有效提高探测灵敏度。The white light interference-based photoacoustic imaging device and method of the present invention utilizes the white light interference to detect the reflected light intensity change of the excitation point of the photoacoustic signal in the sample, and can obtain the photoacoustic signals from different depths of the sample at the same time, and has a certain depth resolution. Compared with the current time domain demodulation method, the detection sensitivity can be effectively improved.

附图说明Description of drawings

图1为本发明的基于白光干涉的光声成像装置的结构原理图;1 is a schematic structural diagram of a photoacoustic imaging device based on white light interference of the present invention;

图2为实施例中的激发光触发信号和光开关触发信号的波形图;2 is a waveform diagram of an excitation light trigger signal and an optical switch trigger signal in an embodiment;

图中,I—光学干涉检测系统,II—光声激发系统,III—扫描系统,1—计算机,2—样品,3—低相干光源,4—光纤隔离器,5—光开关,6—2×2光纤耦合器,7—第一准直器,8—第一透镜,9—第一反射镜,10—光谱仪,11—激发光源,12—第二反射镜,13—第二准直器,14—二向色镜,15—Y向扫描振镜,16—X向扫描振镜,17—第二透镜。In the figure, I—optical interference detection system, II—photoacoustic excitation system, III—scanning system, 1—computer, 2—sample, 3—low coherence light source, 4—fiber isolator, 5—optical switch, 6—2 ×2 fiber coupler, 7—first collimator, 8—first lens, 9—first reflector, 10—spectroscope, 11—excitation light source, 12—second reflector, 13—second collimator , 14—dichroic mirror, 15—Y scanning galvanometer, 16—X scanning galvanometer, 17—second lens.

具体实施方式Detailed ways

下面结合附图和具体实施例对本发明做进一步的详细说明。The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

如图1所示,一种基于白光干涉的光声成像装置,包括光学干涉检测系统I、光声激发系统II、扫描系统III及计算机1;所述光学干涉检测系统I一端与计算机1进行电连接,光学干涉检测系统I另一端与扫描系统III进行光连接;所述光声激发系统II一端与计算机1进行电连接,光声激发系统II另一端与扫描系统III进行光连接;所述扫描系统III与样品2进行光连接。As shown in FIG. 1, a photoacoustic imaging device based on white light interference includes an optical interference detection system I, a photoacoustic excitation system II, a scanning system III and a computer 1; one end of the optical interference detection system I is electrically connected to the computer 1. connected, the other end of the optical interference detection system I is optically connected with the scanning system III; one end of the photoacoustic excitation system II is electrically connected with the computer 1, and the other end of the photoacoustic excitation system II is optically connected with the scanning system III; the scanning System III was optically connected to Sample 2.

所述光学干涉检测系统I包括低相干光源3、光纤隔离器4、光开关5、2×2光纤耦合器6、第一准直器7、第一透镜8、第一反射镜9及光谱仪10;所述低相干光源3通过光纤隔离器4与光开关5进行光连接,光开关5一路与计算机1进行电连接,光开关5另一路与2×2光纤耦合器6通过光纤进行连接;所述第一准直器7一端与2×2光纤耦合器6通过光纤进行连接,第一准直器7另一端发出的激光通过第一透镜8射向第一反射镜9;所述光谱仪10一端与计算机1进行电连接,光谱仪10另一端与2×2光纤耦合器6通过光纤进行连接。The optical interference detection system 1 includes a low-coherence light source 3, an optical fiber isolator 4, an optical switch 5, a 2×2 optical fiber coupler 6, a first collimator 7, a first lens 8, a first reflector 9 and a spectrometer 10 ; The low-coherence light source 3 is optically connected to the optical switch 5 through the optical fiber isolator 4, the optical switch 5 is electrically connected to the computer 1, and the other optical switch 5 is connected to the 2 × 2 optical fiber coupler 6 through optical fibers; One end of the first collimator 7 is connected to the 2×2 fiber coupler 6 through optical fibers, and the laser light emitted by the other end of the first collimator 7 is directed to the first mirror 9 through the first lens 8; one end of the spectrometer 10 is It is electrically connected to the computer 1 , and the other end of the spectrometer 10 is connected to the 2×2 fiber coupler 6 through an optical fiber.

所述光声激发系统II包括激发光源11及第二反射镜12;所述激发光源11与计算机1进行电连接,激发光源11发出的激发光直接射向第二反射镜12。The photoacoustic excitation system II includes an excitation light source 11 and a second reflection mirror 12 ; the excitation light source 11 is electrically connected to the computer 1 , and the excitation light emitted by the excitation light source 11 is directed toward the second reflection mirror 12 .

所述扫描系统III包括第二准直器13、二向色镜14、Y向扫描振镜15、X向扫描振镜16及第二透镜17;所述第二准直器13一端与2×2光纤耦合器6通过光纤进行连接,第二准直器13另一端发出的激光依次通过二向色镜14、Y向扫描振镜15、X向扫描振镜16及第二透镜17射向样品2;由所述第二反射镜12反射的激光直接射向二向色镜14,并依次通过Y向扫描振镜15、X向扫描振镜16及第二透镜17射向样品2。The scanning system III includes a second collimator 13 , a dichroic mirror 14 , a Y-direction scanning galvanometer 15 , an X-direction scanning galvanometer 16 and a second lens 17 ; one end of the second collimator 13 is connected to 2× 2. The optical fiber coupler 6 is connected by an optical fiber, and the laser light emitted from the other end of the second collimator 13 is directed to the sample through the dichroic mirror 14, the Y-direction scanning galvanometer 15, the X-direction scanning galvanometer 16 and the second lens 17 in turn. 2. The laser light reflected by the second reflecting mirror 12 is directed toward the dichroic mirror 14, and then directed toward the sample 2 through the Y-direction scanning galvanometer 15, the X-direction scanning galvanometer 16 and the second lens 17 in sequence.

一种基于白光干涉的光声成像方法,采用了所述的基于白光干涉的光声成像装置,包括如下步骤:A photoacoustic imaging method based on white light interference adopts the photoacoustic imaging device based on white light interference, comprising the following steps:

步骤一:启动低相干光源3,由低相干光源3发出的探测激光依次经光纤隔离器4及光开关5进入2×2光纤耦合器6,再通过2×2光纤耦合器6分两路输出,一路作为参考光依次通过第一准直器7、第一透镜8及第一反射镜9后并原路返回2×2光纤耦合器6,另一路作为样品光依次经过第二准直器13、二向色镜14、Y向扫描振镜15、X向扫描振镜16、第二透镜17及样品2后并原路返回2×2光纤耦合器6;Step 1: Start the low-coherence light source 3, and the detection laser emitted by the low-coherence light source 3 enters the 2×2 fiber coupler 6 through the fiber isolator 4 and the optical switch 5 in turn, and then outputs two channels through the 2×2 fiber coupler 6 , one way as the reference light passes through the first collimator 7, the first lens 8 and the first reflector 9 in turn and returns to the 2×2 fiber coupler 6 in the same way, and the other way as the sample light passes through the second collimator 13 in sequence , the dichroic mirror 14, the Y-direction scanning galvanometer 15, the X-direction scanning galvanometer 16, the second lens 17 and the sample 2 and return to the 2×2 fiber coupler 6 in the same way;

步骤二:当参考光和样品光原路返回2×2光纤耦合器6后,会通过2×2光纤耦合器6直接进入光谱仪10,之后由计算机1进行光谱分析;Step 2: After the reference light and the sample light return to the 2×2 fiber coupler 6, they will directly enter the spectrometer 10 through the 2×2 fiber coupler 6, and then the computer 1 will perform spectral analysis;

步骤三:由计算机1向激发光源11发出触发信号,促使激发光源11启动,当激发光触发信号为高电平时,由激发光源11输出激发光,激发光依次经过第二反射镜12、二向色镜14、Y向扫描振镜15、X向扫描振镜16及第二透镜17射向样品2,并与样品光在样品2内部汇聚在一点上,样品2吸收激光能量后会产生光声压,光声压又会促使样品2内部激发点的光学折射率发生变大,进而使反射光强增大,并产生光声信号;反之,当激发光触发信号为低电平时,激发光源11暂停输出激发光;Step 3: The computer 1 sends a trigger signal to the excitation light source 11 to prompt the excitation light source 11 to start up. When the excitation light trigger signal is at a high level, the excitation light source 11 outputs the excitation light, and the excitation light passes through the second reflecting mirror 12 and the two-way in turn. The chromatic mirror 14, the Y-direction scanning galvanometer 15, the X-direction scanning galvanometer 16 and the second lens 17 are directed towards the sample 2, and converge with the sample light at a point inside the sample 2, and the sample 2 will generate photoacoustics after absorbing the laser energy The photoacoustic pressure will increase the optical refractive index of the excitation point inside the sample 2, thereby increasing the reflected light intensity and generating a photoacoustic signal; on the contrary, when the excitation light trigger signal is low, the excitation light source 11 Pause the output excitation light;

步骤四:由计算机1向光开关5发出触发信号,当光开关触发信号为高电平时,由低相干光源3发出的探测激光可通过光开关5直达样品2;反之,当光开关触发信号为低电平时,由低相干光源3发出的探测激光不可通过光开关5;Step 4: The computer 1 sends a trigger signal to the optical switch 5. When the optical switch trigger signal is at a high level, the detection laser emitted by the low-coherence light source 3 can reach the sample 2 directly through the optical switch 5; on the contrary, when the optical switch trigger signal is When the level is low, the detection laser emitted by the low-coherence light source 3 cannot pass through the optical switch 5;

步骤五:当激发光触发信号处于低电平周期且光开关触发信号处于高电平周期时,将检测到的光声信号记为S0(K);此外,当激发光触发信号和光开关触发信号均处于高电平周期时,将检测到的光声信号记为S1(K);其中,K表示光谱仪10的波数坐标;本实施例中,如图2所示,光开关触发信号的高电平周期与低电平周期的时长比为1:1且交替出现,而激发光触发信号的高电平周期与低电平周期的时长比为1:3且交替出现;Step 5: When the excitation light trigger signal is in the low level period and the optical switch trigger signal is in the high level period, the detected photoacoustic signal is recorded as S0(K); in addition, when the excitation light trigger signal and the optical switch trigger signal are When both are in the high-level period, the detected photoacoustic signal is recorded as S1(K); wherein, K represents the wave number coordinate of the spectrometer 10; in this embodiment, as shown in FIG. The time-length ratio of the flat period to the low-level period is 1:1 and occurs alternately, while the time-length ratio of the high-level period to the low-level period of the excitation light trigger signal is 1:3 and alternately occurs;

步骤六:对光声信号S0(K)和光声信号S1(K)进行预处理,消除直流分量及强度归一化,然后进行快速傅里叶变换,可得到光声信号S0(K)的幅度谱F0(u)以及光声信号S1(K)的幅度谱F1(u),其中u表示频率,频率u和深度z成正比,且z=au,式中,a为比例系数,且比例系数a为已知量,可通过测量一个已知深度的样品2进行确定;最后,通过公式z=au可以得到不同深度的样品2反射光强分布F0(z)和F1(z),其中,F0(z)为激发光源11暂停输出激发光时样品2的反射光强,F1(z)为激发光源11输出激发光时样品2的反射光强;Step 6: Preprocess the photoacoustic signal S0(K) and the photoacoustic signal S1(K), eliminate the DC component and normalize the intensity, and then perform fast Fourier transform to obtain the amplitude of the photoacoustic signal S0(K) The spectrum F0(u) and the amplitude spectrum F1(u) of the photoacoustic signal S1(K), where u represents the frequency, the frequency u is proportional to the depth z, and z=au, where a is the scale factor, and the scale factor a is a known quantity, which can be determined by measuring a sample 2 with a known depth; finally, through the formula z=au, the reflected light intensity distributions F0(z) and F1(z) of sample 2 at different depths can be obtained, where F0 (z) is the reflected light intensity of the sample 2 when the excitation light source 11 stops outputting the excitation light, and F1(z) is the reflected light intensity of the sample 2 when the excitation light source 11 outputs the excitation light;

步骤七:计算样品2不同深度的光声信号P(z),计算公式为P(z)=F1(z)-F0(z);Step 7: Calculate the photoacoustic signal P(z) of sample 2 at different depths, and the calculation formula is P(z)=F1(z)-F0(z);

步骤八:通过X向扫描振镜16和Y向扫描振镜15实现二维扫描,在计算机1中实现二维成像。Step 8: Two-dimensional scanning is realized by the X-direction scanning galvanometer 16 and the Y-direction scanning galvanometer 15 , and two-dimensional imaging is realized in the computer 1 .

实施例中的方案并非用以限制本发明的专利保护范围,凡未脱离本发明所为的等效实施或变更,均包含于本案的专利范围中。The solutions in the embodiments are not intended to limit the scope of the patent protection of the present invention, and all equivalent implementations or modifications that do not depart from the present invention are included in the scope of the patent of this case.

Claims (1)

1. A photoacoustic imaging method based on white light interference adopts a photoacoustic imaging device based on white light interference, and the device comprises an optical interference detection system, a photoacoustic excitation system, a scanning system and a computer; one end of the optical interference detection system is electrically connected with the computer, and the other end of the optical interference detection system is optically connected with the scanning system; one end of the photoacoustic excitation system is electrically connected with the computer, and the other end of the photoacoustic excitation system is optically connected with the scanning system; the scanning system is optically connected with the sample; the optical interference detection system comprises a low-coherence light source, an optical fiber isolator, an optical switch, a 2 x 2 optical fiber coupler, a first collimator, a first lens, a first reflector and a spectrometer; the low-coherence light source is optically connected with the optical switch through the optical fiber isolator, one path of the optical switch is electrically connected with the computer, and the other path of the optical switch is connected with the 2 x 2 optical fiber coupler through the optical fiber; one end of the first collimator is connected with the 2 x 2 optical fiber coupler through an optical fiber, and laser emitted by the other end of the first collimator is emitted to the first reflector through the first lens; one end of the spectrometer is electrically connected with the computer, and the other end of the spectrometer is connected with the 2 x 2 optical fiber coupler through an optical fiber; the photoacoustic excitation system comprises an excitation light source and a second reflecting mirror; the excitation light source is electrically connected with the computer, and the excitation light emitted by the excitation light source directly irradiates the second reflecting mirror; the scanning system comprises a second collimator, a dichroic mirror, a Y-direction scanning galvanometer, an X-direction scanning galvanometer and a second lens; one end of the second collimator is connected with the 2X 2 optical fiber coupler through an optical fiber, and laser emitted by the other end of the second collimator sequentially passes through the dichroic mirror, the Y-direction scanning galvanometer, the X-direction scanning galvanometer and the second lens and is emitted to the sample; the laser reflected by the second reflecting mirror directly irradiates to the dichroic mirror and irradiates to a sample through the Y-direction scanning vibrating mirror, the X-direction scanning vibrating mirror and the second lens in sequence; the method is characterized in that: the method comprises the following steps:
the method comprises the following steps: starting a low-coherence light source, enabling detection laser emitted by the low-coherence light source to enter a 2X 2 optical fiber coupler through an optical fiber isolator and an optical switch in sequence, outputting the detection laser in two paths through the 2X 2 optical fiber coupler, enabling one path of the detection laser to serve as reference light to pass through a first collimator, a first lens and a first reflector in sequence and return to the 2X 2 optical fiber coupler in the original path, enabling the other path of the detection laser to serve as sample light to pass through a second collimator, a dichroic mirror, a Y-direction scanning vibration mirror, an X-direction scanning vibration mirror, a second lens and a sample in sequence and return to the 2X 2 optical fiber coupler in the original path;
step two: when the reference light and the sample light original path return to the 2 x 2 optical fiber coupler, the reference light and the sample light directly enter the spectrometer through the 2 x 2 optical fiber coupler, and then the spectral analysis is carried out by the computer;
step three: the computer sends a trigger signal to the excitation light source to prompt the excitation light source to start, when the trigger signal of the excitation light is at a high level, the excitation light source outputs the excitation light, the excitation light sequentially passes through the second reflecting mirror, the dichroic mirror, the Y-direction scanning vibrating mirror, the X-direction scanning vibrating mirror and the second lens to emit to the sample, and is converged at one point with the sample light inside the sample, the sample can generate light sound pressure after absorbing laser energy, and the light sound pressure can prompt the optical refractive index of an excitation point inside the sample to be increased, so that the reflected light intensity is increased, and a light sound signal is generated; on the contrary, when the exciting light trigger signal is at a low level, the exciting light source stops outputting exciting light;
step four: sending a trigger signal to the optical switch by the computer, and when the trigger signal of the optical switch is at a high level, directly sending detection laser sent by the low-coherence light source to the sample through the optical switch; on the contrary, when the trigger signal of the optical switch is at a low level, the detection laser emitted by the low-coherence light source cannot pass through the optical switch;
step five: when the excitation light trigger signal is in a low level period and the optical switch trigger signal is in a high level period, recording the detected photoacoustic signal as S0 (K); further, when both the excitation light trigger signal and the optical switch trigger signal are in the high level period, the detected photoacoustic signal is denoted as S1 (K); wherein K represents the wave number coordinate of the spectrometer;
step six: preprocessing a photoacoustic signal S0(K) and a photoacoustic signal S1(K), eliminating a direct current component and normalizing the intensity, and then performing fast Fourier transform to obtain an amplitude spectrum F0(u) of the photoacoustic signal S0(K) and an amplitude spectrum F1(u) of the photoacoustic signal S1(K), wherein u represents frequency, the frequency u is in direct proportion to a depth z, and the z is au, wherein a is a proportionality coefficient, and a is a known quantity, which is determined by measuring a sample with a known depth; finally, the reflected light intensity distributions F0(z) and F1(z) of the samples at different depths are obtained by the formula z au, wherein F0(z) is the reflected light intensity of the sample when the excitation light source stops outputting the excitation light, and F1(z) is the reflected light intensity of the sample when the excitation light source outputs the excitation light;
step seven: calculating photoacoustic signals P (z) of different depths of the sample, wherein the calculation formula is P (z) ═ F1(z) -F0 (z);
step eight: two-dimensional scanning is realized through the X-direction scanning galvanometer and the Y-direction scanning galvanometer, and two-dimensional imaging is realized in a computer.
CN202110683906.8A 2021-06-21 2021-06-21 Photoacoustic imaging device and method based on white light interference Active CN113406008B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110683906.8A CN113406008B (en) 2021-06-21 2021-06-21 Photoacoustic imaging device and method based on white light interference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110683906.8A CN113406008B (en) 2021-06-21 2021-06-21 Photoacoustic imaging device and method based on white light interference

Publications (2)

Publication Number Publication Date
CN113406008A CN113406008A (en) 2021-09-17
CN113406008B true CN113406008B (en) 2022-05-24

Family

ID=77681837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110683906.8A Active CN113406008B (en) 2021-06-21 2021-06-21 Photoacoustic imaging device and method based on white light interference

Country Status (1)

Country Link
CN (1) CN113406008B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114264612A (en) * 2022-01-05 2022-04-01 东北大学秦皇岛分校 A PAI-OCT dual-mode imaging system
CN114965285B (en) * 2022-05-11 2025-02-11 重庆邮电大学 Photoacoustic spectroscopy detection system based on direct acoustic wave coupling-multi-path photoacoustic cell
CN115112770B (en) * 2022-07-21 2024-05-03 东北大学秦皇岛分校 Photoacoustic imaging device and method
CN115791632B (en) * 2022-12-09 2024-06-14 东北大学秦皇岛分校 Intravascular ultrasound imaging device and method

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7525661B2 (en) * 2004-02-17 2009-04-28 Andreas Mandelis Laser photo-thermo-acoustic (PTA) frequency swept heterodyned lock-in depth profilometry imaging system
US9833148B2 (en) * 2011-08-11 2017-12-05 University of Washington Through its Center For Commerciallzation Methods and systems for integrated imaging using optical coherence tomography and photoacoustic imaging
CN103163111B (en) * 2013-02-25 2015-11-04 天津大学 A combined fluorescence mesoscopic imaging and OCT early cervical cancer detection system
CN103344569B (en) * 2013-06-21 2015-10-28 中国科学院上海光学精密机械研究所 Polarization complex frequency domain optical coherence tomography method and system
CN105424605B (en) * 2015-11-18 2018-05-01 安徽大学 Photoacoustic spectrum measuring device and method based on low-coherence fiber differential interference non-contact vibration measurement
CN105334262B (en) * 2015-12-04 2017-12-19 东北大学 Noncontact optoacoustic detection method and device based on Through Optical Interference Spectra
CN109363639A (en) * 2018-11-13 2019-02-22 东北大学秦皇岛分校 A high-speed non-contact photoacoustic imaging system based on optical path modulation combined with sensitivity compensation
CN110243763B (en) * 2019-07-02 2020-08-25 东北大学 Non-contact photoacoustic imaging device and method
CN110367941B (en) * 2019-08-20 2022-01-28 东北大学秦皇岛分校 Detection light fusion non-contact photoacoustic-optical coherence tomography dual-mode imaging system
CN112924389A (en) * 2021-01-26 2021-06-08 深圳大学 Multi-mode imaging system and method based on photoacoustic and optical coherence tomography technology

Also Published As

Publication number Publication date
CN113406008A (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN113406008B (en) Photoacoustic imaging device and method based on white light interference
US10548479B2 (en) Optical coherence elastography (OCE) method under acoustic radiation force excitation using OCT doppler variance methods and OCT correlation-based methods
US20080228073A1 (en) System and method for optoacoustic imaging of peripheral tissues
CN110243763B (en) Non-contact photoacoustic imaging device and method
CN108535194B (en) Photoacoustic microscopic imaging system and method based on surface plasma resonance
CN109363639A (en) A high-speed non-contact photoacoustic imaging system based on optical path modulation combined with sensitivity compensation
JP2000517414A (en) Method and apparatus for three-dimensional ultrasonic microscopy using short pulse excitation and three-dimensional ultrasonic microscope used therein
CN115112770B (en) Photoacoustic imaging device and method
Pelivanov et al. A 1 kHz a-scan rate pump-probe laser-ultrasound system for robust inspection of composites
JP2022516551A (en) Equipment and methods for testing subjects
JPH0271147A (en) Non-destructive inspection method and device
CN115517628B (en) A photoacoustic system for detecting circulating tumor cells in freely moving mice
CN113367660B (en) Photoacoustic Doppler flow velocity measuring device and method
CN110367941A (en) A kind of detection light fusion Noncontact optoacoustic-optical coherence tomography double-mode imaging system
CN111466884B (en) Optical coherent elastography device
CN115153453A (en) Quantum-enhanced all-optical photoacoustic signal detection device and method
CN111449629B (en) Optical coherence elastography method and device
CN110169758B (en) An all-optical photoacoustic endoscopic imaging device and method
CN112881297A (en) Speckle interference detection system and method based on photoacoustic cross coupling technology
CN1168980C (en) Method and device for measuring photoacoustic signals in biological tissue using probe ultrasound beam
CN102879335A (en) Portable noninvasive real-time photoacoustic viscoelastic detector
CN108896221B (en) Shock wave signal detection device and method based on Mach-Zehnder interference
Shen et al. Fiber-optic ultrasonic probe based on refractive-index modulation in water
Chitnis et al. An adaptive interferometric sensor for all-optical photoacoustic microscopy
CN205826515U (en) Acoustical signal detector based on surface wave and reflecting light sonomicroscope

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
GR01 Patent grant
GR01 Patent grant
OL01 Intention to license declared
OL01 Intention to license declared