复频域实时重建医学光声图像的研究

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

复频域实时重建医学光声图像的研究

余双春，袁　杰*

出版日期:2016-07-03 发布日期:2016-07-03
作者简介: 南京大学电子科学与工程学院，南京，210093
基金资助:
国家自然科学基金(61201425)，江苏省自然科学基金(BK20131280)

Realtime reconstruction of the medical complex frequency domain photoacoustic image

Yu Shuangchun，Yuan Jie*

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

摘要/Abstract

摘要： 提出将光声成像中的傅里叶域重建方法与分布式并行计算相结合的方法，应用于实际处理医学成像中，实现实时重建医学光声图像的目的．滤波反投影算法是光声图像重建的一个广泛使用的算法，但由于在时域重建的时间复杂度大，并不适合实时重建．傅里叶域重建方法相对于滤波反投影算法具有时间复杂度上的优越性，且产生较少的伪影．并行处理技术，已经广泛应用于各类工程学科，其影响力及范围随着研究深入逐步扩大．分布式并行计算将提高图像处理速度，使实时重建医学图像成为可能．将MATLAB作为平台，实现了对医学图像在复频域的并行加速处理，给出了并行处理与串行处理的测试结果的比较．实验结果表明：傅里叶域重建方法与分布式并行计算的结合，使医学光声成像重建时间大大缩短．

Abstract: In recent years，photoacoustic imaging has developed gradually as a new medical imaging modality.Traditional medical imaging methods，such as Xray Tomography，Optical Coherence Tomography，MRI and ultrasound imaging，cannot meet the requirement of nondestructive and nonionizing method，as well as deeper imaging depth and higher imaging contrast of organ.Speaking of the advantages of photoacoustic effect in biological imaging，it is noninvasive，nonionizing，and easy to perform.It has great prospects for development in the field of medical diagnosis and treatment.The filtered back projection algorithm is a widely used algorithm in photoacoustic image reconstruction.However，because of the time complexity in the process of reconstruction，filtered back projection algorithm is not suitable for realtime reconstruction.In order to solve the above problem，this article is a combination of photoacoustic imaging Fourier domain reconstruction and distributed parallel computing to achieve the purpose of realtime medical photoacoustic image reconstruction.Comparing to the widely used filtered backprojection algorithm，Fourier domain reconstruction methods is advanced in time complexity and has fewer artifacts.The parallel processing technology has also been widely used in various engineering disciplines.Its influence and scope of application have been expanded with further research，emerging in many different areas other than scientific computing.Distributed parallel computing makes realtime reconstruction of medical images possible by accelerating and improving the image processing speed.In this paper，MATLAB is the platform of accelerating the speed of the medical image processing in complex frequency domain.The test results and comparison of parallel processing and serial processing are also presented.The result of experiment indicates that by using Fourier domain reconstruction and distributed parallel computing，the reconstruction time is greatly shortened and the efficiency is improved.

余双春，袁　杰*. 复频域实时重建医学光声图像的研究[J]. 南京大学学报(自然科学版), 2016, 52(3): 536–541.

Yu Shuangchun，Yuan Jie*. Realtime reconstruction of the medical complex frequency domain photoacoustic image[J]. Journal of Nanjing University(Natural Sciences), 2016, 52(3): 536–541.

参考文献

[1]　Hu S，Wang L V．Photoacoustic imaging and characterization of the microvasculature.Journal of biomedical optics，2010，15(1)：011101－011101－15.
[2] Wang L V．Multiscale photoacoustic microscopy and computed tomography.Nature photonics，2009，3(9)：503－509.
[3] Wang L V．Prospects of photoacoustic tomography.Medical Physics，2008，35(12)：5758－5767.
[4] Webb A，Kagadis G C．Introduction to biomedical imaging.Hoboken：Wiley，2003，30(8)：2267.
[5] Wang X，Roberts W W，Carson P L，et al.Photoacoustic tomography：A potential new tool for prostate cancer.Biomedical Optics Express，2010，1(4)：1117－1126.
[6] 陈清明，程祖海，朱海红．脉冲激光在水中激发声脉冲的光声能量转换效率．中国激光，2007，3：341－344.(Chen Q M，Cheng Z H，Zhu H H．Laser acoustic energy conversion efficiency induced by pulse laser in water.Chinese Journal of Lasers，2007，3：341－344.)
[7] Zhang W，Wu Z Y，Yu Q X．Photoacoustic spectroscopy for fast and sensitive ammonia detection.Chinese Optics Letters，2007，11：677－679.
[8] Robert A K，William L K，Daniel R R，et al.Application of thermoacoustic computed tomography to breast imaging.Proc Spie(The International Society for Optical Engineering)，1999，3659：426－430.
[9] Xu M，Wang L V．Pulsedmicrowareinduced thermoacoustic tomography：Filtered backprojection in a circular measurement configuration.Medical Physics，2002，29(8)：1661－1669.
[10] 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.
[11] Köstli K P，Frenz M，Bebie H，et al.Temporal backward projection of optoacoustic pressure transients using Fourier transform methods.Physics in Medicine and Biology，2001，46(7)：1863.
[12] Wang K，Anastasio M A．A simple fourier transformbased reconstruction formula for photoacoustic computed tomography with a circular or spherical measurement geometry.Physics in Medicine and Biology，2012，57(23)：N493.
[13] Wang L．Parallel computing technology review.Information Technology，2012，(10)：112－115.
[14] 沈绪榜，张发存，冯国臣等．计算机体系结构的分类模型．计算机学报，2005，28(11)：1760－1762.(Shen X B，Zhang F C，Feng G C，et al.The classification model of computer architectures.Chinese Journal of Computers，2005，28(11)：1760－1762.)
[15] 卢风顺，宋君强，银福康等．CPU/GPU协同并行计算研究综述．计算机科学，2011，38(3)：5－9.(Lu F S，Song J Q，Yin F K，et al.Survey of CPU/GPU synergetic parallel computing.Computer Science，2011，38(3)：5－9.)
[16] Treeby B E，Varslot T K，Zhang E Z，et al.Automatic sound speed selection in photoacoustic image reconstruction using an autofocus approach.Journal of Biomedical Optics，2011，16(9)：090501－090501－3.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed