Saat gadm
memetakan peta Rusia R
, saya memiliki masalah pada meridian 180 derajat, di mana bagian wilayah Rusia diplot secara terpisah. Anda dapat melihatnya di plot menggunakan kode berikut (itu tepi kiri dan bagian ada di kanan). Adakah yang bisa menyarankan bagaimana menghindari masalah ini?
Ini adalah versi singkat dari kode untuk pertanyaan saya:
require(sp)
rus<-url("http://www.gadm.org/data/rda/RUS_adm1.RData")
print(load(rus))
gadm$regions = as.factor(1:88)
spplot(gadm,"regions")
sumber
sp
menyebutkan bahwa kelas ini mendukung "metode termasuk cetak / tampilkan, plot, subset, [, [[, \ $, nama, redup, ringkasan, ..." Halaman bantuan untuksubset
menunjukkan sintaksisnya. Kuncinya adalah mengetahui bidang mana yang akan digunakan untuk menemukan Moskow! Saya menemukannya dengan melihat nilai-nilai unik dari semua bidang, akhirnya turunlevels(gadm$NAME_1)
. Ini mengarah kespplot(subset(gadm.prj, subset=(NAME_1=="Moskva")),"regions")
, yang melakukan trik.Saya melihat bahwa whuber sudah mengalahkan saya untuk itu, tetapi saya akan memposting jawaban saya.
Jika Anda ingin menyimpan data dalam proyeksi longlat, Anda dapat menambahkan parameter " Pembungkus Longitudinal " ke string Proj4.
Ini akan memusatkan peta Anda di garis bujur 105 derajat timur.
sumber
spTransform
.