I wonder whether the ``interpolation`` argument is redundant with the ``order`` argument. I think ``order=0`` means nearest-neighbor, so we could simply drop the interpolation keyword, and tell users to set ``order=0`` for nearest neighbor.