Récupérer des points spatiaux d'une brique raster

Question

J'essaie d'extraire et de tracer certains points sur une trame. Je l'ai fait ce qui suit:

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) 

Mais, la ligne: coordinates(pt) <- c("X1","X2")

est de donner cette erreur qui affecte le tracé:

Erreur dans coordinates<- (*tmp*, valeur = c ("X1 "," X2 ")): les coordonnées de réglage ne peuvent pas être faites sur les objets spatiaux, où ils ont déjà été réglés .....

Tout conseil pl facilité?

Answer 1

Votre code fonctionne correctement. L'erreur ne se produit que lorsque vous exécutez cette ligne deux fois:

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

Pour tracer ces sur une trame r

plot(r) 
points(pt)