get_acs()设置geometry=TRUE时突发报错求助（tidycensus）

使用tidycensus包的get_acs()函数获取2017-2021年5年期ACS数据，设置geometry=TRUE时突然报错，此前代码运行无异常。

运行代码 #

my_vars <- c(
   total_pop = "B01003_001",
   pop_poverty = "B17001_001"
)

tracts_acs <- get_acs(
      geography = "tract",
      state = "CA",
      variables = my_vars,
      geometry = T) %>% 
   select(GEOID, variable, estimate) %>%
   st_transform(6440)

报错信息 #

获取2017-2021年5年期ACS数据
正在从Census网站下载要素几何数据。如需在后续会话中缓存shapefile，请设置options(tigris_use_cache = TRUE)。
使用州CA的FIPS代码'06'
错误：无法打开"/private/var/folders/zx/5k5b4tzx3gxfj_187kpqzq6w0000gn/T/Rtmpq3ED0D"；源文件可能已损坏或不支持。请查看st_drivers()获取支持的格式列表。
此外，警告信息：
In unzip(file_loc, exdir = tmp) : 从zip文件解压时出错1
错误：几何数据下载失败。请稍后重试或检查Census Bureau网站状态。

解决方法 #

• 启用本地缓存：在代码开头添加options(tigris_use_cache = TRUE)，让系统缓存已下载的shapefile，避免重复下载可能损坏的临时文件，后续会话可直接复用缓存数据。
• 清理临时文件：手动删除报错中提到的临时目录/private/var/folders/zx/5k5b4tzx3gxfj_187kpqzq6w0000gn/T/Rtmpq3ED0D下的所有文件，再重新运行代码。
• 更新依赖包：确保tidycensus、tigris和sf包处于最新版本，运行以下命令更新：

  update.packages(c("tidycensus", "tigris", "sf"))
  

• 分步获取并合并数据：拆分数据获取流程，先拿统计数据，再单独获取地理数据后合并：

  # 第一步：获取ACS统计数据（不加载几何信息）
  acs_stats <- get_acs(
    geography = "tract",
    state = "CA",
    variables = my_vars,
    geometry = FALSE
  ) %>% 
    select(GEOID, variable, estimate)
  
  # 第二步：获取加州tract地理数据
  ca_tract_geo <- tigris::tracts(state = "CA", year = 2021) %>% 
    st_transform(6440)
  
  # 第三步：合并统计数据与地理数据
  tracts_acs <- dplyr::left_join(ca_tract_geo, acs_stats, by = "GEOID")
  

• 检查数据源状态：确认Census Bureau的地理数据服务器运行正常，若服务器临时故障，可稍后重试。

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