أكثر

ناتج نموذج تصحيح الغلاف الجوي لصورة لاندسات 8؟

ناتج نموذج تصحيح الغلاف الجوي لصورة لاندسات 8؟


فيما يلي الإجراء والملاحظة والسؤال / الطلب:

مصدر البيانات: Earth Explorer Landsat 8 Grass GIS الإصدار: 7.0.2 RC1 Windows 7 32 bit OS

لقد كنت أحاول إجراء تصحيح الغلاف الجوي باستخدام مرجع Grass wiki

https://grasswiki.osgeo.org/wiki/Atmospheric_correction

  • الخطوة 01:

    لقد تخطيت استخدام أداة Toar واستخدمت حاسبة الخريطة النقطية لإنشاء:

    أ).الانعكاس = REFLECTANCE_MULT_BAND_X (0.00002) * DN + REFLECTANCE_ADD_BAND_X (-0.1)

    ب).الانعكاس مع تصحيح زاوية الشمس = الانعكاس / الخطيئة (Ɵ)

    هنا الخطيئة (Ɵ) = الخطيئة (49.52284282) = 0.760664796730206

  • الخطوة 02:

    تم توصيل نطاقات الانعكاس التي تم إنشاؤها في الخطوة 01 للنطاق 4 و 3 و 2 amd 8 مع ملف التحكم الخاص بها في i.atcorr.

  • الخطوة 03: نطاق IHS Pan Sharpened 4،3،2 مع النطاق 8 المتولد في الخطوة 02

  • الخطوة 04: يتم إنشاء مركب RGB باستخدام Red = band 4 ، و Green = band 3 ، و Blue = band 2. (جميع النطاقات مشتقة في الخطوة 03)

أرفق الصورة الناتجة للمقارنة. من الواضح أنني فاتني شيء. يجب أن تكون الصورة المعدلة جوًا أكثر وضوحًا. إنها لا تبدو حتى وكأنها صورة RGB.

كما يتم إرفاق ملف البيانات الوصفية ، ومخرجات الأوامر (بحيث يمكنك اتباع خطواتي ، لملاحظة ما إذا كان هناك أي شيء غير عادي) وملف تحكم للرجوع إليه. لقد استخدمت نموذج الغلاف الجوي الاستوائي ونموذج الهباء الجوي الحضري. تم جمع القيم المستخدمة في ملف التحكم من ويكي Atmospheric correction.

هل يمكن لأي شخص أن يخبرني كيف يمكنني الحصول على صورة صحيحة معدلة للغلاف الجوي؟

كما أنني لاحظت أنه لم يعد هناك أي امتداد i.landsat.dehaze متاحًا. هل يمكن أن تخبرني من فضلك كيف قد يكون التجفيف ممكنًا في غياب التمديد؟

ملف البيانات الوصفية ================================================= ================================ GROUP = L1_METADATA_FILE GROUP = METADATA_FILE_INFO ORIGIN = "الصورة مجاملة من هيئة المسح الجيولوجي الأمريكية" REQUEST_ID = "0501510272804_00036" LANDSAT_SCENE_ID = "LC81370442015300LGN00" FILE_DATE = 2015-10-27T08: 52: 24Z STATION_ID = "LGN" PROCESSING_SOFTWARE_VERSION = "LPGS_2.5.1" END_GROUP = METADATA_FILE_INFO مجموعة = PRODUCT_METADATA DATA_TYPE = "L1T" ELEVATION_SOURCE = "GLS2000" OUTPUT_FORMAT = "GEOTIFF" SPACECRAFT_ID = "LANDSAT_8" SENSOR_ID = "OLI_TIRS" WRS_PATH = 137 WRS_ROW = 44 NADIR_OFFNADIR = "NADIR" TARGET_WRS_PATH = 137 TARGET_WRS_ROWEN = 44.085PRODATE_ROWEN = 44 Duration_WRS_ROW = 44 CORNER_UL_LON_PRODUCT = 89.24071 CORNER_UR_LAT_PRODUCT = 24.18946 CORNER_UR_LON_PRODUCT = 91.51825 CORNER_LL_LAT_PRODUCT = 22.01936 CORNER_LL_LON_PRODUCT = 89.29959 CORNER_LR_LAT_RODUCT = 91.29959 CORNER_LR_LAT_RODUCT = 91.29959 CORNER_LR_LAT_RUCT_PRODUCT = 91.51825 ON_X_PRODUCT = 117900.000 CORNER_UL_PROJECTION_Y_PRODUCT = 2676000.000 CORNER_UR_PROJECTION_X_PRODUCT = 349500.000 CORNER_UR_PROJECTION_Y_PRODUCT = 2676000.000 CORNER_LL_PROJECTION_X_PRODUCT = 117900.000 CORNER_LL_PROJECTION_Y_PRODUCT = 2439600.000 CORNER_LR_PROJECTION_X_PRODUCT = 349500.000 CORNER_LR_PROJECTION_Y_PRODUCT = 2439600.000 PANCHROMATIC_LINES = 15761 PANCHROMATIC_SAMPLES = 15441 REFLECTIVE_LINES = 7881 = 7721 REFLECTIVE_SAMPLES THERMAL_LINES = 7881 = 7721 THERMAL_SAMPLES FILE_NAME_BAND_1 = "LC81370442015300LGN00_B1.TIF" FILE_NAME_BAND_2 = "LC81370442015300LGN00_B2.TIF" FILE_NAME_BAND_3 = "LC81370442015300LGN00_B3.TIF" FILE_NAME_BAND_4 = "LC81370442015300LGN00_B4.TIF" FILE_NAME_BAND_5 = "LC81370442015300LGN00_B5.TIF" FILE_NAME_BAND_6 = "LC81370442015300LGN00_B6.TIF" FILE_NAME_BAND_7 = "LC81370442015300LGN00_B7.TIF" FILE_NAME_BAND_8 = "LC81370442015300LGN00_B8.TIF" FILE_NAME_BAND_9 = "LC81370442015300LGN00_B9.TIF" FILE_NAME_BAND_10 = "LC81370442015300LGN00_B10.TIF" FILE_NAME_ BAND_11 = "LC81370442015300LGN00_B11.TIF" FILE_NAME_BAND_QUALITY = "LC81370442015300LGN00_BQA.TIF" METADATA_FILE_NAME = "LC81370442015300LGN00_MTL.txt" BPF_NAME_OLI = "LO8BPF20151027040956_20151027054850.01" BPF_NAME_TIRS = "LT8BPF20151027040602_20151027054456.01" CPF_NAME = "L8CPF20151001_20151231.01" RLUT_FILE_NAME = "L8RLUT20150303_20431231v11.h5" END_GROUP = PRODUCT_METADATA مجموعة = IMAGE_ATTRIBUTES CLOUD_COVER = 4.67 = 4.75 CLOUD_COVER_LAND IMAGE_QUALITY_OLI = 9 = 9 IMAGE_QUALITY_TIRS TIRS_SSM_POSITION_STATUS = "الاسمية" ROLL_ANGLE = -0.001 SUN_AZIMUTH = 150،12621205 SUN_ELEVATION = 49.52284282 EARTH_SUN_DISTANCE = 0.9939432 GROUND_CONTROL_POINTS_VERSION = 2 = 298 GROUND_CONTROL_POINTS_MODEL GEOMETRIC_RMSE_MODEL = 6.823 GEOMETRIC_RMSE_MODEL_Y = 4.812 GEOMETRIC_RMSE_MODEL_X = 4.837 GROUND_CONTROL_POINTS_VERIFY = 88 GEOMETRIC_RMSE_VERIFY = 4.846 END_GROUP = IMAGE_ATTRIBUTES GROUP = MIN_MAX_RADIANCE RADIANCE_MAXIMUM_BAND_1 = 769.35406 RADIANCE_MINIMUM_BAND_1 =- MUM_BAND_2 = 787.82794 RADIANCE_MINIMUM_BAND_2 = -65،05909 RADIANCE_MAXIMUM_BAND_3 = 725.97668 RADIANCE_MINIMUM_BAND_3 = -59،95139 RADIANCE_MAXIMUM_BAND_4 = 612.18445 RADIANCE_MINIMUM_BAND_4 = -50،55439 RADIANCE_MAXIMUM_BAND_5 = 374.62637 RADIANCE_MINIMUM_BAND_5 = -30،93677 RADIANCE_MAXIMUM_BAND_6 = 93.16617 RADIANCE_MINIMUM_BAND_6 = -7،69369 RADIANCE_MAXIMUM_BAND_7 = 31.40199 RADIANCE_MINIMUM_BAND_7 = -2،59319 RADIANCE_MAXIMUM_BAND_8 = 692.82434 RADIANCE_MINIMUM_BAND_8 = -57،21366 RADIANCE_MAXIMUM_BAND_9 = 146.41248 RADIANCE_MINIMUM_BAND_9 = -12،09079 RADIANCE_MAXIMUM_BAND_10 = 22.00180 RADIANCE_MINIMUM_BAND_10 = 0.10033 RADIANCE_MAXIMUM_BAND_11 = 22.00180 RADIANCE_MINIMUM_BAND_11 = 0.10033 END_GROUP = MIN_MAX_RADIANCE مجموعة = MIN_MAX_REFLECTANCE REFLECTANCE_MAXIMUM_BAND_1 = 1.210700 REFLECTANCE_MINIMUM_BAND_1 = -0.099980 REFLECTANCE_MAXIMUM_BAND_2 = 1.210700 REFLECTANCE_MINIMUM_BAND_2 = -0.099980 REFLECTANCE_MAXIMUM_BAND_3 = 1.210700 REFLECTANCE_MINIMUM_BAND_3 = -0.099980 الانعكاس _MAXIMUM_BAND_4 = 1.210700 REFLECTANCE_MINIMUM_BAND_4 = -0.099980 REFLECTANCE_MAXIMUM_BAND_5 = 1.210700 REFLECTANCE_MINIMUM_BAND_5 = -0.099980 REFLECTANCE_MAXIMUM_BAND_6 = 1.210700 REFLECTANCE_MINIMUM_BAND_6 = -0.099980 REFLECTANCE_MAXIMUM_BAND_7 = 1.210700 REFLECTANCE_MINIMUM_BAND_7 = -0.099980 REFLECTANCE_MAXIMUM_BAND_8 = 1.210700 REFLECTANCE_MINIMUM_BAND_8 = -0.099980 REFLECTANCE_MAXIMUM_BAND_9 = 1.210700 REFLECTANCE_MINIMUM_BAND_9 = -0.099980 END_GROUP = MIN_MAX_REFLECTANCE مجموعة = MIN_MAX_PIXEL_VALUE QUANTIZE_CAL_MAX_BAND_1 = 65535 QUANTIZE_CAL_MIN_BAND_1 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_2 QUANTIZE_CAL_MIN_BAND_2 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_3 QUANTIZE_CAL_MIN_BAND_3 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_4 QUANTIZE_CAL_MIN_BAND_4 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_5 QUANTIZE_CAL_MIN_BAND_5 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_6 QUANTIZE_CAL_MIN_BAND_6 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_7 QUANTIZE_CAL_MIN_BAND_7 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_8 QUANTIZE_CAL_MIN_BAND_8 = 1 QUANTIZE_ CAL_MAX_BAND_9 = 65535 QUANTIZE_CAL_MIN_BAND_9 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_10 QUANTIZE_CAL_MIN_BAND_10 = 1 = 65535 QUANTIZE_CAL_MAX_BAND_11 QUANTIZE_CAL_MIN_BAND_11 = 1 = END_GROUP MIN_MAX_PIXEL_VALUE مجموعة = RADIOMETRIC_RESCALING RADIANCE_MULT_BAND_1 = 1.2709E-02 RADIANCE_MULT_BAND_2 = 1.3014E-02 RADIANCE_MULT_BAND_3 = 1.1993E-02 RADIANCE_MULT_BAND_4 = 1.0113E-02 RADIANCE_MULT_BAND_5 = 6.1886E-03 RADIANCE_MULT_BAND_6 = 1.5390E-03 RADIANCE_MULT_BAND_7 = 5.1874E-04 RADIANCE_MULT_BAND_8 = 1.1445E-02 RADIANCE_MULT_BAND_9 = 2.4186E-03 RADIANCE_MULT_BAND_10 = 3.3420E-04 RADIANCE_MULT_BAND_11 = 3.3420E-04 RADIANCE_ADD_BAND_1 = -63،54622 RADIANCE_ADD_BAND_2 = -65،07210 RADIANCE_ADD_BAND_3 = -59،96339 RADIANCE_ADD_BAND_4 = -50،56450 RADIANCE_ADD_BAND_5 = -30،94296 RADIANCE_ADD_BAND_6 = -7،69523 RADIANCE_ADD_BAND_7 = -2،59370 RADIANCE_ADD_BAND_8 = -57،22510 RADIANCE_ADD_BAND_9 = -12،09321 RADIANCE_ADD_BAND_10 = 0.10000 RADIANCE_ADD_BAND_11 = 0.10000 REFLECTANCE_MULT_BAND_1 = 2.0000E-05 RE FLECTANCE_MULT_BAND_2 = 2.0000E-05 REFLECTANCE_MULT_BAND_3 = 2.0000E-05 REFLECTANCE_MULT_BAND_4 = 2.0000E-05 REFLECTANCE_MULT_BAND_5 = 2.0000E-05 REFLECTANCE_MULT_BAND_6 = 2.0000E-05 REFLECTANCE_MULT_BAND_7 = 2.0000E-05 REFLECTANCE_MULT_BAND_8 = 2.0000E-05 REFLECTANCE_MULT_BAND_9 = 2.0000E-05 REFLECTANCE_ADD_BAND_1 = -،100000 REFLECTANCE_ADD_BAND_2 = -0.100000 REFLECTANCE_ADD_BAND_3 = -0.100000 REFLECTANCE_ADD_BAND_4 = -0.100000 REFLECTANCE_ADD_BAND_5 = -0.100000 REFLECTANCE_ADD_BAND_6 = -0.100000 REFLECTANCE_ADD_BAND_7 = -0.100000 REFLECTANCE_ADD_BAND_8 = -0.100000 REFLECTANCE_ADD_BAND_9 = -0.100000 END_GROUP = RADIOMETRIC_RESCALING مجموعة = TIRS_THERMAL_CONSTANTS K1_CONSTANT_BAND_10 = 774.8853 K1_CONSTANT_BAND_11 = 480.8883 K2_CONSTANT_BAND_10 = 1321،0789 K2_CONSTANT_BAND_11 = 1201.1442 END_GROUP = TIRS_THERMAL_CONSTANTS GROUP = PROJECTION_PARAMETERS MAP_PROJECTION = "UTM" DATUM = "WGS84" ELLIPSOID = "WGS84" UTM_ZONE = 46 GRID_CELL_SIZE_PANCHROMIDELL = 30.00 GRID_CELL_SIZE_THERMAL = 30.00 ORIENTATION = "NORTH_UP" RESAMPLING_OPTION = "CUBIC_CONVOLUTION" END_GROUP = PROJECTION_PARAMETERS END_GROUP = L1_METADATA_FILE END Command Output Log ===================== ==================================================== ===== شريط تصحيح الغلاف الجوي 4 (الجمعة 13 نوفمبر 23:54:00 2015) i.atcorr -r - إدخال زائد = TOAR_B4 @ معلمات TOAR = E:  LandCover_Classification  Aeronet  Atcorr_B4.txt الإخراج = ATM_B4 * * ***************************** 6S الإصدار 4.2b ***************** ************* * * هوية الشروط الهندسية * * ------------------------------ - * * ملاحظة لاندسات 8 * * * * شهر: 10 يوم: 27 * * زاوية ذروة الشمس: 40.69 درجة زاوية سمت شمسي: 150.19 درجة * * زاوية عرض ذروة: 0.00 درجة زاوية سمتي: 0.00 درجة * * زاوية تشتت: 139.31 فرق الزاوية السمتي درجة: 150.19 درجة * * * * وصف نموذج الغلاف الجوي * * ----------------------------- * * هوية نموذج الغلاف الجوي * * استوائي (uh2o = 4.12g / cm2 ، uo3 = .247cm-atm ) * * * * هوية نوع الهباء الجوي: * * نموذج الهباء الجوي الحضري * * * * هوية الحالة البصرية: * * الرؤية: 7.54 كم. سميك. 550 نانومتر: 0.5500 * * * * حالة طيفية * * ------------------ * * لاندسات 8 أحمر * * قيمة وظيفة المرشح: * * wl inf = 0.626 mic wl sup = 0.682 ميكروفون * * * * نوع الهدف * * ----------- * * أرضية متجانسة * * انعكاس ثابت على الأطياف 0.000 * * * * وصف الارتفاع المستهدف * * ------- --------------------- * * ضغط الأرض [mb] 1012.54 * * ارتفاع الأرض [كم] 0.004 * * المحتوى الغازي على مستوى الهدف: * * uh2o = 4.110 g / cm2 uo3 = 0.248 cm-atm * * * * تم تفعيل التصحيح الجوي * * -------------------------------- * تصحيح الغلاف الجوي ... اكتمل تصحيح الغلاف الجوي. (الجمعة 13 نوفمبر 23:56:08) انتهى الأمر (دقيقتان و 8 ثوانٍ) نطاق تصحيح الغلاف الجوي 3 (السبت 14 نوفمبر 00:04:55 2015) i.atcorr -r -verbose input = TOAR_B3 @ TOAR parameters = E :  LandCover_Classification  Aeronet  Atcorr_B3.txt الإخراج = ATM_B3 * ******************************** الإصدار 4.2b *** *************************** * * هوية الشروط الهندسية * * ---------------- --------------- * * ملاحظة لاندسات 8 * * * * شهر: 10 يوم: 27 * * زاوية ذروة الشمس: 40.69 درجة زاوية سمتي شمسي: 150.19 درجة * * زاوية عرض ذروة: 0.00 درجة زاوية سمتي: 0.00 درجة * * زاوية التشتت: 139.31 درجة فرق زاوية السمتي: 150.19 درجة * * * * وصف نموذج الغلاف الجوي * * -------------------- --------- * * هوية النموذج الجوي: * * استوائي (uh2o = 4.12g / cm2، uo3 = .247cm-atm) * * * * هوية نوع الهباء الجوي: * * نموذج الهباء الجوي الحضري * * * * هوية الحالة البصرية: * * الرؤية: 7.54 كيلو متر. سميك. 550 نانومتر: 0.5500 * * * * حالة طيفية * * ------------------ * * لاندسات 8 أخضر * * قيمة وظيفة المرشح: * * wl inf = 0.513 mic wl sup = 0.600 ميكروفون * * * * نوع الهدف * * ----------- * * أرضية متجانسة * * انعكاس ثابت على الأطياف 0.000 * * * * وصف الارتفاع المستهدف * * ------- --------------------- * * ضغط الأرض [mb] 1012.54 * * ارتفاع الأرض [كم] 0.004 * * المحتوى الغازي على مستوى الهدف: * * uh2o = 4.110 g / cm2 uo3 = 0.248 cm-atm * * * * تم تفعيل التصحيح الجوي * * -------------------------------- * تصحيح الغلاف الجوي ... اكتمل تصحيح الغلاف الجوي. (السبت 14 نوفمبر 00:06:28 2015) انتهى الأمر (1 دقيقة و 32 ثانية) نطاق تصحيح الغلاف الجوي 2 (السبت 14 نوفمبر 00:12:59 2015) i.atcorr -r --verbose input = TOAR_B2 @ TOAR parameters = E :  LandCover_Classification  Aeronet  Atcorr_B2.txt output = ATM_B2 * ******************************** الإصدار 4.2b *** *************************** * * هوية الشروط الهندسية * * ---------------- --------------- * * مراقبة لاندسات 8 * * * * شهر: 10 يوم: 27 * * زاوية ذروة الشمس: 40.69 درجة زاوية سمتي شمسي: 150.19 درجة * * زاوية عرض ذروة: 0.00 درجة زاوية سمتي: 0.00 درجة * * زاوية التشتت: 139.31 درجة فرق زاوية السمتي: 150.19 درجة * * * * وصف نموذج الغلاف الجوي * * -------------------- --------- * * هوية النموذج الجوي: * * استوائي (uh2o = 4.12g / cm2، uo3 = .247cm-atm) * * * * هوية نوع الهباء الجوي: * * نموذج الهباء الجوي الحضري * * * * هوية الحالة البصرية: * * الرؤية: 7.54 كيلو متر. سميك. 550 نانومتر: 0.5500 * * * * حالة طيفية * * ------------------ * * لاندسات 8 أزرق * * قيمة وظيفة المرشح: * * wl inf = 0.436 mic wl sup = 0.527 ميكروفون * * * * نوع الهدف * * ----------- * * أرضية متجانسة * * انعكاس ثابت على الأطياف 0.000 * * * * وصف الارتفاع المستهدف * * ------- --------------------- * * ضغط الأرض [mb] 1012.54 * * ارتفاع الأرض [كم] 0.004 * * المحتوى الغازي على مستوى الهدف: * * uh2o = 4.110 g / cm2 uo3 = 0.248 cm-atm * * * * تم تفعيل التصحيح الجوي * * -------------------------------- * تصحيح الغلاف الجوي ... اكتمل تصحيح الغلاف الجوي. (السبت 14 نوفمبر 00:14:32 2015) انتهى الأمر (1 دقيقة و 32 ثانية) شريط تصحيح الغلاف الجوي 8 (السبت 14 نوفمبر 00:17:58 2015) i.atcorr -r - الإدخال المفرط = TOAR_B8 @ معلمات TOAR = E :  LandCover_Classification  Aeronet  Atcorr_B8.txt output = ATM_B8 * ******************************** الإصدار 4.2b *** *************************** * * هوية الشروط الهندسية * * ---------------- --------------- * * مراقبة لاندسات 8 * * * * شهر: 10 يوم: 27 * * زاوية ذروة الشمس: 40.69 درجة زاوية سمتي شمسي: 150.19 درجة * * زاوية عرض ذروة: 0.00 درجة زاوية سمتي الرؤية: 0.00 درجة * * زاوية التشتت: 139.31 درجة فرق زاوية السمت: 150.19 درجة * * * * وصف نموذج الغلاف الجوي * * -------------------- --------- * * هوية النموذج الجوي: * * استوائي (uh2o = 4.12g / cm2، uo3 = .247cm-atm) * * * * هوية نوع الهباء الجوي: * * نموذج الهباء الجوي الحضري * * * * هوية الحالة البصرية: * * الرؤية: 7.54 كيلو متر. سميك. 550 نانومتر: 0.5500 * * * * حالة طيفية * * ------------------ * * لاندسات 8 بانكروماتيك * * قيمة وظيفة المرشح: * * wl inf = 0.488 mic wl sup = 0.692 ميكروفون * * * * نوع الهدف * * ----------- * * أرضية متجانسة * * انعكاس ثابت على الأطياف 0.000 * * * * وصف الارتفاع المستهدف * * ------- --------------------- * * ضغط الأرض [mb] 1012.54 * * ارتفاع الأرض [كم] 0.004 * * المحتوى الغازي على مستوى الهدف: * * uh2o = 4.110 g / cm2 uo3 = 0.248 cm-atm * * * * تم تفعيل التصحيح الجوي * * -------------------------------- * تصحيح الغلاف الجوي ... اكتمل تصحيح الغلاف الجوي. (السبت 14 نوفمبر 00:19:53 2015) انتهى الأمر (دقيقة و 54 ثانية) شحذ عموم (السبت 14 نوفمبر 00:23:58 2015) i.pansharpen - أحمر زائد = ATM_B4 @ TOAR أخضر = ATM_B3 @ TOAR blue = ATM_B2 @ TOAR pan = ATM_B8 @ TOAR output = PAN_RGB method = ihs تنفيذ شحذ عموم مع صورة عالية الدقة: 15.000000 باستخدام IHS <-> خوارزمية RGB التحويل إلى مساحة ألوان IHS ... مطابقة الرسم البياني ... التحويل مرة أخرى إلى مساحة ألوان RGB والشحذ ... التعيين جداول ألوان معادلة باللون الرمادي لإخراج الصور ... تم إنشاء خرائط الإخراج التالية التي تم شحذها عمومًا: PAN_RGB_red PAN_RGB_green PAN_RGB_blue لتصور الإخراج ، قم بتشغيل: g.region -p raster = PAN_RGB_red d.rgb r = PAN_RGB_red g = PAN_RGB_green b = PAN_RGB ، اجمع القنوات في خريطة RGB واحدة باستخدام "r.composite". يمكن إعادة موازنة ألوان القناة باستخدام i.colors.enhance. مجموعة  غير موجود بعد. جارٍ الإنشاء… إضافة الخرائط النقطية إلى المجموعة … إضافة خريطة نقطية  لتجميع إضافة خريطة نقطية  لتجميع إضافة خريطة نقطية  للمجموعة (السبت 14 نوفمبر 00:26:31 2015) انتهى الأمر (2 دقيقة و 32 ثانية) مركب (السبت 14 نوفمبر 00:31:31 2015) r.composite - أحمر زائد = PAN_RGB_red @ TOAR green = PAN_RGB_green @ TOAR blue = PAN_RGB_blue @ TOAR output = RGB_PAN.tif إنشاء جدول ألوان لخريطة نقطية للإخراج ... كتابة خريطة نقطية … r.composite كاملة. خريطة نقطية  خلقت. (السبت 14 تشرين الثاني (نوفمبر) 00:32:04 2015) انتهى الأمر (33 ثانية) نموذج ملف التحكم للنطاق 2 =========================== ==================================================== === 18 10 27 4.41 90.373 23.113 1 3 0 0.55 -0.004-1000116

أود أن أقترح اتخاذ خطوة واحدة في كل مرة.

جرب أولاً وحدة i.landsat.toar وقارن القيم مع نتائجك بالخطوة 1.

ثانيًا ، يجب أن تكون حذرًا مع وحدة i.pansharpen حيث أعتقد أنها تقبل بيانات 8 بت فقط ، لذا فإن البيانات العائمة الناتجة عن تصحيح الغلاف الجوي ستعطيك مشكلات.

إذا كنت ترغب في زيادة وضوح البيانات المصححة جوًا ، فحاول استخدام mapcalc (تأكد من أن إعدادات منطقتك ذات دقة 15 م × 15 م). تتم كتابة الصيغ في مستندات i.pansharpen ومتاحة هنا: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

لذا تحقق من نطاقاتك الفردية باستخدام منحدر لوني بتدرج الرمادي ، ثم إذا كانت جميعها تبدو جيدة ، فعليك الانتقال إلى المركب.


شاهد الفيديو: تصحيح مؤثرات الغلاف الجوى فى الصور الفضائية atmospheric correction