Ok I finally got it. The problem was that due to a swapping of the order of lat/lon in gdal v3. A description of the issue is found here https://github.com/OSGeo/gdal/issues/1546 I am attaching an updated version of level1_landsat8.py which should fix the issue. Could you check ? Though it is not qu...

Hello,
It seems indeed that there are cases where lat and lon are swapped. I'm investigating and will keep you updated !
Cheers,
François

Dear Ding, Polymer is an atmospheric correction algorithm to retrieve water colour. It is not designed to work over land. With default processing parameters, the cloud mask triggers over most of the land pixels, and polymer returns bad results when it doesn't. An optional land mask (GSW) can be used...

Hello,
Please find my input and output files for this scene at the following address :
https://we.tl/t-JxBXtlrjxf
Cheers,
François

The command line interface polymer_cli.py is very basic and only supports default arguments.
You can not pass the argument resolution=20 unless you change it at the python level in polymer_cli.py by passing the argument resolution=20 to the object Level1.

A basic example of GSW usage is : from polymer.main import run_atm_corr from polymer.level2_nc import Level2_NETCDF from polymer.level1_msi import Level1_MSI from polymer.gsw import GSW run_atm_corr( Level1_MSI('<L1C product>', landmask=GSW('<gsw tiles storage directory>')) Level2_NETCDF(), ) The ti...

There is a quotation mark error, the argument must be provided outside of the string.

Hello, I tried to reproduce your error with HDF output but couldn't, in my case it didn't fail. Maybe if you could export your anaconda environment, if you are using anaconda, that could be helpful (conda env export > environment.yml). Concerning the error of geolocation in SNAP, this could be due t...

Dear ding, Polymer is an atmospheric correction algorithm to retrieve water colour. It is not designed to work over land. With default processing parameters, the cloud mask triggers over most of the land pixels, and polymer returns bad results when it doesn't. An optional land mask (GSW) can be used...

Dear GHesketh4, I tried the same file but did not get a file full of NaNs. Water reflectance is recovered over the glint area. Which version of Polymer did you use ? LC08_L1TP_011029_20200615_20200625_01_T1_polymer.jpg ding: please open a new thread if you have a precise question that is unrelated t...