 Research
 Open Access
 Published:
Segmentation of abdomen MR images using kernel graph cuts with shape priors
BioMedical Engineering OnLine volumeÂ 12, ArticleÂ number:Â 124 (2013)
Abstract
Background
Abdominal organs segmentation of magnetic resonance (MR) images is an important but challenging task in medical image processing. Especially for abdominal tissues or organs, such as liver and kidney, MR imaging is a very difficult task due to the fact that MR images are affected by intensity inhomogeneity, weak boundary, noise and the presence of similar objects close to each other.
Method
In this study, a novel method for tissue or organ segmentation in abdomen MR imaging is proposed; this method combines kernel graph cuts (KGC) with shape priors. First, the region growing algorithm and morphology operations are used to obtain the initial contour. Second, shape priors are obtained by training the shape templates, which were collected from different human subjects with kernel principle component analysis (KPCA) after the registration between all the shape templates and the initial contour. Finally, a new model is constructed by integrating the shape priors into the kernel graph cuts energy function. The entire process aims to obtain an accurate image segmentation.
Results
The proposed segmentation method has been applied to abdominal organs MR images. The results showed that a satisfying segmentation without boundary leakage and segmentation incorrect can be obtained also in presence of similar tissues. Quantitative experiments were conducted for comparing the proposed segmentation with other three methods: DRLSE, initial erosion contour and KGC without shape priors. The comparison is based on two quantitative performance measurements: the probabilistic rand index (PRI) and the variation of information (VoI). The proposed method has the highest PRI value (0.9912, 0.9983 and 0.9980 for liver, right kidney and left kidney respectively) and the lowest VoI values (1.6193, 0.3205 and 0.3217 for liver, right kidney and left kidney respectively).
Conclusion
The proposed method can overcome boundary leakage. Moreover it can segment liver and kidneys in abdominal MR images without segmentation errors due to the presence of similar tissues. The shape priors based on KPCA was integrated into fully automatic graph cuts algorithm (KGC) to make the segmentation algorithm become more robust and accurate. Furthermore, if a shelter is placed onto the target boundary, the proposed method can still obtain satisfying segmentation results.
Introduction
The recent development of open magnetic resonance imaging has provided new opportunities for next generation imageguided surgical and interventional applications. Imageguided surgery is a standard surgical procedure for abdominal disorders that can reduce surgical trauma and open surgery burdens [1]. However, during surgical planning and surgical navigation based on MR images, there are two problems have to be faced: the shape deformation of the organs and the similarity among abdominal organs. For these reasons, an effective and robust algorithm for abdominal organs segmentation is helpful and very important in imageguided surgery and surgical navigation system [2, 3].
Usually, basic image information, such as intensity and gradient, are used for segmentation; Ostu, kmeans clustering, region growing, etc. are the most widely used algorithms for segmentation. However, they are not suitable for MR segmentation. Because of its weak boundary, intensity inhomogeneity and noise, the segmentation of MR images is considered a complex procedure [4]. For this reason, developing different advanced and intelligent algorithms for MR image segmentation has become a research hotspot over the last few years.
Abdomen MR image segmentation is a challenging task, because majority of tissues in abdomen are soft tissues, the intensities of abdominal tissues are very similar and shape change in a complex way due to respiratory movements [5]. To provide more information about the tissues and organs in abdomen to the doctors, more effective and robust algorithm to segment abdomen MR images are developed. For example, Hassan et al. [6] proposed a novel method to segment liver MR images automatically. This algorithm utilizes artificial neural networks and watershed algorithm. Moreover, Sheng et al. [7] applied a waveletbased kmeans clustering method to segment the human kidney from MRI data set.
Since the level set method can be performed on a fix Cartesian grid with no need to parameterize these objects, in the past decade it was used more and more frequently in image segmentation [8â€“11]. Moreover, the level set method can represent contours with complex topology and change their topology in a natural way. However, the algorithms converge to a local minimum easily, so the results are easily affected by the initial values. Additionally, the execution time may be very long in some applications, especially with large images and multiobject segmentations [12, 13].
Graph cuts techniques have received considerable attention for their global energy optimal advantages. It is used in more and more image segmentation applications for different medical images, such as MR images [14]. However, it requires to choose the object and background seeds interactively, implying a timeconsuming procedure. Kernel graph cuts is a fully automatic algorithm based on graph cuts proposed by Salah et al. [15]. It consists of a multiregion image segmentation algorithm based on graph cuts via kernel mapping of the image data. This algorithm is not successful for abdominal organs segmentation due to the weak boundary and surrounding objects with similar intensities. The graph cuts based on active contours method was proposed by Xu et al. [16]. The method identifies an initial contour around the target and then forms the narrow banded area containing the target via dilation and erosion. The segmentation focuses on image data in the narrow banded area, and separates the abdominal organs with similar intensities. However it cannot solve the weak boundary problem that leads to edge leakage. Integrating the prior information into the segmentation is a popular solution for the weak boundary; the prior information can lead to a more accurate segmentation result according to [17â€“20]. Asem et al. [21] proposed graphs cuts integrated shape priors, and applied this method to kidney segmentation in abdomen MR images. The result of the segmentation shows that it can overcome the edge leakage, however, it uses the probabilistic model, which increases its complexity; moreover, it also needs interactive operations. Chen et al. [21] integrate shape information generated from Active Appearance Model (AAM) into graph cuts for abdominal 3D organ segmentation. However, the long computation time and the dependence of initial location limit its applications. Comparing to the linear PCA (principle component analysis), the nonlinear PCA model performs better on problems with nonlinear deformation [22]. Malcolm et al. [23] proposed a segmentation model that integrates the KPCAbased shape priors into graph cuts, but the model executes iteratively so it is timeconsuming.
In this paper, we propose a new method, which combines kernel graph cuts with KPCA to segment abdominal organs. First, a seed point is chosen inside the target organ, and then the region growing algorithm is used to obtain approximately segmentation. Second, image morphology operations (dilation and erosion) are used to form an initial contour close to the target organ. Third, KPCA is used to obtain the shape priors by training the relevant shape templates set. At last, the shape priors are integrated into kernel graph cuts to make a better segmentation. The main contribution of this paper is the proposal of novel methods that combine kernel graph cuts algorithm with shape priors for abdominal organs extraction. The shape priors based on KPCA help our algorithm to increase segmentation accuracy while intensity is not sufficient to obtain accurate segmentation results.
This paper is organized as follows. Section Methodology describes the methodology which gives a description of kernel graph cuts and shape priors based on KPCA and detailed introduction of our proposed method. Section Experiments presents the experiments results obtained using the proposed novel method. Finally, the discussion about the algorithms is presented and conclusions are drawn in Section Discussion and Section Conclusions, respectively.
Methodology
This section starts by briefly describing the KGC and KPCA.
Kernel graph cuts
Graph cuts algorithm was introduced by Boykov et al. [14] for binary image segmentation application. The purpose is to segment an object from a given image using a set of seeds (object and background) placed by user.
The graph cuts algorithm aims to cast the energybased image segmentation problem into a graph structure global mincut problem. The energy function of graph cuts contains two terms: a regionbased term R(A) and a boundary term B(A), where A stands for an object or background pixel assignment. The regionbased term evaluates the penalty for assigning a particular pixel to a given region. The boundary term evaluates the penalty for assigning two neighboring pixels to different regions. These two terms often weight by 0â€‰â‰¤â€‰Ïµâ€‰â‰¤â€‰1 for relative influence, and the energy function is expressed as follows:
KGC was proposed by Salah et al. [15] for automatic segmentation by mapping image data into high dimension through kernel function. Graph cuts method is a supervised algorithm which requires user intervention for choosing seeds (object & background). The proposed energy function contains two terms: an original kernelinduced data term which evaluates the deviation of the mapped image data and a regularization term expressed as a function of the region indices. The energy function is expressed as follows:
Where E({Î¼_{1}},Î´) measures kernelinduced non Euclidean distances between the observations and the regions parameters Î¼_{ 1 }. Ï† is a nonlinear mapping from the observation space I to a higher dimensional mapped space J, and the radial basis function (gauss function) kernel is used as usual. a is a positive factor. Î´ is an indexing function which assigns each point of the image to a region. lâ€‰âˆˆâ€‰L Is a pixel label in some finite set of labels L. PâˆˆR_{ 1 } is a pixel in each region which is characterized by one label l. R(Î´((p),Î´(q,)) is a smooth regularization function and D is a neighborhood set containing all pairs of neighboring pixels {p,q}.
According to the Mercerâ€™s theorem [24], which states that kernel function can be expressed as a dot product in a highdimensional space, explicitly the mapping Ï† is not available. Instead, the kernel function is as follow:
Substitution of the kernel function gives:
Equation (4) is solved forâˆ¥Ï•(I_{ p })â€‰âˆ’â€‰Ï•(Î¼)âˆ¥^{2} and substituted in (2). Thus, the kernelinduced energy function is given by:
Now, based on the kernelinduced energy function, the graph cuts algorithm can be executed to segment images more efficiently.
Kernel principle component analysis
KPCA is a nonlinear feature extractor performed in the feature space F. The basic idea of the method is to map the data from the input space S to a feature space F via nonlinear map Ï•:Sâ†’F. Because KPCA is able to capture nonlinear features in the data comparing to linear PCA, it can be used more effectively if a preimage of the projection in the feature space is available. Rathi et al. [22] proposed a novel method to reconstruct a unique approximate preimage of a feature vector and applied it for statistical shape analysis.
To form the statistical model of shape space S, the preimage of the projection (in the KPCA space) of a test point xâˆˆS should be found, as shown in Figure 1 from [25]. Let {x_{ 1 }, X_{ n }}âŠ‚S be a set of aligned training shapes represented by binary mask where 1 is object and 0 is background and spread as vectors.
First, the NÃ—N kernel matrix K with the radial basis function (gauss function) has to be computed:
Second, following Eigen decomposition has to be considered :
H is the centering matrix given by H=C\xe2\u02c6\u2019\frac{1}{N}c{c}^{T},C is the NÃ—N identity matrix, câ€‰=â€‰[11â€‰â€¦â€‰1]^{T} is an NÃ—1 vector, Uâ€‰=â€‰[âˆ‚_{1},â€‰â‹¯â€‰âˆ‚_{ N }]^{T} with âˆ‚_{ i }â€‰=â€‰[a_{i 1},â€‰â€¦,â€‰a_{ iN }]^{T} is the matrix containing the eigenvectors and Î£â€‰=â€‰diag(Î»_{1},â€‰â€¦,â€‰Î»_{ N }) contains the corresponding eigen values.
Third, given a point xâˆˆS one can compute its projection PÏ•(x)âˆˆF and a subspace is spanned by the first n eigenvectors given by:
Where \stackrel{\xc2\xaf}{\mathrm{\xcf\u2020}}=\frac{1}{N}{\displaystyle {\mathrm{\xe2\u02c6\u2018}}_{i=1}^{N}\mathrm{\xcf\u2020}}\left({x}_{i}\right); \stackrel{\xcb\u0153}{\mathrm{\xcf\u2020}} is the map centralized by \stackrel{\xcb\u0153}{\mathrm{\xcf\u2020}}\left(x\right)=\mathrm{\xcf\u2020}\left(x\right)\xe2\u02c6\u2019\stackrel{\xc2\xaf}{\mathrm{\xcf\u2020}}. {V}_{k}={\displaystyle \underset{i=1}{\overset{N}{\mathrm{\xe2\u02c6\u2018}}}\frac{{a}_{k\phantom{\rule{0.1em}{0ex}}i}}{\sqrt{{\mathrm{\xce\xbb}}_{k}}}}\stackrel{\xcb\u0153}{\mathrm{\xcf\u2020}}\left({x}_{i}\right) is the k th orthogonal eigenvector of the covariance matrix in F. The projection of test point x in F project onto the k th component by B_{ k }. Then {\mathrm{\xce\xb2}}_{k}=\frac{1}{\sqrt{{\mathrm{\xce\xbb}}_{k}}}{\displaystyle \underset{i=1}{\overset{N}{\mathrm{\xe2\u02c6\u2018}}}{a}_{k\phantom{\rule{0.1em}{0ex}}i}}\stackrel{\xcb\u0153}{k}\left(x,{x}_{i}\right)
Finally, the method in [22] is used to compute the approximate preimage \widehat{x}:
And the distance used in feature space in equation (9) is defined as follows:
where M={\displaystyle \underset{i=1}{\overset{N}{\mathrm{\xe2\u02c6\u2018}}}\frac{1}{{\mathrm{\xce\xbb}}_{i}}}{\xe2\u02c6\u201a}_{k}{{\xe2\u02c6\u201a}_{k}}^{T}, Î³â€‰=â€‰[âˆ‚_{1},â€‰â€¦,â€‰âˆ‚_{ n }]Î², \stackrel{\xcb\u0153}{\mathrm{\xce\xb3}}=\mathrm{\xce\xb3}+\frac{1}{N}\left(1\xe2\u02c6\u2019{c}^{T}\mathrm{\xce\xb3}\right).
The preimage \widehat{x} contains the information that is used to form the shape priors. In following section, the proposed method is presented and summarized.
Proposed method
To mitigate the effect of weak boundary in MR images and to segment abdominal organs from surrounding objects with similar intensities, a novel segmentation method is proposed. The method incorporates KPCA with KGC enlightened by [23]. KPCA is used to form the shape priors based on the statistical model proposed by Rathi et al. [22]. Enlightened by the GCBAC (graph cuts based Active Contours) proposed by Xu et al. [16], the initial dilation and erosion contour idea is introduced into the method to locate the position of shape priors. The framework of the segmentation method is shown in Figure 2.
The proposed segmentation framework can be described in two phases.
Phase I is the presegmentation phase. It includes two procedures to form the image data.
First, to form the initial contour, a seed point is chosen inside the target region manually, and then the region growing algorithm is used to segment the target region approximately. Based on the result of region growing algorithm, a lot of isolated small regions in the target region are not segmented correctly. For this reason, the morphological dilate and erode operation are executed over the presegmented contour using the region growing algorithms. The dilate operation aims to eliminate those isolated small regions and to form a continuous contour around the target region; the erode operation aims to draw the contour near the real target region after the dilate operation. Same morphological structuring element is used for the operation of erode and dilate. As a result, the initial contour is obtained. However, because of the weak boundary in the MR images, a further processing is needed to overcome this problem to reach a more accurate segmentation.
Second step includes obtaining the shape priors. Because of the differences varying from person to person, and different parameters of MR imaging, the shapes of abdominal organs to be segmented differ from each other. Moreover, the respiratory movements make the deformation more complex. Since the deformation is nonlinear, the KPCA is used to train the shape templates determined by experts. Before training the data set, all the shape templates and the initial contour should be aligned. So the image registration is needed here, and only translating, scaling and rotating transforms are taken into consideration during the registration process. Suppose that X_{ i } is the vector of one shape template, one of the shape templates Xj is to be chosen, the fittest parameters can be obtained through translating tj, scaling and rotating transform M(s_{ j },Î¸_{ j }). The target function is obtained by minimizing the error measure Ej:
W is the weight matrix, and the least square method is usually used to solve it.
Commonly, most applications, such as feature extraction and pattern classification, only need the new features generated by KPCA. However, for some other applications, reconstructing the preimage from the KPCA features is needed. In this case KPCA feature is not necessary to describe the deformation patterns; on the contrary, it is required to reconstruct the shapes from the KPCA features. For a Gaussian kernel, the preimage \widehat{x} can be obtained using Equations (610).
Phase II is the segmentation phase. Since the registration is composed of translating, scaling and rotating transform, and the abdominal organs change from patient to patient, the result of registration contour does not represent the real boundary of MR images, and a more accurate segmentation procedure need to be followed. The image data constructs the graph using energy function (5) in KGC, introducing the shape priors \widehat{x} into the data term to overcome the weak boundary. This work aims to use the data term for representing the penalty of pixel attribute to the object or background. Thus, it is assumed that nonuniform shape priors Pp(O) and P_{p}(B) represent penalty of the shape priors attribute to the object or background at a pixel P. A parameter Î·(0â‰¤Î·â‰¤1) is also introduced to represent the weight of relative influence between kernelinduced data term J_{ K } and shape priors, so the new data terms can be written as follows:
Since the preimage has value between 0 and 1, P_{ p }(O) is directly used to represent \widehat{x} and set Ppâ€‰=â€‰( 1 P_{ p }(O)). The smooth term used the original term. At last, based on Equation (5), the new energy function is given by:
Where l must be O or B, which stands for object or background. Thus, the multiregion segmentation of KGC turns to be the binary segmentation of object or background. The new energy function is used to construct graph, compute mincut and get the segmentation.
Experiments
Segmentation is performed on abdominal organs in the abdomen MR images of T1 sequence, and the proposed novel segmentation method is validated using the MATLAB 7.11 program on Windows XP with Lenovo PC with Intel (R) Core (TM) 2 Duo CPU, E7500, and tested using the liver and kidney dataset in abdomen MR images with size of 462â€‰Ã—â€‰310 pixels. All MR images are obtained by Siemens 3.0T MR imaging equipment and all the shape templates are segmented manually by different experts. The size of train set is 30.
Validation
First, the original MR images have to be segmented as shown in Figure 3. Figure 3 (a) is used to segment liver in abdominal organs, and Figure 3 (b) is used to segment both of the two kidneys in abdominal organs.
Second, a seed point is placed inside the liver and kidney region in the MR images, and then the region growing algorithm is used to segment the image approximately. The result of region growing algorithm is shown in Figure 4. From the Figure 4, it can be seen that many isolated small regions are not segmented, and boundary leakage and incorrect segmentation are observed. Many blood vessels exist inside the organ tissues and the noise caused by MR equipment lead to lots of isolated small regions; meanwhile the soft tissues are very similar and the overlaps between different soft tissues lead to weak boundary. So, more work is needed to solve this problem.
Third, the morphological dilate and erode operation are implemented to fix the problem caused by region growing algorithm. The result contour of region growing is dilated and then eroded as demonstrated in Figure 5, the yellow contour is the result of dilation and the blue contour is the result of erosion. The type of morphological structuring element that is used in dilation and erosion is a disk whose radius can be adjusted depending on the result of region growing algorithm. In this experiment, the size of radius is 15 pixels in liver segmentation, and 10 pixels in kidney segmentation. Same morphological structuring element is used in the same MR image during the dilation and the erosion operation. The blue contour is close to the real boundary of abdominal organs as can be seen in Figure 5, but boundary leakage is still observable.
Finally, the shape priors are added into the segmentation to guide the contour more accurate to the real boundary. The size of the training shape template set is 30, and a registration is performed between initial erosion contour (blue contour) and all shape templates (green contour). The result is shown in Figure 6. Figure 7 demonstrates the shape priors which is obtained by training the shape template set using KPCA method. Finally, KGC is combined with shape priors to segment the region inside initial dilation contour (yellow contour) based on function (14). The result of segmentation is shown in Figure 8, from which we can find the boundary leakage is eliminated and the segmentation becomes more correct and accurate.
Quantitative verification
In order to verify the proposed segmentation method, quantitatively experiments are performed to compare our method and other three methods: DRLSE [11], initial erosion contour, KGC in the initial dilation contour without shape priors. The comparison is based on two quantitative performance measures: the probabilistic rand index (PRI) and the variation of information (VoI) [26â€“29].
The PRI counts the fraction of pairs of pixels whose labels are consistent between the computed segmentation and the ground truth. The VoI metric defines the distance between two segmentations like the average conditional entropy. Since the result segmentation boundaries of the proposed method are very close to the ground truth, visible difference between them cannot be identified, for this reason the PRI and VoI are used to quantify the segmentation results.
For each segmentation method, a higher value of PRI and a lower value of VoI imply that the segmentation results are closer to the expert manual segmentation. The statistic data are illustrated in Table 1, and the results of different segmentation algorithms are shown in Figure 9. As can be seen in the Table 1, no matter whether the abdominal organ is liver, right kidney or left kidney, the proposed method has the highest PRI values and the lowest VoI values. As shown in Figure 9, the segmentation with the proposed method has better performance than the other methods. The KGC with shape priors based on KPCA can overcome the boundary leakage and segment every abdominal organ independently without incorrect segmentation of the similar tissues. It indicates that both for liver and for kidney segmentation, the proposed method is better than the other methods.
Parameter adjustment
The value of parameter represents the relative influence of KGC data term and shape priors term. In Figure 8, the value of Î· increases gradually from liver to right kidney. Since the incorrect segmentation in the Figure 9(gi) decreases gradually from liver to right kidney, the KGC data term weighted in the energy function become high. If Î· has a high value, the weight of KGC data term is high too and the weight of shape priors term is small. For different organs different organs, the value of Î· can be adjusted to optimize segmentation performance.
As shown in Figure 10, when =0.65, 0.72 and 0.75, the segmentation of liver, right kidney and left kidney reach the best accuracy. From the third row of Table 1, the value of PRI is 0.8903, 0.9808 and 0.9913, and the value of VoI is 2.1212, 0.5351 and 0.3965 from liver to right kidney when using KGC inside initial dilation contour without shape priors. In another words, the segmentation error decreases from liver to right kidney, as also shown in Figure 9 (gi).
As mentioned above, by adjusting the value of parameter Î·, the segmentation can be optimized. If the incorrect segmentation proportion of KGC inside initial dilation contour without shape priors is small, the weight of KGC data term is high, so that the value of Î· become bigger. On the contrary, if the incorrect segmentation proportion of KGC inside initial dilation contour without shape priors is high, the value of Î· should become smaller. In this way, a satisfying segmentation result can be obtained finally.
Discussion
Because of the noise, weak boundary, intensity inhomogeneity and similar intensities among different abdominal organs, the segmentation of abdominal organs is normally considered a challenging task. Furthermore, the deformation caused by individual difference and respiratory movement makes the segmentation task even more difficult.
KGC is a fully automatic segmentation algorithm based on graph cuts. If the KGC algorithm is used alone to segment abdomen MR image, the segmentation result is not satisfying as shown in the Figure 11 (a). The blood vessels in organ and the overlaps among abdominal organs affect the segmentation, thus, a fully complete boundary of liver or other organs cannot be obtained. To make the segmentation procedure focus on a given organ, the initial contour is produced by region growing algorithm and morphology operation. But this is considered not enough. The segmentation result can be seen from Figure 11(b). For this reason some shape priors has been used to make the segmentation algorithm more robust and accurate. Since KPCA can handle nonlinear deformable information, the shape priors based on KPCA is integrated into KGC. Moreover, if a shelter is placed onto the targetâ€™s boundary, a satisfying segmentation result can be also obtained. This is validated by experiments as shown in Figure 11 (c).
Additionally, to develop reliable prior knowledge, we should choose patients with similar age and weight. In this way their shape prior would guide to the correct segmentation.
On the clinical side, the work assumed that the geometric inaccuracy, such as distortion, does not exist. Because of this, the MR imaging technology has experienced a rapid development to overcome the magnetic field inhomogeneity. Nowadays, the image spatial resolution can be very high and geometric distortion of MR images can be ignored. However, if a more accurate segmentation is needed, geometric distortion should be corrected at first. This will be investigated in future works.
Conclusions
In this paper, a novel method is proposed to segment abdominal organs integrating kernel graph cuts with KPCA shape priors after a series of preprocessing on the abdomen MR images. The morphology operation can eliminate the isolated small region after region growing algorithm. The kernel graph cuts is a fully automatic segmentation algorithm, and it also has global minimization and polynomial time complexity characteristics. The shape priors, generated by preimage projection via KPCA, can handle nonlinear deformation. Experiments on liver and kidney segmentation of abdomen MR image showed that the novel method can obtain satisfying results. The kernel function used in this method is gauss function, and other kernel functions have not been tested yet. Currently, the seed point is obtained manually, but automation will be considered in future works. For the time performance, the algorithm can be parallelized on Graphic processors to achieve higher performance.
Abbreviations
 MR:

Magnetic resonance
 KGC:

Kernel graph cuts
 KPCA:

Kernel principle component analysis
 PRI:

Probabilistic rand index
 VoI:

Variation of information
 DRLSE:

Distance regularized level set evolution
 AAM:

Active appearance model
 PCA:

Principle component analysis
 GCBAC:

Graph cuts based active contours.
References
Tsutsumi N, Tomikawa M, Uemura M, et al.: Imageguided laparoscopic surgery in an open MRI operating theater[J]. Surgical Endoscopy 2013, 27(6):2178â€“2184. 10.1007/s004640122737y
Oliveira DAB, Feitosa RQ, Correia MM: Segmentation of liver, its vessels and lesions from CT images for surgical planning[J]. Biomedical engineering online 2011, 10(1):1â€“23. 10.1186/1475925X101
Peterhans M, Vom Berg A, Dagon B, et al.: A navigation system for open liver surgery: design, workflow and first clinical applications[J]. The International Journal of Medical Robotics and Computer Assisted Surgery 2011, 7(1):7â€“16. 10.1002/rcs.360
Ma Z, Tavares JMRS, Jorge RN, et al.: A review of algorithms for medical image segmentation and their applications to the female pelvic cavity[J]. Comput Methods Biomech Biomed Engin 2010, 13(2):235â€“246. 10.1080/10255840903131878
Artan Y, Haider MA, Yetik IS: Graphbased active contours using shape priors for prostate segmentation with MRI[C]/Biomedical Imaging: From Nano to Macro, 2011. IEEE International Symposium on. IEEE; 2011:1459â€“1462. Chicago, IL; 2011. Date of Conference: March 30 2011April 2 2011
Masoumi H, Behrad A, Pourmina MA, et al.: Automatic liver segmentation in MRI images using an iterative watershed algorithm and artificial neural network[J]. Biomedical Signal Processing and Control 2012, 7(5):429â€“437. 10.1016/j.bspc.2012.01.002
Li S, ZÃ¶llner FG, Merrem AD, et al.: Waveletbased segmentation of renal compartments in DCEMRI of human kidney: initial results in patients and healthy volunteers[J]. Comput Med Imaging Graph 2012, 36(2):108â€“118. 10.1016/j.compmedimag.2011.06.005
Caselles V, Kimmel R, Sapiro G: Geodesic active contours[J]. International journal of computer vision 1997, 22(1):61â€“79. 10.1023/A:1007979827043
Chan TF, Vese LA: Active contours without edges[J]. Image Processing, IEEE Transactions on 2001, 10(2):266â€“277. 10.1109/83.902291
Li C, Kao CY, Gore JC, et al.: Implicit active contours driven by local binary fitting energy[C]. Minneapolis, MN; 2007. Date of Conference: 17â€“22 June 2007
Li C, Xu C, Gui C, et al.: Distance regularized level set evolution and its application to image segmentation[J]. Image Processing, IEEE Transactions on 2010, 19(12):3243â€“3254.
Ben Ayed I, Mitiche A, Belhadj Z: Polarimetric image segmentation via maximumlikelihood approximation and efficient multiphase levelsets[J]. Pattern Analysis and Machine Intelligence, IEEE Transactions on 2006, 28(9):1493â€“1500.
Li C, Huang R, Ding Z, et al.: A level set method for image segmentation in the presence of intensity inhomogeneities with application to MRI[J]. Image Processing, IEEE Transactions on 2011, 20(7):2007â€“2016.
Boykov YY, Jolly MP: Interactive graph cuts for optimal boundary & region segmentation of objects in ND images[C]//Computer Vision, 2001. ICCV 2001. Proceedings. Eighth IEEE International Conference on. IEEE 2001, 1: 105â€“112.
Salah MB, Mitiche A, Ayed IB: Multiregion image segmentation by parametric kernel graph cuts[J]. Image Processing, IEEE Transactions on 2011, 20(2):545â€“557.
Xu N, Ahuja N, Bansal R: Object segmentation using graph cuts based active contours[J]. Comput Vis Image Underst 2007, 107(3):210â€“224. 10.1016/j.cviu.2006.11.004
Uzunbas MG, Zhang S, Pohl KM, et al.: Segmentation of myocardium using deformable regions and graph cuts[C]. Barcelona; Date of Conference: 2â€“5 May 2012
Slabaugh G, Unal G: Graph cuts segmentation using an elliptical shape prior[C]//Image Processing, 2005. Volume 2. II edition. Genoa, Italy: ICIP 2005. IEEE International Conference on. IEEE; 2005:1222â€“1225. September 11â€“14
Wang H, Zhang H, Ray N: Adaptive shape prior in graph cut image segmentation[J]. Pattern Recognition 2013, 46(5):1409â€“1414. 10.1016/j.patcog.2012.11.002
Song Z, Tustison N, Avants B, et al.: Adaptive graph cuts with tissue priors for brain MRI segmentation[C]//Biomedical Imaging: Nano to Macro, 2006. 3rd IEEE International Symposium on. Arlington, VA: IEEE; 2006:762â€“765. Date of Conference: 6â€“9 April 2006
Ali AM, Farag AA, ElBaz AS: Graph cuts framework for kidney segmentation with prior shape constraints[M]//Medical Image Computing and ComputerAssisted Interventionâ€“MICCAI 2007. Volume 4791. Berlin Heidelberg: Springer; 2007:384â€“392.
Chen X, Udupa JK, Bagci U, et al.: Medical image segmentation by combining graph cuts and oriented active appearance models[J]. Image Processing, IEEE Transactions on 2012, 21(4):2035â€“2046.
Rathi Y, Dambreville S, Tannenbaum A: Statistical shape analysis using kernel PCA [J]. San Jose, California, USA: Society for Imaging Science and Technology (IS&T) and SPIE.; 2006.
Malcolm J, Rathi Y, Tannenbaum A: Graph cut segmentation with nonlinear shape priors[C]//Image Processing, 2007. San Antonio, TX: ICIP 2007. IEEE International Conference on. IEEE; 2007:365â€“368. Sept. 16 2007Oct. 19 2007
Muller KR, Mika S, Ratsch G, et al.: An introduction to kernelbased learning algorithms[J]. Neural Networks, IEEE Transactions on 2001, 12(2):181â€“201. 10.1109/72.914517
Dambreville S, Rathi Y: NAMIC Publications Database on Kernel PCA for Segmentation. 2012. http://www.namic.org/Wiki/index.php/Projects:KPCASegmentation
Unnikrishnan R, Pantofaru C, Hebert M: A measure for objective evaluation of image segmentation algorithms[C]//Computer Vision and Pattern RecognitionWorkshops, 2005. San Diego, CA, USA: CVPR Workshops. IEEE Computer Society Conference on. IEEE; 2005:34â€“34. Date of Conference: 25â€“25 June 2005
MeilÄƒ M: Comparing clusterings by the variation of information[M]//Learning theory and kernel machines. Volume 2777. Springer Berlin Heidelberg; 2003:173â€“187.
Yang AY, Wright J, Ma Y, et al.: Unsupervised segmentation of natural images via lossy data compression[J]. Comput Vis Image Underst 2008, 110(2):212â€“225. 10.1016/j.cviu.2007.07.005
Acknowledgments
This work is supported in part by grants from National Natural Science Foundation of China for the Youth(NSFC:61103165), Guangdong Innovative Research Team Program (No. 2011S013) of China, Shenzhen Basic Research for Outstanding Youth(JC201005260248A), Shenzhen Basic Research Program(JCYJ20130401170306812) and National Natural Science Foundation of China (NSFC: 81171402). We thank Shibin Wu for helpful conversations.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare they have no competing interests.
Authorsâ€™ contributions
QL suggested the algorithm for images analyzing and processing, implemented it and analyzed the images. WJ gave the suggestion on algorithm analysis, experiment discussion and manuscript modification. TW, NG, SF and LL performed the acquisition of the abdominal MR images and manuscript discussion. JG and YQ expressed opinions on the evaluation metric of the segmentation results. All authors have read and approved the final manuscript.
Qing Luo, Wenjian Qin contributed equally to this work.
Authorsâ€™ original submitted files for images
Below are the links to the authorsâ€™ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Luo, Q., Qin, W., Wen, T. et al. Segmentation of abdomen MR images using kernel graph cuts with shape priors. BioMed Eng OnLine 12, 124 (2013). https://doi.org/10.1186/1475925X12124
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1475925X12124
Keywords
 Abdominal Organ
 Initial Contour
 Weak Boundary
 Shape Prior
 Shape Template