Snakes reparameterization for noisy images segmentation and targets tracking Idrissi Sidi Yassine, Samir Belfkih. Lycée Tawfik Elhakim Zawiya de Noaceur, route de Marrakech, Casablanca, maroc. Laboratoire Génie des Systèmes, Ecole Nationale des sciences appliquées, Kenitra, maroc. doctorat.yassine@gmail.com samir.belfkih@hotmail.fr Abstract: Active contours or snakes are very popular tools for image segmentation and video tracking. Unfortunately, the utility of this method is limited in the noisy images, because the noise disturbs the evolution of snake, which leads to detect a wrong object. To overcome this problem, in this paper we propose an efficient method of parameterization; this method gives good results and can be applied in objects tracking. I. Introduction: Active contours or snakes are curves or surfaces that move within an image domain to capture desired image features. The first model of parametric active contour was proposed by Kass et al. [1] and named snakes due to the appearance of contour evolution. This first model doesn t have large capture range: we must initialize the snake close to the true boundaries. Thereafter, numerous methods have been developed to improve the capture range problems, including the balloon force [2], distance potential forces [2], multi-resolution methods [4], gradient vector flow (GVF) [5], its generalization: generalized gradient vector flow (GGVF) [6], and vector field convolution (VFC) [7]. The GVF and VFC snakes are popular due to their abilities of attracting the active contour toward object boundary from a sufficiently large distance and the ability of moving the contour into object cavities, and these two properties doesn t exist in other models. However, the later fields are easily influenced by noise and they have limited success in tackling noise. To address this problem, the authors of [8] demonstrate that a better utilization of structural information is more efficient and can significantly improve noise robustness. There are about two decades [9], it has been proven that multi-resolution active contours reduce noise sensitivity and provide faster convergence. To reduce the effect of high noise levels, the authors of [10] introduce the concept of a multi-scale tensor vector field, which helped to significantly improve the results of segmentation. In [11], we proposed to adapt the parameter controlling the degree of smoothness of the gvf field according to the local variance; the qualitative and quantitative results are satisfactory. In the Sobolev active contours [12], the active contour model is reformulate by redefining the notion of gradient in accordance with Sobolev-type inner products, this method represents a noteworthy advance and gives good results for medium noises. This paper is a continuation of efforts to overcome the difficulty of traditional parametric snake s model to attain the true target in the noisy images case. Our method is different from the existing methods, it exploits the following idea: we can avoid the effect of noise if we make a perturbation of the points of the active contours during its evolution. We will see that this idea is also particularly useful in tracking applications. The rest of this paper is as follows: Section 2 briefly recalls the traditional snake and GVF, VFC methods. In Section 3, we explain how and why the noise disturbs the evolution of snakes, and then we expose in 514
detail our idea to overcome this problem. In section 4, we show and we commit an example of successful tracking in real world video sequences by our method of parametrization. Section 5 gives the conclusions. II. Background 1) Parametric active contours An active contour or snake is represented by a parametric curve =,, 0,1 that deforms and move through the image domain to minimize the energy functional: = 1 2 + + 1 Where α and β are pre-determined weighting parameters controlling the smoothness and tautness of the contour, respectively. denotes the external energy that is the issue mostly discussed in active snake models. The external energy is defined on the entire image domain so that the desired features, e.i. edges, would have lower values. It is easy to show that the that minimizes must satisfy the Euler Lagrange equation: Which can be considered as a force balance equation: =0 (2) # $% +# =0 (3) Where # $% = is the internal force that controls the smoothness of the snake and # = &, is the external force that pulls the snake towards the desired image contour. To find a solution of equation (2), the snake is made dynamic by treating as function of time ' as well as-i.e.,'. Namely, we solve: ) (,'= ),',' & + (4),0= * A numerical solution of the above equation can be achieved by discretizing the equation and solving the discrete system [1, 3]. 2) Gradient vector flow In traditional active contours model [1], the external force is based only on edge information. Therefore the external force has a large magnitude only in the immediate vicinity of the boundary. Since this external force has small capture range, the initial snake must be carefully chosen to be close to the true boundaries. Several researchers proposed alternatives solutions to increase the capture range of the external force (see [2,5,6,7]), the most popular and first solution is, the Gradient Vector Flow (GVF) [5] and its generalization, the GGVF [6].GVF is an external force field,,,, constructed by diffusing the edge map gradient vectors (edge force )-,-. away from edges to the homogeneous regions, at the same time keeping the constructed field as close as possible to the edge force near the edges. This is achieved by minimizing the following energy functional: 2 ( µ ) 2 2 2 2 2 2 2 x y x y x y x y 1 E ( U, V ) = ( U + U + V + V ) + ( f + f ) ( U f ) + ( V f ) dxdy 2 (5) 515
Where is a non-negative parameter expressing the degree of smoothness of the field /, 0. The interpretation of this equation is straightforward; the first integrand keeps the field, /, 0 smooth. The second integrand forces the vector field to resemble the initial edge force near the edges (i.e., where the edge force strength is high). Using the calculus of variation, the GVF can be determined by solving the Euler equations: ( x y )( y ) 2 2 2 µ V f + f V f = (6) ( x y )( x ) 2 2 2 µ U f + f U f = (7) Solving (6) and (7) for (U, V) results in gradient vector flow (GVF) that acts as an external force field for the active contour. Once the GVF force field (U, V) is computed via (6) and (7), he is replacing the external force in the snake evolution equations as follow: 0 0 2 4 X = α d X β d X + U ( X, Y ) 2 2 t ds ds Y d 2 X d 4 X = α β + V ( X, Y ) 2 2 t ds ds (8) (9) The GVF field provides a large capture range and the ability to capture the boundary concavities [5, 6]. He is widely used in the segmentation process by active contours. 3) Vector field convolution The vector field convolution (VFC) is an effective external force for active contour model [7]. In VFC, the external field v is generated by convolving the image edge-map f with an isotropic vector field kernel k. The prefixed vector kernel k is defined as: 1,=2, 4, (10) Where 2, is the magnitude of the vector at,, and 4, is the unit vector pointing to the origin (note that the origin of the discrete kernel is located at the center of the matrix 1). Figure 1: An example of the VFC kernel 5, The direction of vectors of 1 at each position, is given as: 516
4, = 6,. 6 ], 7 =8 + (11) Where,, are positions with respect to the kernel s center. The magnitude 2 of kernel s vector element is given by: 2,=7+ 9: or 2, =exp 6>? (12) If we split 1 into its x-component 1, and its y-component 1.,, the field (VFC) is given by:, =-, 1,=-, 1,,-, 1., (13) Once the VFC force field is computed via (13), it replaces the external force in the snake evolution in the same manner as GVF field. In this method (VFC), the field is calculated using a simple convolution of the edge map with a user-defined vector field kernel. This simplicity combined with its robustness to noise makes the VFC a very promising technique. As compared to gradient vector flow [5], VFC yields improved robustness to noise and initialization. It also allows a more flexible force field definition and a reduced computational cost [7, 14]. III. Active contours reparameterization 1) Problem Two major challenges are faced in active contours: the poor capture range and the high sensitivity towards noise. Let us see how the noise affects the evolution: the GVF and VFC fields are calculated from the edge map, if the image is noisy, we will have parasitical contours, which leads to a messy field, preventing a good segmentation. Illustrate this with the following example: the following figure presents a synthetic image generated by the harmonic equation: 7 = + cos 4 E+F and its GVF field. The same synthetic image is again represented with impulsive noise and its GVF field. Figure 2: (a) synthetic image, (b) GVF of (a), (c) (a) with impulse noise, (d) GVF field of (a) 517
We can notice that the GVF field in (b) and the GVF field in (d) are different. The noise changes the direction of several vectors of the field, which makes the field messy, and consequently, the attraction of the active contour towards the desired edges becomes difficult. 2) Solution By inspecting figure 2, we can easily see that if a point of the active contour reached a noise point, it cannot leave it, because all vectors surrounding the noise points head toward these points. Hence, the idea is to shift points of the active contour, so that they no longer coincide with the noise. To use the previous idea, the active contour must be reparameterized after each iteration of deformation, such a parameterization must be uniform (a same distance between successive points of the active contour). We can make this parametrization as follows: take all points of the active contour G $ = $, $ ), for H = 0,,4, $ introduce ' * =0 and ' $ = LM K L where K $ represents the length of each segment G $9 G $, every ' $ represents the cumulative length of the piecewise line that joins the points from G * to G $. Using the two sets of values {' $, $ } and {' $, $ } as interpolation data, we obtain the two splines P and P., with respect to the independent variable t, that interpolate ' and ', respectively. The parametric curve P'=P ',P. ' is called the cumulative parametric spline, and approximates satisfactorily even those curves with large curvature. Moreover, it can also be proved (see [13]) that it is geometrically invariant. In our approach, we choose a fixed distance between successive points of the active contour after each iteration of deformation, this distance should not be large to well approach the shape, at the same time it should not be small because we want to shift these points to prevent noise points. After that, we perturb these new points by a fixed spacing less than the spacing between two successive points of active contour, these points move together and thus their spacing remains the same. To compare results, we present the initial active contour by 4 equidistant points, and then to estimate the accuracy of the results at each iteration, we calculate the number of pixels between the active contour result and the true shape in the image. We use five versions to compare: a) No perturbation. b) A perturbation with a fixe gap = 2, 3. c) A random perturbation at each iteration in the range 0,2, 0,3. d) A perturbation that maximizes the probability of avoiding the points of noise, i-e, have the highest possible external energy. The result is: Figure 3: The number of pixels between the active contour result and the true shape at each iteration. For all cases, the area decrease rapidly from the initial value, but, when the number of iteration becomes large, it levels off. When no perturbation is used, the number of pixels between the snake and the true shape, levels off at a high value. When perturbations were used, the area continues to drop, because there are more possibilities for increasing the energy. We also see that using a fixed gap for the perturbation for all itérations is preferable. 518
IV. Application in target tracking. In this section, we show how the proposed parameterization improved the tracking; the snake was used to track a single leukocyte from video sequences, the background where the leukocyte rolling is very noisy. In the first frame, the target is manually selected, and then it is refined by an active contour algorithm. For the following frames, the predicted in the previous frame result is used to initialize the snake. We see that the proposed algorithm can be used for tracking real-world scenarios and helps to prevent the divergence of the active contour from the target of interest. Note that without the proposed algorithm, the snake is not smooth, develops loops and the tracking result is also flawed. V. Conclusion We have proposed a fast and simple approach to overcome the problem of noise for parametric active contours. We have shown that the proposed algorithm is more robust to noise. Using this algorithm we are also able in the presence of noise to track targets in real-world scenarios. VI. References [1] M. Kass, A. Witkin, and D. Terzopoulos, Snakes: Active contour models, International Journal of Computer Vision, vol. 1, pp. 321-333, 1987. [2] L. D. Cohen, On active contour models and balloons, CVGIP: Image Understand., vol. 53, pp. 211-218, Mar. 1991. [3] L. D. Cohen and I. Cohen, Finite-element methods for active contour models and balloons for 2-D and 3-D images, IEEE Trans. Pattern Anal. Machine Intell., vol. 15, pp. 1131-1147, Nov. 1993 [4] N. Ray, B. Chanda and J. Das, A fast and flexible multiresolution snake with a definite termination criterion, Patter Recognition, vol. 34, pp.1483-1490, 2001. [5] C. Xu and J. L. Prince, Snakes, shapes, and gradient vector flow, IEEE Transactions on Image Processing, vol. 7, pp. 359, 1998. 519
[6] C. Xu, J.L. Prince, Generalized gradient vector external forces for active contours, Signal Processing, vol. 71 pp. 131-139, 1998. [7] B. Li and S. T. Acton, Active contour external force using vector field convolution for image segmentation, IEEE Transactions on Image Processing, vol. 16, pp. 2096-2097, 2007. [8] A. Kumar, A. Wong, A. Mishra, D. A. Clausi, and Paul Fieguth, Tensor vector field based active contours, in Proc. International Conf. on Image Processing, 2011. [9] B. Leroy, I. Herlin, and L. Cohen, Mul ti-resolution algorithms for active contour models, in ICAOS 96, vol. 219, pp. 58 65. Springer Berlin / Heidelberg, 1996. [10] Kumar, A., A. Wong, A. Mishra, D. A. Clausi, and P. Fieguth, "Tensor vector field based active contours", 18th IEEE International Conference on Image Processing (ICIP 2011), Brussels, Belgium, September, 2011. [11] I. Yassine, S. Belfkih Adaptive Weighting Function for Gradient vector Flow using Local Variance, in Proc International Conference on Multimedia Computing and Systems, 2011. [12] Sundaramoorthi, G., Yezzi, A., & Mennucci, A. (2007). Sobolev active contours. International Journal of Computer Vision, 73(3), 345 366. [13] Su B. and Liu D. (1989) Computational Geometry: Curve and Surface Modeling. Academic Press, New York. [14] Scott T. Acton, N. Ray Biomedical Image Analysis: Segmentation, Morgan and Claypool Publishers, 2009. [15] Scott T. Acton, N. Ray Biomedical Image Analysis: Tracking, Morgan and Claypool Publishers, 2008. 520