Trabajo final de la materia "Introducción a las técnicas inteligentes de resolución de problemas de planificación, secuenciación y ejecución". Maestría en Aplicación de Información Espacial (CONAE - FAMAF)
- Notebook con el código explicado: Notas y código.ipynb
- Archivo Readme: README.md
- Imagen NDVI resultante: ndviImage.tif
- .gitignore
Para este trabajo se implementó un programa en Python para el cálculo de NDVI en un cuerpo de agua.
Se implementó a partir de un archivo zipeado del recorte de la zona de estudio de una imagen satelital Landsat 8, en una fecha dada. El archivo zip contiene las bandas ópticas del sensor.
La ejecución se conforma de 5 pasos:
- Carga de las librerías
- Carga de las imágenes
A partir de un archivo .zip que contiene todas la bandas del sensor se procede a su extracción - Exploración de las imágenes
- Lectura de metadatos (dimensiones, sistema de coordenadas de referencia).
- Visualización de las bandas.
- Histogramas de las bandas.
- Cálculo de outliers y enmascaramiento de los mismos
- Cálculo de MNDWI
- A partir de las bandas 3 y 6 se construye el índice
Fórmula de referencia MNDWI = (Verde – SWIR)/(Verde + SWIR) - Construcción de máscara de valores superiores de MNDWI para recortar el área del lago.
- A partir de las bandas 3 y 6 se construye el índice
- Recorte de bandas con la máscara de MNDWI
- Se recortan las bandas 4 y 5 en base a la máscara de MNDWI.
- Cálculo de NDVI y exportación de imagen
- A partir de las bandas 4 y 5 recortadas se calcula el NDVI.
Fórmula de referencia NDVI = (NIR – Rojo)/(NIR + Rojo)
- A partir de las bandas 4 y 5 recortadas se calcula el NDVI.
Para la ejecución demo se utilizó una imagen del 14-05-2023 y se selecciónó como zona de estudio el Embalse Los Molinos. La zona de estudio corresponde a un subset espacial, en forma de cuadrilátero, que contiene el lago.
Construido y ejecutado en Python 3.11.7
- numpy
- matplotlib
- rasterio
- zipfile
- https://eos.com/es/make-an-analysis/ndvi/
- Xu, H. (2006). Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. International Journal of Remote Sensing, 27(14), 3025–3033. https://doi.org/10.1080/01431160600589179
Biol. Manuel Zeballos