PDF version of this entire document

Distances in Image Space

The most straightforward way to measure the distance between images is to treat each image as a vector formed by concatenating the pixel/voxel intensity values, then take the Euclidean distance. It means that each pixel/voxel in one image is compared against its spatially corresponding pixel/voxel in another image. Although this has the merit of simplicity, it does not provide a very well-behaved distance measure since it increases rapidly for quite small image misalignments [#!wang!#]. This observation led us to consider an alternative distance measure, based on the 'shuffle difference', inspired by the 'shuffle transform' [#!Shuffle!#]. If we have two images 61#61 and 62#62, then the shuffle distance between them is defined as

63#63 (11)

where there are 64#64 pixels (or voxels) indexed by 65#65, and 66#66 is the set of pixels in a neighbourhood of radius 67#67 around 65#65.

The idea is illustrated in Figure . Instead of taking the sum-of-squared-differences between corresponding pixels, the minimum absolute difference between each pixel in one image and the values in a neighbourhood around the corresponding pixel is used. This is less sensitive to small misalignments, and provides a better-behaved distance measure. The tolerance for misalignment is dependent on the size of the neighbourhood (67#67), as is illustrated in Figure .

Figure: The calculation of a shuffle difference image

68#68

It should be noted that the shuffle distance as defined above depends on the direction in which it is measured (see Figure ), hence is not a true distance. It is trivial to construct a symmetric shuffle distance, by averaging the distance calculated both ways between a pair of images. However, it was found that the improvement obtained using this was not significant, and did not justify the increased computation time. In what follows, we use the asymmetric shuffle distance.

Figure: Examples of the shuffle difference image: from one image to a second image (left), from the second image to the first (centre), and the symmetrical shuffle distance image (right)

69#69

Figure: A comparison between shuffle distance using varying size neighbourhoods (radius 67#67). Left: original image, right: warped image, centre, from the left: shuffle distance with 70#70(Euclidean), 71#71 and 72#72 pixels.

73#73