R:在插值时处理丢失的数据

时间:2017-12-22 22:22:52

标签: r list dataframe interpolation

我需要用(X,Y,Z)数据的R执行插值:(X,Y)是平面上的点,Z是某个属性。此外,我想使用scatterplot3drglplot3D库将数据绘制为3D表面。作为第一个解决方案,我想使用akima库。数据以csv结构形式出现,其中缺少某些单元格(XY)对或Z属性。这是一个最小的例子:

data<-read.csv(text="lat,lon,T
                      40.075,50.096,84.2
                      40.263,50.549,82.2
                      ,,73.6
                      40.328378,50.587409,84.8
                      39.856,50.435,,
                      40.1602,49.49,61.7")

(在本例中,X和Y为纬度和经度,Y为纬度和经度,T为某一深度的温度)。在某些行中缺少X,Y缺少。来自akima的插值例程不处理缺失值,并且将此数据提供给它会产生错误。

library(akima)
s<-interp(data$x, data$y, data$T)

如何可以动态丢弃丢失的数据,或者执行插值,可能使用不同的例程来教导忽略NA?

0 个答案:

没有答案