Simulation of brain tumors in MR images for evaluation of segmentation efficacy

Obtaining validation data and comparison metrics for segmentation of magnetic resonance images (MRI) are difficult tasks due to the lack of reliable ground truth. This problem is even more evident for images presenting pathology, which can both alter tissue appearance through infiltration and cause geometric distortions. Systems for generating synthetic images with user-defined degradation by noise and intensity inhomogeneity offer the possibility for testing and comparison of segmentation methods. Such systems do not yet offer simulation of sufficiently realistic looking pathology. This paper presents a system that combines physical and statistical modeling to generate synthetic multi-modal 3D brain MRI with tumor and edema, along with the underlying anatomical ground truth, Main emphasis is placed on simulation of the major effects known for tumor MRI, such as contrast enhancement, local distortion of healthy tissue, infiltrating edema adjacent to tumors, destruction and deformation of fiber tracts, and multi-modal MRI contrast of healthy tissue and pathology. The new method synthesizes pathology in multi-modal MRI and diffusion tensor imaging (DTI) by simulating mass effect, warping and destruction of white matter fibers, and infiltration of brain tissues by tumor cells. We generate synthetic contrast enhanced MR images by simulating the accumulation of contrast agent within the brain. The appearance of the the brain tissue and tumor in MRI is simulated by synthesizing texture images from real MR images. The proposed method is able to generate synthetic ground truth and synthesized MR images with tumor and edema that exhibit comparable segmentation challenges to real tumor MRI. Such image data sets will find use in segmentation reliability studies, comparison and validation of different segmentation methods, training and teaching, or even in evaluating standards for tumor size like the RECIST criteria (response evaluation criteria in solid tumors).

[1]  M. Christian,et al.  Measuring response in solid tumors: unidimensional versus bidimensional measurement. , 1999, Journal of the National Cancer Institute.

[2]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.

[3]  Guy Marchal,et al.  Multimodality image registration by maximization of mutual information , 1997, IEEE Transactions on Medical Imaging.

[4]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[5]  R. Wasserman,et al.  A patient-specific in vivo tumor model. , 1996, Mathematical biosciences.

[6]  William J. Schroeder,et al.  The Visualization Toolkit , 2005, The Visualization Handbook.

[7]  James C. Gee,et al.  Spatial transformations of diffusion tensor magnetic resonance images , 2001, IEEE Transactions on Medical Imaging.

[8]  G. Johnson,et al.  Peritumoral diffusion tensor imaging of high-grade gliomas and metastatic brain tumors. , 2003, AJNR. American journal of neuroradiology.

[9]  K Masaoka,et al.  The finite element analysis of brain oedema associated with intracranial meningiomas. , 1990, Acta neurochirurgica. Supplementum.

[10]  V. A. Barker,et al.  Finite element solution of boundary value problems , 1984 .

[11]  Heinz-Otto Peitgen,et al.  A Framework for the Generation of Realistic Brain Tumor Phantoms and Applications , 2004, MICCAI.

[12]  V. Cristini,et al.  Nonlinear simulation of tumor necrosis, neo-vascularization and tissue invasion via an adaptive finite-element/level-set method , 2005, Bulletin of mathematical biology.

[13]  D. Louis Collins,et al.  Twenty New Digital Brain Phantoms for Creation of Validation Image Data Bases , 2006, IEEE Transactions on Medical Imaging.

[14]  Hervé Delingette,et al.  In Silico Tumor Growth: Application to Glioblastomas , 2004, MICCAI.

[15]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[16]  Keith D. Paulsen,et al.  Model-Updated Image-Guided Neurosurgery Using the Finite Element Method: Incorporation of the Falx Cerebri , 1999, MICCAI.

[17]  Guido Gerig,et al.  Effects of Healthy Aging Measured By Intracranial Compartment Volumes Using a Designed MR Brain Database , 2005, MICCAI.

[18]  B Horwitz,et al.  A mathematical model for vasogenic brain edema. , 1978, Advances in neurology.

[19]  M. Gurtin,et al.  An introduction to continuum mechanics , 1981 .

[20]  Lawrence O. Hall,et al.  Automatic tumor segmentation using knowledge-based techniques , 1998, IEEE Transactions on Medical Imaging.

[21]  M. Christian,et al.  [New guidelines to evaluate the response to treatment in solid tumors]. , 2000, Bulletin du cancer.

[22]  Josef Kittler,et al.  Proceedings of the 4th International Conference on Pattern Recognition , 1988 .

[23]  Guido Gerig,et al.  Level-set evolution with region competition: automatic 3-D segmentation of brain tumors , 2002, Object recognition supported by user interaction for service robots.

[24]  John H. Gilmore,et al.  Automatic segmentation of MR images of the developing newborn brain , 2005, Medical Image Anal..

[25]  Ron Kikinis,et al.  Registration of 3D Intraoperative MR Images of the Brain Using a Finite Element Biomechanical Model , 2001, IEEE Trans. Medical Imaging.

[26]  Nicholas Ayache,et al.  Fast and Simple Calculus on Tensors in the Log-Euclidean Framework , 2005, MICCAI.

[27]  Guido Gerig,et al.  Valmet: A New Validation Tool for Assessing and Improving 3D Object Segmentation , 2001, MICCAI.

[28]  Keith D. Paulsen,et al.  Model-updated image guidance: initial clinical experiences with gravity-induced brain deformation , 1999, IEEE Transactions on Medical Imaging.

[29]  Dinggang Shen,et al.  Deformable Registration of Brain Tumor Images Via a Statistical Model of Tumor-Induced Deformation , 2005, MICCAI.

[30]  Christos Davatzikos,et al.  Finite Element Modeling of Brain Tumor Mass-Effect from 3D Medical Images , 2005, MICCAI.

[31]  Guido Gerig,et al.  User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability , 2006, NeuroImage.

[32]  Hervé Delingette,et al.  Realistic simulation of the 3-D growth of brain tumors in MR images coupling diffusion with biomechanical deformation , 2005, IEEE Transactions on Medical Imaging.

[33]  Alan C. Evans,et al.  BrainWeb: Online Interface to a 3D MRI Simulated Brain Database , 1997 .

[34]  Karol Miller,et al.  Biomechanics of the brain for computer-integrated surgery. , 2010, Acta of bioengineering and biomechanics.

[35]  S I Rapoport.  A mathematical model for vasogenic brain edema. , 1978, Journal of theoretical biology.

[36]  William M. Wells,et al.  Simultaneous truth and performance level estimation (STAPLE): an algorithm for the validation of image segmentation , 2004, IEEE Transactions on Medical Imaging.

[37]  S. R. Jammalamadaka,et al.  Directional Statistics, I , 2011 .

[38]  Marc Levoy,et al.  Fast texture synthesis using tree-structured vector quantization , 2000, SIGGRAPH.

[39]  L. R. Dice.  Measures of the Amount of Ecologic Association Between Species , 1945 .

[40]  Guido Gerig,et al.  A brain tumor segmentation framework based on outlier detection , 2004, Medical Image Anal..

[41]  R. Kikinis,et al.  Automated segmentation of MR images of brain tumors. , 2001, Radiology.

[42]  Christos Davatzikos,et al.  Nonlinear elastic registration of brain images with tumor pathology using a biomechanical model [MRI] , 1999, IEEE Transactions on Medical Imaging.

[43]  S. Pizer,et al.  Measuring tortuosity of the intracerebral vasculature from MRA images , 2003, IEEE Transactions on Medical Imaging.