Statistical morphological skull stripping of adult and infant MRI data

Methods of Brain Extraction from Magnetic Resonance Images of Human Head: A Review

Medical images are providing vital information to aid physicians in diagnosing a disease afflicting the organ of a human body. Magnetic resonance imaging is an important imaging modality in capturing the soft tissues of the brain. Segmenting and extracting the brain is essential in studying the structure and pathological condition of brain. There are several methods that are developed for this purpose. Researchers in brain extraction or segmentation need to know the current status of the work that have been done. Such an information is also important for improving the existing method to get more accurate results or to reduce the complexity of the algorithm. In this paper we review the classical methods and convolutional neural network-based deep learning brain extraction methods.

Performance Improvement in Brain Tumor Detection in MRI Images Using a Combination of Evolutionary Algorithms and Active Contour Method

The process of treating brain cancer depends on the experience and knowledge of the physician, which may be associated with eye errors or may vary from person to person. For this reason, it is important to utilize an automatic tumor detection algorithm to assist radiologists and physicians for brain tumor diagnosis. The aim of the present study is to automatically detect the location of the tumor in a brain MRI image with high accuracy. For this end, in the proposed algorithm, first, the skull is separated from the brain using morphological operators. The image is then segmented by six evolutionary algorithms, i.e., Particle Swarm Optimization (PSO), Artificial Bee Colony (ABC), Genetic Algorithm (GA), Differential Evolution (DE), Harmony Search (HS), and Gray Wolf Optimization (GWO), as well as two other frequently-used techniques in the literature, i.e., K-means and Otsu thresholding algorithms. Afterwards, the tumor area is isolated from the brain using the four features extracted from the main tumor. Evaluation of the segmented area revealed that the PSO has the best performance compared with the other approaches. The segmented results of the PSO are then used as the initial curve for the Active contour to precisely specify the tumor boundaries. The proposed algorithm is applied on fifty images with two different types of tumors. Experimental results on T1-weighted brain MRI images show a better performance of the proposed algorithm compared to other evolutionary algorithms, K-means, and Otsu thresholding methods.

Conventional and Deep Learning Methods for Skull Stripping in Brain MRI

Skull stripping in brain magnetic resonance volume has recently been attracting attention due to an increased demand to develop an efficient, accurate, and general algorithm for diverse datasets of the brain. Accurate skull stripping is a critical step for neuroimaging diagnostic systems because neither the inclusion of non-brain tissues nor removal of brain parts can be corrected in subsequent steps, which results in unfixed error through subsequent analysis. The objective of this review article is to give a comprehensive overview of skull stripping approaches, including recent deep learning-based approaches. In this paper, the current methods of skull stripping have been divided into two distinct groups—conventional or classical approaches, and convolutional neural networks or deep learning approaches. The potentials of several methods are emphasized because they can be applied to standard clinical imaging protocols. Finally, current trends and future developments are addressed giving special attention to recent deep learning algorithms.

A Novel Approach for Manual Segmentation of the Amygdala and Hippocampus in Neonate MRI

The gross anatomy of the infant brain at term is fairly similar to that of the adult brain, but structures are immature, and the brain undergoes rapid growth during the first 2 years of life. Neonate magnetic resonance (MR) images have different contrasts compared to adult images, and automated segmentation of brain magnetic resonance imaging (MRI) can thus be considered challenging as less software options are available. Despite this, most anatomical regions are identifiable and thus amenable to manual segmentation. In the current study, we developed a protocol for segmenting the amygdala and hippocampus in T2-weighted neonatal MR images. The participants were 31 healthy infants between 2 and 5 weeks of age. Intra-rater reliability was measured in 12 randomly selected MR images, where 6 MR images were segmented at 1-month intervals between the delineations, and another 6 MR images at 6-month intervals. The protocol was also tested by two independent raters in 20 randomly selected T2-weighted images, and finally with T1 images. Intraclass correlation coefficient (ICC) and Dice similarity coefficient (DSC) for intra-rater, inter-rater, and T1 vs. T2 comparisons were computed. Moreover, manual segmentations were compared to automated segmentations performed by iBEAT toolbox in 10 T2-weighted MR images. The intra-rater reliability was high ICC ≥ 0.91, DSC ≥ 0.89, the inter-rater reliabilities were satisfactory ICC ≥ 0.90, DSC ≥ 0.75 for hippocampus and DSC ≥ 0.52 for amygdalae. Segmentations for T1 vs. T2-weighted images showed high consistency ICC ≥ 0.90, DSC ≥ 0.74. The manual and iBEAT segmentations showed no agreement, DSC ≥ 0.39. In conclusion, there is a clear need to improve and develop the procedures for automated segmentation of infant brain MR images.

Discrete pre-processing step effects in registration-based pipelines, a preliminary volumetric study on T1-weighted images

Pre-processing MRI scans prior to performing volumetric analyses is common practice in MRI studies. As pre-processing steps adjust the voxel intensities, the space in which the scan exists, and the amount of data in the scan, it is possible that the steps have an effect on the volumetric output. To date, studies have compared between and not within pipelines, and so the impact of each step is unknown. This study aims to quantify the effects of pre-processing steps on volumetric measures in T1-weighted scans within a single pipeline. It was our hypothesis that pre-processing steps would significantly impact ROI volume estimations. One hundred fifteen participants from the OASIS dataset were used, where each participant contributed three scans. All scans were then pre-processed using a step-wise pipeline. Bilateral hippocampus, putamen, and middle temporal gyrus volume estimations were assessed following each successive step, and all data were processed by the same pipeline 5 times. Repeated-measures analyses tested for a main effects of pipeline step, scan-rescan (for MRI scanner consistency) and repeated pipeline runs (for algorithmic consistency). A main effect of pipeline step was detected, and interestingly an interaction between pipeline step and ROI exists. No effect for either scan-rescan or repeated pipeline run was detected. We then supply a correction for noise in the data resulting from pre-processing.

Methods on Skull Stripping of MRI Head Scan Images-a Review

The high resolution magnetic resonance (MR) brain images contain some non-brain tissues such as skin, fat, muscle, neck, and eye balls compared to the functional images namely positron emission tomography (PET), single photon emission computed tomography (SPECT), and functional magnetic resonance imaging (fMRI) which usually contain relatively less non-brain tissues. The presence of these non-brain tissues is considered as a major obstacle for automatic brain image segmentation and analysis techniques. Therefore, quantitative morphometric studies of MR brain images often require a preliminary processing to isolate the brain from extra-cranial or non-brain tissues, commonly referred to as skull stripping. This paper describes the available methods on skull stripping and an exploratory review of recent literature on the existing skull stripping methods.

Quantitative evaluation of robust skull stripping and tumor detection applied to axial MR images

To isolate the brain from non-brain tissues using a fully automatic method may be affected by the presence of radio frequency non-homogeneity of MR images (MRI), regional anatomy, MR sequences, and the subjects of the study. In order to automate the brain tumor (Glioblastoma) detection, we proposed a novel approach of skull stripping for axial slices derived from MRI. Then, the brain tumor was detected using multi-level threshold segmentation based on histogram analysis. Skull-stripping method, was applied by adaptive morphological operations approach. This is considered an empirical threshold by calculation of the area of brain tissue, iteratively. It was employed on the registration of non-contrast T1-weighted (T1-WI) and its corresponding fluid attenuated inversion recovery sequence. Then, we used multi-thresholding segmentation (MTS) method which is proposed by Otsu. We calculated the performance metrics based on the similarity coefficients for patients (n = 120) with tumor. The adaptive algorithm of skull stripping and MTS of segmented tumors were achieved efficient in preliminary results with 92 and 80 % of Dice similarity coefficient and 0.3 and 25.8 % of false negative rate, respectively. The adaptive skull stripping algorithm provides robust skull-stripping results, and the tumor area for medical diagnosis was determined by MTS.

Brain extraction in pediatric ADC maps, toward characterizing neuro-development in multi-platform and multi-institution clinical images

Apparent Diffusion Coefficient (ADC) maps can be used to characterize myelination and to detect abnormalities in the developing brain. However, given the normal variation in regional ADC with myelination, detection of abnormalities is difficult when based on visual assessment. Quantitative and automated analysis of pediatric ADC maps is thus desired but requires accurate brain extraction as the first step. Currently, most existing brain extraction methods are optimized for structural T1-weighted MR images of fully myelinated brains. Due to differences in age and image contrast, these approaches do not translate well to pediatric ADC maps. To address this problem, we present a multi-atlas brain extraction framework that has 1) specificity: designed and optimized specifically for pediatric ADC maps; 2) generality: applicable to multi-platform and multi-institution data, and to subjects at various neuro-developmental stages across the first 6 years of life; 3) accuracy: highly accurate compared to expert annotations; and 4) consistency: consistently accurate regardless of sources of data and ages of subjects. We show how we achieve these goals, via optimizing major components in a multi-atlas brain extraction framework, and via developing and evaluating new criteria for its atlas ranking component. Moreover, we demonstrate that these goals can be achieved with a fixed set of atlases and a fixed set of parameters, which opens doors for our optimized framework to be used in large-scale and multi-institution neuro-developmental and clinical studies. In a pilot study, we use this framework in a dataset containing scanner-generated ADC maps from 308 pediatric patients collected during the course of routine clinical care. Our framework leads to successful quantifications of the changes in whole-brain volumes and mean ADC values across the first 6 years of life.

Brain extraction based on locally linear representation-based classification

Brain extraction is an important procedure in brain image analysis. Although numerous brain extraction methods have been presented, enhancing brain extraction methods remains challenging because brain MRI images exhibit complex characteristics, such as anatomical variability and intensity differences across different sequences and scanners. To address this problem, we present a Locally Linear Representation-based Classification (LLRC) method for brain extraction. A novel classification framework is derived by introducing the locally linear representation to the classical classification model. Under this classification framework, a common label fusion approach can be considered as a special case and thoroughly interpreted. Locality is important to calculate fusion weights for LLRC; this factor is also considered to determine that Local Anchor Embedding is more applicable in solving locally linear coefficients compared with other linear representation approaches. Moreover, LLRC supplies a way to learn the optimal classification scores of the training samples in the dictionary to obtain accurate classification. The International Consortium for Brain Mapping and the Alzheimer's Disease Neuroimaging Initiative databases were used to build a training dataset containing 70 scans. To evaluate the proposed method, we used four publicly available datasets (IBSR1, IBSR2, LPBA40, and ADNI3T, with a total of 241 scans). Experimental results demonstrate that the proposed method outperforms the four common brain extraction methods (BET, BSE, GCUT, and ROBEX), and is comparable to the performance of BEaST, while being more accurate on some datasets compared with BEaST.

Knowledge-guided robust MRI brain extraction for diverse large-scale neuroimaging studies on humans and non-human primates

Accurate and robust brain extraction is a critical step in most neuroimaging analysis pipelines. In particular, for the large-scale multi-site neuroimaging studies involving a significant number of subjects with diverse age and diagnostic groups, accurate and robust extraction of the brain automatically and consistently is highly desirable. In this paper, we introduce population-specific probability maps to guide the brain extraction of diverse subject groups, including both healthy and diseased adult human populations, both developing and aging human populations, as well as non-human primates. Specifically, the proposed method combines an atlas-based approach, for coarse skull-stripping, with a deformable-surface-based approach that is guided by local intensity information and population-specific prior information learned from a set of real brain images for more localized refinement. Comprehensive quantitative evaluations were performed on the diverse large-scale populations of ADNI dataset with over 800 subjects (55∼90 years of age, multi-site, various diagnosis groups), OASIS dataset with over 400 subjects (18∼96 years of age, wide age range, various diagnosis groups), and NIH pediatrics dataset with 150 subjects (5∼18 years of age, multi-site, wide age range as a complementary age group to the adult dataset). The results demonstrate that our method consistently yields the best overall results across almost the entire human life span, with only a single set of parameters. To demonstrate its capability to work on non-human primates, the proposed method is further evaluated using a rhesus macaque dataset with 20 subjects. Quantitative comparisons with popularly used state-of-the-art methods, including BET, Two-pass BET, BET-B, BSE, HWA, ROBEX and AFNI, demonstrate that the proposed method performs favorably with superior performance on all testing datasets, indicating its robustness and effectiveness.

Brain extraction using the watershed transform from markers

Isolation of the brain from other tissue types in magnetic resonance (MR) images is an important step in many types of neuro-imaging research using both humans and animal subjects. The importance of brain extraction is well appreciated-numerous approaches have been published and the benefits of good extraction methods to subsequent processing are well known. We describe a tool-the marker based watershed scalper (MBWSS)-for isolating the brain in T1-weighted MR images built using filtering and segmentation components from the Insight Toolkit (ITK) framework. The key elements of MBWSS-the watershed transform from markers and aggressive filtering with large kernels-are techniques that have rarely been used in neuroimaging segmentation applications. MBWSS is able to reliably isolate the brain without expensive preprocessing steps, such as registration to an atlas, and is therefore useful as the first stage of processing pipelines. It is an informative example of the level of accuracy achievable without using priors in the form of atlases, shape models or libraries of examples. We validate the MBWSS using a publicly available dataset, a paediatric cohort, an adolescent cohort, intra-surgical scans and demonstrate flexibility of the approach by modifying the method to extract macaque brains.

Multiscale segmentation of the skull in MR images for MRI-based attenuation correction of combined MR/PET

BACKGROUND AND OBJECTIVE

Combined magnetic resonance/positron emission tomography (MR/PET) is a relatively new, hybrid imaging modality. MR-based attenuation correction often requires segmentation of the bone on MR images. In this study, we present an automatic segmentation method for the skull on MR images for attenuation correction in brain MR/PET applications.

MATERIALS AND METHODS

Our method transforms T1-weighted MR images to the Radon domain and then detects the features of the skull image. In the Radon domain we use a bilateral filter to construct a multiscale image series. For the repeated convolution we increase the spatial smoothing in each scale and make the width of the spatial and range Gaussian function doubled in each scale. Two filters with different kernels along the vertical direction are applied along the scales from the coarse to fine levels. The results from a coarse scale give a mask for the next fine scale and supervise the segmentation in the next fine scale. The use of the multiscale bilateral filtering scheme is to improve the robustness of the method for noise MR images. After combining the two filtered sinograms, the reciprocal binary sinogram of the skull is obtained for the reconstruction of the skull image.

RESULTS

This method has been tested with brain phantom data, simulated brain data, and real MRI data. For real MRI data the Dice overlap ratios are 92.2%±1.9% between our segmentation and manual segmentation.

CONCLUSIONS

The multiscale segmentation method is robust and accurate and can be used for MRI-based attenuation correction in combined MR/PET.

LABEL: pediatric brain extraction using learning-based meta-algorithm

Magnetic resonance imaging of pediatric brain provides valuable information for early brain development studies. Automated brain extraction is challenging due to the small brain size and dynamic change of tissue contrast in the developing brains. In this paper, we propose a novel Learning Algorithm for Brain Extraction and Labeling (LABEL) specially for the pediatric MR brain images. The idea is to perform multiple complementary brain extractions on a given testing image by using a meta-algorithm, including BET and BSE, where the parameters of each run of the meta-algorithm are effectively learned from the training data. Also, the representative subjects are selected as exemplars and used to guide brain extraction of new subjects in different age groups. We further develop a level-set based fusion method to combine multiple brain extractions together with a closed smooth surface for obtaining the final extraction. The proposed method has been extensively evaluated in subjects of three representative age groups, such as neonate (less than 2 months), infant (1-2 years), and child (5-18 years). Experimental results show that, with 45 subjects for training (15 neonates, 15 infant, and 15 children), the proposed method can produce more accurate brain extraction results on 246 testing subjects (75 neonates, 126 infants, and 45 children), i.e., at average Jaccard Index of 0.953, compared to those by BET (0.918), BSE (0.902), ROBEX (0.901), GCUT (0.856), and other fusion methods such as Majority Voting (0.919) and STAPLE (0.941). Along with the largely-improved computational efficiency, the proposed method demonstrates its ability of automated brain extraction for pediatric MR images in a large age range.

An accurate skull stripping method based on simplex meshes and histogram analysis for magnetic resonance images

Skull stripping methods are designed to eliminate the non-brain tissue in magnetic resonance (MR) brain images. Removal of non-brain tissues is a fundamental step in enabling the processing of brain MR images. The aim of this study is to develop an automatic accurate skull stripping method based on deformable models and histogram analysis. A rough-segmentation step is used to find the optimal starting point for the deformation and is based on thresholds and morphological operators. Thresholds are computed using comparisons with an atlas, and modeling by Gaussians. The deformable model is based on a simplex mesh and its deformation is controlled by the image local gray levels and the information obtained on the gray level modeling of the rough-segmentation. Our Simplex Mesh and Histogram Analysis Skull Stripping (SMHASS) method was tested on the following international databases commonly used in scientific articles: BrainWeb, Internet Brain Segmentation Repository (IBSR), and Segmentation Validation Engine (SVE). A comparison was performed against three of the best skull stripping methods previously published: Brain Extraction Tool (BET), Brain Surface Extractor (BSE), and Hybrid Watershed Algorithm (HWA). Performance was measured using the Jaccard index (J) and Dice coefficient (κ). Our method showed the best performance and differences were statistically significant (p<0.05): J=0.904 and κ=0.950 on BrainWeb; J=0.905 and κ=0.950 on IBSR; J=0.946 and κ=0.972 on SVE.

Automatic skull stripping in MRI based on morphological filters and fuzzy c-means segmentation

In this paper a new automatic skull stripping method for T1-weighted MR image of human brain is presented. Skull stripping is a process that allows to separate the brain from the rest of tissues. The proposed method is based on a 2D brain extraction making use of fuzzy c-means segmentation and morphological operators applied on transversal slices. The approach is extended to the 3D case, taking into account the result obtained from the preceding slice to solve the organ splitting problem. The proposed approach is compared with BET (Brain Extraction Tool) implemented in MRIcro software.

Design and construction of a brain phantom to simulate neonatal MR images

This paper presents the design and construction of a 3D digital neonatal neurocranial phantom and its application for the simulation of brain magnetic resonance (MR) images. Commonly used digital brain phantoms (e.g. BrainWeb) are based on the adult brain. With the growing interest in computer-aided methods for neonatal MR image processing, there is a growing demand a digital phantom and brain MR image simulator especially for the neonatal brains. This is due to the pronounced differences between adult and neonatal brains not only in terms of size but also, more importantly, in terms of geometrical proportions and the need to subdivide white matter into two different tissue types in neonates. Therefore the neonatal brain phantom created in the here presented work consists of 9 different tissue types: skin, fat, muscle, skull, dura mater, gray matter, myelinated white matter, nonmyelinated white matter and cerebrospinal fluid. Each voxel has a vector consisting of 9 components, one for each of these nine tissue types. This digital phantom can be used to map simulated magnetic resonance signal intensities resulting in simulated MR images of the newborns head. These images with controlled degradation of the image data present a representative, reproducible data set ideal for development and evaluation of neonatal MRI analysis methods, e.g. segmentation and registration algorithms.

Skull stripping based on region growing for magnetic resonance brain images

In this paper, we propose a new skull stripping method for T1-weighted magnetic resonance (MR) brain images. Skull stripping has played an important role in neuroimage research because it is a basic preliminary step in many clinical applications. The process of skull stripping can be challenging due to the complexity of the human brain, variable parameters of MR scanners, individual characteristics, etc. In this paper, we aim to develop a computationally efficient and robust method. In the proposed algorithm, after eliminating the background voxels with histogram analysis, two seed regions of the brain and non-brain regions were automatically identified using a mask produced by morphological operations. Then we expanded these seed regions with a 2D region growing algorithm based on general brain anatomy information. The proposed algorithm was validated using 56 volumes of human brain data and simulated phantom data with manually segmented masks. It was compared with two popular automated skull stripping methods: the brain surface extractor (BSE) and the brain extraction tool (BET). The experimental results showed that the proposed algorithm produced accurate and stable results against data sets acquired from various MR scanners and effectively addressed difficult problems such as low contrast and large anatomical connections between the brain and surrounding tissues. The proposed method was also robust against noise, RF, and intensity inhomogeneities.

A semi-automatic image segmentation method for extraction of brain volume from in vivo mouse head magnetic resonance imaging using Constraint Level Sets

In vivo magnetic resonance imaging (MRI) of mouse brain has been widely used to non-invasively monitor disease progression and/or therapeutic effects in murine models of human neurodegenerative disease. Segmentation of MRI to differentiate brain from non-brain tissue (usually referred to as brain extraction) is required for many MRI data processing and analysis methods, including coregistration, statistical parametric analysis, and mapping to brain atlas and histology. This paper presents a semi-automatic brain extraction technique based on a level set method with the incorporation of user-defined constraints. The constraints are derived from the prior knowledge of brain anatomy by defining brain boundary on orthogonal planes of the MRI. Constraints are incorporated in the level set method by spatially varying the weighting factors of the internal and external forces and modifying the image gradient (edge) map. Both two-dimensional multislice and three-dimensional versions of the brain extraction technique were developed and applied to MRI data with minimal brain/non-brain contrast T(1)-weighted (T(1)-wt) FLASH and maximized contrast T(2)-weighted (T(2)-wt) RARE. Results were evaluated by calculating the overlap measure (OM) between the automatically segmented and manually traced brain volumes. Results demonstrate that this technique accurately extracts the brain volume (mean OM=94%) and consistently outperformed the region growing method applied to the T(2)-wt RARE MRI (mean OM=81%). This method not only successfully extracts the mouse brain in low and high contrast MRI, but can also be used to segment other organs and tissues.

Ingress of blood-borne macrophages across the blood-brain barrier in murine HIV-1 encephalitis

Blood-borne macrophage ingress into brain in HIV-1 associated neurocognitive disorders governs the tempo of disease. We used superparamagnetic iron-oxide particles loaded into murine bone marrow-derived macrophages (BMM) injected intravenously into HIV-1 encephalitis mice to quantitatively assess BMM entry into diseased brain regions. Magnetic resonance imaging tests were validated by histological coregistration and enhanced image processing. The demonstration of robust BMM migration into areas of focal encephalitis provide 'proof of concept' for the use of MRI to monitor macrophage ingress into brain.