extract_geometry_wms.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "77c03259a4e541438b54625a2d716bdb",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HTML(value='\\n    <style>\\n        .output_scroll {\\n            height: unset !important;\\n            border…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import requests\n",
    "import mercantile\n",
    "from io import StringIO\n",
    "from urllib.parse import urlparse, parse_qsl, urlencode, urlunparse\n",
    "from affine import Affine\n",
    "from shapely.geometry import shape, box, mapping, Polygon\n",
    "from shapely.ops import cascaded_union\n",
    "import xml.etree.ElementTree as ET\n",
    "from rasterio.io import MemoryFile\n",
    "from rasterio import features\n",
    "import numpy as np\n",
    "from rasterio.features import sieve\n",
    "from shapely.ops import transform\n",
    "from pyproj import Transformer, CRS\n",
    "import json\n",
    "import matplotlib.pyplot as plt\n",
    "from descartes.patch import PolygonPatch\n",
    "from scipy.ndimage.filters import gaussian_filter\n",
    "from ipywidgets import interact, interactive, FloatSlider, interact_manual, HTML\n",
    "import warnings\n",
    "warnings.simplefilter(\"ignore\")\n",
    "\n",
    "\n",
    "style = \"\"\"\n",
    "    <style>\n",
    "        .output_scroll {\n",
    "            height: unset !important;\n",
    "            border-radius: unset !important;\n",
    "            -webkit-box-shadow: unset !important;\n",
    "            box-shadow: unset !important;\n",
    "        }\n",
    "    </style>\n",
    "    \"\"\"\n",
    "display(HTML(style))\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Enter URL: E.g. https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox}\n",
      "\n",
      "https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox}\n",
      "\n",
      "Enter Zoom level: E.g. 10\n",
      "\n",
      "10\n",
      "Send request to: https://wms.zh.ch/OGDOrthoZH?service=WMS&request=GetCapabilities&version=1.3.0\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.754097979680026,8.0859375,47.98992166741417,8.4375\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.51720069783939,8.0859375,47.754097979680026,8.4375\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.27922900257082,8.0859375,47.51720069783939,8.4375\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.04018214480666,8.0859375,47.27922900257082,8.4375\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.754097979680026,8.4375,47.98992166741417,8.7890625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.51720069783939,8.4375,47.754097979680026,8.7890625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.27922900257082,8.4375,47.51720069783939,8.7890625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.04018214480666,8.4375,47.27922900257082,8.7890625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.754097979680026,8.7890625,47.98992166741417,9.140625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.51720069783939,8.7890625,47.754097979680026,9.140625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.27922900257082,8.7890625,47.51720069783939,9.140625\n",
      "Send request to https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS=EPSG:4326&WIDTH=512&HEIGHT=512&BBOX=47.04018214480666,8.7890625,47.27922900257082,9.140625\n",
      "Extracted the following geometry:\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD5CAYAAADP2jUWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAuDklEQVR4nO3deZAe93ng9+/T13u/c2IGF0HwBG+CEqyVxKIsU7QsUbSSVVQSlVJcciXFSrK1JcfZaKOtyu4fKVfiSpTSJpujWNraXZcSmS5Z8q5lidJGMq1IlkgDAihRBAgQJAhggMHcxzvvvFf3kz/6HXAIDjDvHJh+j+dTHM68091vPz0zePrXv1NUFWOMMb3FSToAY4wxO8+SvzHG9CBL/sYY04Ms+RtjTA+y5G+MMT3Ikr8xxvQgr9UdRcQFjgJjqvqUiDwHHGpu7gfmVPXwNcccAp5b9a3bgX+qql+90bmGh4f14MGDrYZmelgUKaEqYQRhGNGIlDBS4i7M0sI7xPs5DjgiiMSfA88h8KxsZDrHsWPHplR1V6v7t5z8gS8CJ4EigKp+dmWDiHwFmL/2AFV9DTjc3McFxoBvr3eigwcPcvTo0Q2EZnpBpR5ydrzE5EKVmVKNxUpjJXejkeK5ztWk7TqtJP5YpPENI4yUMIw/V8OI3f1p3nfnIMOF1M27KGO2iYi8tZH9W0r+IrIf+ATwR8AfXrNNgM8Aj6/zNh8BzqrqhgI0JoqUs+MlXjo7Qz2MSHkugecwkPWJ//y2xhHBcQXfBfz4e6rKXKnGXx69xF278xy+bYB8eiNlJWPaW6t/zV8FvgQU1tj2GHBFVc+s8x5PA99oPTRjYi+8OsG5yTJ9GZ9ixt+Rc4oIhYxPTpU3JpZ4Y2KJB27pY7AQkPZdUp5DyndI+e6GnjKMaRfrJn8ReQqYUNVjIvLhNXb5HOskdREJgE8CX77BPs8AzwAcOHBgvbBMj5hbqvHWZJnhfLAtpfyNckQYyAU0wohfnZ+La5lWxaEKvisM5gMOH+xnd386kTiN2ahWSv6PAp8UkSeBNFAUka+r6udFxAM+Bbx3nff4OPALVb1yvR1U9VngWYAjR47YhEMGgNOXFnEdSTyheq7DYP7ddf+qSqQwt1Tn+RPj7O5P81v3j5AO3ASiNKZ163ZnUNUvq+p+VT1IXHXzI1X9fHPzE8ApVb24ztus+3RgzGq1RsRLr0/zyoWFHavq2QwRwXWEfNpjKB9webbC5blK0mEZs66t9mV7Vz2+iOwVke+uep0Ffhv41hbPZXqAqvLW5BLfevEiJy8uMpgPOqZOXURI+Q7np5aSDsWYdW2o+4KqvgC8sOr1F9bY5xLw5KrXZWBoswGa3rG4XOfF12e4MFUmn/YYzAdJh7Rh2cBlbHqZWiOycQKmrVnfNZO4MFJOji3wizdmEYGhhBp3t4PnOjSiBn/x0hiP3jPM3gFrADbtyYomJlG1RsRf/eIyR1+foZD26M92buJfMZCLn1i+//I4r1xYSDgaY9ZmJX+TqOPnZpkp1RjqslG0mcAlUmV8rsKDB/qSDucdVJWpxRonxxYYn11GJO5N5UjcbiGA4/D29xAcpzn1BeC5wkhfmqFCisF8YNVbHcqSv0nMTKnGyYuLDGTbtzfPVniOMFuqMT5XYaQvhZPwE42qcm5yiV++Nc9MqYbnOuRSbnPbykc8L5KGgMazZ6wcS/wtVOH81DKgKFDM+BwYznJob4FCG/fMMu9kyd8kZrkWxqXJDunNs1GB57C43OB7xy9zaG+BDx4aTjSeqcUaf/3rSXKBu23tKqpKrRHx6sUFXrkwz117CjxwSx99XXpD7yaW/E1iUn53VxeICMWsT6UWMl2qJR0OUwtVHBGyqe37Zx93b3VJ+e7VOZjOXFrkttE8jxzsp2g3gbbV3f/6TFtLec7b9QpdzHGEai0kipK72Eo95PTlRdI38YbrOPFUGAP5gDcnSvz8zPRNO5fZOkv+JjEp3+2F3I/vCguVBj87PZXIDaBcbfD9E+PML9XJ7sC0EyvzIY3PVaiH0U0/n9kcq/YxifHduDdJpJp4Y+jNJCIM5wNeu1SiESqH9hXZVUxt68jlehgxU6oRRfFcQ3HDLYSqHD07S6UWMrCDg+YcEVSVC1Nl9g5mSHlOx3fh7TaW/E1iRIRM4BKGiuN1d2IQEYYKAW9Nl3lzcgnPdTi4K8fBXVlG+9L4m+guqapML9Y4e6XEmcslQlWEVTVpzS98TxJpgA18l5+cmiLSuIqvP+cz3Owemk975DMemcDt6ht/O7PkbxIzMV+h1ohI+70xA6YjwkA2Ln2HkXJuYonXxxdxJJ4YznMdfFfwXQffi78OvPjroPk914knkitVGpwcW2BxuY4jQqF5fDvJpzxIxTepMFIWynWmFqqECo7E9yZHhP6sz6G9Be7eu9ZyIeZmseRvdtxStcHxN2c5c7lE2nd7cpCQ67xdGo8ipREp1XpIpRZXg0VRXG0T6apqHHhHKTmbchnMtf+IaBHBcwXPdcheM5YvjJTlWsjPTk8z2p+2LqI7yJK/2TFhpJy+tMjRN2ZRVQbzgT3yE/eSCbp0rMN6XCeu+quFET8/Pc1vPzxqfxM7xJK/2RGzpRovvDrB3FKdvqyP32ZVFCZZxbTH2Mwy5yaXuH0kn3Q4PcGSf5daroW8ObHE6cuLlJYbGz7ec4VCJp5orT/nU8j4ZAOXbMolvYlGupNjCyyWGwx32Rw+ZnuICMWMx0tnZrh1ONcxazh0Mkv+XerYG7OcGotXwcqnN/5rVlWWKiFzS0vUxyMEYaUriQgUMj6DOZ/7buljV3H9hD65UCWT6o2GXbM5Kd9lqVQjjNSS/w6w5N+l6o2IXCruSrc5gudChncfH6nSCJWL08u8MbHEbSM5Hrlt4LqNdWGkzJXr9FtjnjFtw5J/l9KbOHbWESHw4m6IqsrF6TLnJsvcu6/ASF86nqytOTUwArV6BIo15BnTRiz5d6EwUsrVeMbMm01E6MsGhJHy2qVFTo0tEvfglpX/QOnJ7pzGtDNL/h2sWg+5NFshjJRIlTCMB9NcmC4zuVBlaAeH87vNSb2MMZ3Bkn+Hmlyo8je/nqBUjVfdUI1L2SsDajp5HVzTu3phor92Ycm/DanGIz5r9YhaI6LaiD/XGhHlaoPF5QZnLi+Sbi7KYUw3WKjEnQI86+mzIyz5t4EwUn751hxjM8ss10Iq9Ygw0nfW2evbs196zakB2m0uF2M2q1oP0Qh+64GRrl3Zrd1Y8k/YwnKdH786ydRijWzg4jnxYBfrGWN6RSOMWKg0+OhDoza3zw6y5J+QeDHtMj89NQUCgznf6uhNT4lUmS/XAfjAXUPsH8omHFFvseSfgHoj4ugbs5wcW6CY9kj1yJTGxkBc8FmsNKg1Iu7anefwbQObGoVutsZ+4jtouRYytVjl2NlZ5pZqDNmslqbHhJEys1RjT3+a37hz0OZ6SpAl/5tEVSnXQqYXa4zPLnNxZpmF5TpCPDJ2yP7oTQ+q1EP2DqT5nYd3WzVnwiz53wSRKs+fGGdivooQr1qUCTpj4Q1jbqYwUgoZa99qB5b8bwJVmJivWiOuMU2q2hynEsbLO5rEtdxRXERcETkuIt9pvn5ORE40P86JyInrHNcvIt8UkVMiclJEPrBNsbc9S/zGNNu6SjUcR7hnX4GDI7mkQzJsrOT/ReAkUARQ1c+ubBCRrwDz1znunwPPq+qnRSQArD+XMT2iVGkQqfK7793b0roPZue0VPIXkf3AJ4CvrbFNgM8A31hjWxH4EPAvAVS1pqpzW4jXGNMhFit1RODjj+yxxN+GWq32+SrwJSBaY9tjwBVVPbPGttuBSeBfNauMviYiaz7zicgzInJURI5OTk62GJYxpl3VG8qj9wwzaPNPtaV1k7+IPAVMqOqx6+zyOdYo9Td5wHuA/1NVHwGWgP92rR1V9VlVPaKqR3bt2rV+5MaYtrfQHMFr2k8rJf9HgU+KyDngT4HHReTrACLiAZ8CnrvOsReBi6r6YvP1N4lvBsaYLpf2Hc5PlZMOw1zHuslfVb+sqvtV9SDwNPAjVf18c/MTwClVvXidY8eBCyJyqPmtjwCvbj1sY0y7SwcukwtV6o21aotN0rba4fZprqnyEZG9wNdU9cnmt/4h8H83e/q8Afz+Fs9pjOkAjgiq8P/85Dy5tEdfxqM/F9CX9cmlPbIpl1zKw3fFukUnYEPJX1VfAF5Y9foLa+xzCXhy1esTwJFNxmeM6WCD+QBVpR4q06Ua43PxsqPiCAJECoErFLI+/VmfgVxAPuORS3kU0h7pwCY9vFlsqF0LVJXZpbjhypF4vVrHEVxpfnYER2xQlzFrERECL57Timt6fKoqkUKlGnKhXOfslaV4OVLiJR3v2pPnvv1F+m196G1nyX8d9TDipddnOH1pEadZWoFVa43q21+v3BhcR2iEVs9pzHpEBFfADVzSvLOUH0bK6+MlTl9a5JbhLA/d2s9wwebH2i6W/G9gcbnOC7+eZHqxyuA60y+rKko8r4+qkrZZO43ZEtcRBnJxtdH4bIXz05cYKaY5fLCfPQNpmw59i2wR2OsYmynz745eYmG5zlAhte4fmojgSFzq91wH19YhNWZbiAjFrM9QLmChXOf7J8a5OL2cdFgdz0r+QBQpf3d2hoVynXoYUQuV2VKNfMoanIxpFyJCPu1RDyMm5iscGLZpwrbCkj8wNrPMry8skE25OBI33g7mAhwrvRvTdrKByy/fmufSbIX79xe5Y3c+6ZA6Us9X+6gqx8/NXu1znAlcUr5rid+YNpXyXYYLAYvlOq9dWkw6nI7V88n/ynyV6cUaWaveMaZjiAi+1/Ppa0t6+qenqpw4N0fKc6z7mDGmp/R08p9ejEcc5tPW9GGM6S09m/zL1Qa/eHMWz+YVMcb0oJ4r8i5VGrx6cYGTYwsoMJD1kw7JGLMJy7WQ3QPppMPoWD2T/EuVBr++MM+pS4ug0Jf1bSCWMR1qvlyjP+fzvjsHkw6lY/VE8p9arPLd45dRhb6MJX1jOlmp2iDluzzx0Chp33rpbVZP1Pm/enEBR4TBXGCJ35gOF4banCoaao0IVV3/IPMuXV3yDyNlplTjzStL9Oesbt+YblDMeFyYLjM2vXx1Rt2U75D2XTJB/JFNuWQDl1TgEnjO1Y9U87PrWEePrkr+s6Uas0s1JuarTCxUmV2qAeB7js0AaEyXEBGG8m/PmquqhJHSCCPmlkKmF+PXYaQgIPH/rq4RoBrPGLqrGHD//j72DmZ6skaga5L/ci3kO8cuESp4jpDyHfqzviV9Y7qciOC5greB6v9IldlSnR/+6grpwOW+/UXuGM2T66ExP11zpfVGhAJDeVvxxxhzY05zhtB82qPWiDj+5izH35zjwK4s9+4rMtK3/jTuna5rkn9ojT7GmE0IPIfBfIpIlbHpZc5NLlHM+Dx0Sx937sl3bdtA1/T2iSJL/saYzXNE6Mv6DOdTRKHy/702xVI1TDqsm6Zrkn9oyd8Ys03SgYvrxGt7dKuuSf7RqoXUjTFmO3RzL6CuSf6hZX9jzDZSpasXdequ5N+9vydjzA5SVVS1q0v+bdnbJ1K4Ml/Z0DHz5boN8zbGbAuFZp2/Jf8dNV+u8fzxcTb6c7dJnowx20Gbg0W7WVsmf1UYtMFaxpiEqCqu293Jvy3r/Lv7R26MaXeRgue0ZXrcNi2X/EXEBY4CY6r6lIg8Bxxqbu4H5lT18BrHnQMWgRBoqOqR9c5lNffGmCSpKr5nyX/FF4GTQBFAVT+7skFEvgLM3+DY31LVqVZPZO22xpgkRQqeVfuAiOwHPgF8bY1tAnwG+Mb2hmaMMclQVXy3u0v+rV7dV4EvAdEa2x4Drqjqmescq8APROSYiDyz8RCNMWZnWckfEJGngAlVPXadXT7HjUv9j6rqe4CPA/9ARD50nfM8IyJHReRopVKxPvvGmMSoatc3+LZydY8Cn2w23P4p8LiIfB1ARDzgU8Bz1ztYVS81P08A3wbed539nlXVI6p6JJtN07CJ2owxCVEFv9dL/qr6ZVXdr6oHgaeBH6nq55ubnwBOqerFtY4VkZyIFFa+Bj4KvLLeOV1HqDXWqmEyxpibL+qB3j5bvbqnuabKR0T2ish3my9HgZ+IyMvAS8Bfqerz672p5wh1S/7GmIQoNsL3HVT1BeCFVa+/sMY+l4Anm1+/ATy80aBcx7GVuYwxiYkUK/knwXGka5dOM8Z0AOvqmQzXwYb5GmOSI929kAu0afJ3RECw7p7GmMRY8k9IIe1RDy35G2N2noh09Spe0MbJvz/nW3dPY0xirOSfkMFcYMnfGJMMteSfmL5cgFqrrzEmAQrMlGpd3e7Ytsk/n/ZwbFkXY0wCCmmPn5+e5oe/mmBxuZ50ODdF2yb/XMq1cr8xJhGB5zBcCBifW+bbL41x8uICYZfNN9a2yT8TuDiOEHXZD9wY0xlEhL5sQD7t8fMz0/zg5fGuugG0bfIXEUb7UsyV611d72aMaW++6zCUDxifq/DmxFLS4Wybtk3+AB+6bxcjfSlmluwGYIxJjohQSHsce2OGRtgdvRDbOvmnfZcnHhrlwHCGqVKNyG4AxpiEpHyXcjXk9fFS0qFsiw3N6pkE33X4zftGSJ+Z5tSlRQZzQdf3vzXGtKdixufF12eYmK/Sl/MpZnyygUs25ZJNeR2Vm9o++UM82OL9dw+R9l1OnJtjIOfjdfmMe8aY9hN4Dnk8Lk6XeXMiHokkAmg8NiCbcilmfAZyPvfuK1LM+glHfH0dkfwhnuztkdv6yQQuPzszbU8AxphEBJ5DsMZc/6pKI1Lml2pcma/w1mSZjz+ym0KmPW8AHVV8FhHu3V9kKB9Q75JGF2NMdxARfNchm/IYzMU56vsvX2G5FiYd2po6KvkDRJEyX66veec1xph2Ucz4lCoNrsxXkg5lTR2XQUuVBkpzzn9jjGlrStimU9N3XPJfrDRslS9jTEdQaNsq6o5L/nNLNZvt0xjTEVyEaptOTd9xyX9yoUpg3TyNMR1AHKFWb8/k3zFdPVdMLdYQEWqNCAFcV6z+3xjTlhyhbRel6qgitKoyUgxIBy4AoSpTi7WummnPGNMdGmFEpR627bQ0HVXyFxE+dN/I1deqyvFzc/zy3BxDhVSCkRljTKzWiFhcruO5Du+5bYB79hWTDmlNHZX8V4sipVwLmV6s4VobgDEmYZVaSKnaIB24vO+uIe4YzZHy3aTDuq62T/6qyuxSncXlOvPlOjOl2tXXK09TA7n2HD5tjOl+1XrIYqVBIePz2L3DHNyV64i5x9o6+S9VGrz4+jTnp5avruYbeA6+K/RlfWvoNcYkrlQN+eChYe7ancfpoPnG2jb5n5tY4ienplBgMOcjluiNMW0oUkW180Yfte2zyd+dnSHwHPqzlviNMe2rmPb42elpvv3SGOcmljpm3fGWS/4i4gJHgTFVfUpEngMONTf3A3OqeriVY1s5X60RkUu17YOJMcYA8QpfKd9luRby17+eYCAX8Bt3DrJ3IN3WBdeNZNcvAieBIoCqfnZlg4h8BZhv9dhW1EOlg6rPjDE9LhO4ZAKXpWqDH7w8zkgxxXvvGGS0L9WWN4GWqn1EZD/wCeBra2wT4DPANzZ67PXEvXi0LX9gxhhzI7mUx1A+YGG5wfeOX+YHL19haqGadFjv0mqd/1eBLwFrjVN+DLiiqmc2cexVIvKMiBwVkaNTU1MIlviNMZ1JRMin45vA1GKVvzx2iR+9MsFsqZZ0aFetm/xF5ClgQlWPXWeXz3H9Uv96x16lqs+q6hFVPTI4NLTe7sYY0/ZEhGLGZygfMDazzL/9uzF+cnKSpWoj6dBaqvN/FPikiDwJpIGiiHxdVT8vIh7wKeC9Gz32Rid0RMilPebLNfqyQetXY4wxbUhE6M/6RKqcvryIiPDoPcOJxrRuyV9Vv6yq+1X1IPA08KNVyfsJ4JSqXtzEsdclAh8/vJtc2md2qX0ek4wxZiscEQZyAW9NJt8ldKv9/J/mmiofEdkrIt/d4vuSS3t87PBuBvMBM0s1tE1nxjPGmI3wXId6qEwnXP8v7ZhUjxw5okePHgXi/v4/PjnJxenmFA8CGd8l7TvWG8gY05EWK3VcR3j8gVF2FbdnRmIROaaqR1rdv+1HUQWewxMPjlKph8yUakzOV7kwXWZ68e27Zsp3yPhuR82rYYzpXYW0T7na4K9+cZkP3j3EXXvyO16YbfvkvyLtu+wdyLB3IMPDB/upNyJmlmpMLVS5ML3MxHwFVVCUlOeSTbk28Zsxpm1lUx6+F/GTU1MsVRs8ctvAjp6/Y5L/tXzPYbQvzWhfmvtv6SOMlNmlGlOLVc5dWeLKfIXBvC3wYoxpX77rMJgPePmtOe7YnaeY2bnp6dt2YreNch1huJDinr1FPnjPMArWSGyMaXuuE69DfuLc3I6et2uS/2rFjM9ALqBSb8+Fk40xZrW+rM/Z8RIzO9gDqCuTP8ChvQVKlYaV/o0xbc8RIfAc/vqVCd6cKNEIb37BtWPr/Ndz+2ie81NlxmaWGcoH1i3UGNPWipm4B9DfvDpF2nd48EAfd4zmSQc3Zx3gri35B57DRx4c5c7deaYWa4QdssCCMaZ3ZZszgvquw9+dneXPfnaBn5+ZZr5c3/ZzdW3JH+KGlEfvGSab8nj53BwDOb8jFlY2xvS2wHMYygeEkXL60iKnxhbYP5jlwVv7GO1Lb8s5uj4TOiK857Z+PnBoiNly3Z4AjDEdw3XiuYAGcwFX5is8f3x829YG6PrkD/GMevfuK7KrkKLesB5AxpjOsjI1dOA5/PTU1LYUYnsi+a+oh5FNAWGM6Vj5tMfMUo2TF2+0am5reir51xqRrQtsjOlo/VmfY2/ObbkRuKeSfz1UK/kbYzqa5zq4jvCz09NEWxjH1DPJP4qUMIpsZWBjTMcrpj0uzS5zdry06ffomeTfiBRBbLCXMabjiQh9GZ+XXp/ZdCeW3kn+YYQV+40x3SLwHBqRUtvkVBA9lPytf78xpruIwmar/Xsm+Qe+QzZwmS/bgvDGmC4hm5+6vmeSf9p3efKRPRSztiC8MaY7qMa9GDejZ5I/QC7t8bHDu9k/mGG6VNtSNyljjEma5wo//NWVTfX576nkD3EjyYfvH+HAcJbSciPpcIwxZtOKGZ96GPFXv7i84WN7LvlDPFnSvsEM9cjm+THGdLZC2t9UR8aeTP4QVwGJ9f00xnQBz914Luvd5J/yLPUbY7pCtIlZPns4+buII1TrYdKhGGPMlmxmhueeTf4p3+W37t9FqRraDcAY09E203OxZ5M/wP6hLE88OEKpGlKxG4AxpkNtptd6Tyd/iG8Av/3QCOWa3QCMMZ1pM3X+LS/gLiIucBQYU9WnROQ54FBzcz8wp6qHrzkmDfwYSDXP9U1V/WcbjvImmi/XmZiv4rsOtUZE2neTDskYY1pWrjbQTfReaTn5A18ETgJFAFX97MoGEfkKsNa6YlXgcVUtiYgP/EREvqeqP994qBtXa0RMzFdoRArxf0A8F0YjVM6Ml5hcqCJAIePhuz3/IGSM6SALy3VEhCcf2bPhY1tK/iKyH/gE8EfAH16zTYDPAI9fe5zGE+isrDbgNz9u6pwKkSqTC1VeHy/xxpUlokjR5ikFQYlndlZV0oHLYM63Of6NMR1FVZlbqlPI+nzkwRGKGX/D79Fqyf+rwJeAwhrbHgOuqOqZtQ5sVhcdA+4E/ndVfXHDUW7AT09N8fp4Cc8VCmkf15ZtNMZ0kUiVmVKNvYMZfvO+XZuuql63nkNEngImVPXYdXb5HPCN6x2vqmGzLWA/8D4ReeA653lGRI6KyNHJycn1I1/7XFyYXqY/69OfDSzxG2O6Shgp06Uad+8t8MSDo1tqo2ylkvtR4JMicg74U+BxEfk6gIh4wKeA59Z7E1WdA14APnad7c+q6hFVPbJr166Wgr9WpR5Rb0R4VndvjOkSkSrVeshCuc7MUo0jtw/wwbuHtly4XbfaR1W/DHwZQEQ+DPwjVf18c/MTwClVvbjWsSKyC6ir6pyIZJr7//GWIr6BhXLdlmo0xnQkVaUexom+FkZvzz0m0Jf12TOQ4cBwhluGc9tyvo309lnL01xT5SMie4GvqeqTwB7g3zTr/R3gz1T1O1s853XNl+voZsY5G2NMghaX61QbEYWMz97BDLuKKfqyPoWMTz7t3ZQq7A0lf1V9gbjqZuX1F9bY5xLwZPPrXwKPbCXAjZhcrFp3TWNMR1FVqmHEk+/Zw2hfesfO21WZstaIsF6bxphOslQNGe1LM1JM7eh5uyr53zqcpbbJ9SyNMWanqSqVesjhW/t3fLzRVuv828regQy5tMtMqYYI9GdtAJcxpn1VGxH9uYDdAztX3bOiq0r+6cDlP3rffj55ZC9p3930qvbGGLMTlqshtw5ncRIopHZV8gdwHGEgH7B3MEOlZrN0GmPaV6RKNpXMZJJdl/xX7OlPUbdun8aYNiYipBKaSbhrk/9APmXjvYwx7U3A38Ti69uha5N/X9Ynm3KZL9fQzSxzY4wxN1EYKVGkZAIr+W8r1xE+/sgehgopZpbqSYdjjDFXqSozSzUOH+xnqLCz/ftXdG3yB8inPT768G4CV2iEUdLhGGMMALNLdW4fyfHwwf7EYujq5A/xE8BIX5pK3ZK/MSZ588t1BnI+Hzw0nEgXzxVdn/wB9gykqTUs+RtjkrVQqeM5wuMPjhJ4yabfrhrhez1DhZRN9WyMSYyqMluuU0h7PPHgKPl08qk3+Qh2QH/ORzX+Bdh0D8aYnRRpvPrWLUNZHrt3eEurb22nnkj+ad8ln/aoh0rgWfI3xtx81XpIqRqCKvfvL3LkjsG2Wlq2J5I/wJ7+NGfHSwT5IOlQjDFdRlWJFBphxHI9JIyUbMrjvbf3c+uuHMWMn3SI79Izyf+BA32Mz1eYLtUYyPmJtrIbYzrD6qQeRkojVBqrpo0RAQVU45G6uZTHHYMZ7tidZ1cx1dZ5pmeSf1/W53ffu5ejZ2c4NbZIX9ZPvLXdGJOcMNJmQo9oNBO7ojgIKoDGid1rJvViNq4+LmTi2QPSvkvad0j7LinfweuwVQR7JvkDBJ7DB+4eYs9Ahp+emmK5HlJMe9YIbEyXW6l/dyQupSuK5zhkUy79+SBO6mmPXMojHcTJfCWpd+vSsD2V/CGeRe+2kRxDhYAfn5xktlSnP9t+9XHGdLpINdFqD1VlsdKg3ojIpT3ed+cgI8XU1eTuOdLTBb+eS/4rihmfjzwwyp/97YXE/0iN6RaqSrkWslwLryZWaVahQPz07buC7zk37d9cPYxYXG6gwC1DGe7dV2R3fxqnjXratIOeTf4AmcDlluEsl2eX27I13phOoaqUKg2qjYjhQsCjh4bZO5ih1ohYqjQoVRsslOvMlGrMLtWZL9dX7geogiMrNwaHlXuFql5tTF399cr5VNfYDwg84aFb+7hzd56C/bu+rp5O/gCH9hY4P1VOOgxjOlKkysJynUao7BvM8NCt/Yz2pa6W+jOBSyZwGeadM1dGUfyEsHJjmFuqMVN6+6bgOoIjguMIrgiuE384DjirXq/+cByhP+uzbzDTcY2vSej55L+7P03Kc6iHUdc27Biz3cIoTvqRwm0jOR64pchgPmi5Dt1xhHzaI5/2GL3JsZq19Xzydx3h0N4Cr1yYZyBnA8CMuZFGGDG/3ECAQ/sK3LevSNE6THSknk/+ALeN5vjl+Tmb+8d0hFojImrWecPbdd4r/dLX+l789TXfW3lN3AsOVRCQq7XuzX8LqxpsXUd4+NY+7t5bIJey9NHJ7LcH9Gd9BvMBpeWGNRCZtlapxz1pBvIBjsT13/EHiCNXv7dSB+4QV7E4IogD7qp6dEdWbZP4BiDN469+hndsG8oHiS04braXJX/iP+xHDw3z169MMFOq0p8LrOunaUu1RsRtIzkeu3dX0qGYDmctnE1DhRSf/I193LE7z/RizRZ/MW2pEarVsZttYSX/VQLP4dFDw+wbzPC3r01TqjTesQjMWs8CK/Wuq6pFgXjQWL81IJttpkChDRYCMZ2v5b8iEXGBo8CYqj4lIs8Bh5qb+4E5VT18zTG3AH8C7AYi4FlV/efbEPdNE0//kGf/UJZGqERR3CgWRXq1kS1qzvS3MuPf6u+vfP7b16ZphFFX9zdeGWizUidsbj6RuO+8MVu1kSLEF4GTQBFAVT+7skFEvgLMr3FMA/ivVfUXIlIAjonIv1fVV7cQ847wXYettGudHS/x1mSZ/pzfEeMHIo1nOIyaMx2GkRLq2ze9lUbAq083zR4jriOEzSlus0E8y6HdCG4ewZK/2R4tJX8R2Q98Avgj4A+v2SbAZ4DHrz1OVS8Dl5tfL4rISWAf0PbJf6vef/cQg/mAk2MLzC/XSXsuuZTbVolx9ZB81xHSvkPKdykELhnfvTo6Mx24+K7Ew+89h8BtfvYcXEcoVRqMTZc5dWmRmVINBHIpj5RnN4LtpmrJ32yPVkv+XwW+BBTW2PYYcEVVz9zoDUTkIPAI8OIG4utYuZTHI7cN8OCBPi7NVnj14jzjc9V3THIVqVJIezvWdS7SeM7yehhRb8RzmO8ZyPDIwX5GVg3J36h82uPQviJ37y2wuNzgwnSZ1y4tMrNUR5rbbe2ErYsixXHEfpZmW6yb/EXkKWBCVY+JyIfX2OVzwDfWeY888OfAH6jqwnX2eQZ4BuDAgQPrhdUxPNfhwHCWA8NZ6s3BOSvtBONzFV48Pc1SLaQ/u/XVxVYWpqiHby9QsdJXW+FqiXwwH9CX9bl9NM9wofUh+esREYpZn/uzfdy3v8jcUp3zU2VOX15kulTFdYS+jG9PA5vUiJS8rT9htomo6o13EPkfgP+EuP4+TVzn/y1V/byIeMAY8F5VvXid433gO8D3VfV/aSWoI0eO6NGjR1u/ig5WqYUcPzfHa5cW2eyMsyuzGXqOkM94FDMefVmfYsYnm/KuVt+k/Js3je6NRKrMlGq8cn6eN64sUch4pG2g0IaVKg12FVM88ZDNhmPeTUSOqeqRlvdfL/lf8+YfBv6Rqj7VfP0x4Muq+pvX2V+AfwPMqOoftHqeXkr+KxaW64Rh67+L1RxHyDTr5du5VKiqXJgu87evTVOtR/TbWso3pKrxlMjVEEXxXYe/d9cgd+5eq/bV9LqNJv+tdhh+mmuqfERkL/A1VX0SeJT4qeFXInKiucs/UdXvbvG8XacX1hMQEQ4M5xjpS3P8zTlOjS0gxDcv33Xw3Piz08NdR1cWQ6nU40GGhYzHQ7f2sW8ww1AhhWsLkphtsqGS/07pxZJ/L5pcqDJfrrNQrrGw3GBhuRH3PqqHBL5DMd39N0SI22qWKg3qkSLASF+K20fz7OlPU8hYHb9pzU6X/I3ZtF3FFLuKqXd9f75c5/njlylVGuS7dDSrqjJXjufD91zh1l1Zbt2VY6QvZe0hZkd0578s09H6sj4fPbyb7x0fp1xtkO3CqYOrjXjxoMcfHGHYqnNMAqzDsGlLA7mAjz48Sj1Ulmth0uFsu5TnUKmHhJFa4jeJsORv2tZwIcVvPzxKpR5SqXfXDUBESPsuJ87FiwgZs9Ms+Zu2NtqX5iMPjlCqNK7OIdQtcimXifkqc0v1pEMxPciSv2l7+wbjEdLdVv0jzdW0Qiv5mwRY8jcd4cBwlmqXLrDjWZ2/SYAlf9MR1uoS2g1WpsU2ZqdZ8jcdoZj1SXkO9bDLSv9KVy/4Y9qX/dWZjuCIcM++AnPlOsu1sGt6yFjJ3ySl+0bPmK51+GA/u4opXjo7c3WtgJX5qhXFd+JFZjxX8Jz2nuTubWp1/iYRlvxNxxAR9g9l2T+UJYziwV/LtZByNZ4TaK5cZ6FcZ2G5zsJy9K6bQzbw2moVrJXlMR1L/iYBlvxNR3IdIZ/2mnP/vLsxeOXmUK42KNdCSssNTl9eZKpUJeO5ZBNeUnMlrjt25xOLwfQ2S/6mK73z5hC7b3+RS7PLnDg3x+RClcBzKOzQyliNMKLaiD8iVYppn995eBd7BtI3/dzGrMWSv+kZjhNXG+0bzDAxX+Xl83OMTS/juUIxsz0Ly2hzneRqI6LWiK5WO/mew3AhYLQ/zWA+YE9/2nr5mERZ8jc9R0QY7U/z0f7d8fKSF+LlJVfaCERBhbfXPubtr0XjF81P76DES2pmA5fdA2lG+9L053z6sj7ZINlqJmOu1ZaLuYjIJPBWQqcfBqYSOvd2sutoL3Yd7aUbr+NWVd3V6oFtmfyTJCJHN7IaTruy62gvdh3txa7DBnkZY0xPsuRvjDE9yJL/uz2bdADbxK6jvdh1tJeevw6r8zfGmB5kJX9jjOlBPZn8ReS/EpFfi8grIvINEVlzmKWI/IaIhCLy6Z2OsRXrXYeIfFhE5kXkRPPjnyYV63pa+Z00r+dEc7+/SSLO9bTwO/lvVv0+Xmn+fQ0mFe/1tHAdfSLylyLycnO/308q1htp4ToGROTbIvJLEXlJRB5IKtYbEZEvNq/h1yLyB2tsFxH5X0Xk9ea1vGfdN1XVnvoA9gFvApnm6z8DvrDGfi7wI+C7wKeTjnsz1wF8GPhO0rFu07X0A68CB5qvR5KOe7N/W6v2/13gR0nHvcnfxz8B/rj59S5gBgiSjn0T1/E/Af+s+fU9wA+TjnuN63gAeAXIEg/M/X+Bu67Z50nge8RjD98PvLje+/ZkyZ/4B5gREY/4B3ppjX3+IfDnwMROBrZBrVxHp1jvWv5j4Fuqeh5AVdv197KR38nngG/sSFQbt951KFCQeNhynjj5N3Y2xJasdx33AT8EUNVTwEERGd3ZENd1L/BzVS2ragP4G+DvX7PPfwD8icZ+DvSLyJ4bvWnPJX9VHQP+Z+A8cBmYV9UfrN5HRPYR/3D/r52PsDWtXEfTB5qP5t8Tkft3NMgWtXgtdwMDIvKCiBwTkd/b6TjXs4HfCSKSBT5GXMBoKy1ex78gTkqXgF8BX1TVtlpmrcXreBn4FICIvA+4Fdi/k3G24BXgQyIy1Py7eRK45Zp99gEXVr2+2PzedfVc8heRAeK75G3AXiAnIp+/ZrevAv9YVcMdDq9lLV7HL4iHfD8M/G/AX+xokC1q8Vo84L3AJ4DfAf47Ebl7RwNdR4vXseJ3gZ+q6sxOxdeqFq/jd4ATze2HgX8hIsUdDHNdLV7H/0hcqDhB/LR/nDZ7glHVk8AfA/8eeJ74hnVtjGtNHHXDrpw9l/yBJ4A3VXVSVevAt4APXrPPEeBPReQc8Gng/xCR/3BHo1zfutehqguqWmp+/V3AF5HhnQ91Xa38Ti4Cz6vqkqpOAT8GHt7hONfTynWseJr2rfJp5Tp+n7gaTlX1deK69Xt2OM71tPpv5PdV9TDwe8TtF2/ueKTrUNV/qarvUdUPEVexnblml4u882lgP+tUA/di8j8PvF9Ess36yo8AJ1fvoKq3qepBVT0IfBP4L1X1L3Y80htb9zpEZHdz28ojrQNM73ik61v3WoB/CzwmIl7z0ffvrbFP0lq5DkSkD/hN4mtqR61cx/nm92nWkR8C3tjRKNfXyr+RfhEJmi//M+DHqrqww3GuS0RGmp8PEFdTXVtw+HfA7zV7/byfuIrr8o3es+emdFbVF0Xkm8RVIg3ix7xnReQ/b25v23r+1Vq8jk8D/4WINIBl4Gltdg1oJ61ci6qeFJHngV8CEfA1VX0lsaDXsIG/rb8P/EBVl5KJ9MZavI7/HvjXIvIr4iqHf9x8ImsbLV7HvcCfiEhI3JvsP00q3nX8uYgMAXXgH6jq7DXX8V3itoDXgTLxk9kN2QhfY4zpQb1Y7WOMMT3Pkr8xxvQgS/7GGNODLPkbY0wPsuRvjDE9yJK/Mcb0IEv+xhjTgyz5G2NMD/r/AZl9OmJKUo/yAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "GeoJSON output:\n",
      "{\n",
      "    \"type\": \"Polygon\",\n",
      "    \"coordinates\": [\n",
      "        [\n",
      "            [\n",
      "                8.969565987549165,\n",
      "                47.27922900257082\n",
      "            ],\n",
      "            [\n",
      "                8.973652061631208,\n",
      "                47.290534195958585\n",
      "            ],\n",
      "            [\n",
      "                8.991620623576104,\n",
      "                47.29457086136035\n",
      "            ],\n",
      "            [\n",
      "                9.005352031243234,\n",
      "                47.30990758297694\n",
      "            ],\n",
      "            [\n",
      "                9.006038422582417,\n",
      "                47.32756532234069\n",
      "            ],\n",
      "            [\n",
      "                9.000676912362332,\n",
      "                47.333805471963935\n",
      "            ],\n",
      "            [\n",
      "                9.000545281887767,\n",
      "                47.34476201996809\n",
      "            ],\n",
      "            [\n",
      "                8.992305974992021,\n",
      "                47.361493894815304\n",
      "            ],\n",
      "            [\n",
      "                8.975849481486328,\n",
      "                47.37098348005585\n",
      "            ],\n",
      "            [\n",
      "                8.971706445238945,\n",
      "                47.381015667238294\n",
      "            ],\n",
      "            [\n",
      "                8.965216675256876,\n",
      "                47.38771797063729\n",
      "            ],\n",
      "            [\n",
      "                8.963468048272603,\n",
      "                47.39542209737156\n",
      "            ],\n",
      "            [\n",
      "                8.938528199508504,\n",
      "                47.41237663879035\n",
      "            ],\n",
      "            [\n",
      "                8.939326185746767,\n",
      "                47.416388386520254\n",
      "            ],\n",
      "            [\n",
      "                8.953853363988665,\n",
      "                47.425174262141304\n",
      "            ],\n",
      "            [\n",
      "                8.95316675433757,\n",
      "                47.44236734388886\n",
      "            ],\n",
      "            [\n",
      "                8.949046956742931,\n",
      "                47.44794326963111\n",
      "            ],\n",
      "            [\n",
      "                8.930548868582774,\n",
      "                47.4511054833359\n",
      "            ],\n",
      "            [\n",
      "                8.926388586998849,\n",
      "                47.46142306719717\n",
      "            ],\n",
      "            [\n",
      "                8.917004302802026,\n",
      "                47.46922662315857\n",
      "            ],\n",
      "            [\n",
      "                8.918834554399174,\n",
      "                47.53292821138423\n",
      "            ],\n",
      "            [\n",
      "                8.901670907392289,\n",
      "                47.54218188909647\n",
      "            ],\n",
      "            [\n",
      "                8.868865916405799,\n",
      "                47.5455407012233\n",
      "            ],\n",
      "            [\n",
      "                8.867357739355919,\n",
      "                47.54881691500381\n",
      "            ],\n",
      "            [\n",
      "                8.872142537508253,\n",
      "                47.55560542402674\n",
      "            ],\n",
      "            [\n",
      "                8.872142659867928,\n",
      "                47.56994341167184\n",
      "            ],\n",
      "            [\n",
      "                8.868119394653224,\n",
      "                47.57120256929449\n",
      "            ],\n",
      "            [\n",
      "                8.867336123929583,\n",
      "                47.57642142557468\n",
      "            ],\n",
      "            [\n",
      "                8.854291412530415,\n",
      "                47.58428672633104\n",
      "            ],\n",
      "            [\n",
      "                8.825173457573266,\n",
      "                47.589344393666195\n",
      "            ],\n",
      "            [\n",
      "                8.823081055139692,\n",
      "                47.600285347914735\n",
      "            ],\n",
      "            [\n",
      "                8.828766694319269,\n",
      "                47.60715747816466\n",
      "            ],\n",
      "            [\n",
      "                8.833690511430428,\n",
      "                47.60835457512499\n",
      "            ],\n",
      "            [\n",
      "                8.83427306117602,\n",
      "                47.616541290107925\n",
      "            ],\n",
      "            [\n",
      "                8.843991177800724,\n",
      "                47.623623393572316\n",
      "            ],\n",
      "            [\n",
      "                8.848796590242632,\n",
      "                47.63148658461631\n",
      "            ],\n",
      "            [\n",
      "                8.848796612739507,\n",
      "                47.64814005503311\n",
      "            ],\n",
      "            [\n",
      "                8.843990547900123,\n",
      "                47.65831718497481\n",
      "            ],\n",
      "            [\n",
      "                8.823641734826705,\n",
      "                47.66254039218784\n",
      "            ],\n",
      "            [\n",
      "                8.822017537797622,\n",
      "                47.670347142451945\n",
      "            ],\n",
      "            [\n",
      "                8.812406815795528,\n",
      "                47.67636212035885\n",
      "            ],\n",
      "            [\n",
      "                8.802793954267301,\n",
      "                47.67913813016058\n",
      "            ],\n",
      "            [\n",
      "                8.776017259875442,\n",
      "                47.679600798972565\n",
      "            ],\n",
      "            [\n",
      "                8.76846614932927,\n",
      "                47.67543766450524\n",
      "            ],\n",
      "            [\n",
      "                8.767340005051972,\n",
      "                47.669262700555535\n",
      "            ],\n",
      "            [\n",
      "                8.761873545080382,\n",
      "                47.66272562576313\n",
      "            ],\n",
      "            [\n",
      "                8.75061203239404,\n",
      "                47.66248130108188\n",
      "            ],\n",
      "            [\n",
      "                8.740278501788687,\n",
      "                47.6574096621237\n",
      "            ],\n",
      "            [\n",
      "                8.71573568502211,\n",
      "                47.65880177600464\n",
      "            ],\n",
      "            [\n",
      "                8.704599964023423,\n",
      "                47.66341072524623\n",
      "            ],\n",
      "            [\n",
      "                8.700725496287088,\n",
      "                47.668531573251855\n",
      "            ],\n",
      "            [\n",
      "                8.700481282914796,\n",
      "                47.6745113524843\n",
      "            ],\n",
      "            [\n",
      "                8.696626695330584,\n",
      "                47.67663468730201\n",
      "            ],\n",
      "            [\n",
      "                8.696361298558974,\n",
      "                47.68191665828765\n",
      "            ],\n",
      "            [\n",
      "                8.691151000636543,\n",
      "                47.68630977119764\n",
      "            ],\n",
      "            [\n",
      "                8.689494843480851,\n",
      "                47.694409288384705\n",
      "            ],\n",
      "            [\n",
      "                8.684689419988228,\n",
      "                47.697645804335366\n",
      "            ],\n",
      "            [\n",
      "                8.671018924373122,\n",
      "                47.69929920813748\n",
      "            ],\n",
      "            [\n",
      "                8.659970733819469,\n",
      "                47.70736226728731\n",
      "            ],\n",
      "            [\n",
      "                8.646924966059146,\n",
      "                47.708750312887545\n",
      "            ],\n",
      "            [\n",
      "                8.629762198076458,\n",
      "                47.708750290390654\n",
      "            ],\n",
      "            [\n",
      "                8.60916305778468,\n",
      "                47.70458609140903\n",
      "            ],\n",
      "            [\n",
      "                8.603671808211102,\n",
      "                47.70042327279162\n",
      "            ],\n",
      "            [\n",
      "                8.601367483770083,\n",
      "                47.69335705137093\n",
      "            ],\n",
      "            [\n",
      "                8.587192778589452,\n",
      "                47.68191665828765\n",
      "            ],\n",
      "            [\n",
      "                8.584446196558213,\n",
      "                47.67636437824451\n",
      "            ],\n",
      "            [\n",
      "                8.585667596062308,\n",
      "                47.65890700311242\n",
      "            ],\n",
      "            [\n",
      "                8.574146513941008,\n",
      "                47.649991048039595\n",
      "            ],\n",
      "            [\n",
      "                8.572682827085538,\n",
      "                47.63922267561858\n",
      "            ],\n",
      "            [\n",
      "                8.55975282171292,\n",
      "                47.639480904626474\n",
      "            ],\n",
      "            [\n",
      "                8.551481473465307,\n",
      "                47.64536177094886\n",
      "            ],\n",
      "            [\n",
      "                8.532942044754378,\n",
      "                47.64860060097403\n",
      "            ],\n",
      "            [\n",
      "                8.506165575305115,\n",
      "                47.6486005687317\n",
      "            ],\n",
      "            [\n",
      "                8.50204868660155,\n",
      "                47.646751220206575\n",
      "            ],\n",
      "            [\n",
      "                8.501985122495256,\n",
      "                47.63775413387961\n",
      "            ],\n",
      "            [\n",
      "                8.499170330105198,\n",
      "                47.63461844028525\n",
      "            ],\n",
      "            [\n",
      "                8.489000569580991,\n",
      "                47.634257234230425\n",
      "            ],\n",
      "            [\n",
      "                8.484163914905439,\n",
      "                47.62969629085424\n",
      "            ],\n",
      "            [\n",
      "                8.462220869751345,\n",
      "                47.627316979977785\n",
      "            ],\n",
      "            [\n",
      "                8.450694762418347,\n",
      "                47.6162640000497\n",
      "            ],\n",
      "            [\n",
      "                8.437503946577381,\n",
      "                47.61204820199898\n",
      "            ],\n",
      "            [\n",
      "                8.437499999999986,\n",
      "                47.58096763733005\n",
      "            ],\n",
      "            [\n",
      "                8.413468431750433,\n",
      "                47.58197320821402\n",
      "            ],\n",
      "            [\n",
      "                8.404786136868076,\n",
      "                47.573128889437804\n",
      "            ],\n",
      "            [\n",
      "                8.397677826846238,\n",
      "                47.57271951827349\n",
      "            ],\n",
      "            [\n",
      "                8.391955910715069,\n",
      "                47.55694036558283\n",
      "            ],\n",
      "            [\n",
      "                8.381198082923012,\n",
      "                47.55004867353099\n",
      "            ],\n",
      "            [\n",
      "                8.378982594336376,\n",
      "                47.54322725120533\n",
      "            ],\n",
      "            [\n",
      "                8.36746618023274,\n",
      "                47.540332787307364\n",
      "            ],\n",
      "            [\n",
      "                8.366859793700819,\n",
      "                47.51720069783939\n",
      "            ],\n",
      "            [\n",
      "                8.340684897222106,\n",
      "                47.517196651399225\n",
      "            ],\n",
      "            [\n",
      "                8.337940470304243,\n",
      "                47.51394572267918\n",
      "            ],\n",
      "            [\n",
      "                8.344120279874552,\n",
      "                47.4595683837805\n",
      "            ],\n",
      "            [\n",
      "                8.356302754990763,\n",
      "                47.44516240382725\n",
      "            ],\n",
      "            [\n",
      "                8.353733194624247,\n",
      "                47.42610610482178\n",
      "            ],\n",
      "            [\n",
      "                8.360905670965298,\n",
      "                47.4259543246192\n",
      "            ],\n",
      "            [\n",
      "                8.363342285156245,\n",
      "                47.42266893070532\n",
      "            ],\n",
      "            [\n",
      "                8.361527469875877,\n",
      "                47.41860614546694\n",
      "            ],\n",
      "            [\n",
      "                8.347551069946636,\n",
      "                47.4158728732461\n",
      "            ],\n",
      "            [\n",
      "                8.34205917659486,\n",
      "                47.411225725227034\n",
      "            ],\n",
      "            [\n",
      "                8.341373663424736,\n",
      "                47.395428275896236\n",
      "            ],\n",
      "            [\n",
      "                8.34823771545445,\n",
      "                47.38706386570661\n",
      "            ],\n",
      "            [\n",
      "                8.37262998437049,\n",
      "                47.37943996442438\n",
      "            ],\n",
      "            [\n",
      "                8.375019280818673,\n",
      "                47.370329952877604\n",
      "            ],\n",
      "            [\n",
      "                8.384157885186566,\n",
      "                47.359752485081955\n",
      "            ],\n",
      "            [\n",
      "                8.383259013749992,\n",
      "                47.33825981656004\n",
      "            ],\n",
      "            [\n",
      "                8.39035173757096,\n",
      "                47.330285054010304\n",
      "            ],\n",
      "            [\n",
      "                8.391498803538775,\n",
      "                47.319202599966886\n",
      "            ],\n",
      "            [\n",
      "                8.399081519295152,\n",
      "                47.30948158182306\n",
      "            ],\n",
      "            [\n",
      "                8.390006280385407,\n",
      "                47.304916733281985\n",
      "            ],\n",
      "            [\n",
      "                8.377765433545175,\n",
      "                47.30478840431247\n",
      "            ],\n",
      "            [\n",
      "                8.370212543952075,\n",
      "                47.30107151897201\n",
      "            ],\n",
      "            [\n",
      "                8.368152856273147,\n",
      "                47.286202618787065\n",
      "            ],\n",
      "            [\n",
      "                8.36952604130393,\n",
      "                47.280627314550166\n",
      "            ],\n",
      "            [\n",
      "                8.373726248778953,\n",
      "                47.27922900257082\n",
      "            ],\n",
      "            [\n",
      "                8.375705631307687,\n",
      "                47.267559851397785\n",
      "            ],\n",
      "            [\n",
      "                8.379585622563491,\n",
      "                47.26373216605508\n",
      "            ],\n",
      "            [\n",
      "                8.382572417261335,\n",
      "                47.24561471902929\n",
      "            ],\n",
      "            [\n",
      "                8.392932141955834,\n",
      "                47.2358175878443\n",
      "            ],\n",
      "            [\n",
      "                8.398365214921863,\n",
      "                47.21713552479884\n",
      "            ],\n",
      "            [\n",
      "                8.405915938110713,\n",
      "                47.210133458965586\n",
      "            ],\n",
      "            [\n",
      "                8.437499999999986,\n",
      "                47.209280024465684\n",
      "            ],\n",
      "            [\n",
      "                8.4381882769779,\n",
      "                47.20453079823674\n",
      "            ],\n",
      "            [\n",
      "                8.447114061633252,\n",
      "                47.199862065707485\n",
      "            ],\n",
      "            [\n",
      "                8.454926582358215,\n",
      "                47.19962792461046\n",
      "            ],\n",
      "            [\n",
      "                8.461616873778944,\n",
      "                47.19425531484997\n",
      "            ],\n",
      "            [\n",
      "                8.49105619311432,\n",
      "                47.19425934936846\n",
      "            ],\n",
      "            [\n",
      "                8.517836068525108,\n",
      "                47.198461379333665\n",
      "            ],\n",
      "            [\n",
      "                8.524310362081467,\n",
      "                47.20258722555207\n",
      "            ],\n",
      "            [\n",
      "                8.55423435592164,\n",
      "                47.205387161343324\n",
      "            ],\n",
      "            [\n",
      "                8.574835175101269,\n",
      "                47.196147504244585\n",
      "            ],\n",
      "            [\n",
      "                8.583756949053017,\n",
      "                47.18865672305817\n",
      "            ],\n",
      "            [\n",
      "                8.595631044323278,\n",
      "                47.18685450211213\n",
      "            ],\n",
      "            [\n",
      "                8.605355224157181,\n",
      "                47.17352527011053\n",
      "            ],\n",
      "            [\n",
      "                8.60779210796711,\n",
      "                47.1615755379768\n",
      "            ],\n",
      "            [\n",
      "                8.611222944946643,\n",
      "                47.15877573561781\n",
      "            ],\n",
      "            [\n",
      "                8.638048150911477,\n",
      "                47.154408331590474\n",
      "            ],\n",
      "            [\n",
      "                8.653107713977,\n",
      "                47.14523612360122\n",
      "            ],\n",
      "            [\n",
      "                8.686751294812055,\n",
      "                47.14476923520714\n",
      "            ],\n",
      "            [\n",
      "                8.707350053100217,\n",
      "                47.14990485613047\n",
      "            ],\n",
      "            [\n",
      "                8.712841253312329,\n",
      "                47.154105579431295\n",
      "            ],\n",
      "            [\n",
      "                8.713527425203562,\n",
      "                47.1737107696219\n",
      "            ],\n",
      "            [\n",
      "                8.710239273173439,\n",
      "                47.18271277083933\n",
      "            ],\n",
      "            [\n",
      "                8.736186306551122,\n",
      "                47.19285678028009\n",
      "            ],\n",
      "            [\n",
      "                8.738496170162653,\n",
      "                47.218614640986175\n",
      "            ],\n",
      "            [\n",
      "                8.79936115809331,\n",
      "                47.21993826665251\n",
      "            ],\n",
      "            [\n",
      "                8.822368199536603,\n",
      "                47.233922410201224\n",
      "            ],\n",
      "            [\n",
      "                8.84880168370356,\n",
      "                47.22880914613987\n",
      "            ],\n",
      "            [\n",
      "                8.879698682507366,\n",
      "                47.229742922928004\n",
      "            ],\n",
      "            [\n",
      "                8.895225671260674,\n",
      "                47.238087363555294\n",
      "            ],\n",
      "            [\n",
      "                8.94561665125737,\n",
      "                47.244683351538264\n",
      "            ],\n",
      "            [\n",
      "                8.955227698326983,\n",
      "                47.24935115991394\n",
      "            ],\n",
      "            [\n",
      "                8.962093364251043,\n",
      "                47.25682090121835\n",
      "            ],\n",
      "            [\n",
      "                8.96895958852839,\n",
      "                47.267558949214106\n",
      "            ],\n",
      "            [\n",
      "                8.969565987549165,\n",
      "                47.27922900257082\n",
      "            ]\n",
      "        ]\n",
      "    ]\n",
      "}\n"
     ]
    }
   ],
   "source": [
    "print(\"Enter URL: E.g. https://wms.zh.ch/OGDOrthoZH?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&LAYERS=ortho&FORMAT=image/jpeg&CRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox}\")\n",
    "print(\"\")\n",
    "url = input()\n",
    "print(\"\")\n",
    "print(\"Enter Zoom level: E.g. 10\")\n",
    "print(\"\")\n",
    "zoom = int(input())\n",
    "\n",
    "def process_tile(tile, projs):\n",
    "    \"\"\" Extract geometry from area covered by tile \"\"\"\n",
    "    parameters = {}\n",
    "    if 'EPSG:4326' in projs:\n",
    "        parameters['proj'] = 'EPSG:4326'\n",
    "    elif 'EPSG:3857' in projs:\n",
    "        parameters['proj'] = 'EPSG:3857'\n",
    "    else:\n",
    "        for proj in projs:\n",
    "            try:\n",
    "                CRS.from_string(proj)\n",
    "            except Exception as e:\n",
    "                print(str(e))\n",
    "                continue\n",
    "            parameters['proj'] = proj\n",
    "            break\n",
    "    if 'proj' not in parameters:\n",
    "        raise RuntimeError(\"No support for available projections.\")\n",
    "\n",
    "    bounds = mercantile.bounds(*tile)\n",
    "    if parameters['proj'] == 'EPSG:4326' and '1.3.0' in url:\n",
    "        bbox = \",\".join(map(str, (bounds.south, bounds.west, bounds.north, bounds.east)))\n",
    "    elif parameters['proj'] == 'EPSG:4326':\n",
    "        bbox = \",\".join(map(str, (bounds.west, bounds.south, bounds.east, bounds.north)))\n",
    "    else:\n",
    "        proj_transformer = Transformer.from_crs(\"epsg:4326\", parameters['proj'], always_xy=True)\n",
    "        bl = proj_transformer.transform(bounds.west, bounds.south)\n",
    "        tr = proj_transformer.transform(bounds.east, bounds.north)\n",
    "        bbox = \",\".join(map(str, bl + tr))\n",
    "\n",
    "    parameters['bbox'] = bbox\n",
    "\n",
    "    parameters['width'] = 512\n",
    "    parameters['height'] = 512\n",
    "    query_url = url.format(**parameters)\n",
    "    print(\"Send request to\", query_url)\n",
    "\n",
    "    pixel_x = (bounds.east - bounds.west) / parameters['width']\n",
    "    pixel_y = (bounds.north - bounds.south) / parameters['height']\n",
    "    geotransform = (bounds.west, pixel_x, 0, bounds.north, 0, -pixel_y)\n",
    "\n",
    "    r = requests.get(query_url)\n",
    "    data = r.content\n",
    "    with MemoryFile(data) as memfile:\n",
    "        with memfile.open() as dataset:\n",
    "            data = np.zeros(shape=(dataset.height, dataset.width))\n",
    "            for band in range(dataset.count):\n",
    "                data += dataset.read(band + 1)\n",
    "\n",
    "            # Convert extrema to 0\n",
    "            data[data < 2] = 0\n",
    "            data[data > 254] = 0\n",
    "\n",
    "            # Blur data to avoid holes\n",
    "            data = gaussian_filter(data, sigma=3)\n",
    "\n",
    "            # Convert other pixels to 1\n",
    "            data[data > 0] = 1\n",
    "\n",
    "            # Filter small areas\n",
    "            data = sieve(data.astype(np.uint8),\n",
    "                         size=int(dataset.height * dataset.width * 0.001),\n",
    "                         connectivity=8)\n",
    "\n",
    "            # Create shapes for all areas with value = 1\n",
    "            shapes = list(features.shapes(data.astype(np.uint8), transform=Affine.from_gdal(*geotransform)))\n",
    "            buffer_dist = max(pixel_x, pixel_y) * 5\n",
    "            geoms = [shape(s[0]).buffer(buffer_dist).buffer(-buffer_dist)\n",
    "                     for s in shapes if int(s[1]) == 1]\n",
    "            return geoms\n",
    "\n",
    "\n",
    "def parse_wms(xml):\n",
    "    \"\"\" Rudimentary parsing of WMS Layers from GetCapabilites Request\n",
    "        owslib.wms seems to have problems parsing some weird not relevant metadata.\n",
    "        This function aims at only parsing relevant layer metadata\n",
    "    \"\"\"\n",
    "    wms = {}\n",
    "    # Remove prefixes to make parsing easier\n",
    "    # From https://stackoverflow.com/questions/13412496/python-elementtree-module-how-to-ignore-the-namespace-of-xml-files-to-locate-ma\n",
    "    try:\n",
    "        it = ET.iterparse(StringIO(xml))\n",
    "        for _, el in it:\n",
    "            _, _, el.tag = el.tag.rpartition('}')\n",
    "        root = it.root\n",
    "    except:\n",
    "        raise RuntimeError(\"Could not parse XML.\")\n",
    "\n",
    "    root_tag = root.tag.rpartition(\"}\")[-1]\n",
    "    if root_tag in {'ServiceExceptionReport', 'ServiceException'}:\n",
    "        raise RuntimeError(\"WMS service exception\")\n",
    "\n",
    "    if root_tag not in {'WMT_MS_Capabilities', 'WMS_Capabilities'}:\n",
    "        raise RuntimeError(\"No Capabilities Element present: Root tag: {}\".format(root_tag))\n",
    "\n",
    "    if 'version' not in root.attrib:\n",
    "        raise RuntimeError(\"WMS version cannot be identified.\")\n",
    "    version = root.attrib['version']\n",
    "    wms['version'] = version\n",
    "\n",
    "    layers = {}\n",
    "\n",
    "    def parse_layer(element, crs=set(), styles={}, bbox=None):\n",
    "        new_layer = {'CRS': crs,\n",
    "                     'Styles': {},\n",
    "                     'BBOX': bbox}\n",
    "        new_layer['Styles'].update(styles)\n",
    "        for tag in ['Name', 'Title', 'Abstract']:\n",
    "            e = element.find(\"./{}\".format(tag))\n",
    "            if e is not None:\n",
    "                new_layer[e.tag] = e.text\n",
    "        for tag in ['CRS', 'SRS']:\n",
    "            es = element.findall(\"./{}\".format(tag))\n",
    "            for e in es:\n",
    "                new_layer[\"CRS\"].add(e.text.upper())\n",
    "        for tag in ['Style']:\n",
    "            es = element.findall(\"./{}\".format(tag))\n",
    "            for e in es:\n",
    "                new_style = {}\n",
    "                for styletag in ['Title', 'Name']:\n",
    "                    el = e.find(\"./{}\".format(styletag))\n",
    "                    if el is not None:\n",
    "                        new_style[styletag] = el.text\n",
    "                new_layer[\"Styles\"][new_style['Name']] = new_style\n",
    "        # WMS Version 1.3.0\n",
    "        e = element.find(\"./EX_GeographicBoundingBox\")\n",
    "        if e is not None:\n",
    "            bbox = [float(e.find(\"./{}\".format(orient)).text.replace(',', '.'))\n",
    "                    for orient in ['westBoundLongitude',\n",
    "                                   'southBoundLatitude',\n",
    "                                   'eastBoundLongitude',\n",
    "                                   'northBoundLatitude']]\n",
    "            new_layer['BBOX'] = bbox\n",
    "        # WMS Version < 1.3.0\n",
    "        e = element.find(\"./LatLonBoundingBox\")\n",
    "        if e is not None:\n",
    "            bbox = [float(e.attrib[orient].replace(',', '.')) for orient in ['minx', 'miny', 'maxx', 'maxy']]\n",
    "            new_layer['BBOX'] = bbox\n",
    "\n",
    "        if 'Name' in new_layer:\n",
    "            layers[new_layer['Name']] = new_layer\n",
    "\n",
    "        for sl in element.findall(\"./Layer\"):\n",
    "            parse_layer(sl,\n",
    "                        new_layer['CRS'].copy(),\n",
    "                        new_layer['Styles'],\n",
    "                        new_layer['BBOX'])\n",
    "\n",
    "    # Find child layers. CRS and Styles are inherited from parent\n",
    "    top_layers = root.findall(\".//Capability/Layer\")\n",
    "    for top_layer in top_layers:\n",
    "        parse_layer(top_layer)\n",
    "\n",
    "    wms['layers'] = layers\n",
    "\n",
    "    # Parse formats\n",
    "    formats = []\n",
    "    for es in root.findall(\".//Capability/Request/GetMap/Format\"):\n",
    "        formats.append(es.text)\n",
    "    wms['formats'] = formats\n",
    "\n",
    "    # Parse access constraints and fees\n",
    "    constraints = []\n",
    "    for es in root.findall(\".//AccessConstraints\"):\n",
    "        constraints.append(es.text)\n",
    "    fees = []\n",
    "    for es in root.findall(\".//Fees\"):\n",
    "        fees.append(es.text)\n",
    "    wms['Fees'] = fees\n",
    "    wms['AccessConstraints'] = constraints\n",
    "\n",
    "    return wms\n",
    "\n",
    "\n",
    "def get_wms(wms_url):\n",
    "    wms_args = {}\n",
    "    u = urlparse(wms_url)\n",
    "    url_parts = list(u)\n",
    "    for k, v in parse_qsl(u.query):\n",
    "        wms_args[k.lower()] = v\n",
    "\n",
    "    def get_getcapabilitie_url(wms_version=None):\n",
    "\n",
    "        get_capabilities_args = {'service': 'WMS',\n",
    "                                 'request': 'GetCapabilities'}\n",
    "        if wms_version is not None:\n",
    "            get_capabilities_args['version'] = wms_version\n",
    "\n",
    "        # Keep extra arguments, such as map or key\n",
    "        for key in wms_args:\n",
    "            if key not in {'version', 'request', 'layers', 'bbox', 'width', 'height', 'format', 'crs', 'srs'}:\n",
    "                get_capabilities_args[key] = wms_args[key]\n",
    "\n",
    "        url_parts[4] = urlencode(list(get_capabilities_args.items()))\n",
    "        return urlunparse(url_parts)\n",
    "\n",
    "    wms = None\n",
    "    for wmsversion in ['1.3.0', '1.1.1', '1.1.0', '1.0.0', None]:\n",
    "        url = get_getcapabilitie_url(wms_version=wmsversion)\n",
    "        try:\n",
    "            print(\"Send request to: {}\".format(url))\n",
    "            r = requests.get(url)\n",
    "            xml = r.text\n",
    "            wms = parse_wms(xml)\n",
    "            break\n",
    "        except Exception as e:\n",
    "            print(str(e))\n",
    "    if not wms:\n",
    "        print(\"Could not contact WMS server\")\n",
    "    return wms\n",
    "\n",
    "\n",
    "def simplify(geom, distance=200):\n",
    "    \"\"\" Simplify geometry in epsg:3857\"\"\"\n",
    "    transformer = Transformer.from_crs(\"epsg:4326\", \"epsg:3857\")\n",
    "    transformer_back = Transformer.from_crs(\"epsg:3857\", \"epsg:4326\")\n",
    "    geom_3857 = transform(transformer.transform, geom)\n",
    "    geom_3857_simplified = geom_3857.simplify(distance, preserve_topology=False)\n",
    "    geom_simplified = transform(transformer_back.transform, geom_3857_simplified)\n",
    "    return geom_simplified\n",
    "\n",
    "\n",
    "def get_bounding_box(url):\n",
    "    wms = get_wms(url)\n",
    "    wms_args = {}\n",
    "    u = urlparse(url)\n",
    "    for k, v in parse_qsl(u.query):\n",
    "        wms_args[k.lower()] = v\n",
    "    layers = wms_args['layers'].split(',')\n",
    "    bbox = box(*wms['layers'][layers[0]]['BBOX'])\n",
    "    projs = wms['layers'][layers[0]]['CRS']\n",
    "    for layer in layers[1:]:\n",
    "        bbox = cascaded_union([bbox, wms['layers'][layer]['BBOX']])\n",
    "        projs = projs & wms['layers'][layer]['CRS']\n",
    "    return bbox, projs\n",
    "\n",
    "\n",
    "bbox, projs = get_bounding_box(url)\n",
    "\n",
    "geoms = []\n",
    "tiles = list(mercantile.tiles(west=bbox.bounds[0],\n",
    "                              south=bbox.bounds[1],\n",
    "                              east=bbox.bounds[2],\n",
    "                              north=bbox.bounds[3],\n",
    "                              zooms=zoom))\n",
    "\n",
    "for tile in tiles:\n",
    "    tile_geoms = process_tile(tile, projs)\n",
    "    if len(tile_geoms) > 0:\n",
    "        geoms.extend(tile_geoms)\n",
    "\n",
    "geom = cascaded_union(geoms)\n",
    "geom = simplify(geom)\n",
    "\n",
    "\n",
    "def plot_geometry(geom):\n",
    "    bounds = geom.bounds\n",
    "    fig, ax = plt.subplots()\n",
    "    if isinstance(geom, Polygon):\n",
    "        geom = [geom]\n",
    "    for g in geom:\n",
    "        patch = PolygonPatch(g, facecolor='#6699cc', edgecolor='#6699cc',\n",
    "                             alpha=0.5, zorder=2)\n",
    "        ax.add_patch(patch)\n",
    "    ax.set_xlim(bounds[0], bounds[2])\n",
    "    ax.set_ylim(bounds[1], bounds[3])\n",
    "    plt.show()\n",
    "\n",
    "print(\"Extracted the following geometry:\")\n",
    "plot_geometry(geom)\n",
    "\n",
    "print(\"GeoJSON output:\")\n",
    "print(json.dumps(mapping(geom), indent=4))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "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.8.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}