Topology-preserving tissue classification of magnetic resonance brain images

SpringLS: a deformable model representation to provide interoperability between meshes and level sets

A new type of deformable model is presented that merges meshes and level sets into one representation to provide interoperability between methods designed for either. The key idea is to use a constellation of triangular surface elements (springls) to define a level set. A Spring Level Set (SpringLS) can be interpreted as a mesh or level set and used in place of them in many instances. There is no loss of shape information in the transformation from triangle mesh or level set into SpringLS. As examples, we present results for joint segmentation/spherical mapping of a human brain cortex and atlas/non-atlas segmentation of a pelvis.

A compressed sensing approach for MR tissue contrast synthesis

The tissue contrast of a magnetic resonance (MR) neuroimaging data set has a major impact on image analysis tasks like registration and segmentation. It has been one of the core challenges of medical imaging to guarantee the consistency of these tasks regardless of the contrasts of the MR data. Inconsistencies in image analysis are attributable in part to variations in tissue contrast, which in turn arise from operator variations during image acquisition as well as software and hardware differences in the MR scanners. It is also a common problem that images with a desired tissue contrast are completely missing in a given data set for reasons of cost, acquisition time, forgetfulness, or patient comfort. Absence of this data can hamper the detailed, automatic analysis of some or all data sets in a scientific study. A method to synthesize missing MR tissue contrasts from available acquired images using an atlas containing the desired contrast and a patch-based compressed sensing strategy is described. An important application of this general approach is to synthesize a particular tissue contrast from multiple studies using a single atlas, thereby normalizing all data sets into a common intensity space. Experiments on real data, obtained using different scanners and pulse sequences, show improvement in segmentation consistency, which could be extremely valuable in the pooling of multi-site multi-scanner neuroimaging studies.

Topology preserving warping of 3-D binary images according to continuous one-to-one mappings

The estimation of one-to-one mappings is one of the most intensively studied topics in the research field of nonrigid registration. Although the computation of such mappings can be now accurately and efficiently performed, the solutions for using them in the context of binary image deformation is much less satisfactory. In particular, warping a binary image with such transformations may alter its discrete topological properties if common resampling strategies are considered. In order to deal with this issue, this paper proposes a method for warping such images according to continuous and bijective mappings while preserving their discrete topological properties (i.e., their homotopy type). Results obtained in the context of the atlas-based segmentation of complex anatomical structures highlight the advantages of the proposed approach.

Topological correction of brain surface meshes using spherical harmonics

Surface reconstruction methods allow advanced analysis of structural and functional brain data beyond what can be achieved using volumetric images alone. Automated generation of cortical surface meshes from 3D brain MRI often leads to topological defects and geometrical artifacts that must be corrected to permit subsequent analysis. Here, we propose a novel method to repair topological defects using a surface reconstruction that relies on spherical harmonics. First, during reparameterization of the surface using a tiled platonic solid, the original MRI intensity values are used as a basis to select either a "fill" or "cut" operation for each topological defect. We modify the spherical map of the uncorrected brain surface mesh, such that certain triangles are favored while searching for the bounding triangle during reparameterization. Then, a low-pass filtered alternative reconstruction based on spherical harmonics is patched into the reconstructed surface in areas that previously contained defects. Self-intersections are repaired using a local smoothing algorithm that limits the number of affected points to less than 0.1% of the total, and as a last step, all modified points are adjusted based on the T1 intensity. We found that the corrected reconstructions have reduced distance error metrics compared with a "gold standard" surface created by averaging 12 scans of the same brain. Ninety-three percent of the topological defects in a set of 10 scans of control subjects were accurately corrected. The entire process takes 6-8 min of computation time. Further improvements are discussed, especially regarding the use of the T1-weighted image to make corrections.

Simple paradigm for extra-cerebral tissue removal: algorithm and analysis

Extraction of the brain-i.e. cerebrum, cerebellum, and brain stem-from T1-weighted structural magnetic resonance images is an important initial step in neuroimage analysis. Although automatic algorithms are available, their inconsistent handling of the cortical mantle often requires manual interaction, thereby reducing their effectiveness. This paper presents a fully automated brain extraction algorithm that incorporates elastic registration, tissue segmentation, and morphological techniques which are combined by a watershed principle, while paying special attention to the preservation of the boundary between the gray matter and the cerebrospinal fluid. The approach was evaluated by comparison to a manual rater, and compared to several other leading algorithms on a publically available data set of brain images using the Dice coefficient and containment index as performance metrics. The qualitative and quantitative impact of this initial step on subsequent cortical surface generation is also presented. Our experiments demonstrate that our approach is quantitatively better than six other leading algorithms (with statistical significance on modern T1-weighted MR data). We also validated the robustness of the algorithm on a very large data set of over one thousand subjects, and showed that it can replace an experienced manual rater as preprocessing for a cortical surface extraction algorithm with statistically insignificant differences in cortical surface position.

Topological well-composedness and glamorous glue: a digital gluing algorithm for topologically constrained front propagation

We propose a new approach to front propagation algorithms based on a topological variant of well-composedness which contrasts with previous methods based on simple point detection. This provides for a theoretical justification, based on the digital Jordan separation theorem, for digitally "gluing" evolved well-composed objects separated by well-composed curves or surfaces. Additionally, our framework can be extended to more relaxed topologically constrained algorithms based on multisimple points. For both methods this framework has the additional benefit of obviating the requirement for both a user-specified connectivity and a topologically-consistent marching cubes/squares algorithm in meshing the resulting segmentation.

Distance Transforms in Multi Channel MR Image Registration

Deformable registration techniques play vital roles in a variety of medical imaging tasks such as image fusion, segmentation, and post-operative surgery assessment. In recent years, mutual information has become one of the most widely used similarity metrics for medical image registration algorithms. Unfortunately, as a matching criteria, mutual information loses much of its effectiveness when there is poor statistical consistency and a lack of structure. This is especially true in areas of images where the intensity is homogeneous and information is sparse. Here we present a method designed to address this problem by integrating distance transforms of anatomical segmentations as part of a multi-channel mutual information framework within the registration algorithm. Our method was tested by registering real MR brain data and comparing the segmentation of the results against that of the target. Our analysis showed that by integrating distance transforms of the the white matter segmentation into the registration, the overall segmentation of the registration result was closer to the target than when the distance transform was not used.

COMPRESSED SENSING BASED INTENSITY NON-UNIFORMITY CORRECTION

We present a compressed sensing based approach to remove gain field from magnetic resonance (MR) images of the human brain. During image acquisition, the inhomogeneity present in the radio-frequency (RF) coil appears as shading artifact in the intensity image. The inhomogeneity poses problem in any automatic algorithm that uses intensity as a feature. It has been shown that at low field strength, the shading can be assumed to be a smooth field that is composed of low frequency components. Thus most inhomogeneity correction algorithms assume some kind of explicit smoothness criteria on the field. This sometimes limits the performance of the algorithms if the actual inhomogeneity is not smooth, which is the case at higher field strength. We describe a model-free, non-parametric patch-based approach that uses compressed sensing for the correction. We show that these features enable our algorithm to perform comparably with a current state of the art method N3 on images acquired at low field, while outperforming N3 when the image has non-smooth inhomogeneity, such as 7T images.

Multi-parametric neuroimaging reproducibility: a 3-T resource study

Modern MRI image processing methods have yielded quantitative, morphometric, functional, and structural assessments of the human brain. These analyses typically exploit carefully optimized protocols for specific imaging targets. Algorithm investigators have several excellent public data resources to use to test, develop, and optimize their methods. Recently, there has been an increasing focus on combining MRI protocols in multi-parametric studies. Notably, these have included innovative approaches for fusing connectivity inferences with functional and/or anatomical characterizations. Yet, validation of the reproducibility of these interesting and novel methods has been severely hampered by the limited availability of appropriate multi-parametric data. We present an imaging protocol optimized to include state-of-the-art assessment of brain function, structure, micro-architecture, and quantitative parameters within a clinically feasible 60-min protocol on a 3-T MRI scanner. We present scan-rescan reproducibility of these imaging contrasts based on 21 healthy volunteers (11 M/10 F, 22-61 years old). The cortical gray matter, cortical white matter, ventricular cerebrospinal fluid, thalamus, putamen, caudate, cerebellar gray matter, cerebellar white matter, and brainstem were identified with mean volume-wise reproducibility of 3.5%. We tabulate the mean intensity, variability, and reproducibility of each contrast in a region of interest approach, which is essential for prospective study planning and retrospective power analysis considerations. Anatomy was highly consistent on structural acquisition (~1-5% variability), while variation on diffusion and several other quantitative scans was higher (~<10%). Some sequences are particularly variable in specific structures (ASL exhibited variation of 28% in the cerebral white matter) or in thin structures (quantitative T2 varied by up to 73% in the caudate) due, in large part, to variability in automated ROI placement. The richness of the joint distribution of intensities across imaging methods can be best assessed within the context of a particular analysis approach as opposed to a summary table. As such, all imaging data and analysis routines have been made publicly and freely available. This effort provides the neuroimaging community with a resource for optimization of algorithms that exploit the diversity of modern MRI modalities. Additionally, it establishes a baseline for continuing development and optimization of multi-parametric imaging protocols.

Unsupervised MRI segmentation of brain tissues using a local linear model and level set

Real-world magnetic resonance imaging of the brain is affected by intensity nonuniformity (INU) phenomena which makes it difficult to fully automate the segmentation process. This difficult task is accomplished in this work by using a new method with two original features: (1) each brain tissue class is locally modeled using a local linear region representative, which allows us to account for the INU in an implicit way and to more accurately position the region's boundaries; and (2) the region models are embedded in the level set framework, so that the spatial coherence of the segmentation can be controlled in a natural way. Our new method has been tested on the ground-truthed Internet Brain Segmentation Repository (IBSR) database and gave promising results, with Tanimoto indexes ranging from 0.61 to 0.79 for the classification of the white matter and from 0.72 to 0.84 for the gray matter. To our knowledge, this is the first time a region-based level set model has been used to perform the segmentation of real-world MRI brain scans with convincing results.

Longitudinal changes in cortical thickness associated with normal aging

Imaging studies of anatomic changes in regional gray matter volumes and cortical thickness have documented age effects in many brain regions, but the majority of such studies have been cross-sectional investigations of individuals studied at a single point in time. In this study, using serial imaging assessments of participants in the Baltimore Longitudinal Study of Aging (BLSA), we investigate longitudinal changes in cortical thickness during aging in a cohort of 66 older adults (mean age 68.78; sd. 6.6; range 60-84 at baseline) without dementia. We used the Cortical Reconstruction Using Implicit Surface Evolution CRUISE suite of algorithms to automatically generate a reconstruction of the cortical surface and identified twenty gyral based regions of interest per hemisphere. Using mixed effects regression, we investigated longitudinal changes in these regions over a mean follow-up interval of 8 years. The main finding in this study is that age-related decline in cortical thickness is widespread, but shows an anterior-posterior gradient with frontal and parietal regions, in general, exhibiting greater rates of decline than temporal and occipital. There were fewer regions in the right hemisphere showing statistically significant age-associated longitudinal decreases in mean cortical thickness. Males showed greater rates of decline in the middle frontal, inferior parietal, parahippocampal, postcentral, and superior temporal gyri in the left hemisphere, right precuneus and bilaterally in the superior parietal and cingulate regions. Significant nonlinear changes over time were observed in the postcentral, precentral, and orbitofrontal gyri on the left and inferior parietal, cingulate, and orbitofrontal gyri on the right.

Topology-corrected segmentation and local intensity estimates for improved partial volume classification of brain cortex in MRI

In magnetic resonance imaging (MRI), accuracy and precision with which brain structures may be quantified are frequently affected by the partial volume (PV) effect. PV is due to the limited spatial resolution of MRI compared to the size of anatomical structures. Accurate classification of mixed voxels and correct estimation of the proportion of each pure tissue (fractional content) may help to increase the precision of cortical thickness estimation in regions where this measure is particularly difficult, such as deep sulci. The contribution of this work is twofold: on the one hand, we propose a new method to label voxels and compute tissue fractional content, integrating a mechanism for detecting sulci with topology preserving operators. On the other hand, we improve the computation of the fractional content of mixed voxels using local estimation of pure tissue intensity means. Accuracy and precision were assessed using simulated and real MR data and comparison with other existing approaches demonstrated the benefits of our method. Significant improvements in gray matter (GM) classification and cortical thickness estimation were brought by the topology correction. The fractional content root mean squared error diminished by 6.3% (p<0.01) on simulated data. The reproducibility error decreased by 8.8% (p<0.001) and the Jaccard similarity measure increased by 3.5% on real data. Furthermore, compared with manually guided expert segmentations, the similarity measure was improved by 12.0% (p<0.001). Thickness estimation with the proposed method showed a higher reproducibility compared with the measure performed after partial volume classification using other methods.

MR CONTRAST SYNTHESIS FOR LESION SEGMENTATION

The magnetic resonance contrast of a neuroimaging data set has strong impact on the utility of the data in image analysis tasks, such as registration and segmentation. Lengthy acquisition times often prevent routine acquisition of multiple MR contrast images, and opportunities for detailed analysis using these data would seem to be irrevocably lost. This paper describes an example based approach which uses patch matching from a multiple contrast atlas with the intended goal of generating an alternate MR contrast image, thus effectively simulating alternative pulse sequences from one another. In this paper, we deal specifically with Fluid Attenuated Inversion Recovery (FLAIR) sequence generation from T1 and T2 pulse sequences. The applicability of this synthetic FLAIR for estimating white matter lesions segmentation is demonstrated.

Neonatal brain image segmentation in longitudinal MRI studies

In the study of early brain development, tissue segmentation of neonatal brain MR images remains challenging because of the insufficient image quality due to the properties of developing tissues. Among various brain tissue segmentation algorithms, atlas-based brain image segmentation can potentially achieve good segmentation results on neonatal brain images. However, their performances rely on both the quality of the atlas and the spatial correspondence between the atlas and the to-be-segmented image. Moreover, it is difficult to build a population atlas for neonates due to the requirement of a large set of tissue-segmented neonatal brain images. To combat these obstacles, we present a longitudinal neonatal brain image segmentation framework by taking advantage of the longitudinal data acquired at late time-point to build a subject-specific tissue probabilistic atlas. Specifically, tissue segmentation of the neonatal brain is formulated as two iterative steps of bias correction and probabilistic-atlas-based tissue segmentation, along with the longitudinal atlas reconstructed by the late time image of the same subject. The proposed method has been evaluated qualitatively through visual inspection and quantitatively by comparing with manual delineations and two population-atlas-based segmentation methods. Experimental results show that the utilization of a subject-specific probabilistic atlas can substantially improve tissue segmentation of neonatal brain images.

Digital Topology on Adaptive Octree Grids

The theory of digital topology is used in many different image processing and computer graphics algorithms. Most of the existing theories apply to uniform cartesian grids, and they are not readily extensible to new algorithms targeting at adaptive cartesian grids. This article provides a rigorous extension of the classical digital topology framework for adaptive octree grids, including the characterization of adjacency, connected components, and simple points. Motivating examples, proofs of the major propositions, and algorithm pseudocodes are provided.

Online resource for validation of brain segmentation methods

One key issue that must be addressed during the development of image segmentation algorithms is the accuracy of the results they produce. Algorithm developers require this so they can see where methods need to be improved and see how new developments compare with existing ones. Users of algorithms also need to understand the characteristics of algorithms when they select and apply them to their neuroimaging analysis applications. Many metrics have been proposed to characterize error and success rates in segmentation, and several datasets have also been made public for evaluation. Still, the methodologies used in analyzing and reporting these results vary from study to study, so even when studies use the same metrics their numerical results may not necessarily be directly comparable. To address this problem, we developed a web-based resource for evaluating the performance of skull-stripping in T1-weighted MRI. The resource provides both the data to be segmented and an online application that performs a validation study on the data. Users may download the test dataset, segment it using whichever method they wish to assess, and upload their segmentation results to the server. The server computes a series of metrics, displays a detailed report of the validation results, and archives these for future browsing and analysis. We applied this framework to the evaluation of 3 popular skull-stripping algorithms--the Brain Extraction Tool [Smith, S.M., 2002. Fast robust automated brain extraction. Hum. Brain Mapp. 17 (3),143-155 (Nov)], the Hybrid Watershed Algorithm [Ségonne, F., Dale, A.M., Busa, E., Glessner, M., Salat, D., Hahn, H.K., Fischl, B., 2004. A hybrid approach to the skull stripping problem in MRI. NeuroImage 22 (3), 1060-1075 (Jul)], and the Brain Surface Extractor [Shattuck, D.W., Sandor-Leahy, S.R., Schaper, K.A., Rottenberg, D.A., Leahy, R.M., 2001. Magnetic resonance image tissue classification using a partial volume model. NeuroImage 13 (5), 856-876 (May) under several different program settings. Our results show that with proper parameter selection, all 3 algorithms can achieve satisfactory skull-stripping on the test data.

A hybrid geometric-statistical deformable model for automated 3-D segmentation in brain MRI

We present a novel 3-D deformable model-based approach for accurate, robust, and automated tissue segmentation of brain MRI data of single as well as multiple magnetic resonance sequences. The main contribution of this study is that we employ an edge-based geodesic active contour for the segmentation task by integrating both image edge geometry and voxel statistical homogeneity into a novel hybrid geometric-statistical feature to regularize contour convergence and extract complex anatomical structures. We validate the accuracy of the segmentation results on simulated brain MRI scans of both single T1-weighted and multiple T1/T2/PD-weighted sequences. We also demonstrate the robustness of the proposed method when applied to clinical brain MRI scans. When compared to a current state-of-the-art region-based level-set segmentation formulation, our white matter and gray matter segmentation resulted in significantly higher accuracy levels with a mean improvement in Dice similarity indexes of 8.55% ( p < 0.0001) and 10.18% ( p < 0.0001), respectively.

Brain Tissue Segmentation of Neonatal MR Images Using a Longitudinal Subject-specific Probabilistic Atlas

Brain tissue segmentation of neonate MR images is a challenging task in study of early brain development, due to low signal contrast among brain tissues and high intensity variability especially in white matter. Among various brain tissue segmentation algorithms, the atlas-based segmentation techniques can potentially produce reasonable segmentation results on neonatal brain images. However, their performance on the population-based atlas is still limited due to the high variability of brain structures across different individuals. Moreover, it may be impossible to generate a reasonable probabilistic atlas for neonates without tissue segmentation samples. To overcome these limitations, we present a neonatal brain tissue segmentation method by taking advantage of the longitudinal data available in our study to establish a subject-specific probabilistic atlas. In particular, tissue segmentation of the neonatal brain is formulated as two iterative steps of bias correction and probabilistic atlas based tissue segmentation, along with the guidance of brain tissue segmentation resulted from the later time images of the same subject which serve as a subject-specific probabilistic atlas. The proposed method has been evaluated qualitatively through visual inspection and quantitatively by comparing with manual delineation results. Experimental results show that the utilization of a subject-specific probabilistic atlas can substantially improve tissue segmentation of neonatal brain images.

Topology preserving warping of binary images: application to atlas-based skull segmentation

Lots of works have been recently carried out in the field of non-rigid registration to ensure the estimation of one-to-one mappings. However, warping a binary image with such transformations may alter its discrete topological properties if common resampling strategies are considered. This paper proposes an original method for warping a binary image according to some continuous and bijective mapping, while preserving its discrete topological properties. Results obtained in the context of atlas-based segmentation highlight the interest of the approach. Indeed, the method has been successfully applied to the segmentation of skull structures from a database of 15 CT-scans, providing both geometrically and topologically satisfactory results.

Homeomorphic brain image segmentation with topological and statistical atlases

Atlas-based segmentation techniques are often employed to encode anatomical information for the delineation of multiple structures in magnetic resonance images of the brain. One of the primary challenges of these approaches is to efficiently model qualitative and quantitative anatomical knowledge without introducing a strong bias toward certain anatomical preferences when segmenting new images. This paper explores the use of topological information as a prior and proposes a segmentation framework based on both topological and statistical atlases of brain anatomy. Topology can be used to describe continuity of structures, as well as the relationships between structures, and is often a critical component in cortical surface reconstruction and deformation-based morphometry. Our method guarantees strict topological equivalence between the segmented image and the atlas, and relies only weakly on a statistical atlas of shape. Tissue classification and fast marching methods are used to provide a powerful and flexible framework to handle multiple image contrasts, high levels of noise, gain field inhomogeneities, and variable anatomies. The segmentation algorithm has been validated on simulated and real brain image data and made freely available to researchers. Our experiments demonstrate the accuracy and robustness of the method and the limited influence of the statistical atlas.

Statistical and topological atlas based brain image segmentation

This paper presents a new atlas-based segmentation framework for the delineation of major regions in magnetic resonance brain images employing an atlas of the global topological structure as well as a statistical atlas of the regions of interest. A segmentation technique using fast marching methods and tissue classification is proposed that guarantees strict topological equivalence between the segmented image and the atlas. Experimental validation on simulated and real brain images shows that the method is accurate and robust.

Topology correction of segmented medical images using a fast marching algorithm

We present here a new method for correcting the topology of objects segmented from medical images. Whereas previous techniques alter a surface obtained from a binary segmentation of the object, our technique can be applied directly to the image intensities of a probabilistic or fuzzy segmentation, thereby propagating the topology for all isosurfaces of the object. From an analysis of topological changes and critical points in implicit surfaces, we derive a topology propagation algorithm that enforces any desired topology using a fast marching technique. The method has been applied successfully to the correction of the cortical gray matter/white matter interface in segmented brain images and is publicly released as a software plug-in for the MIPAV package.

Digital homeomorphisms in deformable registration

A common goal in deformable registration applications is to produce a spatial transformation that is diffeomorphic, thereby preserving the topology of structures being transformed. Because this constraint is typically enforced only on the continuum, however, topological changes can still occur within discretely sampled images. This work discusses the notion of homeomorphisms in digital images, and how it differs from the diffeomorphic/homeomorphic concepts in continuous spaces commonly used in medical imaging. We review the differences and problems brought by considering functions defined on a discrete grid, and propose a practical criterion for enforcing digital homeomorphisms in the context of atlas-based segmentation.