A Split Window Algorithm for Retrieving Land Surface Temperature from FY-3D MERSI-2 Data

: The thermal infrared (TIR) data from the Medium Resolution Spectral Imager II (MERSI-2) on the Chinese meteorological satellite FY-3D have high spatiotemporal resolution. Although the MERSI-2 land surface temperature (LST) products have good application prospects, there are some deviations in the TIR band radiance from MERSI-2. To accurately retrieve LSTs from MERSI-2, a method based on a cross-calibration model and split window (SW) algorithm is proposed. The method is divided into two parts: cross-calibration and LST retrieval. First, the MODTRAN program is used to simulate the radiation transfer process to obtain MERSI-2 and Moderate Resolution Imaging Spectroradiometer (MODIS) simulation data, establish a cross-calibration model, and then calculate the actual brightness temperature (BT) of the MERSI-2 image. Second, according to the characteristics of the near-infrared (NIR) bands, the atmospheric water vapor content (WVC) is retrieved, and the atmospheric transmittance is calculated. The land surface emissivity is estimated by the NDVI-based threshold method, which ensures that both parameters (transmittance and emissivity) can be acquired simultaneously. The validation shows the following: 1) The average accuracy of our algorithm is 0.42 K when using simulation data; 2) the relative error of our algorithm is 1.37 K when compared with the MODIS LST product (MYD11A1); 3) when compared with ground-measured data, the accuracy of our algorithm is 1.23 K. Sensitivity analysis shows that the SW algorithm is not sensitive to the two main parameters (WVC and emissivity), which also proves that the estimation of LST from MERSI-2 data is feasible. In general, our algorithm exhibits good accuracy and applicability, but it still requires further improvement.


Introduction
Temperature is one of the key physical parameters that characterize the energy transfer between various Earth layers. As an important index that reflects the state of the land surface, land surface temperature (LST) can provide information on factors such as the heat and mass fluxes between the land surface and the atmosphere, vegetation water stress, and soil moisture. Moreover, the acquisition of this information is indispensable in many research fields such as hydrology, environment, ecology,

Data Introduction
FY-3D is the fourth satellite in China's second-generation polar-orbiting meteorological satellite family. Launched in November 2017, FY-3D was placed in operation as the main satellite for China's low-altitude afternoon observations. MERSI-2 is one of the inherited instruments carried by FY-3D, which has been significantly upgraded and functionally integrated based on two imaging instruments (MERSI-1 and VIIRS) on FY-3. Finally, the performance of FY-3D has been significantly improved, with high-precision on-board calibration and lunar calibration capabilities. FY-3D can retrieve quantitative atmospheric, land, and ocean parameters, such as clouds, aerosols, water vapor, land surface characteristics, and ocean color parameters, with high precision. MERSI-2 has 25 channels with a scanning range of ±55.4 • and a sensor spectral coverage of 0.412-12.0 µm, which completely covers the band range of the MERSI-1 remote sensors on FY-3A/B/C. Specifically, the satellite has 15 visible and NIR bands (0.402-0.965 µm), four shortwave infrared bands (1.02/1.23-2.155 µm), two medium-wave infrared bands (3.71-4.1275 µm), one water vapor band (6.95-7.45 µm) and three longwave infrared bands (8. 4-12.5 µm). In terms of product applications, channels 1-7 are mainly used for land and cloud boundaries; channels 8-15 are used for ocean color, plankton and biogeochemical remote sensing; channels 16,17,18,22, and 23 are used for atmospheric water vapor remote sensing; channel 19 is a cloud band; and channels 20, 21, 24, and 25 are used for land, water, and cloud temperature sensing. Table 1 shows the continuity and evolution of the MERSI to MERSI-2 bands and the comparison of MERSI-2 with MODIS, MERSI, and VIIRS. Table 2 shows that the additional MERSI-2 bands include five infrared bands, one cirrus band (1.38 µm), one NIR water vapor band, one band at the wavelength of 7.2 µm, and two TIR bands with 250 m resolution. Generally, MERSI-2 is similar to MODIS and VIIRS in terms of band setting. In summary, MERSI-2 represents a functionally integrated and upgraded version of its first-generation sensors of the same track type and has been greatly improved in terms of data quality, number of bands, and resolution, further expanding the detection range and application fields.
Based on the SW algorithm characteristics, MERSI-2 Band 24 (10.3-11.3 µm) and Band 25 (11.5-12.5 µm) are selected for LST retrieval. Figure 1 shows the spectral response functions of the MERSI-2 TIR bands (Bands 24/25) and the MODIS TIR bands 31 (10.78-11.28 µm) and 32 (11.77-12.27 µm). As seen in Table 2 and Figure 1, TIR band ranges of MERSI-2 and MODIS are similar, and the center wavelengths are close, which indicates the rationality of using MODIS data as one of the cross-calibration and validation references.
include five infrared bands, one cirrus band (1.38 μm), one NIR water vapor band, one band at the wavelength of 7.2 μm, and two TIR bands with 250 m resolution. Generally, MERSI-2 is similar to MODIS and VIIRS in terms of band setting. In summary, MERSI-2 represents a functionally integrated and upgraded version of its first-generation sensors of the same track type and has been greatly improved in terms of data quality, number of bands, and resolution, further expanding the detection range and application fields. Based on the SW algorithm characteristics, MERSI-2 Band 24 (10.3-11.3 μm) and Band 25 (11.5-12.5 μm) are selected for LST retrieval. Figure 1 shows the spectral response functions of the MERSI-2 TIR bands (Bands 24/25) and the MODIS TIR bands 31 (10.78-11.28 μm) and 32 (11.77-12.27 μm). As seen in Table 2 and Figure 1, TIR band ranges of MERSI-2 and MODIS are similar, and the center wavelengths are close, which indicates the rationality of using MODIS data as one of the cross-calibration and validation references.

The Research Approach
The framework and technical methods of this research are described in Figure 2. To further explore the possibility of retrieving LST from MERSI-2 data and provide support for the development of an LST product algorithm, the main approach of this study is as follows: first, we propose a method for retrieving LST using the SW algorithm after cross-calibration of the BT and conduct an example application analysis. Second, in the cross-calibration stage, as MERSI-2 and MODIS data have similar imaging times and TIR band range settings (MERSI-2 Bands 24/25 and MODIS Bands 31/32), the MERSI-2 and MODIS simulation data are obtained by using the MODTRAN program to simulate the radiation transfer process, and the BT relationship model (BT Rel-Model) is established to convert the MODIS BT (simulated value) to MERSI-2. Next, according to the BT Rel-Model, the BT calibration model (BT Cal-Model) is established to correct the BT of the MERSI-2 image to the actual value and then calculate the cross-calibrated BT from MERSI-2. Third, in the LST retrieval stage, according to the characteristics of NIR bands, we can retrieve the water vapor content (WVC) and then calculate the atmospheric transmittance from the WVC in the image. The land surface emissivity is estimated by the NDVI-based threshold method so that the main parameters for retrieving the LST from the same image are determined. Finally, we perform a sensitivity analysis for the main parameters (WVC, atmospheric transmittance, and emissivity) in the SW algorithm to evaluate the applicability and promote the application of the algorithm.

The Research Approach
The framework and technical methods of this research are described in Figure 2. To further explore the possibility of retrieving LST from MERSI-2 data and provide support for the development of an LST product algorithm, the main approach of this study is as follows: first, we propose a method for retrieving LST using the SW algorithm after cross-calibration of the BT and conduct an example application analysis. Second, in the cross-calibration stage, as MERSI-2 and MODIS data have similar imaging times and TIR band range settings (MERSI-2 Bands 24/25 and MODIS Bands 31/32), the MERSI-2 and MODIS simulation data are obtained by using the MODTRAN program to simulate the radiation transfer process, and the BT relationship model (BT Rel-Model) is established to convert the MODIS BT (simulated value) to MERSI-2. Next, according to the BT Rel-Model, the BT calibration model (BT Cal-Model) is established to correct the BT of the MERSI-2 image to the actual value and then calculate the cross-calibrated BT from MERSI-2. Third, in the LST retrieval stage, according to the characteristics of NIR bands, we can retrieve the water vapor content (WVC) and then calculate the atmospheric transmittance from the WVC in the image. The land surface emissivity is estimated by the NDVI-based threshold method so that the main parameters for retrieving the LST from the same image are determined. Finally, we perform a sensitivity analysis for the main parameters (WVC, atmospheric transmittance, and emissivity) in the SW algorithm to evaluate the applicability and promote the application of the algorithm.

The Principle of LST Retrieval
During transmission through the atmosphere to the satellite sensor, land surface radiation is subject to many disturbances, such as those from the land surface type, soil moisture, near-surface

The Principle of LST Retrieval
During transmission through the atmosphere to the satellite sensor, land surface radiation is subject to many disturbances, such as those from the land surface type, soil moisture, near-surface air temperature, and WVC. Therefore, it is a prerequisite to obtain the LST to accurately quantify the mode of interaction between the LST and main interference factors. The derivation of the LST retrieval algorithm is based on the process of land surface heat radiation transfer, which can be quantified as shown in Equation (1) [27]: where B λ (T λ ) denotes the radiance at wavelength λ received by the sensors at the top of the atmosphere; B λ (T s ) denotes the radiation emitted by the ground objects and is often quantified using the Planck function (Equation (2)); ε λ is the land surface emissivity at wavelength λ, and (1 − ε λ ) is the reflectance. In addition, τ λ is the atmospheric transmittance at wavelength λ from the ground to the remote sensor; and L ↑ λ denotes the upwelling path radiances, while L ↓ λ denotes the downwelling path radiances.
where B λ (T s ) indicates the intensity of black radiation (W/m 2 µm); h indicates the Planck constant where T a represents the average temperature of upward radiance of atmosphere, and T ↓ a represents the average temperature of the downward radiance of atmosphere. Qin et al. (2001) [18,22] concluded that replacing T ↓ a with has no substantial effect on the calculation of the equation according to the results of an analysis and comparison. Thus, the land surface heat radiation transfer equation can usually be simplified to Equation (4): The parameters in Equation (4) can be further clarified from left to right: B λ (T λ ) is a known number, and the remaining parameters are unknowns; that is, the equation contains four unknowns.

The SW Algorithm Derivation
Many scholars have proven the feasibility of using the SW algorithm to retrieve LSTs [28,29]. In addition, according to different sensors and research areas, different methods are chosen to simplify and solve the parameters of the thermal radiation transfer equation, to propose different SW algorithms, and implement constant improvements. The general expression of the SW algorithm is as follows: where is a parameter; and represent the BT of channels m and n, respectively; and the units for T s , T m , and T n are absolute temperatures (K). Referring to Equation (5) and continuing to solve Equation (4), further simplification is needed. The traditional algorithm involves performing the Taylor expansion on the Planck function and taking one item for simplification. Many studies have proven that linear simplification of the Planck equation is also applicable [22]. Thus, the correlation analysis between the radiation intensity and temperature for MERSI-2 Bands 24 and 25 is carried out, and the results show that there is an approximate linear relationship between the radiation intensity and the temperature, as indicated by a scatter-plot analysis of the two channels between 273 and 322 K (0-49 • C) (see Figure 3).
The results (Equation (6)) of simplifying the Planck function are substituted into the thermal radiation transfer equation (Equation (4)), and we can obtain Equation (7) to calculate LST by the final simplification. where To obtain LST, it is necessary to obtain only two parameters: land surface emissivity ε λ and atmospheric transmittance τ λ . By simplifying the Planck function with the SW algorithm, we can reduce the initial need to obtain three complex unknown parameters to two unknown parameters, which are easier to obtain and greatly simplify the LST retrieval workload while ensuring high accuracy.

Practical Example Analysis
To perform an application analysis, we selected a MERSI-2 image from August 3, 2018 (UTC 4: 40) and a MODIS image from August 3, 2018 (UTC 4: 20). During the hot summer, the vegetation in the Bohai Sea area of China is prosperous and the ground recognition conditions are sufficient. The MERSI-2 image was radiometrically scaled, reflectivity calibrated, geometrically corrected, and cross-calibrated to obtain the effective BT, reflectivity, and other information. Finally, the SW algorithm was used to analyze the spatial differences in LST and heat in this area.

BT Cross-Calibration
Through our analysis, we found that the BT detected by the MERSI-2 instrument deviated greatly from the actual value. To determine the reasons, we can think from the following two aspects: on one hand, the MERSI-2 sensor may be offset; on the other hand, the TIR data from MERSI-2 need to be further improved via laboratory calibration and field calibrations, which also hinders the development of MERSI-2 TIR data into allow for use in more applications. If no cross-validation is conducted, it is nearly impossible to retrieve LSTs. Furthermore, MODIS LST products are widely used and have been calibrated by NASA for many years, providing an accurate reference object for LST retrieval. To further improve the LST retrieval accuracy, we propose a method using the SW algorithm after cross-calibration of the BT of MERSI-2 TIR bands (Bands 24/25). Accordingly, in the BT cross-calibration stage, two important models are established in this study: the BT relationship model (BT Rel-Model) and the BT calibration model (BT Cal-Model).

BT Relationship Model Establishment
MERSI-2 has imaging times that are similar to those of the MODIS data as well as similar TIR range settings (MERSI-2 Band 24 and MODIS Band 31, MERSI-2 Band 25 and MODIS Band 32). To improve the accuracy of radiation correction, we used the MODTRAN program to simulate the radiation transfer process to obtain MERSI-2 and MODIS simulation data in the mid-latitude region, and the BT Rel-Model was established between the BTs of MODIS and MERSI-2. As one of the key parameters of LST, land surface emissivity is mainly determined by the land surface type and the band intervals of the remote sensor. Thus, at a spatial resolution of 1 km, the pixels in the Bohai Sea area can be divided into soil, vegetation, and water.
More concretely, the first step is to obtain the BTs of the MERSI-2 and MODIS data in TIR bands. We used the MODTRAN to simulate the radiation transfer process and set the temperature (290-320 K, step length is 1 K), WVC (0.4-3.5 g/cm 2 , step length is 1 g/cm 2 ), and land surface emissivity (vegetation, water, and soil) values, for a total of 2976 groups of data. Furthermore, we used the PHP program to use MODTRAN to calculate the BTs of MERSI-2 and MODIS.
The land surface emissivities ε ω , ε ν , and ε s for water, vegetation, and soil, respectively, are usually based on the spectral reflectivity of the ground object, such as the reflectivity spectrum library of ground objects. Therefore, according to the experimental statistics of a variety of soil, vegetation, and water measured reflectivities (as shown in Figure 4) provided by the Jet Propulsion Laboratory (JPL: URL http://speclib.jpl.nasa.gov), we selected their average value to represent the land surface emissivities of various objects, where the input of the relevant parameters can be referred to as the technical roadmap ( Figure 2) and Table 3.
area can be divided into soil, vegetation, and water.
More concretely, the first step is to obtain the BTs of the MERSI-2 and MODIS data in TIR bands. We used the MODTRAN to simulate the radiation transfer process and set the temperature (290-320 K, step length is 1 K), WVC (0.4-3.5 g/ , step length is 1 g/ ), and land surface emissivity (vegetation, water, and soil) values, for a total of 2976 groups of data. Furthermore, we used the PHP program to use MODTRAN to calculate the BTs of MERSI-2 and MODIS.  The land surface emissivities , , and for water, vegetation, and soil, respectively, are usually based on the spectral reflectivity of the ground object, such as the reflectivity spectrum library of ground objects. Therefore, according to the experimental statistics of a variety of soil, vegetation, and water measured reflectivities (as shown in Figure 4) provided by the Jet Propulsion Laboratory (JPL: URL http://speclib.jpl.nasa.gov), we selected their average value to represent the land surface emissivities of various objects, where the input of the relevant parameters can be referred to as the technical roadmap ( Figure 2) and Table 3. Table 3. Sensor characteristics commonly used for retrieving LST.

Sensor
Channel Range Center Frequency  Sampling statistics and regression analysis of the MERSI-2 and MODIS BT simulation results (as shown in Figure 5) and the correlation coefficients are very high (both over 0.998).  The second step is to build the BT Rel-Model using the simulation data, MERSI-2 BT =f (MODIS), and the MERSI-2 BT can be obtained from the MODIS BT. Additionally, we further subdivided the model into Model 1 and Model 2, as shown in Equation (9) The simulation results show that with Model 1, we can approximately convert the BT of MODIS Band 31 to the BT of MERSI-2 Band 24; with Model 2, we can approximately convert the BT of MODIS Band 32 to the BT of MERSI-2 Band 24. In the same way, this law can also be applied to images. For example, with Model 1, we can convert the BT of the MODIS Band 31 image to the BT of the MERSI-2 Band 24 image.

BT Calibration Model Establishment
During the establishment of the BT Cal-Model, we convert the BT of the MODIS image to the actual BT of MERSI-2 according to the BT Cal-Model. Next, the statistical analysis of the BT of the MERSI-2 image and the BT converted from MODIS ( Figure 6) indicated that the BT of the MERSI-2 image is much higher than the real situation, and the correlations are very high (all greater than 0.93). The second step is to build the BT Rel-Model using the simulation data, MERSI-2 BT =f (MODIS), and the MERSI-2 BT can be obtained from the MODIS BT. Additionally, we further subdivided the model into Model 1 and Model 2, as shown in Equation (9): The simulation results show that with Model 1, we can approximately convert the BT of MODIS Band 31 to the BT of MERSI-2 Band 24; with Model 2, we can approximately convert the BT of MODIS Band 32 to the BT of MERSI-2 Band 24. In the same way, this law can also be applied to images. For example, with Model 1, we can convert the BT of the MODIS Band 31 image to the BT of the MERSI-2 Band 24 image.

BT Calibration Model Establishment
During the establishment of the BT Cal-Model, we convert the BT of the MODIS image to the actual BT of MERSI-2 according to the BT Cal-Model. Next, the statistical analysis of the BT of the MERSI-2 image and the BT converted from MODIS ( Figure 6) indicated that the BT of the MERSI-2 image is much higher than the real situation, and the correlations are very high (all greater than 0.93).

BT Calibration Model Establishment
During the establishment of the BT Cal-Model, we convert the BT of the MODIS image to the actual BT of MERSI-2 according to the BT Cal-Model. Next, the statistical analysis of the BT of the MERSI-2 image and the BT converted from MODIS ( Figure 6) indicated that the BT of the MERSI-2 image is much higher than the real situation, and the correlations are very high (all greater than 0.93).
(a) (b) Therefore, we established a validation model that converted the BT of the MERSI-2 image to the actual value, MERSI-2 BT (actual value) = f (MERSI-2 image data), which could be subdivided into Model 3 and Model 4, as shown in Equation (10): Currently, we can use Model 3 to correct the BT of a MERSI-2 Band 24 image to the cross-calibrated value based on its actual value. Similarly, we can use Model 4 to modify the BT of the MERSI-2 Band 25 image to the cross-calibrated value. As shown in Figures 7 and 8, by sampling the statistics, the relationships between the original, actual, and cross-calibrated values of the BT of the MERSI-2 image can be compared to obtain the following conclusion: the original value is higher than the actual value, but after correction by the cross-calibration method, the cross-calibrated value is basically consistent with the actual value. Therefore, we established a validation model that converted the BT of the MERSI-2 image to the actual value, MERSI-2 BT (actual value) =f (MERSI-2 image data), which could be subdivided into Model 3 and Model 4, as shown in Equation (10): Currently, we can use Model 3 to correct the BT of a MERSI-2 Band 24 image to the cross-calibrated value based on its actual value. Similarly, we can use Model 4 to modify the BT of the MERSI-2 Band 25 image to the cross-calibrated value. As shown in Figure 7 and Figure 8, by sampling the statistics, the relationships between the original, actual, and cross-calibrated values of the BT of the MERSI-2 image can be compared to obtain the following conclusion: the original value is higher than the actual value, but after correction by the cross-calibration method, the cross-calibrated value is basically consistent with the actual value.

Estimation of WVC and Atmospheric Transmittance
Atmospheric transmittance is an important factor in the radiation transfer process, and it is also one of the input parameters in the SW algorithm. For technical reasons, atmospheric transmittance is difficult to obtain in real time, so how to accurately and efficiently estimate atmospheric transmittance has become the focus of many scholars. However, since the WVC is the main factor affecting atmospheric transmittance, the atmospheric transmittance can be obtained by establishing

Estimation of WVC and Atmospheric Transmittance
Atmospheric transmittance is an important factor in the radiation transfer process, and it is also one of the input parameters in the SW algorithm. For technical reasons, atmospheric transmittance is difficult to obtain in real time, so how to accurately and efficiently estimate atmospheric transmittance has become the focus of many scholars. However, since the WVC is the main factor affecting atmospheric transmittance, the atmospheric transmittance can be obtained by establishing the relationship between the WVC and atmospheric transmittance. As a general rule, two main methods are used to establish this relationship. The first method is to estimate the WVC by using the measured data from a ground meteorological detection site to estimate the atmospheric transmittance. However, the land surface meteorological stations are unevenly distributed and affected by the topography, so the point-like water vapor data cannot accurately reflect the regional conditions. The second method is to use real-time atmospheric data to obtain the WVC and calculate atmospheric transmittance through radiation transfer process software such as MODTRAN, 6S, or LOWTRAN. With this method, it is also difficult to obtain the atmospheric profile data synchronized with the image to ensure the accuracy of the simulation results. Therefore, it is better to obtain WVC and atmospheric transmittance from the same remote sensing image. . Therefore, we can estimate the WVC directly from the same remote sensing image through NIR bands and then estimate the atmospheric transmittance to meet the accuracy and real-time requirements of parameter estimation.
After performing experiments, Kaufman et al. (1992) [30] found that it is feasible to use the ratio method to estimate the WVC, and many studies have proven the feasibility of this method [22,31]. The two ratio methods are calculated as follows: where C1 = 0.8, C2 = 0.2, and the weight assignment can be adjusted according to the situation; i is the water vapor absorption band; j and k are the atmospheric windows; and ρ is the reflectance of the band. The principle for the above two ratio methods is the same: in NIR bands, the ratio of the water vapor absorption band reflectance to the atmospheric window reflectance is approximated as the atmospheric transmittance. For the relationship between the atmospheric transmittance and the  (13) and (14)) by using the MODTRAN program.
In Equations (13) and (14), for complex surfaces, α = 0.02, β = 0.651, and the other land surface conditions are shown in the study of Kaufman et al. [30]. Thus far, we have estimated the WVC from NIR bands.

Estimating Atmospheric Transmittance from WVC
In the previous section, we calculated the WVC in NIR bands. From this information, we can simulate the atmospheric transmittance of TIR bands with the change in WVC through the MODTRAN program and then establish a model to estimate the atmospheric transmittance from the WVC. Therefore, the changes in the atmospheric transmittance of MERSI-2 Bands 24 and 25 with WVC in the summer midlatitude atmosphere are shown in Figure 9 below.  Figure 9 indicates that with the increase in WVC, the atmospheric transmittance shows a decreasing trend. After linear fitting and polynomial fitting, we selected the results of the cubic polynomial fitting (as shown in Equation (15), which had the highest correlation, to approximate the relationship between the atmospheric transmittance and the WVC in TIR bands: Thus far, it is possible to synchronously retrieve the WVC and the atmospheric transmittance in the same image, thereby greatly improving the accuracy and real-time performance of LST retrieval.

Estimating Land Surface Emissivity by Using the NDVI-Based Threshold Method
As one of the key parameters of LST, land surface emissivity is mainly determined by the land surface type and the band intervals of the remote sensor. The calculation of land surface emissivity based on NDVI values is a simple and effective method that has been applied and validated in many studies [31][32][33][34]. According to Mao's (2005) [22] idea for estimating the land surface emissivity from MODIS images with 1 km spatial resolution, the MERSI-2 image pixels in the Bohai Sea area can also be divided into soil, vegetation, and water. Therefore, Equation (16) is often used to roughly estimate  Figure 9 indicates that with the increase in WVC, the atmospheric transmittance shows a decreasing trend. After linear fitting and polynomial fitting, we selected the results of the cubic polynomial fitting (as shown in Equation (15), which had the highest correlation, to approximate the relationship between the atmospheric transmittance and the WVC in TIR bands: Thus far, it is possible to synchronously retrieve the WVC and the atmospheric transmittance in the same image, thereby greatly improving the accuracy and real-time performance of LST retrieval.

Estimating Land Surface Emissivity by Using the NDVI-Based Threshold Method
As one of the key parameters of LST, land surface emissivity is mainly determined by the land surface type and the band intervals of the remote sensor. The calculation of land surface emissivity based on NDVI values is a simple and effective method that has been applied and validated in many studies [31][32][33][34]. According to Mao's (2005) [22] idea for estimating the land surface emissivity from MODIS images with 1 km spatial resolution, the MERSI-2 image pixels in the Bohai Sea area can also be divided into soil, vegetation, and water. Therefore, Equation (16) is often used to roughly estimate the land surface emissivity of the MERSI-2 mixed pixels: where ε represents the land surface emissivity of the mixed pixels, ε ω , ε ν and ε s are the land surface emissions of water, vegetation, and soil, respectively; P ω , P ν , and (1 − P ω − P ν ) are the proportions of water, vegetation, and soil in the mixed pixels, respectively; and R ω , R ν , and R s are temperature ratios, and their calculation formula is R i = (T i /T) 4 . In addition, Qin (2004) [35] estimated the typical ground temperature ratio to be within the range of 5-45 • C. Therefore, we estimate the average temperature ratios of the three types of land surfaces as R ω = 0.99565, R ν = 0.99240 and R s = 1.00744.
There have been many studies on the NDVI-based threshold method, and the threshold values set for vegetation, water, and soil differ among studies, but the research principles are consistent. If there are obvious vegetation, water, and soil areas in the image (approximated as pure pixels), we can take the average value of NDVI in this area as the threshold, namely, NDVI v , NDVI ω and NDVI s . Therefore, when the NDVI of a mixed pixel is greater than the threshold, we take P i (i = w, v, s) = 1; otherwise, we calculate the ratio of the three features (vegetation, water, and soil) in the mixed pixel. In most studies [22,34,35], the NDVI value of dense and healthy vegetation is greater than 0.5, the NDVI of water is less than 0, and the NDVI value of bare soil is usually in the range of 0.03-0.2. In this study, the NDVI values of obvious water, vegetation, and bare soil areas are statistically analyzed; finally, we obtained NDVI ω = 0, NDVI v = 0.5 and NDVI s = 0.2.
In the mixed pixels, the ratio of water in the near-coast pixels can be estimated by P ω = QDN/DN1, where QDN is the DN of the pixels, and DN1 is the smallest DN of the land pixels. For land pixels, we can determine the vegetation coverage based on the NDVI value, as shown in Equations (17) and (18) [22]: Regarding the land surface emissivity values of water, vegetation, and soil (ε ω , ε ν and ε s ), according to the experimental statistics of the variety of soil, vegetation, and water measured reflectivity values (as shown in Figure 4), we select their average values as the land surface emissivity; that is, for MERSI-2 Band 24, we take ε ω = 0.992, ε ν = 0.9826, and ε s = 0.974, and for MERSI-2 Band 25, we take ε ω = 0.9862, ε ν = 0.987, and ε s = 0.979.

Retrieval Results and Analysis
To date, we have determined three unknown parameters involved in the SW algorithm solution: atmospheric transmittance, land surface emissivity, and BT (cross-calibrated). We substitute these equations into the SW algorithm (Equation (7)) to obtain the LST retrieval results in the Bohai Sea area, as shown in Figure 10

Retrieval Results and Analysis
To date, we have determined three unknown parameters involved in the SW algorithm solution: atmospheric transmittance, land surface emissivity, and BT (cross-calibrated). We substitute these equations into the SW algorithm (Equation (7)) to obtain the LST retrieval results in the Bohai Sea area, as shown in Figure 10. According to the LST retrieval results (Figure 10), we can obtain the following information: (1) The LST ranges from 18.98-56.26°C and the geometric information and temperature values are within the reasonable range for summer LST. (2) The urban heat island effect is prominent, and from According to the LST retrieval results (Figure 10), we can obtain the following information: (1) The LST ranges from 18.98-56.26 • C and the geometric information and temperature values are within the reasonable range for summer LST. (2) The urban heat island effect is prominent, and from north to south, we can clearly identify central cities such as Harbin, Changchun, Shenyang, Dalian, Beijing, and Tianjin. The temperatures of these urban centers are generally higher than those of the surrounding areas, showing the phenomenon of a temperature decrease from the center to the surrounding areas.
(3) We can also identify the northwestern region of Qinhuangdao (such as Chengde) and the dense vegetation area adjacent to the Japanese sea area because the temperatures in these areas are slightly lower than those in the surrounding areas, while the Tongliao and Chifeng areas are sparsely populated, and their temperatures are also the highest.

Validation
Validation of the algorithm is indispensable before the algorithm is promoted and applied. On the whole, there are three main methods for the validation of the algorithm for LST retrieval. The first method is to simulate the atmospheric radiation transfer process through atmospheric modeling software such as LOWTRAN or MODTRAN and evaluate the accuracy of the algorithm based on the obtained data. The second method is to compare the LST retrieval results with mature LST products and then evaluate the relative error. The third method is to use data measured from ground monitoring sites to evaluate the accuracy of the LST retrieval. However, these three evaluation methods have their own limitations: (1) The first method evaluates only the average accuracy and does not fully reflect the real ground thermal radiation transfer process during satellite transit. (2) The second method is very strict in terms of the imaging time of the sensors for comparison and is more suitable for comparison between data with similar imaging times. (3) For the third method, it is difficult to provide observation data synchronized with the satellite, and the point data do not easily reflect the spatial differentiation of the regional temperature. Therefore, to comprehensively evaluate the algorithm and retrieval results proposed in this study, the above three evaluation methods will be used to evaluate the accuracy.

Validation through a Standard Atmosphere Simulation
In this study, the principle of using the MODTRAN program to simulate the radiation transfer process is to input parameters such as spectral bandwidth, frequency, and center wavelength of MERSI-2 Bands 24/25 and set the LST, WVC, land surface type, and land surface emissivity values.
According to this method, the process of thermal radiation conduction at the land surface for the atmospheric conditions in the midlatitude region can be simulated; finally, we can obtain the following six types of output parameters: set LST, WVC, atmospheric transmittance, irradiance, surface type (emissivity of vegetation/water/soil), and BT. Further work involves recalculating the atmospheric transmittance and obtaining the actual value of the LST using the transmittance polynomial estimation model and the SW algorithm derived in this study. By comparing the actual value with the set LST, we can obtain the average accuracy of the SW algorithm, as shown in Figure 11. methods will be used to evaluate the accuracy.

Validation Through a Standard Atmosphere Simulation
In this study, the principle of using the MODTRAN program to simulate the radiation transfer process is to input parameters such as spectral bandwidth, frequency, and center wavelength of MERSI-2 Bands 24/25 and set the LST, WVC, land surface type, and land surface emissivity values. According to this method, the process of thermal radiation conduction at the land surface for the atmospheric conditions in the midlatitude region can be simulated; finally, we can obtain the following six types of output parameters: set LST, WVC, atmospheric transmittance, irradiance, surface type (emissivity of vegetation/water/soil), and BT. Further work involves recalculating the atmospheric transmittance and obtaining the actual value of the LST using the transmittance polynomial estimation model and the SW algorithm derived in this study. By comparing the actual value with the set LST, we can obtain the average accuracy of the SW algorithm, as shown in Figure  11. Figure 11. Technical route of validation through standard atmosphere simulation. Figure 11. Technical route of validation through standard atmosphere simulation.
In Figure 11, Tr is the transmittance. Through simulation, we obtained 2976 sets of data, which can reasonably evaluate the precision of the SW algorithm, and some of the process data and results are shown in Tables 4 and 5.  Generally, the average accuracy of our algorithm is 0.42 K when using simulation data, and the RMS is 0.19. In addition, the error distribution frequency and normal distribution curve for the SW algorithm are shown in Figure 12, which indicates that the SW algorithm proposed in this study has good accuracy and applicability. Generally, the average accuracy of our algorithm is 0.42 K when using simulation data, and the RMS is 0.19. In addition, the error distribution frequency and normal distribution curve for the SW algorithm are shown in Figure 12, which indicates that the SW algorithm proposed in this study has good accuracy and applicability.

Comparison with MODIS LST Products
Although it is difficult to evaluate LST retrieval results based on measured data, Wan et al. (2004) [36] performed extensive work and verified the MODIS LST products with measured data. The results show that the accuracy of the MODIS LST products in the experimental area is within 1 °C. Therefore, we select the MODIS LST products (MYD11A1, as shown in Figure 13) provided by NASA and the results of the MERSI-2 retrieval results using the SW algorithm for relative error analysis.

Comparison with MODIS LST Products
Although it is difficult to evaluate LST retrieval results based on measured data, Wan et al. (2004) [36] performed extensive work and verified the MODIS LST products with measured data. The results show that the accuracy of the MODIS LST products in the experimental area is within 1 • C. Therefore, we select the MODIS LST products (MYD11A1, as shown in Figure 13) provided by NASA and the results of the MERSI-2 retrieval results using the SW algorithm for relative error analysis. Remote Sens. 2017, 9,   As shown in Figure 13, the values from MODIS LST product (MYD11A1) ranges from 13.48-54.58 °C; that is, the retrieval result of MERSI-2 is basically consistent with the spatial heterogeneity of the MODIS LST product at high, medium, and low temperatures. Through our analysis of the LST profile under clear sky conditions in the images (as shown in Figure 14), the LST trends of MERSI-2 and the MODIS LST product are basically the same, while the average temperature error is 1.37 K, and the RMS is 1.17. The imaging times of MODIS and MERSI-2 are close but not the same, and the MODIS LST products have some errors compared with the real LST. Therefore, when we compare the retrieved LST of MERSI-2 with MODIS LST products, the final results are only the spatial analysis of the LSTs and relative errors and are not absolutely precise.

Comparison with Ground-Measured Data
It is not ideal to evaluate the spatial distribution of LST over a large scale by using point-like measured data. On one hand, it is difficult to obtain ground data that is consistent with the resolution of image pixels when satellites are in transit, and on the other hand, it is difficult to obtain parameters such as real-time atmospheric profile data. Therefore, we consider this method only for As shown in Figure 13, the values from MODIS LST product (MYD11A1) ranges from 13.48-54.58 • C; that is, the retrieval result of MERSI-2 is basically consistent with the spatial heterogeneity of the MODIS LST product at high, medium, and low temperatures. Through our analysis of the LST profile under clear sky conditions in the images (as shown in Figure 14), the LST trends of MERSI-2 and the MODIS LST product are basically the same, while the average temperature error is 1.37 K, and the RMS is 1.17. The imaging times of MODIS and MERSI-2 are close but not the same, and the MODIS LST products have some errors compared with the real LST. Therefore, when we compare the retrieved LST of MERSI-2 with MODIS LST products, the final results are only the spatial analysis of the LSTs and relative errors and are not absolutely precise. As shown in Figure 13, the values from MODIS LST product (MYD11A1) ranges from 13.48-54.58 °C; that is, the retrieval result of MERSI-2 is basically consistent with the spatial heterogeneity of the MODIS LST product at high, medium, and low temperatures. Through our analysis of the LST profile under clear sky conditions in the images (as shown in Figure 14), the LST trends of MERSI-2 and the MODIS LST product are basically the same, while the average temperature error is 1.37 K, and the RMS is 1.17. The imaging times of MODIS and MERSI-2 are close but not the same, and the MODIS LST products have some errors compared with the real LST. Therefore, when we compare the retrieved LST of MERSI-2 with MODIS LST products, the final results are only the spatial analysis of the LSTs and relative errors and are not absolutely precise.

Comparison with Ground-Measured Data
It is not ideal to evaluate the spatial distribution of LST over a large scale by using point-like measured data. On one hand, it is difficult to obtain ground data that is consistent with the resolution of image pixels when satellites are in transit, and on the other hand, it is difficult to obtain parameters such as real-time atmospheric profile data. Therefore, we consider this method only for

Comparison with Ground-Measured Data
It is not ideal to evaluate the spatial distribution of LST over a large scale by using point-like measured data. On one hand, it is difficult to obtain ground data that is consistent with the resolution of image pixels when satellites are in transit, and on the other hand, it is difficult to obtain parameters such as real-time atmospheric profile data. Therefore, we consider this method only for evaluating the results and algorithm accuracy of LST retrieval. The ground monitoring station data were obtained from the China Meteorological Administration. After quality control, 11 representative stations with single ground land cover types and flat terrain, which had imaging times close to those of the MERSI-2 data, were selected to evaluate the accuracy of the proposed algorithm. As shown in Figure 15, the evaluation results show that the accuracy of the SW algorithm is 1.23 K, and the RMS is 0.46. In general, the accuracy of the algorithm is relatively high, but the algorithm still requires further improvement in subsequent research.
Remote Sens. 2017, 9, x FOR PEER REVIEW 19 of 25 evaluating the results and algorithm accuracy of LST retrieval. The ground monitoring station data were obtained from the China Meteorological Administration. After quality control, 11 representative stations with single ground land cover types and flat terrain, which had imaging times close to those of the MERSI-2 data, were selected to evaluate the accuracy of the proposed algorithm. As shown in Figure 15, the evaluation results show that the accuracy of the SW algorithm is 1.23 K, and the RMS is 0.46. In general, the accuracy of the algorithm is relatively high, but the algorithm still requires further improvement in subsequent research. Figure 15. SW algorithm-retrieved LST (°C) compared with ground-measured data.

Discussion
Parameter sensitivity analysis is very important in algorithm development and practical applications. In the SW algorithm, emissivity and transmittance are the main influencing factors, and WVC is the factor that influences atmospheric transmittance. Therefore, the sensitivity analysis on the three main influencing parameters, namely, WVC, atmospheric transmittance, and emissivity, was carried out. Equation (19) can approximately evaluate the effect of parameter errors:

Sensitivity Analysis of WVC
In the SW algorithm, the atmospheric transmittance is estimated from WVC, so the sensitivity of the algorithm to the WVC determines the applicability of the algorithm. In this study, the MODTRAN program is used to simulate the radiation transfer process (see Figure 11) to obtain standard atmospheric simulation data for analyzing the sensitivity of WVC. Table 7 shows the partial LST inversion results obtained by changing the WVC error (−90-100%). From Table 6, we can also obtain the change in transmittance with the variation in WVC errors (−90-100%), where the LST error represents the difference between the retrieved LST before and after the change, and the RMS represents the root mean square.

Discussion
Parameter sensitivity analysis is very important in algorithm development and practical applications. In the SW algorithm, emissivity and transmittance are the main influencing factors, and WVC is the factor that influences atmospheric transmittance. Therefore, the sensitivity analysis on the three main influencing parameters, namely, WVC, atmospheric transmittance, and emissivity, was carried out. Equation (19) can approximately evaluate the effect of parameter errors:

Sensitivity Analysis of WVC
In the SW algorithm, the atmospheric transmittance is estimated from WVC, so the sensitivity of the algorithm to the WVC determines the applicability of the algorithm. In this study, the MODTRAN program is used to simulate the radiation transfer process (see Figure 11) to obtain standard atmospheric simulation data for analyzing the sensitivity of WVC. Table 6 shows the partial LST inversion results obtained by changing the WVC error (−90-100%). From Table 6, we can also obtain the change in transmittance with the variation in WVC errors (−90-100%), where the LST error represents the difference between the retrieved LST before and after the change, and the RMS represents the root mean square. From Table 6 and Figure 16, we can obtain the following information: (1) When there is no change in WVC (0%), the transmittances of MERSI-2 Bands 24 and 25 are 0.85 and 0.77, respectively, while the average LST error is 0.38, and the RMS is 0.13. (2) When the WVC error is 60%, the average transmittance errors of MERSI-2 Bands 24 and 25 are 0.11 and 0.15, respectively, and the accuracy of the LST retrieval is the highest (LST error is 0.36 K, and the RMS is 0.14). (3) When the WVC error is between −20% and 100%, the average transmittance error of MERSI-2 Bands 24 and 25 is between 0 and 0.2, while the mean LST error is very small (all less than 0.5 K) and the trend is stable; that is, the algorithm is not sensitive to the WVC.  Table 6 and Figure 16, we can obtain the following information: 1) When there is no change in WVC (0%), the transmittances of MERSI-2 Bands 24 and 25 are 0.85 and 0.77, respectively, while the average LST error is 0.38, and the RMS is 0.13. 2) When the WVC error is 60%, the average transmittance errors of MERSI-2 Bands 24 and 25 are 0.11 and 0.15, respectively, and the accuracy of the LST retrieval is the highest (LST error is 0.36 K, and the RMS is 0.14). 3) When the WVC error is between −20% and 100%, the average transmittance error of MERSI-2 Bands 24 and 25 is between 0 and 0.2, while the mean LST error is very small (all less than 0.5 K) and the trend is stable; that is, the algorithm is not sensitive to the WVC.  Kaufman and Gao (1992) [30] retrieved the WVC from MODIS under clear weather conditions with a maximum inversion error between −13% and +13%. Considering that MERSI-2 and MODIS have similar TIR band settings, the sensitivity analysis of the WVC for MERSI-2 can be analyzed with reference to the sensitivity analysis of MODIS. Finally, we obtain the following results: the algorithm proposed in this study is not sensitive to WVC, which proves that it is feasible to use this algorithm to estimate the transmittance from WVC, and the practicability of the SW algorithm can be further improved.

Sensitivity Analysis of Atmospheric Transmittance
To further analyze the sensitivity of atmospheric transmittance, we conducted a separate study on atmospheric transmittance, as shown in Figure 17. When we simultaneously change the transmittance error (−0.1-0.12) of MERSI-2 Bands 24 and 25 at the same ratio, the LST error ranges from 0.39-1.3 K. When the transmittance error ranges from −0.025-0.02, the LST error is less than 0.5 K. When the transmittance error of MERSI-2 Band 24 remains unchanged, and the transmittance error (−0.01-0.02) of MERSI-2 Band 25 is changed, the LST error ranges from 0.39-0.92 K. Similarly, when we change only the transmittance error (−0.01-0.05) of MERSI-2 Band 24, the corresponding LST error ranges from 0.40-1.02 K. In general, the algorithm is sensitive to the atmospheric transmittance relative to the sensitivity of WVC. The algorithm is not sensitive to WVC, which proves that it is feasible to use the WVC to estimate the atmospheric transmittance and subsequently retrieve the LST. The introduction of WVC will greatly improve the applicability of the algorithm, and in future research, we can use prior knowledge of the sensitivity analysis to further improve the accuracy of the LST retrieval. Kaufman and Gao (1992) [30] retrieved the WVC from MODIS under clear weather conditions with a maximum inversion error between −13% and +13%. Considering that MERSI-2 and MODIS have similar TIR band settings, the sensitivity analysis of the WVC for MERSI-2 can be analyzed with reference to the sensitivity analysis of MODIS. Finally, we obtain the following results: the algorithm proposed in this study is not sensitive to WVC, which proves that it is feasible to use this algorithm to estimate the transmittance from WVC, and the practicability of the SW algorithm can be further improved.

Sensitivity Analysis of Atmospheric Transmittance
To further analyze the sensitivity of atmospheric transmittance, we conducted a separate study on atmospheric transmittance, as shown in Figure 17. When we simultaneously change the transmittance error (−0.1-0.12) of MERSI-2 Bands 24 and 25 at the same ratio, the LST error ranges from 0.39-1.3 K. When the transmittance error ranges from −0.025-0.02, the LST error is less than 0.5 K. When the transmittance error of MERSI-2 Band 24 remains unchanged, and the transmittance error (−0.01-0.02) of MERSI-2 Band 25 is changed, the LST error ranges from 0.39-0.92 K. Similarly, when we change only the transmittance error (−0.01-0.05) of MERSI-2 Band 24, the corresponding LST error ranges from 0.40-1.02 K. In general, the algorithm is sensitive to the atmospheric transmittance relative to the sensitivity of WVC. The algorithm is not sensitive to WVC, which proves that it is feasible to use the WVC to estimate the atmospheric transmittance and subsequently retrieve the LST. The introduction of WVC will greatly improve the applicability of the algorithm, and in future research, we can use prior knowledge of the sensitivity analysis to further improve the accuracy of the LST retrieval.

Sensitivity Analysis of Land Surface Emissivity
The change in emissivity is very important for the retrieval accuracy of the algorithm, and the sensitivity of our algorithm to emissivity was also analyzed, as shown in Figure 18. When we simultaneously change the emissivity error (−0.023-0.008) of MERSI-2 Bands 24 and 25, the corresponding LST error ranges from 0.17-1.09 K; when the emissivity error ranges from −0.013-0, the LST error is less than 0.5 K. When we keep the emissivity error of MERSI-2 Band 24 unchanged, only the transmittance error (−0.01-0.013) of MERSI-2 Band 25 changes, and the corresponding LST error ranges from 0.15-1.33 K. Similarly, when we change only the transmittance error (−0.01-0.008) of MERSI-2 Band 24, the corresponding LST error ranges from 0.17-1.65 K. Through statistical analysis, we found that the emissivities of most features of MERSI-2 Bands 24 and 25 (10.30-11.30

Sensitivity Analysis of Land Surface Emissivity
The change in emissivity is very important for the retrieval accuracy of the algorithm, and the sensitivity of our algorithm to emissivity was also analyzed, as shown in Figure 18. When we simultaneously change the emissivity error (−0.023-0.008) of MERSI-2 Bands 24 and 25, the corresponding LST error ranges from 0.17-1.09 K; when the emissivity error ranges from −0.013-0, the LST error is less than 0.5 K. When we keep the emissivity error of MERSI-2 Band 24 unchanged, only the transmittance error (−0.01-0.013) of MERSI-2 Band 25 changes, and the corresponding LST error ranges from 0.15-1.33 K. Similarly, when we change only the transmittance error (−0.01-0.008) of MERSI-2 Band 24, the corresponding LST error ranges from 0.17-1.65 K. Through statistical analysis, we found that the emissivities of most features of MERSI-2 Bands 24 and 25 (10.30-11.30 µm and 11.50-12.50 µm) are higher than 0.97, and the degree of change is very small. Therefore, the algorithm is proven to not be sensitive to land surface emissivity. μm and 11.50-12.50 μm) are higher than 0.97, and the degree of change is very small. Therefore, the algorithm is proven to not be sensitive to land surface emissivity. The above analysis indicates that when the WVC and land surface emission error are within a certain range, the change in the surface temperature inversion error is within the acceptable range, which shows that the algorithm proposed in this study has good applicability. Furthermore, Figures  16 and 18 show that the land surface emissivity is more sensitive to the estimation of the LST by the SW algorithm than to the WVC. Based on the dynamic changes in NDVI, we determined the land surface emission rate of MERSI-2 data in real time, overcoming the previous limitation that the emission rate did not change over time when land surface classification data were used. Moreover, we consider that at the spatial resolution of 1 km, the visible light bands can be used to divide the features in the image into soil, vegetation, and water. The NDVI is used to determine the proportions of the major land surface types in the mixed pixels. Thus, the land surface emissivity values corresponding to the pure pixels in the image and the mixed pixels are calculated, and finally, the estimation of the land emissivity is accurate at the subpixel level. In fact, in addition to the type of land surface, the land surface emissivity is also affected by land surface fluctuations (especially in mountainous areas) and is related to the geometry of different vegetation types and the location of vegetation in the pixels. Therefore, in future research, we need to pay more attention to the retrieval of the LST of heterogeneous pixel components to continuously improve the retrieval accuracy of LST.

Conclusions
Some deviations exist in TIR band radiance of MERSI-2; therefore, a cross-calibration model was built in this study that fully utilizes MODIS data. Without cross-validation, LST retrievals are nearly impossible. Therefore, a method for retrieving LST using the SW algorithm was proposed after BT cross-calibration, with the aim of providing support for the development of an LST product algorithm. Three methods were used to evaluate the accuracies of the retrieval results and algorithms. Moreover, we performed a detailed analysis of the main parameters that influence the algorithm, which are the WVC, transmittance, and emissivity, and finally, we proved that the algorithm has good applicability. The main conclusions are as follows: (1) Cross-validation is of great significance for the study of LST retrievals. After cross-calibration of the MERSI-2 BT, the scaled value was basically consistent with the actual value, which was not the case for the uncalibrated value. To analyze a practical example, we used the cross-calibrated BT to retrieve the LST of the Bohai Sea area.
(2) In the LST retrieval stage, according to the characteristics of the NIR bands, we can retrieve the WVC and then calculate the atmospheric transmittance from the WVC of each pixel in the image. The above analysis indicates that when the WVC and land surface emission error are within a certain range, the change in the surface temperature inversion error is within the acceptable range, which shows that the algorithm proposed in this study has good applicability. Furthermore, Figures 16  and 18 show that the land surface emissivity is more sensitive to the estimation of the LST by the SW algorithm than to the WVC. Based on the dynamic changes in NDVI, we determined the land surface emission rate of MERSI-2 data in real time, overcoming the previous limitation that the emission rate did not change over time when land surface classification data were used. Moreover, we consider that at the spatial resolution of 1 km, the visible light bands can be used to divide the features in the image into soil, vegetation, and water. The NDVI is used to determine the proportions of the major land surface types in the mixed pixels. Thus, the land surface emissivity values corresponding to the pure pixels in the image and the mixed pixels are calculated, and finally, the estimation of the land emissivity is accurate at the subpixel level. In fact, in addition to the type of land surface, the land surface emissivity is also affected by land surface fluctuations (especially in mountainous areas) and is related to the geometry of different vegetation types and the location of vegetation in the pixels. Therefore, in future research, we need to pay more attention to the retrieval of the LST of heterogeneous pixel components to continuously improve the retrieval accuracy of LST.

Conclusions
Some deviations exist in TIR band radiance of MERSI-2; therefore, a cross-calibration model was built in this study that fully utilizes MODIS data. Without cross-validation, LST retrievals are nearly impossible. Therefore, a method for retrieving LST using the SW algorithm was proposed after BT cross-calibration, with the aim of providing support for the development of an LST product algorithm. Three methods were used to evaluate the accuracies of the retrieval results and algorithms. Moreover, we performed a detailed analysis of the main parameters that influence the algorithm, which are the WVC, transmittance, and emissivity, and finally, we proved that the algorithm has good applicability. The main conclusions are as follows: (1) Cross-validation is of great significance for the study of LST retrievals. After cross-calibration of the MERSI-2 BT, the scaled value was basically consistent with the actual value, which was not the case for the uncalibrated value. To analyze a practical example, we used the cross-calibrated BT to retrieve the LST of the Bohai Sea area.
(2) In the LST retrieval stage, according to the characteristics of the NIR bands, we can retrieve the WVC and then calculate the atmospheric transmittance from the WVC of each pixel in the image. The land surface emissivity is estimated by the NDVI-based threshold method so that the main parameters for retrieving LST from the same image are acquired. Taking the Bohai Sea area as the research area, we used three methods to validate our algorithm. The results show that the average accuracy of the algorithm is 0.42 K when standard atmospheric simulation data are used, and the RMS is 0.19. When our retrieval results are compared with the MODIS LST product, the LST trends of MERSI-2 and MODIS are basically the same, while the average LST error is 1.37 K, and the RMS is 1.17. When our retrieval results are compared with ground-measured data, the evaluation results show that the accuracy of the SW algorithm is 1.23 K, and the RMS is 0.46. In general, the accuracy of the algorithm is relatively high, but the algorithm requires further improvement in subsequent research.
(3) After sensitivity analysis of the main parameters that affect the accuracy of the algorithm, namely, WVC, transmittance, and emissivity, the results show that the algorithm is not sensitive to the WVC and emissivity but is sensitive to transmittance. However, the feasibility of estimating transmittance from WVC is proven. Within a the reasonable interval (−13-13%) of the WVC and when the WVC error is between −20% and 100%, the average transmittance error of MERSI-2 Bands 24 and 25 is between 0 and 0.2, and at this time, the average LST errors are small (all less than 0.5 K), and the trends are stable; that is, the algorithm is not sensitive to WVC. Within a reasonable interval of emissivity change retrievals (generally higher than 0.97), when we change the emissivity error (−0.023-0.008) of MERSI-2 Bands 24 and 25, we can obtain an LST error range of 0.17-1.09 K, which also shows that our algorithm is not sensitive to emissivity. Compared with the interval of error variations for the WVC, our algorithm is slightly sensitive to atmospheric transmittance because when the MERSI-2 Bands 24 and 25 transmittance errors (−0.1-0.12) are changed simultaneously, the LST errors are 0.39-1.3 K. In general, the algorithm for estimating transmittance from atmospheric WVC to retrieve the LST has good applicability.