Estimation of coronary artery movement using a non-rigid registration with global-local structure preservation

Background: At present, coronary artery disease (CAD) is the leading cause of death worldwide. Many studies have shown that CAD is strongly associated with the motion characteristics of the coronary arteries. Although cardiovascular imaging technology has been widely used for the diagnosis of CAD, the motion parameters of the heart and coronary arteries cannot be directly calculated from the images. In this paper, we propose a point set registration method with global and local topology constraints to quantify coronary artery movement. Methods: The global constraint is the motion coherence of the point set which enforces the smoothness of the displacement field. The local linear embedding based topological structure and the local feature descriptor i.e., the 3D shape context, are designed to retain the local structure of the point set. We incorporate these constraints into a maximum likelihood framework and derive an expectation-maximization algorithm to obtain the transformation function between the two point sets. The proposed method was compared with four existing algorithms using simulated data and applied to the real data obtained from 4D CT angiograms. Results: For the simulation data, the proposed method achieves a lower registration error than the comparison algorithms. For the real data, the proposed method shows that, in most cases, the right coronary artery achieves a larger velocity than the left anterior descending and left circumflex branches, and there are three well-defined velocity peaks, during the cardiac cycle for these branches. Conclusion: The proposed approach is feasible and effective in quantifying coronary artery movement and thus adds to the diagnostic power of coronary imaging.