Next Article in Journal
Coupling Random Forest, Allometric Scaling, and Cellular Automata to Predict the Evolution of LULC under Various Shared Socioeconomic Pathways
Previous Article in Journal
Mapping Tree Species Using CNN from Bi-Seasonal High-Resolution Drone Optic and LiDAR Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Development of Dark Hyperspectral Absolute Calibration Model Using Extended Pseudo Invariant Calibration Sites at a Global Scale: Dark EPICS-Global

by
Padam Bahadur Karki
1,2,
Morakot Kaewmanee
2,
Larry Leigh
2,* and
Cibele Teixeira Pinto
2
1
Department of Electrical Engineering and Computer Science, South Dakota State University (SDSU), Brookings, SD 57007, USA
2
Image Processing Laboratory, Department of Electrical Engineering and Computer Science, South Dakota State University (SDSU), Brookings, SD 57006, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(8), 2141; https://doi.org/10.3390/rs15082141
Submission received: 3 March 2023 / Revised: 31 March 2023 / Accepted: 13 April 2023 / Published: 18 April 2023

Abstract

:
This research aimed to develop a novel dark hyperspectral absolute calibration (DAHAC) model using stable dark targets of “Global Cluster-36” (GC-36), one of the clusters from the “300 Class Global Classification”. The stable dark sites were identified from GC-36 called “Dark EPICS-Global” covering the surface types viz. dark rock, volcanic area, and dark sand. The Dark EPICS-Global shows a temporal variation of 0.02 unit reflectance. This work used the Landsat-8 (L8) Operational Land Imager (OLI), Sentinel-2A (S2A) Multispectral Instrument (MSI), and Earth Observing One (EO-1) Hyperion data for the DAHAC model development, where well-calibrated L8 and S2A were used as the reference sensors, while EO-1 Hyperion with a 10 nm spectral resolution was used as a hyperspectral library. The dark hyperspectral dataset (DaHD) was generated by combining the normalized hyperspectral profile of L8 and S2A for the DAHAC model development. The DAHAC model developed in this study takes into account the solar zenith and azimuth angles, as well as the view zenith and azimuth angles in Cartesian coordinates form. This model is capable of predicting TOA reflectance in all existing spectral bands of any sensor. The DAHAC model was then validated with the Landsat-7 (L7), Landsat-9 (L9), and Sentinel-2B (S2B) satellites from their launch dates to March 2022. These satellite sensors vary in terms of their spectral resolution, equatorial crossing time, spatial resolution, etc. The comparison between the DAHAC model and satellite measurements showed an accuracy within 0.01 unit reflectance across the overall spectral band. The proposed DAHAC model uncertainty level was determined using Monte Carlo simulation and found to be 0.04 and 0.05 unit reflectance for the VNIR and SWIR channels, respectively. The DAHAC model double ratio was used as a tool to perform the inter-comparison between two satellites. The sensor inter-comparison results for L8 and L9 showed a 2% difference and 1% for S2A and S2B across all spectral bands.

Graphical Abstract

1. Introduction

Satellite imagery is a means to uncover the Earth’s changes; however, it is critical to have credible data for scientific study, which are obtained by optical satellite sensors. A wide range of applications, including atmospheric physics and geoscience, make use of physical quantities converted from the digital data recorded from satellite images. The reflection or emission of radiation by the Earth’s surface or atmosphere is being used by satellites, allowing the retrieval of the corresponding physical quantities. Although voltage or recorded digital data serve as the sensors’ primary measuring quantity, calibration must be carried out to compare the sensor-derived digital data with incoming radiance in their physical interest. At present, satellite instruments are typically well-designed and calibrated before launch. Though the satellite equipment is advanced and robust, they eventually deteriorate in space owing to factors including temperature, mechanical, electrical, or UV radiation exposure [1]. In order to ensure the quality of the derived variables and products, reliable remote sensing relies on a sensor re-calibration, commonly called radiometric calibration [2].

1.1. Absolute Radiometric Calibration

Data evaluation and the quality of satellite data highly depend on the absolute radiometric calibration of the sensor. During the absolute calibration process, the satellite images are processed and converted into digital numbers or the measured voltage in the form of physical quantities such as at-sensor spectral radiance or the top of atmosphere (TOA) reflectance [3,4]. The absolute calibration approach helps in characterizing and measuring the sensor’s performance from the early pre-launch stage to its on-orbit operation. Calibration can also be performed using different methods, for example pre-launch, onboard, and vicarious calibration. The results of post-launch calibration are compared with pre-launch calibration using several different methods in the lab [5]. After the satellite is placed in its orbit, the sensor characterization and its accuracy measurement are needed throughout its lifetime. Onboard calibrators, such as solar diffuser panels, lamps, etc., can be used to perform radiometric calibration [6,7]. However, onboard calibrators require routine operation and maintenance. They are expensive and also prone to the effects of harsh conditions in the space environment. For these reasons, many satellite sensors do not include onboard calibration. In order to overcome these limitations, pseudo invariant calibration sites (PICSs) were discovered for satellite calibration. The satellite imagery over the PICSs across the Earth’s surface is spatially, temporally, and spectrally stable [8].

1.2. Stable Calibration Sites

Several studies on PICS-based assessment have been performed to identify the stable pixels for sensor radiometric stability for more than 20 years [3]. Spatially uniform PICSs exhibit stable spectral features over time, higher reflectance, and no atmospheric influence on upward radiation [9]. A spatial uniformity of at least 3% and a temporal variability of 1–2% was found by Cosnefroy et al. over twenty desert locations [10]. There are six PICSs for North Africa that the Committee on Earth Observation Satellites (CEOS) recommended showing temporal variability in all bands of 3% or less [11]. According to Helder et al., Libya 4, Libya 1, Algeria 3, Arabia 2, Egypt 2, and Egypt 1 had variabilities of less than 3% and were capable of monitoring long-term trends [8]. Research studies on identifying stable pixels were then conducted by the South Dakota State University Image Processing Laboratory (SDSU IP Lab) to identify the best locations with 3% or less of temporal, spatial, and spectral variability. Less temporal uncertainty was considered in Libya 4, Niger 1, Sudan 1, Niger 2, Egypt 1, and Libya 1. With a temporal variation of less than 3%, Libya 4 stood out as the most-consistently stable site among them and could be widely employed in radiometric calibration work [3,12,13,14]. Shrestha et al. [6] applied the K-means unsupervised classification method, which resulted in identifying 19 clusters representing distinct land surface types in North Africa with pixels having 5% or less of both the spatial and temporal uncertainty. Fajardo et al. expanded the unsupervised classification to 160 clusters at a global scale and identified a cluster with stable bright sites named “global EPICS” [5]. The global EPICS showed a temporal variation of less than 4% for all bands, with some as low as 2.7%. The global classification with 160 clusters was further improved by the author to develop the “300 Class Global Classification” with its 300 clusters around the globe [15]. Past studies, including the research at SDSU IP Lab, focused on absolute calibration using stable bright sites, i.e., PICSs. However, the absolute calibration model developed in the past works neither covers the dynamic range of reflectance measurements nor developed the calibration model using stable dark targets at a global scale. This work proposes an algorithm to build an absolute calibration model using low reflectance dark targets capable of performing satellite calibration in addition to a bright target absolute calibration model. In this study, low reflectance measurements were identified from the stable dark targets of “Global Cluster-36” (GC-36), one of the clusters from the “300 Class Global Classification”. This study used the identified stable dark targets of the GC-36, named “Dark EPICS-Global”.

1.3. Evolution on Development of Absolute Calibration Model

Govaerts et al. [16] developed an absolute calibration model using the PICSs for geostationary satellite sensors. To characterize the atmospheric effects more accurately, an advanced radiative transfer model was developed for the Libya 4 site in 2012. The expanded model included the polarization effects and non-spherical aerosol models and showed accuracy within 3% [17]. In 2010, Helder et al. developed the concept of an empirical absolute calibration model using PICSs and the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) as a reference sensor. The model was validated with Landsat-7 (L7) with an accuracy of 3% and 6% for the visible and short-wave infrared (SWIR) channels, respectively [12]. Based on daily radiance observations over Libya 4 using a geostationary Earth orbit (GEO) sensor, Bhatt et al. developed a desert daily exoatmospheric radiance model (DERM) [18]. The model accuracy remained consistent within 0.4% and 1.9% for Meteosat-8 and Meteosat-7 with the reference Meteosat-9 DERM. Likewise, GOES-11 DERM showed an accuracy of 1% and 3% while calibrating GOES-10 and GOES-15, respectively [19]. In order to minimize the bidirectional reflectance distribution function (BRDF) impact generated by the sensor’s off-nadir observations and seasonal fluctuations owing to solar position shift, Mishra et al. [4] considered a view zenith angle (VZA) for the absolute calibration model. Using Terra MODIS as a reference sensor, this model achieved 3% accuracy and uncertainty within 2% for six spectral bands [4]. In 2017, the coastal aerosol band was included as part of an improvement to the existing absolute calibration model using data from Landsat-8 (L8) images named the refined absolute calibration model. The accuracy of the refined absolute calibration model was within 3% for all spectral bands [20]. Raut et al. [3] expanded the absolute calibration model considering five additional Saharan Desert PICSs and found the model performance within 3% accuracy with 2% precision for three sites viz. Egypt 1, Libya 1, and Sudan 1. However, due to a lack of insufficient reliable Hyperion data, the model performance showed less accuracy for Niger 1 and Niger 2 [3].
In the previous studies, the scene-center-specific spherical angle was used to develop a BRDF model to finalize the absolute calibration model [3,4,12]. The reprocessing of the Landsat archive was performed in 2017, including the geometric error of the 12 m root-mean-squared error and 3% radiometric uncertainty [21]. All Landsat Collection-1 data produced by this method contain files including the Sun illumination and sensor viewing angle coefficients, as well as details about the quality assessment (QA) bands. Using the existing information in the Level-1 product of Landsat-8, Farhad et al. developed a novel four-angle BRDF model using the solar and view geometry angles in the form of Cartesian coordinates. This four-angle BRDF model preserves the data nature to achieve a robust fit. The estimated temporal variation across all the spectral bands was within 1.8% over Libya 4 PICSs [22]. In 2019, the Extended PICS Absolute Calibration (ExPAC) model with solar and view angles was developed using extended pseudo invariant calibration site in North African desert sites (EPICS-NA) as a target and L8 as a reference sensor. The ExPAC model showed a prediction accuracy of 2% after validation with the L7, L8, Sentinel-2A (S2A), and Sentinel-2B (S2B) sensors [23]. Chaity et al. developed an empirical hyperspectral absolute calibration model for the Libya 4 PICS (hyperspectral APICS) using the L8 Operational Land Imager (OLI) and Earth Observing-1 (EO-1) Hyperion sensors [1]. To better represent the angles, the hyperspectral APICS model used pixelwise four-angle BRDF information in Cartesian coordinates. This model had the potential to perform absolute calibration in 1 nm spectral resolution and showed an accuracy and precision of 6% and 4%, respectively, for the off-nadir viewing sensor. Similarly, the model showed an accuracy and precision of 3% for the nadir viewing sensor, respectively [1]. However, the hyperspectral APICS model was constrained for the Libya-4 PICS and showed poor performance for non-Landsat bands, although this could be improved.

1.4. Objectives of the Study

The main objective of this study was to develop a novel dark hyperspectral absolute calibration (DAHAC) model using the L8, S2A, and EO-1 Hyperion measurements over the identified stable dark sites (Dark EPICS-Global) of “Global Cluster-36” (GC-36). L8 OLI and S2A MSI were used as reference sensors and EO-1 Hyperion as the hyperspectral library to convert the multispectral profile of L8 and S2A in the hyperspectral domain. The dark hyperspectral data and their angular information in the Cartesian coordinate form were used to build the four-angle hyperspectral BRDF model, which was further simplified to develop the DAHAC model. Unlike Chaity’s development, the proposed model is unique with its novel algorithm to develop a more accurate model using the “low reflectance measurements” at a global scale. The DAHAC model is capable of predicting the top of atmosphere (TOA) reflectance in all existing hyperspectral bands for any sensors. The DAHAC model’s validation was performed with respect to L7, S2B, and Landsat-9 OLI (L9). Furthermore, the DAHAC model double ratio was estimated to perform sensor inter-comparison between L8 and L9 for Landsats and S2A and S2B for Sentinel-2.
This paper is organized into the following sections: Section 1 illustrates the background and literature review on radiometric calibration along with the absolute calibration model development techniques based on PICSs and its evolution towards EPICS Global. Section 2 outlines the sensor characteristics used in this study, the selection of dark targets for the study area and its validation, the dark hyperspectral dataset generation, and the DAHAC model development process. Section 3 discusses the DAHAC model validation results with respect to several sensors. Section 4 presents the DAHAC model uncertainty estimation process. Section 5 describes the sensor inter-comparison based on the DAHAC model double ratio estimation. Finally, Section 6 summarizes and concludes the paper.

2. Methodology

This section explains the development of the dark hyperspectral absolute calibration (DAHAC) model using reflectance measurements of L8, S2A, and EO-1 over Dark EPICS-Global. Figure 1 shows the overall architecture of the study area selection and sensor data pre-processing for generating the dark hyperspectral dataset required for the DAHAC model’s development. This work started with a thorough analysis of different satellite sensors viz. L7, L8, L9, S2A, S2B, and EO-1 Hyperion. The multispectral sensors (L8 and S2A) were utilized as the reference sensors, and EO-1 was used as a hyperspectral library source. The preliminary selection and identification of many dark targets around the globe from 300 global clusters was the first major step. GC-36 was selected as the most-suitable dark target as its spectral response is closely similar to the volcanic site used for the cross-calibration of S2A and L8 [22]. In this study, the selected dark sites from GC-36 were validated by comparing their TOA reflectance with the overall GC-36 TOA reflectance measurements to provide Dark EPICS-Global. Furthermore, this section also explains the detailed process to generate the normalized hyperspectral profile for two multispectral sensors, L8 and S2A, using EO-1 Hyperion as a hyperspectral library. The procedures to generate the dark hyperspectral dataset for creating the 4-angle BRDF model are also presented in detail. Finally, the process of developing the DAHAC model from the 4-angle BRDF model is also discussed in the sections below.

2.1. Satellite Data

2.1.1. Landsat-7, -8, -9

Since 1972, Landsat missions have continuously collected multispectral data covering the entire planet with a temporal resolution of 16 days. L7, which was launched on 15 April 1999, has eight spectral bands and operates at a mean altitude of 705 km in a Sun-synchronous orbit. Prior to the launch of L8, L7 had a satellite calibration of 5%, making it the most-stable sensor within the Landsat family [7]. However, Scan Line Corrector (SLC) issues with L7 have existed since May 2003 and have resulted in scenes with wedge-shaped data gaps [24]. The radiometric stability of the ETM+ satellite has been evaluated with the help of on-orbit calibration and vicarious measurements. The partial- and full-aperture solar calibration is continuously performed by taking the source as the Sun for onboard calibrators [25]. The official science mission of L7 ended on 6 April 2022. L7 started its extended science mission on 5 May 2022, from a lower orbit of 697 km [26].
L8 is part of the Landsat series, which was launched on 11 February 2013. It has two different instruments: an OLI sensor and the Thermal Infrared Sensor (TIRS). L8 has been producing high-quality image information for Earth observation with spatial resolutions of 30 m for eight separate spectral bands and 15 m for the panchromatic band. The uncertainty associated with the L8 post-launch reflectance calibration is within 2% [7]. The improvement in the Landsat image product was performed in 2018, which provided a Level-1 product with angle information and the quality assessment band. The radiometric calibration was updated in the coastal aerosol and blue band for the Level-1 Collection-2 product [27].
The most recent satellite in the Landsat constellation is Landsat 9, which was launched on 27 September 2021. Landsat 9 is an identical instrument to Landsat 8 with better radiometric resolution. The increased quantization from 12 to 14 bits in L9 enables the sensor to acquire more small variations, resulting in better image quality. The Sun-synchronous orbit of L9 is at an altitude of 705 km and takes 16 days to complete one orbit cycle [28].
Equation (1) shows the conversion of Landsat imagery digital number (DN) values to TOA reflectance [29].
ρ o b s e r v e d = M ρ × Q C a l + A ρ c o s α
where M ρ and A ρ are the multiplicative and additive factors in metadata format, Q C a l represents the calibrated DN value, and α is the solar zenith angle.

2.1.2. Sentinel-2A, -2B

The European Space Agency (ESA) launched S2A on 23 June 2015 and S2B on 7 March 2017 under the Copernicus program. S2A and S2B are positioned in a Sun-synchronous orbit at an altitude of 786 km, phased at 180° relative to one another. These sensors are equipped with push-broom sensor multi-spectral instruments (MSIs) that measure the solar reflectance for 13 different spectral bands with 10m, 20m, and 60m spatial resolutions. The orbital period and combined orbital period for these two satellites are ten days and five days, respectively. The MSI focal plane detectors are divided into 12 different modules to allow these sensors to provide an image with a 290 km swath width at a 20.6° field of view [30,31]. In 2017, the study on the Sentinel-2 satellite showed that the absolute calibration was better than 5% for the overall sensor spectral band [32]. The European Space Agency reported that Sentinel-2 products would soon be upgraded. Sentinel-2A and -2B have undergone multiple processing systems during the course of their lifetimes. The processed data with a baseline of 2.0 or above, from the launch time to October 2022, were used in this study. The TOA reflectance is calculated using Equations (2) and (3).
ρ λ = D N c a l Q
where D N c a l represents the calibrated DN value and Q is a scale factor that is equal to 10,000. This scaling factor accounts for the Earth–Sun distance, exoatmospheric irradiation, and cosine correction. The Sentinel-2 Processing System Version 4.0 has been installed since 25 January 2022. Thus, the TOA reflectance conversion equation is updated as shown in Equation (3).
ρ λ = D N c a l + O f f s e t Q
where offset = −1000 is reported in the metadata image information. The estimated TOA reflectance for S2A is illustrated in Section 2.2.4.

2.1.3. Earth Observing-1 Hyperion

EO-1 Hyperion was a hyperspectral satellite that was part of NASA’s New Millennium Program. This satellite was launched on 21 November 2000 and was discontinued on 20 March 2017 [33]. There were two instruments on board: the Advanced Land Imager (ALI), a hyperspectral imaging spectrometer, and the Linear Etalon Imaging Spectral Array (LEISA) Atmospheric Corrector (LAC). The push-broom hyperspectral sensor covers the 400–2500 nm spectral range and 242 band images (including 196 onboard calibrated bands [34]) with a 10 nm spectral resolution and a 30 m spatial resolution over a 7.7 km swath width. The EO-1 stability analysis for 16 years of operation with varying orbital precession revealed calibration uncertainty for all spectral bands within 5 and 10% [35]. The following Equation (4) is used to convert the EO-1 Hyperion DN values to TOA reflectance.
ρ λ = D N c a l h × π × d 2 E S u n × s i n ϕ × c o s θ
where D N c a l represents the calibrated digital number value, h defines the scale factor, d indicates Earth–Sun distance in A.U., E S u n demonstrates the conversion of calibrated radiance to reflectance derived from the ChKur solar spectrum (ESUN(ChKur)) [36], and ϕ and θ are the Sun elevation and sensor look angle, respectively. The summarized information for the sensors used in this study is shown in Table 1.

2.2. Dark Target Selection

This section explains the steps for identifying dark targets by choosing the stable pixels from the dark target clusters under the “300 Class Global Classification” using the L8 imagery information. This section also describes the zonal mask creation process and filtering technique.

2.2.1. Identifying Stable Dark Target Using Landsat-8

Three different global clusters, Global Cluster-6 (GC-6), Global Cluster-26 (GC-26), and Global Cluster-36 (GC-36), were discovered to be the dark clusters among the 300 clusters from the “300 Class Global Classification”. These clusters include surface types viz. dark rock, dark sand, volcanic area, dark pond, etc. GC-36 was selected as the dark cluster for this study based on the higher pixel counts and pixel density. Since GC-36 is a GeoTIFF format, the Geospatial Data Abstraction Library (GDAL) is used to generate Keyhole Markup Language (KML) files for GC-36. Using a kmL file, L8 WRS-2 system, Google Earth Pro, and Earth Explorer program, an initial visual inspection was performed to choose the WRS-2 paths/rows with the surface type: dark rock, dark sand, volcanic area, and vegetation area on each continent. Initially, 13 distinct WRS-2 paths/rows from the Middle East, South Africa, and North Africa were chosen. However, the surface types viz. dark rock, dark sand, and volcanic area showed a similar spectral response as the Libyan volcanic site that was used for cross-calibration purposes between the L8 and S2A sensors [22].
Based on the spectral characteristics, eight different WRS-2 paths/row were chosen as dark sites for this study. The eight selected WRS-2 paths/rows have pixel counts within the range of 21,726–2,359,755. The 8 WRS-2 paths/rows in the three continents with surface type description are illustrated in Table 2. The selected WRS-2 paths/rows for L8, along with the global map are shown in Figure 2. They are further used to perform temporal stability analysis in Section 2.2.4.

2.2.2. Creation of Zonal Mask for Satellites

Following the selected dark pixels’ location information in Section 2.2.1, a zonal dark pixel mask was generated using the “gdalwarp” function from the GDAL library based on the spatial resolutions of 10 m, 20 m, 30 m, and 60 m and the Universal Transverse Mercator (UTM) zone. To accommodate images positioned in two distinct UTM zones, the zonal mask was extended by 100 km in the east–west direction.

2.2.3. Filtering and Dark Pixel Validation

The primary goal of applying filters is to obtain a clear pixel by removing the clouds from images. The Landsat products with Level 1 Collection 2 (L1C2) are provided with a pixel quality assessment band (BQA) [37]. Based on the BQA information, a binary cloud mask is built, including fill values, dilated cloud, cirrus, cloud, cloud shadow, cloud confidence, and cloud shadow confidence, which is then applied for each L7, L8, and L9 satellite image. In the case of the Sentinel-2 sensor, the provided binary cloud mask is applied for each Sentinel-2 image. Similarly, to remove the cloudy pixels from Hyperion images, Band-123 (1376 nm) image information with almost a 100% absorption feature is used as a binary cloud mask.
Although these cloud binary masks can be employed at the pixel level, we chose to use them first at the scene level to eliminate dates where the region of interest (GC-36 pixels) is substantially obscured by clouds. This is because cloud filters are not perfect and can miss cloud pixels on days that are extremely cloudy. To remove the cloudy image, a potential threshold containing a clear pixel for each image was determined at 40%. In other words, if the number of clear pixels in the cloud binary mask is larger than 40% of the number of pixels in the satellite images, the image is selected for applying the pixel-by-pixel cloud mask; otherwise, the entire scene is discarded.
Even after filtering the cloud in the sensor images, there could be a chance that the cloud shadow pixel has the same pixel values as the clear dark pixel. The clear dark pixels in the sensor images were validated using reference L8 image information. First, the BQA filter with cloud shadow and cloud shadow confidence was employed on the L8 images. The L8 images were also applied to the BQA filter without cloud shadow and cloud shadow confidence. Further, the convolution theory was applied to the resulting two sets of images in order to distinguish the real dark pixels in the images. This process clarified that the pixels associated with each image were clear dark pixels.

2.2.4. Dark Target Data for Absolute Calibration

The cloud-free images were used to estimate the TOA reflectance for each sensor. The estimated TOA reflectance for L8 (1150 images) for dark sites ranges from 0.043 to 0.17, as shown in Figure 3. It can be observed that the selected dark sites were temporally stable with a standard deviation of 0.02 across all the spectral bands. The selected sites (paths/rows) for L8 were validated with respect to GC-36 by comparing the estimated TOA reflectance. GC-36 itself is defined in a way that the pixels included within this class represent identical attributes based on unsupervised K-means clustering. The mean and standard deviation of TOA reflectance for each of the selected sites for each band was compared against the mean and 3 σ standard deviation of GC-36, as shown in Figure 4. The minimum and maximum 3 σ values for GC-36 are shown by red-colored dotted lines. In Figure 4, it can be observed that, for all spectral bands, the estimated mean TOA reflectance for each site is within the boundary line, except for the SWIR-2 band crossing the boundary line within an acceptable range. Therefore, the selected dark sites under GC-36 for the L8 sensor were validated and represented as “Dark EPICS-Global” for this study.
After identifying Dark EPICS-Global using L8 as a reference sensor, the locations of these selected dark sites were further used to investigate the availability of data for the S2A and EO-1 sensors. There are 12 tiles representing the dark sites selected for the S2A sensor. Figure 5 shows the estimated TOA reflectance using 775 images obtained from the 12 stable dark tiles. The TOA reflectance for each tile was within the range of 0.04–0.18, with a standard deviation of 0.02 for all spectral bands. Likewise, among the selected dark sites, three distinct WRS paths/rows, i.e., 168/50, 183/46, and 184/43 were available for the EO-1 sensor. A total of 64 Hyperion images were obtained from the selected stable sites, and their estimated hyperspectral TOA reflectance for the 196 calibrated bands are shown in Figure 6. The average TOA reflectance of all 64 EO-1 Hyperion images with their standard deviation is shown in Figure A1 presented in Appendix B. The estimated TOA reflectance using Dark EPICS-Global for L8, S2A, and EO-1 is called dark target data. The estimated dark target data were further used to generate the dark hyperspectral dataset (DaHD).

2.3. Dark Hyperspectral Dataset

This section explains the dark hyperspectral data generation process to develop the hyperspectral absolute calibration model. L8 and S2A were used as reference sensors since they are well-calibrated and their combination covers a wide range of angular information, as well as non-Landsat spectral bands. L8 and S2A have their TOA reflectance measurements in the multispectral domain, which were converted into the hyperspectral domain. The multispectral to hyperspectral profile conversion was performed using EO-1 Hyperion as a hyperspectral library.
The hyperspectral profile of L8 and S2A was achieved in two steps requiring two inputs, “Input 1” and “Input 2”, as shown in Figure 7. Input 1 consisted of the processed image ( ρ L 8 ), S2A ( ρ S 2 A ), and EO-1 Hyperion ( ρ h ). After going through the process of converting the multispectral domain to the hyperspectral domain as described in Figure 8, the output was the L8 and S2A hyperspectral profile. This hyperspectral profile was validated with L8 and S2A data to see if they represent the measured TOA reflectance. It was found that there was some discrepancy between the hyperspectral profile and the satellite measurements of approximately 1 unit TOA reflectance in the CA, green, and SWIR2 bands, as shown in Figure 9. Therefore, the EO-1 hyperspectral profile ( ρ h c a l ) was modified by applying the relative gain to the Hyperion data as described in Section 2.3.2. “Input 2” consisted of ρ L 8 , ρ S 2 A , and ρ h c a l , used to obtain the final hyperspectral profile for L8 and S2A. The detailed procedures to convert the multispectral to hyperspectral profile are explained in the sections below.

2.3.1. Estimation of Satellite Hyperspectral Profile

Figure 8 shows the detailed conversion process of the multispectral profile into its hyperspectral profile. There were three satellite data: L8 (M = 1150 images), S2A (P = 775 images), and EO-1 (N = 64 images), represented by ( ρ L 8 i = 1 M ), S2A ( ρ S 2 A k = 1 P ), and ( ρ h j = 1 N ), respectively. Using Equation (5), the band-integrated TOA reflectance ( ρ h s a t ) for the satellites viz. L8 ( ρ h L 8 j = 1 N ) and S2A ( ρ h S 2 A j = 1 N ) was estimated by integrating the relative spectral response (RSR) of the reference sensors ( R S R s a t ) with the Hyperion profile ( ρ h ) at each sampled wavelength (dλ), weighted by the corresponding RSR of the reference sensor.
ρ h s a t ( λ ) = λ 1 λ 2 ρ h × R S R s a t d λ λ 1 λ 2 R S R s a t d λ
The TOA reflectance ratio ( R L 8 i = 1 , j ) was estimated between L8 TOA reflectance ( ρ L 8 i = 1 ) for the first image with respect to the N number of EO-1 to L8-integrated TOA reflectance ( ρ h L 8 j = 1 N ) across all the bands (b1…b7). The average TOA reflectance ratio ( A L 8 i = 1 , j = 1 N ) was taken across the bands for the N number of R L 8 i = 1 , j = 1 N . EO-1 to L8-normalized TOA reflectance ρ n o r m L 8 i = 1 , j = 1 N was then computed by multiplying ρ h L 8 j = 1 N with A L 8 i = 1 , j = 1 N . Further, the mean-squared errors ( M S E i = 1 , j = 1 N ) were calculated among ρ n o r m L 8 j = 1 N with respect to ρ L 8 i = 1 using Equation (6) given below:
M S E i j = ( ρ s a t i ρ n o r m s a t i j ) 2 No . of Bands
where ρ s a t and ρ n o r m s a t are the TOA reflectance and normalized TOA reflectance of satellites (“sat” refers to either L8 or S2A). For ρ L 8 i = 1 , the minimum MSE value was observed from the N number of MSEs at instant j = a (say). The final hyperspectral spectrum ( ρ L 8 h y p e r i = 1 ) was then picked from the Hyperion library ( ρ h j = 1 N ) at j = a . Finally, the L8-normalized hyperspectral TOA reflectance ( ρ L 8 h y p i = 1 ) was obtained by further multiplying ρ L 8 h y p e r i = 1 with the average TOA reflectance ratio ( A L 8 i = 1 , j = a ) for the first ρ L 8 i = 1 . This process was repeated for M = 1150 L8-processed images to achieve the M number of L8 normalized hyperspectral TOA reflectance profiles. The same process was applied to the Sentinel-2A data to obtain P = 775 S2A-normalized hyperspectral profiles ( ρ S 2 A h y p k = 1 P ).
The normalized hyperspectral TOA reflectance for both L8 ( ρ L 8 h y p i ) and S2A ( ρ S 2 A h y p k ) was further integrated into the multispectral bands ( ρ L 8 h y p m u l t ) and compared with their observed TOA reflectance ( ρ L 8 i = 1 ) measurements to check the accuracy of their hyperspectral profile. Figure 9 shows the ρ L 8 i = 1 (orange color), ρ L 8 h y p i = 1 (grey color), ρ L 8 h y p m u l t (purple color), and absolute difference between ρ L 8 i = 1 and ρ L 8 h y p m u l t (red colored). It can be observed that there was a discrepancy of 0.0106 TOA reflectance between ρ L 8 i = 1 and ρ L 8 h y p m u l t . This might be because of EO-1 Hyperion’s different band-to-band attributes in their relative domain. This small offset can further be minimized by finding the relative gain and performing relative calibration on the Hyperion sensor, which brought the satellite (L8 and S2A) hyperspectral profile to match its multispectral domain.

2.3.2. Relative Calibration on EO-1 Sensor

In the relative calibration process, the ratios between satellite measurements ( ρ L 8 i and ρ S 2 A k ) and the band-integrated TOA reflectance ( ρ L 8 h y p m u l t and ρ S 2 A h y p m u l t ) were first calculated. Then, the mean ratio across the central wavelengths was taken for L8 and S2A, called the “relative gain” values.
Figure 10 shows the relative gain for L8 (orange color) and S2A (blue color) versus the central wavelength. The relative gain values for L8 and S2A were combined to be the “superspectral relative gain” values (orange color), as shown in Figure 11. In order to obtain the “hyperspectral relative gain” values, the modified Akima cubic Hermite (makima) interpolation was applied to the superspectral relative gains. The interpolated values can be obtained based on a piecewise function of polynomials with at most three degrees, which is the advantage of using the “makima” approach [38]. The interpolation was employed separately for the VNIR and SWIR channels due to the two sensors onboard Hyperion. Before employing the interpolation on the VNIR channel, the average relative gain between L8 and S2A was taken for the “green”, “red”, and “NIR” bands because they are very close to each other. Then, the interpolation was performed to obtain interpolated relative gain values with a 1 nm spectral resolution for the VNIR channels. Similarly, for the SWIR channels, the average relative gain of L8 and S2A was computed and taken as a relative gain with a 1 nm spectral resolution. The interpolated relative gain for VNIR and SWIR channels was further extrapolated and combined in order to cover the entire hyperspectral wavelength region. The EO-1 sensor used in this study has a hyperspectral TOA reflectance profile with a 10 nm spectral resolution. In order to match the EO-1 sensor’s resolution, the combined interpolated relative gain values were integrated to obtain a 10 nm resolution of the hyperspectral relative gain values.
Figure 11 shows the hyperspectral (grey color) and superspectral (orange-color) relative gain plot. Finally, the hyperspectral relative gain was multiplied with the TOA reflectance of the 64 Hyperion processed images ( ρ h j = 1 N ), which completed the relative calibration on the EO-1 sensor, as shown in Figure 12. The output of the relative calibration process, which is the relative calibrated EO-1 Hyperion processed image ( ρ h c a l ), was used as an input parameter for “Input 2” in Figure 7.

2.3.3. Normalized Hyperspectral Profile after Relative Calibration

The processed image information for L8 ( ρ L 8 ), S2A ( ρ S 2 A ) and relative calibrated EO-1 Hyperion-processed image information ( ρ h c a l ) formed a new input data called “Input 2”. Input 2, shown in Figure 7, was now used in the multispectral to hyperspectral profile conversion for L8 and S2A following the same procedures illustrated in Figure 8 and explained in Section 2.3.1. The normalized hyperspectral profile obtained from the conversion process using “Input 2” was then compared with their satellite measurements. Figure 13 shows the L8-normalized hyperspectral profile, ρ L 8 h y p (grey color), along with its band-integrated multispectral form, ρ L 8 h y p m u l t (purple circle). The orange circle represents the observed TOA reflectance ( ρ L 8 i = 1 ). The absolute difference (red bar) was calculated between ρ L 8 i = 1 and band-integrated multispectral TOA reflectance ( ρ L 8 h y p m u l t ), as shown in Figure 13. It can be concluded that performing the relative calibration on the EO-1 sensor reduced the offset and the hyperspectral profile highly matched its observed multispectral domain.
There are M and N numbers of hyperspectral profile for L8 ( ρ L 8 h y p ) and S2A ( ρ S 2 A h y p ), respectively. Figure 14 shows the L8- normalized hyperspectral profile (gray color), the L8 measurements (orange circle), and the relative calibrated Hyperion profile (blue color). The S2A-normalized hyperspectral profile (grey color), the S2A measurements (orange circle), and the relative calibrated Hyperion profile (blue color) are shown in Figure 15. It can be observed that ρ L 8 and ρ S 2 A aligned within the range of normalized hyperspectral profile for both the L8 and S2A sensors. It can also be justified that the normalized hyperspectral profile obtained in this process truly represents satellite data in the multispectral domain.
The DaHD consists of a combination of L8- and S2A-normalized hyperspectral profiles, a total of 1925 (M = 1150, N = 775) profiles. The DaHD, consisting of reflectance measurements in the range of 0.04–0.2, was then used to develop the DAHAC model.

2.4. Dark Hyperspectral Absolute Calibration Model Development

2.4.1. Four-Angle Hyperspectral BRDF Model

Since the Earth’s surface is the non-Lambertian target, solar illumination and sensor viewing geometry can significantly impact the TOA reflectance of a given target. Most often, the BRDF can be used to model this effect [1]. In this study, the DaHD was used to develop the four-angle hyperspectral BRDF model. During this process, the spherical angles associated with the DaHD viz. a solar zenith angle (SZA) ranging from 15–60°, a solar azimuth angle (SAA) from 31–163°, a view zenith angle (VZA) from 0.03–10°, and a view azimuth angle (VAA) from −177–180° were converted into the Cartesian coordinates form using Equations (7)–(10), respectively.
Y 1 = s i n ( S Z A ) × c o s ( S A A )
X 1 = s i n ( S Z A ) × s i n ( S A A )
Y 2 = s i n ( V Z A ) × c o s ( V A A )
X 2 = s i n ( V Z A ) × s i n ( V A A )
where X 1 and Y 1 represent the Cartesian coordinates obtained using SZA and SAA, while X 2 and Y 2 represent the Cartesian coordinates obtained using VZA and VAA.
In order to achieve symmetry with regard to the scattering plane and to provide a robust fit to the four-angle hyperspectral BRDF model, the DaHD along with the angular information in the converted Cartesian coordinates were then mirrored to each quadrant. Consequently, the TOA reflectance became a continuous function of independent variables [22]. Using the new set of the independent variables in the Cartesian coordinates form ( X 1 , Y 1 , X 2 , Y 2 ) and the mirrored DaHD, the four-angle multi-quadratic least-squares regression model was derived as shown in Equation (11). It is important to note that the four-angle hyperspectral BRDF model developed in this study utilizes all the DaHD imagery data to cover all the angular information.
ρ ( λ , X 1 , Y 1 , X 2 , Y 2 ) = β 0 + β 1 X 1 ( λ ) + β 2 Y 1 ( λ ) + β 3 X 2 ( λ ) + β 4 Y 2 ( λ ) + β 5 X 1 Y 1 ( λ ) + β 6 X 1 X 2 ( λ ) + β 7 X 1 Y 2 ( λ ) + β 8 Y 1 X 2 ( λ ) + β 9 Y 1 Y 2 ( λ ) + β 10 X 2 Y 2 ( λ ) + β 11 X 1 2 ( λ ) + β 12 Y 1 2 ( λ ) + β 13 X 2 2 ( λ ) + β 14 Y 2 2 ( λ )
where ρ is the predicted TOA reflectance with respect to the wavelength ( λ ) and four-angle Cartesian coordinates ( X 1 , Y 1 , X 2 , Y 2 ). β 0 represents the BRDF intercept, and the 14 BRDF model coefficients β 1 , β 2 , β 3 , β 4 , β 5 , β 6 , β 7 , β 8 , β 9 , β 10 , β 11 , β 12 , β 13 , and β 14 represent quadratic and linear coefficients.
Once the four-angle BRDF model had been developed, the predicted TOA reflectance was compared against the DaHD, as shown in Figure 16. The TOA reflectance difference was calculated between the DaHD profile and the predicted four-angle BRDF model to determine the level of agreement between them. The distribution of the TOA reflectance differences for all 196 bands is shown in the histogram in Figure 17. It can be observed that the TOA reflectance difference was normally distributed and the majority of the difference lied within a 0.02 TOA reflectance for all 196 hyperspectral bands.

2.4.2. Dark Hyperspectral Absolute Calibration Model

In order to develop the DAHAC model, a hypothesis test was applied to all 15 BRDF model coefficients at the 95% significance level for all 196 hyperspectral bands. The hypothesis test results (model coefficients, estimate, standard error, t-statistics values, p-values, and each coefficient statistical significance) for the band with a 864 nm wavelength are summarized in Table 3.
All of the Hyperion bands were validated by this statistical analysis, excluding a few bands at 942 nm and 1386 nm, which represent water vapor and cirrus with high absorption features. The p-value for the t-statistic hypothesis evaluated whether or not the BRDF coefficients were equal to zero. When the p-value exceeds 0.05, the corresponding BRDF coefficient becomes insignificant and vice versa. There were only seven significant coefficients viz. X 1 X 2 , Y 1 Y 2 , X 1 2 , Y 1 2 , X 2 2 , Y 2 2 , and the intercept, as shown in Table 3. They were used to represent the DAHAC model. It is important to note that the atmospheric parameters were not considered when developing the model in this study since they had a negligible effect on the atmospheric model and the absolute calibration model [39]. The DAHAC model with seven coefficients was derived using the solar zenith and azimuth angles and sensor zenith and azimuth angles, as shown in Equation (12) below:
ρ ( λ , X 1 , Y 1 , X 2 , Y 2 ) D A H A C = B 0 ( λ ) + B 1 ( λ ) X 1 X 2 + B 2 ( λ ) Y 1 Y 2 + B 3 ( λ ) X 1 2 + B 4 ( λ ) Y 1 2 + B 5 ( λ ) X 2 2 + B 6 ( λ ) Y 2 2
where ρ ( λ , X 1 , Y 1 , X 2 , Y 2 ) D A H A C is the model-predicted TOA reflectance with respect to the wavelength ( λ ) and the Cartesian coordinates ( X 1 , Y 1 , X 2 , Y 2 ). The predicted TOA reflectance had a spectral resolution of 10 nm. B 0 represents the intercept, and B 1 , B 2 , B 3 , B 4 , B 5 , and B 6 are the model coefficients corresponding to X 1 X 2 , Y 1 Y 2 , X 1 2 , Y 1 2 , X 2 2 , and Y 2 2 , respectively.
Unlike the hyperspectral model developed by Chaity et al. [1], it should be noted that the DAHAC model does not require any adjustment factor in the model equation to scale the Hyperion spectrum due to low reflectance in all bands. The DAHAC model coefficient values for 196 hyperspectral bands are shown in Figure 18.
The validation of the DAHAC model (7 coefficients) with the four-angle hyperspectral BRDF model (15 coefficients) for all hyperspectral bands is shown in Figure 19a. It can be observed that the plot shows a linear behavior, which means the DAHAC model statistically agrees with the four-angle hyperspectral BRDF model with a slope equal to one and a bias close to zero. Figure 19b shows the TOA reflectance differences between the two models, and the mean absolute error (MAE) calculated using Equation (13) was found to be close to zero.
M A E = j = 1 n | y j y ^ j | n
where y j : four-angle hyperspectral BRDF Model predicted TOA reflectance, y ^ j : DAHAC model-predicted TOA reflectance, and n: number of hyperspectral bands.
This provided significant evidence that the DAHAC model is capable of predicting the DaHD. The DAHAC model is further validated with Landsat missions and Sentinel-2 missions in Section 3. The model performance was evaluated using two metrics: accuracy and precision in unit reflectance due to the low-reflectance target. The accuracy and precision were measured by calculating the mean and standard deviation of the differences between the satellite-observed and model-predicted TOA reflectance using Equations (14) and (15), respectively.
Accuracy = m e a n ( ρ s a t ρ m o d e l )
Precision = s t d ( ρ s a t ρ m o d e l )
where ρ s a t refers to the satellite-observed TOA reflectance and ρ m o d e l represents the DAHAC-model-predicted TOA reflectance, i.e., ρ ( λ , X 1 , Y 1 , X 2 , Y 2 ) D A H A C .

3. Results and Discussion

In this section, the results of the model validation to determine the model accuracy with respect to different sensors are presented and analyzed.

3.1. DAHAC Model Validation Process

In order to evaluate the DAHAC model performance and its accuracy, it was validated with five different sensors viz. L7, L8, L9, S2A, and S2B. During the validation process, the spherical angular information (SZA, SAA, VZA, and VAA) was first transformed into the Cartesian coordinates form for each sensor (L7, L8, L9, S2A, and S2B). Using the converted Cartesian coordinates and the DAHAC model coefficients shown in Figure 18 as the input parameters for the DAHAC model, the hyperspectral TOA reflectance was predicted for each sensor. The hyperspectral TOA reflectance response from the model had a 10 nm resolution, which was then subjected to cubic interpolation to obtain a finer 1 nm spectral resolution. Since multispectral sensors were used for validation in this study, the RSR of each sensor was then used to integrate the interpolated hyperspectral TOA reflectance into the multispectral TOA reflectance. Finally, the model response in the form of multispectral TOA reflectance was compared with the observed TOA reflectance of each sensor used for validation purposes. The results from the DAHAC model validation for L7, L8, L9, S2A, and S2B are presented in the following subsections.

3.1.1. DAHAC Model Validation with Landsat Missions

The DAHAC model-predicted and the observed L8 TOA reflectance for seven bands were compared with respect to the decimal year and SZA, as shown in Figure 20a–g.
The left-hand side plots with respect to the decimal year show the seasonal variation and were well-captured by the model for every band. The right-hand side plots show the decreasing linear trend of the observed TOA reflectance with respect to the SZA, and the DAHAC model-predicted TOA reflectance agreed well with L8 measurements for all spectral bands.
The performance of the DAHAC model on the Landsat sensors (L7, L8, and L9) was evaluated by measuring the difference between the observed TOA reflectance of the Landsat sensors with the DAHAC model’s response. The error distribution for L7, L8, and L9 across the bands are further shown in the violin plot in Figure 21a–c, respectively. The error bar in each violin plot shows the mean and 1 σ standard deviation of the TOA reflectance difference. It is clearly seen from the violin distribution in Figure 21a–c that, for all three Landsat sensors (L7, L8, L9), the majority of TOA reflectance difference fell within ±0.02 for all spectral bands. The average mean difference ranged between 0 and 0.012 for all bands.

3.1.2. DAHAC Model Validation with Sentinel-2 Missions

The DAHAC model validation process for Sentinel-2 missions was the same as for the Landsat missions. Figure 22 shows the comparison between the predicted TOA reflectance and observed S2A measurements for non-Landsat bands (Red-Edge1, Red-Edge2, Red-Edge3, NIR1) with respect to the decimal year and SZA. It can be observed that the seasonal variation in the Sentinel measurements matched very well with the DAHAC-model-predicted TOA reflectance for all non-Landsat bands. Figure 23a,b show that the majority of the TOA reflectance difference fell within ±0.02 TOA reflectance for Sentinel-2 missions (S2A, S2B). The average mean difference ranged within 0 and 0.004 TOA reflectance for all spectral bands.
The hyperspectral BRDF model developed in [1] depends on the hyperspectral cross-scale factor to normalize the Hyperion intercept spectral profile for different spectral bands. The advantage of using the DAHAC model is that it will be independent of calculating and estimating the adjustment factors for scaling the Hyperion spectrum. The dark target data were obtained from all three sensors (L8, S2A, and EO-1) in the range of 0.04–0.2 TOA reflectance for all spectral bands used for the DAHAC model development. Following the DAHAC model specifications for the solar zenith angle (15–60°), solar azimuth angle (31–163°), view zenith angle (0.03–10°), and view azimuth angle (−177–180°), the predicted results showed that the model was able to estimate the TOA reflectance measurements with the highest accuracy of ±0.012 and a precision within 0.02 unit reflectance for all Landsat satellites spectral bands. The model performance for the Sentinel sensors (S2A and S2B) showed an accuracy and precision within ±0.004 and ±0.02 unit reflectance. Additionally, the TOA reflectance predicted using the DAHAC model had a 10 nm spectral resolution, and a finer 1 nm resolution was then obtained by further subjecting the model to cubic interpolation. This implies that the DAHAC model can perform the absolute calibration of any sensor over Dark EPICS-Global with a 1 nm spectral resolution.

4. Uncertainty Analysis

In this section, the uncertainty analysis for the DAHAC model is performed. The sensors viz. L8, S2A, and EO-1 Hyperion used for the DAHAC model’s development have currently been determined to have uncertainties at 2% [7], 4% [39], and 5–10% [35], respectively. It was initially assumed that the DAHAC model’s uncertainty is related to the uncertainty on the relative calibration of the EO-1 sensor, all the steps involved in the model development, the DaHD generation process, and the development of the BRDF model, as described in Section 2.3 and Section 2.4.
The flowchart shown in Figure 24 illustrates the steps followed in the DAHAC model’s uncertainty calculation.
The standard deviation for coefficient estimates of the DAHAC model was calculated using the formula: C o e f f s t d = C I U C I L 2 × t v a l u e , where C I U and C I L represent the upper bound and lower bound of the confidence interval (CI), respectively, and the t-value corresponds to the 95% CI and the degrees of freedom. The initial distribution of the mean and standard deviation of the seven DAHAC model coefficients for 196 bands are tabulated in Table A1 presented in Appendix A. Monte Carlo simulation was then employed to generate multivariate normal random numbers for the model coefficients using the model coefficients and their standard deviation information with the number of iterations set to 100. Using the randomized model coefficients obtained from this simulation and the angular information (SZA, SAA, VZA, VAA) from the DaHD in Cartesian coordinates, the DAHAC model predicted the TOA reflectance at each iteration level. The standard deviation of the predicted TOA reflectance of the DAHAC model was then estimated for all 196 bands. The simulation was performed by increasing the number of iterations at different levels (100, 500, 1000, 1500, 2000, 2500) until the standard deviation of the predicted TOA reflectance of the DAHAC model converged. Once the standard deviation converged to a stationary point, the simulation stopped, and the standard deviation at that point represented the DAHAC model’s uncertainty.
Figure 25 shows the results of the standard deviation of the DAHAC-model-predicted TOA reflectance using the Monte Carlo simulation. The three plots in Figure 25a–c represent different wavelength ranges: (426–905) nm, (912–1507) nm, and (1517–2395) nm, respectively. It can be observed that the standard deviation nearly converged between 1500 and 2500 iterations; the range was stable for all 196 bands. The uncertainty of the DAHAC model was then selected to be the model standard deviation at the 1500 iterations level. Figure 26 shows the DAHAC-model-predicted TOA reflectance and the DAHAC model’s uncertainty associated with each hyperspectral band. This analysis showed that the model uncertainty for the VNIR and SWIR channels was within 0.04 and 0.05 unit reflectance, respectively.

5. The DAHAC Model’s Double Ratio for Sensors Inter-Comparison

The ratio between the DAHAC model’s TOA reflectance and satellite measurements was used to determine the DAHAC model’s double ratio. The DAHAC model’s double ratio was used as a metric to perform the inter-comparison between two satellites across each spectral band. The double ratio approach eliminated any possible bias or error that could be embedded in the DAHAC model. Thus, it can be used as a transfer mechanism to directly compare two satellites, for instance L8 and L9 or S2A and S2B. During this process, the ratio between the model-predicted TOA reflectance and observed TOA reflectance measurements was calculated using Equation (16), called the “reflectance ratio”. Further, the DAHAC model’s double ratio was calculated using Equation (17).
Reflectance Ratio s a t = DAHAC   Model   TOA   Reflectance Observed   TOA   Reflectance
DAHAC   Model   Double   Ratio = Reflectance   Ratio s a t Reflectance   Ratio r e f
where “sat” refers to the L8, L9, S2A, and S2B satellite sensors. “ref” represents the reference sensors viz. L8 and S2A.
The image information from November 2021 to October 2022 was used to perform the inter-comparison between L8 and L9. For S2A and S2B, the data for the inter-comparison analysis were taken from July 2017 to July 2022 over the dark target. In order to prevent the study from being influenced by site variation resulting from various atmospheric conditions, the image information was assumed to be stable within seven days. According to Gross et al., a smaller difference in the view zenith angle between two sensors leads to a reduced impact of several aspects of the BRDF [40]. Using nearby coincident pairs within seven days with a view zenith angle (VZA) difference less than 2°, 46 and 617 pairs were obtained for L8 and L9 and S2A and S2B, respectively. The DAHAC model’s double ratio result showed that L8 and L9 were comparable within 1% for all bands except for the red band within 2%, as shown in Figure 27. The Sentinel-2A and 2-B inter-comparison result showed the agreement within 1% for eleven spectral bands, as shown in Figure 28. In this way, the DAHAC model’s double ratio helped to evaluate the measurement differences across each spectral band for different sensors.
Being a part of the Landsat Cal/Val team, the SDSU IP Lab performed the cross-calibration between L8 and L9, showing the agreement within 0.5–1% [41]. The cross-calibration result between S2A and S2B was within ±2% over the Libya 4 site [42]. The DAHAC model’s double ratio gave similar results as reported in previous studies [15,23,41]. It proved that the DAHAC model’s double ratio can be a complementary tool to be used for satellite calibration.

6. Conclusions

This paper presented the development of the novel DAHAC model using the normalized hyperspectral profile of L8 and S2A over Dark EPICS-Global. Dark EPICS-Global includes different surface types viz. dark rock, volcanic area, and dark sand with a temporal variation of 0.02 unit reflectance, as illustrated in Table 2. The measurements of the L8, S2A, and EO-1 Hyperion sensors were used to generate the dark hyperspectral data for the DAHAC model’s development. L8 and S2A are well-calibrated and were used as reference sensors and the EO-1 Hyperion sensor as the hyperspectral library source. This study proposed an algorithm to obtain the normalized hyperspectral profile for L8 and S2A using the multispectral (L8, S2A) and hyperspectral (EO-1 Hyperion) domain data. The relative calibration of the EO-1 sensor significantly improved the agreement between the hyperspectral profile of L8 and S2A with its multispectral domain. The DaHD was generated by aggregating the normalized hyperspectral profile of L8 and S2A.
The DaHD with TOA reflectance in the range of 0.04–0.2 across all spectral bands was used to develop the four-angle hyperspectral BRDF model. The four-angle hyperspectral BRDF model was further simplified to develop the DAHAC model by considering the significant coefficients in the model parameters. The simplified DAHAC model was leveraged to predict the TOA reflectance with a 10 nm spectral resolution and further subjected to cubic interpolation to achieve a finer 1 nm resolution. Thus, the DAHAC model can perform the absolute calibration of any sensor with a 1 nm spectral resolution over Dark EPICS-Global.
The DAHAC model was validated with several multispectral sensors viz. L7 ETM+, L8, L9, S2A, and S2B. The validation results showed that the model was able to predict the TOA reflectance measurements with the highest accuracy of ±0.012 and a precision within 0.02 unit reflectance for all spectral bands of the Landsat missions (L7, L8, and L9). The model performance for the Sentinel-2 missions (S2A and S2B) showed an accuracy and precision within ±0.004 and ±0.02 unit reflectance. In contrast to the hyperspectral APICS model estimation accuracy within limited bands [1], the DAHAC model can predict the satellite measurements across all the spectral bands, especially for non-Landsat equivalent bands with higher accuracy. This result proved the DAHAC model is a robust tool for performing absolute calibration of any sensors on multispectral and hyperspectral instruments.
A thorough analysis was performed to determine the model’s overall uncertainty. Using the Monte Carlo simulation on the DAHAC model, the DAHAC-model-predicted TOA reflectance and the standard deviation were estimated at different iteration levels. The DAHAC model’s uncertainty was found to be within 0.04 and 0.05 unit reflectance for the VNIR and SWIR channels, respectively.
The sensor inter-comparison for the Landsat and Sentinel-2 missions was performed using the DAHAC model’s double ratio. The agreement between L8 and L9 was within 0.06–2% and 1% for S2A and S2B. These results match closely the results from the L8–L9 cross-calibration (0.5–1%) [41] and the S2A-S2B inter-comparison over Libya 4 (2%) [42].
The DAHAC model is capable of a nadir look angle of up to 10° and a solar zenith angle of up to 60°, however, it would be more challenging for larger viewing angles. When the view zenith angle increases, the RSR shifts towards a shorter wavelength, leading to changes in the sensor’s received radiance information [43]. However, the RSR shift effect was not considered in this study. In addition, the DAHAC model does not consider factors such as atmospheric scatterings, Rayleigh scattering, aerosol optical load, and gas absorption properties before developing the model.
However, increasing the DAHAC model’s accuracy, decreasing its uncertainty, and improving the TOA reflectance consistency are still possible. The calibration variations between L8 and S2A with other sensors and random anomalies at the time of the sensor overpass due to atmospheric conditions were most likely the cause of the reduced accuracy. In general, the model accuracy can be improved by using more data samples with high-quality hyperspectral images for absolute calibration and sensor inter-comparison applications.
In the future, the proposed model’s development algorithm will be applied over stable bright sites to check its reliability and performance accuracy. The ability of the proposed algorithm to perform satellite absolute calibration will be tested for the dynamic range of TOA reflectance measurement. The RSR shift effect will also be considered during the absolute calibration model’s development with larger viewing angles of satellite data.

Author Contributions

Conceptualization, P.B.K., M.K., L.L. and C.T.P.; methodology, P.B.K., M.K. and L.L.; software, P.B.K.; validation, P.B.K.; formal analysis, P.B.K., L.L., M.K. and C.T.P.; investigation, P.B.K.; resources, M.K. and L.L.; writing—original draft preparation, P.B.K.; writing—review and editing, M.K. and L.L.; visualization, P.B.K.; supervision, M.K. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by USGS EROS, Grant Number SA2000371.

Data Availability Statement

Not applicable.

Acknowledgments

Thanks to Juliana Fajardo Rueda, Dinithi Siriwardana Pathiranage, and all the members of the Image Processing Laboratory for their continuous support and assistance. Finally, the authors would like to gratefully acknowledge all the Reviewers for their insightful comments and suggestions for improving the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The initial distribution of mean and standard deviation of the seven coefficients of the DAHAC model are given below.
Table A1. Showing the initial distribution of the mean and standard deviation of the DAHAC model’s coefficients for 196 bands.
Table A1. Showing the initial distribution of the mean and standard deviation of the DAHAC model’s coefficients for 196 bands.
Wavelength (nm) B 0 mean B 0 SD B 1 mean B 1 SD B 2 mean B 2 SD B 3 mean B 3 SD B 4 mean B 4 SD B 5 mean B 5 SD B 6 mean B 6 SD
426.80.1685.90E-040.3081.70E-020.26.90E-03−0.0378.00E-04−0.0293.20E-0317.2434.70E-01−1.6664.30E-02
4370.1555.50E-040.3141.60E-020.1756.40E-03−0.0347.40E-04−0.0233.00E-0316.4534.40E-01−1.5424.00E-02
447.20.1344.90E-040.2721.40E-020.1525.70E-03−0.0296.60E-04−0.0212.70E-0315.83.90E-01−1.4633.60E-02
457.30.134.70E-040.2741.40E-020.1455.50E-03−0.0296.40E-04−0.0222.60E-0315.3263.80E-01−1.3913.50E-02
467.50.1314.80E-040.2811.40E-020.1495.60E-03−0.036.50E-04−0.0252.60E-0316.1893.90E-01−1.4553.50E-02
477.70.1334.90E-040.2821.40E-020.1555.70E-03−0.0336.60E-04−0.0282.70E-0316.3353.90E-01−1.4493.60E-02
487.90.1274.60E-040.2751.30E-020.1495.40E-03−0.0346.20E-04−0.0292.50E-0315.263.70E-01−1.3353.40E-02
4980.1254.60E-040.2721.30E-020.1485.30E-03−0.0356.20E-04−0.032.50E-0314.9983.60E-01−1.2993.40E-02
508.20.1264.50E-040.271.30E-020.1545.30E-03−0.0376.10E-04−0.0342.50E-0314.5613.60E-01−1.2473.40E-02
518.40.1224.40E-040.261.30E-020.155.10E-03−0.0385.90E-04−0.0352.40E-0313.243.50E-01−1.123.20E-02
528.60.1214.20E-040.2521.20E-020.154.90E-03−0.0415.70E-04−0.0372.30E-0311.4793.40E-01−0.9533.10E-02
538.70.124.10E-040.2411.20E-020.1514.80E-03−0.0435.50E-04−0.042.20E-039.9413.30E-01−0.8123.00E-02
548.90.1194.00E-040.2321.20E-020.1534.70E-03−0.0465.50E-04−0.0432.20E-038.6453.20E-01−0.693.00E-02
559.10.1214.10E-040.2251.20E-020.1584.80E-03−0.0495.50E-04−0.0462.20E-037.0543.30E-01−0.5463.00E-02
569.30.124.00E-040.2061.20E-020.164.70E-03−0.0525.40E-04−0.0512.20E-035.0033.20E-01−0.3672.90E-02
579.50.124.00E-040.1891.20E-020.1644.70E-03−0.0565.40E-04−0.0552.20E-032.6533.20E-01−0.1613.00E-02
589.60.1264.30E-040.1911.30E-020.1795.10E-03−0.0625.90E-04−0.0632.40E-032.0913.50E-01−0.1033.20E-02
599.80.1284.60E-040.1811.30E-020.1935.40E-03−0.0666.20E-04−0.0722.50E-031.0913.70E-01−0.0123.40E-02
6100.134.60E-040.1821.30E-020.1835.30E-03−0.0676.20E-04−0.0672.50E-03−0.6643.70E-010.1443.40E-02
620.10.1284.50E-040.1751.30E-020.1815.30E-03−0.0676.10E-04−0.0672.50E-03−1.0253.60E-010.1753.30E-02
630.30.1274.60E-040.1681.30E-020.1825.40E-03−0.0686.20E-04−0.0692.50E-03−2.2413.70E-010.2843.40E-02
640.50.1294.70E-040.1681.40E-020.1865.50E-03−0.076.40E-04−0.0712.60E-03−2.9113.80E-010.3463.50E-02
650.70.1355.10E-040.1691.50E-020.2076.00E-03−0.0766.90E-04−0.0812.80E-03−2.8894.10E-010.353.80E-02
660.90.1315.00E-040.1621.50E-020.25.90E-03−0.0756.80E-04−0.0782.70E-03−3.3824.00E-010.3913.70E-02
6710.1365.10E-040.1691.50E-020.1976.00E-03−0.0776.90E-04−0.0772.80E-03−4.424.10E-010.4853.80E-02
681.20.1395.20E-040.1711.50E-020.2016.10E-03−0.0797.10E-04−0.0792.90E-03−4.9674.20E-010.5373.90E-02
691.40.1325.20E-040.1551.50E-020.1976.10E-03−0.0787.00E-04−0.0792.90E-03−6.1924.20E-010.653.80E-02
701.60.1385.90E-040.1521.70E-020.2296.90E-03−0.0868.00E-04−0.0943.20E-03−6.8584.70E-010.724.40E-02
711.70.1355.90E-040.1481.70E-020.2257.00E-03−0.0868.00E-04−0.0953.30E-03−7.724.80E-010.8014.40E-02
721.90.1367.40E-040.1322.20E-020.2768.70E-03−0.0971.00E-03−0.1224.10E-03−9.3435.90E-010.9765.50E-02
732.10.1386.70E-040.1392.00E-020.2487.80E-03−0.0939.00E-04−0.1073.70E-03−9.9265.30E-011.014.90E-02
742.20.1385.70E-040.161.70E-020.1936.60E-03−0.0857.70E-04−0.0783.10E-03−10.7284.50E-011.0574.20E-02
752.40.1325.50E-040.1561.60E-020.1716.40E-03−0.0827.40E-04−0.073.00E-03−12.1144.40E-011.184.00E-02
762.60.14.70E-040.1031.40E-020.1395.50E-03−0.0686.30E-04−0.0592.60E-03−12.5213.70E-011.2063.50E-02
772.80.1375.80E-040.1561.70E-020.1756.80E-03−0.0877.90E-04−0.0723.20E-03−14.4174.70E-011.3844.30E-02
7830.1415.90E-040.1621.70E-020.1826.90E-03−0.0888.00E-04−0.0753.20E-03−13.8134.70E-011.3384.40E-02
793.10.1396.10E-040.1531.80E-020.1967.10E-03−0.098.20E-04−0.0823.30E-03−13.8094.90E-011.3424.50E-02
803.30.1376.10E-040.1471.80E-020.1977.10E-03−0.098.30E-04−0.0843.30E-03−14.0174.90E-011.3614.50E-02
813.50.136.70E-040.1212.00E-020.2317.90E-03−0.0939.10E-04−0.1033.70E-03−13.6615.40E-011.3545.00E-02
823.60.1266.90E-040.1112.00E-020.2358.10E-03−0.0939.30E-04−0.1073.80E-03−13.8575.50E-011.3725.10E-02
833.80.1296.30E-040.1241.90E-020.2117.40E-03−0.098.50E-04−0.0933.50E-03−14.1595.10E-011.3844.70E-02
8440.1366.00E-040.1511.80E-020.1797.10E-03−0.0898.10E-04−0.0763.30E-03−15.5244.80E-011.4944.40E-02
854.20.1396.10E-040.1611.80E-020.1697.10E-03−0.0898.20E-04−0.073.30E-03−16.4034.90E-011.5744.50E-02
864.40.1366.00E-040.161.80E-020.1577.00E-03−0.0878.10E-04−0.0653.30E-03−16.9834.80E-011.6244.40E-02
874.50.1326.00E-040.1571.70E-020.1497.00E-03−0.0868.10E-04−0.0623.30E-03−17.694.80E-011.6884.40E-02
884.70.1276.00E-040.1521.70E-020.1457.00E-03−0.0858.10E-04−0.0613.30E-03−18.7314.80E-011.7844.40E-02
894.90.1196.30E-040.1191.80E-020.1817.40E-03−0.0878.50E-04−0.0823.50E-03−17.8355.10E-011.7124.70E-02
9050.1086.90E-040.0882.00E-020.2058.10E-03−0.0889.40E-04−0.0963.80E-03−17.9355.50E-011.7375.10E-02
912.50.1641.30E-030.0093.80E-020.4871.50E-02−0.1331.70E-03−0.2047.00E-03−7.5231.00E+000.6569.50E-02
922.50.1461.20E-03−0.0133.40E-020.4281.40E-02−0.1191.60E-03−0.1886.30E-03−10.0739.30E-010.9388.50E-02
932.60.1021.10E-030.0173.30E-020.3491.30E-02−0.1011.50E-03−0.1726.20E-03−13.3699.00E-011.2688.30E-02
942.70.1051.20E-030.0343.40E-020.3631.40E-02−0.1071.60E-03−0.1796.50E-03−13.7969.40E-011.3238.70E-02
952.80.1081.10E-030.0553.30E-020.3521.30E-02−0.1061.50E-03−0.1756.20E-03−12.3929.10E-011.248.40E-02
962.90.1249.90E-040.0982.90E-020.3361.20E-02−0.1081.30E-03−0.1575.40E-03−11.5457.90E-011.1547.30E-02
9730.1318.10E-040.1352.40E-020.2749.50E-03−0.1031.10E-03−0.1264.50E-03−13.7096.50E-011.3736.00E-02
983.10.1376.80E-040.1722.00E-020.2188.00E-03−0.0979.20E-04−0.0953.70E-03−13.9585.50E-011.3655.00E-02
993.20.1416.10E-040.1571.80E-020.1967.10E-03−0.0918.20E-04−0.0793.30E-03−13.7994.90E-011.3044.50E-02
1003.30.1395.90E-040.1541.70E-020.1867.00E-03−0.0898.00E-04−0.0743.30E-03−13.9224.80E-011.314.40E-02
1013.30.145.70E-040.161.70E-020.1776.70E-03−0.0877.70E-04−0.0653.10E-03−13.3664.60E-011.2464.20E-02
1023.40.145.70E-040.1541.70E-020.1756.70E-03−0.0867.70E-04−0.0643.10E-03−13.4164.60E-011.2524.20E-02
1033.40.1325.40E-040.1581.60E-020.1546.30E-03−0.0817.30E-04−0.0573.00E-03−12.9814.30E-011.2284.00E-02
1043.50.1315.30E-040.1551.60E-020.1526.20E-03−0.0797.20E-04−0.0562.90E-03−12.6614.30E-011.1973.90E-02
1053.60.1265.50E-040.1431.60E-020.1456.40E-03−0.0797.40E-04−0.0613.00E-03−14.3964.40E-011.3764.00E-02
1063.70.1235.40E-040.1371.60E-020.1476.40E-03−0.0797.30E-04−0.0633.00E-03−14.3614.30E-011.3734.00E-02
1073.80.1265.70E-040.1361.70E-020.1696.70E-03−0.0857.80E-04−0.0743.10E-03−15.2614.60E-011.4574.20E-02
1083.90.1255.90E-040.1221.70E-020.1886.90E-03−0.0878.00E-04−0.0833.20E-03−14.7194.70E-011.4154.40E-02
10940.126.70E-040.0942.00E-020.2217.90E-03−0.0929.10E-04−0.1023.70E-03−15.1055.40E-011.495.00E-02
1104.10.1088.20E-040.0472.40E-020.2679.60E-03−0.0941.10E-03−0.134.50E-03−14.5966.60E-011.4726.00E-02
1114.10.0738.20E-040.0072.40E-020.2379.70E-03−0.0771.10E-03−0.1284.50E-03−12.4116.60E-011.2866.10E-02
1124.20.0588.40E-040.0232.50E-020.229.80E-03−0.0731.10E-03−0.1254.60E-03−12.5666.70E-011.296.20E-02
1134.30.0699.00E-040.0372.60E-020.2471.10E-02−0.0821.20E-03−0.1344.90E-03−13.3757.20E-011.3656.60E-02
11440.0739.40E-040.0322.70E-020.2611.10E-02−0.0861.30E-03−0.145.10E-03−13.7057.50E-011.4026.90E-02
11540.0929.70E-040.0552.90E-020.2911.10E-02−0.0991.30E-03−0.1475.30E-03−14.3737.80E-011.4527.20E-02
11640.118.30E-040.0772.40E-020.2729.80E-03−0.0991.10E-03−0.134.60E-03−14.2136.70E-011.4246.10E-02
11740.1147.50E-040.0862.20E-020.248.70E-03−0.0961.00E-03−0.1124.10E-03−15.7846.00E-011.5525.50E-02
11840.1147.30E-040.0942.10E-020.2318.50E-03−0.0959.80E-04−0.1084.00E-03−15.9225.80E-011.5665.40E-02
11940.1116.60E-040.1161.90E-020.1977.80E-03−0.099.00E-04−0.0953.60E-03−15.6915.30E-011.5674.90E-02
12050.1116.40E-040.1271.90E-020.1877.50E-03−0.0898.70E-04−0.093.50E-03−15.8775.20E-011.5784.80E-02
12150.1166.10E-040.1471.80E-020.1697.20E-03−0.0898.30E-04−0.0783.40E-03−15.8274.90E-011.5514.50E-02
12250.1175.90E-040.1611.70E-020.1456.90E-03−0.0868.00E-04−0.0653.20E-03−16.1954.70E-011.5774.40E-02
12350.1175.70E-040.1481.70E-020.1286.70E-03−0.0817.70E-04−0.0583.10E-03−16.3864.60E-011.594.20E-02
12450.1155.70E-040.1441.70E-020.126.70E-03−0.087.70E-04−0.0543.10E-03−16.9734.50E-011.6434.20E-02
12550.1145.40E-040.1211.60E-020.1386.40E-03−0.087.30E-04−0.0553.00E-03−16.824.30E-011.5964.00E-02
12650.1075.40E-040.111.60E-020.1416.30E-03−0.0797.30E-04−0.0573.00E-03−16.8774.30E-011.594.00E-02
12750.1145.60E-040.1151.60E-020.1436.50E-03−0.0837.60E-04−0.0583.10E-03−17.6274.50E-011.6584.10E-02
12850.125.70E-040.1131.70E-020.1626.70E-03−0.0867.70E-04−0.0673.10E-03−16.7774.60E-011.5934.20E-02
12950.1135.60E-040.0961.60E-020.1796.60E-03−0.0837.60E-04−0.0773.10E-03−14.1274.50E-011.3714.10E-02
13050.1066.30E-040.0651.90E-020.2157.40E-03−0.0858.60E-04−0.0973.50E-03−13.1015.10E-011.2964.70E-02
13160.17.00E-040.0482.10E-020.2448.20E-03−0.0869.50E-04−0.1133.90E-03−11.745.60E-011.1745.20E-02
13260.0898.00E-040.0172.30E-020.2599.30E-03−0.0841.10E-03−0.1274.40E-03−11.8576.40E-011.25.90E-02
13360.0798.70E-040.022.60E-020.2451.00E-02−0.0841.20E-03−0.1314.80E-03−14.887.00E-011.4946.40E-02
13460.0395.30E-04−0.0041.50E-020.1196.20E-03−0.0467.10E-04−0.0732.90E-03−10.894.20E-011.0883.90E-02
13560.0049.30E-050.0212.70E-03−0.0051.10E-03−0.0051.30E-04−0.0025.10E-04−2.0677.40E-020.2276.90E-03
13660.0023.80E-050.0041.10E-0304.50E-04−0.0025.20E-05−0.0032.10E-04−0.6033.10E-020.082.80E-03
13760.0031.00E-040.0192.90E-03−0.0011.20E-03−0.0051.30E-04−0.0025.50E-04−2.0588.00E-020.2157.30E-03
13860.0051.20E-040.0193.70E-030.0031.50E-03−0.0081.70E-04−0.0066.80E-04−3.1021.00E-010.3199.20E-03
13960.0081.70E-040.0215.10E-030.0162.00E-03−0.0142.30E-04−0.0139.50E-04−5.0061.40E-010.4961.30E-02
14060.0132.60E-040.0187.50E-030.0353.00E-03−0.0223.50E-04−0.0261.40E-03−7.3972.10E-010.7241.90E-02
14160.0244.30E-040.0161.30E-020.0775.00E-03−0.0375.80E-04−0.0522.40E-03−10.5493.40E-011.0383.20E-02
14260.0315.50E-040.021.60E-020.116.40E-03−0.0487.40E-04−0.073.00E-03−12.4124.40E-011.2324.00E-02
14370.0366.30E-040.0171.80E-020.1277.40E-03−0.0558.50E-04−0.083.50E-03−14.7495.00E-011.4454.60E-02
14470.0487.80E-040.0372.30E-020.1769.10E-03−0.071.10E-03−0.1044.30E-03−15.9686.20E-011.5915.80E-02
14570.0669.70E-040.0232.80E-020.2461.10E-02−0.091.30E-03−0.1335.30E-03−18.2627.80E-011.8057.20E-02
14670.0578.90E-040.0122.60E-020.2151.00E-02−0.0811.20E-03−0.1194.90E-03−17.9377.10E-011.7586.60E-02
14770.079.90E-040.0032.90E-020.2621.20E-02−0.0911.30E-03−0.1375.40E-03−18.4777.90E-011.7537.30E-02
14870.0871.00E-030.0293.00E-020.2961.20E-02−0.1021.40E-03−0.1475.60E-03−17.6388.20E-011.6757.50E-02
14970.1028.20E-040.0632.40E-020.2799.60E-03−0.0941.10E-03−0.134.50E-03−12.2376.60E-011.1756.10E-02
15070.1117.20E-040.092.10E-020.2588.50E-03−0.0939.80E-04−0.1154.00E-03−11.7275.80E-011.1175.30E-02
15170.1156.70E-040.0952.00E-020.2327.80E-03−0.0929.00E-04−0.13.70E-03−14.365.30E-011.3254.90E-02
15270.1186.10E-040.1091.80E-020.2057.20E-03−0.0898.30E-04−0.0853.40E-03−15.0824.90E-011.3734.50E-02
15370.125.70E-040.1331.70E-020.1866.70E-03−0.0857.70E-04−0.0763.10E-03−13.8714.50E-011.264.20E-02
15480.1225.70E-040.141.70E-020.1786.60E-03−0.0847.70E-04−0.0713.10E-03−14.1714.50E-011.2814.20E-02
15580.1225.70E-040.1441.70E-020.1716.70E-03−0.0847.70E-04−0.0683.10E-03−14.8594.60E-011.3464.20E-02
15680.1155.50E-040.1291.60E-020.1636.50E-03−0.0817.50E-04−0.0653.00E-03−15.4554.40E-011.3974.10E-02
15780.1155.50E-040.1531.60E-020.1596.40E-03−0.0817.40E-04−0.0643.00E-03−13.6454.40E-011.254.00E-02
15880.1195.50E-040.1551.60E-020.1636.40E-03−0.0817.40E-04−0.0643.00E-03−12.9814.40E-011.1934.00E-02
15980.1145.40E-040.1331.60E-020.1556.40E-03−0.0797.30E-04−0.0613.00E-03−15.1014.30E-011.3594.00E-02
16080.1155.50E-040.1451.60E-020.1586.50E-03−0.0817.50E-04−0.0613.00E-03−14.8184.40E-011.334.10E-02
16180.125.50E-040.1591.60E-020.1576.40E-03−0.0817.40E-04−0.063.00E-03−12.7244.40E-011.1814.00E-02
16280.1195.40E-040.1531.60E-020.1566.30E-03−0.087.30E-04−0.063.00E-03−12.6814.30E-011.1794.00E-02
16380.1185.40E-040.151.60E-020.1596.30E-03−0.087.20E-04−0.0612.90E-03−12.44.30E-011.1454.00E-02
16480.1195.40E-040.1541.60E-020.1656.40E-03−0.0817.40E-04−0.0643.00E-03−12.2434.40E-011.1294.00E-02
16590.1185.30E-040.151.60E-020.1666.20E-03−0.087.20E-04−0.0642.90E-03−11.2194.30E-011.053.90E-02
16690.1195.40E-040.1521.60E-020.1686.30E-03−0.0817.30E-04−0.0652.90E-03−10.9664.30E-011.0254.00E-02
16790.1195.30E-040.1461.60E-020.1676.20E-03−0.0817.20E-04−0.0642.90E-03−11.1244.30E-011.0483.90E-02
16890.1185.40E-040.1341.60E-020.1746.30E-03−0.0827.20E-04−0.0692.90E-03−11.0914.30E-011.0483.90E-02
16990.1165.50E-040.131.60E-020.1736.40E-03−0.0837.40E-04−0.0693.00E-03−12.2444.40E-011.1654.00E-02
17090.1165.60E-040.1281.60E-020.1816.60E-03−0.0857.60E-04−0.0743.10E-03−12.1754.50E-011.1624.10E-02
17190.1135.60E-040.1191.60E-020.1896.50E-03−0.0837.60E-04−0.083.10E-03−11.0974.50E-011.0664.10E-02
17290.1115.70E-040.1111.70E-020.1976.70E-03−0.0837.70E-04−0.0843.10E-03−10.8924.60E-011.0514.20E-02
17390.1116.10E-040.11.80E-020.2187.20E-03−0.0878.30E-04−0.0943.40E-03−10.5974.90E-011.0274.50E-02
17490.1086.20E-040.091.80E-020.2247.30E-03−0.0868.40E-04−0.0973.40E-03−10.3455.00E-011.0034.60E-02
17590.1036.50E-040.0681.90E-020.2387.70E-03−0.0858.80E-04−0.1053.60E-03−9.3245.20E-010.924.80E-02
17690.0967.60E-040.042.20E-020.2659.00E-03−0.0871.00E-03−0.1224.20E-03−9.9486.10E-010.9875.60E-02
17800.0868.60E-040.0182.50E-020.2711.00E-02−0.0881.20E-03−0.1314.70E-03−11.4746.90E-011.1356.30E-02
17900.0728.70E-04−0.0092.50E-020.2541.00E-02−0.0821.20E-03−0.1284.80E-03−12.5357.00E-011.2316.40E-02
18000.0456.60E-04−0.0171.90E-020.1727.70E-03−0.0588.90E-04−0.0923.60E-03−11.2095.30E-011.0744.90E-02
18100.0294.70E-04−0.0131.40E-020.1095.50E-03−0.046.40E-04−0.0622.60E-03−9.5763.80E-010.9013.50E-02
18200.0192.40E-04−0.0056.90E-030.0612.80E-03−0.0213.20E-04−0.0331.30E-03−4.5671.90E-010.4021.70E-02
18300.0075.40E-050.0091.60E-030.0116.30E-04−0.0027.30E-05−0.0033.00E-04−0.2924.30E-02−0.0084.00E-03
18400.0027.10E-050.0242.10E-03−0.0068.40E-040.0019.60E-050.0053.90E-041.2615.70E-02−0.0995.30E-03
18500.0026.80E-050.0232.00E-03−0.0057.90E-040.0019.20E-050.0043.70E-041.2295.40E-02−0.0985.00E-03
18600.0071.10E-040.0023.10E-030.0151.20E-03−0.0081.40E-04−0.0095.80E-04−2.438.50E-020.2297.80E-03
18700.0057.20E-050.0042.10E-030.0088.40E-04−0.0059.70E-05−0.0053.90E-04−1.4765.70E-020.1385.30E-03
18800.011.80E-040.0025.30E-030.032.10E-03−0.0152.50E-04−0.021.00E-03−5.0621.50E-010.4671.30E-02
18910.0091.60E-040.0014.80E-030.0281.90E-03−0.0142.20E-04−0.0179.00E-04−4.5361.30E-010.4171.20E-02
19010.0068.90E-050.022.60E-030.0061.00E-03−0.0071.20E-04−0.0054.90E-04−1.5397.20E-020.1646.60E-03
19110.0091.30E-040.0163.80E-030.0181.50E-03−0.0111.80E-04−0.0137.10E-04−2.6561.00E-010.2729.60E-03
19210.012.00E-040.0065.90E-030.0312.40E-03−0.0162.70E-04−0.0221.10E-03−5.4671.60E-010.5371.50E-02
19310.0152.90E-040.0028.50E-030.0493.40E-03−0.0233.90E-04−0.0351.60E-03−7.6812.30E-010.7482.10E-02
19410.0315.60E-04−0.0161.60E-020.1216.60E-03−0.0467.60E-04−0.073.10E-03−13.1694.50E-011.24.10E-02
19510.0315.30E-04−0.0171.60E-020.1216.20E-03−0.0457.20E-04−0.0672.90E-03−12.2954.20E-011.1163.90E-02
19610.0376.30E-040.031.80E-020.1217.30E-03−0.0558.50E-04−0.0733.40E-03−14.9325.00E-011.4694.60E-02
19710.0648.00E-040.0692.30E-020.1999.40E-03−0.081.10E-03−0.1054.40E-03−15.3146.40E-011.5475.90E-02
19810.0898.60E-040.062.50E-020.2671.00E-02−0.0951.20E-03−0.1274.70E-03−15.0616.90E-011.4266.30E-02
19910.0727.10E-040.0122.10E-020.2078.30E-03−0.0779.60E-04−0.13.90E-03−16.0135.70E-011.4985.20E-02
20020.0233.90E-040.0071.10E-020.0614.50E-03−0.0355.20E-04−0.0342.10E-03−13.1083.10E-011.2192.90E-02
20120.034.90E-040.0191.40E-020.0775.70E-03−0.0466.60E-04−0.0432.70E-03−16.3643.90E-011.5183.60E-02
20220.0747.00E-040.0642.10E-020.1928.20E-03−0.089.50E-04−0.0863.80E-03−17.8965.60E-011.5855.20E-02
20320.1026.80E-040.0952.00E-020.2217.90E-03−0.099.20E-04−0.0933.70E-03−16.1735.40E-011.4275.00E-02
20420.096.40E-040.0741.90E-020.1797.50E-03−0.0848.70E-04−0.083.50E-03−18.7535.10E-011.7344.70E-02
20520.0686.30E-040.0441.90E-020.1427.40E-03−0.0748.60E-04−0.0663.50E-03-21.4925.10E-011.9744.70E-02
20620.0795.90E-040.0511.70E-020.1776.90E-03−0.0738.00E-04−0.0763.20E-03−16.3554.70E-011.4424.30E-02
20720.0946.30E-040.091.80E-020.2017.40E-03−0.0838.50E-04−0.0833.50E-03−15.6825.00E-011.3684.60E-02
20820.1016.30E-040.1081.80E-020.1817.40E-03−0.0888.50E-04−0.0753.40E-03−17.1255.00E-011.5874.60E-02
20920.1076.40E-040.1181.90E-020.1917.50E-03−0.0918.60E-04−0.0793.50E-03−16.115.10E-011.5084.70E-02
21020.1115.60E-040.1331.60E-020.1976.50E-03−0.0847.50E-04−0.0773.00E-03−10.0814.50E-010.9274.10E-02
21130.1125.60E-040.1311.60E-020.1996.60E-03−0.0857.60E-04−0.0783.10E-03−10.2254.50E-010.9354.10E-02
21230.1096.30E-040.1511.80E-020.1767.40E-03−0.0918.50E-04−0.073.40E-03−15.055.00E-011.4164.60E-02
21330.116.20E-040.1541.80E-020.1647.30E-03−0.098.40E-04−0.0623.40E-03−15.3325.00E-011.4394.60E-02
21430.1165.60E-040.1321.60E-020.2016.50E-03−0.0847.50E-04−0.0743.00E-03−10.7124.40E-010.9294.10E-02
21530.1135.80E-040.1231.70E-020.2146.70E-03−0.0867.80E-04−0.0823.20E-03−10.6914.60E-010.9244.20E-02
21630.1056.40E-040.1351.90E-020.1767.50E-03−0.098.60E-04−0.0683.50E-03−16.1485.10E-011.4984.70E-02
21730.1056.40E-040.1321.90E-020.1777.50E-03−0.098.70E-04−0.0693.50E-03−16.2585.10E-011.5084.70E-02
21830.1095.30E-040.1351.50E-020.2096.20E-03−0.0787.10E-04−0.0822.90E-03−5.2354.20E-010.4533.90E-02
21930.1055.00E-040.1261.50E-020.1965.90E-03−0.0746.80E-04−0.0772.70E-03−5.7024.00E-010.4933.70E-02
22030.1015.50E-040.1261.60E-020.1826.40E-03−0.087.40E-04−0.0673.00E-03−11.4994.40E-011.0144.00E-02
22130.1065.50E-040.1261.60E-020.1866.50E-03−0.0817.50E-04−0.0673.00E-03−11.3374.40E-010.9944.10E-02
22240.1064.70E-040.1341.40E-020.1855.50E-03−0.0716.30E-04−0.0672.60E-03−4.7073.80E-010.3893.50E-02
22340.1034.60E-040.1291.30E-020.1795.40E-03−0.076.20E-04−0.0652.50E-03−4.93.70E-010.4073.40E-02
22440.0975.20E-040.1121.50E-020.1516.10E-03−0.0757.00E-04−0.0572.90E-03−11.6244.20E-011.13.80E-02
22540.0945.20E-040.0991.50E-020.156.10E-03−0.0757.10E-04−0.0572.90E-03−12.6224.20E-011.1873.90E-02
22640.0994.80E-040.0881.40E-020.1895.60E-03−0.0726.50E-04−0.0692.60E-03−7.9693.80E-010.6743.50E-02
22740.0994.70E-040.0881.40E-020.1885.50E-03−0.0716.40E-04−0.0682.60E-03−7.6253.80E-010.6363.50E-02
22840.0924.40E-040.1021.30E-020.1545.10E-03−0.0665.90E-04−0.0612.40E-03−6.1843.50E-010.6193.20E-02
22940.0894.60E-040.0941.30E-020.1595.40E-03−0.0686.20E-04−0.0662.50E-03−7.2353.70E-010.7283.40E-02
23040.0964.80E-040.1031.40E-020.1795.60E-03−0.0726.40E-04−0.0722.60E-03−6.053.80E-010.63.50E-02
23140.0914.90E-040.071.40E-020.1845.70E-03−0.0696.60E-04−0.0772.70E-03−6.4433.90E-010.653.60E-02
23240.0965.00E-040.0991.50E-020.2045.90E-03−0.0726.80E-04−0.0832.70E-03−3.7814.00E-010.3773.70E-02
23350.0914.90E-040.0731.40E-020.25.80E-03−0.0686.70E-04−0.0842.70E-03−4.1013.90E-010.4163.60E-02
23450.0865.80E-040.0511.70E-020.226.80E-03−0.0727.90E-04−0.0963.20E-03−5.9354.70E-010.5924.30E-02
23550.096.00E-040.0721.70E-020.2297.00E-03−0.0758.10E-04−0.0993.30E-03−4.634.80E-010.4724.40E-02
23650.0865.90E-040.0371.70E-020.2326.90E-03−0.0718.00E-04−0.1013.20E-03−4.4534.70E-010.4464.40E-02
23750.0846.80E-040.032.00E-020.2537.90E-03−0.0769.20E-04−0.1143.70E-03−5.1125.40E-010.4995.00E-02
23850.0888.30E-040.042.40E-020.2759.70E-03−0.0811.10E-03−0.1284.50E-03−5.5916.60E-010.6036.10E-02
23950.0958.10E-040.0522.40E-020.2879.50E-03−0.0771.10E-03−0.1294.40E-030.0976.50E-010.0856.00E-02

Appendix B

The mean and standard deviation from the estimated EO-1 Hyperion TOA reflectance of 64 images across 196 bands are shown in Figure A1 below. Note that each color code was randomly chosen for each wavelength.
Figure A1. EO-1 Hyperion TOA reflectance mean and standard deviation of 64 images with respect to 196 bands.
Figure A1. EO-1 Hyperion TOA reflectance mean and standard deviation of 64 images with respect to 196 bands.
Remotesensing 15 02141 g0a1

References

  1. Chaity, M.D.; Kaewmanee, M.; Leigh, L.; Teixeira Pinto, C. Hyperspectral Empirical Absolute Calibration Model Using Libya 4 Pseudo Invariant Calibration Site. Remote Sens. 2021, 13, 1538. [Google Scholar] [CrossRef]
  2. Bacour, C.; Briottet, X.; Bréon, F.M.; Viallefont-Robinet, F.; Bouvet, M. Revisiting Pseudo Invariant Calibration Sites (PICS) Over Sand Deserts for Vicarious Calibration of Optical Imagers at 20 km and 100 km Scales. Remote Sens. 2019, 11, 1166. [Google Scholar] [CrossRef]
  3. Raut, B.; Kaewmanee, M.; Angal, A.; Xiong, X.; Helder, D. Empirical Absolute Calibration Model for Multiple Pseudo-Invariant Calibration Sites. Remote Sens. 2019, 11, 1105. [Google Scholar] [CrossRef]
  4. Mishra, N.; Helder, D.; Angal, A.; Choi, J.; Xiong, X. Absolute Calibration of Optical Satellite Sensors Using Libya 4 Pseudo Invariant Calibration Site. Remote Sens. 2014, 6, 1327–1346. [Google Scholar] [CrossRef]
  5. Fajardo Rueda, J.; Leigh, L.; Pinto, C.; Kaewmanee, M.; Helder, D. Classification and Evaluation of Extended PICS (EPICS) on a Global Scale for Calibration and Stability Monitoring of Optical Satellite Sensors. Remote Sens. 2021, 13, 3350. [Google Scholar] [CrossRef]
  6. Shrestha, M.; Leigh, L.; Helder, D. Classification of North Africa for Use as an Extended Pseudo Invariant Calibration Sites (EPICS) for Radiometric Calibration and Stability Monitoring of Optical Satellite Sensors. Remote Sens. 2019, 11, 875. [Google Scholar] [CrossRef]
  7. Markham, B.; Barsi, J.; Kvaran, G.; Ong, L.; Kaita, E.; Biggar, S.; Czapla-Myers, J.; Mishra, N.; Helder, D. Landsat-8 Operational Land Imager Radiometric Calibration and Stability. Remote Sens. 2014, 6, 12275–12308. [Google Scholar] [CrossRef]
  8. Helder, D.; Basnet, B.; Morstad, D. Optimized identification of worldwide radiometric pseudo-invariant calibration sites. Can. J. Remote Sens. 2010, 36, 527–539. [Google Scholar] [CrossRef]
  9. Chander, G.; Angal, A.; Xiong, X.J.; Helder, D.L.; Mishra, N.; Choi, T.J.; Wu, A. Preliminary assessment of several parameters to measure and compare usefulness of the CEOS reference pseudo-invariant calibration sites. Sens. Syst.-Next-Gener. Satell. XIV 2010, 7826, 678–689. [Google Scholar] [CrossRef]
  10. Cosnefroy, H.; Leroy, M.; Briottet, X. Selection and characterization of Saharan and Arabian desert sites for the calibration of optical satellite sensors. Remote Sens. Environ. 1996, 58, 101–114. [Google Scholar] [CrossRef]
  11. USGS EROS Archive—Committee on Earth Observation Satellites (CEOS) Legacy—Calibration/Validation Test Sites. Available online: http://xxx.lanl.gov/abs/https://www.usgs.gov/centers/eros/science/usgs-eros-archive-committee-Earth-observation-satellites-ceos-legacy?qt-science_center_objects=0#qt-science_center_objects (accessed on 5 March 2022).
  12. Helder, D.; Thome, K.J.; Mishra, N.; Chander, G.; Xiong, X.; Angal, A.; Choi, T. Absolute Radiometric Calibration of Landsat Using a Pseudo Invariant Calibration Site. IEEE Trans. Geosci. Remote Sens. 2013, 51, 1360–1369. [Google Scholar] [CrossRef]
  13. Tuli, F.T.Z.; Pinto, C.T.; Angal, A.; Xiong, X.; Helder, D. New Approach for Temporal Stability Evaluation of Pseudo-Invariant Calibration Sites (PICS). Remote Sens. 2019, 11, 1502. [Google Scholar] [CrossRef]
  14. Vuppula, H. Normalization of Pseudo-Invariant Calibration Sites for Increasing the Temporal Resolution and Long-Term Trending. Master’s Thesis, South Dakota State University, Brookings, SD, USA, 2017. [Google Scholar]
  15. Shah, R.; Leigh, L.; Kaewmanee, M.; Pinto, C.T. Validation of Expanded Trend-to-Trend Cross-Calibration Technique and Its Application to Global Scale. Remote Sens. 2022, 14, 6216. [Google Scholar] [CrossRef]
  16. Govaerts, Y.; Clerici, M. Evaluation of Radiative Transfer Simulations Over Bright Desert Calibration Sites. IEEE Trans. Geosci. Remote Sens. 2004, 42, 176–187. [Google Scholar] [CrossRef]
  17. Govaerts, Y.; Adriaensen, S.; Sterckx, S. Optical Sensor CAlibration using simulated radiances over desert sites. In Proceedings of the 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 22–27 July 2012. [Google Scholar] [CrossRef]
  18. Bhatt, R.; Doelling, D.R.; Morstad, D.; Scarino, B.R.; Gopalan, A. Desert-Based Absolute Calibration of Successive Geostationary Visible Sensors Using a Daily Exoatmospheric Radiance Model. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3670–3682. [Google Scholar] [CrossRef]
  19. Bhatt, R.; Doelling, D.R.; Wu, A.; Xiong, X.; Scarino, B.R.; Haney, C.O.; Gopalan, A. Initial Stability Assessment of S-NPP VIIRS Reflective Solar Band Calibration Using Invariant Desert and Deep Convective Cloud Targets. Remote Sens. 2014, 6, 2809–2826. [Google Scholar] [CrossRef]
  20. Kaewmanee, M.; Helder, D. Refined Absolute PICS Calibration Model Over Libya-4 Using 784 Sentinel2A and Landsat 8 Collection-1 Data for Validation. In Proceedings of the PECROA 20, Sioux Falls, SD, USA, 13 November 2017. [Google Scholar]
  21. Landsat Collection-1 Level-1 Product. Available online: http://xxx.lanl.gov/abs/https://docslib.org/doc/7625025/Landsat-collection-1-level-1-product-definition (accessed on 10 March 2022).
  22. Farhad, M.M.; Kaewmanee, M.; Leigh, L.; Helder, D.L. Radiometric Cross Calibration and Validation Using 4 Angle BRDF Model between Landsat 8 and Sentinel 2A. Remote Sens. 2020, 12, 806. [Google Scholar] [CrossRef]
  23. Leigh, L.; Shrestha, M.; Hasan, N.; Kaewmanee, M. Classification of North Africa for Use as an Extended Pseudo Invariant Calibration Site for Radiometric Calibration and Stability Monitoring of Optical Satellite Sensors. In Proceedings of the CALCON 2019, Utah State University, Logan, UT, USA, 19–21 September 2019. [Google Scholar]
  24. Goward, S.N.; Masek, J.G.; Williams, D.L.; Irons, J.R.; Thompson, R. The Landsat 7 mission: Terrestrial research and applications for the 21st century. Remote Sens. Environ. 2001, 78, 3–12. [Google Scholar] [CrossRef]
  25. Andrefouet, S.; Bindschadler, R.; de Colstoun, E.B.; Choate, M.; Chomentowski, W.; Christopherson, J.; Doorn, B.; Hall, D.; Holifield, C.; Howard, S. Preliminary Assessment of the Value Of Landsat-7 Etm+ Data Following Scan Line Corrector Malfunction; EROS Data Center, US Geological Survey: Garretson, SD, USA, 2003.
  26. USGS Landsat-7 Mission. Available online: http://xxx.lanl.gov/abs/https://www.usgs.gov/Landsat-missions/Landsat-7 (accessed on 10 April 2022).
  27. Landsat Collection 2 Level-1 Data. Available online: http://xxx.lanl.gov/abs/https://www.usgs.gov/Landsat-missions/Landsat-collection-2-level-1-data (accessed on 6 March 2022).
  28. Landsat Mission: Landsat-9. Available online: http://xxx.lanl.gov/abs/https://www.usgs.gov/Landsat-missions/Landsat-9 (accessed on 5 August 2022).
  29. Landsat 8 Data Users Handbook. Available online: http://xxx.lanl.gov/abs/https://www.usgs.gov/Landsat-missions/Landsat-8-data-users-handbook (accessed on 2 September 2021).
  30. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  31. Li, J.; Roy, D.P. A Global Analysis of Sentinel-2A, Sentinel-2B and Landsat-8 Data Revisit Intervals and Implications for Terrestrial Monitoring. Remote Sens. 2017, 9, 902. [Google Scholar] [CrossRef]
  32. Revel, C.; Lonjou, V.; Marcq, S.; Desjardins, C.; Fougnie, B.; Luche, C.; Guilleminot, N.; Lacamp, A.S.; Lourme, E.; Miquel, C.; et al. Sentinel-2A and 2B absolute calibration monitoring. Eur. J. Remote Sens. 2019, 52, 122–137. [Google Scholar] [CrossRef]
  33. Ungar, S.; Pearlman, J.; Mendenhall, J.; Reuter, D. Overview of the Earth Observing One (EO-1) mission. IEEE Trans. Geosci. Remote Sens. 2003, 41, 1149–1159. [Google Scholar] [CrossRef]
  34. Jing, X.; Leigh, L.; Helder, D.; Teixeira Pinto, C.; Aaron, D. Lifetime Absolute Calibration of the EO-1 Hyperion Sensor and its Validation. IEEE Trans. Geosci. Remote Sens. 2019, 57, 9466–9475. [Google Scholar] [CrossRef]
  35. Franks, S.; Neigh, C.S.R.; Campbell, P.K.; Sun, G.; Yao, T.; Zhang, Q.; Huemmrich, K.F.; Middleton, E.M.; Ungar, S.G.; Frye, S.W. EO-1 Data Quality and Sensor Stability with Changing Orbital Precession at the End of a 16 Year Mission. Remote Sens. 2017, 9, 412. [Google Scholar] [CrossRef]
  36. Micijevic, E.; Haque, M.O.; Mishra, N. Radiometric calibration updates to the Landsat collection. In Earth Observing Systems XXI; Butler, J.J., Xiong, X.J., Gu, X., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, WA, USA, 2016; Volume 9972, pp. 108–119. [Google Scholar] [CrossRef]
  37. Landsat Collection 2 Quality Assessment Bands. Available online: https://www.usgs.gov/Landsat-missions/Landsat-collection-2-quality-assessment-bands (accessed on 6 March 2022).
  38. Available online: https://www.mathworks.com/help/matlab/ref/interp1.html (accessed on 18 June 2022).
  39. Barsi, J.A.; Alhammoud, B.; Czapla-Myers, J.; Gascon, F.; Haque, M.O.; Kaewmanee, M.; Leigh, L.; Markham, B.L. Sentinel-2A MSI and Landsat-8 OLI radiometric cross comparison over desert sites. Eur. J. Remote Sens. 2018, 51, 822–837. [Google Scholar] [CrossRef]
  40. Gross, G.; Helder, D.; Leigh, L. Extended Cross-Calibration Analysis Using Data from the Landsat 8 and 9 Underfly Event. Remote Sens. 2023, 15, 1788. [Google Scholar] [CrossRef]
  41. Gross, G.; Helder, D.; Begeman, C.; Leigh, L.; Kaewmanee, M.; Shah, R. Initial Cross-Calibration of Landsat 8 and Landsat 9 Using the Simultaneous Underfly Event. Remote Sens. 2022, 14, 2418. [Google Scholar] [CrossRef]
  42. Sterckx, S.; Wolters, E. Radiometric Top-of-Atmosphere Reflectance Consistency Assessment for Landsat 8/OLI, Sentinel-2/MSI, PROBA-V, and DEIMOS-1 over Libya-4 and RadCalNet Calibration Sites. Remote Sens. 2019, 11, 2253. [Google Scholar] [CrossRef]
  43. Cui, Z.; Kerekes, J.P. Impact of Wavelength Shift in Relative Spectral Response at High Angles of Incidence in Landsat-8 Operational Land Imager and Future Landsat Design Concepts. IEEE Trans. Geosci. Remote Sens. 2018, 56, 5873–5883. [Google Scholar] [CrossRef]
Figure 1. Flowchart showing study area selection and sensor pre-processing for the DAHAC model’s development.
Figure 1. Flowchart showing study area selection and sensor pre-processing for the DAHAC model’s development.
Remotesensing 15 02141 g001
Figure 2. GC-36 stable dark target pixels in red color. WRS-2 path/row images footprint for L8 in the different continents: (a) Middle East, (b) South Africa, and (c) North Africa.
Figure 2. GC-36 stable dark target pixels in red color. WRS-2 path/row images footprint for L8 in the different continents: (a) Middle East, (b) South Africa, and (c) North Africa.
Remotesensing 15 02141 g002
Figure 3. L8 dark target TOA reflectance of eight WRS-2 paths/rows for seven spectral bands.
Figure 3. L8 dark target TOA reflectance of eight WRS-2 paths/rows for seven spectral bands.
Remotesensing 15 02141 g003
Figure 4. TOA reflectance comparison between selected dark target (L8) and GC-36 and red dotted lines showing the 3 σ TOA reflectance of GC-36.
Figure 4. TOA reflectance comparison between selected dark target (L8) and GC-36 and red dotted lines showing the 3 σ TOA reflectance of GC-36.
Remotesensing 15 02141 g004
Figure 5. TOA reflectance of twelve tiles for S2A spectral bands ranging from 0.04 to 0.18 TOA reflectance.
Figure 5. TOA reflectance of twelve tiles for S2A spectral bands ranging from 0.04 to 0.18 TOA reflectance.
Remotesensing 15 02141 g005
Figure 6. EO-1 Hyperion hyperspectral TOA reflectance for 64 images obtained from the selected stable dark sites.
Figure 6. EO-1 Hyperion hyperspectral TOA reflectance for 64 images obtained from the selected stable dark sites.
Remotesensing 15 02141 g006
Figure 7. Block diagram showing satellite hyperspectral profile generation process.
Figure 7. Block diagram showing satellite hyperspectral profile generation process.
Remotesensing 15 02141 g007
Figure 8. Overall flowchart showing the Dark Hyperspectral Dataset generation process.
Figure 8. Overall flowchart showing the Dark Hyperspectral Dataset generation process.
Remotesensing 15 02141 g008
Figure 9. Showing L8 normalized hyperspectral TOA reflectance (grey curve), L8 multispectral TOA reflectance (orange circle), L8 normalized hyperspectral integrated TOA reflectance (purple circle), and absolute difference (red bar): before relative calibration.
Figure 9. Showing L8 normalized hyperspectral TOA reflectance (grey curve), L8 multispectral TOA reflectance (orange circle), L8 normalized hyperspectral integrated TOA reflectance (purple circle), and absolute difference (red bar): before relative calibration.
Remotesensing 15 02141 g009
Figure 10. Relative gain for L8 and S2A vs. the central wavelength.
Figure 10. Relative gain for L8 and S2A vs. the central wavelength.
Remotesensing 15 02141 g010
Figure 11. Superspectral and hyperspectral (for calibrated Hyperion wavelengths) relative gain.
Figure 11. Superspectral and hyperspectral (for calibrated Hyperion wavelengths) relative gain.
Remotesensing 15 02141 g011
Figure 12. Block diagram showing the relative calibration on the EO-1 Hyperion-processed image.
Figure 12. Block diagram showing the relative calibration on the EO-1 Hyperion-processed image.
Remotesensing 15 02141 g012
Figure 13. Showing L8 normalized hyperspectral TOA reflectance (grey curve), L8 multispectral TOA reflectance (orange circle), L8 normalized hyperspectral integrated TOA reflectance (purple circle), and absolute difference (red bar): after relative calibration.
Figure 13. Showing L8 normalized hyperspectral TOA reflectance (grey curve), L8 multispectral TOA reflectance (orange circle), L8 normalized hyperspectral integrated TOA reflectance (purple circle), and absolute difference (red bar): after relative calibration.
Remotesensing 15 02141 g013
Figure 14. L8-normalized hyperspectral TOA reflectance (grey line), EO-1 hyperspectral TOA reflectance after relative gain correction (blue line), and L8 multispectral TOA reflectance (orange circle).
Figure 14. L8-normalized hyperspectral TOA reflectance (grey line), EO-1 hyperspectral TOA reflectance after relative gain correction (blue line), and L8 multispectral TOA reflectance (orange circle).
Remotesensing 15 02141 g014
Figure 15. S2A-normalized hyperspectral TOA reflectance (grey line), EO-1 hyperspectral TOA reflectance after relative gain correction (blue line), and S2A multispectral TOA reflectance (orange circle).
Figure 15. S2A-normalized hyperspectral TOA reflectance (grey line), EO-1 hyperspectral TOA reflectance after relative gain correction (blue line), and S2A multispectral TOA reflectance (orange circle).
Remotesensing 15 02141 g015
Figure 16. Dark hyperspectral dataset (grey) vs. 4-angle hyperspectral BRDF-model-predicted TOA reflectance profile (blue).
Figure 16. Dark hyperspectral dataset (grey) vs. 4-angle hyperspectral BRDF-model-predicted TOA reflectance profile (blue).
Remotesensing 15 02141 g016
Figure 17. TOA reflectance differences between DaHD vs. 4-angle hyperspectral BRDF model for 196 bands.
Figure 17. TOA reflectance differences between DaHD vs. 4-angle hyperspectral BRDF model for 196 bands.
Remotesensing 15 02141 g017
Figure 18. DAHAC model coefficients values for 196 hyperspectral bands.
Figure 18. DAHAC model coefficients values for 196 hyperspectral bands.
Remotesensing 15 02141 g018
Figure 19. Showing (a) TOA reflectance of 4-angle hyperspectral BRDF model (15 coefficients) vs. DAHAC model (7 coefficients) and (b) histogram with TOA reflectance difference between 4-angle hyperspectral BRDF model vs. DAHAC model TOA reflectance.
Figure 19. Showing (a) TOA reflectance of 4-angle hyperspectral BRDF model (15 coefficients) vs. DAHAC model (7 coefficients) and (b) histogram with TOA reflectance difference between 4-angle hyperspectral BRDF model vs. DAHAC model TOA reflectance.
Remotesensing 15 02141 g019
Figure 20. DAHAC model validation on L8 across seven spectral bands.
Figure 20. DAHAC model validation on L8 across seven spectral bands.
Remotesensing 15 02141 g020
Figure 21. TOA reflectance difference between Landsats’ observed TOA reflectance and the DAHAC model’s response across seven spectral bands, showing the mean TOA reflectance difference and standard deviation (pink).
Figure 21. TOA reflectance difference between Landsats’ observed TOA reflectance and the DAHAC model’s response across seven spectral bands, showing the mean TOA reflectance difference and standard deviation (pink).
Remotesensing 15 02141 g021
Figure 22. DAHAC model validation on S2A for non-Landsat bands.
Figure 22. DAHAC model validation on S2A for non-Landsat bands.
Remotesensing 15 02141 g022
Figure 23. TOA reflectance difference between S2A’s observed information and the DAHAC model’s response across eleven spectral bands.
Figure 23. TOA reflectance difference between S2A’s observed information and the DAHAC model’s response across eleven spectral bands.
Remotesensing 15 02141 g023
Figure 24. Flowchart showing the calculation of DAHAC model’s uncertainty.
Figure 24. Flowchart showing the calculation of DAHAC model’s uncertainty.
Remotesensing 15 02141 g024
Figure 25. DAHAC model’s uncertainty using Monte Carlo simulation at different iteration levels (100–2500).
Figure 25. DAHAC model’s uncertainty using Monte Carlo simulation at different iteration levels (100–2500).
Remotesensing 15 02141 g025
Figure 26. Mean TOA reflectance of the DAHAC model (blue colored) and DAHAC model’s uncertainty in unit reflectance (orange colored) for a 426–2395 nm wavelength range.
Figure 26. Mean TOA reflectance of the DAHAC model (blue colored) and DAHAC model’s uncertainty in unit reflectance (orange colored) for a 426–2395 nm wavelength range.
Remotesensing 15 02141 g026
Figure 27. Double ratio between L8 and L9. Orange dotted lines show a maximum 2% differences between L8 and L9.
Figure 27. Double ratio between L8 and L9. Orange dotted lines show a maximum 2% differences between L8 and L9.
Remotesensing 15 02141 g027
Figure 28. Double ratio between S2A and S2B. Orange dotted lines show a maximum 1% differences between S2A and S2B.
Figure 28. Double ratio between S2A and S2B. Orange dotted lines show a maximum 1% differences between S2A and S2B.
Remotesensing 15 02141 g028
Table 1. Different sensor information used in the DAHAC model’s development.
Table 1. Different sensor information used in the DAHAC model’s development.
InformationSensors
Landsat-7Landsat-8Landsat-9Sentinel-2ASentinel-2BEO-1 Hyperion
No. of Images76411508477545664
Image Dates2000–20222013–20222021–20222015–20222017–20222001–2017
SZA Range (°)20–5820–6020–6015–5115–5423–77
SAA Range (°)36–16035–16035–15831–16431–16570–145
VZA Range (°)0.10–80.03–50.17–80.04–100.09–100.04–25
VAA Range (°)−94–136−178–180−80–118−89–130−162–135−82–98
No. of Sites8 WRS-2
paths/rows
8 WRS-2
paths/rows
8 WRS-2
paths/rows
12 Tiles12 Tiles3 paths/rows
Table 2. Showing WRS-2 paths/rows, surface types, and pixel count across different continents.
Table 2. Showing WRS-2 paths/rows, surface types, and pixel count across different continents.
ContinentsPaths/RowsSurface DescriptionPixel Count
Middle East159/40dark rock73,544
163/37dark rock36,070
168/51dark rock, dark sand2,359,755
170/42volcanic area710,696
South Africa180/75volcanic area, dark rock40,626
181/73dark rock, dark sand21,726
North Africa183/46volcanic area30,009
184/43volcanic area133,950
Table 3. Showing 4-angle hyperspectral BRDF model coefficients for 864 nm wavelength.
Table 3. Showing 4-angle hyperspectral BRDF model coefficients for 864 nm wavelength.
CoefficientEstimateStandard
Error
t-Statisticp-ValueStatistical
Response
Intercept0.146e-04226.200Significant
X 1 −7.59e-183.30e-04−2.3e-141.0Insignificant
Y 1 −2.71e-173.88e-04−6.98e-141.0Insignificant
X 2 7.33e-180.0061.18e-151.0Insignificant
Y 2 1.50e-180.0027.64e-161.0Insignificant
X 1 Y 1 3.24e-189.59e-043.38e-151.0Insignificant
X 1 X 2 0.160.029.081.27e-19Significant
X 1 Y 2 −4.67e-170.005−8.7e-151.0Insignificant
Y 1 X 2 2.79e-170.021.27e-151.0Insignificant
Y 1 Y 2 0.160.00722.325.42e-107Significant
X 2 Y 2 2.45e-180.054.58e-171.0Insignificant
X 1 2 −0.088.13e-04−107.410Significant
Y 1 2 −0.060.003−19.692.75e-84Significant
X 2 2 16.970.48−35.294.32e-253Significant
Y 2 2 1.620.0436.635.69e-271Significant
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Karki, P.B.; Kaewmanee, M.; Leigh, L.; Pinto, C.T. The Development of Dark Hyperspectral Absolute Calibration Model Using Extended Pseudo Invariant Calibration Sites at a Global Scale: Dark EPICS-Global. Remote Sens. 2023, 15, 2141. https://doi.org/10.3390/rs15082141

AMA Style

Karki PB, Kaewmanee M, Leigh L, Pinto CT. The Development of Dark Hyperspectral Absolute Calibration Model Using Extended Pseudo Invariant Calibration Sites at a Global Scale: Dark EPICS-Global. Remote Sensing. 2023; 15(8):2141. https://doi.org/10.3390/rs15082141

Chicago/Turabian Style

Karki, Padam Bahadur, Morakot Kaewmanee, Larry Leigh, and Cibele Teixeira Pinto. 2023. "The Development of Dark Hyperspectral Absolute Calibration Model Using Extended Pseudo Invariant Calibration Sites at a Global Scale: Dark EPICS-Global" Remote Sensing 15, no. 8: 2141. https://doi.org/10.3390/rs15082141

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop