Covariance-driven retinal image registration initialized from small sets of landmark correspondences