Add stac demo

Author: arbennett

coding_modules/module_7_cloud_datasets/stac_demo.ipynb

@@ -0,0 +1,177 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import xarray as xr\n",
+    "import pystac_client\n",
+    "import planetary_computer\n",
+    "\n",
+    "from dask.distributed import Client, LocalCluster"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "cluster = LocalCluster(n_workers=4, threads_per_worker=1)\n",
+    "client = Client(cluster)\n",
+    "client"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "path = 'gs://gcp-public-data-arco-era5/ar/1959-2022-6h-1440x721.zarr'\n",
+    "era5_ds = xr.open_dataset(path, engine='zarr').chunk('auto')\n",
+    "era5_temperature = era5_ds['2m_temperature']\n",
+    "era5_precipitation = era5_ds['total_precipitation_6hr']"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "era5_precipitation"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "catalog = pystac_client.Client.open(\n",
+    "    \"https://planetarycomputer.microsoft.com/api/stac/v1/\",\n",
+    "    modifier=planetary_computer.sign_inplace,\n",
+    ")\n",
+    "collection = catalog.get_collection(\"cil-gdpcir-cc-by\")\n",
+    "\n",
+    "summary = collection.summaries.to_dict()\n",
+    "available_keys = list(summary.keys())\n",
+    "print(available_keys)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## What physical variables are stored in the dataset?"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "summary['cmip6:variable']"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## What GCMs are used in the dataset?"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "summary['cmip6:source_id']"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "##  What are the SSPs in the dataset?"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "summary['cmip6:experiment_id']"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "## How do we use this information to actually get the data?"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "search = catalog.search(\n",
+    "    collections=['cil-gdpcir-cc-by'],\n",
+    "    query={'cmip6:experiment_id': {'eq': 'historical'}}\n",
+    ")\n",
+    "historic_data = search.item_collection()\n",
+    "historic_data"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "search = catalog.search(\n",
+    "    collections=['cil-gdpcir-cc-by'],\n",
+    "    query={'cmip6:source_id': {'eq': 'GFDL-ESM4'}}\n",
+    ")\n",
+    "gfdl_data = search.item_collection()\n",
+    "gfdl_data"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "init",
+   "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.13"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}