Recent contributions to the body of knowledge on traumatic brain injury (TBI) favor the view that multimodal neuroimaging using structural and functional magnetic resonance imaging (MRI and fMRI, respectively) as well as diffusion tensor imaging (DTI) has excellent potential to identify novel biomarkers and predictors of TBI outcome. This is particularly the case when such methods are appropriately combined with volumetric/morphometric analysis of brain structures and with the exploration of TBI-related changes in brain network properties at the level of the connectome. In this context, our present review summarizes recent developments on the roles of these two techniques in the search for novel structural neuroimaging biomarkers that have TBI outcome prognostication value. The themes being explored cover notable trends in this area of research, including (1) the role of advanced MRI processing methods in the analysis of structural pathology, (2) the use of brain connectomics and network analysis to identify outcome biomarkers, and (3) the application of multivariate statistics to predict outcome using neuroimaging metrics. The goal of the review is to draw the community's attention to these recent advances on TBI outcome prediction methods and to encourage the development of new methodologies whereby structural neuroimaging can be used to identify biomarkers of TBI outcome.
A method for automated location of shape differences in diseased anatomical structures via high resolution biomedical atlases annotated with labels from formal ontologies is described. In particular, a high resolution magnetic resonance image of the myocardium of the human left ventricle was segmented and annotated with structural terms from an extracted subset of the Foundational Model of Anatomy ontology. The atlas was registered to the end systole template of a previous study of left ventricular remodeling in cardiomyopathy using a diffeomorphic registration algorithm. The previous study used thresholding and visual inspection to locate a region of statistical significance which distinguished patients with ischemic cardiomyopathy from those with nonischemic cardiomyopathy. Using semantic technologies and the deformed annotated atlas, this location was more precisely found. Although this study used only a cardiac atlas, it provides a proof-of-concept that ontologically labeled biomedical atlases of any anatomical structure can be used to automate location-based inferences.
In a previous cross-sectional study on baseline data, we demonstrated that the volume of brain white matter hyperintensities (WMH) in the splenium of corpus callosum (SCC) predicted the current mobility function of older persons. The primary aim of this follow-up study was to determine the relation of WMH volume change in SCC (SCC-∆WMH) with change in mobility measures. A secondary aim was to characterize the global and regional progression of WMH. Mobility function and WMH burden were evaluated at baseline and at 2 years in 77 community-dwelling individuals (baseline age, 82 ± 4). Regional WMH in SCC, as well as genu and body of corpus callosum, subregions of corona radiata, and superior longitudinal fasciculus were determined using a white matter parcellation atlas. The total WMH volume increased 3.3 ± 3.5 ml/year, mainly through enlargement. Significant WMH increases were observed in all selected regions, particularly within the corona radiata. While at baseline and follow-up we observed correlations between WMH burden and several measures of mobility, longitudinal change correlated only with change in chair rise (CR). SCC-∆WMH showed the highest correlation (r = -0.413, p = 0.0002) and was the best regional predictor of CR decline (OR = 1.5, r(2) = 0.3). The SCC-∆WMH was more than five times larger in the CR-decline group compared to the no-decline group (p = 0.0003). The SCC-∆WMH (top quartile) showed a higher sensitivity/specificity for CR decline compared to change in total WMH, 63/88% versus 52/84%, respectively. The findings suggest that accrual of WMHs in posterior areas of the brain supporting inter-hemispheric integration and processing of visual-spatial information is a mechanism contributing to age-related mobility deterioration.
The XCEDE (XML-based Clinical and Experimental Data Exchange) XML schema, developed by members of the BIRN (Biomedical Informatics Research Network), provides an extensive metadata hierarchy for storing, describing and documenting the data generated by scientific studies. Currently at version 2.0, the XCEDE schema serves as a specification for the exchange of scientific data between databases, analysis tools, and web services. It provides a structured metadata hierarchy, storing information relevant to various aspects of an experiment (project, subject, protocol, etc.). Each hierarchy level also provides for the storage of data provenance information allowing for a traceable record of processing and/or changes to the underlying data. The schema is extensible to support the needs of various data modalities and to express types of data not originally envisioned by the developers. The latest version of the XCEDE schema and manual are available from http://www.xcede.org/ .
The nonlocal means (NLM) filter has become a popular approach for denoising medical images due to its excellent performance. However, its heavy computational load has been an important shortcoming preventing its use. NLM works by averaging pixels in nonlocal vicinities, weighting them depending on their similarity with the pixel of interest. This similarity is assessed based on the squared difference between corresponding pixels inside local patches centered at the locations compared. Our proposal is to reduce the computational load of this comparison by checking only a subset of salient features associated to the pixels, which suffice to estimate the actual difference as computed in the original NLM approach. The speedup achieved with respect to the original implementation is over one order of magnitude, and, when compared to more recent NLM improvements for MRI denoising, our method is nearly twice as fast. At the same time, we evidence from both synthetic and in vivo experiments that computing of appropriate salient features make the estimation of NLM weights more robust to noise. Consequently, we are able to improve the outcomes achieved with recent state of the art techniques for a wide range of realistic Signal-to-Noise ratio scenarios like diffusion MRI. Finally, the statistical characterization of the features computed allows to get rid of some of the heuristics commonly used for parameter tuning.
We propose a novel probabilistic framework to merge information from diffusion weighted imaging tractography and resting-state functional magnetic resonance imaging correlations to identify connectivity patterns in the brain. In particular, we model the interaction between latent anatomical and functional connectivity and present an intuitive extension to population studies. We employ the EM algorithm to estimate the model parameters by maximizing the data likelihood. The method simultaneously infers the templates of latent connectivity for each population and the differences in connectivity between the groups. We demonstrate our method on a schizophrenia study. Our model identifies significant increases in functional connectivity between the parietal/posterior cingulate region and the frontal lobe and reduced functional connectivity between the parietal/posterior cingulate region and the temporal lobe in schizophrenia. We further establish that our model learns predictive differences between the control and clinical populations, and that combining the two modalities yields better results than considering each one in isolation.
Functional MRI studies have uncovered a number of brain areas that demonstrate highly specific functional patterns. In the case of visual object recognition, small, focal regions have been characterized with selectivity for visual categories such as human faces. In this paper, we develop an algorithm that automatically learns patterns of functional specificity from fMRI data in a group of subjects. The method does not require spatial alignment of functional images from different subjects. The algorithm is based on a generative model that comprises two main layers. At the lower level, we express the functional brain response to each stimulus as a binary activation variable. At the next level, we define a prior over sets of activation variables in all subjects. We use a Hierarchical Dirichlet Process as the prior in order to learn the patterns of functional specificity shared across the group, which we call functional systems, and estimate the number of these systems. Inference based on our model enables automatic discovery and characterization of dominant and consistent functional systems. We apply the method to data from a visual fMRI study comprised of 69 distinct stimulus images. The discovered system activation profiles correspond to selectivity for a number of image categories such as faces, bodies, and scenes. Among systems found by our method, we identify new areas that are deactivated by face stimuli. In empirical comparisons with previously proposed exploratory methods, our results appear superior in capturing the structure in the space of visual categories of stimuli.
For a given cognitive task such as language processing, the location of corresponding functional regions in the brain may vary across subjects relative to anatomy. We present a probabilistic generative model that accounts for such variability as observed in fMRI data. We relate our approach to sparse coding that estimates a basis consisting of functional regions in the brain. Individual fMRI data is represented as a weighted sum of these functional regions that undergo deformations. We demonstrate the proposed method on a language fMRI study. Our method identified activation regions that agree with known literature on language processing and established correspondences among activation regions across subjects, producing more robust group-level effects than anatomical alignment alone.
The overall goal of this research is the design of statistical atlas models that can be created from normal subjects, but may generalize to be applicable to abnormal brains. We present a new style of joint modeling of fMRI, DTI, and structural MRI. Motivated by the fact that a white matter tract and related cortical areas are likely to displace together in the presence of a mass lesion (brain tumor), in this work we propose a rotation and translation invariant model that represents the spatial relationship between fiber tracts and anatomic and functional landmarks. This landmark distance model provides a new basis for representation of fiber tracts and can be used for detection and prediction of fiber tracts based on landmarks. Our results indicate that the measured model is consistent across normal subjects, and thus suitable for atlas building. Our experiments demonstrate that the model is robust to displacement and missing data, and can be successfully applied to a small group of patients with mass lesions.
Least Squares (LS) and its minimum variance counterpart, Weighted Least Squares (WLS), have become very popular when estimating the Diffusion Tensor (DT), to the point that they are the standard in most of the existing software for diffusion MRI. They are based on the linearization of the Stejskal-Tanner equation by means of the logarithmic compression of the diffusion signal. Due to the Rician nature of noise in traditional systems, a certain bias in the estimation is known to exist. This artifact has been made patent through some experimental set-ups, but it is not clear how the distortion translates in the reconstructed DT, and how important it is when compared to the other source of error contributing to the Mean Squared Error (MSE) in the estimate, i.e. the variance. In this paper we propose the analytical characterization of log-Rician noise and its propagation to the components of the DT through power series expansions. We conclude that even in highly noisy scenarios the bias for log-Rician signals remains moderate when compared to the corresponding variance. Yet, with the advent of Parallel Imaging (pMRI), the Rician model is not always valid. We make our analysis extensive to a number of modern acquisition techniques through the study of a more general Non Central-Chi (nc-χ) model. Since WLS techniques were initially designed bearing in mind Rician noise, it is not clear whether or not they still apply to pMRI. An important finding in our work is that the common implementation of WLS is nearly optimal when nc-χ noise is considered. Unfortunately, the bias in the estimation becomes far more important in this case, to the point that it may nearly overwhelm the variance in given situations. Furthermore, we evidence that such bias cannot be removed by increasing the number of acquired gradient directions. A number of experiments have been conducted that corroborate our analytical findings, while in vivo data have been used to test the actual relevance of the bias in the estimation.
Extracting structure of interest from medical images is an important yet tedious work. Due to the image quality, the shape knowledge is widely used for assisting and constraining the segmentation process. In many previous works, shape knowledge was incorporated by first constructing a shape space from training cases, and then constraining the segmentation process to be within the learned shape space. However, such an approach has certain limitations due to the number of variations, eigen-shapemodes, that can be captured in the learned shape space. Moreover, small scale shape variances are usually overwhelmed by those in the large scale, and therefore the local shape information is lost. In this work, we present a multiscale representation for shapes with arbitrary topology, and a fully automatic method to segment the target organ/tissue from medical images using such multiscale shape information and local image features. First, we handle the problem of lacking eigen-shapemodes by providing a multiscale shape representation using the wavelet transform. Consequently, the shape variances existing in the training shapes captured by the statistical learning step are also represented at various scales. Note that by doing so, one can greatly enrich the eigen-shapemodes as well as capture small scale shape changes. Furthermore, in order to make full use of the training information, not only the shape but also the grayscale training images are utilized in a multi-atlas initialization procedure. By combining such initialization with the multiscale shape knowledge, we perform segmentation tests for challenging medical data sets where the target objects have low contrast and sharp corner structures, and demonstrate the statistically significant improvement obtained by employing such multiscale representation, in representing shapes as well as the overall shape based segmentation tasks.
Thalamo-cortical feedback loops play a key role in the processing and coordination of processing and integration of perceptual inputs and outputs, and disruption in this connection has long been hypothesized to contribute significantly to neuropsychological disturbances in schizophrenia. To test this hypothesis, we applied diffusion tensor tractography to 18 patients suffering schizophrenia and 20 control subjects. Fractional anisotropy (FA) was evaluated in the bilateral anterior and posterior limbs of the internal capsule, and correlated with clinical and neurocognitive measures. Patients diagnosed with schizophrenia showed significantly reduced FA bilaterally in the anterior but not the posterior limb of the internal capsule, compared with healthy control subjects. Lower FA correlated with lower scores on tests of declarative episodic memory in the patient group only. These findings suggest that disruptions, bilaterally, in thalamo-cortical connections in schizophrenia may contribute to disease-related impairment in the coordination of mnemonic processes of encoding and retrieval that are vital for efficient learning of new information.
The purpose of this study is to evaluate perfusion indices and pharmacokinetic parameters in solitary pulmonary nodules (SPNs). Thirty patients of 34 enrolled with SPNs (15-30 mm) were evaluated in this study. T1 and T2-weighted structural images and 2D turbo FLASH perfusion images were acquired with shallow free breathing. B-spline nonrigid image registration and optimization by χ² test against pharmacokinetic model curve were performed on dynamic contrast-enhanced MRI. This allowed voxel-by-voxel calculation of k(ep) , the rate constant for tracer transport to and from plasma and the extravascular extracellular space. Mean transit time, time-to-peak, initial slope, and maximum enhancement (E(max) ) were calculated from time-intensity curves fitted to a gamma variate function. After blinded data analysis, correlation with tissue histology from surgical resection or biopsy samples was performed. Histologic evaluation revealed 25 malignant and five benign SPNs. All benign SPNs had k(ep) < 1.0 min⁻¹. Nineteen of 25 (76%) malignant SPNs showed k(ep) > 1.0 min⁻¹. Sensitivity to diagnose malignant SPNs at a cutoff of k(ep) = 1.0 min⁻¹ was 76%, specificity was 100%, positive predictive value was 100%, negative predictive value was 45%, and accuracy was 80%. Of all indices studied, k(ep) was the most significant in differentiating malignant from benign SPNs.
This article reports large item effects in a study of computer-based learning of neuroanatomy. Outcome measures of the efficiency of learning, transfer of learning, and generalization of knowledge diverged by a wide margin across test items, with certain sets of items emerging as particularly difficult to master. In addition, the outcomes of comparisons between instructional methods changed with the difficulty of the items to be learned. More challenging items better differentiated between instructional methods. This set of results is important for two reasons. First, it suggests that instruction may be more efficient if sets of consistently difficult items are the targets of instructional methods particularly suited to them. Second, there is wide variation in the published literature regarding the outcomes of empirical evaluations of computer-based instruction. As a consequence, many questions arise as to the factors that may affect such evaluations. The present article demonstrates that the level of challenge in the material that is presented to learners is an important factor to consider in the evaluation of a computer-based instructional system.
Diffusion tensor imaging (DTI) constitutes the most used paradigm among the diffusion-weighted magnetic resonance imaging (DW-MRI) techniques due to its simplicity and application potential. Recently, real-time estimation in DW-MRI has deserved special attention, with several proposals aiming at the estimation of meaningful diffusion parameters during the repetition time of the acquisition sequence. Specifically focusing on DTI, the underlying model of the noise present in the acquired data is not taken into account, leading to a suboptimal estimation of the diffusion tensor. In this paper, we propose an optimal real-time estimation framework for DTI reconstruction in single-coil acquisitions. By including an online estimation of the time-changing noise variance associated to the acquisition process, the proposed method achieves the sequential best linear unbiased estimator. Results on both synthetic and real data show that our method outperforms those so far proposed, reaching the best performance of the existing proposals by processing a substantially lower number of diffusion images.
Active contours are very popular tools for video tracking and image segmentation. Parameterized contours are used due to their fast evolution and have become the method of choice in the Sobolev context. Unfortunately, these contours are not easily adaptable to topological changes, and they may sometimes develop undesirable loops, resulting in erroneous results. To solve such topological problems, one needs an algorithm for contour self-crossing detection. We propose a simple methodology via simple techniques from differential topology. The detection is accomplished by inspecting the total net change of a given contour's angle, without point sorting and plane sweeping. We discuss the efficient implementation of the algorithm. We also provide algorithms for locating crossings by angle considerations and by plotting the four-connected lines between the discrete contour points. The proposed algorithms can be added to any parametric active-contour model. We show examples of successful tracking in real-world video sequences by Sobolev active contours and the proposed algorithms and provide ideas for further research.
Extracting anatomical and functional significant structures renders one of the important tasks for both the theoretical study of the medical image analysis, and the clinical and practical community. In the past, much work has been dedicated only to the algorithmic development. Nevertheless, for clinical end users, a well designed algorithm with an interactive software is necessary for an algorithm to be utilized in their daily work. Furthermore, the software would better be open sourced in order to be used and validated by not only the authors but also the entire community. Therefore, the contribution of the present work is twofolds: first, we propose a new robust statistics based conformal metric and the conformal area driven multiple active contour framework, to simultaneously extract multiple targets from MR and CT medical imagery in 3D. Second, an open source graphically interactive 3D segmentation tool based on the aforementioned contour evolution is implemented and is publicly available for end users on multiple platforms. In using this software for the segmentation task, the process is initiated by the user drawn strokes (seeds) in the target region in the image. Then, the local robust statistics are used to describe the object features, and such features are learned adaptively from the seeds under a non-parametric estimation scheme. Subsequently, several active contours evolve simultaneously with their interactions being motivated by the principles of action and reaction-this not only guarantees mutual exclusiveness among the contours, but also no longer relies upon the assumption that the multiple objects fill the entire image domain, which was tacitly or explicitly assumed in many previous works. In doing so, the contours interact and converge to equilibrium at the desired positions of the desired multiple objects. Furthermore, with the aim of not only validating the algorithm and the software, but also demonstrating how the tool is to be used, we provide the reader reproducible experiments that demonstrate the capability of the proposed segmentation tool on several public available data sets.
Quantitative analysis has tremendous but mostly unrealized potential in healthcare to support objective and accurate interpretation of the clinical imaging. In 2008, the National Cancer Institute began building the Quantitative Imaging Network (QIN) initiative with the goal of advancing quantitative imaging in the context of personalized therapy and evaluation of treatment response. Computerized analysis is an important component contributing to reproducibility and efficiency of the quantitative imaging techniques. The success of quantitative imaging is contingent on robust analysis methods and software tools to bring these methods from bench to bedside. 3D Slicer is a free open-source software application for medical image computing. As a clinical research tool, 3D Slicer is similar to a radiology workstation that supports versatile visualizations but also provides advanced functionality such as automated segmentation and registration for a variety of application domains. Unlike a typical radiology workstation, 3D Slicer is free and is not tied to specific hardware. As a programming platform, 3D Slicer facilitates translation and evaluation of the new quantitative methods by allowing the biomedical researcher to focus on the implementation of the algorithm and providing abstractions for the common tasks of data communication, visualization and user interface development. Compared to other tools that provide aspects of this functionality, 3D Slicer is fully open source and can be readily extended and redistributed. In addition, 3D Slicer is designed to facilitate the development of new functionality in the form of 3D Slicer extensions. In this paper, we present an overview of 3D Slicer as a platform for prototyping, development and evaluation of image analysis tools for clinical research applications. To illustrate the utility of the platform in the scope of QIN, we discuss several use cases of 3D Slicer by the existing QIN teams, and we elaborate on the future directions that can further facilitate development and validation of imaging biomarkers using 3D Slicer.
The registration of a pair of point sets as well as the estimation of their pointwise correspondences is a challenging and important task in computer vision. In this paper, we present a method to estimate the diffeomorphic deformation, together with the pointwise correspondences, between a pair of point sets. Many of the registration problems are iteratively solved by estimating the correspondence, locally optimizing certain cost functionals over the rigid or similarity or affine transformation group, then estimating the correspondence again, and so on. This type of approach, however, is well-known to be susceptible to suboptimal local solutions. In this paper, we first adopt the perspective of treating the registration as a posterior estimation optimization problem and solve it accordingly via a particle-filtering framework. Second, within such a framework, the diffeomorphic registration is performed to correct the nonlinear deformation of the points. In doing so, we provide a solution less susceptible to local minima. We provide the experimental results, which include challenging medical data sets where the two point sets differ by 180 (°) rotation as well as local deformations, to highlight the algorithm's capability of robustly finding the more globally optimal solution for the registration task.
Individuals with 22q11.2 deletion syndrome (22q11.2DS) evince a 30% incidence of schizophrenia. We compared the white matter (WM) of 22q11.2DS patients without schizophrenia to a group of matched healthy controls using Tract-Based-Spatial-Statistics (TBSS). We found localized reduction of Fractional Anisotropy (FA) and Axial Diffusivity (AD; measure of axonal integrity) in WM underlying the left parietal lobe. No changes in Radial Diffusivity (RD; measure of myelin integrity) were observed. Of note, studies in chronic schizophrenia patients report reduced FA, no changes in AD, and increases in RD in WM. Our findings suggest different WM microstructural pathology in 22q11.2DS than in patients with schizophrenia.