USGS Landsat 9 Level 2, Collection 2, Tier 1

LANDSAT/LC09/C02/T1_L2
Dostępność zbioru danych
2021-10-31T00:00:00Z–2025-08-30T23:07:34.096000Z
Dostawca zbioru danych
Fragment kodu Earth Engine
ee.ImageCollection("LANDSAT/LC09/C02/T1_L2")
Odstęp między kolejnymi wizytami
16 dni
Tagi
cfmask cloud fmask global l9sr landsat lasrc lc09 lst reflectance satellite-imagery sr usgs

Opis

Ten zbiór danych zawiera skorygowane atmosferycznie odbicie światła słonecznego od powierzchni i temperaturę powierzchni gruntu uzyskane na podstawie danych z czujników OLI/TIRS satelity Landsat 9. Obrazy te zawierają 5 pasm widzialnych i bliskiej podczerwieni (VNIR) oraz 2 pasma podczerwieni krótkofalowej (SWIR) przetworzone na ortorektyfikowane odbicie powierzchniowe, a także 1 pasmo podczerwieni termicznej (TIR) przetworzone na ortorektyfikowaną temperaturę powierzchni. Zawierają one również pasma pośrednie używane do obliczania produktów ST, a także pasma kontroli jakości.

Produkty Landsat 9 SR są tworzone za pomocą kodu odbicia powierzchni lądu (LaSRC). Wszystkie produkty ST z kolekcji 2 są tworzone za pomocą algorytmu jednokanałowego opracowanego wspólnie przez Rochester Institute of Technology (RIT) i Jet Propulsion Laboratory (JPL) należące do National Aeronautics and Space Administration (NASA).

Zebrane dane są dzielone na nakładające się „sceny” o wymiarach około 170 km x 183 km przy użyciu standardowej siatki odniesienia.

Niektóre zasoby mają tylko dane SR, w których przypadku pasma ST są obecne, ale puste. W przypadku zasobów z pasmem ST i SR parametr „PROCESSING_LEVEL” jest ustawiony na „L2SP”. W przypadku komponentów zawierających tylko pasma SR wartość parametru „PROCESSING_LEVEL” jest ustawiona na „L2SR”.

Dodatkowa dokumentacja i przykłady użycia

Pliki Landsat Collection 2 są dostępne publicznie w zasobniku Google Cloud Storage na zasadzie płatności przez zamawiającego. Pliki są indeksowane w regularnie aktualizowanej tabeli BigQuery, aby ułatwić analizę: earth-engine-public-data.geo_index.landsat_c2_index.

Uwagi dostawcy danych:

  • Aby można było skutecznie przetworzyć dane produktów na temperaturę powierzchni, muszą one zawierać zarówno dane optyczne, jak i termiczne, ponieważ do czasowego dostosowania produktu ASTER GED do docelowej sceny Landsat wymagany jest wskaźnik ASTER NDVI. Dlatego nie można przetworzyć danych uzyskanych w nocy, aby określić temperaturę powierzchni.

  • W przypadku danych o temperaturze powierzchni w odniesieniu do chmur i prawdopodobnie cieni chmur występuje znany błąd. Charakterystyka tych problemów została udokumentowana przez Cooka i współautorów (2014).

  • ASTER GED zawiera obszary, w których brakuje danych o średniej emisyjności wymaganych do prawidłowego generowania produktów ST. Jeśli brakuje informacji ASTER GED, w tych obszarach będą też brakować dane ST.

  • Zbiór danych ASTER GED jest tworzony na podstawie wszystkich pikseli scen ASTER, które zostały zarejestrowane w latach 2000–2008 przy bezchmurnym niebie. Ten zbiór danych ma globalny zasięg przestrzenny, ale w niektórych obszarach brakuje informacji o średniej emisyjności z powodu ciągłego zanieczyszczenia chmurami w pomiarach ASTER.

  • USGS dodatkowo sprawdza w ASTER GED wartości niefizyczne (emisyjność < 0,6), aby usunąć wszelkie niedoszacowania emisyjności spowodowane niewykrytymi chmurami. W przypadku każdego piksela bez danych wejściowych ASTER GED lub z niefizyczną wartością emisyjności wynikowe produkty Landsat ST mają brakujące piksele. Brakujące piksele Landsat ST będą spójne w czasie (od 1982 r. do dziś) ze względu na statyczny charakter danych klimatycznych ASTER GED. Więcej informacji znajdziesz w artykule landsat-collection-2-surface-temperature-data-gaps-due-missing.

Pasma

Rozmiar piksela
30 metrów

Pasma

Nazwa Jednostki Min. Maks. Skaluj Przesunięcie Rozmiar piksela Długość fali Opis
SR_B1 1 65455 2,75e-05 -0,2 metry 0,435–0,451 μm

Odbicie powierzchniowe pasma 1 (ultraniebieskie, aerozol przybrzeżny)

SR_B2 1 65455 2,75e-05 -0,2 metry 0,452–0,512 μm

Odbicie powierzchniowe pasma 2 (niebieskiego)

SR_B3 1 65455 2,75e-05 -0,2 metry 0,533–0,590 μm

Odbicie powierzchniowe w paśmie 3 (zielonym)

SR_B4 1 65455 2,75e-05 -0,2 metry 0,636–0,673 μm

Odbicie powierzchniowe pasma 4 (czerwonego)

SR_B5 1 65455 2,75e-05 -0,2 metry 0,851–0,879 μm

Odbicie powierzchniowe w paśmie 5 (bliska podczerwień)

SR_B6 1 65455 2,75e-05 -0,2 metry 1,566–1,651 μm

Pasmo 6 (krótkofalowe promieniowanie podczerwone 1) – współczynnik odbicia powierzchni

SR_B7 1 65455 2,75e-05 -0,2 metry 2,107–2,294 μm

Współczynnik odbicia powierzchni w paśmie 7 (krótkofalowe promieniowanie podczerwone 2)

SR_QA_AEROSOL metry Brak

Atrybuty aerozolu

ST_B10 K 0 65535 0,00341802 149 metry 10,60–11,19 μm

Temperatura powierzchni w paśmie 10. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_ATRAN 0 10000 0,0001 metry Brak

Transmisja atmosferyczna. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_CDIST km 0 24000 0,01 metry Brak

Odległość piksela od chmury. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_DRAD W/(m^2*sr*um)/ DN 0 28000 0,001 metry Brak

Downwelled Radiance. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_EMIS 0 10000 0,0001 metry Brak

Emisyjność pasma 10 oszacowana na podstawie danych ASTER GED. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_EMSD 0 10000 0,0001 metry Brak

Odchylenie standardowe emisyjności. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_QA K 0 32767 0,01 metry Brak

Niepewność pomiaru temperatury skóry. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_TRAD W/(m^2*sr*um)/ DN 0 22000 0,001 metry Brak

Pasmo termiczne przekonwertowane na radiację. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

ST_URAD W/(m^2*sr*um)/ DN 0 28000 0,001 metry Brak

Promieniowanie skierowane ku górze. Jeśli parametr „PROCESSING_LEVEL” ma wartość „L2SR”, to pasmo jest w pełni zamaskowane.

QA_PIXEL metry Brak

Atrybuty jakości pikseli wygenerowane na podstawie algorytmu CFMASK.

QA_RADSAT metry Brak

Kontrola jakości nasycenia radiometrycznego

Właściwości obrazu

Właściwości obrazu

Nazwa Typ Opis
ALGORITHM_SOURCE_SURFACE_REFLECTANCE CIĄG ZNAKÓW

Nazwa i wersja algorytmu odbicia powierzchniowego.

ALGORITHM_SOURCE_SURFACE_TEMPERATURE CIĄG ZNAKÓW

Nazwa i wersja algorytmu pomiaru temperatury powierzchni.

CLOUD_COVER LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Procent zachmurzenia (0–100), -1 = nie obliczono.

CLOUD_COVER_LAND LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Procent zachmurzenia nad lądem (0–100), -1 = nie obliczono.

COLLECTION_CATEGORY CIĄG ZNAKÓW

Kategoria kolekcji scen, „T1” lub „T2”.

DATA_SOURCE_AIR_TEMPERATURE CIĄG ZNAKÓW

Źródło danych o temperaturze powietrza.

DATA_SOURCE_ELEVATION CIĄG ZNAKÓW

Źródło danych o wysokości.

DATA_SOURCE_OZONE CIĄG ZNAKÓW

Źródło danych o ozonie.

DATA_SOURCE_PRESSURE CIĄG ZNAKÓW

Źródło danych o ciśnieniu.

DATA_SOURCE_REANALYSIS CIĄG ZNAKÓW

Źródło danych ponownej analizy.

DATA_SOURCE_TIRS_STRAY_LIGHT_CORRECTION CIĄG ZNAKÓW

Źródło danych do korekcji światła rozproszonego TIRS.

DATA_SOURCE_WATER_VAPOR CIĄG ZNAKÓW

Źródło danych o parze wodnej.

DATE_PRODUCT_GENERATED LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Sygnatura czasowa daty wygenerowania produktu.

EARTH_SUN_DISTANCE LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Odległość Ziemi od Słońca (AU).

EPHEMERIS_TYPE CIĄG ZNAKÓW

Identyfikator informujący użytkownika o typie efemerydy orbitalnej: „DEFINITIVE” lub „PREDICTIVE”. Jeśli to pole nie występuje, użytkownik powinien założyć, że jest to „PREDICTIVE”.

GEOMETRIC_RMSE_MODEL LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Łączna średnia kwadratowa błędów (RMSE) reszt geometrycznych (w metrach) w kierunku poprzecznym i podłużnym. Ten parametr występuje tylko wtedy, gdy parametr L1_PROCESSING_LEVEL ma wartość L1TP.

GEOMETRIC_RMSE_MODEL_X LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Średnia kwadratowa błędów (RMSE) reszt geometrycznych (w metrach) mierzonych na punktach kontrolnych używanych do korekty precyzji geometrycznej w kierunku poprzecznym. Ten parametr występuje tylko wtedy, gdy L1_PROCESSING_LEVEL ma wartość L1TP.

GEOMETRIC_RMSE_MODEL_Y LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Średnia kwadratowa błędów (RMSE) reszt geometrycznych (w metrach) zmierzonych na punktach kontrolnych użytych do korekty precyzji geometrycznej w kierunku wzdłuż ścieżki. Ten parametr występuje tylko wtedy, gdy L1_PROCESSING_LEVEL ma wartość L1TP.

GEOMETRIC_RMSE_VERIFY LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Błąd RMSE reszt geometrycznych (w metrach) zmierzony na produkcie skorygowanym pod względem terenu niezależnie za pomocą GLS2000. Ten parametr jest obecny tylko wtedy, gdy poziom przetwarzania L1_PROCESSING_LEVEL to L1TP.

GROUND_CONTROL_POINTS_MODEL LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Liczba punktów kontrolnych użytych w procesie korekty precyzji. Ten parametr jest obecny tylko wtedy, gdy L1_PROCESSING_LEVEL ma wartość L1TP.

GROUND_CONTROL_POINTS_VERIFY LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Liczba punktów kontrolnych na ziemi użytych do weryfikacji produktu skorygowanego pod względem terenu.

GROUND_CONTROL_POINTS_VERSION LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Wersja zbioru danych GCP używana w procesie korekty precyzji. Ten parametr jest obecny tylko wtedy, gdy L1_PROCESSING_LEVEL to L1TP.

IMAGE_QUALITY_OLI PRZ

Jakość obrazu pasm OLI. 1 – najgorsza, 9 – najlepsza, 0 – jakość nie została obliczona. W przypadku Landsata 8 ten parametr jest obniżany w przypadku scen zebranych przy użyciu dolnych 12 bitów z czujnika OLI (TRUNCATION_OLI = "LOWER").

IMAGE_QUALITY_TIRS PRZ

Jakość obrazu pasm TIRS. 1 – najgorsza, 9 – najlepsza, 0 – jakość nie została obliczona. Jest ona również korygowana w dół w przypadku scen przetwarzanych z wartością „SWITCHED” w polu TIRS_SSM_POSITION_STATUS.

L1_DATE_PRODUCT_GENERATED CIĄG ZNAKÓW

Data wygenerowania produktu L1.

L1_LANDSAT_PRODUCT_ID CIĄG ZNAKÓW

Identyfikator produktu Landsat dla odpowiedniego produktu poziomu 1.

L1_PROCESSING_LEVEL CIĄG ZNAKÓW

Poziom przetwarzania odpowiedniego produktu L1.

L1_PROCESSING_SOFTWARE_VERSION CIĄG ZNAKÓW

Wersja oprogramowania do przetwarzania odpowiedniego produktu L1.

LANDSAT_PRODUCT_ID CIĄG ZNAKÓW

Identyfikator produktu Landsat

LANDSAT_SCENE_ID CIĄG ZNAKÓW

Krótki identyfikator sceny Landsat

PROCESSING_LEVEL CIĄG ZNAKÓW

„L2SP”, gdy występują pasma SR i LST, lub „L2SR”, gdy występują tylko pasma SR.

PROCESSING_SOFTWARE_VERSION CIĄG ZNAKÓW

Wersja oprogramowania do przetwarzania, w którym utworzono produkt.

ROLL_ANGLE LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Wartość kąta przechyłu statku kosmicznego w centrum sceny. Wartość przechyłu jest podawana w odniesieniu do układu Yaw Steering Frame (YSF), którego oś X jest wyrównana z wektorem chwilowej prędkości ścieżki naziemnej. Obrót wokół tej osi X odbywa się zgodnie z regułą prawej dłoni: dodatni obrót powoduje, że przyrządy są skierowane na lewo od ścieżki naziemnej, a ujemny obrót powoduje, że przyrządy są skierowane na prawo.

SCENE_CENTER_TIME CIĄG ZNAKÓW

Czas obserwacji w formacie ISO 8601.

SENSOR_ID CIĄG ZNAKÓW

Nazwa czujnika.

SPACECRAFT_ID CIĄG ZNAKÓW

Nazwa statku kosmicznego.

SUN_AZIMUTH LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Kąt azymutu słońca w stopniach dla lokalizacji środka obrazu w momencie jego rejestracji. Wartość dodatnia oznacza kąty na wschód lub zgodnie z ruchem wskazówek zegara od północy. Wartość ujemna oznacza kąty na zachód lub w kierunku przeciwnym do ruchu wskazówek zegara od północy.

SUN_ELEVATION LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Kąt elewacji słońca w stopniach dla lokalizacji środka obrazu w momencie jego rejestracji. Wartość dodatnia oznacza scenę dzienną. Wartość ujemna oznacza scenę nocną. Uwaga: do obliczenia współczynnika odbicia potrzebny jest kąt zenitu słonecznego, który wynosi 90 stopni minus kąt elewacji słońca.

TARGET_WRS_PATH PRZ

Najbliższa ścieżka WRS-2 do środka sceny obrazu Line-ofSight (LOS).

TARGET_WRS_ROW PRZ

Najbliższy wiersz WRS-2 do środka sceny LOS obrazu. Wiersze 880–889 są zarezerwowane dla bieguna północnego, a 990–999 dla bieguna południowego, gdzie WRS-2 nie jest zdefiniowany.

TEMPERATURE_MAXIMUM_BAND_ST_B10 LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Maksymalna osiągalna wartość temperatury dla pasma 10.

TEMPERATURE_MINIMUM_BAND_ST_B10 LICZBA ZMIENNOPRZECINKOWA O PODWÓJNEJ PRECYZJI

Minimalna osiągalna wartość temperatury dla pasma 10.

TIRS_SSM_MODEL CIĄG ZNAKÓW

Wskazuje, w jaki sposób określono pozycję lustra skanującego (SSM) w scenie TIRS z satelity Landsat 8.

  • „FINAL” oznacza ostateczne szacunkowe wartości enkodera wygenerowane po zdarzeniu przełączenia.

  • „ACTUAL” oznacza rzeczywiste wartości enkodera.

TIRS_SSM_POSITION_STATUS CIĄG ZNAKÓW
  • „NOMINAL” oznacza, że SSM działał prawidłowo w przypadku tej sceny.

  • „SWITCHED” oznacza, że SSM zmienił tryb działania w scenie i może mieć problemy z jakością obrazu TIRS, co ma bezpośredni wpływ na wartość IMAGE_QUALITY_TIRS.

  • „ESTIMATED” oznacza, że pozycja SSM została oszacowana, co może być mniej dokładne niż stan „NOMINAL”.

TRUNCATION_OLI CIĄG ZNAKÓW
  • „LOWER” oznacza, że użyto 12 bitów o mniejszej wadze.

  • „UPPER” oznacza, że użyto 12 bitów wyższego rzędu.

Normalny tryb obcinania to „UPPER”.

WRS_PATH PRZ

Numer ścieżki WRS sceny.

WRS_ROW PRZ

Numer wiersza WRS sceny.

Warunki korzystania z usługi

Warunki korzystania z usługi

Zbiory danych Landsat są tworzone przez instytucje federalne, dlatego należą do domeny publicznej i można ich używać, przenosić i reprodukować bez ograniczeń wynikających z praw autorskich.

Należy podać informację o źródle danych USGS, dodając wiersz z cytatem, np. jak w przykładzie poniżej.

(Nazwa produktu, obrazu, zdjęcia lub zbioru danych) dzięki uprzejmości U.S. Geological Survey

Przykład: zdjęcie Landsat-7 udostępnione przez U.S. Geological Survey

Więcej informacji o prawidłowym cytowaniu i uznawaniu produktów USGS znajdziesz w wytycznych dotyczących systemu tożsamości wizualnej USGS.

Odkrywanie za pomocą Earth Engine

Edytor kodu (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)');

Konfiguracja Pythona

Informacje o interfejsie Python API i używaniu geemap do interaktywnego programowania znajdziesz na stronie środowiska 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
Otwórz w edytorze kodu