Ok, now I'm scared. I copied the exact code from your post and got
length(tri) = 1481 and
dim(geodel) = 15773
You are right though, there seem to bee duplicates in my original bm,
but removing them still got diffeent results. I fear there's something
weird going on with my R installation.
Che
I just tried the following:
require(deldir)
require(geometry)
set.seed(42)
bm <- data.frame(x=sample(1:82,800,TRUE),y=sample(1:82,800,TRUE))
del <- deldir(bm)
tri <- triang.list(del)
geodel <- delaunayn(bm)
length(tri) # Got [1] 1481
dim(geodel) # Got [1] 14813
So all seems to be in harmony
Hello again,
I have found further depths of confusion concerning delaunay
triangulations to explore.
Here is the code I'm using to create the confusing results:
bm <- getbm(x) # a data.frame with 2 columns and 800 rows, values are
integers ranging from 1 to 82
del <- deldir(bm) # creating an obj
On 21/05/14 23:34, Raphael Päbst wrote:
I believe you are right. A night of sleep has done wonders for my
understanding of the problem.
Oh sleep it is a gentle thing
Beloved from pole to pole!
Thank you for your patience and help!
Patience? Moi? This must be some new use of the word "pati
I believe you are right. A night of sleep has done wonders for my
understanding of the problem.
Thank you for your patience and help!
Raphael
On 5/21/14, Rolf Turner wrote:
>
>
> It sounds to me as though you are simply getting yourself flummoxed by
> the fact that the different packages produc
It sounds to me as though you are simply getting yourself flummoxed by
the fact that the different packages produce their output in different
formats. The information in the output will be the same (as Boris has
indicated) --- it will just be arranged differently. Learn to interpret
the ou
Thanks for the answer!
I'll post a sample tomorrow, I have however found the following:
triang.list() gives me the coordinates of the triangle's vertices,
while delaunayn() gives me the indices of those coordinates. Thus it
should be more or less simple to convert the output of deldir() into
that
deldir() uses Lee and Schacter's algorithm, while the geometry package is a
(partial) implementation of Barber et. al's Quickhull. Since both algorithms
are correct, they should give the same results for the same data.
How about you post a small input dataset and list the output that you need...
Thank you very much, this looks promising.
I have a follow-up question however, probably due to my thickness when
it comes to the underlying math.
I am translating (as closely as possible) some code that has
originally been written for Mathlab and uses the delaunay() function
there. Now, if I und
Try:
install.packages("deldir")
library(deldir)
?deldir
set.seed(16180)
x <- runif(20); y <- runif(20); window <- c(0,1,0,1)
tess <- deldir(x, y, rw = window)
plot.deldir(tess, wpoints="real", wlines="tess")
Cheers,
Boris
On 2014-05-14, at 8:18 AM, Raphael Päbst wrote:
> Hello everyone!
> I ha
Hello everyone!
I have returned to R after a longish break and am currently working on
a project where I need Delaunay-Triangulations and Voronoi-Diagrams.
If I understood it correctly, the Geometry-Package only offers
functions for the Delaunay-Triangulation at the moment. Is this
correct and if s
11 matches
Mail list logo
