ENH: implement image interpolation