Try this: rinvmod2Ext <- extend(rinvmod2lonlat, with(bbMap2, c(ll.lon, ur.lon, ll.lat, ur.lat)))
contourplot(rinvmod2Ext, zscaleLog=TRUE, margin=FALSE) + layer(grid.raster(gmap2, x=lonCenter2, y=latCenter2, width=width2, height=height2, default.units='native'), under=TRUE) 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/12/10 Agustin Lobo <alobolis...@gmail.com>: > Yes, it does for me as well: I had messed with the parameters of the > map in the original function. > Here I've made a page showing what I do and what I get: > https://dl.dropboxusercontent.com/u/3180464/contourOngmap_log.html > > Is it possible doing the reverse, overlaying the contour on top of the > map instead of putting the map under the contour? > I need the plume in the context of the country. > > Agus > > > Agus > > On Tue, Dec 10, 2013 at 11:51 AM, Oscar Perpiñan > <oscar.perpi...@gmail.com> wrote: >> It works for me using rinvmod2lonlat (long-lat coordinates). >> >> 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/12/10 Agustin Lobo <alobolis...@gmail.com>: >>> and would it possible a similar plot using contourplot() (or >>> levelplot()) overlaid >>> on a gmap? >>> >>> I've tried the equivalent to the >>> stamen.R function that you refer to: >>> >>> contourplot(rinvmod2goog,zscaleLog=TRUE, >>> at=my.at,colorkey=myColorkey,margin=FALSE,add=TRUE) + >>> layer(grid.raster(gmap2, >>> x=lonCenter2, y=latCenter2, >>> width=width2, height=height2, >>> default.units='native'),under=TRUE) >>> >>> and get no error but just the contour. I've tried with my object in >>> pseudomercator and "lon-lat" geographic coordinates: >>> rinvmod2goog <- projectRaster(from=rinvmod2, crs=CRS("+init=epsg:3857")) >>> rinvmod2lonlat <- projectRaster(from=rinvmod2, crs=CRS("+proj=longlat >>> +ellps=WGS84 +no_defs")) >>> >>> both rinvmod2 and and gmap2 display fine independently. >>> >>> Data: >>> https://dl.dropboxusercontent.com/u/3180464/rinvmod2.rda >>> https://dl.dropboxusercontent.com/u/3180464/gmap2.rda >>> (to be used with load()) >>> >>> Thanks >>> >>> Agus >>> >>> On Wed, Nov 20, 2013 at 6:47 PM, Oscar Perpiñan >>> <oscar.perpi...@gmail.com> wrote: >>>> 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.waich...@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@r-project.org >>>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>>> >>>> _______________________________________________ >>>> R-sig-Geo mailing list >>>> R-sig-Geo@r-project.org >>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >> >> _______________________________________________ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo