基于声学散射透镜的低成本光声断层成像方法.pdf
《基于声学散射透镜的低成本光声断层成像方法.pdf》由会员分享,可在线阅读,更多相关《基于声学散射透镜的低成本光声断层成像方法.pdf(8页珍藏版)》请在咨信网上搜索。
1、NATURASCIENCEMar.,20232023年3 月JOURNALCOHNANINNIVERSITYVol.59,No.2第5 9 卷第2 期南京大学学报(自然科学)DOI:10.13232/ki.jnju.2023.02.005基于声学散射透镜的低成本光声断层成像方法何君君1,殷杰1.2*,刘广东1,马倩倩1(1.南京科技职业学院,南京,2 10 0 48;2.南京大学声学研究所,南京,2 10 0 93)摘要:光声断层成像需要使用换能器阵列及多通道系统,硬件成本较高,制约了其广泛应用。为降低成本,设计了一种基于声学散射透镜的光声断层成像方法,利用声学多重散射过程中的多路径效应达到增
2、大探测阵列等效数值孔径的作用,从而减少信号探测过程所需换能器数量.通过数值模拟实验,首先分析了散射透镜格林函数的精确测量方法,随后利用四个散射透镜共同作用实现了高质量光声断层成像,并对影响成像质量的相关因素进行了具体研究,提出的散射透镜无须复杂制作工艺或特殊材料,系统只需四个超声换能器即可实现快速成像,可降低光声断层成像系统的硬件成本。关键词:光声断层成像,低成本,散射透镜,格林函数中图分类号:0 42 9文献标志码:ALow-cost photoacoustic tomography based on acoustic scattering lensesHe Junjun,Yin Jiel.
3、2*,Liu Guangdong,Ma Qianqian(1.Nanjing Polytechnic Institute,Nanjing,210048,China;2.Key Laboratory of Modern Acoustics,MOE,Institute of Acoustics,Nanjing University,Nanjing,210093,China)Abstract:Photoacoustic tomography(PAT)is a hybrid biomedical imaging modality based on photoacoustic(PA)effect,itc
4、ombines optical excitation with acoustic detection.As a no-ionizing and non-invasive technology,it can provide bothstructural and chemical informations in pre-clinical and clinical applications,including tumor diagnosis,drug deliveringmonitoring,little animal whole-body imaging,brain funtional imagi
5、ng and so on.As a result,PAT has become one of thefastest growing imaging method in the past decade.One chanllenge that has continuously attracted the attention of researchersand engineerers in biomedical imaging field is to reduce the cost without degradation of performance.PAT system usuallyinvolv
6、es transducer array and multi-channels system,which lead to a high cost and bring limitations on widely application.Aswe know,with the known Greens function,multi-paths effect of multi-scattering can enlarge equivalent numerical aperture ofthe detecting array,and PA images can be reconstructed with
7、much less transducers by time reversal algorithm.In this paper,we propose a PAT modality based on acoustic scattering lenses,which utilizing the multi-paths effect of multi-scattering.Firstly,a method for accurate of Greens function measuring is verified,quantitative studies indicate that,the accurr
8、acy ofGreens function measuring can be improved by choosing the appropriate coherent gathering zones and increasing the numberof transducers participated in Greens function measuring.Then,we demostrate that high-quality PA images can be achievedwith four scattering lenses,while PA image quality can
9、be improved by controlling the signal length.Scattering lens proposedin this paper needs not precision manufacture or special material,and just four transducers can accomplish quick imaging,withwhich PAT system cost can be reduced.Key words:photoacoustic tomography,low-cost,scattering lens,Greens fu
10、nction基金项目:国家自然科学基金(118 7 42 17),江苏省自然科学基金(BK20181077),江苏高校“青蓝工程”收稿日期:2 0 2 2-12-0 6*通讯联系人,E-mail:232南京大学学报(自然科学)第5 9 卷光声断层成像(PhotoacousticTomography,PAT)是一种基于光声效应的新型非侵人式生物医学成像方法,其具有光学成像对比度高以及声学成像可在深层组织获得高分辨率的特点 1-3 .PAT成像过程中,首先利用非聚焦脉冲激光束照射待测生物组织,生物组织中的光吸收体吸收了光波能量后由于热胀冷缩效应激发出超声波并向周边传播;随后,布置在组织周边的超声换
11、能器接收超声波并将其转化为电信号传输至上位机,上位机通过特定的重构算法形成生物组织内部光吸收系数分布图像.PAT本质上是光学成像,但光学参量的载体是超声波,它打破了光学散射对成像深度的约束,可以在5 6 cm的成像深度获得良好的声学空间分辨率。近年来,PAT在诸多医学领域展现出了广阔的应用前景,如:肿瘤诊断 4-5 、血管疾病诊断 6 、小动物全身成像 7-8 、关节炎诊断 9、药物输运监测 10 以及其他各类医学应用 11-13 早期的PAT系统大多采用单换能器单通道的硬件配置 14,利用步进电机带动换能器沿预先设定的轨道移动,根据合成孔径原理达到阵列成像的效果,这种方式存在鲁棒性较差,扫描
12、时间过长等缺点,难以适应当今临床医学领域对快速成像的需求.近年来,国内外研究机构针对多通道PAT系统开展了相关研究并取得了系列进展 7-8 .然而,多通道系统成本高昂的缺点始终是制约PAT大范围应用的重要因素.为减少换能器及通道数量,国内外学者提出了基于声学透镜的PAT15-16、基于声学腔穴的PAT17等方法,但这些方法亦存在器件需要特殊材质制作、加工难度大、鲁棒性较差等问题.为解决上述问题,本文提出了一种基于声学散射透镜的低成本光声断层成像方法,散射透镜由若干尺寸相同的圆形钢柱构成,单个信号探测换能器固定在散射透镜一侧,成像区位于散射透镜的另一侧,在光声成像的过程中,由于随机散射过程的多路
13、径效应,散射透镜起到了增大换能器等效数值孔径的作用,成像区至换能器的格林函数由超声干涉测量法测得,基于四个散射透镜全包围光吸收体共同作用,利用时间反转算法,可实现高质量的光声断层成像.所提出的散射透镜无须复杂制作工艺或特殊材质,只需四个超声换能器即可实现信号快速探测及成像,进而降低系统硬件成本1理论分析首先考虑如图1所示场景,为避免繁重的数值计算过程,本文采用了2 D模型.多路径效应的产生需要声波的多重散射,因而散射透镜要满足下列条件:(1)散射透镜的厚度要大于声波平均自由程;(2)散射体的尺寸要与声波波长接近.散射透镜核心元件为等距分布的钢制圆柱体,单个超声换能器位于散射透镜左侧ra处,该换
14、能器与散射透镜的相对位置关系保持固定不变.在进行光声成像前,需要在散射透镜右侧放置一个超声换能器阵列,用于测量散射透镜的格林函数.光吸收体位于散射透镜右侧的成像区内,在脉冲激光束的照射下,光吸收体吸收光波能量,并向外辐射光声信号,光声信号的初始声压值与光吸收系数成比例关系,即:q(r,t)A(r)aH(t)/at(1)A(r)为成像区内任一点r处的光吸收系数,H(t)为光脉冲的时域波形,q(r,t)为源,c为声速,光声信号的波动方程为:1 0p(r,t)Vp(r,t)q(r,t)(2)at2部分光声信号经过散射透镜并被位于r处的换能器接收,接收到的光声信号可表示为:pa(rd,t)dtdrq(
15、r,t)G(ra,r,t)(3)R其中,T为信号长度,R为被探测器包围的成像区域,G(rd,r,t)为位置ra与r之间的格林函数.将探测器接收到的信号在时域上反转并发射,反转声场将自适应地汇聚在源点处 18 ,即:p.(r,t)=dtdrpa(ra,t)G(ra,r,t)(4)0R式(4)在均匀介质中可演变为延迟叠加算法(Delay-and-Sum,DAS),通常情况下,利用式(4)重构光声信号初始声压需要换能器阵列分布的曲面包围R,然而,散射透镜的存在增大了探测器的等效数值孔径,部分原本无法被接收到的信号成分由于多路径效应,经多次散射后被换能器接收,使得利用较少换能器进行光声断层成像成为可2
16、33何君君等:基于声学散射透镜的低成本光声断层成像方法第2 期能 19.由式(4)不难看出,成像实现的重点在于格林函数G(rd,r,t)的确定.格林函数的确定可以使用超声干涉测量法,其利用超声信号相关性分析,消除两组信号中共同路径的部分 2 0 .如图1所示,在散射透镜右侧放置一组超声换能器阵列,阵列中的各单元所在位置为rs1,rs2,r s n,散射透镜左侧ra处的换能器向散射透镜发射一个波形为s(t)宽带超声信号,散射透镜右侧换能器阵列各单元接收到的超声信号依次记为u(t),uz(t),,u(t),此时,成像区内某点r至超声换能器r.的格林函数可利用下式近似求得:G(ra,r,w)J_F(
17、a)U(w)G(r,rs,)ds(5)其中 G(ra,r,),G(r,rs,w),U()分别是G(rd,r,t),G(r,rst,t),u,(t)的傅里叶变换形式,i=1,2,n;F(w)=S*(w)为波形滤波器,用于在格林函数重构过程中消除s(t)波形对计算的影响.在实际测量过程中,式(5)可离散化为:G(ra,r,a)ZF(a)U()G(r,rs,)(6)1由于模型中成像区和换能器阵列之间的介质为无损均匀介质,故G(r,rs,t)=(7)4元Si在地震干涉测量学等场合,格林函数的重构需要r1,rs,,rs所在的曲面包围ra及r并对其进行各向同性照射 2 0 .由于应用场景不同,本研究采用单
18、侧照射方式,由相关理论可知,在公式(6)中对于格林函数重构起主要作用的叠加信号来自菲涅尔区内的换能器单元 2 0 1,菲涅尔区如图1绿色虚线所示.然而,由于散射透镜的散射效应,必然会造成信号有效相干叠加区发生变化,需要通过定量研究提高格林函数的计算精度.此外,参与叠加的信号数量也与格林函数计算精度相关,故完成一次测量后可将换能器阵列移位,重复一次测量过程,得到一组新的信号ui(t),u2(t),,u(t),并与第一次测量结果共同代入公式(6)中计算,以提高计算精度。u(t)Mm日u(t)有效叠加区Yd测试点换能器成像区移位u(t)Hui(t)Time(us)散射透镜换能器阵列图1散射透镜结构及
19、格林函数测量过程示意图Fig.1Structure of the scattering lens and the illustra-tion for Greens function measuring2数值模拟实验2.1散射透镜的格林函数计算本文利用COMSOL软件建立了相应的数值模型,如图1所示,散射透镜为矩形,中心位置为(一7 mm,0 m m),尺寸设置为3 mm12mm.散射体为直径0.6mm的圆柱,边界条件为硬边界,相邻散射体中心间距为1mm,散射体分布密度为0.9 2 mm-2.理论上散射体排布方式可采用随机分布或规则分布,只要产生足够的多重散射即可,但在随机分布情况下,有可能出现
20、某一区域散射体分布过于密集或过于稀疏的情况,导致多路径效应只集中在散射透镜的某一部分发生,为避免这一情况,根据经验分析,本文中散射体采用上述布置方式。一个超声换能器位于r(一11mm,0mm)处,一个含有2 3 个等间距分布单元的超声换能器线型阵列位于散射透镜右侧,其轴坐标为4mm,y 轴坐标分布范围为一5.5 5.5 mm.ra处的超声换能器向散射透镜发射中心频率为3 MHz的高斯脉冲信号:s(t)=2元 f8(1-1/f.)exp(一 元 fe(t 1/f)(8)s(t)波形如图1插图所示.换能器阵列中各单元接收到的信号即为u(t),u(t),u(t),波形如图2 a所示,换能器阵列接收到
21、的波阵面如图2b所示,从图中可以看出光声信号在时域上存在多次振荡,这是由于散射透镜对光声信号造成了多重散射,多重散射的出现满足了多路径效应的形成条件.234南京大学学报(自然科学)第5 9 卷1.0(a)u(t)0.50.0-0.5-1.051015202530Time(us)(b)Max2015105Min51015202530Time(s)图2(a)经过散射后的超声脉冲波形;(b)换能器阵列探测到的波阵面Fig.2(a)The US pulse scattered by scattering lens,(b)wavefront picked upby the transducerarray
22、为了研究格林函数测量的精确性,我们在成像区中取一点测试点r(1m m,0 m m),如图1中红点所示.令其发射中心频率为3 MHz波形为s(t)的超声信号,该信号被ra处的超声换能器接收,表示为u(t),可知u(t)=G(rd,rt,t)?s(t),其中G(r d,r t,t)表示r与ra之间实际的格林函数.随后,做r与ra的连线延长线,并与换能器阵列相交,以交点为中心,依次增加交点两侧参与格林函数计算的换能器单元数量,在不同换能器数目的情况下利用公式(6)计算格林函数G(ra,r,t),并令us(t)=G(rd,It,t)s(t),可通过比较us(t)与u(t)波形相似程度确定格林函数计算的
23、精确性,本文选用相关系数比较二者相似程度,相关系数定义为:Couu,(r),u.(r)Ru,(r),u(r)=(9)Varu.(r)Varu(r)使用交点处及左右邻近共5 个换能器单元(=36.87大致为菲涅尔区)数据计算得到的格林函数所对应的u(t)及ut(t)波形如图3 a所示,此时相关系数R=0.839;使用11个换能器数据(=7 9.6 1)计算得到的格林函数所对应的u(t)及u(t)波形如图3 b所示,此时相关系数R=0.935,使用的换能器数目与R值之间的定量关系如图3 c所示.可见,由于散射透镜的出现,对格林函数重构起作用的换能器范围超出了传统的菲涅尔区,信号有效相干叠加区域与重
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 声学 散射 透镜 低成本 断层 成像 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。