SIFT flow code, plus our code