Where did the tumor start? An inverse solver with sparse localization for tumor growth models
From MaRDI portal
Publication:5000579
DOI10.1088/1361-6420/AB649CzbMATH Open1467.92067arXiv1907.06564OpenAlexW3102357792MaRDI QIDQ5000579FDOQ5000579
Authors: Shashank Subramanian, Klaudius Scheufele, Miriam Schulte, George Biros 
Publication date: 14 July 2021
Published in: Inverse Problems (Search for Journal in Brave)
Abstract: We present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction-diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of Magnetic Resonance Imaging (MRI) scans from a single time snapshot drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. The precise time relative to tumor initiation is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint) and that the initial condition magnitude in the maximum norm equals one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over an existing solver that uses a two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the existing solver.
Full work available at URL: https://arxiv.org/abs/1907.06564
Recommendations
- Image-Driven Biophysical Tumor Growth Model Calibration
- Numerical reconstruction of brain tumours
- An inverse problem formulation for parameter estimation of a reaction-diffusion model of low grade gliomas
- Coupling brain-tumor biophysical models and diffeomorphic image registration
- An image-driven parameter estimation problem for a reaction-diffusion glioma growth model with mass effects
Reaction-diffusion equations (35K57) Inverse problems for PDEs (35R30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Pathology, pathophysiology (92C32)
Cites Work
- A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
- Title not available (Why is that?)
- Line search algorithms with guaranteed sufficient decrease
- Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information
- Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit
- Iterative hard thresholding for compressed sensing
- CoSaMP: Iterative signal recovery from incomplete and inaccurate samples
- Title not available (Why is that?)
- An inverse problem formulation for parameter estimation of a reaction-diffusion model of low grade gliomas
- A quasi Tikhonov regularization for a two-dimensional backward heat problem by a fundamental solution
- Bayesian calibration, validation, and uncertainty quantification of diffuse interface models of tumor growth
- On Tikhonov regularization with non-convex sparsity constraints
- Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit
- Heat source identification based on \(\ell_1\) constrained minimization
- Density-dependent quiescence in glioma invasion: instability in a simple reaction–diffusion model for the migration/proliferation dichotomy
- A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type. (Reprint)
- Recovery of sparse signals using OMP and its variants: convergence analysis based on RIP
- A multilayer grow-or-go model for GBM: effects of invasive cells and anti-angiogenesis on growth
- An inverse problem for the recovery of the vascularization of a tumor
- An image-driven parameter estimation problem for a reaction-diffusion glioma growth model with mass effects
- Ill-posed medicine—an introduction to image registration
- An iterative thresholding algorithm for linear inverse problems with multi-constraints and its applications
- Adjoint method for a tumor invasion PDE-constrained optimization problem in 2D using adaptive finite element method
- Numerical reconstruction of brain tumours
- Brain–Tumor Interaction Biophysical Models for Medical Image Registration
- Recovering the source and initial value simultaneously in a parabolic equation
- Simulation of glioblastoma growth using a 3D multispecies tumor model with mass effect
- Methoden zur numerischen Simulation der Progression von Gliomen
- Coupling brain-tumor biophysical models and diffeomorphic image registration
Cited In (10)
- A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design
- Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics
- Learning physics-based models from data: perspectives from inverse problems and model reduction
- An Offline-Online Decomposition Method for Efficient Linear Bayesian Goal-Oriented Optimal Experimental Design: Application to Optimal Sensor Placement
- Mathematical analysis of a model-constrained inverse problem for the reconstruction of early states of prostate cancer growth
- An inverse problem for the recovery of the vascularization of a tumor
- Analyzing the role of the Inf-Sup condition for parameter identification in saddle point problems with application in elasticity imaging
- A heterogeneous patient-specific model of glioblastoma multiforme tumor through an inverse problem
- Large-scale Bayesian optimal experimental design with derivative-informed projected neural network
- A unified Bayesian inversion approach for a class of tumor growth models with different pressure laws
Uses Software
This page was built for publication: Where did the tumor start? An inverse solver with sparse localization for tumor growth models
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5000579)
