如何用R函数将XYZ数据转为插值栅格网格并可视化？

步骤1：安装并加载所需工具包 #

我们用akima包做插值（支持线性或样条插值），首次使用先安装：

install.packages("akima")
library(akima)

步骤2：准备原始数据 #

用你提供的数据框，设置随机种子保证结果可复现：

set.seed(123)
df <- data.frame("X" = c(1:10), "Y" = c(3,6,1,7,8,9,4,2,10,5), "Z" = rnorm(10, 5))

步骤3：生成规则网格 #

先确定X和Y的范围，生成等间距的网格点（步长by可调整，越小网格越精细）：

x_grid <- seq(min(df$X), max(df$X), by = 0.5)
y_grid <- seq(min(df$Y), max(df$Y), by = 0.5)

步骤4：插值填充缺失区域 #

用interp函数根据原始数据对规则网格插值，得到完整的Z值矩阵：

# 线性插值（速度快，适合简单场景）
interp_result <- interp(x = df$X, y = df$Y, z = df$Z, 
                        xo = x_grid, yo = y_grid, 
                        linear = TRUE)

# 若要更平滑的结果，可改用样条插值：
# interp_result <- interp(x = df$X, y = df$Y, z = df$Z, 
#                         xo = x_grid, yo = y_grid, 
#                         linear = FALSE)

步骤5：绘制填充等高线图 #

直接把插值结果传入filled.contour即可：

filled.contour(x = interp_result$x, 
               y = interp_result$y, 
               z = interp_result$z,
               main = "插值后的XYZ数据填充等高线图",
               xlab = "X轴",
               ylab = "Y轴")

备选：用薄板样条插值（更平滑） #

如果需要更平滑的插值效果，可以用fields包的薄板样条插值：

install.packages("fields")
library(fields)

# 生成网格数据框
grid_df <- expand.grid(X = x_grid, Y = y_grid)
# 拟合插值模型
tps_model <- Tps(df[,c("X","Y")], df$Z)
# 预测网格点的Z值
grid_df$Z <- predict(tps_model, grid_df)
# 转换为矩阵格式
z_matrix <- matrix(grid_df$Z, nrow = length(x_grid), byrow = FALSE)
# 绘图
filled.contour(x = x_grid, y = y_grid, z = z_matrix,
               main = "薄板样条插值填充等高线图")

内容的提问来源于stack exchange，提问作者ZT_Geo