scipy.ndimage.zoom interpolation is a bottleneck for custom surfaces

[douglase]

Design solution is to use matched surfaces that don't require interpolation, a poppy solution would be to cache the zoomed surfaces and check that they haven't changed (with a hash?).

To speed up the interpolation itself could implement accelerated versions:
https://github.com/EconForge/interpolation.py