المساعد الشخصي الرقمي

مشاهدة النسخة كاملة : التحويل من DN إلى Reflectance (1 من 2)



مهندس /وضاح عبدالله
02-03-2008, 10:44 PM
الدكتور / وسام
السلام عليكم ورحمة الله وبركاتة

أرجو تعاونك عن كيفية أستخدام Spatial Modeler في برنامج الإرداس لتحويل الصورة من قيم DN إلى Reflectance
وأيضاعن كيفية تطبيق معادلة الخط المستقيم في الموديلير. علما بأني طالب دراسات عليا وأحتاج تطبيق الموديلير
في بحثي.

جزاك الله خيرا

د. وسام محمد
02-04-2008, 08:50 AM
التحويل من القيمة الرقمية للبكسل DN إلى قيمة الأنعكاس الراديومتري
هذه الرسالة رداً على طلب الأخ المهندس وضاح عبد الله حول كيفية استخدام الـ Model Maker في برنامج ERDAS لتحويل قيمة البكسل الرقمية في الصورة DN إلى قيمة الأنعكاس الراديومتري Reflectance وذلك باستخدام طريقة الخط المستقيم.
تتم هذه العملية على مرحلتين:
المرحلة الأولى التحويل من DN إلى Radiance
وتستخدم فيها طريقة الخط التجريبي Empirical Line وتعتمد على المعادلة التالية:
Radiance = ((Lmax-Lmin)* Dn/Qcalmax)+Lmin
حيث Qcalmax هي أقصى قيمة يمكن للـ DN أن تحصل عليها (وتساوي 255 في حالة لاندسات مثلاً) ، Lmax و Lmin القيمة العظمى والصغرى للـ radiance في كل نطاق band ( وتكون في ملف الـ Metadata الخاص بالصورة)
الصورة المبنية فيما يلي تبين ملف metadata يحتوي هذه البيانات بالنسبة لصورة لاندسات TM:
(انظر الصورة رقم 1)
لتنفيذ هذا التحويل نقوم بالخطوات التالية:
1 - افتح برنامج ERDAS ومنه أنقر فوق Modeler ومنها اختار Model Maker لتظهر النافذة التالية:
(أنظر الصورة رقم - 2)
2 - اختار رمز الصورة وقم بسحبه ورميه فوق صفحة الmodeler كما في الشكل التالي:
(أنظر الصورة رقم 3)
3 - قم بالنقر مرتين فوق رمز الصورة لتنفتح نافذة خصائصه وقم بتخصيص اسم الصورة كما هو مبين في الشكل التالي:
(أنظر الصورة رقم 4)
4 - اختار رمز العملية وقم بسحبه والقاءه فوق صفحة النموذج ثم أختر السهم وقم برسمه من الصورة إلى العملية لتحصل على الشكل التالي:
(أنظر الصورة رقم 5)
5 - قم بالنقر مرتين فوق العملية لتنفتح نافذة العمليات ثم قم بإدخال المعالدلة كما هو مبين في الشكل التالي:
(أنظر الصورة رقم 6)
6- قم بأضافة صورة ليتم تفريغ الصورة الناتجة فيها واربطها بسهم مع العملية لتحصل على الشكل التالي:
(أنظر الصورة رقم 7)
7 - قم بالنقر فوق الصورة الناتجة لتقوم بتحديد خصائصها وتأكد أن نوع البيانات Float Single وانه تم استعاد الصفر من الحسابات الأحصائية، كما في الشكل التالي:
(أنظر الصورة رقم 8)
8 - قم بتكرار الخطوات من 4 إلى 7 لجميع النطاقات التي ترغب في تحويلها (أنا سأكتفي بالنطاقات الثلاثة الأولى). تحصل في النهاية على الشكل التالي:
(أنظر الصورة رقم 9)
9 - الخطوة التالية لتجميع الصورة الناتجة في صورة واحدة، أختار عملية وألقها ثم أرسم أسهم من الصور الناتجة إلى العملية ثم أنقر على العملية لتحريرها، اختار من قائمة Functions الـ Data Generation ثم أختار منها العملية STACKLAYERS واكتب أسامي النطاقات كما هو مبين بالشكل التالي
(أنظر الصورة رقم 10)
10 - في النهاية أضف صورة جديدة وسهم من عملية التجميع للصورة حتى يتم تجميع الصورة فيها من المفروض أن يكون الـ model الذي تستخدمه كمال بالشكل التالي:
(أنظر الصورة رقم 11)
11 - لا تنس أن تقوم بحفظ الـ Model ثم قم بتنفيذه عن طريق process>run
جرب هذه الخطوات وسوف أقوم بشرح المرحلة الثانية في رسالة أخرى لضيق الوقت، ولا تنسنا بالدعاء


http://www.gisclub.net/vb/images/uploads/773_3247047a6a61ebd0d2.JPG
الصورة رقم 1

http://www.gisclub.net/vb/images/uploads/773_2961247a6a6435c4c5.JPG
الصورة رقم 2

http://www.gisclub.net/vb/images/uploads/773_3262547a6a65a0292a.JPG
الصورة رقم 3

http://www.gisclub.net/vb/images/uploads/773_62947a6a6b38f3ac.JPG
الصورة رقم 4

http://www.gisclub.net/vb/images/uploads/773_2715447a6a6d088ef9.JPG
الصورة رقم 5

http://www.gisclub.net/vb/images/uploads/773_2960647a6a6f6b3441.JPG
الصورة رقم 6

http://www.gisclub.net/vb/images/uploads/773_2960747a6a70dba337.JPG
الصورة رقم 7

http://www.gisclub.net/vb/images/uploads/773_2961747a6a7207e732.JPG
الصورة رقم 8

http://www.gisclub.net/vb/images/uploads/773_75747a6a738b2bcc.JPG
الصورة رقم 9

http://www.gisclub.net/vb/images/uploads/773_75647a6a753a1ccf.JPG
الصورة رقم 10

حمود سعدد
02-04-2008, 09:00 PM
دكتور وسام وبالنسبة للبندات الحرارية

د. وسام محمد
02-05-2008, 06:50 AM
الأخ حمود كل النطاقات يمكن أن يتم عليهم نفس التحويل، لكنه عادة تستبعد النطاقات الحرارية من التحليل بسبب دقتها المساحية المختلفة عن بقية النطاقات وعدم وجود مكان لها في ترتيب النطاقات الخاصة بدراسة الظواهر الأرضية حيث أن النطاقات الحرارية تمثل الطاقة الكهرومغناطيسية المنبعثة من الأرض وليست المنعكسة عن سطح الأرض.

حمود سعدد
02-05-2008, 12:24 PM
جزاك الله الف خير

ندى صالح
02-05-2008, 03:58 PM
السلام عليكم ورحمة الله وبركاته
سلمت يمينك دكتور وســــــــــــــــــــــــــام على هذا الشرح الوافي
ولكني اسأل هل هذه الطريقه هي الطريقه المتبعه دائما في حساب الانعكاسيه الطيفيه من الصور الفضائيه بالنسبه الى القمر الصناعي لانسات
ولك مني جزيل والامتنان

مهندس /وضاح عبدالله
02-05-2008, 11:03 PM
السلام عليكم ورحمة الله وبركاتة
أشكرك يادكتور وسام وصدق من أسماك وسام النادي وجزاك الله خيرا وجعلة في ميزان حسناتك
شرح أكثر من ممتاز وقمت بتطبيقة والحمدلله وصلت لنتيجة وإن شاء الله سوف ننتظر باقي الشرح للتحويل إلى Reflectance

وكان لي طلب أخر أرجوك وهو أنة بعد الحصول على صورة الريفليكتنس إن شاء الله أنا عندي نقط مساحية ‏بالأحداثيات المترية وأريد أن أستخراج قيم ‏DN‏ لها
أن ستخراج قيم ‏DN‏ لها من صورة الريفليكتنس التي حصلنا عليها.‏

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

د. وسام محمد
02-06-2008, 07:15 AM
الأخ حمود جزانا الله وأياكم بالخير
الأخت ندى هذه الطريقة يمكن أتباعها مع كافة المجسات سواء المحملوة علة متن لاندسات أو غير من المنصات بشرط توفر القيم اللازمة للمعادلة وعادة ما تكون هذه القيم مرفقة بالصورة في صورة ملف يطلق عليه اسم الـ header أو الـ metadata.
الأخ وضاح جزانا الله واياكم بالخير لكن لك يكن هناك داعي لعمل اقتباس لكل الرسالة فقد صار رايك أكثر من واضح يا وضاح.
إن شاء الله أقوم بكتابة الجزء الثاني وبعد ذلك سنقوم معاً بتحديد قيمه الـ DN لنقاطك المحددة.

محمد ابريك كروش
03-07-2008, 01:43 PM
الف شكر للدكتور علي المعلومات القيمة

شاكر احمد
03-10-2008, 09:25 PM
الأستاذ الفاضل الدكتور وسام محمد،
السلام عليكم ورحمة الله وبركاته،
الشكر قليل بحقكم لما تقدموه من خدمات جليلة وإسهامات رائعة لجمهور النادي ولطلبة العلم.
إن سمحتم لدي سؤال عن موضوع ال Change Detection، فإن كانت لدينا صورتين فضائيتين للقمر الصناعي Landsat مأخوذة بالمتحسس ETM+ وملتقطة الأولى في الشهر الرابع والثانية في الشهر السابع من نفس السنة، لذا ينصح بإجراء عملية Radiometric normalization لواحدة منهما لمعالجة الاختلافات الزمنية بينهما والناتج عن اختلاف زاوية ارتفاع الشمس عن الأرض عند التقاط الصورتين وما يؤدي ذلك من اختلاف في الظلال على الأرض ودرجة سطوع الأغطية الأرضية.
رجائي لكم هو: هل تتفضلون في توضيح الخطوات العملية لمعالجة تلك الحالة وكيفية اجراء عملية Radiometric normalization بإستخدام برنامج ERMapper او Erdas Imagine لكي يمكن بعدها من انجاز المعالجة والتحليل لكلا الصورتين.
ولكم منا فائق الإحترام والتقدير.
أخوكم المهندس شاكر احمد

حمود سعدد
03-28-2008, 01:34 AM
السلام وعليكم
الف شكر وتقدير للدكتور وسام ولك الاقئمين على هذا الصرح العلمي وفقكم الله وسدد على طريق الخير خطاكم
دكتور وسام بالنسبة للتحويل من (digital number to surface reflectance) اذا نوفرت القياسات المعملية (spectral reflectance)
هذا ولكم جزيل الشكر والتقدير

د. وسام محمد
03-28-2008, 06:38 PM
أخي حمود سعد
أعتقد أن سؤالك يمكن ان أفهمه بطريقتين في نفس الوقت، فإما أنك تسأل عن طريقة استخدام الـ Model Maker في هذا العمل، أو أنك تسأل عن طريقة العمل بصورة عامة
بالنسبة للـ Model Maker أعتقد أن هذه المشاركة بجزيئها تبين طريقة استخدامه بصورة جيدة
أما عن طريقة ربط الإنعكاس الطيفي المعين معملياً بذلك المعين من المجس، فيمكنم أن تستخدم حساب الإرتباط والإنحدار البسيط لتبين ما إذا كانت هناك علاقة وشكل هذه العلاقة.

حمود سعدد
03-29-2008, 08:14 AM
الف شكر دكتور وسام على الايضاح ما اقصده كما قلت هو معايرة الانعكاس الذي من المجس مع الانعكاس الحقلي زذكرن طريقة (حساب الإرتباط والإنحدار ) دكتور لو تكرمت كيفية عمل هذه الطريقة
وفقك الله وسدد وسدد على طريق الخير خطاك

د. وسام محمد
03-30-2008, 07:38 AM
الأخ حمود
طريقة حساب الارتباط والإنحدار أو Correlation and Regression طريقة بسيطة ويمكن أن تعرف تفصيلتها من أي كتاب أحصاء وهي تعتمد على أنتاج معادلة خطية بسيطة كالتالية
y = mx + c
حيث y المتغير التابع و x المتغير المستقل
بينما
m = (Sigma XY - (Sigma X * Sigma y/n) ) / (Sigma x^2 - (Sigma x)^2/n)
و
c = average y - m * average x
وهذا بعد حساب معامل الارتباط طبعا
برجاء التأكد من المعادلات من أي مرجع إحصاء حيث أنني أكتب معتمداً على ذاكرتي