Совпадение пространственных точек с треугольниками Делоне

У меня есть 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?

0 ответов

Другие вопросы по тегам