Comment puis-je tracer des valeurs sur une carte raster?

Question

J'ai une carte raster uniennes certaines plage de valeurs que j'importé

Je veux superposer des points sur ce qui proviennent d'un fichier CSV dans ce format: dput (droplevels (tête (points, 10))):

points <- structure(list(lat = c(37.423333, 37.423333, 35.896667, 32.834722, 
32.834722, 32.834722, 32.834722, 32.834722, 32.834722, 32.834722 
), lon = c(-122.188333, -122.188333, -121.087222, -116.622222, 
-116.622222, -116.622222, -116.622222, -116.622222, -116.622222, 
-116.622222)), .Names = c("lat", "lon"), row.names = c(NA, 10L 
), class = "data.frame") 

etc.

J'ai essayé fortify(prcp) mais qui a conduit à cette erreur:

Error: ggplot2 doesn't know how to deal with data of class RasterLayer

Comment puis-je le faire?

TRY:

library(raster) 
library(sp) 

coordinates(points) = ~lon+lat 
proj4string(points) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0") 


r <- getData("worldclim",var="bio",res=10) 
r<-r[[12]] #Annual Precipitation 
names(r) <- c("Prec") 
r <- crop(r, extent(-130, -60, 20, 60)) 
spplot(r) + layer(panel.points(x, y, col="green", cex=0.1, pch=1), data=points) 

Mais la sortie est tout vert en quelque sorte, même si les points ne doivent avoir que 9000 points.

Answer 1

Essayez cette solution:

library(raster) 
library(sp) 
library(latticeExtra) 

points <- structure(list(lat = c(37.423333, 37.423333, 35.896667, 32.834722, 
32.834722, 32.834722, 32.834722, 32.834722, 32.834722, 32.834722 
), lon = c(-122.188333, -122.188333, -121.087222, -116.622222, 
-116.622222, -116.622222, -116.622222, -116.622222, -116.622222, 
-116.622222)), .Names = c("lat", "lon"), row.names = c(NA, 10L 
), class = "data.frame") 

coordinates(points) = ~lon+lat 
proj4string(points) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0") 

r <- getData("worldclim",var="bio",res=10) 
r<-r[[12]] #Annual Precipitation 
names(r) <- c("Prec") 
r <- crop(r, extent(-130, -60, 20, 60)) 

spplot(r)+spplot(points, col.regions="green")