EF3 – Model-based Imaging

Project

EF3-11

Quantitative Tissue Pressure Imaging via PDE-Informed Assimilation of MR Data

Project Heads

Alfonso Caiazzo, Karsten Tabelow, Ingolf Sack

Project Members

Cristian Cárcamo Sánchez

Former Members

Felipe Galarce Marín (currently Assistant Professor, PUCV, Valparaiso, Chile)

Project Duration

01.01.2022 − 31.12.2023

Located at

WIAS

Description

Background

Medical imaging of tissues based on MRI enables the non-invasive quantitative characterization of important biomarkers, and it is therefore a pillar of clinical diagnostic of tissue diseases. Biological tissues are by nature poroelastic and the understanding of biphasic aspects such as porosity and interstitial pressure might provide the identification of abnormal pressure increase in the cardiovascular system or cerebro-spinal fluid. To this purpose, the mechanical behavior of the tissue is modeled as a biphasic (poroelastic) material.

Magnetic resonance elastography (MRE) is an imaging technique sensitive to the stiffness properties of tissues. In MRE, the propagation of shear waves in the audible frequency range (10 – 1000 Hz) through the tissue – recorded as a three-dimensional image of the inner displacement field. MRE has been largely used in the context of inverse problems for elastic and viscoelastic tissues. However, to characterize pathological pressure gradients, the data assimilation problem shall take into account (i) a poroelastic underlying physical model and (ii) both the estimation of mechanical parameters and the reconstruction of the physical state (i.e., the IFP – interstitial fluid pressure). In particular, the estimation of pressure gradients is extremely demanding, due to the high dimensionality of the unknown space and to the limited availability and quality of experimental data. In fact, MRE acquisition is practically constrained by the length of examination time, and displacement data are typically available only on a subregion of the tissue of interest and with a limited resolution.

 

Within the MATH+ project EF3-9 we established a new statistical framework to infer on porosity maps from inversion recovery MR data and we developed a pipeline to setup a patient-specific poroelastic model from MPRAGE and IRMRI data. The project EF3-11 targets the quantitative estimation of interstitial pressure from medical imaging. To this purpose, we aim at developing a mathematical framework to solve the data assimilation problem in a physics-informed reduced space, i.e., a manifold of numerical solutions of the poroelastic equations modeling the tissue mechanics.

 

Methods

We tackle data assimilation of MRE data combining (i) a state reconstruction over the whole domain and (ii) a variational projection method for the estimation of pressure gradients. The first step allows to exploit the underlying physical model to extend the available MRE measurement to the full 3D space and avoid the issue of setting internal pressure boundaries.

The displacement reconstruction algorithm will be based on the so-called parametrized- background data-weak (PBDW) method [Maday et al. 2014], an approach designed to reconstruct a physical solution over the whole domain given only partial information. Originally introduced for wave equations, the PBDW has been recently extended in the context of ultrasound images of blood flow [Galarce et al., 2021]. In this approach, one generates a manifold of suitable solutions of an underlying PDE dependent on few physical parameters, by solving the problem on a finite approximation of the parameter space. Then, the reconstruction is sought on a low-dimensional subspace of this PDE- informed manifold, solving an optimization problem depending on the available observations and including additional terms to account for the model bias

 

Project Publications

F. Galarce, A. Caiazzo, J. Mura

Bias and multiscale correction methods for variational state estimation

Abstract: Data assimilation performance can be significantly impacted by biased noise in observations, altering the signal magnitude and introducing fast oscillations or discontinuities when the system lacks smoothness. To mitigate these issues, this paper employs variational state estimation using the so-called parametrized-background data-weak method. This approach relies on a background manifold parametrized by a set of constraints, enabling the state estimation by solving a minimization problem on a reduced-order background model, subject to constraints imposed by the input measurements. The proposed formulation incorporates a novel bias correction mechanism and a manifold decomposition that handles rapid oscillations by treating them as slow-decaying modes based on a two-scale splitting of the classical reconstruction algorithm. The method is validated in different examples, including the assimilation of biased synthetic data, discontinuous signals, and Doppler ultrasound data obtained from experimental measurements.

 

C. Cárcamo, A. Caiazzo, F. Galarce, J. Mura

A stabilized total pressure-formulation of the Biot’s poroelasticity equations in frequency domain: Numerical analysis and applications

Comp. Meth. Appl. Mech. Engnr., Vol 432(A), pp. 117345 (2024) – WIAS Preprint 3101 (2022)

Abstract: This work focuses on the numerical solution of the dynamics of a poroelastic material in the frequency domain. We provide a detailed stability analysis based on the application of the Fredholm alternative in the continuous case, considering a total pressure formulation of the Biot’s equations. In the discrete setting, we propose a stabilized equal order finite element method complemented by an additional pressure stabilization to enhance the robustness of the numerical scheme with respect to the fluid permeability. Utilizing the Fredholm alternative, we extend the well-posedness results to the discrete setting, obtaining theoretical optimal convergence for the case of linear finite elements. We present different numerical experiments to validate the proposed method. First, we consider model problems with known analytic solutions in two and three dimensions. As next, we show that the method is robust for a wide range of permeabilities, including the case of discontinuous coefficients. Lastly, we show the application for the simulation of brain elastography on a realistic brain geometry obtained from medical imaging.

 

F. Galarce, K. Tabelow, J. Polzehl, C. Panagiotis Papanikas, V. Vavourakis, L. Lilaj, I. Sack, A. Caiazzo.

Displacement and pressure reconstruction from magnetic resonance elastography images: Application to an in silico brain model

SIAM J Imaging Sciences, Vol 16(2), pp. 996-1027 (2023)WIAS Preprint 2933 (2022)

Abstract: Magnetic resonance elastography is a motion sensitive image modality that allows to measure in vivo tissue displacement fields in response to mechanical excitations. This paper investigates a data assimilation approach for reconstructing tissue displacement and pressure fields in an in silico brain model from partial elastography data. The data assimilation is based on a parametrized-background data weak methodology, in which the state of the physical system — tissue displacements and pressure fields — is reconstructed from the available data assuming an underlying poroelastic biomechanics model. For this purpose, a physics-informed manifold is built by sampling the space of parameters describing the tissue model close to their physiological ranges, to simulate the corresponding poroelastic problem, and compute a reduced basis. Displacements and pressure reconstruction is sought in a reduced space after solving a minimization problem that encompasses both the structure of the reduced-order model and the available measurements. The proposed pipeline is validated using synthetic data obtained after simulating the poroelastic mechanics on a physiological brain. The numerical experiments demonstrate that the framework can exhibit accurate joint reconstructions of both displacement and pressure fields. The methodology can be formulated for an arbitrary resolution of available displacement data from pertinent images. It can also inherently handle uncertainty on the physical parameters of the mechanical model by enlarging the physics-informed manifold accordingly. Moreover, the framework can be used to characterize, in silico, biomarkers for pathological conditions, by appropriately training the reduced-order model.

News & Events

Workshop on Biophysics-based nodeling and data assimilation in medical imaging (BioPhysMed 2023)

When & Where: WIAS Berlin, August 30 – September 1, 2023

More info here

Computational & Mathematical Biomedical Engineering (CMBE 2022)

When & Where: June 27-29, in Milan (Politecnico)

Minisymposium: Inverse problems in biofluids and -solids

More info here

 

SIAM Imaging Science (IS) Conference

When: March 21-25, Online

Minisymposium: Biomechanical Model-Based MR Imaging, Inverseproblems, and Applications

 

Presentations:

MRE & Mechanics Session:

  • Jing Guo (MRE Group, Charité, Berlin)
  • Felipe Galarce (WIAS Berlin)
  • Federica Caforio (Medical U of Graz)
  • Sebastien Imperiale (INRIA Saclay)

Cardiovascular MRI Session:

  • David Nordsletten (U Michigan)
  • Katerina Skardova (U Prag)
  • Jorge Aguayo (U Chile)
  • Jeremía Garay (U Groningen)

Relevant Publications

N. Jaitner, Y. Safraou, M. Anders, J. Schattenfroh, T. Meyer, B. Huang, J. Jordan, O. Bohem, A. Caiazzo, T. Schaefftler, J. Mura, J. Guo, I. Sack. Noninvasive assessment of portal pressure by combined measurement of volumetric strain and stiffness of in vivo human liver. Acta Biomateralia (2025) (in press, availanle online)

F. Galarce, A. Caiazzo, J. Mura. Bias and multiscale correction methods for variational state estimation. Appl. Math. Mod., Vol 138(A), pp. 115761 (2025)

C. Cárcamo, A. Caiazzo. F. Galarce, J. Mura. A stabilized total pressure-formulation of the Biot’s poroelasticity equations in frequency domain: Numerical analysis and applications Comp. Meth. Appl. Mech. Engnr., Vol 432(A), pp. 117345 (2024) 

F. Galarce, K. Tabelow, J. Polzehl, C.P. Papanikas, V. Vavourakis, L. Lilaj, I. Sack, A. Caiazzo. Displacement and pressure reconstruction from magnetic resonance elastography images: Application to an in silico brain model. SIAM J Imaging Sciences, Vol 16(2), pp. 996-1027 (2023)WIAS Preprint 2933 (2022)

 

F. Galarce, J.-F. Gerbeau, D. Lombardi, and O. Mula. Fast reconstruction of 3d blood flows from doppler ultrasound images and reduced models. Comp. Meth. Appl. Mech. Engnr., 375:113559, 2021.

L. Lilaj, H. Herthum, T. Meyer, M. Shahryari, G. Bertalan, A. Caiazzo, J. Braun, T. Fischer, S. Hirsch, and I. Sack. Inversion recovery MR elastography of the human brain for improved stiffness quantification near fluid-solid boundaries. Mag. Res. Med.  2021 86(5):2552-2561

Y. Maday, A.T. Patera, J.D. Penn, and M. Yano. A Parameterized-Background Data-Weak approach to variational data assimilation: formulation, analysis, and application to acoustics. Int. J. Num. Methods Engnr., 102(5):933–965, 2014

I. Sack and T. Schäffter, editors. Quantification of Biophysical Parameters in Medical Imaging. Springer, 2018

Related Pictures

View of MPRAGE Brain image. The red area indicates the subregion where MRE data are available

Patient-specific computational mesh created from the anatomical images

 

Example of forward simulation: displacement field

Example of forward simulation: pressure field