QGIS: Complemento "Visibility analysis"

Quiero mostraros el complemento de QGIS para realizar estudios de observación a partir de puntos demominado VISIBILITY ANALYSIS, en concreto quería que conocieseis los tres tipos de análisis que puede realizar y sus funciones.

Para ello he realizado un estudio de observación de una capa de puntos sobre un MDE con un radio de visión de 5000 metros, probando los tres tipos de análisis que son:

Binary Viewshed: es el análisis más esencial, esto nos muestra las zonas visibles a partir de cada punto, en un radio de 5000 metros

También vemos que da otro valor a las celdas del raster de salida donde coinciden las zonas visibles de dos puntos

Depth bellow horizon: esto nos muestra un raster con las alturas que debe alcanzar cada punto de ese radio de 5000 metros para que se vea desde el punto de observación

Horizon: nos muestra los puntos visibles del borde del horizonte en el radio de 5000 metros, las podemos ver en color rojo

RELACIÓN ENTRE LAS MIGRACIONES DESDE LAS CCAA Y MADRID Y NÚMERO DE CONTAGIOS

En esta entrada quiero poner en relación los movimientos migratorios entre las CCAA y Madrid y su relación con el número de casos de contagio.

Creo que como geógrafo tenía que investigar un poquito si existe relación entre el número de casos con la movilidad de población entre las CCAA y Madrid, comunidad donde se han registrado el mayor número de contagios, ya que me resulta curioso que comunidades cuya población y actividad económica no es de las mayores de España tengan ese alto número de casos, y me resulta curioso la gran coincidencia que hay en los dos mapas que he elaborado entre la relación migratoria de las CCAA a Madrid y el número de contagios, pudiéndose ser esta una de las causas principales de la propagación en esas zonas.

Esto siempre es una opinión personal y puede ser errónea perfectamente, por ello pongo los mapas para que la gente saque las conclusiones que considere.

Ruta hacia el antiguo Puente de hierro del embalse de Ciudad Real

Esta es una bonita ruta turística hacia un emblema de Ciudad Real, el Puente de hierro, situado sobre el Embalse del Vicario, que unía la línea del ferrocarril Madrid - Badajoz, y que fue inaugurado en junio de 1927 el cual está declarado monumento histórico - artístico.

Más información aquí

Fuente: elsayon.blogspot.com

Haz clic en el mapa para acceder a la ruta

CREAR UN CAMPO AUTONUMÉRICO EN UNA TABLA CON UNA COLUMNA REORDENADA DISTINTA AL FID

Respondiendo a una pregunta de un lector del blog, vamos a ver como crear una columna autonmérica en una tabla donde hemos ordenado otra columna con un criterio distinto al FID, como vemos en este ejemplo:

Como vemos tenemos ordenada la columna área de menor a mayor, por lo que el campo FID se ha descolocado, entonces ahora queremos crear con este orden una columna nueva autonumérica correspondiente a este nuevo orden.

Después de comprobar que con Arcgis y Qgis no podemos realizar esta opción y la programación con python nos da siempre una autonumérica pero con el mismo orden que el FID, lo he realizado de una manera poco ortodoxa pero eficaz, siguiendo los siguientes pasos:

- Abrimos con Microsoft Excel el archivo .dbf del shapefile que contiene la tabla de atributos, y ordenamos la columna "área" de menor a mayor con la opción de Excel "ordenar y filtrar".

- Creamos una nueva columna que denominaremos "orden".

- En su primer campo ponemos un "1" y con el botón derecho del ratón pinchamos en la esquina inferior derecha y sin soltar arrastramos hasta el final de los campos de la tabla.

- Cuando soltamos el botón del ratón aparece un menú en Excel y elegimos en él "rellenar serie".

- Esto nos crea ya esta columna autonumérica que buscábamos con el orden de las áreas de menor a mayor, cuyo FID está desordenado.

- Guardamos la tabla en formato de Excel, quizás nos lo guarde con el nombre de "copia" y el nombre de la tabla.

- Abrimos Arcmap y abrimos el shapefile.

- Ahora abrimos en Arcmap la tabla de Excel que hemos creado, esto se nos abrirá solamente como tabla.

- Abrimos la tabla de atributos del shapefile y en la opciones de la tabla creamos un join con la tabla nueva cuyo campo común va a ser "área".

- Ya tenemos las tablas relacionadas y ya tenemos en la tabla del shapefile la nueva columna denominada "orden", con la nueva clasificación autonumérica basada en las áreas de menor a mayor.

- Ahora guardamos el shapefile con otro nombre mediante Data / export data.

- Abrimos el nuevo shapefile que hemos creado y abrimos la tabla de atributos, que debe tener los atributos del shapefile original más los del join, solamente nos queda eliminar los campos repetidos y dejar los campos originales más el campo "orden" que hemos creado.

Espero haber solucionado tu pregunta con este método, este proceso también se puede realizar con Qgis.

SERIE COMPLEMENTOS QGIS: STATION LINES

En esta ocasión vamos a hablar de un complemento curioso que puede ser útil para algunos trabajos que necesiten dividir algunas líneas de los terrenos en segmentos equidistantes, pudiendo elegir nosotros la distancia entre cada segmento, estamos hablando del complemento Station Lines de Qgis.

En este ejemplo vamos a ver las distintas opciones que tenemos con este complemento, el cual debemos descargarlo e instalarlo desde Administrar e instalar complementos de la pestaña complementos.

En un shapefile de líneas podemos crear segmentos solamente en los bordes de ellas, para ello abrimos el complemento Station Lines, en Layers elegimos el shapefile de lineas en el que queremos trabajar, también podemos seleccionar solamente las líneas a las que queremos aplicarles dicho complemento seleccionándolas en el shapefile y activando la casilla Use only selected features, en Station Choices elegimos Vértices, y damos a aceptar.

Si queremos dividir las líneas en segmentos equidistantes abrimos el complemento Station Lines y ahora en Station Choices elegimos Distance, en el apartado Distance elegimos la distancia entre segmentos que podemos escogerla mediante un campo de nuestra tabla de atributos eligiendo la opción Field o ponerla a mano con la opción Value, en el apartado Side podemos elegir a que lado de la línea se dibujarán los segmentos, derecha, izquierda o centro, en el apartado Length podemos elegir la longitud que tendrá cada segmento, y en el apartado Angle podemos inclinar los segmentos los grados que elijamos, vamos a ver algunos ejemplos.

SERIE COMPLEMENTOS QGIS: RASTER INTERPOLATION

Vamos a ver otro complemento interesante y muy útil de QGIS como es "Raster interpolation", mediante el cual podemos interpolar los valores de un raster a una capa de puntos, obteniendo en cada punto el valor de la celda que se corresponde con el dicho raster.
En este ejemplo voy a utilizar un raster de datos de elevación y un shapefile de puntos aleatorios.
Para ello debemos instalar el complemento "raster interpolation" desde la pestaña complementos / administrar e instalar complementos.
Abrimos el raster del que queremos obtener los datos, en mi caso un MDT, y abrimos el shapefile de puntos a los cuales queremos interpolarlos.

Para poder realizar la interpolación y aplicar el complemento debemos crear un campo en la tabla de atributos del shapefile de tipo numérico y ponerlo en modo edición.
Abrimos el complemento raster interpolation y donde dice Raster layer ponemos nuestro MDT, elegimos la banda, en mi caso es "0", en Destination layer elegimos el shapefile de puntos, en destination field elegimos el campo numérico que hemos creado para que aparezcan los datos, donde pone interpolation method elegimos el método de interpolación por el cual queremos que se realice, aunque en los tres métodos obtenemos unos datos prácticamente iguales, se puede dejar por defecto "bi-linear", también podemos activar la casilla process only selected features para que se realice la interpolación solamente en los puntos que hayamos seleccionado, y damos a "go".

Entonces nos aparecerán los datos de altitud de cada punto en el campo que hemos creado en nuestra tabla de atributos.

SERIE COMPLEMENTOS QGIS: PROFILE TOOL

Vamos a ver otro complemento interesante y útil de QGIS como es "profile tool", con el cual podemos crear un perfil de un terreno a partir de un MDE y una polilínea, ya sea creada directamente o a partir de un shapefile.

Descargamos el complemento profile tool de complementos / administrar e instalar complementos, una vez instalado accedemos a él en la pestaña complementos / profile tool / terrain profile.

Abrimos el MDE que nos interese para crear nuestro perfil topográfico y tenemos dos opciones, o crear una polilínea directamente sobre él seleccionando temporary polyline, o abriendo un shapefile de polilíneas y seleccionando la que nos interese con selected polyline.

Una vez creado el perfil topográfico podemos editarlo un poco en tamaño, posición, y guardarlo en varios formatos como pdf, png, svg.

En la pestaña table podemos copiar en un documento de texto o una hoja de cálculo todas las cotas del perfil topográfico con sus coordenadas.