سؤال حول التصحيح الراديومترى

[شيرين ايوب]

هل من المنطقى انه بعد عمل عمليه التصحيح الراديومترى و التحويل الى ريدينس تظهر قيم بيكسل اقل من صفر بالسالب و اعلى من 255 ؟؟؟؟؟
و اذا كان فما دلاله ذلك؟
و جزاكم الله خيرا

[د. وسام محمد]

هل من الممكن أن ترسلي المعادلة التي أستخدمتيها في عملية التصحيح، ربما أمكنني معرفة موطن الخلل.

[شيرين ايوب]

دكتور هى نفس المعادله التى استخدمتها حضرتك فى الشرح على المنتدى
على الرابط التالى
http://www.gisclub.net/vb/showthread.php?t=2409
ولا اعلم اين الخطأ أرجو إفادتى
و جزاك الله خيرا

[د. وسام محمد]

حسناً إذا كنت تستخدمين ERDAS Imagine 9.0 أو نسخة أحدث فهناك في القائمة Interpreter > Spectral Enhancement > Landsat 7 Reflectance هي أداة ممتازة لعمل التحويل مباشرة بمجرد إدخال القيم الموجودة في الـ Metadata الخاصة بالصور.
أما إذا كنت تريدين استخدام الـ Model Maker فيمكنك ارساله لي لأعرف الخطأ وأصححه.
تقبلي تحياتي.

[شيرين ايوب]

شكرا دكتور
و اقوم بتنفيذه المعادله على الشكل التالي
($n1_l5176038_03819980420_b10 * (193.000 +1.520) / 255) -1.520
و لكن هل يمكن التحويل الى ريفلكتنس قبل التحويل الى ريدينس و هل هذه الإضافه فى البرنامج لصور لاندسات 7 فقط ام هل يمكن استخدمها لكل الأقمار بما فيهم لاندسات1 ؟

[شيرين ايوب]

دكتور وسام
أين الخطأ فى المعادله
مع الشكر

[د. وسام محمد]

الأخت شيرين
هل أنت في مرحلة تحويل الصورة من Digital Number إلى Radiance لأني أحسب هذه المعادلة هي المستخدمة للتحويل من DN إلى Radiance أم أنا مخطاً. عموماً لو أنني غير مخطأ فمعادلتك تعني أن:
Radiance = DN*(Lmax - Lmin)/255+Lmin
حيث DN القيمة العددية للبكسل، وLmax وLmin هما قيمة الـ radiance العظمى والدينا التي يقيسها المجس Sensor في المنطقة الطيفية التي يمثلها هذا النطاق. مقارنة بمعادلتك هذا يعني التالي:
Lmax = 193.000
Lmin = -1.520
وفي هذه الحالة لابد أن تكون معادلتك سليمة يا شيرين وعليك التأكد من أن القيميتين اللتين استخدمتهما كما التعريف.

[شيرين ايوب]

طبقتها أكثر من مره يا دكتور و قد قمت بمراجعته البيانات بدقه شديده و لكن لا ادرى لماذا تظهر القيم بالسالب ؟
و اذا لم يكن الخطا فى تطبيق المعادله فهل يمكن ان يكون هناك خطأ فى المرئيه أم هل للسالب دلاله أم بذلك يعتبر التطبيق خطأ و سيأثر على عملية التصنيف فيما بعد .؟؟؟
شكرا دكتور على إهتمامك
وجزاك الله خيرا

[د. وسام محمد]

هل يمكنك أن ترسلي الـ Metadata أو الـ Header للصورة، ربما تكوني اخترت متغيرات خاطئة؟؟

[شيرين ايوب]

GROUP = L1_METADATA_FILE
GROUP = METADATA_FILE_INFO
ORIGIN = "Image courtesy of the U.S. Geological Survey"
REQUEST_ID = "0800903080099_00001"
PRODUCT_CREATION_TIME = 2009-03-09T05:01:19Z
STATION_ID = "EDC"
LANDSAT5_XBAND = "1"
GROUND_STATION = "XXX"
LPS_PROCESSOR_NUMBER = 0
DATEHOUR_CONTACT_PERIOD = "9021607"
SUBINTERVAL_NUMBER = "01"
END_GROUP = METADATA_FILE_INFO
GROUP = PRODUCT_METADATA
PRODUCT_TYPE = "L1T"
ELEVATION_SOURCE = "GLS2000"
PROCESSING_SOFTWARE = "LPGS_9.4.0"
EPHEMERIS_TYPE = "PREDICTIVE"
SPACECRAFT_ID = "Landsat5"
SENSOR_ID = "TM"
SENSOR_MODE = "SAM"
ACQUISITION_DATE = 1990-08-04
SCENE_CENTER_SCAN_TIME = 07:43:27.8250440Z
WRS_PATH = 176
STARTING_ROW = 38
ENDING_ROW = 38
BAND_COMBINATION = "1234567"
PRODUCT_UL_CORNER_LAT = 32.7002357
PRODUCT_UL_CORNER_LON = 31.0136193
PRODUCT_UR_CORNER_LAT = 32.7144422
PRODUCT_UR_CORNER_LON = 33.6188888
PRODUCT_LL_CORNER_LAT = 30.7580024
PRODUCT_LL_CORNER_LON = 31.0546440
PRODUCT_LR_CORNER_LAT = 30.7711770
PRODUCT_LR_CORNER_LON = 33.6060979
PRODUCT_UL_CORNER_MAPX = 313800.000
PRODUCT_UL_CORNER_MAPY = 3619800.000
PRODUCT_UR_CORNER_MAPX = 558000.000
PRODUCT_UR_CORNER_MAPY = 3619800.000
PRODUCT_LL_CORNER_MAPX = 313800.000
PRODUCT_LL_CORNER_MAPY = 3404400.000
PRODUCT_LR_CORNER_MAPX = 558000.000
PRODUCT_LR_CORNER_MAPY = 3404400.000
PRODUCT_SAMPLES_REF = 8141
PRODUCT_LINES_REF = 7181
PRODUCT_SAMPLES_THM = 4071
PRODUCT_LINES_THM = 3591
BAND1_FILE_NAME = "L5176038_03819900804_B10.TIF"
BAND2_FILE_NAME = "L5176038_03819900804_B20.TIF"
BAND3_FILE_NAME = "L5176038_03819900804_B30.TIF"
BAND4_FILE_NAME = "L5176038_03819900804_B40.TIF"
BAND5_FILE_NAME = "L5176038_03819900804_B50.TIF"
BAND6_FILE_NAME = "L5176038_03819900804_B60.TIF"
BAND7_FILE_NAME = "L5176038_03819900804_B70.TIF"
GCP_FILE_NAME = "L5176038_03819900804_GCP.txt"
METADATA_L1_FILE_NAME = "L5176038_03819900804_MTL.txt"
CPF_FILE_NAME = "L5CPF19900701_19900930_04"
END_GROUP = PRODUCT_METADATA
GROUP = MIN_MAX_RADIANCE
LMAX_BAND1 = 169.000
LMIN_BAND1 = -1.520
LMAX_BAND2 = 333.000
LMIN_BAND2 = -2.840
LMAX_BAND3 = 264.000
LMIN_BAND3 = -1.170
LMAX_BAND4 = 221.000
LMIN_BAND4 = -1.510
LMAX_BAND5 = 30.200
LMIN_BAND5 = -0.370
LMAX_BAND6 = 15.303
LMIN_BAND6 = 1.238
LMAX_BAND7 = 16.500
LMIN_BAND7 = -0.150
END_GROUP = MIN_MAX_RADIANCE
GROUP = MIN_MAX_PIXEL_VALUE
QCALMAX_BAND1 = 255.0
QCALMIN_BAND1 = 1.0
QCALMAX_BAND2 = 255.0
QCALMIN_BAND2 = 1.0
QCALMAX_BAND3 = 255.0
QCALMIN_BAND3 = 1.0
QCALMAX_BAND4 = 255.0
QCALMIN_BAND4 = 1.0
QCALMAX_BAND5 = 255.0
QCALMIN_BAND5 = 1.0
QCALMAX_BAND6 = 255.0
QCALMIN_BAND6 = 1.0
QCALMAX_BAND7 = 255.0
QCALMIN_BAND7 = 1.0
END_GROUP = MIN_MAX_PIXEL_VALUE
GROUP = PRODUCT_PARAMETERS
CORRECTION_METHOD_GAIN_BAND1 = "CPF"
CORRECTION_METHOD_GAIN_BAND2 = "CPF"
CORRECTION_METHOD_GAIN_BAND3 = "CPF"
CORRECTION_METHOD_GAIN_BAND4 = "CPF"
CORRECTION_METHOD_GAIN_BAND5 = "CPF"
CORRECTION_METHOD_GAIN_BAND6 = "IC"
CORRECTION_METHOD_GAIN_BAND7 = "CPF"
CORRECTION_METHOD_BIAS = "IC"
SUN_AZIMUTH = 107.6177898
SUN_ELEVATION = 56.5322914
OUTPUT_FORMAT = "GEOTIFF"
END_GROUP = PRODUCT_PARAMETERS
GROUP = CORRECTIONS_APPLIED
STRIPING_BAND1 = "NONE"
STRIPING_BAND2 = "NONE"
STRIPING_BAND3 = "NONE"
STRIPING_BAND4 = "NONE"
STRIPING_BAND5 = "NONE"
STRIPING_BAND6 = "NONE"
STRIPING_BAND7 = "NONE"
BANDING = "N"
COHERENT_NOISE = "N"
MEMORY_EFFECT = "Y"
SCAN_CORRELATED_SHIFT = "Y"
INOPERABLE_DETECTORS = "N"
DROPPED_LINES = "N"
END_GROUP = CORRECTIONS_APPLIED
GROUP = PROJECTION_PARAMETERS
REFERENCE_DATUM = "WGS84"
REFERENCE_ELLIPSOID = "WGS84"
GRID_CELL_SIZE_THM = 60.000
GRID_CELL_SIZE_REF = 30.000
ORIENTATION = "NUP"
RESAMPLING_OPTION = "CC"
MAP_PROJECTION = "UTM"
END_GROUP = PROJECTION_PARAMETERS
GROUP = UTM_PARAMETERS
ZONE_NUMBER = 36
END_GROUP = UTM_PARAMETERS
END_GROUP = L1_METADATA_FILE
END
دكتور ده ملف ميتا داتا خاص لصوره قمر لاندسات تى ام 1990 ارجو التوضيح وجزاك الله خيرا

[شيرين ايوب]

ممكن حضرتك تقولى الحل فى الحاله هذه ماذا اعمل ؟

[د. وسام محمد]

أنا لا أعرف أين يقع الخطأ، فاختيار Lmax و Lmin سليم، هل فعلتي أي شء بالصورة أي نوع من التحسين مثلاً؟؟

[شيرين ايوب]

لا ابدا يادكتور هذه اول خطوه لى بعد ان قمت بتنزيلها من النت
فما رأى حضرتك فى هذا الأمر ؟؟؟؟

[د. وسام محمد]

إعطيني عنوان الصورة لأجرب تصحيحه، فإن نجحت أرسلت إليك الخطوات تفصيلاً.

[شيرين ايوب]

دكتور الصور دى انا نزلتها من على موقع.gloves
L5176038_03819900804
هل هذا هو المقصود بعنوان الصوره?
مع العلم ان هذه المشكله لم تواجهنى فى هذه الصوره فقط و لكن فى باقى السلسله الزمنيه التى استخدمها لصور لاندسات تى ام او اى تى ام

[شيرين ايوب]

دكتور، هل يوجد حل اخر؟

ــــــــــــــــــــــــــــــــــــــــــــ

ادارة النادي : تم تعديل المشاركة

[د. وسام محمد]

سوف اعمل عليها يوم الاجازة وانتظري الرد مساء الجمعة او السبت

[شيرين ايوب]

شكرا يا دكتور و جزاك الله خيرا

[د. وسام محمد]

انا احاول تنزيل الصورة من الانترنت، حجمها 190 ميجا بايت والخط عندي ليس سليم، سأحاول ثانية

[محمد جميل]

الاخت شرين

ممكن ضغط الصورة وتقسيمها الي عده ملفات لتسهيل عملية ارسالها ع الميل او حتي رفعها الي اي موقع ليتمكن د.وسام من تنزيلها
وافادتنا بالامر

مع الشكر

[شيرين ايوب]

الأخ محمد
حاولت و لكن مع الأسف النت عندى ايضا به مشكله فلا يمكننى سوي فتح المتصفح فقط .
و انا فى أمس الحاجه لحل هذه المشكله

[المهندس الصغير]

السلام عليكم
بخصوص معادلة التصحيح التي كتبتها الاخت شيرين، أود ان اقول انها معادلة صحيحة وان قيم Dnللمرئيات التي تتعامل معها ايضاً صحيحة ويمكن التأكد من ذلك ببساطة وذلك بتعويض Dn=0 أو اي قيمة أقل من 1.52 فإن ناتج المعادلة يكون سالباً. بإمكانكي ان تدققي ذلك.

[شيرين ايوب]

و هل هذا صحيح أن تكون هناك قيم ناتجه بالسالب ؟
أرجو توضيح ما أشرت أليه إذا سمحت
وجزاك الله خيراً