来自轮廓线的高度图与SciPy griddata

时间:2013-01-13 18:01:57

标签: python scipy gis interpolation

我有一个带有高度轮廓线的地形图,并希望使用scipy.interpolate.griddata从中生成高度图(网格,每个单元格的高度为值)。我的想法是将轮廓线的点作为非结构化数据点添加到插值算法中,并询问整个网格的插值(高度)。

由于某种原因,griddata()返回一个网格,该网格主要不在两条轮廓线之间插值,两条轮廓线之间的网格点得到的值较低,而不是两条轮廓线值之间的值。

下面是我的示例20x20图像,其中包含轮廓线(左)和生成的插值网格(1000x1000)。仅来自轮廓线的非零值(非黑色像素)用作插值的输入。 “最近”,“线性”和“立方”插值方法仍然存在同样的问题。任何想法为什么我不能得到我期待的heighmap? Image of contour lines and interpolation result

源代码可以是found here - http://pastebin.com/umNmSrE2

更新:问题已解决。实际上输入图像不好是问题(见下面的评论)。新结果是:

enter image description here

0 个答案:

没有答案