Извлечение пространственных точек растрового кирпича

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 ("X1 "," X2 ")): установочные координаты не могут быть выполнены на пространственных объектах, где они уже установлены .....

Любые советы pl легкость?

Answer 1

Ваш код работает нормально. Ошибка происходит только при запуске этой строки дважды:

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

Чтобы наносить их поверх растра r

plot(r) 
points(pt)