grdwindinversion 0.2.7__tar.gz → 0.3.2__tar.gz

{grdwindinversion-0.2.7 → grdwindinversion-0.3.2}/.github/workflows/build.yml

@@ -4,10 +4,12 @@ on:
   release:
     types:
       - created
+  push:
+    # Triggers on push to any branch
 
 jobs:
   build-and-upload-conda-package:
-    if: true 
+    if: contains(github.event.head_commit.message, '[force-build]') || startsWith(github.ref, 'refs/tags/')
     runs-on: ubuntu-latest
     #${{ github.event.workflow_run.conclusion == 'success' }}
     steps:
@@ -29,6 +31,16 @@ jobs:
           conda activate base
           conda install -y -c conda-forge conda-build anaconda-client
 
+      - name: Set version without 'v'
+        run: |
+          if [[ "${GITHUB_REF}" == refs/tags/* ]]; then
+            VERSION=${GITHUB_REF_NAME#v}
+          else
+            VERSION="0.0.0"  # Dummy version for testing
+          fi
+          echo "VERSION=${VERSION}" >> $GITHUB_ENV
+          echo "GIT_DESCRIBE_TAG=${VERSION}" >> $GITHUB_ENV
+
       - name: Build and Upload Conda package
         env:
           ANACONDA_TOKEN: ${{ secrets.ANACONDA_TOKEN }}

grdwindinversion-0.3.2/.github/workflows/ci.yml

@@ -0,0 +1,105 @@
+name: CI Workflow for grdwindinversion
+
+on:
+    push:
+        branches:
+            - main
+    pull_request:
+        branches:
+            - main
+
+jobs:
+    test:
+        runs-on: ubuntu-latest
+        timeout-minutes: 360 # 6 hours limit for the job
+
+        steps:
+            # Checkout the code
+            - name: Checkout code
+              uses: actions/checkout@v4
+
+            - uses: mamba-org/setup-micromamba@v1
+              with:
+                  micromamba-version: "1.5.9-1" # any version from https://github.com/mamba-org/micromamba-releases
+                  channels: tcevaer, conda-forge
+                  init-shell: bash
+                  post-cleanup: "all"
+
+            - name: Configure Conda channel priority to disabled
+              run: |
+                  conda config --set channel_priority disabled
+
+            - name: Create environment and install tools
+              run: micromamba create -n grdwind_env pytest conda-build boa python=3.10 -y
+
+            - name: Build package
+              run: |
+                  cd recipe
+                  eval "$(micromamba shell hook --shell bash)"
+                  micromamba activate grdwind_env
+                  conda mambabuild .
+
+            # Install the built package into the environment
+            - name: Install the built package
+              run: |
+                  eval "$(micromamba shell hook --shell bash)"
+                  micromamba activate grdwind_env
+                  conda install --use-local grdwindinversion -y
+
+              # Cache the test data if previously downloaded (up to 10 GB limit for the cache)
+              # WARNING : modify the key if the data is modified !!
+            - name: Cache test data
+              uses: actions/cache@v4
+              id: cache
+              with:
+                  path: ./test_data
+                  key: test-data-v3
+                  restore-keys: test-data-v3
+
+              # Download test data if not already cached
+            - name: Download test data
+              if: steps.cache.outputs.cache-hit != 'true' # Only download if cache miss
+              run: |
+                  mkdir -p ./test_data/
+                  wget https://cloud.ifremer.fr/index.php/s/ExLQ2TnYAqozPWE/download -O /tmp/ecmwf.zip
+                  unzip /tmp/ecmwf.zip -d ./test_data/
+                  wget https://cloud.ifremer.fr/index.php/s/kRgdOOPsjoZieZR/download -O /tmp/l1.zip
+                  unzip /tmp/l1.zip -d ./test_data/
+              timeout-minutes: 200 # Adjust depending on the size of your data
+
+              # Set up xsar configuration
+            - name: Setup xsar configuration
+              run: |
+                  mkdir -p ~/.xsar
+                  echo "data_dir: /tmp" > ~/.xsar/config.yaml
+                  echo "auxiliary_dir: ./test_data/auxiliary" >> ~/.xsar/config.yaml
+                  echo "path_dataframe_aux: ./test_data/auxiliary/active_aux.csv" >> ~/.xsar/config.yaml
+
+              # Set up grdwindinversion configuration
+            - name: Setup grdwindinversion configuration
+              run: |
+                  mkdir -p ~/.grdwindinversion
+                  echo "'ecmwf_0100_1h': ./test_data/ECMWF/forecast/hourly/0100deg/netcdf_light/%Y/%j/ECMWF_FORECAST_0100_%Y%m%d%H%M_10U_10V.nc" > ~/.grdwindinversion/data_config.yaml
+                  echo "'ecmwf_0125_1h': ./test_data/ECMWF/0.125deg/1h/forecasts/%Y/%j/ecmwf_%Y%m%d%H%M.nc" >> ~/.grdwindinversion/data_config.yaml
+                  echo "unit_test_s1_product: './test_data/L1/S1A_IW_GRDH_1SDV_20210909T130650_20210909T130715_039605_04AE83_C34F.SAFE'" >> ~/.grdwindinversion/data_config.yaml
+                  echo "unit_test_rcm_product: './test_data/L1/RCM1_OK2767220_PK2769320_1_SCLND_20230930_214014_VV_VH_GRD'" >> ~/.grdwindinversion/data_config.yaml
+                  echo "unit_test_rs2_product: './test_data/L1/RS2_OK141302_PK1242223_DK1208537_SCWA_20220904_093402_VV_VH_SGF'" >> ~/.grdwindinversion/data_config.yaml
+
+                  #echo "'nc_luts_path': ./test_data/GMFS/nc_luts" >> ~/.grdwindinversion/data_config.yaml
+                  #echo "'lut_cmod7_path': './test_data/GMFS/v1.6/GMF_cmod7_official/cmod7_and_python_script'" >> ~/.grdwindinversion/data_config.yaml
+                  #echo "'lut_ms1ahw_path': './test_data/GMFS/v1.6/GMF_cmodms1ahw'" >> ~/.grdwindinversion/data_config.yaml
+
+              # Run the tests
+            - name: Run tests
+              run: |
+                  eval "$(micromamba shell hook --shell bash)"
+                  micromamba activate grdwind_env
+                  pytest
+
+                  # Optionally, upload test artifacts (NetCDF files or logs) if needed
+                  #- name: Upload test artifacts
+                  #  if: failure()  # Only upload on failure
+                  #  uses: actions/upload-artifact@v2
+                  #  with:
+                  #    name: test-output
+                  #    path: ./test_output/

grdwindinversion-0.3.2/PKG-INFO

@@ -0,0 +1,67 @@
+Metadata-Version: 2.1
+Name: grdwindinversion
+Version: 0.3.2
+Summary: Package to perform Wind inversion from GRD Level-1 SAR images
+Author-email: Antoine Grouazel <antoine.grouazel@ifremer.fr>
+License: MIT
+Keywords: xarray,earth-observation,remote-sensing,satellite-imagery,Sentinel-1,RCM,RadarSat2,sar,synthetic-aperture-radar
+Classifier: Development Status :: 4 - Beta
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python
+Classifier: Programming Language :: Python :: 3 :: Only
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Intended Audience :: Science/Research
+Classifier: Topic :: Scientific/Engineering
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+License-File: AUTHORS.rst
+Requires-Dist: xsar
+Requires-Dist: xsarsea
+Requires-Dist: xarray
+Requires-Dist: xarray-datatree
+Requires-Dist: pyyaml
+Requires-Dist: numpy
+Requires-Dist: scipy
+Requires-Dist: fsspec
+Requires-Dist: aiohttp
+
+
+[![Python Version](https://img.shields.io/pypi/pyversions/grdwindinversion.svg)](https://pypi.org/project/grdwindinversion/)
+[![Dependencies Status](https://img.shields.io/badge/dependencies-up%20to%20date-brightgreen.svg)](https://github.com/umr-lops/grdwindinversion/pulls?utf8=%E2%9C%93&q=is%3Apr%20author%3Aapp%2Fdependabot)
+
+Package to perform Wind inversion from GRD Level-1 SAR images
+
+-   Free software: MIT license
+-   Documentation: https://grdwindinversion.readthedocs.io.
+
+## Usage
+
+```python
+
+   SAR_L1-to-L2_wind_processor -h
+   usage: SAR_L1-to-L2_wind_processor [-h] --input_file INPUT_FILE [--config_file CONFIG_FILE] --outputdir OUTPUTDIR [--verbose] [--overwrite]
+
+   Perform inversion from S1(L1-GRD) SAFE, L1-RCM, L1-RS2 ; using xsar/xsarsea tools
+
+   options:
+     -h, --help            show this help message and exit
+     --input_file INPUT_FILE
+                           input file path
+     --config_file CONFIG_FILE
+                           config file path [if not provided will take config file based on input file]
+     --outputdir OUTPUTDIR
+     --verbose
+     --overwrite           overwrite existing .nc files [default is False]
+```
+
+## Features
+
+This Python library (based on `xarray`) allows to perform wind inversion from level-1 GRD (projected magnitude image).
+Mission supported:
+
+-   Sentinel-1
+-   RCM
+-   RadarSat-2

grdwindinversion-0.3.2/README.md

@@ -0,0 +1,37 @@
+
+[![Python Version](https://img.shields.io/pypi/pyversions/grdwindinversion.svg)](https://pypi.org/project/grdwindinversion/)
+[![Dependencies Status](https://img.shields.io/badge/dependencies-up%20to%20date-brightgreen.svg)](https://github.com/umr-lops/grdwindinversion/pulls?utf8=%E2%9C%93&q=is%3Apr%20author%3Aapp%2Fdependabot)
+
+Package to perform Wind inversion from GRD Level-1 SAR images
+
+-   Free software: MIT license
+-   Documentation: https://grdwindinversion.readthedocs.io.
+
+## Usage
+
+```python
+
+   SAR_L1-to-L2_wind_processor -h
+   usage: SAR_L1-to-L2_wind_processor [-h] --input_file INPUT_FILE [--config_file CONFIG_FILE] --outputdir OUTPUTDIR [--verbose] [--overwrite]
+
+   Perform inversion from S1(L1-GRD) SAFE, L1-RCM, L1-RS2 ; using xsar/xsarsea tools
+
+   options:
+     -h, --help            show this help message and exit
+     --input_file INPUT_FILE
+                           input file path
+     --config_file CONFIG_FILE
+                           config file path [if not provided will take config file based on input file]
+     --outputdir OUTPUTDIR
+     --verbose
+     --overwrite           overwrite existing .nc files [default is False]
+```
+
+## Features
+
+This Python library (based on `xarray`) allows to perform wind inversion from level-1 GRD (projected magnitude image).
+Mission supported:
+
+-   Sentinel-1
+-   RCM
+-   RadarSat-2

{grdwindinversion-0.2.7 → grdwindinversion-0.3.2}/docs/conf.py

@@ -17,11 +17,11 @@
 # relative to the documentation root, use os.path.abspath to make it
 # absolute, like shown here.
 #
+import grdwindinversion
 import os
 import sys
 sys.path.insert(0, os.path.abspath('..'))
 
-import grdwindinversion
 
 # -- General configuration ---------------------------------------------
 
@@ -43,6 +43,12 @@ extensions = [
     'sphinxcontrib.mermaid',
     'nbsphinx',
     'jupyter_sphinx',
+    'myst_parser'
+]
+myst_enable_extensions = [
+    "deflist",        # Pour les listes de définitions
+    "linkify",        # Pour rendre les URLs cliquables automatiquement
+    "colon_fence",    # Pour activer des blocs ::: comme Markdown GFM
 ]
 # order by source
 autodoc_member_order = 'bysource'
@@ -69,7 +75,12 @@ templates_path = ['_templates']
 # You can specify multiple suffix as a list of string:
 #
 # source_suffix = ['.rst', '.md']
-source_suffix = '.rst'
+# source_suffix = '.rst'
+
+source_suffix = {
+    '.rst': 'restructuredtext',
+    '.md': 'markdown',
+}
 
 # The master toctree document.
 master_doc = 'index'
@@ -112,8 +123,11 @@ todo_include_todos = False
 # The theme to use for HTML and HTML Help pages.  See the documentation for
 # a list of builtin themes.
 #
-html_theme = 'alabaster'
-#html_theme = 'sphinx_rtd_theme'
+# html_theme = 'alabaster'
+# html_theme = 'sphinx_rtd_theme'
+# html_theme = 'furo'
+html_theme = 'sphinx_rtd_theme'
+
 
 # Theme options are theme-specific and customize the look and feel of a
 # theme further.  For a list of options available for each theme, see the
@@ -125,8 +139,12 @@ html_theme = 'alabaster'
 # relative to this directory. They are copied after the builtin static files,
 # so a file named "default.css" will overwrite the builtin "default.css".
 html_static_path = ['_static']
-#html_style = 'css/grdwindinversion.css'
+html_style = 'css/grdwindinversion.css'
+# html_style = 'css/xsarsea.css'
 
+html_theme_options = {
+    'logo_only': False,
+}
 # -- Options for HTMLHelp output ---------------------------------------
 
 # Output file base name for HTML help builder.
@@ -187,6 +205,3 @@ texinfo_documents = [
      'One line description of project.',
      'Miscellaneous'),
 ]
-
-
-

grdwindinversion-0.3.2/docs/examples/streaks-display.ipynb

@@ -0,0 +1,324 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Compute and display streaks"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# test on an dummy GRD IW product\n",
+    "import xsar,os\n",
+    "import grdwindinversion \n",
+    "filename = xsar.get_test_file('S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE')\n",
+    "outdir = \"/tmp/\"\n",
+    "config_path = os.path.join(os.path.dirname(grdwindinversion.__file__),'config_prod_streaks.yaml')\n",
+    "assert os.path.exists(config_path)\n",
+    "overwrite = True\n",
+    "resolution = '1000m'"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import logging\n",
+    "logger = logging.getLogger('grdwindinversion.gradientFeatures')\n",
+    "logger.setLevel(logging.DEBUG)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## preprocess & load config"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from grdwindinversion.inversion import preprocess "
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "xr_dataset, out_file, config = preprocess(\n",
+    "    filename, outdir, config_path, overwrite, False, resolution)\n",
+    "\n",
+    "model_co = config[\"l2_params\"][\"model_co\"]\n",
+    "model_cross = config[\"l2_params\"][\"model_cross\"]\n",
+    "copol = config[\"l2_params\"][\"copol\"]\n",
+    "crosspol = config[\"l2_params\"][\"crosspol\"]\n",
+    "copol_gmf = config[\"l2_params\"][\"copol_gmf\"]\n",
+    "crosspol_gmf = config[\"l2_params\"][\"crosspol_gmf\"]\n",
+    "dual_pol = config[\"l2_params\"][\"dual_pol\"]\n",
+    "ancillary_name = config[\"ancillary\"]\n",
+    "sensor_longname = config[\"sensor_longname\"]\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## process gradients\n",
+    "- in xr_dataset, it will add variables related to heterogeneity mask an\n",
+    "- it will create xr_dataset_streaks"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from grdwindinversion.inversion import process_gradients\n",
+    "\n",
+    "if config[\"add_gradientsfeatures\"]:\n",
+    "    xr_dataset, xr_dataset_streaks = process_gradients(\n",
+    "        xr_dataset, config)\n",
+    "else:\n",
+    "    xr_dataset_streaks = None"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "xr_dataset_streaks"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## Plotting the results"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "def get_uv_from_dir(wdir):\n",
+    "    \"\"\"\n",
+    "    Get u, v from wind direction\n",
+    "    \n",
+    "    Parameters\n",
+    "    ----------\n",
+    "    wdir : float\n",
+    "        Wind direction in degrees, meteo convention\n",
+    "    \n",
+    "    Returns\n",
+    "    -------\n",
+    "    u : float\n",
+    "        u component of the wind\n",
+    "    v : float\n",
+    "        v component of the wind\n",
+    "    u_norm : float\n",
+    "        Normalized u component of the wind\n",
+    "    v_norm : float\n",
+    "        Normalized v component of the wind  \n",
+    "    \"\"\"\n",
+    "    u = np.cos(np.radians(270-wdir))\n",
+    "    v = np.sin(np.radians(270-wdir))\n",
+    "    sq_mean = np.sqrt(u**2+v**2)\n",
+    "    return u, v, u/sq_mean, v/sq_mean"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import matplotlib.pyplot as plt\n",
+    "import numpy as np\n",
+    "import cartopy.crs as ccrs\n",
+    "from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER\n",
+    "import matplotlib\n",
+    "cmin=0;cmax=30\n",
+    "norm = matplotlib.colors.Normalize(vmin=cmin,vmax=cmax,clip=False)\n",
+    "import os \n",
+    "\n",
+    "#ancillary \n",
+    "ds = xr_dataset\n",
+    "p1 = int(ds.line.size/16)\n",
+    "p2 = int(ds.sample.size/16)\n",
+    "\n",
+    "lons = ds.longitude[::p1,::p2].values.flatten()\n",
+    "lats = ds.latitude[::p1,::p2].values.flatten()\n",
+    "winddir_ancillary = ds.ancillary_wind_direction[::p1,::p2]\n",
+    "u_ancillary, v_ancllary, u_norm_ancillary, v_norm_ancillary = get_uv_from_dir(winddir_ancillary)\n",
+    "\n",
+    "longitude = xr_dataset_streaks['longitude'].values\n",
+    "latitude = xr_dataset_streaks['latitude'].values"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "def plot_wind_from_streaks(streaks_dir, **kwargs):\n",
+    "    \"\"\"\n",
+    "    Plot wind from streaks\n",
+    "    \n",
+    "    Parameters\n",
+    "    ----------\n",
+    "    streaks_dir : xr.Dataset\n",
+    "        Dataset containing streaks direction in degrees in meteo convention\n",
+    "    \n",
+    "    **kwargs : dict\n",
+    "        window_size : int\n",
+    "            Size of the window for smoothing the streaks\n",
+    "        downscale_factor : int\n",
+    "            Downscale factor for smoothing the streaks\n",
+    "        pol : str\n",
+    "            Polarization to plot\n",
+    "        varname : str\n",
+    "            Variable name to plot\n",
+    "    \"\"\"\n",
+    "    \n",
+    "    \n",
+    "    # handle kwargs\n",
+    "    window_size = kwargs.get('window_size', None)\n",
+    "    downscale_factor = kwargs.get('downscale_factor', None)\n",
+    "    pol = kwargs.get('pol', None)\n",
+    "    pol_display = kwargs.get('pol', 'VV')\n",
+    "    varname = kwargs.get('varname', None)\n",
+    "    savefig = kwargs.get('savefig', False)\n",
+    "    savedir = kwargs.get('savedir', './tmp')\n",
+    "    \n",
+    "    # Plot\n",
+    "    fig = plt.figure(figsize=(10, 9))\n",
+    "    ax = plt.axes(projection=ccrs.PlateCarree())\n",
+    "\n",
+    "    # Calculer les composantes u et v à partir des angles streaks_dir\n",
+    "    u, v, u_norm, v_norm = get_uv_from_dir(streaks_dir)\n",
+    "\n",
+    "    # Display streaks direction\n",
+    "    #ax.quiver(longitude, latitude, u_norm, v_norm, edgecolors='k',  norm=norm, pivot= 'mid', scale_units='xy', scale=4., zorder=10, width=0.1/25,transform = ccrs.PlateCarree(), color = 'red', label = 'owiWindStreaks')\n",
+    "    ax.quiver(longitude, latitude, u_norm, v_norm, edgecolors='k',  norm=norm, pivot= 'mid', scale_units='xy', scale=8., zorder=10, width=0.1/25,transform = ccrs.PlateCarree(), color = 'red', label = 'owiWindStreaks')\n",
+    "\n",
+    "    # Display ancillary wind direction\n",
+    "    #ax.quiver(lons, lats, u_norm_ancillary, v_norm_ancillary, edgecolors='k',  norm=norm, pivot= 'mid', scale_units='xy', scale=4., zorder=10, width=0.1/25,transform = ccrs.PlateCarree(), color='blue', label = \"owiAncillaryWindDirection\")\n",
+    "    ax.quiver(lons, lats, u_norm_ancillary, v_norm_ancillary, edgecolors='k',  norm=norm, pivot= 'mid', scale_units='xy', scale=8., zorder=10, width=0.1/25,transform = ccrs.PlateCarree(), color='blue', label = \"owiAncillaryWindDirection\")\n",
+    "\n",
+    "    # Display sigma0 copol\n",
+    "    im = ax.pcolormesh(xr_dataset.longitude,xr_dataset.latitude, xr_dataset.sigma0.sel(pol=pol_display),cmap='gray')\n",
+    "    cbar = plt.colorbar(im, ax =ax)\n",
+    "    cbar.set_label(f'Sigma0, pol = {pol_display}')\n",
+    "    \n",
+    "    # Expand the map by 1 degree in height\n",
+    "    lon_min, lon_max = np.nanmin(lons) - .2, np.nanmax(lons) + .2 \n",
+    "    lat_min, lat_max = np.nanmin(lats) - .2, np.nanmax(lats) + .7\n",
+    "    ax.set_extent([lon_min, lon_max, lat_min, lat_max], crs=ccrs.PlateCarree())\n",
+    "\n",
+    "\n",
+    "    # Configurer les axes\n",
+    "    ax.set_xlabel('Longitude')\n",
+    "    ax.set_ylabel('Latitude')\n",
+    "    ax.set_aspect('equal', adjustable='box')\n",
+    "\n",
+    "    gl = ax.gridlines(draw_labels=True)\n",
+    "    gl.top_labels = False\n",
+    "    gl.right_labels = False\n",
+    "    gl.xformatter = LONGITUDE_FORMATTER\n",
+    "    gl.yformatter = LATITUDE_FORMATTER\n",
+    "    ax.set_title(f\"{os.path.basename(filename)}\\nvarname = {varname}\\nwindow size = {window_size} ; downscale_factor = {downscale_factor} ; pol= {pol}\")\n",
+    "    legend = ax.legend(facecolor='white', loc = 'upper right')\n",
+    "    legend.set_zorder(100)\n",
+    "    if savefig:\n",
+    "        os.makedirs(savedir, exist_ok=True)\n",
+    "        plt.savefig(os.path.join(savedir, f\"wind_streaks_{varname}_{window_size}_{downscale_factor}_{pol}.png\"))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "savefig = False\n",
+    "\n",
+    "varname = 'dir_mean_smooth'\n",
+    "kwargs = {'varname' : varname, 'savefig': savefig, 'savedir': \"./tmp\"}\n",
+    "plot_wind_from_streaks(xr_dataset_streaks[varname], **kwargs)\n",
+    "\n",
+    "varname = 'dir_smooth_mean'\n",
+    "kwargs = {'varname' : varname, 'savefig': savefig, 'savedir': \"./tmp\"}\n",
+    "plot_wind_from_streaks(xr_dataset_streaks[varname], **kwargs)\n",
+    "\n",
+    "\n",
+    "# plot individual solutions \n",
+    "varname = 'dir_smooth'\n",
+    "for downscale_factor in xr_dataset_streaks.downscale_factor.values:\n",
+    "   for pol in xr_dataset_streaks.pol.values:\n",
+    "       for window_size in xr_dataset_streaks.window_size.values:\n",
+    "           streaks_dir = xr_dataset_streaks[varname].sel(\n",
+    "               downscale_factor=downscale_factor, pol=pol, window_size=window_size).values\n",
+    "           kwargs = {'varname' : varname , 'window_size': window_size, 'downscale_factor': downscale_factor, 'pol': pol, 'savefig': savefig, 'savedir': \"./tmp\"}\n",
+    "           plot_wind_from_streaks(streaks_dir, **kwargs)\n",
+    "\n",
+    "# plot one individual solution\n",
+    "# downscale_factor = 1\n",
+    "# pol = 'VV'\n",
+    "# window_size = 3200\n",
+    "\n",
+    "# streaks_dir = xr_dataset_streaks[varname].sel(\n",
+    "#     downscale_factor=downscale_factor, pol=pol, window_size=window_size).values\n",
+    "# kwargs = {'varname' : varname , 'window_size': window_size, 'downscale_factor': downscale_factor, 'pol': pol, 'savefig': True, 'savedir': \"./tmp\"}\n",
+    "# plot_wind_from_streaks(streaks_dir, **kwargs)\n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.11.0"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}