كيفية حساب المساحات المتحولة من صنف الى آخر

[شاكر احمد]

السلام عليكم ورحمة الله وبركاته،

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

لدي لو سمحتوا تساؤل هو كيفية حساب المساحات المتحولة من صنف الى صنف آخر من خلال صورتين فضائيتين مصنفتين الى 5 صنوف من الغطاء الارضي land cover وبإستخدام برنامج ArcGIS.

ارجو ان يكون جوابكم بإسلوب الخطوة خطوة ان امكن مع شكري وتقديري لكم وكل عام وانتم بخير.

[مهندس جيولوجي ابراهيم العمري]

اخي الكريم
السلام عليكم ورحمة الله وبركاته

يمكنك مراجعة الموضوع التالي ففية ان شاء الله الفائدة

مقارنة بين مساحات الاراضى الزراعية (http://www.gisclub.net/vb/showthread.php?t=6512&highlight=%CD%D3%C7%C8+%C7%E1%CA%DB%ED%D1)

بحوث في تطبيقات العلوم الجيومكانية في مراقبة الغطاء النباتي وتدهور الأراضي (http://www.gisclub.net/vb/showthread.php?p=34168)

رسالة ماجستير في كشف التغير change detection تحصلت علي درجة ممتاز (http://www.gisclub.net/vb/showthread.php?t=5796&highlight=%CD%D3%C7%C8+%C7%E1%CA%DB%ED%D1)

كما انصحك بالبحث في المنتدى فهناك الكثير من المواضيع المتعلقة
http://i176.photobucket.com/albums/w185/Abu_Omar/Thanks/EID.gif

[شاكر احمد]

اخي الكريم المهندس الجيولوجي ابراهيم العمري،

كل عام وانتم بخير وشكرا جزيلا لكم على الرد والاشارة الى المواضيع ذات العلاقة بالموضوع.

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

لذلك احببت ان اسأل هنا عسى ان تتكرم علينا بالإجابة والتوضيح بإسلوبك المميز، او ان بعضا من اعضاء المنتدى ممن لهم خبرة في برنامج ArcGIS يفيدنا ويفيد الاخرين و يشرح الطريقة بكيفية استعمال برنامج ArcGIS في حساب المساحات المتحولة من صنف الى صنف آخر من خلال صورتين فضائيتين مصنفتين الى عدد من صنوف استعمالات الاراضي او الغطاء الارضي.

شكرا جزيلا ثانية لك اخي الكريم ولكل من يساهم في الموضوع

[مهندس جيولوجي ابراهيم العمري]

اخي الكريم شاكر احمد

اشكرك أولاً على اسلوبك الجميل في الطرح، واشكرك ثانياً على اهتمامك وإطلاعك على المواضيع التي اشرت اليها في مشاركتي السابقة، فبارك الله فيك،،،

أما بخصوص أستفسارك فهو يحتاج إلى بعض التوضيح حتى اتمكن من الاجابة عنه ومن ذلك:
1- صورتين فضائيتين مصنفتين الى عدد من صنوف استعمالات الاراضي او الغطاء الارضي: هل تعني هنا أنه قد تم التصنيف مسبقاً فأي برنامج تم استخدامة لعمل التصنيف وما هو نوع التصنيف؟

2- في برنامج ArcGIS وفي اعتقادي يمكن الاستفادة من خاصية Change deduction للاجابه عن استفسارك حيث سوف يتم إنشاء طبقة جديدة تحتوي أو تمثل الاراضي المتحولة والتي يمكن بعد ذلك استخراج مساحاتها بكل سهولة، وهذا ملف يشرح طريقة الاستفادة من هذة الخاصية "بحثت عنه سابقاً ووجته الان فقط" وهو موضوع جميل جداً للاستاذ الفاضل رمضان الشافعي على هذا الرابط "كيفية عمل Change deduction بواسطة برنامج ARC GIS (http://www.gisclub.net/vb/showthread.php?t=881)".

3- أعتذر منك عن عدم اتمام الشرح وبصورة مفصلة كما تريد وذلك لعدم توفر البيانات لدي ونظراً لضيق الوقت هذة الايام فلو ارفقت بيانات مع اعطائي للوقت الكافي قد اعد موضوعاً يشرح ذلك.

هذا وتقبل فائق تقديري وكل عام والجميع بالف خير

[شاكر احمد]

اخي الكريم المهندس الجيولوجي ابراهيم العمري،

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

ارجو المعذرة ان لم يكن تساؤلي السابق واضحا بما فيه الكفاية لذا سأحاول هنا ان شاء الله ان اعيد كتابة المسالة بشكل آخر. اذا كانت لدي صورتين فضائيتين (لاندسات) ملتقطتين بتاريخين مختلفين ولتكن 1990 و 2007. بعد ان تم تصنيف هذه الصورتين ببرنامج ايرداس (تصنيف موجه و غير موجه) فكان الناتج لدي هو صورتين مصنفتين لعدد من صنوف استعمالات الأرض land use.

تساؤلي هنا هو كيفية حساب المساحات التي تحولت من صنف محدد في الصورة المصنفة لعام 1990 الى صنف آخر في الصورة المصنفة لعام 2007 ولكل صنف ببرنامج ArcGIS.

هل بالإمكان عمل جدول بتلك المساحات التي تحولت من صنف الى آخر، وهل هناك امكانية بالبرنامج لرسم خارطة بتلك التغيرات كالزيادة مثلا في صنف أو النقصان في صنف آخر. ولتكن صنوف استعمالات الأرض تلك على سبيل المثال هي: مناطق سكنية، مناطق صناعية، حقول زراعية، غابات، مسطحات مائية، أراضي جرداء، اراضي مراعي، اراضي محاصيل.

شكرا لكم وعذرا لإزعاجكم في ايام العيد، وارجو ان تأخذ وقتك الكافي للإجابة، وان احتجت حضرتك للصورتين المصنفتين فأني لااعرف كيف سابعث لك بهما اذ ليس لدي عنوانك الالكتروني.

جعل الله كل ايامكم وايام الجميع اعيادا. راجيا قبول فائق تحياتي وتقديري لك وللجميع.

[مهندس جيولوجي ابراهيم العمري]

يمكنك استخدام احد مواقع مشاركة الملفات والتي من اجملها
http://www.4shared.com
وكذلك موقع
http://www.box.net

[مهندس جيولوجي ابراهيم العمري]

أخي الكريم
السلام عليكم ورحمة الله وبركاته

هذه طريقة حساب المساحة لكل صنف Classes في الصورة المصنفة Classified raster، وهي الطريقة التي اعرفها مع تأكدي أن هناك طرق أخرى لإجراء نفس العملية، قمت بعمل شرح للطريقة على إحدى الصور ولم أكمل تاركاً الباقي لك وذلك للأسباب التالية:

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

اتبع الخطوات التالية:

أولاً قم بإدراج الصورة المراد العمل عليها ثم بعد تحديد نوع المرجع الجغرافي مثلاً ‘UTM1984ZoneN37" وإعداد مشروع العمل كذلك ليتوافق مع المرجع الصحيح للصورة.

ثانياً قم بتحويل الصورة من Raster to Vector كما في الصورة التالية:

http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/1.jpg
http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/2.jpg

ثالثاً قم بتوحيد التصنيفات Classes المتشابهة عن طريق الأداة Merge كما في الصورة التالية:

http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/3a.jpg
http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/3b.jpg
http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/4.jpg

وقم بتكرار هذه الخطوة على كل التصنيفات Classes المتشابهة في الخريطة أو Vector ولتحصل في النهاية على جدول مشابه تماما لما في الصورة الأولي الأصلية Raster.

رابعاً بعد الانتهاء من دمج وتوحيد التصنيفات Classes المتشابهة في الخطوة السابقة نقوم الآن بحساب المساحات لكل تصنيف Classes وذلك عن طريق الأداة Calculate Area في الأدوات كما في الصور التالية:

http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/5.jpg
http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/6.jpg

سوف نجد بعد ذلك انه تم إنشاء طبقة تحتوي على مساحة كل تصنيف وذلك حسب المرجع الجغرافي الذي عليك إعداده كما أشرنا سابقاً، وتظهر النتيجة كما في الصورة التالية:

http://i176.photobucket.com/albums/w185/Abu_Omar/Add%20Point-ArcGIS/AreaCulculate/7.jpg

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

مع تمنياتي لك بالتوفيق،،،

[محمد الحسن الطيب]

جزاك الله خيرا اخي ابراهيم
و مبروك عليك النسخة Arcgis 10

[شاكر احمد]

اخي الكريم المهندس الجيولوجي ابراهيم العمري،

أقدم لك شكري الجزيل لإهتمامك وحرصك على متابعة الموضوع والإجابة على تساؤلات الأعضاء وتقديم الحلول الممكنة لهم وجعل الله سبحانه وتعالى كل ذلك في ميزان حسناتك يوم لا ينفع لا مال و لا بنون.

لقد راجعت اجابتك المسهبة و المفصلة والواضحة اخي الفاضل كما في كل اجاباتك السابقة، وبدا لي من اجابتك ان النتيجة النهائية التي نتحصل عليها من تطبيق طريقتك هي حساب المساحة الصافية لكل صنف مثل الصنف AAA او BBB او CCC في الصورة الجديدة مقارنة بالصورة القديمة، اي بمعنى آخر كم زادت او نقصت المساحة الإجمالية للصنف الواحد بعد ان نطرح المساحة القديمة من المساحة الجديدة. لقد وجدت أن هذا مايمكن الحصول عليه من الطريقة التي تطرقت اليها حضرتك مشكورا.

ولكني كنت اتصور ان هناك طريقة تجد لنا كم تحول من الصنف AAA الى الصنف BBBو CCCو DDDو EEE وكم تحول من الصنف BBB الى AAAو CCCو DDDو EEE وكم تحول من الصنف CCC الى AAAو BBBو DDDو EEEو هكذا لبقية الصنوف، وبحيث تكون المخرجات بهيئة جداول وكذلك خرائط توضح التحول من صنف لآخر. وهذا في الحقيقة هو الذي اسأل وأبحث عنه إن كان ذلك ممكنا بإستخدام البرنامج العملاق ArcGIS.

قد يكون لديك أخي الكريم تعليقا على الموضوع و لا أعرف هل بإمكان ArcGIS ان يقوم بهكذا وظائف؟

أقدم خالص شكري لكم مرة ثانية، راجيا قبول فائق تحياتي وتقديري لكم و للجميع.

[مهندس جيولوجي ابراهيم العمري]

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

[محمد الحسن الطيب]

اخي ابراهيم هل من الممكن ان ترسل لي الشيب فايل لأكمل عليه طلب الاخ إن كان ذلك ممكن
و ربما اكون قد وجدت الحل
mohamedhacen@gmail.com
و الله اعلم

[محمد الحسن الطيب]

هل يمكن الحصول علي الصورتين و ذلك افضل للعملية
اخي شاكر

[مهندس جيولوجي ابراهيم العمري]

ممكن اخي الكريم هذا هو الشيب فايل مع الصور الاصلية على هذا الرابط:

http://www.4shared.com/file/bllnSYLL/CHANGE.html

أو

http://www.box.net/shared/6ad4vdx29f

وبارك الله فيك،،،

[محمد الحسن الطيب]

بارك الله فيك اخي ابراهيم
هل من الممكن الحصول علي المرئية الثانية و جزاك الله كل خير

[محمد الحسن الطيب]

لم انتبه ان الملف الثاني يحتوي علي المرئية الثانية
و جاري العمل و جزاك الله كل خير

[شاكر احمد]

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

ارجو قبول فائق احتراماتي و عظيم تقديري لكم وللجميع.

[مهندس جيولوجي ابراهيم العمري]

بالتوفيق اخي محمد الحسن ولد الشيخ محمد وأرجوا ان تفيدنا بالطريقة في حال وجدتها لتعم الفائدة
وبالتوفيق للجميع
قال تعالي: (وَتَعَاوَنُوا عَلَى الْبِرِّ وَالتَّقْوَى وَلا تَعَاوَنُوا عَلَى الإِثْمِ وَالْعُدْوَانِ وَاتَّقُوا اللَّهَ إِنَّ اللَّهَ شَدِيدُ الْعِقَابِ)

[محمد الحسن الطيب]

http://www.gisclub.net/vb/images/uploads/1362_7974c8eb98071d09.jpg

[محمد الحسن الطيب]

نسيت ان اذكر بأن الاماكن التي تتقاطع مع Bbb هي الزيادات التي تحولت من طبقات 1985 بستثناء Bbb85 فهي بالطبع موجودة ضمن Bbb95 أو كانت من ضمن التقاطع عند اماكن التغيير , و الله اعلم

[محمد الحسن الطيب]

http://www.gisclub.net/vb/images/uploads/1362_25084c8ebdbe92a6a.jpg
هذه الصورة توضح المناطق التي طرأت علي العنصر BBB من كل العناصر الاخري
العنصر BBB بالون الاصفر
العناصر الاخري بالون الازرق الفاتح

[مهندس جيولوجي ابراهيم العمري]

بارك الله فيك اخي محمد الحسن ولد الشيخ محمد
ومبروك عليك الاصدار العاشر ArcGIS10

[محمد الحسن الطيب]

http://www.gisclub.net/vb/images/uploads/1362_79864c8fb9d5bfcb0.jpg

[محمد الحسن الطيب]

http://www.gisclub.net/vb/images/uploads/1362_180044c8fc586829f2.jpg

[شاكر احمد]

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

كما لا انسى جهد المهندس ابراهيم العمري جزاه الله خيرا الذي ساهم في الجزء الاول من الحل. وفقكم الله و ادامكم لكل خير. والسلام عليكم ورحمة الله وبركاته.