Skip to content

explain _repr_svg and print wms/wfs responses directly #240

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
tomkralidis merged 1 commit into from
Jun 21, 2026
+73 −37
Merged

explain _repr_svg and print wms/wfs responses directly #240

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
28 changes: 26 additions & 2 deletions workshop/jupyter/content/notebooks/02-geometry.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
"source": [
"# Geometry\n",
"\n",
"This lesson brings you a bit back to highschool math: (Euclidian) Geometry, \n",
"This lesson brings you a bit back to highschool math: ([Euclidian](https://en.wikipedia.org/wiki/Euclidean_geometry)) Geometry, \n",
"hopefully with more fun as within (geo)spatial IT/GIS we usually deal with\n",
"real-world objects (roads, lakes, forests, etc) abstracted as geometries.\n",
"\n",
Expand Down Expand Up @@ -204,7 +204,31 @@
"source": [
"from shapely.geometry import Polygon\n",
"polygon = Polygon([(0, 0), (3, 0), (3, 4)])\n",
"polygon"
"polygon",
""
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Jupyter (IPython) benefits from Shapely's `geometry._repr_svg_()` method to render any shapely geometry as image.\n",
"Whenever you `log` the value of a shapely geometry, an image (svg) is rendered of the geometry."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"# Explicit use of _repr_svg_() method\n",
"svg = polygon._repr_svg_()\n",
"print(svg)"
]
},
{
Expand Down
12 changes: 8 additions & 4 deletions workshop/jupyter/content/notebooks/09-publishing.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
"source": [
"# Publishing\n",
"\n",
"**NB This Section is optional and not included in the FOSS4G 2024 Workshop**\n",
"**NB This Section is optional and not included in the FOSS4G 2026 Workshop**\n",
"\n",
"It may still be valuable to follow if time left or at home. Mind that the service containers can only be run in a local (Docker)\n",
"environment, i.e. not through the Workshop Binder version in the cloud.\n",
Expand Down Expand Up @@ -58,8 +58,7 @@
"\n",
"Here you will see 10 feature collections listed on the resulting webpage. Feature collections are\n",
"identified by the `\"itemType\": \"feature\"` in the collection definition in the JSON response.\n",
"\n",
"Now let's add the WOUDC station data to our pygeoapi instance.\n"
"\n"
]
},
{
Expand Down Expand Up @@ -116,7 +115,12 @@
"# Get items (paged) in Lakes collection\n",
"lakes = oa_feat.collection('lakes')\n",
"lakes_query = oa_feat.collection_items('lakes')\n",
"lakes_query['features'][0]"
"# Plot first lake\n",
"lake = lakes_query['features'][0]\n",
"print(lake.get('properties').get('name'))\n",
"from shapely.geometry import shape\n",
"lakepolygon = shape(lake.get('geometry'))\n",
"lakepolygon"
]
},
{
Expand Down
70 changes: 39 additions & 31 deletions workshop/jupyter/content/notebooks/10-remote-data.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@
"`OWSLib` enables you to connect to these services as \n",
"a client, mainly to fetch and query data from them. \n",
"Currently WMS, WFS, WCS, CSW, WPS, SOS, WMC and the more recent OGC APIs like \n",
"*OGC API - Features* (formerly called \"WFS version 3\") are supported.\n",
"*OGC API - Features* are supported.\n",
"The list of supported services is growing. \n",
"\n",
"Documentation: https://owslib.readthedocs.io\n"
Expand Down Expand Up @@ -150,7 +150,14 @@
"# Get items (paged) in Lakes collection\n",
"lakes = oa_feat.collection('lakes')\n",
"lakes_query = oa_feat.collection_items('lakes')\n",
"lakes_query['features'][0]"
"lakes_query['features'][0]\n",
"\n",
"# Plot first lake\n",
"lake = lakes_query['features'][0]\n",
"print(lake.get('properties').get('name'))\n",
"from shapely.geometry import shape\n",
"lakepolygon = shape(lake.get('geometry'))\n",
"lakepolygon"
]
},
{
Expand Down Expand Up @@ -258,7 +265,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"Download and save the data (image)"
"Fetch and display the image via WMS"
]
},
{
Expand All @@ -278,30 +285,8 @@
" bbox=[1.0, 50.0, 10.0, 54.0],\n",
" format=\"image/png\")\n",
"\n",
"save_fp = 'test/10-wms.png'\n",
"with open(save_fp, 'wb') as out:\n",
"\tout.write(img.read())"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Display image in Notebook or view in browser via [this link](test/10-wms.png)."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"pycharm": {
"name": "#%%\n"
}
},
"outputs": [],
"source": [
"from IPython.display import Image\n",
"Image(filename=save_fp)"
"from IPython.display import Image, display\n",
"display(Image(data=img.read()))\n"
]
},
{
Expand Down Expand Up @@ -500,7 +485,20 @@
},
"outputs": [],
"source": [
"response = wfs.getfeature(typename='ms:cities', bbox=(15.7500260759, 42.65, 19.59976, 45.2337767604))"
"response = wfs.getfeature(typename='ms:cities', bbox=(15.7500260759, 42.65, 19.59976, 45.2337767604))\n",
"data = response.read()\n",
"\n",
"# Parse the GML response as Shapely geometries\n",
"from io import BytesIO\n",
"from geopandas import read_file\n",
"gdf = read_file(BytesIO(data))\n",
"\n",
"# Combine the points in a MultiPoint\n",
"from shapely import MultiPoint\n",
"mpgeom = MultiPoint(list(gdf.geometry))\n",
"\n",
"# Plot the geometry\n",
"mpgeom"
]
},
{
Expand Down Expand Up @@ -888,7 +886,7 @@
"tags": []
},
"source": [
"Download data from selected layer and write to a file."
"Download data from selected layer and plot."
]
},
{
Expand All @@ -908,9 +906,19 @@
"source": [
"identifier = 'ge:boreholes_cores_rbins'\n",
"features = boreholes_wfs.getfeature([identifier])\n",
"data = features.read()\n",
"\n",
"# Parse the GML response as Shapely geometries\n",
"from io import BytesIO\n",
"from geopandas import read_file\n",
"gdf = read_file(BytesIO(data))\n",
"\n",
"# Combine the points in a MultiPoint\n",
"from shapely import MultiPoint\n",
"mpgeom = MultiPoint(list(gdf.geometry))\n",
"\n",
"with open('test/10-boreholes.gml', 'w', encoding='UTF-8') as out:\n",
" out.write(str(features.read(), 'UTF-8'))"
"# Plot the geometry\n",
"mpgeom"
]
},
{
Expand Down