Algorithm to compute SIFT descriptors