{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Brexit" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import geopandas\n", "import topojson\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "tags": [ "remove-input" ] }, "outputs": [ { "data": { "text/markdown": [ "# Brexit\n", "## Brexit dataset\n", "\n", "This dataset contains results for the Brexit vote at the local authority district, and administrative boundaries.\n", "\n", "## `brexit_vote.csv`\n", "\n", "- `Source`: Electoral Commission\n", "- `URL`\n", "\n", "> [http://www.electoralcommission.org.uk/__data/assets/file/0014/212135/EU-referendum-result-data.csv](http://www.electoralcommission.org.uk/__data/assets/file/0014/212135/EU-referendum-result-data.csv)\n", "\n", "- `Processing`: no processing was required for this dataset, see original source for additional information\n", "\n", "## Local Authority District boundaries\n", "\n", "- `Source`: Office for National Statistics\n", "- `URL`\n", "\n", "> [https://data.gov.uk/dataset/1261032e-c683-4408-a82d-37281d0e250d/local-authority-districts-december-2016-full-clipped-boundaries-in-great-britain](https://data.gov.uk/dataset/1261032e-c683-4408-a82d-37281d0e250d/local-authority-districts-december-2016-full-clipped-boundaries-in-great-britain)\n", "\n", "- `Processing`: no processing was required for this dataset, see original source for additional information\n" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from IPython.display import display_markdown\n", "\n", "display_markdown(open(\"README.md\").read(), raw=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Local Authority District (LAD) geometries\n", "\n", "The original geometries are downloaded from the Office of National Statistics through [`data.gov.uk`](https://data.gov.uk/dataset/1261032e-c683-4408-a82d-37281d0e250d/local-authority-districts-december-2016-full-clipped-boundaries-in-great-britain) and stored locally:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "lads = geopandas.read_file('./Local_Authority_Districts_December_2016_'\n", " 'Generalised_Clipped_Boundaries_in_the_UK_WGS84/'\n", " 'Local_Authority_Districts_December_2016_'\n", " 'Generalised_Clipped_Boundaries_in_the_UK_WGS84.shp'\n", " )" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Topology simplification\n", "\n", "To obtain a much more lightweight set of geometries, we simplify the layer, respecting its topology. First we build the topology:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "topology = topojson.Topology(lads)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Then simplify it:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "simple_uk = topology.toposimplify(.025)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Inspect the resulting geography:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPAAAAD4CAYAAADSBmv/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAABv3UlEQVR4nO2ddXhc2Xn/P+feucMsGDHaMjMtY7LZ4IaxwTbbJk3KTZqmDGnaNCn/kt1wGmpos4HNZiFZBtvrNaPIYhwN89zz++NKY8mSbHktg+z7eR4/Hrz3zGjee8554fsKKSUmJiZLE+VSD8DExOSlYxqwickSxjRgE5MljGnAJiZLGNOATUyWMJaLebLy8nLZ1NR0MU9pYrLkeeGFF8aklBVzPXdRDbipqYndu3dfzFOamCx5hBAn53vOXEKbmCxhTAM2MVnCmAZsYrKEMQ3YxGQJYxqwickSxjRgE5MljGnAJiZLGNOATUyWMBc1kcPEZDHpj6QZiKSxW1TW1fku9XAuCaYBmyxJjg3FedPnn6GoSz5+5wrW1noRQlzqYV10TAM2WTIksgUeOjREyGtHl5JEtsD1y8q4c23VVWm8YBqwyRJBSsnHf7ifXx8dQVMVPv3Gdfz5K1dyw/JyqnyOSz28S4bpxDK5LHnyxCiPHx/lyGAMgH996BhHB2P4HRqfefN6NjX4yBV1fn1slJFY5hKP9tJhzsAmlxXjiSz/+tBxvrOzB4CffuQGACyKQsdoElURaBaFKp+TlVVe8kWdvb0R7lhTdSmHfckwDdjksiCdK/KroyN8/If7uXNtFXUBB83lLprKHIwlMrx9ez1ffboLr0OjpdwFwMtWhy7xqC89pgGbXBbYLIJXr6/mprZyDvZH+fn+QZLZAk+1j/MH/7eXe969hac/dgs5Hcrctlnvl1JelY4s04BNLgsUxXDHeOwa1X4Hb9xUQzJX5I++t49yt41fHBjikSPDvGZ9DWVuG+lcESHArqmEkzkUAX6n9RJ/iouP6cQyuaRkC0Xu39s/4zEFuH1ViNZKN1U+O/2RNL86NkK1z0Glx8Z3dp7k8493YBECKSVBl/WqNF4wZ2CTS4zNoqKpCh/9zovU+u08emSEWDqP226h2uegaywJwEQyR+dokpDXzju2N17iUV8+mDOwySXHrikc6Ivw5PFRtjUHyRR0irrErqm8bFUlPofGZ9+6gfdc24DLZs450zG/DZOLQr6oo6lzzxe3rQxx28oQg5E0P9s/iMuq8o0P7KDMpfGX9x8insnzxSc7uWN1FQ8dHuaPXr4CVbn6HFZzsaAZWAjhF0L8QAhxVAhxRAhxrRAiKIR4WAhxYvL/wIUerMnSJJ7J85O9A2d9XbXfwQdvauHJj99GQ5mT3okk25sCWFSF65eV87oN1dgsKkNXceLG6Sx0Cf0fwINSypXABuAI8GfAo1LK5cCjk/dNTGbhsWu8aUvdgl8/NbuOxHP8/OAghaLOpno/TeVufu/25dT6r97UydM5qwELIbzATcCXAaSUOSllBLgL+Prky74OvP7CDNHkasWqKjzTEaa53MU1LWWXejiXJQuZgVuAUeCrQogXhRBfEkK4gJCUchBg8v/Kud4shLhbCLFbCLF7dHR00QZucuUzksiiCPjq+7bjNp1Xc7IQA7YAm4HPSyk3AUnOYbkspbxXSrlVSrm1omLO7hAmJrPIFXQeOzrKlsYA9UEHlnkcYFc7C7ms9QF9UsrnJ+//AMOAh4UQ1VLKQSFENTByoQZpcvVhtSh89q0b6BhNXJUpkgvlrJc1KeUQ0CuEWDH50O3AYeAnwHsnH3svcP8FGaHJVYsQgmWVnks9jMuahW4sPgp8SwhhBTqB92MY//eEEL8J9ABvuTBDNDExmY8FGbCUci+wdY6nbl/U0ZiYmJwTpmfAxGQJYxqwyWVBJl+81ENYkpgGbHJZ0D6SuNRDWJKYBmxyWbC29uoUZj9fTAM2WTBFXVLU5bzPSyl5un3sIo7IxMxPMzkj4WSO/9vVy7amAKPxLNe2ls2rfpHJ6zit6kUe4dWNOQObnJGAU6NQ1Pnhnn5aKtx47NqM53MFnb++/yAv9kwwkcpR4bEh5fyztMniYs7AJmdECMFHb18OgK5LdnaN47CpDEUy7Ggp44EDQ6yr8xNJ5dlQ50cxC+0vKqYBm8xC1yWKIohn8vzGl3fisqq8Yk0VN7VVcE1rOX9+3wGqvXZubKvg2tYgNX4HNou5dL4UmAa8CMylSZwr6OSL+pLScGofifOdnb38ZN8ALqvKGzbV8qk3rEURgmgqj82iUNQlf3rHCjSLAhK+8Fgn2UKRT7xqFSGv/VJ/hKsOcTH3K1u3bpW7d+++aOe7GlhMQfNPPXCEe5/onPGYRRH81zs2sbHez+PHR8nrkk31ftbW+pBSEk3neezYKLeuqMTn1OY5ssn5IIR4QUo5Vyqz6cRa6ixmqd3pQnFCwGvWVzEUy/DT/QMIoD7goLHMWTq332nlzrVVpvFeIpbO+s7kgvOmzXXcuLycXEHn/V/bhSoEJ0aSOG0aFkXwy4FhNFXwzh2NPHhwkH96w3p8Tg27Zu5/LxWmAZsA8MCBQf78vgOsq/Xxzu0N3Peha1lZ7ZthnP/56Al+sm+ARw4bLU5Mu730mAZsAsDtqyp5W1899zzeyZMnxvjgDS1IKXDaLNQGHDx6eJhfHhriW7+1w3RWXUaYBmwCGC1OPvHKVWiKQkGXvGtHAzU+OwhBbzjFn/3oAF9671bTeC8zTAM2mcGfvGLFjPupXIGmchefftM6qnym8V5umAZsUiKWyQNgsyjYLCrHhuKsqDI0qe7aWHvGQgaTS4NpwJchA5E039vdS1vIQ1vIQ2OZc96+QmeiqMs5ewjt7g5zoD9KmdtGudtKuduGIgR/df9BnukYB4wQks2iUB9w8r+/uYMqn93sR3QZYhrwZUa2UOSzDx3nh3v6So997X1b+cxDx7lzbRUfvmUZf/+zw1gUQWO5i9dtqEZIcNk1FAGxdIFwKkfXWIL/+XUHG+v9NJe7GI5lsKoKT5wYZVf3xFnHISW4bRY+/xubzaXzZYxpwC+RF05OcGQwxju2NyzazBRN5/mNLz3Pgf7ojMc/9sMDjMSzvHJtFd96/iRfe6YbgCqvHasq+OKTXbisKpsaAnz92W6mJ9e9cPLsxjofL18dmlV9ZHJ5YaZSniPHh+P8y4PHeOTIMAAvW1XJl967bVGOXSjq7OuL8F+/auexYzPb0NT47CiKoG8iXXrsVeuqePzYKMnchdGTUhXBtS1lfPO3dlyQ45ssjDOlUpoz8ALpm0jxuYeP8+MX+5nuy6kLOBfl+L3hFF95uotEOs/Rwfis5weis1tqPnBgaFHOPR9FXfLbN7dc0HOYnB+mAc9DIlvghZMT7O4Os6s7zJ6TEXJFfdbrHIukQFHUJT94oY94prAoxzsfqn12Xr+pljdtrjU7I1zmXNEGrOs6irJw762UkkS2wC8ODPGJ+w6cMWzSWObkT+5YwWvWV5/H+Iy622S2QH3Azjc+sJ0Pf2sP4WSObGH2xeJCoqmC126o4U2b67i2pcwszF8iXNEGvBDjzRd1vvp0F1U+Bzua/fzx9/Zz+8oQzWVO2keTc77n1eur+dM7VtAxmuDjP9zPu69pYl3duakqjsQzfOwH+9neHOSRQ8PYrSqRVJ5wMsdNbRVkCzpPnhjlYrgoQl4b//POzWxtCl74k5ksKle0AZ+NRLbAZx86xnOdYdpH4mxt9LOhPsBT7aPUB510j6conDYLKwI+fEsryVyBRLbAX712zUvqXVvpsXPvu7fy2YeO8WJfZIahPnx4+Hw/2oIQAl6+KsQ/vmEdFR7bRTmnyeKyoF+eEKIbiANFoCCl3CqE2Ah8AbADBeDDUsqdF2icF4Sn28f46tPdpfvPdk4wlshz4jSRcZtFodxtY3mli9eur2FlyIOqKqypeelaxgORNDV+B+/YXs93d/USTedf8rHOFVURvGpdNb97aysrq7wX7bwmi8+5TB23Simni/7+C/C3UspfCCFeNXn/lsUc3IVmR1OA/37nJr7yVBe9E2nimZnGG/LaWFPjY1dXmP5Imv5ImsODcTY1BmipcL/k847EMnzwG7soFGF5yH3RjNeqKrxhUy0fuqWVpnLXRTmnyYXlfJbQEpi6fPuAgfMfzoVnynEEsKc3gk1VWFHlZTiWJZM/5Ti6Y02IKq+d7+/uI1/UaalwcUtbBaoiON+U4IDTyhffvZXv7urhC493nd/BFoBDU3nbtnp+++YWqn2OC34+k4vHQg1YAg8JISRwj5TyXuAPgF8KIf4VQ5rnurneKIS4G7gboKGh4bwHfL5M967etjJELJMnV5SkcwWy7UXGEjkADvRFefebG3nk8DBjiRwCqA04+c0bms/r/Jl8Eauq4HNqvHZDLe+/voVDg1G+u7OXn+0fXNAxAk6Noi7x2DUS2cKcM7iqCCo9Nt6wqZbfvKGZMre5x70SWVAmlhCiRko5IISoBB7GaPj9ZuBxKeUPhRBvBe6WUr7sTMc510ys+ZLxF5N9vREePTLMe69r4v891sHukxNEUjm2NQao8jl4+/Z68kVJyGvDaTWud/0TKWrPIYFDSsmnHzzKYCTDh25ppdJj494nOhFC0DGaQFMFUhoNvk7ff58Nq0WhtcJN30SKfFFHl4YiJsDaWi8/++iN53Q8k8uPM2VinXMqpRDib4AE8JeAX0ophaGsFpVSntEjshADTueK9E2kWB7ykC0UL6jecKGoU9Al//bwcb7+bDeFosTvtOJ3avzzm9ZR5rTic1oJuOZuJXI2/uYnh3ji+Cirarz81g3NfOqBI+zqnmB5pbtkqCurPBR0SftIApdVpbnCxcH+2Dmfq7XChdehUe628a4dDVR67Ng0hbF4lsODMQ4NxHjZqhB3rq16SZ/F5NJxXqmUQggXoEgp45O37wD+DmPPezPwGHAbcGIxBuuwqtg0BSnlBTPerrFESYz8gYP9/Gz/IJm8Tshr4+/uWssNy8oXRc+5byJF51iSgWiadbW+Ukhq+ix7dOhU2mQyV3xJxgvQMS1mPRWGUgTo0tgDf+6tG0zjvQJZyK80BNw3KV9qAb4tpXxQCJEA/kMIYQEyTO5zF4OG4IXzkIaTOZrKXDx+fJS1tT5uWV7By1ZWct/eAV69rnrexl3TyRV0MvkiXseZK3XuvqmFWr+Dn+0f5NO/OLpYH+GMqIpgRcjDlsYAyVyBp9vH+NJ7tp010SSSyrGvL4oqBOtqfaZM7BLhrAYspewENszx+FPAlgsxqPNluqcZoGPUmPGKRZ1IOs8P9/QRTmT517esRxECRQjetaNxQcduH0nwW1/fxWA0wydfvYr3XNs072sP9sd4pmOcCo+N8WTuvD7T2djSGGBbUxC/05CAzRV13DYLf/Wa1bMuStFUnt6JFM91jrOvL8q+3gg94dSM1zQEnayr9bGuzkdTmZP+SIaT40m6x1NYFEEklWNVtbf0rz7oIJLKMxbPMprIMhrPMpbIMZbIMjZ5P5LK8/LVIe6+qYUav+kNXwyuyEys6cZ77xMdfOqBo9T6HaRyBSZShsd2bY2X5zrDCKFwx5r5l5bhZI6gy1rqgNAfSTORypMt6PzV/YfY2RXm029aP2c2Vo3fQSSdZzSeXfwPOcm2pgDLKj0c7I/yhcc7So//6StW8K4djSUnYK6g8+MX+3mmY7wk0i5h3nzvnnCKnnCKnx+Y7Rnf0RxkT0+EPT2Rcx7v157p5lvPn+T1G2v59JvWmyof58lVUQ88Es/wTPsoO7sj/OiFPmyaitdhYVtTkH9503osZ5Cr+eA3dmNVFf7qtasIeR3kizqv/I8naZ+2j/3kq1bxwZtakFLy7Z09vGtHI0cGY/x4bz+vWFPFfXv6+fGL/WQKRQq6XLT8ZkXAqmovhwaMfbPfqfGyVSFWVnm4trWMd395JwKwa2rp4lXpsZHKFUlkX3rVk92isLbWx+7zEAsA+OZv7uCG5eXndYyrgauuHjhX0BmJZ6jw2LBZVCo9dl63oQ5NUWgMOtnXF8FlVXnnjsZ5jbdzNMFEKs/v376cP/n+Pu789yf56G3LURXB2hovsXSekcmZ1WlTS+/Z1hAA4Ed7+vjik13c83gnG+v92DSF+HkYzVwIYYSf3rS5Do/dwkQqxzt3NPDroyO8+fPPks7PLvQfiWfZ1hRYkKzOdNbWeHFYVbrHUlhUgUUVbKjzsb8/+pIvSD8/MGAa8HlyRRqw1aLMKrTf1xfhD7+/H6dV5YM3NPHWbQ1UeObWehpLZPmPR0/w0KFh/uLVKxHARCrPP/z8MFaLQiavs6zSXTJgn0PjyGCMrz3dxd/etRaAl60K8eiRETrHkuztjVyQz6kqgsODMQ4PxrCqCuVuK/fvPXtCnCIEmxv8aKqClIYa5XRv+OnU+OxE0nkODpzykA/OITBwLihi8cQQrmauSAOei1g6z1u21NJU7uamZUEmkrl5Ddhjt3BtSxnHhuJ88seHsKoKv31TC48dG+XYsPFDn1pCX9daRrXPwd7eCCPxbKkVyY6WMmr8djrHkjQEnbxrRz33PtGFXVNZV+vjwUPnr6aRm1YznCvqc6p2zMXzXeFZj+1oDs75+BSL2f+o3G3l3962kRuXVyzaMa9WrhoDHk3k+NbzvVhVhaZ3buRlq+d3XNksKm/f3sDbttWzsytM93iSV6+rYUdTgPv2DbK3d4LesKFN9UzHOMmfHuIfXr+WdK7AWDxDucdOtlDkn964npF4huFYlidPjFIbcJDMFhbFeBeTlVWekqd+LgaiGXY0B2k/z/NUee2srPbwL29aT6XZ4WFRuGraizomZ5BcUecj39nLkTl0p05HCMGOljLetq0Bt93CUDzHT/cNoAqBZ5rXOZ0v8tmHjzOWyPFnPzpI11iCfFFSH3SyqtrLr4+OMBTNoAhYXumhpcLF5eR81VSllAM+hQCCTo3mchfran0Mxc5vyQyGislX3rvNNN5F5KqYgUfiGX7vuy+iTWop//HL21hdc+51sO/c0UCN3863nu+he/xU0f3x4QTHhxO8cHKCeKZAfyRNXcBB/0SagWgaVQg+ctsy/u5nh9lPlLs21rKsws1DF6lw/2z0hJM4NLXk9GoqczIQzRBO5QmnXlqpo6YKNtb7AUG2UETXJf/1jk2mVM8ic1UYcKXHzqvXVfPgwSFsFgX9POI4t6yopKXczbUtZfSEU1R6bWys9/Ohb+7hn9+0np1dYTpH4/SMp2gLufnHN6zF79D4wNd3I6URe73vxf7F+3DnSbnbSmuFm+e7wrSUu7BaFLrHkzP21y+FMpeNjpEk4VQOIeBr799uzrwXgKvCgAH++I42CkWdSCpL19jcWlcLpaHMyQduaGYklkFTFQIuKw/94U3oUvLZh45xYiRBhcfGW7fV43da2dUd5h9fv5YD/VEeOjxMrqDTN5EqJZVcbFxWhVXVPtL5IkcGY4wlwjg1haDLygs9E+cdp/Y5NBrLnOzrjbAi5CHo0tjVNc7NbabTarG5Kgw4ksrhtql87m0b+cnefl6/qe6c3t8+EudkOMWm+gBBl5Xd3WH8Tm2G5GrIa+cvfnyAVK7I3921hrduraegS377f3ezq2uCXFHHoghq/A6GYpnznuFeKte0BOkaSyIxyhxr/A7CyRypXPG8EzOmiKYNZRMdODYcZ2WVh/dff3511CZzc8VnYh3sj/IXPz7Ipno/v3VDM2UeK3btpV230rkif/z9vTxwYAhNFfzb2zbymvU1pecHImkqPDY0VSGTL/LpXxxlPJnj5rYKHj8+yjPtY8QzhTn1pS8WWxsDi2aoC2FDnY+vf2D7gopETObmqsvEmmIwmuYff36Em9sqePOWOmqD55Y4MJX/PIXDqvI/79zM0+3jPHJkmNtXhpBScmIkQVvIMyNBf2dXmG/v7CFX0PnpvgFCXtslMd6AU6O10k02X8ShqYSTeeoDDgq6JJrOkzqtLctUCeIUO5qDJdWP0XgWv1OjqcxoltY9nuJsfOqN60zjvYBc0QZc7rTysTtXsLra+5Lac0433umP3bC8vJQC+H+7enixZ4JPv+lUwdbx4Rjf293LLW0VHBqIMRBNMxy7cAUN09nS6KeoG9loAiOrbPc8aZNra7wUdFkqj6zw2BiJZynqRghMUwTPnZbcMRzLlj7L5gb/GQsa7JpSCt+ZXBiu2CV052gCn0Mj6LLOaYiLwd6eCTbU++c8/v6+CM92jPPLQ0P0hFOz4qxgzHZBlxW/00pj0EmZ28qjR0bOqfTQosCmhiC6lAzHMvRH0i/JCVXltSGEOOcUybaQm+PDs5NAblxezsfvXMna2pcuvWticFUuoZPZwnlJvy6EjZOFC+Fkjq6xJH0TKSKpPDe3VfDgwSH290U4NBCbt02KLmEskUNVBLesrGTPyYkFGW9zuQu3zYLTqpLMFtjVPX8K5EIZS2SpegmKlboOm+r9jMQz9EcybGkM8KevWME1LWXnPSaTs3PFGvC6Ov8FP4euSz7/eAf/9vDxGR0chAC/Q6OgS9pCHgYi6XkNU1MFr11fwy8ODNI+kqCpzIkiBD3h2V0hpnDZ1Fk9hM+Xgg61fseM9qULoX00gUURtIU83NRWyT+9cd2ijsvkzFyxBryYSCn55aEh0rkib9h8KgSlKILfvXUZb9xcy0gsy7HhOPt7JxiMZkjldCSSztEk17aWUeaycmIkwZ6eCTJ5nQ11Pt64uY5XrqviJ3sHGIpluH1VJYcGYhwZjM1pvNubg+QLOkeHz54G+lJ4vivMtqYg6cm2MT3h1II0sAu6ZDSR5SazNPCiYxrwGSjqkp/tH+BLT3bxrh313LHmVCfCgUgar0PjfV/Zye6TE/z2TS185NZlPN0+xqNHR/HYLKX63+l6z3UBB197/3aay10MRtP80wNHONAXw+uwoKkK0XSOfHG21bisKvFMniqv/YImsE8tx9dUe1ke8tA5mphzPKcjJdxkJmpcdEwDPg1dlxwbjpPOFSlzW2kLefjq+7dRPk0YvahLvv5MN198srPkZb3niU4OD8b4z7dvpNbv4MGDQ2QLOn6nVqobrvLa+fu71tIbTlHUdT770HFyRZ3xZJb2yWqgrY0BMnmdKq8dh6aSKeiMxjNGbnIyh8tqobnCzcFFXkKfTjSTx2ZRFmS8G+p8fOqN6xZFydPk3DC/8Uk+9/BxdnaNM5bI8clXraK10jWnOmYqV8BptfDy1ZUUdcnT7WMcG44jhOBAf5TX/NfTrKv1sbrGy7amIG1VbrrHkly/rIJtTX6+u6uXa1rK+fbzPUjg+c4w6XwRu6Zw64oKnu0YJ5IuzBl20lSFvb0TVPkcbGrwY1UVElmjeGJKKqfaZ6dQlLx4niICIa+doUmPdNCp0VjuIpMrcmIkzpRPbkXIwx/d0cYrzqApZnJhuWLDSOfKN57t5shgjD+5Y8W8bUim1C6nOkaEk1lyBcmzHWN8a2cPL/ZE5hSJ294c5N/fthGXzYLPobG7O8zTHWP0jKfI65JcXmcwmmZf3/yzqsdmYXWNl2NDMSLp2dI8FkVQ0CU7moMc6IuQyp9/wojPYeHuG1toLHcRzxToGktgs6hsawrislnYVO83q4suAldlGOlcOZM87BRTP9YpJcWgy0ZvOMVf/+QQugRVCIqcMmC3zcLH71zBseE45W4bvRMpvvncSf7r0RNkFpgL7bVbWF3t5eBA9IyKGQVdYrMo7OmZWNCydzpTF4fxZI72kQQrqzzctbEGr0Njd/cEX3/2ZGkbEHRZecWaKjO+e5lgGvA58pN9AwgkqVwRh9XC15/pxmpRCSez6BJaKly8ZUs9z3aOk84VCHntvHN7A+mCTtdokiODMcrcNvojCwvX6BKKEio8dsrd8ozpi7V+B70TZ09vPJ14tsDRIaPoYP2mWvb1RfjnB4/NeXy3zcKffH8fb9xcy+s31VI5jyyRycXBXEKfA1O50d967iT9kTR2TcWuKdg1FZtFYVml0RFh+munbv/3oyc4MBjj+FCc/kj6rLOkXVNQhJiRq1zjs59V98ppVVlR5UEAx4biJHOzlSnnY22tF01RFrx/VhXBXRtr+IfXry01fjNZfMwl9CIxZZDvuubsXRyEEOTyRbL5Ij2RFD0TKT72ipV8+FsvzJugAUYoptZvp2Mkyc7TMqzqg86zGnAqV+TFyfxkq0VhfZ0Pq6pwYiROdI6983S6xpKsqfHNKmgAw4OeK+qEJxNS1tX6eM+1jbx2Q82iCt6ZnBsLMmAhRDcQB4pAYepqIIT4KPARoAD8XEr5sQs0ziWJVVPpi6Q5Ppzgr1+zhocOD2GdlHI9HSHg+tZyHJrCd3b2AkZ7k55winK3jUgqRzSdR1ME+QV2GM8VdPZPOsa2NQXoGEkihDFzNpU5OT6SIJrOE3RaCbqspTDQ9MNvawrwiVetYvNk2mj7SIJ0rjij19LFaANrMjfnMgPfKqUcm7ojhLgVuAtYL6XMTvYOvqrIF412JbeurJwRJ57ihZNhPv9YJzevqMBpU9nSFOSvf3J4xmtURbC9KUjHaIKdXeEZ5YYWRfCmzXU8fnyEgi45OhR/yfW8LpuFcrfG8ZEkK6s85IoSh6aQzAjGk7k5Uz3/7W0beMNp4gfLKmfnl5vGe+k4nyX0h4BPSymzAFLKkcUZ0tJASslHvr2Hoi55y9b6GY8/emSELzxuNAvfUOfjLVvq6BpL8r6v7SLksxMfSSCEoRv13+/cxJ//6EDJyzud2oCD8WSWcretVM20++QE1zQH0YHBSJres+Quq4qg3G1lKJohmSvyxs01bKgLsL8vSjI7d7x5RcjDnWurZhnvFD94oY+usQROq9GeZntz8By+OZPFZKEGLIGHhBASuEdKeS/QBtwohPhHjPaifyKl3HWBxnnZIYTgt25oYlvzqfxfKSUf+uaeku5zldfOF9+zlWg6z9ee7qJQ0Et60n9wextN5U7+9Af76A2n0VSBzaKSK+jkijpW1eg/tLLKw/HhOCuqPDx4cIiN9X5emAwV7WgOMhrPkinoeO0WGoJOHFYVh1Xl5FgKl03l5FhyRg3vr46OUtQh4LTyxs21eO0aE6kcveE0FlVw+8pKrltWjl1TZ3V5HEtk+fgP9vPo0VPXarum8N7rmqhw23jZqhBN5ReuNazJbBZqwNdLKQcml8kPCyGOTr43AFwDbAO+J4Rokae5tYUQdzPZO7ihoWHxRn4ZMGW8Q9EMlR4bPzswSG3AwavWVfFiT4QvvmcrlV47+3sjtI+mmEifErG7YXkZ4WS+ZLwfuXU5//mrE/zWDc3c80QnuaLOmhovQ9EM/RNpmstd/N1daxiOZQi4rDxxfJTnu8I0lzuxWVQsiuDIYIyiNGZQzaJweDCOpgpWV3vpDSeJZ4tkckV6wim6x5IcGYqTK+gIYcSCY5kC33q+B4siWB7ysK7Wy4Z6P4cGYhwdjHF0KD5LwSOT17nn8U4A/vGBI7xsZSWv21hLXcBBrd9Bhcd2weqxTRZowFLKgcn/R4QQ9wHbgT7gR5MGu1MIoQPlwOhp770XuBeMMNIijv2SUCjq7OoKE3TbuH9vPz/bP4jXYeH61nLuecL4ITeWOfn7u9ayttaHlJKilNyyooL20QSZvDETdo4mqZ6sv71xeQUDkTSv31jDda1lfOPZk1zTEuTlq0K85r+fwiIEP57seaQpAs2i4HNorKv1UeG2kcwV6B5PYdNUllW68dgspU4L+aLk8GCMlnIXmxscxHMF9pyMzPhMFkWwqtpbShQp6JKxRJbu8RQ/2tOHXbPQWuHCoamzDBgMYYJtTUEksK8vysNHTs3QVlXhprYK/uddm7BZTG/1YnNWAxZCuABFShmfvH0H8HdAArgNeEwI0QZYgbH5j7T0+eWhIT72g/24bRa8kw3NwAjX3LjMqMQJeW08+oc3YZn8sQohWF9n6EaPxrM4rSrXtJTxhcc7+OJ7tvL1D2zn2Y5xPnrbMr7weAd3/+8LfOoNa3nNhhpsFpVoKs94IssNy8rJFgzj2dsbYTCaYTCawWu30FjmQpeSTQ3+0hLc77QyNG1/2zmWpHMsOSm2PpN80fAir6/zYbeoDMUy9IRTpb7GiWyBg/1R/E4r62q9TKTyWFUFTVXQVIHLZmF/f5T0HMadK+o8cmSYD31zD1/4jS1YLVdNM5CLwkJm4BBw3+QyyAJ8W0r5oBDCCnxFCHEQyAHvPX35fKXx4MEhouk8WxoDrK31kcoVODmeIlfQqQk4ePYTtzEWz5aMdwpFwN+8bg1Pt4/SWObi8491MJ7M8Vvf2M1PP3IDN7dVUNQl25uDfPXpbh45MsKbttTz030D3NxWwfdf6GNv7wSrqr20jyQIee2EPDYkhurm4YEoK6o8SN24n87rqJOzYvd4ktF4FiGgucyFx25hS4OfsWSWfEFSF3SiCjg+HCeWKbCx3k+1z07Iayu1IN3c4Kc3nCbg0jg0EFtQjfDpxNJ5/ur+g3zqDevM/OlFxMzEWiDxTJ4PfmM317SU8bu3LuP9X93FU+2nFhwuq8ov//CmWS0zu8eSNJW76BxN8Jf3H+Tp9nFsFqPb4W/f3IqqCNpHEqXc4kzeSMRoCDpI54t88BsvzBKibwsZKpORdJ62kIdUrsjRobjh4DotxGSzCFZWeekNpwmn5lYFsVkUVtd4SwkgU1zTHCSeLZQaiL8UHFaVFSE3hwdiqIrCn75iBe+/vsncF58DZibWIjCRyvE7N7dyywoj3P23d63hTZ9/hshkdwVdGkkOpxvw15/t5q6Ntbz1nmdLYu6/fXMrf/TyNo4MxviD7+7l2HCc61rL+Ne3bKDG7+DaVkNPKprOsyLkmWHAyyvdDMUy+B1WoulCaZYMODVcc6QzZgvyjFVOYMR2e+bIsT5dkfJccWgqTUEne3snz18s8nc/O8yDh4b4x9evRVGgtcJz5oOYnBFzBn4JSCn53+dO8mzHOGtrfayq9rCp3k/ANTuZI1fQ+dbzJ/nbn55K4PjUG9bxxs213Pgvvy7tMwH+/FUrufum1lnH2N0d5pmOcY4MxvjFwSE2N/jpGU8xNkfyxfQ+v0JwVoXKzQ1+XuyJsJi/Aosi2FDv48RIgtg86ZuqIqjx2/nGB3bQVOY0Z+QzYM7Ai4wQgvdc27SgEkSrReHEyCnZVa/dwg2TcdZPvmoVf/B/e0vPPdsxzhs21c5qPL61KcjWpiAnRuIcGYyxpyeC06qypsZLMlsoVSgJjID9pnp/SV425LVx8AxL4M7RJIowKp4Wi4IuaR9JUum2E0vP3Xe4qEt6w2lu/dfHSmWX717A92kyE9MleBFon2bAb95ST0OZscyu8tmxTfPK7j45wRPHR2e9f4rllR4+eFMLYBQtRNP5GcLpEqMjxIu9EUbiWcrdVo6dRQCvpcK1qMY7RTSdZziWYUXo7EvkRLbA5x/roHAJW84sVcwZ+AIzHMtwbMgwomqfnXdsP5V2eU1LGb/4/Rv571+1o6kKQbeVFVVz9y0u6pLnO8d5pn0cdVIVJOC0zisva7UoeB0WNtT5UBUFCXSMJOaVty13Wwl57VhVhY7RBLHMmSuXzobNorAs5Ca/QOGCgWiGR46McOdaU57nXDAN+ALywskw7//qrpIxjCdyswyjpcLN59628azHets9z5aKGDY3+BmKZhiJz19auLbGy/NdMz3StX47VV47QzHjfcsr3RR1idduYSyRYyyRY0tjACEEbpuFRPalGfGmBj8nhuKzvNpnoi7g4KY2U5b2XDGX0BeQF05OzDDYVTVevPZzu2YOxzL8wXdfnFGBdHgwRn3QiWWeeKrVopTqdqfTH8kgpaQ+4EAICLis7OuLzhjjRDJHUddJZgs0BJ0LWgLbLAorqzxGYUNTkONDcRLnICQARp71157pPqf3mJgz8EumUNT59s4eYuk8H7lt+ZyveevWeu5cU41mEWiqgt2i4LZr53Sehw4Pl9Iop8jkdZ7vCuOerN/VVEFTmYt8QSfgtpIv6PM6robjWQJOjWqvnZ1zhIk6p4WsesIpQl4bm+r9WFRBOlecdVyXTaXO7+To0PmJzWfy+jnN2CYGpgGfhYlkjl8fG+GRI8Ps7YlQ4bVTH3BQH3RSH3Cyod5P91iS2oBjVgdEv9N6Tq01cwV9RqphJl/k8WPzO7Vq/HaODyfY0hAoxWy7w2fXxJpI5Zkgf9bXwcxuhG0hN36nVop9u20Wavz2szrKFspjx0Z4rnPc7Kt0DpgGPMmf33eAQwMx6vwO6gIOnFYLT7WPsuc0qdiBaIZ9c2hG2TWFB3//pvMqp5tuvM90jPHJ+w7OysICI77bWOakxufAZbPMkt65UKhCEEnlUQRsrPdj01T29ixes/CCbrSwMQ144ZgGDBwaiPKdnT1IyZzGuRAyeZ2ecGpR6mGf6xznXV96fs4kjO1NAU6MJOgeS9E9du4KlOfDseE4N7eV0zGaLPUFbixzUtTlOTdFmwurqnB4IEY4mSPoMpuCL4SrzoDHEllODCeQUmJRjRnly091vaSeutNxaConRuK0VLhmpVOeK5/40QFqfMZKIJUr0j5idACs8NhI5XUmUgtb/p6JKq+d0UR2TiH601ld7cWmKcQzeR4/PrPg7OR4CptFocJtQ1F4SY3MtzcHsVsUktkCmbzOX/z4AJ9760ZTLG8BXDUGHE3luf1zjzOWOPUD29YU5NhQ7CXFPLc3BRiJZxlNZFlb4yNX0PnSk11c01JGXeDcxxdJ5fjC4528Zn0162t9PHBgsKQdPZVhFX+JYZ25GI4bSRZHh+JU++xUemycGEmQyRdRhCgpZzaVOckWihwenD+bq8prx2O3nDHjay6squDa1nIODkQRCNbUeBlLZNF1+PGL/bx9+5UlAHEhuCoMOF/U+cXBwRnGC0Ynvu3NwTm9sS6rSo3fgcduwaIoZAtFDvRHWRHycDKcYmf3qb3f9I4Jv/utPfzN69aUih4Wys6uMF94vIORWIbbVlVy/75Tnuf55sjllW4ay5w8cuTc5cikhA/d0sr6Oj+NQSeKIpBS0jWW5LnOcYZjGWr8Dv7fYx1kztKmpcpnP2PXiLkIODV8To2xRJZ3bGtgPJkjnMzSUu7ixEiCH+7Jck1rGU1lpkTPmbgqDDiTL/JnPzow53NHBmJsqPPNqNjRVEG13zEjh1kRhsF0jCVLVUVz0T2e4n1f3cUr1oT4q9euoda/sK73U/pXP3qxn18fm6lokZsnxbA/kuZdOxp4vitMPFPAqiqUua0MzqEdXeaysqzSzfHhOLesqGQgkmY0nqV5cs+ezhX52f4Bav0Ovrurl67RJJl8kRVVnrPOrJnCwmK+0/WmWyvdJDIFxuJZHj06jNeu4dBUusdS+Jwaui751dERPnB984KOfbVyxRtwNJ3nj7+3b97np5alWxsDSCS5giRbKHJ8eGYSvi7B57SSG547Of90fnlomCeOj/GR25bxwRtbzqpEUeN38KMPX8cf/d8+Hjw0hBBw55oq/uyVK1EVwfd29/Ht53tmrSJ+tn+QH37oOk4MJ9hY7+Ou/3kagGtagvSMp0pC8BUeG39/1xqayt1zjuVLT3by2YeP47VbZmwpusaStFa46BhNEnQaF4FcsUg0XaDMZSWdL+KxadgsgpYKN167Rm/41HndNgtraryEk4ZET0uZi6DTyq7ucGllMTytImt1tZdDA1HW1Ph46OAQb9tab7YtPQNXfDnh2+559ozLu6BTw+vQ8Ng1bJrC7u75wyK3rqggnSuiS5BIdB10KcnrOoWiNNISHRoCZmROzaWvPB9SSv7tkRP8aE8fv/j9G/FMS/xIZguGFM+TnTOWtf/8pnW8bZuxX5xI5rjniU6uaw3y/Rf6UQW8en0NN7dVnPEicvtnH6NjdHbICowUzPqgk329UdL5uWfbqe6IYITUNjcEjIqjiRQDkTN3kzgdl1UttYR53YYa/vMdm87p/VcaV3w5YSJb4FvPnURVBLeurKQnnOJgX5QD/dEZzpcKt42GMieRVI5YusBEKkc4lacgJYUFhEKODcXP2toEYH2dDyFgW2OQN2yu5ZVrqxf8WYQQ/NHL23BoKn/y/X3c8+5TfzeXzcIf37GCd+1o5DO/PMb9e/sp6HKGBz3gsvJnr1wJwE1tZ9+HZ/JF/vxHB+Y1XoBKr53nOs+8x53eLqa5zMUzHeNnPfd8TO/ntL8vwhcea+cd2xvxOc8ti+1q4IqYgff1Rnj3l58/ozd5c4OfEyMJ4vO8ZmWV56zpgHUBB1ZVmZFuOB0BLKt0YbWo/M87N1MfdNI1luTYUJxjQzHW1/l52erQgj7Tr44O8wff3cuX3rvtjMLpU/vxlyIWNxhN88Fv7OZg//x73G1NAY4Oxef93uai2md4pU/fhpwrW5sCHOyLkinovOfaRv72dWuuysL/M83AV4QBg/FjHE/kKOiGTvFEMkde10lmi4wncowmzhyf3FTvX1BXvumKFx6byvKQB01V0KVEl5IXTkYQwthTa6oyYyayKIL/fucm7lzgjJwr6OhSXrB46Fef6uJ7u3sJp3KMxrNzitW1lLuwWpRzznVeEfIsOMVydbUXp03lhZMTSAl+h8aqGg/PdoSpCzi4ua2CN2+p5YtPdjEWz3FtaxnXtZaxqSFwVahcXrFLaEONIkm524bHrvHgwSH+32MdM2RqFkr+LMXkQZeV1goXEsmGOh9+p8YzHeOljCSHplIfdBB0aiwPedjfF8VpVWfkDhd0yR99bx+tFW6WL6DK50L+OCOpHF96qqsUa9ZUwbaGAIcGoiSyp5awnWNJNjX4z/n4nWMJtjUF6JtIMRg989/DZVXZ1T2B36mRzBRQFMGRQcP46wIOXjg5Qd9EmscnxQ52dof5j0dP4NAMid6b28q5c201Vb6rr1fxkp2BI6kc//jzI3z/hb7zOs6GOh+Zgl4qup/r+ezk86d/U+tqfbMK6qd0qHwOCyuqPOzsmu0UW1bp5oHfu/GSzh6/ODDI7357z6xZN+iy0lruYvfJCSRGvHZjg5/9vdGSGICmirP2N55iRZWHY0NxNEWAoPS+kNdGU5mLcDJH70SKVdVeCkWdwwMxVlZ76Q2niGUKLKtwU+WzlfbY+YJOQTd8FoWixO/UeL4rjEURvGxViHdf28j1y66suuIragmdL+p84kcH+MFLMNzWChflbhtFKYmn8wzFskTT86clLq90c3I8SW6eH2ulx0aVz04qW6R99NR+b1tTgEMDMTL5Ihvq/dhUQxFjOJahezzFlsYAP/zQdec8/sXm0SPD/P53985ZuL+swjXpvBpHl4YIXZ3fwUgiC9JoBj4az5b0uE5nqrWKBMYTWXrCKSyqQlulh3S+wImRxKz0VbfNwsZ6H12jSQZjGWr9Dvoj6RkXGY9NZUWVlwP9UbLzxONbK1x87M6VvGLNlaHucUUZMBihlq8/080//eLovH/EuVhV5eHESKJ0NffYLdx9Ywuv21jDvU908t1dvRR1ybbGANFMnrF4bl4t5dPZ1OAnlspjUY1i+rn23FP752/91o7LZpY4PhznN7++i4FIZkZetNdhobnMdVZJ2pVVHjx2C9F0npF4trRdmO4rWChnap1aH3DQO5Fme1OAWKaAz6GRzhfpGk2WYvlumwVNFUyk8tg1hZ985AbaFrBVudy54gx4ihPDcX7/u3vPmKd7OnaLQkuFm2tayvi925fNqNdtH4nz74+c4Gf7B895LAv11r5zewOfeuO6cz7+hSSczPG5h4/xf7t6yRclzeUuUrm5W4+eDY/dwvJKN93jqTlVQebjTE7EqefW1niZSOXonxZXFhizfSyTJ5ou4HVYqHDb6BhNsiLk4f6PXL/kiyLOZMBL2oW3POThgd+/kaf/7DZet6FmQe/JFHQOD8ao8dtnFdsvq/TM2d9nOiurPKyt8dJa4SLktTGlaqMIMcN4HZrhwJpCVQQfuXUZf3vXmgV+uotH0GVlQ52/tD+tcFtnbC3s2sJ/JvFMgT09EfyOhcdsr2kOlpqxzcVU5OjgQGyG8YKRJ947kSY6qT+tKQojkxeeY8Nx/ur+g+gvpRfMEmFJe6GnqPU7uH1VJZF0nkQmz1giR88ZlCkCTo3nOsO8cXPdrLrTf3v7Rv7pgSM8cXys5KGt9NhoKXfRE07NCqesqvLgdWhkpmUoWS0KX37vVmyaykAkzWA0zbamIJsaXkKZ0kXi9ZtqebZznB/t6efV62vY1BDgK0934dBUbl8V4r4X+wG4fWUlXodWuu93akZfYk3FaVUnnVJGG5emMidumwWbRaEojaZsYFzcllW6cVhVBiNpXuiZYFNDgJ1dYXwOC9U+B91jSTKT26M9PRGay11zihucTtBlJZMvUuU3FDa/t7uPY0Nx/vEN60rta64kFrSEFkJ0A3GgCBSmT+dCiD8BPgNUSCnP2J3wYqZSfn93L5964AgTqTyqImgMOjkZTrGsws1fvmY1Nyw/+x60N5zin35xhAcPDp1TQ6+2kJuH/vDm8xj9paNzNMH3d/fxpac6sSgKmxr8fPDGFj73yHE+eusyesIpvvB4B/VBJ1ZVKXUlXMied0OdD7umsmeyQfl0/A4NRRGlZXeV12iwNrUHd2oKq2t8JHMFjg7FF1S/rQi4vrWM8WSe48NxPnLrMt6xo4GQd2mFmxYrDnzr6QYqhKgHXg70nMf4Lghv2VrPy1aF2N8fZUtjAJdV5RvPnqQu4OD6ZQuTbKkPOnm+MzzLeO/aWEO528bB/igv9kZmVSfVn2dB/6WkpcKN1SKo8tnZ3BDgU29Yh6oI7m+7HkURZAtFHj8+ypMnTv0UPDZ13hzp6ezri+KYzJN+4eTEjPTLyGnRgKFYBqf11N41lddnyOruOYsAXu2kNNKT7UYijcuq8u+PnuA/f3WC37imkT95xQq85ygweDlyvkvofwM+Bty/CGNZdAIuKze3VZTuv/e6pnN6f76ol/rp3tRWQVvIw4mROO/c3oBlUsAuky/ywskJ/utXJ4ik8vze7ct55RIXJ3/T5nq2NgZZX+efUQnUPhLnC493zir4iGeL7O+LcsuKcsYTebrHkvOKD6QnFTXnq8OezkA0zY7mIMeH40yk8rhsKmtqfPP6KcrdVur8DnrCKfoj6dIWCE7lV3/ng9ew4wrS3FqoAUvgISGEBO6RUt4rhHgd0C+l3Hel5qdqqsKX37dtxmNTnQOnsGsq1y8rZ0tjAJtFuSJydRvKnKX2L9MZjGbQdUmuqM+o7V1e6caiCh47ZszK041uPnIFvdRhYj6m5HOdVpVb2ioYjGZKRr+1MYAuJZqqEM8UGIkbHRv3niXsVR9cuqujuVioAV8vpRwQQlQCDwshjgKfBO442xuFEHcDdwM0NFy5EilLPVRxJuLpPO2jCSKpPIlcgU31fnafnGBHcxCHpvDYaTpZU3vhLY1+7JrKzq5wac+7ZjLv2WVVqQ845k0EcdtUllW46Yuk8dg0DvZHGEueuiAcGoiRLxaZvnsZS8wOW4W8NgpFWcoiu/eJTv7mdZdfJOClsiADllIOTP4/IoS4D7gZaAamZt86YI8QYruUcui0994L3AuGE2sRx25yEZBScvO/PsbmBj+PHBmZ0bI0lSvOaM42He+kFNHT7eMsq3CTKRSJZfLYNKXU0zjotLK9OUg4kUUoAkWI0szeOZpgb1+Um9rKebZjHL/DynWtQdI5naKUODWVQwNR4tkz770DTisnp10kvvncSV67oYYtjZdvROBcOGuATwjhEkJ4pm5jzLq7pJSVUsomKWUT0AdsPt14TZY+qVyRVVUeYuk8ZS4rNywrZ1mFix3NQQ70R2ftdacMui7gpD+SNsJvEykCDg2rqsxwPoVTOXZ2hVEUwYnhBMeG4hwZjHNsKE6+KKnx2TkyYNwuSonPbsWiCvb3RXmuK0yZe3Y/5tMZjBoZZlPx+oIu+Y0vPc9Dh66Mn+pCIvQh4CkhxD5gJ/BzKeWDF3ZYJpcL+aJOwGVlZ/cE48kcT54YI+S1c3jA2GtmC0XW1XpZVuHm5rZyNFWwqtrDsaEYfRNpOseSbG7wc2AghhCCHc1G/6R1tV58DgsCKHfb5izsGIhmaCw39qzhZI5fHBqieyzJ6moPTqtKagH9l6q8dm5dWcHG+kDJiNP5Ir/zzRf48lNdi/Y9XSqWdCqlycWjYzTOAweGiKRyPHBgiMFohoBToy3koajr7D4ZwaLAhjo/qqrQG06hCGgqd/F8Z3hGyGiKGp+dKp+dA31R1tf7SGaLqIrAabWgCMgVdXRdoiqC/X0RdGkYeziZY1O9n40Nfk6Op3js+CgrQx421Pvx2i14HRbSeZ1IKo/fYSHgsiGlxKapfOqBIzMy5t53XRN/9ZrVKPM0irscuGLrgU0uDkPRDP/+8Al+un8QiyK4rrUMv0NDVZWSw2pDnY/9fVE6RhMsD3nwOTQKus7RwThbGgOk80W6x5LEMgUEhmFXuK3EMwXyumQomiGVK87pud7RHKSgG3nPbSE3kXSehjInUhqz6bbGAEUp2ds7wYH+GC6rynWt5QzFMrSPJEjni1gUwV++ZjUP/N6NfPQ7L7K/L0JTmYvv7e6lbyLNPe/egnoZG/F8mDOwyRkZiKT5nW++QP9EuuTJXVl1ykCzeR27piKB9pEEQhhZVdF0HkUIxpM5/A6N5SE3dk0lmS1wfDgxq4TRqgq2NAYo6EYGVSyTp30kQWuFmzKXlbFElp6JNLmCXtrTOqwqdk1lPJFjWaWbk2NJ8vOEpXY0BynqEqtF8NHblrOxPoDDqpItFNnZFUZVBNe1Xh4VYqdjzsAmL5mPfudFVEWQzRdLHSKODsWxqgLvpKFuqg+UGqzV+OzU+Ox0j6e4tqWM0USGrtFkyfNcH3BQ7bPTPZ4kX5RUeGzU+OzYNZVnTxPOK3dZUYTg6WmyRFZVoYhEl4agwpTY3m0rK7j3iVMFEQLY0hggX9SxWRR2doexqgqNZU7e8cXnuXF5OX98xwo21vu5pqUMZYnG700DNjkjn3jlSt731V04bRY2NPhJZotYVUEqV8Rps7C/L8LO7jAVHhv5oo7LZuGZzjDlbitj8QyKgK1NQUbiWfJFnaZyF8+0j6GpCluaA2QKRWwWYyasDzgYimawKLCqxkdRl0bF2LTqzimR+yqvnY5J4f26gIPusRRra7zYNbWkJqIIwdGheKlmPFvQiU2mbOYKOr/3nRdpLHOyrNLNm7fUsaZm6RU7mEvoJUwmX2QskSXgtF5Q8fNvP3eSH+zpYzSeJZrJ01LuZm9vhLaQm3K3DSRkizouq5ET3RNO0Vrh4tmOMOtqvTisFnJFHbtFmTHLuqwqa2t9pbxoTRWsrvbSH0ljVRUGohlURbCx3s+h/mipOqkx6CRTKOK0WvDaLUgk+/uMmvBavwOnVSWSyjOenCnUJwRsaQhgUeDgQHzGMl5TBT/96A2srPJesO/xpXLFFvRfrfxs/wCfe+g4J8OpUiriYgug7+wK87/PneRQfxSXTcXnsPJ81zibGwKMJ3NUuG3YNYWDAzEmkjkKumRLgx8d6BxNsqrKgy4NATq7pmDXVNoq3TN6Sk3RFnLTPpKYYWy1fgfpXLGkiKIqghtayzjQH5ulkvLKtVX84uDZ47pCgMOikJqn19O6Wh/3ffi6Up775YK5B75C6BlPkS0UGY1n6ZtIz8gj/tn+AQ4Pxnjw928s/QCllKXc7KNDMRKZApsbAjNCJvmiztPtY+zunuCFkxP4HBrvu76JrY0BfnFgsKSBffvKSh7941t45PAw6XyBwUiGgUialVVuUtkitsmla66gIzAyqSwWhR3NQdpHEownc7PyxIWAaq8dj11jdbWX8WSOgFPDZdNQhSEkcGIkyUg8Q0PQyYu9kTm1v59qH2NZhYv2M4jTr6/zEUvn503dBDjQH+WeJzr53VuXnf2PcZlgGvASQgj44Z5+Ktw2w0szDV0aXuA3f+FZcgWd37imkXdsr2c4luHXR0d4pmOcEyMJRmIZ2kIesoUi2YJOfcDJI0eGKeqS7c1BHjw0xIOHhlhd7eXd1zZwy/Iynu+e4Injo9z8mV/PW4drtShsqvcTTuboGE0wMpmXPL2tSuI049tUb5QFTnW7mNLEWlfrw6EZiRpem4X+iCTgspa6R9gtSmk5DYYKiILhyS7ocs7abVURZzTeKf7jkRO8fHVoyWhpmUvoJYiu63z24eP8z687gJld/xqDTu5YG2I0luXQQNTIouqawKoqeOyWUihoivqAA59D4+BADK/DQlulZ4awXNCpsbrGy4nhxIwmZHNhUQQtFS4iKUPg7nT8Dgt1ASe6LnE7NPb1TpAtnPr9BZ1WckWd1go3qmJI0PZHUoSTeW5dUUFBl3SPJZlI5VkzqY8VTeVLSSJtITcnRhKEvDYCTmupgEJiXO8WKrK3oc7Hjz58/WUTFzaX0FcYiqLwp69YSX3AyZMnxtjVHea91zXRF07xnV29fPGJUymCWxqtk+9hTvnY3ok01ZMtUGPpAj3hFKurvcQyefom0oRTeZ5qH2dHc5B8USc8T4lghdtGc4WLjpEEE9P2qGUuK41lTlRFcGwozsGBGBZFsKnBj01TyRZOjSmcys1Q9tjaGGB5pQebZuRQr6zy0DvZv+p0Y9xQ5ys5yAq65NDATOmjDfU+VAELkbPe1xct9Ue+3DENeAnz9u0N3L4qxInhOHv7InxnV++M56fLtM7VpHtFyIMu5YzC+pF4lpF4FpdVJeS1UdQlrRVuouk8y0Ke0mv9Do1lITfRVJ7BSJrRRHaGlK6mCNbU+tAUwa7TpGILumRX98ScRf0CWFvjxW0zVgs2i8LOrhg3Li8jfYZG49M10CJzXGT29UbZ1hQoxaPPRI3PviSMF0wDXvJUeGz84uAg//LgsVnP7T45QVvIPW+Tsc5RI3NqdbV3ljRvMldkc0OApzvGGEtM9YKyEPLaGI5lqfbbSWQKM5qgT7G9OcjxoThFXXKwPzKvXlbmtGKEHc1Bnpv2uhqfnWxBwefUODgQx2VVS9pmcNrWocxFW8hw3ulS0hNOzZLF3d8X5cZlZTw3rT55c4OfgUia+oCT8VSO91zTyNu3L526ddOAlzDHh+N89elufjypEDkXqiJYVeUhr0vaRxIoAtomc5UPD8TIFIoETlPmnCKVL85wCMWzBWoDDqp8duLpAo1lTk6OJ2fMjCtDbqN5twQ1kqKow56eCQJObUae89bGAKl8EU0RpfTH6V71creVkM9OoSgJODVsFoXeiTSRVJZrmoMUJfSEkzSVudClnDWz1szRJ0kCJ8Np6oNOyl02FEXQNWrs7TfWB/jaB7afMZ4upSSSzuO1WVBVpeTlLxT1SxZ6Mg14iaHrkkeODPO1Z7rn7MFrtSjUBxzUBZz0TaRKTcJsFsHmBj9dY8mSNK7HplIoSgqT2U0OTWFDnRHL1aUkX9Apc1lZXulmKJYh4LKiqQq7usOsrPLw62OjaKpgTY0Xl83CybEkHodW8lRbVIUN9X5yhSIOq4XDA4Yet9tmYTRhZGZNGe+2piBWVbC10U88U8Rtt/DCyQm2NwWJZ/N0jSVZXe2l1u/ghWmqlhVuG7o+c2ltsyjUBhwUpZwxC6+r8bKnN4KURqw64NRYWeXlPdc18eFbWueUQ4qm8rzmv59kMJJBUxUjlqypxDJ57BZBmdtOOJnFZlFZVulhechoXLe80k1byDNLtnixMb3QS4j2kTgf/MYLM/SR79pYww3LymkIGjpWVV47P9k3wJee7JrVeG0+drQEyBckXrtm7GXjWRqCTuyaws6u8KzeUKuqPPSGUyROWwJvrPOVZqLu8eQsiZvllS48dit7eiZwWFXW1/rQVEF4Uva1oEu2NgY4PhxnbZ2PfEEnlsoT8tlJ5Yqz2q5sawqw5+QEDk3FYbPM6krZVO4k5LEjpZGCmcgUyBV1Ixd7NMHaOj/hVI6iLvmNaxp58+a6WWWFv/2/u/nloeEZj1V6bEjJWVvWAtyyooKPvWIlq2teeoaXmYl1BXCwP8rd39hNmdvGYDTNWCLH797ayp++YuWs1/762Ajv/+qusx5zVbUHi6LMa+ghj+2soaPpNAadNJU7efz4/PLgdQEHfROn1CKFMGR4q7x2Dg1ESeaKCGF4r6u8dg4OxEotWu2aQkPQSa6gE/LaZ+yrV1Z5aJ/W9+pM5/Y5NBTBrNLFT7xyJTevqCDotDKayPKfj56YZbwCQ/Tv5BwxZb9Do8Jjw223EE7k6JlIIaXxGe/aUMMf37HiJYnqmWGkJU4qV6DMbeXJj9+GqgiiqTxPtY/x6vVzNwrP5oszWoBOd/ZMJ5zMnbH/UdBtpT7opHM0OW+Tt5VVHhLZAmUuKxOpPE+3j83y9q4IudFUI+2y2mc3iu7tGgPRND3hNBPJHF67hUy+yPJKNz6HxomRBEeH4mxrChipmdUeOkYSHB9OsKneP8spdnQozrUtQdpHkqVZdToum0pqUj9rvo6U+/si/PLQEMeG4rx1W/0s4wW4flkZ+aKkymvnQF8Ep82C36HhsKkc6IvN0LdeEfIwFMsQTef58d4BHjgwxDu21/PR25cbOeSLgDkDLwF0XZ6zYsTP9w9yzxNGs/OhWGZWBpXXYUFKztqMDYyl6omRxKzwjBDG3nV39ynx+6BToy7oxGuzkNMlE8lcyVO9rcnPru4IK0Iejg0b+/D6gIPGMidjiRyqwoz4bbnbiCF3jCSIpE+N8/Tw047mIMOxDH0TKQq60e6ltcJFoQj7+iKAkUp5fCg+I4MLKNUqu20qL/ZESwY4l+fcYzeMtXfaCmIKiwIbGwKzNLMbgk4Kus5YPEdbyI3faSVf1PnYnSsXLKxnzsBLnJci9zIYTbN/mkayRYFNDQHGEzlG41li6bMbLsCaGi+7uydmNTcHQ51yZ1eYhqCTgFNjX1+U5SEPz3eFWVXlIVfUS+mPioCibtToTt+rDsWMiqNqnwOEYZzxTJ7BaIaxRI6xRI71dT7WOSxMpAo4NMP7u70pQCJbxKEpqAoz0iQjqTyJTIFUrki1z47HbkEIwTUtZTM6SFgUQSSdZ1f3BNc0B9GnXeWm97pyWFXq/HY6RpPzXvDqgy7cVpUdzUEUIUhk80Z++KRBex0WGsqcPNcZJpzM8Y8/P8yPPnz9gv4GZ8I04CuQ7rEkn3+sY8ZjmxoWlsQwnQqPDYsi5jTe6fSEU/SEYUWVp7SfnnJmLa80ls8Bl8ZwNEP7aJL6gIPWCudkna7AogpG4oZ65On5yq0VLo4MxtjcECCezlPucvHY8dEZr6kPzEy60FSBy2bh2GT8e3COLb4i4KblFTQEnWTzOs91hbmmJchznWGEmBnSago6OXJaU7vTCTi1kj52udtKQ9CJzaKU+kF1jyV54MAQfqfGa9ZXUeV10DOemlNA/1wwDfgKpKncxXuubeL4cAyHZiGSzp2xS8J8tJS7KOiSa05LsJiPY9N+5AORNLFMnpDXEK6zqgoBl5VrPDYE8FxnGIdVpdbvmJEMsrnBj82icKAvgsVidDxsrXBzZDBGW8g95wwoMfabLpuKLiVHh+Ls6YlQ7rbSFvKQyRdLcrZrajzoOnSMJWddCIpFaaR51vtnZI8VpWR5pYsTI3NXO21u8M/ohTy1cpiOx27hprZyAg6Nvb0RBtwZPvaKFWf9Ts/G5VX4aLJobGsy0ih/sKePR46MsOfkBGtqvLSF3As+xpHBGALmLOE7G62VbvJFSbXPzpMnxohljKVqQZc82xlGYiRG1E2bPQVGuOfZzjDr6gO0VrgIOq0cG44TzxbYfTLCCz2nDMs92d2hqcxFIlvgQF8Um0UtpY2OTW4XesIp1tX6qPDYODGSxOvQ2FTvp9w9M0abLRS5rrUMRRF47afmtuPDCXyO+eO5e3oidI+nWF83v6JHPFPgieNj3L9vkLqAk6FohoePjJzjtzobcwa+QvnloaEZHmaJ0Y5keeXZDdhts9AWcrO/LzIr9roQNtYZ8d11tT5AGrOqpvLyVZU8etT40QpgeaWHvb1GqiVAoagbM35LkGODMcKpPA6ryppqLwcHDLXJtiqjzO/YUJzGMid9kQxPtY9R47OzptZLOj/7YuN3ajNCZdOVNMcShjjBmlovRwdj7O83UkrLXFauay2jUNQRQiAEbGrwkyvoHBqIzTpHyGvj8ByPOzSFVdVeusaSNJa5yE/6BcaTWRrPc/kMpgFfsZzernOKcDLH8kr3nDnMUySyBbrHU2xuCLCnZ6LUf8hpVVlZ5aGoSw4NxOaMuZa5NApS8nT7qSyxHc1BnjoxxqoqD1ubgsjJBmlTvX87R5PUBx1oFoV8roiUlKqe0rkix4bj3LKigmc6xnlxWmeH9tEEbSEvB/qjDEQzDEQzbGsyPLtTUjwvnOECtK8vym0rKsgUizx5fHRGpZJNUzg+HJ+xFLYocMOyCrY1BSjqRmO1VK5ALJ2nZXKZPxrPUpTQVOYk5LVjUQXPdoyjS5hInRr7lsbAojQcNw34CuSfHzzK/XsH5nzOoooFpfeFkzl2JnNU++wEXVbK3Vae7Rgv7SV3NAfJFQ2ROEUYxxyJZ6n02Ejmiqyr9ZW6Ndo1Y6d2ZCiOpgpq/Q4mUjmWV7jIFHV6w+kZWU0hrw2ndkr6prXCTedoclYf5mxBcqA/yk3Lyzk2HGc4li3pHJS7rWjqmb33bSE3I/EsFlWwtSmIrkuOjcSJpQsMRDJsavCXDLjMZaXaZ5+1b24LuemdSHMybISWLIog6NToHk/RPZ4qrS5ayl1UeG04NJVEpsDf3bX2rH+DhWAa8BVEvqjzmV8e494nOud9jUNT2dOz8GXxVG+hiWRuRkrl1DJ0Xa2PrrEkdk2l0mOb0bh7VbWnlIu9ucGP22ahYzRBldeQnY2mCzQEnWxtCrC/L1oy0JDXjkNTqfDYKOiSbF6fUS4IRpZYXdCBRVHI5Iul7YLNorKmxkuuUOS5zjDLK930TqRK+2JNETSUOQk4rbzQMzEjPi4EBBxG0+8dzYaS5oY6P/v6ItQFHKUVAxhxZSmNlcj0hUhBlyWHVlvITSyTN5q1jSVL8kRffM/W80qtnI5pwFcA/ZE0337+JP+3q4+xs+Tn+p3WBUnLTLGjOcih/uisvOcpOkcTZPJFDvRHqfXb2dIYKBl3ldeOy2Zhd/cEe3oiBF1WBEbsdwojBJWiscxBlddBtqCzrzeCxIjteuwWrNMqfYJODZ/TSoXbxp6eMMsqPWiqYEOdn4FIiifbZ6ZxnhhJsLbGi82isLcvSlO5k2i6UIpPA2xvCiCEYF9fhNZKN7kBQz+swmMDKWkLuWe8HgxZn7kE+qZ/b7tPTszZ//j0HtPnw4IMWAjRDcSBIlCQUm4VQnwGeC2QAzqA90spI4s2MpMzIqXkiRNj/O+zJ/n1sZEzNsqeTuG0yp0zEXBqtI8k5jVeMOqGp+iPZKgLOPE7NCLpPOPJHJ5p3txwMkeNz47Nos5K7zw5nubk+MwMp6nQUFEa/ZGay130R9J0jyfpGkvitRsz+upqH0JAKpufs7Y5kS2gKBpbGvyTFUXZGZI/ugRVGLrRu7qN4gi71UgYUYQo1VNX++zUThb6WzWV7c1BEpnCrPOFPDa6x5Mz/iaqImgLeaj02HAvogTwuRzpVinl9Mvbw8AnpJQFIcQ/A58APr5oIzMpIaXkh3v6KXdZ+cxDx0hmC8QyhRmxx+k0lTkpc9mwTAqwD8cyFHTJ8ko3FkVgUcQZk/6nqPLZS0vghfJ8V7iUhri/L0rQacVts5TkfAaiGZhsjLax3s/x4Tj904TvwKgi8tg0To4nqfLZGYpmiGUKs3S2YpkCG+t87J1MlwQYT2ZZXeXh8LSYdLnHNiPFseK0POQpT7vbppLMFknni4zEsvRH0jMuMlVe+yyvvKYKllW66Rw1ZHFr/Q6EgL6JNNqkyJ6URomjIuC3bmw+p+/zbLzkS4GU8qFpd58D3nz+wzGZTqGok8oXufsbu3muM4xFEWyo92NVBXZNNWphpU5RN+p3FUXgc2i8OBmXPJ2ppe2qag894RTJszTHTmaLC9aRmg+bpuB3aiUDDroMg67wWPn1sVEq3DZWhDwMRtPEMgWayp10j50ae41flOLQIa+NhqCTiVQOVShE0jmc1lM/Yb9Tw6IoHBmOs7nBz/HhBEGXdYbx7mgOMprIEknnaCpzYVEFAoHbpnJsOIFmkTQGnXjtFmoDDiKpPKmc0Zx8X19kVh52vmgIJdgtCk1lLnwOC6m8Tt9EmpvbKnm6fQxNFfze7cu5Y3UV9cHFlepZqAFL4CEhhATukVLee9rzHwD+b1FHdpXyX4+e4IWeCdbX+vjGcyep8Tlonwz5FHRZCovMJ1OzEFQhSC+gt25POFVKL1wIWxsDWBRBIltgZZWHvok0tQEHsXSeDfU+BiMZRuJZwskctX47Tk1hNJElnMyyLOShEE7RXObCYVGJZwv0T6Rx2yysqfHisVnIFnQUIWifzIjSFIFVTbGqyqj6yRWKDKTzSEnJW57IFqj2GfW7jWUunu8K47Eb7UudVhWAcCLL4cFTq4CpENuO5mBJ/GCKnV1htjQGSn8HRRjGkSnoHB2O49AUXrm2mnW1Pr69s4e1tV7+6Q3rWXeGJI/zYaEGfL2UckAIUQk8LIQ4KqV8AkAI8UmgAHxrrjcKIe4G7gZoaFg6WkMXk4FImidPjPLEiTGeOD5KPFPgsWNGuOL0CiCPTWVZyLPgPe/paIqYtTQ8E1KeulhUuG3kikWicxRCWFWjE+GU2EClx0oqV2B39wSqIqjzO0pL4B3NwRktVooSKtxWBiNpfn3sVJhmfZ2PwwPRGZ0UpsIyAHldzqgM2lhvtDid+mhT++ahaIZEtsBQLIvXrrKs0s1oPFvyKjeVLzyhQlMFqmKsJNK5AooQbGzwk8kXsSgKmqoQSeXZ0eLhL169itdvql200sG5WJABSykHJv8fEULcB2wHnhBCvBd4DXC7nKcucXK2vheMcsJFGfUVQiJb4O9/epj/29171teqAtqqPBwfTsxIZjhXHDZ13nrYueibSDESy1LltTMczxB0WllX6+VA/0zHTbnbRqXHVjLglnI3IW+R/kiaxqCT4ViG61qDpHM6RSkJOrVSsoYiYDiWnZWyub8vOisPO1uY3wknpaDKa2cgmsGuKayq8nKwP0JLhVG+6LVbWFHl4WB/dIaOV284zY7mIHtOTsxoT3p0KD7rs7ZVejjYHyM1uYJpq3TzTMf4jHDUGzdWc21rkPV1CysXPB/OasBCCBegSCnjk7fvAP5OCHEnhtPqZinlwuMSJiUsiuDvX7+Wv37dagq65Jn2MX7nm3vmfG1RgteuURdwcHI8hVVV0KUsOaMay5wMRNKlIv7TUYQR+qjy2jmenpmFtbraQyxTKCllhLw2rKpCjd/BkcGY0YB7MvQznswxnsyxpTGAKoxSx7F4lvbRJAPRDNqk5rPdqjKWzLJqMnMrW9B5pmPupfi2prm3A/UBR8l4FQHrav3Ml5sxJQHUWOZkechI/JDAjW2VjCWybG8KoKqGGMLp8rTFyc/ntltKRR92TWF5pZtkrkCd307fpKNNUcWMMQxGM7z7mkYePjzM4OSFw2JRUcTFKTNYyAwcAu6bFPyyAN+WUj4ohGgHbBhLaoDnpJS/c8FGegVi14w9mBWFeCZPcdrvqj7owGW14LJZ2NcboaBLnu8Ks6HOx0kM59BbttSzvy+Cw6rydPvYrGWxVRWsnWxT8mJvhK2NAQajmVKz63i2ABI6RpNUeGw0BZ1U+QwZm+FsluYKF6uqveyaVrA/xQsnDV2r7GnKldcvL2dfb6RkCO0jSawWheYy16yuEGAY5nxpnUIIbl1RQbZQZCKVJ5WbvXRvC7lxWFX29UZpKnNycCBW8s4PRTOoqiCT13HbLGTyBQq64YVO5gqsqfGWeh5rqkLPeIoVIQ9+p0ZymgbX6movK20aNk1lX18Er93ClsYA17eWURdw0lzh4q6NNezunuBt2+qNlqgXibMasJSyE9gwx+NLpwPUZUwym6djJMmf/egA772ukY/fuZKJVI5wMkeZ20oiUyhlTl2/rIz3XNuEy2rBalGwWhQ+/soVfPnJTvom0jPE7jRVsDzkYU9PhDKXUZ/6xImZSQ7Thd9r/Q52doUJ+eyUu600lzspFo2LRlvITX8kPctr7bVbGM4VURXBa9ZVsa7Wz7d29swoXfTYLKyo9tA1T+MxfTJveMrohACfXSPospLJF/n1sVFuaavAqhY5MhhFEYa2lV0zShEfPz6KEHDbigqe7hgjWzASL3rDaaQ0RPLaR5I0ljlRhODwYAy/Q6ModUZiWU6GUzN0uo6PxNlU7+fItNju4cEYmxv8vNgbobncyd03tbC1MchP9w8SSefY1lQPwJbG4OwPeIExM7EuMS6bxvp6P//329fg0FTu3zvAfS/2kcoVjSbaVpWWchdWi8pfvWYNK6pmN926+6ZWblkRYn9fhIMDUQYjGUJeG9/eaeytx5M5VlZ5sFmU0h7SZVWxTFsLGrO7HzhVVG9VBdsaAxwejE0W6BfxOSwlJ5bbZqGs2sY/vH4NX3m6mwcODhLPGMa7vSkACMaSWXZ3T7C+zkdjmZNwKlfSZHZZ1VLmWFOZk1i6QDSdI5LOl4ox7JpCKl8sOZx0acRYLYqgbDKnW0ojJjzVZ+n4cIIbl5XzfHeYbEESTefY05NlRZWHLQ1+9vVFcWgWmsoNkfgp493RHKQnnGRPT4RN9X56wqnSquH4UByrqtA1luKexzv4i4lDk2qWl9YxaxrwZYLHbuTgvnxNiIFImn19UW5dWcHbtzWUmmzNtYQEQ/1iV3eYn+wbYHW1l9FEtlS2t7LKwy0rKmkIOtjZPZW/7OXoUJznOsOsqfEQzxhNuff1RVhf6y1lSeWK0qhAqvbi0BQ6R5PEswXqAg6S2QJ3rK7ihuVl/OH39pVUGkNeG8sq3dg0lSdPjOHQDJmZnV3hkne4eyxFyGsjnSsSyxRYW+ubN71Tl8Ye9fR4dJnbOiOne/fJCW5uqyCazpPMFrBpgraQG5dVpbnc0Lw6OBDj8ECM61vLyRWNWO31rWUlUb7BaIYav5PBaJaBaJrxZI5qn81QzfTZKeowGs+Szhf429e2kM4XCbounId5IZiidlcIjx8b4emOcbrGkqhCcPOKCm5ZUWFoTU3y2LER/vLHByn32Ng7KXAOp8JENotCS7mLqkkdZosieLZznIDTyg3Lyo1wUMBBY5mLV6+vZiia4d1feZ7e8GyRN1VAQ9BF1/jcS+ftTYFSLnGN345Ts9A+msBrt1DmtuG1W7BpRsplMmvoWyWzxRl51Ne1lpHN66iqsZMt6pKusRThZI5Kj42ReJYVIQ9Bl4bECMlV+ewIJOFknoBTo6BL7JrK811hpJQkskVuWVFBPJ1HUxWimTwVbhuFoo7facVmUVBVhdXVHuoCTu5YU7Vof8P5MHWhTUrc92Ifn3rgKPlCkbqgE6fVUopnTnUeODpoKGDcsqKCuzbUcOfaahyTSQ8AE8kcP90/wGNHR9jZNc6ykIdMXp+R9LC9KcC+vui8YZ/GSRH6oi7JF3Uy+SI+h3bGAgFNFWxuCNA1lsRpVVEVQTiZI50rkinobKzz4bRZyOSKOG0qT7WPc8PycvIFI3RVF7CTzuk4Jg2zodyJIhSyhWKpTamuS359bISCDm/YVMsbNtWSLxTRLCo3Li+fs3vDhcY0YJNZPNMxxoe/tWfOTn5WVXDD8nI+8+YNlJ2WhHDP4x187uHjJcO8rrWMo0NxVld7CCdzaBYVTRE4bSpPnEHgfS4cViM99PTSwdPZ0RykoEv2nDTUMm9bUcl4MotNU4mkcug61AXsRNN5sgWJ36lR1CVuu9G9wW5RcDkspLOGeICqCPb0TLC21seJ4QTXtZYxFM1gt6q4bRYePjxMtc/Oq9eGeOPmep7pHGd5yMO6Wt8Fb50CpqysyRxc11rOg79/E0+eGKU3nGIgmsFlVSlz2yhzW9lQ55/x4/ze7l5GYhn+49ETpViz02r0CKrx23mqfWafpu3N5+6RtaoKlV7DgENeo2rH2P4bs54QRs5370SK8Wllk8lcgX19UW5aVk6Fx0b/RJqn2se4blk5ZAqMTbaL2dIYwKGp2CwK3aMpRmOZUqXVDcvKGIhkaKlwkcoWeLE3QmuFi4Sm8uYtdaysctNa4eZfHjxKNFMglulBUxXuvqmFN26uAwwp2oFIelISN0ssnSeazrO1Kcg1LYtXQjgd04CvYqp8dt6ytX7e56dWZ8lsgX/+xVFSOaNG1mPT0FRB0G3lcH+UseTMWTzo1JilJL8Aouk8xaKhgjmayM6qwZ2O36mhywKVHhsgWVfr44n2MZZXuqny2fHaNYSEwwNRlld42FDnJ5bOIybzwDP5IluaAjzfGaa10s2JYaNLYZnLSluVh62NAbIFnWg6x2AkzUgsw/97rJNwMsfGej8WRWFllYcnjo/yVPsYT50YZTSRQ1MUmsqdBJ1WesIpNIvCyqrFKd6fC9OATeZlar/XMZognS+SzuukIxkgw4Y6H08cH8OuKaXkB0PBMs+JkQQv9ERoLnfNiE0vBAkcH0ngmrbnnu+F+aKkP5Kh0mvH59BoCDoZTWQJuqzoUvL4iTE2N/jZ0xPBOa6wvNJjzNTLyxmJZRiMGumd0zPEmspdPDvZ9bGp3InPbuHpyft2i8KO5iAnRhKEkzmjr3DQSZnbykjcCDcVdJ2A08pzXWFq/Hbetb2Bl60OndN3cC6YBryEmOpHezZGYpmS0kVvOD35f4r+SJpllS48dsMrizQadb1+Y+2s+HK2UCSazhNLF0phrCk8dkspfpvJ66U2KdOp8NiILTDnui7goM7vYCKd50B/FCkltX4fQbeVA33ROQsvVMVItazy2YhniqXlfCSVp30kUVr+T5UxpnI6tsnMt2xBL4WkcoVpAu5lzlLlV33AQX84PSM3OlPQ2dUdZntzcPL7zHB4MFYS0oPJ8Nukf2AgkuGRIyPcfXPrgr6Hl4JpwEuInnCKqklFi7k4NBDlzZ9/dkb7kClCHhtVPjsj8eys7oGff6yDlnKX0Thtct9W63dQ5rZyfCjO5sYA2cljWhQjvvrCycicY3BbVfK6xOewlMr+zoZFEezrj84ocTwyGGVLYxCXzUJLuRtNFQzFMlT77KRzRY6PJBhP5nHZLISTOY4OxY2sr5CHEyNxxpM5qrz2GXI84/EMNy4vJ5Mr4NQU8kWJLiXlbiu5gk6Vz06uYHjEB6KZUqvTXEGnqBsa15pFmVVeeWwoxupqL7FMnmqfneGocXFThLHUv5CYBryEaCxznfH5F05O4LFbWFXtYXhSUWKK4XgWv9OKz2mkKZ6u5jEluGbXFFZXe3HbLRwbihPNFHj8+Cgb6oxUwjK3FafVgs0iUIUwpGALOqtqfKRzRY4Mxbi2pWzO5uPz0T2emlV11BbylgocphqUeWwWLIqYkfRxdChOhdtGhdtKQZdMpHLU+B1YFGFkirltbGkMMBLLoCPoD6fQBaTyOssq3OQKOmOJHNpkSKo+6KTSa6cnnJpRYNFS7mIkniExhwiCQ7MwFMsQTuYYjGZ42apK3r6jnjduqqPKZ1/w9/BSMA34CmJllZfPv2szf3H/Ia5rLeOHe/pmLD+nlrqKgLU1Xvom0jNSFpvLXXjtWumHqyqCDXU+rBYFl83CH758GZ9/rLOUYWW1KPicGiPx7Az95XSuSKXHRq3fwfHh+AzdrPmYLvETcGrEMrOX3/FsgdU13hkGvLrai0UVk0vvma+f6i88xaYG/4xSzKKUpefzuixpX4GRJ95a4eLQgFE6qKnKjK1EQ9BJ9aRx9kfSlDs0PnB9E69cV01rxcK7X5wvpgEvUf7hZ4fpGE2wosrLsko3Dx4c4pEjw7xqXRU2i0JPOMXGej+D0TSD0ZlaUro0jHl1tRd3MsdQNE1TmeFwykwrtSvqspSD3BZyc6g/VnrecGoViaTztIXcuG0W2kcSxDIF8kWdkbghHNda4SKSMgTuFGEI1dksChKjjcpILMt4MjejM+BEKs+yCjc1fkPNY3rj7hd7I6UKounZXKezqtozK0NMmfQfaKohMFfhthFJze4b5bFb6BxNEE7lERjbD02FtpAHRQiKupwxQ79qXRX/711bzunvt1iYBrxEOTGSoGssOUPBAuCBA0Mz7q+p9s4yYDCW49N1jqeyqCo9Rhx4aiZeVummzGWdVew+hdduKc1cDs1IMZw+63WMJqn22bmlrZz9/TGODMVLUjXpfJFcQcdjt5DMFkrOoEJRoqqC0XgOwczZOVfQCXlsDMezs3r9TlHhttEW8swQ5At5baUZdGO9v9SpsanMSb4oS84ugKDTSt+EMctLjO3H8KSaiBBwx+pQSeBuNJFlRejChYnOhmnASxBdlzzdPoZFEYS8thk9kKYIOq0ksnmODMW4YVlZSeDNbbdQKEqOzNVzExhPZGkud6EI2N4UJJYxPMOnG68QUO21Uxd0lkTegi4rueJsoxqMZoik8hQnJW2f7wrTUuEqVQHFMwWOZU4tX8tcVu5cW8XqGoVD/TEcVi+aqjAYNTS16gJOGsqcHDxNFaTMZSXo0nBaLYzEMqXuEYlMgUg6XwoDOaxGKaaUkpFYBo9Dm2HAJ8MptjcHiaRyHB9OsLUxgJiUnX3T5jrevr2+5EjUdUnxImYzno5pwEuQ0USWgi5ZWe3l+GkhnKBToz7o5EB/FFURtFa4kRgz68F5mnJNXQCWV7opTgoHTKfKa2Ntja+UCRVN5+mbzN4aiJ4qLuiPZPDYLbOqj1aEPKiKmKGf3DmanLUnnSKazvO93b3c1FbB27fX8+3njaynD93Uyq7uca5tLSedL9IQdNI7kWZfb4RsQWdZpbvUXHyqo4TPobG62ktLhVGk8aujI7SPJPDZtZLhVXkN77yURucKXUryBZ2Xrw7xmze0sCLkZnnIg2sOPWdFEShc/PzoKcxc6CVGUZfc80QH//LgMdbVenFaLeyfFoKZzyg8dsuM3rpOTSHotqEIUdKkmiv8NEWF20aN3z5j2X0mNjUYmU8eu4bLplIoSjL5Ipm8TpXPRqEoiaTzFHU5owhiR3OQKp/R4eHWFZVE03nqg040VTAczdIxmqC10s3e3glqfA4++eOD9E2keNW6amMJX+NjbY2XDfUBdF0iBHPGzrOFIseG4rSPJPj1sVGCTo1nOsbxOSz8zs3LLmjyxbliFjNcITx6ZJhP/+LoLAmaKWWNhqDzjIUA25sCxNIFFEVwZCjGhjo/HSNx1tf5S9lGp1PuttFc7mRtrY9MXic76bw6NhQvhZ7m4vSwUEu5i7dtq+c/Hz3BzSsqEEKwptrL9uYguYLOi70RrKrCB29q4cWeCe7fO8DP9g9iUQSRdI5Xra3m7dsbZuVYxzN5fnV0hNesr5mVcHKujMWzlHsubX3vXJjFDEscXZe872u7eOK0znhTpLIF1tf5sGsqfROpeSVjc0VDu3iKvb0RNFVQ6T31o20LufnwLctorXAbHRLscyciSGmU3X35qS52dU2QK+psawqwrzdCrijJTtsLC2GogkwJo//Lm9bjPu241y0rJ5zM8f3dvTx8eJiHjwzPCAv9ZN8ADqtKtc9OffCUDKzHrnHXxtp5v7tz4XI03rNhzsCXMaPxLGOJLM92jPN3Pzt8xtdO7wQ4F9uaAiXP6+lYLQqfeOVKdAnvvqYRq+XcFBX7JlL8bN8AA9E077u2mVTeUPg4MhjjyGCco0OxkgzO2lof3737mpKg33wcG4rzP79u5+cHBnFoKu/a0cAnXrXqnMZ1pWAuoZcY9+/t51dHR9AlPN85Pqsn0Fw0T2YKzdcuZU2Nh8FodkYGVpnLysZ6P3ff1MKORSh3O9gfpdbvIDBHjWw8k8dmUc/54hBN5fE6LJekkP5y4UwGfHHEa03Oia8+3c39ewfoGEmQnpS2ORtdY0kEsKXhVGK9w6qytSnA6movhwbi1Pod7GgO4p/sgRt0WckWdDY3Lo4A+bJK95zGC8ZS91yNF8Dn1K5q4z0b5h74MuO5zvFSDvPpbSvPhNOq0lbpYXnIjdtuIV8scrA/NqOx18F+o+1IwKmxqcHPieE4d6wJkS/qaOr5X8vPtiw2WXzMGfgy4o+/t4/3fmUno/MsmYWAl606Fd7wTuu9+9r1NaQLOnaLwjMdY/RNpPE6ZjqKpjZLkXQen0PjjjVVDMeyWBTzZ7BUMWfgy4T79/bzy0ND84rAvXlLHXUBB01lrpLczM0rKmgsc+F3aBSKkp8fGCSRzfPEn97KZx46RiSV58blFezsGqc3nC5lSUkJT54Y46O3LeP91ze/pKWtyeWB6cS6xBwdivHX9x+akf1U4bHNmIXL3Vae+8Ttk+LqCyOTLxJL56nw2HjixBj/8chxXuyN8Mq1VTxwYIh/eP1afuOaxkX9LCYXBjMOfJmQzBZ44vgoG+r97O2N8M8PHqV/Il1K1leE4Qi6flk5jx4ZodbvoNJr4+a2inMyXjD2o1N70pvbKmitcPHokRGODMa4ZUUFG+v9i/3xTC4BCzJgIUQ3EAeKQEFKuVUIEcRo6t0EdANvlVLOL+p7FdM9luRvfnqIZzrGyRV0Kj027Jo6K2tKl0YlTlOZi0++ehWvWETR8LqAk/de17RoxzO5PDiXGfhWKeV0LZY/Ax6VUn5aCPFnk/c/vqiju0L4ytNdPHbMaML1ubduoNJj5/mucZ7rHGdvb4R8UWJRBA6ryjt3NJiGZrJgzmcJfRdwy+TtrwOPYRrwnBR0yY3Ly7ltZWVJQ/iG5eUAfPO5k/zFjw9S0CXxaT16TUwWwoKcWEKILmACIxJxj5TyXiFERErpn/aaCSnlrIwAIcTdwN0ADQ0NW06ePLlYY7+kTLXbOF+klLz93mfZ0xPho7ct53dvXXbeSfkmVxbnnUophKiRUg4IISqBh4GPAj9ZiAFP50rxQneNJUhmC6yt9S/K8aSUTKTyF6VNh8nS47xTKaWUA5P/jwD3AduBYSFE9eQJqoGRxRnu5Y/PodG2iDIqQgjTeE1eEmc1YCGESwjhmboN3AEcBH4CvHfyZe8F7r9Qg7zcCLpsZvKDyWXBQpxYIeC+yYRyC/BtKeWDQohdwPeEEL8J9ABvuXDDNDExmYuzGrCUshPYMMfj48DtF2JQJiYmC8NcB5qYLGFMAzYxWcKYBmxisoQxDdjEZAljGrCJyRLGNGATkyWMacAmJkuYi6rIIYQYBRa7mqEcGDvrqy4e5njOjDmeMzPXeBqllBVzvfiiGvCFQAixe75E70uBOZ4zY47nzJzreMwltInJEsY0YBOTJcyVYMD3XuoBnIY5njNjjufMnNN4lvwe2MTkauZKmIFNTK5aTAM2MVnCXBEGLITYKIR4TgixVwixWwix/TIY00eFEMeEEIeEEP9yqccDIIT4EyGEFEKUX+JxfEYIcVQIsV8IcZ8Qwn8JxnDn5N+nfVIW+ZIhhKgXQvxaCHFk8vfy+wt+s5Ryyf8DHgJeOXn7VcBjl3g8twKPALbJ+5WXwXdUD/wSI5Gm/BKP5Q7AMnn7n4F/vsjnV4EOoAWwAvuA1Zfw+6gGNk/e9gDHFzqeK2IGxpC7nVKZ8wEDl3AsAB8CPi2lzEJJDPBS82/AxzjVpPCSIaV8SEpZmLz7HFB3kYewHWiXUnZKKXPAdzF0zi8JUspBKeWeydtx4AhQu5D3XikG/AfAZ4QQvcC/Ap+4tMOhDbhRCPG8EOJxIcS2SzkYIcTrgH4p5b5LOY55+ADwi4t8zlqgd9r9PhZoMBcaIUQTsAl4fiGvXzLNzYQQjwBzNQv6JIY21x9KKX8ohHgr8GXgZZdwPBYgAFwDbMMQ/2uRk2ukSzCeP8dYtl40zjQeKeX9k6/5JFAAvnUxxwbMpZx/yVcmQgg38EPgD6SUC+vufin3Qou4h4hyKqYtgNglHs+DwC3T7ncAFZdoLOswNLu7J/8VMFREqy7xd/Re4FnAeQnOfS3wy2n3PwF84hJ/HxqGj+KPzuV9V8oSegC4efL2bcCJSzgWgB9PjgMhRBuGo+SSVLxIKQ9IKSullE1SyiaM5eJmKeXQpRgPGB5gjD5ar5NSps72+gvALmC5EKJZCGEF3o6hc35JEIZm85eBI1LKz53Le5fMEvosfBD4DyGEBcgw2YvpEvIV4CtCiINADnivnLzMmgDw34ANeHhSb/w5KeXvXKyTSykLQoiPYMx4KvAVKeWhi3X+ObgeeDdwQAixd/KxP5dSPnC2N5qplCYmS5grZQltYnJVYhqwickSxjRgE5MljGnAJiZLGNOATUyWMKYBm5gsYUwDNjFZwvx/DcY/vNrNJj8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "simple_uk.to_gdf().plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And we can write it to a GeoJSON file:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "simple_uk.to_gdf().to_file('./local_authority_districts.geojson', driver=\"GeoJSON\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "\n", "## Download link\n", "\n", "- {download}`[Download the simplified LAD geometries as a GeoJSON] `\n", "- {download}`[Download the EU referendum vote results as a CSV] `" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.8" } }, "nbformat": 4, "nbformat_minor": 4 }