USGS Landsat 9 Level 2, Collection 2, Tier 1

LANDSAT/LC09/C02/T1_L2
זמינות קבוצת הנתונים
2021-10-31T00:00:00Z–2026-05-31T23:02:42.951000Z
מפיק מערך הנתונים
קטע קוד ל-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, וגם רצועות של בקרת איכות.

מוצרי Landsat 9 SR נוצרים באמצעות קוד ההחזרה של פני הקרקע (LaSRC). כל מוצרי ST של Collection 2 נוצרים באמצעות אלגוריתם חד-ערוצי שנוצר בשיתוף פעולה בין Rochester Institute of Technology ‏ (RIT) לבין National Aeronautics and Space Administration ‏ (NASA) Jet Propulsion Laboratory ‏ (JPL).

רצועות של נתונים שנאספו נארזות ב'סצנות' חופפות ששטח כל אחת מהן הוא בערך 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 μm

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

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

טמפרטורת פני השטח של פס 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)

המרחק בין כדור הארץ לשמש (יחידות אסטרונומיות).
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)

מספר נקודות הבקרה הקרקעיות שנעשה בהן שימוש בתהליך תיקון הדיוק. הפרמטר הזה מופיע רק אם 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 למוצר התואם ברמה 1.
L1_PROCESSING_LEVEL מחרוזת

רמת העיבוד של המוצר התואם ברמה 1.
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.

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

  • ‫LOWER מציין שנעשה שימוש ב-12 הסיביות התחתונות.

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

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

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

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

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

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

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

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

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

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

פרטים נוספים על ציטוט נכון של מוצרים של USGS ועל מתן קרדיט ל-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 מפורט מידע על 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
פתיחה ב-Code Editor
קבוצת הנתונים הזו מכילה נתונים על השתקפות פני השטח שעברו תיקון אטמוספרי ועל טמפרטורת פני הקרקע, שנגזרים מהנתונים שנוצרו על ידי חיישני OLI/TIRS של Landsat 9. התמונות האלה מכילות 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-05-31T23:02:42.951000Z
-90 -180 90 180
Google Earth Engine