【发布时间】：2021-06-17 08:11:27
【问题描述】：
我想计算我所在州（北卡罗来纳州）邮政编码的人口密度。我能够从美国人口普查中提取邮政编码人口和多边形，并使用以下代码绘制北卡罗来纳州地图：
library(tidycensus)
geo <- get_acs(geography = 'zcta', # Get zip code-level data
state = 'NC', # for NC
year = 2019, # for 2019
table = 'B01003', # from this US Census table
geometry = TRUE) %>% # and return the geospatial polygons for mapping
dplyr::rename('population' = estimate, 'zipcode' = NAME) %>%
select(-moe) %>%
arrange(zipcode)
p <- tm_shape(geo) +
tm_polygons('population')
p
这会按邮政编码映射人口。为了按邮政编码计算和映射人口密度，我需要每个邮政编码多边形的面积（以英里或平方公里为单位）。我正在努力寻找一种方法来 (a) 从美国人口普查网站获取这些数据，(b) 在其他地方找到它，或者 (c) 使用多边形几何来计算它。
任何建议将不胜感激。
【问题讨论】：
标签： geospatial tmap tidycensus