R sample to count locations in areas

countInAreas.R

##################################################################################
# Import school locations
# Import some areas
# Join counts
# write out

library(readxl)
library(writexl)
library(sf)

HOME <- Sys.getenv("HOME")
setwd(paste(HOME, "\\Documents\\LocalProjects\\SchoolsInAreas", sep=""))
getwd()

# first get the column names
cNames <- read_excel(path="Schools.xlsx", sheet="Analysis of directory", skip=1, n_max=1) %>% names()
# then the raw data
schs <- read_excel(path="Schools.xlsx", sheet="Analysis of directory",  skip=3, col_names=cNames[1:42])

schsSF <- st_as_sf(schs, coords = c("Longitude", "Latitude"), crs=4326)

pAreas = st_read("layers\\area.gdb")
pAreas_wgs84 <- st_transform(pAreas, 4326)

schoolsInAreas <- st_join(schsSF, pAreas_wgs84, join=st_intersects)

write_xlsx(schoolsInAreas, "schoolsInAreas.xlsx")