Python二维多项式拟合等价于IRAF的“surface”

Question

我正在处理一些数据，试图创建一个二维多项式拟合，就像IRAF的surfit（参见here）。我有16个数据点分布在网格模式中（即16个不同x和y坐标下的像素值），需要对它们进行拟合以生成1024x1024阵列。我尝试过很多不同的方法，从astropy.modeling或{}开始，但是没有什么能比IRAF的surfit给出正确的结果。我想这是因为我只使用了16个数据点，但这就是我的全部！结果应该如下所示：

但我得到的看起来更像这样：

或者这个：

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