2.1. Dataset Description
The Dunhuang test site is in the Gobi desert, approximately 35 km west of the city of Dunhuang. The test site was an operational radiometric calibration and validation site for Chinese satellite sensors in 2001, located on the eastern edge of the Kumutage Penniform Desert in Gansu province, Southwest China. The whole target area for vicarious calibration is situated on a stabilized alluvial fan, 30 km × 30 km in size [
10]. The atmosphere is dry, clean, and typically has low levels of aerosol loading, making it beneficial for the calibration experiments; the site was chosen as one of the Committee on Earth Observation Satellites (CEOS) calibration and validation test sites.
The VIRR onboard FY-3C is a heritage instrument from the Multispectral Visible Infrared Scanning Radiometer (MVISR) onboard FY-1C and D. It provides images in 10 spectral bands between 0.44 and 12.5 μm, with a spatial resolution of 1.1 km at nadir, and it includes five visible-near infrared bands, two shortwave bands, and three middle and thermal infrared bands; its data records could be used for vegetation and ocean colour monitoring. VIIRS collects radiometric and imagery data in 22 spectral bands within the visible and infrared region ranging from 0.4 to 12.5 μm, including 16 moderate-resolution (750-m pixels) and five imagery resolution (375-m pixels) bands, plus one panchromatic “Day-Night Band”. The VIIRS spectral data are calibrated and geolocated in ground processing to generate Sensor Data Records (SDRs) [
11,
12]. The solar reflective bands of VIIRS, covering similar wavelength range as MODIS, are also calibrated by the solar diffuser (SD) and lunar observations, with a calibration accuracy of approximately 2% in most bands. Thus, bands M3, M4, I1, and I2 are used to cross-calibrate the corresponding VIRR bands 8, 9, 1, and 2; the characteristics of VIRR and VIIRS visible-near infrared channels are shown in
Table 1 and
Figure 1.
In addition, because the two sensors onboard different platforms have different overpassing times, only clear-sky scenes observed by both sensors are employed during 2014, so that the atmospheric effect resulting from different acquisition time could be reduced as much as possible, and no temporal matching is considered. In this study, a total of 11 scenes of VIRR images over Dunhuang site are acquired and used as the data source for cross-calibration with the corresponding VIIRS calibrated and geolocated SDR data. The dates and viewing geometries of these scenes are shown in
Table 2. Note that some of the pairs have large differences in viewing direction between the two sensors, and this difference could incur larger calibration errors without an accurate Bi-directional Reflectance Distribution Function (BRDF) correction. As an example,
Figure 2 shows the image of VIIRS band I2 (see
Figure 2a) and VIRR band 2 (see
Figure 2b) over the Dunhuang test site on 16 September 2014.
Table 1. The characteristics of VIRR and VIIRS visible-near infrared channels.
Table 1. The characteristics of VIRR and VIIRS visible-near infrared channels. | Band | Centre Wavelength (μm) | Spectral Range (μm) | Spatial Resolution at Nadir (m) |
---|
FY-3C/VIRR | 1 | 0.630 | 0.58–0.68 | 1100 |
2 | 0.865 | 0.84–0.89 | 1100 |
8 | 0.505 | 0.48–0.53 | 1100 |
9 | 0.555 | 0.53–0.58 | 1100 |
NPP/VIIRS | I1 | 0.640 | 0.60–0.68 | 375 |
I2 | 0.865 | 0.85–0.88 | 375 |
M3 | 0.488 | 0.478–0.498 | 750 |
M4 | 0.555 | 0.545–0.565 | 750 |
Figure 1. Spectral response functions of VIIRS and VIRR.
Figure 1. Spectral response functions of VIIRS and VIRR.
Table 2. The dates and viewing geometries of VIIRS and VIRR scenes.
Table 2. The dates and viewing geometries of VIIRS and VIRR scenes.Date | VIIRS | VIRR |
---|
| SZA | SAA | VZA | VAA | SZA | SAA | VZA | VAA |
---|
8 January 2014 | 62.68 | −173.55 | 37.31 | 74.33 | 64.50 | 162.11 | 21.07 | −77.11 |
13 January 2014 | 62.03 | −172.32 | 29.87 | 75.32 | 63.48 | 163.14 | 30.24 | −75.71 |
24 January 2014 | 59.56 | −174.62 | 37.35 | 74.32 | 61.78 | 160.39 | 23.50 | −76.73 |
29 January 2014 | 58.40 | −173.02 | 29.88 | 75.31 | 60.19 | 161.54 | 32.36 | −75.38 |
13 March 2014 | 43.43 | −171.67 | 37.61 | 74.22 | 45.60 | 157.12 | 29.94 | −75.79 |
14 March 2014 | 42.73 | −178.36 | 54.18 | 71.07 | 46.77 | 150.94 | 2.28 | −87.08 |
24 March 2014 | 39.05 | −171.91 | 43.72 | 73.23 | 41.62 | 154.53 | 22.92 | −76.86 |
6 May 2014 | 24.28 | −165.10 | 48.87 | 72.26 | 26.86 | 147.57 | 18.94 | −77.53 |
24 July 2014 | 21.72 | −156.56 | 29.69 | 75.37 | 22.27 | 152.56 | 47.02 | −72.89 |
11 September 2014 | 36.16 | −168.11 | 49.32 | 72.22 | 37.19 | 159.92 | 30.21 | −75.77 |
16 September 2014 | 38.40 | −165.41 | 43.84 | 73.27 | 38.56 | 163.97 | 37.89 | −74.55 |
Figure 2. The images of VIIRS band I2 (a) and VIRR band 2 (b) over the Dunhuang test site on 16 September 2014.
Figure 2. The images of VIIRS band I2 (a) and VIRR band 2 (b) over the Dunhuang test site on 16 September 2014.
2.2. Cross-Calibration Approach
The cross-calibration method involves comparison of the radiance/reflectance measured by the calibrated sensor with that measured by a well-calibrated sensor as a reference. This exercise can be reduced to spatiotemporal coincidences,
i.e., acquisition by both sensors at the same time and with the same viewing geometries (for example, the SNO method). Nevertheless, such coincidences are not very frequent when comparing two sensors with different orbits, altitudes, cycles, and local equatorial crossing time. In these conditions, the reference is, in general, not acquired at the same time for exactly the same spectral range and for the same viewing geometry [
13]. For this reason, some corrections must be applied to take into account these aspects.
To alleviate the impact of viewing geometry on the cross-calibration, the Bidirectional Reflectance Distribution Function (BRDF) characteristics of Dunhuang site measured in 2011 with SVC HR1024 spectrometer are used. Measurements were acquired with the viewing zenith angle scanning from 0° to 70° with a step of 14°, and the relative azimuth angles between the sun and viewing direction varied from 0° to 180°, with a step of 30°. Seven datasets of hemispherical scanning measurements were used with the solar zenith angle ranging from 29° to 52°. Examples of the multi-angle Bidirectional Reflectance Factor (BRF) measurements corresponding to VIIRS bands M3, M4, I1, and I2 are shown in
Figure 3. The figure shows a general increasing trend towards the backward scattering direction for these bands. In this study, first, with the atmospheric parameters at the VIIRS scenes’ acquisition time [the water vapor content (WVC) extracted from the National Centers for Environmental Prediction (NCEP) reanalysis data, the ozone content extracted from the product of Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and the assumed visibility of 40 km (see
Table 3), the surface reflectance of Dunhuang site is derived from VIIRS reflectance at the top of atmosphere (TOA) using the radiative transfer model 6S. Next, the angular effect of surface reflectance is corrected with the BRDF model proposed by Roujean
et al. [
14], fitted with the multi-angle BRF measurements, and the surface reflectance along with VIRR viewing direction could be acquired. Subsequently, based on the atmospheric parameters at the VIRR scene acquisition times, the corresponding TOA reflectance is simulated using the 6S model. In this study, the default solar model in the 6S radiative transfer model is used to characterize the solar irradiance [
15].
Figure 3. Measured Bidirectional Reflectance Factor of the Dunhuang site. (a) VIIRS M3; (b) VIIRS M4; (c) VIIRS I1; and (d) VIIRS I2.
Figure 3. Measured Bidirectional Reflectance Factor of the Dunhuang site. (a) VIIRS M3; (b) VIIRS M4; (c) VIIRS I1; and (d) VIIRS I2.
Table 3. Atmosphere parameters at VIIRS and VIRR acquisition times.
Table 3. Atmosphere parameters at VIIRS and VIRR acquisition times.Date | WVC at VIIRS Scenes Acquisition Time (g·cm−2) | WVC at VIRR Scenes Acquisition Time (g·cm−2) | Ozone (atm-cm) | Visibility (km) |
---|
8 January 2014 | 0.260 | 0.231 | 0.334 | 40.0 |
13 January 2014 | 0.269 | 0.171 | 0.335 |
24 January 2014 | 0.369 | 0.464 | 0.332 |
29 January 2014 | 0.461 | 0.526 | 0.327 |
13 March 2014 | 0.298 | 0.272 | 0.341 |
14 March 2014 | 0.635 | 0.619 | 0.341 |
24 March 2014 | 0.551 | 0.524 | 0.347 |
6 May 2014 | 0.751 | 0.754 | 0.280 |
24 July 2014 | 1.234 | 1.198 | 0.280 |
11 September 2014 | 0.848 | 0.844 | 0.297 |
16 September 2014 | 0.397 | 0.318 | 0.297 |
Different applications and technological developments in Earth observation necessarily require different spectral coverage [
16]. Thus, spectral bands differ significantly among sensors, even for bands designed to observe at the same region of the electromagnetic spectrum; as a result, these sensors yield fundamentally different measurements that are not directly comparable. To remove the effect of spectral characteristics on cross-calibration, the spectral band adjustment factor
k in a given spectral band
i is calculated as a ratio of the TOA reflectance from two sensors with a simulation method. In this method, a series of TOA reflectance for VIIRS and VIRR sensors in a given spectral bands
i are calculated using the 6S model. The band surface reflectances, varying from 0.1 to 0.5 with a step of 0.1, are used as inputs to drive the 6S model, together with the selected atmospheric states, the sensor’s spectral response function, and the illumination and viewing geometries. Subsequently, the corresponding spectral band adjustment factor
is fitted by linear regression method.
After alleviating the impact of viewing geometry and spectral characteristics, the measured TOA reflectance of VIRR in a given band
i could be evaluated with the simulated VIRR reflectance values. The flowchart is shown in
Figure 4:
Figure 4. The flowchart of the cross-calibration procedure.
Figure 4. The flowchart of the cross-calibration procedure.