18.1. Análisis de direcciones. Diagrama de rosas

Análisis de direcciones (diagrama de rosas), a partir de lineamientos mapeados en un GIS.

Introducción: Se explica como hacer un diagrama de rosas para analizar datos direccionales a partir de los lineamientos mapeados en un GIS. Los datos pueden ser de dirección u orientación (datos axiales, ambivalencia de 180°). Además los datos pueden ser ponderados (o no) según su longitud para que los lineamientos de mayor magnitud tengan mayor peso. Se puede utilizar para el estudio de direcciones, estructuras, etc.

Script (link).


0. Datos: Archivo de texto plano (Datos.txt) con 2 columnas (Lon Lat). Cada lineamiento debe estar separado por el símbolo > (Ver como exportar los datos desde Global Mapper y QGis). No es necesario que cada linea este definida por 2 puntos.

1. Procesamiento de datos: Primero es necesario calcular el acimut (en sentido antihorario desde el Norte) y la distancia (en km) de los segmentos. Mapproject permite realizar esto. 

• –Af: permite incluir una nueva columna (3ª) con los datos de acimut. 
• –Gk: agrega una 4ª columna con la longitud del segmento medido en km (entre 2 puntos sucesivos). 
• -o3,2: selecciona las columnas con datos a extraer (4a y 3a, en ese orden espera los datos psrose).

gmt mapproject -fg "Datos.txt" -AF -Gk > "temp_rumbo" -o3,2

Archivo de salida SIN la opción -o3,2. Valores de acimut y distancia en 3a y 4a columna respectivamente.

Archivo de salida CON la opción -o3,2. Valores de distancia y acimut en 1a y 2a columna respectivamente.

2. Filas NaN: Es importante incluir previamente el siguiente comando para que se eviten analizar las filas sin datos (NaN). Caso contrario se distorsiona el diagrama de rosas.

gmtset IO_NAN_RECORDS skip

3. Analizar datos: psrose -I analiza los datos. Es importante que se agreguen los mismos argumentos que se utilizarán para hacer la figura (-D, -T) según el caso.
Debido a la presencia de NaN en las filas, la terminal muestra una advertencia antes de informar los datos de número de datos (n) y acimut promedio (mean az) entre otros parámetros.

gmt psrose "temp_rumbo" -I -D -T

4. Diagrama de Rosas: se dibuja con psrose a partir de los datos procesados (temp_rumbo). Es importante utilizar los argumentos de procesamiento según el tipo de dato.

Argumentos generales y estéticos:

• –R0/1/0/360: región en coordenadas polares (radio y dirección). Es decir primero los 2 valores del radio normalizado (0 a 1) y luego la región (de 0 a 360°).
• –Bx: divisiones del radio.
• –By90g30: divisiones de las direcciones. Anotaciones cada 90° y líneas cada 30°.
• –B+g: color del fondo del gráfico (lightblue).
• -LW,E,S,N: Letras de los puntos cardinales.
• -F: No dibuja la escala.
• –W: estilo de las líneas.
• –G: define el color naranja (orange) de los sectores.   
• –S5cn: define radio de la figura de 5 cm normalizada.

Argumentos de procesamiento:

• –A10: crea sectores cada 10°.
• –D: centra los datos en el sector.
• -T: Los datos son orientaciones en lugar de direcciones (ambivalencia de 180°).
• -Zu: Los datos NO son ponderados por su longitud (todos los segmentos tienen valor de 1). Por defecto SI son ponderados.

gmt psrose "temp_lineamiento.txt" >> %OUT% -R0/1/0/360 -Bx0.5g0.2 -By90g30 -B+glightblue -W0.5p -A10  -S5cn -Gorange -D

 

5. Agregar información: con pstext se agrega la información obtenida del análisis de los datos.

REM    Texto con info
    echo N = 206        | gmt pstext -R1/10/1/10 -JX10 -O -K >> %OUT% -F+cTL -Ya0.375c
    echo Mean Az = 285  | gmt pstext -R1/10/1/10 -JX10 -O -K >> %OUT% -F+cTL 

Diagrama con datos de direcciones (-T) y ponderado por la longitud (sin -Zu).

  

Diagrama con datos de orientación y valor (con -T y -Zu).

Diagrama con datos de direcciones y ponderado por la longitud (sin -T ni -Zu).

Diagrama con datos de direcciones (sin -T) y valor unificado (-Zu).

Más información:
Ver Ejemplo 6 de la galería de GMT para análisis a partir de direcciones.