Radiometric conversion of sentinel 2 images

Optical Toolbox
#14

Thank you very much kiryl23

#15

Hi All,
I am bit confused on the radiometric convertion. Is this on for reflectance value? What is the different with sen2cor approach?
Thanks

#16

Hi Ben,
For this band math, should it be done before using the sen2cor or after?
In my research, I will using S2 and Landsat-8 for data fusion.
Regards

#17

Hi Adhwa,

A lot has changed since I posted that answer. Maybe I can make your life easier. What reflectance product are you trying to achieve and maybe I can answer your question more thoroughly?

Ben

#18

Adhwa,

Not sure if you got my email. Anyway, Sentinel-2A images are in JP2000 format. When importing in a GIS or matrix laboratory, divide each band by the scaling factor (10,000) to retrieve TOA reflectance (rTOA as I have posted in the formula above). Use rTOA to retrieve TOA Radiance in the equation above. Hope this helps. Sorry it took me a while to reply.

Ben

#19

Hi Ben,
Thank you for the information. Actually, I am bit confused on the radiometric convertion and sen2cor. Since I need the reflectance value for the fusion between S2 and L8. I wonder which is more practical to be used ( radiometric convertion and sen2cor) and whether it is the same process or can be integrated.
For example L8, I have already process up to atmospheric correction and convert it to reflectance value using band math. Is it the same process for S2?
Regards

#20

Hi Adhwa,
sen2cor is a plugin designed for transforming Sentinel-2 MSI data on L1C (Top-of-atmosphere (TOA) reflectance) to L2A (Bottom-of-atmosphere (BOA) reflectance). So, when you want to perform atmospheric correction to the Sentinel-2 data in sen2cor, you don’t have to transform it to TOA radiance first. So, basically, you just have to input the image file and let the plugin works.

A different steps is required when you want to perform atmospheric correction with another algorithm (e.g. FLAASH, which is built-in in commercial software such as ENVI). As the algorithm needs TOA radiance data as the input, you have to tranform the data first, which is from TOA reflectance to TOA radiance, using equation given by Igor above (also see Distance of Earth-Sun). Performing atmospheric correction (e.g. FLAASH) allows you to have Surface reflectance data as the output.

#21

Hi Fifauzan,
Thank you for the explanation. I think I know what to do now.
Regards,

#22

Hello @fifauzan ,

Please, can i do the atmospheric correction with the ENVI software with Flaash algorithm

thank you in advance.

#23

Hi @ELHACHIMI ,
Sorry for the late reply. This paper might help: https://www.researchgate.net/publication/307959143_Comparison_of_Physically_and_Image_Based_Atmospheric_Correction_Methods_for_Sentinel-2_Satellite_Imagery?ev=prf_high

1 Like
#24

Hi @fifauzan,
Thank you very much for your help, this is a very important job for me, and this work will help me a lot :slight_smile:.

1 Like
#27

Hello @safah_94 ,
Yes my email is elhachimijaouad@gmail.com

#28

@ELHACHIMI ok thanks. please looked your email

#29

Hello all. . .

Has anyone ever processed Sentinel 2A using ENVI software ???

#30

Maybe you ask this in the support forum of ENVI?
http://harrisgeospatial.com/Support/Forums.aspx

#31

okk. . . thank u @marpet

#32

Thank you very much kiryl23

#33

According to the technical document ‘S2-PDGS-MPC-ATBD-1’, the computation of rTOA follows the equation:
image
where there is an absolute calibration coefficient named A_k;NTDI.
But it seems to be no such constant parameter in your converting equation ?

Would you please tell the reason ?
If this coefficient equals to 1, it should be omitted other than given out in the workflow equation.

Thank you so much !

#34

Or could you please explain why this constant coefficient is omitted when rTOA is converted to LTOA ?
Thank you so much ~ !

#35

Hello!
Unfortunately, certain pages on the network are no longer available.
But there is code and an explanation in the comments.
For example, you can see the calculation available here:

Method ReflectanceToRadianceOp.computeTile() in module:
https://github.com/senbox-org/s2tbx/blob/master/s2tbx-reflectance-to-radiance/src/main/java/org/esa/s2tbx/reflectance2radiance/ReflectanceToRadianceOp.java

or here:

Method L2A_Tables.TOA_refl2rad() in module:
https://github.com/JanBojko/SEN2COR/blob/master/sen2cor/L2A_Tables.py

Everything is explained in the comments, how to calculate TOA Radiance.
Аctual reference to the algorithm (29.12.2023): https://sentinels.copernicus.eu/web/sentinel/technical-guides/sentinel-2-msi/level-1c/algorithm-overview

S2_TOARadiance
S2_TOARadiance1571×732 111 KB