GEOSAT 2 Atmospherically Corrected Images: algorithm validation

: Solar radiation reflected by the Earth´s surface to satellite sensors is modified by its interaction with the atmosphere. The application of atmospheric correction of optical satellite imagery is an essential and needed pre-processing tool for modeling biophysical variables, multi-temporal analysis, and digital classification processes. As a result, true surface reflectance values are obtained without atmosphere influence. To assess this process, GEOSAT (part of the ESA’s Third -Party Mission Programme) performs an optimization of the GEOSAT 2 very high resolution (VHR) multi-spectral imagery adapting the well-known 6S model to the different wavelengths covered by the GEOSAT 2 spectral bands (VHR, PAN). 6S model predicts surface reflectance (BOA) using information from the apparent reflectance (TOA) captured by the satellite sensor and the corresponding atmospheric conditions. To perform the atmospheric correction (AC), both the configuration of the atmosphere at the time of capture and the conditions of scene pointing and luminosity, must be considered. The first is mainly determined by three values: water vapor, ozone, and the number of air-suspended particles (aerosols). For the latter, the geometry of the scene, as well as the respective sun and sensor observation positions are the values to be considered. To validate the resultant GEOSAT 2 AC images, obtained from applying the GEOSAT atmospheric correction algorithm, different common areas between these and Sentinel-2 L2A products have been selected. Then, band-by-band (R, G, B & NIR) operations, such as calculation of the mean square error (RMSE) and a regression analysis were performed. Then, spectral profiles for the three generic land coverages (vegetation, soil and water) were also gathered over the spectral range of GEOSAT 2 and S2 corresponding bands. The outcomes, once analyzed, lead us to conclude that the results obtained by applying the promising GEOSAT AC algorithm are satisfactory and seem to correctly estimate BOA reflectance values for vegetation and water coverages. To extend the study and improve the result ground reflectance values will be required.


Introduction
Radiation leaving the Earth's surface experiences an important interaction with the atmosphere before it is registered by a satellite sensor.There are many ways to compensate for atmospheric contributions to an Earth image observed by an optical satellite sensor.They vary from empirical methods, which modify the brightness of each sensor spectral band, to more robust procedures based on complex and robust physical models describing the atmospheric radiation paths and considering the most contributing atmospheric gases and aerosols.Within the latter, Radiative Transfer models are the most appropriate methods, as they are able to effectively couple and compensate the mentioned atmospheric effects.For several decades, efforts have been made to implement computer codes that simulate and correct for atmospheric disturbances on optical observations.This is the case for the MODTRAN (MODerate resolution atmospheric TRANsmission) [1].Several computer packages have been implemented based on the MODTRAN philosophy, such as the 6S (Second Simulation of a Satellite Signal in the Solar Spectrum) [2], which focusses on correcting for atmospheric effects on airborne and satellite optical images.An important feature of the 6S code is that it can be optimized and customized for specific Earth observing systems, which is the case for GEOSAT 2 imagery.In addition to the already demonstrated good performance of these sophisticated atmospheric correction packages, their application to a particular study and imagery must always be validated in order to assess the uncertainty of the correction process.This study presents results of the optimization of the 6S code for producing GEOSAT 2 level 2 (DE2) products and the respective evaluation of the corrected imagery.For this purpose, processed GE-OSAT 2 DE2 and SENTINEL 2A MSIL2A products with close overpasses have been used to preliminary assess the suitability of the GEOSAT 2 atmospheric correction processor.

Model Description
GEOSAT 2 Surface Reflectance (SR) product is derived from the standard radiance product being processed first to TOA reflectance and then atmospherically corrected to BOA reflectance using the well-known radiative transfer (RT) method 6S [2,3,4].This algorithm has been configured and applied to the different wavelengths covered by the GE-OSAT 2spectral bands (VHR, PAN).For the sake of simplicity and speed we have considered only the bidirectional reflectance (BDR), considering a uniform, or Lambertian, surface.Based on the radiative transfer theory and assuming that the target is a Lambertian surface, the surface reflectance (), in terms of the at-sensor radiance (), can be expressed as follows: The additional parameters of this function: the path radiance (  ), the transmissivity () from surface to satellite, taken into account the transmissivities, caused by both absorption and scattering, the direct solar irradiance (  ) and the diffuse solar irradiance (  ), are all obtained by executing the 6S model.
A typical orthorectified GEOSAT 2 image can have more than 30 million pixels.Thus, considering a pixel-wise 6S model execution results in excessive time demanding.To reduce the time needed to retrieve BOA reflectance products, an optimization approach, based on lookup tables (LUT), has been performed on GEOSAT 2 imagery.Thus, the  3   2  values are derived from the spatial-temporal table defined in [5] which provides us a valid approximation.On the other hand, MERRA-2 service from the Global Modeling and Assimilation Office (of the National Aeronautics and Space Administration -NASA), is used to retrieve closer AOT values to GEOSAT 2image sensing time.Concretely, the hourly (M2T1NXAE [6]) and monthly (M2TMNXAER [7]) aerosol optical thickness at 550 nm are accessed from this NASA service, moreover acquisition geometry angles are taken to center time, and surface elevation is determined through the SRTM 1 Arc-Second Global (Earth Resources Observation and Science (EROS) Center (USGS), 2018) (30m) elevation data from USGS/NASA [8].

Model Validation
To assess the atmospheric correction method for GEOSAT 2 (DE2) products, we utilized atmospherically corrected Sentinel 2 (S2) mission level 2A products.The evaluation procedure involved selecting a common area between both products and resample DE2 products to coincide with the S2 MSI spectral bands resolution.The bands and their spectral ranges are as follows: For Sentinel 2, there are slight differences in the spectral sensitivity of the MSI instruments aboard the S2A and S2B missions., but the spectral response sensitivity functions can be considered practically identical between both missions for the bands considered in this study.Regarding the spectral sensitivity for the DE2 mission bands, no data is available for this study.However, according to the below table (Table 1), the central bands of the DE2 mission are nearly coincident with those of S2, as well as their spectral widths.In this context, for the forthcoming comparison, it is assumed that analogous spectral bands are being evaluated between DE2 and S2.For the validation of the algorithm, several procedures have been considered, • First, the root mean square error (RMSEbi) has been found, band by band, between the values of the product DE2 and S2. • The second procedure relied on a regression analysis between the corresponding bands of both products, fitting a linear function to the two datasets.In this analysis, the regression coefficient (R 2 ) and the fitting error (RMSE) are obtained, where the latter can be interpreted as a deviation from the fitted function.

•
Thirdly, for generic land cover types such as vegetation, soil, and water, spectral profiles are acquired at specific positions within the study area, and these profiles are depicted graphically.

Overal validation: Brazil
Figure 1 displays the location of the DE2 product, acquired on November 25, 2021 (20211125), in relation to the S2A product acquired on the same day (Table 2).The elapsed time between the two is approximately 39 minutes.The relative location between both products is depicted in Figure 1.a, whereas Figure 1.b illustrates the area observed by the S2A product and the atmospheric conditions at the time of observation for both products.According to the provided information, DE2 product data have been acquired with an almost nadiral geometry (roll=3, pitch=0).Nevertheless, it is noted that for the S2 missions, the angle of observation or incidence over the scene can vary between 2 degrees and 12 degrees, depending on the observed zone, being the average angle in this zone between 2.6 o and 8.5 o .DE2 and S2A products were acquired with similar viewing geometries.It is worth to note, that the common observed area in both images lies within the lowest incidence angles.
After cropping out the common areas between both products, two images are obtained with the same number of rows and columns.The resulting images are shown in Figure 2. It should be noted that the product DE2 was resampled to the spatial resolution (10 m) corresponding to the (VIS-NIR) bands of the MSI instrument.It was not possible to select a larger area due to the atmospheric conditions that affect both products.For this area the RMSEbi obtained are shown in the table below: In general, it can be observed that the spectral differences are very low for the first three bands.Although the greatest differences are observed between bands B4_DE2 and B8_S2B, they remain relatively low.In spectral terms, these differences can be considered suitable for all bands in the two images.
In the second procedure, based on a regression analysis between the bands indicated in Table 3, the regression coefficients (R 2 ) and the deviations or RMSEaj of the fit, as shown in Table 4, were obtained.Subsequently, Figure 3 displays the graphs of these regressions.Except for the blue bands (DE2 and S2A), the results obtained from the regression analysis are generally satisfactory.All coefficients are above 0.8 (R 2 ≥ 0.8), except for the two above mentioned blue bands.The results from both methods are consistent with each other.Furthermore, according to Table 3, the bands that fit the best are Band 4 of DE2 and Band 8 of S2A (R 2 = 0.98), which, based on the obtained result, would be practically identical.Given the outcomes observed for Bands 1 and 2 of DE2 and S2A, respectively, it would be advisable to extend the analysis to a new area within the same products.The third evaluation has been conducted through the measurement and comparison of spectral profiles.Three cases have been selected for vegetation and soil, while for water, only two positions could be identified, the graphic representation shown, in the horizontal axe the spectral resolution (λ) in micrometers (µ m), and in the vertical axe the reflectance value between 0,1.The results are presented in the supplementary materials.

4.Results
For the selected analysis area, based on the obtained spectral profiles, the following considerations are made.There are no significant differences detected among the various land covers, and a good agreement is observed among all profiles, including vegetation, soils, and water.The differences typically range between 1% and 2%, which is quite acceptable for this type of comparison.For these three types of covers and depending on the considered spectral band, these differences are practically imperceptible.
Water appears to maintain a satisfactory spectral behavior.However, due to the scarcity of representative water bodies, only two spectral profiles could be measured.As a result, unlike the results obtained for vegetation and soils, these findings cannot be considered definitive.

5.Conclusions
This study shows the preliminary results obtained using the atmospheric correction processor for GEOSAT 2 products based on the 6S code.Although the spectral resolutions between the GEOSAT 2 and the MSI systems are slightly different, the similarity metrics used provide values that confirm the good performance of this processor.Only the blue

Citation:
To be added by editorial staff during production.Academic Editor: Firstname Lastname Published: date Copyright: © 2023 by the authors.Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/license s/by/4.0/).

Table 2 .
Product Id for DE2 and S2A.