plotting.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "e220c39c-9f68-40d8-a7e2-e48fd5134904",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import requests"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "institutional-caribbean",
   "metadata": {},
   "outputs": [],
   "source": [
    "import folium"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "43f6bbfd-713c-4922-8c85-522e39667322",
   "metadata": {},
   "outputs": [],
   "source": [
    "fname = '/mnt/lustre/datadive/Ecology/environmentaldata.csv'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "e7a69f11-0104-4a47-a850-99aa9eb931e2",
   "metadata": {},
   "outputs": [],
   "source": [
    "data = pd.read_csv(fname, low_memory=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "3658b0fc-9dca-44f9-ad1c-7590b50a4e73",
   "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>id</th>\n",
       "      <th>sample.samplingPoint</th>\n",
       "      <th>sample.samplingPoint.notation</th>\n",
       "      <th>sample.samplingPoint.label</th>\n",
       "      <th>sample.sampleDateTime</th>\n",
       "      <th>determinand.label</th>\n",
       "      <th>determinand.definition</th>\n",
       "      <th>determinand.notation</th>\n",
       "      <th>resultQualifier.notation</th>\n",
       "      <th>result</th>\n",
       "      <th>codedResultInterpretation.interpretation</th>\n",
       "      <th>determinand.unit.label</th>\n",
       "      <th>sample.sampledMaterialType.label</th>\n",
       "      <th>sample.isComplianceSample</th>\n",
       "      <th>sample.purpose.label</th>\n",
       "      <th>sample.samplingPoint.easting</th>\n",
       "      <th>sample.samplingPoint.northing</th>\n",
       "      <th>yr</th>\n",
       "      <th>month</th>\n",
       "      <th>day</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>http://environment.data.gov.uk/water-quality/d...</td>\n",
       "      <td>http://environment.data.gov.uk/water-quality/i...</td>\n",
       "      <td>SO-19801A08</td>\n",
       "      <td>APSLEY FARM BOREHOLE A</td>\n",
       "      <td>2000-03-22 10:30:00</td>\n",
       "      <td>BOD ATU</td>\n",
       "      <td>BOD : 5 Day ATU</td>\n",
       "      <td>85</td>\n",
       "      <td>&lt;</td>\n",
       "      <td>3.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>mg/l</td>\n",
       "      <td>UNCODED</td>\n",
       "      <td>False</td>\n",
       "      <td>WASTE MONITORING (AGENCY AUDIT - PERMIT)</td>\n",
       "      <td>442350.0</td>\n",
       "      <td>146570</td>\n",
       "      <td>2000</td>\n",
       "      <td>3</td>\n",
       "      <td>22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>http://environment.data.gov.uk/water-quality/d...</td>\n",
       "      <td>http://environment.data.gov.uk/water-quality/i...</td>\n",
       "      <td>SO-19801A08</td>\n",
       "      <td>APSLEY FARM BOREHOLE A</td>\n",
       "      <td>2000-03-22 10:30:00</td>\n",
       "      <td>Ammonia(N)</td>\n",
       "      <td>Ammoniacal Nitrogen as N</td>\n",
       "      <td>111</td>\n",
       "      <td>&lt;</td>\n",
       "      <td>0.50</td>\n",
       "      <td>NaN</td>\n",
       "      <td>mg/l</td>\n",
       "      <td>UNCODED</td>\n",
       "      <td>False</td>\n",
       "      <td>WASTE MONITORING (AGENCY AUDIT - PERMIT)</td>\n",
       "      <td>442350.0</td>\n",
       "      <td>146570</td>\n",
       "      <td>2000</td>\n",
       "      <td>3</td>\n",
       "      <td>22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>http://environment.data.gov.uk/water-quality/d...</td>\n",
       "      <td>http://environment.data.gov.uk/water-quality/i...</td>\n",
       "      <td>SO-19801A08</td>\n",
       "      <td>APSLEY FARM BOREHOLE A</td>\n",
       "      <td>2000-03-22 10:30:00</td>\n",
       "      <td>Chloride Ion</td>\n",
       "      <td>Chloride</td>\n",
       "      <td>172</td>\n",
       "      <td>NaN</td>\n",
       "      <td>14.10</td>\n",
       "      <td>NaN</td>\n",
       "      <td>mg/l</td>\n",
       "      <td>UNCODED</td>\n",
       "      <td>False</td>\n",
       "      <td>WASTE MONITORING (AGENCY AUDIT - PERMIT)</td>\n",
       "      <td>442350.0</td>\n",
       "      <td>146570</td>\n",
       "      <td>2000</td>\n",
       "      <td>3</td>\n",
       "      <td>22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>http://environment.data.gov.uk/water-quality/d...</td>\n",
       "      <td>http://environment.data.gov.uk/water-quality/i...</td>\n",
       "      <td>SO-19801A08</td>\n",
       "      <td>APSLEY FARM BOREHOLE A</td>\n",
       "      <td>2000-03-22 10:30:00</td>\n",
       "      <td>pH</td>\n",
       "      <td>pH</td>\n",
       "      <td>61</td>\n",
       "      <td>NaN</td>\n",
       "      <td>7.15</td>\n",
       "      <td>NaN</td>\n",
       "      <td>phunits</td>\n",
       "      <td>UNCODED</td>\n",
       "      <td>False</td>\n",
       "      <td>WASTE MONITORING (AGENCY AUDIT - PERMIT)</td>\n",
       "      <td>442350.0</td>\n",
       "      <td>146570</td>\n",
       "      <td>2000</td>\n",
       "      <td>3</td>\n",
       "      <td>22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>http://environment.data.gov.uk/water-quality/d...</td>\n",
       "      <td>http://environment.data.gov.uk/water-quality/i...</td>\n",
       "      <td>SO-19801A08</td>\n",
       "      <td>APSLEY FARM BOREHOLE A</td>\n",
       "      <td>2000-03-22 10:30:00</td>\n",
       "      <td>COD as O2</td>\n",
       "      <td>Chemical Oxygen Demand :- {COD}</td>\n",
       "      <td>92</td>\n",
       "      <td>NaN</td>\n",
       "      <td>27.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>mg/l</td>\n",
       "      <td>UNCODED</td>\n",
       "      <td>False</td>\n",
       "      <td>WASTE MONITORING (AGENCY AUDIT - PERMIT)</td>\n",
       "      <td>442350.0</td>\n",
       "      <td>146570</td>\n",
       "      <td>2000</td>\n",
       "      <td>3</td>\n",
       "      <td>22</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                  id  \\\n",
       "0  http://environment.data.gov.uk/water-quality/d...   \n",
       "1  http://environment.data.gov.uk/water-quality/d...   \n",
       "2  http://environment.data.gov.uk/water-quality/d...   \n",
       "3  http://environment.data.gov.uk/water-quality/d...   \n",
       "4  http://environment.data.gov.uk/water-quality/d...   \n",
       "\n",
       "                                sample.samplingPoint  \\\n",
       "0  http://environment.data.gov.uk/water-quality/i...   \n",
       "1  http://environment.data.gov.uk/water-quality/i...   \n",
       "2  http://environment.data.gov.uk/water-quality/i...   \n",
       "3  http://environment.data.gov.uk/water-quality/i...   \n",
       "4  http://environment.data.gov.uk/water-quality/i...   \n",
       "\n",
       "  sample.samplingPoint.notation sample.samplingPoint.label  \\\n",
       "0                   SO-19801A08     APSLEY FARM BOREHOLE A   \n",
       "1                   SO-19801A08     APSLEY FARM BOREHOLE A   \n",
       "2                   SO-19801A08     APSLEY FARM BOREHOLE A   \n",
       "3                   SO-19801A08     APSLEY FARM BOREHOLE A   \n",
       "4                   SO-19801A08     APSLEY FARM BOREHOLE A   \n",
       "\n",
       "  sample.sampleDateTime determinand.label           determinand.definition  \\\n",
       "0   2000-03-22 10:30:00           BOD ATU                  BOD : 5 Day ATU   \n",
       "1   2000-03-22 10:30:00        Ammonia(N)         Ammoniacal Nitrogen as N   \n",
       "2   2000-03-22 10:30:00      Chloride Ion                         Chloride   \n",
       "3   2000-03-22 10:30:00                pH                               pH   \n",
       "4   2000-03-22 10:30:00         COD as O2  Chemical Oxygen Demand :- {COD}   \n",
       "\n",
       "   determinand.notation resultQualifier.notation  result  \\\n",
       "0                    85                        <    3.00   \n",
       "1                   111                        <    0.50   \n",
       "2                   172                      NaN   14.10   \n",
       "3                    61                      NaN    7.15   \n",
       "4                    92                      NaN   27.00   \n",
       "\n",
       "   codedResultInterpretation.interpretation determinand.unit.label  \\\n",
       "0                                       NaN                   mg/l   \n",
       "1                                       NaN                   mg/l   \n",
       "2                                       NaN                   mg/l   \n",
       "3                                       NaN                phunits   \n",
       "4                                       NaN                   mg/l   \n",
       "\n",
       "  sample.sampledMaterialType.label  sample.isComplianceSample  \\\n",
       "0                          UNCODED                      False   \n",
       "1                          UNCODED                      False   \n",
       "2                          UNCODED                      False   \n",
       "3                          UNCODED                      False   \n",
       "4                          UNCODED                      False   \n",
       "\n",
       "                       sample.purpose.label  sample.samplingPoint.easting  \\\n",
       "0  WASTE MONITORING (AGENCY AUDIT - PERMIT)                      442350.0   \n",
       "1  WASTE MONITORING (AGENCY AUDIT - PERMIT)                      442350.0   \n",
       "2  WASTE MONITORING (AGENCY AUDIT - PERMIT)                      442350.0   \n",
       "3  WASTE MONITORING (AGENCY AUDIT - PERMIT)                      442350.0   \n",
       "4  WASTE MONITORING (AGENCY AUDIT - PERMIT)                      442350.0   \n",
       "\n",
       "   sample.samplingPoint.northing    yr  month  day  \n",
       "0                         146570  2000      3   22  \n",
       "1                         146570  2000      3   22  \n",
       "2                         146570  2000      3   22  \n",
       "3                         146570  2000      3   22  \n",
       "4                         146570  2000      3   22  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "42d44dd6-013e-4899-8731-2ca40291cb59",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['id', 'sample.samplingPoint', 'sample.samplingPoint.notation',\n",
       "       'sample.samplingPoint.label', 'sample.sampleDateTime',\n",
       "       'determinand.label', 'determinand.definition', 'determinand.notation',\n",
       "       'resultQualifier.notation', 'result',\n",
       "       'codedResultInterpretation.interpretation', 'determinand.unit.label',\n",
       "       'sample.sampledMaterialType.label', 'sample.isComplianceSample',\n",
       "       'sample.purpose.label', 'sample.samplingPoint.easting',\n",
       "       'sample.samplingPoint.northing', 'yr', 'month', 'day'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "95ffceb4-c64d-49d7-af86-f89a5e273f3d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "486176"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(data['sample.sampleDateTime'].unique())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "fb9c06d4-829e-4f75-bbb7-94615c7049ff",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'http://environment.data.gov.uk/water-quality/id/sampling-point/SO-19801A08'"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.iloc[0]['sample.samplingPoint']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "1a937afe-e74c-4348-9976-1df83d162508",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'http://environment.data.gov.uk/water-quality/data/measurement/SO-617981-0085'"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.iloc[0].id"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "d4aa3630-9721-44e1-a391-db3834464bc3",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "51.216762 -1.395023\n"
     ]
    }
   ],
   "source": [
    "req = requests.get(data.iloc[0]['sample.samplingPoint'])\n",
    "sample_data = req.json()\n",
    "lat = sample_data['items'][0]['lat']\n",
    "long = sample_data['items'][0]['long']\n",
    "\n",
    "print(lat, long)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "48dff5cc-015c-47c9-ad3d-5e6b814d92b0",
   "metadata": {},
   "outputs": [],
   "source": [
    "from shapely.geometry import Point\n",
    "import geopandas as gpd\n",
    "from geopandas import GeoDataFrame"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "74a294ce-ab35-40b8-a9d3-74f48e25287a",
   "metadata": {},
   "outputs": [],
   "source": [
    "geometry = [Point(xy) for xy in zip(data['sample.samplingPoint.easting'][:10], data['sample.samplingPoint.northing'][:10])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "b2277a5e-4018-450a-ba4a-5a60d6bbd1cc",
   "metadata": {},
   "outputs": [],
   "source": [
    "# The world low-res maps are in EPSG:4236 (lat, long coords)\n",
    "# But the UK regions are in British National Grid\n",
    "# So we'll hold onto both versions\n",
    "gdf = GeoDataFrame(data[:10], geometry=geometry, crs=27700) \n",
    "gdf_4326 = gdf.to_crs(4326)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "referenced-appendix",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Shapefile from ONS\n",
    "# https://geoportal.statistics.gov.uk/datasets/ons::regions-december-2021-en-bfc\n",
    "uk_regions = gpd.read_file('resources/RGN_DEC_2021_EN_BFC.shp')\n",
    "uk_SE = uk_regions[(uk_regions['RGN21NM'] == \"South East\")].geometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "id": "fatal-wellington",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALMAAAD4CAYAAACni9dcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAn+klEQVR4nO2de1yN2f7HP0u5RCO0K100ZVAUpZLTcZvqGA4ZY4aMcZ3jlnvuzJifgxm5xmE4gyHTmHFpXAYNGSrk0lSUSqFIJIpuRGrv/f39oTqiXTv78jw9+3m/XuvV3s9lrc+Tj9V3XZ61GBFBREQINOBagIiIuhDNLCIYRDOLCAbRzCKCQTSziGDQ12ZhEomEbGxstFmkiMCIi4t7TEQm1Z3TqpltbGwQGxurzSJFBAZj7K6ic2KYISIYRDOLCAbRzCKCQTSziGAQzSwiGEQziwgG0cwigkE0czUQER4/fgxtTI8tLS1FRkYGYmJiIJPJNF6ekNHqoEl9YeHChdi0aRMMDAxgb2+Pjh07Vvlpa2sLff3af3VlZWV48OAB7t27h/v37+PevXtVPt+/fx9PnjyBoaEhioqKcPfuXVhaWmrhCYWJaOY3SE5ORlBQEDIzM9GgQQOkpKQgJSUFqampiIiIQGpqKh4+fIgPPvigirnz8vLeMurjx49hZmaGNm3awMrKCm3atIGtrS169+5d+d3MzAzjx4+HRCIRjawiTJtvmri5uRGfh7OJCF5eXhg6dCimTZum8Lrnz5/j5s2blUbPyMiAiYlJpUGtrKxgZWWF1q1b11qDR0REYOzYsbh+/ToMDQ3V/UiCgzEWR0Ru1Z0Ta+bX2LdvHwoKCuDn51fjdU2bNoWzszOcnZ1VKu/ly5eYMmUKNm3aJBpZDYhmLufp06eYP38+Dhw4AD09Pa2UuXbtWtjZ2eGTTz7RSnlCRzRzOcuWLUPfvn3x97//XSvlpaWlYePGjYiLi9NKebqAaGa8avT99NNPSE5O1kp5RIRp06Zh4cKFeP/997VSpi6g8/3MRITp06dj6dKlMDU11UqZBw4cQHZ2Nvz9/bVSnq6g8zWzso0+dVFQUIA5c+YgJCQEDRs21EqZuoJOm/n1Rp8ygyDqYMmSJfDx8dFabK5L6LSZtd3oi4mJwcGDB7UWm+saOmtmbTf6pFIpJk+ejDVr1qBVq1ZaKVPXUMrMjLEMAE8ByABIiciNMeYM4AcATQBIAUwlor80pFOtcNHo27JlC1q0aIFRo0ZppTydhIhqTQAyAEjeOHYKwD/LPw8AEFlbPq6ursQHfv31V3J2dqaysjKtlHfv3j0yNjam1NRUrZQnZADEkgJ/qRJmEIDm5Z+NADxQIS+tUVRUhHnz5iEkJERrjT5/f39MmzYNdnZ2WilPV1H2X5MAnGKMEYBtRLQdgD+AMMbYOrzqr662FcUYmwRgEgBYW1urLFhVli9fjo8++khrjb7Q0FDEx8djz549WilPp1FUZVPVkMKy/KcpgAQAvQFsAvBZ+XFfAKdry4frMCMpKYkkEgk9evRIK+UVFxeTjY0NnTp1Sivl6QKoIcxQagSQiLLKf+YAOAzAHcBYAIfKLwkpP8ZbiING34oVK+Dh4YG+fftqpTxdp9YwgzHWDEADInpa/vkjAMvxKkbuAyASgBeAWxrUqTKHDh1Cfn6+1kb6kpKS8OOPPyIxMVEr5YkoFzObATjMGKu4/lciOskYewbgP4wxfQAlKI+L+cqxY8fg5+enlUafXC7HlClTsHz5crRu3Vrj5Ym8otZ/WSK6DcCpmuNRAFw1IUoTXL58GbNnz9ZKWbt370ZpaSkmTeL1/2/BoRMjgHl5eXjw4AEcHR01XlZubi4WL16MkydP1jjJv6ysDImJiYiOjkZ0dDSys7PRoEEDpRJjrMr3Jk2aYM2aNTr/topOmDk6Ohpubm5aeYNkwYIFGDlyJLp27Vp5jIiQmZlZadzo6GjEx8fDxsYG7u7u8PDwwPvvvw8iglwur1OSSqXw8/NDQECAxp+N7+iEmS9fvoy//e1vGi/n7NmzOHPmDC5duoQzZ85UMS9jDN27d0f37t2xfPlyuLm5oXnz5rVnWgvJycn44IMPYGRkpIYnqN/ojJmnT5+u8XJmzpyJJ0+ewM7ODl27dkX37t0xatQobN68GW3atEF5I1qtxMXFwdW13jRdNIrgzSyXyxEdHY2ff/5Z42UFBATA3Nwcjo6OWpt4L5r5fwjezKmpqTA2NtbKQMmAAQM0XsabxMbGYsiQIVovl48I/h3Ay5cvw8PDg2sZGkEmkyEhIaFKY1OXEbyZL126pJXGHxekpqbCwsJCbPyVI3gzC7lmjo2NhZtbtStV6SSCNnNhYSHu3LmDLl26cC1FI4iNv6oI2swxMTHo2rWrYF/pF81cFUGbWcghhlQqRUJCAlxcXLiWwhsEbWahN/4sLS3VMoooFARrZiLS2jA2F4ghxtsI1sxpaWkwNDSEhYUF11I0gtiT8TaCNbOQQwxArJmrQ7BmFnKIIZVKce3aNXHk7w0EbWah9mSkpKSIjb9qEKSZi4uLcePGDcHWXHFxcWK8XA2CNHNsbCw6d+6Mxo0bcy1FIxw5ckRcErcaBDkFVMghxrFjx3D9+nXs27ePaym8Q5Bmzs7OFuQQ9rNnzzB9+nQEBQWhSZMmXMvhHYIMM+bNm4egoCDcuHGDaylq5ZtvvoGnpye8vLy4lsJLBFkzW1lZYcmSJZg6dSpOnz6tkXfvtE1cXBz27t2LpKQkrqXwFkHWzAAwffp05OfnC2L1TalUiokTJ2LNmjWQSCRcy+EtgjWzvr4+tm3bhgULFiAvL49rOSqxefNmtGzZEqNHj+ZaCq8RrJkBoFu3bhg6dCgWLlzItZR3JjMzE9999x1++OEHQYRLmkTQZgaAb7/9Fn/88QcuXLjAtZQ6Q+U7ufr7+6N9+/Zcy+E9gjezkZERNmzYAD8/P5SVlXEtp04cPHgQt2/fxoIFC7iWUi8QvJkBYNiwYbCyskJgYCDXUpSmsLAQs2bNwrZt29CoUSOu5dQLdMLMjDFs2bIFa9euxZ07d7iWoxSLFy+Gj48PevbsybWUeoMg+5mro23btpg7dy6mT5+O48eP87oxdenSJRw5ckTcybWO6ETNXMHcuXORkZGBw4cPcy1FIWVlZZg0aRICAwPRsmVLruXUK3TKzI0aNcKkSZMQHh7OtRSFrF+/HlZWVhg+fDjXUuodOhNmVGBoaIgXL15wLaNa0tPTsW7dOsTExPA6DOIrOlUzA0DTpk3x/PlzrmW8BRFhypQpWLhwIWxtbbmWUy9RysyMsQzGWCJjLJ4xFvva8RmMsVTGWDJjbI3mZKoPvpr5119/xaNHj+Dv78+1lHpLXcIMTyJ6XPGFMeYJYDAAJyJ6yRjTzk6RKsJHM+fl5WHevHn4/fffBTkPW1uoEmZMAbCKiF4Clbu38h4DAwPexcwLFizA0KFD4e7O601ueY+yZq7YCD6ufGN3AOgAoBdjLJoxdpYx1q26GxljkxhjsYyx2NzcXHVoVgm+1cxxcXHYuXMnDAwMsHPnTvzxxx+4evUqHj58CJlMxrW8eoWyYUZPIsoqDyX+ZIyllt/bCsDfAHQDcIAx1rZ8s+5KiGg7gO0A4ObmRuAYvpnZxsYGmzdvRnZ2NqKiopCdnV2ZCgoKYGJiAnNzc5ibm8PCwqLyc8V3FxcXNGigc+34alHKzK9vBM8Yq9gI/j6AQ+Xm/YsxJgcgAcB99VsDfDOzsbGxwp2wSktL8ejRI2RnZ+PBgweVJo+NjcWDBw+QnJyML7/8EkuXLtWyan6iykbwzwB4AohgjHUA0AjAY8U58QM+xsyKaNSoEdq0aYM2bdpUe/7mzZvo0aMH5s6dq/O7swLKxcxmAKIYYwkA/gIQSkQnAewC0JYxlgRgH4Cxb4YYfESZmjklJQX14FHQoUMHeHp6Yvv27VxL4QdEpLXk6upKXCOTyYgxRnK5vNrz4eHhxBijzZs3a1nZu3HlyhWytLSkkpISrqVoBQCxpMBfOtdyaNCgARo3boySkpK3zuXl5WHs2LH4/vvvsXz5csTFxXGgsG507doVnTt31sqmnbxHkcs1kfhQMxMRtWzZkp48eVLlmFwup88++4z8/f2JiGj//v3Utm1bKigo4EJinTh79iy1a9eOpFIp11I0DsSauSrVxc1BQUG4desWAgICAAC+vr7o168fJk6cyPv4uVevXjA1NcVvv/3GtRRuUeRyTSS+1Mzt27enGzduVH6/efMmSSQSSkpKqnLdixcvyNnZmbZs2aJtiXXm2LFj5OTkpLAtIBQg1sxVeb1mLisrwxdffIF///vfcHBwqHJdkyZNcODAASxduhRXr17lQqrSDBw4EHK5HCdPnuRaCmfopJlf72teunQpzMzMMHXq1Gqvbd++PWbPno2VK1dqU2KdYYxh0aJFlWGSLqKTZq6omSMjI7F7927s2rVL4WR4IsJvv/2GESNGaFll3fH19UVWVla9XCNEHeismbOysjBmzBjs3LkTpqaKZ68eOXIERIRPPvlEewLfEX19fcyfP193a2dFwbQmEl8agL6+vmRmZkYzZ86s8TqZTEadO3emo0ePakmZ6rx48YLMzc0pISGBaykaAWIDsCpNmzaFiYkJVq9eXeN1Bw8eRJMmTeDj46MlZarTpEkT+Pv7Y9WqVVxL0T6KXK6JxJea+erVq3Tnzp0ar5FKpdSxY0c6ceKEdkSpkcLCQjI2Nqa0tDSupagdiDVzVZydnWFjY1PjNfv370eLFi3Qr18/7YhSI82bN4efnx/Wrl3LtRStwkiLo1tubm4UGxtb+4WvQUR4+vQpHj9+XJmkUil8fHw0NildKpXCwcEBW7duhbe3t0bK0DS5ubmws7NDcnIyzM3NuZajNhhjcURU7b5xvFg3Y//+/bh+/XoVw+bm5lZ+bty4MSQSCSQSCUxMTHD//n0cOHAAQUFBGnkB9JdffkHr1q3r9d4hJiYmGDVqFDZs2IA1a+rFi/Mqw4uauX379ujduze6du1axbQSiQTGxsZv7az0/Plz+Pr6gogQEhKCpk2bqk1jWVkZ7O3tsWvXLvTp00dt+XJBZmYmunbtirS0NMEs9VVTzcyLBuC4ceNo69atdWoIlJaW0qhRo6hHjx6Un59fp3tr4vfffye+NFTVwdixY2nFihVcy1Ab4HsDsHfv3jh//nyd7mnYsCF++uknuLm5oU+fPsjOzlaLFg8PD2RkZCAtLU0t+XHNwoULsXnzZl6996gpeGPmc+fOgeoY8jRo0AAbNmyAr68vevbsifT0dJW1mJiYYM6cOVi8eLHKefGBjh07okePHvjxxx+5lqJ5FFXZmkiK/nzL5XIyNzen9PT0d/7z88MPP5CFhQXFx8e/cx4VFBcXk5WVFV28eFHlvPjAX3/9RVZWVvTy5UuupagM+B5mMMbeKdR4ncmTJ2Pjxo346KOPEBUVpZKepk2bYsWKFZg7d26d/1rwkW7dusHe3h6//PIL11I0Ci/MDPwv1FCFYcOGYc+ePfj0008RGhqqUl6jR49GcXExDh06pFI+fGHx4sVYvXq1sFdJUlRlayLV1Etw7do1ateunVr+FEVHR5OZmRkFBwerlE9YWBi1adOG1q1bR5cvX67Xf6blcjm5u7vTb7/9xrUUlUANYQZvzCyTyahly5b04MEDtTz09evXqU2bNrRx40aV8jl06BBNnTqVunTpQs2aNaPevXvT4sWLKTQ0lPLy8tSiVVscPnyYXF1d6/WrVfXCzEREH3/8Me3fv189T01Ed+/eJTs7O1qyZIla/gELCgro5MmT9M0335CnpycZGhqSg4MDTZ48mYKDgyk9PZ3XRpHJZNSpUyc6deoU11LemZrMzIsRwArWrVuHjIwMfP/992orMzc3F//85z/h5uaGLVu2QE9PT215l5WVISEhARcuXMCFCxcQFRUFIkLPnj3Ro0cP9OjRA87Ozrxaczk4OBi7d+/m9b4uNcH7EcAKoqOjqXPnzmr4/1uVoqIi8vLyomHDhml05R+5XE63b9+mn3/+mfz8/Khz585kaGhInp6edPXqVY2VWxdKS0vJ1taWRo8eTRERESSTybiWVCdQX8KM7OxsatGihXqe+g1evHhBn376KfXt25eKi4s1UkZ15OfnU2BgINnb29Pz58+1Vm5N5OTkUGBgIHXu3JlsbGxo6dKlKvXxa5N6Y+ZDhw5R//791fPU1VBSUkJGRkaUmJiosTIU4evrS3PmzNF6uTUhl8vpypUrNHPmTJJIJNSnTx8KCgqip0+fci1NIfXGzP7+/rRy5Ur1PHU1/PTTT9SnTx+N5V8Tjx8/JgsLCzp79iwn5dfGy5cv6dChQ/Txxx+TkZERjRkzhi5cuMC1rLeoN2Z2dXWl8+fPq+ep30AqlVKHDh3o9OnTGslfGY4ePUq2trZUVFTEmQZlePjwIQUGBpKRkRFlZWVxLacKNZmZNyOARUVFSE1NRbdu1W6NojIhISGQSCScTrgfNGgQPvzwQ8yfP58zDcpgZmaG2bNnY/DgwfVr/TpFLtdEqqlmPnHihMZCAJlMRg4ODrx4ObWgoICsra15oaU2jh8/Tj169OBaRhVQH2rm8+fPo1evXhrJ+8iRIzAwMODFy6lGRkbYtWsXJkyYgPz8fK7l1Ejfvn2RkpKCe/fucS1FKXhj5nPnzqF3795qz5eI8O233+Kbb77hzX7U3t7eGDJkCGbOnMm1lBpp1KgRBg8ejJCQEK6lKAUvzFxSUoKrV6/Cw8ND7XmHhoZCLpdj0KBBas9bFVatWoXo6Gjez8obPnw4Dhw4wLUM5VAUf2giKYqZz549S926dVNrbEX0v5liISEhas9bHVy4cIHMzMzo0aNHXEtRSGlpKUkkkloXzdEW4HvMrKkQ49SpU3j27Bk+/fRTteetDv7+979j3LhxmDx58qt+Uh7SsGFDDBkypF6EGkqZmTGWwRhLZIzFM8Zi3zg3lzFGjDHJu4o4d+4cXF1dkZubi8zMTNy4cQPx8fG4dOkSwsPDERoaijNnztRpYjkRYcWKFfj66695vYPpsmXLkJaWhj179nAtRSHDhw/H/v37uZZRK3VZBMaTiKpsWskYa4NXm1xmqiIiKysL48ePR9OmTWFgYFBtysnJwZMnTzBjxgyMHz8eRkZGNeYZGRmJnJwcDB8+XBVpGqdx48YIDg5Gv3794OnpCSsrK64lvUWfPn1w7949pKWloV27dlzLUYyi+OP1BCADgKSa478BcFJ0/s2k6noUly9fphEjRlDLli1p+vTpdPPmTYXXenp6UlBQkErlaZPly5dT3759eTsfesqUKRqdaqAsUHU4G8AdAFcAxAGYVH5sMID/UA1mLz83CUAsgFhra2u1PND9+/fpq6++IhMTExo4cCCdOnWqigmioqLI1taWSktL1VKeNnj58iXp6+vTs2fPuJZSLZGRkeTk5MS1DLWY2bL8pymABAC9AUQDMKJazPx6UvdKQc+fP6cdO3aQo6MjOTg40LZt26i4uJj69etH27ZtU2tZmiYjI4MsLCy4lqEQqVRK5ubmlJqayqmOmsysVMuIiLLKf+YAOAygDwBbAAmMsQwAVgCuMMZa1y3IUQ0DAwNMmDAB165dw3/+8x8cP34czZo1Q1hYGMaOHatNKSqTmpoKe3t7rmUoRE9PD0OHDuV1n3OtZmaMNWOMvVfxGa8afDFEZEpENkRkA+A+ABcieqhRtYo1wtvbG0ePHsWtW7dQUlKCxo0bcyHlnUlJSUHHjh25llEjfO/VUKZmNgMQxRhLAPAXgFAi4u1mc+3atat3Rgbqh5k9PDxQUFCA5ORkrqVUS61mJqLbRORUnhyI6LtqrrGhN7rtROoG38MM4NXafr6+vrwNNfg7mqBj1IeaGXi11+D+/ft5OWIpmpkHPHnyBC9fvqwX2zV0794dJSUlSExM5FrKW4hm5gEVIQZfpqjWBGOssnbmG6KZeUBqamq9CDEq4GuoIZq5DhARsrOzUVpaqtZ8U1JSeN/4ex1XV1cQEa5evcq1lCrwYrep+kJERAT69+8PIoKRkRHMzc3RunVrmJubK0yGhoa15puSkoJJkyZp4QnUA2Osss/ZxcWFazmV8GqtOb6zdu1aZGVlYf369Xjy5Amys7OrTQ8fPqz8rKenV6vhvby8cP78edjZ2XH9iEoTHx+PIUOG4Pbt21qN9Xm/D2B94dq1a/D09ISenh5MTU1hamoKJycnhdcTEYqKiqo1emJiYuUxIkLbtm21+CSq4+TkhEaNGiEmJgbu7u5cywEgmrlOJCQkwN/fX+nrGWMwMjKCkZFRvYqJleH1Xg2+mFlsACpJaWkpbt26hU6dOnEthTeMHDkSv/zyC0pKSriWAkA0s9KkpKTA1tYWBgYGXEvhDfb29nBxceHNK1+imZXk2rVrNcbHusq8efOwfv16yOVyrqWIZlaWhIQEdOnShWsZvMPT0xMGBgY4ceIE11JEMyuLWDNXD2MMc+fOxbp167iWIppZWcSaWTG+vr5IT09HXFwcpzpEMyvBo0ePIJVKYWlpybUUXtKwYUPMmjUL69ev51SHaGYlqKiV68OsNq6YMGECwsLCcPfuXc40iGZWAjFerh0jIyN8+eWXCAgIwMuXLznRIJpZCRISEkQzK8Hs2bNx9epVGBsbw9vbG99++y2ioqLUPstQEaKZleDatWti408JLC0tER0djfv372P27NkoKCjArFmzYGxsjL59++K7777DhQsXNGZucW6GEpSVleHRo0dcy6g3tGjRAj4+PvDx8QEAFBQU4Pz584iMjMSMGTNw69YteHh44MMPP4Snpyfc3NzUsoutOAVUCcLCwjBlyhQkJiaiWbNmXMup9+Tn51eaOzIyEmlpafDw8ICnpyc+/PBDuLq6KjR3TVNARTMryahRo9C6dWteDA4Ijby8vCrmzsjIQHp6Olq1avXWtaKZ1UBubi4cHR3xxx9/wNXVlWs5gkUmk6Fp06YoLCxEkyZN3jpfk5nFBqCSmJiYYO3atZg4cSKkUinXcgTLw4cP0bJly2qNXBuimevA6NGjIZFIsGHDBq6lCJbMzExYW1u/072imesAYww//PADVq9ejfT0dK7lCJK7d+/i/ffff6d7RTPXkbZt22LRokW83lSnPiPWzFrG398f+fn5CA4O5lqK4BBrZi2jr6+PHTt2YMGCBcjJyeFajqAQa2YOcHFxwZgxYzB79myupQiKzMxMsWbmgmXLluHSpUu8eGVIKNy9e1esmbkgJycHHTp0wKpVq7iWIggKCwshlUqrHflTBtHM70BycjJGjx4NV1dXdO3albcrydc3KuLld30JQjRzHYiOjsYnn3wCb29vdOrUCbdv30ZAQADMzMy4liYIVImXAXEKaK0QEU6fPo2AgACkp6djwYIF2Lt3r7gYjAZQJV4GlDRz+V5/TwHIAEiJyI0xthbAIAClANIBfElEBe+shGfI5XIcOXIEAQEBKC4uxqJFizBixAi1zLsVqR5V+pgBvPve2Xi1H6B++efVAFbXlo+6d2jVBKWlpRQUFET29vbUrVs3Onz4MMlkMq5lCZ6wsDCSSCQUHR1d43VQw3bDb5n5jfNDAPxSWz58NnNxcTFt2rSJrK2tydvbm06fPs3bTdmFxv79+8nU1JTOnz9f67XqMPNbG8G/cf4YgFEK7lX7RvDqRiaTkaWlJfn4+NRaM4iol61bt5KlpSUlJCQodb06zPzWRvCvnfsar/bTZrXlw+eaecCAAbRz506uZegMcrmcli1bRh988AGlp6crfV9NZn7XjeDdAYAxNg6AD4CR5QXVW+bMmYPAwEBxJpwWkMvlmDVrFg4dOoSoqCi17RrwrhvBJzHG+gNYAOBjInquFjUc4uXlBT09Pfz5559cSxE0paWlGD16NOLj4xEZGYnWrVurLW9VNoL/HsB7AP5kjMUzxn5QmyoOYIxh9uzZCAwM5FqKoJk0aRKKiooQFhaGFi1aqDXvWvuZieg2gLeW8yGidmpVwgNGjBiBxYsXIykpCY6OjlzLEST6+vpwdXXVyKCTOJz9Go0bN8a0adOwceNGrqUIlkWLFuH7779HUVGR2vMWzfwGfn5+OHjwoLiCkYZo164d+vXrh61bt6o9b9HMbyCRSDB8+HD897//5VqKYPnqq6+wceNGPH+u3n4DcaJRNfj7+8PDwwO3bt2Co6MjHB0d4eDgABsbGzRoIP7/VxUHBwf06NEDO3bswKxZs9SWr7iikQKuX7+OmJgYJCUlISkpCcnJycjLy0PHjh2rGNzR0RGWlpbiQuR15MqVK/j444+Rnp6Oxo0bK32fuDyXmigsLERycjKSk5MrDZ6UlISSkpJKY1f8dHR0hKmpKdeSec2AAQMwePBgTJ48Wel7RDNrmMePH1ca+3WT6+npvWVwBwcHtGzZkmvJvODixYsYNWoUbty4ofTUWtHMHEBEePjwYaXBT5w4UTm6GBkZiT59+nCskB94eXlh3LhxGDNmjFLXiwsncgBjDK1atUJeXh5OnDiBuLg4TJ06FXFxcaKRX2PJkiVYuXIlZDKZynmJvRkaICEhATt37sSvv/4KZ2dn/Otf/8Lvv/8uvmpVDZ6enmjVqhUOHjwIX19flfISa2Y1UVBQgK1bt8LNzQ0+Pj5o0aIFYmJicPr0aXzxxRfcGvnxY+DLLwE19+uqA8YYBg4cqJ6lzhTNDdVE4vN85ndBJpPRmTNnaOTIkWRkZETDhg2jkydPklQq5Vra/8jNfTVtvSLl5nKtqAq7d++m1q1bU1xcnFLXQ9XJ+epKQjKzXC6nCRMmkL29PW3cuJFyeWaSSsaNq2rmceO4VlTJhg0byNramlJSUpS+RzSzBlixYgW5uLjQ06dPuZZSM8XFVc1cXMy1IpLL5fTNN9+QnZ0d3b17t0731mRmMWZ+B4KDg7Fz506EhobC0NCQazk107QpkJsLjBsHFBe/+s4hcrkcM2bMwPHjx3Hu3DmV1sl4E7E3o46cOXMG8+fPR0REhFrfktAoEgkQFMS1CpSVlWHcuHG4f/8+IiIiYGRkpNb8RTPXgaSkJIwYMQIhISHo1KkT13LqFS9evMCwYcPAGMPJkyfFyflc8uDBAwwcOBAbN24UBz3qSGFhIfr3748WLVrg0KFDGuumFM2sBE+fPsXAgQPh5+eHL774gms59YqcnBx4enqiS5cuCA4O1ujyZqKZa6GsrAzDhg2Du7s7Fi1axLWcekVmZiZ69eqFQYMGYdOmTRqfCy6auQaICFOnToWenh62bNkizlmuA6mpqejVqxemTp2KZcuWaeV3JzYAa2DlypW4cuUKzp49C3194f6qZDIZ7ty5A7lcDj09vbdSgwYNqj1ece5No8bFxcHHxwerV69WejacOhDuv5CK7NmzBzt27MClS5f435dcR4gI169fR3h4OMLDw3H27Fk0b94c+vr6kMlkkMvlkMlkVdKbxyq+y+VyMMbeMviePXswePBgrT6XaOZqCA8Px9y5cxEREQFzc3OtlCmTyRAeHg5vb2+1x5ZEhNu3b1eaNzw8HIaGhvDy8oKvry+2bt36zs9JRFWMLZPJoKen9057X6uMoqFBTaT6MJydmJhIJiYmFBERobUy5XI5zZw5kwwNDal37951WkhQEffu3aPg4GAaN24cWVtbk4WFBY0aNYp27dpFd+7cUV00R0Ccm6EcWVlZZG1tTXv27NFquevWrSMHBwd68uQJrV+/noyNjWnLli11WuQ8JyeHDhw4QH5+ftShQwcyNjamoUOH0tatWyk1NVUwa02LZlaCoqIicnZ2pu+++06r5e7du5esrKwoMzOz8lhKSgq5u7uTt7c3ZWRkVHtfQUEBHT16lPz9/alLly5kZGREPj4+FBgYSPHx8YJd7V80cy2UlpZS//79aeLEiVqtwSIiIsjExISuXbv21rmysjIKCAggiURCO3bsoGfPnlFYWBgtWrSI3N3dydDQkP7xj3/QypUr6fLly1RWVqY13VwimrkGKuYl9+/fX6uGSEpKIhMTEzpz5kyN1yUmJpKLiws1adKEevbsSf/3f/9HkZGRVFJSoiWl/KImM+t8b0ZAQADi4uK02peclZWFAQMGYMOGDfDy8qrxWkdHR/z1118oLS0V3yGsBZ028549e7B9+3ZcvHgR7733nlbKLCwsxIABAzBlyhSMHDlSqXv09PREIyuBzpo5IiICc+fORXh4OCwsLJS+Lz09HWVlZWjfvj309PTqVGZpaSk+++wz9OzZEwsXLqyrZJFa0EkzJycn4/PPP8e+ffvg4OCg9H0hISGYOnUqmjdvjpycHDg5OcHFxaUyderUSWGoQkQYP348DA0NsWnTJnGehwbQOTNXzEsODAyEp6enUvcQEdatW4dNmzbh9OnTcHJyQn5+PuLj43HlyhX8+eefWLVqFe7duwdHR8cqBnd0dETjxo3x9ddf49atWwgPD69zjS6iHDq1PNfTp0/Rp08fDB06FF999ZVS90ilUsyYMQMXL15EaGgorKysasw/ISEBV65cqUxpaWmwtbVFaWkpLl68CBMTE3U9jk4irjWHV6YcNGgQrKyssH37dqX+zD979gzDhw+HTCbDgQMH0Lx58zqX++LFCyQmJsLW1lY0shoQ15rDqzj59OnT6Nixo1Lrmj148AC9e/eGhYUFjh079k5GBgADAwO4u7uLRtYCSpmZMZbBGEss3yIttvxYK8bYn4yxW+U/eb1Oq5OTExISEnD8+HG4uroiKipK4bWJiYnw8PCAr68vtm/frtFXfUTUiKLRlNcTqtkIHsAaAIvKPy8CsLq2fPgwAiiXy2nfvn1kaWlJY8aMoYcPH1Y5HxYWRiYmJrR3716OFIrUBDS0CMxgAD+Vf/4JwCcq5KU1GGMYPnw4UlJSYGZmBkdHR2zatAlSqRQ7d+7E6NGjcfDgQXz++edcSxWpI0o1ABljdwDkAyAA24hoO2OsgIhalJ9nAPIrvr9x7yQAkwDA2tra9e7du+pTrwZSUlIwffp0pKenQ19fH6GhobCzs+NalogCVO7NYIxZElEWY8wUwJ8AZgA4+rp5GWP5RFRj3Mx115wiiAgRERHo3Lmz2FDjOTWZWalBEyLKKv+Zwxg7DMAdwCPGmDkRZTPGzAHkqE2xlmGM1TrhR4T/1BozM8aaMcbeq/gM4CMASQCOAhhbftlYAL9rSqSIiDIoUzObAThcPsigD+BXIjrJGIsBcIAxNh7AXQCqreEvIqIitZqZiG4DcKrm+BMA3poQJSLyLujMCKCI8BHNLCIYRDOLCAbRzCKCQTSziGDQ6nxmxlguXnXj1QckAB5zLYLHcPX7eZ+Iqh2m1aqZ6xOMsVhFw6Yi/Pz9iGGGiGAQzSwiGEQzK2Y71wJ4Du9+P2LMLCIYxJpZRDCIZhYRDKKZa4Ax9m/GWFb5W+nxjLEBXGviGsZYf8bYDcZYGmOMVxsjijFzDTDG/g3gGRGt41oLH2CM6QG4CaAvgPsAYgCMIKLrnAorR6yZReqCO4A0IrpNRKUA9uHVW/q8QDRz7UxnjF1jjO3i+0I3WsASwL3Xvt8vP8YLdN7MjLHTjLGkatJgAP8F8AEAZwDZANZzqVWkZnRuSds3IaJ/KHMdY2wHgOMalsN3sgC0ee27VfkxXqDzNXNNlC+hUMEQvHorXZeJAdCeMWbLGGsE4HO8ekufF+h8zVwLaxhjzni1klMGgMmcquEYIpIyxqYDCAOgB2AXESVzLKsSsWtORDCIYYaIYBDNLCIYRDOLCAbRzCKCQTSziGAQzSwiGEQziwiG/wfXEp54eH4aewAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))\n",
    "ax = world[(world.name == \"United Kingdom\")].plot(color='white', edgecolor='black')\n",
    "gdf_4326.plot(ax=ax, color='red', markersize=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "fe017a49-8156-468b-a3c1-a714e419772e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWgAAAD4CAYAAADB9HwiAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABNsElEQVR4nO2dd3hUxfrHP0PoLfQWOoSOtAgoRQGBUEQQVOBSBO6liAUFUeBnryh6VUBAEOFKRxCQ3osigQChgyT0AAEMIYVAyr6/P/Zk3fS2m7NJ5vM882T2PTPvvLObfDM7Z86MEhE0Go1G43rkMTsAjUaj0SSNFmiNRqNxUbRAazQajYuiBVqj0WhcFC3QGo1G46LkNTsAR1OmTBmpXr262WFoNJpcwuHDh++ISFln+M5xAl29enV8fX3NDkOj0eQSlFKXneVbT3FoNBqNi6IFWqPRaFwULdAajUbjomiB1mg0GhdFC7RGo9G4KFqgNRqNxkXRAq3RaDQuihZojSmICH/99Rdbt24lODjY7HA0Gpckxz2oonF97ty5g7e3N1euXCFv3rwEBwczbtw4+vbty6OPPmp2eBqNy6BH0Jos5eLFi3To0IHHHnuMGzducP36dTZs2MCdO3do2bIlc+bMMTtEjcZlUDntRBUvLy/Rj3q7FrGxsQD4+/vTunVrXnvtNd577z2UUrYyIkKfPn24fPkyR48eNStUjSbdKKUOi4iXM3zrEXQ25uLFizbxc1Xeeecd8ubNS968ealXrx5jxozh/fffjyfOAEopChQoQMuWLU2KVKNxPVIVaKVUFaXULqXUaaXUKaXUa4b9faVUoFLKz0jd7epMUkr5K6XOKaW62tm9DZu/UuptO3sNpZSPYV+ulMpv2AsYr/2N69Ud2vtsTs2aNWnSpInZYSRJQEAAixcv5uOPP6ZNmzY8/vjj7Nmzh08//TTZOu3atSM8PDwLo9RoXBwRSTEBFYHmRr4Y8BfQAHgfmJBE+QbAMaAAUAMIANyMFADUBPIbZRoYdVYA/Y38bGCMkX8JmG3k+wPLU4u3RYsWkhuYP3++AALIqlWrzA4nEf369bPF9+DBg1TLBwcHCyBTp07Ngug0GscB+EoqupTRlOoIWkRuiMgRIx8GnAE8UqjyDLBMRB6KyEXAH2hpJH8RuSAiUcAy4Bll/a7bEfjFqL8Q6G3na6GR/wXopBJ+N86lbNiwgS+//JJmzZrRt29fYmJizA4JgKNHj9KlSxd++eUX5syZw4ULFyhQoECq9UqWLBnvp0ajSecctDHF0AzwMUwvK6WOK6XmK6Xi/rI8gKt21a4ZtuTspYEQEYlJYI/ny7h+zyifMK6RSilfpZTv7du309OlbEdERATTpk1j1apVNGrUiL179wLw+uuvmxpXVFQUV69eZc+ePWzbto01a9YwcuRIatSokWYf48ePZ+TIkbRo0YJVq1Y5MVqNJpuQ1qE2UBQ4DDxrvC6PddoiD/AJMN+wzwAG2dX7EehnpHl29sFG2TJYR9Zx9irASSN/Eqhsdy0AKJNSnDl5isPPz0+qVasmgJQvX95mX7VqlQBy9uzZLI/JYrHIiy++aJvOAKRt27YZ9nf69GmZOnWqFCpUSCIiIhwYqUbjHDBzigNAKZUPWAUsFpHVhrAHiUisiFiAuVinMAACDZGNo7JhS87+N1BCKZU3gT2eL+O6u1E+1xAdHY2IcOvWLZo2bUrr1q2JiYnh5s2btjLPPvss33zzDb169eLhw4dZGt+lS5dYsGABV69eJTw8nKioKPbt25dhf/Xr12fixIm4u7tz4sQJB0aq0WQ/0rKKQ2EdBZ8Rka/t7BXtivXBOtoFWAf0N1Zg1AA8gYPAIcDTWLGRH+tNv3XGf6BdWEfYAEOBtXa+hhr5fsBOo3yO5/vvv6dVq1YULVqUevXqsWzZMgCWL19O3ryJHwB99dVXuXfvHk8++WSWivSWLVsAmDBhAkWKFCFfvnyZ9iki/P3336xbty7TvjSabE1qQ2ygLdavrscBPyN1B34GThj2dUBFuzpTsE5HnAO62dm7Y10FEgBMsbPXxCri/sBKoIBhL2i89jeu10wt3pwyxeHt7S3PPvus3Lt3T1599VUB5Ny5c1K+fPlkpzPu3r0rffv2ldatW0tUVFSWxBkbGyuTJ08WQK5cueIQn8ePHxdA/Pz8HOJPo3EmOHGKwylOzUzZWaBjY2Olf//+Urx4cQGkf//+tmsdO3YUQJo0aSKA/Pe//41X986dO7J06VKpW7euAPLee+9lWdx3796VggULSpUqVeT06dOZ8hUbGytt27YVT09PiY2NdVCEGo3zcKZA682SXIipU6eybNkyevbsiYgwZswY27WtW7fSuHFjihQpQkhICO7u7gBMmjSJL774AovFAsALL7zAU089xVNPPZVlcZcoUYK7d+/SqlUrevfuzcmTJzM01REZGUnfvn25fv06mzdvJk8e/aCrJpfjLOU3K2XnEXTTpk2lU6dO6aoDSLFixWT+/PkOiyM2NlYuXLggMTExNpvFYpEXXnhBqlSpIidOnEiyXnBwsFStWlWmT5+eoXZffvlladGihdy4cSND9TUaM8DsVRwa5xMQEICfnx/dunVLV72ZM2cSFhbG8OHDHRKHiJAvXz5q1qxJ3rx5UUqhlCJPnjwcOXKEq1evMm/evCTrFitWjLp16/Lrr7+mu10/Pz+WLVvG4sWLqVChQma7odHkDJyl/Gal7DqCfumll6R27doSHR2d7rqAVK9ePVPtP3z4UObNmyejR4+2rWeeO3eurF69Wry9vWXWrFly9epVadasmQBy//79ePVXrlwpgCilZO7cuelq+9atW1KoUCH57rvvMtUHjcYM0DcJc75AHzlyRIoVK5buerGxsQKIj49Phtv29fWVihUrCiBdu3aV9evXy4YNG8RisSQqe+XKFQHkzJkz8ezTpk2T8uXLJ7Knhbp164qnp2eS7Wk0ro4zBVpPcbgIgYGBNG3aNN31jh8/DpBo+860MG3aNJRSeHl50axZM6Kioti8eTM9evSge/fuSfqsUqUKzz//fLxpjOjoaBYvXkxQUBCFChVKdxyFCxemUKFCBAYGpl5Yo8lFaIF2EbZu3Yq3t3e66zVt2pQ333yTli1bpmurzkOHDvHmm2+ycOFCzp8/z4YNG9K88uLRRx9l+/bttpUj27dv5/Tp0+zcuZNq1aqluw+HDh2iZ8+eNGvWLN11NZocjbOG5mal7DrF8cgjj8iBAwcyXL9IkSKydu3aNJdv3Lix9OzZM0NtRUREiIeHh4wYMcK2Pnv8+PEZ8hXHtWvXBJC7d+9myo9Gk9WgpzhyNlFRURw/fpxGjRpl2EeNGjV45pln0lRWRDhx4gQjRozIUFuFCxdm165d/PjjjwDMnj2bd955J0O+4li71vp0f0RERKb8aDQ5CS3QLoCPjw+1a9emSJEiGfaxZ88eADZv3pym8nXq1KFPnz589913GWrP09OTU6dOsWHDBkaNGmV7cCajxG2w5OGR0lbjGk3uQgu0CzB//nxGjhyZKR9xm+Lv3r071bJKKaZPnw5A//79M9xmgwYNKFiwYIbrx3Hjxg1WrFhh23hJo9FY0QJtImFhYdy4cYNNmzbx9NNPZ8pX4cKFAVJdCRIVFcWvv/5K7969+fzzzylXrlym2nUER48exWKx0KVLF7ND0WhcCi3QJvHbb79RuXJlKlWqRFBQULpOHkkKPz8/ADp37oyIEBwcnKjMsWPHKFCgAAMHDuT999/nlVdeyVSbjqJhw4aA9WlKjUbzD1qgnUS3bt0oXLgwN27cSPL6vn37CA0N5cqVK/j4+KTp3L6U+Ptv6zkGw4YNo169epQuXZru3bsze/Zsdu3axbhx42jatClfffUVERERTJw40TbqNptq1arRunVrNmzYYHYoGo1LoayrRHIOXl5e4uvra3YYNGrUiAcPHnDx4kUWL15M//79eeedd/D19aV58+YsWrSIQoUKcebMmQw9ZJIQEaFnz554e3sTHR3NwIEDWbFiBYcPH8bHxwcvLy+eeeYZ+vbt65K7xLVq1YqAgADu3LljdigaTbpQSh0WES9n+NbbjTqJcuXKMWjQIO7du8eAAQNo0KABV69eZfPmzdSpU4dVq1bRokULh4gzWG/8JRyBvvrqqw7x7WyCg4M5ePCg7dQYjUZjRQu0kxg7diz9+vVj4MCBdO/enSZNmgDWx6vHjx9vcnSuQ2hoKKVLl6Zs2bL07dvX7HA0GpfC9b7r5hD69u3Lhx9+yL179/jjjz8AmDNnjhbnBKxYsQKAI0eOJHnWokaTm0l1DlopVQX4H1Ae6zaUP4jIt0qpL4GngSisZwwOE5EQpVR14AzW8wgBDojIaMNXC2ABUAjYCLwmIqKUKgUsB6oDl4DnReSucWDtt1jPMrwPvCgiR1KK11XmoO358ssv2bhxI7t27TI7FJcjT548NG/eHB8fH9zc3MwOR6NJN86cg07LCDoGGC8iDYDWwFilVANgG9BIRB7BehDsJLs6ASLS1Eij7eyzgP9gPenbE4jbHehtYIeIeAI7jNcA3ezKjjTqZzsiIiLYvXs3oaGhZofiUhw6dAgRYcWKFVqcNZokSFWgReRG3KhVRMKwjo49RGSriMQYxQ4AlVPyo5SqCBQXkQPGBiP/A3obl58BFhr5hQns/zP2JDkAlDD8ZBveffddPvjgAy5dukTx4sXNDselmDlzJi+88AI1a9Y0OxSNxiVJ1xy0MX3RDPBJcGk4sMnudQ2l1FGl1B6lVDvD5gFcsytzzbABlBeRuAXDN7FOp8TVuZpMHfu4RiqlfJVSvrdv305Pl5xKSEgIH330EfPnz8/QNpw5maioKDZu3MiECRPMDkWjcVnSLNBKqaLAKmCciITa2adgnQZZbJhuAFVFpBnwBrBEKZXmoaMxuk7X4mwR+UFEvETEq2zZsump6lRiY2MBHHZeYE7iwYMH3L17Fz8/P6ZPn257rzQazT+kSaCVUvmwivNiEVltZ38R6An8yxBWROShiPxt5A9jvYFYBwgk/jRIZcMGEBQ3dWH8vGXYA4EqydRxeYoWLcqgQYPMDsMlKV68ODNnzmTJkiW8+uqrPP7448TExKReUaPJRaQq0MZKih+BMyLytZ3dG5gI9BKR+3b2skopNyNfE+sNvgvGFEaoUqq14XMIsNaotg4YauSHJrAPUVZaA/fspkJcmr1791KwYEEWLVrEU089ZXY4LsnIkSPZuXMnISEhHDx4kC5durB161Y++eQToqKizA5PozGdtIyg2wCDgY5KKT8jdQdmAMWAbYZttlG+PXBcKeUH/AKMFpG4nXteAuYB/lhH1nHz1p8DnZVS54GnjNdgXYp3wSg/16ifLViyZAmPPfYYq1evZtu2bWaH49K4u7sTFBRE48aNGT9+PB9++CFvv/22XvWiyfXovTicQGxsLHnz5uXLL7/UN8EyQGBgIGPHjmXHjh106tSJRYsWUbRoUbPD0miSxOx10Jp0EBMTQ79+/XB3d8/wkVK5gaioKG7dupXkzUEPDw/WrFmDv78/V65cYeLEiSZEqNGYjxZoB3Pu3DnWrFmDr68vJUuWNDscl+Szzz6jQIEClC9fnscffzzZcuXLl+err75i1qxZXLt2LdlyGk1ORQu0g2nQoAGtWrVi7969ZoficsTGxjJ37lwmT57M999/z2+//cbBgwdTvCHYvn17ANs2qhpNbkILtINRSlGjRg3bBvqaf9iyZQsjR45k9erVjBkzhsBA64rJlEbHbm5uhISEUKpUKSZOnKjXS2tyFVqgHcz9+/fZsmWL3jozCdq0aUPt2rU5f/48AKNGjWLChAmpniz+888/s2/fPr755ht+/fXXrAhVo3EJ9P6ODiY4OJi7d+/q/SWSwN3dnZ9++onnnnuOypUrU6NGDUaOHImXlxeDBw+mXr16zJ07lxYtWtC2bVuUUpw4cYKPP/4YgH//+9/06tXL5F5oNFmHFmgH4+Fh3SrEx8eHVq1amRyN69G2bVv69evHqFGjsFgsnDp1ijfeeAMvL+sqpSZNmvD666/j7u6Om5sbjRo1ws3NjaioKPLly2dy9BpN1qKnOByMUopatWrZNqLXJGb69OmEhYUxefJkhg4dyv/93/8hIty4cYOjR48SHh7OokWLGDBgAMePH2ffvn1anDW5Ei3QTqBKlSq0bdvW7DBcnrfeeovLly/b5qArVKiAUooiRYrQs2dPZsyYoaeLNLkaLdAOJDY2lo8//pgLFy7QrVs3s8NxeZRSXL58mT179pgdikbjkug5aAdx+fJlhg4dSnR0NDt37qRgwYJmh+TyuLm58fvvv/P8888TGxurT1XRaBKgBdpBLFq0iD179nD9+nUqVsxWh76YSuvWrbl16xZhYWGUKFHC7HA0GpdCT3E4iMmTJ/Paa68xbtw4s0PJVri5uTF8+HCaNWvGZ599ZnY4Go1LoXezcyCBgYFUrlwZi8WCdctrTVoQEXbv3k3Hjh3p2bMnFy9e5NtvvyU8PJz8+fPr+XyNS6N3s8smVKpUiQoVKrB9+3azQ8lWKKXo0KEDp06dokmTJuTJk4ennnqKr776iu7duxMUFGR2iBqNKWiBdiBKKUJDQ1N9dFmTNA0aNODjjz/m2LFj3Lp1i7179/Lyyy/blt+FhYWZHaJGk6VogXYwb7/9NuvXrzc7jGyNUoq4w38/+eQTGjduDMDgwYPNDEujyXK0QDuYLl26mB1CjqJ48eIcPXqUQoUK2Xa/02hyC2k5NLaKUmqXUuq0UuqUUuo1w15KKbVNKXXe+FnSsCul1HdKKX+l1HGlVHM7X0ON8ueVUkPt7C2UUieMOt8Zh8om24YrM3LkSBo2bGh2GDkKHx8fypQpo/fY1uQ60jKCjgHGi0gDoDUwVinVAHgb2CEinsAO4zVAN6wneXsCI4FZYBVb4D2gFdASeM9OcGcB/7Gr523Yk2vDZcmTJw+nTp3S+xY7kNq1awPoZXiaXEeqAi0iN0TkiJEPA84AHsAzwEKj2EKgt5F/BvifWDkAlFBKVQS6AttEJFhE7gLbAG/jWnEROSDWNX//S+ArqTZcll27duHu7s4bb7xBTlvCaBblypVjx44dfPTRRwQEBJgdjkaTZaRrDlopVR1oBvgA5UXkhnHpJlDeyHsAV+2qXTNsKdmvJWEnhTZclhIlSnDhwgW2bt1K06ZNuX37ttkh5QgqVqxIwYIFCQkJMTsUjSbLSLNAK6WKAquAcSISan/NGPk6dbiYUhtKqZFKKV+llK/Zgvjnn39Srlw56tWrR0BAAKNHjzY1npxCaGgoDx48wMvLC6UUI0aM0N9QNDmeNAm0UiofVnFeLCKrDXOQMT2B8fOWYQ8EqthVr2zYUrJXTsKeUhvxEJEfRMRLRLzilmeZxb///W+KFSvGmjVriIiI4NixY6bGk1OoVKkSp06d4rfffmPAgAHMnz9fn/StyfGkZRWHAn4EzojI13aX1gFxKzGGAmvt7EOM1RytgXvGNMUWoItSqqRxc7ALsMW4FqqUam20NSSBr6TacEkiIiI4ffo0+/bt4/Dhw6xbt44TJ06YHVaOoU6dOvzwww8sXbqU0aNHU6VKldQraTTZGRFJMQFtsU4tHAf8jNQdKI11ZcV5YDtQyiivgJlAAHAC8LLzNRzwN9IwO7sXcNKoM4N/9ghJso2UUosWLcQsoqKiBJBx48aZFkNOJiAgQABp3Lix2aFoNDYAX0lFlzKa9GZJDmbOnDmMHj1az486gcjISAoXLsz27dvp1KmT2eFoNIDeLClbEXdo7MKFC1MpqUkvR44cAdDHiWlyDVqgHUzcWXovvvgiL730ktnh5CjintCMiIgwORKNJmvQAu0Exo4dyzvvvMPKlSvNDiVHUaJECZ544gn++OMPs0PRaLIELdBO4OHDh3z00Ud8/vnnZoeS42jfvj1r1qwxOwyNJkvQAu0E9Pyz8xgzZgzLli0jLCwMi8VidjgajVPRAu1g/Pz8GDVqFE8//TSDBg0yO5wcR8WKFXnmmWcoXry4XsmhyfFogXYwcVtirlu3jgIFCpgcTc5kyZIl/PXXX+zevVufsqLJ0WiBdjB9+/YF0DeynIynpye9evXivffeMzsUjcZpaIF2MPv37wegadOmREdH6wdWnMjEiRP573//y8OHD80ORaNxClqgHUznzp0BKFq0KIUKFaJ3795apJ1EmzZtAPRJK5ocixZoB1OiRAlKlixJkyZNWLduHTt27ODChQtmh5Vjad++PQsWLDA7DI3GKeQ1O4CcSHBwsC0/fPhwHnvsMYKCgjCOWtQ4iH379rF3715Wr16demGNJhuiR9BO5uOPP+b27dv6sW8nMG/ePAYNGkSfPn3MDkWjcQpaoJ1M8eLF+eyzz5g9ezanTp0yO5wcRYUKFQgPDzc7DI3GaWiBzgLefvttBg0axOuvv252KDmKZs2aERoamnpBjSabogU6i/j+++/Zv38/Xl5evPHGG2aHkyO4e/cu0dHRZoeh0TgNfZMwiyhWrBg9evRgxYoVlCtXzuxwcgTdu3fnpZdeYvfu3Tz55JNmh+OynDp1ioMHD1KhQgWUUpQoUQIfHx8AmjRpgpubG6VLl6ZBgwYmR6pJiBboLGT58uV4eXnh7+9vdig5ggcPHgBQq1YtkyNxXVasWMHYsWPp2rUrt27d4t69e0RFRdGsWTN++uknKlSowM2bNwHrnH6xYsWoWbMmW7ZsoUiRIgwfPpxixYrRt29fmjZtSp48+kt3luKss7TMSmaeSZgWFi1aJNa3XZNZWrRoIQ0aNDA7DIeD9QxQCQwMzJQfi8UidevWlU2bNiV5/cKFCyIiMm3aNAkLC5OffvpJatSoIQsXLpQePXrIs88+a4ulTJkyMmrUKLFYLJmKKSeCE88kTMuhsfOBW8BJO9ty/jlA9hLgZ9irA5F212bb1WmB9RBZf+A7/jkYthSwDevBsNuAkoZdGeX8sR5Y2zwtHXJ1gb5x44YAcvjwYbNDydbcvn1bAPHz8zM7FIcSGRlpE0RAfvnll3jXN2/eLB9//HGafP33v/+1CWxGiasPSMWKFWXUqFEZ9pVTMVug2wPN7QU6wfWvgHflH4FOrtxBoLUhvJuAbob9C+BtI/82MNXIdzfKKaOeT1o65OoCLSIyb948KVWqlMTGxpodSrZl06ZN0rlzZ7PDcDgnTpwQQLZu3SrDhw8XQIYMGSIxMTFy9uxZm1jOmzcvVV9vv/12kiKfFiwWixw9elQCAwNtbcbFdufOHTlw4EBGupcjMVWgJQXhNcTzKuCZSrmKwFm71wOAOUb+HFDRrtw5Iz8HGGBXx1YupZQdBFpEpHTp0vLee++ZHUa25YUXXhBApk6dKsePHzc7HIcQEhJiE8NPP/1URESuXLki9evXt9mLFy8uEydOFECWLFki4eHhEhMTE8/PpUuX4o1800N4eHi8uvbf9pK6Vr9+ffnss8/k1q1bjnkTsiGuLNDt7YMzykUAR4E9QDvD7gVstyvXDlhv5EPs7CruNbAeaGt3bQfglUx8IwFfwLdq1aoOfvudw86dO6V8+fJy7Ngxs0PJljx48ED27dsnTzzxhABSq1YtuXz5stlhZZgjR45IhQoVBJAFCxbI2bNnbddiY2Nl0aJFcujQIRERuXfvXiKhjEtVq1aV0qVLS/ny5SUjfwtJ+WzVqpW88sor0qBBA+natas89thj0qdPH2nbtq00btxYOnfunOmplOyMMwU6s6s4BgBL7V7fAKqKyN9KqRbAGqVUw7Q6ExFRSqV76zcR+QH4AcDLyytbbB3XoUMHXn31VZ5//nlee+01xowZY3ZI2YoCBQrQtm1bdu/eTXBwMKVLl2bbtm2MGDHC7NDSzYMHD2jevDldu3bl/PnzFC1aNN71PHny8K9//cv2unjx4nEDEwAiIyPZv38/QUFBtGjRgpo1a5IvX750xxF3anocJ06c4Pnnn8fd3R1/f3+GDRvG7t27OXToEEePHqVRo0a2snH7zCilaN++PXXq1KFRo0Y8fPiQRx99lCeeeEKvAMkAyv6DTraQUtWxjngb2dnyAoFACxG5lky93cAEo9wuEaln2AcAT4rIKKXUOSN/QylVEdgtInWVUnOM/FKjjq1cSrF6eXmJr69vqn1yBSwWC5s2baJnz57s2bOH9u3bmx1StmXgwIEsXbqUS5cuUa1aNbPDSTMxMTHUqlWLK1eukJa/RWcgIkRFRVGwYEGbTSmVajwJr8fExHDy5EmuX7+Ov78/586dw2KxsGPHDtzd3Tlw4ABubm5O6YOZKKUOi4iXU5ynZZhNElMcgDewJ4GtLOBm5GtiFeZSxuuENwm7G/YviX+T8Asj34P4NwkPpiXW7DIHbc+gQYNk7ty5ZoeRrfHx8Yn3tbxPnz5mh5Qm/Pz8BJDFixdnWZtnz56Vp556SqZMmZJoOuPxxx+Xhg0bStOmTWXmzJnJTqWQjumM6OhoAcTb29tlluk99dRTAki5cuUEkBIlSsj69esz5AuTV3EsxTp1EQ1cA0YY9gXA6ARl+wKnsC6xOwI8bXfNCzgJBAAz+Gf0Xhrr/PJ5YLudoCtgplH+BMnMPydM2U2g79+/L4ULF44356jJGBEREfLmm29KkyZNBJCWLVvKsmXLzA4rReJEMCtJSmxLlSplywcHB4uIyNWrV5MV5zp16qSrzYsXL9rqRkdHO6Nbibh//76cO3dOgoOD5bXXXkuyHx9//LEMGTJEhg4dKoDcuHEj3e2YKtDZLWUXgf7zzz/l66+/loIFCwqQY1YiuAonTpywPWhx/PhxmTdvnjx48MDssBIxYsSILBOtpUuX2oQpb968iW4E2gt0hw4dBIi7JxRvFUlc/pFHHklX3Lt27RJAOnXq5JD+hIWFyffffy9vvfWWdOvWTQYMGCA//vijzJ0713bDNWHq06ePrFy5Uu7evZvI34ABA2TChAnpjkMLdA4T6LiHLAB56aWXZOrUqRIeHm52WDmOzZs3CyDFihUTLy8v6dixY6afznM0kZGRopSSrVu3OrWduAdf4tLGjRvlhx9+SHEKIy4VKFAg3uu9e/dmaKpDxDowAcTf3z/NdXbu3CkLFy6Uhg0byrfffisTJkyQdu3aSf78+aV06dICSMeOHaVmzZrxVpQ8//zzIiISHBycpmcOAClSpEi6+mPU0wKd1pQdBPrBgwfi7u4ugPj4+JgdTo4mIiJCHj58KFFRUfLKK69I27Zt5dixY4nWDpvFp59+KmBdz+0s7AcEcenmzZu263fu3JGGDRumSazjBDkmJibeqPT555+3/S67u7uLh4dHktMF9mu6ExIbGytnzpyRo0ePirEayzaSd3NzE0DGjBkjH330kaxcuTJeHxwBIJMmTcpIPS3QaU3ZQaBFRM6fPy+AXgedhcTGxkqPHj0EkH79+pkdjoiIzJ4926lz0GvWrIknrlFRUUl+vRf5ZxsCQI4ePSoiYnttsVgS3eBLSrw9PT0T2eLWp8c9YANIvXr1ErUfd61BgwY2UV66dKn89ddfDn1PksJisQggu3fvTnddLdA5UKBjY2Ntc3whISFmh5OreP755+OJ4tq1a6VNmzby5ptvZnks69atc5pAx23MBciWLVvSXO/hw4eplokTNEAKFy4stWrVkkKFCgkgefLkSXUEnpC4a99++22a43QkcZ9DVFRUuus6U6D1dqMmkSdPHq5fv07lypXx9PTkwIED1KxZ0+ywcgX58uWjePHiAPz888+8/vrrlC1bloULF/LBBx9QqFChLItlx44dTvM9aNAgAGJjY9P1kEj+/PlTvC4iNn/79u2jXLlyRERE0Lx5c8C6vj85GjVqxNixYylfvjzXrl3jxo0bnD9/HoBp06bx6quvpjlOR3Hp0iV69epF7969M/SAjzPRj/aYSMWKFbl16xb58uWjVq1a+Pj4WL/WaJzK4cOH2bdvHxs2bGDIkCFs2rSJVatWcevWLQoXLsytW7eyLJajR4/y2muvOcX3+PHjAdi7d69D/Xbo0MGWb9u2LXXq1KFZs2aA9Rgy+xFgSEhIvLonT57E09OTBw8e0LBhQwYPHsySJUuIjo62xZtRlFIopejevXuK5f7++286d+5MgwYNUEpRo0YNAFauXJmp9p2Cs4bmZqXsMsWRkLiHBoYPH252KDmaa9euSZEiRSQkJERq164tr732WrzrzzzzjLz++utOjyMiIkLeeustgbTtTJcRRo0aJYDDVwjdvn1b1q5dm6ayBw4ciDe94eh7LsHBwane1Hz00UfFzc1NTp8+LSL/TKcsXrxYKleunKHVKPag56BzvkCLiHTq1Elmz55tdhg5mqFDh8qYMWNk27ZtUqlSJbl9+3a865988okULVpUfH19ZcOGDWn2Gx0dLbdu3ZLQ0NA0lX/00UdtwuCsbWePHz8uOHDdcUb4+uuv44ll3759HeLXYrFIz5494/kuVaqULFu2TJYvXy7NmjVLVbjPnTvnkFi0QOcSgZ4yZYpUqVIlXWtENemjRYsW8sknn0i3bt1kwYIFia4nXC/cv3//VB/GuHXrlnTt2jVenbCwsGTLx8TESMuWLaV3796Z7k9KbNiwIdOjQ0cQF8OiRYsc4u/hw4fxPqO4kXFqfPLJJ/HqOeqGpDMFWs9BuwBXrlxhw4YNfPDBB0RGRtpummgcT5MmTfjtt9/YvXs3bdq0SXQ9bsOg7777jpCQEFauXEm+fPm4fft2sj4XLFjAli1buHXrFr169WLZsmXJHmIbHR1N586dOXjwIIMHD3ZIn5IjtbnYrOLEiRNYLJZ4O/KlFxHBz8+Pxx9/nAIFCgDwwgsvICLUr18/TT7sN9EaOXJk9ti4yVnKb1bKjiPouHkwd3d3adSokctsKJMT+eCDD6RSpUpy8ODBJK/fv39fduzYYZt2iDtR5IUXXkhUNu7EEkCGDRtms9+8eVMAuX79eqI6rVq1kp49e2bZY+dx8WX303uOHDkSb/T7448/pqt+wrnq77//Xl588UWHxIae4si5Ah0UFCRubm5y9OhR6dSpk8tv7pMbadiwobRu3TqeLTY2VpRS0qJFC/n9998T1QFkw4YNMm3aNAkKChIR6xOkQJrWGTuK1atXxxOmPn36ZNsBQGbet4SnwRw9elTKlSvnkM9CC3QOFmhfX1/Jnz9/th/h5GQGDhwogMydO9f2iPiUKVMkb968cu/evSTrDBkyRGrVqhVvDjjuQNi4p/SyCh8fH/nss89sseTWA4sPHz4c7xtFhw4dHLLNqxboHCzQt2/fFg8PD9m5c6fZoWiS4fLly7a9KqZMmSLvvPOOAImW6NmzYMECAeu5gWDdU6JatWoCJHlzUpP1fPTRR1KtWrVM+3GmQKfpRJXsRHY6USWONm3a0KRJE77//nuzQ9GkwNmzZ203pHx9fWnRokWyZUNDQ1myZAkjRozg999/Z+HChVSoUIF27drRo0ePrApZkwKRkZF4eHhw4sQJPDw8MuzHmSeqaIF2AVq2bElgYCCBgYFmh6JJhQMHDlCmTBlq165tdigaBzB69GiqVq3K5MmTM+zDmQKtl9mZzOzZszl06BC9e/c2OxRNGmjdurUW5xxEtWrVmDJlitlhJIseQZtM/fr1qVq1Kps2bdKnHms0WUxoaCju7u7p3lDKHlNH0Eqp+UqpW0qpk3a295VSgUopPyN1t7s2SSnlr5Q6p5Tqamf3Nmz+Sqm37ew1lFI+hn25Uiq/YS9gvPY3rld3WK9dCG9vbx555BEtzhqNCRQvXpx8+fIRHh5udihJkhZVWID1BO+E/FdEmhppI4BSqgHQH2ho1PleKeWmlHLDegBsN6ABMMAoCzDV8FUbuAuMMOwjgLuG/b9GuRzBvXv3ePDgAQAXLlxg2rRphIWFmRyVRpM7KVCgAK46k5CqQIvIXiA4jf6eAZaJyEMRuQj4Ay2N5C8iF0QkClgGPKOUUkBH4Bej/kKgt52vhUb+F6CTUT7b0r9/f5RSlChRgkKFCjFp0iT69esHQJEiRUyOTqPJnbi7uyfaFtVVyMz36peVUseNKZCShs0DuGpX5pphS85eGggRkZgE9ni+jOv3jPKJUEqNVEr5KqV8U9ozwUx2797N8uXLmTZtGnfu3OHy5csEBQUxZMgQpk2bpqc4NBqTKF26NMHBaR2DZi0ZVYVZQC2gKXAD+MpRAWUEEflBRLxExKts2bJmhpIkFouFLl26MHbsWMaPH0/p0qWpWrUq8+fP5++//870RuUajSbjhIeH8/fff5sdRpJkSKBFJEhEYkXEAszFOoUBEAhUsSta2bAlZ/8bKKGUypvAHs+Xcd3dKJ/t2LZtG7Vr1+bbb79NdK1UqVJpd3TnDgwbBvfvOzA6jSZ34+7unrNG0EqpinYv+wBxKzzWAf2NFRg1AE/gIHAI8DRWbOTHeiNxnfGY5C6gn1F/KLDWztdQI98P2CmuOpOfAH9/fxYtWkT79u0pVqwYAwYM4Lnnnsvc9oZ37kDZsrBgARQpYn2t0WgyjauKM5D6obFKqaXAk0AZpdQ14D3gSaVUU6wbj1wCRgGIyCml1ArgNBADjBWRWMPPy8AWwA2YLyKnjCbeApYppT4GjgI/GvYfgZ+VUv5Yb1L2z2xns4KgoCA8PT3x8PCgV69evPLKK5w7d46XX345c47ffDPx659+ypxPjSaXEx4eTkhICE888YTZoSSJflDFwaxevZpx48YREBDg2BOC79+3jpzjiIiAwoUd51+jyYXUqVOHsmXL8scff2TYh37UOxvRpUsXGjRowKBBg7h3757jHBcuDLdvw4svanHWaBzA7t27OX/+PB9++KHZoSSLFmgHU7RoUVavXs2ff/7Jl19+6VjnZcpYpzW0OGs0GSYqKoq5c+cyYsQIPv/8czp16mR2SMmS6hy0Jv0ULlyYJ598koCAALND0Wg0dvz0008MHz6chg0b8sUXX9C3b1+zQ0oRLdBOok6dOpw8eTL1ghqNxuls3bqVrl1tWwNlm79NPcXhJMaOHcuWLVs4dOiQ2aFoNLmay5cvxxPn9evXmxhN+tAC7SRKlizJa6+9xjfffGN2KBpNriQ2NpZJkybRsGFDypcvD0CfPn2y1Yk2eorDiYSHh1OtWjWzw9BocgWXL19m8eLF+Pr6cujQIa5du0arVq3Yvn07rVu3Nju8DKEF2olcvHiRXr16mR2GRpPjidvocvDgwfTr149BgwbRsmVLPDw8yM6bYGqBdiKrV6+mRo0anDlzxnbYqEajcSzXr1+35f/3v/+ZGInj0XPQTmTGjBmEh4fTrl07li5danY4Gk2OQUQIDQ1lypQpthO5XXVP58ygR9BOZOzYsQBUqlSJGTNmMGDAAJMj0miyJ1euXGHz5s1s2LCBK1eucPbsWdupRHnz5uWjjz7C3d3d5CgdjxboLKBSpUoU1k//aTTpwmKx4ObmRokSJQgJCeGZZ57B29ubpk2b0rhxY4oWLWp2iE5HC3QWsG7dOlq2bJl6QY1GA8DSpUsZOHAgAI0bN2br1q0ULFjQ5KiyHj0HnQX89ttvtl82jUaTMitWrGDgwIHUrVsXi8XC3r17c6U4gxboLKF79+68+uqrZoeh0bg8a9eu5YUXXqBu3bqcPXs2Wy+RcwRaoLOA4cOHkz9/frPD0GhcmmPHjtG7d29atWrF2bNnzQ7HJdACnQXcvn2bMmXKmB2GRuOyBAYG0rRpU6ZOncqBAwfMDsdl0AKdBVy+fJl69eqZHYZG47L06NGD9u3bM3HiRLNDcSlSFWil1Hyl1C2l1Ek725dKqbNKqeNKqV+VUiUMe3WlVKRSys9Is+3qtFBKnVBK+SulvlPG5JJSqpRSaptS6rzxs6RhV0Y5f6Od5g7vfRZx+fJl22YtGo0mPgcPHuTq1ats3LjR7FBcjrSMoBcA3gls24BGIvII8Bcwye5agIg0NdJoO/ss4D9YT/r2tPP5NrBDRDyBHcZrgG52ZUca9bMl169fp0CBAuS08x81GkewefNm+vTpQxH7Mzc1QBoEWkT2Yj1V2962VURijJcHgMop+VBKVQSKi8gBsarU/4DexuVngIVGfmEC+//EygGghOEn29GhQweGDBnCCy+8wO3bt80OR6NxKfbu3cuzzz5rdhguiSPmoIcDm+xe11BKHVVK7VFKtTNsHsA1uzLXDBtAeRG5YeRvAuXt6lxNpk624r333iM4OBhfX1/atm1rdjgajUtx6tQpGjdubHYYLkmmBFopNQWIARYbphtAVRFpBrwBLFFKFU+rP2N0ne55AKXUSKWUr1LK1xkj1NjY2ExvxFKyZEl69OhBZGQkFovFMYFpNNmcwMBAbt68SeXKKX4Jz7VkWKCVUi8CPYF/GcKKiDwUkb+N/GEgAKgDBBJ/GqSyYQMIipu6MH7eMuyBQJVk6sRDRH4QES8R8SpbtmxGuxSPjRs30qtXL3bt2kXVqlUpU6YM9+/fz7C/GTNmMGPGDK5evZojd93S5E7++OMPoqKiMlz/5MmTNGrUKNc/kJIcGRJopZQ3MBHoJSL37exllVJuRr4m1ht8F4wpjFClVGtj9cYQYK1RbR0w1MgPTWAfYqzmaA3cs5sKcThHjhzh2WefpWrVqnh4eNCjRw9+++03hgwZwvfff0+3bt2YN29ehm707dy5k1deeYUqVaqwcuVKSpUq5YQeaDRZy/3792nbti0LFizIsI/IyEhq1KjhuKByGiKSYgKWYp26iMY6DzwC8Mc6P+xnpNlG2b7AKcN2BHjazo8XcBLrqHoGoAx7aayrN84D24FShl0BM43yJwCv1GIVEVq0aCHpYceOHTJ69GipUKGCfPbZZ3L69Gnx9/eX48ePCyCjRo0SEZFTp07FTb+Ip6enfPfdd3L37t00tXH16lUBZNeuXemKTZM8gHTs2NHsMHI1tWvXFkBCQ0Mz7OPkyZNSq1YtB0aV9QC+kgZtykhyilMzU3oF+oMPPhBAvv3223j2kJAQKVq0aDwR/uGHH2TdunWyc+dOefzxxwWQdu3ayZo1a2Tv3r1y/vz5ZNt5+eWXBZD169enKz5N0gDy7LPPmh1GrsXHx8c2YDl+/HiG/cTGxkq5cuXkwoULDowua3GmQOf6Jwmfe+45qlSpwqFDh2wbgAO4u7sTFhZGiRIlbLb//Oc/PP3003To0IE//viD0NBQatasybvvvkvv3r3x9PRk9+7dSbbz/vvvA7B//34n9iZ3EHeT1VHnPVosFn1fIJ3ky5fPln/kkUcyPIecJ08eevTowapVqxwVWs7CWcpvVkrvCFpE5O+//5ZWrVpJ27ZtJSwsLN3149i0aZNUqlRJbt++Hc8eGRkpBw4ckNKlS8v9+/cz7F9j5V//+pdt9AZIwYIFZfv27enycebMGfH395fDhw/b/GjSz507d2zvX9myZW358ePHy6JFi+TmzZup+vD19ZXKlStn278N9AjauZQqVYpdu3YBMHPmzAz78fb2ZtiwYXTv3p2goCAA7t69S+XKlWndujVvvvkmhQoVckjMuZl+/frFe/3gwQPCwsLS5aN+/frUrl2bFi1aADB//nyHxZeb2L9/P02aNAHg//7v/9i3bx+DBw/m4sWLDBo0iHr16lGlShX69OnD+++/z6+//sqpU6d4+PChzUeLFi2oVq0ae/bsMasbLkvcjbocg5eXl/j6+mao7qFDh2jfvj0XL16kQoUKGfIhIrz77rvMmjWLAQMGcPjwYerUqUPPnj3p2bNnrt143JGEh4dTrFgx2+uM7BYYExNj+5oeHBxMyZIlHRpjTicmJoZ33nmHefPmMW/ePJ5++mny5Ek83hMRAgICOHLkCH5+fpw4cYK//vqLy5cvU7lyZSpVqkT+/Pnx8/Njz549NGzY0ITeZA6l1GER8XKKc2cNzc1KGZnisGfs2LHSvXt3sVgsmfLj7+8v77zzjsyfP19iYmLSVffTTz9Ntn2LxSLly5eXo0ePZiq+7EpgYGC86Y1///vfGfJjsVhsX8U16SMmJsb2/gcGBmbIx8OHD+Xs2bOye/du2bp1q75JmEwyXVAdnTIr0FFRUVKzZk358ssvxWKxyM2bN6Vjx44yZswY2bhxo5w6dSpT/lOjQ4cOKc6JfvrppwLI0qVLnRqHqxLX/z///DPTvlJ6nx1JXDsJ701kRx4+fCi9evWSzp07y4MHD8wOxyVwpkDrQ2MTkC9fPrZv307z5s25cuUKx48fZ8+ePURGRjJrlnVDPetn4hxS+6o9efJkANq0aZNsGRHJsU9mxU0RPfbYY5n+HBYtWsSgQYPYuXMnHTt2dER4ibDffN5+RVB2JDIyksKFC1OiRAkCAwMpUKCA2SHlePRNwiSoUaMG+/btY/r06Vy/fp2wsDD2799Pw4YNWb58uVPbrl69ui0fGhrKnj17khSiHTt2JOsjT548OVag33jjDQB8fHwy7atdO+teXn379s20r6Q4e/Ysjz32GAD/+9//yJs3e4+HRo4cSd68eQkKCqJw4cJmh5MryN6/MU6kUaNGREdHx/ujunTpEo7a6yM5vv76a1ve3d3dlvfz87PdLQfStG5XKeXU0b6ZOGI0Wq1aNQCGDRuWaV8JuXPnDvXr1wesB6E6as22WWzYsIG1a9cSEBCgz9fMQvQIOgXixDlu6VBERITTBdpisSS5KVPTpk3jjYrHjRuXqMzSpUtz7MgZrEsW43DU5/DII4/E+6foKOwPPc2uW8w+ePCAH374gT59+jB06FCWLFlC1apVzQ4rV6EFOhUePHhA+/btcXNzY86cOTRq1Mip7SmlKFSokO0mgf3TjfZcv3493muLxcLAgQNtr9esWZPjRs/h4eG2/IQJExzi8/jx4w7xk5C2bdvy9NNPA1C6dGmntOEszp07x/jx42nYsCGLFy+mW7dunD17lp49e5odWu7DWXcfzUqZXcWRkMjISAHk9OnTDvWbHiwWiwwcONC2GmDChAni4eFhWxUQGxsbb+lZTubnn3+29XP//v2Z8jV48GABnLoaoUCBAi79mQQHB8uyZcts72nDhg0FkOHDh8tvv/0m0dHRZofo8qCX2Zkn0OHh4QJIcHCwQ/2ml3v37sUT4Lfeektat24tkZGR8tJLL9mu+fv7mxpnVuCof0a7du0SQH744QcHRZaYRYsWCSAjR46Un3/+2WntpAeLxSIbNmyQChUqxHsvAfnmm28y/Y8vt6EF2kSBnjNnjpQuXVpErFuOTpkyRWJjYx3aRlqpUaOG7Q8pNjZWevXqJUAi8c7J7Nixw9bXkJCQTPnKqvetXr16AoiXl5dT20kJPz8/qV+/vnTq1EkAKVmypKxcuVKioqJMiymnoAXaRIGOGwG9/PLLUqlSJQFk3bp1Dm0jPdiPdqKjo6VWrVq2r6VpEZq4KZvdu3enq92pU6em+4lIR7N27VpbP7t16+YQn/b/8JxF/fr1BZCIiAintZGQkJAQmTZtWqIRMiA7d+40bZCRE3GmQOubhKkwYMAA9u/fT7ly5ZgyZQpgXXJkFvbnGebLlw9vb+94a6dTI26zpieffNL2EIVSCqUU586dS1S+V69eKKV46623+Oqrr1i/fj2bNm1ixIgRtno3btwgNjYWgK1bt3Lw4MFM9DB5XnzxRcB6RuTGjRsd4vPPP/8E4NFHH+Xs2bNOOTPyzJkzdOrUKUvWDgcFBTFx4kRq1qzJoUOHmDVrFnv27CEmJsb2R9+hQ4ck983QuCDOUn6zkqNH0Alp0aKFAKZ/NcRuRLRp0yZ55JFHBJBr166luV7CtGnTJunRo4cMGTJEZs6cGW9u2z4ppSRfvnyJ7EuXLrXl4w46mDt3rsOmEOJ8O5LTp08n+348+eSTUqFCBalfv75Uq1ZNtm7dmqE9WgCZPHlymspmdA+YqKgoeffdd6VUqVLy8ssvZ+u9LbIb6CkO1xHor776SgBZuXKlU9tJC3Eb/gASHh5uO7WlS5cuyW6mtG3btngiNH36dBkyZIg89dRTNttjjz2WopADkj9/filZsqR88cUXyZYRcZyo/v777wJI8+bNM+0r4XsQl0aPHm37OX36dJk5c6Z4e3snKle1alWZOnWq3Lx5M02CCkjx4sVTnVaIey/TS0xMjHTu3Fm8vb3l8uXL6a6vyRxaoF1IoONu+ERGRjq1HRGRoKCgVDeitxeOQ4cOydWrV+Xjjz+WChUqyHPPPSfDhw9PsV7cH/SAAQMErJuui4hcunRJAHn66acTCdSsWbMEkKJFi8b7JzFkyBBb/rPPPpOoqCiHCXR0dLTNV0Z28vvrr78S9WPs2LHpOqDBYrHIgQMHZPr06Yl8jRs3Ltk5+rFjx6bpfcjoezVr1ix54oknTL9HkFsxXaCB+cAt4KSdrRSwDethr9uAkoZdAd9hPVj2ONDcrs5Qo/x5YKidvQXWg2H9jboqpTZSSs4W6KJFi8qkSZOc2kYc48aNS/Umlr1A2v9xJ9yW8+uvv45XLzg4WAoWLCiA3L9/X3r06CGAzJgxQ0REwsLCBKwHgkZHR9tuLiZMccKeUnLU2YELFy60+Zw0aZLcvHlT9u/fL3379pXvvvtOnn32WWnZsqXtvMjq1asnGY+vr2+mY7F/399///14/rt27Srh4eHxym/cuFEg+WWQcf48PDzSFcfOnTsFHLO7nyZjuIJAtweaJxDoL4C3jfzbwFQj3x3YZAh1a8BH/hHbC8bPkkY+TtQPGmWVUbdbSm2klJwp0EeOHBFAduzY4bQ27LHfdxeQxo0bJ1lu8uTJAsjBgwcTXYuLGbvRchxxJ5fbp0uXLonIP4Jx69YtWyxxD13Epfv379uWbaWUHLFi4OrVq/L666/b/pGkJdWuXVteeeUVmThxotPXsVssFlmyZEm8pZBx/8BiYmLk5s2bKY6Qk/tHmhIPHz601cvs/uWajGO6QFtjoHoCgT4HVDTyFYFzRn4OMCBhOWAAMMfOPsewVQTO2tlt5ZJrI6XkTIG+fPmyAPLJJ584rY2EpPaU4P79+8XDw8MmpMkxbtw4yZMnT7wnwxLOxYaGhsark1DwqlWrJteuXRNA6tevLyIioaGh8cpERETY8m3atMnQ0rLIyEhZs2aNzU+RIkUEkBIlSsjixYvFx8cnU2dHZgUrVqxI9h/HmDFj5Pr16/HKpyTeSXHw4EEBpF69evLuu+86OnxNOnBVgQ6xy6u418B6oK3dtR2AFzAB+D87+zuGzQvYbmdvB6xPqY0kYhsJ+AK+VatWdfDbHx+wrmLISuxHSn5+fvLXX3/JzZs3JSIiQho1aiQLFy5M1UdMTIy0adNGevbsKUuXLpWAgABZv359POGYPn26rXxAQEAiYVm8eLHtJmlcWrduXbwRa0a5f/++zJgxI0lBGzBggOzatSvDvs3EYrFIUFCQzJ8/39afiRMnStmyZWXVqlW2cukR6LjppmLFiumRswvg8gJtvL4rThRo+zZSSs6eg46bF46bCrh582aqo1dHkNxoDNK+5O/+/fvy9ddfS+/evaV06dLxfAwbNkzAerPL29tb3N3d5fHHH5e33npL2rZtK25ubtK0adNEbY8YMcI2Wk4LN27ckL/++ktERCIiIuLtrQFIu3btZMKECRl+n1yduG8wPj4+UqZMGds3l4kTJ6ZZoOPeKy3OroGrCnSum+IQsU45VK5cWbp37y5vvPGGAFK3bl2ntikicvTo0XhCdu/evUytJDl79qwUKVJEKlWqJMeOHZPIyEgZOXKkdO3aVVatWpXsnK3FYpEVK1ZI//79pXfv3nLv3j0ZPny4La6Eq0bs55/t4x8/fryUKVNGunXrJu+++67cv38/w33Jrly4cEGaN28uJUuWTPRNJjnxjSuzefPmLI5WkxyuKtBfEv8G3hdGvgfxbxIeNOylgItYbxCWNPKljGsJbxJ2T6mNlJKzBVpE4j1avX//fgHk/PnzGfY3efLkRGL7888/y7Bhw2TevHnxVkp8/vnniepv3rw50ci2T58+qY6wbt++bdvPIjIyMlM7l23ZssXWdsuWLcXf39+2muLgwYNy+PDhRDGuXbs2w+3lJC5fvizr16+XCxcu2N6bqlWrJtpSICgoSADZtm2bSZFqksJ0gQaWAjeAaOAaMAIojXX64jyw3U5sFTATCMC6dM7Lzs9wrEvp/IFhdnYv4KRRZwb/LLNLso2UUlYI9JUrV6RPnz6yb98+sVgs4unpKWvWrMmQL/ubatu3b5fFixfbNkECpHfv3vLDDz+Ir69vkoL7448/ChBvNcWYMWMEkI8++ijdX4MtFkum9oyIiIiQTz75JJEYJ1z2RzpuiOU2LBaL7Ny5UypVqiTe3t7y5ptvykcffSStW7eWwYMHmx2eJgGmC3R2Slkh0PbE3bBp27at7fFnd3d3GTp0aJp9xAlW+fLl5YknnpB///vfsnHjxlRHtAcPHpQ8efLIxIkT4/mJiIiQXbt2Sfny5eXrr79O1wMMoaGhcufOnTSXT46ES/iWL18uImJ78EXPn6bOnTt3ZOXKlfLee+/J0KFDZeHChYlW2mjMRwu0Cwu0iMizzz4rCxYskIsXL8rdu3elRYsW8vjjj6eprv1qibRisVhk37590rt3b3nllVdsdl9fXwHkP//5j4hIPIEcMmSIPHz4UCwWS7Lrks+cOePQzevv3r1ra9/f39/28Iunp6fD2tBozEYLtIsLtD0J1y0XKVJERo8eLQEBAUmW37Bhg4D18enUCA0NtT1mXK1aNXn99dfj7Yls/3TbgwcPZN++fbZH05NLmzdvlocPH4qIOGVt8e3bt/XUhiZHowU6Gwl0HJGRkXLw4EH5/vvv5dVXX5UyZcpIvXr1ZPDgwfLzzz/H20UtpRt6e/bskTp16sQTuEaNGiW7KU5CEfzyyy9ttps3b8r58+fl+vXrsn37dqlVq5btcejnnnvOadMO9t8SEj6godFkd7RAZ0OBTkhUVJRs2LBBmjRpkmhEmZQw3r9/X4YNGyZly5aVOXPmJDn6TQpAmjRpEs925MiRZOOyWCy2G42//fZbpvqo0eRGnCnQetfuLCJfvnx0794dPz+/RBvjh4WFMWnSJNsG+OXKlaNw4cKEh4dz/vx5Ro0alchfsWLFEtn+/vtvAGbNmhXP3qxZs2TjUkoxfPhwAJ5++mlCQ0PT3TeNRuMctECbQJ06dYiIiLC97tevH6tXr7a97tSpEytWrGD58uWsWrUqXt2QkBBiY2OTFNJt27YBUKtWrQzHNnv2bIefKKLRaDJG3HrjHIOXl5f4+vqaHUaauX37NuXKlWPnzp107NgRsB5rpZSylRGReK9T4sMPP2Ty5MnkzZs3XXHY+69RowYBAQFpblOjyc0opQ6LiJczfOsRtMkMHjwYgLZt29psCYUxPUL57rvvplucAaKjo1m7di1vv/02Fy9e5OHDh+n2odFoHIsWaJM4deoU3bt3Z8uWLfz222/ky5eP5557jlOnTpkST968edm8eTM///wzv//+OwULFjQlDo1G8w/pH2ppMs3x48eZOHEi9+/f58cff2ThwoVMmDCBc+fOcezYMZo3b85XX31FpUqVnB7L7t27Wbx4MXv27OH8+fMEBARQs2ZNp7er0WhSRwt0FnPu3Dn69evH+fPnAQgODua5557j+eefp3Tp0pQoUYJffvmFDh06cODAAUqWLOmUOMLCwmjSpAkXL17kq6++4oUXXiAiIkKLs0bjQmiBzmLq1q3LqVOnOHPmDDVq1EhyuVzz5s0JCQmhVatWbN++napVqzo0hqioKLp27UrTpk05duxYkjFoNBrz0XPQJpAvXz4eeeSRFIVx+vTpjBgxgnbt2vHhhx9y69Yth7X/1ltvcefOHVatWqXFWaNxYbRAuyhubm689dZbLFmyhKtXr+Lp6Ym3tzenT5/OtO9vvvmGCxcu6GV0Go2LowXaxWnTpg1z587l5s2bdO3alXbt2vHGG28QGxubKb/Tpk1zUIQajcZZaIHOJhQqVIjXX3+d8+fPc/ToUQYOHEh0dHSG/Xl7ezswOo1G4wy0QGczSpUqxaZNmwgJCaFjx46k96nJuP06Klas6IzwNBqNA9ECnQ0pWLAgq1evZsiQIXTv3p158+alqV5ERARlypRBKYW7u7uTo9RoNJklwwKtlKqrlPKzS6FKqXFKqfeVUoF29u52dSYppfyVUueUUl3t7N6GzV8p9badvYZSysewL1dK5c94V3MWRYoU4T//+Q9//PEHn3zyCR999FGS5USE+fPno5SiaNGiAPj4+GRlqBqNJoNkeB20iJwDmgIopdyAQOBXYBjwXxGJdxdKKdUA6A80BCoB25VSdYzLM4HOWA+kPaSUWicip4Gphq9lSqnZWA+rjb+XZi7H09OT33//nQ4dOlC9enUGDx7MvXv3OHToED4+PsyePZuwsDAmTJiAh4cH48aNMztkjUaTRhz1oEonIEBELqewdOsZYJmIPAQuKqX8gZbGNX8RuQCglFoGPKOUOgN0BAYaZRYC76MFOhEeHh7MnDmTLl268MYbb3Dnzh08PT2pW7cuU6dOZcCAAXpJnUaTDXGUQPcHltq9flkpNQTwBcaLyF3AAzhgV+aaYQO4msDeCigNhIhITBLl46GUGgmMBBz+1F12oXPnzvz5558cO3aMXr166ZuAGk0OINM3CY154V7ASsM0C6iFdfrjBvBVZttIDRH5QUS8RMSrbNmyzm7OZWndujWjRo3S4qzR5BAcMYLuBhwRkSCAuJ8ASqm5wHrjZSBQxa5eZcNGMva/gRJKqbzGKNq+vEaj0eR4HLHMbgB20xtKKfvhWx/gpJFfB/RXShVQStUAPIGDwCHA01ixkR/rdMk64zDGXUA/o/5QYK0D4tVoNJpsQaZG0EqpIlhXX9ifavqFUqop1pOnL8VdE5FTSqkVwGkgBhgrIrGGn5eBLYAbMF9E4natfwtYppT6GDgK/JiZeDUajSY7oc8k1Gg0mkygzyTUaDSaXIgWaI1Go3FRtEBrNBqNi6IFWqPRaFyUHHeTUCl1G7ichqJlgDtODscMcmq/QPctO5JT+wX/9K2aiDjlCbkcJ9BpRSnl66w7r2aSU/sFum/ZkZzaL8iavukpDo1Go3FRtEBrNBqNi5KbBfoHswNwEjm1X6D7lh3Jqf2CLOhbrp2D1mg0GlcnN4+gNRqNxqXRAq3RaDQuSrYUaKWUm1LqqFJqvfE6ycNlja1Nlxt2H6VUdTsfLnmAbRJ9W6CUumh3CG9Tw66UUt8Z8RxXSjW38zFUKXXeSEPt7C2UUieMOt8p4xwspVQppdQ2o/w2pVRJJ/TrktG2n1LKN6V2s1PfkumX0w9OTul324F9K6GU+kUpdVYpdUYp9VhO+MxS6JvrfW4iku0S8AawBFhvvF4B9Dfys4ExRv4lYLaR7w8sN/INgGNAAaAGEIB1q1M3I18TyG+UaZBSG1nQtwVAvyTKdQc2AQpoDfgY9lLABeNnSSNf0rh20CirjLrdDPsXwNtG/m1gqhP6dQkok8CWZLvZqW/J9Ot9YEISZR32e5fc77aD+7YQ+LeRzw+UyAmfWQp9c7nPzeEC4+yE9WSVHVgPlF1vfLh3gLzG9ceALUZ+C/CYkc9rlFPAJGCSnc8tRj1bXcM+yUjJtuHMvhm2BSQt0HOAAXavzwEVsR6gMCdhOePaWTu7rVxcXSNfETjnhL5dIrGQJdludupbMv16n6T/0B32e5fc77YD++UOXEzoM4d8Zsn1zeU+t+w4xfENMBGwGK9TOlzWA+NAWuP6PaO8zZ6gTnL2NB9gm0m+IX7f4vjE+Nr4X6VUAcOW3j54GPmEdoDyInLDyN8EymeyH0khwFal1GFlPeQ3pXazU9+S6hdYD04+rpSab/cV3ZG/d8n9bjuKGsBt4CdlnXKbp6wHdOSEzyy5voGLfW7ZSqCVUj2BWyJy2OxYHE0KfZsE1AMexfo18S1nxiHWf+3OWHvZVkSaYz3DcqxSqn0WtevsNpLqV5YfnOwE8gLNgVki0gyIwDrdYCMbf2bJ9c3lPrdsJdBAG6CXUuoSsAzrVMC3GIfLGmXsD5e1HVRrXHfHehhtcgfYJme3HWCbRBuOIlHflFKLROSGWHkI/AS0TNi3NPYh0MgntAMEKeMsSePnLUd2DEBEAo2ft4BfsfYjuXazTd+S6peIBIlIrIhYgLlk/DNL6fcuud9tR3ENuCYiPsbrX7CKWrb/zEimb674uWUrgRaRSSJSWUSqY51g3yki/yL5w2XXGa8xru80/iO73AG2yfRtkN0vqgJ6E/8Q3iHG3fPWwD3ja+EWoItSqqTxFa0L1vmvG0CoUqq14WtIMu+Tw/umlCqilCoWlzdiOplCu9mib8n1S2XNwcnJ/W47BBG5CVxVStU1TJ2wniearT+zlPrmkp+bIyffszIBT/LPjbSaxhvmD6wEChj2gsZrf+N6Tbv6U7DegT2HcffYsHcH/jKuTbGzJ9lGFvRtJ3DC+GVZBBQ17AqYacR5AvCyqz/ciNMfGGZn9zL8BAAz+OdJ0tJYb06eB7YDpRzcn5pY73AfA07Fva/JtZtd+pZCv3424j6O9Q+yoqN/71L63XZg/5oCvkY/1mBdhZGtP7NU+uZyn5t+1Fuj0WhclGw1xaHRaDS5CS3QGo1G46JogdZoNBoXRQu0RqPRuChaoDUajcZF0QKt0Wg0LooWaI1Go3FR/h/8rkzMZrp04QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ax = uk_SE.plot(color='white', edgecolor='black')\n",
    "gdf.plot(ax=ax, color='red', markersize=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "26293199-135e-4813-8fb4-9f220f1fdc4f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "id                                          http://environment.data.gov.uk/water-quality/d...\n",
       "sample.samplingPoint                        http://environment.data.gov.uk/water-quality/i...\n",
       "sample.samplingPoint.notation                                                     SO-19801A08\n",
       "sample.samplingPoint.label                                             APSLEY FARM BOREHOLE A\n",
       "sample.sampleDateTime                                                     2000-03-22 10:30:00\n",
       "determinand.label                                                                     BOD ATU\n",
       "determinand.definition                                                        BOD : 5 Day ATU\n",
       "determinand.notation                                                                       85\n",
       "resultQualifier.notation                                                                    <\n",
       "result                                                                                    3.0\n",
       "codedResultInterpretation.interpretation                                                  NaN\n",
       "determinand.unit.label                                                                   mg/l\n",
       "sample.sampledMaterialType.label                                                      UNCODED\n",
       "sample.isComplianceSample                                                               False\n",
       "sample.purpose.label                                 WASTE MONITORING (AGENCY AUDIT - PERMIT)\n",
       "sample.samplingPoint.easting                                                         442350.0\n",
       "sample.samplingPoint.northing                                                          146570\n",
       "yr                                                                                       2000\n",
       "month                                                                                       3\n",
       "day                                                                                        22\n",
       "geometry                                       POINT (-1.3950297049434972 51.216766219183306)\n",
       "Name: 0, dtype: object"
      ]
     },
     "execution_count": 66,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf.iloc[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "id": "3d2c5de6-f301-4eb5-b53f-89d4a2be8d62",
   "metadata": {},
   "outputs": [],
   "source": [
    "df = data.groupby(by='sample.samplingPoint.label').head(1)\n",
    "geometry = [Point(xy) for xy in zip(df['sample.samplingPoint.easting'][:],\n",
    "                                    df['sample.samplingPoint.northing'][:])]\n",
    "map_locations = GeoDataFrame(df['sample.samplingPoint.label'][:], geometry=geometry, crs=27700).to_crs(4326)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "id": "flying-smile",
   "metadata": {
    "scrolled": true
   },
   "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 src=\"about:blank\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" data-html=%3C%21DOCTYPE%20html%3E%0A%3Chead%3E%20%20%20%20%0A%20%20%20%20%3Cmeta%20http-equiv%3D%22content-type%22%20content%3D%22text/html%3B%20charset%3DUTF-8%22%20/%3E%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%3Cscript%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20L_NO_TOUCH%20%3D%20false%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20L_DISABLE_3D%20%3D%20false%3B%0A%20%20%20%20%20%20%20%20%3C/script%3E%0A%20%20%20%20%0A%20%20%20%20%3Cstyle%3Ehtml%2C%20body%20%7Bwidth%3A%20100%25%3Bheight%3A%20100%25%3Bmargin%3A%200%3Bpadding%3A%200%3B%7D%3C/style%3E%0A%20%20%20%20%3Cstyle%3E%23map%20%7Bposition%3Aabsolute%3Btop%3A0%3Bbottom%3A0%3Bright%3A0%3Bleft%3A0%3B%7D%3C/style%3E%0A%20%20%20%20%3Cscript%20src%3D%22https%3A//cdn.jsdelivr.net/npm/leaflet%401.6.0/dist/leaflet.js%22%3E%3C/script%3E%0A%20%20%20%20%3Cscript%20src%3D%22https%3A//code.jquery.com/jquery-1.12.4.min.js%22%3E%3C/script%3E%0A%20%20%20%20%3Cscript%20src%3D%22https%3A//maxcdn.bootstrapcdn.com/bootstrap/3.2.0/js/bootstrap.min.js%22%3E%3C/script%3E%0A%20%20%20%20%3Cscript%20src%3D%22https%3A//cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js%22%3E%3C/script%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//cdn.jsdelivr.net/npm/leaflet%401.6.0/dist/leaflet.css%22/%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//maxcdn.bootstrapcdn.com/bootstrap/3.2.0/css/bootstrap.min.css%22/%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//maxcdn.bootstrapcdn.com/bootstrap/3.2.0/css/bootstrap-theme.min.css%22/%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//maxcdn.bootstrapcdn.com/font-awesome/4.6.3/css/font-awesome.min.css%22/%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css%22/%3E%0A%20%20%20%20%3Clink%20rel%3D%22stylesheet%22%20href%3D%22https%3A//cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css%22/%3E%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cmeta%20name%3D%22viewport%22%20content%3D%22width%3Ddevice-width%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20initial-scale%3D1.0%2C%20maximum-scale%3D1.0%2C%20user-scalable%3Dno%22%20/%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cstyle%3E%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23map_20998f52f0014fe5bdf4336ead0d2db7%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20position%3A%20relative%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20width%3A%20100.0%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20height%3A%20100.0%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20left%3A%200.0%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20top%3A%200.0%25%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%3C/style%3E%0A%20%20%20%20%20%20%20%20%0A%3C/head%3E%0A%3Cbody%3E%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%3Cdiv%20class%3D%22folium-map%22%20id%3D%22map_20998f52f0014fe5bdf4336ead0d2db7%22%20%3E%3C/div%3E%0A%20%20%20%20%20%20%20%20%0A%3C/body%3E%0A%3Cscript%3E%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20map_20998f52f0014fe5bdf4336ead0d2db7%20%3D%20L.map%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22map_20998f52f0014fe5bdf4336ead0d2db7%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20center%3A%20%5B50.82%2C%20-1.01%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20crs%3A%20L.CRS.EPSG3857%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zoom%3A%2010%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20zoomControl%3A%20true%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20preferCanvas%3A%20false%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29%3B%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20tile_layer_6c34ab717fb44cd0806ef11df99c7b60%20%3D%20L.tileLayer%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22https%3A//%7Bs%7D.tile.openstreetmap.org/%7Bz%7D/%7Bx%7D/%7By%7D.png%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22attribution%22%3A%20%22Data%20by%20%5Cu0026copy%3B%20%5Cu003ca%20href%3D%5C%22http%3A//openstreetmap.org%5C%22%5Cu003eOpenStreetMap%5Cu003c/a%5Cu003e%2C%20under%20%5Cu003ca%20href%3D%5C%22http%3A//www.openstreetmap.org/copyright%5C%22%5Cu003eODbL%5Cu003c/a%5Cu003e.%22%2C%20%22detectRetina%22%3A%20false%2C%20%22maxNativeZoom%22%3A%2018%2C%20%22maxZoom%22%3A%2018%2C%20%22minZoom%22%3A%200%2C%20%22noWrap%22%3A%20false%2C%20%22opacity%22%3A%201%2C%20%22subdomains%22%3A%20%22abc%22%2C%20%22tms%22%3A%20false%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_5eaaaae2a542491c8de9302ca109d631%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.83378852730109%2C%20-0.9817906946486339%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_56febd90a4ec46af9d4d8ea82bb106d1%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_c52b04e08e364869953b2c141418b601%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_c52b04e08e364869953b2c141418b601%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20BRIDGE%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_56febd90a4ec46af9d4d8ea82bb106d1.setContent%28html_c52b04e08e364869953b2c141418b601%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_5eaaaae2a542491c8de9302ca109d631.bindPopup%28popup_56febd90a4ec46af9d4d8ea82bb106d1%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_938251ef9faf49d8a6aa4ca2ac154b2a%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.8294403545175%2C%20-0.9989243345818196%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_45e14ebbd54945d2b2d2297ae68fd63a%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_87941088c09d480eb1ea824df8f6b281%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_87941088c09d480eb1ea824df8f6b281%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20CHAN%20SOUTHMOOR%20POINT%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_45e14ebbd54945d2b2d2297ae68fd63a.setContent%28html_87941088c09d480eb1ea824df8f6b281%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_938251ef9faf49d8a6aa4ca2ac154b2a.bindPopup%28popup_45e14ebbd54945d2b2d2297ae68fd63a%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_58e40676704143399729a240ed22f2f7%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.825369960659515%2C%20-0.9961717708797697%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_d425bcfd84c4419ab338c2d36a725f6d%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_9e2294621b8f423ab649aebf9de961c4%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_9e2294621b8f423ab649aebf9de961c4%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20CHANNEL%20THE%20DEEPS%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_d425bcfd84c4419ab338c2d36a725f6d.setContent%28html_9e2294621b8f423ab649aebf9de961c4%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_58e40676704143399729a240ed22f2f7.bindPopup%28popup_d425bcfd84c4419ab338c2d36a725f6d%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_36696c4a5f844c68a50cb12188dd1f7e%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.79280810321534%2C%20-1.0270898893757867%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_1264b96f770842c4907f414bc866e69b%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_e8924bd8f59146c7ba924853e96e853a%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_e8924bd8f59146c7ba924853e96e853a%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20HARBOUR%20MOUTH%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_1264b96f770842c4907f414bc866e69b.setContent%28html_e8924bd8f59146c7ba924853e96e853a%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_36696c4a5f844c68a50cb12188dd1f7e.bindPopup%28popup_1264b96f770842c4907f414bc866e69b%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_cab9f40875f44c56870525f55c218f6d%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.82566164409545%2C%20-0.9987211691148886%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_92537ad18d81436eb36fbb5722584163%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_f002a1777107466ab7e800b484a11dd7%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_f002a1777107466ab7e800b484a11dd7%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20HARBOUR-SFW%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_92537ad18d81436eb36fbb5722584163.setContent%28html_f002a1777107466ab7e800b484a11dd7%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_cab9f40875f44c56870525f55c218f6d.bindPopup%28popup_92537ad18d81436eb36fbb5722584163%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_46f33c61b0be4b118c83b1c5b893af31%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.80801832606203%2C%20-1.0178326655699979%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_ea56fb872a874dbeaffb13d3a856dfb9%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_e72ff5e5f1aa4ca7ab7fb91a7ca75288%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_e72ff5e5f1aa4ca7ab7fb91a7ca75288%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ETHE%20RAFTS%20LANGSTONE%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_ea56fb872a874dbeaffb13d3a856dfb9.setContent%28html_e72ff5e5f1aa4ca7ab7fb91a7ca75288%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_46f33c61b0be4b118c83b1c5b893af31.bindPopup%28popup_ea56fb872a874dbeaffb13d3a856dfb9%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_e34c80509e5947de967a98355d8b4746%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.81329273275425%2C%20-1.0035277722931057%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_95cc1b481b3741aa9b58723936f02c52%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_6869e1ee3804410cb31a20aa15e071c7%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_6869e1ee3804410cb31a20aa15e071c7%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3EALPHA%20BUOY%20LANGSTONE%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_95cc1b481b3741aa9b58723936f02c52.setContent%28html_6869e1ee3804410cb31a20aa15e071c7%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_e34c80509e5947de967a98355d8b4746.bindPopup%28popup_95cc1b481b3741aa9b58723936f02c52%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_15719a135e5a4364b5651b8cb7273bd2%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.83160778860656%2C%20-1.0428975191776155%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_b9ff3a1f6f9f4caebdff15a03a377c21%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_a543c16fca9c415e84d7c28aaaa87b79%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_a543c16fca9c415e84d7c28aaaa87b79%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3EA2030%20BRIDGE%2C%20LANGSTONE%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_b9ff3a1f6f9f4caebdff15a03a377c21.setContent%28html_a543c16fca9c415e84d7c28aaaa87b79%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_15719a135e5a4364b5651b8cb7273bd2.bindPopup%28popup_b9ff3a1f6f9f4caebdff15a03a377c21%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_66232914e0d44aaeadb5e1ab5cf62efb%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.82032203539504%2C%20-1.037449629213488%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_16be5f71222e4b7bbd5e235e05767005%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_ac4c9f65e0a44652ae0c38838fb8bed3%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_ac4c9f65e0a44652ae0c38838fb8bed3%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3EBROOM%20CHANNEL%20LANGSTONE%20HARBR%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_16be5f71222e4b7bbd5e235e05767005.setContent%28html_ac4c9f65e0a44652ae0c38838fb8bed3%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_66232914e0d44aaeadb5e1ab5cf62efb.bindPopup%28popup_16be5f71222e4b7bbd5e235e05767005%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_2e20ea7a437e4f6ca184107474f4a42a%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.80625583306303%2C%20-1.0221272615526298%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_15f5d04fe4fc459da48e1018870fb01e%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_c10d9c52590a477f948b2379129af09c%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_c10d9c52590a477f948b2379129af09c%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ENW%20SINAH%20BUOY%2C%20LANGSTONE%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_15f5d04fe4fc459da48e1018870fb01e.setContent%28html_c10d9c52590a477f948b2379129af09c%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_2e20ea7a437e4f6ca184107474f4a42a.bindPopup%28popup_15f5d04fe4fc459da48e1018870fb01e%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_91f5f086458c4b6aa54972d37f432294%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.78659306073791%2C%20-1.0367236125316468%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_d185034a62224efb9f4cedd2415b171e%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_3aeda91d02d34ec59a92b0d6e82d1ff1%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_3aeda91d02d34ec59a92b0d6e82d1ff1%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3EDUMP-INV%20LANGSTONE%20SED-%2ACLOSED%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_d185034a62224efb9f4cedd2415b171e.setContent%28html_3aeda91d02d34ec59a92b0d6e82d1ff1%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_91f5f086458c4b6aa54972d37f432294.bindPopup%28popup_d185034a62224efb9f4cedd2415b171e%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_bbb9717b1d444c2e80c3abb133b9c0e2%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.78659306073791%2C%20-1.0367236125316468%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_94975ab087d245b9b1fb3597e5e3f59f%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_b82b17cec3914286bd204b00484208ab%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_b82b17cec3914286bd204b00484208ab%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3EDUMP-INV%20LANGSTONE%20WEED-CLOSED%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_94975ab087d245b9b1fb3597e5e3f59f.setContent%28html_b82b17cec3914286bd204b00484208ab%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_bbb9717b1d444c2e80c3abb133b9c0e2.bindPopup%28popup_94975ab087d245b9b1fb3597e5e3f59f%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_709fa8bc7c0e40b8afe5372420b35d5c%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.809071101154316%2C%20-1.0351544537222745%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_d834ca1b719b428bbc22ade46ce7d558%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_01300651c0524e4b99799551fe8cfeb9%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_01300651c0524e4b99799551fe8cfeb9%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20HARBOUR%20CSEMP%20MUSSELS%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_d834ca1b719b428bbc22ade46ce7d558.setContent%28html_01300651c0524e4b99799551fe8cfeb9%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_709fa8bc7c0e40b8afe5372420b35d5c.bindPopup%28popup_d834ca1b719b428bbc22ade46ce7d558%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20marker_dc91e9c8d5cb49a38b83977b9ecdca45%20%3D%20L.marker%28%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B50.8178727312132%2C%20-1.0133673131025998%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%29.addTo%28map_20998f52f0014fe5bdf4336ead0d2db7%29%3B%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%20%20%20%20%20%20%20%20var%20popup_9c140df0707844e8afdffbbdfcff68c7%20%3D%20L.popup%28%7B%22maxWidth%22%3A%20%22100%25%22%7D%29%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20var%20html_9ac202bc088f46ff834a8cd8b0ca1a55%20%3D%20%24%28%60%3Cdiv%20id%3D%22html_9ac202bc088f46ff834a8cd8b0ca1a55%22%20style%3D%22width%3A%20100.0%25%3B%20height%3A%20100.0%25%3B%22%3ELANGSTONE%20HARBOUR%20WFD%20BI%3C/div%3E%60%29%5B0%5D%3B%0A%20%20%20%20%20%20%20%20%20%20%20%20popup_9c140df0707844e8afdffbbdfcff68c7.setContent%28html_9ac202bc088f46ff834a8cd8b0ca1a55%29%3B%0A%20%20%20%20%20%20%20%20%0A%0A%20%20%20%20%20%20%20%20marker_dc91e9c8d5cb49a38b83977b9ecdca45.bindPopup%28popup_9c140df0707844e8afdffbbdfcff68c7%29%0A%20%20%20%20%20%20%20%20%3B%0A%0A%20%20%20%20%20%20%20%20%0A%20%20%20%20%0A%3C/script%3E onload=\"this.contentDocument.open();this.contentDocument.write(    decodeURIComponent(this.getAttribute('data-html')));this.contentDocument.close();\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7fe90fcef5b0>"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "langstone = [loc for loc in map_locations['sample.samplingPoint.label'] if 'LANGSTONE' in loc]\n",
    "langstone_locations = map_locations[map_locations['sample.samplingPoint.label'].isin(langstone)]\n",
    "m = folium.Map(location=[50.82, -1.01], tiles=\"OpenStreetMap\", zoom_start=10)\n",
    "for name, geom in zip(langstone_locations['sample.samplingPoint.label'][:], langstone_locations['geometry'][:]):\n",
    "   folium.Marker(\n",
    "      location=[geom.y, geom.x],\n",
    "      popup=name,\n",
    "   ).add_to(m)\n",
    "m"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "id": "swedish-botswana",
   "metadata": {},
   "outputs": [],
   "source": [
    "m.save('langstone-sites-map.html')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "incomplete-weather",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (Data Languages Lecture)",
   "language": "python",
   "name": "datalanguages"
  },
  "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.9"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {},
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}