我是GIS和R的新手，我正在做一个将犯罪数据映射到城市中特定社区的项目。数据文件为CSV格式，具有事件ID，经度和纬度。邻域多边形由shapefile定义。这是我在R中的代码：
require(sp)
require(rgdal)
require(maps)
crime<-read.csv("crime sample.csv")
crime<-subset(crime,!is.na(crime$Latitude))
coordinates(crime)<-c("Longitude","Latitude")
Neigh <- readOGR(".", "Neighborhoods_2012b")
proj4string(crime)<-proj4string(Neigh)
inside.Neigh <- !is.na(over(crime, as(Neigh, "SpatialPolygons")))
但是，结果表明没有点落在任何多边形内。
proj4string(Neigh)
给
[1] "+proj=tmerc +lat_0=36.66666666666666 +lon_0=-88.33333333333333 +k=0.999975
+x_0=300000 +y_0=0 +datum=NAD83 +units=us-ft +no_defs +ellps=GRS80 +towgs84=0,0,0
使用此链接的数据文件可用： file
谢谢你的帮助。