Promedio de polilíneas a partir de un número de líneas

Question

Tengo varias polilíneas con diferente longitud, dirección y curvatura. Quiero crear un promedio o conjunto de las polilíneas.

Lo que hice previamente en ArcGIS 10.1 es convertir las polilíneas en puntos en un tamaño de celda (25 km) que es bastante mayor que el área de estudio (10). De este modo, pude "promediar" los puntos. A continuación, volví a convertir las celdas ráster en líneas. En cierto modo se han "promediado", pero las líneas son muy rectas y forman pequeñas cajas en la intersección siguiendo el tamaño de la celda.

He leído algunas soluciones usando R. Me interesa que me den sugerencias para hacer esto en R.

Answer 1

Utiliza el paquete en mi repositorio de github: https://gitlab.com/b-rowlingson/meanline

Este es un ejemplo de uso. Primero haz algunos datos - es una lista de matrices de dos columnas de coordenadas x-y en las columnas:

 > xylist = list(cbind(1:10, rnorm(10)), cbind(1:10, rnorm(10)), cbind(1:10, rnorm(10)))

Entonces:

 > xyaverage = average_lines(xylist)

Aquí están las líneas y la línea media generada por lo anterior:

> plot(do.call(rbind,xylist), type="n")
> lapply(xylist,lines)
> lines(xyaverage, col="red")