Совпадение пространственных точек с треугольниками Делоне
У меня есть 2 файла SpatialPointsDataFrame. Я вычисляю треугольники Делоне из первого шейп-файла (shp1
) как в коде ниже. Затем я хочу назначить каждую точку из моего второго шейп-файла (shp2
) к одному из треугольников Делоне, удерживаемых внутри объекта vtess
,
vtess <- deldir(shp1$coords.x1,shp1$coords.x2)
Я могу назначить каждого shp2
укажите на мозаику мозаики следующим образом:
tl <- tile.list(vtess)
shp2$tile <- 99
for (i in 1:length(shp2$tile)){
shp2$tile[i] <- which.tile(shp2$coords.x1[i], shp2$coords.x2[i], tl)
}
Тем не менее, пока я могу найти список треугольников, используя triang.list(vtess)
, здесь нет which.triang
команда, поэтому невозможно назначить аналогичным образом shp2
указывает на треугольники.
Я предполагаю, что мне нужно (1) создать SpatialPolygon из треугольников, а затем (2) перекрыть shp2
указывает на новый файл многоугольника. Это самый эффективный способ? И если да, то как мне создать шейп-файл многоугольника из треугольников, содержащихся внутри vtess
?