How to overlay two spatial plots?

[matteodefelice]

I have gridded data and a set of stations and I would like to plot them together with spatialPlot. Is it possible?

[miturbide]

Hi Matteo,
spatialPlot uses spplot (package sp), therefore, you could convert your station data into a SpatialPoints* object and use argument sp.layout like this:

library(visualizeR)
library(sp)
stations <- SpatialPoints(getCoordinates(VALUE_Iberia_pr))
spatialPlot(climatology(EOBS_Iberia_pr), sp.layout = list(stations, first = FALSE, pch = 2, col = "black"))

[matteodefelice]

Thanks!