diff --git a/docs/workshop_latest.rst b/docs/workshop_latest.rst index cf17b7dd..2689fd5a 100644 --- a/docs/workshop_latest.rst +++ b/docs/workshop_latest.rst @@ -547,9 +547,9 @@ L2 satellite data already included in the package. >>> from wavy.satellite_module import satellite_class as sc >>> path_to_files= '/home/patrikb/wavy/tests/data/CCIv3_20Hz/' - >>> sd = '2019-3-24 18' - >>> ed = '2019-3-24 20' - >>> region = 'NordicSeas' + >>> sd = '2019-3-24 10' + >>> ed = '2019-3-24 11' + >>> region = 'Sulafj' >>> sco = sc(sd=sd, ed=ed, >>> nID='L2_20Hz_s3a', name='s3a', >>> region=region) @@ -559,16 +559,12 @@ You can have a first look at the data: .. code-block:: python3 - >>> sco.quicklook(a=True) + >>> sco.quicklook(a=True, land_mask_resolution='f') .. image:: ./oslo_ws24_filters_1.png :scale: 80 - -There could be some land interactions with the satellite track. It is possible to -filter it using *sco.filter_landMask()* function. Additionnally points closer than -a certain distance to the coast can be defined using *sco.filter_distance_to_coast()* -with *llim* argument. Finally, values under threshold (resp. over) can also be discarder using the -*sco.apply_limits()* method with *llim* argument (resp. *ulim*). + +There could be some land interactions with the satellite track. It is possible to filter it using *sco.filter_landMask()* function. Additionnally points closer than a certain distance to the coast can be defined using *sco.filter_distance_to_coast()* with *llim* argument. If specifics (range gate resolution and height over ground) about the satellite intsrument are known one can compute a possible foot print land interaction based on the size of the pulse limited footprint size byt adding *.filter_footprint_land_interaction()*. Finally, values under threshold (resp. over) can also be discarder using the *sco.apply_limits()* method with *llim* argument (resp. *ulim*). The corresponding code is the following: @@ -576,7 +572,8 @@ The corresponding code is the following: >>> sco_filter = sco.filter_landMask()\ ... .filter_distance_to_coast(llim=650)\ - ... .apply_limits(llim=0.5) + ... .filter_footprint_land_interaction()\ + ... .apply_limits(llim=0.1) .. image:: ./oslo_ws24_filters_2.png @@ -586,20 +583,23 @@ Some additional despiking method can be applied with *sco.despike_blockStd()*. .. code-block:: python3 - >>> sco = sco.despike_blockStd(slider=20, sigma=2, chunk_min=5) + >>> sco = sco.despike_blockStd(slider=20, sigma=2, chunk_min=5, sampling_rate_Hz=20) .. image:: ./oslo_ws24_filters_3.png :scale: 80 + +In order to show the computed x-track pulse-limited footprint size one can add: + +.. code-block:: python3 + + >>> sco = sco.despike_blockStd(slider=20, sigma=2, chunk_min=5, sampling_rate_Hz=20) -Another operation you can perform is smoothing the time serie, using a running mean -with the *sco.filter_runmean* method: +Another operation you can perform is smoothing the time serie, using a running mean with the *sco.filter_runmean* method: .. code-block:: python3 - >>> sco = sco.filter_runmean(window=11, chunk_min=5, sampling_rate_Hz=20) + >>> sco = sco.filter_runmean(window=5, chunk_min=5, sampling_rate_Hz=20) -.. image:: ./oslo_ws24_filters_4.png - :scale: 80 10. Saving data to netcdf ######################### diff --git a/wavy/quicklookmod.py b/wavy/quicklookmod.py index a369a0c6..d36037fb 100755 --- a/wavy/quicklookmod.py +++ b/wavy/quicklookmod.py @@ -75,7 +75,8 @@ def quicklook(self, a=False, projection=None, **kwargs): vmax = kwargs.get('vmax', np.nanmax(plot_var)) levels = kwargs.get('levels', - np.arange(vmin, vmax, .5)) + np.arange(0, vmax+.5, .5)) + #np.arange(vmin, vmax, .5)) cflevels = kwargs.get('cflevels', levels) clevels = kwargs.get('clevels', levels)