grdwindinversion 0.2.3.post9__tar.gz → 0.2.3.post11__tar.gz

{grdwindinversion-0.2.3.post9/grdwindinversion.egg-info → grdwindinversion-0.2.3.post11}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: grdwindinversion
-Version: 0.2.3.post9
+Version: 0.2.3.post11
 Summary: Package to perform Wind inversion from GRD Level-1 SAR images
 License: MIT
 Classifier: Development Status :: 2 - Pre-Alpha

grdwindinversion-0.2.3.post11/docs/examples/wind-inversion-from-grd.ipynb

@@ -0,0 +1,196 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "e578297a-d2dc-4cb4-b016-1181649d227b",
+   "metadata": {},
+   "source": [
+    "## Example of Wind inversion"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "fa19ad8e-ebaa-430e-bd91-c21d517e7e88",
+   "metadata": {},
+   "source": [
+    "### define a product GRD Sentinel-1"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "id": "ef483f15-4690-460b-b688-fa966bab8cd8",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import xsar"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "id": "e04fa1b0-eb49-4e45-af0f-0be798e26f95",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "/home/vincelhx/Documents/autoentreprise/IFREMER/libs/fork_xsar_vinc/xsar/src/xsar/xsar.py:212: UserWarning: Downloading https://cyclobs.ifremer.fr/static/sarwing_datarmor/xsardata/S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE.zip\n",
+      "  warnings.warn(\"Downloading %s\" % file_url)\n"
+     ]
+    },
+    {
+     "ename": "ClientResponseError",
+     "evalue": "503, message='Service Unavailable', url=URL('https://cyclobs.ifremer.fr/static/sarwing_datarmor/xsardata/S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE.zip')",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mClientResponseError\u001b[0m                       Traceback (most recent call last)",
+      "Cell \u001b[0;32mIn[6], line 2\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[38;5;66;03m# test on an dummy GRD IW product\u001b[39;00m\n\u001b[0;32m----> 2\u001b[0m input_file \u001b[38;5;241m=\u001b[39m \u001b[43mxsar\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mget_test_file\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mS1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m)\u001b[49m\n",
+      "File \u001b[0;32m~/Documents/autoentreprise/IFREMER/libs/fork_xsar_vinc/xsar/src/xsar/xsar.py:213\u001b[0m, in \u001b[0;36mget_test_file\u001b[0;34m(fname)\u001b[0m\n\u001b[1;32m    211\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mexists(os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mjoin(res_path, fname)):\n\u001b[1;32m    212\u001b[0m     warnings\u001b[38;5;241m.\u001b[39mwarn(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mDownloading \u001b[39m\u001b[38;5;132;01m%s\u001b[39;00m\u001b[38;5;124m\"\u001b[39m \u001b[38;5;241m%\u001b[39m file_url)\n\u001b[0;32m--> 213\u001b[0m     local_file \u001b[38;5;241m=\u001b[39m \u001b[43murl_get\u001b[49m\u001b[43m(\u001b[49m\u001b[43mfile_url\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    214\u001b[0m     warnings\u001b[38;5;241m.\u001b[39mwarn(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mUnzipping \u001b[39m\u001b[38;5;132;01m%s\u001b[39;00m\u001b[38;5;124m\"\u001b[39m \u001b[38;5;241m%\u001b[39m os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mjoin(res_path, fname))\n\u001b[1;32m    215\u001b[0m     \u001b[38;5;28;01mwith\u001b[39;00m zipfile\u001b[38;5;241m.\u001b[39mZipFile(local_file, \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mr\u001b[39m\u001b[38;5;124m'\u001b[39m) \u001b[38;5;28;01mas\u001b[39;00m zip_ref:\n",
+      "File \u001b[0;32m~/Documents/autoentreprise/IFREMER/libs/fork_xsar_vinc/xsar/src/xsar/utils.py:664\u001b[0m, in \u001b[0;36murl_get\u001b[0;34m(url, cache_dir)\u001b[0m\n\u001b[1;32m    640\u001b[0m \u001b[38;5;250m\u001b[39m\u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m    641\u001b[0m \u001b[38;5;124;03mGet fil from url, using caching.\u001b[39;00m\n\u001b[1;32m    642\u001b[0m \n\u001b[0;32m   (...)\u001b[0m\n\u001b[1;32m    660\u001b[0m \u001b[38;5;124;03mDue to fsspec, the returned filename won't match the remote one.\u001b[39;00m\n\u001b[1;32m    661\u001b[0m \u001b[38;5;124;03m\"\"\"\u001b[39;00m\n\u001b[1;32m    663\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;124m'\u001b[39m\u001b[38;5;124m://\u001b[39m\u001b[38;5;124m'\u001b[39m \u001b[38;5;129;01min\u001b[39;00m url:\n\u001b[0;32m--> 664\u001b[0m     \u001b[38;5;28;01mwith\u001b[39;00m fsspec\u001b[38;5;241m.\u001b[39mopen(\n\u001b[1;32m    665\u001b[0m             \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mfilecache::\u001b[39m\u001b[38;5;132;01m%s\u001b[39;00m\u001b[38;5;124m'\u001b[39m \u001b[38;5;241m%\u001b[39m url,\n\u001b[1;32m    666\u001b[0m             https\u001b[38;5;241m=\u001b[39m{\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mclient_kwargs\u001b[39m\u001b[38;5;124m'\u001b[39m: {\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mtimeout\u001b[39m\u001b[38;5;124m'\u001b[39m: aiohttp\u001b[38;5;241m.\u001b[39mClientTimeout(total\u001b[38;5;241m=\u001b[39m\u001b[38;5;241m3600\u001b[39m)}},\n\u001b[1;32m    667\u001b[0m             filecache\u001b[38;5;241m=\u001b[39m{\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mcache_storage\u001b[39m\u001b[38;5;124m'\u001b[39m: os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mjoin(os\u001b[38;5;241m.\u001b[39mpath\u001b[38;5;241m.\u001b[39mjoin(config[\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mdata_dir\u001b[39m\u001b[38;5;124m'\u001b[39m], \u001b[38;5;124m'\u001b[39m\u001b[38;5;124mfsspec_cache\u001b[39m\u001b[38;5;124m'\u001b[39m))}\n\u001b[1;32m    668\u001b[0m     ) \u001b[38;5;28;01mas\u001b[39;00m f:\n\u001b[1;32m    669\u001b[0m         fname \u001b[38;5;241m=\u001b[39m f\u001b[38;5;241m.\u001b[39mname\n\u001b[1;32m    670\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/core.py:102\u001b[0m, in \u001b[0;36mOpenFile.__enter__\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m     99\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21m__enter__\u001b[39m(\u001b[38;5;28mself\u001b[39m):\n\u001b[1;32m    100\u001b[0m     mode \u001b[38;5;241m=\u001b[39m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mmode\u001b[38;5;241m.\u001b[39mreplace(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mt\u001b[39m\u001b[38;5;124m\"\u001b[39m, \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m\"\u001b[39m)\u001b[38;5;241m.\u001b[39mreplace(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mb\u001b[39m\u001b[38;5;124m\"\u001b[39m, \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m\"\u001b[39m) \u001b[38;5;241m+\u001b[39m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mb\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[0;32m--> 102\u001b[0m     f \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mfs\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mopen\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mpath\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mmode\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmode\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    104\u001b[0m     \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mfobjects \u001b[38;5;241m=\u001b[39m [f]\n\u001b[1;32m    106\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mcompression \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/implementations/cached.py:444\u001b[0m, in \u001b[0;36mCachingFileSystem.__getattribute__.<locals>.<lambda>\u001b[0;34m(*args, **kw)\u001b[0m\n\u001b[1;32m    408\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21m__getattribute__\u001b[39m(\u001b[38;5;28mself\u001b[39m, item):\n\u001b[1;32m    409\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m item \u001b[38;5;129;01min\u001b[39;00m [\n\u001b[1;32m    410\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mload_cache\u001b[39m\u001b[38;5;124m\"\u001b[39m,\n\u001b[1;32m    411\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m_open\u001b[39m\u001b[38;5;124m\"\u001b[39m,\n\u001b[0;32m   (...)\u001b[0m\n\u001b[1;32m    442\u001b[0m         \u001b[38;5;66;03m# all the methods defined in this class. Note `open` here, since\u001b[39;00m\n\u001b[1;32m    443\u001b[0m         \u001b[38;5;66;03m# it calls `_open`, but is actually in superclass\u001b[39;00m\n\u001b[0;32m--> 444\u001b[0m         \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;01mlambda\u001b[39;00m \u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkw: \u001b[38;5;28;43mgetattr\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mtype\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mitem\u001b[49m\u001b[43m)\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[38;5;21;43m__get__\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m    445\u001b[0m \u001b[43m            \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43margs\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkw\u001b[49m\n\u001b[1;32m    446\u001b[0m \u001b[43m        \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    447\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m item \u001b[38;5;129;01min\u001b[39;00m [\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m__reduce_ex__\u001b[39m\u001b[38;5;124m\"\u001b[39m]:\n\u001b[1;32m    448\u001b[0m         \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mAttributeError\u001b[39;00m\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/spec.py:1151\u001b[0m, in \u001b[0;36mAbstractFileSystem.open\u001b[0;34m(self, path, mode, block_size, cache_options, compression, **kwargs)\u001b[0m\n\u001b[1;32m   1149\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n\u001b[1;32m   1150\u001b[0m     ac \u001b[38;5;241m=\u001b[39m kwargs\u001b[38;5;241m.\u001b[39mpop(\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mautocommit\u001b[39m\u001b[38;5;124m\"\u001b[39m, \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_intrans)\n\u001b[0;32m-> 1151\u001b[0m     f \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_open\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m   1152\u001b[0m \u001b[43m        \u001b[49m\u001b[43mpath\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1153\u001b[0m \u001b[43m        \u001b[49m\u001b[43mmode\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mmode\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1154\u001b[0m \u001b[43m        \u001b[49m\u001b[43mblock_size\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mblock_size\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1155\u001b[0m \u001b[43m        \u001b[49m\u001b[43mautocommit\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mac\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1156\u001b[0m \u001b[43m        \u001b[49m\u001b[43mcache_options\u001b[49m\u001b[38;5;241;43m=\u001b[39;49m\u001b[43mcache_options\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1157\u001b[0m \u001b[43m        \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m,\u001b[49m\n\u001b[1;32m   1158\u001b[0m \u001b[43m    \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m   1159\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m compression \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[1;32m   1160\u001b[0m         \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01mfsspec\u001b[39;00m\u001b[38;5;21;01m.\u001b[39;00m\u001b[38;5;21;01mcompression\u001b[39;00m \u001b[38;5;28;01mimport\u001b[39;00m compr\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/implementations/cached.py:444\u001b[0m, in \u001b[0;36mCachingFileSystem.__getattribute__.<locals>.<lambda>\u001b[0;34m(*args, **kw)\u001b[0m\n\u001b[1;32m    408\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21m__getattribute__\u001b[39m(\u001b[38;5;28mself\u001b[39m, item):\n\u001b[1;32m    409\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m item \u001b[38;5;129;01min\u001b[39;00m [\n\u001b[1;32m    410\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mload_cache\u001b[39m\u001b[38;5;124m\"\u001b[39m,\n\u001b[1;32m    411\u001b[0m         \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m_open\u001b[39m\u001b[38;5;124m\"\u001b[39m,\n\u001b[0;32m   (...)\u001b[0m\n\u001b[1;32m    442\u001b[0m         \u001b[38;5;66;03m# all the methods defined in this class. Note `open` here, since\u001b[39;00m\n\u001b[1;32m    443\u001b[0m         \u001b[38;5;66;03m# it calls `_open`, but is actually in superclass\u001b[39;00m\n\u001b[0;32m--> 444\u001b[0m         \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;01mlambda\u001b[39;00m \u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkw: \u001b[38;5;28;43mgetattr\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mtype\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mitem\u001b[49m\u001b[43m)\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[38;5;21;43m__get__\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m(\u001b[49m\n\u001b[1;32m    445\u001b[0m \u001b[43m            \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43margs\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkw\u001b[49m\n\u001b[1;32m    446\u001b[0m \u001b[43m        \u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    447\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m item \u001b[38;5;129;01min\u001b[39;00m [\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124m__reduce_ex__\u001b[39m\u001b[38;5;124m\"\u001b[39m]:\n\u001b[1;32m    448\u001b[0m         \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mAttributeError\u001b[39;00m\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/implementations/cached.py:695\u001b[0m, in \u001b[0;36mWholeFileCacheFileSystem._open\u001b[0;34m(self, path, mode, **kwargs)\u001b[0m\n\u001b[1;32m    693\u001b[0m             f2\u001b[38;5;241m.\u001b[39mwrite(data)\n\u001b[1;32m    694\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n\u001b[0;32m--> 695\u001b[0m     \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mfs\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mget\u001b[49m\u001b[43m(\u001b[49m\u001b[43mpath\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mfn\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    696\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39msave_cache()\n\u001b[1;32m    697\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_open(path, mode)\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/asyn.py:115\u001b[0m, in \u001b[0;36msync_wrapper.<locals>.wrapper\u001b[0;34m(*args, **kwargs)\u001b[0m\n\u001b[1;32m    112\u001b[0m \u001b[38;5;129m@functools\u001b[39m\u001b[38;5;241m.\u001b[39mwraps(func)\n\u001b[1;32m    113\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mwrapper\u001b[39m(\u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs):\n\u001b[1;32m    114\u001b[0m     \u001b[38;5;28mself\u001b[39m \u001b[38;5;241m=\u001b[39m obj \u001b[38;5;129;01mor\u001b[39;00m args[\u001b[38;5;241m0\u001b[39m]\n\u001b[0;32m--> 115\u001b[0m     \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[43msync\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mloop\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mfunc\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43margs\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[38;5;241;43m*\u001b[39;49m\u001b[43mkwargs\u001b[49m\u001b[43m)\u001b[49m\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/asyn.py:100\u001b[0m, in \u001b[0;36msync\u001b[0;34m(loop, func, timeout, *args, **kwargs)\u001b[0m\n\u001b[1;32m     98\u001b[0m     \u001b[38;5;28;01mraise\u001b[39;00m FSTimeoutError \u001b[38;5;28;01mfrom\u001b[39;00m \u001b[38;5;21;01mreturn_result\u001b[39;00m\n\u001b[1;32m     99\u001b[0m \u001b[38;5;28;01melif\u001b[39;00m \u001b[38;5;28misinstance\u001b[39m(return_result, \u001b[38;5;167;01mBaseException\u001b[39;00m):\n\u001b[0;32m--> 100\u001b[0m     \u001b[38;5;28;01mraise\u001b[39;00m return_result\n\u001b[1;32m    101\u001b[0m \u001b[38;5;28;01melse\u001b[39;00m:\n\u001b[1;32m    102\u001b[0m     \u001b[38;5;28;01mreturn\u001b[39;00m return_result\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/asyn.py:55\u001b[0m, in \u001b[0;36m_runner\u001b[0;34m(event, coro, result, timeout)\u001b[0m\n\u001b[1;32m     53\u001b[0m     coro \u001b[38;5;241m=\u001b[39m asyncio\u001b[38;5;241m.\u001b[39mwait_for(coro, timeout\u001b[38;5;241m=\u001b[39mtimeout)\n\u001b[1;32m     54\u001b[0m \u001b[38;5;28;01mtry\u001b[39;00m:\n\u001b[0;32m---> 55\u001b[0m     result[\u001b[38;5;241m0\u001b[39m] \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mawait\u001b[39;00m coro\n\u001b[1;32m     56\u001b[0m \u001b[38;5;28;01mexcept\u001b[39;00m \u001b[38;5;167;01mException\u001b[39;00m \u001b[38;5;28;01mas\u001b[39;00m ex:\n\u001b[1;32m     57\u001b[0m     result[\u001b[38;5;241m0\u001b[39m] \u001b[38;5;241m=\u001b[39m ex\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/asyn.py:562\u001b[0m, in \u001b[0;36mAsyncFileSystem._get\u001b[0;34m(self, rpath, lpath, recursive, callback, **kwargs)\u001b[0m\n\u001b[1;32m    560\u001b[0m     callback\u001b[38;5;241m.\u001b[39mbranch(rpath, lpath, kwargs)\n\u001b[1;32m    561\u001b[0m     coros\u001b[38;5;241m.\u001b[39mappend(\u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39m_get_file(rpath, lpath, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs))\n\u001b[0;32m--> 562\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;01mawait\u001b[39;00m _run_coros_in_chunks(\n\u001b[1;32m    563\u001b[0m     coros, batch_size\u001b[38;5;241m=\u001b[39mbatch_size, callback\u001b[38;5;241m=\u001b[39mcallback\n\u001b[1;32m    564\u001b[0m )\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/asyn.py:246\u001b[0m, in \u001b[0;36m_run_coros_in_chunks\u001b[0;34m(coros, batch_size, callback, timeout, return_exceptions, nofiles)\u001b[0m\n\u001b[1;32m    240\u001b[0m     \u001b[38;5;28;01mif\u001b[39;00m callback \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m _DEFAULT_CALLBACK:\n\u001b[1;32m    241\u001b[0m         [\n\u001b[1;32m    242\u001b[0m             t\u001b[38;5;241m.\u001b[39madd_done_callback(\u001b[38;5;28;01mlambda\u001b[39;00m \u001b[38;5;241m*\u001b[39m_, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39m__: callback\u001b[38;5;241m.\u001b[39mrelative_update(\u001b[38;5;241m1\u001b[39m))\n\u001b[1;32m    243\u001b[0m             \u001b[38;5;28;01mfor\u001b[39;00m t \u001b[38;5;129;01min\u001b[39;00m chunk\n\u001b[1;32m    244\u001b[0m         ]\n\u001b[1;32m    245\u001b[0m     results\u001b[38;5;241m.\u001b[39mextend(\n\u001b[0;32m--> 246\u001b[0m         \u001b[38;5;28;01mawait\u001b[39;00m asyncio\u001b[38;5;241m.\u001b[39mgather(\u001b[38;5;241m*\u001b[39mchunk, return_exceptions\u001b[38;5;241m=\u001b[39mreturn_exceptions),\n\u001b[1;32m    247\u001b[0m     )\n\u001b[1;32m    248\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m results\n",
+      "File \u001b[0;32m~/miniconda3/envs/xsar_N3_local/lib/python3.10/asyncio/tasks.py:408\u001b[0m, in \u001b[0;36mwait_for\u001b[0;34m(fut, timeout)\u001b[0m\n\u001b[1;32m    405\u001b[0m loop \u001b[38;5;241m=\u001b[39m events\u001b[38;5;241m.\u001b[39mget_running_loop()\n\u001b[1;32m    407\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m timeout \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m:\n\u001b[0;32m--> 408\u001b[0m     \u001b[38;5;28;01mreturn\u001b[39;00m \u001b[38;5;28;01mawait\u001b[39;00m fut\n\u001b[1;32m    410\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m timeout \u001b[38;5;241m<\u001b[39m\u001b[38;5;241m=\u001b[39m \u001b[38;5;241m0\u001b[39m:\n\u001b[1;32m    411\u001b[0m     fut \u001b[38;5;241m=\u001b[39m ensure_future(fut, loop\u001b[38;5;241m=\u001b[39mloop)\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/implementations/http.py:248\u001b[0m, in \u001b[0;36mHTTPFileSystem._get_file\u001b[0;34m(self, rpath, lpath, chunk_size, callback, **kwargs)\u001b[0m\n\u001b[1;32m    245\u001b[0m     size \u001b[38;5;241m=\u001b[39m \u001b[38;5;28;01mNone\u001b[39;00m\n\u001b[1;32m    247\u001b[0m callback\u001b[38;5;241m.\u001b[39mset_size(size)\n\u001b[0;32m--> 248\u001b[0m \u001b[38;5;28;43mself\u001b[39;49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43m_raise_not_found_for_status\u001b[49m\u001b[43m(\u001b[49m\u001b[43mr\u001b[49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[43mrpath\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m    249\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m isfilelike(lpath):\n\u001b[1;32m    250\u001b[0m     outfile \u001b[38;5;241m=\u001b[39m lpath\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/fsspec/implementations/http.py:214\u001b[0m, in \u001b[0;36mHTTPFileSystem._raise_not_found_for_status\u001b[0;34m(self, response, url)\u001b[0m\n\u001b[1;32m    212\u001b[0m \u001b[38;5;28;01mif\u001b[39;00m response\u001b[38;5;241m.\u001b[39mstatus \u001b[38;5;241m==\u001b[39m \u001b[38;5;241m404\u001b[39m:\n\u001b[1;32m    213\u001b[0m     \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mFileNotFoundError\u001b[39;00m(url)\n\u001b[0;32m--> 214\u001b[0m \u001b[43mresponse\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mraise_for_status\u001b[49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\n",
+      "File \u001b[0;32m~/.local/lib/python3.10/site-packages/aiohttp/client_reqrep.py:1005\u001b[0m, in \u001b[0;36mClientResponse.raise_for_status\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m   1003\u001b[0m \u001b[38;5;28;01massert\u001b[39;00m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mreason \u001b[38;5;129;01mis\u001b[39;00m \u001b[38;5;129;01mnot\u001b[39;00m \u001b[38;5;28;01mNone\u001b[39;00m\n\u001b[1;32m   1004\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mrelease()\n\u001b[0;32m-> 1005\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m ClientResponseError(\n\u001b[1;32m   1006\u001b[0m     \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mrequest_info,\n\u001b[1;32m   1007\u001b[0m     \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mhistory,\n\u001b[1;32m   1008\u001b[0m     status\u001b[38;5;241m=\u001b[39m\u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mstatus,\n\u001b[1;32m   1009\u001b[0m     message\u001b[38;5;241m=\u001b[39m\u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mreason,\n\u001b[1;32m   1010\u001b[0m     headers\u001b[38;5;241m=\u001b[39m\u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mheaders,\n\u001b[1;32m   1011\u001b[0m )\n",
+      "\u001b[0;31mClientResponseError\u001b[0m: 503, message='Service Unavailable', url=URL('https://cyclobs.ifremer.fr/static/sarwing_datarmor/xsardata/S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE.zip')"
+     ]
+    }
+   ],
+   "source": [
+    "# test on an dummy GRD IW product\n",
+    "input_file = xsar.get_test_file('S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE')\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "f4886386-b19a-4c43-a790-192f4866a692",
+   "metadata": {},
+   "source": [
+    "### perform wind inversion"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d3fda3a2-2ee6-4401-be8f-c0fab8344a49",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from grdwindinversion.inversion import makeL2 \n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "c8924754-34f5-4a27-a598-492022b3b3f5",
+   "metadata": {},
+   "source": [
+    "### define config file"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "f95ee426-0b46-4d39-ae90-241e8bbaa332",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "True\n"
+     ]
+    }
+   ],
+   "source": [
+    "import grdwindinversion\n",
+    "import os\n",
+    "config_file_s1 = os.path.join(os.path.dirname(grdwindinversion.__file__),'config_prod.yaml')\n",
+    "grdwindinversion.__file__\n",
+    "print(os.path.exists(config_file_s1))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "3cbe649d-44d7-4a7f-9041-137092a61e07",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "out_folder = '/tmp/test_windinversion'"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "b90c78bb-8643-4d64-a463-28fcb68d9b38",
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "NameError",
+     "evalue": "name 'input_file' is not defined",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
+      "Cell \u001b[0;32mIn[4], line 1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m out_file,outds \u001b[38;5;241m=\u001b[39m makeL2(\u001b[43minput_file\u001b[49m,out_folder, config_path\u001b[38;5;241m=\u001b[39mconfig_file_s1,overwrite\u001b[38;5;241m=\u001b[39m\u001b[38;5;28;01mTrue\u001b[39;00m)\n\u001b[1;32m      2\u001b[0m outds\n",
+      "\u001b[0;31mNameError\u001b[0m: name 'input_file' is not defined"
+     ]
+    }
+   ],
+   "source": [
+    "out_file,outds = makeL2(input_file,out_folder, config_path=config_file_s1,overwrite=True)\n",
+    "outds"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "91137cfe",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "9a83becc",
+   "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.10.12"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

{grdwindinversion-0.2.3.post9 → grdwindinversion-0.2.3.post11}/grdwindinversion/inversion.py

@@ -219,12 +219,12 @@ def inverse(dual_pol, inc, sigma0, sigma0_dual, ancillary_wind, dsig_cr, model_v
     Parameters
     ----------
     dual_pol: bool
-        True if dualpol, False if singlepol 
+        True if dualpol, False if singlepol
     inc: xarray.DataArray
         incidence angle
     sigma0: xarray.DataArray
         sigma0 to be inverted
-    sigma0_dual: xarray.DataArray 
+    sigma0_dual: xarray.DataArray
         sigma0 to be inverted for dualpol
     ancillary_wind=: xarray.DataArray (numpy.complex28)
         ancillary wind
@@ -232,10 +232,10 @@ def inverse(dual_pol, inc, sigma0, sigma0_dual, ancillary_wind, dsig_cr, model_v
     dsig_cr=: float or xarray.DataArray
         parameters used for
 
-            | `Jsig_cr=((sigma0_gmf - sigma0) / dsig_cr) ** 2`    
-    model_vv=: str 
+            | `Jsig_cr=((sigma0_gmf - sigma0) / dsig_cr) ** 2`
+    model_vv=: str
         model to use for VV or HH polarization.
-    model_vh=: str 
+    model_vh=: str
         model to use for VH or HV polarization.
 
     Returns
@@ -245,8 +245,8 @@ def inverse(dual_pol, inc, sigma0, sigma0_dual, ancillary_wind, dsig_cr, model_v
 
     See Also
     --------
-    xsarsea documentation 
-    https://cyclobs.ifremer.fr/static/sarwing_datarmor/xsarsea/examples/windspeed_inversion.html    
+    xsarsea documentation
+    https://cyclobs.ifremer.fr/static/sarwing_datarmor/xsarsea/examples/windspeed_inversion.html
     """
     logging.debug("inversion")
 
@@ -306,7 +306,7 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
     Returns
     -------
     xarray.Dataset
-        final dataset 
+        final dataset
     dict
         encoding dict
 
@@ -321,17 +321,33 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
         'elevation': 'owiElevationAngle',
         'ground_heading': 'owiHeading',
         'land_mask': 'owiLandFlag',
-        'offboresight': 'owiOffBoresightAngle',
         'mask': 'owiMask',
         'windspeed_co': 'owiWindSpeed_co',
-        'windspeed_cross': 'owiWindSpeed_cross',
-        'windspeed_dual': 'owiWindSpeed',
         'winddir_co': 'owiWindDirection_co',
-        'winddir_cross': 'owiWindDirection_cross',
-        'winddir_dual': 'owiWindDirection',
         'ancillary_wind_speed': 'owiAncillaryWindSpeed',
         'ancillary_wind_direction': 'owiAncillaryWindDirection',
     })
+
+    if "offboresight" in xr_dataset:
+        xr_dataset = xr_dataset.rename(
+            {"offboresight": "owiOffboresightAngle"})
+
+    xr_dataset.owiLon.attrs["units"] = "degrees_east"
+    xr_dataset.owiLon.attrs["long_name"] = "Longitude at wind cell center"
+    xr_dataset.owiLon.attrs["standard_name"] = "longitude"
+
+    xr_dataset.owiLat.attrs["units"] = "degrees_north"
+    xr_dataset.owiLat.attrs["long_name"] = "Latitude at wind cell center"
+    xr_dataset.owiLat.attrs["standard_name"] = "latitude"
+
+    xr_dataset.owiIncidenceAngle.attrs["units"] = "degrees"
+    xr_dataset.owiIncidenceAngle.attrs["long_name"] = "Incidence angle at wind cell center"
+    xr_dataset.owiIncidenceAngle.attrs["standard_name"] = "incidence"
+
+    xr_dataset.owiElevationAngle.attrs["units"] = "degrees"
+    xr_dataset.owiElevationAngle.attrs["long_name"] = "Elevation angle at wind cell center"
+    xr_dataset.owiElevationAngle.attrs["standard_name"] = "elevation"
+
     xr_dataset['owiNrcs'] = xr_dataset['sigma0_ocean'].sel(pol=copol)
     xr_dataset.owiNrcs.attrs = xr_dataset.sigma0_ocean.attrs
     xr_dataset.owiNrcs.attrs['units'] = 'm^2 / m^2'
@@ -341,7 +357,8 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
     # NESZ & DSIG
     xr_dataset = xr_dataset.assign(
         owiNesz=(['line', 'sample'], xr_dataset.nesz.sel(pol=copol).values))
-    xr_dataset.owiNesz.attrs = xr_dataset.nesz.attrs
+    xr_dataset.owiNesz.attrs['units'] = 'm^2 / m^2'
+    xr_dataset.owiNesz.attrs['long_name'] = 'Noise Equivalent SigmaNaught'
 
     xr_dataset['owiNrcs_no_noise_correction'] = xr_dataset['sigma0_ocean_raw'].sel(
         pol=copol)
@@ -358,6 +375,9 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
             'swath_number_flag': 'owiSwathNumberFlag'
         })
 
+        xr_dataset["owiSwathNumber"].attrs["standart_name"] = "swath number"
+
+    # sigma0_raw__corrected cross
     if "sigma0_raw__corrected" in xr_dataset:
         xr_dataset['owiNrcs_no_noise_correction_recalibrated'] = xr_dataset['sigma0_raw__corrected'].sel(
             pol=copol)
@@ -374,19 +394,25 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
 
         xr_dataset = xr_dataset.rename({
             'dsig_cross': 'owiDsig_cross',
-            'nesz_cross_final': 'owiNesz_cross_final'
+            'nesz_cross_final': 'owiNesz_cross_final',
+            'winddir_cross': 'owiWindDirection_cross',
+            'winddir_dual': 'owiWindDirection',
+            'windspeed_cross': 'owiWindSpeed_cross',
+            'windspeed_dual': 'owiWindSpeed',
         })
-
+        # nrcs cross
         xr_dataset['owiNrcs_cross'] = xr_dataset['sigma0_ocean'].sel(
             pol=crosspol)
         xr_dataset.owiNrcs_cross.attrs['units'] = 'm^2 / m^2'
         xr_dataset.owiNrcs_cross.attrs['long_name'] = 'Normalized Radar Cross Section'
         xr_dataset.owiNrcs_cross.attrs['definition'] = 'owiNrcs_cross_no_noise_correction - owiNesz_cross'
 
+        # nesz cross
         xr_dataset = xr_dataset.assign(owiNesz_cross=(
             ['line', 'sample'], xr_dataset.nesz.sel(pol=crosspol).values))  # no flattening
-        xr_dataset.owiNesz_cross.attrs = xr_dataset.nesz.attrs
-        # unused
+        xr_dataset.owiNesz_cross.attrs['units'] = 'm^2 / m^2'
+        xr_dataset.owiNesz_cross.attrs['long_name'] = 'Noise Equivalent SigmaNaught'
+
         xr_dataset['owiNrcs_cross_no_noise_correction'] = xr_dataset['sigma0_ocean_raw'].sel(
             pol=crosspol)
 
@@ -394,6 +420,7 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
         xr_dataset.owiNrcs_cross_no_noise_correction.attrs[
             'long_name'] = 'Normalized Radar Cross Section, no noise correction applied'
 
+        #  sigma0_raw__corrected cross
         if "sigma0_raw__corrected" in xr_dataset:
             xr_dataset['owiNrcs_cross_no_noise_correction_recalibrated'] = xr_dataset['sigma0_raw__corrected'].sel(
                 pol=crosspol)
@@ -406,23 +433,7 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
 
             xr_dataset.owiNrcs_cross.attrs['definition'] = 'owiNrcs_cross_no_noise_correction_recalibrated - owiNesz_cross'
 
-    xr_dataset["owiWindSpeed_co"].attrs["comment"] = xr_dataset["owiWindSpeed_co"].attrs["comment"].replace(
-        "wind speed and direction", "wind speed")
-
-    xr_dataset["owiWindDirection_co"].attrs["comment"] = "wind direction in meteorological convention, 0=North, 90=East"
-
-    if dual_pol:
-        xr_dataset["owiWindSpeed"].attrs["comment"] = xr_dataset["owiWindSpeed"].attrs["comment"].replace(
-            "wind speed and direction", "wind speed")
-        xr_dataset["owiWindDirection"].attrs["comment"] = "wind direction in meteorological convention, 0=North, 90=East"
-
-        xr_dataset["owiWindSpeed_cross"].attrs['comment'] = "wind speed inverted from model %s (%s)" % (
-            crosspol_gmf, crosspol)
-
-        xr_dataset["owiWindDirection_cross"].attrs["comment"] = "wind direction in meteorological convention, 0=North, 90=East, copied from dualpol"
-
-        xr_dataset.owiWindSpeed_cross.attrs['model'] = crosspol_gmf
-        xr_dataset.owiWindSpeed_cross.attrs['units'] = 'm/s'
+    #  other variables
 
     xr_dataset['owiWindQuality'] = xr.full_like(xr_dataset.owiNrcs, 0)
     xr_dataset['owiWindQuality'].attrs[
@@ -431,7 +442,7 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
     xr_dataset['owiWindQuality'].attrs['flag_values'] = np.array([
         0, 1, 2, 3])
     xr_dataset['owiWindQuality'].attrs['flag_meanings'] = "good medium low poor"
-    xr_dataset['owiWindQuality'].attrs['comment'] = 'not done yet'
+    xr_dataset['owiWindQuality'].attrs['comment'] = 'NOT COMPUTED YET'
 
     xr_dataset['owiWindFilter'] = xr.full_like(xr_dataset.owiNrcs, 0)
     xr_dataset['owiWindFilter'].attrs['long_name'] = "Quality flag taking into account the local heterogeneity"
@@ -440,7 +451,7 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
         0, 1, 2, 3])
     xr_dataset['owiWindFilter'].attrs[
         'flag_meanings'] = "homogeneous_NRCS, heterogeneous_from_co-polarization_NRCS, heterogeneous_from_cross-polarization_NRCS, heterogeneous_from_dual-polarization_NRCS"
-    xr_dataset['owiWindFilter'].attrs['comment'] = 'not done yet'
+    xr_dataset['owiWindFilter'].attrs['comment'] = 'NOT COMPUTED YET'
 
     xr_dataset = xr_dataset.rename(
         {"line": "owiAzSize", "sample": "owiRaSize"})
@@ -455,50 +466,6 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
 
     xr_dataset.compute()
 
-    for var in ['footprint', 'multidataset', 'rawDataStartTime', 'specialHandlingRequired']:
-        if var in xr_dataset.attrs:
-            xr_dataset.attrs[var] = str(xr_dataset.attrs[var])
-        if "approx_transform" in xr_dataset.attrs:
-            del xr_dataset.attrs["approx_transform"]
-
-    xr_dataset.attrs["TITLE"] = "Sentinel-1 OWI Component"
-    xr_dataset.attrs["missionPhase"] = "Test"
-    xr_dataset.attrs["acquisitionStation"] = "/"
-    xr_dataset.attrs["softwareVersion"] = "/"
-    xr_dataset.attrs["pythonVersion"] = str(
-        sys.version_info.major)+'.'+str(sys.version_info.minor)
-    xr_dataset.attrs["polarisationRatio"] = "/"
-    xr_dataset.attrs["l2ProcessingUtcTime"] = datetime.datetime.now().strftime(
-        "%Y-%m-%dT%H:%M:%SZ")
-    xr_dataset.attrs["processingCenter"] = "/"
-    try:
-        xr_dataset.attrs["firstMeasurementTime"] = datetime.datetime.strptime(xr_dataset.attrs['start_date'],
-                                                                              "%Y-%m-%d %H:%M:%S.%f").strftime(
-            "%Y-%m-%dT%H:%M:%SZ")
-        xr_dataset.attrs["lastMeasurementTime"] = datetime.datetime.strptime(xr_dataset.attrs['stop_date'],
-                                                                             "%Y-%m-%d %H:%M:%S.%f").strftime(
-            "%Y-%m-%dT%H:%M:%SZ")
-    except:
-        xr_dataset.attrs["firstMeasurementTime"] = datetime.datetime.strptime(xr_dataset.attrs['start_date'],
-                                                                              "%Y-%m-%d %H:%M:%S").strftime(
-            "%Y-%m-%dT%H:%M:%SZ")
-        xr_dataset.attrs["lastMeasurementTime"] = datetime.datetime.strptime(xr_dataset.attrs['stop_date'],
-                                                                             "%Y-%m-%d %H:%M:%S").strftime(
-            "%Y-%m-%dT%H:%M:%SZ")
-    xr_dataset.attrs["clmSource"] = "/"
-    xr_dataset.attrs["bathySource"] = "/"
-    xr_dataset.attrs['oswAlgorithmName'] = 'grdwindinversion'
-    xr_dataset.attrs["owiAlgorithmVersion"] = grdwindinversion.__version__
-    xr_dataset.attrs["gmf"] = config['GMF_'+copol_gmf+'_NAME'] + \
-        ", " + config["GMF_"+crosspol_gmf+"_NAME"]
-    xr_dataset.attrs["iceSource"] = "/"
-    xr_dataset.attrs["owiNoiseCorrection"] = "False"
-    xr_dataset.attrs["inversionTabGMF"] = config['GMF_'+copol_gmf +
-                                                 '_NAME'] + ", " + config["GMF_"+crosspol_gmf+"_NAME"]
-    xr_dataset.attrs["wnf_3km_average"] = "/"
-    xr_dataset.attrs["owiWindSpeedSrc"] = "owiWindSpeed"
-    xr_dataset.attrs["owiWindDirectionSrc"] = "/"
-
     table_fillValue = {
         "owiWindQuality": -1,
         "owiHeading": 9999.99,
@@ -523,9 +490,6 @@ def makeL2asOwi(xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf,
             else:
                 encoding[var].update({'_FillValue': None})
 
-    xr_dataset.attrs["xsar_version"] = xsar.__version__
-    xr_dataset.attrs["xsarsea_version"] = xsarsea.__version__
-
     return xr_dataset, encoding
 
 
@@ -611,21 +575,7 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
                 xr_dataset = xsar_dataset.datatree['measurement'].to_dataset()
 
         xr_dataset = xr_dataset.rename(map_model)
-        # add attributes
         xr_dataset.attrs = xsar_dataset.dataset.attrs
-        xr_dataset.attrs['L1_path'] = xr_dataset.attrs.pop('name')
-        xr_dataset.attrs["sourceProduct"] = sensor
-        xr_dataset.attrs["missionName"] = sensor_longname
-        if ((recalibration) & ("SENTINEL" in sensor_longname)):
-            xr_dataset.attrs["path_aux_pp1_new"] = os.path.basename(os.path.dirname(
-                os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_pp1_new'])))
-            xr_dataset.attrs["path_aux_cal_new"] = os.path.basename(os.path.dirname(
-                os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_cal_new'])))
-
-            xr_dataset.attrs["path_aux_pp1_old"] = os.path.basename(os.path.dirname(
-                os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_pp1_old'])))
-            xr_dataset.attrs["path_aux_cal_old"] = os.path.basename(os.path.dirname(
-                os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_cal_old'])))
 
     except Exception as e:
         logging.info('%s', traceback.format_exc())
@@ -663,8 +613,33 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
     variables = list(set(xr_dataset) - set(black_list))
     xr_dataset = xr_dataset[variables]
 
-    # TODO Better land mask
-    # xr_dataset.land_mask.values = cv2.dilate(xr_dataset['land_mask'].values.astype('uint8'),np.ones((3,3),np.uint8),iterations = 3)
+    #  lon/lat
+    xr_dataset.longitude.attrs["units"] = "degrees_east"
+    xr_dataset.longitude.attrs["long_name"] = "Longitude at wind cell center"
+    xr_dataset.longitude.attrs["standard_name"] = "longitude"
+
+    xr_dataset.latitude.attrs["units"] = "degrees_north"
+    xr_dataset.latitude.attrs["long_name"] = "Latitude at wind cell center"
+    xr_dataset.latitude.attrs["standard_name"] = "latitude"
+
+    #  incidence
+    xr_dataset.incidence.attrs["units"] = "degrees"
+    xr_dataset.incidence.attrs["long_name"] = "Incidence angle at wind cell center"
+    xr_dataset.incidence.attrs["standard_name"] = "incidence"
+
+    #  elevation
+    xr_dataset.elevation.attrs["units"] = "degrees"
+    xr_dataset.elevation.attrs["long_name"] = "Elevation angle at wind cell center"
+    xr_dataset.elevation.attrs["standard_name"] = "elevation"
+
+    # offboresight
+    # TOREMOVE
+    if "offboresight" in xr_dataset:
+        xr_dataset.offboresight.attrs["units"] = "degrees"
+        xr_dataset.offboresight.attrs["long_name"] = "Offboresight angle at wind cell center"
+        xr_dataset.elevation.attrs["standard_name"] = "offboresight"
+
+    # masks (no ice / no_valid)
     xr_dataset.land_mask.values = binary_dilation(xr_dataset['land_mask'].values.astype('uint8'),
                                                   structure=np.ones((3, 3), np.uint8), iterations=3)
     xr_dataset.land_mask.attrs['long_name'] = 'Mask of data'
@@ -672,8 +647,6 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
     xr_dataset.land_mask.attrs['flag_values'] = np.array([0, 1])
     xr_dataset.land_mask.attrs['flag_meanings'] = 'valid no_valid'
 
-    # MASK
-    # Careful : in sarwing process sometimes there are 2 & 3. Not made here
     logging.debug("mask is a copy of land_mask")
 
     xr_dataset['mask'] = xr.DataArray(xr_dataset.land_mask)
@@ -683,7 +656,7 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
     xr_dataset.mask.attrs['flag_values'] = np.array([0, 1, 2, 3])
     xr_dataset.mask.attrs['flag_meanings'] = 'valid land ice no_valid'
 
-    # ANCILLARY
+    # ancillary
     xr_dataset['ancillary_wind_direction'] = (
         90. - np.rad2deg(np.arctan2(xr_dataset.model_V10, xr_dataset.model_U10)) + 180) % 360
 
@@ -715,13 +688,15 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
         1].strip()
     xr_dataset = xr_dataset.drop_vars(['model_U10', 'model_V10'])
 
-    # NRCS & NESZ
+    # nrcs processing
     xr_dataset['sigma0_ocean'] = xr.where(xr_dataset['mask'], np.nan,
                                           xr_dataset['sigma0'].compute()).transpose(*xr_dataset['sigma0'].dims)
     xr_dataset['sigma0_ocean'] = xr.where(
         xr_dataset['sigma0_ocean'] <= 0, np.nan, xr_dataset['sigma0_ocean'])
 
     xr_dataset['sigma0_ocean'].attrs = xr_dataset['sigma0'].attrs
+    #  we forced it to nan
+    xr_dataset['sigma0_ocean'].attrs['comment'] = "clipped, no values <=0"
 
     xr_dataset['sigma0_ocean_raw'] = xr.where(xr_dataset['mask'], np.nan,
                                               xr_dataset['sigma0_raw'].compute()).transpose(*xr_dataset['sigma0_raw'].dims)
@@ -729,21 +704,24 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
         xr_dataset['sigma0_ocean_raw'] <= 0, np.nan, xr_dataset['sigma0_ocean_raw'])
     xr_dataset['sigma0_ocean_raw'].attrs = xr_dataset['sigma0_raw'].attrs
 
+    # processing
     if dual_pol:
-
         if config["apply_flattening"]:
             xr_dataset = xr_dataset.assign(nesz_cross_final=(
                 ['line', 'sample'], windspeed.nesz_flattening(xr_dataset.nesz.sel(pol=crosspol), xr_dataset.incidence)))
             xr_dataset['nesz_cross_final'].attrs[
                 "comment"] = 'nesz has been flattened using windspeed.nesz_flattening'
+
         else:
             xr_dataset = xr_dataset.assign(
                 nesz_cross_final=(['line', 'sample'], xr_dataset.nesz.sel(pol=crosspol).values))
             xr_dataset['nesz_cross_final'].attrs["comment"] = 'nesz has not been flattened'
-        # dsig
 
-        sigma0_ocean_cross = xr_dataset['sigma0_ocean'].sel(pol=crosspol)
+        xr_dataset.nesz_cross_final.attrs['units'] = 'm^2 / m^2'
+        xr_dataset.nesz_cross_final.attrs['long_name'] = 'Noise Equivalent SigmaNaught'
 
+        # dsig
+        sigma0_ocean_cross = xr_dataset['sigma0_ocean'].sel(pol=crosspol)
         xr_dataset["dsig_cross"] = windspeed.get_dsig(config["dsig_"+crosspol_gmf+"_NAME"], xr_dataset.incidence,
                                                       sigma0_ocean_cross, xr_dataset.nesz_cross_final)
 
@@ -753,38 +731,148 @@ def makeL2(filename, out_folder, config_path, overwrite=False, generateCSV=True,
         sigma0_ocean_cross = None
         dsig_cross = 0.1  # default value set in xsarsea
 
+    model_vv = config["GMF_"+copol_gmf+"_NAME"]
+    model_vh = config["GMF_"+crosspol_gmf+"_NAME"]
+
     wind_co, wind_dual, windspeed_cr = inverse(dual_pol,
-                                               inc=xr_dataset.incidence,
+                                               inc=xr_dataset['incidence'],
                                                sigma0=xr_dataset['sigma0_ocean'].sel(
                                                    pol=copol),
                                                sigma0_dual=sigma0_ocean_cross,
                                                ancillary_wind=xr_dataset['ancillary_wind'],
                                                dsig_cr=dsig_cross,
-                                               model_vv=config["GMF_" +
-                                                               copol_gmf+"_NAME"],
-                                               model_vh=config["GMF_"+crosspol_gmf+"_NAME"])
+                                               model_vv=model_vv,
+                                               model_vh=model_vh)
 
-    #  get windspeeds
+    # windspeed_co
     xr_dataset['windspeed_co'] = np.abs(wind_co)
-    xr_dataset['windspeed_dual'] = np.abs(wind_dual)
+    xr_dataset["windspeed_co"].attrs["units"] = "m.s⁻1"
+    xr_dataset["windspeed_co"].attrs["long_name"] = "Wind speed inverted from model %s (%s)" % (
+        model_vv, copol)
+    xr_dataset["windspeed_co"].attrs["standart_name"] = "wind_speed"
+    xr_dataset["windspeed_co"].attrs["model"] = wind_co.attrs["model"]
+    del xr_dataset["windspeed_co"].attrs['comment']
+
+    # winddir_co
+    xr_dataset['winddir_co'] = (
+        90 - (np.angle(-np.conj(wind_co), deg=True)) + xr_dataset.ground_heading) % 360
+    xr_dataset["winddir_co"].attrs["units"] = "degrees_north"
+    xr_dataset["winddir_co"].attrs["long_name"] = "Wind direction in meteorological convention, 0=North, 90=East, inverted from model %s (%s)" % (
+        model_vv, copol)
+    xr_dataset["winddir_co"].attrs["standart_name"] = "wind_direction"
+    xr_dataset["winddir_co"].attrs["model"] = wind_co.attrs["model"]
 
+    # windspeed_dual / windspeed_cr / /winddir_dual / winddir_cr
     if dual_pol:
+        xr_dataset['windspeed_dual'] = np.abs(wind_dual)
+        xr_dataset["windspeed_dual"].attrs["units"] = "m.s⁻1"
+        xr_dataset["windspeed_dual"].attrs["long_name"] = "Wind speed inverted from model %s (%s) & %s (%s)" % (
+            model_vv, copol, model_vh, crosspol)
+        xr_dataset["windspeed_dual"].attrs["standart_name"] = "wind_speed"
+        xr_dataset["windspeed_dual"].attrs["model"] = wind_dual.attrs["model"]
+        del xr_dataset["windspeed_dual"].attrs['comment']
+
+        xr_dataset['winddir_dual'] = (
+            90 - (np.angle(-np.conj(wind_dual), deg=True)) + xr_dataset.ground_heading) % 360
+        xr_dataset["winddir_dual"].attrs["units"] = "degrees_north"
+        xr_dataset["winddir_dual"].attrs["long_name"] = "Wind direction in meteorological convention, 0=North, 90=East inverted from model %s (%s) & %s (%s)" % (
+            model_vv, copol, model_vh, crosspol)
+        xr_dataset["winddir_dual"].attrs["standart_name"] = "wind_direction"
+        xr_dataset["winddir_dual"].attrs["model"] = wind_dual.attrs["model"]
+
         xr_dataset = xr_dataset.assign(
             windspeed_cross=(['line', 'sample'], windspeed_cr))
-    else:
-        xr_dataset['windspeed_cross'] = windspeed_cr
-
-    #  get winddirections
-    xr_dataset['winddir_co'] = (
-        90 - (np.angle(-np.conj(wind_co), deg=True)) + xr_dataset.ground_heading) % 360
-
-    xr_dataset['winddir_dual'] = (
-        90 - (np.angle(-np.conj(wind_dual), deg=True)) + xr_dataset.ground_heading) % 360
-    xr_dataset['winddir_cross'] = xr_dataset['winddir_dual'].copy()
+        xr_dataset["windspeed_cross"].attrs["units"] = "m.s⁻1"
+        xr_dataset["windspeed_cross"].attrs["long_name"] = "Wind Speed inverted from model %s (%s)" % (
+            model_vh, crosspol)
+        xr_dataset["windspeed_cross"].attrs["standart_name"] = "wind_speed"
+        xr_dataset["windspeed_cross"].attrs["model"] = "%s" % (model_vh)
+
+        xr_dataset['winddir_cross'] = xr_dataset['winddir_dual'].copy()
+        xr_dataset["winddir_cross"].attrs["units"] = "degrees_north"
+        xr_dataset["winddir_cross"].attrs["long_name"] = "Wind direction in meteorological convention, 0=North, 90=East, copied from dualpol"
+        xr_dataset["winddir_cross"].attrs["standart_name"] = "wind_direction"
+        xr_dataset["winddir_cross"].attrs["model"] = "No model used ; content is a copy of dualpol wind direction"
 
     xr_dataset, encoding = makeL2asOwi(
         xr_dataset, dual_pol, copol, crosspol, copol_gmf, crosspol_gmf, config)
 
+    #  add attributes
+    firstMeasurementTime = None
+    lastMeasurementTime = None
+    try:
+        firstMeasurementTime = datetime.datetime.strptime(xr_dataset.attrs['start_date'],
+                                                          "%Y-%m-%d %H:%M:%S.%f").strftime(
+            "%Y-%m-%dT%H:%M:%SZ")
+        lastMeasurementTime = datetime.datetime.strptime(xr_dataset.attrs['stop_date'],
+                                                         "%Y-%m-%d %H:%M:%S.%f").strftime(
+            "%Y-%m-%dT%H:%M:%SZ")
+    except:
+        firstMeasurementTime = datetime.datetime.strptime(xr_dataset.attrs['start_date'],
+                                                          "%Y-%m-%d %H:%M:%S").strftime(
+            "%Y-%m-%dT%H:%M:%SZ")
+        lastMeasurementTime = datetime.datetime.strptime(xr_dataset.attrs['stop_date'],
+                                                         "%Y-%m-%d %H:%M:%S").strftime(
+            "%Y-%m-%dT%H:%M:%SZ")
+
+    attrs = {
+        "TITLE": "Sentinel-1 OWI Component",
+        "productOwner": "IFREMER",
+        "sourceProduct": (xr_dataset.attrs["safe"] if "safe" in xr_dataset.attrs else os.path.basename(xr_dataset.attrs["product_path"])),
+        "sourceProduct_fullpath": xr_dataset.attrs.pop('name'),
+        "missionName": sensor_longname,
+        "missionPhase":  "Operational",
+        "polarisation": xr_dataset.attrs['pols'],
+        "acquisitionStation": '',
+        "xsar_version": xsar.__version__,
+        "xsarsea_version": xsarsea.__version__,
+        "pythonVersion": f"{sys.version_info.major}.{sys.version_info.minor}.{sys.version_info.micro}",
+        "polarisationRatio": "/",
+        "l2ProcessingUtcTime": datetime.datetime.now().strftime("%Y-%m-%dT%H:%M:%SZ"),
+        "processingCenter": "IFREMER",
+        "firstMeasurementTime": firstMeasurementTime,
+        "lastMeasurementTime": lastMeasurementTime,
+        "clmSource": "/",
+        "bathySource": "/",
+        "oswAlgorithmName": "grdwindinversion",
+        "owiAlgorithmVersion": grdwindinversion.__version__,
+        "gmf": config['GMF_'+copol_gmf+'_NAME'] + ", " + config["GMF_"+crosspol_gmf+"_NAME"],
+        "iceSource": "/",
+        "owiNoiseCorrection": "True",
+        "inversionTabGMF": config['GMF_'+copol_gmf +
+                                  '_NAME'] + ", " + config["GMF_"+crosspol_gmf+"_NAME"],
+        "wnf_3km_average": "False",
+        "owiWindSpeedSrc": "owiWindSpeed",
+        "owiWindDirectionSrc": "/",
+        "ancillary_source": xr_dataset.attrs['ancillary_source']
+    }
+
+    if ((recalibration) & ("SENTINEL" in sensor_longname)):
+        attrs["path_aux_pp1_new"] = os.path.basename(os.path.dirname(
+            os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_pp1_new'])))
+        attrs["path_aux_cal_new"] = os.path.basename(os.path.dirname(
+            os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_cal_new'])))
+
+        attrs["path_aux_pp1_old"] = os.path.basename(os.path.dirname(
+            os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_pp1_old'])))
+        attrs["path_aux_cal_old"] = os.path.basename(os.path.dirname(
+            os.path.dirname(xsar_dataset.datatree['recalibration'].attrs['path_aux_cal_old'])))
+
+    # new one to match convention
+    _S1_added_attrs = ["product", "ipf", "multi_dataset", "footprint",
+                       "coverage", "orbit_pass", "platform_heading"]
+    _RS2_added_attrs = ["passDirection", "swath", "footprint", "coverage"]
+    _RCM_added_attrs = ["swath", "footprint", "coverage", "productId",]
+
+    for sup_attr in _S1_added_attrs + _RS2_added_attrs + _RCM_added_attrs:
+        if sup_attr in xr_dataset.attrs:
+            attrs[sup_attr] = xr_dataset.attrs[sup_attr]
+    for var in ['footprint', 'multidataset']:
+        if var in attrs:
+            attrs[var] = str(attrs[var])
+
+    xr_dataset.attrs = attrs
+
     os.makedirs(os.path.dirname(out_file), exist_ok=True)
 
     xr_dataset.to_netcdf(out_file, mode="w", encoding=encoding)

{grdwindinversion-0.2.3.post9 → grdwindinversion-0.2.3.post11/grdwindinversion.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: grdwindinversion
-Version: 0.2.3.post9
+Version: 0.2.3.post11
 Summary: Package to perform Wind inversion from GRD Level-1 SAR images
 License: MIT
 Classifier: Development Status :: 2 - Pre-Alpha

grdwindinversion-0.2.3.post9/docs/examples/wind-inversion-from-grd.ipynb

@@ -1,134 +0,0 @@
-{
- "cells": [
-  {
-   "cell_type": "markdown",
-   "id": "e578297a-d2dc-4cb4-b016-1181649d227b",
-   "metadata": {},
-   "source": [
-    "## Example of Wind inversion"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "fa19ad8e-ebaa-430e-bd91-c21d517e7e88",
-   "metadata": {},
-   "source": [
-    "### define a product GRD Sentinel-1"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "ef483f15-4690-460b-b688-fa966bab8cd8",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "import xsar"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "e04fa1b0-eb49-4e45-af0f-0be798e26f95",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "# test on an dummy GRD IW product\n",
-    "input_file = xsar.get_test_file('S1A_IW_GRDH_1SDV_20170907T103020_20170907T103045_018268_01EB76_Z010.SAFE')\n"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "f4886386-b19a-4c43-a790-192f4866a692",
-   "metadata": {},
-   "source": [
-    "### perform wind inversion"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "d3fda3a2-2ee6-4401-be8f-c0fab8344a49",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "from grdwindinversion.inversion import makeL2 \n",
-    "\n"
-   ]
-  },
-  {
-   "cell_type": "markdown",
-   "id": "c8924754-34f5-4a27-a598-492022b3b3f5",
-   "metadata": {},
-   "source": [
-    "### define config file"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "f95ee426-0b46-4d39-ae90-241e8bbaa332",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "import grdwindinversion\n",
-    "import os\n",
-    "config_file_s1 = os.path.join(os.path.dirname(grdwindinversion.__file__),'config_S1.yaml')\n",
-    "grdwindinversion.__file__\n",
-    "print(os.path.exists(config_file_s1))"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "3cbe649d-44d7-4a7f-9041-137092a61e07",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "out_folder = '/tmp/test_windinversion'"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "27bc671b-d78c-4d70-b898-eac4093253e0",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "makeL2?"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "id": "b90c78bb-8643-4d64-a463-28fcb68d9b38",
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "out_file,outds = makeL2(input_file,out_folder, config_path=config_file_s1,overwrite=True)\n",
-    "outds"
-   ]
-  }
- ],
- "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.9.15"
-  }
- },
- "nbformat": 4,
- "nbformat_minor": 5
-}