Efficient 3D rigid-body registration of micro-MR and micro-CT trabecular bone images

Registration of 3D images acquired from different imaging modalities such as micro-magnetic resonance imaging (µMRI) and micro-computed tomography (µCT) are of interest in a number of medical imaging applications. Most general-purpose multimodality registration algorithms tend to be computationally intensive and do not take advantage of the shape of the imaging volume. Multimodality trabecular bone (TB) images of cylindrical cores, for example, tend to be misaligned along and around the axial direction more than that around other directions. Additionally, TB images acquired by µMRI can differ substantially from those acquired by µCT due to apparent trabecular thickening from magnetic susceptibility boundary effects and non-linear intensity correspondence. However, they share very similar contrast characteristics since the images essentially represent a binary tomographic system. The directional misalignment and the fundamental similarities of the two types of images can be exploited to achieve fast 3D registration. Here we present an intensity cross-correlation based 3D registration algorithm for registering 3D specimen images from cylindrical cores of cadaveric TB acquired by µMRI and µCT in the context of finite-element modeling to assess the bone's mechanical constants. The algorithm achieves the desired registration by first coarsely approximating the three translational and three rotational parameters required to align the µMR images to the µCT scan coordinate frame and fine-tuning the parameters in the neighborhood of the approximate solution. The algorithm described here is suitable for 3D rigid-body image registration applications where through-plane rotations are known to be relatively small. The accuracy of the technique is constrained by the image resolution and in-plane angular increments used.