Simulation-based joint estimation of body deformation and elasticity parameters for medical image analysis.

Estimation of tissue stiffness is an important means of noninvasive cancer detection. Existing elasticity reconstruction methods usually depend on a dense displacement field (inferred from ultrasound orMR images) and known external forces.Many imaging modalities, however, cannot provide details within an organ and therefore cannot provide such a displacement field. Furthermore, force exertion and measurement can be difficult for some internal organs, making boundary forces another missing parameter. We propose a general method for estimating elasticity and boundary forces automatically using an iterative optimization framework, given the desired (target) output surface. During the optimization, the input model is deformed by the simulator, and an objective function based on the distance between the deformed surface and the target surface is minimized numerically. The optimization framework does not depend on a particular simulation method and is therefore suitable for different physical models. We show a positive correlation between clinical prostate cancer stage (a clinical measure of severity) and the recovered elasticity of the organ. Since the surface correspondence is established, our method also provides a non-rigid image registration, where the quality of the deformation fields is guaranteed, as they are computed using a physics-based simulation.

Integration of a Multigrid ODE solver into an open medical simulation framework.

In this paper, we present the implementation of a Multigrid ODE solver in SOFA framework. By combining the stability advantage of coarse meshes and the transient detail preserving virtue of fine meshes, Multigrid ODE solver computes more efficiently than classic ODE solvers based on a single level discretization. With the ever wider adoption of the SOFA framework in many surgical simulation projects, introducing this Multigrid ODE solver into SOFA's pool of ODE solvers shall benefit the entire community. This contribution potentially has broad ramifications in the surgical simulation research community, given that in a single-resolution system, a constitutively realistic interactive tissue response, which presupposes large elements, is in direct conflict with the need to represent clinically relevant critical tissues in the simulation, which are typically be comprised of small elements.

Neurosurgery Simulation Using Non-linear Finite Element Modeling and Haptic Interaction.

Real-time surgical simulation is becoming an important component of surgical training. To meet the real-time requirement, however, the accuracy of the biomechancial modeling of soft tissue is often compromised due to computing resource constraints. Furthermore, haptic integration presents an additional challenge with its requirement for a high update rate. As a result, most real-time surgical simulation systems employ a linear elasticity model, simplified numerical methods such as the boundary element method or spring-particle systems, and coarse volumetric meshes. However, these systems are not clinically realistic. We present here an ongoing work aimed at developing an efficient and physically realistic neurosurgery simulator using a non-linear finite element method (FEM) with haptic interaction. Real-time finite element analysis is achieved by utilizing the total Lagrangian explicit dynamic (TLED) formulation and GPU acceleration of per-node and per-element operations. We employ a virtual coupling method for separating deformable body simulation and collision detection from haptic rendering, which needs to be updated at a much higher rate than the visual simulation. The system provides accurate biomechancial modeling of soft tissue while retaining a real-time performance with haptic interaction. However, our experiments showed that the stability of the simulator depends heavily on the material property of the tissue and the speed of colliding objects. Hence, additional efforts including dynamic relaxation are required to improve the stability of the system.

Image estimation from marker locations for dose calculation in prostate radiation therapy.

Tracking implanted markers in the prostate during each radiation treatment delivery provides an accurate approximation of prostate location, which enables the use of higher daily doses with tighter margins of the treatment beams and thus improves the efficiency of the radiotherapy. However, the lack of 3D image data with such a technique prevents calculation of delivered dose as required for adaptive planning. We propose to use a reference statistical shape model generated from the planning image and a deformed version of the reference model fitted to the implanted marker locations during treatment to estimate a regionally dense deformation from the planning space to the treatment space. Our method provides a means of estimating the treatment image by mapping planning image data to treatment space via the deformation field and therefore enables the calculation of dose distributions with marker tracking techniques during each treatment delivery.

Physically-based validation of deformable medical image registration.

We propose a new approach for validating deformable image registration algorithms. Since difference images do not necessarily reflect the 3D correspondence of organs, we propose to use the deformation fields generated in our FEM-based simulations to assess the displacement resulted from other registration methods. Unlike traditional FEM-based registration methods, the boundary condition for the target organ is not given explicitly. Instead it is driven by inter-organ contact forces generated by boundary conditions on surrounding organs to reduce the uncertainty induced by geometry-based surface matching. To validate our system, real CT images of the male pelvis are analyzed, and the prostate can be reasonably registered without matching its surface to the image. Several registration methods are then evaluated using our system.