INERTIAELLIPSE Inertia ellipse of a set of points.
Note: Deprecated! Use equivalentEllipse instead.
ELL = inertiaEllipse(PTS);
where PTS is a N*2 array containing coordinates of N points, computes
the inertia ellipse of the set of points.
The result has the form:
ELL = [XC YC A B THETA],
with XC and YC being the center of mass of the point set, A and B are
the lengths of the inertia ellipse (see below), and THETA is the angle
of the main inertia axis with the horizontal (counted in degrees
between 0 and 180).
A and B are the standard deviations of the point coordinates when
ellipse is aligned with the principal axes.
Example
pts = randn(100, 2);
pts = transformPoint(pts, createScaling(5, 2));
pts = transformPoint(pts, createRotation(pi/6));
pts = transformPoint(pts, createTranslation(3, 4));
ell = inertiaEllipse(pts);
figure(1); clf; hold on;
drawPoint(pts);
drawEllipse(ell, 'linewidth', 2, 'color', 'r');
See also
equivalentEllipse
Package: matgeom