El análisis de deforestación y talas ilegales en series temporales puede ser identificado rápidamente jugando con series multitemporales de deforestación en Google Earth Engine. Los recursos radar, con aportación diaria de datos, ofrece series temporales de datos que te pueden ayudar a comparar mosaicos temporales para identificar las zonas territoriales afectadas y las teselas forestales arrasadas de manera antrópica.
Los datos de Sentinel 1 son una interesante vía para los análisis de deforestación en Google Earth Engine sin necesidad de descargar, de forma individual, amplias series de datos diarios. Sólo necesitas entender un poco la dinámica de datos rádar y los filtros necesarios para rescatar las bandas de la colección radar de Sentinel. Puedes descargar como referencia este script de análisis multitemporal para localizar las masas forestales por medio de una simbología de colores.
Desde la colección de datos rádar de Sentinel 1 necesitarás hacer una selección de la banda en la que estés interesado. Deberás jugar con el insturmento IW (Interferometric Wide Swath) y acceder a las bandas VV, o VH, disponibles en la colección de imágenes rádar pudiendo emplear alguna de ellas. La banda VH te ayudará a identificar de manera más nítida presencia de masas de agua entre las zonas forestales. De igual forma, puedes jugar con una o dos de las órbitas (ascendente y descendente) del satélite y empezar a filtrar los datos de la colección que mostrarán la deforestación en Google Earth Engine.
Para identificar las tasas de cambio en las masas forestales puedes hacer comparativas entre dos momentos temporales jugando con efecto swipe o jugar con un sencillo análisis multitemporal de tres momentos clave. Filtrando los tres momentos temporales de análisis y componiendo el mosaico a través de la delimitación del ROI territorial podrás extraer los datos de tu zona de interés.
Una sencilla composición de las tres bandas temporales generará una imagen RGB que te ayudará a representar, a través de una rampa de colores, los lugares en los que hay presencia o ausencia de datos en cada momento temporal.
Para la interpretación de los datos puedes trabajar una sencilla leyenda sobre el visor, haciendo una clasificación de datos de presencia, ausencia o presencia parcial de masa forestal en cada momento temporal. Sólo necesitarás definir las etiquetas de la leyenda y asignarle la simbología tradicional generada en las imágenes multitemporales.
Adicionalmente puedes incorporar en el visor la representación de datos radar para cada uno de los momentos temporales y observar e interpretar las coberturas forestales a través de las bandas.
Una vez tengas montado tu script, puedes correrlo y visualizar los datos de deforestación en Google Earth Engine.
Si los datos rádar de Sentinel 1 no cubren tus momentos temporales de análisis, dispones de alternativas vinculadas con los análisis de deforestación en Google Earth Engine trabajando con las colecciones del Hansen Global Forest Change. En este caso dispones de un repertorio de datos, generados por Landsat 7, y basados en pérdidas de masa forestal a través de imágenes ópticas con las que analizar índices de vegetación como el NDVI y periodos históricos de 18 años evaluando ganancias y pérdidas de masas forestales con cobertura mundial.
Si quieres evaluar pérdidas de masas forestales a través de incendios en Earth Engine, también puedes jugar con los análisis multitemporales combinando índices de calcinación y severidad de incendios a través de este otro script de GEE.
Dispones de otros scripts temáticos que podrás descargar desde esta entrada.
