Procrustes analysis

In statistics, Procrustes analysis is a form of statistical shape analysis used to analyse the distribution of a set of shapes. The name Procrustes refers to a bandit from Greek mythology who made his victims fit his bed either by stretching their limbs or cutting them off.

Background

Here we just consider objects made up from a finite number $k$ of points in$n$ dimensions, these points are called landmark points.

The shape of object can be considered as a member of an equivalence class formed by removing the translational, rotational and scaling components.

For example, translational components can be removed from an objectby translating the object so that the mean of all the points lies at the origin.Likewise the scale component can be removed by scaling the object so thatthe sum of the squared distances from the points to the origin is 1.

Mathematically: take $k$ points in two dimensions, say

:$((x\_1,y\_1),(x\_2,y\_2),dots,(x\_k,y\_k)),$.

The mean of these points is where

:.

Now translate these points so that themean is translated to the origin ,giving the point .Likewise scale can be removed by finding the size of the object

:

and dividing the points by the scale giving points . Other methods for removing the scale are also used.

Removing the rotational component is more complex. Consider two objects with scale and translation removed, let the points of these be $((x\_1,y\_1),ldots)$, $((w\_1,z\_1),ldots)$. Fix one of these and rotate the other around the origin so that the sum of the squared distances between the points is minimised. A rotation by angle $heta\; ,!$ gives $(u\_1,v\_1)\; =\; (cos\; heta\; w\_1-sin\; heta\; z\_1,sin\; heta\; w\_1+cos\; heta\; z\_1)\; ,!$. The Procrustes distance is

:$d=sqrt\{(u\_1-x\_1)^2+(v\_1-y\_1)^2+cdots\}$,

the distance can be minimised by using a least squares technique to find the angle $heta$ which gives the minimum distance.

Variations

There are many ways of representing the shape of an object. The shape of object can be considered as a member of an equivalence class formed by taking the set of all sets of "k" points in "n" dimensions, that is R^kn and factoring out the set of all translations, rotations and scalings. A particular representation of shape is found by choosing a particular representation of the equivalence class. This will give a manifold of dimension "kn"-4. Procrustes is one method of doing this with particular statistical justification. Bookstein obtains a representation of shape by fixing the position of two pointscalled the bases line. One point will be fixed at the origin and the other at (1,0)the remaining points form the "Bookstein" coordinates.

It is also common to consider "shape and scale" that is with translational and rotational components removed.

Examples

Shape analysis is used in biological data to identify the variations of anatomical features characterised by landmark data, for example in considering the shape of jaw bones. [ [http://www.rps.psu.edu/mar94/goodall.html "Exploring Space Shape"] by Nancy Marie Brown, Research/Penn State, Vol. 15, no. 1, March 1994]

One study by David George Kendall examined the triangles formed by standing stones to deduce if these were often arranged in straight lines. The shape of a triangle can be represented as a point on the sphere, and the distribution of all shapes can be though of a distribution over the sphere.The sample distribution from the standing stones was compared with the theoretical distribution to show that the occurrence of straight lines was no more than average. [ [http://links.jstor.org/sici?sici=0883-4237(198905)4%3A2%3C87%3AASOTST%3E2.0.CO%3B2-P "A Survey of the Statistical Theory of Shape"] , by David G. Kendall, Statistical Science, Vol. 4, No. 2 (May, 1989), pp. 87-99]

See also

*Active shape model
*Alignments of random points
*Biometrics
*Generalized procrustes analysis
*Image registration
*Kent distribution
*Morphometrics
*Orthogonal Procrustes problem
*Procrustes

References

* F.L. Bookstein, "Morphometric tools for landmark data", Cambridge University Press, (1991).
* J.C. Gower, G.B. Dijksterhuis, "Procrustes Problems", Oxford University Press (2004).
* K.V. Mardia, I.L.Dryden, "Statistical Shape Analysis", Wiley, Chichester, (1998).