د. وسام محمد
02-07-2008, 08:23 AM
إذا كنت تقرأ هذا المقال فأنا أفترض أنك قد قرأت المقال السابق له في نفس الموضوع، لكن إذا لم تكن قرأت هذا المقال فأنصحك أن تقوم بقراءته على هذا الرابط.
http://www.gisclub.net/vb/showthread.php?t=2409
الجزء الثاني: التحويل من الـ Radiance إلى الـ Reflectance
لإجراء هذا التحويل نستخدم معادلة جديدة هي
Reflectance = (pi * radiance * d * d) / (ESUN * cos(theta))
حيث pi هي ثابت بقيمة 3.14 و d المسافة بين الشمس والأرض لحظة التصوير ، theta هي زاوية ميل الشعاع الشمسي الشاقط على الأرض، وأخيراً ESUN قيمة الأشعاع الشمسي المتوسط خارج الغلاف الجوي.
بالنسبة لـ d و theta يمكن الحصول عليهم من ملف الـ header الخاص بالصورة.
أما ESUN فهي المعضلة الكبرى حيث أنه يقوم على حسابها المحطة التي تستقبل الصور وتصدر في منشورات خاصة على سبيل المثال الكتيب المعنون بـ Modtran User Manual الصادر عن معمل أبحاث القوات الجوية الأمريكية وله أصداراتان أحدهما قديمة تعود إلى أواخر التسعينات والثانية حديثة صادرة عام 2006. وهذا الكتيب يضم البيانات الخاصة بالأقمار الصناعية الأمريكية Landsat 4,5 والفرنسية SPOT 1,2,3 والهندية IRS1,2,3
وعموما الجدول المبين في الصورة التالية يقدم قيم ESUN للمجس TM المحمول على متن القمر الصناعي Landsat.
http://www.gisclub.net/vb/images/uploads/773_1094347aa9454ca2f9.JPG
الآن لنبدأ العمل مع ERDAS حيث سوف نقوم بفتح نموذج جديد لنبدأ العمل عليه.
1 - أضف صورة إلى النموذج الجديد ثم قم بتخصيص هذه الصورة لتصبحة صورة الـ Radiance التي قمنا بإنتاجها في المقال السابق، عندئذ من المفروض أن يكون شكل نموذجك كما بالصورة التالية:
http://www.gisclub.net/vb/images/uploads/773_2022747aa948bac69f.JPG
2 - قم بأضافة عملية إلى النموذج (رمز الدائرة) ثم قم بربطها بواسطة سهم مع صورة الـ radiance ثم أنقر عليها مرتين لينفتح نافذة الـ Function Definition وأكتب فيها المعادلة المبينة في الشكل التالي:
http://www.gisclub.net/vb/images/uploads/773_1396447aa94ad36920.JPG
لا تنسى أني حصلت على قيم d و theta من ملف الـ heard أو الـ metadata الخاص بالصورة التي أستخدمها، أما قيمة ESUN فقد حصلت عليها من الجدول السابق. كما أنني أعمل على النطاق الأول فقط.
3 - أضف صورة جديدة وأجعلها تستقبل البيانات الناتجة عن العملية السابقة لا تنسى أن تجعل بيانات هذه الصورة من النوع Single Float.
4 - كرر الخطوة رقم 2 و 3 لكل النطاقات التي تريد أن تنفذ عليها عملية التحويل (لا تنسى انا سوف اكتفي فقط بالنطاقات الثلاثة الأولى) من المفترض أن تحصل على الشكل التالي في النموذج الخاص بك.
http://www.gisclub.net/vb/images/uploads/773_1481947aa94e7a86c0.JPG
5 - أضف عملية جديدة ثم ارسم أسهم من جميع النطاقات إلى هذه العملية، أنقر مرتين على العملية لتنفتح نافذة Function Definition أختار من العمليات Data Generation ومنها أختار Staklayers ثم أضف النطاقات. من المفترض أن يكون الشكل لديك كالتالي.
http://www.gisclub.net/vb/images/uploads/773_1481047aa950a1b18e.JPG
6 - أضف صورة ليتم أختزان نتاج العملية السابقة فيها ولا تنسى أن تكون من النوع Float Singe. من المفترض أن يكون نموذجك كما بالشكل التالي:
http://www.gisclub.net/vb/images/uploads/773_3126947aa953ab3c22.JPG
7 - قم بحفظ النموذج من خلال File>Save as ثم قم بتنفيذه من خلال Process>Run.
لا تحرمنا من دعائكم بالخير
http://www.gisclub.net/vb/showthread.php?t=2409
الجزء الثاني: التحويل من الـ Radiance إلى الـ Reflectance
لإجراء هذا التحويل نستخدم معادلة جديدة هي
Reflectance = (pi * radiance * d * d) / (ESUN * cos(theta))
حيث pi هي ثابت بقيمة 3.14 و d المسافة بين الشمس والأرض لحظة التصوير ، theta هي زاوية ميل الشعاع الشمسي الشاقط على الأرض، وأخيراً ESUN قيمة الأشعاع الشمسي المتوسط خارج الغلاف الجوي.
بالنسبة لـ d و theta يمكن الحصول عليهم من ملف الـ header الخاص بالصورة.
أما ESUN فهي المعضلة الكبرى حيث أنه يقوم على حسابها المحطة التي تستقبل الصور وتصدر في منشورات خاصة على سبيل المثال الكتيب المعنون بـ Modtran User Manual الصادر عن معمل أبحاث القوات الجوية الأمريكية وله أصداراتان أحدهما قديمة تعود إلى أواخر التسعينات والثانية حديثة صادرة عام 2006. وهذا الكتيب يضم البيانات الخاصة بالأقمار الصناعية الأمريكية Landsat 4,5 والفرنسية SPOT 1,2,3 والهندية IRS1,2,3
وعموما الجدول المبين في الصورة التالية يقدم قيم ESUN للمجس TM المحمول على متن القمر الصناعي Landsat.
http://www.gisclub.net/vb/images/uploads/773_1094347aa9454ca2f9.JPG
الآن لنبدأ العمل مع ERDAS حيث سوف نقوم بفتح نموذج جديد لنبدأ العمل عليه.
1 - أضف صورة إلى النموذج الجديد ثم قم بتخصيص هذه الصورة لتصبحة صورة الـ Radiance التي قمنا بإنتاجها في المقال السابق، عندئذ من المفروض أن يكون شكل نموذجك كما بالصورة التالية:
http://www.gisclub.net/vb/images/uploads/773_2022747aa948bac69f.JPG
2 - قم بأضافة عملية إلى النموذج (رمز الدائرة) ثم قم بربطها بواسطة سهم مع صورة الـ radiance ثم أنقر عليها مرتين لينفتح نافذة الـ Function Definition وأكتب فيها المعادلة المبينة في الشكل التالي:
http://www.gisclub.net/vb/images/uploads/773_1396447aa94ad36920.JPG
لا تنسى أني حصلت على قيم d و theta من ملف الـ heard أو الـ metadata الخاص بالصورة التي أستخدمها، أما قيمة ESUN فقد حصلت عليها من الجدول السابق. كما أنني أعمل على النطاق الأول فقط.
3 - أضف صورة جديدة وأجعلها تستقبل البيانات الناتجة عن العملية السابقة لا تنسى أن تجعل بيانات هذه الصورة من النوع Single Float.
4 - كرر الخطوة رقم 2 و 3 لكل النطاقات التي تريد أن تنفذ عليها عملية التحويل (لا تنسى انا سوف اكتفي فقط بالنطاقات الثلاثة الأولى) من المفترض أن تحصل على الشكل التالي في النموذج الخاص بك.
http://www.gisclub.net/vb/images/uploads/773_1481947aa94e7a86c0.JPG
5 - أضف عملية جديدة ثم ارسم أسهم من جميع النطاقات إلى هذه العملية، أنقر مرتين على العملية لتنفتح نافذة Function Definition أختار من العمليات Data Generation ومنها أختار Staklayers ثم أضف النطاقات. من المفترض أن يكون الشكل لديك كالتالي.
http://www.gisclub.net/vb/images/uploads/773_1481047aa950a1b18e.JPG
6 - أضف صورة ليتم أختزان نتاج العملية السابقة فيها ولا تنسى أن تكون من النوع Float Singe. من المفترض أن يكون نموذجك كما بالشكل التالي:
http://www.gisclub.net/vb/images/uploads/773_3126947aa953ab3c22.JPG
7 - قم بحفظ النموذج من خلال File>Save as ثم قم بتنفيذه من خلال Process>Run.
لا تحرمنا من دعائكم بالخير