使用R spacetime加载时空数据集

时间:2016-06-09 14:44:18

标签: r geospatial shapefile

使用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中的列匹配?

0 个答案:

没有答案