Radiometric Variations of On-Orbit FORMOSAT-5 RSI from Vicarious and Cross-Calibration Measurements

: A new Taiwanese satellite, FORMOSAT-5 (FS-5), with a payload remote sensing instrument (RSI) was launched in August 2017 to continue the mission of its predecessor FORMOSAT-2 (FS-2). Similar to FS-2, the RSI provides 2-m resolution panchromatic and 4-m resolution multi-spectral images as the primary payload on FS-5. However, the radiometric properties of the optical sensor may vary, based on the environment and time after being launched into the space. Thus, maintaining the radiometric quality of FS-5 RSI imagery is essential and signiﬁcant to scientiﬁc research and further applications. Therefore, the objective of this study aimed at the on-orbit absolute radiometric assessment and calibration of on-orbit FS-5 RSI observations. Two renowned approaches, vicarious calibrations and cross-calibrations, were conducted at two calibration sites that employ a stable atmosphere and high surface reﬂectance, namely, Alkali Lake and Railroad Valley Playa in North America. For cross-calibrations, the Landsat-8 Operational Land Imager (LS-8 OLI) was selected as the reference. A second simulation of the satellite signal in the solar spectrum (6S) radiative transfer model was performed to compute the surface reﬂectance, atmospheric e ﬀ ects, and path radiance for the radiometric intensity at the top of the atmosphere. Results of vicarious calibrations from 11 ﬁeld experiments demonstrated high consistency with those of seven case examinations of cross-calibration in terms of physical gain in spectra, implying that the practicality of the proposed approaches is high. Moreover, the multi-temporal results illustrated that RSI decay in optical sensitivity was evident after launch. The variation in the calibration coe ﬃ cient of each band showed no obvious consistency (6%–24%) in 2017, but it tended to be stable at the order of 3%–5% of variation in most spectral bands during 2018. The results strongly suggest that periodical calibration is required and essential for further scientiﬁc applications.


Introduction
FORMOSAT-5 (FS-5), launched in August 2017, is the first indigenously developed Earth observation satellite by the National Space Organization of Taiwan. FS-5 is the follow up satellite of FS-2 with a 2-day revisiting orbit. FS-5 carries a Taiwanese, domestic-made remote sensor, remote sensing instrument (FS-5 RSI). FS-5 RSI has a similar spectral bandwidth and radiometric properties to FS-2 RSI, but it has a higher radiometric resolution (8-bit to 12-bit) and spatial resolution in multi-spectral To completely evaluate and calibrate the variations in radiometric parameters of FS-5 RSI postlaunch, this study intends to examine the performance of vicarious [8,12,[19][20][21][22][23][24] and cross-calibrations [13,[24][25][26][27][28][29][30]. Vicarious calibrations will be performed using FS-5 RSI imagery based on simultaneous measurements from ground-based instruments including the effects of the spectral band difference effect (SBDE), viewing angles, and atmosphere parameters. On the other hand, the Landsat-8 calibration and validation team (CVT) has continuously monitored the stability of the instrument with onboard vicarious calibration and pseudo invariant calibration site (PICS) based calibration methods [10,11]. An extensive study of the LS-8 OLI instrument and its preflight calibration has been described in previous literature [17,31,32]. Therefore, the well-calibrated OLI sensor was selected to serve as the cross-calibration reference in this study. The original designs of FS-2 and FS-5 did not include on-board calibration systems such as an internal calibrator (IC), full aperture solar calibrator (FASC), partial aperture solar calibrator (PASC), and solar diffuser [21,33,34]. Consequently, the objective of this study aimed to construct a viable scheme of operational radiometric calibration onorbit for FS-5 RSI by integrating the advantages of vicarious and cross approaches. The results of radiometric calibration could assist the National Space Organization in maintaining the radiometric quality of FS-5 RSI imagery.

Calibration Sites
Proper calibration sites are required for on-orbit absolute radiometric calibration. To minimize the environmental and BRDF effects, desert areas are usually preferable due to the clear weather with low aerosol loading and the near-Lambertian surfaces [18,26,35]. In addition, the high surface reflectivity of these areas enhances the signal-to-noise ratio (SNR) more effectively. Thus, the Second, vicarious calibration uses natural reference targets for instrument calibration coincident with a satellite overpass. This method has been applied to estimate absolute radiance calibrations of FS-2 RSI using an airport runway as the calibration site on Dongsha Island [5]. With the input of in situ atmospheric parameters and surface reflectance, the TOA (top of atmosphere) radiance of each spectral band is computed by the LOWTRAN-7 radiation transfer model. Compared to the coincident digital number recorded by the FS-2 RSI sensor, results indicated there were obvious variations in physical gains (approximately 20%) in the period after launch in September 2004. In addition, Liu et al. [15] conducted radiometric calibrations of the FORMOSAT-2 RSI based on ground data collected by the advanced spaceborne thermal emission and reflection radiometer (ASTER) calibration team at Ivanpah Playa and Railroad Valley in the United States from 2006 to 2008 [6]. The vicarious calibration results showed that the degradations of sensitivity were 8% to 12% (visible), 9.6% (PAN), and 2.2% (near-infrared), respectively. However, several prerequisites needed to be considered such as homogeneous and bright surfaces, atmospheric stability, and cloud cover effects of the calibration site. Therefore, vicarious calibration is rather complex, time consuming, and laborious, although it is practical in that it gives precise results. Third, the cross-calibration approach compares the radiation delivered by sensors on different satellites, which is generally applied for monitoring and identifying the radiometric performance of satellite optical sensors. Cross-calibration is typically done through comparisons of band-to-band data, which have overlapping spectra and similar spectral response functions of different sensors. For instance, the spectral band difference effect (SBDE) of Landsat-8 OLI (Operational Land Imager) and Landsat-7 ETM+ (Enhanced Thematic Mapper Plus) sensors was investigated by Mishra et al. [16], and the results showed that Earth Observing 1 (EO-1) Hyperion images could be used to compensate the relative spectral response differences. In addition, by using the advanced spaceborne thermal emission and reflection radiometer (ASTER) onboard Terra as a reference, Kamei et al. [17] presented results on the cross-calibration of FS-2 RSI based on the sites at La Crau, Taklimakan Desert, Ivanpah Playa, Alkali Lake, and Railroad Valley Playa. In their study, TOA radiance was derived from the second simulation of a satellite signal in the solar spectrum (6S) model including the effects of water vapor (H 2 O), ozone (O 3 ), and aerosols. Non-uniform reflection characteristics of the surface were also taken into account using the MODIS bidirectional reflectance distribution function (BRDF). The results of the BRDF corrections at Ivanpah Playa and Railroad Valley Playa sites were within 5% when compared to the preflight values from 2012. However, major error sources from this approach include observation time dissimilarity, spectral response function differences in conjunction with spectral surface reflectance and atmospheric optical depth, and observation area difference [18].
To completely evaluate and calibrate the variations in radiometric parameters of FS-5 RSI post-launch, this study intends to examine the performance of vicarious [8,12,[19][20][21][22][23][24] and cross-calibrations [13,[24][25][26][27][28][29][30]. Vicarious calibrations will be performed using FS-5 RSI imagery based on simultaneous measurements from ground-based instruments including the effects of the spectral band difference effect (SBDE), viewing angles, and atmosphere parameters. On the other hand, the Landsat-8 calibration and validation team (CVT) has continuously monitored the stability of the instrument with onboard vicarious calibration and pseudo invariant calibration site (PICS) based calibration methods [10,11]. An extensive study of the LS-8 OLI instrument and its preflight calibration has been described in previous literature [17,31,32]. Therefore, the well-calibrated OLI sensor was selected to serve as the cross-calibration reference in this study. The original designs of FS-2 and FS-5 did not include on-board calibration systems such as an internal calibrator (IC), full aperture solar calibrator (FASC), partial aperture solar calibrator (PASC), and solar diffuser [21,33,34]. Consequently, the objective of this study aimed to construct a viable scheme of operational radiometric calibration on-orbit for FS-5 RSI by integrating the advantages of vicarious and cross approaches. The results of radiometric calibration could assist the National Space Organization in maintaining the radiometric quality of FS-5 RSI imagery.

Calibration Sites
Proper calibration sites are required for on-orbit absolute radiometric calibration. To minimize the environmental and BRDF effects, desert areas are usually preferable due to the clear weather with low aerosol loading and the near-Lambertian surfaces [18,26,35]. In addition, the high surface reflectivity of these areas enhances the signal-to-noise ratio (SNR) more effectively. Thus, the calibration sites at Alkali Lake and Railroad Valley Playa, located in the Great Basin Desert of the USA, which are regularly used for radiometric calibrations of various satellite sensors due to their higher spatial uniformity and temporal stability [9][10][11]18,[23][24][25][26]30,[36][37][38], were selected in this study. Figures 2 and 3 illustrate subsets of the FS-5 RSI Level 1A images covering the corresponding sites.  The surface reflectance of the calibration sites for vicarious calibrations measured not only the reflectance at the sensor viewing angle, but also the reflectance at different zenith and azimuth angles to obtain the complete reflectivity information of the calibration sites, as shown in Figures 3 and 4, where most surface reflectivity of the calibration sites were near-Lambertian at the same zenith angle along with azimuth angles. In this study, this measurement provided BRDF and SBDE corrections for cross-calibration in the difference of observing geometry and spectral bandwidth between FS-5 and L-8.     The surface reflectance of the calibration sites for vicarious calibrations measured not only the reflectance at the sensor viewing angle, but also the reflectance at different zenith and azimuth angles to obtain the complete reflectivity information of the calibration sites, as shown in Figures 3 and 4, where most surface reflectivity of the calibration sites were near-Lambertian at the same zenith angle along with azimuth angles. In this study, this measurement provided BRDF and SBDE corrections for cross-calibration in the difference of observing geometry and spectral bandwidth between FS-5 and L-8.   The surface reflectance of the calibration sites for vicarious calibrations measured not only the reflectance at the sensor viewing angle, but also the reflectance at different zenith and azimuth angles to obtain the complete reflectivity information of the calibration sites, as shown in Figures 3 and 4, where most surface reflectivity of the calibration sites were near-Lambertian at the same zenith angle along with azimuth angles. In this study, this measurement provided BRDF and SBDE corrections for cross-calibration in the difference of observing geometry and spectral bandwidth between FS-5 and L-8.

Satellite Datasets
Eighteen cloudless scenes acquired by FS-5 RSI during the period from November 9, 2017 to September 25, 2018 were collected. Eleven of them were examined with vicarious calibration, as shown in Table 1, and seven sets of FS-5 RSI images associated with concurrent LS-8 OLI images were used for cross-calibrations, as presented in Table 2. Level 1A products of FS-5 RSI were prepared for both radiometric calibrations, and LS-8 OLI level 1T products were taken as the reference imagery for cross-calibration. The Level 1T imagery of OLI has been radiometrically calibrated, geometrically corrected, and coordinately projected [18]. The solar zenith angle (SZA), solar azimuth angle (SAA), view zenith angle (VZA), and view azimuth angle (VAA) for both sensors are also included in Tables  1 and 2.

Satellite Datasets
Eighteen cloudless scenes acquired by FS-5 RSI during the period from November 9, 2017 to September 25, 2018 were collected. Eleven of them were examined with vicarious calibration, as shown in Table 1, and seven sets of FS-5 RSI images associated with concurrent LS-8 OLI images were used for cross-calibrations, as presented in Table 2. Level 1A products of FS-5 RSI were prepared for both radiometric calibrations, and LS-8 OLI level 1T products were taken as the reference imagery for cross-calibration. The Level 1T imagery of OLI has been radiometrically calibrated, geometrically corrected, and coordinately projected [18]. The solar zenith angle (SZA), solar azimuth angle (SAA), view zenith angle (VZA), and view azimuth angle (VAA) for both sensors are also included in Tables 1  and 2.

Atmospheric Parameters
For visible and NIR regions of the EM wave spectrum, the effects of atmospheric scattering and absorption on radiation transfer were significant, particularly the aerosol optical properties. Most atmosphere parameters are built-up in the radiative transfer code (RTC), except for aerosol optical properties. The aerosol profile and aerosol optical depth at 550 nm (AOD550) are required for TOA radiance calculations in the RTC. In this study, the aerosol profile was set in desert areas at the CNES (National Centre for Space Studies) as a reference list of calibration sites [36]. The AOD550 was retrieved from the AERONET (AErosol RObotic NETwork) station at the Alkali Lake and Railroad Valley Playa sites. For calibration sites without the AERONET station, AOD550 could be derived from the satellite products such as Terra MODIS (Moderate Resolution Imaging Spectroradiometer) (MOD04_L2) [39] and CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) (CAL_LID_L2_05kmAPro) products [40].
In addition to aerosols, other atmospheric materials were also included in radiative transfer such as the absorption of ozone and water vapor. In the study, the data of ozone and water vapor were collected from on-site measurements of a sun photometer and are listed in Table 3.

Radiation Transfer of Satellite Observations
Under the assumption of the Lambertian ground surface in the 6S model simulation, TOA reflectance is determined by the TOA radiance and solar constant, as shown in Equation (1) [15,20]: where ρ * is the TOA reflectance observed by the satellite's sensor; µ s is the cosine of the solar zenith angle; µ v is the cosine of the sensor viewing zenith angle; ρ a is the path radiance of atmosphere; ϕ is the relative azimuth angle between the sun and the sensor; τ is the atmosphere optical depth; ρ is the ground surface reflectance; T is the total transmittance from the sun to the ground surface; s is the atmosphere albedo; ρ is the mean surface reflectance; and t d is the diffuse transmittance on the surface-to-sensor path. In this study, the input parameters for the RTC simulation were surface reflectance, atmospheric properties (water vapor, ozone, aerosol vertical profile, and optical depth and type), sensor geometry (SZA, SAA, VZA, and VAA), and spectral response function [41][42][43].

Atmospheric Effects
Desert regions are usually selected as calibration sites because of their atmospheric stability and clear sky conditions, which can reduce the uncertainties and complexity in calculating the radiation transfer. However, there are still some atmospheric parameters that affect radiative transfer in visible and NIR regions including aerosol, ozone, and water vapor content. Thus, radiometric calibrations are usually implemented with measurements of atmospheric components. When there is a lack of in situ measurements of atmospheric parameters, the atmospheric model (US Standard 1962) built-in 6S code and satellite aerosol products are employed to provide synchronous observed atmospheric data. The US Standard 1962 atmospheric profile is a model that defines values and profiles of atmospheric temperature, pressure, and content (Table 4). However, the atmospheric model may significantly affect the simulation output of transported radiance ( Figure 5) when the actual atmospheric state differs greatly from the empirical atmospheric model, which results in errors for subsequent radiometric calibration procedures. In order to assess the uncertainty of simulated at-sensor radiance in vicarious calibration (reflectance based on FS-2 RSI bandwidth) and cross-calibration (reflectance based on LS-8 OLI bandwidth), which are affected by atmospheric components, the coefficient of variation (CV) is quantified as where σ(R i,a ) and R i,a are the standard deviation and average of at-sensor radiance simulated under a atmospheric conditions in the i band, respectively.  SBDE is a prerequisite for cross-calibration. Previous studies have shown that for the same target area, the difference of observed TOA radiance can reach as high as 25% by different sensors with similar spectral bands. By computing the relative difference of at-sensor radiance of vicarious

Spectral Band Difference Effect (SBDE)
SBDE is a prerequisite for cross-calibration. Previous studies have shown that for the same target area, the difference of observed TOA radiance can reach as high as 25% by different sensors with similar spectral bands. By computing the relative difference of at-sensor radiance of vicarious calibration (reflectance based on FS-5 RSI bandwidth) and cross-calibration (reflectance based on LS-8 OLI bandwidth), the SBDE for the same atmosphere condition is defined as follows: where R cross,i,a is the at-sensor radiance for band i under atmosphere condition a based on the cross calibration, and R vica.,i,a is the at-sensor radiance for band i under atmosphere condition a based on the vicarious calibration.

Environmental Effect
Experiments were also set up to determine the environmental effect in AERONET of visible band AOD (aerosol optical depth) products at 550 nm in various target sizes centered at calibration sites. The reflectance of the calibration site (ROC) and the reflectance of the environment (ROE) are fixed and different values. We altered the radius centered at the calibration site under the same atmosphere parameters and geometry conditions.

Relationship between Radiation and Digital Number of FS-5 RSI
Regarding the digital number (DN) output of the FS-5 RSI, conversion between the detected signals to the spectral at-sensor radiance in each band can be described from the following equation [15], where C(R) is the output signal (DN); K 0 is the mean calibration coefficient, depending on the band; G sel is the dynamic range gain equaling to 2 j−1 2 where j is a gain number (for FS-5 RSI, j = 1, 2, 4); η is the relative response of each detector in the CMOS array; R is the mean spectral radiance (W/m 2 /sr/µm) observed by FS-5 RSI; N is the noise; and C(0) is the offset. Considering the high SNR requirement (92 in the PAN band and over 140 in the MS bands) [6], the noise was ignorable. For the FS-5 product level 1A or above, C(0) was equal to 0, and η was equal to 1.

Vicarious Calibration
In this study, vicarious calibration is a reflectance-based method, which relies on ground surface reflectance measurements. The in situ ground surface reflectance at Alkali Lake and Railroad Valley were measured with a FieldSpec 3 portable spectroradiometer (made by Analytical Spectral Devices, Inc.) in the spectral range from 350 nm to 2500 nm. A SPECTRALON (made by Spectra Vista Corporation) standard reflection panel was used as the estimated ideal diffuse standard surface. Then, the in situ ground surface reflectance was measured simultaneously with the overpass of FS-5 RSI. Based on the ground surface reflectance measurements and atmospheric parameters, at-sensor radiance of FS-5 RSI was simulated by the RTC model (6S). Then, the mean calibration coefficient K 0 was calculated from the grey level of the pixel, at-sensor radiance, and corresponding gain number based on Equation (4). The workflow of this procedure is shown in Figure 6. In order to facilitate the comparison Remote Sens. 2019, 11, 2634 9 of 18 of calibration results of this study with the latest instrument calibration values (i.e., dim files provided by NSPO), the relative difference between calibration values is represented as the following: where K 0.i,vica is the radiometric conversion coefficient computed from vicarious calibration in band i, and K 0.i,dim is the radiometric conversion coefficient provided under dim files in band i.
Corporation) standard reflection panel was used as the estimated ideal diffuse standard surface. Then, the in situ ground surface reflectance was measured simultaneously with the overpass of FS-5 RSI. Based on the ground surface reflectance measurements and atmospheric parameters, at-sensor radiance of FS-5 RSI was simulated by the RTC model (6S). Then, the mean calibration coefficient K0 was calculated from the grey level of the pixel, at-sensor radiance, and corresponding gain number based on Equation (4). The workflow of this procedure is shown in Figure 6. In order to facilitate the comparison of calibration results of this study with the latest instrument calibration values (i.e., dim files provided by NSPO), the relative difference between calibration values is represented as the following: . , * 100% where . , is the radiometric conversion coefficient computed from vicarious calibration in band i, and . , is the radiometric conversion coefficient provided under dim files in band i.

Cross-Calibration
Cross-calibration is generally employed for post-launch satellite sensors [44]. Nevertheless, coincidences are rare when comparing the two satellite sensors with different cycles, orbits, swaths, spatial, and local overpass times, and so on [1,45]. Thus, certain corrections have to be exerted on these aspects. First, the selected pairs of images for FS-5 RSI and LS-8 OLI were geometrically corrected to match with each other. This step is to maintain the consistency across all images and to avoid any discrepancies during the subsequent process. Second, FS-5 RSI images were cropped to the desired size within the areas of LS-8 OLI images. Then, the images were resampled to 25 by 25 grids

Cross-Calibration
Cross-calibration is generally employed for post-launch satellite sensors [44]. Nevertheless, coincidences are rare when comparing the two satellite sensors with different cycles, orbits, swaths, spatial, and local overpass times, and so on [1,45]. Thus, certain corrections have to be exerted on these aspects. First, the selected pairs of images for FS-5 RSI and LS-8 OLI were geometrically corrected to match with each other. This step is to maintain the consistency across all images and to avoid any discrepancies during the subsequent process. Second, FS-5 RSI images were cropped to the desired size within the areas of LS-8 OLI images. Then, the images were resampled to 25 by 25 grids [16]. Finally, the brightness grids were selected for cross-calibration. The selected grids of LS-8 OLI images (at-sensor radiance) were used as input in RTC (6S model) with the corresponding geometric conditions, atmosphere parameters, and spectral response function (SRF) to remove the atmosphere effect and obtain the spectral surface reflectance of LS-8 OLI. The corresponding spectral surface reflectance (Re f l i ) of each sensor can be described as the following equation: where λ Max (i) is the maximum wavelength of the Full Width at Half Maximum (FWHM) of band i; λ min (i) is the minimum wavelength of the FWHM of band i; and Refl(λ) and SRF(λ) are the surface reflectance and spectral response at wavelength λ, respectively.
To properly correct the viewing geometry on cross calibration, the difference in the spectral bandwidth and BRDF effect were also considered. Despite the fact that the calibration site surface was near-Lambertian, the BRDF effect remained significant. The ground based in situ measurements of surface reflectance were used to quantify the influence for the correction of these two effects. The difference (δ band ) of spectral radiance between FS-5 RSI and LS-8 OLI in the same observing geometry can be estimated by Equation (7) below, where Refl OLI and Refl RSI are the mean surface reflectance within the spectral bandwidth of LS-8 OLI and FS-5 RSI, respectively. The BRDF effect on surface reflectance between nadir and off-nadir (≥ 50 0 ), δ viewing , is estimated by the following equation, Re f l nadir * 100 (8) where Refl off−nadir and Refl nadir are the mean surface reflectance in off-nadir and in nadir, respectively.
The coefficient (C) for the BRDF and spectral corrections is then derived from the combination of Equation (7) and Equation (8) and rearranged as, where Refl RSI, off−nadir and Refl OLI, nadir are the mean surface reflectance of FS-2 RSI and LS-8 OLI, respectively. The flow chart of the cross-calibration is shown in Figure 7. The angular effect of surface reflectance was corrected with the BRDF according to Equation (9). Consequently, based on geometric conditions, the SRF of FS-5 RSI, and the same atmospheric parameters at the LS-8 OLI and FS-5 RSI scene acquisition times, the corresponding at-sensor radiance of FS-5 RSI was simulated using the RTC model.

Sensitivity Study: At-Sensor Radiance/Reflectance Simulation
The FS-5 RSI at-sensor radiance for all bands was calculated by 6S simulation. Errors in atmospheric parameters (water vapor, ozone content, and AOD) may create potential uncertainties in RTC simulation. Therefore, the sensitivity of simulated FS-5 RSI at-sensor radiance was tested by changing the content of the control variables (ozone, water vapor, and AOD). The test results ( Figure  8) showed that the maximum differences of CV values in at-sensor radiance calculations were

Sensitivity Study: At-Sensor Radiance/Reflectance Simulation
The FS-5 RSI at-sensor radiance for all bands was calculated by 6S simulation. Errors in atmospheric parameters (water vapor, ozone content, and AOD) may create potential uncertainties in RTC simulation. Therefore, the sensitivity of simulated FS-5 RSI at-sensor radiance was tested by changing the content of the control variables (ozone, water vapor, and AOD). The test results (Figure 8) showed that the maximum differences of CV values in at-sensor radiance calculations were approximately 4%, 2%, and 4% for the different control variables, separately. Most CV values had similar tendencies in Alkali Lake and Railroad Valley Playa, except for the AOD sensitivity test. This phenomenon may have been caused by the lower reflectance in Railroad Valley Playa. Therefore, the scattering effect caused by aerosol particles was reduced. The SBDE test in Alkali Lake and Railroad Valley Playa between FS-5 RSI and LS-8 OLI is shown in Figure 9. The relative errors caused from SBDE were approximately 1% in the MS bands and 4% in the PAN bands. The variations of SBDE in MS bands due to different atmospheric parameter contents were lower than 0.5%. Therefore, different atmospheric conditions had limited impacts on SBDE. On the other hand, the relative errors in PAN bands were significant because of the huge difference between the RSF of FS-5 RSI and L-8 OLI. This resulted in the underestimation of the atsensor radiance based on cross-calibration. The SBDE test in Alkali Lake and Railroad Valley Playa between FS-5 RSI and LS-8 OLI is shown in Figure 9. The relative errors caused from SBDE were approximately 1% in the MS bands and 4% in the PAN bands. The variations of SBDE in MS bands due to different atmospheric parameter contents were lower than 0.5%. Therefore, different atmospheric conditions had limited impacts on SBDE. On the other hand, the relative errors in PAN bands were significant because of the huge difference between the RSF of FS-5 RSI and L-8 OLI. This resulted in the underestimation of the at-sensor radiance based on cross-calibration. The environmental effect on at-sensor reflectance in different AOD550 and sizes of the calibration site are shown in Figure 10. It includes two cases: ROC > ROE and ROC < ROE. For both cases, the sizes of the calibration sites are important. The larger calibration site tended to have a lower environmental effect. As Figure 10 shows, the calibration site with a radius larger than 10 km can be considered a homogeneous surface. In addition, AOD550 played a vital role in environmental effects. The consistency of at-sensor reflectance between the different sizes of the calibration sites was determined by AOD550 loading. The higher AOD550 loading indicated a higher atmospheric effect. The environmental effect on at-sensor reflectance in different AOD550 and sizes of the calibration site are shown in Figure 10. It includes two cases: ROC > ROE and ROC < ROE. For both cases, the sizes of the calibration sites are important. The larger calibration site tended to have a lower environmental effect. As Figure 10 shows, the calibration site with a radius larger than 10 km can be considered a homogeneous surface. In addition, AOD550 played a vital role in environmental effects. The consistency of at-sensor reflectance between the different sizes of the calibration sites was determined by AOD550 loading. The higher AOD550 loading indicated a higher atmospheric effect. Remote Sens. 2019, 9, x FOR PEER REVIEW 14 of 20  Figure 11 displays the results of the vicarious calibration at two calibration sites, Alkali Lake and Railroad Valley Playa, in 2017 and 2018. The variation in calibration coefficient (K0) of each spectral band is evident, which could be mainly caused by the unstable optical system in different environments between the surface and space after launch. Most of the spectral K0s increased in 2017 (except Dec. 06), indicating the decreasing sensitivity to the optical property, which was similar to the behavior of FS-2 RSI in 2006 [19]. In 2018, the K0s tended toward stable and a slight increase excluding the red band. Table 5 summarizes the relative differences of K0 in 2018 for the calibration sites at Alkali Lake and Railroad Valley Playa. According to the previous study, the uncertainty of the absolute radiometric calibration was around ±5%, which was particularly somewhat worse right after launch [33]. Hence, there seemed to be great results for both sites. However, for panchromatic channels, the relative difference was up to 6%.

Vicarious Calibration
It is worth noting that some of the K0 values from Alkali Lake were relatively smaller than those of Railroad Valley Playa such as the results on Dec. 06, 2017 and Sep. 06, 2018, as shown in Figure 11. The main factor could be the obvious environmental effect caused by the very small bright playa with a dark surrounding in the region of Alkali Lake (see also Figure 2a). The radiance observed at the top of the atmosphere may be reduced because the distinct radiances between the bright playa scattered out and the dark surrounding scattered into, which could result in the smaller value and variation of  Figure 11 displays the results of the vicarious calibration at two calibration sites, Alkali Lake and Railroad Valley Playa, in 2017 and 2018. The variation in calibration coefficient (K 0 ) of each spectral band is evident, which could be mainly caused by the unstable optical system in different environments between the surface and space after launch. Most of the spectral K 0 s increased in 2017 (except Dec. 06), indicating the decreasing sensitivity to the optical property, which was similar to the behavior of FS-2 RSI in 2006 [19]. In 2018, the K 0 s tended toward stable and a slight increase excluding the red band. Table 5 summarizes the relative differences of K 0 in 2018 for the calibration sites at Alkali Lake and Railroad Valley Playa. According to the previous study, the uncertainty of the absolute radiometric calibration was around ±5%, which was particularly somewhat worse right after launch [33]. Hence, there seemed to be great results for both sites. However, for panchromatic channels, the relative difference was up to 6%.

Vicarious Calibration
It is worth noting that some of the K 0 values from Alkali Lake were relatively smaller than those of Railroad Valley Playa such as the results on Dec. 06, 2017 and Sep. 06, 2018, as shown in Figure 11. The main factor could be the obvious environmental effect caused by the very small bright playa with a dark surrounding in the region of Alkali Lake (see also Figure 2a). The radiance observed at the top of the atmosphere may be reduced because the distinct radiances between the bright playa scattered out and the dark surrounding scattered into, which could result in the smaller value and variation of K 0 s. Therefore, sincerely considering the environmental effect on the Alkali Lake calibration site is suggested. K0s. Therefore, sincerely considering the environmental effect on the Alkali Lake calibration site is suggested. Figure 11. The FS-5 RSI calibration coefficient (K0) of vicarious calibration starting at preflight. The black, blue, green, red, and magenta lines stand for the Pan, Blue, Green, Red, and NIR (Near-infrared) spectral bands, respectively. In Figure 12, most of the 0 values retrieved from cross-calibration had the same tendencies. The relative differences of 0 values from cross-calibration results are listed in Table 7. The variation of K0 may be due to different viewing angles, surface reflectance characteristics, and observation durations between the FS-5 RSI and LS-8 OLI satellites. In addition, both sensors had a large difference in bandwidth between the near-infrared and panchromatic bands; this difference in bandwidth may also have affected the calibration results. Figure 11. The FS-5 RSI calibration coefficient (K 0 ) of vicarious calibration starting at preflight. The black, blue, green, red, and magenta lines stand for the Pan, Blue, Green, Red, and NIR (Near-infrared) spectral bands, respectively.

Cross Calibration
In Figure 12, most of the K 0 . values retrieved from cross-calibration had the same tendencies. The relative differences of K 0 values from cross-calibration results are listed in Table 6. The variation of K 0 may be due to different viewing angles, surface reflectance characteristics, and observation durations between the FS-5 RSI and LS-8 OLI satellites. In addition, both sensors had a large difference in bandwidth between the near-infrared and panchromatic bands; this difference in bandwidth may also have affected the calibration results.  Figure 12. Same as Figure 11, but for the cross-calibration. The standard deviations of the calibration coefficients (K0) are summarized in Table 7. The spectral K0 retrieved from the vicarious calibration and cross calibration had the same tendency, although the magnitudes between the two methods were different. The higher standard deviation of vicarious calibration could have been affected from the error source of the ground reflectance measurement and non-Lambertian ground characteristic, which is suggested as 3.3% [46]. The results of radiometric calibration demonstrated the proposed procedures were acceptable and practical. In comparison to FS-2 RSI, the K0 changed suddenly in the early years, but remained stable for the later years. The phenomenon was the same as the other satellite sensor. The calibration coefficients will affect the radiometric quality of the image, thus periodical calibration is necessary and very important for further scientific research and applications.   Figure 11, but for the cross-calibration.
The standard deviations of the calibration coefficients (K 0 ) are summarized in Table 7. The spectral K 0 retrieved from the vicarious calibration and cross calibration had the same tendency, although the magnitudes between the two methods were different. The higher standard deviation of vicarious calibration could have been affected from the error source of the ground reflectance measurement and non-Lambertian ground characteristic, which is suggested as 3.3% [46]. The results of radiometric calibration demonstrated the proposed procedures were acceptable and practical. In comparison to FS-2 RSI, the K 0 changed suddenly in the early years, but remained stable for the later years. The phenomenon was the same as the other satellite sensor. The calibration coefficients will affect the radiometric quality of the image, thus periodical calibration is necessary and very important for further scientific research and applications.

Conclusions
In this study, two kinds of on-orbit absolute radiometric approaches-vicarious calibration and cross-calibration-were examined for a new optical sensor RSI onboard FORMOSAT-5. The consistent results indicate that the proposed approaches in on-orbit radiometric calibration are highly practical. According to the results, the findings of this study can be concluded below.

1.
Atmospheric variables are essential for the accurate radiometric calibration of on-orbit sensors including aerosol property, ozone, and water vapor, especially for the large viewing angle of the observations. For example, the effect of improper ozone content on radiometric calibration could reach 6%-10% uncertainty in the Red, Green, and Pan spectral bands (Figure 5a). Thus, on-orbit cross calibration is recommended rather than vicarious calibration when atmospheric measurements are lacking.

2.
Both vicarious and cross-calibration methods show a decaying optical sensitivity, which resulted in rapid changes (6%-24% in radiometric coefficient) during the first year after launch. Most calibration results remained stable, around 3%-5% relative difference, in the second year.
To maintain image quality, further research is suggested to keep track of FORMOSAT-5 RSI radiometric variation.