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.