R地图 - 在地图上绘制点时过滤器坐标输出

我有一个事件位置（经度，纬度）的大型数据框。看到下面的示例：R地图 - 在地图上绘制点时过滤器坐标输出

events <- structure(list(lat = c(45.464944, 40.559207, 45.956775, 44.782831, 45.810287, 35.913357, 43.423855, 45.2359, 45.526025, 41.91371, 46.340833, 40.696482, 42.367164, 41.913701, 41.89167, 46.046206, 41.108316, 45.514132, 45.688118, 37.090387, 43.446555, 41.913712, 46.614833, 45.368825, 41.892168), lon = c(9.163453, 8.321587, 12.983347, 10.886471, 9.077844, 10.560439, 6.768272, 9.23176, 9.375761, 12.466994, 6.889444, 17.316925, 13.352248, 12.466992, 12.516713, 14.497758, 16.871019, 12.056176, 9.176543, 15.293506, 6.77119, 12.466993, 15.194904, 11.110711, 12.516085)), .Names = c("lat", "lon"), row.names = c(NA, 25L), class = "data.frame")

我想只绘制发生在特定国家（例如，意大利）的事件。

以下地块的所有事件：

library(rworldmap) 
library(maps) 
par(mar=c(0,0,0,0)) 
plot(getMap(resolution="low"), xlim = c(14,14), ylim = c(36.8,47)) 
points(x=events$lon, y=events$lat, pch = 19, col ='red' , cex = 0.6)

我怎样才能筛选出属于国家的边界以外的事件？预先感谢您的帮助和支持。

来源

2015-02-23 Sal

看看'SP :: over'。你可以用它来过滤只有特定多边形中的点。 – hrbrmstr 2015-02-23 11:47:49

分配多边形映射到对象：

m = getMap(resolution="low")

转换事件为sp对象 - 第一经度，纬度，然后：

pts = SpatialPoints(events[c(2,1)])

分配CRS：

proj4string(pts) = proj4string(m)

情节点在任何多边形内：

内意大利

points(pts[m,], pch = 3)

选择点：

it = m[m$ISO_A2 == "IT",] 
points(pts[it,], col = 'blue')

（这使用sp::over，但幕后）

来源

2015-02-23 14:28:14

谢谢。我花了整个下午尝试'sp :: over'，但总是错过了相关的参数，我最终无处可去。这一步一步的解决方案简单而高效。太好了！ – Sal 2015-02-23 16:08:53

'over'在'sp'的CRAN页面上有一个小插曲，以补充它的文档有点简洁。主要的用例是空间选择（参见上文）和聚合，两者都可以重复使用。 – 2015-02-23 16:54:40

只是想对@EdzerPebesma说+1 ......这个星期花了很多Google搜索终于偶然发现了这个简洁的答案，这对于像我这样的R新手来说也是有意义的。 – 2015-08-19 20:46:12

R地图 - 在地图上绘制点时过滤器坐标输出

回答

相关问题