A method of aligning a pair of images with (101) a first image and a second image, wherein said images comprise a plurality of intensities corresponding to a domain of points in a D-dimensional space includes identifying (102) feature points on both images using the same criteria, computing (103) a feature vector for each feature point, measuring a feature dissimilarity (104) for each pair of feature vectors, wherein a first feature vector of each pair is associated with a first feature point on the first image, and a second feature vector of each pair is associated with a second feature point on the second image. A correspondence mapping (105) for each pair of feature points is determined using the feature dissimilarity associated with each feature point pair, and an image transformation (106) is defined to align (108) the second image with the first image using one or more pairs of feature points that are least dissimilar.