一种基于光声成像选取最佳声速组的研究

南京大学学报(自然科学版) ›› 2016, Vol. 52 ›› Issue (3): 528–535.

一种基于光声成像选取最佳声速组的研究

叶　濛，袁　杰*

出版日期:2016-07-03 发布日期:2016-07-03
作者简介: 南京大学电子科学与工程学院，南京，210093

The study of selecting the optimum SOS group based on PA imaging

Ye Meng，Yuan Jie*

Online:2016-07-03 Published:2016-07-03
About author:School of Electronic Science and Engineering，Nanjing University，Nanjing，210093，China

摘要/Abstract

摘要： 光声成像是一种新兴的医学影像技术，具有图像对比度高、分辨率高、传递信息量大的特点，在生物医学领域有着广泛的应用．为了达到更好的重建效果，需要补偿声速在生物组织中传播的不均匀性．光声成像的时域反投影重建算法能够改变系统声速以及修正接收光声信号的时间，使得重建图像的质量与传播声速的取值密切相关，从而在修改声速的同时得到相应的成像结果．使用时域反投影算法重建图像，根据更新的声速组的搜索范围和搜索步长，采用不同的图像聚焦程度的评价方法，自适应地选取最佳声速组．研究中，对选取最佳声速的过程进行了多组实验，通过逐渐地缩小声速组的搜索范围和搜索步长，根据相应聚焦程度的评价方法，选取出最佳声速组．最终该组的声速值会收敛于正确的声速值附近，实验结果中的最大误差率都小于0.5%.

Abstract: Photoacoustic(PA)imaging and sensing，an emerging hybrid technology involving both light and sound，has excellent sensitivity to the chemical and physiological information in biological samples and has been explored for potential application in the diagnosis osteoporosis，breast cancer，arthritis and relevant disease.Research on photoacoustic tomography(PAT)got prosperous development for its being promisingly characterized with noninvasive and nonionizing.PA imaging has a wide range of applications in the biomedical field，which with the characteristics of good contrast，great resolution and transmitting of vast information amount.In order to achieve better reconstruction effect，we need to compensate the speed of inhomogeneous sound propagation in biological tissue.This paper analyzes the basic PA imaging reconstruction algorithm and indicates the superiority of the timedomain back projection algorithm，which not only changes system speed of sound(SOS)but also corrects the transmission time of acoustic pressure signal.The SOS in biomedical tissue typically used 1540 m?s－1，despite the fact that in soft tissues this property is known to range from approximately 1450 m?s－1 in fat to over 1600 m?s－1 in muscle.These 5% SOS errors will significantly reduce the contrast and spatial resolution of PAT image.The timedomain back projection reconstruction algorithm is better than frequencydomain reconstruction algorithm and delayandsum method.According to the updated searching range and steplength，we can adopt different degree assessment criteria of focus to choose the optimum SOS group adaptively.Through diminishing the searching range and steplength of the group，we make the choice procedure to be convergence.The performance of the choice the optimum SOS group procedure is illustrated by simulation.Furthermore，the performances of the selected speeds are also compared in the maximum error ratio to the correct speed group.Simulation and experiments verify the choice procedure method and the correct SOS so as to build better focused PAT images.

叶　濛，袁　杰*. 一种基于光声成像选取最佳声速组的研究[J]. 南京大学学报(自然科学版), 2016, 52(3): 528–535.

Ye Meng，Yuan Jie*. The study of selecting the optimum SOS group based on PA imaging[J]. Journal of Nanjing University(Natural Sciences), 2016, 52(3): 528–535.

参考文献

[1] Xu M, Wang L V. Photoacoustic imaging in biomedicine. Review of Scientific Instruments, 2006, 77(4): 041101.
[2] 张弛. 光声成像的图像重建算法研究. 硕士学位论文. 上海: 复旦大学, 2009. (Zhang C. The study of image reconstruction algorithm based on photoacoustic imaging. Master Dissertation. Shanhai: Fudan University, 2009. )
[3] 宋智源. 光声成像技术的重建算法与实验研究. 硕士学位论文. 天津: 天津大学, 2009. (Song Z Y. Research on the reconstruction algorithm and experiments of photoacoustic imaging technology. Master Dissertation. Tianjin: Tianjin University, 2009. )
[4] 张弛, 汪源源. 声速不均匀介质的光声成像重建算法. 光学学报, 2008, 28(12): 2296－2301. (Zhang C, Wang Y Y. Photoacoustic imaging reconstruction algorithm for speed of sound inhomogeneous medium. Journal of Optics, 2008, 28(12): 2296－2301. )
[5] Kruger R A, Liu P, Appledorn C R. Photoacoustic ultrasound(PAUS)reconstruction tomography. Medical Physics, 1995, 22(10): 1605－1609.
[6] 杨宁利, 张淑仪, 孙利等. 脉冲光声量热法研究氧合血红蛋白的光解反应. 南京大学学报(自然科学), 2007, 43(1): 73－78. (Yang N L, Zhang S Y, Sun L, et al. Photodissociation of oxyhemoglobins by pulsed photoacoustic calorimetry. Journal of Nanjing University(Natural Sciences), 2007, 43(1): 73－78. )
[7] Xu M, Wang L V. Timedomain reconstruction for thermoacoustic tomography in a spherical geometry. IEEE Transactions on Medical Imaging, 2002, 21(7): 814－822.
[8] Köstli K P, Beard P C. Twodimensional photoacoustic imaging by use of Fouriertransform image reconstruction and a detector with an anisotropic response. Applied Optics, 2003, 42(10): 1899－1908.
[9] Anastasio M A, Zhang J, Pan X, et al. Halftime image reconstruction in thermoacoustic tomography. IEEE Transactions on Medical Imaging, 2005, 24(2): 199－210.
[10] Haltmeier M, Schuster T, Scherzer O. Filtered backprojection for thermoacoustic computed tomography in spherical geometry. Mathematical Methods in the Applied Sciences, 2005, 28(16): 1919－1937.
[11] Kunyansky L A. Explicit inversion formulae for the spherical mean Radon transform. Inverse Problems, 2007, 23(1): 373.
[12] 杜恒, 刘东权. 超声系统声速优化模拟与算法. 中国生物医学工程进展—2007中国生物医学工程联合学术年会论文集(上册). 北京: 中国生物医学工程学会, 2007, 656－659. (Du H, Liu D Q. Optimization and algorithm of sound velocity in ultrasonic system. In: Progress in Chinese Society of Biomedical Engineering－2007 Proceedings of the annual conference of the Chinese Society of Biomedical Engineering Association(1). Beijing: Beijing Chinese Society of Biomedical Engineering Institute, 2007, 656－659. )
[13] He H, Liu D C. Sound speed optimization based on acoustic point spread function. In: The 3rd International Conference on Bioinformatics and Biomedical Engineering. Beijing: IEEE, 2009, 1－4.
[14] 李晓瑛, 刘东权. 基于图像对比度分析的自适应超声声速优化. 生物医学工程学杂志, 2012, 28(6): 1094－1097. (Li X Y, Liu D Q. An adaptive ultrasound sound speed optimization based on image contrast analysis. Journal of Biomedical Engineering, 2012, 28(6): 1094－1097. )
[15] Cho M H, Kang L H, Kim J S, et al. An efficient sound speed estimation method to enhance image resolution in ultrasound imaging. Ultrasonics, 2009, 49(8): 774－778.
[16] Sakaguchi T, Hirano T, Watanabe Y, et al. Inner head acoustic field for boneconducted sound calculated by finitedifference timedomain method. Japanese Journal of Applied Physics, 2002, 41(5S): 3604.
[17] Riad S M. The deconvolution problem: an overview. Proceedings of the IEEE, 1986, 74(1): 82－85.
[18] Bennia A, Riad S M. An optimization technique for iterative frequencydomain deconvolution. IEEE Transactions on Instrumentation and Measurement, 1990, 39(2): 358－362.
[19] 朱孔凤, 姜威, 王端芳等. 一种新的图像清晰度评价函数. 红外与激光工程, 2005, 34(4): 464－468. (Zhu K F, Jiang W, Wang D F, et al. New kind of clarityevaluation function of image. Infrared and Laser Engineering, 2005, 34(4): 464－468. )

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed