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