library(orsifronts)
## Loading required package: sp
library(raster)
library(rgeos)
## rgeos version: 0.3-23, (SVN revision 546)
## GEOS runtime version: 3.5.1-CAPI-1.9.1 r4246
## Linking to sp version: 1.2-4
## Polygon checking: TRUE
ex <- extent(10, 130, -80, -30)
p <- rgeos::gPolygonize(rgeos::gNode(rgeos::gUnion(raster::intersect(orsifronts, as(ex, "SpatialPolygons")), as(ex, "SpatialLines"))))
## Warning in raster::intersect(orsifronts, as(ex, "SpatialPolygons")): non
## identical CRS
## Warning in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td,
## unaryUnion_if_byid_false, : spgeom1 and spgeom2 have different proj4
## strings
d <- data.frame(region = seq(length(p)))
pd <- SpatialPolygonsDataFrame(p, d, match.ID = FALSE)
plot(pd, col = rainbow(nrow(pd)), asp = "")
plot(pd, col = grey(seq(0.1, 0.9, length = nrow(pd))), asp = "")
text(coordinates(pd), lab = pd$region, col = rainbow(nrow(pd)))
extract(pd, cbind(c(90, 110, 120), c(-70, -60, -30)))
## point.ID poly.ID region
## 1 1 4 4
## 2 2 5 5
## 3 3 2 2