如何在R语言Terra包中对SpatRaster对象的缺失值进行插值？

问题背景 #

你通过样本点坐标创建并赋值了一个SpatRaster对象，现在需要对其中标记为0的缺失值进行插值处理。首先需要明确：初始创建的rast对象默认像元值为0，你手动赋值的是有效像元，剩余的0属于缺失值，需要先转为标准缺失值标记NA：

library(terra)
# 将栅格中值为0的像元转为缺失值NA
raster_layer[raster_layer == 0] <- NA

原代码错误原因 #

你之前的代码存在两个核心问题：

• 混用了terra和raster包的对象体系：raster_layer是terra的SpatRaster对象，转矩阵后用raster()创建新对象会因包不兼容报错
• approx()是一维线性插值函数，仅适合单维度序列数据，无法适配栅格的空间插值场景

可行解决方案 #

方案1：基于原始样本点的空间插值（推荐） #

既然你有包含坐标和值的样本数据集Test，直接用这些点做空间插值更精准，无需先手动填充栅格：

# 将Test转换为SpatVector空间点对象（替换x、y为你数据中的实际坐标列名）
test_points <- vect(Test, geom = c("x", "y"))

# 克里金插值（适合空间自相关强的数据）
r_interp_krige <- interpolate(raster_layer, test_points, method = "krige")

# 反距离权重插值（操作简单，适合多数场景）
r_interp_idw <- interpolate(raster_layer, test_points, method = "idw")

# 查看插值结果
plot(r_interp_krige)

方案2：行列维度线性插值 #

如果仅需要在栅格的行/列方向做简单线性补值，可结合terra::app实现：

# 按列方向插值
interp_col <- app(raster_layer, function(x) {
  approx(seq_along(x), x, xout = seq_along(x), na.rm = FALSE)$y
})

# 按行方向插值（先转置处理，再转置回原方向）
interp_row <- t(app(t(raster_layer), function(x) {
  approx(seq_along(x), x, xout = seq_along(x), na.rm = FALSE)$y
}))

# 结合行列插值结果取均值
r_interp_comb <- (interp_col + interp_row) / 2

plot(r_interp_comb)

方案3：邻近像元填充 #

如果需要快速用最近的有效像元填充缺失值，可使用terra::fill：

# 用最近邻像元填充缺失值
r_fill_near <- fill(raster_layer, method = "near")
plot(r_fill_near)

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