Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Med Imaging**|**v.12; 2012**|**PMC3368740

Formats

Article sections

- Abstract
- Background
- Methods
- Results
- Discussion
- Conclusion
- Competing interests
- Authors' contributions
- Pre-publication history
- References

Authors

Related links

BMC Med Imaging. 2012; 12: 10.

Published online 2012 May 16. doi: 10.1186/1471-2342-12-10

PMCID: PMC3368740

Mohammad-Reza Nazem-Zadeh,^{1,}^{2,}^{3} Sona Saksena,^{4} Abbas Babajani-Fermi,^{4,}^{5} Quan Jiang,^{3} Hamid Soltanian-Zadeh,^{}^{1,}^{4,}^{6} Mark Rosenblum,^{5} Tom Mikkelsen,^{7} and Rajan Jain^{4,}^{7}

Mohammad-Reza Nazem-Zadeh: ude.hcimu.dem@azmezanm; Sona Saksena: ude.hfh.dar@sanos; Abbas Babajani-Fermi: ude.hfh.dar@bsabba; Quan Jiang: ude.hfh.dar@1gnaijQ; Hamid Soltanian-Zadeh: ude.hfh.dar@sdimah; Mark Rosenblum: gro.shfh@1bnesorm; Tom Mikkelsen: gro.shfh@1lekkimt; Rajan Jain: ude.hfh.dar@jnajar

Received 2011 August 8; Accepted 2012 May 16.

Copyright ©2012 Nazem-Zadeh et al; licensee 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.

This article has been cited by other articles in PMC.

This paper presents a three-dimensional (3D) method for segmenting corpus callosum in normal subjects and brain cancer patients with glioblastoma.

Nineteen patients with histologically confirmed treatment naïve glioblastoma and eleven normal control subjects underwent DTI on a 3T scanner. Based on the information inherent in diffusion tensors, a similarity measure was proposed and used in the proposed algorithm. In this algorithm, diffusion pattern of corpus callosum was used as prior information. Subsequently, corpus callosum was automatically divided into Witelson subdivisions. We simulated the potential rotation of corpus callosum under tumor pressure and studied the reproducibility of the proposed segmentation method in such cases.

Dice coefficients, estimated to compare automatic and manual segmentation results for Witelson subdivisions, ranged from 94% to 98% for control subjects and from 81% to 95% for tumor patients, illustrating closeness of automatic and manual segmentations. Studying the effect of corpus callosum rotation by different Euler angles showed that although segmentation results were more sensitive to azimuth and elevation than skew, rotations caused by brain tumors do not have major effects on the segmentation results.

The proposed method and similarity measure segment corpus callosum by propagating a hyper-surface inside the structure (resulting in high sensitivity), without penetrating into neighboring fiber bundles (resulting in high specificity).

Corpus callosum is the largest inter-hemispheric fiber bundle in the human brain [1,2]. Most of the fibers interconnect homologue cortical areas in roughly mirror image sites but a large number of the fibers have heterotypic connections ending in asymmetrical areas [3]. Previous studies have mainly investigated effects of various pathologies on the corpus callosum [4-7]. However, a fully automated, fast, and accurate method for segmenting corpus callosum without penetrating into irrelevant neighboring structures, using data acquired in routine clinical protocols, is still lacking.

Previously, image processing methods have been proposed for segmenting corpus callosum in anatomical magnetic resonance images (MRI) [8-10]. These methods rely on intensity information of two-dimensional images and their results may need pruning. Recently, attention has been oriented towards diffusion tensor imaging (DTI) to segment white matter tracts of the brain [11,12]. Although the tensor model fails to describe higher order anisotropies in heterogeneous areas where more than one fiber population exists, it is practically useful for extracting major white matter tracts, particularly the ones with predominant diffusivity pattern such as corpus callosum. When using DTI data, the fiber bundles can be extracted by: a) clustering of fibers resulting from tractography into fiber bundles [13-19]; or b) segmenting fiber bundles via hyper-surface propagation based on local properties of diffusion tensor, diffusion signal, or orientation distribution function (ODF) [20-27]. Since clustering methods rely on the tractography results, they do not work properly if the tractography results are inaccurate. On the other hand, segmentation methods based on hyper-surface propagation do not use the tractography results and are thus more robust to noise.

In a region-based segmentation framework, a similarity measure between successive tensors is typically used. Some of the hyper-surface propagating methods in the literature concentrated on scalar quantities derived from the tensor data which do not reflect complete tensor information [20]. Other methods benefit from the entire information contained in the DTI data [21-30]. Wang and Vemuri [22] proposed a statistical level-set segmentation method. However, the tensors derived in this framework are not necessarily positive semi-definite, leading to inappropriate results especially when consecutive tensors are much different. Metrics like Kullback-Leibler divergence and J-divergence [23,24] have also been proposed. One of the most promising methods is introduced by Jonasson et al. [25]. They defined a new similarity measure called normalized tensor scalar product (NTSP). Comparing NTSP with other similarity measures, they demonstrated superiority of their proposed measure. To segment brain structures like thalamic nuclei, they modified their framework to favor the propagation of multiple hyper surfaces without overlapping [26]. Lenglet et al. [27] defined a dissimilarity measure and statistics between tensors based on the Riemannian distances. Although improving the segmentation results, this approach is computationally expensive. Defining a Log-Euclidean distance, another metric was defined by Arsigny et al. [28] which has lower computational burden. Weldeselassie and Hamarneh [29] used their proposed similarity measure in an energy minimization framework. Awate et al. [30] used the similarity measure in a Markov random field framework.

In terms of quantitative evaluation of diffusion parameters, previous studies have compared DTI-based indices in normal appearing white matter and corpus callosum in multiple sclerosis [4], stroke [5], schizophrenia [6], and Huntington's [7] and also studied the DTI methods to assess corpus callosum regions across the human lifespan [31]. For segmenting corpus callosum and its subdivisions in these studies, however, two-dimensional (2D) methods were applied and DTI-based indices compared in the mid-sagittal plane. However, without recruiting a three dimensional (3D) method to segment the whole corpus callosum and its subdivisions, the extracted quantities may be inaccurate.

Since the tensor model is not capable of describing heterogeneous diffusion behavior in the crossing fiber bundles, some studies used High Angular Resolution Diffusion Imaging (HARDI) data to segment specific bundles [32-38]. However, the HARDI data is not widely acquired in clinical centers and hence, the tensor-based methods are still of more practical use in clinical research.

In this manuscript, we present a 3D method to segment corpus callosum. Based on tensor and anisotropy values of neighboring voxels, a similarity measure is proposed and used as a speed function in the proposed level-set method. In this method, the principal diffusion direction (PDD) and prior information about the diffusivity pattern in corpus callosum are used to avoid inclusion of neighboring fiber bundles. Then, the Witelson subdivisions of corpus callosum are automatically identified [39]. The idea of using diffusivity pattern in corpus callosum has been used by Lee et al. [40]. However, they performed a 2D segmentation on the mid-sagittal plane. Moreover, since their method uses the left-right component of PDD to delineate corpus callosum boundaries, it is not applicable in more lateral sagittal planes, where corpus callosum connects to minor and major forceps and considerable anterior-posterior component of PDD exists.

The proposed segmentation method can be identically used for segmenting corpus callosum of control subjects as well as patients with glioblastoma. However, since we use the geometric information of the diffusivity pattern in corpus callosum and the glioblastoma tumor may change the original shape and diffusivity pattern of corpus callosum, we validate the reproducibility of the proposed method through realistic simulations of corpus callosum rotations under glioblastoma tumor pressure. We study the effects of rotation by different Euler angles (azimuth, elevation and skew) quantitatively and demonstrate that even in extreme cases with large rotations, brain tumors do not have major effects on the segmentation results generated by our proposed method. We apply the method to the DTI data of normal subjects and brain cancer patients with glioblastoma and show its superiority to some of the previously published methods in the literature.

This study is approved by the institutional review board and is compliant with the Health Insurance Portability and Accountability Act (HIPAA). Between February 2006 and December 2008, 19 patients (8 males, 11 females; mean age 60.7 years) with treatment naïve glioblastoma underwent MRI with DTI using a 3 T scanner at our institution. Based on anatomical MRI findings, the patients were divided into two groups: Group 1 including patients with tumors not infiltrating corpus callosum (n = 12); and Group 2 including patients with tumors infiltrating corpus callosum (n = 7). For the control group, 11 patients (5 males, 6 females; mean age 48 years) who underwent brain MRI for nonspecific headache or single idiopathic seizure with normal MRI were included.

All the patients underwent both conventional MRI and DTI on a 3 T scanner (Excite HD, GE Medical Systems, Milwaukee, WI) using an 8-channel head coil. Diffusion weighted images (DWIs) were acquired in 25 diffusion gradient directions. The reconstructed DWIs have intra-slice resolution of 256 × 256 with voxel size of 0.98 × 0.98 × 2.5 mm. To have the same step size in each direction for the front propagation in the level-set method and to avoid extensive computation, we interpolated the data into 128 × 128 grid with cubical, 1.9 × 1.9 × 1.9 mm voxels.

We use a level-set method that takes into account several properties when formulating the segmentation problem [41,42]. The method smoothes the propagating hyper-surface automatically and leads to a regularized segmentation. Among the advantages of the method, its property of generalizing from 2D to 3D and higher dimensions and automatic splitting and merging of the surfaces are notable. In addition, it simplifies calculation of geometric quantities needed in the proposed method such as normal to surface and curvature.

For bundle segmentation in areas with high level of similarity in diffusion, the hyper-surface needs to grow in the direction normal to the surface. Moreover, the resulting fiber bundles must be smooth. Consequently, we use a level-set whose growing hyper-surface is the zero level-set of the following function which is a reduced form of the Hamilton-Jacobi partial differential equation:

$${D}_{t}\phi \left(r,t\right)+F\left(r,t\right)\left|\right|\nabla \phi \left(r,t\right)\left|\right|-\kappa \left(r,t\right)\left|\right|\nabla \phi \left(r,t\right)\left|\right|=0$$

(1)

where r ^{n }is the state space, : ^{n }× → is the level-set function, D_{t } is the partial derivative of with respect to the time variable *t*, = D_{r } is the gradient of with respect to the state space variables, and F(*r, t*) is the speed in the direction normal to the surface, extracted from the spherical harmonic coefficients of the neighboring voxels. The sign · stands for the magnitude operator. The curvature κ(*r, t*) is used to fulfill the smoothness constraint. To calculate the first order spatial partial derivative (*r, t*), we use the 5th order of the upwind method [41].

Corpus callosum is a commissural fiber bundle with a specific diffusion pattern which can be coded as prior knowledge in the segmentation framework. Using this information, we prevent the hyper-surface from propagating into adjacent white matter structures such as cingulum, tapetum, minor and major forceps, and tracts of the corona radiata. Although dissimilarity among tensors helps in this case, smooth and gradual transition in shape and direction of the DTI tensors from corpus callosum to minor and major forceps makes the segmentation difficult.

We define a similarity measure between every voxel on the propagating hyper-surface and its neighbors in the propagation direction, based on tensor and anisotropy values. This similarity measure is used as the speed term F(*r, t*) in Equation (1). The hyper-surface propagation in each step depends only on the speed function of the boundary voxels. This speed function term moves the hyper-surface to fill the whole fiber bundle, while the regularizing curvature term in Equation (1) is in charge of smoothing corpus callosum without changing its real structure.

Using the fact that the diffusivity in corpus callosum is perpendicular to the mid-sagittal plane of the brain, we consider a threshold (*PDD _{x }Threshold*) for the x-component (left-right component) of the PDD (

In the proposed algorithm, we consider two more thresholds on collinearity of the PDD vectors (*Collinearity_Threshold*) and similarity of fractional anisotropy (FA) values (*FA_Threshold*) in the neighboring voxels.

The proposed segmentation steps are as follows:

1. Select the initial seeds in corpus callosum in mid-sagittal plane manually.

2. Initiate the hyper-surface as the congregation of small spheres around the seed points.

3. Do until convergence

• For each point *r *on the hyper-surface at step *t*:

a) Calculate the normal direction to the surface.

b) Calculate the 26-neighborhood and keep the neighbors *n _{r }*for

c) If *PDD*(*r*). *PDD*(*n _{r}*) >

$$\text{Then}\phantom{\rule{1em}{0ex}}F\left(r,t\right)=\sum _{{n}_{r}}FA\left(r\right).FA\left({n}_{r}\right).\frac{tr\left[D{\left(r\right)}^{*}D\left({n}_{r}\right)\right]}{tr{\left[D\left(r\right)\right]}^{*}tr\left[D\left({n}_{r}\right)\right]}$$

(2)

• Threshold *F*(*r, t*) with *F _ Threshold *to diminish the effect of negligible speeds.

• Use the resultant speed in the level-set framework.

4. Extract the zero level-set as the segmented corpus callosum.

In Equation (2), *tr(.) *is the matrix trace and *D(r) *is the tensor at point *r*.

After segmenting corpus callosum, Witelson subdivisions of corpus callosum are automatically extracted [39]. First, the critical point between genu and rostrum of corpus callosum is calculated, where the curvature of the structure boundary in mid-sagittal plane changes. Then, the segmented corpus callosum is automatically subdivided into Witelson subdivisions in the mid-sagittal plane: rostrum, genu, rostral body, anterior mid-body, posterior mid-body, isthmus and splenium [39]. Moreover, the user can visualize and confirm the calculated mid-sagittal plane and the critical joining point between genu and rostrum. The critical point can be selected manually if close supervision is preferred or needed.

We implemented the proposed algorithm in MATLB R2008a using a PC with Intel^{® }Core™ 2Duo CPU (E8400@ 3.00 GHz, 3.00 GHz) and 4 GByte RAM and 64 bit VISTA operating system.

To evaluate quality of the segmentation results, the dice correctness measure [43] is calculated using the following relation:

$$Correctness=\frac{N\left({S}_{a}\cap {S}_{y}\right)}{\left[N\left({S}_{a}\right)+N\left({S}_{y}\right)\right]/2}$$

(3)

where S_{a }and *S _{r }*are the automatic and manual (reference) segmentation results, respectively, and

We tested sensitivity of corpus callosum Witelson segments to *PDD _{x}_Threshold *for a normal subject. Figures Figures11 and and22 show the number of True-Positives, the number of False-Positives, and the Dice correctness measure over a range of

Figure Figure33 shows the number of True-Positives, the number of False-Positives, and the Dice correctness measure over a range of *Collinearity _Threshold *for the Witelson segments of the Corpus Callosum of a normal subject. It is shown that when *Collinearity _Threshold *values are in a specific range (here 0.65 to 0.75), the Dice correctness measure is maximized. However, the sensitivity of the segmentation method to *Collinearity _Threshold *is less than its sensitivity to *PDDx_Threshold*. We found the following set of parameters optimal for segmenting corpus callosum: PDD_{x}_Threshold = 0.55, Collinearity _Threshold = 0.7, *FA_Threshold *= 0.1, and *F_Threshold *= 0.05.

The tensor-based method proposed by Jonasson et al. [32] with NTPP similarity measure is quite sensitive to the speed threshold. If the speed threshold is selected high enough to prevent the front from propagating into the neighboring structures, corpus callosum is not segmented entirely (low sensitivity). However, if the speed threshold is chosen low enough to segment the entire corpus callosum, the front propagates into irrelevant fiber structures such as superior longitudinal fasciculus, cingulum, minor forceps, and tracts of corona-radiata (low specificity). Figure Figure4a4a shows the results of our implementation of their method. In this figure, we tuned the speed threshold so that the whole structure is segmented. Figure Figure4b4b demonstrates high sensitivity (segmenting almost the whole structure) and high specificity (without major penetration into the neighbouring structures) of our proposed method compared to the Jonasson's method.

We have also compared the results of the proposed method with those of the Jonasson's method quantitatively. Table Table1 shows1 shows the average Dice measures for different Witelson subdivisions of corpus callosum for the control subjects. Note that the best performance is obtained in the genu area with least penetration inside adjacent fiber bundles and the worst performance is obtained in the posterior mid-body subdivision with penetration mainly inside cingulum and tracts of corona-radiata. For the tumor patients, the performance of the Jonasson's method was poor and thus we do not present them.

Table Table2 shows2 shows the average Dice measure for different Witelson subdivisions of corpus callosum in the control subjects, Group 1 patients (tumor not infiltrating corpus callosum), and Group 2 patients (tumor infiltrating corpus callosum). Dice measures ranging from 90% to 98% indicate that the automatic segmentation results generated for the control subjects and the Group 1 patients are in excellent agreement with the manual segmentation results. For Group 2 patients where the tumors infiltrated corpus callosum, the Dice measures ranging from 81% to 92% indicate good agreements of the automatic and manual segmentation results.

Performance of the proposed method can be ascertained by visualizing the segmentation results (Figures (Figures5,5, ,6,6, ,7).7). Figure Figure55 shows the segmentation results of the proposed method for 11 control subjects. Figure Figure66 shows the results for 12 Group 1 patients where the glioblastoma tumor does not infiltrate corpus callosum. These 12 patients were selected based on the severity of the effect of tumor on the rotation of corpus callosum. Note that although these cases include major geometric deviations from the normal state, the proposed method has successfully segmented corpus callosum and its subdivisions. Figure Figure77 shows the segmentation results for 7 Group 2 patients where the glioblastoma tumor has infiltrated corpus callosum. This figure illustrates that the proposed method has successfully segmented corpus callosum and its subdivisions in these extreme cases. Note that in all cases, the segmented corpus callosum does not penetrate into the adjacent fiber bundles such as cingulum, tapetum, minor and major forceps, superior longitudinal fasciculus, or tracts of corona radiata.

We have used the geometric information of the diffusivity pattern in corpus callosum to prevent the front from penetrating into the neighboring structures. However, a tumor may change the shape and diffusivity pattern of corpus callosum from its original shape and diffusivity pattern.

To study the rotational effect of a tumor on the segmentation process, we rotate the segmented corpus callosum and its neighbors up to 5 voxels in different Euler angles (*azimuth*-*elevation*-*skew*), the tensor (T), and the principal diffusion direction (PDD) of each voxel for all control subjects.

$$R\left(\phi ,\theta ,\psi \right)={R}_{x}\left(\phi \right){R}_{y}\left(\theta \right){R}_{z}\left(\psi \right)$$

(4)

where R(*,θ,ψ*) is the rotation matrix with respect to the spherical Euler angles of *,θ *and *ψ *respectively. R _{x }(), *R _{y }*(

$$\left({x}^{\prime},{y}^{\prime},{z}^{\prime}\right)=R{\left(\phi ,\theta ,\psi \right)}^{*}\left(x,y,z\right)$$

(5)

where (x, y, z) is the coordinates for a voxel within a 5-voxel neighborhood of corpus callosum, while (x', y'z') is the rotated voxel coordinates.

For each of the rotated coordinates, we calculate the tensor (T) and the principal diffusion direction (PDD) by:

$$T\left({x}^{\prime},{y}^{\prime},{z}^{\prime}\right)={\text{R}}^{\mathsf{\text{t}}}\mathsf{\text{(}}\varphi \mathsf{\text{,}}\theta \mathsf{\text{,}}\psi \text{)*T(x,y,z)*R(}\varphi \mathsf{\text{,}}\theta \mathsf{\text{,}}\psi \mathsf{\text{)}}$$

(6)

$$PDD\left({x}^{\prime},{y}^{\prime},{z}^{\prime}\right)=\text{R(}\varphi \mathsf{\text{,}}\theta \mathsf{\text{,}}\psi \text{)*PDD(x,y,z)}$$

(7)

The means and standard deviations of the Dice measures for the subdivisions of the corpus callosum rotated in different Euler angles (5, 10, 15, 20, 25, and 30 degrees) in 11 control subjects are shown in Table Table3.3. The top, middle, and bottom rows represent the segmentation results of the rotated corpus callosum and its neighbors under different azimuth angles (around *x *axis), elevation angles (around *y *axis), and skew angles (around *z *axis), respectively. As can be seen from the Tables Tables22 and and3,3, the rotation of corpus callosum does not have a major effect on the segmentation results. It means that although the segmentation algorithm relies on the geometric information of corpus callosum and its diffusion pattern, it segments a corpus callosum deformed under the potential tumor pressure.

Means and standard deviations of the Dice measures in 11 normal subjects where corpus callosum subdivisions are rotated under different Euler angles

We used the Wilcoxon two sample tests to compare sensitivity of the outer subdivisions (rostrum, genu, and splenium) and the inner subdivisions (rostral body, anterior and posterior mid-body, and isthmus) of corpus callosum with respect to the rotation. For two arrays *A *and *B*, the Wilcoxon test performs a paired two-sided signed rank test of the null hypothesis that data in the vector *A-B *come from a symmetric distribution with zero median. P-values less than 0.01 are considered statistically significant. From the p-values in Table Table4,4, it can be inferred that the outer subdivisions have significantly lower average Dice measures than the inner subdivision for all azimuth and elevation rotation angles (*p*-values less than 0.0004). For skew angles, however, the outer subdivisions do not show significantly lower average Dice measures than the inner subdivision (5 out of 6 p-values are more than 0.01). Considering the overall diffusivity inside corpus callosum perpendicular to its main axis, the overall diffusivity pattern changes more dramatically around the *x *and *y *axes (azimuth and elevation) than around the *z *axis (skew). Therefore, the segmentation results are more sensitive to rotations in the azimuth and elevation directions. Figures Figures8,8, ,99 and and1010 show the segmentation results of the proposed method for the corpus callosum and its neighbors rotated 30° in the azimuth, elevation, and skew directions, respectively.

The *p*-value of Wilcoxon two-sample tests between outer (rostrum, genu, and splenium) and inner (rostral body, anterior and posterior mid-body, and isthmus) subdivisions of corpus callosum

The proposed method is more accurate than the tensor-based method that uses the NTSP similarity measure [25]. Comparing with the manual segmentation by experts, we demonstrated accuracy of our proposed method even when the tumor infiltrates corpus callosum.

The proposed method works in 3D. Extraction of diffusion indices such as mean diffusivity and fractional anisotropy over the entire corpus callosum generates a reliable quantification of the structure that cannot be achieved by an analysis of the mid-sagittal plane only [4-7].

An optimal set of parameters is chosen to segment the entire corpus callosum (with high sensitivity) without penetrating into adjacent fiber structures (with high specificity). We have used *Collinearity_Threshold *to prevent the front from penetrating into cingulum. The threshold for the x-component of the principal diffusion direction (*PDD _{x}_Threshold*) is used to prevent the front from propagating into major and minor forceps.

The proposed method defines Witelson subdivisions of corpus callosum automatically.

The proposed method has successfully segmented corpus callosum and its subdivisions in diffusions MRI data of normal subjects and brain tumor patients.

Depending on size, shape, type and proximity to the corpus callosum, the glioblastoma tumor may cause rotation, shrinkage or more severe disruptions like tearing of the corpus callosum fibers. Amongst the mentioned effects, we evaluated the linear rotational effect of tumor on corpus callosum. With this simplification and without loss of generality, we showed that the change in diffusivity pattern due to a tumor does not change the segmentation accuracy dramatically. However, if the effect of tumor infiltration inside the corpus callosum is sever and changes the structure dramatically, it may be impossible to segment the structure entirely and accurately.

The proposed segmentation method is seed-based and needs the initial seed points for hyper-surface to propagate. This requires the operators to define the seeds manually using anatomical landmarks. To automate the process, the initial seeds may be obtained from fiber atlases of control subjects [34,44]. However, when segmenting fiber bundles of pathological cases, abnormalities like tumors may change the fiber bundles and thus conventional atlas registration methods may not be applicable. To solve this problem, Zacharaki et al. [45] proposed a method for transferring structural and functional information from neuro-anatomical brain atlases into individual patient's data. Application of such methods in our case of segmenting corpus callosum in patients with tumors would be still in question.

Since the tensor model is not capable of describing heterogeneous diffusion behavior in crossing fiber bundles, some studies segmented the desired fiber bundles using HARDI data. Generalizing the level-set method presented in [22] to the HARDI data, Descoteaux and Deriche [35] applied a region-based statistical surface evolution to the image of the ODFs to find coherent white matter fiber bundles. This is equivalent to the maximization of a posteriori probability which obtains the desired segmentation for the observed ODFs and presumes Gaussian distributions in different partitions of the Q-ball images. Although their method appropriately propagates through the regions of fiber crossings, it propagates the hyper-surface inside all of the crossing fiber bundles and segments them as a whole, not individually. They applied their method to extract corpus callosum and tracts of corona-radiata. However, their method is sensitive to initial seeds and suffers from limited anatomical knowledge of the operator who defines the seeds. In another study, a k-means algorithm has been employed to find the clusters using Euclidean distance as dissimilarity measure [36].

In a different study, a Position-Orientation Space (POS) is introduced by combining the geometric space with the spherical ODF space from the HARDI data, where two crossing fiber populations with different orientations in the spatial domain are resolved by applying a front propagation method in the level-set framework in a five-dimensional space [32]. Hagmann et al. [36] performed fiber bundle segmentation in POS based on the Markov Random Fields (MRF). In a similar study, McGraw et al. [37] proposed a mixture of von Mises-Fisher distributions to model the ODF again in MRF. Assuming that the spatial relationships are modeled by the MRF, this method estimates a hidden random field of fiber bundles from the observed ODF profiles using a Maximum a Posteriori (MAP) formulation. However, dimensional reinforcement of the problem causes disadvantages such as increasing the computational cost.

Using Spherical Harmonic Coefficients (SHC) as features of functions on the sphere [38], a method has been proposed for fiber bundle segmentation [33]. However, without masking the speed term with a measure of anisotropy (such as FA) that has low values in the crossing areas, the growing hyper-surface may penetrate into irrelevant fiber bundles that have common areas with the bundle of interest. In another work [34], the authors proposed an atlas-based method introducing a novel similarity based on PDD's and spherical harmonic coefficients. Integrating PDD's into the framework, along with a proper PDD-selection algorithm, leads to the segmentation of most of important fiber bundles in the brain without penetration into irrelevant fiber bundles that have common areas with the bundle of interest. Using an atlas [44] to find the initial seeds for the fiber bundles, the proposed method overcomes limitations of the semi-automatic methods that suffer from limited anatomical knowledge and subjectivity of the operator who defines the seed voxels. Note that since the number of diffusion measurements in the tensor data is not adequate for fitting SHC and estimating more than one PDD, we used the same methods originally applied on the HARDI data.

Although the success of recent studies in segmenting fiber bundles in the HARDI data is promising, such protocol in not being widely used in clinical centers. On the other hand, diffusion tensor imaging is clinically feasible and thus the tensor-based methods are of more interest.

In this paper, we have proposed a 3D method based on a new DTI similarity measure to segment corpus callosum and determine its subdivisions. The method propagates a hyper-surface within corpus callosum without penetration into the neighboring fiber bundles. Segmentation of corpus callosum, the largest commissural fiber bundle in the brain, makes it possible to quantify various diffusion characteristics in its subdivisions, opening new perspectives for monitoring disease evolution or prognosis.

MRI and other data of the tumor patients originally acquired for patient care are retrospectively and anonymously used in this research to train, test, and evaluate the proposed methods. This usage was reviewed and approved by the Institutional Review Board (IRB) committee of Henry Ford Hospital, Detroit, Michigan, USA. Details of the data and results are described in the manuscript.

The authors declare that they have no competing interests.

MN implemented the algorithms presented in the manuscript, tested the algorithms, and prepared the first draft of the manuscript. RJ participated in the design and coordination of the study, helped with the testing and evaluation of the proposed method, and revised the manuscript. SS contributed to the manual segmentation, revision of the manuscript, and final evaluation of the results. AB performed the statistics on the segmentation results. HS contributed to the revision of the manuscript and analysis of the results. MR, QJ, TM contributed to the image acquisition, selection of the cases, and classifying them into different groups. All authors read and approved the final manuscript.

The pre-publication history for this paper can be accessed here:

- de Lacoste MC, Kirkpatrick JB, Ross ED. Topography of the human corpus callosum. J Neuropathol Exp Neurol. 1985;44:578–591. doi: 10.1097/00005072-198511000-00004. [PubMed] [Cross Ref]
- Clarke S, Miklossy J. Occipital cortex in man: organization of callosal connections, related myelo- and cytoarchitecture, and putative boundaries of functional visual areas. J Comp Neurol. 1990;298:188–214. doi: 10.1002/cne.902980205. [PubMed] [Cross Ref]
- Dorion A, Chantôme M, Hasboun D. et al. Hemispheric asymmetry and corpus callosum morphometry: a magnetic resonance imaging study. Neuroscience Research. 2000;36:9–13. doi: 10.1016/S0168-0102(99)00102-9. [PubMed] [Cross Ref]
- Hasan K, Gupta R, Santos R, Wolinsky J, Narayana P. Diffusion tensor fractional anisotropy of the normal-appearing seven segments of the corpus callosum in healthy adults and relapsing-remitting multiple sclerosis patients. Journal of Magnetic Resonance Imaging. 2005;21:735–743. doi: 10.1002/jmri.20296. [PubMed] [Cross Ref]
- Gupta R, Saksena S, Hasan K. et al. Focal Wallerian degeneration of the corpus callosum in large middle cerebral artery stroke: serial diffusion tensor imaging. Journal of Magnetic Resonance Imaging. 2006;24:549–555. [PubMed]
- Kubicki M, Styner M, Bouix S. et al. Reduced interhemispheric connectivity in schizophrenia-tractography based segmentation of the corpus callosum. Schizophrenia research. 2008;106:125–131. doi: 10.1016/j.schres.2008.08.027. [PMC free article] [PubMed] [Cross Ref]
- Rosas HD, Lee SY, Bender AC, Zaleta AK, Vangel M, Yu P, Fischl B, Pappu V, Onorato C, Cha JH, Salat DH, Hersch SM. Altered white matter microstructure in the corpus callosum in Huntington's disease: Implications for cortical. Neuroimage. 2010;49:2995–3004. doi: 10.1016/j.neuroimage.2009.10.015. [PubMed] [Cross Ref]
- Lundervold A, Duta N.
*Taxt T*. Jain A: Model-guided segmentation of corpus callosum in MR images. Citeseer. 1999. pp. 231–237. - Lee C, Huh S, Ketter T, Unser M. Automated segmentation of the corpus callosum in midsagittal brain magnetic resonance images. Optical Engineering. 2000;39:924–935. doi: 10.1117/1.602449. [Cross Ref]
- Brejl M, Sonka M. Object localization and border detection criteria design in edge-based image segmentation: automated learning from examples. IEEE Transactions on Medical Imaging. 2000;19:973–985. doi: 10.1109/42.887613. [PubMed] [Cross Ref]
- Basser P, Mattiello J, Le Bihan D. MR diffusion tensor spectroscopy and imaging. Biophys. 1994;66:259–267. doi: 10.1016/S0006-3495(94)80775-1. [PubMed] [Cross Ref]
- Basser P, Mattiello J, Le Bihan D. Estimation of the effective self-diffusion tensor from the NMR spin echo. J Magn Reson. 1994;103:247–254. doi: 10.1006/jmrb.1994.1037. [PubMed] [Cross Ref]
- Brun A, Knutsson H, Park H, Shenton M, Westin C. Clustering fiber traces using normalized cuts. MICCAI LNCS. 2004;3216:368–375. [PMC free article] [PubMed]
- Maddah M, Grimson W, Warfield S, Wells W. A unified framework for clustering quantitative analysis of white matter fiber tracts. Med Imag Analysis. 2008;12:191–202. doi: 10.1016/j.media.2007.10.003. [PMC free article] [PubMed] [Cross Ref]
- Xu Q, Anderson A, Gore J, Ding Z. Unified Bundling and Registration of Brain White Matter Fibers. IEEE Transactions on Medical Imaging. 2009;28(9):1399–1411. [PMC free article] [PubMed]
- Voineskos A, Rajji T, Lobaugh N, Age-related decline in white matter tract integrity and cognitive performance: A DTI tractography and structural equation modeling study. Neurobiology of Aging. 2010. in press doi:10.1016/j.neurobiolaging.2010.02.009. [PMC free article] [PubMed]
- Wassermann D, Bloy L, Kanterakis E, Verma R, Deriche R. Unsupervised white matter fiber clustering and tract probability map generation: Applications of a Gaussian process framework for white matter fibers. Neuroimage. 2010;51(1):228–241. doi: 10.1016/j.neuroimage.2010.01.004. [PMC free article] [PubMed] [Cross Ref]
- Li H, Xue Z, Guo L, Liu T, Hunter J, Wong S. A hybrid approach to automatic clustering of white matter fibers. Neuroimage. 2010;49:1249–1258. doi: 10.1016/j.neuroimage.2009.08.017. [PubMed] [Cross Ref]
- Eckstein I, Shattuck D, Stein J. et al. Active fibers: Matching deformable tract templates to diffusion tensor images. Neuroimage. 2009;47:82–89. [PMC free article] [PubMed]
- Zhukov L, Museth K, Breen D, Whitaker R, Barr A. Level set modeling and segmentation of DT-MRI brain data. J Electron Imag. 2003;12:125–133. doi: 10.1117/1.1527628. [Cross Ref]
- Vemuri B, Chen Y, McGraw T, Wang Z, Mareci T. Fiber tract mapping from diffusion tensor MRI. Proc IEEE Workshop on Variational and Level Set Methods in Computer Vision. 2001. pp. 81–88.
- Wang Z, Vemuri B. Tensor field segmentation using region based active contour model. ECCV LNCS. 2004;3024:304–315.
- Rousson M, Lenglet C, Deriche R. Level set and region based surface propagation for diffusion tensor MRI segmentation. Computer vision and mathematical methods in medical and biomedical image analysis. 2004. pp. 123–134.
- Feddern C, Weickert J, Burgeth B. Level set methods for tensor valued images. 2003. pp. 65–72. Geometric and Level Set Methods in Computer Vision: In Proceedings of the Second IEEE Workshop on Variational.
- Jonasson L, Bresson X, Hagmann P, Cuisenaire O, Meuli R, Thiran J. White matter fiber tract segmentation in DT-MRI using geometric flows. Med Imag Analysis. 2005;3:223–236. [PubMed]
- Jonasson L, Hagmann P, Pollo C, Bresson X, Wilson C, Meuli R, Thiran J. A level set method for segmentation of the thalamus and its nuclei in DT-MRI. Signal Processing. 2007;87:309–321. doi: 10.1016/j.sigpro.2005.12.017. [Cross Ref]
- Lenglet C, Rousson M, Deriche R. A statistical framework for DTI segmentation. IEEE Trans Med Imag. 2006;25:675–700.
- Arsigny V, Fillard P, Pennec X, Ayache N. Log-Euclidean metrics for fast and simple calculus on diffusion tensors. Magnetic Resonance in Medicine. 2006;56:411–421. doi: 10.1002/mrm.20965. [PubMed] [Cross Ref]
- Weldeselassie Y, Hamarneh G. DT-MRI segmentation using graph cuts. Citeseer: 65121 K-65121.
- Awate S, Hui Z, Gee J. A fuzzy, nonparametric segmentation framework for DTI and MRI analysis: With applications to DTI-tract extraction. IEEE Transactions on Medical Imaging. 2007;26:1525–1536. [PubMed]
- Hasan K, Ewing-Cobbs L, Kramer L, Fletcher J, Narayana P. Diffusion tensor quantification of the macrostructure and microstructure of human midsagittal corpus callosum across the lifespan. NMR in Biomedicine. 2008;21:1094–1101. doi: 10.1002/nbm.1286. [PMC free article] [PubMed] [Cross Ref]
- Jonasson L, Bresson X, Thiran J, Wedeen V, Hagmann P. Representing diffusion MRI in 5-D simplifies regularization and segmentation of white matter tracts. IEEE Trans On Med Imaging. 2007;26:1547–1554. [PubMed]
- Nazem-Zadeh M, Davoodi-Bojd E, Soltanian-Zadeh H. Level set fiber bundle segmentation using spherical harmonic coefficients. Computerized Medical Imaging and Graphics. 2010;34:192–202. doi: 10.1016/j.compmedimag.2009.09.003. [PubMed] [Cross Ref]
- Nazem-Zadeh M, Davoodi-Bojd E, Soltanian-Zadeh H. Atlas-Based Fiber Bundle Segmentation Using Principal Diffusion Directions and Spherical Harmonic Coefficients. NeuroImage. 2011;54:246–164. Supplement 1. [PubMed]
- Descoteaux M, Deriche R. High angular resolution diffusion MRI segmentation using region-based statistical surface evolution. Journal of Mathematical Imaging and Vision. 2008;33(2):239–252.
- Hagmann P, Jonasson L, Deffieux T, Meuli R, Thiran J, Wedeen V. Fibertract segmentation in position orientation space from high angular resolution diffusion MRI. Neuroimage. 2006;32:665–675. doi: 10.1016/j.neuroimage.2006.02.043. [PubMed] [Cross Ref]
- McGraw T, Vemuri B, Yezierski R, Mareci T. Segmentation of High Angular Resolution Diffusion MRI modeled as a field of von mises-fisher mixtures. Proc European Conf on Computer Vision. 2006. pp. 463–475.
- Hess C, Mukherjee P, Han E, Xu D, Vigneron D. Q-ball reconstruction of multimodal fiber orientations using the spherical harmonic basis. Magn Reson Med. 2006;56:104–117. doi: 10.1002/mrm.20931. [PubMed] [Cross Ref]
- Witelson S. Hand and sex differences in the isthmus and genu of the human corpus callosum: a postmortem morphological study. Brain. 1989;112:799. doi: 10.1093/brain/112.3.799. [PubMed] [Cross Ref]
- Lee S, Cheng J, Chen C, Tseng W. An automatic segmentation approach for boundary delineation of corpus callosum based on cell competition. Proc Engineering in Medicine and Biology Society Conf of the IEEE. 2008. pp. 5514–5517. [PubMed]
- Osher S, Sethian J. Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulations. J Comput Phys. 1988;79:12–49. doi: 10.1016/0021-9991(88)90002-2. [Cross Ref]
- Osher S, Paragios N. Chapter: Level set methods. Geometric level set methods in imaging, vision, and graphics. Springer-Verlag New York Inc. 2003. pp. 3–20.
- Dice LR. Measures of the amount of ecologic association between species. Ecology. 1945;26:297–302. doi: 10.2307/1932409. [Cross Ref]
- Wakana S, Caprihan A, Panzenboeck M. et al. Reproducibility of quantitative tractography methods applied to cerebral white matter. Neuroimage. 2007;36:630–644. doi: 10.1016/j.neuroimage.2007.02.049. [PMC free article] [PubMed] [Cross Ref]
- Zacharaki E, Hogea C, Shen D, Biros G, Davatzikos C. Non-diffeomorphic registration of brain tumor images by simulating tissue loss and tumor growth. Neuroimage. 2009;46:762–774. doi: 10.1016/j.neuroimage.2009.01.051. [PMC free article] [PubMed] [Cross Ref]

Articles from BMC Medical Imaging are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |