Hello! Here is some code that I'm working on: part of a larger project
cat("ready to krige\n") st.u <- krigeST(index.index~1,w,newdata=pred.grd,modelList=stv) print(str(st.u)) plotKML(st.u,"var1.pred") And here is the output: ready to krige Formal class 'STFDF' [package "spacetime"] with 4 slots ..@ data :'data.frame': 9075 obs. of 1 variable: .. ..$ var1.pred: num [1:9075] 104 104 103 103 103 ... ..@ sp :Formal class 'SpatialPixels' [package "sp"] with 5 slots .. .. ..@ grid :Formal class 'GridTopology' [package "sp"] with 3 slots .. .. .. .. ..@ cellcentre.offset: Named num [1:2] -106 26 .. .. .. .. .. ..- attr(*, "names")= chr [1:2] "x1" "x2" .. .. .. .. ..@ cellsize : Named num [1:2] 0.38 0.38 .. .. .. .. .. ..- attr(*, "names")= chr [1:2] "x1" "x2" .. .. .. .. ..@ cells.dim : Named int [1:2] 33 25 .. .. .. .. .. ..- attr(*, "names")= chr [1:2] "x1" "x2" .. .. ..@ grid.index : int [1:825] 793 794 795 796 797 798 799 800 801 802 ... .. .. ..@ coords : num [1:825, 1:2] -106 -106 -106 -105 -105 ... .. .. .. ..- attr(*, "dimnames")=List of 2 .. .. .. .. ..$ : NULL .. .. .. .. ..$ : chr [1:2] "x1" "x2" .. .. ..@ bbox : num [1:2, 1:2] -106.7 25.8 -94.1 35.3 .. .. .. ..- attr(*, "dimnames")=List of 2 .. .. .. .. ..$ : chr [1:2] "x1" "x2" .. .. .. .. ..$ : chr [1:2] "min" "max" .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots .. .. .. .. ..@ projargs: chr "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0" ..@ time :An xts object on 2001-01-01/2011-01-01 containing: Data: int [1:11, 1] 1 2 3 4 5 6 7 8 9 10 ... - attr(*, "dimnames")=List of 2 ..$ : NULL ..$ : chr "timeIndex" Indexed by objects of class: [Date] TZ: UTC xts Attributes: NULL ..@ endTime: POSIXct[1:11], format: "2002-01-01" "2003-01-01" ... NULL Plotting the first variable on the list KML file opened for writing... Temporal support not indicated... using 'dtime=0'. Writing to KML... Closing var1.pred.kml sh: 1: gnome-open: not found > My question is: why am I getting the "dtime=0" message when I have values for @time and @endTime? Is it because one is xts and the other is POSIX, please? Thanks, Erin [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo