working configuration:
R-3.1.1
WIN-64
goal: want to perfrom ploygon clipping, giving bounding box values of my
study area and extract the area from world map.
I followed:
> clip.extent <- as(extent(76.3700, 31.7439, 78.6541, 33.2653),
> "SpatialPolygons")
> proj4string(clip.extent) <- CRS(proj4string(WorldMap))> LSMap <-
> gIntersection(WorldMap, clip.extent, byid = TRUE)
Error :
In RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_not_poly,
"rgeos_intersection") :
spgeom1 and spgeom2 have different proj4 strings
Any help or assistance would be highly obliged
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
