In this paper, a novel multiresolution algorithm for registering multimodal images, using an adaptive Monte Carlo scheme is presented. At each iteration, random solution candidates are generated from a multidimensional solution space of possible geometric transformations, using an adaptive sampling approach. The generated solution candidates are evaluated based on the Pearson type-VII error between the phase moments of the images to determine the solution candidate with the lowest error residual. The multidimensional sampling distribution is refined with each iteration to produce increasingly more plausible solution candidates for the optimal alignment between the images. The proposed algorithm is efficient, robust to local optima, and does not require manual initialization or prior information about the images. Experimental results based on various real-world medical images show that the proposed method is capable of achieving higher registration accuracy than existing multimodal regi...