Check out also function joinPolys in the PBSmapping package (http://www.inside-r.org/packages/cran/pbsmapping/docs/joinPolys)
- - - A. Márcia BARBOSA, post-doctoral fellow CIBIO/InBIO - University of Évora (Portugal) http://modtools.wordpress.com/barbosa 2016-02-18 10:04 GMT+00:00 Loïc Dutrieux <loic.dutri...@wur.nl>: > Hi, > > Have you tried gIntersection from the rgeos package. > > Given 3 SpatialPolygon* objects you can probably do something like: > > sp12 <- gIntersection(sp1, sp2) > sp123 <- gIntersection(sp12, sp3) > > Cheers, > Loïc > > > On 02/18/2016 10:25 AM, John Wasige wrote: > >> Dear all, >> >> I have three shapefiles that I want to intersect. Any body with an idea on >> how to do the intersection? Shapefiles below; >> library(rgdal) >> library(rgeos) >> library(maptools) >> library(sp) >> ###shape1 >> Rain_poly_merge <- readOGR(dsn="G:/EPIC_CUP/Rain_poly_merge.shp", >> layer="Rain_poly_merge") >> Rain_poly_merge <- spTransform(Rain_poly_merge, CRS("+init=epsg:4326")) >> >> ###shape2 >> Rain_poly_merge <- readOGR(dsn="G:/EPIC_CUP/yield_maize_average2002.shp", >> layer="Rain_poly_merge") >> Rain_poly_merge <- spTransform(Rain_poly_merge, CRS("+init=epsg:4326")) >> >> ###shape3 >> Rain_poly_merge <- readOGR(dsn="G:/EPIC_CUP/yield_maize_average2008.shp", >> layer="Rain_poly_merge") >> Rain_poly_merge <- spTransform(Rain_poly_merge, CRS("+init=epsg:4326")) >> >> [[alternative HTML version deleted]] >> >> _______________________________________________ >> R-sig-Geo mailing list >> R-sig-Geo@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >> >> > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo