For non-rigid registration, the objects in medical images are usually treated as a single deformable body with homogeneous stiffness distribution. However, this assumption is invalid for certain parts of the human body, where bony structures move rigidly, while the others may deform. In this paper, we introduce a novel registration technique that models local rigidity of pre-identified rigid structures as well as global non-rigidity in the transformation field using triangular B-splines. In

more »

... rast to the conventional registration method based on tensor-product B-splines, our approach recovers local rigid transformation with fewer degrees of freedom (DOFs), and accurately simulates sharp features (C 0 continuity) along the interface between deformable regions and rigid structures, because of the unique advantages offered by triangular B-splines, such as flexible triangular domain, local control and space-varying smoothness modeling. The accurate matching of the source image with the target one is accomplished through the use of a variational framework, in which a composite energy, measuring the image dissimilarity and enforcing local rigidity and global smoothness, is minimized subject to pre-defined point-based constraints. The algorithm is tested on both synthetic and real 2D images for its applicability. The experimental results show that, by accurately modeling sharp features using triangular B-splines, the deformable regions in the vicinity of rigid structures are less constrained by the global smoothness regularization and therefore contribute extra flexibility to the optimization process. Consequently, the registration quality is improved considerably.