چکیده:
پیشینه و هدف در طول دو دهه اخیر نیاز شدید به اطلاعات دمای سطح زمین جهت مطالعات محیطی و فعالیت های مدیریتی و برنامه ریزی، برآورد دمای سطح زمین را به یکی از موضوعات مهم علمی تبدیل کرده است. از سویی دیگر روش های مختلفی جهت تخمین دمای سطح زمین ارائه شده است که هرکدام نتایج متفاوتی را برای مناطق مختلف در پی داشته است. در این پژوهش الگوریتم هایی که در مطالعات مختلف هرکدام نتایج قابل قبولی داشته، انتخاب و مورد ارزیابی قرارگرفته است. در حوزه مطالعات حرارتی آنچه بهعنوان یک نقص اساسی در پایش دمای سطح زمین به شمار میآید، نبود ایستگاه های هواشناسی کافی جهت آگاهی از مقادیر دمایی در نقاط فاقد ایستگاه و محدودیت اطلاعاتی در تهیه داده های دمایی به خصوص برای مناطق وسیع است. منطقه موردمطالعه نیز با این کمبود رو به رو است و این محدودیت، اهمیت موضوع انتخاب شده برای این پژوهش جهت تخمین دمای سطح زمین با استفاده از فناوری سنجشازدور را بیشتر نمایان می سازد. هدف از این تحقیق، تخمین دمای سطح شهرستان اردبیل و ارزیابی دقت چهار الگوریتم تک کاناله، تک پنجره بهبودیافته، رابطه معکوس تابع پلانک و معادله انتقال تابش، مقایسه دقت دو ماهواره لندست 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 قابل پیش بینی بود.
Background and Objective Over 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 mono-window, 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: single-channel
algorithms, An Improved mono-window, 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 single-channel
algorithms, the inversion of Planck’s function and
RTE, the first method of emission and for the An
Improved Mono-Window 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 single-channel
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 single-channel algorithms, An Improved
mono-window, 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 single-channel algorithm was more accurate than
the other methods. Comparison of the results of the
single-channel 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 single-channel algorithm, the An Improved Mono-
Window 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 out-of-town 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.