我目前正在使用“Akima”插值例程来进行二维线性插值。我目前正在尝试通过排除不良数据点和依赖于它们的插值来尽可能地进行线性插值。我不想做任何样条拟合,只是线性插值。
我可以想出两种使用现有的 akima 包来做到这一点的方法;
- 通过将二维数据集划分为没有丢失数据点的有效子集,然后对每个子集进行插值,然后合并结果。
- 或者将缺失值设置为无意义的值(在我的例子中为-1.0),然后将结果标记为任何插值 NA。不幸的是,插值节点的索引似乎没有返回,所以我必须自己找到这些节点,在这种情况下我应该编写自己的例程。
每一个都有点痛苦,我确信一定有更好的方法,或者必须有一个包来执行上述操作之一,因为我确信这是许多人遇到的常见问题。
非常感谢任何有关替代插值例程或使用 akima interp 的方法的建议。
鲍勃
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)