Hi Thank you for your useful info.
However for following statement, > voronoiShapefile <- function() > { -----X8--------------X8-----------X8-- > for (i in seq(along=polys)) > { > pcrds <- > cbind(tritest.vp[[8]][,1],tritest.vp[[8]][,2] ) I think all "8"s needed to "i" like this. cbind(tritest.vp[[i]][,1],tritest.vp[[i]][,2] ) BTW, any method finding neightbour voronoi'cells' IDs? Regards.
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo