Fast computation of mutual information in the frequency domain with applications to global multimodal image alignment

TL;DR

This paper introduces a fast, frequency domain-based algorithm for computing mutual information across all displacements, enabling efficient global multimodal image alignment with significant speed improvements over previous methods.

Contribution

The authors develop an asymptotically faster algorithm for mutual information computation in the frequency domain, facilitating global multimodal image alignment.

Findings

Achieves high success rates in rigid transformation recovery across diverse datasets.

Outperforms existing local and deep learning-based methods in accuracy.

GPU implementation yields speed-ups of 100 to over 10,000 times.

Abstract

Multimodal image alignment is the process of finding spatial correspondences between images formed by different imaging techniques or under different conditions, to facilitate heterogeneous data fusion and correlative analysis. The information-theoretic concept of mutual information (MI) is widely used as a similarity measure to guide multimodal alignment processes, where most works have focused on local maximization of MI that typically works well only for small displacements; this points to a need for global maximization of MI, which has previously been computationally infeasible due to the high run-time complexity of existing algorithms. We propose an efficient algorithm for computing MI for all discrete displacements (formalized as the cross-mutual information function (CMIF)), which is based on cross-correlation computed in the frequency domain. We show that the algorithm is equivalent to a direct method while asymptotically superior in terms of run-time. Furthermore, we propose a method for multimodal image alignment for transformation models with few degrees of freedom (e.g. rigid) based on the proposed CMIF-algorithm. We evaluate the efficacy of the proposed method on three distinct benchmark datasets, of aerial images, cytological images, and histological images, and we observe excellent success-rates (in recovering known rigid transformations), overall outperforming alternative methods, including local optimization of MI as well as several recent deep learning-based approaches. We also evaluate the run-times of a GPU implementation of the proposed algorithm and observe speed-ups from 100 to more than 10,000 times for realistic image sizes compared to a GPU implementation of a direct method. Code is shared as open-source at \url{github.com/MIDA-group/globalign}.