Raw Message

Message-ID: <CAJgdCD4rMt1A9CoZqhfMyJp2Kidwn4Fqja6KuqXnjysPcbkAtw@mail.gmail.com>
Date: 2016-02-18T09:25:02Z
From: John Wasige
Subject: Intersect shapefiles

?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]]