简化SF对象的多边形

如何在sf不引入间隙和裂缝的情况下简化多边形？

以shapefile为例，我将使用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()手册页中指出的内容：

将几何体转换为另一种类型：简化或显式转换

和：

争论：性格 目标类型，如果缺少，则尝试简化；当x为sfg类型（即单个几何）时，则需要指定to。

当我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()

rgeos::gSimplify()当使用以下参数指定时，@sk的建议可以进行拓扑感知的简化（即在不创建条的情况下进行简化）：

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的程序包，可以将sf对象转换为sp-对于栅格/多边形交互，我可以这样做。因此，您可以执行以下操作：

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