CN113933245A - 一种基于单波长透射式光声显微镜的双组分定量成像法 - Google Patents
一种基于单波长透射式光声显微镜的双组分定量成像法 Download PDFInfo
- Publication number
- CN113933245A CN113933245A CN202110973416.1A CN202110973416A CN113933245A CN 113933245 A CN113933245 A CN 113933245A CN 202110973416 A CN202110973416 A CN 202110973416A CN 113933245 A CN113933245 A CN 113933245A
- Authority
- CN
- China
- Prior art keywords
- component
- sample
- photoacoustic
- target sample
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 45
- 230000005540 biological transmission Effects 0.000 title claims abstract description 22
- 230000031700 light absorption Effects 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000000149 argon plasma sintering Methods 0.000 claims abstract description 15
- 230000003287 optical effect Effects 0.000 claims abstract description 4
- 230000008878 coupling Effects 0.000 claims description 11
- 238000010168 coupling process Methods 0.000 claims description 11
- 238000005859 coupling reaction Methods 0.000 claims description 11
- 239000007788 liquid Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 7
- 238000010895 photoacoustic effect Methods 0.000 claims description 5
- 238000002604 ultrasonography Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000006096 absorbing agent Substances 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 239000000203 mixture Substances 0.000 abstract description 12
- 230000005284 excitation Effects 0.000 abstract description 10
- 230000007547 defect Effects 0.000 abstract description 2
- 230000001575 pathological effect Effects 0.000 abstract description 2
- 239000008280 blood Substances 0.000 description 24
- 210000004369 blood Anatomy 0.000 description 24
- 239000008267 milk Substances 0.000 description 19
- 235000013336 milk Nutrition 0.000 description 18
- 210000004080 milk Anatomy 0.000 description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 238000000386 microscopy Methods 0.000 description 4
- 238000004445 quantitative analysis Methods 0.000 description 4
- 235000020247 cow milk Nutrition 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 2
- 108010054147 Hemoglobins Proteins 0.000 description 2
- 102000001554 Hemoglobins Human genes 0.000 description 2
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical compound O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 description 2
- 108010064719 Oxyhemoglobins Proteins 0.000 description 2
- 241001494479 Pecora Species 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 108010002255 deoxyhemoglobin Proteins 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 239000012466 permeate Substances 0.000 description 2
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000006059 cover glass Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008557 oxygen metabolism Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/01—Arrangements or apparatus for facilitating the optical investigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/01—Arrangements or apparatus for facilitating the optical investigation
- G01N2021/0106—General arrangement of respective parts
- G01N2021/0112—Apparatus in one mechanical, optical or electronic block
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)
- Spectroscopy & Molecular Physics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于单波长透射式光声显微镜的双组份定量成像法，当脉冲激光照射样品后，由于光声效应激发光声信号，超声换能器可以检测到两个携带有样本不同信息的光声信号。通过分析两个光声信号的大小，并根据样品中不同组分光吸收和光散射特性的差异，从而实现对样品中的双组分成分的定量成像。本发明提出的方法，基于光学分辨率的透射式光声显微镜，可以实现高对比度的双组分混合物含量的定量成像；该方法无需不同波长或不同能量激光的多次差异性照射，避免了多个激发源和多次激发导致的许多缺点，在揭示基本生理和病理现象方面具有较高的易用性和实用性。
Description
技术领域
本发明涉及一种基于单波长透射式光声显微镜的双组分定量成像法，是一种利用透射模式光声显微镜在单波长单次照射下实现双组分定量成像的方法。
背景技术
对生化样品中不同组分的定量分析和成像具有重要的应用价值。例如，组织中多种发色团浓度的定量信息可以对许多基本的生理和病理现象提供重要的信息；通过血管中氧合血红蛋白和脱氧血红蛋白的定量成像，还能获得的血氧饱和度图像，从而反映组织的氧代谢状况。
光声显微镜是一种新兴的基于光声效应的生物医学成像技术。当激光脉冲照射样品时，入射光被样品中的发色团(如血红蛋白、脂类、黑色素、胶原、DNA/RNA等)吸收，导致局部温度升高和热弹性膨胀，进而产生超声波，这就是光声效应。光声效应产生的超声波又被称为光声信号。光声信号的强度与样品的光吸收系数成正比。由于生物组织中各种生物分子组份具有特定的光吸收特性，因此，光声显微镜可以有效地提供生物组织功能和分子图像。
光声显微镜对组织中多种成分的成像显示出巨大潜力。基于双波长激发的光声显微镜已经成功地测量了含氧血红蛋白和脱氧血红蛋白的浓度，并进一步提取血氧饱和度。此外，基于分子弛豫效应的光声显微镜，使用不同强度的激光来区分氧合血红蛋白和脱氧血红蛋白。然而，这些方法通常都需要多波长激发、多次激发，因此需要多台激光器，或者波长可调谐激光器，从而实现不同波长或不同能量的差异性照射，增加了系统的复杂性和成本。而且，多次照射会降低测量精度和生物安全性。总之，多组分的定量成像仍然是一个十分有吸引力的挑战。
发明内容
发明目的：为了克服以往多组分定量成像需要多个激发源的局限性，实现在单波长单次照射条件下的双组分定量成像，本发明提供一种基于单波长透射式光声显微镜的双组分定量成像法，采用单波长透射模式光声显微镜扫描样品，在每点激光激发后，聚焦超声换能器可以同时采集到两个光声信号，这两个信号的幅度与样品的光吸收和光散射特性有关，通过分析这两个信号的强度，并且结合样品中两个组分光吸收和光散射特性的不同，可以实现对样品中双组分的量化并成像。
技术方案：为实现上述目的，本发明采用的技术方案为：
一种基于单波长透射式光声显微镜的双组份定量成像法，包括如下步骤：
步骤1：将双组份的目标样本置于样本片上，将样本片安装在单波长透射式光声显微镜上，即先将样本片固定在共聚焦皿底部，再向共聚焦皿内注入耦合液体；球形聚焦超声换能器安置在共聚焦皿上方，物镜安装在共聚焦皿下方，球形聚焦超声换能器与物镜同轴共聚焦对准；脉冲激光通过物镜聚焦在目标样本上，球形聚焦超声换能器表面浸入共聚焦皿中的耦合液体内，以实现良好的声耦合；
步骤2：激光器产生单波长的脉冲激光，脉冲激光通过物镜聚焦在目标样本上；脉冲激光与目标样本相互作用，产生光吸收和光散射；
步骤3：由于光吸收和光散射，目标样本吸收部分脉冲激光的能量，该部分能量经光声效应产生光声信号SII；光声信号SII在共聚焦皿中的耦合液体内传播，球形聚焦超声换能器在tII时刻接收到光声信号SII；
步骤4：目标样本未吸收的剩余部分脉冲激光能量在共聚焦皿中的耦合液体内传播，于tI时刻在球形聚焦超声换能器表面产生光声信号SI，tI时刻早于tII时刻；
步骤5：针对目标样本的每个采样点，球形聚焦超声换能器采集到两个在时序上分离的光声信号，将这两个光声信号转换为电信号，经放大、数字化后导入计算机；
步骤6：使用二维电移台承载样本片移动，对目标样本进行逐点扫描；
步骤7：分析目标样本每个采样点上的两个光声信号，结合目标样本中每种组份的光吸收和光散射特性，定量求解目标样本各组份含量；
步骤8：重复步骤7，以各组份浓度为成像参数，可以得到目标样本中各组份含量的分布图像。
具体的，脉冲激光经过目标样本的光吸收和光散射后，到达球形聚焦超声换能器表面的脉冲激光被减弱，进而引起光声信号SI强度的降低，光声信号SI幅度的降低量与目标样本的光吸收和光散射成正比；目标样本中的光吸收体引起光声信号SII强度的增加，光声信号SII幅度的增加量与目标样本的光吸收成正比；因此，通过分析两个光声信号幅度的大小，结合目标样本中两个组份的光吸收和光散射特性，能够对两个组份的含量进行量化，这是单波长单次照射实现双组分定量成像的基本原理；
设目标样本的两个组份分别为组份A和组份B，在样本片(x,y)位置处两组组份的含量分别为ma(x,y)和mb(x,y)，组份A和组份B的吸收系数分别是εa和εb，组份A和组份B的散射系数分别是σa和σb；
对于光声信号SI，令s1表示目标样本导致的光声信号SI强度的降低量(即无目标样本区域与有目标样本区域光声信号SI的幅度差(如图1(b)所示))，存在如下关系：
(εa+σa)·ma+(εb+σb)·mb＝α1·s1 (1)
其中：α1是与脉冲激光强度、系统增益和球形聚焦超声换能器表面的光吸收特性有光的系统参数；
对于光声信号SII，令s2表示目标样本导致的光声信号SII强度的增加量(即无目标样本区域与有目标样本区域光声信号SII的幅度差(如图1(b)所示))，存在如下关系：
εa·ma+εb·mb＝α2·s2 (2)
其中：α2是与脉冲激光强度、系统增益和球形聚焦超声换能器的响应有关的系统参数；
联立方程(1)和(2)，得到如下方程组：
其中：k11＝(εa+σa)/α1，k12＝(εb+σb)/α1，k21＝εa/α2，k22＝εb/α2。
对于方程(3)中的系数k11、k12、k21、k22，可以根据已知样品组份的光吸收和光散射系数等计算确定；也可以对已知浓度的校准样本进行校准实验，从而确定方程(3)中的系数k11、k12、k21、k22，具体包括如下步骤：
制备n个校准样本，校准样本与目标样本均为含有组份A和组份B的双组份样本，且校准样本中的各组份含量已知，分别记为ma (1)、mb (1)、……、ma (n)、mb (n)；ma (i)和mb (i)分别表示第i个校准样本中组份A和组份B的含量；
将校准样本安装在单波长透射式光声显微镜上，s1 (1)，s2 (1)，……，s1 (n)，s2 (n)为从这些校准样检测到的信号；s1 (i)和s2 (i)分别表示第i个校准样本中光声信号SI强度的降低量和光声信号SII强度的增加量；
可得到如下方程：
M*K＝S (4)
其中：
采用最小二乘法求解方程(4)，从而估计系数k11、k12、k21、k22，即K＝(MTM)-1MTS；由于存在四个未知系数，因而至少需要两个校准样本(即n≥2)才能确定四个未知系数。
具体的，根据系数k11、k12、k21、k22，以及单波长透射式光声显微镜测量获得的信号幅度s1和s2，求解的方程(3)，可得样本中两种组份的含量ma和mb：
进一步计算得到样本中各组分的浓度ca和cb：
其中：ca和cb分别为样本中组份A和组份B的浓度。
有益效果：本发明提供的基于单波长透射式光声显微镜的双组份定量成像法，相对于现有技术，具有如下优势：1、本发明在每次单波长激光单次照射样品后，可以通过同一个球形聚焦超声换能器实现两个光声信号的采集，从而得到混合物中双组分的两幅定量图像；2、本发明对样品中双组分含量进行量化成像的整个过程，无需不同波长或不同能量激光的多次差异性照射，避免了多个激发源和多次激发会引发的许多缺点；3、本发明对双组分含量量化的成像信号由同一个换能器接收，故两幅图像的位置完全匹配度，易于配准叠加。可以通过叠加得到复合图，进一步分析得到更多信息。
附图说明
图1(a)为本发明基于单波长透射式光声显微镜的双组分定量成像系统原理图；
图1(b)为球形聚焦超声换能器表面接收到的光声信号的信号图，s1表示光声信号SI因为目标样本引起的幅度减小量，s2表示光声信号SII因为目标样本导致的幅度增大量；
图2为采用本发明对以血液和牛奶的二元混合物为目标样本的定量成像实验组图：
(a)为目标样本的照片，标记为“1”、“2”、“3”、“4”的四个目标样本中，血液与牛奶含量配比分别为1.0:0.0、0.667:0.333、0.333:0.667、0.0:1.0，标记为“5”的目标样本是水，不含有任何血液或牛奶成分；
(b)是以光声信号SI峰值为成像参量绘制的图像；
(c)是以光声信号SII峰值为成像参量绘制的图像；
(d)是牛奶含量的定量图像；
(e)是血液含量的定量图像；
图3(a)为对图2中牛奶含量的实验结果和实际浓度的定量分析，厚度表示目标样本的厚度，目标样本厚度对定量成像结果几乎没有影响；
图3(b)为对图2中血液含量的实验结果和实际浓度的定量分析，厚度表示目标样本的厚度，目标样本厚度对定量成像结果几乎没有影响；
图4为沿两个发色团边界处的浓度梯度成像组图，各组分在交界处相互扩散、渗透：
(a)为牛奶的浓度图像；
(b)为目标样本的照片；
(c)为血液的浓度图像；
(d)为(a)和(c)的复合图，血液和牛奶分别用伪彩色红色和绿色表示；
(e)为各组分浓度沿图(d)中垂线的分布；
(f)为各组分浓度沿图(d)横线的分布。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1(a)为基于单波长透射式光声显微镜的双组分定量成像系统，图1(b)为每次激发后，同一换能器接收到的两个信号，样品点和参考点(没有样品的区域)之间的幅度差异用s1和s2量化表示。用具有不同光吸收和光散射特性的二元混合物对成像方法进行验证与应用，下面结合实例对本发明做出进一步的说明。
步骤1：实验样品制备。用无菌脱纤维绵羊血和牛奶制备四组二元混合物样品，四组样品中血液和牛奶的体积比分别为3:0、2:1、1:2和0:3，即血液(牛奶)含量为1.0(0.0)、0.667(0.333)、0.333(0.667)、0.0(1.0)。将四组混合物分别注入位于共聚焦皿底、厚度约为0.16mm的四个矩形微管中(图2(a)所示样品1-4)。将水注入最右侧的矩形微管中，作为对照组(图2(a)所示样品5)。共聚焦皿中盛水用作超声耦合。
步骤2：波长532nm、脉冲宽度约为8ns、重复频率10kHz的脉冲激光，通过共焦皿底下、数值孔径为0.4的物镜聚焦在样品上。光电二极管监测激光强度，并触发采集卡记录信号。聚焦超声换能器安置在共聚焦皿上方，并与下方物镜保持同轴共聚焦。该聚焦超声换能器的中心频率22.72MHz，-6dB相对带宽77.36％，有效直径3mm，焦距6.83mm。换能器表面浸入共聚焦皿中的液体内，以实现良好的声耦合。实验装置示意图如图1(a)所示。
步骤3：超声换能器接收每个采样点产生的两个在时序上分离的声信号SI和SII，转换为电信号，经由小信号放大器以46dB的增益放大，并以250MHz的采样频率数字化后，保存到计算机中。
步骤4：二维电移台以扫描间隔50μm承载样品移动，进行逐点扫描。
步骤5：提取信号SII的峰-峰值。图2(c)给出了以SII峰-峰值作为成像参数的图像。
步骤6：提取信号SI的峰-峰值。图2(b)给出了以SI峰-峰值作为成像参数的图像。
步骤7：根据每个采样点检测到的两个信号大小，计算求得到混合物中两个组分(牛奶和血液)的含量，并绘制图像，如图2(d)和(e)所示。
步骤71：逐点计算样品s1值。将无样品区域的信号SI的幅度减去有样品区域信号SI的幅度，两者之差即为样品该点的s1值。
步骤72：逐点计算样品s2值。将有样品区域的信号SII的幅度减去无样品区域信号SII的幅度，两者之差即为样品该点的s1值。
步骤73：利用两个校准样估计方程组(3)的系数。
配比两个校准样，其血液-牛奶体积比已知，分别为0:1和1:0，即ma (1)＝0.0，mb (1)＝1.0，ma (2)＝1.0，mb (2)＝0.0。并采用步骤71-72的方法获得两个校准样品的s1 (1)，s2 (1)，s1 (2)，s2 (2)。代入公式(4)可以求得方程组(3)的系数k11、k12、k21、k22。
步骤74：目标样品未知的ma和mb可以由公式(6)求得：
进一步计算得到样本中各组分的浓度ca和cb：
ca＝ma/(ma+mb)，cb＝mb/(ma+mb)
步骤75：采用步骤74-75的方法，逐点计算目标样品的空间每点的各组分浓度值ma(x,y)和mb(x,y)，并分别以牛奶和血液浓度作为成像参数，绘制混合物中两个组分(牛奶和血液)的定量图像。如图2(d)和(e)所示，牛奶的含量从样品1到样品4逐渐增加(图2(d))，而血液的含量在四个混合物样品中逐一减少(图2(e))。此外，图像显示，第五个样品不含任何血液和牛奶。综上，成像结果与5个样品中牛奶和血液的实际含量趋势相符。
图3给出了对图像进行定量分析的结果。从图2(d)和(e)中提取每个混合物样品中各组分的平均浓度。四个二元混合物样品中，牛奶浓度分别为-0.094±0.143，0.405±0.077，0.661±0.044，0.950±0.020，血液浓度分别为1.094±0.143，0.595±0.077，0.339±0.044，0.050±0.02。这些数值与四个样品的实际浓度相符。另外，考虑到实际样品不可能完全均匀，并且样品的厚度会影响成像性能，我们在另外四个厚度加倍的样品上重复上述实验。计算得到牛奶浓度为-0.031±0.080，0.398±0.070，0.692±0.034，0.953±0.019，血液浓度为1.031±0.080，0.602±0.070，0.308±0.034，0.047±0.019。可以看出，这些值仍然正确地反映了样品的实际浓度。此外，结合上述数据及图3(a)和图3(b)可以看出，不同厚度样品得到的各组分浓度没有显著差异。
此外，本发明还可以用于两个发色团边界处的浓度梯度成像，如图4所示，样品中所有组分(牛奶、血液、水)都相互扩散、渗透，在交界处形成了渐变的浓度梯度。所获得的图像准确地反映了这种浓度梯度，进一步验证了所提出方法的实用性。
步骤1：将无菌脱纤维绵羊血、牛奶和水分别点样在共聚焦皿底，并用盖玻片盖住。从图4(b)样品照片中可见，各组分在交界处相互扩散、渗透。
步骤2：在新样品上重复案例一种的成像步骤，得到牛奶和血液的浓度图像，绘制图4(a)和(c)。
步骤3：将图4(a)和(c)配准后叠加，得到一个复合图，即图4(d)。图4(d)中，血液和牛奶分别用伪彩色红色和绿色表示。
步骤4：提取图4(d)中任意垂线和横线上各组分的浓度，分别绘制在图4(e)和(f)中。如图4(e)和(f)所示，样品中各组分的空间分布与实际情况基本一致。
以上实验结果均表明，基于单波长透射式光声显微镜的双组分定量成像方法，可以通过分析每次单波长激光照射后，同一个球形聚焦超声换能器采集到的两个光声信号大小，结合样品中两个组分光吸收和光散射特性的差异，进而对样品中双组分含量进行量化并成像。本发明提出的方法只需要单波长激光单次照射即可实现双组分定量成像，避免了以往多组分定量成像需要多个激发源的局限性。因此，本发明方法在单波长激光单次照射的条件下实现了双组分的定量成像，具有较高的易用性和实用性。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种基于单波长透射式光声显微镜的双组份定量成像法，其特征在于：包括如下步骤：
步骤1：将双组份的目标样本置于样本片上，将样本片安装在单波长透射式光声显微镜上，即先将样本片固定在共聚焦皿底部，再向共聚焦皿内注入耦合液体；球形聚焦超声换能器安置在共聚焦皿上方，物镜安装在共聚焦皿下方，球形聚焦超声换能器与物镜同轴共聚焦对准；脉冲激光通过物镜聚焦在目标样本上，球形聚焦超声换能器表面浸入共聚焦皿中的耦合液体内；
步骤2：激光器产生单波长的脉冲激光，脉冲激光通过物镜聚焦在目标样本上；脉冲激光与目标样本相互作用，产生光吸收和光散射；
步骤3：由于光吸收和光散射，目标样本吸收部分脉冲激光的能量，该部分能量经光声效应产生光声信号SII；光声信号SII在共聚焦皿中的耦合液体内传播，球形聚焦超声换能器在tII时刻接收到光声信号SII；
步骤4：目标样本未吸收的剩余部分脉冲激光能量在共聚焦皿中的耦合液体内传播，于tI时刻在球形聚焦超声换能器表面产生光声信号SI，tI时刻早于tII时刻；
步骤5：针对目标样本的每个采样点，球形聚焦超声换能器采集到两个在时序上分离的光声信号，将这两个光声信号转换为电信号，经放大、数字化后导入计算机；
步骤6：使用二维电移台承载样本片移动，对目标样本进行逐点扫描；
步骤7：分析目标样本每个采样点上的两个光声信号，结合目标样本中每种组份的光吸收和光散射特性，定量求解目标样本两种组份的含量；
步骤8：重复步骤7，以各组份浓度为成像参数，可以得到目标样本中各组份含量的分布图像。
2.根据权利要求1所述的基于单波长透射式光声显微镜的双组份定量成像法，其特征在于：脉冲激光经过目标样本的光吸收和光散射后，到达球形聚焦超声换能器表面的脉冲激光被减弱，进而引起光声信号SI强度的降低，光声信号SI幅度的降低量与目标样本的光吸收和光散射成正比；目标样本中的光吸收体引起光声信号SII强度的增加，光声信号SII幅度的增加量与目标样本的光吸收成正比；
设目标样本的两个组份分别为组份A和组份B，在样本片(x,y)位置处两组组份的含量分别为ma(x,y)和mb(x,y)，组份A和组份B的吸收系数分别是εa和εb，组份A和组份B的散射系数分别是σa和σb；
对于光声信号SI，令s1表示目标样本导致的光声信号SI强度的降低量，存在如下关系：
(εa+σa)·ma+(εb+σb)·mb＝α1·s1 (1)
其中：α1是与脉冲激光强度、系统增益和球形聚焦超声换能器表面的光吸收特性有光的系统参数；
对于光声信号SII，令s2表示目标样本导致的光声信号SII强度的增加量，存在如下关系：
εa·ma+εb·mb＝α2·s2 (2)
其中：α2是与脉冲激光强度、系统增益和球形聚焦超声换能器的响应有关的系统参数；
联立方程(1)和(2)，得到如下方程组：
其中：k11＝(εa+σa)/α1，k12＝(εb+σb)/α1，k21＝εa/α2，k22＝εb/α2。
3.根据权利要求2所述的基于单波长透射式光声显微镜的双组份定量成像法，其特征在于：对于方程(3)中的系数k11、k12、k21、k22，根据已知样品组份的光吸收和光散射系数计算确定；或者对已知浓度的校准样本进行校准实验，从而确定方程(3)中的系数k11、k12、k21、k22，具体包括如下步骤：
制备n个校准样本，校准样本与目标样本均为含有组份A和组份B的双组份样本，且校准样本中的各组份含量已知，分别记为ma (1)、mb (1)、……、ma (n)、mb (n)；ma (i)和mb (i)分别表示第i个校准样本中组份A和组份B的含量；
将校准样本安装在单波长透射式光声显微镜上，s1 (1)，s2 (1)，……，s1 (n)，s2 (n)为从这些校准样检测到的信号；s1 (i)和s2 (i)分别表示第i个校准样本中光声信号SI强度的降低量和光声信号SII强度的增加量；
可得到如下方程：
M*K＝S (4)
其中：
采用最小二乘法求解方程(4)，从而估计系数k11、k12、k21、k22，即K＝(MTM)-1MTS。
4.根据权利要求2所述的基于单波长透射式光声显微镜的双组份定量成像法，其特征在于：根据系数k11、k12、k21、k22，以及单波长透射式光声显微镜测量获得的信号幅度s1和s2，求解的方程(3)，可得样本中两种组份的含量ma和mb：
进一步计算得到样本中各组分的浓度ca和cb：
其中：ca和cb分别为样本中组份A和组份B的浓度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110973416.1A CN113933245B (zh) | 2021-08-24 | 2021-08-24 | 一种基于单波长透射式光声显微镜的双组分定量成像法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110973416.1A CN113933245B (zh) | 2021-08-24 | 2021-08-24 | 一种基于单波长透射式光声显微镜的双组分定量成像法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113933245A true CN113933245A (zh) | 2022-01-14 |
| CN113933245B CN113933245B (zh) | 2023-06-06 |
Family
ID=79274463
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110973416.1A Active CN113933245B (zh) | 2021-08-24 | 2021-08-24 | 一种基于单波长透射式光声显微镜的双组分定量成像法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113933245B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118641480A (zh) * | 2024-08-14 | 2024-09-13 | 之江实验室 | 一种光声蛋白筛选方法和光声蛋白筛选设备 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007054751A1 (en) * | 2005-11-14 | 2007-05-18 | Szegedi Tudományegyetem | Method and apparatus for the selective photoacoustic detection of contaminants in a gaseous sample by making use of spectrally distant exciting wavelengths |
| US20090054763A1 (en) * | 2006-01-19 | 2009-02-26 | The Regents Of The University Of Michigan | System and method for spectroscopic photoacoustic tomography |
| EP2148183A1 (en) * | 2008-07-25 | 2010-01-27 | Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH) | Imaging device and method for optoacoustic imaging of small animals |
| CN101696931A (zh) * | 2009-10-15 | 2010-04-21 | 清华大学深圳研究生院 | 物质浓度无损检测方法及其装置 |
| US20140142404A1 (en) * | 2008-10-23 | 2014-05-22 | The Washington University | Single-cell label-free photoacoustic flowoxigraphy in vivo |
| US20160305914A1 (en) * | 2013-11-19 | 2016-10-20 | Washington University | Systems and methods of grueneisen-relaxation photoacoustic microscopy and photoacoustic wavefront shaping |
| CN106725348A (zh) * | 2017-02-27 | 2017-05-31 | 集美大学 | 一种吸收体吸收系数测量和同时光声成像的无损检测方法及其装置 |
| US20190192008A1 (en) * | 2016-06-30 | 2019-06-27 | Deutsches Krebsforschungszentrum Stiftung des öffentlichen Rechts | Machine Learning-Based Quantitative Photoacoustic Tomography (PAT) |
| CN111693465A (zh) * | 2020-07-15 | 2020-09-22 | 南京大学 | 一种同时获得光吸收和光散射双对比度的显微成像法 |
| CN112683808A (zh) * | 2020-12-15 | 2021-04-20 | 电子科技大学 | 一种基于光声谱多组分装置与信号处理方法 |
-
2021
- 2021-08-24 CN CN202110973416.1A patent/CN113933245B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007054751A1 (en) * | 2005-11-14 | 2007-05-18 | Szegedi Tudományegyetem | Method and apparatus for the selective photoacoustic detection of contaminants in a gaseous sample by making use of spectrally distant exciting wavelengths |
| US20090054763A1 (en) * | 2006-01-19 | 2009-02-26 | The Regents Of The University Of Michigan | System and method for spectroscopic photoacoustic tomography |
| EP2148183A1 (en) * | 2008-07-25 | 2010-01-27 | Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH) | Imaging device and method for optoacoustic imaging of small animals |
| US20140142404A1 (en) * | 2008-10-23 | 2014-05-22 | The Washington University | Single-cell label-free photoacoustic flowoxigraphy in vivo |
| CN101696931A (zh) * | 2009-10-15 | 2010-04-21 | 清华大学深圳研究生院 | 物质浓度无损检测方法及其装置 |
| US20160305914A1 (en) * | 2013-11-19 | 2016-10-20 | Washington University | Systems and methods of grueneisen-relaxation photoacoustic microscopy and photoacoustic wavefront shaping |
| US20190192008A1 (en) * | 2016-06-30 | 2019-06-27 | Deutsches Krebsforschungszentrum Stiftung des öffentlichen Rechts | Machine Learning-Based Quantitative Photoacoustic Tomography (PAT) |
| CN106725348A (zh) * | 2017-02-27 | 2017-05-31 | 集美大学 | 一种吸收体吸收系数测量和同时光声成像的无损检测方法及其装置 |
| CN111693465A (zh) * | 2020-07-15 | 2020-09-22 | 南京大学 | 一种同时获得光吸收和光散射双对比度的显微成像法 |
| CN112683808A (zh) * | 2020-12-15 | 2021-04-20 | 电子科技大学 | 一种基于光声谱多组分装置与信号处理方法 |
Non-Patent Citations (6)
| Title |
|---|
| AMOS DANIELLI ET AL.: "Single-wavelength functional photoacoustic microscopy in biological tissue", 《OPTICS LETTERS》 * |
| 孙正等: "定量光声层析成像的研究进展", 《发光学报》 * |
| 尚姗姗 等: "基于光谱编码的血管内光声组分成像", 《激光生物学报》 * |
| 燕康;张建辉;张帆;周瑛子;龚小竞;: "高精度光声光谱定量成像系统的设计与应用", 机械工程学报 * |
| 陶超;殷杰;刘晓峻;: "生物组织光声成像技术综述", 数据采集与处理 * |
| 马瑶等: "单波长光声光谱血氧饱和度检测技术研究（英文）", 《光谱学与光谱分析》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118641480A (zh) * | 2024-08-14 | 2024-09-13 | 之江实验室 | 一种光声蛋白筛选方法和光声蛋白筛选设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113933245B (zh) | 2023-06-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Strohm et al. | High frequency label-free photoacoustic microscopy of single cells | |
| JP5197779B2 (ja) | 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法 | |
| Strohm et al. | Single cell photoacoustic microscopy: a review | |
| JP3433508B2 (ja) | 散乱吸収体計測方法及び散乱吸収体計測装置 | |
| Strohm et al. | High resolution ultrasound and photoacoustic imaging of single cells | |
| Hysi et al. | Photoacoustic ultrasound spectroscopy for assessing red blood cell aggregation and oxygenation | |
| EP0592200B1 (en) | Apparatus for measuring optical information in scattering medium and method therefor | |
| JP6366367B2 (ja) | 被検体情報取得装置、被検体情報取得装置の制御方法、および、プログラム | |
| US20090069674A1 (en) | Measurement apparatus | |
| Estrada et al. | Virtual craniotomy for high-resolution optoacoustic brain microscopy | |
| JP2009538418A (ja) | 光音響撮像方法 | |
| CN110361357A (zh) | 一种用于皮肤检测的单阵元光声谱信号获取系统及方法 | |
| JP6461288B2 (ja) | 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法 | |
| WO2018043193A1 (ja) | 情報取得装置および信号処理方法 | |
| Pilatou et al. | Analysis of three-dimensional photoacoustic imaging of a vascular tree in vitro | |
| CN113933245B (zh) | 一种基于单波长透射式光声显微镜的双组分定量成像法 | |
| US20060187533A1 (en) | Method and device for time-gated amplification of photons | |
| CN111693465B (zh) | 一种同时获得光吸收和光散射双对比度的显微成像法 | |
| Trunina et al. | Study of water diffusion in human dentin by optical coherent tomography | |
| Pleitez et al. | Mid-infrared optoacoustic microscopy with label-free chemical contrast in living cells and tissues | |
| Zhang | Submicron-resolution photoacoustic microscopy of endogenous light-absorbing biomolecules | |
| Paltauf | Dual-wavelength optoacoustic imaging | |
| Lengenfelder | Remote photoacoustic sensing using speckle-analysis for biomedical imaging | |
| Weernink | Light fluence marker for quantitative photoacoustic imaging | |
| Abdalsalam | Optical Property Estimation in Multi-Layered Tissues Using Structured Interrogation Frequency-Domain Near-Infrared Spectroscopy |
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 |