python中等高线的高阶插值

Question

你们中有谁知道高阶插值方法(Catmull-Rom样条，三次插值等)在Python中创建2D轮廓？

Matplotlib和OpenCV分别提供了函数measure.find_contours()，contours()和findContours()，但它们都是基于线性插值(也称为行进平方)的，我正在研究一些在Python语言中精度更高的东西，最好是Python.任何指针都将受到高度的感谢。

https://www.dropbox.com/s/orgr2yqhbbk2xnr/test.PNG

在上图中，我试图从f(x，y)=x^3+y^3的标量域中提取iso-value 25。我正在寻找6个点，比线性插值给出的6个红点具有更高的精度。

Answer 1

对于非结构化2D数据(或三角化数据)，您可能会对以下类感兴趣：

http://matplotlib.org/api/tri_api.html?highlight=cubictriinterpolator#matplotlib.tri.CubicTriInterpolator

其从用户定义的Triangulation和三角测量节点处定义的字段提供Clough-Tocher (三次)插值器。还可以通过帮助器类UniformTriRefiner使用它

http://matplotlib.org/api/tri_api.html?highlight=refine_field#matplotlib.tri.UniformTriRefiner.refine_field http://matplotlib.org/mpl_examples/pylab_examples/tricontour_smooth_user.png

然而，自适应插值的选择当然取决于您的数据集。