National Alliance for Medical Image Computing

The Publication Database hosted by SPL

All Publications | Upload | Advanced Search | Gallery View | Download Statistics | Help | Import | Log in

3D Nonrigid Registration via Optimal Mass Transport on the GPU

Institution:
1Georgia Institute of Technology, Atlanta, GA, 30332, USA
2Emory University, Atlanta, GA, 30332, USA
Publisher:
Elsevier Science
Publication Date:
Dec-2009
Journal:
Med Image Anal
Volume Number:
13
Issue Number:
6
Pages:
931-40
Citation:
Med Image Anal. 2009 Dec;13(6):931-40.
PubMed ID:
19135403
PMCID:
PMC2811327
Keywords:
Non-rigid registration, Optimal Mass Transport, Monge-Kantorovich, Variational Methods, Multigrid, GPU
Appears in Collections:
NAC, NA-MIC
Sponsors:
NFS
P41 RR13218 RR/NCRR NIH HHS/United States
U54 EB005149 EB/NIBIB NIH HHS/United States
Generated Citation:
ur Rehman T., Haber E., Pryor G., Melonakos J., Tannenbaum A. 3D Nonrigid Registration via Optimal Mass Transport on the GPU. Med Image Anal. 2009 Dec;13(6):931-40. PMID: 19135403. PMCID: PMC2811327.
Downloaded: 731 times. [view map]
Paper: Download, View online
Export citation:

In this paper we present a new computationally efficient numerical scheme for the minimizing flow approach for optimal mass transport (OMT) with applications to non-rigid 3D image registration. The approach utilizes all of the gray-scale data in both images, and the optimal mapping from image A to image B is the inverse of the optimal mapping from B to A. Further, no landmarks need to be specified, and the minimizer of the distance functional involved is unique. Our implementation also employs multigrid, and parallel methodologies on a consumer graphics processing unit (GPU) for fast computation. Although computing the optimal map has been shown to be computationally expensive in the past, we show that our approach is orders of magnitude faster then previous work and is capable of finding transport maps with optimality measures (mean curl) previously unattainable by other works (which directly influences the accuracy of registration). We give results where the algorithm was used to compute non-rigid registrations of 3D synthetic data as well as intra-patient pre-operative and post-operative 3D brain MRI datasets.

Additional Material
1 File (212.791kB)
Rehman-MEDIA2008-fig4.jpg (212.791kB)