Sampled Limited Memory Methods for Massive Linear Inverse Problems

12/17/2019
by   Julianne Chung, et al.
0

In many modern imaging applications the desire to reconstruct high resolution images, coupled with the abundance of data from acquisition using ultra-fast detectors, have led to new challenges in image reconstruction. A main challenge is that the resulting linear inverse problems are massive. The size of the forward model matrix exceeds the storage capabilities of computer memory, or the observational dataset is enormous and not available all at once. Row-action methods that iterate over samples of rows can be used to approximate the solution while avoiding memory and data availability constraints. However, their overall convergence can be slow. In this paper, we introduce a sampled limited memory row-action method for linear least squares problems, where an approximation of the global curvature of the underlying least squares problem is used to speed up the initial convergence and to improve the accuracy of iterates. We show that this limited memory method is a generalization of the damped block Kaczmarz method, and we prove linear convergence of the expectation of the iterates and of the error norm up to a convergence horizon. Numerical experiments demonstrate the benefits of these sampled limited memory row-action methods for massive 2D and 3D inverse problems in tomography applications.

READ FULL TEXT

page 14

page 16

research
01/06/2023

Immersed Boundary Method for the Complete Electrode Model in Electrical Impedance Tomography

We propose an immersed boundary scheme for the numerical resolution of t...
research
07/21/2022

Convergence analysis of multi-step one-shot methods for linear inverse problems

In this work we are interested in general linear inverse problems where ...
research
02/07/2019

Speeding up scaled gradient projection methods using deep neural networks for inverse problems in image processing

Conventional optimization based methods have utilized forward models wit...
research
03/30/2020

A Wavelet Based Sparse Row-Action Method for Image Reconstruction in Magnetic Particle Imaging

Magnetic Particle Imaging (MPI) is a preclinical imaging technique capab...
research
08/17/2022

Translation invariant diagonal frame decomposition of inverse problems and their regularization

Solving inverse problems is central to a variety of important applicatio...
research
12/13/2022

On the Polyak momentum variants of the greedy deterministic single and multiple row-action methods

For solving a consistent system of linear equations, the classical row-a...
research
05/02/2021

Bridging and Improving Theoretical and Computational Electric Impedance Tomography via Data Completion

In computational PDE-based inverse problems, a finite amount of data is ...

Please sign up or login with your details

Forgot password? Click here to reset