R——缺失值的二维数据集插值方法

我目前正在使用“Akima”interp例程来进行2d线性插值。我目前正试图通过排除坏的基准点和依赖于它们的插值，尽可能地进行线性插值。我不想做任何样条曲线拟合，只是线性插值

使用现有的akima软件包，我可以想出两种方法来实现这一点

通过将2d数据集划分为不缺少数据点的有效子集，然后对每个子集进行插值，然后合并结果

或者将缺少的值设置为无意义值（在我的例子中为-1.0），然后在任何插值为NA的地方标记结果。不幸的是，插值节点的索引似乎没有返回，因此我必须自己找到这些节点，在这种情况下，我应该编写自己的例程

每一个都有一点痛苦，我相信一定有更好的方法或者一定有一个包来做上面的事情，因为我相信这是很多人都遇到过的常见问题

对于使用akima interp的替代插值例程或方法的任何建议，我们将不胜感激。

---

鲍勃

你看了

Amelia

套餐了吗？

你看了

Amelia

套餐了吗？

你抱怨“疼痛”。我们在回答一个没有可行示例的问题时会有什么痛苦呢？谢谢您的回答，但这还不是一个具体的代码问题。这是一个一般方法问题，也是一个与一般问题类别的建议相关的问题。您是否检查了空间任务视图？那里有很多2D插值……我认为很多插值例程都假设数据集中没有NAs（我知道gstat有）。因此，我不确定是否有现成的解决方案适合您。您关于“排除坏数据点”的说明很可能隐藏了大量潜在的灾难。如果你知道有不好的值，你怎么能开始插值呢？首先，你对坏值的标准是什么？通常，人们可能首先将数据拟合到函数，然后拒绝，比如说3西格玛数据。不管你喜欢与否，除非你提供一个样本和你想要的结果，否则你不会在这里得到有用的回复。我们在回答一个没有可行示例的问题时会有什么痛苦呢？谢谢您的回答，但这还不是一个具体的代码问题。这是一个一般方法问题，也是一个与一般问题类别的建议相关的问题。您是否检查了空间任务视图？那里有很多2D插值……我认为很多插值例程都假设数据集中没有NAs（我知道gstat有）。因此，我不确定是否有现成的解决方案适合您。您关于“排除坏数据点”的说明很可能隐藏了大量潜在的灾难。如果你知道有不好的值，你怎么能开始插值呢？首先，你对坏值的标准是什么？通常，人们可能首先将数据拟合到函数，然后拒绝，比如说3西格玛数据。不管你喜欢与否，除非你提供一个样本和你想要的结果，否则你在这里不会得到一个有用的回答。你能在回答中更完整一点吗？请解释为什么Amelia软件包是相关的。你能在回答中更完整一点吗？请解释为什么Amelia软件包是相关的。