The identification of tumors in the internal organs of chest, abdomen, and pelvis anatomic regions can be performed with the analysis of Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) data. The contrast agent is accumulated differently by pathologic and healthy tissues and that results in a temporally varying contrast in an image series. The internal organs are also subject to potentially extensive movements mainly due to breathing, heart beat, and peristalsis. This contributes to making the analysis of DCE-MRI datasets challenging as well as time consuming. To address this problem we propose a novel pairwise non-rigid registration method with a Non-Parametric Bayesian Registration (NParBR) formulation. The NParBR method uses a Bayesian formulation that assumes a model for the effect of the distortion on the joint intensity statistics, a non-parametric prior for the restored statistics, and also applies a spatial regularization for the estimated registration with Gaussian filtering. A minimally biased intra-dataset atlas is computed for each dataset and used as reference for the registration of the time series. The time series registration method has been tested with 20 datasets of liver, lungs, intestines, and prostate. It has been compared to the B-Splines and to the SyN methods with results that demonstrate that the proposed method improves both accuracy and efficiency.

A brain MRI protocol typically includes several imaging contrasts that can provide complementary information by highlighting different tissue properties. The acquired datasets often need to be co-registered or placed in a standard anatomic space before any further processing. Current registration methods particularly for multicontrast data are computationally very intensive, their resolution is lower than that of the images, and their distance metric and its optimization can be limiting. In this work a novel and effective non-rigid registration method is proposed that is based on the restoration of the joint statistics of pairs of such images. The registration is performed with the deconvolution of the joint statistics with an adaptive Wiener filter. The deconvolved statistics are forced back to the spatial domain to estimate a preliminary registration. The spatial transformation is also regularized with Gaussian spatial smoothing. The registration method has been compared with the B-Splines method implemented in 3DSlicer and with the SyN method implemented in the ANTs toolkit. The validation has been performed with a simulated Shepp-Logan phantom, a BrainWeb phantom, the real data of the NIREP database, and real multi-contrast datasets of healthy volunteers. The proposed method has shown improved comparative accuracy as well as analytical efficiency.