我正在尝试在点数据和面数据之间执行空间连接。
我的csv文件A中有指示事件空间坐标的数据,还有另一个文件shapefile B,该文件以多边形形式包含区域的边界。
head(A)
month longitude latitude lsoa_code crime_type
1 2014-09 -1.550626 53.59740 E01007359 Anti-social behaviour
2 2014-09 -1.550626 53.59740 E01007359 Public order
3 2014-09 -1.865236 53.93678 E01010646 Anti-social behaviour
head(B@data)
code name altname
0 E05004934 Longfield, New Barn and Southfleet <NA>
1 E05000448 Lewisham Central <NA>
2 E05003149 Hawcoat <NA>
我想将犯罪数据A添加到我的shapefile B中,以映射在我的区域A中发生的犯罪事件。不幸的是,由于code
A中的代码所引用的单位不同于B中的代码,因此我无法执行属性连接。
我已经阅读了许多教程和帖子,但找不到答案。我试过了:
joined = over(A, B)
和overlay
,但没有实现我想要的功能。
有没有办法直接进行这种连接,还是需要从A到另一种格式的中间转换?
从概念上讲,我想选择属于code
B区域的A点(类似于“基于ArcGIS中的空间位置进行连接”)。
有人遇到这个问题并解决了吗?
point.in.polygon
是,这是否将保留变量month
和crime_type
。你知道吗
point.in.poly
最终选择了属于相关多边形的那些点。谢谢。
point.in.polygon()
包裹了sp
吗?