En el artículo precedente se produjo en una hoja de cálculo una serie de coordenadas que correspondían a una rotación. Se cargaron en QGIS como capa de texto delimitado por espacios y se sugirió su conversión a vectorial, tipo polígono, con la extensión Points2One. Indagando sobre la posibilidad de hacer dicha operación directamente desde GRASS-QGIS, se descubrió que ello es posible con v.in.ascci. El módulo v.in.ascci admite dos formatos, tipo puntos y estándar. A pesar de que el archivo ascii puede cargarse para que se convierta directamente en vectorial tipo polígono, se considerarán los diferentes tipos de archivo ascii para convertir a vectorial punto, línea y polígono; respectivamente.
Para puntos, el archivo que inicialmente era éste (artículo precedente):
punto x y 1 639731.837577387 1005200.22069874 2 627197.968621817 1025465.79515429 3 622084.779997606 1046398.11172247 4 632550.434564209 1070843.95514829 5 645006.625153522 1089710.65383974 6 665151.262694107 1100229.68306209 7 679345.843517489 1100894.25445434 8 687711.304539127 1093719.20775055 9 697918.301157753 1085220.36713669 10 714927.261776934 1070449.16524505 11 712889.707408732 1046609.59547414 12 688962.613106175 1010131.42230662 13 667530.346401376 993320.069534544 14 652028.174142708 994452.767756846 15 639731.837577387 1005200.22069874
se convierte en:
639731.837577387|1005200.22069874 627197.968621817|1025465.79515429 622084.779997606|1046398.11172247 632550.434564209|1070843.95514829 645006.625153522|1089710.65383974 665151.262694107|1100229.68306209 679345.843517489|1100894.25445434 687711.304539127|1093719.20775055 697918.301157753|1085220.36713669 714927.261776934|1070449.16524505 712889.707408732|1046609.59547414 688962.613106175|1010131.42230662 667530.346401376|993320.069534544 652028.174142708|994452.767756846 639731.837577387|1005200.22069874
y se graba como rotado1.txt. Se ejecuta, en cónsola de GRASS-QGIS, el siguiente comando:
v.in.ascii input="/home/zeito/TEXTOS/rotado1.txt" out=rotado1
La visualización del vectorial se encuentra en la siguiente imagen:
Para incorporar el vectorial tipo línea el archivo siguiente, rotado2.txt (modo estándar), se encabeza con VERTI: para señalar el cuerpo de puntos y L 15 1 que indica que es un vectorial tipo línea con 15 puntos y una sola categoría. Termina en 1 1 (capa 1, categoría 1).
VERTI: L 15 1 639731.837577387 1005200.22069874 627197.968621817 1025465.79515429 622084.779997606 1046398.11172247 632550.434564209 1070843.95514829 645006.625153522 1089710.65383974 665151.262694107 1100229.68306209 679345.843517489 1100894.25445434 687711.304539127 1093719.20775055 697918.301157753 1085220.36713669 714927.261776934 1070449.16524505 712889.707408732 1046609.59547414 688962.613106175 1010131.42230662 667530.346401376 993320.069534544 652028.174142708 994452.767756846 639731.837577387 1005200.22069874 1 1
Se ejecuta, en cónsola de GRASS-QGIS, el siguiente comando:
v.in.ascii format=standard input="/home/zeito/TEXTOS/rotado2.txt" out=rotado2
cuya visualización es la siguiente:
Para importar un vectorial tipo polígono (modificador B en la cabecera del archivo) se requieren los puntos que constituyen el vectorial y su centroide (modificador C en la cabecera del archivo) dispuestos de esta forma:
VERTI: B 15 639731.837577387 1005200.22069874 627197.968621817 1025465.79515429 622084.779997606 1046398.11172247 632550.434564209 1070843.95514829 645006.625153522 1089710.65383974 665151.262694107 1100229.68306209 679345.843517489 1100894.25445434 687711.304539127 1093719.20775055 697918.301157753 1085220.36713669 714927.261776934 1070449.16524505 712889.707408732 1046609.59547414 688962.613106175 1010131.42230662 667530.346401376 993320.069534544 652028.174142708 994452.767756846 639731.837577387 1005200.22069874 C 1 1 667773.000000000 1049331.00000000 1 1
donde el 1 1 al final del archivo significa capa 1, categoría 1. La sintaxis del comando para importarlo en cónsola de GRASS-QGIS es:
v.in.ascii format=standard input="/home/zeito/TEXTOS/rotado3.txt" out=rotado3
y su visualización posterior:
Para finalizar, un ejemplo que incluye dos polilíneas con dos categorías:
VERTI: L 7 1 639731.837577387 1005200.22069874 627197.968621817 1025465.79515429 622084.779997606 1046398.11172247 632550.434564209 1070843.95514829 645006.625153522 1089710.65383974 665151.262694107 1100229.68306209 679345.843517489 1100894.25445434 1 1 L 8 1 687711.304539127 1093719.20775055 697918.301157753 1085220.36713669 714927.261776934 1070449.16524505 712889.707408732 1046609.59547414 688962.613106175 1010131.42230662 667530.346401376 993320.069534544 652028.174142708 994452.767756846 639731.837577387 1005200.22069874 1 2
Su visualización es similar al segundo caso con la salvedad de que su tabla atributiva contiene dos registros.
Pingback: Conversion de un archivo vectorial de puntos en un vectorial de línea usando un script de python (GRASS-QGIS) |