## Mean shift: A robust approach toward feature space analysis (2002)

### Cached

### Download Links

- [nichol.as]
- [nichol.as]
- [www.caip.rutgers.edu]
- [www.comaniciu.net]
- [pages.cs.wisc.edu]
- [www.cs.jhu.edu]
- [www.cs.cmu.edu]
- [pages.cs.wisc.edu]
- [www.cs.cmu.edu]
- [www.caip.rutgers.edu]
- [courses.csail.mit.edu]
- [www.cs.rutgers.edu]
- [www.ics.uci.edu]
- [www.wisdom.weizmann.ac.il]
- [kowon.dongseo.ac.kr]
- [www.cs.unr.edu]
- [www.cse.unr.edu]
- [www.cse.unr.edu]
- [www.ics.uci.edu]
- [pages.cs.wisc.edu]
- DBLP

### Other Repositories/Bibliography

Venue: | IEEE Transactions on Pattern Analysis and Machine Intelligence |

Citations: | 1469 - 34 self |

### BibTeX

@ARTICLE{Comaniciu02meanshift:,

author = {Dorin Comaniciu and Peter Meer},

title = {Mean shift: A robust approach toward feature space analysis},

journal = {IEEE Transactions on Pattern Analysis and Machine Intelligence},

year = {2002},

volume = {24},

pages = {603--619}

}

### Years of Citing Articles

### OpenURL

### Abstract

Abstract A general nonparametric technique is proposed for the analysis of a complex multimodal featurespace and to delineate arbitrarily shaped clusters in it. The basic computational module of the technique is an old pattern recognition procedure, the mean shift. We prove for discrete data the convergence ofa recursive mean shift procedure to the nearest stationary point of the underlying density function and thus its utility in detecting the modes of the density. The equivalence of the mean shift procedureto the Nadaraya-Watson estimator from kernel regression and the robust M-estimators of location is also established. Algorithms for two low-level vision tasks, discontinuity preserving smoothing andimage segmentation are described as applications. In these algorithms the only user set parameter is the resolution of the analysis, and either gray level or color images are accepted as input. Extensiveexperimental results illustrate their excellent performance.

### Citations

3921 | Pattern Classification and Scene Analysis - Duda, Hart - 1973 |

2649 | Introduction to Statistical Pattern Recognition”, 2nd edition - Fukunaga - 1990 |

2162 | Density Estimation for Statistics and Data Analysis. Monographson Statistics and Applied Probability - Silverman - 1986 |

2152 |
Algorithms for Clustering Data
- Jain, Dubes
- 1988
(Show Context)
Citation Context ... an objective function that expresses the quality of the decomposition (i.e.,the index of cluster validity). The objective function typically compares the inter- versus intra-cluster variability [30],=-=[28]-=- or evaluates the isolation and connectivity of the delineated clusters [43]. . Finally,since in most of the cases the decomposition is task dependent,top-down information provided by the user or by a... |

1328 |
Finding Groups in Data: An Introduction to Cluster Analysis
- Kaufman, Rousseuw
- 1990
(Show Context)
Citation Context ...mizes an objective function that expresses the quality of the decomposition (i.e.,the index of cluster validity). The objective function typically compares the inter- versus intra-cluster variability =-=[30]-=-,[28] or evaluates the isolation and connectivity of the delineated clusters [43]. . Finally,since in most of the cases the decomposition is task dependent,top-down information provided by the user or... |

1123 | Pfinder: Real-time tracking of the human body
- Wren, Azarbayejani
- 1997
(Show Context)
Citation Context ...4c. The segmentation step does not add a significant overhead to the filtering process. The region representation used by the mean shift segmentation is similar to the blob representation employed in =-=[64]-=-. However, while the blob has a parametric description (multivariate Gaussians in both spatial and color domain), the partition generated by the mean shift is characterized by a nonparametric model. A... |

740 | Nonlinear Programming. Athena Scientific - Bertsekas - 1999 |

667 | Statistical pattern recognition: a review - Jain, Duin, et al. - 2009 |

581 | P.: Real-time tracking of non-rigid objects using mean shift
- Comaniciu, Ramesh, et al.
(Show Context)
Citation Context ... already discusses a simple example. However,by introducing adequate objective functions,the optimization problem can acquire physical meaning in the context of a computer vision task. For example,in =-=[14]-=-,by defining the distance between the distributions of the model and a candidate of the target,nonrigid objects were tracked in an image sequence under severe distortions. The distance was defined at ... |

399 | Kernel Smoothing - M, Jones - 1995 |

385 |
Mean shift, mode seeking, and clustering
- Cheng
- 1995
(Show Context)
Citation Context ...eature space [25, 60, 63]. Our approach to mode detection and clustering is based on the mean shift procedure, proposed in 1975 by Fukunaga and Hostetler [21] and largely forgotten till Cheng’s paper =-=[7]-=- rekindled the interest in it. In spite of its excellent qualities, the mean shift procedure does not seem to be known in the statistical literature. While the book [54, Sec.6.2.2] discusses [21], the... |

365 | Multivariate Density Estimation - Scott - 1992 |

360 | Non-parametric model for background subtraction
- Elgammal, Davis
- 2000
(Show Context)
Citation Context ...ametric toolbox developed in this paper is suitable for a large variety of computer vision tasks where parametric models are less adequate, for example, modeling the background in visual surveillance =-=[18]-=-. The complete solution toward autonomous image segmentation is to combine a bandwidth selection technique (like the ones discussed in Section 3.1) with top-down task related high level information. I... |

296 |
The estimation of the gradient of a density function, with applications in pattern recognition
- Fukunaga, Hostetler
- 1975
(Show Context)
Citation Context ...ineated based on the local structure of the feature space [25],[60],[63]. Our approach to mode detection and clustering is based on the mean shift procedure,proposed in 1975 by Fukunaga and Hostetler =-=[21]-=- and largely forgotten until Cheng's paper [7] rekindled interest in it. In spite of its excellent qualities,the mean shift procedure does not seem to be known in statistical literature. While the boo... |

278 | Robust anisotropic diffusion
- Black, Sapiro, et al.
- 1998
(Show Context)
Citation Context ...stopping criterion and after a sufficiently large number of iterations, the processed image collapses into a flat surface. The connection between anisotropic diffusion and M-estimators is analyzed in =-=[5]-=-. A recently proposed noniterative discontinuity preserving smoothing technique is the bilateral filtering [59]. The relation between bilateral filtering and diffusion based techniques was analyzed in... |

209 | Smoothing Methods in Statistics - Simonoff - 1996 |

187 | Robust analysis of feature spaces: color image segmentation
- Comaniciu, Meer
- 1997
(Show Context)
Citation Context ...y a linear mapping property [65,p.166]. Our first image segmentation algorithm was a straightforward application of the feature space analysis technique to an L*u*v* representation of the color image =-=[11]-=-. The modularity of the segmentation algorithm enabled its integration by other groups to a large variety of applications like image retrieval [1],face tracking [6],object-based video coding for MPEG-... |

171 |
A reliable data-based bandwidth selection method for kernel density estimation
- Sheather, Jones
- 1991
(Show Context)
Citation Context ...s on the Laplacian of the unknown density being estimated,and its performance is not well understood [62,p. 108]. For the univariate case,a reliable method for bandwidth selection is the plug-in rule =-=[53]-=-,which was proven to be superior to leastsquares cross-validation and biased cross-validation [42],[55,p. 46]. Its only assumption is the smoothness of the underlying density. . The second bandwidth s... |

157 | Mean shift analysis and applications
- Comaniciu, Meer
- 1999
(Show Context)
Citation Context ...Psi j=1;2::: and n ^fh;K(j)o j=1;2::: converge, and n ^fh;K(j)oj=1;2::: is also monotonically increasing. The proof is given in the Appendix. The theorem generalizes the result derived differently in =-=[13]-=-, where K was the Epanechnikov kernel, and G the uniform kernel. The theorem remains valid when each data point xi is associated with a nonnegative weight wi. An example of nonconvergence when the ker... |

144 | An adaptive clustering algorithm for image segmentation - Pappas - 1992 |

98 | The variable bandwidth mean shift and data-driven scale selection
- Comaniciu, Ramesh, et al.
- 2001
(Show Context)
Citation Context ...nally, since in most of the cases the decomposition is task dependent, top-down information provided by the user or by an upper-level module can be used to control the kernel bandwidth. We present in =-=[15]-=- a detailed analysis of the bandwidth selection problem. To solve the difficulties generated by the narrow peaks and the tails of the underlying density, two locally adaptive solutions are proposed. O... |

97 | Multiscale image segmentation by integrated edge and region detection
- Tabb, Ahuja
(Show Context)
Citation Context ...ature,we will mention only some whose basic processing relies on the joint domain. In each case,a vector field is defined over the sampling lattice of the image. The attraction force field defined in =-=[57]-=- is computed at each pixel as a vector sum of pairwise affinities between the current pixel and all other pixels,with similarity measured in both spatial and range domains. The region boundaries are t... |

82 |
Adaptive noise smoothing filter for images with signal-dependent noise
- Kuan, Sawchuk, et al.
- 1985
(Show Context)
Citation Context ...ther hand, adaptively reduce the amount of smoothing near abrupt changes in the local structure,i.e.,edges. There are a large variety of approaches to achieve this goal,from adaptive Wiener filtering =-=[31]-=-,to implementing isotropic [50] and anisotropic [44] local diffusion processes, a topic which recently received renewed interest [19],[37], [56]. The diffusion-based techniques,however,do not have a s... |

64 | Distribution Free Decomposition of Multivariate Data
- Comaniciu, Meer
- 1999
(Show Context)
Citation Context ...ks retained for the final classification are marked with red dots. 14can be reliably supported by input domain information. The multimodal feature space analysis technique was discussed in detail in =-=[12]-=-. It was shown experimentally that for a synthetic, bimodal normal distribution the technique achieves a classification error similar to the optimal Bayesian classifier. The behavior of this feature s... |

55 | Inference of Surfaces, 3D Curves, and Junctions From Sparse, Noisy, 3D Data
- Guy, Medioni
- 1997
(Show Context)
Citation Context ...ch are based on in situ optimization. Under this paradigm the solution is obtained by using the input domain to define the optimization problem. The in situ optimization is a very powerful method. In =-=[23]-=- and [58] each input data point was associated with a local field (voting kernel) to produce a more dense structure from where the sought information (salient features, the hyperplane representing the... |

51 | Edge flow: A framework of boundary detection and image segmentation
- Ma, Manjunath
- 1997
(Show Context)
Citation Context ...tor computed at that pixel and projected into the spatial domain. However,in contrast to [57],the mean shift procedure moves in the direction of this vector,away from the boundaries. The edge flow in =-=[34]-=- is obtained at each location for a given set of directions as the magnitude of the gradient of a smoothed image. The boundaries are detected at image locations which encounter two opposite directions... |

50 | Parametric and non-parametric unsupervised cluster analysis
- Roberts
- 1997
(Show Context)
Citation Context ...rs arising from the dominant colors, and a decomposition of the space into elliptical tiles will introduce severe artifacts. Enforcing a Gaussian mixture model over such data is doomed to fail, e.g., =-=[49]-=-, and even the use of a robust approach with contaminated Gaussian densities [67] cannot be satisfactory for such complex cases. Note also that the mixture models require the number of clusters as a p... |

49 | Cluster-based Probability Model and its Application to Image and Texture Processing
- Popat, Picard
- 1997
(Show Context)
Citation Context ...67] cannot be satisfactory for such complex cases. Note also that the mixture models require the number of clusters as a parameter,which raises its own challenges. For example,the method described in =-=[45]-=- proposes several different ways to determine this number. Arbitrarily structured feature spaces can be analyzed only by nonparametric methods since these methods do not have embedded assumptions. Num... |

47 |
Computer vision face tracking as a component of a perceptual user interface, in
- Bradski
- 1998
(Show Context)
Citation Context ... representation of the color image [11]. The modularity of the segmentation algorithm enabled its integration by other groups to a large variety of applications like image retrieval [1],face tracking =-=[6]-=-,object-based video coding for MPEG-4 [22],shapedetectionandrecognition[33],andtextureanalysis [47],to mention only a few. However,since the feature space analysis can be applied unchanged to moderate... |

45 | Comparison of data-driven bandwidth selectors
- Park, Marron
- 1990
(Show Context)
Citation Context ... [62, p.108]. For the univariate case a reliable method for bandwidth selection is the plug-in rule [53], which was proven to be superior to least squares cross validation and biased cross-validation =-=[42]-=-, [55, p.46]. Its only assumption is the smoothness of the underlying density. ffl The second bandwidth selection technique is related to the stability of the decomposition. The bandwidth is taken as ... |

43 | Nonparametric multivariate density estimation: a comparative study
- Hwang, Lay, et al.
- 1994
(Show Context)
Citation Context ...never the feature space has more than (say) six dimensions,the analysis should be approached carefully. Employing projection pursuit,in which the density is analyzed along lower dimensional cuts,e.g.,=-=[27]-=-,is a possibility. To conclude,the mean shift procedure is a valuable computational module whose versatility can make it an important component of any computer vision toolbox. APPENDIX Proof of Theore... |

42 | Deformable shape detection and description via model-based region grouping
- Sclaroff, Liu
- 2001
(Show Context)
Citation Context ...n algorithm enabled its integration by other groups to a large variety of applications like image retrieval [1],face tracking [6],object-based video coding for MPEG-4 [22],shapedetectionandrecognition=-=[33]-=-,andtextureanalysis [47],to mention only a few. However,since the feature space analysis can be applied unchanged to moderately higher dimensional spaces (see Section 5),we subsequently also incorpora... |

36 |
Robust Statistical Procedures
- Huber
- 1996
(Show Context)
Citation Context ...than practical importance. 2.6 Relation to Location M-estimators The M-estimators are a family of robust techniques which can handle data in the presence of severe contaminations, i.e., outliers. See =-=[26]-=-, [32] for introductory surveys. In our context only the problem of location estimation has to be considered. 11Given the data and the scale , will define , the location estimator as (30) where, is a... |

35 |
Edge-preserving smoothers for image processing (with discussion
- CHU, GLAD, et al.
- 1998
(Show Context)
Citation Context ...with the requirements to be satisfied by the objective function . The relation between location M-estimators and kernel density estimation is not well investigated in the statistical literature, only =-=[9]-=- discusses it in the context of an edge preserving smoothing technique. 3 Robust Analysis of Feature Spaces Multimodality and arbitrarily shaped clusters are the defining properties of a real feature ... |

33 |
Robust Regression
- Li
- 1985
(Show Context)
Citation Context ...ractical importance. 2.6 Relation to Location M-estimators The M-estimators are a family of robust techniques which can handle data in the presence of severe contaminations, i.e., outliers. See [26], =-=[32]-=- for introductory surveys. In our context only the problem of location estimation has to be considered. 11Given the data and the scale , will define , the location estimator as (30) where, is a symme... |

32 |
Bilateral Filtering for Gray and Color
- Tomasi, Manduchi
- 1998
(Show Context)
Citation Context ...lat surface. The connection between anisotropic diffusion and M-estimators is analyzed in [5]. A recently proposed noniterative discontinuity preserving smoothing technique is the bilateral filtering =-=[59]-=-. The relation between bilateral filtering and diffusion-based techniques was analyzed in [3]. The bilateral filters also work in the joint spatial-range domain. The data is independently weighted in ... |

29 | Geodesic Active Contours for Supervised Texture Segmentation - PARAGIOS, DERICHE - 1999 |

28 |
Detection of diffuse and specular interface reflections and inter-reflections by color image segmentation
- Bajcsy, Lee, et al.
- 1996
(Show Context)
Citation Context ...4. Optional: Eliminate spatial regions containing less than M pixels. The cluster delineation step can be refined according to a priori information and,thus,physics-based segmentation algorithms,e.g.,=-=[2]-=-,[35],can be incorporated. Since this process is performed on region adjacency graphs,hierarchical techniques like [36] can provide significant speed-up. The effect of the cluster delineation step is ... |

23 |
Transform for Line Recognition : Complexity of Evidence Accumulation
- RISSE, Hough
- 1989
(Show Context)
Citation Context ... the space. The problem of color representation will be discussed in Section 4,but the employed parameterization has to be carefully examined even in a simple case like the Hough space of lines,e.g., =-=[48]-=-,[61]. The presence of a Mahalanobis metric can be accommodated by an adequate choice of the bandwidth matrix (2). In practice,however,it is preferable to have assured that the metric of the feature s... |

22 | Bilateral Filtering and Anisotropic Diffusion: Towards a Unified Viewpoint,” technical report, HP Laboratories
- Barash
- 2000
(Show Context)
Citation Context ... A recently proposed noniterative discontinuity preserving smoothing technique is the bilateral filtering [59]. The relation between bilateral filtering and diffusion-based techniques was analyzed in =-=[3]-=-. The bilateral filters also work in the joint spatial-range domain. The data is independently weighted in the two domains and the center pixel is computed as the weighted average of the window. The f... |

15 |
Data sharpening as a prelude to density estimation
- CHOI, P
- 1999
(Show Context)
Citation Context ...to be known in the statistical literature. While the book [54, Sec.6.2.2] discusses [21], the advantages of employing a mean shift type procedure in density estimation were only recently rediscovered =-=[8]-=-. As will be proven in the sequel a computational module based on the mean shift procedure is an extremely versatile tool for feature space analysis and can provide reliable solutions for many vision ... |

14 | A clustering method based on the estimation of the probability density function and on the skeleton by influence zones, Pattern Recognition Lett
- Herbin, Bonnet, et al.
- 1996
(Show Context)
Citation Context ...ma of the p.d.f., that is,to the modes of the unknown density. Once the location of a mode is determined,the cluster associated with it is delineated based on the local structure of the feature space =-=[25]-=-,[60],[63]. Our approach to mode detection and clustering is based on the mean shift procedure,proposed in 1975 by Fukunaga and Hostetler [21] and largely forgotten until Cheng's paper [7] rekindled i... |

12 |
Adaptive nonlocal filtering: A fast alternative to anisotropic diffusion for image enhancement
- Fischl, Schwartz
- 1999
(Show Context)
Citation Context ... of approaches to achieve this goal,from adaptive Wiener filtering [31],to implementing isotropic [50] and anisotropic [44] local diffusion processes, a topic which recently received renewed interest =-=[19]-=-,[37], [56]. The diffusion-based techniques,however,do not have a straightforward stopping criterion and,after a sufficiently large number of iterations,the processed image collapses into a flat surfa... |

12 | Segmentation and interpretation of multicolored objects with highlights, Comput. Vision Image Understand
- Maxwell
- 2000
(Show Context)
Citation Context ...ptional: Eliminate spatial regions containing less than M pixels. The cluster delineation step can be refined according to a priori information and,thus,physics-based segmentation algorithms,e.g.,[2],=-=[35]-=-,can be incorporated. Since this process is performed on region adjacency graphs,hierarchical techniques like [36] can provide significant speed-up. The effect of the cluster delineation step is shown... |

11 |
Finding Groups in Data. An Introduction to Cluster Analysis
- Kauffman, Rousseeuw
- 1990
(Show Context)
Citation Context ...izes an objective function that expresses the quality of the decomposition (i.e., the index of cluster validity). The objective function typically compares the inter- versus intra-cluster variability =-=[30, 28]-=- or evaluates the isolation and connectivity of the delineated clusters [43]. 15Finally, since in most of the cases the decomposition is task dependent, top-down information provided by the user or b... |

10 |
New interpretation and improvement of the nonlinear anisotropic diffusion for image enhancement
- Monteil, Beghdadi
- 1999
(Show Context)
Citation Context ...pproaches to achieve this goal,from adaptive Wiener filtering [31],to implementing isotropic [50] and anisotropic [44] local diffusion processes, a topic which recently received renewed interest [19],=-=[37]-=-, [56]. The diffusion-based techniques,however,do not have a straightforward stopping criterion and,after a sufficiently large number of iterations,the processed image collapses into a flat surface. T... |

7 | ªStatistical Pattern Recognition: A Review,º - Jain, Duin, et al. - 2000 |

7 | Nonparametric Robust Method for Computer Vision - Comaniciu - 2002 |

6 |
ªPfinder: Real-time Tracking of the Human Body,º
- Wren, Azarbayejani, et al.
- 1997
(Show Context)
Citation Context ...4c. The segmentation step does not add a significant overhead to the filtering process. The region representation used by the mean shift segmentation is similar to the blob representation employed in =-=[64]-=-. However,while the blob has a parametric description (multivariate Gaussians in both spatial and color domain),the partition generated by the mean shift is characterized by a nonparametric model. An ... |

5 |
The cascaded hough transform as an aid
- Tuytelaars
- 1998
(Show Context)
Citation Context ...space. The problem of color representation will be discussed in Section 4,but the employed parameterization has to be carefully examined even in a simple case like the Hough space of lines,e.g., [48],=-=[61]-=-. The presence of a Mahalanobis metric can be accommodated by an adequate choice of the bandwidth matrix (2). In practice,however,it is preferable to have assured that the metric of the feature space ... |