تخمین دمای سطح اراضی اردبیل با استفاده از تصاویر لندست و ارزیابی دقت روش های برآورد دمای سطح زمین با داده های میدانی

پیشینه و هدف در طول دو دهه اخیر نیاز شدید به اطلاعات دمای سطح زمین جهت مطالعات محیطی و فعالیت های مدیریتی و برنامه ریزی، برآورد دمای سطح زمین را به یکی از موضوعات مهم علمی تبدیل کرده است. از سویی دیگر روش های مختلفی جهت تخمین دمای سطح زمین ارائه ‌شده است که هرکدام نتایج متفاوتی را برای مناطق مختلف در پی داشته است. در این پژوهش الگوریتم هایی که در مطالعات مختلف هرکدام نتایج قابل قبولی داشته، انتخاب و مورد ارزیابی قرارگرفته است. در حوزه مطالعات حرارتی آنچه به‌عنوان یک نقص اساسی در پایش دمای سطح زمین به شمار می‌آید، نبود ایستگاه‌ های هواشناسی کافی جهت آگاهی از مقادیر دمایی در نقاط فاقد ایستگاه و محدودیت اطلاعاتی در تهیه داده های دمایی به خصوص برای مناطق وسیع است. منطقه موردمطالعه نیز با این کمبود رو به رو است و این محدودیت، اهمیت موضوع انتخاب ‌شده برای این پژوهش جهت تخمین دمای سطح زمین با استفاده از فناوری سنجش‌ازدور را بیشتر نمایان می سازد. هدف از این تحقیق، تخمین دمای سطح شهرستان اردبیل و ارزیابی دقت چهار الگوریتم تک کاناله، تک پنجره بهبودیافته، رابطه معکوس تابع پلانک و معادله انتقال تابش، مقایسه دقت دو ماهواره لندست 5 و لندست 8 در برآورد دمای سطح زمین. مواد و روش هادر این پژوهش از سه نوع داده استفاده‌شده است؛ تصاویر ماهواره لندست 5 و 8، داده‌ های دو ایستگاه هواشناسی، و داده ‌های زمینی برداشت‌شده با دماسنج دیجیتالی. تصاویر مورداستفاده از دو ماهواره لندست 5 و لندست 8 بافاصله زمانی 19 ساله انتخاب ‌شده است. داده های هواشناسی مورداستفاده نیز از دو ایستگاه سینوپتیک موجود در محدوده موردمطالعه اخذ گردید. علاوه بر دمای سطح زمین، داده های رطوبت نسبی، حداقل دما و حداکثر دمای 24 ساعت نیز در دو تاریخ مدنظر اخذ گردید، همچنین دونقطه از منطقه موردمطالعه انتخاب و دمای سطح زمین در موقعیت این دو ایستگاه هم‌زمان با عبور ماهواره با استفاده از دو دماسنج دیجیتالی ثبت شد. جهت مدل‌سازی تابش و میزان انتقال اتمسفری از نرم‌افزار محاسبه‌گر تحت وب modtran استفاده‌شده است. توان تشعشعی با دو روش گسیلمندی بر اساس شاخص ndvi و گسیلمندی بر اساس حدآستانه گذاری ndvi و دمای سطح زمین با چهار الگوریتم تک کانال، تک پنجره بهبودیافته، رابطه معکوس تابع پلانک و معادله انتقال تابشی  با استفاده از باند 6 لندست 5 و باند 10 لندست 8 در نرم افزارmatlab  برای دو سال 2000 و 2019 کدنویسی گردید. درنهایت دقت الگوریتم ها با استفاده از داده های دمای سطح ایستگاه سینوپتیک و نمونه برداری میدانی مورد ارزیابی قرار گرفت. نتایج و بحث نتایج نشان داد  که برای سه الگوریتم تک کانال، رابطه معکوس تابع پلانک و rte، روش اول گسیلمندی و برای الگوریتم تک پنجره بهبودیافته روش دوم گسیلمندی از دقت بالاتری برخوردار بوده است. داده های دمای سطح اخذ شده از ایستگاه های هواشناسی در سال 2000 ازلحاظ زمانی 12 دقیقه اختلاف و برای سال 2019 اختلاف 4 دقیقه ای بازمان عبور ماهواره دارد. ایستگاه اول هواشناسی تا حدودی در محدوده شهری واقع‌شده است. نتایج نشان داد که مهم‌ترین عامل بیشتر بودن اختلاف ایستگاه اول با lst برآورد شده در مقایسه با ایستگاه دوم همین عامل باشد، چراکه ناهمگونی پیکسل ها و تغییرات زیاد سطوح در محدوده شهری باعث تداخل ارزش پیکسل ها و به دنبال آن احتمال بروز خطا در برآورد دمای سطح در محدوده انسان ساز شهری را بالا  می برد. برای ایستگاه زمینی نیز دونقطه با محیطی همگن و خارج از محدوده شهری با کاربری کشاورزی (یونجه) و کاربری بایر که محصول آن برداشت‌ شده بود، انتخاب و دمای سطح آن‌ها هم‌زمان با عبور ماهواره اندازه گیری شد. نتایج خروجی تخمین دمای سطح زمین با دو ایستگاه سینوپتیک و دو ایستگاه زمینی مورد مقایسه و ارزیابی قرار گرفت. در هر دو تاریخ الگوریتم تک کانال کمترین اختلاف را با ایستگاه های ثبت دما نشان داد. نتیجه گیری در این پژوهش با استفاده از تصاویر ماهواره لندست 5 و لندست 8 چهار الگوریتم برآورد دمای سطح زمین شامل روش های تک کانال، تک پنجره بهبودیافته، رابطه معکوس تابع پلانک و ربطه انتقال تابش کدنویسی و نقشه های دمای سطح زمین شهرستان اردبیل برای دو سال 2000 و 2019 در محیط نرم افزار متلب کدنویسی و استخراج گردید. باند 6 ماهواره لندست 5 برای سال 2000 و از باند 10 ماهواره لندست 8 به دلیل مقدار نویز کمتر نسبت به باند 11 و نزدیکی به مقدار 9.66 که بیشترین تابش زمین برای سال 2019 استفاده شد. مقایسه نقشه های دمای سطح حاصل از الگوریتم ها با ایستگاه های سینوپتیک و زمینی نشان داد که در هر دو سال 2000 و 2019 الگوریتم تک کانال دقت بیشتری نسبت به بقیه روش ها داشته است. مقایسه نتایج روش تک کانال، نشان از اختلاف 2.5+ و 2 با ایستگاه های 1 و 2 برای سال 2000 و اختلاف دمای 1.3+، 0.9+، 1 و 0.9 به ترتیب با ایستگاه های 1، 2، 3 و 4 برای سال 2019 را نشان می دهد. استفاده مستقیم از ضرایب انتقال‌پذیری اتمسفر در فرآیند روش تک کانال، در بالا بودن دقت این روش موثر بوده است. ازنظر دقت بعد از الگوریتم تک کانال، به ترتیب روش تک پنجره بهبودیافته، الگوریتم rte و درنهایت الگوریتم رابطه معکوس تابع پلانک قرار گرفتند. نتایج مقایسه خروجی هر چهار الگوریتم با داده های ایستگاه های 1، 2، 3 و 4، نشان از دقت بالاتر ایستگاه های زمینی برداشت‌شده با دماسنج دیجیتالی نسبت به داده های ایستگاه های هواشناسی دارد، ازجمله دلایل آن می توان به قرارگیری ایستگاه های هواشناسی (به خصوص station_1) در محدوده شهری با توجه به ناهمگن بودن محیط شهری و امکان تداخل پیکسلی و تداخل دمایی کاربری ها اشاره کرد، درحالی‌که ایستگاه های زمینی از محدوده خارج از شهر و از محیطی با پیکسل های همگن (بایر و کشاورزی) انتخاب گردید. همچنین نتایج هر چهار الگوریتم مستخرج از تصویر لندست 8 در مقایسه با نتایج چهار الگوریتم حاصله از تصویر لندست 5، دقت بیشتری را نشان می دهد و با توجه به بهبود توان تفکیک مکانی سنجنده tirs نسبت به tm، دقت بیشتر خروجی های سنجنده tirs قابل پیش بینی بود.

Estimation of Ardabil land surface temperature using Landsat images and accuracy assessment of land surface temperature estimation methods with ground truth data

Background and ObjectiveOver the past two decades, the intense need for land surface temperature information for environmental studies and management and planning activities has made estimating the land surface temperature one of the most important scientific topics. On the other hand, different methods have been proposed to estimate the land surface temperature, each of which has resulted in different results for different regions. In this study, the algorithms that have had acceptable results in different studies have been selected and evaluated. In the field of thermal studies, what is considered as a major defect in monitoring the land surface temperature is the lack of sufficient meteorological stations to know the temperature values in places without stations and information limitations in preparing temperature data, especially for large areas. The study area is also facing this shortage, and this limitation further highlights the importance of the topic selected for this study to estimate the surface temperature using remote sensing technology. Verification and validation of results obtained from estimating the land surface temperature are other basic and discussed topics in thermal studies. The purpose of this study is an estimation of temperature in Ardabil city and evaluate the accuracy of the four singlechannel algorithms, the improved monowindow, the Planck's inversion function method and the radiative transfer equation (RTE) method, to compare the accuracy of the two Landsat 5 and Landsat 8 satellites in estimating the land surface temperature. Materials and Methods Three types of data have been used in this study; Landsat 5 and 8 satellite images, data of two meteorological stations and ground data harvested with a digital thermometer. The images used are from the two satellites Landsat 5 and Landsat 8 with a time interval of 19 years. The meteorological data used were obtained from two synoptic stations in the study area. In addition to land surface temperature, relative humidity, minimum temperature and maximum temperature data of 24 hours were also obtained on two dates. Also, two points of the study area were selected and land surface temperature in the position of these two stations simultaneously with the satellite Recorded from two digital thermometers. MODTRAN web version calculator software version 6 has been used to model the radiation and the amount of atmospheric transmission. Emissivity with two methods of LSE methods based on NDVI and LSE NDVI Thresholds Method and land surface temperature with four algorithms: singlechannel algorithms, An Improved monowindow, inversion of Planck’s function and radiative transfer equation using band 6 Landsat 5 and band 10 Landsat 8 bands. It was coded in MATLAB software for 2000 and 2019. Finally, the accuracy of the algorithms was evaluated using synoptic station surface temperature data and field sampling. Results and Discussion The collected data and results are analyzed and while presenting the output maps, the accuracy of the methods with terrestrial and meteorological data as well as the accuracy of Landsat 5 and Landsat 8 satellites in estimating the land surface temperature has been compared and evaluated. The results showed that for the three singlechannel algorithms, the inversion of Planck’s function and RTE, the first method of emission and for the An Improved MonoWindow algorithm, the second method of emission had a higher accuracy. Land surface temperature data obtained from meteorological stations in 2000 differ by 12 minutes in terms of time and by 2019 differ by 4 minutes in terms of satellite transit time. The first meteorological station is located somewhat within the city limits and according to the results, it seems that the most important factor is the greater difference between the data of the first station and the estimated LST compared to the second station is the same factor because the heterogeneity of pixels and large changes in levels in urban areas interfere with a pixel value. And subsequently increases the likelihood of errors in estimating surface temperature within the urban anthropogenic range. For the ground station, two points with a homogeneous environment and outside the urban area with agricultural use (alfalfa) and barren use of the harvested product were selected and their surface temperature was measured at the same time as the satellite. The output results of land surface temperature estimation were compared and evaluated with two synoptic stations and two ground stations. In both histories, the singlechannel algorithm showed the least difference with the temperature recording stations. Conclusion In this research, using Landsat 5 and Landsat 8 satellite images, four algorithms for estimating the land surface temperature of the earth, including singlechannel algorithms, An Improved monowindow, inversion of Planck’s function and radiative transfer equation and land surface temperature maps of Ardabil city for two 2000 and 2019 were coded and extracted in MATLAB software environment. The band 6 Landsat 5 satellite was used for 2000 and the band 10 Landsat 8 satellite was used for 2019 due to less noise than the 11th band and the proximity of 9.66 (which is the highest radiation in this range). Comparison of land surface temperature maps obtained by the algorithms with synoptic and ground stations showed that in both 2000 and 2019, the singlechannel algorithm was more accurate than the other methods. Comparison of the results of the singlechannel method with the stations shows a difference of  +2.5 and 2 with stations 1 and 2 for the year 2000 and a temperature difference of  +3.3, +0.9, 1 and 0.9. Shows stations 1, 2, 3 and 4 for 2019, respectively. It seems that the direct use of atmospheric transmittance coefficients in the singlechannel method process has been effective in the high accuracy of this method. In terms of accuracy, after the singlechannel algorithm, the An Improved MonoWindow method, the RTE algorithm, and finally the Planck function inverse correlation algorithm were placed, respectively. The results of comparing the output of all four algorithms with the data of stations 1, 2, 3 and 4, show that the ground stations harvested with a digital thermometer are more accurate than the data of meteorological stations. One of the reasons for this is the location of meteorological stations (especially, Station_1) in the urban area due to the heterogeneity of the urban environment and the possibility of pixel interference and temperature interference of land uses, while ground stations from the outoftown area. And was selected from an environment with homogeneous pixels (barren and agricultural). Also, the results of all four algorithms extracted from the Landsat 8 image show more accuracy compared to the results of the four algorithms obtained from the Landsat 5 image, and due to the improved spatial resolution of the TIRS sensor compared to the TM, the TIRS sensor output is more accurate, It was predictable.