Descargar DEM de 30 metros

En la red existe un gran numero de fuentes de datos disponibles, pero en este articulo me voy a centrar en una que me ha dejado sorprendido y que veo que tiene poca difusión.

Son los modelos de elevaciones GDEM.

Estos modelos de elevaciones han sido desarrollado en conjunto y puesto a disposición del público por el Ministerio de Economía, Comercio e Industria (METI) de Japón y la Administración Nacional Aeronáutica y del Espacio (NASA). Se generan a partir de datos recogidos por el satélite ASTER. Ya se encuentran en la versión 2 y esto implica un mejoramiento constante del producto entregado.

Las características principales se encuentran en la pagina pero aquí las traduzco.

• Escenas de un Arco de grado
• Resolución del pixel 30mx30m
• Coordenadas Geográficas, Datum WGS84
• Formato Geotiff, con niveles digitales de 16 bits enteros positivos
• Niveles digitales especiales, como océano y lagos con valor 0 y sin datos -9999.

Las aplicaciones que se pueden hacer son muchas, enumerarlas seria difícil, asi que les dejo un link que muestra algunas.

Y como descargo los DEM de 30 metros???....entramos a la pagina oficial del proyecto. 

http://gdem.ersdac.jspacesystems.or.jp/

Nos registramos (primera vez en el sitio) y después nos logeamos.

Para buscar la zona de interés, pulsamos la opción Search (menú de la izquierda), nos debe mostrar la pagina donde podemos seleccionar los datos a descargar.

podemos mover el mapa para centrarlo a nuestro sitio, luego dar doble click para hacer un acercamiento, y repetimos este doble click hasta que tengamos una visual que nos permita decidir cuales son los archivos a descargar. Picamos el botón Grid , para activar la grilla de selección y luego el botón Start para poder seleccionar sobre el mapa...

Como podrán notar pueden seleccionar todas las que necesiten y/o utilizar las otras opciones para buscar y descargar los datos. Pulsamos el boton Next, donde se listan todos los DEM a descargar,

Nuevamente pulsamos Next , seleccionamos la categoría en la que utilizaríamos los datos, (yo escogí agricultura), y pulsamos Agree (aceptamos los términos) y en la nueva ventana pulsamos Download.

Esperamos a que se procese la información y guardamos en disco.

Estos modelos de elevaciones tiene mayor resolución espacial que los SRTM, por consiguiente las información generada se acerca mas a lo que vemos en la realidad. 

Por ultimo, espero hacer algunos ejercicios para demostrar su poder y les dejo una captura con gvSIG y QGIS

Éxitos...

La trampa del Hardware en los receptores GPS

Saludos a tod@s.

En esta entrada no traigo un manual sobre un programa o algo parecido. Este articulo busca advertir sobre un problema que me ha pasado.

Imagínate la siguiente situación. Compras un súper radio para escuchar emisoras AM/FM. El manual dice que sirve para eso, el vendedor también  lo indica,  tu emocionado, lo compras. Cuando recibes el radio, revisas el manual, lo prendes, sintonizas una emisora AM, y todo ok, pero cuando tratas de sintonizar una emisora FM!!!!...Salta ese mensaje que no te esperabas, "por favor inserte el serial de activación de esta función".

Ese mensaje "por favor inserte el serial de activación de esta función", es lo que llamo "La Trampa del Hardware".

La verdad me parece algo muy triste, confuso y engañoso, que se ve mucho en los receptores GPS (desde ahora GPS). Se compra un GPS que supuestamente sirven para tomar datos L1 raw, L1/L2, L1/GLONASS, RTK, etc, pero que solo seria posible con los seriales adicionales para poder usar esas características. Engañoso porque no es claro desde el principio y sobre todo porque un equipo de estos, vale unos USD $ 3.000, y los seriales para activar las otras funciones valen mas o menos USD $ 1.000, ojo!!!, cada serial. 

Me pregunto, quien quiere compra un GPS, que supuestamente sirve para todo, pero no lo puedo utilizar sino compro los seriales adicionales. Aclaro que no estoy en contra de que esas características tengan su valor, sea cual sea,  lo que me molesta de todo esto es que vendan en GPS por un valor considerable, y que la activaciones cuesten tanto como el equipo. A demás no es claro este tema de la activaciones de funciones por seriales en los manuales o antes de la venta, pereciera que su primer interés es vender el equipo, y después como se hizo semejante inversión casi quedas obligado a comprar los seriales para utilizar el GPS al maximo. 

Para finalizar esta reflexión, te recomiendo que documentes bien el equipo que quieres comprar (sea cual sea), pregunta por todas la funciones que trae, si están activas o necesitas de seriales o programas especiales que solo lo vende quien desarrolla el equipo, pide que prueben el equipo con todas las funciones, recuerda que si es por tu concepto que compran una herramienta, y luego te das cuenta que necesitas comprar cosas adicionales, no vas a quedar bien parado, también ten presente, que "La Trampa del Hardware", es muy  común en la mayoría de los equipos para trabajos especializados.

Espero que muy pronto podamos tener opciones al estilo gvSIG, Qgis, Mapserver, Geoserver, Postgres+Postgis, OpenLayers, Spatialite, Grass, OrfeoToolbox, en cuanto a dispositivos o herramientas para tomas de datos, en otras palabras lo que por ahí, ya se conoce como el OpenHardware...

Espero les sirva esta información.

Clasificación No Supervisa gvSIG y Condicional en la Calculadora Raster

En este ejemplo se hace una clasificación No-Supervisada con gvSIG y Sextante y después solo se extrae la cobertura de interés por ejemplo los cuerpos de agua.
Un poco de teoría. Existen diferentes métodos o técnicas que permiten extraer los elementos presentes en una imagen. Para nuestro caso hablo de las cobertura en una imagen de satélite. Esos métodos de clasificación se agrupan en dos: supervisados y no supervisados, a nivel general.

Clasificación supervisada[1]

Este tipo de clasificación supone que se tiene un conocimiento de la zona, osea somos capaces de proporcionar o indicar cuales son la muestras (pixeles en diferentes bandas) para entrenar el algoritmo. Luego del entrenamiento el clasificador es capaz de determinar a que grupo o clase pertenecen el resto de elementos (pixeles), y así generar una capa que representa, por ejemplos las coberturas de una zona.

Clasificación no supervisada[1]

A diferencia de la supervisada se supone que no contamos con conocimiento de la zona, por lo que no podríamos entrenar los algoritmos para la clasificación supervisada. En general en este tipo de clasificación parte de la propiedades estadísticas y de vecindad que tiene la información (pixeles en diferentes bandas) y se van agrupando formando clases de manera iterativa. La clasificación no supervisada se la suele llamar también clustering(agrupamiento).

Veamos los pasos.

1- Abrimos gvSIG y cargamos nuestra imagen de satélite. Para este ejercicio utilizo una imagen Landsat, realizo una composición RGB(432).

2- Abrimos Sextante y buscamos Herramientas de análisis para capas raster/Clasificación no supervisada(clustering)...

3- La configuramos así, en la pestaña Parámetros recomiendo seleccionar tres bandas, para este ejercicio la 4, la 3 y la 2. En numero de clases tomo 10 pero dependiendo de la imagen pueden ser menos o mas. Solo tener en cuenta que es mejor que el clasificador agrupe lo mas acertado posible y esto se logra con un buen trio de bandas y suficientes clases. En otras palabras si mi me interesan tres clases de la imagen es mejor clasificarla en 7 y después escoger la clases que me forma mis elementos de interés.

4- en la pestaña Parámetros Región de Análisis la opción que mejor resultado me  proporcionado es  utilizar la extensión de otra capa, después de esto pulsamos el botón aceptar.

5- Una muestra del resultado.

Entendiendo un poco, los que tenemos es una capa resultado que agrupa en 10 clases los pixeles de las bandas 4,3 y 2 que participan en el proceso. A cada clase le asigna un numero como identificador y sobre los cuales podemos hacer diferentes operaciones.
Hagamos el ejercicio de extraer los cuerpos de agua con ayuda de la calculadora raster. procedemos así,
1- Teniendo nuestra capa resultado activa, seleccionamos el botón de información (botón con letra i de color blanco sobre un circulo azul) y picamos lo que consideramos como cuerpo de agua. así podremos saber cual es el valor del pixel asignado a esta clase. Para nuestro ejemplo ahora sabemos que los cuerpos de agua tiene el valor 1

2- En Sextante buscamos Herramientas de calculo para capas raster/calculadora de mapas...

3- Escribimos la sintaxis para extraer el valor. Esto realmente se puede hacer de otras formas pero quiero que veas la alternativa con el condicional if(), básicamente funciona así, if(condición, proceso 1, proceso 2), el proceso 1 se ejecuta si la condición es verdadera de lo contrario (osea si la condición es falsa) se ejecuta el proceso 2. Nuestro condicional quedaría como , if(Agregados Band 1 == 1, 1, 0). Puedes anidar varios condicionales.

4- Podemos observar el resultado. De una forma muy rápida se extrajo el cuerpo de agua, aunque en el resultado muestra otras zona que seguramente no lo son, esto se puede mejorar con otra combinación de bandas, o aumentando el numero de clases, también puedes vectorizar y eliminar los polígonos que no necesitas.

Saludos y espero que les sirva...
Tarea: Trata de hacerlo con mas coberturas anidando por ejemplo dos condicionales, te explico como.
if(condición1, proceso 1, if(condición2, proceso 3, proceso 4)proceso 2), el proceso 1 se ejecuta si la condición1 es verdadera de lo contrario (osea si la condición1 es falsa) se  evalúa la condición2, si esta es verdadera se ejecuta el proceso 3 y si la condición 2 es falsa se ejecuta el proceso 4.
Dejo un vídeo para mayor facilidad

[1] http://es.wikipedia.org/wiki/Reconocimiento_de_patrones

Pequeño Resumen de Datos Libres

Cuando estamos realizando una investigación para la universidad o laboramos en una entidad publica o privada, y necesitamos datos, es cuando por lo general empiezan los problemas. Valen mucho, son restringidos, no han aprobado el presupuesto, en fin, obstáculo tras obstáculo. Te invito a que esos inconvenientes no te detengan y menos en esta época donde existen instituciones que han puesto información muy valiosa para el uso del publico. Esta es una lista que puedes ir construyendo haciendo un buen uso de los buscadores, sin embargo acá te dejo un resumen que me gustan mucho porque los considero básicos a para cualquier trabajo.

• Si necesitas los nombres geográficos de tu área de interés o de un país completo, esta es la fuente indicada, descargas los archivos en (http://www.geonames.org/), vienen codificados por pais y formato de texto. Recuerda que te puedes ayudar de estas entradas para saber como utilizarlos en cualquier programa.
• Para los limites administrativos de departamentos y municipios (términos para Colombia) puedes decargas los SHP desde esta pagina (http://www.gadm.org/): Sin embargo hay que aclarar que la información a escala municipal esta muy generalizada, pero es mejor que nada.
• Es posible que necesites determinar las coberturas de una zona, comparar imágenes de diferentes años, extraer o actualizar las vías y zonas urbanas o simplemente generar una mapa con una imagen de satélite como fondo. Si una imagen Landsat cumple con tus requerimientos puedes mirar aquí para que sepas como descargarlas y aquí como armarlas.
• Uno de los datos mas costoso en general en los SIG, es la altura sobre todo de grandes extensiones, pero en esta época, también los hay de acceso libre. Estos Modelos de elevaciones los puedes bajar como se indica aquí, y a partir de ellos puedes determinar cuencas, red hidrico, capas de pendiente, de aspecto, modelos de sombras, de visibilidad, o un gran ETC!!!.

Y si quieres un gran resumen de datos a nivel mundial te dejo esta pagina, mejor imposible... http://freegisdata.rtwilson.com/.

Pero bueno tenemos datos y como los trabajamos, que programa compramos??. Me atrevo asegurar que si eres una persona curiosa con ganas de aprender y de tener alternativas te encantara esta otra lista de programas, la mayoría libre que puedes utilizar con diferentes propósitos. Si estas desesperado porque la lista es larga simplemente te comento que puedes usar gvSIG y/o QGIS para que inicies.

Este pequeño resumen este pensado para trabajos puntuales, cabe aclarar que si tu idea es un algo mas corporativo también existen alternativas libres que se pueden combinar con algo de que menciono arriba. Finalizo  diciendo que los limites son los datos en si (querer genera un mapa de coberturas a escala 1:10.000 con una imagen Landsat) y nuestros conocimientos, el primero es insalvable, toca cambiar la fuente de información, el segundo lo podemos mejorar, solo depende de las ganas, de la disposición, de la comunicación con otras personas.

Exitos!!!

Instalar gvSIG 1.11 en Ubuntu 12.04 de 64 Bits

Saludos a todos

Esto para mi fue un dolor de cabeza pero bueno, se pudo solucionar. Y lo comparto para que le ahorre tiempo alguien.

Paso a explicar lo que me pasaba.
- Lo primero era descargar el instalador de gvSIG, los cual no fue mayor problema.
- Asignarle a ese instalador permisos de ejecución que tampoco fue inconveniente.
- Ejecutar el instalador...aqui aparecía el primer problema, este era el error que mostraba la consola.

/home/humberto/Descargas/gvSIG-1_11-1305-final-lin-i586-withjre-j1_5.bin: 12: /home/humberto/Descargas/gvSIG-1_11-1305-final-lin-i586-withjre-j1_5.bin: ./launcher-Linux: not found

/home/humberto/Descargas/gvSIG-1_11-1305-final-lin-i586-withjre-j1_5.bin: 21: /home/humberto/Descargas/gvSIG-1_11-1305-final-lin-i586-withjre-j1_5.bin: ./launcher-Linux: not found

Lo que pude entender era que faltaban unos paquetes de compilación (eso creo) que se son los siguientes:

sudo apt-get install g++-multilib
sudo apt-get install gcc-multilib

Aquí desapareció el primer error pero surgió otro. Lo muestro a continuación.

error while loading shared libraries: libXi.so.6

Este problema surgía porque estaba instalando una maquina java de 32 Bits que estaba buscando librerías de 32 Bits y mi instalación era de 64 Bits. La solución instalar las librerías de 32 Bits para el sistema, usamos el siguiente comando:

 sudo apt-get install ia32-libs

Nuevamente ejecutamos nuestro instalador....y el resto ya lo conocen.....

Espero les sirva y como dije al principio les ahorre tiempo.
Ref primer error
Ref segundo error

Organizar datos de Coordenadas geográficas Parte II

ACLARACIÓN: iniciando el desarrollo de esta parte me di cuenta que gvSIG solo entiende el separador de decimales como punto. Aquí puedes hacer varias cosas pero te menciono dos, o cambias la configuración regional del equipo y exportas la tabla, o al resultado del archivo CSV cambias las comas por puntos con la herramienta de Buscar y Remplazar de cualquier editor de texto (no de tablas) y salvas.

Para la parte II haré  el ejercicio de cargar la tabla con gvSIG y QGIS.
Mis datos de tabla están en WGS84.  Guardamos la tabla  en formato CSV, delimitado por punto y coma.

En este punto hemos exportado la tabla a un formato de intercambio. Procedemos a cargarla en gvSIG, desde el Gestor de proyectos, picamos el elemento Tabla y pulsamos el botón Nuevo, en la nueva ventana, pulsamos el botón Añadir, buscamos la tabla teniendo presente que archivos tipo este en CSV, pulsamos Abrir,  la tabla queda cargada y finalizamos esta parte pulsando el botón aceptar.

La tabla se despliega de forma automática.

Teniendo una Vista abierta, nos vamos al menú Vista/Añadir capa de eventos...y asigamos los campos de Latitud y Longitud donde corresponden.

Pulsamo Aceptar y los puntos deben aparecen el vista...

Ya aquí si queremos podemos guardar esta información como Shape file, realizar operaciones, consultas etc.
Ahora con QGIS los hacemos de la siguiente forma. Vamos al menú Capa/Añadir capa de texto delimitado.

En el cuadro de dialogo que aparece buscamos nuestro archivo de texto, lo seleccionamos, seleccionamos el delimitador, que para este caso es punto y coma,  configuramos los campos que corresponden a la Latitud y la Longitud, y si todo esta bien, pulsamos Aceptar.

Y si todo esta bien nos debe mostrar algo como esto.

Y si lo que quieres es comparar aquí te dejo una ultima imagen.

Por ultimo les dejos un vídeo.

Organizar datos de Coordenadas geográficas Parte I

Saludos a tod@s.

He notado que se presentan inconvenientes a la hora de poder utilizar tablas que poseen coordenadas en algún programa se sistemas de información geográfica.

Voy a dar unas pautas básicas para organizar y utilizar esta información sobre todo cuando son coordenadas geográficas aunque pueden aplicarse a coordenadas planas.

1 – La latitud y la longitud deben estar en grados decimales. Osea 4°15'2” a 4,2833333333.

2 – La latitud debe estar es una columna y al longitud en otra.

3 – La latitud norte y sur se diferencian en el signo. Positivo norte y negativo sur, 10°20'6” S a -10,4333333333.

4 – La longitud tiene el mismo principio. Este positivo y oeste negativo.

5 – Verificar como entiende tu programa el separador de decimales como punto o coma, sobre todo cuando el programa SIG es capaz de leer los datos desde la hoja de calculo o un DBF.

5 – Por ultimo, los campo de latitud y longitud deben se numéricos. Aunque parezca obvio, se puede cumplir todas características anteriores y el campo estar declarado como texto, sencillamente ningún programa (por lo menos yo no lo conozco) podrá usar las columnas para dibujar las coordenadas de la tabla.

Teniendo presente estas reglas, cualquier información que no las cumpla debe ser estandarizada con ayuda de algún programa, yo recomendaría una hoja de calculo, aunque se puede hacer también con la herramientas para tabla de los SIG.

vamos a los ejemplos.

 

Como pueden observar en la figura de arriba, las columnas están separadas por punto y coma, las coordenadas se encuentran en grados minutos y segundos y ademas la latitud y longitud se mezcla con letras y símbolos. Para que puedan ser usada en algún programa SIG hay que transformar esta información, veamos como.

Utilizando LibreOffice – Calc (también se puede hacer con Excel) importamos el archivo,

Nos queda algo como esto.

Cuando ya tenemos la tabla importada debemos poder convertir los valores de grados, minutos y segundo a grados decimales. Primero utilizamos las herramientas de texto para separar la latitud en tres columnas, una de grados otra de minutos y otra de segundos. Igual se hace con la longitud.

La función utilizada se llama MID(en Excel se llama EXTRAE), la cual extrae de una cadena de texto, una porción partiendo de una posición inicial y una cantidad de caracteres.

Como podrán notar en algunas columnas hay caracteres y números, esto se puede solucionar de manera rápida con la herramienta Buscar y Remplazar, sobretodo si son muchos datos. Y finalmente la información queda así.

Ahora bien solo nos queda pasar de grados minutos y segundos a grados decimales. La formula general es: grados decimales = ((segundos/60) + minutos)/60 + grados, aplicando esta formula en nuestra tabla tenemos.

Tenemos la información estandarizada, ahora a probar con algún programa SIG. Tenga en cuenta que lo mas probable es que tengas que guardar la información en otro formato, como dbf o csv. Voy a hacer el ejercicio con gvSIG y Qgis dos programas excelentes, y la tabla en formato csv. Espera la parte II
A continuación le dejo un vídeo de lo explicado y que también puedes ver en Youtube

gvSIG y Calculadora Raster de Sextante

Hace algunos días un compañero de la lista de usuarios gvSIG tenia una necesidad, trate de ayudarlo dándole las indicaciones pero se mi hizo difícil. Así que hice un vídeo para que me entendiera mejor. El vídeo consiste en extraer una porción de información de un archivo raster y luego se  vectoriza.

Saludos y espero que les sirva

Cambiar de UTF-8 a Latin1 en Postgres 8.4 con pgAdmin III

La idea es mostrar otra alternativa a las que normalmente se encuentran buscando en la web. Empecemos.

Primero abrimos nuestro  pgAdmin III y buscados dentro de las BD existentes una llamada postgres, entramos. Selecionamos Catalogos/PostgreSQL/Tablas/pg_database, la abrimos. En la columna encoding, se encuentran los códigos de cada una de la BD creadas.

 Para este ejemplo 6 corresponde a UTF-8 y 8 a Latin1. Si seleccionamos la BD que nos interesa desde el Explorador de Objetos se muestra en el panel SQL el tipo de encoding.

Para proceder a cambiar el encoding simplemente vamos la tabla pg_database, en la columna encoding de la BD de interés  y cambiamos 6 por 8, guardamos los cambios.

Si queremos nuevamente podemos ver la BD en el Explorador de objetos, confirmaremos los cambios.

Para finalizar les dejo una tabla con todos los valores que se pueden tener.
Aclaro que todo esto fue hecho con Postgres 8.4.

Saludos y espero que les sirva.

CODIGO	NOMBRE	DESCRIPCION
0	SQL_ASCII	ASCII
1	EUC_JP	Japanese EUC
2	EUC_CN	Chinese EUC
3	EUC_KR	Korean EUC
4	JOHAB	Korean EUC (Hangle base)
5	EUC_TW	Taiwan EUC
6	UNICODE	Unicode (UTF-8)
7	MULE_INTERNAL	Mule internal code
8	LATIN1	ISO 8859-1
9	LATIN2	ISO 8859-2
10	LATIN3	ISO 8859-3
11	LATIN4	ISO 8859-4
12	LATIN5	ISO 8859-9
13	LATIN6	ISO 8859-10
14	LATIN7	ISO 8859-13
15	LATIN8	ISO 8859-14
16	LATIN9	ISO 8859-15
17	LATIN10	ISO 8859-16
18	ISO_8859_5	ISO 8859-5
19	ISO_8859_6	ISO 8859-6
20	ISO_8859_7	ISO 8859-7
21	ISO_8859_8	ISO 8859-8/ECMA 121 (Latin/Hebrew)
22	KOI8	KOI8-R(U)
23	ALT	Windows CP866
24	WIN874	Windows CP874 (Thai)
25	WIN1250	Windows CP1250
26	WIN	Windows CP1251
27	WIN1256	Windows CP1256 (Arabic)
28	TCVN	TCVN-5712/Windows CP1258 (Vietnamese)

ARMAR IMAGENES SEPARADAS EN BANDAS EN UNA SOLA CON gvSIG

Hace mucho tiempo no escribo nada en este blog que inicie y deje votado... pero bueno quiero volver y devolver algo de lo que he aprendido de otros que muy amablemente colocan información grandiosa en sus paginas.

Quise regresar con un tip, de esos que ahorran tiempo y aportando algo al uso del gvSIG.

Si por algún motivo necesitas armar una imagen de satélite que descargaste y tienes las bandas por separado, este es un procedimiento para unirlas en un solo archivo. Aclaro que solo lo he hecho para bandas en formato tif, para formato Raw puedes mirar aqui

Después de abierto gvsig, crea una vista y configura el sistema de referencia al de la imagen que descargarte. Por lo general para Landsat que es la del ejemplo corresponde WGS 84 con UTM y la respectiva zona. Esto lo aconsejo para evitar reproyecciones al vuelo que pueden demorar el procedimiento.

Procedamos a cargar la primera banda de la imagen

 Luego  nos vamos a la propiedades de la imagen y seleccionamos la pestaña Bandas, como se muestra en la imagen solo aparece una banda. 

Pulsamos el botón Añadir se abre un cuadro de dialogo, seleccionamos las bandas restantes y tenemos cuidado de hacerlo en orden para que así queden en el archivo final.

Ya tenemos casi todo listo, la gran ventaja de este procedimiento es que en este punto ya podemos hacer combinaciones de bandas, realces, etc.

pulsar el botón Aplicar y luego Cerrar del dialogo propiedades. Nos debe quedar algo como esto.

Para hacer permanente esta unión de bandas click derecho sobre la capa y escojemos Salvar como...no es necesario modificar nada, le damos el nombre de salida y pulsamos Guardar...esperamos y lo agregamos a la vista para verificar.

Espero les sirva esta información...

USGS LIBERA TODA SU COLECCIÓN DE IMÁGENES LANDSAT!!!

Es una gran noticia, el Servicio Geológico de los Estados Unidos liberara toda su colección de imágenes Landsat. Este es un aporte incalculable, datos desde los 70's.
En enero del año 2009 el USGS planea tener a disposicion toda su coleccion de imágenes Landsat. Desde Landsat-1, hasta Landsat-7. Los tiempos establecidos son los siguientes, para Julio de 2008 todas las nuevas aquisiciones de Landsat-7, Septiembre de 2008 todos los datos de Landsat-7, Diciembre de 2008 todos los datos de Landsat-5, Enero de 2009 el resto de datos Landsat.
Las imagenes pueden ser previsualizadas y descargadas en http://glovis.usgs.gov/

Al tocar la zona de interes es lanzado el visor que permite descargar la informacion....

consulte el documento oficial en este aqui

ARMADO DE IMAGENES LANDSAT GRATUITAS CON SOFTWARE LIBRE

Con un programa Open source es posible armar el archivo que tiene las bandas que componen la imagen de satélite. Este procedimiento sera explicado para landsat pero aplica para cualquier imagen que deba ser armada desde cero.

Lo primero y quizas lo mas importante es ubicar el archivo que trae los metadatos de la imagen. En otras palabras, las dimensiones en filas, columnas, tamaño del pixel, resolucion radiometrica, fecha de la toma, entre otros. Por lo general es un archivo de texto. Como la extension del archivo no siempre es igual, la forma mas facil de indentificarlo dentro de todo el grupo es observar lo siguiente.

Para imagenes Landsat MSS

Para imgenes Landsat TM
Para imagenes Landsat ETM+
Después de descargar la información revisamos cual archivo contiene los metadatos. El posible archivo que quedo por descarte trate de abrilo con un editor de texto. Al abrirlo podemos encontrar algo como esto para una imagen Landsat MSS.

Pero que programa utilizamos. Este ejercicio se puede realizar con muchos programas algunos muy fáciles, como ERDAS o ENVI. Para este caso utilizaremos un programa Open Source, o libre se llama FWTools 2.1.0. Este programa permite crear las imágenes en formato por ejemplo TIFF o IMG entre muchas opciones. Descargar e instalar el programa, ese es el siguiente paso. Ejecute el programa le aparecerá el siguiente visualizador.

En el menú File/Open Raw Image aparecerá un cuadro de dialogo que permitirá hacer el importe de cada una de las bandas descargadas de www.landcover.org

Este cuadro debe ser llenado con unos parámetros básicos, el resto de parámetros pueden quedar por defecto. Image width e Image height, se toman del archivo de metadatos, quedando algo asi.

Pulsamos el botón Save y colocamos al archivo de salida un nombre apropiado que nos permita identificar que banda es. Puede ser B1 o BANDA1. No es una regla solo es una sugerencia. Luego pulsamos el botón OK. Este paso lo repetimos como bandas queramos unir en un solo archivo.

Para este ejemplo solo utlizare tres bandas pero recuerden que se pueden unir más. En el programa nos vamos al menú Image/Compose… en el cuadro de dialogo pulsando el botón Load File buscamos los archivos producidos en el paso anterior. Escogemos b1.vrt. Si todos los archivos se encuentran en el mismo directorio serán cargados automáticamente. Seleccionamos las bandas con el mouse, y pulsamos el botón Add, las bandas deben pasar al recuadro Output. Asegúrese de tenerlas ordenadas. Luego pulse el botón Current View, debe mostrase la imagen en el visualizador.

Teniendo la imagen en el visualizador procedemos a exportarla el formato que necesitamos, como por ejemplo GeoTIFF. El menú File/Export… permite llevar la imagen a muchos formatos. En el recuadro Output coloque el nombre del el archivo de salida acompañado de la extensión, para este caso tif en las opciones Output Format podemos escoger la que necesitemos, dejamos Gtiff, Output Resolution dejamos Full, para que la imagen quede con el máximo de información. Por ahora el resto de parámetros pueden quedar por defecto. Pulsamos el botón Export y tenemos una imagen armada esperando ser procesada e interpretada con este o con cualquier otro programa.

Solo tenga presente algo. Algunas veces las bandas estan en formato TIFF. Para este caso solo es necesario hacer los dos ultimos pasos.