给定一组点位置(在NYC中为拾取位置)和限定感兴趣区域的多边形(在NYC中为邮政编码),我试图了解每个点是否在每个多边形内,如果是,则应确定哪个领取地点位于邮政编码所在的位置。
daf <- read.csv(file="Oct2308.csv", header=TRUE, sep=",")
sodo <- readOGR("ZIP_CODE_040114.shp", layer="ZIP_CODE_040114")
proj4string(sodo)
extent(sodo)
utm18nCRS <- crs(sodo)
utm18nCRS
plot.sodo<- SpatialPointsDataFrame(daf[,12:13],daf,proj4string = utm18nCRS)
crs(plot.sodo)
plot.sodo$PickUp_Zipcode <- over(plot.sodo,sodo)$ZIPCODE
write.csv(plot.sodo,file="~/Desktop/RFiles/ZipCodes/rfiles/unitedplotsodo.csv")
我的结果似乎是NA:
PickUp_Zipcode
NA
NA
NA
NA
NA