Encuentre sublíneas, dadas las compensaciones inicial y final, a través de la línea de comando
Aug 19 2020
Tengo un archivo CSV que se ve esencialmente así:
id,id_wkt,start,end,wkt
11,11,0,43,"LINESTRING(3.3623676 51.3132714,3.3626663 51.3130068,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.3643705 51.3117777)"
12,11,43,217,"LINESTRING(3.3623676 51.3132714,3.3626663 51.3130068,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.3643705 51.3117777)"
Consiste en registros con una cadena lineal WKT. Algunos de los WKT aparecen varias veces y la geometría real debe obtenerse de la cadena lineal utilizando las compensaciones de inicio y fin (longitud en metros a lo largo de la línea) en sus respectivas columnas.
Me gustaría extraer cadenas de sublíneas en función de estas compensaciones de longitud. ¿Es esto posible a través de una herramienta de línea de comandos como OGR?
Respuestas
Socci Aug 18 2020 at 22:08
Después de algunas investigaciones, de hecho encontré una solución, usando ogr2ogr
:
ogr2ogr output.csv -lco GEOMETRY=AS_WKT input.csv -dialect SQLITE -sql "SELECT id, id_wkt, ST_Line_Substring(GEOMETRY, start / GeodesicLength(GEOMETRY), end / GeodesicLength(GEOMETRY)) FROM input"
Da
WKT,id,id_wkt
"LINESTRING (3.3623676 51.3132714,3.3626663 51.3130068,3.36274075608485 51.3129464307024)","11","11"
"LINESTRING (3.36274075608485 51.3129464307024,3.363068 51.3126811,3.3636323 51.3122504,3.3640754 51.3119614,3.36436609040842 51.3117804449745)","12","11"