Rank Minimization-Based Regularization Method for Sparse-View Photoacoustic Image Reconstruction

Abstract

The photoacoustic tomography (PAT) is a new biomedical imaging modality. It has great advantages in early diagnosis of the human disease and accurate monitoring of disease progression. In photoacoustic imaging, when a beam of short-pulsed laser illuminates the biological tissue, the photoacoustic effect leads to the emergence of acoustic waves in the tissue. The initial acoustic pressure in the tissue reveals the structures of the tissue. The purpose of the PAT reconstruction problem is to obtain the initial acoustic pressure in the tissue from the collected photoacoustic signal information. In this paper, we propose a rank minimization-based regularization model for the sparse-view photoacoustic image reconstruction problem. We design a proximal alternating iterative algorithm to solve the model and the convergence of the algorithm is demonstrated by utilizing the Kudyka-Łojasiewicz theory. The experimental results show that the proposed method is competitive with the existing state-of-the-art PAT reconstruction methods in terms of both reconstructed quantities and visual effects for the sparse-view PAT reconstruction problem.

Keywords

Photoacoustic tomography (PAT) reconstruction / Rank minimization / Proximal alternating iterative method / Kudyka-Łojasiewicz theory / 68U10 / 94A08

References

[1] ArridgeS, BeardP, BetckeM, CoxB, HuynhN, LuckaF, OgunladeO, ZhangE. Accelerated high-resolution photoacoustic tomography via compressed sensing. Phys. Med. Biol., 2016, 61248908

[2] AttouchH, BolteJ, RedontP, SoubeyranA. Proximal alternating minimization and projection methods for nonconvex problems: an approach based on the Kurdyka-Łojasiewicz inequality. Math. Oper. Res., 2010, 35(2): 438-457

[3] Attouch, H., Bolte, J., Svaiter, B.F.: Convergence of descent methods for semi-algebraic and tame problems: proximal algorithms, forward-backward splitting, and regularized Gauss-Seidel methods. Math. Program. 137(1/2), 91–129 (2013)

[4] Bai, Z.-Z.: Sharp error bounds of some Krylov subspace methods for non-Hermitian linear systems. Appl. Math. Comput. 109(2/3) 273–285 (2000)

[5] BaiZ-Z. Motivations and realizations of Krylov subspace methods for large sparse linear systems. J. Comput. Appl. Math., 2015, 283: 71-78

[6] BaiZ-Z, PanJ-YMatrix Analysis and Computations, 2021, Philadelphia. Society for Industrial and Applied Mathematics.

[7] BaiZ-Z, WangL. On convergence rates of Kaczmarz-type methods with different selection rules of working rows. Appl. Numer. Math., 2023, 186: 289-319

[8] Bai, Z.-Z.: Wang, L., Muratova, G.V.: On relaxed greedy randomized augmented Kaczmarz methods for solving large sparse inconsistent linear systems. East Asian J. Appl. Math. 12(2), 323–332 (2022)

[9] BaiZ-Z, WuW-T. On greedy randomized Kaczmarz method for solving large sparse linear systems. SIAM J. Sci. Comput., 2018, 40(1): 592-606

[10] BaiZ-Z, WuW-T. On relaxed greedy randomized Kaczmarz methods for solving large sparse linear systems. Appl. Math. Lett., 2018, 83: 21-26

[11] BaiZ-Z, WuW-T. On greedy randomized coordinate descent methods for solving large linear least-squares problems. Numer. Linear Algebra, 2019, 2642237

[12] BaiZ-Z, WuW-T. On greedy randomized augmented Kaczmarz method for solving large sparse inconsistent linear systems. SIAM J. Sci. Comput., 2021, 43(6): 3892-3911

[13] BaiZ-Z, WuW-T. Randomized Kaczmarz iteration methods: algorithmic extensions and convergence theory. Jpn. J. Ind. Appl. Math., 2023, 40(3): 1421-1443

[14] BarryJR, LeeEA, MesserschmittDGDigital Communication, 2012, Berlin. Springer.

[15] BellAG. On the production and reproduction of sound by light. Am. J. Sci., 1880, 3(118): 305-324

[16] BolteJ, DaniilidisA, LewisA. The Łojasiewicz inequality for nonsmooth subanalytic functions with applications to subgradient dynamical systems. SIAM J. Optimiz., 2007, 17(4): 1205-1223

[17] BolteJ, SabachS, TeboulleM. Proximal alternating linearized minimization for nonconvex and nonsmooth problems. Math. Program., 2014, 146(1/2): 459-494

[18] Dean-BenXL, NtziachristosV, RazanskyD. Acceleration of optoacoustic model-based reconstruction using angular image discretization. IEEE Trans. Med. Imaging, 2012, 31(5): 1154-1162

[19] DingL, Dean-BenXL, RazanskyD. Real-time model-based inversion in cross-sectional optoacoustic tomography. IEEE Trans. Med. Imaging, 2016, 35(8): 1883-1891

[20] GengT, SunG, XuY, HeJF. Truncated nuclear norm minimization based group sparse representation for image restoration. SIAM J. Imaging Sci., 2018, 11(3): 1878-1897

[21] GolubGH, Van LoanCFMatrix Computations, 2013, Baltimore. Johns Hopkins University Press.

[22] HanY, DingL, BenXLD, RazanskyD, PrakashJ, NtziachristosV. Three-dimensional optoacoustic reconstruction using fast sparse representation. Opt. Lett., 2017, 42(5): 979-982

[23] HuangYM, YanHY. Weighted nuclear norm minimization-based regularization method for image restoration. Commun. Appl. Math. Comput., 2021, 3: 371-389

[24] HuangYM, YanHY, WenYW, YangX. Rank minimization with applications to image noise removal. Inform. Sci., 2018, 429: 147-163

[25] Jansen, K., van-Soest, G., van-der-Steen, A.F.W.: Intravascular photoacoustic imaging: a new tool for vulnerable plaque identification. Ultrasound Med. Biol. 40(6), 1037–1048 (2014)

[26] JohnMJ, BarhumiI. Fast and efficient PAT image reconstruction algorithms: a comparative performance analysis. Signal Process., 2022, 201108691

[27] KrugerRA, LiuP, FangYR. Photoacoustic ultrasound (PAUS)—reconstruction tomography. Med. Phys., 1995, 2(10): 1605-1609

[28] KunyanskyLA. Explicit inversion formulae for the spherical mean Radon transform. Inverse Probl., 2007, 23(1): 373-383

[29] LauferJG, ZhangEZ, TreebyBE, CoxBT, BeardPC, JohnsonP, PedleyB. In vivo preclinical photoacoustic imaging of tumor vasculature development and therapy. J. Biomed. Opt., 2012, 175056016

[30] LiX, QiL, ZhangS, HuangSX, ChenWF. Model-based optoacoustic tomography image reconstruction with non-local and sparsity regularizations. IEEE Access, 2019, 7: 102136-102148

[31] LiuF, GongX, WangLV, GuanJ, MengJ. Dictionary learning sparse-sampling reconstruction method for in-vivo 3D photoacoustic computed tomography. Biomed. Opt. Express, 2019, 10(4): 1660-1677

[32] LiuXX, LuJ, ShenLX, XuC, XuYS. Multiplicative noise removal: nonlocal low-rank model and its proximal alternating reweighted minimization algorithm. SIAM J. Imaging Sci., 2020, 13(3): 1595-1629

[33] MordukhovichBSVariational Analysis and Generalized Differentiation, 2006, Berlin. I. Basic Theory. Springer.

[34] ShawCB, PrakashJ, PramanikM. Least squares QR-based decomposition provides an efficient way of computing optimal regularization parameter in photoacoustic tomography. J. Biomed. Opt., 2013, 18(8): 080501-080501

[35] SivasubramanianK, PeriyasamyV, PramanikM. Non-invasive sentinel lymph node mapping and needle guidance using clinical handheld photoacoustic imaging system in small animal. J. Biophotonics, 2018, 111201700061

[36] TangJ, ColemanJE, DaiX, JiangH. Wearable 3-D photoacoustic tomography for functional brain imaging in behaving rats. Sci. Rep., 2016, 6125470

[37] WangZ, BovikAC, SheikhHR, SimoncelliEP. Image quality assessment: from error visibility to structural similarity. IEEE Trans. Image Process., 2004, 13: 600-612

[38] WangZX, WangHQ, RenSL. Research on ADMM reconstruction algorithm of photoacoustic tomography with limited sampling data. IEEE Access, 2021, 9: 113631-113641

[39] XiaJ, ChatniMR, MaslovK, GuoZJ, WangK, AnastasioM, WangLV. Whole-body ring-shaped confocal photoacoustic computed tomography of small animals in vivo. J. Biomed. Opt., 2012, 175050506

[40] XieY, QuY, TaoD, WuW, YuanQ, ZhangW. Hyperspectral image restoration via iteratively regularized weighted Schatten $ p $-norm minimization. IEEE Trans. Geosci Remote Sens., 2016, 54(8): 4642-4659

[41] XuMH, WangLV. Time-domain reconstruction for thermoacoustic tomography in a spherical geometry. IEEE Trans. Med. Imaging, 2002, 21(7): 814-822

[42] XuMH, WangLV. Universal back-projection algorithm for photoacoustic computed tomography. Phys. Rev. E, 2005, 711016706

[43] XuMH, XuY, WangLV. Time-domain reconstruction algorithms and numerical simulations for thermoacoustic tomography in various geometries. IEEE Trans. Biomed. Eng., 2003, 50(9): 1086-1099

[44] XuY, FengD, WangLV. Exact frequency-domain reconstruction for thermoacoustic tomography. I. Planar geometry. IEEE Trans. Med. Imaging, 2002, 21(7): 823-828

[45] Xu, Y., Xu, M.H., Wang, L.V.: Exact frequency-domain reconstruction for thermoacoustic tomography. II. Cylindrical geometry. IEEE Trans. Med. Imaging 21(7), 829–833 (2002)

[46] YamagaI, Kawaguchi-SakitaN, AsaoY. Vascular branching point counts using photoacoustic imaging in the superficial layer of the breast: a potential biomarker for breast cancer. Photoacoustics, 2018, 11: 6-13

[47] YanHY, HuangYM, YuYC. A matrix rank minimization-based regularization method for image restoration. Digit. Signal Process., 2022, 130103694

[48] YanHY, ZhengZ. Image cartoon-texture decomposition by a generalized non-convex low-rank minimization method. J Franklin Inst., 2024, 361(2): 796-815

[49] YuY, PengJ, YueS. A new nonconvex approach to low-rank matrix completion with application to image inpainting. Multidimens. Syst. Signal Process., 2019, 30: 145-174

[50] ZhangY, WangYY, ZhangC. Total variation based gradient descent algorithm for sparse-view photoacoustic image reconstruction. Ultrasonics, 2012, 52(8): 1046-1055

RIGHTS & PERMISSIONS

Shanghai University