An Iterative Boundary Random Walks Algorithm for Interactive Image Segmentation
Abstract
The interactive image segmentation algorithm can provide an intelligent ways to understand the intention of user input. Many interactive methods have the problem of that ask for large number of user input. To efficient produce intuitive segmentation under limited user input is important for industrial application. In this paper, we reveal a positive feedback system on image segmentation to show the pixels of selflearning. Two approaches, iterative random walks and boundary random walks, are proposed for segmentation potential, which is the key step in feedback system. Experiment results on image segmentation indicates that proposed algorithms can obtain more efficient input to random walks. And higher segmentation performance can be obtained by applying the iterative boundary random walks algorithm.
I Introduction
Interactive image segmentation with simple user input has been an important research topic in the field of image analysis recently [1, 2, 3, 4, 5]. There are some successful applications, such as video surveillance [6], image retrieval [7], medical image analysis [8], and object detection [9], etc. To segment a foreground object out from its surrounding background, the user should give a certain user inputs for interactive algorithms. A large number of user inputs are expensive and tedious to acquire, especially in processing video sequences. This paper attempts to interactively cut out a desired object with less user inputs. Our goal is to develop iterative segmentation algorithm that acquire more reliable input, i.e. semiseeds, according to segmentation results.
Ia Related Work
Generally, interactive image segmentation algorithm produces a binary map, i.e., a pixel belongs to either foreground or background. The approaches for interactive segmentation can roughly be classified into two categories according to the type of user inputs, and they are briefly reviewed as follows.
1) The user input is an area that close to the desired cutting contour. Active contour [10] and intelligent scissor [11] are two earliest interactive image segmentation methods. They both utilize the boundary properties of image. The active contour algorithm requires to place input near the desired boundary, and intelligent scissor algorithm requires to place input along the desired contour of the foreground object. The intelligent scissor algorithm is difficult to be used in processing images with low contrast or noisy boundaries. The problem of active contour algorithm is likely to be involved in local minimum. Recently, graph cut based methods [12], like LazySnapping algorithm [13] and GrabCut method [14], were developed by utilizing the boundary and regional properties. The LazySnapping algorithm models an image as a graph where each node represents a pixel and two nodes are connected with a weighted edge. The GrabCut method models the foreground and background pixels according to the Gaussian Mixture Models. However, GrabCut is difficult to process image whose foreground and background share similar color distribution. Another popular approach, SIOX algorithm [15], is derived from color signature. It works well with noise and videos, but depends heavily upon the foreground and background color distribution.
2) The user input is the labeling of some pixels which belong to background or foreground. Random walks algorithm [16] efficiently segment the image with only two types of user input, background seeds and foreground seeds. After setting the background seeds and foreground seeds, each pixel in image can be classified by the probability that a random walker starting from one pixel first arrives at one of the foreground seeds. The performance of Random walks algorithm is sensitive to the positions and number of the seeds. Constraints random walks [1] adds two more types of user inputs to reflect the user’s intention. By solving the constraints harmonic functions, it can achieve better image segmentation performance than the random walks algorithm. However, it cannot handle transparent or semitransparent boundaries such as semilucent hair.
IB Our Work
In fact, many efficient interactive image segmentation algorithms can provide intelligent ways to understand the intention of user input. However, less of that considers image segmentation under less user input. To efficient produce intuitive segmentation under less input, this paper reveals the positive feedback system on image segmentation to show the pixels of selflearning. With the positive feedback system, image is initially segmented by basic random walks algorithm using the original user inputs. After that, the segmentation results are subsequently used with the original user inputs to resegment the image iteratively. This process is repeated until convergence is reached.
In this paper, we proposed a segmentation potential to learn the connection of user input and segmentation results. The segmentation potential is close related with the misclassified possibility. It can indicate that how to select the segmentation results to enlarge the input. Two methods, iterative random walks and boundary random walks, are developed for segmentation potential. Iterative random walks can add more background and foreground seeds based on segmentation results. Boundary random walks can directly modify the probability value of boundary seed. Lastly, the iterative boundary random walks algorithm is proposed by combined iterative random walks and boundary random walks.
The contributions of this paper are summarized as follows.

A positive feedback system is proposed for image segmentation. The pixels of selflearning can be achieved within feedback system.

The proposed algorithms can alleviate limited user input problems by utilizing useful information from the segmentation results.

Two segmentation algorithms, boundary random walks and iterative random walks, are proposed to work for segmentation potential. Higher segmentation results can be efficiently obtained for the proposed algorithms as verified on image segmentation.
The remainder of this paper is organized as follows. Section II reviews the basic random walks algorithm and reveals the positive feedback system on image segmentation. Section III details the boundary random walks and iterative random walks. Some experiment results obtained by the proposed algorithms are provided in Section IV and conclusion are given in Section V.
Ii Brief Review on Random Walks Algorithm for Image Segementation
To make this paper be selfcontained, we briefly review the random walk algorithm for image segmentation in this section. An image can be modeled as a graph as shown in Fig. 1. Each node of the graph represents a pixel and only the neighboring nodes are connected with undirected edges shown in Fig. 1(a). Let denotes a set of vertices and denotes a set of edges bounded by vertices and . The graph can be represented by . The weight of edge is defined as , and the degree of node is defined as . In applications of image segmentation, the edge weight could be defined as
(1) 
where is a scaling factor, and are the gray values corresponding to vertices and respectively.
As illustrated in Fig. 1(b), the input to random walks algorithm [16] are the marked foreground seeds and marked background seed . By defining of each vertice as the probability that a random walker starts from vertices and arrives at the first time before reaching , we have for all the background seeds and for all the foreground seeds . The probability for the unmarked vertices can be calculated by solving the following optimization problem
(2) 
Differentiating the objective function of (2) to and setting the derivative to zero, we have
(3) 
In practice, the harmonic function (3) is difficult to be solved. An important work [16] is to transform the above harmonic equation to be a linear equation
(4) 
where the vector represents the probabilities of the marked seeds (background seeds and foreground seeds), and is the probability vector of unmarked vertices. and are submatrix of the Laplacian matrix L
(5) 
where the element of L is defined as
(6) 
There are many efficient methods available for solving the linear equation (4). After solving (4), the foreground object is segmented as the set of pixels whose probability are greater than 0.5, as shown in Fig. 1(b)(c). In practice, the probability of a vertex is more close to 0.5, the corresponding pixel is more possible to be misclassified. Fig. 1(c) indicate that the pixels whose probability within a range almost located at the boundary of segmentated object. These pixels, named boundary seeds , contain many misclassified pixels. It is possible to exploit the information contained in to improve image segmentation performance.
Iii Proposed Methods for Image Segmentation
After analyzing the conventional random walk algorithm for image segmentation, we firstly propose a new index, called segmentation potential, to reflect the performance of segmentation. Two new methods, iterative random walks and boundary random walks, are proposed to enhance the segmentation performance by reducing the segmentation potential. By combining the advantages of the above two new methods, an iterative boundary random walks algorithm is also proposed to further improve segmentation performance.
Iiia Segmentation Potential
The random walks method can be regarded as a semisupervised learning problem [17], where background and foreground seeds are defined as the labeled samples. Despite the unlabeled samples are used to learn the classification model, a large number of sufficient labeled samples are expensive to obtain. In the case of limited training samples, if we know the samples which have been correctly classified, these samples could be used to enlarge the set of training samples. The retrained classifier is possible to be more precise [18].
Since usually we have no knowledge on the true labels of unlabel samples, it is important to design a scheme with which we can integrate the information on segmentation results, background/foreground seeds, boundary seeds into segmentation process. A potential way is to involve the concept of feedback system as shown in Fig. 2. Before introducing the feedback system in Fig. 2, we first propose a new index, called segmentation potential, to evaluate the stability of segmentation. High stability means that the image could be reliably segmented and the probabilities of pixels are close to or . Hence, we define the segmentation potential as
(7) 
where is probability of pixel after random walk processing. It is clear that the lower the segmentation potential , the more reliable the image to be correctly segmented. In this paper, we design new scheme on segmentation aiming to reduce the segmentation potential in processing.
Turning back to the scheme in Fig. 2, the input of more foreground/background seeds is able to reduce the segmentation potential and the more boundary seeds found in processing may cause the higher value of segmentation potential. High value of segmentation potential means potential poor performance on image segmentation. Using the segmentation results in further processing is able to enlarge the set of foreground/background seeds. From the relationship presented in Fig. 2, we are able to devise some new iterative algorithms to improve the performance of segmentation.
IiiB Iterative Random Walks
The pixels with probability close to or could be reliably classified as foreground/background seeds. In the first proposed iterative method, we select the segmented pixel whose probability value far away as new background/foreground seed
(8) 
where is a threshold. Since there are a large number of pixels could be selected as new background/foreground seeds, a part of these seeds based on (8) could been randomly selected and merged into the input sets and .
With the selected new foreground/background seeds, the basic random walks is modified to repeat iteratively. The details of proposed iterative random walk (IRW) algorithm is presented in Algorithm 1. The segmented result of the previous iteration could be used to enlarge the input marked sets of the next iteration. After some iterations, a large number of new background/foreground seeds can be obtained. In Fig. 3(a), the illustration of iterative random walks algorithm in 3 iterations is shown. With the autoselected background/foreground seeds involved in processing, more stable segmentation results with less number of boundary seeds and low segmentation potential could be obtained.
IiiC Boundary Random Walks
The segmentation potential plays an important role in segmentation processing as shown in Fig. 2. In practice, the higher the segmentation potential, the higher probability that pixels misclassified. The IRW algorithm reduces segmentation potential by adding more background and foreground seeds into processing. In this Section, we propose another way to reduce the segmentation potential by directly removing the boundary seeds using a modified random walk algorithm.
As we have discussed, the boundary pixels have probabilities near 0.5. These pixels also cause the increasing of segmentation potential. In order to reduce the number of boundary pixels, we modify the random walk algorithm by introducing another optimization object, i.e., the reducing of segmentation potential. Hence, the new algorithm, called boudary random walk algorithm, is formualted as the following multiple objective optimization probem
(9) 
For easy processing, the above optimization problem could be reformulated as
(10) 
where is a tradeoff factor whose range is , to guarantee the convexity of (10).
Differentiating the object function (10) with respect to each for and setting the derivative equal to zero, we obtain the following equations
(11) 
Similar to (4), the above equations can be transformed to the following linear equations
(12) 
where the vector e is a binary vector whose th element is defined as
(13) 
Compared with (4), the difference is that the new Laplacian matrix in (12) is constructed by . , are submatrices of defined as (4).
IiiD Iterative Boundary Random Walks Algorithm
Both the iterative random walks and boundary random walks algorithms can efficiently reduce segmentation potential. Hence, better performance on segmentation is achieved. To further improve segmentation performance, we can use the boundary random walks approach to replace the basis random walks approach in IRW algorithm and get a new method called iterative boundary random walks (IBRW) algorithm (The details of IBRW is given in Algorithm 2). In Fig. 3(b), the illustration of IBRW is presented. Compared with IRW in Fig. 3(a), the IBRW could converges quickly and reduces segmentation potential more efficiently than IRW.
Iv Numerical Results and Discussions
In order to evaluate the performance of proposed methods, in this section, the proposed methods are tested on some image datasets, including the Berkeley dataset [19] and the MSRC dataset [20].
Iva Experimental Study for IRW Algorithm
We first verify the performance improvement by adding more background and foreground seeds into random walks algorithm in each iteration, i.e., the IRW algorithm. The key merit of IRW is to mitigate the limited training sample problem by utilizing the previously segmented results. Hence, it is important to obtain reliable background and foreground seeds (semiseeds) based on probability distribution to improve segmentation results iteratively
Fig. 4 shows the image segmentation results of IRW algorithm in 3 iterations. The first column of Fig. 4 is the results of basic random walks, i.e., the initialization of IRW algorithm. The next 3 rows are the results of IRW. The IRW automatically select the semiseeds by (8). From the second row of Fig. 4, i.e., the probability maps of IRW in each iteration, it can see that the pixels near the left shoulder have higher probability values after the semiseeds are involved into the IRW. This is the positive feedback effect of semiseeds, i.e., we can find more semiseeds from probability maps of lower segmentation potential, and in return, more semiseeds lead to probability distribution of lower segmentation potential. Compared with the results in Fig. 4(a)(d), after adding 70 semiseeds into random walk algorithm, the number of misclassified pixels is reduced from 623 to 584. The boundary pixels whose probability values close to 0.5 are reduced from 205 to 141. More details of IRW in segmentation process are shown in the third and fourth rows of Fig. 4. We also provide the zoomin segmentation results (left shoulder). From the zoomin results, it is clear that the misclassified pixels are iteratively corrected by adding more semiseeds. Those experimental results verify that the introducing of semseeds into random walk algorithm could significantly improve its performance.
IvB Experimental Study for IBRW Algorithm
We also carried out some experiments to evaluate the performance of BRW and IBRW algorithms. In Fig. 5, the BRW is compared with the basic random walk algorithm. To conduct a fair comparison, same background and foreground seeds are used to initiate the random walk and BRW algorithm. From Fig. 5(b), it is clear that the number of boundary seeds is much smaller that that of the random walk algorithm. The probability map of BRW also has lower segmentation potential. These experimental results reveal the effectiveness of the proposed BRW algorithm.
In Fig. 6, more results are presented to compare the performance of IBRW and RW. In the initialization of both algorithms, the same input foreground/background seeds are used as shown in Fig. 6(a). After several iterations, as shown in Fig. 6(d), the sets of foreground/background seeds for IBRW are enlarged with adding of semiseeds which are automatically selected by IBRW. It can be seen that the segmentation results of IBRW outperforms the random walks algorithm. Moreover, from the probability maps of IBRW and RW (Fig. 6(b) and Fig. 6(e)), it is clear that the segmentation potential has been efficiently reduced.
In Fig. 7, the performance comparison between IRW and IBRW is also illustrated in terms of boundary pixels and error pixels. From the results shown in Fig. 7, it is clear that both IRW and IBRW have decreased number of boundary and error pixels in each iteration. Since the IBRW uses the modified boundary random walk algorithm (10), IBRW has more faster convergence speed than IRW. From Fig. 7, we can also find that the IBRW can converge in about 3 iterations generally.
We further compare the performance of IBRW with some stateofart algorithms including GrabCut, LazySnapping, constrained random walks on 3 different images. The input of GrabCut algorithm is a rectangle covering the target object as shown in the top row of Fig. 8. For the other algorithms, the same background/foreground strokes are used as inputs. Compared with the other popular image segmentation methods, we can find that proposed IBRW achieves better segmentation performance.
In the following experiment, 50 images form MSRC dataset are tested by the proposed IBRW algorithm and more stateofart algorithms, like LazySnapping [13], GrabCut [14], SIOX [15], Random walks [16] and Constrained random walks [1]. All of these algorithms use the exactly same inputs of background/foreground strokes. Table I shows the error rates of all algorithms. The error rate is defined as the ratio between number of misclassified pixels and number of unclassified pixels, where misclassified pixels exclude those from the unclassified region [20]. For the other stateofart algorithms, we directly quote the best results from the literature [1], [13][16]. It can be seen that proposed algorithm outperforms the other stateofart algorithms and achieves very low error rates. Meaning that you have different input for these experiments.
IvC Sensitivity Analysis of Parameters
There are three parameters used in the proposed algorithm: the threshold for selection of semiseeds, the parameter for selection of boundary seeds and the tradeoff factor for BRW. It is important to study the sensitivity of proposed IBRW to these parameters. In this Section, we conduct some experiments on investigating the sensitivity of the proposed methods on Berkeley dataset.
The threshold plays a critical role in selection of background/foreground semiseeds. In the first experiment, the parameter is selected as discrete candidates in range . As shown in Fig. 9(a), we can find that IRW has low error rate when the threshold in the range of . However, if the threshold is less than 0.4, the error rate is less sensitive to the parameter. Hence, we suggest to use = 0.1 in experiments. If the threshold is too small, the selection of new seeds will involve in local field. If the threshold is too big, the selection of new seeds will unreliable.
The tradeoff factor controls the compromise between the weighted difference of probabilities and the segmentation potential. We expect that the probabilities of boundary seeds could be far away 1/2. Hence, it should use larger . However, to guarantee the problem (10) to be a convex optimization problem, the factor is limited in the range . Fig. 9 (b) shows that the IBRW algorithm performance is deteriorate when the factor is out of range .
Lastly, we evaluate the sensitivity of IBRW versus different combinations of and . Generally, large value of is encouraged to emphasize the selection of boundary seeds. The IBRW algorithm achieves bad performance for large ¦Ä values, since many nonboundary pixels are regarded as boundary seeds. As shown in Fig. 9 (c), there are many parameter combination near the optimal performance on area of small ¦Ë and ¦Ä. We fixed as (0.005, 0.1) in our experiments.
V Conclusion
The interactive image segmentation algorithm under limited user input play an important role in industrial application. It can avoid the expensive and tedious user input. To let machine more intelligent understand the intention of limited user input, the pixels of selflearning must present on the image segmentation. In this study, the most important work is to reveal positive feedback system on image segmentation. It provides a way to extend the basic random walks algorithm. We use the iterative boundary random walks algorithms, which it is combined iterative random walks with boundary random walks approach, to reduce segmentation potential. Experiment results on image segmentation shows that proposed algorithms can obtain more efficient input. And higher segmentation performance can be obtained by applying the IBRW algorithm. The proposed methods can also be applied to many other pattern recognitions field, e.g. limited training sample. Future work could focus on the selection of semiseeds, instead of the threshold method in IBRW.
References
 [1] W. Yang, J. Cai, J. Zheng, and J. Luo, “Userfriendly interactive image segmentation through unified combinatorial user inputs,” IEEE Trans. Image Process., vol. 19, no. 9, pp. 2470–2479, 2010.
 [2] L. Garcia Ugarriza, E. Saber, S. R. Vantaram, V. Amuso, M. Shaw, and R. Bhaskar, “Automatic image segmentation by dynamic region growth and multiresolution merging,” IEEE Trans. Image Process., vol. 18, no. 10, pp. 2275–2288, 2009.
 [3] J. Sourati, D. Erdogmus, J. Dy, and D. Brooks, “Accelerated learningbased interactive image segmentation using pairwise constraints,” IEEE Trans. Image Process., vol. 23, no. 7, pp. 3057–3070, July 2014.
 [4] B. Ham, D. Min, and K. Sohn, “A generalized random walk with restart and its application in depth upsampling and interactive segmentation,” IEEE Trans. Image Process., vol. 22, no. 7, pp. 2574–2588, July 2013.
 [5] T. N. A. Nguyen, J. Cai, J. Zhang, and J. Zheng, “Robust interactive image segmentation using convex active contours,” IEEE Trans. Image Process., vol. 21, no. 8, pp. 3734–3743, Aug 2012.
 [6] A. Meghdadi and P. Irani, “Interactive exploration of surveillance video through action shot summarization and trajectory visualization,” IEEE Trans. Vis. Comput. Gr., vol. 19, no. 12, pp. 2119–2128, Dec 2013.
 [7] J. Bach, S. Paul, and R. Jain, “A visual information management system for the interactive retrieval of faces,” IEEE Trans. Knowl. Data En., vol. 5, no. 4, pp. 619–628, Aug 1993.
 [8] D. L. Pham, C. Xu, and J. L. Prince, “Current methods in medical image segmentation1,” Annu. Rev. Biomed. Eng., vol. 2, no. 1, pp. 315–337, 2000.
 [9] J. Ye and G. Xu, “Geometric flow approach for regionbased image segmentation,” IEEE Trans. Image Process., vol. 21, no. 12, pp. 4735–4745, Dec 2012.
 [10] M. Kass, A. Witkin, and D. Terzopoulos, “Snakes: Active contour models,” Int. J. Comput. Vis., vol. 1, no. 4, pp. 321–331, 1988.
 [11] E. N. Mortensen and W. A. Barrett, “Interactive segmentation with intelligent scissors,” Graph. Models Image process., vol. 60, no. 5, pp. 349–384, 1998.
 [12] Y. Boykov and M.P. Jolly, “Interactive graph cuts for optimal boundary amp; region segmentation of objects in nd images,” in Proc. IEEE Int. Conf. Computer Vision, vol. 1, 2001, pp. 105–112.
 [13] Y. Li, J. Sun, C.K. Tang, and H.Y. Shum, “Lazy snapping,” ACM Siggraph, vol. 23, no. 3, pp. 303–308, Aug. 2004.
 [14] C. Rother, V. Kolmogorov, and A. Blake, “”grabcut¡± – interactive foreground extraction using iterated graph cuts,” ACM Siggraph, vol. 23, pp. 309–314, 2004.
 [15] G. Friedland, K. Jantz, and R. Rojas, “Siox: simple interactive object extraction in still images,” in Proc. IEEE Int. Symp. Multimedia, Dec 2005, p. 7.
 [16] L. Grady, “Random walks for image segmentation,” IEEE Trans. Pattern Anal. Mach. Intell., vol. 28, no. 11, pp. 1768–1783, 2006.
 [17] O. Chapelle, B. Schlkopf, and A. Zien, SemiSupervised Learning, 1st ed. The MIT Press, 2010.
 [18] Q. Jackson and D. Landgrebe, “An adaptive classifier design for highdimensional data analysis with a limited training data set,” IEEE Trans. Geosci. Remote. Sens., vol. 39, no. 12, pp. 2664–2679.
 [19] D. Martin, C. Fowlkes, D. Tal, and J. Malik, “A database of human segmented natural images and its application to evaluating segmentation algorithms and measuring ecological statistics,” in Proc. IEEE Int. Conf. Computer Vision, vol. 2, July 2001, pp. 416–423.
 [20] A. Blake, C. Rother, M. Brown, P. Perez, and P. Torr, “Interactive image segmentation using an adaptive gmmrf model,” in Proc. Eur. Conf. Computer Vision. SpringerVerlag, May 2004.