使用zctas()获取北卡2020年ZCTA代码报错，如何解决及获取普查区块边界？

你遇到的报错是因为tigris包的zctas()函数仅在2000和2010年支持按州筛选ZCTA数据，2020年的ZCTA无法通过该方式按州获取。但要获取2020年北卡罗来纳州的人口普查区块边界，直接用tigris包的blocks()函数就能实现，全程无需手动下载或处理Shapefile。

实现代码 #

首先加载依赖包：

library(tidyverse)
library(tigris)
library(sf)

调用blocks()函数获取全州区块数据：

# 获取北卡罗来纳州2020年人口普查区块边界
nc_blocks_2020 <- blocks(state = "NC", year = 2020)

• 参数说明：state = "NC"指定北卡罗来纳州（也可以用州FIPS代码"37"），year = 2020指定2020年普查数据，函数会自动完成数据下载并导入为sf格式的空间对象。

如果需要筛选特定县的区块，比如达勒姆县（FIPS代码063），可以用dplyr过滤：

# 筛选达勒姆县的2020年普查区块
nc_durham_blocks_2020 <- nc_blocks_2020 %>%
  filter(COUNTYFP == "063")

内容的提问来源于stack exchange，提问作者alocklear