# NOT RUN {
run <- FALSE
if (require(maptools)) run <- TRUE
if (run) {
nc1 <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
proj4string=CRS("+proj=longlat +datum=NAD27"))
}
if (run) {
print(system.time(oU <- gUnionCascaded(nc1)))
if (version_GEOS0() >= "3.8.0") {
print(system.time(oCU <- gCoverageUnion(nc1)))
}
}
# }
Run the code above in your browser using DataLab