Calcular orientación o dirección de líneas en QGIS

Al hacer adqusición de datos geofísicos en el mar, la embarcación debe realizar una serie de "líneas de adquisición" las cuáles por lo general son un conjunto de rutas paralelas entre si y otras perpendículares. En muchos casos se planifican para que sean norte-sur y este-oeste. Cada una de estas líneas es indentificada para facilitar el acceso a los datos asociados a ellas.

Recientemente recibí unas líneas de adqusición que formaban una cuadrícula, con variaciones debido a que las embarcaciones no pueden mantener un rumbo exacto durante su navegación. Estas líneas tenían un inconveniente, en realidad eran un conjunto muy grande de segmentos, cada segmento era un objeto espacial. La tarea que debía hacer era que cada línea se convirtiera en un objeto único en lugar de un conjunto de objetos. La primero que se me ocurrió hacer fue seleccionar las líneas de acuerdo a su orientación para luego unirlas. Sin embargo, primero debía saber la dirección de cada línea o segmenteo. Esto lo pude hacer incorporando un campo nuevo en el que se introdujo automáticamente el valor de la orientación o dirección de las líneas. La solución la obtuve aquí y a continuación les dejo la expresión que deben usar en el campo nuevo para optener la dirección de las líneas.

case
when yat(-1)-yat(0) < 0 or yat(-1)-yat(0) > 0 then 
(atan((xat(-1)-xat(0))/(yat(-1)-yat(0)))) * 180/3.14159 + 
(180 *
(((yat(-1)-yat(0)) < 0) + 
(((xat(-1)-xat(0)) < 0 AND (yat(-1) - yat(0)) >0)*2)
))
when ((yat(-1)-yat(0)) = 0 and (xat(-1) - xat(0)) >0) then 90
when ((yat(-1)-yat(0)) = 0 and (xat(-1) - xat(0)) <0) then 270
end

 Posteriormente publicaré el resto del procedimiento con el que solucioné mi problema. 

Crear base de datos espaciales local Spatialite con QGIS

Trabajar con bases de datos espaciales locales tiene grandes ventajas, una de las que más me gusta es la de poder contener en un solo archivo una gran cantidad de capas, a diferencia de los comúnmente usados shapefiles, que por cada capa existen diversos archivos.

A continuación explicaré lo sencillo que es crear una base de datos local Spatialite en QGIS.

Primero debes verificar que tienes el Panel del Explorador activo.

Si este panel no se encuentra activo, puedes habilitarlo haciendo clic derecho sobre la barra de herramientas y selecciona Panel del Explorador. Otra vía para hacerlo es ir a la barra de herramientas y seleccionar en el menú la opción “Ver” y luego “Páneles”.

Una vez activo el Panel de Explorador nos deplazamos hasta el final con el cursor hasta que podamos observar diferentes opciones de conexión.

Se nos mostrará la opción Spatialite. Al hacer clic derecho se tendremos la opción de conectar a una base de datos ya existente y la de crear una base de datos nueva. Bastará con hacer clic sobre la última opción, seleccionar la ubicación y el nombre de la base de datos nueva. De esta forma tendremos nuestra base de datos lista para importar nuestros archivos vectoriales.

Algoritmos de procesamiento vectoriales desaparecen después de actualizar a QGIS 2.16. Cómo solucionarlo.

Luego de actualizar de QGIS 2.14 a 2.16 desaparecieron los algoritmos de procesamietos vectorial del menú. Para solucionarlo en Linux desinstalé la aplicación y borré las carpetas /home/miusuario/.qgis2 y /usr/share/qgis/ . Luego de borrar, reinstalé la aplicación y pude tenera acceso a las herramientas del menú vectorial. Este error tiene que ver con que el plugin ftools ha sido añadido directamente como herramientas de procesamiento en la versión 2.16.

Saga GIS dejó de correr en Ubuntu

Siempre estoy probando software y haciendo experimentos con ellos. No soy computista, soy biólogo marino, pero me siento cómodo con los retos en informática y computación.

Recientemente tras algunas actualizaciones de algunos paquetes SAGA GIS dejó de funcionar, no iniciaba y obtenía un error "segmentation fault (core dumped)" en español Violación de segmento (core 'generado'). Esto me sucedía con Ubuntu 14.04 y Saga 2.2.4

El problema se presentaba porque entre tantas actualizaciones la biblioteca GDAL se actualizó a la versión 1.11 y SAGA GIS trabaja con la versión 1.10. Así que mediante Synaptic forcé la instalación de la la versión 1.10.1 y todo volvió a la normalidad. Creo que la actualización de GDAL se hizo con la instalación de GRASS 7,así que me quedé con la versión 6.4.3. Luego investigaré si puedo tener las dos versiones de GDAL.

Despliegue de datos en 3D usando QGIS

Para este ejercicio empleamos QGIS 2.6 y datos de la cuenca del río Canoabo, Venezuela. Los datos incluían un Modelo Digital de Elevación, una imagen Landsat, capa de vialidad.

Lo primero a realizar es cargar los datos. Luego se debe crear un nuevo Mapset en Grass. Prestar atención en la selección de la extensión del del Mapset, en este caso seleccione Venezuela de los predeterminados.

Haga clic sobre grass tolos y busque las herramientas de importar raster y vector desde la vista de QGIS.

Se ejecutaron los módulos designando el nombre de las capas que se crearán y luego se cargaron en el proyecto.

Busque la la herramienta nviz en las herramientas de GRASS. Seleccionamos el MDE como capa para determinar la elevacion y para este caso se seleccionó como capa de solape la de vialidad.

A partir de allí podremos visualizar en 3D los datos de la cuenca del Río Canoabo.