Plotting un ráster detrás de un archivo de forma
¿Cómo puedo trazar un objeto "raster" detrás de un objeto shapefile? Ambos trazan bien por sí mismos, pero los puntos no se trazan sobre el ráster:
require(rgdal)
require(maptools)
require(raster)
myproj = "+proj=utm +zone=12 +north +ellps=WGS84 +units=m"
shp = readShapeSpatial(fn.shp, proj4string = CRS(myproj))
ras = raster(fn.tif)
plot(ras)
plot(shp, bg="transparent", add=TRUE)