13.  Convertimos nuestra lista de números a una cadena de caracteres con la función stream str, luego creamos un bucle que nos genera el valor de NDVI para cada año (desde 1999 hasta 2019).
Utilizamos  ee. ImageCollectionpara acceder a todas las series temporales relacionadas con Landsat 7. Más específicamente, utilizaremos el producto llamado LANDSAT / LE07 / C01 / T1_SR para calcular NDVI.
Todas las imágenes existentes en los conjuntos de datos mencionados anteriormente se almacenarán en la variable llamada 17collection.
También filtramos todas las imágenes por área (usando filterBounds) y por año (usando filterDate). Como creamos una imagen NDVI anual, cambiamos el año varias veces (21 veces, de 1999 a 2019).
Ejecutamos el composite que es una imagen  generada a partir de muchas otras imágenes y agrego el resultado a la lista que cree en el inicio bands=( ).    
La composición es un proceso muy común en la teledetección que se refiere al proceso de combinar imágenes superpuestas espacialmente en una sola basada en una función de agregación.
En este caso, crearemos una sola imagen que contenga el valor máximo de NDVI por píxel. Es decir. El valor máximo de NDVI en el píxel A podría ocurrir en la fecha x, mientras que ese valor máximo en el píxel B podría ocurrir en la fecha. Esto significa que la fecha de la imagen sería un compuesto de todas las imágenes combinadas, Esto es posible gracias a la función  qualityMosaic.