lunes, 30 de mayo de 2016

Mapa del mes: EUSKADI relieve PACK.

Este mayo que acaba hemos preparado un mapa del mes distinto: un pack completo de Euskadi bajo el tema común del relieve. En dicho lote hemos incluido los principales productos que solemos usar en cartografía cuando trabajamos con relieves de una zona, y los hemos agrupado en un archivo RAR.

Su contenido es el siguiente (148 megas):


- ASTER_euskadi_38m.asc: modelo de elevaciones de Euskadi en formato asc con paso de malla 38,2 metros.

- Euskadi_areas.shp: polígonos vectoriales representando las alturas de Euskadi en rangos de 20 metros..

- Euskadi_curvas.shp: lineas vectoriales representando las curvas de nivel de Euskadi con rango de 20 metros.

- Relieve_color.tif: imagen tif representando el relieve de Euskadi con tintas hipsográficas por alturas.

- Sombras.tif: archivo tif conteniendo el relieve sombreado (hillshade) de Euskadi.

Vamos a mostrar y comentar un poco más en detalle cada uno de los productos contenidos en el lote.

ASTER_euskadi_38m.asc
A partir de los datos de relieve de ASTER GDEM v2 abiertos en Global Mapper desde su conexión online, hemos encuadrado la zona de Euskadi y exportado en formato asc a la mayor resolución que permite este proceso (38,2 metros). El área es cuadrada y cubre la totalidad de Euskadi y zonas limítrofes, pues es una faena exportar justo por los límites de la Comunidad para los usuarios cercanos a las fronteras. Reproyectado a UTM zona 30 ETRS89.
Archivo de relieve con diferentes representaciones por defecto según el programa: a la izquierda abierto en QGIS en escala de grises sin sombreado; a la derecha abierto en Global Mapper con sombreado y coloreado de alturas.
Euskadi_areas.shp
Desde el mismo Global Mapper y a partir del relieve anterior hemos extraído los polígonos por alturas del relieve con un rango de 20 metros. En realidad se trata de un proceso similar a extraer las curvas de nivel, pero generando polígonos cerrados en vez de lineas, lo que nos da muchas posibilidades a la hora de trabajar con ellos.
Detalle de zona con los polígonos por alturas.
Euskadi_curvas.shp
Prácticamente igual que el anterior, esta vez se trata de las curvas de nivel de Euskadi, extraídas cada 20 metros. Procesos ambos muy sencillos, que Global Mapper ejecuta con enorme fluidez como ya hemos explicado en algún artículo. Tanto las áreas como las curvas están sacadas en UTM zona 30 ETRS89.
Vista general del archivo de curvas de nivel.
Relieve_color.tif
Abriendo el archivo de relieve asc en QGIS hemos generado un relieve coloreado con tintas hipsográficas representando las diferentes altitudes de la comunidad. Empleando la gama de colores standard con variaciones cada 200 metros, el resultado es el siguiente:
Creemos que la imagen tif puede dar juego como fondo de cartografía, o como base para diferentes productos.

Sombras.tif
También desde QGIS hemos realizado el clásico proceso hillshade, o relieve de sombras, sobre el archivo de relieve asc original. Con los parámetros por defecto del programa, este archivo es imprescindible a la hora de aplicar sombreados sobre otras cartografías, dándoles esa sensación de volumen que las distingue de los mapas planos.
Obviamente, la resoluciones del archivo de color y el de sombras son la misma del asc original (38,2 metros), apropiadas para trabajos a escalas intermedias. También están explicados en algún artículo del blog estos procesos.

Esperamos que este lote os pueda ser de alguna utilidad al trabajar sobre la zona de Euskadi. Lo tenéis disponible para descarga en la sección Mapa del Mes. En la siguiente imagen se aprecia la generosa cobertura de terreno que abarca este pack de relieve respecto a la comunidad. Saludos!!
Área cubierta por los archivos respecto al territorio de la Comunidad de Euskadi.

jueves, 26 de mayo de 2016

SASPlanet: novedades mayo.

Ya nos va costando encontrar cosas interesantes, pero ahí vamos y nos ponemos desde ya a buscar las siguientes para la próxima. Las novedades en mapas de este mes de mayo para nuestros lotes personalizados de SASPlanet disponibles en la sección DESCARGAS son las siguientes:

CATALUNYA

- TOPOtotal 2. Nueva versión del conjunto de mapas topográficos de Catalunya que comprende todas las escalas y que están mejor escalados y se visualizan mejor que en el TOPOtotal que había hasta ahora. Captura a zoom 19:
- ORTOtotal. Del mismo WMS que el anterior, multiescala continuo con ortofotos de Catalunya. Captura a zoom 19.

SPAIN

- Montes. Servicio WMS del MAGRAMA con la propiedad de los Montes Públicos. Si os es más útil como CAPA cambiad la configuración a asLayer=1. Con Leyenda. Captura a zoom 15.


GALICIA

- Cartas Náuticas de las costas gallegas provenientes del visor de Intecmar. Captura zoom 15. Thanks Pascal!




Para el lote SASPlanetINTERNATIONAL, las novedades son:

CHILE

- IGM50K. Mapas del Instituto Geográfico Militar a escala 1:50000. Abarcan todo el país. Captura a zoom 15.
- Raster50K. Topográfico a escala 50k de todo Chile. Similar al anterior. Captura a zoom 16.


COLOMBIA


- BaseRelive. Mapa base con relieve de Colombia proveniente del IGAC. Captura a zoom 14.
- Imagenes. Mosaico de imágenes del IGAC, contiene tanto Landsat como ortofotos 1:2000 en algunas zonas. Carga algo lenta y con errores. Captura a zoom 14 donde se aprecia área de imágenes Landsat y área de ortofotos.


CAMERÚN

- GeoCameroun. Topográfico de sorprendente nivel para ser zona africana. Captura a zoom 16.
- GeoCameroun OmbrageCAPA con sombreado de relieve de Camerún. Ideal para combinar con el mapa anterior. Captura a zoom 14. Thanks Pascal.

TURKEY

- PergoRoads. De la empresa Pergo nos llegan estos mapas de Turquía, bastante detallados y con estilo OSM. Thanks jeanbon. Captura a zoom 14.


lunes, 23 de mayo de 2016

INE: fuente inagotable de datos.

A menudo repetimos que para elaborar una buena cartografía temática, para que nuestros mapas reflejen fielmente realidades, es imprescindible contar con los datos adecuados. Y también a menudo decimos que conseguir dichos datos no es tarea fácil entre la maraña de información que nos rodea, muchas veces caótica y escasa de criterio.
Por suerte siempre tendremos el INE. 
El Instituto Nacional de Estadística (INE) es como un oasis en medio del desierto, o mejor aún, como un remanso en medio de la jungla. En sus páginas se esconden los datos más insospechados, las tablas más inesperadas; pero también los principales indicadores de un país y todos aquellos datos básicos para una buena cartografía por temas. Sin duda un indispensable a tener siempre en cuenta cuando buscamos información que cartografiar.
Sería absurdo hacer un resumen de la información disponible en el INE, pues es tanta y tan variada que nunca daríamos una idea de lo que allí se puede encontrar. La información se estructura en algunas grandes categorías que nos orientan el rumbo, pero os recomiendo que paséis un rato navegando por sus páginas para haceros una idea aproximada de todos los datos a nuestro alcance.
Portada web del INE con las principales áreas de búsqueda de datos.
Una vez que hayamos alcanzado la tabla de datos de nuestro interés, el INE nos ofrece la posibilidad de descargarla en varios formatos habituales: Excel (extensiones XLS y XLSx), CSV (separado por tabuladores, coma o punto y coma) y Pc-Axis (formato propietario del INE para el que nos ofrecen un software gratuito que no he probado, pero que tiene buena pinta y parece que permite elaborar cartografía de manera sencilla combinando los datos con vectoriales).
Ejemplo de tabla INE (ocupación de alojamientos rurales 2014 por provincias y meses) y acceso al botón descarga
con la ventana de elección de los diferentes formatos disponibles.
Descargada la tabla, podremos aplicar en QGIS su importación (como vimos en este artículo) y/o su unión con vectoriales que los representen (como vimos en este otro artículo). He probado a descargar la misma tabla en formatos xls, xlsx y csv de punto y coma, y con QGIS 2.14 Essen es posible abrirlos todos (los excel como capa vectorial y los csv importando capa de texto delimitado). Eso si, difícilmente podremos evitar una previa manipulación de la tabla para adaptar correctamente su estructura y poder tenerla en QGIS perfectamente utilizable, pues es normal que contengan distintas cabeceras que debemos suprimir o que debamos retocar los nombres de la columna en el campo común de unión.
Tabla CSV del ejemplo anterior con los parámetros ajustados para una correcta importación.
Otro ejemplo con el CSV de la producción editorial por CCAA y temas.
Y no se nos ocurre mejor forma de terminar este breve artículo que haciendo unos mapas con la cantidad de tocayos que tengo repartidos por España: nada menos que 303.247 Javieres. Y con una media de edad de 30,9 años.
Quién la pillara.... Saludos a todos!!

miércoles, 18 de mayo de 2016

Fotogramas de vuelos, nueva descarga CNIG.

El IGN ha añadido un nuevo producto a la lista del Centro de Descargas: fotogramas de vuelos. Se trata de los fotogramas pertenecientes a los conocidos vuelos históricos 1929-1930 Ruiz de Alda, 1973-1986 Interministerial, 1980-1986 Nacional, 1989-1991 Costas, y los vuelos PNOA de 2012, 2013 y 2014.
Y decimos vuelos conocidos porque todos ellos ya se encontraban en la Fototeca Digital, desde donde tenemos la posibilidad de visualizar e imprimir estas imágenes, incluidas las de las dos series de vuelos americanos (1945-1946 y 1956-1957) que de momento no parecen estar presentes en este nuevo producto.
En este caso se trata de la descarga en formato ECW de los fotogramas individualizados, de diferentes años, escalas y sistemas de referencia, resultado de escanear el negativo analógico del fotograma, y que tampoco están escaneados en su totalidad ni cubren la totalidad del territorio. Los fotogramas correspondientes a los vuelos PNOA son digitales a partir de 2007, y cubren las diferentes zonas (a diferentes resoluciones) según el plan de vuelos anual que puede consultarse en este mapa.
Ejemplo de fotograma de año 1983.
Los fotogramas están georeferenciados, aunque en diferentes proyecciones y de una forma aproximada, con lo que no debemos esperar en ningún caso la precisión de las ortofotos PNOA, y debemos tomarlos más como imágenes históricas que como un producto de exactitud fotogramétrica. En una prueba con los fotogramas del Vuelo Interministerial encontramos por ejemplo el siguiente desfase en la calibración:
Los fotogramas nos aparecen desplazados unos 260 metros dirección Sureste respecto al raster IGN.
Estas desviaciones y las características de las imágenes de cada vuelo están perfectamente explicadas en el documento pdf de la información auxiliar del producto.
Reseñar que la nomenclatura de las imágenes es diferente para cada vuelo, cosa muy lógica teniendo en cuenta que los fotogramas tienen diferentes extensiones y resoluciones. Por tanto, y dada la enorme cantidad de imágenes disponibles, nos será tremendamente complicado adivinar qué fotograma concreto cubre la zona concreta que nos interese, y aunque refinemos mucho la búsqueda, normalmente nos veremos obligados a descargar bastantes imágenes para acertar el tiro.
Diferentes tamaños, resoluciones y terreno cubierto por cada fotograma dependiendo del vuelo al que pertenezca. De izda. a dcha: PNOA2014, Vuelo Interministerial y Vuelo Nacional.
Las características de las imágenes de los vuelos históricos (excluido PNOA) son las siguientes:

Enhorabuena una vez más a nuestro IGN por seguir aportando a la biblioteca pública común material de enorme valor y calidad. Un saludo.

jueves, 12 de mayo de 2016

QGIS: uniendo tablas de datos.

A menudo queremos representar en un mapa algún archivo vectorial pero nos encontramos con que su tabla de atributos es escasa en datos o no contiene aquellos que necesitamos. Y también a menudo sucede que esa información que necesitamos la tenemos en una base de datos aparte. Es entonces cuando debemos fusionar ambos archivos de alguna forma que nos permita disponer de toda la información en aquel archivo con el que deseamos trabajar. En este caso práctico vamos a ver lo que para mi es uno de los procesos más básicos e imprescindibles a la hora de crear nuestros mapas: enriquecer la información de un archivo vectorial con una tabla de datos externa, lo que nos ampliará las posibilidades de representación.
Vamos a utilizar para ello los polígonos de los municipios de España en formato vectorial shp (incluido en la descarga Líneas Límite Municipales) y la base de datos del Nomenclátor Geográfico de Municipios y Entidades de Población en formato mdb Access, ambos disponibles en la sección Equipamiento Geográfico de Referencia Nacional del Centro de Descargas del IGN.

Una vez descargados ambos archivos, los abrimos en QGIS. La base de datos mdb, al constar de varias tablas, nos preguntará cuál de ellas queremos abrir en QGIS: con abrir la tabla municipios nos basta para este caso.
El vectorial shp de polígonos municipales que hemos utilizado es el correspondiente al datum ETRS89.
Ambos abiertos en QGIS y consultando sus correspondientes tablas de atributos tenemos lo siguiente:
Apreciamos rápidamente que ambas tablas comparten como columna de valores comunes el nombre del municipio (llamado NOMBRE en la tabla del vectorial, y NOMBRE_ACTUAL en la tabla mdb). Podríamos haber tenido suerte y que el campo del Código INE fuera también idéntico, pero inexplicablemente está formado por distinto número de dígitos en cada tabla con lo que no es utilizable.
A continuación debemos unir ambas tablas. Para ello abrimos las propiedades de la capa vectorial (que es a la que queremos añadir los campos) y en la pestaña Uniones creamos una nueva unión según se ve en la siguiente imagen:
Creamos una nueva unión desde el símbolo +. En Unir capa obviamente seleccionamos la tabla mdb; en Unir campo y Campo objetivo seleccionamos el campo común de cada una de las dos tablas. Podemos elegir si unir la tabla entera o sólo algunas de las columnas (zona verde) y también el prefijo que distinguirá a las nuevas columnas añadidas (en azul). Pulsamos en Aceptar.
La nueva unión queda añadida a la lista de uniones mostrando sus parámetros básicos. Aplicamos y aceptamos.
Ahora nuestra capa vectorial de polígonos cuenta en su tabla de atributos, además de con sus datos originales, con los nuevos campos que le hemos añadido desde la tabla de datos mdb.
Columnas o campos nuevos añadidos (en amarillo) con el nombre de campo predefinido en el prefijo (verde).
Por lo que, y esto es lo importante, cuando vayamos a modificar el estilo de la capa (graduado, categorizado, etc...) las nuevas columnas de datos están ya disponibles para poder utilizarse en la simbología de capa.
La unión de datos hasta este momento es sólo virtual y se pierde si cerramos alguna tabla, por lo que conviene guardar la capa vectorial como un nuevo archivo shp que mantenga la nueva estructura de la tabla de datos de forma permanente. Simplemente botón derecho sobre la capa y Guardar como...

¿Qué hemos conseguido con esto? Que donde había una tabla de datos con escasa información que apenas podíamos clasificar por ningún dato relevante, tenemos ahora algunos de los datos más importantes para cualquier simbología que queramos aplicar a los municipios: población, superficie, perímetro, altitud...

Jugando ahora con nuestra nueva tabla de datos integrada en el vectorial y el magnífico QGIS nos será muy fácil elaborar interesantes mapas sobre la base de los municipios españoles. Para muestra tres ejemplos:


El origen y formato de la tabla de datos que queremos unir al archivo vectorial puede ser diverso, pero sirva este artículo de idea global sobre las posibilidades de este proceso. Un saludo.

lunes, 9 de mayo de 2016

CORINE Land Cover.

Si el mes pasado hablábamos del SIOSE, el producto estrella en cuanto a ocupación del suelo del IGN, no podemos dejar de rematar la faena hablando de CORINE Land Cover (CLC), el otro producto destacable dentro de esta temática en nuestro imprescindible Centro de Descargas.
Lo primero que se nos ocurre pensar es ¿qué diferencias tiene con el SIOSE? Podéis echar un vistazo a la web del IGN sobre el asunto, o también a la web de la EEA para una completa información. Nosotros, en nuestro estilo de andar por casa, resumiríamos las principales diferencias así:
  • La escala de referencia de CORINE es 1:100000 frente al 1:25000 del SIOSE.
  • La unidad mínima cartogafiable por defecto son 25 hectáreas (0,5 a 2 hectáreas en SIOSE).
  • La clasificación del terreno y sus niveles de jerarquía son más simplificados que en SIOSE.
Digamos que, básicamente, CORINE Land Cover es un proyecto a nivel europeo que trabaja la información a mayores escalas, por lo que es un producto más simplificado tanto en polígonos como en categorías de datos. Como solemos hacer, una vez ofrecidos los enlaces para interesados, vamos a lo práctico. Un primer vistazo a las posibles descargas del producto nos arroja estos resultados:
Junto a los archivos de metadatos de cada serie (naranja), lo cierto es que aparecen varias descargas de varios productos y actualizaciones (amarillo) y el producto de 2006 (verde) que da la impresión de ser el único unificado e ideal para descargar.
Las descargas en este producto, a diferencia del SIOSE que son por Comunidades Autónomas, son a nivel nacional. Descargado y descomprimido el archivo ZIP correspondiente a 2006 este es su contenido:


- Carpeta con simbología (códigos de colores RGB empleados en la clasificación en formato mdb Access).

- Vectoriales con los polígonos para las Islas Canarias (amarillo).



- Vectoriales con los polígonos para toda España. Canarias incluido (verde).


- Diversos archivos TXT y PDF con información del producto y su metodología, que también contienen los imprescindibles códigos numéricos de clasificación del terreno (Anexo del pdf).

La verdad que este tipo de productos tienen un abordaje complicado, casi siempre debido a la codificación numérica de los diferentes tipos de coberturas. Abrimos el vectorial de toda España -CLC06_es.shp- en QGIS para hacernos la siempre recomendable visión global (542 megas, no es poca cosa):
CORINE Land Cover abierto en QGIS con detalle de zona.
Un vistazo siempre igual de recomendable a la tabla de atributos del archivo vectorial nos proporciona información muy importante; en este caso los campos que contiene dicha tabla son los siguientes:
155770 polígonos en la tabla. Area_ha: extensión en hectáreas del polígono. Shape_leng: longitud del vectorial del polígono (equivalente a su perímetro en metros). Shape_Area: extensión del polígono, esta vez en metros cuadrados. code_06: código identificativo de la ocupación del suelo. ID: identificador único de cada polígono.
Aunque todos los campos tienen su importancia, obviamente el campo code_06 es vital. En él se contiene el código que clasifica la cobertura del suelo para dicho polígono y que podemos encontramos en el anexo del PDF de la descarga (44 clases en total). Clasificamos y categorizamos una zona al azar para verlo mejor:
La utilidad de CORINE Land Cover es evidente, sobre todo a escalas intermedias, y su facilidad de uso en comparación al SIOSE toda una ventaja. Así por ejemplo, en pocos minutos, sobre el sombreado de relieve del Mapa Base del IGN, activamos las coberturas que nos interesen y conseguimos un interesante mapa de ejemplo con abundante y utílisima información:
Coberturas con su código CORINE sobre relieve de una zona de Huesca.
No olvidemos para terminar que este servicio está disponible vía WMS (también WFS), y en su correspondiente capa con su correspondiente leyenda dependiente de la escala de visualización, en la misma URL que el SIOSE: http://www.ign.es/wms-inspire/ocupacion-suelo? y que a través de una petición GetFeatureInfo podemos acceder a los datos básicos de cada polígono. También contiene capas con las variaciones históricas de la cobertura de los polígonos en cada actualización CORINE.
Consulta en Global Mapper de información sobre polígono desde el servicio WMS de CORINE.
Si os interesan los aspectos más técnicos sobre el futuro de los productos relacionados con la ocupación del suelo, se han publicado guías al respecto (final de página bajo el epígrafe Implementación española) de cómo los datos deben ajustarse a la directiva INSPIRE para fomentar la máxima interoperabilidad entre organismos. Saludos y a disfrutar de este extraordinario producto, para mi básico en la elaboración de cartografía topográfica, pues nos permite aportar un útil e interesante contenido extra.

jueves, 5 de mayo de 2016

Nuevos productos IGN: lineas de costa.

Nuestro IGN sigue haciendo los deberes y publicando material cartográfico con una buena frecuencia. En este caso, y en colaboración con el Instituto Hidrográfico de la Marina, pone a nuestra disposición tres nuevos productos en la sección Equipamiento Geográfico de Referencia Nacional.
Se trata de la Linea de Costa, las Lineas de Base Rectas y la Zona Económica Exclusiva en el Mediterráneo. Todos ellos archivos vectoriales en formato shape SHP de los que tenemos cumplida información y detalles en el interior de la correspondiente descarga. Son archivos a priori muy sencillos, pero "cuyo fin es proporcionar una infraestructura básica de datos geográficos al sector SIG español. Destinado a todo tipo de usuarios (empresas, Administración), para todo tipo de aplicaciones y sobre todo tipo de plataformas (sistemas vectoriales, ráster, orientados a objetos, etcétera)."
Hemos representado con QGIS los tres productos en un mapa rápido y sencillo, para que tengáis una idea general de los lineas de que se trata. Sobre la ortofoto de la isla de Ibiza, este sería su contenido:
Hemos dejado los nombres originales de los vectoriales, de tal forma que SLCONS serían las líneas de las construcciones artificiales, DEPCNT la linea de bajamar, COALNE la linea de costa natural y lbr las lineas de base rectas.
Cada vectorial consta de su tabla de datos o atributos, más simple o más completa, con los datos necesarios.
Tabla de atributos del vectorial Lineas de Base Rectas.
Y este sería el aspecto general, con la linea y puntos de la zona económica exclusiva del Mediterráneo incluida.
En definitiva, unos datos de referencia que completan aún más el equipamiento geográfico nacional y que cumplen su función de infraestructura básica espacial para trabajos SIG. Se agradece que nuestro IGN los ofrezca en un repositorio centralizado, evitando mucha búsqueda inútil y pérdidas de tiempo. Un saludo.