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. התמונות האלה מכילות 5 פסי אור נראה ואינפרה-אדום קרוב (VNIR) ו-2 פסי אינפרה-אדום בגלים קצרים (SWIR) שעברו עיבוד להאחדת קנה מידה של השתקפות פני השטח, ופס אחד של אינפרה-אדום תרמי (TIR) שעבר עיבוד להאחדת קנה מידה של טמפרטורת פני השטח. הם מכילים גם רצועות ביניים שמשמשות לחישוב מוצרי ה-ST, וגם רצועות QA.

מוצרי Landsat 9 SR נוצרים באמצעות קוד ההחזרה של פני הקרקע (LaSRC). כל מוצרי Collection 2 ST נוצרים באמצעות אלגוריתם של ערוץ יחיד שנוצר במשותף על ידי המכון הטכנולוגי של רוצ'סטר (RIT) והמעבדה להנעה סילונית (JPL) של מנהל האווירונאוטיקה והחלל הלאומי (NASA).

רצועות של נתונים שנאספו נארזות ב"סצנות" חופפות שמכסות כ-170 ק"מ על 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 μm

החזרת אור מפני השטח ברצועה 2 (כחול)
SR_B3 1 65455 2.75e-05 -0.2 ‫30 מטרים ‎0.533-0.590 μm

מידת ההחזרה של פני השטח ברצועה 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 μm

החזרת אור מפני השטח ברצועה 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 ק"מ 0 24000 0.01 ‫30 מטרים ללא

המרחק של הפיקסל מהענן. אם הערך של 'PROCESSING_LEVEL' הוא 'L2SR', הרצועה הזו מוסתרת לחלוטין.
ST_DRAD W/(m^2*sr*um)/ 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 W/(m^2*sr*um)/ DN 0 22000 0.001 ‫30 מטרים ללא

פס תרמי שהומר לבהירות קרינה. אם הערך של PROCESSING_LEVEL מוגדר כ-L2SR, הרצועה הזו מוסתרת לגמרי.
ST_URAD W/(m^2*sr*um)/ DN 0 28000 0.001 ‫30 מטרים ללא

קרינה עולה. אם הערך של PROCESSING_LEVEL מוגדר כ-L2SR, הרצועה הזו מוסתרת לגמרי.
QA_PIXEL ‫30 מטרים ללא

מאפייני איכות הפיקסלים שנוצרו מאלגוריתם CFMASK.
QA_RADSAT ‫30 מטרים ללא

בקרת איכות של רוויה רדיומטרית

מאפייני תמונה

מאפייני תמונה

שם סוג תיאור
ALGORITHM_SOURCE_SURFACE_REFLECTANCE מחרוזת

השם והגרסה של האלגוריתם של החזרת האור מהמשטח.
ALGORITHM_SOURCE_SURFACE_TEMPERATURE מחרוזת

השם והגרסה של האלגוריתם לחישוב טמפרטורת פני השטח.
CLOUD_COVER נקודה צפה בדיוק כפול (Double Precision Floating Point)

אחוז כיסוי העננים (0-100), ‎-1 = לא מחושב.
CLOUD_COVER_LAND נקודה צפה בדיוק כפול (Double Precision Floating Point)

אחוז כיסוי העננים מעל היבשה (0-100), ‎-1 = לא מחושב.
COLLECTION_CATEGORY מחרוזת

קטגוריית אוסף הסצנות, "T1" או "T2".
DATA_SOURCE_AIR_TEMPERATURE מחרוזת

מקור הנתונים של טמפרטורת האוויר.
DATA_SOURCE_ELEVATION מחרוזת

מקור נתוני הגובה.
DATA_SOURCE_OZONE מחרוזת

מקור נתוני האוזון.
DATA_SOURCE_PRESSURE מחרוזת

מקור נתוני הלחץ.
DATA_SOURCE_REANALYSIS מחרוזת

מקור נתונים לניתוח מחדש.
DATA_SOURCE_TIRS_STRAY_LIGHT_CORRECTION מחרוזת

מקור נתונים לתיקון אור תועה של TIRS.
DATA_SOURCE_WATER_VAPOR מחרוזת

מקור הנתונים של אדי המים.
DATE_PRODUCT_GENERATED נקודה צפה בדיוק כפול (Double Precision Floating Point)

חותמת הזמן של התאריך שבו המוצר נוצר.
EARTH_SUN_DISTANCE נקודה צפה בדיוק כפול (Double Precision Floating Point)

המרחק בין כדור הארץ לשמש (AU).
EPHEMERIS_TYPE מחרוזת

מזהה שמודיע למשתמש על סוג האפמריס של המסלול שבו נעשה שימוש: "DEFINITIVE" או "PREDICTIVE". אם השדה לא קיים, המשתמש צריך להניח שהערך הוא PREDICTIVE.
GEOMETRIC_RMSE_MODEL נקודה צפה בדיוק כפול (Double Precision Floating Point)

ה-RMSE (שגיאת השורש הממוצע הריבועי) המשולב של השאריות הגיאומטריות (במטרים) בכיוונים של תנועה לרוחב המסלול ותנועה לאורך המסלול. הפרמטר הזה מופיע רק אם L1_PROCESSING_LEVEL הוא L1TP.
GEOMETRIC_RMSE_MODEL_X נקודה צפה בדיוק כפול (Double Precision Floating Point)

‫RMSE (שורש ממוצע הריבועים) של השאריות הגיאומטריות (במטרים) שנמדדו ב-GCP (נקודות בקרה קרקעיות) ששימשו לתיקון הדיוק הגיאומטרי בכיוון הרוחבי. הפרמטר הזה קיים רק אם הערך של L1_PROCESSING_LEVEL הוא L1TP.
GEOMETRIC_RMSE_MODEL_Y נקודה צפה בדיוק כפול (Double Precision Floating Point)

‫RMSE (שורש ממוצע הריבועים) של השאריות הגיאומטריות (במטרים) שנמדדו ב-GCP (נקודות בקרה קרקעיות) ששימשו לתיקון הדיוק הגיאומטרי בכיוון לאורך המסלול. הפרמטר הזה קיים רק אם הערך של L1_PROCESSING_LEVEL הוא L1TP.
GEOMETRIC_RMSE_VERIFY נקודה צפה בדיוק כפול (Double Precision Floating Point)

שגיאת השורש הממוצע הריבועי (RMSE) של השאריות הגיאומטריות (במטרים) שנמדדו במוצר עם תיקון פני השטח באופן עצמאי באמצעות GLS2000. הפרמטר הזה קיים רק אם הערך של L1_PROCESSING_LEVEL הוא L1TP.
GROUND_CONTROL_POINTS_MODEL נקודה צפה בדיוק כפול (Double Precision Floating Point)

מספר נקודות ה-GCP שנעשה בהן שימוש בתהליך תיקון הדיוק. הפרמטר הזה קיים רק אם הערך של L1_PROCESSING_LEVEL הוא L1TP.
GROUND_CONTROL_POINTS_VERIFY נקודה צפה בדיוק כפול (Double Precision Floating Point)

מספר נקודות הבקרה הקרקעיות (GCP) שנעשה בהן שימוש באימות של המוצר עם תיקון פני השטח.
GROUND_CONTROL_POINTS_VERSION נקודה צפה בדיוק כפול (Double Precision Floating Point)

גרסת מערך הנתונים של GCP שמשמשת בתהליך תיקון הדיוק. הפרמטר הזה קיים רק אם הערך של L1_PROCESSING_LEVEL הוא L1TP.
IMAGE_QUALITY_OLI INT

איכות התמונה של רצועות ה-OLI. ‫1 = הכי גרוע, 9 = הכי טוב, 0 = האיכות לא מחושבת. ב-Landsat 8, הפרמטר הזה מותאם כלפי מטה לסצנות שנאספו באמצעות 12 הביטים התחתונים מחיישן OLI (TRUNCATION_OLI = "LOWER").
IMAGE_QUALITY_TIRS INT

איכות התמונה של פסי ה-TIRS. ‫1 = הכי גרוע, 9 = הכי טוב, 0 = האיכות לא מחושבת. הוא גם מותאם כלפי מטה לסצנות שעובדו עם "SWITCHED" לערך TIRS_SSM_POSITION_STATUS.
L1_DATE_PRODUCT_GENERATED מחרוזת

תאריך יצירת המוצר של מוצר L1 התואם.
L1_LANDSAT_PRODUCT_ID מחרוזת

מזהה מוצר Landsat למוצר L1 המתאים.
L1_PROCESSING_LEVEL מחרוזת

רמת העיבוד של מוצר L1 התואם.
L1_PROCESSING_SOFTWARE_VERSION מחרוזת

גרסת תוכנת העיבוד של מוצר L1 התואם.
LANDSAT_PRODUCT_ID מחרוזת

מזהה המוצר של Landsat
LANDSAT_SCENE_ID מחרוזת

מזהה קצר של סצנת Landsat
PROCESSING_LEVEL מחרוזת

‫L2SP אם יש גם פס SR וגם פס LST, או L2SR אם יש רק פס SR.
PROCESSING_SOFTWARE_VERSION מחרוזת

גרסת תוכנת העיבוד שבאמצעותה נוצר המוצר.
ROLL_ANGLE נקודה צפה בדיוק כפול (Double Precision Floating Point)

זווית הגלגול של החללית במרכז הסצנה. ערך הגלגול ניתן בהתייחסות למסגרת ההיגוי של הסבסוב (YSF), שציר ה-X שלה מיושר עם וקטור המהירות של מסלול הקרקע המיידי. הטיה סביב ציר ה-X מתבצעת לפי כלל יד ימין: הטיה חיובית גורמת למכשירים להצביע שמאלה ממסלול הקרקע, והטיה שלילית גורמת למכשיר להצביע ימינה.
SCENE_CENTER_TIME מחרוזת

זמן התצפיות, כמו במחרוזת ISO 8601.
SENSOR_ID מחרוזת

שם החיישן.
SPACECRAFT_ID מחרוזת

שם החללית.
SUN_AZIMUTH נקודה צפה בדיוק כפול (Double Precision Floating Point)

זווית האזימוט של השמש במעלות במיקום מרכז התמונה בזמן הצילום של מרכז התמונה. ערך חיובי מציין זוויות מזרחה או בכיוון השעון מצפון. ערך שלילי מציין זוויות לכיוון מערב או נגד כיוון השעון מהצפון.
SUN_ELEVATION נקודה צפה בדיוק כפול (Double Precision Floating Point)

זווית הגובה של השמש במעלות במיקום של מרכז התמונה בזמן הצילום של מרכז התמונה. ערך חיובי מציין סצנה שמתרחשת ביום. ערך שלילי מציין סצנה שמתרחשת בלילה. הערה: כדי לחשב את ההחזרה, צריך את זווית השמש בנקודת הזניט, שהיא 90 פחות זווית הגובה של השמש.
TARGET_WRS_PATH INT

נתיב ה-WRS-2 הקרוב ביותר למרכז הסצנה של קו הראייה (LOS) בתמונה.
TARGET_WRS_ROW INT

השורה הקרובה ביותר של WRS-2 למרכז הסצנה של LOS בתמונה. השורות 880 עד 889 שמורות לקוטב הצפוני, והשורות 990 עד 999 שמורות לקוטב הדרומי, שבהם לא מוגדר WRS-2.
TEMPERATURE_MAXIMUM_BAND_ST_B10 נקודה צפה בדיוק כפול (Double Precision Floating Point)

ערך הטמפרטורה המקסימלי שאפשר להשיג עבור רצועה 10.
TEMPERATURE_MINIMUM_BAND_ST_B10 נקודה צפה בדיוק כפול (Double Precision Floating Point)

ערך הטמפרטורה המינימלי שאפשר להשיג עבור פס 10.
TIRS_SSM_MODEL מחרוזת

מציין איך נקבע המיקום של Landsat 8 TIRS Scene Select Mirror‏ (SSM).

  • הערך FINAL מציין את הערכים המשוערים הסופיים של המקודד שנוצרו אחרי אירוע ההחלפה.

  • ‫ACTUAL: מציין את הערכים בפועל של המקודד.
TIRS_SSM_POSITION_STATUS מחרוזת

  • ‫NOMINAL מציין שה-SSM פעל בצורה תקינה בסצנה הזו.

  • ‫SWITCHED מציין שמצב הפעולה של ה-SSM השתנה בסצנה, ויכול להיות שיש בעיות באיכות התמונה של ה-TIRS, שמשפיעות ישירות על הערך IMAGE_QUALITY_TIRS.

  • הסטטוס ESTIMATED מציין שהמיקום של ה-SSM משוער, ויכול להיות שהוא לא מדויק כמו הסטטוס NOMINAL.
TRUNCATION_OLI מחרוזת

  • הערך "LOWER" מציין שנעשה שימוש ב-12 הביטים הנמוכים.

  • ‫UPPER מציין שנעשה שימוש ב-12 הביטים העליונים.

מצב החיתוך הרגיל הוא 'UPPER'.
WRS_PATH INT

מספר הנתיב של הסצנה ב-WRS.
WRS_ROW INT

מספר השורה של הסצנה ב-WRS.

תנאים והגבלות

תנאים והגבלות

מערכי הנתונים של Landsat הם נתונים שנוצרו על ידי הממשל הפדרלי, ולכן הם נחלת הכלל ואפשר להשתמש בהם, להעביר או לשכפל אותם ללא הגבלת זכויות יוצרים.

כדי לתת קרדיט ל-USGS כמקור הנתונים, צריך להוסיף שורת טקסט עם ציטוט, כמו בדוגמה הבאה.

(שם המוצר, התמונה, הצילום או מערך הנתונים) באדיבות הסקר הגאולוגי של ארצות הברית

דוגמה: תמונה של Landsat-7 באדיבות הסקר הגאולוגי של ארה"ב

בהנחיות של USGS בנושא מערכת הזהות החזותית מפורט מידע נוסף על ציטוט נכון ומתן קרדיט למוצרים של USGS.

סיור עם פלטפורמת Earth Engine

Code Editor (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 בשפת Python ועל השימוש ב-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
פתיחה ב-Code Editor
קבוצת הנתונים הזו מכילה נתונים על החזרת אור מפני השטח וטמפרטורת פני השטח, שתוקנו מבחינה אטמוספרית ונגזרו מהנתונים שהופקו על ידי חיישני Landsat 9 OLI/TIRS. התמונות האלה מכילות 5 פסים של אור נראה ואינפרה-אדום קרוב (VNIR) ו-2 פסים של אינפרה-אדום בגלים קצרים (SWIR) שעברו עיבוד להחזרת פני שטח שעברה האחדת קנה מידה, ופס אחד של אינפרה-אדום תרמי (TIR) שעבר עיבוד ל…
LANDSAT/LC09/C02/T1_L2, cfmask,cloud,fmask,global,l9sr,landsat,lasrc,lc09,lst,reflectance,satellite-imagery,sr,usgs
2021-10-31T00:00:00Z/2026-10-11T23:19:19.950000Z
-90 -180 90 180
Google Earth Engine