Упростить полигоны объекта sf

Как мне упростить sfмногоугольник без введения пробелов и осколков?

С шейп-файлом, например, я бы использовал 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

Я пробовал, sf::st_cast()который, из man-страниц, заявляет:

Приведите геометрию к другому типу: либо упростите, либо приведите явно

и:

для аргументации: характер; тип цели, если отсутствует, попытка упрощения; когда x имеет тип sfg (т. е. одиночную геометрию), то необходимо указать.

Когда я toпропал без вести, это не сработало, как ожидалось (я знал, что это слишком хорошо, чтобы быть правдой!):

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

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

В настоящее время я открываю файл rgdal::readOGR(), упрощаю его, сохраняю и загружаю снова sf.

Есть ли способ лучше?

rgeos::gSimplify()

Предложение @sk rgeos::gSimplify()может делать топологически упрощенные упрощения (то есть упрощает без создания осколков), если указано со следующими аргументами:

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

gSimplify@dataоднако, не сохраняет фрейм, поэтому мы должны создать его заново:

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

И это действительно приводит к меньшему размеру файла (можно настроить tolаргумент, чтобы уменьшить его), и я подтвердил, что это не создало проскальзываний, изучив его в QGIS.

object_size(regions_gSimplify)
# ~8MB

Так что, хотя это верная альтернатива, у rmapshaper::ms_simplify()меня все еще есть та же проблема, а именно, что она не работает с 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"

Ответ @obrl_soil также может быть применен gSimplify(), просто используйте его вместо ms_simplify().

Вы можете привести объект sf к sp для пакетов, которые еще не поддерживают sf - я делаю это справедливо для взаимодействий растров / полигонов. Так что вы могли бы сделать:

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