多边形上的点的叠加在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