Calculating the 2D motion of lumbar vertebrae using splines

In this study we investigate the use of splines and the ICP method [Besl, P., McKay, N., 1992. A method for registration of 3d shapes. IEEE Transactions on Pattern Analysis and Machine Intelligence 14, 239-256.] for calculating the transformation parameters for a rigid body undergoing planar motion parallel to the image plane. We demonstrate the efficacy of the method by estimating the finite centre of rotation and angle of rotation from lateral flexion/extension radiographs of the lumbar spine. In an in vitro error study, the method displayed an average error of rotation of 0.44 +/- 0.45 degrees, and an average error in FCR calculation of 7.6 +/- 8.5 mm. The method was shown to be superior to that of Crisco et al. [Two-dimensional rigid-body kinematics using image contour registration. Journal of Biomechanics 28(1), 119-124.] and Brinckmann et al. [Quantification of overload injuries of the thoracolumbar spine in persons exposed to heavy physical exertions or vibration at the workplace: Part i - the shape of vertebrae and intervertebral discs - study of a yound, healthy population and a middle-aged control group. Clinical Biomechanics Supplement 1, S5-S83.] for the tests performed here. In general, we believe the use of splines to represent planar shapes to be superior to using digitised curves or landmarks for several reasons. First, with appropriate software, splines require less effort to define and are a compact representation, with most vertebra outlines using less than 30 control points. Second, splines are inherently sub-pixel representations of curves, even if the control points are limited to pixel resolutions. Third, there is a well-defined method (the ICP algorithm) for registering shapes represented as splines. Finally, like digitised curves, splines are able to represent a large class of shapes with little effort, but reduce potential segmentation errors from two dimensions (parallel and perpendicular to the image gradient) to just one (parallel to the image gradient). We have developed an application for performing all the necessary computations which can be downloaded from http://www.claritysmart.com.