Este tutorial está basado en uno que en alguna ocasión nos resultó muy útil para armar un mosaico de ortofotos (ver entrada original por Vehrka de Geomaticblog); para este tutorial trabajarémos con algunos Modelos Digitales de Elevación (MDE) del INEGI  utilizando GDAL, Kosmo y posteriormente Sextante.

antes de comenzar es necesario tenga instalado FWTools, si ya lo tiene prosiga a realizar los siguientes pasos, sino revise la entrada para instalar FWTools en Windows

1.- Descargue de la página de INEGI los MDE a utilizar

Del sitio del INEGI descargue los MDE de las cartas 1:50,000 E13B11 y E13B12, esto descargará los siguientes archivos:

MDE-n200000s194500e1044000o1050000.zip
MDE-n200000s194500e1042000o1044000.zip

Descomprímalos en la carpeta de su elección y obtendra los siguientes archivos:

MDE-n200000s194500e1042000o1044000.bil
MDE-n200000s194500e1042000o1044000.blw
MDE-n200000s194500e1042000o1044000.hdr
MDE-n200000s194500e1044000o1050000.bil
MDE-n200000s194500e1044000o1050000.blw
MDE-n200000s194500e1044000o1050000.hdr

2.- Revisar los archivos .bil

Kosmo no lee archivos .bil, por lo tanto para validar que los archivos descargados están correctos utilizarémos el siguiente comando para ambos archivos .bil:

gdalinfo MDE-n200000s194500e1042000o1044000.bil

Esto arrojará el siguiente resultado:

Driver: EHdr/ESRI .hdr Labelled
Files: MDE-n200000s194500e1042000o1044000.bil
MDE-n200000s194500e1042000o1044000.hdr
Size is 1200, 900
Coordinate System is `’
Origin = (-104.666666666666671,20.000000000000000)
Pixel Size = (0.000277777777778,-0.000277777777778)
Corner Coordinates:
Upper Left  (-104.6666667,  20.0000000)
Lower Left  (-104.6666667,  19.7500000)
Upper Right (-104.3333333,  20.0000000)
Lower Right (-104.3333333,  19.7500000)
Center      (-104.5000000,  19.8750000)
Band 1 Block=1200×1, ColorInterp=Undefined
NoData Value=0

De aquí podemos conocer algunas cosas sobre nuestro archivo, como el tamaño 1200 x 900, que NO tiene ningún sistema de coordenadas asignado, el tamaño del pixel y las coordenadas extremas, INEGI nos dice que la información se proporciona en coordenadas geográficas, las unidades de Z están en metros, el datum es ITRF92 época 1988.0 y el elipsoide GRS80.

3.- Mosaiqueo de los MDE

Hasta este momento Usted no ha visto en Kosmo los MDE, esto porque todo el tiempo hemos trabajado con el formato .bil, en un par de pasos más podrá ya visualizarlos mosaiqueados en Kosmo.

*Actualización:

El programa gdal_merge.py en la versión que descargaron genera el siguiente error:

File “C:\ARCHIV~1\FWTOOL~1.7\bin\gdal_merge.py”, line 463, in ?
sys.exit(main())
File “C:\ARCHIV~1\FWTOOL~1.7\bin\gdal_merge.py”, line 433, in main
if quiet == 0 and verbose == 0:
nboundLocalError: local variable ‘quiet’ referenced before assignment

Para solucionar este problema tiene que descargar el archivo gdal_merge.py corregido que se encuentra en la página http://lists.osgeo.org/pipermail/gdal-dev/2010-March/024106.html o en  la caja de Archivos  que se encuentra en la GeoComuna y reemplazarlo en la carpeta bin (C:\Archivos de programa\FWTools2.4.7\bin) , una vez reemplazado el archivo gdal_merge.py el error no debe aparecer.

Para generar el mosaico con los dos archivos .bil, utilizaremos el comando gdal_merge de la siguiente manera;

gdal_merge -o <archivo_salida.bil> <mde1.bil> <mde1.bil>

gdal_merge -o mosaico_4326.bil MDE-n200000s194500e1042000o1044000.bil MDE-n200000s194500e1044000o1050000.bil

Al archivo generado mosaico_4326.bil sólo falta transformarlo a GeoTiff y para ello utilizarémos el comando gdal_translate

gdal_translate –of gtiff <imagen_entrada.bil> <imagen_salida.tif>

gdal_translate –of gtiff mosaico_4326.bil mosaico_4326.tif

Ahora abra mosaico_4326.tif en Kosmo y listo!

Si requiere este mosaico en UTM es necesario reproyectarlo, para ello utilizarémos el comando gdalwarp

gdalwarp -s_srs EPSG:4326 -t_srs EPSG:32613 mosaico_4326.tif mosaico_32613.tif

Para esta entrada colaboraron:

Carlos Ruíz, Jesús Medina, Olivia Contreras y Juan José Del Toro