i.atcorrgrass - عبر الإنترنت في السحابة

برنامج:

اسم

أنا أتكور - يقوم بتصحيح الغلاف الجوي باستخدام خوارزمية 6S.
6S - المحاكاة الثانية لإشارة القمر الصناعي في الطيف الشمسي.

الكلمات الرئيسية

الصور وتصحيح الغلاف الجوي

موجز

أنا أتكور
أنا أتكور --مساعدة
أنا أتكور [-إراب] إدخال=الاسم [نطاق=الحد الأدنى - الحد الأقصى] [ارتفاع=الاسم] [رؤية=الاسم]
المعلمات=الاسم الناتج=الاسم [تجديد=الحد الأدنى - الحد الأقصى] [-اعادة الكتابة] [-مساعدة] [-مطنب]
[-هدوء] [-ui]

الأعلام:
-i
إخراج الخريطة النقطية كعدد صحيح

-r
تم تحويل الخريطة النقطية المدخلة إلى انعكاس (الافتراضي هو التألق)

-a
المدخلات من ETM + صورة تم التقاطها بعد 1 يوليو 2000

-b
المدخلات من ETM + صورة تم التقاطها قبل 1 يوليو 2000

--الكتابة فوق
السماح لملفات الإخراج بالكتابة فوق الملفات الموجودة

--مساعدة
طباعة ملخص الاستخدام

- الإسراف
إخراج وحدة مطول

--هادئ
إخراج وحدة هادئة

--ui
فرض إطلاق مربع حوار واجهة المستخدم الرسومية

المعلمات:
إدخال=الاسم [مطلوب]
اسم الخريطة النقطية المدخلة

نطاق=الحد الأدنى - الحد الأقصى
نطاق الإدخال
الافتراضي: 0,255

ارتفاع=الاسم
اسم الخريطة النقطية لإدخال الارتفاع (بالمتر)

رؤية=الاسم
اسم الخريطة النقطية لرؤية المدخلات (بالكيلومترات)

المعلمات=الاسم [مطلوب]
اسم ملف نصي الإدخال مع معلمات 6S

الناتج=الاسم [مطلوب]
اسم لخريطة نقطية الإخراج

تجديد=الحد الأدنى - الحد الأقصى
إعادة مقياس خريطة الناتج النقطية
الافتراضي: 0,255

الوصف

أنا أتكور يقوم بتصحيح الغلاف الجوي على الخريطة النقطية للإدخال باستخدام خوارزمية 6S
(الثاني محاكاة of الأقمار الصناعية سيجنل in القادم الطاقة الشمسية طيف). خوارزمية مفصلة
الوصف متاح في مرفق حوسبة علوم انعكاس سطح الأرض
موقع الكتروني.

مهم ملاحظة: حالياًّ منطقة إعدادات . تم تجاهله! تم تعديل المنطقة لتغطية
إدخال الخريطة النقطية قبل إجراء التصحيح الجوي. الإعدادات السابقة هي
استعادته بعد ذلك. يخبرنا هذا العلم أنا أتكور لمحاولة تسريع العمليات الحسابية. لكن،
سيؤدي هذا الخيار إلى زيادة متطلبات الذاكرة.

إذا علم -r يتم استخدام البيانات النقطية المدخلة على أنها الانعكاس. خلاف ذلك ، المدخلات
يتم التعامل مع البيانات النقطية على أنها إشعاع القيم ويتم تحويلها إلى انعكاس في
أنا أتكور مدة العرض. البيانات الناتجة هي دائما انعكاس.

لاحظ أنه يجب تحديد وقت تجاوز القمر الصناعي بتوقيت غرينتش (GMT).

مثال على معلمات 6S:
8 - الشروط الهندسية = Landsat ETM +
2 19 13.00 -47.410 -20.234 - شهر يوم hh.ddd خط عرض ("hh.ddd" بالساعات العشرية بتوقيت جرينتش)
1 - وضع الغلاف الجوي = استوائي
1- نموذج الهباء الجوي = قاري
15 - الرؤية [كم] (تركيز نموذج الهباء الجوي)
-0.600 - متوسط ​​الارتفاع المستهدف فوق مستوى سطح البحر [كم] (هنا 600 متر فوق سطح البحر)
-1000 - ارتفاع المستشعر (هنا ، مستشعر على متن قمر صناعي)
64 - النطاق الرابع من ETM + Landsat 4
إذا لم يكن الموضع متاحًا في خطوط الطول والعرض (WGS84) ، فإن ملف m.proj تحويل
يمكن استخدام الوحدة النمطية لإعادة الإسقاط من إسقاط مختلف.

6S CODE معامل خيارات

A. هندسي الشروط
رمز الوصف أيقونة

1 ميتيوسات تدخل الملاحظة الشهر واليوم والساعة العشرية (التوقيت العالمي hh.ddd) n. ل
العمود ، ن. على الخط. (المقياس الكامل 5000 * 2500)

2 يذهب الشرق تدخل الملاحظة الشهر واليوم والساعة العشرية (التوقيت العالمي hh.ddd) n. ل
العمود ، ن. على الخط. (مقياس كامل 17000 * 12000) ج

3 يذهب غرب تدخل الملاحظة الشهر واليوم والساعة العشرية (التوقيت العالمي hh.ddd) n. ل
العمود ، ن. على الخط. (المقياس الكامل 17000 * 12000)

4 avhrr (PM noaa) أدخل الشهر واليوم والساعة العشرية (التوقيت العالمي hh.ddd) n. ل
العمود (1-2048) ، xlonan ، hna تعطي الطول (xlonan) والجسر
ساعة (hna) عند العقدة الصاعدة عند خط الاستواء

5 avhrr (AM noaa) أدخل الشهر واليوم والساعة العشرية (التوقيت العالمي hh.ddd) n. ل
العمود (1-2048) ، xlonan ، hna تعطي الطول (xlonan) والجسر
ساعة (hna) عند العقدة الصاعدة عند خط الاستواء

6 المهندس (بقعة) أدخل month ، day ، hh.ddd ، long. ، lat. *

7 tm (لاندسات) أدخل month، day، hh.ddd، long.، lat. *

8 etm + (لاندسات 7) أدخل month، day، hh.ddd، long.، lat. *

9 ليس (IRS 1C) أدخل month، day، hh.ddd، long.، lat. *

10 الأسطر النجمية أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

11 أفنير أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

12 ikonos أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

13 RapidEye أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

14 VGT1 (سبوت4) أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

15 VGT2 (سبوت5) أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

16 الرؤية الكونية 2 أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

17 كويك بيرد أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

18 لاندسات 8 أدخل الشهر ، اليوم ، hh.ddd ، long. ، lat. *

* نوت: بالنسبة لتجارب HRV و TM و ETM + و LISS و ASTER ، فإن خطوط الطول والعرض هي
إحداثيات مركز المشهد. يجب أن يكون خط العرض> 0 في نصف الكرة الشمالي و <0 بالنسبة إلى
الجنوب. يجب أن يكون خط الطول> 0 في نصف الكرة الشرقي و <0 في حالة الغرب.

B. جوي نموذج
رمز معنى

0 لا امتصاص غازي

1 استوائي

2 خط الوسط الصيفي

3 شتاء خط الوسط

4 الصيف شبه القطبية

5 شتاء شبه قطبي

6 الولايات المتحدة القياسية 62

7 حدد نموذج الغلاف الجوي الخاص بك كمجموعة مما يلي
5 معلمات لكل قياس: ضغط الارتفاع [km]
[mb] درجة الحرارة [k] كثافة الماء [g / m2] كثافة o3 [g / m3]
على سبيل المثال: يوجد قياس واحد للمسبار اللاسلكي لكل منهما
ارتفاع من 0-25 كم في خطوة 1 كم ، قياس واحد لكل منها
ارتفاع 25-50 كم في خطوة 5 كم ، واثنان واحد
قياسات الارتفاعات 70 كم و 100 كم. هذا يجعل 34
القياسات. في هذه الحالة ، هناك 34 * 5 قيم يجب إدخالها.

8 تحديد نموذج الغلاف الجوي الخاص بك مع تقديم قيم
محتوى بخار الماء والأوزون: uw [g / cm2] uo3 [cm-atm]
الملف الشخصي مأخوذ منا 62.

C. الهباء الجوي نموذج
رمز معنى أيقونة

0 لا الهباء الجوي

1 نموذج قاري

2 النموذج البحري

3 نموذج حضري

4 نموذج شيتل للهباء الصحراوي الخلفية

5 ـ حرق الكتلة الحيوية

6 نموذج الستراتوسفير

7 حدد النموذج الخاص بك أدخل النسبة المئوية الحجمية لكل مكون: c(1) =
٪ حجمي يشبه الغبار c(2) = النسبة المئوية الحجمية للذوبان في الماء
c(3) = النسبة المئوية الحجمية للمحيطات c(4) =٪ الحجمي من السخام الكل
القيم بين 0 و 1.

8 حدد نموذجك الخاص وظيفة توزيع الحجم: سجل متعدد الوسائط عادي (حتى 4
أساليب).

9 تحديد نموذجك الخاص وظيفة توزيع الحجم: جاما المعدلة.

10 حدد نموذجك الخاص وظيفة توزيع الحجم: Junge Power-Law.

11 حدد النموذج الخاص بك قياسات مقياس الضوء الشمسي ، 50 قيمة كحد أقصى ، تم إدخالها على النحو التالي: r
و d V / d (logr) حيث r هو نصف القطر [ميكرون] ، V هو
الحجم ، d V / d (logr) [cm3 / cm2 / micron]. متبوعًا بـ: nr
و n لكل طول موجي حيث nr و n على التوالي
الجزء الحقيقي والخيالي من معامل الانكسار.

D. بخاخ الأيروسول من التركيز نموذج (الرؤية)
إذا كان لديك تقدير لمعلمة رؤية الأرصاد الجوية v ، فأدخل مباشرة إلى
قيمة v [km] (العمق البصري للهباء الجوي (AOD) سيتم حسابه من الهباء الجوي القياسي
حساب تعريفي).

إذا كان لديك تقدير للعمق البصري للهباء الجوي ، فأدخل 0 للرؤية وفي a
أدخل السطر التالي في العمق البصري للهباء عند 550 نانومتر (iaer تعني 'i' للإدخال و
"aerosol") ، على سبيل المثال:
0 - الرؤية
0.112 - العمق البصري للهباء الجوي 550 نانومتر

ملاحظة: إذا كانت قيمة iaer تساوي 0 ، أدخل -1 للرؤية.

E. الهدف ارتفاع (xps) ، مدخل بطاقة الذاكرة : نعم المنصة (إكس بي بي)
الارتفاع المستهدف (xps ، بالسالب [km]): xps> = 0 يعني أن الهدف عند مستوى سطح البحر.
بخلاف ذلك ، تعبر xps عن ارتفاع الهدف (على سبيل المثال ، متوسط ​​الارتفاع) بـ [km] ، معطى
كقيمة سالبة

منصة المستشعر (xpp ، بالسالب [km] أو -1000):
xpp = -1000 تعني أن المستشعر موجود على متن القمر الصناعي.
xpp = 0 يعني أن المستشعر على مستوى الأرض.
-100 <xpp <0 يحدد ارتفاع المستشعر المعبر عنه بـ [km] ؛ هذا الارتفاع
معطى نسبي إلى القادم الهدف الارتفاع كقيمة سالبة.

لمحاكاة الطائرات فقط (xpp لا يساوي 0 ولا يساوي -1000): puw ، po3
(محتوى بخار الماء ومحتوى الأوزون بين الطائرة والسطح)
تايرب (سماكة الهباء الجوي عند 550 نانومتر بين الطائرة والسطح)

في حالة عدم توفر هذه البيانات ، أدخل قيمًا سالبة لها جميعًا. puw ، سوف po3 بعد ذلك
من ملف التعريف القياسي us62 وفقًا للقيم الموجودة على الأرض
مستوى. سيتم حساب taerp وفقًا لملف تعريف أسي يبلغ 2 كم للهباء الجوي.

F. المستشعرات فرقة
هناك احتمالان: إما تحديد الظروف الطيفية الخاصة بك (الرموز -2 ، -1 ، 0 ،
أو 1) أو اختر رمزًا يشير إلى نطاق أحد الأقمار الصناعية المحددة مسبقًا.

حدد الظروف الطيفية الخاصة بك:

رمز معنى

-2 أدخل wlinf ، wlsup. ستكون وظيفة المرشح مساوية لـ 1
على النطاق بأكمله (مثل iwave = 0) ولكن الإخراج خطوة بخطوة
ستتم طباعتها.

-1 أدخل wl (أحادي الشرط ، الامتصاص الغازي مشمول).

0 أدخل wlinf ، wlsup. وظيفة التصفية ستكون مساوية لـ
1 على النطاق بأكمله.

1 أدخل وظائف wlinf و wlsup ووظائف مرشح المستخدم (lambda) بواسطة
خطوة 0.0025 ميكرومتر.

نطاقات القمر الصناعي المحددة مسبقًا:

رمز معنى

2 ميتيوسات النطاق مقابل (0.350-1.110)

3 يذهب الشرق النطاق مقابل (0.490-0.900)

4 يذهب النطاق الغربي مقابل (0.490-0.900)

5 avhrr (نوا6) النطاق 1 (0.550-0.750)

6 avhrr (noaa6) الفرقة 2 (0.690-1.120)

7 avhrr (نوا7) النطاق 1 (0.500-0.800)

8 avhrr (noaa7) الفرقة 2 (0.640-1.170)

9 avhrr (نوا8) النطاق 1 (0.540-1.010)

10 avhrr (noaa8) الفرقة 2 (0.680-1.120)

11 avhrr (نوا9) النطاق 1 (0.530-0.810)

12 avhrr (noaa9) الفرقة 1 (0.680-1.170)

13 avhrr (نوا10) النطاق 1 (0.530-0.780)

14 avhrr (noaa10) الفرقة 2 (0.600-1.190)

15 avhrr (نوا11) النطاق 1 (0.540-0.820)

16 avhrr (noaa11) الفرقة 2 (0.600-1.120)

17 الهريفي1 (بقعة 1) النطاق 1 (0.470-0.650)

18 ساعة 1 (بقعة 1) النطاق 2 (0.600-0.720)

19 ساعة 1 (بقعة 1) النطاق 3 (0.730-0.930)

20 hrv1 (Spot1) عموم الفرقة (0.470-0.790)

21 الهريفي2 (بقعة 1) النطاق 1 (0.470-0.650)

22 ساعة 2 (بقعة 1) النطاق 2 (0.590-0.730)

23 ساعة 2 (بقعة 1) النطاق 3 (0.740-0.940)

24 hrv2 (Spot1) عموم الفرقة (0.470-0.790)

25 tm (لاندسات 5) النطاق 1 (0.430-0.560)

26 tm (لاندسات 5) النطاق 2 (0.500-0.650)

27 tm (لاندسات 5) النطاق 3 (0.580-0.740)

28 tm (لاندسات 5) النطاق 4 (0.730-0.950)

29 tm (لاندسات 5) النطاق 5 (1.5025-1.890)

30 tm (لاندسات 5) النطاق 7 (1.950-2.410)

31 MSS (لاندسات 5) النطاق 1 (0.475-0.640)

32 مللي ثانية (لاندسات 5) الفرقة 2 (0.580-0.750)

33 مللي ثانية (لاندسات 5) الفرقة 3 (0.655-0.855)

34 مللي ثانية (لاندسات 5) الفرقة 4 (0.785-1.100)

35 ماس (ER2) النطاق 1 (0.5025-0.5875)

36 MAS (ER2) الفرقة 2 (0.6075-0.7000)

37 MAS (ER2) الفرقة 3 (0.8300-0.9125)

38 MAS (ER2) الفرقة 4 (0.9000-0.9975)

39 MAS (ER2) الفرقة 5 (1.8200-1.9575)

40 MAS (ER2) الفرقة 6 (2.0950-2.1925)

41 MAS (ER2) الفرقة 7 (3.5800-3.8700)

42 MODIS النطاق 1 (0.6100-0.6850)

43 نطاق MODIS 2 (0.8200-0.9025)

44 نطاق MODIS 3 (0.4500-0.4825)

45 نطاق MODIS 4 (0.5400-0.5700)

46 نطاق MODIS 5 (1.2150-1.2700)

47 نطاق MODIS 6 (1.6000-1.6650)

48 نطاق MODIS 7 (2.0575-2.1825)

49 avhrr (نوا12) النطاق 1 (0.500-1.000)

50 avhrr (noaa12) الفرقة 2 (0.650-1.120)

51 avhrr (نوا14) النطاق 1 (0.500-1.110)

52 avhrr (noaa14) الفرقة 2 (0.680-1.100)

53 بولدر النطاق 1 (0.4125-0.4775)

54 POLDER band 2 (غير قطبي) (0.4100-0.5225)

55 POLDER band 3 (غير قطبي) (0.5325-0.5950)

56 نطاق بولدر 4 P1 (0.6300-0.7025)

57 POLDER band 5 (غير قطبي) (0.7450-0.7800)

58 POLDER band 6 (غير قطبي) (0.7000-0.8300)

59 نطاق بولدر 7 P1 (0.8100-0.9200)

60 POLDER band 8 (غير قطبي) (0.8650-0.9400)

61 etm + (لاندسات 7) النطاق 1 (0.435-0.520)

62 etm + (لاندسات 7) النطاق 2 (0.506-0.621)

63 etm + (لاندسات 7) النطاق 3 (0.622-0.702)

64 etm + (لاندسات 7) النطاق 4 (0.751-0.911)

65 etm + (لاندسات 7) النطاق 5 (1.512-1.792)

66 etm + (لاندسات 7) النطاق 7 (2.020-2.380)

67 etm + (لاندسات 7) النطاق 8 (0.504-0.909)

68 ليس (آي آر سي 1C) النطاق 2 (0.502-0.620)

69 liss (IRC 1C) النطاق 3 (0.612-0.700)

70 liss (IRC 1C) النطاق 4 (0.752-0.880)

71 liss (IRC 1C) النطاق 5 (1.452-1.760)

72 الأسطر النجمية النطاق 1 (0.480-0.645)

73 أستر باند 2 (0.588-0.733)

74 أستر باند 3N (0.723-0.913)

75 أستر باند 4 (1.530-1.750)

76 أستر باند 5 (2.103-2.285)

77 أستر باند 6 (2.105-2.298)

78 أستر باند 7 (2.200-2.393)

79 أستر باند 8 (2.248-2.475)

80 أستر باند 9 (2.295-2.538)

81 أفنير النطاق 1 (0.390-0.550)

82 فرقة أفنير 2 (0.485-0.695)

83 فرقة أفنير 3 (0.545-0.745)

84 فرقة أفنير 4 (0.700-0.925)

85 ikonos الشريط الأخضر (0.350-1.035)

86 شريط إيكونوس أحمر (0.350-1.035)

87 نطاق إيكونوس نير (0.350-1.035)

88 RapidEye الشريط الأزرق (0.438-0.513)

89 شريط RapidEye الأخضر (0.463-0.594)

90 شريط أحمر RapidEye (0.624-0.690)

91 فرقة RapidEye RedEdge (0.500-0.737)

92 نطاق RapidEye NIR (0.520-0.862)

93 VGT1 (سبوت4) النطاق 0 (0.400-0.500)

94 VGT1 (SPOT4) النطاق 2 (0.580-0.782)

95 VGT1 (SPOT4) النطاق 3 (0.700-1.030)

96 VGT1 (SPOT4) نطاق MIR (1.450-1.800)

97 VGT2 (سبوت5) النطاق 0 (0.400-0.550)

98 VGT2 (SPOT5) النطاق 2 (0.580-0.780)

99 VGT2 (SPOT5) النطاق 3 (0.700-1.000)

100 VGT2 (SPOT5) نطاق MIR (1.450-1.800)

101 فرقة WorldView 2 Panchromatic (0.447-0.808)

102 الشريط الأزرق الساحلي WorldView 2 (0.396-0.458)

103 WorldView 2 النطاق الأزرق (0.442-0.515)

104 WorldView 2 النطاق الأخضر (0.506-0.586)

105 شريط أصفر WorldView 2 (0.584-0.632)

106 WorldView 2 شريط أحمر (0.624-0.694)

107 WorldView 2 نطاق الحافة الحمراء (0.699-0.749)

108 نطاق WorldView 2 NIR1 (0.765-0.901)

109 نطاق WorldView 2 NIR2 (0.856-0.1043)

110 كويك بيرد فرقة بانكروماتيك (0.405-1.053)

111 فرقة كويك بيرد الزرقاء (0.430-0.545)

112 شريط QuickBird Green (0.466-0.620)

113 شريط كويك بيرد الأحمر (0.590-0.710)

114 نطاق QuickBird NIR1 (0.715-0.918)

115 لاندسات 8 شريط الهباء الساحلي (0.427 نانومتر - 0.459 نانومتر)

116 Landsat 8 Blue Band (436 نانومتر - 527 نانومتر)

117 لاندسات 8 جرين باند (512 نانومتر - 610 نانومتر)

118 لاندسات 8 ريد باند (625 نانومتر - 691 نانومتر)

119 نطاق لاندسات 8 بانكروماتيك (488 نانومتر - 692 نانومتر)

نطاق 120 Landsat 8 NIR (829nm-900nm)

121 فرقة لاندسات 8 سيروس (1340 نانومتر -1409 XNUMX نانومتر)

122 فرقة لاندسات 8 SWIR1 (1515 نانومتر - 1697 نانومتر)

123 فرقة لاندسات 8 SWIR2 (2037 نانومتر - 2355 نانومتر)

أمثلة

جوي تصحيح of a اندسات 7 قناة
يعتمد المثال على مجموعة بيانات نموذجية في ولاية كارولينا الشمالية (توقيت غرينتش -5 ساعات). أولا وضعنا
المنطقة الحسابية لخريطة القمر الصناعي ، على سبيل المثال القناة 4:
ز المنطقة النقطية = lsat7_2002_40 -p
من المهم التحقق من البيانات الوصفية المتاحة لموقع الشمس الذي يجب أن يكون
المحددة لتصحيح الغلاف الجوي. أحد الخيارات هو التحقق من وقت تجاوز القمر الصناعي
مع موقع الشمس كما ورد في البيانات الوصفية. بالنسبة لمجموعة بيانات نموذج ولاية كارولينا الشمالية ، فهم
تم تخزينها أيضًا لكل قناة ويمكن استرجاعها على النحو التالي:
r.info lsat7_2002_40
في هذه الحالة ، لدينا: SUN_AZIMUTH = 120.8810347 ، SUN_ELEVATION = 64.7730999.

إذا كانت البيانات الوصفية لموقع الشمس غير متوفرة ، فيمكننا أيضًا حسابها من الممر العلوي
الوقت على النحو التالي (ص قناع الشمس يستخدم SOLPOS):
r.sunmask-height = الارتفاع للخارج = السنة الوهمية = 2002 شهر = 5 أيام = 24 ساعة = 10 دقائق = 42 ثانية = 7 منطقة زمنية = -5
# .. التقارير: سمت الشمس: 121.342461 ، زاوية الشمس فوق هورتس (تصحيح الانكسار): 65.396652
إذا كان وقت العبور غير معروف ، فاستخدم توقع تجاوز القمر الصناعي.

تحويل DN (الرقم الرقمي = قيم البكسل) إلى إشعاع في أعلى الغلاف الجوي (TOA) ، باستخدام
الصيغة
Lλ = ((LMAXλ - LMINλ) / (QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ
أين:

· Lλ = الإشعاع الطيفي عند فتحة المستشعر بالواط / (متر مربع * ستير
* µm) ، الإشعاع الظاهر كما يراه جهاز استشعار الأقمار الصناعية ؛

· QCAL = قيمة البكسل المُعايرة الكمية في DN ؛

· LMINλ = الإشعاع الطيفي الذي يتم تحجيمه إلى QCALMIN بالواط / (متر
مربع * ستير * ميكرومتر) ؛

· LMAXλ = الإشعاع الطيفي الذي يتم قياسه إلى QCALMAX بالواط / (متر
مربع * ستير * ميكرومتر) ؛

· QCALMIN = الحد الأدنى لقيمة البكسل المُعايرة الكمية (المقابلة لـ
LMINλ) في DN ؛

· QCALMAX = الحد الأقصى لقيمة بكسل المعايرة الكمية (المقابلة لـ
LMAXλ) في DN = 255.
LMINλ و LMAXλ هما الإشعاع المرتبط بالحد الأدنى والأقصى لقيمة DN ،
ويتم الإبلاغ عنها في ملف البيانات الوصفية لكل صورة ، أو في الجدول 1. مكاسب عالية أو منخفضة
يتم أيضًا الإبلاغ عن الكسب في ملف البيانات الوصفية لكل صورة من صور القمر الصناعي. الحد الأدنى لقيمة DN
(QCALMIN) هو 1 لصور Landsat ETM + (انظر كتيب Landsat ، انظر الفصل 11) ، و
الحد الأقصى لقيمة DN (QCALMAX) هو 255. QCAL هي قيمة DN لكل بكسل منفصل في
صورة لاندسات.

نستخرج المعاملات ونطبقها للحصول على خريطة التألق:
CHAN = 4
r.info lsat7_2002 _ $ {CHAN} 0 -h | tr '\ n' '' | sed 's + ++ g' | tr ':' '\ n' | grep "LMIN_BAND $ {CHAN} \ | LMAX_BAND $ {CHAN}"
LMAX_BAND4=241.100,p016r035_7x20020524.met
LMIN_BAND4=-5.100,p016r035_7x20020524.met
QCALMAX_BAND4=255.0,p016r035_7x20020524.met
QCALMIN_BAND4=1.0,p016r035_7x20020524.met
التحويل إلى الإشعاع (يتم إجراء هذا الحساب للنطاق 4 ، والنطاقات الأخرى ، و
يجب استبدال الأرقام المكتوبة بخط مائل بالقيم المرتبطة بها):
r.mapcalc "lsat7_2002_40_rad = ((241.1 - (-5.1)) / (255.0 - 1.0)) * (lsat7_2002_40 - 1.0) + (-5.1)"
# البحث عن متوسط ​​الارتفاع (الهدف فوق مستوى سطح البحر ، يُستخدم كقيمة تهيئة في ملف التحكم)
r.univar الارتفاع
قم بإنشاء ملف تحكم "icnd.txt" للقناة 4 (NIR) ، بناءً على البيانات الوصفية. للجسر
الوقت ، نحتاج إلى تحديد الساعات العشرية:
10:42:07 بالتوقيت المحلي NC = 10.70 ساعة عشرية (الدقائق العشرية: 42 * 100/60) وهو
15.70 بتوقيت جرينتش:
8 - الشروط الهندسية = Landsat ETM +
5 24 15.70 -78.691 35.749 - شهر يوم hh.ddd خط عرض ("hh.ddd" بالساعات العشرية بتوقيت جرينتش)
2 - وضع الغلاف الجوي = خط العرض المتوسط ​​في الصيف
1- نموذج الهباء الجوي = قاري
50 - الرؤية [كم] (تركيز نموذج الهباء الجوي)
-0.110 - متوسط ​​الارتفاع المستهدف فوق مستوى سطح البحر [كم]
-1000 - حساس على متن القمر الصناعي
64 - النطاق الرابع من ETM + Landsat 4
أخيرًا ، قم بتشغيل التصحيح الجوي (-r لخريطة إدخال الانعكاس ؛ -a للتاريخ> يوليو
2000):
i.atcorr -r -a lsat7_2002_40_rad Elev = معلمات الارتفاع = إخراج icnd_lsat4.txt = lsat7_2002_40_atcorr
لاحظ أنه تتم قراءة قيمة الارتفاع من ملف "icnd_lsat4.txt" في البداية إلى
حساب التحويل الأولي. من الضروري إعطاء قيمة يمكن أن تكون هي الوسيلة
قيمة نموذج الارتفاع. لتصحيح الغلاف الجوي ثم الارتفاع النقطي
يتم استخدام القيم من الخريطة.

لاحظ أن العملية حسابية مكثفة.
لاحظ أيضًا أن أنا أتكور تقارير زاوية ارتفاع الشمس فوق الأفق بدلاً من الشمس
زاوية الذروة.

متبقي وثائق قضايا

1. يجب شرح تأثير وأهمية قيمة الرؤية أو الخريطة أيضًا
كيفية الحصول على تقدير للرؤية أو العمق البصري للهباء عند 550 نانومتر.

استخدم i.atcorrgrass عبر الإنترنت باستخدام خدمات onworks.net