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.92067OpenAlexW3102357792MaRDI QIDQ5000579FDOQ5000579
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.
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
Cites Work
- scientific article; zbMATH DE number 42093 (Why is no real title available?)
- scientific article; zbMATH DE number 5060482 (Why is no real title available?)
- A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
- A multilayer grow-or-go model for GBM: effects of invasive cells and anti-angiogenesis on growth
- A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type. (Reprint)
- A quasi Tikhonov regularization for a two-dimensional backward heat problem by a fundamental solution
- Adjoint method for a tumor invasion PDE-constrained optimization problem in 2D using adaptive finite element method
- An image-driven parameter estimation problem for a reaction-diffusion glioma growth model with mass effects
- An inverse problem for the recovery of the vascularization of a tumor
- An inverse problem formulation for parameter estimation of a reaction-diffusion model of low grade gliomas
- An iterative thresholding algorithm for linear inverse problems with multi-constraints and its applications
- Bayesian calibration, validation, and uncertainty quantification of diffuse interface models of tumor growth
- Brain–Tumor Interaction Biophysical Models for Medical Image Registration
- CoSaMP: Iterative signal recovery from incomplete and inaccurate samples
- Coupling brain-tumor biophysical models and diffeomorphic image registration
- Density-dependent quiescence in glioma invasion: instability in a simple reaction–diffusion model for the migration/proliferation dichotomy
- Heat source identification based on \(\ell_1\) constrained minimization
- Ill-posed medicine—an introduction to image registration
- Iterative hard thresholding for compressed sensing
- Line search algorithms with guaranteed sufficient decrease
- Methoden zur numerischen Simulation der Progression von Gliomen
- Numerical reconstruction of brain tumours
- On Tikhonov regularization with non-convex sparsity constraints
- Recovering the source and initial value simultaneously in a parabolic equation
- Recovery of sparse signals using OMP and its variants: convergence analysis based on RIP
- Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information
- Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit
- Simulation of glioblastoma growth using a 3D multispecies tumor model with mass effect
- Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit
Cited In (14)
- Coupling brain-tumor biophysical models and diffeomorphic image registration
- A Fast and Scalable Computational Framework for Large-Scale High-Dimensional Bayesian Optimal Experimental Design
- Analyzing the role of the inf-sup condition for parameter identification in saddle point problems with application in elasticity imaging
- An inverse problem formulation for parameter estimation of a reaction-diffusion model of low grade gliomas
- Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics
- Joint state-parameter estimation for tumor growth model
- 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
- A heterogeneous patient-specific model of glioblastoma multiforme tumor through an inverse problem
- Numerical reconstruction of brain tumours
- 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
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)
