Quarto output

Author: MartinSchobben

notebooks/courses/microwave-remote-sensing.ipynb

@@ -8,7 +8,7 @@
     "# Microwave Remote Sensing\n",
     "\n",
     "\n",
-    "This course at the TU Wien teaches course, students are able to read, visualize and analyze Synthetic Aperture Radar (SAR) data. This will aid interpretation of SAR data based upon a physical understanding of sensing principles and the interaction of microwaves with natural objects.\n",
+    "This course at the TU Wien teaches students to read, visualize and analyze Synthetic Aperture Radar (SAR) data. This will aid interpretation of SAR data based upon a physical understanding of sensing principles and the interaction of microwaves with natural objects.\n",
     "\n",
     ":::{note}\n",
     "These notebooks contain interactive elements. The full interactive elements can only be viewed on Binder by clicking on the Binder badge or 🚀 button.\n",

notebooks/courses/microwave-remote-sensing.yml

@@ -7,25 +7,20 @@ dependencies:
   - cmcrameri
   - dask
   - datashader
-  - eomaps
   - folium
   - graphviz
   - holoviews
   - huggingface_hub
   - hvplot
   - intake
   - intake-xarray
-  - ipywidgets
-  - ipympl
   - jupyter
   - jupyter_bokeh
   - mamba
   - matplotlib
   - nodejs
-  - netcdf4
   - numpy
   - odc-stac
-  - pip
   - pyproj
   - pystac-client
   - python=3.11
@@ -35,8 +30,5 @@ dependencies:
   - scipy
   - seaborn
   - snaphu
-  - stackstac
   - xarray
-  - zarr
-  - pip:
-    - graphviz
+  - zarr=2.18.4

notebooks/courses/microwave-remote-sensing/02_in_class_exercise.ipynb

@@ -65,7 +65,7 @@
     "lonmin, lonmax = 16, 17\n",
     "bounds = (lonmin, latmin, lonmax, latmax)\n",
     "\n",
-    "time_range = \"2022-01-01/2022-12-30\"\n",
+    "time_range = \"2022-07-01/2022-07-31\"\n",
     "\n",
     "items = (\n",
     "    pystac_client.Client.open(\"https://stac.eodc.eu/api/v1\")\n",

notebooks/courses/microwave-remote-sensing/03_in_class_exercise.ipynb

@@ -50,7 +50,7 @@
    "source": [
     "url = \"https://huggingface.co/datasets/martinschobben/microwave-remote-sensing/resolve/main/microwave-remote-sensing.yml\"\n",
     "cat = intake.open_catalog(url)\n",
-    "gtc_dc = cat.gtc.read()\n",
+    "gtc_dc = cat.gtc.read().compute()\n",
     "gtc_dc"
    ]
   },
@@ -157,7 +157,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "coef_dc = cat.coef.read()\n",
+    "coef_dc = cat.coef.read().compute()\n",
     "coef_dc.hvplot.image(x=\"x\", y=\"y\", robust=True, data_aspect=1, cmap=\"Greys_r\",\n",
     "                     groupby=\"band\", rasterize=True).\\\n",
     "                        opts(frame_height=600, framewise=False, aspect=\"equal\")"
@@ -223,7 +223,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "rtc_dc = cat.rtc.read()\n",
+    "rtc_dc = cat.rtc.read().compute()\n",
     "\n",
     "# add to existing cube\n",
     "rtc_dc = xr.concat([coef_dc, rtc_dc], dim=\"band\")\n",

notebooks/courses/microwave-remote-sensing/08_in_class_exercise.ipynb

@@ -96,7 +96,7 @@
     "val_range = dict(vmin=0, vmax=255, cmap=\"gray\")\n",
     "\n",
     "for i, ds in enumerate(datasets):\n",
-    "    im = abs(ds.vv).plot(ax=ax[i], add_colorbar=False, **val_range)\n",
+    "    im = ds.intensity.plot(ax=ax[i], add_colorbar=False, **val_range)\n",
     "    ax[i].tick_params(axis=\"both\", which=\"major\")\n",
     "\n",
     "cbar = fig.colorbar(im, ax=ax, orientation=\"horizontal\", shrink=0.9, pad=0.2)\n",
@@ -120,18 +120,14 @@
    "outputs": [],
    "source": [
     "# Compute the intensity and phase from complex data\n",
-    "intensity = abs(iw2_ds.vv)\n",
-    "phase = xr.apply_ufunc(\n",
-    "    np.angle, iw2_ds.vv, dask=\"parallelized\", output_dtypes=[float]\n",
-    ")\n",
     "cmap_hls = sns.hls_palette(as_cmap=True)\n",
     "\n",
     "fig, axes = plt.subplots(1, 2, figsize=(15, 6))\n",
     "\n",
-    "intensity.plot(ax=axes[0], cmap=\"gray\", robust=True)\n",
+    "ds.intensity.plot(ax=axes[0], cmap=\"gray\", robust=True)\n",
     "axes[0].set_title(\"Intensity Measurement of IW2\")\n",
     "\n",
-    "phase.plot(ax=axes[1], cmap=cmap_hls)\n",
+    "ds.phase.plot(ax=axes[1], cmap=cmap_hls)\n",
     "axes[1].set_title(\"Phase Measurement of IW2\")\n",
     "\n",
     "plt.tight_layout()"