simplify_compare.R

library(sf)

## Linking to GEOS 3.5.1, GDAL 2.1.2, proj.4 4.9.3

nc <- read_sf(system.file("shape/nc.shp", package="sf"), quiet = TRUE)
plot(st_geometry(st_simplify(nc, preserveTopology = TRUE, dTolerance = .07)), col = "grey")

## Warning in st_simplify.sfc(st_geometry(x), preserveTopology, dTolerance):
## st_simplify does not correctly simplify longitude/latitude data, dTolerance
## needs to be in decimal degrees

plot(st_geometry(nc), add = TRUE, col = "transparent", border = "firebrick")

library(rmapshaper)
library(sp)
plot(rmapshaper::ms_simplify(as(nc, "Spatial"), keep_shapes = TRUE, keep = 0.01), col = "grey")
plot(st_geometry(nc), add = TRUE, col = "transparent", border = "firebrick")