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.
