Coupling brain-tumor biophysical models and diffeomorphic image registration

DOI10.1016/J.CMA.2018.12.008zbMATH Open1440.92043arXiv1710.06420OpenAlexW2766956330WikidataQ92133231 ScholiaQ92133231MaRDI QIDQ1987828FDOQ1987828

Christos Davatzikos, Amir Gholami, Klaudius Scheufele, Andreas Mang, Miriam Schulte, George Biros

Publication date: 16 April 2020

Published in: Computer Methods in Applied Mechanics and Engineering (Search for Journal in Brave)

Abstract: We present the SIBIA (Scalable Integrated Biophysics-based Image Analysis) framework for joint image registration and biophysical inversion and we apply it to analyse MR images of glioblastomas (primary brain tumors). In particular, we consider the following problem. Given the segmentation of a normal brain MRI and the segmentation of a cancer patient MRI, we wish to determine tumor growth parameters and a registration map so that if we "grow a tumor" (using our tumor model) in the normal segmented image and then register it to the patient segmented image, then the registration mismatch is as small as possible. We call this "the coupled problem" because it two-way couples the biophysical inversion and registration problems. In the image registration step we solve a large-deformation diffeomorphic registration problem parameterized by an Eulerian velocity field. In the biophysical inversion step we estimate parameters in a reaction-diffusion tumor growth model that is formulated as a partial differential equation. In this paper, our contributions are the introduction of a PDE-constrained optimization formulation of the coupled problem, the derivation of the optimality conditions, and the derivation of a Picard iterative scheme for the solution of the coupled problem. In addition, we perform several tests to experimentally assess the performance of our method on synthetic and clinical datasets. We demonstrate the convergence of the SIBIA optimization solver in different usage scenarios. We demonstrate that we can accurately solve the coupled problem in three dimensions ($256^3$ resolution) in a few minutes using 11 dual-x86 nodes. Also, we demonstrate that, with our coupled approach, we can successfully register normal MRI to tumor-bearing MRI while obtaining Dice coefficients that match those achieved when registering of normal-to-normal MRI.

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

zbMATH Keywords

tumor growthparallel algorithmsadjoint-based methodsatlas registrationbiophysically constrained diffeomorphic image registration

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) Biomedical imaging and signal processing (92C55) Image processing (compression, reconstruction, etc.) in information and communication theory (94A08) Reaction-diffusion equations (35K57) Newton-type methods (49M15) Parallel algorithms in computer science (68W10)

Cites Work

• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• FAIR
• Computing large deformation metric mappings via geodesic flows of diffeomorphisms
• Optimization with PDE Constraints
• Computational Optimization of Systems Governed by Partial Differential Equations
• Truncated-Newton algorithms for large-scale unconstrained optimization
• An inverse problem formulation for parameter estimation of a reaction-diffusion model of low grade gliomas
• Patient-specific mathematical neuro-oncology: using a simple proliferation and invasion tumor model to inform clinical practice
• Choosing the Forcing Terms in an Inexact Newton Method
• Algorithms for PDE-constrained optimization
• Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations
• An inverse problem for the recovery of the vascularization of a tumor
• A mathematical method for parameter estimation in a tumor growth model
• An image-driven parameter estimation problem for a reaction-diffusion glioma growth model with mass effects
• Large-scale PDE-constrained optimization
• Selection, calibration, and validation of models of tumor growth
• Adjoint method for a tumor invasion PDE-constrained optimization problem in 2D using adaptive finite element method
• Selection and validation of predictive models of radiation effects on tumor growth based on noninvasive imaging data
• Brain–Tumor Interaction Biophysical Models for Medical Image Registration
• Constrained $H^1$-Regularization Schemes for Diffeomorphic Image Registration
• An Inexact Newton--Krylov Algorithm for Constrained Diffeomorphic Image Registration
• A fully coupled space-time multiscale modeling framework for predicting tumor growth
• 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 In (5)

• An operator-splitting approach for variational optimal control formulations for diffeomorphic shape matching
• Where did the tumor start? An inverse solver with sparse localization for tumor growth models
• CLAIRE: A Distributed-Memory Solver for Constrained Large Deformation Diffeomorphic Image Registration
• Image-Driven Biophysical Tumor Growth Model Calibration
• Analyzing the role of the Inf-Sup condition for parameter identification in saddle point problems with application in elasticity imaging

Uses Software

• TAO
• FAIR.m

This page was built for publication: Coupling brain-tumor biophysical models and diffeomorphic image registration

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