[INFO] Introduction of additional Radiometric Offset in PB04.00 products

Dear friends,

I write this to underscore the information in the Tables here
Processing Baseline - Sentinel Online and in the L1C and L2A Data Product Quality Reports - Sentinel Online

Along with several other evolutions, Processing Baseline 04.00 L1C (TOA) and L2A (BOA) products (TTO 25th of January, 2022) now contain an Offset in the metadata. This offset has been included to accomodate ‘noisy’ pixels over dark surfaces, that can sometimes present negative radiometric values at L1B and L1C.

Prior to PB 04.00, these negative values were clamped to a predefined range (e.g. [1-32767]), thereby leading to a loss of information. In order to avoid this, and following an analysis by expertise, a radiometric offset has been introduced which shifts the range. The conclusions of the analysis resulted in the following specifications for the Offset:

  • At L1B, RADIO_ADD_OFFSET = -100
  • At L1C, RADIO_ADD_OFFSET = -1000
  • At L2A, RADIO_ADD_OFFSET = -1000

SNAP is User-configurable for this change; the Radiometric Offset can be added in SNAP by selection of the relevant radio button (Tools > Options > Sentinel-2 Reader > Add negative radiometric offset (L1C -L2A) ):

image

Additional information can be found in the [Sentinel-2 Products Specification Document (PSD)](Sentinel-2 Products Specification Document (PSD) - Sentinel Online .)

I hope this information is of benefit to you.

Cheers

Jan

OPT-MPC S2 Technical Manager

4 Likes

Hi @Jan,
I posted a question yesterday (Clamp negative reflectance to 0 (NODATA) for Baseline 04.00 with new OFFSET correction?) that perhaps you could help me with (pert of it repeated below).

While comparing images between the current Baseline 04.00 and earlier Baseline versions I subtract the OFFSET value Baseline 04.00 images depending on product. This sometimes leads to, as expected, negative values in the resulting raster.

My question is: What is the correct way to present these values if clamping them to a floor of 0? Should they be regarded as NODATA (value 0) - Implying missing or faulty data? Or should they be viewed as a minimum value (like 1)? The reasoning being that these have a registered value that is out of bounds, but assumed to be approximately 0.

Thanks!

Hi @Eddie
I would suggest that they should be regarded as NO DATA (value 0). This would at least flag them as unreliable. This is based on my understanding of negative reflectance values being linked to the Atmospheric Correction processing step. And thus they are an unwanted contribution to the output. So making them 0 identifies them as non-nominal, and helps serve as an avoidance measure. Or a flag of doubt. Giving them a value of 1 would potentially attribute them some relevance, and may serve to confuse.

I hope this helps.

Cheers

Jan

OPT-MPC S2 Technical Manager

Thanks for your swift response @Jan,
thats great!
Is this perhaps something that could be clarified in the “Sentinel-2 Products Specification Document” or some related documentation?

Hi @Eddie
No worries. We are currently undertaking a wholesale update of Sentinel Online, so I will suggest that the content is updated to reflect these changes.

As far as the PSD goes, it’s really more a description of the Product levels and their packaging (and associated schemas and metadata), so the specifics of the Radiometric Offset aren’t in that remit. But I’ll mention it to the L2A expertise to see if it can be included in the ATBD .

Cheers

Jan

1 Like

Hi Jan, I could not find this information elsewhere.

Is the RADIO_ADD_OFFSET value always the same for all bands or is it variable meaning it is mandatory to check the metadata for every scene and extract the exavt value.

Cheers

Hello @alesinar

The RADIO_ADD_OFFSET is fixed:

At L1B, RADIO_ADD_OFFSET = -100
At L1C, RADIO_ADD_OFFSET = -1000
At L2A, RADIO_ADD_OFFSET = -1000

In the L1C product it is available from the DS metadata at <Radiometric_Offset_List>
In the L2A product it is available from the DS metadata at <BOA_ADD_OFFSET_VALUES_LIST>

Page 444 of the [PSD] identifies its use in converting DN to physical values.

I hope this helps you.

Cheers

Jan

Great, thank you!

1 Like

Hello @Jan. I am working with using an NDVI threshold we identified via fieldwork to map bare ground. To be consistent with the previous years (prior to v0400 offset) I need to get the post-v0400 images to work with our threshold. I have tried loading a L2A scene into SNAP with the negative radiometric offset ticked and then exported the bands i need (4 & 8a) as TIFs but unfortunately it does not make any difference to the values in the images. What am I doing wrong here?

All the best,

Hello @tyhayward
Can you let me know the product details please? I’ll download it and try from this end.

Thanks in advance

Jan

Hi @Jan thanks for getting back.
The product i attempted this within SNAP was:
S2B_MSIL2A_20230119T112309_N0509_R037_T30UWC_20230119T125155

It would be good if you could share how this methodology is applied as I think that is where I might be going wrong.

I have since done a simple band calculation of -1000 to bands B04 & B8a and then changed any negative values to 0. Doing this I have achieved more reasonable NDVI responses but I am still unsure as to if this was correct or not. This was performed in ArcGIS PRO.

All the best,

Hi @Jan did you manage to take a look at this? It would be a great help for the work that we do to get a functioning methodology in place.

all the best,

Hello @tyhayward

I’m sorry for not getting back to you sooner. I have looked at profile plots over in your AOI for a couple of L2As and they both exhibit similar values. One thing I did not however, is that the built up areas are sometimes rendered as cloud shadow in the Scene Classification Layer (SCL) and this is also seen in the WVP quality layer too:

Cheers

Jan