Two successive operations are undertaken in order to obtain a Level-1C product from Level-1B products:
Resampling comprises five main steps:
The masks of saturated areas, no-data and technical quality are also resampled in the target geometry. The resampling can process all or some spectral bands (all by default).
Tiling ModuleThe goal of this module is to select the list of predefined tiles which will be produced. The Earth's surface is split into several adjacent tiles in a given representation (UTM/WGS84). Some tiles may overlap between them. Each tile is defined by:
A list of tiles which intersect the footprint defined in the metadata is determined. If the intersection is void, the empty tiles are set to no-data by the resampling process.
The auxiliary data geometries are computed by this module.
The resampling function is intended to obtain radiometric information in the target geometry. This process is in two steps:
The geometric transformation allows linking of the points of the target image with the points of the initial image. For SENTINEL-2, a resampling grid is used and computed for each spectral band and each detector. For each tile, there are 13 by 12 resampling grids.
Figure 1: Resampling Grid
The values of the resampling grid are computed for the footprint of the detector plus a small margin by extrapolation of the viewing directions. For each point of the orthoimage, the coordinates of the corresponding pixel in the native image are calculated by bi-linear interpolation of the four adjacent nodes. The grid computation is divided into three steps:
Also, the resampling grid provides complementary information:
A point on the target image can have no homologous point in the initial image. To avoid this, the initial image is projected on the target geometry to determine its footprint. Only pixels included in the interest area of the target geometry are interpolated. The footprint of two detectors within the same spectral band may also overlap. The useful area of each detector is cropped at the column defined in the GIPP inter-detector of the overlapping area for each line (the middle column by default).
Interpolation estimates the radiance values of the target point, knowing the radiance of neighbouring pixels. A linear algorithm is used for this operation with B-spline functions.
TOA Reflectance Computation
The numeric digital counts (CN) of each pixel image (i,j) and each spectral band (k) are converted in TOA reflectance (ρ). This conversion takes into account the equivalent extra-terrestrial solar spectrum (Es), the incoming solar direction defined by its zenith angle (θs) for each pixel of the image and the absolute calibration (Ak) of the instrument MSI.
The conversion equation is:
Equation 1: Top of Atmosphere conversion
Equation 2: Earth Sun distance
t is the Julian Day corresponding to the acquisition date (reference day: 01/01/1950).
0.01673 is the Earth orbit eccentricity.
0.0172 is the Earth angular velocity (radians/day).
The parameters Ak and Es are provided by the GIPP and are also included in the ancillary data of the Level-1 products.
The sun zenith angles are determined at this level too. A sun angle grid is computed by regularly down-sampling the target geometry (Level-1C tile). The cosine of the zenith angle θs is defined at each point of the grid using the ground coordinates and the datation of the corresponding pixel acquisition. The azimuth angle is not processed here.
Reflectance, often between 0 and 1, is converted into integer values, in order to preserve the dynamic range of the data by applying a fixed coefficient (1000 by default).
For Level-1C products, the noise model is also adapted to the new range of TOA reflectance radiometric values. The parameters (Ak,Es) of the Level-1B noise model are corrected using Equation 1 by replacing CNk with the Level-1B parameters Ak1B and Es1B.
Computation of the Masks
The last operation to obtain the final Level-1C product is computation of the cloud masks and the land/water mask for a tile in the Level-1C geometry. The cloud mask specifies cloud type: cirrus or opaque cloud. The computation of these masks is explained and detailed in the Masks section.
Level-1B technical masks (defective, no_data pixels, etc) are then reprojected to the Level-1C geometry.