Dear r sig geo users,
I trying to add location points (longitude-time) onto a hovmoller plot
(rasterVis). These are locations of objects at given longitude and time.
An example of what I want to do here
library("rasterVis")
library("lattice")
r <- raster(nrows=10, ncols=10)
r <- setValues(r, 1:ncell(r))
r1 <- raster(nrows=10, ncols=10)
r1 <- setValues(r1, 1:ncell(r))
r2 <- raster(nrows=10, ncols=10)
r2 <- setValues(r2, 1:ncell(r))
St=stack(r,r1,r2)
idx <- seq(as.Date('2008-01-15'), as.Date('2008-1-17'), 'day')
SISmm <- setZ(St, idx)
hov<-hovmoller(SISmm, contour=FALSE, panel=panel.levelplot.raster,
yscale.components=yscale.
raster.subticks,interpolate=TRUE,
par.settings=RdBuTheme)
#location data in decimal degrees. these are the points I want to add onto
the hovmoller plot
x<-c(50.00000,25.00000,0.00000)
date<-as.Date('2008-01-15', '2008-01-15','2008-01-17')
loc<-data.frame(x,date)
#attempt to overlay the location onto the hov plot. This is clearly not the
correct approach but I am still at the very bottom of the R learning
curve...
plot(hov, position=c(0,0,1,1))
plot(loc,position=c(0,0,1,1),newpage=F)
Thankful for any pointers or suggestions.
Regards,
Don
[[alternative HTML version deleted]]