Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

How to calculate the area of polygon overlap in R?

Does anyone know how to calculate the area in common between 2 or more polygons in R? I would like to have the output of such a calculation be the coordinates of a new polygon for that area of overlap. Cheers

like image 651
Marc in the box Avatar asked Apr 02 '12 13:04

Marc in the box


1 Answers

EDIT: these days I would use the 'intersect', 'cover', 'erase', 'union' and related functions in the 'raster' package. They do the hard work to keep the top-level object and attributes.

ORIG: You could use the rgeos package with its gIntersection function. Successive calls between pairs and resulting intersections will get you there. See

library(rgeos)
?gIntersection

You will need to get into the structure of "SpatialPolygons" in the sp package to get the final coordinates. See the vignette("sp").

like image 75
mdsumner Avatar answered Oct 20 '22 17:10

mdsumner