Compliance boundary conditions for patient-specific deformation simulation using the finite element method

Estimation of boundary conditions for patient-specific liver simulation during augmented surgery

PURPOSE

Augmented reality can improve the outcome of hepatic surgeries, assuming an accurate liver model is available to estimate the position of internal structures. While researchers have proposed patient-specific liver simulations, very few have addressed the question of boundary conditions. Resulting mainly from ligaments attached to the liver, they are not visible in preoperative images, yet play a key role in the computation of the deformation.

METHOD

We propose to estimate both the location and stiffness of ligaments by using a combination of a statistical atlas, numerical simulation, and Bayesian inference. Ligaments are modeled as polynomial springs connected to a liver finite element model. They are initialized using an anatomical atlas and stiffness properties taken from the literature. These characteristics are then corrected using a reduced-order unscented Kalman filter based on observations taken from the laparoscopic image stream.

RESULTS

Our approach is evaluated using synthetic data and phantom data. By relying on a simplified representation of the ligaments to speed up computation times, it is not estimating the true characteristics of ligaments. However, results show that our estimation of the boundary conditions still improves the accuracy of the simulation by 75% when compared to typical methods involving Dirichlet boundary conditions.

CONCLUSION

By estimating patient-specific boundary conditions, using tracked liver motion from RGB-D data, our approach significantly improves the accuracy of the liver model. The method inherently handles noisy observations, a substantial feature in the context of augmented reality.

Robust Reconstruction of Elasticity Using Ultrasound Imaging and Multi-Frequency Excitations

Biomedical parameters of tissue can be important indicators for clinical diagnosis. One such parameter that reflects tissue stiffness is elasticity, the imaging of which is called elastography. In this paper, we use displacements from harmonic excitations to solve the inverse problem of elasticity based on a finite-element method (FEM) formulation. This leads to iterative solution of nonlinear and nonconvex problems. In this paper, we show the importance and selection of viable initializations in numerical simulation studies and propose techniques for the fusion of multiple initializations for ideal reconstructions of unknown tissue as well as combining information from excitations at multiple frequencies. Results show that our method leads up to 76% decrease in root-mean-squared error (RMSE) and 9.9 dB increase in contrast-to-noise ratio (CNR) in simulations with noise, when compared to conventional iterative FEM without multiple initializations and frequencies. As the wave patterns in individually selected frequencies may introduce artifacts, a joint inverse-problem solution of multi-frequency excitations is introduced as a robust solution, where CNR improvements of up to 11.9 dB are observed. We also present the methods on a tissue-mimicking gelatin phantom study using mechanical excitation and ultrafast plane-wave ultrasound imaging, where the RMSE was improved by up to 51%. An experiment of ablation via heating an ex-vivo bovine liver shows that reconstruction artifacts are reduced with our proposed method.