Adapted to 2019 services

0f135d33 · Luís de Sousa · 750d0ba8 · 0f135d33 · 750d0ba8 · 0f135d33
Commit 0f135d33 authored 5 years ago by Luís de Sousa
--- a/WCS-getExtent.ipynb
+++ b/WCS-getExtent.ipynb
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "This time a connection is made to the service for soil pH:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from owslib.wcs import WebCoverageService\n",
+    "wcs = WebCoverageService('http://maps.isric.org/mapserv?map=/map/phh2o.map', version='1.0.0')"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "A bounding box broadly matching Senegal:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "bbox = (-1784000, 1356000, -1140000, 1863000)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Now fetch the coverage for Senegal:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "response = wcs.getCoverage(\n",
+    "    identifier='phh20_0-5cm_mean', \n",
+    "    crs='urn:ogc:def:crs:EPSG::152160',\n",
+    "    bbox=bbox, \n",
+    "    resx=1, resy=1, \n",
+    "    format='GeoTIFF')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "with open('./data/Senegal_pH_0-5_mean.tif', 'wb') as file:\n",
+    "    file.write(response.read())"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "With the data on the client side some regular interaction can be started with a library like resterIO:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 14,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "RasterioIOError",
+     "evalue": "'./data/Senegal_CEC0m.tif' not recognized as a supported file format.",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mCPLE_OpenFailedError\u001b[0m                      Traceback (most recent call last)",
+      "\u001b[0;32mrasterio/_base.pyx\u001b[0m in \u001b[0;36mrasterio._base.DatasetBase.__init__\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;32mrasterio/_shim.pyx\u001b[0m in \u001b[0;36mrasterio._shim.open_dataset\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;32mrasterio/_err.pyx\u001b[0m in \u001b[0;36mrasterio._err.exc_wrap_pointer\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;31mCPLE_OpenFailedError\u001b[0m: './data/Senegal_CEC0m.tif' not recognized as a supported file format.",
+      "\nDuring handling of the above exception, another exception occurred:\n",
+      "\u001b[0;31mRasterioIOError\u001b[0m                           Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-14-3c4dada614d1>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mrasterio\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mCEC\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mrasterio\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"./data/Senegal_CEC0m.tif\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"GTiff\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      3\u001b[0m \u001b[0mshow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mCEC\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtitle\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Cation Exchange Capacity at 0 metres in Senegal'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcmap\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'gist_ncar'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/usr/lib/python3/dist-packages/rasterio/env.py\u001b[0m in \u001b[0;36mwrapper\u001b[0;34m(*args, **kwds)\u001b[0m\n\u001b[1;32m    428\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    429\u001b[0m         \u001b[0;32mwith\u001b[0m \u001b[0menv_ctor\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0msession\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msession\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 430\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    431\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    432\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0mwrapper\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/usr/lib/python3/dist-packages/rasterio/__init__.py\u001b[0m in \u001b[0;36mopen\u001b[0;34m(fp, mode, driver, width, height, count, crs, transform, dtype, nodata, sharing, **kwargs)\u001b[0m\n\u001b[1;32m    214\u001b[0m         \u001b[0;31m# None.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    215\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mmode\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'r'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 216\u001b[0;31m             \u001b[0ms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mDatasetReader\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdriver\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msharing\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msharing\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    217\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0mmode\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'r+'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    218\u001b[0m             \u001b[0ms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_writer_for_path\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdriver\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msharing\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msharing\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32mrasterio/_base.pyx\u001b[0m in \u001b[0;36mrasterio._base.DatasetBase.__init__\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;31mRasterioIOError\u001b[0m: './data/Senegal_CEC0m.tif' not recognized as a supported file format."
+     ]
+    }
+   ],
+   "source": [
+    "import rasterio\n",
+    "CEC = rasterio.open(\"./data/Senegal_CEC0m.tif\", driver=\"GTiff\")\n",
+    "show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Next steps:\n",
+    "https://publicwiki.deltares.nl/display/OET/WCS+primer\n",
+    "https://geoscripting-wur.github.io/PythonRaster/"
+   ]
+  }
+ ],
+ "metadata": {
+  "@webio": {
+   "lastCommId": "C7CD903630FF4971A278510DC0E27447",
+   "lastKernelId": "5e40268c-f755-4e83-aec6-935c11d843f6"
+  },
+  "kernelspec": {
+   "display_name": "Python 3",
+   "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.6.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
+%% Cell type:markdown id: tags:
+
+This time a connection is made to the service for soil pH:
+
+%% Cell type:code id: tags:
+
+``` python
+from owslib.wcs import WebCoverageService
+wcs = WebCoverageService('http://maps.isric.org/mapserv?map=/map/phh2o.map', version='1.0.0')
+```
+
+%% Cell type:markdown id: tags:
+
+A bounding box broadly matching Senegal:
+
+%% Cell type:code id: tags:
+
+``` python
+bbox = (-1784000, 1356000, -1140000, 1863000)
+```
+
+%% Cell type:markdown id: tags:
+
+Now fetch the coverage for Senegal:
+
+%% Cell type:code id: tags:
+
+``` python
+response = wcs.getCoverage(
+    identifier='phh20_0-5cm_mean',
+    crs='urn:ogc:def:crs:EPSG::152160',
+    bbox=bbox,
+    resx=1, resy=1,
+    format='GeoTIFF')
+```
+
+%% Cell type:code id: tags:
+
+``` python
+with open('./data/Senegal_pH_0-5_mean.tif', 'wb') as file:
+    file.write(response.read())
+```
+
+%% Cell type:markdown id: tags:
+
+With the data on the client side some regular interaction can be started with a library like resterIO:
+
+%% Cell type:code id: tags:
+
+``` python
+import rasterio
+CEC = rasterio.open("./data/Senegal_CEC0m.tif", driver="GTiff")
+show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')
+```
+
+%% Output
+
+    ---------------------------------------------------------------------------
+    CPLE_OpenFailedError                      Traceback (most recent call last)
+    rasterio/_base.pyx in rasterio._base.DatasetBase.__init__()
+    rasterio/_shim.pyx in rasterio._shim.open_dataset()
+    rasterio/_err.pyx in rasterio._err.exc_wrap_pointer()
+    CPLE_OpenFailedError: './data/Senegal_CEC0m.tif' not recognized as a supported file format.
+
+During handling of the above exception, another exception occurred:
+    RasterioIOError                           Traceback (most recent call last)
+    <ipython-input-14-3c4dada614d1> in <module>()
+          1 import rasterio
+    ----> 2 CEC = rasterio.open("./data/Senegal_CEC0m.tif", driver="GTiff")
+          3 show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')
+    /usr/lib/python3/dist-packages/rasterio/env.py in wrapper(*args, **kwds)
+        428
+        429         with env_ctor(session=session):
+    --> 430             return f(*args, **kwds)
+        431
+        432     return wrapper
+    /usr/lib/python3/dist-packages/rasterio/__init__.py in open(fp, mode, driver, width, height, count, crs, transform, dtype, nodata, sharing, **kwargs)
+        214         # None.
+        215         if mode == 'r':
+    --> 216             s = DatasetReader(path, driver=driver, sharing=sharing, **kwargs)
+        217         elif mode == 'r+':
+        218             s = get_writer_for_path(path)(path, mode, driver=driver, sharing=sharing, **kwargs)
+    rasterio/_base.pyx in rasterio._base.DatasetBase.__init__()
+    RasterioIOError: './data/Senegal_CEC0m.tif' not recognized as a supported file format.
+
+%% Cell type:markdown id: tags:
+
+Next steps:
+https://publicwiki.deltares.nl/display/OET/WCS+primer
+https://geoscripting-wur.github.io/PythonRaster/
--- a/WCS.ipynb
+++ b/WCS.ipynb
--- a/WCS_basics.ipynb
+++ b/WCS_basics.ipynb
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "First load the `WebCoverageService` class from the OWSLib and create a connection to a service:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 16,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "from owslib.wcs import WebCoverageService\n",
+    "wcs = WebCoverageService('http://maps.isric.org/mapserv?map=/map/bdod.map', version='1.0.0')"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "A bounding box broadly matching Senegal:"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The absense of errors means that a connection was sucessfully established. To get more information about this service, start by identifying the operations available from this service:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "metadata": {
+    "scrolled": true
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "['GetCapabilities', 'DescribeCoverage', 'GetCoverage']\n"
+     ]
+    }
+   ],
+   "source": [
+    "print([op.name for op in wcs.operations])"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The full list of coverages available from this service is in the `contents` property:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "metadata": {
+    "scrolled": true
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "['bdod_15-30cm_Q0.95', 'bdod_0-5cm_mean', 'bdod_5-15cm_Q0.95', 'bdod_15-30cm_Q0.5', 'bdod_30-60cm_Q0.05', 'bdod_100-200cm_Q0.05', 'bdod_15-30cm_Q0.05', 'bdod_100-200cm_Q0.95', 'bdod_0-5cm_Q0.5', 'bdod_0-5cm_Q0.05', 'bdod_5-15cm_Q0.5', 'bdod_100-200cm_mean', 'bdod_60-100cm_Q0.95', 'bdod_5-15cm_Q0.05', 'bdod_30-60cm_mean', 'bdod_5-15cm_mean', 'bdod_100-200cm_Q0.5', 'bdod_30-60cm_Q0.95', 'bdod_15-30cm_mean', 'bdod_0-5cm_Q0.95', 'bdod_60-100cm_Q0.05', 'bdod_60-100cm_Q0.5', 'bdod_60-100cm_mean', 'bdod_30-60cm_Q0.5']\n"
+     ]
+    }
+   ],
+   "source": [
+    "print(list(wcs.contents))"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "That is a large set of coverages, but it is easy to filter the dictionary. For instance, to get the name of all  coverages for the 0 cm to 5 cm depth interval:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "['bdod_0-5cm_mean', 'bdod_0-5cm_Q0.5', 'bdod_0-5cm_Q0.05', 'bdod_0-5cm_Q0.95']\n"
+     ]
+    }
+   ],
+   "source": [
+    "names = [k for k in wcs.contents.keys() if k.startswith(\"bdod_0-5cm\")]\n",
+    "print(names)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Or to search for all the coverages reporting the median prediction:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 10,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "['bdod_15-30cm_Q0.5', 'bdod_0-5cm_Q0.5', 'bdod_5-15cm_Q0.5', 'bdod_100-200cm_Q0.5', 'bdod_60-100cm_Q0.5', 'bdod_30-60cm_Q0.5']\n"
+     ]
+    }
+   ],
+   "source": [
+    "q0_5_covs = [k for k in wcs.contents.keys() if k.find(\"Q0.5\") != -1]\n",
+    "print(q0_5_covs)"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "These are the SoilGrids predictions for bulk density for the six standard depths defined in the GlobalSoilMap specifications.\n",
+    "\n",
+    "The details for one of these coverages can be inspected using the identifiers above:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 14,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "[urn:ogc:def:crs:EPSG::152160, urn:ogc:def:crs:EPSG::152160]"
+      ]
+     },
+     "execution_count": 14,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "bdod_5_15_median = wcs.contents['bdod_5-15cm_Q0.5']\n",
+    "bdod_5_15_median.supportedCRS"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 15,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "['GEOTIFF_16']"
+      ]
+     },
+     "execution_count": 15,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "bdod_5_15_median.supportedFormats"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "bbox = (-17.1, 12.65, -13.15, 14.95) "
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Now fetch the coverage for Senegal:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 17,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "response = wcs.getCoverage(\n",
+    "    identifier='phh20_0-5cm_mean', \n",
+    "    crs='urn:ogc:def:crs:EPSG::152160',\n",
+    "    bbox=bbox, \n",
+    "    resx=1, resy=1, \n",
+    "    format='GeoTIFF')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 18,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "with open('./data/Senegal_CEC0m.tif', 'wb') as file:\n",
+    "    file.write(response.read())"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "With the data on the client side some regular interaction can be started with a library like resterIO:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 14,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "RasterioIOError",
+     "evalue": "'./data/Senegal_CEC0m.tif' not recognized as a supported file format.",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mCPLE_OpenFailedError\u001b[0m                      Traceback (most recent call last)",
+      "\u001b[0;32mrasterio/_base.pyx\u001b[0m in \u001b[0;36mrasterio._base.DatasetBase.__init__\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;32mrasterio/_shim.pyx\u001b[0m in \u001b[0;36mrasterio._shim.open_dataset\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;32mrasterio/_err.pyx\u001b[0m in \u001b[0;36mrasterio._err.exc_wrap_pointer\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;31mCPLE_OpenFailedError\u001b[0m: './data/Senegal_CEC0m.tif' not recognized as a supported file format.",
+      "\nDuring handling of the above exception, another exception occurred:\n",
+      "\u001b[0;31mRasterioIOError\u001b[0m                           Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-14-3c4dada614d1>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mrasterio\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mCEC\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mrasterio\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"./data/Senegal_CEC0m.tif\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"GTiff\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      3\u001b[0m \u001b[0mshow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mCEC\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtitle\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Cation Exchange Capacity at 0 metres in Senegal'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcmap\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'gist_ncar'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/usr/lib/python3/dist-packages/rasterio/env.py\u001b[0m in \u001b[0;36mwrapper\u001b[0;34m(*args, **kwds)\u001b[0m\n\u001b[1;32m    428\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    429\u001b[0m         \u001b[0;32mwith\u001b[0m \u001b[0menv_ctor\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0msession\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msession\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 430\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    431\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    432\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0mwrapper\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/usr/lib/python3/dist-packages/rasterio/__init__.py\u001b[0m in \u001b[0;36mopen\u001b[0;34m(fp, mode, driver, width, height, count, crs, transform, dtype, nodata, sharing, **kwargs)\u001b[0m\n\u001b[1;32m    214\u001b[0m         \u001b[0;31m# None.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    215\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mmode\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'r'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 216\u001b[0;31m             \u001b[0ms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mDatasetReader\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdriver\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msharing\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msharing\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    217\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0mmode\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'r+'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    218\u001b[0m             \u001b[0ms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mget_writer_for_path\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdriver\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mdriver\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msharing\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0msharing\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32mrasterio/_base.pyx\u001b[0m in \u001b[0;36mrasterio._base.DatasetBase.__init__\u001b[0;34m()\u001b[0m\n",
+      "\u001b[0;31mRasterioIOError\u001b[0m: './data/Senegal_CEC0m.tif' not recognized as a supported file format."
+     ]
+    }
+   ],
+   "source": [
+    "import rasterio\n",
+    "CEC = rasterio.open(\"./data/Senegal_CEC0m.tif\", driver=\"GTiff\")\n",
+    "show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Next steps:\n",
+    "https://publicwiki.deltares.nl/display/OET/WCS+primer\n",
+    "https://geoscripting-wur.github.io/PythonRaster/"
+   ]
+  }
+ ],
+ "metadata": {
+  "@webio": {
+   "lastCommId": "15E983DED5BA45E28DEE0958687E8D50",
+   "lastKernelId": "9ef91650-4d26-43d0-a331-f95652b81618"
+  },
+  "kernelspec": {
+   "display_name": "Python 3",
+   "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.6.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
+%% Cell type:markdown id: tags:
+
+First load the `WebCoverageService` class from the OWSLib and create a connection to a service:
+
+%% Cell type:code id: tags:
+
+``` python
+from owslib.wcs import WebCoverageService
+wcs = WebCoverageService('http://maps.isric.org/mapserv?map=/map/bdod.map', version='1.0.0')
+```
+
+%% Cell type:markdown id: tags:
+
+A bounding box broadly matching Senegal:
+
+%% Cell type:markdown id: tags:
+
+The absense of errors means that a connection was sucessfully established. To get more information about this service, start by identifying the operations available from this service:
+
+%% Cell type:code id: tags:
+
+``` python
+print([op.name for op in wcs.operations])
+```
+
+%% Output
+
+    ['GetCapabilities', 'DescribeCoverage', 'GetCoverage']
+
+%% Cell type:markdown id: tags:
+
+The full list of coverages available from this service is in the `contents` property:
+
+%% Cell type:code id: tags:
+
+``` python
+print(list(wcs.contents))
+```
+
+%% Output
+
+    ['bdod_15-30cm_Q0.95', 'bdod_0-5cm_mean', 'bdod_5-15cm_Q0.95', 'bdod_15-30cm_Q0.5', 'bdod_30-60cm_Q0.05', 'bdod_100-200cm_Q0.05', 'bdod_15-30cm_Q0.05', 'bdod_100-200cm_Q0.95', 'bdod_0-5cm_Q0.5', 'bdod_0-5cm_Q0.05', 'bdod_5-15cm_Q0.5', 'bdod_100-200cm_mean', 'bdod_60-100cm_Q0.95', 'bdod_5-15cm_Q0.05', 'bdod_30-60cm_mean', 'bdod_5-15cm_mean', 'bdod_100-200cm_Q0.5', 'bdod_30-60cm_Q0.95', 'bdod_15-30cm_mean', 'bdod_0-5cm_Q0.95', 'bdod_60-100cm_Q0.05', 'bdod_60-100cm_Q0.5', 'bdod_60-100cm_mean', 'bdod_30-60cm_Q0.5']
+
+%% Cell type:markdown id: tags:
+
+That is a large set of coverages, but it is easy to filter the dictionary. For instance, to get the name of all  coverages for the 0 cm to 5 cm depth interval:
+
+%% Cell type:code id: tags:
+
+``` python
+names = [k for k in wcs.contents.keys() if k.startswith("bdod_0-5cm")]
+print(names)
+```
+
+%% Output
+
+    ['bdod_0-5cm_mean', 'bdod_0-5cm_Q0.5', 'bdod_0-5cm_Q0.05', 'bdod_0-5cm_Q0.95']
+
+%% Cell type:markdown id: tags:
+
+Or to search for all the coverages reporting the median prediction:
+
+%% Cell type:code id: tags:
+
+``` python
+q0_5_covs = [k for k in wcs.contents.keys() if k.find("Q0.5") != -1]
+print(q0_5_covs)
+```
+
+%% Output
+
+    ['bdod_15-30cm_Q0.5', 'bdod_0-5cm_Q0.5', 'bdod_5-15cm_Q0.5', 'bdod_100-200cm_Q0.5', 'bdod_60-100cm_Q0.5', 'bdod_30-60cm_Q0.5']
+
+%% Cell type:markdown id: tags:
+
+These are the SoilGrids predictions for bulk density for the six standard depths defined in the GlobalSoilMap specifications.
+
+The details for one of these coverages can be inspected using the identifiers above:
+
+%% Cell type:code id: tags:
+
+``` python
+bdod_5_15_median = wcs.contents['bdod_5-15cm_Q0.5']
+bdod_5_15_median.supportedCRS
+```
+
+%% Output
+
+    [urn:ogc:def:crs:EPSG::152160, urn:ogc:def:crs:EPSG::152160]
+
+%% Cell type:code id: tags:
+
+``` python
+bdod_5_15_median.supportedFormats
+```
+
+%% Output
+
+    ['GEOTIFF_16']
+
+%% Cell type:code id: tags:
+
+``` python
+bbox = (-17.1, 12.65, -13.15, 14.95)
+```
+
+%% Cell type:markdown id: tags:
+
+Now fetch the coverage for Senegal:
+
+%% Cell type:code id: tags:
+
+``` python
+response = wcs.getCoverage(
+    identifier='phh20_0-5cm_mean',
+    crs='urn:ogc:def:crs:EPSG::152160',
+    bbox=bbox,
+    resx=1, resy=1,
+    format='GeoTIFF')
+```
+
+%% Cell type:code id: tags:
+
+``` python
+with open('./data/Senegal_CEC0m.tif', 'wb') as file:
+    file.write(response.read())
+```
+
+%% Cell type:markdown id: tags:
+
+With the data on the client side some regular interaction can be started with a library like resterIO:
+
+%% Cell type:code id: tags:
+
+``` python
+import rasterio
+CEC = rasterio.open("./data/Senegal_CEC0m.tif", driver="GTiff")
+show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')
+```
+
+%% Output
+
+    ---------------------------------------------------------------------------
+    CPLE_OpenFailedError                      Traceback (most recent call last)
+    rasterio/_base.pyx in rasterio._base.DatasetBase.__init__()
+    rasterio/_shim.pyx in rasterio._shim.open_dataset()
+    rasterio/_err.pyx in rasterio._err.exc_wrap_pointer()
+    CPLE_OpenFailedError: './data/Senegal_CEC0m.tif' not recognized as a supported file format.
+
+During handling of the above exception, another exception occurred:
+    RasterioIOError                           Traceback (most recent call last)
+    <ipython-input-14-3c4dada614d1> in <module>()
+          1 import rasterio
+    ----> 2 CEC = rasterio.open("./data/Senegal_CEC0m.tif", driver="GTiff")
+          3 show(CEC, title='Cation Exchange Capacity at 0 metres in Senegal', cmap='gist_ncar')
+    /usr/lib/python3/dist-packages/rasterio/env.py in wrapper(*args, **kwds)
+        428
+        429         with env_ctor(session=session):
+    --> 430             return f(*args, **kwds)
+        431
+        432     return wrapper
+    /usr/lib/python3/dist-packages/rasterio/__init__.py in open(fp, mode, driver, width, height, count, crs, transform, dtype, nodata, sharing, **kwargs)
+        214         # None.
+        215         if mode == 'r':
+    --> 216             s = DatasetReader(path, driver=driver, sharing=sharing, **kwargs)
+        217         elif mode == 'r+':
+        218             s = get_writer_for_path(path)(path, mode, driver=driver, sharing=sharing, **kwargs)
+    rasterio/_base.pyx in rasterio._base.DatasetBase.__init__()
+    RasterioIOError: './data/Senegal_CEC0m.tif' not recognized as a supported file format.
+
+%% Cell type:markdown id: tags:
+
+Next steps:
+https://publicwiki.deltares.nl/display/OET/WCS+primer
+https://geoscripting-wur.github.io/PythonRaster/