2016-07-14 8 views
0

Ich versuche, einige Punkte in einem Raster zu extrahieren und zu zeichnen. Ich habe folgendes:Abrufen räumlicher Punkte von Rasterstein

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) 

Aber die Zeile: coordinates(pt) <- c("X1","X2")

wird diesen Fehler geben, die das Plotten beeinflusst:

Fehler in coordinates<- (*tmp*, Wert = c ("X1 "," X2 ")): Setzen von Koordinaten kann nicht auf Räumliche Objekte erfolgen, wo sie bereits eingestellt wurden.

Alle Ratschläge pl Leichtigkeit?

Antwort

0

Ihr Code funktioniert gut. Der Fehler tritt nur auf, wenn Sie diese Zeile laufen zweimal:

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

diese r eines Rasters auf zeichnen

plot(r) 
points(pt) 
+0

Das sind Punkte von Raster extrahiert werden, aber ich weiß nicht, wie plotten sie auf dem Raster .. – NAmo

+0

Verwenden Sie 'plot' gefolgt von' points' – RobertH