A parallel approach for alignment of multi-modal gridbased data

Multi-modal registration is still a big challenge i n image processing. In this article we present a va riation of the wellknown fast Fourier transform (fft) accelerated meth ods for finding the alignment between two datasets, i.e. the rigid transformation consisting of a rotation as well as a translation mapping all regions in both datasets belonging together in an optimal manner. Our method can be applied to such multi-modal registration problems as computer tomography (CT)/positron emission tomography (PET) matching, or matching CT data with data obtained by magnetic resonance imaging (MRI). We reformulate the alignment problem into an optimization problem concerning a metric measure. The particular form of the proposed objective function can be exploited to fft-accelerate the translational part of the align ment-problem. Thus the reduced problem can be solved for the three rem aining degrees of freedom of the rotatory part usin g standard optimizers, such as downhill-simplex or Powell’s me thod. A further advantage of our approach is the straight forward parallelization of the objective function’ s computation. Our implementation on a graphic processing unit (GPU) yielded a speedup factor between 5 and 25 depending on the size of the data. The results show, that the application of a GPU can be highly rewarding for all fft accelera ted algorithms.

[1]  Paul N. Swarztrauber,et al.  Multiprocessor FFTs , 1987, Parallel Comput..

[2]  H. Levine Medical Imaging , 2010, Annals of Biomedical Engineering.

[3]  M Thelen,et al.  Imaging of the lungs using 3he MRI: Preliminary clinical experience in 18 patients with and without lung disease , 1997, Journal of magnetic resonance imaging : JMRI.

[4]  Frieder Enzmann,et al.  Evaluation of positron-emission-tomography for visualisation of migration processes in geomaterials , 2008 .

[5]  Paul N. Swarztrauber,et al.  Ordered Fast Fourier Transforms on a Massively Parallel Hypercube Multiprocessor , 1991, J. Parallel Distributed Comput..

[6]  Pablo Chacón,et al.  Fast rotational matching of rigid bodies by fast Fourier transform acceleration of five degrees of freedom. , 2003, Acta crystallographica. Section D, Biological crystallography.

[7]  David R. Haynor,et al.  PET-CT image registration in the chest using free-form deformations , 2003, IEEE Transactions on Medical Imaging.

[8]  Vipin Kumar,et al.  The Scalability of FFT on Parallel Computers , 1993, IEEE Trans. Parallel Distributed Syst..

[9]  B. N. Chatterji,et al.  An FFT-based technique for translation, rotation, and scale-invariant image registration , 1996, IEEE Trans. Image Process..

[10]  K Cowtan,et al.  Fast Fourier feature recognition. , 2001, Acta crystallographica. Section D, Biological crystallography.

[11]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[12]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[13]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[14]  Paul A. Viola,et al.  Alignment by Maximization of Mutual Information , 1997, International Journal of Computer Vision.

[15]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[16]  C. Morandi,et al.  Registration of Translated and Rotated Images Using Finite Fourier Transforms , 1987, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[17]  Satish K. Tripathi,et al.  Parallel and Distributed Computing Handbook , 1995 .

[18]  William H. Press,et al.  Numerical recipes in C , 2002 .

[19]  Wei Pan,et al.  An Adaptable-Multilayer Fractional Fourier Transform Approach for Image Registration , 2009, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[20]  Eldad Haber,et al.  Intensity Gradient Based Registration and Fusion of Multi-modal Images , 2006, MICCAI.

[21]  Guy Marchal,et al.  Multimodality image registration by maximization of mutual information , 1997, IEEE Transactions on Medical Imaging.

[22]  Mark J. Harris,et al.  Parallel Prefix Sum (Scan) with CUDA , 2011 .

[23]  Claudia Leopold,et al.  Parallel and distributed computing , 2000 .