插值空间数据

Question

我是一位有用的助手，可以翻译您的文本。

我有一个非常大的数据集，包括经纬度和平均降水量(appt)。我想插值一些额外的数据点，以便创建选择州的平均降水量的ggplot图。但我不确定如何插值数据。使用interp函数的结果产生了NA。

如何在给定经度/纬度时插值额外的数据点（约5,000 appt）？

为什么使用interp函数时会得到NA？

dput：

data <- structure(list(longitude = c(-91.40953, -90.30213, -89.26907, 
  -93.32515, -94.13632), latitude = c(33.403216, 34.505369, 32.438327, 
  36.283347, 39.28965), APPT = c(90.4899996121724, 50.4899996121724, 
  30.4899996121724, 70.4899996121724, 93.4899996121724)), .Names = c("longitude", 
  "latitude", "APPT"), row.names = c(NA, 5L), class = "data.frame")

当前代码：

library(akima)

sp_data <- interp(x = data$longitude,
                  y = data$latitude,
                  z = data$APPT)

dInterp <- data.frame(expand.grid(x = sp_data$x,
                                  y = sp_data$y), z = c(sp_data$z))

Current Output : head(dInterp)

          x        y  z
1 -94.13632 32.43833 NA
2 -94.01152 32.43833 NA
3 -93.88672 32.43833 NA
4 -93.76192 32.43833 NA
5 -93.63711 32.43833 NA
6 -93.51231 32.43833 NA

- Vedda

---

我不确定您确切地需要什么，但我曾经回答过另一个类似输入数据（经度、纬度、雨）的问题，该数据进行了插值处理，然后使用ggplot绘制了图表; 请参见这里。 - Henrik

@Henrik 感谢您的帖子。我已经查看并成功绘制了等高线图。我的目标是向数据集中添加额外的值，以便在地图上绘制时填补许多空洞...我认为 interp() 是我需要的函数。 - Vedda

interp 只进行插值，因此网格上任何超出输入数据的点都将为 NA。要进行外推，请设置 extrap=TRUE。 - kdauria

Answer 1

问题在于akima::interp没有填充每个条目。因此，当您查看数据的唯一“角落”时，您只会看到NA值。您需要“向下滚动”以查看插值值。它只填充存在数据的区域：

library(akima) # should have been part of the question:

 contourplot(z ~ x+y, data=dInterp)

73   -90.14268 32.61400       NA
74   -90.01788 32.61400       NA
75   -89.89308 32.61400       NA
76   -89.76828 32.61400       NA
77   -89.64347 32.61400 40.94645
78   -89.51867 32.61400 37.13339
79   -89.39387 32.61400 33.32033
80   -89.26907 32.61400       NA
81   -94.13632 32.78968       NA
82   -94.01152 32.78968       NA

并且：

112  -90.26748 32.78968       NA
113  -90.14268 32.78968       NA
114  -90.01788 32.78968 51.40291
115  -89.89308 32.78968 47.58984
116  -89.76828 32.78968 43.77678
117  -89.64347 32.78968 39.96372
118  -89.51867 32.78968 36.15065
119  -89.39387 32.78968       NA
120  -89.26907 32.78968       NA
121  -94.13632 32.96535       NA

还有许多值的聚类在该包含1600行的数据框中。

- IRTFM

---

谢谢！在大数据文件中试验后，它似乎正在生成超出纬度/经度范围的点。 我只是想为给定的纬度/经度创建额外的点，以便在地图上绘制时填补空洞。 我是否使用了正确的函数？ - Vedda