Least Median Squares Estimation and Estimators
This page contains a description of the least median squares estimation
technique with a sample application. It also describes possible uses of
the technique in computer vision.
Introduction
Least squares estimation is a technique used
to find paramaters for a given equation which is the best fit for a set
of data points. The idea is to minimise squares of the
offsets ("the residuals") of the points from the curve. The sum of the
squares of the offsets is used instead of the offset absolute values
because this allows the residuals to be treated as a continuous
differentiable quantity. However, because squares of the offsets are
used, outlying points can have a disproportionate effect on the fit, a
property which may or may not be desirable depending on the problem at
hand.
Least Mean squares tries to minimise the error over all of the
data points in the set. In the case of computer vision (due to digital
photography) the data in question is almost always noisy. This means
that there is a large number of outliers and least mean squares often
gives a poor fit to the data. One solution to this is to use least median squares estimation.
Least Median Squares Estimation
It is very difficult to define the whole algorithm for least median
squares as a mathamatical formula for a whole data set so an algorithm
that generates a least median of squares solution
for a data set is given below:
1. Pick m random sets of points with size p from the data set where p is the number of paramaters in the equation being solved.
2. for each subset use a method such as Least mean of squares to find a solution for the paramaters from that data set.
3.for each set of paramaters p obtained calculate the median of the squared residuals M with respect to the whole data set:
4. The solution is pj for which Mj is minimal among all mMj 's.
Least median of squares estimation is robust to outliers due to its
high breakdown value of 50%. This is the fraction of outliers that can
be tolerated whilst still returning a good solution. However the high
breakdown value means that LMedS does not cope well with Gaussian
noise.
Examples of Use
LMedS like many estimators can be applied
to problems in computer vision. Possible areas for use are fitting a
line to a set of points from an image, fitting a group of lines to a
set of points (such as a model), finding the optical flow in an image,
one area of particular interest is applying it to the discovery of the
Epipolar geometry of a pair of images to be used in stereo vision.
In this case the paramaters that the algorithm is estimating
are the rotation between the coordinate system of the two images, the
translation vector between the coordinate systems and the scale factor.
The data points that it is using are the pixel values from the image
that may be filtered in some way first.
The result of applying LMedS to the problem generates the following Epipolar Lines
Links
Equations from: http://www-sop.inria.fr/robotvis/personnel/zzhang/Publis/Tutorial-Estim/node25.html
Images from: A Highly Robust Regressor And Its Application In Computer Vision, E.Izquierdo
Further Information: http://www-sop.inria.fr/robotvis/personnel/zzhang/Publis/Tutorial-Estim/node25.html
Further Information:Least median of squares: a suitable
objective function for stock assessment models?, Kyle W. Shertzer and
Michael H. Prager, Available At: Link