我想对一个大矩阵（52017455 x 150）的每一行的NA值做一个简单的线性插值，通常我可以通过使用带有apply函数的zoo包中的na.approx函数来做到这一点，但是，这太慢了，我想知道是否有 Rcpp 或 data.table 函数或 packge 来执行此操作？一般来说，我的数据是遥感时空网格，我尝试在光栅包中使用 approxNA 函数，但速度不够快。
x=rnorm(120)
x[c(3,8,16,22)]=NA
data=data.frame(t(x))
data=data[rep(1,52017455),]
result=apply(data,1,zoo::na.approx)