CN104873175A - 扩散光学层析和光声层析联合测量***及方法 - Google Patents
扩散光学层析和光声层析联合测量***及方法 Download PDFInfo
- Publication number
- CN104873175A CN104873175A CN201510346284.4A CN201510346284A CN104873175A CN 104873175 A CN104873175 A CN 104873175A CN 201510346284 A CN201510346284 A CN 201510346284A CN 104873175 A CN104873175 A CN 104873175A
- Authority
- CN
- China
- Prior art keywords
- chromatography
- optoacoustic
- signal
- optical
- detection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 84
- 238000005259 measurement Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000003325 tomography Methods 0.000 title abstract description 14
- 238000003384 imaging method Methods 0.000 claims abstract description 56
- 238000001514 detection method Methods 0.000 claims abstract description 44
- 238000004587 chromatography analysis Methods 0.000 claims description 50
- 238000004458 analytical method Methods 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 239000000835 fiber Substances 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 239000007788 liquid Substances 0.000 claims description 8
- 239000004065 semiconductor Substances 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 239000013307 optical fiber Substances 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000002902 bimodal effect Effects 0.000 claims description 2
- 230000005284 excitation Effects 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 3
- 238000013480 data collection Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 13
- 238000010521 absorption reaction Methods 0.000 description 9
- 210000004369 blood Anatomy 0.000 description 8
- 239000008280 blood Substances 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 108010054147 Hemoglobins Proteins 0.000 description 5
- 102000001554 Hemoglobins Human genes 0.000 description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 5
- 229910052760 oxygen Inorganic materials 0.000 description 5
- 239000001301 oxygen Substances 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- 108010064719 Oxyhemoglobins Proteins 0.000 description 2
- 238000003759 clinical diagnosis Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000013481 data capture Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009543 diffuse optical tomography Methods 0.000 description 2
- MOFVSTNWEDAEEK-UHFFFAOYSA-M indocyanine green Chemical compound [Na+].[O-]S(=O)(=O)CCCCN1C2=CC=C3C=CC=CC3=C2C(C)(C)C1=CC=CC=CC=CC1=[N+](CCCCS([O-])(=O)=O)C2=CC=C(C=CC=C3)C3=C2C1(C)C MOFVSTNWEDAEEK-UHFFFAOYSA-M 0.000 description 2
- 229960004657 indocyanine green Drugs 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000001320 near-infrared absorption spectroscopy Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 108010022579 ATP dependent 26S protease Proteins 0.000 description 1
- 238000004497 NIR spectroscopy Methods 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 239000007864 aqueous solution Substances 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000003822 epoxy resin Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229940028435 intralipid Drugs 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000010895 photoacoustic effect Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 229920000647 polyepoxide Polymers 0.000 description 1
- 229920002981 polyvinylidene fluoride Polymers 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000013334 tissue model Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
- A61B6/035—Mechanical aspects of CT
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- Acoustics & Sound (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种扩散光学层析和光声层析的联合测量***及方法,该***包含扩散光学层析测量部分(1)、仿CT空间扫描***(2)、光声层析测量部分(3)及工业控制计算机(4);其中:所述扩散光学层析测量部分(1)和所述光声层析测量部分(3)共同连接仿CT空间扫描***(2)形成一个完整的多模态测量***;所述工业控制计算机(4)中的***自动控制软件用于实现扩散光学信号和光声信号的自动化探测和数据采集。与现有技术相比,本发明克服现有单一模态成像技术的不足,弥补现有技术在空间分辨率和定量性上无法兼得的缺陷,使成像方法兼具有光声层析技术的空间分辨率和扩散光学层析技术的定量精度,继而实现高质量的光学参数成像功能。
Description
技术领域
本发明涉及光学成像技术在生物医学方向上的应用领域,特别涉及一种扩散光学层析和光声层析联合测量***、多模态测量方法。
背景技术
随着生物、医学、光电技术及计算机技术的发展,各种先进医学影像设备不断呈现,为临床诊断提供了多种模式的医学图像。由于成像机理不同和人体组织结构的高度复杂性,单一模式的医学图像不能提供医生所需要的足够信息,不同模态的医学图像提供的信息是不能覆盖、相互补充的。因此,在临床诊断和治疗中,对于多模态医学成像技术的发展有着十分重要的意义。
随着多种模式医学成像手段的发展,现代医学研究发展方向也逐渐地发生转变,即从通过传统的解剖学成像方法对已形成的疾病进行诊断,发展至利用医学功能成像技术揭示疾病的形成机制。光学成像作为一种功能成像手段,因在灵敏度、特异性、对比度、定量性和辐射安全性等关键指标上具有其它模态无法兼有的综合优势而倍受关注,已在生物医学研究领域得到广泛应用。尤其随着近红外组织光谱技术(near-infrared spectroscopy,NIRS)的兴起,针对活体组织的光学功能检测技术得到了迅速发展。其中,建立在NIRS基础上的扩散光学层析成像技术(diffuse opticaltomography,DOT)是针对组织体生化功能信息的无创检测技术,可以实现活体组织深部生理和病理指标相关信息的三维实时定量获取。但是,受到组织体高散射特性的影响,该类纯光学成像技术在成像深度和空间分辨率上无法兼得,其空间分辨率被限制在成像组织尺寸的5%-10%;基于光声效应的光声层析成像技术(photoacoustictomography,PAT)方法有效地结合了光学成像和超声成像的优点,是一种以超声作为媒介的生物光子成像方法。该方法采用脉冲激光照射生物组织,组织吸收光能量后产生超声(光声)信号,通过测量携有生物组织光学吸收特征信息的光声信号,能够重建出组织体的光能吸收分布图像。由于超声波低衰减和高穿透性等特点,大幅度降低了组织体散射效应的影响,PAT方法能够对深层组织进行高空间分辨率和高对比度的结构和功能成像。但是,传统PAT技术通过重建初始声压计算组织体光能吸收密度的空间分布,并假设组织体的吸收系数与光能吸收密度成线性关系,最终获得吸收系数重建图像。当面向深层组织成像时,由于复杂生物组织体内吸收和散射分布不均匀性,这种线性假设关系与实际情况不符,组织体内光能分布特征成为了影响光学参数重建定量性的关键问题。
因此,寻找一种既保有光学功能成像的高对比度优势,又同时满足成像定量性和高空间分辨率要求的医学成像方法,是生物医学成像领域亟待解决的问题。换而言之,该方法需要具有DOT成像技术的高定量性以及PAT成像技术的高空间分辨率特性,即保有两种成像方式所长,互相补充不足,从而得到较为理想的定量性和空间分辨率的光学功能图像。把不同模态成像技术提供的互补信息综合起来,不仅可以在一幅图像上同时表达多幅图像的信息,还可以通过理论学上的创新形成一种全新的多模态成像方法,那么此方法所提供的图像可以为医生做出确切诊断提供愈来愈详细和精确的生理病理信息,也能为医生对病灶的观察和疾病的诊断提供更直观、全面和清晰的判断依据,提高疾病的检出率。
发明内容
为了克服上述现有技术的问题,本发明提出了一种扩散光学层析和光声层析联合测量***及方法,利用扩散光学层析和光声层析的联合测量的手段,实现多模态测量,最终可获得目标体的光学参数在二维或三维空间上的分布图像。
本发明提出一种扩散光学层析和光声层析的联合测量***,该***包含扩散光学层析测量部分(1)、仿CT空间扫描***(2)、光声层析测量部分(3)及工业控制计算机(4);其中:所述扩散光学层析测量部分(1)和所述光声层析测量部分(3)共同连接仿CT空间扫描***(2)形成一个完整的多模态测量***;所述工业控制计算机(4)中的***自动控制软件、数据处理和图像重建软件用于实现扩散光学信号和光声信号的自动化探测、数据采集及预处理和图像重建;其中,自动控制软件用于实现对扩散光学层析和光声层析联合测量***各组成部分的时序控制,以实现测量过程的自动化;数据处理和图像重建软件用于对联合测量***测量数据进行预处理,并采用图像重建算法对扩散光测量数据与光声测量数据进行图像重建,从而得到目标体扩散光学层析和光声层析双模态联合重建图像;
所述扩散光学层析测量部分(1)进一步还包括光学接收***(11)、稳态光源***(12)、滤光***(13)和检测***(14);所述光声层析测量部分(3)进一步还包括光声探测***(31)、脉冲光源***(32)和数据采集***(33)。
所述仿CT空间扫描***包括水箱(21)、浸入水箱中设置的电控升降台和旋转台(22)、成像腔(23)以及用于挂载第一超声换能器(311)、第二超声换能器(312)以实现信号环形扫描的电控环形探测旋转台(24),目标体(5)置于圆柱形的成像腔(23)内,成像腔(23)维持圆柱形状且能保证超声信号有效传播;成像腔(23)固定于电控升降台和旋转台(22)上,通过对升降间距和旋转角度的控制,调整不同成像面(6)上不同光源入射角度,即实现空间采样率的任意设定;所述水箱(21)中承载的水作为超声耦合液,将整个扫描***浸于耦合液中,用于声学信号的传播及减少超声信号的衰减。
所述扩散光学层析测量部分(1)还包括光学接收***(11)、稳态光源***(12)、电控滤光***(13)和检测***(14);其中:
所述光学接收***(11)由目标体(5)对侧的同探测平面位置上等角度间隔均匀分布8根探测光纤组成,采用与成像腔同轴的半圆形光纤支架固定;
所述稳态光源***(12)包括两个波长分别为λ1和λ2的第一稳态半导体激光器(121)、第二稳态半导体激光器(122),且连接同轴入射光纤;入射光经过准直器(123)后垂直入射至目标体(5)表面;分别接入到8:1光开关(111)的不同通道内,8:1光开关(111)出射光经准直器(112)至电控滤光***(13);
所述电控滤光***(13)包括滤光轮控制器(131)和电控滤光轮(132),出射光通过电控滤光***入射至检测***(14);
所述检测***(14)包括光电倍增管(141)与计数模块(142),所述光电倍增管(141)接收光信号并转换成电信号,由计数模块(142)进行处理与记录。
所述光声层析测量部分(3)还包括光声探测***(31)、脉冲光源***(32)和数据采集***(33);其中:
所述脉冲光源***(32)采用脉冲激光器(321)的分布式多束入射方式来实现高能量脉冲入射光对目标体(5)表面的均匀照射,连接一组双臂光纤束(322),光纤束传导激光分别均匀照射成像面的1/4半圆域;
所述光声探测***(31)由挂载于电控环形探测旋转台(24)上的第一超声换能器(311)和第二超声换能器(312)组成,相对入射光源位于反射测量面;其探头的环形扫描位置由电控环形探测旋转台(24)来调整;目标体(5)出射的光声信号由第一超声换能器(311)和第二超声换能器(312)探测后转化成电信号,传送至数据采集***(33);
所述数据采集***(33)包括前置放大器(331)与高速数据采集卡(332)。
本发明还提出一种扩散光学层析和光声层析的联合测量方法,该方法包括以下步骤:
步骤(1)、将选定波长的稳态光源***连接到同轴入射光纤,入射光纤经过准直器后垂直入射至目标体表面,在其对侧的同平面位置上等角度间隔均匀分布8根探测光纤,分别接入到多通道光开关的不同通道内,通过8:1光开关的切换获得不同探测位置的出射光信号;8:1光开关出射光经准直器准直后通过电控滤光轮后入射至光电倍增管,光电倍增管接收光信号并转换成电信号由计数模块进行处理与记录。
步骤(2)、脉冲激光器作为光声激发光源产生脉冲激光,由双臂光纤束将激发光传送至成像面的不同入射位置,垂直入射至待测目标体表面,由脉冲激光器激发产生的光声信号经过超声耦合液,传送到超声换能器被转换成电信号,电信号经过前置放大器后,由高速数据采集卡进行数据采集,最终将采集的数据传送给工业控制计算机进行处理;通过电控环形探测旋转台调整超声换能器的环形扫描位置,重复以上的测量步骤,最终可以得到不同角度投影下待测目标体出射的光声信号;
步骤(3)、通过步骤(1)获得单个光源入射角度下的扩散光学信号,然后再通过步骤(2)采集到待测目标体中激发出的光声信号;利用工业控制计算机中的***自动控制软件来调节仿CT空间扫描***中的电控旋转台,从而获得新的光源入射角度,重复步骤(1)和(2),获得一个二维成像面上任意光源入射角度下的扩散光学信号和光声信号;再利用***自动控制软件来调节电控升降台和旋转台,继续重复以上步骤,便可以获得待测目标体任意成像面上的信息,将所有成像面上的信息组合起来便得到了整个目标体的三维信息;将获取的信息传送给数据处理及图像重建软件,最终可获得目标体的光学参数在二维或三维空间上的分布图像。
与现有技术相比,本发明克服现有单一模态成像技术的不足,弥补现有技术在空间分辨率和定量性上无法兼得的缺陷,使成像方法兼具有光声层析技术的空间分辨率和扩散光学层析技术的定量精度,继而实现高质量的光学参数成像功能。
附图说明
图1为扩散光学层析和光声层析的联合测量***框图;其中附图标记:1、扩散光学层析测量部分;2、仿CT空间扫描***;3、光声层析测量部分;4、工业控制计算机;11、光学接收***;12、稳态光源***;13、电控滤光***;14、检测***;31、光声探测***;32、脉冲光源***;33、数据采集***;
图2为双波长双频宽扩散光学层析测量部分和光声层析联合测量部分组成的联合测量***结构示意图;其中附图标记:111、8:1光开关;112、准直器;121、第一稳态半导体激光器;122、第二稳态半导体激光器;131、滤光轮控制器和132、电控滤光轮;141、光电倍增管;142、计数模块;311、第一超声换能器;312、第二超声换能器;321、脉冲激光器;322、双臂光纤束;331、前置放大器;332、高速数据采集卡;5、目标体;
图3为仿CT空间扫描***的结构示意图;其中附图标记:21、水箱;22、电控升降台和旋转台;23、成像腔;24、电控环形探测旋转台;6、成像面。
具体实施方式
下面将结合附图对本发明的具体实施方式进行详细描述,这些实施方式若存在示例性的内容,不应解释成对本发明的限制。
实施例1:
如图2所示,本发明的具体实施案例1为双波长双频宽扩散光学层析测量部分和光声层析联合测量部分组成的联合测量***结构。
扩散光学层析测量部分实施过程(如图1和图2所示):两个波长分别为λ1和λ2的第一稳态半导体激光器(121)、第二稳态半导体激光器(122)(平均功率<10mW)连接芯径62.5μm的同轴入射光纤,入射光经过准直器(123)后垂直入射至目标体(5)表面(设该入射角度为0°),由其对面101.25°~258.75°的同平面位置上等角度间隔均匀分布8根芯径500μm的探测光纤构组成的光学接收***(11)接收(采用与成像腔同轴的半圆形光纤支架固定),分别接入到8:1光开关(111)的不同通道内(通过光开关的切换获得不同探测位置的光信号),光开关出射光经准直器(112)至电控滤光***,电控滤光***包括滤光轮控制器(131)和实现不同滤光片切换功能的电控滤光轮(132)(载有OD=1~5的不同衰减强度的中性密度滤光片以衰减过强扩散光信号使得探测***工作在线性区域),出射光通过电控滤光***入射至检测***,检测***包括光电倍增管(PMT)(141)与计数模块(142)(光子计数器),PMT接收光信号并转换成电信号由计数模块(142)(采用可编程逻辑技术实现计数门宽和积分时间可调式设计)进行处理与记录。光声层析测量部分实施过程(如图1和图2所示):脉冲光源***采用分布式多束入射方式来实现高能量入射光对目标体(5)表面的均匀照射。脉冲光源***(32)采用脉冲激光器(321)(调Q开关Nd:YAG激光器(波长532nm)泵浦的OPO激光器,波长可调谐范围600~900nm),拟选择与DOT测量部分相近的两个波长,脉冲宽度≤10ns,重复频率10/20Hz,脉冲能量~50mJ。脉冲光源***连接一组双臂光纤束(322),光纤束传导激光分别均匀照射成像面的1/4半圆域。光声探测***(31)由两个宽频超声换能器(311)、(312)(采用水浸非聚集超声探头)组成;两个超声换能器(311)、(312)相对入射光源位于反射测量面,以保证近光源区域光声信号测量的准确性,并结合仿CT扫描模式实现目标体整体探测区域光声信号的有效覆盖。每个超声换能器的接收面直径约为6mm,中心频率分别为2.25MHz和10MHz,通频带宽分别为1.35MHz和7.28MHz。为了得到相同光源下的多角度超声投影,利用电控环形探测旋转台(24)来调整超声探头的扫描位置,调整方向范围分别为成像面的第一象限和第四象限内(设DOT部分的稳态光源入射角度为0°),通过选择不同角度间距调整相同光源入射角度下的扫描投影数。目标体(5)出射的光声信号由两个超声换能器(311)、(312)探测后转化成电信号,再传送至数据采集***(33)。数据采集***包括多通道前置信号放大器(331)及多通道高速数据采集卡(332)(设置采样频率为100MHz);超声换能器组将接收到的光声信号送入前置信号放大器的不同通道中,再通过多通道高速数据采集卡(332)采集不同通道的信号,最终将数据传送到工业控制计算等待处理。前置信号放大器的频率响应和数据采集卡的采样频率根据超声换能器的主频来进行选择。
图3为仿CT空间扫描***的结构示意图:待测目标体置于圆柱形的成像腔(23)内,成像腔(23)拟采用环氧树脂材料和PVDF膜共同制成,既能够维持圆柱形状又能保证超声信号的有效传播,其直径约为2-3cm,壁厚≤1mm。成像腔(23)固定于电控升降台和旋转台(22)上,通过对升降间距和旋转角度的控制,调整不同成像面(6)上不同光源入射角度,即采用此仿CT的扫描模式来实现空间采样率的任意设定。此外,为了有利于声学信号的传播及减少超声信号的衰减,整个扫描测量结构浸于水箱(21)中。利用电控升降台和旋转台来实现三维空间上的多角度扫描,通过调整扫描间隔和范围来自由控制空间采样密度。水箱中承载的水作为超声耦合液,将整个扫描***浸于耦合液中,从而有利于光声信号传播以及减少超声信号的衰减。
实施例2:
利用实施例1所述***进行生物组织模型的扩散光学层析和光声层析联合测量。该***可实现对病灶区域高对比度和高空间分辨率成像,能够对肿瘤新生血管和血氧情况进行监测,以及内源或外源特异性标志物的在体观测。举例说明,在活体小动物实验中采用吲哚菁绿(indocyanine green,ICG)作为外源吸收对比剂时,工作波长λ1根据ICG的吸收光谱峰值波长设置为785nm。此外,为有效获取组织体内血氧信息,***工作波长λ2根据血红蛋白等消光点(约805nm)拟设置为830nm。测量时,将待测目标体置于圆柱型成像腔中,成像腔与目标体之间的空隙采用仿组织体匹配溶液(拟采用Intralipid水溶液加入India Ink混合制成)进行填充。参照实施例1中扩散光学层析测量部分和光声层析测量部分的实施步骤,测量过程首先通过扩散光学层析测量部分获得工作波长为λ1的单个光源入射角度下的扩散光学信号;再通过光声层析测量部分采集目标体激发出的光声信号;接着利用工业控制计算机中的***自动控制软件来调节仿CT空间扫描***中的电控旋转台,从而获得新的光源入射角度;重复测量步骤,可分别获得一组二维成像断层上的扩散光学信号和光声信号;然后在同一成像断层完成工作波长为λ2的扩散光学层析测量和光声层析测量,工作波长为λ2的测量过程与λ1波长下的测量过程相同;再利用***自动控制软件来调节电控升降台,继续重复以上步骤,便可以获得目标体其他成像断层上的信息,将所有断层信息整合得到目标体的三维信息。将获取的信息传送给数据处理及图像重建软件,最终可获得目标体分别在λ1和λ2工作波长下的光学参数的二维或三维空间上分布图像。本实施例由采用双波长测量方式分别获得在λ1和λ2工作波长下的光学吸收参数和分布,继而进行血氧饱和度的计算。血氧饱和度SaO2是指血液中氧合血红蛋白的容量占全部可结合的血红蛋白的百分比,血红蛋白具有氧合血红蛋白HbO2和还原血红蛋白Hb两种形式。由于在该双波长下生物组织中的主要光学吸收物质为血红蛋白,则光学吸收参数与HbO2和Hb浓度关系式可表示为
其中和分别为血液中HbO2的摩尔吸光系数和浓度;εHb和CHb分别是血液中Hb的摩尔吸光系数和浓度。由公式(1)可得
即可分别求得HbO2和Hb的浓度,从而可以计算出组织体内血氧饱和度,具体公式为
Claims (5)
1.一种扩散光学层析和光声层析的联合测量***,其特征在于,该***包含扩散光学层析测量部分(1)、仿CT空间扫描***(2)、光声层析测量部分(3)及工业控制计算机(4);其中:所述扩散光学层析测量部分(1)和所述光声层析测量部分(3)共同连接仿CT空间扫描***(2)形成一个完整的多模态测量***;所述工业控制计算机(4)中的***自动控制软件、数据处理和图像重建软件用于实现扩散光学信号和光声信号的自动化探测、数据采集及预处理和图像重建;其中,自动控制软件用于实现对扩散光学层析和光声层析联合测量***各组成部分的时序控制,以实现测量过程的自动化;数据处理和图像重建软件用于对联合测量***测量数据进行预处理,并采用图像重建算法对扩散光测量数据与光声测量数据进行图像重建,从而得到目标体扩散光学层析和光声层析双模态联合重建图像;
所述扩散光学层析测量部分(1)进一步还包括光学接收***(11)、稳态光源***(12)、滤光***(13)和检测***(14);所述光声层析测量部分(3)进一步还包括光声探测***(31)、脉冲光源***(32)和数据采集***(33)。
2.如权利要求1所述的扩散光学层析和光声层析的联合测量***,其特征在于,所述仿CT空间扫描***包括水箱(21)、浸入水箱中设置的电控升降台和旋转台(22)、成像腔(23)以及用于挂载第一超声换能器(311)、第二超声换能器(312)以实现信号环形扫描的电控环形探测旋转台(24),目标体(5)置于圆柱形的成像腔(23)内,成像腔(23)维持圆柱形状且能保证超声信号有效传播;成像腔(23)固定于电控升降台和旋转台(22)上,通过对升降间距和旋转角度的控制,调整不同成像面(6)上不同光源入射角度,即实现空间采样率的任意设定;所述水箱(21)中承载的水作为超声耦合液,将整个扫描***浸于耦合液中,用于声学信号的传播及减少超声信号的衰减。
3.如权利要求2所述的扩散光学层析和光声层析的联合测量***,其特征在于,所述扩散光学层析测量部分(1)还包括光学接收***(11)、稳态光源***(12)、电控滤光***(13)和检测***(14);其中:
所述光学接收***(11)由目标体(5)对侧的同探测平面位置上等角度间隔均匀分布8根探测光纤组成,采用与成像腔同轴的半圆形光纤支架固定;
所述稳态光源***(12)包括两个波长分别为λ1和λ2的第一稳态半导体激光器(121)、第二稳态半导体激光器(122),且连接同轴入射光纤;入射光经过准直器(123)后垂直入射至目标体(5)表面;分别接入到8:1光开关(111)的不同通道内,8:1光开关(111)出射光经准直器(112)至电控滤光***(13);
所述电控滤光***(13)包括滤光轮控制器(131)和电控滤光轮(132),出射光通过电控滤光***入射至检测***(14);
所述检测***(14)包括光电倍增管(141)与计数模块(142),所述光电倍增管(141)接收光信号并转换成电信号,由计数模块(142)进行处理与记录。
4.如权利要求2所述的扩散光学层析和光声层析的联合测量***,其特征在于,所述光声层析测量部分(3)还包括光声探测***(31)、脉冲光源***(32)和数据采集***(33);其中:
所述脉冲光源***(32)采用脉冲激光器(321)的分布式多束入射方式来实现高能量脉冲入射光对目标体(5)表面的均匀照射,连接一组双臂光纤束(322),光纤束传导激光分别均匀照射成像面的1/4半圆域;
所述光声探测***(31)由挂载于电控环形探测旋转台(24)上的第一超声换能器(311)和第二超声换能器(312)组成,相对入射光源位于反射测量面;其探头的环形扫描位置由电控环形探测旋转台(24)来调整;目标体(5)出射的光声信号由第一超声换能器(311)和第二超声换能器(312)探测后转化成电信号,传送至数据采集***(33);
所述数据采集***(33)包括前置放大器(331)与高速数据采集卡(332)。
5.一种扩散光学层析和光声层析的联合测量方法,其特征在于,该方法包括以下步骤:
步骤(1)、将选定波长的稳态光源***连接到同轴入射光纤,入射光纤经过准直器后垂直入射至目标体表面,在其对侧的同平面位置上等角度间隔均匀分布8根探测光纤,分别接入到多通道光开关的不同通道内,通过8:1光开关的切换获得不同探测位置的出射光信号;8:1光开关出射光经准直器准直后通过电控滤光轮后入射至光电倍增管,光电倍增管接收光信号并转换成电信号由计数模块进行处理与记录。
步骤(2)、脉冲激光器作为光声激发光源产生脉冲激光,由双臂光纤束将激发光传送至成像面的不同入射位置,垂直入射至待测目标体表面,由脉冲激光器激发产生的光声信号经过超声耦合液,传送到超声换能器被转换成电信号,电信号经过前置放大器后,由高速数据采集卡进行数据采集,最终将采集的数据传送给工业控制计算机进行处理;通过电控环形探测旋转台调整超声换能器的环形扫描位置,重复以上的测量步骤,最终可以得到不同角度投影下待测目标体出射的光声信号;
步骤(3)、通过步骤(1)获得单个光源入射角度下的扩散光学信号,然后再通过步骤(2)采集到待测目标体中激发出的光声信号;利用工业控制计算机中的***自动控制软件来调节仿CT空间扫描***中的电控旋转台,从而获得新的光源入射角度,重复步骤(1)和(2),获得一个二维成像面上任意光源入射角度下的扩散光学信号和光声信号;再利用***自动控制软件来调节电控升降台和旋转台,继续重复以上步骤,便可以获得待测目标体任意成像面上的信息,将所有成像面上的信息组合起来便得到了整个目标体的三维信息;将获取的信息传送给数据处理及图像重建软件,最终可获得目标体的光学参数在二维或三维空间上的分布图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510346284.4A CN104873175A (zh) | 2015-06-19 | 2015-06-19 | 扩散光学层析和光声层析联合测量***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510346284.4A CN104873175A (zh) | 2015-06-19 | 2015-06-19 | 扩散光学层析和光声层析联合测量***及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104873175A true CN104873175A (zh) | 2015-09-02 |
Family
ID=53941058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510346284.4A Pending CN104873175A (zh) | 2015-06-19 | 2015-06-19 | 扩散光学层析和光声层析联合测量***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104873175A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107485408A (zh) * | 2017-09-13 | 2017-12-19 | 北京大学 | 一种实现xct和光声成像的双模成像***及其成像方法 |
CN108366736A (zh) * | 2015-12-21 | 2018-08-03 | 皇家飞利浦有限公司 | 用于组织状况测量的设备 |
CN108553080A (zh) * | 2017-12-26 | 2018-09-21 | 天津大学 | 一种面向小动物皮下肿瘤的立体扫描光声介观成像*** |
CN109187423A (zh) * | 2018-09-10 | 2019-01-11 | 西安电子科技大学 | 一种扩散光学层析成像*** |
CN111053534A (zh) * | 2019-12-17 | 2020-04-24 | 电子科技大学 | 一种基于荧光与光声层析成像的动态脑功能双模态成像*** |
CN111466875A (zh) * | 2020-03-12 | 2020-07-31 | 西安电子科技大学 | 一种旋转式扩散光学成像*** |
CN114027802A (zh) * | 2022-01-10 | 2022-02-11 | 杭州纬壹医疗科技有限公司 | 一种扩散光学层析成像*** |
CN115844331A (zh) * | 2022-12-02 | 2023-03-28 | 天津大学 | 一种多角度的光声层析成像***及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080123083A1 (en) * | 2006-11-29 | 2008-05-29 | The Regents Of The University Of Michigan | System and Method for Photoacoustic Guided Diffuse Optical Imaging |
EP2036487A2 (en) * | 2007-09-12 | 2009-03-18 | Canon Kabushiki Kaisha | Measurement method and measurement apparatus |
-
2015
- 2015-06-19 CN CN201510346284.4A patent/CN104873175A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080123083A1 (en) * | 2006-11-29 | 2008-05-29 | The Regents Of The University Of Michigan | System and Method for Photoacoustic Guided Diffuse Optical Imaging |
EP2036487A2 (en) * | 2007-09-12 | 2009-03-18 | Canon Kabushiki Kaisha | Measurement method and measurement apparatus |
Non-Patent Citations (5)
Title |
---|
ANG LI: "Reconstructing chromospheres concentration images directly by continuous-wave diffuse optical tomography", 《OPTICS LETTERS》 * |
CHEN XU ET AL: "Indocyanine green enhanced co-registered diffuse optical tomography and photoacoustic tomography", 《JOURNAL OF BIOMEDICAL OPTICS》 * |
DAVID A. BOAS ET AL: "Imaging the Body with Diffuse Optical Tomography", 《 IEEE SIGNAL PROCESSING MAGAZINE》 * |
DAVID A. BOAS ET AL: "Volumetric diffuse optical tomography of brain activity", 《OPTICS LETTERS》 * |
YIHAN WANG ET AL: "Photoacoustic tomography guided diffuse optical tomography for small-animal model", 《PROC. OF SPIE》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108366736A (zh) * | 2015-12-21 | 2018-08-03 | 皇家飞利浦有限公司 | 用于组织状况测量的设备 |
CN107485408A (zh) * | 2017-09-13 | 2017-12-19 | 北京大学 | 一种实现xct和光声成像的双模成像***及其成像方法 |
CN107485408B (zh) * | 2017-09-13 | 2021-01-26 | 北京大学 | 一种实现xct和光声成像的双模成像***及其成像方法 |
CN108553080A (zh) * | 2017-12-26 | 2018-09-21 | 天津大学 | 一种面向小动物皮下肿瘤的立体扫描光声介观成像*** |
CN109187423A (zh) * | 2018-09-10 | 2019-01-11 | 西安电子科技大学 | 一种扩散光学层析成像*** |
CN111053534A (zh) * | 2019-12-17 | 2020-04-24 | 电子科技大学 | 一种基于荧光与光声层析成像的动态脑功能双模态成像*** |
CN111466875A (zh) * | 2020-03-12 | 2020-07-31 | 西安电子科技大学 | 一种旋转式扩散光学成像*** |
CN111466875B (zh) * | 2020-03-12 | 2021-08-24 | 西安电子科技大学 | 一种旋转式扩散光学成像*** |
CN114027802A (zh) * | 2022-01-10 | 2022-02-11 | 杭州纬壹医疗科技有限公司 | 一种扩散光学层析成像*** |
CN115844331A (zh) * | 2022-12-02 | 2023-03-28 | 天津大学 | 一种多角度的光声层析成像***及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20220054017A1 (en) | Laser Optoacoustic Ultrasonic Imaging System (LOUIS) and Methods of Use | |
CN104873175A (zh) | 扩散光学层析和光声层析联合测量***及方法 | |
Wang et al. | Photoacoustic microscopy and computed tomography: from bench to bedside | |
Nyayapathi et al. | Dual scan mammoscope (DSM)—a new portable photoacoustic breast imaging system with scanning in craniocaudal plane | |
JP5541662B2 (ja) | 被検体情報取得装置およびその制御方法 | |
US20090054763A1 (en) | System and method for spectroscopic photoacoustic tomography | |
US20100094134A1 (en) | Method and apparatus for medical imaging using near-infrared optical tomography combined with photoacoustic and ultrasound guidance | |
CN101612034B (zh) | 重构混浊介质光学参数的时间分辨测量***及方法 | |
US20180078143A1 (en) | Photo-acoustic imaging apparatus and methods of operation | |
JP2012179348A (ja) | 音響波取得装置 | |
US20130190594A1 (en) | Scanning Optoacoustic Imaging System with High Resolution and Improved Signal Collection Efficiency | |
JP2007267848A (ja) | 腫瘍検査装置 | |
Wang et al. | Noninvasive reflection mode photoacoustic imaging through infant skull toward imaging of neonatal brains | |
CN105596011B (zh) | 一种无创血糖检测装置 | |
Thompson et al. | Laser-induced ultrasound transmitters for large-volume ultrasound tomography | |
CN115397337A (zh) | 血红蛋白浓度测定***、经***探针、附件以及血红蛋白浓度测定方法 | |
Sinha et al. | Differentiation between malignant and normal human thyroid tissue using frequency analysis of multispectral photoacoustic images | |
CN211131008U (zh) | 一种基于多光谱超声调制的乳腺肿瘤多模态成像装置 | |
Xia et al. | Recent progress in multimodal photoacoustic tomography | |
TW201629490A (zh) | 陣列式近場光學高散射材料檢測方法 | |
JP2019111435A (ja) | 情報取得装置 | |
CN211962001U (zh) | 一种高光谱深层三维散射光影像仪 | |
Saxena et al. | From bench to market: commercialization of photoacoustic imaging | |
US20220183564A1 (en) | Photoacoustic breast imaging system and method | |
Haisch et al. | OPUS–Optoacoustic imaging combined with conventional ultrasound for breast cancer detection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150902 |
|
RJ01 | Rejection of invention patent application after publication |