如何在R中进行回归插值？

Question

我有不同气象站的气候数据。我想使用回归而不是克里金或IDW来插入这些数据。任何人都可以告诉你如何在R中做到这一点？

Answer 1

一些示例数据（您应该将问题包括在内）

library(raster)
r <- raster(system.file("external/test.grd", package="raster"))
ra <- aggregate(r, 10)
d <- na.omit(data.frame(v=values(ra), xyFromCell(ra, 1:ncell(ra))))

适合模特

m <- glm(v ~ ., data=d)

预测到光栅

p <- raster(r)
p <- interpolate(p, m)

删除不需要的区域

p <- mask(p, r)