我有一个地形图的高度等值线,并想从它生成一个高度图(网格,每个单元格有高度值)与scipy.interpolate.griddata. 我的想法是将等高线的点作为非结构化数据点添加到插值算法中,并请求整个网格的插值值(高度)。在
由于某些原因,griddata()返回的网格大多不是在两条等高线之间插值的,两条等高线之间的网格点得到的值很低,而不是两条等高线值之间的值。在
下面是我的示例20x20图像,其中包含等高线(左)和生成的插值网格(1000x1000)。只有来自等高线的非零值(非黑色像素)被用作插值的输入。“最近”、“线性”和“三次”插值方法也存在同样的问题。你知道为什么我不能得到我想要的海图吗?
源代码可以是found here - http://pastebin.com/umNmSrE2
更新:问题已解决。实际上,错误的输入图像是问题所在(见下面的注释)。新的结果是:
目前没有回答
相关问题 更多 >
编程相关推荐