Mark: > My conclusion is that I need a > nearest-neighbour type interpolation routine. José: > Try Radial Basis Functions. This was my thought as well. Scipy has a handy implementation: http://www.scipy.org/Cookbook/RadialBasisFunctions Zach