En el SIyAH-NA se han adoptado las bandas 3 (VIS) y 7 (SWIR) del sensor MODIS (platafomas TERRA y AQUA, resolución espacial 500 m), debido a su asociación a la presencia de agua en superficie. Además, resulta una combinación VIS/SWIR muy sensible a la presencia de agua en comparación con las otras obtenidas mediante combinación de bandas VIS y NIR del sensor MODIS (en efecto, es la que produce mayor realce de las zonas anegadas). Esto es, si asumimos que los valores NDWI>=0 están asociados píxeles 'complementamente inundados', la combinación de bandas 3 y 7 permite detectar una mayor cantidad que otras parejas VIS/SWIR (disminuye la omisión y se garantiza una baja comisión por el carácter conservador del método de umbral). Y esto reviste valor considerando que, dada la baja resolución espacial de los datos MODIS y el tipo de asociación entre presencia de agua y NDWI, el carácter rígido de una clasificación por umbral NDWI suele conducir a estimaciones de valor de área saturada por debajo de las que pudieran regitrarse in situ, una vez definida una unidad de agregación espacial (e.g. cuencas vertientes, sistemas palustres). Así, fue selecionada tal combinación de bandas para la generación de un proucto NDWI (NDWI37 MOD09A1), útil para la generación de un estimado asociado al valor de área saturada para una o distintas unidades de agregación espacial.
En principio, se generó la capa PROD_ndwi37_YYYYMMDD (correspondiente al producto PROD de MODIS, fecha YYYMMDD). Luego, a partir de su procesamiento masivo, hemos podido advertir cierta 'volatilidad' en el índice construido con bandas 3 y 7. Particularmente, hemos detectado dos problemas. En primer lugar, la generación de valores nulos al aplicar el índice seleccionado. Esto afectaba, sobre todo, píxeles de ambientes palustres y en planicies aluviales. Por lo general, estos píxeles nulos se encontraban flanqueados de píxeles con NDWI37>=0. Así, se tacleó el problema asignando a cada 'píxel problema' el valor de la media local luego de aplicar un filtro móvil (ventana) de tamaño 7x7 (píxeles). En segundo lugar, detectamos valores 'falsos positivos' en algunas áreas palustres (e.g. 'Pampa Arenosa'). El efecto de no corregir esto conduce a un incremento en el sesgo negativo de los estimados del valor de área saturada (tanto como su distribució espacial), por aumento en el error de omisión, haciendo uso de la técnica de umbral NDWI fijado en 0. En consecuencia, para los 'falsos positivos' se optó por un reemplazo por los valores correspondientes al ndwi45 (combinación bandas 4 y 5), filtrados mediante una ventana de media movil, tamaño 7x7 (píxeles), considerando los píxeles NDWI37 'vecinos' no nulos). Los primeros experimetos mostraron una mejora en los estimados, permitiendo incorporar al análisis píxelescuyos valres debieran ser positivos, al menos bajo el supuesto de una alta autcorrelación espacial.
La capa resultante se denomina PROD_ndwi37c_YYYYMMDD y se presenta en formato Geotiff, con rango de valores válido entre 0 - 2000, siendo 9999 el valor para los datos declarados finalmente nulos (i.e. cobertura nubosa , grandes centros urbanos, nulos originales que no se ha podido rellenar). Para transformar los datos al dominio original del índice [-1,1] se debe realizar:
X/1000-1=Y
Siendo 0<=X<=2000 el valor de la capa PROD_ndwi37_YYYYMMDD (UInt16) e Y su valor correspondiente a PROD_ndwi37_YYYYMMDD.La capa refiere al método de relleno utilizado para reemplazar valores nulos 'originales' de la capa PROD_ndwi37_YYYYMMDD correspondiente al producto PROD, fecha YYYMMDD, y generar PROD_ndwi37c_YYYYMMDD.
La capa resultante se denomina CDP_MODIS_metodo37_YYYMMDD y se presenta en formato Geotiff, con rango de valores válido entre 0 y 4 correspondiente a la variable nominal 'método de obtención valor ndwi'. Luego, han de considerarse las siguientes etiquetas para cada valor:
0 = Dato correspondiente a capa ndwi37 (sin tratamiento)
1 = Dato nulo en capa ndwi37, reemplazado por media local obtenida sobre la base de aplicación de filro 7 x 7
2 = Falso positivo en capa ndwi37, reemplazado por media local combinndo ndwi45 (píxel problema) y ndwi37 (píxeles vecinos) obtenida sobre la base de aplicación de filro 7 x 7.
3 = Nubes y grandes centros urbanos. Excluidos del análisis
4 = Datos nulos sin rellenar
Los productos generados son de carácter operativo. Esto es, se obtienen a partir de la implementación de rutinas de análisis automatizadas sobre la base de la captura sistemática de datos provistos por LPDAAAC-USGS, EEUU (lpdaac.usgs.gov). Luego, las aplicaciones y análisis que pudieran desarrollarse sobre estos están siempre sujetas a validación por parte de los usuarios, tomando en responsabilidad propia las conclusiones derivadas. El INA no se responsabiliza del uso indebido de los productos ofrecidos.
Fecha de publicación: 2017-01-02
Contacto: Leandro Giordano (lgiordano@ina.gob.ar)