Accuracy Evaluation and Parameter Analysis of Land Surface Temperature Inversion Algorithm for Landsat-8 Data

,


Introduction
Land surface temperature (LST) is an essential parameter in many research fields, such as ecology, climatology, urban thermal environment, urban heat island, and hydrology [1][2][3].Especially in recent years, more attention has been paid to carbon dioxide emissions [4,5], polar climate warming and glacier melting [6,7], urban heat island (UHI) [8][9][10], urban heatwave and human living health [11], forest fires [12], drought monitoring [13], soil evapotranspiration [14], and precipitation and surface runoff [15,16].e traditional LST measurement method is based on the temperature measurement data of meteorological stations for further spatial interpolation, which is a classical pointbased measurement method [17].Although this method has high accuracy, the point-based temperature measurement method cannot represent the regional temperature.e timeliness of the method is insufficient.With the development and research of space-based infrared remote sensing, it is more and more popular to use satellite remote sensing to retrieve regional and global surface temperatures.
In many satellite observation systems equipped with thermal infrared sensors, Landsat satellites have become one of the most widely used satellites due to their observation history of more than 40 years.Since its launch on February 11, 2013, the newest generation of Landsat-8 has been in orbit for more than eight years.It has sent back a large amount of Earth observation data at a rate of 550 images per day.Together with Landsat-7 ETM+, it forms a repeated observation period with an interval of 8 days, which provides precious data for monitoring ecological environment changes on the Earth's surface.Because of its advantages of high spatial resolution, large scanning field of view, and accessible data, researchers have applied it to further scientific research and ecological construction.
Since the launch of Landsat-8, various processing documents and parameters have been posted on the Landsat.Many researchers have developed a variety of LST inversion algorithms based on Landsat-8.
e main LST inversion algorithms include the Radiative Transfer Equation (RTE) [18,19], Single Channel (SC) algorithm [20,21], Mono Window (MW) algorithm [22], and Split Window (SW) algorithm [21,[23][24][25].In these surface temperature inversion algorithms, multiple basic parameters are needed for input [26].e multiple parameters are estimated variables rather than standard ones.Accurate and reliable inversion algorithms should calculate the accuracy of the LST data used in the study to make the use of the product possible.If the inversion algorithm has a significant systematic error or a significant instability error, the algorithm or product will not be popular among users [27].
In this study, LST inversion, spatial analysis, and multiscale analysis methods were considered for Landsat data representing four seasons.e commonly used LST inversion algorithms are evaluated.First of all, we designed nine remote sensing inversion schemes for LST.We fitted the LST calculated by the satellite with the measured temperature values at the ground stations at the same time.Based on the same scene remote sensing data, the accuracies of the nine inversion algorithms were compared and analyzed.e accuracy of each LST inversion algorithm on the time scale was analyzed separately on the images of different seasons.
en, the inversion accuracies of all the LST inversion methods on the regional scale (inside city and outside city) of the same surface features are compared and analyzed.Finally, we test the dependence of each inversion algorithm on parameters by equal step size.e aims of this study are as follows: (1) evaluate the errors of multiple inversion algorithms on the same scene data, (2) evaluate the error of each inversion algorithm in seasons (spring, summer, autumn, and winter), (3) evaluate the regional error of each inversion algorithm, and (4) perform dependence intensity and sensitivity analysis of the evaluation algorithm on parameters.

Study Areas.
e study area was selected as an area covered by Landsat-8 data (orbit number 123, 032).e longitude and latitude of the data center are 40.3°N and 116.7 °E.It mainly covers Beijing, the capital of China.e climate type in the study area is semihumid and semiarid monsoon climate in the warm temperate zone, with four distinct seasons, hot and rainy in summer and cold and dry in winter.It provides a natural, convenient condition for studying the difference of inversion algorithms on a seasonal scale.

Data Compilation.
is study contained four pieces of Landsat-8 data based on USGS (United States Geological Survey, https://www.usgs.gov/).Landsat-8 has 11 bands, of which seven multispectral bands have a spatial resolution of 30 m. rough resampling, the spatial resolution of the two thermal infrared bands is also 30 m. e revisit period is 16 days for global coverage.Since the surface reflectance is affected by the atmospheric conditions at the time of image acquisition and it is also a potential factor for studying the surface thermal model, when screening the data, we select the imaging data under good weather condition, which was a clear day with less than 2% cloud cover, and use the thermal band data to retrieve the LST conduct research.In order to investigate the accuracy of the Landsat-8 land surface temperature inversion algorithm in season, the imaging dates of January 31, 2017 (winter), May 7, 2017 (spring), July 10, 2017 (summer), and October 30, 2017 (autumn), were selected.All images were named by the 8-digit number of years, months, and day of the collection date for the convenience of recording.For example, the image on May 7, 2017, as shown in Figure 1, was named 2017-05-07.
is study also collected and sorted the surface temperature data of 20 meteorological stations in the study area.
e basic information of the meteorological stations is shown in Table 1.
e positions of these meteorological stations in the study area are shown in the red dots in Figure 1.
e temperature data were collected from the China Meteorological Administration and measured simultaneously as remote sensing images.To study the accuracy of the inversion algorithm of the LST inside and outside the city, points of interest (POIs) were established for four kinds of ground objects (water, vegetation, soil, and building), and 20 POIs were established for each ground object.Pure pixels were selected for each point of interest from Landsat-8 multispectral data.High-resolution Google Earth images also examined the selection of all POIs.As shown in Figure 1, those outside the city are represented by green dots, while yellow dots represent those inside the city.

Multispectral Image Processing.
Using the Landsat-8 website, the user can convert the Digital Numbers (DNs) of the image to the spectral radiance of the top of the atmosphere (TOA).By applying the equation where L λ (Watts•m −2 •srad −1 µm −1 ) refers to the TOA spectral radiance, M L is the multiplicative rescaling factor of the corresponding band, it can be obtained from the header file of the data with the field "Radiance_Mult_Band_x � ," x stands for band number.Q cal is the pixel values, and A L is the offset of the data.e multispectral bands of Landsat-8 data can also be converted to the reflectivity of the top of the atmosphere by the following equation: where ρ λ ′ is the reflectance of the top of the atmosphere at band λ without correction of the solar angle, M ρ is the reflectivity adjustment factor of band λ, and it can be obtained from the header file of the data with the field where θ z is the solar zenith angle of the image center, and solar altitude angle θ E can also be used, but sine function should be used.We used the solar altitude angle provided in the image for our calculations.
In this study, the COST model atmospheric correction algorithm of Pat and Chavez [28] is adopted to conduct Landsat-8 atmospheric correction.COST requires few parameters, and the correction is mainly carried out with the data of the image itself [29] by the following equation: where Q h is the modification value of atmospheric influence, which can be obtained by the darkest pixel method, τ � cos[(90 − θ z )π/180], and τ is the atmospheric transmittance estimated based on θ z .After COST atmospheric correction, NDVI (normalized difference vegetation index) is calculated by the following equation: where ρ nir COST is the reflectance image of the near infrared band and ρ red COST is the reflectance image of the red band.

ermal Infrared Image
Processing. e brightness temperature (BT) can be calculated using the following equation: where T is the BT; K 1 and K 2 for the Landsat-8 band 10 are 774.89Watts•m −2 •srad −1 µm −1 and 1321.08K, respectively.L λ is the thermal infrared band after radiation calibration and COST atmospheric correction.

Estimation of Land Surface Emissivity.
In this study, we applied two LSE (LSE1 and LSE2) estimation methods; both methods are based on NDVI.e first method LSE1 refers to the NDVI threshold method proposed by Sobrino [30] calculated by the following equation: where ε is the LSE and P v is calculated according to NDVI using the following equation: where NDVI soil is the NDVI value of bare soil or areas without vegetation coverage and NDVI veg is the NDVI value of pixels completely covered by vegetation, namely, the NDVI value of pure vegetation pixels.Empirical values NDVI veg � 0.6 and NDVI soil � 0.05 were taken; that is, when the NDVI of a pixel is greater than 0.60, P v value is 1.0, and when NDVI is less than 0.05, P v value is 0. e other method LSE2 is to firstly divide the surface into water, natural surface, and urban area according to the method proposed by Qin et al. [31] and then calculate the LSE for the three types of surface separately.e calculation is by the following equation: However, the concept of natural surfaces and urban areas is vague and cannot be quantitatively calculated.Based on remote sensing images of Beijing and its surrounding areas, we conducted local experiments and found that when NDVI is less than 0.05, it corresponds to water; when NDVI is greater than or equal to 0.05 or less than or equal to 0.6, it corresponds to urban areas; and when NDVI is greater than 0.6, it corresponds to the natural surface.e calculation is in the following equation:  [32,33].

RTE Inversion Methods.
e RTE is a direct method for LST retrieval using a single TIR band. is can be given by the following equation: where is the brightness value of TIRS received by satellite sensor and λ is the tenth band of TIRS.At present, the accuracy of the eleventh band calibration of TIRS is still not quantified, and USGS still does not encourage the use of eleventh band for relevant surface temperature inversion, so we only consider the tenth band in the RTE and SC.ε refers to the LSE, B(T s ) is the blackbody radiance energy, and B(T s ) at a temperature of T s is calculated by the inversion of the following equation: and, eventually, T s can be obtained from the inversion of Planck's law as in the following equation: 3.4.2.SC Inversion Methods.Jiménez-Muñoz et al. [20] proposed a revised version of SC for LST retrieval using Landsat TIR data.Concerning the SC, T s is obtained from the following equations: where b c is equal to 1320 K for Landsat-8 band 10 and L and T are obtained by equation ( 1) and equation ( 6), respectively.ψ 1 , ψ 2 , ψ 3 is a function of water vapor content (w) expressed in the following equation: When w ≥ 3g • cm − 2 , the accuracy of LST will be significantly affected, so Jiménez-Muñoz suggests using equation ( 18) to calculate ψ 1 , ψ 2 , ψ 3 :

MW Inversion Methods. Mono Window (MW) algorithm was originally developed for Landsat TM data.
Later, some scholars evaluated the ability of MW to process Landsat ETM + data and Landsat-8 data, which can be used in the following equation: where T a is the effective mean atmospheric temperature and a and b values for the Landsat-8 band 10 are −67.355and 0.458606, respectively.
where T 10 and T 11 are the brightness temperatures of band 10 and band 11 and ε m is the average LSE of band 10 and band 11.Δε is the LSE interpolation of band 10 and band 11.
In this study, the LSE of main ground objects was obtained through ASTER spectral library, and the estimated value of three NDVI segments in equation ( 10) was matched through the calculation of the average value of LSE of relevant ground objects.Table 3  Rozenstein et al. [25] proposed an SW algorithm for inversion of land surface temperature using Landsat-8 TIR data based on the MW algorithm of QIN [23].e algorithm is expressed by the following equations:

Advances in Meteorology
where the parameters determined by atmospheric transmittance and LSE, and a i , b i are, respectively, the regression coefficients of band 10 and band 11 determined according to different temperature ranges; this is shown in Table 4.In this study, parameters with temperature ranging from 0 °C to 40 °C were selected.ε i is the LSE of band 10 and band 11. τ i (θ) is the atmospheric transmittance corresponding to the solar zenith angle.Rozenstein et al. provided the relationship between atmospheric transmittance and atmospheric water vapor content, as shown in Table 5; namely, atmospheric transmittance was estimated by atmospheric water vapor content, which was also the calculation method in this paper.In this study, the parameters determined by mid-latitude summer are selected.

Validation Method and Performance Metrics.
In this study, based on four commonly used LST inversion methods and different parameters, a total of nine schemes are designed, which are shown in Table 6.
e T-based technique, examined in this research, is a direct way of comparing the satellite-derived LST with measured temperatures at meteorological stations [34][35][36][37][38].
e capability of the T-based technique depends mostly on the accuracy of the measured temperatures at meteorological stations and pixel scale effect.
e major benefit of the T-based method is that the accuracy of LST algorithm can be accurately evaluated [39].In this study, satellite-derived LST and station-based LST were analyzed considering the performance metrics such as Root Mean Square Error (RMSE) and average Bias [40].e formulas of these metrics are as follows: where T L8 and T Station are the Landsat-8-derived LST and station-based LST, respectively, and n refers to the number of pieces of data; in this study, n is equal to 20.
Inputs of parameters such as LSE, up radiance, down radiance, atmospheric transmittance, atmospheric water and gas content, and atmospheric average temperature are needed in the inversion algorithm.Each of the above parameters is not directly measured by the instrument but has been estimated and processed.e errors of these inputs of parameters will affect the accuracy of LST inversion.Sensitivity analysis is an application of how the errors (numerical, statistical, or otherwise) of the model output are divided and allocated to different sources of uncertainty in the model input [41].e following equation is utilized: where x represents one of the input parameters and Δx is the step size of the parameter setting.In the study species, atmospheric transmittance (τ), emissivity, and NDVI were set as 0.01, upwelling radiance and downwelling radiance were set as 0.1 unit, water vapor content (w) was set as 0.1 unit, mean surface temperature was set as 1 unit, and S e was the LST difference calculated for each increase in step size; T s (x + Δx) and T s (x) refer to the LST calculated for "x + Δx" and "x," respectively.

Results and Analysis
In order to accurately evaluate the accuracy of the LST inversion algorithm, measured data from 20 meteorological stations in the study area were used to evaluate the stability e results are shown in Figure 2. By comparing the inversion results of each algorithm with the measured data, the statistical Root Mean Square Error (RMSE) and average Bias are shown in Table 7 of Appendix.It can be obtained from Figure 2 and Appendix's Table 7 that the inversion results of the LST algorithm are entirely consistent with the graph trend of the measured data of meteorological stations.e high elevation meteorological sites (e.g., Foyeding, Yanqing, Tanghekou, and Zhaitang) were low.In the graphic display, the line with the measured data of adjacent stations shows a concave shape.All LST inversion algorithms also have the same characteristics.In terms of seasonal scale, the inversion results of the nine algorithms in winter have the highest matching degree with the measured data of meteorological stations.e values of RMSE and Bias were relatively small.Followed by the data in autumn, the matching degrees in spring and summer were inferior.e values of RMSE and Bias were relatively large.In winter, the inversion results of LST1, LST2, LST3, LST4, and LST5 were the closest to the measured data of meteorological stations, with difference of less than 2 °C, followed by the inversion results of LST6, which were 3 °C higher than the measured data, and the difference between the inversion results of LST7, LST8, and LST9 and the measured data was about 7 °C.In spring and summer, almost all the algorithm inversion results were higher than the measured temperature.At this time, the closest to the measured temperature were LST1, LST2, LST4, and LST5, with difference of about 5 °C, followed by LST3 inversion results, about 10 °C.e inversion results of LST8 and LST9 were 20 °C higher than the measured results, and the reliability of the results will be seriously questioned.In autumn, the inversion results of LST1, LST2, LST4, and LST5 were slightly lower than the measured data by 2 °C, the inversion results of LST3 were 5 °C higher than the measured data, and the inversion results of LST6, LST7, LST8, and LST9 were 8 °C higher than the measured data.
From the perspective of the algorithm inversion results, the inversion results of LST1, LST2, LST4, and LST5 were the closest to the measured data of meteorological stations, with the inversion error within ±5 °C, followed by the inversion results of LST3 at about 10 °C.e worst are the inversion results of LST6, LST7, LST8, and LST9, with the inversion error around 20 °C.In general, in the season with low temperature, the error of the inversion results of the algorithm was smaller than that of the measured data.With the increase of temperature, the inversion results of all the inversion algorithms were higher than those of the measured data.LST1, LST2, LST4, and LST5 have a relatively small error, followed by LST3 LST6, LST7; the inversion errors of LST8 and LST9 were relatively large.
We performed linear fitting between nine LST inversion results and measured temperatures in four seasons.e fitting results are shown in Figure 3.According to the statistical fitting results, it can also be seen that, in winter, all the relevant results were relatively ideal, with R 2 greater than 0.73.Except for the two SW algorithms LST6 and LST7, the fitting slopes of the other algorithms were all less than 1. e fitting results in autumn show that R 2 was around 0.5, and all the fitting slopes were less than 1.However, the relevant intercept data are relatively large.In spring and summer, the proper slopes were more significant than 1. e absolute value of the fitting intercept was also relatively large.It is important to note that the fitting results of the LST7 algorithm in spring, summer, and autumn show no significant phenomenon; that is to say, the error between the inversion results of the LST7 algorithm and the measured results is unacceptable.
Overall, the inversion results of the RTE and SC algorithm calculated by LSE1 and LSE2 have an excellent fitting effect with the measured data; the inversion results of the SC algorithm are calculated by w. e inversion result of MW was higher compared to the measured temperature.e inversion effect of SW was not particularly ideal, which may be caused by the instability of the radiation calibration of the 11th band of TIRS.

Results on Different POIs.
It is well known that intensive human activities affect the surrounding environment.Industrial production and emissions affect atmospheric transmittance, water vapor content, and upward and downward radiation.All LST inversion algorithms need these parameters to support the calculation.erefore, it is necessary to explore whether the inversion results of similar surface features with the same reflectivity in densely populated cities and outside cities close to the natural surface will be affected by various inversion algorithms and to what extent.e four selected ground objects (20 POIs for each ground object) were analyzed inside and outside the city.Considering the statistical error, 10%-90% of the points were selected as the upper and lower whiskers of the box map, and the upper and lower boundaries of the box represented 25%-75% of the values, respectively.e result is shown in Figure 4.
As shown in Figure 4, in the same season, the same LST inversion algorithm results show that the height of the box map inside the city was less than the height of the box map outside the city.e standard deviation of inversion results of all LST inversion algorithms on POI inside the city was less than the standard deviation of inversion results outside the city.On the one hand, the main reason may be the urban heat island effect and on the other hand it may be the influence of intensive human activities inside the city on the parameters of the inversion algorithm.In addition, in a relatively narrow space, using the same estimated parameters for calculation, the inversion results are relatively concentrated.However, the outer space of the city is relatively broad, which belongs to the natural surface.ere are differences in estimated atmospheric parameters, and the inversion results will be different in different outer areas.As a result, the inversion results of land surface temperature are relatively scattered.
Both inside and outside the city, the height of the water box chart is minimal.e POI inversion temperature difference between inside the city and outside the city is about 1 °C in autumn.e average inversion temperature of POI  inside the city is about 3 °C higher than that outside the city in other seasons.It shows that the standard deviation of all LST inversion algorithms on POI of the water body is minimal, which also conforms to the characteristics of the water body being relatively stable, and LSE is close to the black body.Surprisingly, the height of the box chart of vegetation in all categories is considerable.Among all the POI inversion temperatures of vegetation outside the city in spring, the POI inversion temperature difference of vegetation outside the city is about 5 °C (this value does not count LST7 because it is precarious).
e POI inversion temperature difference of other vegetations is 8 °C.It indicates that the standard deviation of all LST inversion algorithms on vegetation POI is substantial, indicating that inversion algorithms are more sensitive to vegetation inversion, which will be explained in detail in the next section.In addition, it is also possible that when selecting pure pixels of vegetation, the vegetation types of pixels are different.Some pixels are grassland, some are woodland, and some are shrubs.Although these are all vegetation pixels, their reflectivity is inconsistent, which may lead to specific errors.
e box graph of soil POI temperature retrieved by the SW (LST6 and LST7) algorithm is relatively large; in particular, the LST7 algorithm is very obvious.e height of the box graph of soil POI temperature retrieved by other algorithms is relatively tiny, with a difference of about 4 °C outside the city and 3 °C inside the city.It shows that the standard deviation of all LST inversion algorithms on soil POI is minimal.It is worth noting that, in the selection of soil POI in this paper, the soil is dry and bright areas were selected instead of the dark soil with high moisture content to avoid the error caused by the significant difference in reflectivity.
Except for LST6 and LST7, the height of the box graph of building POI inversion by other LST inversion algorithms is small, with a difference of about 5 °C outside the city and 3 °C inside the city, indicating that the standard deviation of inversion results of all LST inversion algorithms on building POI is also small.

Sensitivity Analysis of Model Parameters.
is paper adopts the control variable method to conduct sensitivity analysis for different parameters to explore the sensitivity of various LST inversion algorithms to parameters.First, assume all the DN values on the image are the same (in this study, take the average of the band after atmospheric correction; band 10 is set to 10.4 and band 11 is set to 9.4).
en select a group of parameters as variables, and the rest as invariants to control the variation of variable parameters for sensitivity analysis.Finally, the sensitivity of the parameters was analyzed.As the data with good weather conditions were selected in the remote sensing image screening previously, the interval of atmospheric transmittance is set as 0.5∼1.0 with a step size of 0.01 unit, and the interval of surface reflectance is set as 0.9∼1.0 with a step size of 0.01 unit.e interval of DNVI is set as −0.2∼1.0 with a step size of 0.01 units.
e interval of upward and downward radiation is set as 0∼5, and the step size is set as 0.1 units; the interval of atmospheric water vapor content is set as 0-2.5, and the step size is set as 0.1 units; the interval of average surface temperature is set as 0-100, and the step size is set as 1 unit.e influence of each parameter on the results in different algorithms was obtained by controlling variables.For the comparability of the results, LST was normalized, as shown in Figure 5.
It can be seen from Figure 5(a) that the influence of atmospheric transmittance on the inversion algorithms of RTE, SC, and MW is a logarithmic function.With the increase of atmospheric transmittance, the inversion of land surface temperature gradually decreases.e atmospheric transmittance has a minor effect on the curvature of RTE (LST1 and LST2) and a more significant effect on the curvature of SC (LST4 and LST5) and MW (LST8 and LST9).Figure 5(b) shows that the influence of specific emissivity on all algorithms is linear.As the specific emissivity increases, the land surface temperature decreases linearly.When the emissivity increases by 0.1 units, the inversion surface temperature decreases by 1 unit.Figure 5(c) shows that the influence of NDVI on LST1, LST3, LST4, and LST8 is piecewise linear.Segmentation is because, during parameter setting, specific emissivity is segmenting according to NDVI mandatory threshold.Specific emissivity is set as a fixed value for parts with NDVI less than 0.05 and DNVI greater than 0.6.e inversion temperature is also a fixed value.For the part with DNVI greater than 0.05 and less than 0.6, the specific emissivity changed with the change of NDVI.e inversion temperature decreased linearly with the increase of NDVI.e influence of NDVI on LST2, LST5, and LST9 is also segmenting.However, when NDVI is more than 0.05, the influence of NDVI on LST presents a quadratic function; when NDVI is equal to 0.4, the function has a minimum value; when NDVI is more than 0.05 and less than 0.4, the inversion LST decreases the NDVI.When the NDVI is more than 0.4, the inversion LST increases the NDVI.
In Figure 5(d), the influence of uplink radiation on the SC (LST4 and LST5) algorithm is linear.With the increase of uplink radiation, the land surface temperature inversion by the SC algorithm decreases linearly.
e influence of upline radiation on the RTE (LST1 and LST2) algorithm is logarithmic.However, the curvature of the curve is minimal.With the increase of upline radiation, the surface temperature retrieved by the RTE algorithm also decreases.As shown in Figure 5(e), the influence of downward radiation on the SC and RTE algorithm is linear.As downward radiation increases by 1 unit, the inversion land surface temperature decreases by 0.2 units.As shown in Figure 5(f ), the influence of atmospheric water vapor content on LST3 and LST6 is linear.As the atmospheric water vapor content increases by 0.1 units, the surface temperature retrieved by LST3 and LST6 increases by about 0.3 units.e influence of atmospheric water vapor content on LST7 is nonlinear.With the increase of atmospheric water vapor, the LST retrieved by LST7 also increases, but the saturated phenomenon appears after 2.5. is result is consistent with the result of Jiménez et al. (2010).In Figure 5(g), the average land surface temperature parameters required in the MW algorithm are presented.It can be seen that, with the increase of the average land surface temperature by 1 unit, the inversion of the land surface temperature by MW decreases by 0.01 units.

Discussion
NASA's ACPC obtains the atmospheric parameters required for the LST inversion method based on the MODTRAN radiation transfer code.It is not possible to find atmospheric profiles in situ (radiosonde data, etc.) at any time and any place.erefore, even if this use (simulating atmospheric profiling information with ACPC) affects the accuracy of the method, it is clear from our results and the literature [8,10,14,18,20,[22][23][24][29][30][31][32][33]] that NASA's ACPC provides an excellent and effective simulation.e development of the LST inversion method has its source of error, as it involves parameterization steps for inverting the coefficients and estimating some initial parameters.erefore, it is necessary to compare the errors of various algorithms in the system.In general, the inversion results of the RTE algorithm and SC algorithm calculated by LSE1 and LSE2 have an excellent fitting effect with the measured data.e inversion results of MW are higher compared to the measured temperature.e inversion effect of SW is not particularly ideal, which may be due to the instability of radiation calibration at band 11 of TIRS.Hanqiu [42] used SC (LST3 and LST4) and SW (LST8 and LST9) algorithms to invert the LST in Fuzhou on August 4, 2013, and concluded that the SC result was better than the SW result.e inversion error of LST4 was −2.77 °C.However, since the atmospheric water vapor content parameter was more than 4 on that day, the error of LST3 is −17.44 °C, and the inversion errors of LST8 and LST9 are −13.61°C and −5.75 °C, respectively.Wang [41] compared the sensitivity analysis of the three algorithms, SC, MW, and SW, and found that SW had the lowest sensitivity to the input parameter error.Asia Siddiqui et al. [19] selected the data of two landscapes in January and April to study the heat island in Bangalore, India, based on the RTE method and found that the relationship between LST and NDVI was firstly positively correlated and then negatively correlated.e various LST inversion algorithms are most stable in winter, and the stability becomes worse in other seasons.Aliihsan Sekertekin and Stefania Bonafoni [40] used Landsat-5, Landsat-7, and Landsat-8 data to retrieve surface temperature from rural areas and found that RTE and SC presented similar results regardless of the season, while MW differed from RTE and SC for all seasons, especially in summer.
e parameter error had little influence on the results in the humid environment.MW and SC methods are more sensitive to the error of input parameters, especially under hot and humid conditions. is result is consistent with our inversion results in urban areas.On the one hand, with the increase of temperature, the properties of ground objects (specific heat capacity and specific emissivity) have distinct differences, leading to a significant difference in reflectivity and then resulting in a significant difference in the inversion of surface temperature.On the other hand, considering the errors of surface temperature inversion of similar surface features POI from inside and outside the city, the increase of temperature will have a particular influence on the parameters of the LST inversion algorithm, which may lead to errors in the inversion results.e stability of various inversion algorithms for temperature inversion inside and outside the city is tested.e results show that the standard deviations of inversion results for POIs of LST inversion algorithms in urban areas are smaller than those of inversion results outside the city.On the one hand, the main reason may be that the intensive human activities in cities have a particular influence on land surface temperature inversion parameters.On the other hand, the inversion results of the same reflectivity characteristics will be relatively concentrated when the same estimated parameters are used for calculation in a very narrow space inside the city.However, the space outside the city is relatively broad and belongs to the natural surface, so the estimated atmospheric parameters in different regions will be different, which may lead to the relatively scattered inversion results of land surface temperature.
It is essential to analyze the sensitivity/uncertainty of the LST inversion algorithm by considering all the input parameters.Sekertekin and Bonafoni's [43] simultaneous ground-based LST measurements were collected from Atmospheric Radiation Measurement (ARM) and Surface Radiation Budget Observing Network (SURFRAD) stations, located at different rural environments of the United States to examine the efficiency of different LST algorithms, using both daytime and nighttime Landsat-8 data and in situ measurements.Concerning the in situ sensitivity results, the effects on LST of the uncertainty of the downwelling and upwelling radiance were almost identical in daytime and nighttime.e accuracy of the LST retrieval methods for daytime Landsat-8 data varied between 2.17 K Root Mean Square Error (RMSE) and 5.47 K RMSE considering all LST methods and LSE models.MW with two different LSE models presented the best results for the daytime.Wang et al. [44] presented an improved Mono Window (IMW) algorithm for LST retrieval from the Landsat-8 TIRS band 10 data, and sensitivity analysis conducted for the IMW algorithm revealed that the possible error in estimating the required atmospheric water vapor content has the most significant impact on the probable LST estimation error.is result is similar to our inversion results of LST3, LST6, and LST7.Although from our results and literature we can see clearly that NASA ACPC provides a satisfying and effective simulation, in our precise and fuzzy provided coordinates, respectively, phenology, area type (urban and rural), and satellite altitude parameters such as atmospheric profile parameters calculation, calculation results by ACPC difference are not very big; however, when we use the estimated atmospheric profile parameters for further calculation, there will be error propagation.Sensitivity analysis of the LST inversion algorithm can help researchers to understand the source of error and reduce the error as much as possible.

Conclusions
In all seasons, the inversion results of the nine LST algorithms are in complete agreement with the graph trends of the measured data of all meteorological stations.In winter, the inversion results of the nine algorithms have the highest matching degree with the measured data of meteorological stations, and the RMSE and Bias values are relatively small.e second is autumn's LST1, LST2, LST4, and LST5 inversion results of the nearest meteorological station, with inversion error of ±5 °C; the following inversion result of LST3 is about 10 °C; the worst are the inversion results of LST6, LST7, LST8 and LST9, with inversion error of about 20 °C.In the low-temperature season, the error of the algorithm inversion results is less than that of the measured data.With the increase of temperature, the inversion results of each inversion algorithm are all higher compared to the measured data.e inversion errors of LST1, LST2, LST4, and LST5 are small, followed by LST3, LST6, LST7, LST8, and LST9.
According to the relevant results of the nine LST inversion results and the measured temperature, the fitting results in winter are all relatively ideal, with R 2 greater than 0.73.In other seasons, the fitting effect is general, and R 2 is about 0.5.In general, the inversion results of the RTE algorithm and SC algorithm calculated by LSE1 and LSE2 have an excellent fitting effect with the measured data, followed by the inversion results of the SC algorithm calculated by w.By analyzing the stability of various inversion algorithms inside and outside the city, it can be seen that the standard deviation of inversion results of each LST inversion algorithm for POI inside the city is less than that outside the city.
e influence of atmospheric transmittance on RTE, SC, and MW inversion algorithms is a logarithmic function.With the increase of atmospheric transmittance, the inversion surface temperature decreases gradually.
e influence of emissivity on all algorithms is linear.e inversion of LST decreases by 1 for every 0.1 increase of emissivity.e influence of NDVI on LST1, LST3, LST4, and LST8 is piecewise linear, and the influence on LST2, LST5, and LST9 is in the form of a quadratic function.When NDVI is equal to 0.4, the function has a minimum value.When NDVI is more than 0.05 and less than 0.4, the inversion land surface temperature decreases.When NDVI is more than 0.4, the inversion LST increases; this is the opposite of the previous results.e effect of uplink radiation on the SC (LST4 and LST5) algorithm is linear.With the increase of uplink radiation, the LST retrieved by the SC algorithm decreases linearly.e influence of upline radiation on the RTE (LST1 and LST2) algorithm is logarithmic.With the increase of upline radiation, the land surface temperature retrieved by the RTE algorithm decreases.e influence of downslope radiation on SC and RTE algorithms is linear.With the increase of downslope radiation by 1 unit, the inversion LST decreases by 0.2 units.e surface temperature retrieved by LST3 and LST6 increases by about 0.3 units for every 0.1 unit increase of atmospheric water vapor content.e influence of atmospheric water vapor content on LST7 is nonlinear.With the increase of atmospheric water vapor, the surface temperature retrieved by LST7 also increases, but, after 2.5 units, it shows saturation phenomenon.

Figure 1 :
Figure 1: Study areas.e remote sensing data of "LC08_L1TP_123032_20170507_20170515_01_T1," displayed in true color, is named 2017-05-07, with red dots representing the location of the meteorological station, green dots representing the points of interest outside the city, and yellow dots representing the points of interest inside the city.

Figure 3 :Figure 4 :
Figure 3: LST inversion results are fitted with the measured temperature.

Figure 5 :
Figure 5: Sensitivity analysis of parameters of various LST inversion algorithms.e vertical axis represents the normalized LST, and the horizontal axis represents the parameters of each variable.e units of parameters are given in Section 3.

Table 1 :
Basic information of meteorological stations in the study area.," where x stands for band number.A ρ is the reflectivity adjustment parameter of band λ, the header file of the data with the field "Reflectance_Add_Band_x � ."ρ λ ′ can be further corrected by equation (3) into atmospheric top reflectivity ρ λ : (Zhang et al. 2016)rical equations to estimate T a through air temperature (T o ), since it is an essential parameter of MW(Zhang et al. 2016).In this study, T a was computed for the mid-latitude summer and winter region and T o was obtained from the mean of all meteorological stations.
[21]4.SW Inversion Methods.Although the USGS does not encourage users to use band 11 calculations, some researchers have developed a Split Window algorithm for Landsat-8 that references other satellite two-channel algorithms.Jiménez-Muñoz et al.[21]proposed a revised version of SW for LST retrieval using Landsat TIRS data.e algorithm is simple to calculate.After the brightness temperature of band 10 and band 11 of TIR is obtained, it can be calculated by the following equation:

Table 2 :
e linear equations for the calculation of the effective mean atmospheric temperature (T a ) from the near-surface air temperature (T o ).

Table 3 :
Estimation of the associated LSE.

Table 4 :
a i and b i over different temperature ranges.

Table 5 :
Relationship between atmospheric transmittance and atmospheric water vapor content.

Table 6 :
Different LST inversion methods.LST9 of the inversion results of different LST inversion algorithms in different seasons.In addition, considering the difference of inversion results of the same reflectivity ground objects by different LST inversion algorithms, we designed two groups of experimental areas inside the city and outside the city.Each experimental area has four types of ground objects, respectively, and each type of ground object has 20 POIs.Finally, we analyzed the dependence of different LST inversion algorithms on parameters by equal step size.

Table 7 :
Root Mean Square Error (RMSE), average bias, slope, intercept, and R 2 of LST inversion and field measured temperature.