Hi,

I'm using the package 'tripack' to look at the Voronoi tesselations of point patterns. I want to calculate the area of the intersection of each polygon with the unit square (in which all of my points lie). So far I have been using:

voronoi.area(voronoi.mosaic(runif(10),runif(10)))

This returns NA for edge polygons which are not bounded. However, their intersection with some bounded area is bounded. So, I'd like to know how to get the area of the intersection of each polygon with some defined area (ie. the unit square). Is this possible?

Cheers,

Jeremy
-- 
j.silver@ugrad.unimelb.edu.au


