Skip to content
Prev 16159 / 29559 Next

plotKML and time stamp

Juliane, try

xc = round(runif(10,-1,1)+27.2, 2)
yc = round(runif(10,-1,1)+-82.5, 2)
xy = cbind(xc, yc)
tc = Sys.time() + 3600 * (1:10)
df = data.frame(z1 = round(5 + rnorm(10), 2), z2 = 20:29)

library(sp)
xy.sp = SpatialPoints(xy)
proj4string(xy.sp) = CRS("+proj=longlat +datum=WGS84")

library(spacetime)
ST = STIDF(xy.sp, tc, df)

library(plotKML)
plotKML(ST)

see also

vignette(spacetime)
On 09/30/2012 03:39 AM, Struve,Juliane wrote: