{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Basic plotting of geospatial data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Ensuring high resolution plots from mac device\n",
    "%config InlineBackend.figure_format = 'retina'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Plotting vector data\n",
    "Vector data are data represented by a series of vertices, sometimes connected by edges. You typically have 4 types:\n",
    "1. Points. Single vertices (e.g. a latitude and longitude coordinate pair that represents the location of an ATM machine)\n",
    "2. Lines. Sets of vertices that are connected by edges (e.g. 3 latitude/longitude pairs that define the location of  a road).\n",
    "3. Polygons. Sets of vertices connected by edges where the last vertex and the first vertex are the same (e.g. 4 latitude/longitude pairs that define the outline of a building)\n",
    "4. Multi-polygons. A collection of polygons (e.g. a university represented as a collection of polygon buildings or a building with an inner courtyard defined by an outer and inner polygon)\n",
    "\n",
    "In this example, we'll start by walking through how to plotpolygons representing a grid accross the state of Connecticut."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "import geopandas as gpd\n",
    "import rasterio\n",
    "from rasterio.plot import show\n",
    "from rasterio.features import rasterize\n",
    "\n",
    "plt.close('all')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABNMAAAOKCAYAAACxtxb5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAABYlAAAWJQFJUiTwAACAR0lEQVR4nOz9e7S1Z1kf+n8vRA6igUWrP62KQGoQqWybcNCkcoiYDcihSrLJzxYp/IrHggId1F1EDMreRTAi4KFVUdQ9mmAcO1REiwIRSkoYJFWGEgFDMhBFEV7llIOE3L8/5lxxMt+53vdZ653PWvd85uczxhoz65nPfc3v9fJS4tV530+11gIAAAAAnNwdjjoAAAAAAGwKwzQAAAAAGMgwDQAAAAAGMkwDAAAAgIEM0wAAAABgIMM0AAAAABjIMA0AAAAABjJMAwAAAICBDNMAAAAAYCDDNAAAAAAYyDANAAAAAAa641EH4MSq6vokpyW54YijAAAAAEzFvZN8orV2n/0uNEzr32l3vetd73n/+9//nkcdBAAAAGAKrr322tx0000HWmuY1r8b7n//+9/z6quvPuocAAAAAJNw1lln5ZprrrnhIGudmQYAAAAAAxmmAQAAAMBAhmkAAAAAMJBhGgAAAAAMZJgGAAAAAAMZpgEAAADAQIZpAAAAADCQYRoAAAAADGSYBgAAAAADGaYBAAAAwECGaQAAAAAwkGEaAAAAAAxkmAYAAAAAAxmmAQAAAMBAhmkAAAAAMJBhGgAAAAAMZJgGAAAAAAMZpgEAAADAQIZpAAAAADCQYRoAAAAADGSYBgAAAAADGaYBAAAAwECGaQAAAAAwkGEaAAAAAAxkmAYAAAAAA61lmFYzT6+qd1TVJ6vqxqr6X1X1rKr6vKV7711V7QQ/l5zgc55aVe+sqk9V1cer6oqqetwJ7r9rVV1UVe+tqpur6iNV9dqquv8J1nxFVb26qv6yqm6pqhuq6uVVtXOCNWdX1Ruq6ti893dX1Q8u9w4AAADAZrvjmuq8JslTknwkyaVJPp3kUUl+OsnDquqC1lpbWvNHSS5fUeuPV31AVb0syXOTfCjJLyS5U5ILk/xWVT2ztfaqpfvvnOT3kpyT5F3zLF+Z5IIk31pV57bWrlpac3qSK5N8SZLXJfnTJA9J8gNJHl1V57TWPra05olJfjPJzfPejyV5fJKfmn/2Bav6AQAAAGDznPIwrar+ZWaDtOuTPKS19tH59c9P8tokT0ry1CS/srT0D1trPzrwM87ObJB2XZIHt9b+dn79pUmuTvKyqnp9a+2GhWXPyWyYdVmSJ7fWbpuvuTSzId6rq+rrdq/P/Wxmg7RntdZeufD5Fyd5dpIXJ/meheunZTbY+2ySR7TW3jW//oIkb05yflVd2Frb89t2AAAAAGyOdWzz/Pb560/uDtKSpLX2mSQvmP/6zFP8jN0B1ot3B2nzz7ghyc8kuXOSp+1er6paWPO8xYFZa+11Sd6W5GuTPHxhzX2TnJdkt+aiF2b2bbunVNXdFq6fn+SLk1yyO0ibf8bNSX54/uv37q9VAAAAAHq1jmHal85fP7Divd1rZ1bVPZbe+ydV9d1V9R/nrw88wWecO3/93RXv/c7SPUlyepJ7JXlfa+36gWt2//mNS99WS2vtk0nenuQLknzDwFxvTXJjkrPnW04BAAAA2HDrODNt99to91nx3n0X/vlrkrxj4fdvmf/crqquSPLU1toHF67dLcmXJ/lUa+3DKz7j/fPXMxau3W/++r49Mh90zXnzNW862ZrW2q1VdX2SB2T253DtHnWTJFV19R5vfc2J1gEAAABweNbxzbTXz1+fU1X33L1YVXdMctHCfbtPw7wxyY8lOWt+bSez7ZZvSfKIJG9a2kp59/nrx/f4/N3r99iQNQAAAABsqHV8M+2SJP86yWOSvKeq/ltmA7NHZbbd8v1JvjqzQ/rTWvtIkh9ZqvHWqjovyf9I8tAk/zazp2/ux/LTQk+kelvTWjtrZYHZN9bO3MdnAgAAADCSU/5m2vx8sSck+fdJ/iqzJ3s+PcmHkvyLJB+b3/qRk9S5Nckvzn992MJbu9/uuntWW/XtsJOtOe0I1wAAAACwodaxzTOttVtbaz/ZWvv61tpdW2untdYeneQ9Sb4+yU1J/mRAqb+Zv96+zbO19ukkf5HkC6vqy1as+er56+K5Ze+dv56R1UZfM9/mep8kt2b1wxkAAAAA2DBrGaadwFOS3CXJa1trnxlw/+6TMpeHT2+evz56xZrHLN2TJNcl+WCSM6pq1YMRVq15y/z1vKr6nD+XqvqiJOdkNhRcfIjCiXI9LLOnf17ZWrtlxfsAAAAAbJi1DNOq6rQV1x6c5D8l+VSSFy1cf2hV3WnF/ecmefb8119fevvn56/Pr6qdhTX3TvL9SW5J8su711trbWHNTywOx6rqiUm+KbNvzf3BwprrkrwxyW7NRRdl9m25X51/U27XZZk9zfTCqnrQwmfcJcmPz3/9ueVegZmq2rqfe97znif/gwEAAKBb63gAQZL8XlXdlOSPk3wyyQOSPDazIde3t9YWv2n2kiQPqKorMjtXLUkemOTc+T+/oLV25WLx1tqVVXVxkuckeXdVXZbkTkmenOSeSZ7ZWrthKdPFSR6X5PwkV1XVm5LcK8kFmT0g4enz894WfV+SK5O8oqq+Ocm1mT0Q4ZGZbe98/lKuT1TVMzIbql1RVZckOZbZGXL3m1+/9AR/brD1ZrPvmar6nN9P5VoPtfa6BwAAgM21rmHaZUkuzOypnndN8peZPUzgP60Ycv1akm9L8uDMtlt+fpK/TvLaJK9qrb1t1Qe01p5bVe9O8u+SfFeS25Jck+SlrbXXr7j/lqp6VJIfSvIdmX3r7RNJLk/ywtbae1asuW7+DbMXZbZ187FJPpzkFUkuaq0dW7Hm8qp6eGaDtidltq31zzIb/L2iLf9f0gAAAABsrDLr6VtVXX3mmWeeefXVVx91FFi7Hr9Nts5aQ9cBAABwuM4666xcc80117TWztrv2nV9Mw3gQJa3Pa7aBnnQaz3UGrqOfu3s7OTYseO+mAwAAGwpwzTgSPX2bbJ11pJhOhkAAAB2reVpngAAAACwDQzTAAAAAGAgwzQAAAAAGMiZacCR6vGhAeusJcM0M0BvPCgDAODwGKYBR6q3g+rXWUsGGabcjwz9ZNi9BgDA4bDNEwAAAAAGMkwDAAAAgIEM0wAAAABgIGemAUeqx4Pq11lLBhnGqiWDDEOuncw6H1zg3Lb1OpX/bPxnAZvLA2VgMximAUeq14O8p3IwuQz9ZJhaPzL0k+FUa63TVP5Me8lwKnrsRwYZptrPujMA/bPNEwAAAAAGMkwDAAAAgIEM0wAAAABgIGemAUdq+VyIVedEHPRaD7VkkGGsWjLIsI5a6z6bZ0p/pr1mGGpT+pFhezOss9Y2ZIDD4iEYwximAUeqtwNk11lLBhmm3I8M/WQ41VrrNJU/Uxmm2Y8M/WSYWj8y9JNhav1M4d8Npso2TwAAAAAYyDANAAAAAAayzRM4UstfI171teKDXuuhlgwyjFVLBhnWUeugWzn2WjelP1MZ1l9LBhnGqiWDDGPV2sYMB/13g21jmAYcqR7OBJjaGQcy9Jlhav3I0E+GU611UFP+M5Vhmv3I0E+GqfUjQz8ZptbPpv27wTaxzRMAAAAABjJMAwAAAICBDNMAAAAAYCBnpgFHanlP/qo9+ge91kMtGWQYq5YMMqyj1qmcizL1P1MZ1l9LBhnGqiWDDGPV2sYMp/LvBtvEMA04Uj0csDm1A0Nl6DPD1PqRoZ8Mp1rroKb8ZyrDNPuRoZ8MU+tHhn4yTK2fTft3g21imycAAAAADGSYBgAAAAADGaYBAAAAwEDOTAOO1PKe/FV79A96rYdaMsgwVi0ZZFhHrVM5F2Xqf6YyrL+WDDKMVUsGGcaqtY0ZTuXfDbaJYRpwpHo4YHNqB4bK0GeGqfUjQz8ZTrXWQU35z1SGafYjQz8ZptaPDP1kmFo/m/bvBtvENk8AAAAAGMgwDQAAAAAGss0TOBR7fV14+fqq+w56rYdaMsgwVi0ZZBiz1hBT/zOVYf21ZJBhrFoyyDBWrW3McNB/L9g2hmnAoen1TICpnXEgQ58ZptaPDP1kmFo/MvSTYWr9yNBPhqn1I0M/GabWz1Fl4ORs8wQAAACAgQzTAAAAAGAgwzQAAAAAGMiZacChWbX/fvnakHs2qZYMMoxVSwYZxqolgwxj1ZJBhrFqySDDWLW2McOqezieYRpwaHo9YHNqB4bK0GeGqfUjQz8ZptaPDP1kmFo/MvSTYWr9yNBPhqn1c1QZODnbPAEAAABgIMM0AAAAABjIMA0AAAAABnJmGnBoVu2/X7425J5NqiWDDGPVkkGGsWrJIMNYtWSQYaxaMsgwVq1tzLDqHo5nmAYcml4P2JzagaEy9Jlhav3I0E+GqfUjQz8ZptaPDP1kmFo/MvSTYWr9HFUGTs42TwAAAAAYyDANAAAAAAYyTAMAAACAgZyZBhyaVfvvl68NuWeTaskgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVWsbM6y6h+MZpgGHptcDNqd2YKgMfWaYWj8y9JNhav3I0E+GqfUjQz8ZptaPDP1kmFo/R5WBk7PNEwAAAAAGMkwDAAAAgIFs8wQOzaqvDC9fG3LPJtWSQYaxaskgw1i1ZJBhrFoyyDBWLRlkGKuWDMnOzk6OHTt23D3bzjANODS9ngkwtTMOZOgzw9T6kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+uklA8ezzRMAAAAABjJMAwAAAICBDNMAAAAAYCBnpgGHZtV+++VrQ+7ZpFoyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYfX7GKYBh2gTDthcZy0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj+9ZOB4tnkCAAAAwECGaQAAAAAwkGEaAAAAAAzkzDTg0Kzab798bcg9m1RLBhnGqiWDDGPVkkGGsWrJIMNYtWSQYaxaMqx+H8M04BBtwgGb66wlgwxT7keGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+unlwwczzZPAAAAABjIMA0AAAAABjJMAwAAAICBnJkGHJpV++2Xrw25Z5NqySDDWLVkkGGsWjLIMFYtGWQYq5YMMoxVS4bV72OYBhyiTThgc521ZJBhyv3I0E+GqfUjQz8ZptaPDP1kmFo/MvSTYWr99JKB49nmCQAAAAADGaYBAAAAwEC2eQKHZtVXhJevDblnk2rJIMNYtWSQYaxaMsgwVi0ZZBirlgwyjFVLhtXvY5gGHKJNOBNgnbVkkGHK/cjQT4ap9SNDPxmm1o8M/WSYWj8y9JNhav30koHj2eYJAAAAAAMZpgEAAADAQIZpAAAAADCQM9OAQ7Nqv/3ytSH3bFItGWQYq5YMMoxVSwYZxqolgwxj1ZJBhrFqybD6fQzTgEO0CQdsrrOWDDJMuR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK2fXjJwPNs8AQAAAGAgwzQAAAAAGGgtw7SaeXpVvaOqPllVN1bV/6qqZ1XV5+2x5uyqekNVHZvf/+6q+sG97p+veWpVvbOqPlVVH6+qK6rqcSe4/65VdVFVvbeqbq6qj1TVa6vq/idY8xVV9eqq+suquqWqbqiql1fVzgnW7LsXAAAAADZPLe+PPVCRql9N8pQkH0nyW0k+neRRSb42yW8muaAtfFBVPXF+/eYklyY5luTxSe6X5LLW2gUrPuNlSZ6b5ENJLktypyQXJrlnkme21l61dP+dk7wpyTlJ3pXkzUm+MskFSf4+ybmttauW1pye5MokX5LkdUn+NMlDkjwyyXuTnNNa+9jSmn33sh9VdfWZZ5555tVXX30qZeDI2WsPAACwWaoqt91221HHGMVZZ52Va6655prW2ln7XXvKDyCoqn+Z2SDt+iQPaa19dH7985O8NsmTkjw1ya/Mr5+W5BeSfDbJI1pr75pff0FmA6/zq+rC1tolC59xdmaDtOuSPLi19rfz6y9NcnWSl1XV61trNyxEe05mg7TLkjy5tXbbfM2lSS5P8uqq+rrd63M/m9kg7VmttVcufP7FSZ6d5MVJvmfh+r57gW22CQdsrrOWDDJMuR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK2fXjJwvHVs8/z2+etPtvkgLUlaa59J8oL5r89cuP/8JF+c5JI2Hz7N7785yQ/Pf/3epc/YHWC9uM0HafM1NyT5mSR3TvK03es1+097d83z2sLArLX2uiRvy+xbcw9fWHPfJOcl2a256IWZfdvuKVV1t1PsBQAAAIANtY5h2pfOXz+w4r3da2dW1T3m/3zu/PV3V9z/1iQ3Jjm7Zts0d51oze8s3ZMkpye5V5L3tdauH7hm95/f2D7322pprX0yyduTfEGSbxiYa69eVqqqq1f9JPmak60FAAAA4HCc8jbPJLvfRrvPivfuu/DPX5PkHZmdJZYk71u+ubV2a1Vdn+QB87XXzr8J9uVJPtVa+/CKz3j//PWMhWt7fsYprjlvvuZNJ1uzqpc96sLWWPUV4eVrQ+7ZpFoyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYfX7rGeY9vok/98kz6mqS1prx5Kkqu6Y5KKF+3bmr3efv358j3q71+9xwPt7X7NS2+PAu5p9O+3Mk62HTbAJZwKss5YMMky5Hxn6yTC1fmToJ8PU+pGhnwxT60eGfjJMrZ9eMnC8dQzTLknyr5M8Jsl7quq/Zba98VGZbbd8f5KvzuyQ/iF2/5NqJ7zrePu5/yCfcVhrAAAAAOjUKZ+Z1mbniz0hyb9P8leZPdnz6Uk+lORfJPnY/NaPzF93v6119z1KnrZ038nuv/vSfQf5jMNcAwAAAMCGWscDCNJau7W19pOtta9vrd21tXZaa+3RSd6T5OuT3JTkT+a3v3f+esZynZptDb1Pklszf3hBa+3TSf4iyRdW1Zet+Pivnr8unlu252cc1ppVvQAAAACw2daxzfNEnpLkLkle01r7zPzam5P8qySPTvJfl+5/WGZPzHxra+2Whetvntd6dJJfXlrzmIV7dl2X5INJzqiq+7Tjn+i5as1b5q/nVdUd2sITPavqi5Kck9lQ8B1LufbbC2ytVfvtl68NuWeTaskgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVUuG1e+zpm+mVdVpK649OMl/SvKpJC9aeOuyzJ4AemFVPWjh/rsk+fH5rz+3VO7n56/Pr6qdhTX3TvL9SW7JwpCtzU7L213zE1V1h4U1T0zyTZl9a+4PFtZcl+SNSXZrLrooyd2S/Or8m3Kn0gtsrdba5/wsXxtyz9BrPdSSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTQ4bda3yudX0z7feq6qYkf5zkk0kekOSxmQ25vr21dvs2x9baJ6rqGZkNoq6oqkuSHMvs3LX7za9fuli8tXZlVV2c5DlJ3l1VlyW5U5InJ7lnkme21m5YynRxksclOT/JVVX1piT3SnJBZg9IeHpb+PbZ3PcluTLJK6rqm5Ncm+ShSR6Z2fbO5y/l2ncvAAAAAGyutXwzLbOh0Rdl9lTP5yT5uiS/mOQBrbX/vnxza+3yJA9P8tYkT0ryzCSfma+9sK0YfbbWnpvk32T2kIPvSvKdmZ3D9vjW2qtW3H9LZk8UfVGSeyR5dpJvSXJ5kge31q5asea6JA9K8iuZDdGem9kTSV+R5Btbax9bsWbfvQAAAACwmdbyzbTW2kuTvHSfa96e2bfX9rPmNUles4/7b0rywvnP0DV/nuRp+8y1714AAAAA2DxjP4AA4HabcMDmumvJIMNYtWSQYaxaMsgwVi0ZZBirlgwyjFVLhtXvY5gGHKLlXc9V9TnXln8/lWs91JJBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9dNLBo63rjPTAAAAAGDyDNMAAAAAYCDDNAAAAAAYyJlpwKFZtd9++dqQezaplgwyjFVLBhnGqiWDDGPVkkGGsWrJIMNYtWRY/T6GacAh2oQDNtdZSwYZptyPDP1kmFo/MvSTYWr9yNBPhqn1I0M/GabWTy8ZOJ5tngAAAAAwkGEaAAAAAAxkmydwaFZ9RXj52pB7NqmWDDKMVUsGGcaqJYMMY9WSQYaxaskgw1i1ZFj9PoZpwCHahDMB1llLBhmm3I8M/WSYWj8y9JNhav3I0E+GqfUjQz8ZptZPLxk4nm2eAAAAADCQYRoAAAAADGSYBgAAAAADOTMNWLu99tWvur58bcg9m1RLBhnGqiWDDGPVkkGGsWrJIMNYtWSQYaxaMqx+H8M0YCSbesDmOmvJIMOU+5GhnwxT60eGfjJMrR8Z+skwtX5k6CfD1PrpJQPHs80TAAAAAAYyTAMAAACAgQzTAAAAAGAgZ6YBo1i1t37ItYOu67WWDDKMVUsGGcaqJYMMY9WSQYaxaskgw1i1ZFj9PoZpwEg29YDNddaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTSwaOZ5snAAAAAAxkmAYAAAAAA9nmCYxi1deBh1w76Lpea8kgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVUuG1e9jmAaMZFPPBFhnLRlkmHI/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaP71k4Hi2eQIAAADAQIZpAAAAADCQYRoAAAAADOTMNGAUq/bWD7l20HW91pJBhrFqySDDWLVkkGGsWjLIMFYtGWQYq5YMq9/HMA0YyaYesLnOWjLIMOV+ZOgnw9T6kaGfDFPrR4Z+MkytHxn6yTC1fnrJwPFs8wQAAACAgQzTAAAAAGAgwzQAAAAAGMiZacAoVu2tH3LtoOt6rSWDDGPVkkGGsWrJIMNYtWSQYaxaMsgwVi0ZVr+PYRowkk09YHOdtWSQYcr9yNBPhqn1I0M/GabWjwz9ZJhaPzL0k2Fq/fSSgePZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz04BRrNpbP+TaQdf1WksGGcaqJYMMY9WSQYaxaskgw1i1ZJBhrFoyrH4fwzRgJJt6wOY6a8kgw5T7kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+uklA8ezzRMAAAAABjJMAwAAAICBbPMERrHq68BDrh10Xa+1ZJBhrFoyyDBWLRlkGKuWDDKMVUsGGcaqJcPq9zFMA0ayqWcCrLOWDDJMuR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK2fXjJwPNs8AQAAAGAgwzQAAAAAGMgwDQAAAAAGcmYaMIpVe+uHXDvoul5rySDDWLVkkGGsWjLIMFYtGWQYq5YMMoxVS4bV72OYBoxkUw/YXGctGWSYcj8y9JNhav3I0E+GqfUjQz8ZptaPDP1kmFo/vWTgeLZ5AgAAAMBAhmkAAAAAMJBhGgAAAAAM5Mw0YBSr9tYPuXbQdb3WkkGGsWrJIMNYtWSQYaxaMsgwVi0ZZBirlgyr38cwDRjJph6wuc5aMsgw5X5k6CfD1PqRoZ8MU+tHhn4yTK0fGfrJMLV+esnA8WzzBAAAAICBDNMAAAAAYCDbPIFRrPo68JBrB13Xay0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWrJsPp9DNOAkWzqmQDrrCWDDFPuR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT66eXDBzPNk8AAAAAGMgwDQAAAAAGMkwDAAAAgIGcmQaMYtXe+iHXDrqu11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYfX7GKYBI9nUAzbXWUsGGabcjwz9ZJhaPzL0k2Fq/cjQT4ap9SNDPxmm1k8vGTiebZ4AAAAAMJBhGgAAAAAMZJgGAAAAAAM5Mw0Yxaq99UOuHXRdr7VkkGGsWjLIMFYtGWQYq5YMMoxVSwYZxqolw+r3MUwDRrKpB2yus5YMMky5Hxn6yTC1fmToJ8PU+pGhnwxT60eGfjJMrZ9eMnA82zwBAAAAYCDDNAAAAAAYyDANAAAAAAZyZhowilV764dcO+i6XmvJIMNYtWSQYaxaMsgwVi0ZZBirlgwyjFVLhtXvY5gGjGRTD9hcZy0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj+9ZOB4tnkCAAAAwECGaQAAAAAwkG2ewClb9dXfg16bWi0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWrJsPp9DNOANZnKmQDrrCWDDFPuR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT66eXDBzPNk8AAAAAGMgwDQAAAAAGWtswraq+tareWFUfqqqbquoDVfUbVfWNS/fdu6raCX4uOcFnPLWq3llVn6qqj1fVFVX1uBPcf9equqiq3ltVN1fVR6rqtVV1/xOs+YqqenVV/WVV3VJVN1TVy6tq5wRrzq6qN1TVsaq6sareXVU/WFWfd7I/NwAAAAA2x1rOTKuqlyR5XpKPJbk8yUeT/NMkT0zypKr6ztbary8t+6P5vcv+eI/PeFmS5yb5UJJfSHKnJBcm+a2qemZr7VVL9985ye8lOSfJu5L8dJKvTHJBkm+tqnNba1ctrTk9yZVJviTJ65L8aZKHJPmBJI+uqnNaax9bWvPEJL+Z5OYklyY5luTxSX5q/tkXrOoHpmZ5L/2qvfVDrh10Xa+1ZJBhrFoyyDBWLRlkGKuWDDKMVUsGGcaqJcPq91nDMK2qvjTJv0/y10ke2Fr7yMJ7j0zy5iQvSrI8TPvD1tqPDvyMszMbpF2X5MGttb+dX39pkquTvKyqXt9au2Fh2XMyG2ZdluTJrbXb5msuzWyI9+qq+rrd63M/m9kg7VmttVcufP7FSZ6d5MVJvmfh+mmZDfY+m+QRrbV3za+/YN73+VV1YWttz2/bwVRM5YDNddaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTSwaOt45tnl81r3NVWxikJUlr7S1JPpnki0/xM3YHWC9u80HavP4NSX4myZ2TPG33es3+095d87y2MDBrrb0uyduSfG2Shy+suW+S85Ls1lz0wiSfTvKUqrrbwvXzM+vtkjYfpM0/4+YkPzz/9Xv31yoAAAAAvVrHMO39Sf4+yUOq6h8vvlFVD0vyRUl+f8W6f1JV311V/3H++sATfMa589ffXfHe7yzdkySnJ7lXkve11q4fuGb3n9/YPvfbammtfTLJ25N8QZJvGJjrrUluTHJ2zbacAgAAALDhTnmbZ2vtWFX9hyQXJ3lPVV2e2dlppyd5Qmbnln33iqXfMv+5XVVdkeSprbUPLly7W5IvT/Kp1tqHV9R5//z1jIVr95u/vm+P2Addc958zZtOtqa1dmtVXZ/kAUnum+TaPeomSarq6j3e+poTrQMAAADg8KzlAQSttZdX1Q1JXp3kGQtv/VmSX1na/nljkh/L7NyyD8yvPTDJjyZ5ZJI3VdXXt9Y+PX/v7vPXj+/x8bvX77Fwrec1MEnLe+lX7a0fcu2g63qtJYMMY9WSQYaxaskgw1i1ZJBhrFoyyDBWLRlWv8/6nub5vCT/V5JXJHlVkr/K7BtV/3eS/6dmw7HnJcl8sPYjSyXeWlXnJfkfSR6a5N9m9vTN/Wgnv+UfIve2prV21soCs2+snbmPz4QjMZUDNtdZSwYZptyPDP1kmFo/MvSTYWr9yNBPhqn1I0M/GabWTy8ZON4pn5lWVY9I8pIk/6219pzW2gdaaze21q5J8m1J/iLJc2t2wP+eWmu3JvnF+a8PW3jr4/PXu++x9O5L9w1Zc9oRrgEAAABgQ63jAQSPm7++ZfmN1tqNSd45/5x/PqDW38xfb39iZptt9/yLJF9YVV+2Ys1Xz18Xzy177/z1jKw2+pqqumOS+yS5Nf+wnRUAAACADbaOYdrukyq/eI/3d6///YBau0/KXB4+vXn++ugVax6zdE+SXJfkg0nOqKr7DFyzOww8r6o+58+lqr4oyTlJbkryjoG5HpbZ0z+vbK3dsuJ9AAAAADZNa+2UfpL8H5mdCfZXSb586b3HJLktsyHUP5pfe2iSO62oc26Sm+e1zl567+z59T9LsrNw/d6ZPTn05iT3Xlrzf87X/EaSOyxcf+L8+p8sXp+/99/n7z1z6frF8+s/v3T9tMy+TXdLkgctXL9Lkivnay48xT/fq88888wGPZv/Xffjx48fP378+PHjx48fPxP/2dnZOer/E3QtzjzzzJbk6naAWc06HkBwWZLfT/KoJNdW1f+b2WDt/pltAa0kP9Ra+9j8/pckeUBVXZHkQ/NrD8xsmJYkL2itXbn4Aa21K6vq4iTPSfLuqrosyZ2SPDnJPTMbft2wlOvi+eefn+SqqnpTknsluSCzJ4o+vbV229Ka78tsCPaKqvrmJNdmNvx7ZGbbO5+/lOsTVfWM+Z/BFVV1SZJjSZ6Q5H7z65ee8E8PJqJN5IDNddaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTc4Ztd8rDtNbabVX12CTfn+TCzB468AWZDZXekOQVrbU3Liz5tfk9D87sm2ufn+Svk7w2yataa2/b43OeW1XvTvLvknxXZt94uybJS1trr19x/y1V9agkP5TkO5I8O8knklye5IWttfesWHNdVT0oyYsy27r52CQfzuwppRe11o6tWHN5VT08s0HbkzL7VtqfZTb4e0Vb/lsHAAAAwMZaxzfT0lr7TJKXz39Odu8vJfmlA37Oa5K8Zh/335TkhfOfoWv+PMnT9pnr7ZkN3gAAAACYsLUM0wCWv+q76qu/Q64ddF2vtWSQYaxaMsgwVi0ZZBirlgwyjFVLBhnGqiXD8DrbxjANWIupnAmwzloyyDDlfmToJ8PU+pGhnwxT60eGfjJMrR8Z+skwtX56zrDt7nDUAQAAAABgUximAQAAAMBAhmkAAAAAMJAz04C1WN43v2of/ZBrB13Xay0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWrJMLzOtjFMA9ZiKgdsrrOWDDJMuR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK2fnjNsO9s8AQAAAGAgwzQAAAAAGMgwDQAAAAAGcmYasBbL++ZX7aMfcu2g63qtJYMMY9WSQYaxaskgw1i1ZJBhrFoyyDBWLRmG19k2hmnAWkzlgM111pJBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9dNzhm1nmycAAAAADGSYBgAAAAAD2eYJrMXyV31XffV3yLWDruu1lgwyjFVLBhnGqiWDDGPVkkGGsWrJIMNYtWQYXmfbGKYBazGVMwHWWUsGGabcjwz9ZJhaPzL0k2Fq/cjQT4ap9SNDPxmm1k/PGbadbZ4AAAAAMJBhGgAAAAAMZJgGAAAAAAM5Mw1Yi+V986v20Q+5dtB1vdaSQYaxaskgw1i1ZJBhrFoyyDBWLRlkGKuWDMPrbBvDNGAtpnLA5jprySDDlPuRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T66TnDtrPNEwAAAAAGMkwDAAAAgIEM0wAAAABgIGemAWuxvG9+1T76IdcOuq7XWjLIMFYtGWQYq5YMMoxVSwYZxqolgwxj1ZJheJ1tY5gGrMVUDthcZy0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj89Z9h2tnkCAAAAwECGaQAAAAAwkGEaAAAAAAzkzDRgLZb3za/aRz/k2kHX9VpLBhnGqiWDDGPVkkGGsWrJIMNYtWSQYaxaMgyvs20M04C1mMoBm+usJYMMU+5Hhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrp+cM2842TwAAAAAYyDANAAAAAAayzRNYi+Wv+q766u+Qawdd12stGWQYq5YMMoxVSwYZxqolgwxj1ZJBhrFqyTC8zrYxTAPWYipnAqyzlgwyTLkfGfrJMLV+ZOgnw9T6kaGfDFPrR4Z+Mkytn54zbDvbPAEAAABgIMM0AAAAABjIMA0AAAAABnJmGrAWy/vmV+2jH3LtoOt6rSWDDGPVkkGGsWrJIMNYtWSQYaxaMsgwVi0ZhtfZNoZpwFpM5YDNddaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTc4ZtZ5snAAAAAAxkmAYAAAAAAxmmAQAAAMBAzkwD9mWv/fHL11fdN+TaQdf1WksGGcaqJYMMY9WSQYaxaskgw1i1ZJBhrFoyDK+zbQzTgH2b8gGbU+tHhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrR4Z+Mkytn54zbDvbPAEAAABgIMM0AAAAABjINk9g31Z9rXf52pB71rmu11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYXidbWOYBuzblM8EmFo/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaPzL0k2Fq/fScYdvZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz04B9W7VHfvnakHvWua7XWjLIMFYtGWQYq5YMMoxVSwYZxqolgwxj1ZJheJ1tY5gG7NuUD9icWj8y9JNhav3I0E+GqfUjQz8ZptaPDP1kmFo/MvSTYWr99Jxh29nmCQAAAAADGaYBAAAAwECGaQAAAAAwkDPTgH1btUd++dqQe9a5rtdaMsgwVi0ZZBirlgwyjFVLBhnGqiWDDGPVkmF4nW1jmAbs25QP2JxaPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj8y9JNhav30nGHb2eYJAAAAAAMZpgEAAADAQIZpAAAAADCQM9OAfVu1R3752pB71rmu11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYXidbWOYBuzblA/YnFo/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaPzL0k2Fq/fScYdvZ5gkAAAAAAxmmAQAAAMBAtnkC+7bqa73L14bcs851vdaSQYaxaskgw1i1ZJBhrFoyyDBWLRlkGKuWDMPrbBvDNGDfpnwmwNT6kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT66fnDNvONk8AAAAAGMgwDQAAAAAGMkwDAAAAgIGcmQbs26o98svXhtyzznW91pJBhrFqySDDWLVkkGGsWjLIMFYtGWQYq5YMw+tsG8M0YN+mfMDm1PqRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrp+cM2842TwAAAAAYyDANAAAAAAYyTAMAAACAgZyZBuzbqj3yy9eG3LPOdb3WkkGGsWrJIMNYtWSQYaxaMsgwVi0ZZBirlgzD62wbwzRg36Z8wObU+pGhnwxT60eGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+un5wzbzjZPAAAAABjIMA0AAAAABlrbNs+q+tYkP5Dka5P8oyQfTnJ1kotba/9zxf1nJ/nhJN+Q5C5J/izJq5O8srX22T0+46lJvn/+GZ9N8r+SvKy19vo97r9rkh9KcmGSr0ryiSRXJHlha+3aPdZ8RZIXJXn0Qh+XJ7motfa3e6zZdy+wyVZ9rXf52pB71rmu11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYXidbbOWYVpVvSTJ85J8LLPB00eT/NMkT0zypKr6ztbary/c/8Qkv5nk5iSXJjmW5PFJfirJOUkuWPEZL0vy3CQfSvILSe6U2ZDst6rqma21Vy3df+ckvzev964kP53kK+e1v7Wqzm2tXbW05vQkVyb5kiSvS/KnSR6S2ZDw0VV1TmvtY0tr9t0LbLopnwkwtX5k6CfD1PqRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T66TnDtjvlYVpVfWmSf5/kr5M8sLX2kYX3HpnkzZl90+vX59dOy2wY9tkkj2itvWt+/QXze8+vqgtba5cs1Dk7s0HadUke3ObfEKuql2b27beXVdXrW2s3LER7TmbDrMuSPLm1dtt8zaWZDfxeXVVft3t97mczG6Q9q7X2yoXPvzjJs5O8OMn3LFzfdy8AAAAAbK51nJn2VfM6Vy0O0pKktfaWJJ9M8sULl8+f/37J7vBpfu/NmW2VTJLvXfqM3QHWi9vCVsv58Oxnktw5ydN2r9dsTLq75nmLA7PW2uuSvC2zraIPX1hz3yTnJdmtueiFST6d5ClVdbdT7AUAAACADbWOYdr7k/x9kodU1T9efKOqHpbki5L8/sLlc+evv7ui1luT3Jjk7Jpt0xyy5neW7kmS05PcK8n7WmvXD1yz+89vXPq2Wlprn0zy9iRfkNm5aENy7dXLSlV19aqfJF9zsrUAAAAAHI5T3ubZWjtWVf8hycVJ3lNVl2d2dtrpSZ6Q2bll372w5H7z1/etqHVrVV2f5AFJ7pvk2vk3wb48yadaax9eEeH989czhnzGKa45b77mTSdbs6qXPerCxlm1R3752pB71rmu11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYXidbbOWBxC01l5eVTdk9gTLZyy89WdJfmVp++fd568f36Pc7vV7HPD+3tes1Fo7a9X1mn077cyTrYfDNOUDNqfWjwz9ZJhaPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj89Z9h269jmmap6XmYH/f9KZt9Iu1uSs5J8IMn/U1U/sZ9y89d2wruOt5/7D/IZh7UGAAAAgE6d8jCtqh6R5CVJ/ltr7TmttQ+01m5srV2T5NuS/EWS59bsgP/kH76tdfc9Sp62dN/J7r/70n0H+YzDXAMAAADAhlrHN9MeN399y/IbrbUbk7xz/jn/fH75vfPXM5bvr6o7JrlPklsz+1ZbWmufzmwg94VV9WUrPv+r56+L55bt+RmHtWZVLwAAAABstnWcmbb7pMov3uP93et/P399c5J/leTRSf7r0r0Py+yJmW9trd2ycP3NSZ4yX/PLS2ses3DPruuSfDDJGVV1n3b8Ez1XrdkdBp5XVXdoC0/0rKovSnJOkpuSvGMp1357gY23ao/88rUh96xzXa+1ZJBhrFoyyDBWLRlkGKuWDDKMVUsGGcaqJcPwOttmHd9Me9v89buq6ssX36iqx2Q2hLo5yZXzy5cl+WiSC6vqQQv33iXJj89//bmlz/j5+evzq2pnYc29k3x/kluyMGRrs9Pxdtf8RFXdYWHNE5N8U5L3JPmDhTXXJXljkt2aiy7K7By4X51/U27XQXqBjdda+5yf5WtD7lnnul5rySDDlPuRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T66TXD7vVtto5vpl2W5PeTPCrJtVX1/yb5qyT3z2wLaCX5odbax5KktfaJqnrGfN0VVXVJkmNJnpDkfvPrly5+QGvtyqq6OMlzkry7qi5LcqckT05yzyTPbK3dsJTr4vnnn5/kqqp6U5J7JbkgyY1Jnt4Wvn02932ZDf1eUVXfnOTaJA9N8sjMtnc+fynXvnsBAAAAYHOd8jfT5gOpxyZ5dmbf9vq2JM9N8g1J3pDkf2+t/fTSmsuTPDzJW5M8Kckzk3wms2HZhW3FmLO19twk/yazQd13JfnOJH+S5PGttVetuP+WzAZ8L0pyj3m+b0lyeZIHt9auWrHmuiQPyuyppA+d93F6klck+cbdgeCp9gIAAADAZlrHN9PSWvtMkpfPf4aueXtmQ7j9fM5rkrxmH/fflOSF85+ha/48ydP2mWvfvQAAAACwedYyTAO2y9QP2FxnLRlkGKuWDDKMVUsGGcaqJYMMY9WSQYaxaskwvM62MUwD9m1593JVfc615d+HXjvoul5rySDDlPuRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T66TnDtlvH0zwBAAAAYCsYpgEAAADAQLZ5Avu26mu9y9eG3LPOdb3WkkGGsWrJIMNYtWSQYaxaMsgwVi0ZZBirlgzD62wbwzRg36Z8JsDU+pGhnwxT60eGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+un5wzbzjZPAAAAABjIMA0AAAAABjJMAwAAAICBnJkG7NuqPfLL14bcs851vdaSQYaxaskgw1i1ZJBhrFoyyDBWLRlkGKuWDMPrbBvDNGDfpnzA5tT6kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT66fnDNvONk8AAAAAGMgwDQAAAAAGMkwDAAAAgIGcmQbs26o98svXhtyzznW91pJBhrFqySDDWLVkkGGsWjLIMFYtGWQYq5YMw+tsG8M0YN+mfMDm1PqRoZ8MU+tHhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrp+cM2842TwAAAAAYyDANAAAAAAayzRPY015f3111ffnakHvWua7XWjLIMFYtGWQYq5YMMoxVSwYZxqolgwxj1ZJheJ1tY5gGnFAP+/E3tZYMMky5Hxn6yTC1fmToJ8PU+pGhnwxT60eGfjJMrZ+eM2w72zwBAAAAYCDDNAAAAAAYyDANAAAAAAZyZhpwQqv2ww+5dtjreq0lgwxj1ZJBhrFqySDDWLVkkGGsWjLIMFYtGYbX2TaGacAJ9XC45abWkkGGKfcjQz8ZptaPDP1kmFo/MvSTYWr9yNBPhqn103OGbWebJwAAAAAMZJgGAAAAAAMZpgEAAADAQM5MA05o1X74IdcOe12vtWSQYaxaMsgwVi0ZZBirlgwyjFVLBhnGqiXD8DrbxjANOKEeDrfc1FoyyDDlfmToJ8PU+pGhnwxT60eGfjJMrR8Z+skwtX56zrDtbPMEAAAAgIEM0wAAAABgIMM0AAAAABjImWnACa3aDz/k2mGv67WWDDKMVUsGGcaqJYMMY9WSQYaxaskgw1i1ZBheZ9sYpgEn1MPhlptaSwYZptyPDP1kmFo/MvSTYWr9yNBPhqn1I0M/GabWT88Ztp1tngAAAAAwkGEaAAAAAAxkmydwQqu+wjvk2mGv67WWDDKMVUsGGcaqJYMMY9WSQYaxaskgw1i1ZBheZ9sYpgEn1MN+/E2tJYMMU+5Hhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrp+cM2842TwAAAAAYyDANAAAAAAYyTAMAAACAgZyZBpzQqv3wQ64d9rpea8kgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVUuG4XW2jWEacEI9HG65qbVkkGHK/cjQT4ap9SNDPxmm1o8M/WSYWj8y9JNhav30nGHb2eYJAAAAAAMZpgEAAADAQIZpAAAAADCQM9OAE1q1H37ItcNe12stGWQYq5YMMoxVSwYZxqolgwxj1ZJBhrFqyTC8zrYxTANOqIfDLTe1lgwyTLkfGfrJMLV+ZOgnw9T6kaGfDFPrR4Z+Mkytn54zbDvbPAEAAABgIMM0AAAAABjIMA0AAAAABnJmGnBCq/bDD7l22Ot6rSWDDGPVkkGGsWrJIMNYtWSQYaxaMsgwVi0ZhtfZNoZpwAn1cLjlptaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTc4ZtZ5snAAAAAAxkmAYAAAAAA9nmCZzQqq/wDrl22Ot6rSWDDGPVkkGGsWrJIMNYtWSQYaxaMsgwVi0ZhtfZNoZpwAn1sB9/U2vJIMOU+5GhnwxT60eGfjJMrR8Z+skwtX5k6CfD1PrpOcO2s80TAAAAAAYyTAMAAACAgQzTAAAAAGAgZ6YBJ7RqP/yQa4e9rtdaMsgwVi0ZZBirlgwyjFVLBhnGqiWDDGPVkmF4nW1jmAacUA+HW25qLRlkmHI/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaPz1n2Ha2eQIAAADAQIZpAAAAADCQYRoAAAAADOTMNOCEVu2HH3LtsNf1WksGGcaqJYMMY9WSQYaxaskgw1i1ZJBhrFoyDK+zbQzTgBPq4XDLTa0lgwxT7keGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+un5wzbzjZPAAAAABjIMA0AAAAABrLNEzihVV/hHXLtsNf1WksGGcaqJYMMY9WSQYaxaskgw1i1ZJBhrFoyDK+zbQzTgBPqYT/+ptaSQYYp9yNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTc4ZtZ5snAAAAAAxkmAYAAAAAAxmmAQAAAMBAzkwDkuy9733V9SHXDntdr7VkkGGsWjLIMFYtGWQYq5YMMoxVSwYZxqolw/A628YwDbhdr4dbbmotGWSYcj8y9JNhav3I0E+GqfUjQz8ZptaPDP1kmFo/PWfYdrZ5AgAAAMBAhmkAAAAAMNApD9Oq6t9UVTvJz2cX7r/3Se695ASf9dSqemdVfaqqPl5VV1TV405w/12r6qKqem9V3VxVH6mq11bV/U+w5iuq6tVV9ZdVdUtV3VBVL6+qnROsObuq3lBVx6rqxqp6d1X9YFV93pA/QwAAAAA2wzrOTPvDJBft8d43JTk3ye+seO+Pkly+4vofrypUVS9L8twkH0ryC0nulOTCJL9VVc9srb1q6f47J/m9JOckeVeSn07ylUkuSPKtVXVua+2qpTWnJ7kyyZckeV2SP03ykCQ/kOTRVXVOa+1jS2uemOQ3k9yc5NIkx5I8PslPzT/7glX9QI9W7X0/6LXDXtdrLRlkGKuWDDKMVUsGGcaqJYMMY9WSQYaxaskwvM62OeVhWmvtDzMbqB2nqv7n/B//y4q3/7C19qNDPqOqzs5skHZdkge31v52fv2lSa5O8rKqen1r7YaFZc/JbJh1WZInt9Zum6+5NLMh3qur6ut2r8/9bGaDtGe11l658PkXJ3l2khcn+Z6F66dlNtj7bJJHtNbeNb/+giRvTnJ+VV3YWtvz23bQk14Pt9zUWjLIMOV+ZOgnw9T6kaGfDFPrR4Z+MkytHxn6yTC1fnrOsO1GOzOtqv5Zkm9I8hdJfvsUy+0OsF7c5oO0JGmz4dnPJLlzkqctfHYtrHleWxiYtdZel+RtSb42ycMX1tw3yXlJdmsuemGSTyd5SlXdbeH6+Um+OMklbT5Im3/GzUl+eP7r9+6vVQAAAAB6NeYDCL57/vpLrbXPrnj/n1TVd1fVf5y/PvAEtc6dv/7uivd+Z+meJDk9yb2SvK+1dv3ANbv//Mb2ud9WS2vtk0nenuQLMhsQDsn11iQ3Jjm7ZltOAQAAANhw6zgz7ThVddck/zrJbUl+cY/bvmX+s7juiiRPba19cOHa3ZJ8eZJPtdY+vKLO++evZyxcu9/89X17fPZB15w3X/Omk61prd1aVdcneUCS+ya5do+6SZKqunqPt77mROsAAAAAODyjDNOS/B9J7pHkt1trf7703o1Jfiyzc8s+ML/2wCQ/muSRSd5UVV/fWvv0/L27z18/vsdn7V6/x8K1ntdAt1btfT/otcNe12stGWQYq5YMMoxVSwYZxqolgwxj1ZJBhrFqyTC8zrYZa5j2XfPX/7z8RmvtI0l+ZOnyW6vqvCT/I8lDk/zbzJ6+uR/t5Lfcbvc/+W7WtNbOWllg9o21M/fxmXBgvR5uuam1ZJBhyv3I0E+GqfUjQz8ZptaPDP1kmFo/MvSTYWr99Jxh2639zLSq+tokZyf5UJI3DF3XWrs1/7Al9GELb318/nr3PZbefem+IWtOO8I1AAAAAGyoMR5AcLIHD5zI38xfb39iZptt9/yLJF9YVV+2Ys1Xz18Xzy177/z1jKw2+pqqumOS+yS5Nf+wnRUAAACADbbWbZ5VdZckT8nswQO/dIASu0/KXB4+vXle99FJfnnpvccs3LPruiQfTHJGVd2nHf9Ez1Vr3jJ/Pa+q7tAWnuhZVV+U5JwkNyV5x1KufzXP9V+XPuNhmT39862ttVsCG2DV13UPeu2w1/VaSwYZxqolgwxj1ZJBhrFqySDDWLVkkGGsWjIMr7Nt1n1m2gVJdpK8vh3/4IEkSVU9NMn/aq39/dL1c5M8e/7rry8t+/nMhmnPr6rLW2t/O19z7yTfn+SWLAzZWmutqn4+yf+V5Ceq6sm7w7GqemKSb0ryniR/sLDmuqp6Y2ZP7Pz+JK9c+PyLMvu23H9u//BghCS5LMlLklxYVa9srb1r/hl3SfLj83t+btWfA/So1/34m1pLBhmm3I8M/WSYWj8y9JNhav3I0E+GqfUjQz8ZptZPzxm23bqHabsPHvgvJ7jnJUkeUFVXZHauWjJ7mue5839+QWvtysUFrbUrq+riJM9J8u6quizJnZI8Ock9kzyztXbD0udcnORxSc5PclVVvSnJvTIb+N2Y5Olt4dtnc9+X5Mokr6iqb05ybWYPRHhkZts7n7+U6xNV9YzMhmpXVNUlSY4leUKS+82vX3qCPwsAAAAANsjahmlVdf8k/yInf/DAryX5tiQPzmy75ecn+eskr03yqtba21Ytaq09t6reneTfZTa0uy3JNUle2lp7/Yr7b6mqRyX5oSTfkdm33j6R5PIkL2ytvWfFmuuq6kFJXpTZ1s3HJvlwklckuai1dmzFmsur6uGZDdqelOQuSf4ss8HfK9ryCBcAAACAjbW2YVpr7dokJ/2uX2vtl3Kw89TSWntNktfs4/6bkrxw/jN0zZ8nedo+c709s8EbAAAAABO27m2ewAZbtff9oNcOe12vtWSQYaxaMsgwVi0ZZBirlgwyjFVLBhnGqiXD8DrbxjANuF2vh1tuai0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj89Z9h2dzjqAAAAAACwKQzTAAAAAGAgwzQAAAAAGMiZacDtVu19P+i1w17Xay0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWrJMLzOtjFMA27X6+GWm1pLBhmm3I8M/WSYWj8y9JNhav3I0E+GqfUjQz8ZptZPzxm2nW2eAAAAADCQYRoAAAAADGSbJ3C7VV/XPei1w17Xay0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWrJMLzOtjFMA27X6378Ta0lgwxT7keGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+un5wzbzjZPAAAAABjIMA0AAAAABjJMAwAAAICBnJkG3G7V3veDXjvsdb3WkkGGsWrJIMNYtWSQYaxaMsgwVi0ZZBirlgzD62wbwzTgdr0ebrmptWSQYcr9yNBPhqn1I0M/GabWjwz9ZJhaPzL0k2Fq/fScYdvZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz04Dbrdr7ftBrh72u11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WSYXidbWOYBtyu18MtN7WWDDJMuR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK2fnjNsO9s8AQAAAGAgwzQAAAAAGMgwDQAAAAAGcmYacLtVe98Peu2w1/VaSwYZxqolgwxj1ZJBhrFqySDDWLVkkGGsWjIMr7NtDNOA2/V6uOWm1pJBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9dNzhm1nmycAAAAADGSYBgAAAAAD2eYJ3G7V13UPeu2w1/VaSwYZxqolgwxj1ZJBhrFqySDDWLVkkGGsWjIMr7NtDNOA2/W6H39Ta8kgw5T7kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+uk5w7azzRMAAAAABjJMAwAAAICBDNMAAAAAYCBnpgG3W7X3/aDXDntdr7VkkGGsWjLIMFYtGWQYq5YMMoxVSwYZxqolw/A628YwDbhdr4dbbmotGWSYcj8y9JNhav3I0E+GqfUjQz8ZptaPDP1kmFo/PWfYdrZ5AgAAAMBAhmkAAAAAMJBhGgAAAAAM5Mw04Har9r4f9Nphr+u1lgwyjFVLBhnGqiWDDGPVkkGGsWrJIMNYtWQYXmfbGKYBt+v1cMtNrSWDDFPuR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT66fnDNvONk8AAAAAGMgwDQAAAAAGss0TuN2qr+se9Nphr+u1lgwyjFVLBhnGqiWDDGPVkkGGsWrJIMNYtWQYXmfbGKYBt+t1P/6m1pJBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9dNzhm1nmycAAAAADGSYBgAAAAADGaYBAAAAwEDOTANut2rv+0GvHfa6XmvJIMNYtWSQYaxaMsgwVi0ZZBirlgwyjFVLhuF1to1hGnC7Xg+33NRaMsgw5X5k6CfD1PqRoZ8MU+tHhn4yTK0fGfrJMLV+es6w7WzzBAAAAICBDNMAAAAAYCDDNAAAAAAYyJlpsKVW7XNf57XDXtdrLRlkGKuWDDKMVUsGGcaqJYMMY9WSQYaxaskwvM62MUyDLbYph1tuai0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj89Z9h2tnkCAAAAwECGaQAAAAAwkGEaAAAAAAzkzDTYYst73VftfT/otcNe12stGWQYq5YMMoxVSwYZxqolgwxj1ZJBhrFqyTC8zrYxTIMttimHW25qLRlkmHI/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaPz1n2Ha2eQIAAADAQIZpAAAAADCQbZ6wxZa/nrvq67oHvXbY63qtJYMMY9WSQYaxaskgw1i1ZJBhrFoyyDBWralnqKqTXlt1z87OznGftW0M02CLbcp+/E2tJYMMU+5Hhn4yTK0fGfrJMLV+ZOgnw9T6kaGfDFPrZ+wMHJxtngAAAAAwkGEaAAAAAAxkmAYAAAAAAzkzDbZY1XYdsHkUtWSQYaxaMsgwVi0ZZBirlgwyjFVLBhnGqjWlDB4asF6GabDFtu2AzW3sR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT60eGfjJMrZ91ZmC9bPMEAAAAgIEM0wAAAABgIMM0AAAAABjImWmwxaqme8BmL7VkkGGsWjLIMFYtGWQYq5YMMoxVSwYZxqq1yRlW3cP6GKbBFpvyAZs91JJBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9XOqGRiPbZ4AAAAAMJBhGgAAAAAMZJsnbLHlr/6u+irwQa8d9rpea8kgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVWtKGXZ2dnLs2LHj7uFgDNNgi03tTIDeaskgw5T7kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+ll3BtbHNk8AAAAAGMgwDQAAAAAGOuVhWlX9m6pqJ/n57Ip1Z1fVG6rqWFXdWFXvrqofrKrPO8FnPbWq3llVn6qqj1fVFVX1uBPcf9equqiq3ltVN1fVR6rqtVV1/xOs+YqqenVV/WVV3VJVN1TVy6tq5wRr9t0LAAAAAJunlvfV7rtA1dcn+Zd7vP1NSc5N8tuttcctrHlikt9McnOSS5McS/L4JPdLcllr7YIVn/OyJM9N8qEklyW5U5ILk9wzyTNba69auv/OSd6U5Jwk70ry5iRfmeSCJH+f5NzW2lVLa05PcmWSL0nyuiR/muQhSR6Z5L1JzmmtfWxpzb572Y+quvrMM8888+qrrz6VMnAce+YBAAC2Q1XltttuO+oYXTnrrLNyzTXXXNNaO2u/a0/5AQSttT9M8oer3quq/zn/x/+ycO20JL+Q5LNJHtFae9f8+gsyG3idX1UXttYuWVhzdmaDtOuSPLi19rfz6y9NcnWSl1XV61trNyx8/HMyG6RdluTJrbXb5msuTXJ5kldX1dftXp/72cwGac9qrb1y4fMvTvLsJC9O8j2n0gv0ZOoHbB51LRlkmHI/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaP+vOwPqMdmZaVf2zJN+Q5C+S/PbCW+cn+eIkl7T58ClJWms3J/nh+a/fu1Rud4D14jYfpM3X3JDkZ5LcOcnTFj67FtY8ry0MzFprr0vytiRfm+ThC2vum+S8JLs1F70wyaeTPKWq7naKvQAAAACwocZ8AMF3z19/qbW2eGbaufPX312x5q1Jbkxyds22aQ5Z8ztL9yTJ6UnuleR9rbXrB67Z/ec3ts/9tlpaa59M8vYkX5DZgHBIrr16AQAAAGBDjTJMq6q7JvnXSW5L8otLb99v/vq+5XWttVuTXJ/Z9tP7zmvdLcmXJ/lUa+3DKz7u/fPXM4Z8xmGtWdXLiVTV1at+knzNydYCAAAAcDhO+cy0PfwfSe6R2YMH/nzpvbvPXz++x9rd6/c44P29r4FuLO+bX7WP/qDXDntdr7VkkGGsWjLIMFYtGWQYq5YMMoxVSwYZxqo1pQyr3ufgxhqmfdf89T8fYO3uf8LthHcdbz/3H+QzRl3T9nh6RM2+nXbmPj4TBpv6AZtHXUsGGabcjwz9ZJhaPzL0k2Fq/cjQT4ap9SNDPxmm1s+6M7A+a9/mWVVfm+TsJB9K8oYVt3x8/nr3PUqctnTfye6/+9J9B/mMw1wDAAAAwIYa48y0vR48sOu989czlt+oqjsmuU+SW5N8IElaa5/O7ImgX1hVX7ai3lfPXxfPLdvzMw5rzapeAAAAANhsax2mVdVdkjwlswcP/NIet715/vroFe89LLMnZl7ZWrtl4JrHLN2TJNcl+WCSM6rqPgPXvGX+el5Vfc6fS1V9UZJzktyU5B0Dc+3VCwAAAACbqrW2tp/MBmktyW+d4J7TkvxNkluSPGjh+l2SXDlff+HSmrPn1/8syc7C9Xsn+ViSm5Pce2nN/zlf8xtJ7rBw/Ynz63+yeH3+3n+fv/fMpesXz6///Kn2coA/06vPPPPMBus2//vpx48fP378+PHjx48fP34m/lNVR/1/gnbnzDPPbEmubgeY1az7AQS7Dx74L3vd0Fr7RFU9I8llSa6oqkuSHEvyhCT3m1+/dGnNlVV1cZLnJHl3VV2W5E5JnpzknpkNv25Y+qiLkzwuyflJrqqqNyW5V5ILktyY5OmttduW1nxfZkOwV1TVNye5NslDkzwys+2dzz/VXqAnbeIHbB51LRlkmHI/MvSTYWr9yNBPhqn1I0M/GabWjwz9ZJhaP+vOwPqsbZtnVd0/yb/I3g8euF1r7fIkD0/y1iRPSvLMJJ/JbFh2YVv+WzBb89wk/ybJX2U2tPvOzL5d9vjW2qtW3H9LkkcleVGSeyR5dpJvSXJ5kge31q5asea6JA9K8iuZDdGem+T0JK9I8o2ttY+toxcAAAAANtPavpnWWrs2yeBRZ2vt7Ukeu8/PeE2S1+zj/puSvHD+M3TNnyd52j5z7bsXAAAAADbPurd5Ahtk+au+q776e9Brh72u11oyyDBWLRlkGKuWDDKMVUsGGcaqJYMMY9WaUoZV73NwhmmwxaZ+JsBR15JBhin3I0M/GabWjwz9ZJhaPzL0k2Fq/cjQT4ap9bPuDKzP2s5MAwAAAICpM0wDAAAAgIEM0wAAAABgIGemwRZb3je/ah/9Qa8d9rpea8kgw1i1ZJBhrFoyyDBWLRlkGKuWDDKMVWtKGVa9z8EZpsEWm/oBm0ddSwYZptyPDP1kmFo/MvSTYWr9yNBPhqn1I0M/GabWz7ozsD62eQIAAADAQIZpAAAAADCQYRoAAAAADOTMNNhiy/vmV+2jP+i1w17Xay0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWpNKcOq9zk4wzTYYlM/YPOoa8kgw5T7kaGfDFPrR4Z+MkytHxn6yTC1fmToJ8PU+ll3BtbHNk8AAAAAGMgwDQAAAAAGMkwDAAAAgIGcmQZbbHnf/Kp99Ae9dtjreq0lgwxj1ZJBhrFqySDDWLVkkGGsWjLIMFatKWVY9T4HZ5gGW2zqB2wedS0ZZJhyPzL0k2Fq/cjQT4ap9SNDPxmm1o8M/WSYWj/rzsD62OYJAAAAAAMZpgEAAADAQLZ5whbY6yu9y9dX3XfQa4e9rtdaMsgwVi0ZZBirlgwyjFVLBhnGqiWDDGPVmlKGVe9zcIZpsCW28UyAo64lgwxT7keGfjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+tn3RlYH9s8AQAAAGAgwzQAAAAAGMgwDQAAAAAGcmYabIlVe+SXrw25Z+xaPWRYZy0ZZBirlgwyjFVLBhnGqiWDDGPVkkGGsWpNKcPOzs5x73NwhmmwJbbxgM2jriWDDFPuR4Z+MkytHxn6yTC1fmToJ8PU+pGhnwxT62edGVgv2zwBAAAAYCDDNAAAAAAYyDANAAAAAAZyZhpsiartOmCzl1oyyDBWLRlkGKuWDDKMVUsGGcaqJYMM+61VVUeeYex1HjgwLsM02BLbdsBmD7VkkGHK/cjQT4ap9SNDPxmm1o8M/WSYWj8y9JNhaK1VeuznVDIwLts8AQAAAGAgwzQAAAAAGMg2T9gSVdM+E6DXWjLIMFYtGWQYq5YMMoxVSwYZxqolgwwnurbq/VV2dnaOW9dDP0PucT7a4TNMgy1x1OcSjLmu11oyyDDlfmToJ8PU+pGhnwxT60eGfjJMrR8Z+smwV60hjh07dtJ7ev0z5fDZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz02BLrDorYPnakHvGrtVDhnXWkkGGsWrJIMNYtWSQYaxaMsgwVi0ZZDjRtXUezr/8kILdzxqrn1W1PWygD4ZpsCU24cDQHjJMrR8Z+skwtX5k6CfD1PqRoZ8MU+tHhn4yTK0fGfrJsNe1dRnykIJVGXrIznrZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz02BLVG3GgaE9ZFhnLRlkGKuWDDKMVUsGGcaqJYMMY9WSYe/fD3LtoOtOJfvytVPJ0MMB/Xs9qGBZj9kZxjANtsQmHBjaQ4ap9SNDPxmm1o8M/WSYWj8y9JNhav3I0E+GqfXTa4bDts4/h0039EEFbC7bPAEAAABgIMM0AAAAABjIMA0AAAAABnJmGmyJqs08OHXTa8kgw1i1ZJBhrFoyyDBWLRlkGKvWNmao+tyD93s4uH6vQ/eHXOshP+yHYRpsiU09OHWTa8kgw5T7kaGfDFPrR4Z+MkytHxn6yTC1fo4iQ48cus82sc0TAAAAAAYyTAMAAACAgWzzhC1R1d9ZD71mWGctGWQYq5YMMoxVSwYZxqolgwxj1TqVdYvXln8/lQyrag25tuoe54lBfwzTYEv0eNZDjxmm1o8M/WSYWj8y9JNhav3I0E+GqfUjQz8ZeulniINmAKbNNk8AAAAAGMgwDQAAAAAGMkwDAAAAgIGcmQZbYGdnJ1XTPkC211oyyDBWLRlkGKuWDDKMVUsGGcaqddB7hlj179FVn/uQAA8IgO1jmAZb4NixY8ddq5rWAbI91pJBhin3I0M/GabWjwz9ZJhaPzL0k+Go+jmIVf8eDWCbJwAAAAAMZJgGAAAAAAMZpgEAAADAQM5Mgy22fHbEqrMkDnrtsNf1WksGGcaqJYMMY9WSQYaxaskgw/Lvi9eWfx96beg6DwkA1skwDbbYJh8guwm1ZJBhyv3I0E+GqfUjQz8ZptaPDP1kANh0tnkCAAAAwECGaQAAAAAwkG2esKV2dnZS1ef5GUedYZ21ZJBhrFoyyDBWLRlkGKuWDH1kqDrY2WSrrg1d57wyYGoM02BLHTt27HN+rzr68zN6yDC1fmToJ8PU+pGhnwxT60eGfjJMrR8ZAFgX2zwBAAAAYCDDNAAAAAAYyDANAAAAAAZyZhqQZPUDCZIc+Nphr+u1lgwyjFVLBhnGqiWDDGPVOuoMVes7eP+ws59KLYf/A6yfYRqQ5PgHEiSzfxnbxkN5p9SPDP1kmFo/MvSTYWr9yNBPhin2sy6b9OcAwPrZ5gkAAAAAAxmmAQAAAMBAhmkAAAAAMJAz04ATqtquw4nXXUsGGcaqJYMMY9WSYb0Zqg52wP1B1+2Va0p/pqdyzzqtK/s6/3Mdu2cAZgzTgBPaxsOJp9SPDP1kmFo/MvSTYWr9TC3DUZj6n+l+1o1pXdmHGvrnAMD4bPMEAAAAgIEM0wAAAABgoLUO06rqm6rqN6vqw1V1y/z1jVX12IV77l1V7QQ/l5yg/lOr6p1V9amq+nhVXVFVjzvB/Xetqouq6r1VdXNVfaSqXltV9z/Bmq+oqldX1V/Oe7ihql5eVTsnWHN2Vb2hqo5V1Y1V9e6q+sGq+rwhf24AAAAAbIa1nZlWVT+c5MeSfDTJ65N8OMk/TvLPkzwiyRuWlvxRkstXlPrjPeq/LMlzk3woyS8kuVOSC5P8VlU9s7X2qqX775zk95Kck+RdSX46yVcmuSDJt1bVua21q5bWnJ7kyiRfkuR1Sf40yUOS/ECSR1fVOa21jy2teWKS30xyc5JLkxxL8vgkPzX/7AtW9QObYGdnZ+XZG0OuHfa6XmvJIMNYtWSQYaxaMqx/3WFb9b/fm/JnWrW+hzPs7Oz5/xd+yvb6Mz7ZtVPJtPyZh90zAP9gLcO0qrogs0Ha7yf59tbaJ5fe//wVy/6wtfajA+ufndkg7bokD26t/e38+kuTXJ3kZVX1+tbaDQvLnpPZMOuyJE9urd02X3NpZkO8V1fV1+1en/vZzAZpz2qtvXLh8y9O8uwkL07yPQvXT8tssPfZJI9orb1rfv0FSd6c5PyqurC1tue37aBnx44dO+5a7XHY7SYfTjzlfmToJ8PU+pGhnwxT62cKGY7a8v9+b9Kf6aZY9e9IU/xMAFY75W2eVXWHJC9JcmOS71gepCVJa+0zp/gxuwOsF+8O0uZ1b0jyM0nunORpC5lqYc3zFgdmrbXXJXlbkq9N8vCFNfdNcl6S3ZqLXpjk00meUlV3W7h+fpIvTnLJ7iBt/hk3J/nh+a/fu79WAQAAAOjVOr6ZdnaS+2T2DbC/rapvTfLPMtv2+M7W2v/cY90/qarvTvKPknwsyf9srb17j3vPnb/+7or3fifJC+b3vHB+7fQk90ryvtba9Xus+ab5mrcsfcYbl76tltbaJ6vq7ZkN274hyZsG5HprZgPGs6vqzq21W/boLUlSVVfv8dbXnGgdAAAAAIdnHcO0B89f/zrJNUm+bvHNqnprkvNba3+ztO5b5j+L916R5KmttQ8uXLtbki9P8qnW2odXfP77569nLFy73/z1fXtkPuia8+Zrdodpe65prd1aVdcneUCS+ya5do+6sFGGnqM25J51ruu1lgwyjFVLBhnGqiXD+tcdtVM5A3XVtXWt6/HPCgCGWMcw7Uvmr9+T5Pokj0pyVZKvSvKTSf73JL+R2UMIktm3tX4ss3PLPjC/9sAkP5rkkUneVFVf31r79Py9u89fP77H5+9ev8fCtZ7XrNRaO2vV9fk31s482Xo4LEPOURt6Lsq61vVaSwYZptyPDP1kmFo/U8jQm4Oegbrq2rrXAcAmOuUz05J83vy1MvsG2ptaa59qrf1Jkm/L7OmbD6+qb0yS1tpHWms/0lq7prX2d/Oft2b2ra+rkvzTJP/2ADn2c3rp7v9y97gGAAAAgE6tY5i2+0CAD7TW/mjxjdbaTUn++/zXh5yoSGvt1iS/OP/1YQtv7X676+5ZbdW3w0625rQjXAMAAADAhlrHMO2989e/2+P93WHbXQfU2j1X7fYnZs63e/5Fki+sqi9bsear56+L55btZjojq42+pqrumNmDGW7NP2xnBQAAAGCDrePMtLdmNjD66qq6U2vt75fe/2fz1xsG1PqG+evy8OnNSZ6S5NFJfnnpvccs3LPruiQfTHJGVd1nxRM9V63ZfarneVV1h8UnelbVFyU5J8lNSd6xlOtfzXP916XPeFiSL0jy1pM9yROmYMyDmze5lgwyjFVLBhmWf1+8tvz70GsHXXcq/awzw17XTpbhoNn3urYJVj2U4KB/R3avn+za8u87Ozv7zg0APTjlb6a11j6a5NLMtjr+yOJ7VfUtmT2A4ONJfnd+7aFVdaflOlV1bpJnz3/99aW3f37++vyq2llYc+8k35/kliwM2drsdNPdNT9RVXdYWPPEJN+U5D1J/mBhzXVJ3phkt+aiizL7ttyvLjwYIUkuS/LRJBdW1YMWPuMuSX58/uvPLfcKU9Rau/1n+feh1w66rtdaMsgw5X5k6CdDay233XbbCX8feu2g606ln3VmWL52VP+5boJjx46t7e/IQX9WPRgBADbBOr6ZliTPSfLQzIZdD0vyzsye5vltST6b5Bmttb+b3/uSJA+oqisyezhBMnua57nzf35Ba+3KxeKttSur6uL557y7qi5LcqckT05yzyTPbK3dsJTp4iSPS3J+kquq6k1J7pXkgsyeKPr0tvDts7nvS3JlkldU1TcnuXbe1yMz2975/KVcn6iqZ2Q2VLuiqi5JcizJE5Lcb3790hP+yQEAAACwMdYyTGutfaSqHprkhzMboH1Dkk8m+e0k/3drbXFr5K/N73lwZtstPz/JXyd5bZJXtdbetsdnPLeq3p3k3yX5riS3JbkmyUtba69fcf8tVfWoJD+U5Dsy+9bbJ5JcnuSFrbX3rFhz3fwbZi/KbOvmY5N8OMkrklzUWjvu//ustXZ5VT08s0Hbk5LcJcmfZTb4e0XbpP8vSgAAAABOaF3fTMt80PSc+c+J7vulJL90wM94TZLX7OP+m5K8cP4zdM2fJ3naPnO9PbPBGwAAAAATtrZhGrC99jrEeNmmHOS9zloyyDBWLRn6yLDOw/KHruv10Pbl/y3oJXsPfycBgGkxTANO2fIBwlV13CHMQ64ddF2vtWSQYcr9yMCyXg+TP4q/WwDAtJ3y0zwBAAAAYFsYpgEAAADAQLZ5AqPY1POP1l1LBhnGqrWNGapWnyl2mBl6Pa+MvfXw3wsAYFoM04BRbPL5R1M6z0mGfjJMrZ+jyLCshwz0z5lpAMC62eYJAAAAAAMZpgEAAADAQIZpAAAAADCQM9OAUWzKgeZj15JBhrFqTT3DkIP+d3Z2VtY5zAz0r4f/XgAA02KYBoxiUw40n9oB7TL0m2Fq/YydYYhjx47tew3bxwMIAIB1s80TAAAAAAYyTAMAAACAgQzTAAAAAGAgZ6YB+zL0LJgeDzQ/iloyyDBWrSllcNA/Y+rhvxcAwLQYpgH7tqkHmk/tgHYZNivD1PpZZwYYkwcQAADrZpsnAAAAAAxkmAYAAAAAAxmmAQAAAMBAzkwDkuzvjJdNONC8l1oyyDBWrU3O4EwpDpO/3wDAuhmmAbeb0oHmPdSSQYYp9+OAdjaFv98AwLrZ5gkAAAAAAxmmAQAAAMBAtnnCllrneUurrm3yOVDrrCWDDGPVOuwMVfU515Z/38+1nZ2d4z4PxrCzs3Ogv6f+fgMAJ2KYBltsamcw9VZLBhmm1A9somPHjh11BABggmzzBAAAAICBDNMAAAAAYCDDNAAAAAAYyJlpMDGrDhIfeu82Hqo+di0ZZBir1pgZ9vP/jgAAwLYxTIMJ2pQDzXvMMLV+ZOgnw6b1AwAArGabJwAAAAAMZJgGAAAAAAMZpgEAAADAQM5Mgw0y9ByjHg8036QM66wlgwxj1Rp6z8murbpnZ2fnuNoAAMCMYRpsmE0+0HwTMkytHxn6yXAU/QAAAOtnmycAAAAADGSYBgAAAAAD2eYJndrrfLSjPoNpGzKss5YMMoxVyzlnAABwNAzToGM9nsE09QxT60eGfjKM3Q8AAHA4bPMEAAAAgIEM0wAAAABgIMM0AAAAABjImWlwBPZ6uMCQ+6Z+QHsPGdZZSwYZxqg19P8NAQAA1s8wDY7IJhxovo0ZptaPDP1kWHc/AADA0bDNEwAAAAAGMkwDAAAAgIEM0wAAAABgIGemwRrt5xyj3g40l2GcWjLIMEYtZ6YBAMDRMUyDNdvkA81lmF4/MvSTYd39AAAAR8M2TwAAAAAYyDANAAAAAAYyTAMAAACAgZyZBgMNPaNoUw80l2GcWjLIMEYtZ6YBAMDRMUyDfZjy4egyTLMfGfrJsO5+AACAo2GbJwAAAAAMZJgGAAAAAAPZ5gkr7LWFaurnOcmw/loyHE2GqjrujLGT3bPXtYNm2E/W/a7b2dk57n0AAOBwGKbBHrbtPCcZptnPtmZYpx76GbtHAABgONs8AQAAAGAgwzQAAAAAGMgwDQAAAAAGcmYaW2evhwsMuW9KB7TLcDi1ZBg/w9D/Th/Uzs7O2h5msHxt6DoPHAAAgH4YprGVejtMXIZ+Mkytn23JMKZjx46NWh8AANgstnkCAAAAwECGaQAAAAAwkGEaAAAAAAzkzDQmYz/nJm3bAe0yHF0tGdabweH8AADAUTNMY1J6PRxdhs3JMLV+ppYBAADgqNnmCQAAAAADGaYBAAAAwECGaQAAAAAwkDPT2Fg9HKq+zloy9JNhnbVkWP86AACAo2SYxkbr7XB0GTY/w9T6mUIGAACAntjmCQAAAAADGaYBAAAAwEC2edKd/WzrOurznGSYZoZ11pJh/esAAACOkmEaXdqk85xkmFaGqfUzhQwAAAA9sc0TAAAAAAYyTAMAAACAgdY6TKuqb6qq36yqD1fVLfPXN1bVY1fce3ZVvaGqjlXVjVX17qr6war6vBPUf2pVvbOqPlVVH6+qK6rqcSe4/65VdVFVvbeqbq6qj1TVa6vq/idY8xVV9eqq+st5DzdU1curaucEa/bdCwAAAACbZ21nplXVDyf5sSQfTfL6JB9O8o+T/PMkj0jyhoV7n5jkN5PcnOTSJMeSPD7JTyU5J8kFK+q/LMlzk3woyS8kuVOSC5P8VlU9s7X2qqX775zk9+b13pXkp5N85bz2t1bVua21q5bWnJ7kyiRfkuR1Sf40yUOS/ECSR1fVOa21jy2t2Xcv22zo+Uc9HHIuw/ZmWGetTc5QVYeaYdXn7ezs+f+PAQAAcCTWMkyrqgsyG6T9fpJvb619cun9z1/459MyG4Z9NskjWmvvml9/QZI3Jzm/qi5srV2ysObszAZp1yV5cGvtb+fXX5rk6iQvq6rXt9ZuWPjY52Q2zLosyZNba7fN11ya5PIkr66qr9u9PvezmQ3SntVae+XC51+c5NlJXpzke06lFzK5w9FlmFaGqfVzKhmWrTMDAADApjrlbZ5VdYckL0lyY5LvWB6kJUlr7TMLv56f5IuTXLI7fJrfc3OSH57/+r1LJXYHWC/eHaTN19yQ5GeS3DnJ0xYy1cKa5y0OzFprr0vytiRfm+ThC2vum+S8JLs1F70wyaeTPKWq7naKvQAAAACwodZxZtrZSe6T2TbOv62qb62q/1BVP1BV37ji/nPnr7+74r23ZjaUO7tm2zSHrPmdpXuS5PQk90ryvtba9QPX7P7zG5e+rZb5gPDtSb4gyTcMzLVXLwAAAABsqHVs83zw/PWvk1yT5OsW36yqtyY5v7X2N/NL95u/vm+5UGvt1qq6PskDktw3ybXzb4J9eZJPtdY+vOLz3z9/PWPh2p6fcYprzpuvedPJ1qzqZY+6SZKqunqPt77mROsAAAAAODzrGKZ9yfz1e5Jcn+RRSa5K8lVJfjLJ/57kNzJ7CEGS3H3++vE96u1ev8cB7+99zdab0gHtMkwzwzprHXWG2uMhAsvXhhz0v7OzM6jW8jUPEQAAAKZkHcO0z5u/VmbfQPuj+e9/UlXfltm3th5eVd/YWvufA+rt/l9g+z2tej/3H+QzRl3TWjtrZYHZN9bO3Mdndm9KB7TLML0MU+xnXY4dOzZKXQAAgE2yjjPTdh8I8IGFQVqSpLV2U5L/Pv/1IfPX3W9r3T2rnbZ038nuX/XtsP1+xmGuAQAAAGBDrWOY9t7569/t8f7usO2uS/efsXxjVd0xs4cZ3JrkA0nSWvt0kr9I8oVV9WUr6n/1/HXx3LI9P+Ow1qzqBQAAAIDNto5h2lszGxh9dVXdacX7/2z+esP89c3z10evuPdhmT0x88rW2i0L10+05jFL9yTJdUk+mOSMqrrPwDVvmb+eV1Wf8+dSVV+U5JwkNyV5x8Bce/Wy9XbPU9o9U2nx91O5dtjrZJhmhk3qZ9WP88kAAADGdcrDtNbaR5NcmtlWxx9ZfK+qviWzBxB8PMnvzi9fluSjSS6sqgct3HuXJD8+//Xnlj7m5+evz6+qnYU1907y/UluSfLLC5nawpqfqIXhWFU9Mck3JXlPkj9YWHNdkjcm2a256KIkd0vyq/Nvyu06SC9br7V2+8/y76dy7bDXyTDNDJvUz2233XbcNeeaAQAAjGsdDyBIkuckeWhmw66HJXlnZk/z/LYkn03yjNba3yVJa+0TVfWMzAZRV1TVJUmOJXlCkvvNr1+6WLy1dmVVXTz/nHdX1WVJ7pTkyUnumeSZrbUbljJdnORxSc5PclVVvSnJvZJckOTGJE9vrd22tOb7klyZ5BVV9c1Jrp339cjMtnc+fynXvnsBAAAAYHOtY5tnWmsfyWzo9FNJvjLJs5Kcm+S3k3xTa+03lu6/PMnDM9si+qQkz0zymcyGZRe23a9hfO6a5yb5N0n+Ksl3JfnOJH+S5PGttVetuP+WJI9K8qIk90jy7CTfkuTyJA9urV21Ys11SR6U5Ffm/Tw3yelJXpHkG1trH1uxZt+9AAAAALCZyqynb1V19Zlnnnnm1VdffdRR1qKqsvh3bvn3U7l22OtkmGaGTeoHAACAgznrrLNyzTXXXNNaO2u/a9e1zRMG2dnZuf0w9V3Lv5/KtcNeJ8M0M6yz1royrHofAACAw2eYxqFaPhx9G749JMNmZei5HwAAAI7eWs5MAwAAAIBtYJgGAAAAAAMZpgEAAADAQM5M40iteiBBsrmHxMswjQzrrOUBBAAAANNimMaRWn4gQbL5h8TLsNkZeu4HAACAo2ebJwAAAAAMZJgGAAAAAAMZpgEAAADAQM5Mo0ubeki8DNPIsM5a68qws7Nz3PsAAAAcPsM0urTJh8TLsNkZeu0HAACAPtjmCQAAAAADGaYBAAAAwEC2edKdnZ2djT3XSoZpZFhnrYOuc0YaAABAnwzT6M6xY8eOu7Yp51rJsPkZeukHAACAPtnmCQAAAAADGaYBAAAAwECGaQAAAAAwkDPT2Bg9HhIvwzQzrLPWqt89bAAAAGBzGaaxMXo8JF6G6WUYux8AAAA2m22eAAAAADCQYRoAAAAADGSYBgAAAAADOTONjbCzszP5Q+9l6CfDumqteh8AAIDNZpjGRjh27Nhx16Z06L0M/WRYdz8AAABMi22eAAAAADCQYRoAAAAADGSbJxtr1Tlqm3q2lgx9ZThorar6nGs7OzvHrQEAAGCzGaaxsZbPUdv0s7Vk6CPDqdQCAABg+mzzBAAAAICBDNMAAAAAYCDDNAAAAAAYyJlpTMqmHFQvw/gZlh8GsO5aHi4AAACwnQzTmJRNOahehvEzrLLOWgAAAGwn2zwBAAAAYCDDNAAAAAAYyDANAAAAAAZyZhqTsbOzc+SH3h/Fum3McNAHAqz6O+LhAgAAAOyHYRqTcezYseOubePB+9uQ4aBW/R0BAACA/bDNEwAAAAAGMkwDAAAAgIEM0wAAAABgIGemMXlTOnh/GzN4GAAAAAA9MUxj8qZ08P42ZgAAAICe2OYJAAAAAAMZpgEAAADAQLZ5Mmk7Ozsbe1bYJmeoqpW1T3bN+WgAAAD0zjCNSTt27Nhx1zblrLBNzgAAAABTZZsnAAAAAAxkmAYAAAAAAxmmAQAAAMBAzkxj6+x1EP6q+w5y7agP/+9lHQAAAEyRYRpbp7U26cP/jyoDAAAAbAPbPAEAAABgIMM0AAAAABjIMA0AAAAABnJmGltpSof/H7TW8oMY9noww5BrOzs7x30eAAAATJFhGltpaof/H6QWAAAAsH+2eQIAAADAQIZpAAAAADCQbZ5spaM+r+ywM6xaAwAAAOyfYRpb6ajPKzuKM9MAAACAU2ebJwAAAAAMZJgGAAAAAAMZpgEAAADAQM5MYytt6gMIquqk11bds7Ozc1xtAAAAYP8M09hKm/oAAgAAAOBo2eYJAAAAAAMZpgEAAADAQIZpAAAAADCQM9PYOnsd4r/qvoNcW9c6Dw0AAACA/himsXVaaxvxAAIAAACgP7Z5AgAAAMBAhmkAAAAAMJBhGgAAAAAM5Mw0ts7Ozk53DyDwsAEAAADYDIZpbJ1jx459zu89PIAAAAAA2Ay2eQIAAADAQGsZplXVDVXV9vj5q6V7732Ce1tVXXKCz3lqVb2zqj5VVR+vqiuq6nEnuP+uVXVRVb23qm6uqo9U1Wur6v4nWPMVVfXqqvrLqrpl3tvLq2rPfXhVdXZVvaGqjlXVjVX17qr6war6vJP92QEAAACwOda5zfPjSV6+4vqn9rj/j5JcvuL6H6+6uapeluS5ST6U5BeS3CnJhUl+q6qe2Vp71dL9d07ye0nOSfKuJD+d5CuTXJDkW6vq3NbaVUtrTk9yZZIvSfK6JH+a5CFJfiDJo6vqnNbax5bWPDHJbya5OcmlSY4leXySn5p/9gV79E8nVp2hlqw+C+1k14bcs/uZAAAAwOZZ5zDt71prP7qP+/9w6P1VdXZmg7Trkjy4tfa38+svTXJ1kpdV1etbazcsLHtOZsOsy5I8ubV223zNpZkN8V5dVV+3e33uZzMbpD2rtfbKhc+/OMmzk7w4yfcsXD8ts8HeZ5M8orX2rvn1FyR5c5Lzq+rC1tqe37bj6C2foQYAAACwl005M213gPXi3UFaksyHZz+T5M5JnrZ7vWZfA9pd87zFgVlr7XVJ3pbka5M8fGHNfZOcl2S35qIXJvl0kqdU1d0Wrp+f5IuTXLI7SJt/xs1Jfnj+6/fur1UAAAAAerXOYdqdq+pfV9V/rKofqKpHnuTMsH9SVd89v/+7q+qBJ7j33Pnr765473eW7kmS05PcK8n7WmvXD1yz+89vXPq2Wlprn0zy9iRfkOQbBuZ6a5Ibk5w933IKAAAAwIZb5zbPL03ya0vXrq+qp7XW/mDF/d8y/7ldVV2R5KmttQ8uXLtbki9P8qnW2odX1Hn//PWMhWv3m7++b4+sB11z3nzNm062prV2a1Vdn+QBSe6b5No96iZJqurqPd76mhOtAwAAAODwrOubab+c5JszG6jdLcnXJfnPSe6d5Heq6n9buPfGJD+W5KwkO/Ofhyd5S5JHJHnT0lbKu89fP77HZ+9ev8eGrAEAAABgQ63lm2mttYuWLv1xku+pqk9l9uCAH03ybfN7P5LkR5buf2tVnZfkfyR5aJJ/m9nTN/cVYx/37j5asZs1rbWzVhaYfWPtzH18JgAAAAAjGfsBBD8/f33YyW5srd2a5BdX3L/77a67Z7VV3w472ZrTjnANAAAAABtq7GHaR+avdzvhXf/gb5bvb619OslfJPnCqvqyFWu+ev66eG7Ze+evZ2S10ddU1R2T3CfJrUk+sEdNAAAAADbI2MO0b5y/Dh0m7T4pc/n+N89fH71izWOW7kmS65J8MMkZVXWfgWveMn89r6o+58+lqr4oyTlJbkryjoG5HpbZ0z+vbK3dsuJ9AAAAADbMKQ/TquoBVXXPFde/Ksmr5r/++sL1h1bVnVbcf26SZy/fP7e7XfT5VbWzsObeSb4/yS2ZPQQhSdJaawtrfmJxOFZVT0zyTUnek+QPFtZcl+SNmT004fuXPv+izL4t96vzb8rtuizJR5NcWFUPWviMuyT58fmvP7fcKwAAAACbaR0PILggyQ9V1VuSXJ/kk0lOT/KtSe6S5A1JXrZw/0uSPKCqrkjyofm1ByY5d/7PL2itXbn4Aa21K6vq4iTPSfLuqrosyZ2SPDnJPZM8s7V2w1Kui5M8Lsn5Sa6qqjcludc8741Jnt5au21pzfcluTLJK6rqm5Ncm9kDER6Z2fbO5y/l+kRVPSOzodoVVXVJkmNJnpDkfvPrl+75JwcAAADARlnHMO0tmQ2O/nlm2zrvluTvMnsy568l+bX5N8V2/VpmT/Z8cGbbLT8/yV8neW2SV7XW3rbqQ1prz62qdyf5d0m+K8ltSa5J8tLW2utX3H9LVT0qyQ8l+Y7MvvX2iSSXJ3lha+09K9ZcN/+G2Ysy27r52CQfTvKKJBe11o6tWHN5VT08s0HbkzIbIP5ZZoO/Vyz1DgAAAMAGK7OevlXV1WeeeeaZV1999VFHAQAAAJiEs846K9dcc801rbWz9rt27AcQAAAAAMBkGKYBAAAAwECGaQAAAAAwkGEaAAAAAAxkmAYAAAAAAxmmAQAAAMBAhmkAAAAAMJBhGgAAAAAMZJgGAAAAAAMZpgEAAADAQIZpAAAAADCQYRoAAAAADGSYBgAAAAADGaYBAAAAwECGaQAAAAAwkGEaAAAAAAxkmAYAAAAAAxmmAQAAAMBAhmkAAAAAMJBhGgAAAAAMZJgGAAAAAAMZpgEAAADAQIZpAAAAADBQtdaOOgMnUFUfu+td73rP+9///kcdBQAAAGASrr322tx0003HWmv/aL9rDdM6V1XXJzktyQ1HHAXG9jXz1z890hRwePydZ9v4O8+28XeebeLvO5vo3kk+0Vq7z34XGqYBXaiqq5OktXbWUWeBw+DvPNvG33m2jb/zbBN/39k2zkwDAAAAgIEM0wAAAABgIMM0AAAAABjIMA0AAAAABjJMAwAAAICBPM0TAAAAAAbyzTQAAAAAGMgwDQAAAAAGMkwDAAAAgIEM0wAAAABgIMM0AAAAABjIMA0AAAAABjJMAwAAAICBDNMAAAAAYCDDNGA0VfUVVfXqqvrLqrqlqm6oqpdX1c4+63xrVb2xqj5UVTdV1Qeq6jeq6hvHyg4HsY6/8zXz9Kp6R1V9sqpurKr/VVXPqqrPGzM/DFVV51fVK6vqbVX1iapqVfXrB6y1lv+tgDGt6+/8Ov+7A2Nax9/VqvpHVfVvq+r/rao/m/97/Mer6n9U1f+vqswj2Fh3POoAwDRV1elJrkzyJUlel+RPkzwkyQ8keXRVndNa+9iAOi9J8rwkH0tyeZKPJvmnSZ6Y5ElV9Z2tNf8SypFb19/5JK9J8pQkH0lyaZJPJ3lUkp9O8rCquqC11kZoAfbjh5P8b0k+leRDSb7mIEXW+N8bGNta/s6vsQ6MbR1/Vy9I8nNJPpzkLUk+mOT/k+Tbk/xiksf49xo2lWEaMJafzez/OHpWa+2Vuxer6uIkz07y4iTfc6ICVfWlSf59kr9O8sDW2kcW3ntkkjcneVESwzR6sI6/8/8ys0Ha9Uke0lr76Pz65yd5bZInJXlqkl9Zf3zYl2dn9n9c/VmSh2f2fyQdxCn/9wYOybr+zq+rDoxtHX9X35fkCUl+u7V22+7FqvqPSd6Z2b/XfHuS3zzltHDIfK0SWLuqum+S85LckORnlt5+YWbftHlKVd3tJKW+KrP/d+qqxUFakrTW3pLkk0m+eB2Z4VSs8e/8t89ff3J3kJYkrbXPJHnB/NdnnnJgOEWttbe01t5/Kt8mWON/b2B06/g7v846MLZ1/F1trb25tfZbi4O0+fW/SvLz818fcQox4cgYpgFjOHf++sYV/+P5ySRvT/IFSb7hJHXen+Tvkzykqv7x4htV9bAkX5Tk99eSGE7Nuv7Of+n89QMr3tu9dmZV3eOAOaEn6/rvDQCb5zPz11uPNAUckGEaMIb7zV/ft8f775+/nnGiIq21Y0n+Q2ZnK7ynqv5LVf3fVfXaJG9M8ntJvnsNeeFUreXvfGZnAibJfVa8d9+Ff3bGDlOwrv/eALBBquqOSb5z/uvvHmUWOCjDNGAMd5+/fnyP93ev3+NkhVprL89s69sdkzwjyQ9ldpjpnyf5leXtn3BE1vV3/vXz1+dU1T13L87/pfOihfs85ZApWNv/VgCwUf5Tkn+W5A2ttf9+1GHgIAzTgKNQ89eTnsFQVc9LcllmB66fnuRuSc7KbMvb/1NVPzFSRlinoX/nL0nyO5n9Xd/9NubLk/xhksfmH76p89kRMkJvBv9vBQCboaqeleS5mT29+SlHHAcOzDANGMPutwnuvsf7py3dt1JVPSLJS5L8t9bac1prH2it3dhauybJtyX5iyTPnR9iDUdpLX/n5+dGPSGzp9j+VWb/kvn0zJ6m9S+SfGx+q29kMgVr+e8NAP//9u43VJOqjgP496ep7UKmliZUsGGBpkv2x5KicrHSwheLRhhWikRIvsjI6o3hmm9C0DDrRdA/TAgiMsEKBdtEkKLCKM1Eq9WiJFYtNd3NrV8vZq7crrt3Z32e670rnw9cDjNznjNn4MzdZ7935pz9Q1VdmOTqJL9Lsmmc0gX2S8I0YCXcM5Z7mufmNWO5p3lyFpwxls9Yiru7n8iwpPYBSV6/rx2EOZvXmE937+ruK7v7xO5e192HdvfpGb54npjkySR3zdphWAPmdt8AsLZV1UVJvpzkzgxB2oOr2yOYjTANWAkL4dd7qur/fs9U1YuSvC1DIPCzvbRzyFgeuYfjC/v//Ww6CXM0rzG/nA8neWGS73b3U3urDPuB5+K+AWCVVdVnk3wxw7QVm8x5zPOBMA2Yu+7+Q4bVNjckuXDJ4csyzHt2bXf/K0mq6qCqOraqjllS97ax/FhVvXzxgap6b4b/aO1Icvt8rwD2zRzHfKrq0N3sOynDZL2PJ/n8fHsPK2tP431f7xvYXyz3Ox6ej/byveZzGb7D/CrJqd29/RkNwH6ous3pCszf+I/p7UmOSnJDkruTvCXJpgyv7Ly1ux8a625I8qck93f3hkVtHJDkpiTvSvJYkuszzCN1XIZXQCvJRd199XNyUbCMeYz58djPMzyNc2eGcX98hsUHdiY506pXrAVVtTnJ5nHz6CSnZVgYZuGPINu7++Kx7obsebxPvm9gNc1xzE9uB1bTPMZ8VZ2bYRGx/yS5JrufA3Nbd39rzt2HFSdMA1ZMVb0yw1M0pyd5SZK/JflBkssWTzi6ly+dB2V4YuHsJK9Nsj7JwxnmS/tSd9+80tcBU81pzH86w3g/Jsm6JH/NECp/obu3rfQ1wBRVtSXJpctUeXpsLzfex+OT7htYTfMa8/vSDqymeYz5CW0kya3dfcqz7ymsDmEaAAAAAExkzjQAAAAAmEiYBgAAAAATCdMAAAAAYCJhGgAAAABMJEwDAAAAgImEaQAAAAAwkTANAAAAACYSpgEAAADARMI0AAAAAJhImAYAAAAAEwnTAAAAAFhxVfX+qrqmqm6rqkerqqvquhU4z8aquraq/lxVO6vq71V1a1V9ZB7tv2AejQAAAADAXlyS5HVJHk/ylyTHzvsEVXVekq8leSLJjUm2JTksyQlJ3pfk2lnPIUwDAAAA4LnwyQwh2n1J3plk6zwbr6qTMwRpdyY5vbsfXHL8oHmcx2ueAAAAAKy47t7a3fd2d0/9TFV9sKq2VtUjVbWjqu6uqkuq6pDdVL8iyYFJPrQ0SBvP/9QM3X+aJ9MAAAAAWHOq6utJzs/wNNv3k/wjyclJLk9yalW9u7t3jXVfkeTtSX6Z5K6q2pTkjUk6ya+TbO3u/86jX8I0AAAAANaUce6z85Ncn+Sc7n5y0bEtSS5NcmGSq8fdJ43lvUl+kuSUJU3+tqrO7O77Zu2b1zwBAAAAWGs+kWRXkvMXB2mjy5M8lOScRfuOGssPJDkuyZlJXpzk1Um+nWRjkh9W1cGzdsyTaQAAAACsGVW1PsOqn9uTXFRVu6u2M0NotuDAReVHu/vGcfvRqjp3rPumJGcl+c4s/ROmAQAAALCWHJ6kkhyZ4XXOKR4Zy51JfrT4QHd3Vd2QIUx7c2YM07zmCQAAAMBa8s+xvKO7a7mfRZ+5Zywf28NCAwth27pZOydMAwAAAGDN6O7Hk9yV5PiqOmLix36T4bXQl1bVy3Zz/ISx3DZr/4RpAAAAAKw1VyU5OMk3quqwpQer6vCqesPCdnfvSvLVcfOKqjpgUd2NSc7LsKDB92btWHX3rG0AAAAAwLKqanOSzePm0UlOS/LHJLeN+7Z398WL6n8lyceTPJzkpiQPJDkiyauSvCPJN7v7gkX11ye5JcnJSe5I8tMM866dleH1zk9191UzX4cwDQAAAICVVlVbsvyCAvd394YlnzkjyQUZFg44LEOw9kCSm5Nc192/X1J/fZLPJDk7Q+i2I8kvklzZ3T+ey3UI0wAAAABgGnOmAQAAAMBEwjQAAAAAmEiYBgAAAAATCdMAAAAAYCJhGgAAAABMJEwDAAAAgImEaQAAAAAwkTANAAAAACYSpgEAAADARMI0AAAAAJhImAYAAAAAEwnTAAAAAGAiYRoAAAAATCRMAwAAAICJhGkAAAAAMJEwDQAAAAAm+h8snh3T6Y0tBQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 720x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 453,
       "width": 617
      },
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(num='main', figsize=(10,10))\n",
    "\n",
    "data_path = \"./data/DeliveryGrid2016/DeliveryGrid.shp\"\n",
    "\n",
    "# Plot the grid of all tiles in Connecticut\n",
    "full_data = gpd.read_file(data_path)\n",
    "ctgrid = full_data.to_crs(epsg=6434).copy()\n",
    "ctgrid.plot(ax=ax, color='none', linewidth=0.5, edgecolor='black')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsEAAAI+CAYAAABddde4AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAABYlAAAWJQFJUiTwAABkwElEQVR4nO39eZBka17f97+f5zlrrrX2du9l5s6gmbgjmEECsWM2g5AIFiFhkAQ2WGAIjGRD/IQREQ6sCGOJTSCb1QNWIAkzwhIYJNAgYSFkgm00gGEYQAxz5259e601t7M+vz9O5a3uulXd1d2VndmVn1dEck6e5TnfzK7LfOrUc57HeO8REREREVkmdt4FiIiIiIg8bgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSCc66QWPM80AP+NBZty0iIiIicoc3Anve+2cf9MQzD8FAL03Tteeee25tBm2LiIiIiADwB3/wB4zH44c6dxYh+EPPPffc2nvf+94ZNC0iIiIi0vjoj/5ofuu3futDD3Ou+gSLiIiIyNJRCBYRERGRpaMQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSUQgWERERkaWjECwiIiIiS0chWERERESWjkKwiIiIiCydYN4FLIqyhH/9r8H7k1/wcPvmvX+Ra3uSa3+Qcy9ehB//cdjYeLifTxERETlbCsFAVcFP/RR8yZfMuxI5z37pl+CLv3jeVYiIiAgoBANNCJ4G4E/4BDCmWTfm8DV9f3R55/6TzjntMfdq+6RjF6Xt8/i57nfMadv+0R+F3/xNKApERERkQSgEA2F4uP5t33Z3iBF5VP1+syzL+dYhIiIih/RgHE3ojaJmXXfr5Kw51ywVgkVERBaHQvCBaQhWUJGzNg3B+gVLRERkcag7xIE4hsFAQUXudtJoD9P1k445brt+wRIREVkcCsEHpneCf+RHIEnuH2zuXJ++P27bw4Smo+v3aneW1z7Ltud57Ydt+6wpBIuIiCwOheADGxvw6qvwr/7VvCuRRXTnyA/3Gk3ipH2rq/Cpnzq/+kVERORuCsEH/s//E37hF5r1Bw07i75/kWt7EmoXERGR80ch+MBHfETzEhEREZHzT6NDiIiIiMjS0Z3gA1/91fCLv/j67Sf9SfxBts/qWF3vybheuw3f+Z3wjncc346IiIg8fgrBwN/6W82oECKz8mM/Bv/gH8y7ChEREZlSCAa2t5vlV3wFfPZn373vuOGyHnQIrZOOf9xtP2j7avv0x5/Uxq/9WvPQ5c7Og11DREREZkshGGi1muXKCly+PNdS5Jy5fbtZKgSLiIgsFj0YR9NnE2AymW8dcv50Os1yd3e+dYiIiMjdFIJRCJbZmYZg3QkWERFZLArBHIbg8Xi+dcj5oxAsIiKymNQnmMMQ/IEPwLvffbh9+rCT969/8OnOfXcu77fvaHv3a/es2zttu6e55tFrn1W7T+pnOW57njfr6g4hIiKyWBSCgdXVZvne9zYvkbPW6827AhEREbmTQjDw+Z8P3/ANh0/yTyc8OG55r32Pcuys21+kWpax/T/35xAREZEFohBMc5dOExmIiIiILA89GCciIiIiS0d3gg9cuwbf/M2H7497AOp++0573LK08STU+CS04Rx80zfBl3wJIiIickYUgmkC8Bd9UTPFrcgi+oEfUAgWERE5SwrBwGBwGID/h//h+GOmDzkd9/7ovtMcd9pzzuJa9zpP13r0a93rvEf9d37+efhf/hcYjU5uR0RERB6cQjBgD3pGX74Mn/M5861F5E7Tn02FYBERkbOlB+M4DBp1Pd86RI6K42ap2QxFRETOlkIwhyH46MNJIvOmECwiIjIbCsHoTrAsrmkIVncIERGRs6U+wRyG4KKAF15o7ggfHb5q+v647ccNdXWaNh7kWF3v8X2GRfm83kNVNe91J1hERORsKQTTjMMKsLsLX/EVcy1F5Fi9XhOK7zXahIiIiJyeQjBw4QJ8+ZfDe97TvDfmMGxM14++n/W+eV57kfYtal2Pe9/HfqwCsIiIyFlSCKYJF//4H8+7ChERERF5XPRgnIiIiIgsHd0JPvChD8G/+BenfyjrXstFPkfXm/8587zeR3wEvPOd6lohIiKiEAy88gp87dfCL/zCvCsRma3f+A34238b3vrWeVciIiIyXwrBwIsvHgbgL/7iu++SHb1jducDS/faf5pjTtr/INc5yzYe9dxl+K6epFqP1vwd39EMAbiz8/o2RURElo1CMLCy0iyfeQa+7uvmWorIzKyvNyF4d3felYiIiMyfHowD+v1mORzOtw6RWWq3m6VCsIiIiEIwcBiCB4P51iEyS9MQvLc33zpEREQWgUIw0Go1s8bleTN1ssh5pDvBIiIih9QnmOahoX4ftrbg7/99CMPjh5eautdwVY9y3Fm3cZbHzet7WOTPcBZtPM7vYTxulgrBIiIiCsGvecMbmhD87/7dvCsRmR1r4SM/ct5ViIiIzJ9C8IGf+in4pV86fH+vIbHuN1zWabadRRtPWruqbf7txvFhH3gREZFlphB84I1vhK/8ynlXISIiIiKPgx6MExEREZGlozvBB55/Ht797mb9Xg8cnWbf42zjSWt3kWt70tp9mDbe8hb4ru9q+gaLiIgsM4Vg4EMfgq/5Gvi3/3belYjM1s/9HHz1V8Nzz827EhERkflSCAZ+93cPA/Dnf/7d++710NGDHne/B5ce9rgHrW0e1zzuOH2nj3bN44671zW/67uaEVCce33bIiIiy0YhmMM/DX/8x8M3fMN8axGZlWn4TZL51iEiIrII1DOQwxB83EQDIudFnjdLhWARERGFYOAwBNf1fOsQmaVpCE7T+dYhIiKyCNQdAt0JlsXn/eGrru9enrTt6PYsa9rSnWARERGFYGA+d4JPG1xO2jat9VGOPemcWR170mc8izYW+XpnVfNZsBYC/VcvIiKiEAyHIfj974e/8TfOPvwct03kYVj7+pcxp9/2BV9w/MgTIiIiy0YhGHjTm5on5ycT+OAHH991HzTEnJdjF7m2RTj2uO3GKLyKiIicJYVg4I1vhJdfhmvXHl/IU6ARERERmR+F4AOXLjUvERERETn/FIKPGAyabhHTfrt39vc9bttpjnmUbWp/8a65rO0D/JW/Al/2ZYiIiDzxFILvMJnAl385/N//97wrEVlMv/mbCsEiInI+KATf4Ru/8TAA93qH/XbvfCjpUbbd+f64bY/a/uPYdr/P8CR8rof5DI9j2yJ/htEI/v7fP5xwQ0RE5EmnEHyH6UxaX/u18CVfMt9aRBbJ3l6z1KyKIiJyXth5F7BIpiF4OrOWiDTswf+nUAgWEZHzQiH4DtPpZBWCRe427SKhECwiIueFQvAdpneC1e9R5G66EywiIueN+gTf4SxD8J1DTB3ddtz7445XW2prUdqa/jehECwiIueFQvAdpiH4X/5L+Pmfv3couHPbcftEziOrvx2JiMg5oRB8h4/7OFhZgZ0dKMuzafPoMFen3faw+9SW2prldT73cxERETkXFILv8La3wa1bUFWH2x42bIiIiIjI4lIIPsK55iUiIiIi55dC8IEbN2B7u1k/OlvXg64/aW2IiIiILBuFYODmTfjKr2wehltWZx3GH/a889TGk17/0TbSFP7234bnnkNEROSJpxAMfOd3Hgbgp58+ecioB91+dP+jbj+Lmk4ayUIjXchpJAn8wA/MuwoREZFHpxDM4bBPX/3V8Nf+2nxredweJUg/CaF/EWuaZa2zqul3fgd+7udgMkFERORcUAhmuWfDUt9gOY2ybELwMv43IiIi55OGvucwBKsrgMjxpr8oKQSLiMh5oRDMct8JFjkNhWARETlv1B2CwxA8mcDe3vHH3KvLwMN0JzjpnLO8zsO0dZaf5XFd56z/beT1FIJFROS8UQjmMAS/613NS5aXftk53nQWRXUZEhGR80IhGPjsz4Yf/VHY3z9+//2GFTvtdp3zeM95GA9Tw7IwBj7xE+ddhYiIyNlQCAY+/uPhhRfmXYU8LosQ0J/Ec6yFVuvk80RERJ4kCsGydNSHWERERBSC5dyZTv5Q102odW7eFYmIiMiiUQg+8IEPwPPPN8Gprg9D1Enr99uv8+Z33p2SBH7yJ+HzPm8+P1ciIiKymBSCgT/6I/iyL4P/+B/nXYmcFWubUDyZwP/7/yoEi4iIyN0UgoF3vvMwAH/MxzR/Qrf2cErh495Pt8Hd++513knHP2obj3Lt03zGebbxsN8ZwD/7Z/BDP6RRHUREROT1FIKB0ahZ/nf/HXzhF861FDlDmuBBRERETqJpkzkMwUky3zrkbGk6bBERETmJQjCHITiO51uHnC3dCRYREZGTqDsEuhN8lqb9b6fDlJ3FtodtZzhslgrBIiIicpRCMIch+Ed+pHmY6s4Hqabhalbbzioo3m/b4/hMi0ohWERERI5SCAYuXGiWH/zgfOs4b46OEjHtnnCW2+53bKsFf+kvPb7PLCIiIk8GhWDgH/0j+NqvPZxh7HGEs1mfP8+aRERERBadQjCQpvBpnzbvKkRERETkcdHoECIiIiKydHQnGNjZgS/90tMde5o/9+uYx3PMaY776I+Gv/N3TteWiIiILI+lD8G7u/AZnwG//dvzrkRm4V/8C/jKr4RLl+ZdiYiIiCySpQ/BZXkYgL/9208+7jRDgJ3VMae1aDUtWj3/8B/C9vbhEHgiIiIiU0sfgqdT63Y68LEfO99a5Gz9H/9HE4LzfN6ViIiIyKJZ+gfjpiF4kSd7kIcTHPyKl2XzrUNEREQWj0LwwTegWcXOnzBslroTLCIiIkcpBOtO8Lk1vROsECwiIiJHqU/wQQjOMvj6r7/3sXcOx3V0aK57DdV1v2NPe+48rnm/Gh7HNe917L3Oe+WVZqkQLCIiIkctfQiOY7hyBa5ehd///XlXI7OwujrvCkRERGTRLH0IthZ+7/fgD//wcNvRrhH3er8Ix6qmk98/9RS84x2IiIiI3GXpQzDA2hp84ifOuwoREREReVyW/sE4EREREVk+uhN8h1/9Vfjv//vD99OHrI4un5R9i1jTvPYlCfytvwVvehMiIiIiCsFTv/qr8EmfNO8qZJaqCv63/23eVYiIiMgiUAg+sLfXLLtd+PZvP93DWMeNLXx0372OOcu2Tvuw2Lzamuf38Ad/AO9+N+zsvP48ERERWU4KwQem4wW/9a3w3HPzrUXOVq/XhODxeN6ViIiIyKLQg3EHNH3y+RXHzVIhWERERKYUgg9o+uTzK4qapUKwiIiITCkEH9Cd4PNLd4JFRETkKPUJPjANwTduwE/+ZLN+58NWRx+4OulBrJMe4jqujUdte97XPG3b865tNLp7KSIiIqIQfKDTaZavvgo/+IPzrUVm48qVeVcgIiIii0Ih+MCf+TPwPd8DL7zQvDfm7okXjnv/qMfNsm3VcPd7a+GjPxoRERERQCH4NcbcPVuciIiIiJxfCsFL7M7+s9NXXT/Y9sdxztNPw8bGvL8tEREROU8Ugu/wK78C3/qtT0YwfNRzniRpCi+/DGtr865EREREzguF4AO//uvwKZ8y7yrm487+s0f71U5HzXiQfdbeu2/ucftO6s/70kvN0GZXryoEi4iIyNlRCD6wvd0snYNv//bjA9o8w+FZ7DuprkX2FV9x+LCiiIiIyFlRCD4wnVXs7W/XKAKL6EkI7CIiIvLk0IxxB6aziuX5fOsQERERkdlTCD4wvRNclvOtQ46nO8EiIiJyltQd4sD0TnCWNS+4O3gd14/2uEkaRERERGTxKQQfmN4J/tCH4HM+59HaOm1YPvrg2tHti97G47jWq6++/lgRERGRR6UQfODNb4aP+Rh43/sOx9OFu8fXvfP9dP04dx4jjy6O4eLFeVchIiIi54lC8IEogve85+HPv19YPm2o1nGv3/fssxojWERERM6WQvAZUb9gERERkSeHRocQERERkaWjECwiIiIiS0chWERERESWjkKwiIiIiCwdhWARERERWToKwSIiIiKydBSCRURERGTpKASLiIiIyNKZyWQZf3L7P/EX3/lJ7O4P6PW6WANpnEANeE8QBgxGY6ra4/Hs7Q9JWy3wFdloyKX1VSLnML7GGsf29jZVWdBOWwwnGTv7+/S7bbqdNhfWVxkMB6Rx1KwPBtTe0+11MMZgnGUwGBMnKZNJRm0d3oXcmlhGOYRRRBQF7O3us7a6wrXrN7l05QK7u3t02x18XTGZTNjYWOeVV17l8pVL3Ly5xUq/R1mXTMYZKyt9rl27xeUrF7h+7QarayuMhkOcs7TSFjdu3ObS5Qtcu3aD9Y1VRsMRBnjX5/+TWXz9IiIiInIfMwnBVVVhDcRRxO1bW4ShpWy1aSUpgbU4IJ9MCKMEXIDHUNYVgbOYMGCU55g4xpclznrKqiIKQzCQlwX4mjAMSOKIIHAkcUhdlwRhRLe/wv7+HmEU4QLX1OIsRZFjrMUZSNOQKKiprKWqKlppyDiMCMOClac6WDvm6ae6hC6grgxhuEqeZ2y8YRVPRnvV0WrVVJXHdVOMLUnXLAkDLr5hhaqqMK2EMAjAwEbYwQU5mx+2grNAp4Wv61l89SIiIiJyCjMJwUVZ8sqr17HWUpQFadonDEPAEwaOqizottsMRhNu3bzF1u4O7U6bVrtFFBj2hgOqMqedtIjDiDRt4ayn02rhwpBtA91WGwOURUGn3QIDRVVye2uLdisFYDweU3kPeAaDIXle0Ouv4KuCwJdQ1FDVeJvTiQLKYkRgDQawRY2pHaExuLKkFVjyPMNaS5xarC2ovccFMJkMadmSiIwYGGZDojiCErJxhi8LWv0+1jShe/oSERERkfmYSQiu65qtnV1W11aJ0gQXhHQ6bUxVU1cVg8GIrKwYZgVlVROEAUVRMBoPMWlCJ4moqgrnLGEY0Om0yMZjhqMxZVESh+FBkO5R5jmFhTAKD0Iq9Fd6DPaH5EVOq93GB2CcI4wgSRM67TaT8ZDEWrKixBiDC0KiOMHXNdZYkjim9p66arpDtFstwjCiqkrCMCAvKsBT+xrrLN1eB+9hkud4YwmjmKqq6fQTqqoiimK897iwuQNsjbpji4iIiMzLTEJwGIZc3LzI1u42rW7KpMjZ2d0lMo6qLCmKkts7u2wPxzgX0un0CJxlMNpj7MeExrPS6RK6kCLPqYqCMAjxdU02GZPlOWVZEgUBUWiZ7I5ZWVnBOkMQBEwmGaPJmCIviJMEY8BZQ+BiwjCkqiswUHlPmrb40CvX+f0/+T2KahpQTfNBjMHapuYwCKh81YTwKKKqaipfEwQBVVlgnaWqaoqqghqSOKIqK8BQViVRFBE4B6Zp31nL/++jZvHti4iIiMj9zCQEW+tYXVtjMB5R5CW3RtvcrG6xvrJKO0lIoojVfo/SG3b2h1Rlha8r2mmLPB9h8Fhr8HjKosBhcM4SxQF5FpDnOXVdU1QldQ1h6DDWcO36dTY21jAWWmkKaUpd18RxRBA6rAmw1mAwZFlOEATUtsZ7y/v+6APktaEG8B7vPcYYMM3xxhh8XeEP1o3x1P7gAxsPHjDgPc3x2OkbAJxruoF4PNbY6WYRERERmYOZhGCAcZ7T7nTA1/i65vbtHa5ev8nG2jqXN9aIgoinL14kTQeMJhnD4T6bm6vEYYvQOeqqIstyrIeqKqgrqApDp5USBAGjyYQsywFPN2iRxDGtNKXX7VFVFUEY4KxjMBiQ5TlFURIETS611uGCAOssZZHT67SIwojaW1bXN3HWHXTHcJRVhbUW8E2wBaq6xllDUZY4a6nq6iAYG7xvknFVeqyBGo81hrpuQrU1HAZsEREREZmLmYRg7z0eKMqCfqdFMcl56tJFxlnB/nDI86+8yubaCj4vCENH2yWU1YQ0SYkCS2jBlDUusJRlQRpHOF9jvMfXMBiMGBc5SZJQFTlRFFBWFa1Wivc1ZVkCUJkKjCVNAqIoZjTKGGcZ8UFXhSAIcIFjrd9hpdvFtzr8qXd8Ahv9PtdvbvHUUxe5evUG6+srjEYjjHV0Wi1u3tziwqV1rr16g4sXNtgfDAkCRxCE7O0NuHxhjZeuXuepK5e4eWubtBXjq5rROGdjfYXbt25z8cL6LL56ERERETmFmTydVdUVQdQ8aDaaZKRpShQE9Notkihia2eX9/2nD/D8Ky9za3uLKAqJwpAgCHHOUeQFYRhSliU1ntEkb4Y/63YZZRmj8ZiqqhiOhk03g6qkrAqstdQegjBgb7CPsbC62qeqKgaDAdYayqLEWosLQ8qywlhLYC1XLmxgrGE8npDnBeNRc6d5MpmQ5zmTLGcyycmynOFoTFXWjMcZeVExmWRkWUFVVYwnGUVRMR43yyxvzimKkuFoTF6UDAZjJlkxi69eRERERE5hJneCjbEMRyMmWUY2HtNKEibjMWurq3TSlEGacHN7DJOMjbV18J6LFzZJo5CymFDXNWVVkeclYRzhjMEDedkMLdZqJVjnyIuCJInodlKMhzCMuHX7NhcvbtLtdJlMMqihLEo8ntpXJEnUPNB28CCbr2u8qbh8YZ1bL77KxdUuK92UfNKl30rIVrqstFOSwGIwxHFCudKhl0Zs9JulKVOsc03bnYRWHLDWTWlFlpV2Qhg2fZGpa9qxY6Wb0opm1hNFRERERO5jZt0h9nb3GY2GtFsp40nGYDCg02lz8cIGzhra7ZSbt2/jfUVR5Kyv9kkiRx0YxnVNGARY4yjKgk6/R10WjPOcLG9CcqfTot9pU/maJEkIo4jBYEiWFWSTnCIvMNbQ6bSp9muiJGZvb4D3nqoocM414dQ4qrLiqUsb/PH12wShwzlHkjRBO4wCjLPYwOG9IQgDgigkjAKiNCaOY/KyInCOKI5IWiVxktDpdklbKUVZEUUheI8xjk6nTZ7lpGkyi69eRERERE5hZiE4DCOiuMQFjv3hkKIsuL21RRhEXLiwQTpoQVVTFSXJShMm67IA71nrdcmznNJDGMfkk4z9/R2oPJ20DcaTRBEGQ5GXbO3sUNUlw3FG2koJXEgV1IBnNBoTxzGD0YiyrPA0Ywpnk3EzRrGDvCiIwoBOK+HatZtUazWvvHoL6+Dq1Wtsbm4yzib42pB3C1599TpJFHDt2nUssL27TxAEtNOUV2/cxBl48eWreDy3bm2TtpLmuJ09al/x8svXKIpyFl+9iIiIiJzCzELwaDIhDALGk4yqyJsh0bKC7b09auPYvLBJ2krY399jdXWFdpry/POvEpqaaKVHJ42pPZQ1XLtxk7woWO+vkEQxRZE1I0P4mv39fQpfEQSOOE3J8pxJlpHnOXEUUtU1gQnY29sny0uKsgDfjNU7Gk/odMJm9Iey5A2XLrBbQhoYOpGhHUA3tLRDQ1BZqtrTCjy9yJE4aAWQGE83AGxF6jydEBLnaTtPK4DUeRJT4axj4jyJqWkFhsT5+36PIiIiIjIbM+uY6pzBuoBiPCaNE3zlWV3ts7834MUXX+CDz3+AVpoSOEu/12MwHDIcjWjFIXlZ0kkTYuvIa6gNYCxxnGKsgdIxHI1opzHtdptxPiGKI+IowjmHx1PkeTNFc1VRFBWttEUUVc0DbLWnqmuMaR6UC5wjdPBhm32idh/vDReTFZIULnz4hWa0i25EEqfkec7TH75JUWSsvmGVbjumWgsoyxJvDJdaq6Sp5eKfukgSBzzT3cDTDIlWbDQPAL5p/WnQSMEiIiIiczOjENxMNjGZNKM4GBOy0u8z2NslOQirt7d3GIzG4GuejSLyg1nZrHUMsxxnR3RbKTs7+/iyZDAcUa7X1DWUvqaqm1Ejkiig1+0QBSE3b96i021T1yl5UZHSjOUbBGEz49vBxBl1VRG4pj+wtQbnHFEQUeMJyjFZltF2jqCuKSfjZkY4Y4Gc1BrAEruK4XiIi2pCYymNb8Yn9hUtoApqImso67KZ2KPI6UYhlppslJPn+Wy+ehERERG5rxnNGNfMluaso9vpUk5GBEFAHKe0Wyl1XbG6tsYLL71MXmQMBiNeeOEFNtf7BGFIYG0zRNlozPraGp12m1F2MB5wUeJpJqEIw6CZ9AJDnheMxmOKsiAIQuq6Ik1iqqJkMhmTFwWVh0k24cLGGsUko6oqrLFkk4y8KDCBIy/L5g510EzY0eu2qeqaKIooiwrwBM5hbEy336XIcwaDIWXl6bRTwlYKHsbjMUmaEgeOosgJwpQwCKmqkrzMaXdas/jqRUREROQUZtMnGIijGGst1lnWL18kjUKMh7oq2dvbI05LLly4wHg0ZjwcsrOzQ5pEtJOkeXgtjIgCS5wm5MMxnTii9p6sbPoCJ1HYTGlsLUkSU/qaVqsJlu0koQrDZhSGJGaSZcRRRJ3ltFop1hiiMKTX6VDVniiK2RlM+L33fYCy9oRRSJJEBM5RVzV1XeMxzXBq3lN7j3OGKIoYjzMmk4xuv4cvK5yzWGMJwoAo2sJXFS50BGGAr5uRKSZ5E8D52Fl8+yIiIiJyP7PpDuGhriv2B3vEUczm2iqvvnqdTpoSOUe/36fwnjRNKfJmZIanr1yiqEr2BvtsrvRwgSUIQoosJ4wC8DAYDtnf2wdjqOuKJImpKsdoklH6CheEXL64SVUUVHVFUZb0ox5pWrE/GOCcBQ+GZtxh6xxV3UxaYazjt973n8hqjzdN4AWoPVjTfCaMb8IwnqZThAFjm6mSzfSgpn1jml6/dV1jplMl2+a63vume8VXzuTbFxEREZH7mEkIrqqKIA7pBl0219cYDYd4XzczvUUhLgioPezvD3CBwxhHr9tje28HbyAvctLAEUYhrVaL3Z1dDJ7haIQBqrqZPrnIS4qyxAUh+6MhSRqzXpTEUYDPa9qdNkVREEUhvZUe+4MhTO/m4smyDOcsta9YW12h3+uyPSlY27h08DlqnLWUZQWmCfYA3jfhFu+pDvZ5X4NvhmUzxgCmCcd4/EE4bt43DXg9GCciIiIyNzMbIu3mzS36vS57u7tEYUAQOiyOsvZURUmv36OqPbe3tpm4nCJvpkYejUYE3tOJUwyG/f0h40lGu5VijCWO42bYNaCuPVleMM728Qass9zc2qLfadNppUzG49cmp9jZ2yNKE4IwADzZZIL3NVUFURQBsLm2iisqnvuoP8dKp8v1m1s8/cxlXnr5GhcubDAYDDDGkaYxW1s7XLq4zquv3uDShU12dvdI04gayMY5aytdbtzeYXNzndtb27TTBF/XZHlOr9/hxo2tWXz1IiIiInIKdlYNl1XJrds3ubV1C6jotFu0WwlVWeICS1WU3L51G+8hLyr2RyOu37hJNsmIopgsL/DGMJlMmjun3tPudOh0uwD42mOsw7mQKI6JophWmgLmoE9vTFlW7O3u0em06PY6VFVFVVaUZYmzFl97oihswnBZ8OY3PEUchkyygiwvGI0mFEXJaDwmm+SMxzmTLCfLS4bjMVVVsz8ckxcVw3HGOC8pS8/+MKMoa/b2RmRFxWA4YTwpmWQFg9GErKjY2RvM6qsXERERkfuY2Z1g8MRRCL5mMhmThiEWy1q/A9ZRFAV1WRJECaPJmKIoCIOACxc2GY5GRIFjkuUEgWv68hpDUVaUVUm322YyyfBYioNgG0YhWV6Q5xNWe22iMGTM+GDq45ggyAmCChcEGGspqpKyqsgmGUWZA5b1lQ7++YqVbotup81af0waBay0W3TSGF+VGGvppDHjNCYJLd00Jg4trSggsoYwMLQiR+QsaWRxvqYVBaShpa6hChyhNSShm8VXLyIiIiKnMJMQbIxhPB7R2Vylk8bYwJEXOWvdPnXV9MXNi5J+vwc2YGtnu7lDHAXNrGydNt6AtxAEcfMYW9X0x7UHk2HU3oOpm5EgqoOH5OoKPBR5gXOO8TgjiBzjyYTJJKMsKnxVsb+/z7VXr7G+sY4vSqy1FHnF+kqfOAyJo4goDGi1W4RRRLvTotVKDh5oM6TthLVqlbTVZmN9jV63DRiSNMIaQ6vVpttp4TGsrPZJ05g0Sajrg4k72gmhC2fx1YuIiIjIKcxsxrgojuBgZrbBcAhJQtXpUGU5aRBR156d4ZCiybYkrYTIGcLI0eu1qfKcdqdNGsYUk4y8HtNqpVy7cZ2NtVXKoiIvCpyzdNotsjzHOHDOURQlL7x4lZ3dXSpfUdWezY0Ntq9dB+tYX+3x1JVLVJU/GGbN0Ou1aLW7bKz0uX7jNpNeyfUbNzGh5eqNWxjr2N3dxzlHa5Jw49YWla955ZVr1MDNW7dJWylhELC7P+Tixiovvvgql65c5NbtLTqdNr4sGY0nrK6v8MrLr8JnzOrbFxEREZF7mUmfYI/nwuYmLgwpy5o4SdgbDri1dYs4DvBVQVXkUDcPqGV5Rp5NWFtbIYoc1kGSxATOETiLc4YsyxgOB8RhwObmJr6uKctmNraqrPDeUxYlzjrG45yXrr7KcDzG2QBnLLu7e/iD0RzCMMQ5R+1r8DRTHgP4mjc/dZG2rUhtTSc0pNbTsp7UeFJTkzpP5Ktm3XoSWx28alrOE5uKyJfEBhJbk9ia2NTNduOJbXNMYv0svnoREREROYUZ9QmGGze2KKucdrtFWdf4umZvPKKXpgS1pSoL0jhkMBiST8a0Oyl7g30C1ybotWiFMVu3b9HvdohcSKfdZjgcEAUhr77yajMRh7V4ajCeVitlPBlRe8/e/gBjDL12l0sXN8jznLKu2NvfJ0lbzfsiAwxREjEcDvF1RVWVfPgzF7kwLjA25KnuKmHoWX2mTxTXXEibGe9cEJCtR4RhyJXWJYIg4g0rl4iikKqqyDdS6tpz5a2XwRqe7qxRe0+aJBR5ThSFvKn/zCy+ehERERE5hZn1CZ7kJVVVUZb7eN+i200Zjgbc2ttlLelwe2sbGwRcWl8lDiCrS/b3dklCyLKETpTQSlKqymMctNIEA+zu7jZBsyyp6prQhlR1QTEeMx5PSJKYOI5Jo5AkjpqH9CyMJuNmTGIM+SQniUPyosIZg6UZAzgILGVZsJYGZPkEqgntOMXZnMHeDkkcY7IcMkPoPa62MJ5gw+Zhu3JSE8cRIcDBc28uCMkmE2pfQTbG1DVRFRE7jRMsIiIiMi8zuxOc5xmBs/R6XZxtxvTFOiZ5yb4fN7PFlQWdJKabbHJzb4dRBa12DMZT1yVxEOFsgK8hK3OMsVS1Z29/n0leYJyjHQREYY0HijynlcQEztJpt2mlEbvbO8TthLryrPZXsAZWV/v4qiQIasIgaNr1kJcVURThq4ogCOivdPG1p9fv0mqnVHVFHEUURYH3zdTOzjWBufZNALbOEsTBwVBwAd7TdOlIIjCG8WjUTKxhFIJFRERE5mVGD8b5Zjxf4xmPxwTOEgQWYx1hFGFwWFsRBQFRFLDW72EDuLZTUlclzjiiMCRyIaENGe4PKIuC2nuqqiaOY/K65ur1m1y6tIn3NUVesL7WxzlDr9MhjSLiOCJwlgpPVXmiICRJU7KsIA4dBo91ljwruH79FluDEWEU0UwGVxM4RxhG5HlBEARY04x/PBiMaLXbVHXNeDwG20zHbA0URUmrlWCNwbrm6y2KgqqsCKKQyWSCs7aZfe4TZvPti4iIiMi9zSQEO+fw3jPc36PdbtHr9yjyjPFoRCdMaPXa+LJmMNynPjjHenjrG59lZ3+PyWDCJMzxzhO2Q6w1xHFEXlTNo3yVJ4kjrly+wNb2Dpc2Nrh9e4vhcEzgoNftsHlhg7qqCENLXpSAYTIeUXlPK40xZU2Zl+RlSRRFZKXnP/zm/0dtHLXnYFpjj689gQuaqZZ9jffw2ozHpnmwzntL4By19/i6pvYVxhjsQTj2HoyxzVTLB1MsW90IFhEREZmbmYTgsizZ393FAHVds7c3oN9tk7ucrd1d+u02oXUkaQrAaDJhZaVPYB3PXLxCXdcMByP2hvvs7OwSOsfe3h6rKys4a3FxRD9JiZKYVrtDv9vFYLh+8zppKyGOIwajIUkUkaYpUVRRlBXWGKg9zgVUdYkNmnGJi7qmlaSEQcS48iTtLniDP4jCmINBNHzdrHuPxxMEzSQezgVYY7HWYZ2jqioMzV3jJhR78Ia69hgzTdEaHUJERERkXmY2TnAQRs2UxFVJXmfsm2aq48pXDLIJviiJnCUrWzjbjNVblCW7+wN63Tb9XpsyjxmNxoCh3WoBHExc0SKJYsaTnOvXr7N1+zb4mvW1VazxYGBre4tLFy+Ac0RBTKdTU3lPGKeUZY0BirIiL3LG4wkrvQ4rvS6Rt3zcp34W3XaPre1dLl7Y4Nr1W2xsrDAcDHGu6cKxtzdgY3Odq9dvcXFjjb29IUEYEoeGre1dLlzc5OorN7lwcY2bt7bodDvUVTNEW6/f4dq1m7P66kVERETkPmYyTjAYjAFLzaUL6/R7HXxV0m63sMaytbtLkMRkVUXhK4q6bvrXTnL2BoOm/y2W0AVc3Fhnrd+j3+0CHmMcrVab3b0B167f4PbWNlEUEUbhayNI5FlOmrZI04SyKInCsLlb6wIwljwvGE8yiqKkrj0uaH4XuLS5RhiGDMclVeUZDCbkpWd/OGGSVwwnBaNJwSSv2B9m5GXN/v6QoqwYjsaMJxmjSc5wnJEXFfuDEaNJwWiSM84qxlnF/nBCllfsDcaz+epFRERE5L5mEoKNgcAZkjggG49Y6XYIA8dqv8PmxgrOGqqqIstLbu/sU3ooqhrrLCv9HmEQYoxh4+IF2v0+obOEgaMqS6IgJJtkYAx7gxFguX5zi6vXbjIaT0jTFlVZU5U12zu7mINrOeeYTDKKqumva/BY40nimCgIMb7m2Q+7gsETBwZnDWkcElpIQkdgII0CktAROUsrDomcoR1FJGFAOwrptRK6acJKu0UahfQ7Ke04oBNHtCNHJwloxyGRhXY8s5vwIiIiInIfMxoizZMmIb1umzgM8HWNs1CVOYEzdNotMJ4ojRjnObd2dtno9+glEYEx+KrChSHbN28ShhHjrKAsKpxxzZi+ZUU+yalrTxhHZHnB+CBs7+8PmUwyXn71OmkSsbG+wrNveAO+9hgDo8GAbiulrpph26o8ByzOWvqtFrGzXNpco9ftk4SO3koPg6ff75JlGXgIwoA0juh1O1SXNmilMax1ca4ZL7jbbZNEESsrXXqdDnhPp9OmKErotQ9Gj1CfYBEREZF5mUkItgZ63YROJ8FgMB6iKGQ8GdHtdojDiNFoTFYUOBMwHI0JraUTr7K62jwgNx4O6bRbTXiuK6I4oqwqfF1hnaPbbvGWZz+M4XjC9Vu38VUzTFppDC4I8TXUJdy6tYvhZTqtpitGGBqCIKAqC6oqxzoHBiyGTpqyvtrnhVdeZX0l59VrN3nqqZIXX77KlSsX2d3Zw1hHEofcur3NlSuXePnlqzz91CW2tnZIWynOGnZ29sgqz0svXqPGcu3Va6ysrTCeZJRFyfr6Ci+/cm0WX72IiIiInMJMQnAUhVjg1u0t2mlM4ELW1tbIy+YhtNA5AufIsoK8yhnuDYmcY3t/QFEUrHW7rLRbBIFjOB6R5xNybyjKgjSJsMYwziZUZUloDSudFq04ZDgagTG00xaTSU47bZEXOVtbu1B7gijA2oC6rnHOYaxtHtarcwyGMIq5srFGYTyxrUisJ/AFERWxrYlthXOGkJKIgsRWxJQktiYxFS1bYy3kztOyNS1XE5uCxNVEFNigZlKWpK4mppzFVy8iIiIipzCzGeNacQtf10RhRFVWbO/sYazHWRiNJ0RBROgC1npdAhxFUbK9u0cUrlFUFXlRYg3cunmLGksSJ1RVxWA4JAwLRpMMfE1ZVownEzqdNs4aVtfWuH17GxcE3N7ewRhw1hCEId57wjA46CNsscZQ+ZogCPC1pywL3vTUJmGrTxCGPN27QBA4LqWbxLGjaPWwFqIgIFuLSJKATbdGGJVcjNuUB/2Nn11bJYocV9oXiOOIN/QvUtVNn+eqqEiSiKc7T83iqxcRERGRU5jROMEVZV5isdRVRVVVFHlOnITEYUTkHBjHcD+j33U8dfkSVZnjbDPGcFGUxGsJVZGTpjGYgN29Idu7O7TbKd12h3GWMRxNqL1nMBhQ42mnKePRmLr2TX/fuiIOHIGz+Lqm1W3jrSUMA+q6wDqLMabpqzudvKLMYbxLPvbgK+I0pQ24Mie2UFcVrTil07IU+ZC2zclH+yRRSIWHusbmE4qJPxgnOSCJIrIiw+ee0BhMNsFk2Sy+ehERERE5hdnMGGctdVkTRhEYSJKQOE4YT8Z0220mwxFF6RkPR7xaXOfShQ2eunyBfDIijUM6aYuyLFldWaEoCl599QZF5cmLkmowwmDhYBa5ZhKKmv3BgDwvKKsd6tqT5xkbq31aSYwxEIYBRV5Q1DXdTgdnHUEAWV5gncMYg68rwtBhDFRF3qz7mjBwlFWBtZZ2GmONIc9z8DVJHOEchEFEUVVAfTBT3EG9ZUkJVFVBEIbNlx460qQ/i69eRERERE5hJiE4jiLiJOH5l1+mt9qm12pTlSUWQ5lXpGFMtr9Hv90mbbVIkpitnT2iwODwRL0+xjdTFkdhTL/bZTAa46yhyAuG43EztFnUPCwXOof3Tf9h5xxxEGKjkDgKabcShoNR0/e3KvDW4WuY5BU3d3aY5AXBwZBsVVWDGeG9J0li4hJ2RwOCIKSqKqy11PUQay15nlNWNWXtD6ZnDvDAJMtot1LyomA8zsjLgqIoD9o2+LomCBxRGMFHz+LbFxEREZH7mUkIzvKCl169Tpgk3N7apyw8K70OnVZClmU461jvdXBBxN5ozGB/SJbnRJEjDmBztU+athiPxxRlQRxHTLKMbiulrD3jIqesKrJBE1CnsxCnSUIQONIwIo4jAmdpH4wKAR5rDGEcEycp125v8zP/5j8wqWo8lsP+EFD7GmctURhSFAUY39x99geTHfsa4yxVVeExGG/wBvx0KmR/MDWyacYoBg6nXgYMHoOHr5zFty8iIiIi9zOjGRs83hgmkxxnQ7K8pK4hzwsi57DO4ZzF1DVpHPHy9ZuknQ7D4Yg6tbx88xp5f4VOlGBrqH1FEDg67RbjvODm7i7OOdI4JokTyrLZP55k1FVNaStMntNZWWE8nhAlIb1+l9E4x7iAPM/YXFsjiiKG4wxjw4PcajAcxGFrqXCY0OFrD9ZiDoKyOcjLQQDGOOq6xhgoq2bEB++bY7yHwPnmPONeC+vOHqyIiIiIyFzMZpxgZymrZhQEak/ta8ZZRhq16XY7UJXNw2l5CUVOK46bqeucpd/rgYGiqtjbH5CEzZBomxc2eOHFlwmCgCiKuL2zx7PPPEVRNiNGWGsJw5A0ikjikLKsyLIcazwuMIRhgMkKwtDhD7pQrPR6lEHOJ3765xBFMYPhmNWVHtvbe1zYXGVvb584SbDWMhqNWVvtc+PWNpsbq2zv7NFupfjaMxqPabVSrl6/yfraKlu3t+n3umR5jgHa7RbbO/tsbqyytbNLv9chz/JZfPUiIiIicgozmTYZ73GmJh9PSOKQdpqQTcavXc0fhNzKGPqrPd749GW6acyFtVV6nS7OBtSlJ0lSkjimKCtevXYD76GqDoZdq2oG4wzrAtbWNkiTFr4GX9eMxtlBP9wK72u6nXYz1TIQhhFBGGAxPHP5IkmSkFUBBCnjHGzYZjjxlASMspq8spS1ZTAqKHHsDXMq7xhOSiYlZBWMs4oax3BSUxEyyTxZZShqyyivqUzI3jCj8pbhuCAvoaxm8s2LiIiIyCnMbNrkjZU+ZVVRVhWXLl1gf39AmiaUvgLvGUwmOJo+s2+4dIVer8vOYIDBEhCwubpKhKEomru5nXYbZwPyoqSdJji3wWA4wbmIwdYuxkAnjUjTtJneGKi9x7qATrdDUdY4X1LVNWVVEQchm2sr/NHVW5i6xOExvmqGOMPjy5LQQuBo7iYbCAyE1oKviQNHaE3TvQFwBiJnsNSEzmDqktCBL2ssFcHB+ZGzWJpJNURERERkPmYzRJpzXLq4SVmW7A4G7GzvcPnyRUbDISZspix21pJPcmztefnqq2xubLLS7ZMXGe20hfNNAJ5kE1ZWViiLkk6rzdUbN7EYkiji5Ws3iJKEwXhMFFjS0FGWFe00pSgL8DUrK33SNMXmBUXpsdZS5DnOQL/TIokCPvyNl0nSFiu9Fqv9Hr1WTLebMl7rY6whcI7VboeVXhsHrK52aSchYRAB0G0l9HodgsCxttpjpZ2QpgnGQJEXtNop3TRhpdcmiRztdou61q1gERERkXmZSQguypIszymKHO9rjDOMxyOcM8RpxHhUUpUV7TQhDSJ85anqmsFwD2MN5WRMK7CkUUQ3jpu+td6TFwV7e3vkVY3H0e10uX7jJp1Wi26akkQRaRLjrKGuS8ASRyG7u3sY48jzHBenB+MCW1Z6HdppzB/9yQv0e32u3dzizc9+GB/84MtcuXKB7Z09gigkiUK2bm/zYVcu8Ccfepk3vuEKN25tk6YtQufY3tnl8qUN/uSFl3jm6ctcf/Um/V4XYw2DwZC19VVeeuUaz1za5PrNm/RXV7C1HowTERERmZeZ3Qne3t0jTUNc4EiiGO8rkjSlKkrqsiKJAmJjcd6Ttto4Y8iyCZ5miLSd/SFF2jzA1tw5NngPnU6b27v7TPKc0DniwHFhfZVuK8HUFfgKsBRFTq/Xpabpf+yswbngYGgHw3gyofbwzKULVMbj6oKAioAK5wsCagIqQiyRdThfYX1FaGsC4wmoiSw0Y0iUOCoi6wmpiGxNaGoArC8P2qlwvnlFpj4cTk1EREREHrvZ3AkuCiZ5QRSHhKEjDBxxGGCALMuIwoDAQzsOiVyEMZa93V2CwNLt9oiiiNDG5FmGxdBKEyZZxv54wt5oROWhKAvSNCFwbQJnqcqSOLQkcUSeZThr6a/2KH1NMc6Jo4SsqAlbTV/hLJtgreWtH3YZH/eJk5TL/YQkgdazGySJYzPpkMQRdV2xGfXpryRspldwgWXVdQjCkMA5iv4aaRqyEV8kTVMupRu4IMBaQ7mREoYBl9JLWGfZbK0RRxGV7gSLiIiIzM2Mxgk2JHGMpxkNt9fpUpVFM2qEdZRZQa/dYq3TZTAc4yJLFAbU1lNVOeAwxtHvNHeI8TWTLGN7b49bO7sEUULaSljtdciy5i5xXuQ4G9DqrNHtthkOh/i6pgY67TatVotbW3vNzG/G473H+5rAVNhyHzMaYcdjhrvNdMmUAa6qoAibh9nqknp3iKs9RZETuwBbWcKwOa7a3SGoSqo8JPA15M1nj4MASnCTSTNSRRLjfMFkMpnNVy8iIiIi9zWj7hCWOEkp8hxrYWt7j9V+B1/VBC5g8+Iak5099gdjqtpjvGfzwiaj8RDvK5yByXhEHQS004Q8y9je3eP29g5V3YywEEUhvV4HKs+t21uUVYmzMds72zxz5TJQgzU4oNPrQAVJEuEC00zHHEVgfHPXOCiJgoDIVMRpcDC9cUUchzhrwXts0Mw6V1uPCx3GQhg4wBPFIWVZklMTBJaqbu5UYyAMIpx1BK6ZJKOqKjCeVhrP4qsXERERkVOYzWQZNiDt9DDDIeDp9Ve4dfMmezu7rPS6TPYmbK72uH57h36/T+ChKmu67Q7j0YjJaMJoNKKqS565fJndvX0mWUYrjcmKisHeNs50GexHxFFEK40ZjAqcM3hf4ano9XsHQdQyGY1xrpleuSgLjGlmbHPGkKQp4AmCgDAImZQl3hi8NQyLkiAwDIdDirIZ7i2JY7KswBiw1pLl+cH0yDVVVZPlOdZZ6qrGG0NRlM31fIU1hqIsD6ZeNvAps/j2RUREROR+ZhKCy7pmkvvmkbEyp6qg11thtb/K3u4+LnRs7w0xNmA4nhCEIdYYRvsjjPF4bwiCZka5STYhCiNaSULbOYIoYmtnmyQOiUKL9xXtNGKld4Fep401hslkQhwnxHGCcQG3b95mku/SW13BWMskm2CoiAJHZEKKosTjMTbkX/+7X+Xq7R2Kuqaum/7DZVnjDdR1jbW2mUYZ8NQHn9gCBmOgquuD6ZcNnulMyR7jofk/zTjKAHztLL59EREREbmfmYTg2ntKHHHaodjb5ubtbS6trTIZDNjZ2WV1tc9Kr0ddFkTOMh5nDKIRgQWqmt3dPeq6pNWK6XU7lJOMssipbRMsV3s9otBhvMc5w/pan07aIk0TBoMh42xMVZU4G5ENR7RaKUEUsr2zS5i2CQKHryuquiYvC8bjEZ12mzQOCcKI3dEEbIQxjtr7gxElLFioMRgH4DEAGIwFQ9NdwhqLsQ68b16YgyRcY2xzhq/r4742EREREXlMZhOCq4qd3R3e+qc+nCybEAWWvcGQapLR7rTAwDibMBqMaKcxK90OxhparRRf5iShYzCaUOaG8WBE4By1h/EkZ2d/n163QxykJHFIWRTs7uzQbbUYDgaAZ6XfAyxpkrC/N8AFrunqkOUQRvR6HfJx0z0hz3OCIMD75oG3zY01wpdv8Emf8dmkaZcsz2mlCXv7I1ZX+mxv7dDrdxlPmtEl0jRhb3/AxvoqN2/eZm1tldG4GZ0iiAIGgxFrK11evXaTC5tr7O8NSFsJ1phZfPUiIiIicgozmbzXWst4OOaVl69SVZ6d3T129gZEaUpR1gzHY7yBbq9Dp9shSRLCKKTbadNptXjTs2/g6cuXSKKIqq7I8wKPZzQe42nC6mAwIHABFzcvsrm2TlWVdLod4jgmimIGgyEvvPACaRqRtlLiVkoYOKqipMxzOq0WURjg62YWuaqqqOua9dUecRxhwxRvY8o6IEw6DCc1YdJiUhq8i6hxlJXD24jBqAQXMRhXFLVjlFVMipqqtuwPCzABg1EGNmA4KZnkdXO3WERERETmYkYPxll6vT55XgDNg2C1MYyzoulOYBxZluODgLJsSjB4tre3CQ20Wym9Tpsim1DXNWsrq1jniOKY3f19WmlC4AxpmmKMIQhDXr76Mk8//RQr/RUm4zFRHDEYDwnCgLIo2drZJooiiroZuK32HucC4jiiKIrX2lpfWyF0TVcLX1VN14W6Au+p6wpLhaXp6uB91fT89U2fX2tMMwycMViafsHOQtN1wuPrmsBZrKbKEBEREZmrGY0T3ITaVpowHO6B95RVxWA4JAwDijznwtoqZZ5T1zW7e3s406aTpgxGQ7qdlDiMabdaeDxVVdJup4RxhD8YMziJUnzt2R+OGI1G7O4NSG5vNUEzDEniiI31DcCQFwWdbpvBeEwQBThjCZxjMp5ggCgMiMKQrCgJjCFNIi5urNBu98nynHY7JQoDNtZWaIWONE2oVntUVUWSJrSjkLW1HpGFfr9DlnfxtSeOI1Y6CWurXZLA0e/3WO12CMOAMApn9dWLiIiIyH3MJAQ3j5HVlPmYdhJjafoJT19BGGCtoygLWklEGkd4Dzs7u3Q7KTYImxnhWinWGMIwIAgjBjdvEYYhe3sDgl4fvKHMS8bjCcY5bty8TZ5nxHHEU1eu4AJHXlW02i22dpth1tbWu8RxhMFjLTjrmmfY8ITO0k4CVropL796nW6vZG9vj/XNdV586SpFXXPt6nVWV1coqpLhcMTmxjqvvHydN1QVL770CpevXGB/OKIuPd1Oi1ev3eSNb3yaF1+6yuXLF9m6vUO7FZOkCXzkLL59EREREbmfmfQJrquKKp+QT0ZEzpCGjr3tLaq6JIojqA3bu7tUtaeqm24JzjiCIGQyyfngh17k9vYOSZLiXEgYxozHE/KyIAgC3vqWP8Xe/j57+/tkWUZRlNR1E2STJGZzY50gComSBDzs7u2T5zlgKMsCXzezxeE9eVlgrSFwAWEY4gw8+8xTRHgcFZaagJoIj60rAlPjTI3zFc6XGF9gfYGjmX0uMB7nSxwljoqACqocVxcE1AS+hKokUpdgERERkbmZzZ1gA4PhAIxv7uIaQ0nNaDIhjmKcg+FwTBQF1Hv7rHTajLOc8WRMXdU4W5OVJc6FrPR65GXFaDJp+hZXnu2dXcqqpigr4iii0+2yN9jHBY6qromiGOsC6qrprxs4R+QCqtqT5wW5y+i0EsDgDh5QC4OAuqowVLz1mUuUYRdjI4peQLudcCG+SJJEXEg2CMOALHewEhMnEZvBOv2ViNVgkyhyXO70yLKcMHBsput0uxHr4SZJGnApXaXGE4Yz+f1DRERERE5hNiHYWowLyLIJN25v0W51GGUlbReRlyWRCwjSgNFoTBgGYB02gPf98Z8QRhHOwMWNPkkc0e20wXjquqbf6zMYjLh9a4sgcNRYbu7sEkUhYRwRBZbV1TXSVpu6rpouDlFEkqbUO/vUBoIwoKoqyqokDMPmYTdjscYQRGEzOYeFts2ofU5RlRS3t0iMoRyVtKII4y22aAK4qzMiV2AGGSvOUWRDwjCkGu+RlyVxHGOHGUlZUO/uUuR5M3OdxgoWERERmZuZ9Ql+6qkrDEcjyqpifX2TD33wT6i9p/aGrCxoJwnGGIajMXuDEaauSJK0ueubT6hvV2yurRJEIXEYkOcZWZYxHI/Jy4Ld/X3SJOXmzjbdbpvQWd7wzGXW1tcJwoj93W1cHBGlKfk4w1pD6AJcEBI5Q10216uqCmOgKEsA2p0Ok0mGtZ5WFFMGljoNsdaSF3kzYxweAoNzlrIs6bUTirLCBo7UBxhrqcqYsrQkSYK1liRMKcsS40ua2eVm9kyiiIiIiNzHTJJYGAT0koA6g7jVwVQ5gQHqmvFoyHg8xq27ZviyvGB7Z5d+u8WF9U3GWcbO/g6jbMTNvW0+LL9AEnUJw5CsKIijkCgMSaKIIAipfM1wMiLqpHhgOBrhjMEFjjIvSNoxw3zCeDym1++zP5qwurlOXVXNhG4uoPL+tambbRwRuQAPVMZSlDXjyQRjDcZaLAbvDVme0+21KGoYDidgDbZs2rTOMshLoihmVNR4X1HXnqqum57GVYV1mixDREREZF5mEoIDZ7m8vkJoPC4I2dsf0koi8rzEWoc1hp29XS6srdJttZoH6aqafqfNhbU+YVCzMwJvPMNsRCuKKKqSSTYhCEPyg1nc6qpkbaWH2fO004RWmlDkGUUcESUpxXBEXVWUZXPdLMtpHYwHnGUZL1y9zn/4j78HzmKshYMZkp11lFVFWVbUHvKyPJgFucY5h/eeqqoJw5CyLKkPujYYa/C+uROOgapuRgP2NXhq6oN9Ho8xlv/1s2fx7YuIiIjI/czmb/LGYKiJQtfcUa1KrDF4asLAka6tUpU5Kyt9rPfUeYHxNUWRc3GjT7f7Bq7dvoGhYnt7h5aLcMZQFiV7+0PANCE0cBg8Vy5sUNVV04e4221q8BAGjvFoTJq2MCaj9gbvPYPhCPD0+j12BiOGeQ7G4g8muMCDN+C9xxpHDRhjaGY6LsE3n5Gs5GBn00XCM025GGuo6wqsazbVd0yP4Wswmi5DREREZF5mEoLLsqAsC9ppwtb2LnVdMx6NqKoaghBjoB0nDPf2cQY6aUI7beOrqnkwLXRcWd+klUbURd5MgBElRFHEOMtZW+1TFDlJFBFGAXESkqYJ3U6Hsq6JnWNyMDVz7WucsZRlxWA4Zv3iBYyvyCY5vU6HbqdNbFd429s/hjhJGQxH9DptxpMJSZriPVR1TZpEDAZjet0OO7t7dLsd8rzAWYMLAsbjMd1Om529fbqdFqPhmCgKwRjGk5xOK2Fnd5+VXofBcEQYRbP46kVERETkFGY2TlcQGHa3t8mzjPFwyFq/z5ve8CxvfPoZ3vGn38Zar48F0jghTRPyPCdtpXgseVFhMFR5RWhDjIciz7HWkCYpURTR7XRIkhhnLVVRMBqM+OAHnmf75haDwYi8KHBhxGAwJM+bB9qcs5RFQV3X1HWNwbO5tkqr3cWHXbrrlyHo0tu4DGGXdv8CrrWCD9ok3XVyH9Fe2aD0IVF7BR8kVDYhSLqMS0vcXWNSOKL2CrVLIEhxSZfSRHRX1skqS6e/RmVjwnZvVl+9iIiIiNzHTEKwOegeYCy00phuK+HihU2efuoyH/mOj2Cwt8twsEeRTdjf3+f27W3yqiYrK2o8Vd2EYA+UVQXGkOUZ/mDcX1/XlHlBWVZUZU1Vesqiotfp0kpb3Lq9xXgyoS4r+v0+1oV4bwjCEGMNzjmstYRhyOpKH++nD69VVL6i9v5gGLWKqiwP9tVNXQd3l6GmKkuom3XwTcW+whqaOn2zzVcVRVlhvKc6ONdriDQRERGRuZlJd4i3bD7HP/9r//7kAz51Fld9SJ/4EOd8yj32fdpD1iEiIiIij42mLRMRERGRpaMQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSUQgWERERkaWjECwiIiIiS0chWERERESWjkKwiIiIiCwdhWARERERWToKwSIiIiKydBSCRURERGTpKASLiIiIyNJRCBYRERGRpaMQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSUQgWERERkaWjECwiIiIiS0chWERERESWjkKwiIiIiCwdhWARERERWToKwSIiIiKydBSCRURERGTpKASLiIiIyNJRCBYRERGRpaMQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0gnkXICIicp699BL86q8+ejvGqI3z2Majnr+xAZ/0SY/WxrJSCBYREZmhv/AX4Pd/f95VyHn2278NH/VR867iyaMQLCIiMkMvvNAsP+VTwLmHb8f7s6nnUZ1FHYvwWRblczxKG9O/MPzdvws//dOPXsuyUQgWERGZEe9hOGzWv/VbHy0Eixz1Mz8D3/u9cOHCvCt5MunBOBERkRmZTJogHIYKwHL2dnaapULww1EIFhERmZHpXeAkmW8dcj5NQ/Dm5lzLeGKpO4SIiMiMTEMwwK/8SjMSwJ0ve3Arytp7b7vfOXfue9BzTjrvfudMXzI/29vNUneCH45CsIiIyIxUVbPc34f/8X+cby2z8jAB+mHC/1n/4vC4rzeLX1Q++MHmHIXgh6MQLCIiMiPPPgvf8i3w/vc3fYPr+vXL47bdbznrc05zzNSd6zIfTz017wqeTArBIiIiM2IMfNu3zbuK2fD+8LWoYX5WvwAs0jlvfWvzkgenECwiIiIP7M4+wRr5Qp5EGh1CRERERJaOQrCIiIiILB11hxAREXnMfv7n4Tu+A8ry9fuOm0b3pKl1F2m7almsGv/xP4ZP+qTj90lDIVhEROQx+77vg1/+5XlXIefZJ38y/MZvwMd+7LwrWVwKwSIiIo/Z9O7d13wN/Ok/ffrzTpqc4kG3P0j7s2z7Qa95Fm0/aPtP4nf+dV8He3savu5+FIJFREQes+nECG94A3zkR863Fjl/1tebENxqzbuSxaYH40RERB6zaQg+qT+nyKPIsmaZpvOtY9EpBIuIiDxm0xCsP1fLLOR5s1QIvjd1hxAREXnMpiH4Z38W3vOe5v20v+d0Eoo730+PP8t9R18nbX/QfdPP9zD7HtdnXYR99/suHsVk0iwVgu9NIVhEROQxW19vlu95z3zrkMX1KKF+MGjeKwTfm0KwiIjIY/b3/h58/Mc3f7b2vukW4f3xr1nsm8c19Rnvve+oR+0q82f/rELw/SgEi4iIPGabm/BVXzXvKmTRnGWwXls7m64V55lCsIiIiMgCOKs+wXI6CsEiIiJztr8P29sn9/s8btvjOlbkvFIIFhERmaMXX4TnnoPRaN6VnOxxhu5ZH7vItc3qc3ze58F/8V88+L/7eacQLCIiMkd/9EdNAA5D6PcPt9/Zv/NRt925vN+2O9s52p48mX7+5xWCj6MQLCIiMkfTUQDe8Q74zu+cby1HzSuIH7ftQc45i3of9TMsQr15Dt/1XYfjBsvdFIJFRETmqKqapV3AOVzVL/jJlmVNCNbMhMdbwP/kRERElsc0oCxiCJYn2/QXGIXg4+lOsIiIyBwpBC+304z9+7DbiqK5hkLw8RSCRURE5mjaHWIwgN///cPAcmegmde2eV9/3jU9jms9Dm95y+O5zpNGIVhERGQB/O7vwtd//byrkHmY9r22tnlN189im3PNL1fyegrBIiIic/Spn9qM43r9+tmGH217crbJfCgEi4iIzNHaGvzsz867CpHlo274IiIiIrJ0FIJFREREZOmoO4SIiMgC+bVfg6/6KhgO795+XN/Ro9uW5ZhFr2/Rjjn6/m/+TfiiL3p9O8tGIVhERGSB/MzPwPvfP+8q5Dz79/8efuqn4C/9pXlXMl8KwSIiIgtketfuS74EvvALD7d7//pjj2477pj7nXNW7Zym3Vm1o+/m9G187/fCiy/C1av3v/55pxAsIiKyQKYzx7XbcOnSfGuR8+ftb29CsHPzrmT+9GCciIjIApmGE011K7MwvTOsEKwQLCIislCmd4JP8+dykQc1nabbKgGqO4SIiMgimYaT69fvnu72Qfqm3uvYB+1DqnPP17kvv9wsdSdYIVhERGShBAf/y/zudzcvkVlI03lXMH8KwSIiIgvkL//lZgir/f27t99r7NcHHSdWbS1GW/OqY3MTPvdzWXoKwSIiIgvkrW+Ff/Nv5l2FyPmnbtEiIiIisnQUgkVERERk6ag7hIiIyILwHj7v8w4fiDPm8HXW72fR5uO4xpNa9yJ9N5/92fCf/+csPYVgERGRBVEU8HM/N+8q5Lz7zu+E/+f/gc/4jHlXMl8KwSIiIguiKJplFDVh2PvXjwF70rbp+tHlac+fRZsnnXeWbS7r57xfmydd40d/tJkw41d/VSFYIVhERGRBTENwGB6OFyxylvb34Sd+Qj9foAfjREREFkaeN0sFFJmVum6WmjFOIVhERGRh3HknWGQWqqpZKgQrBIuIiCyM6Z1gBRSZFd0JPqQ/uIiIiCyI6Z3g69fhz//5+dbyJLnzQTC5N90JPqQQLCIisiCeeQaefRaef/7wrrDIWet04OM+bt5VzJ9CsIiIyIJIU/jAByDL5l3Jk2c6EYTcXxDo4UtQCBYREVko1jZhWERmSyFYRETkHKlr2Nl5/XS6x72svfd+kfNMIVhEROQc+czPhH//78+uvdOE5dMes2zHLVptn/d58NEffXY/G086hWAREZFz5Nd+rVm2283y6LS5x72fDpt1dMrd6bbpiALyZHvnO+Fnf1ZBeEohWERE5ByxBzMA/PN/Dkny8O0cDct1fe8gfdzrLI49GtBPOvYs61v0az5MbT/2Y3D1Krz3vQrBUwrBIiIi58i0L++jjp17tF+wxpV9sm1twb/8l4chWjRjnIiIyLkyvROsCSTkTtNfaBSCDykEi4iInCPTEKywI3fSz8XrqTuEiIjIOTINO+9+9+F4w3d2bbhzuYzb7uzicZ63HaU7wa+nECwiInKOtFrNOMHf//3zrkTm6WgoPvpgnSgEi4iInCs/+IPNA1D3G63gtNse9fxF2rZo9ZzlZznquO3tNvy5P/f6Y5eVQrCIiMg58vmf37xk+dwvRFsLgZLfa/RViIiIiJwDR/tBy71pdAgRERERWTq6EywiInIOvetd8A/+QfMg1EmjJmj9bNfnff3Trn/VV8HKCktPIVhEROQc+v7vh/e8Z95VyCJ6/nn4e38Per15VzJfCsEiIiLn2Dd+I3z4hzfrd44WcPShqUXYfnT/k7r9NJ91Xt/NT/4k/MAPwDd/s0KwQrCIiMg5NJ0045ln4Lnn5luLLI5f/EXY2tIoEaAH40RERM6laQg+bgxZWV5V1SwVghWCRUREzqVpCNYMYXInheBD+gpERETOoWkI/oVfgN/93Wb9zvFjTxpL9rixZo879qT9pz32Qa5/lm3N4rM8yLHzrjXPm2UYHl/HMlEIFhEROYf6/Wb5b//tfOuQxWOMQjAoBIuIiJxL3/Ed8Gf/LJTl8SMH3Ll+v/0Pcqz2z3f/aY79z/4ziGOWnkKwiIjIOfSmN8G3fMu8qxBZXHowTkRERESWjkKwiIiIiCwddYcQERE5p3791+EP//Dhzz9pJAOd++Sf+7mfC93uw1/3PFAIFhEROYdu3IBP/uTDcWFF7vR3/y58wzcsdxBWCBYRETmHrl5tAnCvB5/wCQ/XxqPMNjevmeoe9rrL9Fl/6ZfgW78V/upfVQgWERGRc2Zvr1m+8Y3wzd8811JkwXzgA/DSS5pNUA/GiYiInEPTEJym861DFs+0H7FCsIiIiJw70xDcbs+3Dlk8CsENdYcQERE5h6Yh+Ld+C77+65t1Yw5fs3x/3D5rZ3/dx/3+Sf1cWda8VwgWERGRc+fy5Wa5s9O8RI56lKHZzgOFYBERkXPoC74A3v9+2N5uRhGY16uu53t9fb7jXx/+4fDcc/P+KZ0vhWAREZFzatlDjsi96ME4EREREVk6CsEiIiIisnTUHUJEROSc8h5+//ebkSKmfUGn26fLk9YfZZvaWex2vuM7EBSCRUREzq0f/3H48i+fdxWyaOoavuu75l3F/CkEi4iInFMvvdQs19fh0qVm/ei4sVNHtx3d9zDnHD33Qds7i3Ye9tqnaW+Wn39Wdf/v/zt893fDd36nhkhTCBYRETmnkqRZfvqnw3/73863FlkMP/IjzZ3gugbn5l3NfOnBOBERkXMqjptlns+3Dlkc07u/yz5bHCgEi4iInFvTO8EKwTJlD5KfQrC6Q4iIiJxb0zvBW1vwgQ8067Pqj/qo5zxqn+RZ1nye6E7wIYVgERGRcypNm+Vv/mbzkodzlg8LzvuXielfBRSCFYJFRETOrU/7NPjzfx6uXTubcWcfZP/jvNZZ7z/q6LFPuje/+bCrzDJTCBYRETmn1tbg3e+edxVPpvMc+K9c0cgQoBAsIiIi8jrnvW+waHQIEREREVlCuhMsIiIid/EefuEX4Pr15v1xD2cd3XbS+lkeuwxtPY7r9nqHI4csM4VgERERucsv/zL8hb8w7ypkVn70R+HLvgyiaN6VzJdCsIiIiNzlT/6kWT79NLztbXePinDcg1Z3rh8dZeFhjz06EsO9rvso9ZxVWw9T7zw+2+4u/I2/ARcvwud+LktNIVhERETucvt2s/yET4Cv+7r51iJn65u/GX7jN+ZdxWLQg3EiIiJyl2kI7vXmW4fILCkEi4iIyF22tpqlQrCcZ+oOISIiIneZ3gn2Hm7cON2IA6cZteC4KYQf5Bg5O/pOFYJFRETkiO3tZvm939u8FsVpgvTR484qgD/qLwDzvPadbf+n/4QcUAgWERGRu3zpl8Lzz0NVHY4qcNxUvMe9znr/nY7bJg/nypV5VzB/CsEiIiJyl6/5mua1KE4bou+1T/sP91+5Au94x2z+rZ4kCsEiIiKy0NQvWGZBo0OIiIiIyNLRnWARERG5y//1f8H3fd/r774efXjrfuunOfekbcY0f7q/13nHHXO/7fdy2nZO2+bR80/T1knHnPbco+tHj0mSZsKMt7/95LaWhUKwiIiIvOaFF+Cv/tXmoTg5n27cgF/8xXlXMX8KwSIiIvKa973vMAC/6U2v33+v0RlO2nfn9tMc86jn3uv9WbRxvzaPu8ZJx9+vtoc576R9ed68n06GsuwUgkVEROQ1V682S2vhb/7N+dYiZ+sHf7AZJziO513JYtCDcSIiIvKa69ebpXPzrUPOXlE0S4XghkKwiIiIvEYh+PyadnNRCG6oO4SIiIi85ubNZuk9/MZvPFpf2Hv1cX2Uc+/X1nmq7SxruXWrWSoENxSCRURE5DWjUbPMMnjXu+Zbi8zGxsa8K1gM6g4hIiIir/mmb4KP+Ih5VyGz9Nf/+rwrWAy6EywiIiKv+eRPhl/6Jfg3/+buP+kf/fP+Sdse9Hi1+3jbfeYZ+LRPQ1AIFhERkSM2NuCv/bV5VyEyW+oOISIiIiJLR3eCRURE5C57e7CzA3Xd/Bl9urxz/aTloh2zqHXNq/bVVfjhH266RSw7hWARERF5zWgEf+fvwA/8wLwrkVl597vhq7963lXMn0KwiIiIvOaHfugwAF+8CMbc/2UPOlfeuW7t7I691znHnf8wx56mrtMc+zCf+6w+w9HXj/0Y/Mf/qHGCpxSCRURE5DVp2iw/7/PgG79xvrXI2Zr+27bb861jUejBOBEREXlNq9UsJ5P51iFnbzxulgrBDYVgERERec00BOf5fOuQszf9xWb6b7zsFIJFRETkNdM/metO8PmTZc1Sd4Ib6hMsIiIir1m0O8FHh/o6bvive22/c1tVzb+9edb80kvNd6o7wQ2FYBEREXnN9E7wH/8xfNM3zT/sydlKErh8ed5VLAaFYBEREXnNh31YM8TWYADvec+8q2lYe/fLuSdv26LU85a3wMrKvP9FF4NCsIiIiLzmqafg934PPvShxQhw07FvRc6aQrCIiIjc5W1va14i55lCsIiIiLzOzg787M8evp/OjHav9dMed57WF6WOB1m3Fv7Mn4EwZKkpBIuIiMhdtrfh274Nvvu7512JzMqXfin8xE/Mu4r5UggWERGRu7zyymEA/qzPunuf9/dfP83xD3vOWdexbMcPBs1QadeuHX/eMlEIFhERkWO98Y3wLd8y7yrkLP3yL8P/9D/B2tq8K5k/PXMpIiIidznp7qI8+fb2muXq6nzrWAQKwSIiIiJLYn+/WepOsEKwiIiIyNJQCD6kPsEiIiJyrFdega/+6sP3R4cGO7ptun7a4+537GmPu9+1Z9Hm47r2WbUzXf7hHzZLhWCFYBERETniyhVIUxiP4QMfmHc1MgvPPjvvCuZPIVhERETusr4OL7wAL7/cvJ8+KHe/YbpOu+1hzlE7Z3fO5iZ85mey9BSCRURE5HU2N5uXyHmlB+NEREREZOnoTrCIiIic6Kd/Gn73d5t1Y17/oNVxy0c95nFcY9nr6Pfh7W8/+WG+ZaAQLCIiIsf66Z+GL/qieVchs/ITPwFf+qXzrmJ+FIJFRETkWD/9083y0iX4rM863H70Yat7PYh1Fvum7x9239FjHnTfo3yORfxsV6/CrVvNw4/LTCFYREREjmUPnhz6r/4r+JzPmW8tcnZ++IfhXe+adxXzpwfjRERE5FjTEFzX861DZBYUgkVERORY0xB89M/pIueBQrCIiIgcS3eC5TxTn2ARERE51jQE/9N/Cv/qXx0Os3XnC45/Pz33Xsfer52zaOOsajnLdk772R+kjQdp55VXTv43XyYKwSIiInKst7ylWd640bzkfAmWPAUu+ccXERGRk3zDNzSjQoxGh0Nted90j7jz/XGvZT5mEWs6ekyvB1/8xfP+CZsvhWARERE5ljHwtrfNuwqR2dCDcSIiIiKydHQnWEREZAlMJpBlzfqdD0o96PJhzxVZNArBIiIi51xVNf17f+iH5l3JowXwk5azaPO81+IcfOM3wpd8CUtLIVhEROSc+5//58MA3Go1y+kEGEcnwji6/UH2n3TOvc6X+fn+71cIFhERkXOs222Wf/kvw9d//eO//qME68e1f5FrO+3+0577R38E3/M96qqiECwiInLO9XrNcjSaz/XVP3ixDAbN0rn51jFvGh1CRETknJveCZ5XCJbFMp0GWyFYREREzrVpCB6P51uHLAaF4Ia6Q4iIiJxz0+4QL74IP/7jzfrRvqTH9R896ZiTzjupjXm3fb92H/Xaj6vts7p2UTTbFIJFRETkXLtwoVleuwY/8iPzrUUWx8d93LwrmC+FYBERkXPuLW+Bd74T/viPXz+O7HGve+1/lHNnvV+1nX5/GML6+qP9XD3pFIJFRESWwFd91bwrEFksejBORERERJaO7gSLiIgsia2t5gXH/5n8pPXTHjfvdZEHoRAsIiKyBG7dgm/4Bvin/3TelczeooTyRV+/cqUZLWQ6esiyUQgWERFZAt///YcB+MqVew8j9qRtP+q4KYTl9X7rt+BXfgX+4l+cdyXzoRAsIiKyBMqyWX7lV8J/+V/Ot5ZZOW7cXW0/fvv3fR/83u9BVbG0FIJFRESWwDQEn+cJEu78k7/cW6fTLJf5jrlGhxAREVkC0xBs9b/8wuEvCwrBIiIicq4tw51geXDLHILVHUJERGQJTPt+DgbwyivN+nHDix3tTjDPY47r2rBoxzypdCdYIVhERGQpTEPwP/knzUvOxqOE6ON+IXhcxwyHzbpCsIiIiJxrf+WvwC//MozHJ48ccNzyXvvO6pgn4RonOe1xi6jVgre9bd5VzI9CsIiIyBL49E+H971v3lWcD4sa1B/0Gp1OE4SXlUKwiIiIyAM4j32El5FCsIiIyBL7rd9qlif1Xz1p2t1ZHzOPaz6uumQxKASLiIgsqfe+Fz7mY+ZdxfJatHB+5/LLvgy++7sf7fMtOoVgERGRJTUNwG9+cxN8HuZBsTun4j3NOY/yUNpJ0wEfPWfR6jnJcbUuine+UyFYREREzqG6Plz/kR+ZXx3LZN6/AJzmlwPv4Yu/GPb3m2H1zvPkKgrBIiIiS2gagjWN8uPzpPQNbrWacYT39mB1dd7VzI5+9EVERJaQQrCcpNNpljs7cy1j5vSjLyIisoSmM8gt+l1JefyWJQSrO4SIiMgSmt4Jrmt417ua9aMjB9xvZIE73x+37UGOOe46j/P6T0KNj+v67XazVAgWERGRc8c5CEMoCvjhH553NbKItrfnXcFsKQSLiIgsoSRp7gD/+q8fjgpw3Atmu/9xXEM1Pvj+ixfh7W9/9J+zRaYQLCIisqS+6Iual8gy0oNxIiIiIrJ0FIJFREREZOmoO4SIiIjc5Xd+B/7ZPzt8f9JMY8e9f5Bj1faT1ba18DVfA5/2aZwLCsEiIiLymve+F77wC+Hll+ddiSyia9cUgkVEROQc+p3fOQzAX/VVh9uPTqpxr0k27nfs/SboODqu7WmPfdBrzevco/sf5dz7tfUg597r/Jdfhn/0jyDL7t3+k0QhWERERF4zmTTLL/gC+Ot/fb61yOJ4//ubEDydZOU80INxIiIi8pppCI6i+dYhi2V6h1ghWERERM4lhWA5jj1IjOcpBKs7hIiIiLxm2ufzt38bfuAHmjuA5+E1DXGPsm7t4R3Rh12fxetxOI93ghWCRURE5DWtVrN8//ublyy+x/FLRFk211IIFhERkXPp674O1tZgMGgCj/fN60HWH/d5y3rtqen2x+Ftb3s813kcFIJFRETkNb0e/Df/zbyrkNOYht/HFcAB3vzm+X7ms6QQLCIiIvIEurNPsHPzreVJpNEhRERERGTp6E6wiIiI3NP3fA987/cevj9phrLj1h/kWLX95LUdhvCN3wif/Mk8cRSCRURE5ET/8B82IUfkJEWhECwiIiLnzN5es/ysz4L/+r++e9/REQqOWz/u+Ac5dtbHq5b7H3/SsS+8AD/8w3D79vH7F51CsIiIiJxoOlHEhQtw6dJ8a5HFculSE4J3d+ddycPRg3EiIiJyomkIflzj0MqTo9Npljs7cy3joSkEi4iIyImmIfg8zRQmZ0MhWERERM4t3QmWkyRJ8/MxGkGez7uaB6c+wSIiInKiaQj+mZ+BX/zF44+5X0B+lAB9r3Nndd1HafdJ+6yPet3pXwh2dpp+408ShWARERE50Ud9VDMb2WTSvESOetvbYG1t3lU8OIVgEREROdFnfibcutX8yfte7pxI4UH26dzFOPdR2o7jw78YPEkUgkVEROSeVlaal8h58gTmdhERERGRR6MQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSUQgWERERkaWjECwiIiIiS0chWERERESWjkKwiIiIiCwdhWARERERWTrGe3+2DRpzO03Tteeee+5M2xURERERudMf/MEfMB6Pt7z36w967ixC8PNAD/jQmTYsIiIiInK3NwJ73vtnH/TEMw/BIiIiIiKLTn2CRURERGTpKASLiIiIyNJRCBYRERGRpaMQLCIiIiJLRyFYRERERJaOQrCIiIiILB2FYBERERFZOgrBIiIiIrJ0FIJFREREZOkoBIuIiIjI0lEIFhEREZGloxAsIiIiIktHIVhERERElo5CsIiIiIgsHYVgEREREVk6CsEiIiIisnQUgkVERERk6SgEi4iIiMjSUQgWERERkaXz/wfnoNKduDL4wgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 287,
       "width": 352
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(num='main', figsize=(6,6))\n",
    "\n",
    "# %%\n",
    "# Plot the image transforming the coordinates from pixel coordinates to EPSG6434 coordinates\n",
    "# data = rasterio.open('data/images/180700_se.tif')\n",
    "# show(data.read([1,2,3]).astype('float')/65536, transform=data.transform, ax=ax)\n",
    "data = rasterio.open('data/images/cropped.tif')\n",
    "show(data.read([1,2,3]).astype('float')/255, transform=data.transform, ax=ax)\n",
    "\n",
    "# %%\n",
    "# PLOT SEGMENTATION DATA\n",
    "raw_segmentation = gpd.read_file(\"data/annotations/segmentation.geojson\")\n",
    "segmentation = raw_segmentation.to_crs(epsg=6434).copy()\n",
    "segmentation = segmentation.assign(label=0) # Sets default label to 0\n",
    "# Convert Other, TP, FP, and FN to numerical values\n",
    "for (i, entry) in enumerate(segmentation.iterrows()):\n",
    "    if entry[1].default == 'TP':\n",
    "        segmentation.at[i,'label'] = 1\n",
    "    elif entry[1].default == 'FP':\n",
    "        segmentation.at[i,'label'] = 2\n",
    "    elif entry[1].default == 'FN':\n",
    "        segmentation.at[i,'label'] = 3\n",
    "# segmentation.plot(ax=ax, linewidth=1, ec=(0,0,1,1), fc=(0,0,1,0.25))\n",
    "\n",
    "# Need to feed in pairs of geometries and the corresponding values to use\n",
    "shapes = []\n",
    "for (i, entry) in enumerate(segmentation.iterrows()):\n",
    "    shapes.append((segmentation.at[i,'geometry'], segmentation.at[i,'label']))\n",
    "\n",
    "\n",
    "result = rasterize(shapes,out_shape=data.shape, transform=data.transform)\n",
    "masked_labels = np.ma.masked_where(result == 0, result)\n",
    "colors = [\n",
    "          (0.01,0.91,0.04),\n",
    "          (1,0,0),\n",
    "          (1,0.78,0.2)]\n",
    "cmap = mpl.colors.LinearSegmentedColormap.from_list(\n",
    "    'Custom cmap', colors, N=len(colors))\n",
    "show(masked_labels, transform=data.transform, ax=ax, alpha=0.4, cmap=cmap, interpolation='none')\n",
    "\n",
    "ax.set_xticks([])\n",
    "ax.set_yticks([])\n",
    "\n",
    "# Plot manual annotations\n",
    "raw_annotations = gpd.read_file(\"data/annotations/annotations.geojson\")\n",
    "annotations = raw_annotations.to_crs(epsg=6434).copy()\n",
    "annotations.plot(ax=ax, linewidth=1, edgecolor=(0,0,1,1), facecolor=(0,0,1,0.25))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "text/plain": [
       "<Figure size 432x288 with 0 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "9cc1658958d109e4d0e85844406ea2e7d6468316002f2d5e050fa2244cc4c3ff"
  },
  "kernelspec": {
   "display_name": "Python 3.9.6 64-bit ('ds': conda)",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.6"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {},
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}