Hi,

I have a map tile from RgoogleMaps that I want to use as a background image for 
an spplot of gridded data, where the data will be plotted with an alpha value < 
1.  I am unable to plot them together. 

I was going to set up an example with meuse.grid, but I couldn't get CRS() to 
work with the meuse spec provided in the rgdal documentation.  So allow me to 
show as simply as I can with my code:

# background image
backgr <- list('grid.raster', myGoogleMapTile, x=x.ftCenter, y=y.ftCenter,
               width=backgr.width.ft, height=backgr.height.ft, 
default.units='native')
my.layout <- list(backgr)  

# data frame df has x, y, and z columns
coordinates(df) <- c("x", "y")
spdf <- as(df, "SpatialPixelsDataFrame")

spplot(spdf)  # plots df only, as expected
spplot(spdf, sp.layout=my.layout)  # only the background image and colorkey show

Thinking that the second call of spplot() is plotting the background image over 
the gridded data, I tried using first=T as shown below, but both ways generated 
an error message.  

backgr <- list('grid.raster', myGoogleMapTile, x=x.ftCenter, y=y.ftCenter,
               width=backgr.width.ft, height=backgr.height.ft, 
default.units='native', first=T)
my.layout <- list(backgr, first=T)

What am I missing here?

Thanks,
Scott Waichler
Pacific Northwest National Laboratory
Richland, WA, USA

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to