# An automated and simple method for brain MR image extraction

- Haiyan Zhang†
^{1}, - Jiafeng Liu†
^{1}, - Zixin Zhu†
^{2}and - Haiyun Li
^{1}Email author

**10**:81

https://doi.org/10.1186/1475-925X-10-81

© Zhang et al; licensee BioMed Central Ltd. 2011

**Received: **16 June 2011

**Accepted: **13 September 2011

**Published: **13 September 2011

## Abstract

### Background

The extraction of brain tissue from magnetic resonance head images, is an important image processing step for the analyses of neuroimage data. The authors have developed an automated and simple brain extraction method using an improved geometric active contour model.

### Methods

The method uses an improved geometric active contour model which can not only solve the boundary leakage problem but also is less sensitive to intensity inhomogeneity. The method defines the initial function as a binary level set function to improve computational efficiency. The method is applied to both our data and Internet brain MR data provided by the Internet Brain Segmentation Repository.

### Results

The results obtained from our method are compared with manual segmentation results using multiple indices. In addition, the method is compared to two popular methods, Brain extraction tool and Model-based Level Set.

### Conclusions

The proposed method can provide automated and accurate brain extraction result with high efficiency.

## Background

Brain extraction which segments magnetic resonance (MR) head images into brain and non-brain region is often required for analyses of neuroimage data. Accurate and automated brain extraction plays an important role in the analyses because brain region should be isolated before other processing algorithms such as tissue classification, registration or cortical surface reconstruction can be made [1–3]. For example, functional images such as Functional magnetic resonance image (FMRI) and Positron Emission Tomography (PET) image usually contain few non-brain tissues, whereas high resolution MR images often contain some non-brain tissues(i.e., skin, fat, muscle, etc.), and if the non-brain tissues of MR images can be accurately removed beforehand, the registration robustness will be improved greatly [2]. Furthermore, as a pre-processing step, brain extraction is usually performed before a full segmentation of the brain region into grey matter (GM), white matter (WM), and cerebrospinal fluid (CSF), so that the segmentation problem can be simplified [4, 5]. On the other hand, brain extraction is also a difficult and time-consuming pre-processing step performed in neuroimage analysis due to the complexity of human brain anatomy and weak boundaries between brain and non-brain tissues.

Ma et al. [6] pointed out that researchers should combine the application background with practical requirements to design a proper algorithm for a medical image segmentation task. Although many brain extraction algorithms (BEAs) have been proposed to accurately segment brain from non-brain tissues, their segmentation quality varies greatly and has important influence on the results of subsequent image analysis. Boesen et al. [7] compared the performance of four BEAs and concluded that the brain extraction tool (BET) and the brain surface extractor (BSE) was significantly faster than the statistical parametric mapping(SPM) or Minneapolis consensus strip(McStrip). Compared to two different manual strip-masks, however, McStrip outperformed BET, SPM and BSE based on the Correct Boundary and Pertinent Boundary criteria and misclassified the least number of brain voxels. These popular methods have both their advantages and weaknesses, and none of them can be accurate and robust enough for large-scale neuroimage analysis [7, 8]. Subsequent research work are aimed at developing fully automated, accurate and robust BEAs for MR images. To facilitate large-scale neuroimage analysis, Zhuang et al. [9] developed a new automatic BEA called the model based level set method (MLS) which can provide robust performance for large-scale neuroimage analysis. As the number of subjects increases and the real-time image processing for clinical application develops, the need for fully automated, simple and fast brain extraction algorithms will become critical. In this paper, we first proposed a new method satisfying the requirement of both fully automated brain extraction and accurate brain extraction result; then we summarized the experimental results, evaluation, and comparison of our method to BET and MLS; finally, we discussed the advantages and disadvantages of our method for brain MR image extraction.

## Methods

We developed an accurate and simple brain extraction method using an improved geometric active contour model (GAC) which can not only solve the boundary leakage problem but also is less sensitive to intensity inhomogeneity. The proposed brain extraction method comprises three major steps: image intensity parameters are first estimated and a binary image of the head is calculated for the following segmentation procedures. Then the initial contour is automatically determined within the brain region. Finally, the proposed geometric active contour model is applied to extract the brain region on each of the slices.

### The proposed GAC model with a new local region-based signed pressure force function for brain MR image extraction

*ϕ(x, t)*be a 3-D scalar function whose zero level set defines the geometric active contour. The traditional geometric active contour formulation is written below [10, 11]:

*k*is the curvature and

*v*

_{ 0 }is a constant, and

*c(k+v*

_{ 0 }

*)*determines the overall evolution speed of level sets

*ϕ*along their normal direction. The curvature

*k*has the effect of smoothing the contour, while

*v*

_{ 0 }has the effect of shrinking or expanding contour at a constant speed. The speed of contour evolution is coupled with the image data through a multiplicative stopping term

*c*. This scheme works well for objects that have good contrast, but when it is used for detecting boundaries between brain and non/brain tissues, this contour tends to leak through the boundary or fail to reach the desired boundary. To solve the so-called boundary leaking problem, Zhang et al. [13] proposed a region-based active contour model with a region-based signed pressure force (SPF) function which can efficiently stop the contours at weak or blurred edges. The SPF function is also called the region function which is derived from the image and has values in the range [-1, 1]. The region function modulates the sign of the pressure forces using region information so that the contour shrinks when it is outside the object of interest or expands when it is inside the object. The model constructed by Zhang et al. [13] only used the image statistical information of the entire region inside and outside the contour, which can't successfully segment brain MR images with intensity inhomogeneity. However, intensity inhomogeneity often occurs in MR images due to technical limitations or motion artifacts and may cause difficulties in image segmentation. Li et al. [14] proposed a region-based local binary fitting (LBF) model which utilized the local intensity information in local regions to solve the problem caused by intensity inhomogeneity. To solve the boundary leakage and intensity inhomogeneity which were both common in brain MR image extraction, we constructed a new local region-based SPF function which utilized two smooth constants that locally approximated the image intensities inside and outside the contour as follows:

*G*

_{ σ }is the Gaussian kernel with standard deviation

*σ*and * denotes the convolution operator. In Eq.(3), We take the maximum absolute value of the numerator as the denominator, which can ensure the Eq.(3) has values in the range [-1,1], so Eq.(3) can serve as a SPF function.

*f*

_{ 1 }and

*f*

_{ 2 }are two smooth constants defined as follows:

*H(ϕ)*is the Heaviside function which is generally approximated by a smooth function

*H*

_{ ε }defined as follows:

*ε*is a positive constant.

*f*

_{ 1 }and

*f*

_{ 2 }of Eq.(4) can be regarded as the weighted averages of image intensities in a Gaussian window inside and outside the contour, respectively. Therefore, the proposed SPF function can not only utilize the locally statistical information inside and outside the contour to control the evolution but also use the smoothing effect of Gaussian filter, which is thus less sensitive to noise and intensity inhomogeneity. Substituting the SPF function in Eq.(3) for the edge potential function

*c(x)*in Eq.(1), the level set formulation of the proposed model for brain extraction is as follows:

The proposed geometric active contour model utilizes the new local region-based SPF function to solve the boundary leaking problem and intensity inhomogeneity which traditional geometric active contour models fail to solve. So our model not only works well for objects that have good contrast but also for objects with weak boundary such as brain surface.

### Estimation of image intensity parameters and binary image of the head

*t*

_{ 1 }

*, t*

_{ 2 }] is determined to ignore the voxels with unusual intensities, where

*t*

_{ 1 }

*, t*

_{ 2 }are the intensity values in the histogram chosen such that the accumulated number of voxels reaches 2% and 98%, respectively. Subsequently, a threshold

*t*is chosen empirically to separate the background and the skin or muscle tissue that covers the head which can be calculated as:

*T*c is a constant different for axial, coronal and sagittal orientation. The value of

*T*

_{ c }is determined according to image data empirically and sometimes several local thresholds may be needed to solve leakage problem. With the threshold

*t*, we can get the binary image of head as shown in Figure 1.(d), which is useful for the determination of initial contour.

### Automatic initialization

*ϕ*, is initialized as a signed distance function

*ϕ*

_{ 0 }for computational efficiency and re-initialization has been extensively used as a numerical remedy for maintaining stable curve evolution [12]. However, in practice, the re-initialization process is complicated and expensive. In this work, the initial function,

*ϕ*

_{ 0 }is defined as a binary level set function as in the work of Lie [16]. Such initialization is very simple to implement in practice and can greatly improve computational efficiency. The initializing method is as follows:

- 1.
Find the most left and right voxels automatically by searching all the voxels with the maximum intensity of the binary image of the head.

- 2.
Compute the distance

*d*_{ l-r }between the most left and right voxels. - 3.
Approximate the shape of the brain in axial as a square with its length set to be the distance

*d*_{ l-r }, and the shape of the brain in coronal or sagittal orientation as a rectangle with a 8:7 and 4:3 ratio of width (equal to the distance*d*_{ l-r }) to height, respectively. - 4.
Initialize the zero level set curve as a circle and the initial circle is positioned at the center of the square or rectangle with its radius, r, equal to one third of the length of the square side(for axial orientation) or one third of the width of the rectangle(for coronal and sagittal orientation).

- 5.To simplify the initialization of the level set function, we use a binary level set function as in the work of Lie [16]. Each level set function can only take two values at convergence, then the initial function
*ϕ*_{ 0 }is defined as:${\varphi}_{0}\left(x,\phantom{\rule{2.77695pt}{0ex}}y\right)=\left\{\begin{array}{c}-\rho \phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\left(x,\phantom{\rule{2.77695pt}{0ex}}y\right)\in \mathsf{\text{int}}\phantom{\rule{1em}{0ex}}\left({\Omega}_{1}\right)\\ \rho \phantom{\rule{1em}{0ex}}\phantom{\rule{1em}{0ex}}\phantom{\rule{2.77695pt}{0ex}}\phantom{\rule{2.77695pt}{0ex}}\left(x,\phantom{\rule{2.77695pt}{0ex}}y\right)\in out\left({\Omega}_{1}\right)\end{array}\right.\phantom{\rule{1em}{0ex}}$(8)

where *ρ* is a positive constant, *ϕ*
_{
0
} *(x, y)* denotes *ϕ(x, y, t)* at t = 0, Ω is the enclosing interface. With such initialization of the level set function, not only the re-initialization procedure is completely eliminated, but also the level set function *ϕ* is no longer required to be initialized as a signed distance function.

### Correction of leakage through weak boundaries based on local thresholds estimation

*Fcurv*to prevent high curvature, and lastly segmented the same slice again. However, the correction step is not completely effective and also time-consuming as shown in Figure 3.(h). In this work, we try to solve the problem with a higher efficiency. So we propose an method to correct effectively the leakage of weak boundaries based on local thresholds estimation. Our method to prevent such leakage comprises two major steps. First, we confirm if leakage occurs or not during segmenting process. By observing, we find such leakage through weak boundaries often occurs at some points with high curvature (shown in Figure 3.(c)), so we search such points and, if the number of these points reaches a pre-set value which is determined according to image data empirically, we estimate that segmentation leakage will occur, so the segmentation will be stopped immediately. Then, a new segmentation for the same slice will start again based on the local thresholds estimation. In local thresholds estimation, the brain region is divided into several parts with two different thresholds which is used to separate the background and other non-brain tissues as mentioned previously. As illustrated by Figure 3.(d), the thresholds of part I and part II are the same but higher than those of the other parts. The reason to make such improvement lies in the fact that, in some areas, the CSF is often thinner than other parts, leading to similar intensities of non-brain tissues with brain tissues and the leakage through weak boundaries is liable to occur. As a result, these parts are separated with a higher threshold, which can reveal more details of the weak boundary and lower the risk of the segmentation leakage. Experimental results show that local thresholds estimation can solve the boundary leaking problem effectively as shown in Figure 3.(f). Figure 3.(h), (i), and 3(j) also list the segmentation results with leakage of MLS, our segmentation results without leakage, and expert segmentation results, respectively.

## Results and Discussion

^{2}and the slice thickness is 1 mm. We used the following default setting of the parameters: σ = 1.5, time step Δ

*t*= 1,

*ρ*= 1,

*v*

_{ 0 }= 5. Figure 2 shows the final results on four sample images displayed in three orientations. To measure the extraction accuracy of our algorithm, 10 normal MRI brain data sets and the corresponding manual segmentations were obtained from the Internet Brain Segmentation Repository (IBSR) developed by the Centre for Morphometric Analysis (CMA) at Massachusetts General Hospital. Each volume has around 65 coronal slices, with 256 × 256 pixels per slice. The slice resolution is 1.02 × 1.04 mm

^{2}, and slice thickness is 3.1 mm. Figure 4 shows the final results of our method on eight sample volumes as well as manual segmentation results displayed in coronal orientation. We also computed the

*sensitivity, specificity, Jaccard*index,

*Dice*index and the

*FP_Rate*of our segmentation results using the manual segmentation results provided by the IBSR (shown in Table 1).

Performance comparison of BET, MLS and the proposed method for multiple indices using the IBSR data sets

Method | Sensitivity | Specificity | Jaccard | Dice | FP_rate |
---|---|---|---|---|---|

BET | 0.999(0.001) | 0.982(0.005) | 0.896(0.045) | 0.945(0.026) | 0.115(0.063) |

MLS | 0.982(0.03) | 0.991(0.008) | 0.925(0.041) | 0.961(0.022) | 0.069(0.055) |

Our method | 0.973(0.01) | 0.993(0.003) | 0.923(0.022) | 0.960(0.012) | 0.05(0.022) |

*FP_Rate*and

*Specificity*. Second, in our experiments, BET is always conservative and often includes some non-brain tissues, leading to the best sensitivity but worst specificity and

*FP_Rate*coefficients. The reason BET had such performance, perhaps because it was more important to avoid removing brain tissue than to remove all non-brain tissues for clinical application. Third, MLS is superior to the our method with regard to sensitivity and has similar performance on Jaccard and Dice indices. Generally speaking, when compared to BET, our method does not include too many non-brain tissues and need not tune many input parameter, so it is accurate and simple to use. Compared to MLS, due to our automatic initialization method, our method is more efficient. So as an automated and simple brain extraction tool, our method can accurately extract brain tissue with high efficiency.

## Conclusions

We proposed an automated and simple brain extraction method using an improved geometric active contour model. Our method has the following advantages over existing brain extraction algorithms: first, Our method uses a binary level set function to eliminate the expensive re-initialization of the existing brain extraction algorithms, it is thus more efficient. Second, the method not only utilizes the image statistical information to construct a new local region-based SPF function, but also corrects the leakage through extremely weak boundaries based on local thresholds estimation, thus can successfully segment brain tissue with weak boundaries. Third, the initial contour can be automatically set inside the brain with sufficiently large radius to improve the automation and the efficiency of the brain extraction. Last but most importantly, our method is very simple and easy to use. No preprocessing step is needed and all the results can be obtained using the original, noisy MR data. Thus, the proposed method can extract brain tissue with high efficiency and full automation compare to two other methods. However, our method was tested using normal adult MRI brain data sets only, and larger sample data sets including different age groups and abnormal anatomy structures such as tumor are needed in order to further test our method as a fully automated, simple and robust method for brain extraction.

## Notes

## Declarations

### Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 30670576) and Scientific Research Key Program of Beijing Municipal Commission of Education (Grant No. kz200810025011). The authors would also like to thank the anonymous reviewers for their constructive comments.

## Authors’ Affiliations

## References

- Shattuck DW, Sandor-leahy SR, Schaper KA, Rottenberg DA, Leahy RM:
**Magnetic resonance image tissue classification using a partial volume model.***NeuroImage*2001,**13:**856–876.View ArticleGoogle Scholar - Smith SM:
**Fast robust automated brain extraction.***Hum Brain Mapp*2002,**17:**143–155. 10.1002/hbm.10062View ArticleGoogle Scholar - Ségonna F, Dale AM, Busa E, Glessner M, Salat D, Hahn HK:
**A hybrid approach to the skull stripping problem in MRI.***NeuroImage*2004,**22:**1060–1075. 10.1016/j.neuroimage.2004.03.032View ArticleGoogle Scholar - Dale AM, Fischl B, Sereno MI:
**Cortical surface-based analysis I: segmentation and surface reconstruction.***NeuroImage*1999,**9:**179–194. 10.1006/nimg.1998.0395View ArticleGoogle Scholar - Zhang YY, Brady M, Smith S:
**Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm.***IEEE Trans Med Imaging*2001,**20:**45–57. 10.1109/42.906424View ArticleGoogle Scholar - Ma Z, Tavares JMRS, Jorge RN, Mascarenhas T:
**Review of Algorithms for Medical Image Segmentation and their Applications to the Female Pelvic Cavity.***Comput Meth Biomech Biomed Eng*2010,**13**(2):235–246. 10.1080/10255840903131878View ArticleGoogle Scholar - Boesen K, Rehm K, Schaper K, Stoltzner S, Woods R, Lϋders E, Rottenberg D:
**Quantitative comparison of four brain extraction algorithms.***NeuroImage*2004,**22**(3):1255–1261. 10.1016/j.neuroimage.2004.03.010View ArticleGoogle Scholar - Rehm K, Schaper K, Anderson J, Woods R, Stoltzner S, Rottenberg D:
**Putting our heads together: a consensus approach to brain/non-brain segmentation in T1-weighted MR volumes.***NeuroImage*2004,**22**(3):1262–1270. 10.1016/j.neuroimage.2004.03.011View ArticleGoogle Scholar - Zhuang AH, Valentino DJ, Toga A:
**Skull-stripping magnetic resonance brain images using a model-based level set.***NeuroImage*2006,**32:**79–92. 10.1016/j.neuroimage.2006.03.019View ArticleGoogle Scholar - Caselles V, Catte F, Coll T, Dibos F:
**A geometric model for active Contours in image processing.***Numerische Mathematik*1993,**66**(1):1–31. 10.1007/BF01385685MATHMathSciNetView ArticleGoogle Scholar - Malladi R, Sethian J, Vemuri B:
**Shape modeling with front propagation: A Level Set Approach.***IEEE Trans Pattern Anal Machine Intell*1995,**17**(2):158–175. 10.1109/34.368173View ArticleGoogle Scholar - Osher S, Fedkiw R:
*Level set methods and dynamic implicit surfaces.*New York, Springer-Verlag; 2002.Google Scholar - Zhang KH, Zhang L, Song HH:
**Active contours with selective local or global segmentation: A new formulation and level set method.***Image and Vision Computing*2010,**28:**668–676. 10.1016/j.imavis.2009.10.009View ArticleGoogle Scholar - Li CM, Kao CY, Gore JC, Ding ZH:
**Minimization of region-scalable fitting energy for image segmentation.***IEEE Trans Image Process*2008,**17:**1940–1949.MathSciNetView ArticleGoogle Scholar - Liu JX, Chen YS, Chen LF:
**Accurate and robust extraction of brain regions using a deformable model based on radial basis functions.***Journal of Neuroscience Methods*2009,**183:**255–266. 10.1016/j.jneumeth.2009.05.011View ArticleGoogle Scholar - Lie J, Lysaker M, Tai XC:
**A binary level set model and some applications to Mumford-Shah image segmentation.***IEEE Trans Image Process*2006,**15**(5):1171–1181.MATHView ArticleGoogle Scholar - Huh S, Ketter TA, Sohn KH, Lee C:
**Automated cerebrum segmentation from three-dimensional sagittal brain MR images.***Comput Biol Med*2002,**32:**311–328. 10.1016/S0010-4825(02)00023-9View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.