Spectral clustering for TRUS images
© Mohamed and Salama; licensee BioMed Central Ltd. 2007
Received: 13 October 2006
Accepted: 15 March 2007
Published: 15 March 2007
Identifying the location and the volume of the prostate is important for ultrasound-guided prostate brachytherapy. Prostate volume is also important for prostate cancer diagnosis. Manual outlining of the prostate border is able to determine the prostate volume accurately, however, it is time consuming and tedious. Therefore, a number of investigations have been devoted to designing algorithms that are suitable for segmenting the prostate boundary in ultrasound images. The most popular method is the deformable model (snakes), a method that involves designing an energy function and then optimizing this function. The snakes algorithm usually requires either an initial contour or some points on the prostate boundary to be estimated close enough to the original boundary which is considered a drawback to this powerful method.
The proposed spectral clustering segmentation algorithm is built on a totally different foundation that doesn't involve any function design or optimization. It also doesn't need any contour or any points on the boundary to be estimated. The proposed algorithm depends mainly on graph theory techniques.
Spectral clustering is used in this paper for both prostate gland segmentation from the background and internal gland segmentation. The obtained segmented images were compared to the expert radiologist segmented images. The proposed algorithm obtained excellent gland segmentation results with 93% average overlap areas. It is also able to internally segment the gland where the segmentation showed consistency with the cancerous regions identified by the expert radiologist.
The proposed spectral clustering segmentation algorithm obtained fast excellent estimates that can give rough prostate volume and location as well as internal gland segmentation without any user interaction.
Prostate disease has a foremost impact on the quality of life of elderly men. Benign enlargement of the prostate frequently causes bladder outlet difficulty [1, 2]. Malignant diseases of the prostate are considered a significant cause of death. Moreover, there has been an increase in the reported prevalence of prostate cancer. This is believed to be due to increased prostate awareness, prostate specific antigen (PSA) screening, and imaging techniques such as TransRectal UltraSound (TRUS).
Prostate Cancer is typically diagnosed by conduction a biopsy operation. Prostate Cancer diagnosis is usually aided by using TRUS images. This procedure is two manifold, first the prostate boundary is detected, Second the prostate tissue should be segmented and/or classified to different regions. The work done in this paper focuses on the TRUS image segmentation. The gland is first segmented from the background, second the outlined gland is further segmented into different regions.
The first segmentation process (prostate boundary segmentation from TRUS images) is crucial for some major applications in prostate disease decisions such as aiding diagnosis [3, 4] and treatment planning [5, 6]. Prostate Brachtherapy involves accurately insertion of radioactive materials (seeds) into the gland according to a predetermined plan  in which the gland location and boundary should be determined.
The purpose of the second segmentation process is mainly Computer Aided Diagnosis (CAD). The internal segmentation process involves segmenting the prostate gland into regions where different regions that represent different tissue textures in the gland are identified. The region segmentation highlights the different gland regions for either the purpose of feature analysis, or to augment the radiologist decision in highlighting the suspicious regions.
The internal segmentation system proposed in this paper can serve as a preprocessing stage for any CAD system. In the earlier work features were constructed for the whole image, where the whole image is divided into small squares and features are then constructed from each of those squares [7, 8]. This is considered a time consuming process. Therefore internally segmenting the prostate and studying only the highly suspicious regions is expected to be more accurate and efficient.
The proposed segmentation method has a different foundation than the previously proposed segmentation systems as it relies mainly on graph theory techniques. On the other hand, the older segmentation methods are mainly: edge base segmentation, texture based segmentation and model based segmentation. Each of these methods requires prior knowledge such as a seed point or three points on the boundary. In the proposed TRUS image segmentation method, spectral clustering treats the image as a weighted undirected graph and finds the segmentation result by obtaining the minimum cut of this weighted graph based on the graph theory methods without any prior knowledge to either a seed point or any point on the boundary.
The TRUS images used in this work are obtained from University of Western Ontario and are derived from Aloka 2000 ultrasound machine using a broadband 7 MHz linear transducer and a field of view of approximately 6 cm. A set of 29 radiologist identified TRUS images were used for this study.
2. Related work
The prostate boundaries are typically identified from TRUS images. Although manual outlining of the prostate border enables the prostate volume to be determined accurately , it is time consuming and tiresome. Moreover, since, in TRUS images' quality is not very good, therefore, traditional edge detectors are unable to extract the correct boundaries. Therefore, a number of investigations have been devoted to design automatic or semi-automatic methods that are suitable for segmenting the prostate boundary from ultrasound images.
A segmentation method that depends on clustering each pixel of an ultrasound image was introduced in . In this method, along with the relative position of the pixel, four energy measures were used to determine the cluster that the pixel belongs to. A drawback of this method is that the number of clusters is not predictable for a particular image; therefore, the prostate might be represented by disconnected regions.
Artificial Neural Networks (ANNs) was also introduced as a method for prostate segmentation in transrectal ultrasound images . This method segments images to prostate and non-prostate regions. Three neural network architectures have been proposed. This method needs lots of training data in order to train the ANN; moreover the training process is lengthy and computationally expensive.
Active contours were introduced in  and are used since then as one of the main methods for prostate boundary detection. The active contours model is used in  for prostate boundary detection where constraints were imposed on the model's deformation according to a predefined model shape. In this method, one-dimensional wavelet transform was applied on the radial function of both the prior model and the deformed model. While it was demonstrated that this method detected the prostate boundary accurately for typical gland shapes, the dependence of the statistically derived prior model has limited its ability to segment the prostate with atypical shape.
In an attempt to enhance the active contours method a cubic spline interpolation technique was used in  to identify an initial contour based on four user-defined points. Then, Lobregt's discrete dynamic contour (DDC) model  was used to refine the boundary. This method was shown to be effective if the initial contour was defined accurately, however, the result was less satisfactory for segmenting an irregular boundary that could not be accurately approximated by the initial contour, and further human intervention was required under this condition.
Another semi-automatic segmentation algorithm based on the dyadic wavelet transform and the discrete dynamic contours was used in . In this method first a spline interpolation is used to determine the initial contour based on four user-defined initial points. Then the discrete dynamic contour refines the initial contour based on the approximate coefficients and the wavelet coefficients generated using the dyadic wavelet transform. A selection rule is used as well to choose the best contour.
A common deformable model was also used in  to segment the prostate in transrectal ultrasound images. The new enhancement was the use of a Gabor filter bank in both multiple scales and multiple orientations to characterize the prostate boundaries. The Gabor features are then reconstructed to be invariant to the rotation of the ultrasound probe. Then, the segmentation is obtained by minimizing the energy function of the prostate shape model. The model focuses on the similarity of different Gabor features at different deformation stages using a multiresolution technique.
Another deformable models based research presented an approach where model initialization and constraining model evolution are based on prior knowledge about the prostate shape . The prostate shape has been modeled using deformable superellipses.
Deformable models were also used in  for automatic segmentation of trans-abdominal ultrasound images of the prostate. In this method a filter is used to enhance the contours without changing the information in the image. Adaptive morphological and median filtering are employed to detect the noise-containing regions and smooth these areas. Then a heuristic optimization algorithm begins to search for the contour initialized from a prostate model.
All the active contours based methods depend mainly on the initial points set by the user as well as the initial contour generation. Most of the research is focused on changing the number of initial points or changing the method used to obtain the initial contour. Optimizing the energy function is also an area of research using the active contour method.
Another segmentation technique to extract prostate contours from Transrectal Ultrasound images using Sticks filter to reduce the speckle was proposed in . Equi-spaced radii were then projected from an arbitrary seed point inside the prostate cavity towards its boundary. Candidate edge points are then obtained along each radius which include the edge points and some false points. This approach is dependent on the choice of the seed point which might mislead the prostate contour extraction method.
The sticks method was used in another algorithm for prostate boundary detection in . The algorithm provides prostate edge detection as a visual guide to the observer using edge delineation. It is then followed by manual editing. This edge detection algorithm contains three stages. First, the sticks algorithm is used to enhance contrast and reduce speckle in the image. Second, the resulting image is smoothed using an anisotropic diffusion filter. Finally, some basic prior knowledge of the prostate, such as shape and echo pattern, is used to detect the most probable edges which indicate the prostate shape. In the last stage, the information is then integrated by using a manual linking procedure on the detected edges. The drawback of this method is that it depends on prior knowledge of some prostate features, a limitation that makes it limited to typical prostate gland shapes and echo patterns. Moreover it needs manual editing at the final stage to obtain the prostate boundaries.
Most of the above-mentioned methods depended on statistical estimation for initialization and some of these initialization methods were not accurate enough . Some other methods depend on choosing the right seed point, otherwise the algorithm will not converge to the right boundary. While all the above methods depend on solving optimization problems that are parameter sensitive and time consuming.
Generally, prostate segmentation methods have limitations when the image contains shadows with similar gray level and texture attached to the prostate. In these cases the segmentation error may increase considerably. Another problem may be the lack of sufficient number of training images if a learning technique is used. Algorithms based on active contours have been quite successfully implemented with the major drawback that they depend on user interaction to determine the seed points or the initial snake.
Based on the previous literature review of the existing methods, a new approach should ideally be:
Independent on user interaction as user interaction (e.g. defining seed points or initial contours) has drawbacks such as time consumption, human error etc.
Independent on training images where training images is typically difficult to obtain, especially if the samples should be prepared by an expert. Hence, sample-based learning should be avoided.
Independent on noise level where the approach must be robust with respect to the presence of noise and shadow.
The proposed segmentation algorithm in this paper has a totally different basis as it doesn't depend on the inherited snakes' algorithm. This will get rid of designing the energy function, optimizing it and accurately selecting the seed points or initial contour points.
3. Spectral clustering
The Human Visual System (HVS) can effectively identify objects in a scene and can often segment the scene to coherent segments or clusters. There has been a tremendous amount of research done to achieve the same level of performance obtained by the HVS. Various methods have been introduced in literature to segment ultrasound images such as Kalman filters  and statistical shape models in which the prostate is segmented from the background in order to determine the gland volume. However most of these methods require large amount of human interaction. The Gabor filters were introduced in this field as a method to internally segment the prostate gland that was already manually segmented from the background  in which a Gabor filter was designed to automatically and accurately segment the TRUS images.
Spectral Clustering methods are applied in this paper in order to give a fast segmentation results that don't require either filter design or any human interaction. Spectral Clustering have been introduced for data clustering and was applied in different fields. The spectral clustering usually represents the data by a weighted graph and the eigenvectors of the affinity (or similarity) matrix of this graph are used for the segmentation . In the problem of image segmentation the image pixels are considered as the data points as shown in .
3.1. Graph based image segmentation
Given an image I, a graph G = (V, E, W) is constructed with the pixels represented by the graph nodes V, and the pixels within a distance V ≤ G r are connected by a graph edge E. The weight W(i, j) measures the likelihood of pixel i and j being in the same cluster. Partitioning of this graph represents the image segmentation [24–26].
Assigning weights to graph edges
The pair-wise pixel affinity graph determines the segmentation accuracy. Therefore, as recommended in  two simple local grouping cues are used which are the intensity and contours.
1. Grey Level Intensity: neighboring pixels with close intensity are most likely to be in the same region.
Where X i and I i represent pixel location and intensity.
Connecting pixels considering only intensity and location usually gives bad segmentation due to the texture that is present in the TRUS images. Therefore the principal image contours (edges) are also considered for the segmentation of TRUS images.
2. Dominant Contours: the image edges are considered useful when the neighboring regions have the same cutter. The affinity between two pixels is calculated by measuring the image edges between them.
Where line(i, j) is a straight line joining pixels i and j and Edge(x) is the edge strength at location(x).
The two cues are combined in this work in the form:
Spectral clustering segmentation algorithms
In , a clustering algorithm based on thresholding the largest eigenvector of the affinity matrix was suggested. While in  the authors have argued for using a totally different eigenvector for solving these types of segmentation problems. Rather than examining the first eigenvector of W they examined the generalized eigenvectors. Let D be the degree matrix of W:
The generalized eigenvector y i is a solution to:
(D — W)y i = λ i Dy i (5)
Solving the generalized eigenvector minimizes the Ncut which in turn produce the optimum segmentation as proved in . In this case the generalized eigenvector corresponding to the second smallest eigenvalue was used. Thresholding this second eigenvector to obtain the segmentation result was suggested. This method is adopted for the application of TRUS image segmentation and it yields to a segmentation that minimizes the normalized cut:
where: A ∪ B = V and A ∩ B = 0, and
Therefore the solution to the segmentation problem minimizes the affinity between groups normalized by the affinity within the same group. In this work, the spectral clustering is used for the first time for TRUS image segmentation using the approach proposed in .
4. Proposed lgorithm implementation
The edge map of the TRUS image is obtained using Canny edge detection method
The Affinity matrix is created using equation 3.
The eigenvectors are calculated and reshaped to be shown in the figures.
Eigen vector Discretization: In the best case scenario the second smallest eigenvector should take on two discrete values and the signs of the values can tell how to partition the graph. The second smallest eigenvector obtained in our case is a continuous vector; therefore it needs to be discretized in order to find a splitting point to obtain the clustering. In this work the splitting point that minimizes Ncut that is shown in equation (6) is chosen.
5. Validation methods
Distance δ= Average Euclidean distance (in pixels) between the algorithm-based segmentation and the manual segmentation. For each pixel the distance is defined as the shortest Euclidean distance between that pixel and the pixels located on the other contour.
Area Difference AD : where S Manual is the area of the manual segmentation and S Algorithm is the area of the algorithm-based segmentation.
Area Overlap AO:
Since the main purpose of this paper is to introduce the new method and prove its concepts, one validation method is used which is the AO. The AO can be considered as a good representation for the algorithm success in segmenting the prostate as it measures the area of the gland that the algorithm could capture. The authors realize that more images are needed for the investigation to be generalized.
6. Experimental results
In this section, some results that show the correlation between the desired segmentation and the eigenvectors of the affinity matrix corresponding to TRUS images already segmented (either from the background or into cancerous and non cancerous regions) by an expert radiologist.
6.1. Segmenting the prostate from the background
6.2. Internal regions segmentation
The spectral clustering image segmentation algorithm is mainly used in this work for internal gland segmentation. However it was applied in the previous subsection for the prostate gland segmentation from the TRUS image and it obtained high segmentation accuracy which proves its capability of dealing with TRUS images effectively. Therefore, it is used in this sub-section for ROI (Region of Interest) segmentation from the manually segmented prostate.
The proposed algorithm in this work is faster than the Gabor multi-resolution analysis that was used earlier for prostate ROI identification  on the expense of giving a rough estimate of the internal segmentation than the earlier presented work. The algorithm proposed in this paper can then be used for suspicious regions estimation in an online application. It can also be used to support the decisions obtained using the feature analysis methods especially if the later contradicts the radiologist's decision.
Typical prostate shape
Unusual prostate shape
From the above results it is clear that the spectral clustering can be used as a preliminary estimate for the cancerous regions during the imaging procedure. Moreover it can be used as a support for the decision making either by the radiologist or by the other CAD systems. It can also estimate the prostate volume, location and size from the TRUS image effectively. The second smallest eigenvector proved to be consistent with the radiologist manual segmentation for the two different segmentation problems tackled in this paper
In this paper a novel technique is proposed to segment the prostate using TRUS images. The new strategy that is introduced in this work is based on the spectral clustering algorithm. Spectral clustering has the benefit of being built on a totally different foundation that doesn't include any contour or seed point estimation. The proposed spectral clustering segmentation method is inspired from the graph theory techniques. The idea of the spectral clustering depends mainly on treating the image as a weighted graph and searched for the minimum cut of that graph. This idea has a totally different prospective than the well known prostate segmentation methods such as the snakes. The proposed spectral clustering segmentation method is accurate, simple to implement and doesn't involve any energy functions to be built or optimized. From the results and analysis shown in this paper, it can be concluded that spectral clustering can be considered as a new advance for prostate segmentation that proved its ability to accurately segment the gland from the background for typical as well as atypical prostate shapes. Moreover, it is also clear from the results obtained in this work that the spectral clustering method is able to roughly segment the cancerous regions that proved consistency with the regions identified by the doctor. The proposed method is able to recognize regions regardless of the prostate shape and the spatial location of the cancer within the gland.
In conclusion, the algorithm obtained fast excellent estimates that can give rough prostate volume as well as cancerous regions segmentation which can be used for online application.
- Lee C, Kozlowski JM, Grayhack JT: "Etiology of benign prostatic hyperplasia". Urol Clin North Amer 1995, 22: 237–246.Google Scholar
- Barry MJ: "The epidemiology and natural history of benign prostatic hyperplasia". Curr Opin Urol 1994, 4: 3–6. 10.1097/00042307-199401000-00002View ArticleGoogle Scholar
- Catalona WJ, Beiser JA, Smith DS: "Serum free prostate specific antigen and prostate specific antigen density measurements for predicting cancer in men with prior negative prostatic biopsies". J Urol 2167, 158: 2162–1997. 10.1016/S0022-5347(01)68187-4View ArticleGoogle Scholar
- Presti JC Jr, Hovey R, Bhargava V, Carroll PR, Shinohara K, Moul JW: "Prospective evaluation of prostate specific antigen and prostate specific antigen density in detection of carcinoma of prostate: ethnic variations". J Urol 1997, 157: 907–912. 10.1016/S0022-5347(01)65080-8View ArticleGoogle Scholar
- Grimm PD, Blasko JC, Ragde H, Sylvester J, Clarke D: "Does brachytherapy have a role in the treatment of prostate cancer". Hemat Oncol Clin North Amer 1996, 10: 653–673. 10.1016/S0889-8588(05)70359-2View ArticleGoogle Scholar
- Hill CR, ter Haar GR: "High intensity focused ultrasound potential for cancer treatment". Br J Radiol 1303, 68: 1296–1995.View ArticleGoogle Scholar
- Bassat O, Sun Z, Mestas JL, Gimenez G: "Texture Analysis Of Ultrasound Images of the Prostate By Means of Co-occurrence Matrices". Ultrasonic Imaging 1993, 15: 218–237. 10.1006/uimg.1993.1014View ArticleGoogle Scholar
- Ulrich Scheipers, Helmut Ermert, Hans-Joerg Sommerfeld, Miguel Garcia-Schürmann: "Ultrasonic multifeature tissue characterization for prostate diagnostics". Ultrasound in Medicine and Biology 2003,29(8):1137–1149. 10.1016/S0301-5629(03)00062-0View ArticleGoogle Scholar
- Tong S, Cardinal HN, Downey DB, Fenster A: Analysis of linear, area, and volume distortion in 3D ultrasound imaging. Ultrasound Med Biol 1998, 24: 355–73. 10.1016/S0301-5629(97)00268-8View ArticleGoogle Scholar
- Richard WD, Keen CG: Automated texture-based segmentation of ultrasound images of the prostate. Computerized Medical Imaging and Graphics 1996,20(3):131–140. 10.1016/0895-6111(96)00048-1View ArticleGoogle Scholar
- Prater JS, Richard WD: Segmenting ultrasound images of the prostrate using neural networks. Ultrasound Imaging 1992, 14: 159–185. 10.1016/0161-7346(92)90005-GView ArticleGoogle Scholar
- Kass M, Witkin A, Terzopoulos D: "Snakes: active contour models". Int J Comput Vision 1987, 1: 321–331. 10.1007/BF00133570View ArticleGoogle Scholar
- Knoll C, Alcaniz M, Grau V, Monserrat C, Juan MC: "Outlining of the prostate using snakes with shape restrictions based on the wavelet transform". Pattern Recognition 1999, 32: 1767–1781. 10.1016/S0031-3203(98)00177-0View ArticleGoogle Scholar
- Ladak HM, Mao F, Wang Y, Downey DB, Steinman DA, Fenster A: Prostate boundary segmentation from 2D ultrasound images. Medical Physics 2000, 27: 1777–1788. 10.1118/1.1286722View ArticleGoogle Scholar
- Lobregt S, Viergever M: A discrete dynamic contour model. IEEE Trans Med Imaging 1995, 14: 12–24. 10.1109/42.370398View ArticleGoogle Scholar
- Bernard Chiu, George H Freeman, Salama MMA, Aaron Fenster: "Prostate segmentation algorithm using dyadic wavelet transform and discrete dynamic contour". Phys Med Biol 49(21):4943–4960. 7 November 2004Google Scholar
- Shen D, Zhan Y, Davatzikos C: "Segmentation of the prostate boundaries from ultrasound images using statistical shape model". IEEE Transactions on Medical Imaging 22(2003):539–551.Google Scholar
- Gong L, Pathak SD, Haynor DR, Cho PS, Kim Y: Parametric Shape Modeling Using Deformable Superellipses for Prostate Segmentation. IEEE Transactions on Medical Imaging 2004,23(3):340–349. 10.1109/TMI.2004.824237View ArticleGoogle Scholar
- Betrounia N, Vermandela M, Pasquierc D, Maoucheb S, Rousseaua J: Segmentation of abdominal ultrasound images of the prostate using a priori information and an adapted noise filter. Computerized Medical Imaging and Graphics 2005, 29: 43–51. 10.1016/j.compmedimag.2004.07.007View ArticleGoogle Scholar
- Pathak SD, Chalana V, Haynor DR, Kim Y: Edge-guided boundary delineation in prostate ultrasound images. IEEE Transactions on Medical Imaging 2000, 19: 1211–1219. 10.1109/42.897813View ArticleGoogle Scholar
- Abolmaesumi P, Sirouspour M: "An interacting multiple model probabilistic data association filter for cavity boundary extraction from ultrasound images.". IEEE Transactions on Medical Imaging 23(2004):772–784.Google Scholar
- Mohamed SS, Salama MMA, Kamel M, El-Sadaany EF, Rizkalla K, Chin J: "Prostate cancer multi-feature analysis using TRUS images". Physics in Medicine and Biology 50(15):N175-N185. 7 August 2005 10.1088/0031-9155/50/15/N02Google Scholar
- Ng AY, Jordan M, Weiss Y: "On spectral clustering: Analysis and an algorithm". In Advances in Neural Information Processing Systems 14. MIT Press; 2002:849–856.Google Scholar
- Shi J, Malik J: "Normalized cuts and image segmentation". IEEE Transactions on Pattern Analysis and Machine Intelligence 22(2000):888–905.Google Scholar
- Timothee Cour, Florence Benezit, Jianbo Shi: "Spectral segmentation with multiscale graph decomposition". IEEE International Conference on Computer Vision and Pattern Recognition (CVPR) 2005.Google Scholar
- Perona P, Freeman WT: "A factorization approach to grouping". In Proc ECCV Edited by: Burkardt H, Neumann B. 1998, 655–670.Google Scholar
- Shi J, Malik J: "Normalized cuts and image segmentation". IEEE Conf. on Computer Vision and Pattern Recognition 1997, 731–737.Google Scholar
- Archip N, Rohling R, Cooperberg P, Tahmasebpour H, Warfield SK: "Spectral Clustering Algorithms for Ultrasound Image Segmentation". MICCAI, LNCS 3750 2005, 862–869.Google Scholar
- Sahba Farhang, Tizhoosh HamidR, Salama MagdyM: A coarse-to-fine approach to prostate boundary segmentation in ultrasound images. BioMedical Engineering OnLine 2005,4(4):58. 10.1186/1475-925X-4-58View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.