jmt2080ad / polylinesplitter Goto Github PK

Really useful code! However, when I try splitting up a polyline into very small segments, I get an error. For example:

rivsplit <- splitLines(riv, dist = 1000) # this is OK
rivsplit <- splitLines(riv, dist = 10) # this results in an error

riv is a SpatialLine object with one feature. It is about 33 km long.
The error occurs at line 125 of polylineSplitter, at split(xydf,dist)

"Error in data.frame(x = x, y = y, end = end) : arguments imply differing number of rows: 7055, 7056"

Hi, nice set of functions. However, data on the sp object input is deleted. I'd like to keep the original data within the output, ie. duplicating rows.

I tried this, where Name is the original ID, however, it throws the folllowing Error object 'Name' not found:

splitLines2<- function(spobj, dist, start = T, sf = F){ xydf<-coordBuild(spobj) if (start == F){ xydf<-xydf[rev(rownames(xydf)),] } spoints <- split(xydf, dist) linelist <- list() lineslist <- list() id <- 1 Name <- Name if(!sf) { j <- 1 for(i in 1:(nrow(spoints)-1)){ linelist[j] <- Line(spoints[c(i, i + 1), c(1:2)]) j = j + 1 if(spoints[i+1,3] == 1){ lineslist[id]<-Lines(linelist, ID = id) id = id+1 Name = Name linelist<-list() j = 1 } } return(SpatialLinesDataFrame(SpatialLines(lineslist), data = data.frame(id = 0:(length(lineslist)-1), Name = Name))) } else { start <- 1 for(i in 1:(nrow(spoints)-1)){ if(spoints[i+1,3] == 1){ lineslist[[id]] <- sf::st_linestring(as.matrix(spoints[c(start:(i + 1)), c(1:2)], ncol = 2)) id <- id + 1 Name <- Name start <- i + 1 } } return(sf::st_sf(id = 1:length(lineslist), Name = Name, geom = sf::st_sfc(lineslist))) } }