Image-analysis methods play an important role in helping detect brain changes in diagnosis of AD. In the paper, the symmetric logdomain diffeomorphic demons algorithm, with the properties of symmetry and invertibility, is used to compute the pair-wise registration, whose deformation field is then used to calculate the Riemannian distance between them. The spectral embedding algorithm is performed based on the Riemannian distance matrix to project images onto a low-dimensional space where each image is represented as a point and its neighboring points correspond to images of high anatomical similarity. Finally, the quick shift clustering method is employed in the embedded space to partition the dataset into subgroups.