USGS Landsat 9 Level 2, Collection 2, Tier 1

LANDSAT/LC09/C02/T1_L2
Доступность наборов данных
2021-10-31T00:00:00Z–2026-10-09T23:34:52.889000Z
Поставщик набора данных
Фрагмент кода Earth Engine
ee.ImageCollection("LANDSAT/LC09/C02/T1_L2")
Интервал повторного посещения
16 дней
Теги
cfmask cloud fmask global l9sr landsat lasrc lc09 lst reflectance satellite-imagery sr usgs

Описание

Этот набор данных содержит скорректированные с учетом атмосферных условий значения отражательной способности поверхности и температуры поверхности суши, полученные на основе данных, полученных датчиками Landsat 9 OLI/TIRS. Эти изображения содержат пять видимых и ближних инфракрасных (VNIR) диапазонов и два коротковолновых инфракрасных (SWIR) диапазона, обработанных для ортотрансформированного отражения поверхности, и один тепловой инфракрасный (TIR) диапазон, обработанный для ортотрансформированной температуры поверхности. Они также содержат промежуточные каналы, используемые при расчете продуктов ST, а также каналы контроля качества.

Продукты Landsat 9 SR создаются с помощью кода отражения поверхности Земли (LaSRC). Все продукты Collection 2 ST созданы с помощью одноканального алгоритма, разработанного совместно Рочестерским технологическим институтом (RIT) и Лабораторией реактивного движения (JPL) Национального управления по аэронавтике и исследованию космического пространства (NASA).

Собранные данные упаковываются в перекрывающиеся "сцены" размером примерно 170 км x 183 км с использованием стандартизированной опорной сетки.

Некоторые объекты содержат только данные SR, и в этом случае полосы ST присутствуют, но пусты. Для объектов с диапазонами ST и SR свойство PROCESSING_LEVEL имеет значение L2SP. Для объектов, в которых есть только каналы SR, для параметра PROCESSING_LEVEL задано значение L2SR.

Дополнительная документация и примеры использования

Файлы Landsat Collection 2 общедоступны в сегменте Google Cloud Storage на условиях оплаты отправителем. Файлы индексируются в регулярно обновляемой таблице BigQuery, что упрощает анализ: earth-engine-public-data.geo_index.landsat_c2_index.

Примечания поставщика данных:

  • Чтобы успешно преобразовать данные в температуру поверхности, продукты данных должны содержать как оптические, так и тепловые данные, поскольку ASTER NDVI требуется для временной корректировки продукта ASTER GED в целевой сцене Landsat. Поэтому ночные снимки не могут быть обработаны для определения температуры поверхности.

  • Известна ошибка в извлечении температуры поверхности, связанная с облаками и, возможно, тенями от облаков. Описание этих проблем приведено в работе Cook et al., (2014).

  • ASTER GED содержит области, в которых отсутствуют средние данные об излучательной способности, необходимые для успешного создания продуктов ST. Если в файле ASTER GED отсутствует информация, в этих областях не будет данных ST.

  • Набор данных ASTER GED создан на основе всех пикселей сцен ASTER, полученных с 2000 по 2008 год при безоблачном небе. Хотя этот набор данных имеет глобальный пространственный охват, в нем отсутствуют данные о средней излучательной способности для некоторых областей из-за постоянного загрязнения облаками измерений ASTER.

  • Специалисты USGS дополнительно проверяют нефизические значения (коэффициент излучения < 0,6) в ASTER GED, чтобы удалить любые заниженные значения коэффициента излучения, вызванные незамеченными облаками. Для любого заданного пикселя без входных данных ASTER GED или нефизического значения излучательной способности в результирующих продуктах Landsat ST отсутствуют пиксели. Отсутствующие пиксели Landsat ST будут согласованными во времени (с 1982 года по настоящее время) с учетом статического характера средних климатологических данных ASTER GED. Подробнее landsat-collection-2-surface-temperature-data-gaps-due-missing…

Диапазоны частот

Диапазоны частот

Размер пикселя: 30 метров (все диапазоны)

Название Единицы Мин. Макс. Масштаб Смещение Размер пикселя Длина волны Описание
SR_B1 1 65455 2,75e-05 -0,2 30 метров 0,435–0,451 мкм

Отражательная способность поверхности в диапазоне 1 (ультрасиний, прибрежный аэрозоль)

SR_B2 1 65455 2,75e-05 -0,2 30 метров 0,452–0,512 мкм

Отражательная способность поверхности в диапазоне 2 (синий)

SR_B3 1 65455 2,75e-05 -0,2 30 метров 0,533–0,590 мкм

Поверхностное отражение в диапазоне 3 (зеленый)

SR_B4 1 65455 2,75e-05 -0,2 30 метров 0,636–0,673 мкм

Отражательная способность поверхности в диапазоне 4 (красный)

SR_B5 1 65455 2,75e-05 -0,2 30 метров 0,851–0,879 мкм

Поверхностное отражение в диапазоне 5 (ближний инфракрасный)

SR_B6 1 65455 2,75e-05 -0,2 30 метров 1,566–1,651 мкм

Отражательная способность поверхности в диапазоне 6 (коротковолновый инфракрасный 1)

SR_B7 1 65455 2,75e-05 -0,2 30 метров 2,107–2,294 мкм

Отражательная способность поверхности в диапазоне 7 (коротковолновый инфракрасный 2)

SR_QA_AEROSOL 30 метров Нет

Атрибуты аэрозолей

ST_B10 K 0 65535 0,00341802 149 30 метров 10,60–11,19 мкм

Температура поверхности (канал 10). Если для параметра PROCESSING_LEVEL задано значение L2SR, этот диапазон полностью маскируется.

ST_ATRAN 0 10000 0,0001 30 метров Нет

Пропускание атмосферы. Если для параметра PROCESSING_LEVEL задано значение L2SR, этот диапазон полностью маскируется.

ST_CDIST km 0 24000 0,01 30 метров Нет

Расстояние в пикселях до облака. Если для параметра PROCESSING_LEVEL задано значение L2SR, этот диапазон полностью маскируется.

ST_DRAD Вт/(м^2*ср*мкм)/ DN 0 28000 0,001 30 метров Нет

Излучение, направленное вниз. Если для PROCESSING_LEVEL задано значение L2SR, эта полоса полностью замаскирована.

ST_EMIS 0 10000 0,0001 30 метров Нет

Коэффициент излучения канала 10, рассчитанный на основе данных ASTER GED. Если для PROCESSING_LEVEL задано значение L2SR, эта полоса полностью замаскирована.

ST_EMSD 0 10000 0,0001 30 метров Нет

Стандартное отклонение излучательной способности. Если для PROCESSING_LEVEL задано значение L2SR, эта полоса полностью замаскирована.

ST_QA K 0 32767 0,01 30 метров Нет

Неопределенность диапазона температуры поверхности. Если для параметра PROCESSING_LEVEL задано значение L2SR, этот диапазон полностью маскируется.

ST_TRAD Вт/(м^2*ср*мкм)/ DN 0 22000 0,001 30 метров Нет

Термальный диапазон преобразован в излучение. Если для параметра PROCESSING_LEVEL задано значение L2SR, этот диапазон полностью маскируется.

ST_URAD Вт/(м^2*ср*мкм)/ DN 0 28000 0,001 30 метров Нет

Исходящее излучение. Если для PROCESSING_LEVEL задано значение L2SR, эта полоса полностью замаскирована.

QA_PIXEL 30 метров Нет

Атрибуты качества пикселей, сгенерированные на основе алгоритма CFMASK.

QA_RADSAT 30 метров Нет

Контроль качества радиометрического насыщения

Свойства изображения

Свойства изображения

Название Тип Описание
ALGORITHM_SOURCE_SURFACE_REFLECTANCE STRING

Название и версия алгоритма отражения поверхности.

ALGORITHM_SOURCE_SURFACE_TEMPERATURE STRING

Название и версия алгоритма измерения температуры поверхности.

CLOUD_COVER DOUBLE

Процент облачности (0–100), -1 = не рассчитано.

CLOUD_COVER_LAND DOUBLE

Процент облачного покрова над сушей (0–100), -1 = не вычисляется.

COLLECTION_CATEGORY STRING

Категория подборки сцен: "T1" или "T2".

DATA_SOURCE_AIR_TEMPERATURE STRING

Источник данных о температуре воздуха.

DATA_SOURCE_ELEVATION STRING

Источник данных о высоте.

DATA_SOURCE_OZONE STRING

Источник данных об озоне.

DATA_SOURCE_PRESSURE STRING

Источник данных о давлении.

DATA_SOURCE_REANALYSIS STRING

Источник данных для повторного анализа.

DATA_SOURCE_TIRS_STRAY_LIGHT_CORRECTION STRING

Источник данных для коррекции паразитной засветки TIRS.

DATA_SOURCE_WATER_VAPOR STRING

Источник данных о водяном паре.

DATE_PRODUCT_GENERATED DOUBLE

Временная метка даты создания товара.

EARTH_SUN_DISTANCE DOUBLE

Расстояние от Земли до Солнца (а. е.).

EPHEMERIS_TYPE STRING

Идентификатор, который сообщает пользователю, какой тип эфемерид использовался: "DEFINITIVE" (окончательный) или "PREDICTIVE" (прогнозный). Если поле отсутствует, пользователь должен предположить, что используется режим "Прогнозирование".

GEOMETRIC_RMSE_MODEL DOUBLE

Объединенная среднеквадратическая ошибка геометрических невязок (в метрах) в поперечном и продольном направлениях. Этот параметр присутствует, только если L1_PROCESSING_LEVEL имеет значение L1TP.

GEOMETRIC_RMSE_MODEL_X DOUBLE

Среднеквадратическая ошибка геометрических невязок (в метрах), измеренная на опорных точках, которые использовались для коррекции геометрической точности в направлении, перпендикулярном траектории полета. Этот параметр присутствует, только если для параметра L1_PROCESSING_LEVEL задано значение L1TP.

GEOMETRIC_RMSE_MODEL_Y DOUBLE

Среднеквадратическая ошибка геометрических невязок (в метрах), измеренная на опорных точках, которые использовались для коррекции геометрической точности в направлении вдоль траектории. Этот параметр присутствует, только если L1_PROCESSING_LEVEL имеет значение L1TP.

GEOMETRIC_RMSE_VERIFY DOUBLE

Среднеквадратическая ошибка геометрических невязок (в метрах), измеренная на продукте с коррекцией рельефа независимо с помощью GLS2000. Этот параметр присутствует только в том случае, если для параметра L1_PROCESSING_LEVEL задано значение L1TP.

GROUND_CONTROL_POINTS_MODEL DOUBLE

Количество опорных точек, использованных при коррекции точности. Этот параметр присутствует только в том случае, если для параметра L1_PROCESSING_LEVEL задано значение L1TP.

GROUND_CONTROL_POINTS_VERIFY DOUBLE

Количество опорных точек, использованных при проверке продукта с коррекцией рельефа.

GROUND_CONTROL_POINTS_VERSION DOUBLE

Версия набора данных GCP, использованная в процессе корректировки точности. Этот параметр присутствует только в том случае, если для L1_PROCESSING_LEVEL задано значение L1TP.

IMAGE_QUALITY_OLI ПЕР

Качество изображений, полученных с помощью каналов OLI. 1 – наихудшее качество, 9 – наилучшее, 0 – качество не рассчитано. Для Landsat 8 этот параметр уменьшается для сцен, собранных с использованием нижних 12 бит от датчика OLI (TRUNCATION_OLI = "LOWER").

IMAGE_QUALITY_TIRS ПЕР

Качество изображений в диапазонах TIRS. 1 – наихудшее, 9 – наилучшее, 0 – качество не рассчитано. Она также корректируется в меньшую сторону для сцен, обработанных со значением SWITCHED для TIRS_SSM_POSITION_STATUS.

L1_DATE_PRODUCT_GENERATED STRING

Дата создания товара, соответствующего товару уровня 1.

L1_LANDSAT_PRODUCT_ID STRING

Идентификатор продукта Landsat для соответствующего продукта уровня 1.

L1_PROCESSING_LEVEL STRING

Уровень обработки для соответствующего продукта L1.

L1_PROCESSING_SOFTWARE_VERSION STRING

Версия ПО для обработки соответствующего продукта L1.

LANDSAT_PRODUCT_ID STRING

Идентификатор продукта Landsat

LANDSAT_SCENE_ID STRING

Короткий идентификатор сцены Landsat

PROCESSING_LEVEL STRING

"L2SP", если присутствуют диапазоны SR и LST, или "L2SR", если присутствуют только диапазоны SR.

PROCESSING_SOFTWARE_VERSION STRING

Версия ПО, в которой был создан продукт.

ROLL_ANGLE DOUBLE

Угол крена космического аппарата в центре сцены. Значение крена задается в системе координат Yaw Steering Frame (YSF), ось X которой совпадает с вектором мгновенной скорости движения по земле. Крен вокруг оси X определяется по правилу правой руки: положительный крен приводит к тому, что инструменты указывают влево от наземного трека, а отрицательный – вправо.

SCENE_CENTER_TIME STRING

Время наблюдения в виде строки ISO 8601.

SENSOR_ID STRING

Название датчика.

SPACECRAFT_ID STRING

Название космического корабля.

SUN_AZIMUTH DOUBLE

Азимут солнца в градусах для центра изображения в момент съемки. Положительное значение указывает на угол к востоку или по часовой стрелке от севера. Отрицательное значение указывает на угол к западу или против часовой стрелки от севера.

SUN_ELEVATION DOUBLE

Угол возвышения солнца в градусах для центральной точки изображения в момент съемки. Положительное значение указывает на сцену, снятую днем. Отрицательное значение указывает на ночную сцену. Примечание. Для расчета отражательной способности требуется зенитный угол солнца, который равен 90 – угол возвышения солнца.

TARGET_WRS_PATH ПЕР

Ближайший путь WRS-2 к центру сцены изображения в пределах прямой видимости (LOS).

TARGET_WRS_ROW ПЕР

Ближайшая строка WRS-2 к центру сцены LOS изображения. Строки 880–889 зарезервированы для Северного полюса, а 990–999 – для Южного полюса, где система WRS-2 не определена.

TEMPERATURE_MAXIMUM_BAND_ST_B10 DOUBLE

Максимально достижимое значение температуры для диапазона 10.

TEMPERATURE_MINIMUM_BAND_ST_B10 DOUBLE

Минимальное достижимое значение температуры для диапазона 10.

TIRS_SSM_MODEL STRING

Указывает, как было определено положение зеркала выбора сцены (SSM) Landsat 8 TIRS.

  • "FINAL" – окончательные расчетные значения кодировщика, сгенерированные после события переключения.

  • "ФАКТИЧЕСКИЕ" – фактические значения кодировщика.

TIRS_SSM_POSITION_STATUS STRING
  • "NOMINAL" означает, что SSM работал нормально для этой сцены.

  • "SWITCHED" означает, что SSM переключил режим работы в сцене и может иметь проблемы с качеством изображения TIRS, которые напрямую влияют на значение IMAGE_QUALITY_TIRS.

  • "ОРИЕНТИРОВОЧНОЕ" означает, что положение SSM было определено приблизительно и может быть не таким точным, как при статусе "НОМИНАЛЬНОЕ".

TRUNCATION_OLI STRING
  • "LOWER" – использовались нижние 12 бит.

  • "UPPER" означает, что использовались 12 старших битов.

По умолчанию используется режим UPPER.

WRS_PATH ПЕР

Номер пути сцены в системе WRS.

WRS_ROW ПЕР

Номер строки сцены в WRS.

Условия использования

Условия использования

Наборы данных Landsat созданы федеральным правительством, поэтому относятся к общественному достоянию и могут использоваться, передаваться и воспроизводиться без ограничений авторского права.

При использовании данных USGS необходимо указать источник, добавив строку с цитатой, как в примере ниже.

(название продукта, изображения, фотографии или набора данных), предоставлено Геологической службой США

Пример: "Изображение Landsat-7 предоставлено Геологической службой США".

Подробную информацию о том, как правильно ссылаться на продукты Геологической службы США и упоминать их, можно найти в Руководстве по фирменному стилю Геологической службы США.

Исследуйте данные с помощью Earth Engine

Редактор кода (JavaScript)

var dataset = ee.ImageCollection('LANDSAT/LC09/C02/T1_L2')
    .filterDate('2022-01-01', '2022-02-01');

// Applies scaling factors.
function applyScaleFactors(image) {
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  var thermalBands = image.select('ST_B.*').multiply(0.00341802).add(149.0);
  return image.addBands(opticalBands, null, true)
              .addBands(thermalBands, null, true);
}

dataset = dataset.map(applyScaleFactors);

var visualization = {
  bands: ['SR_B4', 'SR_B3', 'SR_B2'],
  min: 0.0,
  max: 0.3,
};

Map.setCenter(-114.2579, 38.9275, 8);

Map.addLayer(dataset, visualization, 'True Color (432)');

Настройка Python

Информацию о Python API и использовании geemap для интерактивной разработки можно найти на странице Среда Python.

import ee
import geemap.core as geemap

Colab (Python)

dataset = ee.ImageCollection('LANDSAT/LC09/C02/T1_L2').filterDate(
    '2022-01-01', '2022-02-01'
)


# Applies scaling factors.
def apply_scale_factors(image):
  optical_bands = image.select('SR_B.').multiply(0.0000275).add(-0.2)
  thermal_bands = image.select('ST_B.*').multiply(0.00341802).add(149.0)
  return image.addBands(optical_bands, None, True).addBands(
      thermal_bands, None, True
  )


dataset = dataset.map(apply_scale_factors)

visualization = {
    'bands': ['SR_B4', 'SR_B3', 'SR_B2'],
    'min': 0.0,
    'max': 0.3,
}

m = geemap.Map()
m.set_center(-114.2579, 38.9275, 8)
m.add_layer(dataset, visualization, 'True Color (432)')
m
Открыть в редакторе кода