Ao traçar o gadmmapa da Rússia R, tenho um problema no meridiano de 180 graus, onde parte da região russa é plotada separadamente. Você pode vê-lo no gráfico usando o código a seguir (é a borda esquerda e a parte está à direita). Alguém pode sugerir como evitar esse problema?
Aqui está uma versão curta do código para minha pergunta:
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")
Respostas:
Projete os dados. Este exemplo usa uma projeção estereográfica polar centralizada perto do meio da Rússia (longitude 105 graus leste):
require(rgdal)
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "regions")
A região do problema está no canto superior direito. Você pode ver que essa região foi dividida em duas partes e, em virtude da projeção, remontada em uma área contígua: o limite comum é visível. É assim que está nos dados: é explicitamente representado como dois recursos separados. Para remover a divisão visível, você precisaria mesclar as duas partes (uma operação muito mais fácil com um GIS do que com R).
spmenciona que esta classe suporta "métodos incluindo impressão / exibição, plotagem, subconjunto, [, [[, \ $, nomes, dim, resumo, ..." A página de ajuda para subsetmostra sua sintaxe. O truque é saber qual campo usar para encontrar Moscou! Eu o encontrei olhando para valores únicos de todos os campos, finalmente chegando levels(gadm$NAME_1). Isso levou a spplot(subset(gadm.prj, subset=(NAME_1=="Moskva")),"regions"), o que faz o truque.
Vejo que whuber já me venceu, mas vou postar minha resposta de qualquer maneira.
Se você deseja manter os dados na projeção de longlat, pode adicionar o parâmetro " Longitudinal Wrapping " à string Proj4.
library(rgdal)
p4 <- projection(gadm)
p4_new <- paste(p4, "+lon_wrap=105")
gadm_new <- spTransform(gadm, CRS(p4_new))
spplot(gadm_new, "regions")Isso centralizará seu mapa na longitude 105 graus leste.
spTransform.