{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Example of DOV search methods for interpretations (informele hydrogeologische stratigrafie)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[![Binder](https://mybinder.org/badge_logo.svg)](https://mybinder.org/v2/gh/DOV-Vlaanderen/pydov/master?filepath=docs%2Fnotebooks%2Fsearch_informele_hydrostratigrafie.ipynb)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Use cases explained below\n",
    "* Get 'informele hydrogeologische stratigrafie' in a bounding box\n",
    "* Get 'informele hydrogeologische stratigrafie' with specific properties\n",
    "* Get 'informele hydrogeologische stratigrafie' in a bounding box based on specific properties\n",
    "* Select 'informele hydrogeologische stratigrafie' in a municipality and return date\n",
    "* Get 'informele hydrogeologische stratigrafie' based on fields not available in the standard output dataframe\n",
    "* Get 'informele hydrogeologische stratigrafie' data, returning fields not available in the standard output dataframe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import inspect, sys"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# check pydov path\n",
    "import pydov"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get information about the datatype 'Informele hydrogeologische stratigrafie'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pydov.search.interpretaties import InformeleHydrogeologischeStratigrafieSearch\n",
    "itp = InformeleHydrogeologischeStratigrafieSearch()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A description is provided for the 'Informele stratigrafie' datatype:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'Een informele hydrogeologische stratigrafie van een boring geeft, op basis van de (gecodeerde) lithologie, een indeling weer naar de al dan niet watervoerende eigenschappen van een bepaald beschreven diepte-interval. \\r\\nHet verschil met de hydrogeologische stratigrafie is dat voor de informele hydrogeologische stratigrafie geen algemene standaard gebruikt wordt.  Bij de hydrogeologische stratigrafie wordt de Hydrogeologische Codering van de Ondergrond van Vlaanderen (HCOV). Deze kan beschouwd worden als de officiele hydrogeologische codering voor het Vlaams Gewest.'"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "itp.get_description()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The different fields that are available for objects of the 'Informele hydrogeologische stratigrafie' datatype can be requested with the get_fields() method:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pkey_interpretatie\n",
      "Type_proef\n",
      "Proefnummer\n",
      "pkey_boring\n",
      "x\n",
      "y\n",
      "start_interpretatie_mtaw\n",
      "diepte_tot_m\n",
      "gemeente\n",
      "Auteurs\n",
      "Datum\n",
      "Opdrachten\n",
      "betrouwbaarheid_interpretatie\n",
      "Geldig_van\n",
      "Geldig_tot\n",
      "eerste_invoer\n",
      "geom\n",
      "diepte_laag_van\n",
      "diepte_laag_tot\n",
      "beschrijving\n"
     ]
    }
   ],
   "source": [
    "fields = itp.get_fields()\n",
    "\n",
    "# print available fields\n",
    "for f in fields.values():\n",
    "    print(f['name'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can get more information of a field by requesting it from the fields dictionary:\n",
    "* *name*: name of the field\n",
    "* *definition*: definition of this field\n",
    "* *cost*: currently this is either 1 or 10, depending on the datasource of the field. It is an indication of the expected time it will take to retrieve this field in the output dataframe.\n",
    "* *notnull*: whether the field is mandatory or not\n",
    "* *type*: datatype of the values of this field\n",
    "\n",
    "Alternatively, you can list all the fields and their details by inspecting the `get_fields()` output or the search instance itself in a notebook:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        \n",
       "                <div class=\"classname code\">\n",
       "                    pydov.search.interpretaties.InformeleHydrogeologischeStratigrafieSearch\n",
       "                </div>\n",
       "            <div class=\"description\"><p>Een informele hydrogeologische stratigrafie van een boring geeft, op basis van de (gecodeerde) lithologie, een indeling weer naar de al dan niet watervoerende eigenschappen van een bepaald beschreven diepte-interval. \r\n",
       "Het verschil met de hydrogeologische stratigrafie is dat voor de informele hydrogeologische stratigrafie geen algemene standaard gebruikt wordt.  Bij de hydrogeologische stratigrafie wordt de Hydrogeologische Codering van de Ondergrond van Vlaanderen (HCOV). Deze kan beschouwd worden als de officiele hydrogeologische codering voor het Vlaams Gewest.</p></div><div class=\"fields\"><div>\n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>pkey_interpretatie</b> - URL die verwijst naar de gegevens van deze informele hydrogeologische stratigrafie op de website. Voeg '.xml' toe om een XML voorstelling van deze gegevens te verkrijgen.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Type_proef</b> - Het type van de proef waarvan de informele hydrogeologische stratigrafie een interpretatie is. Bij een informele hydrogeologische stratigrafie zal dit altijd 'Boring' zijn.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li><li>codelist:</li>\n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        \n",
       "                <div class=\"classname code\">\n",
       "                    pydov.util.codelists.FeatureCatalogueValues\n",
       "                </div>\n",
       "            \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Boring</b> - Boring - <i>De informele hydrogeologische stratigrafie is een interpretatie van een boring.</i></p>\n",
       "            </div>\n",
       "        \n",
       "            </div>\n",
       "        </ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Proefnummer</b> - Het proefnummer van de boring waarvan de informele hydrogeologische stratigrafie een interpretatie is.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>pkey_boring</b> - URL die verwijst naar de gegevens van de boring op de website. Voeg '.xml' toe om een XML voorstelling van deze gegevens te verkrijgen.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>x</b> - De x-coördinaat van de boring in het Lambert72 coördinaatsysteem (in meter, EPSG:31370).</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>y</b> - De y-coördinaat van de boring in het Lambert72 coördinaatsysteem (in meter, EPSG:31370).</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>start_interpretatie_mtaw</b> - De hoogte van het aanvangspeil van de boring in het TAW stelsel (in meter).</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>diepte_tot_m</b> - De maximumdiepte van de informele hydrogeologische stratigrafie in meter. Dit is de 'tot' van de diepste geïnterpreteerde laag.</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>gemeente</b> - De gemeente waarin de boring van deze interpretatie gelegen is.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Auteurs</b> - De auteur of auteurs van de informele hydrogeologische stratigrafie.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Datum</b> - De datum waarop de informele hydrogeologische stratigrafie werd beschreven.</p><ul><li>type: <span class=\"code small\">date</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Opdrachten</b> - De opdracht(en) waaraan de informele hydrogeologische stratigrafie gekoppeld is.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>betrouwbaarheid_interpretatie</b> - De betrouwbaarheid van het geheel van de informele hydrogeologische stratigrafie (goed, twijfelachtig of onbekend). Deze wordt toegekend door een DOV-medewerker verantwoordelijk voor de ingave.</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Geldig_van</b> - Het tijdstip vanaf wanneer de informele hydrogeologische stratigrafie geldig is.</p><ul><li>type: <span class=\"code small\">date</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>Geldig_tot</b> - Het tijdstip tot wanneer de informele hydrogeologische stratigrafie geldig is.</p><ul><li>type: <span class=\"code small\">date</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>eerste_invoer</b> - Het tijdstip waarop deze interpretatie voor het eerst in DOV ingevoerd werd.</p><ul><li>type: <span class=\"code small\">datetime</span></li><li>notnull: <span class=\"code small\">True</span></li><li>query: <span class=\"code small\">True</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>geom</b> - None</p><ul><li>type: <span class=\"code small\">geometry</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">False</span></li><li>cost: <span class=\"code small\">1</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>diepte_laag_van</b> - Diepte van de bovenkant van de laag informele hydrostratigrafie in meter.</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">False</span></li><li>cost: <span class=\"code small\">10</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>diepte_laag_tot</b> - Diepte van de onderkant van de laag informele hydrostratigrafie in meter.</p><ul><li>type: <span class=\"code small\">float</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">False</span></li><li>cost: <span class=\"code small\">10</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        \n",
       "            <style type=\"text/css\">\n",
       "                div.pydov {\n",
       "                    background-color: rgba(170,170,170,0.05);\n",
       "                    padding: 10px;\n",
       "                    padding-left: 20px;\n",
       "                    border-left: 1px solid #fee439;\n",
       "                    border-radius: 10px;\n",
       "                    margin: 10px 0;\n",
       "                    max-height: 50vh;\n",
       "                    overflow: auto;\n",
       "                }\n",
       "\n",
       "                .code {\n",
       "                    font-family: monospace;\n",
       "                }\n",
       "\n",
       "                .small {\n",
       "                    font-size: 0.8rem;\n",
       "                }\n",
       "            </style>\n",
       "            <div class=\"pydov\">\n",
       "        <p><b>beschrijving</b> - Benoeming van de eenheid van de laag informele hydrostratigrafie in vrije tekst (onbeperkt in lengte).</p><ul><li>type: <span class=\"code small\">string</span></li><li>notnull: <span class=\"code small\">False</span></li><li>query: <span class=\"code small\">False</span></li><li>cost: <span class=\"code small\">10</span></li><li>multivalue: <span class=\"code small\">False</span></li></ul>\n",
       "            </div>\n",
       "        </div></p></div>\n",
       "            </div>\n",
       "        "
      ],
      "text/plain": [
       "<pydov.search.interpretaties.InformeleHydrogeologischeStratigrafieSearch at 0x7f9498f8ba10>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "itp"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Example use cases"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get 'informele hydrogeologische stratigrafie' in a bounding box"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Get data for all the 'informele hydrogeologisch stratigrafie' interpretations that are geographically located within the bounds of the specified box.\n",
    "\n",
    "The coordinates are in the Belgian Lambert72 (EPSG:31370) coordinate system and are given in the order of lower left x, lower left y, upper right x, upper right y."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n",
      "[000/001] .\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pkey_interpretatie</th>\n",
       "      <th>pkey_boring</th>\n",
       "      <th>betrouwbaarheid_interpretatie</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>start_interpretatie_mtaw</th>\n",
       "      <th>diepte_laag_van</th>\n",
       "      <th>diepte_laag_tot</th>\n",
       "      <th>beschrijving</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>QUARTAIR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>2.5</td>\n",
       "      <td>13.6</td>\n",
       "      <td>FORMATIE VAN MALDEGEM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>13.6</td>\n",
       "      <td>25.5</td>\n",
       "      <td>FORMATIE VAN MALDEGEM (Lid van Wemmel)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>25.5</td>\n",
       "      <td>33.2</td>\n",
       "      <td>FORMATIE VAN LEDE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>33.2</td>\n",
       "      <td>67.3</td>\n",
       "      <td>FORMATIE VAN BRUSSEL</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                  pkey_interpretatie  \\\n",
       "0  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "1  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "2  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "3  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "4  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "\n",
       "                                         pkey_boring  \\\n",
       "0  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "1  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "2  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "3  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "4  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "\n",
       "  betrouwbaarheid_interpretatie         x         y  start_interpretatie_mtaw  \\\n",
       "0                          goed  172148.5  182587.5                     10.94   \n",
       "1                          goed  172148.5  182587.5                     10.94   \n",
       "2                          goed  172148.5  182587.5                     10.94   \n",
       "3                          goed  172148.5  182587.5                     10.94   \n",
       "4                          goed  172148.5  182587.5                     10.94   \n",
       "\n",
       "   diepte_laag_van  diepte_laag_tot                            beschrijving  \n",
       "0              0.0              2.5                                QUARTAIR  \n",
       "1              2.5             13.6                   FORMATIE VAN MALDEGEM  \n",
       "2             13.6             25.5  FORMATIE VAN MALDEGEM (Lid van Wemmel)  \n",
       "3             25.5             33.2                       FORMATIE VAN LEDE  \n",
       "4             33.2             67.3                    FORMATIE VAN BRUSSEL  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from pydov.util.location import Within, Box\n",
    "\n",
    "df = itp.search(location=Within(Box(100000, 100000, 200000, 200000, epsg=31370)))\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The dataframe contains one 'informele stratigrafie' interpretation where three layers ('laag') were identified. The available data are flattened to represent unique attributes per row of the dataframe.\n",
    "\n",
    "Using the *pkey_interpretatie* field one can request the details of this interpretation in a webbrowser:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "https://www.dov.vlaanderen.be/data/interpretatie/1981-217464\n"
     ]
    }
   ],
   "source": [
    "for pkey_interpretatie in set(df.pkey_interpretatie):\n",
    "    print(pkey_interpretatie)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get 'informele hydrogeologische stratigrafie' with specific properties"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next to querying interpretations based on their geographic location within a bounding box, we can also search for interpretations matching a specific set of properties. For this we can build a query using a combination of the 'InformeleHydrogeologischeStratigrafie' fields and operators provided by the WFS protocol.\n",
    "\n",
    "A list of possible operators can be found below:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['PropertyIsBetween',\n",
       " 'PropertyIsEqualTo',\n",
       " 'PropertyIsGreaterThan',\n",
       " 'PropertyIsGreaterThanOrEqualTo',\n",
       " 'PropertyIsLessThan',\n",
       " 'PropertyIsLessThanOrEqualTo',\n",
       " 'PropertyIsLike',\n",
       " 'PropertyIsNotEqualTo',\n",
       " 'PropertyIsNull',\n",
       " 'SortProperty']"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "[i for i,j in inspect.getmembers(sys.modules['owslib.fes2'], inspect.isclass) if 'Property' in i]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this example we build a query using the *PropertyIsEqualTo* operator to find all interpretations that are within the community (gemeente) of 'Herstappe':"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n",
      "[000/001] c\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pkey_interpretatie</th>\n",
       "      <th>pkey_boring</th>\n",
       "      <th>betrouwbaarheid_interpretatie</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>start_interpretatie_mtaw</th>\n",
       "      <th>diepte_laag_van</th>\n",
       "      <th>diepte_laag_tot</th>\n",
       "      <th>beschrijving</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>QUARTAIR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>2.5</td>\n",
       "      <td>13.6</td>\n",
       "      <td>FORMATIE VAN MALDEGEM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>13.6</td>\n",
       "      <td>25.5</td>\n",
       "      <td>FORMATIE VAN MALDEGEM (Lid van Wemmel)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>25.5</td>\n",
       "      <td>33.2</td>\n",
       "      <td>FORMATIE VAN LEDE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>33.2</td>\n",
       "      <td>67.3</td>\n",
       "      <td>FORMATIE VAN BRUSSEL</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                  pkey_interpretatie  \\\n",
       "0  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "1  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "2  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "3  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "4  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "\n",
       "                                         pkey_boring  \\\n",
       "0  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "1  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "2  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "3  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "4  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "\n",
       "  betrouwbaarheid_interpretatie         x         y  start_interpretatie_mtaw  \\\n",
       "0                          goed  172148.5  182587.5                     10.94   \n",
       "1                          goed  172148.5  182587.5                     10.94   \n",
       "2                          goed  172148.5  182587.5                     10.94   \n",
       "3                          goed  172148.5  182587.5                     10.94   \n",
       "4                          goed  172148.5  182587.5                     10.94   \n",
       "\n",
       "   diepte_laag_van  diepte_laag_tot                            beschrijving  \n",
       "0              0.0              2.5                                QUARTAIR  \n",
       "1              2.5             13.6                   FORMATIE VAN MALDEGEM  \n",
       "2             13.6             25.5  FORMATIE VAN MALDEGEM (Lid van Wemmel)  \n",
       "3             25.5             33.2                       FORMATIE VAN LEDE  \n",
       "4             33.2             67.3                    FORMATIE VAN BRUSSEL  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from owslib.fes2 import PropertyIsEqualTo\n",
    "\n",
    "query = PropertyIsEqualTo(propertyname='gemeente',\n",
    "                          literal='Rotselaar')\n",
    "df = itp.search(query=query)\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Once again we can use the *pkey_interpretatie* as a permanent link to the information of these interpretations:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "https://www.dov.vlaanderen.be/data/interpretatie/1981-217464\n"
     ]
    }
   ],
   "source": [
    "for pkey_interpretatie in set(df.pkey_interpretatie):\n",
    "    print(pkey_interpretatie)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get 'informele hydrogeologische stratigrafie' in a bounding box based on specific properties"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can combine a query on attributes with a query on geographic location to get the interpretations within a bounding box that have specific properties.\n",
    "\n",
    "The following example requests the interpretations of boreholes only, within the given bounding box.\n",
    "\n",
    "(Note that the datatype of the *literal* parameter should be a string, regardless of the datatype of this field in the output dataframe.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n",
      "[000/001] c\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pkey_interpretatie</th>\n",
       "      <th>pkey_boring</th>\n",
       "      <th>betrouwbaarheid_interpretatie</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>start_interpretatie_mtaw</th>\n",
       "      <th>diepte_laag_van</th>\n",
       "      <th>diepte_laag_tot</th>\n",
       "      <th>beschrijving</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.5</td>\n",
       "      <td>QUARTAIR</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>2.5</td>\n",
       "      <td>13.6</td>\n",
       "      <td>FORMATIE VAN MALDEGEM</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>13.6</td>\n",
       "      <td>25.5</td>\n",
       "      <td>FORMATIE VAN MALDEGEM (Lid van Wemmel)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>25.5</td>\n",
       "      <td>33.2</td>\n",
       "      <td>FORMATIE VAN LEDE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>goed</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>33.2</td>\n",
       "      <td>67.3</td>\n",
       "      <td>FORMATIE VAN BRUSSEL</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                  pkey_interpretatie  \\\n",
       "0  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "1  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "2  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "3  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "4  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "\n",
       "                                         pkey_boring  \\\n",
       "0  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "1  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "2  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "3  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "4  https://www.dov.vlaanderen.be/data/boring/1981...   \n",
       "\n",
       "  betrouwbaarheid_interpretatie         x         y  start_interpretatie_mtaw  \\\n",
       "0                          goed  172148.5  182587.5                     10.94   \n",
       "1                          goed  172148.5  182587.5                     10.94   \n",
       "2                          goed  172148.5  182587.5                     10.94   \n",
       "3                          goed  172148.5  182587.5                     10.94   \n",
       "4                          goed  172148.5  182587.5                     10.94   \n",
       "\n",
       "   diepte_laag_van  diepte_laag_tot                            beschrijving  \n",
       "0              0.0              2.5                                QUARTAIR  \n",
       "1              2.5             13.6                   FORMATIE VAN MALDEGEM  \n",
       "2             13.6             25.5  FORMATIE VAN MALDEGEM (Lid van Wemmel)  \n",
       "3             25.5             33.2                       FORMATIE VAN LEDE  \n",
       "4             33.2             67.3                    FORMATIE VAN BRUSSEL  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from owslib.fes2 import PropertyIsEqualTo\n",
    "\n",
    "query = PropertyIsEqualTo(\n",
    "            propertyname='Type_proef',\n",
    "            literal='Boring')\n",
    "\n",
    "df = itp.search(\n",
    "    location=Within(Box(100000, 100000, 200000, 200000, epsg=31370)),\n",
    "    query=query\n",
    "    )\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can look at one of the interpretations in a webbrowser using its *pkey_interpretatie*:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "https://www.dov.vlaanderen.be/data/interpretatie/1981-217464\n"
     ]
    }
   ],
   "source": [
    "for pkey_interpretatie in set(df.pkey_interpretatie):\n",
    "    print(pkey_interpretatie)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Select 'informele hydrogeologische stratigrafie' in a municipality and return date"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can limit the columns in the output dataframe by specifying the *return_fields* parameter in our search.\n",
    "\n",
    "In this example we query all the 'informele stratigrafie' interpretations in the community of Rotselaar and return their date:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Datum</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1981-01-01</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Datum\n",
       "0  1981-01-01"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "query = PropertyIsEqualTo(propertyname='gemeente',\n",
    "                          literal='Rotselaar')\n",
    "df = itp.search(query=query,\n",
    "                return_fields=('Datum',))\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Datum</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>unique</th>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>top</th>\n",
       "      <td>1981-01-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>freq</th>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "             Datum\n",
       "count            1\n",
       "unique           1\n",
       "top     1981-01-01\n",
       "freq             1"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.describe()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get 'informele hydrogeologische stratigrafie' based on fields not available in the standard output dataframe"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To keep the output dataframe size acceptable, not all available WFS fields are included in the standard output. However, one can use this information to select interpretations as illustrated below.\n",
    "\n",
    "For example, make a selection of the interpretations in the municipality of Rotselaar, before 1/1/1990:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pkey_interpretatie</th>\n",
       "      <th>Datum</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>1981-01-01</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                  pkey_interpretatie       Datum\n",
       "0  https://www.dov.vlaanderen.be/data/interpretat...  1981-01-01"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from owslib.fes2 import And, PropertyIsEqualTo, PropertyIsLessThan\n",
    "\n",
    "query = And([PropertyIsEqualTo(propertyname='gemeente',\n",
    "                               literal='Rotselaar'),\n",
    "             PropertyIsLessThan(propertyname='Datum', \n",
    "                                 literal='1990-01-01')]\n",
    "            )\n",
    "df = itp.search(query=query,\n",
    "                return_fields=('pkey_interpretatie', 'Datum'))\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get 'informele hydrogeologische stratigrafie' data, returning fields not available in the standard output dataframe"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As denoted in the previous example, not all available fields are available in the default output frame to keep its size limited. However, you can request any available field by including it in the *return_fields* parameter of the search:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[000/001] .\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pkey_interpretatie</th>\n",
       "      <th>pkey_boring</th>\n",
       "      <th>x</th>\n",
       "      <th>y</th>\n",
       "      <th>start_interpretatie_mtaw</th>\n",
       "      <th>gemeente</th>\n",
       "      <th>Auteurs</th>\n",
       "      <th>Proefnummer</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>https://www.dov.vlaanderen.be/data/interpretat...</td>\n",
       "      <td>https://www.dov.vlaanderen.be/data/boring/1981...</td>\n",
       "      <td>172148.5</td>\n",
       "      <td>182587.5</td>\n",
       "      <td>10.94</td>\n",
       "      <td>Rotselaar</td>\n",
       "      <td>Onbekend - Vlaamse Maatschappij voor Watervoor...</td>\n",
       "      <td>B/3005-016</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                  pkey_interpretatie  \\\n",
       "0  https://www.dov.vlaanderen.be/data/interpretat...   \n",
       "\n",
       "                                         pkey_boring         x         y  \\\n",
       "0  https://www.dov.vlaanderen.be/data/boring/1981...  172148.5  182587.5   \n",
       "\n",
       "   start_interpretatie_mtaw   gemeente  \\\n",
       "0                     10.94  Rotselaar   \n",
       "\n",
       "                                             Auteurs Proefnummer  \n",
       "0  Onbekend - Vlaamse Maatschappij voor Watervoor...  B/3005-016  "
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "query = PropertyIsEqualTo(\n",
    "            propertyname='gemeente',\n",
    "            literal='Rotselaar')\n",
    "\n",
    "df = itp.search(query=query,\n",
    "                return_fields=('pkey_interpretatie', 'pkey_boring',\n",
    "                               'x', 'y', 'start_interpretatie_mtaw', 'gemeente', 'Auteurs', 'Proefnummer'))\n",
    "\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Visualize results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Using Folium, we can display the results of our search on a map."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "# import the necessary modules (not included in the requirements of pydov!)\n",
    "import folium\n",
    "from folium.plugins import MarkerCluster\n",
    "from pyproj import Transformer"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "# convert the coordinates to lat/lon for folium\n",
    "def convert_latlon(x1, y1):\n",
    "    transformer = Transformer.from_crs(\"epsg:31370\", \"epsg:4326\", always_xy=True)\n",
    "    x2,y2 = transformer.transform(x1, y1)\n",
    "    return x2, y2\n",
    "\n",
    "df['lon'], df['lat'] = zip(*map(convert_latlon, df['x'], df['y'])) \n",
    "# convert to list\n",
    "loclist = df[['lat', 'lon']].values.tolist()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><span style=\"color:#565656\">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc=\"&lt;!DOCTYPE html&gt;\n",
       "&lt;html&gt;\n",
       "&lt;head&gt;\n",
       "    \n",
       "    &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap-glyphicons.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;\n",
       "    \n",
       "            &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,\n",
       "                initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;\n",
       "            &lt;style&gt;\n",
       "                #map_d3f3d1ee3d71dd10bc2da953fa4e9030 {\n",
       "                    position: relative;\n",
       "                    width: 100.0%;\n",
       "                    height: 100.0%;\n",
       "                    left: 0.0%;\n",
       "                    top: 0.0%;\n",
       "                }\n",
       "                .leaflet-container { font-size: 1rem; }\n",
       "            &lt;/style&gt;\n",
       "\n",
       "            &lt;style&gt;html, body {\n",
       "                width: 100%;\n",
       "                height: 100%;\n",
       "                margin: 0;\n",
       "                padding: 0;\n",
       "            }\n",
       "            &lt;/style&gt;\n",
       "\n",
       "            &lt;style&gt;#map {\n",
       "                position:absolute;\n",
       "                top:0;\n",
       "                bottom:0;\n",
       "                right:0;\n",
       "                left:0;\n",
       "                }\n",
       "            &lt;/style&gt;\n",
       "\n",
       "            &lt;script&gt;\n",
       "                L_NO_TOUCH = false;\n",
       "                L_DISABLE_3D = false;\n",
       "            &lt;/script&gt;\n",
       "\n",
       "        \n",
       "    &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/leaflet.markercluster/1.1.0/leaflet.markercluster.js&quot;&gt;&lt;/script&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/leaflet.markercluster/1.1.0/MarkerCluster.css&quot;/&gt;\n",
       "    &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/leaflet.markercluster/1.1.0/MarkerCluster.Default.css&quot;/&gt;\n",
       "&lt;/head&gt;\n",
       "&lt;body&gt;\n",
       "    \n",
       "    \n",
       "            &lt;div class=&quot;folium-map&quot; id=&quot;map_d3f3d1ee3d71dd10bc2da953fa4e9030&quot; &gt;&lt;/div&gt;\n",
       "        \n",
       "&lt;/body&gt;\n",
       "&lt;script&gt;\n",
       "    \n",
       "    \n",
       "            var map_d3f3d1ee3d71dd10bc2da953fa4e9030 = L.map(\n",
       "                &quot;map_d3f3d1ee3d71dd10bc2da953fa4e9030&quot;,\n",
       "                {\n",
       "                    center: [50.953140668280824, 4.683962451078404],\n",
       "                    crs: L.CRS.EPSG3857,\n",
       "                    ...{\n",
       "  &quot;zoom&quot;: 12,\n",
       "  &quot;zoomControl&quot;: true,\n",
       "  &quot;preferCanvas&quot;: false,\n",
       "}\n",
       "\n",
       "                }\n",
       "            );\n",
       "\n",
       "            \n",
       "\n",
       "        \n",
       "    \n",
       "            var tile_layer_169c35779694a2c671f795cda2228f21 = L.tileLayer(\n",
       "                &quot;https://tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,\n",
       "                {\n",
       "  &quot;minZoom&quot;: 0,\n",
       "  &quot;maxZoom&quot;: 19,\n",
       "  &quot;maxNativeZoom&quot;: 19,\n",
       "  &quot;noWrap&quot;: false,\n",
       "  &quot;attribution&quot;: &quot;\\u0026copy; \\u003ca href=\\&quot;https://www.openstreetmap.org/copyright\\&quot;\\u003eOpenStreetMap\\u003c/a\\u003e contributors&quot;,\n",
       "  &quot;subdomains&quot;: &quot;abc&quot;,\n",
       "  &quot;detectRetina&quot;: false,\n",
       "  &quot;tms&quot;: false,\n",
       "  &quot;opacity&quot;: 1,\n",
       "}\n",
       "\n",
       "            );\n",
       "        \n",
       "    \n",
       "            tile_layer_169c35779694a2c671f795cda2228f21.addTo(map_d3f3d1ee3d71dd10bc2da953fa4e9030);\n",
       "        \n",
       "    \n",
       "            var marker_cluster_c33c0e13952f967a6379710ebb2c8933 = L.markerClusterGroup(\n",
       "                {\n",
       "}\n",
       "            );\n",
       "        \n",
       "    \n",
       "            var marker_2b16f847e5f81f97d4706ba614098f72 = L.marker(\n",
       "                [50.953140668280824, 4.683962451078404],\n",
       "                {\n",
       "}\n",
       "            ).addTo(marker_cluster_c33c0e13952f967a6379710ebb2c8933);\n",
       "        \n",
       "    \n",
       "        var popup_85c4b3532fdb2e10ed68b00a727f9050 = L.popup({\n",
       "  &quot;maxWidth&quot;: &quot;100%&quot;,\n",
       "});\n",
       "\n",
       "        \n",
       "            \n",
       "                var html_76e6a38d1d88bd22370f5562dbd582bf = $(`&lt;div id=&quot;html_76e6a38d1d88bd22370f5562dbd582bf&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;B/3005-016&lt;/div&gt;`)[0];\n",
       "                popup_85c4b3532fdb2e10ed68b00a727f9050.setContent(html_76e6a38d1d88bd22370f5562dbd582bf);\n",
       "            \n",
       "        \n",
       "\n",
       "        marker_2b16f847e5f81f97d4706ba614098f72.bindPopup(popup_85c4b3532fdb2e10ed68b00a727f9050)\n",
       "        ;\n",
       "\n",
       "        \n",
       "    \n",
       "    \n",
       "            marker_cluster_c33c0e13952f967a6379710ebb2c8933.addTo(map_d3f3d1ee3d71dd10bc2da953fa4e9030);\n",
       "        \n",
       "&lt;/script&gt;\n",
       "&lt;/html&gt;\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7f94982567b0>"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# initialize the Folium map on the centre of the selected locations, play with the zoom until ok\n",
    "fmap = folium.Map(location=[df['lat'].mean(), df['lon'].mean()], zoom_start=12)\n",
    "marker_cluster = MarkerCluster().add_to(fmap)\n",
    "for loc in range(0, len(loclist)):\n",
    "    folium.Marker(loclist[loc], popup=df['Proefnummer'][loc]).add_to(marker_cluster)\n",
    "fmap\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": ".venv (3.13.5)",
   "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.13.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}