R - Spatstat - Поиск ближайшего соседа по идентификатору с использованием таблицы данных

У меня есть два больших фрейма данных, которые называются intersections(представляющие пересечения уличной системы) и users(представляющие пользователей сети) следующим образом:

intersectionsимеет три колонки: x, yи label_street. Они соответственно представляют положение перекрестка в квадратном окне наблюдения (скажем [0,5] x [0,5]) и улицу, на которой он расположен. Вот пример:

intersections <- data.frame(x=c(0.147674, 0.235356, 0.095337, 0.147674), y=c(0.132956, 0.150813, 0.087345, 0.132956), label_street = c(5,6,5,6))

head(intersections)

            x        y label_street
1    0.147674 0.132956            5
2    0.235356 0.150813            6
3    0.095337 0.087345            5
4    0.147674 0.132956            6

Поскольку перекресток находится на пересечении нескольких улиц, каждая (x,y)комбинация в intersectionsтаблице появляется как минимум дважды, но с разными label_street(например, строки 1 и 4 в предыдущем примере). label_streetНе может быть номер строки (именно поэтому он начинается в 5 в моем примере).

users4 колонки: x, y, label_street, ID. Они соответственно представляют положение пользователя, улицу, на которой он расположен, и уникальны для IDкаждого пользователя. В этом фрейме данных нет дубликатов, так как пользователь находится на уникальной улице и имеет уникальный ID. Вот пример ( IDи label_streetможет не быть номером строки)

users <- data.frame(x = c(0.20428152, 0.17840619, 0.12964668, 0.20423856, 0.19349761, 0.10861251), y = c(0.14448448, 0.13921481, 0.11724543, 0.14447573, 0.14228827, 0.09891443), label_street = c(6,6,5,6,6,5), ID = c(2703, 3460, 4325, 12506, 19753, 21282))


head(users)
              x          y label_street      ID
1    0.20428152 0.14448448            6    2703
2    0.17840619 0.13921481            6    3460
3    0.12964668 0.11724543            5    4325
4    0.20423856 0.14447573            6   12506
5    0.19349761 0.14228827            6   19753
6    0.10861251 0.09891443            5   21282

То , что я хочу сделать , это следующее: для каждой точки (x,y)из intersections, получить IDи расстояние до его ближайшего соседа разделяющего же street_label по прибытиюusers

У меня есть рабочее решение, использующее spatstatфункцию nncrossпоиска ближайшего соседа и plyrфункцию adplyдля работы с данными.

Мое рабочее решение выглядит следующим образом:

1) Напишите определяемую пользователем функцию, которая получает идентификатор и расстояние до ближайшего соседа строки в таблице запросов.

 NN <- function(row,query){
 df <- row
 window <- c(0,5,0,5)   #Need this to convert to ppp objects and compute NN distance using nncross
 NN <- nncross(as.ppp(row[,1:2],window),as.ppp(query[,1:2],window))
 df$NN.ID <- query$ID[NN$which]
 df$dist <- NN$dist
 return(df)
}

2) Примените эту определяемую пользователем функцию построчно к моим «пересечениям» фреймов данных, при этом запрос будет подмножеством пользователей, использующих тот же street_label, что и строка:

result <- adply(intersections, 1, function(row) NN(row, users[users$label_street == row$label_street, ])

Результат на примере следующий:

head(result)
           x           y    label_street     NN.ID         NN.dist
1   0.147674    0.132956               5      4325      0.02391247
2   0.235356    0.150813               6      2703      0.03171236
3   0.095337    0.087345               5     21282      0.01760940
4   0.147674    0.132956               6      3460      0.03136304

Поскольку мои реальные фреймы данных будут огромными, я думаю, что вычисление матриц расстояний для просмотра ближайшего соседа не будет эффективным и adplyбудет медленным. Есть у кого-нибудь идея data.tableподобного решения? Я только сейчас об основах data.tableи всегда считал это очень эффективным по сравнению с plyr.