Skip to content

Instantly share code, notes, and snippets.

@moi90
Created August 26, 2022 14:32
Show Gist options
  • Save moi90/28539fa82e55399a254d8e6938eefe56 to your computer and use it in GitHub Desktop.
Save moi90/28539fa82e55399a254d8e6938eefe56 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 50,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"\n",
"from scipy.ndimage import distance_transform_edt\n",
"from skimage.morphology import binary_dilation, disk, binary_erosion\n",
"import matplotlib.pyplot as plt\n",
"import skimage.data\n",
"from skimage.util import img_as_ubyte\n",
"from skimage.filters import threshold_otsu\n",
"import numpy as np"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7fbc8d4589d0>"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPMAAAD8CAYAAACioJLqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABgpElEQVR4nO29d5gkVbn4/zlV1XlyTrs7u7OzibALLMsCIgoqUUBMiChGVEzoNd/7+xruNVwxRy4KCCogAgoqWQVEYJe0u7A5787s5Dzd06Gqzu+Pqu7p7unJPT2B+jzPPDNT4dSp8J7wnjcIKSUODg7zH2W2K+Dg4JAdHGF2cFggOMLs4LBAcITZwWGB4Aizg8MCwRFmB4cFQs6FWQhxvhBitxBinxDiS7m+voPDQkXkcp1ZCKECe4A3Ak3Ac8C7pJQ7clYJB4cFSq575g3APinlASllFLgTuDTHdXBwWJBoOb5eLXA06f8m4LT0g4QQ1wDXAKiop/gpyE3tco0QCFUlWuJBC0tEf8jarKpIw5jlyjmkI1waGCYA0jSHt7vdSE2F0JD1v6KgF/swPODpiIKuIwGkRGgqIEBKwtVevJ06MhwBwCzyo/SGki4oEEKkXCtMkKiMiEz1y7UwZ6rEiHG+lPJG4EaAAlEiTxPnznS9ZgV55jq03Uc5dM1KjBMHqb/iZdSKcoQQ6K1txN60Hu/Bboy9B2a7qvMaJRBAKSpEP9aCtrgO6XZh7D2AWlyMubwOXtqJkhdAFBZgdnYjI5ZwKfn5GD09KOvWECvyIhVB62kepAr193Zi7NiDedZJqKEoansf+tEm1JXLkYebENJH32tWEAsolG3ugo4eRJ6f8LIyvLtbkUNDCC2AqXVjiqBV0T5GSohM3bZJ/n3U+8y1MDcBi5L+rwOO5bgOcwLhcqNuP4jR10/hQZPuEwEpMbu6Ufx+ALSBGCIcnd2KLgCEqiJ1HaREP9I0vEMRiO37MXUdo68fevuSThKYIauXNLfuRAWQksXP+jHDEQzTGjkp/3oJCej2acbufdYf4TB5f9xkbYuX2QHawcOJY+nqzup95nrO/BzQKIRYKoRwA1cA9+e4DrOO4vVy7JPrWfH3ILFzT0aNSPSoitA0UFWoKgdAbHoFGZ2YMCv5+ailJSheLyiq9eMAgNHfj9HWbv0jpfUDGF3dCYElXREsZaKHTj7HDIXAnJtToJwKs5RSBz4BPAzsBO6SUm7PZR1mGyU/n10/PJEPfPABPlf+OL2fGqTglS58r/iQuo5aVkrzRZUAqA1LMLu6JySYSkE+sq4SpagQbXEtaknRDN+Jw1wjp0tTU2Ehzpm1ulqGVlWhBXXUbfswg0GU/HzMgYGEUkzqOorfP9xzODhgzZn7ZfecUIA5AHpTM66mZgDiekpzYMD6Q0prfgeOIE8DtayU6PFLMDWBGjNxH+hAP9qU8VitqpLY0iqkW0EqgmihhqkKCl9qI9RYhn9bE3pL64jzhMeDsWENylNbrGG4ECOH6xNFUUGaqecnj8jS92W6j6ld2cFh8ih+PzKmgzStpbekj1N4PJZ+IL5NURGKSDRsammJ9XHrOkZPD1pdLTISxezpSRwTP0+rKOPYWxtY8o79HOwpobawj50vL2Llr/MwAh5QBEOVHqQCgaYhBiq9tJ6uEisyQJGsaDhGX8TL4SdrCJ8YoqbUx9FDp+Lu1FAiAj1PohcYICR5FUEqvrqInpAPnztGVNcYCHmIdPsyPwT7HE0x6esJIMMq3tIhGso7iRganYMBNNVEVUzqC7vRhNXcH+ovIRhxE/nUM6M+X0eYFyJpggAgNC31o58sSb2O4vVaZakqQgiU4iJkLIYcDCLqqtHLC9D6rDVXvdCHGowgwjHaz6rA22ciTNBCJq6+KMIwUQcjdK0vRZgSX4fOQJ2LQJtBqFxFMUCNSrpOEBhLh3C7dfS9q/Cs6mNpSTfVPhebW5bT1+dHhjSqlnShCkmJ1kR9Xhcbig9xsv8QS5b2wCVQqkoUIF9xo6DQaQxx1PAQEDp1GniFhoaKKhRYl3T/J0z90WWTDQWdo+5zhHkBorhdEF+OsVHLStFb24aFWlGt4xQFUVvF4bdVoedJlKhgyd/6MXwuXF3BhADrxX4QAqkJWk/yEQuAtEeB0WITrSbEJY0v89aiu1nuCtNhCEuoFBgwJQNSo0jRMST0mS7qNHsqAXQZghpN4BUaYamTJzz0m2H8igsNFdM2RVAQlpCdkXbDdaP3VsP42RMLcubtn0OpD/L2lS/RFingmeZ63tO4mbcVvIQLNx7hmubTnz0cBdg0UfLzATCDoVHnNcmKLLWgAFQVo6dnzHKFy42SF0DqOqK6gmhtISJmYnpU3G1BREc3RGPIJdVEynxITWGoTAUJBYfCNL/WDwLUMBie4XKHqkz8iwZYU9HKcfktqMJkwPCSr4Z5fd4OipQIfmGQrwi6zdQ6lSjgEgpeoc3Lj375P99Pw1VbQEprdGEvP6kFBcj6GkKLC+ht0BhYZhJY3M/GmsO8sfgVTvU0U6m68Svu2b4FNpx3lOe3hh0F2KRQVBSvB+H3QVEBemUhuk/F9NireSa4+2MM/H/99A15MV4swt8q6Vkj8bYr6AGJ4QHXoEAqoEZA90uiFToioiCKooh4UV1uPF3Dyo5IqcHZG3ZwdtF2eg0/Z/n/znKXQUyaeIVKm6GzPVph7fP9hXJVQ0HBI6zX2WUOUax4cYnJrjUrwPBcr2wBLVUfjA1S/+vhqYIZDif2Gf39sK0f7zaowvoBOOpy85uiDdxSU0ZocT69yzQGGg0qGzo5t3oPr83bxSp3D9WqbwrPOvu8+nrm5Pmkolr2toDi8yLy8xlcW8ORy0zWNh5lffFhVnuP0eDqoE7T8Qo18dIMKekzo+QrGj7hpsUI0WeqVKomZWoge/V1yAqv2XY5eRcfmZ7eII4QKB4PSnERRlUpfavyGVykMLhMZ8mydl5fuYfX5O1mpauPStVHTBpZ69VfXT2z3aOmIAQiL4C5uJIj5+XDun58nigl/iFOLTkMQL33KKs8xzjZHSZP8aYV6rF/kssk5QXVaXnUZf9uHLLAoBnGvKUCqWfJxl1KzHAYs6UVWlopeIlhVyAheNZXxOaS8zGqihlYGkAbkjS/XuGFd/6QQmUULXcWmPPCLDxu5Knr0ANWVU1NEKrQCFUKpGoNXWOFJijgLg+xpqqVk4uOoiT5b7gUnTXeZta6O6kcd0iULsgO851vdWyg6MEd5MQIU0rMUMjSkTQ1k/e8tbnxiXz+5/VncH3VSzN26TkvzLUrurn9zp/jT1K4eIRmaTUnTV72KuYwLzCkyZ/veQ2L+p+e1XqYAwP89b4zuP4jMyfMcz4GmF8olKkB/Io78TM1QXZ4NfLnYBFLf98829UAYOndXbwQmTkvOEcqHBY0n3/0XegHD892NQAwdu7jA1uvnrHyHWF2WLDsiQVp/O3QbFdjGNPAf1chIXNmemdHmB0WLO/c8kHYPLc8bIsf3c9t/UtnpGxHmB0WJD1GiLzfFM65QAJGRwfXP/LmGSnbEWaHBclnm88j7+GXZ7saGWn4Y5gj+mDWy3WE2WHBEZExttx2wpz1B1c27+C6Q2/JfrlZL9HBYZa5sXc51X/cO9vVGBUZi3Lg7kZiMrtTgCkLsxBikRDin0KInUKI7UKIT9vbvyaEaBZCbLF/Lkw658t2WprdQojzsnEDDg7p/OzPF2J0dMx2Ncak9r6jPDKUXRv+6fTMOvAfUsrVwEbg40KINfa+H0op19k/DwDY+64AjgPOB35hp6txcMgajw8pNNw2twUZQD/SxHWbr8hqmVMWZilli5TyRfvvAaxom7VjnHIpcKeUMiKlPAjsw0pX4+CQNd7/xAcw9uyf7WqMj5RU3+Wmz8zeOnhW5sxCiHrgJGCTvekTQohtQoibhRDF9rZMqWkyCr8Q4hohxPNCiOc7uubW0oLD3GV/bJDGm/SpB9XLMXn/2MUPutZnrbxpC7MQIg+4B7hOStkP/BJowIqg1AJ8P35ohtMzPnUp5Y1SyvVSyvXlpc5I3GFiXLXjapRNr8x2NSaM0d/PH+57bdbKm5YwCyFcWIL8eynlvQBSyjYppSGlNIFfMTyUdlLTOMwYPUYI9ddl2Qk+kEOW3tPLtmh4/AMnwHS02QK4CdgppfxB0vbqpMPeAsSbyvuBK4QQHiHEUqAR2DzV6zs4JPOllnPJf2j+9MpxzJf3cM2Oq7JS1nT8mc8E3gO8LITYYm/7CvAuIcQ6rCH0IeAjAFLK7UKIu4AdWJrwj0uZ5YU2h1clMWmw+baTqAjOrs/ylDAN5J3lRNbGph0kccrCLKV8iszz4AfGOOebwDenek0Hh0zc1LeYmrv3M78G2MOUPryfW7+8hGsKpzfrdCzAHOY1hjT54T2XoLe2zXZVpozR1s53Hpu+84UjzA7zmoeG/Cy/ef7rURv+EJm284UjzA7zmk8+fPWciSQyHdTndvK5I5dOqwxHmB3mLdujQ6y4JTjb1cgKMhJhz10rp+V84Qizw7zlbc9dg3xx52xXI2vU3t80LecLR5gd5iXtRpCqX3vmXCSR6aAfPson/33llM93hNlhXvK+fW/H8/jcjCQyZaRk8R9UOo2pTR0cYXaYdwyaYbpvWoyMRGa7KlnH98QO/qf97Cmd6wizw7zjCy2vo/hPC6xXtjGDQR758wYMaY5/cBqOMDvMK2LS4N+3n4wZXBha7EzU393BC9HJ6wIcYXaYV/y8t4Ha2/fl/LpqUWHOrmXsPsCHtr530uc5wuwwbzCkyQ13X4jR1p7T66qVFTT/pga1tCQ3FzQNAn+YfOYLR5gd5g2PDvlY9tuWnF+3+V3LefyUm+h9w4qcXbP4kT3c2De56znC7DBvuPaB92HsO5jTayr5+Zx91XMUq34G39UHSm4i3xhd3fzswfMndY4jzA7zgs2RGCtv6s/5dXvffBzfqvoXALesvRXlhNz1zg13Bdkfm7jzhSPMDvOC9z73Acxtu3J6TcXrpfBDR8lTvACc4nGz/51FObu+eGk3H9//zgkf7wizw5ynRR+k+iZPzqNuDly0lt813pWy7b0X/zNnmm0Zi9J+1+IJO184wuww53n/vnfm3HRTaBrKNe2UqamOD9eVbKHngtU5q0f1Xw7zp+DEtOjTjc55SAjxsp2G5nl7W4kQ4lEhxF77d3HS8U56mjlMyIxyU18VjY+/jzU/v5ZVv/4Y5+64hLsGZy5B+HgMmmH6fr0o56abkXPXcfvq347Ynqd4GXj7AIhMEbOyj958jC8+8Y4JHTudgH5xXi+l7Ez6/0vA36WU3xFCfMn+/4tp6WlqgMeEECucoH6zz7ZomKu3vg/vnUUU/+MAy9q2DO9UVG6tPpP/vXgJl1z7BF8q2zrtwHOT4SutZ1H81x3k8iMRmsbgJ/up0/Iy7v/Z2ju4vv7NOQuK0HC7QdN5g6PWJ85MDLMvBW61/74VuCxpu5OeZo6xLRrmE5/+FJVvO0DBHc+ONMgwDfTmY5Tf9Bybzq3m+N9+KqspVcYiJg2e+N2pGP251WIbp5/AXSfcPOr+M70xmi4ZKxNTdtGe3s5nJxCFZLrCLIFHhBAvCCGusbdVSilbwMpHBVTY2530NHOMPnOIq7/3WXz3P4eMjT2MVoqLkaEhln/tJdY98Kmc1O+mvsXU3ZnjvFGKyuGPGyx1jd4LuoTK2iteQXg8OamSjETYde/KcZ0vpivMZ0opTwYuwMoCOVauDSc9zRzj/QcuoerXL05IS2x0dWOGQpiRCGu+cYzjn303g2Z2MjFkvJ40+fEfLs151E1xyhru23jDuMd9veYBODF3a86L7m3ioSH/mMdMS5illMfs3+3An7CGzW3xrBb27/i4zUlPM4eISYODdzZihocFUi0tofXTZzB02YaRCp54RA8p0ZuaWXTlfk66/TNTctWbCH8OFrHs1qYZKXtUhGD3Rz2sdo8tNABLXXkcvGzsOWw2mUgUkumkpwkIIfLjfwNvwkpFcz9wtX3Y1cB99t9Oepo5xCNDAaofSbVzDm1s4NHPXc+3f3AD6vKlY55vhsOs+MEB3nPo3Bmp3xfvvxL90JEZKXs0lLWr+cM5v5zw8Rec/xyKf3zBzwpSsuQOFZ3RG8/p9MyVwFNCiK1YQvk3KeVDwHeANwoh9gJvtP9HSrkdiKeneQgnPc2s0WkE+cY33j/Cztn395f51JE387e+dYjQ+ENovbWNpu82siXLy0Z7YkGW3zm9GNKTRgh2fTSPDZ6Ja+q/WP44sdNWzWClUvH+awfNsfxR908n2foBKeVa++c4O/UMUsouKeW5UspG+3d30jnflFI2SClXSikfnOq1HabHpdvfQ8ldL2Xc98r9q4iYGuFV1QiXe9yyfPc/xzVfu27aAdyTuWLrB3IedVNd0cD/veGWSZ1TreVx4LLxn1G2MINBBrtHHwk4FmCvMn7UU0/BF9wpc+U4ZjhMzVMhzincwadvuIM9vzoBbVHd2AVKSfFvN3P2Q58hImPTrl+7EaTwFwU5j7q5+6NlvMk/+fp/4PWP5zRwgbtr9FGQI8yvImLS4LafXTCmw4I6GOG6zVfwH/dejdQFbecvQi0rHbtg02D5bTH2xaafuu2jBy/Dm2PTTW3JIr590R1TOve6ki30vyF3Q+2x8k9nwwLMYR4RKR7HDPFAE43frIFjhzF6egAmZH2l/Hsrb73ts2z/4M9RxdT6iJAZ5diNDRSGn53S+VNl34freGugh6n0bXmKl/a3DbH0HpFzR5B0nJ75VcRBPUz5lrGHkubAAMb23QlBnijy9BPxrZ3cOel85tjZlDy4Z1plTBatqpLPXH7/lBsggJ+tvwNtyaLxD5xhHGF+FeEXEkWfud7j1rW/AeBgbHDSBiUhM8qLv1yH0dU9/sFZ5Mh7Gvhg4fSWwF7vC9P85nF0CznAGWa/SjCkyW29p+A93DsjTguuwx3c37+Oy/75OmqfkHSvUokWSL5z+e95a974ttW39S+l4qGDOU2YrhYXc/lVT+AS07MydAmVFe/cTf+NnlkNzO/0zAscQ5rsiQVZ/eT7+febV2DsnpkwtWZXN3fdeg4rvvASBc8cpv72Jnwdgh996V3c1l825rkxafCjuy9Bb2mdkbqNRusVq/hK2ZaslPW9xfdhnpw7RVgmnJ55AfP1jjXc/bvXUfdwNw3796HPYOB4Mxym+vtPI8ESSiEo2l9F3xKN333gItRb/sq787synrstarD07p4xbJuyj1pQwBkfeDFr7pyLtTwOXO6n4ZmsFDclnJ55AWJIkzfufDObL2mg5vqnMbftyn0GCCnJ294JArS9x/jB998xaoCDD217L+Yre3Nave5L1vCd6ieyWuYnLnwwd7G1M+AI8wLkzXsuxnVlBP3w0fEPnkk6uijaryMHBqj4/TYu2vn2EYe0G0GKf56XWyMRRUW9qj0RqC9bXFO4J6extdNxhHmBcddgIfLTBTnP+pAJo7ePvJdbMcNhzGAQ+YMKWtLMPq87cjGep3bktF5i7Sp+uer2rJfrV9z0vn0wZ7G103GEeQGxMxriJ1+5Iuchacci2fPJ+49tfOXYBYn/IzLGnltWYYZCOa3T/ncVcKI7u71ynN+ccgvqqoYZKXs8HGFeQFz67EfJ+/MLs12NUZGRCP/+x/GJ/2/sXU7FfblNAqcWF/OhCx6bsfJPcascunxs7f1M4QjzAuF/OlfR+JX+MW135wKLHovSbljKuJ/87UKMjo6cXr/3vJV8vHjmbL9VoXDBW55FyR/dVXGmcIR5AbA9OsRjXzoL/cCh2a7KCITLjbZ0SeJ/96ZdXN/xGjZHYjTemltrLxSVyLu7s674Suf/Vf6b0OtyF1s7jiPM85yIjHH5bz+L58HnZ7sqGZGxaEpIWjMY5J5tJ/OBl96HsTO3Q2zl+EZuOn5kLOxsU6j4OHJZLlfNLRxhnudcfeg8ln33lVn32JkMJU+5Kb05kHOf5f1XFrMuRxE1v3zGA6jl5Tm5VhxHmOcxERlj380rMQcGZrsqk6L89q34Hnoxp9dUy0r58MWP5Ox6V+Ufovfc3Gq1pxPQb6Wdlib+0y+EuE4I8TUhRHPS9guTznHS02SRx4byqXw4t0HvsoEZCuVcUddxyQo+VZy7JTu/4qbjkpkLRZyJ6cQA2y2lXCelXAecAoSwwu0C/DC+T0r5AEBaeprzgV8IMU13lVc53zt4Hvqx3DonzEeEy0311QdzmlYH4PMnPZLToXa2htnnAvullGMl33HS02SZY5tqcj7vnJesW8nnFj2U88telX+I/rOX5ex62RLmK4DkIEqfEEJsE0LcnJQF0klPk0UMaVKU26Ac85Z97wxwuif335FfcXPszdMPcjhRpi3MQgg3cAnwR3vTL4EGYB3QAnw/fmiG0530NFNkUEbIPzp7jvDzCV+7QpsxxONDCn3mED1GiIOx3MTl/tbGP6HV1uTkWtnwZ74AeFFK2QYQ/w0ghPgV8Ff7Xyc9TRZp0sHdNDNRQxYadb98mSt3fxb/kSB9q/KJBQRDF/Tz2Kn/h0sINvzzk2x+/U9HJFbPBpcE2vjFqXX4mmf+U8+GML+LpCG2EKI6ngUSeAtWyhqw0tPcLoT4AVZ+Zic9zTTYFa2Enr7Zrsa8wBwYwHffZqQQFO32YQ4NIV5Yw3nnfIFwuUSrDxLOwjr9tmiYGtWgTA0waIYZMHU6TI3OEzQW/Xn69zEe0xJmIYQfKwXNR5I2f1cIsQ5rCH0ovk9KuV0IEU9Po+Okp5kW24YWIQdzHHBgviNlwkNLbtlB7W4fZjiC8dq1vEH5GM+f8auMpp4xaYwbJywiY3x+/9v41rJ7uaWvgc5YPn986jSqGzsIl5kIz8zHBxNyjlsOrV/rlZsfnv0wpnONTzSfxr7TjTnvWDEv2Hgihk+j+doYmmYw2O3npMbDvLZ0L8vc7VzkH5xUKN7t0SHyFQMF6DRc3N23nhfee1xWXFM3yb/TL7szBj93YoDNUw4MliKNlvEPdABAaFpKwydcbpSli5B+DwzFiOW58D6Rj+kCd4kk2qDRHClmlaeFFiNEk+5DESateiHdRh6LXF3sj1ZyYWAPj4WWcbbvAEtdeRjSZFe0ku98590U7wqhdQwgWzswB2beYMUR5nnK4e5i6qSjP8yE0DTMDccRK3Dj6osSqvHSfL5J3QMKasRE9ym0nCVQKsIYugI9bhASTxcEjkkKDpu0HqqnZ2AJW1rWoYaiqF2WyawcCkM0hqytRBkM8YeGC/BsO8LNZ15G23qVvKNQ9tIgJZutyH65nEc6wjxPCfXPrBvfvEZViZR60EIGrqYu8rd2s+pBMMPWnNUDNN4nkIYxqoNKcq5FCSPjefdY0US1Q0cwAN+fO6j/c7ZvZHI4wjxP0dpzl0p0viEjEbx/sRZKRtMoyNx7KM44jtfUPMXfMk4COIc5hxIIILSZ6z+dnnkeEpMG/vYF2LXMRxQV4dKseXowiFZXixwaov/1jYSLFHpXglGo4y0Oc279Hjb98lRKbp6ZSPmOMM9DQjKKvy13Nr8OqShrV9NzQiHBGgVXv2SoUuBvkZT9ejN6TQlaex/BSpVQtcTfIrjkjc+y3NvG15+8lCWtM7eU6AjzHGN7dIhfdZ3Fxrz9vD2vK+P6Zpth4mkZcEw5Zwlz2y6Kd3so0nVruUsk5Wbe/DI6UP0PD92nlFK0a4AXbsjnqbM3UudT8O/tmLH35ghzjjGkyXMRyUavSqcRpFjxYSJ514HzePmJRkp2SJAw8GkvV+Q/lbEMFxKpOeqOWUNKzHA45f90jJ17Kdy5FykEUko8Dz5nbZ/BajnCPAPEpEGPGaZY8fJgKD9hQbQ9OsQ1u96N+HU5fKiDzs2VRGpjqL0aRTsFDX/abeUnFoJ/nLMeFmcW5gGpIUKOx9S8IIcWlo4wZ5mf9iyhLVZAmWuA7YO1PH/rWj59og6KZPltMQqe34XUj8B9CgHjsLVGYr/wRKstJY2/ibH9jUMc5/bN2r04zC8cYZ4khjQ5oodY6soDLAN7j3BxRB/kQKyAGlcPP/vzheQdgaqHmqg4uolKRSA8HsxQCHOCLbWy6RWu2/8OrlvyKOXqAA/2r+VjJc/RYSh0GHmIIadndkjlVSXMvx8o5WTPUY4Z+bzWG83oCTNohvEIF9uiBqd4LMOMFn2QMtWHS6i0GCH2xoo59/H3UFgU4sqG57jz4Cn07ith8cMGG779HMv+0IO5bVfCYEGaTNohQhoG+3fU8FvvGWzeV0/p4x5uX3Y2H7v8QSKmCzmQG+d6h/nDghZmQ5oMySh5ipdOI8jFgSFMCV7Rh2Ib7IVMS6hdQuWG3lp+cselnPXml3ixfRFPn3QHCiIhyH3mEBe9+GECdxay6qGd6Mct5Z6lb8Qbkax8bDci4OefPz6dskPbJ11XtbyczguXU7J9AOXgMYyePqQqafpBI6ueOojR3kGp5uKv/zgHJWoi+rdk+Wk5zHcWnEo0Jg02R2L8O2xySA9x8Y4r2B8b5IDuJk94eDhUy2LNj4kkZEZ5JSZQ7IhG7y04SCxPsusbJ9B5tIj/aNnITf11vPvgm/iv9hM45xufJby9iHCRgLIS1G378XYbuAZNhN+HHByk+DfPYPT3A6AWFSJc7tFTfAqRyElk9vah+0Ht6MPsH0TxenB3qbgGDYylVYQvOhVx3HLcbYO4DrVbyyEODknMG3/mZAfxPnMIv3DTaQyRr2h8qeVs3lG6ia8fuISBiIdITEM3FNx/L6T68U7Kbmpl05F6vr7uL+yLVPLR4hd4OLSYelcHV/3jGn549p18Z+/53LD69/xP00UEzwuhFOSz87u11FT0UuAJk+eK0Pm1enzbmwmdWEfvtYO47y6m9IHdhDY2EPMr9C1TESYs+vVO0DR2/2cDpVsEJXe8iNA0hM9L9/mNtJ8mcfcoRMoNPnzmE/zx4EkU+YdYlNdDW6iAmkAfi33dnOQ/TI3Wg1cY7I1WUKX18sjACfzp1rOp+enzyFh0lt+Ow5SJN/ATiK4qXG4G3nIy+X96kWejD47qzzznhXn5CX75lXtO4t62k9nxfD3aoiDKVqs3C1cZeKpC1P1Ypa/BR9GeENr2g1BbSc+6UtrOlLxxwzb+fc9J/PojP+XKf30YOaTha9KofjpCy+kexMl9VP/Ug+lW+O8bbuSfg2v4+xfOQhuM4eoOEflJmIPNZVRX9nJZ3Vbu/eYb6bo0hNnsR1ZEEG0e6k5sJXpTFa6gSetpKkbDEKLJR/0pTRw4VkZVeR+X1W2lUuvjHP8hOg0XBoIbO87Gp8a4rOhFipQhPrPvHRw6UMHaVUf44/IHeCECS7QhqrU8YtLght5l3HjrRQytHWL5B3bMeOQKBwslEEBGY+M2npOJJqL4/al5qZMNT9LL1TSUZUsw9uwfMzjBnBdmX9UiecrpnyJweBD5yh7EquUI08TYsQfhciPcLmQkQttHNqCFJGV/2IrweDAa62jfkE9/g4npN1AHVUq3CPqXCcKLo6y8IcLua92obpOGnxoce20eyhk9BPcX4ulUkAoM1ercev6N3Nm1kSXeLhRh8vJALZ+tepSbu16DLlVqPL0sdneywt1GiRqmw/Dx3NAyri7YgUdouISKiUlY6mwKF/Dtgxdi/qSSI28xISbwHXVR8WKMpnNV1JCg4XcdHL68Au00qy5L1h7jxsbb2Rsr5YeH34j5jXKkIlCf2Dq/Y2YrKkIRsxMpRVEn9ewUvx8ZjU6srpMseyIITUMpLMDo6WOT8cjUhVkIcTNwMdAupTze3lYC/AGox4rz9Q4pZY+978vAB7GWTT8lpXzY3n4K8BvABzwAfFpOoCUp1MrlRu+Fw62Y3YKpBQVIw0ApyEdvbUNbsgg5MIgMhjAjERSfD6Wqgh2fq6DuMYG7T8fUBJ0nuhlcHeV9pzxNU7iYoO5mS0stp9UdJmaqrMlrYW+oghda69B1ld+ccgtHY6X8a2AF971wElq3xtvP+zevzd/Nn7pP5tFdq3nP2k1ckL+NKx+4FlEU5W3HvcT7ip+mSLGcIcISLrz5C6gRWPK7wxgdnfS882S8PQa+Y0EwgR37EG43Mhq1huQuDWPFYrS2XvZ+tI6l9wURhoSXdi6IUEHpkT9mtIzRer1MDYoQqGVlhNcuxv34VmtfmoAqXi9mJDJqmZMRZrW8fMI5qpVAgGcG75+WML8WGARuSxLm7wLdUsrvCCG+BBRLKb9op6C5AytTRQ3wGLBCSmkIITYDnwaexRLmn0gpHxzvBgoCNfL0xg/C/qMIrwfh9YKiMLi2BqlBX71GLA9ieRLTK5FFMUrLBjindg/n5m+nQLHM7pa7wnhtrXWHEaFI0Xgl6mKlK0Kx6ucTzaex6wvHoQ3G6F0ZoOsEQcF+cAUlRR88yp6XF7Hicy+ilBTR9O7lHP/WnShCsv33a/B2mxTuC8KW3UTPPoGuEzxECyFSYqCVh3FtC1B/+1FkMERwYwPev262Xro0QSjWB2XYH4BQUNwuUFWEqmKGQqhlpbResoyqh5sgEgVNAynRcxC+dabou2ojJX/bjdHXP+xcPM63qOTnDyfJU1TUVQ0YO/ei5OWlJM9Ti4sxenqSTlRTjHPGRAgQCvL0E1A270DGovS+53R83Tr+Z/djdPdkz6rLvlZG4R+lAZr2MFsIUQ/8NUmYdwOvk1K2CCGqgcellCvtXhkp5bft4x4GvobVe/9TSrnK3v4u+/yPpF8rHc+yWnnajVeytKCL0wv3s6lvGS1DBRxX2MIny/5Freof4YxwMDbIUlcel+09j617FuPq1Njynh+zJyY5we0iInU8QuPbXWu4ZdvpqJqB2Btg2R96EFKCbtC1oZyivUG0tj7aX1dD+f22qSWgFhRgNixC7R3E7OwG00SuqAegf2U+rqBJYHcXQjeI1hShbtqBNAzUvABDZ6zE3R1G3XU4ofW2H7L1W0q0ulqi9eVoL+5B1FRi7DtI73tPx3BDxZ/3ILxezL7+eZf9ERi+T6EkGjPF7UIE/KBp1vq5qlrD2pjdY9ofu/B4kNGo9YxqazC7ujHD4RFz1RHz0WSSek6ttgbp92LsPTB2lV1ulLwA5sAAwuNBKS9FP5SUsC8ueGPMe8clvff3+0FVR7zjmQjoVxmPjW0LdIW9vRar540TT0ETs/9O354RIcQ1wDUApTVu+oa8PNm2nNPX7+cfL64hcFiDezxcfOnpfO9jv2KVu4cBU2W5S+MTTa/jqQfWcu07/sbBexsIeKF0h84z7/Dx4Yc+xLoTDrD9yeWces5Onju6mJp73RQ8fQgpJZHjFqE9vgXF7aJoz36r9wNKDh5OMZA3+vtRdh/EiOkgTYSmoXb3I3v68JU04nlhH+ZQGMXvRzlwCCkEwu3G6O/Hv7MVGY4Qb0QTQ8Wkj0BvakZtbcM0DNT2LpT8fEq29WFu3YkhpfVBFRTAfBRmKYc/XCHQFtVAJIre2obi9SIKC6zpUiSCuqIBNBVjh5WHR/F5EVUV6IePgktLPEO1pBi9rSMhDGYohLZkEQiRKnSQIjB6SxtKwLI3iAtsSo8er3IsitEbs3rsYBAZiVj6GlWxRoslxUifB2P77sQ5E50CqAUFGIPBEb2zGQqhBCYXlD/bRiOjpaCZcGoasNLTADcCrDzBK59d/1sUFDqNIe5ZeYzFp/Twr1OW4XkafnTRpTRfUIFUIVwm8bUJlv5+Lw/+YgW1ZR2Y+w6jlhbzv++5ilU7dhHJz2e50kTbo0tZvq8No60dXdcRLjeeo3mIijLMgUHUkmLMvv7RX4gQoAhkTCLA+sAA7752ZG0VyrF2ZG0l2B9HvJcJL6+ge7WHmr81IzTNGmZmegb2dc1QCKW0BGPLDmuHojJ48TryH9kx2uObswhNQ5rDjZhaVWkFxwuHUVc0WD1kW3vieGPP/sTfit/quTGtxiBZSPWWkZkwzbYOlKqKMTXRituV6PmE24UoLky8L2tjUk8rJZZyw3o3QtMww1GIxhCDweFpko00jAnNn80xtN9mMJgyYhMeD4yRJXaqwtwWz1xhD7Pjb2C0FDRN9t/p28flSEc5Twz5OcM7wM+7T6f5kcWIRwqpWOJjoBZo66D6Z4cQxzVivrwHTMPqRYWAzk4A9NY2RGubtb23D8XrRTl6DD3pQctYFITA6OpBxqIoimLNW8lsimkGhwPQS10kXrzZ0YXo6UMU5CO6ekmPB+J+dheV3fXI3n6wPwDF6x12qVNUzDNPZP/b3az6xj6Mzi5kaGi4AGliuAViHhqNJD9HaVojEMBSOgll9CGqECgVZSN72TEww2HMw0fHPEYpLsK0GwIzGMQ8kJZUILk+9rqw4nFhRiLD92IameOJSQkTyPGQaSkrZZQgZSLUkOLzzogw3w9cDXzH/n1f0vYRKWhsBdiAEGIjsAl4L/DTiVxIicGnf/thELDo4SC1zzyNBHwvWWpx0+NBGgZtZxajrt9AxR+3W3NRKdGW1SM9LmjvgphuDY8DAcyhMJiGNfcdCicUMPGeQPF6JzUfjb9YoWkoxUWWYip5Ppz0Ys1gELWpHRmJJOZ1Znj4pQuXhjoYZdEjroQQm8FQSlnFj+5NnW/Pd6S0PtykXih9/whBHu3YtPMUvzWMzjSHNrqtXlgJBFKyXWRC2M4ywqVBmgCqpSVjKsYmm81CeD2pt2F/X0bv2OmIxhVmIcQdwOuAMiFEE/BVLCG+SwjxQeAI8HYYNwXNxxhemnrQ/hkXxZAs+7/9lnJkaAgDS9hQVcxgkO4rTqb8scNEiiGWDxXu4YTa+oFDVqtfWgJV5dDfn9KjGoNB1LwA2OcYnV3WXKhhCSTNfyCt9xwFqespGmbh8YCRmnVCLS2xFGmjDL8Ujwdj6y78+/wInxdCIdSCPIzePqtuLg2js2sij27ukzwMVVS0yvLEkFloWkrwu8Szt0dA3e/bSOm2fuQLY9vBjyWgcQEzh8Z+r2C929GmXEZPn63QG6UntqcWIxqNURqk9I5EuNxIPTaucm3OG40UiBJ5mjgXecZa1FAMc8sO1NISgqcvJ/DcIRgKI3V9zJeWCSU/HxmNoixbDM1tqUqIDHMd8zXrUJ5+GbW4EBkMWa2nUDB6exPzmRGtrxCo+fmJXlTJz0e4XQmt+MhKqShrGmHfIcxw2BpFBBdwPilbMIXLjYxFUctKMfsGkHoMbXGdpWneuXe4AWR43q0uW4yx/9DwBz7G/FSrqkRvbUMtKsToH8y5sY3QNGsOLRSEqibm78LjSQQCTGAHCEzR5CsqitcDpsmzsYfoi3VknGPNC0cL4XKjbN6O3GkNg4XXS94rrci+fnBpkxZkAKWoEKV+EcbuA+DxWL03oC5fOuJlq+XlqM++YvWK3T3WfGz5IkRhPsJtOVKkC7Li9ydaYa22Bm1RHTIaxRwt2ZvdSouOboTPCkgwQpDjo4z5iLD0CmplheV8Amg11QCoZSWWrmBgMLEerB9pghZLFZP8zCzNvwkDwZSeSigicxhbRU2MZIzevkkLsrp8KWpBQUp5KbelacPbFDVjHZJXK5KH0DISGfmOpWnpUkxjuK6mgTQMzGgMqY9e/zkvzMLlwjhtDdFz1tHysVNQVzci8/zoh45gRmMc/OQq1MZlo5+vaajl5Yn/lfx862MyTTo3lqNVlCH8XogPoTq7rWM0zVK8eL0Qi4JQ6Lry5OFyDrWAolgtaCYNiKomfkeXVaJXF1trp9FR7HultF6ermdcHrFuRrFa6fx81NWNYz63uYZwuy3hNYzE89KPtaAWFKC3tlkfrK1YEpqGWlKc0PSPGPFIiZlp/qiqINKE2pxGcj0hoKfPsii0vdu0mqoUjzWp6xNvIEwjoy5GuNyI9cfbBcqM9RWqNfIQ2ujZKOe8MMcK3bj3t+I5NkjNP7oRPf3Io8dQi4tBmhTsTxpmpWl41YIC1EW1yHDYUk4dvwo5NISMRTH7Byh/ohkZDqMfPpr4cMzgEMIOShAP3GYOBhFuF+X/ah0O8dPZZc3JzcwpTuIvTUYiKP96CTa/bCvCho9VVzeOaMlHHYKDpanv6LAahPb5NW+WkQhGVTHmoqrhj1VKzEiEyAXrk4xJLPNKs38QrarS2uRyW42qjXC5LQFPEyppm1iOJ7zC47HWoTOg+P3W3Na2zjK6upExPfE+9eZjaJUV1pzek6qowjRQ8gKJUUh8BDLus4lF4eW9ifsfUV9NA9PuMIzR46XP+TlzXnGdPL3/jBHLBHGb2lEX54VAqKplJK/rVoJtzTU8X9E0y9JoAlpG4fEgY/rwXGY6c67ktcvR5nnZsCiagyheL0p1JfrBwynbk/UN8fnz8M60ZyBE4j0qfr+VP2oMXUdGhEAtKso4AtLqF9N9eg0Fdzw7Yl9ChzHGddTVjZYOpr9/xHFqacnYtguAVl2Vcd08/lzGsgCb8z2z0htK/fghZeiUWBbyeNCW1Q+fl5cHJ63G6O+35tRSJtaS4+elC3LycDwZGYlYiotodNKCHO9dEv/blknx+8h8QZn6Ow3F683Ygs91xJI69PKCEZZNSlFh4u+4ICteL+pxK1F8aQEN4+8RW1M9xvtI7xmFpllTsvhSWNL2eM+vHzpC4R+fz1xg45IUAVXy81Pn04A83DysHU8W5OJiet+4AnVx3ejBKshsAKOWlaIU5o96Tpw5L8wpjCVIhoF0JS1lDAygHutCLS1BW1RH5IJTAat1VJM+nmTkwEDWhES43JY5YWtbynb90JFp97bSMC3F2zwhLrzm/kOoe5ustfbycqLnrbdMNhNOJmnDZiEmZ9aYLDzl5dZyThLSlAjDTBnSC5fbskpTbFFQ1JTGJT6U1upqMbfssOetmjUcLytB1tekXsMwLMFLE1ijp4eCP71kjUomO7LTdSgtHvewBRMDTOo6tKW6kuktbSBNVEUlsKsdHTAHBi1LmgyMtY48qagetqnnWEI7WUOCEfWZR4EJ4j2Vkp9v6SaEwOjowL/Hj/R5oKgAJTSENAyE241SUoR+pAn5ytQTlMtQaMTzV/MCiSG+8HgIn3si/mf2gWEMG+GYRoreIv6c9abmhMFRwojjaPPICxsGRDOnDkp+Z2pZqR0jfRSvqeQie/sQo62CJDG/euZxSFjIJA9jpcTo6k68RBmJjGtJk05csz32QSK1NTbGfkEyEkGsPz51KKioYw7BEufOt3BB9lqpCPitJSR7VCF7+5FHWzAOHLF0EYal7dUPH0U5YeWIKcqkLmkv+WjL6lGOXwVYRkJCc1lD9WiUwM52jJ6ehCAnFFqj+T77/VZHIASctBq1qnKEAlPqOqK0GLUgL7FNq66ytOHxb0hRMfsGUDwe1LwA6vKlw0tco3xnE9HIz29hHuXDj5vwJUhu+aYwjE73asp8kEQowxY9E3n48sWdyFg08dEqSdZrCwV1RYOljQ74Mbt7kKZMaHyNnh7L88iloZaVpDwz0dYNfmu+LFQ1oSdQAgHr71HefUrjKARSEYgjllWeWlyYWNNGSsz2TrS6Yee9UZcNAaSJ2dMzrKjbvh903XLdTEM/3GQZp9j1iaysYfCNa9BqrSG5EvBbDbKqYgwGMQ83DzckGXKLTZQ5r82OW4BNmKmGo0nXMksTxeMZ14QzayRrPidid5x+zhxHLShAlBSl2FirRYX0nr+a4qeOWmawtgZfLS0ZYbIa1+YKj8ea05oSs78fGdNRAn7LA2209flEIcPvOL4KkgjJM9aSYDpJz10tKkwZ6SVPn+Jmx/pJjWjP7xrxLakFBZih0KS+1Xmtzc7IWL1rmpGA4vejVVeNW6San6QtNA1Li2obfoym5c4myZZdWmUFygkrh3favZBaVmr9LraVIdKc8FrmbCE0zZor9/cjB4OJ3gkhEEWFFPzpJczOLhS7Z1I8HoRr5AglrsySkQgMhS2lkGGgHN9oKT+HhhLlJn6whrjq8qV2ITKpvKTVkLhxUHK9PZ6Rz9YOjaxVV1r193hGTNmEOjxUNm0LL+3FPRk7BWNgYITr5HSYl8IsNBdaXW3qfGWUYZeMRjH7+lOGU5lI90Iyw5Hhee9U5qhjzH3VosKE8Cr5+SherxUHyj5HhsMoXcP1iX8gZjyLhTRRKytQCwtQK2e+oZkO8d7PMvSI0nnOEsCyCDPz/InRjxmN2evG4eHlGSGs95a21pyISqLriLZu6x3HhUVKFJ8PtawM4fGgt7Qij7WlVysF4fMh4sYe8XcQiaTqJhQVtaIcltYiA/bwP4PpppmseLN/x82N1eLi1I4ozYgoY90ymaiOwrwUZhmLWtrFQnuNb7ThpqKiFBYgAgFk/+gujemtMmCZ3oWtaCEZFWb2i9eWLkmY+oHdE524CrnRMs/T6mqHTUjjRQ8GLU8bLPdGeVwD0fPWW8b0YM234stsQqCUFFnB/iIRa4kkOITw+zAHg8M+wXMYo6cHpTAfpaiQsscOWhtNSazcP7z0ZBrDzgU2akU5RmVR6qgJay02tL7e0o2UFiXs4wGrx1y6CBkOJ4a7Kb2fbT+dbL1ltLZjdHWj5uejFuSlLoUlKVPN7l6EbkJ7pzXnnoC3VcpzsJ1yJsNkeu55KcxxEvOc0eaN8WWGsiKoLBu1nNH0BtqiuuEGIwm1qBBl7WrLTheQ4Yj1cdimiCIUQd26z66kQd9FxyXWLtWCAqSuozYuTdRR7R7E+/RuzKGhRMMkk0xKjbb2hEZdP9ZqNWYHD9ueOHPbeES43KiVFZhLqtCPNhFtqLLWaIsK0f61DWlKpL2UI2PRhFAKlxs5GES+sH1ErDTF68U1EEMpKbZCLif1XsLtRvQNpgy7U5aEigtBKNawXrEtBGNWXDFjMGi7mib1hnGfaCGQegxj517MIcs8WI335hNlKvqpSZwzr4V5ooiBEPJoUmCTtBeQcc1WCMzSgmFBTT63phKls49oQwWEreFYPNAcgLHvIGYohHC5Mbp7cA8YUGw1CqHXrLQ+DI+W6AH0g4eHDfDthsnYvS+lOmYoRDySZ3wJQ6usmNPGI8LlRhqGFaTvFet+XLubMcMRus5fbgmhNFEryob1ADZSt4Q1HSUvD5Gfj/b8LvQmO/RSf3/iuclIBKOjM8X+O3kObXT3WNr0ilIUtytVgx1/9mkjscTQ2f6JO4QYg8FJC+iw1lqkWI+l33/8mHhDAgzP/UdhwRiNjGWEMaWhqJSWxU/6dVQVY+dekBK1smj0NclT1qD0WR+BuzeGCFmtuffRrdZIQElrR23DfuHWEteRMX3YvNHnszSftp29Vr8Yo7k1kQRgrqLmBRClxcieXoxIJBEjuuzJJiukkpRWkLy6arSAH1waRlML8pRVyH0j35s5MIASi6FUV0Jnd0YvpBFmuhXlCEVBb2kdNsv1e1BKS6zheHUFoqMbGRqaVIQZq2Ed6RuQ7i+Q/H/CQUQoVmOQdF/pxP0K4g2Sse/gmPVZMMKs+LwYE/2op7EcF38pamkJHG5Ft5dQhKoO7yvIg55BzKPH4IRG1Jf2oA8NgZRoVZWYwRCiqQ3TDrCuBPyYwRCK14OorUK2diCqKxBdPRg9fajxoX4olAgja3b1ILye1JBCcwjh8SCEsHrNDCGO4gEQUVSkYWK+sgulugrZP2CFtY0ZlptjWjQSoQhLYRZ31ogrxxQVtbQEGQyO8G837CCB1lRFQfi8GNt2JeKzqYoKisAcHESrrkrt2dOXDNM14pm8nNzu4fNFhmXSDPHBxowzN8ElyIW3zjwLpHv6xG135eIaRFvXiIwFwuW2QrSqtj+0oljRIb1eKxZ0craEeAB3LGVabHEZ4pmXLcFfughz36G51zOPtU6erpnONKJSVLtxjGGeuRb1uZ1j3+NEPnbbi04amV1WEwpKaVq22qYx4j7iDWlyzC+1uNgyBx0YmFYnMVHz3mmtMwshbhZCtAshXknadr0QYpcQYpsQ4k9CiCJ7e70QYkgIscX+uSHpnFOEEC8LIfYJIX4i5mN4SXttcTykbgUPNLfvxujosNY6k+ZHMha1hlVJ1j5GTx968zFriSXpo9CqKy0LMSkZWl2F60ArwcvXEzttFcb23XNPkGHsJZd0e+niIiv1bdyhYekStMpyhEtDaC5c2w9nKiUV00Bd0ZDwQ864nCOlFW6ooT6j44aSF7Aa17RgA8LtHk67a/f4MhpLvDujr3/agpxMXJE6soLjm/lORAH2G+D8tG2PAsdLKU8E9gBfTtq3X0q5zv75aNL2X2IFtm+0f9LLnLvYyxlabQ3K4pHr1SNspeM21vYLlkNDmCuXjDgvbrEUOW2F1fsKYRmG2MoR4fEgi/LB1mz7tllDU19bFM+BjjmryZ7M2qiMWKGUEo2SaSJjMUs/EItay1r2B54e5C+F9k470qpMGSKn2BeYBvJYm6WcTA8K0dc/clXDVnalH2sODFgNSGmJZas9FUFOE87EMlqSIjWl0bG9tcYscrxrSimfBLrTtj0ipYwP8p8lNSb2COzY2gVSymfsZHG3AZeNd+1cMpYllf66ddYQq6cX84illEmJ/UTqB6wE/KglRQktqtHbh3zuZauF33hiyvVkJIJ3yyEUj8cKqxO1ojBavroSY+c+5MCgFR20ohijqhT3wXZLy66O31rPChOpl73EZA4MJALjgzWXTjflNGzlkFJcnKq9T3r+MhxBq0hdfhSqSs9rUiOKJHzbzVQBFIqlOdYW1Q3bHcQ14KOYiZqDwcSy2ojbs5e+RvSoE/VlV1R7BDAcsWQ8s89sLE19gNSwuUuFEC8JIZ4QQpxlb6tlkulphBDPCyGejzHOMHICw48JDVGWL8nc0wmB9s8XLUWOYaAUFKCWl49YElLyAtaaakGBFXao3DbPtIdjQtOQK5eg9luGBvHeXOo6RmcXSlEhsiCQWFOVuq3JNg1kOALL6xFt3SiRGHpTM0Zb+5zNBjmhob8dkikRoC+OoqI2LiP41tPSGksXIm84h1R8WUtZu9pKa5MXQEaiKQZAMhal6FErtU08TJC6fGnCSCX5mtK0AhYYre0p1mSZiPeYcSsx4XInIm3G66bk51kNRJrzTLoBTNzcNR1LUy5SonmOx7SEWQjxn1jxsX9vb2oBFkspTwI+ixUQvwAmn55GSrleSrnexRhzVCGGPZXGIr0VTHdXBOjsTdgHp1UmIZDSMBEBnxXEwM7SGC/P6B+k4wOnEN5gBdoz9x1GLSoaLkbXkc+/ksibBNYHFrcBNoMha8krXkWPh8F3bETx+xF5AZSuXoz2DoydexMtdXzYOR1XwTlBstCYBsa+gxQ8eWBYyIWSMJRJEG8kj3ViRiIYnV1WT5kWI8vyGRYYzS3WtMUeZieI273nBezwUCMFJ30Ulh5RU+oxyz7c7u3jDXQ81nbycDmxHGXbCkhTjnCXVbxey5gmEpnYKCd+3oSPTEMIcTVW3uZ3x/MsSykjUsou++8XgP3ACqaRnmbU63s8E3Y1TEerrSE9oqbZ3485ypBJLSkC7N5UN6yllMHBRBlCcyFUlaq/HcH9xMvDw0TDsDIc2tZbwuVOiXIiAz4rdndBHkpZiR3k3hZUVaXwH3tRSksgEkkYN5hnrkW1o06ottFIejSTmUCcdBzmWSfN+HUAyxqro2N47uh2jRg1qWUl1vqx2zWsm4hFRwij8HgQJ69JhIkyh4ZQi4oSCkmtosyaj5aVWMou29VSLSocttNOj8BpG3PEy1B8PktLnmHoLA1jeDTh8Yw4RvF5UYqLUrZZSlDTshaMj3ImoB+ZkjALIc4HvghcIqUMJW0vF0Ko9t/LsBRdB+yMkQNCiI22Fvu9DKe0mfz1PR7E6oapnm4ZkaQNoWQkMuo8JnkOpzc1J8z/kuddcXtxaRiJXEFSysRwXC0pRi0rSbEuihufYBhIO8dV3KRR+H1Wb9PRidHbh9HVbZkqmjIhvHrzsSnFDJ8KoSUBjr3GN/6BM0C6hl9oGkZdOUZHF0Z75/CBSeGAwHZi8fsT1meAVY7blbB51lvbQFEtK7xgCFYuRaupHjPGthpPZmBPicZ8B8ma/aQeOK6tN4PBzHm2pZx0CqKJJFtPpKcB2rDS03wZ8ADxr/xZKeVHhRBvBb6BNfQ2gK9KKf9il7Oe1PQ0n5QTWOTOuM5sW9Ck5O3N0RKNWlycMJhXvF6kYab0BtqiOoyyQnhlb2K7WlycUGgl7I/HWxdNCl6IoqItrrV8gWfLh9nOqpCrxmNMRolwqlZWQH5g2FJKUa00sG73sJND2nNND7onvJ5Uo5EJoi2rRz90dMLvVfF5E2ahk2ET/6Df7MrYTc8fo5GkEKsJAwD7QQy+YyOFTxxIWPqMSzYFIr2sZEODpHrG0+HIaDSl7oli/H5rxOG13PbSnd7nUyCCnJP2zJONPMTSRcj9h1P8iVOMPiorMNo7hufgcYWTEGj1ixMBGBMpZqQcfpdJHciIEMFTZZz3vDCCEySFWE2fKxf+c/+4MbdSyKJQCJeWuh6YNKxSS4oRa5ZbXkMDA7ZdrkKm0DAypmP09aO3dVhzZ59vRGYGh1FIHsomm1vGdMx43q78/IRBRvvlK62QuxtPtDqA+PswjeG46qpqZeEUSsLTLeGfHH+XyVVIHp1NRyE5jfc8f4R5DIyOjgllRtSWLkn1XbaNQeKRRNSiwokbPMSVI9HoqENPo6sbufMAZnfv8GknrrSWVOw4ZcLjsVp6exkq/kHpLa1ZjULxaiJhThuLJow+4j2p2lCPr9PEPHAYNlsZJNW8QJLmXNj24gZmV7eVkH0onLAZGHkxS1mmFhRY9vpp4ZWFyz0yldAEAzdOlnkrzOnBx5NJLNinIfsHrOFW0sK94vdDWZGlVBkn28DwBURCwy3cbivi4yika1jNbbstn2e7N1erKiyNbMaT5/YUaK6SnhU0OeGBeego/j9vTtFQx1cKlEDAit5SkEfcj1m4reGzkpeX2d1USsyg5TttdHWPVKzqMeTh5lSjIu/YJsGjdirjaLTnz5w5jenMUdSCgkQK10TydUAtLJhcNIi0ufurgbizwXwmPv9VPB5wuUa6HyYH/hvjO0tRvI6SSih+vuL1pjrQTBAlELDcbvuteOObzMcWwJwZUs33pqpsUNQUZ3YzGEwMb0XAn9LLxnv4UdOoTnGde94iBEqaySSQsCMfcbid1WO0smYNW29hRmOWdV36blUdM4uGWlzMiDS+aZ5gcZRli4GRy2sTxQyFJrxENb+EORuMoWBIrCHbqGWlVquYKQxrNj7GOeooMSpSpoTKVfLzURuXWbbTaXG2EGI4F3E8YkYS6WaNuSSun9AqylCSgtUnEEpqKKM0jJ6eMb+jFCHv7hvRqGXKIKkWFY5oQNTycqtzmWDo5fklzDnW6OrNx2Y2e8Qcn+Ikk8kRxRwYwNh7wLqPJLNDxWtZ55m2a6Bwuwm/9riUDzp9eidcbkvgZ0AxNBp6e2eioRYeD8qJq1CLizNakiXXc1JeYQMDCSOjRBmahli5FPW4pHDK1RXWEmb8Gdh+7OnnjsX8EuZZZIQx/GQEcTQf23mEEvCNOZJIZHrQNITfZ60c2L2xjETwHe0fnsIIgUjLZClcmu3RZCbKmfGRS0oqWIn5yt5xA+mbG9ZMKu5apnjZZjQGew5B63DQCtnUatkg+LyJxOoTWaFJxhHmDKQMi2yb6mmF58nl3Hq0JZRpYvT2jd6AKWqi51YKC5ChIeRgKCV0rjx4NBHIXm1cRmxNXYpjSyKwXtzOeiIpgbKAWlxs2QoU5o878hOahranOfF81YICy87fjvI5YaSJGQ5bJrr2M4iPYsy+fpSkYA2TeZfzVps9ZWbAkiqRhHsOIFxu1JrK4RhbM3ahkdrbeDqWcZ+FHcsr7g01btlp70zx+62wxFn4doXHkxIQYPwTkuav8VRImUIRCYFaUpywNEsh2aQ0frgd9C/+LcVdI82BAYSmJXr4hWEBli1mYN49W6F7RmiQhUC4XehHp+WQNjGSbaNtzHA4VZBHmf9qNVUp7obp9xEfxsbdQ9NXE8wM6VqniowvF9lCmuLuaBsVsfHEYbsG29pM8fuHc5plqouUliCnY7vtqmmKN9VOoRR/fkp+PqIw34rbribVZ4ye+tUnzDPArCxPjaYoUpQR7p0zihxOXj6icRml4dSbmlOEPqUxVFTU2urEdqUgD3NxxfD+CcZhmwzaojp6r9o4vEGaVty24kLLwOel3SnLQ2pBAeaQnU86XReS/F6S1qoT1ddcVm7mtNjcenNLyv9GXz/6oSMYXd0THvU5wjxfSUrnohYUWD1HsgZ5JhRIcT/fJBSPx0r2ZptMThvTQD9wyLqcy40MhlDbelN7xiybufafWsv/fu3/OPb5MxI5vIyOYU13SmMjRCKMUYouZJRwQFbqnCQtfiya2f88veHL1BCaxsg4ZUk4wjwGcz3DouKz7cxVdUTeIysTQpZfr5QYfekJ9sKYkYjl251l/YuM2QnhXBokBSjI9kjI1xKm3/TyzKd/wMEflyBKii2hXt04oufV6hdbzhcVaQn70pfa7NGDjETGnYYNXLEReea6YSViIJA5w8U4OMI8BjO6xpwF4sMvo6dnRF2Trdyyyijzw2xeK914wjjajGE7qyRHY5kMyTG6Rux7Zis/ueYKru86hR1n/I6lf2ghtrgcY+fekQ1HNIZQxLjutuPqUZJGTb3LFbTeoYSNvhkMjp9rOgOOMC8AMlkUzWdGxtiynSJsN9gJNbJpOgXhdo+5Pqz+80Weu2INVx16HT+ueYZrb7uH0OWnjZiq6M3Hpj0yUPLzLTdJu46LHhlAasqks0qOKHdaZzvMCZK9gl7VJAlwelRMc2BgXAcRY+deei73sPrJ93OBv4cbfvAjWq87PfsGP6ZJdHk1Wm211Qhvfhlz686Mo5vJrF8vPGGeK/bO6fVQVI59/gy6P3D63KnjQiNJGDJZXk0EvbWNhg/t48RbPkWlavLoZ6/nwP+cmjmH9xgIj8eaX2eqZjCI8tQW9KNNGZV5QtOsZAiME18sjammp/maEKI5KQ3NhUn7vmynoNkthDgvaXtu0tNkQQkTf5DZqodwuen80Ab+/Invol/aM+Hg9WN6HS0wRvVMmyxj2HanW/aNhhkMsvTrL3D2Lz9PWEpevOqH7PrxCRnjW8dJ969XltSBMcYSYdzSzczwvapqQkcwGaaangbgh0lpaB4AEEKsAa4AjrPP+UU8WifzLT1Nlswihaax9/qT+dN/Xc8b/vZZqj/YMeE512SM7OcDY83rM3qmTYWxvJmSn2d6zqu0JTcZi1L3nU287f99ngM67Lr4F+z6yYpRtczpborG3gNWzzuF+qZEip1EVJIppacZg0uBO+342QeBfcCG+ZCeJpmEgfs0l3aEptF2zQaeeev3+a/mC1n9n3snbTy/kJjwvN52oZwwWWh00404ADANim99lms/92n+Gixlz5tu5OCNdSMEP0EGg5GU3WP4SI/JBI2ApvO1fsLOAnmzECLeXNUCyUbB8TQ0M5eeZqbIwnJL+LyT+OF/3MBfBhto/Y+lU1pueFViZk67Ohpy44lTF5RxC5cE7t7E9756JXtiUbad8Rt2/aQh8/XGc9SYhLcVMGyTPsMukL8EGoB1WClpvm9vHy0Nzcykp5kBFK93ch4wo6AtWcSp//087UY+v7vuYsTTW7NQO4dMKM/vnPFQRgV3buLDX/oMz0RUdp77f+z6/nGTXg6cdGOeKYTzGExJmKWUbVJKQ0ppAr8CNti7moDktHvxNDRZT08zU5jh8PQ+DEVFXd2I73dDfLj0KX7y5StwP/x89iroMIKc6BakJP/OZ/nPz32ErVHYcvGPOfSfJ+cumMIE7m+q6Wmqk/59CxDXdN8PXCGE8AghlmIpujZnOz3NXEYJ+On8Hnxr8X1c+d+fI3Dv5tmu0pxktAiqcx3/vZv4yPc/zYBp8Nj7rqf9Y6fNdpUSTGRp6g7gGWClEKJJCPFB4Lv2MtM24PXAZwCklNuBu4AdwEPAx6WU8bHCx4BfYynF9pOaBnZBIDSNff91PA+feCtv/u3nKL3p2QWljc4myvJ6K171PKTiF5t4w28+T5GiccvnfkjoLXNDoF99wQlmkP4rN3Lvd77H2U9fS8MH9s37kLQOoyM8Hnb/4gT2nP9/3NS3mPvecRbmK7usfXaggSmXnZQKJx0nOEEOUFc38v7/734+33QxjZ9pdwR5gSMjEVZ/bj9nbrmCDxcexfXzXtTSEhSv10pgN52ypxgyacELcy4cEJRAgM7vQak6SMdnF6O3tM74NR1mH6Onh/Jrw3yl/WTuWf43dn29EanrmVO05oAFL8wz7oAgBPv/60R+d/xv+N9vvRue3Taz13OYU+iHj7L5y6eyParz70u+T9d7Tp21uix4YZ5pouet5+53/ZDz//JZim9zNNevRtwPP897fv4Z/IrKJ7/4R9Q1K2alHo4wTwO1tIT6r+3i5q7XsOq/D056kd9hgSAltT9+ntOevob3FnQS+lGUvqs2TusbUPLzJ32+I8zTQFZX8PiW1ez82OqRkSdsBYZWVZmIreywcJGxKMv+K8QNvbX8/bh7Oeuzm1B8vqmXN4VQwo4wTwPzlV2s+Ohm2PzyqMcYtWVWDKtcRszMMWppidNYAcae/fzmW2+m3Qjx+fKnaH/P2omHOEp7fhldI8crYtJnzDPUNSusnjGbKCra0iUTOlS+sN1KtjbH1/Ong9k/mNUYYDO9AjGTqW+K7nyeMx//FBVqAC7uQqgTFLGJROccp8Fc8MIsDxzB7OnNapltnziNqju6suKQsRDIduDDsVYglPx8ouefitq4bOrlx9dxZyAHmNR1Vv13L3cNFvLgulvov2Rd9gofp8Fc8MI8bceJNOSZ6/j5dT/j8b2NVooUh4mRpZ5QDg3h396CzEbWjhnKAWbsPcC3f/Ju/ELl1C8+n53INRNgwQtzNlG8XuQ3uug1/az43tSSZ6cwQ0neFjJS19GPNk05xte0meD7qrplK1fuv4zrqzZx+MMrxz8hCzjCPAmaP34yd634A5+94/1WNMXpIrMfOH7OskDuc6IaajMYZOBbdTTpQ3zx6rvQliR5Bo/SIGRUlk2isXeEeYKoRYW87l3P8ZPu9TT8bP9sV8dhlkifso0V5M/96EtcuPljvDu/nV3X1Y4rmGppcWoCdphw8EdwhHnCHPnIcfx31ZPcc+vrxs1m4PDqQYYjoyvRTIO6n2jsikX49kV3DIfeHSW2nN7ahrFzX2r5k5jTO8I8AdTycj703gf4attZ1P36lfFPcJhfTGONXMaiYwqc8u9tvOXZjwCw+9qaUXtnrboKra52dI31BOroCPMEOPyhRq4q2M5Tvzx1REjV+Yo8fe2UkpMtRIRLmzmjF9Og4Ts6D3afyM/ecjNqQz1Ic+T82O1CjqHUU3zecQ1QHGEeB7WslA+++yHevutKym9fIEH5FJVQjdfxubZJiVM9A5hbd7K9q4o3+obYf3VVImdWMvrhoxnDMMeD6wuvZ1wDFEeYx6HjzSt4S/42+H7Fgvr4A/dscvJT5Qop8f+siHYjxH+87T7U8vLxz7GJjwSNru5xl+Ommp7mD0mpaQ4JIbbY2+uFEENJ+25IOic36WmyiPB4qLz6EJdv+RCefy4gP2W7F1K83pmLN/0qQy0upv3aMxi6bAOhy09j8O2npZilev+xjbfteC9X5R+i/dLlmQuZpkhMKT2NlPKd8dQ0wD3AvUm79yelrflo0vb5lZ4GME5dzf9bcj9Fv8hbmL3Y8npEbdVs12JuMknBktEo/g4DKQQ9jSq6T0FGh4fSMhJB3lLBYV3npA9vy6yvmOZa/LTS09i96zuAO8YqY76lpwFACPa9X+NHx96E9187Zrs2M4I4cgzz8ATyIeWaueCBNUnBMoNB8v64Cf+fNlF7/SZ8XTrK2tUIlxvF70dbsoi+dwywzOXiZ3WPk/cXhf53bURd3Zi1Kk/XyvwsoE1KuTdp21IhxEtAP/BfUsp/MYX0NFi9OF5mx5lBbVzGt19zDz/85hUUhZ6ZlTrMNHNWMz+DyqgUhEis+QqXhlBVhNuF8HpB05D5fqRbI1bsQ/ep6AEFUxNIBULlCkKCNiTpb4C8QxAuF4TLTCpXdBCKDPL6pVu57eXTqCrr48rFz7EzWMOluy/j3hV/4lf1fyX0vwbn/fgLVO/cO3Y9J8h0hfldpPbKLcBiKWWXEOIU4M9CiOOYQnoa4EawQu1Os46TR1HZ9Ykynh1soOTurSxcT+Q5jhDDPWQiFauCUATC50PJC4BiCaNZXACqZeseKfNheFUihQqxgCBSLFBiEK6QaMsHCIfceP1RKgsHKPcNArDI10OFu5sybYBaVw9eEaNGG8ArJIWKSpMO+2OlHIsV87b8Pdw72MiRSCl/O3wcF9bt47G7NhAtkPiPKRSeEKajJ5/bd62n9m4XwcpKbtEvpuBwlOaz3VylXcT+exsJtJrU/XMf2Wq6pizMQggNuBw4Jb5NShkBK9OblPIFIcR+YAXzKD0NgFZRxg8u+D1fvOs91C/QXjmXCE2zftxucLsQeQHMojwApEvF8KgMLvISqlLQQhKpwFCFIFwfgSEVbUDFrA0j+9ysWN3E3q2LuO5ND7LSc4wB00fQdHOm7xBhqfL7ntM4v3AbrXohD3efwNdrHqJOy+POgWLO8B3l2XAtYdPFOf5D1GlWHXqMEHt1Fye5FVwieYgf4N9hk7+El/Boxxp+WH8Pn/nre3lyw26237Uab7ek7EiYvbtLWSz3QX4A2dSC8eRKlglJuNxL3lP7yHe5kIYJFSUsfsRP5AYXNb1bMYPBrAkyTK9nfgOwS0qZGD4LIcqBbimlIYRYhqXoOiCl7BZCDAghNgKbsNLT/HQ6FZ9J+l5TzzLXfdTfNzjbVZl1EoLo84GmIfIDSK8bqSnESv0MlbkIVaiEy8DwSfR8E0+HSmRZmKKiIG7NoL6wm9V5rdS5uylSQ9S7OqlUo5iAVwhcCPIUDy6hYkgT1R76/i3k5edHz+HAk/VUlPVx39m/5c2vvBczX+eaon14hIuvdqygUB2i13Tz9ic/SlFxkCfaltPSXoSMKvxX9UP0GCGe7N/I/+46j4FBHxuXHuRM3yG2RcMcipVwOLqEC/O2E5Eag2aEAWlyd/+JHBgq5x9/PoX8I5Ku88N8hrdS+Sz0fHmIauN5UAQyEiFh/9VmP7NntoKi4lcEhq4Pjyo6OlCA7Dtd2tcdL6OFnZ7mdUCZXd2vSilvEkL8BnhWSpm8/PRW4Bt2fQ372L/Y+9ZjacZ9WKlpPiknkE4j1xkthKbRdNdKVpW3MfiGgRFabG3JImQ4AgV5mAePzIg/bNZRVGs+6PXA0lrC1XloQZ1QtQdvl44eUBkqUYkWCGJ5EC43EVVh3G6dE6uPsTKvjaWeDorUEIu0bvxCZ2uklvP8zQkhjDNohrm+6xS+Xr4dgNv6y3hb3jH8ykjrpU4jyDFD5US3l53REAHFxAVsjZby1551PPTkSRTtFHRt0FHzYqyobmfwx3WEC1XK3neYSyq3cuNPL0GNQOdGnconVCJFCmpE4u0x0b2CWJ6ge61BwW4NBFQ9O4gyGKVzQzFSAdMFnl6JkOAaNHEFdbTeCGpnH7KvHzMSQcZ0FJ8Xcyicu/n8KIyV0cJJT5OGOOU4fnjPr7j85s+x+BtPj9ivVVUidZ3Qqcvwb96P0TXRPPTTrJc9TBVuFyJgzRXN4jykqhIt9xEp0ggXCcJlgmiBRC/WUfNjFBeEWFLYTaknyOpAC28M7KRctbQAecJFDAMXKh6hoQqF9x85i49UPM5G70iN8sHYIB2mhw0eF38LeQmbbi4L9PLviMKA6WWtu5OAULil73hKtEHeV9DO1zvW8NXy4dWATiPIA8ElxKTGP3tWsb+vlI8tfYKvbbqExdXdhH5fTdnmTkRPP0ZnF0pxMSLgI9xQjrstaGngGxahdvaBqmC2d2KGQgiPB6Gq1nKQUFItrOJ5jhcAYwlzdmOmLAC6jy+gVJXUPJV5XVlvtcZSnge7JjffEQKhuUARlp1tfj64NKQiiFUVIiREi1wMVmkE6wRqGILLY2iBGGXFA6wqbqfe10Wdu5tGTytuDGq0IRSgSNHwCFfafG80hv1xnwzDa73De27sq+HxHSv5Ud0jieMMabI9FqXb8HNX9xvY21/O2eV7Wezu5OsPvo2XXv9vfv/cRpSgihIWSJdEDQnc/YJvLDEo2aJwx6LXodnGcyW7dPJfbidcX4rnaA8lQ2HujJzEyoEdoKp4godSnqvR0QEdoB06MqyIfGn7iKGqjERG16iOIchj5XUalxw2EkoggIzGYIwITY4wJxF862lc9NknUIBwmYu89AMU1VrC0DTrt8+HzA+gl+YxVO3FVAEhCFYphKolhkdi5hm4CyPUlfayurCNgBahzNXJqb6DFClDXLf3ndyy8mccM/z896E3c9vyP3AgVsByV39CQZMZlU5D4BUqPjG2Ab4hTXQMft23jLfn78KQkhjw4tDx/K2viBItyIDh5a83nUVVl8lli66gubOIWI8HNahStBv8nQZ5O7pwDwR5eqiSZ3z1rOjZwouBclYObgPDGHXKkR40Rwe0A4eyqvyZKuqiWmRXz4SX6dSy0mEb6hz29jISsRqdMXCEGdCW1dOzoYqPf/2PfPPlC1BXmCz59G5ae9dz9BwXKKCXxVjbcJQ1BS0Uu4IUqkPUuzqo1fopV00KFTeKbYMzXg/59Y41bPCE8SseDjeX0tfo4j1/uZaale0UKirNejGv81n9ULJCqN0I0mtCiQLdJjw9tJSv//MyRFRB+g2ICaqWdtHdH0CaglifB99RDW+3xPAIyreE+aP/fLwdYdBN1L4gsrMblEJkOEJl2J5W3CFYKg+PqHfKp9TbZ/3OQfge4XKjFOZndESYLrKrByYRAMCM33eOmYhu5lU/ZxYuN/6/F3JHwwN4hAtDmvyidyk/evEcvrnhz1yR3zPqufEe75mwh7/1reV/KjfjEa4Rx70QibJIixGVEhX4+KG3cGywkGuXPc5Pvv92BusEy37fxuFve4nFVGKDbpACNJO87Z6EkqZ4dxRtMIrpVnF1BaGzF6OzE6QcTiOavDabI9SCghk1QBEuN0pRoTXkzjKK14tSVore1Jz1smcCZ848BlKPsa2pFs9ySwhVobDG24Q0FIrUEPcH/QyYPk72HGWZy0XIjNFqQKkqMWyh+cBfP0nDH8P86JeD1Lm7ORwp49mepexpK8ft1nE9WET+UZ1wsYprSJJ3YICSQ838oegMylpepDSmY5gGi670gCmRemxMgVRgxBA10XLPQuOcUZAV1Qr8n4X6yFh0RgQZrOitomBhOJu86ntmgMhFpxK6tpe3Lt6CIRVuve8clt3dy753FxE4Kqh5qBXp99BzQiG+Dh3XQAzdb7WDiiFxbzmI0dePWlJk+aoOhTHDM+sjOx8wX7MOddOOrMfVzjbC5Z7zdYzjLE1NAOFyo/i8SF13BDFbZHnIP5eFTl2zAmI6xt4DM3odZ5g9AWQsijHBD0XJz8ccGJjhGi0A5nhHkU1Mnwtllu/XEeZJYGVzVDBzZCjikMpc7ZXByik222O5OR82yCyeO8qJ6LJK0FTEklG9Nx3mOfM5f5jTM08UIRC6iX5w5PqrwwJBCJTyUszD8zPW25zvmZWe4GxXwUJKyxvGYVyEx5P17Io5QUr0w0dnuxZTZs4Ls8M8xDCmlCzcYXrMw+bTYa4zL9xCFyBOz+zgsEBwhNnBYYHgCLODwwLBEWaHVxULOVneRNLTLBJC/FMIsVMIsV0I8Wl7e4kQ4lEhxF77d3HSOV+209DsFkKcl7R93qWocVhYmINZXOpU1HEzM+aSifTMOvAfUsrVwEbg40KINcCXgL9LKRuBv9v/Y++7AjgOKwXNL4RIeOvPuxQ1DguLbJqEaotqUGsqs1bedJlIepoWKeWL9t8DwE6sbBSXArfah93KcLqZS4E7pZQRKeVBYB+wYV6mqHFwGAP98NE5ZWQyqTmzEKIeOAkr9nWllLIFLIEHKuzDaoHkO4ynoplwihohxDVCiOeFEM/HWIAJ2xwcZoAJC7MQIg8r4+N1UsqxYsSMlopmwilqpJQ3SinXSynXu/BkOsQhizhpXRcGE7IAE0K4sAT591LKePrWNiFEtZSyxR5Ct9vbm4BFSafHU9FMKUXNAD2Dj8m7d0+knvOMMqBztisBQPYSd8yde8ouc+m+loy2Y1xhtjXONwE7pZQ/SNp1P3A18B37931J228XQvwAqMFSdG22U9ZMJUXNbinl+gkcN68QQjy/0O5rId4TzJ/7mkjPfCbwHuBlIcQWe9tXsIT4LiHEB4EjwNsBpJTbhRB3ATuwNOEfl1LG/bY/RmqKmgezcxsODg5zPgbYfGkVJ8tCvK+FeE8wf+5rPliA3TjbFZghFuJ9LcR7gnlyX3O+Z3ZwcJgY86FndnBwmACOMDs4LBDmrDALIc63HTX2CSG+NNv1mSxCiEO2U8kWIcTz9rZJO6fMNkKIm4UQ7UKIV5K2zXsnm1Hu62tCiGb7nW0RQlyYtG/u35eUcs79ACqwH1gGuIGtwJrZrtck7+EQUJa27bvAl+y/vwT8r/33GvsePcBS+97V2b4Hu26vBU4GXpnOfQCbgdOxLAEfBC6Yg/f1NeBzGY6dF/c1V3vmDcA+KeUBKWUUuBPLgWO+MynnlNxXbyRSyieB9Kj/897JZpT7Go15cV9zVZhHc9aYT0jgESHEC0KIa+xtk3VOmavMmJPNHOATQoht9jA8Pn2YF/c1V4V5wk4Zc5gzpZQnAxdg+YC/doxjF8L9QhacbGaZXwINwDqgBfi+vX1e3NdcFebRnDXmDVLKY/bvduBPWMPmNntoxgSdU+Yqk72PKTnZ5BopZZuU0pBSmsCvGJ7qzIv7mqvC/BzQKIRYKoRwY0UuuX+W6zRhhBABIUR+/G/gTcArDDunwEjnlCuEEB4hxFJs55Tc1npSTOo+7KH4gBBio63tfW/SOXOGeANl8xasdwbz5b5mU6M4jrbxQmAPlubwP2e7PpOs+zIs7edWYHu8/kApVoilvfbvkqRz/tO+193MsqY37V7uwBpyxrB6og9O5T6A9VjCsR/4Gbb14Ry7r98CLwPbsAS4ej7dl2PO6eCwQJirw2wHB4dJ4gizg8MCwRFmB4cFgiPMDg4LBEeYHRwWCI4wOzgsEBxhdnBYIPz/Gawa2WfQNVEAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"image = img_as_ubyte(skimage.data.eagle())\n",
"thr = threshold_otsu(image)\n",
"image = image > thr\n",
"plt.imshow(image)"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {},
"outputs": [],
"source": [
"def footprint_dilation(img, radius):\n",
" footprint = disk(radius)\n",
" return binary_dilation(img, footprint)\n",
"\n",
"def footprint_erosion(img, radius):\n",
" footprint = disk(radius)\n",
" return binary_erosion(img, footprint)\n",
"\n",
"def isotropic_dilation(img, radius):\n",
" dist = distance_transform_edt(~img)\n",
" return dist < radius\n",
"\n",
"def isotropic_erosion(img, radius):\n",
" dist = distance_transform_edt(img)\n",
" return dist > radius"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 13.3 s, sys: 8.58 s, total: 21.9 s\n",
"Wall time: 21.9 s\n"
]
}
],
"source": [
"%%time\n",
"\n",
"dilated_footprint = footprint_dilation(image, 75)"
]
},
{
"cell_type": "code",
"execution_count": 41,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 224 ms, sys: 2.33 s, total: 2.56 s\n",
"Wall time: 2.55 s\n"
]
}
],
"source": [
"%%time\n",
"\n",
"dilated_isotropic = isotropic_dilation(image, 75)"
]
},
{
"cell_type": "code",
"execution_count": 46,
"metadata": {},
"outputs": [],
"source": [
"def iou(a,b):\n",
" return (a & b).sum() / (a | b).sum()"
]
},
{
"cell_type": "code",
"execution_count": 42,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7fbc8b4d64a0>"
]
},
"execution_count": 42,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABBsAAAFpCAYAAADdkCCpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAACZYklEQVR4nOzdd3wUZf4H8M8zsy09gSSQCiF0FBAQxN7Fil3sZz2xoqind97PO9vZz3YWPMthx957V0CK0nsnIZSEkLbJlpnn90c2YQNpuzu7s+Xzfr14AVtmvmmfzHznmecRUkoQERERERERERlFMbsAIiIiIiIiIoovbDYQERERERERkaHYbCAiIiIiIiIiQ7HZQERERERERESGYrOBiIiIiIiIiAzFZgMRERERERERGSrizQYhxAQhxEohxBohxG2R3j8RUaJjDhMRmY9ZTETxTkgpI7czIVQAqwAcA6AMwFwA50opl0WsCCKiBMYcJiIyH7OYiBJBpEc2jAWwRkq5TkrpBvAmgIkRroGIKJExh4mIzMcsJqK4F+lmQwGAzX7/L/M9RkREkcEcJiIyH7OYiOKeJcL7E+08ttd9HEKIKwFcCQAq1NHJSA93XTHDVZyMfTMrzS4jYhZX58C+uaHNY0JVAasFEAK6VYHi0aFbFQivDqFLwOuF9GqAAITVCt1uhVQFpAIIXUIqAloSIC0SQpWQWvO3pWLRoSgSFqHDpnohALh0CyABq6IhWXVDAqj3OqBDINPiRA/VBZXzrCak+YtclVLKHLPrCAJzOETM4bY5vBfpl8Pw5bDDCqn6vcQvh/ekWjVYFb1btaVZm5CjNjKHE1QM5zAQRBYLm220tVduuOuiWCUkIAWgSEBv79vL/7XY+7vN/7GWbbX8G/D7P5r3IdH1fighuDeXdZjFkW42lAEo8vt/IYAte75ISjkNwDQASBc95DhxVGSqi3KWvsW49KsfcEZqktmlRMwLNb3x7pmHQdQ2oOLEIriOqcVlg2bh6NRlSBFe9FAkduoS2SqwU5NwShU/Owdh2pqDYbNo+PvATzDOXoU0xQYFCnToUKDAKtSud94lnnwlMjVvzUazawgSczgEiZ7D2yYUof6Yelw59FccnboMDqHt9XqPL4efXH44ku0e3DbwCxyWVIE0xdb6GuYwGSGGcxgIIovtxUUyf+qUiBRHUa7lHN+vYSAtEsIrWv9ufVzZ3S/Yu51FFLoNU27uMIsj3WyYC2CAEKIEQDmASQDOi3ANUUtYLFCysiBsVsgkO0SjC7KhAVpNLSAlth9ZgFNSqgEYcYAWGy7L2IqeH3wKTSo4KaUKdmH1PeNofU2W79OR4buwNcxWjqvGvOW3lRS/fyfO546oA8zhTvjncAv/HK48pAAnpVQBsHa8kTjTnRze0zBbOa4a/5rfIykdvpYoQTGLu0mqEkJrewW9ebSqSQVFg3aaBi0NBv9GA5DgnycyXUSbDVJKrxDiWgBfovms70Up5dJI1hBt1PR0uEf1x6Zj7Rh44AacmzcHPdV6pCmNqNOTsKypAK+u2x/yi56oO9hp0JWg2HJqSr3vX4lzcE8ULszhvSlpafCOKMXmY5LR95CNOL9gNnqq9a3PL2sqwIurxsP2dTpqD2n0O9lOHMxhImMxizsgmkfrC4ndJ9Ryj6v0Ok+giWJFpEc2QEr5GYDPIr3fiBECSlIShMUC0TMLekYKGvqkoj5fhSdVwJsEeFMldJuEbtNx8tg/cHvuE8izpO6xIQWACxOS1+GmHuvgGa1Bhw4e6BFRqBI5h+Uecww09gKOPXEubs59CsV75XCzCcnrcNMB6+AZxxwmIuPEfRZ3Qlpl85wtdq15/qxGFYqz+YLano0E//+zyUAUWyLebIgrQkAt7Yuq8b2xc19Ay3WjT0ElxvTchAJ7NfZP+gE5aiPyVRWpSsfDTYH2D3D9NY9oSLxRDUREnRIClr7FqDowD1UjBHOYiMgs7U062NFLPQLCI4BGpfWtQPffT0Sxgc2GFkLAUpCPhhH52DnYCq8D0O0SiltA8QKpZTq8DoGm7OY4lALQ9qvDE6PexGFJzg6G1SrgfapERN0kBCy9ctEwuhg7B1uh7xGrqWU6rM7mI1FntoKmngLafnV4ZNTbODqpjjlMRGSUNisT+FaO0UTryALdoUMkac0jEpp8k2axUUBEe0iMZoOiQqgqlNQUCLttr6c9pXlYfakFDxzyNo5L3ooMZe9Zxj1SgwIBVbS3vBaH1BIRdcqXw+0/J6CPHozVl1nxz4M+wMTUDzrM4ebbGAAL1D3ymDlMRNRdUpV7D9SSzSMO9n7x3o8rTQrgUnavckBE1I6obzYIIaA4HFCyMgH/A1UhoPXOgmZXUdfHAakKaFZAWoCG/OZlX6QAtOIm9M/bgT5pVdgnZR2KbFV77WOobRsGWluufLW/nFkiTsxIRASgeQ4ER9tbEERaGoTDDq13FnSrAt2uoq7Q3mbleHeaQEOBbJPD7bEKHVfmPIfhtpZ9dJbDzGIiSmzSJiF9Z/nCo+yex0AAepIGxalCKgAUCZnqBbwKlHp1r23AsscECFJAeNXmEQrdaSKw0UBEXYj6ZkPSYInhbzRhfOovSBau1sdVoWOAtQYOIZClOKD4HeG2P/qgMxxiS0TUEftgYJ/XXW0e2zd5DQqs1a05rADI2COLgUDyuLP5FIiIqIW06M13iAGAJgHdl7sSEE1+kyzqAmJX+6O+lEYFuzdCRBQeUd9sKLQ68UCvBR082/WEXkREFJo+1gY81PuPDp5lDhMRRVLLqg3t4WoNRBRN2NIkIiIiIiIiIkOx2UBEREREREREhmKzgYiIiIiIiIgMxWYDERERERERERmKzQYiIiIiIiIiMhSbDURERERERERkKDYbiIiIiIiIiMhQbDYQERERERERkaHYbCAiIiIiIiIiQ7HZQERERERERESGYrOBiIiIiIiIiAzFZgMRERERERERGYrNBiIiIiIiIiIyFJsNRERERERERGQoNhuIKCCVWgOcutvsMoiIEhZzmIiIYoHF7AKIKHbMd7lx3e1TYWnSUdPXgroBGnr1q8T+OZswJnU9Rjs2o49FIFVxmF0qEVFcYg4TEVGsYLOBiLrtvDmXo++MuYCuIQlAbwAQAqvtdqxJ2w9v5h4FZ590VA+wonaQF9nFu3Bo3hocmbEMg62VKLTYYRdWkz8KIqLYxRwmIqJYwWYDEXXLek898l+0AbrW9gkpoTc1AU1NwI4dsC9tPvhtOQBelpSEFVlHQuudhbrSNNSUKGgY5Eb/PttwWM5qjElej31sVeilJsEqVBM+MiKi2MAcJiKiWMJmAxF1y1VrJsH+w2LIQN4kJXSnE7rTCZRvQep8ILXlOUXFryk5mJVRCk9xNnYOTUZDoYCrrwvDS8owNmsDjkxdhv7WJmQoDh4AE1HCYw4TEVEsYbOBiLrk1N2ofqUIWa4y4zaqa9Dr6qDX1UGUlaPnTKCn76lGRcUvqb3wS69hcBdmorbYjhv+OgPnp1UZt38iohjCHCYioljDZgMRdenZXYOR/flaaF2/1Bi6Bq22FqithboayAJwx/gzcP4p0yJVARFRVGEOExFRrOHSl0TUKU3q+O9bE6Bt225qHSXv6ajw1ptaAxGRGZjDREQUi9hsIKJOfepMRd83t5pdBuy/LMPUzSebXQYRUcQxh4mIKBax2UBEnbrhx/OgrVlvdhnQnU4sfn8INKmbXQoRUUQxh4mIKBax2UBEHVrudmLAix5ABjT3edgUfbQNPzRxfXgiShzMYSIiilVsNhBRh85beCmUucvNLqOVtnodJs893+wyiIgihjlMRESxKuhmgxCiSAjxvRBiuRBiqRDiBt/j/xBClAshFvj+nOD3ntuFEGuEECuFEMcZ8QEQUXhUag1IeSUD0uM2u5TdpETOu0mo0RvNriRqMIuJ4hdzODYwh4mI2hfK0pdeAFOllL8LIdIAzBdCfO177t9Syof9XyyEGApgEoBhAPIBfCOEGCiljNgqTkTUfTeVHY+0zxYj2u7MzfhmFR7YMQ739VpkdinRgllMFKeYwzGDOUxE1I6gRzZIKSuklL/7/l0HYDmAgk7eMhHAm1JKl5RyPYA1AMYGu38iCh+P1LDo9X2gO51ml7IXrboa7356kNllRA1mMVF8Yg7HDuYwEVH7DJmzQQjRF8B+AH7zPXStEGKREOJFIUSW77ECAJv93laGDoJYCHGlEGKeEGLejio2eYki7dld/VDw3gazy+hQ308asNwdfQfgZjMyi5nDROZiDsemcB4Ta/UN4SqbiCgsQm42CCFSAbwLYIqUshbAMwBKAYwEUAHgkZaXtvP2dqdWllJOk1KOkVKOyemphloiEQVAkzqe+OQEeMu3mF1Kh8T8FZi6/kyzy4gqRmcxc5jIPMzh2BTuY2I1NcX4oomIwiikZoMQwormUH1NSvkeAEgpt0kpNSmlDuB57B4WVgagyO/thQCi97coUYL6ucmCAS9Xml1Gp6THjYp3+sLD21sBMIuJ4g1zOPYwh4mI9hbKahQCwAsAlkspH/V7PM/vZacBWOL790cAJgkh7EKIEgADAMwJdv9EFB6X/HQJtBVrzC6jS/kfbcJHDVldvzDOMYuJ4g9zOLYwh4mI2hfKahQHAbgQwGIhxALfY38FcK4QYiSah4NtAPBnAJBSLhVCzACwDM2z9l7DWXeJostaTz36v6ABst3RnFHFW1aOW2aehTOOfcHsUszGLCaKI8zhmMQcJiJqR9DNBinlL2j/nrPPOnnPvQDuDXafRBReFy+/CKm/LWv/xtFoIyX6vKmg4sh65FlSza7GNMxiovjCHI49zGEiovYZshoFEcW+as0J/eVcSI/b7FK6zfHzMty97WizyyAiMgRzmIiI4gmbDUQEAPj71iOQ+cVys8sIiN7QgB/fHwVN6maXQkQUMuYwERHFEzYbiAgu6cGvL4+GtqvG7FIC1uf9HZjtMrsKIqLQMIeJiCjesNlARHirLg/5H240u4ygaCvX4aqFF5hdBhFRSJjDREQUb9hsIEpwmtRx18dnwVtWbnYpwdE1pL2Zjnq9yexKiIiCwhwmIqJ4xGYDUYKb7QIG/K/a7DJCkvntajyza5jZZRARBYU5TERE8YjNBqIEd+Evl0NfujIs21bT06EOGQBl5FDoh+0HZfhgWPr1hZKcbOh+tMoqTPvsWEO3SUQUKRHJ4eGDW/9Y+vWF4nAYuh/mMBFR90hbTCxubAiL2QUQkXnWe+rRd7oApHGhpzgcaDxqX2w6R8PlI3/FCWlfI03xIE0RqNMlyrVUzKgahy9+PAD932qAnL8M0LWQ91v6dj2Wn+3EEJuxjQwionAKZw5vPF3isv1/wUm+HG5RrqXif9sPxk8/jUfp2/XMYSKiMJIKIPwX7PGK3f9u+Wec9h/YbCBKYBP/uAL5Pyw0Jt+EgOeY0RC3bseMQY8hV03xPWH3/QFyVaDUChxa8Btw3m9YcIYLp/80GYPvq4W2ck1o+1+wEtevPQdfD/k4tO0QEUWQ4Tl81CiI23Z0mMOAL4eLfwEu+AULzmIOExGFk9A7+X+cNhla8DYKogRVqTUg8+U0SK835G0Jqw3lt47Ho8/9B98O/cjvALdzI+12rDvmRZzy/iw0nTQWEKLrN3VAetzYOaMQHhn61TkiokgwOocrbhqP+6ZNCyqHj373d7hO3J85TEQUKoHdIxYCJOPs7DzOPhwi6q5L152BlK+WhLwdYbFgzb9G4etrHsRIu73rN7Tjqsxy3P3ENNROGhdSLb0/24x367ND2gYRUaQYncNf3vAgDnCoQW3jph7rcNeTzzOHiYhCJRH0iIU9R0HEOjYbiBKQU3djy6sl0J3O0DYkBMqnjMXPZz+MPEtqSJs61AH85a5XoR0+KuhteMvKcfuPZ4ZUBxFRJBiWwwAqrjUuh2/+5+uQB44IehvMYSJKVEaPSoiHUQ5x8CEQUaAeqhqJ3PdCn/lcP3gknr/6yZAPcFucmlKPHvdshJqZEdwGpES/t3RUeOsNqYeIKFwMy+FD9sOT1z9tWA6fkVqLzAfLmMNERAEyelRCPIxyYLOBKMFoUseMtw6HVrUzpO0Iqw2VNzcGPWS3Iy/3/Rxbzx0a9Ptts1fgji0TDKyIiMhYRufwocauYskcJiIiQ7DZQJRgPnamo8+MipC34zpyON4b+V8DKmorWbHhgEv/gJKWFtT79YYGzHlvODQZB+1gIopLzGEiIkoEbDYQJZipn14Abc36kLYhrDbsuroeJVZjhu3u6a7e36LhyCFBv7/o/a34oclqYEVERMZhDhMRUSJgs4EogcxxeTDg1dDvo9XGDcU7Ybia1iJbTcGWc9xBL8Gmrd2Iq+edb3BVREShYw4TEVGiYLOBKIGcP+sKyPlLQ9uIEFhzmYrSMF1Na/Ho/jNg6Vsc3Jt1DblvJaFGbzS2KCKiEDGHiYgoUbDZQJQg1nrqUfSiCsggF/71UQeW4plDXjWoqo4dn1yHzacXBP3+tO9X4dGqMQZWREQUGuYwERElEjYbiBLEmQsuh+2HhSFvZ/VlOZiQ7DKgos5ZhYrhZyyDkpIS1Pu16mq8+cmhBldFRBQ85jARESUSNhuIEkCl1oDkVzMgvd6QtmMpLMCUkz4xqKquPVb0KRqOHhb0+0veqcUid5OBFRERBYc5TEREiYbNBqIE8LeKo5Hx5fKQt7P+4j64MmND6AV1U7aagvKzPEFPUIbFKzF5xXnGFkVEFATmMBERJRo2G4jinFN3Y95/R0KrrQ1pO5aCfEw+71NYhWpQZd3zyLi3g56gTHq9cM/oBZf0GFwVEVH3MYeZw0REiYjNBqI491DVSPR6d2XI29l4YV9clbnOgIoCc2JyTUgTlOV+uRGv1hYZWBERUWCYw8xhIqJExGYDURzzSA0z3jwcWtXOkLajZmXh6LPmRPxqGtA8QdnBk36HkpYW1Pu95Vtwzw+nGFwVEVH3MIeZw0REiYrNBqI49n5DD/R9oyzk7Ww9ZzDu6zXTgIqCc1fvb9F48OCg39/vHQ1l3noDKyIi6h7mcDPmMBFR4mGzgSiO3f7ZufBu2BTSNtSsLBx8+TwkKzaDqgpctpqCjRODnJwMgO3Xpbhp00QDKyIi6h7mcDPmMBFR4mGzgShOzXF5MODV0K8ibTt7MO7v/asBFYXm9sM+gdorN6j36k1NWPnuIGhSN7gqIqKOMYd3Yw4TESUeNhuI4tQ5302GnL80pG0oyckouWC1qVfTWpyXtg47j+4X9PsLPq7A141JBlZERNQ55nBbzGEiosTCZgNRHFrlacCA/3oAKUPaTs3Jw/Fs3w8Nqio0qYoDnkk7ISyWoN6vrd2Aa+eca3BVRETtYw7vjTlMRJRY2GwgikMXLb0YytzQr6b1mLwR2WqKQVWFbto+r0IMGxDcm6VE3jt2VGtOY4siImoHc7gdzGEiooQSUrNBCLFBCLFYCLFACDHP91gPIcTXQojVvr+z/F5/uxBijRBipRDiuFCLJ6K9VWoNUF/OhvR6Q9pOwzH74OXSdwyqyhgjbRZsPj6r6xd2IO2b5bhvx0EGVhQdmMVE0YU53DHmcOvrmcNEFPeMGNlwhJRypJRyjO//twH4Vko5AMC3vv9DCDEUwCQAwwBMAPC0ECYsFk0U5/5WcTQyvlgW0jaExYKaS+qi6moaAKhCQd8J6yGswd27rNXW4uOPxxtcVdRgFhNFCeZwx5jDzGEiShzhuI1iIoD/+f79PwCn+j3+ppTSJaVcD2ANgLFh2D9RwvJIDb+9sh+02tqQtuM+fAQ+2O95g6oy1n193weGBzmEF0DJ+7sw3+U2sKKoxSwmMgFzuGvMYeYwESWGUJsNEsBXQoj5QogrfY/1klJWAIDv75Y1kgoAbPZ7b5nvsb0IIa4UQswTQszbUaWFWCJR4niyegAK3lkX0jaExYLKa50osaYaVJWxhtscWH9aetDvl0tW4colFxhYUVQwPIuZw0TBYQ53jTkc3DGxVt8QptKJiMIj1GbDQVLKUQCOB3CNEOLQTl4r2nms3SmapZTTpJRjpJRjcnpyVBlRd2hSx/MzJsBbsTWk7XgOG4H3o/RqWouDj14MJTk5qPdKrxe2N3qgXm8yuCpTGZ7FzGGiwDGHu4c5HNwxsZoaXbfUEBF1JaRmg5Ryi+/v7QDeR/MQsG1CiDwA8P293ffyMgBFfm8vBLAllP0T0W4fNGSi5I3QDnChqCj/sxulUXo1rcVf876ANjL4IbxZX67CEztHGFiRuZjFRNGBOdx9zGHmMBHFv6CbDUKIFCFEWsu/ARwLYAmAjwBc7HvZxQBaFof+CMAkIYRdCFECYACAOcHun4jauuXrc6GtDnHo7qghmL7/iwZVFD6l1lSsPTMp6PdrVTvx8mdHGliReZjFRNGDOdx9zGHmMBHFP0sI7+0F4H0hRMt2XpdSfiGEmAtghhDiMgCbAJwFAFLKpUKIGQCWAfACuEZKyRuBiQwwx+XBwJca2h+D2V1CYOVkB8barUaVFVbnHDETC9Izg56EreSDBiyf5MQQW3DDgKMIs5goCjCHA8ccZg4TUXwTUob0azHsxoxwyDlfFnX9QqIENuTXC1F8zjJAD/5YRRk+GHd88DoOcoRjkRrjVXjrcd4VU2D7cl5Q7xdWG5yfFOCnfd8P6v1q3pr5fsubxTXmMFHXmMOBYw4Hxl5cJPOnTjG7DCKiNjZMubnDLI6N32ZE1KG1nnrkP2sL6QAXAFZfmIUD7AYVFQF5llRsmCQB0d48W12THjca3+oNp54Qy68RURgxh5nDRES0NzYbiGLcZSsvgPWnxSFtQx3QDw9MfA2qiK1IuOfAD2DJzwv6/bmfrcNLtaUGVkREiYg5zBwmIqK9xdZvNCJqo1JrgHdaL0hPaFeFVl/RC6em7DKmqAg6LaUC2yb0Cfr93q3b8OiXJxpYERElGuYwc5iIiNrHZgNRDJuy+USkfRba1TRLXm9MOemTmLuaBgDJig3K6ZUQluDnuu33bhPWe+oNrIqIEglzmDlMRETti73fakQEAPBIDUtfHwrd6QxpOxv+1A9XZmwwpigTPDX0dYjB/YN+vzp3OW7YcKaBFRFRomAON2MOExFRe9hsIIpRM+pzkf/hppC2oebk4E/nfwmrUA2qKvJG21RsOqVH0O+XLhc2v9UPLukxsCoiSgTM4WbMYSIiag+bDUQx6h8fng3v5rKQtrHttP64JnO5QRWZQxUKjjhtPpS0tKC3kf/eWrxay6UdiSgwzOFmzGEiImoPmw1EMWh2k4b+r1aHtA01MwPjrvgDyYrNoKrM84/e36Hx0MFBv9+7dRvu+fEUAysionjHHG6LOUxERHtis4EoBp373Z+hL14Z0ja2nzUUD+b9YExBJstWU7DxTD3otd4BYMB0Fxa5mwysiojiGXO4LeYwERHtic0GohizytOA/i97ASmD3oaSlob9Ll+EVMVhYGXmunP8x7D07hX0+8XsJTjt58kGVkRE8Yo53D7mMBER+WOzgSjGnPnH5VBnLwlpG7UThuJf+V8ZVFF0ODu1DNtOKAl+A7qGQffV49ldBcYVRURx6fT5VzKH22FIDt9di4d2lhpXFBERBUUa0Clgs4EohtTrTUh9MwPS6w16G8JigedPO5GtphhYmfmSFRtSz6mAsAZ/77O2fDXevPkE/MRRvETUgRq9EZmvpTKH22FIDq9ai09vOZI5TERdMuJkmDom9NC3wS8RUQx5pGoUMr9eFdpGRgzCtKGvGlNQlPnPwDeA4QNC2ob983n4661/xrv16QZVRUTx5LGq0Uj/YXVoG2EOd8r+RXMOv1mXZVBVRBSPjDgZpvBis4EoRnikhrdfPxxa1c6QtrP2rHSMtNsNqiq6DLMlYd0ZITYJpETKO7/h2UtPx0mrjke9zstrRNSMOdw1I3P4pQtPYg4TEZkslBEkbDYQxYh7Koej+MXQrqap2T1xyQnfGVRRdDrt+FlQMzNC3o7yywLopzTgiP+7EaeuPg7rPfXQJFvoRImMOdw9RuUw5ixuk8PL3U7mMBFRhO01giSARYcshlZCRGGhSR1vfXAY+uyYGdJ2dk4YgOt7fAogfmY/39Nt2TNx8tE3IuWd30LellZbix4vzoLrdQeuGnk1to9ORW2pDj1JB3Br6MUSUcxgDndfuHL4xmFXYNsBGagt1X2Xy24OeftERBQYKQDhW4xJqp2vysRmAyUEj9SwxuPCNw1DMDFtKYotqWaXFJBPnano9+pWaCFsQ1gs0CbtjKtl1tqTpSaj8mwnUt4VIS1L509vagJmL0LubCDX99gmQ7ZMlDhacvjz+n1wUtpiDLTG1uSIzOHuC1sOz1+K3Pm7c3iDIVsmSlAtV6eN+RGlBOI/0kFonQ9z4G0UFLc0qWOTtx5XbD4Io/59HW6aeBk+P7QUh31xo9mlBeyGH8+DtmZ9SNsQg/vjP8NeM6ii6PbEqDehDuhndhlECU+TOlZ5GnDZpoNbc/ibQ/rguC+mmF1awJjDgWEOE0WZPc8JJdhooLDjyAaKSfV6E0b9ciU81Q4MGlSOHnYnMm2NqHSlYNn23vAsS0fqZqDnIifURWuQ3zATLU24wU/3xg9HKTg8KTbu+9zkrUfpa3rIV4c2nNkDo22qQVVFt6OSXLjxjF4o/Ndas0shilv1ehP2/epaoEnBoMHlSLftnsTPP4dzfq+HWLoW+U6/HH6mAT8czRyOZ8xhosiRVgnh2d1NkBYJ4Q3gxnqiMGGzgWLS+LmXovTSNdCdTkhFxU5FYKeqAloNCrSdbQ4I9zyU1Rcux+TpV2H+FY8hWQl+LfBIuWz1ubDOWrbXxxEINTMDp0ycCVUkxmAmVSg4+LQ/sOmpNOh1dWaXQxSXxs+9FIOvXdaaw7XK7gNb/xxu7+KZvnA5rn7xKsyZ/GhM3FLAHA4cc5gocvwbDUAHQ9s76JVKhUtIUvgkxm88iitPVvdB0R0adKez+QFdg/R6IV0uSK+36ytPUqLvo4uxz3dXRf2s1vV6E2r/V9h8r2oIao8ejNuyQ5vULNbcnfcNnIcNMbsMorj00M5SFN3m3juHfX+6k8N9/r0Aw7+7mjkcx5jDRGEm0P7KAAEMwmKjgcKJzQaKKZ86HXj/xmOgLV0Z0nb0ujoMubUcY3+fBJf0GFSd8W4sPwo9P1oW2kYUFTsnNSBLTTamqBiRq6Zg06n8DUpktE+dDnx+0xHQVq4JaTu608kcjnPMYaIw47wLFOXYbKCY8ejOfnj0yvNh/WqeIdvzbt2GXudvxeBPrsZ2rcGQbRrJIzXMfW0EtF01IW1HHdIfz4+ablBVseXmA7+E2iu36xcSUbfcUzk4LDk85INrUeGtN2SbRmIOh445TGQ8GaYzuHBtlxIXv6Uo6i11N6L0u0vwzen7wfLdfEO3rdXWYtC1C3Hy327G1eUHoF4PbZiskWbU56Lgw80hb2fNBT1wkCMxf9QvSV+L6iM5GzpRqFpyeObpQ8KSwwNvmI9T77iFORyHLk5fjdqDS8wugyiuiDCNZuAtFfHNjGZSYv7mo6jmkRrWe+rxZHUflHx8BaZceDX6X7QI2qrwzGgtPW5kTp+FDcfYcdTtN+Ld+vSw7CcQmtTxz/fPhndjaAe5anZPXHTi9wZVFXuSFRt2nOQCBGdkJgpESw7fVzmobQ6HuPRjR6TX25rDR982BU9W9wnLfgLBHDZGquLA1tPdzGGiEOx1khihWyekhfdoxBMzmklsNlBUqdac2OfFa3HVWZPx+aGlGPjnuVB+/gPQtbDvW9tVg8xXZuH+VRPCvq+u/OpS0H96VcjbqTxpIG7qsciAimLXPWM+gKUg3+wyiGLGdq0Bw6ddh6vOmoxfDsuPeA5nvDobL609IOz76gpz2DjMYaIg+Xp0Zo044PKZFCo2GyiqjPtlMkru+R2Ysxha1U5TamhoMn85zKt+vwDaitBGcgirDbZJ22Jiec9wOi11O7aeUGx2GUQxQZM6xv9wLfrcN685h6urTamjyW01Zb/+mMPGOS11O3YczRwmCoRUwMkfyXCRvpWCzQaKGo9V98WAv9VAulym1qEvNfc2imrNiR5vpoR8FdFzyL54fWhiTkjmzy6sqDvCySG8RN3weHV/DP57FaTHbWodcglzOJ7YhRVVR0XPXBxEsSBc8zIQRRKbDRQVtmsNePXR4+Fdt8HsUpA3023q6hS3bjkG6V8tD3k7685RUGxJNaCi2HfzyK9h4WzoRJ0q89bjjUeOg3fDJrNLQcGPTaauTsEcNt5No7+BpXcvs8sgin7C9ycKmg1cnSL+RLqJxW8hMp0mdYz76gZkTzd2hvNgOX5ZjtvLjzNl3x6p4be3RkCrrQ1pO5a83rjj0I8Nqir2XZC2AVVHczZ0oo5oUschn9+Inq/MNbsUAIB1zkrcscWc+XOYw+FxSfpaVB/GHCbqlihoNABcnSIuycg2kdhsINPdsGU8ht5ZYfqw3RZ6QwMWvLQvnHrk63m3PhtF75WFvJ3yM/vhgvTQl2uLF8mKDTtPbOStFEQduHjjkRh65yZIr9fsUgAwh+NRsmLD9olNzGGiTkg1SroMFNciOboh6GaDEGKQEGKB359aIcQUIcQ/hBDlfo+f4Pee24UQa4QQK4UQ5lw6pqjyQ6OClTcOgbes3OxS2uj1xlKcvuq0iO/39u/OCnkIs5KSgv3OXwy7MH+CtWhy16iPYMnPM7sMwzGLKVRfOO3YfktfeLduM7uUNnq9sRQTV54R8f0yh8Pn1v2+istb2pjDZAgBCE1EzaiGNtgjjCsygl/PoJsNUsqVUsqRUsqRAEYDcAJ43/f0v1uek1J+BgBCiKEAJgEYBmACgKeFEGpI1VNMq9ebcN0zV0H5ZYHZpexFq60F/tIjomu9L3c7MWB66JNjekcNxN35nxtQUXw5JWUbKo+M3NczUpjFFIpqzYm/PHUZxK8LzC5lL1ptLSxT02Iyh7Xh/fF/eczhPV2QtgGVx/QzuwzDMYcp3kXy5NRIUgEbJe2I5O0xRt1GcRSAtVLKjZ28ZiKAN6WULinlegBrAIw1aP8Ugw6adwkKn15gdhkdknMX452/HIc367Iisr/T5vwZytylIW9n+5gkFHJCsr0kKza4z6wGlLg+nmMWU0AOnXc58p6Jjvly2qMvWBaTObztgBSUWJnDe0pWbGg8Y1e830rBHKbARcmEkB2J1bkbhI6o/ryaKkIxbFSzYRKAN/z+f60QYpEQ4kUhRMsRQgEA/5sXy3yPUQJy6m6kvp4B3ek0u5ROOT6ZgxcvnYi/bBsJTYYvaZe7nSh82mrI/dJ1fWP0N0IE3DHkM1h65ZhdRjgxi6nbnLob6W+km77ccFdacvja8nHM4TjwlyFfxfuqFMxhClwMnRBzhQoKRMjfLkIIG4BTALzte+gZAKUARgKoAPBIy0vbeXu7P1pCiCuFEPOEEPN2VIW2xjVFpzkuBzLnbTW7jG4Rvy7AklMKMeCdq/GV02r4wa4mdZzw+RSoPy00ZHvSwZ+ZjhydvA11Y4vNLiMsjM5i5nD8i7UcXn9adthy2CM1nPjJjYblcMTXFoshJ6WUMYeDPCbW6s1blpuoRdhHORh4xV1amMUditCnxoje1PEAfpdSbgMAKeU2KaUmpdQBPI/dw8LKABT5va8QwJb2NiilnCalHCOlHJPTM66HPCesV3YcCG1zu1/+qOTdXIYBU37DY6echn2mXYuz1x2FHxqNae0eu/xUDP3nRkAP/YROWCywZkT3VUozZShJKDs2bofvGprFzOH4F3M5XFbemsP7PnstTlp1vHE5vOx0DP7HGsNyWPSIjtWVolGGkoTyw+P20mhYj4nV1JQwlk4UnKge6aDH7TGfISKx+okR3x7nwm+4mBDCf7r30wAs8f37IwCThBB2IUQJgAEA5hiwf4pBC3bkQ3o9ZpcRGCmhLVuF4n/ORO0Rdbh/0vk4e91RQW9OkzpOXX0c7NfYDJkFXnE4sOZf++O7A58OeVvx7MKDf4GaFZn7vyOMWUwBieUcLrp7JrzHVLXmcLAjHVpyOOlaK7TKqpDLa8nhHw99MuRtxbMLjvyZOcwcpjhh+EgHA89/Y3WuiUgRWvibMSE1G4QQyQCOAfCe38MPCiEWCyEWATgCwI0AIKVcCmAGgGUAvgBwjZSSY3MT1M7yTEDG7tAm6XFDzl2MhnOTMOTXC1GpBTa0cbvWgCE/XQLP+RZoK9eEXpCiYtONozD/3Ec5OWQXbuo5DzXHDDK7DEMxiykYcZPDkxwY9OOlQeXwwO8ug+c8lTkcYTf1nIfaowaaXYahmMNEwY9yiOrREfEuzP0GSyhvllI6AfTc47ELO3n9vQDuDWWfFB9S1of0rRc1vJvL0PfCSpx05lQUX7Uajxd/iLxODjLLvPW4fsNpqHimFP3e+wNegyZmqz9jDD666kFkKDzA7UqGkoQtx2gYOMPsSozDLKZgxE0Ol5VjwKVVQeXwAANz2HnqGLzz54eRoSQbsr14lqEkofxoiYHvmF2JcZjDRMGPJDB0BELLyXPs9tIjSorwTjMUH0caFFM8UkP6hvgZ16Q3NSHj1dmoey8Z5x16IzYfraJoRAX6pVUhxeJCg9eOtbXZKFuYh+KvPbD/sgzpztmGZaCa3RMH3DYXpVxmrdvGD12DnVYbpIf3VVNiiuccPv/gKdh8pBWFo7egX9ruWyPW1mZjy/w8FH3jhm3WcsNzeOzf5mKIjY2G7hq/72rmMFE8E74TWd+vGqns8W+J8DQE2GQISLhvNWGzgSKuWm9C8rb4O7jQnU7YvpiL0i8AKCrKVRVCVSE1N+zaZpTqG5pfZ/B+KyYNwvu9HwdgNXjL8WtA6nbMsWbwIJcSVjznsPWreej3FVpzuIVd24wS5nDUYA4TxTnZ9oq5/0kt51JIHLxDhiKuQZew7moyu4zw0jVIjxt6U1PzgZQBM5y3R3E4UHrOKtgFD3ADUWyrgkhOMrsMItPs0gHrTqfZZYSXL4db/oQrh4XdjpKzVzOHA6Tw8iMRkemkVYZ1VQo2GyjiPBCAzpamIQaV4J9FH5tdRcw5JHktdk4YYHYZRKZZ7MqDsr3a7DLighhQgv8r+sTsMmLO2RnzsGvicLPLIDKV7tChO3hMbBiudBkw4REQYVwilM0Gijge5Bqnep8MDLTazC4j5gy0puCW/3sdYr9hZpdCZAqPVCG9XrPLiAvVIzIxzMa7UgM1xJaMa//xNpSRQ80uhcg0ikuB4uLpmGE4YCrq8LubIq5Od0B6Ymxt9yjV1EOBVahdv5D2ckZqLdbfpkJYeJJAiadJWmN62ctowhwO3vlpVVh3u4U5TIkrXJMkEgUijN+DbDZQxP28ayD02nqzy4gLriyzK4htH4x7FvpYjm6gxMMcpmjBHCYiI/B2lOjEZgNFnEu3AJKBYARbrdkVxLYhtmSsmeQwuwyiiCt3ZkBq4ZkwkSgQQ2zJWHsWJ+ylxCQtHNZgFN6OYoAwTN3ArwpF3MJt+bxX2CCWRgmNjZuQnHLQfCjJyWaXQRRRa7fmhG11hkSjupnDoTru4AVQHGz8ElEI2LeJSmw2UMQ1NnJCQ6P0WNGEKr3R7DJi2llZc6D0zjW7DKKI4jGZcXqsaMJ2Lc6XEQ2z83vOhJLf2+wyiCJOeLl8AkWRMBwcsNlAEads4nBJo1iXbsJbtZzJOxT5qhNaVorZZRBFFHPYONZFGzC9Zj+zy4hpzGEiCpkAl76MQmw2UMQ5tjMJjKJVVuGxr483u4yY1ku1oTGfB7mUWJjDxtGqqzHt26PMLiOmMYeJKGRc2SMqsdlAEeWRGlK28t5WIw18cRd+aOSPcrBUIeBJ4eePEgdz2HjM4dBJfvqIiExn9KSljHaKqHrdhaRKTg5pJH3xSlw548/wSE72FgwFCjRrF1d5Ba8CU/xgDhtPX7wSV771Z7ikx+xSiIiIgiZVY5sNFkO3RtSFOqnDtrOJo5yMJCX6P7wSh448G78MfxuqYA8xEB6pwV7btlFjKSxA1eFF2DEasBU2INnhQpPbCleTDfalSej9mwv45h2TKiYKDXM4DKRE/weX4eDh52P2fm8yhw1gKSxA9UFF2HZAcw638LgtsC9NQv7PjcBPzGEi2kPL9SH+kgtKe0uISgUQvgGR0iIhNNH6+e1qVBqbDRRRWzU7lBoneA3eWFrVTvS4Og1nvHw83u3/OQ90A+CBhqphFvRZNQBQFKy8Igt3TngH56R9ALuw7v2GI4DtkxuQVxj5WomMsMKdBXVrNTi2wVjarhrkXpfJHA6CBxqqhlrQZ2kJYLe15vCZqe8jWWlnBasjgLKr6tGnKPK1Ehml5Qqy0Dh60lBsMoROoM3nUfjdebnnCiqii7sy2WygiFIhOSQ9TLzrNwKTemOfJy7GzAOeQ5aabHZJMSFDScKCa57E/MsBN1Qc2rrUezuNBp9clROZUezyQAW8bDWEQ0sOD3rkUsw95GnmcDdlKEn4/drH8dsVzbm7O4c7Xiq70JIa/sKIwqll9QRegQsbqUo2c0zGZgNF1FYtHaK23uwywk7t2QOeffqgaqgDnhQB3QqobiC1TEfmoiroq9dDhuFg31uxFX0v3IVjzp+KsVf9gX/lfY8MpftL3HmkBqtQ232uWnOiTurYqtmxxZsFTSr4uW4gajxJmFdRhIbyNNx4xBe4LmujUR9OxFiFigMcXb+OiGKHmpkBz779ULWPA97k3QebqWU6Mn/fDn39prDl8IBLg8vhlrl3Osrh7VoDmqTECncWnNKOlU15WOPMhVcqMZ/DdmH1azIQxb89rxDHrT2ukvsPyd/zOcN3HeZGQ8vHIm0S8Iour/LHDAO/Jmw2EBlIHViK1Zfm4pITvsMlmR8hV01uM5TWJT1Y6AZuXHkO9Om5yPxwMfSGhk62GDi9qQk9X5iFDe9m4IRjb0TFSW6cNHQxxqSuR29LDXLVemzXUrHVm4EVjfn4aP2+cK1JhzfHjT+NmoU7spe0O/z3x6Zc3PjteejzsUTKsu2QzkbImloAHhSlbsOOUzLQ57gdhn4sRGS8ZU0FkM5Gs8sIm5YcvuD4H3FZ1qfIayeHf3NZcfuq0yOTw8dMwZbjNJw8ciHGpK7f67WtObwqHVCA04+Zhbty5+51G5cmdbxdNxgP/3Q8St7RkLRmB+TOXdAbmyBUBUUpW1F5EnOYKFboDh3CIxLuyruZJ+R7zj0gLRLCK9o2fgSgJ2kQTSpEB8tptryv+QWiudHg/7owN1FiiZAyuj8TY0Y45JwveVNevPjU6cDTBx8G79ZtZpdiKCUtDZuu3Rf3XjodJyfXdute3WrNiUmrzob8v2yIXxcEt9/kZCjpaYDdBmg6pMcDfVcNpMvV9nUOB0RKMkRqCqTDDtHkgmxyQTjs8BT2xJapHrw3ehoGWju/PUCTOua7Nby+czw+/mkMbNUKcg6qwHlFc3Bm2ipkJ9DtBWremvlSyjFm1xEJzOH4cm35OKw5UEJ63GaXYiglORmbrx+Jey4LPIfxtx7A7EVB71dJT2v9v5Sy0xxuw2IJOodf2HEovv15BOxVCjIP3Yrziufi9NTlyEug2wsSKYcBwF5cJPOnTjG7DDKInuaFcKkQ7sRqNnTKgJN0qaDDJkGb13XUbAiU761S+O03wZoNG6bc3GEWc2QDRdRWTyakJ76WBlN75WLVowVYdNhjvomsujcpWJaajC+HfIJP/+fAvX/7E9Jm/AZ01PwTAordDlitEHYb6g7tj/KJHpy2zwIclr4AKYoLmlRQpyfhw6qR+GX+SJS854X9j3UQyUnwFOdA2nbXVX5wLvKOKMM5BfNwUNJaDLMlAei6UaAKBWPtCsbmzcNj58zb49nEaTQQUXRRe+Vi5cMF+P2IR3y3LASWw+9OT8cjfz+vezkMQCQloe6IgSg7Xsfpo+bjsPQFrS/zSAs+rBqJX2fth76feNrkMJTmo1JpEagYn4QeR1TgguLfgsvhwlnAubP2eDZxGg1EMc+rAPEy7N4Mfif0/nMziD0ivN15G0TzLRbNoxMQWnPA9742+02gRkNX2GygiFrYUAS93tjhqmZS09Ox6ZkcrBz3AlTR8URWnTkxuQkF9z+Oa+QNyPh2FWSTq/k+Yl0CIwZiw8npyNx/Ow7pvRZZlmpkWJw4L+3jdiYe0wB4cEbqT0Cfn1B2Sj0+rB+CImsV9rdvh9VvYs4sJcnvql/353QgIoo2LTm8atwLUEVweXZGai1KfTmc/uliAGiTw5uPTUfawc05rAiJXtbtuCTj0w4ngGzN4TPb5rDqy2EFzGGiRKc0csWavXQ0t0PLY36NA+nrE0DuMTfDHif67d6mwsZAxLDZQBG1qaEHoFWZXYYxFBUr7hmC5eOegtreEokBGGm345kHH8Psxn6o9KRhfWM2tjWl4Z4+/8VwW3szZnU9w3mhJRXXZG72/Y9Xu4ioWZkzE5Dxk8Or/j4MS8c9YVgOf3LHCOgQ2NjYkzlMRNRNhq/8oEhA32OZRb/tx81kjHGOzQaKqDVV2Sjwxsd8Da7jRuGLiY/CLoy5fWC4zYHhti2+/630/c2puYnIWPGWw5+d/bCxOZy9co9HmcNERF0xeqLLhFmtI85x/A5FVMP2+LivX9jtcF67q8uJvIiIog1zmIiIiCKBzQaiQCkqNt42Gp8Nf9nsSoiIAmbdpZpdQuiYw0QU4/Qk3gdA8Y/NBoooa4ar6xdFucrLx+KTSx9MqGUeiYiiyc6LmcNEFOPsmtkVhEzyTJK6wG8RiiirzWt2CSFRhg/GX25+HaVWTvRFRGQGZfhgTL3tTeYwEcU0ocb+MgicpJG6wmYDRZSI5blehMCKG1JxdmqN2ZUQESUmXw5PSqs2uxIiIiLqApsNFFE2S+yObPAeOQrvH/kfs8sgIgqJNy92b2fTDx2Jt458xuwyiIhCJp2xuyigVADE8gVEihg2GyiiJhQth5ISe/fYqunp6PnPDRhpt5tdChFRSCYNnxezOZx17yaMtVvNLoWIKGRKY+yehgkdQOzfBUIRELvf5RSTDkpbBSUt9u6z3X7WMLzY91OzyyAiChlzmIiIiCKhy2aDEOJFIcR2IcQSv8d6CCG+FkKs9v2d5ffc7UKINUKIlUKI4/weHy2EWOx77gkhYvrufQpSpuIEkpPMLiMgwm5H/kXrkao4zC6FEhizmBKZsNvR+8INzGEyFXOYjCRVCRkHk0QSdaY7IxteBjBhj8duA/CtlHIAgG99/4cQYiiASQCG+d7ztBCiZUHvZwBcCWCA78+e26QE0MfihNYzzewyArNPfzzc912zqyB6GcxiMkCKcAN2m9llBEQMKsF9fd83uwyil8EcJoNIi4S0xVizQYBzNVBAumw2SCl/ArBzj4cnAvif79//A3Cq3+NvSildUsr1ANYAGCuEyAOQLqWcJaWUAKb7vYcoqlXul44SC6+mkbmYxWSU/tZaaLmZZpcRkB1jszDEyrkayFzMYTKS4lJib94GCc7VQAEJ9ju8l5SyAgB8f+f6Hi8AsNnvdWW+xwp8/97z8XYJIa4UQswTQszbUaUFWSJFI4cQ8KbE1gFjXR/A2noxgiiqhC2LmcMUTZjDFMUidkys1TcYWjiZTDSPbiCKZ0a309obWCM7ebxdUsppUsoxUsoxOT15cBFP0hQbmrJjq9mgJfEXAcWckLOYOUxRhcN2KfYYfkyspsbeKjLUMalISIdudhkBYXOEAhVss2GbbxgYfH9v9z1eBqDI73WFALb4Hi9s53FKMDs1F1LKGs0uIyCKl0e5FLWYxURE5mIOU1CEJqDUx1gzX+cxMQUm2GbDRwAu9v37YgAf+j0+SQhhF0KUoHnSmzm+YWV1QogDfDPuXuT3HkogTgkoTV6zywiI2sRgpajFLKaAuSQgPLF1awybvhTFmMOUMERsDcSgKGDp6gVCiDcAHA4gWwhRBuBOAPcDmCGEuAzAJgBnAYCUcqkQYgaAZQC8AK6RUrYc0UxG8yy+SQA+9/2hBLPCkw1l207EVFZxxBhFAWYxGYU5TBQc5jARUWC6bDZIKc/t4KmjOnj9vQDubefxeQD2Cag6ijuz6gdAr60zu4xuU3vl4sIzvjW7DCJmMRmGOUwUHOYwEVFgYmy9FYp1H6wbDt3pNLuMbqs8vhRTeiw2uwwiIsPMWLEf9IbYmdW+cgJzmIjij1Q5ZIviH5sNFFEje5dDSUoyu4xu2zWhAcmKzewyiIgMM6ZoMxSHw+wyum3nsU3MYSKKO0KLrbloJM8aKQj8tqGI+mfBJ9D3KTW7jG5RkpNx2qBFZpdBRGSovxd8CjksdnL4jKF/mF0GEZHxBGJqWV9ODknBYLOBIqqvJRmVI1PNLqNbRHIShidvNrsMIiJDDbQ6sGN0utlldAtzmIjilgQnv6W4x2YDRZQqFDjzY6ONu3PCAByTvMnsMoiIDBVLObzrqAE4LGmj2WUQESU03kJBweK3DkWcZov+Nq6lIB+n3Po9ctUUs0shIkpIloJ8nPS371FoiY3RcEREgYqVSSJ5CwUFi80GijjFE/1X1DZc1Be39FxmdhlERAmLOUxE8S7WJokkChSbDRRx1nqzK+icmp6OCWfOhlWoZpdCRBQWlihf+ZI5TEREFPvYbKCIqtEb0XOJx+wyOuUaMwB/yf3Z7DKIiMKiWnMie3F057BnRCluzPnJ7DKIiBJbjK2YQdGHzQaKqCVuO5JXV5pdRqe2j7JzrgYiilsL3alIXr7V7DI6tXVcEudqIKL4Fgsn8lwxg0LEZgNF1H42Lzack2d2GZ2qL/WaXQIRUdiMsTux8dwis8voFHOYiOKeBGS0NxuIQsRmA0VUsmLDY5c8D2X4YLNLaZeamYFzD5htdhlERGGTqjjw6OVRnMPp6Tj7gDlml0FEFHZc5YHiHZsNFHHHJnuw4qp0s8tol3u/UkzJnmV2GUREYRXNOezZrxTX9fzF7DKIiBJbLNzmQVGPzQYyxT1HvQtLYYHZZexlx0gH52sgooTwzyPfg6Ug3+wy9rJ9FOdrICIyHedrIAOw2UCmODN1K7ZNKDa7jL005jBViSgxnJ5ahu3H9jG7jL0wh4ko4XAEAcUpNhvIFHZhRdo5WyCsNrNLacPTQzO7BCKiiEhVHEg5twLCYjG7FCKixBaFPVY9iRNKUOjYbCDTPND/HSj9o+uqmlrPHwkiShx39f8ASkn0jTIjIkooUTg/gnBHWUEUk3hmRabZxypRs08Ps8toQ21isBJR4hhjc6NmZI7ZZbTB2dmJKOFE4fwIQuMxMYWOzQYyTbJigzNXNbuMNpK2C2iSR7pElBiYw0RERBQubDaQqbxRtPCDkpKC/metgir4Y0FEZAbmMBElKmmNsqENRAbgb3MylcVpdgW77ThvOJ7v+5HZZRARJSzmMBElKuHhbQsUf9hsINNoUofaGB1dXM+xY/Dwbc8hS002uxQioojRpA61KTpyWDt8FB64bRpzmIjIRJJnh2QgrndFpmmUbqSVecwuA2pWFnr9cxUO5xI/RJRgoiqH/7UGRyVx+WEiIjNxkl4yEntXZC5h/pCxrZMG49mib8wug4go4jRIiCgY2MAcJqKEZ/4hMZHh2Gwg0yhQoDnMTVYlORmDL1iBZMVmah1ERImKOUxEBEg1Cjq/RAbjbRRkGg80WJzhHaslLBaovXLhKc6B7lDhyrTCscMFS50LYuMWeIeV4P6ipwCkhrUOIqKEpaiw5PWCpzgHUt3dYG7JYW1gMe4ufAbMYSJKZMIbwQtwAtAdOpTG3dedpYLmkW7seZCB2Gwg0ziEBU1ZFljDsG01JwfbTu+P9DO34Ma+X2N/+3ZYhUCysMIpPdihCXxSvy8W10kUW3iAS0RktJYctp26HX/p/yXG2bfC6nfr3A5NYEbNGKysV1BqZQ4TUQITgG7XoTSFb9C5VPzmY5CAcPntS7DRQOHBZgOZxiM12Hd5Dd2msNtRedEoTLjmF/w95zPYRUsrY/eBbDJsyFaBIT3WAj3WGrp/IqJEJ6w2VP5pdKc5DADZKnBnzjIgZ1nkiyQiiiYSYW00AHtP/Njm/2wyUJhEfbPBC06JGq+cUoPFadzM42pmBpY/PABzj3sY2WoKEJYxE0SJhzkcv1xSh9pocA7fPwi/nfgwcpnDRESmaTOSgcgkUT9B5MrtvVDmrTe7DAqDBl0C0phWqpqejk0vFmLF8c/4Gg1EZJSVW3thE3M4LoUlh0/+j6/RQEREZmGjgaJB1DcbdAuQqUT9AAwKQok1FVtu9EBYQvv6CosFyx8cjPnjXvYbrktERtFtQA/mcFwyModX3DOEOUxEZDKuakHRpMtmgxDiRSHEdiHEEr/HHhJCrBBCLBJCvC+EyPQ93lcI0SiEWOD786zfe0YLIRYLIdYIIZ4QQnRrytWhPbYhVXEE8aFRLPhozHPwHjw8pG3UnjEGP5/wKA9wKa6ZmcXM4fhmVA7/dOojzGGKa2YfExN1h9D47UTRozsjG14GMGGPx74GsI+UcjiAVQBu93turZRypO/PVX6PPwPgSgADfH/23Ga71OgffEEhKLWmYu2k4K+oKWlpKL5+FQq5ogTFv5dhUhYzh+NbqTUVa89hDhN1w8sw8ZiYEgD7BBRnujyClFL+BGDnHo99JaVsWUZgNoDCzrYhhMgDkC6lnCWllACmAzg1qIop7vxp/C9Q09ODem/90UPxTPEnBldEFH2YxRRO5x4wG0paWlDvZQ5TomAOExEFxojLVZcC+Nzv/yVCiD+EED8KIQ7xPVYAoMzvNWW+x4hwWsbvQGHvwN8oBCrOciFLTTa+KKLYwyymoE3KmgNRlBf4G4VA2Sle5jBRM+YwhSaE6RY4VwNFo5BmhBJC/A2AF8BrvocqABRLKauEEKMBfCCEGIb2BwV1+BMhhLgSzcPLUFzAScniXY7ihZYe+P3gam4O7hj1WRgqIoot4chi5nBiCTqHe/bATWO/CUNFRLElEsfEalaWsUVT/BAAJO/BoOgT9MgGIcTFAE4CcL5vGBiklC4pZZXv3/MBrAUwEM1dW/9hZYUAtnS0bSnlNCnlGCnlmJyearAlUoxIVlS4M+0Bv88zuAAnpKwPQ0VEsSNcWcwcpu7wDCvGOenLzC6DyFSROiZWU7mkbNwTgG4PYs1KyaUuKToF1WwQQkwA8BcAp0gpnX6P5wghVN+/+6F50pt1UsoKAHVCiAN8M+5eBODDkKunuODUNVhr3AG/r6avAz2VpDBURBQbmMVkNuYwJTrmMBlKAoqLkzJT/OhybKwQ4g0AhwPIFkKUAbgTzTPt2gF87VutZ7Zvlt1DAdwlhPAC0ABcJaVsmUhnMppn8U1C8/1s/ve0UQJzSkBp8gZ8m5qXtwhTAmEWUzg5JSBcWuC3C3PULiUQ5jBFK2mREF4GMkWfLpsNUspz23n4hQ5e+y6Adzt4bh6AfQKqjhLCt86BUDZvhRbg+7xJAqpg95cSA7OYwunrhkFB5TBRImEOU7Ris4GiFc/UyHSPLzsCWtXOrl+4B8dOCU3yBjUiolA9uuiooHKYiIiMI4M8M1OaeEpH0YnfmWSqas2JzLfSABn4cj3JO7yol64wVEVElDiqNSdy3ktmDhMRmYyTPFK8YbOBTDV500lI/3RxUO9NWbYd81ypBldERJRYQsrhJVvxY2NPgysiIqKA8S4KikJsNpBp1nrqseWh/tAbGoJ6v7a5HHesOtXYooiIEshSd2NIOezdvAX/WnO8wVURESUYIxoFgQ9OIwo7NhvIFPV6E0545RYkfTQ/6G1IrxfaW7lw6oEvm0lElOhq9Eac+fJUJH04N/iN6BpzmIiIiNrFZgOZQoOEbVgNnKeOgeJwBL2d7A9X4LatBxlYGRFRYnBJHcq+NXCeNpY5TERkFt7+QHGMzQYyRYaShMXjXse7jz8K76c50I4YBYjA01arrsbPL+8Pl/SEoUoioviVq6Zg6fjXWnNYP2y/oLbDHCYiCoEEb4GguMVmA5kqV03B10M+xj0vPI+tU8ZDWG2BbUBRsWuYFxao4SmQiCjOteTwbS9Ox/arD2QOExGZSQBSZfeB4gObDRQVDnCo+Pamh7D23tFQUlK6/8YxQ/HWhP9AFfxWJiIKxVFJGj677UHmMBGRmSQgNN5bQfGBRwYUNbLVFCw873GsvWN4t6+srfpTEsbarWGujIgoMeRZUnfnsMXSrfcwh4mIiKg9bDZQVElWbPjlgodRedHoLl+rZvfEXw7/NAJVEREljmTFhm/Ofwg7z9+/y9eqWVm4+tBvI1AVEVECEujeBJICkDyroyjEb0uKOrlqCi646XNYCgs6fZ1ncBEmpq6MUFVERImj2JKK827tRg4P64ML0hdGqCoiogTT3ckjJSD0cBdDFDg2GygqXZO5Fhsu6tPpa2r7JSFXTY5QRUREiYU5TERERKFgs4GiklWoOOaMOVDT0zt8jTcpggURESUYq1Bx1OlzoaSldfgaDtslIjIfV6+gaMXDBIpaN+f8AM/wfh0+r9kEZz8nIgqjG3O+h7ZvxzlMREQR0NXcDVx5mKIUz9QoauWpyagc0fHwXFudhCZ5gxoRUbgUW5KxY78AlsEkIiLjdTF3g3BzqUyKTmw2UNRShYKGoo6TNaXCg1q9KYIVERElFuYwERERBYvNBtrLIncTavRGs8sAAJx7/E9Qc3LafS552Vb83JQd4YqIiMJvvsuNCm+92WUA6CKHl2zBl87OV6wgIopZsbKkZHeXyCSKsFj48QmaR2pwSY/ZZcSUTd56TPzyOugyOiaa+Wv2Amw9s3+7z3nLK3DrH2dEuCIiCgRzOHDrPfU469Pr4IiSOWk6zeEtFbhzwckRroiIKDKkiJElJbu7RCZRhEXHkUyYWIUKu7CaXUZM+bh+CI4ftRhZUbKUmV1YkXvWJgirbe8ndQ29piehLEqu/hHR3pjDgfusYQgmjFsYVTmccfoWCItl7yelRO9XHMxhIopL0dZoiIlRFkR++C1LbVyTuRlPF8w2u4w2bunzBdSC3u0+l/T1QkxceGmEKyIiCp9ozOHbSz6DWpjf7nPMYSKiyIi25gdRV6K+2SA5JijhjbLVoalf+/cLS5cLPR5IxnoPr6oRhQtzmPa31zCHiYiIKCBR32xYvjUXRy07Be/Wp2OTtx5lvj9rPfV4uTYX02rav9JC8cMqFGiOjr9VLbua0MRxZURhs3xrLsYvPKM1h1v+rPXU44Wa3rivcpDZJVKYWYUCzd5xzlp31KNOtnObBREZqrPDHWlhYzjR8XCYok3UHxlYtjfAOmELXsw5EM/n94S0qoAuobg1QEqM/t8SAFvMLpPCSIECqXY8xa7weOGJ/r4ZUcyybG9A+okb8GKv5hyGaP55FF4dkBJDX1xpcoVERImhzTB6gbaTAupcjiARSKXj2yl4mwVFm6hvNgCA9HrhrdgKVGxtfUwHsP7+8fgkZyFiYIAGhUCHDsXTcbdeuNzYqSWj+buCiMJC1/bKYQlfDveeB+Zw/BOdXDRlDhOZYI+fSZ5oJgZ+nSmWxOzRoaVPEe48bQbUKFkajMJno1fCsaWuw+f1rdvxXvXoCFZERABzOJF0J4dn7BwbwYqIiKhdHOBCUSRmjxC3nFyEs1O3m10GRcAbu8YC68s7fF5vasIns0dFsCIiAoCK4wuZwwni+apDgLWbO3xeb2rCl7+OjFxBRNSK9+knju58rSWbDRRFYjOehID7iBpYhWp2JRQBb359MPS6jq+oAUDpmy6s8jREqCIighBwHlnPHE4Qn3y3P/SGzjO23ztNzGEiE3R2ixPFl+7cQsHbLCiaxGSzQUlKwsR+i80ugyLg1yYdA16r6fJ16uwlOH3+lRGoiIiA5hw+fcBCs8ugCGAOE0U5NhuIKErFZLNBWCzIs3V94BMsj9TQ7+tLMWTa1figITVs+6HOaVLHhZ9Nhr5weZevlV4vCu5X8JXTGoHKiCjcOeySHpR8egVz2GSB5nDhvxR81JAcgcqIKFKkym5GTOLtFBQFYrLZAFWFXfGEbfMeqaHoHQv6PrIYq129wrYf6tx1Ww7EkPvLANnNX3JzFuOvD1yOtZ768BZGRGHnkRr6vgvmsMku2XQ4hty7sds5LOcuxt3/upg5TBQpAuE/qeQkAKbjvBwUq2LzW1fT4NRtYdt8y1KLSkoyiqw7w7Yf6tizuwqw7spSeMs6nhiyPT1fmIOJT9+Kpe7GMFVGRC00XjaJa8/uKsD2q4ualzwNQI+XmcNEESMR9tsohN48ukFaOcLBLEHNw8AvF0WBLpsNQogXhRDbhRBL/B77hxCiXAixwPfnBL/nbhdCrBFCrBRCHOf3+GghxGLfc08IIYI+StUbmzBnV0mwb+/SFk2DfZsT0utFk+Sw/EibWjEK7196FOQfSwN/s66h4MFZuPT/bsK0mnzjiyMySbRlcaRyWKgKHMIbtv3Q3jSp49rycSHn8OW334hndxUYXyCRSaIthyNJ6ALQor7M+GPAp5yjIshM3fn2exnAhHYe/7eUcqTvz2cAIIQYCmASgGG+9zwtROtU5c8AuBLAAN+f9rbZLdLjxqzFA4J9e5c+rBsOZX0ZXCP64ujkdWHbD7W13lOP/t9fguXn9AFmLwp+Q1Iic/osfDBxPPq9fRXerU+HR2rGFUpkjpcRRVksPW7M/n1gMG/tltd2jYWydjM8I/phYloQJ7wUlPWeegz64TKsP6tXyDmc/sZsfHTqAd3KYU1y+nSKCS8jinK4RThPJlu2LQVXOYhV/LqRmSxdvUBK+ZMQom83tzcRwJtSSheA9UKINQDGCiE2AEiXUs4CACHEdACnAvg8mKIBoPhToPLEBmSrKcFuokNP/3IUBu6agy2H2FFo4cRk4eTU3filKQXXzjsXhS/a0P/7RdA8bkO2ra1aiwFT1uHFx47APScUoHZ8I8aVbECatQl2xYtdniTM2dwHlt/T0Nhbx5VHfYsz0/9AocUOxa8Px6X9KBpEYxYXfSlReWp4cnj6rIMwsHYOyg5LQjFzOKycuhvfNGbi5vlntuawNww5/K+j81F9WBPGlWyAXW0eraJJgTmb+8A2Ow0NxTrOPfxXXJI1C4UWe+s2FCjMYYoK0ZjDQHhPJlu2zRPWCBFoe/tDGG6FkKqEaBmlIto2klqekwq/5mSMLpsNnbhWCHERgHkApkopqwEUAJjt95oy32Me37/3fLxdQogr0dzxhQPtz2qd/MNyTNl8Il7t+0MIH0L7Tt9/Hn65cDwc+3G+hnCp8NbjwG9vQMEnFmT8sQ0lG5YCumZ8pkoJ7/qNyP3PRuT+B9hpsWCnqkIIFdLrRB/v7iVUf8zMx/f9x6G+bwo0q4DQgbpiBd9f/1BYTqaIDBKWLO5WDv8YgRwezRwOl5YcLvzYgvS55SgpD28O93x+I3o+35zD/vxz+Pf0LMwbeBXq+6a0zklXV6ziyxseRB6bThS9InJMrGZlGVw2JSLhfzuMbHunRstz/o0GNh4oFME2G54BcDea+213A3gEwKVo/84i2cnj7ZJSTgMwDQDSRY92X6fX1WH1M+NRed+nhp8IPpL3O6rv+wV1UgfAgxujzW7ScMWzt2LQo3MgvV5E8m5s6fUCXm+733zarhpgXg1S5u1+LLNnD/zngv1xZ86yiNVIFICwZTFzOL791ARc9+StGPSkSTncAa22Fpi3pE0OZ2Rm4NmLxuGfObydhqJSxI6J7cVFnPIvDrU5mTfjK9zFPtlooFAEdZeXlHKblFKTUuoAngcw1vdUGYAiv5cWAtjie7ywncdDkjXjdxzw8zVhudczS03m0N1uqtacqNacXb5uvaceo+adg79fcjnyH5rV6QFntNCqduL1Tw7j/cQUlaIhi7Nm/I79v7meOWyyQHP47ksuQe/HYiSHd9XgzU8OZQ5TVIqGHG6tRWUvIhbxZJ7iWVAjG4QQeVLKCt9/TwPQMivvRwBeF0I8CiAfzZPezJFSakKIOiHEAQB+A3ARgCdDKx2QLhcG3rgFJ7x+Cj4b/BFUwelWI82pu3HA/6Yid76O8tM8OGOfP3BM+hI4FA90qaBWd+D9qtH4ce5QlL7jRu7MpZAG3QscKf1fqMCbZ+bg/LQqs0shaiMasli6XBhy22acUMwcNku93oQD/3szshdpKDtBx5mj5+GY9CVQfUewDbod71eNxk8zh6Hfhy7mMJGBoiGHW2vhahGmaxml0GZehA5eQ5QIumw2CCHeAHA4gGwhRBmAOwEcLoQYieaBNxsA/BkApJRLhRAzACwD4AVwjZSt009PRvMsvklongQn6Ilw/GnbtsNyWTEOeGoSZu73BieRiiCP1DDs68kY/K9F0BsaMOB9YElyMpb1OBawqICUgNsDbWc1Brh+AxCbS/56123Ag8+eg5OmPoIMJcnscihBRXMWt+TwmMfPw+wxr8AuuGRwpHikhn2/vgaDH1oA3enEQP8cbqFp0HZWo7+r+fbxWM3hh/5zDo67hXPokHmiOYfbForY/EGPA60Tavo3Gvb4esR8o4HfXxQAIWV0f7ekix5ynDiqy9epPXtgzc2DMO3s53B4Uqz/FIeXJnWcs+5YzF9WgrknPBbUgZtTd2PYV1djyK3roVXG/5UmxeHAiseGY9XJz7ChRQAANW/NfCnlGLPriIRu53B6OtbcNgzTzmEOd8U/h3874THkBpHDNXojRn55HYb8JTFyWNjtWPnYCKw6hTlMzRIph4HmORvyp07p+oU8GSSiCNow5eYOszhumg0AACGgDB+M3s+V4aXin8NbWAw7e91RaJjkgLajEluvHI1brn0Lk1J3dHv48xdOO6a+eBmKH18IvaEhzNVGhrDaoPQpgHNQNtxpKjSbgOqWSNrhgWPlVni3bIWalYHlD/XD4mOfQqriMLtkMlkiHeQyh4135tqj0XiuHdq2Hdh++WhMnTIj4By++b+Xoeix36E3NYW52sgQFgvU4kI0DMmBO213MyFphweOZeXwbt0GtUcWc5haJVIOAwE0GygoQd/eEOfNnc5uCSECOm82hLL0ZfSREqJsGwYkbze7kqg1oz4D1bcVQSlbAADo9eRMvPbFYfjH5b1wzUmf45KM5e3eKlCvN+Gd+mLc9fVpGPRcDQqXzEQ8XLdUMzOw/YyhSJu0Bf8sfR8jbI1IFrbW52v0JvzS1Au3/XEacl9LwpBb1mBC/nn4Zfh7JlZtnK+cVnxUPQrX536HEoujzdVCp+5Gpe7mBH0UGCkhNm5BafIOsyuJWq/V9UTd1DygrHnJx5xnZuG1b5pzePKJX+KijMXtjjir0RvxTl0J7v1mIgY9V4OCOMvh5LO34s7+H2GMvX6vHP7SWYy7F57QmsOH5VyE+aNnmFi1cb5yWvHuzjG4Pvc7DLTamMNknDg/CTaCf4Mh6Nsb4u1z3NJX8H1cnd0SEpfzT/DnxlDxNbIBgHbEKLw+/cmghqTGu0XuJky+dQpSZ8xu93m1Vy4axvbFjpEWNOVokDYJ4RawVavIXqQhY/ZmeMsNmTDZfIoK56ljMPgvS/Dvgm+7dYWswluP05b8CUJIzBrxbgSK3NuvTTrG2qUhQ4grvPU48V+3oNf0RZCD+2LnPumo7SdgH1GNE/ssxQ53GlbX5OCJgW9iuI1XEPeUSFfUAs7hw0fhlelPII8nSHvpModzcuDcvy+2j7aiKUdrfdxWrSL3Dy/SftsEb8XWSJUbXkKgceL+GHT70oTN4TJvPU696xbkvL4QGNgXVSMyUNdHQB9WjzMGLsAOdxoWVubjv8NeYQ63I5FyGODIhojzO+lsc3WfJ6Nxaa/GiQCksvvrHpeNFYMkzsgGAJX7OtCTk/i167IlFyH7owUdXgnTtm2H4+PtKPq4/eejf4G07hF2OzZPHY2P/vwgSq2pALp3AJdnScXske/AqbsB2Lp8vZE0qePElSej+uVinHzz97gje0VI23Pqbhz22i3o9/xc6JoGoSjYcYgHDx86A8cnVyJZ8X18RUBHnx+P1DC9tgBPrDwCnvlZSBlXyUlaCQBQOdyBXDXZ7DKiUpc5vGMH7J/tQNFn7T8fVzl802jM+PMjGGZLQqzk8GGLz4T3lVyc9JcfDMnhI1+7Bf1e9uWwVcWOQzz41yHv4qSUit3Nl27msDYrC0mHMIeJDOM/qaP/1f14bDTEWANFqhJCF4bWvFcjQbb9une30cDbTtqKu2aDJwVceq0Dn494CYdPuQUF9880uxTTCIsFa+8ahTnnPYwsNbirrq0n4hF2VdGPGHr3NpRYHACCP5Cs8NbjsFdvQem9iyDS0rD6tsF47ewnMNbesoJA9z4+BQIOxQPXwiyUPr8W8lkdQx67HEsPn8bVCBIcc7hjn494CUdefQvyHmUON+dwcE0ps3L4kj4zccC96zHQakMoOVzmrcdRr9yC0vsWQSQnY83tw/DqJP8c7l7zpU0O/28t5IvMYaLu4kmhnxhqNITNHrePBL0Zfk+1EXdHg540/rR0JFtNweNXPAcxZh+zSzFNxTVj8f25DwV9gGsWVSg4NaUeA60pQV+x0qSOF2p646R7b0G//5sLJbsHyl7qjWUXPOV3gBtYTeenVWHx5U+iz8c18AwtxMDJazH0uytRr8fHhHVERstWU/DvaxI7h7dfsX/M5vBlGVsxzJYUcg6f9o9bUHKnL4enF2DxRU+EnMO9P2iAd2ABc5haSRuPiTuT6CeFUo3d7w+hGTOqQfqfCUsY03QR2N24oPhrNljrwvvVXe52YrsWuyswHJWkYdV1Ngi73exSIk7sNwz3XfciChPsPnKn7sZDO0sx5H/X4L0TD0D2c7OgZGZgx9MOLBob+nBbq1DxdMFs/GnaR6g7aggGXrUKo3+9EprkjW0UHszh2CX2G4a7p76UcDlcoze2yeEeL86CkpHemsOhjkKwChUvFP+Ci1/8uDmH/7wSI36YzBxOcMLNMx7qWKI3W4AwzcGwR9Milps6Roi7ZoO9Wob1l+vEWZNx+o034YilE7HKE5sHuz8d8QR6fJeMNf8+APKgkRCWuLubZi/CYsGam204MTnxrvTcs2MMfjhlX5T8dRa86zZAWCxYfk8pfh35pqH7OT+tCpf9631oIweg9PqtmLj6REO3T7Ej0jnc0b6i+UTrpyOeQOZ3KVj3wHjmcAK4a9tB+OHkfdrk8Ip/DAhbDuv7lGLQ1HLmMFE3ybg7I6J2tdNfCevXXrCpE3dHN9mLG1GuOcO2TJR1cQpS3pkJ8bEd1465BuWHJcMzvB490p1wWLxweqzYsT0dBw1ai1f7/hCWGkJVaEnF6yXfAyXfY9XpDTjj9yvQ+xEblF8XAlG+Okmw3IePwJcHPw4gsa6mAcAdOfMw7NYDMOi6LZAeN5qO2w/fHf8orML4z8Wf0rfjqwc3YtcpAvqf0zH19VF4JO93w/dD0S1nQQM2eJ2+CViN55/D1+03GeWHp8I9sh6ZaY2tr9mxPR0Di7fhs8EfReX8EYWWVLxZ8h1Q8h1WTWrO4byHbRC/LjC7tLBJ5By+u9ds7HPT/hh4Yzmk14um4/bDt6c8ErYc/uSRCjhP1ZnDic6ge9CD2Zf/zP3SKiE80X3CxVUGEoTce1UJfu3DK/qOwEJkWbwOj+44PGzbTz6wEsJqg3S5IH5dgML7ZqJk0iJknLgW9uM2IuukNejztoLr874OWw1GGmhNweJxr+Pv01/Glqnj4/PqmqKi7HJP2E58ol2yYsPMEx5F9aTREFYbmq6pRkkYPxcv9f0Km64YjG2H52BU6saw7Yeil7J4LR7fcUTYtu+fw5i9CAX3N+dw1omrW//0ewX4W8knUdlo2FNLDv91+nRsveFA5nAcSlZs+GniI6g+d/+I5PBr/T7HpssHY9thzOGEZtQ96EHsq7XRoEoIb3Q3GtrDkQ7xK6LNhfi8hhuQuPtR0urq8N3mAWHb/vR9XoY+ZsjeT0gJSAlLfh72v2deUBM9melQB/DV9Q9i7T37x92BrtqvGM/s/5rZZZgqz5KK/a//He4jhuOlodPDui+7sOLDyQ/imzsewflpVWHdF0UnvaEB32wYFLbtT9/nZcjRgzt83lKQjyEPLMGh3ZvQP2ocnqTji6nxmcOW4gL8e8wMs8swVaEvhz2H7huRHH578sP4/O8PM4cTnNknzUYvTxgpcXu1m5MXUoTFXbNBTUvDpH7hGy44xJaMjSd0MIO2EFh5YzHu7zU/bPsPpzxLKr4/7yFsuWEsIOIniWpG9cJ4e2PXL4xz9+f9CMftWzDEFv4Z4EutqTE30zwZR01PxzkDwpvDm4/p4KqwL4cfz58Vtv2HU0sOb716rNmlGGrX2Hwc4ag1uwzT3Z/3I3D7jojk8BBbMnLVlLDvh6JbuE+au2xmxGCjwZ/ZzRrDRXK0S5SLu69tlIq/T3Nhb5yUvjCsu8gYVdnuybjr+DF478zHYmLYbkcKLal46pqnIQ8cYXYphtk5WDFtTfZokqo48Nmgz8wugxJBYW+clhHee8RTD9jRbg67jx0dFzn82A3PQh400uxSDMMcbpaqOPDt0I/MLoPIMB02M+LkmlXcjnAgfm0jJHaPxjqg261IE96w7mNs7iYIS9vbJCx5vTH67vkYbouxcbvtONQBbJnqiYphvIrDASU5GUpKChSHI6gRF7qdLVyiSJI2i2k5vO+9C+Mihw9P0qMvh5OTg89hK3OYiMh0cdIEothh/lGMwbyZdqQo4f1JyrQ6IWzpkB5362NrJ5fgo96fIV76N++NnoZrx14NMTO8o0T2IgQsJX1Qdko+3AfV4ZiSFciy1sEqNNRpDvxY0R81c3LR5+MaYOFKSG/XJzQcJkWUGNZOLsFHeczhkPlyeMuEfDQeXodj+61AuqV5qWePVPFjRX/U/ZKLoi+7n8NEFHm6XYfwKGG9gitV2f68DOwvGmbP1RMgsHvVD1UGtLSiFIDg14YiKO6aDdbqJjToElDDt49GzQpoWuv/LXm9MeXM6FxeLVgDrSlYd0YSSmdGbp/qoP5YOTkbj5z4Kk5MroFVtPNF7LUAGAksuNiFSXMvR/HjavOBeCdLdiputnGJ4l285vCGU5JREq05fEX3c1jozGGiSFNc4c/DWJ0AMlB7nfBH0F779V/1I4BGQ7vbIgqz+Dkq8xENTdih28O6j12eZEht909r2dn9cEnGhrDu0wyjxq2GsEbgHltFxY6rxmPShz9g7dnP4tSU+vYPcP2MtNux4uBXcOcrL2LzX8d3Wmf6OgmP1Dp8noiMFYkcbtTa3kIRrzk88uBVEcvhnZcGn8NbbmEOEyWkBGg0AJE9SZdq8J/UUN5LFA5x12yQW7bh7erwzuL947r+rbdQqNk9cfIlP8MuYmupy+44Ped3KD0yw7oPYbdjwz/H4qO/PoSL0isDfv9BDgW/XvUwVj4xEkpK+7Nu95xfjaVuDvMlihRZvhWvVh0Y1n38uHYAc9ggLTn87j+Cz+Gfrus8h7PnVGK+K9RKiShgHFQUcwIdrRDQe7n0JUVY3DUb9IYGvPtr+JoNHqkh7dfdS1aVnz8If88J76zrZklTG8M7OZmiouzG0fj1Tw+j0NLBMnbdkKUmY+nJT2HtHcPbrVeuXIdrVp4bSqVxQ5M6NMkxdBReutOJT38ZHbbte6SGzB92TwIZzzmcrLggwrkUsRAov6E5h4uNyOG/DQeUvUdE6Gs2Yuqqs0KpNG4whymecF6sGMOlLynC4jIi+r/pwnK3Myzbfr+hB/I/LQMAWAoLcPqlP8Tl1TQAUCGDmnW8u5pOHI23r3oY2QasA56s2PDj+Q+h+tz993pOetyQL+WiUmsIeT+x7oGqIfhX1VCzy6AEUPp2U1hzuNfX5QDiP4fDzXXCGLwy+d+G5fCX5z+EGuZwp+7cMQL3VO5jdhmUKMJ8YpkwcwBEYkSAkdvnCIau8fMTEXHZbFB+W4JTf7vK8O1qUsftn5wL74ZNAIC1V/TBHdlLDN9PtGiSViBMs4yrPXug39+WY4gtuesXd1OeJRUTbv4Jak7OXs9lfLwIpy690LB9xaoXvz4C7/z3SLikx+xSKM4xh42hh/HXtNqzB0r+vgIj7cbNr1FiTcVxt/7MHO7Em18fjPdfOJw5TBRLIjEiwMjtcwQDRYm4bDZIrxf9/uXFu/Xphm73yV39MOjpbQAAtVcurjnz07ia+XxP3+waBn1XTVi2vfHPg/Fs0XeGb/ev2Quw5dwBez2uO51I/XsKptdmG77PWLHA5UL/NxuQP30pzlpzstnlUJyTXi/63ecx/Gfu8er+zGGDmJHDabc58EJNb8P3GSsWuFzoP6Me+dOXYuLKU80uhxIIb3cgf1IBr+y3NGMS/fMQZnEbPfqCZXj4zvOwyN1kyPbmu9yYcfcEaGvWAwC2nNMfV2auMWTb0erTxftCbzLm8+dP2WcwHrn0hbAMe7YLKwaevRKKw7HXc3LuYvz3L6fj28YwrosaxU7/eTLk/GXQdtXAdUsOnqzuY3ZJFOf0hcvx7J1nGprD7915TGLl8PwRMZnD/c5cDdHOiAl9wTL877ZTmMO7aoCpmcxhihijb3dItOZFzK7y0MGJtNDBkQ8UEXEdFelv/obzn7oJS92NIW1nlacBlzw5BWlv/QYAUFJSsP/5CyN+j3C15sRXTmtEhl5Wag0ofj883x4r/5yBCcnhm5b8jsJPgf59230u6cM5uOfqS/CfXUVh2380eqGmNwb/qw7QfUvPzVmMd285Fs/uKjC3MIp7aTOMzeGU9+YAMDeHv3DaUa8b3wDYU6XWgOJPw7PtcOfw34o+hdJFDj9Qtffoh3j2WHXfNjks/1jKHKaYlTBzNfiEskJEtEu4xtGeHy+bLmEV399eUiL/kVm4/K834oOG4GbZ/sJpxwX/dzPyHvsNkM3fjc4jh+He/K+MrLRbxvxwDR6bcBKOmHodBk6fjJKPrkT/H/4U9MfWmSvXn4rk75Yavl1LQT7+75j3Dd+uv0FWFbuGZ3b4vO3Lefj0tHEo/fYSw664RrMFLhde+MepELUNsBQWQB1YCmXkUNhqPHjylYlY5eGEbRRGvhy+4rYb8Vpdz6A28W59Oi68I3py+MljJuC4G29ozeGSj6+IrRzO6x32HN7HJlA9IqvD521fzsOPZwxPqBx+587jgB07YSksaM7iYYNgq/Hgv/8+BTPqM8wukRJI0CeX0Xi+7VeTtPGscS/d+JQkXONI7v4ZkKps/beepENPSrBPRgSEcV3DKCEl0l+fjWdXTsQdfwM+HP0cSq1dHxSu99Tj9AWXI/deK7LmzNr9hKJi12V1yDVg5u5AvFbXE4PuqoG2Zj3S1qxH2pu7n7vjtotw7LWPIVmxGbKvMm89tj1VitSG2YZsz1/N+CKcnfoBAGNqbY9dWFF7Rh0yZlggO5jgUlu5Bv0vVnHzAX/GmvMcuPiQn3Fmxnz0UnWoENAgUaNLfNEwBE8vPxR/3+czTEqrDlvN4ZSheHDMX3/G2JS1cAgP0pQmZCpuZCpAslCRqkT2e5kSkJRIe3M2Xlk5Af/6e+A53OseKzLnRk8Oe9dvROr6jUh9e/dzRufwpjDmcO0BfXBa6vsAkgzfdotu5fCqtbtz+FwHzjr4N1yUNQu91N0HezW6xId1w/H88oNw5/BPYjqHD/rbbyi9eztKbc1zjuSoDeitar7nbQC4ogpFRtAnl9F4Lu9Xk3BHYzekEwJR+TmVSgffIy2f3iisOWCqBHQBSAHh+3iES2n9NxlHSBndn9V00UOOE0cZsi1hscBz6AisO1fBFeN+wsT0BeihaHAIBU1Sxy5dwUd1I/DsnMNQ8oaE7cfFkB53222M2Qf3v/2CobN3d8UjNQyfdh2K/zmz3efV9HTsfDMXv46YEfJEaS7pweCPr8Gga3/v8AAxFBvuHo+Vlz1j+Hb3tNztxJRJVwGzF3Xr9Wp6OmTffHizkqDZVaguDWq9G8qGCmg7q+E8dSxue+h/ODE5/q/AUfeoeWvmSynHmF1HJIQlh89WcdH4X3F2xjz0ULTW53fpCmbUjMHLcw5CyVt6zOSwkpaGLa8U4vf9XzMkh4e8cx0GTp2XUDmspKUBJQXwZu1ugqj1bijrtkDbtYs5THtJpBwGAHtxkcyfOsXYje55Aimw98R50X2qsJcOT5aD2ZYq4/oWiq4Y+bmk+LVhys0dZnH8j2zwI71eWL6bj4HfAb+k98LPxZdAS7PDm2aFpc4DtcEFsakCA3fNa359O9tYe1ZaRA9wAeDGLQei5D8roXXwvFZbi55/TsclMw7HS8U/BH2g65IejJlzMYbcugJamJa8dPfs6KMw1hBbMtaekYLSji4KKs2TkykO39dSVaHUNcICwFOYBk+qDQ1DHdAPyoDmAJx5Om7+4ywcMf4Fw65cEiUi/xyem94Tc3w53EJtcEFs2IKBtXObX9/ONqIxh/W6OhReV2dYDg/++7K4yOH1p6SgpIscBpqzWFgsEL4cbixMg1QFZJYVtQcNgW4FpACun30ujjj8OeYwkVH2CFlpkRAe0e5zscLIk+N4aTR0p2nQ3muE3vy4TNKgNKiAAPRUL4RTjZvPDYVXQjUb/Gm1tcCSWgjsHrjYVTapOTm4YMKPYa6srXq9CbOfHYWelbM6fZ13cxmqJhXhsOfOxLf7vhXwpGnVmhOjvr4eQ6aubf7chIlICc/Bc3tOOnIuZv5pPBryBCAAZ7EXwqFBsWkYWrAVNsWLwuRdsCteDHSsQ6ltO/ItdchXVViFCgvUdk4YeIBLZBT/HG4RDzl8wH8m4ceRrwZ8QlypNWD/z6dgyF/WhDWHI+mE4+Zi7oIDUF+gwOsApAVo7OuGsOgoztuJHo7mOWMGpO3YK4cV37RSe38emcNE4dLaaIgDIV2Vj8AtDtImI3brh9CbG0kAILwCepoXSr0FUgAQzaM3pFWHcO3dKBc6IJy+5rAElLqEPX2kIPC7JQDl5w/AbdmfIZL3Vd5YfhRy317a4dU0f96Nm5E+KQMjr78Bd5z/Fial7ujy6ppHanioaijeefpIDH75D2hhWGLNn3RG7lvusbx58Nz7G6wikCXWOH8BUTSLhRzOOS8do6dMCSqHB70wH9oet43Essfy5sH171kBNsCZw0QUupBGOERgVEek55gQ3t37a2kYNM9R0Py40k6joVWMjnIh87HZ0E2Wgnycdel3EV1mzam78cd/h6NnbedX0/xpu2pQfNdMvP7qIbj33HyMP3kRrsj9EYOsLliFAg0SHqljoTsd/958LMrfKUH+O2uRs3VWl1cUOyOsNgi1OaREWhqEww49KxVamgOuLBsas1V4HcBBI5aHsJfABdZoIKJoFjM5XFvbJoeHn7ACU/K+xiCrC4poPqjzSB2/uzLxwIbjsfODwtYcDvV4riWLRUoyRHIyZGoSvD1SIC0CdUV2uFMFDhoZ2RyO9PKkRETt6Wykg6HzPCi+k3iJ3SMkfH9Lm+82FZ68U4LostkghHgRwEkAtksp9/E99haAQb6XZALYJaUcKYToC2A5gJW+52ZLKa/yvWc0gJfRPP31ZwBukNE+O6WftVf2xUc9P0EkVwvVoaPqAA96fV0M74ZNAb3Xu24Diu7dgC2POHBX0VlwFWVBtykQuoTi1mFfvwPalm3I9VTAKwQgBJSkJEAIQFGgZGYAVgv0ZAc8OcnwJqmoK7RAcwg0FEjoNglbYQOSHS4oAtgnuwJ59hoAwLDk1SiwVKOXWo8cVcIhVCSJ5mGvoU6cRpSomMUm5/BXRfBu3BzQe1tyeNfDdtxV3JzDLfeMCE3Cvn4HrOVbkOvd1JrDAKDY7YCqQqSmQKQkQ1pUuAsyWnNYKkBDAeBNaf6yKb2bkJ7mBAAM6FGJfsmV6OvYhIG2rchUGlFoab59LUNxQIFgDhMFiTkc29o0E/aY/HLPRoP/xJAtjYg9GxJSQestCHq6F3ApUFwKZJIG6RXN/xZos8IBGw2UaLozsuFlAE8BmN7ygJTynJZ/CyEeAVDj9/q1UsqR7WznGQBXApiN5mCdAODzgCs2gaUgHzed/UHED9BSFQfWTJiGu8bsi5nX7A/llwUdv1gIQCitk2wBgOiRCWm3wZ2XDs2uQHMoqCtUIRXAeUIhLM4iuHpoyOy7CxlJTRifvR6pqguq0LFf0mykK01IU9zIVyWsQkGSsAX4OUgO6eMnojZeRgJnsdk5fOvIMVh0w8hu5TDgN+GhL4el1QLdKiAVgYY8CzypAhheBCmK0FCsI31ANVLtzbdP7J+9ET2tDehnX4a+1kpYhYZ+Fm+QOWz3/SEiA7yMBM7haBDICIQ2Iwx8dIcOpUmBVCSgNJ/862lewKs0L33o27ZM1oBGtfnWA9E8NGHPZRGbX+trEtfuPqVSGtQ9XuNXAxsNlGC6bDZIKX/ydWf3IoQQAM4GcGRn2xBC5AFIl1LO8v1/OoBTESPBWnFyH1yc/iHMWANbFQr+mbMUM14swwMPnwd7rURdoQJpAdzpEt40HdKhIzW7AQUZNRiaUYG+jioAwFDHr8hRG5CvanCI5sm2ApuwTEU412Enou5L9CzedoK5OfxI3u+Y8eJaPPDweUjeoaO2WIXui1NXTwlPT29rDgPAPplbUGzf2ZrDVujItwioELALawC3eKm+P5wUkchsiZ7DACIycWKnu/cfWaBKQDY3AXSHDqGJNvMgSLsO7PFYywSYQhNomYinvQkP/R9rXXWBjQKigIU6Z8MhALZJKVf7PVYihPgDQC2AO6SUPwMoAFDm95oy32PtEkJcieaOLxxmXx0XAk1H15l+z+nZqTWYeOcTAAK9/9URnoKIKJoYnsXRlsO1RzUwh4komoX9mFjNyjK86Fjmv/Si0rj3iK/2HuNyjUSRFWqz4VwAb/j9vwJAsZSyync/2gdCiGEA2vvJ7rA/KKWcBmAaAKSLHqb2EdWMdJw3cJ6ZJbQy+0CbiKKW4VkcbTl80dA5ZpbQijlMRB0I+zGxvbjI/Gvr5ldARDEk6GaDEMIC4HQAo1sek1K6ALh8/54vhFgLYCCau7aFfm8vBLAl2H1HUuXEobiuxxfg7QREFI0SIYuZw0QUzRIhh4mIghHKTFtHA1ghpWwdCiaEyBGi+UZUIUQ/AAMArJNSVgCoE0Ic4Lun7SIAH4aw77ATFgvqzjkA5936OTIUHuASUdSK2yxmDhNRjIjbHCYiCkV3lr58A8DhALKFEGUA7pRSvgBgEtoOFwOAQwHcJYTwonnalauklDt9z03G7mV+PkcUT4Sjpqdj+YOD8fMJD6PQkmp2OURECZfFano6Vtw9FD+dxhwmouiQaDncwn8ZSCKiQIhoX9Y3XfSQ48RREdufmpmBTS8WYv64l3lvLhF1SM1bM19KOcbsOiKBOUxE0SiRchhonrMhf+oUs8sgImpjw5SbO8ziyC5YHuWE3Y7lDw3kAS4RkUmE1Ybl9w9iDhMRmY2DGYgoRGw2+Km8aBTmTniMB7hERCapvHg0fjvx38xhIiKzRffgZyKKAWw2+Kg5OZhwzS/IVlPMLoWIKCGpOTmYcO0vyGUOExEREcU8Nht8tp7RH3/P+d3sMoiIEhZzmIiIiCh+sNmA5uXVep5VxmG7RERmUVQkn7qNOUxEZDLJswMiMkiXS1+aTeuZgoYjxiGlvBHq6jJoVTu7flOA1N69cF2frw3fLhFRPNCzUtBwlC+HV26GVl1t+D4s+b1xa+kXhm+XiChuRGjCRqFHZj9EFP+ivtkwuGA7vn3iKWzTXPiwbhge/fk4DJjugpi9BNA1Q/bRNLA3xtmrAPA+YSKiPQ0s2p3Dr+4ajed/PQwDXjE+hw907ABzmIioA34TNkqFTQEiin5RP1BKQMAurCi2pOK6rI1Yf8o0PPz6c1j7ynCI/YYZso+aEhuylCRDtkVE0c8jNXzltMIlPWaXEhMUvxz+a/ZKrJ/IHCai0Hikhg8aUpnDQWKjgYiMEs5bp6K+2dCe4TYH1hzxEq6e8T4q/zweUNSQtudNFlBFTH4qiCgAHqlhWk0+hj93Hf59xpn4b00/s0uKWS05fMVbH2PnpaHncGz+NiKiQPnn8LRTjscT1YPNLiluGHHCIFWud0mUkMJ0m1bU30bRmVNSnBj214dwUu9bUXz3b0EP53VnGFwYEZnKqbuxU3djm2ZDlZaCla58vLJhLNxf5KDgg00o3jwTOoBHvj4R15z1rNnlxrQzUmsx8s6HcFJRaDlMRPGlXm/CLt2LNZ50uKWKla58/Hf1gVC+zkL+R805rAF45ttjcMtZa80uNy4YMdpBOnSIhhCbx0QUdfxvvdrzNqxwjpSK6WYDAJRaU/HZpQ/ijO23IOeZWUFtI3mbhCZ1jm4gimGa1HF1+UH4es5w5P0kkLrJCbXaCVHvhHQ2IqtmHaCvhtfvPYNerMFPJwOHOkwrOy605PCZFbcge1pwOUxEsa81h2eNQP5PEinljbBUVANeDdLZiN67VgBSMofDTaDN/A6BUNhoIIofLaMVJAAhWx+I5G1YcXF2XWJNxaM3Pwd54Iig3m+rlfCCV+OIYt3cl0ZiwLW/IXXGbGD2Imgr18BbvqV59YR2rrjri1bgsrcnwyP58x+qEmsq7r31ReYwUYKb+9JIDLhhNlLe/Q2YvQjejZt357Dc+wxYX7QCl82YzLkbiIiMJtHaeBRahJaz2UNcNBsA4PAkHVW3N0JYbQG/N3NJNVZ6eJBLFMtUoaA20CkYpMSAh1bhiMVnhaWmRDMh2RV8Di/eiYXuMBRFRBETdA7fvwyHL5oUlpoSUqjTLphzTkJEZhMw/Oc/5m+j8PfO8Bdx+SE3wPLd/IDeJ9duxINbJuDVvj+EpzAik2hSxwqPCw3Sgs2entjozoZTt+HnHf3h0iwQQqJXch2e6/MpMuJgJQDdEfi4MK2yChlXp+GS1w/Bf4t+5O1UIQo2h/W1G/HoluPwZsl3YaqMyBwu6cESt4QGgcVNRaj0pqFOc+DXHf0gZfNRXarNhZdL30G2GgdLvyqBn+lqu2rQ45os5nC04ByRFK98txhJVTZf6fe/zQC+x3URFz8DrR9jQG8yvo64ajaUWFOxbpKCgQEeq+pNTfjtp/0ANhsozrxd3xMvXXgSLJV1QGMTZIMTUtOAxgrYZfOJeW1mJu77bjwe6LXA3GKNEGQ31rtuA7afVYjxz56DH0a8hmQl8Cvz1KzEmop1Z6kB57B0ufDHd6OAy9hsoPjyRl0B3jrnKCh1jZC1dZCNTYCUsDVuan2NlpqKB348GA/1/sPESs3VksP7PXk+Zo/5H3PYCEHO3RDUSQpRtBOAFL6pC3yN3j1/PuKl0WCEPSeRDFbctY7PHzsbSkrgVwZKPm7EKk9DGCoiAlZ5GjClYgwmrT8SI+aci4E/XYT+r03GPrPPD+t9qju86VBXbIR33QZ4K7ZCq62F3tDQPH+BlICU0Kqr8cEn48NWQyQpjcFHmndzGXpO2opxT0zBk9V9OI9DCM4cNzeoHO7zBXOYwqclh49ZfjL2/e08DPzpIgx4dTJKv70krDlcpyVBrC+Hd90GaJVV0BsaoDudrRkMKaHX1eGjTw8IWw0RpQd/kurdXIaC8zcxh03GRgNFgrTINheJpDXMZ/ly98lzhyfRcdRoCPXn2KhJJOOu2XBKxu8QBb0Dfp8yZxkuWPIn4wsiAnDcF1Ow8iALdh1Wi96nLkfJpEUovWUWim93YYk7fMn25Y6hzQe1XSh9dQc+akgOWx2RoEkd6etC24ZeV4f8B2fii+P2wZgHr8OEFSfig4ZUrPfUo8xbj03eeix3O/FkdR9jio5TZ2XNYQ5T1GnJYeWYLcg/bRlKJi1Cv1tnYdA9dcxhgxiSww0NbXL40MWnteZwyx/mcDfF0YkTxR/hbTuKQHjY5IpmMsiuQVzdRgEAPRUXpMMa8Pukx42UJzKw/Dknhthi+5d9InFJDzxSwzbNix6Kgiapo04KbPBkYoeWDk0K1Ggp+HVXKZZsy4OzPBX23k7MH/9CxIZozm7SMGiaE3pT017Pyc1b8FLlIRhd8FtY9r1iXh+Ueiu6fJ22ai2mvnsxTrzwPzF7r+wGrxM58+sN2Za3rBy9Hy8H/mPB870OhZ6TCd2mQmg6hEuDKN8G4FtD9hWPmMOJxT+HNb/LVBs8mdjizQLQfHU/anN4Yxke33o0pvf5KSz7XjWzL0q6k8Mr1+DmGM/hVZ4m9JpVAyMuiLXksPDLYSl8y7Z5mMNEnZEWCehi99XpEJZDjSSjhu63u+1Euj2o5ZYRAz+XbbYVwPdT3DUbMhQBd05KUB+Y7es/cPLbU7Ho3Md5r2CEaFJHo3TDI3VYhYIdmhc1uhW10o6FjX2gQWD2rn6oakrB2vIcyEYVSqMKe6UCSyOQtlmHtV5DUlkd3DkpUJs0WHY1QtQ7IZ1NzbcLaBr0+loUeHcCAJSUFOz31xvwynlPYKw98BOiQD++8z6+BgN+b7+ZoDud+PbTA6FdMcvwg8u1nnqUfLD3gfWe1IGl2HBOLxxz1O+G7j/S/l52MpQlaw05yG0hvV54y7cA5Vua/4+Y+F0d02xf/4FTX5uKPy5mDkdKSw5v0TQ49d2/PZe687HTm4ptnnTM21kMTSrNOey0IHmTBYq7+d5X/xwW2u6fEFHvhKz1NQClvlcOj75lCv538eNhz2GP1DrP4aYmzP1yFLQrfghLDvf5rLHL16kDS7Hh7F444LAlhu4/0u4uPxFYtcHQbe6ZwwBzuLukKiGtEkpTcN/XCXVyFkWkKgEVgLe5WSCtEsIjmhsIsu3weD1Zg3ArzaME/LdhlRAaAHfXX79wnuAHKpx1BPK9HE2fk6BI39wUYdx+d8Vds6FGl7BVOoM74dA1DLhzIUYWXoHlh74Us1cWIs2pu+GUHmz0WrHWk4NVTXlQIOFQPHDqNiyqLUB5fQa2bO4J4Wr+nAq3QOpGBZZGibTNXqhuHbpVgb2yEUptI4THC1lVDQDQG+uheKsxQG7usAYdu7+Zu7rDU29oQN87ZuOOTy/H2Kd/xz25i0P/JHTg0eoBGHL/BnjbWVu8Rcnr2/Dp+ak4JaXrYbaBOHfJJejx27IO80DtlYvVN5bioTNewcnJtb7v99j8nndJD5a+PQS9nTPNLoVCpWvod+9CjOzHHA5Eeznsb0FNIeav7AvRpAI6kLJJhepqfk51NedwUkUDhNO1+001vollNQ2ycSsUqXeZw92lNzSg+J+zcMdX4c/hB6qGYfB966BFYw7n5GD11P7MYQoLoYmQmgVsNBjI70qwnuaFaFQhvAJ6ZvN8Mcqu3U1XaZGATYeQavMcKL73tff1UJxqu7vbaw6rTk4OI3ZSvcfKD9EsphsNUSbumg11uhXC7Q36/brTif5TKjBm2nmYPeYV2EV4r7iYSZM6VKG0+dslvajQ3HBLBdu0VGzX0jC/oQS/7eiL2iY7qivSodarSC5XAAFY6yUcO3Uk7fDAtqUGorYeel3zlSyhqs0Hqe46pHp3YqDs+kZOia6bBYaQEmLmQnxVPjisB7m/VvWHd9uOTl+jrVqLm9++GMdf/BSsov1fGoFa6m5E6uMZkB53u88rwwcj+7kt+LC4ZZ+xeXDb4v7KESh4Y01kvneoS3W6FcIT/Fcj0XJ4Ty05vMGbAQCo0lLxzvYxWLMzG1IK1GzOgOISSClTYHE2H7V1lMMtpLsWAz3zOq0l4sdWvhz+smxIWHN4XnUfaDsqO31NOHJ4kbupyxx2PFWFZaXMYaJo1jKyAAB0u97cxPGKNle/pU1CdGMUAQAodbtPv/ybDK2PuRTAtTsPWkctxMBJeqeiuP6YH8kQxeKu2fBZ3XCgfFtI29C2bUfexS4Mvv8azD3x3zGz7rVHavBIDTV684FNkwQ2etNRqzuwqLEYP+/oj7JdGXDWJCF9sQ32aglnbwFrPaBbAPsuCXuNjpRN9RBeHcrOOsDjgb6rBnb3ZuToGnI62X/MHWAIAZsa3qpzHPXYbLVAujrfT/9HV2Hi+JPx2aDPQt5njd6I0964CSVfz273eTl+BI7774+YkrUBzeP0YluFtx6fPXIYMrfNMrsU8vmwdj/DcnjIvddi9imPIjeGcrhSa2zNQ48EVnp6wi1VLGosxoebhsPpsrbJ4T215LC6vca3EQ+0qmrkeKoAALmd7D/mcjiK9H9oBY7d/3R8P+zDkLdVozfi7P9NRZ+v28+leMvh9Z56fP7QocjY1v7vHYoDMXLPv7+9TiD9Pgb/BkILPaX5lgSpyOYT/tbtSAjfZXnFreweaeC37S4bDTHwuQv3Cbe0yr0mhYxELd25HShuGg1ROHok7poNLyw8EP1rQ1+nWttVg0E3LMRxv9+M0679Hn/pudSwqx3BqPDW48fGImz29MDiugKs2NkL1TUpkACUTUlQnQJpGyXstTqStrsAXUJxe6FW1QFNLugNTqChAoV6WZvtZrWzr5b70uPl565DUsLlDe+PwI6mVEhP51fUAECrrIJ6ZT9MfXMUHskLfu6Eer0JI7++FoP/+Qf0doYMqwP6of+Ty3wHuLHPqbtx2Gu3oN8b86MpVxPey/MPxMDazq+id4e2qwYDp/yBExbcjFOvi84crtyRDtsmW/Pa3DqQvkEiY20jFLfvtF/XoW6tbp47psGJHvWr0aOT4fwtJIDgx+hRMLTqaiRdnWVcDv/r93Zz2NK3GCVPrIirHD7utVtQ8vocs0uhPRh68hjlv2SlKoEMD1BrbR0JIO064FKa713fo/72Vj5QGpp/vwi0fc6/8RDtn4dQhPuEO5DVJgyd2LCjRkMMNtC6FIUfT1w1G7ZrDSh607gPSbpcyJ42CzM/LMb+5xyJvmesxU1FX2KErRGpwh6xe4k1qeOg92/G4HvWAh439PoGZGlrkNXFASsPVrtWuTkTGBW+7W+uzUC23N6t12qr12H5pH44YVpvfDjw44BPqiq89Tj43Zsx+I7F7c64Lqw2rLk7FZ/lzw1ou9GqUmvAuHemYuA9i6B3MEyZIm+71oA+7xp3n6/0uNHz+VmY+VFzDvc+dSP+2vfTqM9h/0eZw52Lxhw+4ukifDH0nYBv4dnkrceRb92Cwf9Y1GEOr7qvBz7K/wCxftsE4JfD9y6CrnNcTbQJ61XqKBt2LjQB7Gw7qXAg8xZQdArp+6yrZgK/HyIirpoN56+ahKRvFxt+RV7bth29ntiOpufseKDgVDgHZKOmxApvioAzT4e1TsFzFz+NQx0G79jn79tHYshDm+Hd0fm9/xS47LkqPKdoYbtaWr+oJ7Llqm6/Xlu9DuqZPbDPjdfivrNfw6kpu7o8mXLqbty29SDM+fdoDHhrHnRv+6c2jceNxHcHPgogNZAPISp96nTg9meuxoAn5rHREGXCncN42oIHiv1yOLm5seFNAp69LHw5fGPFOOZwmOTOVOE6xRO2uTkafs8OOIeTzs7C8KnXB5XDpW/M6fDEu/G4kfj+4EehivjI4b89fjUGPMscTkTtnQDqqRqEU02oJkc8CfvnNoiRBCHVk0DNhGj+uYibZsMidxO8D/aC0tTxTNmhki4XvOs2wLZuQ9u5C4TAn3pchdVnPm34VbZ369Mxd8poKGWh3xoSsxQVlj6FqB/WCw15KjSbgOqWsDZIZC6rBVZvhN7QENSmc78tw/u39cDZqTUGF918havw28APwLSqneh7xyy88L+jcce5vXHgCYtwae5PGGRthOob2tckdfzm6o2H1h4L91u9kPPBCmRUz+4wV4XVhsbJ1Si0xP4BriZ13Pj2JSh5dBZkN4ajU+TMd7nDn8Neb8c5nBu+HF4+ZRiUsgWGbjemKCosRfmo3zcPDXm7m7PWBonMpTXAyvXtXsnvjuyfyvBhQ3bYcrjou66XntyTVl3dJodHHLcC1+V9s1cOz2zKx30rJ0B5r2drDnckLnP4KeZwVDJpeLhSH/7b3KL1hCqi/OeeMPAkMxY+t519vNF8wh1u0fxxi2j/JTFmhEPO+bKo09dUag047JlbUPivWYBJH4+lIB+lH+7AEwYOUX+hpjde/ttEJL/X/trg8U5xOFB34gjUXViLR/edgTH2+tZh05rU4YWGdR4PHt56LOa8NxzFr21oXos7EEJg9VNjse605wyv/5TVE+A5blfQB+AtFIcDoigfnvwMaFYFQgdUlwbrhu3Qtm2H7GAkgz+x/7545O1pGGZLCqmWaPHsrgJ8fNx+8G4u6/rFYfKNfGe+lHKMaQVEUHdz+PCnbkHBgybmcF5v5L7fgJeKfzZsm9Nq8vHK305O+ByuvbAW9wz7EIclVSFV2Fuf90LDSo+GhyuOxR/v7xO3OSzsdijFBfDkZ0AK30RxHh3WDdvhrdgGdOMWAjF6GB589wUMt4Vp+E2EMYcjz15cJPOnTun6hVFwL3oin/iFXWdf3yj42kdK0N9jCfQ5ipQNU27uMItj/4ZBAMs8KSj6qta0A1wA8JZvwaqrBmFaTX7I23JJDyatPxJvXXZcwh7gyvEj0PBhHt577FEsHPsGjkrSkKEktV6xVIUCu7BiiC0ZLxT/gt9veBJHfbkCVVeMh7Dbu9i6/44kBj+9Cx80GHulqVJrwI5n+4Z8gAsAelMTtNXroPz4B6zfzIflu/kQvy6At3xLtxoNAFB+RFrcNBoA4IqMzdhwQbHZZZCfJe40FH1RbW4OV2xFxbV9DMvh09ccg/cuOiphc1jsv29rDi8a+wZOSXG25nDLH7uwYrjNgel9ftqdw5ePh7AEMHAyBnJYulytOaz+8DvUH35vzeHuNBoAYMvhGRhmtXX9whjBHKbOGNpoEACMmwoo9uz5sXMeAgBsZsWKuGg2HOoAdt3VBCU52dQ65LwleOeyY/DXbcOD3sZrdT2x/6M3oPYUQMxcaGB1MUII1J99AK6d/jZ+2vf9bi93ZxUqbuqxDl/838NYff9+UNPTu71LbelK3PnURdiuBXcrRnsOmX0VMt6LnltfGkc6zS7BUKpQcPBpf5j+M0+7HZ6ko/JezfSvSUsOT60IfsbBF2p6Y/9Hb0DTaRrk3MUGVhcjfDl86WsfBZXDH//fQ1j94JiAc/iuRy80LIc1qWP8z9dEVQ47RzsjNqFpJKhCwdiJi6E44mOkRlyJphPOUJoE/sv4RdPHFCuC/NxL1dhPtoy22IvH76UobsZF25c/aJ/vOx3VpwV/km8UMXMhFpxVitLvLsEqT/cPmiq89Rg17xy8es5xyHtkJrSqnWGsMnrVnzUO/3ngcZySEtzJcbaagsVnPYEVTwwI6AAo7z9zcOBbU+HUQ5/k6rz1R6DfLTWQLlfI2zJKjwzjGindpcnwtpyvz/0O6M+ratHkqxEvR00OrzijEKXfBpfD75x1OHP4gceDnkMhz5KKxWcHnsO5z881LIcv3HAUBv6lkjkc5hy+sffXQP++Yd0HUcIL9uQ4yPd1uFRkkDgCIQKiuIESN82GLDUZPS7fFNgQ+jDR1qzHgD8txuRLrkfpt5fg3fp0VGt7nzzX6034ymnFsFnnY9I1NyH3rA3QFywzoeLoIMePwF/unY6RIX4NkxUb/jj6Kaz7+36A6F5gSq8XA/65FCNevyGkK2tnrzsK1ZN7wbthU9DbMJwQSLZ6wroLTeqo0Rsx3+XG9Vv2x4BXJ+PN+pyu3xiC/lYLagdnhnUfFJhoymHvhk0YcMnuHJ5em91hDn/QkNo2h5esMKHi6CD23xc33POGYTm84dZRQeVwhbc+qP1qUseZa49uzmET5xIwg38OX7bp4Ijk8CCritqhmWHdB8W4PU+COouDQG4XiHcG3zpi9GiFWBbJkRZRN6rDBF3eVCmEKAIwHUBvADqAaVLKx4UQPQC8BaAvgA0AzpZSVvveczuAywBoAK6XUn7pe3w0gJcBJAH4DMAN0sAZKv/Tbwb+PPqaqLj9QHq9sHw3H/2/F3ix4BA8NSyvdZk2qQLWeom0zRpSl2xF4aYVgK4ldKYqycloumtX0CMa9pShJOHt8/+NKd9fA+s387v1Hr2uDqV/nYvjNt2Cv1z3BialVXd7f2Xeehz203UY/LdK6Buj7ERFSjg94VlSrsWNFeOw5K8j4NhcA7mxHAOSVkI7JbxjuixQ4UmK4nFjBoqlHH6s39u4fvhVQBTcfuCfw2/1Go2XRpyKmhIrpNL8fSOkL4cXbkFhOXNYSU5G0321hq0KkaEk4a1LH8WUmYHn8Elrb8EtN/1/e+ceI1dVx/HPr/to6RsohWVLbYFWLYFYqAQV0ESjQJBifaREpYhJNbFKgyZUG41/+YzEGAwoAQGtVgg0VCOmUJ4iim1pKbCUtraWbrdbBLSLFLoz+/OPOdvObvcx07l37rnd7yeZ7N0z9/Gdc2c+O/fsueesrNrDFz98He9e3kFPe2QergN9PLzjZWaNq4+HC2NGhochXy6G+Adp9EYHA+seOe+hikn4j9HR9laoaRBGiLLBqJ6fiZg/f/WikhGcCsDX3X2DmU0A1pvZg8A1wFp3/4GZLQOWATeY2RxgIXAWcCrwkJnNdvcicDOwGPgbJbFeAjyQ1IuZ2TSeve8bR8tfk9pjArhT2N1O8+52Bvr/QmXD+x377L/8HO5/10+Ayu4NroRzmsfw3+u6OOmxZrzCOcC9UGDqTX/lrrUfYfmXTuTbl9zHVRPaB5z/veg9PN99kMUvfI7mX5zArD9tpBDpXOP/fnVCqvv/4wtnM2vNOnqHSWs8fjKTG9IfJ6LnmJm8d1hy4+EzGo9jz0UTaEluYp7acaewt5PmvZ3y8BDE5OEpv3iKux4/7OFPjd/F+FFH3pLR38Oz5eFDHm4YpzFtUiA3LoZIL3TKLj6tYIM+J+LgqN9DOpeH6O1VkvTtKUccJ9LGxWG/qrt7B9ARlrvMrA1oBeYDHwqr3Qk8CtwQyle6+9vADjPbBpxvZjuBie7+FICZ3QVcScJi7ZpzkJYkdyhSxxobOXj1a0ypcBCyavjdOb9iyflfwZ7cWNV2xbatnLl0K/dMey+3fHABr5wHja1v0tRUpLu7gYP/a2bi5mZantjP5I0v4oVtUXt1wjNjKH64J5XByYrew6Sn+14E+Ftv8Z/iWCDdBoeGOK8pEidPHm6wUXSdJQ/nDWts5K3PR+rh1nnc+oEF7L2g5OFeDv6vmYmbmjn1sf8w+bmXovfwxH+k5+FuL3L8XzK6fSnmSk+YPLlYiDQov6A94uI2J1NKpn1R7g3ep2Eh7UaGQ8eJsKEBKuvZcAgzmwHMBf4OnByki7t3mNnUsForpVbaXnaHsu6w3L88URqO0/+o8saomdP56ZyVpDGEyOymceyYfxynP3l02xd2tzNpRTuTVgz8fA6cCsCpD7/O+q8VOX908nX85wNjaflzB30mf+su0FmYBPw78eP18nrPAca3j5DWhjLy4OFRoyubClDEQ8NprXz/3atS2ffspnH869LjmHG0Hm7fw/i793Dm3QM/H+n3qyNoeeQ11l+fjofXHBjHyQ/vlYfrSB5cLHJIyrcf9L8QrpbyC9ojLm5z8qU47YvyejUu5IWK/+KZ2XjgXmCpu+8fatUBynyI8oGOtdjM1pnZuldere5L6yjLyTtdHKLr7JOY25xeI9EZ83ZVN+f7MUjP5pf4woZrUtn3koeuprh9Z5+yYlcX9+6am8rxenn4wKmMaWtP9RixkRcPi/yxf+4pXDgmvdkSZr5fHu55YWvdPXz3zvNSOV4v979xBqOffznVY8RIVi4uvlHdZ1SD0+WQlKf51IVwRqRc7T4q/WMcLRVpyMyaKEl1hbvfF4o7zawlPN8C7Avlu4HTyjafBuwJ5dMGKD8Cd/+lu89z93knndhQ6WsBoFiQWfNGV2sDY0c1p7b/j05tY9TYEX7vak+RU34+hraDyd7WsPZAA++6uQv6j2nlzhuPTU112rVvPr2Awt7O1PYfG3nycF7+uyEOIw/XgQw8/PYjU1L18PfWX0px3yup7T9GsnRxw/jqbnOKtVu1ECMNT7khwHqI9rvXsFfmZmbAbUCbu99Y9tRqYFFYXgTcX1a+0MxGm9lMYBbwdOhe1mVmF4R9Xl22TWI07qp8Tm8RB93pjpnFpIY3oWlk/0cNoPHRjVz2p6WJ7a/oPSz59Zfo2dQ24PPTV3Xy6FvpzILRUXiD6b9pPPLL9TGKPCzyTpOpdwyUPHz5qusTawB427v56u2De3jaH+ThJMmbi0X+8OZ0P08x/we83iQ9HehQPYlGcsOfDTfLjpldCDwBbObwrZHfonSP2t3AdGAX8Gl3fy1ssxy4ltKovUvd/YFQPo/D0/w8AHx1uGl+zKwL2HIUr61eTCHNGyKTIfaMylc7sWc8FvO9w93TncQ+IA9XxLH4HqsnseeD+DPGng/iz1htvrp5GOTiCoj9/QXxZ1S+2ok947GYb1AXD9vYkDVmts7d52WdYzBizwfxZ1S+2ok9o/LlmzzUT+wZla92Ys8Yez6IP2Ps+bIm9vqJPR/En1H5aif2jCMtnwY4EEIIIYQQQgghRKKosUEIIYQQQgghhBCJkofGhl9mHWAYYs8H8WdUvtqJPaPy5Zs81E/sGZWvdmLPGHs+iD9j7PmyJvb6iT0fxJ9R+Won9owjKl/0YzYIIYQQQgghhBAiX+ShZ4MQQgghhBBCCCFyRLSNDWZ2iZltMbNtZrYsowynmdkjZtZmZs+b2XWh/Ltm1m5mG8PjsrJtvhkybzGzj9Up504z2xyyrAtlJ5jZg2a2Nfw8PouMZvbOsnraaGb7zWxp1nVoZreb2T4ze66srOo6M7PzQt1vM7Ofhfmy08r3YzN70cyeNbNVZjY5lM8wswNldXlLRvmqPqdp5Rsi4+/L8u00s42hvO51mBfk4oozRuvhcLzoXBy7h4fIKBfXlk8erhJ5uKqc0bo4Rg+HY0TtYnk4tYz1cbG7R/cAGoDtwOlAM7AJmJNBjhbg3LA8AXgJmAN8F/jGAOvPCVlHAzPDa2ioQ86dwJR+ZT8CloXlZcAPs8xYdl73Au/Iug6Bi4FzgedqqTPgaeB9gFGaJ/vSFPN9FGgMyz8syzejfL1++6lnvqrPaVr5BsvY7/mfAN/Jqg7z8EAuribjTnLg4bLzmrmLB/FINB4eIqNcXEO+fs/Lw8PXoTxcXc6d5MDFROLhcJyoXTxIPnm4xoz9nk/NxbH2bDgf2Obu/3T3g8BKYH69Q7h7h7tvCMtdQBvQOsQm84GV7v62u+8AtlF6LVkwH7gzLN8JXFlWnlXGDwPb3f1fQ6xTl3zu/jjw2gDHrrjOzKwFmOjuT3npE3hX2TaJ53P3Ne5eCL/+DZg21D7qnW8I6l5/w2UMLbGfAX431D7SzpgD5OLaiNHDEImLY/fwYBnl4mTyycMVIw/XTowujsLDEL+L5eF0M6bt4lgbG1qBl8t+383QQksdM5sBzAX+HoqWhK47t5d1LcoqtwNrzGy9mS0OZSe7eweU/kAAUzPOCLCQvm/kmOoQqq+z1rDcv7weXEupRbGXmWb2jJk9ZmYXhbIs8lVzTrOsv4uATnffWlYWSx3GhFxcOXnxMMTt4jx5GOTiWpCHKyNrXxxBxB6G/Lg4Zg9DvlwsD9dGqi6OtbFhoPs/vO4pAmY2HrgXWOru+4GbgTOA9wAdlLqeQHa5P+Du5wKXAl8xs4uHWDeTjGbWDFwB3BOKYqvDoRgsU1Z1uRwoACtCUQcw3d3nAtcDvzWziRnkq/acZnmur6LvH/lY6jA2onr9kbs4eg9Drl0cnUfk4pqRhysjqtcfuYchBy7OsYchMo/Iw4mQqotjbWzYDZxW9vs0YE8WQcysiZJUV7j7fQDu3unuRXfvAW7lcJemTHK7+57wcx+wKuTpDN1deru97MsyIyXpb3D3zpA1qjoMVFtnu+nbbSv1rGa2CLgc+GzowkToivVqWF5P6f6v2fXOdxTntO71B2BmjcAC4Pe9ZbHUYYTIxRWSEw9D/C6O3sMhm1xcA/JwVcjDVZATF8fuYciBi+Xh2qmHi2NtbPgHMMvMZobWv4XA6nqHCPew3Aa0ufuNZeUtZat9Augd2XM1sNDMRpvZTGAWpYE00sw4zswm9C5TGjDluZBlUVhtEXB/VhkDfVrNYqrDMqqqs9CtrMvMLgjvlavLtkkcM7sEuAG4wt3fLCs/ycwawvLpId8/M8hX1Tmtd74yPgK86O6HuoLFUocRIhdXli8vHob4XRy1h0EuTgh5uHLk4coz5sXFsXu499jRulgeToz0XewJj2ia1AO4jNJIt9uB5RlluJBS95BngY3hcRnwa2BzKF8NtJRtszxk3kIdRkumNDrxpvB4vreugBOBtcDW8POEDDOOBV4FJpWVZVqHlCTfAXRTaqn74tHUGTCPkkC2AzcBlmK+bZTu8+p9L94S1v1kOPebgA3AxzPKV/U5TSvfYBlD+R3Al/utW/c6zMsDubiSfNF7OBwzKhcP4pFoPDxERrm4hnyh/A7k4WrqUR6uLGP0LiYyD4djRO3iQfLJwzVmDOV3kLKLLWwohBBCCCGEEEIIkQix3kYhhBBCCCGEEEKInKLGBiGEEEIIIYQQQiSKGhuEEEIIIYQQQgiRKGpsEEIIIYQQQgghRKKosUEIIYQQQgghhBCJosYGIYQQQgghhBBCJIoaG4QQQgghhBBCCJEoamwQQgghhBBCCCFEovwfogZgb6PLyS0AAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 1296x432 with 3 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, axes = plt.subplots(ncols=3, figsize=(18,6))\n",
"\n",
"axes[0].imshow(dilated_footprint)\n",
"axes[1].imshow(dilated_isotropic)\n",
"axes[2].imshow(dilated_footprint.astype(int)-dilated_isotropic.astype(int), vmin=-1, vmax=1)\n",
"\n",
"print(\"IoU:\", iou(dilated_footprint, dilated_isotropic))"
]
},
{
"cell_type": "code",
"execution_count": 51,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 10 s, sys: 8.32 s, total: 18.3 s\n",
"Wall time: 18.3 s\n"
]
}
],
"source": [
"%%time\n",
"\n",
"eroded_footprint = footprint_erosion(image, 75)"
]
},
{
"cell_type": "code",
"execution_count": 52,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 220 ms, sys: 432 ms, total: 652 ms\n",
"Wall time: 650 ms\n"
]
}
],
"source": [
"%%time\n",
"\n",
"eroded_isotropic = isotropic_erosion(image, 75)"
]
},
{
"cell_type": "code",
"execution_count": 53,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"IoU: 1.0\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABBsAAAFpCAYAAADdkCCpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAA3WUlEQVR4nO3deZzcdZ3n8fenqo900umkQ+6LHAQE1AGJiHiMLoMgHoCOY1ARjzHCwkhGmRWPVWd33R3vXXYUBpUNKIIRZIgKyCHKKEdIQiAkIUnnbtIkIQd0p9PVdXz2j/4FKkkfVd2/qt+vql/Px6MfXf2tX1W9uzq+KT/1+/3K3F0AAAAAAABhSUQdAAAAAAAAVBeGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAEJV9mGDmZ1vZuvNrMXMri334wPAcEcPA0D06GIA1c7cvXwPZpaUtEHSuZJaJT0p6RJ3X1u2EAAwjNHDABA9uhjAcFDuPRvOlNTi7pvdvVvS7ZIuLHMGABjO6GEAiB5dDKDqlXvYME3SjryfW4M1AEB50MMAED26GEDVqynz41kva8ccx2FmCyUtlKSkkmeMVFOpcwFAUdq1/0V3nxB1jkGghwFUhQruYWkQXWx1dWfUTppY6lwAUJTuHa19dnG5hw2tkmbk/Txd0s6jN3L3GyXdKElNNs7fZOeUJx0AFOhBv2Nb1BkGiR4GUBUquIelQXRx/cwZPvULi8oSDgAKtXXRNX12cbkPo3hS0jwzm21mdZIWSFpa5gwAMJzRwwAQPboYQNUr654N7p4xs6sk/V5SUtJN7r6mnBkAYDijhwEgenQxgOGg3IdRyN3vkXRPuR8XANCDHgaA6NHFAKpduQ+jAAAAAAAAVY5hAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCVRN1AACIlJkSI0dKiYSsvk7KZGTNYyWzV7fJ5eT7DsizWclduUOHJI8sMQBUF3oYAKoSwwYAw1rN9GkafXunTm/aofE17WrPjdAbGx5WrWVf2abLa7W8c47SntRLmQY99uJs6W8iDA0AVYQeBoDqxLABiFLeuzaWTMpzLuWyUiIpS5isoUFKp+WZTM9GyaQSjaNkI0Yo19wkb6hV57QGdY9KKDU2oUk3LpenuyP6ZSpT5vk2bT4wR7fP/kPeakJHH2X29hGbX/1h0jNKliUdgJKjhyNHDwNAdWLYAPQnkVRyTJNyc6dp15lNOvDajOpfTCqRNh2akVaiIaNcOqkR2+qU9waMkl1SMiV1j5G6pmRlTT0vPD1nqttWr3SjSxNSOm5chzLZhNxN847bo92do/X8i2N16rQ21SUzOnV0m/ZnRmrnoTGSpOa6Tr256VmNS3boNXV7NMJc4xN1Sprp5Luu0kRe4BYvl5WWjFfnX3VrZKIu6jQAjkYPVz96GACqEsMGVL5EUsnmMbIRI+QN9cqOH62aFw5I6cwxm/rY0cqOrpckpRtr1Tm5VqkxCXVOPvLAz+7xWdWN69LM8fv17smrdUHjbzW3pkFJK/ycqlnPFbV98Ua9cmlVKqV5Pz9UwseqbuPv26Rbvjhbl499PuooQGWih+nhIaKHAaD6MGxA6ZjJ6nreofDubslffSFpNTVKjB0jq6mR6o96FyORUHpqs7Ijknp5Zp26xpsOvq5LjU09L+K6DtUp+0KDcvU5NU7p0PjGg/rwtOWaVfeixiY6Nb3mkFrSTery2mMiTa15SRMSPS9+R1hCjYl6JWQFvBgdNcD1xyrtC9wjffzpT2jKsjVle7xqk921W99++L26/OJ/izoKEC56uOjbDBY9PDT0MABUH4YN1Sr/DM6SEvX1PcebJpPHvOA85naWeGWX1czoOuVqE2qfXquO6aauyRmNm3FAL7ePlFob1LjdlDzkGt2aVsfUWmVGmtpn59R04n7V12b0+uN2SpKe2TtH+9tHKnWwTpMmvaSpjS/pPRNWa2yyUxNqXj4iQp2ymlPbpVqZmhIjinyxmJTUqOk1OUmpXq4fUcR9VYbd2YMa9+PGnt1QMWizf53R9vd1aGZNY9RRUC3oYdHDKAY9DADVhWFDqSWSSjaOkupqpeOalWtqUGpcvbzm1RehNQezqtvVLjuUkhIJeX3wTlA2J+vskqfTsoaeXVNTU5q0+w316pibkfJex1q3SUkp0WXKjspp9OR2uZtqEjmNG9Wpk8fsUmvncWqu79T2jknasadZ6Zfqj8xal9O0qfvUWJfSuyat0fsaf6MxCVOtJdRo9X2+2Mx6Th2e0kirU631cbqm6cU8aUkN5h2s4Wrh5g+o4eE1ykUdpMLV/WWNrtl+oZbMeSjqKAgbPdyDHi4Zejgc9DAAVBeGDSFLjBolmzVd+05v1p7zU3rdjJ16y3EtmlTzkk4b0arRltG4ZFK1eedQ7vC09mQTeiHbqBGW1thEzztBaU9oe6ZZXV6r45Idmpzs0KRkTs2J4o5ZHbzC3llIWkJjrKHEWdCbzly3dtw6R+M7H4s6SsXLdXVp3dI3KHv1A2Xd9Rrho4dRTvRweOhhAKguDBsKYaZEfb1sTJOsrk7ZiWMkM3VNbFB3U1IHpyTUPUbqnntIF53ytC477madVJtUvR19rGrvu46OVJ0mJqVTlVXPxzy9+oLxtPrOvC15IYkj3XDgNZp0V4vYcTccM37dpvsWjtR7RnZFHQVHo4cRU/RwuOhhAKgeDBsOC06idfiEWbkJY3Xg5CbtO8U08vX79dfTW/S20Y9rdOKQXlO3X5I0JpHsY5fV6jseFfGT9Zx+fMf5mrnn0aijVI3spq363OOX6D3/6f9FHWV4oodRYejh8NHDAFA9Km/YkHfCrUR9vZRMKtE8VqpJymuSUjKp9MTR8mTPdt1jatQ5MalUs6lzck4jjm9XXc2x7z/UJHN68+Qtaq49oDePWq45tft0fE1dL++KSYXu1gqU0n2HRmr2bbt5Ny1M7pq8tF4d7+hSY4L/s9onehiQRA+XBD0MAFUj9sOG9ORR2vKPb1amMafG6S9r9IhXz2x9+nHPa3L9Xs0fuVyjE4dUa1nVKaepNRnVBmftqrWEGqxuEMf+cWIsxNtVf/yYTtywIuoYVWfMQxt03b6/0pfHr486SmzQw0Dv6OHSoIcBoDrEfthw6sQ9Wvbx6wvYMhF8SVJ9fxsCFe+Z7i6d+JN+PjoPg5bdu0833f9OffkjvMg9jB4GjkUPlw49DADVgVP9AhXosqc/IVu+NuoYVWvukoNa19058IYAhi16uLToYQCofAwbgAqzP9upUT8bI89koo5StezpDfrHzR+KOgaAmKKHS48eBoDKN+hhg5nNMLOHzWydma0xs6uD9W+Y2fNmtir4uiDvNl8ysxYzW29m54XxCwDDzdWt52v0vc9GHaOqeSql3UtmKuXpqKMMiC4Gyo8eLj16GAAq31D2bMhI+oK7nyzpLElXmtkpwXU/cPfTgq97JCm4boGkUyWdL+lHZsd8VhmAfnTmurX6Z69V7uDBqKNUvSm/2aY7OiZHHaMQdDFQRvRw+dDDAFDZBj1scPc2d18ZXG6XtE7StH5ucqGk29095e5bJLVIOnOwjw8MRz88cLKm/Gpj1DGGhczONn31zxdHHWNAdDFQXvRw+dDDAFDZQjlng5nNknS6pCeCpavM7Bkzu8nMmoO1aZJ25N2sVX0UsZktNLPlZrZ8z14+vRqQpKzndNOS85TdsyfqKMODu2bf5mrLdESdpGBhdjE9DByLHi6zYd7Dwf290sXZDvamAVBZhjxsMLNGSXdKWuTuL0u6XtJcSadJapP0vcOb9nLzXj8vyt1vdPf57j5/wnHsVQZI0gOHGjT79l1RxxhW6h97Tt944dyoYxQk7C6mh4Fj0cPlN5x7WDqyi5ONo8IPDQAlNKRhg5nVqqdUb3X3X0uSu+9y96y75yT9WK/uFtYqaUbezadL2jmUxweGkyv+42PKbtwcdYxhJXfwoP581+nKei7qKP2ii4HyoIfLjx4GgMo1lE+jMEk/lbTO3b+ftz4lb7OLJR0+XfNSSQvMrN7MZkuaJ2nZYB8fGE42pTt0wk+zkvf6xgdK6Pile/WXVHw/JZguBsqDHo4OPQwAlalmCLd9i6RLJa02s1XB2pclXWJmp6lnd7Ctkj4rSe6+xsyWSFqrnrP2XunuHAgMFOBvV/29Jj2xtvd9LFFS2XUtunzlx7T27J9HHaUvdDFQBvRwdOhhAKhMgx42uPuf1fsxZ/f0c5tvSvrmYB8TGI5aMx0ae/1oebo76ijDUy6rcbeP0ktnHdKYREPUaY5BFwOlRw9HjB4GgIoU333SAEiS/vOWv9WIR9ZEHWNYa3p4o/7vvtOjjgEgIvRw9OhhAKg8DBuAGOvIdWn3j2cp19kZdZRhLbt3n265551RxwAQAXo4HuhhAKg8DBuAGPvqrrdq3O/WRx0Dkmb/+0Gt6+b/bADDDT0cH/QwAFQWhg1ATKU9qwd/daay+/dHHQWSbMVzumLDR6KOAaCM6OF4oYcBoLIwbABi6oYDc3T8z7ZGHQMBT3er45dTlPJ01FEAlAk9HC/0MABUFoYNQExd99sLlHl+Z9QxkGfS77frl+1TBt4QQFWgh+OHHgaAysGwAYihx7uyOuFWdtuNm0zr8/r6Hz8QdQwAZUAPxxM9DACVg2EDEEMfeWShcqs5IVkczb09o+2ZjqhjACgxeji+6GEAqAwMG4CY2ZLu0JybXXKPOgp6UbvsOV27431RxwBQQvRwvNHDAFAZGDYAMXPRU59RzZ+fjToG+pDr7NSaX50cdQwAJUQPxxs9DACVgWEDECO7swfV/ONGebo76ijox7RbN0YdAUCJ0MOVgR4GgPhj2ADEyKLt71XDw2uijoEBZPfsiToCgBKhhysDPQwA8cewAYiJzly3nvv5a5Tr7Iw6CgAMS/QwAADhYdgAxMSXd52tybc/F3UMABi26GEAAMLDsAGIgbRn9eCvzlR2P5/pDgBRoIcBAAgXwwYgBn760kwd/7OtUccAgGGLHgYAIFwMG4CIZT2n797zPmWe3xl1FAAYluhhAADCVxN1ACAMKU9rVzalg7mE1nZPliTNrd2jccm0xiVq1JgYEXHCvv2us1En3rhH2aiDAMAQ0MMAACAfwwZUpLRn9fChEfrW1ndr26qpmrhcamrpkHVlZHv2STmXJjQr11CrgzMbtfuMhJpO26tLZi3Xx8as1sTkqKh/hVdc/YeP6sQNT0YdAwCKQg8DAID+MGxAxejMdeuRrtH62vr3q+sPEzT97jbVbN6uubntkiQPvl4RfAb3yBXSrLskmenBscdr6dnnaNv7TJ8++xEtGrcq0nfbVqS6ddJPuo7MDQAxRQ8DAIBCMWxA7K1IdeuKtR+V3zFeE/+4U83bNku5jcXv7uqu7P79qv/dkzrxd9Kjx03V0ve9UycuXKcfzrxXYxINpYjfr08+fZmmrFhb9scFgGLQwwAAoFgMGxBby1Jpffj3V+qkfzuocavXyzMblQnx/rN796l58WPae+donfvBz+vMK1fqB1MfVa0lQ3yUvr2YPagxNzdJOY4SBhBP9DAAABgsPo0CsfTZ1jfry5/6rE684kn5U2vkmTBf3h4p196u5sWPafN7x+g1v7xSm9IdJXusfFdtf58af7+6LI8FAMWihwEAwFAwbEDsvHPNhdrx0SlKPrxS8vIdRZvdtVsnXPOkPvLVf9LvOkt7/HBnrlsbfnaScp2dJX0cABgMehgAAAwVwwbEymdb36yGKxLKbtwcTYBcVmN/9ph+8JmP6Dv75pbsYf55z5ma/MvnSnb/ADBY9DAAAAgDwwbExn2d9dqy6ERlW7ZEHUXJh1fqoUvP0g8PzAj9vtOe1W+XnK3s/v2h3zcADAU9DAAAwsKwAbGwPdOha6/7tOzRp6OO8gp/ao3uXHSe/ngo3P+ZLD3YrFm3Px/qfQLAUNHDAAAgTAwbELmUp3XO7f+kSf/6RNRRjlF7/3J97fN/r5++NDm0+/zibz6izJZtod0fAAwVPQwAAMLGsAGRu2TTBZr3v56L7UePNdy9TEsuO1dfaHvDkO9rWSqtE37eHkIqAAgPPQwAAMLGsAGRast06MXvzY7/cbPLVmvN35885GOHP/zHy+Wr1oYUCgCGjh4GAAClwLABkfrwuks18r74HB/cH39qje74x/P0eNfg3vlbkerWST/qLuvHyAHAQOhhAABQCkMaNpjZVjNbbWarzGx5sDbOzB4ws43B9+a87b9kZi1mtt7MzhtqeFS2F7MH5TdMlKdSUUcpWN39K3TZbVcp7cW/0P3Eqk/IV/BuGsJHF2Ow6GEgHPQwABwrjD0b3unup7n7/ODnayU95O7zJD0U/CwzO0XSAkmnSjpf0o/MLBnC46NCfe2Fc9T4YIW96HPX3Os26arn31rUzTalOzTh+pGxPR4aVYEuRtHoYSBU9DAA5CnFYRQXSro5uHyzpIvy1m9395S7b5HUIunMEjw+KsD+bKee+MnpyrVX3km6srt2a/3XX6sVqe6Cb3PuHz+nuj+sKl0o4Fh0MfpFDwMlRw8DGNaGOmxwSfeb2QozWxisTXL3NkkKvk8M1qdJ2pF329Zg7RhmttDMlpvZ8j17eQeiGr1jxac18eanoo4xaPX3Ldcn/3WR2jIdA267KpXSvB+m5ZlMGZJhmAq9i+nh6kcPA6Eq+WvibMfBEkUHgNIY6rDhLe7+BknvlnSlmb29n22tl7Vez9Dk7je6+3x3nz/hOPYqqzbrujs14X83KNfVFXWUwXPX1P+zTG/99TXqzPX/ztpn1lwqW7muTMEwTIXexfRwdaOHgdCV/DVxsnFUGDkBoGyGNGxw953B992S7lLPLmC7zGyKJAXfdwebt0rK/7yq6ZJ2DuXxUXk6cl26+OZrlHykMs583h/PZHTS19bqjU98UlnP9bpN1nPK3Dued9NQUnQxikEPA+GjhwHgWIMeNpjZKDMbffiypHdJelbSUkmXBZtdJunu4PJSSQvMrN7MZkuaJ2nZYB8flenNT35Ks7+1qmpO0JV9+WXN+tx+/ZcX5vd6/W86mzTl921lToXhhC5GsehhIFz0MAD0rmYIt50k6S4zO3w/v3D3+8zsSUlLzOzTkrZL+pAkufsaM1siaa2kjKQr3QfxuVWoWLe2H6cZX80q29kZdZRQZZ7fqUe/e5b2f/vPak6OPOK6a5Z+THNbHo8oGYYJuhgFo4eBkqCHAaAXgx42uPtmSX/Vy/peSef0cZtvSvrmYB8TlSvtWf3Pmz+s6WsejTpKSYy9+xm9+9KP609/dZvqrVZSzwnJTvhFe+8HYQIhoYtRKHoYKA16GAB6V4qPvgSOcWv7FM26ZVvUMUom19mpcZfu05u+dbW+vOv1Snlal676pPT0+qijAYAkehgAAJQXwwaUxT//8SJlWp+POkZJZffu06TrHtWqcyfojd+7WtO/4ZyQDEBs0MMAAKCchnLOBqAgKU9r8iPDZ66VfXGvpnz/UfV+XnQAKD96GAAAlNvweeWBSGQ9px/uP0nNy/dEHQUAhiV6GAAARIFhA0pme6ZDJ/zus3ro/a9XdsOmqOMAwLBDDwMAgKhwGAVK4o+HEvriP/+TTvrFk8pwvCwAlB09DAAAosSwAaG74cA0Lfnc+Rr74GN83BgARIAeBgAAUWPYgFAtapuvNVe/VrV/XhF1FAAYluhhAAAQBwwbEJrv7JurDZfOUWLtqqijAMCwRA8DAIC44ASRCMV9nfX6/T/8tbJrN0QdBTFTM32aamZMjzoGUPXoYfSFHgYARIFhA4Zse6ZDX/vvn1Ly4ZVRR0EMrV80Qy0LZ0QdA6hq9DD6Qw8DAKLAsAFD8mL2oN61+J807tYno46CGEqeOFffv/AWvfuCJ5Vsaoo6DlCV6GH0hx4GAESFYQMGLeVpnfnvn9esb66U87Fq6MWmSyfqPSM79OWJf9KhN58YdRyg6tDDGAg9DACICsMGDNrFGy7Ua/7rc/JUKuooiKGaaVN1xQfuVdISmpgcpa1/a1FHAqoOPYz+0MMAgCgxbMCgbEl36NC/TFX2wEtRR0FM7bx4lq4cu+mVn7/110tUczzHDANhoYcxEHoYABAlhg0YlIue+ozq//BM1DEQU8nmZr3tE0+q1pKvrF08ap92vocXuUBY6GH0hx4GAESNYQOKdsvL4zXlfyTl6e6ooyCmXvjwa/Qvk/9yxFqtJTX+gztktXURpQKqBz2MgdDDAICoMWxAUdoyHbruux+SL3826iiIqeTYMTr70ys1MnHsi9nvzLlDduoJEaQCqgc9jIHQwwCAOGDYgIJlPaezH1ik8TeviDoKYuzAeSfrv09+uNfrTquv1+YPjSlzIqB60MMoBD0MAIgDhg0o2OWtb9MpX21lt130KTFypMZevl3NyZF9brPgvY8o2dxcxlRA9aCHMRB6GAAQFwwbUJBnurvU8vVTlGl7IeooiLGO816nW074Vb/bLBr3pNrfwWe9A8Wih1EIehgAEBcMG1CQC+/7nOp+vzzqGIgzM+1ecEjjk6P63aw5OVJtH0pJxue9A8WghzEgehgAECMMGzCgR7qk11zfLrlHHQUxlpw7S98/Y0lB2/6PM+5WzczpJU4EVA96GIWghwEAccKwAQO67P6Fyj3zXNQxEHPb/m6Kzm/oLGjbDza+qNaL+ax3oFD0MApBDwMA4oRhA/q1IX1Q825J8W4a+pVsatK5H1ympBVWKbWW1Mkfek6JESNKnAyofPQwCkEPAwDihmED+nVlywIlnlwXdQzEXOdbTtKXJ/6pqNt8e8ZSdZ99aokSAdWDHkYh6GEAQNwwbECfsp7T7t/M4CPWMKBdb6wt+jYzaxq1eQEVBPSHHkah6GEAQNzwXxj0aUV3VtPufzHqGKgAs767Shd//vP6QMu5Snm64Nv9r7++QzUzOEEZ0Bd6GIWihwEAccOwAX36t13vlG/eHnUMVIBcZ6caf/WEut7frdfd8jntzh5U2rNqy3Roe6ZDW9Idast0HPMC+IONL2rnhTMjSg3EHz2MQtHDAIC4qYk6AOLrDytP0Yldy6KOgQqSPfCS5vzXJ3Xh6i8oU29qXt+pRCoj5SSvTejlE0Zp15ukz/6nh/QPzWs0MlGn0e9rk91QI89koo4PxA49jGLRwwCAuGDYgF6lPK1Jj7LjC4rnmYyafvH4qz/nXde0TGr6hfSnCbN0+4fO1bs+86hybrK6Ol7kAkehhzFY9DAAIA4G/SrGzE4ys1V5Xy+b2SIz+4aZPZ+3fkHebb5kZi1mtt7MzgvnV0ApbE6n1fzMgahjoEpl9+zRxB89qtXnT1TTRw8o11nY58LjWHRx9aKHUUr0cHjoYQDo3aD3bHD39ZJOkyQzS0p6XtJdkj4p6Qfu/t387c3sFEkLJJ0qaaqkB83sRHfPDjYDSuf3B0+RbW+LOgaqXHbX7qgjVDy6uHrRwygHenjo6GEA6F1Y+2eeI2mTu2/rZ5sLJd3u7il33yKpRdKZIT0+Qrb8pVnKHTwUdQwAxaGLqwg9DFQkehgAAmENGxZIui3v56vM7Bkzu8nMmoO1aZJ25G3TGqwhhp7ZPYXPdQcqD11cRehhoCLRwwAQGPKwwczqJL1f0q+CpeslzVXP7mRtkr53eNNebu69rMnMFprZcjNbvmcve5RF4eDmMVFHAFCEsLuYHo4ePQxUllK/Js52HAw3MACUWBh7Nrxb0kp33yVJ7r7L3bPunpP0Y726W1irpBl5t5suaWdvd+juN7r7fHefP+G4ZAgRUazRWzgDOlBhQu1iejh69DBQcUr6mjjZOKqE0QEgfGG8krlEebuLmdmUvOsulvRscHmppAVmVm9msyXNk8SHh8dQytMauTsXdQwAxaGLqwg9DFQkehgA8gz60ygkycxGSjpX0mfzlr9tZqepZ3ewrYevc/c1ZrZE0lpJGUlXctbdeOrMpTVyN8cJA5WCLq4+9DBQWehhADjWkIYN7t4p6bij1i7tZ/tvSvrmUB4Tpbcjm1DdC+3iv3pAZaCLqw89DFQWehgAjsUBoThGZ65W1sU7agAQFXoYAABUOoYNOMaTh+bI9x+IOgYADFv0MAAAqHRDOowClW1/tlM/f/lk3bLlTXrxhSZZZ1I+Iqe63TWafXBF1PEAoOrRwwAAoFoxbBiGNqQP6oMrP6OxPxut0X/aqHH7NmqcH/nxzr1+2DMAIBT0MAAAqHYMG6pI2rP6S1etfrbnbD3z4lS91N6guvqMGkekdObEbXp4xzzpL2M15bFOTXviWXkmw8nHACBE9DAAAEAPhg0VIOs5rUl36+a9Z+vBHSepK1WrbDahESPSeuPU7frw+GW6ZdfZWnfbyZpyX5tyO3ZqXGqDxh2+AzNtqKvT1O51UvDOGe+YAUDh6GEAAIDiMGyIsbRn9cMDc3XDne/W7DsPSBu3aXLnuiO2aauv1/+d/C7l9h3QxPZHe3+HzF2eSpUjMgBUFXoYAABgcBg2xNAz3V1auPZjSt89QZN/u03HP/+Ycn1s66mUMtt2lDUfAFQ7ehgAAGBoGDaUWdZz2p3t1K5srXZnG9WU6NLoRLeOrzHtyWb0N4/8g074YVZjl6+VZ1qUiTowAFQZehgAAKD0GDbkeTF7UJszdWrPjVDWE2pKdGlObZcmJkcN6X6zntOGdJc+89zH1HHPZE1Y2ana3e2yQyl5fa1UV6uOeWOV6HbNe+ApeSbDsbwAhiV6GAAAoDoM+2HD9kyHPr/tIq2790RNWpZSw+a9Ujp4H6u2RqmZ47Rr/gh1nNytt5zcomum/F6n1ddLkjpz3VqXlg7kGrTy0Czd0/ZaHehsUMfBEWoc1aX2jgaNWN2gEXtd49Z0qvHJdRqV3ixJxxzT27C25zsvbgEMN/QwAABA9Rl2w4bOXLeWdEzXD1veoQPPHqfpD6c14k/PanrXo5J0zO6yyc1bNfWPPZdfrK/Xta/9lDb9XZOyjTlN/YPU9Oxe2cFD8pfbVdfRqom5rCYGt5141H3xAhYA6GEAAIDhYNgMG/ZnO3Xl9vdow/97jSb9drPG7dqocb5Bkvo86dfRPJWSVqzRnBWvrvH56ABQGHoYAABg+KjKYUPas0p7VstSI3T16g8r80SzJjydVsN/PKfj2h/jZF8AUGL0MAAAwPBW0cOGzly37uiYqp9se5t2bB0v60oq2WUau14atSujkS37NXnjBinX875Xoe+cAQAKQw8DAACgN7EfNmxLj9JXd79Of9kzRx2peh3qrlUmk1B65yiNX2kaf/9mNezerhNzW465LbvWAsDQ0cMAAAAoVuyHDannpBVnj1bdoe0a58ee2otdcQGgtOhhAAAAFCv2wwa5K9fZGXUKABi+6GEAAAAUKRF1AAAAAAAAUF0YNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCxbABAAAAAACEimEDAAAAAAAI1YDDBjO7ycx2m9mzeWvjzOwBM9sYfG/Ou+5LZtZiZuvN7Ly89TPMbHVw3XVmZuH/OgBQnehiAIgWPQwAxSlkz4bFks4/au1aSQ+5+zxJDwU/y8xOkbRA0qnBbX5kZsngNtdLWihpXvB19H0CAPq2WHQxAERpsehhACjYgMMGd39E0r6jli+UdHNw+WZJF+Wt3+7uKXffIqlF0plmNkVSk7s/5u4u6Za82wAABkAXA0C06GEAKM5gz9kwyd3bJCn4PjFYnyZpR952rcHatODy0eu9MrOFZrbczJanlRpkRACoeiXrYnoYAApSttfE2Y6DoQYHgFIL+wSRvR1z5v2s98rdb3T3+e4+v1b1oYUDgGFiyF1MDwPAkIT+mjjZOCq0cABQDoMdNuwKdgNT8H13sN4qaUbedtMl7QzWp/eyDgAYPLoYAKJFDwNAHwY7bFgq6bLg8mWS7s5bX2Bm9WY2Wz0nvVkW7FbWbmZnBWfc/XjebQAAg0MXA0C06GEA6EPNQBuY2W2S3iFpvJm1Svq6pH+RtMTMPi1pu6QPSZK7rzGzJZLWSspIutLds8FdXaGes/g2SLo3+AIAFIAuBoBo0cMAUBzrORFufDXZOH+TnRN1DAA4woN+xwp3nx91jnKghwHE0XDqYUmqnznDp35hUdQxAOAIWxdd02cXh32CSAAAAAAAMMwxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCxbABAAAAAACEimEDAAAAAAAIFcMGAAAAAAAQKoYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoWLYAAAAAAAAQsWwAQAAAAAAhIphAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUAw4bzOwmM9ttZs/mrX3HzJ4zs2fM7C4zGxuszzKzQ2a2Kvi6Ie82Z5jZajNrMbPrzMxK8hsBQBWiiwEgWvQwABSnkD0bFks6/6i1ByS91t1fL2mDpC/lXbfJ3U8Lvi7PW79e0kJJ84Kvo+8TANC3xaKLASBKi0UPA0DBBhw2uPsjkvYdtXa/u2eCHx+XNL2/+zCzKZKa3P0xd3dJt0i6aFCJAWAYoosBIFr0MAAUJ4xzNnxK0r15P882s6fM7E9m9rZgbZqk1rxtWoM1AEA46GIAiBY9DAB5aoZyYzP7iqSMpFuDpTZJM919r5mdIenfzexUSb0di+b93O9C9exephEaOZSIAFD1StHF9DAAFK4cr4mTzc3hhgaAEhv0sMHMLpP0XknnBLuByd1TklLB5RVmtknSieqZ2ubvVjZd0s6+7tvdb5R0oyQ12bg+CxgAhrtSdTE9DACFKddr4vqZM+hiABVlUIdRmNn5kr4o6f3u3pm3PsHMksHlOeo56c1md2+T1G5mZwVn3P24pLuHnB4AhjG6GACiRQ8DQN8G3LPBzG6T9A5J482sVdLX1XOm3XpJDwSf1vN4cJbdt0v6b2aWkZSVdLm7Hz6RzhXqOYtvg3qOZ8s/pg0A0A+6GACiRQ8DQHEGHDa4+yW9LP+0j23vlHRnH9ctl/TaotIBACTRxQAQNXoYAIoTxqdRAAAAAAAAvIJhAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCxbABAAAAAACEimEDAAAAAAAIFcMGAAAAAAAQKoYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoWLYAAAAAAAAQsWwAQAAAAAAhIphAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAEI14LDBzG4ys91m9mze2jfM7HkzWxV8XZB33ZfMrMXM1pvZeXnrZ5jZ6uC668zMwv91AKA60cUAEC16GACKU8ieDYslnd/L+g/c/bTg6x5JMrNTJC2QdGpwmx+ZWTLY/npJCyXNC756u08AQO8Wiy4GgCgtFj0MAAUbcNjg7o9I2lfg/V0o6XZ3T7n7Fkktks40symSmtz9MXd3SbdIumiQmQFg2KGLASBa9DAAFGco52y4ysyeCXYpaw7WpknakbdNa7A2Lbh89HqvzGyhmS03s+VppYYQEQCqXkm6mB4GgIKV5TVxtuNg2LkBoKQGO2y4XtJcSadJapP0vWC9t2POvJ/1Xrn7je4+393n16p+kBEBoOqVrIvpYQAoSNleEycbRw0xKgCU16CGDe6+y92z7p6T9GNJZwZXtUqakbfpdEk7g/XpvawDAAaJLgaAaNHDANC3QQ0bguPNDrtY0uGz8i6VtMDM6s1stnpOerPM3dsktZvZWcEZdz8u6e4h5AaAYY8uBoBo0cMA0LeagTYws9skvUPSeDNrlfR1Se8ws9PUs9vXVkmflSR3X2NmSyStlZSRdKW7Z4O7ukI9Z/FtkHRv8AUAKABdDADRoocBoDjWcyLc+Gqycf4mOyfqGABwhAf9jhXuPj/qHOVADwOIo+HUw5JUP3OGT/3CoqhjAMARti66ps8uHsqnUQAAAAAAAByDYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCxbABAAAAAACEimEDAAAAAAAIFcMGAAAAAAAQKoYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoWLYAAAAAAAAQsWwAQAAAAAAhIphAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChGnDYYGY3mdluM3s2b+2XZrYq+NpqZquC9VlmdijvuhvybnOGma02sxYzu87MrCS/EQBUIboYAKJFDwNAcWoK2GaxpH+VdMvhBXf/8OHLZvY9SS/lbb/J3U/r5X6ul7RQ0uOS7pF0vqR7i04MAMPTYtHFABClxaKHAaBgA+7Z4O6PSNrX23XBJPbvJN3W332Y2RRJTe7+mLu7ekr6oqLTAsAwRRcDQLToYQAozlDP2fA2SbvcfWPe2mwze8rM/mRmbwvWpklqzdumNVjrlZktNLPlZrY8rdQQIwJA1Qu9i+lhAChKyV8TZzsOhp8aAEqokMMo+nOJjpzgtkma6e57zewMSf9uZqdK6u1YNO/rTt39Rkk3SlKTjetzOwCApBJ0MT0MAEUp+Wvi+pkz6GIAFWXQwwYzq5H0AUlnHF5z95TU8xaYu68ws02STlTP1HZ63s2nS9o52McGAPSgiwEgWvQwAPRuKIdR/I2k59z9lV3BzGyCmSWDy3MkzZO02d3bJLWb2VnBMW0fl3T3EB4bANCDLgaAaNHDANCLQj768jZJj0k6ycxazezTwVULdOxJcN4u6Rkze1rSHZIud/fDJ9K5QtJPJLVI2iTOugsABaOLASBa9DAAFGfAwyjc/ZI+1j/Ry9qdku7sY/vlkl5bZD4AgOhiAIgaPQwAxRnqp1EAAAAAAAAcgWEDAAAAAAAIFcMGAAAAAAAQKoYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoWLYAAAAAAAAQsWwAQAAAAAAhIphAwAAAAAACBXDBgAAAAAAECqGDQAAAAAAIFQMGwAAAAAAQKgYNgAAAAAAgFAxbAAAAAAAAKFi2AAAAAAAAELFsAEAAAAAAISKYQMAAAAAAAgVwwYAAAAAABAqhg0AAAAAACBUDBsAAAAAAECoGDYAAAAAAIBQMWwAAAAAAAChYtgAAAAAAABCxbABAAAAAACEimEDAAAAAAAIFcMGAAAAAAAQKoYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoRpw2GBmM8zsYTNbZ2ZrzOzqYH2cmT1gZhuD7815t/mSmbWY2XozOy9v/QwzWx1cd52ZWWl+LQCoHvQwAESPLgaA4hSyZ0NG0hfc/WRJZ0m60sxOkXStpIfcfZ6kh4KfFVy3QNKpks6X9CMzSwb3db2khZLmBV/nh/i7AEC1oocBIHp0MQAUYcBhg7u3ufvK4HK7pHWSpkm6UNLNwWY3S7oouHyhpNvdPeXuWyS1SDrTzKZIanL3x9zdJd2SdxsAQB/oYQCIHl0MAMUp6pwNZjZL0umSnpA0yd3bpJ7ylTQx2GyapB15N2sN1qYFl49eBwAUiB4GgOjRxQAwsIKHDWbWKOlOSYvc/eX+Nu1lzftZ7+2xFprZcjNbnlaq0IgAUNXoYQCIXlRdnO04WHxYAIhQQcMGM6tVT6ne6u6/DpZ3BbuBKfi+O1hvlTQj7+bTJe0M1qf3sn4Md7/R3ee7+/xa1Rf6uwBA1aKHASB6UXZxsnFUeL8IAJRBIZ9GYZJ+Kmmdu38/76qlki4LLl8m6e689QVmVm9ms9Vz0ptlwW5l7WZ2VnCfH8+7DQCgD/QwAESPLgaA4ljPeWn62cDsrZL+Q9JqSblg+cvqOUZtiaSZkrZL+pC77wtu8xVJn1LPWXsXufu9wfp8SYslNUi6V9I/+AABzKxd0vpB/G7lMl7Si1GHGEDcM5Jv6OKesRrzHe/uE0oR5mj0cEGq8d9YOcU9nxT/jHHPJ8U/Y7H5ytbDEl1cgLj/+5Lin5F8Qxf3jNWYr88uHnDYEDUzW+7u86PO0Ze455Pin5F8Qxf3jOSrbJXw/MQ9I/mGLu4Z455Pin/GuOeLWtyfn7jnk+KfkXxDF/eMwy1fUZ9GAQAAAAAAMBCGDQAAAAAAIFSVMGy4MeoAA4h7Pin+Gck3dHHPSL7KVgnPT9wzkm/o4p4x7vmk+GeMe76oxf35iXs+Kf4ZyTd0cc84rPLF/pwNAAAAAACgslTCng0AAAAAAKCCxHbYYGbnm9l6M2sxs2sjyjDDzB42s3VmtsbMrg7Wv2Fmz5vZquDrgrzbfCnIvN7MzitTzq1mtjrIsjxYG2dmD5jZxuB7cxQZzeykvOdplZm9bGaLon4OzewmM9ttZs/mrRX9nJnZGcFz32Jm1wWfl12qfN8xs+fM7Bkzu8vMxgbrs8zsUN5zeUNE+Yr+m5YqXz8Zf5mXb6uZrQrWy/4cVgq6uOCMse3h4PFi18Vx7+F+MtLFQ8tHDxeJHi4qZ2y7OI49HDxGrLuYHi5ZxvJ0sbvH7ktSUtImSXMk1Ul6WtIpEeSYIukNweXRkjZIOkXSNyRd08v2pwRZ6yXNDn6HZBlybpU0/qi1b0u6Nrh8raRvRZkx7+/6gqTjo34OJb1d0hskPTuU50zSMklvlmTq+Zzsd5cw37sk1QSXv5WXb1b+dkfdTznzFf03LVW+vjIedf33JH0tquewEr5EFxeTcasqoIfz/q6Rd3EfPRKbHu4nI108hHxHXU8PD/wc0sPF5dyqCuhixaSHg8eJdRf3kY8eHmLGo64vWRfHdc+GMyW1uPtmd++WdLukC8sdwt3b3H1lcLld0jpJ0/q5yYWSbnf3lLtvkdSint8lChdKujm4fLOki/LWo8p4jqRN7r6tn23Kks/dH5G0r5fHLvg5M7Mpkprc/THv+V/gLXm3CT2fu9/v7pngx8clTe/vPsqdrx9lf/4GyhhMYv9O0m393UepM1YAunho4tjDUky6OO493FdGujicfPRwwejhoYtjF8eih6X4dzE9XNqMpe7iuA4bpknakfdzq/ovtJIzs1mSTpf0RLB0VbDrzk15uxZFldsl3W9mK8xsYbA2yd3bpJ7/QEiaGHFGSVqgI/8hx+k5lIp/zqYFl49eL4dPqWeieNhsM3vKzP5kZm8L1qLIV8zfNMrn722Sdrn7xry1uDyHcUIXF65SeliKdxdXUg9LdPFQ0MOFibovjhHjHpYqp4vj3MNSZXUxPTw0Je3iuA4bejv+w8ueImBmjZLulLTI3V+WdL2kuZJOk9Smnl1PpOhyv8Xd3yDp3ZKuNLO397NtJBnNrE7S+yX9KliK23PYn74yRfVcfkVSRtKtwVKbpJnufrqkz0v6hZk1RZCv2L9plH/rS3Tkf+Tj8hzGTax+/5h3cex7WKroLo5dj9DFQ0YPFyZWv3/Me1iqgC6u4B6WYtYj9HAoStrFcR02tEqakffzdEk7owhiZrXqKdVb3f3XkuTuu9w96+45ST/Wq7s0RZLb3XcG33dLuivIsyvY3eXwbi+7o8yontJf6e67gqyxeg4DxT5nrTpyt62SZzWzyyS9V9JHg12YFOyKtTe4vEI9x3+dWO58g/iblv35kyQzq5H0AUm/PLwWl+cwhujiAlVID0vx7+LY93CQjS4eAnq4KPRwESqki+Pew1IFdDE9PHTl6OK4DhuelDTPzGYH078FkpaWO0RwDMtPJa1z9+/nrU/J2+xiSYfP7LlU0gIzqzez2ZLmqedEGqXMOMrMRh++rJ4TpjwbZLks2OwySXdHlTFwxNQsTs9hnqKes2C3snYzOyv4t/LxvNuEzszOl/RFSe9398689QlmlgwuzwnybY4gX1F/03Lny/M3kp5z91d2BYvLcxhDdHFh+Sqlh6X4d3Gse1iii0NCDxeOHi48Y6V0cdx7+PBjx7aL6eHQlL6LPeQzmob1JekC9ZzpdpOkr0SU4a3q2T3kGUmrgq8LJP1M0upgfamkKXm3+UqQeb3KcLZk9Zyd+Onga83h50rScZIekrQx+D4uwowjJe2VNCZvLdLnUD0l3yYprZ5J3acH85xJmq+eAtkk6V8lWQnztajnOK/D/xZvCLb9YPC3f1rSSknviyhf0X/TUuXrK2OwvljS5UdtW/bnsFK+RBcXki/2PRw8Zqy6uI8eiU0P95ORLh5CvmB9sejhYp5HeriwjLHvYsWsh4PHiHUX95GPHh5ixmB9sUrcxRbcEAAAAAAAIBRxPYwCAAAAAABUKIYNAAAAAAAgVAwbAAAAAABAqBg2AAAAAACAUDFsAAAAAAAAoWLYAAAAAAAAQsWwAQAAAAAAhIphAwAAAAAACNX/B2/IqwCXwhs4AAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 1296x432 with 3 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, axes = plt.subplots(ncols=3, figsize=(18,6))\n",
"\n",
"axes[0].imshow(eroded_footprint)\n",
"axes[1].imshow(eroded_isotropic)\n",
"axes[2].imshow(eroded_footprint.astype(int)-eroded_isotropic.astype(int), vmin=-1, vmax=1)\n",
"\n",
"print(\"IoU:\", iou(eroded_footprint, eroded_isotropic))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.10.4 (conda)",
"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.10.4"
},
"orig_nbformat": 4,
"vscode": {
"interpreter": {
"hash": "e6c6e9b29087cbfec9dabff591da7a30e404354c2c5afe71b287f8ce14117625"
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment