Simplificar polígonos do objeto sf

Como simplifico um sfpolígono sem introduzir espaços e lascas?

Com um shapefile, por exemplo, eu usaria rmapshaper::ms_simplify():

library("pryr")
library("rgdal")
library("rmapshaper")

download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip",
              destfile = "regions.zip")
unzip("regions.zip")
regions <- readOGR(".", "england_gor_2011")
object_size(regions)
# ~13MB

regions <- ms_simplify(regions)
object_size(regions)
# < 1MB

Eu tentei o sf::st_cast()que, nas páginas de manual, afirma:

Transmitir geometria para outro tipo: simplificar ou converter explicitamente

e:

argumentar: personagem; tipo de destino, se ausente, a simplificação é tentada; quando x é do tipo sfg (ou seja, uma única geometria), deve ser especificado.

Quando deixei tocomo falta, isso não funcionou como esperado (eu sabia que era bom demais para ser verdade!):

library("sf")
regions <- sf::read_sf("england_gor_2011.shp")
object_size(regions)
# ~13MB

regions <- sf::st_cast(regions)
object_size(regions)
# Still 13MB

Atualmente, estou abrindo o arquivo rgdal::readOGR(), simplificando-o, salvando-o e carregando-o novamente com sf.

Existe uma maneira melhor?

rgeos::gSimplify()

A sugestão de @sk rgeos::gSimplify()pode fazer simplificações com reconhecimento topológico (isto é, simplifica sem criar lascas) quando especificada com os seguintes argumentos:

library("rgeos")
regions_gSimplify <- gSimplify(regions, tol = 0.05, topologyPreserve = TRUE)

gSimplifynão preserva o @dataquadro, portanto, devemos recriá-lo:

regions_df <- regions@data
regions_gSimplify <- sp::SpatialPolygonsDataFrame(regions_gSimplify, regions_df)

E isso realmente resulta em um tamanho de arquivo menor (pode alterar o tolargumento para diminuí-lo) e eu confirmei que isso não havia criado lascas ao examiná-lo no QGIS.

object_size(regions_gSimplify)
# ~8MB

Portanto, embora essa seja uma alternativa válida, rmapshaper::ms_simplify()ainda tenho o mesmo problema, a saber, que não funciona com sf:

regions_sf <- sf::read_sf("england_gor_2011.shp")
object_size(regions_sf)

regions_gSimplify <- gSimplify(regions_sf, topologyPreserve = TRUE, tol = 0.05)
# Error in gSimplify(regions_sf, topologyPreserve = TRUE, tol = 0.05) : 
# no slot of name "proj4string" for this object of class "sf"

A resposta do @obrl_soil também pode ser aplicada gSimplify(), basta usá-la no lugar de ms_simplify().

Você pode converter um objeto sf em sp, para pacotes que ainda não suportam sf - eu faço isso bastante para interações raster / polígono. Então você poderia fazer:

simplepolys <- rmapshaper::ms_simplify(input = as(sfobj, 'Spatial')) %>%
  st_as_sf()