1. مقدمه
جنگلها اکوسیستمهایی هستند که منابع زیستی را حفظ میکنند و طیف وسیعی از محصولات و خدمات را برای رفع نیازهای انسان ارائه میکنند [ 1 ، 2 ]. همانطور که گزارش شده است، 31 درصد از کل سطح زمین توسط جنگل ها پوشیده شده است که نقش مهمی در حفظ اکوسیستم ها و تامین غذا، علوفه و سرپناه دارند. جنگل ها همچنین در کنترل تخریب زمین با به حداقل رساندن فرسایش خاک و بیابان زایی نقش موثری دارند. علاوه بر این، حتی تنوع زیستی در یک اکوسیستم منطقه ای به طور مستقیم تحت تأثیر پوشش جنگلی اطراف قرار می گیرد [ 3 ]. زیست توده جنگلی به عنوان مخازن کربن در نظر گرفته می شود. جنگل ها به طور مستقیم نقش مهمی در چرخه کربن و آب، تولید اولیه ناخالص و خالص ( GPP و NPP ) دارند.، تغییرات آب و هوایی و گرم شدن کره زمین، تشعشع و تعادل محیطی، حتی کیفیت هوا، و فعالیت های انسانی [ 1 ، 4 ]. ارزیابی دقیق زیست توده جنگل در بالای زمین ( AGB ) برای ارزیابی منابع جنگلی، ارزیابی چرخه کربن و مدیریت پایدار جنگل بسیار مهم است [ 1 ، 2 ، 3 ، 4 ]. بنابراین، ارزیابی AGB جنگل برای شناسایی تغییرات احتمالی آب و هوا به دلیل جنگلزدایی و تعادل اکولوژیکی جهانی، و همچنین ترویج تکامل بیولوژیکی جهانی، جانشینی جامعه، و فرآیندهای کاهش ضروری است [ 1 ، 5 ، 6 ]]. در سطح جهانی، جنگل به سرعت در حال زوال ممکن است تهدید بزرگی نه تنها برای منابع جنگلی بلکه برای تعادل اکولوژیکی آینده باشد [ 1 ]. بنابراین، ارزیابی و پایش دقیق AGB جنگل برای مدیریت پایدار و کاهش، و همچنین تشخیص سناریو تغییرات آب و هوایی به دلیل جنگلزدایی مهم و حیاتی است [ 1 ، 5 ، 6 ]. جنگلها پتانسیل زیادی برای سوخت نسل بعدی با برداشت محصولات جانبی (مانند چوب برس، الوار، هیزم) و تامین بخش زیادی از مواد تجدیدپذیر (مانند میوههای خوراکی، آجیل، کاغذ، گیاهان دارویی) دارند که میتوانند به انرژیهای تجدیدپذیر تبدیل شوند. 1 ، 2 ، 3 ، 4, 7 , 8 , 9 ]. ارزیابی زیست توده جنگل برای ارزیابی منابع انرژی تجدیدپذیر و برنامه ریزی آینده برای حفظ تعادل در محیط زیست، اقتصاد و امنیت انرژی بسیار مهم است [ 1 ، 2 ، 3 ، 4 ].]. بنابراین، ارزیابی و پایش زیست توده جنگل برای ارزیابی منابع زیستی و اکوسیستم زمینی بسیار مهم است. جنگل ها به عنوان یک مخزن کربن در نظر گرفته می شوند. جنگل ها از طریق فتوسنتز کربن را از جو جذب کرده و در بافت گیاه ذخیره می کنند. در سطح منطقه ای و همچنین جهانی، برآورد زیست توده پوشش گیاهی جنگل یک تمرین مهم در تعیین ذخیره کربن در جزء درخت غالب و محاسبه چرخه کربن است. سپس در زیست توده جنگل (یعنی تنه ها، شاخه ها، ریشه ها و برگ ها)، در مواد آلی مرده (آشغال و چوب خشک) و در خاک رسوب می کند. این فرآیند جذب و رسوب کربن توسط پوشش گیاهی جنگل به عنوان ترسیب کربن شناخته می شود که می تواند به مدیریت پایدار منابع طبیعی کمک کند. برای کاهش CO 2انتشار گازهای گلخانه ای از طریق ترسیب کربن، این محرک تغییرات آب و هوایی توجه قابل توجهی را به اکوسیستم های جنگلی به عنوان یک گزینه مناسب جلب کرده است [ 10 ].
در سطح جهانی، سیستمهای مشاهده جهانی پیشرفته و سریع و تجزیه و تحلیل منابع طبیعی جغرافیایی با استفاده از تکنیکهای سنجش از دور در حال توسعه هستند. سنجش از دور تکنیکی پیشرفته، سریع، مقرون به صرفه و سازگار با محیط زیست برای پایش و نقشه برداری بسیاری از منابع طبیعی جغرافیایی و اقلیم از جمله جنگل، دینامیک آب پوشش گیاهی و نظارت بر تغییرات و ارزیابی زیست توده است [ 1 ، 6 ، 11 ، 12 ، 13 ، 14 ]. امکانات پیشرفته و باکیفیت جمعآوری و توزیع دادهها توسط سازمانهای فضایی مختلف، از جمله اداره ملی هوانوردی و فضایی ( ناسا )، مرکز علوم از راه دور کانادا (کانادا) ارائه میشود.CCRS ) و آژانس فضایی اروپا ( ESA ) [ 1 ، 15 ]. به طور کلی، تکنیکهای سنجش از دور، زیست توده جنگل را مستقیماً استخراج نمیکنند، اما به دست آوردن پارامترهایی برای تخمین زیست توده جنگل، مانند ارتفاع جنگل، شاخص سطح برگ، کسری از پوشش گیاهی ( FVC )، یا تولید اولیه خالص را امکانپذیر میسازند. امروزه به دلیل دقت نسبتاً بالاتر نسبت به سیستم های معمولی، سیستم های سنجش از دور عملکرد سریع و مقرون به صرفه در ارزیابی زیست توده جنگل و ارزیابی تغییرات توجه زیادی را به خود جلب کرده اند. روش سنجش از دور معمولاً زیست توده و تشخیص/اطلاعات تغییر را بهجای مستقیم با بازیابی پارامترهای سنجش از دور ارزیابی میکند. این پارامترهای استخراج شده به ارزیابی کمک می کنندنقشه برداری و نظارت AGB ، جنگل و پوشش گیاهی با رویکردهای یکپارچه سازی، وارونگی و مدل سازی [ 1 ، 6 ، 13 ، 14 ]. علاوه بر این، دادههای باکیفیت موجود از حسگرهای پیشرفته روی ماهوارههایی که اخیراً توسط آژانسهای فضایی مختلف به فضا پرتاب شدهاند، در مدلسازی اکوسیستمها به شیوهای نوآورانه و همچنین برای استخراج سیستمهای موجودی جدید مبتنی بر خشکی برای تخمین ارتفاع تاج جنگل، زیست توده، ذخایر کربن، و تغییرات مکانی و زمانی با دقت نسبتاً بالا [ 6 ، 13 ، 14 ]]. در دسترس بودن تصاویر Sentinel-2 با وضوح فضایی بالا (10 تا 60 متر) از برنامه کوپرنیک فرصت جدیدی برای مطالعه و ارزیابی منابع زیستی جنگلی و ذخایر کربن طبیعی ارائه می دهد. در حال حاضر، مجموعه دادههای چندگانه سنجش از راه دور از طریق منابع غیرفعال و فعال از پلتفرمهای مختلف از جمله دادههای فضابردی، هوایی و زمینی به دست میآیند. این داده ها با استفاده از رویکردهای پردازش نوآورانه با مهارت های پیشرفته آماری، ریاضی و یادگیری ماشینی برای کاربردهای ذکر شده استفاده می شوند [ 11 ، 16 ، 17 ]. معمولاً دادههای حسگرهای فعال دارای نسبت سیگنال به نویز بالا ( S:N ) و نفوذ بهتر در مقایسه با سنسورهای غیرفعال هستند. بنابراین، سیستم های فعال مانند رادار دیافراگم مصنوعی (SAR )، تشخیص نور و محدوده ( LIDAR ) می تواند به تاج پوشش جنگل ها، پوشش های گیاهی و پوشش های زمین نفوذ کند. در نتیجه، سنجش از دور فعال به ابزاری قدرتمند برای اندازهگیریهای ساختاری و حجمی عمودی جنگلها و زیستتوده تبدیل شده است. داده های حسگرهای غیرفعال پتانسیل زیادی برای شناسایی، طبقه بندی و نقشه برداری از پوشش جنگلی، انواع و طبقه بندی پوشش گیاهی دارند [ 14 ، 18 ، 19 ، 20 ].
از دهه گذشته، ارزیابی منابع جنگلی و بررسی اتمسفر به طور قابل ملاحظه ای با تکنیک های سنجش از دور تقویت شد [ 11 ]. این تکنیک ها اغلب برای تخمین جنگل AGB با استفاده از مجموعه داده های مختلف استفاده می شوند. دادههای مایکروویو مانند SAR برای ویژگیهای بافت و استخراج اطلاعات ضریب پراکندگی پسزمینه مفید هستند و دادههای نوری مانند Sentinel-2 برای استخراج اطلاعات شاخصهای گیاهی و شاخص سطح برگ ( LAI ) استفاده میشوند [ 1 ، 14 ، 21 ]. به عنوان مثال، Sentinel-1 به عنوان داده SAR Cتصاویر نوری باند و Sentinel-2 با وضوح فضایی بالا (10 تا 60 متر) به طور رایگان از برنامه کوپرنیک اروپا برای نقشه برداری و نظارت موثر پوشش گیاهی و جنگل، تخمین زیست توده و ارزیابی منابع طبیعی در دسترس هستند [ 1 ، 11 ]. به طور مشابه، تعداد زیادی مجموعه داده سنجش از راه دور وجود دارد که به طور رایگان از ناسا ، ESA و سایر آژانسهای فضایی برای کاربردهای فوق در دسترس هستند. روشهای مدلسازی و رگرسیون مختلف بین دادههای مایکروویو سنجش از دور (یعنی Sentinel-1- قطبشهای عمودی-عمودی ( VV ) و عمودی-افقی ( VH )) و دادههای نوری (یعنی Sentinel-2- LAI و شاخص گیاهی تفاوت نرمالشده (NDVI )) و داده های میدانی در حال حاضر برای ارزیابی زیست توده و اعتبارسنجی حقیقت زمینی در حال توسعه هستند [ 22 ]. برای تخمین ساختار بیوفیزیکی، دادههای مایکروویو Sentinel-1 برای استخراج تخمین ارتفاع درخت مبتنی بر InSAR از یک جنگل تپهای استفاده شد. این مطالعه کاهش عدم قطعیت تخمین ارتفاع درخت را حتی از یک زمین مواج نشان داد [ 23 ]. تصاویر نوری پتانسیل زیادی برای نظارت چند زمانی، بهبود طبقهبندی نوع جنگل با وضوح فضایی بالا، و زمان دسترسی مجدد سریع فراهم میکنند [ 24 ]. به عنوان مثال، تصاویر سنجش از دور نوری یعنی. Sentinel-2 توسط محققان مختلف در طبقه بندی، نقشه برداری و پایش جنگل استفاده شده است [ 18 ], 19 , 20 , 25 ]. رویکردها و آزمایشهای مختلفی برای استخراج پارامترهای مختلف در تخمین AGB جنگل با استفاده از مجموعه دادههای سنجش از دور نوری و مایکروویو، از جمله ارتفاع درخت، قطر درخت در ارتفاع سینه ( DBH ) و حجم درخت در دسترس است [ 25 ، 26 ]. مدلهای تجربی و معادلات آلومتریک مختلفی برای تخمین AGB جنگل از دادههای میدانی در سطح نمودار با استفاده از تکنیکهای اندازهگیری دستی، به عنوان مثال، DBH درخت ، ارتفاع، و گاهی اوقات چگالی چوب برای بهبود دقت ایجاد شدهاند [ 27 ].]. مدلهای رابطهای مختلفی بین سطح پایه یا قطر درخت، پارامترهای ارتفاع وزنی (همچنین به عنوان ارتفاع لوری شناخته میشود) و AGB برای ارزیابی دقیقتر زیست توده [ 28 ، 29 ، 30 ] و AGB جنگل مستقیم بهدستآمده از SAR و LiDAR وجود دارد. داده ها با استفاده از مدل های فوق از مناطق بزرگ [ 7 ، 31 ، 32 ، 33 ]. به طور مشابه، مطالعات سنجش از دور نوری نشان می دهد که شاخص های LAI ، FVC و پوشش گیاهی پتانسیل قابل توجهی برای ارزیابی AGB جنگل دارند [ 34 ]. راFVC و LAI در ارزیابی زیست توده برگ از طریق یک تابع تناسب مفید هستند [ 35 ]. معمولاً AGB با استفاده از داده های میدانی در سطح نمودار مجموع شاخ و برگ و زیست توده چوب را محاسبه می کند [ 36 ]. برای محاسبه LAI یا FVC ، چندین شاخص پوشش گیاهی آزمایش شده است [ 37 ، 38 ، 39 ]، که هنوز برای انواع جنگل های خاص چالش برانگیز هستند. علاوه بر این، دقت AGBتخمین از دادههای سنجش از دور نوری و مایکروویو توسط سطح دقت پردازش دادهها و استخراج پارامتر، از جمله عوامل اکولوژیکی، زیستی، و توپوگرافی محدود شده است [ 6 ، 13 ، 14 ، 22 ، 23 ، 24 ، 25 ، 26 ، 40 ، 41 . ، 42 ]. در این مطالعه، یک مدل تجربی مبتنی بر AGB میدانی و FVC از سنجش از دور نوری برای سادهسازی و ارزیابی سریعتر AGB جنگل کار شده است. اطلاعات میدانی (ارتفاع درخت، درخت DBHتعداد درختان در سطح کرت (30 در 30 متر)، تراکم ساقه درخت، و FVC برای محاسبه AGB میدانی از یک سایت آزمایشی در جنگل حفاظتشده توندی، ایالت جارکند در شرق هند جمعآوری شد.
2. مواد و روشها
برای محل آزمایش در جنگل حفاظت شده توندی، AGB جنگل با استفاده از تصاویر Sentinel-2 ارزیابی شد. شرح مختصری از محل آزمایش، جزئیات داده ها و روش پیشنهادی به ترتیب در ادامه آمده است.
2.1. منطقه مطالعه
منطقه مورد مطالعه 144 کیلومتر مربع را در داخل جنگل حفاظت شده توندی در ناحیه دنباد، ایالت جارکند، هند پوشش می دهد. این ایالت هندی به دلیل شرایط اقلیمی و فیزیوگرافی متنوع یکی از غنی ترین ایالت ها از نظر تنوع زیستی است. از جهت شمال غربی به شمال شرقی گسترش می یابد و تپه های دوموندا دارای دو سر (بالا) در حاشیه شرقی آن هستند.
در نیمه شمالی این ولسوالی چند تپه مجزا با ابعاد مختلف پراکنده شده است. این منطقه یک جنگل برگریز معمولی است که مساحت کل جنگلی آن 23605 کیلومتر مربع است که حدود 29.61٪ از کل مساحت زمین جارخند را تشکیل می دهد. در هند، جنگل های حفاظت شده و حفاظت شده توسط دولت های ایالتی مربوطه اعلام می شود. در جنگلهای حفاظتشده، همه فعالیتها ممنوع است مگر اینکه اجازه داده شود. در جنگل های حفاظت شده برعکس است، یعنی فعالیت های خاص مانند شکار، چرا و غیره مجاز است مگر اینکه ممنوع شود. از مجموع اراضی جنگلی جارخند، 28/81 درصد در زیر پوشش جنگل های حفاظت شده و 58/18 درصد در رده جنگل های حفاظت شده قرار دارد.
از نظر جغرافیایی، جارکند دارای دو دسته جنگلی اصلی است، یعنی جنگلهای خشک منطقه استوایی و جنگلهای مرطوب منطقه استوایی [ 43 ]. محل آزمایش در رده جنگل های خشک منطقه گرمسیری قرار می گیرد و با Shorea robusta (نام محلی Sal) به عنوان غالب ترین گونه درخت پوشیده شده است، در محدوده 23°45’40” تا 24°05’50” شمالی و 85 درجه و 57 دقیقه قرار دارد. 30 اینچ تا 86°35’55 اینچ شرقی ( شکل 1 a). جنگلهای حفاظتشده توندی در بخش جنگلهای دنباد در زیر گروه 5-B، یعنی جنگلهای برگریز خشک استوایی شمالی قرار دارند [ 8 ]. پوشش جنگلی ذخیره توندی حاوی محصول تقریباً خالصی از نهالها و قطبهای Shorea robusta است. محل آزمایش توسط تپه های کوچک با توپوگرافی موج دار تحت سلطه است. شوره روبوستاغالب بوده و طبقات بالای جنگل را با ارتفاع متوسط تا 12 متر در محل آزمایش تشکیل می دهد. هوموس به جز در مناطق دور افتاده که پوشش جنگلی به خوبی محافظت می شود وجود ندارد. این منطقه دارای آب و هوای معمولی از نوع موسمی با سه فصل مشخص است: گرم، بارانی و زمستان. رطوبت هوا در فصل بارندگی بسیار زیاد و در هوای گرم بسیار کم است. حداکثر دما از 45 درجه سانتیگراد بالاتر می رود و باد گرمی به نام “لو” اغلب از آوریل تا ژوئن تا زمان شروع موسمی می وزد. رعد و برق معمولاً در ماه می رخ می دهد که با کاهش موقتی چند درجه همراه است. به دلیل صنعتی شدن و فعالیت های معدنی سنگین، ذرات معلق فراوان و ذرات گرد و غبار در این فصل در جو یافت می شود. بادهای موسمی معمولاً در اواسط ژوئن رخ می دهد و تا پایان سپتامبر ادامه می یابد.
از نظر زمین شناسی، توندی شامل سنگ پرکامبرین است، سنگ کشور شیست کوارتز-فلدسپاتی است که توسط برخی اجسام آذرین، به عنوان مثال، متادولوریت ها و آمفیبولیت ها به شدت نفوذ کرده است. تپه و تپه های برجسته در این منطقه عمدتاً شامل سنگ های نفوذی تیره رنگ، سخت و فشرده (عمدتا آمفیبولیت و متادولوریت) است. خاک و آبرفت ها عمدتاً از گنیس و شیست کوارتز فلدسپاتیک به دست می آیند. تشکیل خاک در منطقه جنگلی کم عمق تا بسیار کم عمق است و عمق در دشت تا زمین های مواج به ندرت از 60 سانتی متر (24 اینچ) تجاوز می کند. به طور کلی لومی قرمز با جیب های سفالی است. فرسایش رخ می دهد و با شیب متفاوت است. خاک سطحی سخت و فشرده است [ 44 ].
2.2. مجموعه داده ها
محصول تصویر چند طیفی Sentinel-2 ( MSI ) سطح-1C که در 4 دسامبر 2019 به دست آمد و منطقه مورد مطالعه را پوشش میدهد به عنوان دادههای آزمایش سنجش از راه دور نوری انتخاب شده است. شکل 1b تصویر ترکیبی با رنگ کاذب ( FCC ) از Sentinel-2 ( L1C ) را نشان می دهد که محل آزمایش انتخاب شده در منطقه مورد مطالعه را در بر می گیرد. وضوح فضایی بین 10 تا 60 متر و وضوح طیفی بین 15 تا 180 نانومتر برای باندهای مختلف تصویر جمع آوری شده Sentinel-2 متفاوت است ( جدول 1 ).
هنگامی که دادههای Sentinel-2 به وضوح پیکسل 30 متر کاهش مییابند (مثلاً Landsat 8) از باندهای تصویر تقریباً معادل استفاده میکنند، Sentinel-2 از سایر مجموعههای داده فضایی موجود بهتر عمل میکند. علاوه بر این، Sentinel-2 شامل نوار لبه قرمز با کیفیت بالا است. پیشبینی جنگل AGB با استفاده از Sentinel-2 دقت بهتری را در مقایسه با Landsat 8 تأیید کرد [ 45]. بررسی میدانی برای جمعآوری اطلاعات جنگلهای زمینی برای توسعه روش مبتنی بر سنجش از دور و اعتبارسنجی نتایج انجام شد. به منظور پوشش کامل تنوع جنگلی منطقه مورد مطالعه، مکان های حقیقت زمینی بر اساس دسترسی در جنگل و بدون توجه به دسته های ارتفاع درخت به طور تصادفی انتخاب شدند. علاوه بر این، نمونههای زمینی از مکانهای مختلف با چنین الگویی جمعآوری شد تا کل تنوع جنگلی منطقه مورد مطالعه را پوشش دهد.
قطعه های نمونه انتخاب شده نمایشی از کل منطقه مورد مطالعه است. نمونه های زمین از سه پروفیل مختلف که با فاصله نزدیک به 1 کیلومتر از هم جدا شده بودند جمع آوری شد. فاصله بین نمونه ها نیز تقریباً 1 کیلومتر تعیین شد. نمونه های زمینی به صورت تصادفی از این پروفیل ها بر اساس دسترسی به جنگل با جهت های مستقل جمع آوری شد. جمع آوری داده های میدانی با استفاده از یک سیستم موقعیت یاب جهانی ( GPS )، ارتفاع سنج و تجهیزات میدانی جانبی انجام شد. نمونه های زمینی شامل ارتفاع درخت، DBH درخت و تعداد درختان در قطعه نمونه، تراکم ساقه درخت و FVCاز 22 مکان مختلف با اندازه قطعه تقریبی 30 متر × 30 متر (تقریباً 0.09 هکتار)، از سایت مورد مطالعه جمع آوری شد. نمونههای مزرعهای از 22 مکان شامل تقریباً 550 درخت منفرد بود و میانگین 25 درخت در هر قطعه در طبقات مختلف ارتفاعی اندازهگیری شد. میانگین ارتفاع درخت اندازه گیری شده برای هر مکان مزرعه به عنوان ارتفاع درخت محل برای مطالعه بیشتر استفاده شد. در FVC زمینتخمین در هر مکان نمونه، ابتدا تعداد کل درختان و شعاع پوشش هر درخت محاسبه شد. پوشش سطحی تک درختان دایره ای در نظر گرفته شد. بنابراین، سطح پوشش هر درخت با استفاده از شعاع پوشش آنها محاسبه شد. کل سطح پوشش درخت در هر مکان نمونه با جمع کردن سطح پوشش درختان جداگانه برآورد شد. در نهایت، FVC میدانی با در نظر گرفتن نسبت مساحت تحت پوشش درختان به کل مساحت محل نمونه محاسبه شد. در نتیجه، دادههای زمین و FVC برای محاسبه AGB زمینی و اعتبارسنجی AGB مدلسازیشده از دادههای Sentinel-2 استفاده شد. جزئیات بیشتر از داده های زمینی و مکان ها در جدول 2 ارائه شده است. توالی نمونه های زمین ارائه شده در جدول 2 ترتیب جمع آوری داده ها بدون توجه به ترکیب پروفایل است. بنابراین، توالی نمونه ارائه شده در جدول 2 مکان نمونه تصادفی را بدون توجه به الگوی فضایی نمایه نشان می دهد.
2.3. روش شناسی
روش شناسی در دو مرحله ارائه شده است. در مرحله اول، یک مدل کلی توسعهیافته برای ارزیابی جنگل AGB از تصاویر سنجش از دور نوری ارائه شده است، در حالی که فاز دوم آزمایش رویکرد توسعهیافته ویژه تصویر Sentinel-2، شاخصهای پوشش گیاهی و پارامترهای آماری مورد استفاده در این مطالعه را نشان میدهد. .
مرحله I. روش استاندارد برای تخمین جنگل AGB با استفاده از داده های سنجش از راه دور نوری.
روش تخمین AGB جنگل در شش مرحله زیر بیان شد و نمودار جریان روش در شکل 2 نشان داده شده است .
مرحله 1: پیش پردازش: در این مرحله، تصویر سنجش از دور نوری برای کاهش نویز، تصحیح و کالیبراسیون رادیومتری، تصحیح جوی و نمونهبرداری مجدد فضایی از پیش پردازش شد. نویز در اکثر تصاویر ماهواره ای به دلیل نقص سنسور و گرمایش باقی می ماند و تأثیرات محیطی باعث کاهش کیفیت داده ها می شود. از این رو، حذف محتوای نویز در یک تصویر برای افزایش کیفیت داده ضروری است. تصاویر سنجش از دور توزیع شده معمولاً حاوی عدد دیجیتال هستند که یک مقدار تابش کالیبره شده با بهره و افست خاص حسگر است. بنابراین، تصحیح و کالیبراسیون رادیومتریک برای تبدیل اعداد دیجیتال به تشعشع مورد نیاز است.
علاوه بر این، تصحیح اتمسفر برای کاهش اثرات پراکندگی و جذب اتمسفر در دادههای تابشی و تبدیل تابش به دادههای بازتابی انجام میشود. تصویر Sentinel-2 دارای وضوح مکانی متفاوت در باندهای طیفی مختلف است. نمونهگیری مجدد فضایی برای دستیابی به وضوح فضایی واحد از وضوحهای فضایی مختلف در باندهای طیفی مختلف استفاده شد. در پیش پردازش، کاهش نویز و تصحیح رادیومتریک با استفاده از الگوریتمهای مبتنی بر آمار تصویر که با استفاده از بهره و افست حسگر توسعه یافتهاند، انجام شد. یک روش ترکیبی تفریق میدان مسطح و جسم تاریک برای تصحیح جوی و تبدیل یک تصویر رادیومتری به یک تصویر بازتابی توسعه داده شد.
مرحله 2: محاسبه FVC جنگل – در این مرحله، FVC جنگل در هر سطح پیکسل از تصاویر سنجش از دور نوری با استفاده از مدل FVC پیشنهاد شده توسط ژانگ و همکاران محاسبه شد. [ 46 ] به شرح زیر:
فرض کنید اطلاعات طیفی پیکسل پوشش گیاهی خالص و پیکسل خاک خالص به صورت S v و Ss نمایش داده شوند. سپس اطلاعات طیفی ( S ) یک تصویر در هر پیکسل از اطلاعات طیفی پوشش گیاهی ( S veg ) و پوشش غیر گیاهی ( So ) تشکیل شده است . اگر FVC جنگل در یک پیکسل به صورت Fc نمایش داده شود ، آنگاه کسری از پوشش غیر گیاهی باید 1- Fc باشد . بنابراین، از نظر ریاضی، FVC جنگل در سطح پیکسل را می توان به صورت زیر محاسبه کرد:
با استفاده از معادله (4)، FVC از هر دو تصویر سنجش از دور نوری و داده های میدانی محاسبه می شود.
مرحله 3: تخمین AGB جنگل میدانی – در این مرحله، AGB جنگل میدانی از داده های میدانی جمع آوری شده محاسبه شد. جنگل مزرعه AGB در سطح کرت انفرادی بر حسب کیلوگرم (کیلوگرم) با استفاده از مدل تجربی ساده شده زیر پیشنهاد شده توسط چاو و همکاران به دست آمد. [ 27 ] به شرح زیر:
که در آن μ وزن مخصوص (g cm – 3 ) ساقه درخت است، D DBH ( cm)، H ارتفاع درخت (m) و κ ثابت برای نوع جنگل است.
مرحله 4: توسعه مدل بین FVC و AGB — در این مرحله، مدلسازی رگرسیون خطی برای یافتن بهترین منحنی برازش بین AGB جنگل زمینی محاسبهشده از دادههای میدانی با استفاده از مرحله 3 و تصویر جنگل FVC بهدستآمده با استفاده از مرحله 2 انجام میشود.
بهترین مدل خطی مناسب بین AGB جنگل میدان و FVC جنگل تصویر مربوطه به صورت زیر نمایش داده می شود:
که در آن α و β به ترتیب افزایش و افست برای FVC جنگل تصویری هستند که برای محاسبه AGB جنگل استفاده می شود .
مرحله 5: تولید تصویر جنگلی AGB با استفاده از تصویر FVC جنگل – فرض کنید تصویر FVC جنگل به صورت FVC ( i , j ) و تصویر جنگل AGB به صورت AGB ( i , j ) نشان داده شده است، جایی که ( i , j ) محل پیکسل سپس جنگل AGB پیکسلی بر حسب کیلوگرم به صورت زیر محاسبه می شود:
مرحله 6: ارزیابی دقت – دقت AGB جنگل محاسبه شده از تصویر سنجش از دور نوری با اعتبارسنجی زمینی ارزیابی شد. اعتبارسنجی زمین با مقایسه مقدار AGB بین میدان محاسبهشده و تخمین تصویر در مکانهای مختلف زمین انجام شد. دقت نتیجه در درصد به صورت زیر نمایش داده می شود:
که در آن N تعداد کل نقاط اعتبارسنجی است، و n تعداد نقاط اعتبار سنجی توافق شده با آن است که برابر یا کمتر از یک انحراف استاندارد بین اندازه گیری های تصویر و میدان است.
در نهایت، AGB جنگل تعمیمیافته بهدستآمده از سطح کرت فردی بر حسب کیلوگرم به تن در هکتار (tha -1 ) استاندارد میشود.
فاز دوم. ارزیابی مدل جنگلی AGB با استفاده از داده های Sentinel-2.
برای آزمایش و تأیید اثربخشی روش پیشنهادی، آزمایش در روش متوالی بر روی Sentinel-2 بهعنوان یک تصویر سنجش از راه دور نوری برای محاسبه AGB جنگل از سایت تست جنگل رزرو شده توندی، Jharkhand، هند استفاده شد. مراحل متوالی در محاسبه AGB جنگل در زیر آورده شده است.
پس از پیش پردازش، نمونه برداری مجدد فضایی در ابتدا برای رساندن همه باندهای طیفی به وضوح فضایی 10 متر انجام شد (به عنوان مثال، باندهای 2-4 با تفکیک مکانی 10 متر و باندهای 5-7، 8a با وضوح فضایی 20 متر) و سپس نمونه برداری مجدد به وضوح فضایی 30 متر، یعنی تقریباً برابر با وضوح مکانی داده های میدانی است.
استخراج اطلاعات طیفی مستقیم پوشش گیاهی و غیر گیاهی (خاک) برای محاسبه FVC جنگل مشکل است. بنابراین، شاخص های مختلف پوشش گیاهی برای تعیین کمیت اطلاعات پوشش گیاهی در پیکسل یک تصویر استفاده شد. در آزمایش ما، برای تعیین کمیت مقدار پوشش گیاهی در یک پیکسل Sentinel-2، چهار شاخص پوشش گیاهی مورد ارزیابی قرار گرفت: شاخص گیاهی تفاوت نرمال شده ( NDVI )، شاخص گیاهی اصلاح شده ( MVI )، شاخص گیاهی تعدیل شده با خاک ( SAVI )، و خاک اصلاح شده. شاخص گیاهی تنظیم شده ( MSAVI ). NDVI ، MVI ، SAVI ، و MSAVIتصاویر از تصویر Sentinel-2 با استفاده از معادلات زیر تولید شدند:
که در آن L ضریب تنظیم پسزمینه سایهبان است ( مقدار L برابر 0.5 در دادههای بازتابی) و NIR و RED مقادیر بازتابی در باندهای مادون قرمز و قرمز نزدیک تصویر سنجش از راه دور چند باندی هستند.
در آزمایش ما، NIR (باند 8) و RED (باند 4) با دقت از تصویر بازتابی Sentinel-2 برای محاسبه این شاخصهای پوشش گیاهی انتخاب شدند. بیشترین مقدار شاخص پوشش گیاهی در یک تصویر به عنوان پیکسل پوشش گیاهی خالص در نظر گرفته شد و کمترین مقدار به عنوان پیکسل خاک خالص تعیین شد. در نتیجه، پیکسل های خالص برای اطلاعات پوشش گیاهی و غیر گیاهی (خاک) برای محاسبه تصاویر FVC جنگل از شاخص های پوشش گیاهی محاسبه شده با استفاده از رابطه (4) شناسایی شدند. بهترین شاخص پوشش گیاهی در بین شاخص های پوشش گیاهی برای محاسبه AGB جنگل بر اساس ضریب همبستگی ( R ) و میانگین خطای مطلق ( MAE ) انتخاب شد.) بین FVC جنگل زمینی و FVC جنگل تصویری . R درجه تشابه آماری بین یک جفت متغیر است، در حالی که MAE اندازه گیری خطاهای بین مشاهدات زوجی است. R و MAE بین FVC زمین و FVC از Sentinel-2 توسط NDVI ، MVI ، SAVI ، و MSAVI مربوط به مکانهای نمونه زمین محاسبه شد. بر اساس بالاترین همبستگی و کمترین معیار MAE ، تصویر FVC مبتنی بر MSAVI در AGB انتخاب شد.محاسبه علاوه بر این، خط روند با رگرسیون خطی بین زمین و FVC مبتنی بر MSAVI برازش شد که در معادله (13) نشان داده شده است:
در جنگل حفاظت شده تندی، گونه درختی غالب Shorea robusta (سال) است. بر این اساس، بر اساس [ 27 ]، جنگل زمینی AGB با استفاده از معادله زیر از داده های میدانی محاسبه شد ( جدول 2 ):
که در آن AGB جنگل زیست توده بالای زمین (کیلوگرم)، μ (0.667) وزن مخصوص (تراکم ساقه درخت) (g cm – 3 )، D DBH ( cm)، H ارتفاع درخت (m) و ثابت 0.509 به دست آمده از [ 27 ]. AGB جنگل زمینی محاسبه شده از داده های میدانی و تصویر انتخابی جنگل FVC در جدول 5 نشان داده شده است.
یک خط روند بین AGB جنگل زمینی محاسبه شده و FVC جنگل بازیابی شده توسط رگرسیون خطی به دست می آید. مدل زیر (معادله (15)) از 15 نمونه تمرین زمینی (جدول 5) و از 7 نمونه زمین باقی مانده برای اعتبارسنجی نتایج به دست آمده از مدل توسعه یافته استفاده شد (جدول 6).
دقت جنگل AGB بر اساس Sentinel-2 با مقایسه هفت اندازهگیری میدانی باقیمانده با مکانهای مربوطه ارزیابی شد. لطفاً توجه داشته باشید که نمونههایی که برای آموزش انتخاب شدهاند برای تأیید اعتبار انتخاب نشدهاند تا از تعصب در ارزیابی دقت جلوگیری شود. بنابراین، نمونههای آموزشی و اعتبارسنجی قبل از مدلسازی روند از هم جدا شدند. برای پردازش داده ها و اجرای روش پیشنهادی و آزمایش در ارزیابی جنگل AGB ، عمدتاً از نرم افزار برنامه نویسی زبان داده های تعاملی ( IDL ) و محیط برای تجسم تصاویر ( ENVI ) استفاده شد.
3. نتایج و بحث
نتایج ارزیابی جنگل AGB از منطقه مورد مطالعه با استفاده از روش شناسی و تجزیه و تحلیل تجربی شرح داده شده در بخش 2 در زیر ارائه شده است.
تصاویر تخمینی FVC جنگل بر اساس NDVI ، MVI ، SAVI و MSAVI به ترتیب در شکل 3 نشان داده شده است. FVC زمین و تصویر محاسبه شده از داده های مزرعه و Sentinel-2 به ترتیب برای هر شاخص پوشش گیاهی استخراج شده از مکان های مربوطه، در جدول 3 نشان داده شده است. R و MAE بین FVC زمین و FVC برآورد شده با استفاده از شاخص های پوشش گیاهی از Sentinel-2 مورد تجزیه و تحلیل قرار گرفت. مقادیر تخمینی R و MAE در جدول 4 نشان داده شده است. نتایج نشان می دهد که FVC تصویر مبتنی بر MSAVI بیشترین همبستگی و کمترین MAE را ارائه می دهد. بر اساس مدل تجربی AGB [ 27 ]، AGB جنگل زمینی محاسبه شده از داده های میدانی و تصویر انتخابی جنگل FVC در جدول 5 نشان داده شده است.
طبق نتایج حاصل از آزمایش با استفاده از داده های ماهواره ای همراه با داده های اولیه زمینی، FVC جنگل مبتنی بر تصویر از طریق MSAVI محاسبه شد . این بالاترین همبستگی و کمترین میانگین خطای مطلق را با توجه به FVC زمین در مقایسه با سایر شاخص های پوشش گیاهی مورد استفاده در این مطالعه نشان داد ( جدول 4 ؛ شکل 3 ). این مقایسه ها نشان داد که MSAVI بهترین شاخص پوشش گیاهی برای محاسبه FVC جنگل مبتنی بر تصویر بود . علاوه بر این، FVC جنگل مبتنی بر تصویر یک رابطه خطی را با مقدار همبستگی بالای 0.9995 نشان داد که با استفاده از معادله (13) به دست آمده است.
مدل تجربی (معادله (14)) در محاسبه AGB جنگل زمینی مفید بود . علاوه بر این، مدل خطی بین AGB جنگل زمینی و FVC جنگل مبتنی بر تصویر انتخاب شده مقدار همبستگی خوبی 0.75 را نشان داد (معادله (15)). مدل مشتق شده به FVC جنگل مبتنی بر تصویر برای محاسبه AGB مبتنی بر تصویر منطقه مورد مطالعه اعمال شد ( شکل 4 ). این مدل با استفاده از 15 نمونه از 22 نمونه جمع آوری شده از میدان آموزش داده شد. دقت AGB جنگل مبتنی بر Sentinel-2با مقایسه هفت اندازه گیری میدانی باقی مانده با مکان های مربوطه ارزیابی شد. انحراف معیار 3166.50 کیلوگرم در هر قطعه بود که از برآورد جنگل AGB از هر 22 مکان محاسبه شد. نتیجه نشان می دهد که تمام ارزیابی بین AGB جنگل به دست آمده از تصویر و میدان محاسبه شده در یک انحراف استاندارد توافق شده است. ارزیابی فوق بر اساس یک گونه درختی، یعنی Shorea robusta انجام شد. این گونه درختی تقریباً 90 درصد در منطقه مورد مطالعه غالب است.
در این مطالعه، ما دیگر گونههای درختی را نادیده گرفتهایم تا از پیچیدگی روش توسعهیافته جلوگیری کنیم. بنابراین، می توان در نظر گرفت که این مطالعه تقریباً بیش از 90 درصد دقت را ارائه کرده است. عدم قطعیت نتیجه نهایی با استفاده از آزمون مجذور کای مورد آزمایش قرار گرفت. مقدار آزمون کای اسکوئر 0.19 بود، یعنی 19 درصد عدم قطعیت در برآورد.
در نهایت، مدل توسعه یافته برای محاسبه AGB جنگلی منطقه مورد مطالعه با استفاده از تصویر Sentinel-2 ( شکل 4 ) استفاده شد. بر اساس نقشه AGB جنگلی تولید شده، کل AGB جنگل محاسبه شده منطقه مورد مطالعه 7,914,875,916.57 کیلوگرم (7,914,875.92 تن) در 170,704,500 متر مربع ( 17,070.45 هکتار) برآورد شده است. در نتیجه، منطقه مورد مطالعه دارای جنگل AGB به طور متوسط 463.66 تن در هکتار است.
علاوه بر این، کارایی روش با قاضی و همکاران مقایسه شد. مدل [ 47 ] در برآوردهای جنگلی AGB . قاضی و همکاران مدل برآورد جنگل AGB یک رابطه خطی بین ناحیه پیش بینی شده تاج ( CPA ) است که از یک تصویر سنجش از دور نوری و AGB مبتنی بر میدان بدست می آید .
مدل فوق برای Sentinel-2 و AGB زمینی به دست آمده از منطقه مورد مطالعه به دست آمده است که در معادله (16) نشان داده شده است:
به همین ترتیب، اثربخشی روش توسعه یافته با Pandit و همکاران مقایسه شد. مدل [ 9 ]. این مدل بر اساس رگرسیون تصادفی جنگل برای تخمین AGB جنگل توسعه یافته است. این مدلها برای اعتبارسنجی جنگل محاسبهشده AGB اعمال شدند که در جدول 5 و جدول 6 نشان داده شده است .
نتایج نشان می دهد که پنج با قاضی و همکاران. مدل و شش توسط Pandit و همکاران. مدل از میان هفت نقطه اعتبارسنجی که در یک انحراف استاندارد با تخمین های زمینی مطابقت دارند. توجه داشته باشید که انحراف استاندارد جنگل AGB مورد استفاده برای مقایسه از کل داده های زمین محاسبه شد. علاوه بر این، قدرت رویکرد پیشنهادی با استفاده از ضریب همبستگی و میانگین خطای مطلق مورد آزمایش قرار گرفت.
نتایج به دست آمده با این رویکرد، قاضی و همکاران. [ 47 ] مدل، و Pandit و همکاران. مدل [ 9 ] با مقایسه نتایج آنها با برآوردهای AGB جنگل های زمینی مورد ارزیابی قرار گرفت . همبستگی بین نتایج حاصل از تصویر Sentinel-2 و داده های زمینی با رویکرد فوق در مقایسه با قاضی و همکاران (0.9937) بیشتر بود. مدل (0.9796) و Pandit و همکاران. مدل (0.9884). در حالی که میانگین خطای مطلق بین نتایج حاصل از داده های تصویری و زمینی با رویکرد پیشنهادی کمتر (97/582) از قاضی و همکاران بود. مدل (18428.41) و Pandit و همکاران. مدل (1509.64). نتایج نشان داد که رویکرد توسعه یافته در مقایسه با قاضی و همکاران برتر است. [ 47 ] و Pandit و همکاران. [ 9 ] مدل ها. جنگلارزیابی AGB توسط قاضی و همکاران. مدل توسط CPA مبتنی بر NDVI انجام شد در حالی که تخمین AGB با رویکرد پیشنهادی توسط FVC مبتنی بر MSAVI به دست آمد . به طور مشابه، یک مدل رگرسیون تصادفی مبتنی بر جنگل برای ارزیابی AGB جنگل توسط Pandit و همکاران توسعه داده شد. [ 9 ]. این مدل با استفاده از الگوریتم جنگل تصادفی رابطه ای بین شاخص پوشش گیاهی مشتق شده از لبه قرمز و AGB زمین ایجاد کرد [ 9 ].
این نتایج همچنین نشان داد که انتخاب شاخصهای پوشش گیاهی نیز ممکن است پارامترهای مهمی در تخمین AGB جنگل با دقت باشد. ممکن است دلیلی وجود داشته باشد که MSAVI در مقایسه با NDVI برای تشخیص پوشش گیاهی در جنگل دقیق تر است. علاوه بر این، تخمین دقیق CPA در مقایسه با تخمین FVC از دادههای سنجش از راه دور نوری مشکل است. علاوه بر این، Pandit و همکاران. مدل نشان داد که نه تنها انتخاب الگوریتم قوی ارزشمند است، بلکه شناسایی و ارزیابی دقیق پوشش گیاهی نیز در جنگل AGB حیاتی است.توسعه مدل با این حال، مدل توسعهیافته انتخاب شاخص پوشش گیاهی را برای تعیین کمیت پوشش جنگلی قبل از توسعه مدل جنگلی AGB در نظر گرفت. رویکرد توسعهیافته مستقیم و ساده برای ارزیابی دقیق جنگل AGB است. علاوه بر این، یک رویکرد AGB جنگل تعمیمیافته ارائه شده است به طوری که رویکرد توسعهیافته میتواند برای سایر دادههای سنجش از راه دور نوری برای تخمین جنگل AGB از مکانهای مطالعه مختلف اعمال شود. به طور کلی، این نتایج و تجزیه و تحلیل از منطقه مورد مطالعه نشان می دهد که Sentinel-2 دارای پتانسیل قابل توجهی برای تخمین AGB جنگل است. بر اساس رویکرد توسعهیافته و تصویر Sentinel-2، میانگین جنگل AGB66/463 تن در هکتار از منطقه مورد مطالعه به دست آمد که به طور موثر با برآوردهای مبتنی بر مشاهدات زمینی مطابقت داشت.
یکی از اهداف کلیدی این کار، توسعه روشی جدید برای ارزیابی دقیق، دقیق و سریع زیست توده جنگلی روی زمین از دادههای سنجش از دور نوری ماهوارهای بود. برای این منظور از داده های ماهواره Sentinel-2 استفاده شد. FVC زمین از تصاویر سنجش از دور نوری در هر مکان نمونه زمین استخراج شد. جدای از شاخصهای مختلف پوشش گیاهی موجود، روش جدیدی برای ارزیابی دقیق، دقیق و سریع و اعتبار سنجی تخمین زیست توده جنگلهای بالای زمین با استفاده از سنجش از دور نوری در منطقه جنگلی حفاظتشده توندی، هند توسعه داده شده است. این منطقه جنگلی به طور معمول نشان دهنده یک نوع جنگل نیمه گرمسیری است که به وفور در منطقه ای که مطالعه برای آن انجام شده است یافت می شود. اعتبار سنجی زمینی و مقایسه با CPA توسعه یافتهتایید شدند. این روش نتایج بهتری را نسبت به هر روش شناخته شده دیگری که تاکنون ارائه شده است نشان داد. با استفاده از این روش، AGB 463 تن در هکتار برای منطقه مورد مطالعه به دست آمد و در ادامه، داده های به دست آمده از ماهواره نیز با داده های تخمینی زمین مقایسه شد. در نهایت، با استفاده از مدل توسعه یافته، سطح دقت مطلوبی به دست آمد. این منجر به درک و درک ارزش نتیجه تحقیق از این کار می شود.
این مطالعه همچنین نشان داد که متغیرهای بیوفیزیکی انواع پوشش گیاهی بهدستآمده از Sentinel-2 در محاسبه AGB جنگل با استفاده از رگرسیون محلی مفید هستند. با این حال، این مطالعه فرض کرد که کل منطقه مورد مطالعه از یک گونه درختی تشکیل شده است، در حالی که در واقعیت جنگل ها از گونه های درختی متعدد با مقادیر جغرافیایی متفاوت در شرایط زمین مواج تشکیل شده اند. با این حال، در مقایسه با مجموعه دادههای موجود و تکنیکهای موجود، روش پیشنهادی در ارزیابی سریع AGB جنگل بسیار مفید است. علاوه بر این، در دسترس بودن داده های سنجش از دور چندطیفی نوری با کیفیت بالا در مقیاس جهانی، رایگان یا با هزینه کم از منابع متعدد، ممکن است به جامعه جهانی برای جنگل AGB کمک کند.برآورد کردن.
در نهایت، این نتایج منجر به این نتیجه می شود که مدل به دست آمده به سرعت AGB جنگل را از تصاویر سنجش از دور نوری به طور مستقیم با دقت بالا تولید می کند. استفاده از دادههای سنجش از دور نوری در مقایسه با دادههای سنجش از دور مایکروویو، که برای استخراج پارامترهای بیوفیزیکی گیاهی در برآورد AGB به پردازش دادههای پیچیده نیاز دارند، سریع و ساده است. علاوه بر این، روشهای قبلی [ 9 ، 47 ] و دادههای مایکروویو در محاسبات جنگلی AGB به عوامل متعددی مانند عوامل اکولوژیکی، زیستی و توپوگرافی، از جمله ارتفاع واقعی درخت، قطر، تراکم ساقه، و شکل اطلاعات تاج درخت بستگی دارد.
4. نتیجه گیری
هدف از این تحقیق توسعه یک روش قوی و مقرون به صرفه برای ارزیابی سریع AGB جنگل بود. در ابتدا، FVC زمینی با استفاده از کسری از پوشش زمین توسط درختان و تصویر FVC با استفاده از شاخص پوشش گیاهی از دادههای سنجش از دور نوری محاسبه شد. در مرحله بعد، جنگل زمینی AGB از داده های زمین (تراکم درخت، ارتفاع درخت و DBH درخت ) با استفاده از یک مدل تجربی [ 27 ] استخراج شد. متعاقباً، یک مدل رابطهای بین جنگلهای زمینی AGB و سنجش از دور FVC با ترکیب این اطلاعات با استفاده از رویکرد برازش منحنی خطی ایجاد شد.
برای تخمین AGB جنگل ، مدل ابداع شده بر روی تصویر FVC به دست آمده از تصویر سنجش از دور نوری اعمال شد. اثربخشی رویکرد توسعهیافته با استفاده از تصاویر Sentinel-2 بهعنوان دادههای سنجش از راه دور نوری آزمایشی ارزیابی شد. منطقه جنگلی حفاظت شده توندی، واقع در ایالت جارکند در شرق هند، به عنوان مکان آزمایشی برای این مطالعه انتخاب شد. نتایج بهدستآمده با رویکرد توسعهیافته تطابق بسیار خوبی با نتایج بهدستآمده توسط برآوردهای زمینی نشان داد. استفاده از روش توسعهیافته در جنگل حفاظتشده توندی نشان داد که محل آزمایش شامل جنگل AGB استبه طور متوسط 463 تن در هکتار. نتیجه به دست آمده توسط رویکرد پیشنهادی به خوبی با نتایج گزارش شده در مطالعات قبلی از منطقه مورد مطالعه مطابقت دارد. نتیجه بیش از 90% تطابق با اندازه گیری زمین را نشان داد.
این مطالعه هنوز دارای یک دامنه بهبود آینده است که از اطلاعات توپوگرافی در مدل توسعهیافته استفاده میکند.
بدون دیدگاه