Publications by Year: 2014

Sridharan R, Dalca AV, Golland P. An Interactive Visualization Tool for Nipype Medical Image Computing Pipelines. Int Conf Med Image Comput Comput Assist Interv. Workshop on Interactive Medical Image Computing. 2014;17.
Dalca AV, Sridharan R, Rost NS, Golland P. tipiX: Rapid Visualization of Large Image Collections. Int Conf Med Image Comput Comput Assist Interv. Workshop on Interactive Medical Image Computing. 2014;14.
Zhu L, Kolesov I, Gao Y, Kikinis R, Tannenbaum A. An Effective Interactive Medical Image Segmentation Method using Fast GrowCut. Int Conf Med Image Comput Comput Assist Interv. Workshop on Interactive Methods. 2014;17 (WS).Abstract
Segmentation of anatomical structures in medical imagery is a key step in a variety of clinical applications. Designing a generic, automated method that works for various structures and imaging modalities is a daunting task. In this paper, we present an effective interactive segmentation method that reformulates the GrowCut algorithm as a clustering problem and computes a fast, approximate solution. The method is further improved by using an efficient updating scheme requiring only local computations when new user input becomes available, making it applicable to high resolution images. The algorithm may easily be included as a user-oriented software module in any number of available medical imaging/image processing platforms such as 3D Slicer. The efficiency and effectiveness of the algorithm are demonstrated through tests on several challenging data sets where it is also compared to standard GrowCut.
Zhu MICCAI WS 2014
Kolesov I, Zhu L, Karasev P, Tannenbaum A. A Control Framework for Interactive Deformable Image Registration. Int Conf Med Image Comput Comput Assist Interv. Workshop on Interactive Methods. 2014;17 (WS). Kolesov MICCAI 2014
Farhat N. Tutorial: Preparing Data for 3-D Printing using 3D Slicer. 2014.Abstract
This tutorial demonstrates how to prepare data for 3D printing using the open source software 3D Slicer. The following topics are highlighted in the tutorial: introduction to the 3D Slicer interface, loading data into 3D Slicer, volume rendering, cropping image volumes, creating label maps, creating surface models, and saving data in file formats appropriate for 3D printing.
3D Printing
Kahn T, Jolesz FA, Lewin JS. Proceedings of the 10th Interventional MRI Symposium. 10th Interventional MRI Symposium. 2014;10 :1-85. 2014 iMRI Symposium Proceedings
Kapur T, Tempany CM, Jolesz FA. Proceedings of the 7th Image Guided Therapy Workshop. Image Guided Therapy Workshop. 2014;7 :1-60. 2014 IGT Workshop Proceedings
Chenevert TL, Malyarenko DI, Newitt D, Li X, Jayatilake M, Tudorica A, Fedorov A, Kikinis R, Liu TT, Muzi M, et al. Errors in Quantitative Image Analysis due to Platform-Dependent Image Scaling. Transl Oncol. 2014;7 (1) :65-71.Abstract
PURPOSE: To evaluate the ability of various software (SW) tools used for quantitative image analysis to properly account for source-specific image scaling employed by magnetic resonance imaging manufacturers. METHODS: A series of gadoteridol-doped distilled water solutions (0%, 0.5%, 1%, and 2% volume concentrations) was prepared for manual substitution into one (of three) phantom compartments to create "variable signal," whereas the other two compartments (containing mineral oil and 0.25% gadoteriol) were held unchanged. Pseudodynamic images were acquired over multiple series using four scanners such that the histogram of pixel intensities varied enough to provoke variable image scaling from series to series. Additional diffusion-weighted images were acquired of an ice-water phantom to generate scanner-specific apparent diffusion coefficient (ADC) maps. The resulting pseudodynamic images and ADC maps were analyzed by eight centers of the Quantitative Imaging Network using 16 different SW tools to measure compartment-specific region-of-interest intensity. RESULTS: Images generated by one of the scanners appeared to have additional intensity scaling that was not accounted for by the majority of tested quantitative image analysis SW tools. Incorrect image scaling leads to intensity measurement bias near 100%, compared to nonscaled images. CONCLUSION: Corrective actions for image scaling are suggested for manufacturers and quantitative imaging community.
Savadjiev P, Whitford TJ, Hough ME, Clemm von Hohenberg C, Bouix S, Westin C-F, Shenton ME, Crow TJ, A J, Kubicki M. Sexually Dimorphic White Matter Geometry Abnormalities in Adolescent Onset Schizophrenia. Cereb Cortex. 2014;24 (5) :1389-96.Abstract

The normal human brain is characterized by a pattern of gross anatomical asymmetry. This pattern, known as the "torque", is associated with a sexual dimorphism: The male brain tends to be more asymmetric than that of the female. This fact, along with well-known sex differences in brain development (faster in females) and onset of psychosis (earlier with worse outcome in males), has led to the theory that schizophrenia is a disorder in which sex-dependent abnormalities in the development of brain torque, the correlate of the capacity for language, cause alterations in interhemispheric connectivity, which are causally related to psychosis (Crow TJ, Paez P, Chance SE. 2007. Callosal misconnectivity and the sex difference in psychosis. Int Rev Psychiatry. 19(4):449-457.). To provide evidence toward this theory, we analyze the geometry of interhemispheric white matter connections in adolescent-onset schizophrenia, with a particular focus on sex, using a recently introduced framework for white matter geometry computation in diffusion tensor imaging data (Savadjiev P, Kindlmann GL, Bouix S, Shenton ME, Westin CF. 2010. Local white geometry from diffusion tensor gradients. Neuroimage. 49(4):3175-3186.). Our results reveal a pattern of sex-dependent white matter geometry abnormalities that conform to the predictions of Crow's torque theory and correlate with the severity of patients' symptoms. To the best of our knowledge, this is the first study to associate geometrical differences in white matter connectivity with torque in schizophrenia.

Nguyen T, Cevidanes L, Paniagua B, Zhu H, Koerich L, De Clerck H. Use of shape correspondence analysis to quantify skeletal changes associated with bone-anchored Class III correction. Angle Orthod. 2014;84 (2) :329-36.Abstract
OBJECTIVE: To evaluate the three-dimensional (3D) skeletal changes in the mandibles of Class III patients treated with bone-anchored maxillary protraction using shape correspondence analysis. MATERIAL AND METHOD: Twenty-five consecutive patients with skeletal Class III who were between the ages of 9 and 13 years (mean age, 11.10 ± 1.1 years) were treated using Class III intermaxillary elastics and bilateral miniplates (two in the infrazygomatic crests of the maxilla and two in the anterior mandible). Cone-beam computed tomography (CBCT) was performed for each patient before initial loading (T1) and at 1 year out (T2). From the CBCT scans, 3D models were generated, registered on the anterior cranial base, and analyzed using 3D linear distances and vectors between corresponding point-based surfaces. RESULTS: Bone-anchored traction produced anteroposterior and vertical skeletal changes in the mandible. The novel application of Shape correspondence analysis showed vectors of mean (± standard deviation) distal displacement of the posterior ramus of 3.6 ± 1.4 mm, while the chin displaced backward by 0.5 ± 3.92 mm. The lower border of the mandible at the menton region was displaced downward by 2.6 ± 1.2 mm, and the lower border at the gonial region moved downward by 3.6 ± 1.4 mm. There was a downward and backward displacement around the gonial region with a mean closure of the gonial angle by 2.1°. The condyles were displaced distally by a mean of 2.6 ± 1.5 mm, and there were three distinct patterns for displacement: 44% backward, 40% backward and downward, and 16% backward and upward. CONCLUSION: This treatment approach induces favorable control of the mandibular growth pattern and can be used to treat patients with components of mandibular prognathism.
Veeraraghavan H, Miller JV. Faceted visualization of three dimensional neuroanatomy by combining ontology with faceted search. Neuroinformatics. 2014;12 (2) :245-59.Abstract
In this work, we present a faceted-search based approach for visualization of anatomy by combining a three dimensional digital atlas with an anatomy ontology. Specifically, our approach provides a drill-down search interface that exposes the relevant pieces of information (obtained by searching the ontology) for a user query. Hence, the user can produce visualizations starting with minimally specified queries. Furthermore, by automatically translating the user queries into the controlled terminology our approach eliminates the need for the user to use controlled terminology. We demonstrate the scalability of our approach using an abdominal atlas and the same ontology. We implemented our visualization tool on the opensource 3D Slicer software. We present results of our visualization approach by combining a modified Foundational Model of Anatomy (FMA) ontology with the Surgical Planning Laboratory (SPL) Brain 3D digital atlas, and geometric models specific to patients computed using the SPL brain tumor dataset.
Rathi Y, Pasternak O, Savadjiev P, Michailovich OV, Bouix S, Kubicki M, Westin C-F, Makris N, Shenton ME. Gray Matter Alterations in Early Aging: A Diffusion Magnetic Resonance Imaging Study. Hum Brain Mapp. 2014;35 (8) :3841-56.Abstract

Many studies have observed altered neurofunctional and structural organization in the aging brain. These observations from functional neuroimaging studies show a shift in brain activity from the posterior to the anterior regions with aging (PASA model), as well as a decrease in cortical thickness, which is more pronounced in the frontal lobe followed by the parietal, occipital, and temporal lobes (retrogenesis model). However, very little work has been done using diffusion MRI (dMRI) with respect to examining the structural tissue alterations underlying these neurofunctional changes in the gray matter. Thus, for the first time, we propose to examine gray matter changes using diffusion MRI in the context of aging. In this work, we propose a novel dMRI based measure of gray matter "heterogeneity" that elucidates these functional and structural models (PASA and retrogenesis) of aging from the viewpoint of diffusion MRI. In a cohort of 85 subjects (all males, ages 15-55 years), we show very high correlation between age and "heterogeneity" (a measure of structural layout of tissue in a region-of-interest) in specific brain regions. We examine gray matter alterations by grouping brain regions into anatomical lobes as well as functional zones. Our findings from dMRI data connects the functional and structural domains and confirms the "retrogenesis" hypothesis of gray matter alterations while lending support to the neurofunctional PASA model of aging in addition to showing the preservation of paralimbic areas during healthy aging.

Helmer KG, Pasternak O, Fredman E, Preciado RI, Koerte IK, Sasaki T, Mayinger M, Johnson AM, Holmes JD, Forwell LA, et al. Hockey Concussion Education Project, Part 1. Susceptibility-weighted imaging study in male and female ice hockey players over a single season. J Neurosurg. 2014;120 (4) :864-72.Abstract
OBJECT: Concussion, or mild traumatic brain injury (mTBI), is a commonly occurring sports-related injury, especially in contact sports such as hockey. Cerebral microbleeds (CMBs), which appear as small, hypointense lesions on T₂*-weighted images, can result from TBI. The authors use susceptibility-weighted imaging (SWI) to automatically detect small hypointensities that may be subtle signs of chronic and acute damage due to both subconcussive and concussive injury. The goal was to investigate how the burden of these hypointensities changes over time, over a playing season, and postconcussion, in comparison with subjects who did not suffer a medically observed and diagnosed concussion. METHODS: Images were obtained in 45 university-level adult male and female ice hockey players before and after a single Canadian Interuniversity Sports season. In addition, 11 subjects (5 men and 6 women) underwent imaging at 72 hours, 2 weeks, and 2 months after concussion. To identify subtle changes in brain tissue and potential CMBs, nonvessel clusters of hypointensities on SWI were automatically identified, and a hypointensity burden index was calculated for all subjects at the beginning of the season (BOS), the end of the season (EOS), and at postconcussion time points (where applicable). RESULTS: A statistically significant increase in the hypointensity burden, relative to the BOS, was observed for male subjects with concussions at the 2-week postconcussion time point. A smaller, nonsignificant rise in the burden for female subjects with concussions was also observed within the same time period. There were no significant changes in burden for nonconcussed subjects of either sex between the BOS and EOS time points. However, there was a statistically significant difference in the burden between male and female subjects in the nonconcussed group at both the BOS and EOS time points, with males having a higher burden. CONCLUSIONS: This method extends the utility of SWI from the enhancement and detection of larger (> 5 mm) CMBs, which are often observed in more severe cases of TBI, to cases involving smaller lesions in which visual detection of injury is difficult. The hypointensity burden metric proposed here shows statistically significant changes over time in the male subjects. A smaller, nonsignificant increase in the burden metric was observed in the female subjects.
Pasternak O, Koerte IK, Bouix S, Fredman E, Sasaki T, Mayinger M, Helmer KG, Johnson AM, Holmes JD, Forwell LA, et al. Hockey Concussion Education Project, Part 2. Microstructural white matter alterations in acutely concussed ice hockey players: a longitudinal free-water MRI study. J Neurosurg. 2014;120 (4) :873-81.Abstract
OBJECT: Concussion is a common injury in ice hockey and a health problem for the general population. Traumatic axonal injury has been associated with concussions (also referred to as mild traumatic brain injuries), yet the pathological course that leads from injury to recovery or to long-term sequelae is still not known. This study investigated the longitudinal course of concussion by comparing diffusion MRI (dMRI) scans of the brains of ice hockey players before and after a concussion. METHODS: The 2011-2012 Hockey Concussion Education Project followed 45 university-level ice hockey players (both male and female) during a single Canadian Interuniversity Sports season. Of these, 38 players had usable dMRI scans obtained in the preseason. During the season, 11 players suffered a concussion, and 7 of these 11 players had usable dMRI scans that were taken within 72 hours of injury. To analyze the data, the authors performed free-water imaging, which reflects an increase in specificity over other dMRI analysis methods by identifying alterations that occur in the extracellular space compared with those that occur in proximity to cellular tissue in the white matter. They used an individualized approach to identify alterations that are spatially heterogeneous, as is expected in concussions. RESULTS: Paired comparison of the concussed players before and after injury revealed a statistically significant (p < 0.05) common pattern of reduced free-water volume and reduced axial and radial diffusivities following elimination of free-water. These free-water-corrected measures are less affected by partial volumes containing extracellular water and are therefore more specific to processes that occur within the brain tissue. Fractional anisotropy was significantly increased, but this change was no longer significant following the free-water elimination. CONCLUSIONS: Concussion during ice hockey games results in microstructural alterations that are detectable using dMRI. The alterations that the authors found suggest decreased extracellular space and decreased diffusivities in white matter tissue. This finding might be explained by swelling and/or by increased cellularity of glia cells. Even though these findings in and of themselves cannot determine whether the observed microstructural alterations are related to long-term pathology or persistent symptoms, they are important nonetheless because they establish a clearer picture of how the brain responds to concussion.
Sasaki T, Pasternak O, Mayinger M, Muehlmann M, Savadjiev P, Bouix S, Kubicki M, Fredman E, Dahlben B, Helmer KG, et al. Hockey Concussion Education Project, Part 3. White matter microstructure in ice hockey players with a history of concussion: a diffusion tensor imaging study. J Neurosurg. 2014;120 (4) :882-90.Abstract
OBJECT: The aim of this study was to examine the brain's white matter microstructure by using MR diffusion tensor imaging (DTI) in ice hockey players with a history of clinically symptomatic concussion compared with players without a history of concussion. METHODS: Sixteen players with a history of concussion (concussed group; mean age 21.7 ± 1.5 years; 6 female) and 18 players without a history of concussion (nonconcussed group; mean age 21.3 ± 1.8 years, 10 female) underwent 3-T DTI at the end of the 2011-2012 Canadian Interuniversity Sports ice hockey season. Tract-based spatial statistics (TBSS) was used to test for group differences in fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD), and the measure "trace," or mean diffusivity. Cognitive evaluation was performed using the Immediate Postconcussion Assessment and Cognitive Test (ImPACT) and the Sport Concussion Assessment Tool-2 (SCAT2). RESULTS: TBSS revealed a significant increase in FA and AD, and a significant decrease in RD and trace in several brain regions in the concussed group, compared with the nonconcussed group (p < 0.05). The regions with increased FA and decreased RD and trace included the right posterior limb of the internal capsule, the right corona radiata, and the right temporal lobe. Increased AD was observed in a small area in the left corona radiata. The DTI measures correlated with neither the ImPACT nor the SCAT2 scores. CONCLUSIONS: The results of the current study indicate that a history of concussion may result in alterations of the brain's white matter microstructure in ice hockey players. Increased FA based on decreased RD may reflect neuroinflammatory or neuroplastic processes of the brain responding to brain trauma. Future studies are needed that include a longitudinal analysis of the brain's structure and function following a concussion to elucidate further the complex time course of DTI changes and their clinical meaning.
Garlapati RR, Roy A, Joldes GR, Wittek A, Mostayed A, Doyle B, Warfield SK, Kikinis R, Knuckey N, Bunt S, et al. More accurate neuronavigation data provided by biomechanical modeling instead of rigid registration. J Neurosurg. 2014;120 (6) :1477-83.Abstract
It is possible to improve neuronavigation during image-guided surgery by warping the high-quality preoperative brain images so that they correspond with the current intraoperative configuration of the brain. In this paper, the accuracy of registration results obtained using comprehensive biomechanical models is compared with the accuracy of rigid registration, the technology currently available to patients. This comparison allows investigation into whether biomechanical modeling provides good-quality image data for neuronavigation for a larger proportion of patients than rigid registration. Preoperative images for 33 neurosurgery cases were warped onto their respective intraoperative configurations using both the biomechanics-based method and rigid registration. The Hausdorff distance-based evaluation process, which measures the difference between images, was used to quantify the performance of both registration methods. A statistical test for difference in proportions was conducted to evaluate the null hypothesis that the proportion of patients for whom improved neuronavigation can be achieved is the same for rigid and biomechanics-based registration. The null hypothesis was confidently rejected (p < 10(-4)). Even the modified hypothesis that fewer than 25% of patients would benefit from the use of biomechanics-based registration was rejected at a significance level of 5% (p = 0.02). The biomechanics-based method proved particularly effective in cases demonstrating large craniotomy-induced brain deformations. The outcome of this analysis suggests that nonlinear biomechanics-based methods are beneficial to a large proportion of patients and can be considered for use in the operating theater as a possible means of improving neuronavigation and surgical outcomes.
Asami T, Hyuk Lee S, Bouix S, Rathi Y, Whitford TJ, Niznikiewicz M, Nestor P, McCarley RW, Shenton ME, Kubicki M. Cerebral white matter abnormalities and their associations with negative but not positive symptoms of schizophrenia. Psychiatry Res. 2014;222 (1-2) :52-9.Abstract
Although diffusion tensor imaging (DTI) studies have reported fractional anisotropy (FA) abnormalities in multiple white matter (WM) regions in schizophrenia, relationship between abnormal FA and negative symptoms has not been fully explored. DTI data were acquired from twenty-four patients with chronic schizophrenia and twenty-five healthy controls. Regional brain abnormalities were evaluated by conducting FA comparisons in the cerebral and each lobar WMs between groups. Focal abnormalities were also evaluated with a voxel-wise tract specific method. Associations between structural WM changes and negative symptoms were assessed using the Scale for the Assessment of Negative Symptoms (SANS). The patient group showed decreased FA in the cerebrum, especially in the frontal lobe, compared with controls. A voxel-wise analysis showed FA decreases in almost all WM tracts in schizophrenia. Correlation analyses demonstrated negative relationships between FA in the cerebrum, particularly in the left hemisphere, and SANS global and global rating scores (Anhedonia-Asociality, Attention, and Affective-Flattening), and also associations between FA of left frontal lobe and SANS global score, Anhedonia-Asociality, and Attention. This study demonstrates that patients with chronic schizophrenia evince widespread cerebral FA abnormalities and that these abnormalities, especially in the left hemisphere, are associated with negative symptoms.
Appia V, Yezzi A, Arepalli C, Faber T, Stillman A, Tannenbaum A. A complete system for automatic extraction of left ventricular myocardium from CT images using shape segmentation and contour evolution. IEEE Trans Image Process. 2014;23 (3) :1340-51.Abstract
The left ventricular myocardium plays a key role in the entire circulation system and an automatic delineation of the myocardium is a prerequisite for most of the subsequent functional analysis. In this paper, we present a complete system for an automatic segmentation of the left ventricular myocardium from cardiac computed tomography (CT) images using the shape information from images to be segmented. The system follows a coarse-to-fine strategy by first localizing the left ventricle and then deforming the myocardial surfaces of the left ventricle to refine the segmentation. In particular, the blood pool of a CT image is extracted and represented as a triangulated surface. Then, the left ventricle is localized as a salient component on this surface using geometric and anatomical characteristics. After that, the myocardial surfaces are initialized from the localization result and evolved by applying forces from the image intensities with a constraint based on the initial myocardial surface locations. The proposed framework has been validated on 34-human and 12-pig CT images, and the robustness and accuracy are demonstrated.
Naaz F, Chariker JH, Pani JR. Computer-Based Learning: Graphical Integration of Whole and Sectional Neuroanatomy Improves Long-Term Retention. Cogn Instr. 2014;32 (1) :44-64.Abstract
A study was conducted to test the hypothesis that instruction with graphically integrated representations of whole and sectional neuroanatomy is especially effective for learning to recognize neural structures in sectional imagery (such as MRI images). Neuroanatomy was taught to two groups of participants using computer graphical models of the human brain. Both groups learned whole anatomy first with a three-dimensional model of the brain. One group then learned sectional anatomy using two-dimensional sectional representations, with the expectation that there would be transfer of learning from whole to sectional anatomy. The second group learned sectional anatomy by moving a virtual cutting plane through the three-dimensional model. In tests of long-term retention of sectional neuroanatomy, the group with graphically integrated representation recognized more neural structures that were known to be challenging to learn. This study demonstrates the use of graphical representation to facilitate a more elaborated (deeper) understanding of complex spatial relations.
Crabb MG, Davidson JL, Little R, Wright P, Morgan AR, Miller CA, Naish JH, Parker GJM, Kikinis R, McCann H, et al. Mutual Information as a Measure of Image Quality for 3-D Dynamic Lung Imaging with EIT. Physiol Meas. 2014;35 (5) :863-79.Abstract
We report on a pilot study of dynamic lung electrical impedance tomography (EIT) at the University of Manchester. Low-noise EIT data at 100 frames per second were obtained from healthy male subjects during controlled breathing, followed by magnetic resonance imaging (MRI) subsequently used for spatial validation of the EIT reconstruction. The torso surface in the MR image and electrode positions obtained using MRI fiducial markers informed the construction of a 3D finite element model extruded along the caudal-distal axis of the subject. Small changes in the boundary that occur during respiration were accounted for by incorporating the sensitivity with respect to boundary shape into a robust temporal difference reconstruction algorithm. EIT and MRI images were co-registered using the open source medical imaging software, 3D Slicer. A quantitative comparison of quality of different EIT reconstructions was achieved through calculation of the mutual information with a lung-segmented MR image. EIT reconstructions using a linear shape correction algorithm reduced boundary image artefacts, yielding better contrast of the lungs, and had 10% greater mutual information compared with a standard linear EIT reconstruction.