Page 1 of 1
About the OLCI process - land masking
Posted: Fri Nov 02, 2018 11:46 am
by liugers11
Hi,
Thank you for providing the use of Polymer. I have one question in the use of Polymer for processing OLCI. What is the difference when landsat mask selecting "None" and "default"?
best wishes
Ge
Re: About the OLCI process
Posted: Mon Nov 05, 2018 8:53 am
by fsteinmetz
Dear Ge,
If you use
landmask=None, there is no landmask applied at all, Polymer will run also over land pixels.
With
landmask="default", it uses the landmask embedded in OLCI Level1, but this will mask out inland waters.
The last option,
landmask=GSW(), uses an external land mask (JRC global surface water - it will download it on the fly) and is recommended if you are interested in inland waters. For example, for OLCI:
Code: Select all
from polymer.gsw import GSW
run_atm_corr(
Level1_OLCI(..., landmask=GSW(agg=8)), # GSW native resolution is ~55m ; agg=8 aggregate to a resolution
Level2_NETCDF() # of about 440m to match OLCI full resolution
)
There is some documentation in the Python docstrings, in the python files themselves. We will try to improve the documentation in the near future.
Re: About the OLCI process
Posted: Sat Dec 22, 2018 11:44 am
by liugers11
Thank you for your answers. It really helps.