检索栅格砖的空间点

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*，值= c（＆＃34; X1＆＃34;，＆＃34; X2＆＃34;））：   设置坐标不能在已经设置的空间对象上完成.....

请问任何建议？

Answer 1

您的代码运行正常。只有在您运行此行两次时才会出现错误：

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

将这些绘制在栅格r

之上

plot(r)
points(pt)