Dear list:



I am trying to calculate the intersection area between two irregular
polygons (see example data below).



set.seed(1234)

theta <- seq(0, 2 * pi, length=(100))

poly1<-  cbind(c(0 + 1 * cos(theta) + rnorm(100, sd=0.1)), c(0 + 2 *
sin(theta)))

poly2<-  cbind(c(0 + 2 * cos(theta) ), c(-1 + 1.5 * sin(theta)+ rnorm(100,
sd=0.1)))



plot(x, y, type = "n", , xlim=c(-5,5), ylim=c(-5,5))

polygon(poly1)

polygon(poly2)





My data correspond to geographical positions, therefore I cannot generate it
using equations as I did in the example.

Are there any straightforward algorithm to calculate the intersection area
between polygons poly1 and poly2?



Thank you so much in advance



Jonas

	[[alternative HTML version deleted]]

