arvi1000 · March 5, 2021 05:03 · arvi1000 · Mar 5, 2021
diff --git a/district_zip_overlap.R b/district_zip_overlap.R
 options(tigris_use_cache = TRUE)

 library(tigris)
 library(sf)

 # get shapefiles (coarse resolution)
 sd_sf <- school_districts(state = 'CA', cb = T)
 zip_sf <- zctas(cb=T, state='CA')

 # which zips overlap OUSD?
 ousd_sf <- sd_sf[grep('Oakland', sd_sf$NAME), ] 
 oakland_zip_indexes <- which(sapply(st_overlaps(zip_sf, ousd_sf), length) > 0)
 ousd_zips <- zip_sf[oakland_zip_indexes, ]

 # plot
 library(ggplot2)
 ggplot() +
  geom_sf(data=ousd_zips, color='grey', fill='white', alpha=0.5) +
  geom_sf(data=ousd_sf, color='red', fill=NA) +
  theme_void()

 # print the zips
 ousd_zips$ZCTA5CE10
	options(tigris_use_cache = TRUE)

	library(tigris)
	library(sf)

	# get shapefiles (coarse resolution)
	sd_sf <- school_districts(state = 'CA', cb = T)
	zip_sf <- zctas(cb=T, state='CA')

	# which zips overlap OUSD?
	ousd_sf <- sd_sf[grep('Oakland', sd_sf$NAME), ]
	oakland_zip_indexes <- which(sapply(st_overlaps(zip_sf, ousd_sf), length) > 0)
	ousd_zips <- zip_sf[oakland_zip_indexes, ]

	# plot
	library(ggplot2)
	ggplot() +
	geom_sf(data=ousd_zips, color='grey', fill='white', alpha=0.5) +
	geom_sf(data=ousd_sf, color='red', fill=NA) +
	theme_void()

	# print the zips
	ousd_zips$ZCTA5CE10