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]] ______________________________________________ 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.