# A Modified Brain MR Image Segmentation and Bias Field Estimation Model Based on Local and Global Information.

1. IntroductionAs an important medical treatment technique, magnetic resonance imaging (MRI) has greatly enhanced the efficiency of the doctor's diagnosis and avoided the numerous anatomical surgeries or diagnostic laparotomies. In various applications of MR image, brain MR image has played a vital role in the detection and diagnosis of brain diseases. In practice, however, brain MR image has some deficiencies such as noise interference, intensity inhomogeneity, low contrast, and the partial volume effect of brain tissue. Therefore, it is very difficult to segment accurately brain MR images. Now, most of the image segmentation methods assume that the intensity distribution of MR images is uniform, only from the perspective of antinoise performance to segment MR image, which is bound to cause inaccurate segmentation result. For this reason, to achieve higher segmentation accuracy, it is necessary to estimate the bias field of MR images.

Currently, fuzzy C-means (FCM) clustering technique has been widely applied in image segmentation, and it has the advantages of unsupervised segmentation, simple arithmetic, and fast speed of computation [1-3]. Meanwhile, the algorithm also has shortcomings. It only considers single pixel information in the image while ignoring the influence of neighborhood spatial pixel, which is sensitive to noise. Therefore, many scholars put forward a series of improved algorithms [4-9]. Ahmed et al. [4] presented an algorithm that incorporated the neighboring spatial information of a pixel in local window into the objective function of the FCM algorithm to improve antinoise capability (called BCFCM). Cai et al. [5] modified BCFCM algorithm (called FGFCM) where the gray level value and local spatial information of adjacent pixels are all similarly considered in the measurement. Krinidis et al. [6] proposed a clustering algorithm called FLICM (fuzzy local information C-means), which introduces an adaptive control factor in its objective function without trial-and-error experiments, as well as obtaining better antinoise performance and image segmentation accuracy. However, the aforementioned algorithms only considered the local neighborhood information of the image without taking into account the global structural information of an image. As a result, the FCM algorithm with local spatial constraints generally cannot get satisfactory segmentation results.

In recent years, a few FCM algorithms with global constrained methods were presented [10-12]. Zhao et al. [11] proposed an efficient fuzzy clustering scheme, in which a nonlocal constraint item is introduced into the objective function of the improved FCM, and the global structure information of the image plays a very important role in the process of image segmentation. Feng et al. [12] proposed a nonlocal FCM clustering method with edge preservation, and the detail information of the image can be well preserved. The antinoise ability of these algorithms has been greatly improved for image segmentation. However, these algorithms ignored the correction of intensity nonuniformity in MR image.

In order to get rid of the influence of bias field, many studies have been launched in this respect by some scholars [13-20]. Wells et al. [13] suggested an adaptive estimation method of bias field using the expectation maximization model. The model can automatically segment each brain tissue, but it needs to accurately know the specific information of each tissue and results in high computational complexity. Sled et al. [14] proposed a nonparametric nonuniform intensity normalization (N3) approach in MR data. The method is iterative and seeks the smooth multiplicative field that maximizes the high frequency content of the distribution of tissue intensity, and it requires no a priori knowledge and can be applied to almost any MR image. Tustison et al. [15] improved N3 algorithm by replacing the B-spline smoothing strategy used in the original N3 framework with an advantageous alternative and modified the iterative optimization scheme to improve convergence performance (called N4ITK). N3 and N4ITK are available to the public, so they have been applied widely in medical image processing. However, N3 and N4ITK are only designed for bias correction, and the parameters estimated for B-spline are very time-consuming; they are generally used for image segmentation [16]. Li et al. [17] presented a scheme of bias field estimation and image segmentation based on coherent local intensity clustering (CLIC); the Gaussian kernel function is incorporated into the weighting of the local neighborhood in the model to measure the spatial distance between neighbor pixels in FCM objective function. Although CLIC model can correct the bias field, this algorithm only utilizes grayscale distribution information in local neighborhood without considering global information of the whole image. Therefore, CLIC algorithm is sensitive to noise and has lower segmentation accuracy.

Through the analysis mentioned above, it can be seen that FCM and its variants only use the difference between the central pixel and its neighboring pixels to calculate the local similarity measure. When both the central pixel and some of the neighboring pixels are abnormal pixels, it may fail to analyze the impact of each neighboring pixel on the local similarity measure exhaustively. In this condition, the local spatial information derived from the image may play a negative role in guiding the noisy image segmentation. However, the nonlocal information can take advantage of the high degree of redundancy in the image. In other words, every pixel in the image can find a set of samples with a similar neighborhood configuration to it. Then, the pixel under consideration could be influenced by the weighted averaging over these samples. If nonlocal adaptive spatial constraint term is introduced, then the objective function, the local neighborhood information, and global structure information of an image can be comprehensively utilized. Therefore, a novel model of simultaneous MRI brain tissue segmentation and bias field estimation is proposed in this paper, which includes two improvements. Firstly, an improved FCM objective function is constructed including local spatial constraint term and nonlocal constraint term, which can significantly enhance the antinoise performance and the segmentation accuracy. Secondly, the bias field was coupled with the model as a multiplicative additional field, and thus bias field can be effectively estimated from MR images to reduce the impact of intensity inhomogeneity for image segmentation.

2. Related Work

2.1. FCM Clustering with Spatial Constraints. In view of some drawbacks of standard FCM algorithm, a modified scheme is proposed by Ahmed et al. [4]. The objective function introduces a regularization term with spatial neighborhood constraints in the local window, and the new objective function is

[mathematical expression not reproducible], (1)

where [x.sub.k] is the gray value of the kth pixel, [v.sub.i] is the ith cluster center, [u.sub.ik] is the fuzzy membership degree of the kth pixel belonging to the ith cluster and follows the constraint [[summation].sup.c.sub.i=1] [N.sub.[OMEGA]] is the number of neighborhood pixels, [x.sub.r] is the gray value of the neighboring pixels in local window, and [[OMEGA].sub.k] represents the set of neighbors that exist in a window around [x.sub.k]. The parameter m is a weighting exponent controlling classification results, and [alpha] is a key parameter to control the balance between the first term and the second term in (1). By minimizing (1), [u.sub.ik] and [v.sub.i] can be derived as follows:

[mathematical expression not reproducible]. (2)

2.2. FCM Clustering with Improved Fuzzy Partitions. FCM clustering with spatial constraints and its variants have better segmentation performance than the standard FCM algorithm. However, they still have some disadvantages. For example, their antinoise ability is not strong enough, and the execution efficiency of the algorithm is low, and the adaptive ability of related parameters is poor. Zhu et al. [10] proposed a more effective FCM clustering model (called GIFP_FCM), in which a novel energy function is built by introducing a new membership constraint item; the convergence speed and segmentation accuracy of GIFP_FCM algorithm get more significant enhancement than FCM. The corresponding objective function of GIFP_FCM is as follows:

[mathematical expression not reproducible], (3)

where [a.sub.k] is an adjusted parameter for controlling the clustering performance of GIFP_FCM, [a.sub.k] = [alpha] min {[[parallel] [x.sub.k] - [v.sub.i] [parallel].sup.2], i [member of] {1, 2, ..., c}}, and 0 < [alpha] < 1. By minimizing (3) using Lagrangian optimization, the fuzzy membership [u.sub.ik] and clustering center [v.sub.i] are derived as follows:

[mathematical expression not reproducible]. (4)

2.3. Bias Field Model. In practical application, affected by radio frequency coil, magnetic resonance equipment, and the variability of different brain tissue, intensity inhomogeneity generally appears in MR image. As a result, bias field correction is a very important task for MR image segmentation.

In MR image, bias field presents smooth variation of gray level of the pixels of the same tissue in spatial domain. Therefore, bias field can be regarded as a multiplicative component of an MR image. Let X be the observed image and [X.sub.0] be the true image; b and N denote bias field and noise in the image, respectively. The model is as follows:

X = b[X.sub.0] + N. (5)

2.4. CLIC Model. In CLIC model [17], a Gaussian kernel function is first used in local neighborhood by standard FCM algorithm to segment the MR image and then extended to the whole image. In the small neighborhood of an image, bias field can be approximated as constant; the pixel gray level can also be regarded as a constant. In this case, the bias field can be estimated when ignoring the additional noise in (5). The clustering criterion function of CLIC is

[mathematical expression not reproducible], (6)

where b(r) is the bias field, K(r - k) is the weight of a truncated Gaussian kernel assigned to the intensity [x.sub.k], and the weighting function K is written as follows:

[mathematical expression not reproducible], (7)

where [sigma] is the standard deviation of the Gaussian kernel, [rho] is the radius of the adjacent window, and a is a constant to normalize the Gaussian kernel.

In CLIC model, image segmentation and bias field estimation can be completed simultaneously, but this model also has some drawbacks. First of all, as for the weighting of spatial measure, Gaussian kernel is only related to the local neighborhood of current target pixel, without taking into account the texture structure of the whole image. As a result, some of the pixels will get the error classification. Second, the model cannot effectively eliminate the influence of noise, because CLIC is modeled based on the standard FCM. Finally, a large number of nuclear convolution calculations result in high computational complexity in CLIC.

3. The Proposed Method

During MR brain image segmentation, the local neighborhood information and global structure information of the image should be used for ensuring the accuracy of the segmented results and robustness to noise. In this section, we reconstructed the objective function based on CLIC and GIFP_FCM algorithms, and a nonlocal regularization term is added in the modified model. In fact, the variation pattern of MR image pixels has certain regularity; that is, the brain tissue region has a strong global similarity, and it should be treated with nonlocal constraint.

3.1. Nonlocal Spatial Constraint. Given a noisy MR image, the similarity between the kth and lth pixels is measured by their gray value vectors x([N.sub.k]) and x([N.sub.l]), where x([N.sub.k]) and x([N.sub.l]) represent the square neighborhood of radius s and centered at pixels [x.sub.k] and [x.sub.l], respectively. The similarity of gray value vector is determined by Gaussian weighted Euclidean distance [[parallel] x([N.sub.k]) - x([N.sub.l]) [parallel].sup.2.sub.2,[sigma]]; the expression is as follows:

[mathematical expression not reproducible], (8)

where [sigma[ denotes the standard deviation of the Gaussian kernel, x([N.sub.k]) denotes a window of radius s around [x.sub.k], [x.sup.(p)]([N.sub.k]) is the pth element in this window, and [[sigma].sup.(p)] is defined as follows:

[[sigma].sup.(p)] = [s.summation over (v=max(d,1))] 1 / [(2V +1).sup.2] s, (9)

where [y.sub.p] = mod(p, (2s + 1)) and [z.sub.p] = floor(p, (2s + 1)) + 1. ([y.sub.p], [z.sub.p]) denote the coordinates of the pth element in the square window x([N.sub.i]), d = max([absolute value of ([y.sub.p] - s - 1)], [absolute value of ([z.sub.p] - s - 1)]). Therefore, the similarity measure [S.sub.kl] between the kth and lth voxels in the MR image is as follows:

[S.sub.kl] = 1 / [Z.sub.k] e (- [[parallel] x ([N.sub.k]) - x ([N.sub.l]) [parallel].sup.2.sub.2,[omega]]/h), (10)

where h is a constant proportional to the noise deviation [sigma]; it controls the decay of the similarity measure [S.sub.kl] and satisfies 0 < [S.sub.kl] < 1 and [mathematical expression not reproducible] represents a search window of radius u around the fcth voxel. [Z.sub.k] denotes a normalized constant; it is defined as follows:

[mathematical expression not reproducible]. (11)

Therefore, for the kth pixel of the image, the nonlocal average method obtains its estimate value [[bar.x].sub.k] utilizing a weighted average method, and the expression is as follows:

[[bar.x].sub.k] = [summation over (l[member of][w.sup.u.sub.k])] [S.sub.kl] [x.sub.l]. (12)

3.2. Objective Function. The modified objective function incorporating local operation and nonlocal operation is as follows:

[mathematical expression not reproducible], (13)

where [u.sub.ik] is the membership degree of fuzzy clustering, [v.sub.i] is the ith cluster center, [[OMEGA].sub.k] represents a neighborhood window of radius r and centered at [x.sub.k], and [[bar.x].sub.k] is the nonlocal calculated value of the kth pixel, and the definition of K(r-k) and b(r) is the same as (6). [[beta].sub.k] is the weighting parameter to balance local coupling item and nonlocal regularization item of (13); the definition of parameter [[beta].sub.k] is

[[beta].sub.k] = max ([w.sup.u.sub.k]) - min ([w.sup.u.sub.k]) / max ([w.sup.u.sub.k]). (14)

Theorem 1. Assume [[summation].sup.c.sub.i=1] [u.sub.ik] = 1, 0 [less than or equal to] [u.sub.ik] [less than or equal to] 1, and m > 1.By an optimization using Lagrange multiplier method in (13), [u.sub.ik] and [v.sub.i] can be derived as follows:

[mathematical expression not reproducible], (15)

[mathematical expression not reproducible]. (16)

Proof. According to Lagrange multiplier method, (13) can be converted to unconstrained optimization problem:

[mathematical expression not reproducible]. (17)

Compute the partial derivative of L with respect to [u.sub.ij] and [[lambda].sub.k], respectively, and let [partial derivative]L / [partial derivative][u.sub.ik] = 0 and [partial derivative]L / [partial derivative][[lambda].sub.k] = 0; that is,

[mathematical expression not reproducible], (18)

[partial derivative]L / [partial derivative][[lambda].sub.k] = 1 [c.summation over (i=1)][u.sub.ik] = 0. (19)

From (18), we obtain

[mathematical expression not reproducible]. (20)

By substituting (20) into (18), we obtain

[mathematical expression not reproducible]. (21)

By substituting (21) into (20), we obtain

[mathematical expression not reproducible]. (22)

Similarly, let

[partial derivative]L / [partial derivative][v.sub.k] = 0; that is,

[mathematical expression not reproducible]. (23)

From (23), we obtain

[mathematical expression not reproducible]. (24)

The theorem is completely proved.

3.3. Bias Field Estimation. To estimate the bias field [b.sub.k] of the image in (17), we adopt the same analysis method and take the partial derivative of L with respect to [b.sub.k]. Let [partial derivative]L / [partial derivative][b.sub.k] = 0; is,

[mathematical expression not reproducible]. (25)

From (25), [b.sub.k] can be obtained:

[mathematical expression not reproducible]. (26)

3.4. The Algorithm Flow. The program flow of the proposed algorithm can be summarized as follows.

Step 1. The following are given: the number of clusters c, the exponent of fuzziness m, the radius of local window s, the radius of search window U, bias field [b.sub.0] = 1,andstopcriterion

Step 2. Get the fuzzy cluster prototypes [V.sup.(1)] = {[v.sup.(1).sub.1]) , [v.sup.(1).sub.2], ..., [v.sup.(1).sub.c]} using C-means algorithm and set the iteration initial value t = 0.

Step 3. Compute nonlocal similarity measure [S.sub.kl] using (10) and then obtain the nonlocal weighted average [[bar.x].sub.k] for the kth pixel using (12).

Step 4. Compute and update the membership degree [u.sup.(t).sub.ik] by (15).

Step 5. Compute and update the clustering prototypes [v.sup.(t).sub.i] by (15).

Step 6. Get the bias field estimate [b.sup.(t).sub.k] by (26).

Step 7. If max [parallel] [V.sup.(t+1)] - [V.sup.(t)] [parallel] < [epsilon], then output results; otherwise, go to Step 3 and set t = t + 1.

4. Experimental Results

To validate the validity of the proposed algorithm, several bias field estimation and image segmentation algorithms are taken as comparative methods. In all experiments, the parameter settings are as follows: the stop criterion [epsilon] = 0.001 and m = 2 and the radius of two windows r = 10 and s = 3, respectively.

4.1. Bias Field Correction. First of all, the proposed method is applied to the 1.5T- and 3T-weighted brain MR images, and Figure 1 illustrates the experimental results. Column (a) shows three MR images with different level intensity nonuniformity (INU), column (b) shows the estimated bias field, column (c) shows the corrected images, and column (d) shows the segmentation results. In column (d), the brain tissue of corrected image becomes very homogeneous, and the cerebrospinal fluid (CSF), white matter (WM), and grey matter (GM) can be also clearly identified.

In addition, Figure 2 compares results by BCFCM [4], N3 [14], N4ITK [15], and the proposed algorithm on MR images. The original images, estimated bias fields, and bias corrected images are shown in the first, second, and third columns of Figures 2(a)-2(d), respectively. The histogram of the original MR image and the histograms of the bias corrected images by BCFCM, N3, N4ITK, and our methods are shown in Figure 3. In the histograms, the right two significant peaks correspond to the GM and WM, respectively. The peak of the cerebrospinal fluid (CSF) is not distinct since its volume is relatively small. From the histograms of the bias corrected images recovered by BCFCM, N3, N4ITK, and our methods, we see that the histograms of specific tissues approximately satisfy Gaussian distribution but with significantly different variances. These results validate that our model is more consistent with the intensity distribution of the image with intensity inhomogeneity than other algorithms.

4.2. Antinoise Ability Analysis. In the second experiment, the images from BrainWeb simulated MR image database [21] are applied to analyze the antinoise ability of five algorithms (FCM, BCFCM [4], GIFP_FCM [10], CLIC [17], and our method). These simulated brain images are of T1-weighted contrast with a 1mm slice thickness, 15% Rician noise, and no intensity inhomogeneities. The noisy MR images are segmented into four classifications: CSF, WM, GM, and the background. Before segmentation, the extracranial tissues were removed. The segmentation results using the five methods are shown in Figure 4, respectively.

Figure 4(a) is the MR image with 15% (l = 15) Rician noise. Figures 4(b)-4(f) show the corresponding experimental results by the five methods, respectively. Figure 4(g) shows the ground truth. Figure 4 shows that the ability of our method is excellent at the preservation of image detail and robustness to noise, while the other four compared algorithms are relatively poor.

In order to quantitatively analyze the antinoise performance of the five algorithms, 8 MR images (the level of Rician noise ranges from 5% to 20%) are selected as the experimental samples. The statistical results (average values) of the Jaccard similarity (JS) [22] values of GM, WM, and CSF are shown in Table 1. JS was used for comparison and quantitative evaluation. Hence,

JS = [absolute value of ([A.sub.i] [intersection] [B.sub.i])] / [absolute value of ([A.sub.i] [union] [B.sub.i])], (27)

where [A.sub.i] denotes the set of pixels belonging to the ith class identified by the clustering algorithm, while [B.sub.i] denotes the set of pixels belonging to the ith class in the ground truth. As a fuzzy similarity measure, the larger the JS value, the better the clustering performance. It can be seen from the experimental results in Table 1 that, with the increase of noise level of MR image, JS values of all algorithms are reduced. However, our method has higher values than the other four algorithms, and it illustrates that the proposed algorithm has better clustering performance and stronger robustness to noise.

4.3. Brain Tissue Segmentation of Noisy MR Image. In the third experiment, three MR images are selected for brain tissue segmentation by five algorithms, and these images are corrupted by 12% Rician noise and 40% intensity inhomogeneity. The noisy MR images are shown in Figure 5(a), and the corresponding segmentation results of FCM, BCFCM, GIFP_FCM, CLIC, and the proposed algorithm are shown in Figures 5(a)-5(f), respectively. Figure 5(g) shows the ground truth of the three test images. Figure 5 shows that the proposed algorithm not only can estimate bias field but also can effectively ensure segmentation accuracy of three MR images.

To evaluate the segmentation accuracy with the variety of bias fields, we tested 10 MR brain images using three bias field estimation algorithms: BCFCM, CLIC, and the proposed method. The JS values of WM, GM, and CSF are compared when the INU level of bias fields varied from 10% to 80% for the MR images with noise 9%, and the comparison results are illustrated in Figure 6, respectively. It can be observed from Figure 6 by quantitative comparison that JS values of three algorithms are progressively smaller with the increase of the INU level, and our proposed model can preserve more image detail than the other two methods.

5. Conclusion

The traditional FCM algorithms with local information to suppress noise have some limitations in image segmentation. Due to the nonlocal weighted measure being unable to reflect the distance of real pixel points to the clustering center, a robust MRI brain tissue segmentation and bias field estimation model based on local and global information is proposed in this paper. During the process of calculating similarity measure, the algorithm utilized the nonlocal information to adjust parameters and reduce the difficulty of parameter setting; meanwhile, bias fields of MR brain image with intensity inhomogeneity is corrected. The segmentation results of the experiment showed that the proposed method presents stronger robustness to noise and accurately estimates the bias field, and more detailed structure of the MR brain image can be efficiently preserved.

http://dx.doi.org/10.1155/2016/9871529

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant no. 81301297), the National Key Scientific Instrument and Equipment Development Projects of China (2012YQ04014001 and 2012YQ04014010), Natural Science Foundation of Heilongjiang Province (F201241), the Fundamental Research Funds for the Central Universities (HEUCF160412), and Developing Key Laboratory of Sensing Technology and Systems in Cold Region of Heilongjiang Province and Ministry of Education, Heilongjiang University, China (Grant no. P201410).

References

[1] S. Naz, H. Majeed, and H. Irshad, "Image segmentation using fuzzy clustering: a survey," in Proceedings of the 6th International Conference on Emerging Technologies (ICET '10), pp. 181-186, IEEE, Islamabad, Pakistan, October 2010.

[2] K. S. Tan and N. A. M. Isa, "Color image segmentation using histogram thresholding Fuzzy C-means hybrid approach," Pattern Recognition, vol. 44, no. 1, pp. 1-15, 2011.

[3] M. Gong, Z. Zhou, and J. Ma, "Change detection in synthetic aperture radar images based on image fusion and fuzzy clustering," IEEE Transactions on Image Processing, vol. 21, no. 4, pp. 2141-2151, 2012.

[4] M. N. Ahmed, S. M. Yamany, N. Mohamed, A. A. Farag, and T. Moriarty, "A modified fuzzy C-means algorithm for bias field estimation and segmentation of MRI data," IEEE Transactions on Medical Imaging, vol. 21, no. 3, pp. 193-199, 2002.

[5] W. Cai, S. Chen, and D. Zhang, "Fast and robust fuzzy c-means clustering algorithms incorporating local information for image segmentation," Pattern Recognition, vol. 40, no. 3, pp. 825-838, 2007.

[6] S. Krinidis and V. Chatzis, "A robust fuzzy local information C-means clustering algorithm," IEEE Transactions on Image Processing, vol. 19, no. 5, pp. 1328-1337, 2010.

[7] I. Despotovic, E. Vansteenkiste, and W. Philips, "Spatially coherent fuzzy clustering for accurate and noise-robust image segmentation," IEEE Signal Processing Letters, vol. 20, no. 4, pp. 295-298, 2013.

[8] M. Gong, Y. Liang, J. Shi, W. Ma, and J. Ma, "Fuzzy C-means clustering with local information and kernel metric for image segmentation," IEEE Transactions on Image Processing, vol. 22, no. 2, pp. 573-584, 2013.

[9] Z. Zhao, L. Cheng, and G. Cheng, "Neighbourhood weighted fuzzy c-means clustering algorithm for image segmentation," IET Image Processing, vol. 8, no. 3, pp. 150-161, 2014.

[10] L. Zhu, F.-L. Chung, and S. Wang, "Generalized fuzzy C-means clustering algorithm with improved fuzzy partitions," IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 39, no. 3, pp. 578-591, 2009.

[11] F. Zhao, L. Jiao, H. Liu, and X. Gao, "A novel fuzzy clustering algorithm with non local adaptive spatial constraint for image segmentation," Signal Processing, vol. 91, no. 4, pp. 988-999, 2011.

[12] J. Feng, L. C. Jiao, X. Zhang, M. Gong, and T. Sun, "Robust nonlocal fuzzy c-means algorithm with edge preservation for SAR image segmentation," Signal Processing, vol. 93, no. 2, pp. 487-499, 2013.

[13] W. M. Wells III, W. E. L. Crimson, R. Kikinis, and F. A. Jolesz, "Adaptive segmentation of MRI data," IEEE Transactions on Medical Imaging, vol. 15, no. 4, pp. 429-442,1996.

[14] J. G. Sled, A. P. Zijdenbos, and A. C. Evans, "A nonparametric method for automatic correction of intensity nonuniformity in mri data," IEEE Transactions on Medical Imaging, vol. 17, no. 1, pp. 87-97, 1998.

[15] N. J. Tustison, B. B. Avants, P. A. Cook et al., "N4ITK: improved N3 bias correction," IEEE Transactions on Medical Imaging, vol. 29, no. 6, pp. 1310-1320, 2010.

[16] K. Zhang, Q. Liu, H. Song, and X. Li, "A variational approach to simultaneous image segmentation and bias correction," IEEE Transactions on Cybernetics, vol. 45, no. 8, pp. 1426-1437, 2015.

[17] C. Li, C. Xu, A. Anderson et al., "MRI tissue classification and bias field estimation based on coherent local intensity clustering: a unified energy minimization framework," in Proceedings of the 21st International Conference on Information Processing in Medical Imaging (IPMI '09), pp. 289-299, Williamsburg, Va, USA, July 2009.

[18] S. Wang, X. Ji, J. Zhang et al., "Improved coupled model for MR images segmentation and bias restoration," Journal of Image and Graphics, vol. 17, no. 9, pp. 1175-1180, 2012.

[19] J. Zhang, H. Yang, Y. Chen, L. Fang, and T. Zhan, "Brain MR image segmentation and bias correction model based on non-local information," Journal of Computer-Aided Design & Computer Graphics, vol. 25, no. 4, pp. 526-540, 2013.

[20] C. Li, J. C. Gore, and C. Davatzikos, "Multiplicative intrinsic component optimization (MICO) for MRI bias field estimation and tissue segmentation," Magnetic Resonance Imaging, vol. 32, no. 7, pp. 913-923, 2014.

[21] C. A. Cocosco, V. Kollokian, R. K.-S. Kwan et al., "BrainWeb: Online Interface to a 3D MRI Simulated Brain Database," 2011, http://brainweb.bic.mni.mcgill.ca/brainweb/.

[22] F. Masulli and A. Schenone, "A fuzzy clustering based segmentation system as support to diagnosis in medical imaging," Artificial Intelligence in Medicine, vol. 16, no. 2, pp. 129-147, 1999.

Wang Cong, (1) Jianhua Song, (1,2) Kuan Luan, (1) Hong Liang, (1) Lei Wang, (1) Xingcheng Ma, (1) and Jin Li (1)

(1) College of Automation, Harbin Engineering University, Harbin 150001, China

(2) Electronic Engineering College, Heilongjiang University, Harbin 150080, China

Correspondence should be addressed to Jin Li; lijin@hrbeu.edu.cn

Received 1 February 2016; Revised 29 June 2016; Accepted 27 July 2016

Academic Editor: Rong Chen

Caption: Figure 1: Bias field estimation and corresponding segmentation results. (a) Three MR images, (b) the estimated bias field, (c) the images of removed bias field, and (d) the segmented images.

Caption: Figure 2: Comparison of the estimated bias fields of several algorithms. (a) BCFCM, (b) N3, (c) N4ITK, and (d) the proposed method.

Caption: Figure 3: Histograms of original image and bias corrected images in Figure 2. (a) original image, (b) BCFCM, (c) N3, (d) N4ITK, and (e) our method.

Caption: Figure 4: Antinoise ability comparison of the five algorithms. (a) Noisy image, (b) FCM, (c) BCFCM, (d) GIFP.FCM, (e) CLIC, (f) the proposed algorithm, and (g) ground truth.

Caption: Figure 5: Experimental results for threeMR images. (a) NoisyMR images, (b) FCM, (c) BCFCM, (d) GIFP FCM, (e) CLIC, (f) the proposed method, and (g) ground truth.

Caption: Figure 6: JS comparison of the brain tissue segmentation with three algorithms: (a) JS of WM segmentation, (b) JS of GM segmentation, and (c) JS of CSF segmentation.

Table 1: Comparison of JS on simulated MR images with different levels of noise. Algorithm Tissues 5% 10% 15% 20% FCM WM 0.8373 0.7458 0.6403 0.5773 GM 0.7189 0.5325 0.4211 0.3021 CSF 0.8359 0.6751 0.5238 0.4001 BCFCM WM 0.9022 0.8563 0.7459 0.6832 GM 0.8754 0.7681 0.5765 0.4308 CSF 0.8868 0.8354 0.6869 0.5369 GIFP.FCM WM 0.9491 0.8993 0.8357 0.7365 GM 0.9011 0.8377 0.7324 0.5793 CSF 0.9268 0.8634 0.7957 0.7263 CLIC WM 0.9007 0.8557 0.7461 0.6989 GM 0.8766 0.7706 0.6013 0.5426 CSF 0.8993 0.8457 0.6956 0.5536 The proposed WM 0.9513 0.9320 0.8697 0.7833 algorithm GM 0.9295 0.8723 0.7838 0.6619 CSF 0.9457 0.8803 0.8329 0.7937

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Cong, Wang; Song, Jianhua; Luan, Kuan; Liang, Hong; Wang, Lei; Ma, Xingcheng; Li, Jin |

Publication: | Computational and Mathematical Methods in Medicine |

Date: | Jan 1, 2016 |

Words: | 5030 |

Previous Article: | Fluid-Structure Simulations of a Ruptured Intracranial Aneurysm: Constant versus Patient-Specific Wall Thickness. |

Next Article: | Posterior Estimates of Dynamic Constants in HIV Transmission Modeling. |

Topics: |