Saya ingin menggeneralisasi SpatialPolygonsDataFrame
dalam R untuk membuat plot lebih efisien. Poligon saya menghasilkan plot pdf besar yang lambat dimuat dan sulit diintegrasikan ke dalam dokumen.
Saya mencoba menggunakan gSimplify
fungsi dari rgeos
paket, tetapi sayangnya itu tidak menghemat topologi poligon.
Apakah ada alternatif yang bisa saya gunakan?
r
topology
generalization
yellowcap
sumber
sumber
gSimplify()
memiliki parameter opsionaltopologyPreserve
yang standarnya adalahFALSE
. Sudahkah Anda mencoba mengaturnyaTRUE
?Jawaban:
Ada diskusi tentang ini di r-sig-geo . Untuk jawaban yang pasti Anda harus bertanya di sana, karena ada orang yang tahu wawasan spasial R.
Tapi, Anda juga bisa melakukan ini di aplikasi desktop GIS (ekspor bentuk menggunakan perintah writeOGR dari rgdal atau writePolyShape () dari maptools) seperti QuantumGIS, GRASS atau SAGA.
Untuk QuantumGIS gunakan Vector / Geometry Tools / Sederhanakan geometri (Saya telah menguji dan tidak mempertahankan topologi, tetapi diterapkan ke Rumania admin1 shapefile, terlihat baik dengan 2000 node).
Untuk penggunaan GRASS v.generalisasi (baca manual untuk info tentang algoritme, ada beberapa).
Untuk SAGA, Anda harus terlebih dahulu mengubah poligon menjadi garis (Bentuk - Garis / Konversi Poligon ke Garis), kemudian menyederhanakan garis (Bentuk Garis / Penyederhanaan Garis), dan akhirnya mengubah garis menjadi poligon (Bentuk - Poligon / Mengubah Garis ke Poligon). Saya telah menguji ini dan hasilnya memiliki topologi yang dipertahankan.
Ceria dan semoga sukses,
sumber
Sekarang ada
rmapshaper
paket luar biasa dari Andy Teucher yang mencakup fungsi penyederhanaan yang "melakukan penyederhanaan multi-poligon yang disadari secara topologi".Dari repo github mereka, contoh penggunaan:
dimana
states_sp
sebuahspatialPolygons*
benda.Lihat paket README untuk informasi lebih lanjut: https://github.com/ateucher/rmapshaper
sumber
itu akan membantu:
gSimplify(..., topologyPreserve=T)
dariGEOS
paketsumber
gSimplify
tidak benar-benar mempertahankan topologi selama beberapa poligon. Saya pikir ini lebih berhubungan dengan topologi single atau multipolygons. Untuk shapefile dengan banyak poligon, topologi tidak dilestarikan.