使用R

Question

如何计算通过多边形的线数，给出以下内容。

1. 我从SpatialPolygonDataFrame访问多边形（例如SpatialPolygon @polygons [[1]]）
2. 这些行是SpatialLines

Here is a sample picture

Answer 1

假设spdf是SpatialPolygonDataFrame类的对象，sl是SpatialLines类的对象

尝试

libray(rgeos) overlap <- over(sl, spdf) sum(!is.na(overlap))

或直接

libray(rgeos) overlap <- length(sl[spdf,])