In this paper, a novel method for simultaneously registering multiple images acquired from different imaging modalities is presented. The optimal alignment is computed as the set of transformations that minimize the dispersion of the multi-dimensional joint phase moment distribution. Dispersion is measured as the cumulative quadratic orthogonal distance between samples from the joint phase moment distribution and the corresponding multi-dimensional fitted hyperline. The proposed method is designed to be computationally efficient, robust to signal non-homogeneities and noise, and maintains internal consistency amongst all images. Experimental results using real-world medical and remote sensing images show that the proposed method achieves a high level of registration accuracy when simultaneously registering multiple multimodal images.