scripts/exploratory-analysis-mdned.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "be69aa8d",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import seaborn as sns\n",
    "import os"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d3bb35ae",
   "metadata": {},
   "outputs": [],
   "source": [
    "dir_geoparsed_results = os.path.dirname(os.getcwd())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "98abd712",
   "metadata": {},
   "source": [
    "## Reading Toponym Resolution Results of Edinburgh Geoparser"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "e370fc45",
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>lgl_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Alexandria</td>\n",
       "      <td>31.3113</td>\n",
       "      <td>-92.4451</td>\n",
       "      <td>Alexandria 31.3113,-92.4451</td>\n",
       "      <td>Alexandria</td>\n",
       "      <td>31.3113 , -92.4451</td>\n",
       "      <td>31.2928 , -92.4592</td>\n",
       "      <td>Alexandria 31.3113,-92.4451</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Rapides Parish</td>\n",
       "      <td>31.1669</td>\n",
       "      <td>-92.4835</td>\n",
       "      <td>Rapides Parish 31.1669,-92.4835</td>\n",
       "      <td>Rapides Parish</td>\n",
       "      <td>31.1669 , -92.4835</td>\n",
       "      <td>31.2 , -92.54</td>\n",
       "      <td>Rapides Parish 31.1669,-92.4835</td>\n",
       "      <td>[(31.36351, -92.58681)]</td>\n",
       "      <td>18.717774</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Avoyelles Parish</td>\n",
       "      <td>31.0669</td>\n",
       "      <td>-91.9668</td>\n",
       "      <td>Avoyelles Parish 31.0669,-91.9668</td>\n",
       "      <td>Avoyelles Parish</td>\n",
       "      <td>31.0669 , -91.9668</td>\n",
       "      <td>31.07 , -92.0</td>\n",
       "      <td>Avoyelles Parish 31.0669,-91.9668</td>\n",
       "      <td>[(31.07624, -92.00138), (31.07624, -92.00138)]</td>\n",
       "      <td>0.706195</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Cottonport</td>\n",
       "      <td>30.9841</td>\n",
       "      <td>-92.0535</td>\n",
       "      <td>Cottonport 30.9841,-92.0535</td>\n",
       "      <td>Cottonport</td>\n",
       "      <td>30.9841 , -92.0535</td>\n",
       "      <td>30.98408 , -92.05346</td>\n",
       "      <td>Cottonport 30.9841,-92.0535</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>490</td>\n",
       "      <td>Pineville</td>\n",
       "      <td>31.3224</td>\n",
       "      <td>-92.4343</td>\n",
       "      <td>Pineville 31.3224,-92.4343</td>\n",
       "      <td>Pineville</td>\n",
       "      <td>31.3224 , -92.4343</td>\n",
       "      <td>31.3224 , -92.4343</td>\n",
       "      <td>Pineville 31.3224,-92.4343</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0              name      lat      lon  \\\n",
       "0           0        Alexandria  31.3113 -92.4451   \n",
       "1           1    Rapides Parish  31.1669 -92.4835   \n",
       "2           2  Avoyelles Parish  31.0669 -91.9668   \n",
       "3           3        Cottonport  30.9841 -92.0535   \n",
       "4         490         Pineville  31.3224 -92.4343   \n",
       "\n",
       "                             lgl_key            toponym             old_coord  \\\n",
       "0        Alexandria 31.3113,-92.4451        Alexandria    31.3113 , -92.4451    \n",
       "1    Rapides Parish 31.1669,-92.4835    Rapides Parish    31.1669 , -92.4835    \n",
       "2  Avoyelles Parish 31.0669,-91.9668  Avoyelles Parish    31.0669 , -91.9668    \n",
       "3        Cottonport 30.9841,-92.0535        Cottonport    30.9841 , -92.0535    \n",
       "4         Pineville 31.3224,-92.4343         Pineville    31.3224 , -92.4343    \n",
       "\n",
       "               new_coord                            mlg_key  \\\n",
       "0     31.2928 , -92.4592        Alexandria 31.3113,-92.4451   \n",
       "1          31.2 , -92.54    Rapides Parish 31.1669,-92.4835   \n",
       "2          31.07 , -92.0  Avoyelles Parish 31.0669,-91.9668   \n",
       "3   30.98408 , -92.05346        Cottonport 30.9841,-92.0535   \n",
       "4     31.3224 , -92.4343         Pineville 31.3224,-92.4343   \n",
       "\n",
       "                        geocoded_coordinates_list  median_error_distance  \n",
       "0                                              []                    NaN  \n",
       "1                         [(31.36351, -92.58681)]              18.717774  \n",
       "2  [(31.07624, -92.00138), (31.07624, -92.00138)]               0.706195  \n",
       "3                                              []                    NaN  \n",
       "4                                              []                    NaN  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lgl_edinburgh = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/lgl_geocoded_results_Edinburgh_Geoparser.csv')\n",
    "df_lgl_edinburgh.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "b4c28cec",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAAF3CAYAAAAGv5jnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAogElEQVR4nO3debwcVZ338c/XhLDvEYQkSFBwWAYBI5sDBkEFZh7Q0YdFXCJIQERFEYbFUcB5VBYXFJRNQBhWETUiCqKgI8omypKwhRggAYWw6LAn8Hv+OOeSTtF9u25u97n3dr7v16tft+vU6apfVXf/bvWpU6cUEZiZWRmvGeoAzMyWJE66ZmYFOemamRXkpGtmVpCTrplZQU66ZmYFOenaoEi6TtLsdmUDXOYUSSFp8uCi664c47lDHYeNLE66NUmanL9kn6tZfxlJB0n6taTHJM2X9JSkmyUdL+mfKvX7Es37u7MFr4qvb32tHjNLxGGJpDdLOlPSvZKekfS8pAckXS7pQ5KWGuoYrTNGD3UAvUjSesAVwIbAb4BvAI8AKwCbAfsCn5O0TkTMHao4s28BNzcp/99BLPNdgAbx+iWKpKOALwFPApcAdwIvAuOBnYDzgH8BDhiqGK1znHQ7TNKywM+ANwD/HhE/alJnGeAzwHC4HPB/IuKyTi4wIl7s5PK6LR9FjoqI54dg3VOA/wf8CnhfRPy9UuVYSZsBbyscWtdJErB8RDzdi+trxc0Lnfcx4J+AE5slXICIeD4ivhIRD5cNbfFJWjX//J2Xf/5eJ+ktLeq2bOeVtLakiyQ9KelZSVdJ2qDFakdLOib/zH5B0u2S9mqyvqZtq83ahvPyQtLGkr4uaQ7wPLB1nr+upB9K+kd+/ETSxBz7dS22dxtJv8n75XFJZ0laocU2Nb5uDPBl0q+KPZokXAAi4s8RcWqT1+8k6ercbPV83j8HtljXeyRdn2N8Oj/fvUm92fm92iI3jT0t6QlJ35e0RqXuipL+S9KN+XPxgqSZkr4qablK3b7muSmSPiFpBmm/fy7P31jSDyTNzcv5q6RrJf1rZTlLSzpK0vS8zU9J+qmkzQeyvqHkI93O62uTPWtIo6hvRUljm5Q/FxHPwCtHglcBbwXOB24gNZNcAzw+gHUtD/w2v/4oYCLwaeAnkjaJiJcq9Y/Pr/lOnv4ocJGkZSLi3AGst5kLgOeAr5F+cTwiaXXgf4A1gdOAu4DtgGtzHM1sRmpKOge4EJgM7Ae8DExtE8PbgLWA8yLiiYEEL2lqjvEG0pHyM8A7ge9KekNEHNZQ9yDgVOBu4LhcPAX4saQDIuKMyuLHk468fwhcBmxBahKbJOmtEfFsrjeOdJDxw7ztC4C3A4cDmwPvbhL6IcDqwJnAX4GH8n7/dZ5/GvAAMBaYBGxF+uXY9zn8BbAt6XN4CrAysD9wvaTtI+KWdutrElNZEeFHjQfpyxTA59rUexz4e5PyUaQPUuNj2Yb5U/Ly319oe/rW1+pxSkPdqbns2MoyDsnlsyvl17UoC+DwSvlhufzdTWJ7AFi5oXzlXPZEZd8FcG4/2zi5oeyYXHYdMLpS/4Q8b58W5ddVyoOUXLeqlP8MmA+s0OY9+GRexmeazFux8llZvWHeWqSjtgubvO5k4CVgvTy9KvA0MBNYqaHeSsD9pKPsVRrKZ+eYDqkst6857IiGsjHAUk1i+FKuu2WT788TwBqV+rvleXu02V99Mby7Ur4S8GDj+9Pf+ob64eaFzlsJ+EeT8g2BxyqPTxSMq5XjSEdI1ce3G+q8h/RF/lrltd+l+ba28jLpxF2jviOc9ZvU/240/OTOz08jJZLJA1hvM9+MiAWVsv9DOuF5UaX8pH6W84eIuLFS9mvSr8h128SwUv7bbB+ew6KflQca5r0fWBr4nqSxjQ/gp6Rmw51y3XeSjtK/FRGvrCc//xbp5O5OLOofLPx10ec7ufy9Dct4MSLmA0ganZugxpJ+AUE6Sq06LyIerZT1vce7SFqp+oIGHyQdrf+xss1jgF8C/6J0TqXd+oaUmxc67x8s/DI1+gvpCwDwZvr/Ig9Ibj+stiE+EfVOaN0REde0qbMe8EjjlxYgIl6QNIuUBOt4OF59sqqveWL1JvXvalI2oyGmwbi3SdlE4KaIeLmxMCIelfRUi+XMalLW3zY16tufzT4vx5L+wUD6Z/eGhnkb5r/9vW9r5r8T89/pTer0lVX35azqZ6fhvV6kbm66OBDYmFefI2r2uXjVfo+I30g6j/TLZB9JN5O27ZKImNFQdUNgWdI/oVbGsmgTQrP3eUg56XbencD2kiZGxF/6CiO1j14DIKl6hDVYnwO+WCnbgfQTejipttk26kYXs/4+38/2M28gBrNNd+a/m1VnRMQdwB0Akp5ssdwPk47Mm2n2z6CjJH2W9A/hatJR88Okrm7jgHNpfqK+6X6PiI9IOhHYhdSOfihwtKRDIuKUvlWS9sln+wmrmpA79T53jJNu510GbE86wXB0oXWeB/yuUnZbB5c/C3iXpJUaj3YlLU068qkmhU7ZEPhJpWyjhpj6PAGs1uT1Az0ang28UdJrGo9281n7VQa4rDquJ53cea+kz0ZE3ZOS9+W/82r8SunbTxuTTo41arYvAdaTNKbxaLfhvb67od6HSPtsl8r+2rntFjQREXeS/hGdKGkV4Ebgq5JOjdRQex/wWuDX1V8jI4nbdDvvLNIH8zBJ721Rp6NHdRExKyKuqTw6mQh/QjoReGil/OM0/2ncKR+XtHLfRH5+IPAU6aKTPvcC2zR2U5K0Kqm3w0D8lHSSau9KeVe6GeWkdhTppNkljdtaUf28XAq8QOrDW23DRNLKOUlCaut8BvikpBUb6qxIOpH3dK7TaCXgoErZQbn8xw1lL5FOVr0Sn6TRwBEttqMpSatJWiQXRcRTpCa55YBlcvF5wOtocaQrac1m5cONj3QHbkelixuq5kXEaRHxXO5beAVwuVLfzqtJRzQrkfrw7kn6wDbrvvI+VS4Rzu6PiOoJnk7YrsX2AFyQjzDOIfVg+IKkicAfSF2C/i/pDHi3PkfzgBslnZOnPwqsA3wsFnZbgtR16L+BX0s6n3RUuj/p5NPrBrC+44EPAOdI2pL0z3M7UheleXThYpaIOEfS2qQz/vdL6rsibT7pH8C7SFej3dnwmjmSPk76B39X3uYHSEeB/0w68bkRqQfJU5IOJ3UZu1EL+zNPAd4IHBCv7h98P/BFSZsAfwTeQuoydjeLngi9DPgK8HNJl5M+3x/IsQ/Eh4HPSPoRqZfFfFLXs3cDl0bEc7neyaTzIidKegfphOU/SJ+JHUk9OnYY4LrLG+ruEyPlwcIuKK0ed1fqL0vqnXAt6Qs7n3SEdgupC9KbKvWntFn+Lzq8Pe3WFzR0qSL9fP8e6STRM6T24km07h7WtiyXr5vXdUyT2HYinVB6kHRkdwfwgRbbcxgp8bxAOgG3L/13GVu3xXImApeTulL9g3SUPzG/h1dW6tbuqlbj/diMlETvI7VDPp+3+0eks/bNuma9Lc9/lNSW+nD+vB0KLFOp+17g9/m9eyY/f0+TZc7O79UWpKT2DKn56HxgzUrdUcCRpET5Qt7/J5Caharv6eRcNqXFtn8/L+eZvN9vy9uxdKXuaOBTpEvX+7blPlK/63fVWd9QP5QDNLMWcuf9ecDpEdH0iq9eoXQl4eyImDzEofQst+maNWjWRsrCNspq26fZgBVr05V0NvBvwKMRsUmT+SK12exK+nk1JSJuLRWfWXalpAeAW0kHJTuSPre/Z9GTSGaLpeSR7rlAf11JdiFdlbQ+6aTNdwvEZFZ1Bekk4ZdI7ZMbk/qi7hyvHhvCbMCKtulKWhe4osWR7umka6cvytP3kE5EtOr8bWY24gynLmPjWLQL1Zxc9qqkm0dYmgqw0UYbvWX69GZXOJqZDZmWffFH5Im0iDgjIiZFxKRll2123sPMbHgaTkl3LjChYXp8LjMz6xnDKelOAz6sZGvSmLRuzzWznlKyy9hFpKtExirdIuWLwFIAEXEacCWpu9hMUpexgV43b2Y27BVLuhFRHUSkOj8YHoN6m5l1zXBqXjAz63lOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkV5KRrZlaQk66ZWUFOumZmBTnpmpkVpIgY6hgGZcyYMTF//vwBvWbt8ROY+9CDXYrIzAy1nDHSk66k2PP03w/oNZccsC0jfbvNbFhrmXTdvGBmVpCTrplZQU66ZmYFOemamRXkpGtmVpCTrplZQU66ZmYFOemamRXkpGtmVpCTrplZQU66ZmYFFU26knaWdI+kmZKOaDJ/HUnXSvqTpNsl7VoyPjOzbiuWdCWNAk4FdgE2AvaWtFGl2ueBSyNic2Av4Dul4jMzK6Hkke6WwMyImBURLwIXA7tX6gSwUn6+MvBwwfjMzLpudMF1jQMeapieA2xVqXMMcLWkTwLLAzuVCc3MrIzhdiJtb+DciBgP7AqcL+lVMUqaKukWSbcUj9DMbBBKJt25wISG6fG5rNF+wKUAEfEHYBlgbHVBEXFGREyKiElditXMrCtKJt2bgfUlTZQ0hnSibFqlzoPAjgCSNiQl3ccKxmhm1lXFkm5ELAAOBq4C7iL1Upgu6ThJu+VqhwL7S7oNuAiYEr6vjpn1kJIn0oiIK4ErK2VfaHg+A3hbyZjMzEoabifSzMx6mpOumVlBTrpmZgU56ZqZFeSka2ZWkJOumVlBTrpmZgU56ZqZFeSka2ZWkJOumVlBtZKupKKXC5uZ9aq6R7qPSDopj/xlZmaLqW7SPQrYFrhT0h8k7SdphS7GZWbWk2ol3Yg4MyK2BTYBfgf8F+no92xJHhXMzKymAZ1Ii4i7IuIw0l0fjgI+APxW0t2SDmx2ax0zM1toQCfI8h0f/h3YF3gH6aj3e8DawH8Ck0l3hDAzsyZqJV1JW5AS7d7AfOA84OCIuLehzhWAbxRpZtaPuke6NwNXA1OBn+Rb71TNBi7uUFxmZj2pbtJdLyIe6K9CRDwDfHTwIZmZ9a66J77WlfT2aqGkt0vavsMxmZn1rLpJ9xvAqk3KV8rzzMyshrpJ903AbU3K78zzzMyshrpJ9zlgrSbl44AXOxeOmVlvq5t0rwKOl/RKE4Ok1YCv5HlmZlZD3d4LnwN+C8yWdHsu2xR4FNizG4GZmfWiWkk3Ih6R9GZgH2CzXPx94MKIeLZLsZmZ9ZzalwHn5HpmF2MxM+t5tZOupPHA9sAaVNqCI+LrHY7LzKwn1R17YR/gbGAB8BgQDbMDcNI1M6uh7pHuccDXgP+MiJe6GI+ZWU+r22VsTeAsJ1wzs8Gpm3SvBLbqZiBmZkuCus0LvyRdHLExcAdpTN1XRMTlnQ7MzKwX1U26p+e/RzWZF8CozoRjZtbb6l4c4XufmZl1gJOpmVlBtZKukoMkTZf0rKT1cvkRkvbobohmZr2j7pHup4HPA2cAaiifCxzc6aDMzHpV3aR7ILB/RJxMuiqtz63Axh2PysysR9VNuq8n3SWiaj6wbOfCMTPrbXWT7ixgiybluwIzOheOmVlvq9tP9yTgFEnLkdp0t5H0IeBwYN9uBWdm1mvq9tM9R9Jo4MvAcsD5wMPApyLiki7GZ2bWUwYyiPmZwJmSxgKviYhHuxeWmVlvqp10+0TEvG4EYma2JKg7iPkdLDpw+SIiYtOORWRm1sPqHuleVpleinSDyrcBp3YyIDOzXlb3RNqxzcolHUbqw2tmZjUMdsCby0m3ZTczsxoGm3S3B57tRCBmZkuCuifSplWLgLWAzYGmTQ9mZvZqdU+kPV6ZfhmYDhwVEVd3NiQzs95V90TaR7sdiJnZkqDonSMk7SzpHkkzJR3Ros4ekmbkAdMvLBmfmVm31W3T/Qv9XBzRKCLWa7GMUaQ+ve8E5gA3S5oWETMa6qwPHAm8LSKelLRGnXWamY0Uddt0zwU+C9wE/CGXbQNsCXwdeLrGMrYEZkbELABJFwO7s+jQkPsDp0bEkwAe38HMek3dpDsROD4ivtxYKOlIYOOI+GCNZYwDHmqYngNsVamzQV7u9aTbuh8TEb+oLkjSVGBqzdjNzIaNukn332k+iPkPSM0BnTIaWB+YDIwHfivpnyPiqcZKEXEG6X5tSKrV7GFmNhzUPZH2DCkRVk2m/sURc4EJDdPjc1mjOcC0iJgfEX8B7iUlYTOznlD3SPcbwKmSJgE35LKtgY8Ax9Rcxs3A+pImkpLtXsAHKnV+DOwNnJPH7d2AdKsgM7OeULef7gmSZpNuxb5HLr4L+EhEXFpzGQskHQxcRWqvPTsipks6DrglIqblee+SNAN4CTgsIqoXZpiZjViKGNlNopJiz9N/P6DXXHLAtoz07TazYU2tZtS+OELSMpLeL+k/JK2Sy94gabUOBGhmtkSoe3HEG4FrgBWAVUi9Fp4CPp6nP9aV6MzMekzdI91vAlcDawLPNZRPA3bocExmZj2rbu+FbYGtI+IlaZGmigeBtTselZlZjxrIgDdLNSlbB/h7h2IxM+t5dZPu1aSxF/qEpJVIA5j/rONRmZn1qLrNC4cCv5Z0D7AMcAnwRuBvLOy3a2ZmbdS9OGKupM1IV4ttQTpCPgO4ICKe6++1Zma2UNukK2kp0uhgO0bE2cDZXY/KzKxHtW3TjYj5wHxqDmJuZmat1T2R9m3gSEl124DNzKyJukl0O+DtwFxJd5KGenxFROzW6cDMzHpR3aQ7D/hhNwMxM1sS9Jt0Je0G/Ny3YDcz64x2bbo/Albtm5B0h6QJ/dQ3M7N+tEu61TEh16X55cBmZlbDQMZeMDOzQWqXdINF++dWp83MbADa9V4Q8BtJC/L0csDPJb3YWCkiNu1GcGZmvaZd0j22Mu1uY2Zmg9Bv0o2IatI1M7NB8Ik0M7OCnHTNzApy0jUzK8hJ18ysICddM7OCWvZekPThuguJiPM6E46ZWW/rr8vYqZXpMaRxF17O068h3VHiBcBJ18yshpbNCxGxYt8D2Au4nTSY+TL5sR3wZ+ADBeI0M+sJddt0TwI+FRHXR8SC/LgeOAT4WteiMzPrMXWT7rpUbtGTPQus07FozMx6XN2keyPwLUnj+gry828AN3QjMDOzXlQ36e4HrA7MljRb0mxgNrAGsH93QjMz6z21bkwZEfdL2hR4J/BPufgu4JqI8Pi6ZmY11b0bMDm5Xp0fZma2GGpfkSbpIEnTJT0rab1cdoSkPboXnplZb6mVdCUdAnweOINFb1Y5Fzi482GZmfWmuke6BwL7R8TJwIKG8luBjTselZlZj6qbdF8P3NmkfD6wbOfCMTPrbXWT7ixgiybluwIzOheOmVlvq9t74STgFEnLkdp0t5H0IeBwYN9uBWdm1mvq9tM9R9Jo4Muk27CfDzxMGo/hki7GZ2bWUwbST/dM4ExJY4HXRMSj3QvLzKw31e0y9gVJOwBExLy+hCtpeUlf6GaAZma9pO6JtGOAX0j6ZKV8BeCLHY3IzKyHDeQeafsBx0r6jqRR3QrIzKyXDSTpXg1sA7yDdNS7cndCMjPrXXWTbgBExD3A1rnsJmCDbgRlZtar6ibdV8ZbiIingJ2BXwK/6EJMZmY9q26XsWOBp/smIuIl4GBJfwS270ZgZma9qO7FEce2KD8HOKejEZmZ9bCWSVfSZ4HvRMTz+XkrERHf6HxoZma9p78j3U8C3weez89bCdINKtuStDNwMjAKOCsivtqi3vuAy4C3RsQtdZZtZjYStEy6ETGx2fPFlfv2nkq6z9oc4GZJ0yJiRqXeisCnSXcgNjPrKQPppztYWwIzI2JWRLwIXAzs3qTel4DjSUfYZmY9pV2bbi0R8fUa1cYBDzVMzwG2qqxzC2BCRPxM0mH9xDYVmFo3PjOz4aJdm24dAdRJuv2S9Jq8nCltVxhxBul+bUjyLeDNbMSo1abbIXOBCQ3T43NZnxWBTYDrJAG8DpgmaTefTDOzXlGyTfdmYH1JEyWNAfYCpvXNjIi/R8TYiFg3ItYFbgCccM2sp9QexFzSqsAuwDrAmMZ5EXFcu9dHxAJJBwNXkbqMnR0R0yUdB9wSEdP6X4KZ2chXK+lK2hr4GfAC8FpSs8BaeXo20DbpAkTElcCVlbKmg6BHxOQ6yzQzG0nqNi+cCFxA6oHwPGl4x3WAW0jdu8zMrIa6SXdT4JSICOAlYOmI+BvwH6S7SpiZWQ11k+6LDc//Brw+P38aWLujEZmZ9bC6J9JuBd4K3AtcB/yXpDWBDwK3dyc0M7PeU/dI92jg4fz888BjwLeBVfGVYWZmtdUdT/eWhuePkbqOmZnZAJW8OMLMbIlXt5/uqqReCjsAa1BJ1hGxRscjMzPrQXVPpJ0HbEwa1Pxv5LsDm5nZwNRNupOBt0fErV2Mxcys59Vt071/AHXNzKyFuon008BXJL0533bHzMwWQ93mhZnAsqSLJMjj3b4iIpyIzcxqqJt0LwJWBj6FT6SZmS22ukl3ErBlRNzZzWDMzHpd3TbdGcBK3QzEzGxJUDfpfh74uqSdJK0pabXGRzcDNDPrJXWbF/ru9nA1i7bnKk/7RJqZWQ11k+4OXY3CzGwJ0TbpSloKOAH4cETc0/2QzMx6V9s23YiYD0zE3cTMzAat7om07wP7dzMQM7MlQd023eWBfSS9E/gj8EzjzIj4VKcDMzPrRXWT7obkS4CB9Srz3OxgZlZT3dv1uPeCmVkH1D3SBUDSMsAbSUe390fE812JysysR9U6kSZpKUknAk8CtwF3AE9KOiF3KTMzsxrqHukeD+wNHAj8LpdtB3yFlLg/1/nQzMx6T92k+wFg34i4sqHsfkmPAWfhpGtmVkvdfrork27ZU3U/sErHojEz63F1k+5tpAHMqz4N/Llj0ZiZ9bi6zQuHA1dK2gm4IZdtDawN7NKNwMzMelGtI92I+C2wAXAZsEJ+/AB4U0T8rr/XmpnZQrX76UbEw8DRXYzFzKzn9Zt0694VIiKe6Ew4Zma9rd2R7jzaj60QNZZjZma0T5b9jbmwM6n3woLOhWNm1tv6TboR8ZtqmaTNgRNJV6SdDnypO6GZmfWeuv10kTRR0oXATcDjwEYR8amIeKxr0ZmZ9Zi2SVfS6pJOBu4GXgdsGxF7RkSzK9TMzKwf/SZdSUeTLvV9O7B7RLwjIm4uEpmZWQ9qdyLtS8BzwBzgIEkHNasUEbt1OjAzs17ULumeh2/HY2bWMe16L0wpFIeZ2RKhdu8FMzMbPCddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK6ho0pW0s6R7JM2UdEST+Z+VNEPS7ZJ+Jen1JeMzM+u2YklX0ijgVNLdgzcC9pa0UaXan4BJEbEp6SaYJ5SKz8yshJJHulsCMyNiVkS8CFwM7N5YISKujYhn8+QNwPiC8ZmZdV3JpDsOeKhhek4ua2U/4OddjcjMrLBheUNJSR8EJpHG8W02fyowtWhQZmYdUDLpzgUmNEyPz2WLkLQTcDTw9oh4odmCIuIM4Ixc30NPmtmIUbJ54WZg/XyvtTHAXsC0xgr5ppenA7tFxKMFYzMzK6JY0o2IBcDBwFXAXcClETFd0nGS+u48cSKwAvADSX+WNK3F4szMRqSibboRcSVwZaXsCw3PdyoZj5lZab4izcysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK8hJ18ysICddM7OCnHTNzApy0jUzK2jJTLqvGY2kAT3GTVhnqKM2sx4weqgDGBIvL2DP038/oJdccsC2XQrGzJYkS+aRrpnZEHHSNTMryEnXzKwgJ10zs4KcdM3MCnLSNTMryEm3rsXo2+v+vWZWtWT2010ci9G3F9y/18wW5SNdM7OCnHTNzApy0jUzK6ho0pW0s6R7JM2UdEST+UtLuiTPv1HSuiXj64pCg+uMm7DOsF3P6DHL9NxJyFL723pPsRNpkkYBpwLvBOYAN0uaFhEzGqrtBzwZEW+UtBdwPLBnqRi7otDgOg/PeWhYr2exTkJ+fHskDeg1o5Zampfmv9D11wAeNMkWS8neC1sCMyNiFoCki4HdgcakuztwTH5+GXCKJEVEFIxz6OWj4yXeYv7DKvUaK2fchHV4eM5DA3rN2uMnMPehB7sU0eJTqXwm6f3AzhHxsTz9IWCriDi4oc6duc6cPH1/rjOvsqypwNQ8uQlwZ4FNaGcsMK9trd6PAYZHHMMhBnAcwy0GKBPHvIjYudmMEdlPNyLOAM4AkHRLREwa4pCGRRzDIYbhEsdwiMFxDL8YhkMcJU+kzQUmNEyPz2VN60gaDawMPF4kOjOzAkom3ZuB9SVNlDQG2AuYVqkzDfhIfv5+4NdLXHuumfW0Ys0LEbFA0sHAVcAo4OyImC7pOOCWiJgGfA84X9JM4AlSYm7njK4FPTDDIY7hEAMMjziGQwzgOBoNhxhgiOModiLNzMx8RZqZWVFOumZmBY3opNvusuJBLnuCpGslzZA0XdKnc/lqkn4p6b78d9VcLknfyrHcLmmLhmV9JNe/T9JHWq2zn1hGSfqTpCvy9MR8mfTMfNn0mFze8jJqSUfm8nskvXsxYlhF0mWS7pZ0l6RthmhffCa/H3dKukjSMt3eH5LOlvRo7kfeV9axbZf0Fkl35Nd8S2p+ZUyLOE7M78ntkn4kaZV229jqe9NqP9aJo2HeoZJC0thu7o9WMUj6ZN4f0yWd0O19sVgiYkQ+SCfj7gfWA8YAtwEbdXD5awFb5OcrAvcCGwEnAEfk8iOA4/PzXYGfAwK2Bm7M5asBs/LfVfPzVQcYy2eBC4Er8vSlwF75+WnAx/Pzg4DT8vO9gEvy843y/lkamJj326gBxvB94GP5+RhgldL7AhgH/AVYtmE/TOn2/gC2B7YA7mwo69i2Azflusqv3WUAcbwLGJ2fH98QR9NtpJ/vTav9WCeOXD6BdKL8AWBsN/dHi32xA3ANsHSeXqPb+2KxckunFlT6AWwDXNUwfSRwZBfX9xPSuBH3AGvlsrWAe/Lz04G9G+rfk+fvDZzeUL5IvRrrHQ/8CngHcEX+IM5r+KK9sh/yB36b/Hx0rqfqvmmsVzOGlUnJTpXy0vtiHPBQ/qKOzvvj3SX2B7Bu5QvekW3P8+5uKF+kXrs4KvPeC1zQ7PvQt420+N7097mqGwfp0v03A7NZmHS7tj+avCeXAjs1qdfVfTHQx0huXuj7AvaZk8s6Lv8s3Ry4EVgzIh7Js/4KrNkmnsHG+U3gcODlPL068FRELGiyvFfWlef/PdcfbAwTgceAc5SaOc6StDyF90VEzAVOAh4EHiFt3x8pvz+gc9s+Lj8fTCx99iUdGS5OHP19rtqStDswNyJuq8wquT82ALbLzQK/kfTWxYxhUPuinZGcdIuQtALwQ+CQiPhH47xI/wa71udO0r8Bj0bEH7u1jppGk37KfTciNgeeIf2kfkW39wVAbjfdnfRPYG1geaDp9e0lldj2diQdDSwALhiCdS8HHAV8ofS6K0aTfgVtDRwGXNqqfXwojeSkW+ey4kGRtBQp4V4QEZfn4r9JWivPXwt4tE08g4nzbcBukmYDF5OaGE4GVlG6TLq6vFaXUQ92X80B5kTEjXn6MlISLrkvAHYC/hIRj0XEfOBy0j4qvT+gc9s+Nz9f7FgkTQH+Ddgn/wNYnDgep/V+bOcNpH+Et+XP6njgVkmvW4w4BrM/5gCXR3IT6dfh2MWIYTD7or1OtVOUfpD+q80ivdl9jeAbd3D5As4DvlkpP5FFT6CckJ//K4ueMLgpl69Gag9dNT/+Aqy2GPFMZuGJtB+waCP/Qfn5J1j0xNGl+fnGLHoiYRYDP5H2P8Cb8vNj8n4oui+ArYDpwHJ52d8HPllif/Dq9sOObTuvPnG06wDi2Jk0POprK/WabiP9fG9a7cc6cVTmzWZhm27X9keTfXEgcFx+vgGp6UDd3hcD/i53akFD8SCdGb2XdAby6A4v+19IPxlvB/6cH7uS2nt+BdxHOlPa90ERaZD2+4E7gEkNy9oXmJkfH13MeCazMOmulz+YM/OHo+9s7TJ5emaev17D64/Osd1Di7Pjbda/GXBL3h8/zl+U4vsCOBa4mzSc5/n5i9TV/QFcRGpDnk86mtqvk9sOTMrbcz9wCpUTlm3imElKLn2f0dPabSMtvjet9mOdOCrzZ7Mw6XZlf7TYF2OA/86vvRV4R7f3xeI8fBmwmVlBI7lN18xsxHHSNTMryEnXzKwgJ10zs4KcdM3MCnLStRFB0imSrhvqOKokrSrpb5LekKcnN46y1cH1fELSTzu5TBsaTrrWdZLOzYnoe03mHZ/nXdGF9V6Xl119XNxQp7H8WUmzJF0oabuaqzkKuDIi7u90/BVnAW8ZQFw2TDnpWikPAXvkgXKAVy7N/TBpAJtuOYc0clXj44BKnf1z+YakTvYvAr+RdFh/C85jDnyMdG+/roqIF0jDe36q2+uy7nLStVJuJ129tUdD2b8CzwPXNVZUGrT9JElP5sc3SZdtNta5TtJ3JH1Z0rw8oPVJkqqf6Wcj4q+Vx98rdZ7K5Q9ExLURMQX4KvAVSW/sZ5t2JV21eH2rCkoDqf9I0q2S1pC0bj6q3iuPhPVcHrltU0mbSPq9pGck/U7SxMrippHG4liun5hsmHPStZK+R7r0s8++pCPR6mWRh5KOPg8gjWU6CtinyfL2IY2stS1wMHAIsGeHYv0a6fvxnn7qbAf8MVpc1ilpJeAXpHEGJkfEow2zjyUNOr458BTpstZvky5X3ZJ0CfO3Kou8hTRewDYD2xQbTpx0raQLgUmS1s8jUO0MnNuk3iGkAWQujYi7gU+TxqytmhERX4iIeyPiUuBaYMdKnamSnq48DmoXaEQ8Tho5bL1+qr0eeLjFvDVyPP8LvDsqw4ICX4+IK/P2fY10d4Nv5yPt6aQxB3aoxPQsaUzgddvFb8PX6PZVzDojIp6U9CPSEe5TwHUR8WDjkKeSVia1r/6h4XUvS7qRRYfhg9Rk0ehhUrJrdAnpqLLRYzVDFv2Pk7ss8LcW864iDbryvlg4GHajxtj7lnFHpWx5ScvlZNvnubxeG6GcdK20s0lDMj7N4Ae9nl+ZDl796+3vETFzoAvOXb5eSxr6r5V5pNHWmrmC1H69CWn0r6rG2KOfsur2rEb9fxo2DLl5wUr7Fal3wFjSEJGLyCe5HiGNpwqkO8qS2jlLOpQ0CPaP+6nzJ1KzQDP/SRqH9VeSNutEQLkv8DKkI2gboXyka0VFREjalDRG6gstqp0MHCnpXtJP7oNITQ6PtKjfn+Vy+3GjFyPiiYbpVXKdMaS7IHyE1JXt8Db9b68Cjpe0em4DXkREHJ3/YVwjacd49f3DBmo7YFZE3DfI5dgQctK14iLif9tU+RrwOtIFAZAGK7+A1I92oD6aH42uJw1S3+fM/PcFUmK/gdTb4Lf9LTgi7pB0E+muFKe2qHNUTry/krQj6UTY4tq7IVYboTyIudkgSNqZdGS+UUS81MX1bEJqmtmgST9jG0Hcpms2CBHxC9JR7vh2dQdpbeDDTrgjn490zcwK8pGumVlBTrpmZgU56ZqZFeSka2ZWkJOumVlBTrpmZgX9f0q1KZpV6NgLAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_lgl_edinburgh, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('LGL - Edinburgh Geoparser', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "09769175",
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>wikipedia_name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>page</th>\n",
       "      <th>geovirus_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Pandi</td>\n",
       "      <td>Pandi, Bulacan</td>\n",
       "      <td>14.87</td>\n",
       "      <td>120.95</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Pandi,_Bulacan</td>\n",
       "      <td>Pandi, Bulacan 14.87,120.95</td>\n",
       "      <td>Pandi, Bulacan</td>\n",
       "      <td>14.87 , 120.95</td>\n",
       "      <td>14.866666666667 , 120.95</td>\n",
       "      <td>Pandi, Bulacan 14.87,120.95</td>\n",
       "      <td>[(14.865, 120.95722)]</td>\n",
       "      <td>0.797781</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>15.47</td>\n",
       "      <td>120.75</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Central_Luzon</td>\n",
       "      <td>Central Luzon 15.47,120.75</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>15.47 , 120.75</td>\n",
       "      <td>15.466666666667 , 120.75</td>\n",
       "      <td>Central Luzon 15.47,120.75</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>15.00</td>\n",
       "      <td>121.08</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Bulacan</td>\n",
       "      <td>Bulacan 15,121.08</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>15 , 121.08</td>\n",
       "      <td>15.0 , 121.08333333333</td>\n",
       "      <td>Bulacan 15,121.08</td>\n",
       "      <td>[(14.79945, 120.94436)]</td>\n",
       "      <td>26.838576</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Manaoag</td>\n",
       "      <td>Manaoag, Pangasinan</td>\n",
       "      <td>16.04</td>\n",
       "      <td>120.48</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Manaoag,_Pangasinan</td>\n",
       "      <td>Manaoag, Pangasinan 16.04,120.48</td>\n",
       "      <td>Manaoag, Pangasinan</td>\n",
       "      <td>16.04 , 120.48</td>\n",
       "      <td>16.043889 , 120.485556</td>\n",
       "      <td>Manaoag, Pangasinan 16.04,120.48</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>15.92</td>\n",
       "      <td>120.33</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Pangasinan</td>\n",
       "      <td>Pangasinan 15.92,120.33</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>15.92 , 120.33</td>\n",
       "      <td>15.916666666667 , 120.33333333333</td>\n",
       "      <td>Pangasinan 15.92,120.33</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0           name       wikipedia_name    lat     lon  \\\n",
       "0           0          Pandi       Pandi, Bulacan  14.87  120.95   \n",
       "1           1  Central Luzon        Central Luzon  15.47  120.75   \n",
       "2           2        Bulacan              Bulacan  15.00  121.08   \n",
       "3           3        Manaoag  Manaoag, Pangasinan  16.04  120.48   \n",
       "4           4     Pangasinan           Pangasinan  15.92  120.33   \n",
       "\n",
       "                                                page  \\\n",
       "0       https://en.wikipedia.org/wiki/Pandi,_Bulacan   \n",
       "1        https://en.wikipedia.org/wiki/Central_Luzon   \n",
       "2              https://en.wikipedia.org/wiki/Bulacan   \n",
       "3  https://en.wikipedia.org/wiki/Manaoag,_Pangasinan   \n",
       "4           https://en.wikipedia.org/wiki/Pangasinan   \n",
       "\n",
       "                       geovirus_key               toponym         old_coord  \\\n",
       "0       Pandi, Bulacan 14.87,120.95       Pandi, Bulacan    14.87 , 120.95    \n",
       "1        Central Luzon 15.47,120.75        Central Luzon    15.47 , 120.75    \n",
       "2                 Bulacan 15,121.08              Bulacan       15 , 121.08    \n",
       "3  Manaoag, Pangasinan 16.04,120.48  Manaoag, Pangasinan    16.04 , 120.48    \n",
       "4           Pangasinan 15.92,120.33           Pangasinan    15.92 , 120.33    \n",
       "\n",
       "                            new_coord                           mlg_key  \\\n",
       "0            14.866666666667 , 120.95       Pandi, Bulacan 14.87,120.95   \n",
       "1            15.466666666667 , 120.75        Central Luzon 15.47,120.75   \n",
       "2              15.0 , 121.08333333333                 Bulacan 15,121.08   \n",
       "3              16.043889 , 120.485556  Manaoag, Pangasinan 16.04,120.48   \n",
       "4   15.916666666667 , 120.33333333333           Pangasinan 15.92,120.33   \n",
       "\n",
       "  geocoded_coordinates_list  median_error_distance  \n",
       "0     [(14.865, 120.95722)]               0.797781  \n",
       "1                        []                    NaN  \n",
       "2   [(14.79945, 120.94436)]              26.838576  \n",
       "3                        []                    NaN  \n",
       "4                        []                    NaN  "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_geovirus_edinburgh = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/geovirus_geocoded_results_Edinburgh_Geoparser.csv')\n",
    "df_geovirus_edinburgh.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "e77b05f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAF3CAYAAACSWTEVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAApXElEQVR4nO3debgcVZ3/8feHhBD2LbIlQRIFhmUQMLIpEAQUcR746fgTEJUdFVBQB4ZNB9BR2VQUHNk32RE1IgqCRkZkR1kSCIYQJQHZQdkT/M4f5zTpVLrvrc7t7lu5/Xk9Tz+3+9Tpqm+f7v7e6lOnTikiMDOzwbfYYAdgZmaJE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCFXlKSJkkLS3oMdS1U0apN2tJOkmZImtyHEjpF0XH6daw12LNY5Qz4hSxop6SBJv5H0tKQ5kl6QdKekEyX9SxdiuCp/mTbuo44kPZpjW7LTMXVbTnrRx+2Tgx1jr6jCd8IaGz7YAXSSpPHAtcB6wO+A7wBPAMsAGwP7Av8hac2ImN3BUM4FPgbsAxzapM52wFrAmRHxqqSbgSWBOR2Mq9tmAUc1WXbLQq5zKLZTx1ToO2ENDNmEnPcyfwG8A/hoRPykQZ2RwBeBTp+ueAPwGLCnpMMj4o0GdfbJf88FiIh/Aq/1t2JJApaOiJfaFWwHvRgRP2rnCsu2U5VIWjYi/jEI263Sd6Lrut3uC7W9iBiSN+DzpA/Vfy/Ec5cHTgSmA68DTwOXAeMb1B0FnEFKuG/kv2cAKxfqnZDj+ViDdSwHvALcX1c2Mdffu1EZcDAwNcd3XF4+E5jcYP2N1jUSOA6Ylrf9AnA/cHKH3o+ZwAMt1N8V+CMp2T4GfA3Ysa82adJO+wBTcjv9BTiiSWyTgU2B3wAvAc8BFwKrFOoel9e9VrP1FMoCuADYHvh9XvfkuuWfy+/B68CfgUNy3AFMbLDddYFvkH5tvA7cC+y8iH4naq9zh/z6/pLXfR+we4P1fgC4ApgBvJo/szcA2zaoOzm/H+OBq/P7GXnZYsBheTv/AP6e34NzgcUL65kA/AR4Jsc2DTgGGF52e63chuweMqmLAOCcVp4kaXngD8CawHmkL/PqwEHA7ZImRMRfCnXfmeveA2xC+pK9X9JmMe8/5PnAsaQEcXVhs7uTfnafWzLMw4CVgbOBv5E+8K06g/Tz9CLg26RfS2sD71+IdZU1TNKoJsuerX1bJH0E+DHpA34CMJfUbh9ucXufBVYltesLwCeBEyXNiohLC3XHADfl7V5NSs77AhMkvSciXmlx2/UmAP9Oer8urBVK+k/gW6TPzVHAUsDhpGTXzIWk7plTgBGkz8JPJa0TETP7iaNq34maE4GlgR/kx/sAl0kaGREX1NXbG1iJ9JmdBYwG9gdukrRdRPxvYb3LkLplbiEl0VVy+TGkz9XPgR8CbwLjgF2AJcjdX5I+DFxD+id0KinJbpmfuzHw/0tur7xO7A1V4QY8S/qJXCwfRvoPXn9bsm75aaT/vu8qPO/tpP+kF9SV/TfpP/xBhboH5/KvFcpvIiWX1Qvlt5L++46qK5tI8z2/5yjsueXlMym/h/wccF0X34+ZOYZmt1F1789fSXsk9e2xPGkPqpU95MeB5evKlyIlu1ubxHZYobz20/3IurLjaH0POYAdCuUr5c/ZfcDIuvLVgBdpvod8LaC68vfk8m8uat8J5u0h/6XwPtXe6+cKcSzdIPZV82flukL55Lzurzd4zj3A1H7aaiRpZ+dmFtwbrn0uJpbZXiu3oTzKYjnSh6VoPdKXsv52MLzVH7sn6U2YLWlU7Qa8DNxG+tlU85H8/LMK2zgzl3+kUH4u6cP/6VpBPqK9BTApIp4p+douioinStZt5kVgA0kbDnA9rZhJ6nZodHsx13k3MBY4v749IuJF0t5MK87Pz6ut4xXSe7h2g7p/Z94eWs0PcnnxfWzVvRFxY6FsR9KX/n8i4q0+8Ij4G3BJH+s6LXIGyPXvJHWDNHpNRVX8TkBqg/r3qfZer0j651orf7l2X9IyklYm7d3eDmze5DWf0qDsRWC0pPc1eQ6k92dV0i/bFQqv+7pc5wMNntdoe6UN5S6Lv5M+gEWPkhob4F3M34BvI3UFfIDmPxv/WXd/HHBXRMytrxARcyU9TPrZW+8a0k/nfUg/0yD9LIb0866sh1uo28xhwMXA/ZJmAL8l/YT7eaQDZU1JWq1Q9EZEPFdimy83SExF4/Pfhxosm1piG/VmNCh7lvQeL1A3CgdbI+L13DbjG9RvRaP3a1z+O63BskZlNa28pqIqficAHmxQVnuv32p7Se8g7YF/EFihUD9Y0NMR8UKD8qOBnwL/K+lx0t7tL4Cr6z4D6+W/fX0vVy25vdKGckJ+ANhG0riIeLRWmP/L3gggaW7hOcp/b2RewmybiHhN0qXAQZK2Iv1n/xSpP+z6FlbVrD+z0YcSGrzPEfGzfJLBzsC2pAMr+5E+pDsUk1PBE4XHv6NuT6ZC3uzAOpu1MTT/Pg2k/7mo2WtSk/J6lftOlCVpGdJe+tLAd0kHoP9B+mdwFI2PfTRs94i4NSf3D5KGm24HfAI4VtL78s5F7XUfDvypSViPl9leK4ZyQr4a2IbU6X9Myec8TdqDXa7EnhykvZV1JQ2v3yOQNBxYh8Z7M+eSDobsQ+pHXI101LvPvdKSnsvrLGq4h5c/eD8CfpR/mn4LOII0wuGqPrazY+Hx862H2lStzRqdnLB+G7dTNF7SiPp/RJKWILVd/d567ZfASqQumFrdkaQDXdNLbq/23HVJIzvqrVs66tZU9TuxHvCzQlntva7V3x5YA9g3Is6vryjp6yXimk+kYaI/zjckHUQ60L0fcDJpxAuU+1XXNkO5D/kc0hfp8HzUvpH59ipyUrwE2EzSxxo+Qao/cvpT0k+6/QvVDsjlPyk+PyLuIf3H3Y15Bzpa6a7oy8PAv0gaXRfvEnk71JUNk7RCIa4gDTODxkm9vu6NhdvdbYk+uZv0i2Gf+hEZkpYjjZrolOVI/yjrHZTLf1pXVut+2KFQ94u09n36NelA7udyMgfe6g7as4X1tKKS3wlSGyxft77lSe/1C6RfXzDvl8F88Un6AM37jxtqMtLnnvy39tm/HngKOFLSAt8HSUtKWraV7ZYxZPeQI53t9mHSUelrlOYquIF05HQ50h7YbqQ3un7Y2DHAe4ErJV1JOmjxBumI8s6khLF3rnsSaejLGZI2JSW0TUj/Zafl5Y2cC3wf2Il0VL7RXsPCOJ00hO5GST8kDYv6FAv+lFoWeELSpBzzU6S+v8+R9nZ/3qZ4ipbv4xTp+yPi3oh4U9IXgSuBOySdTRqZsi+pr3TNDsX2CPBf+SDn3aSDi/uSEtj36urdSHpvT8gHlR4F3kc6MFv2oCwR8ayk40ljim+R9CPSKJADSUl/Am0+OaPC34lnSMPnanu++5De5/1j3nDD3+c4T81dbbNIQ88+Req++NcWmuJBSbeRugwfJ/2yOTC/pstzW70s6dOkfzDTJJ1H+vWzQm6nj5IOUE5uYbv9G8gQjUXhRhrfezDpoNUzpDGGLwB3kT4c6zZ4zlLAV0hv9KukvqoHSeNINy/UfRvpaPysvO5ZpJ8+o/qIacW83gA+1aTOREoM8WrwvL1IH/w3SMniCFL/2lvPIyXqbwJ3kJLc66Sf0OcBa3fofZhJ38Pevl6o/1HSL4nXGcCJIQ3iuIDCgH0WPDHkZdI/pouBVRusYx3gV8w7oeZK0pjYmTQ5MaSPdjmYlIDrTwypncCxWV2942hhuN2i8J1g/hNDjicNd3w9b+MTDWLYKLf783n7k4Gtm7ynk4GZTV7/kaT+6KfqPl9XAZs2qLshqVtvNuk79SRpnPVXgJXKbK+Vm/LKzKwiJH2flJhXjzQMbkhSmqHvfGC7iJg8uNFUw1DuQzartPq+47qy1Unj1B8YysnYGutaQpZ0nqSnJD3QZLkkfU/SdEn35f4ns6FsoqT7JR0v6QBJ/02am2IZ0s9q6zHdPKh3Aemg00VNln+IdLbR2qSjpv9Di0dPzRYx00kHEw8gnXzxGqkf95vRxaFWVh1d7UPOR0evjYgFTteVdCbpwMRl+fE00rnixZMQzMyGpCoNexvN/ENtarM5LZCQJR1IGqbC+uuv/+4pU6Z0JUAzszZpeGblInlQLyLOiogJETFhySWH3NWOzKxHVSkhzybN8lUzJpeZmfWEKiXkScCn82iLLUjztrr/2Mx6Rtf6kCVdRjqDapSkWcB/AYsDRMQPSXOM7kw68vwK864xZ2bWE7qWkCNij36WB4VJcMzMekmVuizMzHqaE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRSgiBjuGARkxYkTMmTNnoZ+/xpixzH7sr22MyMysX2pYuKgnZEmx25l/WOjnX/GZrVjU28DMFjkNE7K7LMzMKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpyQzcwqwgnZzKwiupqQJe0kaZqk6ZKObLB8TUm/lfRHSfdJ2rmb8ZmZDaauJWRJw4AzgA8B6wN7SFq/UO1Y4MqI2ATYHfhBt+IzMxts3dxD3gyYHhEzIuIN4HJg10KdAJbL95cHHu9ifGZmg2p4F7c1Gnis7vEsYPNCneOAGyR9Hlga2KE7oZmZDb6qHdTbA7ggIsYAOwMXS1ogRkkHSrpL0l1dj9DMrEO6mZBnA2PrHo/JZfX2A64EiIhbgZHAqOKKIuKsiJgQERM6FKuZWdd1MyHfCawtaZykEaSDdpMKdf4KbA8gaT1SQn66izGamQ2ariXkiJgLHAJcDzxIGk0xRdIJknbJ1b4MHCDpXuAyYO+IiG7FaGY2mLp5UI+IuA64rlD21br7U4H3djMmM7OqqNpBPTOznuWEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRpRKypK6eYm1m1ovK7iE/IemUPAObmZl1QNmEfDSwFfCApFsl7SdpmQ7GZWbWc0ol5Ig4OyK2AjYEfg98nbTXfJ4kz85mZtYGLR3Ui4gHI+Jw0tU+jgY+Adws6SFJn210uSUzMyunpYN1+UofHwX2Bd5P2ls+F1gD+AowkXQlEDMza1GphCxpU1IS3gOYA1wEHBIRD9fVuRbwRUfNzBZS2T3kO4EbgAOBn+XLMRXNBC5vU1xmZj2nbEIeHxF/6atCRLwM7DPwkMzMelPZg3BrSdq2WChpW0nbtDkmM7OeVDYhfwdYsUH5cnmZmZkNUNmEvC5wb4PyB/IyMzMboLIJ+VVg9Qblo4E32heOmVnvKpuQrwdOlPRWt4WklYBv5mVmZjZAZUdZ/AdwMzBT0n25bCPgKWC3TgRmZtZrSiXkiHhC0ruAPYGNc/GFwKUR8UqHYjMz6ymlT53OiffsDsZiZtbTSidkSWOAbYBVKPQ9R8S32xyXmVnPKTuXxZ7AecBc4Gkg6hYH4IRsZjZAZfeQTwBOBb4SEW92MB4zs55VdtjbqsA5TsZmZp1TNiFfB2zeyUDMzHpd2S6LX5NODNkAuJ80J/JbIuKadgdmZtZryibkM/PfoxssC2BYe8IxM+tdZU8M8bXyzMw6zInWzKwiSiVkJQdJmiLpFUnjc/mRkj7e2RDNzHpD2T3kQ4FjgbMA1ZXPBg5pd1BmZr2obEL+LHBARJxGOluv5h5gg7ZHZWbWg8om5LeTrg5SNAdYsn3hmJn1rrIJeQawaYPynYGp7QvHzKx3lR2HfApwuqSlSH3IW0r6FHAEsG+ngjMz6yVlxyGfL2k48A1gKeBi4HHgCxFxRQfjMzPrGa1MUH82cLakUcBiEfFU58IyM+s9pRNyTUQ804lAzMx6XdkJ6u9n/knp5xMRG7UtIjOzHlV2D/nqwuPFSRc7fS9wRjsDMjPrVWUP6h3fqFzS4aQxymZmNkADnVzoGmDPdgRiZtbrBpqQtwFeaUcgZma9ruxBvUnFImB1YBOgYXeGmZm1puxBvWcLj/8JTAGOjogb2huSmVlvKntQb59OB2Jm1uu6esUQSTtJmiZpuqQjm9T5uKSpeTL8S7sZn5nZYCrbh/wofZwYUi8ixjdZxzDSmOUdgVnAnZImRcTUujprA0cB742I5yWtUmabZmZDQdk+5AuALwF3ALfmsi2BzYBvAy+VWMdmwPSImAEg6XJgV+afvvMA4IyIeB7A82WYWS8pm5DHASdGxDfqCyUdBWwQEZ8ssY7RwGN1j2cBmxfqrJPXewswDDguIn5VXJGkA4EDS8ZuZrZIKJuQP0rjCeqvInUxtMtwYG1gIjAGuFnSv0bEC/WVIuIs0vX9kFSqK8XMrOrKHtR7mZQkiyZS/sSQ2cDYusdjclm9WcCkiJgTEY8CD5MStJnZkFd2D/k7wBmSJgC35bItgL2A40qu405gbUnjSIl4d+AThTo/BfYAzs/zLq9DunyUmdmQV3Yc8kmSZgKHAh/PxQ8Ce0XElSXXMVfSIcD1pP7h8yJiiqQTgLsiYlJe9gFJU4E3gcMjonhSipnZkKSIRbsLVlLsduYfFvr5V3xmKxb1NjCzRY4aFZY+MUTSSEkfk/SfklbIZe+QtFKbAjQz62llTwx5J3AjsAywAml0xQvA5/Lj/TsSnZlZDym7h/xd4AZgVeDVuvJJwHZtjsnMrCeVHWWxFbBFRLwpzdf18VdgjbZHZWbWg1qZXGjxBmVrAi+2KRYzs55WNiHfQJrLoiYkLUeanP4XbY/KzKwHle2y+DLwG0nTgJHAFcA7gSeZNy7ZzMwGoOyJIbMlbUw6i25T0p71WcAlEfFqX881M7Ny+k3IkhYnzdK2fUScB5zX8ajMzHpQv33IETEHmEPJCerNzGzhlD2o933gKEll+5zNzKxFZRPs1sC2wGxJD5Cm43xLROzS7sDMzHpN2YT8DPDjTgZiZtbr+kzIknYBfhkR+3QpHjOzntVfH/JPgBVrDyTdL2lsH/XNzGwh9ZeQi3N2rkXjU6jNzGyAWpnLwszMOqi/hBzMP/64+NjMzNqkv1EWAn4naW5+vBTwS0lv1FeKiI06EZyZWS/pLyEfX3jsoW9mZh3SZ0KOiGJCNjOzDvFBPTOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4poOspC0qfLriQiLmpPOGZmvauvYW9nFB6PIM1j8c/8eDHSlUReB5yQzcwGqGmXRUQsW7sBuwP3kSaqH5lvWwN/Aj7RhTjNzIa8sn3IpwBfiIhbImJuvt0CHAac2rHozMx6SNmEvBaFyzZlrwBrti0aM7MeVjYh3w58T9LoWkG+/x3gtk4EZmbWa8om5P2AlYGZkmZKmgnMBFYBDuhMaGZmvaXURU4j4hFJGwE7Av+Six8EbowIz49sZtYGZa86TU68N+SbmZm1Wekz9SQdJGmKpFckjc9lR0r6eOfCMzPrHaUSsqTDgGOBs5j/wqezgUPaH5aZWe8pu4f8WeCAiDgNmFtXfg+wQdujMjPrQWUT8tuBBxqUzwGWbF84Zma9q2xCngFs2qB8Z2Bq+8IxM+tdZUdZnAKcLmkpUh/ylpI+BRwB7Nup4MzMeknZccjnSxoOfANYCrgYeJw0v8UVHYzPzKxntDIO+WzgbEmjgMUi4qnOhWVm1nvKDnv7qqTtACLimVoylrS0pK92MkAzs15R9qDeccCvJH2+UL4M8F9tjcjMrEe1ck29/YDjJf1A0rBOBWRm1qtaScg3AFsC7yftLS/fmZDMzHpT2YQcABExDdgil90BrNOJoMzMelHZhPzW/BUR8QKwE/Br4FcdiMnMrCeVHfZ2PPBS7UFEvAkcIuluYJtOBGZm1mvKnhhyfJPy84Hz2xqRmVmPapqQJX0J+EFEvJbvNxMR8Z32h2Zm1lv62kP+PHAh8Fq+30yQLnbaL0k7AacBw4BzIuJbTer9O3A18J6IuKvMus3MFnVNE3JEjGt0f2HlsctnkK7LNwu4U9KkiJhaqLcscCjpStdmZj2jlXHIA7UZMD0iZkTEG8DlwK4N6n0NOJG0Z25m1jP660MuJSK+XaLaaOCxusezgM0L29wUGBsRv5B0eB+xHQgcWDY+M7NFQX99yGUEUCYh90nSYnk9e/e7wYizSNf3Q1IMdNtmZlVQqg+5TWYDY+sej8llNcsCGwKTJQGsBkyStIsP7JlZL+hmH/KdwNqSxkkaAewOTKotjIgXI2JURKwVEWsBtwFOxmbWM0pPUC9pReBDwJrAiPplEXFCf8+PiLmSDgGuJw17Oy8ipkg6AbgrIib1vQYzs6GtVEKWtAXwC+B14G2krobV8+OZQL8JGSAirgOuK5Q1nOA+IiaWWaeZ2VBRtsviZOAS0kiJ10hTcK4J3EUaomZmZgNUNiFvBJweEQG8CSwREU8C/0m6moiZmQ1Q2YT8Rt39J4G35/svAWu0NSIzsx5V9qDePcB7gIeBycDXJa0KfBK4rzOhmZn1lrJ7yMcAj+f7xwJPA98HVsRnzJmZtUXZ+ZDvqrv/NGn4m5mZtVE3TwwxM7M+lB2HvCJpNMV2wCoUEnlErNL2yMzMekzZg3oXARuQJqx/knwVajMza5+yCXkisG1E3NPBWMzMelrZPuRHWqhrZmYLoWySPRT4pqR35UsxmZlZm5XtspgOLEk6QYQ8X/FbIsJJ2sxsgMom5MuA5YEv4IN6ZmYdUTYhTwA2i4gHOhmMmVkvK9uHPBVYrpOBmJn1urIJ+Vjg25J2kLSqpJXqb50M0MysV5Ttsqhd5eMG5u8/Vn7sg3pmZgNUNiFv19EozMys/4QsaXHgJODTETGt8yGZmfWmfvuQI2IOMA4PdTMz66iyB/UuBA7oZCBmZr2ubB/y0sCeknYE7gZerl8YEV9od2BmZr2mbEJej3zaNDC+sMxdGWZmbVD2Ek4eZWFm1mFl95ABkDQSeCdpr/iRiHitI1GZmfWgUgf1JC0u6WTgeeBe4H7geUkn5WFxZmY2QGX3kE8E9gA+C/w+l20NfJOU1P+j/aGZmfWWsgn5E8C+EXFdXdkjkp4GzsEJ2cxswMqOQ16edBmnokeAFdoWjZlZDyubkO8lTU5fdCjwp7ZFY2bWw8p2WRwBXCdpB+C2XLYFsAbwoU4EZmbWa0rtIUfEzcA6wNXAMvl2FbBuRPy+r+eamVk5pcchR8TjwDEdjMXMrKf1mZDLXg0kIp5rTzhmZr2rvz3kZ+h/rooosR4zM+tHf4m0rzksdiKNspjbvnDMzHpXnwk5In5XLJO0CXAy6Uy9M4GvdSY0M7PeUnYcMpLGSboUuAN4Flg/Ir4QEU93LDozsx7Sb0KWtLKk04CHgNWArSJit4hodOaemZktpD4TsqRjSKdHbwvsGhHvj4g7uxKZmVmP6e+g3teAV4FZwEGSDmpUKSJ2aXdgZma9pr+EfBG+RJOZWVf0N8pi7y7FYWbW80qPsjAzs85yQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzq4iuJmRJO0maJmm6pCMbLP+SpKmS7pN0k6S3dzM+M7PB1LWELGkYcAbpKtXrA3tIWr9Q7Y/AhIjYiHRB1ZO6FZ+Z2WDr5h7yZsD0iJgREW8AlwO71leIiN9GxCv54W3AmC7GZ2Y2qLqZkEcDj9U9npXLmtkP+GVHIzIzq5BKXpxU0ieBCaR5mBstPxA4sKtBmZl1WDcT8mxgbN3jMblsPpJ2AI4Bto2I1xutKCLOAs7K9T09qJkNCd3ssrgTWDtfm28EsDswqb5CvoDqmcAuEfFUF2MzMxt0XUvIETEXOAS4HngQuDIipkg6QVLtiiMnA8sAV0n6k6RJTVZnZjbkdLUPOSKuA64rlH217v4O3YzHzKxKfKaemVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBPyYsORNKDb6LFrDvarMLMhYPhgBzDo/jmX3c78w4BWccVntmpTMGbWy7yHbGZWEU7IZmYV4YTcDgPshx4+YqT7sc2su33IknYCTgOGAedExLcKy5cALgLeDTwL7BYRM7sZ40IZYD/0FZ/Zyv3YZta9PWRJw4AzgA8B6wN7SFq/UG0/4PmIeCfwHeDEbsXX60aPXXPQ99AHGsNAf2n4V4YNtm7uIW8GTI+IGQCSLgd2BabW1dkVOC7fvxo4XZIiIroYZ096fNZjA9vL/9w2SBpwHIP5S6Mdr2GNMWOZ/dhfB7SORd3osWvy+KzHBrSOwW7HwXoN6lauk/QxYKeI2D8//hSweUQcUlfngVxnVn78SK7zTGFdBwIH5ocbAg904SWUMQp4pt9anec4FlSVWKoSB1QnlqrEAd2L5ZmI2KlYuEiOQ46Is4CzACTdFRETBjkkoDqxOI4FVSWWqsQB1YmlKnHA4MfSzVEWs4GxdY/H5LKGdSQNB5YnHdwzMxvyupmQ7wTWljRO0ghgd2BSoc4kYK98/2PAb9x/bGa9omtdFhExV9IhwPWkYW/nRcQUSScAd0XEJOBc4GJJ04HnSEm7P2d1LOjWVSUWx7GgqsRSlTigOrFUJQ4Y5Fi6dlDPzMz65jP1zMwqwgnZzKwiFumELGknSdMkTZd0ZAfWP1bSbyVNlTRF0qG5fCVJv5b05/x3xVwuSd/L8dwnadO6de2V6/9Z0l7NttlPPMMk/VHStfnxOEm35+1dkQ+WImmJ/Hh6Xr5W3TqOyuXTJH1wIeNYQdLVkh6S9KCkLQejTSR9Mb8vD0i6TNLIbrWJpPMkPZXHztfK2tYGkt4t6f78nO9Jjc9YaRLHyfm9uU/STySt0N9rbfZdataeZWOpW/ZlSSFp1GC0SS7/fG6XKZJO6kabtCwiFskb6cDgI8B4YARwL7B+m7exOrBpvr8s8DDptO+TgCNz+ZHAifn+zsAvAQFbALfn8pWAGfnvivn+igsRz5eAS4Fr8+Mrgd3z/R8Cn8v3DwJ+mO/vDlyR76+f22kJYFxuv2ELEceFwP75/ghghW63CTAaeBRYsq4t9u5WmwDbAJsCD9SVta0NgDtyXeXnfqiFOD4ADM/3T6yLo+FrpY/vUrP2LBtLLh9LOpj/F2DUILXJdsCNwBL58SrdaJOWv1vtWlG3b8CWwPV1j48CjurwNn8G7AhMA1bPZasD0/L9M4E96upPy8v3AM6sK5+vXsltjwFuAt4PXJs/lM/UffHeao/84d8y3x+e66nYRvX1WohjeVIiVKG8q21CSsiP5S/u8NwmH+xmmwBrFb70bWmDvOyhuvL56vUXR2HZR4BLGn1Haq+VJt+lvj5jrcRCmgbhXcBM5iXkrrYJKYnu0KBex9ukldui3GVR+0LWzMplHZF/4m4C3A6sGhFP5EV/A1btJ6Z2xPpd4Ajgn/nxysALETG3wTrf2l5e/mKu3444xgFPA+crdZ+cI2lputwmETEbOAX4K/AE6TXezeC0SU272mB0vt+OmPYl7U0uTBx9fcZKkbQrMDsi7i0s6nabrANsnbsafifpPQsZx4DbpC+LckLuGknLAD8GDouIv9cvi/RvsqNjByX9G/BURNzdye2UNJz0c/B/ImIT4GXSz/O3dKlNViRNRjUOWANYGlhgboDB0o026I+kY4C5wCWDtP2lgKOBrw7G9guGk35NbQEcDlzZrA96MC3KCbnMqdgDJmlxUjK+JCKuycVPSlo9L18deKqfmAYa63uBXSTNBC4ndVucBqygdIp5cZ3NTkFvR5vNAmZFxO358dWkBN3tNtkBeDQino6IOcA1pHYajDapaVcbzM73FzomSXsD/wbsmf85LEwcz9K8Pct4B+kf5r35szsGuEfSagsRy0DbZBZwTSR3kH5pjlqIOAbaJn1rV99Ht2+k/3gzSG94rdN9gzZvQ6QJ879bKD+Z+Q/enJTvf5j5D1TckctXIvW7rphvjwIrLWRME5l3UO8q5j+4cFC+fzDzH8C6Mt/fgPkPYMxg4Q7q/S+wbr5/XG6PrrYJsDkwBVgqr/tC4PPdbBMW7KdsWxuw4AGsnVuIYyfStLZvK9Rr+Frp47vUrD3LxlJYNpN5fcjdbpPPAifk++uQuiPUjTZp6bvVrhUNxo10pPZh0tHQYzqw/veRfnbeB/wp33Ym9SPdBPyZdOS29oERaRL+R4D7gQl169oXmJ5v+wwgponMS8jj84d0ev6Q1I4gj8yPp+fl4+uef0yObxpNjlKXiGFj4K7cLj/NX5yutwlwPPAQafrVi/OXqittAlxG6rueQ9r72q+dbQBMyK/rEeB0CgdR+4ljOinh1D6zP+zvtdLku9SsPcvGUlg+k3kJudttMgL4UX7+PcD7u9Emrd586rSZWUUsyn3IZmZDihOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsi2SJJ0uafJgx1EkaUVJT0p6R348sX6WszZu52BJP2/nOm3wOSFb10m6ICepcxssOzEvu7YD252c1128XV5Xp778FUkzJF0qaeuSmzkauC4iHml3/AXnAO9uIS5bBDgh22B5DPh4npgIeOuU5k+TJgzqlPNJM4fV3z5TqHNALl+PdFLBG8DvJB3e14rz3A37k64N2VER8TppKtYvdHpb1j1OyDZY7iOd0fbxurIPA68Bk+srKk3Mf4qk5/Ptu6TTW+vrTJb0A0nfkPRMnqD8FEnFz/grEfG3wu3FQp0XcvlfIuK3EbE38C3gm5Le2cdr2pl0ZuctzSooTZT/E0n3SFpF0lp5b3z3PAvZq3kWvY0kbSjpD5JelvR7SeMKq5tEmuNkqT5iskWIE7INpnNJp8nW7Evagy2ePvpl0l7rZ0jzzw4D9mywvj1Js5ttBRwCHAbs1qZYTyV9X/5fH3W2Bu6OJqe/SloO+BVpvoaJEfFU3eLjSZPJbwK8QDr99/uk03o3I536/b3CKu8izbmwZWsvxarKCdkG06XABElr5xnAdgIuaFDvMNJEPVdGxEPAoaT5houmRsRXI+LhiLgS+C2wfaHOgZJeKtwO6i/QiHiWNHvb+D6qvR14vMmyVXI8/wA+GIVpXIFvR8R1+fWdSrqSxffzHvoU0twN2xVieoU0r/Na/cVvi4bh/Vcx64yIeF7ST0h7xi8AkyPir/XT1EpantSfe2vd8/4p6Xbmnx4RUjdIvcdJibDeFaS90XpPlwxZ9D3H8ZLAk02WXU+a1ObfY97k5vXqY6+t4/5C2dKSlsqJuObVvF0bApyQbbCdR5o68yUGPpH5nMLjYMFfgS9GxPRWV5yHrb2NNCVjM8+QZr5r5FpSf/mGpBnYiupjjz7Kiq9nJcr/Q7GKc5eFDbabSKMYRpGm8pxPPuD2BGkeXCBdsZjUr9pNXyZNav7TPur8kdTV0MhXSHPn3iRp43YElMc6jyTtedsQ4D1kG1QREZI2Is1t+3qTaqcBR0l6mPQz/iBSN8YTTer3ZancX13vjYh4ru7xCrnOCNJVL/YiDcc7op/xxdcDJ0paOfc5zycijsn/TG6UtH0seJ25Vm0NzIiIPw9wPVYRTsg26CLiH/1UORVYjXQyBKTJ6C8hjRNu1T75Vu8W0sUIas7Of18nJf3bSKMibu5rxRFxv6Q7SFckOaNJnaNzUr5J0vakg3ILa4+6WG0I8AT1Zm0kaSfSHv36EfFmB7ezIam7Z50G46htEeU+ZLM2iohfkfaOx/RXd4DWAD7tZDy0eA/ZzKwivIdsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV8X8EzCAHHzp5JQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_geovirus_edinburgh, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('GeoVirus - Edinburgh Geoparser', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "a1dc7924",
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>wikipedia_name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>page</th>\n",
       "      <th>wiktor_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, California</td>\n",
       "      <td>34.9514</td>\n",
       "      <td>-120.4333</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Ca...</td>\n",
       "      <td>Santa Maria, California 34.9514,-120.4333</td>\n",
       "      <td>Santa Maria, California</td>\n",
       "      <td>34.9514 , -120.4333</td>\n",
       "      <td>34.951388888889 , -120.43333333333</td>\n",
       "      <td>Santa Maria, California 34.9514,-120.4333</td>\n",
       "      <td>[(34.95303, -120.43572), (34.95303, -120.43572...</td>\n",
       "      <td>0.283927</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Bulacan</td>\n",
       "      <td>14.8208</td>\n",
       "      <td>120.9636</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Bu...</td>\n",
       "      <td>Santa Maria, Bulacan 14.8208,120.9636</td>\n",
       "      <td>Santa Maria, Bulacan</td>\n",
       "      <td>14.8208 , 120.9636</td>\n",
       "      <td>14.818333333333 , 120.95633333333</td>\n",
       "      <td>Santa Maria, Bulacan 14.8208,120.9636</td>\n",
       "      <td>[(14.83333, 121.0), (14.83333, 121.0), (14.833...</td>\n",
       "      <td>4.981280</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Ilocos Sur</td>\n",
       "      <td>17.3585</td>\n",
       "      <td>120.4981</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Il...</td>\n",
       "      <td>Santa Maria, Ilocos Sur 17.3585,120.4981</td>\n",
       "      <td>Santa Maria, Ilocos Sur</td>\n",
       "      <td>17.3585 , 120.4981</td>\n",
       "      <td>17.366666666667 , 120.48333333333</td>\n",
       "      <td>Santa Maria, Ilocos Sur 17.3585,120.4981</td>\n",
       "      <td>[(16.9518, 120.4445), (16.9518, 120.4445), (16...</td>\n",
       "      <td>46.315270</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Romblon</td>\n",
       "      <td>12.4417</td>\n",
       "      <td>122.1042</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Ro...</td>\n",
       "      <td>Santa Maria, Romblon 12.4417,122.1042</td>\n",
       "      <td>Santa Maria, Romblon</td>\n",
       "      <td>12.4417 , 122.1042</td>\n",
       "      <td>12.441667 , 122.104167</td>\n",
       "      <td>Santa Maria, Romblon 12.4417,122.1042</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>Santa Cruz</td>\n",
       "      <td>Santa Cruz, Chile</td>\n",
       "      <td>-34.6372</td>\n",
       "      <td>-71.3631</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Cruz%2C_Chile</td>\n",
       "      <td>Santa Cruz, Chile -34.6372,-71.3631</td>\n",
       "      <td>Santa Cruz, Chile</td>\n",
       "      <td>-34.6372 , -71.3631</td>\n",
       "      <td>-34.633333333333 , -71.366666666667</td>\n",
       "      <td>Santa Cruz, Chile -34.6372,-71.3631</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0         name           wikipedia_name      lat       lon  \\\n",
       "0           0  Santa Maria  Santa Maria, California  34.9514 -120.4333   \n",
       "1           1  Santa Maria     Santa Maria, Bulacan  14.8208  120.9636   \n",
       "2           2  Santa Maria  Santa Maria, Ilocos Sur  17.3585  120.4981   \n",
       "3           3  Santa Maria     Santa Maria, Romblon  12.4417  122.1042   \n",
       "4           4   Santa Cruz        Santa Cruz, Chile -34.6372  -71.3631   \n",
       "\n",
       "                                                page  \\\n",
       "0  http://en.wikipedia.org/wiki/Santa_Maria%2C_Ca...   \n",
       "1  http://en.wikipedia.org/wiki/Santa_Maria%2C_Bu...   \n",
       "2  http://en.wikipedia.org/wiki/Santa_Maria%2C_Il...   \n",
       "3  http://en.wikipedia.org/wiki/Santa_Maria%2C_Ro...   \n",
       "4   http://en.wikipedia.org/wiki/Santa_Cruz%2C_Chile   \n",
       "\n",
       "                                  wiktor_key                   toponym  \\\n",
       "0  Santa Maria, California 34.9514,-120.4333  Santa Maria, California    \n",
       "1      Santa Maria, Bulacan 14.8208,120.9636     Santa Maria, Bulacan    \n",
       "2   Santa Maria, Ilocos Sur 17.3585,120.4981  Santa Maria, Ilocos Sur    \n",
       "3      Santa Maria, Romblon 12.4417,122.1042     Santa Maria, Romblon    \n",
       "4        Santa Cruz, Chile -34.6372,-71.3631        Santa Cruz, Chile    \n",
       "\n",
       "               old_coord                             new_coord  \\\n",
       "0   34.9514 , -120.4333     34.951388888889 , -120.43333333333   \n",
       "1    14.8208 , 120.9636      14.818333333333 , 120.95633333333   \n",
       "2    17.3585 , 120.4981      17.366666666667 , 120.48333333333   \n",
       "3    12.4417 , 122.1042                 12.441667 , 122.104167   \n",
       "4   -34.6372 , -71.3631    -34.633333333333 , -71.366666666667   \n",
       "\n",
       "                                     mlg_key  \\\n",
       "0  Santa Maria, California 34.9514,-120.4333   \n",
       "1      Santa Maria, Bulacan 14.8208,120.9636   \n",
       "2   Santa Maria, Ilocos Sur 17.3585,120.4981   \n",
       "3      Santa Maria, Romblon 12.4417,122.1042   \n",
       "4        Santa Cruz, Chile -34.6372,-71.3631   \n",
       "\n",
       "                           geocoded_coordinates_list  median_error_distance  \n",
       "0  [(34.95303, -120.43572), (34.95303, -120.43572...               0.283927  \n",
       "1  [(14.83333, 121.0), (14.83333, 121.0), (14.833...               4.981280  \n",
       "2  [(16.9518, 120.4445), (16.9518, 120.4445), (16...              46.315270  \n",
       "3                                                 []                    NaN  \n",
       "4                                                 []                    NaN  "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_wiktor_edinburgh = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/wiktor_geocoded_results_Edinburgh_Geoparser.csv')\n",
    "df_wiktor_edinburgh.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "d649dcfc",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAF3CAYAAACSWTEVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAArMklEQVR4nO3debgcVbnv8e/PhAABGcOcAEHAQ+CgQGQUCIIa0QvH4TKISpRB5TB4HBk8yOBFGRxQ8DAoIMgoBzEqCKIElDkgU4BAiJEkIARkkDmB9/6xVpNKp3vv2nt3966kf5/n6Wd3r1pd9XZ177dXr1q1ShGBmZkNvrcNdgBmZpY4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE3ILSJohaVJvZdYzSetKCknH9FTWj/VOkjSjBSG2jaQJ+XWOG+xYbPB0bUKWND7/A3y7wbKt87LXJA1vsPz3kt6UNGIA268lmlK3kusc1+C5L0q6S9J/SRra33hLbHtSL6/hm+3ati1I0hBJn5J0laR/SHpd0r8k3SPpx5LeM9gxWmNt+wddBPwFmAeMa7Bsp7xsGLAtcF1tQU5q7wXuj4inc/E7gb6eYTMH+HRd2ceAjwInAA/2cX1FFwNXAQJWBz4DfB/YCDhwAOvtzWvA/k2W3d3Pdf4dWJr0flgvciPhSmA74E7gf4CZwFLAGGB34GBJ20bELYMVpzXWtQk5Il6UdAfwHknDI+LlwuJxwB+AzfL96wrL3gMsA0wqrOu1fmz/JeAXxTJJ65MS8h8iYlKj55V0V0S8tW5JPwEeAvaXdFREzBnAunsyr7jdVoh0KumrrVxnu0l6e0T8axC2K+CXpGR8SESc1qDOl4DPkr48FyuSlgbmRkRHvrzbsb2u7bLIrie1grerFeQW8HbADfm2U91zxhWeW3tOqf5iSaMlTZX0uKR3lQ1S0ghJp0uamX9+zsyPVy7z/Jz8byW1mN9RdrvtJOm9km6S9IqkJyWdBizboF6P/cqSPiLpDkmvSnpC0snNumYkrSfp15Kel/SCpF9JWq+uTtO+3EZ90bX3XtJmkq6R9Dxwb2H5x3NXwauSHpP0LUm75G1MaBDm2yR9VdKjucvsYUn79rQvCz5C+nxe1CgZA0TEvIg4OyLuqnsdS0o6UtKUHOtzkn4jabMG+2EZSd8pxPgPSedLWqeuXq0LbYKkQ/JreTX/PaTBereUdF5e/nLuZrlJ0kcb1D0vr3sVSedIehJ4CRiZl39G0u35dbwkabqkCyWtUreeDSRdkD87r+f382RJy/Rle63StS3k7HrgSOa3iGF+C/gG4AXgVEnL5KRGrht5eWmSNid1IzwLbBMRfy/5vOWBm4H1gXOAu0gt9y8C75O0ZcnWWC0R/7MvcfeVmverP1drSUjaivSr41/AicBzwF7A+X3c3K7AQcAZpH2zO/BV0j4+oa5u7VfNbcARwAb5uVtL2iwi/tHHbRetDfyJ1Dr9X/IXi6Q9Sd1HjwLHkrpd9gX+Tw/rOoHURXMmqRX7ReA8SdMi4qZe4vhE/vvTvgQvaQng96TuuQuA04DlgQOAmyTtEBGTC3WvITVaLge+R9qXXwQ+IGlsRMyq28QhpK6zM0nv+d7AjyStFBHHFup9FPg34DJSV9XKpP11haR9IuKiBuH/AfgHcDzpPX5R0qeBnwN/Bo4GXgFGkT4vq5K6C5G0Bel9ey7HNht4F3AosJ2kHSNibm/ba7pj+yMiuvZG+uC/BtxcKDuC9KEZSupzDeADednQ/AbcXbeeGcCkZmXA+/M6bwZW7iGeY/L2xhXK/l8uO6iu7n/m8uMLZeNy2dHACGAV4N+B03P5bW3cl5PyNprdxhbq3gy8DmxYKBsG3J7rHlMoX7eHspeAdQvlAu4HnmgS2w/ryj+ay88olE2ofw/q1jOjwfscwP515UNJ/+BPAisWypcFpufnTGiw3b8Cwwrla+XP6MUl3oM78zpWbLBsRN1tmcKy/8rP+2Ddc5YDHqPw2SYl6QBOqqv74Vx+QYPP47+AkQ3e67l15cs0iHs4MBV4oK78vLzuXzR4zhWkxtTQXvbXPaSuvLc3+VxMKLO9Vt66ussiIl4htZjGFn6ijCMl6HkR8SDwFPO7KWqt5+spSdKngN+Rvol3john+hjmR0nf6GfVlZ+Zyxf6OUdqjc0hxX4vqSV4BakF2U6vkr58Gt2mAkhaFdgG+HVEPFx7YkS8Dvygj9u7MiJmFNYRpPdmdUkLdX8A3y0+iIhf5bj+o4/brfdP4Ny6si2ANYHzIuLZwjZfJLXom/lJ3he1+rOBh0mt0N4sl/++UCzM+2JO3e3EQpVPkRLTnUrdYyPyL51hpBbhe5X6SyF93t4EvlPcRkT8jnTgdndJ9Xnlwii0mgvv9VAKvxZi/q9QJA1X6pIbTvrf2UjScizslAZlz+fnfViSGixH0r8DmwIXAUvWve6/kL7sP1Byey3T7V0WkP6Btyd96P5I+ilW/Ll7I/P7kcflv5NKrnsLYAfST7yPRcQb/YhvNDA56g4cRMQ8SQ8Dmzd4zlmkn89LkFrI3yD1dZU6OJb72YYUit6IcgcC34iI63qpU+uzfajBsgfKxFcwvUFZ7QtvZRb8OflcNO6WeBD4j7puqb56tMF7Ozr/ndqgfqOymmavaZ0G5fVqiXg5UrdNzSukL0WA1ag7mEz6Jbg0+ad8EyNIozVGA48Xv2QKpgDvznWfKpQ3GjFUe6/f6sPPX9bfJjUcVm3wnBWo+7IhfVnVO4H0f3cl8IykG4CrgUtjfvfeRvnvsfnWyGoNyhptr2WckFNCPpqUbF9gfv9xzQ3AD3IrYxypdXBjyXU/QvpZthMwntRS7oRHConxakl/IX3rn0Hqq+3NHSyYAP5O6iaomp6+4Bq2jEroafhis/+Xl5uU90ez11Tm9dxP+oJ+N4VfcfnL4jpIB0SbrPs+4Ms9rLtdI3NSAKkley0pUZ4KTCa1dN8gjQr5JA0GIcSCo6NqZY9IGgPsnG87AmcDx+b+8EeZvz+/R+o/b2ShL51G22slJ2S4hdRy3ImUkF8hJaSaG0j7aRyp9Xx3k9ZBIy8Au5He8Csk7RERv+5jfNOBd0oaWmwlK40k2JDGLaoFRMTNki4APiPpRxFxcy9P2YfUYqp5pY8x9+Rv+e+/NVg2poXbqbeCpNUbtJI3Ap4qtI5rBz1XarCO0aQv2DJm5L/vbLCsUVkrXE4ac74/fehWIzUcVgH+FBFv9lJ3OjBe0goR8VzdsjGkz/zTdeUbsbDae137/G5KOqB2XER8q1hRUrOx7U1FGop6Vb4haVdSg+jLpOMvj+SqZX7VdUxX9yHDW2/cLaTuhY8AtxT78EitjmeAr1E3/rjk+l8g9UXdBvxS0sf7GOKVpH+W+g/lAbn8VyXXczyptXFcbxUj4qaIuK5w6+3ofmkR8SRpCN7ukjaslUsaRjq41E6HFx/k4VTvJO3jmtpP0l3q6u5N6hMuazLwBDBB0oqF9SwLfKEP6+mL35IaEJ+UdHCTOo1a2ueTRkE0bCFLKv50v5KUN+r35YdIo38mNkjq+0gaWahbe6/fyDHD/F8GC8QnaRMaHydpqslIn9owv9oX7V9J/9tfUN3Qx7yOoZIafSm3lVvIyfWkFvK2wALfzhERkv7M/AM/fWl51NbxYv7A/ga4RNKnIuLSkk8/Cfi/wOl56NxfSR/8/Uh9kSeVjGGapEtI/xzbR8Sf+/o6ShiaD2I2Mr3QMv8y6YvtJkmnM3/YWzs/j08DH5O0Zt52bdjbk6TRLQBExFRJ1wGfzz+j7yZ1AXwUmEbql+9V7uP/KnAhcLukn5GGvU0gfcGPpu9nd/a2zZD0CVLS/LHSOOffkPp+h5OGPu6Rq88oPPVUUh/zyZLeRzqI9gJpON/OzP8FCWm0wb7AN3L3x42kIZm1fXlkg9AeBm6TdAZpxMUnSQfIj4+ImbnOg6Q+6K8rTVcwlfQL8POk7pQt+rArrpX0HGnY20xS3/ME8igQeGtffTq/1nslnZO3Pzy/no+RRlyd14ftDlw7h3AsKjfSqdC14Vk7Nlj+pbxsHrB8g+Uz6GHYW6FsadIBvnnAPg3WcwwNhlyRWsI/AWaRfjLPIg1lG1FXb1x+/lebvM6NSC2R69uwDycV9mGj2y/q6u9AGv72Kukf+XRgE/o27O2YBnHU9uG6dbHNIB1A+jUp2fwr31+/wTpWJx0UfYF0YPDqvO8m0XjY26Qe9sv/JY10eY00hOxbzB9WtUeh3oRG730x/j68F0NJp+Vfnfft3Px67837eWyT5xxK6q57Kd8eIX2hfKCu7jKkURbTScMXnyIlunWafB4n5HU/kvfDI8BhDWJYJ+/3OaR++dvzvmr0np5HHljTYD0HMH+88OukXypXATs12eYZ+X18nfRleWd+faPKbK+VN+WNmVmHSPoKafjUNhFx62DH0y5KZzteD3w2Is4b1GAWEV3fh2zWLpKGSRpSV7Ys6aDSM8zv1zQDOpiQ8/nfT0m6v8lySfqRpGmS7s39pWaLsvWAaUrzPhwo6VukboPRwDdjwYPHZh09qHce6Rz5ZvMVfIh0oGUDYCvStIFbdSQys/aYQxpRsg/pRId5pANUh0fEZYMZmFVTR/uQ81HZ30bEJg2WnUk6OHJxfjyVdIDjiY4FaGY2iKo07G0t0hCVmlm5bKGELOlA8kTrY8aM2WLKlCkdCdDMrEUannm5SB7Ui4izImJsRIxdeumle3+CmdkioEoJeTZpztKakbnMzKwrVCkhTyTNtSBJWwPPu//YzLpJx/qQJV1MOnNnhKRZpDOWlgCIiDNIZ9LsSjo99WXSDE9mZl2jYwk5IvbuZXmQBsybmXWlKnVZmJl1NSdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4roaEKWNF7SVEnTJB3eYPnakq6X9FdJ90ratZPxmZkNpo4lZElDgNOBDwFjgL0ljamr9k3gsojYDNgL+Emn4jMzG2ydbCFvCUyLiOkR8TpwCbB7XZ0Alsv3lwce72B8ZmaDamgHt7UWMLPweBawVV2dY4BrJR0CLAPs0pnQzMwGX9UO6u0NnBcRI4FdgQskLRSjpAMlTZY0ec6cOR0P0sysHTqZkGcDowqPR+ayov2AywAi4hZgKWBE/Yoi4qyIGBsRY1dZZZU2hWtm1lmdTMh3ABtIGi1pGOmg3cS6Oo8BOwNI2oiUkN0ENrOu0LGEHBHzgIOBa4AHSaMppkg6TtJuudpXgAMk3QNcDEyIiOhUjGZmg0mLer4bO3ZsTJ48ebDDMDPrCzUqrNpBPTOzruWEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRpRKypE7Om2xm1pXKtpCfkHRKnoHNzMzaoGxCPhLYFrhf0i2S9pO0bBvjMjPrOqUSckScHRHbApsAfwG+TWo1nyNpu3YGaGbWLfp0UC8iHoyIr5Gu9nEk8EngRkkPSfpCo8stmZlZOX06WJev9PEx4HPA+0it5Z8BawL/DYwjXQnEzMz6qNQE9ZI2JyXhvYG5wPnATyPi4UKdjYHJEbF0m2JtaNiwYTF37tx+P3/NkaOYPfOxFkZkZtarhhPUl20h3wFcCxwI/DpfjqneDOCSfoU2AHPnzmXPM2/u9/Mv/fy2LYzGzKz/yibk9SLi7z1ViIiXgM8OPCQzs+5U9iDcupJ2rC+UtKOkHVock5lZVyqbkH8ArNigfLm8zMzMBqhsQn4ncE+D8vvzMjMzG6CyCfkVYI0G5WsBr7cuHDOz7lU2IV8DnCjprW4LSSsB38nLzMxsgMqOsvgqcCMwQ9K9uWxT4Clgz3YEZmbWbUol5Ih4QtK7gH2Ad+finwMXRcTLbYrNzKyrlD51Oifes9sYi5lZVyudkCWNBHYAVqWu7zkivt/iuMzMuk6phCxpH+AcYB4wByhOgBGAE7KZ2QCVbSEfB3wP+O+IeKON8ZiZda2yw95WI83u5mRsZtYmZRPyVcBW7QzEzKzble2y+APpxJCNgftIcyK/JSKuaHVgZmbdpmxCPjP/PbLBsgCGtCYcM7PuVfbEEF8rz8yszZxozcwqolRCVnKQpCmSXpa0Xi4/XNIe7Q3RzKw7lG0hHwZ8EziLBS/ONxs4uNVBmZl1o7IJ+QvAARFxKulsvZq7gI1bHpWZWRcqm5DXIV0dpN5cYOnWhWNm1r3KJuTpwOYNyncFHmhdOGZm3avsOORTgNMkDSf1IW8j6dPA14HPtSs4M7NuUnYc8rmShgInAMOBC4DHgUMj4tI2xmdm1jX6MkH92cDZkkYAb4uIp9oXlplZ9ymdkGsi4ul2BGJm1u3KTlB/HwtOSr+AiNi0ZRGZmXWpsi3ky+seL0G62Ol2wOmtDMjMrFuVPah3bKNySV8jjVE2M7MBGujkQlcA+7QiEDOzbjfQhLwD8HIrAjEz63ZlD+pNrC8C1gA2Axp2Z5iZWd+UPaj3TN3jN4EpwJERcW1rQzIz605lD+p9tt2BmJl1u45eMUTSeElTJU2TdHiTOntIeiBPhn9RJ+MzMxtMZfuQ/0YPJ4YURcR6TdYxhDRm+f3ALOAOSRMj4oFCnQ2AI4DtIuJZSauW2aaZ2eKgbB/yecCXgduBW3LZNsCWwPeBF0usY0tgWkRMB5B0CbA7C07feQBwekQ8C+D5Msysm5RNyKOBEyPihGKhpCOAjSPiUyXWsRYws/B4FrBVXZ0N83pvAoYAx0TE7+tXJOlA4MCSsZuZLRLKJuSP0XiC+l+SuhhaZSiwATAOGAncKOnfI+K5YqWIOIt0fT8klepKMTOrurIH9V4iJcl64yh/YshsYFTh8chcVjQLmBgRcyPib8DDpARtZrbYK9tC/gFwuqSxwK25bGtgX+CYkuu4A9hA0mhSIt4L+GRdnSuBvYFz87zLG5IuH2VmttgrOw75JEkzgMOAPXLxg8C+EXFZyXXMk3QwcA2pf/iciJgi6ThgckRMzMs+IOkB4A3gaxFRf1KKmdliSRGLdhespNjzzJv7/fxLP78ti/o+MLNFjhoVlj4xRNJSkj4h6RuSVshl75C0UosCNDPramVPDFkfuA5YFliBNLriOeCL+fH+bYnOzKyLlG0h/xC4FlgNeKVQPhHYqcUxmZl1pbKjLLYFto6IN6QFuj4eA9ZseVRmZl2oL5MLLdGgbG3g+RbFYmbW1com5GtJc1nUhKTlSJPT/67lUZmZdaGyXRZfAf4kaSqwFHApsD7wJPPHJZuZ2QCUPTFktqR3k86i25zUsj4LuDAiXunpuWZmVk6vCVnSEqRZ2naOiHOAc9oelZlZF+q1Dzki5gJzKTlBvZmZ9U/Zg3o/Bo6QVLbP2czM+qhsgt0e2BGYLel+0nScb4mI3VodmJlZtymbkJ8G/redgZiZdbseE7Kk3YCrI+KzHYrHzKxr9daH/CtgxdoDSfdJGtVDfTMz66feEnL9nJ3r0vgUajMzG6C+zGVhZmZt1FtCDhYcf1z/2MzMWqS3URYCbpA0Lz8eDlwt6fVipYjYtB3BmZl1k94S8rF1jz30zcysTXpMyBFRn5DNzKxNfFDPzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpqOspD0mbIriYjzWxOOmVn36mnY2+l1j4eR5rF4Mz9+G+lKIq8BTshmZgPUtMsiIt5euwF7AfeSJqpfKt+2B+4GPtmBOM3MFntl+5BPAQ6NiJsiYl6+3QR8Cfhe26IzM+siZRPyutRdtil7GVi7ZdGYmXWxsgn5NuBHktaqFeT7PwBubUdgZmbdpmxC3g9YGZghaYakGcAMYFXggPaEZmbWXUpd5DQiHpW0KfB+4N9y8YPAdRHh+ZHNzFqg7FWnyYn32nwzM7MWK32mnqSDJE2R9LKk9XLZ4ZL2aF94Zmbdo1RClvQl4JvAWSx44dPZwMGtD8vMrPuUbSF/ATggIk4F5hXK7wI2bnlUZmZdqGxCXge4v0H5XGDp1oVjZta9yibk6cDmDcp3BR5oXThmZt2r7CiLU4DTJA0n9SFvI+nTwNeBz7UrODOzblJ2HPK5koYCJwDDgQuAx0nzW1zaxvjMzLpGX8Yhnw2cLWkE8LaIeKp9YZmZdZ+yw96OlrQTQEQ8XUvGkpaRdHQ7AzQz6xZlD+odA/xe0iF15csC32ppRGZmXaov19TbDzhW0k8kDWlXQGZm3aovCflaYBvgfaTW8vLtCcnMrDuVTcgBEBFTga1z2e3Ahu0IysysG5VNyG/NXxERzwHjgT8Av29DTGZmXanssLdjgRdrDyLiDeBgSXcCO7QjMDOzblP2xJBjm5SfC5zb0ojMzLpU04Qs6cvATyLi1Xy/mYiIH7Q+NDOz7tJTC/kQ4OfAq/l+M0G62GmvJI0HTgWGAD+NiO82qfdx4HLgPRExucy6zcwWdU0TckSMbnS/v/LY5dNJ1+WbBdwhaWJEPFBX7+3AYaQrXZuZdY2+jEMeqC2BaRExPSJeBy4Bdm9Q73jgRFLL3Mysa/TWh1xKRHy/RLW1gJmFx7OAreq2uTkwKiJ+J+lrPcR2IHBg2fjMzBYFvfUhlxFAmYTcI0lvy+uZ0OsGI84iXd8PSTHQbZuZVUGpPuQWmQ2MKjwemctq3g5sAkySBLA6MFHSbj6wZ2bdoJN9yHcAG0gaLWkYsBcwsbYwIp6PiBERsW5ErAvcCjgZm1nXKD1BvaQVgQ8BawPDissi4rjenh8R8yQdDFxDGvZ2TkRMkXQcMDkiJva8BjOzxVuphCxpa+B3wGvAKqSuhjXy4xlArwkZICKuAq6qK2s4wX1EjCuzTjOzxUXZLouTgQtJIyVeJU3BuTYwmTREzczMBqhsQt4UOC0iAngDWDIingS+QbqaiJmZDVDZhPx64f6TwDr5/ovAmi2NyMysS5U9qHcX8B7gYWAS8G1JqwGfAu5tT2hmZt2lbAv5KODxfP+bwBzgx8CK+Iw5M7OWKDsf8uTC/Tmk4W9mZtZCnTwxxMzMelB2HPKKpNEUOwGrUpfII2LVlkdmZtZlyh7UOx/YmDRh/ZPkq1CbmVnrlE3I44AdI+KuNsZiZtbVyvYhP9qHumZm1g9lk+xhwHckvStfisnMzFqsbJfFNGBp0gki5PmK3xIRTtJmZgNUNiFfDCwPHIoP6pmZtUXZhDwW2DIi7m9nMGZm3axsH/IDwHLtDMTMrNuVTcjfBL4vaRdJq0laqXhrZ4BmZt2ibJdF7Sof17Jg/7HyYx/UMzMboLIJeae2RmFmZr0nZElLACcBn4mIqe0PycysO/XahxwRc4HReKibmVlblT2o93PggHYGYmbW7cr2IS8D7CPp/cCdwEvFhRFxaKsDMzPrNmUT8kbk06aB9eqWuSvDzKwFyl7CyaMszMzarGwLGQBJSwHrk1rFj0bEq22JysysC5U6qCdpCUknA88C9wD3Ac9KOikPizMzswEq20I+Edgb+ALwl1y2PfAdUlL/autDMzPrLmUT8ieBz0XEVYWyRyXNAX6KE7KZ2YCVHYe8POkyTvUeBVZoWTRmZl2sbEK+hzQ5fb3DgLtbFo2ZWRcr22XxdeAqSbsAt+ayrYE1gQ+1IzAzs25TqoUcETcCGwKXA8vm2y+Bd0bEX3p6rpmZlVN6HHJEPA4c1cZYzMy6Wo8JuezVQCLin60Jx8yse/XWQn6a3ueqiBLrMTOzXvSWSHuaw2I8aZTFvNaFY2bWvXpMyBFxQ32ZpM2Ak0ln6p0JHN+e0MzMukvZcchIGi3pIuB24BlgTEQcGhFz2hadmVkX6TUhS1pZ0qnAQ8DqwLYRsWdENDpzz8zM+qnHhCzpKNLp0TsCu0fE+yLijo5EZmbWZXo7qHc88AowCzhI0kGNKkXEbq0OzMys2/SWkM/Hl2gyM+uI3kZZTOhQHGZmXa/0KAszM2svJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4roaEKWNF7SVEnTJB3eYPmXJT0g6V5Jf5S0TifjMzMbTB1LyJKGAKeTrlI9Bthb0pi6an8FxkbEpqQLqp7UqfjMzAZbJ1vIWwLTImJ6RLwOXALsXqwQEddHxMv54a3AyA7GZ2Y2qDqZkNcCZhYez8plzewHXN3WiMzMKqSSFyeV9ClgLGke5kbLDwQO7GhQZmZt1smEPBsYVXg8MpctQNIuwFHAjhHxWqMVRcRZwFm5vqcHNbPFQie7LO4ANsjX5hsG7AVMLFbIF1A9E9gtIp7qYGxmZoOuYwk5IuYBBwPXAA8Cl0XEFEnHSapdceRkYFngl5LuljSxyerMzBY7He1DjoirgKvqyo4u3N+lk/GYmVWJz9QzM6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4pwQn7bUCT167bWqLUHO3ozW4wMHewABt2b89jzzJv79dRLP79ti4Mxs27mFrKZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQT8kAM4Cw/n+lnZvV8pt5ADOAsP/CZfma2ILeQzcwqwgnZzKwinJDNzCrCCdn6Za1Ra/uAplmL+aCe9cvjs2b6gKZZi7mFPJg8bM7MCtxCHkyDOGxurVFr8/ismf1+vpm1nhNyl3KXg1n1OCEvynKXh5ktHjqakCWNB04FhgA/jYjv1i1fEjgf2AJ4BtgzImZ0MsZFyqJ8PcABfJmsOXIUs2c+1u9ND7S7ZqDbt85bVN7zjiVkSUOA04H3A7OAOyRNjIgHCtX2A56NiPUl7QWcCOzZqRitgwbxy2TA3TVf3KHfXyZDlliSN+a+1u9tD/T53fplMpjvOZTf751sIW8JTIuI6QCSLgF2B4oJeXfgmHz/cuA0SYqI6GCcVnWD3VUzwC+Tgfbdd+uXyUCfPyADPQBft9+bpTR1KtdJ+gQwPiL2z48/DWwVEQcX6tyf68zKjx/NdZ6uW9eBwIH54SbA/R14CWWMAJ7utVb7OY6FVSWWqsQB1YmlKnFA52J5OiLG1xcukgf1IuIs4CwASZMjYuwghwRUJxbHsbCqxFKVOKA6sVQlDhj8WDp5YshsYFTh8chc1rCOpKHA8qSDe2Zmi71OJuQ7gA0kjZY0DNgLmFhXZyKwb77/CeBP7j82s27RsS6LiJgn6WDgGtKwt3MiYoqk44DJETER+BlwgaRpwD9JSbs3Z7Ut6L6rSiyOY2FViaUqcUB1YqlKHDDIsXTsoJ6ZmfXMkwuZmVWEE7KZWUUs0glZ0nhJUyVNk3R4G9Y/StL1kh6QNEXSYbl8JUl/kPRI/rtiLpekH+V47pW0eWFd++b6j0jat9k2e4lniKS/Svptfjxa0m15e5fmg6VIWjI/npaXr1tYxxG5fKqkD/YzjhUkXS7pIUkPStpmMPaJpP/K78v9ki6WtFSn9omkcyQ9lcfO18patg8kbSHpvvycH0mNz+ZoEsfJ+b25V9KvJK3Q22tt9r/UbH+WjaWw7CuSQtKIwdgnufyQvF+mSDqpE/ukzyJikbyRDgw+CqwHDAPuAca0eBtrAJvn+28HHgbGACcBh+fyw4ET8/1dgasBAVsDt+XylYDp+e+K+f6K/Yjny8BFwG/z48uAvfL9M4Av5vsHAWfk+3sBl+b7Y/J+WhIYnfffkH7E8XNg/3x/GLBCp/cJsBbwN2Dpwr6Y0Kl9AuwAbA7cXyhr2T4Abs91lZ/7oT7E8QFgaL5/YiGOhq+VHv6Xmu3PsrHk8lGkg/l/B0YM0j7ZCbgOWDI/XrUT+6TP/1utWlGnb8A2wDWFx0cAR7R5m78mzcUxFVgjl60BTM33zwT2LtSfmpfvDZxZKF+gXsltjwT+CLwP+G3+UD5d+Md7a3/kD/82+f7QXE/1+6hYrw9xLE9KhKor7+g+ISXkmfkfd2jeJx/s5D4B1q37p2/JPsjLHiqUL1Cvtzjqln0UuLDR/0jttdLkf6mnz1hfYiFNg/AuYAbzE3JH9wkpie7SoF7b90lfbotyl0XtH7JmVi5ri/wTdzPgNmC1iHgiL/oHsFovMbUi1h8CXwfezI9XBp6LiHkN1vnW9vLy53P9VsQxGpgDnKvUffJTScvQ4X0SEbOBU4DHgCdIr/FOBmef1LRqH6yV77cips+RWpP9iaOnz1gpknYHZkfEPXWLOr1PNgS2z10NN0h6Tz/jGPA+6cminJA7RtKywP8CX4qIF4rLIn1NtnXsoKSPAE9FxJ3t3E5JQ0k/B/8nIjYDXiL9PH9Lh/bJiqTJqEYDawLLAAvNDTBYOrEPeiPpKGAecOEgbX84cCRw9GBsv85Q0q+prYGvAZc164MeTItyQi5zKvaASVqClIwvjIgrcvGTktbIy9cAnuolpoHGuh2wm6QZwCWkbotTgRWUTjGvX2ezU9Bbsc9mAbMi4rb8+HJSgu70PtkF+FtEzImIucAVpP00GPukplX7YHa+3++YJE0APgLsk78c+hPHMzTfn2W8g/SFeU/+7I4E7pK0ej9iGeg+mQVcEcntpF+aI/oRx0D3Sc9a1ffR6RvpG2866Q2vdbpv3OJtiDRh/g/ryk9mwYM3J+X7H2bBAxW35/KVSP2uK+bb34CV+hnTOOYf1PslCx5cOCjf/08WPIB1Wb6/MQsewJhO/w7q/Rl4Z75/TN4fHd0nwFbAFGB4XvfPgUM6uU9YuJ+yZfuAhQ9g7dqHOMaTprVdpa5ew9dKD/9LzfZn2Vjqls1gfh9yp/fJF4Dj8v0NSd0R6sQ+6dP/VqtWNBg30pHah0lHQ49qw/rfS/rZeS9wd77tSupH+iPwCOnIbe0DI9Ik/I8C9wFjC+v6HDAt3z47gJjGMT8hr5c/pNPyh6R2BHmp/HhaXr5e4flH5fim0uQodYkY3g1MzvvlyvyP0/F9AhwLPESafvWC/E/VkX0CXEzqu55Lan3t18p9AIzNr+tR4DTqDqL2Esc0UsKpfWbP6O210uR/qdn+LBtL3fIZzE/Ind4nw4Bf5OffBbyvE/ukrzefOm1mVhGLch+ymdlixQnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2RZJkk6TNGmw46gnaUVJT0p6R348rjjLWQu385+SftPKddrgc0K2jpN0Xk5SP2uw7MS87Ldt2O6kvO762yWFOsXylyVNl3SRpO1LbuZI4KqIeLTV8df5KbBFH+KyRYATsg2WmcAeeWIi4K1Tmj9DmjCoXc4lzRxWvH2+rs4BuXwj0kkFrwM3SPpaTyvOczfsT7o2ZFtFxGukqVgPbfe2rHOckG2w3Es6o22PQtmHgVeBScWKShPznyLp2Xz7Ien01mKdSZJ+IukESU/nCcpPkVT/GX85Iv5Rd3u+rs5zufzvEXF9REwAvgt8R9L6PbymXUlndt7UrILSRPm/knSXpFUlrZtb43vlWcheybPobSppE0k3S3pJ0l8kja5b3UTSHCfDe4jJFiFOyDaYfkY6Tbbmc6QWbP3po18htVo/T5p/dgiwT4P17UOa3Wxb4GDgS8CeLYr1e6T/l//ooc72wJ3R5PRXScsBvyfN1zAuIp4qLD6WNJn8ZsBzpNN/f0w6rXdL0qnfP6pb5WTSnAvb9O2lWFU5IdtguggYK2mDPAPYeOC8BvW+RJqo57KIeAg4jDTfcL0HIuLoiHg4Ii4Drgd2rqtzoKQX624H9RZoRDxDmr1tvR6qrQM83mTZqjmefwEfjLppXIHvR8RV+fV9j3Qlix/nFvoU0twNO9XF9DJpXud1e4vfFg1De69i1h4R8aykX5Faxs8BkyLiseI0tZKWJ/Xn3lJ43puSbmPB6REhdYMUPU5KhEWXklqjRXNKhix6nuN4aeDJJsuuIU1q8/GYP7l5UTH22jruqytbRtLwnIhrXsnbtcWAE7INtnNIU2e+yMAnMp9b9zhY+Ffg8xExra8rzsPWViFNydjM06SZ7xr5Lam/fBPSDGz1irFHD2X1r2clyn+hWMW5y8IG2x9JoxhGkKbyXEA+4PYEaR5cIF2xmNSv2klfIU1qfmUPdf5K6mpo5L9Jc+f+UdK7WxFQHuu8FKnlbYsBt5BtUEVESNqUNLfta02qnQocIelh0s/4g0jdGE80qd+T4bm/uuj1iPhn4fEKuc4w0lUv9iUNx/t6L+OLrwFOlLRy7nNeQEQclb9MrpO0cyx8nbm+2h6YHhGPDHA9VhFOyDboIuJfvVT5HrA66WQISJPRX0gaJ9xXn823optIFyOoOTv/fY2U9G8ljYq4sacVR8R9km4nXZHk9CZ1jsxJ+Y+SdiYdlOuvvQux2mLAE9SbtZCk8aQW/ZiIeKON29mE1N2zYYNx1LaIch+yWQtFxO9JreORvdUdoDWBzzgZL17cQjYzqwi3kM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCL+P6lQXjHTJr3QAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_wiktor_edinburgh, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('WikToR - Edinburgh Geoparser', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bb7aeb85",
   "metadata": {},
   "source": [
    "## Reading Toponym Resolution Results of CamCoder"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "a189e35b",
   "metadata": {
    "scrolled": true
   },
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>lgl_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Alexandria</td>\n",
       "      <td>31.3113</td>\n",
       "      <td>-92.4451</td>\n",
       "      <td>Alexandria 31.3113,-92.4451</td>\n",
       "      <td>Alexandria</td>\n",
       "      <td>31.3113 , -92.4451</td>\n",
       "      <td>31.2928 , -92.4592</td>\n",
       "      <td>Alexandria 31.3113,-92.4451</td>\n",
       "      <td>[(31.31129, -92.44514), (31.31129, -92.44514),...</td>\n",
       "      <td>671.350717</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Rapides Parish</td>\n",
       "      <td>31.1669</td>\n",
       "      <td>-92.4835</td>\n",
       "      <td>Rapides Parish 31.1669,-92.4835</td>\n",
       "      <td>Rapides Parish</td>\n",
       "      <td>31.1669 , -92.4835</td>\n",
       "      <td>31.2 , -92.54</td>\n",
       "      <td>Rapides Parish 31.1669,-92.4835</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Avoyelles Parish</td>\n",
       "      <td>31.0669</td>\n",
       "      <td>-91.9668</td>\n",
       "      <td>Avoyelles Parish 31.0669,-91.9668</td>\n",
       "      <td>Avoyelles Parish</td>\n",
       "      <td>31.0669 , -91.9668</td>\n",
       "      <td>31.07 , -92.0</td>\n",
       "      <td>Avoyelles Parish 31.0669,-91.9668</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Cottonport</td>\n",
       "      <td>30.9841</td>\n",
       "      <td>-92.0535</td>\n",
       "      <td>Cottonport 30.9841,-92.0535</td>\n",
       "      <td>Cottonport</td>\n",
       "      <td>30.9841 , -92.0535</td>\n",
       "      <td>30.98408 , -92.05346</td>\n",
       "      <td>Cottonport 30.9841,-92.0535</td>\n",
       "      <td>[(30.98408, -92.05346)]</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>490</td>\n",
       "      <td>Pineville</td>\n",
       "      <td>31.3224</td>\n",
       "      <td>-92.4343</td>\n",
       "      <td>Pineville 31.3224,-92.4343</td>\n",
       "      <td>Pineville</td>\n",
       "      <td>31.3224 , -92.4343</td>\n",
       "      <td>31.3224 , -92.4343</td>\n",
       "      <td>Pineville 31.3224,-92.4343</td>\n",
       "      <td>[(42.15536, -75.1635)]</td>\n",
       "      <td>1947.840891</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0              name      lat      lon  \\\n",
       "0           0        Alexandria  31.3113 -92.4451   \n",
       "1           1    Rapides Parish  31.1669 -92.4835   \n",
       "2           2  Avoyelles Parish  31.0669 -91.9668   \n",
       "3           3        Cottonport  30.9841 -92.0535   \n",
       "4         490         Pineville  31.3224 -92.4343   \n",
       "\n",
       "                             lgl_key            toponym             old_coord  \\\n",
       "0        Alexandria 31.3113,-92.4451        Alexandria    31.3113 , -92.4451    \n",
       "1    Rapides Parish 31.1669,-92.4835    Rapides Parish    31.1669 , -92.4835    \n",
       "2  Avoyelles Parish 31.0669,-91.9668  Avoyelles Parish    31.0669 , -91.9668    \n",
       "3        Cottonport 30.9841,-92.0535        Cottonport    30.9841 , -92.0535    \n",
       "4         Pineville 31.3224,-92.4343         Pineville    31.3224 , -92.4343    \n",
       "\n",
       "               new_coord                            mlg_key  \\\n",
       "0     31.2928 , -92.4592        Alexandria 31.3113,-92.4451   \n",
       "1          31.2 , -92.54    Rapides Parish 31.1669,-92.4835   \n",
       "2          31.07 , -92.0  Avoyelles Parish 31.0669,-91.9668   \n",
       "3   30.98408 , -92.05346        Cottonport 30.9841,-92.0535   \n",
       "4     31.3224 , -92.4343         Pineville 31.3224,-92.4343   \n",
       "\n",
       "                           geocoded_coordinates_list  median_error_distance  \n",
       "0  [(31.31129, -92.44514), (31.31129, -92.44514),...             671.350717  \n",
       "1                                                 []                    NaN  \n",
       "2                                                 []                    NaN  \n",
       "3                            [(30.98408, -92.05346)]               0.000000  \n",
       "4                             [(42.15536, -75.1635)]            1947.840891  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lgl_camcoder = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/lgl_geocoded_results_CamCoder.csv')\n",
    "df_lgl_camcoder.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "1ed8305f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVkAAAF3CAYAAAAPVDidAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAk4ElEQVR4nO3debQcVbn38e/PhCFhDEYUkiCJgBJ4UTAyeYGgqBHvguvwMogDggRFQK4Cl0G5gHepTCoIKgEB8UUGETUXoyBckAsIJKIMCQSTECEBIUFAZhJ43j/2Pkmn6XNOddK7T7r5fdaqdap27a5+qs7p51Tv2rVLEYGZmZXxhoEOwMysmznJmpkV5CRrZlaQk6yZWUFOsmZmBTnJmpkV5CRr1kUk7S8pJI0f6FgscZJ9HZA0Pn/wjqxYf3VJh0j6H0kLJC2S9JSkqZJOkfSOuvo9H+xPlNmDXuMcJumEHNdTkl6WNE/SLyR9TJLaGU8vMQ6VdISk/5X0j3wsH5M0JR+3wQMdo5XlX7AtQ9IY4Gpgc+APwHeBR4E1gXcBBwBHStooIuYPYJzbAr8G1gcmA5cA/wRGALsDvwC+BPxgAGPcBPgNsBlwHfAtYCEp5t2AC4GxwNEDFaOV5yRrS0gaQkoKbwM+FhG/bFBndeDfgQG7VVDSW4D/BlYHdomIm+uqfEPSh4BhbQ8uy8fyamAM8PGIuKquyimS3gO8p+3BNUnSIGC1iHh+oGPpRG4usFqfB94BnNYowQJExIsR8a2IeKS9oS3jKNLZ4H80SLAARMQ1EXFZz7KkvSVNlvSQpJckLZT0K0lb1b9W0lxJN0p6p6TrJD0r6XFJZ0ganJtTTpc0X9KLkm6StHndZj4PvB04o0GC7YlxakQsc6Yt6d8k3SLpufy+t0jas9HrJR0k6f68P7MkHQE0bCKRtE5u6pmV6y+QdGn+5lJbr6fpZzdJX5c0G3gR2KvRdq1/PpO1Wj1tqucPaBT9+zjwMvCTJl5zKPAEMAn4O+lsfSJwi6RtIuKvdfVHAr8HLgeuBD4IfAVYDGwBDAG+DQwHjgR+JWnziHg1v77nWE6qGqCkQ4BzgPuBk3Px/nnbB0fEpJq6R5Cacu4CjgOG5jgeb7DddYBbgY2AC4DpwAbAIcDtksZFxN/qXnY6sApwHqkZZmbV/bA6EeGpyydgPOnr/ZH91HsCeLpB+SBSMqmdhtSs3z9v/xNt2Je18nvd3eTr1mhQtjnwEvCDuvK5+T3+b135n4BXSW3Bqik/PNf/UH/Hso/4hgHPArOAtWvK1wZmA88A6+aydYHngBnA0Jq6I/M2AhhfU34m8ALwzrr3fCspgV7U4Hc5s3bbnpZ/cnOB1Vqb9KGrtzmwoG76UhvjqrV2/tkozl5FxHMAStaWNJy0HzOB7Rq8ZH5E/Lyu7GbS1/HvR85I2f/mn5vWxflMEyF+AFgDOCsiluxbnj+LdOFxt1z8QdKZ6zlR004aEfNIFwCXyD0s9gNuAuZLGt4zkRL1bXl79X4YboNtCTcXWK1/sjSJ1XqQlAQA3kn6KtkSktYkJZBa/4iIl3t5SU8CWqvJ99ka+AbprH6NutUPNnhJo7Ine1nXU/7GmrJ/Nhnj6PxzeoN1PWVj6n7e36DujLrlN+W4Pkj6p9LIqw3KHuilrjXJSdZq3QvsLGl0RCxJJPks8DoASYtb/J5HAv9ZV7YrcGOjyhHxjKS/Ae+QNCQiXujvDSRtRDqT+ycp0c4kncUF8D1em+QBXuljk72tq73o1HMsx0TEnP5iLKgnpuuAU5p4nc9iW8RJ1mpdCexMujJ+fJve82LS1/Bad/XzmqtI3cg+TbULSx8lJdI9IuKG2hWS3khql221X7D0WB5XoX5PIt4CuL5u3di6Oj0/39FH3R4LgKdI7bzXVYjDWsxtslbrfNJX0KMkfbSXOi29iyoi5kTEdXXTk/287FRS8jhV0g4Ng5Q+KGmfvNhz5qm6OgcBb1mR+PtwPumM+cg+umC9O/cogNST4TngMElr1dRZCziMdEHr9zV1XwC+JGloTd2RwCdr3yNSb4dLgG17uyNP0vrN755V5TPZ15f355sJ6i2MiB9FxAuSPkLqRH+VpBuBa0ldntYmnTntTUpaDzfYzsfrb7nNZkfEpS3ZAyAi/i7pX0lX+W+W9CuWNgdsCEwA/gX4Yn7Jb0lff38q6WxSG+p7SXeGzabA5yAins8x/obUBetaUnJ8gtROuivwIdI/DCLiKUlHk7pw3S7poryp/YFNgIMj4ulc90lJXye1jd8q6WLShbAvAH8Ftq4L5/i8v1dIuoJ0setlUu+C3Um9JvZv8SGwHgPdvcFT+YmlXbh6m+6vqz+E1HvgBtJtoItIXzmnkZLC2+vq79/P9n9XaL/WI7XnTgOeJiWOeaRmjz3q6u5MapZ4Ju/Lb4AtSW2/c+vqzgVubPB+J+b92biufONcfmKD1wwlNW3cTErui4DH8vt/GhhUV/+jpD6tz+XpVuDfetn/g0lnyy+Run4dAXyOui5cNXF8HbiHdBb8DHAfqR/sdg1+l+Mbvaen5iflA2tmZgW4TdbMrKC2JVlJF+T7v+/tZb0knZXvrb5b0jbtis3MrJR2nsleRLog0ZsPk+6Y2ZR0T/kP2xCTmVlRbUuyEXET8I8+quwJXBzJbcC6kjZoT3RmZmWsTF24RrBst6B5uezR+oqSJpLOdhk7duy7p09vdCeimVnb9Np/vCMvfEXEpIgYFxHjhgwZMtDhmJn1amVKsvOBUTXLI3OZmVnHWpmS7GTgM7mXwfaksThf01RgZtZJ2tYmK+lS0p1HwyXNI92pswpARPwImEK6xW8W6RbIz7UrNjOzUtqWZCNi337WBwM3ELSZWRErU3OBmVnXcZI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgtqaZCVNkDRT0ixJxzRYv5GkGyT9WdLdknZvZ3xmZq2miGjPG0mDgAeADwDzgKnAvhExo6bOJODPEfFDSWOBKRGxcV/bXXXVVWPRokVNx7PhyFHMf/ihpl9nZtaAelsxuI1BbAvMiog5AJIuA/YEZtTUCWDtPL8O8Eh/G120aBF7n3tr08FcfvCOTb/GzKxZ7UyyI4CHa5bnAdvV1TkRuFbSYcAawG7tCc3MrIyV7cLXvsBFETES2B34qaTXxChpoqRpkqa1PUIzsya0M8nOB0bVLI/MZbUOBK4AiIg/AqsDw+s3FBGTImJcRIwrFKuZWUu0M8lOBTaVNFrSqsA+wOS6Og8B7weQtDkpyS5oY4xmZi3VtiQbEYuBQ4FrgPuAKyJiuqSTJe2Rq30VOEjSXcClwP7Rru4PZmYFtPPCFxExBZhSV3ZCzfwM4L3tjMnMrKSV7cKXmVlXcZI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKcpI1MyuoUpKV1Nbbb83MukXVM9lHJZ2eR8YyM7OKqibZ44AdgXsl/VHSgZLWLBiXmVlXqJRkI+K8iNgR2BK4Gfgv0tntBZI8apaZWS+auvAVEfdFxFGkpxocB3wSuEnS/ZK+0OhRMWZmr2dNXdDKTzT4GHAA8D7SWe2PgQ2BrwPjSU88MDMzKiZZSduQEuu+wCLgYuDQiHigps7VgB9saGZWo+qZ7FTgWmAi8Ov8KJl6c4HLWhSXmVlXqJpkx0TE3/qqEBHPAZ9b8ZDMzLpH1QtVG0vapb5Q0i6Sdm5xTGZmXaNqkv0uMKxB+dp5nZmZNVA1yb4duKtB+b15nZmZNVA1yb4AbNCgfATwcuvCMTPrLlWT7DXAKZKWNBlIWg/4Vl5nZmYNVO1dcCRwEzBX0t25bCvgcWDvEoGZmXWDSkk2Ih6V9E5gP+BdufgnwM8i4vlCsZmZdbzKt9XmZHpewVjMzLpO5SQraSSwM7A+dW25EfGdFsdlZtYVqo5dsB9wAbAYWABEzeoAnGTNzBqoeiZ7MnAG8PWIeKVgPGZmXaVqF643A+c7wZqZNadqkp0CbFcyEDOzblS1ueD3pJsRtgDuIY0pu0REXNXqwMzMukHVJHtu/nlcg3UBDGpNOGZm3aXqzQh+dpeZ2XJw8jQzK6hSklVyiKTpkp6XNCaXHyNpr7Ihmpl1rqpnsl8GvgZMAlRTPh84tNVBmZl1i6pJ9gvAQRFxJumurx53Alu0PCozsy5RNcm+lfQUhHqLgCGtC8fMrLtUTbJzgG0alO8OzGhdOGZm3aVqP9nTgbMlDSW1ye4g6dPA0cABpYIzM+t0VfvJXihpMPBNYCjwU+AR4PCIuLxgfGZmHa2ZQbvPA86TNBx4Q0Q8Xi4sM7PuUDnJ9oiIhSUCMTPrRlUH7b6HZQfqXkZEbNWyiMzMukjVM9kr65ZXIT1Q8b3AOa0MyMysm1S98HVSo3JJR5H60JqZWQMrOkDMVaTHhJuZWQMrmmR3Bp5vRSBmZt2o6oWvyfVFwAbA1kDDpgQzM6t+4euJuuVXgenAcRFxbWtDMjPrHlUvfH2udCBmZt2orU9GkDRB0kxJsyQd00udvSTNyAOE/6yd8ZmZtVrVNtkH6eNmhFoRMaaXbQwi9an9ADAPmCppckTMqKmzKXAs8N6IeFLS+lXe08xsZVW1TfYi4CvAHcAfc9kOwLbAd4BnK2xjW2BWRMwBkHQZsCfLDpV4EHBORDwJ4PERzKzTVU2yo4FTIuKbtYWSjgW2iIhPVdjGCODhmuV5wHZ1dTbL272F9JjxEyPid/UbkjQRmFgxdjOzAVM1yX6MxoN2/5z09b5VBgObAuOBkcBNkv5PRDxVWykiJpGeN4akSs0YZmYDoeqFr+dIia/eeKrfjDAfGFWzPDKX1ZoHTI6IRRHxIPAAKemamXWkqmey3wXOkTQOuC2XbQ98Fjix4jamAptKGk1KrvsAn6yr8ytgX+DCPG7tZqRH35iZdaSq/WRPlTSX9GjwvXLxfcBnI+KKittYLOlQ4BpSe+sFETFd0snAtIiYnNd9UNIM4BXgqIiovxHCzKxjKKKzmzQlxd7n3tr06y4/eEc6fd/NbKWh3lZUvhlB0uqSPiHpPyStm8veJmm9FgRoZtaVqt6MsAlwHbAmsC6pV8FTwBfz8ueLRGdm1uGqnsl+D7gWeDPwQk35ZGDXFsdkZtY1qvYu2BHYPiJekZZpengI2LDlUZmZdYlmBohZpUHZRsDTLYrFzKzrVE2y15LGLugRktYmDdj9m5ZHZWbWJao2F3wV+B9JM4HVgcuBTYDHWNpv1szM6lS9GWG+pHeR7sbahnQGPAm4JCJe6Ou1ZmavZ/0mWUmrkEbPen9EXABcUDwqM7Mu0W+bbEQsAhZRcdBuMzNbquqFr+8Dx0qq2oZrZmZUv/C1E7ALMF/SvaShD5eIiD1aHZiZWTeommQXAr8oGYiZWTfqM8lK2gP4rR8Jbma2fPprk/0lMKxnQdI9kkb1Ud/MzGr0l2Trx0jcmMa315qZWQPNjF1gZmZN6i/JBsv2j61fNjOzPvTXu0DAHyQtzstDgd9Kerm2UkRsVSI4M7NO11+SPalu2d24zMya0GeSjYj6JGtmZk3whS8zs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4J67V0g6TNVNxIRF7cmHDOz7tJXF65z6pZXJY1b8GpefgPpiQkvAU6yZmYN9NpcEBFr9UzAPsDdpMG7V8/TTsBfgE+2IU4zs45UtU32dODwiLglIhbn6RbgCOCMYtGZmXW4qkl2Y+oeOZM9D2zUsmjMzLpM1SR7O3CWpBE9BXn+u8BtJQIzM+sGVZPsgcAbgbmS5kqaC8wF1gcOKhOamVnnq/QgxYiYLWkr4APAO3LxfcB1EeHxZc3MelH1abXkZHptnszMrILKd3xJOkTSdEnPSxqTy46RtFe58MzMOlulJCvpCOBrwCSWfbjifODQ1odlZtYdqp7JfgE4KCLOBBbXlN8JbNHyqMzMukTVJPtW4N4G5YuAIa0Lx8ysu1RNsnOAbRqU7w7MaF04ZmbdpWrvgtOBsyUNJbXJ7iDp08DRwAGlgjMz63RV+8leKGkw8E3SY8F/CjxCGs/g8oLxmZl1tGb6yZ4HnCdpOPCGiHi8XFhmZt2haheuEyTtChARC3sSrKQ1JJ1QMkAzs05W9cLXicDvJB1WV74m8J8tjcjMrIs084yvA4GTJP1A0qBSAZmZdZNmkuy1wA7A+0hnteuUCcnMrHtUTbIBEBEzge1z2R3AZiWCMjPrFlWT7JLxCiLiKWAC8HvgdwViMjPrGlW7cJ0EPNuzEBGvAIdK+hOwc4nAzMy6QdWbEU7qpfxC4MKWRmRm1kV6TbKSvgL8ICJezPO9iYj4butDMzPrfH2dyR4G/AR4Mc/3JkgPVOyXpAnAmcAg4PyI+HYv9T4OXAm8JyKmVdm2mdnKqNckGxGjG80vr9y39hzSc8LmAVMlTY6IGXX11gK+THpCrplZR2umn+yK2haYFRFzIuJl4DJgzwb1vgGcQjqDNjPraP21yVYSEd+pUG0E8HDN8jxgu7r33AYYFRG/kXRUH7FNBCZWjc/MbKD01yZbRQBVkmyfJL0hb2f/ft8wYhLpeWNI8iPJzWylValNtkXmA6Nqlkfmsh5rAVsCN0oCeAswWdIevvhlZp2qnW2yU4FNJY2WtCqwDzC5Z2VEPB0RwyNi44jYGLgNcII1s45WedBuScOADwMbAavWrouIk/t7fUQslnQocA2pC9cFETFd0snAtIiY3PcWzMw6T6UkK2l74DfAS8CbSF/zN8jLc4F+kyxAREwBptSVNRz0OyLGV9mmmdnKrGpzwWnAJaQeAi+ShjvcCJhG6m5lZmYNVE2yWwFnR0QArwCrRcRjwH+QnppgZmYNVE2yL9fMPwa8Nc8/C2zY0ojMzLpI1QtfdwLvAR4AbgT+S9KbgU8Bd5cJzcys81U9kz0eeCTPfw1YAHwfGIbvvDIz61XV8WSn1cwvIHXlMjOzfrTzZgQzs9edqv1kh5F6EewKrE9dco6I9VsemZlZF6h64etiYAvSIN6PkZ9ea2ZmfauaZMcDu0TEnQVjMTPrOlXbZGc3UdfMzLKqifPLwLckvTM/RsbMzCqo2lwwCxhCuimBPN7rEhHhxGtm1kDVJHspsA5wOL7wZWZWWdUkOw7YNiLuLRmMmVm3qdomOwNYu2QgZmbdqGqS/RrwHUm7SXqzpPVqp5IBmpl1sqrNBT1PM7iWZdtjlZd94cvMrIGqSXbXolGYmXWpfpOspFWAU4HPRMTM8iGZmXWPfttkI2IRMBp32zIza1rVC18/AQ4qGYiZWTeq2ia7BrCfpA8AfwKeq10ZEYe3OjAzs25QNcluTr6lFhhTt87NCGZmvaj6+Bn3LjAzWw5Vz2QBkLQ6sAnp7HV2RLxYJCozsy5R6cKXpFUknQY8CdwF3AM8KenU3MXLzMwaqHomewqwL/AF4OZcthPwLVKiPrL1oZmZdb6qSfaTwAERMaWmbLakBcD5OMmamTVUtZ/sOqRH0NSbDazbsmjMzLpM1SR7F2nA7npfBv7SsmjMzLpM1eaCo4EpknYDbstl2wMbAh8uEZiZWTeodCYbETcBmwFXAmvm6efA2yPi5r5ea2b2ela5n2xEPAIcXzAWM7Ou02eSrfrUg4j4R2vCMTPrLv2dyS6k/7EJosJ2zMxel/pLjn2NWTCB1LtgcevCMTPrLn0m2Yj4Q32ZpK2B00h3fJ0LfKNMaGZmna9qP1kkjZb0M+AO4AlgbEQcHhELikVnZtbh+k2ykt4o6UzgfuAtwI4RsXdENLoDzMzMavSZZCUdT7p1dhdgz4h4X0RMbUtkZmZdoL8LX98AXgDmAYdIOqRRpYjYo9WBmZl1g/6S7MX48TJmZsutv94F+7cpDjOzrlS5d4GZmTXPSdbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgpxkzcwKamuSlTRB0kxJsyQd02D9VyTNkHS3pOslvbWd8ZmZtVrbkqykQcA5pKfbjgX2lTS2rtqfgXERsRXpoY2ntis+M7MS2nkmuy0wKyLmRMTLwGXAnrUVIuKGiHg+L94GjGxjfGZmLdfOJDsCeLhmeV4u682BwG+LRmRmVthK+QBESZ8CxpHGsW20fiIwsa1BmZkth3Ym2fnAqJrlkblsGZJ2A44HdomIlxptKCImAZNyfQ/FaGYrrXY2F0wFNs3PClsV2AeYXFshP6TxXGCPiHi8jbGZmRXRtiQbEYuBQ4FrgPuAKyJiuqSTJfU8WeE0YE3g55L+ImlyL5szM+sIbW2TjYgpwJS6shNq5ndrZzxmZqX5ji8zs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMryEnWzKwgJ1kzs4KcZM3MCnKSNTMr6PWbZN8wGElNTyNGbTTQkZtZBxk80AEMmFcXs/e5tzb9sssP3rFAMGbWrV6/Z7JmZm3gJGtmVpCTrJlZQU6yZmYFOcmamRXkJGtmVpCTrJlZQU6yZmYFOcmamRXkJGtmVpCTrJlZQU6yZmYFOck2y6N3mVkTXr+jcC0vj95lZk3wmayZWUFOsmZmBTnJmpkV5CRrZlaQk6yZWUFOsu3SIV2/RozaqCPiNOsUbe3CJWkCcCYwCDg/Ir5dt3414GLg3cATwN4RMbedMRazvF2/vrgzkpp+3YYjRzH/4Yeaft0j8x7uiDjNOkXbkqykQcA5wAeAecBUSZMjYkZNtQOBJyNiE0n7AKcAe7crxpVSp/TLbfM/kUGrrMYri15q+nVO6q9PI0ZtxCPzHl6u167o30w7z2S3BWZFxBwASZcBewK1SXZP4MQ8fyVwtiRFRLQxzu6QmydWeivwT6QTknq7/xksbzJp9/u1+7gs7zc0WPETFrUrf0n6BDAhIj6flz8NbBcRh9bUuTfXmZeXZ+c6C+u2NRGYmBe3BO5twy70ZTiwsN9a3R8DrBxxOAbH0O4YFkbEhEYrOvK22oiYBEwCkDQtIsYNZDyOYeWKwzE4hpUphnb2LpgPjKpZHpnLGtaRNBhYh3QBzMysI7UzyU4FNpU0WtKqwD7A5Lo6k4HP5vlPAP/j9lgz62Rtay6IiMWSDgWuIXXhuiAipks6GZgWEZOBHwM/lTQL+AcpEfdnUrGgq3MMS60McTiGxDEkAxpD2y58mZm9HvmOLzOzgpxkzcwK6ugkK2mCpJmSZkk6poXbHSXpBkkzJE2X9OVcvp6k30v6a/45LJdL0lk5jrslbVOzrc/m+n+V9Nne3rOPWAZJ+rOkq/PyaEm35/e6PF9ERNJqeXlWXr9xzTaOzeUzJX1oOWJYV9KVku6XdJ+kHdp9LCT9e/5d3CvpUkmrlz4Wki6Q9Hjuv91T1rL9lvRuSffk15wlvfYuiV5iOC3/Lu6W9EtJ6/a3f719Vno7hlXiqFn3VUkhaXi7j0UuPywfj+mSTi19LJoWER05kS6ezQbGAKsCdwFjW7TtDYBt8vxawAPAWOBU4JhcfgxwSp7fHfgtIGB74PZcvh4wJ/8clueHNRnLV4CfAVfn5SuAffL8j4Av5vlDgB/l+X2Ay/P82HxsVgNG52M2qMkYfgJ8Ps+vCqzbzmMBjAAeBIbUHIP9Sx8LYGdgG+DemrKW7TdwR66r/NoPV4zhg8DgPH9KTQwN948+Piu9HcMqceTyUaSL2X8Dhg/AsdgVuA5YLS+vX/pYNJ1PWrGRgZiAHYBrapaPBY4t9F6/Jo25MBPYIJdtAMzM8+cC+9bUn5nX7wucW1O+TL0K7zsSuB54H3B1/gNcWPMBW3IM8h/6Dnl+cK6n+uNSW69iDOuQEpzqytt2LEhJ9uH84Rycj8WH2nEsgI3rPtQt2e+87v6a8mXq9RVD3bqPApc0+gz07B+9fFb6+nuqGgfp9vd3AnNZmmTbdixIiXG3BvWKHotmpk5uLuj54PWYl8taKn/V3Bq4HXhzRDyaV/0deHM/saxojN8DjgZezctvBJ6KiMUNtrfkvfL6p3P9FY1hNLAAuFCp2eJ8SWvQxmMREfOB04GHgEdJ+/Yn2n8soHX7PSLPr0gsAAeQzvyWJ4a+/p76JWlPYH5E3FW3qp3HYjNgp/w1/w+S3rOcMazQsehLJyfZ4iStCfwCOCIi/lm7LtK/u2L93yT9K/B4RPyp1HtUNJj0Fe2HEbE18Bzpa/ISbTgWw0iDB40GNgTWABreJ95Opfe7P5KOBxYDlwzAew8FjgNOaPd71xlM+oazPXAUcEWj9tyB1MlJtsptustN0iqkBHtJRFyVix+TtEFevwHweD+xrEiM7wX2kDQXuIzUZHAmsK7SLcf12+vtluQVPU7zgHkRcXtevpKUdNt5LHYDHoyIBRGxCLiKdHzafSygdfs9P88vVyyS9gf+FdgvJ/vlieEJej+G/Xkb6Z/eXflvdCRwp6S3LEccK3Is5gFXRXIH6Vvf8OWIYUWORd9a0eYwEBPpP9gc0i+6pwF7ixZtW6TBw79XV34ay170ODXPf4RlG/rvyOXrkdozh+XpQWC95YhnPEsvfP2cZRvnD8nzX2LZiz1X5PktWPYCwByav/D1v8Db8/yJ+Ti07VgA2wHTgaF5uz8BDmvHseC1bYAt229ee7Fn94oxTCANEfqmunoN948+Piu9HcMqcdStm8vSNtl2HosvACfn+c1ITQEqfSya+vy0YiMDNZGuYj5Aulp4fAu3+y+kr4F3A3/J0+6kdpvrgb+Srmj2/IGINCD5bOAeYFzNtg4AZuXpc8sZz3iWJtkx+Q9yVv6j6LmqunpenpXXj6l5/fE5tpk0uGpb4f3fBUzLx+NX+QPS1mMBnATcTxrW8qf5w1P0WACXktqAF5HOmA5s5X4D4/L+zAbOpu7iYh8xzCIlk56/zR/1t3/08lnp7RhWiaNu/VyWJtl2HotVgf+XX3sn8L7Sx6LZybfVmpkV1MltsmZmKz0nWTOzgpxkzcwKcpI1MyvISdbMrCAnWVupSTpb0o0DHUc9ScMkPSbpbXl5fO1IVC18ny9J+u9WbtPay0nWipF0UU48P26w7pS87uoC73tj3nb9dFlNndry5yXNkfQzSTtVfJvjgCkRMbvV8dc5H3h3E3HZSsZJ1kp7GNgrDyoDLLnV9TOkAV9KuZA0ulPtdHBdnYNy+eakju0vA3+QdFRfG8737X+e9Ey6oiLiJdJQl4eXfi8rw0nWSrubdHfUXjVlHwFeBG6srag0QPnpkp7M0/dIt0LW1rlR0g8kfVPSwjyI8+mS6v+Wn4+Iv9dNT9fVeSqX/y0iboiI/YFvA9+StEkf+7Q76Y7AW3qroDRw+C8l3SlpfUkb57PmffJoUS/kUc22krSlpFslPSfpZkmj6zY3mTSOxdA+YrKVlJOstcOPSbdT9jiAdKZZf7vhV0lnlweTxvMcBOzXYHv7kUaf2hE4FDgC2LtFsZ5B+lz8Wx91dgL+FL3cLilpbeB3pHv1x0fE4zWrTyINtL018BTpVtHvk24B3ZZ0S/BZdZucRrrnfofmdsVWBk6y1g4/A8ZJ2jSP0jQBuKhBvSNIA65cERH3A18mjdlab0ZEnBARD0TEFcANwPvr6kyU9GzddEh/gUbEE6SRtcb0Ue2twCO9rFs/x/MM8KGoGyIT+E5ETMn7dwZpBP/v5zPp6aT79neti+l50pi4G/cXv618BvdfxWzFRMSTkn5JOoN9CrgxIh6qHfZT0jqk9tE/1rzuVUm3s+zQdJCaIGo9QkputS4nnTXWWlAxZNH3OLFDgMd6WXcNaaCSj8fSAaBr1cbes4176srWkDQ0J9ceL+T3tQ7jJGvtcgFpiMJnWfGBnhfVLQev/Vb2dETManbDuQvWm0jD4fVmIWkkskauJrU/b0kaIatebezRR1n9/qxH9X8SthJxc4G1y/Wkq/fDScMlLiNflHqUNKYokJ56SmqnbKevkgZ+/lUfdf5M+prfyNdJY5FeL+ldrQgo98VdnXSGbB3GZ7LWFhERkrYijRP6Ui/VzgSOlfQA6Sv0IaQmhEd7qd+Xobn9t9bLEfGPmuV1c51VSSP9f5bUtezofvq/XgOcIumNuQ13GRFxfP4HcZ2k98drn4HVrJ2AORHx1xXcjg0AJ1lrm4h4pp8qZwBvIXXAhzQ49yWkfqzN+lyeat1CGpC9x3n550ukRH4bqTfATX1tOCLukXQH6akL5/RS57icaK+X9H7ShavltW9NrNZhPGi32XKQNIF05j02Il4p+D5bkppaNmvQz9c6gNtkzZZDRPyOdBY7sr+6K2hD4DNOsJ3LZ7JmZgX5TNbMrCAnWTOzgpxkzcwKcpI1MyvISdbMrCAnWTOzgv4/lrEJF6LmeIMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_lgl_camcoder, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('LGL - CamCoder', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "38dcb857",
   "metadata": {
    "scrolled": true
   },
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>wikipedia_name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>page</th>\n",
       "      <th>geovirus_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Pandi</td>\n",
       "      <td>Pandi, Bulacan</td>\n",
       "      <td>14.87</td>\n",
       "      <td>120.95</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Pandi,_Bulacan</td>\n",
       "      <td>Pandi, Bulacan 14.87,120.95</td>\n",
       "      <td>Pandi, Bulacan</td>\n",
       "      <td>14.87 , 120.95</td>\n",
       "      <td>14.866666666667 , 120.95</td>\n",
       "      <td>Pandi, Bulacan 14.87,120.95</td>\n",
       "      <td>[(14.88722, 120.85722), (14.88722, 120.85722),...</td>\n",
       "      <td>10.229429</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>15.47</td>\n",
       "      <td>120.75</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Central_Luzon</td>\n",
       "      <td>Central Luzon 15.47,120.75</td>\n",
       "      <td>Central Luzon</td>\n",
       "      <td>15.47 , 120.75</td>\n",
       "      <td>15.466666666667 , 120.75</td>\n",
       "      <td>Central Luzon 15.47,120.75</td>\n",
       "      <td>[(15.30714, 120.79762)]</td>\n",
       "      <td>18.458642</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>15.00</td>\n",
       "      <td>121.08</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Bulacan</td>\n",
       "      <td>Bulacan 15,121.08</td>\n",
       "      <td>Bulacan</td>\n",
       "      <td>15 , 121.08</td>\n",
       "      <td>15.0 , 121.08333333333</td>\n",
       "      <td>Bulacan 15,121.08</td>\n",
       "      <td>[(15.0, 121.08333), (15.0, 121.08333), (15.0, ...</td>\n",
       "      <td>0.000358</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Manaoag</td>\n",
       "      <td>Manaoag, Pangasinan</td>\n",
       "      <td>16.04</td>\n",
       "      <td>120.48</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Manaoag,_Pangasinan</td>\n",
       "      <td>Manaoag, Pangasinan 16.04,120.48</td>\n",
       "      <td>Manaoag, Pangasinan</td>\n",
       "      <td>16.04 , 120.48</td>\n",
       "      <td>16.043889 , 120.485556</td>\n",
       "      <td>Manaoag, Pangasinan 16.04,120.48</td>\n",
       "      <td>[(16.03333, 120.48333), (16.03333, 120.48333)]</td>\n",
       "      <td>1.197965</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>15.92</td>\n",
       "      <td>120.33</td>\n",
       "      <td>https://en.wikipedia.org/wiki/Pangasinan</td>\n",
       "      <td>Pangasinan 15.92,120.33</td>\n",
       "      <td>Pangasinan</td>\n",
       "      <td>15.92 , 120.33</td>\n",
       "      <td>15.916666666667 , 120.33333333333</td>\n",
       "      <td>Pangasinan 15.92,120.33</td>\n",
       "      <td>[(15.91667, 120.33333), (15.91667, 120.33333)]</td>\n",
       "      <td>0.000514</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0           name       wikipedia_name    lat     lon  \\\n",
       "0           0          Pandi       Pandi, Bulacan  14.87  120.95   \n",
       "1           1  Central Luzon        Central Luzon  15.47  120.75   \n",
       "2           2        Bulacan              Bulacan  15.00  121.08   \n",
       "3           3        Manaoag  Manaoag, Pangasinan  16.04  120.48   \n",
       "4           4     Pangasinan           Pangasinan  15.92  120.33   \n",
       "\n",
       "                                                page  \\\n",
       "0       https://en.wikipedia.org/wiki/Pandi,_Bulacan   \n",
       "1        https://en.wikipedia.org/wiki/Central_Luzon   \n",
       "2              https://en.wikipedia.org/wiki/Bulacan   \n",
       "3  https://en.wikipedia.org/wiki/Manaoag,_Pangasinan   \n",
       "4           https://en.wikipedia.org/wiki/Pangasinan   \n",
       "\n",
       "                       geovirus_key               toponym         old_coord  \\\n",
       "0       Pandi, Bulacan 14.87,120.95       Pandi, Bulacan    14.87 , 120.95    \n",
       "1        Central Luzon 15.47,120.75        Central Luzon    15.47 , 120.75    \n",
       "2                 Bulacan 15,121.08              Bulacan       15 , 121.08    \n",
       "3  Manaoag, Pangasinan 16.04,120.48  Manaoag, Pangasinan    16.04 , 120.48    \n",
       "4           Pangasinan 15.92,120.33           Pangasinan    15.92 , 120.33    \n",
       "\n",
       "                            new_coord                           mlg_key  \\\n",
       "0            14.866666666667 , 120.95       Pandi, Bulacan 14.87,120.95   \n",
       "1            15.466666666667 , 120.75        Central Luzon 15.47,120.75   \n",
       "2              15.0 , 121.08333333333                 Bulacan 15,121.08   \n",
       "3              16.043889 , 120.485556  Manaoag, Pangasinan 16.04,120.48   \n",
       "4   15.916666666667 , 120.33333333333           Pangasinan 15.92,120.33   \n",
       "\n",
       "                           geocoded_coordinates_list  median_error_distance  \n",
       "0  [(14.88722, 120.85722), (14.88722, 120.85722),...              10.229429  \n",
       "1                            [(15.30714, 120.79762)]              18.458642  \n",
       "2  [(15.0, 121.08333), (15.0, 121.08333), (15.0, ...               0.000358  \n",
       "3     [(16.03333, 120.48333), (16.03333, 120.48333)]               1.197965  \n",
       "4     [(15.91667, 120.33333), (15.91667, 120.33333)]               0.000514  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_geovirus_camcoder = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/geovirus_geocoded_results_CamCoder.csv')\n",
    "df_geovirus_camcoder.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "6deb4a33",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAF3CAYAAACSWTEVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAnVUlEQVR4nO3debgcVZnH8e+PhLAvYRVIkERACQwKRjYHiIBOxHlgdBgWGXYIiCiMDMiiDOC4sIuCSlgFQYgMakQUBhQZF5aAsgQIhBAhgSGEbZQtCb7zxzkXOpXue+ve20sl/fs8Tz+36/Sp6reru997+tSpU4oIzMys85bqdABmZpY4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7K1jKRxkkLSgZ2OxVpL0kxJt3c6jsWdE/JiQNKyko6U9CtJz0uaL+llSfdIOkPS+9oQw49ycv1AL3Uk6ckc23KtjqlTJO0saZKkpyW9Kekvku6V9FVJIzodH4CksZKukDRD0uuSXpX0kKRvtuPzYgMjnxhSbZJGAzcCmwC/AW4BngVWBD4A7AasBqwfEbNbGMd44BfAtyLi6AZ1dgJuAy6KiCMkLQUMA+ZHxFutiq1d8uu5CDgU+DNwHfA46TV+EPhnYF5ErNWxIAFJ/wH8BzAXuAZ4mNT42hTYA1gTGB4Rf2nic84EZkbEuGZtsytFhG8VvQHLAY8A84BPNqizLHAisG6LY1kKeIr0JR/WoM5VQAAf6ue2BazY6f1dIs7T8+u7pt4+AFYFzutwjAfnGH8FrNLgM/UNYOUmP+9M4PYmb3OlTr/nbX//Oh2Ab728OfC5/OX66gDWXQU4A5gOvAk8D/wQGF2n7hrAhcDTOfk/nZdXL9TrSUh71NnGysBrwIM1ZeNy/QPrlQGfJbXe3gROzY/X/WI32NaywKnAtPzcLwMPAme14L1YC3g9x7dsyXXWBc4B/gS8BLyRX+8XgSGFugfm17czcAqpBf46cBewTa6zI/Bb4FXSr6QvF7YxLJf/BVizH6+t1Puf644EJgGvAP8H/Ax4Ty/v2y6kX3Uv59f/AHBEnXozgduBLYCb8/af7PR3sN23oViV7ZH/XtKflSStAvweWB+4DJgKrAMcCdwlaWxE/LlQd8Nc9z7Sl+IzwE6Stop3ftpeDnwJOAi4vvC0e5NaX5eWDPMYYHXgYuB/SUmgvy4ktQivBM4FhgIbATsNYFt9+QTpH8CVEfFGyXU2Bz4F/Bh4AlgaGE9qoY4GDq+zzjeAIcD5pAR7LHCLpP1J+3YicDWwJ3C6pCcj4gd53Q8D7wKuiojnywTYn/df0qrAHaSk/D3SP5cdgV+T3vvitifkencCXyX9I/ko8F1J74mI4wqrrE9q2f8I+C9St1x36fR/BN8a34AXgFfqlA8htWpqb8vVPH4+qXX1/sJ67ya1aq6oKfsqqWV2ZKHuZ3P5VwrltwELgHUK5X8gtXTXqCkbR+MW8ovAWnVe20zKt5BfBG5q03txTn7+T/VjneXIx2kK5VcBb9XuQ95pId9HTXcI6RhBAPOBsTXlPa3hP9SU9fyi+kI/Yiz9/gNfy2UHFep+M5ffXlO2DqlFfE2d5zw/v/7RNWUz8zYObcf7WdWbR1lU28qkBFq0CakLovb2WUgjHYB9SS2Z2ZLW6LmRWih3Ah+r2dYn8/oTC89xUS7/ZKH8UtI/hP17CvJR+22AyRExt+RruzIi5pSs28grwKaSNhvkdspYOf+t937UFRGvR842koZJWi2/DzeT+uTH1lntuxExr2b5f/LfuyJiSs225wF3k34RDDhG+vf+/xPwHOkXSa0z6mx3D2AZ4NLaz2B+/T8jvf5dCuu8SPoV1rXcZVFt/8c7X7JaT5J++gG8Hzi75rE1SV0BHyN9oer5W839UcCUiFhQWyEiFkh6DNiysO4NpP7Ag3jni3hw/ntZoxdSx2P9qNvIMaTW5oOSZpB+Ov8M+FlE/K23FSW9q1A0LyJe7GWVniS3UtngJA0FTiD989qQdPCy1vA6q82oXYiIl9L/WJ6sU/cl0ns94Bjp3/s/GrgnCiNmIuJZSS8XtrtJ/ntrL8+9dmH5ieK2u40TcrU9BOwgaVREvP2FjIhXyR90SQsK6/R86W+lfstlUCLiDUnXAEdK2o500Gk/YBap5VfWa42eokH5Ip/ViPippA2AXUl9mbsAhwD/I2mXQkuz6NnC8m9I3SKNPJT/bkHqEy7jXFI3wnWkroE5pK6HLUnvTb1fqI0SUplEVRtjp/V8Dvdn0X3dY0ZhudFnoms4IVfb9cAOpHGvJ5dc53lSC3bliOitddJjBvBeSUNrW0m5dbcxi35pIHVbHElqJa9GOpD01b5apSW9mLdZNLpe5dyq/QHwg9xd8w3geGB30sGhRj5aWH6pj7h+TuoT3U/SVyPizT7qQ/pHdUdE7F1bKGnDEusOxO9IB0j/SdLqEfFCiXX68/7PADaSNKS2JStpHdKQv1qP579zS34ODZ+pV3WXAI8Cx0kq9uX2WOhncE6KVwNbSdqj7gpS7YkLPyF1cxxaqHZYLl+kNRgR95GGcu3FOwd/+tNd0ZvHgPdJWq8m3mXy81BTNiQf9a+NK4A/5sV6Sb227q2F27191J8DnAVsAFwiaVixjqSVJZ1XU/QWhfdH0grAv/X2XAOVfxGcTOqyuE7SIl0X+azPr0nq6Qr7CeXf/5+Suhn2L9T9Yp1wJpEO8p5W76xNSavk99VquIVcYRHxuqRPkM7UuyHPFXALqRW0MvA+UlJ8i4WHjZ1MGgI1SdIk0oG8eaRRFrsC95KO6gOcCfwLcKGkLUkJbQvST/9p+fF6LgW+TRrGdXtE1GtJD8QFpCF0t0r6Hmk0wX4s+nN2JeBZSZNzzHNI/aGfIbV2f9akeGqdSho9cCjw95KuJY3zHkY6a/JfSPu5J+FeDxwu6TpSF9LapP72Mi3XAYmIyySNJJ2pNz13L/WcqbdJjnEt4Ot5lf68/2cCnwYulvRB0nDKccC2pBOGauOYJekzpEbFI5KuIo2tXhP4O9IBwjGk0RXWo9PDPHzr+0YaPvVZ0kGruaR+yJeBKaQvyXvrrLM88GXSiRKvk04WeIQ07nfrQt01ge+Q+oHn578XUjOErc72h+ftBrBfgzrj6OXEkF62fQApGcwjHcw6njS2+O31SEnw66SRBi+QWmMzSS31jVr8fuxC6g6ZlWP8C+mf3FdYeCjb8qRW9Z9J3R2Pkw7y7VxnvxyYy8bVeb6gZqhiTfkV5B8GdR4bC3w/7783yCftkPq1NyrULf3+k8YKX086gFjmxJAPk1rZc/K+eiZ/jo+l5gSbRut3281zWZiZVYT7kM3MKqJtCVnSZZLmSHqoweOS9C1J0yU9kPuzzMy6RjtbyFeQDgA18nHSWUcbAROA77YhJjOzymhbQo6IO0hjTBvZnXQ6bUTEncCqeXyjmVlXqNKwt/VYeOjWrFy2yFk+eRapCQBjxoz54NSpU9sSoJlZkxRPowcW04N6ETExIsZGxNjllltirxRkZl2mSgl5Nmme1R4jcpmZWVeoUkKeDOyfR1tsQ5oHuNGkJGZmS5y29SFL+iHpLK01JM0indq5NEBEfA+4iXRa73TSWUUHtSs2M7MqaFtCjoh9+ng8KEwgY2bWTarUZWFm1tWckM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCrCCdnMrCKckM3MKsIJ2cysIpyQzcwqwgnZzKwinJDNzCpCaU6fxdewYcNi/vz5A15/3REjmf30U02MyMysT3WvGLLYJ2RJsddFvx/w+tcdvh2L+z4ws8XOknMJJzOzJZETsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRTghm5lVhBOymVlFOCGbmVWEE7KZWUU4IZuZVYQTsplZRbQ1IUsaL2mapOmSTqjz+PqSfi3pj5IekLRrO+MzM+uktiVkSUOAC4GPA2OAfSSNKVT7EjApIrYA9ga+0674zMw6rZ0t5K2A6RExIyLmAdcCuxfqBLByvr8K8Ewb4zMz66ihbXyu9YCna5ZnAVsX6pwK3CLpc8AKwC7tCc3MrPOqdlBvH+CKiBgB7ApcJWmRGCVNkDRF0pS2R2hm1iLtTMizgZE1yyNyWa1DgEkAEfEHYFlgjeKGImJiRIyNiLEtitXMrO3amZDvATaSNErSMNJBu8mFOk8BOwNI2oSUkJ9vY4xmZh3TtoQcEQuAo4CbgUdIoymmSjpd0m652rHAYZLuB34IHBgR0a4Yzcw6qZ0H9YiIm4CbCmWn1Nx/GPhwO2MyM6uKqh3UMzPrWk7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhVRKiFLausp1mZm3ahsC/lZSWfnGdjMzKwFyibkk4DtgIck/UHSIZJWbGFcZmZdp1RCjoiLI2I7YDPgt8B/klrNl0ny7GxmZk3Qr4N6EfFIRBxHutrHScCngTskPSrpiHqXWzIzs3L6dbAuX+njU8DBwE6k1vKlwLrAl4FxpCuBmJlZP5VKyJK2JCXhfYD5wJXAURHxWE2dGwFfdNTMbIDKtpDvAW4BJgA/zZdjKpoJXNukuMzMuk7ZhDw6Iv7cW4WIeBU4aPAhmZl1p7IH4TaQtGOxUNKOknZockxmZl2pbEI+Dxhep3zl/JiZmQ1S2YT8XuD+OuUP5cfMzGyQyibk14F16pSvB8xrXjhmZt2rbEK+GThD0tvdFpJWA76eHzMzs0EqO8ri34E7gJmSHshlmwNzgL1aEZiZWbcplZAj4llJ7wf2BT6Qi78PXBMRr7UoNjOzrlL61OmceC9uYSxmZl2tdEKWNALYAViLQt9zRJzb5LjMzLpO2bks9gUuAxYAzwNR83AATshmZoNUtoV8OnAO8OWIeKuF8ZiZda2yw97WBi5xMjYza52yCfkmYOtWBmJm1u3Kdln8N+nEkE2BB0lzIr8tIm5odmBmZt2mbEK+KP89qc5jAQxpTjhmZt2r7IkhvlaemVmLOdGamVVEqYSs5EhJUyW9Jml0Lj9B0p6tDdHMrDuUbSEfDXwJmAiopnw2cFSzgzIz60ZlE/IRwGERcT7pbL0e9wGbNj0qM7MuVDYhv5t0dZCi+cByzQvHzKx7lU3IM4At65TvCjzcvHDMzLpX2XHIZwMXSFqe1Ie8raT9gOOBg1sVnJlZNyk7DvlySUOBrwHLA1cBzwCfj4jrWhifmVnX6M8E9RcDF0taA1gqIua0Liwzs+5TOiH3iIi5rQjEzKzblZ2g/kEWnpR+IRGxedMiMjPrUmVbyNcXlpcmXez0w8CFzQzIzKxblT2od1q9cknHkcYom5nZIA12cqEbgH2bEYiZWbcbbELeAXitGYGYmXW7sgf1JheLgHWALYC63RlmZtY/ZQ/qvVBY/hswFTgpIm5pbkhmZt2p7EG9g1odiJlZt2vrFUMkjZc0TdJ0SSc0qLOnpIfzZPjXtDM+M7NOKtuH/CS9nBhSKyJGN9jGENKY5Y8Cs4B7JE2OiIdr6mwEnAh8OCJekrRWmec0M1sSlO1DvgL4AnA38Idcti2wFXAu8NcS29gKmB4RMwAkXQvszsLTdx4GXBgRLwF4vgwz6yZlE/Io4IyI+FptoaQTgU0j4l9LbGM94Oma5VnA1oU6G+ft/g4YApwaEb8sbkjSBGBCydjNzBYLZRPyp6g/Qf2PSF0MzTIU2AgYB4wA7pD0dxHxcm2liJhIur4fkkp1pZiZVV3Zg3qvkpJk0TjKnxgyGxhZszwil9WaBUyOiPkR8STwGClBm5kt8cq2kM8DLpQ0Frgzl20DHACcWnIb9wAbSRpFSsR7A58u1PkJsA9weZ53eWPS5aPMzJZ4ZcchnylpJnA0sGcufgQ4ICImldzGAklHATeT+ocvi4ipkk4HpkTE5PzYxyQ9DLwFHBcRxZNSzMyWSIpYvLtgJcVeF/1+wOtfd/h2LO77wMwWO6pXWPrEEEnLStpD0hclrZrL3iNptSYFaGbW1cqeGLIhcCuwIrAqaXTFy8Bn8vKhLYnOzKyLlG0hfxO4BVgbeL2mfDLwkSbHZGbWlcqOstgO2CYi3pIW6vp4Cli36VGZmXWh/kwutHSdsvWBV5oUi5lZVyubkG8hzWXRIyStTJqc/udNj8rMrAuV7bI4FviVpGnAssB1wIbAc7wzLtnMzAah7IkhsyV9gHQW3ZaklvVE4OqIeL23dc3MrJw+E7KkpUmztO0cEZcBl7U8KjOzLtRnH3JEzAfmU3KCejMzG5iyB/W+DZwoqWyfs5mZ9VPZBLs9sCMwW9JDpOk43xYRuzU7MDOzblM2Ic8F/quVgZiZdbteE7Kk3YBfRMRBbYrHzKxr9dWH/GNgeM+CpAcljeylvpmZDVBfCbk4Z+cG1D+F2szMBqk/c1mYmVkL9ZWQg4XHHxeXzcysSfoaZSHgN5IW5OXlgV9ImldbKSI2b0VwZmbdpK+EfFph2UPfzMxapNeEHBHFhGxmZi3ig3pmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV0XCUhaT9y24kIq5sTjhmZt2rt2FvFxaWh5HmsfhbXl6KdCWRNwEnZDOzQWrYZRERK/XcgL2BB0gT1S+bb9sDfwI+3YY4zcyWeGX7kM8GPh8Rv4uIBfn2O+AY4JyWRWdm1kXKJuQNKFy2KXsNWL9p0ZiZdbGyCfku4FuS1uspyPfPA+5sRWBmZt2mbEI+BFgdmClppqSZwExgLeCw1oRmZtZdSl3kNCKekLQ58FHgfbn4EeDWiPD8yGZmTVD2qtPkxHtLvpmZWZOVPlNP0pGSpkp6TdLoXHaCpD1bF56ZWfcolZAlHQN8CZjIwhc+nQ0c1fywzMy6T9kW8hHAYRFxPrCgpvw+YNOmR2Vm1oXKJuR3Aw/VKZ8PLNe8cMzMulfZhDwD2LJO+a7Aw80Lx8yse5UdZXE2cIGk5Ul9yNtK2g84Hji4VcGZmXWTsuOQL5c0FPgasDxwFfAMaX6L61oYn5lZ1+jPOOSLgYslrQEsFRFzWheWmVn3KTvs7RRJHwGIiLk9yVjSCpJOaWWAZmbdouxBvVOBX0r6XKF8ReA/mhqRmVmX6s819Q4BTpP0HUlDWhWQmVm36k9CvgXYFtiJ1FpepTUhmZl1p7IJOQAiYhqwTS67G9i4FUGZmXWjsgn57fkrIuJlYDzw38AvWxCTmVlXKjvs7TTgrz0LEfEWcJSke4EdWhGYmVm3KXtiyGkNyi8HLm9qRGZmXaphQpb0BeA7EfFGvt9IRMR5zQ/NzKy79NZC/hzwfeCNfL+RIF3stE+SxgPnA0OASyLiGw3q/TNwPfChiJhSZttmZou7hgk5IkbVuz9QeezyhaTr8s0C7pE0OSIeLtRbCTiadKVrM7Ou0Z9xyIO1FTA9ImZExDzgWmD3OvW+ApxBapmbmXWNvvqQS4mIc0tUWw94umZ5FrB14Tm3BEZGxM8lHddLbBOACWXjMzNbHPTVh1xGAGUScq8kLZW3c2CfTxgxkXR9PyTFYJ/bzKwKSvUhN8lsYGTN8ohc1mMlYDPgdkkA7wImS9rNB/bMrBu0sw/5HmAjSaMkDQP2Bib3PBgRr0TEGhGxQURsANwJOBmbWdcoPUG9pOHAx4H1gWG1j0XE6X2tHxELJB0F3Ewa9nZZREyVdDowJSIm974FM7MlW6mELGkb4OfAm8CapK6GdfLyTKDPhAwQETcBNxXK6k5wHxHjymzTzGxJUbbL4izgatJIiTdIU3CuD0whDVEzM7NBKpuQNwcuiIgA3gKWiYjngC+SriZiZmaDVDYhz6u5/xzw7nz/r8C6TY3IzKxLlT2odx/wIeAx4HbgPyWtDfwr8EBrQjMz6y5lW8gnA8/k+18Cnge+DQzHZ8yZmTVF2fmQp9Tcf540/M3MzJqonSeGmJlZL8qOQx5OGk3xEWAtCok8ItZqemRmZl2m7EG9K4FNSRPWP0e+CrWZmTVP2YQ8DtgxIu5rYSxmZl2tbB/yE/2oa2ZmA1A2yR4NfF3S+/OlmMzMrMnKdllMB5YjnSBCnq/4bRHhJG1mNkhlE/IPgVWAz+ODemZmLVE2IY8FtoqIh1oZjJlZNyvbh/wwsHIrAzEz63ZlE/KXgHMl7SJpbUmr1d5aGaCZWbco22XRc5WPW1i4/1h52Qf1zMwGqWxC/khLozAzs74TsqSlgTOB/SNiWutDMjPrTn32IUfEfGAUHupmZtZSZQ/qfR84rJWBmJl1u7J9yCsA+0r6KHAv8GrtgxHx+WYHZmbWbcom5E3Ip00DowuPuSvDzKwJyl7CyaMszMxarGwLGQBJywIbklrFT0TEGy2JysysC5U6qCdpaUlnAS8B9wMPAi9JOjMPizMzs0Eq20I+A9gHOAL4bS7bHvg6Kan/e/NDMzPrLmUT8qeBgyPippqyJyQ9D1yCE7KZ2aCVHYe8CukyTkVPAKs2LRozsy5WNiHfT5qcvuho4E9Ni8bMrIuV7bI4HrhJ0i7AnblsG2Bd4OOtCMzMrNuUaiFHxB3AxsD1wIr59iPgvRHx297WNTOzckqPQ46IZ4CTWxiLmVlX6zUhl70aSES82JxwzMy6V18t5Ln0PVdFlNiOmZn1oa9E2tscFuNJoywWNC8cM7Pu1WtCjojfFMskbQGcRTpT7yLgK60Jzcysu5Qdh4ykUZKuAe4GXgDGRMTnI+L5lkVnZtZF+kzIklaXdD7wKPAuYLuI2Csi6p25Z2ZmA9RrQpZ0Mun06B2B3SNip4i4py2RmZl1mb4O6n0FeB2YBRwp6ch6lSJit2YHZmbWbfpKyFfiSzSZmbVFX6MsDmxTHGZmXa/0KAszM2stJ2Qzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6sIJ2Qzs4poa0KWNF7SNEnTJZ1Q5/EvSHpY0gOSbpP07nbGZ2bWSW1LyJKGABeSrlI9BthH0phCtT8CYyNic9IFVc9sV3xmZp3WzhbyVsD0iJgREfOAa4HdaytExK8j4rW8eCcwoo3xmZl1VDsT8nrA0zXLs3JZI4cAv2hpRGZmFVLJi5NK+ldgLGke5nqPTwAmtDUoM7MWa2dCng2MrFkekcsWImkX4GRgx4h4s96GImIiMDHX9/SgZrZEaGeXxT3ARvnafMOAvYHJtRXyBVQvAnaLiDltjM3MrOPalpAjYgFwFHAz8AgwKSKmSjpdUs8VR84CVgR+JOlPkiY32JyZ2RKnrX3IEXETcFOh7JSa+7u0Mx4zsyrxmXpmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGxmVhFOyGZmFeGEbGZWEU7IZmYV4YRsZlYRTshLDUXSgG7rjVy/09Gb2RJkaKcD6Li/LWCvi34/oFWvO3y7JgdjZt3MLWQzs4pwQjYzqwgn5MEYRP+z+6DNrKitfciSxgPnA0OASyLiG4XHlwGuBD4IvADsFREz2xljvwyi/xncB21mC2tbC1nSEOBC4OPAGGAfSWMK1Q4BXoqIDYHzgDPaFV9HuIVtZjXa2ULeCpgeETMAJF0L7A48XFNnd+DUfP964AJJiohoY5ztM9gW9md2QNKA1l13xEhmP/3UgJ97vZHr88yspwe8/mCef7DPPWTpZXhr/psDXr+TsQ/2fbPOKL7vjVKa2pXrJO0BjI+IQ/PyfsDWEXFUTZ2Hcp1ZefmJXGduYVsTgAl5cTPgoTa8hDLWAOb2Wav1HMeiqhJLVeKA6sRSlTigfbHMjYjxxcLFchxyREwEJgJImhIRYzscElCdWBzHoqoSS1XigOrEUpU4oPOxtHOUxWxgZM3yiFxWt46kocAqpIN7ZmZLvHYm5HuAjSSNkjQM2BuYXKgzGTgg398D+NUS239sZlbQti6LiFgg6SjgZtKwt8siYqqk04EpETEZuBS4StJ04EVS0u7LxJYF3X9VicVxLKoqsVQlDqhOLFWJAzocS9sO6pmZWe98pp6ZWUU4IZuZVcRinZAljZc0TdJ0SSe0YPsjJf1a0sOSpko6OpevJum/JT2e/w7P5ZL0rRzPA5K2rNnWAbn+45IOaPScfcQzRNIfJd2Yl0dJuis/33X5YCmSlsnL0/PjG9Rs48RcPk3SPwwwjlUlXS/pUUmPSNq2E/tE0r/l9+UhST+UtGy79omkyyTNyWPne8qatg8kfVDSg3mdb0n1zwBqEMdZ+b15QNKPJa3a12tt9F1qtD/LxlLz2LGSQtIandgnufxzeb9MlXRmO/ZJv0XEYnkjHRh8AhgNDAPuB8Y0+TnWAbbM91cCHiOd9n0mcEIuPwE4I9/fFfgFIGAb4K5cvhowI/8dnu8PH0A8XwCuAW7My5OAvfP97wGfyfePBL6X7+8NXJfvj8n7aRlgVN5/QwYQx/eBQ/P9YcCq7d4nwHrAk8ByNfviwHbtE2AHYEvgoZqypu0D4O5cV3ndj/cjjo8BQ/P9M2riqPta6eW71Gh/lo0ll48kHcz/M7BGh/bJR4BbgWXy8lrt2Cf9/m41a0PtvgHbAjfXLJ8InNji5/wp8FFgGrBOLlsHmJbvXwTsU1N/Wn58H+CimvKF6pV87hHAbcBOwI35Qzm35ov39v7IH/5t8/2huZ6K+6i2Xj/iWIWUCFUob+s+ISXkp/MXd2jeJ//Qzn0CbFD40jdlH+THHq0pX6heX3EUHvskcHW970jPa6XBd6m3z1h/YiFNg/B+YCbvJOS27hNSEt2lTr2W75P+3BbnLoueL2SPWbmsJfJP3C2Au4C1I+LZ/ND/Amv3EVMzYv0mcDzwt7y8OvByRCyos823ny8//kqu34w4RgHPA5crdZ9cImkF2rxPImI2cDbwFPAs6TXeS2f2SY9m7YP18v1mxHQwqTU5kDh6+4yVIml3YHZE3F94qN37ZGNg+9zV8BtJHxpgHIPeJ71ZnBNy20haEfgv4JiI+L/axyL9m2zp2EFJ/wjMiYh7W/k8JQ0l/Rz8bkRsAbxK+nn+tjbtk+GkyahGAesCKwCLzA3QKe3YB32RdDKwALi6Q8+/PHAScEonnr9gKOnX1DbAccCkRn3QnbQ4J+Qyp2IPmqSlScn46oi4IRc/J2md/Pg6wJw+YhpsrB8GdpM0E7iW1G1xPrCq0inmxW02OgW9GftsFjArIu7Ky9eTEnS798kuwJMR8XxEzAduIO2nTuyTHs3aB7Pz/QHHJOlA4B+BffM/h4HE8QKN92cZ7yH9w7w/f3ZHAPdJetcAYhnsPpkF3BDJ3aRfmmsMII7B7pPeNavvo9030n+8GaQ3vKfTfdMmP4dIE+Z/s1B+FgsfvDkz3/8ECx+ouDuXr0bqdx2eb08Cqw0wpnG8c1DvRyx8cOHIfP+zLHwAa1K+vykLH8CYwcAO6v0P8N58/9S8P9q6T4CtganA8nnb3wc+1859wqL9lE3bByx6AGvXfsQxnjSt7ZqFenVfK718lxrtz7KxFB6byTt9yO3eJ0cAp+f7G5O6I9SOfdKv71azNtSJG+lI7WOko6Ent2D7f0/62fkA8Kd825XUj3Qb8DjpyG3PB0akSfifAB4ExtZs62Bger4dNIiYxvFOQh6dP6TT84ek5wjysnl5en58dM36J+f4ptHgKHWJGD4ATMn75Sf5i9P2fQKcBjxKmn71qvylass+AX5I6rueT2p9HdLMfQCMza/rCeACCgdR+4hjOinh9Hxmv9fXa6XBd6nR/iwbS+HxmbyTkNu9T4YBP8jr3wfs1I590t+bT502M6uIxbkP2cxsieKEbGZWEU7IZmYV4YRsZlYRTshmZhXhhGyLJUkXSLq903EUSRou6TlJ78nL42pnOWvi83xW0s+auU3rPCdkaztJV+QkdWmdx87Ij93Ygue9PW+7eLu2pk5t+WuSZki6RtL2JZ/mJOCmiHii2fEXXAJ8sB9x2WLACdk65WlgzzwxEfD2Kc37kyYMapXLSTOH1d4OL9Q5LJdvQjqpYB7wG0nH9bbhPHfDoaRrQ7ZURLxJmor1861+LmsfJ2TrlAdIZ7TtWVP2CeAN4PbaikoT858t6aV8+ybp9NbaOrdL+o6kr0mamycoP1tS8TP+WkT8b+H2SqHOy7n8zxHx64g4EPgG8HVJG/bymnYlndn5u0YVlCbK/7Gk+yStJWmD3BrfO89C9nqeRW9zSZtJ+r2kVyX9VtKowuYmk+Y4Wb6XmGwx4oRsnXQp6TTZHgeTWrDF00ePJbVaDyfNPzsE2LfO9vYlzW62HXAUcAywV5NiPYf0ffmnXupsD9wbDU5/lbQy8EvSfA3jImJOzcOnkSaT3wJ4mXT677dJp/VuRTr1+1uFTU4hzbmwbf9eilWVE7J10jXAWEkb5RnAxgNX1Kl3DGminkkR8ShwNGm+4aKHI+KUiHgsIiYBvwZ2LtSZIOmvhduRfQUaES+QZm8b3Uu1dwPPNHhsrRzPX4B/iMI0rsC5EXFTfn3nkK5k8e3cQp9KmrvhI4WYXiPN67xBX/Hb4mFo31XMWiMiXpL0Y1LL+GXg9oh4qnaaWkmrkPpz/1Cz3t8k3cXC0yNC6gap9QwpEda6jtQarfV8yZBF73McLwc81+Cxm0mT2vxzvDO5ea3a2Hu28WChbAVJy+dE3OP1/Ly2BHBCtk67jDR15l8Z/ETm8wvLwaK/Al+JiOn93XAetrYmaUrGRuaSZr6r50ZSf/lmpBnYimpjj17Kiq9nNcr/Q7GKc5eFddptpFEMa5Cm8lxIPuD2LGkeXCBdsZjUr9pOx5ImNf9JL3X+SOpqqOfLpLlzb5P0gWYElMc6L0tqedsSwC1k66iICEmbk+a2fbNBtfOBEyU9RvoZfySpG+PZBvV7s3zur641LyJerFleNdcZRrrqxQGk4XjH9zG++GbgDEmr5z7nhUTEyfmfya2Sdo5FrzPXX9sDMyLi8UFuxyrCCdk6LiL+0keVc4B3kU6GgDQZ/dWkccL9dVC+1fod6WIEPS7Of98kJf07SaMi7uhtwxHxoKS7SVckubBBnZNyUr5N0s6kg3IDtU9NrLYE8AT1Zk0kaTypRT8mIt5q4fNsRuru2bjOOGpbTLkP2ayJIuKXpNbxiL7qDtK6wP5OxksWt5DNzCrCLWQzs4pwQjYzqwgnZDOzinBCNjOrCCdkM7OKcEI2M6uI/wcL04/7YDbZrAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_geovirus_camcoder, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('GeoVirus - CamCoder', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "2d14d026",
   "metadata": {
    "scrolled": true
   },
   "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>Unnamed: 0</th>\n",
       "      <th>name</th>\n",
       "      <th>wikipedia_name</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>page</th>\n",
       "      <th>wiktor_key</th>\n",
       "      <th>toponym</th>\n",
       "      <th>old_coord</th>\n",
       "      <th>new_coord</th>\n",
       "      <th>mlg_key</th>\n",
       "      <th>geocoded_coordinates_list</th>\n",
       "      <th>median_error_distance</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, California</td>\n",
       "      <td>34.9514</td>\n",
       "      <td>-120.4333</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Ca...</td>\n",
       "      <td>Santa Maria, California 34.9514,-120.4333</td>\n",
       "      <td>Santa Maria, California</td>\n",
       "      <td>34.9514 , -120.4333</td>\n",
       "      <td>34.951388888889 , -120.43333333333</td>\n",
       "      <td>Santa Maria, California 34.9514,-120.4333</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Bulacan</td>\n",
       "      <td>14.8208</td>\n",
       "      <td>120.9636</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Bu...</td>\n",
       "      <td>Santa Maria, Bulacan 14.8208,120.9636</td>\n",
       "      <td>Santa Maria, Bulacan</td>\n",
       "      <td>14.8208 , 120.9636</td>\n",
       "      <td>14.818333333333 , 120.95633333333</td>\n",
       "      <td>Santa Maria, Bulacan 14.8208,120.9636</td>\n",
       "      <td>[(16.15687, 119.92664)]</td>\n",
       "      <td>185.275483</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Ilocos Sur</td>\n",
       "      <td>17.3585</td>\n",
       "      <td>120.4981</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Il...</td>\n",
       "      <td>Santa Maria, Ilocos Sur 17.3585,120.4981</td>\n",
       "      <td>Santa Maria, Ilocos Sur</td>\n",
       "      <td>17.3585 , 120.4981</td>\n",
       "      <td>17.366666666667 , 120.48333333333</td>\n",
       "      <td>Santa Maria, Ilocos Sur 17.3585,120.4981</td>\n",
       "      <td>[(17.3906, 120.4587), (17.3906, 120.4587), (17...</td>\n",
       "      <td>3.730377</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>3</td>\n",
       "      <td>Santa Maria</td>\n",
       "      <td>Santa Maria, Romblon</td>\n",
       "      <td>12.4417</td>\n",
       "      <td>122.1042</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Maria%2C_Ro...</td>\n",
       "      <td>Santa Maria, Romblon 12.4417,122.1042</td>\n",
       "      <td>Santa Maria, Romblon</td>\n",
       "      <td>12.4417 , 122.1042</td>\n",
       "      <td>12.441667 , 122.104167</td>\n",
       "      <td>Santa Maria, Romblon 12.4417,122.1042</td>\n",
       "      <td>[]</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>4</td>\n",
       "      <td>Santa Cruz</td>\n",
       "      <td>Santa Cruz, Chile</td>\n",
       "      <td>-34.6372</td>\n",
       "      <td>-71.3631</td>\n",
       "      <td>http://en.wikipedia.org/wiki/Santa_Cruz%2C_Chile</td>\n",
       "      <td>Santa Cruz, Chile -34.6372,-71.3631</td>\n",
       "      <td>Santa Cruz, Chile</td>\n",
       "      <td>-34.6372 , -71.3631</td>\n",
       "      <td>-34.633333333333 , -71.366666666667</td>\n",
       "      <td>Santa Cruz, Chile -34.6372,-71.3631</td>\n",
       "      <td>[(-34.63881, -71.36576), (-34.63881, -71.36576...</td>\n",
       "      <td>0.614602</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Unnamed: 0         name           wikipedia_name      lat       lon  \\\n",
       "0           0  Santa Maria  Santa Maria, California  34.9514 -120.4333   \n",
       "1           1  Santa Maria     Santa Maria, Bulacan  14.8208  120.9636   \n",
       "2           2  Santa Maria  Santa Maria, Ilocos Sur  17.3585  120.4981   \n",
       "3           3  Santa Maria     Santa Maria, Romblon  12.4417  122.1042   \n",
       "4           4   Santa Cruz        Santa Cruz, Chile -34.6372  -71.3631   \n",
       "\n",
       "                                                page  \\\n",
       "0  http://en.wikipedia.org/wiki/Santa_Maria%2C_Ca...   \n",
       "1  http://en.wikipedia.org/wiki/Santa_Maria%2C_Bu...   \n",
       "2  http://en.wikipedia.org/wiki/Santa_Maria%2C_Il...   \n",
       "3  http://en.wikipedia.org/wiki/Santa_Maria%2C_Ro...   \n",
       "4   http://en.wikipedia.org/wiki/Santa_Cruz%2C_Chile   \n",
       "\n",
       "                                  wiktor_key                   toponym  \\\n",
       "0  Santa Maria, California 34.9514,-120.4333  Santa Maria, California    \n",
       "1      Santa Maria, Bulacan 14.8208,120.9636     Santa Maria, Bulacan    \n",
       "2   Santa Maria, Ilocos Sur 17.3585,120.4981  Santa Maria, Ilocos Sur    \n",
       "3      Santa Maria, Romblon 12.4417,122.1042     Santa Maria, Romblon    \n",
       "4        Santa Cruz, Chile -34.6372,-71.3631        Santa Cruz, Chile    \n",
       "\n",
       "               old_coord                             new_coord  \\\n",
       "0   34.9514 , -120.4333     34.951388888889 , -120.43333333333   \n",
       "1    14.8208 , 120.9636      14.818333333333 , 120.95633333333   \n",
       "2    17.3585 , 120.4981      17.366666666667 , 120.48333333333   \n",
       "3    12.4417 , 122.1042                 12.441667 , 122.104167   \n",
       "4   -34.6372 , -71.3631    -34.633333333333 , -71.366666666667   \n",
       "\n",
       "                                     mlg_key  \\\n",
       "0  Santa Maria, California 34.9514,-120.4333   \n",
       "1      Santa Maria, Bulacan 14.8208,120.9636   \n",
       "2   Santa Maria, Ilocos Sur 17.3585,120.4981   \n",
       "3      Santa Maria, Romblon 12.4417,122.1042   \n",
       "4        Santa Cruz, Chile -34.6372,-71.3631   \n",
       "\n",
       "                           geocoded_coordinates_list  median_error_distance  \n",
       "0                                                 []                    NaN  \n",
       "1                            [(16.15687, 119.92664)]             185.275483  \n",
       "2  [(17.3906, 120.4587), (17.3906, 120.4587), (17...               3.730377  \n",
       "3                                                 []                    NaN  \n",
       "4  [(-34.63881, -71.36576), (-34.63881, -71.36576...               0.614602  "
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_wiktor_camcoder = pd.read_csv(dir_geoparsed_results+'/geoparsed-results/wiktor_geocoded_results_CamCoder.csv')\n",
    "df_wiktor_camcoder.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "5f87da79",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0.0, 17500.0)"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWEAAAF3CAYAAAB0cPpRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAm/UlEQVR4nO3debwcVZn/8c/XhMgOgbAnQhBQAoOCYR8kKGJAfzDqDIsoIEtUZBsdkEUR0NEBUQQFJWhAGBGQAc0oAsIQkB1EliQQCCFKArKDsifw/P4455JK0/d23U53V+7t7/v16tftOnWq+um6t597+tSpU4oIzMysGu+oOgAzs27mJGxmViEnYTOzCjkJm5lVyEnYzKxCTsJmZhVyEu5ikmZLmtKozLqTpCmSZlcdx2DnJDzASBovKSR9q866rfK61yQtXWf9lZLelDRiEV5/nfwapR4l9zmuzrYvSrpL0r9LGtpsvP14X2MlnSdplqRXJL0kaaqkH0h6b7tfvwxJG0g6S9IDOb5XJD0oaaKkzauOz5rT9j9ua7kbgfnAuDrrdsjrhgHbANf0rMiJ7J+BqRHxdC5+D9Dfq3WeAj5bU/ZJ4BPAt4H7+7m/ol8CVwACVgf2Ab4PbAhMWIT99knSN4BvAE8DFwLTSQ2UjYA9gEMkDY+If7QrhhIxHgD8GHiVdJzuJv2uNwA+BRwkaaOImF5VjNYcJ+EBJiJelHQHsLmkpSPi5cLqccAfgE3z82sK6zYHlgGmFPb1WhOv/xLw38UySeuRkvAfImJKve1Kuisi3tq3pLOAB4ADJR0XEU8twr7rkrQ/cAJwHfCJiHihZv1RpAStVr92WZJ2BCaS/jl8NCIeq1l/DHBoFbH1l6Tlqvxntjhyd8TAdB2ptbttT0Fu6W4LXJ8fO9RsM66wbc82pfp/JY2WNEPSY5LeVzZISSMknSnpUUmv559nSlq5zPY54d9KSoDvLvu6/YhvGPCfwIvAHrUJOMfwSkQcHRF/z9ssJ+lbkm6T9HTu+pkp6b9qu4AK3Sz7STo4H8NXJd0n6eO5zj/lbqK/S3pG0hmSlqgJ4+R8DPaoTcA5xvkRcVqxFSxpGUnfkfRwjvFvks6XtHad4zBc0jn5/byU+4I/0MdxGyvp8sL7nyHpuNpuo54+ZUnrSrpU0rPA33vbb7dyS3hgug44lgUtX1jQ0r2e9Id+uqRlciIj1428vjRJm5G6CJ4Dto6Iv5TcbgXgZmA9YBJwF6mF/kXgQ5K2KNki6km+z/Yn7pK2JXV7XNCPVvZawIHA/5C6LuYD2wNHkd7fR+ts8yVgOPBTUnfCYcDlkv4NOIfUvfBrYCdSi/ZJ4FuQ/gECmwF/LNvVkJP4Vfn9XQp8D1ifdOx3kjQ2IubU1N0cuID0T+/9pG9Rz9TZ98eAy4CZeb/PAlsDJ+Xt/q1mk2VJf3M3AccBq5Z5D10lIvwYYA9gKeA14OZC2THAP0j/WDckJdyd8rqhpNbe3TX7mQ1M6a0M+Eje583Ayn3Ec0J+vXGFsv/MZQfX1P1SLv9moWxcLjseGAGsAvwTcGYuv61Nx/HQvP8v92ObYcASdcq/mfe1RZ33NRdYoVC+SS5/E/hkzX7+BDxeWP5/ue4Z/YjxoLzNKTXlH8vlFxTKJuSyE2vqHpHLZxfKlgT+BtwADK2p/+91/gam5LJvVf2ZWZwf7o4YgCLiFeA2YKykZXLxOFJSnh8R95NaU+Pyup5W8nWUJOkzwO+A/wM+HBFvaxU18AnSSbyJNeVn5/JP1NnmxLzuSeBe4GBSq2u3fr52Wcvnn6W/IkfE6xExD1IXUP4qP4IF/e9b1tnsvCh0dUTEvfk1H4uIy2rq3gisLmnZZmMkHds3ge/UxP470gm93ST1fPb/BXiD1Kot+nGd1/wIsBpwLrBi7m4akd//FbnOTnXiObUfsXcdJ+GB6zpgCeCfa/qDe9zAgn7hcfnnlJL7/gBwPnAtqaX2ShPxjQZmRMT8YmFefhBYt842E0kf9F2Ar5K+6o4kfYVvSNIqklYvPFZpsElPklmuzP4Lr3OwpHtJ30aeJf3jmJJXD6+zyaw6Zc8Bj/RSDtDTb95MjKNJCf65Ouum5X31DFNcl9TyXijhRjppWxv3hvnnJNJ7Lj4eyOtWq9nmqYh4vh+xdx33CQ9c15G+vo8jfVB7+oN7XA+clltU40gtoxtK7vshYB4piY8ntYg74aGI6GlR/l7SjaSW4U+APUtsfwdQPPH0F2CdPupPzT83LRugpC+TWo1XA2cAjwGvk/qKz6N+w+aNXnbXWzksGI3R7xjbqCemI0kt6npqTxy+XLeWvcVJeOC6hdRC3IGUhF8hJaEe15N+v+NIreS7e2kZ1fN3YFfgSuAySbtHxG/6Gd8s4D2ShhZbw7nVvgH1W4cLiYibJV0A7CPpjIi4ucEme5P6y3s0asHfROrj/BdJK5fscvksqd9854h4s6dQ0vgS2/ZbRDwi6c/AtpLeGxEPNNwoHdvxklas0wodQ/r9Pl2ou5Ok5YutYUnvJLWSi38zD+WfLxX+WdoicnfEAJW/Lt5C6jr4OHBLRLxeqDKVdHb7SGrGB5fc/99J/Xu3Ab+S9Kl+hvhr0gm2A2vKD8rll5fczzdJLcaTGlWMiJsi4prC46YG9V8nnbFfDrhY0tu+8ktaUtK3JfX0zb5BOtmkQp2hwNEl308zvpp/XiRp9ToxDpF0hKQxuejXpM/20TX1dia1qCcX/oH8BhgCfKVmt19kQX90j6tI/fVHS1qpThxL1TuG1je3hAe260gt4W1IFxS8JSJC0h9JJ1566vZLpAtDdgb+l5QAPhMRF5fc/BTScKUz8zC3P5MSwAHAjLy+TAwzJV0E7C1pu4j4Y3/fR4P9T5I0inT8ZkoqXjG3YX4Pq7LgJNel+fnvJV1GSlSfJnXftEVE/EHSBNLJshmSilfMrUe6Yu7dwMZ5k/OAfYGvSlqH1A21HulE5xOk4Y09ziWNkDg+D4e7hfR7+jfgYQo5IiJekrQPKcnPkDSJNFRtReC9LLhyckrr3n0XqHp4hh/NP0iXIUd+bF9n/RF53XwKQ6QK62fTxxC1QtlSpFbQfGDvOvs5gZrhSbl8FeAsYA4pSc0hDTsbUVNvXN7+P3p5nxuSWqDXtfFYjgV+TjpZ9iqpL/M+0mXT6xfqDSENB5xJOjH3F9I/lJ5hgSfUeV/7lTn2NcdynTrr3kNKxA/m+F4l/UM7G9i0pu4ypH8Ws0h91k+SxgGvXWe/KwE/I31zeomURMfmn7Pr1N+YdNXk3LzvJ0jDGL8OrFSoV3d7PxZ+KB8sMzOrgPuEzcwq1LEkLGmSpCclTe1lvfJ18zMl3Zv7Ec3MBrVOtoTPI4057c3OpOvb1yedKPhxB2IyM6tUx5JwRNxA35Ow7AacH8mtpMsi1+hMdGZm1VichqitBTxaWJ6Tyx6vrZiH60wAGDNmzAemTZvWkQDNzJrU63zUA/LEXERMjIixETF2qaWWaryBmdlianFKwnOBUYXlkbnMzGzQWpyS8GTSHAGStBXwQkS8rSvCzGww6VifcL7UchwwQtIc0mWiSwBExE9I85HuQroS6WXgc52KzcysKh1LwhGxV4P1QbrrgplZ11icuiPMzLqOk7CZWYWchM3MKuQkbGZWISdhM7MKOQmbmVXISdjMrEJOwmZmFXISNjOrkJOwmVmFnITNzCrkJGxmViEnYTOzCjkJm5lVyEnYzKxCTsJmZhVyEjYzq5CTsJlZhZyEzcwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQk7CZmYVchI2M6uQk7CZWYWchM3MKuQkbGZWISdhM7MKOQmbmVXISdjMrEJOwmZmFXISNjOrkJOwmVmFnITNzCrkJGxmViEnYTOzCikiqo5hkQwbNizmzZvX1LZrjhzF3Ef/2uKIzMzeRr2uGOhJWFLscfbNTW178ee3YaC/fzMbEHpNwu6OMDOrkJOwmVmFnITNzCrkJGxmViEnYTOzCjkJm5lVyEnYzKxCTsJmZhVyEjYzq5CTsJlZhZyEzcwq5CRsZlYhJ2Ezswp1NAlLGi9phqSZko6us/5dkq6T9GdJ90rapZPxmZl1WseSsKQhwJnAzsAYYC9JY2qqfQ24JCI2BfYEzupUfGZmVehkS3gLYGZEzIqI14GLgN1q6gSwfH6+AvBYB+MzM+u4oR18rbWARwvLc4Ata+qcAFwt6VBgGWDHzoRmZlaNxe3E3F7AeRExEtgFuEDS22KUNEHSnZLu7HiEZmYt1MkkPBcYVVgemcuKDgAuAYiIW4AlgRG1O4qIiRExNiLGtilWM7OO6GQSvgNYX9JoScNIJ94m19T5K/BhAEkbkpLwUx2M0cysozqWhCNiPnAIcBVwP2kUxDRJJ0naNVf7CnCQpHuAXwL7he/EaWaDWCdPzBERVwBX1JQdX3g+Hdi2kzGZmVVpcTsxZ2bWVZyEzcwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQk7CZmYVchI2M6uQk7CZWYWchM3MKlQqCUvq6OXNZmbdomxL+HFJp+aZzczMrEXKJuFjgW2AqZJukXSApGXbGJeZWVcolYQj4pyI2AbYGLgR+BapdTxJkmc9MzNrUr9OzEXE/RFxJOmuGMcCnwZukPSApC/UuxWRmZn1rl8n3PIdMT4J7A98iNQq/hmwJvB1YBzpjhlmZlZCqSQsaTNS4t0LmAecDxwSEQ8W6vwW8I03zcz6oWxL+A7gamAC8Jt8q6Jas4GLWhSXmVlXKJuE142Iv/RVISJeAj636CGZmXWPsifS1pG0fW2hpO0lfbDFMZmZdY2ySfg0YHid8uXzOjMza0LZJPwe4J465VPzOjMza0LZJPwKsEad8rWA11sXjplZdymbhK8CTpb0VpeEpJWA7+R1ZmbWhLKjI/4DuAGYLeneXLYJ8CSwRzsCMzPrBqWScEQ8Lul9wN7A+3Pxz4ELI+LlNsVmZjbolb5sOSfbc9oYi5lZ1ymdhCWNBD4IrEpNX3JEfL/FcZmZdYWyc0fsDUwC5gNPAVFYHYCTsJlZE8q2hE8Cvgd8PSLeaGM8ZmZdpewQtdWAnzoBm5m1VtkkfAWwZTsDMTPrRmW7I/5AulhjI+A+0pzCb4mIy1odmJlZNyibhM/OP4+tsy6AIa0Jx8ysu5S9WMP3jjMzawMnVzOzCpVKwkoOljRN0suS1s3lR0vavb0hmpkNXmVbwocDXwMmAiqUzwUOaXVQZmbdomwS/gJwUEScTrpqrsddwEYtj8rMrEuUTcJrk+6iUWsesFTrwjEz6y5lk/AsYLM65bsA01sXjplZdyk7TvhU4EeSlib1CW8t6bPAUcD+7QrOzGywKztO+FxJQ4FvA0sDFwCPAYdFxMVtjM/MbFDrz6Tu5wDnSBoBvCMinmxfWGZm3aF0Eu4REU+3IxAzs25UdlL3+1h4IveFRMQmLYvIzKyLlG0JX1qzvATphp/bAme2MiAzs25S9sTcifXKJR1JGkNsZmZNWNQJfC4D9m5FIGZm3WhRk/AHgZdbEYiZWTcqe2Jucm0RsAawKVC3q8LMzBore2LumZrlN4FpwLERcXVrQzIz6x5lT8x9rt2BmJl1o47eWUPSeEkzJM2UdHQvdXaXND1PIH9hJ+MzM+u0sn3Cj9DHxRpFEbFuL/sYQhpT/BFgDnCHpMkRMb1QZ33gGGDbiHhO0qplXtPMbKAq2yd8HvBl4Hbglly2NbAF8H3gxRL72AKYGRGzACRdBOzGwlNhHgScGRHPAXh+CjMb7Mom4dHAyRHx7WKhpGOAjSLiMyX2sRbwaGF5DrBlTZ0N8n5vAoYAJ0TElbU7kjQBmFAydjOzxVbZJPxJ6k/q/itS90GrDAXWB8YBI4EbJP1TRDxfrBQRE0n3u0NSqW4SM7PFUdkTcy+REmOtcZS/WGMuMKqwPDKXFc0BJkfEvIh4BHiQlJTNzAalsi3h04AzJY0Fbs1lWwH7AieU3McdwPqSRpOS757Ap2vq/BrYCzg3z1u8AenWSmZmg1LZccKnSJoNHA7snovvB/aNiEtK7mO+pEOAq0j9vZMiYpqkk4A7I2JyXreTpOnAG8CREVF7oYiZ2aChiIHdpSop9jj75qa2vfjz2zDQ37+ZDQjqbUXpizUkLSnpXyV9VdKKuezdklZqQYBmZl2p7MUa6wHXAMsCK5JGRTwPfDEvH9iW6MzMBrmyLeEfAFcDqwGvFMonAzu0OCYzs65RdnTENsBWEfGGtFDXxl+BNVselZlZl+jPBD5L1Cl7F/BCi2IxM+s6ZZPw1aS5I3qEpOVJE7r/ruVRmZl1ibLdEV8B/k/SDGBJ4GJgPeAJFowbNjOzfip7scZcSe8nXc22GakFPRH4RUS80te2ZmbWu4ZJWNISpNnPPhwRk4BJbY/KzKxLNOwTjoh5wDxKTupuZmbllT0x90PgGEll+5DNzKyEskl1O2B7YK6kqaSpLd8SEbu2OjAzs25QNgk/DfxPOwMxM+tGfSZhSbsCv/ct783M2qNRn/DlwPCeBUn3SRrVR30zM+uHRkm4dg7Mdah/+bKZmTWhP3NHmJlZizVKwsHC44Nrl83MbBE0Gh0h4HpJ8/Py0sDvJb1erBQRm7QjODOzwa5REj6xZtnD1MzMWqjPJBwRtUnYzMxayCfmzMwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQr2OjpC0T9mdRMT5rQnHzKy79DVE7cya5WGkeSPezMvvIN1x4zXASdjMrAm9dkdExHI9D2BP4F7S5O5L5sd2wN3ApzsQp5nZoFS2T/hU4LCIuCki5ufHTcARwPfaFp2Z2SBXNgmvQ80tjbKXgXe1LBozsy5TNgnfBpwhaa2egvz8NODWdgRmZtYNyibhA4CVgdmSZkuaDcwGVgUOak9oZmaDX6kbfUbEw5I2AT4CvDcX3w9cExGeX9jMrEll77ZMTrZX54eZmbVA6SvmJB0saZqklyWtm8uOlrR7+8IzMxvcSiVhSUcAXwMmsvDNP+cCh7Q+LDOz7lC2JfwF4KCIOB2YXyi/C9io5VGZmXWJskl4bWBqnfJ5wFKtC8fMrLuUTcKzgM3qlO8CTG9dOGZm3aXs6IhTgR9JWprUJ7y1pM8CRwH7tys4M7PBruw44XMlDQW+Tbrt/QXAY6T5JC5uY3xmZoNaf8YJnwOcI2kE8I6IeLJ9YZmZdYeyQ9SOl7QDQEQ83ZOAJS0j6fh2BmhmNpiVPTF3AnClpENrypcFvtHSiMzMukh/7jF3AHCipLMkDWlXQGZm3aQ/SfhqYGvgQ6RW8QrtCcnMrHuUTcIBEBEzgK1y2e3ABu0IysysW5RNwm/NFxERzwPjgT8AV7YhJjOzrlF2iNqJwIs9CxHxBnCIpD8BH2xHYGZm3aDsxRon9lJ+LnBuSyMyM+sivSZhSV8GzoqIV/Pz3kREnNb60MzMBr++WsKHAj8HXs3PexOkG342JGk8cDowBPhpRPxXL/U+BVwKbB4Rd5bZt5nZQNRrEo6I0fWeNyuPLT6TdJ+6OcAdkiZHxPSaessBh5Pu8GxmNqj1Z5zwotoCmBkRsyLideAiYLc69b4JnExqgZuZDWqN+oRLiYjvl6i2FvBoYXkOsGXNa24GjIqI30k6so/YJgATysZnZra4atQnXEYAZZJwnyS9I+9nv4YvGDGRdL87JMWivraZWVVK9Qm3yFxgVGF5ZC7rsRywMTBFEsDqwGRJu/rknJkNVp3sE74DWF/SaEnDgD2ByT0rI+KFiBgREetExDrArYATsJkNaqUndZc0HNgZeBcwrLguIk5qtH1EzJd0CHAVaYjapIiYJukk4M6ImNz3HszMBp9SSVjSVsDvgNeAVUjdCGvk5dlAwyQMEBFXAFfUlNWdFD4ixpXZp5nZQFa2O+K7wC9IIxxeJU1n+S7gTtJwMjMza0LZJLwJ8KOICOAN4J0R8QTwVdJdN8zMrAllk/DrhedPAGvn5y8Ca7Y0IjOzLlL2xNxdwObAg8AU4FuSVgM+A9zbntDMzAa/si3h44DH8vOvAU8BPwSG4yvXzMyaVnY+4TsLz58iDVUzM7NF1MmLNczMrEbZccLDSaMgdgBWpSZ5R8SqLY/MzKwLlD0xdz6wEWmS9yfId182M7NFUzYJjwO2j4i72hiLmVnXKdsn/HA/6pqZWUllE+vhwHckvS/fpsjMzFqgbHfETGAp0kUb5Pl+3xIRTsxmZk0om4R/CawAHIZPzJmZtUzZJDwW2CIiprYzGDOzblO2T3g6sHw7AzEz60Zlk/DXgO9L2lHSapJWKj7aGaCZ2WBWtjui524YV7Nwf7Dysk/MmZk1oWwS3qGtUZiZdamGSVjSEsApwD4RMaP9IZmZdY+GfcIRMQ8YjYelmZm1XNkTcz8HDmpnIGZm3ahsn/AywN6SPgL8CXipuDIiDmt1YGZm3aBsEt6QfMkysG7NOndTmJk1qeztjTw6wsysDcq2hAGQtCSwHqn1+3BEvNqWqMzMukSpE3OSlpD0XeA54B7gPuA5SafkIWxmZtaEsi3hk4G9gC8AN+ay7YDvkBL5f7Q+NDOzwa9sEv40sH9EXFEoe1jSU8BPcRI2M2tK2XHCK5BucVTrYWDFlkVjZtZlyibhe0gTutc6HLi7ZdGYmXWZst0RRwFXSNoRuDWXbQWsCezcjsDMzLpBqZZwRNwAbABcCiybH78C3hMRN/a1rZmZ9a70OOGIeAw4ro2xmJl1nT6TcNm7ZkTEs60Jx8ysuzRqCT9N47khosR+zMysjkbJs685I8aTRkfMb104Zmbdpc8kHBHX15ZJ2hT4LumKubOBb7YnNDOzwa/sOGEkjZZ0IXA78AwwJiIOi4in2hadmdkg1zAJS1pZ0unAA8DqwDYRsUdE1LuCzszM+qHPJCzpONKlydsDu0XEhyLijo5EZmbWBRqdmPsm8AowBzhY0sH1KkXErq0OzMysGzRKwufj2xeZmbVNo9ER+3UoDjOzrlR6dISZmbWek7CZWYWchM3MKuQkbGZWISdhM7MKOQmbmVXISdjMrEJOwmZmFepoEpY0XtIMSTMlHV1n/ZclTZd0r6RrJa3dyfjMzDqtY0lY0hDgTNLdmccAe0kaU1Ptz8DYiNiEdFPRUzoVn5lZFTrZEt4CmBkRsyLideAiYLdihYi4LiJezou3AiM7GJ+ZWcd1MgmvBTxaWJ6Ty3pzAPD7tkZkZlaxxfIGnZI+A4wlzWNcb/0EYEJHgzIza4NOJuG5wKjC8shcthBJOwLHAdtHxGv1dhQRE4GJub6n2jSzAauT3RF3AOvne9UNA/YEJhcr5JuIng3sGhFPdjA2M7NKdCwJR8R84BDgKuB+4JKImCbpJEk9d+b4LrAs8CtJd0ua3MvuzMwGhY72CUfEFcAVNWXHF57v2Ml4zMyq5ivmzMwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQk7CZmYVchI2M6uQk7CZWYWchM3MKuQkbGZWISdhM7MKOQmbmVXISdjMrEJOwmZmFXISNjOrkJOwmVmFnITNzCrkJGxmViEnYTOzCjkJm5lVyEnYzKxCTsJmZhVyEjYzq5CTsJlZhZyEzcwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQk7CZmYVchI2M6tQdyfhdwxFUlOPtUa9q+rozWwQGFp1AJV6cz57nH1zU5te/PltWhyMmXWj7m4Jm5lVzEnYzKxCTsJmZhVyEjYzq5CTsJlZhZyEzcwq5CRsZlYhJ2Ezswo5CZuZVchJ2MysQk7CZmYVchI2M6uQk7CZWYWchM3MKuQkbGZWISdhM7MKOQmbmVXISbhZTd4aybdFMrOi7r690aJo8tZIvi2SmRW5JWxmVqGOJmFJ4yXNkDRT0tF11r9T0sV5/W2S1ulkfGZmndaxJCxpCHAmsDMwBthL0piaagcAz0XEesBpwMmdiq9jmuxLdn+y2eDUyT7hLYCZETELQNJFwG7A9EKd3YAT8vNLgR9JUkREB+Nsryb7kgEu/uIHkdTUtkOWeCdvzHttwGw70OKtats1R45i7qN/beo1F8Vao97FY3MebWrbKmJenONVp/KbpH8FxkfEgXn5s8CWEXFIoc7UXGdOXn4413m6Zl8TgAl5cWNgagfeQiMjgKcb1mo/x7Ewx7Ewx7GwTsXxdESMr7diQI6OiIiJwEQASXdGxNiKQ3IcjsNxOI6mdPLE3FxgVGF5ZC6rW0fSUGAF4JmORGdmVoFOJuE7gPUljZY0DNgTmFxTZzKwb37+r8D/Dar+YDOzGh3rjoiI+ZIOAa4ChgCTImKapJOAOyNiMvAz4AJJM4FnSYm6kYltC7p/HMfCHMfCHMfCHEfWsRNzZmb2dr5izsysQk7CZmYVGtBJuNFl0Iu471GSrpM0XdI0SYfn8pUk/UHSQ/nn8FwuSWfkWO6VtFlhX/vm+g9J2re312wQzxBJf5b027w8Ol/aPTNf6j0sl/d66bekY3L5DEkfbSKGFSVdKukBSfdL2rqK4yHp3/PvZKqkX0paslPHQ9IkSU/mMe09ZS07BpI+IOm+vM0ZUv2rc3qJ47v5d3OvpMslrdjovfb2GerteJaJo7DuK5JC0ogqjkcuPzQfk2mSTmn38WhKRAzIB+nk3sPAusAw4B5gTAv3vwawWX6+HPAg6XLrU4Cjc/nRwMn5+S7A7wEBWwG35fKVgFn55/D8fHgT8XwZuBD4bV6+BNgzP/8J8MX8/GDgJ/n5nsDF+fmYfIzeCYzOx25IP2P4OXBgfj4MWLHTxwNYC3gEWKpwHPbr1PEAPghsBkwtlLXsGAC357rK2+7cjzh2Aobm5ycX4qj7XunjM9Tb8SwTRy4fRToJ/xdgREXHYwfgGuCdeXnVdh+PpnJNq3bU6QewNXBVYfkY4Jg2vt5vgI8AM4A1ctkawIz8/Gxgr0L9GXn9XsDZhfKF6pV87ZHAtcCHgN/mP8inCx+4t45F/sPfOj8fmuup9vgU65WMYQVS8lNNeUePBykJP5o/sEPz8fhoJ48HsE7Nh70lxyCve6BQvlC9RnHUrPsE8It6n42e90ovn6G+/r7KxkGaduB9wGwWJOGOHg9S4tyxTr22Ho/+PgZyd0TPh7HHnFzWcvkr7KbAbcBqEfF4XvU3YLUG8bQizh8ARwFv5uWVgecjYn6dfb71enn9C7n+osYxGngKOFepW+Snkpahw8cjIuYCpwJ/BR4nvb8/0fnjUdSqY7BWft6KmPYntRybiaOvv6+GJO0GzI2Ie2pWdfp4bABsl7sRrpe0eZNxLNLxaGQgJ+GOkLQs8D/AERHx9+K6SP8W2zrGT9LHgScj4k/tfJ0ShpK+7v04IjYFXiJ99X5Lh47HcNJET6OBNYFlgLrX5FehE8egEUnHAfOBX1Tw2ksDxwLHd/q16xhK+sa0FXAkcElvfcpVGshJuMxl0ItE0hKkBPyLiLgsFz8haY28fg3gyQbxLGqc2wK7SpoNXETqkjgdWFHp0u7affZ26feixjEHmBMRt+XlS0lJudPHY0fgkYh4KiLmAZeRjlGnj0dRq47B3Py86Zgk7Qd8HNg7/0NoJo5n6P14NvJu0j/Ie/Lf7EjgLkmrNxHHoh6POcBlkdxO+iY5ook4FuV4NNaqfo1OP0j/5WaRfuE9negbtXD/As4HflBT/l0WPglzSn7+MRY+6XB7Ll+J1Jc6PD8eAVZqMqZxLDgx9ysWPlFwcH7+JRY+EXVJfr4RC5+MmEX/T8z9EXhPfn5CPhYdPR7AlsA0YOm8758Dh3byePD2vseWHQPefiJql37EMZ40NewqNfXqvlf6+Az1djzLxFGzbjYL+oQ7fTy+AJyUn29A6mpQu49Hvz/XrdpRFQ/S2dYHSWc0j2vxvv+Z9LXyXuDu/NiF1D90LfAQ6cxrzx+LSJPWPwzcB4wt7Gt/YGZ+fG4RYhrHgiS8bv4DnZn/QHrOAC+Zl2fm9esWtj8uxzeDXs4yN3j99wN35mPy6/yB6fjxAE4EHiBNYXpB/jB15HgAvyT1Rc8jtbQOaOUxAMbm9/Uw8CNqToQ2iGMmKdH0/L3+pNF7pZfPUG/Hs0wcNetnsyAJd/p4DAP+O29/F/Chdh+PZh6+bNnMrEIDuU/YzGzAcxI2M6uQk7CZWYWchM3MKuQkbGZWISdhG1Ak/UjSlKrjqCVpuKQnJL07L48rziDWwtf5kqT/beU+rVpOwtYxks7LielnddadnNf9tg2vOyXvu/ZxUaFOsfxlSbMkXShpu5IvcyxwRUQ83Or4a/wU+EA/4rLFnJOwddqjwO558h/grcuJ9yFNytMu55Jm5So+Pl9T56BcviFpsP/rwPWSjuxrx3m+hANJ90hsq4h4jTSl6WHtfi3rDCdh67R7SVeW7V4o+xjwKjClWFFpIvtTJT2XHz8gXV5arDNF0lmSvi3p6Tyx96mSav+2X46Iv9U8Xqip83wu/0tEXBcR+wH/BXxH0np9vKddSFdX3tRbBaXJ5S+XdJekVSWtk1vde+YZvl7Js9NtImljSTdLeknSjZJG1+xuMmk+kaX7iMkGCCdhq8LPSJep9tif1FKtvXzzK6TW6edJc7gOAfaus7+9SbOGbQMcAhwB7NGiWL9H+pz8Sx91tgP+FL1cfippeeBK0hwJ4yLiycLqE0kTsG8KPE+6/PaHpMtqtyBddn1GzS7vJM1zsHX/3ootjpyErQoXAmMlrZ9n1xoPnFen3hGkyXAuiYgHgMNJ8/XWmh4Rx0fEgxFxCXAd8OGaOhMkvVjzOLhRoBHxDGlWtHX7qLY28Fgv61bN8fwD+GjUTIcKfD8irsjv73ukuz78MLfEp5HmS9ihJqaXSfMir9Moflv8DW1cxay1IuI5SZeTWsDPA1Mi4q/FqV4lrUDqn72lsN2bkm5j4ekGIXVxFD1GSn5FF5NanUVPlQxZ9D1H8FLAE72su4o0ecynYsGk4EXF2Hv2cV9N2TKSls7Jt8cr+XVtgHMStqpMIk1D+SKLPgH4vJrl4O3f8l6IiJn93XEeYrYKaYrD3jxNmlGunt+S+r83Js1sVqsYe/RRVvt+VqL8PxFbjLk7wqpyLWn0wQjStJgLySfNHifNJQuku/WS+kk76SukycB/3UedP5O6Eer5Omn+2Wslvb8VAeWxyEuSWtg2wLklbJWIiJC0CWl+2Nd6qXY6cIykB0lf0Q8mdVE83kv9viyd+5+LXo+IZwvLK+Y6w0h3iNiXNHTuqAbjf68CTpa0cu5DXkhEHJf/gVwj6cPx9nuv9dd2wKyIeGgR92OLASdhq0xE/KNBle8Bq5MuUIA0gfsvSON4++tz+VF0E2ny/h7n5J+vkRL9raTRDDf0teOIuE/S7aQ7d5zZS51jcyK+VtKHSSfWmrVXIVYb4Dypu1kLSBpParmPiYg32vg6G5O6cjaoM87ZBiD3CZu1QERcSWoFj2xUdxGtCezjBDx4uCVsZlYht4TNzCrkJGxmViEnYTOzCjkJm5lVyEnYzKxCTsJmZhX6/0pH/K6JbStaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "sns.displot(df_wiktor_camcoder, x=\"median_error_distance\",stat=\"probability\",bins=20)\n",
    "plt.xlabel('MdnED (km)', fontsize=14)\n",
    "plt.ylabel('Normalized Frequency', fontsize=14)\n",
    "plt.title('WikToR - CamCoder', fontsize = 18)\n",
    "plt.ylim(0,1)\n",
    "plt.xlim(0,17500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5ffdc1aa",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}