Cómo acceder al catálogo de imágenes de GEE con información de imágenes ya procesadas anteriormente

Hace algún tiempo bajé el set completo de unas imágenes de Landsat 5 y por daño irreversible en el disco duro perdí toda la información allí disponible. Lo que pude rescatar fue lo que ya había almacenado en la nube (mi cuenta en dropbox). Entre ellas, unas imágenes ya cortadas con una máscara vectorial pero para las cuales preservé el código completo del nombre asignado a ellas en las imágenes originales. La banda 3 correspondía a la siguiente: LT50380322011235PAC01_B3.

El código anterior se refiere a una imágen Lansat 5 para la cual recuerdo que en el nombre se debía especificar el path/row y la fecha de adquisición correspondiente por parte del satélite. En la imagen siguiente se muestra la máscara vectorial que usé en esa oportunidad y que corresponde al resultado del último post. Por tanto, en la zona oscura de la API de GEE marqué un punto arbitrario con el cursor y copié las coordenadas desplegadas en la pestaña del Inspector.

gee3

Las coordenadas obtenidas en el Inspector fueron usadas en la web siguiente:

https://landsat.usgs.gov/wrs-2-pathrow-latitudelongitude-converter

para obtener el path/row respectivo; tal como se visualiza en la imagen siguiente:

gee4

Por tanto, los 6 dígitos 038032 a continuación de LT5 hacen referencia al path/row de la imagen. La fecha está en una notación, 2011235, donde los primeros cuatro dígitos corresponden al año 2011 y los tres siguientes al día juliano. Para convertir día juliano en día calendario se empleó el script python siguiente:

def leap_year(year):
    return  ( year % 4 == 0 and year % 100 !=0 ) or year % 400 == 0

year = int(raw_input("year = ? ")) 
leap = leap_year(year)
 
days = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31]
 
if leap == True:
    days[1] = 29
 
day = int(raw_input("julian day = ? ")) 
sum = 0
 
for k, item in enumerate(days):
    sum += item
     
    if sum >= day:
        dd = item - (sum - day)
         
        if k+1 < 10 or dd < 10:
            mm = str(k+1).zfill(2)
            dd = str(dd).zfill(2)
        else:
            mm = str(k+1)
            dd = str(dd)
 
        formatted_day = str(year) + "-" + mm + "-" +  dd
         
        break
 
print formatted_day

La ejecución del script anterior permite determinar que el día juliano 235 corresponde a la fecha '2011-08-23' que será usada en el método 'filterDate' de GEE.

Para ubicar la imageCollection en “Search places y datasets…” filtramos por Landsat y ubicamos la Landsat 5 TM collection; tal como se tiene en la imagen siguiente:

gee5

Corroboramos la metadata e importamos a GEE.

gee6

Completamos la variable importada con la información correspondiente al script a continuación y ejecutamos.

var imageCollection = ee.ImageCollection("LANDSAT/LT05/C01/T1")
  .filterDate('2011-08-23', '2011-08-24')
  .filterBounds(ee.Geometry.Point(-112.038, 40.1725));
  
print (imageCollection);

Map.addLayer(imageCollection,
             {bands: ['B1']},
             'images');

El resultado es el siguiente:

gee7

Se añadió a la vista sólo la banda 1 pero la imageCollection contiene 8 bandas. El ‘print’ en el código del script permite visualizar en la Console toda la metadata correspondiente a estas imágenes; tal como se presenta a continuación:

gee8

Si se usa el Inspector se puede observar que la ImageCollection LT05_038032_20110823 corresponde, tal como se esperaba, a una serie de sólo 8 imágenes adquiridas el día 2011-08-23 a las 17:56:58.

Anuncios
Esta entrada fue publicada en Google Earth Engine. Guarda el enlace permanente.

Una respuesta a Cómo acceder al catálogo de imágenes de GEE con información de imágenes ya procesadas anteriormente

  1. Pingback: Cómo acceder en GEE a imágenes individuales ya transformadas (radiancia, reflectancia, temperatura de brillo) empleando el código de nombre de las imágenes Landsat 5 | El Blog de José Guerrero

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión /  Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión /  Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión /  Cambiar )

Conectando a %s