Bagaimana cara menyederhanakan sf
poligon tanpa memasukkan celah dan sliver?
Dengan shapefile, misalnya, saya akan menggunakan 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
Saya sudah mencoba sf::st_cast()
yang, dari halaman manual, menyatakan:
Transmisikan geometri ke tipe lain: sederhanakan, atau gulirkan secara eksplisit
dan:
argumen: karakter; tipe target, jika hilang, penyederhanaan dicoba; ketika x adalah tipe sfg (yaitu, satu geometri) maka perlu ditentukan.
Namun ketika saya pergi to
sebagai hilang, ini tidak berfungsi seperti yang diharapkan (saya tahu itu terlalu bagus untuk menjadi kenyataan!):
library("sf")
regions <- sf::read_sf("england_gor_2011.shp")
object_size(regions)
# ~13MB
regions <- sf::st_cast(regions)
object_size(regions)
# Still 13MB
Saat ini saya sedang membuka file dengan rgdal::readOGR()
, menyederhanakannya, menyimpan ini, kemudian memuatnya lagi dengan sf
.
Apakah ada cara yang lebih baik?
rgeos::gSimplify()
Saran @sk tentang rgeos::gSimplify()
dapat melakukan penyederhanaan sadar-topologi (yaitu menyederhanakan tanpa membuat sliver) ketika ditentukan dengan argumen berikut:
library("rgeos")
regions_gSimplify <- gSimplify(regions, tol = 0.05, topologyPreserve = TRUE)
gSimplify
tidak menyimpan @data
bingkai, jadi kita harus membuatnya kembali:
regions_df <- regions@data
regions_gSimplify <- sp::SpatialPolygonsDataFrame(regions_gSimplify, regions_df)
Dan ini memang menghasilkan ukuran file yang lebih kecil (dapat mengubah tol
argumen untuk membuatnya lebih kecil) dan saya mengkonfirmasi ini belum membuat sliver dengan memeriksanya di QGIS.
object_size(regions_gSimplify)
# ~8MB
Jadi, meskipun ini merupakan alternatif yang valid untuk rmapshaper::ms_simplify()
saya masih memiliki masalah yang sama, yaitu bahwa itu tidak bekerja dengan 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"
Jawaban @obrl_soil juga dapat diterapkan gSimplify()
, cukup gunakan sebagai ganti ms_simplify()
.
st_simplify
seharusnya melakukan itu? (belum menggunakannya)
st_simplify
, terima kasih telah menunjukkannya. Saya lebih suka algoritma yang rmapshaper::ms_simplify
defaultnya melebihi semua yang lain yang saya coba sejauh ini, tapi saya akan bermain dengan opsi baru (pembaruan: whoa lanjutkan dengan hati-hati, preserveTopology = TRUE
pasti belum bekerja dengan benar)
regions
) tetapi di luar itu tidak lagi mempertahankan topologi. Ketika itu pecah pada titik tertentu saya akan mengatakan itu bukan perilaku yang diinginkan