在R中的空间数据帧中查找最近的网格值

我正在学习空间数据帧。我创建了一个简单的示例，它基于一组初始点，使用反向距离加权（IDW）在网格上插值。现在我想获得每个初始点的插值。我认为这应该是相对容易的，但很难找到一种优雅的方式。换句话说，我想在“geog2.idw”中为40个点（x，y）中的每个点找到最近的网格值

代码

#根据Brian J.Knaus创建的示例进行调整
#包裹
图书馆（gstat）
要求（gstat）
需要（格子）
需要（光栅）
需要（sp）
##### ##### ##### ##### #####
#生成数据。
种子（1）
x您可以将
x
和
y
转换为
空间点
：
空间点（cbind（x，y））
并在
上使用
功能：

pts <- SpatialPoints(cbind(x, y))
pts.idw <- over(pts, geog2.idw["var1.pred"])

# Writing extracted values on the map:
image(geog2.idw, xlim=c(-125, -110), ylim=c(40, 47))
contour(geog2.idw, add=T, lwd=0.5)
map("state", lwd=1, add=T)
map("county", lty=3, lwd=0.5, add=T)
text(pts, round(pts.idw[,"var1.pred"],4), col = "blue")

pts我收到一个错误：“错误：找不到函数“map””您可以将
x
和
y
转换为
空间点
：
空间点（cbind（x，y））
并在
上使用函数
。
pts <- SpatialPoints(cbind(x, y))
pts.idw <- over(pts, geog2.idw["var1.pred"])

# Writing extracted values on the map:
image(geog2.idw, xlim=c(-125, -110), ylim=c(40, 47))
contour(geog2.idw, add=T, lwd=0.5)
map("state", lwd=1, add=T)
map("county", lty=3, lwd=0.5, add=T)
text(pts, round(pts.idw[,"var1.pred"],4), col = "blue")