R 根据包含点的多边形（NYC civic地理空间数据）为点添加标签

我有纽约5449棵树的经纬度，以及55个不同邻里列表区域（NTA）的形状文件。每个NTA在shapefile中都有一个唯一的NTA代码，我需要在long/lat表中附加第三列，告诉我每个树属于哪个NTA（如果有的话）

我已经在stackoverflow上使用其他多边形中的点线程取得了一些进展，尤其是，但在尝试使用gContains时仍然会出错，并且不知道如何检查/标记不同多边形的每棵树（我猜是某种sapply还是for loop？）

下面是我的代码。可以在以下位置找到数据/形状文件：

库（rgdal）
图书馆（rgeos）
图书馆（GG2）
#导入数据
setwd（“”）
xy简单地尝试一下：

ShapeFile <- readShapeSpatial("Shapefile.shp")
points <- data.frame(long=xy$INTPTLON10, 
                         lat =xy$INTPTLAT10,
                         stringsAsFactors=F)
dimnames(points)[[1]] <- seq(1, length(xy$INTPTLON10), 1)
points <- SpatialPoints(points)

df <- over(points, ShapeFile)

ShapeFile我试过这么做，但它只给了我NAs？xy$nrow=超过（pts，空间多边形(nyc@polygons)，returnlist=TRUE）尝试1）不带任何参数-简单地说是空间点和空间多边形。2） 在我的应用程序中，我总是为坐标指定name-dimension，如
dimnames（coords）[[1]]我在这里尝试了代码，但我似乎得到了满表的NAs？原始xy表中有一些树缺少lon/lat，但我已将它们更改为0s…不确定这是否是导致NAs的原因？
ShapeFile <- readShapeSpatial("Shapefile.shp")
points <- data.frame(long=xy$INTPTLON10, 
                         lat =xy$INTPTLAT10,
                         stringsAsFactors=F)
dimnames(points)[[1]] <- seq(1, length(xy$INTPTLON10), 1)
points <- SpatialPoints(points)

df <- over(points, ShapeFile)