# modified from spacetime:
library(sp)
library(spacetime)
sp <- cbind(x = c(0,0,1), y = c(0,1,1))
row.names(sp) <- paste("point", 1:nrow(sp), sep="")
sp <- SpatialPoints(sp)
time <- as.POSIXct("2010-08-05") + 3600 * (10:12)
x <- STI(sp, time)
st_as_sftime(x)
# convert a Track object from package trajectories to an sftime object
library(trajectories)
x1_Track <- trajectories::rTrack(n = 100)
x1_Track@data$speed <- sort(rnorm(length(x1_Track)))
x1_sftime <- st_as_sftime(x1_Track)
# convert a Tracks object from package trajectories to an sftime object
x2_Tracks <- trajectories::rTracks(m = 6)
x2_sftime <- st_as_sftime(x2_Tracks)
# convert a TracksCollection object from package trajectories to an sftime object
x3_TracksCollection <- trajectories::rTracksCollection(p = 2, m = 3, n = 50)
x3_sftime <- st_as_sftime(x3_TracksCollection)
# convert an sftime object to an sftime object
st_as_sftime(x3_sftime)
# convert an sf object to an sftime object
g <- st_sfc(st_point(c(1, 2)), st_point(c(1, 3)), st_point(c(2, 3)),
st_point(c(2, 1)), st_point(c(3, 1)))
x4_sf <- st_sf(a = 1:5, g, time = Sys.time() + 1:5)
x4_sftime <- st_as_sftime(x4_sf)
# convert a Tracks object from package trajectories to an sftime object
x5_stars <- stars::read_stars(system.file("nc/bcsd_obs_1999.nc", package = "stars"))
x5_sftime <- st_as_sftime(x5_stars, time_column_name = "time")
# this requires some thought to not accidentally drop time dimensions. For
# example, setting `merge = TRUE` will drop the time dimension and thus throw
# an error:
if (FALSE) {
x5_sftime <- st_as_sftime(x5_stars, merge = TRUE, time_column_name = "time")
}
# convert a data frame to an sftime object
x5_df <-
data.frame(a = 1:5, g, time = Sys.time() + 1:5, stringsAsFactors = FALSE)
x5_sftime <- st_as_sftime(x5_df)
# convert a ppp object to an sftime object (modified from the sf package)
if (require(spatstat.geom)) {
st_as_sftime(gorillas, time_column_name = "date")
}
# convert a psp object to an sftime object (modified from the spatstat.geom
# package)
if (require(spatstat.geom)) {
# modified from spatstat.geom:
x_psp <-
psp(
runif(10), runif(10), runif(10), runif(10), window=owin(),
marks = data.frame(time = Sys.time() + 1:10)
)
st_as_sftime(x_psp, time_column_name = "time")
}
# convert an lpp object to an sftime object (modified from the
# spatstat.linnet package)
if (require(spatstat.geom) && require(spatstat.linnet)) {
# modified from spatstat.linnet:
# letter 'A'
v <- spatstat.geom::ppp(x=(-2):2, y=3*c(0,1,2,1,0), c(-3,3), c(-1,7))
edg <- cbind(1:4, 2:5)
edg <- rbind(edg, c(2,4))
letterA <- spatstat.linnet::linnet(v, edges=edg)
# points on letter A
xx <-
spatstat.geom::ppp(
x=c(-1.5,0,0.5,1.5), y=c(1.5,3,4.5,1.5),
marks = data.frame(time = Sys.time() + 1:4, a = 1:4),
window = spatstat.geom::owin(
xrange = range(c(-1.5,0,0.5,1.5)),
yrange = range(c(1.5,3,4.5,1.5)))
)
x_lpp <- spatstat.linnet::lpp(xx, letterA)
# convert to sftime
st_as_sftime(x_lpp, time_column_name = "time")
}
# convert an sftrack object to an sftime object (modified from sftrack)
if (require(sftrack)) {
# get an sftrack object
data("raccoon")
raccoon$timestamp <- as.POSIXct(raccoon$timestamp, "EST")
burstz <-
list(id = raccoon$animal_id, month = as.POSIXlt(raccoon$timestamp)$mon)
x_sftrack <-
as_sftrack(raccoon,
group = burstz, time = "timestamp",
error = NA, coords = c("longitude", "latitude")
)
# convert to sftime
st_as_sftime(x_sftrack)
}
# convert an sftraj object to an sftime object (modified from sftrack)
if (require(sftrack)) {
# get an sftrack object
data("raccoon")
raccoon$timestamp <- as.POSIXct(raccoon$timestamp, "EST")
burstz <-
list(id = raccoon$animal_id, month = as.POSIXlt(raccoon$timestamp)$mon)
x_sftraj <-
as_sftraj(raccoon,
time = "timestamp",
error = NA, coords = c("longitude", "latitude"),
group = burstz
)
# convert to sftime
st_as_sftime(x_sftraj)
}
# convert a cubble_df object from package cubble to an sftime object
if (requireNamespace("cubble", quietly = TRUE, versionCheck = "0.3.0")) {
# get a cubble_df object
data("climate_aus", package = "cubble")
# convert to sftime
climate_aus_sftime <-
st_as_sftime(climate_aus[1:4, ])
climate_aus_sftime <-
st_as_sftime(cubble::face_temporal(climate_aus)[1:4, ])
}
Run the code above in your browser using DataLab