如何解决将人口普查区分配给相交的学校边界
我正在尝试将人口普查区分配到学校出勤范围。
我有德克萨斯州学校的数据,我正在使用 NCES 学校出勤边界 (SAB) 数据集来查找位于这些 SAB 内的人口普查区域。
library(sf)
library(tigris)
library(dplyr) # for pipes
library(ggplot2)
# first: pull School Attendance Boundary Data from NCES: https://data-nces.opendata.arcgis.com/datasets/school-attendance-boundary-survey-2015-2016-1
# this query filters to TX
query <- "https://nces.ed.gov/opengis/rest/services/K12_School_Locations/SABS_1516/MapServer/0/query?where=stAbbrev%20%3D%20%27TX%27&outFields=*&outSR=4326&f=json"
tx_sabs <- st_read(query,stringsAsFactors = FALSE) # this takes a minute
# SABs are larger than census tracts
ggplot(tx_sabs) + geom_sf()
# grab all tx census tracts
tx_tracts <- tigris::tracts(state = "48",cb = T) %>%
sf::st_as_sf() %>%
select(STATEFP,COUNTYFP,tract = GEOID)
ggplot(tx_tracts) + geom_sf()
# set to same CRS
census_crs <- st_crs(tx_tracts)
# convert SAB crs to census CRS (4269)
tx_sabs <- tx_sabs %>%
st_transform(crs = census_crs)
# merge all census tracts within SABs into SAB
tx_sab_tracts <- sf::st_join(tx_sabs,tx_tracts,join = st_intersects,left = TRUE)
# tx_sab_tracts %>%
# ggplot() + geom_sf(aes(geometry = geometry)) +
# theme(legend.title = element_blank())
所以这里的 st_join
是将人口普查区编号分配给它与之相交的每个 SAB。由于左连接,我无法将人口普查区域形状与 SAB 形状进行比较,以查看重叠。所以我想我的问题是我怎么能确定它按照我想要的方式工作?我担心人口普查区被排除在 SAB 之外。
谢谢!
(编辑提及dplyr)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。