R:为什么 st_join 会给出无效时间错误?

Neb*_*oyd 5 left-join geospatial nearest-neighbor r-sf

我正在尝试SpatialPointsDataFrames使用sf::st_join(). 两个文件都已使用转换,st_as_sf()但是当我尝试连接时出现错误

rep(seq_len(nrow(x)),lengths(i)) 错误:'times' 参数无效

在这一点上,我尝试交换 x 和 y 参数,并调整参数的无数变体,但似乎没有任何效果。我已经检查了 的帮助文件sf::st_join(),但没有看到任何关于times参数的信息?所以我不确定它来自哪里以及为什么它不断抛出这个错误......

下面是我的数据集的一个示例,它使用进一步向下的代码产生相同的错误

> head(sf.eSPDF[[1]])
Simple feature collection with 6 features and 8 fields
geometry type:  POINT
dimension:      XY
bbox:           xmin: 35.9699 ymin: -3.74514 xmax: 35.97065 ymax: -3.74474
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs
# A tibble: 6 x 9
  TIME                ELEVATION LATITUDE LONGITUDE DATE       V1                  V2                  Survey            geometry
  <dttm>              <chr>        <dbl>     <dbl> <date>     <dttm>              <dttm>               <dbl>         <POINT [°]>
1 2012-01-20 07:26:05 1018 m       -3.74      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1 (35.97047 -3.74474)
2 2012-01-20 07:27:35 1018 m       -3.74      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1 (35.97057 -3.74486)
3 2012-01-20 07:27:39 1019 m       -3.74      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1  (35.9706 -3.74489)
4 2012-01-20 07:27:47 1020 m       -3.74      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1 (35.97065 -3.74489)
5 2012-01-20 07:28:05 1020 m       -3.74      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1 (35.97035 -3.74498)
6 2012-01-20 07:28:26 1019 m       -3.75      36.0 2012-01-20 2012-01-20 00:00:00 2012-01-31 00:00:00      1  (35.9699 -3.74514)

Run Code Online (Sandbox Code Playgroud)
> head(sf.plt.centr)
Simple feature collection with 6 features and 1 field
geometry type:  POINT
dimension:      XY
bbox:           xmin: 35.75955 ymin: -3.91594 xmax: 36.0933 ymax: -3.401
epsg (SRID):    4326
proj4string:    +proj=longlat +datum=WGS84 +no_defs
    PairID                  geometry
1        1   POINT (36.0933 -3.6731)
42      92 POINT (36.02593 -3.91594)
83     215 POINT (36.06496 -3.75837)
124    225   POINT (35.83156 -3.401)
165    251 POINT (35.75955 -3.54388)
206      2 POINT (36.08752 -3.69128)

Run Code Online (Sandbox Code Playgroud)

下面是我用来检查工作解决方案的代码


sf.eSPDF<-lapply(eSPDF, function(x){
  st_as_sf(as(x, "SpatialPointsDataFrame"))
})


sf.plt.centr<-st_as_sf(as(plt.centr, "SpatialPointsDataFrame"))

x1<-head(sf.eSPDF[[1]])
x2<-head(sf.plt.centr)

check<-st_join(x1, x2, join=st_nn, maxdist = Inf, returnDist = T, progress = TRUE)



Run Code Online (Sandbox Code Playgroud)

如您所见,我要加入的文件是列表中的一个对象。该列表中的所有对象都具有与给定示例相同的结构。最终我想得到一个代码,将sf.plt.centr文件连接到列表中的每个文件。就像是...


big.join<-lapply(sf.eSPDF, function(x){
  st_join('['(x), sf.plt.centr, st_nn, maxdist = Inf, returnDist = T, progress = TRUE)
}


Run Code Online (Sandbox Code Playgroud)