# Evapotranspiración con datos MODIS desde Google Earth Engine

Es este ejemplo se muestra un script para obtener datos de ***Evapotranspiración (ET)*** donde se usa el conjunto de datos evapotranspiración/flujo de calor latente [MOD16A2 versión 6]
(https://developers.google.com/earth-engine/datasets/catalog/MODIS_006_MOD16A2) que es un producto compuesto de 8 días con una resolución espacial de 500 metros/pixel y contiene información desde el 2001.

![image.png](https://cdn.hashnode.com/res/hashnode/image/upload/v1662094460175/xXmdJ29mi.png align="left")

El algoritmo utilizado para la recopilación de productos de datos ***MOD16*** se basa en la lógica de la ecuación de Penman-Monteith, incluye entradas de datos de reanálisis diarios junto con productos de datos de detección remota MODIS, como la dinámica de las propiedades de la vegetación, el albedo y la cobertura terrestre. 

En los siguientes pasos vamos a mostrar como calcular la ET acumulada en 2021 para México.

1. Se importan algunos script de la herramientas VICAL, útil  para la visualización. 
```
var imp = require('users/InifapCenidRaspa/VICAL:Exportaciones');
var St = require('users/InifapCenidRaspa/VICAL:Style');
```

2. Se importa la colección de imagenes y, se define el intervalo (fecha inicial y final) y la banda a utilizar (**ET**). 
```
var collection = ee.ImageCollection("MODIS/006/MOD16A2")
  .filterDate('2021-01-01', '2021-12-31')
  .select('ET');
```
Las bandas con las que cuenta este Dataset se muestran en la siguiente imagen, de estas bandas, los valores de píxel para las dos capas de evapotranspiración (ET y PET) son la suma de los ocho días dentro del período compuesto.

![image.png](https://cdn.hashnode.com/res/hashnode/image/upload/v1662094970990/sgj4gLn-d.png align="left")


5.  Definimos la región, que este caso seria México, por tanto usamos el archivo disponible en GEE de *Polígonos de límite internacional*. Para cambiar el país se debe cambiar las dos letras **MX** por el código de país, para conocer el este código puede visitar la dirección [https://en.wikipedia.org/wiki/List_of_FIPS_country_codes](https://en.wikipedia.org/wiki/List_of_FIPS_country_codes). Por ejemplo para Perú seria [**PE**]
```
var region = ee.FeatureCollection("USDOS/LSIB_SIMPLE/2017")
              .filter(ee.Filter.eq('country_co', 'MX'));// 'MX' de México
```


5.  Calculamos la ET y asignamos el estadísticos, que en este caso seria **SUMA**. Posteriormente recortamos con respecto a nuestra región de interés.  
Después del paso anterior multiplicamos los valores obtenidos por 0.1 que es la escala con la que nos da GEE la ET de este DataSet.
```
var composite =collection.sum().clip(region);
var composite =composite.multiply(0.1);
```

5. Definimos la paleta de visualización y  la leyenda que se mostrara en el mapa
```
var paletaCl = {min: 0, max: 1300, palette:St.ETO};
var panL = imp.Leyenda('ET (mm)', paletaCl);
```
6. Finalmente, se pida a GEE que muestre el mapa con los valores de ET 
```
Map.centerObject(region)//centrar a la region
Map.addLayer(composite, paletaCl, ' Suma_ET');
Map.add(panL);// leyenda
```
7. Resultado 
En la siguiente Figura se muestra el resultado de la implementación del código. 
![image.png](https://cdn.hashnode.com/res/hashnode/image/upload/v1662096075631/mrKVGbSHa.png align="left")

Revisa el código: [https://code.earthengine.google.com/2fb626db113b2d2efad82d05060c1bef](https://code.earthengine.google.com/2fb626db113b2d2efad82d05060c1bef)

si quiere exportar el resultado en formato Raster puede escribir el siguiente codigo
```
Export.image.toDrive({
  image: composite,
  description: 'ET_Export',
  folder: 'earthengine',
  fileNamePrefix: 'ET',
  scale: 500,
  maxPixels: 1e10
})
```