Перемещение Аляски и Гавайев на тематическую карту США с ggplot2

Ответ 1

Здесь, как это сделать, проецируя и преобразуя. Вам понадобится:

require(maptools)
require(rgdal)

fixup <- function(usa,alaskaFix,hawaiiFix){

  alaska=usa[usa$STATE_NAME=="Alaska",]
  alaska = fix1(alaska,alaskaFix)
  proj4string(alaska) <- proj4string(usa)

  hawaii = usa[usa$STATE_NAME=="Hawaii",]
  hawaii = fix1(hawaii,hawaiiFix)
  proj4string(hawaii) <- proj4string(usa)

  usa = usa[! usa$STATE_NAME %in% c("Alaska","Hawaii"),]
  usa = rbind(usa,alaska,hawaii)

  return(usa)

}

fix1 <- function(object,params){
  r=params[1];scale=params[2];shift=params[3:4]
  object = elide(object,rotate=r)
  size = max(apply(bbox(object),1,diff))/scale
  object = elide(object,scale=size)
  object = elide(object,shift=shift)
  object
}

Затем прочитайте в своем шейп файле. Используйте rgdal:

us = readOGR(dsn = "states_21basic",layer="states")

Теперь перейдите в равную область и запустите функцию fixup:

usAEA = spTransform(us,CRS("+init=epsg:2163"))
usfix = fixup(usAEA,c(-35,1.5,-2800000,-2600000),c(-35,1,6800000,-1600000))
plot(usfix)

Параметры - это повороты, масштабирование, сдвиг x и y для Аляски и Гавайев соответственно, и они были получены методом проб и ошибок. Тщательно их настройте. Даже изменение параметра шкалы Гавайи до 0.99999 отослало его с планеты из-за большого числа вовлеченных.

Если вы хотите вернуть это значение в lat-long:

usfixLL = spTransform(usfix,CRS("+init=epsg:4326"))
plot(usfixLL)

Но я не уверен, что вам нужно использовать преобразования в ggplot, так как мы сделали это с spTransform.

Теперь вы можете перейти через бизнес ggplot2 fortify. Я не уверен, имеет ли это значение для вас, но обратите внимание, что порядок состояний отличается от версии usfix. Аляска и Гавайи теперь являются последними двумя состояниями.

Ответ 2

Я опубликовал пакет fiftystater R на CRAN, чтобы предоставить простое решение. Он основан на шагах от ответа Spacedman (будет ссылка, но недостаточная репутация) и будет опубликован как кадр данных формы ggplot2::geom_map с именем fifty_states, чтобы удалить необходимость установки зависимостей, отследить файл исходной формы или настроить elide значения.

library(ggplot2)
library(mapproj)
library(fiftystater)

crimes <- data.frame(state = tolower(rownames(USArrests)), USArrests)

p <- ggplot(crimes, aes(map_id = state)) + 
  geom_map(aes(fill = Assault), map = fifty_states) + 
  expand_limits(x = fifty_states$long, y = fifty_states$lat) +
  coord_map()
p

^{карта пятидесяти состояний}

Плотность шума может быть очищена обычным способом, а в пакете также есть функция fifty_states_inset_boxes, чтобы добавить вставки:

p + scale_x_continuous(breaks = NULL) + 
  scale_y_continuous(breaks = NULL) +
  labs(x = "", y = "") +
  theme(panel.background = element_blank()) +
  fifty_states_inset_boxes()

^{пятьдесят состояний с вставками}

Ответ 3

Как только вы начнете перемещаться вокруг, вы можете просто представить Аляску и Гавайи как квадратные коробки где-нибудь в Мексиканском заливе. Это имело бы дополнительное преимущество, позволяя рассказать, какой цвет был на Гавайях.

Возможно, вы могли бы просто пройти весь путь и использовать искаженную систему, в которой каждое государство имеет равную площадь, тогда вы могли бы увидеть Род-Айленд.

Примеры картограмм США на изображениях Google показывают вид вещей. Не знаю, сколько из них имеют шейп файлы или данные с ними.

Вы действительно заинтересованы в относительных размерах состояний или хотите получить представление, которое позволяет людям узнать, какое значение для состояния?