Hello,
I am processing L1B imagery to L2 using an atmospheric correction workflow based on Py6S (for optical data). However, after the correction, I observe color artifacts in the output images, as seen in the attached example. Also, the input data had the same artifacts.
My pipeline includes:
- Extracting solar irradiance and atmospheric parameters (water vapor, ozone, AOD) from NetCDF data
TOA to Radiance Conversion using the metadata
SixS Radiative Transfer Model for surface reflectance calculation
Manual assignment of CRS and transform if missing
Potential issues I suspect:
- Incorrect aerosol profile (currently using
AeroProfile.Maritime) - Geometric misalignment in resampling
- Incorrect atmospheric input values
Has anyone faced similar issues with Py6S-based correction? Any suggestions on debugging the reflectance artifacts would be highly appreciated.
Thanks in advance!