Corrected format in getCoverage example

7fb0285d · Luís de Sousa · a33b5699 · 7fb0285d
Commit 7fb0285d authored 4 years ago by Luís de Sousa
--- a/02-WCS-getExtent.ipynb
+++ b/02-WCS-getExtent.ipynb
@@ -14,7 +14,7 @@
  },
  {
   "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
@@ -31,7 +31,7 @@
  },
  {
   "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
@@ -47,16 +47,29 @@
  },
  {
   "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 3,
   "metadata": {},
-   "outputs": [],
+   "outputs": [
+    {
+     "ename": "AttributeError",
+     "evalue": "'WebCoverageService_1_0_0' object has no attribute 'supportedFormats'",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-3-c0e26f5657aa>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      4\u001b[0m     \u001b[0mbbox\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mbbox\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      5\u001b[0m     \u001b[0mresx\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m250\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mresy\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m250\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 6\u001b[0;31m     \u001b[0mformat\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mwcs\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msupportedFormats\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;31m# GeoTiff 16 bit integer\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      7\u001b[0m )\n",
+      "\u001b[0;31mAttributeError\u001b[0m: 'WebCoverageService_1_0_0' object has no attribute 'supportedFormats'"
+     ]
+    }
+   ],
   "source": [
    "response = wcs.getCoverage(\n",
    "    identifier='phh2o_0-5cm_mean', \n",
    "    crs='urn:ogc:def:crs:EPSG::152160',\n",
    "    bbox=bbox, \n",
    "    resx=250, resy=250, \n",
-    "    format='GEOTIFF_16')"
+    "    format=''\n",
+    ")"
   ]
  },
  {
@@ -90,7 +103,7 @@
   "outputs": [],
   "source": [
    "import rasterio\n",
-    "CEC = rasterio.open(\"./data/Senegal_pH_0-5_mean.tif\", driver=\"GTiff\")"
+    "ph = rasterio.open(\"./data/Senegal_pH_0-5_mean.tif\", driver=\"GTiff\")"
   ]
  },
  {
@@ -108,8 +121,15 @@
   "source": [
    "from rasterio import plot\n",
    "%matplotlib inline\n",
-    "plot.show(CEC, title='Mean pH between 0 and 5 cm deep in Senegal', cmap='gist_ncar')"
+    "plot.show(ph, title='Mean pH between 0 and 5 cm deep in Senegal', cmap='gist_ncar')"
   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
  }
 ],
 "metadata": {

 %% Cell type:markdown id: tags:

 Obtaining a map segment from a WCS
 ====================

 Obtaining a map segment is another operation greatly simplified by the [OWSLib](https://geopython.github.io/OWSLib/) package. This case is conducted from a pre-established bounding box.

 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:

 The `getCoverage` method can now be used to fetch the map segment within the bounding box. Note the other parameters, [Section 1](01-WCS-basics.ipynb) showed how to obtain them.

 %% Cell type:code id: tags:

 ``` python
 response = wcs.getCoverage(
    identifier='phh2o_0-5cm_mean',
    crs='urn:ogc:def:crs:EPSG::152160',
    bbox=bbox,
    resx=250, resy=250,
-    format='GEOTIFF_16')
+    format=''
+)
 ```

+%% Output
+
+    ---------------------------------------------------------------------------
+    AttributeError                            Traceback (most recent call last)
+    <ipython-input-3-c0e26f5657aa> in <module>()
+          4     bbox=bbox,
+          5     resx=250, resy=250,
+    ----> 6     format=wcs.supportedFormats[0] # GeoTiff 16 bit integer
+          7 )
+    AttributeError: 'WebCoverageService_1_0_0' object has no attribute 'supportedFormats'
+
 %% Cell type:markdown id: tags:

 Now fetch the coverage for Senegal and save it to disk:

 %% 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 [rasterIO](https://rasterio.readthedocs.io/). First open the file from disk:

 %% Cell type:code id: tags:

 ``` python
 import rasterio
-CEC = rasterio.open("./data/Senegal_pH_0-5_mean.tif", driver="GTiff")
+ph = rasterio.open("./data/Senegal_pH_0-5_mean.tif", driver="GTiff")
 ```

 %% Cell type:markdown id: tags:

 Finally, use the `plot` class to plot the map:

 %% Cell type:code id: tags:

 ``` python
 from rasterio import plot
 %matplotlib inline
-plot.show(CEC, title='Mean pH between 0 and 5 cm deep in Senegal', cmap='gist_ncar')
+plot.show(ph, title='Mean pH between 0 and 5 cm deep in Senegal', cmap='gist_ncar')
+```
+
+%% Cell type:code id: tags:
+
+``` python
 ```