Recovery Assessment via Geometric Metamorphosis

Many traumatic brain injury (TBI) cases involve blunt force damage resulting from the brain hitting the skull due to rapid acceleration/deceleration, e.g., during car accidents. The damage exists as lesions with hemorrhagic and nonhemmorhagic components. In order to plan treatment for TBI, patients typically receive an MRI scan on initial presentation in the clinic, and another one four-to-eight months later. The clinical goal is to distinguish areas in these scans that represent permanent damage from those that represent transient effects (see Figure 1). Determining the regions where the hemorrhage has receded can be particularly useful in predicting long-term outcome. However, the geometry of the lesion, the deformation of the brain due to swelling, the presence of the lesion, and the infiltration of the hemorrhage into the brain can cause drastic changes between these scans and confound the determination of correspondence and hemorrhage recession.

Figure 1. MRI scans from a TBI patient. Left: Initial scan. Right: Rigidly registered scan acquired eight months later. The TBI lesion is circled in both. Notice the presence of brain deformations, the change in the lesion’s shape, and the conversion of some tissue from a diseased state (hemorrhagic) to a healthy state. Traditional registration methods fail to align these data: the mapping from one to the other cannot be fully explained by a spatial transform.

Using data from TBI patients from the University of California at Los Angeles, Kitware and the University of North Carolina at Chapel Hill have been developing registration algorithms that can help distinguish permanent damage from transient effects in longitudinal MRI scans.

Our methods can distinguish global tissue deformations, local changes in the geometry of a lesion, and local changes in the composition of the tissue and the lesion. We refer to this collection of changes as “geometric metamorphosis.” These changes are also evident in tumor cases where post-treatment assessment requires determination of changes in tumor geometry, tumor infiltration, scarring, and overall brain morphology. Similarly, in stroke cases, there is a clinical need to predict chronic changes in blood perfusion from acute scans.

Working as a team, we have produced a paper that details our algorithms and provides extensive experimental results. That paper has been accepted for publication at the Medical Image Computing and Computer-Assisted Intervention (MICCAI) conference to be held in Toronto, Canada, September 18-22, 2011. [Neithammer2011]

Related Works
Most image transformation models will fail in the presence of lesions that infiltrate or recede from healthy tissue and that deform the healthy tissue, as is common with tumors and other lesions. Direct application of a classical, deformable registration method will likely produce unrealistic deformation estimates. To overcome this limitation, deformable registration methods with weak and strong models of appearance change have also been proposed [3, 6]; proposed changes include models of tumor growth, for example. However, these methods are application-specific and degrade when their tumor growth model does not match the changes in the specific patient being registered.

Others [4, 7, 1] attempt to pre-mask areas in and around each lesion that cannot be matched. However, areas within these regions do not then contribute to the computation of the deformable registration. Therefore, registration in and around the lesion, where it may be clinically most important, may be poor. Our method explicitly includes a geometric model of the pathology so that its deformations can be explicitly captured in conjunction with the deformations of the underlying image.

The basis of our work is that a registration method should be able to distinguish image appearance changes arising from the composition of (a) background deformations of the image with (b) foreground deformations of an embedded geometric object, e.g., a tumor. We model these transformations through a large displacement diffeomorphic metric mapping (LDDMM) [5]. LDDMM, however, is an inexact matching registration algorithm that only allows for the deformation of the source image and not for a change of its appearance.

In order to explicitly model appearance changes, we augment the LDDMM registration model with an extra control to model foreground deformation of a geometric model, which induces image appearance changes through an image composition model.

Specifically, we define the geometric metamorphosis problem as the minimization of:

where I0 is the source image that is being warped to the target image I1 and I(1) represents the warped source image. T1 is the geometric model of the lesion in I0; T2 is the geometric model of the lesion in I1; and Iτ(1) represents the warped T1. v is the foreground (non-lesion) deformation field. vτ is the background (lesion) deformation field. Note that the geometric model T1 and its image Iτ are subject to both deformations, whereas the source image is only subjected to the background deformation. L is a differential operator that controls the spatial regularity of these deformation fields. σ controls the weighting of the image match term. w ∈ (0, 1) controls the trade-off between background and foreground deformations. Sim denotes a similarity measure of choice; we use the L2 distance measure. Two similarity terms are used to assure matching of (a) the regions that correspond in both images and (b) the geometric models. Ic(•, •, •) denotes the image composition model. Additional detail is given in the full version of our paper [2].

[1] Brett, M., Lee, A., Rorden, C., Ashburner, J.: Spatial normalization of brain images with focal lesions using cost function
     masking. Neuroimage 14(2), 486-500(2001)
[2] Niethammer, M., Hart, G., Pace, D., Vespa, P., Irimia, A., Van Horn, J., and Aylward, S., Geometric Metamorphosis. To
     Appear, Medical Image Computing and Computer-Assisted Interventions, Toronto, September, 2011
[3] Hogea, C., Davatzikos, C., Biros, G.: An image-driven parameter estimation problem for a reaction-diffusion glioma growth
     model with mass effects. Journal of Mathematical Biology 56(6), 793-825 (2008)
[4] Lamecker, H., Pennec, X.: Atlas to image-with-tumor registration based on demons and deformation inpainting. In: Proc.
     MICCAI Workshop on Computational Imaging Biomarkers for Tumors – From Qualitative to Quantitative (CIBT’2010) (2010)
[5] Miller, M.I., Younes, L.: Group actions, homeomorphisms, and matching: A general framework. International Journal of
     Computer Vision 41(1/2), 61-84 (2001) [6] Prastawa, M., Bullitt, E., Gerig, G.: Simulation of brain tumors in MR images for
     evaluation of segmentation efficacy. Medical image analysis 13(2), 297-311 (2009)
[7] Stefanescu, R., Commowick, O., Malandain, G., Bondiau, P., Ayache, N., Pennec, X.: Non-rigid atlas to subject registration
     with pathologies for conformal brain radiotherapy. In: Barillot,C., Haynor, D.R., Hellier, P. (eds.) Medical Image Computing
     and Computer-Assisted Intervention, MICCAI 2004. pp. 704-711 (2004)

Stephen Aylward is Kitware’s Director of Medical Imaging Research. He is also an adjunct professor at UNC, an Associate Editor of IEEE transactions on Medical Imaging, and PI for a wide range of NIH-funded research efforts.



Gabe Hart joined Kitware in June, 2010 as an R&D engineer in the medical imaging group. Prior to joining, he studied medical image analysis at UNC, Chapel Hill. He splits his time between the ITKVideo and SimpleITK groups in the ITK v4.0 team and the NAMIC algorithms research team.



Marc Niethammer is an Assistant Professor at the University of North Carolina in the Department of Computer Science and the Biomedical Research Imaging Center. He received his Ph.D. in Electrical and Computer Engineering from the Georgia Institute of Technology and post-doctoral training at Harvard Medical School/Brigham and Women’s Hospital.

Danielle Pace is a Research and Development Engineer in Kitware’s medical imaging team. She joined Kitware’s North Carolina office in July 2010, and works on both research in deformable image registration and software engineering for 3D Slicer and Kitware’s commercial consulting projects.



Jack Van Horn is an Assistant Professor of Neurology at UCLA. He completed his doctoral work at the University of London, he has over 20 years of experience in the field of human neuroimaging, neurology, cognitive neuroscience, and large-scale data analysis and visualization.


Questions or comments are always welcome!