BTN100: la base topográfica de referencia.

Con motivo de la actualización este mes de julio de 2015 de la información contenida en la BTN100 hemos pensado que es un buen momento para profundizar un poco más en ella y desgranar tanto su contenido como la utilidad e importancia que tiene en el entramado cartográfico nacional. Con esta actualización el IGN además asegura que "con el nuevo protocolo de actualización continuo de BTN100 podrás tener cada tres meses todas las novedades: nuevas inauguraciones de infraestructuras, cambios en la denominación de municipios y poblaciones, además de todas las incidencias que nos vayáis reportando." Sin duda un motivo más para congratularse.

Como vemos en la captura de arriba, en la información básica de su cabecera en el Centro de Descargas, se nos dan ya varios datos fundamentales:
1.- El formato en el que se ofrece es shapefile, es decir archivos vectoriales .shp.
2.- Contiene 56 capas de información geográfica, con datos topográficos y temáticos. Os pongo captura del archivo shape de cada capa en la que se ve su nombre y el tamaño del mismo (aunque esto es un poco indiferente pues los shp ya sabéis que van acompañados de archivos adjuntos .dbf y .prj que son sus bases de datos y localización).

3.- Podemos descargar un archivo continuo para toda España (con todas las capas, que está por encima de 700 megas) o archivos independientes agrupados por temáticas.

Descarga de la BTN100 en la que se observan las 7 diferentes temáticas así como su tamaño. El zip llamado BTN100 sería el archivo completo con todas las temáticas.

4.- El sistema geodésico de referencia es ETRS89 con coordenadas geográficas de latitud y longitud. La escala continua es de 1:100000 (resolución de 20 metros).

En el enlace Información auxiliar BTN100 podemos encontrar varios documentos en PDF con información técnica del producto (fechas, especificaciones, catálogo de objetos,...) así como el archivo xml de los metadatos. Así mismo en este enlace del IGN tenemos algo más de información sobre la BTN100.
Ya que la BTN100 está diseñada para su explotación mediante Sistemas de Información Geográfica (SIG) y la producción de productos cartográficos, abramos alguna de sus capas shape en QGIS para ver qué aspecto tienen y qué información contienen.

Capa de polígonos de los Núcleos de Población abierta en QGIS. En su tabla de atributos, con 26913 elementos, vemos la información de la que se le ha dotado (nombre y población como datos más relevantes o útiles).

Otro ejemplo de la tabla de atributos de la capa de puntos de los vértices geodésicos (11141 elementos) mostrando sus columnas de información tremendamente útiles: nombre, coordenadas, altitud,...

La BTN100 nos proporciona de esta forma toda la información básica y los elementos necesarios para elaborar cartografía a escalas intermedias (para escalas más detalladas contamos con la BTN25). Obviamente los archivos vectoriales abiertos en bruto en un programa ofrecerán el aspecto de un conglomerado de líneas, puntos y polígonos sin ninguna simbología, estilo ni etiquetado.

BTN100 abierta en bruto en Global Mapper. El Control Center nos muestra las 56 capas de que consta y el número de elementos de cada una. Captura sobre la ciudad de Madrid a escala 1:100000.

Pero con un manejo básico de QGIS y en pocos minutos podremos tener cartografiada una zona, simplemente añadiendo los elementos principales del territorio y dotándoles de una simbología, estilo y etiquetado básicos. Ya luego es cuestión de las necesidades y gustos de cada caso añadir los elementos necesarios, buscar la escala idónea, añadir algún fondo, etc, etc... hasta conseguir el mapa que buscamos.

Captura a escala 1:100000 con simplemente algunas capas principales de la BTN100 modificadas.

De gran importancia es también la cantidad y calidad de los datos que conformen las tablas de atributos de las diferentes capas, pues ello nos permitirá hacer análisis de datos geográficos con mayor facilidad utilizando algoritmos y herramientas de procesado que todo buen SIG contiene (y que yo apenas se manejar). Pero teniendo los datos y con simple simbología de capas podemos hacer en minutos sencillos estudios demográficos...

...curiosas representaciones de los vértices geodésicos por altura...

...un mapa de calor con la capa de puntos Aeropuertos...

...o un buffer de 100 kms en la capa de puntos Estaciones de GPS para ver su solapamiento.

La BTN100 es pues una base de datos geográfica actualizada que resulta de obligada referencia a la hora de utilizar recursos geográficos a las escalas adecuadas. Nuestra imaginación, conocimientos y necesidades harán el resto a la hora de configurar el análisis o resultado cartográfico que deseemos.

SASPlanetINTERNATIONAL: nueva actualización.

Pasamos a enumerar las últimas incorporaciones de mapas de esta nueva revisión de nuestro lote SASPlanetINTERNATIONAL de mapas del mundo disponible en la sección DESCARGAS. Como sabéis este lote contiene mapas de países que no son España, ni sus vecinos Andorra, Portugal y Francia, los cuales están incluidos en la descarga SASPlanet.Release.141212+Mapas3.0

GUATEMALA

Gracias al amigo del blog Neko que nos pasó la dirección de una IDE de Guatemala hemos cubierto la cartografía principal para SASPlanet de este país. Hay multitud de capas más, como podéis ver en su archivo de capacidades (por si queremos añadir alguna otra que nos interese):

http://ide.segeplan.gob.gt/geoserver/ows?service=wms&request=getcapabilities

- Hojas Cartográficas: mapa topográfico tipo escaneado del año 2010. Captura a zoom 16.

- Mapa BASE: bonito mapa topográfico tipo vectorial. Captura a zoom 17.

- ORTOFOTO: ortofoto de todo el país. Comienza con imágenes Landsat y a partir del zoom 14 se visualiza la ortofoto en detalle.

- Carreteras/Veredas: muestra carreteras y veredas. Añadida como CAPA pues tiene mayor utilidad.

Ortofoto de Guatemala con capa de carreteras y veredas superpuesta. Zoom 18.

AUSTRALIA

- Topográfico: Precioso y completísimo topográfico de Australia. Una pena que desaparezca a zoom 14.

- Mapa BASE: mapa base topográfico de Australia hasta a nivel callejero. Captura a zoom 16.

- TOPO 250K: topográfico a escala 250k. Captura a zoom 13 (a mayor zoom se pìxela ya demasiado).

CANADÁ

- Toporama: Interesante mapa con la topografía básica de Canadá. Desaparece a zoom 20. Ejemplo zoom 16.

MARRUECOS

- AGADIRtopo: Dentro de la escasez de cartografía marroquí encontramos esta verdadera joya topográfica. Lástima que sólo cubra la zona de Agadir.

IRLANDA

- Mapa Base: Completo multimapa que comienza con la vista del relieve del terreno, a zoom 13 se superpone el topográfico base y a zoom 17 entra la ortofoto (2005). Captura a zoom 15.

ISLANDIA

Gracias a nuestro genial amigo jeanbon del foro en francés ttMaps incluimos dos mapas de Islandia.

- IsAtlasVect: es una bonita mezcla de topográfico y vectorial. Ejemplo a zoom 14.

- Isatlas100r: es el clásico mapa topográfico estilo escaneado que se pixela a zooms muy altos. Ejemplo zoom 14.

SUIZA

Gracias a Pascal, amigo del blog, tenemos dos mapas de Suiza (los anteriores habían dejado de funcionar). Ambos mapas incluyen también el Principado de Liechtenstein.

- Swisstopo-swissttlm3D: precioso mapa topográfico de gran calidad con datos vectoriales sobre relieve sombreado y sensación 3D. Con una alta precisión e integración de la tercera dimensión, swissTLM3D es el conjunto de datos vectoriales 3D más preciso y más completo de Suiza. Ejemplo a zoom 17.

- Swisstopo-Cartes: mapa multiescala topográfico básico de Suiza en el que el mapa va ganando detalle a mayor zoom. Ejemplo a zoom 17.

De momento estos son los mapas que hemos podido añadir al lote INTERNATIONAL para SASPlanet. Ya sabéis que si conocéis algún servicio wms que os interese podéis comentarlo e intentamos añadirlo a la próxima revisión de mapas. Saludos.

SASPlanet: nueva revisión.

Nueva revisión de nuestro lote SASPlanet disponible en la sección DESCARGAS. Estas son las novedades o modificaciones efectuadas en esta ocasión en nuestro repaso por la magnífica cartografía de esta España nuestra. Saludos.

CANARIAS

- Actualizado Topo Integrado 2010 por el de 2011. Los hay de 2012 y 2013 pero de islas distintas en servicios wms distintos.

- Modificada URL de la orto express 2013.
- Modificada URL del Hipsométrico (antes elevaciones).
- Añadido precioso modelo del terreno con datos LIDAR de 2011-2012. Resolución: 2,5 m/pixel.

.

SPAIN

- Añadidas dos nuevas capas del PNOA histórico (2009 y 2010).

- Modificada url de 1ª edición MTN50. Ahora consta de dos capas: una muestra el numero de hoja y es enlace a las distintas ediciones del centro de descargas (esta función no está disponible en SASPlanet, si en otros visores de servicios wms). La otra es el mosaico con las hojas de la primera edición MTN50.

El número de hoja no molesta demasiado, pero si queremos extraer la imagen limpia del mapa bastará con eliminar la capa "Enlace" en nuestro archivo params.txt correspondiente a este mapa.

Visualización en Global Mapper con la posibilidad de enlace directo al Centro de Descargas
 para descarga de las versiones disponibles de cada hoja MTN50.

- Actualizado SIOSE al año 2011. Dos niveles de zoom. Visible a partir de zoom 12; con Leyenda nueva.

- Revisada CAPA de CARTOCIUDAD con nuevos parámetros.

CASTILLA-LA MANCHA

- Añadidas CAPAS de carreteras y de nombres (ésta sólo actualizada y más completa en una zona concreta de la Comunidad).

Capas de nombres y carreteras sobre fondo PNOA. Zoom 15.

- Añadido Mapa Base. En realidad es un magnífico mapa base de ESRI que cubre toda España, pero es la única forma que hemos encontrado de tener alguna cartografía de Castilla-La Mancha,  pues no conocemos que tenga ni un solo servicio wms público. Lamentable y un caso único en España.

ANDALUCIA

- Actualizada la Ortofoto 2010-11 por la de 2013 (0,5 metros/pixel).
- DEA100. Mapa de los Datos Espaciales de Andalucía para escalas intermedias.

- Toporaster10: Sustituye al MTA10 ya que integra toda la topografía armonizada de Andalucía. Lo hemos configurado como CAPA por las razones que explico abajo.

- Sombras SIOSE: es un sombreado orográfico con los usos del suelo de Andalucía. Perfecto para fondos de mapas; de hecho, el producto final Toporaster10 está formado por la combinación del vectorial más este fondo sombreado. Desaparece a zoom 18.

Existe un servicio WMTS que ofrece el Toporaster10 completo (vectoriales más fondo sombreado); su dirección es http://www.ideandalucia.es/geowebcache/service/wmts? Yo no estoy seguro si SASPlanet puede visualizar servicios WMTS (creo que si), en cualquier caso desconozco cómo configurarlos. Si alguien sabe cómo, podría verse el Toporaster completo con su precioso fondo de sombras en el mismo mapa, tal y como os muestro en estas capturas conseguidas superponiendo la capa sobre el fondo. Una preciosidad.

CASTILLA Y LEÓN

- Añadido Topográfico Urbano 1:1000. Magnífico mapa visible a partir de zoom 20. Si el mapa sale algo deformado cambiad la proyección a geográfica WGS84 en el menú View.

CATALUNYA

- Añadida la cartografía topográfica urbana a escalas 1:1000 y 1:2000. Visible a partir de zoom 20 en los zooms anteriores se visualiza un esquema con las zonas disponibles. Ejemplo a zoom 21.

ARAGÓN

Modificados los parámetros del servidor pues había dejado de funcionar la carga de los cuatro mapas disponibles.

MDT: de Global Mapper a CompeGPS.

Recientemente hablamos en este artículo sobre los datos de relieve ASTER GDEM y la forma de conseguirlos desde su página web. Hoy vamos a trabajar un poco con ellos desde Global Mapper ya que, aparte de poder abrir los archivos físicos descargados, tiene la característica de tener una conexión predeterminada con el servicio web de ASTER. La idea es que, en vez de tener que descargar previamente los datos desde su página web, más luego manipularlos y recortarlos a nuestro gusto, podamos elegir directamente la porción de terreno que más nos interese sobre la superficie de los datos y descargarlos, sin necesidad de hacer la descarga primero de los archivos desde su página web.

Además, una vez conseguida la zona que queramos, importaremos el relieve en CompeGPS para salvarlo en su formato propio .cdem (que no se si es mejor o peor pero desde luego en relación calidad/tamaño es insuperable).

Lo más complicado con lo que me he topado en este proceso es el asunto de las proyecciones. En esa locura difícil de entender para todos me resultaba casi imposible que Compe me colocara el DEM en su sitio (siempre teniendo como referencia el servicio wms raster del IGN). Tras muchas pruebas y dudando aún de la fiabilidad del resultado, estos fueron mis pasos:

1. Abrimos en Global Mapper la conexión web predeterminada a los datos ASTER GDEM.

2. Con nuestra zona de interés encuadrada en pantalla, seleccionamos la zona que vamos a exportar. Puede hacerse cualquier polígono con la herramienta vectorial o utilizar algún archivo vectorial que ya tengamos. En este caso hemos cargado el vectorial shape con los recintos de las Comunidades Autónomas de España (disponible en muchos lugares, por ejemplo, en la página del Equipamiento Geográfico de Referencia Nacional del Centro de Descargas del CNIG). Para este ejemplo lo hicimos con la Comunidad Autónoma de Euskadi.

Datos ASTER encuadrados sobre la Península, con la capa de recintos de Comunidades Autónomas cargada y la Comunidad de Euskadi seleccionada. Para seleccionarla podemos usar la herramienta Info o la herramienta Digitizer.

3. Exportamos los datos seleccionados al formato .asc (ya que CompeGPS en principio lo importa sin problemas). La versión web de ASTER nos da una resolución algo menor a la de los archivos descargados de su web como vemos en la zona naranja (parece ser que no sirve los datos a la misma calidad). Por supuesto en este paso el parámetro de la resolución podríamos modificarlo a nuestro gusto y conseguir de esta forma un MDT de 40, 60, 80 metros o la resolución que le pongamos. Muy importante marcar la casilla Crop to Selected Area Feature(s) -cortar al área seleccionada- en la pestaña Export Bounds si no queremos que nos exporte todo el planeta.

Global Mapper nos pedirá un nombre y una ubicación donde exportar nuestro archivo .asc. El proceso puede alargarse más o menos dependiendo de la velocidad de conexión, la resolución puesta, la disponibilidad del servidor, la potencia de vuestro pc, etc...

4. Global Mapper nos genera nuestro archivo de relieve con la proyección del original, que en este caso es la Mercator de Google (EPSG:3857), y que a mi me ha costado tanto que CompeGPS me la coloque en su sitio. Así que he optado por reproyectarlo a una más habitual para nosotros. Cerramos todo lo que tengamos abierto en Global Mapper y abrimos nuestro recién creado archivo .asc. La manera correcta de corregir la proyección es a través del menú Tools/Configure... y en la pestaña Projection modificamos los parámetros como se ve en la imagen (Proyección UTM, zona 30 y Datum ETRS89. He usado estos porque son los utilizados por el IGN en los MDT que nos ofrece para descarga).

Reproyectamos nuestro archivo de relieve con los parámetros indicados. Aplicamos y OK.

Por supuesto una vez cambiada la proyección, exportamos el archivo de nuevo sobrescribiendo el que teníamos para modificarlo.

La resolución del MDT nos ha mejorado (no me preguntéis por qué, supongo que efectos de la reproyección). Esta vez no tenemos nada que seleccionar, así que dejamos marcado All Loaded Data para que exporte todo lo que tenemos en pantalla (nuestro MDT de Euskadi).

Cerramos todo de nuevo y abrimos nuestro MDT reproyectado. Queremos comprobar su precisión sobre el servicio WMS de cartografía raster del IGN antes de irnos a CompeGPS.

Nuestro archivo EUSKADI.asc reproyectado (en amarillo su proyección y resolución actuales) abierto junto con el servicio wms raster del IGN sobre el puerto de Bilbao. La exactitud es impresionante.

5. Importamos nuestro MDT en formato .asc en Compe. Y automáticamente se nos pregunta en qué formato queremos guardarlo (elegimos cdem) y en qué proyección (rellenamos con los datos de nuestro MDT: UTM zona 30 ETRS89).

CompeGPS comienza a generar nuestro MDT en su formato .cdem y al terminar lo muestra abierto en pantalla.

Nuestro archivo abierto ya en formato CDEM y con su ventana de propiedades mostrando resolución y proyección.

6. Hacemos la misma comprobación de exactitud sobre el wms del IGN con la herramienta transparencia.

La precisión del MDT vuelve a ser perfecta.

Y, para terminar, nos damos el lujo y el homenaje de hacer una vista en 3D sobre la ortofoto con nuestro cdem.

Ya tenemos nuestro MDT en formato propio cdem de CompeGPS de todo Euskadi a una resolución según se nos dice de 28 metros. Y el cdem ocupa sólo 7,41 megas!!! Desde luego una relación de tamaño impresionante en comparación con otros formatos (algún truco me sospecho...).

Notas mentales:
- CompeGPS me transforma a CDEM sin problema los archivos .asc de tamaño relativo (Rioja, 95 megas; Euskadi, 150 megas; Navarra, 238 megas), pero por ejemplo con tamaños como el de Canarias (364 megas) y superiores me da error a la hora de fabricar su CDEM. Desconozco si es problema de mi PC o un bug de Compe. Para este artículo he utilizado las versiones Global Mapper 16.1 y CompeGPS LAND 7.6.6.
- El lío de las proyecciones es tan tremendo que supongo que muchos estaréis igual.... Yo he hecho la prueba con estas y no parece ir mal, pero no estaría de más que algún entendido en el asunto hiciera un sabio comentario que nos colocara en nuestro sitio.