
Hello, I have been spending some time on Google looking for an
intersection of polygons and polylines. I have indeed a
SpatialPolygonsDataFrame (administrative boundaries) and thousands of
SpatialLines (outputs GCM's hurricanes tracks simulations),
the SpatialPolygonsDataFrame have an attribute (names of the polygon), I
want to find the intersection of the polygon and each polyline, and
create a new column for each polyline in the attribute of the polygon,
and fill 1 if the polyline hit the polygon, 0 ifelse.
This is a little example to express the case

library(sp)
P <- SpatialPolygons(list(Polygons(list(Polygon(cbind(c(1,1.5,2.8,2,1),
c(3,4,3.5,2.3,3)))), "p1"), 
Polygons(list(Polygon(cbind(c(1.5,3.5,4,3.5,2.8,1.5),
c(4,4.3,3.5,3.8,3.5,4)))), "p2"),
Polygons(list(Polygon(cbind(c(2,2.8,3.5,3,2.5,2),c(2.3,3.5,3.8,2,1,2.3))
)), "p3"),
Polygons(list(Polygon(cbind(c(3,3.5,4,3.8,3),c(2,3.8,3.5,1.5,2)))),
"p4")), 1:4)
nom<-LETTERS[1:4]
attr <- data.frame(nom, row.names = sapply(slot(P,
"polygons"),function(x) slot(x, "ID"))) 
polyg<- SpatialPolygonsDataFrame(P, attr)

dat<-as(polyg,"data.frame")
nom
p1 A
p2 B
p3 C
p4 D

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")))

#plot(polyg)
#plot(L,add=T)

The result that I want is like this

nom l1
p1 A 0 
p2 B 0
p3 C 1
p4 D 1

Thanks you!




<span id=m2wTl><font face="Arial, Helvetica, sans-serif" size="2" style="font-size:13.5px">_______________________________________________________________<BR>Obtenez votre e-mail GRATUIT avec 10 Mo sur <a href="http://www.mail2france.com" target="new">http://www.mail2france.com</a></font></span>
	[[alternative HTML version deleted]]

