Hi Ravi,
After performing 'gIntersects' try merging this result to your dataframe,
then select the TRUE values?
Jesse
cz_zip <- gIntersects(camapzip_temp,camap_base, byid=TRUE)
camap_base@data<-cbind(camap_base@data, czip)
camap2<-camap_base[(camap_base@data[,/???/] %in% c("TRUE)),] # Insert your
column number for czip in place of ???
plot(camap2)
head(camap2@data)
-----
--------------------------------------------------------
Jesse D Berman, PhD
Yale University
School of Forestry and Environmental Studies
Post-Doc Fellow
--
View this message in context:
http://r-sig-geo.2731867.n2.nabble.com/Clipping-a-Map-tp7585156p7585162.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
