Publications by Year: 2007

2007

RATIONALE AND OBJECTIVES: To perform a retrospective, quantitative assessment of the anatomic relationship between intra-axial, supratentorial, primary brain tumors, and adjacent white matter fiber tracts based on anatomic and diffusion tensor magnetic resonance imaging (MRI). We hypothesized that white matter infiltration may be common among different types of tumor. MATERIAL AND METHODS: Preoperative, anatomic (T1- and T2-weighted), and LINESCAN diffusion tensor MRI were obtained in 12 patients harboring supratentorial gliomas (World Health Organization [WHO] Grades II and III). The two imaging modalities were rigidly registered. The tumors were manually segmented from the T1- and T2-weighted MRI, and their volume calculated. A three-dimensional tractography was performed in each case. A second segmentation and volume measurement was performed on the tumor regions intersecting adjacent white matter fiber tracts. Statistical methods included summary statistics to examine the fraction of tumor volume infiltrating adjacent white matter. RESULTS: There were five patients with low-grade oligodendroglioma (WHO Grade II), one with low-grade mixed oligoastrocytoma (WHO Grade II), one with ganglioglioma, two with low-grade astrocytoma (WHO Grade II), and three with anaplastic astrocytoma (WHO Grade III). We identified white matter tracts infiltrated by tumor in all 12 cases. The median tumor volume (+/- standard deviation) in our patient population was 42.5 +/- 28.9 mL. The median tumor volume (+/- standard deviation) infiltrating white matter fiber tracts was 5.2 +/- 9.9 mL. The median percentage of tumor volume infiltrating white matter fiber tracts was 21.4% +/- 9.7%. CONCLUSIONS: The information provided by diffusion tensor imaging combined with anatomic MRI might be useful for neurosurgical planning and intraoperative guidance. Our results confirm previous reports that extensive white matter infiltration by primary brain tumors is a common occurrence. However, prospective, large population studies are required to definitively clarify this issue, and how infiltration relates to histologic tumor type, tumor size, and location.
Dimaio SP, Pieper S, Chinzei K, Hata N, Haker SJ, Kacher DF, Fichtinger G, Tempany CM, Kikinis R. Robot-assisted Needle Placement in Open MRI: System Architecture, Integration and Validation. Comput Aided Surg. 2007;12(1):15–24.
In prostate cancer treatment, there is a move toward targeted interventions for biopsy and therapy, which has precipitated the need for precise image-guided methods for needle placement. This paper describes an integrated system for planning and performing percutaneous procedures with robotic assistance under MRI guidance. A graphical planning interface allows the physician to specify the set of desired needle trajectories, based on anatomical structures and lesions observed in the patient’s registered pre-operative and pre-procedural MR images, immediately prior to the intervention in an open-bore MRI scanner. All image-space coordinates are automatically computed, and are used to position a needle guide by means of an MRI-compatible robotic manipulator, thus avoiding the limitations of the traditional fixed needle template. Automatic alignment of real-time intra-operative images aids visualization of the needle as it is manually inserted through the guide. Results from in-scanner phantom experiments are provided.
Sierra R, Dimaio SP, Wada J, Hata N, ekely G abor S, Kikinis R, Jolesz FA. Patient Specific Simulation and Navigation of Ventriculoscopic Interventions. Stud Health Technol Inform. 2007;125:433–5.
In this paper a comprehensive framework for pre-operative planning, procedural skill training, and intraoperative navigation is presented. The goal of this system is to integrate surgical simulation with surgical planning in order to improve the individual treatment of patients. Various surgical approaches and new, more complex procedures can be assessed using a safe and objective platform that will allow the physicians to explore and discuss possible risks and benefits prior to the intervention. A simulation environment extends the pre-operative planning in a natural way, as it allows for direct evaluation of the surgical approach envisioned for each case. In addition, by providing intraoperative navigation based on this simulation, surgeons can carry out the previously optimized plan with higher precision and greater confidence.
Von Spiczak J, Samset E, DiMaio S, Reitmayr G, Schmalstieg D, Burghart C, Kikinis R. Device connectivity for image-guided medical applications. Stud Health Technol Inform. 2007;125:482–4.
The integration of medical devices with software applications is crucial for image-guided medical applications. This work describes a general device interface that has been designed for high-frequency streaming of multi-modal events, thus providing maximum performance and flexibility for such applications. Several sample applications and performance tests are provided to demonstrate the usability of the concept.
Nain D, Haker S, Bobick A, Tannenbaum A. Multiscale 3-D shape representation and segmentation using spherical wavelets. IEEE Trans Med Imaging. 2007;26(4):598–618.
This paper presents a novel multiscale shape representation and segmentation algorithm based on the spherical wavelet transform. This work is motivated by the need to compactly and accurately encode variations at multiple scales in the shape representation in order to drive the segmentation and shape analysis of deep brain structures, such as the caudate nucleus or the hippocampus. Our proposed shape representation can be optimized to compactly encode shape variations in a population at the needed scale and spatial locations, enabling the construction of more descriptive, nonglobal, nonuniform shape probability priors to be included in the segmentation and shape analysis framework. In particular, this representation addresses the shortcomings of techniques that learn a global shape prior at a single scale of analysis and cannot represent fine, local variations in a population of shapes in the presence of a limited dataset. Specifically, our technique defines a multiscale parametric model of surfaces belonging to the same population using a compact set of spherical wavelets targeted to that population. We further refine the shape representation by separating into groups wavelet coefficients that describe independent global and/or local biological variations in the population, using spectral graph partitioning. We then learn a prior probability distribution induced over each group to explicitly encode these variations at different scales and spatial locations. Based on this representation, we derive a parametric active surface evolution using the multiscale prior coefficients as parameters for our optimization procedure to naturally include the prior for segmentation. Additionally, the optimization method can be applied in a coarse-to-fine manner. We apply our algorithm to two different brain structures, the caudate nucleus and the hippocampus, of interest in the study of schizophrenia. We show: 1) a reconstruction task of a test set to validate the expressiveness of our multiscale prior and 2) a segmentation task. In the reconstruction task, our results show that for a given training set size, our algorithm significantly improves the approximation of shapes in a testing set over the Point Distribution Model, which tends to oversmooth data. In the segmentation task, our validation shows our algorithm is computationally efficient and outperforms the Active Shape Model algorithm, by capturing finer shape details.
Melonakos J, Gao Y, Tannenbaum A. Tissue Tracking: Applications for Brain MRI Classification. Proc SPIE Int Soc Opt Eng. 2007;6512.
Bayesian classification methods have been extensively used in a variety of image processing applications, including medical image analysis. The basic procedure is to combine data-driven knowledge in the likelihood terms with clinical knowledge in the prior terms to classify an image into a pre-determined number of classes. In many applications, it is difficult to construct meaningful priors and, hence, homogeneous priors are assumed. In this paper, we show how expectation-maximization weights and neighboring posterior probabilities may be combined to make intuitive use of the Bayesian priors. Drawing upon insights from computer vision tracking algorithms, we cast the problem in a tissue tracking framework. We show results of our algorithm on the classification of gray and white matter along with surrounding cerebral spinal fluid in brain MRI scans. We show results of our algorithm on 20 brain MRI datasets along with validation against expert manual segmentations.
Yu P, Grant E, Qi Y, Han X, Segonne F, Pienaar R, Busa E, Pacheco J, Makris N, Buckner RL, Golland P, Fischl B. Cortical surface shape analysis based on spherical wavelets. IEEE Trans Med Imaging. 2007;26(4):582–97.
In vivo quantification of neuroanatomical shape variations is possible due to recent advances in medical imaging and has proven useful in the study of neuropathology and neurodevelopment. In this paper, we apply a spherical wavelet transformation to extract shape features of cortical surfaces reconstructed from magnetic resonance images (MRIs) of a set of subjects. The spherical wavelet transformation can characterize the underlying functions in a local fashion in both space and frequency, in contrast to spherical harmonics that have a global basis set. We perform principal component analysis (PCA) on these wavelet shape features to study patterns of shape variation within normal population from coarse to fine resolution. In addition, we study the development of cortical folding in newborns using the Gompertz model in the wavelet domain, which allows us to characterize the order of development of large-scale and finer folding patterns independently. Given a limited amount of training data, we use a regularization framework to estimate the parameters of the Gompertz model to improve the prediction performance on new data. We develop an efficient method to estimate this regularized Gompertz model based on the Broyden-Fletcher-Goldfarb-Shannon (BFGS) approximation. Promising results are presented using both PCA and the folding development model in the wavelet domain. The cortical folding development model provides quantitative anatomic information regarding macroscopic cortical folding development and may be of potential use as a biomarker for early diagnosis of neurologic deficits in newborns.
Martin-Fernandez M, Alberola-López C, Ruiz-Alzola J, Westin CF. Sequential anisotropic Wiener filtering applied to 3D MRI data. Magn Reson Imaging. 2007;25(2):278–92.
We present three different sequential Wiener filters, namely, isotropic, orientation and anisotropic. The first one is similar to the classical Wiener filter in the sense that it uses an isotropic neighborhood to estimate its parameters. Here we present a sequential version of it. The orientation Wiener filter uses oriented neighborhoods to estimate the structure orientation present at each voxel, giving rise to a modified estimator of the parameters. Finally, the anisotropic Wiener filter combines both approaches adaptively so that the appropriate approach is locally selected. Several synthetic experiments are presented showing the performance of the filters with respect to their parameters. A mean square error analysis is performed using a publicly available magnetic resonance imaging (MRI) brain phantom and a comparison with other filtering approaches is carried out. In addition, results from filtering real MRI data are presented.
Onitsuka T, McCarley RW, Kuroki N, Dickey CC, Kubicki M, Demeo SS, Frumin M, Kikinis R, Jolesz FA, Shenton ME. Occipital lobe gray matter volume in male patients with chronic schizophrenia: A quantitative MRI study. Schizophr Res. 2007;92(1-3):197–206.
Schizophrenia is characterized by deficits in cognition as well as visual perception. There have, however, been few magnetic resonance imaging (MRI) studies of the occipital lobe as an anatomically defined region of interest in schizophrenia. To examine whether or not patients with chronic schizophrenia show occipital lobe volume abnormalities, we measured gray matter volumes for both the primary visual area (PVA) and the visual association areas (VAA) using MRI based neuroanatomical landmarks and three-dimensional information. PVA and VAA gray matter volumes were measured using high-spatial resolution MRI in 25 male patients diagnosed with chronic schizophrenia and in 28 male normal controls. Chronic schizophrenia patients showed reduced bilateral VAA gray matter volume (11%), compared with normal controls, whereas patients showed no group difference in PVA gray matter volume. These results suggest that reduced bilateral VAA may be a neurobiological substrate of some of the deficits observed in early visual processing in schizophrenia.