cloudnetpy 1.59.2__py3-none-any.whl → 1.59.4__py3-none-any.whl

cloudnetpy/plotting/plot_meta.py

@@ -18,6 +18,9 @@ class PlotMeta(NamedTuple):
         log_scale: Whether to plot data values in a logarithmic scale.
         moving_average: Whether to plot a moving average in a 1d plot.
         contour: Whether to plot contours on top of a filled colormap.
+        zero_line: Whether to plot a zero line in a 1d plot.
+        time_smoothing_duration: The duration of the time smoothing window
+            (in 2d plots) in minutes.
     """
 
     cmap: str = "viridis"
@@ -26,6 +29,8 @@ class PlotMeta(NamedTuple):
     log_scale: bool = False
     moving_average: bool = True
     contour: bool = False
+    zero_line: bool = False
+    time_smoothing_duration: int = 0
 
 
 _COLORS = {
@@ -122,6 +127,9 @@ _CLABEL = {
 }
 
 
+_MWR_SINGLE_SMOOTHING = 10
+_MWR_MULTI_SMOOTHING = 30
+
 ATTRIBUTES = {
     "rain-radar": {
         "rainfall_rate": PlotMeta(
@@ -134,25 +142,30 @@ ATTRIBUTES = {
             cmap="coolwarm",
             plot_range=(223.15, 323.15),
             contour=True,
+            time_smoothing_duration=_MWR_SINGLE_SMOOTHING,
         ),
         "potential_temperature": PlotMeta(
             cmap="coolwarm",
             plot_range=(260, 320),
             contour=True,
+            time_smoothing_duration=_MWR_SINGLE_SMOOTHING,
         ),
         "equivalent_potential_temperature": PlotMeta(
             cmap="coolwarm",
             plot_range=(260, 320),
             contour=True,
+            time_smoothing_duration=_MWR_SINGLE_SMOOTHING,
         ),
         "relative_humidity": PlotMeta(
             plot_range=(0, 120),
             contour=True,
+            time_smoothing_duration=_MWR_SINGLE_SMOOTHING,
         ),
         "absolute_humidity": PlotMeta(
             plot_range=(1e-4, 1e-2),
             log_scale=True,
             contour=True,
+            time_smoothing_duration=_MWR_SINGLE_SMOOTHING,
         ),
     },
     "mwr-multi": {
@@ -160,18 +173,24 @@ ATTRIBUTES = {
             cmap="coolwarm",
             plot_range=(223.15, 323.15),
             contour=True,
+            time_smoothing_duration=_MWR_MULTI_SMOOTHING,
         ),
         "potential_temperature": PlotMeta(
             cmap="coolwarm",
             plot_range=(260, 320),
             contour=True,
+            time_smoothing_duration=_MWR_MULTI_SMOOTHING,
         ),
         "equivalent_potential_temperature": PlotMeta(
-            cmap="coolwarm", plot_range=(260, 320), contour=True
+            cmap="coolwarm",
+            plot_range=(260, 320),
+            contour=True,
+            time_smoothing_duration=_MWR_MULTI_SMOOTHING,
         ),
         "relative_humidity": PlotMeta(
             plot_range=(0, 120),
             contour=True,
+            time_smoothing_duration=_MWR_MULTI_SMOOTHING,
         ),
     },
     "fallback": {
@@ -521,5 +540,8 @@ ATTRIBUTES = {
         "pia": PlotMeta(
             plot_range=(0, 3),
         ),
+        "lwp": PlotMeta(
+            zero_line=True,
+        ),
     },
 }

cloudnetpy/plotting/plotting.py

@@ -17,8 +17,10 @@ from matplotlib.ticker import AutoMinorLocator
 from matplotlib.transforms import Affine2D, Bbox
 from mpl_toolkits.axes_grid1 import make_axes_locatable
 from numpy import ma, ndarray
+from scipy.ndimage import uniform_filter
 
 from cloudnetpy.exceptions import PlottingError
+from cloudnetpy.instruments.ceilometer import calc_sigma_units
 from cloudnetpy.plotting.plot_meta import ATTRIBUTES, PlotMeta
 
 
@@ -480,6 +482,14 @@ class Plot2D(Plot):
 
         alt = self._screen_data_by_max_y(figure_data)
 
+        if (duration := self._plot_meta.time_smoothing_duration) > 0:
+            sigma_units = calc_sigma_units(
+                figure_data.time, alt * 1e3, sigma_minutes=duration, sigma_metres=0
+            )
+            valid_time_ind = ~np.all(self._data.mask, axis=1)
+            smoothed_data = uniform_filter(self._data[valid_time_ind, :], sigma_units)
+            self._data[valid_time_ind, :] = smoothed_data
+
         image = self._ax.pcolorfast(
             figure_data.time_including_gaps,
             alt,
@@ -539,6 +549,8 @@ class Plot1D(Plot):
         )
         if self._plot_meta.moving_average:
             self._plot_moving_average(figure_data)
+        if self._plot_meta.zero_line:
+            self._ax.axhline(0, color="black", alpha=0.5)
         self._fill_between_data_gaps(figure_data)
         self.sub_plot.set_yax(ylabel=units, y_limits=self._get_y_limits())
         pos = self._ax.get_position()

cloudnetpy/version.py

@@ -1,4 +1,4 @@
 MAJOR = 1
 MINOR = 59
-PATCH = 2
+PATCH = 4
 __version__ = f"{MAJOR}.{MINOR}.{PATCH}"

{cloudnetpy-1.59.2.dist-info → cloudnetpy-1.59.4.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: cloudnetpy
-Version: 1.59.2
+Version: 1.59.4
 Summary: Python package for Cloudnet processing
 Author: Simo Tukiainen
 License: MIT License

{cloudnetpy-1.59.2.dist-info → cloudnetpy-1.59.4.dist-info}/RECORD

@@ -8,7 +8,7 @@ cloudnetpy/metadata.py,sha256=v_VDo2vbdTxB0zIsfP69IcrwSKiRlLpsGdq6JPI4CoA,5306
 cloudnetpy/output.py,sha256=WoVTNuxni0DUr163vZ-_mDr1brXhY15XSlGMrq9Aoqw,14700
 cloudnetpy/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 cloudnetpy/utils.py,sha256=0TlHm71YtSrKXBsRKctitnhQrvZPE-ulEVeAQW-oK58,27398
-cloudnetpy/version.py,sha256=3kvolgwQsNaJOi84o81Nihx0L61gFF-c1CaOq-zOKT8,72
+cloudnetpy/version.py,sha256=Ek-ltBpir522jGt9XlS2hHrcIB47W2_qNJ8Umx09pAA,72
 cloudnetpy/categorize/__init__.py,sha256=gP5q3Vis1y9u9OWgA_idlbjfWXYN_S0IBSWdwBhL_uU,69
 cloudnetpy/categorize/atmos.py,sha256=fWW8ye_8HZASRAiYwURFKWzcGOYIA2RKeVxCq0lVOuM,12389
 cloudnetpy/categorize/atmos_utils.py,sha256=wndpwJxc2-QnNTkV8tc8I11Vs_WkNz9sVMX1fuGgUC4,3777
@@ -92,8 +92,8 @@ cloudnetpy/model_evaluation/tests/unit/test_plotting.py,sha256=h9V8JKmrO4v9bOvv-
 cloudnetpy/model_evaluation/tests/unit/test_statistical_methods.py,sha256=Ra3r4V0qbqkpDuaTYvEIbaasl0nZ5gmTLR4eGC0glBQ,9724
 cloudnetpy/model_evaluation/tests/unit/test_tools.py,sha256=Ia_VrLdV2NstX5gbx_3AZTOAlrgLAy_xFZ8fHYVX0xI,3817
 cloudnetpy/plotting/__init__.py,sha256=lg9Smn4BI0dVBgnDLC3JVJ4GmwoSnO-qoSd4ApvwV6Y,107
-cloudnetpy/plotting/plot_meta.py,sha256=HkG60l_UFOCy80mu5xoDlGJgeYUcWa1y_HfugaFWZWQ,15074
-cloudnetpy/plotting/plotting.py,sha256=j1OsTDlAhS78w9CUNfPVCH2vhlA3cl3cjC2wlCywqio,31788
+cloudnetpy/plotting/plot_meta.py,sha256=cLdCZrhbP-gaobS_zjcf8d2xVALzl7zh2qpttxCHyrg,15983
+cloudnetpy/plotting/plotting.py,sha256=LwTAHAdu3tFc-rnspTx3jfvNvmK54A-bRE8B3ty2-L8,32408
 cloudnetpy/products/__init__.py,sha256=2hRb5HG9hNrxH1if5laJkLeFeaZCd5W1q3hh4ewsX0E,273
 cloudnetpy/products/classification.py,sha256=0E9OUGR3uLCsS1nORwQu0SqW0_8uX7n6LlRcVhtzKw4,7845
 cloudnetpy/products/der.py,sha256=mam6jWV7A2h8V5WC3DIeFp6ou7UD1JOw9r7h2B0su-s,12403
@@ -107,8 +107,8 @@ cloudnetpy/products/mie_lu_tables.nc,sha256=It4fYpqJXlqOgL8jeZ-PxGzP08PMrELIDVe5
 cloudnetpy/products/mwr_tools.py,sha256=PRm5aCULccUehU-Byk55wYhhEHseMjoAjGBu5TSyHao,4621
 cloudnetpy/products/product_tools.py,sha256=rhx_Ru9FLlQqCNM-awoiHx18-Aq1eBwL9LiUaQoJs6k,10412
 docs/source/conf.py,sha256=IKiFWw6xhUd8NrCg0q7l596Ck1d61XWeVjIFHVSG9Og,1490
-cloudnetpy-1.59.2.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
-cloudnetpy-1.59.2.dist-info/METADATA,sha256=nHUt1jpzNXmkUrk6jbbNUtx-wCnDKT0tVbNnX25O09U,5784
-cloudnetpy-1.59.2.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
-cloudnetpy-1.59.2.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
-cloudnetpy-1.59.2.dist-info/RECORD,,
+cloudnetpy-1.59.4.dist-info/LICENSE,sha256=wcZF72bdaoG9XugpyE95Juo7lBQOwLuTKBOhhtANZMM,1094
+cloudnetpy-1.59.4.dist-info/METADATA,sha256=G2Zctni6oT3Nr7TbNTuMDMcTnCAvBdH9Qn5IWS32GtE,5784
+cloudnetpy-1.59.4.dist-info/WHEEL,sha256=oiQVh_5PnQM0E3gPdiz09WCNmwiHDMaGer_elqB3coM,92
+cloudnetpy-1.59.4.dist-info/top_level.txt,sha256=ibSPWRr6ojS1i11rtBFz2_gkIe68mggj7aeswYfaOo0,16
+cloudnetpy-1.59.4.dist-info/RECORD,,