Mostrando entradas con la etiqueta cdem. Mostrar todas las entradas
Mostrando entradas con la etiqueta cdem. Mostrar todas las entradas

lunes, 29 de junio de 2015

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.

jueves, 26 de septiembre de 2013

Modelos Digitales del Terreno.

Llega la hora de hablar de los Modelos Digitales del Terreno (MDT en adelante), cuyo acrónimo en inglés es DEM (Digital Elevation Model). Por cualquiera de las dos formas los conoceréis.
Como sabéis ya hemos conseguido los mapas de la zona que nos interese; ahora es importante hacerse con los MDT de la misma zona, pues querremos representarlos en tres dimensiones para darle el máximo realismo que podamos a la representación virtual del terreno. Para esto ya necesitaremos un software específico que sea capaz de colocar juntos los mapas y los MDT, y sea capaz también de representarlos. En la página de Software básico tenéis mi recomendación: CompeGPS Land. Un programa todo en uno con el que pocas cosas más necesitaremos.
Lo primero e ideal es tener un archivo MDT base que nos cubra toda la Península; dicho archivo lo tenéis en la sección Descargas en formato propio de Compe (.cdem) y abierto en el programa tiene este aspecto. (Siempre que los archivos estén en los formatos comunes que maneja CompeGPS el programa los trata como mapas, por lo que abrir un relieve es el mismo proceso que abrir un mapa).


Y este otro aspecto, en una zona montañosa, en tres dimensiones.

Bien, lo importante de un MDT es básicamente su resolución, esto es, la distancia entre puntos que contienen altura. A menor distancia, más resolución y realismo tendrá el relieve. El propio CompeGPS al preguntarle por las propiedades del archivo nos da la resolución de esa malla que forma el MDT. En este caso nos dice que el tamaño de celda es de 75,1 x 92,8 metros.


Este archivo será suficiente para la inmensa mayoría de nuestros propósitos, pero vamos a aprender también a buscarnos y hacernos nuestros propios MDT. En la página Catálogo de productos del IGN se nos ofrecen tres opciones de descarga para MDT's: 
- Modelo digital del terreno con paso de malla de 200 m, con distribución por provincias.
- Modelo digital del terreno con paso de malla de 25 m, con la misma distribución de hojas que el MTN50.
- Modelo digital del terreno con paso de malla de 5 m, con la misma distribución de hojas que el MTN50.

Vamos a tirar por el camino de enmedio y hacernos con el MDT de 25 metros de la hoja que nos interesa que, por seguir con el ejemplo, será la 157 del MTN50. Tras la consabida búsqueda avanzada en el centro de descargas del IGN, tendremos delante listo para descargar el archivo que buscamos.

Como veis el formato de archivo que nos ofrece el IGN para estos archivos es .asc. No siendo el propio de CompeGPS ( es .cdem como dijimos) el programa es capaz de importar este formato y guardarlo en el suyo propio. Hay multitud de formatos de archivo MDT, tampoco os volváis locos al principio porque con el tiempo es cuestión de controlar algún conversor e ir guardándolos en el formato que más nos interese dependiendo del sotware que utilicemos.
Abrimos CompeGPS, menú Mapas/Importar Relieve 3D.

Elegimos el archivo descargado (en este caso al estar la hoja 157 en zona de huso 29, el zip del IGN contiene dos versiones del archivo .asc, una con huso 29 y otra con huso 30) e inmediatamente CompeGPS nos pregunta en qué formato queremos guardar el archivo .asc en la misma carpeta. Nos da a elegir entre .hdr (otro formato clasico de MDT's) y .cdem. Elegimos este último por ser el propio de CompeGPS y con el que mejor trabaja.

Por último CompeGPS nos pregunta los datos geográficos del archivo; podéis dejarlo todo como en la imagen sólo teniendo cuidado de que la casilla zona corresponda con el huso del MDT que hayáis importado. En este caso yo estoy importando el .asc con huso 29, pues pongo 29 en la casilla zona.

Tras unos segundos en los que el programa ejecuta la conversión, nos aparece el MDT de la hoja 157 del MTN50, con proyección en huso 29 y una resolución de paso de malla de 25 x 25 metros.

Y en la carpeta que contenga el archivo .asc tendréis también otro con idéntico nombre y extensión .cdem. Por cierto el archivo .asc ya no sería necesario y se puede borrar; CompeGPS ya sólo necesitará en el futuro el archivo .cdem para funcionar.
Deciros por último que si la zona que buscamos coincide entre dos hojas del MTN, o entre tres, podemos seguir este proceso descargando los MDT´s correspondientes, importando en CompeGPS, y una vez abiertos todos los que necesitemos en pantalla el programa nos da la opción de unir los MDT en uno solo. Lo guardamos como .cdem y listo.
Un poco más de información sobre los MDT aquí.