CT to MR registration of complex deformations in the knee joint through dual quaternion interpolation of rigid transforms.

Purpose: To develop a method that enables CT to MR image registration of complex deformations typically encountered in rotating joints such as the knee joint. Methods: We propose a workflow, denoted Quaternion Interpolated Registration (QIR), consisting of three steps, which uses prior knowledge of tissue properties to initialise deformable registration. First, the rigid skeletal components were individually registered. Next, the deformation of soft tissue was estimated using a dual quaternion-based interpolation method. Finally, the registration was fine-tuned with a rigidity-constrained deformable registration step. The method was applied to paired CT and MR images of the knee of 92 patients. It was compared to registration using B-Splines (BS) and B-Splines with a Rigidity Penalty (BSRP). Registration accuracy was evaluated using the Dice Similarity Coefficient (DSC), Mean Absolute Surface Distance (MASD) and 95th percentile Hausdorff Distance (HD95) on bone, and DSC on muscle and fat. To evaluate the rigidity of bone in the registration, the Jacobian Determinant (JD) was calculated. Results: QIR achieved improved results with 0.86, 0.76 mm and 1.88 mm on the DSC, MASD and HD95 metrics on bone, compared to 0.84, 1.40 mm and 4.99 mm for the BS method and 0.84, 1.40 mm and 3.56 mm for the BSRP method. The average DSC of muscle and fat was 0.78 and 0.86 for the QIR, 0.76 and 0.84 for both BS and BSRP. Comparison of the median and interquartile ranges (IQR) of the JD indicated that the QIR (1.00 median, 0.03 IQR) resulted in higher rigidity in the bone compared to the BS (0.98 median, 0.19 IQR) and BSRP (1.00 median, 0.05 IQR). Conclusion: This study showed that QIR could improve the outcome of complex registration problems, encountered in joints involving rigid and non-rigid bodies such as occur in the knee, as compared to a conventional registration approach.