I was told, among other things, to look for previous messages by Oscar that would explain how to do this. Here is Oscar's solution that I found to work for the case where you want to plot two rasters that have different extents. In my case, gmap is a satellite image obtained with RgoogleMaps, and z are the results I want to plot over the background image. I made the alpha values of z less than 1 so that the background image would show through. library(rasterVis) levelplot(z~x*y) + layer_(grid.raster(gmap, y=latCenter, width=width, height=height, default.units='native')) Scott Waichler
image background for spplot of gridded data
2 messages · Waichler, Scott R, Oscar Perpiñan
Hi, For the sake of completeness: just a few days ago I received a comment to a post about spplot and ggmap (http://procomun.wordpress.com/2013/04/24/stamen-maps-with-spplot/). There I was using a SpatialPointsDataFrame object. Because of that comment I posted another example with SpatialPolygonsDataFrame: https://gist.github.com/oscarperpinan/7482848 Best, Oscar. ----------------------------------------------------------------- Oscar Perpi??n Lamigueiro Grupo de Sistemas Fotovoltaicos (IES-UPM) Dpto. Ingenier?a El?ctrica (ETSIDI-UPM) URL: http://oscarperpinan.github.io Twitter: @oscarperpinan 2013/11/20 Waichler, Scott R <Scott.Waichler at pnnl.gov>:
I was told, among other things, to look for previous messages by Oscar that would explain how to do this. Here is Oscar's solution that I found to work for the case where you want to plot two rasters that have different extents. In my case, gmap is a satellite image obtained with RgoogleMaps, and z are the results I want to plot over the background image. I made the alpha values of z less than 1 so that the background image would show through. library(rasterVis) levelplot(z~x*y) + layer_(grid.raster(gmap, y=latCenter, width=width, height=height, default.units='native')) Scott Waichler
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo