المعايرة الراديومترية والتصحيحات الجوية لصور الأقمار الصناعية

في     هذا الدرس  سأقوم بتقديم  موجز نظري عن عن   التصحيحات الجوية والمعايرة   لصور بعض

الأقمار الصناعية وساقوم    بعمل المعايرة بشكل يدوي باستخدام حاسبة   الصور في أي برنامج ( ENVI,Qgis,arcmap)

لفهم    الطريقة التي   تستخدمها البرمجيات لمعايرة الصور .

ومن ثم حساب المؤشر  NDVI

(Band NIR-  BAND red)/ (Band NIR+  BAND red)

 للبيانات  قبل المعايرة والتصحيحات الجوية  DN

للبيانات   بعد حساب الانعكاس أعلى الغلاف الجوي  TOA

للبيانات   التي تمثل الانعكاس  على سطح الأرض BOA

ومقارنة النتائج مع بعضها   البعض لإظهار دور التصحيحات والمعايرة في  الحصول على نتائج دقيقة

توضيح  بعض المفاهيم

المعايرة الراديومترية    والتصحيحات الجوية : عمليات   تفيد في تحسن جودة البيانات و تحسن امكانية تفسير الصور  من خلال معايرة قيم البيكسل للحصول قيم انعكاس تمثل القيمة الحقيقة  للبيكسل علي سطح الأرض حيث تكون القيمة المسجلة في حساسات القمر الصناعي (DN) والطائرات   مختلفة عن القيمة الحقيقة المنبعثة من سطح الأرض (BOA)

Digital Number -DN

وهي  البيانات الخام  التي تمثل الطاقة التي  يتم ملاحظتها وقياسها على الحساس مباشرة وعلى الرغم  من أن DN ترتبط بقيم الانعكاس السطحي إلا أنها ليست متماثلة ويجب    معايرة هذه البيانات للحصول على قيمradiance (الإشعاعية ) او الانعكاس على قمة الغلاف الجوي   reflectance top of Atmosphere ومن ثم إجراء التصحيحات الجوية للحصول على قيم الانعكاس الحقيقة علي سطح الارض  TOP من البيانات الخام المكتسبة مباشرة من قبل الحساس حيث لا نستطيع استخدامها مباشرة لحساب المؤشرات وإجراء التحليلات والمقارنة حيث سنحصل على قيم  خاطئة وكذلك الصورة بدون معايرة تكون باهتة وصعبة التفسير,وتتم المعايرة للحصول على radiance أوreflectance of top Atmosphere           بتطبيق بعض المعادلات باستخدام البيانات Metadata لبيانات القمر الصناعي الموجودة ضمن الملف MTL file

مثال   عن هذه البيانات   بيانات القمر لاندسات  landsat -level-1 من المرحلة الاولي   

TOA Radiance and TOA reflectance

حيث   يتم التحويل من نمط البيانات الخام  DN إلى الانعكاس اعلا الغلاف الجوي TOA  او الاشعاعية أعلى الغلاف الجوي radiance  حسب ماتقضي الحاجة وسيتم التحويل حسب المعادلات  التي سيتم شرحها أدناه وهي الطريقة التي تتبعها البرامج   في معايرة البيانات

Conversion   from DN to TOA Radiance

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

واحدة الإشعاعية في النظام الدولي للوحدات هي واط لكل متر مربع لكل ستراديان (W·sr−1·m−2

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

  نقوم بتطبيق المعادلة التالية

Lλ = MLQcal + AL

Lλ     

    قيم الاشعاعية  radince

 مدخلات   المعالدلة  

AL   

       (RADIANCE_ADD_BAND_x (رقم  الباند )

ML       

  (RADIANCE_MULT_BAND__(رقم  الباند )

Qcal     

  pixel values -DN  البيانات الخام المسجلة على الحساس

القيم السابقة   التي تستخدم لحساب   الإشعاعية أعلى الغلاف الجوي  نحصل عليها من الملف MTL المرفق مع بيانات القمر الصناعي اثناء التحميل

Conversion from DN  to TOA Reflectance

للحصول  على قيمة الانعكاس في  أعلى الغلاف الجوي يجب  تطبيق المعادلة التالية  ولدينا الصيغتين التاليتين حيث  يتم استخدام احيانا زاوية ارتفاع  الشمس  أو الزاوية   المحصورة بين  الشمال و الشعاع الواصل بين الشمس  و الموقع

قيم المعادلة    يتم الحصول عليها  من MTL file المرفق   البيانات اثناء التحميل

Lλ = (MLQcal + AL)/ sin (θSE ) or Lλ = (MLQcal + AL)/ cos (θSZ )

 مدخلات   المعالدلة  

AL   

 (RADIANCE_ADD_BAND_x (رقم  الباند )

ML       

(RADIANCE_MULT_BAND__(رقم  الباند )

Qcal     

  pixel values (DN)  البيانات الخام المسجلة على الحساس

θSE

زاوية ارتفاع الشمس

θSz

الزاوية    المحصورة بين اتجاه الشمال    و الشعاع الواصل بين الشمس والموقع

قيم المعادلة    يتم الحصول عليها  من MTL file المرفق  مع البيانات اثناء التحميل

surface reflectance BOA

 

 Conversion from TOA Reflectance to surface reflectance

بعد حساب  قيم الانعكاس  اعلا الغلاف الجوي  TOA  تنتهي مرحلة المعايرة ,وتبدأ مرحلة التصحيحات الجوية   للحصول على قيم الانعكاس على سطح الأرض BOA .   وهذا يتطلب معرفة   الظروف الجوية للموقع .  وزمن التقاط الصورة وغيرها   من المعلومات

Atmospheric Correction Models

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

Dark Object Subtraction Method

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

بعض  البرمجيات التي   تساعد في إجراء التصحيحات الجوية للبيانات

برنامج السناب من خلال الإضافة sen2core    يقوم بتصحيح بيانات القمر سينتينيل 2 من المستوى الأول level 1-C وفي الرابط التالي  شرح لطريقة التنصيب وطريقة إجراء التصحيحات الجوية علي بيانات القمر سينتينيل 2  

او من خلال الخوارزمية Quick atmospheric  correction في برنامج الاينفي

من  خلال الإضافة   sime classification   لبرنامج ال Qgis بتطبيق الخوارزميةDOS1

حيث سيتم شرحها في المثال  العملي التالي

 

فكرة   المشروع

سيتم حساب   قيم المؤشر NDVI    على بيانات لم تخضع لعملية المعايرة والتصحيحات    ومن ثم علي بيانات تم معايرتها وتحويلها إلى انعكاس أعلى الغلاف الجوي  TOA ومن ثم حساب قيمة المؤشر بعد معايرة البيانات وتصحيحها BOA
لمعرفة   كيف تؤثر  المعايرة والتصحيحات على  النتائج

 في ال 31من شهر   مايو 2003 تعرض مصحح    خط المسح Scan Line Corrector -SLC

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

اولا وسنقوم بمعالجة الصورة بإزالة هذه الفجوات باستخدام     برنامج الأنفي
ثانيا    معايرة هذه البيانات   تحويلها الي إشعاعية radaince  وانعكاس TOA أعلى الغلاف الجوي  بشكل يدوي باستخدام المعادلات المشروحة   سابقا

ثالثا وحساب الانعكاس على سطح الارض BOA   مباشرة من قيم الحساس DN

باستخدام  الخوارزمية DOS1       المستخدمة في برنامج Qgis عن طريق   الايضافة

semi  classification

تحميل البيانات وتصحيح الفجوات

نقوم   بتحميل البيانات   بيانات القمر الصناعي landsat 7

LE07_L1TP_173037_20040914_20170119_01_T

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

• 

وسيتم    تصحيحها بنسخ  الاضافة landsat gap fill التي قد  لا تكون ضمن صندوق ادوات البرنامج  

لذلك يجب إضافتها   انظر هنا

نحدد   الباند التي يتم تصحيحها     وفي مثالنا سنقوم بتصحيح الباند رقم 4   والباند رقم 3

 ملاحظة :  تستطيع تصحيح   كل الصورة دون الحاجة إلى  تصحيح كل باند علي حدا

• 

بعد   اجراء التصحيح     نلاحظ اختفاء الفجوات

• 

ونعيد       الخطوة   السابقة للحزمة    الثالثة    ايضا

ثم  نقوم بحساب   قيمة المؤشر NDVI   من الحزمتين
ملاحظة باستخدام   البيانات DN ستكون   النتائج غير صحيحة للمؤشرات  

• 

حساب الانعكاس أعلى الغلاف الجوي TOA من نمط بيانات الحساسDN

نقوم   بالتحويل إلى  انعكاس    أعلى الغلاف الجوي   باستخدام برنامج ال Qgis

باستخدام   المعادلة المشروحة سابقا     باستخدام حاسبة الصورة

   بتطبيق المعادلة التالية :

L_λ = (M_LQ_cal + A_{L)/ sin} (θ_SE )

نقوم بفتح الملف  MTL    الذي يحتوي علي الميتاداتا

• 

 حيث    نلاحظ ان الملف  يحتوي على زاوية ارتفاع الشمس والمسافة بين الأرض  والشمس وغيرها من المعلومات كما توضح الصور التالية

• 

وكذلك القيم الاخرى التي تعوض في المعادلة السابقة

• 

ونقوم     بتعويض هذه   القيم في  المعادلة السابقة  مع الانتباه الي ان زاوية ارتفاع الشمس تحسب  بالراديان

• 

نعيد الخطوات نفسها للباند الاخرى    

• 

بعد الحصول على  جميع الباندات بصيغة انعكاس أعلى الغلاف الجوي   نقوم بحساب قيمة المؤشرNDVI

• 

ملاحظة

ونستطيع   الحصول على  الانعكاس والاشعاعية  أعلى الغلاف مباشرة دون   الحاجة لتطبيق المعادلات السابقة
 باستخدام بعض البرمجيات كالاينفي

نلاحظ   الفرق بين  قيم NDVI  المحسوب من البيانات  الخام ND وقيمة NDVI     المحسوبة من تحويل البيانات إلى   انعكاس أعلى الغلاف الجوي TOA

نلاحظ   ان الفرق كبير   حيث يصل إلى الضعف  أحيانا كما تظهر الصورة التالية  حيث يظهر المؤشر بدون معايرة البيانات     قيمة منخفضة 0.2 تدل على غطاء نباتي فقير  بينما تعطي قيمة الموشر بعد معايرة البيانات   القيمة 0.41 التي تدل على وجود غطاء نباتي

• 

تصحيح البيانات باستخدام الخوارزمية DOS1 للحصول علي الانعكاس علي سطح الارض BOA

الان سأستعرض طريقة   ضمن الاضافة semi classification   تقوم بتصحيح البيانات مباشرة من DN الحصول  على انعكاس  على سطح الأرض:

حيث   يتم حساب     الانعكاس على سطح الارض   بالمعادلة التالية

ρ=[π∗(Lλ−Lp)∗d2]/(ESUNλ∗cosθs)

حيث   يدخل في هذه المعادلة المسافة  بين الارض والشمس ومسار الاشعة عبر الغلاف الجوي  والقيمة

ESUNλ   التي تعطي بجداول  خاصة وهناك   عدة خوارزميات   من   DOS      والطريقة المتبعة بالبرنامج هية  Dos1   التي    لها شروط بارومترات خاصة نعوضها  بالمعادلة العامة للحصول  في المعادلة أعلاه.

  لمزيد من المعلومات راجع الرابط   التالي

نقوم    بفتح برنامج ال Qgis    ومن الاضافة   semi classification   نقوم    باضافة الحزم  وكذلك اضافة الملف MTL    كما تظهر الصورة  التالية

• 

    ونلاحظ    أيضا بأن القيم السابقة يتم   استيرادها من الملف MTL   وتملأ  الجداول   بشكل اوتماتيكي

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

• 

 بعد   انتهاء التعديل  نلاحظ ان البرنامج    يقوم بتغيير   اسم    البيانات   ساستعرض    مثلا   القناة   السادسة 
حيث  يتم   ايضافة   الاحقة RT_

RT_LE07_L1TP_173037_20040914_20170119_01_T1_B6_VCID_1

   ومن البيانات   المعدلة     ساقوم  بحساب   الموشر    NDVI   من البيانات   التي   تكون  عبارة عن انعكاس  حقيقي   من علي  سطح الارض  

وبحسب الموشر NDVI من البيانات المعدلة

• 

بمقارنة  النتائج نلاحظ  ان قيمة المؤشر   المحسوبة من البيانات الماخوذة من  الحساس    مباشر digital  numbers تعطي نتائج     خاطئة بينما استخدام بيانات الانعكاس أعلى الغلاف الجوي TOA وبيانات   الانعكاس على سطح الأرض BOA تعطي نتائج لحد ما متقاربة .لذلك عند القيام بأي  مشروع استشعار يجب فحص البيانات المستخدمة ومعرفة اي مستوى تكون البيانات المتاحة    وهل تحتاج الى  عملية تعديل ومعياره   وبشكل خاص عند استخدام سلاسل  زمنية من الصور

• 

ملاحظة

تستطيع   تحميل صور القمر  سينتينيل 2 من المستوى الثاني  مباشرة level-2/A ولكن هذه الصور
متوفرة   بداء من 2 أيار 2017 وتكون من النمط BOA

كاتب المقال:Martin Ealya

نشر البيانات المكانية باستخدام الجيوسيرفير

  المهارات المكتسبة من الدرس

طريقة تنصيب الجيوسيرفير
طريقة نشر البيانات المكانية عن طريق الجيوسيرفير
 اجراء بعض الاستعلامات  علي البيانات  في متصفح الويب باستخدام  الفلتر CQL
 تصدير   البيانات الي نظام المعلومات الجغرافية  ك WFS

قبل البدء بنشر البيانات المكانية عن طريق الجيوسيرفير ساقوم بتوضيح بعض المصطلحات الهامة حول الجيوسيرفير والنشر عبر الويب

WMS

عبارة عن  استجابة لطلب خريطة من  قاعدة بيانات جيومكانية  او اكثر .  والطلب يتضمن  الامتداد الجغرافي   و والطبقات المطلوب  تحميلها  والاستجابة عبارة عن ملف    HTML   يحمل ملف راستر   صورة يتم عرضها بالمتصفح او    اي نظام معلومات جغرافية

WFS  

 وهي    نفس  مبدء الخدمة السابقة ولكن الاستجابة تكون  عبارة عن طبقة شعاعية  يمكن   تعديلها   واجراء التحليلات عليها  من قبل العميل باستخدام نظام العلومات الجغرافية

Geography Markup Language –GML

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

OGC

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

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

الكوارث …..الخ

واصبحت   منذ عام     2007   عضو  في   موسسة Wide Web Consortiums -W3C
رابط الصفحة

  http://www.opengeospatial.org

W3C

هو   عبارة عن مؤسسة  دولية تعمل علي  تطوير الويب    وعمل معايير  موحدة   تأسست في 1 أكتوبر 1994 في مختبر MIT لعلوم الكمبيوتر في كامبريدج ، ماساتشوستس.
رابط     الصفحة
https://www.w3.org/participate/

OpenLayers

هي   عبارة عن مكتبة    تمكن من    وضع    خرائط  ديناميكية  في صفحات الانترنيت  بالاضافة   الى تحميل
  الطبقات    والمتجهات  علي هذه  الخرائط    وهي مفتوحة المصدر  تماما  وتسخدم في لغة البرمجة جافا   سكريبت
رابط الصفحة
https://openlayers.org/

مقالة تعليمية عنها من موقعنا

https://thatgisguy.com/2017/10/07/web-map/

GoServer

 هو خادم برمجيات مفتوح  المصدر   مكتوب بلغة جافا   تتيح للمستخدم  مشاركة  وتحرير  البيانات الجيومكانية ونستطيع   نشر مخطط او خريطة  باستخدام    معايير     وضوابط       OGC

 وباعتباره  مشروع  مفتوح المصدر   تم تطويره   واختباره  ودعمه   من  قبل مجموعة   متنوعة من الافراد   والمؤسسات  في  جميع  انحاء العالم  والبرنامج   هو المرجع   لخدمات   OGCالتي تقوم   بنشر  خدمة WMS   ,WFS عن    طريق برنامج الجيوسيرفير   .

وقد بداء البرنامج في بادئ الامر   في عام  2001     من قبل    مشروع التخطيط المفتوح  TOPP  والبرنامج عبارة عن بعض الادوات  تساعد المواطنين في مشاركة الحكومة  في عملية التخطيط  للمدن وتطوير  المناطق الحضرية  وجعل الحكومات  اكثر  شفافية  في التخطيط  من خلال  الامكانية  علي نشر المخطاطات   واجراء  الاستعلامات عليها

وكانت   هذه الفكرة   هي النواة  لتطوير الجيوسيرفير  وجعله  مشابه  لشبكة الانترنت   في تبادل وتحميل  والبحث عن المعلومات الجغرافية.

  وقد قام المطورون    للجيوسيرفير    بتأسيس مشروع  Geotools    عبارة   عن ادوات جافا GIS-Java

 قدرة الجيوسيفير    على نشر خدمتي  WMS  ,WFS   بالاضافة الي  ربطه ايضا بمشاريع   اخرى   مثل Postgis  ,openlayer   وغيرها عزز من فعالية   الجيوسيرفير  

حيث يمكن    للجيوسيرفير  قراءة  البيانات  من اكثر 12  مصدر  للبيانات  المكانية  ونشرها بصيغ   مختلفة .

تحميل   الجيوسيرفير

 ملاحظة   هامة اثناء  خطوات    التنصيب سيتطلب منك ادخال   اسم المستخدم وكلمة المرور   قم      بحفظ كلاهما بكتابتهما على ورقة خارجية  لانك لن تستطيع  الدخول   الى  الجيوسيرفر   مجددا في حال نسيتهما    

نذهب   الي  الصفحة الرسمية لتحميل الجيوسيرفير

قم بتحميل ملف التنصيب حسب نظام التشغيل الخاص بك

 نقوم  باختيار   اسم المستخدم وكلمة المرور 
ومن ثم ادخال : 8080 

بعد اكتمال تنصيب البرنامج يجب   البدء باقلاع السيرفير حيث  عند تشغيل الجهاز يجب القيام دائما بهذه العملية

بعد  الانتهاء من الاقلاع   نقوم بفتح 
 geoserver web page Admin
حيث  تكون ضمن المتصفح الافتراضي  علي الجهاز خاصتك

   و نقوم بادخال كلمة السر  التي قمنا بانشاءها   خلال تنصيب البرنامج

البداء    بعملية    نشر البيانات  الجفرافية

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

يمكنك  تحميل البيانات  من هنا

1-workspace

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

نقوم بعمل ورك سبيس تحمل الاسم فرانكفورت

نلاحظ   ظهورها في القائمة المنسدلة

 في   حال فتحنا المسار التالي
C:\Program Files (x86)\GeoServer 2.14.1\data_dir\workspaces\FRANKFURT
 نلاحظ تشكل مجلد يحمل اسم الورك سبيس الذي قمنا بانشاءه   كما توضح الصورة التالية حيث تستطيع لاحقا الوصل الى  كل الملفات التي ستضمنها   والتي سنقوم باضافتها بالخطوات  التالية

2-stores

بعدها  نذهب الي القائمة  Store   حيث   سنحدد    نمط البيانات  التي  سنقوم باضافتها    ونشرها  لاحقا . في   مثالنا  فقط shapefil    حيث   سنقوم    بنشر  اربع طبقات
طبقة  تمثل الابنية السكنية  عبارة عن مضلعات
طبقة تمثل    الابنية الحكومية عبارة عن مضلعات
طبقة   تمثل  الشوارع   عبارة  عن نقط
طبقة تمثل الشوارع  عبارة عن خطوط 

 نقوم باختيار اضافة طبقة كما في الصورة ادناه

نقوم   باضافة   الطبقة   من جهاز  الكمبيوتر   الى    الجيوسيرفير  من خلال الضغط على المستعرض وادخال الطبقة

3-layers

نذهب الي التبويب طبقة لاظهار الطبقة التي قمنا بانشائها في الخطوة السابقة

بعد انشاء الطبقة  نقوم بفتح  خصائص الطبقة  لضبط  حدود الطبقة   ونظام الاسقاط  و والعرض  وغيرها من الصفات

نلاحظ  في الصورة ادناه ان نظام الاسقاط المستخدم  يعبر عنه برقم   25832
EPSG :25832
بعدها نقوم بحساب امتداد الطبقة  ويتم حسابها بشكل  اوتوماتيكي  بالضغط  على  الخيار الذي تظهره الصورة التالية

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

بعد    ضبط  الامتداد   الطبقة    تصبح جاهزة للنشر  ولكن بدون  تنسيق   حيث   تظهر بالتنسيق الافتراضي  وسنقوم   بالخطوات اللاحقة للتحكم بالعرض  

4-layer preview

سنقوم بعرض   الطبقة     بالاعداد ت  الافتراضية   من  خلال الذهاب الي التبويب  عرض  الطبقة

ساقوم بعرضها  علي شكل open layer    حيث  تستطيع عرضها باستخدام انماط اخرى kml,GML   و…….غيرها ايضا

نلاحظ  ظهور طبقة   الشوارع   في المتصفح   مباشرة  مع امكانية التكبير والتصغير والتحكم بالعرض   وتطبيق  فلتر لعرض النتائج التي يرغب بها المستخدم   

CQL filter

يستطيع المستخدم اجراء العديد من الاستعلامات من خلال هذا الفلتر حيث يمكن من عرض الطرقات التي طولها اكبر من 400 متر في هذا القسم من فرانكفورت

• 

5-styles

 يمكن  الجيوسيرفير من تصدير لبيانات  بصيغ مختلفة   

 من التبويب  نمط  نستطيع التحكم    بنمط العرض   حيث تظهر القائمة الانماط المتاحة.   وهي الخيارات الافتراضية الموجودة في الجيوسيرفير  وسنقوم   بخطوات لاحقة   بانتاج النمط الخاص بنا

عند اختيار اي نمط  من الانماط السابقة   نلاحظ انه     مكتوب علي شكل  كود xml

سنقوم   بتصدير النمط  الخاص بنا  الي  الجيوسيرفير       حيث  سنقوم   بتصدير   الطبقة  الي قنا بنشرها Qgis   علي   شكل   WFS  ومن  ثم  تصنيف     هذه  الطبقة   فنلاحظ   هنا   فعالية الجيوسيرفير  في
نشر البيانات الجغرافية وامكانية اجراء   التحليلات   علي البيانات   ضمن نظم المعلومات الجغرافية

نقوم بنسخ  الجزء  من  الرابط   الذي تظهره الصورة  ادناه

Add WFS

    نضغط   علي Add WFS layer   من   القائمة    اضافة   طبقة   او من خلال  الرمز الذي يظهر  في القائمة  علي يسار  الشاشة

فتظهر   نافذة     نضغط     اولا  علي الخيار edit  من   ثم    تظهر نافذة جديدة
نقوم بنسخ الرابط الذي قمنا بنسخه    في الحقل URL

نضغط      بعدها    علي     connect    فنلاحظ   يتم تحميل الطبقة  في   النافذة   التي بالاسفل      نقوم   بتاشيرها  ونضغط   على  اضافة

نلاحظ يتم  اضافة طبقة  الشوارع  كخدمة  WFS     التي  تملك  نفس  خصائص shapefile   ونقوم   بالدخول الي خصائص     الطبقة   من   التبويب  Style و نقوم بعمل تصنيف  كما تظهر الصور ادناه مع  الخطوات    ومن

ومن ثم نقوم بتصدير  التنسيق بالصيغة SLD

SLD: هو  كود  XML محدد بواسطة OGC لتحديد مظهر طبقات الخريطة. ،  يمكن وصف مظهر الطبقة بغض  النظر  عن العناصر المرسومة ضمن  الطبقة   ولكن   لنوع محدد  (    نقطي او شعاعي او مضلع ).

نقوم   بحفظ   التنسيق   وفق   الرابط التالي ضمن    ملفات الجيوسيرفير

C:\Program Files (x86)\GeoServer 2.14.1\data_dir\styles

نقوم باضافة التنسيق الجديد الي التنسيقات الافتراضية المخزنة بالجيوسيرفير

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

 بعد الضغط  علي تحديث  يتم تحويل النمط SLD  الي كود XML  يظهر في النافذة التي بالاسفل  وبعدها نضغط علي  حفظ فيصبح من الانماط   المعتمدة ضمن الجيوسيرفير  

• 

الان نعيد  نشر الطبقة السابقة  مجددا ولكن   باضافة العرض الي قمنا بانتاجه  في الخطوات السابقة   بالدخول    الي
 layer التبويب

• 

 ندخل علي القائمة نشر كما تظهر الصورة ادناه

• 

• 

نقوم بعرض الطبقة من جديد بعد اضافة التنسيق الجديد لطبقة الشوارع

• 

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

• 

نعيد الخطوات السابقة للطبقات الثلاث الباقية ونقوم بنشرها

طبقة المباني الحكومية

• 

طبقة الابنية السكنية

• 

طبقة المدارس

• 

5-layers Groups

  في هذه المرحلة نقوم   با ضافة جميع    الطبقات   في طبقة   واحدة  ومن   ثم اعاد    العرض  لهذه المجموعة معا

نذهب     الي التبويب add groups layer      ونقوم باضافة   مجموعة طبقات جديدة

• 

ندخل اسم المجموعة التي سنقوم بانشاءها

• 

ونقوم ايضا باضافة   نظام الاسقاط للمجموعة     من ثم حساب الامتداد 

 

• 

نقوم   باضافة  الطبقات السابقة  الي المجموعة  كما   تظهر الصورة  ادناه

• 

    بعد ان يتم اضافة   كل الطبقات   نقوم  بالتحكم  بترتيب الطبقات  حيث    تكون   الطبقة  في الاعلى    بترتيب القائمة    هي الطبقة الاساس بالعرض   

• 

بعد الانتهاء من ترتيب الطبقات نضغط على حفظ ونقوم بعرض المجموعة معا بالذهاب الي التبويب عرض

• 

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

Martin Ealya : كاتب المقال

مراقبة تغيرات منسوب مياه نهر النيل عبر تحليل صور القمر سينتينيل 1

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

حتى تستطيع ان تقارن النتائج وتكتشف فعالية كل قمر

المهارات المكتسبة من الدرس

• • التعرف على برنامج سناب

 

• • تعلم على افضل طريقة لقص الصورة أثناء مشاريع المراقبة للحصول على نفس المشهد

 

• • الخطوات الأساسية لمعالجة الصورة

 

• التعامل مع البيانات التي تمت معالجتها في برامج متنوعة

دراسة لتوسع مدينة القاهرة بواسطة صور الأقمار الصناعية وبرنامج QGIS

في هذا الدرس سنقوم بعملية تصنيف باستخدام خوارزميات مختلفة باستخدام بيانات مختلفة باستخدام اداة تصنيف جديدة باستخدام برنامج الـ QGIS باستخدام الإضافة Semi-Automatic Classification تحت الرابط التالي

وهو ملحق ل Qgis يسمح بالتصنيف شبه التلقائي لصور الاستشعار عن بعد. كما أنه يوفر العديد من الأدوات لتحميل الصور المجانية (لاندسات، سنتينل-2، سنتينل-3، أستير، موديس)، والمعالجة المسبقة للصور(انشاء موزاييك وعرض الصور انشاء تراكيب من القنوات ، المعالجة اللاحقة للتصنيفات،

وسأقوم في هذا الدرس بعمل مشروع مراقبة لمدينة القاهرة ودراسة توسع المدينة من عام 1984 حتى عام 2017 باستخدام أرشيف القمر لاند سات

مراقبة التغيرات المكانية عبر الزمن باستخدام القنوات اللونية

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

مقالة تعليمية: كيفية حساب فيضان نهر النيل بواسطة برنامج QGIS وصور الأقمار الصناعية

مقدمة

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

سلسلة تعليمية – QGIS – الجزء الرابع – تحليل شريحة البيتزا

في الدرس التالي سنقوم بعمل مدخل بسيط للاستشعار وكل ما تحتاجه هو Orfeo Toolbox وهي مكتبة مفتوحة لمعالجة صور الاستشعار عن بعد من وكالة علوم الفضاء الفرنسية CNES و يمكن تنصيبها على منصات مختلفة.

سلسلة تعليمية – QGIS – الجزء الثالث – تحليل شبكات الطرق في مدينة المالكية في سورية وإيجاد أفضل مسار.

مقدمة

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

سلسلة تعليمية – QGIS – الجزء الثاني – حساب مساحة الأشجار المثمرة المغمورة خلال الفيضان؟

مقدمة

بعد التعرف على أهمية نظام المعلومات الجغرافية في المقال السابق: سلسلة تعليمية – QGIS – الجزء الأول – كيف تختار الموقع الأفضل لمنزلك؟. دعونا نخطو خطوة الي الأمام ونتعرف علي فعالية البرنامج في بعض التحليلات الهندسية كإيجاد مساحات التداخل بين طبقات مختلفة وإظهار الصفات الهندسية لعناصر هذه الطبقات وهي من المهمات الشائعة التي تواجها من خلال عملك في مجال الـ GIS (مثلا: توسع طريق او توسيع مخطط تنظيمي ) ففي حالة توسيع طريق فيجب عليك حساب المساحات الذي يجب تملكها لتوسيع هذا الطريق ونوع هذه الاراضي التي سيتم تملكها ومساحتها .

في الخريطة أعلاه مثال مشابه لهذه الحالات حيث تظهر الخريطة فيضان نهر وغمره لمناطق شجرية مغطاة باشجار متساقطة الأوراق وأشجار دائمة الخضرة. في هذا المثال سنقوم بحساب تداخل فيضان نهر Koyukuk River في الولايات المتحدة الأميركية مع المناطق المحيطة فيه وحساب مساحة المناطق المغمورة حسب نوع التغطية النباتية وسنحسب نسبة الارضي المعرضة للفيضان و المغطاة بأشجار متساقطة الأوراق ونسبة المناطق المعرضة للفيضان والغطاة بنباتات دائمة الخضرة

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

هذا المثال والبيانات بناءا على أمثلة موجودة في الكتاب المرفق.

البيانات:

• طبقة الانهار الكبيرة .
• طبقة الاشجار

خطوات العمل

نقوم بتحميل طبقة الانهار Majrivers بالضغط على الرمز الذي يظهر عليه مؤشر الماوس ونقوم بتحميل طبقة الانهار

Metadata

نتأكد ان اسقاط المشروع متوافق مع نظام الطبقة وننتبه الى ان الواحدات المستخدمة في الطبقة هي القدم (القسم الذي تم تاشيره باللون الازرق بموشر الماوس )

2-نقوم باختيار نهر Koyukuk River باستخدام select by Attribute

وذلك من processing tools

بان نكتب بالفلتر select by Attribute

processing tools: select by Attribute

فتظهر النافذة التالية ونقوم بوضع طبقة الانهار ونختار

الصفة التي سيتم الاختيار بموجبها وهي بمثالنا DESCRIPTION

نختار العملية هنا مساوة

في الحقل value نكتب القيمة التي نبحث عنها

فيتم تاشير جميع العناصر التي تمثل نهر Koyukuk River وعددها 506 عنصر

نقوم بعمل Fixed distance buffer للعناصر المحددة من خلال الذهاب الى

 Vector:Geoprocessing Tools: Fixed distance buffer

وندخل المسافة 5000 قدم

ونفعل Dissolve result

سنقوم بهذه الخطوة بحساب المساحة والمحيط لل Buffer بتشكيله حول النهر والذي يمثل منطقة الفيضان باستخدام اداة هامة هي Export/Add geometry columns وهي تفيد لحساب الخصائص الهندسية مباشرة لا مضلع موجود بالطبقة

ونذهب الى processing tools ونكتب بالفلتر Export/Add geometry columns 

فتظهر النافذة التالية

calculate using يساوي layer CRS حتى يقوم بحساب المساحة وفق لنظام إسقاط الطبقة ونستطيع حسابها وفقا لنظام اسقاط المشروع أيضا

نقوم بفتح open Attribute Table للطبقة الناتجة ولطبقة الـ Buffer

بالذهاب الي open Attribute Table لطبقة ال Buffer

ولطبقة Add geom info

ونقارن

فنلاحظ ظهور عمودين للبيانات في الطبقة الناتج احدهما يمثل المساحة والاخر يمثل المحيط. الصورة التالية تظهر البيانات الجديدة و البيانات القديمة ونلاحظ انه قد تم اضافة حقل للمساحة وحقل للمحيط.

نقوم بهذه المرحلة بإضافة طبقة الغابات

نقوم بحساب التقاطع بين طبقة add geom info وطبقة الغابات

vector:Geoprocessing Tools: intersection

ونقوم بفحص البيانات الناتجة عن تقاطع هذين الطبقتين من خلال فتح open Attribute table لطبقة التقاطع فنلاحظ أن البيانات هي كالتالي

حيث لدينا خمس عناصر ويظهر لدينا العمود VEGDESC بأن العنصر الأول والثاني أشجار متساقطة والعنصر الثالث والرابع والخامس أشجار دائمة الخضرة

لذلك سنقوم بعمل مضلع واحد فقط لكل صفة من صفات العمود VEGDESC فسنقوم بتوحيد المضلعات التي لها صفة واحدة من خلال الذهاب الى:

Vector:Geoprocessing Tools:Dissolve

فتظهر النافذة التالية

ومن المهم هنا ازالة تفعيل Dissolve all

واختيار الصفة التي سيتم توحيد العنصر وفقا لها وهي في مثالنا VEGDESC

وذلك من خلال سحبها الي الخانة اليمين من الجدول (selected) كما في الصور التالية

بعد ذلك تفحص البيانات الناتجة عن هذه العملية بالذهاب الي open Attribute table للطبقة

فنلاحظ انه لدينا الان فقط مضلعين احدهما للأشجار المتساقطة والاخر للأشجار المختلطة ومساحة كل منطقة

اخيرا نقوم بحساب نسبة المساحة المغطاة بكل نوع من انواع الشجر الي المساحة الكلية (مساحة النهر + مساحة الفيضان

بالذهاب الى

 Advanced python calculator :processing Toolbox

ونقوم بوضع اسم العمود الذي سيخرج نسبة نوع الشجر

Result file name  وقد قمت بتسميته Martin

Field type = float

ونكتب الفورمولا التالية

value = $geom.area()/ <area>

حيث التابع:

 value = $geom.area()/ <area>

يمثل مساحة المنطقة المشجرة للاشجار المتساقطة والمساحة المشجرة للاشجار المختلطة

يمثل المساحة الكلية

فنحصل علي النتائج التالية

بان 18.9% نسبة الاشجار المتساقطة الاوراق في منطقة الفيضان

6.9 % نسبة الأشجار دائمة الخضرة

نسبة الأشجار الكلية في منطقة الفيضان علي جانبي النهر 18.9+6.9 وتساوي 25.8%

وهذا ماتظهره الصورة التالية بوضع شفافية 50% لمنطقة الفيضان

ملاحظة

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

سأقوم لاحقا باتتمة المشروع واضافة ملف اتمتة المشروع عندما نتطرق لدروس نمذجة المشاريع وعمل geo-algorithms

لذلك لاتنسى الاشتراك في صفحة الفيسبوك That GIS Guy – نظم المعلومات الجغرافية ومشاركة  هذه المقالة مع أصدقائك المختصين. وفي حال لديك أي تعليقات أو أسئلة لاتتردد حيث أنك تفيدنا في تحسين المقالات.

الكاتب: المهندس مارتن إيليا

المراجع:

QGIS: Becoming a GIS Power User – Page 106 

سلسلة تعليمية – QGIS – الجزء الأول – كيف تختار الموقع الأفضل لمنزلك؟

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

هذا المثال والبيانات بناءا على أمثلة موجودة في الكتاب المرفق

المتطلبات:

• برنامج QGIS وتستطيع الحصول عليه من هذا الرابط: الرابط
•  بيانات ملائمة لأجراء هذه التحليل وسنقوم بدروس لاحقة بشرح انواع البيانات وطرق عرضها.

الشروط

• موقع قريب من بحيرة حتى 1000 متر
• موقع قريب من مدرسة ابتدائية حتى 500 متر
• موقع قريب من جامعة حتى 2000 متر
• موقع بعيد عن الشارع الرئيسي حتى 1000 متر

خطوات العمل:

نقوم بفتح البرنامج وإضافات البيانات. البيانات هي عبارة عن ثلاث طبقات:

• طبقة للمؤسسات التعليمية في نورث كارولينا في أميركا schools_wake.shp.
•   طبقة للبحيرات lakes.shp
• طبقة للشوارع الرئيسية  roadsmajor.shp

ويتم فتحها من القائمة Layer : add layer: add vector layer

طبعا عند القيام بأي مشروع قم بفحص البيانات ماهو نظام الاسقاط المستخدم وماهي الـ metadata للبيانات حيث يظهر لديك نوع البيانات و روابط التخزين وعدد العناصر والامتداد للطبقة والوحدات المستخدمة

لا تنسى ضبط نظام الاسقاط للمشروع EPSG  بحيث يتوافق مع نظام اسقاط البيانات وهو في مثالنا :NAD38(HARN)/north carolina EPSG 3358

وكذلك عليك تفحص البيانات وماذا تحتوي هذه الطبقات من بيانات مثلا طبقة المؤسسات التعليمية school_wake. نقوم باختيار الطبقة من نافذة layer panal نضغط بالزر اليميني ونختار open attribute table

نلاحظ ان الجدول يحتوي علي بيانات المؤسسات التعليمة الاسم والكو د ونوع المؤسسة …. الخ وتكون عبارة عن ارقام او كودات مرتبطة مع جدول يشرح دلالة هذه الارقام والكودات. ليس من الضروري فهم جميع محتويات هذه البيانات في هذه المرحلة.

نقوم لان بالتعرف علي طريقة اختيار عناصر محددة من طبقة البحيرات. ندخل علي قائمة الـ processing toolbox ونختار Select by expression. نضغط عليها فتظهر نافذة نقوم باختيار الطبقة وندخل صفات العناصر التي نريد ادخلها وهنا سأختار البحيرات الكبيرة التي مساحتها أكبر من 100 هكتار والتي لها نمط lake/pond.

ثم  ندخل التعبير التالي في حقل Expresion

"AREA">1000000 AND "FTYPE"= 'LAKE/POND'

ملاحظة:

QGIS حساس لعلامات التنصيص لذلك قبل المساواة (“) وبعد المساواة (‘)

نقوم بالضغط موافق فتتشر جميع البحيرات التي مساحتها اكبر من 100 هكتار ولها نمط lake/pond وعدد العناصر الناتجة 81 عنصرا من أصل 15279.

بعد ذلك نقوم بإنشاء مجال Buffer حول هذه البحيرات حيث سيكون ضمن هذا النطاق الموقع الذي اريد اختياره نذهب الى قائمة vector:Geoprocessing Tools :fixed distance buffer 

نقوم باختيار الطبقة وادخال المسافة 1000 متر واختيار dissolve result  ومكان حفظ الطبقة.

فنحصل علي الخريطة التالية:

نقوم الان باختيار المدارس الابتدائية من طبقة school_wake حيث نقوم باختيار Select by expression من خلال نافذة Processing toolbox.  نكتب التعبير التالي:

"GLEVEL" = 'E'

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

نقوم باختيار الجامعات الان من طبقة school_wake:

Select by expression:

"GLEVEL" = 'H'

ونقوم بعدها بتوليد Buffer حول هذه العناصر لكن 2000 متر

ومن ثم نقوم بعمل Buffer حول  طبقة roadmajor.shp بمسافة 500 متر

فنحصل علي الخريطة التالية

الخطوة التالية ستكون ايجاد المنطقة المطلوبة وهي عبارة عن المنطقة الي تكون قريبة من البحيرات وقريبة من الجامعات وقريبة من المدارس الابتدائية وبعيدة عن الطرقات الرئيسية لذلك نأخذ تقاطع هذه المناطق جميعا

مرحلة اولى

سنأخذ تقاطع المناطق المحيطة بالمدراس الابتدائية والجامعات

ثم نقوم بإيجاد تقاطع الطبقة الناتجة عن التقاطع السابق مع المنطقة المحيطة بالبحيرات. نقوم باجراء تقاطع اخر

مرحلة ثانية

ناخذ فرق تقاطع المناطق الثلاث (جامعات بحيرات ومدارس ابتدائية ) عن الطرقات العامة

ثم نأخذ الفرق بين ناتج التقاطعين السابقين والـ Buffer المحيط بالطرقات

الخريطة النهائية الناتجة

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

لاتنسى الاشتراك في صفحة الفيسبوك  That GIS Guy – نظم المعلومات الجغرافية ومشاركة  هذه المقالة مع أصدقائك المختصين. وفي حال لديك أي تعليقات أو أسئلة لاتتردد حيث أنك تفيدنا في تحسين المقالات.

الكاتب: المهندس مارتن إيليا

Martin Elia

المراجع:

هذا المثال والبيانات بناءا على أمثلة موجودة في كتاب:

QGIS 2 Cookbook – Page 126