检索光栅砖的空间点

Retrieve spatial points of raster brick

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

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")) 错误: 无法在已设置坐标的空间对象上设置坐标.....

有什么建议吗?

您的代码运行良好。仅当您 运行 此行两次时才会发生此错误:

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

将这些绘制在光栅上 r

plot(r)
points(pt)