Importación y conversión directa de archivos de texto a vectoriales tipo punto, línea o polígono con v.in.ascii en GRASS-QGIS

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.

About these ads
Esta entrada fue publicada en GRASS, QGIS, SIG, Software Libre. Guarda el enlace permanente.

Una respuesta a Importación y conversión directa de archivos de texto a vectoriales tipo punto, línea o polígono con v.in.ascii en GRASS-QGIS

  1. Pingback: Conversion de un archivo vectorial de puntos en un vectorial de línea usando un script de python (GRASS-QGIS) |

Deja un comentario

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s