## Not run:
# if (require(maps)) {
# world <- map("world", fill=TRUE, col="transparent", plot=FALSE)
# worldSpP <- map2SpatialPolygons(world, world$names, CRS("+proj=longlat +ellps=WGS84"))
# worldSpP <- worldSpP[-grep("Antarctica", row.names(worldSpP)),]
# # incomplete polygons
# worldSpP <- worldSpP[-grep("Ghana", row.names(worldSpP)),]
# # self-intersection mouth of Volta
# worldSpP <- worldSpP[-grep("UK:Great Britain", row.names(worldSpP)),]
# # self-intersection Humber estuary
# worldSpPr <- recenter(worldSpP)
# plot(worldSpPr)
# title("Pacific view without polygon splitting")
# worldSpPnr <- nowrapRecenter(worldSpP)
# plot(worldSpPnr)
# title("Pacific view with polygon splitting")
# }
# ## End(Not run)
if (rgeosStatus()) {
crds <- matrix(c(-1, 1, 1, -1, 50, 50, 52, 52), ncol=2)
rcrds <- rbind(crds, crds[1,])
SR <- SpatialPolygons(list(Polygons(list(Polygon(rcrds)), ID="r1")),
proj4string=CRS("+proj=longlat +ellps=WGS84"))
bbox(SR)
SRr <- recenter(SR)
bbox(SRr)
SRnr <- nowrapRecenter(SR)
bbox(SRnr)
}
Run the code above in your browser using DataCamp Workspace