## CPM: A Deformable Model for Shape Recovery and Segmentation Based on Charged Particles (2004)

### Cached

### Download Links

- [www.cs.rug.nl]
- [www.cs.rug.nl]
- [www.cs.rug.nl]
- [www.cs.rug.nl]
- DBLP

### Other Repositories/Bibliography

Venue: | IEEE TRANS. PATTERN ANAL. MACHINE INTELL |

Citations: | 21 - 3 self |

### BibTeX

@ARTICLE{Jalba04cpm:a,

author = {Andrei C. Jalba and Michael H. F. Wilkinson and Jos B. T. M. Roerdink},

title = {CPM: A Deformable Model for Shape Recovery and Segmentation Based on Charged Particles },

journal = {IEEE TRANS. PATTERN ANAL. MACHINE INTELL},

year = {2004},

volume = {26},

number = {10},

pages = {2004}

}

### OpenURL

### Abstract

A novel, physically motivated deformable model for shape recovery and segmentation is presented. The model, referred to as the charged-particle model (CPM), is inspired by classical electrodynamics and is based on a simulation of charged particles moving in an electrostatic field. The charges are attracted towards the contours of the objects of interest by an electrostatic field, whose sources are computed based on the gradient-magnitude image. The electric field plays the same role as the potential forces in the snake model, while internal interactions are modeled by repulsive Coulomb forces. We demonstrate the flexibility and potential of the model in a wide variety of settings: shape recovery using manual initialization, automatic segmentation, and skeleton computation. We perform a comparative analysis of the proposed model with the active contour model and show that specific problems of the latter are surmounted by our model. The model is easily extendable to 3D and copes well with noisy images.

### Citations

3087 | Snakes: Active contour models
- Kass, Witkin, et al.
- 1987
(Show Context)
Citation Context ...tion, shape recovery, skeleton. æ 1 INTRODUCTION AN important goal in computer vision is to recover the shapes of the objects of interest from visual data. Deformable models introduced by Kass et al. =-=[1]-=- and generalized to the 3D case by Terzopoulos et al. [2] have become extremely popular, offering a unified framework which combines knowledge from geometry, physics, and approximation theory [3]. The... |

1174 |
The Design and Analysis of Spatial Data Structures
- Samet
- 1990
(Show Context)
Citation Context ...t methods used to compute all constituent parts of the model. 3.5.1 Computation of the Coulomb Force An efficient approach to approximate the Coulomb force between free particles is to use a k-d tree =-=[42]-=- data structure to partition space, such that responses to queries about the k nearest neighbors of a particle can be obtained in logarithmic time. At each time step, after the tree is constructed, th... |

1084 | Geodesic active contours
- Caselles, Kimmel, et al.
- 1997
(Show Context)
Citation Context ...he underlying shapes. Several methods have been proposed to address this limitation. Malladi et al. [20] proposed a geometric formulation of the active contour model, based on level sets. As shown in =-=[21]-=-, [4], this method suffers from boundary leakage in the vicinity of blurred edges. In an attempt to overcome this problem, Caselles et al. [21] reformulated the snake evolution as an optimization of t... |

720 | Elastically Deformable Models
- TERZOPOULOS, PLATT, et al.
- 1987
(Show Context)
Citation Context ...mation theory [3]. They have been extensively used in shape recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling [5] and computer animation =-=[6]-=- to texture segmentation [7] and object tracking [8], [9]. A popular deformable model is the snake model [1], describing a closed parametric curve that deforms dynamically and moves towards the desire... |

631 | Shape modeling with front propagation: A level set approach
- Malladi, Sethian, et al.
- 1995
(Show Context)
Citation Context ...filtering for smoothing. All aforementioned deformable models cannot handle topological changes of the underlying shapes. Several methods have been proposed to address this limitation. Malladi et al. =-=[20]-=- proposed a geometric formulation of the active contour model, based on level sets. As shown in [21], [4], this method suffers from boundary leakage in the vicinity of blurred edges. In an attempt to ... |

459 | On active contour models and balloons
- Cohen
- 1991
(Show Context)
Citation Context ...ensive discussion, see [18], [21], [26], [27]. 1. Curve collapse: In the absence of image energy, the curve tends to collapse [28]. One may add an internal pressure term forcing the contour to expand =-=[29]-=-, but this requires the initial contour to be placed completely inside the target object. Moreover, as pointed out in [18], the snake may get outside the contour of the target object in places with we... |

355 | A new Voronoi-based surface reconstruction algorithm
- Amenta, Bern, et al.
- 1998
(Show Context)
Citation Context ... Early efforts relied on heuristics and assumed that the input points are sampled uniformly from the unknown surface. Major breakthroughs were achieved by Hoppe et al. [39] and later by Amenta et al. =-=[40]-=- ([41] for curve reconstruction), who designed algorithms that guarantee correct reconstruction even with nonuniform (but sufficiently dense) samples, and which require the selection of a single param... |

349 |
Computer Simulation Using Particles
- Hockney, Eastwood
(Show Context)
Citation Context ...ducible to OðN log NÞ. 3.5.2 Computation of the Electric Field The electrostatic field (see (7)) is computed using the socalled “particle-particle particle-mesh” (PPPM) method from molecular dynamics =-=[43]-=-. The field is precomputed, and its values are kept in two matrices (with the same sizes as the input image I), one for each component of the 2D vector field. The basic idea of PPPM is to split the Co... |

289 |
Using dynamic programming for solving variational problems in vision
- Amini, Weymouth, et al.
- 1990
(Show Context)
Citation Context ...in an energy functional to be minimized. The internal forces keep the contour smooth, while the external forces attract the snake towards lines, edges, or other low-level image features. Amini et al. =-=[10]-=- pointed out some shortcomings of the original algorithm and proposed an improvement based on discrete dynamic programming. Although their approach is more stable and allows the inclusion of hard cons... |

278 |
Finite-Element Methods for Active Contour Models and Balloons for 2-D and 3-D Images
- Cohen
- 1993
(Show Context)
Citation Context ...field; result. 3. Convergence into boundary concavities: If the external force is weak, the snake is not pulled inside concavities and will stop advancing. Although there are several approaches [29], =-=[33]-=- which try to solve this problem, most of them do not give satisfactory results [18]. The GVF method was shown to remedy this problem to some extent [18], [34], but it still does not work well for thi... |

260 |
Numerical Recipes in C –The art of Scientific Computing, Cambridge Universe
- Press
- 1992
(Show Context)
Citation Context ...LES 1325 time step tjþ1 tj þ t the force F ðriðtÞÞ acting on each particle pi is computed according to (12), in order to obtain the acceleration ai of the particle. Then, one can use Euler’s method =-=[38]-=- to advance the current velocity and position of the particle over the time step. Unfortunately, this method may result in large errors, which accumulate in time. Therefore, we have developed an effic... |

243 |
A Fast Algorithm for Active Contours and Curvature Estimation
- Williams
- 1992
(Show Context)
Citation Context ...inclusion of hard constraints in the energy functional, the method is time consuming and needs careful setting of parameters. A fast, greedy active contour algorithm was proposed by Williams and Shah =-=[11]-=-. Their method retains the improvements of Amini’s algorithm, but is more than an order of magnitude faster. Leymarie and Levine [12] presented a detailed analysis of the snake model, emphasizing its ... |

200 | Surface modeling with oriented particle systems
- Szeliski, Tonnesen
- 1992
(Show Context)
Citation Context ... by Reeves [24] to model natural phenomena such as fire and waterfalls. In these models, particles move under the influence of force fields but do not interact with each other. Szeliski and Tonnensen =-=[25]-=- proposed a system of oriented, interacting particles to model deformable surfaces. Each particle is considered as a surface element and has an associated rotation matrix. The authors define coplanar ... |

154 |
Dynamic programming for detecting, tracking, and matching deformable contours
- Geiger, Gupta, et al.
- 1995
(Show Context)
Citation Context ...pposed to stop when the global energy minimum is reached. However, the snake can be trapped in local minima, and global minimization techniques such as simulated annealing [35] or dynamic programming =-=[36]-=- need to be used. 5. Topological changes: Splitting and merging of contours can be addressed by considering the active contour (or surface) as the zero level set of a higherdimensional function, cf. M... |

151 |
Tracking deformable objects in the plane using an active contour model. Anonymous
- Leymarie, Levine
- 1993
(Show Context)
Citation Context ...pe recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling [5] and computer animation [6] to texture segmentation [7] and object tracking [8], =-=[9]-=-. A popular deformable model is the snake model [1], describing a closed parametric curve that deforms dynamically and moves towards the desired image features under the influence of internal and exte... |

134 |
A geometric model for active contours
- Caselles, Catté, et al.
- 1993
(Show Context)
Citation Context ...s: Splitting and merging of contours can be addressed by considering the active contour (or surface) as the zero level set of a higherdimensional function, cf. Malladi et al. [20] and Caselles et al. =-=[37]-=-. The evolving curve C is embedded as the zero level set of a 2D scalar function ðx; tÞ, satisfying the evolution equation @ @t ^Fiðx; yÞðw0 þ w1 ðx; yÞÞjrj; ð2Þ where w0 and w1 are constants, is th... |

123 | Scale-space theory: A basic tool for analysing structures at different scales. Journal of applied statistics
- Lindeberg
- 1994
(Show Context)
Citation Context ...rticle positions) obtained using the projection (rescaling) of Csþ1 to level s as initialization. 3. Do reconstruction at scale s 0. For the purpose of simplicity we use a standard Gaussian pyramid =-=[45]-=- with three levels, although other pyramids (e.g., wavelet pyramids) may be used as well. The main advantage of this approach is a reduction of the computing costs. Initial convergence leads to a roug... |

91 | Generalized gradient vector flow external forces for active contours
- Xu, Prince
- 1998
(Show Context)
Citation Context ...hough there are several approaches [29], [33] which try to solve this problem, most of them do not give satisfactory results [18]. The GVF method was shown to remedy this problem to some extent [18], =-=[34]-=-, but it still does not work well for thin and long boundary concavities. As the example of Fig. 1 shows, the GVF snake stops at the “entrance” of the concavity because there is no attractive force to... |

79 |
Gradient Vector Flow
- Xu, Prince, et al.
- 1998
(Show Context)
Citation Context ...tour into segments, each of which defines the contour locally. In [17], an active contour model based on the B-spline representation and multiple-stage energy minimization was proposed. Xu and Prince =-=[18]-=- proposed a new external force, which they called gradient vector flow (GVF), and showed that GVF has a large capture range and is able to move the snake into boundary concavities. An accurate and hig... |

73 | Topology adaptive deformable surfaces for medical image volume segmentation
- McInerney, Terzopoulos
- 1999
(Show Context)
Citation Context ...nerney and Terzopoulos 0162-8828/04/$20.00 ß 2004 IEEE Published by the IEEE Computer SocietyJALBA ET AL.: CPM: A DEFORMABLE MODEL FOR SHAPE RECOVERY AND SEGMENTATION BASED ON CHARGED PARTICLES 1321 =-=[22]-=- proposed a new class of deformable models by introducing an affine cell image decomposition. This mechanism induces an iterative reparameterization that enables parametric deformable surfaces to modi... |

70 | Area and length minimizing flows for shape segmentation
- Siddiqi, Lauzière, et al.
- 1998
(Show Context)
Citation Context ...convolution. 2.2 Problems, Variations, and Additions There are several problems related to the traditional snake model, which we briefly mention here. For a more extensive discussion, see [18], [21], =-=[26]-=-, [27]. 1. Curve collapse: In the absence of image energy, the curve tends to collapse [28]. One may add an internal pressure term forcing the contour to expand [29], but this requires the initial con... |

66 |
Particle systems—a technique for modeling a class of fuzzy objects
- Reeves
- 1983
(Show Context)
Citation Context ...ss models in the context of surface reconstruction. A completely different approach to shape deformation is based on particle systems. Such systems have been introduced in computer graphics by Reeves =-=[24]-=- to model natural phenomena such as fire and waterfalls. In these models, particles move under the influence of force fields but do not interact with each other. Szeliski and Tonnensen [25] proposed a... |

58 |
ªDeformable Models in Medical Image Analysis: A
- McInerney, Terzopolous
- 1996
(Show Context)
Citation Context ... al. [1] and generalized to the 3D case by Terzopoulos et al. [2] have become extremely popular, offering a unified framework which combines knowledge from geometry, physics, and approximation theory =-=[3]-=-. They have been extensively used in shape recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling [5] and computer animation [6] to texture se... |

56 | A bayesian approach to dynamic contours through stochastic sampling and simulated annealing
- Storvik
- 1994
(Show Context)
Citation Context ...erion: The deformation is supposed to stop when the global energy minimum is reached. However, the snake can be trapped in local minima, and global minimization techniques such as simulated annealing =-=[35]-=- or dynamic programming [36] need to be used. 5. Topological changes: Splitting and merging of contours can be addressed by considering the active contour (or surface) as the zero level set of a highe... |

41 | Adaptive meshes and shells: Irregular triangulation, discontinuities and hierarchical subdivision
- Vasilescu, Terzopoulos
- 1992
(Show Context)
Citation Context ... introducing an affine cell image decomposition. This mechanism induces an iterative reparameterization that enables parametric deformable surfaces to modify their topology. Vasilescu and Terzopoulos =-=[23]-=- used spring-mass models in the context of surface reconstruction. A completely different approach to shape deformation is based on particle systems. Such systems have been introduced in computer grap... |

36 | Shape recovery algorithms using level sets in 2-D/3-D medical imagery: A state-of-the-art review
- Suri, Liu, et al.
(Show Context)
Citation Context ...ly popular, offering a unified framework which combines knowledge from geometry, physics, and approximation theory [3]. They have been extensively used in shape recovery and medical imaging (see [3], =-=[4]-=- for recent surveys). Other applications range from geometric modeling [5] and computer animation [6] to texture segmentation [7] and object tracking [8], [9]. A popular deformable model is the snake ... |

35 | On the relationship between parametric and geometric active contours
- Xu, Jr, et al.
(Show Context)
Citation Context ...ution. 2.2 Problems, Variations, and Additions There are several problems related to the traditional snake model, which we briefly mention here. For a more extensive discussion, see [18], [21], [26], =-=[27]-=-. 1. Curve collapse: In the absence of image energy, the curve tends to collapse [28]. One may add an internal pressure term forcing the contour to expand [29], but this requires the initial contour t... |

34 |
Discrete smooth interpolation in geometric modeling
- Mallet
- 1992
(Show Context)
Citation Context ...metry, physics, and approximation theory [3]. They have been extensively used in shape recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling =-=[5]-=- and computer animation [6] to texture segmentation [7] and object tracking [8], [9]. A popular deformable model is the snake model [1], describing a closed parametric curve that deforms dynamically a... |

32 | ªBrownian Strings: Segmentating Images with Stochastically Deformable Contours,º
- Grzeszczuk, Levin
- 1997
(Show Context)
Citation Context ...g Number TPAMI-0436-1203. and shortcomings, and proposed an improved termination criterion and a method based on a discrete scale-space representation useful for object tracking. Grzeszczuk and Levin =-=[13]-=- controlled the evolution of the active contour by a simulated annealing process which causes the contour to settle into the global minimum of an image-derived energy function. Peterfreund [14] propos... |

31 |
Multi-resolution algorithm for active contour models
- Leroy, Herliin, et al.
- 1996
(Show Context)
Citation Context ...nsitivity to initialization: The initial contour should be close to the true boundary in order to be attracted. Several techniques have been proposed to rectify this problem, e.g., multiscale methods =-=[32]-=-, [19], pressure forces and distance potentials [29], and, more recently, gradient vector flow (GVF) snakes [18]. Basically, the GVF approach replaces the classical external force with a vector field ... |

29 | Geodesic Active Contours for Supervised Texture Segmentation
- PARAGIOS, DERICHE
- 1999
(Show Context)
Citation Context ... been extensively used in shape recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling [5] and computer animation [6] to texture segmentation =-=[7]-=- and object tracking [8], [9]. A popular deformable model is the snake model [1], describing a closed parametric curve that deforms dynamically and moves towards the desired image features under the i... |

19 | Constraints on Deformable Models: Recovering 3D - Terzopoulos, Witkins, et al. - 1988 |

18 |
Region based strategies for active contour models
- Ronfard
- 1994
(Show Context)
Citation Context ...ut in [18], the snake may get outside the contour of the target object in places with weak response of the gradient operator (boundary leakage). Although the dynamic pressure models proposed in [30], =-=[31]-=- were designed to solve the latter problem, the first still remains. 2. Sensitivity to initialization: The initial contour should be close to the true boundary in order to be attracted. Several techni... |

15 |
The Velocity Snake: Deformable Contour for Tracking in Spatio-Velocity Space
- Peterfreund
- 1999
(Show Context)
Citation Context ... Levin [13] controlled the evolution of the active contour by a simulated annealing process which causes the contour to settle into the global minimum of an image-derived energy function. Peterfreund =-=[14]-=- proposed the velocity snake for boundary tracking of nonrigid objects, by applying velocity control to the evolving contour. Ngoi and Jia [15] presented an active contour model for colour region extr... |

15 | Rethinking classical internal forces for active contour models
- Perrin, Smith
- 2001
(Show Context)
Citation Context ...the traditional snake model, which we briefly mention here. For a more extensive discussion, see [18], [21], [26], [27]. 1. Curve collapse: In the absence of image energy, the curve tends to collapse =-=[28]-=-. One may add an internal pressure term forcing the contour to expand [29], but this requires the initial contour to be placed completely inside the target object. Moreover, as pointed out in [18], th... |

12 |
ªSurface Reconstruction from Unorganized
- Hoppe, DeRose, et al.
- 1992
(Show Context)
Citation Context ...blem studied now for many years. Early efforts relied on heuristics and assumed that the input points are sampled uniformly from the unknown surface. Major breakthroughs were achieved by Hoppe et al. =-=[39]-=- and later by Amenta et al. [40] ([41] for curve reconstruction), who designed algorithms that guarantee correct reconstruction even with nonuniform (but sufficiently dense) samples, and which require... |

10 |
Dynamic-programming for detecting, tracking, and matching deformable contours
- Geiger, Gupta, et al.
- 1995
(Show Context)
Citation Context ...n shape recovery and medical imaging (see [3], [4] for recent surveys). Other applications range from geometric modeling [5] and computer animation [6] to texture segmentation [7] and object tracking =-=[8]-=-, [9]. A popular deformable model is the snake model [1], describing a closed parametric curve that deforms dynamically and moves towards the desired image features under the influence of internal and... |

6 |
Segmented Snake for Contour Detection
- Wong, Yuen, et al.
- 1998
(Show Context)
Citation Context ...ary tracking of nonrigid objects, by applying velocity control to the evolving contour. Ngoi and Jia [15] presented an active contour model for colour region extraction in natural scenes. Wong et al. =-=[16]-=- developed the segmented snake model and demonstrated its ability to handle objects with sharp corners. Their method involves a recursive splitand-merge procedure that divides a contour into segments,... |

3 |
An Active Contour Model for Colour Region Extraction
- Ngoi, Jia
- 1999
(Show Context)
Citation Context ...minimum of an image-derived energy function. Peterfreund [14] proposed the velocity snake for boundary tracking of nonrigid objects, by applying velocity control to the evolving contour. Ngoi and Jia =-=[15]-=- presented an active contour model for colour region extraction in natural scenes. Wong et al. [16] developed the segmented snake model and demonstrated its ability to handle objects with sharp corner... |

3 |
The Crust and the -Skeleton
- Amenta, Bern, et al.
- 1998
(Show Context)
Citation Context ... efforts relied on heuristics and assumed that the input points are sampled uniformly from the unknown surface. Major breakthroughs were achieved by Hoppe et al. [39] and later by Amenta et al. [40] (=-=[41]-=- for curve reconstruction), who designed algorithms that guarantee correct reconstruction even with nonuniform (but sufficiently dense) samples, and which require the selection of a single parameter–t... |

1 |
Fast Active Contour Convergence through Curvature Scale Space Filtering
- Mokhtarian, Mohanna
- 2001
(Show Context)
Citation Context ...range and is able to move the snake into boundary concavities. An accurate and high-speed active contour model based on a reformulation of the internal energy was introduced by Mokhtarian and Mohanna =-=[19]-=-, by removing the curvature part and using curvature scalespace filtering for smoothing. All aforementioned deformable models cannot handle topological changes of the underlying shapes. Several method... |

1 |
Active Region Models for Segmenting
- Ivins, Porrill
- 1994
(Show Context)
Citation Context ...nted out in [18], the snake may get outside the contour of the target object in places with weak response of the gradient operator (boundary leakage). Although the dynamic pressure models proposed in =-=[30]-=-, [31] were designed to solve the latter problem, the first still remains. 2. Sensitivity to initialization: The initial contour should be close to the true boundary in order to be attracted. Several ... |

1 |
An Iterative PPPM Method for Simulating Coulombic
- Beckers, Lowe, et al.
- 1998
(Show Context)
Citation Context ...similarly. In the second step, the charges are spread over a larger neighborhood of grid points, in order to produce a smooth total charge distribution. This step is implemented using the approach in =-=[44]-=-, in which the charges are spread by a diffusion process. The method proceeds by solving Poisson’s equation, (18), on the mesh. Then, the mesh-energy term is computed as a weighted sum over the same g... |