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
