R: создание карты выбранных канадских провинций и штатов США
Я пытаюсь создать карту выбранных канадских провинций / территорий и отдельных штатов США. Пока что самые хорошие карты, по-видимому, сгенерированы с помощью данных GADM:http://www.gadm.org/
Однако я не смог нанести на карту США и Канаду одну и ту же карту или нанести на карту только отдельные провинции / территории и штаты. Например, меня интересуют Аляска, Юкон, СЗТ, Британская Колумбия, Альберта, Монтана и другие.
Кроме того, карта США, кажется, разделена вдоль международной линии дат.
Может кто-нибудь, пожалуйста, помогите мне:
plot the aforementioned provinces/territories and states on a single map avoid having the U.S. split along the International dateline overlay a latitude-longitude grid select a specific projection, maybe the polyconic.Возможно, spplot не позволяет пользователям указывать прогнозы. Я не увидел опцию выбора проекции на странице справки spplot. Я знаю, как выбирать проекции с помощью функции карты в пакете карт, но эти карты выглядели не так хорошо, и я не мог построить желаемое подмножество провинций / территорий и штатов с этой функцией.
Я не знаю, как начать добавлять сетку широта-долгота. Тем не менее, раздел 3.2 файла "sp.pdf" кажется, чтобы обратиться к теме.
Ниже приведен код, который я придумал до сих пор. Я загрузил каждый связанный с картой пакет, на который я наткнулся, и закомментировал данные GADM, за исключением границ провинций / территорий или штатов.
К сожалению, пока мне удалось построить только карты Канады или США.
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)
Спасибо за любую помощь. До сих пор я не добился прогресса с шагами 2 - 4 выше. Возможно, я прошу слишком много. Возможно, я должен просто переключиться на ArcGIS и попробовать это программное обеспечение.
Я прочитал этот пост StackOverflow:
Можно ли использовать R для ГИС?
РЕДАКТИРОВАТЬ
Теперь я позаимствовал электронную копию «Прикладного анализа пространственных данных у R». Bevand et al. (2008) и загрузил (или обнаружил) связанный R-код и данные с веб-сайта книги:
Я также нашел здесь красивый R-код, связанный с ГИС:
https://sites.google.com/site/rodriguezsanchezf/news/usingrasagis
Если и когда я узнаю, как достичь желаемых целей, я опубликую решения здесь. Хотя я могу в конечном итоге перейти к ArcGIS, если не смогу достичь цели в R.