On Fri, May 20, 2011 at 10:53 AM, Rija FANIRIANTSOA
<[email protected]> wrote:
>
> L<-SpatialLines(list(Lines(list(Line(cbind(c(7,6,5,4.2,3,2.8,1),
> c(4,3.5,3,2.8,2.5,1.8,1)))),"l1")))
> require(rgeos)
> gIntersects(L,polyg,byid=TRUE)
l1
p1 FALSE
p2 FALSE
p3 TRUE
p4 TRUE
- that gives you T/F for intersections of L with the polygons in
polyg. Job almost done!
Barry
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
