HYGEOS support forum

Hello,

I am checking bitmask, and although I checked different polymer-corrected NC files, my code doesn't show me any land pixels using land mask flag. I got cloud base pixels, clean water pixels, and case2 pixels, but no land pixels. Should I integrate the scene classification (SCL) masks from L2A Sentinel-2 products?

I appreciate any help you can provide.

Hello,

Yes, Polymer provides a landmask through the Global Surface Water dataset.
You can use it by providing a GSW instance to Level1_MSI through the 'landmask' argument :

from polymer.gsw import GSW
run_atm_corr(Level1_MSI(..., landmask=GSW(directory='/path/to/gsw_data/')), ...)

Make sure that the pointed directory exists and is empty, and the landmask data will be downloaded automatically.

Kind regards,
François