De acuerdo a ESRI la reflectancia aparente permite ajustar los valores del brillo de imágenes (DN) de algunos sensores de satélite, en base a la elevación del sol, la fecha de adquisición y las propiedades del sensor (ganar/ajustes de influencia para cada banda), muy similar a una corrección atmosférica. Este proceso realiza dos correcciones, una basada en la configuración de ganancia, los valores de brillo originales se vuelven a crear a partir de los valores de la imagen al revertir las ecuaciones de ganancia, y la otra corrección tiene que ver con las diferencias en el ángulo del sol y el brillo. Los valores resultantes generalmente son más bajos que los valores de entrada.
ArcGIS 10.1 ejecuta esta función de forma automáticamente para los sensores Landsat, IKONOS y QuickBird, que pueden ser descargados gratuitamente desde EarthExplorer, y que contengan información sobre su fecha de adquisición, elevación del sol para el dataset, ganancia, influencia y el prejuicio e irradiación del sol de cada banda.
Para este ejemplo se ha descargado una imagen ETM Landsat desde EarthExplorer, con todas sus bandas y toda la información necesaria para la corrección incluida en los metadatos clave (somo sugerencia se puede reemplazar la extensión del archivo MTL, de txt a xml).
Desde ArcMap se agrega el archivo XML que contiene los metadatos de la composición multiespectral ETM Landsat, en Pyramid resampling technique seleccionar Bilinear Interpolation, luego aceptar la creación de pirámides.
Una vez obtenida la composición multiespectral, dirigirse a las propiedades, en la pestaña Functions, en Stretch Function, insertar la función Apparent Reflectance Function, y ajustar/aceptar.
Una vez realizado este proceso se cuenta con el ajuste de la reflectancia aparente , posteriormente se puede realizar una corrección radiométrica, clasificación de la imagen, balance de color o creación de mosaicos, es importante señalar que ArcGIS 10 dispone diferentes funciones automatizadas para archivos ráster, para mayor claridad del procedimiento se puede visualizar el videotutorial, también incluye como determinar de forma automática el índice diferencial de vegetación normalizado (NDVI).
NOTA: Debido algunos problemas que han tenido los usuarios para obtener su composición RGB, producido por problemas de lectura, estructura o ausencia del archivo MTLold.txt (cuando se modifica la extensión a *.xlm o *.met), como otra opción, existe la posibilidad de hacerlo manualmente de la siguiente forma:
- Activar Imagine Analysis ubicado en el menú Windows.
- Añadir todas las bandas dentro de ArcMap.
- En Imagine Analysis seleccionar las bandas a combinar, ejemplo RGB 123.
- En Processing ejecutar Composite Bands.
- Luego ya se puede aplicar el procedimiento sugerido en esta entrada.
Excelente aporte, me gustaría saber sobre corrección atmosférica y cálculo de otros índices.
Gracias, espero tener un tiempo libre para preparar algo así, pero si cuentas con ENVI o Geomatica 2013, poseen herramientas para hacerlo de forma sencilla.
Hola. Tengo instalada la versión 10.1, pero no me acepta los xml. Desde el Add Data no me aparecen los xml y si intento pasarlo directamente al layer me marca error. ¿Sabes por qué hace esto?
Gracias
De donde descargaste la imagen, asegúrate que la imagen sea Landsat L7 ETM+ (la estructura de los archivos TIF y TXT, debe ser similar como se muestra en la captura), mira si el archivo de texto MTL no lo reconoce directamente, le cambias la extensión a XML.
A mi me pasa lo mismo, pudiste solucionarlo.
Franz muy buen aporte, facilita mucho el calculo de la reflectancia aparente, puesto que Envi los valores del gain and offset se deben ingresar manualmente. Gracias.
Una pregunta esto también sirve para imágenes rapideye, puesto que tengo el archivo .xml pero al llevarlo al arcmap este me deniega el acceso.
Tengo entendido que funciona con las imágenes MTL de 7 bandas, no he probado con el resto.
Hola, sin duda un gran aporte, quisiera saber si funciona de la misma manera con las imágenes Landsat 8 OLI ya que he intentado subir el archivo .xml pero me lo impide.
Gracias por la ayuda
Hola, en muchas ocasiones tu blog me ha sido de gran ayuda, quisiera preguntarte una duda que tengo aunque no es muy referente a este tema espero que puedas ayudarme. En la imagen final que muestras, la captura de pantalla de arcmap, en tu imagen de los bordos, izquierdo y derecho se ven dos franjas de colores azul y roja, yo descargué unas imágenes de glovis Landsat 5 TM del año 2000 ya que debo hacer una clasificación de vegetación, al momento de unir las bandas también me salen, mi problema es que como utilizo 4 imágenes para mi zona, hago un mosaico en ERDAS para unirlas pero estas franjas aparecen en mi mosaico final y alteran la clasificación ya que esta franja atraviesa toda mi zona de estudio, mi pregunta es; ¿como puedo eliminar esta franja?
Muchas gracias y espero que puedas ayudarme
Puedes hacer un recorte (solo si cuentas con una área de solape extensa) antes de unir. También puedes cambiar el método de fusión en Erdas.
Tengo una imagen landsat ETM+ con nivel de procesamiento L1, las bandas vienen con valores de reflectancia muy altos (positivos y negativos) como las puedo trabajar en IDRISI?