On the Computation of the Average of Spatial Displacements

Abstract

Many applications in biomechanics and medical imaging call for the analysis of the kinematic errors in a group of patients statistically using the average displacement and the standard deviations from the average. This paper studies the problem of computing the average displacement from a set of given spatial displacements using three types of parametric representations: Euler angles and translation vectors, unit quaternions and translation vectors, and dual quaternions. It has been shown that the use of Euclidean norm in the space of unit quaternions reduces the problem to that of computing the average for each quaternion component separately and independently. While the resulting algorithm is simple, the change of the sign of a unit quaternion could lead to an incorrect result. A novel kinematic measure based on dual quaternions is introduced to capture the separation between two spatial displacement. This kinematic measure is then used to formulate a constrained least squares minimization problem. It has been shown that the problem decomposes into that of finding the optimal translation vector and the optimal unit quaternion. The former is simply the centroid of the set of given translation vectors and the latter can be obtained as the eigenvector corresponding to the least eigenvalue of a 4 × 4 positive definite symmetric matrix. It is found that the weight factor used in combining rotations and translations in the formulation does not play a role in the final outcome. Examples are provided to show the comparisons of these methods.