R: Почему st_join выдает ошибку неверного времени? - PullRequest
Новая
       43

R: Почему st_join выдает ошибку неверного времени?

0 голосов
/ 09 июля 2019

Я пытаюсь присоединиться к 2 SpatialPointsDataFrames путем анализа ближайшего соседа, используя sf::st_join().Оба файла были конвертированы с использованием st_as_sf(), но когда я пытаюсь соединиться, я получаю ошибку

Ошибка в rep (seq_len (nrow (x)), length (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)


> 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)

Ниже приведен код, который я использую для проверки рабочих решений 


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)

Как видите, файл, к которому я хочу присоединиться, является объектом в списке.Все объекты в этом списке имеют структуру, идентичную приведенному примеру.В конце концов я хочу получить код, который присоединяет файл sf.plt.centr к каждому из файлов в списке.Что-то вроде ... 


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