1. Introduction
Radiometric calibration is to convert the digital number (DN) value of a satellite sensor to a physical characterization, such as apparent radiance or apparent reflectance, in the visible and near-infrared sensors. Radiometric calibration is a critical step and the basis of the quantification of remote sensing application. The radiometric calibration not only characterizes the status and degradation with the time of the satellite sensors after launch, but also establishes the physical relationship between multi-source remote sensing data and expands the data source of quantitative remote sensing. Therefore, high precision radiometric calibration plays an important role in the development of quantitative remote sensing.
The calibration techniques can be divided into prelaunch calibration and on-orbit (post-launch) calibration. There are five types of on-orbit calibration methods—inter-calibration [
1,
2], solar/lunar calibration [
3,
4], vicarious calibration [
5,
6], scene calibration [
7,
8,
9], and cross-calibration [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19]. Cross-calibration is an important and widely used on-orbit calibration method, which can be used for most optical sensors. The cross-calibration uses high-precision reference satellite data to achieve radiometric calibration of the sensors to be calibrated; thus, it can be used with high frequency and low cost. Furthermore, cross-calibration could be used for historical image.
The research on cross-calibration is going on for more than 20 years. Teillet et al. described a methodology that used airborne hyperspectral imagery as reference data to calibrate multiple satellite sensors [
11]. Thome et al. used the Railroad Valley Playa test site in Nevada to cross-calibrate the Earth Observation-1 Advanced Land Imager (EO-1 ALI), Terra MODIS, and IKONOS with respect to the Landsat 7 Enhanced Thematic Mapper Plus (ETM+) [
12]. Chander et al. presented the results from the cross-calibration of the ETM+ and ALI image pairs using two approaches [
13]. One approach was based on image statistics, and the other was based on vicarious calibration. Thome et al. performed cross-calibrations between Landsat 5 TM and Landsat 7 ETM+ on multiple dates [
14]. Teillet et al. presented the results of TM and ETM+ during the tandem-orbit configuration in 1999 [
15]. Eight tandem image pairs were processed in cross-calibration. Angal et al. evaluated the radiometric calibration agreement between MODIS and ETM+ over an African pseudo-invariant calibration site, Libya 4. The uncertainties of the BRDF, spectral response difference, and water vapor column were accounted for in this cross-calibration [
16].
The cross-calibration method is also used for Chinese satellites by researchers, in virtue of its unique characteristics and precision of the results. GF-1 satellite is the first satellite of the Chinese high-resolution satellite constellation, launched on 26 April 2013, and equipped with two 2 m resolution panchromatic/8 m resolution Multispectral (PMS) cameras and four 16 m resolution wide-field-of-view (WFV) cameras. The China Centre for Resource Satellite Data and Application (CRESDA, the official operational management agency of GF-1 satellite) is responsible for the calibration of GF-1 using site calibration method and the publication of calibration coefficients. Since the launch, many researchers have invested substantial effort in the radiometric cross-calibration of GF-1 sensors. Han et al. (2014b) carried a cross calibration between GF-1 WFV sensor, Shijian 9A satellite (SJ-9A) PMS, and Terra MODIS sensor [
17]. Yang L et al. (2015) cross calibrated the GF-1 satellite with Landsat 8 OLI sensor image [
18]. Li et al. (2016) proposed a cross-calibration method to cross calibrate the GF-1 satellite WFV sensor, taking Landsat 8 OLI sensor as the reference [
19]. The USGS spectral library was used for spectral band adjustment factors (SBAF) correction.
The sensors with similar spatial resolution are reasonably selected as reference for the radiometric cross-calibration of high-resolution, narrow-width satellites, such as GF-1 PMS and Landsat 8 OLI. However, due to the long revisit period and low frequency of satellites passing through the calibration site, it is difficult to obtain reference images with the same imaging time and nadir observation as the satellites to be cross-calibrated.
Under these circumstances, there are two ways to accomplish on-orbit cross-calibration for high-resolution and narrow-swath sensors. The first one is to use satellite images with the same imaging date but large view-angle and low spatial resolution as reference. The other is to select the satellite images with different imaging dates but nadir observation as reference sensor, such as Landsat 8 OLI. Gao et al. carried out the research on this issue. Landsat 8 and MODIS were both selected as reference sensors to cross-calibrate the GF-1 PMS sensor at Dunhuang and Golmud test sites, respectively [
20]. The results show that the calibration accuracy based on Landsat 8 is better than that based on MODIS calibration, which indicates that, at large observation angles, the error caused by the difference of observation angle is larger than that caused by the difference of imaging date.
However, in this paper, both the two test sites were assumed to be Lambertian, without considering the difference caused by the change of surface directional reflection. Some researchers have developed some new cross-calibration methods by BRDF (Bidirectional Reflectance Distribution Function) correction. As to the cross-calibration of GF-1 satellite, Feng et al. has taken Landsat-8 OLI as the reference sensor to cross calibrate the GF-1 satellite WFV sensor [
21]. The large observation angle difference between Landsat 8 OLI image and GF-1 WFV images were corrected using Terra MODIS bi-directional reflectance distribution function (BRDF) products. Zhong et al. developed a new cross-calibration technique for HJ-1satellite CCD (Charge-coupled Device) sensor [
22]. In this method, the Landsat ETM+ and Advanced Space-borne Thermal Emission and Reflection Radiometer global digital elevation model (ASTER GDEM) product are used to develop a BRDF model of a desert site. Yang A. et al. updated the method proposed by Zhong et al. and used for the GF-1 satellite WFV sensor [
23]. The Landsat 8 and the DEM extracted by the three-line camera sensor (TLC) onboard ZiYuan-3 are replaced to develop BRDF model of a desert site. Yang A. et al. also applied this technique to calibrate the GF-4 satellite PMS sensor [
24]. It can be seen that the BRDF models in prior studies are based on ground measured data or existing BRDF products. Moreover, the MODIS BRDF product and ground reflectance BRDF model are based on the given solar angle, which has some drawbacks and limitations. Therefore, as the absence of ground reflectance and without using BRDF products, how to construct a TOA BRDF model of calibration site and correct the TOA reflectance with a large observation angle (such as MODIS images) to verify the accuracy of calibration results is the main purpose of this paper.
The structure of this paper is as follows: GF-1 PMS, MODIS, Landsat 8 OLI, the test sites, and datasets are described in
Section 2.
Section 3 introduces the methods of this paper, including BRDF model construction method and radiometric cross-calibration method. The BRDF model construction results, the radiometric cross-calibration results based on BRDF model and the verification of the accuracy of the method proposed in this paper are illustrated in
Section 4.
Section 5 discusses the uncertainty caused by ground measured data, BRDF model, etc., and the total uncertainty of cross-calibration using MODIS as reference sensor are given. The conclusions of our research are made in
Section 6.
6. Conclusions
An innovative cross-radiation calibration method is proposed in this paper, which is used to cross-calibrate GF-1 PMS camera at Dunhuang and Golmud test site, respectively, taking MODIS and Landsat 8 OLI as reference sensors.
Firstly, a BRDF model is constructed based on the TOA reflectance of MODIS in long time-series sunny day, which are extracted by the TOA BT and the variation coefficient of TOA reflectance. This process is an important step of this paper. Then the TOA reflectance of MODIS images with large imaging angles is corrected by BRDF model. After that, the radiometric cross-calibration of PMS is carried out by using BRDF-corrected MODIS and OLI images, and the radiometric calibration coefficients are calculated by the least square method.
Compared with previous studies, it is found that the accuracy of cross-radiation calibration has been significantly improved after directional correction of large observation angle MODIS images using the BRDF model constructed in this paper. Furthermore, the study also finds that although ground measured parameters of Dunhuang test site are used as substitution for the Golmud test site in cross-calibration, the calibration coefficients obtained are in good agreement with that of the Dunhuang site. By comparing with the site calibration coefficients and analyzing the uncertainty of this paper, the total calibration uncertainties of PMS using MODIS as reference sensor are about 2%–6%, which are similar with the accuracy of vicarious calibration. The result is promising, and in the future, the robustness and universality of the technique proposed in this paper should also be assessed by applying it to other radiometric calibration test sites and similar optical satellite sensors. In addition, as CRESDA has adjusted the gain of GF-1 PMS1 on 25 September 2016, this paper only studies the PMS1 data from the launch of GF-1 satellite on 26 April 2013 to 25 September 2016. The radiation performance and calibration accuracy of the gain-adjusted GF-1 PMS1 sensor will be studied in the follow-up study.