Mean Shift Analysis and Applications
Dorin Comaniciu Peter Meer
Department of Electrical and Computer Engineering
Rutgers University, Piscataway, NJ 08854-8058, USA
fcomanici, [email protected]
to the average of the data p oints within. Details in the Abstract
image are preserved due to the nonparametric charac-
A nonparametric estimator of density gradient, the
ter of the analysis which do es not assume a priori any
mean shift, is employed in the joint, spatial-range
particular structure for the data.
value domain of gray level and color images for dis-
The pap er is organized as follows. Section 2 dis-
continuity preserving ltering and image segmentation.
cusses the estimation of the density gradient and de-
Prop erties of the mean shift are reviewed and its con-
nes the mean shift vector. The convergence of the
vergence on lattices is proven. The prop osed ltering
mean shift pro cedure is proven in Section 3 for discrete
metho d asso ciates with each pixel in the image the clos-
data. Section 4 de nes the pro cessing principle in the
est lo cal mo de in the density distribution of the joint do-
joint spatial-range domain. Mean shift ltering is ex-
main. Segmentation into a piecewise constant structure
plained and ltering examples are given in Section 5.
requires only one more step, fusion of the regions asso-
The prop osed mean shift segmentation is intro duced
ciated with nearby mo des. The prop osed technique has
and analyzed in Section 6.
two parameters controlling the resolution in the spatial
and range domains. Since convergence is guaranteed,
2 Density Gradient Estimation
the technique does not require the intervention of the
Let fx g b e an arbitrary set of n p oints in the
i
i=1:::n
user to stop the ltering at the desired image quality.
d
d-dimensional Euclidean space R . The multivariate
Several examples, for gray and color images, show the
kernel density estimate obtained with kernel K x and
versatilityof the metho d and compare favorably with
window radius h, computed in the p oint x is de ned as
results describ ed in the literature for the same images.
[12, p.76]
n
X
1 Intro duction
1 x x
i
^
f x = K : 1
d
Low level computer vision tasks are misleadingly dif-
nh h
i=1
cult and often yield unreliable results, since the em-
ployed techniques rely up on the correct choice by the
The optimum kernel yielding minimum mean inte-
user of the tuning parameter values. Today,itisanac-
grated square error MISE is the Epanechnikovkernel
cepted fact in the vision community that the execution