В настоящее время я пытаюсь получить список почтовых индексов, которые пересекаются с районами сообщества в Нью-Йорке. Я пытаюсь сделать это в R, используя объединение двух файлов GeoJSON или двух файлов shp, но в обоих случаях я сталкиваюсь с проблемами неидентичного CRS. Я не уверен, как обойти это.
Одной из основных проблем, с которыми я столкнулся, может быть то, что почтовые индексы и районы сообщества не сопоставляются 1-к-1, поэтому почтовые индексы будут появляться несколько раз по районам сообщества.
Я провел много поисков и попробовал следующее:
-функция 'rbind'
-использование spTransform с proj4string вторичного элемента GeoJSON
-функция «над»
-функция «укрепить»
- преобразование в фреймы данных
nycz <- readOGR("Zips.geojson")
nyccd <- readOGR("Community Districts.geojson")
nycz1 <- spTransform(nycz, CRS(proj4string(nyccd)))
p <- SpatialPointsDataFrame(nycz, data.frame(id=1:71))
q <- SpatialPointsDataFrame(nyccd, data.frame(id=1:71))
proj4string(q)
proj4string(p)
plot(new)
plot(p, col="red", add=TRUE)
res <- over(q,p)
table(res$boro_cd)
Мне нужна таблица почтовых индексов по районам сообщества.