similar to: Voronoi

Okay, been working with tripack, seems the most mature package for this. Got it to work well with their test data set - data(tritest). When i tried random numbers to explore further, i am getting some results that don't reconcile. example run this: library(tripack) y <- runif(100) x <- runif(100) vm <- voronoi.mosaic(x,y) plot(vm) par(new=T) plot(x,y,col='blue') when

There's Renka's STRIPACK, and TRIPACK, respectively, ACM TOMS Algorithms 772 and 751, and there's the R package "deldir" which does the Delaunay for a plane, but does anyone have or know of the tessellation in R for a sphere? Also, is there a standard indexing scheme for Delaunay facets, and perhaps of edges in such facets? I'd expect that to be a publication reference,

Dear R-users, Does anyone know how to do weighted Voronoi diagrams (Dirichlet tesselation) in R? To be more specific, I have a set of coordinates for tree locations on a plot, and I'm looking for a way to do the tesselation so that the polygon size for each tree depends on the size of the subject tree, and the size of its neighbors. So, the location of the bisection between two trees

Hi all, I am running R 2.5.0 under Windows XP Media Center Edition. Here's a problem that's been stumping me for a few days now, and I can't find anything useful in the archives. I am using voronoi.mosaic (tripack package) to create proximity polygons for a study of vegetation competition and dynamics. The points lists are read in from a file for each plot, then 8 duplicates

Hello list, I'm just getting started with using R - I have been trying over the past day or so to work out a method for generating voronoi polygons for PostGIS using SQL. I was able to put together a procedure which works relatively well, but is somewhat inefficient. Someone on the PostGIS list pointed me to the deldir() function in R, for which I can export a text file with x/y coordinates

Hi, I've been looking for functions that can do delaunay tesselation and generate voronoi cells. I came across deldir and tripack but both seem to be restricted to 2D points. Are there any packages that can do a tesselation in N dimensions? I know that Matlab and Mathematica use the qhull package to provide functions for this. Does anybody know of any R packages that do this (maybe by calling

Hello. I'd like to convert voronois object created by tripack to polygons to use them in GIS(Geographic Information Systems) software. I tried to create voronoi objects by using following code. > library(tripack) > x <- rnorm(10) > y <- rnorm(10) > plot(x,y) > v<-voronoi.mosaic(x,y) > plot(v) But from here, I could not create polygons. Of course, I know

Is there any pure r code to do delaunay or voronoi diagrams? Thanks! --------------------------------- [[alternative HTML version deleted]]

I would like to create a voronoi map that can show me point values like the ones in colum street_cle. I was able to do it in tileplot(latticeExtra) like the ones attached, and with the following commands: xy=list(x=x2$X, y=x2$Y) xy$street_cle <- x2$street_cle tileplot(street_cle~ x*y, xy) But this is a trellis. Is it possible to be converted to a shapefile or image so I can plot it using

Hi, I have a monthly time series with missing values and I would use stl function to identify seasonality. I tried all settings of na.action but the result is the same: stl(tm245,s.window=11, na.action=na.pass) Error in stl(tm245, s.window = 11, na.action = na.pass) : NA/NaN/Inf in foreign function call (arg 1) Can you help me? Thanks Andrea Toreti [[alternative HTML version

Hi all, Given a number of points (x,y) in a plane, I'd like to plot a map of polygons, so that 1) each polygon contains exactly one point 2) the polygon defines the area for which this specific point is closer than any other point. It's a bit like a map of areas "influenced" by that point, and it's obviously a matter of intersecting the perpendicular bisectors

There seems to be an error in the X11 display code: try: plot(c(-30000,-1,0,1,30000),c(30000,1,0,1,30000),type="l") points(c(-30000,-1,0,1,30000),c(30000,1,0,1,30000)) it gives a big "V" as expexted. Now zoom in using small xlim and ylim: plot(c(-30000,-1,0,1,30000),c(30000,1,0,1,30000),type="l", xlim=c(-5,5),ylim=c(0,5))

hi all members, please, i need you help... now a i´m working with veronoi polygons in a area with projections, but i need cut the polygons left. On other words, i need cut the polygons in the worked area. R help say that use the command voronoi.findrejectsites, but in this command i need put the numbers, any way...this command not cut!! do you can help me? Thank you for help me! José Bustos

Dear Help List, My name is Matt Oliver. I have been using R for about a year and find it very helpful. However, I have a need for a function that I cannot find. I am not very good at programming so I thought I would ask the group. I have an irregular grid of data (x = Longitude, y = Latitude). Each pair of my x,y has a categorical value. Obviously linear or any other numerically based

hello, I want to do a test for spatial correlation. I tried it with geary.test() but I don't understand the required input. x= a numeric vector the same length as the neighbours list in listw (my sampled data, I assume) listw= a listw object created for example by nb2listw (well when I check nb2listw() I get to "neighbours - an object of class nb" - but I couldn't figure

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. --27464147-1968795515-1164715405=:5670 Content-Type: TEXT/PLAIN; charset=iso-8859-1; format=flowed Content-Transfer-Encoding: 8BIT You need to re-install packages when you update R, as it says in the rw-FAQ. Use update.packages(checkBuilt=TRUE) to

Full_Name: Larry Layne Version: 2.4.0 OS: MS Windows 2000 Submission from: (NULL) (129.24.91.249) All of the following libraries loaded just fine in version 2.3.1: > utils:::menuInstallLocal() package 'tripack' successfully unpacked and MD5 sums checked updating HTML package descriptions > library(tripack) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable

Full_Name: Michael R. Allen Version: 1.01 OS: linux redhat Submission from: (NULL) (149.149.32.46) Hello, Thought you might be interested. I tried to load the package tripack needed by the sgeostat package. After some careful searching, three packages were found which had the same error. Here is the output. Mike > library(acepack) Error in dyn.load(x, as.logical(local),

Hello, I've been debugging a few different tessellation shader issues with nouveau, but let's start small. I see this issue on my GK208 with high frequency, and I *think* I've seen it once or twice on my GF108, but it's exceedingly rare, if it does happen. I don't have a GK10x to test on, unfortunately, but I assume it'll have the same issue as the GK208. The issue is

One additional observation that I just made is that on GK208, the blob apparently doesn't use the result of S2R Rx, SR_INVOCATION_ID wholesale in TCS. It either passes it through a I2I.S32.S32 Rx, |Rx| (i.e. absolute value), or even more paradoxically, shl 2; shr 2; which removes the top *2* bits, rather than just the top 1. However I see no such behaviour on GF108. I'm going to test out