GSW landmask with MSI

[gbourdin]

Hi François, hi everyone,
I am trying to apply the GSW landmask to MSI scenes with the command below with polymer-4.17beta2 but I see no difference with the scene I process without using the GSW landmask. The GSW netcdf are downloaded in the correct directory and I get no error when the image is processed but I also noticed no output saying the GSW files are being downloaded. In fact, there is no output whatsoever while the image is being processed besides a warning that does not look related:

Code: Select all

site-packages/pyproj/crs/crs.py:141: FutureWarning: '+init=<authority>:<code>' syntax is deprecated. '<authority>:<code>' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6
  in_crs_string = _prepare_from_proj_string(in_crs_string)

I have tried with or without inputing a custom directory for GSW but nothing change and I still get land pixels processed.
Anyone had a chance to use the GSW mask on MSI with the last version of polymer (polymer-4.17beta2)?
Best,
Guillaume

Code: Select all

run_atm_corr(
    Level1_MSI('/home/S2B_MSIL1C_20190714T152819_N0208_R111_T22WDA_20190714T185950.SAFE/GRANULE/L1C_T22WDA_A012292_20190714T153316',
               resolution='20',
               ancillary=Ancillary_ERA5(directory=anc_dir),
               landmask=GSW(directory=gswdirw),
               sline=startline,
               eline=endline,
               scol=startcol,
               ecol=endcol),
    Level2_NETCDF(outdir='/home',
                  ext='S2B_MSIL2_20190714T152819_N0208_R111_T22WDA_20190714T185950.nc',
                  overwrite=force),
    multiprocessing=0
)