我需要用(X,Y,Z)数据的R执行插值:(X,Y)是平面上的点,Z是某个属性。此外,我想使用scatterplot3d
或rgl
或plot3D
库将数据绘制为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?