我在圆柱坐标（R，z，theta）中有一个规则的网格。在每个网格点，我有一个该点密度的值。我正在寻找关于如何插入密度值以获得3D网格单元内特定点的值的建议，使用其周围的8个网格点的值。是否有一个R包将使这个简单？
我试图寻找答案，但大多数似乎与'网格'绘图包有关。
例如，一些样本数据：R <- c(1,2,3)
z <- c(1,2,3)
th <- c(1,2,3)
dens <- array(rep(1, 3*3*3), dim=c(3, 3, 3))
dens[1,2,1] <- 2
如何在（R，z，th）=（1.5,1.5,1.5）得到dens的值？