$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
تم تصميم الطريقة شبه الآلية لرسم خرائط وتصنيف مجتمعات النباتات لتوليد معلومات دقيقة عن الغطاء الأرضي للأراضي الرطبة الساحلية، التي تهيمن عليها نباتات عشبية منخفضة النمو غير متجانسة وتتميز بالميكروتوبوغرافيا؛ انظر مثالا على مثل هذا الموطن في الصورة في الشكل 1. تم إجراء دراسة حالة للأراضي الرطبة الساحلية الشمالية في منطقة البلطيق في جزيرة هيوما، إستونيا. في الموقع، تم أخذ عينات من أربع مجتمعات نباتية لتدريب نموذج التصنيف: الرواد المفتوح (OP)، الشاطئ السفلي (LS)، الشاطئ العلوي (الولايات المتحدة)، والعشب الطويل (TG). يوفر الشكل 2 نظرة عامة موجزة على المنهجية المنفذة في أربع مراحل مميزة. على النقيض من ذلك، يوضح الشكل 3 دليلا مخططا مفصلا يحدد الخطوات المرقمة والقابلة للتنفيذ مقترنة بأسهم لتمثيل تدفق البيانات، وبالتالي يتم تلخيص المنهجية عبر خطوات ومراحل البروتوكول. المخرجات النهائية لهذه الطريقة، باستخدام حساس متعدد الأطياف ومستشعر RGB، هي خرائط مجتمعات نباتية متعددة الأطياف وRGB، ويتم تخزينها كطبقات في مشروع نظم المعلومات الجغرافية مع جميع مجموعات البيانات الوسيطة. تم تقييم جودة مجموعات البيانات الجغرافية المكانية عند الانتهاء من كل مرحلة بروتوكول، حيث تعتمد جودة النتائج في المراحل التالية على جودة نتائج المرحلة السابقة، والتي تبلغ ذروتها في المرحلة النهائية، حيث تشكل المدخلات مجموعا لجميع الخطوات السابقة.
فيما يلي ملخص نتائج كل من المراحل الأربع للبروتوكول:
في المرحلة الأولى، أجري المسح الجوي للطائرات بدون طيار على ارتفاع 120 مترا، مما أدى إلى مسافة عينة أرضية (GSD) تبلغ حوالي 10 سم/بكسل. أنتج المسح الجوي صورا متعددة الأطياف وصور RGB. بعد معالجة الحركية اللاحقة (PPK)، تم تصحيح العلامات الجغرافية لكل صورة طائرة بدون طيار، وتم تغليف الصور بفسيفساء باستخدام برامج تصوير فوتوغرامترية احترافية لإنتاج خرائط مسح جوي. أسفرت المعالجة الضوئية للصور الجوية متعددة الأطياف عن ملفات نقطية ذات انعكاس أحمر، أخضر، NIR، وأحمر الحواف.
تمت معالجة الصور الجوية متعددة الأطياف وRGB كمشاريع منفصلة بسبب اختلاف دقة الصور الفوتوغرافية بشكل متبيع. انظر الشكل 4 لإعدادات فسيفساء صور RGB. الشكل 4A. خيارات المعالجة من خلال اختيار البلاطات المدمجة وتوليد خرائط مخرجات موحدة وسلسة. الشكل 4B. اختيار حسابات مؤقتة للمؤشر من صور كاميرا RGB أنتج أربع خرائط نقطية: واحدة لكل نطاق (أحمر، أخضر، وأزرق)، ونموذج DSM (نموذج السطح الرقمي). الشكل 4C. يوضح توليد نقاط ربط تلقائية لمطابقة البكسلات الشائعة عبر الصور المتداخلة خلال خطوة المعالجة الأولى. مع 3 مؤشرات RGB (أحمر، أخضر، أزرق) وDSM، أنتج برنامج الفوتوغرامتري نقطيا رماديا مولدا بالبرمجيات: تمثيل أحادي اللون أحادي القناة يجمع ويجمع بين النطاقات المختلفة (مثل الأحمر، الأخضر، الأزرق)، مع ظهور قيم بكسل منخفضة باللون الأسود والقيم العالية باللون الأبيض. صياغة هذا المؤشر الرمادي محددة للبرمجيات وليست موحدة عبر منصات التصوير الفوتوغراميتي.
أنتجت برمجيات التصوير الفوتوغرافي تقارير جودة تصف نتائج المعالجة لكلا المجموعتين من مشاريع تصوير الفسيفساء الجوية: متعددة الأطياف وRGB. توفر هذه التقارير مؤشرات رئيسية للتحقق من الدقة. مؤشرات النتائج عالية الجودة هي انخفاض خطأ الإسقاط: كلما اقتربت من 0 كان ذلك أفضل (أي شيء أقل من 1.0 مقبول عموما). على سبيل المثال، تعني خطوط "مطابقة" كثيفة تطابقات قوية بين معظم الصور. تم تقييم خرائط المخرجات بصريا ووجدت خالية من الفراغات والفروقات البصرية الأخرى. التناقضات البصرية هي علامة على احتمال حدوث أخطاء في الانعكاس والفسيفساء، بالإضافة إلى وجود مشاكل في شبكة المرجع الجغرافي. يستخدم تقرير الفوتوغراميترية، إلى جانب كتيبات البرمجيات المحددة وموارد التدريب، لتحديد المشكلة، وإيجاد حل، وإعادة معالجة الصور الجوية.
خدمت خرائط المسح الجوي للطائرات بدون طيار كمجموعة بيانات أساسية لإعادة حساب المعلومات الطيفية المعقدة إلى مؤشرات بكسلات موحدة (أي مؤشرات النباتات السادسة التي تولدها المرحلة الثالثة)، مما أتاح استخدامها كمدخلات منظمة لتصنيف التعلم الآلي (في المرحلة النهائية من البروتوكول).
تم جمع الملاحظات الفيزيائية الحيوية في الموقع كمجموعة بيانات من الحقيقة الأرضية. في موقع دراسة الحالة (الأراضي الرطبة الساحلية في البلطيق بوريل)، تم تحديد وأخذ عينات من أربع مجتمعات رئيسية من النباتات الوعائية: الرواد المفتوح (OP)، الشاطئ السفلي (LS)، الشاطئ العلوي (الولايات المتحدة)، والعشب الطويل (TG). باستخدام طريقة أخذ عينات الطبقات الرباعية التي تم تحديدها في المرحلة الثانية من البروتوكول، تم تحديد مجتمع النباتات السائدة من خلال وجود أنواع مؤشرات: ساليكورنيا أوروبيا وسويادا ماريتيما؛ LS بنسبة عالية من Glaux maritima وJuncus geradii؛ الولايات المتحدة حسب أنواع المؤشر: Festuca rubra و Leontodon autumnalis؛ بينما يتميز TG بوجود أعشاب طويلة بشكل ملحوظ: Deschampsia cespitosa، Elytrigia repens، وMolinia caerulea. تم جمع ما لا يقل عن 10 عينات من الرباعيات لكل مجتمع نباتي. تضمنت عينة كل ربع ملاحظات حيوية فيزيائية كالتالي: تركيب أنواع النباتات (٪)، رطوبة التربة (٪)، ارتفاع النباتات بالسنتيمترات، ومواقع إحداثيات الربع الدقيقة (quadrat). تم رقمنة هذه المتغيرات كمجموعة بيانات مجمعة، ثم استوردت إلى مشروع نظم المعلومات الجغرافية، وأضيفت إلى ملف شكل بيانات GNSS المربعة العينة، وحفظها كحزمة جغرافية (.gpkg). مثل هذا الملف المستقل عن المنصة، والذي يحتوي على ملاحظات بيوفيزيائية لمجتمعات النباتات الرطبة وبيانات إحداثيات دقيقة للعينات، يخدم غرضا مزدوجا: تقييم ونمذجة البيئية البيئية والبيئية بشكل أعمق. والأهم من ذلك، من أجل التنفيذ الناجح لهذا البروتوكول، يعمل ملف GeoPackage هذا كبيانات تدريب نموذجي وتحقق في عملية التعلم الآلي المحددة في المرحلة الرابعة. ترميز مجتمع النباتات: OP، LS، US، TG؛ تم استخدامها كتسميات تدريبية في عملية التعلم الآلي. الدقة (تسميات عالية الدقة، تدقيق بشري، تنسيق متسق) والبنية (تنسيق نظيف ومستطيل (صفوف/أعمدة)) ضرورية لإنتاج مجموعة بيانات صحة أرضية عالية الجودة. يجب أن تكون معدلات الخطأ شبه الصفرية (مثل الأخطاء المطبعية) لضمان اتساق قواعد التصنيف. على سبيل المثال، إذا كان "OP" هو الحقيقة الأساسية، فإن التأكد من عدم تصنيفها بأي طريقة أخرى (مثل "op" أو "O.P.") أمر ضروري لتجنب التكرار.
لتنفيذ المرحلة الثالثة من البروتوكول، تم استخدام مجموعات البيانات من المسح الجوي الذي أجري في المرحلة الأولى لإعادة حساب بيانات المسح الجوي الخام إلى مؤشرات النباتات (VIs). أولا، لتقليل حجم مجموعة البيانات وتبسيط عملية الحوسبة، يتم "قص" مخرجات الفوتوغرامتري لحذف بيانات البكسل غير ذات الصلة أو منخفضة الجودة خارج حدود موقع الدراسة (انظر الشكل 5 كتوضيح للإجراءات التي تم تنفيذها في برنامج نظم المعلومات الجغرافية).
ثم تم حساب 19 مقياس متعدد الأطياف باستخدام خرائط الانعكاس الحمراء والخضراء والنورية النارية والحواف الحمراء. تم إنشاء ملفات النقطي متعددة الأطياف VI التالية وأدرجتها في مشروع نظم المعلومات الجغرافية: مؤشر دات 4 (Datt4)41، مؤشر النباتات المحسنة (EVI)42، مؤشر الكلوروفيل الأخضر (GCI)43، مؤشر فرق الخضرة (GDI)44، مؤشرات الفرق المعمم للنباتات (GDVI)45، مؤشر نسبة الغطاء تحت الأحمر الأخضر (GIPVI)46، مؤشر الفرق الأخضر النباتي المعتدل (GNDVI)47، مؤشر الفرق الأخضر-الأحمر (GRDI)44، مؤشر النباتات الخضراء-الحمراء (GRVI)45، مؤشر الفرق المعدل للنباتات (mNDVI)48، مؤشر النباتات المعدلحسب التربة 49 (MSAVI: انظر الشكل 6. لقطة شاشة لحاسبة نقطية في برنامج نظم المعلومات الجغرافية مع إدخال صيغة MSAVI)، النسبة البسيطة المعدلة للحافة الحمراء (MSRred_edge)50، مؤشر الفرق الطبيعي للماء (NDWI)51، مؤشر الفرق النباتي المعدل (NDVI)52، مؤشر الحافة الحمراء NDVI (NDVIre)47، مؤشر النباتات المثلثة ذات الحافة الحمراء (RTVIcore)50، مؤشر النباتات المعدلة للتربة (SAVI) 53،54، النسبة البسيطة (SR)55، نسبة الحافة الحمراء البسيطة (SRre)47.
علاوة على ذلك، تم حساب وحفظ ما مجموعه 23 مؤشر VI قائم على RGB كملفات نقطية في مشروع نظم المعلومات الجغرافية لمزيد من التحليل المكاني: مؤشر النسبة الأزرق–الأخضر (BGRI)57، مؤشر السطوع (BI)58، مؤشر الغطاء النباتي الأزرق الواسع النطاق الديناميكي (BRVI)59، مؤشر اللون للغطاء النباتي (CIVE)60، مؤشر الرؤية الخضراء المحسنة للنباتات (EGVI)61، مؤشر الغطاء النباتي المعزز (ERVI)62، مؤشرات الخضراء الزائدة ExG و ExGI63، مؤشر التركيبة (COM)64، الإحداثيات اللونية الخضراء (GCC)34،65، مؤشر الأوراق الخضراء (GLI)66، مؤشر الغطاء النباتي الأخضر-الأحمر (GRVI)67، مؤشر النباتات الخضراء الحمراء المعدل (MGRVI)68، مؤشر الفرق الأخضر-الأزرق المعدل (NGBDI)69، مؤشر الفرق الأخضر-الأحمر المعدل (NGRDI)69، مؤشر نسبة الأحمر-الأخضر-الأزرق ومؤشر الغطاء النباتي الأحمر-الأخضر-الأزرق (RGBRI70 ) ومؤشر RGBVI68)، مؤشر النسبة الأحمر إلى الأخضر (RGRI)57، مؤشر التربة المعدلة للنباتات المرئية (SAVI)53، مؤشر الخضرة المثلثية (TGI)69، مؤشر المقاومة الجوية المرئية (VARI)47، المؤشر النباتي (VEG)71، مؤشر ووبيكي (WI)63.
تشمل المخرجات النهائية خريطتين منفصلتين لتصنيف مجتمع محطات متعددة الأطياف ومحطة RGB محفوظتين كبيانات GeoTIFF ضمن مشروع نظم المعلومات الجغرافية، إلى جانب جميع البيانات الجغرافية المكانية السابقة. تم استخدام 19 مقياس VI متعدد الأطياف، إلى جانب مجموعة بيانات الحقيقة الأرضية، كمتنبئين على مستوى البكسل في سير عمل التعلم الآلي (مثل خوارزمية التردد الراديوي) لتصنيف مجتمعات النباتات الساحلية الشمالية (OP, LS, US, TG) وتوليد خريطة لمجتمع النباتات متعددة الأطياف. وأيضا متنبعات مبنية على RGB (25 مؤشر RGB VI، مؤشر رمادي،و56 مؤشرا أحمر وأزرق وأخضر وDSM مصحح)، إلى جانب مجموعة بيانات الحقيقة الأرضية، لتصنيف وتوليد خريطة مجتمعية نباتية قائمة على RGB (انظر الشكل 7). المخرجات النهائية التي تم رسمها في وحدة تحكم R؛ الشكل 7A يمثل خرائط متعددة الأطياف، والشكل 7B يمثل خرائط مجتمعات النباتات القائمة على RGB).
لتوليد المخرجات النهائية، تم استخدام بيئة R (المرحلة 4 من البروتوكول) لبرمجة وتنفيذ سير عمل التعلم الآلي (مثل خوارزمية التردد الراديوي)، مما أتاح تصور المخرجات الوسيطة، مثل تصنيفات أداء مؤشر النباتات (VI) (الشكل 8). عند رسم خرائط وتصنيف مجتمعات النباتات، قدم وحدة تحكم R رسائل أخطاء وتفسيرات في الوقت الحقيقي، مما سهل استكشاف الأخطاء وضمان التقدم بسلاسة في كل خطوة.
قامت هذه الدراسة بتقييم نموذجين لخوارزمية الترددات الراديوية لخرائط تصنيف مجتمعات النباتات. حقق النموذج الأول للترددات الراديوية، الذي استخدم المؤشرات المتعددة الأطياف كمجموعة بيانات إدخال، دقة 92.34٪ مع درجة F1 بلغت 0.915 في مجموعة بيانات التحقق. هذه نتيجة صحيحة، حيث تتراوح درجات الفورمولا 1 من 0 (الأسوأ) إلى 1 (الكمال)، مع الدرجات الأعلى التي تشير إلى أداء متفوق. تمثل درجة F1 في نموذج الترددات الراديوية المتوسط التوافقي للدقة والاسترجاع. هذا المقياس يوازن فعليا بين الإيجابيات الكاذبة والسلبيات، أي أنه لا يسمح لنوع من الأخطاء بأن يتفوق على الآخر، وهو أمر أساسي لمجموعات البيانات غير المتوازنة، وغالبا ما يتم احتساب متوسط بين الفئات أو الإبلاغ عنه في كل فئة لتقييم الأداء العام للنموذج78,79. بلغ معدل خطأ OOB في النموذج متعدد الأطياف 7.75٪، مما يشير إلى خطر متوسط لسوء التصنيف. تراوحت درجات F1 على مستوى الصف في مجموعة بيانات التحقق من 0.845 في TG إلى 0.984 في مجتمعات النباتات السائدة، مما يشير إلى تباين متوسط في الأداء عبر المجتمعات.
أما النموذج الثاني للراديو الراديوي الذي تم تقييمه فكان مبنيا على الرؤية البصرية والDSM المستمدة من المسح الجوي RGB، محققا دقة بلغت 98.89٪ مع درجة F1 إجمالية بلغت 0.987. كان معدل خطأ OOB فقط 1.14٪، مما يعكس موثوقية عالية جدا. كانت درجات F1 على مستوى الصف مرتفعة باستمرار (0.986–0.993)، مما يشير إلى تصنيف قوي لجميع المجتمعات.
أظهر مخطط الأهمية المتغيرة في R (الشكل 8A) قدرات كل من الورشة البصرية التنبؤية لمجتمعات النباتات الرطبة الساحلية في البوريال. كانت أكثر الظروفان البصرية متعددة الأطياف تأثيرا كمؤشرات هي mNDVI وSR، والتي تشير إلى حساسية قوية لخضرة المظلة والكتلة الحيوية. وبالمثل، كانت GRDI وNDVI وNDVIre وSRre أيضا مرتبة عالية، مما يعكس دور مؤشرات الفروق ذات الحافة الحمراء والتطبيع في تمييز أنواع النباتات. على النقيض من ذلك، يظهر مؤشر RGB ومخطط الأهمية المتغيرة في DSM أن نموذج الارتفاع الهيكلي المستخرج من DSM كمتنبئ هو عامل رئيسي لدقة التصنيف (الشكل 8B). تشمل المتنبئات المهمة الأخرى TGI، والتدرج الرمادي، وVARI، التي تستفيد من تباين انعكاس RGB. لذا، تشير هذه النتائج إلى أن المؤشرات المشتقة من DSM وRGB تحسن بشكل كبير موثوقية التصنيف وتقلل من خطر سوء التصنيف.
تم حفظ مجموعة البيانات الجغرافية المكانية بالكامل كمشروع نظم المعلومات الجغرافية، وهو مجلد يشمل خرائط مسح الطائرات بدون طيار وخرائط الراستر البصرية، وحزمة المراقبة الحيوية GeoPackage (الحقيقة الأرضية)، وخرائط مجتمع النباتات؛ كلها مرتبطة كمشروع QGIS (ملف .qgz). المشروع مناسب للاستخدام غير المتصل ويتيح التوافق بين منصات نظم المعلومات الجغرافية.

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

الشكل 2: مخطط سير العمل المفاهيمية. (1) المسح الجوي القائم على الطائرات بدون طيار، (2) جمع بيانات الحقيقة الأرضية والمرجعية الجغرافية، (3) حساب مؤشر النباتات، و(4) التصنيف الخاضع للإشراف باستخدام الغابة العشوائية في R. توضح الخرائط الفعلية من تنفيذ البروتوكول هذه الطريقة ذات المراحل الأربع. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.

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

الشكل 4: مثال على برنامج التصوير الضوئي. (A–C) لقطة شاشة لمشروع فسيفساء صور RGB: (أ) خيارات المعالجة ل DSM وorthomosaics. (ب) معالجة خيارات المؤشرات الحمراء والخضراء والزرقاء. (ج) عرض تقدم المعالجة، بما في ذلك الصور الجوية وزوايا الكاميرا. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 5. قص ملفات النقطي على الخرائط في برنامج نظم المعلومات الجغرافية. لقطة شاشة لمساحة عمل QGIS مع طبقات مرئية وعلامات تبويب مفتوحة بالتتابع توضح كيفية قص البيانات غير الضرورية: استخراج نقطي-قص نقطي بحسب Extend ("extent" في هذه الحالة تشير إلى شكل مضلع مرسوم كما هو مرغوب ويستخدم لإخفاء البكسلات التي تسقط خارج حدوده). يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.

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

الشكل 7. بيئة الحوسبة الإحصائية لصور شاشة لإصدار R تظهر النتائج النهائية المخططة في وحدة تحكم R (RStudio). (أ) خرائط مجتمعات نباتية متعددة الأطياف و(ب) خرائط مجتمعات نباتية تعتمد على RGB. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 8. بيئة الحوسبة الإحصائية للقطات شاشة R التي توضح تصنيفات الأداء VI. (أ، ب) مشاهد لوحدة تحكم R: (أ) لقطة شاشة عند استخدام الشبكات المتعددة الأطياف كمجموعة بيانات إدخال، و(ب) لقطة شاشة عند استخدام مجموعة بيانات RGB+DSM. تتضمن كل من العرضين (A وB) مخططات ذات أهمية متغيرة تسمى "final_rf" (اسم كائن مولد بواسطة السكريبت). يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
الملف الإضافي 1: معادلات مؤشر النباتات. قائمة بمعادلات مؤشر النباتات (VI) المعتمدة على الأربعة الطيف وRGB المستخدمة في هذا البروتوكول، بما في ذلك صيغ تطبيع نطاق RGB قبل حساب RGB VI. يرجى الضغط هنا لتحميل هذا الملف.
الملف التكميلي 2: سكريبت تصنيف والتحقق من RF. يرجى الضغط هنا لتحميل هذا الملف.