Re: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-26 Thread Michael Sumner
-Original Message- > From: Glenn Stauffer [mailto:gestauf...@gmail.com] > Sent: Friday, January 26, 2018 12:29 PM > To: 'Ege Rubak' <ru...@math.aau.dk> > Subject: RE: [R-sig-Geo] Splitting spatial lines at self intersections > (line > crossings) > > I think I am

Re: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-26 Thread Glenn Stauffer
Message- From: Glenn Stauffer [mailto:gestauf...@gmail.com] Sent: Friday, January 26, 2018 12:29 PM To: 'Ege Rubak' <ru...@math.aau.dk> Subject: RE: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings) I think I am really close to a solution using the RQGIS p

Re: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-26 Thread Marcelino de la Cruz Rot
You should load before library(maptools) and the coercion is done. Cheers Marcelino El 26/01/2018 a las 13:12, Michael Sumner escribió: Ege: I don't believe that this code is enough to get coercion between sp and spatstat: pts <- cbind(c(120:123,121,125),c(100,100,104,102,99,98)) pt2 <-

Re: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-26 Thread Michael Sumner
Ege: I don't believe that this code is enough to get coercion between sp and spatstat: pts <- cbind(c(120:123,121,125),c(100,100,104,102,99,98)) pt2 <- cbind(c(124,124,123,118,124,125),c(100,97,100,104,106,110)) library(sp) library(spatstat) L1 <- SpatialLines(list(Lines(list(Line(pts)), "X")))

Re: [R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-25 Thread Ege Rubak
The `psp` class in spatstat consists of individual line segments and `selfcrossing.psp` checks whether each individual line intersects one of the other lines, which happens at all the points in your plot. If instead you treat each of the two line sequences as a `psp` you can check find the

[R-sig-Geo] Splitting spatial lines at self intersections (line crossings)

2018-01-25 Thread Glenn Stauffer
I have a Spatial Lines object I would like to split at every point where the line self-intersects (crosses or touches itself), or anywhere lines touch each other, if there are multiple lines. I've tried using spatstat to convert the SpatialLines to a psp object and then use the selfcrossing.psp