Как установить максимальную длину стороны треугольника в триангуляции Делоне?

Как удалить расстояния из триангуляции Делоне, которые больше, чем мне нужно?

Пример данных:

x<-rep(1:12, c(2,2,7,9,10,5,4,6,10,10,9,4))
y<-c(1,2,1,2,1:3,5:8,1:9,1:10,2,7:10,8:11,7:12,3:12,3:12,4:12,5,8:10)
x_plus<-seq(0.2:0.8, by=0.1)
x<-x+sample(x_plus, 78, replace=TRUE)
y<-y+sample(x_plus, 78, replace=TRUE)

Построение карты:

plot(x,y)

Триангуляция Делоне с tri.mesh() - package(tripack)

my.triangles<-tri.mesh(x,y)
plot(my.triangles, do.points=FALSE, lwd=0.2)
points(x,y, col = "blue", pch=20)

Как я могу извлечь только более короткие расстояния? Мне не нужны те большие, вы наверняка знаете, какие расстояния я имею в виду. Есть ли какой-то аргумент для этого в функции tri.mesh()? Или это можно было сделать после него?

Сохраняются ли расстояния в этом объекте?

my.triangles

triangulation nodes with neigbours:
node: (x,y): neighbours
1: (1.4,1.7) [5]: 2 3 4 11 12
2: (2,3) [6]: 1 4 7 8 9 11
3: (3,1.8) [4]: 1 4 5 12
.
.
.
76: (12.4,8.8) [5]: 68 69 70 75 77
77: (12.9,9.9) [6]: 70 71 72 75 76 78
78: (13,11) [4]: 72 73 74 77
number of nodes: 78
number of arcs: 221
number of boundary nodes: 10
boundary nodes:  1 11 12 45 56 66 74 75 77 78
number of triangles: 144
number of constraints: 0

r extract delaunay
person Ladislav Naďo    schedule 16.07.2013    source источник
comment
Не могли бы вы сделать ваш пример воспроизводимым? Пакет tripolt недоступен в CRAN.   -  person QuantIbex    schedule 16.07.2013
comment
извините, это посылка: tripack   -  person Ladislav Naďo    schedule 16.07.2013
comment
Каковы конкретные критерии для удаления треугольника/ребра?   -  person QuantIbex    schedule 16.07.2013

Ответы (1)


arrow_upward
6
arrow_downward

Документация (?tri) предполагает, что сегменты имеют форму r$tlist[k] -- r$tlist[r$tlptr[k]]: вы можете удалить те, которые находятся за пределами некоторого порога.

r <- tri.mesh(x,y)
k <- seq_len( r$tlnew - 1 )
i <- r$tlist[k]
j <- r$tlist[r$tlptr[k]]
keep <- i > 0
i <- abs( i[ keep ] )
j <- abs( j[ keep ] )
plot( x, y )
segments( r$x[i], r$y[i], r$x[j], r$y[j], col="grey" )
distances <- sqrt( ( r$x[i] - r$x[j] ) ^ 2 + ( r$y[i] - r$y[j] ) ^ 2 )
threshold <- 2.5  # Choose the threshold manually
i <- i[ distances < threshold ]
j <- j[ distances < threshold ]
segments( r$x[i], r$y[i], r$x[j], r$y[j], lwd = 2 )

меньшая триангуляция

person Vincent Zoonekynd    schedule 16.07.2013