使用R包spacetime
,我将创建一个类STFDF
的数据集。我在m个位置(n x m)观察了n次。这些地点是美国各州。
我设法加载了n次:
"2000-01-01 EST" "2000-02-01 EST" "2000-03-01 EST" "2000-04-01 EST" ...
一个有n列和m行的矩阵。我还将美国各州加载到SpatialPolygonsDataFrame
:
df = loadMyDataFrame( ... ) # my US dataset
times = extractTimes( df )
matrix = extractMatrix( df )
# load US states from a shapefile
library(maptools)
crswgs84=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
# class: SpatialPolygonsDataFrame
states <- readShapePoly("data/cb_2015_us_state_500k/cb_2015_us_state_500k",proj4string=crswgs84)
library(spacetime)
STFDF(states, times, data=matrix)
这样做,我得到:
Error: nrow(object@data) == length(object@sp) * nrow(object@time) is not TRUE
显然,states
未正确排序和过滤。如何确保states
中的行与matrix
中的列匹配?