Orbital Lifetime (2008–2017) Radiometric Calibration and Evaluation of the HJ-1B IRS Thermal Infrared Band

: The infrared sensor (IRS) is a payload on the HJ-1B satellite and includes a thermal infrared band (B08). In order to obtain radiometric calibration coe ﬃ cients and evaluate annual change, this study performed an analysis covering its 10-year orbital lifetime (2008–2017). The cross-calibration of IRS B08 with MODIS was performed using near-simultaneous images over Lake Qinghai, China. The results reveal that the radiometric response of IRS B08 notably changed during its orbital lifetime from year-to-year. The o ﬀ sets ﬂuctuated more than the gain. The top-of-atmosphere (TOA) radiance obtained by calibration coe ﬃ cients in this study was generally in agreement with those obtained by onboard calibrator, within an error range of ± 4.00% from 2008 to 2012. The percent di ﬀ erence compared with ﬁeld validation was within 1.63%. The di ﬀ erence between IRS and MODIS radiance over ﬁeld validation sites was within ± 5%. Approximately a 1% di ﬀ erence occurred between the TOA temperature of IRS and MODIS. The radiometric response of IRS B08 continuously decreased from 2008 to 2013, whereas it ﬂuctuated from 2014 to 2017. Moreover, the DN ﬂuctuated more when the at-aperture radiance was low, although it was more stable at higher radiance. These indicate that the cross-calibration coefficients calculated in this


Introduction
Research utilizing derived thermal infrared data has become a hot topic in recent decades. The value of thermal data archives will continue to grow since they can be applied to many areas, such as the monitoring of evapotranspiration [1,2], the estimation of soil moisture [3], the detection of the geothermal system [4,5] and mineral mapping [6,7]. In order to take proper advantage of thermal infrared data, a critical step is radiometric calibration, which converts observed digital numbers (DNs) to physical radiance quantities [8]. One of the virtual purposes of this process is to provide calibration coefficients to consistently convert DNs into top-of-atmosphere (TOA) radiances. The converted radiance that explicitly embodies the reflective or emissive properties from the terrestrial surface is usually taken as the most important preliminary input in quantitative applications.
In general, there are four main methods used to obtain calibration coefficients: preflight calibration, in-flight calibration, vicarious calibration and cross-calibration. Since the in-orbit environment differs from the simulated environment in the laboratory, post-flight calibration is required even if the preflight

Reference Data
The cross-calibration method requires relatively high calibration accuracy of the reference sensor. The MODIS is an important sensor on the Terra that observes the Earth's surface every 1-2 days. The MODIS instrument has 36 bands, and its thermal infrared (TIR) range is 3-15μm. MODIS band 31 (10.78-11.28 μm) and band 32 (11.77-12.27 μm) overlapped with IRS band 08 (10.50-12.50 μm); a comparison of their channel response functions is presented in Figure 2. In addition, the MODIS TIR bands are perfectly calibrated through the onboard blackbody, and the uncertainty of absolute radiometric accuracy for bands 31 and 32 is 2% [33,34]. Since both are morning satellites the pass time difference between the Terra and the HJ-1B over Lake Qinghai is about 1 h-and sometimes only a few minutes. Therefore, this study selected MODIS B31 and B32 as the reference channels for the IRS B08 cross-calibration. Table 1 compares the thermal band features of both instruments.

Reference Data
The cross-calibration method requires relatively high calibration accuracy of the reference sensor. The MODIS is an important sensor on the Terra that observes the Earth's surface every 1-2 days. The MODIS instrument has 36 bands, and its thermal infrared (TIR) range is 3-15µm. MODIS band 31 (10.78-11.28 µm) and band 32 (11.77-12.27 µm) overlapped with IRS band 08 (10.50-12.50 µm); a comparison of their channel response functions is presented in Figure 2. In addition, the MODIS TIR bands are perfectly calibrated through the onboard blackbody, and the uncertainty of absolute radiometric accuracy for bands 31 and 32 is 2% [33,34]. Since both are morning satellites the pass time difference between the Terra and the HJ-1B over Lake Qinghai is about 1 h-and sometimes only a few minutes. Therefore, this study selected MODIS B31 and B32 as the reference channels for the IRS B08 cross-calibration. Table 1 compares the thermal band features of both instruments.   Near-simultaneous IRS and MODIS images over Lake Qinghai were selected and matched as one-to-one experimental data pairs. Excluding images dominated by cloud or ice, 208 image pairs were chosen. Based on the accuracy distribution characteristics of the twin-channel difference model used in this study (introduced in the section below), 56 pairs with MODIS-view zenith angles >50° were discarded. The remaining 152 pairs were then used to perform the radiometric calibration. The images pairs were further subdivided into three groups based on their corresponding MODIS-view zenith angle, i.e., within 30°, within 40° and within 50°. The annual effective data volumes are graphed in Figure 3a

Image Pairs Matching and Statistics
Near-simultaneous IRS and MODIS images over Lake Qinghai were selected and matched as one-to-one experimental data pairs. Excluding images dominated by cloud or ice, 208 image pairs were chosen. Based on the accuracy distribution characteristics of the twin-channel difference model used in this study (introduced in the section below), 56 pairs with MODIS-view zenith angles >50 • were discarded. The remaining 152 pairs were then used to perform the radiometric calibration. The images pairs were further subdivided into three groups based on their corresponding MODIS-view zenith angle, i.e., within 30 • , within 40 • and within 50 • . The annual effective data volumes are graphed in Figure 3a and the distribution of the MODIS-view zenith angle is plotted in Figure 3b.

Image Preprocessing
Image preprocessing includes fine geometric correction and resampling. The Landsat thematic mapper (TM) image orthorectification method has been successfully employed. Therefore, the IRS and MODIS images were geometrically corrected based on the fine geometric correction data of the Landsat TM (LT05_L1TP_133034_20070822). Since the spatial resolution of IRS B08 is 300 m but are 1 km for MODIS B31 and B32 bands, IRS B08 data were resampled to MODIS 1 km resolution using cubic convolution.
A uniform and stable subarea on the lake is the preferable region of interest (ROI) for performing cross-calibration. During the selection of the ROI, the lake surface was gridded into approximate 20 km × 20 km squares [15], as illustrated in Figure 4. The squares containing both water and land were discarded in order to ensure feature uniformity in the target area. Each square labeled with row and

Image Preprocessing
Image preprocessing includes fine geometric correction and resampling. The Landsat thematic mapper (TM) image orthorectification method has been successfully employed. Therefore, the IRS and MODIS images were geometrically corrected based on the fine geometric correction data of the Landsat TM (LT05_L1TP_133034_20070822). Since the spatial resolution of IRS B08 is 300 m but are 1 km for MODIS B31 and B32 bands, IRS B08 data were resampled to MODIS 1 km resolution using cubic convolution.

of 15
A uniform and stable subarea on the lake is the preferable region of interest (ROI) for performing cross-calibration. During the selection of the ROI, the lake surface was gridded into approximate 20 km × 20 km squares [15], as illustrated in Figure 4. The squares containing both water and land were discarded in order to ensure feature uniformity in the target area. Each square labeled with row and column numbers in Figure 4 was iterated in order to mask and extract the radiance from MODIS B31 and B32, after which the radiance variance was calculated for each. The blue grid (6,4) in Figure 4 corresponds to the smallest variance, indicating that the average thermal radiative characteristics in this grid are more uniform and stable than those in other grids. Therefore, this grid square was selected as the final ROI used to perform radiometric cross-calibration in this study.

Image Preprocessing
Image preprocessing includes fine geometric correction and resampling. The Landsat thematic mapper (TM) image orthorectification method has been successfully employed. Therefore, the IRS and MODIS images were geometrically corrected based on the fine geometric correction data of the Landsat TM (LT05_L1TP_133034_20070822). Since the spatial resolution of IRS B08 is 300 m but are 1 km for MODIS B31 and B32 bands, IRS B08 data were resampled to MODIS 1 km resolution using cubic convolution.
A uniform and stable subarea on the lake is the preferable region of interest (ROI) for performing cross-calibration. During the selection of the ROI, the lake surface was gridded into approximate 20 km × 20 km squares [15], as illustrated in Figure 4. The squares containing both water and land were discarded in order to ensure feature uniformity in the target area. Each square labeled with row and column numbers in Figure 4 was iterated in order to mask and extract the radiance from MODIS B31 and B32, after which the radiance variance was calculated for each. The blue grid (6,4) in Figure 4 corresponds to the smallest variance, indicating that the average thermal radiative characteristics in this grid are more uniform and stable than those in other grids. Therefore, this grid square was selected as the final ROI used to perform radiometric cross-calibration in this study.

TOA Radiance Calibration and Coefficients Regression
The twin-channel difference model, which is a simplified physical model, was established based on radiative transfer model [16]. This model can avoid the effects of atmospheric temperature and density on the calibration process without additional observation parameters and is an ideal high-frequency and high-precision calibration algorithm for in-orbit sensors. Therefore, it was applied here to cross calibrate the IRS B08 with the MODIS B31 and B32. The model is expressed by different combination of the TOA temperature of reference channels, as: where T H08 is the TOA temperature of IRS B08 and T M31 , T M32 are the TOA temperatures of MODIS B31, B32, respectively. The coefficients a 1 , a 2 , a 3 can be obtained from the cross coefficients' lookup table with corresponding observed geometric conditions. Figure 5 illustrates the entire process for deducing radiometric coefficients via cross-calibration. The four main steps are described below.

1.
Obtain the DN of IRS B08, the radiance of MODIS B31 and B32 and the view zenith angles of ROI area from near-simultaneous image pairs; 2.
Convert the TOA radiance of MODIS B31, B32 into TOA temperature according to "temperatureradiance" lookup table and obtain cross coefficients, a i , according to "cross coefficient" lookup table; 3.
Calculate the TOA temperature of IRS B08 using Equation (1) and convert to TOA radiance through the inverse of the Planck function or the "temperature-radiance" lookup table.

4.
Perform coefficients regression. The data obtained from the above process (DN and TOA radiance of IRS B08) were regressed using a linear equation as: Remote Sens. 2020, 12, 2362 where DN is the digital number, L is TOA radiance, and b and g are the radiometric calibration coefficients of the expected offset and gain, respectively. The unit of g is W −1 ·m 2 ·sr·µm, and b is unitless.
3. Calculate the TOA temperature of IRS B08 using Equation (1) and convert to TOA radiance through the inverse of the Planck function or the "temperature-radiance" lookup table.
4. Perform coefficients regression. The data obtained from the above process (DN and TOA radiance of IRS B08) were regressed using a linear equation as: where DN is the digital number, L is TOA radiance, and b and g are the radiometric calibration coefficients of the expected offset and gain, respectively. The unit of g is • • • , and b is unitless. Figure 5. Flow chart of cross-calibration by twin-channel difference model.

Field Experiment
The field validation techniques estimated the TOA radiance over a selected test site. Next, the simulated radiance was convolved with the relative channel spectral response function to obtain the predicted radiance at sensor. Finally, the predicted radiance at sensor was compared with the

Field Experiment
The field validation techniques estimated the TOA radiance over a selected test site. Next, the simulated radiance was convolved with the relative channel spectral response function to obtain the predicted radiance at sensor. Finally, the predicted radiance at sensor was compared with the measured radiance. The TOA radiance values were simulated with the help of the radiation propagation code MODTRAN [35] using field measurements as inputs. Measured radiance was tested with a Cimel CE312 radiometer, which was accurately calibrated by a highly accurate field portable blackbody (Mikron M340) before and after the satellite overpass. Water was used as the target. The experiment was conducted concurrent with the satellite transit time.
The temperature-based method for radiation propagation to the sensor can be expressed as: where L TOA is the TOA radiance, τ is the atmospheric transmittance from ground to sensor, ε(λ) is the spectral emissivity of the observed target, B(λ, T) is the Planck function associated with a target at temperature T, L ↓ (λ) is the spectral downwelling radiance, L ↑ (λ) is the channel upwelling radiance, R(λ) is the relative channel spectral response function of the sensor (IRS in this study) and L obs is the channel radiance assuming the sensor was located on the ground. In the radiance-based method, L obs is calculated using Equation (4): where M is the parameter that corresponds to the ground instrument (the Cimel CE312 in this study).
The atmospheric profile data used by MODTRAN to simulate τ, L ↑ and L ↓ were from atmospheric sounding balloons.
The field experiment sites used for this study are listed in Table 2. They are distributed in the southeast coastal area of China, from temperate to subtropical. The tested water temperature ranged from 280 to 295 K.

Annual Calibration Coefficients
The radiometric calibration coefficients of the IRS B08 from 2008 to 2017 obtained in this study are listed in Table 3, while, the official coefficients provided by the CCRSDA from 2008 to 2012 are listed in Table 3. Table 3. Radiometric calibration coefficients of IRS B08 from 2008 to 2017.

Cross-Calibration Coefficients Official Coefficients
Year

Coefficients Validation
The official coefficients were obtained by on-board calibrators, which were used to validate the calibration results obtained in this study. The DNs of water in the images were approximately 400,

On-board Validation
The official coefficients were obtained by on-board calibrators, which were used to validate the calibration results obtained in this study. The DNs of water in the images were approximately 400, with most ranging from 300 to 500. The radiance differences between the results obtained from the official coefficients and this study are shown in Figure 7. It can be seen that all the percent differences of radiance were within ±4.00%. The best fitting result exhibited a 0.16% difference, while the worst fitting result was −3.64%. These results indicate that the cross-calibration coefficients calculated in this study are reliable.

On-board Validation
The official coefficients were obtained by on-board calibrators, which were used to validate the calibration results obtained in this study. The DNs of water in the images were approximately 400, with most ranging from 300 to 500. The radiance differences between the results obtained from the official coefficients and this study are shown in Figure 7. It can be seen that all the percent differences of radiance were within ±4.00%. The best fitting result exhibited a 0.16% difference, while the worst fitting result was −3.64%. These results indicate that the cross-calibration coefficients calculated in this study are reliable.

Vicarious Comparison
Based on the field experiments, the TOA radiance was simulated with the help of MODTRAN. The comparison of the radiance calculated using the cross-calibration coefficients and the field campaign measurements is presented in Figure 8a. They were close to the 1:1 line and the total error was within 0.11 W·m −2 ·sr −1 ·µm −1 , corresponding to approximately 1.0 K of TOA temperature. The percent difference was within 1.63%. The image in 2014 was also calibrated using the calibration coefficients of the past few years. The radiance differences comparison between the cross-calibration and the field campaign is plotted in Figure 8b. These results indicated that not only the accuracy of the calibration coefficients in this study was reliable, but also the performance of the sensor changed with time.

Vicarious Comparison
Based on the field experiments, the TOA radiance was simulated with the help of MODTRAN. The comparison of the radiance calculated using the cross-calibration coefficients and the field campaign measurements is presented in Figure 8a. They were close to the 1:1 line and the total error was within 0.11 W • m • sr • μm , corresponding to approximately 1.0 K of TOA temperature. The percent difference was within 1.63%. The image in 2014 was also calibrated using the calibration coefficients of the past few years. The radiance differences comparison between the cross-calibration and the field campaign is plotted in Figure 8b. These results indicated that not only the accuracy of the calibration coefficients in this study was reliable, but also the performance of the sensor changed with time.

Image Cross-Comparison
The radiance obtained from images of the IRS B08 and the MODIS B31 and B32 over the field experiment sites was compared; the results from these comparisons are plotted in Figure 9. Each data point on the plots corresponds to a sampling points in the field experiments.

Image Cross-Comparison
The radiance obtained from images of the IRS B08 and the MODIS B31 and B32 over the field experiment sites was compared; the results from these comparisons are plotted in Figure 9. Each data point on the plots corresponds to a sampling points in the field experiments. Figure 9a summarizes the radiance comparison results for the IRS B08 vs. MODIS B31. The data points do not lie precisely on the 1:1 line. In fact, all data points lie below the 1:1 line, indicating that IRS B08 was underestimating radiance in comparison to the MODIS B31. In the percent difference plots shown in Figure 9b, IRS B08 radiance estimates are as much as 5% lower. Figure 9c summarizes the radiance comparison results for the IRS B08 vs. MODIS B32. The percent difference plots presented in Figure 9d show that all data points lie above the zero line, indicating that the IRS B08 appears to estimate a 1% to 4% higher radiance than the MODIS B32.

Relatively Low/Normal/High Radiance Variation
According to the applicable conditions of the twin-channel difference model, coefficients a could be applied to all images in which the TOA temperature ranged from 240 to 310 K. When the temperature was 240 K, the corresponding radiance was approximately 3.0 W • m • sr • μm . When the temperature was 310 K, the corresponding radiance was approximately 10.0 W • m • sr • μm . The water DNs in the selected images were within 500, corresponding to an at-aperture radiation of approximately 8.5 W • m • sr • μm . The corresponding DN values from low radiance (3.0 W • m • sr • μm ) to high radiance (10.0 W • m • sr • μm ) were computed using the crosscalibration coefficients obtained in this study, and the results are presented in Figure 10. These results indicate that the DN at low radiance increased slightly from 2008 to 2013 and fluctuated after that time. It was more stable at normal radiance, however. When the radiance was 8.0 W • m • sr • Figure 9. Comparison of radiance obtained from images of both IRS and MODIS sensors. IRS radiances were calculated using the calibration coefficients of this study. The MODIS radiances were obtained from MODIS L1B products. (a) IRS B08 vs. MODIS B31 radiance; (b) percent difference relative to MODIS B31; (c) IRS B08 vs. MODIS B32 radiance; (d) percent difference relative to MODIS B32.

Relatively Low/Normal/High Radiance Variation
According to the applicable conditions of the twin-channel difference model, coefficients a i could be applied to all images in which the TOA temperature ranged from 240 to 310 K. When the temperature was 240 K, the corresponding radiance was approximately 3.0 W·m −2 ·sr −1 ·µm −1 . When the temperature was 310 K, the corresponding radiance was approximately 10.0 W·m −2 ·sr −1 ·µm −1 . The water DNs in the selected images were within 500, corresponding to an at-aperture radiation of approximately 8.5 W·m −2 ·sr −1 ·µm −1 . The corresponding DN values from low radiance (3.0 W·m −2 ·sr −1 ·µm −1 ) to high radiance (10.0 W·m −2 ·sr −1 ·µm −1 ) were computed using the cross-calibration coefficients obtained in this study, and the results are presented in Figure 10. These results indicate that the DN at low radiance increased slightly from 2008 to 2013 and fluctuated after that time. It was more stable at normal radiance, however. When the radiance was 8.

Channel Response Variation
As shown in Figure 11a, the TOA radiance range in the selected grids was 5.0-8.0 W • m • sr • μm . Therefore, we chose 5.0-8.0 W • m • sr • μm as the input radiance interval to analyze the response of the sensor. The corresponding DNs were calculated with Equation (2) using the coefficients obtained in this study. The change in the DN range is presented in Figure 11b and the attenuation percentages are shown in Figure 11c. The larger the DN range, the higher the sensitivity of the sensor. It can clearly be seen from Figure 11

Channel Response Variation
As shown in Figure 11a, the TOA radiance range in the selected grids was 5.0-8.0 W·m −2 ·sr −1 ·µm −1 . Therefore, we chose 5.0-8.0 W·m −2 ·sr −1 ·µm −1 as the input radiance interval to analyze the response of the sensor. The corresponding DNs were calculated with Equation (2) using the coefficients obtained in this study. The change in the DN range is presented in Figure 11b and the attenuation percentages are shown in Figure 11c. The larger the DN range, the higher the sensitivity of the sensor. It can clearly be seen from Figure 11

Channel Response Variation
As shown in Figure 11a, the TOA radiance range in the selected grids was 5.0-8.0 W • m • sr • μm . Therefore, we chose 5.0-8.0 W • m • sr • μm as the input radiance interval to analyze the response of the sensor. The corresponding DNs were calculated with Equation (2) using the coefficients obtained in this study. The change in the DN range is presented in Figure 11b and the attenuation percentages are shown in Figure 11c. The larger the DN range, the higher the sensitivity of the sensor. It can clearly be seen from Figure 11  In summary, the radiometric response of the IRS B08 changed significantly during its orbital lifetime. A negative number indicates that the radiometric response was lower than the previous year. In other words, under the same radiance input interval, the DN range became smaller. The radiometric response exhibited a downward trend year-by-year from 2009 to 2013, with the largest decrease occurring in 2010. The radiometric response increased by +5.23% in 2014, however and fluctuated after that time.

Discussion
Hundreds of near-simultaneous images of Lake Qinghai were collected from 2008 to 2017. This provided a long-term sequence of images for the thermal infrared cross-calibration of the instrument. Temperature linearization between the IRS B08 and MODIS B31 was performed before crosscalibration, as shown in Figure 12a. They manifested a very close linear correlation. All data points lie very near the 1:1 line, indicating excellent temperature agreement between the two sensors. The percent difference plots (Figure 12b) suggest a difference of approximately 1% relative to IRS. Generally, in Figure 12c, the MODIS B31 displays a higher radiance than the IRS B08, while the MODIS B32 has a lower radiance than the IRS B08. Radiance and temperature distribution as a function of the MODIS and IRS lifetimes are shown in Figure 12d. In different years, each channel exhibited sinusoidal trends. Therefore, they were effective in further analysis. In summary, the radiometric response of the IRS B08 changed significantly during its orbital lifetime. A negative number indicates that the radiometric response was lower than the previous year. In other words, under the same radiance input interval, the DN range became smaller. The radiometric response exhibited a downward trend year-by-year from 2009 to 2013, with the largest decrease occurring in 2010. The radiometric response increased by +5.23% in 2014, however and fluctuated after that time.

Discussion
Hundreds of near-simultaneous images of Lake Qinghai were collected from 2008 to 2017. This provided a long-term sequence of images for the thermal infrared cross-calibration of the instrument. Temperature linearization between the IRS B08 and MODIS B31 was performed before cross-calibration, as shown in Figure 12a. They manifested a very close linear correlation. All data points lie very near the 1:1 line, indicating excellent temperature agreement between the two sensors. The percent difference plots (Figure 12b) suggest a difference of approximately 1% relative to IRS. Generally, in Figure 12c, the MODIS B31 displays a higher radiance than the IRS B08, while the MODIS B32 has a lower radiance than the IRS B08. Radiance and temperature distribution as a function of the MODIS and IRS lifetimes are shown in Figure 12d. In different years, each channel exhibited sinusoidal trends. Therefore, they were effective in further analysis. In summary, the radiometric response of the IRS B08 changed significantly during its orbital lifetime. A negative number indicates that the radiometric response was lower than the previous year. In other words, under the same radiance input interval, the DN range became smaller. The radiometric response exhibited a downward trend year-by-year from 2009 to 2013, with the largest decrease occurring in 2010. The radiometric response increased by +5.23% in 2014, however and fluctuated after that time.

Discussion
Hundreds of near-simultaneous images of Lake Qinghai were collected from 2008 to 2017. This provided a long-term sequence of images for the thermal infrared cross-calibration of the instrument. Temperature linearization between the IRS B08 and MODIS B31 was performed before crosscalibration, as shown in Figure 12a. They manifested a very close linear correlation. All data points lie very near the 1:1 line, indicating excellent temperature agreement between the two sensors. The percent difference plots (Figure 12b) suggest a difference of approximately 1% relative to IRS. Generally, in Figure 12c, the MODIS B31 displays a higher radiance than the IRS B08, while the MODIS B32 has a lower radiance than the IRS B08. Radiance and temperature distribution as a function of the MODIS and IRS lifetimes are shown in Figure 12d. In different years, each channel exhibited sinusoidal trends. Therefore, they were effective in further analysis. In Section 2.1.3, we discarded image pairs with zenith angles >50° in the MODIS view. We selected images from 2011 in order to verify the validity of this choice, since 2011 had the largest number of effective image pairs in this study. As shown in Table 4, calibration accuracy is related to zenith angle. When images with a zenith angle > 50° were included, the worst fitting result occurred. Therefore, it was effective to filter the images based on the zenith angle of the MODIS view. As described in Section 3.3, the radiometric response of the IRS B08 at normal at-aperture radiance was relatively stable. Therefore, we regarded it as a fixed regression point (L = 8.0 W • m • sr • μm , DN = 460). Combined with another point obtained by cross-calibration from each image pair, the linear regression coefficients could be calculated. The gains and offsets are shown in Figure  13a, b. Both gains and offsets were more dispersed in summer than in other seasons. The average surface temperature of Lake Qinghai in summer is approximately 289.15 K, and the highest temperature is 295.45 K, corresponding to the normal radiance range discussed in Section 3.3. The radiation response in this range was stable, so the gap between the DN and TOA radiation derived from MODIS in summer and the fixed point was relatively small. Therefore, the slope of the fitted line was small. A slight change in slope will cause a large change in the reciprocal, which is the gain. As a result, gains in summer exhibited more volatility. Correspondingly, the offsets also fluctuated more in summer than in other seasons. However, gains and offsets in winter were more stable, due to the fact that DN changed more when the temperature (radiance) was low. In Section 2.1.3, we discarded image pairs with zenith angles >50 • in the MODIS view. We selected images from 2011 in order to verify the validity of this choice, since 2011 had the largest number of effective image pairs in this study. As shown in Table 4, calibration accuracy is related to zenith angle. When images with a zenith angle > 50 • were included, the worst fitting result occurred. Therefore, it was effective to filter the images based on the zenith angle of the MODIS view. As described in Section 3.3, the radiometric response of the IRS B08 at normal at-aperture radiance was relatively stable. Therefore, we regarded it as a fixed regression point (L = 8.0 W·m −2 ·sr −1 ·µm −1 , DN = 460). Combined with another point obtained by cross-calibration from each image pair, the linear regression coefficients could be calculated. The gains and offsets are shown in Figure 13a,b. Both gains and offsets were more dispersed in summer than in other seasons. The average surface temperature of Lake Qinghai in summer is approximately 289.15 K, and the highest temperature is 295.45 K, corresponding to the normal radiance range discussed in Section 3.3. The radiation response in this range was stable, so the gap between the DN and TOA radiation derived from MODIS in summer and the fixed point was relatively small. Therefore, the slope of the fitted line was small. A slight change in slope will cause a large change in the reciprocal, which is the gain. As a result, gains in summer exhibited more volatility. Correspondingly, the offsets also fluctuated more in summer than in other seasons. However, gains and offsets in winter were more stable, due to the fact that DN changed more when the temperature (radiance) was low.
Overall, the cross-calibration method is a cost-effective technique that is particularly tailored for long-term historical data. The attenuation characteristics of the sensor during its lifetime can be further analyzed based on the calibration coefficients. The coefficients of the IRS B08, especially from 2013 to 2017, provides an accurate guarantee for the quantitative application of IRS data. These data can be further applied to waterline extraction, temperature retrieval and other areas. The research results have fully enhanced the value of the HJ-1B satellite data. Overall, the cross-calibration method is a cost-effective technique that is particularly tailored for long-term historical data. The attenuation characteristics of the sensor during its lifetime can be further analyzed based on the calibration coefficients. The coefficients of the IRS B08, especially from 2013 to 2017, provides an accurate guarantee for the quantitative application of IRS data. These data can be further applied to waterline extraction, temperature retrieval and other areas. The research results have fully enhanced the value of the HJ-1B satellite data.

Conclusions
Cross-calibration between the IRS B08 and MODIS B31, B32 was performed using nearsimultaneous image pairs of Lake Qinghai over a 10-year lifetime (2008-2017). All images impacted by cloud or ice were discarded. In addition, since the calibration accuracy decreased significantly when images corresponding to MODIS zenith angles >50° were included, those images were excluded as well. The remaining 152 pairs were used to perform radiometric calibration. After various corrections, a very close linear correlation was found between the TOA radiances of IRS B08 and the DNs extracted from HJ-1B IRS images. The cross-calibration coefficients, gains and offsets, were obtained from 2008 to 2017, which compensated for the missing coefficients from 2013 to 2017. The percent difference was within ±4.00% when compared to the on-board calibrator results. Compared with the field campaign results, the percent difference was within 1.63%. The IRS B08 radiance is lower than that of MODIS B31 and higher than MODIS B32 radiance. However, the TOA temperature between the two sensors are highly consistent with approximately 1% difference. The calibration coefficients obtained in this study are highly accurate and can thus be used in the radiometric calibration of HJ-1B IRS images. The percent changes of gains were within ±6.5%, although the percent changes of b were more scattered, indicating that the fluctuations of offsets were greater than those of gains. For the same radiance, the corresponding DN varied in different years.
The DN values changed more in the low radiance range but were relatively stable in the normal radiance range. Therefore, the radiance in the normal range of L = 8.0 W • m • sr • μm was regarded as a fixed point corresponding to DN = 460. Combined with another point obtained based on each image, cross-coefficients were calculated. The coefficients manifested different change characteristics in different seasons. The gain and offset values in summer were more volatile than in other seasons, while the winter values were the most stable. As a result, the DN range of the sensor differed for the same radiance input interval, meaning that the radiometric response of the IRS changed during its orbital lifetime. The response of the IRS decreased from 2008 to 2013 but fluctuated from 2014 to 2017.
With the calibration coefficients provided here, the IRS thermal archive calibrated the radiance to acceptable levels. SST and LST maps can be operationally produced from the IRS radiance data,

Conclusions
Cross-calibration between the IRS B08 and MODIS B31, B32 was performed using nearsimultaneous image pairs of Lake Qinghai over a 10-year lifetime (2008-2017). All images impacted by cloud or ice were discarded. In addition, since the calibration accuracy decreased significantly when images corresponding to MODIS zenith angles >50 • were included, those images were excluded as well. The remaining 152 pairs were used to perform radiometric calibration. After various corrections, a very close linear correlation was found between the TOA radiances of IRS B08 and the DNs extracted from HJ-1B IRS images. The cross-calibration coefficients, gains and offsets, were obtained from 2008 to 2017, which compensated for the missing coefficients from 2013 to 2017. The percent difference was within ±4.00% when compared to the on-board calibrator results. Compared with the field campaign results, the percent difference was within 1.63%. The IRS B08 radiance is lower than that of MODIS B31 and higher than MODIS B32 radiance. However, the TOA temperature between the two sensors are highly consistent with approximately 1% difference. The calibration coefficients obtained in this study are highly accurate and can thus be used in the radiometric calibration of HJ-1B IRS images. The percent changes of gains were within ±6.5%, although the percent changes of b were more scattered, indicating that the fluctuations of offsets were greater than those of gains. For the same radiance, the corresponding DN varied in different years. The DN values changed more in the low radiance range but were relatively stable in the normal radiance range. Therefore, the radiance in the normal range of L = 8.0 W·m −2 ·sr −1 ·µm −1 was regarded as a fixed point corresponding to DN = 460. Combined with another point obtained based on each image, cross-coefficients were calculated. The coefficients manifested different change characteristics in different seasons. The gain and offset values in summer were more volatile than in other seasons, while the winter values were the most stable. As a result, the DN range of the sensor differed for the same radiance input interval, meaning that the radiometric response of the IRS changed during its orbital lifetime. The response of the IRS decreased from 2008 to 2013 but fluctuated from 2014 to 2017.
With the calibration coefficients provided here, the IRS thermal archive calibrated the radiance to acceptable levels. SST and LST maps can be operationally produced from the IRS radiance data, and these data can also be used in various applications, such as waterline extraction, temperature retrieval and other research areas.