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] (developers.google.com/earth-engine/datasets..) 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.

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.

1. 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. 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
   

1. 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);
   

2. 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);
   

3. 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
   

4. Resultado En la siguiente Figura se muestra el resultado de la implementación del código.

Revisa el código: 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
})