搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于聚焦声场模型的光声层析成像时间延迟快速校正反投影方法

罗晓飞 王波 彭宽 肖嘉莹

引用本文:
Citation:

基于聚焦声场模型的光声层析成像时间延迟快速校正反投影方法

罗晓飞, 王波, 彭宽, 肖嘉莹

Back-projection method with fast time-delay correction for photoacoustic tomography reconstruction based on a focused sound field model

Luo Xiao-Fei, Wang Bo, Peng Kuan, Xiao Jia-Ying
PDF
HTML
导出引用
  • 光声层析成像是一种发展迅速的成像技术, 其可提供生物组织的结构和功能信息, 结合了光学成像高光学对比度与声学成像高穿透深度的优点. 然而, 由于现有的反投影成像算法通常将围绕目标扫描的超声换能器等效为一个点探测器, 导致非中心成像区域图像的切向模糊, 严重影响了图像质量. 本文提出一种新的光声层析成像算法, 其采用聚焦声场等效模型, 可以快速有效地克服换能器孔径效应所造成的声场畸变, 恢复非中心成像区域的切向分辨率. 仿真结果表明, 该方法对直径5 mm, 距离旋转中心6 mm的目标, 切向分辨率提升至少达2倍. 实验结果表明, 该方法可以有效地恢复边缘图像的切向模糊, 使得复杂目标的微小结构能被清晰探测. 这种新方法为传统的反投影方法提供了一种有价值的替代选择, 对基于圆/球扫描的光声层析成像系统的设计具有重要的指导作用.
    Photoacoustic tomography is a rapidly developing imaging technology, which can provide structural and functional information of biological tissues. It can integrate the advantage of high optical contrast from optical imaging with those from high penetration depth from ultrasound imaging. However, the existing back-projection algorithm in photoacoustic tomography equates the ultrasonic transducer scanning around the target to a point detector, which leads to a notable tangential blur in the eccentric imaging regions, and thus seriously degrading the image quality. In this paper, we propose a novel photoacoustic tomography reconstruction algorithm, which employs a focused sound field equivalent model to overcome the sound field distortion caused by the transducer’ finite aperture effect and can quickly and effectively restore the elongated tangential resolution in the eccentric imaging regions. Simulation results show that for the target with a diameter of 5 mm and a distance of 6 mm from the rotation center, the tangential resolution is improved by at least twice. Experimental results show that this method can effectively restore the image tangential blur in the off-center regions, where the tiny structures of complex targets can be detected. This new method provides a valuable alternative to the conventional back-projection method and plays an important guiding role in the design of photoacoustic tomography systems based on circle/sphere scanning.
      通信作者: 肖嘉莹, jiayingxiao@csu.edu.cn
    • 基金项目: 湖南省应急厅科技创新专项(批准号: 2021_QYC_10050_26366)、湖南省高新技术产业科技创新引领计划(批准号: 2020SK2003)和中南大学创新驱动项目团队项目(批准号: 2020CX004)资助的课题
      Corresponding author: Xiao Jia-Ying, jiayingxiao@csu.edu.cn
    • Funds: Project supported by the Emergency Science and Technology Project of Hunan Province (Grant No. 2021_QYC_10050_26366), the Department of Science and Technology of Hunan Province, High-tech Industry Science and Technology Innovation Leading Program (Grant No. 2020SK2003), Central South University, Innovation Driven Program Team Project (Grant No. 2020CX004).
    [1]

    Wang L V, Hu S 2012 Science 335 1458Google Scholar

    [2]

    Wang L V, Yao J J 2016 Nat. Methods 13 627Google Scholar

    [3]

    Zhou Y Y, Liang S Y, Li M S, Liu C B, Lai P X, Wang L D 2020 J. Biophotonics 13 e201960229

    [4]

    Zhou Y Y, Cao F, Li H H, Huang X Z, Wei D S, Wang L D, Lai P X 2020 Biomed. Opt. Express 11 2394Google Scholar

    [5]

    Liu Y Y, Fan X, Li Y Y, Qu L L, Qin H Y, Cao Y N, Li H T 2015 Prog. Chem. 27 1459

    [6]

    Yang X M, Maurudis A, Gamelin J, Aguirre A, Zhu Q, Wang L V 2009 J. Biomed. Opt. 14 71770K

    [7]

    Xia J, Wang L H V 2014 Ieee T. Bio-Med. Eng. 61 1380Google Scholar

    [8]

    Menke J 2015 Eur. Radiol. 25 2205Google Scholar

    [9]

    Xiao J Y, Yao L, Sun Y, Sobel E S, He J S, Jiang H B 2010 Opt. Express 18 14359Google Scholar

    [10]

    Wang B, Xiao J Y, Jiang H B 2014 J Neural. Eng. 11 046013Google Scholar

    [11]

    Li M C, Liu C B, Gong X J, Zheng R Q, Bai Y Y, Xing M Y, Du X M, Liu X Y, Zeng J, Lin R Q, Zhou H C, Wang S J, Lu G M, Zhu W, Fang C H, Song L 2018 Biomed. Opt. Express 9 1408Google Scholar

    [12]

    Zhou Y, Yao J J, Wang L H V 2016 J. Biomed. Opt. 21 061007Google Scholar

    [13]

    Burgholzer P, Bauer-Marschallinger J, Grun H, Haltmeier M, Paltauf G 2007 Inverse Probl. 23 S65Google Scholar

    [14]

    Xu M H, Wang L H V 2005 Phys. Rev. E 71 016706Google Scholar

    [15]

    Roitner H, Haltmeier M, Nuster R, O'Leary D P, Berer T, Paltauf G, Grun H, Burgholzer P 2014 J. Biomed. Opt. 19 056011Google Scholar

    [16]

    Li M L, Tseng Y C, Cheng C C 2010 Opt. Express 18 26285Google Scholar

    [17]

    Xu M H, Wang L V 2003 Phys. Rev. E 67 056605Google Scholar

    [18]

    Li C H, Wang L V 2008 Appl. Phys. Lett. 93 033902Google Scholar

    [19]

    Drozdov G, Rosenthal A 2017 Ieee T. Med. Imaging. 36 301Google Scholar

    [20]

    Zhang C, Wang Y Y 2008 J. Opt. Soc. Am. A. 25 2436Google Scholar

    [21]

    Cai D, Li Z F, Li Y, Guo Z D, Chen S L 2017 Opt. Express 25 1421Google Scholar

    [22]

    Luo X F, Xiao J Y, Wang C C, Wang B 2021 Photonics 8 356Google Scholar

  • 图 1  基于二维环形扫描的PAT系统

    Fig. 1.  2D circular-scanning-based PAT system.

    图 2  不同重建模型的反投影算法示意图 (a)点状换能器模型; (b)无限大小换能器模型; (c)聚焦声场换能器模型

    Fig. 2.  Schematics of the back-projection algorithm of different reconstruction models: (a) The model for point-like ultrasonic transducer; (b) the model for ultrasonic transducer with infinite element size; (c) the model for ultrasonic transducer with focused sound field.

    图 3  不同重建模型下换能器中心频率变化的仿真结果 (a)—(d)当中心频率分别为1 MHz, 3 MHz, 5 MHz和10 MHz时, 用传统点模型方法得到的结果; (e)—(h)无限大平面模型的相应结果; (i)—(l)聚焦声场模型的相应结果. 所采用的平面换能器直径为6 mm, 7个点的位置分别为–6 mm, –4 mm, –2 mm, 0 mm, 2 mm, 4 mm和6 mm. 在模拟数据中加入5%的噪声, 所有图像的振幅都归一化为1

    Fig. 3.  Simulation results with different reconstruction methods when the transducer central frequency varies: (a)–(d) The results with the conventional point-like method when the central frequency are 1 MHz, 3 MHz, 5 MHz and 10 MHz, respectively; (e)–(h) the corresponding results with the model of the infinite element size; (i)–(l) the corresponding results with the model of focused sound field. The diameter of the employed planar transducer was 6 mm, and the locations of the seven points are –6 mm, –4 mm, –2 mm, 0 mm, 2 mm, 4 mm, and 6 mm, respectively. The simulated data were added with 5% noise, and the amplitudes of all the images were normalized to 1 for comparison.

    图 4  由3 MHz中心频率超声换能器模拟得到的4个目标的横向剖面, 红线、蓝线和绿线分别表示点状换能器模型、无限大小换能器模型和聚焦声场换能器模型的结果 (a) 0 mm; (b) 2 mm; (c) 4 mm; (d) 6 mm

    Fig. 4.  The simulated lateral profiles for the four targets obtained from 3 MHz central frequency ultrasonic transducer. The locations of the four points are 0 mm (a), 2 mm (b), 4 mm (c), and 6 mm (d), respectively. The red, blue and green lines represent the results from the conventional point-like model, the infinite element size model, and the focused sound field model, respectively.

    图 5  PAT系统和超声换能器的传输特性示意图 (a)二维PAT系统图; (b)超声换能器光声脉冲响应信号; (c)超声换能器光声频谱图

    Fig. 5.  Schematic of the PAT system and transmit characteristic of the transducer: (a) Schematic of the 2D PAT system; (b) the photoacoustic impulse waveform of the transducer; (c)spectral diagram of the transducer.

    图 6  不同重建算法对点目标成像的实验结果 (a)仿体照片; (b)8个点目标的切向分辨率分布图; (c)点状模型BP算法; (d)无限大小模型BP算法; (e)聚焦声场模型BP算法

    Fig. 6.  The experimental results of point target imaging with different reconstruction algorithms: (a) The photograph of phantom; (b) tangential resolution distributions of eight point targets; (c) the BP method of point-like model; (d) the BP method of infinite element size model; (e) the BP method of focused sound field model.

    图 7  不同重建算法对小鼠模型成像的实验结果 (a)小鼠脑图片; (b)点状模型BP算法; (c)无限大小模型BP算法; (d) 聚焦声场模型BP算法

    Fig. 7.  The experimental results of mouse brain imaging with different reconstruction algorithms: (a) The photograph of the mouse brain; (b) the BP method of point-like model; (c) the BP method of infinite element size model; (d) the BP method of focused sound field model.

    表 1  Table 1 不同算法和换能器频率下偏离中心6 mm目标的FWHM和SNR

    Table 1.  The FWHM and SNR results of 6 mm off-center target with different methods and transducer frequencies.

    中心频率点状无限大小聚焦声场
    1 MHz
    FWHM/mm1.851.050.85
    SNR/dB36.0626.1736.80
    3 MHz
    FWHM/mm1.251.040.42
    SNR/dB28.7827.2933.87
    5 MHz
    FWHM/mm0.930.820.38
    SNR/dB26.4224.0728.09
    10 MHz
    FWHM/mm0.450.410.25
    SNR/dB22.9016.7525.31
    下载: 导出CSV
  • [1]

    Wang L V, Hu S 2012 Science 335 1458Google Scholar

    [2]

    Wang L V, Yao J J 2016 Nat. Methods 13 627Google Scholar

    [3]

    Zhou Y Y, Liang S Y, Li M S, Liu C B, Lai P X, Wang L D 2020 J. Biophotonics 13 e201960229

    [4]

    Zhou Y Y, Cao F, Li H H, Huang X Z, Wei D S, Wang L D, Lai P X 2020 Biomed. Opt. Express 11 2394Google Scholar

    [5]

    Liu Y Y, Fan X, Li Y Y, Qu L L, Qin H Y, Cao Y N, Li H T 2015 Prog. Chem. 27 1459

    [6]

    Yang X M, Maurudis A, Gamelin J, Aguirre A, Zhu Q, Wang L V 2009 J. Biomed. Opt. 14 71770K

    [7]

    Xia J, Wang L H V 2014 Ieee T. Bio-Med. Eng. 61 1380Google Scholar

    [8]

    Menke J 2015 Eur. Radiol. 25 2205Google Scholar

    [9]

    Xiao J Y, Yao L, Sun Y, Sobel E S, He J S, Jiang H B 2010 Opt. Express 18 14359Google Scholar

    [10]

    Wang B, Xiao J Y, Jiang H B 2014 J Neural. Eng. 11 046013Google Scholar

    [11]

    Li M C, Liu C B, Gong X J, Zheng R Q, Bai Y Y, Xing M Y, Du X M, Liu X Y, Zeng J, Lin R Q, Zhou H C, Wang S J, Lu G M, Zhu W, Fang C H, Song L 2018 Biomed. Opt. Express 9 1408Google Scholar

    [12]

    Zhou Y, Yao J J, Wang L H V 2016 J. Biomed. Opt. 21 061007Google Scholar

    [13]

    Burgholzer P, Bauer-Marschallinger J, Grun H, Haltmeier M, Paltauf G 2007 Inverse Probl. 23 S65Google Scholar

    [14]

    Xu M H, Wang L H V 2005 Phys. Rev. E 71 016706Google Scholar

    [15]

    Roitner H, Haltmeier M, Nuster R, O'Leary D P, Berer T, Paltauf G, Grun H, Burgholzer P 2014 J. Biomed. Opt. 19 056011Google Scholar

    [16]

    Li M L, Tseng Y C, Cheng C C 2010 Opt. Express 18 26285Google Scholar

    [17]

    Xu M H, Wang L V 2003 Phys. Rev. E 67 056605Google Scholar

    [18]

    Li C H, Wang L V 2008 Appl. Phys. Lett. 93 033902Google Scholar

    [19]

    Drozdov G, Rosenthal A 2017 Ieee T. Med. Imaging. 36 301Google Scholar

    [20]

    Zhang C, Wang Y Y 2008 J. Opt. Soc. Am. A. 25 2436Google Scholar

    [21]

    Cai D, Li Z F, Li Y, Guo Z D, Chen S L 2017 Opt. Express 25 1421Google Scholar

    [22]

    Luo X F, Xiao J Y, Wang C C, Wang B 2021 Photonics 8 356Google Scholar

  • [1] 向鹏程, 蔡聪波, 王杰超, 蔡淑惠, 陈忠. 基于深度神经网络的时空编码磁共振成像超分辨率重建方法. 物理学报, 2022, 71(5): 058702. doi: 10.7498/aps.71.20211754
    [2] 郁钧瑾, 郭星奕, 隋怡晖, 宋剑平, 他得安, 梅永丰, 许凯亮. 超分辨率超快超声脊髓微血管成像方法. 物理学报, 2022, 71(17): 174302. doi: 10.7498/aps.71.20220629
    [3] 隋怡晖, 郭星奕, 郁钧瑾, Alexander A. Solovev, 他得安, 许凯亮. 生成对抗网络加速超分辨率超声定位显微成像方法研究. 物理学报, 2022, 71(22): 224301. doi: 10.7498/aps.71.20220954
    [4] 吕浩昌, 赵云驰, 杨光, 董博闻, 祁杰, 张静言, 朱照照, 孙阳, 于广华, 姜勇, 魏红祥, 王晶, 陆俊, 王志宏, 蔡建旺, 沈保根, 杨峰, 张申金, 王守国. 基于深紫外激光-光发射电子显微技术的高分辨率磁畴成像. 物理学报, 2020, 69(9): 096801. doi: 10.7498/aps.69.20200083
    [5] 高强, 李小秋, 周志鹏, 孙磊. 基于分形谐振器的远场超分辨率扫描成像. 物理学报, 2019, 68(24): 244102. doi: 10.7498/aps.68.20190620
    [6] 席晓琦, 韩玉, 李磊, 闫镔. 螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法. 物理学报, 2019, 68(8): 088701. doi: 10.7498/aps.68.20190055
    [7] 高强, 王晓华, 王秉中. 基于宽带立体超透镜的远场超分辨率成像. 物理学报, 2018, 67(9): 094101. doi: 10.7498/aps.67.20172608
    [8] 唐弢, 赵晨, 陈志彦, 李鹏, 丁志华. 超高分辨光学相干层析成像技术与材料检测应用. 物理学报, 2015, 64(17): 174201. doi: 10.7498/aps.64.174201
    [9] 李杰, 朱京平, 张云尧, 刘宏, 侯洵. 光谱分辨率可调的新型干涉成像光谱技术研究. 物理学报, 2013, 62(2): 024205. doi: 10.7498/aps.62.024205
    [10] 张峰, 闫镔, 汪先超, 江桦, 魏星. 半覆盖锥束CT中扁平物体的高效反投影滤波重建. 物理学报, 2013, 62(16): 168702. doi: 10.7498/aps.62.168702
    [11] 张文喜, 相里斌, 孔新新, 李杨, 伍洲, 周志盛. 相干场成像技术分辨率研究. 物理学报, 2013, 62(16): 164203. doi: 10.7498/aps.62.164203
    [12] 王芳, 赵星, 杨勇, 方志良, 袁小聪. 基于人眼视觉的集成成像三维显示分辨率的比较. 物理学报, 2012, 61(8): 084212. doi: 10.7498/aps.61.084212
    [13] 卢婧, 李昊, 何毅, 史国华, 张雨东. 超分辨率活体人眼视网膜共焦扫描成像系统. 物理学报, 2011, 60(3): 034207. doi: 10.7498/aps.60.034207
    [14] 代秋声, 漆玉金. 针孔单光子发射计算机断层成像的空间分辨率研究. 物理学报, 2010, 59(2): 1357-1365. doi: 10.7498/aps.59.1357
    [15] 张小正, 毕传兴, 徐亮, 陈心昭. 基于波叠加法的近场声全息空间分辨率增强方法. 物理学报, 2010, 59(8): 5564-5571. doi: 10.7498/aps.59.5564
    [16] 吴丹, 陶超, 刘晓峻. 有限方位扫描的光声断层成像分辨率研究. 物理学报, 2010, 59(8): 5845-5850. doi: 10.7498/aps.59.5845
    [17] 赵贵敏, 陆明珠, 万明习, 方莉. 高分辨率扇形阵列超声激发振动声成像研究. 物理学报, 2009, 58(9): 6596-6603. doi: 10.7498/aps.58.6596
    [18] 向良忠, 邢达, 郭华, 杨思华. 高分辨率快速数字化光声CT乳腺肿瘤成像. 物理学报, 2009, 58(7): 4610-4617. doi: 10.7498/aps.58.4610
    [19] 陈湛旭, 唐志列, 万 巍, 何永恒. 基于声透镜成像系统的光声层析成像. 物理学报, 2006, 55(8): 4365-4370. doi: 10.7498/aps.55.4365
    [20] 刘 力, 邓小强, 王桂英, 徐至展. 改善共焦系统轴向分辨率的位相型光瞳滤波器. 物理学报, 2001, 50(1): 48-51. doi: 10.7498/aps.50.48
计量
  • 文章访问数:  4649
  • PDF下载量:  67
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-10-31
  • 修回日期:  2021-12-06
  • 上网日期:  2022-01-26
  • 刊出日期:  2022-04-05

/

返回文章
返回