An exemplary embodiment of the present invention includes a method of registering images. The method includes: for each image, determining an optimum intensity threshold set from a plurality of intensity threshold sets that maximizes a variance between classes of each set, segmenting each image using the corresponding determined optimum intensity threshold set, generating mutual information from a joint histogram of at least two of the segmented images, and registering the at least two images using the mutual information. The joint histogram may be generated using a geometry shader of a graphical processing unit.
Multilevel thresholding for mutual information based registration and image registration using a GPU