USGS Landsat 9 Level 2, Collection 2, Tier 1

LANDSAT/LC09/C02/T1_L2
資料集可用性
2021-10-31T00:00:00Z–2025-08-30T23:07:34.096000Z
資料集供應來源
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 感應器產生的資料推導出的地面溫度。這些影像包含 5 個可見光和近紅外線 (VNIR) 波段,以及 2 個經過正射校正的短波紅外線 (SWIR) 波段,還有 1 個經過正射校正的熱紅外線 (TIR) 波段。這些檔案也包含用於計算 ST 產品的中間頻帶,以及 QA 頻帶。

Landsat 9 SR 產品是使用地表反射率程式碼 (LaSRC) 建立。所有 Collection 2 ST 產品都是使用單一管道演算法建立,該演算法是由羅徹斯特理工學院 (RIT) 和美國國家航空暨太空總署 (NASA) 噴射推進實驗室 (JPL) 共同建立。

收集到的資料條會封裝成重疊的「場景」,使用標準化參考格線,涵蓋約 170 公里 x 183 公里的範圍。

部分資產只有 SR 資料,因此會顯示 ST 頻段,但內容為空白。 如果資產同時具有 ST 和 SR 頻帶,「PROCESSING_LEVEL」會設為「L2SP」。如果素材資源只有 SR 波段,「PROCESSING_LEVEL」會設為「L2SR」。

其他說明文件和使用範例。

Landsat Collection 2 檔案已公開發布在 Google Cloud Storage bucket,使用者可依要求者付費模式存取。這些檔案會編入定期更新的 BigQuery 資料表,方便您進行分析:earth-engine-public-data.geo_index.landsat_c2_index

資料供應商注意事項:

  • 資料產品必須包含光學和熱能資料,才能順利處理為地表溫度,因為系統需要 ASTER NDVI 才能暫時將 ASTER GED 產品調整為目標 Landsat 場景。因此,系統無法處理夜間的溫度資料,並顯示溫度。

  • 在雲層和可能的雲影方面,地表溫度擷取作業存在已知錯誤。Cook 等人(2014)

  • ASTER GED 包含缺少平均放射率資料的區域,而這些資料是成功生成 ST 產品的必要條件。如果缺少 ASTER GED 資訊,這些區域就會缺少 ST 資料。

  • ASTER GED 資料集是根據 2000 年至 2008 年間取得的 ASTER 場景,以及所有晴空像素所建立。雖然這個資料集涵蓋全球空間範圍,但由於 ASTER 測量結果持續受到雲層汙染,部分區域缺少平均放射率資訊。

  • USGS 會進一步篩選 ASTER GED 中的非實體值 (放射率 < 0.6),移除因未偵測到雲層而低估的放射率。如果任何像素沒有 ASTER GED 輸入或非實體放射率值,產生的 Landsat ST 產品就會缺少像素。由於 ASTER GED 平均氣候資料的靜態性質,缺少 Landsat ST 像素的情況在 1982 年至今的期間內會保持一致。詳情請參閱 landsat-collection-2-surface-temperature-data-gaps-due-missing

頻帶

Pixel 大小
30 公尺

頻帶

名稱 單位 最小值 最大值 規模 偏移 像素大小 波長 說明
SR_B1 1 65455 2.75e-05 -0.2 公尺 0.435 至 0.451 微米

波段 1 (超藍色、沿海氣膠) 地表反射率

SR_B2 1 65455 2.75e-05 -0.2 公尺 0.452 至 0.512 微米

Band 2 (藍色) 表面反射率

SR_B3 1 65455 2.75e-05 -0.2 公尺 0.533 到 0.590 μm

第 3 波段 (綠色) 表面反射率

SR_B4 1 65455 2.75e-05 -0.2 公尺 0.636-0.673 μm

第 4 波段 (紅色) 表面反射率

SR_B5 1 65455 2.75e-05 -0.2 公尺 0.851 至 0.879 微米

第 5 波段 (近紅外線) 地表反射率

SR_B6 1 65455 2.75e-05 -0.2 公尺 1.566 至 1.651 微米

波段 6 (短波紅外線 1) 地表反射率

SR_B7 1 65455 2.75e-05 -0.2 公尺 2.107 至 2.294 微米

第 7 波段 (短波紅外線 2) 地表反射率

SR_QA_AEROSOL 公尺

噴霧屬性

ST_B10 K 0 65535 0.00341802 149 公尺 10.60-11.19 μm

頻帶 10 表面溫度。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_ATRAN 0 10000 0.0001 公尺

大氣透射率。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_CDIST km 0 24000 0.01 公尺

像素與雲端的距離。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_DRAD W/(m^2*sr*um)/ DN 0 28000 0.001 公尺

下行輻射。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_EMIS 0 10000 0.0001 公尺

根據 ASTER GED 估算波段 10 的放射率。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_EMSD 0 10000 0.0001 公尺

發射率標準差。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_QA K 0 32767 0.01 公尺

地表溫度頻帶的不確定度。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_TRAD W/(m^2*sr*um)/ DN 0 22000 0.001 公尺

熱感應頻帶轉換為輻射。 如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

ST_URAD W/(m^2*sr*um)/ DN 0 28000 0.001 公尺

上湧輻射。如果「PROCESSING_LEVEL」設為「L2SR」,這個波段會完全遮蓋。

QA_PIXEL 公尺

由 CFMASK 演算法產生的像素品質屬性。

QA_RADSAT 公尺

放射飽和度品質保證

圖片屬性

圖片屬性

名稱 類型 說明
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

用於通知使用者所用軌道星曆類型的 ID:「DEFINITIVE」或「PREDICTIVE」。如果沒有這個欄位,使用者應假設為「PREDICTIVE」。

GEOMETRIC_RMSE_MODEL DOUBLE

幾何殘差的綜合均方根誤差 (RMSE),以公尺為單位,包括沿軌和跨軌方向。只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

GEOMETRIC_RMSE_MODEL_X DOUBLE

在用於幾何精確度校正的地面控制點 (GCP) 上,以幾何殘差 (公尺) 測量的均方根誤差 (RMSE),方向為跨軌方向。只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

GEOMETRIC_RMSE_MODEL_Y DOUBLE

在沿軌方向的幾何精度校正中,以 GCP (地面控制點) 測量的幾何殘差均方根誤差 (公尺)。只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

GEOMETRIC_RMSE_VERIFY DOUBLE

使用 GLS2000 獨立測量,並在經過地形校正的產品上測得的幾何殘差 RMSE (公尺)。只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

GROUND_CONTROL_POINTS_MODEL DOUBLE

精確度校正程序中使用的 GCP 數量。 只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

GROUND_CONTROL_POINTS_VERIFY DOUBLE

用於驗證地形校正產品的地面控制點 (GCP) 數量。

GROUND_CONTROL_POINTS_VERSION DOUBLE

用於精確度修正程序的 GCP 資料集版本。只有在 L1_PROCESSING_LEVEL 為 L1TP 時,才會出現這個參數。

IMAGE_QUALITY_OLI INT

OLI 波段的影像品質。 1 = 最差,9 = 最佳,0 = 未計算品質。 如果是 Landsat 8,這個參數會針對使用 OLI 感應器較低 12 位元收集的場景向下調整 (TRUNCATION_OLI =「LOWER」)。

IMAGE_QUALITY_TIRS INT

TIRS 波段的影像品質。 1 = 最差,9 = 最佳,0 = 未計算品質。 如果以「SWITCHED」處理 TIRS_SSM_POSITION_STATUS 值,系統也會調降該值。

L1_DATE_PRODUCT_GENERATED STRING

對應 L1 產品的產品生成日期。

L1_LANDSAT_PRODUCT_ID STRING

對應 L1 產品的 Landsat 產品 ID。

L1_PROCESSING_LEVEL STRING

對應 L1 產品的處理層級。

L1_PROCESSING_SOFTWARE_VERSION STRING

對應 L1 產品的處理軟體版本。

LANDSAT_PRODUCT_ID STRING

Landsat 產品 ID

LANDSAT_SCENE_ID STRING

Landsat 場景 ID 簡稱

PROCESSING_LEVEL STRING

如果同時存在 SR 和 LST 頻帶,則為「L2SP」;如果只有 SR 頻帶,則為「L2SR」。

PROCESSING_SOFTWARE_VERSION STRING

建立產品的處理軟體版本。

ROLL_ANGLE DOUBLE

場景中心的太空船滾動角度。側傾值以偏擺轉向架構 (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 INT

最接近影像視線 (LOS) 場景中心的 WRS-2 路徑。

TARGET_WRS_ROW INT

最接近影像 LOS 場景中心的 WRS-2 列。第 880 到 889 列保留給北極,第 990 到 999 列保留給南極,這些區域未定義 WRS-2。

TEMPERATURE_MAXIMUM_BAND_ST_B10 DOUBLE

Band 10 可達到的最高溫度值。

TEMPERATURE_MINIMUM_BAND_ST_B10 DOUBLE

第 10 級錶帶可達到的最低溫度值。

TIRS_SSM_MODEL STRING

指出如何判斷 Landsat 8 TIRS 場景選取鏡 (SSM) 位置。

  • 「FINAL」表示在切換事件後產生的最終預估編碼器值。

  • 「ACTUAL」表示實際編碼器值。

TIRS_SSM_POSITION_STATUS STRING
  • 「NOMINAL」表示 SSM 在這個場景中正常運作。

  • 「SWITCHED」表示 SSM 在場景中切換了運作模式,可能導致 TIRS 影像品質問題,直接影響 IMAGE_QUALITY_TIRS 值。

  • 「預估」表示系統預估的 SSM 位置,可能不如「名義」狀態準確。

TRUNCATION_OLI STRING
  • 「LOWER」表示使用了較低的 12 位元。

  • 「UPPER」表示使用了上方 12 位元。

一般截斷模式為「UPPER」。

WRS_PATH INT

場景的 WRS 路徑編號。

WRS_ROW INT

場景的 WRS 列號。

使用條款

使用條款

Landsat 資料集是由聯邦政府建立的資料,因此屬於公共領域,可使用、轉移或複製,不受著作權限制。

請提供 USGS 做為資料來源的確認或出處,方法是加入一行文字引用,如下列範例所示。

(產品、圖片、相片或資料集名稱) 由美國地質調查局提供

示例:Landsat-7 影像由美國地質調查局提供

如要進一步瞭解如何正確引用及註明 USGS 產品出處,請參閱 USGS 外觀形象系統指南

使用 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 環境頁面,瞭解 Python API 和如何使用 geemap 進行互動式開發。

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
在程式碼編輯器中開啟