Stochastic EM methods with variance reduction for penalised PET reconstructions

DOI10.1088/1361-6420/AC2D74zbMATH Open1527.35495arXiv2106.03172MaRDI QIDQ6087345FDOQ6087345

Authors: Z. Kereta, Simon Arridge, K. Thielemans, Bangti Jin 

Publication date: 15 November 2023

Published in: Inverse Problems (Search for Journal in Brave)

Abstract: Expectation-maximization (EM) is a popular and well-established method for image reconstruction in positron emission tomography (PET) but it often suffers from slow convergence. Ordered subset EM (OSEM) is an effective reconstruction algorithm that provides significant acceleration during initial iterations, but it has been observed to enter a limit cycle. In this work, we investigate two classes of algorithms for accelerating OSEM based on variance reduction for penalised PET reconstructions. The first is a stochastic variance reduced EM algorithm, termed as SVREM, an extension of the classical EM to the stochastic context, by combining classical OSEM with insights from variance reduction techniques for gradient descent. The second views OSEM as a preconditioned stochastic gradient ascent, and applies variance reduction techniques, i.e., SAGA and SVRG, to estimate the update direction. We present several numerical experiments to illustrate the efficiency and accuracy of the approaches. The numerical results show that these approaches significantly outperform existing OSEM type methods for penalised PET reconstructions, and hold great potential.

Full work available at URL: https://arxiv.org/abs/2106.03172

Recommendations

• EM-TV Methods for Inverse Problems with Poisson Noise
• The ordered subsets mirror descent optimization method with applications to tomography
• Tomographic reconstruction from Poisson distributed data: a fast and convergent EM-TV dual approach
• ADMM-EM method for \(L_1\)-norm regularized weighted least squares PET reconstruction
• Convergence results for scaled gradient algorithms in positron emission tomography

zbMATH Keywords

expectation maximizationstochastic gradientvariance reductionpositron emission tomographyordered subset expectation maximisation

Mathematics Subject Classification ID

Inverse problems for PDEs (35R30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32)

Cited In (5)

• An analysis of stochastic variance reduced gradient for linear inverse problems ^*
• Title not available (Why is that?)
• On the Convergence of Stochastic Gradient Descent for Linear Inverse Problems in Banach Spaces
• Tomographic reconstruction from Poisson distributed data: a fast and convergent EM-TV dual approach
• An approach to positron emission tomography based on penalized cross-entropy minimization.

This page was built for publication: Stochastic EM methods with variance reduction for penalised PET reconstructions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6087345)