Integration of biomechanical models into image registration in the presence of large deformations

Prone-to-supine breast image registration has potential application in the fields of surgical and radiotherapy planning, and image guided interventions. However, breast image registration of three-dimensional images acquired in different patient positions is a challenging problem, due to large deformations induced to the soft breast tissue caused by the change in gravity loading. Biomechanical modelling is a promising tool to predict gravity induced deformations, however such simulations alone are unlikely to produce good alignment due to inter-patient variability and image acquisition related influences on the breast shape. This thesis presents a symmetric, biomechanical simulation based registration framework which aligns images in a central, stress-free configuration. Soft tissue is modelled as a neo-Hookean material and external forces are considered as the main source of deformation in the original images. The framework successively applies image derived forces directly into the unloading simulation in place of a subsequent image registration step. This results in a biomechanically constrained deformation. Using a finite difference scheme enables simulations to be performed directly in the image space. Motion constrained boundary conditions have been incorporated which can capture tangential motion of membranes and fasciae. The accuracy of the approach is assessed by measuring the target registration error (TRE) using nine prone MRI and supine CT image pairs, one prone-supine CT image pair, and four prone-supine MRI image pairs. The registration reduced the combined mean TRE for all clinical data sets from initially 69.7mm to 5.6mm. Prone-supine image pairs might not always be available in the clinical breast cancer workflow, especially prior to surgery. Hence an alternative surface driven registration methodology was also developed that incorporates biomechanical simulations, material parameter optimisation, and constrained surface matching. For three prone MR images and corresponding supine CT-derived surfaces a final mean TRE of 10.0mm was measured

3-D lung deformation and function from respiratory-gated 4-D x-ray CT images : application to radiation treatment planning.

Many lung diseases or injuries can cause biomechanical or material property changes that can alter lung function. While the mechanical changes associated with the change of the material properties originate at a regional level, they remain largely asymptomatic and are invisible to global measures of lung function until they have advanced significantly and have aggregated. In the realm of external beam radiation therapy of patients suffering from lung cancer, determination of patterns of pre- and post-treatment motion, and measures of regional and global lung elasticity and function are clinically relevant. In this dissertation, we demonstrate that 4-D CT derived ventilation images, including mechanical strain, provide an accurate and physiologically relevant assessment of regional pulmonary function which may be incorporated into the treatment planning process. Our contributions are as follows: (i) A new volumetric deformable image registration technique based on 3-D optical flow (MOFID) has been designed and implemented which permits the possibility of enforcing physical constraints on the numerical solutions for computing motion field from respiratory-gated 4-D CT thoracic images. The proposed optical flow framework is an accurate motion model for the thoracic CT registration problem. (ii) A large displacement landmark-base elastic registration method has been devised for thoracic CT volumetric image sets containing large deformations or changes, as encountered for example in registration of pre-treatment and post-treatment images or multi-modality registration. (iii) Based on deformation maps from MOFIO, a novel framework for regional quantification of mechanical strain as an index of lung functionality has been formulated for measurement of regional pulmonary function. (iv) In a cohort consisting of seven patients with non-small cell lung cancer, validation of physiologic accuracy of the 4-0 CT derived quantitative images including Jacobian metric of ventilation, Vjac, and principal strains, (V?1, V?2, V?3, has been performed through correlation of the derived measures with SPECT ventilation and perfusion scans. The statistical correlations with SPECT have shown that the maximum principal strain pulmonary function map derived from MOFIO, outperforms all previously established ventilation metrics from 40-CT. It is hypothesized that use of CT -derived ventilation images in the treatment planning process will help predict and prevent pulmonary toxicity due to radiation treatment. It is also hypothesized that measures of regional and global lung elasticity and function obtained during the course of treatment may be used to adapt radiation treatment. Having objective methods with which to assess pre-treatment global and regional lung function and biomechanical properties, the radiation treatment dose can potentially be escalated to improve tumor response and local control