Double-stage delay-multiply-and-sum (DS-DMAS) is an algorithm proposed for photoacoustic image reconstruction. The DS-DMAS algorithm offers a higher contrast than conventional delay-and-sum and delay-multiply and-sum but at the expense of higher computational complexity. Here, we utilized a compute unified device architecture (CUDA) graphics processing unit (GPU) parallel computation approach to address the high complexity of the DS-DMAS for photoacoustic image reconstruction generated from a commercial light-emitting diode (LED)–based photoacoustic scanner. In comparison with a single-threaded central processing unit (CPU), the GPU approach increased speeds by nearly 140-fold for 1024 × 1024 pixel image; there was no decrease in accuracy. The proposed implementation makes it possible to reconstruct photoacoustic images with frame rates of 250, 125, and 83.3 when the images are 64 × 64, 128 × 128, and 256 × 256, respectively. Thus, DS-DMAS can be efficiently used in clinical devices when coupled with CUDA GPU parallel computation.
WangXPangYKuGXieXStoicaGWangLV. Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain. Nat Biotechnol. 2003;21(7):803-6.
2.
WangXXieXKuGWangLVStoicaG. Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. J Biomed Opt. 2006;11(2):024015.
3.
WangJLinCYMooreCJhunjhunwalaAJokerstJV. Switchable photoacoustic intensity of methylene blue via sodium dodecyl sulfate micellization. Langmuir. 2018;34(1):359-65.
4.
XuMWangLV. Photoacoustic imaging in biomedicine. Rev Sci Instrum. 2006;77(4):041101.
5.
JokerstJVThangarajMKempenPJSinclairRGambhirSS. Photoacoustic imaging of mesenchymal stem cells in living mice via silica-coated gold nanorods. ACS Nano. 2012;6(7):5920-30.
6.
KimTLemasterJEChenFLiJJokerstJV. Photoacoustic imaging of human mesenchymal stem cells labeled with Prussian blue-poly(l-lysine) nanocomplexes. ACS Nano. 2017;11(9):9022-32.
7.
WangLVHuS. Photoacoustic tomography: in vivo imaging from organelles to organs. Science. 2012;335(6075):1458-62.
8.
GuoBLiJZmudaHSheplakM. Multifrequency microwave-induced thermal acoustic imaging for breast cancer detection. IEEE Trans Biomed Eng. 2007;54(11):2000-10.
9.
PramanikMKuGLiCWangLV. Design and evaluation of a novel breast cancer detection system combining both thermoacoustic (TA) and photoacoustic (PA) tomography. Med Phys. 2008;35(6 Pt. 1):2218-23.
10.
MehrmohammadiMJoon YoonSYeagerDEmelianovSY. Photoacoustic imaging for cancer detection and staging. Curr Mol Imaging. 2013;2(1):89-105.
11.
HaririAWangJKimYJhunjhunwalaAChaoDLJokerstJV. In vivo photoacoustic imaging of chorioretinal oxygen gradients. J Biomed Opt. 2018;23(3):036005.
12.
de La ZerdaAPaulusYMTeedRBodapatiSDollbergYKhuri-YakubBTBlumenkranzMSMoshfeghiDMGambhirSS. Photoacoustic ocular imaging. Opt Lett. 2010;35(3):270-2.
13.
PuKShuhendlerAJJokerstJVMeiJGambhirSSBaoZRaoJ. Semiconducting polymer nanoparticles as photoacoustic molecular imaging probes in living mice. Nat Nanotechnol. 2014;9(3):233-9.
14.
TuraniZFatemizadehEBlumettiTDaveluySMoraesAFChenWMehreganDAndersenPENasiriavanakiM. Optical Radiomic Signatures Derived from Optical Coherence Tomography Images Improve Identification of Melanoma. Cancer Research. 2019;79(8):2021-2030.
15.
YaoJXiaJMaslovKINasiriavanakiMTsytsarevVDemchenkoAVWangLV. Noninvasive photoacoustic computed tomography of mouse brain metabolism in vivo. NeuroImage. 2013;64:257-66.
16.
NasiriavanakiMXiaJWanHBauerAQCulverJPWangLV. High-resolution photoacoustic tomography of resting-state functional connectivity in the mouse brain. Proc Natl Acad Sci U S A. 2014;111(1):21-6.
17.
LiMLOhJTXieXKuGWangWLiCLunguGStoicaGWangLV. Simultaneous molecular and hypoxia imaging of brain tumors in vivo using spectroscopic photoacoustic tomography. Proc IEEE. 2008;96(3):481-9.
18.
ValluruKSWilsonKEWillmannJK. Photoacoustic imaging in oncology: translational preclinical and early clinical experience. Radiology. 2016;280(2):332-49.
19.
SongWWeiQLiuWYiJSheibaniNFawziAALinsenmeierRAJiaoSZhangHF. A combined method to quantify the retinal metabolic rate of oxygen using photoacoustic ophthalmoscopy and optical coherence tomography. Sci Rep. 2014;4:6525.
20.
TaruttisAHerzogERazanskyDNtziachristosV. Real-time imaging of cardiovascular dynamics and circulating gold nanorods with multispectral optoacoustic tomography. Opt Express. 2010;18(19):19592-602.
21.
Van de SompelDSasportasLSJokerstJVGambhirSS. Comparison of deconvolution filters for photoacoustic tomography. PLoS ONE. 2016;11(3):e0152597.
22.
WangLV. Multiscale photoacoustic microscopy and computed tomography. Nat Photonics. 2009;3(9):503-9.
23.
MastandunoMAGambhirSS. Quantitative photoacoustic image reconstruction improves accuracy in deep tissue structures. Biomed Opt Express. 2016;7(10):3811-25.
XuMXuYWangLV. Time-domain reconstruction algorithms and numerical simulations for thermoacoustic tomography in various geometries. IEEE Trans Biomed Eng. 2003;50(9):1086-99.
27.
MozaffarzadehMPeriyasamyVPramanikMMakkiabadiB. An efficient nonlinear beamformer based on p^{th} root of detected signals for linear-array photoacoustic tomography: application to sentinel lymph node imaging, Arxiv:1805.09913, 2018. Available from https://arxiv.org/abs/1805.09913.
28.
UpputuriPKPramanikM. Recent advances toward preclinical and clinical translation of photoacoustic tomography: a review. J Biomed Opt. 2016;22(4):041006.
29.
SivasubramanianKPeriyasamyVDienzoRAPramanikM. Hand-held, clinical dual mode ultrasound—photoacoustic imaging of rat urinary bladder and its applications. J Biophotonics. 2018;11(5):e201700317.
30.
ZengYXingDWangYYinBChenQ. Photoacoustic and ultrasonic coimage with a linear transducer array. Opt Lett. 2004;29(15):1760-2.
31.
MercepEJengGMorscherSLiPCRazanskyD. Hybrid optoacoustic tomography and pulse-echo ultrasonography using concave arrays. IEEE T Ultrason Ferr. 2015;62(9):1651-61.
32.
YinBXingDWangYZengYTanYChenQ. Fast photoacoustic imaging system based on 320-element linear transducer array. Phys Med Biol. 2004;49(7):1339-46.
33.
ZhouQJiXXingD. Full-field 3D photoacoustic imaging based on plane transducer array and spatial phase-controlled algorithm. Med Phys. 2011;38(3):1561-6.
34.
MozaffarzadehMMahloojifarAOroojiMKratkiewiczKAdabiSNasiriavanakiM. Linear-array photoacoustic imaging using minimum variance-based delay multiply and sum adaptive beamforming algorithm. J Biomed Opt. 2018;23(2):026002.
35.
ParkSKarpioukABAglyamovSREmelianovSY. Adaptive beamforming for photoacoustic imaging using linear array transducer. In: Ultrasonics Symposium, Beijing, China, 2-5 November 2008-21 March2009, pp. 1088-91.
36.
AlshayaAHarputSMoubarkAMCowellDMJMcLaughlanJFreearS. Spatial resolution and contrast enhancement in photoacoustic imaging with filter delay multiply and sum beamforming technique. In: Ultrasonics Symposium, Tours, France, 18-21 September-3 November2016, pp. 1-4.
37.
MozaffarzadehMSadeghiMMahloojifarAOroojiM. Double-stage delay multiply and sum beamforming algorithm applied to ultrasound medical imaging. Ultrasound Med Biol. 2018;44(3):677-86.
38.
MozaffarzadehMMahloojifarAOroojiMAdabiSNasiriavanakiM. Double-stage delay multiply and sum beamforming algorithm: application to linear-array photoacoustic imaging. IEEE Trans Biomed Eng. 2018;65(1):31-42.
39.
ParidarRMozaffarzadehMMehrmohammadiMOroojiM. Photoacoustic image formation based on sparse regularization of minimum variance beamformer. Biomed Opt Express. 2018;9(6):2544-61.
40.
MatroneGSavoiaASCalianoGMagenesM. The delay multiply and sum beamforming algorithm in ultrasound B-mode medical imaging. IEEE Trans Med Imaging. 2015;34(4):940-9.
41.
MokdadAAzmiPMokariNMoltafetMGhaffari-MiabM. Cross-layer energy efficient resource allocation in PD-NOMA based H-CRANs: implementation via GPU. IEEE Trans Mob Computing. 2019;18:1246-59. doi:10.1109/TMC.2018.2860985.
42.
Masumnia-BishehKGhaffari-MiabMZakeriB. Evaluation of different approximations for correlation coefficients in stochastic FDTD to estimate SAR variance in a human head model. IEEE Trans Electromagn C. 2017;59(2):509-17. doi:10.1109/TEMC.2016.2614128.
43.
PratxGXingL. GPU computing in medical physics: a review. Med Phys. 2011;38(5):2685-97.
44.
RostamiSRMGhaffari-MiabM. Finite difference generated transient potentials of open-layered media by parallel computing using OPENMP, MPI, OPENACC, and CUDA. IEEE Trans Antenn Propag. 2019:doi:10.1109/TAP.2019.2920253 (epub ahead of print).
45.
KangHLeeSWLeeESKimSHLeeTG. Real-time GPU-accelerated processing and volumetric display for wide-field laser-scanning optical-resolution photoacoustic microscopy. Biomed Opt Express. 2015;6(12):4650-60.
46.
ShanTQiJJiangMJiangH. GPU-based acceleration and mesh optimization of finite-element-method-based quantitative photoacoustic tomography: a step towards clinical applications. Appl Opt. 2017;56(15):4426-32.
47.
DingLDeán-BenXLRazanskyD. Real-time model-based inversion in cross-sectional optoacoustic tomography. IEEE Trans Med Imaging. 2016;35(8):1883-91.
48.
ShanT. Accelerated time domain quantitative photoacoustic tomography (TD-qPAT) based on graphic processing units (GPU) for clinical application of breast cancer imaging. In: Optical Tomography and Spectroscopy, paper OW4D.3, Optical Society of America, 3-6April2018. Available from https://www.osapublishing.org/abstract.cfm?uri=OTS-2018-OW4D.3.
49.
LiuSFengXGaoFJinHZhangRLuoYZhengY. GPU-accelerated two dimensional synthetic aperture focusing for photoacoustic microscopy. APL Photon. 2018;3(2):026101.
50.
WangKHuangCKaoYJChouCYOraevskyAAAnastasioMA. Accelerating image reconstruction in three-dimensional optoacoustic tomography on graphics processing units. Med Phys. 2013;40(2): 023301.
51.
LutzweilerCDeán-BenXLRazanskyD. Expediting model-based optoacoustic reconstructions with tomographic symmetries. Med Phys. 2014;41(1):013302.
52.
ArridgeSBeardPBetckeMCoxBHuynhNLuckaFOgunladeOZhangE. Accelerated high-resolution photoacoustic tomography via compressed sensing. Phys Med Biol. 2016;61(24):8908.
53.
WenTLiLZhuQQinWGuJYangFXieY. GPU-accelerated kernel regression reconstruction for freehand 3D ultrasound imaging. Ultrason Imaging. 2017;39(4):240-59.
GamelinJMaurudisAAguirreAHuangFGuoPWangLVZhuQ. A real-time photoacoustic tomography system for small animals. Opt Express. 2009;17(13):10489-98.
56.
TaruttisANtziachristosV. Advances in real-time multispectral optoacoustic imaging and its applications. Nat Photon. 2015;9(4):219-27.
57.
LiMLiuCGongXZhengRBaiYXingMDuXLiuXZengJLinRZhouHWangSLuGZhuWFangCSongL. Linear array-based real-time photoacoustic imaging system with a compact coaxial excitation handheld probe for noninvasive sentinel lymph node mapping. Biomed Opt Express. 2018;9(4):1408-22.
58.
ArnalBWeiCWPerezCNguyenTMLombardoMPelivanovIPozzoLDO’DonnellM. Sono-photoacoustic imaging of gold nanoemulsions: part II. real time imaging. Photoacoustics. 2015;3(1):11-9.
CuiHYangX. Real-time monitoring of high-intensity focused ultrasound ablations with photoacoustic technique: an in vitro study. Med Phys. 2011;38(10):5345-50.
61.
SongLKimCMaslovKShungKKWangLV. High-speed dynamic 3D photoacoustic imaging of sentinel lymph node in a murine model using an ultrasound array. Med Phys. 2009;36(8):3724-9.
HaririALemasterJWangJJeevarathinamASChaoDLJokerstJV. The characterization of an economic and portable LED-based photoacoustic imaging system to facilitate molecular imaging. Photoacoustics. 2018;9:10-20.
64.
RostamiSRMGhaffari-MiabM. Fast computation of finite difference generated time-domain green’s functions of layered media using OpenAcc on graphics processors. In: 2017 Iranian Conference on Electrical Engineering, pp. 1596-99. doi:10.1109/IranianCEE.2017.7985300.
65.
RuetschGFaticaM. CUDA FORTRAN for Scientists and Engineers. Santa Clara, CA: NVIDIA Corporation; 2011.
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
