lunes, 19 de febrero de 2018

SASPlanet, algunos mapas nuevos.

Hacía mucho que no dábamos un repaso a nuestra colección de mapas para SASPlanet. Lo cierto es que nos es imposible mantener controlados tantos mapas, y nos llevaría un tiempo que no tenemos comprobar aquellos que siguen funcionando o aquellos que por cualquier cambio o baja en los servicios WMS han dejado de funcionar. Como siempre pedimos vuestra colaboración y si algún mapa os interesa dejad un comentario y miraremos de intentar solucionarlo.
Aún así nos gusta de vez en cuando introducir en nuestro lote aquellas novedades que nos parecen muy interesantes y útiles. Luego ya es cosa personal de cada uno colocarlos en su colección (o no), o configurarlos de la forma que estime oportuno. En el blog hay información sobre todo ello. Y para aquellos para los que Google satélite es lo único importante en el mundo la última versión que conocemos a fecha de hoy es la 762. En esta ocasión los elegidos son los siguientes:

SPAIN

- PNOA Histórico: Completamos esta colección de CAPAS en el menú Spain con 4 nuevas incorporaciones: las imágenes del vuelo PNOA 2016, el vuelo OLISTAT (1997-1998), el vuelo Nacional_1981-1986 y el vuelo Interministerial_1973-1986. De todas ellas tenéis información en nuestros artículos sobre el WMS en cuestión.
Ejemplo de un par de vuelos del WMS ortofotos históricas activados como capas.
- Naútica IHM: se trata de las Cartas Naúticas de la Armada, servidas por el Instituto Hidrográfico de la Marina (IHM). Hemos incluido sus cuatro niveles de detalle, ya que están en servicios WMS individuales imposibles de ofrecer agrupados en uno. Nos parecen una maravilla para todos los amantes del mar y la navegación.
 Ejemplo sobre la ría de Camariñas (Galicia) a zoom 16 de SASPlanet del nivel de detalle 2 (arriba) y del nivel de detalle 4 (abajo).


CASTILLA Y LEÓN (CyL)

- PRUG: incorporada como capa en CyL la cartografía del Plan Rector de Uso y Gestión (PRUG) del Parque Nacional de Guadarrama, que actualmente se encuentra en sus últimos días de alegaciones. Se acompaña Leyenda. Ejemplo a zoom 16 de SASPlanet.


NAVARRA

- Orto MA: ortofoto Máxima Actualidad del servicio navarro que ofrece siempre las últimas imágenes disponibles (en este caso la reciente orto de 2017, 25 cms de resolución). Captura a zoom 19 de SASPlanet.

- MTNa_5k_2014: nuevo mapa topográfico 1:5000 ráster de Navarra a partir de su BTA vectorial. Un salto cualitativo frente al anterior y ya muy envejecido 5k (también presente en nuestro SASPlanet). Captura a zoom 19 de la misma zona de la orto anterior.


VALENCIA

- SENTINEL 2017: imagen Sentinel 12/11/2017 de la Comunidad Valenciana en color verdadero (RGB) y 10 metros de resolución. Captura a zoom 15.

- Orto 1956-57: ortofoto 50 cms B/N de la Comunidad Valenciana, año 1956. Captura a zoom 17.

- Orto 2017: ortofoto de la Comunidad Valenciana a 25 cms, año 2017. Captura a zoom 17.


BALEARES

- Relieve: modelo de relieve sombreado de las Islas Baleares. Captura a zoom 13.

Por lo demás hemos aprovechado para corregir algunos mapas que mostraban mal o ningún funcionamiento (otros han quedado para otra ocasión), así como para suprimir algunos que no hemos sido capaces de recuperar o que han sido absorbidos por otros servicios (por ejemplo, CartoCiudad absorbido por el Mapa Base).
Como siempre nuestro lote SASPlanet con mapas incorporados listo para funcionar, disponible en la sección DESCARGAS del blog. Saludos y hasta la próxima revisión.

lunes, 12 de febrero de 2018

Caso práctico LIDAR-QGIS: el campamento romano de Mérida.

Hace unos días dimos cuenta en nuestro Facebook de una impactante noticia: la localización el año pasado gracias a los datos LIDAR de un campamento legionario (castra) en las orillas del río Lácara, un afluente del Guadiana, cerca de Mérida. Podéis leer el artículo original del trabajo para más detalles. Y no es el único caso de descubrimientos arqueológicos usando estos datos.
Concretamente la imagen superior que acompaña el estudio, y en la que se observa la ubicación del campamento nos dejó absolutamente intrigados: en la imagen a la ortofoto PNOA de la zona, sin rastro visible ni para los ojos más avezados. En la imagen b el sombreado a partir de los datos LIDAR mostrando claramente el perímetro del supuesto campamento romano.
Obviamente la intriga y la curiosidad nos pudo una vez más. Y aunque no tenemos los conocimientos suficientes para hacer estudios semejantes, ni tampoco sabemos los procesos detallados empleados en la imagen (se mencionan los filtros Laplacian y Gaussian), no perdíamos nada por trastear un poco con los datos LIDAR de la zona.

Para todo el proceso vamos a utilizar las herramientas LAStools para datos LIDAR, compatibles con QGIS, y que podéis integrar en el programa siguiendo este artículo.

Esto es lo que hicimos a nuestra manera amateur:
1.- Descargar el archivo LIDAR de la zona desde el Centro de Descargas. Tras buscar la ubicación, descargamos su cuadrícula LIDAR 2x2 kms. Concretamente se trata del archivo PNOA_2010_LOTE9_EXT_722-4328_ORT-CLA-COL.LAZ (13,6 megas). La extensión .laz es la clásica de los datos LIDAR comprimidos. Confirmamos que el archivo es el correcto con un vistazo rápido en la primera herramienta LAStool: lasview.
Habilitamos la caja de herramientas de QGIS, y dentro de las Herramientas para datos LIDAR, y del conjunto LAStools, ejecutamos lasview. Simplemente cargamos nuestro archivo LAZ en la ruta y ejecutamos Run para arrancar el visor.
La nueva ventana de visualización efectivamente nos confirma que el archivo cubre la zona que necesitamos. La cerramos.
2.- Extraemos el DEM desde los datos LIDAR. Para ello vamos a ejecutar la utilidad blastAdem de LAStools, que genera un archivo de elevaciones estándar desde la nube de puntos LIDAR.
Ejecutamos la herramienta blastAdem. En la ventana de configuración cargamos nuestro archivo LAZ. En las opciones de filtadro (verde) es MUY IMPORTANTE seleccionar los puntos class 2, que son los catalogados como ground en LIDAR. De esta forma eliminamos vegetación, edificios, ruido, etc... y trabajamos solamente con los puntos clasificados como "tierra". En la categoría Producto podríamos seleccionar directamente hillshade ahorrando un paso pero lo dejamos por defecto para así mostraros algún proceso más de QGIS.
Tan importante es filtrar los puntos ground que si no hacemos filtrado el DEM resultante es este:
Y aplicando el filtro de puntos ground (class 2) el resultado del DEM es este otro:

3.- Generando sombreado hillshade sobre el DEM. La estructura del posible campamento ya se intuye en la imagen, simplemente haciendo un DEM filtrado sobre los puntos ground LIDAR, limpios del resto de clasificaciones. Digamos de manera tosca que sin filtrado tenemos un modelo de elevaciones y con filtrado un modelo de superficie.
Ahora vamos a generar el modelo de sombras para darle el aspecto 3D necesario que remarque el relieve. Desde el menú Ráster/Análisis seleccionamos MDT (modelos de terreno):


En la ventana de configuración de los parámetros ya tenemos por defecto cargada la capa DEM (es la única en el proyecto -en naranja-). Seleccionamos una ruta y formato de salida para el hillshade generado (yo he escogido formato tif -en verde-). 


En la pestaña Modo -amarillo- escogemos obviamente Mapa de sombras (Hillshade). 


Y los parámetros marcados en rosa deciden la exageración vertical -yo he escogido 2 para remarcarlo un poco más- la escala y la orientación y altura de la luz solar con la que se van a generar las sombras del terreno.






Obviamente con todas estas opciones os recomendamos probar y observar los diferentes resultados, pues son muy distintos de aspecto dependiendo los parámetros escogidos. Dejamos marcada la opción Cargar en la vista del mapa para que la capa tif se añada al proyecto.



Una vez todo colocado pulsamos Aceptar y en breves segundos se nos genera el mapa de sombras.

Ahí tenemos nuestro campamento romano perfectamente marcado donde parecía no haber nada... Alucinante.
El mapa de sombras también podemos hacerlo directamente sobre la capa DEM, abriendo sus propiedades, sin necesidad de generar nueva capa (aunque yo recomiendo con capa aparte). En la pestaña Estilo disponemos de todas las opciones para probar.

Y hasta aquí esta magia de los procesamientos de imagen que no deja de sorprendernos con nuevas aplicaciones. Aunque de modo seguramente menos profesional, pero también con software completamente gratuito, lo contentos que nos hemos quedado haciendo aparecer el cuadradito romano... Saludos y hasta la próxima.

miércoles, 7 de febrero de 2018

Nuevas capas en WMS ortofotos históricas.

Hace un año escribíamos el artículo referente al impresionante servicio WMS de ortofotos históricas del IGN, imprescindible siempre que se hable del historial de las imágenes aéreas de España. Ahora este producto, de una utilidad y demanda incuestionables, vuelve a actualizarse para incluir tres nuevas capas de información: los vuelos PNOA 2016, Nacional (1981-1986) e Interministerial (1978-86).
Recordando que la URL de acceso es http://www.ign.es/wms/pnoa-historico? vamos a analizarlos y dar algunos detalles de estas nuevas incorporaciones.
Listado de capas del servicio WMS. En verde capas añadidas en esta actualización.
1.- PNOA 2016.
Se incluye el vuelo del año 2016 correspondiente al Plan Nacional de Ortofotografía Aérea (PNOA), con resoluciones de entre 25 y 50 cms, que abarca las Comunidades Autónomas de Andalucía, Murcia y Extremadura. Si estás interesado en los detalles técnicos de este vuelo puedes consultar sus especificaciones.
Además de la capa global que nos permite conocer los distintos vuelos PNOA en cada punto, la capa de cada vuelo cuenta con la capacidad GetFeatureInfo que nos permite consultar la resolución y la fecha de la ortoimagen.

2.- Nacional (1981-1986).
Sus metadatos nos dicen textualmente: Vuelo fotogramétrico en blanco y negro de ámbito nacional realizado por encargo del Instituto Geográfico y Catastral (actual Instituto Geográfico Nacional). Fechas de vuelo de 1980 a 1986. Fotogramas en blanco y negro a escala 1:30.000. Tamaño de pixel (GSD) de 75 cm y 45 cm, según zonas. Sistema geodésico de referencia ED50 en la Península, Islas Baleares, Ceuta y Melilla, y REGCAN95 en las Islas Canarias y proyección UTM en el huso correspondiente. Las Comunidades Autónomas de Extremadura y Baleares se ofrecen en ETRS89 y UTM en el huso 30 y 31, respectivamente.
En esta actualización se han generado las coberturas de las Comunidades Autónomas de Andalucía, Extremadura, La Rioja, País Vasco y de las Islas Baleares. La capa cuenta también con GetFeatureInfo.

3.- Interministerial (1978-86).
Citando al Blog IDEE, el Vuelo Interministerial se realizó entre noviembre de 1976 y julio de 1986 entre el Ministerio de Agricultura (65,25 %), MOPU (25,07 %), Hacienda (8,34 %) y el Ministerio del Aire (1,34 %) a una escala aproximadamente de 1:18 000 y constituye el primer caso de colaboración entre organismos para la realización de vuelos fotogramétricos
En el siguiente enlace del Mapama disponemos del mapa con las zonas responsabilidad de cada cual:
Las resoluciones varían entre 0,25 y 0,5 m en blanco y negro. Se ha generado ya la cobertura de las Comunidades Autónomas de Andalucía, Extremadura, Castilla y León, La Rioja y País Vasco y más adelante se irán incluyendo el resto de autonomías. Esta capa cuenta también con la habilidad GetFeutureInfo.

Se prevé que las ortofotos de estos vuelos estén pronto disponibles para descarga en la sección correspondiente del Centro de Descargas del CNIG. También pueden consultarse los fotogramas, imprimirlos o pedir certificados en la Fototeca Digital.
Muchas gracias de nuevo a nuestro IGN por seguir ampliando contenidos con regularidad y profesionalidad. Nosotros ya lo tenemos todo colocadito y disponible en SASPlanet.

martes, 6 de febrero de 2018

¿Las bicicletas no son para Guadarrama?

Aunque en varias ocasiones hemos mencionado ya este asunto, la rabiosa actualidad obliga a dedicarle un poco más de espacio y tiempo. Os pongo en antecedentes:
  • El pasado 27 de diciembre se 2017 se publicó el plan Rector de Uso y Gestión (PRUG) del Parque Nacional de la Sierra de Guadarrama para someterlo a información pública y permanecerá expuesto durante el próximo mes y medio para presentar alegaciones.
  • El texto del PRUG regulará la actividad forestal, la ganadería, la practica deportiva de senderismo, escalada y ciclismo.
  • El borrador del PRUG establece severas limitaciones en los recorridos ciclables dentro del Parque Nacional, limitándolos a un listado de caminos recogidos en el Anexo IV.
Así las cosas, el mundo ciclista echa humo y la polémica está servida en foros y comunidades diversas ante la más que probable limitación de uso de la bicicleta dentro del Parque.

Desde Cartografía Digital entendemos ambas posturas y la polémica y dificultades que siempre entraña regular las actividades de un Parque Nacional. Y, siendo conscientes de que la barra libre que imperaba hasta ahora es totalmente incompatible con la debida protección de un espacio natural más masificado de lo habitual dada su situación geográfica, comprendemos también que la comunidad ciclista se ha visto señalada como centro de las prohibiciones. Creemos que ello viene dado fundamentalmente por dos cuestiones:
  • La diferenciación clara con el mundo senderista. Si bien las limitaciones del PRUG también alcanzan a los senderistas tanto en cuanto a zonas prohibidas como en cuanto a número de individuos en los grupos organizados (Art. 45.a), para el mundo ciclista se establecen, además, limitaciones físicas en cuanto a recorridos permitidos (Art. 45.d.1).
  • Esta clara diferenciación entre las dos actividades está en el centro de la polémica. ¿Es más incompatible con la buena conservación del Parque la actividad ciclista que la senderista? ¿En base a qué razones? ¿Erosión de la cubierta terrestre (Art.38)?
La polémica está pues bien servida. El texto completo del Decreto podéis encontrarlo AQUÍ, incluidos los Anexos y cartografía de referencia. Os animamos a leerlo con atención y extraer vuestras propias conclusiones, más allá de las muy interesadas opiniones que cada sector pueda lanzar.

En Cartografía Digital hemos querido colaborar a una mejor comprensión del estado de las cosas, y como lo nuestro son los mapas, lo hemos hecho con ellos. Dado que la claridad de la cartografía contenida en el PRUG puede ser algo difícil de interpretar para el común de los usuarios, hemos elaborado un par de mapas web y los hemos puesto a disposición de grupos ciclistas, para que a su vez la generalidad de la comunidad tenga una visión más amigable de los elementos afectados.

El primero lo encontramos en Alfonso y Amigos, un grupo de aficionados a la bicicleta que han alojado nuestro mapa en su web. Sobre la base de la cartografía topográfica del IGN tenemos superpuesta la cartografía PRUG con todos sus elementos (zonas de protección, recorridos autorizados, zonas de escalada y vivac). A esto le hemos añadido la posibilidad de visualizar sobre ortofoto, y una capa con los recorridos previstos que nos permite obtener información extra de los mismos (nombre del tramo, zona, anchura,...).
PINCHA EN EL MAPA PARA ACCEDER AL VISOR
El segundo es aún más completo y lo tienen alojado los amigos de Getafe Bike en su web. Con la posibilidad de cartografía IGN u ortofoto como mapas base, el mapa cuenta también con las opciones de activación, transparencia o descarga de capas, medición de líneas o áreas sobre el mapa, tablas conteniendo listado de elementos (recorridos, escalada,vivac) para facilitar selección y búsqueda, ventanas pop up con detalles, ... Y una utilidad que nos parece interesante para la preparación por ejemplo de buenas alegaciones: la posibilidad de crear líneas, puntos o áreas en una capa nueva, de forma que se puedan dibujar sobre el terreno aquellos elementos que se considere necesarios para redondear una mejor interpretación del PRUG, y ser así susceptibles de utilización directamente en formato digital.
PINCHA EN EL MAPA PARA ACCEDER AL VISOR
En definitiva, y más allá de opiniones interesadas o ignorantes, una herramienta práctica con la que poder trabajar en lo que verdaderamente importa: generar alegaciones con fundamento e intentar que el resultado final del PRUG deje lo más contentas posible a todas las partes interesadas en la defensa de un bien común. El tiempo se acaba como para perderlo en filosofías.

lunes, 29 de enero de 2018

Caso práctico QGIS: nuevos mapas del Atlas Nacional de España.


Hace unos días dimos cuenta en nuestra página de Facebook de la publicación de un nuevo producto en el Centro de Descargas del CNIG. Englobado bajo la categoría Cartografía para Sistemas de Información Geográfica, el producto lleva por nombre Mapas temáticos del ANE.
En los próximos meses tendrá lugar la publicación de una nueva edición del ANE, denominada España en mapas y perteneciente a la serie Compendios del Atlas Nacional de España (ANE). Todos los mapas de esta publicación estarán incluidos en este nuevo producto y para cada mapa se podrán obtener los datos en formato SHP (siempre que las condiciones de licencia lo permitan), los metadatos en formato XML y PDF, y una versión en JPG y en PDF. Se espera que todos los mapas puedan estar a descarga en el primer trimestre del próximo año 2018.

Por el momento ya están disponibles para descarga los primeros 10 mapas de los más de 800 con que contará la colección.
Podemos descargar de cada uno el JPG, el PDF o un ZIP conteniendo ambos, más los archivos SHP utilizados en su composición (no presentes en todos por desgracia, ya que desde el punto de vista compositivo es lo más interesante). Son todos ellos archivos de poco peso y por tanto tremendamente sencillos de manejar y manipular. Por citar un ejemplo para que os hagáis una idea del contenido, el ZIP del mapa Espana_Elecciones-al-Parlamento-Europeo_2014_mapa_15866_spa incluye los metadatos en formato PDF y XML, el propio mapa en formatos JPG y PDF, y la carpeta llamada datos donde se incluyen archivos útiles para su elaboración (en este caso unos SHP con las provincias de España y archivo Excel con los datos numéricos de las elecciones).
Mapa final y contenido de la carpeta correspondiente.
Hemos comprobado que los PDF de algunos mapas son en realidad extraordinarios GeoPDF, que no sólo están georeferenciados, sino que permiten cargar el mapa desglosado en múltiples capas tanto ráster como vectoriales. Muy buena noticia.
PDF del mapa demarcaciones hidrográficas abierto en formato vectorial y con todo su contenido desglosado por capas.
En Cartografía Digital, aparte de este breve análisis-presentación del producto y de agradecer a nuestro IGN la publicación de una colección tan interesante y prometedora, lo primero que hemos pensado al ver estos mapas es que son un material perfecto para poder hacer uno de nuestros casos prácticos. Para ello vamos a utilizar el material del mapa Espana_Usos-de-agua-reutilizada_2013_mapa_15172_spa y vamos a trabajarlo en QGIS para conseguir un producto de aspecto similar al mapa final presentado (mapa con representación en diagramas).

Si os interesa aprender un poco más del manejo de QGIS, os animo a descargar también vosotros los archivos e ir siguiendo nuestro artículo como ejercicio. Comenzamos!!




Cargamos en QGIS los archivos SHP y la hoja Excel disponibles (para este caso práctico y a efectos de visibilidad en pantalla sólo hemos utilizado el SHP de la Península, descartando el SHP de Canarias):
Para una mayor comodidad posterior renombramos las capas con nombres identificativos que nos ayuden, en este caso con el nombre correspondiente a cada variable de cada hoja y que podemos encontrar en su propia tabla de atributos o en el mismo archivo Excel:
Ejemplo de capas renombradas para mayor comodidad con el contenido de cada una.
El siguiente paso es realizar una unión de datos de las diferentes tablas Excel sobre la capa vectorial SHP cuyo procedimiento ya explicamos en este artículo. Para ello necesitamos saber el campo común que todas las tablas comparten y el campo de datos que queremos añadir a la tabla común. Todo esto lo descubrimos en las tablas de atributos:
Comprobamos que la columna con valor común que comparten las tablas es id_ccaa en la tabla SHP y Código unidad territorial en las tablas Excel.
Procedemos a hacer la unión de datos sobre la capa CCAA; para ello entramos en las Propiedades de la capa y en la pestaña Uniones. Pulsamos el símbolo + para crear una nueva unión:
En la ventana de unión seleccionamos los datos para una correcta unión: la capa que vamos a unir (en este caso nuestra primera hoja Excel Agricultura Variable_8713), el campo de unión y el campo objetivo de ambas tablas (ya hemos visto más arriba cuáles eran) y el campo que queremos unir (como podemos ver en el Excel o las tablas el dato numérico que nos interesa representar es la columna Valor del dato). Podemos añadir también un prefijo personalizado para identificar mejor las columnas pero, al haber renombrado previamente las capas, el que sale por defecto nos puede servir.
Cuando todo esté listo pulsamos Aceptar y nuestra primera unión queda creada. Ahora repetimos el proceso con el resto de capas Excel, pues nuestro objetivo es añadir los datos de todas las hojas Excel al vectorial SHP. Una vez hecho todas las uniones se muestran en la ventana, desde donde podremos eliminar o modificar la que nos interese:
Terminadas de hacer todas las uniones que necesitamos, aplicamos y aceptamos.
Es momento ahora de comprobar que nuestras uniones han sido bien hechas y la tabla de atributos de nuestra capa vectorial SHP contiene las nuevas columnas con los nuevos datos unidos:
Nótese que Ceuta y Melilla aparecen como NULL al no coincidir el campo de unión con las tablas Excel.
Perfecto. Ya disponemos de la base correcta para comenzar la representación de la información.
Como hemos visto en el mapa que pretendemos reproducir, los datos están presentados en forma de típicos diagramas de queso, en los que cada porción representa un porcentaje del total. En las Propiedades de nuestra capa CCAA pinchamos la pestaña Diagramas:
1. En el selector de diagramas elegimos "Gráfico de queso". 2. Comenzamos configurando la primera sección Atributos. 3. Seleccionamos los atributos a representar y los añadimos a la ventana Atributos asignados. Obviamos el valor 8712 ya que es el valor del agua total reutilizada y como vemos en el Mapa ANE se utiliza para determinar la proporción por tamaños de cada queso. 4. Haciendo doble click en cada campo podemos seleccionar un color y una leyenda para cada atributo. También podemos modificar el orden de representación en la esfera horaria del quesito arrastrando cada atributo para modificar su posición en la lista.
Podemos usar el botón Aplicar para ir viendo sobre el lienzo los cambios mientras seguimos configurando secciones:
En la sección Apariencia podemos 1. modificar la transparencia global de los quesitos, su color y anchura de línea y el ángulo de inicio (a partir del cual comienza la representación de las porciones según el sentido horario). 2. Decidir si la visibilidad de los quesos depende de algún atributo o de la escala de visualización (escogemos el valor 8712 para que aquellas CCA sin datos globales no representen su quesito). Pulsamos de nuevo Aplicar para ver los cambios.
A estas alturas el gráfico debe de tener el aspecto de la siguiente imagen (nótese que Ceuta y Melilla representan quesito vacío debido a que las hojas Excel muestran ambas ciudades unidas y el la tabla SHP figuran independientes, por lo que la unión de datos del principio no ha sido correcta y deberíamos modificarlo manualmente):
Continuamos configurando la sección Tamaño:
1. Nuestra primera elección es si dar a todos los quesos un tamaño fijo o un tamaño escalado. 2. Como queremos que el tamaño de los quesos represente la proporción total de agua reutilizada seleccionamos dicho valor en la pestaña Atributo. Establecemos el valor máximo para el diámetro más grande a representar (el botón Encontrar lo completa automáticamente), y un tamaño escalado para controlar desproporciones. Por último, y para evitar esas grandes desproporciones, podemos establecer un tamaño mínimo para los quesos más pequeños a fin de que se visualicen con más facilidad. Pulsamos Aplicar para comprobar los cambios.
Elegimos como Ubicación de los quesitos el centroide de cada CCAA.
Aplicamos y aceptamos todos los cambios y ya tendríamos la base del mapa, con sus datos representados y la simbología básica aplicada. Un poco de trabajo en el diseñador de impresión de QGIS para rematar el mapa y nuestro caso práctico estaría terminado. Saludos!