USGS Landsat 9 Level 2, Collection 2, Tier 1

LANDSAT/LC09/C02/T1_L2
Доступность набора данных
2021-10-31T00:00:00Z–2026-10-11T23:19:19.950000Z
Поставщик набора данных
Фрагмент кода 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. Поэтому ночные снимки не могут быть обработаны для определения температуры поверхности.

  • Известная ошибка в определении температуры поверхности относительно облаков и, возможно, их теней. These issues have been documented by Cook et al., (2014).

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

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

  • Геологическая служба США дополнительно проверяет нефизические значения (излучательная способность < 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 (прогнозируемый). Если поле отсутствует, пользователь должен выбрать 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 к центру изображения в пределах прямой видимости.

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" – окончательные расчетные значения кодировщика, сгенерированные после события переключения.

  • Значение ACTUAL указывает на фактические значения кодера.

TIRS_SSM_POSITION_STATUS STRING
  • NOMINAL – SSM работал нормально.

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

  • ESTIMATED – положение SSM было оценено, что может быть не так точно, как при статусе NOMINAL.

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

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

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

WRS_PATH ПЕР

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

WRS_ROW ПЕР

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

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

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

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

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

(Product, Image, Photograph, or Dataset Name) courtesy of the U.S. Geological Survey

Пример: "Изображение 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
Открыть в редакторе кода