Proceso para obtener la reflectancia del raster de Landsat 8 Archivos a utilizar: -
Raster de Landsat 8 con todas las bandas, Herramientas de cálculo para Landsat 8 agregadas al ArcToolbox:
Se agrega ambas -
; (están en sidweb en un rar)
Espacio de trabajo, puede ser un cantón o una provincia en shapefile o vertor
1. Se corta el raster de manera que se obtenga solamente el espacio de trabajo que nos interesa con la herramienta agregada al ArcToolbox “Pan-Sharpened”.
Se escoge la carpeta completa donde están todas las bandas del LC08
Se busca el Mtl dentro de la misma carpeta.
El área es el shapefile o lugar que tomamos para trabajar.
El nombre que asignemos será un sufijo al nombre con el que ya está nombrada cada banda.
este es el mensaje de que se ha terminado el procedo, no se exportarán los archivos porque ya se guardan automáticamente dentro de la carpeta de LC08
2. Se obtiene la radiansa de cada una de las bandas, para ello se usa otra de las herramientas que se agregó anteriormente “Radiance with atmospheric correction”:
Similar en Workspace se ubica la carpeta del LC08 Mientras que en Meta data… se toma el archivo en txt que fue generado por el Pan-sharpened. Por otro lado el nombre que asignemos seguirá siendo un sufijo.
3. Se obtiene la Reflectancia de la misma forma con la herramienta de “Reflectancie with atmospheric correction”
Completamente igual a obtener la radiansa
Es posible observar cómo se agregaron todos los archivos deseados:
4. Probar combinaciones en la composición generadas por el Pan – sharpened: Notese que estas combinaciones esta dado con las bandas aun no corregidas
El color verde representa la vegetación en el lugar
En esta última se observa que el color rojo corresponde a la vegetación 5. Se Obtiene el NDVI (Índice de Vegetación de Diferencia Normalizada) Con respecto a las bandas de reflectancia correspondientes a las Bandas de Nir Y Red que en el caso del LT8 son las bandas 5 y 4 respectivamente
Se añade las reflectancias de las bandas 4 y 5 y usando raster calulator se realiza la siguiente operación: NDVI= (NIR-R)/(NIR+R)
Este último raster representa la imagen de la vegetación en esa área; y para obtener el área específica se debe transformar ese raster a un polígono o vector:
Para saber el valor neto de la vegeatación se debe borrar los valores que tienen cero ya que no representan mayor parte de la vegetación. Se edita el Shapefile que se creó
Donde del resultante es posible conocer el área de vegetación total que existe en el espacio de estudio. 6. Composición de bandas se utilizará para comparar con la de otra fecha en función de la reflectancia antes calculada.:
Se pueden realizar un sin número de combinaciones entre las bandas a mostrar y el color, con el fin de visualizar de mejor manera la vegetación en el territorio.
El archivo que se crea es temporal, asique si lo requieren es necesario exportarlo a alguna dirección.
El orden al inicio es el orden que se las seleccionó Si se altera ese orden se obtiene: