USGS Landsat 7 Level 2, Collection 2, Tier 1

LANDSAT/LE07/C02/T1_L2
זמינות קבוצת הנתונים
1999-05-28T01:02:17Z–2024-01-19T01:39:58.042000Z
יוצר מערך הנתונים
קטע קוד ל-Earth Engine
ee.ImageCollection("LANDSAT/LE07/C02/T1_L2")
מרווח הזמן בין כל אחזור של נתונים עדכניים
16 ימים
תגים
cfmask cloud etm fmask global landsat lasrc le07 lst reflectance satellite-imagery sr usgs

תיאור

קבוצת הנתונים הזו מכילה נתוני השתקפות מפני השטח וטמפרטורת פני השטח, שעברו תיקון אטמוספרי והופקו מהנתונים שנוצרו על ידי חיישן Landsat 7 ETM+. התמונות האלה מכילות 4 פסי תדרים של אור נראה ואינפרה-אדום קרוב (VNIR) ו-2 פסי תדרים של אינפרה-אדום בגלים קצרים (SWIR) שעברו עיבוד להאחדת קנה מידה של השתקפות פני השטח, ופס תדרים אחד של אינפרה-אדום תרמי (TIR) שעבר עיבוד להאחדת קנה מידה של טמפרטורת פני השטח. הם מכילים גם רצועות ביניים שמשמשות לחישוב מוצרי ה-ST, וגם רצועות QA.

מוצרי Landsat 7 SR נוצרים באמצעות אלגוריתם Landsat Ecosystem Disturbance Adaptive Processing System ‏ (LEDAPS) (גרסה 3.4.0). כל מוצרי Collection 2 ST נוצרים באמצעות אלגוריתם חד-ערוצי שנוצר בשיתוף פעולה בין 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.

  • בנוסף, ב-ASTER GED,‏ USGS מסנן ערכים לא פיזיים (מקדם פליטה < 0.6) כדי להסיר כל הערכת חסר של מקדם הפליטה בגלל עננים שלא זוהו. בכל פיקסל שאין לו קלט ASTER GED או ערך פליטה לא פיזיקלי, מוצרי Landsat ST שמתקבלים כוללים פיקסלים חסרים. הפיקסלים החסרים של Landsat ST יהיו עקביים לאורך זמן (1982 עד היום) בהתחשב באופי הסטטי של נתוני האקלים הממוצעים של ASTER GED. מידע נוסף זמין במאמר בנושא landsat-collection-2-surface-temperature-data-gaps-due-missing

שימו לב שהמסלול של Landsat 7 סטה לזמן רכישה מוקדם יותר מאז 2017.

תחום תדרים

גודל הפיקסל
‫30 מטרים

תחום תדרים

שם יחידות מינימום מקסימום קנה מידה היסט גודל הפיקסל אורך גל תיאור
SR_B1 1 65455 2.75e-05 ‫0.2- מטרים ‫0.45-0.52 מיקרומטר

מידת ההחזרה של פני השטח ברצועה 1 (כחול)

SR_B2 1 65455 2.75e-05 ‫0.2- מטרים ‫0.52-0.60 μm

מידת ההחזרה של פני השטח ברצועה 2 (ירוק)

SR_B3 1 65455 2.75e-05 ‫0.2- מטרים ‫0.63-0.69 μm

מידת ההחזרה של פני השטח ברצועה 3 (אדום)

SR_B4 1 65455 2.75e-05 ‫0.2- מטרים ‫0.77-0.90 μm

מידת ההחזרה של פני השטח ברצועה 4 (אינפרא-אדום קרוב)

SR_B5 1 65455 2.75e-05 ‫0.2- מטרים ‫1.55-1.75 מיקרומטר

מידת ההחזרה של פני השטח ברצועה 5 (תת-אדום בגלים קצרים 1)

SR_B7 1 65455 2.75e-05 ‫0.2- מטרים ‫2.08-2.35 מיקרומטר

מידת ההחזרה של פני השטח ברצועה 7 (תת-אדום בגלים קצרים 2)

SR_ATMOS_OPACITY 0 10000 0.001 מטרים ללא

פרשנות כללית של אטימות האטמוספירה שנוצרה על ידי LEDAPS ומבוססת על הקרינה שנצפתה מעל צמחייה כהה וצפופה (DDV) בסצנה. הפרשנות הכללית של אטימות האטמוספירה היא שערכים (אחרי החלת קנה מידה של 0.001) שקטנים מ-0.1 הם ברורים, ערכים בין 0.1 ל-0.3 הם ממוצעים, וערכים שגדולים מ-0.3 מצביעים על אובך או על מצבי עננות אחרים. ערכי SR מפיקסלים עם אטימות אטמוספרית גבוהה יהיו פחות מהימנים, במיוחד בתנאים של זווית זניתית סולארית גבוהה. הפס SR_ATMOS_OPACITY מיועד למשתמשים מתקדמים ולבדיקת איכות המוצר, והוא לא אומת. במקום זאת, מומלץ לרוב המשתמשים להשתמש במידע על פס QA_PIXEL כדי להבחין בין ענן לבין מכשיר.

SR_CLOUD_QA מטרים ללא

Cloud Quality Assessment

ST_B6 K 0 65535 0.00341802 149 מטרים ‫‎10.40-12.50 μm

טמפרטורת פני השטח ברצועה 6. אם הערך של 'PROCESSING_LEVEL' מוגדר כ-'L2SR', הרצועה הזו מוסתרת באופן מלא.

ST_ATRAN 0 10000 0.0001 מטרים ללא

העברה אטמוספרית. אם הערך של 'PROCESSING_LEVEL' מוגדר כ-'L2SR', הרצועה הזו מוסתרת באופן מלא.

ST_CDIST ק"מ 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. אם הערך של '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 מטרים ללא

Upwelled Radiance. אם הערך של 'PROCESSING_LEVEL' מוגדר כ-'L2SR', הרצועה הזו מוסתרת באופן מלא.

QA_PIXEL מטרים ללא

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

QA_RADSAT מטרים ללא

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

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

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

שם סוג תיאור
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_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.

GROUND_CONTROL_POINTS_MODEL נקודה צפה בדיוק כפול (Double Precision Floating Point)

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

GROUND_CONTROL_POINTS_VERSION נקודה צפה בדיוק כפול (Double Precision Floating Point)

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

IMAGE_QUALITY INT

איכות התמונה המורכבת של הפסים. ‫0 = הכי גרוע, 9 = הכי טוב, ‎-1 = האיכות לא חושבה או הוערכה.

L1_DATE_PRODUCT_GENERATED מחרוזת

תאריך יצירת המוצר עבור מוצר L1 התואם.

L1_LANDSAT_PRODUCT_ID מחרוזת

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

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 מחרוזת

גרסת התוכנה לעיבוד שיצרה את המוצר.

SCENE_CENTER_TIME מחרוזת

השעה של התצפיות בפורמט מחרוזת ISO 8601.

SENSOR_ID מחרוזת

שם החיישן.

SPACECRAFT_ID מחרוזת

שם החללית.

SUN_AZIMUTH נקודה צפה בדיוק כפול (Double Precision Floating Point)

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

SUN_ELEVATION נקודה צפה בדיוק כפול (Double Precision Floating Point)

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

TEMPERATURE_MAXIMUM_BAND_ST_B6 נקודה צפה בדיוק כפול (Double Precision Floating Point)

ערך הטמפרטורה המקסימלי שאפשר להשיג עבור פס 6.

TEMPERATURE_MINIMUM_BAND_ST_B6 נקודה צפה בדיוק כפול (Double Precision Floating Point)

ערך הטמפרטורה המינימלי שאפשר להשיג בתדר 6.

WRS_PATH INT

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

WRS_ROW INT

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

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

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

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

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

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

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

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

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

Code Editor (JavaScript)

var dataset = ee.ImageCollection('LANDSAT/LE07/C02/T1_L2')
    .filterDate('2017-06-01', '2017-07-01');

// Applies scaling factors.
function applyScaleFactors(image) {
  var opticalBands = image.select('SR_B.').multiply(0.0000275).add(-0.2);
  var thermalBand = image.select('ST_B6').multiply(0.00341802).add(149.0);
  return image.addBands(opticalBands, null, true)
              .addBands(thermalBand, null, true);
}

dataset = dataset.map(applyScaleFactors);

var visualization = {
  bands: ['SR_B3', 'SR_B2', 'SR_B1'],
  min: 0.0,
  max: 0.3,
};

Map.setCenter(-114.2579, 38.9275, 8);

Map.addLayer(dataset, visualization, 'True Color (321)');

הגדרה של Python

בדף סביבת Python מפורט מידע על Python API ועל השימוש ב-geemap לפיתוח אינטראקטיבי.

import ee
import geemap.core as geemap

Colab (Python)

dataset = ee.ImageCollection('LANDSAT/LE07/C02/T1_L2').filterDate(
    '2017-06-01', '2017-07-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_B6').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_B3', 'SR_B2', 'SR_B1'],
    'min': 0.0,
    'max': 0.3,
}

m = geemap.Map()
m.set_center(-114.2579, 38.9275, 8)
m.add_layer(dataset, visualization, 'True Color (321)')
m
פתיחה ב-Code Editor