Карты Тигра R: проблемы с картами Гуама и Марианских островов

У меня проблемы с файлами формы Тигра для Гуама и Марианских островов:

  1. Карта Тигра Гуам не содержит информации о муниципалитетах, однако многие наборы данных переписи включают муниципалитеты Гуама. Документация Тигра включает ссылку на муниципалитеты. Кто-нибудь может порекомендовать другой источник для полигональных / шейп-файлов муниципалитета Гуама?
  2. Марианские острова не загружаются правильно, используя функцию округов в Тигре. Я включил мой обходной путь, но я бы предпочел иметь доступ к этой информации напрямую.

,

library(tigris)
library(ggplot2)
library(ggthemes)
library(dplyr)

## Working Puerto Rico Map ----
# Downloads map
PR_map <- counties(72, cb = T) 
# Converts shapefile for ggplot2
PR_map <- ggplot2::fortify(PR_map, region = "COUNTYFP")
# Adds random data
PR_map <- merge(PR_map, data.frame(id = unique(PR_map$id), value = rep(1:4, (length(unique(PR_map$id)) %/% 4 + 4))[1:length(unique(PR_map$id))]), by = "id")
# Plots map
ggplot(PR_map)+geom_map(map = PR_map, aes(x=long, y=lat, fill = value, map_id=id), color="white", size=0.25)+ coord_map("mercator")+ theme_map()

## Broken Guam Map (No municipalities)----
# Downloads map
GU_map <- counties(66, cb = T) 
# Converts shapefile for ggplot2
GU_map <- ggplot2::fortify(GU_map, region = "COUNTYFP")
# Adds random data
GU_map <- merge(GU_map, data.frame(id = unique(GU_map$id), value = rep(1:4, (length(unique(GU_map$id)) %/% 4 + 4))[1:length(unique(GU_map$id))]), by = "id")
# Plots map
ggplot(GU_map)+geom_map(map = GU_map, aes(x=long, y=lat, fill = value, map_id=id), color="white", size=0.25)+ coord_map("mercator")+ theme_map()

## Broken Call for Mariana Islands 
# Does not download Mariana Islands, instead it downloads the entire United States and Territories 
USA_map <- counties(69, cb = T) 

# My work around is to subset to the appropriate region:
MI_map <- subset(USA_map, USA_map$STATEFP == "69")

# Converts shapefile for ggplot2
MI_map <- ggplot2::fortify(MI_map, region = "COUNTYFP")
# Adds random data
MI_map <- merge(MI_map, data.frame(id = unique(MI_map$id), value = rep(1:4, (length(unique(MI_map$id)) %/% 4 + 4))[1:length(unique(MI_map$id))]), by = "id")
# Plots map (The hex turns the map sideways for better viewing.)
ggplot(MI_map)+geom_map(map = MI_map, aes(x=long, y=lat, fill = value, map_id=id), color="white", size=0.25)+ coord_map("hex")+ theme_map()

1 ответ

  1. Решение вопроса Гуама заключается в использовании функции "места".
  2. Ошибка устранена. Кайл Уокер, сопровождающий Tigris, обновил приложение в ответ на этот пост и связанный твит. Смотрите твит.
Другие вопросы по тегам