پیوندهای مفید
اخبار و اعلانات
آمار
| تعداد نشریات | 43 |
| تعداد شمارهها | 1,686 |
| تعداد مقالات | 13,791 |
| تعداد مشاهده مقاله | 32,417,891 |
| تعداد دریافت فایل اصل مقاله | 12,801,028 |
پایش و برآورد وسعت مناطق گرفتار حریق بخشی از کوهستان زاگرس با استفاده از تصاویر ماهوارهای لندست | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| جغرافیا و برنامه ریزی محیطی | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| مقاله 4، دوره 33، شماره 4 - شماره پیاپی 88، دی 1401، صفحه 49-62 اصل مقاله (1023.45 K) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| نوع مقاله: مقاله پژوهشی | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| شناسه دیجیتال (DOI): 10.22108/gep.2022.131560.1470 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| نویسندگان | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| موسی عابدینی* 1؛ مریم محمدزاده شیشه گران2؛ احسان قلعه3 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 1استاد گروه جغرافیای طبیعی دانشگاه محقق اردبیلی، اردبیل، ایران | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 2دانشجوی دکتری ژئومورفولوژی گرایش ژئومورفولوژی دانشگاه محقق اردبیلی، اردبیل، ایران | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 3دانشجوی دکتری ژئومورفولوژی دانشگاه محقق اردبیلی، اردبیل، ایران | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| چکیده | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| جنگلها ازنظر تأمین تعادل اجتماعی و زیستمحیطی جایگاه فوقالعاده مهمی در اکوسیستم دارد. بزرگترین خطر برای چنین جنگلهایی، آتشسوزی است. هرساله سطوح زیادی از جنگلهای منطقۀ زاگرس در غرب ایران، دستخوش آتشسوزی میشود و از بین میرود؛ بنابراین برآورد شکلگیری و ویژگیهای رفتاری آتش برای مقابله با آن بسیار مهم است. بهرهگیری از دادهها و تصاویر ماهوارهای مناسب و در دسترس، اطلاعات مفیدی را دربارۀ شرایط پیش و پس از آتشسوزی عرصههای جنگلی فراهم میکند. این پژوهش، با هدف برآورد مساحت مناطق گرفتار حریق و آگاهی از کارایی، قابلیت دادههای ماهوارۀ لندست و شاخصهای NBR و dNBR در تشخیص، ارزیابی و تهیۀ نقشه جنگلهای سوختهشدۀ منطقۀ زاگرس انجام شد. برای این منظور، پس از تهیۀ تصاویر لازم، آتشسوزیهای اتفاقافتاده در خردادماه سال 1399 در جنگلهای زاگرس با بهرهگیری از تکنیکهای سنجش از دور و انجام پردازشهای لازم بررسی شد. نتایج بهدستآمده نشاندهندۀ آن بود که شاخصهای NBR و dNBR اطلاعات مناسبی را درخصوص تأثیر آتشسوزی و روند تغییرات آن در اختیار قرار میدهد؛ همچنین 13685 هکتار از جنگلهای زاگرس در این آتشسوزی طعمۀ حریق شدهاند. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| کلیدواژهها | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| آتشسوزی؛ شاخصهای NBR و dNBR؛ طبقهبندی؛ زاگرس | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| اصل مقاله | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
مقدمه به دلیل قرارگرفتن کشور ایران در کمربند خشک کرۀ زمین و ناحیۀ پرفشار جنب حارهای، شرایط جوی برای وقوع حوادث غیرمترقبه بهطور کامل فراهم است. مطابق بررسیهای انجامشده، کشور ما جزو 10 کشور حادثهخیز دنیاست. یکی از این حوادث که به وفور در ایران اتفاق میافتد، پدیدۀ آتشسوزی جنگلها و مراتع است. جنگلها، بخش مهمی از سیستم زیستبوم زمین و منبع بزرگی برای مقاصد مختلف ازجمله مخزن ژنتیک، مخزن آب، منبع کربن و یک سرچشمۀ ذخیرۀ انرژی در طبیعت است. جنگلها در بهترکردن محیطزیست و در تعادل نگهداشتن آن نقش اساسی دارد و همزمان برای ارتباط با توسعۀ مناسب در اقتصاد اجتماعی یک منبع مهم طبیعی است (سرکارگر اردکانی و همکاران، 1388: 26)؛ اما این منابع عظیم با آتشسوزیها به خطر افتاده است. هرساله در مناطق مختلف، هزاران هکتار از جنگلها طعمۀ حریق میشود. آتشسوزی جنگل با منشأ طبیعی یا انسانی، اثرهای زیانبار و ویرانگری به شکل مستقیم یا غیرمستقیم بر زندگی بشر به جای میگذارد. علاوه بر تخریب محیطزیست و آلودگی آن، موجب نابودی ذخایر چوب، احشام، زمینهای کشاورزی و چرایی، ساختمانها و جان و مال انسانها میشود و پیامدهای اقتصادی، اجتماعی و روانی زیادی در پی دارد (Niu & Zhai, 2012: 724). براساس گزارش سازمان جنگلها، مراتع و آبخیزداری کشور، عرصههای منابع طبیعی ایران شامل جنگل، مرتع، بیابان، بیشهزار و درختچهزارها، معادل 48/83 درصد از مساحت کشور است که از این بین، سهم جنگلها و مراتع 82/61 درصد است و عمدۀ آتشسوزی ایران در این عرصهها اتفاق میافتد. مهمترین عواملی که زمینۀ آتشسوزی را در جنگلها و مراتع فراهم میکند، کاهش رطوبت نسبی هوا و در دسترس بودن مواد سوختی است. مادۀ سوختی، از خشکشدن شاخ و برگ گیاهان به دلیل خشکسالی یا یخبندان در فصل زمستان فراهم میشود؛ همچنین در فصل پاییز که سطح جنگل را برگهای خشکشده میپوشاند، گرمشدن هوا زمینۀ مناسبی را برای آتشسوزی مواد سوختی فراهم میکند و انباشت مواد سوختی در جنگل، خطر آتشسوزی را افزایش میدهد. نگهداری درختان خشکشده در فصولی که خطر آتشسوزی وجود دارد نیز بر ریسک حریق میافزاید (Niu & Zhai, 2012: 725). از طرف دیگر، بیاحتیاطی ساکنان یا مسافران و نبودِ امکان دستیابی به مناطق آتشگرفته، باد، توپوگرافی و شکل زمین و مدیرت آتش از عوامل اثرگذار در شدت و گسترش آتشسوزی است (محمدی، 2008: 31). آتش، یک پدیدۀ جهانی است که از مؤثرترین فاکتورها در اکوسیستمهای خاکی کرۀ زمین به شمار میرود. به دلیل قرارگرفتن کشور ایران در کمربند خشک کرۀ زمین و ناحیۀ پرفشار جنب حارهای، شرایط جوی لازم برای وقوع آتشسوزی در جنگلها و مراتع فراهم است و گاهی اوقات خسارات زیادی را به جنگلها و مراتع وارد میکند. در راستای نگهداری از جنگلها و حفظ منابع طبیعی، راهبردها و راهکارهای متفاوتی در دنیا استفاده میشود (Kaufman et al, 2003: 1767). برای مبارزه با آتشسوزی جنگل باید همزمان رفتار آتش و نیز نحوۀ مبارزه با آن را در یک دنیای مجازی، قبل، حین و بعد از وقوع آتشسوزی بازسازی و رفتار گسترش آتش، استراتژیها و روشهای مبارزه با آن را بررسی کرد. سنجندههای ماهوارههای سنجش از دور با قدرت تفکیک مکانی و زمانی مناسب، اطلاعات مهمی را برای هشدارهای زودهنگام آتشسوزی در دسترس قرار میدهد؛ همچنین فرصتهایی را برای آنالیز کیفی جنگل و دیگر اکوسیستمها در مقیاسهای مکانی، جغرافیایی، تکنولوژی و ژئومکانی فراهم میکند و اطلاعات و ابزارهای ضروری را برای ایجاد یک نقشۀ حساسیت آتشسوزی از طریق شناسایی، کلاسهبندی و تهیۀ نقشۀ خطر آتشسوزی ایجاد میکند (Labib & Harris 2018: 234). با تکیه بر پتانسیلهای یادشده از دادههای سنجش از دور در زمینۀ تهیۀ نقشۀ خطر آتشسوزی جنگل و بررسی عوامل مؤثر در ایجاد و گسترش آتشسوزی، پژوهشهایی در داخل و خارج از کشور انجام شده است. Quintano et al (2018) شدت سوختگی را با استفاده از سه شاخص dNBR، RdNBR و نسبت سوختگی نسبیشده (RBR) را با استفاده از تصاویر ماهوارهای سنتینل 2 و لندست 8 بررسی کردند. نتایج نشاندهندۀ آن بود که شاخص RdNBR مبتنی بر ماهوارۀ Landsat به بیشترین دقت رسیده است. Garcia-Llamas et al (2019) آتشسوزی جنگلی را بررسی کردند که در منطقۀ کوهستانی کابرر (Cabrere) در شمال غربی اسپانیا در آگوست 2017 با استفاده از دادههای ماهوارهای Landsat-8، Sentinel-2 و Deimos-1 رخ داد. آنها شاخصهای انعکاسی، حرارتی و ترکیبی را روی تصاویر اعمال کردند و به این نتیجه رسیدند که Sentinel-2 نشاندهندۀ نتایج بهتری است و شاخصهای مبتنی بر NBR همبستگی بیشتری را برای تجزیهوتحلیل شدت سوختگی، پوشش گیاهی و خاک ارائه میکند. Chen et al (2020) در پژوهشی، شدت آتشسوزی را با اعمال ده شاخص بر روی دادههای Landsat شناسایی کردند و نتیجه گرفتند که شاخص dGEMI و dNBR قویترین رویکرد را در منطقۀ موردمطالعه ارائه میدهد. Kurnaz et al (2019) در پژوهشی، هفت شاخص dNBR، RBR، RdNBR، dNDVI، dGNDVI، BAI و GEMI را برای شناسایی مناطق گرفتار حریق استفاده کردند که از بین شاخصهای مورداستفاده، شاخص dNBR و RdNBR بیشترین دقت و عملکرد را در شناسایی مناطق گرفتار حریق داشت. هدف از پژوهش حاضر، کشف آتش و تعیین وسعت مناطق گرفتار حریق در آتشسوزی اتفاقافتاده در خردادماه سال 1399 در کوهستان زاگرس شامل بخشهایی از استانهای کرمانشاه، ایلام و کردستان با استفاده از تصاویر ماهوارهای لندست 8 است.
روششناسی پژوهش پهنۀ رویشی زاگرس شامل رشتهکوههای زاگرس، وسیعترین و اصلیترین رویشگاه گونههای مختلف بلوط در ایران بوده است و به همین دلیل این منطقه اهمیت بسیار ویژهای دارد. جنگلهای بلوط منطقۀ زاگرس از منتهیالیه شمال غربی ایران استان آذربایجان غربی آغاز شده است و تا غرب ایران (استانهای کردستان، کرمانشاه، ایلام، چهارمحال و بختیاری، خوزستان و لرستان و جنوب غرب ایران استانهای کهگیلویه، بویر احمد و فارس) ادامه دارد. جنگلهای زاگرس با گستردگی در ۱۱ استان کشور با ۶ میلیون هکتار مساحت، ۴۰ درصد جنگلهای ایران را تشکیل میدهد که حدود ۷۰ درصد تیپ گونههای جنگلی زاگرس را بلوطها شامل میشود. پوشش گیاهی در نواحی غربی شامل گیاهان سردسیری و در نواحی شرقی شامل گیاهان گرمسیری است. بخش عمدهی پوشش گیاهی زیستگاه جنگلی زاگرس را درختان بلوط تشکیل میدهد. در ارتفاعات بالاتر از مرز جنگل (حدود 2300 متر) علفزارها و بوتهزارهای انبوه مشاهده میشود. بارندگی این ناحیه از ۶۰۰ میلیمتر در قسمتهای شمالی تا ۳۰۰ میلیمتر در جنوب غرب در نوسان است. اقلیم منطقه نیز مدیترانهای نیمه خشک با زمستانهای سرد است. فرم غالب جنگلشناسی در زاگرس، جنگلهای شاخه زاد بوده است و این دلالت بر تخریب ناشی از حضور طولانی مدت جوامع انسانی دارد. زاگرس شامل زاگرس شمالی، میانی و جنوبی است که تفاوت آنها بیشتر در قسمت شمالی- جنوبی بهلحاظ تفاوت در چینخوردگیها مشهود است. زاگرس مرکزی اغلب کوهستانی و در دامنۀ ارتفاعی ۸۳۰ تا ۴۴۱۶ متر از سطح دریا قرار گرفته است. میانگین ارتفاع محدوده ۲۳۳۲ متر است. ۳۴ درصد از محدوده، ارتفاع ۲۴۰۰– ۲۱۰۰ متر و کمتر از یک درصد آن، ارتفاع کمتر از ۱۰۰۰ متر دارد. متأسفانه به دلیل رشد جمعیت در قرن اخیر و نیاز جنگلنشینان به سوخت، منبع غذایی برای دامها و زمینهای کشاورزی و بهرهبرداری غیراصولی از این جنگلها، بهکلی چهرۀ آنها را دگرگون کرده است و بهجز در نقاط معدودی، نشانی از انبوهی سابق این جنگلها نیست.
شکل (1) موقعیت جغرافیایی منطقۀ موردمطالعه(منبع: نویسندگان، 1400) Figure (1) Geographical location of the studied area (source: authors, 1400)
در این پژوهش، از تصاویر ماهوارۀ لندست 8 استفاده و این تصاویر از سایت USGS (سازمان زمینشناسی ایالات متحده) اخذ شد. لندست 8، هشتمین ماهواره از این مجموعه است و اطمینان حاصل میشود که مهمترین نقش برنامۀ لندست در تعیین و نظارت بر منابع لازم برای معیشت انسان مانند غذا، آب و جنگلها ادامه خواهد داشت. از تصویر ماهوارۀ لندست 8، سنجندۀ OLI برای استخراج نقشۀ کاربری اراضی و از سنجندۀ TIRS برای استخراج دمای سطح زمین و مناطق دچار حریق استفاده شد. جدول (1) اطلاعات تصاویر ماهوارۀ لندست 5 و 8. (منبع: متادیتا تصاویر) Table (1) information of Landsat 5 and 8 satellite images. (Source: image metadata)
در این پژوهش، آتشسوزیهای رخداده در استانهای کرمانشاه، ایلام و کردستان از آتشسوزیهای جنگلها و مراتع زاگرس در خردادماه 1399 انتخاب شد. آتشسوزی در استان کرمانشاه از روز پنجشنبه 8 خردادماه 1399 آغاز شد و به علت اینکه به موقع مهار نشد، در 9 خردادماه 1399 به محدودۀ استانهای ایلام و کردستان سرایت کرد و گسترش یافت. دادههای استفادهشده با توجه به گستردگی آتشسوزیها در طول این دورۀ زمانی و محدودیت دسترسی به تصاویر ماهوارهای بهصورت موردی همراه با اطلاعات توصیفی است. ابتدا نواحی دارای آتشسوزی شناسایی و سپس بهصورت موردی تصاویر ماهوارهای لندست 8 برای پیش و پس از این دو آتشسوزی تهیه شد. به دلیل شرایط متفاوت اخذ تصاویر سنجش از دور ازجمله حضور اتمسفر در فاصلۀ بین سنجنده و زمین و ...، خطاهایی در تصویر به وجود میآید؛ به همین دلیل، قبل از استفاده از این تصاویر، برخی پیشپردازشها (Pre-Processing) ازجمله تصحیح اتمسفری، تصحیح رادیومتریکی و توپوگرافی برای بهبود کیفیت تصاویر انجام شد
محاسبۀ دمای سطح زمین بسیاری از روشهای بازیابی دمای سطح زمین از پژوهشهای سنجش از دور ارائهشده از سوی چند کانال مادونقرمز دادههای حرارتی استفاده میکند و این روش بهطور تقریبی به سه گروه تقسیم میشود: الگوریتم تک کانال، الگوریتم چند کانال به عنوان مثال، الگوریتم پنجره مجزا و روش چند زمانه. پژوهشها نشاندهندۀ آن است که نقش سنجش از دور حرارتی در مطالعه و برآورد دمای سطح زمین بسیار حائز اهمیت است. دادههای سنجش از دور حرارتی، امکان بررسی و تخمین دمای سطح زمین و همچنین جزایر حرارتی را فراهم میکند. از بین دادههای حرارتی سنجش از دور، تصاویر ماهوارهای جایگاه خاصی دارد و تصاویری همچون لندست، استر، مودیس و ... در برآورد دمای سطح زمین به دلیل داشتن باندهای حرارتی مفید واقع میشود. دمای سطح زمین بهعنوان نمایهای از شدت گرما، از عناصر اساسی شناخت آبوهواست (درویشی و همکاران، 1398: 20). در سنجش از دور، دمای سطح زمین که بالای جو با سنجندهها به دست میآید، دمای روشنایی نامیده میشود. این دما به علت پدیدههای جذب و پخش جوی به نظر کمتر از دمای واقعی در سطح زمین است. برای استخراج دمای سطح زمین، روشهای مختلفی توسعه داده شده است که هرکدام از این روشها مزایا و معایب خاص خود را دارد. روشهای مختلف استخراج دمای سطح زمین عبارت است از: روشهای تک کانالی، توازن انرژی برای سطح زمین (سبال)، روش دفتر علوم لندست، روش استفان بولتزمن و روش پنجرۀ مجزا (کاکه ممی و همکاران، 1399: 30). در این پژوهش، از روش پنجرۀ مجزا استفاده شد که مراحل آن بهصورت زیر است. برآورد دمای روشنایی دمای روشنایی، تشعشع الکترومغناطیسی انتقالی رو به بالا از لایۀ بالای اتمسفر زمین است که با استفاده از ثابتهای حرارتی ارائهشده در فایل مرجع تصویر رادیانس طیفی به دمای روشنایی تبدیل شد. فرایند تبدیل حرارتی با تبدیل مقادیر ارزش رقومی باندهای حرارتی خام باندهای حرارتی سنجندۀ TIRS به رادیانس طیفی لایۀ بالای اتمسفر انجام شد و سپس با استفاده از رابطۀ (1) (رابطۀ پلانک) دمای روشنایی محاسبه شد. درواقع در این مرحله، واریانس طیفی به دمای جسم سیاه در سطح سنجنده تبدیل شد (Aboelnour & Engel, 2018: 60).
که: T= دمای روشنایی ماهواره برحسب کلوین؛ K2- K1= ثابت حرارتی برای هر دو باند حرارتی موجود درمتادیتا؛ Lλ= درخشندگی جو بالاست.
برآورد قابلیت انتشار سطح زمین به دلیل تفاوت قابلیت انتشار یک مادۀ واقعی در مقایسه با جسم سیاه، بسته به نوع مادۀ تشکیلدهندۀ پدیدههای مختلف، در طول موجهای مختلف ناحیۀ مادونقرمز حرارتی با یکدیگر، پدیدههای مختلف با وجود دمای یکسان، با استفاده از این پارامتر تشخیص داده میشود. هنگامی که قابلیت انتشار انرژی از اجسام واقعی در طول موجهای مختلف اندازهگیری شود، مقادیر متفاوتی حاصل خواهد شد که این امر وابسته به طول موج است (ویسی و همکاران، 1395: 31). برای محاسبۀ دمای سطح زمین برآورد قابلیت انتشار در سطح زمین ضروری است که این معیار بهصورت جداگانه در باندهای حرارتی سنجندۀ TIRS ماهوارۀ لندست 8 (باند 10 و 11) محاسبه شود که رابطۀ آن بهصورت زیر است (Gondwe et al, 2018: 20). رابطۀ (2) که: LSE= قابلیت انتشار سطح زمین؛ Es-v= انتشار گیاه و خاک برای باندهای حرارتی طبق جدول (2)؛ FVC= نسبت پوشش گیاهی است. جدول (2) قابلیت انتشار برای گیاه و خاک Table (2) Emissivity for plant and soil
محاسبۀ نهایی دمای سطح زمین (LST) با روش پنجرۀ مجزا این الگوریتم ازجمله روشهای مؤثر در برآورد دمای سطح زمین است و نسبتبه سایر روشها صحت بیشتری دارد (Sobrino et al, 2004: 436). ویژگی مهم این روش، حذف اثرهای جوی است. با توجه به اینکه این الگوریتم، اطلاعات دقیقی دربارۀ نمایۀ جوی کسب ماهوارهای نیاز ندارد، برای بازیابی LST از دو سنسور چند طیفی و حرارتی استفاده میکند که از رابطۀ زیر به دست میآید (Wan & Dozier, 1996, p. 895).
که: TB10-TB11= دمای روشنایی؛C0…C6 = ضرایب الگوریتم دیدنی در جدول (3)؛ W= مقدار بخار آب موجود در جو (رابطۀ 4)؛ Δε= اختلاف قابلیت انتشار سطح زمین؛ ε= میانگین قابلیت انتشار سطح زمین است. رابطۀ (4) W= C0 + C1 * (Tj / Ti) + C2* (Tj/ Ti)2 که: Tj / Ti = میانگین دمای روشنایی؛ C0 = 674/9-؛ C1 = 653/0؛ C2 = 087/9 است.
جدول (3) مقادیر ضرایب Table (3) values of coefficients
شاخص نسبت سوختگی نرمالشده ((Normalized Burn Ratio شاخص نسبت سوختگی نرمالشده (NBR) یکی از شاخصهای پرکاربرد در مطالعات مشاهدۀ آتش است که با مطالعات سنجش از دور انجام میشود. شاخص NBR، برای برجستهکردن مناطق سوزاندهشده در مناطق بزرگ آتشسوزی طراحی شده است. مقادیر شاخص NBR، از مادونقرمز نزدیک (NIR) و باندهای مادونقرمز موج کوتاه (SWIR) با تناسب تفاوتها و مجموع آنها اعمال میشود. شاخص NBR، در مقادیر بین 1- و 1+ بیان میشود. از آنجایی که کاهش پوشش گیاهی در فرایند پس از آتشسوزی رخ میدهد، بازتاب نوار مادونقرمز نزدیک بهشدت کاهش مییابد. پوشش گیاهی سالم بهطور کلی ارزش بالایی دارد؛ درحالیکه زمین خالی یا مناطق تازه سوخته، ارزش شاخص پایینی دارند. معادلۀ شاخص نسبت سوختگی نرمالشده در رابطۀ 5 آورده شده است (Novelli et al, 2016: 405). رابطۀ (5)
شاخص اختلاف نسبت سوختگی نرمالشده (Differenced Normalized Burn Ratio) شاخص اختلاف نسبت سوختگی نرمالشده (dNBR)، اختلاف نواحی سوخته را با استفاده از تفاوت NBR محاسبهشده از تصاویر قبل و بعد از رویداد تعیین میکند. از آنجایی که مقادیر خاکستری در شاخص NBR قبل از آتشسوزی بزرگتر از مقادیر خاکستری در مناطق سوخته پس از آتشسوزی است، نتیجه گرفته میشود که آتشسوزی در مناطق با مقادیر خاکستری مثبت رخ داده است؛ درحالیکه مناطق با مقادیر منفی بر آتشسوزی اثر نمیگذارد. (Parks et al, 2014: 1831). رابطۀ (6) طبقهبندی نظارتشده یکی از روشهای قدرتمند در ارتباط با تولید نقشههای موضوعی، استفاده از روشهای طبقهبندی پیکسل مبنا از طریق چند باند طیفی متفاوت است. بهصورت کلی، ساختار الگوریتم طبقهبندی پیکسل مبنا در دو مرحله خلاصه میشود: خوشههای طیفی ایجادشده در فضای ویژگی و اختصاص داده شدن هر یک از پیکسلها به یک کلاس (Hussaina et al, 2013: 94). طبقهبندی نظارتشده بهطور معمول بهدنبال وجه تشابهی با پیکسلهای معلوم هر کلاس است. این پیکسلهای معلوم را اشیای نمونۀ تعلیمی مینامند. برای انجام طبقهبندی نظارتشده، انتخاب نمونههای تعلیمی لازم از کلاسهای مختلف با شناخت از منطقه، استفاده از اطلاعات جانبی و بر مبنای DN پیکسلهای تصویر رنگی تصاویر ماهوارۀ لندست صورت گرفت. روش بیشترین احتمال (Maximum Likelihood) از میان روشهای طبقهبندی نظارتشده، تاکنون بهعنوان دقیقترین و پراستفادهترین روشها شناخته شده است. در روش مذکور، فرض بر این است که همۀ مناطق آموزشی دارای پراکنش نرمال است. درحقیقت نمونههای کلاسهای آموزشی باید معرف آن کلاس باشد (مختاری و نجفی، 1394: 38).
شکل (2) فلوچارت مراحل تحقیق (منبع: نویسندگان، 1400) Figure (2) flowchart of the research steps (source: authors, 1400) یافتههای پژوهش و تجزیهوتحلیل آنها نقشههای خروجی دمای سطح زمین با استفاده از تصویر ماهوارۀ لندست 8 سنجندۀ TIRS باند 10 در شکل (3) نشان داده شده است. با توجه به نقشههای استخراجشده، دمای سطح زمین برای تصویر قبل از آتشسوزی کمترین و بیشترین دما بهترتیب 6 و 32 درجۀ سانتیگراد و برای تصویر بعد از آتشسوزی بهترتیب 6 و 41 درجۀ سانتیگراد است که افزایش دمای 9 درجۀ سانتیگرادی در بازۀ زمانی قبل و بعد از آتشسوزی حاکی از وجود آتشسوزی در منطقۀ موردمطالعه است.
شکل (3) نقشۀ دمای سطح زمین قبل از آتشسوزی (سمت راست) و بعد از آتشسوزی (سمت چپ) Figure (3) Earth surface temperature map before the fire (right side) and after the fire (left side)
پس از انتخاب الگوریتمهای شناسایی آتشسوزی، به پیشپردازش (تصحیحات هندسی، رادیومتریک و اتمسفری)، پردازش (اعمال روابط و فرمولهای آتشسوزی برای شناسایی آتشسوزی) و پسپردازش (ارزیابی و صحتسنجی نتایج آتشسوزیهای شناساییشده با الگوریتمهای آتشسوزی) تصاویر اقدام شد. نتایج نهایی آتشسوزیهای شناساییشده برای بخشهایی از رشتهکوه زاگرس از سوی الگوریتمهای آتشسوزی نشاندهندۀ آن است که آتشسوزیهای شناساییشده با شاخص NBR طبق شکل (4) در منطقۀ موردمطالعه بهصورت پراکنده و با شدت زیادی رخ داده است.
شکل (4) نقشۀ شاخص NBR قبل از آتشسوزی (سمت راست) و بعد از آتشسوزی (سمت چپ) Figure (4) NBR index map before the fire (right side) and after the fire (left side) از تفاوت بین NBR قبل از آتشسوزی و بعد از آتشسوزی تصاویر برای محاسبۀ dNBR استفاده میشود که برای تخمین شدت سوختگی مورداستفاده قرار میگیرد. مقادیر dNBR از حالتی به حالت دیگر متفاوت است؛ بنابراین اگر تفسیر در موارد خاص برای دستیابی به بهترین نتایج باشد، باید از طریق ارزیابی میدانی انجام شود؛ با این حال، سازمان زمینشناسی آمریکا یک طبقهبندی را برای تفسیر شدت سوختگی پیشنهاد کرد؛ درنتیجه برای برآورد شدت سوختگی در شاخص dNBR شش درجۀ مختلف شدت آتش تعریف میشود و به هر یک رنگی اختصاص داده میشود (جدول 4). جدول (4) طبقهبندی شدت سوختگی شاخص dNBR Table (4) classification of dNBR index burn severity
با توجه به نقشۀ dNBR حاصل از تصاویر ماهوارهای (شکل 5)، مشخص شد که آتشسوزی اتفاقافتاده در منطقۀ موردمطالعه بهصورت پراکنده و با شدت زیاد از شمال غرب آغاز شده و بهسمت جنوب شرق ادامه پیدا کرده است که با شدت و ضعف مختلف، بر منطقۀ وسیعی اثر گذاشته است.
شکل (5) نقشۀ شاخص dNBR منطقۀ موردمطالعه Figure (5) dNBR index map of the studied area ارزیابی نتایج بهدستآمده، امری ضروری است؛ بهطوریکه اگر نتایج ارزیابی نشوند، بدین معنی است که تفسیر کامل نیست. با توجه به هدف اصلی پژوهش (تشخیص مناطق گرفتار حریق)، پس از تفکیک و تشخیص مناطق آتشسوزی با استفاده از شاخصهای مورداستفاده، اقدام به طبقهبندی و تفکیک مناطق سوخته از سایر مناطق شد. در این راستا، از روش پیکسل پایه و الگوریتم بیشترین احتمال برای طبقهبندی استفاده شد که از 3 طبقۀ مناطق گرفتار حریق، مناطق مسکونی و سایر کاربریها در طبقهبندی استفاده شد. با استفاده از نتایج بهدستآمده از نقشۀ طبقهبندی (شکل 6)، وسعت مناطق گرفتار حریق در منطقۀ موردمطالعه، 13685 هکتار برآورد شد.
شکل (6) نقشۀ طبقهبندیشدۀ مناطق گرفتار حریق Figure (6) classified map of the fire affected areas
نتیجهگیری بهرهگیری از دانش سنجش از دور و تصاویر ماهوارهای در ارزیابی و مدیریت پایدار منابع طبیعی و بهتبع آن بررسی چالشها و تهدیدات منابع جنگلی، به دلیل ویژگیهای این مناطق و همچنین کارایی تکنیکهای سنجش از دوری، بهصورت گسترده در حال توسعه و تکامل است. امروزه بررسی قابلیت تصاویر ماهوارههای مختلف، در تشخیص و ارزیابی جنبههای گوناگون علوم جنگلداری از سوی محققان، در بسیاری از مناطق رویشی دنیا انجام شده و اطلاعات بسیار خوبی را در اختیار قرار داده است. خوشبختانه در این ارتباط در دهههای اخیر، تصاویر سنجندههای مختلف لندست بهآسانی و برای دورههای مختلف در دسترس بوده است و در صورتی که از قابلیت مناسبی در تشخیص مناطق جنگلی سوختهشده برخوردار باشد، بهراحتی در تعیین سطوح آتشسوزیهای مناطق مختلف مؤثر است و اطلاعات دقیقی را در اختیار قرار میدهد. این پژوهش با هدف بهکارگیری قابلیت و توان تصاویر ماهوارۀ لندست، شاخص NBR و dNBR در تفکیک و شناسایی عرصههای جنگلی سوختهشده از سایر پدیدهها، در بخشی از جنگلهای زاگرس، طراحی و اجرا شد. بررسی تغییرات پس از آتشسوزی و آگاهی از شدت سوختگی با استفاده از شاخص NBR تجزیهوتحلیل شد. بهعنوان نتیجهگیری کلی گفته میشود که با توجه به صحت نتایج طبقهبندی در تشخیص عرصههای سوختهشده از سایر مناطق با اعتماد بالا، بر توان و قابلیت مناسب تصاویر ماهوارۀ لندست و شاخص NBR و dNBR در تفکیک و تمایز مناطق جنگلی سوختهشده تأکید شده است. براساس نتایج نقشههای شدت آتش بهدستآمده از این پژوهش، ادعا میشود که توان تفکیک مکانی و طیفی تصاویر لندست بهخوبی پاسخگوی تهیۀ آمار و اطلاعات صحیح از عرصههای آتشسوزی و درنتیجه تهیۀ نقشۀ مناطق جنگلی سوختهشده در محدودۀ جنگلهای زاگرس ایران است؛ درنتیجه با توجه به نتایج بهدستآمده از نقشههای دمای سطح زمین مربوط به قبل و بعد از آتشسوزی که حاکی از افزایش 9 درجۀ سانتیگرادی دما در منطقۀ موردمطالعه و همچنین با توجه به استخراج مناطق دچار حریق توسط شاخصهای NBR و dNBR و استخراج مناطق دچار حریق با طبقهبندی به روش نظارتشده (الگوریتم بیشترین احتمال) و همخوانی و همسویی این مناطق با یکدیگر، وسعت مناطق دچار حریق در منطقۀ موردمطالعه 13685 هکتار برآورد میشود. ادب و همکاران(2012) در پژوهشی به مدلسازی خطر آتشسوزی در جنگلهای گلستان به کمک تصاویر مودیس و سیستم اطلاعات جغرافیایی توجه کردند. متغیرهای مورداستفاده براساس ضریب آتشسوزی به طبقاتی تقسیمشده با تجزیهوتحلیل دادهها در محیط GIS به سه منطقۀ خطر کم، متوسط، زیاد برای منطقه به دست آمد. نجفی و همکاران (1394) پژوهشی را با عنوان مدلسازی و تهیۀ نقشۀ خطر وقوع آتشسوزی جنگل و با استفاده از سنجش از دور در منطقۀ باغ شادی یزد بررسی کردند. در این پژوهش، از مدل تحلیل سلسلهمراتبی و معیارهای فیزیولوژی شامل شیب، ارتفاع دما و بارش و ... استفاده شد؛ همچنین از شاخصهای پوشش گیاهی نیز استفاده شد و در آخر منطقه به طبقات پرخطر و خطرناک همپوشانی زیادی با نقاط آتشگرفتۀ گذشته داشته است. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| مراجع | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
منابع
درویشی، شادمان، سلیمانی، کریم، رشیدپور، مصطفی (1398). تأثیر شاخصهای گیاهی و خصوصیات سطح شهری بر تغییرات دمای سطح زمین (مطالعۀ موردی: شهرستان سنندج). سنجش از دور و سامانۀ اطلاعات جغرافیایی در منابع طبیعی، دورۀ 10، شمارۀ 1، ص 35-17.
سرکارگراردکانی، علی، ولدان زوج، محمدجواد، منصوریان، علی (1388). تحلیل فضایی نیروی آتشسوزی مناطق مختلف کشور با استفاده از RS و GIS. محیط شناسی، دورۀ 35، شمارۀ 52، ص 34-25.
کاکه ممی، آزاد، قربانی، اردوان، اصغری سراسکانرود، صیاد، قلعه، احسان، غفاری، سحر (1399). بررسی رابطۀ تغییرات کاربری اراضی و پوشش گیاهی با دمای سطح زمین در شهرستان نمین. سنجش از دور و سامانۀ اطلاعات جغرافیایی در منابع طبیعی، دورۀ 11، شمارۀ 2، ص 48-27.
مختاری، محمدحسین، نجفی، احمد (1394). مقایسۀ روشهای طبقهبندی ماشین بردار پشتیبان و شبکۀ عصبی مصنوعی در استخراج کاربریهای اراضی از تصاویر ماهوارهای لندست TM. مجلۀ علوم آب و خاک، دورۀ ۱۹، شمارۀ ۷۲، ص 45-35.
ویسی، شادمان، ناصری، عبدعلی، حمزه، سعید، مرادی، پوریا (1395). برآورد دمای مزارع نیشکر با استفاده از الگوریتم پنجرۀ مجزا و تصاویر سنجندۀ OLI ماهوارۀ لندست 8. سنجش از دور و سامانۀ اطلاعات جغرافیایی در منابع طبیعی، دورۀ 7، شمارۀ 1، صص 40-27.
Aboelnour, M., & Engel, B. (2018). Application of remote sensing techniques and Geographic Information Systems to Analyze Land Surface Temperature in Response to Land Use/Land Cover Change in Greater Cairo Region, Egypt. Journal of Geographic Information System, 10, 57-88.
Chen, Y., Lara, M.J., & Hu, F.S. (2020). A robust visible near-infrared index for fire severity mapping in Arctic tundra ecosystems. ISPRS Journal of Photogrammetry and Remote Sensing, 159, 101–113.
Gondwe, S.V.C., Muchena, R., & Boys, J. (2018). Detecting land Use and Land Cover and Land Surface Temperature Change in Lilongwe City, Malawi. Journal of Remote Sensing and GIS, 9(2), 17-26.
Hussaina, M., Chen, D., Cheng, A., Wei, H., & Stenley, D. (2013). Change Detection from Remotely Sensed Images: From Pixel based to Object-based Approaches. Journal of Photogrammetry and Remote Sensing, 80, 91–106.
Kaufman, Y.J., Ichoku, C., & Giglio, L. (2003). Fire and Smoke Observation from the Earth Observation System MODIS Instrument-Products, Validation, and Operational Use. International Journal of Remote Sensing, 24(8), 1765-1781.
Kurnaz, B., Bayik, C., & Abdikan, S. (2019). Determination of Forest Fire Area by Using Satellite Images: Muğla Case. 3rd International Conference on Advanced Engineering Technologies, 19-21.
Labib, S.M., & Harris, A. (2018). The potentials of Sentinel-2 and LandSat-8 data in green infrastructure extraction, using object based image analysis (OBIA) method. European Journal of Remote Sensing, 51(1), 231-240.
Niu, R., & Zhai, P. (2012). Study on forest fire danger over Northern China during the recent 50 years. Climatic Change, 111(3-4), 723-736.
Novelli, A., Aguilar, M.A., Nemmaoui, A., Auilar, F.J., & Tarantino, E. (2016). Performance evaluation of object based greenhouse detection fromSentinel-2 MSI and Landsat 8 OLI data: A case study from Almería(Spain). International Journal of Applied Earth Observation and Geoinformation, 52, 403–411.
Parks, S.A., Dillon, G.K., & Miller, C. (2014). A New Metric for Quantifying Burn Severity: The Relativized Burn Ratio. Remote Sensing, 6, 1827-1844.
Quintano, C., Fernández-Manso, A., & Fernández-Manso, O. (2018). Combination of Landsat and Sentinel-2 MSI data for initial assessing of burn severity. International Journal Applied Earth Observation Geoinformation, 64, 221–225.
Sobrino, J.A., Jiménez-Muñoz, J.C., & Paolini, L. (2004). Land surface temperature retrieval from LANDSAT TM 5. Journal of Remote Sensing of Environment, 90(4), 434-440.
Wan, Z., & Dozier, J. (1996). A generalized split-window algorithm for retrieving land-surface temperature from space. Journal of Geoscience and Remote Sensing, 34(4), 892-905.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
آمار تعداد مشاهده مقاله: 1,123 تعداد دریافت فایل اصل مقاله: 502 |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||