{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas as gpd\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "chile = gpd.read_file('data/division_comunal_geo_ide_1/', encoding='utf-8')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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>NOM_REG</th>\n",
       "      <th>NOM_PROV</th>\n",
       "      <th>NOM_COM</th>\n",
       "      <th>SHAPE_LENG</th>\n",
       "      <th>DIS_ELEC</th>\n",
       "      <th>CIR_SENA</th>\n",
       "      <th>COD_COMUNA</th>\n",
       "      <th>SHAPE_Le_1</th>\n",
       "      <th>SHAPE_Area</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Región de La Araucanía</td>\n",
       "      <td>Malleco</td>\n",
       "      <td>Ercilla</td>\n",
       "      <td>170675.595438</td>\n",
       "      <td>48</td>\n",
       "      <td>14</td>\n",
       "      <td>9204</td>\n",
       "      <td>170675.595706</td>\n",
       "      <td>4.979025e+08</td>\n",
       "      <td>POLYGON ((-72.5231622316139 -37.94412233799464...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Región de La Araucanía</td>\n",
       "      <td>Malleco</td>\n",
       "      <td>Lonquimay</td>\n",
       "      <td>346907.195437</td>\n",
       "      <td>49</td>\n",
       "      <td>14</td>\n",
       "      <td>9205</td>\n",
       "      <td>346907.195899</td>\n",
       "      <td>3.932854e+09</td>\n",
       "      <td>POLYGON ((-71.34595698743479 -38.0886599063557...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Región de La Araucanía</td>\n",
       "      <td>Malleco</td>\n",
       "      <td>Los Sauces</td>\n",
       "      <td>169844.122286</td>\n",
       "      <td>48</td>\n",
       "      <td>14</td>\n",
       "      <td>9206</td>\n",
       "      <td>169844.122239</td>\n",
       "      <td>8.519730e+08</td>\n",
       "      <td>POLYGON ((-72.91458811671856 -37.8395074881169...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Región de La Araucanía</td>\n",
       "      <td>Malleco</td>\n",
       "      <td>Lumaco</td>\n",
       "      <td>227705.219204</td>\n",
       "      <td>48</td>\n",
       "      <td>14</td>\n",
       "      <td>9207</td>\n",
       "      <td>227705.218806</td>\n",
       "      <td>1.112643e+09</td>\n",
       "      <td>POLYGON ((-72.88079192476494 -38.0750506676126...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Región de La Araucanía</td>\n",
       "      <td>Malleco</td>\n",
       "      <td>Purén</td>\n",
       "      <td>137030.480474</td>\n",
       "      <td>48</td>\n",
       "      <td>14</td>\n",
       "      <td>9208</td>\n",
       "      <td>137030.480353</td>\n",
       "      <td>4.654760e+08</td>\n",
       "      <td>POLYGON ((-73.03073936806268 -37.8402500719441...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                  NOM_REG NOM_PROV     NOM_COM     SHAPE_LENG  DIS_ELEC  \\\n",
       "0  Región de La Araucanía  Malleco     Ercilla  170675.595438        48   \n",
       "1  Región de La Araucanía  Malleco   Lonquimay  346907.195437        49   \n",
       "2  Región de La Araucanía  Malleco  Los Sauces  169844.122286        48   \n",
       "3  Región de La Araucanía  Malleco      Lumaco  227705.219204        48   \n",
       "4  Región de La Araucanía  Malleco       Purén  137030.480474        48   \n",
       "\n",
       "   CIR_SENA  COD_COMUNA     SHAPE_Le_1    SHAPE_Area  \\\n",
       "0        14        9204  170675.595706  4.979025e+08   \n",
       "1        14        9205  346907.195899  3.932854e+09   \n",
       "2        14        9206  169844.122239  8.519730e+08   \n",
       "3        14        9207  227705.218806  1.112643e+09   \n",
       "4        14        9208  137030.480353  4.654760e+08   \n",
       "\n",
       "                                            geometry  \n",
       "0  POLYGON ((-72.5231622316139 -37.94412233799464...  \n",
       "1  POLYGON ((-71.34595698743479 -38.0886599063557...  \n",
       "2  POLYGON ((-72.91458811671856 -37.8395074881169...  \n",
       "3  POLYGON ((-72.88079192476494 -38.0750506676126...  \n",
       "4  POLYGON ((-73.03073936806268 -37.8402500719441...  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "chile.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio = chile[chile['NOM_REG'] == 'Región del Bío-Bío']\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f8b18c75390>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAToAAAD4CAYAAAB8FSpXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9d5ylaVnn/b2fk3OonKs6T093T/dMdw0MwgAiAhJU2FUQAZUdWJdXF8O6LKsIiqK8LrsvuEoQUZYgIDAvOI4MaWZgQufp6Zwq56qTc3ju/eOpOl3hnKoTq6qb5/v59OdTdc4T7qqu5zrXfYXfJaSU6Ojo6NzJKFu9AB0dHZ1Goxs6HR2dOx7d0Ono6Nzx6IZOR0fnjkc3dDo6Onc8xq1eQDU0NzfL/v7+rV6Gjo7ONuDUqVPzUsqW9Y65LQ1df38/J0+e3Opl6OjobAOEECMbHaNvXXV0dO54dEOno6Nzx6MbOh0dnTse3dDp6Ojc8eiGTkdH546nJkMnhPgTIcQ5IcRZIcR3hBCdy9576eLrF4QQj5c4f0AI8awQ4roQ4p+EEOZa1qOjo6NTjFo9uo9KKQ9JKQ8D3wb+CEAI4QX+N/B6KeXdwL8rcf5fAB+TUu4CgsBv1LgeHR0dnTXUZOiklJFl3zqAJc2ntwBfl1KOLh43u/pcIYQAXg58bfGlfwB+vpb16Ojo6BSj5oJhIcSHgbcBYeBliy/vAUxCiB8CLuB/SSn/cdWpTUBISplb/H4c6FrnPg8BDwH09vbWumydBhOMZ/ji8VGeGwuxr8PNO188gNtqKnqslBLtc09HpzGIjYQ3hRDfBdqLvPV+KeXDy457H2CVUn5ACPEJ4Cjw04ANeBr4OSnl1WXHNwPPLG5bEUL0AP8qpTyw0aKPHj0q9c6I7c3nnx7mDx++UPjeazfxGy8aYCaaQpVwX6+PN97XTTav8q7Pn+IPXrWPve2urVuwzm2LEOKUlPLoesds6NFJKV9R5v2+ADwCfADNO1uQUsaBuBDiCeAe4Oqy4xcArxDCuOjVdQMTZd5LZ5vT1+RY8X0okeWvHrv13//l46N88+wEb39hP9+/PEs8nePjbz5Cq9u62UvV+Qmg1qzr7mXfvgG4vPj1w8BPCSGMQgg7cD9wafm5UnMlfwC8afGlty+ep3MH8JWTY+u+r0p48to8v/+15wB4diiALuqv0yhqzbp+RAhxXghxDngl8NsAUspLwKPAOeA48Bkp5XkAIcQjy8pQ/gD4HSHEdbSY3d/VuB6dbUAsnePb56bKOjaYyBa+/tn/+QTv+vxJrkxH1jlDR6dyakpGSCnfuM57HwU+WuT11yz7+iYwWMsadLYXubzKv52frurcUCLLv12YwWYy8PYH+snkVHa1OjkxHOTZoQXe9+q7OD8ZJp7O8eLd66ry6Ois4LaUadLZnqiq5LUf/xGXp6M1XWcqnOLtnz1OJJVb8fpjF2dIZfP851fs4cW7S5yso1MEvQVMp27Mx9I1GznQ4nXJbJ5drc4Vr48Hk5gMCm8Z1MuLdCpDN3Q6daPVbWX3KuNUKQYBx/p9HOj0EIhlONrnW/H+VDjFt58vL/6no7OEbuh06sqLdjVXfe5gvx+fw8yJ4SBnxkLkVJWTI8E1xz19Y6GWJdbEbDS1ZffWqR49RqdTV7p9tg2P8TvM5FSVdDZPs9NKh9fKeCDJ8eHAiuNWx+gArCaFtz/QV7f1VsrpkRDfPDPBXR1ufvsVeqDwdkE3dDqbyn19Pk6NBDEoAgFMhJJMhJJln28zGfjWc5MYhOD8ZJgur51j/b5NayFLZnM8ekHLKuuta7cPuqHTqSsL8cy671+b1ZIVebW68uBgIstf/+AGf/2DG4XXDnS5OdDp4c9/8eAKw/NPJ0Y52OVlf6e7qnsVYyqsbV1/fH2ed/7DSR56yQ7u39FUt+vrNAbd0OnUDVWVfOfC+jV0kWSOPr8Nr8MMEp4bD9d83/MTEc5PRLCaDDy4t4VcXvLVk2N85+IMPX4bf/rzB3lwT+11dyeGA/zNDzUDG03n+N7lWZ68Ns8H33A3u1qdHOv313wPncawYVP/dkRv6t+e/Pkjl/jkEzcrOudYv4/nxsNkcmqDVgV3d7r5l996cU3XuDQV4d9/8mmiReKGAA/uaeHv3n4Uo0HP72025TT16/8rOnXje5fXyA5uyMnhIAe7PA1YzS2clto3Lh/81oWSRg7g8atz/LdvPF/zfXQag7511akb93R7uT4bq+gcibblrQd9TXbe8UA/N+fi3JyPcXkqyot3N/NrLxpY97xPfP8a12dj/OoL+7i3t3hiYzaS3vD+jzw/zdB8nL944yF2tNRWT6hTX3RDp1M3XrSriX8+PV7xeQal9syl22rEZzfz6gMdtHs0qSdVlSgbXPup6/P8v9/R5KO+eXYSp8XILxzp4o9ff3dhXcPzccaCiQ3XEEvn6Gty6EZuG6JvXXXqxpFe38YHFeHsWBBjlcZuoNnBkR4vkVSOs2MhXvmxx3n4rCZrqCiCbF7l00/c5D/+n1P8wdfOEU/f2n7mVcnHvnt1xfVi6Ryff2aE+/70Md77T2cZDybI5lWy+fK8zq+dGufCZGNjjjqVo3t0OnWjz2/HbjaQyOQrOm9ni4srM9X1yLa4LIQTt0paIqkcH//+dV64owmDIvjQty/y8NnJwvuXZ6J87d0vxGRQ+MT3r3NieG3nBWhKKt84M8F8LE2z01LRmt7wiR/z9Pt+mkxexWsz4SgRI/zgty6wq9XJK+5qw2kxYjcb9Lq8BqFnXXXqyis/9jhXZyqL0/X5bbR7bJwZC1XsCR3r9xU1VjaTAYkklV17PZfVyL8/2s3f/Wi4vPU12RlZ2HjrupyjfT6mIynMRoUdzQ7MRoX3/9x+sjmVJ6/N8YMrc3z/8ixWk4KqgsWkcLDLw+d+bRCzUd9oVUJdpNR1dCrBajJUfM5IIMlIIMmeNmdFRtJnNzEdLt57msyW9iq7fTY+++Phsu8zHkxytN/HRDDBVHjjpASwokf35lwcgB9emVvj7S4Z4kxe5akbC3zp+Chvf6C/7LXplIdu6HTqRiCeqVqmqctrqyhje7jHw2w0zViw/PaxY/0+cnnJWDBBJRuZvCo5ORykw2MttJvF0zlCiSyBeBqv3UQsnSe1aFxLxfPK2dJ/6NsXeWBnE7vb9EFB9UQ3dDp14wvPjFQdhG9ymivqeRUIJkOVKYmksyrnJqrvxJgKpwotYMtJLvPy7GYD9/Z6yKsSISCnSs6Mhsq+R16VvPFvnuK9P7OHdzzQr8fs6oRu6HTqQjiZ5bM/Hqr4vD6/DZvZyLkKW8EW4hkElD1Qx2U1cmOusthhNSQyeZ4dWqnCcv+Af81r6xFJ5fjgty5yd6eHwQG9rawe6FFPnbrwNz+8sWLQzUY4zAaMCliMSlXb3dFAggNd5Tfr7+9wE68wG1wvnh0KVGWwfutLZ/jgty7w8NmJuhVV/6SiGzqdmnni6hx/+7jW7G4zKdzb66XXb2dwwE+xjdc93R4S2TwtLiuZMuvTihFO5nBZykt+XJzc2slix4cCDFbY9D8dSfH3Px7mt798ll/5zLMkt8hQ3wnohk6nZr52aoy9bU46PBYOdHk4PRpiNJDg+FCA+/p8OMwGdjQ7GBzw0+2zcW4ijJRazKtSuSYhtNKNwX4/drOBds/GQp/AGqmmPr+NIz0eDnd7aXNVVidXLceHAxzpqa6v9+mbCzz0+ZOcGA5wO5aEbTV6HZ1O1Ugp+dC3L/LlE2Prehs2s6Hk+7tanYwF4qRz6/8dNjvN9Dc5CCWzFffTLnH/gJ9APEM8nWNyWVKhx2fDYzORUyXTkRSZbJ5Ekfq7euCyGDEYBKEKtvmrecEOPy/f18qr7u6gt8lex9U1jktTEYSAfe310wZcopw6Ot3Q6VTNszcX+KVPPVPzdXa2OLixWGtWjENdHoYDcSLJ0uoh5XJPt6csDbxWl4X+ZgenRwLUu5trf4eby9MRag27mQ0KD+5t4UCnmy6fndff07ntio1T2Twf++5VPv3ETVSpteyZDQqff+cgrS5rXe6hyzTpNJTHr85xoNPN4ICfQ10ejvR46PVX7mE0rWqx8tiMOM23Ym8Wk1IXIwdQNGhYhNlomuNDAe7pqa5/dz0uTkXWTDerhkxe5bGLMzx5fZ7f++pz/Nz/9yT/cm6qavXmRvCVk2N88vGbBaM+NB/nykyUn/rIDxiaL/3hVm/08hKdqnn65gLniwT5j/X7EAhSuTzPT4SxmkpvXYHCexaj4J4eH9PhFM1OM6cX68+UOtSS9fptuKymistYTo0EaXaa6fBYsZmN5FXJzblYRRnmYhwfDnKsX5ufsdwuHer2kMnluTxd/vZcWbTe12Zj/KcvnsZuNvCPvz6I126m2WnGazcX4nqbUZeXzOSZiaRIZPJ8/PvXix6TVVVaNyk2Crqh06mBiRJdCct7T1tdFgLxNEf7fEVHFwI8PxGmr8lONJXj+GK9WTCRwWMz4raauDRVm9z64W4vZ8dDQPkFycuZj2WYj90SDjAbtcxyKqtycar6bO6J4SC9fjttbgvzsQwOs4Er01EcFiNHer08NxZad3vb5bXhtBjXTE9LZPK86W+fBhblqxxm7GYj3T4b//kVu5kKpRgJJPh3R7txW01Vr78YV2eivPv/nCq0vZXCZFDq8gFWLrqh06kKVZWkywhezUa1roHRQIJ97U6cFlNRg7e6aX53qxOTQeHZoQD72l1EqmwtA23rW08yOZXToyGaHGa6fTbGK2hDW81oIMFoYOXPns5lCMQzNDvNWE0G5qJpun02zEaFS1Pa7+G+Ph/nJ8IbdpNEUrnC2MhLUxEeuzhTeO8zT97k0287yoE6Kjx/5smbGxo50H6Hz42HeMEmDRbSDZ1OVfzFo5cJJ8vfvs1G08xG01iNCp0e64qsZzFOL2ubujwd5UivF6MiMBsUoqlc2a1cDrOBC3WuoTMoggNdbm7OxVkIrj/1rBaWe5FLyZrBfh8GRXBiOEiuxljcVDjF6z/xI775n17EoW5vRefORdNI5JqEQiVF2X2bmDHWDZ1ORSQzed7ymWdQi2TrFcGGmcR0Ti06mHojlvpFj/b5sJnLKxLe1eIgp0qGK5RY2ohev53nx8M4rUbcNmP9EiVlcHw4iMUo8DvMBW+5FlQJk6FkRYbuh1dm+Z2vPEev386X/sMLsJkN/ODKLF96dpTvLPMY12Nfu6tinb9aqMmnF0L8iRDinBDirBDiO0KIzmXvvXTx9QtCiMdLnP85IcTQ4nFnhRCHa1mPTuN5+uY8Z0ZDPDcWxms3cajbw2C/j8EBP+0eK81O87rnS2BHs6Pq+4eSWcaCCfa1b6zuoSii7kYOtMyhy2oilsqxr63+dWEb4bCYNpSIr4QvPDta9rGjCwne/43zBOIZzo6FePOnn+F7l2b4wMMXyjZyoHnp/3xqfNOKn2uqoxNCuKWUkcWvfwvYL6V8txDCCzwFvEpKOSqEaJVSrhkRJYT4HPBtKeXXKrmvXke3NcTSOd7x2eMlkwoGoWX1NtpSrZeYKIXVpLCv3c2FyXBBBqmvyU48naPZacFkEBgNSsHz293mxGszlVQQrpXDPR7momkmKlRQKQeB5h1LSnvIg/0+jtfxZ/v8bwzy4t0bz7799BM3+fAjl+p2X7vZwOO//zJaasjANryObsnILeLglpjEW4CvSylHF4+rfA6ezrbDbFDW3W5IYG8ZntZ4MEGXr7zWrSUOdnk4OxZaofU2spBgPqZp4D0/EeHMaAiPzcRAs4NrMzFCFcQQK2VoPkGHt7KfoVwGB/zkJXjtZo4s9g0f6vbQ36T1D5sMoqhyci385aNXNjzm66fH62rkQMsQp9YRSa0XNaejhBAfFkKMAb8C/NHiy3sAnxDih0KIU0KIt61ziQ8vbn8/JoTYvE27TsVsVHWvSlYMnynFdERLSnhs5ZU2HO3zcbpMDzCczBYKUa/NxMoyvNUQTmZJ1bnJ/li/j33troKkUyCe4cxi3/C58TDDC1r/8JFeH4FEhnpWZzw/EeafT60/wa1aUdX1UITm1TWaDbeuQojvAu1F3nq/lPLhZce9D7BKKT8ghPgEcBT4acAGPA38nJRyxcglIUQHMA2YgU8BN6SUHyqxjoeAhwB6e3vvGxkZKe8n1Kkb7/r8Sf7tQuk4TJ/fTrvHypWZaNm9nPcP+Dk3Hi4pfb40xnAkUHmsrb/J3pAYHWjlL9eK9NwODvgIJbJ47VpxcinPy6gIXFYju9tc5PIqQ/PxsouQ7+3VEgdLmek2lwWXzcj12do7DV6ww89AsxOf3YTHZuKN93XT7LSgqpIXfuR7zJQx37YSdjQ7+NJDL6DNXX07WF1mRkgpX1Hm/b4APAJ8ABgHFqSUcSAuhHgCuAdYYeiklFOLX6aFEH8P/N466/gUmjHk6NGj26fH5Q4lncvzycdvcrDbw8v2tvLXP7i+rpFrdVmYDCcrNkjPDgXY2eKgw6NJNl2aihBdlpW1mgxVGTnQatScFiOxMrzMSlnK/B7u8WJQBIqARCbH8aFbnmev34bdbMRtM5HK5ommcoQTGRwWI2PBJG6rqVAgXQljwSR9TXaO9mtdJOPBJDPRNH6HmZ0tjpriks/cDPDMzVtrCsQzvO81d3F8OFB3Iwfwt796X01GrlxqKi8RQuyWUl5b/PYNwOXFrx8GPiGEMKJ5a/cDHytyfoeUckpofSk/D5yvZT069WN0IcH/eOwq73/NXbxsb+uGD2Qwnim7j3Q1N+bi3JiLc7TPRyaXZ6DZgddmQpWShXj1dWotTgszdSjBWI3ZqOCxGrm318uFyXBJ5ZXRQPFi3kAii0Ab1ViNEZ+Lppkr8nMF4hlaXRa8dlNN6ijLeXYoQDiR5c//9fLGB1fB0zcW2LMJ8zFqraP7iBBiL6ACI8C7AaSUl4QQjwLnFt/7jJTyPIAQ4hHgnVLKSeALQogWtEfk7NL5OluPxah5LKdHg/z6507w/AYFullVcrjby4WpcNnDnleTUyX39vp5+uZCVeevpt1rrcrQtbktOCxGZiNp9ne4Fw24BARSSuxmA09cnV8j424yiLJ/9mMDvhXeX724PB3laL+P0YVEXerszo6FuOdD36nDyopzciS4KVPPajJ0Uso3rvPeR4GPFnn9Ncu+fnkt99dpHP/yvBZVePTCdNkTs86Oh9jR7CCazhX1OErhshjY1+EubLnqFVurRn78loyTtv7VfaSlaHKYcVmNtLot5PISo6KgIgnGM/jsWm2hogjyeUkgkSFWRdF0uZwcDuIwG+jz2xgp4VVuB/wOM//lZ/duyr30zgidFUgp+eQTN/mLRy8vfl/Z+VPhZNmdC0sMtDhXxJXmomna3VZa3RYuT0Uqllu3mxR2tDjLmhLW7rYyHdGOu3/Av6bvtBxcViN9TXZOj4ZKGOjNkyMq3DGTx5bJc6jLw0wk1ZAtfK20uiyb1h2h69HprEAIwa++oI93VLmdSGVVAvHy40M++1rppHgmz3QkxbnxcFUN53s73JyfjKyJ7zU5zCsKUw91ewjE0+xrdzHQ7ODZoUDRcYYbEU3lMBq236M0H8twbiJMIpuvi/5dvbk8HeWpG/Obci/do9Mhm1cxLXtQHRYjv3Ssh889NVzxtSTaJ3VOVVGEWNGYvpp97S5sZgPBdeaenhsPs7/DjdNqZHQhznQZmT/zsp9lb5sLq0nBajIwH0szNB/n3l4v2bzKyEKcTF7WpT5sOyt1R1M5hhfi3NfrYzKcREpZ1u9xM1gtutoodEP3E8ToQoJgIsPdnW6MBgVVlfzZI5f4f16+G499pUfitZv43Z/Zw189drXE1UqzFAQfHPAVDF1/k51WtxUpJUII5qKpsgxMTpUFzTe/w0SHx8JUuPRDqghtelav347ZqHBlZu09TlcwUNpvN9HTZGdsIUFgnUzmieEgO1scNDnM3JiL15QtbgTLNfXMRoWDXW7MBoWrs7EV5Tybyd2dbna0VN/3XAm6ofsJ4tELU/zZI5f5+cOd/M9fPsKjF6b5P8+O8F9fvW/NsR0eG+95+S7imXxhlGElCLT43j3dHmajaRwWY1U1Y8sJxLNYTQpH+3zcnIsTSKw0Jl67iR3NDhZiGcaDiaJ9oi6rkS6vdYWCb6/fTovLgiIgmMjS5DATTWURQhBKZgsCBv1NduZjaWLp4sXNS2Uy+9qddPlsTIaS63q0W0Ump/L8hPbhsb/DTTKb31RZ8yVuzMX48bV5Xn2wo+H30g3dTwg/ujbP3/1oCKAQh+r123nvK/aUjC8JIRgLVpf9lKxUGq40QVGKVFbl5EgQn91Ek9OM327GoAjcVhNnRoMlvTVFwLF+rQvj8nSMuzpcOC1Gzk9G1ohfFhP/DiWyhBJZjvR4OTO2vke4ZEQVAXvanFydqW5q2WZwcSqC2ahwrN9HLq9yZqw2NedKSGVVvv38lG7odOrHXR0uFha9i9feo6lpNTnNvOvBnSXPefrGAo+en67L/dPZPP1NdqwmzeDVGhcLJrJ0ea0YDaKgulsMbTtpIZbOFXpIgXXPWY/zk2E6PNaykhaqBK9tfdmq7UAmpxY+lAYH/MxF06RzeaSUSElD43nPjYX4+ulxfvHe7obdA3RD9xPDmdFQQT6p2aE9fB3rDH9eiKV5zxdP122i1EQoxf0D/oKxOdjlZiGWQZUUyjuqueZGwexUTi27Fq4csnlJXpVldx/cmItxqNvD+YlwzeMNN4Ni4YVj/VrG9upMlHCdRUbHg0k++m9XSGbzvPZQZ9lCD5Wy/XLiOnVHSslXTo4Vvm/3bNxb+K/np+seUF/ec/r8RITJcAqDImpSr3BZ1/+s7izjZ62U2Wh63balfe0uBvt9DDQ7MCqioOBbi+DoVnJiOMiJ4SCKEDWVqbzq7nZetreFJsdKL3cqnOL93zjfkJ7kJXRD9xNAIJ4pbBXb3Va6ytBR+6VjPXU3EtdnY6wWxp0IJXGYjRzr99Gxwf0UAff1+ji2qGjsd5h59maAF+0qPWAlU+/p04uMLiQwFOntdVoMhJJZjg8HGZqPk1VV5mOaGm8okSmIANyOBBNZTo4EGRzwV3xuq8vC37z1Xv7+1wb53VcW74Z4frz8bHil6FvXO5xAPMN//qezTC5Oi3rvz+zGXcb2IJ1TMRkV9rU7AYHNpBQC1QZFVLWlzauSXv/a9q65WJq5WJr7B/yF2Nexfh/RVA6b2cCZ0RAui5H9ne4VcTbQYnCnhgMM9vuRaJPJFmIZWt0WDIoglc1zb6+XXF4STecYmo+XHWNbj5loitYiogH7OzwrtsrLi6cDiSyBRAibSeFQn48zY6GGGeJGkq1izccG/IWZsoe6bxWB72xxcE+3l5/Z38bL97XVbY2r0Q3dHc4/nxrnyWu3qs89NlMhIbAe//DU8JoRhMf6fSQzeabCKdI5taKthsdmotNrXTcJ8OxQgL3tLtLZvDZha3Hr7Lebkcg1Rg40Kaa7Oz1r4nClxgAe7HJzYTLC4IC/pnIXKcFlM60xdKdGAtw/4OfEcKBkTC6ZVQtjHJ0WIxemIusO+N5ujAUT9PntFSmvPHFljlQ2j9Vk4O5ON5/7tWOYDAoP7GzalKHauqG7w3n7A/04LEZ+eGWW71yc4b9/8zzjwSQv3NnEzhZnUaN3YTLMx79/bc3ry8tF+vy2sgxdu9uK2agVJ5eT6Uxn8xgVsSI+uLpebjnZvCx4q+WwVD92YijA4R4v8XSO6UiKeDq3brKg22sjp6orMpB++9qMal7CQjyNyaBsOPd2KZxwqNuzpg1uOzMfy6C4BGaDKLsPOZrOEU5msZoMCCF46d7WBq9yJTUNx9kq9OE4laGqmq7bsQ9/d8XrioC+Jge7Wp10eKy0uqzMxVI8fWNhw9ovp8WAxWgoGCSPzcTedpemGCIoTA85NRJcI2e0Ecuzs+VgNxvo8tm4VkO9Wn+TXZtRKiCRzmG3GEAK0vk8ioRrc3HyqqTbZ8NrNxFN5QjE08xGtZ9/6UducVpIZHIVzTcFzVuejabXeNHbFbNRIZdXy84kv3mwlz//xYMNWUtdFIZ1bj/OjAZ55maAd71kB4qiTeX6wP+/VtNUldrovqH5OPf1evnH0fLl6WPpPB6bmRfvbiadyzMbSdfc+bBEpYYxkclzbSbGPT0enquy4HV4IbGuNNS9vV5Oj4ZWyKe3uizsa3cxFU4RTWXp8tpocpiZj2Xo8Bq4XkRqvRRL3nKlRn6rkFJWWC6ztQ6VbujuMM6OhfjdrzzHzfk4Tou2LX3qxgI/ulZaJaLU+EGzUcFiUNjZ6lwxZnAJi0nhR9fWClDWQn+TvSqD2emxkq7zZKzlFEu+zEbTK8Qtx4NJxoPaNrrLZ6Pdbam42DaUyFQk4LlVKEJgMVJSXXk5foeZ//aauzZhVaXRy0vuIGajKd79+VPMxbSH6w8fvsAfPnyBfz0/TSqX554eD8f6fexqdRbKPPa0ObEYlUImrNVl4d5eLwc63TQ5tFkHZ8dC7G1zYVsVz2tymOv+OR1KZqsqGu322xsypWqJ58bD7Gp1ln38RDCJyaBwX68Pr738n+fKTAyLQeGensrlqTaTdE7lULe3rGP7muy4rI0pBC4X3aO7g7CbjXzm7Ufp9Nr4rS+d4UfXb3lx2bxcsa1zmDVV30A8w49vLOCzm9jX7uTydKyoBPf5yQg7WxwsxDM0Oy3k8qo2ycpiJJXL180DCSWyhW1iJVSjJlwO3V4bnV4b05FkxQZ4LJhkLJjEZBDsaHbQ5DQzUobEeSyT58JEmMEBP5OhpDZUp4EzaqtlZCGB325aV9UFyhuB2Wh0Q3cH4bQYC0KVf/S6/bzyY0+UPDYvJXlVLahWBBPZDcft3ZiL0+qyFGJPwwsJnBYDdrOxrg9iNSKWeSlxWY11lRwyKYJ2j7VQurK0La2UbF5ycz7Ozfk4itAUQ1S5vg5eTl3ZjrW7zUkqk8ds1LK52ZxKt9+OIrSZrPUeaF0Os9E0HZ6N9eTu7d160U/d0N2h7G518me/cICvnBxDEaKgPzcXTTMeTHCoy1tVD+hqb8RlNaEooq6GrpqqqjOjIQb7/Tw/ESJZp4e+y2fj/J1wv90AACAASURBVOQtL7geTqMqKejr9TXZyebVguR7f5Mdk1EpZI+Xl28Uyygv1fC5bUZ2tzoLpTOVsLNFm7g2GkxWNOcDtDDHetqAoO0cXrSrueJ11Rvd0N2hCCF4y/19xNI5/uyRlaPqzEaFemXBpsIp7unxMFGlt1OMZ4cCHOv3VTyf9PiwVht3dgMZpXLYjNq2kYUENpOB+/p8BOIZcqrK8EyMTq+VFqcFn8PMTDiJRKzr/UWSOW7MxtjX7lr3uE6PlR6/HVVKtBpdweXpCDfm4rS7LRzocnO+AmM5G03T6bUWnc3x6gPt/PJgLxajwgt2lG7R2yx0Q3eHc6x/bV9iJqdyfDhYc3cAaLV4jch2Lin2+uxmnhsPlR0DHJqP0+W1leyMKBersT76eRuRzOY5NRLEZBC0uq0c7fNhUAQXJkLEMrd+r0d6vGTyKhcmixuiRFbl8nQUn91Er99OXkquTEfp8du1ATRScnEyUrJ0ZTqSJpLMoojKPNd2t5VEJr9CyaXZaeaDr7+b1k0YTF0uuqG7w1k+DGY1FybCDPb7EUIzfpenIxVt+4wK3N1Vfe3aRkyHU7isRnr8dm7OlaeAG05m2d3qZDaaqjpBsrvNWTdpp2anmTa3FYMiiCSzRFKaQejy2rCZjahSks2rXJyMMBFMlvSMz4yF6PRaubvTXdLYwVKs9db/x825eNm/u0RWZbBf68Et93d3ejTEQLODVCZParETZD6W4cRwkJ871MHJ4QBeu4ldrY0fUr0euqG7w3FaSv8XxzN5JLIwSPlwjxezQSGSynBtJobJoHCgy8O12eI6ZEd6K99eVkKXz8bZKozoyZEgL9rZRCKTJ6eqWIyaosjQXIxynl/XOr+zUnhsJna2OEjnVDI5FZ/DRCKdZzSYKGqYKpmUtsRkKMVkKMX+Djd2i4HnKjBI5XJ8OMiRXi/nxsNlCzcMzcc52ufj+lys4Nm99ytn+avHrhSM7Jf+wwt44c6t28Lqhu4OR9lAEmh4IVFoX1qKbfU32enw2khl85weDbKnzUU4uTb2c2Y0iN9hquqhXbHGEtul1XV7lRBL5xYHUd/CYTaws8XJ+cnSIpgdHsuGUumrMRkEvX5bxSUx1bKUzPDZTfT47QzNxYnWsYTjzGiIA53uNQXR63FyJMjdne6Cocvk1IKRe83Bdl6wo3Jpp3qiG7o7nI0kdeaiaQaaHSuGo6xuhTIvyjWNLCRWbG1zKuxqddUU5+vz21iIZ2l2mvHazQgBubzKbCRdXfp1EUsRIxnP5Dk3EabVZWGg2UEmlyeayjETTRNN5fDaTQw0OzAZFDo8NkYCCaaLyDnZTAb8DjNum5FIMofLaqwq41krS9tUk0FsuKWtlPOTEXa3Oklm82WX7FyairCr1bmi9c2oCH7zpbs2RaFkPXRDd4cTSGTY0eJgKpQimS3eaO53mNedAuW2mvjR9TBNDjOdXhM3lsV8Lk1FsBg3VuoohcloIJZOakooqwxssenyu1qdTISSHOrycHwoQLPTUugEWc56YhWrPRWX1ciuViddHmvBuI0GtFjZkpSSEFo9nNmgcGMuxkQoycTmOHAbks1LLkxGONTl4dJ0pG7b2WuzMfwOM8bFfumNUKU2U3d5EfFfvulQVUPI641u6O5whucT3JyLYxDa6MFoKsdoILHiD9e4anu7u9VJNq8yvJBgT5uzMNNhIZ4hksrS47NpFf+KtmXz2y08eb26ievNDjO5JjvpnEo0lVtX+sluUgrewrNDAQaaHVhNygpD53eY2d3qrEjQMprKEU3FuD4bw2QQ2sPqMNHhtXJmk7aj9eDcRLjuogBGRXC418uFiQjJbJ5mp5mBZgeqKgmncmuECy5ORejx2UAIAvEMlk3KXm+EbujucJbkqfOSQszKaTbQ6rEWYiiZ/C2joAhAgMmo4LGZCCeyKzyrbF4Sz+TZ2+bi5nyMC5NRIEqba63abjmcHQvR5bMVFH/b3RZ6/Q4kEoGmZKwomtDlQjzN0Pwtr28umubuTjf3D/iREhKZXGHalyKoqswkm5dk83nimTzTkRS7W53YzQbSObWhvbTbleXer8dmJJzIFhJQR3qK97qOBZOYjQr9TXa+eXYCs1HBYTFwfCjAfCzNf/+5/WWJv9YT3dDd4fzw6tya12KZPOZlYpY3ZmPYTArJrIrRoBCKZwklMxhLdDwE4hkCqwbneOwm5mJpBLC7zcV0JFV0Sla3z0b7YrnFs0MBLCYD48tmx05H0iUVP35qV/MKQ7d6hOFyVKltSUsps5TLkiyTySDY0+bEbFA4X8dYWL1RG6gvuTrzbjKWbtXL5FRMBoVALM1fPHoJv8NCNq/isBgbusZS6IbuDmZ0IVGyur/FaSlkS50WIztbHTx5bYG9bS4uL8Z5Kon1RFM5FAEHu72cGQ1xqNuD357j5mLsz6AIWl0WWpxmTo4EEWjxL4/NxMkya9aSmcoyi5eno+vWEVZCNi+5OhPDZTGyp82Jx2ZiNpJiJFC/jpB6cGI4yI4WR9m1czXda7GD5WyJMpdrK7a12nr+8dcHsZs33+zUdEchxJ8AbwBUYBZ4h5RyUgjx+8CvLLvHXUCLlDKw6vwB4MtAE3AK+FUpZX1n7P0E84XjpYU0r8zEONLrJZdXeX4iwnQkxYt3NRFJ5Qp/tC6LkUxeLSvRMBVO4bGZCrLm58bDCKFdY2+7i+cnQkyFU4UtqqTyIdanR0O0uS3MVKDxNhdNs7/DXSjJqJVoOkd05paXt7/Dhd1sXFFDttnsanUUtpeRZI7QOtLz9USiGdYOjxWb2aDFghUtA3xpamVSxGUx4nWYePHurel7rUlKXQjhllJGFr/+LWC/lPLdq455HfBeKeXLi5z/FeDrUsovCyH+FnhOSvk3G933dpRS//sfD7G/w839m9T3F4hneMlf/mDd4L7TYiCWXpmJNRkER3p9ZHMqZ8ZCHOn11hSQb3VZMBmUmluyQPMAr8xEqeZPdnDAz5mRINkGyTnZzQZanJaKBsbUitWosKvVyfnJCDaTgWQ2z55WB5F0vmhZTCN510t28LYH+vHZTdjNRv71+Sn+6eQYqoSX7W3hHQ/0N6zEpOFS6ktGbhEHxTvF3wx8qcjiBPBy4C2LL/0D8MfAhobuduTxq3NcnYlumqH7q+9c2XB4zWojB9oWbakuzmUx1FLKxr52F8F4pi5GDrSYXLWfy8eHAliMCkd7vMxEUozVUYQANDn3yXCSg10eLEYFVUpyecl4KLkinrlsnEZNDA74OTMaLMQLl0qHrGYjs7HN3xR96smbPHZxhlcfbOc9L9vNqw608+qDHZu+jlLUvFkWQnwYeBsQBl626j078CrgPUVObQJCUsqlp3Ec6FrnPg8BDwH09vbWuuxN59pMjB6ffdPu99SNhZqv0eKyVl3tLwRYjAr1/BCvtlZv+fknR4LYTAY6PFZcViNem5mRhXhVGePVZPOS5yfWxkSdFiNWk0IsncOkKPQ22XFajGRyKrFMDgFMR1L47VpfbF6VKEKAkCC1zPPZsWChfe1on29NkbZBaEkgm0nzLG0mw4rZtX1+Gw6LkaszURoxSlZKuDkf569/cIO//sENDnV7+MI7799yZeElNty6CiG+C7QXeev9UsqHlx33PsAqpfzAstd+CXirlPJ1Ra7bDDwjpdy1+H0P8K9SygMbLfp227oG4hnu/ZPHEAJ++VgPv/nSXfT4G2v0vn56nN/5ynNVn18PBZAljvX7CCezG04WWw+H2UAmrzZslkKX14bdbFgVQN8+eGwmdrQ4iCSzDC8kVvShtrostLotKySWenw2zWhKyfWZWKFFrNlppstrW9Me1wjedF83H33ToYZ3RZSzda3buEMhRC/wyHJDJYT4BvBVKeUXixwvgDmgXUqZE0K8EPhjKeXPbnSv28nQ5fIqH/r2Rf7x6VuJAUXAK+5q409+/gBtDZKyuT4b5dRIEJvZyF8+erkqddy+Jnvdxu+ZDAIpZdXeRD0TCuvR32Sny2fl/ESkqJDBdsBlNdLfZCeUzJLM5JlfZ6taTIpLEbCzxbkpRv1//fJh3nC45EatLjQ8RieE2C2lXJp0/Abg8rL3PMCDwFuLnSullEKIHwBvQsu8vh14uNixtzNGg8K7H9zJP50YI51T2dni4GCXh7FgkuuzsYYZul2troI0zoXJMJ98/GZF5wu0uFO92NHs5MpMLQW3m1N7NbyQoNlpIZtTcVmMtHusSNAG3EgtTmgxKoQSWYwGwXwszc4WJyaDFpdL51QMiuD8eBgW28aqRRHQ7LRgNip0eKxk8iqhRJb5aLrs3trheS0TutwDXPqypUT7XD351BM3G27oyqHWGN1HhBB70cpLRoDlGddfAL4jpVxR0COEeAR4p5RyEvgD4MtCiD8FzgB/V+N6tiWdXht/+Nr9/Pdvnue9P7OH1x7qBLSpXQuxNE3O+tR6leJon59PcsvQ+ewmdre5uD4bW1P4e+sYc8XS2qWwmw1MhmvbBicy+YpFIatGaNpsoBLdwOsxKqJoHPNgl5uZSJpun41cXi7OiIis8WjNBkGXz47XbiKfl5iNCoqATF5yc+7WoKJq51WYDUpRuaVrszGO9HgabuiuzcbI5rXi4a2k1qzrG9d573PA54q8/pplX98EBmtZw+3C6+7p5OPfv8bP7G8rvNbq2hwF1muz0cX7Wejx25gMpTg+FFjTF2kzaX+MPX47ZqOBQB3qsdxWIx0eW43enOZp1VrqshF2s4G9bU5OVqCxV6rZfcnjWi4eYDMZ6G+2YTcZsJoNpDJ5RgKJdQUVaiWcyvLAziYyeZXh+fiKbW4mL2l2mle85rYZGWh2FIaC14rPbuLJa3O8bG/rliqY1C1Gt5ncTjG65YwFEg1PQiznwmSYh/7xFMFEBqfFyFwsvaI841CXh0Q2x/VZ7UHrb7IXJJr2d7i4OFW9cbKZDOzv1MYp1utB7vJamSgyn6BWXFYju1qczMfSdS872WoUAQYhyKqy6EwJoyJwWrXhOpeno+zvcBc+/Do9Vjq9NoSAyVCypt/9Z99xlJfva9v4wCooJ0anD7DeRDbTyAF89eQ405EUiUye2Wh6TQ3auYkwDosRu0lBADbzrUZrRxUqu8vx2E2cGgnWzciZDYIOj60u11rCaTEy2O/HqAjOjIXuOCMHcLTfVyiSno2k1wzhzqmS0GKjvstqXOHhT4ZTnBwJcmI4iJTQ4al+B/LE1Xmy+c0fybiE3ut6B/PHr7+beDrHV0+NlzzmubEwgwM+5qJpLi3z4Got4zDUeZvitZuJpXNVTQcrRSqb5+xYsDBScDtgMoi6ldB0eq0FmXzQtAmNBq3nuJhycLFpXoX3wilaXRaO9vlQFMG1mWhhDvDuVic+h5lMTmUmksJhMZLK5lfEFT/31DCZvMof/Ow+PPbNr63TDd0djJSSVx1oX9fQAagqK1RBQBPULPVAlEOLy1y3OjxYKRd0b68XgeD0WLDqTgnQ9PlObRO9uRaXhYEmB1dno5gNChajgt9hLgy6rsb4+e3mNcZrNprmaJ+vqv/XlZJNJk2QlNXN+xqKgPv6fBgVgUTL4p+fCPPhRy7xl286VPG9a0U3dHcwH/zWRT731PCGx8WLqIKkcyq9fnvVhu7GbJxun63qbGEp+vw2jIrg4lR1Pa/L2ZQM7gYYFMHhHg835uJrJo8tbaU9NiM9PjsWk4Gr09FC8a/LYqTFZSGrqkipJVN8djPhZBan1UiuyFaxXsKcGw0sVyWcWiaPZTEqHO7x8us/1V/zvatBN3QNQlXlhoNpGk05HlW311YyjnZhMqKJLVZROJvI5tnncdXd0OVUSOVU3Dbjhr28G7GentpmIYBLk5HFcpbihJM5wkkti7uzxUG/WVNWvjwVKchglcOxfl9d1YfL5ZX72/jYLx2uOe5bC7qhqzNSSmajad7+2eP8x5fu3LJiyb/54XXOjAY51u8rTGUXwOXplRX/QoFUiYcsmc3jtlno8dkrEpu0mQ0c7PJwogEPlTarQTOerS4LqVye3a2uFd5DOezvcK8Q/Nwqdi1mO8vlRg06c6kSM0Maza+8oG9LjRzohq7u/MNTw/zpv1wip0q+dmqcF+1q1ialbyIXJyP8xaNXANa0B/kdZg53ezm7KLHutpqA0l7XTCTNXDTN/QN+zk+EiW/QLdHls7GrxUEkmWt4L8PStroav9lhMXBxqrZSlYFmO00OC4oQhJKZinp597a5MBkFk5uY6d0KwUuA0U2UriqFbujqyPGhAP/jsauFItInr83zzM2FQifEZvDdizP8zlfOlnw/EM/Q16SVubS7LWWNyFOlNozGoAg6vVZanBasJgMToeSKrendnW6G5mI8Hrzlcc1G0xgVQZPTTCSVI7loKJeyd7WMSgRtYM6lqQhCUHXMzm5S1t06DjQ7aHKYMSiCE8OBQmwvls4TTd0qwj3a5+PmfLxkt0mPz0YklWNfu2tLtpDpdX7GRvKZJ2/yqy/o25J7L6EbujoQSmT48okxPvXETSKrZmD67WaklA2vCs/mVT722FX+9w9vbHjsTDhFs8OMQaksRpVXZWFaPGjbP0UkCw9+OqeuMBh+hzYxKhDPrJi94LObmQgliKZyNZdTyMVrHuta27y+HpFklnu6PFjNBs5PhGl3W8mpKn1NmkKIlJKRQIKDXR6uzcYKccw2t2WxZ9mJ22LkxLIEwsmRIK0uC4MDfi1GKwAhyOZURhbiheTCVhg5AHWT+oVXMxpI8Px4mIPdWzf2UO+MqAP/fGqc3/1qaUmkHc0O/vC1++ny2QjEM1ycjNDpteKzm1mIZ3jFXW2YawiMjwUS/N5Xnyv7ARICXrijqS6adfvaXVhNBmYjKbp9di5Nhbmrw01ehVOjG8fNBNDfbF9T3lIpd3W4VtQBFmN/hxvQlJWFEHUxOD67if5mx20xFrHeoxArYU+bk++898GGXLvh6iU6Ggvx9Uswbs7H+bXPnSj5/ntfsYd3Pbij4hFweVXypeOjfPTfrmyY7getpafdY+XiZKQuRg5Wzn2YXBR6PL5Y0Ht3p5sbszFS62gzSSCcqM6zG+z3IdE8SYG2dczkVdLZPB6beXE8Ypxev51MTuXydKQuJSVmo0ImpyKElkw4Oxqqaeu8WRRr7m8kdrOB3a1OHtzbylvv31qxXN3Q1UhelXzqiaGarvGx715ld5uT15QhPZ3O5Tk/EeZ7l2b51rlJxiqYQpXM5hmej69reOrJhckIXV4buby6roLvrlatidxu1jwtCSC1+bHT4RQLq2JerS4LLS5LwaBuxHCdNPWWONLjZSGWxmY2FLo0Dna5y5ZO2grWKyNqBC6rkf/yqn1bHptbQjd0NTIWSDBfB6mbP3r4PD+6Pk9/k51mpwWDIlClJJNTmQqnuDQV4dx4mPlYuuqYVjCR5a4OF06riXa3dU2BaiOYCCUZ7PcXNXQDzY7FhEWq5NZVCGhymBlocZBI54in8/gcJs6ONV4htxTBRIbrq8o8np+IbJo4aDW0ui1Vy+JXQzSVw23dPuZl+6zkNuXfLkzX5TrzsQxffHa0Ltdaj6U4ltWkMNDs2JRP+eVyT4MDfvKqZDaSYmg+vuH9pYSFeIaFeAYhYLB/6+JMoLWNlZIhnw4n6ypBX0/OjoVqHuZdLs1OC29/YR+vv2fzqg02Qjd0NbLVhZDV0Om11jS/oVKuz8bY2+bCYzfWVE4iF8tcDvd4OTu2+cH/g13udWctBBJZTKlcIYa3nVAlDe3UMS5eO6dKPv7mI7xw5+ZMuyuXre+BuY1Z8kxuN5ocm1vADNDusaxQ0qiFs2Mh7u311uVa5WJQREGtYz2yqmRXi3PD47aCRlZY9PrtXPnTV/OD33vptjNyoBu6mpiNprCaK8uUbge2oqQokcnTXoOe2WpuzMVodprrdr2NONrnK7tv9+JUhF0tDjy27THqb4lyMvPVItE+DAaaHQ27Ry3ohq4GOjw21O0ggVEBg/3+ivpWK8WxzPAvCXn2N9kJxjP0+OonnBlO5jAIQVcdr1kKu9nAqZHKttzX5+L4HY0zxEubUL/dhMW48ZZ0f4e7oeGKeDrXUENaK7dfgGkbMbQ4sPd2IliHORDrcaDLw2ggTqvbytXpKIMDPhLpPOcnI7giqbrGr2aiaTyLMw4amVSxmw0c6PIgpWQ+Vr40fIvLUtO6jAq0u21k8yo+hxmbyYBhMRZ2YTKC0SAIJLI4LUYOdbsWjZ9ACG1amcNiLBifkYXGJp1mo2m+9dwkb90m5SSr0Q1dlcTSOX7vq8+RbIAiRDFt/9XsanUyuhBHrXBWaiMMnckg8DvMzETSSGAqnGYqrJWTLI/LRdN57u50l9VfWy7hZA6fvbFb2PlYhvmY5tE5KghVnB0Ncv+An+uzsRW1gCZFcFenm1xeRQiBffGakWQOBHisJoYX4izEM4wvZnCL1iEuOlCxdK5uqsu1cGI4sG0Nnb51rRKH2cB7Xrar7tcdHPBhMxkwGcSKGFSby1LYonR4rEgpaXZZUBSFLm95sa9DXZ51hx2Xy46WlXGYwz1ewoksgwM+IhtsX0yG+mf+hhcS2CrsKqmWeCaP01LevTJ5ybNDAXwOM11eKwLY1+6kyWXh3HiYi1NRLkxGmA6nyKsSj82Ix2biykyU2Wh60zsZauXkcJBvnpnY6mUURffoqkQIwdfPTGA2KGTqNPRjV6uTc+NhjIpCj8/OzcXhw71+G+FkliaHBZ/DRCqrrtAly+Qlx/p9hXWpUmIQ2tBiRWiJgGaXuabxdSaD4ECnB6NBcH4izKEubSaoqsqCN1FOVjWRybOnzcnNuXjJUYGVsKvFgd9p4WID446r2dvu4tRI+eUt12djWI2CQ90eJNDiVOj22jAogtmo1qa2fDDPYL+fYCLDjbnYtlBBLpeJUJLP/niI193TWdhibxf0pv4qCSUyWE0GTo0E+eLxUf7l3FTV1+ryWmlyWAglsw3V7qp2hGGX10ZbnSvr+/w2vHYzFqOBvJSEk5nC2MVyOdTl4cJUZNM9n3oO6FmPuzpcJDJ5JoLJunwobAaf/41BXry7ZVPvqTf1NxDvYlzoRbuaeWBnE+99xR7e8ulnyp6x4DAbUBTBzmYniWyOcxONb2mqtnXMYzPVvX1oJJCkdVkb2u5WJ3d3upmLpgknM6Rz66/1SK+Xc+PhLdneldKbqzdLXSx2k8LhXi8GIbgxF6tL+KFRRFO1yds3Ct3Q1QEhBLtanfyvXz7CWz7zzIYqFgLNeEyGUwWl382gmrouoyIYDTQmY6cs0+hbPknKZBAc7fcxH03T5rYSz+QwKQoGRaAIQTqfZzKU3LIYVjSVq2lCWqUksionFz1Ig9BKRTJ5FZvJwPOb8AFZCe/54mn+9q338cCuZnJ5teAQbDW6oasjL9zZxFvv7+Pzz4yse9zRTdr61IrdbKDbZ2tY/VWpMpNsXhYe7Horj9SD2WiaJoe5rjNYyyUvWSEcsKfNSS4vKxqS00hUCQ99/hSAZoj/+JUYDVuf89z6FdxhNG1Qrd/sMJPYYO5Co5gKV9Zs3ultnJGD4mMWbxd6/fZNN3LFuDoT4+Z8nMEB/1YvZQ3JbH5Te6rXQzd0debdD+6ky1u8Wv9Al5tgMovFqGzJH6bJoLCj2cGBLje7Wjfux2xqYGX/XR2ubfMQVMLhHi/39nqJbTMjfXwowMGurZMqL8Wj56tP0tUTfetaZ6wmA7//s3v51BM3mItmCCbS2EwGBlqcBOIZ8qrk9GhoSwzd6m1gh8eK0SDo8FiJp/M4LUYyeRWDIrg0GWGugTGoWqTjtxKTQRBKZkFqsTIh4PJUhG3g3HFtJkqP31aRGGuj+dqpcR7c28p9fb4tXYdeXtIAQokML/jz75HKqliMiib1vUpqe3erE5/DTCCe4frs9vNsXrSriR9fr4/cejG2Smqpnhzp8ZLO5asq2WkUJkVwsNtTyJIbFIHTYqTHb8NuNiIXxVxXy001Mt6oCPj0247y03e1NeT6DS8vEUL8CfAGQAVmgXdIKSeFEL8P/Mqye9wFtEgpA6vO/xzwILD0W3+HlLL0rL7bhHAyuzguWptnAGvnCSzPMvodJna3uta0Cm0FJoNWaNzojOJ202urhkx+exk50GSiTo+GONbvQxHaeMZwMkt4YmXHysFON1Jome+5aJq5aKphcvBL4zIbZejKodb9w0ellIeklIeBbwN/BCCl/KiU8vDi6+8DHl9t5Jbx+0vH3glGDjTjVkkPbCCe5dmhAN3+xitxbMQ9PV7cNlNNXRTloCha0/vtyl0drkI/73bkxHCQZ4cCJTsrUjmV8xOaPP9UOEVObezwnK3ulKjJo5NSLjf/Dig6OPLNwJdquc/tRl+THbfVuGbG63q0u61c2WLv4J5uD6eGg5sy/fP8RASBNinMaTGSyuaJpnLbpkxiPZYGda839Hq74yjSr3txKsrhHk9D5nG0bvGHWs3JCCHEh4G3oW0/X7bqPTvwKuA961ziw0KIPwK+B/xXKWXRj0khxEPAQwC9vVs7Om0j0jkVp6UyQ9fbZGd6i9WK07nNHXEsYY2SyUCzg2AiQ6gMNd/NxKgI7upwF7Tp6rXz7vLacFmNpHMqTQ5zIZabl5JkJo8QEIhl1p2iVg2lvLdsnfq2l2M2Krx8X2vdr1sJGxo6IcR3gfYib71fSvmwlPL9wPuFEO9DM2gfWHbM64Afr7NtfR8wDZiBTwF/AHyo2IFSyk8tHsPRo0e3dQblqevzhRmn5dDsNHNpG0yPujYbo81lqftDVS47mh10eK0NzfZWg8dmostrq3sXwmC/f8UktlLadXazgbs7XVyYrJ/HbyshNzUaSNZ9iM4DO5voa9pa5eENY3RSyldIKQ8U+ffwqkO/ALxx1Wu/zDrbVinllNRIA38PDFb6A2w3pJQ8eW2+onPsZuO26BHMq5L+TZbCtpkUYbxhwwAAIABJREFUDnV7ONTtocVl4fxEGK99e0iQW00KR3q92ExKxWMMD3S5Odrn41CXh2LKVIMD/rLHTSYyea5Mx7i/RElSNZU6p0eCdBaRto+mcpwcCXK031eXuFpfk52HXrKj5uvUSq1Z191SymuL374BuLzsPQ9aRvWt65zfIaWcEkII4OeB87WsZ6v57sUZPvLo5YrLRUwGgQLc1+8jlc0zFU6hqivHBG4GBkU0REh0PfKq5NyqUoe97e6y5zM0ih3NDpwWI2fKFDNodpppcpixmgzEM3kuTUUL20OPzcS+dhcXJsIoBsFd7e6KPficqmnb7Wp14rGZEELrmVYlnJ8IYzEK+psdpLJ5bCYDdosRJEhk0XbDnMq6g8xPDgc50uvlzGgIl8VINF3dB/F9fT7u7tj6QuZaY3QfEULsRSsvGQHevey9XwC+I6Vc4Y8LIR4B3imlnAS+IIRoQfs/O7vq/NuGK9NR/uX5Kf72hzcq1qbr9dvwO8y4rabCH6TFqGBURN3VeDfirg7XGqPTaDJ5SYfHytSyrX5qMTa1UYmnySAwKIL9HW7OjoXqqt22lBTpb7IzF00TL9G21+qy0N/k4NxEqKSqSDiZLcyiVQQ1zaUt9SGaZm28E7T42OCAn/PjoTXJE6tpfVfw6lSEg10eLkyGGez3cbyK/uyvn57g7k4Pv/FTAxWfW0/0guE68PSNBX7lM89U9aCtFw9RhPZ+NX9glTLQ7GB4Pr6pyQgAs0FwuMe3Zht3/8DaQdVmo0KPT/tgmAwlmYumySwWud7X6yWTV7EYNfmrYDzDXCzNnlYX5yZCpGrIkLa7rbR7rJwdC+G0GPDZzbS5rZwZC902KsAui5E97S5UKTEqgng6V7IG8Fi/b7F4WK7IwN4/4Oe58cp/ly/f18pn33GsluWvi65Ht0nc1eHCYjSsu+3b1eLg+tzKYLPDbFhXaFOVcGIkyP4OF06riXAiy5WZ+peg7GhxYFJEXY1cX5OdSDK74SzUTF6SL/Jhe3EqQofHSqvLgsVkIJ3NYzEZkFKSymrZ4cyySv5Tq7aYZoPAZjZwfDjAkR4PZ8fDG3qIxdjf4cZmUjAYFAb7/VyfizIWTK5QBL4diKZznCozwXB5Olo0ZvzsUID+JjsWk4ErG8w0Wc4zNxdQVdnQAdoboRu6Grg+G+P5iRDfODO5YWzLYzNzpMfEmTGtz/X4UIC97a4NBS2lpPDJazEq3D/gR0Khb9aoCG7OxXBYjbS7rSiL4ozltPMIAcf6/ZwYClRl5GwmA3/+iwe5ORdjeCGBx2Yik1N5YFcT9/b6+NC3L/LYxZkNr3NlOoLdbFih6hJN5Yimciu2tMs52udjMlQ6s53JSzJJ7WE9MxZmV6sDn92MlHBlpviDXIzpSGrThDa3C8Z1DNLwQqIsQYjl/OFr92+pkQPd0NXEN86Mlz3uMJ3LMxpIaFusXB6vvXLV3nROLRrfMRkEsVSOq0ktfmMxKhzqdmMzGRiej5csF1nycMo1chajwkMv2YHHZsKgCO7r83Go21vy+I+/+Qiv+/iPVrS7FWNvu4vzFbYeVfrgaDLtmkftd5hxWY3sanXy/HgYj81EIpsnWSwOdxuGdmql1WVd1xOfi6YKH9alMBsVbCYDD71kB28e3Pq6V93Q1UC5GTmA85MRzAbB1dkYsXSurs/Pau8tnVMLSYXBAX9JQ/fAziZS2Ty/+dKdnJ+M/N/2zjw+zrLc+9979kkmmckkzb42DemSrkn3FqEsraAsAoIoq4pUwBffV0QPKirKexC348GjclRU5HhUFOH0ICDbEVRauu+1W9I2TdOm2fdZ7vPHM4lZJsns8yS5v5/PfDJ58jzPXJ87M9fc931d1+/iz38/xwU5Dm5ZUUJ3v4/5BU7S7WZSrSZONHdTlplKcWZKyHbZzEY2fWoNbxw8y1uHm3hh12nynDaWFGfQ2N7L346dZ05ueliNZkATm2zujDzXbmCG1tzZT5rNxPmu/qB7ggsKnQkPzuiBlAm6nLX1eNlyvJmqgvSgX1BXLsjj2x9cCIDVlJjubBOhHF2ESCkpyUzhr0dDV/jo90kWFKbFNBlzIrbXtbCszE1DW88o+Z7qkgzuuWgWNosRv1/y58PnKMywMys7bdR9yiLMr7OajGyoyuOyubl89eqqYblZ9a091Lf08NGfvRtW+sKZtl4qstMYmKFFSt2Q/dHNx5spn5HKqZaeQSGGfm9yBFITRZrVyJy8dIQQnOvUVJO17ZDQEvOCObGVMzP5lxsX6UJVeCgq6hoFvR4f73n8DRrbQ59dlGam0NjRF3yZFEesJgN+v8QzJEpoNxt56f61Sc9a/+yzu/jN1lPjnpNmM5FuM7OhKpe/Hm2isa0vLnmGAwoeGSlmWrs9QZf17oC81mQjzWaiLCsVKbX/fUt3/4TbCuNhNQkqstPYG0hryXJYePLWGpYUJ1Z7TkVd44zVZCDVYkLLYgqN8139ZKdZqUtwL4Q+r3/Y8swg4KZlRUl3cn6/HJwVm42CXKeNk809WIwGzEbBqllZHGho50vvm8vl87RKxBf3NPDJZ7bHxZ6OXi9z89KxmAy0dLfiSjFzQU4ap1q66e734bSZaersoyLHQd357kkjN5VuN9E/ZEsjFvR5JX0+P4uLXJgC5R9FrtC3NhKJcnRRIITgzjVlfOEPoRV0OO0m+n0y4U5ugKEF24uLM3j4/fOSYsdQvH7Jc59czdbaZlaVZ2E0Cp56+zg31BThl5J8l53W7v5h3aSWl7kxCOLS3HlAhdlk0BSYsxyWYZvuA2IDhxs7cadaKC9MHUz0znfZsJuNw5qL64XZuenjBg8iZaicl9EgBh2e3lCOLkouDkOVwSgM5DvNSfsgSDTl4G/dsEg3UuYWk4EZaVbeOz9v8Nh9l1QMO2dky7yMFAuLizNCzguLBK8fGtp6ae3uH6wxHRmsaO7qp7mrnxJ3ChazgSONnUg0GaeuPi/NXf3ku+wcPtuZ0MTiRUUuOno9uFIs9Hq8pFhM7E9Ahc0D6yvJiGOfkWhQji5KHBYTWQ4rTSFEAZu7+2nu7qci2zHu3ojdrHWvr8h2kGIx0tbjiUkjmcxUCz++LX4Z6onCYBCsLs+Mq6MboMejpfSYjVo6zb7TbaMqA+pGJH3vO92O1SSQUku+zU6zUpKZwta6loRkq7R09wdWDYn9Qk2z6ded6NeySUKPx4c5jOm6027iyLnRTstuMTI/34nH58dk1KoUOnu9g8uiBQVOdkcpE9Te66XX48Vm1v7tm4+d53xXPyaD4LK5OQihz2VHMC6fl8v3Xj+SsNfz+CTb6lqoyHGQZjXR0t3P8aaxtyD6vP/waGc7+jjb0UdlThqd/V7qW3rId9nITrNhMggMBkGvx8fe+raol+OxeJ9EysPP76Orz8vH187U3XtJOboo2d/QFtaypCI7eHrJvPx0ttQ2YzIQVNRxd30b5TO07P6OPm/QEhyrycDc/HTMRgO9Hh/d/T7cqRbae/oBwZbjzfx2Wz2Li1y8dbiJx18+OPjBWj8vh69eXUVO+mjpHj3yb28mzskNZeie1AU5DlKtJvadbg8pKHGosQOzQVDgslHf2juqsiPSwvmhNHYkT7zV65c8+uJBbl1Zis2sj/y5AZSji5IVMzO5ZE42v9pyMuJ7mI2C/YFv4fE+L9renrYcqcpPp6Gtd1gznfIZqRMmMX9xjMDJy/saubgym5t0kMUeCnoQ5xzYTihw2bBbTKTbTBgNgkNnOsZUl/b4JfVjlK7Vnu9mTl4aByKU1F9c7GJvkmZzQ3m3tpm1FTOSbcYw9LEjPYlJsZh49Nr53L6qNKTz23pGl9YsKHCF3X9g7+l2+n1+FhY6WVbmZmlpxmCia6S8uPdMVNcnkiuGBC+STX1rL0fOdrL9RCvv1rZEnFx9tqOPAw0dLCtzj6kAPB47TrQyMyu8OtR48MlntvP9N47ERZY9UpSjiwFCCL5w5RxWlWeOe97yMveoIMTiIhfbTkS2XOno9bLrVBtbjjfzbm1L1NHcDfOCKebrkztWl4X85ZJoDjS0D6bARMKW483Mzh1dnTIRhS4757uSP9Pt6PXy+MuHePTFA8k2ZRDl6GKEyWjg+zcvocAVvGXhrOxU3g0ind0aZIaXaMxGwaPXzueDNYXJNiUsHtwwmw8v199Su9+nqQHnj/FeCIWdJ1pxp1oocaewqMhFdUkGleM4v7l56bT1eMYU/0w0QmglhnpBOboYkpFq4QcfWTJKa39xsYv61t5REbW5eWljNkRJJA6rieqSDN3VJ06E3WLkq1dXsbBobAWVZBJNH5ABKa665m52nmxlW10Lh850UBpEVGFJsQub2UCvjmpzr1lUwJU62l6YXO/sScCCQhdfv6Zq2DGPzx+0tjXU4ul409Lt4bGXDjIZ656NBsFzG1dx68qSZJsyiu7+2Dc8yh4RFV9amsH2E61sP9EakgZhonhuRz2PbFJL1ynNjUuL+MCSAvJdNk2vP4iUjc1siLjhSDw4dKZDF5HMSDAYBBsvKmfmjOTW7Q5FAJmpsW/afPRsJwsKnBRm2JnhsLJLxzJSL+yq51iQnNFkoBxdHBBC8PD755JuM49ZXzi/wKmLZesAH11TNmq2MJnIc9q5/9ILxlXHTSQS6Or3srzMHVObznf1s7u+jXyXnXyXTdeiAk2d/Vz67f/hjqe2JNsU5ejihdNu4bHrFgTtjZnvsrEjwkhrPFhbkcWHV+hvUz9crlqYz8fWJr+H6AAdvd7BFoXF7tipeszJ01on6nk2B9q2Qr7LzupZWck2RTm6eLKwyMX9QwrUTQbB8jI3gvETg0NlWWkGWQ4L1RPof1VkO0gdIy9rzaws7r/0At0owUbLnWtK+fSlF1DgsutmdnfwTAfnO/tYVuomz2kLq2QwGCkW07g9WfXCl98/l7cfXKeLLx8lvBlnPD4/13z/L3h8frw+OdgvNFYsK3NT29RFisVIv88/WFY0kJrQ2eflyLlOpNTO7fP48PklXf1eHr12ASsnyP2bzEipNW/+4h/2xqV7WqS4U8zMyknjyNnOiAU8g0m/64ncdBubPrWGPq9/zJSrWBGK8KZydAlg18lWbntqy6CWWSwozLCT5bCy86RW8jWgzza/QMunMggxqK02wOzcNNbPy0UC763KZU5eeszs0TO9Hh+PbNrPM5tPkG4zIYku9SNWmAxalL6zzxu2Os28/DRSrWberW3Wdf8em9nA7zeuZm5+/N5rytHpiM/8dhfPbhtfLjxUFhQ6OdDQHnY6wab71lBV4IyJDZMNn19yoKGdqgInvR4f2+taONXSQ2VuGk+/Uxez/02kLChwYjQI6lt7OBtG9Huiblx6QAj4zOWV3HPxrDjdX0mp64bPrq/kxT0Nw3qXRkJNSQaHGjvCdnIz0qyURliDORUwGsSgk7eZjawaskG+oNDJslI3X920n84kpfwMlVaqyHHg8fpHzciDseV4M5U5DmrPd0dd6xwvpITHXz5EvsvGtYuTU32jghEJIjvdxr3rovtGW1qawda6loiWXec6+rjo8Tc53Tq5OswnAiEEH1xaxOZ/uoRLwlCMjheHGztJnaDl4FAONXZSkeOgMictaJRfLxw9m7x0KuXoEsgdq8rIjSJX7Wy7FrmLtIawqbNP13lXySbVauLxGxYyv8AZdWQ0Wvad7qAi2xGyHXvr2znU2MGsGQ7cqeY4Wxc+QsD11cmrpVaOLoHYLUbuiWJWV9fczZbaZrxRyN/o+RtfD7hTLTx/z2p+cedyMpPc/+DYuU5K3CnUlGSwtDSDmpIMlpe5yU4bu+KiqbOP5i4Pi4v0tRcrJbx56GzSXl/t0SWYG2uK+OGbR6kPYwm5tFSbwXl9EqNBcDTCspqNF5VTmBHfUP9UwGAQrCzP5JVPX0h9aw8f/flWmjr7Eh7d9Ek4MkRsdYDFxS68fj/NXaOj+INCrDqTMgeoyAlfeipWqBldgrGYDHxsbVlY12ytbaGrz8eOk61srWuhJcI0FYfVpDstfz2T6bCyoNDFK/dfyJxc/aTi7DjRikAMzvDSrKPnK3VNXbpJmB7g9qe28PO/1nKuo49eT2KVVqJ2dEKIR4QQu4UQO4UQrwgh8gPHnUKI/xJC7BJC7BNC3DHG9dVCiD1CiCNCiO+JafBJvGlpMU576PsoEmho62FZqZt5EeQjGQ2C9y3I464Lk5+hPhnJSLXwIZ3p3p3v6mdrXQubjzcHjaa3dHuoLslI+l7jUDw+ycMv7GPp11/lhh/+LaGvHYsZ3eNSygVSykXAJuBLgeP3APullAuBi4BvCSGCbXr8APg4UBF4bIiBTbrGbjGGLSvU0u1hS20z1hD7sdaUZPCBxQX8x8eXs/vhy3ni5iWYJ5nenJ64obqQe+OUBxYNZVmpHAnSOlOi9aFdWOgiI0V/wYlI5eYjJep3vpRyqAZRKtoYE/iZFpihOYBmYFhehBAiD0iXUr4jtczlXwDXRGvTZGDjReURLS1C3SdaU5HFt29cxKryLFKDLG0U4WEzG7l1VQkPbpjNe6tyxwzqlGWlcufqMn63cSWrZ8W3vM4gtEfPOMvAga2OyhwH7kBwRQ+zvJf3neFrm/Yn7PVi8gkQQnwduBVoAy4OHH4CeAE4DaQBN0opR4YLC4ChKemnAseCvcZdwF0AxcX6WkZEQorFxMWzs/nT/sawrjt0pp2SzJRAg+Lg5KRbWRboLq+IHdlpNjZeVA5odbS/fKeO322vJzfdxoqZbjIdVtZWZOFK0RzKMx9bwUd+vJm3jzTFxR6jEGSkWAilUfWhxk4EWlvNXo8v6v4i0dLn9UfdwzYcQnJ0QohXgWCdUx6SUj4vpXwIeEgI8XngXuBhYD2wE1gHlAN/EkK8NWIGGDJSyieBJ0ErAYvkHnrjqoX5YTu6bo+f1u5+FhQ4sZmNbAnSh+L2VWWsKk++NM5URgjBLStLuWVl6bjnfeP6BXz61zvjUoDv8UuaOkMvF5PAvtPtGIUWlE129WeRO3EZACEtXaWUl0opq4I8nh9x6jPAdYHndwC/lxpHgOPA7BHn1wNDswgLA8emBetmZ2MJcc9tKG09XnbXt7GltpkZI3KqSjJT+IQKOuiGfJedH3ykOmhkNBgzs1K5eXkxv75rBf/3sgtYVuoeN1MkOz18FeMshzXpTg7gX18/kjD5/qiXrkKICinl4cCvVwMHA89PAJcAbwkhcoBK4NjQa6WUDUKIdiHECmAz2vL3X6O1abKQajWxvMzNW4cjX9qUZqYMk0C3mgx6TKGa1qRajVxfU8jTf6vDG2S9ZjMbuHN1GTcvLyY7zTb45bd8ZiafuqSCbXXNfOLpbUE7fJ1u7WVhkZPmzn5OtoSWm3mus4/lZW78UnKqpYeGtuANteONQYDXLxOyZxi1eokQ4ndoTswP1AF3SynrA2kmPwPy0CT0/1lK+cvANTsDUVqEEDWB8+zAH4H75ARGTUb1krH4wZtHeeylgxOfGKDEnUJ7r4eWbg8Wo1aovv1E67BzKnPSeOjKObx2oJHPbpitghE6YfepVh747e5Bbbwsh5WvXDWPy+bmTDizb2jr4bafbhlTzmlJsWvU+yAUstOsdPZ5oxabiJSlpRlcNjeH21eVRbS6ASXTNCnYfOw8Nz75TsjnF2bY6fX4KMlMpbmrj+NN4ytcfOqSCm5ZUYKUklSrSTm9JLOtroXvvvp3TrX08Is7l1EUhsS61+fnrqe38frB4aVUJW47dc2RizXMzUvn8NnwFXFiyVUL8/nehxZHdK2SaZoEzC8MrybRbjZyqqUn5EbF33vtMN97TdtZWFTk4rMbKlWgIolUl2Tw9EeXI6UMu0rFZDQE7T2RlWaLytHtb2inpkRTxgFNHbgix0FpZiqHGjvYebKVfq8fp91MWwwarpsMYtQS/oVdp+nx+FhSnDEY2Y4lytElGbvZqDUf9oRWqH/4bCcFGXbqQ9yPGcrOk63c+pMtXF9dSIHLzn1D+lkoEkskBUBSyqDR2zNtPSwtzWBPfVvI76ORHD7bgdEguLAii0euqaIwQ3Oofr8c3PMVQvB//nMHz+88jcNq4s7VpfikpPZ8N4cbO+jq8+Hx+cl0WDnQ8I/kinynjbn5Th66cg5n23sxGgQf/fnWUU7zT/sb2VrbzO2rSrGP0eMkUpSjSzJCCGY4rJxs6cFsEOS6bLR2e4JqzjmsRubmpQ9+80aC1y/5z3dPkmoxKkc3yRhQSR5JfWsvHX1e8pz2iFtotvV4efTaKm5ePrxixzAiMfrWlSXML3CyoSp30BkOMOAUezw+nvzzMRrb+5ibl8ZNy4oHq3IGKiI23beGLz2/lzcOncNiNLCwyEme087MGan0eHzK0U1F5uWnYzUbcNotbKtrwWY2sLzMzfGmLgoz7NSe1968xe4UttTGpk1iV7+PvfVt01ZafTIyUnHYbBSD+2qVOWm8G+V7I885cV5bdYmb6pLgyegDTjHFYuL+Sy8Y9z5F7hR+fNtSXj3QyPIy92CSdbxQjk4H1J7v5lRLD0cCCqy9Hv/gEmVo/4AZjtiJZjrtZoozY9drVBF/huq5FWXY6e73UZ7twOPzR+3kDEIrc0skRoNg/bxgdQixRzk6HdDY3hvS3kpKDKfzbz94MWk2/RV7K4Jztr2Xb75yaPB3m8XIyZYezseo4uJDy4qndOtLJWeRZLr6vCHry8Wyccu/vHp44pMUuiE73cbvNq4abFEZyxRbi9EQl0innlCOLsmE01g5I4bS3pt2N8TsXorEMC/fyY8+Uk2hy87hINJMkVCamRJQnp7a2xhq6ZpkToWRJrLleDOLilyDTavD4eLKGRS7U9h3up3a81386q4VYd9DkVwa23v4/388gMVsiKpW1WTQup5dMjubdbOzp4XqtHJ0SeZcGM2KQcuFG4jIhtPouKbUPdhA2O+Xo9IGFPrn4JlO/rj3TFT3uHJ+Ht+5cVHE5VaTFeXoksz+0+GrVg1EZIvcdrIdVrZNUOMohFYVMYBycpOT37x7MqLrbllRwnXVhdS39HDF/NxpMYMbiXJ0SaS738uf9kf+DX2yuYeTzT3ML0hnT/1wh2kxGlg3O5uNF5VTmpmKU4dy2orwuK66gNcPnh1XUTgYtee7WFTkGvZlN91Qji6J/NsbR2kPUgERLu29XoxCa48H8OQt1Vx4wYyE50Up4su62Tn8/pOruPFHfwv6vrGbjVxfXchVi/L5yVvHeWnfGYrdKRMm704HlKNLMH1eH9946RA7TrREJKszEiHAZjKS57JT39rDuspsVpZnKic3RZmTl86SkgzePHRu1N/uXFPKA+s1bdulpW6OnuskN92mFGtQji7hPPDb3byw63TM7jc7N40DDVqKygPrKwcDDoqpx6v7G/nua39nb/3ofd0P1hRy78XDa5fLZzgSZZrumV6hFx1wqmV8/bhwKHDZkVKSk25l/bycKZ3ZroCV5ZmkmEfPTWZlO/ji++bGvBB+KqFmdAmmICMlJkvWApcdd6qFh66cQ01JBibVs3XKk2o18dQdS7n1p1vYNkTBZmGhS5XzTYD6dCSY21eVxuQ+eU4bT3xoMStmZionN41ItZp49u6VPHTFHBYGoqg5ETTImW5M+U+Ix+enuSs0Nd5EsLDQSVVBelT3qCnJ4BvXL6Akwd3OFfpACMHHL5zJT26r4frqQj5x4dSuU40FU37peqatl43PbKOn38cD62ezoSoxsjBjYTIaeG9VXtAN5VC468KZPLhh9pid4hXThyyHlW/esDDZZkwKpvyMrsidwheunMvRc13c/cttPPriATy+2Om6hYvH5+dMhO3lvnLVPP7pijnKySkUYTLlHR3A5mP/0Ox68s/H+MZLBxPWOHckrx1o5Ol36sK+bm1FFrfFaH9PoZhuTAtHd9OyomG///tbx3nspUNJmdk5rJFFx/7f5ZUxtkShmD5MC0f398aOUd3rf/g/R6n52qts2h275N2J6PX4+PXWyAqzP/bzrdzzzHZe2RedeoVCMR2ZFo5ubcUMnr17FZYRaRhtPR5KMxMTuZRScs33/8J/RVgV0dTZx3/vaeCBZ3fT2q2fKLJCMRmYFo4OtMbB5dnDS2I+sKQgYV2whBCYjKEHEVbOzGRmVuoo5zy/wBn3jkkKxVRjyqeXDHDsXCfLy9ysmOnmuR31eH2Sh66Yk1AbKnPSJ0wrMRoE1y4u4HPvnU2Ww8qZtl4eem4Prx3UOkAVZ6bQ7/VPO+FEhSIaRLKij9FQU1Mjt27dGtY1Xp8fIQQGAX4JfikHm+omik//eifP7agf8+/zC5z86JZq8l3D+2t293u56xfbePtIEwAfWFzAtUsKWFsxI672KhSTASHENillzXjnTJtpgclowGgQCCEwGkTCnRxAlmPsJefKmZn84Z7Vo5wcaA2Bf3J7DTcvLwbg9zvq+eQz23nn2Pm42apQTCWmjaPTA3bL8J0Cg9CUJz6wpID71s0aNxHYajKy8T3lg9HjBzfMZsVMpVaiUIRCVHt0QohHgKsBP3AWuF1KeVoI4QR+CRQHXuObUsqnglz/JpAHDLTCulxKeXbkeVMFv3/4NsEn3lPOgxtmh3z9uc4+BCCB2qYuvD6/KuhXKEIg2k/J41LKBVLKRcAm4EuB4/cA+6WUC4GLgG8JIcZat31YSrko8JiyTg5g/TytznZhkYvHr1/AdUsKwrr+j3saGPCVP377OAfPhN4TVqGYzkQ1o5NSDg0hpqJNNgj8TBNauyEH0AzErs38JKWqIJ1XPn0h7lQLWY7wpXXuXFPG3vp2Wns8ZDks9CexZlehmExEnV4ihPg6cCvQBlwcOPwE8AJwGkgDbpRSjvWpfEoI4QN+B3xNTsYwcIgIIbggJy3i6/OcdtV4WqGIgAmXrkKIV4UQe4M8rgaQUj4kpSwCngHuDVy2HtgJ5AOLgCeEEMFE2D4spZwPrA3ECqHQAAAE/ElEQVQ8bhnHjruEEFuFEFvPnRvdGEShUCjGYkJHJ6W8VEpZFeTx/IhTnwGuCzy/A/i91DgCHAdG7bpLKesDPzuA/wCWjWPHk1LKGillzYwZKn9MoVCETlTBCCHE0LZDVwMHA89PAJcEzskBKoFjI641CSGyAs/NwPuAvdHYo1AoFMGIdo/un4UQlWjpJXXA3YHjjwA/E0LsAQTwoJSyCUAIsTMQpbUCLwecnBF4Ffj3KO1RKBSKUUQbdb1ujOOngcvH+NuiwM8uoDqa11coFIpQUNmmCoViyqMcnUKhmPIoR6dQKKY8ytEpFIopz6TUoxNCnEOL8o5HFtCUAHMiQc+2gbIvWvRsn55tg8jsK5FSjptcOykdXSgIIbZOJMaXLPRsGyj7okXP9unZNoiffWrpqlAopjzK0SkUiinPVHZ0TybbgHHQs22g7IsWPdunZ9sgTvZN2T06hUKhGGAqz+gUCoUCUI5OoVBMAyatoxNC/FoIsTPwqBVC7Bzx92IhRKcQ4jNjXF8mhNgshDgSuNfYvQhjaJ8QYtmQ47uEENeOcf3PhBDHh5y7SGf2JWv8LhNCbBNC7An8XDfG9V8WQtQPuccVOrLNLYT4kxDicOBnRqxsm8C+TCHEG4HPxRPjXB+3sYuRfeGPn5Ry0j+AbwFfGnHsWeC3wGfGuOY3wE2B5z8ENibCPiAFMAWe56F1TzMFueZnwPWJHr8w7EvW+C0G8gPPq4D6Ma758lj/ex3Y9g3gc4HnnwMeS5B9qcAaNDm1J8a5JiFjF4V9YY/fpJ3RDSCEEMAHgV8NOXYNmqrxvnGuWYfmDAF+DlyTCPuklN1SyoFGQTb+0VAoKURiX5LHb4fUZMBA+//ahRDhdxpKrm1Xo40ZJHbsuqSUbwO98Xi9cInCvrDHb9I7OrReE41SysMAQggH8CDwlXGuyQRah3ygTwHh9R6M0L6AjcuFEPuAPcDdQ+wYydeFELuFEN+J44c5EvuSOn5DuA7YLqXsG+PaewPj99NYLw+jtC1HStkQeH4GyImDbRPZNxHxHjuI3L6wx0/Xjk5M0JgnwIcYMptDm3Z/R0rZqVP7kFJullLOA5YCnxdC2ILc/vNofTaWAm40560n+6ImUvsC184DHgM+McbtfwCUozVnakBbIunFtkGktv4Ke1YfjX0hENXYJcC+QUIev0Ssw+O4vjcBjUDhkGNvAbWBRytaT9l7R1wn0AqHB/aiVgIvJ8K+IOe8DtRMcJ+LgE16sS/Z4wcUAn8HVod4n1Jgr15sAw4BeYHnecChRP5vgdsZZw8s3mMXrX2RjJ+uZ3QhcClwUEp5auCAlHKtlLJUSlkKfBd4VEo5LIIjtRF6A7g+cOg2YGRXs7jYJ7RopSnwvARt1lY78kIhRF7gp0Dbg4hH46CI7Evy+LmA/0bbjP7LWBcOjF+Aa4n9+EVsG1rP49sCzxM2dqGSgLGDKOwjkvGLtadO5AMtMnn3OH//MkOiR8CL/CMqNhPYAhxBi85aE2EfWu/afWh9b7cD14xh3+toe2R7gV8CDp3Zl6zx+wLQFbBv4JEd+NuPCcw+gacD47c78MHI05FtmcBrwGG0plDuRIxd4Hgt2iqnE21vdW6ixy4G9oU9fqoETKFQTHkm+9JVoVAoJkQ5OoVCMeVRjk6hUEx5lKNTKBRTHuXoFArFlEc5OoVCMeVRjk6hUEx5/hdaVbvjTrpKlwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "biobio.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio = biobio[biobio['NOM_COM'] != 'Lebu']\n",
    "biobio = biobio[biobio['NOM_COM'] != 'Coronel']\n",
    "biobio = biobio[biobio['NOM_COM'] != 'Talcahuano']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.axes._subplots.AxesSubplot at 0x7f8b189f1f28>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASMAAAD4CAYAAACnroB1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOy9d3zkd33n//x8p/eiLq2k0fZdr7d4d7VgcCgxcMARJyFHCgmQBDgSkks5uMRxAvmF40LiPI6QXwgJCZzhcEhCQkkxmJjYxrht731XvUvTe/l+7o+vNB5JI2mqJHvn+ddqvu0j7cx7Pu/2egspJQ0aNGiw0SgbvYAGDRo0gIYxatCgwSahYYwaNGiwKWgYowYNGmwKGsaoQYMGmwL9Ri+gEpqbm6XP59voZTRo0KAETp06NSulbFnrvJelMfL5fJw8eXKjl9GgQYMSEEIMlXJew01r0KDBpqBhjBo0aLApaBijBg0abAoaxqhBgwabgoYxatCgwaagKmMkhPiEEOK8EOKsEOK7QojOgmOvn3/9khDi6RWu7xNCvCiEuCmE+HshhLGa9TRo0ODlS7U7o4ellPullAeBfwU+BiCEcAN/AfyIlPIu4L+scP0fAZ+WUm4HAsAvVrmeBg0avEypyhhJKcMFP9qABT2SnwG+LqUcnj9veum1QggBvBH4x/mXvgT8aDXradCgwcuXqosehRCfBN4DhIA3zL+8EzAIIZ4CHMBnpJRfXnJpExCUUmbnfx4FulZ5zgeBDwL09PRUu+wGwEwkxVdeGOLaZIS9nU7ef18fVmPxt4SUEu37o0GD+iDWElcTQjwBtBc59JCU8lsF5z0ImKWUHxdC/DlwBPhhwAI8D7xdSnm94Pxm4IV5Fw0hRDfwbSnlvrUWfeTIEdmowK6ezz55k4cfv5b/uclm5P33bWU0EEeVkvt2tPC2uzuIpbL8yt+e5vd/5C56m2wbuOIGL0eEEKeklEfWOm/NnZGU8v4Sn/ko8BjwcbRdzpyUMgbEhBDfBw4A1wvOnwPcQgj9/O5oCzBW4rMa1ADfEsMyF0vzR9+5mv/5aydH+caZMX7kQCdPXptB/vMlPv2ug3hsjTxDg9pTbTZtR8GPDwAL7+RvAa8VQuiFEFbgGHCl8FqpbcmeBH5i/qX3zl/XYB3I5lT+6fTo6ueokn+/PMXHvnURgB/cmEVRGq5ag/pQbTbtU0KIi0KI88CbgV8DkFJeAb4DnAeOA38jpbwIIIR4rKAE4LeA3xRC3ESLIX2hyvU0KJGZaIr/uLosr1CUQDwDaMbpDX/yFL/61TPcnonWc3kN7kCqCmBLKd+5yrGHgYeLvP62gn/fBvqrWUOD8snmVB6/OFnRtf5Ymn85N47VoOMn+7tRVcm2FjtPX5/h6mSEj7x5J6eGAkjgVVubarvwBq9oXpYSIg0qJ5tTef2fPMVoIFHVfW7PRnnfF48TTmYXvf5vF8bJ5iT//c27qrp/gzuPRjvIHcbgXLxqQwRwYjBAKpujx2td9PqIP4HdpOcnDm+p+hkN7iwaxugOY1uLjXanuap7GHSCfp+XvR1OYqks9/S4Fx2/ORPlyRLjUQ0aLNAwRncYQgju3V55LOdYnxeHWc/xQT9nRkJkciqnh4OLzpESXrg9V+1SK2Y6ktywZzeonEbM6A6k22Nd85xWh4l4OktOlTTZTHS4zQzMxnhxwL/ovKUxIwCnWc9P929clfyJgQBfOzXCa7c38/77tm7YOhqUR2NndAey1kDzfp+X6UiKVFYllVUZDSY4MRhgNpou6f4mg45vnBnj9kyUfzo1yunhQPWLLoPZaIqnrs1wbTJCY3z7y4fGzugOZC6aWvX4SCAOQCZX2Qd5JpLiM9+7wWe+dyP/2oFuN/f0uPn4O+7Kvyal5MvPD/FDO1voa65dm8lkWHPT/v3KFB/48il+9Y3bOdDtXuOqBhtNwxjdYaSz6prFjhOhJL1eK16bkYyqcnEsvOr5pXBuJMi5kSBGncKrtzWRzKj83Ylhnro2w/ZWO594YB+v3lZ9XdLT12f40nODAATjGZ64MsVzt2b5xAP78DVbOdzrrfoZDerDmo2ym5FGo2zlPPj1C3z1+HBZ1xzu9XBhLEQ6q9ZpVfDqrU189YOvquoep4cDvPuvXySRyRU9/vb9HfzZTx1C12hpWVdKbZRtxIzuML53Zarsa04PBdjT7qjDal7CZqp+k/5737y4oiEC+LfzE/yvx66seLzBxtJw0+4wDnS7+ffL5RkkCeTU2uygd7U5eNfRbm7PRLk1E+XGVJQ37m7lF17bt/LzpeThx6/hj6V519FuDnW7i2orTYVXj4UBfO3kCJfHwzz8X/azpYSsYoP1o2GM7jBeu725bGMErCi6VipCgN2kp8lu5IGDnTTbTQCoqlxTCeDxS1P8xVO3APi7EyM4zHrefayXj75lV97lujgWIhBfO9sXTmbZ3eFoGKJNSMNNu8M41FNZVunyRKjiZ25rsXG4x0MkmeW5W3O8+dPf57uXtEZdRREkMzk+++RN/uv/PcnHvnWRZIGrlcmp/FlBVg4gkszyl0/f4tj/eoL//g/nmImkSOfUkndvjzw3yOBsjGyufjGwBuXT2BndYexsc6AIKNfr6mu2caHCrFqz3USwYNfij6X58ydvcsTnJZNTeegbF3miIJZ1czrKV37xGIoi+KNvX+XyRPHnzkbT/NPpUUKJNDajruT1SAn3/++nufj/vYXxYJw2lwmTfvn1Ukoe+uZFDvd4ePW2Jrw2Iya90pDfrRONbNodyKv/8HtMhMprmehrttFiN3FmJFB2/dGxPu+yym3Q3LZ0ViVdZIfishj4sUOdPPLcUEnP2NXm4NpUpKx13butiZvTUbw2I11uCxajjo+/4y4C8TTP3Zzle1eneebGLBaDjqyqYjbo6Pd5+ev3HGmIzJVBzWRnG7zyMOnL984HZmMMzMbo9loY8Zfe9d/qMDE0Fyt6LJpa3kqyQG+TtWRDBJqkyVGfh6G5ONORtQPZAM/d0vrnpiMprk5qhuw/rk4TTy/OyC1k6DK5LN+7Os03z47x4/c0VAlqTcMY3WGMBRMM++MVXbu12cbt2eKGpRiHezyMBeNMlpDlWuBIr4ecKst6DmjV4icGA2zxWDjq86AIQTiZIZzIMhdN0WQ3EopnSOdUpITMCn7qUkNUjIe+cZEjvV56mhpB8FrSMEZ3GF96brDseNECHqsRKN1IpHK5sgwRaLulhV1KJYwGEkX1msaCSYTQ4kUOk55DnU5yqkQIyGRVzo2WHqBPZHI88Nkf8OBb9/Cuo90Vr7XBYhrG6A5iJpLi0RdKd30W2NFmRycEp8pseE1mystWNduNZe+IymEhPBpJZTleEMNSBPT7PBwfLP33C8Qz/PbXz7O308m+Lletl3pH0kjt3yFIKfnTJ64TK8ENWcBqULDoBZlsrqLdys3pKDvb7CWf39dsq2vLyUqoEo4PBujvK69vTZXwoa+c4pP/dpnHLkw0FAKqpGGM7hAeuzDJoy9qPWkui55D3W58TVaOrfABPNzrIZWTuKwmEunKDUQmJzGUkHkyKKIsV6keHB/wc6TXU9Y1o4EEf/3MAL/86Gk+8OWTpLKlG/sGi2kYozsAKSXfODPGjjY7nW4z21sdnBkJMjgX58UBP/0+Lzajjh1tdvr7vHS4zJwaCpBTJZPhJNkyg0w6RXDU5+FYnxeDTtDuXlvmNqNKDm55qSBTAL4mKwe73RzsdtNsX5/BkaeHA+zrclZ07RNXpvnwo6c5NxJs7JIqoFFn9Aonp0p+5xsX+OaZMVKruEAWo47ECi7cng4H1yYjawa+O1xmutwW5mJpBiqM/bxmexOToSSxVC6vSwSwtcWGzagjnZVMRZKkMzniZcakSsVrM5LK5MpyaZfyQztb+KEdzbzt7g463ZYarq5+nB8NYjfp2dpSumtdCqXWGTWM0SucJy5P8f4vV/+36nKbGQuuXCh5T4+bG1MRIqnq3ZTDvR5ODa0dTG53menxWjkx4F9TvbJcDmxx1cRtNBsUXr+zlR1tdna2OXjrvnb0us3lkERTWf7k8Ws88twgQkCv14rLYuCRn++vySjzhoRIAwCeuTnLvi4n/X1e7u5yck+Pm64Kvqm7ljSWNtuNGHQvxYIMOqUmhgg0F60UJkNJjg/42b+l9tmsc6MhDlfYx1dIMqPynUuT/ODmLL/61TM88Nln+ffLU6g1UkGoBV96bpBH5gXppNTGWZ0bDfGaP/oPpsLrN9ygkdp/hfPszVluTi8eRS3QdK4REE9nuTgWXtVNA/JNpVaDwr4tbkYDcQ52uzkxnw6vxUerr9mGUa9wsoRdUSHnRkO0Oky0u8yYDQq5nOT6dJRIkWEB5XBqOEh/n3dRGQBozcaheKbkMgRFaMYa4NJ4mA98+SQui4Gv/OIxrCYdzTYTLqshH2daj963eDrLeDBJKpvjL5++VfScrCppmVdXWA8axugVjJSS0cDyamsJHB986QPW5jAxG01xsNvN2ZHgsvMBTg8H8TVZCSYy+Q9nMqNiNSi0uSycH6lOdP9Ir6dsI1TIdCS1qA3EYtBxpFdTCii3Z62Q4wN++pptNNmNzEXTmA0KF8dCeG1GDmxxcX40tKoh7vVa0enEMoMWSmR4x5//AACP1YDTYsBi0LG73cEHfmgro4EEY4EEP93fg6WMJuBSOD8a5MN/e3rNth6zXqm5+7saDWP0CianSrIlNLVOzX+Ip8IJdrXZsRr1nClilAbnFhu2rc02FCE4Pqi5SueriLHUejOQyOQ4ORSg1WGixW5iZo0hBKux0JdXyFQ4xVQ4RZvThE4RzEbT9DXZkEiuT2k70X6fhzMjwTUbiwPxDIF4BoCrkxG+eXY8f+xLzw/yhfceYXtr7ZQ2/+LJWyX1F4aTWa5Ohrmrc32KOhvG6BWKlJKPfetiWWn5iVCKiVAKh0lHk83IXGx1sbLCnczFsRD39LhRFIFBEYSTWS6NlyY54rUaqjJkxTDpFXa22RmYiRGtIiu2FoXqktemIgihqRQoAl4Y8FNtfmhoLs5b/vQZnvjN15U9QWU6nESnCJoKXC0pJfFVpHkLMeqUiuKLldIIYL8CiSQzvPUzz+S/oQspRfkinVMJJzNlPVOVmit3cjBAMqNiM5XmWuxud2A16VctO6iEvmYbl8bD6HQCe4lrqQVSwosDfi6MhnCaDTW5Z06VTARLV0oA+PaFCe7/30/zq189ky/E/M7FSX7+kRN8//pMSfc40O3CUaPfoRSqMkZCiE8IIc4LIc4KIb4rhOgsOPb6+dcvCSGeXuH6R4QQA/PnnRVCHKxmPQ00/uPqNFcnI5wcCuCxGjiwxUW/z0u/z0ub04x9DfH7VFayrYpak5lIislQih2ta98jns4VbWytlhvTUaxGPdFUjt3tlRUxVoPDol+UbayWr7xYek/hjakIH/vnS4TnlTXf98UT/MfVKT72rYs8da00QwRwYjDAdy5OrlsBZ7Vu2sNSyt8DEEL8N+BjwIeEEG7gL4D/JKUcFkK0rnKPj0op/7HKdTSYJ5zM8LmnXsqOaPGIl1wgu0lHtIQUvNlQ/veU1ahjT7uTsyMBclLbhfmarISTWTpcZvTz27KF+p3d7Q4sBl3FkiarkVMlu9odjAcTVQXGV0IAOgVyavFM4kQoRX+fl9noclG5SnjswiRnhgMc6lm9XUVKyb+cG2emIJj//O05nr89V9FzP/y3p3Ga9Tz34A+v+SVWLVXtjKSUhUEBGy/9v/wM8HUp5fD8eatPDWxQM4w6hRbHyunYVFYtaccyG0nTusp9irGvy8WpYc0Qgea6Dc7F8cfSXBoPc240xLnREB6rgb5mK1cnI0VVHmvFwGyUbk99Yh5HfR6yKjQ7TBzq0fr89ne58DVZ6fdpbTCZGv9uDz9+bc1zvvz8EH/2Hzdr+txwMrsueuFVx4yEEJ8UQowA70bbGQHsBDxCiKeEEKeEEO9Z5RafnHf1Pi2EWL+ihlcoJr2Cusq2OpOTJQ0xHA0m8NqMOEr8Nuz3eTlRRFq2GIF4hoFZbTd0aTyMr04iZeFEllCiulqjpfT7vOxoteflRmYiKc4Ma31+58dCDM7FOT7o51C3h0iZcbe1eOH2HI/PDzJYiauT1U//XYpJr+TrpOrJmu80IcQTQHuRQw9JKb8lpXwIeEgI8SDwK8DH5+97GPhhwAI8L4R4QUp5fck9HgQmASPweeC3gD9YYR0fBD4I0NPTU8KvduehqpKf++KLPHtz5S15X7MNp1mPy2IglFj9w7IgG9Lf5+XsSHBFeY9mu5GxYKKimpTd7Y6qxNRWY1e7o2hGr7/PSyiexmkxcHaV1LteEbitBra12EnnVAZmY4vqs1ZDIheNd+pwmTEbdBX37IG20/yv//cUr93eTE+TFa/ViMOs56eP9eA0G0hlc/zb+YmK778Sfc22mu/yirGmMZJS3l/ivR4FHkMzRqPAnJQyBsSEEN8HDgCLjJGUcuEvlxJC/B/gI6us4/NoBosjR45snlr6DSSWyvJX37/NvduaeNXWJv748WurGqJuj4WhueU1M2txfMDPrjYHLQ4T6ZzK5fHwIv1qu0m/rAapVG7PRNEromxlgFJY0Po+3OtBAIoQhAqKNkH7oFkMOuxmPclMjmA8TSyVw2zQMRZM0Gw3Fh0msBajgQRdbgtHfB7GAon8AIQmm5HtrfaK7rnAD27OQoEnllUlH37Ddr5/fZZwlVXnSxECvvC+o7it9VdNqCoiJYTYIaVcGGr1AHB1/t/fAv5cCKFH2/UcAz5d5PoOKeWE0OrffxS4WM167jSuTob5s+/doNNl5lVbmzg+sHqQMhhPVyw5e20qwrWpCIe63WRyOba12HCaDeSkXBQsLZcuj7Wq3cJK2Iw6rEYdh7rdnB8Jrqh5vdqzFQFOS2Wp7YlQsugElrlYms50DodZX3W7ygIv3J7jp4528yclxJTKRUp48fbcugwgqDY8/ikhxC5ABYaADwFIKa8IIb4DnJ8/9jdSyosAQojHgPdLKceBR4UQLWjJibML1zcojYVZX8/e0mIJS3vQlhJJ5Tjc6+bsSKjicdWKgMO93vxkjWppshkrMkbtLjMWgw5/LM2uNjt5gWshUFWJ1ajj+zdml12nBZZL+92P+Jb3pdWCC2Mh+vu83JyO4l+jsLQUnrkxy+H/+UQNVlack0OBzW+MpJTvXOXYw8DDRV5/W8G/31jN8+9kpJQ8dkHzcv/1/HjJlb6nhoLsaLMzHU6tGTMqZCF2cmoogCKgw2ViIlT5jmiB1YLtK7FUYqRU7eo2hwmjXqHDZSaraoF8VZWEEpm8G6JThFb0mcgQqIGhWInjA35cFkPN/o71otVh4jfu37kuz2q0g7wMyc6PfF6YP1/u53k2kiq7kK3dac4bAFVCNJmj1WFii8fC+dEg5RZQO8x6epusa+7mFp69ILTW3+dleIU5bKvhsRpoc5k5PxpipGiRZf0GAaxEKJHBYjSzf4uL0UCiJrukWtPtteK0rI+ZaLSDvAzRKYL3/9BWfvJIZWNysjm1rEBnp9u8LOMVSWWZjqQ4PRxkb0f5jZS72x1cHAsvW0erw4S3QNDrcI+HcCLN7nYHviYrxwf8ZY8/Aq2cwGxYv7aQUpkMJTk/GkKVkkM10E+qNaeGApwbWR9t8sbO6GWIEAKn2cA7D2/h70+OlH19JKUFUBdqR1b7Rt7X6USnCMZXUXm8OhXhrk4ndpOe2zOxkjrkRYGE2t4OJwadwGzQMRlKMhaMc7jHTTKrcn06Qjyj1iT9nytzLPd6EoxnGA8mONzrYXgujiJeUlPYaJrWSX+8YYxexrQ7zXz4Ddv47JPFxbFWYyGTc6zPy4sDfoSAviYbzQ4TUkoUIRgPJbhYQud9Oqvm63naneY1JTv0imBwLsa2FhuqhMsTy59xari4rlIxWh0mOlxmBudiqxY5nh4OsKPVjsdq5MpkuGbZrFqxIEsCWmvN3V0uDDrB5fEwyQ0Y4QRafG69OvcbxuhlTE+TlY+8eRfBeCY/hqgcdEILhO/vcjE1LzdRbfZoMpzEatRx1OfhxnSUYHxxkNxrM7DFYyWcyHBrpnicxm010OYwLxJF6/Va820u4WQGj9VIKJFBr9O0hBbUHrd4DATj6aL9dxKtgRa0HZ8iBKPBOP5YbSula0E8nePCmOYe7e1w4o+lKnJPq+XqRJgXB/y8bmdL3Z/VMEYvc4QQFTea5uTiTFQpbSKlEE/nODEYoN1pptkmcNs0vWy7Sc+poQD+WPEYhEERHOr1cG4kyLWpCPs6nViMOs6Phhjyxxla4/dcUHo84vNwco0M28KOz6AIfE3Wios214PLE2HMBoX+Pi/JTK7m2k+rEUvn+M7FyYYxarA237syxTNF6mkqwahX6PZasBn1qFIW1UMqh8lwkt4mK0LAlYmVYz47Wu14bEZC8cXV0aW4iMW4OBbCazXgj6+948mokia7aVMbI9Akfhf+Nsf6vEyGkqRzKorQ+g2n6xhfenFgjn87P8Hb93fU7RnQMEYva0YDcX7978/W7H6Dc3H6fd58/9XBbjcToQRSUvGbfWguzj1rZInmYqm8+1QLkhmVZruJdFYtSeVx1B9nX6eTSxPhqpUZ14NirST9fV5UVXJtKlLzWNjtmRh/8t1rhBIZfvRQ56Keu1rSSO2/jPmXcxM1f+MVNkSeHQkyFU5hM+nzfV6VYFvjzetrKk9OtRRGAwl2d6wsqrav00m/z0uP14JRrzASSHCo202Ptz4KAvXm+ICfk0MBzAYdh9fQPFoJIeDt+zt4/a4WnObF/2cDszF+95sXVmyWrgUNY/Qy5r339i6qyakFxcbvDMzG8FgNHPV51tQ40iuCw72e/Hhrp0XPi4N+XtXnXfGaZImazOWy0m7OazUwHkpyfNDPsD9BNJUllMhwejhIIp3j7i5XSfK8m5GZSIrTIwGO+so3SFubbXz2Z+7hkZ/v58Nv2L7suCrhcoWucyk03LSXMcmMmh9vIyXYzfp8lbROERX1n6WyOTrd5mV1RZNhLZtz1OfJf8j7+7yEExmMeoXzoyHcVgM7Wx3LZDb2dji5OBbUZrUhiaVzBOMZOlxmhIBEOsc9PW5SWZWcKrk6GaHDZS7aaFoOE6EEHqshP3ljgb75tpYFCo/PRFPMRFPYjDru6nJyZnjt6R6bDSnLr8oH7f9zgbsLBmPuanOwf4uLN+1t49jWplossSiN8dYvY/7wsSv81fdv53/WCTjY4yGnqkwEkwQSmbK21V6bEa/NuGqLhiKgt8mKXlGYjiTzdT3NduN8T9dyt9Fi0LG91caFsdK+VRfGHvX7PCX3na3EFreF0SVi9gZFcE+vpyQZj7s6nZgNOi6OhWo+NKCedLjMKIJVR5Ivpdlu4oUH34hepyCl5LuXp2i2Gzncu/KuthRKHW/d2Bm9TDkx6OeLzw4sei0nWfSNv6vNUdIAww6XGZNeIZ7OrdkrpkpNJ0gRyqICw9noylXciUyOYBlNuQup69PDAe7pcRNKZJiJpNZsYdnisZDJqosql7s8y41RZr45thQdpYVizsO9bk4NlV6IudFMhJJ0eywogpJlY+ZiKZJZFbtOQQjBW+4qpqlYPxrG6GVEKpvj1GCA716e4vvXZ9Z0H2aiKZxmff5D7LUZ2NHqWOS+SQmnhsvbfVydjHKkzJhENJllW4ttxULHYmRVbfzRAttabDTZTAgBsXRWy+pISTonUYTk4ngEo06bl+a2GAgmMkyGljfFbnFbGJiNlSXodmpIG3U9HkzUZZpJPYiksmXpV33wvq11F91fjYYx2iB+cGOWWzNR3nuvr+RrPvzoGZ64MlXy+f5Ymh6vhf1bXKRzKuOBZFUKg4UoZY6ADcQzBBMZ7up0ljzccSm3ZmKrGrO7u1xcGAstqo/qdJvZ2WZnKpwimsrS6TbjsRhRFK07bq1CykIW6nyO+jycqNJ9XA/KjcBsdMCmYYw2gBduz/Fb/3SesWACm0lPLJXFYdbzpr1tKw7N+8IPBooaIotBh04RbG+1c340uOyb0GbS84NVpGgrYXe7vaK2kW6PNT9QsB4Y9csN5NJA/Ig/wQjazqav2Uqz3biqi1mMRDpXcYJgvSnVTetyW9ZNt2glGsZonRnxx/nwo6fzgeWPfO1c/pjFoOONe1q5d1sTB7a42dOhdcw/e3OWsyMB7u5ycmEsTKfLTLvLTCqbYzaSxh9Lc3YkyIFuF5fHw4vct3oUqM1G05gNCslMeQHddpe5LsqJC5wdDtLttZQ0Rx5gYDZOX7OVbq+VWzPRosH3YlwcD+Oy6On2WCuuEl8PQokM/X2lqVVua7VjMW6sxErDGK0zDrOer7z/GE12I7/0ldOLAs6JjDbdYWHCg8Os51M/fjcPP36Nwbk4LXYTO1rt3JiOMl4k7X1uJMTudjtjwSRtTjPZnIpZr+A064mnczUTvZ+Nptnf5eL82Pr1SK1GX7ONZruRyVCSDmfpxgiYH5kUx6RX2N6qxZoGZmPMrSF0FkpkSWYiHOvzMhpIEElmai6GXwvGAgkcJh2RNQZ3xlIbv/aGMVpn3FZjXuL0d962m3d+7vkVz1VVyee/fzvfN7VQA7MaVyejtDpM+azY4Fwcj9WArsYTOMwVfIuqqsRm1BEroUWjVAyKoMVuytc2Lc2clUoqq+b/ZnpFcHeXi2Q2x41V+vNSWZmPwQlgZ5udeDqHSa/M10ypdLq13rxzq4xEqidjwURJSo1rteysBw1jtIEc6nbzsf+8l2+eHUOvCPQ6hZwqmY4kGQsk2NPhrGg089LKY5tJj16nVDXFYymV1KedHApwrM/LmeEA6Rp8MBUBTQ4T58deyrjVomwuq8q8fMe2FhuRZDb/N93WorWuLATSTXpBKiuRULSxeEHf2msz0u40cXmVhuGV2NXmwG7WcXM6VpZuOWg1YUNrNAE7zXpeVcdixlJpGKMNRFEUfuG1ffhjaf78ycUjiS0GpWaSHqOBBHs7nDU1RicGAyXHIwp5ccBfVUatkEPd7rJE2Crh1kwMh0nPoW430VSWWCrLeCjJFo+FFrsJl0XPdCRNJqeu2uzrj6XJqWrezV6JbqWo2cMAACAASURBVI+FTrclb+yFEFwaCxFN59jittDhWi4BvBpjgQReqxF/fLnb+a4jW3jT3na8NiOHeyvrZ6slDWO0CTja54UnF7+WyKi8OOCv6AO/FIMi6pLFOj7gZ0erHYdZz7nR0scfTYaSa6pBloIos7ygUiKpLGdGghh0gnaXmSO9HnSK4MJokHhBEP9wr5tYKreisQglsoQSUbw2I71ei+YazsTo8Vrx2gwgBedGAysMDNBc0KYyexGzqmRrq43MRG5R3KjLbeF33rZnXYYzlkrDGG0CVms+vTEV4ZjPA0KQyqpcGg+VFXsw6gQ72oqPea4F05EkFqONbo+lZE2guViafp8HfyxFpd7avs7KXNhitDlNtDhMCDQ53kA8g0En6HRbMBt0qKokk9P+9iP+xIoB8lNDQXq9ljVHdvtj6UW646VMSFlgLpbmWJ+Xk0OBko3/ycEAO9vsxKaj+TT/RCjB2ZEgr9vZwvO35+h0WfA11149oRwaxmgTsFrVayCeIaNKTs9XSd/T40anCPyxNMNzcYSAA91uroyHi2r3HOzx1DWd7rYaK1IePD4Y4FV9XtI5lVRWa/gNxNMMzsZKMlDlFl2CVoHua7KRzKhkcioeq4F4OsfgbCyvPV1IufVHAEPzhmpflxOzXseZkWDN65FeHPCXXXh5fSrKUZ+HKxNhoqkcqoRf+spp2pyasJwQ8C+/8lr2dZU/6aVWNIzRJkBZIzZU6M4stEdsb7XT5DCRUyUnBgP4mqxEi+xMzg4HsBh0JKqU6VipeK7ZblwzQLoSoWRmmQKk1aijx2tddWfha7KWXVZgMehodZgXtZfUk4vzTcFeq4HeJhtXJyNV/x8UcmIwwP4uFyOB+DJVgtWuOdjt5uyI9jdIZHL53ey7DndzV+fK+k/rQcMYbQIya3SDj/gTtDtNiwTZl27tF6REbs1EF7lx6ZzkyJbqXJrtLTZGgwk6XBbcVgNITYRtOpKqqoeg2I4wns7lJUS6PRYyOZVIKstkMEk0naPZbqTDZQagzWlmYDZWVLfIatThsRlxmPTE0llMOqUm447KxR/P4I8HsRi0nrlqpXwLOT8WYm+Hk2RWJVFiucSFsRDdHsuiuJTFoOMDP7R13WJwK9EwRpsAVUr2dDi5PRNdUaaizWledTpEs93EMzdnaXeaMelF3l0AGJyLldW9vRSzQUcyozJQRHitmCHo9lqYi6bZ1+Xi9JAfl8VYtIgwu4o/NhFKLtIzclkMbG+10OU2MzgbYyKUzH+r7+1wYDXqEUIL2BoUhcsTIcY2UUNrIqNyfSrKoW4354q07VTK5YkwHS5zycYop0pNbzyZyVec/9lPH2J7q702C6qChjHaBGxtsfPtX7uPZCbH45cm+dxTt7g+FVn0hjXpFxcZ7m53EEtlGQkk2N3uYCSgfTAnw0nMBoU2h4mpSAqjTtDuNLOz1cFztyvrUbOZtFHUqaxKJJFZtWixxW7KB3hPDvrxNdtQYJExanWY2NZiI7ZGVXAhoUSGUCLDzekoRr2CXqfQ6jTSZDeu28TTWnBmJFiTDGkhQmjzzRY0l9qc2thxgWA2mlqWWDg/GmLrfLA6nMhWJSlcSxrGaBNhNuh44GAXDxzsIhBL887PPZeXgVUL/CGDTpDM5HBZDcTSOabCyUVxg2RGxWXR4kpDc7F8/5TDrK9IM/viWAiPTWu3AOhym+nyWJFSIhBkcioGnYIqJZFkNh/jUiUEYml2tDk4ZjchkcRSOQKxNM/f9mPQiYoaVdNZlTQQTycYDybY1ebAZFCIp7PcnC5domTDqHEh9ngwmW8I9lgNzEZeGgZ5pNdTNMt5ezaGxaijw2XmH0+NksqqGPUK50aCzEZT/P477lozlllrGsZok2Iz6ZkKv+SmDM3F0Suaxo9Jr2WepsJJsqosmuovnE66QI/XyqXxMAadYFuLndFAvOiww94mbWCiIrShjnaznomCNouxYHJFBcHXbm9eJOgWWDJ+qJBMTtLmNNPttXKmwsCyKsk/z6hX2NPhQCCKTqndNNTxM15qMBs09YG+JitjgTh/+NgVmh0mMjkVt6W4ckS92Rz7swbLMOoV/v+fOZT/2W7S85a7tLlVO1rthBJZEhm1rJqjeDqHbr4U4OpkhJ1tDnxNL03D0CmCLo8Fp1nPycEAp4f87O1w0uu1lVwPlC6zuPLSeJi5ClLoxZ+tcmUiwkQowY5WO0d6PWxZp9HM5XB8wE/nfBC+3pwdCXB4vkizGJcnIpwaDnJ7NsbxAT9nhoP8tx/ese67IqhyZySE+ATwAKAC08D7pJTjQoiPAu8ueMYeoEVK6V9yfR/wd0ATcAr4OSllbd6ZrwBeu72F+/e08oH7tnJsaxPJTA5fk5UTBYL3LouBRCZXktb1wGwMr83IwHxf1enhIIrQDN3eTifnRgKMBRL5wG9WpewdxunhAC6LoaweqmF/vGwVyNUIxDP5HYJRr7C3w4nNpOPqZO1nipXK9lY7k6EEOkUQTWbLkuGthqyqSRF3eywgtMysXhHc1enkwlhoUVzSYdazxWPlUIWjjqqlKkF+IYRTShme//d/A/ZKKT+05Jx3AL8hpXxjkev/Afi6lPLvhBB/CZyTUn5urefeSYL8wXh6Wcl+JJnh/zw7yM42O2/e286jx4f5vW9erPgZnW4z2RpNJa1GWuRYn5fTQwEydRItc1kM2Ey6ZYJr9cRu1NHltXJtMoLFqNNco2YroUR2URX2evDrP7yDdx3txmM1YjHq+PrpUb51dhyAt+5r56f6e+ry3FIF+Ws2HUQI8SDQI6X8pSWv/y3wpJTyr5e8LoAZoF1KmRVCvBr4fSnlW9Z61p1kjNZiJpLiwa9fKEuOtpC7Op1MBBMljYIuhUoC0oWYDQp3d7kYDSSqHlVUDJNeYVe7A4NOi1CksjlGAwmCNfr9Czni8xSVDjnic3NtIkpknTWEdIqg12vlnYe38Auv6cNsUNaltmjdjJEQ4pPAe4AQ8AYp5UzBMSswCmwv4qI1Ay9IKbfP/9wNfFtKuW+F53wQ+CBAT0/P4aGhoarW/UpASsl9f/xkxQLxRr3CnnYHI/54zYyRw6SvyYfMZtThMBtwWw04LQYGZmJVN9auht2kx2xQiCS1VHdvkw2LQUcmlyOazs2P+E7SbDPR7DChqnJ+0KNACC1edX40mI+tFY4JX8CoE2xtseOy6AklssxF04t+J1+TFYtRx7XJSM3qkFbj1Vub+Jv3HsFWZxH+mhkjIcQTQLGZJQ9JKb9VcN6DgFlK+fGC134S+Fkp5TuK3LcsY1RIY2ekcW4kyAOffbYm9+rv8+KPpapKjXuthpoZtWL0NVsRiKJTbzcDXpvW+hFKZLi9JP7V4TLjshgWVYH3NVtpsmuG7cpkJF+42Oky47UZ10XS9hde08fH3rG3rs+o2dw0KeX9JT7zUeAx4OMFr/0U8NUVzp8D3EIIvZQyC2wBxkp8VgNgb6eTLreFsQrVDQs5PuDHYa7uG7Lba8Ufr18BoiYRqwWDW+1Gzo+FipYmbBT+WAZ/TCtR8FoNbPFamYmkyOTUZRXloP0+A7NxjvR6FlVQj4eSzERT9Hotiyrp68EXnx3gvp3NvGFXa12fUwrVZtN2SClvzP/4AHC14JgLeB3ws8WulVJKIcSTwE+gZdTeC3yr2LkNipNTJZFk7XYiviZbXuGwItYpG3xzOorX6gGpuYUdbi1N7pyfrBJNZTHoBOFEFr1OMBNNsaPFgU4nUKUkk1XR6xTODAfQKaIqOVi9Ah6rCbNRod1pJp1VCSYyTIdTJasZFNaTLZDJSWwmA15rtq67TYC/evrWy98YAZ8SQuxCS+0PAYWZtB8DviulXLRfFUI8BrxfSjkO/Bbwd0KI/wmcAb5Q5XpeMYwFEzzy7AC/8oYduKzFi9BGA/GaicC7rQZulDB9djVSGbWqHriyECIvmRJZo/lUr4iigyr3b3ExGUrS47WSzqrkpFantHT9Rp2gp8mKw2Qgp0qMegVFkO/XW4j7lDMIoBCTobie+OWJMPs6nXU3RhvRQFyMqoyRlPKdqxx7BHikyOtvK/j3baC/mjW8UpgMJWl3mRnxx3nmxix/9r0bTIaT9DbZ+NlX9ebPCyczKEJwcSxUs11Rs92Iy2LgVry6WMzVyUjdBxw6zHp8TdayertWGkSwsHMpLGmwGnV0uS3YTDpMeh3xdI6B2Vhd20z80RTH+ryoUnJ9KrqoRkunCDxWw6LKaq/VQLfXSqyEceSl4DDrefbmLPdua9rQzv1GO8gGc2LQz6999QyhRAa7Wb+sheMHN2Y54vOwu92JlJK3/ukz+RjRge7qhLAcZj072xxMh5M1KzisV+2M22pga7ONiVCSC2P1C+zG07lVNarrQSKr5qeMaDVILxmec6MhDDqB12Zke6udS2MhtrfZOT6gGfxuj4U2pxlFCIbmYkxVUCs24k/w7r95ka996NUc9Xlr80tVQKMdZIP5+xMjTIST8w2vy99I37k0yV89fZt4Oks6p+It0EA266sbutdsN3FqaGXN5XIx65UVJ+JWisdqoN/nJZtTOT0crEvt0Uazv0BdMZrM0eO1LjqeyUn8sTTHB/y4rIa8IQIYCSQ4ORTg+KAfs1G3qoTxWjx1bRp1A6fkNnZGG8yf/JcDRJNZvnNpcsVzvnFmjEQ6x/XpyKKUcbVyprWaPrKA22ogmclxuNezaDhlNaSzKicG/Rs+B74Qg666oHch21ps+V0RaKqeXW7LMtdsgdWqx4fm4nS6XxoYcHk8lBfh393uwGkxkMrmmAwl8ViNhBLp/CglgM8+eYtURuXX37RzVSnketEwRhtMTpW89e72VY0RgD+eWla7cnUygtOiL3ks81Jq3Z09GU7lBeCO9HpQpaxa5nV7q51zFWhs14P2efXJq5MRPFYdJoNCs81IKiu5NhWp6MvBWeT/YCyY4FCPm0AFf7ulciI2ozZFZWmQeiqcwqATecMl0ZKhL9ye40///Tq/+5/rW3tUjIYx2mA++o/n+PrptcurZiPLYzHRVJZ+n4fjFQaMB2ZjdLhMi74da8HWZhtCsEzfuhIW2jY2ErNBYU+Hk5vT0XxwfqHhdiGD1mQz0uY0YTPpuTgeztcNuSwGvDYDORVyUuIw6ecbidPz6pTLd6e1El9bS04kk5OL5IgtBh1HfB5+5lh9etTW4o42RqoqN0QqoZDxEgoWtzbb8kqOS7k2FcVqUBbN7yqVaCrL3c2umhojgeZaZXMSp1lfshzqStTalawEvaJwcTS0agPvXCydV7Pc1eZArxOaYRoNMjBbetaz1iqQpfLAwU7++Cf2L1MUXU/uSGOkqpLxUIL3fPE4D71tDz+8p23d1yCl5DPfu87N6Sj983PRQPswL0wQXWC1wrxQIkOX20KPWV9WvYjDrGd3u6PmaXiJNmxwYeZ9m8NELJ1lW0v57tbdXS5uzaxvZqsYW1tsZY1julZFvVaozjVFK/Hee30baojgDjVGf/HUTT79xA1yquTvT4xwpNe7YmFhvTg+4OdPn9BGWi/tcm+2G/E12/K9SWtFIsaCmkZNv8/LhbEgiTV2Sb4mK71NNgJFRh7XmoVUcyVzzox6pSoFANACxB6rEUXAbCy9LO62GrvbHYCmt7ReOC0b85Ec8ce5Z4N0jBa444zRk9em+cunb+eDjd+9PMW7R7XJmuvFP58b57f/6fyKx2ej6fx0T1+TtaTCtqwqOT7oRz+v1thiN2LU6xieiy2aKnKw2821yTCDc3F04iVd7IVallAiQ3LemPX7vCCo2m3wWA1cHA8hBJQjEqEvcNEsBmVVI7u91Y7HakAgFnXLx9M5piOaqJoioN/n4fp0dEXJkG6vhVgqy/ZWx4a4SzVS9Cmbv37mNg8c7NqYh89zxxijuWiKR18c5q+fuU10icSF22LQxOXrXH2azOT41Lev8shzg2ueOxNJ4bHo0evKW1NWlYvUGg90uxYZo9loKv+hzknY4rHgMBuYiaTyo4iMOoHHamRoLoak+lR2Jqti1ivs7nCW5RYm0lkOdLsw6XRcmQjT5tTaMXzNNoKxtOYSBuLc3eXmykQ4b7TbXWZSmRzbWuzYTXpeHNCmoqhSm2Tb5Tazs80xP1AAbXR4JseIP54PSG+EIQLNfd8IbkxFuT6lSRFvFDUTV1tPKpEQeeTZAX7/Xy6veHx3u4OH3r6HVoeZuViKy+NhfE02HGY9c7E0/+mu9qqC3bdmovzmP5zj3Ehp6Voh4N6tTTx7q7LxQoXs7XBi0Aumwil6vFauTYbZ1eYkk8txpoQxPwZF0OI0Va2QeGCLa8240f4uF1lVYjPpEIJFBX6V0uow0e4yVzSGe70ppoO0Xhzu9fBPv3Rvze9bMwmRVwprxR6uTkb4uS8cX/H47759D+95tQ9jmTOmMjmVLz03yKf//fqq88YW6PFaabIZuDIRqYkhgsU61gvjhhbe8Ae7XVwaC6+aKcqoEikp283SCTji03qukpkcAjjq85Cc1+x2W41MhZOMBBL0eC0kMioXxkM1cVVMeoVUVs27rbc3QSC8FDK58rOi1WA36dneauf+Pa38dJ1kZ0vljjBGyUyOL5XgGq3G//y3K+xud/LaHc0lPe/8aIgnrkzxr+fGGS+jhSGezjITSZIsQWC/FpwdCdHbZCWSyOJfJaDd47HQbDdiMegQQiuSQ0IsnWU8mFhW09LltuAw6xdVF6/GglZRrTjY7SaUSKNTBGeGgwi03e9m6VAvhq/Jys11NJoui4Hf/5G9/NihLev2zNW4I4zRjana6A1/5GvnuH9vK71eG80OI4rQ9HFSGZXxUJLL4yEujIWYjaYrbtWYjaa5q9NJJJmlzWmqawf8AkNzmsCXf2i5MdreaqPJZmIilFwxqySElgHs8Vrzk0qsRn112khVsnS9Em33W8spJLWmxWEqOnCxXoQSGVwbNCOtGHeEMfrOpYma3GcynOQrLwzX5F6rcWk+pW836ej2WirWySmH+LwLqcy7VjlVMhFMcHM6xk1W//BKqRnR2Whau7638qrwWnBXp4NL48V3QPF0jnaXOe+ubibODge5p8dddQtNKbQ4TPzy67fxxt3rX2O3EneEMaq34Hg92Npi43IN2ilK5fJEGF+zlWa7qapM0kLGqpbNsuWw1od5IpTEYdavnwhcGWTq3BFg0In87/yF9x5h/xZ33Z5VCS+/T2mZpLNqzSaWricbMWLY57Xx1PWZtU8sgVNDAfZvca1rBsukV0qalBJJZkvK7G0Eso6hwr0dTv7xl+5lPJigt8lWvwdVyMZ3IdaZqXASm3Fjy9wrIbYBQvOJdBZ3DSvRJ4LJdY1JHNjiLnkQ5bnRELvb7Vg22XtjOlo/91GrGVM2pSGCO8AYdXut65aZqhVHej1V9TetRaFWjWVef3l7q52pSIptzfaaPWcmmsJh1tPmrFzwq1S8NkPZ9TlXJ6N1NZYLDpfXasRQQvHqgS2uusYHQ4kMsXUeHFkOr3g37dJ4iEeeHdzoZZTFVKS+wdW9HU7GgnGa7CZuz0Q56vMQTGQYnIsTTWVrGk8ZDSRoshnZ4rFUPGyyFIx6Hcf6HKhSMh1Oljzip9tjqSqYbdAJ2uangjTZjVgNeq0eC+29Z1AU/PE0LouBXW0OJC9V+sdSGewmA/6YVoJQzy8g0LKmT1yZ2vC2j5V4RRujUDzDR752nnQdCslKqVnZ1WbnxnQURUCpmzODThCog460Ua/gthiYjqSQSMaCScbmK6oLywdmo+map7/nYml2O+rbZjAZSuaNSjnSq5fGQvT3ebk2GVmkPW3UCfZ0OEllVYw6BZNBcyKC8Qw6ReA067k5EyOUyOSNbDEXMYn2Hx9KZDassrqQE4P+TWuMXtFumtOi50Ov21rz+x7r82I36zHpxSJN6jaHCZNe+9brdJuJzaeRdYqC12pc6XaLuLvLVZPBhFtbFscF7u5yEU9lOdbnXbNbvx5SErdmoqzX4InZaIpSk1LxjMrxAT+dbjNtThOKgLs6nTgtBs6Nhrg6GeH8WIiZSIqcKvHYjDjNBi6Nh/HHKq8n2yieuTHLdy6uriq6Ubyid0YA/3x2HJ0iavam2d3u4NRQAKdFT6vTzIhfk+/obbIyHU7R6jRjN+mJprKL/H+PTeGoT5NoEEIgpUQR2roUIYilM7Q4zFyuolDQatSxo82OQadwYTSYnwsmIJ9mL6UiWgjJ9lY7t2eiNXHXdrXZcVmNnB0Jkl6H+J0qS+uDK+TKRASbQWH/FheqCp1uC74mG4oimAolGJyLLypIPNbnZTaa4vZsbMM67SthaC7Ol54b5C13tW3oWKJivKIbZWejKZxmA8/enOXRF4d44sp0xc/c1qI1zU6FU3WdUFFpyrnHa8VjNdQ0Xb29xYbdrMegU8ipkploquwA66EeN2fWoYhvKcf6vCW3olTD3V1OgvEM46Hky2aX9I1fvpdD66hd1GiURRvFA/CG3a28bmcLlyfC/NwXXlxTG3gBm1GHXifY2mwnlMxwtoQO92pZaeDgaggBRr2oed3MzZkYh3tfGsq4t8PBvk4nU+EUwUR6TVmRI72eRRrL60mkRpN212JhhtuCcqYQgutTkRX1kjYDSyV0NguvaGNUiKII9nW5ePgnDvD+L6+9q9IJMBkU/LEMZ0qU/agFNlP58RqzQcdQjRtNFyiMvRRWhFsMOg50O/FH07Q4TMRSWUx6TfZDEZDOSW7PblynfDSVxWsz1m2o5FIiyWzeaBt0gr0dWmZPpyj59p7Nwvu/dJLPv+cIR32evHrCZuCOMUYL3L+3jR892Mk3z46vet5h38YIo5c7IMxlMdBkM3J7tj7Nn/EVZE8SmRwn5z989Xp2NQz743S6zAjK/pNWTSYnFxnu3e0OoqlsXUsbyiGVVXnvFzW5nCabkVO/96YNXpHGKzqbthJN9tVTvx1OU13S66VwvczRyls8lroag/Vyd+pBp9uyKYY/Xp2MMBpIcKxv40ZHr8RcLM3IOmp8r8YdaYx+40078azQ9vD2/R08+dE38F9ft4333etb34UBHquRbS029nU56V0y5rgYlbh1pXJgi2tdxehrxT09bg71uAkmNlfc5sUBP7s2UNZ1JTZLqv+Oc9NAa4f4xi+/hl//+7NMhZPMRFJYjDpet7OF//GW3ZgNOn7i8JYN+SAOFOxyhIBOlxm9TqHDZSaWymI16cnmVBQhuDYZZqJKKdjVqNUI5/VGAJFUFoEWdBdCbJq4zXgosekkTL56YphjW70b3sX/ik7tr4VW4wPJjIpBJ1CEWCTh8MLtOb53ZYrnbs1tmjdzIa/ZVhuN7JW4q9O5KX/vcjjc6yGWym4qhUezXmFXuyOf/dQpAodZT4/XilmvQyJJZHJcHFv8t692MMJqGHSCL//CMV69ranm916X1L4Q4hPAA4AKTAPvk1KOCyE+Cry74Bl7gBYppX/J9Y8ArwMWctLvk1KerWZN5bAwrXSlzu1XbW3iVVubkFIyHUnxuadu8S/nxvOTQzeKhabLsRKm0VZDMrP+ygG1JhhPbzplx2RW5fxoiH6fF4nkxGCAYDxDML64NOOeHjeZnESVEn8szXQ4ya52B9fqYFgzOcnJQX9djFGpVBszelhKuV9KeRD4V+BjAFLKh6WUB+dffxB4eqkhKuCjC+eupyEqByG0ZsiPv2Pvho5yWeBQjwerUVd3iVKHybCo3eXlxt1dTmajtRvdXUsk2lCE1WSFU1mVC2MhLo2HmQglyUkwG+oX5tWVORar1lS1M5JSFu4jbRTPov408NVqnrNZuD4V3RD1wkIO93rWreTg7GgQndDGB1mMOpKZHKH57v7Nzv4tLi6Ph0puUN6MFIugnBsJ1c19bnWYa37Pcqg6gC2E+CTwHjRX6w1LjlmB/wT8yiq3+KQQ4mPA94DfllIW/SoTQnwQ+CBAT8/6j1SRUhuFXQ8FgHLIrPOnKyfh/JJ+uR1tdiaCyU1XyWvSK+xqc2A26jg56K+ZDEqP14LFoCOrSjxWY16oKJuTJDNZhBDMRtLM1HgXllWL/18blNrvjuwmPfeVMPmmnqxpjIQQTwDtRQ49JKX8lpTyIeAhIcSDaEbn4wXnvAN4dhUX7UFgEjACnwd+C/iDYidKKT8/fw5HjhxZ96j70FycvztRfzH+tbg1E13XyuKl7Gi10+YwbZralAWa7Ua8NuMyw1kNAm0nurilpXj8aUGvqJaaRJ4VKqOH/LGa9/y9YXcrbc5NvjOSUt5f4r0eBR5jsTH6KVZx0aSUC2M7UkKI/wN8pMRnrTtzsfSK1cjrSSyd465OF8dj61cdbjfp2NqiKUCa9QqXJyI4zQaSmY2Px1iNOvZ0OBmei3N9qryC0f1dLox6hWQmx6Xx8LIYw9EypruGEhni6Sz9fcUr9/VK6ZpWC5wdCeC1GfDHFtdLBeIZAvOTRM6MBKtWDdjRaufnX+Or7iY1oNps2g4p5Y35Hx8ArhYcc6Flyn52les7pJQTQtMy+FHgYjXrqSfJTA6LQWFfl4tEOsd0JEU6p657Q6RBJ9bdPTLolGXC+kd6PSXrTdeLHa12jHql5Dheq8OE12bEpFcIJTKLptc22Yxsa7VxYTSEyaBjZ5uDy2XGZTI5yfEBPzvb7DjNWlGtEFrz85WJMCa9oK/ZRiydw27SYTVqH79sTuVUkV1OKisRYuUvwNPDQY76tEZmh0lf8WzA/j4vu9s3PjFTbczoU0KIXWip/SHgQwXHfgz4rpRy0b5WCPEY8H4p5TjwqBCiBW1HfHbJ9ZuGa5MRvvLCEDvbHfnsh1GvaPUibXaulfmNXA37ulzrLskRiGvD/gqVENUSv44NOoFep7Cn3VHzeWA35ltndrbZGZqLkcoWX1Obw0Rvs43TQ4EVDehcLM3cQBohtNR7NUmC1XZoF4sYOItBx1Gfh3MjQdIFdUQCbeeXzKy8TsL3xgAAIABJREFUpbo9HWVvh5MrE+GKVRIefXGYA91u3nWku+xra8kdXfRYKp///i3+12NXix7TK4ID3e51ybJt1Hhmu0lzhZamoRe+lQsx6RV8zTZcZj3D/jiz0VTePTnS6yGVzWHQKRh0CtORJJFklq0tdk4PBSqST1mgy2PBazVwYSyM06LHZTHQYjdxdiS46eajrYTLYmBHm51cTqLXCUKJzIqGrb/PSzqroqpyUZysv8/LmeFA2cWRDxzs5DM/daiq9a9EQ8+oAs4MB5aJTs1GU/zl07dXvCarSs4MB7ir04nNpGcumqpLkd3ONvuK2ZVK2d5qZyKYILZGLGwlGdxb0zFaHSa63BaMeoVUNodRr6BK7ZpMTi6Kkyz91rYYFPSK4PiAvyrto/1bXOiEwKAT9Pd5uToRZsSfWJdJvLUklMjklRBWQydY0eAcH/CzvcWGhLLeh8/VsZK/VBrGqIDPPnmTn31VL6/d3swfP36N33nbHv7HP55fM3OlypdGUtuMOvp9XhCSuWgaKTUtpVvTUVxWA20OE0LAzeloSQFNvSK4p8dd8bhoh0nPp965n0vjIUYDCbw2I4l0jvt2NnNgi5sHv36BH9ycXfM+A7MxjHplkWysf15LeyXX51C3e9Vq9USB+3FyKMCuNgcOix4p4cpEuKSEgW7+b7uWQX0lIQG9opDJFf+db87E2NfpLOueH3/H3hqsrDoaxqiAHzu0hd/8h3N84oF9XB4Pc+gPvluyKuQCsXSuaAbGoBOEE5l8wNti1LG31Y7FoOPGVDT/wV6KKmVZMhh2k57339eH3aTJxR71ednb6eTt+zuKnv+5n72Ht37mmTW1dvqabZwtU2ROX2ZFb2FavNluRAH6WuxcnQzjshiIJLOkllhwVUoyG1z7td6oElqdJoZWKT6dDCfp93lW/RIz6hUcJj2/8sbt/Of9nfVYalk0YkZLmIumeN3DT61rxmqldDDAfTuaiSaz3Lu9idNDQZ6/PcfdXS7edbSbRDrL3V1unBY9VqOewdkYu9oddLotZT0/lc3x+KUpvn99hscvTrLFa+Vwr5uB2RgnBgIc7HGXHdDd2WYnns5VLSjW12xjNpoikswWjVEd7vFwanhjq+I3gn1dzmWNtMXY2WYvGnf6ySPdfOJH9wGaUaonjZhRhTTZTdy3o5lvr6PGy7mRIP0+L2PBxLLm13u3NfOLr+3DqNdE8Z+8Os32Vju+5uUjivuKvFYKJr2OHznQydvv7uCP3rk/30AMMOKPMzQX4wNfPkWijMbZmUiK3iZr1caoUFLlxGCA3e12bk7H8sHuudjG1zrVE4/VwI5WB0LAdCRJi8Ocl5AphWITc+/f08Yf/vjdixQqNgMNY1SE997r4+nrM+tW5JjKqhwf9OO2GpbJpH7uqZv8+D1dtDnN6BTB/Xvb6rYOXZE3Z7fXSrfXyut3taxpoN1WAw6znjfvbef5W3PcLFO1shSuTkbz3/YdLvOKfXJNNuOGqytUgsdqoNNtQa8ITHqFmWh6kds/UKbW+ZXx8KIsbKfLzG++aeemM0TQMEZF2dpso8luJL7O2ZhgPEN/QdWvThG8717fhpfpZ3Mqz9/Wsi0GnaDdpc2LW5i0+uqtTVyfivD/2jvv+LiqM+9/zxTNjKRR773ZcrdsyRK2ARtjSkh2TSghhE7oLZtNIXnZD5td4t1NeAMJL9lkKUsJkM0SUngdNgkQehI33I3BxpItWbJlWb1rZs7+MVeyNJqRps+1dL6fz3x0dXXPncfHd5455Xl+zz9tWMSauZkA/PfWRr75yu6I2OOSkspsO0k2Ey1dg6QnxFGRlUhDWx9DThcpNjPHuwZZmJfEgeM9Z0wJoSy7hfa+4bAmwfYOOxHCXTLKpDmgYEfQkUY5Iw8a2/tZ/8g7kxZKo8V4DaFz52Tw9xdWxsSO8Thckte/uoZtDe2cPScDp0vy/F+O8MXaQlwuyEm20tk/PKHKRCR1cQ61uqduZqMgN9lCTpJ1Qo200U2Cfc3dZNktFKXFj4UNFKa5Rx2BjjCiQVF6fESi2j8aVxzAYjJErbJvoChn5IWitPix6N5oYjIIEHDBgmwe2rAIS4QXFv3FajZiNRv5zOLTO3L3nT9nwjWe5W5yk60RD9IccUpauoboHXJSq4ndey60t/YM0dozRGlGAkaD4FBr75gsijvf0EGW3conrT1RrQxbXZxKZ/8wqfFx9A87iI8zsc+PBelQefBvFmA1R043PRSUM/KgMC2eP/zduVzw6DtTBo3ZLSb6hh0sLkjGYjTS1jfE4RCDHR0uSXlmIo9eVRXSffSAyWjgrLL0qESM9ww62FLfjtVkoKowhX3NXZMCAscvhI/KotjMBhxOSUd/D7nJVvJTbFErOtnY3q+NgqKrQmm3ei9EoQeUM/JCY0e/V0eUEGdkYX4ygyNOTAaBEO6Q/V2tXRgEzM+1TxgSB8OJ7kGklGN10D841Eb3wAgWs4F18yK3eB0JLl6Uw7N/boja+w06XOxs7GRBbhIWk4H2/uEpY3HGB122dA3S0jXIwjw7bb3DnOgeGisZbjYaMBgEA8NO9oRBomR5UUrY8/T85Rsv76J7YIRrzyqOyftPhXJGXvjjvhNez8/NtrOlvt2rHIRLuufmFZkJpCbE0dE/4nU3KSHOSGWOHbPRQP+wk4ERJ2kJcXRpQY9//vQUr+5qpjg9gbcOtPKjNw+Otf3cklw2XrqYZB9llvSElJIfv3UoJu+9v+X0dGd+rp04k4H9zd1+5Wvta+7BYnKvRR1t7+eoR3hVOEp2Byp1Ek6GHC7+edN+rqkrGvvC0wvKGXlhR+Pkh81uMbFX+1acam370Mk+0EZVSwuSOdLeP0FmpCIrcdpvxa/8l3cp8E27W/jMolyf0dR6QgjByRhLjMDpxdviNBsmo4FkmxmjQbDnWJfPbPghh3styhvNXQPMyUoMek2xpjj2QZrDDhe7mrqoKoxtaSJP9LFCqiNGnC4OevnmKs1IYCTALeJdTV0IoKowhdrSNFaUpPpM+/CX1/a2TH+RTvjsYv04zSPtA3x6so8Pj3aytaEj6O3t5s5BDrb2sqIkNagNhm1HOijyozhnpLnpmS089d5hXYU9KGfkwfd/f2DSt15NcWrQcqYd/SPsbOxkS727EkSomeQXL/SmAKxP7j1/Dpctz4+1GV6pP9k3tgMXDFsbOoKqFFOWkTBBFypWdPSP8N3ffTRhGSDWqGnaOD441MaT79VPOn+sM/YxKTazke9eukhXow1/+M7fLgTgVx8ei7ElExkVUAtFT3xfcxfpCXGkJphJtJgwGtwKkr4iz5cUJHOotVcX8sXgDqpdVqSfqZpyRhovbT7KQ5v2Tzq/tCB5rPJnLEm2makuTtVlGP9UJFnNbLx0Mfubu3VV1XWU7hBGKS6pKUR6OLOcJCvHuyeWr64pTsUlZUD5fZHmmroizqvMirUZY6hpGu6dn5e3N3p9UIYc+nh4jncP8ugbn8TajKCwxRnZdO/ZupuyxRlFSOqSvihIm6iaUFuaxrYjHXx4NHTx/HDy/F+O8Ojr+nmmZr0zOniih0df/8SrrnRCnJGeQX04I4DdTV26WG8IBpPRwH3r5pAfoLxJJDEZBBmJ4a+Y23iqn6UFyeQkWylIsbHNzwojseCX25to6oj9MgTMcmc0OOLk1ue38difvMfDVObYI17PPhDuWlvuVRLiTKEkI4F71lWgl5lm/4gLp0tSW5oaVptO9Ayx91gXeSlWkmxmXWtwH+scYM3Db3Pvz3fE2pTZ7YysZiP/ctlir4mDpRnxMYuS9caFC7K5bHlBrM0Imatri7i6NvoVgX3R0T/ClvoOFuQmkZscPnWEBXlu8bPxAZh6xGwU5KfYWFkWucRmf5nVzgjc4mW3nF069nt8nJG60jT6fYjQB0pdaRoZiXEsLUie8rq52YnYzJP/O4SA8yoz+cr6OV71hs5E7l03hzvXlpObbB2TtYg1e5u76R1yF2HMslswByiZ64nNbJqgF65XHr5iKe9+8zy+VBf7LwglOwsMDDv57GPvEWcy0DfkoDFEdUJPVpSkcqxzAIFAImnudO+0ZNkt5CRb6R92urPJDUIr5+NixOmif8jJj66uYkmBfrZfw42UkvcPtfHgb/dNSGaNNRmJcZRlJPLxiZ6g1+lqS1PZUq9fSdyyjAReuXMVgw4nucmRW8vzV3ZWOSONDw61cctzWyckT4ZKaUYCiRbTWHKlQbhVHCuz7Yw4XQw5XJNkWZcWprB2biYGIfjskhwqsmJf6TMa9A45ePA3e/nVjmPYrSaQBF0hNZyYtbp4Hf3DAZegWpyfRHycaYLWkh5Jspr49d2rKddKmIcb5YyC4Lbnt/HH/d6TZANleVEKu5q6Ag63f+cbaylO16cSX6QZTcVZkJfEwLCTD492cORUP1WFKTz25kF+vy96uuSeGAQszndPtVu6BgMSQZuq4IJeMBoED35uATesKgn7vZUzCoIjp/q44JF3GQ6h9I0AakpSp0zE9EVxultLSa/iV7FESslzf27ge7//WBeBg/Ny7HQPjoxNuadjUV4SH50BErhPXl/DBWHWWffXGc36BezxFKcncPO4xexgqC1NY2tDR8COCODIqX7WPvw2nSEm085EhBDcuLqUzQ+cT10IOWXh4sDxHrIS/d9929vczaK8JMozE3S9EVHfFjt5E+WMPLhzbTkpIegFtfUOUVuSFnTOz/HuQd1/e8aSJKuZR6+qYl6OPeYf6p1NnSzItfsdo7SrqYvDJ/soTLWRlqC/eDGzUbChKnZR8soZeZBsM3PrOWVBt//0ZB9bGtr9rmvljVDazgbyUmy8dt85PHFddcyDQD850UNFViI1xanUlqRSU5xKXWkaSVbvaZ8S6B500N43Qk1xanSNnYYRp+QDP0qdRwqVKOuFG1eV8PT79QFlc9eWpCIBp8stGXvwRHBJod+8uDKkkdlswWAQnD8/m9e/ei7NXYNc9/Rmegajv/vmcHlXbqwuTuWTEz1ebRp9rvT4nVORFZkdNX9QIyMvJFhM3BjgrsLWhg7a+4b58Ggn24900B3kB8NuMelODlTPZCVZqSpM4c2vraFMR/XAth/pwGY2UlPiHil5E2LTS07YeK746V/45fYmWnsGo54kHrIzEkI8JITYLYTYKYT4oxAiTzufLIT4/0KIXUKIfUKIm3y0rxZC7BFCHBJCPCZ08km8fmUxtgB2tSTQNeAuwjgvJ/DYoDiTgSurC/hSnf6E0s8EsuxWrqjRV7pMa88Q2xo62FzfTomXcI2TPUPUlqbpaoQ07HDx9Zd3UbvxTW56ZmtU3zscI6OHpZRLpJRVwCbgQe383cB+KeVSYC3wAyGEtxTpnwC3AnO018VhsClkUuLjuGpFYUBt2rRSxHFG/7p1ZVk6ly3P55U7V7LzwQt4+MqlMV+UPZO5YWUJN68ObTc0ElRm2/n05OSpnMPlrvNWVZiiy6l5tCvPhuyMpJTjMwETOF0qXgJ2baSTCLQDE+YuQohcIElK+VfpDnh6Hrg0VJvCxT3rKoJq5/Qzdmv9gmwe+UIV1cVpxMep5btQSbCYuH1NGV9dP5f187N97nJVZtu5fU0ZL9+xkuURVjo0GwX9I44pdZN2HO2ks3+EhXlJYwvyoebGhYNf7zjGI1HUOwrLJ0AIsRG4HugCztNOPw68CjQDduAqKaVn8E0+0DTu9ybtnLf3uA24DaCoKDpJfRmJFupK0wIO5z/a3k9uspWWLt8BcQWpNt1VZ5gJZCdZ+cp6d7Vbl0vy9Pv1bNrdTEFqPKsq0km0mFg3L2usmOErd67i0h9/EDE1T5NBkJlo8Uv7fF9zNwbhDpB0Scn+EGvwhUr/sJNoukS/nJEQ4g3AmxL8A1LK30opHwAeEEJ8G7gH+EfgImAnsA4oB14XQrznMZLyGynlE8AT4I7ADuYewbChKj9gZ9Qz6MBsNLAkPxlbnNFr+zvWlFOts63dmYbBILj13DJuPdd3qIYQgseuXsZ9P98REYc0MOIKSPPaJd0BknYfoQHRJpqVTPyapkkp10spF3l5/dbj0heBy7Xjm4BfSTeHgHpgnsf1x4Dxq44F2jndcNHC7KAWGNv7htl9rIutDe3YLRMfrIV5SVyjA8kGhZvi9AR+el2132t9c7ISuXFVCb+47SzuO3/OtFO9YGKh0uLDr0AZDNGsHhKy+xVCzJFSjlq8ATigHR8FzgfeE0JkA5XA4fFtpZQtQohuIcRZwGbcU73/F6pN4SQ90cKSghR2NQYntOaSMCd7YuHGYOptKSKL3WrmCysKeOGvR73+PSHOyB1ryrm8uoBMuwWz5rjqytL5+wvm8v7BNu5+6UOvciPHuwZYVpRCS9cgx6eYuo+nuWuAutI0XFJS39ZHW2/sUoRcLhmVQhAhJ8oKIV7B7WhcwBHgDinlMW2L/1kgF3f+6L9JKV/Q2uzUdt8QQtRo19mA/wHuldMYFalEWV98//cH+Pe3P/X7+rKMBE72DNIz5MRuMVKRZWeHhzNbUpDMty6exx/2Hefbl8xXybE6YVtDO19/eRcNp9wxQLnJVh7asIi1lZmYphk5HTnVx3VPb+Fou/f4oWArzeSn2GjtGfSrPHckWFWezgULsrnurOJp+8AbKms/jLyx/wS3PO/f+5kMkJZgwSAE+ak2TvYM+Xw4R3ngkvl8fnk+TpckyWrGFqccUyz54FAbP37rEO19wzx3cy3ZSf4nxA47XFz79OZJkiGV2Yl87CVS218W5yex51hsJWyvqSti4+cXB9zOX2ekj1UynbMof2rJ2PE4XJBpt7CvuXtS7SxfbHztIza+9hHgzvq//+J5anE7hqyuyGB1RQZSyoCj4eNMBgpT4yc5o1Bz6PYc62ZZYcrYCLsoLZ7CNBtzsuzsbOxkf0s3I04XSVZzWCrImAyTyzi9uPko7X3D1JWmcWME4rmUM/KD1AAzrPc1d5OWYKa9L/CHYkt9O1968q9cXl1ASXo8t51bHvA9FOEhmGQAh9PFVi+liVq6B6kpSWVXY2fQ061mrVLNxQtz+OcNC8nSRmxOlxyLqRJCcPOzW/nTgVbSEuK4tq4Ip5Qcau2loa2f7sERBO41so/H5U/mp9hYVpTC1y6s5FTvEAMjTm7/2fZJO4H/s/c4u5u6uOas4rF1s3ChnJEfWExGkm0mugYcWEyCrCQrHb3D9HrZsk2ymZifkxSS1OiQw8VLm4+Sk2RVzugMo3fI4XVa3tg+QN+gg0y7xW9BNk9O9Q3xoy9WTZL58Izav/WcMlZXZHDJ4pxJ2tajjqtnyMET7xymvX+YJfnJfKGmcGyRejTyetO9Z/N/fr2Hvx5ux2IyUFWYQk6ylTlZiQyOOJUzihXLilJp6hggLd7MlgZ3EuSo0H6W3cKRsQVPW9g0j493D1Lf1hf1sHxF8Bz2KCpgNoqxkdCcbHtIz4bDhV9FMFeWp7Oy3HvpoVHHlWQ18/WLKqe8T1lmIj/7ch1vfnSCVRUZJFkjm7KinJGfHD7ZR0vXAIda3Q/WwIiTrQ3uyg/jv+kKUsMXJJadZAlrLS9F5Hn745Njx+WZCXT2j1Cemcigwxnyl5TZKIiLcliI2Wjg4kW5UXkv5Yz8wOWStHQN+DXXt3ipfRYMVrOBd75xntryP4M4cqqPf3/rdHVim9nIp319nOoLz0j5y2eXzeiyVSr6zg9ae4b8XnT0N6htOgZHXPzHO4env1ChG4rTE3j5jpVjukrBalp5w24xcfPZJWG7nx5RzsgPPg5AtdGfOb2//G5Pc9jupYgOy4pSeezqZVRkJUwbX+YvZRkJ3Hf+HLLsM3vKrqZpfhCIIt/m+nZqilPZdiTwSqKXLM4hLSGOPU1dtPYM8dzNtQHfQxFbGtv7+dGbBxl2hBZMbDIIblhVwrlzM1kzNzNM1ukb5Yz8oLXb/4J9ANuOdFBbkkbDqb6Aiv2dVZbO9StLgOjlAynCy66mTl4PsRDoldUF/Otli4NKvTiTUc7IDz5qCTwMf4sW+FaaEU9KfBw7jk6daGsQTNA3Uo7ozOQXWxuDanfX2nLWzcuivW+YCxZkz0oddOWMpqGzf5h3Pjk5/YU+qG/rB/qZm504qYqExWTgooU53HZuGcXp8WOCX4ozl6tWFLL5cHvAVYkbTvVRUxL74pSxRDmjaXjk9U8YcgRf7noUz3zkl26po7okFYtJbd3PJD63JI+8FBvXPrXZq6haosXEVSsKuWRxDj984yDvHWyjIiuRu88LTuJ4JqGckRcGhp089Lv9fNTSPe30yh+MBoHJKMi2W2jrG+bCBdksL1aOaKayrDCF+blJbPeyiXHPugruWONO8Xn+5lo+OdFLSUa8ehZQzsgrd724nbc+Dn5q5smivKQxHZuHNizkOm2RWjHz+M2OYzz53mH2NU9eZ7xx1cTqJUIIKoMoazVTmV3L9X7S1DG9eLq/FKfH0zvkID/FxmcW5VBdPLvXBWY6ayszvZabWlKQzP0Xz4t6OseZhBoZeSEvxcbB1uCFsEYpy0ggJd7MA5+dT1VhqqqJNgtIiY/j+ZtrueapzRNGR8uLUpVo3jQoN+2Fm1aXhOU++ak2fnzNcqqL05QjmkWkxMfx6j1n842LKpmfmwQQkFrkbEU5Iy+sKs8IuW776op0/vWyxZP0ZBSzA6NBcPd5FTxz4wouW5Yfti+4mYxyRl6IMxm4cKG3MnH+8dX1c3nhy3VhlRNRnJnkJFt55Koqpb7gB2rNyAvDDhcn/NSv9uQHVy7l8uqC6S9UKBQTUCMjL/xm5zF+vSPwWpKXLM5RjkihCBLljLyQaAl8wGgyCP5u/dwIWKNQzA6UM/JgYNjJK9ubAm7ncEmufWoz9/18B++GkMumUMxWlDMah8Pp4qIfvsubB1qDat/aM8Sru5r52su7GByZnJekUCh8o5zROExGA4EoN6yuSKcsI2FSDFFNcaraPVEoAkTtpnlQmW0fKzvkC7NRcEV1Id/6zDySbWYa2/u5/5Xd/PnTU4C72ueI0xX2ulIKxUxGOSMPXJ5aHx7Ulabx2NXLJkTUFqbF85Nrq7npmS18eLST/3j3MF0DI1xeXcCKWa5Ro1D4i/rq9iAj0eLzbxctzOalW8/yGtqfbDPzwi11XFqVB8B/bW3k9p9tZ1dj6BIkCsVsQDkjDzyTGY0GwdzsRK6sLuC2c8unzDGLjzNx+5rT5ai/87cLWVo4c+tcKRThJKRpmhDiIWAD4AJagRullM1CiGTgBaBIe4//K6V8xkv7t4FcYFSz40IpZXBbWWHC5Zo4TfvahXO5a63/KnyneofHjg+f7FXC+gqFn4Q6MnpYSrlESlkFbAIe1M7fDeyXUi4F1gI/EELE+bjHNVLKKu0VU0cEcJGWk3ZWWRoPX7GESwIs7btp9+laZz9842BYtZEUiplMSCMjKeV4ObsEYHRYIQG7cJc4SATagfCV14wgK8vTee2+c8hLsZIS78t/+uautRUcbO1lYNhJpt1C/8gZ8c9WKGJOyLtpQoiNwPVAF3Cedvpx4FWgGbADV0kpfanaPyOEcAKvAN+VcprtrAgjhGBBXlLQ7YvS43nlzlVhtEihmB1MO00TQrwhhNjr5bUBQEr5gJSyEHgRuEdrdhGwE8gDqoDHhRDePuHXSCkXA+dor+umsOM2IcQ2IcS2kydVuoVCMdOY1hlJKddLKRd5ef3W49IXgcu145uAX0k3h4B6YJ6Xex/TfvYALwE+6zlLKZ+QUtZIKWsyM2dHuV+FYjYR0gK2EGLOuF83AAe046PA+do12UAlcNijrUkIkaEdm4HPAXtDsUehUJy5hLpm9G9CiErcW/tHgDu08w8Bzwoh9gACuF9K2QYghNip7b5ZgD9ojsgIvAE8GaI9CoXiDCXU3bTLfZxvBi708bcq7WcfUB3K+ysUipmDisBWKBS6QDkjhUKhC5QzUigUukA5I4VCoQtEjAOeg0IIcRL37l0kyQDaIvweoaDsCw0926dn2yBw+4qllNMGB56RzigaCCG2SSlrYm2HL5R9oaFn+/RsG0TOPjVNUygUukA5I4VCoQuUM/LNE7E2YBqUfaGhZ/v0bBtEyD61ZqRQKHSBGhkpFApdoJyRQqHQBbPaGQkhfiGE2Km9GoQQO7XztePO7xJCfN5H+2eFEPXjrq3SmX2lQojNQohD2r0C19EN3LYLhBDbhRB7tJ/rfLT/jhDi2Lh7XBIu28JkX5oQ4nUhxEHtZ2qU7EsXQrwlhOgVQjw+RftY9Z+/9gXef1JK9XKvm/0AeFA7jgdM2nEu7sonJi9tngWu0LF9/w18UTv+KXBnFGxbBuRpx4uAYz7afAf4egz6zl/7vg98Szv+FvC9KNmXAJyNW47n8SnaxKr//LUv4P6b1SOjUbTCAV8Afg4gpeyXUo4q6Vs5XWggJgRjn9ZmHfBL7dRzwKVRsG2HdEvIAOwDbEII35UxI0wI9m3A3WcQob7zYV+flPJ9YDAS7xcoIdgXcP8pZ+TmHOCElPLg6AkhRJ0QYh+wB7hj3Iffk41CiN1CiEcj+KELxr50oHPc+SYgPxq2jeNy4EMp5ZCPtvdoffef4Z4GhcG+bClli3Z8HMiOgX3TEev+m4qA+2/GOyMxTUEBjavRPP8oUsrNUsqFwArg20KIyTWt4du4tb1XAGnA/TqzLySCtU1ruxD4HnC7j9v/BCjHXbChBfdUQE/2jSHdc42AR8eh2OcHMe2/QPC7/6Ix59TzC7fa5QmgYIpr/gTUTHOftcAmvdiHW+63jdNrSyuBP0TDNqAA+ARY7ed9SoC90eo7f+wDPgZyteNc4ONo/t8CNzLFmkws+88f+4Lpvxk/MvKD9cABKWXT6Anh3oUyacfFuEc/DZ4NhRC52k+Be04ciYICQdkn3U/BW8AV2qkbAM+KLpGwLQX4He7Fyw98NRztO43PE72+88s+3HUQBREMAAAAvElEQVT/btCOI9F3Xu3zl1j1XwAE3n/h9qZn2gv3jtgdHueuw724uRP4ELh03N9e4/RuzJ9wr9nsBV4AEnVmXxmwBTgEvAxYomDbPwB9mm2jryztb0+hjeCAn2l9t1t7cHOj1Hf+2pcOvAkcxF0sIi0a9mnnG3BXYe7Fvda3QC/9F4B9AfefSgdRKBS6QE3TFAqFLlDOSKFQ6ALljBQKhS5QzkihUOgC5YwUCoUuUM5IoVDoAuWMFAqFLvhfj4WpTLWcmksAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "biobio.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['Mulchén', 'Nacimiento', 'Negrete', 'Quilaco', 'Quilleco',\n",
       "       'San Rosendo', 'Santa Bárbara', 'Tucapel', 'Yumbel', 'Concepción',\n",
       "       'Florida', 'Hualpén', 'Hualqui', 'Lota', 'Penco', 'Santa Juana',\n",
       "       'Tomé', 'Bulnes', 'Chillán', 'Chillán Viejo', 'Cobquecura',\n",
       "       'Coelemu', 'Chiguayante', 'San Pedro de la Paz', 'El Carmen',\n",
       "       'Ninhue', 'Ñiquén', 'Pemuco', 'Pinto', 'Portezuelo', 'Quillón',\n",
       "       'Quirihue', 'Ranquil', 'San Carlos', 'San Fabián', 'San Ignacio',\n",
       "       'San Nicolás', 'Treguaco', 'Yungay', 'Coihueco', 'Arauco',\n",
       "       'Cañete', 'Contulmo', 'Curanilahue', 'Los Alamos', 'Tirúa',\n",
       "       'Alto Biobío', 'Antuco', 'Cabrero', 'Laja', 'Los Angeles'],\n",
       "      dtype=object)"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biobio['NOM_COM'].unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[:, 'centroid'] = biobio.loc[:,'geometry'].centroid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "def assign_x(centroid):\n",
    "    return float(centroid.x)\n",
    "\n",
    "def assign_y(centroid):\n",
    "    return float(centroid.y)\n",
    "biobio.loc[:, 'centroid_x'] = biobio['centroid'].apply(assign_x)\n",
    "biobio.loc[:, 'centroid_y'] = biobio['centroid'].apply(assign_y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.drop('centroid',  axis=1, inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "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>NOM_REG</th>\n",
       "      <th>NOM_PROV</th>\n",
       "      <th>NOM_COM</th>\n",
       "      <th>SHAPE_LENG</th>\n",
       "      <th>DIS_ELEC</th>\n",
       "      <th>CIR_SENA</th>\n",
       "      <th>COD_COMUNA</th>\n",
       "      <th>SHAPE_Le_1</th>\n",
       "      <th>SHAPE_Area</th>\n",
       "      <th>geometry</th>\n",
       "      <th>centroid_x</th>\n",
       "      <th>centroid_y</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>120</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Mulchén</td>\n",
       "      <td>311947.908393</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8305</td>\n",
       "      <td>311947.908133</td>\n",
       "      <td>1.922636e+09</td>\n",
       "      <td>POLYGON ((-72.43334352332262 -37.5619680689156...</td>\n",
       "      <td>-72.097023</td>\n",
       "      <td>-37.838358</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>121</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Nacimiento</td>\n",
       "      <td>191416.344754</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8306</td>\n",
       "      <td>191416.344324</td>\n",
       "      <td>9.117771e+08</td>\n",
       "      <td>POLYGON ((-72.72013571650028 -37.3025229857047...</td>\n",
       "      <td>-72.823564</td>\n",
       "      <td>-37.485523</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>122</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Negrete</td>\n",
       "      <td>77213.548250</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8307</td>\n",
       "      <td>77213.548309</td>\n",
       "      <td>1.556746e+08</td>\n",
       "      <td>POLYGON ((-72.59553424428385 -37.5310507766004...</td>\n",
       "      <td>-72.576349</td>\n",
       "      <td>-37.608089</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilaco</td>\n",
       "      <td>282388.766764</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8308</td>\n",
       "      <td>282388.767338</td>\n",
       "      <td>1.126256e+09</td>\n",
       "      <td>POLYGON ((-72.07607421008562 -37.6530209035402...</td>\n",
       "      <td>-71.704976</td>\n",
       "      <td>-37.960748</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>124</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilleco</td>\n",
       "      <td>240420.559462</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8309</td>\n",
       "      <td>240420.559976</td>\n",
       "      <td>1.121994e+09</td>\n",
       "      <td>POLYGON ((-72.03043817152037 -37.2290063644018...</td>\n",
       "      <td>-71.861352</td>\n",
       "      <td>-37.436655</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                NOM_REG NOM_PROV     NOM_COM     SHAPE_LENG  DIS_ELEC  \\\n",
       "120  Región del Bío-Bío  Bío-Bío     Mulchén  311947.908393        47   \n",
       "121  Región del Bío-Bío  Bío-Bío  Nacimiento  191416.344754        47   \n",
       "122  Región del Bío-Bío  Bío-Bío     Negrete   77213.548250        47   \n",
       "123  Región del Bío-Bío  Bío-Bío     Quilaco  282388.766764        47   \n",
       "124  Región del Bío-Bío  Bío-Bío    Quilleco  240420.559462        47   \n",
       "\n",
       "     CIR_SENA  COD_COMUNA     SHAPE_Le_1    SHAPE_Area  \\\n",
       "120        13        8305  311947.908133  1.922636e+09   \n",
       "121        13        8306  191416.344324  9.117771e+08   \n",
       "122        13        8307   77213.548309  1.556746e+08   \n",
       "123        13        8308  282388.767338  1.126256e+09   \n",
       "124        13        8309  240420.559976  1.121994e+09   \n",
       "\n",
       "                                              geometry  centroid_x  centroid_y  \n",
       "120  POLYGON ((-72.43334352332262 -37.5619680689156...  -72.097023  -37.838358  \n",
       "121  POLYGON ((-72.72013571650028 -37.3025229857047...  -72.823564  -37.485523  \n",
       "122  POLYGON ((-72.59553424428385 -37.5310507766004...  -72.576349  -37.608089  \n",
       "123  POLYGON ((-72.07607421008562 -37.6530209035402...  -71.704976  -37.960748  \n",
       "124  POLYGON ((-72.03043817152037 -37.2290063644018...  -71.861352  -37.436655  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biobio.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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>NOM_REG</th>\n",
       "      <th>NOM_PROV</th>\n",
       "      <th>NOM_COM</th>\n",
       "      <th>SHAPE_LENG</th>\n",
       "      <th>DIS_ELEC</th>\n",
       "      <th>CIR_SENA</th>\n",
       "      <th>COD_COMUNA</th>\n",
       "      <th>SHAPE_Le_1</th>\n",
       "      <th>SHAPE_Area</th>\n",
       "      <th>geometry</th>\n",
       "      <th>centroid_x</th>\n",
       "      <th>centroid_y</th>\n",
       "      <th>id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Mulchén</td>\n",
       "      <td>311947.908393</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8305</td>\n",
       "      <td>311947.908133</td>\n",
       "      <td>1.922636e+09</td>\n",
       "      <td>POLYGON ((-72.43334352332262 -37.5619680689156...</td>\n",
       "      <td>-72.097023</td>\n",
       "      <td>-37.838358</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Nacimiento</td>\n",
       "      <td>191416.344754</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8306</td>\n",
       "      <td>191416.344324</td>\n",
       "      <td>9.117771e+08</td>\n",
       "      <td>POLYGON ((-72.72013571650028 -37.3025229857047...</td>\n",
       "      <td>-72.823564</td>\n",
       "      <td>-37.485523</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Negrete</td>\n",
       "      <td>77213.548250</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8307</td>\n",
       "      <td>77213.548309</td>\n",
       "      <td>1.556746e+08</td>\n",
       "      <td>POLYGON ((-72.59553424428385 -37.5310507766004...</td>\n",
       "      <td>-72.576349</td>\n",
       "      <td>-37.608089</td>\n",
       "      <td>2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilaco</td>\n",
       "      <td>282388.766764</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8308</td>\n",
       "      <td>282388.767338</td>\n",
       "      <td>1.126256e+09</td>\n",
       "      <td>POLYGON ((-72.07607421008562 -37.6530209035402...</td>\n",
       "      <td>-71.704976</td>\n",
       "      <td>-37.960748</td>\n",
       "      <td>3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilleco</td>\n",
       "      <td>240420.559462</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8309</td>\n",
       "      <td>240420.559976</td>\n",
       "      <td>1.121994e+09</td>\n",
       "      <td>POLYGON ((-72.03043817152037 -37.2290063644018...</td>\n",
       "      <td>-71.861352</td>\n",
       "      <td>-37.436655</td>\n",
       "      <td>4</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              NOM_REG NOM_PROV     NOM_COM     SHAPE_LENG  DIS_ELEC  CIR_SENA  \\\n",
       "0  Región del Bío-Bío  Bío-Bío     Mulchén  311947.908393        47        13   \n",
       "1  Región del Bío-Bío  Bío-Bío  Nacimiento  191416.344754        47        13   \n",
       "2  Región del Bío-Bío  Bío-Bío     Negrete   77213.548250        47        13   \n",
       "3  Región del Bío-Bío  Bío-Bío     Quilaco  282388.766764        47        13   \n",
       "4  Región del Bío-Bío  Bío-Bío    Quilleco  240420.559462        47        13   \n",
       "\n",
       "   COD_COMUNA     SHAPE_Le_1    SHAPE_Area  \\\n",
       "0        8305  311947.908133  1.922636e+09   \n",
       "1        8306  191416.344324  9.117771e+08   \n",
       "2        8307   77213.548309  1.556746e+08   \n",
       "3        8308  282388.767338  1.126256e+09   \n",
       "4        8309  240420.559976  1.121994e+09   \n",
       "\n",
       "                                            geometry  centroid_x  centroid_y  \\\n",
       "0  POLYGON ((-72.43334352332262 -37.5619680689156...  -72.097023  -37.838358   \n",
       "1  POLYGON ((-72.72013571650028 -37.3025229857047...  -72.823564  -37.485523   \n",
       "2  POLYGON ((-72.59553424428385 -37.5310507766004...  -72.576349  -37.608089   \n",
       "3  POLYGON ((-72.07607421008562 -37.6530209035402...  -71.704976  -37.960748   \n",
       "4  POLYGON ((-72.03043817152037 -37.2290063644018...  -71.861352  -37.436655   \n",
       "\n",
       "   id  \n",
       "0   0  \n",
       "1   1  \n",
       "2   2  \n",
       "3   3  \n",
       "4   4  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biobio.reset_index(drop=True, inplace = True)\n",
    "biobio.loc[:, 'id'] = biobio.index.values\n",
    "biobio.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>NOM_REG</th>\n",
       "      <th>NOM_PROV</th>\n",
       "      <th>NOM_COM</th>\n",
       "      <th>SHAPE_LENG</th>\n",
       "      <th>DIS_ELEC</th>\n",
       "      <th>CIR_SENA</th>\n",
       "      <th>COD_COMUNA</th>\n",
       "      <th>SHAPE_Le_1</th>\n",
       "      <th>SHAPE_Area</th>\n",
       "      <th>geometry</th>\n",
       "      <th>centroid_x</th>\n",
       "      <th>centroid_y</th>\n",
       "      <th>id</th>\n",
       "      <th>centroid_manual_x</th>\n",
       "      <th>centroid_manual_y</th>\n",
       "      <th>add_line</th>\n",
       "      <th>line_loc</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Mulchén</td>\n",
       "      <td>311947.908393</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8305</td>\n",
       "      <td>311947.908133</td>\n",
       "      <td>1.922636e+09</td>\n",
       "      <td>POLYGON ((-72.43334352332262 -37.5619680689156...</td>\n",
       "      <td>-72.097023</td>\n",
       "      <td>-37.838358</td>\n",
       "      <td>0</td>\n",
       "      <td>-72.097023</td>\n",
       "      <td>-37.838358</td>\n",
       "      <td>False</td>\n",
       "      <td>down</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Nacimiento</td>\n",
       "      <td>191416.344754</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8306</td>\n",
       "      <td>191416.344324</td>\n",
       "      <td>9.117771e+08</td>\n",
       "      <td>POLYGON ((-72.72013571650028 -37.3025229857047...</td>\n",
       "      <td>-72.823564</td>\n",
       "      <td>-37.485523</td>\n",
       "      <td>1</td>\n",
       "      <td>-72.823564</td>\n",
       "      <td>-37.485523</td>\n",
       "      <td>False</td>\n",
       "      <td>down</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Negrete</td>\n",
       "      <td>77213.548250</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8307</td>\n",
       "      <td>77213.548309</td>\n",
       "      <td>1.556746e+08</td>\n",
       "      <td>POLYGON ((-72.59553424428385 -37.5310507766004...</td>\n",
       "      <td>-72.576349</td>\n",
       "      <td>-37.608089</td>\n",
       "      <td>2</td>\n",
       "      <td>-72.576349</td>\n",
       "      <td>-37.608089</td>\n",
       "      <td>False</td>\n",
       "      <td>down</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilaco</td>\n",
       "      <td>282388.766764</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8308</td>\n",
       "      <td>282388.767338</td>\n",
       "      <td>1.126256e+09</td>\n",
       "      <td>POLYGON ((-72.07607421008562 -37.6530209035402...</td>\n",
       "      <td>-71.704976</td>\n",
       "      <td>-37.960748</td>\n",
       "      <td>3</td>\n",
       "      <td>-71.704976</td>\n",
       "      <td>-37.960748</td>\n",
       "      <td>False</td>\n",
       "      <td>down</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Región del Bío-Bío</td>\n",
       "      <td>Bío-Bío</td>\n",
       "      <td>Quilleco</td>\n",
       "      <td>240420.559462</td>\n",
       "      <td>47</td>\n",
       "      <td>13</td>\n",
       "      <td>8309</td>\n",
       "      <td>240420.559976</td>\n",
       "      <td>1.121994e+09</td>\n",
       "      <td>POLYGON ((-72.03043817152037 -37.2290063644018...</td>\n",
       "      <td>-71.861352</td>\n",
       "      <td>-37.436655</td>\n",
       "      <td>4</td>\n",
       "      <td>-71.861352</td>\n",
       "      <td>-37.436655</td>\n",
       "      <td>False</td>\n",
       "      <td>down</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "              NOM_REG NOM_PROV     NOM_COM     SHAPE_LENG  DIS_ELEC  CIR_SENA  \\\n",
       "0  Región del Bío-Bío  Bío-Bío     Mulchén  311947.908393        47        13   \n",
       "1  Región del Bío-Bío  Bío-Bío  Nacimiento  191416.344754        47        13   \n",
       "2  Región del Bío-Bío  Bío-Bío     Negrete   77213.548250        47        13   \n",
       "3  Región del Bío-Bío  Bío-Bío     Quilaco  282388.766764        47        13   \n",
       "4  Región del Bío-Bío  Bío-Bío    Quilleco  240420.559462        47        13   \n",
       "\n",
       "   COD_COMUNA     SHAPE_Le_1    SHAPE_Area  \\\n",
       "0        8305  311947.908133  1.922636e+09   \n",
       "1        8306  191416.344324  9.117771e+08   \n",
       "2        8307   77213.548309  1.556746e+08   \n",
       "3        8308  282388.767338  1.126256e+09   \n",
       "4        8309  240420.559976  1.121994e+09   \n",
       "\n",
       "                                            geometry  centroid_x  centroid_y  \\\n",
       "0  POLYGON ((-72.43334352332262 -37.5619680689156...  -72.097023  -37.838358   \n",
       "1  POLYGON ((-72.72013571650028 -37.3025229857047...  -72.823564  -37.485523   \n",
       "2  POLYGON ((-72.59553424428385 -37.5310507766004...  -72.576349  -37.608089   \n",
       "3  POLYGON ((-72.07607421008562 -37.6530209035402...  -71.704976  -37.960748   \n",
       "4  POLYGON ((-72.03043817152037 -37.2290063644018...  -71.861352  -37.436655   \n",
       "\n",
       "   id  centroid_manual_x  centroid_manual_y  add_line line_loc  \n",
       "0   0         -72.097023         -37.838358     False     down  \n",
       "1   1         -72.823564         -37.485523     False     down  \n",
       "2   2         -72.576349         -37.608089     False     down  \n",
       "3   3         -71.704976         -37.960748     False     down  \n",
       "4   4         -71.861352         -37.436655     False     down  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biobio.loc[:, 'centroid_manual_x'] = biobio.loc[:, 'centroid_x']\n",
    "biobio.loc[:, 'centroid_manual_y'] = biobio.loc[:, 'centroid_y']\n",
    "biobio.loc[:, 'add_line'] = False\n",
    "biobio.loc[:, 'line_loc'] = 'down'\n",
    "biobio.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# MANUAL CENTROIDS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[13, 'centroid_manual_x'] = -73.19470720619215\n",
    "biobio.loc[13, 'centroid_manual_y'] = -37.13201571832683\n",
    "\n",
    "biobio.loc[13, 'add_line'] = False\n",
    "biobio.loc[13, 'line_loc'] = 'right'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[26, 'centroid_manual_x'] = -72.06804310767384\n",
    "biobio.loc[26, 'centroid_manual_y'] = -36.2231461440657"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[31, 'centroid_manual_x'] = -72.53895343646452\n",
    "biobio.loc[31, 'centroid_manual_y'] = -36.1627659414129"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[20, 'centroid_manual_x'] = -72.78541117863536\n",
    "biobio.loc[20, 'centroid_manual_y'] = -36.30121608933768"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[21, 'centroid_manual_x'] = -72.89983798750039\n",
    "biobio.loc[21, 'centroid_manual_y'] = -36.46066139037599"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[18, 'centroid_manual_x'] = -72.02403279657192\n",
    "biobio.loc[18, 'centroid_manual_y'] = -36.56854387519827"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[22, 'centroid_manual_x'] = -73.11328813065357\n",
    "biobio.loc[22, 'centroid_manual_y'] = -36.98802202370379\n",
    "\n",
    "biobio.loc[22, 'add_line'] = True\n",
    "biobio.loc[22, 'line_loc'] = 'up'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[23, 'centroid_manual_x'] = -73.25632164173484\n",
    "biobio.loc[23, 'centroid_manual_y'] = -36.854320539197886\n",
    "\n",
    "biobio.loc[23, 'add_line'] = True\n",
    "biobio.loc[23, 'line_loc'] = 'right'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[11, 'centroid_manual_x'] = -73.26732421951033\n",
    "biobio.loc[11, 'centroid_manual_y'] = -36.72038483392178\n",
    "\n",
    "biobio.loc[11, 'add_line'] = True\n",
    "biobio.loc[11, 'line_loc'] = 'right'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[14, 'centroid_manual_x'] = -73.07587936621692\n",
    "biobio.loc[14, 'centroid_manual_y'] = -36.65685992884627\n",
    "\n",
    "biobio.loc[14, 'add_line'] = True\n",
    "biobio.loc[14, 'line_loc'] = 'right'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[49, 'centroid_manual_y'] = -37.36303894687814"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[24, 'centroid_manual_x'] = -71.80591122640529"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.loc[32, 'centroid_manual_y'] = -36.6956869399083"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio_geo_ori = biobio.loc[:, 'geometry'].copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio['geometry'] = biobio_geo_ori.simplify(0.001, preserve_topology=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio.to_file('data/biobio.geojson', driver='GeoJSON')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "biobio_simple = biobio.copy()\n",
    "biobio_simple['geometry'] = biobio_simple['geometry'].simplify(0.02, preserve_topology=True)\n",
    "biobio_simple.to_file('data/biobio_simplified.geojson', driver='GeoJSON')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "geoPropDataExploration",
   "language": "python",
   "name": "geopropdataexploration"
  },
  "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.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}