Temperatura de brillo a partir de imágenes Landsat 8

En otros posts he escrito acerca de la determinación de la temperatura de brillo usando la banda térmica (b6) de imágenes Landsat 5. Hoy voy a describir cómo hacer lo mismo pero usando imágenes Landsat 8. Los recursos que pone a la disposición el USGS (U.S. Geological Survey) para el proceso de downloading permite a los usuarios registrados ubicar fácilmente el dataset completo (11 bandas más el archivo *.txt de la metadata) si se cuenta con un shapefile de la zona de interés. Por otra parte, este archivo de metadata es muy explícito en cuanto a los parámetros que hay que seleccionar e, incluso, existen tutoriales bastante ilustrativos acerca de como hay que proceder publicados por el propio USGS y el Blog de ESRI.

Como los tutoriales anteriores están en inglés vamos a describir un procedimiento usando ArcGis 10.1. Sin embargo, la temperatura de brillo se determinó primero para la banda 10 del set LC80380322013144LGN00 (Utah, USA) mediante la acción combinada del raster calculator de QGIS y gdal_cal.py en GNU/Linux Debian (el script de python se usó porque la acualización a Chugiak me imposibilitó el uso de Commander para acceder a la función logaritmo natural necesaria en la fórmula). Voy a usar está imagen como calibración del procedimiento en ArcGis.

Si se leen el tutorial de Arcgis no se si perderán mucho tiempo tratando de descifrar la sistematización allí planteada usando el ‘Raster Function Template Editor’. No obstante, eso sale mucho más fácil en un sólo paso con el Raster Calculator de ArcGis. En Geoprocessing -> Search For Tools (ventana que aparece a la derecha) ubican el Raster Calculator y en el espacio para la fórmula colocan la siguiente (que pueden usar como plantilla para otros casos):

1321.08/Ln((774.89/((0.0003342*"LC80380322013144LGN00_B10.TIF")+0.1))+1)

la cual es válida sólo para las bandas 10. En la línea anterior ya se incluye la metadata correspondiente para transformar los DN en Radianzas que son necesarias para el cálculo de la temperatura de brillo.La imagen de como luce el Raster Calculator es la siguiente:

landsat8_2

Después de OK, se obtiene una imagen dataset de ArcGis cuyo formato incluye 6 archivos *.adf y uno *.xml. Como mi QGIS soporta este tipo de formato y cuenta con el plugin ‘Value Tool’ examiné ambas imágenes simultáneamente, la producida por ArcGis y la obtenida en Debian, y los resultados fueron idénticos; tal como se observa en la siguiente imagen:

landsat8

Los valores que se reflejan en ‘Value Tool’ corresponden a un punto arbitrario en una de las zonas más calientes de la imagen.

Esta entrada fue publicada en ArcGis10.1, Teledetección. Guarda el enlace permanente.

3 respuestas a Temperatura de brillo a partir de imágenes Landsat 8

  1. Pingback: Temperatura de brillo a partir de imágen...

  2. Hola Jose! te consulto, porque estoy haciendo mapas térmicos usando las bandas 10 y 11 del Landsat 8 y me surgen dos problemas: el primero es que no logro que las referencias del raster de Tb me de datos de temperaturas en las diferentes zonas del mapa. Por otro lado, al querer ponerle color al raster, para que me identifique zonas frias y calientes, solo me tira el mapa con diferentes tonalidades de verde. Le cambio el rango de datos y aún así no es coherente. Cualquier ayuda gracias!

    • Los problemas surgen posiblemente porque no estás aplicando adecuadamente la fórmula. Para las bandas 10 tiene que ser similar a esta en el ráster calculator:

      1321.08/Ln((774.89/((0.0003342*"LC80380322013144LGN00_B10.TIF")+0.1))+1)
      

      Para las bandas 11 varía un poco porque las constantes son diferentes.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s