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()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를 지원하지 않는 패키지의 경우 sf 객체를 sp로 캐스트 할 수 있습니다. 래스터 / 폴리곤 상호 작용을위한 적절한 비트입니다. 그래서 당신은 할 수 있습니다 :

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