[R-sig-Geo] plotKML and time stamp
Edzer Pebesma
edzer.pebesma at uni-muenster.de
Mon Oct 1 22:55:01 CEST 2012
On 10/01/2012 05:47 PM, Struve,Juliane wrote:
> Edzer,
>
> Thank you for your help, this works nicely. When xy.sp are the points and tc is the date/time stamp, why is df needed ?
>
You can create an STI object with just the space-time points, leaving
out the df, but plotKML does not (yet) have a method to display STI
objects -- it wants STIDF with attributes. To work around this, you
could try giving the attribute a constant value, if you don't have sth
that varies.
> Regards,
>
> Juliane
>
>
> -----Original Message-----
> From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Edzer Pebesma
> Sent: Sunday, September 30, 2012 6:06 AM
> To: r-sig-geo at r-project.org
> Subject: Re: [R-sig-Geo] 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:
>> Dear list,
>>
>> After plotting this toy example on Google Earth I would like to add a date/time stamp to it so that the locations will be shown in sequence using the time slider. I am not quite sure how to approach this, can I add a date/time to my spatial points and then submit it to plotKML ? What format does the date/time stamp have to be in to be plotted on Google Earth ?
>>
>> Many thanks for any suggestions.
>>
>> Regards,
>>
>> Juliane
>>
>>
>> xc = round(runif(10,-1,1)+27.2, 2)
>> yc = round(runif(10,-1,1)+-82.5, 2)
>> xy = cbind(xc, yc)
>> xy.sp = SpatialPoints(xy)
>> df = data.frame(z1 = round(5 + rnorm(10), 2), z2 = 20:29) SP.spdf =
>> SpatialPointsDataFrame(xy.sp, df)
>> plot(SP.spdf)
>> #add a projection
>> proj4string(SP.spdf) = CRS("+proj=longlat +datum=WGS84")
>> plotKML(SP.spdf,kmz = FALSE)
>>
>> [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>
> --
> Edzer Pebesma
> Institute for Geoinformatics (ifgi), University of Münster Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
> http://www.52north.org/geostatistics e.pebesma at wwu.de
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics e.pebesma at wwu.de
More information about the R-sig-Geo
mailing list