Abrir EarthExplorer

Ir a https://earthexplorer.usgs.gov/. Es el portal oficial del USGS para descarga manual de imágenes. Para descargas programáticas masivas existe la USGS Machine-to-Machine API; para visualización rápida sin descarga, LandsatLook.

Registro e inicio de sesión

Hacer clic en Register. La cuenta es gratuita pero obligatoria para descargar. Una vez creada, iniciar sesión en el portal.

Definir el área de búsqueda

EarthExplorer ofrece varias opciones para acotar la búsqueda: (1) buscar por nombre de lugar, (2) hacer zoom en el mapa interactivo, (3) cargar un archivo KML o un Shapefile comprimido como ZIP.

Recomendado: cargar un Shapefile con el polígono o caja delimitadora del AOI. Antes de cargarlo, comprimir todos los archivos del shapefile (.shp, .shx, .dbf, .prj) en un único .zip.

Selección de Shapefile y carga del ZIP en EarthExplorer.

Selección del dataset

Pasar a la pestaña Data Sets. Navegar a Landsat → Landsat Collection 2 → Landsat Collection 2 Level-2.

Aquí es importante elegir según fechas de lanzamiento y tiempo en órbita:

Para series 2013–presente conviene combinar Landsat 8 y 9 (frecuencia efectiva de 8 días).

Filtros de búsqueda

Antes de pasar a resultados, aplicar filtros críticos:

• Rango de fechas: ajustar según la ventana de interés (estación seca, antes/después de un evento, fenología de cultivo).
• Cloud Cover (Land): usar Cloud Cover Land y no Cloud Cover Scene. El segundo incluye nubes sobre océano que no afectan tu estudio.
• Day/Night: para Landsat 8/9 las pasadas son diurnas; mantener por defecto.

Umbral típico: nubosidad < 10 % para análisis directo; < 30 % si se va a aplicar máscara de nubes con QA_PIXEL.

Elegir la imagen y descargar

En Results aparecen las escenas que cumplen los filtros. Revisar la previsualización antes de descargar para confirmar que el AOI tenga buena visibilidad (sin nubes, sin sombras de cumulonimbo, sin humo). Iniciar la descarga con el ícono de flecha verde.

Productos disponibles para descarga

El cuadro de descarga ofrece varios productos. Los dos más usados:

• Landsat Collection 2 Level-2 Surface Reflectance Bundle: bandas SR + QA + metadatos (~700 MB).
• Landsat Collection 2 Level-2 Surface Temperature Bundle: banda ST + auxiliares térmicas + QA (~300 MB).

Si solo se trabaja con índices ópticos, el bundle SR basta. Si el estudio incluye temperatura, balance energético o detección de calor urbano, descargar también el bundle ST.

Cargar las bandas en QGIS

Abrir QGIS LTR (versión Long-Term Release) y arrastrar todas las bandas multiespectrales y la térmica al panel de capas. Para Landsat 8/9 son:

Aplicar escala a las bandas multiespectrales

Usar la Calculadora Raster con la fórmula:

"SR_B4@1" * 0.0000275 + (-0.2)

Repetir para cada banda (B1 a B7). Para automatizar, conviene usar el algoritmo Raster calculator (GDAL) en lote desde el procesado, o un script PyQGIS:

bands = ['SR_B1','SR_B2','SR_B3','SR_B4','SR_B5','SR_B6','SR_B7']
for b in bands:
    processing.run("gdal:rastercalculator", {
        'INPUT_A': f'{path}/{escena}_{b}.TIF',
        'BAND_A': 1,
        'FORMULA': 'A*0.0000275-0.2',
        'OUTPUT': f'{out}/{escena}_{b}_SR.tif'
    })

Banda térmica: escala y conversión a Celsius

Para la banda térmica ST_B10, aplicar el factor específico y luego restar 273.15 si se desea Celsius:

("ST_B10@1" * 0.00341802 + 149.0) - 273.15

Ojo: el resultado solo tiene sentido sobre superficies opacas (suelo, vegetación, agua). Sobre nubes o píxeles saturados los valores son irreales (50–80 °C en zonas frías), por eso el siguiente paso lógico es enmascarar con QA_PIXEL.

Composición multibanda

Una vez escaladas las bandas SR, generar un raster multibanda con todas. En QGIS la herramienta más rápida es Construir raster virtual (gdalbuildvrt) con las 7 bandas SR escaladas, marcando la opción Place each input file into a separate band.

Reproyección al UTM de la zona de estudio

Las escenas Landsat C2L2 vienen ya en proyección UTM con datum WGS84, pero la zona corresponde a la del centro de la escena. Si tu AOI está en otra zona UTM (caso típico en zonas fronterizas entre dos husos), reproyectar es obligatorio.

Herramienta: Warp (reproject) de GDAL.

Método de remuestreo según el tipo de banda:

Usar nearest neighbor sobre bandas de calidad evita inventar valores QA inexistentes; usarlo sobre SR introduce escalonado visible.

Recorte por AOI

Último paso del flujo tradicional: recortar el raster al polígono del área de estudio.

Herramienta: Clip raster by mask layer. Configurar:

• Capa de máscara: shapefile del AOI.
• Marcar Match the extent of the clipped raster to the extent of the mask layer para que el raster de salida tenga el bounding box del AOI.
• Activar Keep resolution of input raster para no degradar el píxel.
• Compresión LZW o DEFLATE en Profile para ahorrar disco.