image background for spplot of gridded data

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 <alobolistas at 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.perpinan at 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 <alobolistas at 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.perpinan at 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.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