Python 2D多项式拟合等同于IRAF的'surfit'

Question

我正在处理一些试图创建二维多项式拟合的数据，就像IRAF的surfit一样（参见here）。我有16个数据点以网格图案（即16个不同的x坐标和y坐标的像素值）分布，需要适合生成1024x1024阵列。我尝试过一系列不同的方法，从astropy.modeling或scipy.interpolate开始，但与IRAF的surfit相比，没有任何方法可以提供正确的结果。我想这是因为我只使用了16个数据点，但这就是我的全部！结果应如下所示：

但我得到的更像是这样：

或这个：

如果您对如何最好地完成此任务有任何建议，我将非常感谢您的意见！谢谢。