PDE-constrained optimization in medical image analysis

Abstract

PDE-constrained optimization problems find many applications in medical image analysis, for example, neuroimaging, cardiovascular imaging, and oncological imaging. We review related literature and give examples on the formulation, discretization, and numerical solution of PDE-constrained optimization problems for medical imaging. We discuss three examples. The first one is image registration. The second one is data assimilation for brain tumor patients, and the third one data assimilation in cardiovascular imaging. The image registration problem is a classical task in medical image analysis and seeks to find pointwise correspondences between two or more images. The data assimilation problems use a PDE-constrained formulation to link a biophysical model to patient-specific data obtained from medical images. The associated optimality systems turn out to be sets of nonlinear, multicomponent PDEs that are challenging to solve in an efficient way. The ultimate goal of our work is the design of inversion methods that integrate complementary data, and rigorously follow mathematical and physical principles, in an attempt to support clinical decision making. This requires reliable, high-fidelity algorithms with a short time-to-solution. This task is complicated by model and data uncertainties, and by the fact that PDE-constrained optimization problems are ill-posed in nature, and in general yield high-dimensional, severely ill-conditioned systems after discretization. These features make regularization, effective preconditioners, and iterative solvers that, in many cases, have to be implemented on distributed-memory architectures to be practical, a prerequisite. We showcase state-of-the-art techniques in scientific computing to tackle these challenges.

References in corpus (7)

• Implementation and evaluation of various demons deformable image registration algorithms on GPU
• A coupled surface-Cahn--Hilliard bulk-diffusion system modeling lipid raft formation in cell membranes
• An inverse problem formulation for parameter estimation of a reaction diffusion model of low grade gliomas
• Constrained $H^1$-regularization schemes for diffeomorphic image registration
• An inexact Newton-Krylov algorithm for constrained diffeomorphic image registration
• A Lagrangian Gauss-Newton-Krylov Solver for Mass- and Intensity-Preserving Diffeomorphic Image Registration
• A Semi-Lagrangian two-level preconditioned Newton-Krylov solver for constrained diffeomorphic image registration

Cited by in corpus (7)

• Integrated Biophysical Modeling and Image Analysis: Application to Neuro-Oncology
• Optimizing a DIscrete Loss (ODIL) to solve forward and inverse problems for partial differential equations using machine learning tools
• CLAIRE: A distributed-memory solver for constrained large deformation diffeomorphic image registration
• Image-Driven Biophysical Tumor Growth Model Calibration
• CUQIpy: II. Computational uncertainty quantification for PDE-based inverse problems in Python
• Multi-Node Multi-GPU Diffeomorphic Image Registration for Large-Scale Imaging Problems
• Consistency of sample-based stationary points for infinite-dimensional stochastic optimization