Ich versuche eine Karte von ausgewählten kanadischen Provinzen / Territorien und ausgewählten US-Bundesstaaten zu erstellen. Bisher scheinen die schönsten Karten die zu sein, die mit GADM-Daten erstellt wurden:http://www.gadm.org/

Es war mir jedoch nicht möglich, die USA und Kanada auf derselben Karte zu zeichnen oder nur ausgewählte Provinzen / Territorien und Bundesstaaten. Ich interessiere mich beispielsweise für Alaska, Yukon, NWT, Britisch-Kolumbien, Alberta und Montana.

Außerdem scheint die US-Karte entlang der internationalen Datenlinie geteilt zu sein.

Kann mir bitte jemand helfen:

Zeichnen Sie die oben genannten Provinzen / Territorien und Bundesstaaten auf einer einzigen KarteVermeiden Sie eine Spaltung der USA entlang der internationalen DatenlinieÜberlagern eines Breiten- und LängengradgittersWählen Sie eine bestimmte Projektion, vielleicht die polykonische.

Möglicherweise erlaubt spplot Benutzern nicht, Projektionen anzugeben. Auf der Hilfeseite von spplot wurde keine Option zum Auswählen einer Projektion angezeigt. Ich weiß, wie man Projektionen mit der Kartenfunktion im Kartenpaket auswählt, aber diese Karten sahen nicht so gut aus und ich konnte auch nicht die gewünschte Teilmenge von Provinzen / Territorien und Staaten mit dieser Funktion zeichnen.

Ich kann nicht anfangen, ein Breitengrad-Längengrad-Gitter hinzuzufügen. Abschnitt 3.2 der Datei 'sp.pdf' scheint sich jedoch mit dem Thema zu befassen.

Unten ist der Code, den ich mir bisher ausgedacht habe. Ich habe jedes kartenbezogene Paket geladen, auf das ich gestoßen bin, und GADM-Daten auskommentiert, mit Ausnahme von Landes- / Gebiets- oder Staatsgrenzen.

Leider habe ich es bisher nur geschafft, Karten von Kanada oder den USA zu zeichnen.

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)

Danke für jede Hilfe. Bisher habe ich mit den obigen Schritten 2 - 4 keine Fortschritte gemacht. Vielleicht bitte ich um zu viel. Vielleicht sollte ich einfach zu ArcGIS wechseln und diese Software ausprobieren.

Ich habe diesen StackOverflow-Beitrag gelesen:

Kann R für GIS verwendet werden?

BEARBEITEN

Ich habe jetzt eine elektronische Kopie von "Applied Spatial Data Analysis with R" ausgeliehen. Bevand et al. (2008) und den zugehörigen R-Code und die zugehörigen Daten von der Website des Buches heruntergeladen (oder gefunden):

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

Ich habe hier auch ein paar gut aussehende GIS-bezogene R-Codes gefunden:

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

Wenn ich lerne, wie ich die gewünschten Ziele erreichen kann, werde ich hier Lösungen veröffentlichen. Obwohl ich möglicherweise irgendwann zu ArcGIS wechseln werde, wenn ich die Ziele in R nicht erreichen kann.