объект фрейма данных для использования в качестве объекта as.li nnet

Question

У меня есть объект data.frame, который я могу легко преобразовать в имя_пространства_d_frame , а затем преобразовать его в spatiallinesdataframe , но затем, когда я попытался покрыть как. li nnet он не читает отметки

         X       Y roadID
1 177321.3 3378163      1
2 177321.4 3378168      1
3 177321.4 3378168      1
4 177321.5 3378177      1
5 177321.5 3378186      1
6 177321.5 3378195      1

, затем я делаю этот data.frame в SpatialPointsDataFrame

coordinates(roaDF1) <- c("X","Y")
proj4string(roaDF1)=proj4string(trtrtt)

class       : SpatialPointsDataFrame 
features    : 100412 
extent      : 143516.4, 213981, 3353367, 3399153  (xmin, xmax, ymin, ymax)
crs         : +proj=utm +zone=17 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
variables   : 1
names       : roadID 
min values  :      1 
max values  :    347 

, затем используя SpatialPointsDataFrame для преобразовать SpatialLinesDataFrame (нашел этот код в сети)

LineXX <- lapply(split(roaDF1, roaDF1$roadID), function(x) Lines(list(Line(coordinates(x))), x$roadID[1L]))

linesXY <- SpatialLines(LineXX)
data <- data.frame(roadID = unique(roaDF1$roadID))
rownames(data) <- data$roadID
lxy <- SpatialLinesDataFrame(linesXY, data)
proj4string(lxy)=proj4string(trtrtt)

теперь lxy выглядит следующим образом

> lxy
class       : SpatialLinesDataFrame 
features    : 347 
extent      : 143516.4, 213981, 3353367, 3399153  (xmin, xmax, ymin, ymax)
crs         : +proj=utm +zone=17 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
variables   : 1
names       : roadID 
min values  :      1 
max values  :    347 

наконец-то я хочу согласовать этот lxy с объектом li nnet *

> W3.orig=as.linnet(lxy)
Warning messages:
1: In spatstat::linnet(vertices = V, edges = edges, sparse = TRUE) :
  edge list should not join a vertex to itself; ignored
2: In as.linnet.SpatialLines(lxy) :
  Internal error: could not map data frame to lines

это не потому, что roaDF1 имеет 100421 функций, когда я делаю его 20 КБ, все равно выдает ту же ошибку.

Любая помощь?

Answer 1

Это происходит потому, что данные координат линии включают в себя повторяющиеся вершины (т. Е. Две последовательные строки данных координат идентичны).

Существует предупреждение о «соединении вершины с самим собой», что является важной подсказкой.

Хотя в сообщении говорится, что это было «проигнорировано» (т. Е. Повторяющиеся вершины были удалены), код для обработки меток не смог обработать это, и выдает второе предупреждение и сдается.

Я напишу исправление для as.linnet.SpatialLines для обработки повторяющихся вершин. Вы все равно получите предупреждение!

PS: исправление для maptools::as.linnet.SpatialLines было отправлено сопровождающим maptools, и я ожидаю, что оно будет включено в следующий выпуск.