Problema con el mapa trazado en R

Question

Mientras que el trazado de la gadm mapa de Rusia en R, tengo un problema en 180 grados de meridiano, donde parte de la región rusa se representa por separado. Se puede ver en el gráfico utilizando el siguiente código (es el borde izquierdo y la parte que está a la derecha). Puede alguien sugerir cómo evitar este problema? Gracias de antemano.

Aquí es una versión corta del código de mi pregunta:

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")

Answer 1

Veo que whuber ya me pegaba a él, pero voy a publicar mi respuesta de todos modos.

Si desea conservar los datos en la proyección longlat, puede Agregar el parámetro "Envoltura Longitudinal" a la cadena de Proj4.

library(rgdal)

p4 <- projection(gadm)
p4_new <- paste(p4, "+lon_wrap=105")

gadm_new <- spTransform(gadm, CRS(p4_new))
spplot(gadm_new, "regions")

Esto centrará el mapa en 105 grados este de longitud.

Answer 2

Proyecto los datos. En este ejemplo se utiliza una proyección estereográfica polar centrada cerca de la mitad de Rusia (longitud de 105 grados este):

require(rgdal)
proj4.str <- CRS("+init=epsg:3413 +lon_0=105")
gadm.prj <- spTransform(gadm, proj4.str)
spplot(gadm.prj, "regions")

El problema de la región es en la parte superior derecha. Usted puede ver que esta región ha sido dividida en dos partes y, en virtud de la proyección, a montar en un área contigua: el límite común es visible. Es de esa manera en los datos: es explícitamente representados como dos elementos independientes. Para quitar el visible split, usted tendría que combinar las dos piezas (una operación mucho más fácil de hacer con un SIG que con R).