[英]over with points and polygons in R: get the name of the polygons
在我的data
我有一个带纬度/经度的信号列表。
我有一个形状文件,我用readOGR()
导入了它,我称之为polygons
。
随着代码
data$inside.polygons <- !is.na(over(data, as(polygons, "SpatialPolygons")))
我在一个新的变量data
称为inside.polygons
。 它是一个逻辑变量,描述信号是否在多边形内 (TRUE) 或不在 (FALSE) 内。
是否可以添加带有多边形名称名称的新列?
我创建了一个新表
polygons.table <- data.frame(polygons)
我在 $Polygon.name 中输入了每个多边形的名称
> str(polygons.table)
'data.frame': 233 obs. of 6 variables:
$ Country : Factor w/ 9 levels "Denmark","Estonia",..: 9 9 9 9 9 9 9 4 9 9 ...
$ Polygon.name: Factor w/ 237 levels "Aalborg","Aalborg Portland",..: 114 115 69 192 193 8 237 231 230 224 ...
$ Shape_Leng: num 0.0339 0.0209 0.0399 0.1628 0.1343 ...
$ Shape_Area: num 5.64e-05 2.26e-05 4.22e-05 5.25e-04 5.30e-04 ...
$ LOCodes : Factor w/ 193 levels "DEBOF","DEFLF",..: NA NA 155 184 184 137 193 15 191 192 ...
$ Links : Factor w/ 26 levels "http://eng.port-bronka.ru/",..: NA NA NA NA NA NA NA NA NA NA ...
如何将 Polygon.name 添加到多边形内的每个信号?
谢谢!
得到它:
inside.polygon2 <- over(data, polygons[,"Polygon.name"])
然后我将它添加到我的数据中
data$polygon.name <- inside.polygon2$Polygon.name
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.