Apuntes de Qgis 3.0 PDF
Apuntes de Qgis 3.0 PDF
Apuntes de Qgis 3.0 PDF
0 GIRONA: NOVEDADES Y
DESCARGA. LA REVOLUCIÓN DEL
SOFTWARE LIBRE ESTÁ AQUÍ
Para todos aquellos que aún no conozcan QGIS y que comienzan a interesarse en este potente
software comentar que se trata de un Sistema de Información Geográfica de Software Libre, siendo
un proyecto de la Open Source Geospatial Foundation (OSGeo).
Desde el inicio del proyecto en el año 2002 hasta la actualidad, QGISha experimentado una
evolución que le ha llevado a ser un programa GIS que cubre la mayor parte de las necesidades en
visualización de datos.
QGIS se encuentra bajo la General Public License (GPL), lo cual permite al usuario modificar su
código fuente y garantizar la existencia y acceso a un programa GIS gratuito.
Descubre por qué QGIS es el SIG líder de Código Abierto para escritorio.
Desde 2015, el conjunto de desarrolladores del proyecto, vienen trabajando en la siguiente versión:
la 3.0, cuyo lanzamiento se programó para el pasado mes de diciembre de 2017. Sin embargo,
alcanzada dicha fecha, éste no se produjo y se retrasó hasta el 23 de febrero de 2018. Y ahora sí,
por fin, ya tenemos disponible la nueva versión QGIS 3.0
El motivo de que esta versión haya sufrido esta serie de retrasos en su lanzamiento, se debe a la
complejidad de llevar a cabo la actualización a la librería Qt5 y Python 3. Un gran impacto necesario,
no sólo para la interfaz del software, sino también para que el conjunto de plugins sea compatible
con la nueva API.
1
En este sentido, no todos los plugins disponibles en la versión anterior van a estarlo en la versión 3.0,
se irán incorporando paulatinamente (aunque algunos no se actualizarán, esto depende del
desarrollador de cada plugin).
Detrás del proyecto hay una comunidad, formada por un comité directivo, un equipo de
desarrolladores y los propios usuarios del software. Las contribuciones de esta comunidad pueden
ser en forma de código, informes o correcciones de errores, documentación, promoción y
apoyo. MappingGIS apoya el proyecto como patrocinadores desde noviembre de 2014.
Índice [cerrar]
• QGIS 3.0 Girona
• Descarga de QGIS 3.0 Girona
• ¿Por qué en esta ocasión la versión 3.0 se llama Girona?
• Novedades de QGIS 3.0
o Interfaz de usuario
o Simbología
o Etiquetado
o Diagramas
o Representación
o Digitalización
o Gestión de datos
o Formularios y widgets
o Leyenda de capa
o Compositor de mapas
o Herramientas de análisis
o Procesamiento
o Opciones de aplicación y proyecto
o General
o Proveedores de datos
o Plugins
o QGIS Server
o Programabilidad
Según la hoja de ruta, hasta finales de año y más concretamente, hasta el 26 de octubre de 2018, la
versión LTR seguirá siendo QGIS 2.18 Las Palmas de Gran Canaria.
QGIS 3.0 es una LR (versión de lanzamiento regular) que inicia una nueva fase de desarrollo.
Después del lanzamiento de la versión 2.16 comenzó el trabajo en el nuevo lanzamiento previsto
para la versión 3.0, lo que implicó un cambio en el calendario de lanzamientos.
La versión 3.0 implica muchos cambios de gran alcance, por lo que no está claro por adelantado
cuánto tiempo llevará “acomodar” el conjunto de plugins que se venían disfrutando en la versión 2.x
2
DESCARGA DE QGIS 3.0 GIRONA
Puedes descargar QGIS 3.0 Girona desde la página oficial de descargas, seleccionando el
instalador acorde a la versión de tu sistema operativo.
3
¿POR QUÉ EN ESTA OCASIÓN LA VERSIÓN
3.0 SE LLAMA GIRONA?
QGIS 30. se denomina Girona debido a que esta fue la ciudad anfitriona de la decimoquinta reunión
de desarrolladores en mayo de 2016, justo después de la II Conferencia Internacional de Usuarios y
Desarrolladores de QGIS y las Jornadas SIG Libre
Esta es la decimotercera versión que lleva un nombre de una localidad terrestre, que desde la versión
1.6 adoptaron nombres terrenales: 1.6 Copiapó, 1.7 Wrocław, 1.8 Lisboa, 2.0 Dufour,
2.2 Valmiera, 2.4 Chugiak, 2.6 Brighton, 2.8 Wien, 2.10 Pisa, 2.12 Lyon, 2.14 Essen, 2.16
Nødebo y 2.18 Las Palmas.
El motivo de este modo de designar las versiones de QGIS tiene como objetivo reducir los riesgos
legales que tienen las marcas.
Interfaz de usuario:
2. Mejora en el zoom con la ruleta del ratón en equipos Mac para un uso más fino.
4
4. Se incluye la posibilidad de alternar la visibilidad de los paneles abiertos en la interfaz mediante
el menú Ver/Alternar visibilidad del panel, lo que permite una mejor visualización de los datos.
5
5. Nueva barra de localizadores.
6. Capacidad de guardar el lienzo como imagen o, incluso, como pdf.
6
7. Nuevo cuadro de diálogo para seleccionar color de fondo en la flecha de norte.
Simbología
7
13. Nuevo panel de estilos para administrar la simbología.
14. Mejoras en el Administrador de estilos. Por ejemplo, como añadir a favoritos y crear grupos
de símbolos.
8
15. Soporte para modificar el color y la transparencia de un conjunto de valores al mismo
tiempo.
Etiquetado
Diagramas
20. Propiedades definibles para los datos. Se ha agregado soporte para más propiedades de
datos en diagramas, tales como:
• Distancia
• Prioridad
• Índice Z
• Mostrar
• Ángulo de inicio de diagrama
9
Representación
21. Soporte para “capas en vivo”. Es posible generar capas de forma automática en un intervalo
de tiempo que determinemos. Estas capas son excelentes para actualizar datos que cambian
continuamente.
22. Soporte para agrupar los puntos cercanos de una capa en un solo símbolo.
10
23. Representación en cuadrícula para desplazamiento de puntos.
Digitalización
Gestión de datos
11
28. Revisión de metadatos.
29. Función Pan a entidad desde la tabla de atributos.
12
Formularios y widgets
13
34. Agregar acciones de ámbito de capa.
35. Soporte para mostrar los valores de campo en modo autocompletar dentro del formulario.
36. Control para zoom y flash a entidades en la selección dentro del cuadro de diálogo.
37. Función de selección entre grupo de valores en campos numéricos dentro del formulario.
Leyenda de capa
Compositor de mapas
40. Control de la representación de las líneas del grid. Tanto verticales como horizontales.
41. Creación de una nueva composición a partir de una plantilla qpt, que podemos arrastrar y
soltar a la interfaz de QGIS.
42. Control del espaciado entre las etiquetas de los elementos de la leyenda de mapa.
43. Posibilidad de establecer diferentes SRC para los elementos del mapa.
14
44. Nueva función para limitar los ángulos de una línea. Tanto al dibujar una polilínea como un
polígono.
Herramientas de análisis
45. Posibilidad de ejecutar herramientas de procesamiento con los valores únicos de un
ráster.
Procesamiento
• Aspect (Orientación)
• Slope (Pendiente)
• Ruggedness (Rugosidad)
• Hillshade (Modelo de sombras)
• Relief (Relieve)
15
53. Nuevo algoritmo para extracción de nodos.
54. Integración del complemento de Estadísticas zonales dentro de la Caja de herramientas
de procesado.
55. Nuevo algoritmo para consultas SQL SpatiaLIte.
56. Nuevo algoritmo para extender líneas.
57. Nuevo algoritmo para extraer geometrías basado en una expresión.
58. Nuevo algoritmo para importar dentro de bases de datos SpatiaLite.
59. Nuevo algoritmo de interpolación (IDW / TIN).
60. Nuevo algoritmo de cálculo de geometría por expresión.
61.Nuevo algoritmo para el cálculo del punto más lejano al límite de un polígono.
62. Soporte para extraer atributos con valores nulos o no nulos.
63. Nuevo algoritmo para eliminar geometrías.
64. El plugin heatmap ha sido añadido a la caja de herramientas de procesado.
65. Sustitución del plugin RoadGraph por un conjunto de algoritmos para análisis de redes
en la caja de herramientas de procesado.
66. Posibilidad de exportar un modelo de procesado a pdf o svg.
16
71. Mejoras en la ventana del modelador de procesado. Por ejemplo, el control de zoom en la
barra de herramientas de la ventana, personalización de los datos de entrada y la ubicación de los
paneles de la ficha algoritmos.
17
73. Se ha eliminado el proveedor TauDEM de la Caja de herramientas de procesado y se
implementa como plugin.
74. Nuevo algoritmo para descarga de archivos desde URL.
General
77. Se han eliminado plugins para integrarlos en el bloque de procesado. Por ejemplo, el
conversor de dxf a shp o estadísticas zonales.
78. Posibilidad de configurar la ubicación de los archivos de ayuda de QGIS.
79. Se han eliminado los proveedores OTB y LiDAR de la Caja de herramientas de procesado.
Proveedores de datos
80. Soporte para GeoPackage y para todos los formatos ráster editables admitidos por GDAL.
81. Nuevo cuadro de diálogo para añadir capas unificado.
18
87. Integración de GeoNode. GeoNode es un sistema de gestión de contenidos geoespacialesde
código abierto que facilita su publicación en web.
Plugins
19
QGIS Server
91. Se ha revisado completamente QGIS Server. Todos los servicios (WMS GetMap, WFS
GetFeature, GetLegendGraphics, WCS, GetPrint, etc.) han sido reescritos parcialmente. Algunos,
como WMS, han sido completamente reescritos.
Se han incluido a demás las siguientes características:
20
Programabilidad
97. Desaparece la aplicación QGIS Browser para gestión de archivos, ahora toda la gestión se
realiza desde el panel Navegador.
21
Novedades en la tabla de atributos de QGIS 3
La tabla de atributos permite administrar los datos SIG, además permiten el análisis, acceso y
actualización de la información almacenada. En general, una tabla consta de filas y columnas. Cada
fila representa un objeto espacial, y cada columna almacena la información del objeto espacial.
Para acceder a la Tabla de Atributos de una capa activa, seleccione una capa y haga clic con el
botón derecho, luego en el menú desplegable seleccione Abrir tabla de atributos. También se puede
dirigir al menú Capa > Abrir tabla de atributos o pulsar la tecla F6.
En el tope de la Tabla de Atributos, se encuentra una barra descriptiva que muestra el nombre de la
Capa seguido del total de registros, cuantos están seleccionados y filtrados.
22
¿Cuál es la diferencia entre seleccionar y filtrar?
La selección resalta con un color los registros en el mapa y la tabla, el filtrado muestra los registros
que cumplen un criterio y oculta el resto.
Barra de Herramientas
Seguidamente de la barra descriptiva esta la barra de herramientas cuyos botones realizan las
siguientes funcionalidades:
1. Activa el modo de edición, permite crear nuevos campos, introducir, modificar y borrar información.
2. Activa el modo multiedición. Nuevo en QGIS 3.
3. Guardar los cambios.
4. Refresca la tabla para visualizar cambios realizados.
5. Adicionar registros.
6. Borrar el registro seleccionado.
7. Cortar los registros seleccionados.
8. Copia los registros seleccionados.
9. Pega los datos a un nuevo registro.
10. Selecciona registros utilizando el Dialogo de Expresiones.
11. Seleccionar todos los registros.
12. Invertir la selección.
13. Remover toda la selección. No afecta los registros solo los deja sin seleccionar.
14. Selecciona o filtra los registros utilizando formularios. Nuevo en QGIS 3.
15. Colocar los registros seleccionados en la parte superior.
16. Desplazar el mapa a las entidades seleccionadas.
17. Acercar el mapa a las filas seleccionadas. Si tenemos un registro seleccionado acerca el mapa a dichas
entidades espaciales.
18. Borra un campo y todos los datos almacenados en la columna.
19. Crea un nuevo campo.
20. Abre la calculadora de campos.
21. Formato condicional. Esta herramienta ya existía en versiones anteriores, en QGIS 3 se encuentra alineada
con el resto de botones.
Modo multiedición: permite editar los datos de múltiples registros sin necesidad de utilizar la
calculadora de campos, funciona en edición y modo formulario, es útil cuando deseas colocar el
mismo dato en varios registros, simplemente se seleccionan los registros haciendo clic en la caja de
chequeo ubicada al lado izquierdo. Al escribir en uno de los campos se reescriben los datos en todos
los seleccionados
23
Selecciona o filtra registros utilizando formularios: al igual que la herramienta anterior funciona
en modo formulario, facilita la búsqueda, selección y filtrado de registros por atributos, especialmente
la ejecución de consultas condicionales complejas. En el ejemplo de la figura se seleccionan los
estados de México con una superficie mayor de 10.000 Km2 y una población menor o igual a
1.000.000 de habitantes. En el caso de los campos de texto cuenta con autocompletación para
facilitar la búsqueda y selección.
24
Modos de trabajo:
En QGIS 3 podemos gestionar las tablas en dos modalidades, Tabla y Formulario, para alternar entre
ellas, hacer clic en los botones ubicados en la esquina inferior derecha.
Modo Tabla
25
Modo formulario
Toma más relevancia en QGIS 3, nuevas funcionalidades trabajan solo en esta modalidad
Existen otras operaciones que podemos realizar con las tablas, visualizar solo los registros de las
entidades visibles en el mapa, recién creadas o editadas, vincular tablas de atributos con otros datos,
calcular propiedades geométricas de las entidades, calcular estadísticas e incluso ejecutar scripts
Python. Más adelante iremos cubriendo estos tópicos.
26
OPERACIONES BÁSICAS CON TABLAS EN
QGIS 3
Es posible realizar operaciones en las tablas de atributos de QGIS 3, por ejemplo consultas, cálculo de
propiedades geométricas, concatenación y extracción de textos, operaciones matemáticas y con
fechas.
Las operaciones en las tablas de atributos se ejecutan en: en la Barra de Cálculo de Campos,
ubicada inmediatamente debajo de la barra de herramientas de la tabla de atributos o con
la Calculadora de campos.
Calculadora de Campos:
La calculadora de campos está compuesta por tres paneles, al medio muestra las expresiones que
son las funciones que se puede utilizar (en este artículo se introdujo el uso de expresiones para crear
simbología), en el panel derecho se escribe la operación y el panel izquierdo muestra una
breve explicación con ejemplos de lo que hace la expresión seleccionada. En la parte superior se
configuran las opciones de salida, ya sea en un nuevo campo, campo virtual o actualizar un campo
existente.
27
Operaciones sobre campos de texto
Los campos de texto son denominados Cadenas, haciendo referencia a cadenas de caracteres,
la Calculadora de Campos cuenta con expresiones para hacer operaciones con texto, concatenar,
extraer caracteres o frases, convertir a mayúsculas o minúsculas, convertir números a textos o
viceversa, entre otros.
1.- Rescribir atributos: si deseamos asignar un texto o valor numérico a varios registros, primero
seleccionamos los registros y luego se escribe el valor entre apostrofes.
2.- Concatenar campos: podemos atributos de texto, si uno de los campos es numérico
automáticamente es convertido en texto. En el siguiente ejemplo se tienen los datos de una red de
abastecimiento, un campo con el diámetro en pulgadas (decimal) y otro con el material (texto), se
creara un nuevo campo que los una denominada “Tipo”, para ello se utiliza la función “concat”:
concat(“Material”, ‘ de ‘ , “Diametro”, ‘ ” ‘)
28
1.- Densidad de población: Para calcular la densidad de población para 1990 se utiliza la siguiente
expresión “Pop1990” / “Area”, se configura un nuevo campo denominado “Densidad90”, tipo decimal
con 12 caracteres y 2 decimales:
2.- Incremento de la población: Para conocer el incremento de población entre los años 1990 y 1997,
sencillamente realizamos la resta de los campos “Pop1997” – “Pop1990”.
1.- Calcular las coordenadas de puntos: Abra la calculadora de campos, en la caja de dialogo que se
despliega deja la opción por defecto, “Crear un Nuevo Campo”, en la pestaña “Nombre” introduzca
“Este”. En la parte central de la caja de dialogo busque las funciones, seleccione Geometría y luego
doble clic en $X, la función aparece escrita en el panel de la izquierda y en la parte inferior de la tabla
podemos ver la vista preliminar del resultado. Se da clic sobre el botón Ok y se crea el campo. Repita
el procedimiento utilizando la función $Y para la coordenada Norte
29
¿Porque se utiliza esta función?
El signo $ significa que la función se ejecutara sobre las entidades de la capa (sobre la cual se está
trabajando), por lo cual $X devuelve la coordenada X de cada entidad de la capa.
2.- Calcular la longitud de un línea: En este caso en la Calculadora de Campos selecciona la función
$lenght y le coloca un nombre apropiado al campo, por ejemplo: Longitud.
30
Nota: también es posible consultar las propiedades geométricas de una entidad sin tener que realizar
una operación en la tabla de atributos, sencillamente utilizando el botón identificar se despliega una
caja de dialogo que nos muestra los atributos, luego clic en la pestaña “derivados”.
31
AGREGAR CAPAS DE GOOGLE MAPS EN QGIS
3
Es verdad que con el lanzamiento de QGIS 3 estamos contentos. Puesto que trae muchas
características y actualizaciones como la edición mejorada de geometría, símbolos, diseño, capa 3D,
etc. Pero, un problema para los usuarios impacientes que no nos gusta esperar, es la actualización
de los plugins. Aunque muchos ya están en QGIS 3, pero desafortunadamente no todos, incluyendo
uno de mis favoritos OpenLayers, que puede ser utilizado para añadir Google Maps basemap en un
mapa de QGIS como Google Satellite, terreno y mapa de carreteras. ¿Hay alguna otra manera para
que podamos seguir utilizando Google Map en QGIS? La respuesta es Sí y la mostraré en este
tutorial. Con la ventaja es que si se pueden imprimir.
Google Maps Layer está alojando en algún lugar del servidor de Google y envía la ficha al usuario
que lo solicite. Técnicamente se llama Tile Map Service (TMS). Por lo tanto, sólo tenemos que
encontrar el TMS que utiliza Google para usar las capas de Google Maps. Los Google Maps Layer
TMS que utiliza NextGIS se enumeran a continuación:
Para añadir el servicio de mapas de Google Maps en QGIS 3 simplemente, dirigirse dentro del panel
Navegador a XYZ Tiles. Ahora haga clic en una nueva conexión (New Connection), con el botón
derecho del ratón haga clic en XYZ Tiles.
Indique el nombre, por ejemplo: Google Satellite. En la URL copia y pega uno de los TMS de Google
Maps que se listan arriba. Ajuste el nivel máximo de zoom (19). Acepta todos los cambios.
32
Ahora ya se cuenta con la conexión de Google Satellite en XYZ Tile de QGIS 3. Para visualizar
dentro QGIS simplemente hacer doble clic o arrastrar el servicio XYZ Tile al panel de capas.
33
De la misma manera se puede agregar el resto de TMS descritos al inicio de este tutorial, por
ejemplo, Google Maps.
Con todo, para mayor facilidad del procedimiento a continuación comparto el vídeo.
34
CREANDO UN BUFFER DINÁMICO EN QGIS 3
En SIG el termino Buffer se refiere a un análisis/consulta de área de influencia, en capas vectoriales
se utiliza para definir el área a una distancia fija de una entidad, algunos ejemplos prácticos:
determinar las viviendas a una distancia de 200 metros de un hidrante, representar la zona de
protección de un cauce según la normativa. Los Buffers pueden crearse con múltiples anillos e
incluso a un solo lado.
Buffer dinámico
Utilizando las opciones de Simbología es posible crear un Buffer dinámico en QGIS 3, es decir, un
área de influencia que cambie bajo ciertas condiciones.
Ejemplo practico:
Tipo de análisis/simbología: representar un Buffer dinámico alrededor de un curso de agua.
Capas utilizadas: una capa de líneas que representa el cauce del río y una capa de polígonos
Condiciones:
• El espesor del Buffer será 40 veces menor a la longitud del cauce. De esta manera se expandirá cuando se
alargue la entidad y se retraerá en caso contrario.
• El Buffer se creara cuando el cauce tenga una longitud mayor de 1000 metros.
• Resalte la porción del buffer que intersecte un polígono especifico.
Procedimiento:
1.- En QGIS 3 cargar las capas:
35
Nota: las capas deben tener el mismo sistema de coordenadas, los análisis se realizan en el sistema de
proyección de la capa.
2.- En la figura puede observar que la capa de Polígonos llamada “PropiedadB“, tiene una única
entidad, al identificar dicho polígono, cuenta con tres campos, el campo “Nombre” tiene como atributo
“Finca Azusal”.
3.- En el “Panel estilo de capas”, seleccione la capa de línea que representa el curso de agua (en el
ejemplo “hidroB“), luego en la pestaña Simbología, se mantiene el estilo por defecto, Símbolo Sencillo
(Single Symbol), luego clic al botón con icono de cruz de color verde para adicionar una nueva capa
de símbolo.
36
5.- En la capa de símbolo recién creada y en la pestaña “Tipo de Capa del Símbolo” seleccione
“Geometry Generator”, luego en la pestaña “Tipo de Geometría” seleccione “Poligono/Multipoligono”,
esta es la entidad que debe generar el código que escribiremos más adelante.
6.- Debajo de la pestaña “Tipo de Geometría” hay una caja de texto donde escribir el código que
definirá el comportamiento de nuestra Capa de Símbolo, clic sobre el botón con una “E” (símbolo de
Épsilon), esto abre el Dialogo de Expresiones.
Escriba el siguiente código:
8.- Para que el buffer reaccione al interceptar el polígono verde que se encuentra al sur del cauce,
seleccionamos “Relleno” o “Simple fill” de la capa “Geometry Generator” recién creada, luego repetimos
el proceso anterior, dando clic al botón con el signo en cruz de color verde creamos una nueva Capa
de Símbolos, luego “Geometry Generator” > “Tipo de Geometría” > “Poligono/Multipoligono”.
Despliegue el dialogo de Expresiones y coloque el siguiente código:
intersection($geometry, geometry(get_feature(
‘PropiedadB_62be45a2_718e_46d9_855f_afd886692054’, ‘Nombre’, ‘Finca Azusal’ )))
37
Explicación del código:
• intersection(geometria, geometria): determina la intersección de dos geometrías pasadas como
argumentos.
• geometry(get_feature(Layer, Campo, Atributo): devuelve una geometría al especificar la capa (Layer), el
campo y el atributo que le corresponde. El nombre de la capa se obtiene del Dialogo de Expresiones, en
“Capas de Mapas”, en el caso del ejercicio se hace doble clic sobre PropiedadB y se escribe el código que
está en la expresión
•
38
El resultado se muestra en la imagen animada:
39
EPIDEMIOLOGÍA ESPACIAL CON QGIS 3
En esta oportunidad, se utiliza el análisis espacial para evaluar una epidemia de cólera, recrearemos
el caso de John Snow (1854), Físico Londinense padre de la epidemiología.
Datos:
Extracto del Mapa de Londres de 1984, con la ubicación de las muertes por cólera, disponible
enWikipedia. En este enlace puedes descargar los datos para realizar el ejercicio.
El mapa presenta una barra de escala gráfica en Pies, realizando las transformaciones a metros y
colocando unas coordenadas arbitrarias se referencio a escala.
Luego, se trazaron 558 puntos de casos de cólera (estrellas de color rojo), también los sitios de
tomas de agua (capa “Pozos”, círculos de color verde).
40
1.- Evaluando la distribución de los casos:
La primera interrogante es: ¿los casos ocurren de forma aleatoria o siguen un patrón espacial?
Para empezar, ejecute en menú Vectorial > Herramientas de Análisis > Análisis de Vecinos más Cercanos
Posteriormente, en la caja de dialogo se introduce como Capa de entrada la capa que contiene la
ubicación de los casos. El Panel Result Viewer muestra el resultado, un archivo html.
Interpretación:
El análisis de Vecino más cercano, compara la distribución con una distribución aleatoria normal
(Gaussiana), constituye una prueba de hipótesis.
41
Devuelve los siguientes parámetros
• Índice del Vecino más cercano: si el valor es 1 la distribución es aleatoria, valor mayor de 1 es dispersa,
menor de 1 es agrupada (concentrada en ciertas zonas).
Como el valor obtenido es de 0,647479, el patrón es concentrado.
• Z-Score: -15,93059 valor de significancia estadística, define que tan alejados estamos del centro de la
distribución normal, valores mayores a 2,56 y menores de -2,56 rechazan la hipótesis de aleatoriedad.
• Distancia media observada: 961998 distancia media entre los puntos.
• Distancia media esperada: 7.663555 distancia media si los valores tuvieran una distribución aleatoria.
•
Por consiguiente, para demostrarlo ubico los casos y los pozos que abastecían la zona, en total, se
encuentran 11 pozos
42
2.2- Análisis de densidad de casos por pozo
Para realizar el conteo de cuantos casos ocurrieron alrededor de cada pozo, diríjase a menú Vectorial
> Herramientas de Análisis > Contar Puntos en Polígono:
Produce un duplicado de la capa de Thiessen, con un campo denominado Numpoints con el total de
casos de cólera.
Luego, se representa en color degradado en función de la relación del número de puntos entre el
área ocupada:
43
El resultado es evidente, el área de influencia del pozo N° 8 concentra el 60,57% de los casos.
Seleccione menú Procesos > Caja de Herramientas luego en la Caja de Herramientas de Procesos >
Análisis Vector> Distancia al eje más Próximo (Lineas).
Posteriormente, en Capa de Puntos de Origen se seleccionan los “Casos” y en Capa de Puntos Destino
“Pozos”
El resultado produce una capa de línea que une cada Caso a un Pozo.
44
Al seleccionar las líneas unidas al pozo N° 8 obtenemos el mismo resultado del análisis
anterior, 60,57% de los casos.
2.4- Mapa de calor
La herramienta requiere un radio alrededor del cual determinar la cantidad de puntos, luego, asignar
un valor de densidad.
El resultado es un raster de densidad que asigna valores mayores donde hay mayor agrupación.
Finalmente, se obtiene un resultado más llamativo, en la imagen podemos compararlo con los
polígonos de Thiessen.
45
La Herramienta Mapa de Calor cuenta con opciones de configuración, especialmente en Opciones
Avanzadas, las cuales, se explicaran en posteriores artículos.
Este es un ejercicio muy práctico para el análisis espacial, también disponible el proceso en Trainings
Manual de QGIS.
Conclusión:
Se comprueban los resultados de John Snow, existió una evidente relación entre el pozo N° 8 y la
epidemia de Coléra.
46
VISUALIZACIÓN 3D EN QGIS 3
QGIS 3.0 viene con vista de capas 3D que nos permite visualizar datos SIG en 3D. De esta manera
podemos obtener una visualización más realista de los datos que contienen información de elevación
o altura como el Modelo Digital de Elevación (DEM). Es realmente genial visualizar los datos DEM
en 3D para ver la diferencia de elevación del terreno de la superficie terrestre, en lugar de verlo en
una superficie plana.
2.- Los datos contienen información de elevación en cada píxel. Se representan en una sola banda
gris. Pero se puede mejorar su apariencia. Para ello, abra las propiedades de la capa DEM. En Style
Option (opciones de estilo) cambia el tipo de renderizado a Singleband pseudocolor. A continuación,
seleccione una rampa de color.
47
48
3.- Para visualizar los datos DEM en 3D. Lo hacemos en 3D Map View. Puede encontrarla en el
menú View (Ver).
4.- De repente podría ser que se muestre un error, simplemente aceptamos. Esto se debe a que la
vista de mapa 3D no es compatible con el sistema de coordenadas no proyectado. En este caso, los
datos DEM tienen un sistema de coordenadas geográficas. No tenemos que cambiar el sistema de
coordenadas de datos DEM, sino sólo cambiar el sistema de coordenadas de la vista del mapa. Para
ello, haga clic en la ventana QGIS inferior derecha con el CRS actual.
49
6.- Después de cambiar el sistema de coordenadas. Vuelva a abrir la nueva vista de mapa 3D. La
ventana de visualización del mapa en 3D aparecerá plana, no se ve ninguna vista en 3D. No se
preocupe, intente rotar los datos DEM con el botón de aplazar en el teclado + clic izquierdo del ratón.
Mueva el ratón y verá que los datos DEM rotan.
50
7.- Parece prometedor. Pero aun así no se puede ver una buena superficie del terreno. Para ello,
haga clic en el botón Configure (configurar). La ventana de configuración aparecerá, luego fije la
Elevación a los datos DEM y cambie la Escala vertical a un cierto número, por ejemplo 2. La escala
vertical es una escala para dar elevación al terreno. Una escala más alta elevará más el terreno.
Ahora se puede jugar con los controles de zoom y giros.
51
8.- Ahora puede agregar una capa de imágenes satelitales. Para el ejercicio se ha añadido la capa
de imágenes satelitales mundiales de ESRI web map server. QGIS superpondrá automáticamente
las imágenes satelitales por encima de los datos DEM. Finalmente disfrutar de la visualización 3D
en QGIS 3.
52
ETIQUETADO INTELIGENTE DE PUNTOS EN
QGIS 3
El manejo eficiente de capas con grandes cantidades de entidades presenta dificultades requiriendo
el uso de técnicas de etiquetado inteligente. Una ingente cantidad de entidades satura y ralentiza el
despliegue en pantalla.
QGIS 3 cuenta con un conjunto de herramientas que permiten al usuario generar mapas impresos y
en pantalla, comprensibles, atractivos, estéticos y que cumplan normas elementales de cartografía.
Caso práctico: despliegue y etiquetado de puntos de una sísmica:
Se cuenta con una capa shapefile que contiene 30.154 puntos procedentes de sísmica.
2.- Para mejorar la visualización se regulará la cantidad de puntos visibles de acuerdo a la escala,
en los siguientes rangos:
53
• Escalas < 1:60.000 se mostraran cada 25 puntos
• Escala entre 1:60.000 y 1:25.000, se mostrarán cada 20 puntos
• Escala entre 1:25.000 y 1:10.000, se mostrarán cada 5 puntos
• Escalas mayores a 1:10.000 se muestran todos los puntos
3.- Clic botón derecho sobre la capa de puntos, luego Propiedades > Simbología > Basado en
Reglas (Rule-Based). Para adicionar la primera regla haga clic sobre el botón con signo de suma
ubicado en la parte inferior izquierda.
4.- En Etiqueta se da un nombre a la primera regla: “Cada 50”, en Filtro seleccione la letra E para
abrir el Dialogo de expresiones y se escribe la expresión:$id % 25 = 0
54
Explicación de la expresión:
La expresión itera sobre cada uno de los elementos de la capa verificando si cumplen con la
condición, en caso de ser verdadero la entidad será visible
• $id: variable del sistema que nos arroja el “id” de la entidad actual. El id es un campo que asigna un numero
secuencial a cada entidad de una capa de forma automática.
• %: operador matemático que devuelve el resto de una división.
• $id % 25 = 0: los números múltiplos de 25 son los que cumplen la condición.
5.- Se definen los rangos de escala, adicionalmente puede configurar las propiedades de estilo del
marcador. Se repite el proceso para los demás rangos:
55
8.- Se da un nombre a la Regla (“Cada 200”), en la pestaña Filtro se escribe la Expresión: $id % 200
= 0, la cual solo mostrara las etiquetas de los puntos cuyos “id” sean múltiplos de 200. Se define el
intervalo de escala correspondiente (< de 1:60.000).
9.- Se añaden las reglas para cada 50, 10 y 5 puntos.
10.- En la siguiente imagen se puede observar el resultado.
56
ANÁLISIS ESPACIAL INTERACTIVO EN QGIS 3
QGIS 3 permite la integración de los estilos (simbología) con el análisis espacial, las capacidades
son amplias, en esta oportunidad se muestra un análisis espacial con simbología sencilla, posteriores
artículos profundizan en el tema.
Tipo de análisis/simbología: línea que represente el trayecto más corto de una entidad de línea a un
punto específico.
Procedimiento:
2.- Active el “Panel de estilo de capas”, esta ventana permite administrar las propiedades de todas las
capas. Específicamente, trabajara sobre una capa de línea que representa una ruta dentro de una
ciudad, el objetivo es graficar una línea que represente el trayecto más corto a un punto cuyas
coordenadas especificaremos. En el ejemplo se ubico un punto en pantalla con las
coordenadas Este: 465823 y Norte: 1112780, acorde a la Proyección UTM Zona 19 (SRC 2202).
57
3.- Seleccionando la capa de línea que representa la ruta, se configura el estilo de simbología
como Símbolo Sencillo (Single Symbol), y se adiciona una capa de símbolo al dar clic al botón con un
ícono tipo cruz de color verde.
58
4.- En la capa de símbolo recién creada, en la pestaña “Tipo de Capa del Símbolo” elije “Geometry
Generator”.
Luego en la pestaña “Tipo de Geometría” seleccione “Cadena de Línea/de Multilinea”, este es el tipo de
entidad que debe producir el código que escribirá más adelante.
59
5.- Debajo de la pestaña “Tipo de Geometría” se encuentra una caja de texto donde escribir el código
que definirá el comportamiento de la Capa de Símbolo, haga clic sobre el botón con una “E” (símbolo
de Épsilon), esto abrirá el Dialogo de Expresiones.
60
Si edita la linea y esta seleccionada la opción “Actualización en Vivo” la linea que representa el trayecto
más corto se actualizara automáticamente, si crea nuevas lineas tomaran el mismo comportamiento.
Esta es una introducción, podemos hacer otras cosas, por ejemplo determinar el trayecto a una
entidad de una capa con un atributo dado, resaltar los puntos de unión, o crear una etiqueta que
muestre la longitud del trayecto.
61
PERFIL ESPECTRAL DE IMÁGENES DE
SATÉLITE EN QGIS 3
Una de las mayores dificultades de generar un mapa de uso de la tierra es delimitar las diferentes
coberturas. Establecer límites entre tipos de uso o vegetación es un proceso muchas veces subjetivo.
Existen diversas herramientas que pueden ayudar en el proceso de interpretación entre ellas realizar
un perfil espectral.
Perfil espectral
Un perfil espectral es la representación del comportamiento espectral a lo largo de un trayecto.
Por lo tanto, es un proceso similar a crear un perfil topográfico, en vez de la elevación del terreno se
grafican los valores de los pixeles.
QGIS 3 no cuenta con herramientas para crear un perfil espectral, para ello podemos utilizar el
complemento Profile Terrain
Procedimiento:
62
4.- Seleccione la imagen en el Panel de Capas, en la caja de dialogo de Profile Terrain clic en Add
Layer, aparece un cuadro de dialogo que permite adicionar las bandas.
En la imagen del ejemplo se utiliza una imagen Spot5 que cuenta con cuatro bandas, rojo, verde,
Infrarrojo cercano (IR) Infrarrojo Medio (IM), se asigna un color a cada banda.
5.- Seguidamente, trace una línea sobre la imagen y automáticamente se generan las curvas
espectrales del trayecto.
63
En consecuencia, el perfil muestra un cambio drástico de las bandas del IC e IM (color violeta y rojo
oscuro) al pasar el trayecto por un cuerpo de agua
64
El perfil espectral muestra el alza de los valores en el IC indicativo de vegetación alta y los valores
más bajos en zonas con suelo sin vegetación.
65
La imagen muestra el comportamiento del NDVI en la zona evaluada anteriormente.
Al colocar el cursor sobre el perfil se muestra el valor del NDVI, se trazan ejes paralelos que facilitan
la comparación y representa dicho punto en el mapa.
66
En consecuencia, este procedimiento puede ayudar a evaluar zonas aparentemente homogéneas,
como: zonas boscosas o áreas de cultivo:
Por ejemplo, en la imagen se puede notar el comportamiento del NDVI en una zona de cultivo, el
punto más bajo refleja la presencia del canal de riego.
67
PREPARAR UN MAPA PARA IMPRIMIR O
EXPORTAR EN QGIS 3
Generalmente el trabajo en un SIG tiene como fin producir uno o más mapas impresos
(composiciones), cumpliendo los criterios del proyecto.
Confeccionar un mapa es una ciencia y un arte, es un proceso que debe cumplir reglas cartográficas,
criterios estéticos y de diseño.
QGIS 3 nos ofrece varias herramientas para generar mapas con estilo profesional.
Pasos para confeccionar un mapa para impresión en QGIS 3
Preparando la Vista del Mapa
La organización y configuración de las capas en la Vista del Mapa, facilitara la preparación del mapa
para impresión.
Algunas recomendaciones:
• Agrupar las capas que pueden formar parte de la misma categoría en la leyenda, por ejemplo agrupar en
Vialidad las capas: caminos, ferrocarril.
• Configurar adecuadamente los SRC de las capas y el proyecto.
• Renombrar las capas con un título adecuado para la leyenda.
• Remover las capas innecesarias.
• Utilizar mapas bases para mejorar la estética de nuestro mapa.
Mapas Bases:
Utilizando el PlugIn QuickMapServices u OpenLayers, contando con conexión de internet, se
puede incluir una capa base que nos dará un contexto y permitirá prescindir de otras capas.
• Diríjase al menú Proyecto > New Print Layout; o simplemente utilice la combinación de teclas Ctrl +
P.
• Utilizando el asistente, menú Proyecto > Layout Manager despliega una caja de dialogo donde se puede
administrar las composiciones de mapa, de clic en Crear.
68
En ambas opciones aparece un pequeño dialogo solicitando el nombre de la composición.
En menú Diseño > Page Setup podemos definir las dimensiones que tendrá nuestro mapa
impreso.
Adicionando la Vista del Mapa
El diseño de un mapa consiste en insertar elementos sobre la composición y administrar sus
propiedades, tales como la vista del mapa, barra de escala, tablas, entre otros.
Para añadir la Vista del Mapa, en menú Añadir Elemento > Añadir Mapa, o en el botón con el
mismo ícono en la barra de herramientas.
69
Posteriormente, manteniendo presionado el botón del mouse, luego arrastrándolo define la posición
y extensión que tendrá el mapa en el papel.
70
Definiendo la escala y sistema de coordenadas
Si queremos ajustar nuestra vista a una escala precisa, podemos editar las propiedades del
elemento, sencillamente introducimos el valor de escala.
Seguidamente, puede configurar el sistema de coordenadas en SRC, en el ejemplo se utilizó
UTM/WGS84 Zona 17 Norte.
Creando una ubicación relativa:
Una ubicación relativa es una vista a menor escala del área de estudio, el objetivo es transmitir al
lector el contexto nacional y regional. Para crearla:
71
Nota: con esta opción no es necesario bloquear el Mapa central, cada Mapa puede tener diferente
posición y escala
Rejilla de Coordenadas:
1.- Seleccione el elemento Mapa.
2.- Luego, en el Panel Propiedades del Elemento, expanda las opciones, seleccione Cuadricula.
3.- Para crear una cuadricula, de clic al botón con el signo de sumar de color verde.
72
4.- Para configurar la cuadricula QGIS 3 ha añadido el botón “Modificar Cuadricula”, seleccione y
se despliegan las opciones.
5.- Configure el espaciamiento en los ejes X y Y, el tipo de línea, letra, decimales, posición de las
coordenadas y demás parámetros de la rejilla de coordenadas.
73
Finalmente, en la imagen se muestra el resultado.
74
4.- Incremente la Anchura de la Cabeza de Flecha.
Leyenda:
1.- En menú Añadir Elemento> Añadir Leyenda.
2.- Se dibuja un rectángulo donde queremos ubicar la leyenda.
3.- En el Panel de Propiedades desactiva la opción Auto Actualizar y luego edita la leyenda,
dando clic al signo menos de color rojo, remueva las capas innecesarias.
75
4.- Utilizando el botón con forma de lápiz podemos editar el nombre de la capa o categoría.
5.- En la opción Columnas podemos optimizar el espacio ocupado por la leyenda dividiéndola en
columnas.
76
Barra de escala menú Añadir elementos> Añadir Barra de Escala.
De igual forma textos y gráficos.
Finalmente, el menú Diseño tiene las opciones para imprimir y/o exportar el mapa en diferentes
formatos, como imagen, PDF, SVG.
77
OPERACIONES ESENCIALES EN LA
CALCULADORA RASTER DE QGIS 3
La Calculadora raster de QGIS 3 es un módulo que permite realizar múltiples operaciones sobre
archivos raster, entre ellas matemáticas, lógicas, reclasificaciones y reasignaciones.
Existen diversos tipos de archivos raster, producto de interpolación, rasterizaciones, modelos digitales
de elevación (MDE) e incluso imágenes de satélite.
QGIS 3 permite realizar operaciones sobre archivos raster de múltiples formatos, en el ejemplo se
demuestra la aplicación sobre archivos grd (Surfer), Geotiff (imagen landsat) y hgt (datos de
elevación Mision SRTM)
Calculadora Raster
Para acceder a la Calculadora Raster seleccione menú Raster > Calculadora Raster.
En la parte superior izquierda, en “Bandas Raster” aparecen las capas raster previamente cargadas
en el mapa, las cuales serán insumos para realizar los diferentes procesos. En la parte superior
derecha están los parámetros para configurar el archivo de salida, nombre, sistema de referencia
cartográfico, extensión, entre otros
78
2.- Diferencia entre precipitación anual y evaporación anual
En el siguiente caso se cuenta con un raster de precipitación y otro de evaporación anual, generados
de la interpolación de estaciones meteorológicas. La diferencia entre ellos es un indicativo del tipo
de clima.
En la imagen se puede observar que el raster de precipitación (color verde-azul) cubre mayor
extensión que el de evaporación (rojo-blanco), el número de estaciones del primero (círculos de color
blanco) es mayor y más amplio.
79
En este caso se escribe la ecuación: "PrepA@1" - "EvapA@1"
3.- Condicionales
Es posible extraer zonas que cumplen con dos o más condiciones, por ejemplo: precipitaciones >
1000 mm y elevaciones >= a 1000 m.s.n.m.
80
Aplicación en imágenes de satélite (Raster multibanda)
Sin embargo, de igual forma como se procesan raster productos de interpolación o rasterización, se
pueden realizar sobre archivos de imágenes multibanda.
Previamente cargada una imagen multibanda la calculadora lista cada banda por separado
indicándola con un @
81
2.- Condicional. Extrayendo un espejo de agua
Como resultado, en el NDVI los valores menores a uno (01) corresponden a los cuerpos de agua,
realizamos un condicional para seleccionarlos: "ndvi@1" < 0
Para extraer y calcular la superficie del espejo de agua podemos convertir el raster booleano a
vectorial en:
82
RECLASIFICACIÓN DE UN RÁSTER EN QGIS 3
Las operaciones de reasignación o reclasificación, consisten en la modificación total o parcial de los
valores de un ráster. La reclasificación permite agrupar en grupos los valores actuales de un ráster.
Estas operaciones pueden ser de utilidad en análisis de los diferentes tipos de ráster, Modelos
Digitales de Elevación, ráster discretos, imágenes de satélite.
Ejemplos en ráster de variables continúas:
Denominamos raster continuos aquellos que representan una variable continua en el espacio como
elevación, precipitación, ruido, entre otros, también comprende imágenes de satélite o producto
generados.
83
En el menú Procesos active la Caja de Herramientas, luego ubique GRASS > r.reclass, en la caja de
dialogo se selecciona el MDE, en la caja de texto Reclass Rules colocamos la expresión que define
los intervalos:
• 0 thru 100 = 1
• 100 thru 200 = 2
• 200 thru 300 = 3
• 300 thru 400 = 4
• 400 thru 600 = 5
• 600 thru 1000 = 6
• 1000 thru 1600 = 7
• 1600 thru 2000 = 8
84
Se ejecuta la herramienta y genera el ráster de salida, que por defecto tiene el
nombre Reclassified. Este procedimiento facilitara la cuantificación de la superficie de cada clase y la
representación cartográfica.
85
Para ello se introduce en la Calculadora Raster la expresión:
("Pendiente@1"<= 25) * 25 + (("Pendiente@1" > 25) AND ("Pendiente@1" <=
50)) * 50 + (("Pendiente@1"> 50) AND ("Pendiente@1" <= 75)) * 75 +
("Pendiente@1" > 75) * 100
Explicación: el operador de * asigna un valor a los pixels que cumplen la condición
Como el ráster generado cuenta con solo 4 categorías se puede aplicar un estilo por Único Valor.
86
Generando estadísticas: para cuantificar la superficie ocupada por cada rango de pendiente
seleccione:
Caja de Herramientas de Procesos > Análisis Raster > Informe de valores únicos capa raster.
87
Los valores del NDVI son relacionados con la presencia o ausencia de vegetación, valores bajos
generalmente coinciden con suelos desnudos o cuerpos de agua.
Comparando los valores del NDVI con la imagen de satélite se definen los rangos
• < 0 agua
• 0 a 0,3 suelos desnudos zonas urbanas
• 0,3 a 0,65 herbazales cultivos
• 0,65 a 1 cultivos densos bosques
Los valores del NDVI son de tipo punto flotante (Float) con 14 decimales, la Calculadora Raster tiene
la limitación que no realiza transformaciones a enteros.
Antes de realizar la reclasificación se multiplica el NDVI por una constante, en este caso se utilizó
1000:
88
Finalmente configuramos el estilo y calculamos las estadísticas
89
VISUALIZACIÓN AVANZADA EN 3D USANDO
QGIS
QGIS es una aplicación muy poderosa en lo que a SIG libre se refiere. Una de sus potencialidades
es la posibilidad de desarrollar complementos (plugins) que permitan realizar tareas específicas,
como vistas 3D. Estos plugins son desarrollados por la comunidad y son herramientas flexibles que
permiten incrementar las de ya por sí elevadas capacidades de QGIS.
Si bien es cierto que QGIS 3.0 incluye soporte directo para visualización 3D, uno de los plugins más
completos para visualización tridimensional es Qgis2threejs desarrollado por Minoru Akagi. Este
plugin permite exportar datos de terreno (elevación), imagen y datos vectoriales al navegador web.
Los datos exportados se pueden visualizar (e interactuar) en cualquier navegador web con soporte
WebGL.
Requerimientos: QGIS instalado, el plugin qgis2threejs instalado (se puede obtener del sitio web
indicado o a través del gestor de complementos de QGIS), navegador web con soporte WebGL (Yo
uso Chrome u Opera, aunque funciona en Firefox y Edge).
Desarrollo del ejercicio
En QGIS he cargado 3 capas vectoriales y 2 capas raster. Las capas vectoriales corresponden a el
límite de una cuenca (rojo), una red hidrográfica pequeña (azul) y unos puntos que representan
estaciones pluviométricas en la región de estudio (verde). Las capas raster corresponden a un
modelo digital de elevación y una ortofoto.
90
Instalado el plugin aparece en la barra de herramientas de QGIS . Al hacer clic sobre el botón
se abre el cuadro de diálogo del complemento.
Los valores de altura en la capa MDE son los que se representarán en la salida del complemento,
es decir es la información que dará la percepción de relieve, cualquier error o artefacto en esta capa
generará errores en la visualización de la misma.
En las opciones de tipo de visualización (Display type), podemos escoger 4 opciones, la primera
(Map canvas image) desplegará sobre el relieve generado por el MDE seleccionado, lo que en ese
momento tengamos desplegado en nuestro visor de QGIS. Por esta razón es importante planificar
de antemano que queremos representar. La segunda opción despliega sobre el relieve una capa que
seleccionamos, la tercera opción un archivo de imagen y la cuarta opción únicamente despliega el
relieve con un color sólido escogido.
Una vez seleccionadas las opciones en el cuadro de diálogo, hacer clic en Run. Se abrirá una nueva
pestaña en el navegador predeterminado, en el cual aparecerá la visualización 3D de nuestro
modelo. Adicionalmente se puede guardar el archivo html que genera qgis2threejs para poder
acceder al mismo posteriormente.
91
En las siguientes imágenes se puede observar la diferencia entre la visualización únicamente del
MDE y las capas vectoriales y en la segunda imagen el mismo relieve, pero recubierto con la ortofoto
de la zona de estudio.
92
Opciones avanzadas de Qgis2threejs
El procedimiento anterior, se podría decir es la base del plugin. Sin embargo, sus capacidades son
un poco mayores. El plugin ofrece la capacidad de personalizar las capas vectoriales para su
representación en el modelo tridimensional. En la imagen siguiente, al escoger la capa de estaciones
se puede personalizar la representación de los símbolos de dicha capa en la visualización 3D, en
este caso, se escoge cilindros para representar las estaciones, los mismos que serán de color rojo.
93
Adicionalmente se puede configurar la exageración vertical de la visualización, herramienta que es
particularmente útil en áreas en donde el contraste de alturas no es muy grande y se necesita realzar
ciertas características del terreno.
94
La visualización resultante será:
En la que claramente se puede apreciar que el relieve es algo exagerado. Sin embargo, permite una
mejor visualización de las formas del terreno.
Conclusiones y recomendaciones
95
EVALUANDO PARÁMETROS DE SUELOS
INTERPOLADOS EN QGIS 3
Una de las tareas más llamativas al trabajar con SIG es representar variables biofísicas en 2D/3D,
integrarlas y analizarlas
Generalmente, se parte de datos puntuales recolectados utilizando un sensor, luego son
espacializados mediante métodos de interpolación.
Seguramente han asistido a una conferencia donde el ponente enfoca con un apuntador sobre una
serie de mapas y señala “noten que tan similar son los mapas”, pero, ¿Qué determina la similaridad
entre ellos? (adaptado de Berry, 1999).
El problema
¿Que tan similar o diferente es un mapa de otro?,¿es confiable emitir conclusiones solo con el
análisis visual de los mapas raster?
96
Caso de estudio:
Existen varios escenarios y los métodos aplicables pueden ser distintos:
En este caso veremos algunas estrategias para evaluar cualitativa y cuantitativamente variables de
una zona, específicamente datos interpolados de un estudio de suelos.
Datos disponibles
Para este ejercicio utilizaremos un conjunto de puntos de un análisis de suelos de una zona de
cultivo, denominada SoilSamp, los datos son tomados del tutorial del programa ArcView3.2 (descargar
datos del ejercicio).
Ademas, se cuenta con el modelo digital de terreno del área de estudio
En total consisten en 74 puntos con análisis de los parámetros Potasio (K), Fosforo (P), Biomasa,
pH, Materia orgánica y conductividad eléctrica
97
Estadísticas
QGIS 3 cuenta con herramientas que permiten extraer estadísticas sencillas de las capas (ver el
articulo estadísticas en QGIS 3).
Posteriormente, aprovechar el potencial de la visión geográfica para evaluar las estadísticas por
zonas. Por ejemplo seleccionar los puntos ubicados a menor altura, extraer las estadísticas y
compararlas con el resto del área.
La selección puede realizarse manualmente, para áreas pequeñas y homogéneas, en caso contrario
se recomienda seguir el procedimiento de Reclasificación del MDE > Vectorización > Consulta
espacial > Estadísticas de lo seleccionado.
98
Como en el caso anterior, también, puede graficar puntos seleccionados
El histograma permite ver como se distribuyen los datos y el gráfico de dispersión contrastar dos
variables
En la imagen, el histograma bimodal del Fósforo (Soil_P) y en menor medida del Potasio (Soil_K)
puede indicar la presencia de dos poblaciones, por ejemplo: dos tipos de suelo o tipos de manejo.
Por otro lado, el gráfico de dispersión Potasio/Materia Orgánica presenta la mejor correlación.
99
Analizando las variables interpoladas
Para ello menú Raster > Miscelánea > Construir Raster Virtual, en la caja de dialogo en capas
de entrada seleccione las variables.
El orden en que son seleccionadas serán numeradas y organizadas en el ráster resultante.
Luego, la composición a color (RGB), en la imagen se visualiza Fósforo en el cañón Rojo, Potasio
en el cañón verde y Biomasa azul.
100
Puede realizar la consulta por datos o por grafo, para más detalle consulte el artículo análisis espectral
y multitemporal.
Nota: las variables tienen distintas escalas de valores, es necesario normalizarlas para poder realizar
comparaciones.
En próximos artículos se profundiza el análisis combinando los gráficos con las estadísticas y la
representación espacial, determinando correlaciones, evaluando la precisión de la interpolación.
101
CREACIÓN DE SHAPEFILES EN QGIS 3
QGIS 3 Girona es la versión más reciente de este potente software libre, con notables mejoras en la
administración, visualización, y consulta de archivos. En esta oportunidad se aborda la creación de
archivos vectoriales, resaltando las diferencias con las versiones anteriores.
Archivos vectoriales:
El antiguo Shapefile es el archivo geoespacial de mayor uso en el ámbito SIG, sin embargo, su
sustitución por formatos más versátiles y completos ha sido anunciado hace tiempo, específicamente
por bases de datos espaciales, ESRI ha priorizado el uso de las Geodatabases y ahora QGIS hace lo
propio con SQLite(mediante la extensión Spatial Lite).
En las versiones anteriores de QGIS el archivo vectorial Shapefile era el protagonista seguido
por Spatial Lite (SQLite), la nueva versión de QGIS marca un cambio elevando el estatus dentro de
su interfaz al formato de bases de datos Geopackage (implementación de Spatial Lite). Es importante
señalar que igualmente podemos crear capas Spatial Lite en QGIS 3.
Las bases de datos espaciales superan las limitaciones del Shapefile al permitir almacenar y
gestionar en un solo archivo tablas, datos Raster y Vectoriales. La comparación Spatial
Lite y Geopakage es un tema amplio que abordaremos en otro momento.
¿Cómo crear archivos Shapefile en QGIS 3?
Existen varias alternativas para crear archivos Shapefiles, se puede usar la barra de menús
siguiendo la siguiente secuencia:
Menú Capa > Crear capa > Nueva capa de archivo shape
Otra forma es utilizando las barras de herramientas, en este caso el “Administrador de Capas > Nueva
capa de archivo Shape..”
102
La alternativa más directa es mediante barra del “Administrador de fuentes de datos”
Definiendo las características del Shapefile
Independientemente del modo que seleccionado, se despliega una caja de diálogo en la cual se
puede configurar cada uno de sus parámetros.
• Nombre de archivo.- Define el nombre y ruta donde se almacena en nuevo fichero shapefile.
• Codificación del archivo.- Por defecto aparece System o UTF-8, se recomienda mantenerla para no tener
problemas con incompatibilidad de caracteres en las tablas de atributos.
• Tipo de geometría.- Es importante porque permite definir la nueva capa vectorial, sea punto, línea o
polígono.
• Se tiene la opción de incluir la tercera dimensión “Z”, con clic a esta opción puede añadir y editar los
valores de altura de los vértices de las entidades.
• También existe la opción para añadir valores “M”, útiles para la geocodificación.
• Para definir el sistema de coordenadas geográficas o proyectadas, en el botón Seleccionar SRC se elige
entre los sistemas utilizados recientemente, o con clic sobre el mismo para buscar y seleccionar el sistema
de preferencia.
• En esta misma caja de dialogo se puede diseñar la tabla de atributos, es decir, definir los campos. Por
defecto se crea la tabla con solo el campo por defecto Id, en caso de no agregar ningún otro campo.
• Es posible crear cuatro tipos de campos de acuerdo al tipo de datos: texto, números enteros, números
decimales y datos de fecha.
• Clic en Aceptar / OK para finalizar con la creación del nuevo shapefile.
• En las versiones QGGIS 2.x se puede crear una capa Shapefile mediante la combinación de las teclas
Ctrl+Shift+N, en la versión QGIS 3 la misma secuencia permite crear un Geopakage.
• QGIS 3 incorpora los botones para crear Shapefiles y Geopakages en la barra de “Administrador de fuentes
de datos”
• QGIS 3 da opciones para asignar las dimensiones Z y M.
103
DIBUJANDO CON PRECISIÓN EN QGIS
En esta oportunidad continuamos mostrándoles más herramientas de dibujo preciso en QGIS, para
ello utilizaremos el PlugIn QAD, en el cual llevaremos a cabo tres tareas comunes en AutoCAD como
son:
• Extender una línea hasta el borde de otra entidad
• Cortar o dividir una línea en un punto específico.
• Modificar la longitud de una línea
1.- Extender una línea
Una de las operaciones más comunes en AutoDesk es la posibilidad de extender un segmento de
línea hasta el borde de una entidad cercana.
¿Qué condiciones necesitamos para que sea realmente útil? Que la herramienta nos permita trabajar
con versatilidad y sin restricciones, por ejemplo, que reconozca elementos de la propia capa y de
otras capas, que trabaje con diferentes tipos de geometrías
Primero instalamos y/o activamos el complemento QAD:
104
Activando la herramienta Alarga.
En la siguiente imagen podemos ver unas polilíneas de color rojo, queremos extender el vértice de
la línea en donde está el círculo azul hasta la que está encerrada en el polígono verde. La capa de
línea que vamos a extender debe estar edición.
Primero seleccionamos (con la herramienta de selección de QGIS) la línea “barrera” (hasta donde
deseamos extender la línea), el siguiente paso es activar la herramienta Alargar o tecleamos
“Alargar” en la línea de comandos, luego hacemos clic sobre el segmento que deseamos extender:
105
Caso 2. Extendiendo segmentos de línea hasta el borde de un polígono
En el siguiente caso extenderemos la línea hasta el borde de un polígono, veremos otras formas de
seleccionar y activar la herramienta
Estando activo el complemento QAD podemos seleccionar el polígono simplemente haciendo clic
sobre él, lo cual cambia el cursor a forma de cruz desplazamos para formar un cuadrado de
selección, al hacer clic queda seleccionado el polígono.
Si hacemos clic botón derecho aparece una caja de dialogo que muestra el comando utilizado
recientemente, seleccionamos “Alargar”
106
Ahora hacemos clic sobre la línea a extender:
Para esta tarea podemos utilizar dos herramientas Dividir y Parte, la primera corta la polilínea en el
punto especificado y en el proceso tiende a remover una parte del segmento (Parte) . La segunda
(Divide) separa el segmento de línea en el punto especificado manteniendo la continuidad del
trazado. En las figuras pueden ver ejemplos de ambos casos: Dividir.
107
Parte:
Consultar la longitud de una línea: al igual que el comando Length de Autocad al utilizar la
herramienta Longitud y hacer clic sobre una línea nos arroja el valor de la longitud en unidades del
sistema de proyección en que este configurado el proyecto, por ejemplo si estamos trabajando en
una proyección UTM nuestras unidades serian metros.
108
Dinámicamente: al activar el comando Longitud aparecen varias opciones en la línea de comandos de
QAD, una de ellas es Dinámica si aceptamos esta opción la línea se proyecta siguiendo su dirección
y podemos modificarla interactivamente.
Incremento: cómo podemos inferir de su nombre esta opción permite incrementar la longitud de una
línea a un valor especificado, es decir, adicionamos (si el valor es positivo) o restamos (si el valor es
negativo) una cantidad especificada a la extensión de la línea.
109
Pueden explorar las otras opciones como son: Porcentaje y Total, esta última permite especificar el
valor de la longitud de la línea, es decir, indiferentemente de la extensión de la línea podemos
especificar su valor, por ejemplo, mide 99 metros y queremos que mida 100 metros exactos.
Existen muchas otras herramientas en este interesante complemento que veremos más adelante
110
Dibujando en QGIS con mayor precisión que
AutoCAD
El título de este artículo es un reto, si has utilizado los programas de Autodesk supongo que podrías
pensar “es una estrategia publicitaria”, “una exageración”, los invito a continuar leyendo y verán como
QGISsiendo software libre cumple el desafió de brindarnos una interfaz intuitiva y profesional para
dibujar con precisión en 2D similar a la de Autodesk. Es más…
¡Si deseas dibujar como lo haces en Autodesk, si vas realizar trazados considerando ángulos y
distancias, o si te has propuesto llevar a cabo un catastro con un SIG libre, tienes que ver esto!!
Para cumplir los objetivos de este artículo nos centraremos en un solo Plugin, utilizando muy poco
de las herramientas de dibujo propias de QGIS.
Es pertinente aclarar que existen variados Plugins que mejoran las capacidades de dibujo de QGIS,
entre ellos puedo nombrar: CADDIGITIZE y CADTOOLS.
Sin embargo,ninguno de ellos es tan versátil y completo como el que veremos a continuación, tema
central del artículo.
Plugin QAD:
El objetivo de este Plugin es proveer a QGIS de una interfaz lo más similar a AutoCAD. Dos aspectos
muy positivos:
Las herramientas de Qad pueden ser accedidas a través de la barra de botones que se activa con
el complemento, tal como se resaltan en la imagen que se muestra seguidamente.
111
O pueden accederlas en submenú “Qad” que aparece en el menú “Vectorial“:
NOTA IMPORTANTE: la linea de comandos es opcional, pueden realizar todo su trabajo prescindiendo de
la misma.
2.- Botón derecho: Otra característica similar a AutoCAD es que haciendo un clic derecho podemos
acceder a opciones de forma intuitiva mientras trabajamos, me explico, mientras trazo una polilínea
con clic derecho aparece una caja de dialogo que entre otras cosas me permite terminar el trazado
(Enter), cerrar la polilínea (Cerrar), escoger entre las coordenadas más recientes para colocar el
siguiente vértice, definir la longitud, trazar un arco, entre otros.
112
Incluso, sin estar dibujando pero aún activas las herramientas Qad, si desplazamos el cursor sobre
la pantalla y hacemos clic izquierdo se despliega una caja de dialogo que permite escoger entre las
herramientas utilizadas recientemente o entrar en opciones para mayores configuraciones.
3.- Modo de trabajo: Qad hace mucho más que brindarnos herramientas de dibujo, crea un nuevo
ambiente de trabajo mientras este activo, eso incluye: altera la forma, color y desplazamiento del
cursor, da las opciones con clic derecho ya comentadas, ajuste de precisión (snap), incluso el
comportamiento de las entidades al seleccionarlas, todo de forma similar a AutoCAD.
113
4.-Ajuste de precisión (Osnap): los comandos de precisión son la característica principal de AutoCAD.
¿Qué es un comando de ajuste de precisión?
Es el que nos permite unir entidades con precisión mientras dibujamos, por ejemplo el vértice de una
línea o polígono, con un punto, directamente a una porción de una línea o en forma perpendicular,
paralela, o unirlo a una esquina de una entidad, entre otros. Los comandos de ajuste de precisión
son denominado “Snapping”
Si seleccionamos el botón anteriormente descrito se nos despliega una caja de dialogo para
configurar los diversos Snap y otros parámetros de dibujo:
Y lo más interesante es que este es solo una pequeña parte del gran conjunto de herramientas y
configuraciones disponibles, tenemos también el modo de rastreo polar, herramientas para acotar,
escalar objetos, simetría, buffer, offset, crear matrices, dividir y extender, y mucho más. Simplemente
es necesario jugar con todas la herramientas disponibles y podrás descubrir la infinidad de cosas se
puede hacer.
114
MÚLTIPLES CANVAS EN QGIS 3
Como sabemos, hace poco ha sido lanzada la nueva versión de QGIS 3 Girona, con muchas
funcionalidades y utilidades que pretenden hacernos más sencillo nuestro trabajo con este programa.
Entre estas funcionalidades, una de las que más me gustan es la incorporación de múltiples vistas
de mapa o canvas.
Aquí veremos que se nos abre una ventana que podremos incrustar entre nuestros paneles.
En este sentido dispondremos de algunas de algunas funciones que nos harán más útil su
comportamiento:
Sincronizar con la ventana principal: Con esta característica veremos que ambos mapas se
moverán cuando movamos la ventana principal.
115
Sincronizar escala: Igual que podemos permitir que se mueva en X e Y, podremos hacerlo en Z,
moviéndose ambos mapas a la vez.
Si utilizamos la segunda opción será muy importante ponerle un factor a la escala para que una de
ellas tenga más zoom que la otra, de lo contrario no sería muy útil tener dos vistas iguales.
Entre otras características, vemos que también podemos rotar el mapa, mostrar las anotaciones,
o incluso mostrar un pequeño recuadro de la extensión de la vista principal.
¿Y tú, sabías de la existencia de esta nueva funcionalidad de QGIS 3? Si quieres formarte en QGIS
3, puedes acceder a nuestro Curso de Introducción a los GIS actualizado de QGIS 3.
116
TUTORIAL GIS: CÓMO CONECTAR POSTGIS
DESDE QGIS
En este tutorial aprenderemos a conectar PostGIS desde QGIS, ambas herramientas líderes en el
sofware libre no comercial.
Dentro del mundo del software libre, PostgreSQL y su extensión espacial (PostGIS) se han
convertido en el sistema gestor de bases de datos (SGDB) más extensamente utilizado.
PostGIS permite el uso de objetos GIS(Geographic information systems) incluyendo soporte para
índices GiST basados en R-Tree, y funciones básicas para el análisis de objetos GIS. Con PostGIS
podemos usar todos los objetos que aparecen en la especificación OpenGIS como puntos, líneas,
polígonos, multilíneas, multipuntos, y colecciones geométricas. Además, en los últimos años ha
desarrollado una extensión ráster, consiguiendo así situarse como una de las bases de datos
espaciales libres más completas y potentes que existen actualmente.
PostGis está implementado utilizando como base PostgreSQL, por lo que hereda de forma
automática todas las características de éste, así mismo PostGIS soporta tipos de datos espaciales,
índices espaciales y tiene cientos de funciones espaciales y permite importar y exportar datos a
través de varias funciones propias. Todo esto favorece el poder conectar PostGIS con QGIS, no
cabe duda.
El uso de PostGIS nos permite realizar múltiples operaciones de análisis sin llevar implícito el uso
de un software SIG, lo que suma un aliciente más al conjunto de innumerables ventajas (software
libre, compatibilidad con estándares de OGC, …), que tiene PostGIS.
Sin embargo otra característica de PostGIS es que es operable con otros software SIG de escritorio
como pueden ser OpenJump, Kosmo o Qgis.
117
QGIS y Postgis – Cómo conectar PostGIS
A pesar de que con PostGIS podemos realizar algunas de las principales tareas que realizamos con
un gis como realizar selecciones complejas, utilizar operadores de extracción o la realización de
buffer, muchas veces es preferible trabajar o editar nuestra información cartográfica desde un
software gis. A continuación mostraremos como es de sencillo conectar PostGIS con Qgis y trabajar
la información almacenada en nuestra base de datos desde QGis (Debemos tener una bases de
datos espacial creada en PostGIS).
118
3 . Creada la conexión procedemos a cargar una nueva capa:
4. Nos aparecerá un nuevo cuadro donde seleccionar las capas que queramos insertar.
5. Seleccionamos y ya tendríamos cargadas las capas en Qgis para poder trabajar con
ellas.
119
Y de esta forma tan sencilla hemos aprendido a conectar PostGIS con QGIS. Si quieres conocer
más de cerca de PostGIS, importar ficheros, y trabajar con diferentes tipos de información, te
invitamos a que visites nuestra web de formación y descubras nuestro nuevo Curso Superior en
Gestión de Bases de Datos Espaciales: PostgreSQL/PostGIS.
120
CÓMO TRABAJAR CON RASTER EN QGIS
(HILLSHADES, CURVAS DE NIVEL, PERFILES
TOPOGRÁFICOS Y 3D)
En este tutorial vamos a explicar paso a paso:
QGIS tiene fantásticas herramientas y complementos para el procesamiento de datos del terreno o
elevaciones y para su visualización. Ante la pregunta:
La respuesta es que no necesitas realizar ningún tipo de conversión previa. QGIS lee perfectamente
los archivos .asc (Arc/Info ASCII Grid) gracias al uso de la librería GDAL para leer y escribir
formatos de datos raster, incluyendo Arc/Info Binary Grid, GeoTIFF, Erdas Imagine y muchos más.
La librería GDAL soporta más de 100 formatos raster.
Índice
• Obtener datos del terreno
• Añadiendo datos raster en QGIS
• Convertir a otro formato raster
• Extraer las curvas de nivel
• Crear un hillshade (mapa de sombras)
• Mostrando datos del terreno en 3D
• Crear perfiles topográficos con QGIS
• Cómo utilizar la calculadora ráster en QGIS
121
Para poder realizar la descarga es necesario registrarse en el sitio. Una vez realizada la
descarga tendremos un archivo llamado MDT200-LEON-H30.asc
Una vez añadida la capa veremos los datos del terreno en QGIS. Cada píxel del raster representa la
elevación media en metros en ese lugar. Los píxeles oscuros representan las zonas de baja altitud
y píxeles más claros representan áreas con alta altitud.
122
Haciendo uso de la librería GDAL, QGIS es capaz de convertir prácticamente cualquier formato en
otro, la lista es muy amplia (.SID, . TIFF, . IMG, HDR.ADF, .DT0, .DT1, .DT2, . JPG, . PNG, .BMP,
.TER, .ASC, .SQLITE, .E00, .MBTILES, etc). Consulta la lista de formatos raster de GDAL.
Si no ves el menú Ráster, debes habilitar el plugin GdalTools desde Complementos > Administrar e
instalar complementos… > Instalado.
Clic en Aceptar.
123
Una vez finalizado el proceso, verá las curvas de nivel cargadas en QGIS. Cada línea de esta capa
representa una elevación. Todos los puntos a lo largo de estas líneas estarían a la misma altura. Y
cuanto más próximas estén las líneas más pronunciada será la pendiente.
Como puedes ver en la imagen superior, la herramienta Raster > Análisis >MDT (Modelos del
terreno) también nos permite crear mapas de:
• Pendientes.
• Orientación.
• Escabrosidad.
124
Jugando con la transparencia del hillshade y cambiando la rampa de colores del MDT original este
es el resultado:
125
MOSTRANDO DATOS DEL TERRENO EN 3D
Basta con introducir el raster que contiene los datos de elevación, el complemento utiliza los valores
de elevaciones y los combina con la imagen de QGIS para crear un archivo HTML. la imagen 3D es
dinámica, podemos ampliar, alejar y movernos por el mapa:
126
CREAR PERFILES TOPOGRÁFICOS CON
QGIS
Para crear perfiles topográficos en QGIS necesitamos utilizar el complemento Profile Tool
Como resultado vemos que no hemos trazado demasiado bien el curso del río.
Para finalizar podemos exportar el perfil a PDF, PNG, SVG o PS.
• En esta entrada puedes ver cómo crear perfiles topográficos con ArcGIS.
127
dichos datos. Los resultados de dichas operaciones se muestran en una nueva capa ráster en un
formato admitido por la librería GDAL.
En este tutorial te mostramos cómo modificar los valores de un ráster a partir de una máscara con la
calculadora ráster de QGIS.
128
CONSTRUYENDO UN MODELO DIGITAL DE
ELEVACIÓN QGIS 3
MDE en QGIS3
Sin embargo, cuando se desea modelar una porción de la superficie de la tierra el problema
es diferente, el relieve posee formas y detalles conocidos que se deben reproducir.
El presente artículo forma parte de una secuencia, posterior a la creación del MDE (descargar
ejercicios de práctica), se incorporará un dique, se estimará el espejo de agua y el calculo del
volumen embalsado.
▪ Puntos de elevación.
▪ Líneas de estructura: líneas continuas que agregan información a la superficie.
▪ Líneas de ruptura : producen cambios marcados en la superficie, por ejemplo para representar
ríos, vialidades, líneas de costa.
129
Caso de estudio:
Generación de un MDE de una pequeña superficie de relieve escarpado dentro de una finca
agropecuaria.
Procedimiento:
Cargue los archivos vectoriales en la vista del mapa.
Puede notar que tenemos entidades de puntos, líneas (hidrografía, curvas de nivel cimas y
laderas) y polígonos (área de estudio).
Para acceder a las herramientas de interpolación diríjase a Menú Procesos > Caja de
Herramientas de Procesos, en el Panel seleccione la pestaña Interpolación,
luego Interpolación TIN.
130
Modelo Digital con puntos de elevación
131
En la pestaña Capa Vector seleccione las capas de puntos, en Atributo del punto
seleccionamos el campo Elevation. Luego clic al botón con el signo de suma de color verde
para agregar la variable.
QGIS 3 ofrece dos métodos para generar el TIN, Lineal y Clough-Toucher (cubica), la
segunda produce una superficie más suave.
En la pestaña Extensión clic en el botón de la esquina > Usar extensión de la capa del lienzo.
132
En la Caja de dialogo Select Extent seleccione la capa Area.
Puede configurar el tamaño del pixel en las opciones Número de Filas y Número de
Columnas, un tamaño de pixel más pequeño genera una superficie mejor definida.
Para más detalles sobre como ajustar las dimensiones del pixel consulte el
articulo Interpolación espacial en QGIS 3.
Finalmente puede definir los archivos de salida.
En Interpolado el nombre y ubicación de la superficie
Luego en Triangulación el vectorial con la red de triángulos.
133
Modelo Digital con puntos y líneas de ruptura
Nuevamente, despliegue la herramienta de interpolación TIN y configure las capas.
Introduzca la capa Puntos de igual forma que el caso anterior.
Luego, seleccione en Pestaña Capa Vector > curvas de nivel y en Atributo el campo
Elevation, adiciónelo a la lista, en Tipo > Lineas de Ruptura.
Seguidamente, para la capa Cimas-Laderas seleccione la opción “Usar coordenadas z para
interpolación”, en Tipo > Línea de Estructura.
Para la capa Hidro seleccione la opción “Usar coordenadas z para interpolación”,
en Tipo > Línea de Ruptura.
134
Configure las demás opciones de igual forma que el caso anterior.
Finalmente, compare ambas superficies.
135
INTEGRANDO UN DIQUE EN UN MDE QGIS 3
¿Cómo modificar un MDT?
En este caso desarrollamos un caso práctico utilizando herramientas disponibles en QGIS 3.
Se mostrara como alterar un Modelo Digital de Elevación (MDE) utilizando la calculadora Raster,
específicamente como incorporar un dique.
No se pretende reemplazar programas hidrológicos más adecuados, sino desarrollar una aproximación
utilizando QGIS 3.
Área de estudio:
Corresponde a una unidad de producción pecuaria (finca) donde se desarrolla cría de ganado bovino
y aves.
La finca está ubicada en una zona de piedemonte en el sur hacia colinas altas al norte, es decir, una
tendencia de zonas planas en el sur a escarpadas en el norte.
Caso de estudio:
Las actividades productivas requieren un abastecimiento continuo de agua, se evalúa la ubicación de
un pequeño dique aprovechando el relieve escarpado.
Capas y variables:
Se utilizaran las capas y el procedimiento para generar el modelo digital de terreno explicado en el
artículo “Como generar un Modelo Digital de Elevación QGIS 3”
Se suministra el sitio donde se ubicará el dique.
Procedimiento:
1.- Genere el MDE con un tamaño de pixel de 2,5 de ancho con 2,5 de alto, cargue la capa de puntos
Ubi_dique.
136
En la imagen el diamante cyan representa la ubicación del dique.
2.- Genere perfiles para decidir la dirección y longitud del dique utilizando el complemento Profile
Tool (Consulte el articulo Crear el perfil de un río en QGIS).
En la posición seleccionada la máxima longitud será de 120 metros en la cresta, la máxima elevación
cota 280 m.s.n.m.*
137
Seguidamente, creé un Shapefile de polígono y dibuje la forma del polígono, defina como nombre
DIQUE.
3.- Posteriormente convertirá el dique a formato Raster, para ello menú Raster > Conversión >
Rasterizar (vectorial a raster).
En la pestaña Capa de entrada introduzca la capa de polígono del dique (DIQUE), en A fixed value
to Burn coloque 1, en horizontal y vertical resolución coloque 2,5 (la resolución del MDE).
138
Esto genera un Raster binario de nombre Rasterizado con valor 1 que corresponden al dique
y 0 el resto.
Explicación:
▪ ( “Rasterizado@1” > 0 ) * 280 selecciona los pixeles del dique (con valor 1) luego les asigna el
valor 280 que es la cota máxima del dique.
▪ + ( “Rasterizado@1” = 0 ) * “MDE@1” el signo + concatena otra condicional. Luego selecciona
los pixeles con valor igual a 0 y les asigna el valor de elevación del MDE.
5.- Posteriormente, genere los perfiles y constate la incorporación del dique en el MDE.
139
6.- En la vista de planta y 3D se puede constatar la forma del dique.
140
Afinando el resultado
La ecuación aplicada presenta el inconveniente que asigna el valor de elevación 280 a todo el trayecto
del dique, tanto a zonas con elevaciones menores como mayores de 280.
7.- Se puede afinar el procedimiento incrementando la condicional y con el nombre de salida MDE-
dique2:
(( “Rasterizado@1” > 0 ) AND ( “MDE@1” <= 280 )) * 280 + (( “Rasterizado@1” > 0 ) AND (
“MDE@1” > 280 ))* “MDE@1” + ( “Rasterizado@1” = 0 ) * “MDE@1”
Explicación:
▪ (( “Rasterizado@1” > 0 ) AND ( “MDE@1” <= 280 )) * 280 selecciona los pixels del dique y
asigna la elevación de 280 solo aquellos en los que el MDE la elevación es menor a 280.
▪ (( “Rasterizado@1” > 0 ) AND ( “MDE@1” > 280 ))* “MDE@1” los pixels del dique en los
cuales el MDE tiene elevaciones mayores de 280, le asigna la elevación del MDE.
141
▪ ( “Rasterizado@1” = 0 ) * “MDE@1” los pixels distintos al dique le asigna la elevación del
MDE.
142
CUENCA Y ESPEJO DE AGUA DE UN EMBALSE EN
QGIS 3
Estimación del área embalsada
Los modelos digitales de terreno (MDT) son una forma de representación del relieve y al mismo
tiempo, un insumo para la extracción de múltiples parámetros.
Un MDT es una variable de entrada para procesos complejos como los hidrológicos, análisis de
visibilidad, geomorfológicos, insolación, etc.
En los artículos “como generar un modelo digital de elevación con QGIS 3” y “Como insertar un
dique en un MDE en QGIS 3”, encontrará los procedimientos iniciales así como la data para
reproducir lo planteado en la presente entrada.
Sin embargo, para fines de evaluación de variables debe respetarse la correspondencia entre la calidad
de los datos de entrada y la resolución del MDE
.
Las capas utilizadas para construir el MDE provienen de una base cartográfica 1:25.000, es decir, 1
cm equivale a 250 metros.
143
Se genera el MDE con un tamaño de pixel de 10 metros, para ello configure las dimensiones:
144
Luego, convierta el polígono a formato raster, para ello, menú Raster > Conversión > Rasterice
(vectorial a Raster).
145
2.- Delimitación automática:
En menú Procesos > Caja de Herramientas > SAGA > Terrain Analisys – Hidrology
Procedemos a corregir el MDE, seleccione Fill Sinks, en la pestaña DEM seleccione MDE_10.
Este procedimiento llenará las depresiones que pudieran evitar que se determine adecuadamente las
zonas que contribuyen.
Las demás opciones por defecto clic en Run.
146
El resultado es un raster binario, el área que drena al dique tendrá valor de 100, compárela con la que
realizo manualmente.
Espejo de agua:
Posteriormente se determina el espejo de agua considerando la máxima capacidad de llenado.
En el artículo “Como insertar un dique en un MDE en QGIS 3” diseñamos el dique y le asignamos
como altura máxima 280 m.s.n.m. (aliviadero).
Para determinar el espejo de agua sencillamente extraemos del MDE corregido, dentro del área de
drenaje, las zonas con elevaciones menores a 280 m.s.n.m.
147
En la imagen el resultado:
148
Genera un archivo html donde podemos ver la superficie ocupada por el espejo de agua (56.100 m2
o 5,61 ha), específicamente los pixeles con valor 1.
149
Posteriormente podemos realizar representaciones en tres dimensiones:
150
CREAR UNA COLECCIÓN DE MAPAS PDF EN QGIS
3
La nueva versión de QGIS ha incorporado mejoras en las salidas para impresión incluyendo la opción
para crear reportes con nuestros mapas y tablas de atributos.
La posibilidad de automatizar la preparación de múltiples mapas para impresión ya existía en
versiones anteriores de QGIS con la herramienta Atlas, sin embargo, en la versión QGIS3 el
módulo Report amplía las opciones mientras nos facilita el trabajo. En ArcGIS se llama Data Driven
Pages.
Caso práctico: preparar para impresión mapas agrupados en un documento PDF con su
respectiva portada en no más de 5 minutos:
Contenido del documento: queremos preparar un documento en el que se muestre un mapa de
ubicación del estado Cojedes, Venezuela; seguidamente los mapas de cada una de sus parroquias.
1.- Prepara una vista, en ella carga las diferentes capas Shapefile, como son: poligonal del estado,
hidrografía, vialidad, municipios y parroquias.
2- Configura el etiquetado de las capas municipios y parroquias por sus respectivos nombres.
3.- En el menú Proyecto selecciona New Report aparece una caja de dialogo que me permite darle un
nombre al reporte que voy a crear, seguidamente se abre la ventana de Report que es similar a la de
Layouts, con la diferencia del panel derecho llamado Report Organizer.
151
4.- Estructura el documento en el Report Organizer y llena de contenido las páginas con las
herramientas de Layout. Selecciona la opción Include Report Header y luego le da clic a editar, esto
crea la primera página del documento.
152
5.- Lista la Portada, adiciona paginas! inicia con una página estática, para ello clic al botón con el
signo de suma y luego a la opción Static Layout Section, crea una nueva página.
Igual que en un Layout (composición de mapas) inserta el mapa de la vista, al terminar asegurarse de
seleccionar en las propiedades del mapa Bloquear Capas.
153
6.- Ahora crea otra página para mostrar un mapa con la división político territorial por parroquias,
para ello, prepara un mapa en la vista. Luego entra al asistente de Report y selecciona
en Report, repetir el procedimiento anterior:
7.- En la siguiente sección genera automáticamente los mapas de las 13 parroquias Primero,
asegurarse bloquear las capas del mapa creado anteriormente. Hacer nuevamente clic en el botón con
el signo de suma pero esta vez escoge la opción Field Group Section. Aparece un nuevo elemento en
el Panel Report llamado Grupo y una serie de botones. Antes de continuar prepara el mapa con las
parroquias
154
8.- La configuración del conjunto de mapas es similar a configurar la herramienta Atlas, en la opción
Capa escogemos Parroquia y en Campo Parroquia (el campo que contiene el nombre de cada
parroquia). Luego selecciona Include Body y esto crea la página patrón que se repetirá para cada
Parroquia, da clic en Editar y configura el Layout:
9.- Aquí configura un mapa similar a un Layout, la escala, flecha norte, rejilla de coordenadas. Lo
IMPORTANTE es que configurar la opción Controlled by Report, selecciona la casilla Margen
alrededor de Objeto
155
10.- Finalmente en la barra de herramientas y exporta el reporte como PDF, genera un documento
con tantos mapas como entidades (Parroquias) tenga la capa.
156
Al abrir el archivo PDF creado se puede ver que contiene un total de 18 páginas, 17 de ellas son
mapas:
Con un poco de dedicación se puede generar vistosos documentos, insertando textos, marcos,
imágenes, tablas, incluso falta probar si se puede añadir una vista 3D.
157
CONVERTIR PUNTOS, LÍNEAS Y POLÍGONOS KML
A SHP EN QGIS
Si has realizado tu trabajo en Google Earth, con tan solo unos sencillos pasos se puede pasar los
ficheros KML a shapefiles. No se requiere instalar plugins adicionales, sino simplemente trabajar con
la instalación existente de QGIS.
Primero se requiere tener los puntos, líneas o polígonos en Google Earth, y guardarlo en formato kml
o kmz.
▪ Hacer clic en Añadir capa vectorial del menú Capa > añadir capa, en la ventana emergente
hacer clic en Explorar.
▪ Explorar la ubicación del archivo, en caso de no mostrarse en la ventana, en el tipo de formato
seleccionar kml, kmz y abrir el archivo.
▪ Una vez cargado el polígono hacer clic derecho sobre la capa y seleccionar la opción Guardar
como.
▪ Posteriormente en la opción Formato se selecciona Archivo shape de ESRI. En la actual
ventana asignar un nombre y directorio para el archivo, en SRC se coloca el sistema de
proyección.
▪ Automáticamente se carga el shapefile como una nueva capa al proyecto actual de QGIS.
▪ Este mismo procedimiento se puede realizar la puntos, líneas y polígonos (kml, kmz) creadas
en Google Earth.
▪ Los nuevos shapefiles generados por este proceso se los puede abrir con cualquier aplicación
que soporte este formato.
158
159
EXCEL Y QGIS: IMPORTAR Y
EXPORTAR DATOS.
¿Cómo importar y exportar datos de excel a
QGIS?
En este artículo vamos a ver la manera tanto de importar datos de excel a
QGIS, como las opciones con las que contamos para exportarlos a formato de
datos.
Antes que nada, supongamos que tenemos una tabla sencilla, de puntos, por
ejemplo una hoja de excel con las coordenadas y los códigos de algunos
aeropuertos.
Importar Excel a QGIS.
Lo primero que habrá que hacer es pasar nuestros datos de excel a un formato
de texto plano. CSV es el más recomendado para esto, por tanto, desde excel
exportaremos los datos a CSV. A mí me gusta separarlos por “;” (semicolon), ya
que por coma puede darse inconsistencia en algunos datos.
Exportar a CSV
Tras esto, la herramienta de importación de CSV de QGIS nos permitirá
seleccionar el archivo, y desde esta misma ventana podremos decirle el nombre
de las columnas con las coordenadas, la separación, u otras opciones.
160
Importar CSV
Solo tendremos que darle a la capa una proyección, y tendremos disponible
estos datos para realizar análisis espaciales.
Aeropuertos
¿Y cómo los exporto?
Después de haber recibido muchas consultas al respecto, he visto que hay
muchos artículos sobre cómo importar archivos de excel, pero pocos sobre cómo
exportarlos.
Podríamos encontrarnos con el caso de necesitar tener los datos de una capa en
una tabla de excel, por lo que QGIS permite también de manera bastante sencilla
esta acción.
Antes de nada, debemos ir a nuestra capa en la tabla de contenidos y con el
botón derecho seleccionar guardar como.
161
Guardar cómo
Opciones
¿Y si no quiero exportar todas las columnas?
Desde la versión 2.16, podemos realizar esta acción de una manera sencilla sin
recurrir a ningún plugin como antes. Ahora, bajo la pestaña de Seleccionar
campos para exportar, podemos elegirlos sencillamente marcando el checkbox:
162
Cómo convertir un dxf en shapefile en GRASS-QGIS
El formato de archivos dxf, acrónimo del inglés Drawing Exchange Format, se usa para
dibujos de diseño asistido por computadora pensando en la interoperabilidad entre los
archivos *.dwg usados por el programa AutoCAD y el resto de programas del mercado. Si
se pretende abrir un archivo dxf en QGIS ello no es posible. Para probar pueden bajar una
muestra de un dxf de aquí:
http://justcad.com/downloads.html
En la imagen anterior se observan dos capas, una de líneas y otra de puntos, que se guardan
de manera independiente como shapefiles (formato de ESRI) usando el menú contextual de
la ventana de Capas (botón derecho del ratón sobre la capa). Se
nombraron frame_dxf.shp y points_dxf.shp; respectivamente.
163
en ambiente de QGIS (no en GRASS-QGIS).
164
IMPORTAR COORDENADAS DE
EXCEL A QGIS Y CREAR
POLÍGONOS
Una de las rutinas más comunes en el uso de Sistemas de Información Geográfica es la
construcción de capas espaciales a partir de la información proveniente de campo. Sea
que esto represente coordenadas, vértices de parcelas o retícula de elevaciones, la
información suele venir en ficheros separados por coma u hojas de cálculo de Excel.
1. El archivo de coordenadas geográficas en Excel.
2.
En este caso, estoy intentando importar los asentamientos humanos de la República de
Cuba, el que he bajado de diva-GIS, que por cierto es uno de los mejores sitios para
descargar datos geográficos de cualquier país. Como se puede ver, las columnas B y C
contienen la información relativa a la latitud y longitud en forma de coordenadas
geográficas.
165
Vectorial > XY tools > OpenExcele file as attribute table or
Point layer
En caso que el fichero esté guardado con extensión .xlsx el navegador no lo mostrará, pues
filtra solo los archivos con extensión .xls. Non es problema, podemos aplicar técnicas
antiguas de DOS y escribir en el cambio de nombre, el filtro: *.* (asterisco punto asterisco)
y hacemos Enter; esto va a permitir que se vean todos los archivos en esa
localización. Igual pudimos haber escrito *.xls y hubiera filtrado solo los archivos de
extensión .xls.
166
Luego nos aparece un panel en el que debemos indicar cuál columna equivale a la coordenada en
X, en este caso elegimos la columna de longitud, la columna de latitud para la coordenada Y.
Y allí lo tenemos. La consulta muestra que se ha guardado la capa con los datos que contenía el
fichero de asentamientos humanos de Cuba, que incluye el nombre, la latitud, la longitud, la
clasificación y la provincia administrativa.
167
3. Crear polígonos a partir de coordenadas
En caso que, queremos no solo importar los vértices sino que también se cree un polígono
en el orden de esas coordenadas, podemos utilizar el plugin Points2One. Este plugin
permite identificar como se llamará la capa de destino, si lo que importaremos se construirá
como líneas o como polígono.
168
4. Cómo importar coordenadas de Excel a otros
programas CAD/GIS.
Como recordarán, este proceso lo hemos hecho con muchos otros programas. Tan simple como
QGIS, pocos. Pero aquí aparece cómo hacerlo con AutoCAD, Microstation, Manifold
GIS, AutoCAD Civil 3D, Google Earth.
169
QGIS 3.6 NOOSA:
NOVEDADES, DESCARGA E
INSTALACIÓN
Para todos aquellos que aún no conozcan QGIS y que comienzan a interesarse en este
potente software comentar que se trata de un Sistema de Información Geográfica de
Software Libre, siendo un proyecto de la Open Source Geospatial Foundation (OSGeo).
Desde el inicio del proyecto en el año 2002 hasta la actualidad, QGIS ha experimentado
una evolución que le ha llevado a ser un programa GIS que cubre la mayor parte de las
necesidades en visualización de datos.
QGIS se encuentra bajo la General Public License (GPL), lo cual permite al usuario
modificar su código fuente y garantizar la existencia y acceso a un programa GIS gratuito.
Detrás del proyecto hay una comunidad, formada por un comité directivo, un equipo de
desarrolladores y los propios usuarios del software. Las contribuciones de esta comunidad
pueden ser en forma de código, informes o correcciones de errores, documentación,
promoción y apoyo. MappingGIS apoya el proyecto como patrocinadores desde noviembre
de 2014.
Descubre por qué QGIS es el SIG líder de Código Abierto para escritorio.
Índice [cerrar]
• QGIS 3.6 Noosa
• Descarga de QGIS 3.6 Noosa
• ¿Por qué QGIS 3.6 en esta ocasión se llama Noosa?
• Novedades de QGIS 3.6 Noosa
o Expresiones
o Herramientas de mapa
o Interfaz de usuario
o Simbología
o Digitalización
o Gestión de datos
o Formularios y widgets
o Diseñador de mapas
o Herramientas de análisis
o Procesamiento
170
o Opciones de aplicación y proyecto
o Navegador
o General
o Proveedores de datos
o QGIS Server
o Programabilidad
o Características 3D
171
QGIS 3.6 NOOSA
La versión 3.6 es una versión LR, es decir, inicia una nueva fase de desarrollo y, a su vez,
la nueva versión LTR (Long-Term_Release), pasa a ser la 3.4.5
Según la hora de ruta, el 21 de junio de 2019 se lanzará la próxima versión LR, la 3.8
172
DESCARGA DE QGIS 3.6
NOOSA
Puedes descargar QGIS 3.6 Noosa desde la página oficial de descargas del proyecto,
seleccionando el instalador acorde a la versión de tu sistema operativo.
173
¿POR QUÉ QGIS 3.6 EN ESTA
OCASIÓN SE LLAMA NOOSA?
QGIS 3.6 se denomina Noosa, nombre decidido en la 20ª reunión de desarrolladores entre
el 22 y 26 de agosto de 2018 en Zanzíbar.
174
Esta es la duodécima versión que lleva un nombre de una localidad terrestre, desde las
versiones 0.8.1 a 1.5 fueron denominadas con nombres de las lunas de Júpiter y Saturno,
por ejemplo, Titán, Io, Tethys, etc. Desde la versión 1.8 se adoptaron nombres terrenales:
Wrocław, Lisboa, 2.0 Dufour, 2.2 Valmiera , 2.4 Chugiak, 2.6 Brighton, 2.8 Wien, 2.10 Pisa,
2.12 Lyon, 2.14 Essen, 2.16 Nødebo, 2.18 Las Palmas, 3.0 Girona, 3.2 Bonn y
3.4 Madeira.
175
NOVEDADES DE QGIS 3.6
NOOSA
Esta nueva versión de QGIS incorpora una serie de mejoras y novedades que continúan
completando el nuevo proyecto QGIS 3.x.
Expresiones
Herramientas de mapa
3.Etiqueta de título en mapa. Se ha añadido una nueva opción que permite a los usuarios
agregar de forma rápido un título a las composiciones de mapa, que se presenta con un
fondo con transparencia y color personalizables.
176
8. Modo capa actual como predeterminado. En la herramienta de vértices el modo “todas
las capas” solía ser el predeterminado, en la actualidad se ha modificado por el modo “capa
actual”.
9. Herramienta de identificación de objetos espaciales compatible con capas de
malla.
Interfaz de usuario
10. Nuevo tema para la interfaz de usuario. Se ha añadido el tema Blend of Gray para la
interfaz de usuario (mezcla de grises).
Simbología
13. Manejo mejorado de las bibliotecas de estilos .xml. El panel del navegador muestra
ahora bibliotecas de estilos xml de QGIS, que pueden utilizarse con drag and drop.
14. Nueva regla Forzar orientación a mano derecha en la representación de
polígonos. Esta nueva opción está disponible en el botón Avanzado para símbolos de
relleno. La corrección de orientación se aplica sólo al renderizar, y la geometría original de
la entidad no se modifica.
177
15. Nueva opción para renderizar exclusivamente anillos exteriores o interiores en
líneas simples o de marcador.
16. Nuevo marcador de imagen ráster. Se ha añadido un nuevo tipo de simbología
a QGIS que permite al usuario utilizar un archivo ráster para mostrarlo como marcador. La
simbología de este modo permite la personalización de ancho, alto, relación de imagen,
rotación y opacidad.
17. El relleno ráster puede tener imágenes configuradas con acceso a url’s remotas o
imágenes incrustadas.
178
18. Posibilidad de utilizar valores de representación para simbología de tipo clasificado.
19. Fusión de categorías en el modo de simbología categorizada. Los usuarios podrán
seleccionar varias categorías y agruparlas en una única lo que permite un estilo más simple de
representación cuando se presentan un elevado número de categorías.
20. Opción para mostrar los colores del proyecto en los botones de datos definidos
por el usuario. El menú Colores contiene todos los colores definidos como parte del
esquema de color el proyecto actual (que se define a través de las propiedades del
proyecto).
179
• Nuevo modo “Vista de lista”
• Nueva pestaña “Todos” que muestra todos los símbolos coincidentes.
Digitalización
Gestión de datos
Formularios y widgets
25. Nuevo widget en formularios para campos binarios tipo blob. Ofrece una etiqueta
que muestra si el campo tipo blob está vacío o no.
Diseñador de mapas
180
30. Los elementos de diseño pueden bloquear etiquetas de mapa. Esta función permite
que se marquen como “bloqueadores” para etiquetas de mapa algunos elementos de la
composición, lo que evita que las etiquetas se coloquen por debajo de los mismos.
31. Advertencias sobre las exportaciones. Al exportar diseños de impresión, QGIS ahora
comprueba automáticamente el diseño en busca de errores comunes y advierte a los
usuarios de cuando se encuentran:
32. Control sobre la posición de apilamiento de las extensiones de vista general del
mapa.
33. Nuevas variables de expresiones para elementos de leyenda. Además, los símbolos
representados como parte de una leyenda ahora usan el contexto de expresión completo
disponible para la leyenda.
34. Expresiones dentro del texto del elemento de leyenda. Esta función permite que las
expresiones se incrusten directamente dentro del texto del elemento de leyenda.
Herramientas de análisis
Procesamiento
181
37. Explotar algoritmo HStore. Este algoritmo crea una copia de una capa de entrada y
agrega un nuevo campo para cada clave única que se encuentra en un campo de
tipo HStore.
38. Forzar regla “a mano derecha”. Este nuevo algoritmo obliga a las geometrías de tipo
poligonal a respetar la regla “a mano derecha”, en la que el área que está delimitada por un
polígono está a la derecha del límite.
182
39. Nuevo algoritmo para extraer campos binarios.
40. Nuevo algoritmo para cortar líneas por longitud máxima. Este algoritmo toma una
capa de entrada y divide cada entidad en varias partes, de modo que ninguna parte es más
larga que la longitud máxima especificada.
183
41. Nuevo tipo de parámetro para configuración de autenticación. El
parámetro QgsProcessingParameterAuthConfig permite la selección de las
configuraciones de autenticación disponibles (y la creación de otras nuevas).
42. Nuevas opciones de remuestreo y formato.
43. Algoritmo ráster de estadísticas zonales. Este algoritmo calcula estadísticas para los
valores de una capa ráster, categorizados por zonas definidas en otra capa ráster.
184
46. Nuevo algoritmo para calcular superficies volumétricas. Se ha portado del bloque
de geoprocesamiento de SAGA como algoritmo nativo de QGIS.
47. Exportar modelo como script de Python. Se ha incluido esta opción que ya estaba
disponible en la versión 2.x de QGIS.
48. Borrar opciones de configuración desde la pestaña Avanzado. Este cambio agrega
un nuevo menú que se muestra en la pestaña Avanzado en el cuadro de diálogo de
configuración.
185
Navegador
49. Nueva opción Exportar a archivo para capas ráster y vectoriales dentro del panel
del navegador.
186
50. Propiedades nativas de archivos y carpetas. Nueva opción para acceder a las
propiedades de un archivo o directorio con clic derecho de ratón en el panel del navegador.
51. Vista previa de capas y atributos. QGIS 3.6 recupera la funcionalidad para
previsualizar capas y sus atributos, disponible accediendo a las propiedades de la capa
desde el navegador.
187
54. Abrir directorio en path. De este modo, podemos abrir una nueva ventana con el
terminal cmd en ese directorio.
General
Proveedores de datos
58. Información de servicio abierta para las capas de ArcGIS Feature Server. Las
capas de ArcGIS Feature Server cuentan con una nueva opción de menú que permite abrir
directamente la página web de información de servicio para la capa seleccionada.
59. Soporte binario blob para proveedores de OGR. En lugar de convertir los campos
binarios en cadenas truncadas, estos valores son recuperables como sus contenidos
binarios originales.
60. Marcador de imágenes ArcGIS Feature Server y soporte de relleno de imágenes.
188
61. Soporte de etiquetado de ArcGIS Feature Server.
62. Soporte booleano y campo binario para capas de memoria.
63. Soporte JSON para GeoPackage.
64. Opción de guardar selección de clave principal en el proveedor PostgreSQL.
65. Soporte de valores únicos de varios campos de ArcGIS Feature Server.
QGIS Server
Programabilidad
68. Nueva clase QgsExifTools. Contiene utilidades para recuperar la geoetiqueta de las
imágenes y para configurar la geoetiqueta de una imagen.
69. Iterar sobre partes de un QgsGeometry. Esta nueva API de PyQGIS permite una
iteración sencilla de todas las partes de una geometría, independientemente de su tipo.
70. Mejora en la API QgsLineString de PyQGIS.
71. Mejoras en la API de PyQGIS Geometry Collection.
72. Creación de algoritmos de procesamiento de forma más sencilla a través de @alg
decorator.
73. Bloqueo de solicitudes de red.
74. Controles de validación personalizados en las exportaciones de diseño.
75. API para el cálculo de líneas geodésicas.
Características 3D
76. Sombreado del terreno. Nuevas opciones para que el usuario elija cómo debe de
renderizarse el terreno:
189
78. Procesador 3D basado en reglas.
190
CREANDO MARCADORES
ANIMADOS CON LEAFLET
(SVG, GIFS Y FONT
AWESOME)
Leaflet es una excelente librería JavaScript empleada para publicar mapas elegantes, con
un “aire” moderno de forma rápida y eficaz. Pero además es capaz de integrar recursos
para hacer los mapas más vistosos buscando llamar la atención del usuario. Por ejemplo
con Leaflet podemos animar los marcadores.
En este blog ya tratamos anteriormente cómo insertar iconos de bootstrap en Leaflet, una
de cuyas opciones es la de hacer girar al marcador. Ahora nos centraremos en crear
animaciones utilizando imágenes SVG, GIF y con Font Awesome.
Utilizar un formato vectorial tiene además otras ventajas, como permitirnos crear efectos y
animaciones. Las imágenes SVG se pueden editar con programas open
sourcecomo INKSCAPE.
191
Disponemos también de muchas páginas web que nos proporcionan ejemplos y
explicaciones sobre el uso de SVG. Por ejemplo un rectángulo se forma de la siguiente
manera:
Como vemos se compone de unas etiquetas html <svg> en cuyo interior definimos el
rectángulo (rect) mediante sus dimensiones: largo y ancho y un estilo.
Pasos a seguir
El rectángulo que hemos definido antes podemos utilizarlo como un marcador en Leaflet.
Lo primero es crear un nuevo icono mediante una instancia a L.Icon
Como vemos, hemos introducido algunos cambios en la definición del rectángulo que son
los que nos proporcionan la animación. Estamos modificando la transparencia del marcador
en un periodo de 5 segundos. Esto hará que el marcador vaya “difuminándose” hasta
192
desaparecer al cabo de 5 segundos, y a continuación se vuelve a hacer visible, en un bucle
continuo.
Una vez completada la generación del icono SVG podemos crear el marcador e insertarlo en
el mapa.
<!DOCTYPE html>
<head>
<meta charset="utf-8">
<script src="https://cdnjs.cloudflare.com/ajax/libs/d3/3.5.5/d3.min.js"></script>
<link rel="stylesheet" href="https://unpkg.com/leaflet@1.4.0/dist/leaflet.css" />
<script src="https://unpkg.com/leaflet@1.4.0/dist/leaflet.js"></script>
<style>
body { margin:0;position:fixed;top:0;right:0;bottom:0;left:0; }
#map {
height: 100%;
}
193
.leaflet-marker-icon{
width: 200px;
height: 200px;
}
</style>
</head>
<body>
<div id="map"></div>
<script>
var map = L.map('map').setView([40.965, -5.664], 14);
L.tileLayer('https://{s}.tile.osm.org/{z}/{x}/{y}.png', {
attribution: '© <a href="https://osm.org/copyright">OpenStreetMap</a> contributors'
}).addTo(map);
var CustomIcon = L.Icon.extend({
options: {
iconSize: [120, 90],
iconAnchor: [22, 94],
popupAnchor: [-3, -76]
}
});
var svgrect = "<svg xmlns='https://www.w3.org/2000/svg'><rect x='0' y='0' width='50' height='50' fill='red'><animate
attributeType='CSS' attributeName='opacity' from='1' to='0' dur='5s' repeatCount='indefinite' /></rect></svg>";
Utilizar GIF en Leaflet no tiene nada de particular, podemos utilizarlo igual que cualquier otro
formato de imagen. Un ejemplo es el siguiente:
<!DOCTYPE html>
<html>
<head><title>GeoJson y Leaflet</title>
<meta charset="UTF-8">
<link rel="stylesheet" href="https://unpkg.com/leaflet@1.4.0/dist/leaflet.css" />
194
<script src="https://unpkg.com/leaflet@1.4.0/dist/leaflet.js"></script>
<style>
#map {
width: 700px;
height: 600px; }
#div1 {
font-size:48px;
}
</style>
</head>
<body>
<div id="map"></div>
<script>
var map = L.map('map',{center: [36.8021480 , -5.1393441],zoom: 7});
L.tileLayer('https://{s}.tile.osm.org/{z}/{x}/{y}.png').addTo(map);
var geojson = [{
"type": "Feature",
"geometry": {
"type": "Point",
"coordinates": [-5.9731700, 37.3828300]
},
"properties": {
"name": "Sevilla",
"title": "Andalucía"
}
},{
"type": "Feature",
"geometry": {
"type": "Point",
"coordinates": [ -4.4203400 ,36.7201600]
},
"properties": {
"name": "Málaga",
"title": "Andalucia"
}
}];
El código anterior utiliza un geoJSON para crear y asignar propiedades a dos marcadores.
La creación del icono, como vemos, es igual que en el ejemplo que pusimos para SVG, es
decir: realizamos una instancia a L.Icon y luego definimos la ruta (URL) de la imagen. En
esta ocasión se trata de una bandera de Andalucía flameada por el viento.
195
FONT AWESOME – EL CONJUNTO DE
ICONOS Y HERRAMIENTAS MÁS
POPULARES DE LA WEB
De entre la gran cantidad de recursos web que nos proporcionan iconos, uno de los más
importantes es sin duda Font Awesome.
Podemos utilizar un marcador para crear una función con JavaScript que defina una
animación, de la siguiente forma:
function chargebattery() {
var a;
a = document.getElementById("div1");
a.innerHTML = "";
setTimeout(function () {
a.innerHTML = "";
}, 200);
setTimeout(function () {
196
a.innerHTML = "";
}, 400);
setTimeout(function () {
a.innerHTML = "";
}, 600);
setTimeout(function () {
a.innerHTML = "";
}, 800);
};
Mediante setTimeout hacemos que se ejecute una función en un tiempo determinado medido
en milisegundos. De esta forma podemos modificar un div en el tiempo de tal forma que se
crea una animación. En el ejemplo que estamos utilizando sería la carga de una batería. Como
paso previo hemos definido un div utilizando la clase de leaflet L.divIcon.
setInterval(chargebattery, 1000);
Y creamos el marcador:
197
R Y GIS: QUÉ ES R Y SU
RELACIÓN CON LOS SIG
R es un entorno de software libre y lenguaje de programación para realizar análisis
estadístico. Se trata de uno de los lenguajes de programación más utilizados en
investigación por la comunidad estadística.
El proyecto data de 1997 y se encuentra bajo la Licencia Pública General (GPL por sus
siglas en inglés) del proyecto colaborativo de software libre GNU.
R Y RSTUDIO
La interfaz del núcleo R es una ventana de línea de comandos que proporciona una
excelente flexibilidad y control, pero no es fácil trabajar de este modo.
198
R Y GIS
La relación entre R y los Sistemas de Información Geográfica es larga. R se ha utilizado
para el análisis de datos espaciales desde los primeros días de R. Iniciado por Roger
Bivand, los esfuerzos estructurados de la comunidad comenzaron en 2003 en un taller
previo a la conferencia durante el DSC2003.
199
PAQUETES DE R PARA TRABAJAR CON
DATOS ESPACIALES
Gracias al código abierto se han desarrollado paquetes modulares (actualmente + de
13.500), que son complementos para temas específicos por la comunidad de
desarrolladores.
• Todas las funciones y métodos del paquete sf que emplean datos espaciales tienen
el prefijo st_, que se refiere a espacial y temporal (al igual que en PostGIS).
• Los objetos geográficos son data.frames o tibbles con una columna de geometrías.
• Representa de forma nativa en R los 17 tipos de objetos geográficos simples para
todas las dimensiones (XY, XYZ, XYM, XYZM)
• Interfaces a GEOS para soportar el modelo topológico DE9-IM.
• Se conecta a GDAL, es compatible con todas las opciones de controlador, columnas
de fecha (Date) y fecha y hora (DateTime) (POSIXct) y transformaciones del sistema
de coordenadas de referencia a través de PROJ.
• Utiliza serializaciones WKB (well-know-binary conocidas escritas en C++/Rcpp para
fast I/O con GDAL y GEOS.
• Lee y escribe directamente en bases de datos espaciales como PostGIS usando DBI.
200
(Illustration (c) by Allison Horst)
El paquete sp
Clases y métodos para datos espaciales; el documento de clases donde reside la
información espacial, para datos 2D o 3D. S. Las funciones que se proporcionan son:
imprimir datos como mapas, selección espacial, métodos para recuperar coordenadas,
subconjuntos, impresión, resumen, etc.
Sin embargo el paquete sf es más fácil de aprender y reemplazará paulatinamente a sp, por
lo que se recomienda utilizar el paquete sf. Ya hay páginas web para realizar la migración.
El paquete raster
El paquete raster permite manejar datos ráster con el objeto de clase RasterLayer. El raster
se obtiene a partir de una matriz que contiene los valores de cada celda.
201
El paquete stars
Muchos de los datos raster disponibles en la actualidad son dinámicos y vienen como
series temporales de ráster para diferentes temas.
Los datos raster recientes, como las imágenes de satélite, modelos climáticos o los datos
de pronóstico del tiempo, que no caben en el almacenamiento local necesitan otro tipo de
tratamiento. Existe un paquete para el análisis raster, denominado stars que:
library(tidyverse)
#> ── Attaching packages ────────────────────────────────── tidyverse 1.2.1 ──
202
Fuente: https://keen-swartz-3146c4.netlify.com/intro.html#a-first-map
203
CONFIGURAR R EN QGIS 3 MEDIANTE EL
PLUGIN PROCESSING R PROVIDER
En QGIS 3 existe un complemento denominado Processing R Provider, creado por North
Road.
Una vez instalado podremos configurarlo. Para ello debemos ir al menú superior
Configuración > Opciones…
204
A partir de ese momento si abrimos la caja de herramientas de Procesos (menú superior
Procesos > Caja de herramientas) veremos los algoritmos de R.
LEAFLET Y R
Leaflet es una de las bibliotecas de JavaScript de código abierto más populares para mapas
interactivos. Lo utilizan sitios web que van desde The New York Times hasta GitHub, así
como especialistas en SIG como OpenStreetMap, Mapbox o CARTO.
205
Lo dejamos aquí, aunque podríamos hacer un articulo mucho más largo y será mejor
continuar ampliando el tema con artículos más específicos. Como veis las posibilidades de
trabajo con R y GIS son abundantes, ya que ese trata de una librería bastante madura y
con mucho futuro.
206
CÓMO RECLASIFICAR UNA
IMAGEN EN QGIS
Debido a la importancia y excelentes resultados que el software GRASS proporciona en el
tratamiento de la información, QGIS también ha incorporado gran parte de los algoritmos
del mismo para poder explotar datos cargados directamente a la interfaz de QGIS.
Conviene señalar que no todos los algoritmos disponibles con el proveedor GRASS lo están
en la Caja de herramientas de procesos (ya que en algunos casos se requiere que se
encuentren dentro de la estructura de datos nativa de GRASS GIS).
Si quieres saber más acerca de GRASS y sus utilidades, puedes visitar nuestro
artículo Qué podemos hacer con GRASS GIS 7.
207
QUÉ ES LA RECLASIFICACIÓN RÁSTER
La reclasificación (o también
llamada reasignación) ráster supone una
modificación de los valores de una imagen de forma total o parcial, agrupando
los valores que presenta la imagen de entrada en una clasificación de rangos definida por
el usuario.
Este tipo de operaciones son útiles sobre todo cuando queremos establecer una
clasificación, por ejemplo, de rangos de pendientes para conocer la orografía de nuestra
zona de estudio, modelos de orientaciones, altitud, índices de vegetación, etc…
Es decir, podemos reclasificar imágenes que representan tanto variables continuas como
pueden ser la elevación o la precipitación, por ejemplo, como variables discretas, por
ejemplo, los índices de vegetación.
208
Por defecto, QGIS carga la Caja de herramientas en modo simplificado, mostrando
únicamente los algoritmos propios de QGIS. El modo avanzado contiene características y
algoritmos adicionales que requieren la comprensión de sus propias aplicaciones, por lo
que son un tema avanzado.
209
Podríamos realizar una clasificación de la pendiente en los siguientes rangos:
input_categories=output_category [label]
Una vez creado el archivo, lo guardamos en nuestro pc como Documento de texto (*.txt)
y codificación ANSI.
EL MÓDULO R.RECLASS
Una vez contamos con la imagen ráster que deseamos reclasificar y con el documento que
alberga las reglas de reclasificación, podemos continuar el proceso. Es el momento de
configurar la herramienta, podemos utilizar el buscador de la Caja de herramientas de
procesado para seleccionar el algoritmo r.reclass
210
Introduciríamos los siguientes parámetros en la herramienta:
211
Puedes emplear la herramienta Identificar objetos espaciales para conocer la
clasificación de cada uno de los píxeles.
212
CÓMO CONVERTIR UN PDF EN
SHP, GEOPACKAGE O
CUALQUIER OTRO FORMATO
VECTORIAL
En esta entrada vamos a explicar el proceso para convertir un archivo PDF en SHP o en
cualquier otro formato vectorial (GeoPackage, GeoJSON, GML, etc).
Hemos creado un mapa muy sencillo (contiene únicamente una capa de ríos) y lo hemos
exportado a PDF. Este es su aspecto:
213
Comenzamos:
• Opción 1. Utilizar un conversor online como Convert PDF to Autocad Free Online, el
proceso consta de dos pasos:
• Seleccionar el PDF a convertir.
• Introducir el correo electrónico para recibir el enlace de descarga del archivo.
El tiempo de espera hasta que recibamos el email es de 30 minutos.
Aidecad tiene una versión de prueba que permite realizar 20 conversiones gratuitas.
• Paso 1. Clic en el botón “Add PDF Files “Add Folder” para cargar los archivos PDF
de origen.
• Paso 2. Elegimos el formato de salida y seleccionamos la carpeta de salida para
guardar los archivos resultantes.
• Paso 3. Clic en el botón “Convert” para completar la conversión.
214
Hemos probado las dos opciones mencionadas y ambas funciona bien, el DXF generado
es satisfactorio.
El primer paso es añadir los archivos DXF a la interfaz. Añadimos el DXF usando el botón
“Añadir capa vectorial”:
215
A continuación podemos guardar las entidades en otro formato haciendo clic en el menú
superior Capa… Capa > Guardar como…
216
Pulsamos Aceptar y la nueva capa se añade al proyecto.
Una vez añadida la capa en formato SHP, editando la capa podemos limpiarla, eliminando
las líneas que no nos interesen.
Vector Bender es un complemento para QGIS que permite distorsionar una capa
vectorial para adaptarla, de la misma manera que el georreferenciador funciona para
los ráster.
217
GEOWE: PLATAFORMA GIS LIBRE
PARA CREAR Y EDITAR TUS
PROPIOS MAPAS EN LA WEB
La iniciativa GeoWE es un proyecto software OpenSource, que nació en el año 2015,
orientado a la edición de Información Geográfica en Web. Surgió por la motivación de un
equipo de desarrolladores con ganas de crear las herramientas necesarias para
la visualización y edición de datos espaciales, y ofrecer al mundo SIG una nueva alternativa
para el desarrollo de aplicaciones web orientadas a la edición de datos geográficos.
GeoWE es una de las 10 aplicaciones más importantes para publicar mapas en la nube. En
la actualidad ofrece una amplia gama de servicios y/o aplicaciones:
¿QUÉ ES GEOWE?
GeoWE es un Sistema de Información Geográfica en la Web que tiene como objetivo
principal el llevar la potencialidad de la edición de datos geográficos de los SIG de
Escritorio a la Web. Es un proyecto software orientado a la creación y edición de
Información Geográfica, cuya culminación toma la forma de una aplicación cliente accesible
desde diversos dispositivos.
218
Objetivo
El objetivo principal es constituir una plataforma SIG Web de uso libre para todos los
usuarios de la comunidad y de código abierto. La plataforma está orientada a ofrecer las
funciones y herramientas necesarias para realizar edición avanzada de geodatos, de la
manera más directa y sencilla para el usuario final.
Tecnología
• Framework Google Web ToolKit: para maximizar la compatibilidad con los distintos
navegadores.
• Librería de componentes reutilizables GXT basada en GWT: para obtener una
interfaz de usuario más rica y completa.
• Librería OpenLayers: para la gestión y representación de la información geoespacial.
• Framework Errai basado en GWT: para la construcción de aplicaciones web
complejas.
• Biblioteca JTS (Java Topology Suite): para operaciones espaciales y de
geoprocesamiento.
219
• La interoperabilidad permite integrar datos abiertos sin problemas.
220
ALGUNOS CASOS DE USO
Tal y cómo indican en su web, con el SIG WEB de GeoWE se pueden resolver necesidades
reales de diferentes sectores y temáticas:
• Urbanismo.
• Deporte y ocio.
• Agricultura.
• Emergencia y Cooperación.
• Logística y Gestión de activos.
• Etc.
CONCLUSIONES Y AGRADECIMIENTO
Cómo siempre la mejor forma de conocer una aplicación es probarla y trastear con ella.
Nosotros lo hemos hecho y estamos muy satisfechos con los resultados obtenidos. Es una
herramienta muy intuitiva y fácil de manejar que ofrece algunas de las funcionalidades más
usadas y demandadas en el mundo GIS y estamos seguros que en próximas
actualizaciones nos ofrecerá nuevas mejoras que harán mucho más productivo nuestro
trabajo.
Por nuestra parte, queremos agradecer al equipo de GeoWE el gran trabajo que están
realizando y gracias por brindarnos la oportunidad de poder utilizar la aplicación de forma
libre y abierta. Sin duda alguna, el proyecto GeoWE constituye un ejemplo de que cualquier
objetivo se puede conseguir con iniciativa, motivación y perseverancia.
221
222