Robust registration of thoracic computed tomography (CT) images is strongly impacted by motion during acquisition, high-density objects, and noise, particularly in lower-dose acquisitions. Despite the enhanced registration speed achieved by popular deep learning (DL) methods, their robustness is often neglected. This study aimed to develop a robust thoracic CT image registration algorithm to address the aforementioned issues. A novel, anatomical structure-aware hierarchical registration. By this method, employing a divide-and-conquer approach, dissimilarity metrics, and regularization terms are selected for different regions based on their distinct image features and motion patterns. These terms are then innovatively reconstructed using the Welsch's function, which allows control over the penalty distribution on the loss values. Subsequently, a novel Welsch parameter update strategy is designed for the task of thoracic CT image registration, enabling dynamic sparsity in registration from coarse to fine levels to accommodate various levels of noise and sliding motion. Moreover, the majorization-minimization (MM) algorithm is used to handle the Welsch terms by constructing surrogate functions based on the current variable values for variable update, thereby reducing the complexity of optimization. Experimental results on publicly available deformable image registration lab four-dimensional CT (DIR-Lab 4DCT) and chronic obstructive pulmonary disease (COPD) datasets with and without noise, showed that our proposed method achieves comparable performance to state-of-the-art methods in noise-free scenarios [1.14 and 1.19 mm compared to 1.14 and 1.35 mm target registration errors (TREs)], while demonstrating superior robustness in the presence of noise (1.78 and 2.38 mm compared to 2.00 and 3.31 mm TREs). Ablation studies also validated the effectiveness of each component in the method. A novel and robust algorithm for thoracic CT image registration has been proposed, which has significant potential for valuable clinical applications, including surgical quantitative imaging.
Read full abstract