检索光栅砖的空间点

Question

我正在尝试在栅格上提取和绘制一些点。 我做了以下事情：

library(sp)
pt1 <- data.frame(cbind(4110000,40000))####almost completely degraded
pt2 <- data.frame(cbind(4100000,0))#####green
pt3 <-  data.frame(cbind(4140000,55000))####slightly green
pt <- rbind(pt1, pt2, pt3)
coordinates(pt) <- c("X1","X2")

points(pt)

但是，这行： coordinates(pt) <- c("X1","X2")

给出了影响绘图的错误：

coordinates<-(*tmp*, value = c("X1", "X2")) 中的错误： 无法在已设置的空间对象上设置坐标.....

有什么建议吗？

Answer 1

您的代码运行良好。该错误仅在您运行此行两次时发生：

coordinates(pt) <- c("X1","X2")

在栅格上绘制这些r

plot(r)
points(pt)