Estoy intentando crear un mapa de las provincias / territorios canadienses seleccionados y los estados seleccionados de los Estados Unidos. Hasta ahora, los mapas más bonitos parecen ser los generados con datos GADM:http://www.gadm.org/

Sin embargo, no he podido trazar los EE. UU. Y Canadá en el mismo mapa o parcela solo provincias / territorios y estados seleccionados. Por ejemplo, estoy interesado en Alaska, Yukon, NWT, British Columbia, Alberta y Montana, entre otros.

Además, el mapa de los Estados Unidos parece estar dividido a lo largo de la línea de datos internacional.

¿Puede alguien ayudarme por favor a:

trazar las provincias / territorios y estados mencionados en un solo mapaevitar que los Estados Unidos se dividan a lo largo de la línea de datos internacionalsuperponer una cuadrícula de latitud-longitudSeleccione una proyección específica, tal vez la policónica.

Tal vez spplot no permita que los usuarios especifiquen proyecciones. No vi una opción para seleccionar una proyección en la página de ayuda de spplot. Sé cómo seleccionar proyecciones con la función de mapa en el paquete de mapas, pero esos mapas no parecían ser tan agradables y tampoco pude trazar el subconjunto deseado de provincias / territorios y estados con esa función.

No sé cómo empezar a agregar una cuadrícula de latitud-longitud. Sin embargo, la Sección 3.2 del archivo 'sp.pdf' parece tratar el tema.

A continuación se muestra el código que he encontrado hasta ahora. He cargado todos los paquetes relacionados con el mapa que he encontrado y comentado los datos de GADM, excepto los límites provinciales / territoriales o estatales.

Desafortunadamente, hasta ahora solo he logrado trazar mapas de Canadá o los EE. UU.

library(maps)
library(mapproj)
library(mapdata)
library(rgeos)
library(maptools)
library(sp)
library(raster)
library(rgdal)

# can0<-getData('GADM', country="CAN", level=0) # Canada
  can1<-getData('GADM', country="CAN", level=1) # provinces
# can2<-getData('GADM', country="CAN", level=2) # counties

plot(can1)    
spplot(can1, "NAME_1") # colors the provinces and provides
                       # a color-coded legend for them
can1$NAME_1            # returns names of provinces/territories
# us0 <- getData('GADM', country="USA", level=0)
  us1 <- getData('GADM', country="USA", level=1)
# us2 <- getData('GADM', country="USA", level=2)
plot(us1)              # state boundaries split at 
                       # the dateline
us1$NAME_1             # returns names of the states + DC
spplot(us1, "ID_1")
spplot(us1, "NAME_1")  # color codes states and
                       # provides their names
#
# Here attempting unsuccessfully to combine U.S. and Canada on one map.
# Attempts at selecting given states or provinces have been unsuccessful.
#
plot(us1,can1)
us.can1 <- rbind(us1,can1)

Gracias por cualquier ayuda. Hasta ahora no he progresado con los pasos 2 a 4 anteriores. Tal vez estoy pidiendo demasiado. Quizás debería simplemente cambiar a ArcGIS y probar ese software.

He leído este post de StackOverflow:

¿Se puede usar R para GIS?

EDITAR

Ahora he tomado prestada una copia electrónica de 'Análisis de datos espaciales aplicados con R' Bevand et al. (2008) y descargado (o ubicado) el código R asociado y los datos del sitio web del libro:

http://www.asdar-book.org/

También encontré aquí un código R relacionado con GIS:

https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis

Si y cuando aprenda a cumplir los objetivos deseados, publicaré soluciones aquí. Aunque eventualmente me puedo mudar a ArcGIS si no puedo lograr los objetivos en R.