Análisis Espacio Temporal de Distribución de Terremotos con QGIS, PyQGIS y TimeManager

Los fenómenos de la naturaleza deben ser analizados en su magnitud, pero también en el tiempo y la ubicación, tanto para su comprensión como para su correlación con otros fenómenos. QGIS es un software de Sistemas de Información Geográfica de código libre que nos permite la representación de puntos, líneas y polígono en el espacio. Si bien QGIS esta enfocado en el análisis espacial, existe un complemento llamado TimeManager que complementa la representación temporal de los datos.

QGIS viene con una serie de funciones, herramientas y complementos. Para funcionalidades específicas o particulares se pueden utilizar códigos en Python (PyQGIS) que aceleran el procesamiento y análisis espacial. Este tutorial muestra el procedimiento completo de representación espacio temporal de datos mundiales de terremotos registrados en servidores de la USGS. El tutorial incluye scripts en Python, filtrados (querys), y asignación de estilos.

Link de base de datos de terremotos:

https://earthquake.usgs.gov/earthquakes/feed/v1.0/csv.php

 

Tutorial

 

Scripts

Aquí el script para cargar la tabla:

ruta="file:///C:/Users/Saul/Documents/Ih_AnalisisEspacioTemporaldeTerremotosconQGISyPyQGIS/all_month.csv"
campos="?type=csv&xField=longitude&yField=latitude"
rutatotal=ruta +campos
print(rutatotal)
vlayer = QgsVectorLayer(rutatotal, 'Terremotos', "delimitedtext")

Este el código para formatear la fecha:

from qgis.core import *
from qgis.gui import *

@qgsfunction(args='auto', group='Custom')
def mifuncion(dtg, feature, parent):
    dtg = str(dtg)
    #mi hora debe ser 2016-01-07
    fechahora = dtg[:4] + '-' + dtg[5:7] + '-' + dtg[8:10]
    #fechahora = dtg
    return fechahora

 

Datos de entrada

Descargue los datos para este tutorial en este enlace.