Several fixes and refactorings to the interpolation methods

Fixes #388 (closed), unfortunately the implementation of interpolateNANs with limit is not None is quite a bit slower, then the incorrect previous one. Any ideas to speed things up are highly appreciated, I couldn't find a better solution despite a few hours of trying.