如何在R中設置Delaunay三角剖分中的三角形邊的最大長度？

如何去除比我需要的Delaunay三角剖分的距離？如何在R中設置Delaunay三角剖分中的三角形邊的最大長度？

實施例的數據：

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)

Delaunay三角與tri.mesh（） - 包（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

來源

2013-07-16 Ladislav Naďo

能否請您讓您的例子可以重現？看起來''tripolt'看起來並不適合CRAN。 – QuantIbex

對不起，它是包裝：三件套 –

去除三角形/邊緣的具體標準是什麼？ – QuantIbex

的文檔（?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)

smaller triangulation

來源

2013-07-16 11:08:21

如何在R中設置Delaunay三角剖分中的三角形邊的最大長度？

回答

相關問題