inertia.m

find best fitting plane through x,y,z points and calculate quality metrics (coplanarity/colinearity)

The best fit plane is assumed to pass through the point's centre of mass:
i.e. cent = [mean(points(:,1)) mean(points(:,2)) mean(points(:,3))]
INPUT: points - nx3 xyz vertex list
OUTPUT: Dir - directional cosine corresponding to the pole to the best fit plane
M - vertex coplanarity
K - vertex colinearity
USAGE: [Dir,M,K] = inertia(points);
REFERENCES: Fernández, O., 2005. Obtaining a best fitting plane through 3D georeferenced data. Journal of Structural Geology, 27, 855–858
COPYRIGHT: Thomas Seers 03/06/2014