[R-sig-Geo] Area left in shape file

Barry Rowlingson b.rowlingson at lancaster.ac.uk
Tue Sep 11 00:46:06 CEST 2012


On Mon, Sep 10, 2012 at 11:31 PM, Juan Tomas Sayago
<juantomas.sayago at gmail.com> wrote:
> Dear list,
> I have two shapes with polygons, and need to find the area in one of them
> that is not covered by the other. What do you recommend?

 You can do operations with polygons using the rgeos package - its
probably as simple as something like:

 A = readOGR(".","shape1.shp")
 B = readOGR(".","shape2.shp")

 C = gDifference(A,B)

then get the area of C.

Barry



More information about the R-sig-Geo mailing list