[R-sig-Geo] plotKML and time stamp

Tomislav Hengl hengl at spatial-analyst.net
Fri Oct 5 10:30:38 CEST 2012 

On 04/10/2012 15:00, Struve,Juliane wrote:
> Hello,
>
> Thanks a lot for replying. I can't run this script, the last line gives me an error:
>
> KML file header opened for parsing...
> Parsing to KML...
> Closing  HRtemp08_jan.kml
> Error in kml_compress(file.name = file.name) :
>    KMZ generation failed. Wrong command passed to 'zip = ... option'.
>
> What do I need to change ?

kmz = FALSE

Looks like your 'system' has no zip option. See:

 > Sys.getenv("R_ZIPCMD", "zip")

https://r-forge.r-project.org/scm/viewvc.php/pkg/R/kml_compress.R?view=markup&root=plotkml

>
> Regards,
>
> Juliane
>
> sessionInfo()
>
> R version 2.15.0 (2012-03-30)
> Platform: i386-pc-mingw32/i386 (32-bit)
>
> locale:
> [1] LC_COLLATE=English_United States.1252
> [2] LC_CTYPE=English_United States.1252
> [3] LC_MONETARY=English_United States.1252
> [4] LC_NUMERIC=C
> [5] LC_TIME=English_United States.1252
>
> attached base packages:
> [1] splines   stats     graphics  grDevices utils     datasets  methods
> [8] base
>
> other attached packages:
>   [1] plotKML_0.2-5      dismo_0.7-23       scales_0.2.2       aqp_1.2-7
>   [5] gstat_1.0-14       spacetime_0.8-0    xts_0.8-6          zoo_1.7-7
>   [9] Hmisc_3.9-3        survival_2.36-12   RColorBrewer_1.0-5 stringr_0.6.1
> [13] colorRamps_2.3     plotrix_3.4-5      colorspace_1.1-1   plyr_1.7.1
> [17] raster_2.0-12      pixmap_0.4-11      XML_3.9-4.1        rgdal_0.7-19
> [21] sp_0.9-99
>
> loaded via a namespace (and not attached):
> [1] cluster_1.14.2  dichromat_1.2-4 grid_2.15.0     labeling_0.1
> [5] lattice_0.20-6  munsell_0.4     reshape_0.8.4   tools_2.15.0
>
>
> -----Original Message-----
> From: Tomislav Hengl [mailto:hengl at spatial-analyst.net]
> Sent: Wednesday, October 03, 2012 7:41 AM
> To: Struve,Juliane
> Cc: Edzer Pebesma; r-sig-geo at r-project.org
> Subject: Re: [R-sig-Geo] plotKML and time stamp
>
>
> try:
>
>   > library(plotKML)
>   > data(HRtemp08)
> # format the time column:
>   > HRtemp08$ctime <- as.POSIXct(HRtemp08$DATE, format="%Y-%m-%dT%H:%M:%SZ")
> ## Not run: # create a STIDF object:
>   > library(spacetime)
>   > sp <- SpatialPoints(HRtemp08[,c("Lon","Lat")])
>   > proj4string(sp) <- CRS("+proj=longlat +datum=WGS84")
>   > HRtemp08.st <- STIDF(sp, time = HRtemp08$ctime, data =
> HRtemp08[,c("NAME","TEMP")])
> # write to a KML file:
>   > HRtemp08_jan <- HRtemp08.st[1:500]
>   > shape = "http://maps.google.com/mapfiles/kml/pal2/icon18.png"
>   > kml(HRtemp08_jan, dtime = 24*3600, labels = TEMP, shape = shape,
> colour = TEMP, colour_scale=rep("#FFFF00", 2), kmz = TRUE)
>
> See also:
>
> http://plotkml.r-forge.r-project.org/tutorial.php#sp_classes
> http://plotkml.r-forge.r-project.org/layer.STIDF.html
>
> Note that, when using 'kml' method, you need to open the KML file yourself.
>
> T. Hengl
>
> On 02/10/2012 15:01, Struve,Juliane wrote:
>> Hello,
>>
>> The attribute I am interested in is a site number. When I use the
>> site number as attribute data the size of the bubble is proportional
>> to the number, which does not really make sense given that the data
>> are just site numbers. But when I make it a constant I lose the site
>> number information. How can I get a plot that shows all the locations
>> with the same size bubble but displays also the site number
>> information, say like a label ? It seems that the attribute data need
>> to be numeric.
>>
>> Thank you for your time,
>>
>> Juliane
>>
>> -----Original Message----- From: Edzer Pebesma
>> [mailto:edzer.pebesma at uni-muenster.de] Sent: Monday, October 01, 2012
>> 4:55 PM To: Struve,Juliane; r-sig-geo at r-project.org Subject: Re:
>> [R-sig-Geo] plotKML and time stamp
>>
>>
>>
>> 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
>>
>> _______________________________________________ R-sig-Geo mailing
>> list R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>

More information about the R-sig-Geo mailing list