在R中绘制人口普查Shapefile出现地理偏移问题求助

我已获取美国街区级人口普查数据，使用sf和rgdal包手动将其与街区Shapefile合并。尽管从读取Shapefile到合并的全过程中CRS保持不变，但绘制生成的sf对象时地图出现偏移——东北部部分街区明显超出城市范围，落入密歇根湖中。

数据来源说明 #

• Shapefile：美国人口普查局TIGER/Line数据库的2020年街区图层，使用的伊利诺伊州对应文件为tl_rd22_17_tabblock20.zip
• 人口普查数据：2020年人口普查P1表，针对伊利诺伊州库克县（FIPS代码17031）的街区级人口数据

最小可复现代码 #

library(ggplot2)
library(sf)
library(rgdal)
library(data.table)

dt_census <- fread(path_to_census_data)

blck_shp <- readOGR(dsn = block_shapefile_filepath, layer = "tl_rd22_17tabblock20")

blck_sf <- st_as_sf(blck_shp)

dt <- merge(blck_sf, dt_census, by.x = "GEOID20", by.y = "GEOID")

mini_dt <- dt[, c("POP20", "geometry")]
mini_dt$POP20 <- as.numeric(mini_dt$POP20)

plot(mini_dt)

问题现象 #

运行上述代码后生成的地图出现偏移，东北部部分街区明显超出城市范围，落入密歇根湖中：

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