Hi,

Try this,

d <- data.frame(x=runif(1e3, 0, 30), y=runif(1e3, 0, 30))
d$z = (d$x - 15)^2 + (d$y - 15)^2

library(spatstat)
library(maptools)

W <- ripras(df, shape="rectangle")
W <- owin(c(0, 30), c(0, 30))
X <- as.ppp(d, W=W)
Y <- dirichlet(X)
Z <- as(Y, "SpatialPolygons")
plot(Z, col=grey(d$z/max(d$z)))

and also panel.voronoi in latticeExtra.

Unfortunately I do not know of a solution that uses more efficient
algorithms for computing the Dirichlet tessellation and extracting
tiles than those relying on deldir. The Triangle package (r-forge)
looks promising.

HTH,

baptiste


On 5 September 2011 06:26, awesolow <aweso...@andrew.cmu.edu> wrote:
> Hi,
>
> I have a set of x, y points (longitude/latitude) along with a z value
> representing an attribute at each point.  I want to create a
> Voronoi/Dirichlet tesselation of these points coloring each tile by the z
> value.  I tried searching for a way to solve this and it was suggested to
> use the dirichlet() command to get the correct coloring.  However, my
> coloring is not correct.
>
> Any thoughts?
>
> Thanks in advance.
>
> --
> View this message in context: 
> http://r.789695.n4.nabble.com/Coloring-Dirichlet-Tiles-tp3789746p3789746.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to