多边形上的点的叠加在R中返回NA

多边形上的点的叠加在R中返回NA,r,polygon,spatial,R,Polygon,Spatial,编辑:根据以前的答案,用更具体的信息更新问题 很抱歉不够清晰,但这是我的第一篇文章。我有一组点,我想用over函数覆盖到加利福尼亚州海湾地区的县,这样我就可以得到一个包含点的数据框和一个新列,显示每个点所属的邻域。我看过其他几个相关的问题,但他们的答案似乎不适用于我使用的数据。我使用complete.cases删除带有NA的行。然后我将这些点转换为空间点数据帧。这些点的投影从NA开始,我将它们设置为“+init=epsg:3857”(谷歌地图/开放街道地图默认投影)。然后将投影点转换为多边形,即

编辑:根据以前的答案,用更具体的信息更新问题

很抱歉不够清晰,但这是我的第一篇文章。我有一组点,我想用
over
函数覆盖到加利福尼亚州海湾地区的县,这样我就可以得到一个包含点的数据框和一个新列,显示每个点所属的邻域。我看过其他几个相关的问题,但他们的答案似乎不适用于我使用的数据。我使用
complete.cases
删除带有
NA
的行。然后我将这些点转换为
空间点数据帧
。这些点的投影从
NA
开始,我将它们设置为
“+init=epsg:3857”
(谷歌地图/开放街道地图默认投影)。然后将投影点转换为多边形,即:
“+proj=lcc+lat_1=37.066667+lat_2=38.43333333+lat_0=36.5+lon_0=-120.5+x_0=2000000+y_0=500000.0000000001+datum=NAD83+units=us ft+no_defs+ellps=GRS80+towgs84=0,0,0”


我使用
result如果你在检查投影时得到
NA
作为你的分数,这意味着你没有正确设置投影。谢谢-我将此作为其他问题的解决方案,但他们的解决方案总是使用proj4string(分数)
id  lat         lon
1   37.97102    -122.300752
2   37.809869   -122.253598
3   37.782851   -122.45351
4   37.808497   -122.263405
5   37.767335   -122.425399
6   37.79215    -122.423548
7   39.78373    -100.445882
8   37.79215    -122.423548
9   39.78373    -100.445882
10  37.870594   -122.265938
11  39.78373    -100.445882
12  37.870594   -122.265938
13  39.78373    -100.445882
14  37.96979    -121.769799
15  NA          -100.445882