disdrodb 0.1.0__py3-none-any.whl → 0.1.1__py3-none-any.whl

disdrodb/utils/dataframe.py

@@ -0,0 +1,342 @@
+#!/usr/bin/env python3
+
+# -----------------------------------------------------------------------------.
+# Copyright (c) 2021-2023 DISDRODB developers
+#
+# This program is free software: you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation, either version 3 of the License, or
+# (at your option) any later version.
+#
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with this program.  If not, see <http://www.gnu.org/licenses/>.
+# -----------------------------------------------------------------------------.
+"""Dataframe utilities."""
+import numpy as np
+import pandas as pd
+
+
+def log_arange(start, stop, log_step=0.1, base=10):
+    """
+    Return numbers spaced evenly on a log scale (similar to np.arange but in log space).
+
+    Parameters
+    ----------
+    start : float
+        The starting value of the sequence (must be > 0).
+    stop : float
+        The end value of the sequence (must be > 0).
+    log_step : float
+        The step size in log-space (default is 0.1).
+    base : float
+        The logarithmic base (default is 10).
+
+    Returns
+    -------
+    np.ndarray
+        Array of values spaced in log scale.
+    """
+    if start <= 0 or stop <= 0:
+        raise ValueError("Both start and stop must be > 0 for log spacing.")
+
+    log_start = np.log(start) / np.log(base)
+    log_stop = np.log(stop) / np.log(base)
+
+    log_values = np.arange(log_start, log_stop, log_step)
+    return base**log_values
+
+
+def compute_1d_histogram(df, column, variables=None, bins=10, labels=None, prefix_name=True, include_quantiles=False):
+    """Compute conditional univariate statistics.
+
+    Parameters
+    ----------
+    df : pandas.DataFrame
+        Input dataframe
+    column : str
+        Column name to be binned.
+    variables : str or list, optional
+        Column names for which conditional statistics will be computed.
+        If None, only counts are computed.
+    bins : int or array-like
+        Number of bins or bin edges.
+    labels : array-like, optional
+        Labels for the column bins. If None, uses bin centers.
+
+    Returns
+    -------
+    pandas.DataFrame
+    """
+    # Copy data
+    df = df.copy()
+
+    # Ensure `variables` is a list of variables
+    # - If no variable specified, create dummy variable
+    if variables is None:
+        variables = ["dummy"]
+        df["dummy"] = np.ones(df[column].shape)
+        variables_specified = False
+    elif isinstance(variables, str):
+        variables = [variables]
+        variables_specified = True
+    elif isinstance(variables, list):
+        variables_specified = True
+    else:
+        raise TypeError("`variables` must be a string, list of strings, or None.")
+    variables = np.unique(variables)
+
+    # Handle column binning
+    if isinstance(bins, int):
+        bins = np.linspace(df[column].min(), df[column].max(), bins + 1)
+
+    # Drop rows where any of the key columns have NaN
+    df = df.dropna(subset=[column, *variables])
+
+    if len(df) == 0:
+        raise ValueError("No valid data points after removing NaN values")
+
+    # Create binned columns with explicit handling of out-of-bounds values
+    df[f"{column}_binned"] = pd.cut(df[column], bins=bins, include_lowest=True)
+
+    # Create complete IntervalIndex for both dimensions
+    intervals = df[f"{column}_binned"].cat.categories
+
+    # Create IntervalIndex with all possible combinations
+    full_index = pd.Index(intervals, name=f"{column}_binned")
+
+    # Define grouping object
+    df_grouped = df.groupby([f"{column}_binned"], observed=False)
+
+    # Compute statistics for specified variables
+    variables_stats = []
+    for i, var in enumerate(variables):
+        # Prepare prefix
+        prefix = f"{var}_" if prefix_name and variables_specified else ""
+
+        # Define statistics to compute
+        if variables_specified:
+            # Compute quantiles
+            quantiles = [0.01, 0.05, 0.10, 0.25, 0.50, 0.75, 0.90, 0.95, 0.99]
+            df_stats_quantiles = df_grouped[var].quantile(quantiles).unstack(level=-1)
+            df_stats_quantiles.columns = [f"{prefix}Q{int(q*100)}" for q in df_stats_quantiles.columns]
+            df_stats_quantiles = df_stats_quantiles.rename(
+                columns={
+                    f"{prefix}Q50": f"{prefix}median",
+                },
+            )
+            # Define other stats to compute
+            list_stats = [
+                (f"{prefix}std", "std"),
+                (f"{prefix}min", "min"),
+                (f"{prefix}max", "max"),
+                (f"{prefix}mad", lambda s: np.median(np.abs(s - np.median(s)))),
+            ]
+            if i == 0:
+                list_stats.append(("count", "count"))
+        else:
+            list_stats = [("count", "count")]
+
+        # Compute statistics
+        df_stats = df_grouped[var].agg(list_stats)
+
+        # Compute other variable statistics
+        if variables_specified:
+            df_stats[f"{prefix}range"] = df_stats[f"{prefix}max"] - df_stats[f"{prefix}min"]
+            df_stats[f"{prefix}iqr"] = df_stats_quantiles[f"{prefix}Q75"] - df_stats_quantiles[f"{prefix}Q25"]
+            df_stats[f"{prefix}ipr80"] = df_stats_quantiles[f"{prefix}Q90"] - df_stats_quantiles[f"{prefix}Q10"]
+            df_stats[f"{prefix}ipr90"] = df_stats_quantiles[f"{prefix}Q95"] - df_stats_quantiles[f"{prefix}Q5"]
+            df_stats[f"{prefix}ipr98"] = df_stats_quantiles[f"{prefix}Q99"] - df_stats_quantiles[f"{prefix}Q1"]
+            if include_quantiles:
+                df_stats = pd.concat((df_stats, df_stats_quantiles), axis=1)
+            else:
+                df_stats[f"{prefix}median"] = df_stats_quantiles[f"{prefix}median"]
+        variables_stats.append(df_stats)
+
+    # Combine all statistics into a single DataFrame
+    df_stats = pd.concat(variables_stats, axis=1)
+
+    # Reindex to include all interval combinations
+    df_stats = df_stats.reindex(full_index)
+
+    # Determine bin centers
+    centers = intervals.mid
+
+    # Use provided labels if available
+    coords = labels if labels is not None else centers
+
+    # Reset index and add coordinates/labels
+    df_stats = df_stats.reset_index()
+    df_stats[f"{column}"] = pd.Categorical(df_stats[f"{column}_binned"].map(dict(zip(intervals, coords, strict=False))))
+    df_stats = df_stats.drop(columns=f"{column}_binned")
+
+    return df_stats
+
+
+def compute_2d_histogram(
+    df,
+    x,
+    y,
+    variables=None,
+    x_bins=10,
+    y_bins=10,
+    x_labels=None,
+    y_labels=None,
+    prefix_name=True,
+    include_quantiles=False,
+):
+    """Compute conditional bivariate statistics.
+
+    Parameters
+    ----------
+    df : pandas.DataFrame
+        Input dataframe
+    x : str
+        Column name for x-axis binning (will be rounded to integers)
+    y : str
+        Column name for y-axis binning
+    variables : str or list, optional
+        Column names for which statistics will be computed.
+        If None, only counts are computed.
+    x_bins : int or array-like
+        Number of bins or bin edges for x
+    y_bins : int or array-like
+        Number of bins or bin edges for y
+    x_labels : array-like, optional
+        Labels for x bins. If None, uses bin centers
+    y_labels : array-like, optional
+        Labels for y bins. If None, uses bin centers
+
+    Returns
+    -------
+    xarray.Dataset
+        Dataset with dimensions corresponding to binned variables and
+        data variables for each statistic
+    """
+    # # If polars, cast to pandas
+    # if isinstance(df, pl.DataFrame):
+    #     df = df.to_pandas()
+
+    # Copy data
+    df = df.copy()
+
+    # Ensure `variables` is a list of variables
+    # - If no variable specified, create dummy variable
+    if variables is None:
+        variables = ["dummy"]
+        df["dummy"] = np.ones(df[x].shape)
+        variables_specified = False
+    elif isinstance(variables, str):
+        variables = [variables]
+        variables_specified = True
+    elif isinstance(variables, list):
+        variables_specified = True
+    else:
+        raise TypeError("`variables` must be a string, list of strings, or None.")
+    variables = np.unique(variables)
+
+    # Handle x-axis binning
+    if isinstance(x_bins, int):
+        x_bins = np.linspace(df[x].min(), df[x].max(), x_bins + 1)
+    # Handle y-axis binning
+    if isinstance(y_bins, int):
+        y_bins = np.linspace(df[y].min(), df[y].max(), y_bins + 1)
+
+    # Drop rows where any of the key columns have NaN
+    df = df.dropna(subset=[x, y, *variables])
+
+    if len(df) == 0:
+        raise ValueError("No valid data points after removing NaN values")
+
+    # Create binned columns with explicit handling of out-of-bounds values
+    df[f"{x}_binned"] = pd.cut(df[x], bins=x_bins, include_lowest=True)
+    df[f"{y}_binned"] = pd.cut(df[y], bins=y_bins, include_lowest=True)
+
+    # Create complete IntervalIndex for both dimensions
+    x_intervals = df[f"{x}_binned"].cat.categories
+    y_intervals = df[f"{y}_binned"].cat.categories
+
+    # Create MultiIndex with all possible combinations
+    full_index = pd.MultiIndex.from_product([x_intervals, y_intervals], names=[f"{x}_binned", f"{y}_binned"])
+
+    # Define grouping object
+    df_grouped = df.groupby([f"{x}_binned", f"{y}_binned"], observed=False)
+
+    # Compute statistics for specified variables
+    variables_stats = []
+    for i, var in enumerate(variables):
+        # Prepare prefix
+        prefix = f"{var}_" if prefix_name and variables_specified else ""
+
+        # Define statistics to compute
+        if variables_specified:
+            # Compute quantiles
+            quantiles = [0.01, 0.05, 0.10, 0.25, 0.50, 0.75, 0.90, 0.95, 0.99]
+            df_stats_quantiles = df_grouped[var].quantile(quantiles).unstack(level=-1)
+            df_stats_quantiles.columns = [f"{prefix}Q{int(q*100)}" for q in df_stats_quantiles.columns]
+            df_stats_quantiles = df_stats_quantiles.rename(
+                columns={
+                    f"{prefix}Q50": f"{prefix}median",
+                },
+            )
+            # Define other stats to compute
+            list_stats = [
+                (f"{prefix}std", "std"),
+                (f"{prefix}min", "min"),
+                (f"{prefix}max", "max"),
+                (f"{prefix}mad", lambda s: np.median(np.abs(s - np.median(s)))),
+            ]
+            if i == 0:
+                list_stats.append(("count", "count"))
+        else:
+            list_stats = [("count", "count")]
+
+        # Compute statistics
+        df_stats = df_grouped[var].agg(list_stats)
+
+        # Compute other variable statistics
+        if variables_specified:
+            df_stats[f"{prefix}range"] = df_stats[f"{prefix}max"] - df_stats[f"{prefix}min"]
+            df_stats[f"{prefix}iqr"] = df_stats_quantiles[f"{prefix}Q75"] - df_stats_quantiles[f"{prefix}Q25"]
+            df_stats[f"{prefix}ipr80"] = df_stats_quantiles[f"{prefix}Q90"] - df_stats_quantiles[f"{prefix}Q10"]
+            df_stats[f"{prefix}ipr90"] = df_stats_quantiles[f"{prefix}Q95"] - df_stats_quantiles[f"{prefix}Q5"]
+            df_stats[f"{prefix}ipr98"] = df_stats_quantiles[f"{prefix}Q99"] - df_stats_quantiles[f"{prefix}Q1"]
+            if include_quantiles:
+                df_stats = pd.concat((df_stats, df_stats_quantiles), axis=1)
+            else:
+                df_stats[f"{prefix}median"] = df_stats_quantiles[f"{prefix}median"]
+        variables_stats.append(df_stats)
+
+    # Combine all statistics into a single DataFrame
+    df_stats = pd.concat(variables_stats, axis=1)
+
+    # Reindex to include all interval combinations
+    df_stats = df_stats.reindex(full_index)
+
+    # Determine coordinates
+    x_centers = x_intervals.mid
+    y_centers = y_intervals.mid
+
+    # Use provided labels if available
+    x_coords = x_labels if x_labels is not None else x_centers
+    y_coords = y_labels if y_labels is not None else y_centers
+
+    # Reset index and set new coordinates
+    df_stats = df_stats.reset_index()
+    df_stats[f"{x}"] = pd.Categorical(df_stats[f"{x}_binned"].map(dict(zip(x_intervals, x_coords, strict=False))))
+    df_stats[f"{y}"] = pd.Categorical(df_stats[f"{y}_binned"].map(dict(zip(y_intervals, y_coords, strict=False))))
+
+    # Set new MultiIndex with coordinates
+    df_stats = df_stats.set_index([f"{x}", f"{y}"])
+    df_stats = df_stats.drop(columns=[f"{x}_binned", f"{y}_binned"])
+
+    # Convert to dataset
+    ds = df_stats.to_xarray()
+
+    # Transpose arrays
+    ds = ds.transpose(y, x)
+    return ds

{disdrodb-0.1.0.dist-info → disdrodb-0.1.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: disdrodb
-Version: 0.1.0
+Version: 0.1.1
 Summary: disdrodb provides tools to download, standardize, share and analyze global disdrometer data.
 Author: Gionata Ghiggi
 Project-URL: homepage, https://github.com/ltelab/disdrodb
@@ -47,7 +47,7 @@ Requires-Dist: twine; extra == "dev"
 Requires-Dist: loghub; extra == "dev"
 Dynamic: license-file
 
-# 📦 DISDRODB - A package to standardize, process and analyze global disdrometer data.
+# 📦 DISDRODB - A package to standardize, process and analyze global disdrometer data
 
 |                   |                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   |
 | ----------------- | --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- |
@@ -60,17 +60,16 @@ Dynamic: license-file
 | Linting           | [![Black](https://img.shields.io/badge/code%20style-black-000000.svg?style=flat)](https://github.com/psf/black) [![Ruff](https://img.shields.io/endpoint?url=https://raw.githubusercontent.com/astral-sh/ruff/main/assets/badge/v2.json&style=flat)](https://github.com/astral-sh/ruff) [![Codespell](https://img.shields.io/badge/Codespell-enabled-brightgreen?style=flat)](https://github.com/codespell-project/codespell)                                                                                                                                                                                                                                                     |
 | Code Coverage     | [![Coveralls](https://coveralls.io/repos/github/ltelab/disdrodb/badge.svg?branch=main&style=flat)](https://coveralls.io/github/ltelab/disdrodb?branch=main) [![Codecov](https://codecov.io/gh/ltelab/disdrodb/branch/main/graph/badge.svg?style=flat)](https://codecov.io/gh/ltelab/disdrodb)                                                                                                                                                                                                                                                                                                                                                                                     |
 | Code Quality      | [![Codefactor](https://www.codefactor.io/repository/github/ltelab/disdrodb/badge?style=flat)](https://www.codefactor.io/repository/github/ltelab/disdrodb) [![Codebeat](https://codebeat.co/badges/14ff831b-f064-4bdd-a2e2-72ffdf28a35a?style=flat)](https://codebeat.co/projects/github-com-ltelab-disdrodb-main) [![Codacy](https://app.codacy.com/project/badge/Grade/d823c50a7ad14268bd347b5aba384623?style=flat)](https://app.codacy.com/gh/ltelab/disdrodb/dashboard?utm_source=gh&utm_medium=referral&utm_content=&utm_campaign=Badge_grade) [![Codescene](https://codescene.io/projects/36773/status-badges/code-health?style=flat)](https://codescene.io/projects/36773) |
-| License           | [![License](https://img.shields.io/github/license/ltelab/disdrodb?style=flat)](https://github.com/ltelab/disdrodb/blob/main/LICENSE) |
-| Community         | [![Slack](https://img.shields.io/badge/Slack-disdrodb-green.svg?logo=slack&style=flat)](https://join.slack.com/t/disdrodbworkspace/shared_invite/zt-25l4mvgo7-cfBdXalzlWGd4Pt7H~FqoA) [![GitHub Discussions](https://img.shields.io/badge/GitHub-Discussions-green?logo=github&style=flat)](https://github.com/ltelab/disdrodb/discussions)  |
-| Citation          | [![DOI](https://zenodo.org/badge/429018433.svg?style=flat)](https://zenodo.org/doi/10.5281/zenodo.7680581) |
+| License           | [![License](https://img.shields.io/github/license/ltelab/disdrodb?style=flat)](https://github.com/ltelab/disdrodb/blob/main/LICENSE)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              |
+| Community         | [![Slack](https://img.shields.io/badge/Slack-disdrodb-green.svg?logo=slack&style=flat)](https://join.slack.com/t/disdrodbworkspace/shared_invite/zt-25l4mvgo7-cfBdXalzlWGd4Pt7H~FqoA) [![GitHub Discussions](https://img.shields.io/badge/GitHub-Discussions-green?logo=github&style=flat)](https://github.com/ltelab/disdrodb/discussions)                                                                                                                                                                                                                                                                                                                                       |
+| Citation          | [![DOI](https://zenodo.org/badge/429018433.svg?style=flat)](https://zenodo.org/doi/10.5281/zenodo.7680581)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        |
 
 [**Slack**](https://join.slack.com/t/disdrodbworkspace/shared_invite/zt-25l4mvgo7-cfBdXalzlWGd4Pt7H~FqoA) | [**Documentation**](https://disdrodb.readthedocs.io/en/latest/)
 
 DISDRODB is part of an international joint effort to index, collect and homogenize drop size distribution (DSD) data from around the world.
 
 The DISDRODB project also aims to establish a global standard for sharing disdrometer observations.
-Built on FAIR data principles and Climate & Forecast (CF) conventions, DISDRODB standards facilitates the processing, analysis, and visualization of disdrometer data.
-
+Built on FAIR data principles and Climate & Forecast (CF) conventions, DISDRODB standards facilitate the processing, analysis and visualization of disdrometer data.
 
 ## ℹ️ Software Overview
 
@@ -87,13 +86,13 @@ The software enables you to:
 - Compute empirical and model-based drop size distribution parameters and derive geophysical and polarimetric radar variables of interest (DISDRODB L2 product)
 
 Currently, the DISDRODB Working Group is finalizing the development of the L1 and L2 scientific products.
-If you have ideas, algorithms, data or expertise to share, or you want to contribute to the future DISDRODB products, do not hesitate to get in touch !!!
+If you have ideas, algorithms, data, or expertise to share, or you want to contribute to the future DISDRODB products, do not hesitate to get in touch!!!
 
-Join the [**DISDRODB Slack Workspace**](https://join.slack.com/t/disdrodbworkspace/shared_invite/zt-25l4mvgo7-cfBdXalzlWGd4Pt7H~FqoA) to meet the DISDRODB Community !
+Join the [**DISDRODB Slack Workspace**](https://join.slack.com/t/disdrodbworkspace/shared_invite/zt-25l4mvgo7-cfBdXalzlWGd4Pt7H~FqoA) to meet the DISDRODB Community!
 
 ## 🚀 Quick Start
 
-You're about to create your very own DISDRODB Local Data Archive.
+Create your own DISDRODB Local Data Archive.
 
 ### 📚 Download the DISDRODB Metadata Archive
 
@@ -102,76 +101,62 @@ The DISDRODB Metadata Archive is a collection of metadata files that describe th
 To download the DISDRODB Metadata Archive, navigate to the desired directory and run:
 
 ```bash
-
 git clone https://github.com/ltelab/DISDRODB-METADATA.git
-
 ```
 
-This command will download a directory called `DISDRODB-METADATA`.
-
-Alternatively, you can download a static version of the DISDRODB Metadata Archive using:
+Or download a static snapshot without using git:
 
 ```bash
-disdrodb_download_metadata_archive /custom/path/where/to/download/DISDRODB-METADATA
+disdrodb_download_metadata_archive /path/to/DISDRODB-METADATA
 ```
 
 ### 📚 Define the DISDRODB Configuration File
 
-The disdrodb software needs to know where the local DISDRODB Metadata Archive
-is stored on your local machine, as well as where you want to download the raw stations data
-as well as where to save the DISDRODB products you will generate.
-
-Within the `disdrodb` package, the argument `metadata_archive_dir` refers to the base directory of your local DISDRODB Metadata Archive, while `data_archive_dir` refers to the base directory of your local DISDRODB Data Archive.
+The disdrodb software requires to know two directories:
 
-The `metadata_archive_dir` path corresponds to the `DISDRODB` directory within the `DISDRODB-METADATA` archive.
-The `data_archive_dir` path corresponds to `DISDRODB` directory of choice where
-all DISDRODB products will be saved.
+- `metadata_archive_dir`: the base of your local DISDRODB Metadata Archive
+- `data_archive_dir`: the base of your local DISDRODB Data Archive
 
-To facilitate the creation of the DISDRODB Configuration File, you can adapt and run in python the following code snippet.
-Please note that on Windows these paths must end with `\DISDRODB`, while on macOS/Linux they must end with `/DISDRODB`.
+On Windows, paths must end with `\DISDRODB`; on macOS/Linux, they must end with `/DISDRODB`.
 
 ```python
 import disdrodb
 
 metadata_archive_dir = "/<path_to>/DISDRODB-METADATA/DISDRODB"
-data_archive_dir = "/<path_of_choice_to_the_local_data_archive>/DISDRODB"
+data_archive_dir = "/<path_to>/DISDRODB"
 disdrodb.define_configs(
     metadata_archive_dir=metadata_archive_dir, data_archive_dir=data_archive_dir
 )
 ```
 
-By running this command, the disdrodb software will write a `.config_disdrodb.yml` file into your home directory (i.e. `~/.config_disdrodb.yml`)
-that will be used as default configuration file when running the disdrodb software.
+This creates a `.config_disdrodb.yml` file in your home directory (e.g., `~/.config_disdrodb.yml`).
 
-If you **now close your python session and reopen a new one**, if you will run the following code snippet, you
-should get the `metadata_archive_dir` and `data_archive_dir` paths you just defined in the DISDRODB Configuration File:
+To verify the configuration, open a new Python session and run:
 
 ```python
 import disdrodb
 
-print("DISDRODB Metadata Archive Directory: ", disdrodb.get_metadata_archive_dir())
-print("DISDRODB Data Archive Directory: ", disdrodb.get_data_archive_dir())
+print("Metadata Archive Directory:", disdrodb.get_metadata_archive_dir())
+print("Data Archive Directory:", disdrodb.get_data_archive_dir())
 ```
 
-You can also verify and print the default DISDRODB Metadata Archive and Data Archive directories
-by typing the following command in the terminal:
+Or in the shell:
 
 ```bash
-
-   disdrodb_data_archive_directory
-   disdrodb_metadata_archive_directory
-
+disdrodb_metadata_archive_directory
+disdrodb_data_archive_directory
 ```
 
 ### 📥 Download the DISDRODB Raw Data Archive
 
-To download all data stored into the DISDRODB Decentralized Data Archive, you just have to run the following command:
+To download all data stored into the DISDRODB Decentralized Data Archive,
+you just have to run the following command:
 
 ```bash
 disdrodb_download_archive
 ```
 
-If you aim to download data from a specific source (e.g., EPFL), run:
+To download from a specific source (e.g., EPFL):
 
 ```bash
 disdrodb_download_archive --data-sources EPFL
@@ -185,7 +170,7 @@ To open the local DISDRODB Data Archive directory, type:
 disdrodb_open_data_archive
 ```
 
-### 💫 Transform the raw data to standardized netCDF files
+### 💫 Transform Raw Data to Standardized netCDFs
 
 If you want to convert all stations raw data into standardized netCDF4 files, run the following command in the terminal:
 
@@ -195,40 +180,28 @@ disdrodb_run_l0
 
 Type `disdrodb_run_l0 --help` to see further options.
 
-### 💫 Generate the DISDRODB L1 and L2 products
+### 💫 Generate DISDRODB L1 and p2 products
 
-To generate the DISDRODB L1 and L2 products, run the following commands in the terminal:
+To generate DISDRODB L1 and L2 products, run the following commands in the terminal:
 
 ```bash
 disdrodb_run_l1
-```
-```bash
 disdrodb_run_l2e
-```
-```bash
 disdrodb_run_l2m
 ```
 
-### 💫 Open analysis-ready DISDRODB products
+### 💫 Analyze Analysis‐Ready Products
 
-The disdrodb software `open_dataset` function enable to lazy open all station files of
-a DISDRODB product into a `xarray.Dataset` (or `pandas.DataFrame` for the DISDRODB L0A product).
+The software’s `open_dataset` function **lazily** opens all station files of a given product:
 
 ```python
 import disdrodb
 
-# Define station arguments
-data_source = "EPFL"
-campaign_name = "HYMEX_LTE_SOP3"
-station_name = "10"
-
-# Open all station files of a specific product
 ds = disdrodb.open_dataset(
     product="L0C",
-    # Station arguments
-    data_source=data_source,
-    campaign_name=campaign_name,
-    station_name=station_name,
+    data_source="EPFL",
+    campaign_name="HYMEX_LTE_SOP3",
+    station_name="10",
 )
 ds
 ```