Conseguir capas de mapas útiles (Actualizado)

Existen infinidad de lugares donde conseguir capas de mapas que puedan ser útiles. En este post pretendo hacer una selección de los lugares más útiles para obtener capas útiles para la Ecología, el medio ambiente, los espacios protegidos y la ordenación del territorio, aunque de los mismos sitios hay información útil para otras muchas áreas de conocimiento.

Terrasit

Dentro de la Comunidad Valenciana, la web de referencia para conseguir mapas oficiales es Terrasit. Para poder bajarse los mapas es necesario registrarse. El procedimiento está explicado en el siguiente post:

Otros organismos oficiales

Por desgracia la información oficial está dispersa en las diferentes webs de los diferentes organismos oficiales y no existe un lugar centralizado donde conocer qué información está disponible en cada uno. No queda más remedio que bucear dentro de cada uno de ellos.

Asociaciones, Empresas y Organismos Internacionales.

Otras referencias

 

¿Conoces algún otro sitio de donde bajarse capas útiles?

 

Modificar registros de una capa vectorial en Qgis. Ejemplo para extraer el polígono de un municipio.

Para el ejemplo nos bajaremos la capa de municipios de la web de ArcGis.

Municipios ArcGisDando en “Descargar” se baja el archivo “Municipios_ETRS89_30N.zip”. Antes que nada hay que descomprimir el archivo en una carpeta que contendrá una serie de archivos que se llaman todos igual con diferentes extensiones:

instantánea6Cada uno de estos archivos son diferentes datos sobre la misma capa. Todos ellos juntos forman la capa vectorial. El archivo de referencia de la capa es el de formato “shape” y extensión “.shp”, mientras que el archivo con extensión .dbf es el que contiene la tabla de atributos que se podría abrir con cuidado en una hoja de cálculo, aunque no es aconsejable hacerlo directamente.

Abrir capa vectorial

Para abrir la capa vectorial en Qgis ir a Capa>Añadir capa>Añadir capa vectorial…

añadir capa vectorial Elegir el archivo de la capa en formato shape, elegir la codificación adecuada (normalmente UTF-8  o ISO-8859-1) y pulsar Abrir. Si no se elige la codificación adecuada, las ñ y los  acentos saldrán cómo símbolos extraños.

Una vez abierta la capa con todos los municipios de  España, pulsar con el botón derecho y abrir la tabla de atributos.

instantánea9En la tabla de atributos buscar el municipio que  queremos guardar. Para poder seleccionarlo debemos  activar la edición, lo seleccionamos y pulsamos a “invertir selección” de forma que quedan seleccionados todos los municipios menos el nuestro.

instantánea11Damos a borrar, desmarcamos editar tabla y guardamos la nueva tabla de atributos.

Guardar Capa

Una vez modificada la capa es importante guardarla. No sirve con guardar el proyecto de Qgis, es necesario guardar cada una de las capas del proyecto. Lo ideal es tenerlas ordenadas en una carpeta donde esté metido todo el proyecto de forma ordenada.

Para guardar la capa le damos encima con el botón derecho, damos a guardar como y le ponemos un nombre (en el ejemplo gatova.shp).

instantánea12Ya tenemos la superficie del municipio de Gátova separada del resto.

Guardar el proyecto

Si tenemos más capas y/o queremos guardar los colores y resto de cambios que hayamos hecho en “Propiedades”, es necesario guardar el proyecto. El proyecto guarda las rutas a cada una de las capas y los colores y otras propiedades de cada una de ellas.

Así como las capas en formato shape se pueden abrir en cualquier programa de GIS (ArcGis, gvSig, Qgis,…), el proyecto se guarda en un formato diferente para cada programa, por lo que no se puede abrir en un programa diferente del que se ha utilizado para crearlo.

Cómo instalar Qgis (actualizado)

Instalar Qgis es tan fácil como seguir los pasos de su página oficial:

https://www.qgis.org/en/site/forusers/alldownloads.html

Aunque conviene tener en cuenta que existen diferentes versiones y no siempre la última es la mejor. De hecho, para trabajar, la mejor suele ser la última versión LTR (long term release) que es la última versión estable.

INSTALAR QGIS EN LINUX

En los repositorios suele haber versiones antiguas, aunque funcionales. Es importante instalar además de qgis, el paquete de python (qgis-python) para que funcionen las herramientas y plugins. Para un uso básico probablemente tendrás suficiente con el que salga en el repositorio oficial, pero a veces el que aparece está bastante anticuado y hay que añadir el repositorio de qgis.

Para instalar desde el repositorio de ubuntu desde consola:

sudo apt-get install qgis python-qgis

En varias webs aparece el qgis-plugin-grass que puedes probar a instalarlo también, pero mejor por separado.

INSTALAR LA ÚLTIMA VERSIÓN LTS DE QGIS EN UBUNTU 14.04 LTS

Para tener la última versión LTS es necesario añadir repositorios, actualizar y después instalar. Está todo explicado en el enlace de arriba, pero para simplificar:

1. Editamos el archivo “/etc/apt/sources.list” para lo que hace  falta ser root. Se puede hacer con nano, gedit*, leafpad, kate o cualquier editor de  texto plano:

sudo gedit /etc/apt/sources.list

Agregamos las siguientes líneas al final del texto:
deb http://qgis.org/debian-ltr trusty main
deb http://qgis.org/ubuntugis-ltr trusty main
deb http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu trusty main

cambiando la versión de ubuntu a la que tengas instalada.

2. Activamos la clave de los repositorios. Volvemos a entrar a la consola y escribimos lo siguiente para activar la clave:
gpg --keyserver keyserver.ubuntu.com --recv DD45F6C3
gpg --export --armor DD45F6C3 | sudo apt-key add -

Si la firma da problemas prueba con esto.

3. Y ya sólo queda actualizar e instalar:
sudo apt-get update
sudo apt-get install qgis python-qgis

Una vez instalado correctamente, el siguiente paso es instalar los plugins que vayas a necesitar.

WEBS DONDE SOLUCIONAR PROBLEMAS DE QGIS

 

Referencias sobre botánica

A pesar de que mucha de la información es difícil de encontrar y por tanto de utilizar, sí que existen recursos gratuitos de fácil acceso que nos pueden ayudar a entender la botánica:

Geobotánica y vegetación

Muy aconsejables los apuntes de Francisco José Alcaraz Ariza, muy ordenados y actualizados (no sólo los temas de teoría sino toda la web):
http://www.um.es/docencia/geobotanica/teoria.html

Sobre la vegetación de España por zonas, una buena referencia clásica son los 18 volúmenes de Itinera Geobotánica:

https://floramontiberica.wordpress.com/itinera-geobotanica/

Sobre el estado actual de la situación de los mapas de vegetación:

Sainz et al 2010 La Cartografía sintética de los paisajes vegetales españoles: Una asignatura pendiente en Geobotánica. Ecología 23: 249-272.

Otra aproximación diferente son los Atlas de Idoneidad de especies leñosas de la Península Ibérica.  http://www.opengis.uab.cat/IdoneitatPI/index.html

Y para bajarse los mapas en shp:

Rivas-Martinez 1987 Mapas Series Vegetación: http://www.magrama.gob.es/es/biodiversidad/servicios/banco-datos-naturaleza/informacion-disponible/memoria_mapa_series_veg_descargas.aspx

En la Comunidad Valenciana:
Costa 1982. Pisos bioclimáticos y series de vegetación en el área valenciana. Cuadernos de geografía, ISSN 0210-086X, Nº 31, 1982 , 129-144

A escala europea:

Visor del European Forest Data Center y sus descargas.

Zona andina:

Atlas de los Andes  del Norte y del Centro

Flora

Para búsquedas sobre flora, lo mejor es ir directamente a las fuentes:

y los enlaces incluídos en ellos, de los que destacaría:

y las Bases de datos de biodiversidad:

Pero si sólo os interesan los árboles, lo mejor es  una App para el móvil:

  • ArbolApp. Una App para el móvil con los árboles de la Península Ibérica, Baleares y Canárias (ArbolappCanarias) editada por el CSIC. También tiene versión web en los enlaces anteriores.

Otros

En la parte normativa es importante tener en cuenta las listas-patrón de especies y hábitats publicadas por el Ministerio, con el objetivo de unificar su nomenclatura.

Para buscar más documentos de botánica, recomiendo la web de Jose Luis Benito, que tiene recopilados todos los enlaces importantes:

Y ya puestos, os dejo los enlaces a algunos libros interesantes de conservación de flora:

Si buscáis algo no gratuito igual os sirve este enlace.

Si conocéis más recursos gratuitos que merezcan la pena, podéis ponerlos en comentarios y los añadiré.

Automatizar el paso de GPX a SHP desde R

Desde la lógica más pura del trabajo con datos en R, lo mejor es tener los datos originales sin modificar, ya que siempre que se modifican es una posibilidad nueva de error. Para poder trabajar con esos datos originales lo que hacemos es convertirlos en una matriz de datos con los que podamos trabajar.

En nuestro caso tenemos en una carpeta los datos GPX tal cual salen del GPS y una tabla CSV con el mismo nombre que contiene los datos recogidos en la libreta (Especie, nombre parcela, otros datos) sobre los mismos puntos que contiene el GPX. Es importante que el número de registros en uno y en otro sean los mismos y que estén en el mismo orden, aunque no es necesario que se llamen igual, es interesante hacerlo para poder comprobar que se han unido bien.

Pasos a seguir:

  1. Instalar rgdal
  2. Elegir directorios de trabajo
  3. Importar datos GPX en R
  4. Importar datos CSV en R
  5. Juntar todos los datos en una única SpatialPointsDataFrame
  6. Corregir errores de transcripción
  7. Eliminar repetidos
  8. Exportar SHP

Instalar rgdal

Para importar datos GPX en R necesitaremos instalar la librería rgdal. Al intentar instalarla puede ocurrir que falte gdal-config, por lo que habrá que instalarlo primero en la consola de Linux:

sudo apt-get install libgdal1-dev libproj-dev

y en la consola de R:

install.packages("rgdal", dep=T)         #rgdal package
library(rgdal)

Una vez instalada, sólo es necesario correr la última línea de código para volver a activar la librería las siguientes veces.

Elegir directorios de trabajo

Lo primero que vamos a hacer es fijar los directorios de trabajo:

#   Directorio de entrada de datos
DirDatosBrutos <- "Directorio donde están los datos GPX y CSV/"
#   Directorios para exportar datos
DirDatosSHPwp  <- "Directorio donde se querran guardar los archivos shp de puntos/"
DirDatosSHPtrk <- "Directorio donde se querrán guardar los shp de tracks/"

Observar que al principio no hay barra y al final hay una barra inclinada (en Mac o en Windows puede ser diferente). Para comprobar que están bien elegidos, podemos utilizar:

list.files()
list.files(DirDatosBrutos)
list.files(DirDatosBrutos, ".csv")
list.files(DirDatosSHPwp, ".shp")
list.files(DirDatosSHPtrk, ".shp")

Importar datos GPX en R

Para importar los datos GPX y convertirlos en variables espaciales de puntos (waypoints) o lineas (tracks) utilizaremos readOGR()

#=============================================================================== 
#          LEER DATOS GPX y csv y crear listas
#=============================================================================== 
GPXinNombres <- list.files(DirDatosBrutos, ".GPX")  # Lista nombres gpx importar
CSVinNombres <- list.files(DirDatosBrutos, ".csv")  # Lista nombres csv importar
#
#                         GPXwp           Lista de SpatialPointsDataFrame PUNTOS
GPXwp <- list()
for (i in GPXinNombres){
    nombre <- paste(DirDatosBrutos,i,sep="")
    GPXwp[[i]] <- readOGR(dsn = nombre, layer="waypoints")
}
#
#                       GPXtrk             Lista de SpatialLinesDataFrame TRACKS
GPXtrk <- list()
for (i in GPXinNombres){
  nombre <- paste(DirDatosBrutos,i,sep="")
  GPXtrk[[i]] <- readOGR(dsn = nombre, layer="tracks")
}
#                       COMPROBAR
for (i in 1:length(GPXwp)){
   plot(GPXwp[[i]])
   plot(GPXtrk[[i]],add=T)
}
# Debería salir un dibujo con los traks y waypoints de cada archivo gpx

Importar datos CSV en R

Lo ideal es que tengan el mismo nombre que los correspondientes GPX, pero si no, al menos deberían estar en el mismo orden alfabético. También los nombres de las variables deberían ser en todas las mismas. En el ejemplo son ID, Tipo, Lugar Parcela y Fecha, pero habrá que poner en cada caso los datos tomados en la libreta.

#                         CSVwp     Lista de matrices de datos espaciales
CSVwp <- list()
for (i in CSVinNombres){
nombre <- paste(DirDatosBrutos,i,sep="")
CSVwp[[i]] <- read.table(nombre,header=T,sep=",")      #,sep="\t")     #   Leer tabla
names(CSVwp[[i]]) <- c("ID","Tipo","Lugar","Parcela","Fecha")          #   Arreglar nombres variables
}
#

Limpiar variables del GPX

Podemos ver todas las variables introducidas en cada una de las tablas utilizando str() o head(). En el ejemplo sólo nos quedaremos con las variables 5, 1 y 2, “nombre”, “ele” y “time”, respectivamente, eliminando todas las demás ya que están vacías o no tienen información relevante. Comprobar siempre antes, ya que pueden variar en función del modelo de gps. Para eliminarlas de toda la lista de tablas a la vez, utilizaremos for().

#===============================================================================
#          LIMPIAR GPXs
#===============================================================================
#          Seleccionar las variables útiles y crear variable VarUtiles
head(GPXwp[[1]])
VarUtiles <- c(5,1,2)                           # Puede cambiar según modelo GPS
for(i in 1:length(GPXwp)){
GPXwp[[i]] <- GPXwp[[i]][,VarUtiles]
}
head(GPXwp[[1]])                                             # Comprobar cambios
#

Juntar las variables de los GPX con las de los CSV

Un detalle importante en este paso es que los nombres de las variables en el CSV no deben coincidir con las variables que se hayan conservado del GPX. Si es así debería cambiarse el nombre de una de las dos. En este momento tenemos una lista de gpx con tres variables cada uno y una lista de csv con 4 variables cada uno. Podrían ser más y no pasaría nada mientras estén en el mismo orden y lógicamente el número de puntos en los gpx coincida con el número de lineas de datos de cada csv. Si es así, esto debería funcionar.

#===============================================================================
#          JUNTAR GPX Y CSV
#===============================================================================
for(i in 1:length(GPXwp)){
for (j in names(CSVwp[[i]])){
GPXwp[[i]][[j]] <- CSVwp[[i]][[j]]
}
}
#                       COMPROBAR si se han añadido
head(GPXwp[[1]])

Juntar todos los datos en una única SpatialPointsDataFrame

para juntar todos los datos en una única tabla de datos espaciales “SpatialPointsDataFrame” creamos un objeto introduciremos la primera tabla y luego con for() vamos añadiendo las demás utilizando rbind() para añadir cada vez más filas. Más tarde quitaremos los elementos duplicados por lo que no debemos  preocuparnos si los archivos solapan o si se tomaron puntos repetidos.

WP <- GPXwp[[1]] # Crea tabla GPXWP
for(i in 2:length(GPXwp)){
WP <- rbind(WP,GPXwp[[i]]) # Une resto de tablas a GPXWP
}
# COMPROBAR
WP

Corregir errores de transcripción

Además de comprobar que se han introducido todos los casos es interesante comprobar también si existen errores de transcripción que haya que corregir. Para las variables cuantitativas suele ser interesante ver los valores máximos y mínimos. Las cualitativas se pueden comprobar viendo los niveles de cada factor.

levels (WP$Tipo)                                 # Debería haber sólo 4: A,H,R,S

En nuestro ejemplo resulta que en la variable Tipo se han introducido en algunos casos “Ref” en vez de “R” y “X”  en vez de “H”. Para corregirlo

levels(WP$Tipo)[levels(WP$Tipo)=="REF"] <- "R"   # Cambia REF por R
levels(WP$Tipo)[levels(WP$Tipo)=="X"] <- "H"     # Cambia X por H
levels (WP$Tipo)                                 # Debería haber sólo 4: A,H,R,S

Eliminar repetidos

Los puntos que coinciden en coordenadas y en Tipo se consideran puntos repetidos y por tanto serán eliminados para que en cada coordenada sólo pueda haber un punto de cada tipo.

#===============================================================================
#          ELIMINAR REPETIDOS
#===============================================================================
coordTipo <- cbind(coordinates(WP),WP$Tipo)           # Variables para repetidos
dupl <- duplicated (coordTipo)                        # casos duplicados
#                       COMPROBAR
# length(subset(dupl,dupl==TRUE))                     # Numero casos duplicados
# length(subset(dupl,dupl==FALSE))                    # Numero casos sin duplicar
#
str(WP)
WP <- WP[!dupl,]                                      # ELIMINAR CASOS REPETIDOS
str(WP)
#

Exportar SHP

Ahora ya podemos exportar un shp con todos los puntos de todos los gpx y con los datos del csv incluidos en la tabla de atributos utilizando la función writeOGR().

writeOGR(WP, dsn = DirDatosSHPwp, "WP-Definitivo", driver = "ESRI Shapefile", overwrite_layer = T)  # Crea WP.shp DEFINITIVO

También se puede exportar en otros formatos (GPX, CSV, etc.). Si además queremos exportar una capa de shp para cada uno de los archivos de traks de cada gpx, se puede hacer así

for(i in 1:length(GPXtrk)){
writeOGR(GPXtrk[[i]],DirDatosSHPtrk,names(GPXtrk)[i],driver="ESRI Shapefile",overwrite_layer = T)
}

No he conseguido (tampoco lo he necesitado por ahora) juntar todas las “SpatialLinesDataFrame” de las tracks de los gpx en una sólo ni eliminar los segmentos repetidos. Si en algún momento veo cómo hacerlo lo añadiría aquí. Si sabes cómo hacerlo puedes comentarlo.

Cómo pasar los datos de tu GPS a capas de GIS en shape utilizando Qgis

La mayor parte de la información cartográfica pública se publica en formato shp, por lo que suele ser la mejor forma de editar, pasar a limpio y almacenar nuestros datos de distribución espacial. Los datos que recogemos con el receptor GPS se guardan normalmente (o al menos se pueden exportar) en formato GPX. Los datos de cada punto se suelen meter en libreta y de ahí a una tabla (o directamente en una tablet), pero no suele ser cómodo hacerlo en el mismo aparato, por lo que luego hay que juntar los datos. Vamos a ver cómo juntarlos en formato shp para poder utilizar otras capas de información geográfica y en futuros posts veré la forma de pasar los datos shp a R para poder realizar análisis de distribución espacial.

Si lo único que se quiere es visualizar o modificar los puntos o tracks y volver a guardarlos en gpx o volver a meterlos en el gps, lo mejor sería no complicarse la vida y utilizar Qlandkartegt o su nueva versión Qmapshack.

En este caso vamos a complicarlo un poco más que eso. Partimos de la base de que tenemos un archivo gpx con puntos y tracks y una tabla con los datos de cada punto en .csv (se puede exportar a .csv desde ods, xls o cualquier otro formato). Vamos a juntar ambos con Qgis y guardarlo todo en formato shp, o incluso volver a guardarlo como gpx. Pero mejor, paso a paso.

Si tienes muchos archivos gpx y csv, igual es mejor automatizar el proceso y juntar todos los archivos en un único shp utilizando R.

Primer paso: Abrir gpx en Qgis

Se pueden abrir los archivos gpx en Qgis clicando Capa>Añadir capa>Añadir capa vectorial…

qgis1

Se selecciona el archivo y se abren dos capas, la de waypoints (puntos) y la de tracks (líneas). Obviamente si se quieren guardar el resto de datos también se puede, pero no suelen ser interesantes.

qgis2 qgis3 qgis4

Paso 2: Comprobar que está todo en su sitio

Para comprobar que está todo en su sitio la mejor opción es utilizar alguno de los plugins de mapas web aconsejados (QuickMapServices, OpenLayers) y abrir de fondo algún mapa (osm, gmaps) para ver si los puntos están colocados donde se supone que deben estar. Si no está donde debe lo más probable es que haya que ajustar la proyección.

qgis5

Paso 3: limpiar tabla de atributos

Antes de añadir nuevos datos y marcas en la tabla de atributos es interesante ver qué es lo que hay ahí dentro y quitar todo lo que no sea útil. De todas las variables que introduce el gps en la tabla de atributos,  en nuestro caso la única información  útil es el nombre de cada waypoint, ya que es lo que vamos a utilizar para identificarlos y unir los datos de la tabla csv. Para quitar el resto de variables de la tabla de atributos hay que darle primero a editar y luego a eliminar columnas.

Importante: hay que haber guardado previamente la capa para poder editarla. Para ello: botón derecho>guardar como

qgis6 qgis7

 

 

 

 

 

 

 

En otros casos puede interesar conservar “ele” elevaciones, o “time” el momento en que se realizó el waypoint.

Un detalle importante es fijarse en que el nombre de los waypoints en este caso empiezan en 005. Los  ceros a la izquierda cuentan en el nombre, por lo que tendremos que ponerlos igual en la tabla csv. Para ello es importante ver en formato de número y añadir números a la izquierda antes de exportarlo como csv.

Paso 4: añadir datos a los waypoints a partir de una tabla de datos csv

Lo más importante que tiene que tener la tabla es una de las variables igual a alguna de las variables de la tabla de atributos de la capa. En el ejemplo la variable ID del archivo .csv es el nombre del punto, escrito exactamente igual que la variable “nombre” de la tabla de atributos, ceros a la izquierda incluidos.

qgis8

Para añadir los datos utilizamos el plugin MMQgis, que ya describí anteriormente.

qgis9qgis10

Importante fijarse en cada dato que te pide el formulario:

  • elegir el archivo csv,
  • poner el nombre de la variable,
  • la capa vectorial
  • y el nombre de la misma variable en la capa, para que los datos se unan adecuadamente.
  • También es interesante poner un nombre de archivo shape con el nombre de la nueva capa (ejemplo.shp)
  • y de los datos que no se han unido por si acaso (si todo va bien debería quedar vacío).

 

qgis11qgis12qgis13qgis14

Y como siempre que se trabaja con GIS, que no se te olvide guardar, en dos fases:

  • Cada capa por separado (formato shp, botón derecho y guardar como) para guardar los cambios en las tablas y situación de los puntos.
  • Y el proyecto (formato qgs) para guardar las localizaciones de cada capa, el orden y los colores elegidos.

Paso 5, o cerrando el bucle: Volver a guardar de shp a gpx

En caso necesario podéis volver a guardar la capa que queráis como gpx para volver a introducirla en vuestro gps, pero con los nuevos datos. Para ello botón_derecho, guardar _como,  gpx.

 

Utilizar la capa de Open Street Maps (OpenLayerPlugin) para comprobar el sistema de coordenadas en Qgis

Como ya había comentado en un post anterior, el plugin de OpenLayerPlugin es muy útil para situarnos donde estamos y comprobar que nuestras capas de mapas están encuadradas correctamente en el territorio. Por ejemplo, si tenemos mal puesto el sistema de coordenadas, puede ocurrir que el municipio de Gátova aparezca totalmente descolocado (al lado de Menorca, por ejemplo), un poco descolocado como en el mapa de debajo a la izquierda.

instantánea13

14-15

Sólo con el sistema de coordenadas correcto la figura se verá perfectamente alineada con los límites de la capa de OpenStreetMaps.

En definitiva, este plugin sirve para situarse y ver si nuestros datos están situados correctamente.

Plugins imprescindibles en Qgis

Después de instalar el Qgis, es importante echar un vistazo a los plugins. Una de las grandes ventajas de Qgis es su gran versatilidad que consigue a partir de los plugins que añaden nuevas características al programa.

Para instalar un plugin basta con ir a Complementos>Administrar e instalar complementos…

Qgis Instalar complemento

Y una vez en el administrador de complementos buscar, marcar e instalar los que uno quiera.

Qgis administrador complementos

Existen plugins casi para cualquier cosa, aunque la dificultad está en saber cuáles son realmente útiles. Un primer vistazo a la página de plugins de qgis, ordenada por los mas votados (Most voted) o los más bajados (Top downloads) nos dice ya cuales son los más utilizados.

Qgis most downloaded plugins

De estos, los que más vamos a utilizar son:

OpenLayers Plugin: Permite abrir y utilizar capas de OpenStreetMap, Google Maps, Bing Maps, MapQuest y Apple Maps. Estas capas son muy útiles para comprobar que estamos utilizando la proyección correcta en nuestras capas, entre otras cosas.

OpenLayer tiene el problema reconocido de que no funciona bien al imprimir el mapa, por lo que cuando se va a imprimir es mejor utilizar los mapas de QuickMapServices, que funcionan perfectamente, aunque sólo utiliza mapas libres (no google).

mmqgis: Nos permite exportar la tabla de atributos como CSV y luego volver a importarla, de forma que podemos trabajar datos de la tabla en cualquier hoja de cálculo (Calc, Excel,..) de forma mucho más cómoda que en el editor de tablas.

Table Manager: Nos permite ordenar las columnas en la tabla de atributos.

Más información:

¿Has utilizado algún otro plugin que consideres imprescindible? Ponme la información en un comentario y lo probaré.

Qgis o gvSig

Difícil decisión. Personalmente siempre he encontrado a Qgis más amigable con el usuario y a gvSig más difícil de utilizar. Sin embargo el hecho de que la documentación de gvSig sea en español le da algo de ventaja. Probablemente por la misma razón Qgis es el más utilizado en organismos internacionales, mientras que gvSig es el más utilizado en las administraciones españolas.

No parece fácil decidir:

Para instalar Qgis:

http://garmendia.blogs.upv.es/instalar-qgis/

Para instalar gvSig:

http://garmendia.blogs.upv.es/instalar-gvsig/

Instalar gvSig

gvSig tiene documentación en español, por lo que no debería ser difícil de instalar:

http://www.gvsig.com/es/productos/gvsig-desktop/descargas

Me sorprende que todavía no se encuentre en ningún repositorio oficial, pero así es, hay que bajarse el paquete binario de la web. Si se intenta instalar directamente no funcionará debido a que es un programa que corre sobre JAVA, lo que significa que tienes que tener java instalado previamente.

 

INSTALAR GVSIG EN UBUNTU 14.04 LTS

  1. Instalar Java
  2. Bajarse gvsig
  3. Instalar gvsig

1.Instalar Java

sudo apt-get update
sudo apt-get install openjdk-7-jdk:i386 libidn11:i386 libldap-2.4-2:i386 libjpeg62:i386 libproj0:i386 libgeos-c1:i386 libssl0.9.8:i386http://garmendia.blogs.upv.es/wp-admin/post.php?post=39&action=edit

2.Bajarse gvsig

Ahora sí, elegir desde:

http://www.gvsig.com/es/productos/gvsig-desktop/descargas

3. Instalar gvsig

Dar permiso de ejecutar como programa y correr el programa de instalación.

Seguir los pasos de la documentación.

Manual de gvSig:

http://downloads.gvsig.org/download/web/es/build/html/index.html

 

 

Calculadora de campos de Qgis para poner valores a partir de cualidades (ejemplo Aptitud)

La tabla de atributos de Qgis todavía es una de las asignaturas pendientes en su desarrollo. No permite arrastrar o copia y pegar muchos valores a la vez, de forma que cuando tenemos muchos registros se complica realizar operaciones. Una forma posible es abrir el archivo .dbf en una hoja de cálculo y luego tener cuidado de volver a ponerla en su sitio. Se puede hacer pero no resulta lo más aconsejable. En gvsig se permite exportar la tabla en .CSV y luego juntar campos, lo que facilita esta operación. mmqgis es un pluging de Qgis que permite exportar la tabla de atributos como CSV y luego volver a importarla, de forma que podemos trabajar datos de la tabla en cualquier hoja de cálculo (Calc, Excel,..) de forma mucho más cómoda que en el editor de tablas. Aquí se explican otros plugins interesantes como Table Manager, que nos permite ordenar las columnas en la tabla de atributos.

En cualquier casos en la mayoría de los casos se puede conseguir casi cualquier cosa utilizando la calculadora de campos, que cada vez es una herramienta más versatil.

Para ver todos los detalles de la calculadora de campo, vienen muy bien explicados aquí y aquí, por lo que no me entretendré en eso. Sino que directamente paso a ver cómo

PONER VALORES DE APTITUD A PARTIR DE CUALIDADES

Para ello utilizaré la capa Siose 2011 del municipio de Gátova, bajada de Terrasit (aquí explico cómo) y simplificada eliminando las columnas que no necesito (4 columnas x 193 registros). En la figura los primeros 39 registros:

Screenshot - 031215 - 13:38:35

La idea es crear una columna con valores numéricos para cada tipo de vegetación según la considere más o menos apta para una actividad concreta.

Para ello editamos la tabla (botón de más a la izquierda) y abrimos la calculadora de campos (botón de más a la derecha). También se pueden ordenar los registros por orden alfabético presionando sobre el nombre.

Los valores de aptitud que queremos poner a cada valor de NIVEL 1 son:

Arbolado forestal 2
Matorral 2
Pastizal 1
Cultivos 1
Coberturas artificiales 0
Terrenos naturales sin vegetación 0

Obviamente hay muchas formas de obtener el mismo resultado, voy a exponer las dos que más me gustan, que son con ILIKE y con CASE_ELSE.

A mi me gusta más la primera porque es más fácil detectar errores y las columnas extra me permiten realizar cálculos desde la hoja de cálculo (sólo tengo que copiar el archivo .dbf en otra carpeta y abrirlo.

Detalles a tener en cuenta en las operaciones:

  • las variables se ponen entre comillas dobles
  • El texto se pone entre comillas simples
  • Los valores numéricos se ponen tal cual
  • Conviene escribir la fórmula en un editor de texto por si no sale bien poder corregirla usando “Actualizar campo existente” y para luego saber cómo se ha hecho (en Qgis no se guarda)·

Utilizando ILIKE

En realidad utilizar LIKE, ILIKE o = es casi lo mismo. Para valores numéricos es lo mismo, aunque en ese caso lo mejor es utilizar =, <,>. Para texto ILIKE tiene la ventaja de no discriminar entre mayúsculas y minúsculas.

La mejor opción desde mi punto de vista es dividir el problema en problemas menores fácilmente revisables, haciendo una columna para Aptitud 2 (A2), otra para Aptitud 1 (A1) y si llegara a hace falta otra para A0.

Para el caso de A2 se puede hacer con OR:

"NIVEL1" ILIKE 'Arbolado forestal'
OR
"NIVEL1" ILIKE 'Matorral'

Para la A1 se hace de la misma manera, cambiando los valores por ‘Pastizal’ y ‘Cultivos’.

Luego, una vez creadas las dos columnas, y comprobado que no queda ningún valor sin utilizar, simplemente hay que operar con ellas multiplicando a cada una por su valor, con lo que nos queda una columna con los valores de aptitud de la capa Siose 11, utilizando los valores de la variable NIVEL1 (Apt_S11N1):

"A1"+2*"A2"

En realidad la columna A0 no hace falta en este caso, pero en ocasiones es útil para comprobar que no nos dejamos nada.

Utilizando CASE_ELSE

Si tengo claro que no voy a tener que corregir valores ni a operar desde hoja de cálculo, puedo hacerlo directamente con el Condicional CASE ELSE

CASE WHEN condition THEN result ELSE result END

Screenshot - 031215 - 15:32:06es decir:

CASE
    WHEN
        "NIVEL1" ILIKE 'Arbolado forestal'
        OR
        "NIVEL1" ILIKE 'Matorral'
    THEN
2
    WHEN
        "NIVEL1" ILIKE 'Pastizal'
        OR
        "NIVEL1" ILIKE 'Cultivos'
    THEN
1
    ELSE
0
END

Suele ser buena idea hacerlo de dos formas diferentes y comprobar que los datos que salen son los mismos, sobre todo cuando hay tantas filas que es imposible comprobarlo de forma manual.

IMPORTANTE: Escribir la fórmula en un editor de texto o copiarla y pegarla en un editor de texto antes de dar a aceptar, de forma que si no sale perfecta la podamos modificar y volver a actualizar el campo.

Referencias

http://docs.qgis.org/2.2/pt_BR/docs/user_manual/working_with_vector/field_calculator.html
http://mappinggis.com/2014/05/como-usar-la-calculadora-de-campos-de-qgis/
http://garmendia.blogs.upv.es/bajarse-capas-terrasit/

 

 

Bajarse capas shp y raster de Terrasit (Comunidad Valenciana)

Para poder bajarse las capas de mapas de Terrasit es necesario registrarse previamente. Sólo permite bajarse mapas de cada municipio por separado, por lo que es necesario ir municipio a municipio. Lo que sí permite es ver las capas en su visor web antes de bajarlas.

Además, gran parte de la información ambiental que no está disponible en formato shape ni raster, está disponible en wms, por lo que se puede al menos consultar. Vamos por partes:

  1. REGISTRARSE
    1. Ir a la web de Terrasit: http://terrasit.gva.es/es
    2. Entrar en la pestaña >descargas>descargas>acceso a registro>Crear cuenta nueva
    3. Introducir datos personales para registrarse. Tienes derecho a mentir (como siempre que te pidan datos personales), pero es importante acordarse del usuario y la contraseña ya que luego las necesitarás.
  2. INICIAR SESIÓN
    1. Ir a la web de Terrasit: http://terrasit.gva.es/es
    2. Entrar en la pestaña >descargas>descargas>acceso a registro>Iniciar sesión y meter el nombre de usuario y contraseña.
  3. DESCARGAR CARTOGRAFÍA
    1. Ir a la web de Terrasit: http://terrasit.gva.es/es
    2. Entrar en la pestaña >descargas>cartografía (si no has iniciado sesión no sale)
    3. Seleccionar municipio u hoja a descargar. Las hojas son ortofotos por lo que no nos sirven para combinar los datos, utilizaremos los datos de un municipio (ejemplo: Gátova).

Terrasit descargas gátova

nos bajamos las capas que vayamos a utilizar presionando en la flecha de la derecha. Hay que acordarse de la proyección para luego ponerla bien en Qgis. Algunos de los mapas que pueden resultarnos útiles son:

– Siose 2011 (Usos del suelo y tipos de vegetación)

– Infraestructura verde (Lics, cuevas, parques y parajes, etc.)

– Planeamiento (propiedad y protección suelo)

– Cartografía tridimensional (pendientes, orientaciones, etc)

Si se está haciendo un trabajo de un municipio concreto merece la pena bajarse todas y ver qué información contiene cada una de ellas. Hay capas vectoriales y capas raster.

Todas esas capas y algunas más las podemos visualizar previamente por wms o en el visor de terrasit.

Para conseguir mapas de otros orígenes, ver:

http://garmendia.blogs.upv.es/conseguir-capas-de-mapas/