{ "cells": [ { "cell_type": "markdown", "metadata": { "id": "1kl60n4rCSsz" }, "source": [ "# Create water occurrence image from yearly satellite images\n", "\n", "*Written by Men Vuthy, 2022*\n", "\n", "---\n", "\n", "\n", "\n" ] }, { "cell_type": "markdown", "metadata": { "id": "NMzDLF2lGXYA" }, "source": [ "### Objective\n", "\n", "\n", "\n", "* The objective is to create a water occurrence image from Landsat and Sentinel satellite images. The data contains multiple images of dry season in different years.\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Data folder:\n", "\n", "![img](img-markdown/dataset.png \"Title\")" ] }, { "cell_type": "markdown", "metadata": { "id": "MvbiGrb4CfMP" }, "source": [ "### Installation" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "l-w6XsesCdiV", "outputId": "af7776ee-655c-46d0-ba54-21e0e1fef2f6" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n", "Requirement already satisfied: geopandas in /usr/local/lib/python3.7/dist-packages (0.10.2)\n", "Requirement already satisfied: pandas>=0.25.0 in /usr/local/lib/python3.7/dist-packages (from geopandas) (1.3.5)\n", "Requirement already satisfied: fiona>=1.8 in /usr/local/lib/python3.7/dist-packages (from geopandas) (1.8.21)\n", "Requirement already satisfied: pyproj>=2.2.0 in /usr/local/lib/python3.7/dist-packages (from geopandas) (3.2.1)\n", "Requirement already satisfied: shapely>=1.6 in /usr/local/lib/python3.7/dist-packages (from geopandas) (1.8.4)\n", "Requirement already satisfied: munch in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (2.5.0)\n", "Requirement already satisfied: click-plugins>=1.0 in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (1.1.1)\n", "Requirement already satisfied: certifi in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (2022.6.15)\n", "Requirement already satisfied: six>=1.7 in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (1.15.0)\n", "Requirement already satisfied: setuptools in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (57.4.0)\n", "Requirement already satisfied: click>=4.0 in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (7.1.2)\n", "Requirement already satisfied: attrs>=17 in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (22.1.0)\n", "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.7/dist-packages (from fiona>=1.8->geopandas) (0.7.2)\n", "Requirement already satisfied: pytz>=2017.3 in /usr/local/lib/python3.7/dist-packages (from pandas>=0.25.0->geopandas) (2022.2.1)\n", "Requirement already satisfied: numpy>=1.17.3 in /usr/local/lib/python3.7/dist-packages (from pandas>=0.25.0->geopandas) (1.21.6)\n", "Requirement already satisfied: python-dateutil>=2.7.3 in /usr/local/lib/python3.7/dist-packages (from pandas>=0.25.0->geopandas) (2.8.2)\n", "Looking in indexes: https://pypi.org/simple, https://us-python.pkg.dev/colab-wheels/public/simple/\n", "Requirement already satisfied: rasterio in /usr/local/lib/python3.7/dist-packages (1.2.10)\n", "Requirement already satisfied: numpy in /usr/local/lib/python3.7/dist-packages (from rasterio) (1.21.6)\n", "Requirement already satisfied: click-plugins in /usr/local/lib/python3.7/dist-packages (from rasterio) (1.1.1)\n", "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.7/dist-packages (from rasterio) (0.7.2)\n", "Requirement already satisfied: affine in /usr/local/lib/python3.7/dist-packages (from rasterio) (2.3.1)\n", "Requirement already satisfied: attrs in /usr/local/lib/python3.7/dist-packages (from rasterio) (22.1.0)\n", "Requirement already satisfied: setuptools in /usr/local/lib/python3.7/dist-packages (from rasterio) (57.4.0)\n", "Requirement already satisfied: snuggs>=1.4.1 in /usr/local/lib/python3.7/dist-packages (from rasterio) (1.4.7)\n", "Requirement already satisfied: click>=4.0 in /usr/local/lib/python3.7/dist-packages (from rasterio) (7.1.2)\n", "Requirement already satisfied: certifi in /usr/local/lib/python3.7/dist-packages (from rasterio) (2022.6.15)\n", "Requirement already satisfied: pyparsing>=2.1.6 in /usr/local/lib/python3.7/dist-packages (from snuggs>=1.4.1->rasterio) (3.0.9)\n" ] } ], "source": [ "!pip install geopandas\n", "!pip install rasterio" ] }, { "cell_type": "markdown", "metadata": { "id": "6CHUB7Bq31u6" }, "source": [ "### Code" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "lR57ejgvCPLm", "outputId": "240eb667-8c46-4526-af17-98d5961e41f9" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/content/drive/MyDrive/Colab Notebooks/cheata\n" ] } ], "source": [ "cd /content/drive/MyDrive/Colab Notebooks/cheata" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "id": "doUBsXcECllF" }, "outputs": [], "source": [ "# Import all necessary modules\n", "import os\n", "import glob\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "\n", "import rasterio\n", "import geopandas as gpd" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "5BT6qnHxC5GG", "outputId": "ecd75dcc-c267-4cad-a1b2-ae582c8bf741" }, "outputs": [ { "data": { "text/plain": [ "['tiff/201003051.tif',\n", " 'tiff/201102041.tif',\n", " 'tiff/201202071.tif',\n", " 'tiff/201302251.tif',\n", " 'tiff/201403241.tif',\n", " 'tiff/201502071.tif',\n", " 'tiff/201604091.tif',\n", " 'tiff/201703151.tif',\n", " 'tiff/201803101.tif',\n", " 'tiff/201902131.tif',\n", " 'tiff/202003091.tif']" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Read data\n", "# File and folder paths\n", "file_path = \"tiff\"\n", "\n", "# Make a search criteria to select the ndvi files\n", "q = os.path.join(file_path, \"*.tif\")\n", "\n", "# sorted files by name\n", "paths = sorted(glob.glob(q)) \n", "\n", "paths" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "Tf97BY_bDYzF", "outputId": "69d88156-98a1-4cca-90cc-6c2a62877919" }, "outputs": [ { "data": { "text/plain": [ "[,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Open all files in folder\n", "image_source = []\n", "\n", "for file in paths:\n", " image = rasterio.open(file)\n", " image_source.append(image)\n", "\n", "image_source" ] }, { "cell_type": "markdown", "metadata": { "id": "PTCSc2MjVAJ1" }, "source": [ "**Create water occurrence for Landsat images**" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "id": "6Xwkw40eDt8G" }, "outputs": [], "source": [ "image_binary_landsat = []\n", "\n", "for i in range(len(image_source)-5):\n", " # read all image into numpy array\n", " img = image_source[i].read()\n", "\n", " # make binary image, 0=no water, 1= water\n", " img_binary = np.where(img != 0, 1, img)\n", "\n", " # Append to list\n", " image_binary_landsat.append(img_binary)\n", "\n", "# Sum array to create water occurrence\n", "landsat = sum(image_binary_landsat)" ] }, { "cell_type": "markdown", "metadata": { "id": "EajDJ5xaF1uO" }, "source": [ "*Visualize result*" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 419 }, "id": "RpOgEaBLDlDN", "outputId": "24fcea44-43a2-436b-a77a-d936f7e7ad6b" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc4AAAGSCAYAAABqqiE8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3xc13Xo+9+eM71gBp2oBHsTuwpVbNmSrWbJvcXdSVzucxzn+SaOffPetW8S35duK3GKdeMiO4llR26xbMmSZatXkpLYCTb0jgEGM5g+s98fe0YASAAESIAckOvLDz4Ezpw5c2YGmHX23muvrbTWCCGEEGJubBf7BIQQQoilRAKnEEIIMQ8SOIUQQoh5kMAphBBCzIMETiGEEGIeJHAKIYQQ8yCBUwghhJgHCZxiSVFKfUEp9eBp247NsO29ZznW65RSXYtwjiGl1D8rpfqUUnGl1H6l1EcX+nGEEBeHBE6x1DwBXKeUsgCUUnWAA9h+2rbVhX0XjVLKPs02J/ArYDlwLRAE/gj4C6XUZxfzfOZwbmdsE0LMnwROsdS8iAmU2wo/vwb4DXD0tG0ntNY9SqmPKqUOK6WiSqmTSqlPACilfMCDQL1SKlb4qldK2ZRSn1dKnVBKDSulfqCUqijcp0UppZVSv6OU6gB+Pc35fRBoBt6ltT6ltc5orR8Cfh/4U6VUWeFYTUqpHymlBguP87XiAZRSH5t0zoeUUjsK27VSavWk/b6tlPrzwvevU0p1KaX+WCnVB3xLKfUlpdT9Sql/U0qNAR9RSgWVUt9QSvUqpbqVUn8+6YLjI0qpp5RSf6OUGlFKnVJK3T7p8SqUUt9SSvUUbv/JpNvuVEq9rJQaVUo9o5Taci5vrhBLgQROsaRordPA88BrC5teCzwJPHXatmJrcwC4EygDPgp8RSm1Q2s9DtwO9Git/YWvHuDTwFuBG4F6YAT4x9NO40ZgA3DrNKf4RuDBwvEn+yHgBq4tBKoHgHagBWgA7gNQSr0L+BLwocI5vxkYnsNLA7AMqMC0dj9e2PYW4H4gBPw78G0gi2mRbwduAX530jGuwVyEVAF/BXxDKaUKt30X8AKbgBrgK4Vz3g58E/gEUAl8HfgvpZRrjuctxJIigVMsRY8zESRfgwmcT5627XEArfXPtdYntPE48HDh9pl8EvgTrXWX1jqFCWLvPK2b80ta63GtdWKa+1cBvadv1FpngaHC7VdjgvIfFY6T1Fo/Vdj1d4G/0lq/WDjn41rr9tlfjlflgS9qrVOTzu1ZrfVPtNZ5TCC+A/iDwuMOYILf5LHgdq31/9Fa54B7gTqgttD9fTvwSa31SKEl/XjhPh8Hvq61fl5rndNa3wukgF1zPG8hlhQZ8xBL0RPApwpdqNVa62NKqX7g3sK2Kwr7UOhq/CKwFnOh6AX2z3Ls5cCPlVL5SdtyQO2knztnuf8QJthMUQi8VYXbmzABKjvN/ZuAE7McfzaDWuvkadsmn+tyTDd370QjEttp+/QVv9Faxwv7+TEt2bDWemSax10OfFgp9elJ25yYiwMhLjnS4hRL0bOYpJuPAU8DaK3HgJ7Cth6t9alCV+EPgb8BarXWIeAXQDFqTLc0UCdwu9Y6NOnLrbXunrTPbEsK/Qq4vTCGOtk7MK2w5wqP0TxDsk4nsGqGY8cxgb9o2Wm3T3dek7d1Fs6hatJzK9Nab5rh8U4/rwqlVGiG27582mvm1Vp/bw7HFWLJkcAplpxCN+Ru4LOYLtqipwrbiuObTsAFDALZQuvzlkn79wOVSqngpG3/AnxZKbUcQClVrZR6yzxO77tAF/CfhWQih1LqVuDvMV28EeAFTHfuXyilfEopt1Lq+sL9/xX4Q6XUTmWsLp4L8DLwPqWUpZS6DTPWOmda615MV/XfKqXKColQq5RSZz1O4b4PAv+klCovPK9i1/j/AT6plLqmcM4+pdSblFKB+ZyfEEuFBE6xVD2OSVB5atK2JwvbngDQWkcx2aw/wCT5vA/4r+LOWusjwPeAk4Vs0Hrg7sI+DyulopgW4jVzPanCuOgbMK2w54Ex4O8w46Z/XdgnB9yFSdDpwATa9xRu+0/gy8B/AFHgJ5huUoDPFO43Cry/cNt8fQhzQXEI85rczzRdyzP4IJABjmCSrv6gcM67MS39rxWOeRz4yDmcmxBLgpKFrIUQQoi5kxanEEIIMQ8lFTiVUrcppY4qpY4rpT5/sc9HCCGEOF3JdNUWJoW3YiaQd2EqxPyW1vrQRT0xIYQQYpJSanFeDRzXWp8sVIe5D1P1RAghhCgZpRQ4G5g6EbursE0IIYQoGUuucpBS6uMU63Aq+07lLr+4JySEEAIAnY6iswl19j2XtlIKnN2YcmNFjYVtU2it7wHuAbB5a7Rr3bsvzNkJIYSYVeroDy72KVwQpdRV+yKwRim1Qpk1Dd/LpMnqQgghRCkomRan1jqrlPo94JeABXxTa33wIp+WEEIIMUXJBE4ArfUvMEW4hRBCiJJUSl21QgghRMmTwCmEEELMgwROIYQQYh4kcAohhBDzIIFTCCGEmAcJnEIIIcQ8SOAUQggh5kECpxBCCDEPEjiFEEKIeZDAKYQQ4rKmlAoppe5XSh1RSh1WSl072/4lVXJPCCGEuAjuBh7SWr+zsMiId7adJXAKIYS4bCmlgsBrgY8AaK3TQHq2+0hXrRBCiMvZCmAQ+JZS6iWl1L8qpXyz3UFanEIIIUreLbfepoeHhs7pvnv37jkIJCdtukdrfU/hezuwA/i01vp5pdTdwOeB/3em40ngFEIIUfKGh4Z4+vnd53Rfj0MltdZXznBzF9CltX6+8PP9mMA5I+mqFUIIUfI0oPW5fc16XK37gE6l1LrCppuBQ7PdR1qcQgghlgTNWaLgufs08O+FjNqTwEdn21kCpxBCiKVhkeKm1vplYKau3DNI4BRCCLEkLFp7c55kjFMIIYSYB2lxCiGEWBLOluhzoUjgFEIIsSQsYnLQvEjgFEIIsTSURtyUwCmEEKL0aUombkrgFEIIsQTMoZjBhSJZtUIIIcQ8SItTCCHEkiDJQUIIIcQ8SFetEEIIsQRJi1MIIUTJK66OUgokcAohhFgCdMmMcUpXrRBCCDEP0uIUQgixJEhXrRBCCDEPJRI3JXAKIYRYIkokckrgFEIIsSSUSnKQBE4hhBAlT5dQrVoJnEIIIZaEEombMh1FCCGEmA9pcQohhFgaSqTJKYFTCCHEkiDJQUIIIcQ8SHKQEEIIMQ8lEjclcAohhCh9pbQ6imTVCnG5sFlQ0QCByot9JkIsadLiFOJS5QmA3Wm+t7sgl6Fh+zZCITeRSJKuVw7CYNtFPUUh5qc0mpwSOIW4VCWioJT53u4Eh5v+rgHKy5dz41VNxLfU8fizpwjvfgoyqYt7rkKcjVQOEkJcEMVPmkwKMimyR1/kwFAfmcyVXL2plre9cT0P+z10HmiFkR5Ixi7u+QoxixKJmxI4hbisKAXJKCPhOLeuraTW42ZlhYtXttXz0K+PMvby05DLXuyzFGJa0uIUQlxYbj+VO3bxR+/dys0rqsnk8vTFkrzjinq2147jd9t5wOdm4PknIRW/2GcrxBmkAIIQYmG4/ThaNpJJJGGwHSw7ONyQz5mvZAwcLt7022/jT29bT8jroHM4zt3PtPH07k4aGoJ89Z1buWtdNTsbfTyys5EHfr6PbOcRM04qhJhCAqcQS5llZ9d738y/fWgnubwmPJ4hm8sDMJrMcCISI5PTbKkOsrzSSySe4ba/e4LRkTg3X9fCjVc388Of7OV3Y2n+/eO7CLkcbLguyGtWl3P3T2vpePQh6boVpaM0GpwSOIVYyjwbr+Yf3rGFg71j3Lu3m3KfC7ddYVOK5nInO2pCVPiddI0m+PB39/D8YwehYz9ozX2nNvDlP7qdsVs388vvPcJv5TX/+X9dx1giw46aEF/90A4+a1O0PfKgBE9x0WlKJm5K4BRiyXK4+Og7d5DM5Hj3//4lqcMvgsMFlgPKqsFm4SgLkolGITpkviZnV3Qf5k/udvPX//0mTlx3Fcf3n+TTP/LzD2/f8uoud39wB3/osDj2wn7oP3ERnqQQEyQ5SAhx7iw7NG3kY1c18ae/Okbq6F4znllM6olHAMj0KnB6zkz28VdAJgnt+/izb5Xx2utX4XRavLynnQ9FUnz7I1cS8jlQCr72gR3cv6OBH/7yCOEXHwd/pblv4TGEuDC0JAcJIeaprJp1r7+BocEow33DfOi9u4gls/z4J3sgm57+PsqGe81Wkkcm7aMUwfVbiEVi5DqPMHZ0Pw/09HLnu24gGHTx8u423nb3U3z9t6+iIeQhEs/w/s11XNno59+uqCOdyRGLpenuDDPa0WlaoqXSFBCXthL5NZPAKcRSYNl57yffxqevbUEpODESY1W5ny/8/DAMdc18v3yO5OEXp45Rak3klWdNl26x4EEszG8er+HPP7GLoNfJ08+e5H1feYxrr27hrVtrWFceQGtNXbmXlVUeAi6LsWQjTxyt44WnQuSO7V7c5y9ECZHAKcRSYDnY0RTgsz/Zz9hYivJyD11dEQZ6hk3x9tlMDppKmdZhoZLQZOOtr/DQwdV85oYVeN12nt3bw69+dZAHf/wcNSubCYY8dLcP4va6sewWXp+LdDpn7mx3ztzqFWKBlEiDUwKnEEtCLkPPWJqA18nxY4Mc/sVD0LAOh9Mxv7FGu9ME0nzuzNsSUR755X5aqn3csq6C8WSWEz4n8XiaXC5PeDhGPBIl3n584j5KmZZroMqU7BNisUitWiEucd6gSaBZqOLpuSyHuseoCbqJDEWgohFfwEculyPjr4Bw99yOc5bzybft5+v3xFh75UY+evMKVtX4cdoVlT4Hu9siPPm8JlddAQoiLz87cbxiMXkhFomZjlIakVMCpxALyRcyVXu8QYiFITN4/se0WdB8BXdsruaBfYOmoo/lIFUWJJvJgjdkWn2jvSaD1lNmAul0rcqzyaah7zitvx7kC8/sIbhiJWvXL+P3b1rJ2++o48mNlTx0eJi3bq7hf97np+ORn5tmQKk0BcSlbRF/zZRSbUAUyAFZrfWVM+27aIFTKdUEfAeoxTzde7TWdyulKoDvAy1AG/BurfWIUkoBdwN3AHHgI1rrvYt1fkKct9PH9ZwerPrVuL1uxjtOmnmT53v8mhW86R3X8eGrGvmbXx3nld2nwOUDfyXZRALiYziqlmELBkiXVaDzGpfXRcrlM2ttZtPzD6A2y9xnqIPIUAcvvggf/GGQVa9/Hd/92DWsryjjnt2dfOKOtfzJSwfOXNNTKTPdJTp8fs9fiNNcgMuz12utz/qHa1vEE8gC/11rvRHYBXxKKbUR+DzwqNZ6DfBo4WeA24E1ha+PA/+8iOcmxPk7PRnG6cHutDN+/ADEhqGsxrQA58vtB5cXa8UWbnnrLn776ib+n/sP8ML3f0rq2EsmICejkE6CL4TNMn/G3oCXQHkAAFcwiGPlZmwrt0HtKghUzv3xLTtUL4fy+olt8QgnHnyA9/3LswCsrfHywCt9vO2DbzD7T6a1aXUHa00ArWoGl3f+r4MQpyl2bsz3a6EtWotTa90L9Ba+jyqlDgMNwFuA1xV2uxd4DPjjwvbvaK018JxSKqSUqiscR4iLzxMguGkn9U0VxGJplILxWJLh1lYY6sC3aiO+Mh9R20Zy2RxaazI9pyCdmNvxlYKWbbzlLTtwWDbKvA4yuTz/93f30vHK4alFDIY6zP/BWvCvAKCyJoTf72JwIIrNsjE6NGqCqquaXLaCbHtqogU6uRXqC+FfuxmXx0U+l8fpdhLuD5NJ+6eeXz5H2y8f4K6BCH/5yV2Ewwl8bofplj69dZkYA5uFb80VeHweRgZqyB3fUwiqLtNqToxJKT9RSjTwsFJKA1/XWt8z044XZIxTKdUCbAeeB2onBcM+TFcumKDaOeluXYVtEjhFaUjGiLzyLJHu5dh8AWw2G8GqIHWbr8Dv34FSitb9bZDP4fJ7ySQzs3fXOj0mgahwSbzsxtv43qeux24pnu4K872nO2g7OWiScGaa6pFNmT93Bf3dQwwoRXl1OS6XRdSyyOVMgLTsFrYVm0wwj41DLgPjIzA2+Orjl1f4sSxFPJ4hMzIMAyfPfDytie57jh/ubeFzd63jG8904l+9kdjLT029tM+mwXKQz+Upr/ARKPNwqvuYyQBWNqhZDh0Hz+FNEJez80gOqlJKTZ5sfM80gfEGrXW3UqoGeEQpdURr/cR0B1v0wKmU8gM/BP5Aaz2mJmXfaa11IbrP53gfx3TlgsM/+85CLKTi/MfeVvJAXimGB2onlvACsDsJNDSQTqXRiXGzzWYVpoFkprawQsvM8cZHqLv+9bznDav4xYkBfvpMJz2dQ2TSGXKZnLnPTP1N2TRaa4p/V1prhnqH8Pg82B12tNbk82a1FGVTKBSuoOnOTeUyJnDGI8T2Pcvx3haUrww9FjZdzbM85qM/eYbPvPb9fOZ1K/h7BU8f8U9dgszhBruTxNAQoW3NXL++ml+6b+bwLx81wbP3uHk9hJiPc+92HZot2QdAa91d+H9AKfVj4GrgwgdOpZQDEzT/XWv9o8Lm/mIXrFKqDhgobO8GmibdvbGwbYrCVcI9ADZvjaTyiYunfj3eygrixw+A04NjWTOZZIryqjKG+kaobGli1Bdgw9YW/uebN3ByNM6X7nmG5KEXTKAdbAN/BdVXv4b/9uZ1fOvhEwz2jZBJFwKKhuxA5+zJPXYnky9GU5EoDJwkY7MgWEvdxnWMhcfIZrPmeJksuXA/+IKmxVmUy0Lf8bl/Lg2c4os/O8znb1/LlSsqeLqiEfqOTVwYJGOmRR0Lk8ttY1udnw13rOFbQTfPPX4Y2vdJJq6Yt8X6jVFK+QBbYVjRB9wC/OlM+y9mVq0CvgEc1lr/3aSb/gv4MPAXhf9/Omn77yml7gOuASIyvilKllK4gkH8IT9/89VPUOtxU+V1MZ7O4rJsHAqP8ZOX+7lr6zbWl5exfyiC1ppf/q87+ebebdz7jUfA7mLr9Rv5+E0tPLB/iIHesJleAqAhNdgHkf7pH9/uNEk3TjeppJlLaVkWjA1MjKmmxhkoK8fpdpof+7vN8bJpGO6c/rhzZbMIhxN88P97hDvv2Mwff/oWvvWLZgaefnSiwEIyBp4APq+TKo+L0VSGT93Ywl3blvG1n9bR+8TDMsYp5kyzqNdatcCPCxehduA/tNYPzbSz0ot0JkqpG4Angf1AvrD5f2DGOX8ANAPtmOko4UKg/RpwG2Y6yke11rMWwLR5a7Rr3bsX5fyFmFV1C5/41JsYjKY4cmqEbDaP1hq/38l1G2r40LYGfC47J4djfO4H+znyxPPgL6eycRnf/r0baAx5iKdy2C3FV59u43vf/hV4/FgOB3aHnVRfp6nEM9PfpydgqvUMdUzpJkapqUUOnB6oaABvADoPLVxBBocL15rtpI6/Am4///jXH2FVMMC7//pRxvY8bvZRCtfGa7jvC7fSHh3HYVOUOR0ARDNZ7nuhhyd+9qwpEl+cAjPDY6G1lPRbAlJHf0A+PrAo1TA2bd2hv/+LaXtOz2pzY2DP2bpq52Mxs2qfAmZ6AW+eZn8NfGqxzkeIheRvbCISz/DQI4eJt75iknwsB6QTvORw8Y/167A5HOSHCq08rSHczXDnQe46foKb3/4aaoIeTvVH2f3MMRhqh1yWnM0i5/RAanz2y+tEdOqYIkwfWNIJ6Dt+5vbzlU2boJlOQDrBH3zlcX74P29n/aZ6XnjJAp2HZWu467bN2G2KL35zN9l0lo1bGnnftY1UeZz87vXNrKgNcO+9T5jnO9o7fWDPps3UmHhkoii9EBeRVA4S4hw0tVRzvGeM+KEXJwJWsdsxk4L2fa92s0yhNQx38ui37p/oUp0cIItdnKVO6ynTbDKte/j2ni3cvnUZ6fTb6WgPs3JVNRvqfPzevXuInDwBw528uE/x4lPbeOfbr+SW9RXcvLqc9Ptv4Hvff84kJE0XOAsXHfgrzALdY4Ozt1DFpUlTMlXeJXAKcQ5qq32ks/lzywzNZS+dsT2nB+rXQv9JfvTdR9n0ubfw1XduJZXL83h7mL/65nOmaEPx4kJrOPUS9//TcU69/Q62rapieYWbm2/byqPfHzyzFT1ZLGyKQ4SWmcApLdDLTonETQmcQpyLeDLLW3bU8YzDPfcCB5cim4Xd5UQ3riN38mX+7E/u4c+qV2Bzu8mH+2C0b/qWYTrBnmdaicXS3LyjgR0tIXjPG3j0J8/MvjB2JonVtJ5c3ykTPH0hGB9d3OcoSkapJGIvZsk9IS5NDhdKwWubK3Gt2Xqxz+biSsbIDvXhD/pZ9po3QqgOug6SP75n9kLzuSy07+PoE8/zfOsgLeUu3n9VPb/z8Vuhbu3UNUaVMi1NXwjn2p20rG0w80QTY+D0ysoslxF9jv8WmrQ4hZiPYC3XvPUm/vnd2zgwECGbXoQuV29wfmtsXmzhLqLlNXgDXq6+9WoO7qsjm86SOrFvokygzTJdrP5yGOo0zy+fg5Ee9jz4FHWVt/L2rTXcuCJI/Wdv5xcv9fLK7pM4nA6u3rWK69dU4rQrbIUY+S+JDD3PRMz0G3+FafVPLkkoLk0l0uKUwCnEPFzxxtfw9+/YwoPH+vnq/QfIjZ7nCijTSY0v/DEXk9bko6OMDDpZsaKc5Str6e4cJrVstSl0AFBej69pOZZlMQbQOTbR7xbu5pFHD5HO5njXjjrWVnpZ84ZVpF6/AktBLJ1jf1+cruE4vUMxdm2o5Yu/tZkfb6jhkYcPmG7bYK0p5ZcaN0H5cu4+F4tOAqcQlt1Md1C2V1tBr3YxKmUKkidjEFrG/377JroiCe5/votYJLY43YRLMXEoNkw2VMHRI/3ovCbS1QMen1kVJRWHZBSPz0Mw5GWs7eQZg1Wpg8/xcNsRjhy/gXe9fiUBl0UinSerNb94oYvDT+0101VcPvr7d6FvXME7ti/jXduXcXIkyS/29vDS04coa2xkfGyc3LFZp4CLJaiEkmolcIrLnL8CVdvyanUdgHSoGh3uMwXadd5MgfCVU791MyGXk499ZzfHD3WS6zq6cAUFFsuFSp4ZHyU72MNwz0nwV2L5y/AGvKRcm0n3tEEmSU1tGePjaYIrVxHxl0PnwaljoOOjdLy4l79t7aRpTROhkBubTdHZNjgxxzOXpffgYb47HOHYtatYURugvszBG7YuY11TiPbBGPv2djAu01UuSaWSHCSBU1zeMkl0XpPP5V9d19LpdkJdMym3H4basZVV4Pa6+Z071vJo2xBHXzoGHQcu8onP0eQkm8UWLpSWjkewQjupqA6STGboHxlGVTUQDLqIxVI0La8iWO6nw2bBqZemHmOkB0Z66Ow8SGfzZqobqs374XCbwJnPQbibuN3J3pecHPS48HqdrFweYl1dGbm8j2fHxs0Fj7jELE6iz7mQwCkub94QpOJk+nqgrAZnRRX5fJ5sdMxU83F6cDgdaK2xKcWvDw/iLguQtDuXRgm409fJXEwuL7gDMD5CLmtae6GQh9yqFSxvqcBpt6io8GCzKYJBNy3rGmlLxSemn0xuIeay0LGfqO8qKmrKoaLRZNFqDf4K7P4yxsJjBEIBoiNRBnrDxLY2URPywGh/6TRNxMIqkbdVpqOIy4Zt9U5YvsX84PKa1lhizIxtJmMQ6SfdfsQEzZ5CN6zDDcqsZ/m9x9u4eWM1ZeVlpsC6mMpTRmjdRtB5csO9hAcj9PeO0rKikps219Jc7SMaTZPJ5LHZFB6Pg2VrV7Lmjjtxb7zavCeT5bIk21vJZnJsuGqdyZ4FGB8xy6WlU0SOHSZ+ZC/jh/dycF8XL7zcM3XVF3FJ0ef4tdCkxSkuD0qRH4+asm5gkoHcPtPKiRQWcy5WoZlcjWZskExVA3aHndYDbfhuWYndcQG7P+citMyUobvYY3rZNONj468mWEXtToiF2d1Zy8jIBrauryaf1+RyGptN4XRa3Hp9C+PJLH6/E9vWZvY82wrjY2bJtXwOosMMdPSSSVfhW7mB8cN7wLKjbAqGu8z4bc0K8JaRjCeJjcbAVy5ZtZciXTodCRI4xeVBa+htnfh5sA08ZaYVkxibuu/kAGTZsdls2Gw2yqvLWVfhx7JspZMUZLMmxjGLGawXiWv5WpRSZMpqTGH7gVPmhniEE4/20N+znU1bGqkJecjnNalsHp/Twm7ZaKj2s67OT12lj56hGHsfzpiVXwD6TxJ1uahYVkG8aQP+oJ+y8gB9eU2u96SpTpSMYWtYia/MR8rRQDYdnz0pSpKHxHmQwCkuT7msqX0aC8+6m/+Kq/nEu7ez++QIwyMJ9g9EcTotU8mmFMY5C8ky2CxTN/ZiBE6lwGZRXVdJIODiFJAcG5jaPFCK2MHdPD8wRLCulqaWKoaH4xw9PszYaIz1G+so85olxzKZvJmPWXx9dZ5sJMy410V5TTk1tUESiQxur5ts81pS8QREw6RTabyBQndvwxqyp/aVzgWOOG+m27U0mpwSOIWYiVK8/rVrWFPp5snDWeprfPzg2U6GB8dwl4dIZptMYkspKHRrXhTl9biWNeFy2enpCpM8+NyZfWrF4u3dh4l0Hyay15qS+Tq+/B2kM3leeLmboZ4hk7QVXGb2GTgJ/ScYt1lUbVlHb3eY5HjSBMz+k1DZDP5yXG4Xo4OjWHbLZEhbjpkDp7Q2l6bSiJsSOIWYkdbs2d/HmzZVcfLEkCl4ACYxJZcHl+cin+A8+CtMq3RscGGPqxRWZR25bI4Tew+b+ZZzGYiaXGCiajnXbKzlVH+UdDJNvv2AuT1UB+n4xL6xYQZ7w1TWlpPP5UmNx01g7DoITg++HdfhcDlIJ9OkhwdMV3wmuTQLSohplUjclMApLjNKzSvDYLBnCJ/D/Jno4v00ZKJR6DmyGGe4sArdqGSSJkN4EeSGe02X97kk5ASqqFi9inxec6pjlNEjB0wwLBR2r9uymXQqzXBXH3QdJt7WisO5CYfTASO9E8u6ZdOEjxycOG5iDKpXmC71wbaL36UuFoQkBwlxodmdZhrJdMtWWXZTXP207s5AeYBoJmtamAnixnIAACAASURBVAAaUr0dJvmlVP6KZ2M5TGtztG9xxj+Li0yfK4ebdetrOd43RnhwzARflxfKaihbVkMg4KJieYiNG5fx5M8SMNpLNpslWOEnHKgyx6ioh5HCEmaTdR+C+vXmvZXAeUkolTFOmccpLh+5jJmzqab5tc/nph33qqouI5XNk4wnyWVzpHraTHdkqYyROT2zVwfKps8MKKXA7oTmK/A2NHP0cD/tHaNEjuyDyiZsTRsgUM6yBjNvM53OcqptBKusHLRmvK+PSDhmWpV2J2s2rwCHa/rHcblNy3Om24U4BxI4xeVDa+g7Nn3QmzyPc5KOkxNBJ9txxAShUmppphOlE8TnI1BFzfIGbJaNsZExOvYdMa+/y00+FoGuI/R0DmGzKdLpPG63w5Te8wbBcpBKpvA2tmBV1hGNpkxNXrtzatF9uxObzYYrGIDKpov3XMXCKZEKCNJVKy4vswU9ZTtjDDQ5NEiZ047b5yZeuxJiIxPzC8X8FeadeptWkMvlSMaTZNsPmzFYgHAvAO6124iNjNGpNZZl4XQ7Taasrxx6W0mGvVC7koq6KvK5PLUt9fT3n5zokrU7IVCFzWbD4XSQ8gbmPb4tSk+pvHsSOIUommYczLdsGaOpLNl0FpfPQ87lJDvat3BjZm7/tC3ds1LKJPsspQo5dieb3/ZmPnTjcq5tqKQnmuBj//AkEZ2fyHwd7YPyejKpDMvXNTPUP0I+n8ftceJyO3F73QzHwqbrNZMmNhojm8niK/NNfazyenbecg0+j4NXXunC4XSQsTtlXucSpoF8iVz4SFetEDPxhXjPXVv4zhPtpFMmUFpWodDAQjmXoAmm5bQQ43ZKme7PioZFHwf0bbqaez6wkzetq6O6zMUNq6q4+bWrCWy6cmqd2pEecqNDJOJpXB4XNXXlhEIe0qkMHq+LwKp1MNwJfcdID/WTT6dIp9L4N2wHT8AcI9LP3ueO0dsfI9I/RCadubArxYhFUSI9tdLiFGJadidV265iZ6OPHz985NWpKLlsrnSWrFqIdTa1BpuFrWIZebffFLdfJOMHnuMD91RSVeXlylWVfP71q/nIzgb2tw7hcF5LNpNlrLfPdJlH+omEIyilSI4nTXdrMkXCYaeusZJoW73J5h1qB38FvvptJjgWpRPo9gMcy2ahp9UkhpVIa0WcI6lVK0SJs+w4XQ4cNtvE/E0gm82eeytxMfgrIB6ZPkHI7pwxW3iKWJj8iUVcUcSym67YXJYTv3mME6Fl7HmuDK3hozsa+Y+P7yKZyZHLa2KZLM92jfIX9zxFKhLFFQyQy+XIJcxzyGayRCKJiTmp5fXUbLqC4b5hM5908pSbfO7MxbKFWAASOIWYTipOz9FTRNMbUUqZ+ZvDJbACyelmq7XrKzddl33Hz36cyZfyC1mDVyloWA99J8x4bDwC8QhZ4J++0smLd70Oh8PCYdlY31BGpdfBjroAO3at5sXHz1wsXOc1A8dOmTJ8hXN1OO0Eq4KEh3rOfH+UMhm3Nmti+bilNC4sJpGFrIUoTZNXzUiMYdnghl0t2C0bT+/uZKD1hCnlFo9c3PM8G6VMoXSnZ/7ZpE7PwgVOraH7yPRl78ZH2P1MK3osDE4PjwPEIyzbupXfvmMdLz59lFR3mwl4gQpQNpxejynqoLV5Xv4QkfAYwYoyc5GQrZwoYpHPwYrttKxrZLBvhGwmS2posHTqC4t5MclBF/ssDEkOEqLIZkHTponkn3A3f/C/fkImm+ezN6zgup2NuCqrYNlK0/1YynzlUFZtys3Nd2BooS8KZqoVm0mhT+yFZBRXKIgrFMTbtIKN66ppDDn50Ht3mazjkR7oOADt+0j3tkM+a1rFWkN0hNhojO5TPWB3TDxXyw5OD063k3WrKrnh+tVcc+1qHOWVC/vcxAWlz/HfQivxv34hLqB8zoyRBapM1iZA33Eeui/FVStC7N5fKIYwNlz6hcNj4TPXGS1V8Qi5TA7LblHfXE1LjZ8nT0ToCcdxhCrIjAVNQpbDjaqsx+VxkQ2UkY2OmfdssN0k/+RzE5WULAd4ysjn8+x+qYv6hhBut51MKiPzOZewUnnbJHAKMdlgG5TXT/xsd9K8cys2pRjoLkxrcAcu2unNS7E7s1Q+bWaSy5Id6CSbTnC8t4PjD49AJsmuD/8Wu25Yy/59PhLjE+OSWmt8ZT4qV9UxOjJObLSSdGQUxgZxNa8x3bEDJ02w1c3Eo3FstnJisTTByiCRxKq5jfuKkiNjnEKUolx2StFyx+ptvOemlVR67aTbj5hW3NnmA9avM92dF7tGrNNj5mieTxH2C6X4WhW7iW0Wu589hsPpIJ+fOv3H4/NQ11BOT2eYZCJpgmZhGk3q+D6zHFlhYe9sOk12PEZlaDV5renpCmMFQuQG7aXfayBKloxxCnG6SZmZmcFeeiJpAi47rpYNZvrH2T5w+0/Mnu16oSRj5x80Q8tMVuqFls+R7T5BYnCATCpDLpMjn8sTCAXYddVyaqt9jEfHSXUcM2uMFmvUpuJm5ZrqFjPPteswRPrJFla38fg82B12M3Y6k1Ifv75MFZODzuVroclviBCz0Xk6BmMkVpSZWqmBqjOWHjvDpdSSifTP/z4OF6p5E5ZlkT2+d+oUkfl0HcfCZo6pUuSVDbxBVl1zK+lsnuMnw2Q7W00LtaIBXD4zLq2UqYBUPG+lIBnjhedO0LK6jqHeocL2WdoMl9L7d4kpla5aaXGKy9ccSsxZlXVsbgrxxe+8TKLjhMnwnLwCx6VOz7Nci8NFYPM1NKyow7JbhfmT6tVlxKaMH8/nHPI5SIxx5FAvea1xOCwTNB0uVHmtOX5x32zafFl2yrbswnfF1STHohz5xYNkewrzP6X83pJU/HWc79dCk8ApLk9KQc0s00psFt4t1/PlT9/Ikwf6GNz/CoyPmBanp+zCnutSYdmhvJ5oVyddj/+K1Il95uLE5YPqFtZtXwOBwnSQc6mLm8sS2f04Tz9xFKfTgoYNULcGh9MB0aGJ/VJxXGu3g6+cVCKFy+PCGwqaFVjGR8nGYmbx67LqhXne4sI4x6C5GIFTumrF5anYMnH7p6/5WrOCbVe2sLHKTz6vzdzBbBoGTl34c10KQsvMBUUyBun4mauQ9LZy9OdtE631YgtxvquVaE3m5H6OpjLgdIPOk+7vmjr3NJ8jdewVcLjIpDOEO/tg8FThUzQH3YfNeczWXStKjgbyJdJVK4FTXL6G2mcuoZfPkc3mcdoscjltCgqcbWzzclQMQJmkCV6zlbPLpk1LUylIRM/9Mf0V5gLGZp35nhTHUNMJSCfI97eZYH56gC4GUbGklMrMKrnkEpevXHbmv0S7k8oyNz9rHaD1QLsZ21wKHC6obLowmaFOD1QtN/+Pj84eNL1Blt96F//jy/+NZTfedu6P6QkQXLUGKhvNV+Mm83wdLvAGcW28ZuoYdHTYBM1iC1eIBSCBU4ii2lUQNIkmli9Ac5WXX7/UQ35koPSKu09HKaqvudHUdV3szFBPwIwRjg3MabWYumuu56O3ruaxI4MM9Qyddf9plVWDp4xIRydkM1gOB3avF9eyBjNe7fSY/ryymqn3szuxr9x6eSV1XZLOteCelNwTlzNvkNAVOxh98bFFGvF34mpcgdu7AZvNRsdQHJtNmUCUjM5/PO5C05rB3c8sfpD3BqF6uSk6MMdi8L1P/oovvfCM6aI9l/PzhQitv4LR462vFmnPKQU2i6zLZ/ZJjZM6NE3RiWzaTIsplX4+cc5K5S2UwCmWjniE0b1PL9xfj2U36zpmklgrt1FRW4FSiuoaP1pDR+8Ym1dV0nFqkGissCJKqfzlzuT01p/TY8Yf53PeM821dLhwrdmOsimSJw7M70KiOEXkPIweOWCKHRRpbVrWcylKvxR6DMSsJDlIiHO1UMtdFVTu2MUVm+rYu6ed6poAbrfFyrog79teR0OZh8fbh/h+JGYm15d60JzOfIMmmKo7Iz1TA6MvBF5TQSgZDp9fcs+5mC7zWVx2FvNPUCllAbuBbq31nbPtK2Oc4vKVyzL80vO0Hh/mqqtaCIXcaA3tA1EeODrESDzDt3/Tho7HSmPCfLEqznzG6s7lk2awbWrQVAr/2s3YK6pJHXkRelvnf8wLybKbMdjT2NdeZRKJxJKlz/Frjj4DHJ7LjhI4xeUtk6K3vReH3UYmk8PhsEinc5zoi7J/KMLhVwrzNi/2nD/LDs2bcazaem7Vd+bj9GCrNbEjL5NtP2y6PEu95e0Nnvl+OVzs2LV6Yq1VISZRSjUCbwL+dS77S1etuHwpBd4gHp+HvqFxUqkcZWUu3rGrkXKPnUePhtGRodLopvWU4fJ7zeWz5bjwj5/LmozjOWTQXnQ2+5ldu5kUL/zbD2Ssc4nTi/d3+FXgc8Cc1gyUwCkuPy4vrtVbaGypxWZTVFZ6+cObV3H//n5au0a55+etdO07BMOduDdcha2+ifj+Zy5u8MznSCVS0H9q7uexEGtx2iwTbNIJ04W7FMxUmF6C5pJmkoPOWZVSavekn+/RWt8DoJS6ExjQWu9RSr1uLgeTwCkuHw4Xq2+9nXfcsJz2cIKGcjcbq71kteb3v7WHniMnzALHkz5gk0f3mvq0NSvB7YOuQxdv9YxsZu5JMjbLdO+e6xQam2W6PN1+GOo4t2MIsZD0ebU4h7TWV85w2/XAm5VSdwBuoEwp9W9a6w/MdDAJnOLy4Pbzpo+9i5U1fl48NUI4kqC1S/Gs28HJk2H6nv719EEmkzJrWgZruea2q6kJXcnP7nv81bmEi3GeuHzg8ppFs3NZfGs2k8vlSI6Ozf04+dz5tbCWrcbuLyM73GcSks53DqvNMtV7nJ7SWKtULEmL0eejtf4C8AWAQovzD2cLmiCBU1zKlMJavZOyijKuubIZu2Vjb1uYWMxMaUmn88TjCfoOHTl7YIj08/wPHuDKd7yJO99zIw/cGzu3tSqL7M6pwa2QKVu2aSehygDJRAaf30U8nqa/o9+0+qqbTeCZrbTd+bLsEKzF5vWT7e8wwXshpgBVt7By5yZ+6/Ur+dr3XyZy7LCpOnQhur8XostalIRFHOOcFwmc4tJTqNe6ftcW3nVDMxVeO48eCbN7fx/rVlfSEU4wOmxab9H+wbmP3SWi7P737+NYdxXXvPUmwuEEx37xwPy7bpUisGUX6VSaVDhsWmOJMRjpYaxvgHgsbhZptuwm0WVswARajZljOT5igvZifIhUNHLFjVfSfnKAaNi+cEE6Pspd1y2nL5ohOhq9cNmtDpcpwTfceWEeTyxpWuvHgMfOtp8ETnHpcPtZfuPr2bahBq/LTjiW4p6fHUEpxXDfMJmOo3QfqjWBKBU33aHpxPwCkNZk+rt4/oFeGrdvwb56O9mOI3MvCGDZoboFrbVZ+spmTc1WHWon25c+MxhnUtB9xLT+lq02mbWLEQyGOznw4KipqBQ9x5qy04kOc/eff+vVak0rrr2aU4/0LHhBizNkUkunQL+Y1XkmBy0oCZyiNBUDytk+9AKVuJrWkBoJ894Pvo6haIpoIsPJ7gj9/VHC/WHTekuMmQA5uTzbuX5oZ1P4Vm5gqC9MRU0FI2oDmWN7z97y9IXwrNyIw+lgrKMdchloXA89x0zgrGiAsipoe2X6+xdbf72tZlUST2DhK/jkc4tXpWfSVJZTTz5lguhiB044+1hvMXNYlDgtXbVCzMqyY6usIz9b4Kxo4Jb3voHGSh9j8QzhWIr2rgiBgJOOtiEiB1+aWx3T+RofZbzTrAeZ7DxJzZZtDCQ3QXzMtAKnC6B2J2RSJMIjJAZOToyp5vNQvwa6j0JqHNL+sz++1hDuglAdVDZD10GzvVhRaK4fLsFaU7w+lz3/ZKL5KqUEIcsOOi/joEtAqbxFEjhFacqkyB/fM/PtTg/Xv+11tJ4M8+STxwiEAqxbV43f7+Tgyx2kWvcuXmsmnzOZtgWjQ6NsvmoNh19pJ1vsPg1UQlkNvqpKkvEkuVP7TIuxe1JFL6XAZsPlcZGqbDJTYeba2stlweWlcX0LXd2HzCdK3Vpz/7kkLdksVFUDenTQBOHLyektzFJf9Ua8qkTipgROUSLmk/lo2SFQhd/tYDQco3llLU6nxb5XuhiPjJM9te/CdAEWpDuPsf/UQfOYWoMvhKtpDcqmsOwWgfIAoycmfVArZZJjaldCOk3OXkgCCi0zXdNzeR2UApeXvs5CZqpSpnU913FJm4XO5UyiUalcxi8UmwVNmyCXg96jEz0A/grTSo8Nw+g0y4+JkqaBfIn8rkrgFKUhUGU++OcS8CoaCba04LTb8Pg8dJ4aIJfNkRoeNB+KF7oFMbnbsaoZe2UtLo+LsbaTJJ0e0w1YbOEU17JUCkb7YbSPbHGe5GwXD6e3kiqbePvbd9I9NM7uvCZ3fI9Zcmsu3a02C6qaoe/E0iihNx+BSu786Ft457Za0rk8//DwBvY/cxAcThpWNTIyOEq8//jU+8gY55JRGmFTAqcoFTpvVq4Yap89ycbhAn85q9bU8OIrvYwOjWJZFpl0BgZOXvzW02gf2WyasYi/0HrMm9al3VlokebxBgNk0hkyxfHXYqCf7dz9FaZ16PKZpcLC3fzo3ofBX47L5yG/cjsVtRUkx5OMt74ye9KQzpsWV6kFTbd/9nNSamKf6V6rqmbe8oE3Mp7K8rlv7yUVT5GMJ6lorsfjdZOIJ4n3dp95YRWounBzSsUlQQKnKA3RYfM1E6cH38ad1NRVkMnk+NBrmvnXR9sYGRwhk86Q75tHDdfFlE2f2Q04uTWTy5JOpsmOjZiAarNMUtDZzr24gHNxQWitTWD2l5Ma6INIP8MjFYRWr2X1zTcTHooyHh0n1d56ZiKO1qUXNJ0eqGmB/pNmqtA0fFuu5823beKnDx4gvu/piRtsFrRs4Yufeh190TRf/9fHzMVBb6up8VtcsLx40XG68ylkIS6c8yu5t6AkcIrSZ9nNtA0gMjLOx992Bd96rJ0DT70MTjeEe0orS/N0xRaOzYLG9eaPPxXH1byWbCZLbmRg7vVgi1NS7E7U8itwup2kLQs92gvhbkZf6Ga0sonNN11DQ00T+49U0v3MEyYYLXbmaLDWdEVnkq+WC5xzgFYKhrumFlxwuMwxvEHsDavYuKWReCrLlVev4IlC4LSvvYp3vnU7N60p559/fYqXHnjU/C4o28QFSy4LFVX4G5uIdXdNTE2a/NjFceIS+WAWZ5J5nELMh9NDsCJIqDLATVc18cTRIV558DHzAbmUxqbyORjqwla/AntdM6l4Agbbz20uplLYHebP1+l2kqpfN1GAfriT/T8d4sT6Hbh9bta+8Q14PHZOHR9g7OWnF6dIfXk9TVfuQGtNfDxJPBonlUih2/ad+XiTF+IuBqp0YuICQymwO9lw55twu+3Y7TZ8HgdD4QS94ThXLC+n8rO/SyKVZV1dgIf3dHPfNx+auvybPu33ItKPa81qYqcHzapm0xrtbTVdtpZDCiaUsFK5rpHAKUqXNwhVTTSvb+GL776CjkiK1oE4z/zwVxNdl0tNOkEmNg4jhbG2YgurrNp8aCfGTLflHJ6fmhSAXD4vqUCVyap1uCEZI77vaeJOD6PNG2lYUcfmbY083bdm6pSYs7HsptiCUmYKTrFogTc0kcHrr6Dlmp28ZkcD46ksHqed1q5Rjh3tZ7TDmlj8utCl6g/6Ka8KkkykUTZFMOihrMyF1pru7jEGuwfJJ5McfvmUmV4TC5vHTCfA4eLgxqtweV0kx5M81HFibuUHM0mGW1tNvd9kbCJIJ2NY5TXkPAHzOJZ8JJayyyarVillAbuBbq31nUqpFcB9QCWwB/ig1jqtlHIB3wF2AsPAe7TWbYt9fqKE2J1mGsFQB3VXXkNZmQu3286G5eU8fjLCd/7xRxNJHDYLa9V27E47qe42M054tpZbMcHmYv7xpRNTp0h4g9jqV5NPpyAaNq2e2GljvTbL/D+5dZ1Nk04kcfk8pt5tYswkzkSHUA1r0d2tJkCkE+SP76FzoJrx2Bbq160gUl3J+LED5vWarcVu2am5/mb+9qM72VoX4lD/GH3xJK0DSaKpLD3hOPFUliuay1lf4+Ev79tPZDhCdX0l1dU+mpZXsnL1W4hEUpw42AZac9MbN7O82s9vdncSDHnxeOwMDMTobh+kvrkav9/FsM1GvljU4XSZFOOvPMX4fF93rcHlo7q+msHUSug5YrbFwuROjk28H0upB+Myo7m8WpyfAQ4DZYWf/xL4itb6PqXUvwC/A/xz4f8RrfVqpdR7C/u95wKcn1gszVfg8HrJDHSb1snZpomU13PTLZvJ5jbxhk3V/NV3dpOIJejri1JV5TMB5dWuuDy5VJKyihrKt20nl80Rj8UZ7zhlHqvYypn8QRiPXPy/vMnn4wux/o03MdAXYWRwBD02MH03as0KE3CLU3XGBs3z6DpMKrTMtJRShVBS2YRSCn36tJ6xQcJ7nsFavglvwMv6N97EqWM9pA4+d+bjBSpxL19L04partpUSzST5ZETA3z/+W6Gh+PU1PhQSmHZFJUBN/vawvzHTzoZP/AceIN0DfTgvW4n2WyO7rYBVm9o4LqbrqCvL8rwaIKnn2oldeqwCfZFWnN8vzK9DMHa83uNZ6JsxGPxMws+XKz1VcWStaiBUynVCLwJ+DLwWWX6lm4C3lfY5V7gS5jA+ZbC9wD3A19TSildKmlUYv7CPVSt2kVoQyOHXzgMfcdeTfbAKqy8UcygtOxsvekqXruukl8fGuQbDx0nk87gdDtxux0cP9JtWlTFyjpaQ8cBRjowxdoDVXjqGvEvX4nDuY58Lk9iPEF6bMzUhI2F8axYT+LQCyXzQVm59Ur8fift4wnsDjsZX7kJgjpvVvQoZnvaXShfGToem0i2sVkmkE5OKiqseZkPF9bQrF8L4xETKLSGdILc8T1EgSPtjay5bifHBlogk8SxrBmn24ndbqemLkQ6nSMSSfDwkye5/6djr2YzZtoPc6LYAlY200swubUfHYZYmPHxzVRVeUn1dnCw7TDO5etRSpHL5ciGB6cvhai1eX8tB/hCC1sz12ZB10HGuwsJQG7/RIayWDIul67arwKfAwKFnyuBUa118ZOrC2gofN8AdAJorbNKqUhh/ymlUJRSHwc+DoBjDnU9xcUTCxMZjlBV5SVYV8u4z092uA971TIsyyKXzZHNmDqpNcvraKzx8x+PtdF+otfcX0NiPEH3qT4zfWOmurOpOKQ6SAx3mvG98npsgRA2y2YSRtIJcHmprqug46ijZAJnJpVh96N7TYAszPFk2WqzSLZln8jyHO5A57LmeXgCZjy0vM6MOU7OIC2WAXR6wFeOy+8n7/WSKR5/ctbocCfHns4TXLmG5hVVjI+nGeqPEBkMExkaMcfsPWb2LQZrb9C8dpO7xItTPfwVE5nNWtPT1ktfh81k2MYjpE/smygEcbYPv8SYuXDI52aeszkfdqd5zcLdrx7L1riOfDY7UVDf7py4CCiRD2dxGg35EnlrFi1wKqXuBAa01nsKq2ovCK31PcA9ADZvTYm8jGIKm4Vvy7WMH9pDPBYnmcya8nMOC2tZA06XE8tuYbNs/PXvXMXJkST3PXaK3zzealo2GtMy6W2D6BDZuX6QFVpV9J8gP2iRV2oiSCaidDz7LJTXm2CaTReKKVRe+CzKwiLWY8cOTa2W5Pab0nvKNqUW7pRAVdGIIxAwCUbKBpYDx+ptZDpbJ1pogSpcy8z1qM1mM5mj+Zx5bSbPMQ13E4n0s/8lTIusotEEbYcL+8otZG3W1IuV6S5cfOUmcE7u/lSKsooyIvt3TwTd+azrmUmZ9yi0zFRZ6jtxfuuCKmW6t4uJP7ks+ZMvT83uVcr0TEjQLFmXyxjn9cCblVJ3AG7MGOfdQEgpZS+0OhuB4idEN9AEdCml7EAQkyQklhqlyKQy4PbjdDnpbh8knzczsCzL4o2vX8sHtzcwkEjydPsYP3zoEKlECn/QT6DMw6nD7dB1eP7daJZ99iSP6LD58C/eZjnwNLaQGB9ZuAWbz0YpVt5yG93tA6ROHJj6HJOxiZVOZrivzeMl09s+0Y2rlAmiDveUXdPJNE63ExS4KqtIDQ9Nn6lbfL1y2YmkJK3JDvefGSgtu2nNJqITrVebZQK+1ua2dMJ0F+fy5/+ajvaZFmuw1gS1udbxPV0mZc4zUDUxv7T4O1BWXRgykJZm6dPkS6To3qIFTq31F4AvABRanH+otX6/Uuo/gXdiMms/DPy0cJf/Kvz8bOH2X8v45hKVy5I+/DzYLHLZ5WitsVk2nC4nn3z3NkJuO5//z330dA4Ti8TQ2Syk4nh8HpQCZVkmucUTMIFl8vy+ySb/etid5gP8bIUQJnfTJmMk9j+zMM95rtx+Th3tQsdjJvu159jc53FqTX6oe2rWrdam1ODk12a4E23ZwVULxc2ewNkzRostVqXAW3ZmDddJBQ3Uyu3o/jazfqjNMuPMZTWmZKKyER09S8buXCVjWPWrySUS59czkM+ZpDF/xdSA7vLRdM01dO47bFrboqSVSkS4GJOW/hi4Tyn158BLwDcK278BfFcpdRwIA++9COcmFpLNIhfux1HXSCAU4FNv38RANMM/3fcC2UwWndfoTBpGeiEWZqT/BCNOj/mAdvsJXrGTTCpDvKfTLIMVKMflcaHzmrrmakbDMUb7h83Ymd0BXUcu9jOemd1pujWD1a8WBdDe4Izl5WY03aoek+eCgnk9igk8GjOOHBmY18O4An5S/oozW6mFTy7d3WoCbKLQFevywWCbaXV6g9hsNvJnqz0LJgnIX2las4NtZ2Ze57LkBrtwLGsmU9lkphPNdKFhs0z3ez43fQa33Wm6twsCO15rFgl44cUzpwAJMYsLEji11o8BjxW+PwlcPc0+SeBdF+J8xAXiCeCqPT2HAgAAIABJREFUbWDrzhY+cG0j/z97bx4cSXqed/7yqMo6UUDhPrsb6O7pa2Z6DpIznOEhjnhJsinZWsryWruSQtZatuW1JceuZMeuQrG7YTm0u7ZiZWslW2sdYUmUSMmUKNKkRryGwznYMz09faMbDTTus4C6K6sy89s/3kpU4exCNzCD4dQTUQFUIY8vMwvf873X8/7hq7N8+4uvyKRmRYQEuo5y7NlnuHtrBi+9Ii4z3ZDEoit1/TiVQmWXKQHoBhNFm1hrjLaeDpyKQ1dvG6HHhpkcX6KQLeBO3RD3nlOWLM2+k7T2dLL2xsu1ps2+8PpB1e5pmkzU4bi4Ut0KzI7WrN79brIdS2K1tdXeK7CXFyReuJdrrOrpEm0TK227Zb5PiJtJzClDfhWv0fOVi1ApEUwOUHb6YHF86zZ2ng986BTdrY/xl8/fIDN+e2MM2Ef7IFZnD3a+AHff3P5c5aI8l9ZeQpGQLCoyi+LG3elamzgUkLZib/coBE2ZjCYOBh1DHHvPeVZXsjzU38q/+fOb3HnpwlZ3W7lEajlDz1A3s5dXxKIIhu9dKjB+kZymCSEFLIq5E7QkW7DCFgDZtj60SAxli+VhhsKitNM+yJPPPc5AZ5QvfPFNnImre7f6GkHAgr6H0AMBvEJOCOEgSx/MIFgR3IoQlpOuiiLcDxkYJqpUlPIhpTbGju8Fz5W608xSYzq1FRtWZykrJQS2HTSdj57pIBkJ8PyLlixGQjEh9libEGF6CbLL2OUilBpwfceTLL3yjZplaleTrTZL9TVxqPBuKUdp4t0Iw6Tr9GksyyTWEsHxPO68+Mr2XSgKadIrEXJrOVib21rusBv8LNpyEefWBVL+Z9XEFeW7KwHHc1mtfn7hvyxzIRiWbcMtkk1aSIPn1NR0dEOyOrPLG4lVN2qv7ZJfIgmM/hO4xaKQ5upirY5yv1HfZivaBtkVnFJuYwLU/cIMQOdRsVYT3eIdaOSYvpW9l+4rfqeXnVBI86t/eImAFWDp8iUhRs+FUha9rYuAFcAuV2O4jXQ6qYpHbLieunriw1Ku1MQmqMPjEGgSZxMHgkq5QqlUobU1zN2l3EaVmHqsSUsst76bxf2g/j/K/33z8fzPNycQ+ZqrVlRcyIDeO0I4GiafbYfVBfncdcAwMeKtWGGLwvysuIHzq+KGTXRz8skzrCxnWbk8jjd9wHGzWBJC8SrZOxASEYji6Bs1S6rrmMQO90IGShEIh6gYhggsNEqaIPejvnH3PiF18SVo7RVi7jwm1ml2GW95Bq9nCC0cRe0leWi78bX21GK6zU4phw5NV20T391IDpDsiOO6HsGgzvkjbbwY7xDrZTso9fa6yLaJOXrjl8j7iSSeW5tIlYebT1KwCxsnXzMImSVGv/iXb109YCguVjpIRmwpR3ElWKutzC6DacmCYC/xVKdMpWSL5RhJCGHtViazGXUJOPsC33W/MFbLtLbzQtJOWcpzHPveko73IsNyUUg/UL1nh7lV3bsUh6XQokmcTewvAhZtw8OEwwHS6RJ9HTE6Y1Ur5J2E3Sw014H2QZnIffJ8O6TblibkZ2uPZKam56VWtZgVV7Ly7k14oZiQUCwp11UpQWsvRiCAW8iLy3WlwV6hVZgtCRznaG18D4pyUQjccyEUJ9A9QKVckczd9IosEO7lGg5YIn6xNLEzeZaLsuCob3HWRBPboEmcTewvOo7Q1ZPAsgxSi2vcbQ3REgmI1dAINtcOHkY4ZWmgHbAOJrGoY0isnUbjhLF2gvEY5YW67ddFHsydZe78GGm5CJFWQp1d2EWbQDBAeeZOTdXIKctx6rWCt0OiG4oZdF1HT7Tj5Vb2RzIPIBDC7B7CWV0SC7NcwDODMsZG3NC6URPHj7ZKbLte59cwpZxndVaeq1JNHdtDiMPiqt1nn0oT73os3uHmd64xPZ0mGAqSShWJWiaJwYHG9g9YtRrEtxObxRbq0dKJHgiIq/QgsDzZOGnqBka4mjHsL07CdeOq1sQSS27d1wjIIYbPQyiKFbawwhblyZtCMvXE4To10jRM6D0J/adr5+w6RrDvKMaRs2i6RsAKoPefXD/HA6OwJtrDkRZJAHIdjJ5jYok2ArtQu6dWlMTRYfn9yCN0feDjsrjwidUpi7u9iUMFP8Z5P6/9RtPibGJ/4Towe5OF7DKRY6doTcbI2w7HT3bz2oVdYkyRBHQeoedYP57rsfjKN986GbydsFNMbHkSz4/XWhEhj+zbVEDvubhry7iGWXMvtvbK7z7xtXTVCvwjCbkuuwBdx+h9aATT1Jm+M0tuLVdzz/o1rtuh5wQEgjB1tWbt6QaJ9gTZtSxKKVzHxVuc3D+rzXMpryxCrHW9mbamadKYulGr08fqLOmrGTmO8shnqy3ZmhbmoYf6bpfca+JdjuwKhdtX8Ho/SDpf5uRAK6+FW7YmqZhB0A1Cw2fo6msnEgmiabDY1idlHPcTa4q3PziR+YSpG2LBFdIYx5/gw8+d4dLlOZZf/CsZe89x6SLydmZh+rHEWFJckH79pQ9f1B6E5DV9PYN4aWYJTddQK7O4vhD8bslaVgQ9HMHLpTe4g+N9fdhFm9LkbYzuI7hL07UeoXuFb+1vkVSMsj5vxjtAo6ZQtNmFHI7Lddr5GqlqmpTY6AZmS5t03Jm8Qn5zCPeQJKA0sQnvhu4oTTRBMUM+V2KwI0p7xETrHUGNvV77e7SVnvc8TTgcIBoNUqm4zM+uimpN/WS/V+xXNqSmgRUhOPQQruPy3EfPMtAepevZYf7wtZjE9Ep56BmBlen97R/ZCHRDyFLTobAmiT3bxTPr72O9ZF9mCYYfpVIoCOGmpms1m9tZ+9X6VW9pWlzqPuIdlEtl7EwW8qu46WrcNBS/vxiw76qvtyJDcZk107Piek504ty9LucZOC3XXE+S8Q6sjm6cioM7fVMSpvya1+wyjnZi/9zITbxlOCxrmiZxNrH/iCWlafTCHKFwkIqrGF3Io5Y3yaRFWkmvpFHJForFCivzK1TuVgvTHyQLd7/+u5IDBDp7MQMm8bY4l64u8PXFVcqzE0JG6QUwg5iJpLTg2gzDFBI4qAxN3YDWbtGLdey9E7dTxlmelevwCXc7yTsf7YNisSpPyNW/z+kF7EBI5O6mnbpjPICkYCgmZGcGZQGjG1J6E26pWbK+dKJuEKgubrzxS0KeTlmEERZna4uAUq6WrLQwJudo4h0DiXEeDuZsEmcT+wvDhHgHQ8PdrLbGKdsOL745RyJh1YQGfKzNUUwvUPRVeA5bNm12mYpholo7cF0Xz/WojF+RibnrGMzehEgCTdfWlXs2wHUOVoXGKW9VwNkrttN83Qm+S3izuk4oRqC9C9dxH9xlbZjoxx4lEAxgz05AJIEZb8HJZqqNtCuYR04TjUdJGwGJs87dotJ7AiNQtSCtCIRbRI0qlxLydZ2N4/Lc/dcKbuJdgyZxNrE/8LVDQzEwDG5fmyIUCUlD6orDTz53ntdefhxn+pZMWEod7lo5TSM48gjltVWc6VGcrqNgV+v7QjFRqqk2eq6MXWqcvPY7Fvp2LDbqSTMch7Y+NE1DN3Uc31K8X0QSmAFTSLilAysWAw2MtjYRrZ8Z52f+0Sf5vhMdvDC5yv/1Oxb2tZdh4hKuGZR723FEvmOzNw6Pb6+JfUEzxtnEdw80DWPwFIZpUB6/CqUcrhkk39IJhTTu0TNcnS9w4swg15dma103HmSCPWjEkpSLJbGynDJMXKr9rZiRXpSJbqn7a3QBEEmg9x2Xov1ysTFdVRCyjSXfvszd3aAbsDRBOROV5/qgyK9RHntTrjnahl3IgBVF80X6PYdv31jk0t1VZuay6w3S17WF0wswc11cytG2vcW7q+dsKgYdXhyWdVCTOJvYF7izt3FdpxZ/csriBrQiBK0gn/n8Gzz+3mGZ4KKtu8fS3i74JRhmUIjt7pvb/6cqJXV+yT7RqW2kbCZgiUvTCkBXDxW7gpdZ3N+ZoJH+l/sNP65aseVeJLrFIvXvyW5lLdvBz+yFDYlFCqQ5t1K89pk/E5LUDVH68ffzXCm9WZmqtRDbC0IxjN5h3FtN4jyMUEodmhhnUwChiQeHUutaqRuIwDCh/yHyaxlKV18mnS7RNnJcMlGrYuoNI9G9t326jjWuVgQyCQ+eheHHMEcehUjr7qQWbiEQiWAePSdiAB1Du4smKAVmEDuVkiSWlbnGSVOpe1ubwTB0HW3sePsBTZOXnwAFMk7PlXsRDMvnHUP7d07XkeM7ZRGJWBwX69JHakY6p/jjibdvzP69F3yijrbu35ib2FcodX+v/UZDFqemaWPAy8ALwAtKqT0oPjfxrkXXMNglWLwDSnH7xgzReFTKHszg3koVCmt7SrQJJDupGAGYG733xroB7YOYQYmtOdkMlO8xNiMgcb2AjpFMSuKQnZeyifrG0dHWqih5QSZ6TcNxj9bE2fcLmiatwN4qycJYUvRxA0FYnqol2mSWwHNJPPo+0itphk8NcuebqbfO/VmfWVwvAtEIPFcaoDdxaHFYFK8btTjPAL8JtAO/qmnamKZpf3Zww2riHQ1Ng9YezHhCYpnVGKDneoyc6BRlm73WaFbsexNCncVXmbgmBfPhXWTxEt0i/D1wmvd+8mmchSnUnYuSLXuv0o7lu5RnJ6TmFNANnTMf/TA9Z09D94hYYiCEUr9A2GuvynrsZD0Zppxj9tb+k2ZLZ+1a6pFfk2c4P7Y1O7WYIZ/Jk+xK8tH3DvLQcx96a0o/ukc498M/TOezHyP88PuxRs5J+cpe4JS3lydsook6NBrjdIFK9acHLFZfTTSxFYluaOvFmbgiRflVVAoF2ltCnH7/o1z/wvT+nzeSqDWiLuUgNVsrcg/HhbALaaiU+MiPfYp4OEDedphfyjN+Z6WW7dsIlIL0AireAaEgKEinS3R3xxgcfJzXXorjTd+skaRuCEmHW6TbSO9JSWRptPYy2Q+tXXD38lbL2yfLg4hv5lLbk7Hn7tz9xHVwbn6HlG7wm9PzxLva0Qcewrt75WAzqZ0ySikeOddDT1uEb78xy7yuU7z6SmMLimirKBLtV1eXJvYV78Q6zgxwGfi/gf+glDqE6X1NHBr4TasjCVirS9CYucHV0aMkEiGRPpsb3d8AhFMWyyYQkliXYfLY932IS98Zwxt7vdZ+qv8UP/xED7bj8fXRVWanV0ldevX+tHEdG3vVgWKWmXyBmasFOkaOMTDcy5xpiOUbbsHs6kfXdTRNQ8XP4LkezmbSjCXlGrYjwHg7eN7OyUoHhQexYD0Xpq+SzXZj9gzhHbR4/8oUVz83xdVYkuDgCZLdSYaGu7k509NYvaoVla43B1l728QD4ZDwZsPE+aPAs8A/BH5K07RvA99USv31gY2siXcu7MLW+GW1hVU0GqS7I0L8409w584xFi683LjVda/4nefyvr/zKX7uuRGuLOb49u0UbVGLya42VsYQ6ylgYYYsfv2/jqHrGoahkZq4e/+lMfO3a9e3OgvA8sIYHDuPFbIwh07TM9SFYeisLmco22VQ4MxObixHiSQw+0ZwyptKX/xjFzLSqqvnuMQRS7maazratlFK77BBeQwO9zKlwBn9zv4eWzfElVxvyeZSlG9dZNl9FMPoany2taKQbWbUHma8o+o4lVKfBz6vadop4JPAPwX+J2APaYtNvKsRDBM6dpq5mRRT44s4FYdEe4Ke9zzN/Ne/dO/9IwmGnv0Ak996YWfFl5Yunhhu53tOdPHMcAdHW0P84u9dxAyYEltULhgB2jrbsG2HuallUf3Zj3rS+slZKcitYcT6UEoxMz7HwHAfp8/1UShUmJtNY6+2SkmLWxESrJRw7l7dfmGgFGQW8QIWVlsbthWR/ewiLN+VbTar+RwmZJcZf/kC7afPsDLTen+avpomCWWbXb1dx8RTkJqpxWKrknvOrQvM3A037knIrbz9HXma2BGSIXs4mLPRrNrPAY8CY0hm7X8HvHKA42riuw12gVK+hBWyyK+ugadQyRZCIXP3ST/eTnDgOAErwKnj7aRXHyOfyeOMXthqSYRi/PFXbpKzXa5MrPD6t66J+HprryTWlIugGzgVh+mJBeyVZSHhzNL+X6/nYpgGruuiCjmmxmbIZdr44U+cJnq+l6XcCEvpIul8mRtX50in0jB/R2Jsvnasrx9rmJhHzmIExNWp6Rq6GcJduANtfQQ6eqj4/STvhb3WVe6GYIOkpBSszrJyqSDx70Zjyf5Yw3GItnHy/U8w+p2rtUzpYFj+lktBSyda5yCapuFl1ySz2a3Uur00gsMoMNHEBhyExalpWgj4JmAhnPhZpdQv7bZPo67afwVcVGqnXkNNNNEAxi+SntDWyaDt8ZOMvnRxd0tJ09F0jVAkRHs8xJlzfUxOppm5bWzdb/YGy6uz/N7rLwvp+HV/ShF56DE0TcMwDVYXV6UV2H4k02iaJP2szm6coFMzZP32XZ4Lrks+k+c//M4LoDzMcJigFcR1XFqSLbT3tOMkEziOQ7mjh6AVJD81LpZyMYNTLqMbIcoL05BL4bZ0glNGb+2kUq40HovsHpHkl/2wrPba7Nnvm9kIkWkawROPUb79hix8civMTi1L+YumSRw7FMMMmhhDJ7HCFqFoiMXLl2XhEUkIsRbSjT1nTZPEre28Gb7I/GG16Jt4UNjAR5RSOU3TAsC3NE37klLq5Z12aJQ4LwH/SNO0D1bffwP4f5VSzTbpTewN/qQZ78CyDMjcIzk7s4R9bRnt7FMkIgEGO2IsLRVkYtxsIShVm/h8sYTek2DnKaTWsFri5MeuyTb7VbaxU4X15qzTcJxyuGU9ScWpvgCWNA2CYayRRzhxZgClFJqmURxo55lH+yiWHUan1rh5dVqsSl8AAPDWliS+uZsVqWmYJ58kFAmRm6+Syr4Q5x6JRKlabLZ7pKYwVCnVRP7Xe3wGKGdzEr/VNDAtckspmLsJ0TbM/hEMwwANKnYFe2WJTKRFOrfEO2qKTp5b681pRdfj0FsQSUBb//bEqVSTNA8JDsLiVOL/9VdXgepr1zM1Spy/UT3Yv6++/7HqZz+192E28a5HwCJ69Dh37yw1JoKgFKXZSZ7/TjddXVFc1yN+4izZN17cmQCTAzJxGoZMfDPXsP3EymibWEtGQCZKP7nEz/qMJGQCbjQWt9NkXI9idudYqlJgF7CvvcyVySR6zzEM0yAQDFA83cVIZ4T/4b1D/MnpLn7/T3UKd2/XEovyq0IMO7kZdQPz+OPEEjEq5QpWMom9WH57G29rOugGVlKaUaPA8zxRVCqVpEVaegFmrskzKWbEwly+u05gnuuhaZrs45e59J6Etl6olEXAItYuzzkQEhGOnSzPgAWdR5qJQYccigeKcXZomnah7v1vKaV+y3+jaZoBvAYcB/6dUmrXUGSjxPkepdSjde+/qmnapR23bqKJ3dB9nKMj3Vz98td23iYUE91b/x9ldZbbX/oCt0E+07TdrUY/DrYytfFzX8jbrUCsDZY84g89Qvb6RYm/abpYQLnU1jG8Fcil8G6n8JDC6T+5+gpEEvzmmUcZOtqOYRpiUZULQsSBELQPrHdtoZDeaE32HCcSj5DP5NeJxkx24mjawesFW5GtwhX+c5sbxdbPYLXEcR2JB+uWjmvoOKm8bGMGITmA0dKGuzwri4RoKygPb34cr5Tb2O90YUwWPLoB4Rb0RDsBKyCqTpomf98OTlnINteMcR52PIBy0LJS6smd/lgNQ57XNK0V+DNN084ppa7stH3DAgiapo0opcYANE0bRsQQmmhib9A0Oo70kUoVZSLcDlaEM9/3ca795ZdqFmnXsFgcfsLMPWCefI80ad5cv1fd9xM/8hHGJtcYvxnk6aeOkTs/yO3byyxeekNIqe8EuC6szh1M8lCjcB3IrrD66tdYvWAIsWua/KwiGAqijp5B13Vp3Fx149LaQ7Srk8zcvFjF4Rax3jQd6+hp7FgrlAq1Ti37IU5gBiXTVSn0aBzPcart2Epy/HgHpOfF2vcc7InrYOdxom3ymW8h2oX1TFp3aVrG7bvFOwYxw2GcO2/K8QKW3I/8alXeMA/FDN7sbWwrKkTqlHf/3pSLh7vNXRPAwWfVKqXWNE37GvAJ4IGJ858DX9M07Q7iXDkC/MQDj7KJdx/CLfT2JRi9Or3zRGZFufHmRG0iM4PE+vooJZKEY2Gys7MiZrCLK/XJ959gYqKD+W/M1izUYBjsAnq8lfcca+VfPneCf/o5i1LZxXU9SSzJLovVMnNzb1q6B4V4O0bPMQZHelldzpJeSkli03oHkTzlbA6rpSot6DfZ1g0o5chfebUWO6xzFdujr4sakesIWe1XLWgsidXWJt1f8lnMeAIjEsYuFGXcZlBijbnUxvP5i5Nk/4Y453pWrP8sCmmYLeKE4zL2zCJ0HEGLxAgmk9j5ojS39lwhW59wd0MwjNHagZtbOdyt7t7tOCDBdk3TOoFKlTTDwEeBf73bPvckzqrv91HgBPBQ9eObSqnm8qyJvSFg0fbwE5TLLvati7XP/bIG3ZBXLoVXL/XW0slHv+chro2tYJo6M6bB2nIcJt7c0V378he+RXzoqLgvi9kNk683cZl//dsmfzDcg+N43FrO8PgTQxBPiuVVP1G/zUieO89H3j/M95/t4Cs3U0yv5Llxo4ela1ehXBTRhMwqZlAaPbtd/dIs3C7snjDklLe6atsHxV2512sPhtezi80W6clpmAZebgVndRanZ4R4MkE22y7i9n4Xlc3QDUh0wfQ1ee+XozibphrPRe8+ijd5Tb43c6OonuMQ6sC0gjgDZ9DM2tSmMinxVuw064biuPPjTdI85FAcWFuxXuB3q1ynA3+slPrCbjvckziVUq6maT+qlPo3wJv7M84m3nXQDRg4QzwR4eaFG+uTevTRZzn76CCv/v4fiaXnbiMrVymxsFbk9HCSXMnBdRVrExO7xzhXZ8lmlmqWWSwpllW1Y0b59huM3akqEYVb+PpqRuJcAevQkCZA6vWX+OzMAt++OERqcZXS+PX15BqsKI8/dZxYOEBPW5ivvTLF3N05uUblyTXvpStJNEHL0WEyd8fFxb3b/W3pXK859TNcXcfFmbuDo+ly/lxKnuX0dQI9HyBxbIT0pCmfl3JC1PWZx7qBbuh4yQEh2J6RmnXd0inXbZiQX8VbnKyVw+gGBEOUS2XU/B1o6UIRgqW70D0s+/htznqOy34zN2qZstnlt6ajTBOHEkqpN4HH9rJPo67aFzVN+3XgM0C+7oSv7+VkTbyL0XeKUDTE5MUrGybL/OQdxqIhmYR3iiVmV/j2l1+n4+gAXd0tKKUIdfVQ2i25xXU2lBAc+8Cz9PTEuPzGFLmLL2yMZ1WW3t445m6wCzB1ldmpq7UFRbIfq3eQ46f6Afjkw5185uUZ1pbXpI1b30khjY5BiRc2kvULMHmFzIwJ/afoPHWKpddf2VmlyQxiJRKQSEjwBnCy6e2zezWN1aVV1NxY1fXqCVHVSw4COGW85Rlx5SoFC3eELINhcdsGwxCOSby0tVsaiVfKkjE7dVXqB/ySpFBMkr1mb1Z7hB4Rt20pL/vXE6X/uy9heEjUaZrYisPyaBolzvPVn7+86fOP7ONYmvhuRaIbSllK129KiYAP3RBFmTfy9y5Sn7nO8tocH/vZH+Eb35lG1/WN2rUdQzJhbhfTMoNYlkkiEiQ3cwBdWQ4a9dfTfxorkUDXdVxXcfP6HP/q1iK5tRzlxVnJIg5Ycl9yq1Ja0yhxgiw2ZkeJnR3myN/6JK+/fKvmEvVhiNVor8XWtX8BAolWKnb/1oSsio2ybXnG/vPSjQ2dc9Zh5+Xlny9gCemlF8RSLGaqwhYKKxZGqRBO+Bze+KXaQkkpWTAU0rWsXj/Lerc4brRNCNdPrmri0OEdJbmH6NP+beBo3T6H4wqaOPzwRQ6qmqsAmEES599P+soFccNVY5s7QtMwB05QcTxKRZvC5FhtEo4kaBs5TiaVwb3zxtZidbfC6JW73PjW0jt7Uowl6T81zOLMEna+wMTtOZRS2EUbFsZrmba+Nb3ZomsUTpnxrz5P4OMf5/u//1G++e046UuvCJkFLPQj5/BW5ta72zi+Rmwotr1QQLxdutXUo3tE3PXpBTluKLYlgan+Woz+E7gzt6Q+d20O3Ar2vJQaBbsHKHcelZitGZTjZZb23tTbKYtV28ShxSHhzYaJ878Aa8DrgL9MPCSX0MShx2YBdACnTHq86rrzFWR2Q6Kbrv5ObtxdZWV0dKNr1YqSz+RxM6vbT9xWFG9h4p2f/KFpzIzPSsPqUo5SS6e4LIu57TOMH6S8wi4w+vxXOfYTP8Rv/uNn+dzlE0yv5FlYyFGpeEylV8Sqq5Rqz9TvhboZ+bWtLly3QqC9C9XWKX1b3Qr0noDZUan/rCNRLRJD0zW03hGCoSBeLEpldly+A1aE8lRRZPmgRpzBsHg69tJbs5QTa/ntFIdoYkco9c7rxzmglPrEgY6kiXcXQjGpu9utIXI98qsszS6zMOWJu61+gl6bo2wGa51CNqOUO5gmz281sisbCShzwLHZQpqvf+0Gzxxv46MPJXmi7zjpQoVvTq1w7XwfsZBJyNT5ynemuXNjCm9pWrJyd2q0negWlaPVWcgsUTECVVUgS0iylK/VXNYtpJRSGKaBYRo4ZQc3vVLzTvgt7JQnnotIoiYCcT8Nqb8bviffxTgctNk4cX5b07SHlVKXD3Q0Tbx7sJuFuV23lGAYXdclthcMbyy1UOr+Jsm9wjAbFmD4boF97WX+118ak4VOtFWuv5ARogrFSQ70MHy8UzqS3Ms1bOdlv4q97p6PHjnO3/7+h/mDz13AWVnYaMH6WBjH7hgS+USnIkTZPigLr0Iaqo2rNV3DXqyLYSb7RfQhvbC37OImmrgHGiXOZ4Ef1zRtHFFy05YsAAAgAElEQVSS1xBt3EcObGRNfHdjt9jTdu5Wp4zrulKjuFOm50HjIIS+rUitXtEvg/HFGmLtMvH78UFfYq6Q3r/WYI1gF8s2NX+L4eN/i1hfH7mVqd3HVW/5r81DeoG8M8yN6UGiLVEKWg+V7HJNNUg35N6AZAn7+sLFjNwXz5X7F+/AqTh4lcrGpgG5lGxv57eOpYl3JN6JyUFNNHFwMIMyMe70j+GUceYnt4/lHeYmzruh5zhtR4Yo5os4FQdn7NJ6gkpw+GE0XUN5irJbbRvWMUisrYVcOgd3Lt77+G8FPBelFG0dCXJDZ2HyauOkrhQs3uHVzy5D11FIL4k1qmnQ0iVu17lR8U609aEne0QEPleNr2aXZfv5W3i+XF/99+d+3bUg7mMzKER/SCbrdzsUB9Md5X7QEHEqpXYIHjXRxD6hfRBS0xsbONfDdTa6AgMWRNtoPXkawzRYeemv93+C8+NlB4WKvS6+vvnzcqGIFYtQtsswf2u9h2ku0324LCjlsbpapKMjwtykiROw9mYNKyUWpOfWhAja+gh0D1C5c7m2TWoGr76sZnN/zP1O6GntrbYhm2m88XYTB47DYnHq996kiSbeAiyMyUQ4cEbccJvJpB66Ab0n6H3scf7NP3iKwaG2gxnTfvSs3A2ZRcrjV7HHr4o17bslvaq4vEJE2P3JwnUk67PRdmdvBZTizrW7nDuWpKu/UyzFemgatPZIYtBOiLahBy2xGKtu6opd2Vrn6Yu81/fH9BtWb5dBrGm7f492QzSBEQ7X5BqbOBSo/wrs5bXfaBJnE28f/EnPR7yDYCQshejbwZ8EPRemr+O6Lg93Jfg77x8SCyGW3DpB18fJ9oqDjiNWqqIAxaxY0/UZnZkl7HQazECtT+hhRSnPy5fnOfNQpyj7+Ii3w9DDWP1Ha43Fd4CXz4rQQSCE2ZqEVIOiDYFQTXoPqi3F4nLu3pPiyYC9E+jyFO7cnVoc9bA/g3cDlELd52u/0STOJt4eRBJovSN87Kc+LRmWAMqjnFqWuNR2X/a+U2K9VLddHJvkx377VV66s0bXuYelGfHmVmV+/83N2M4aaR/c3TJ6qzE3Kq3N+k7dv+X0VmBpgvGbM7REAvSPDMhnsSTBgeNYsSph+gk/2yGXAjtPsKMb48hZDNNA6+i/97PQDRGA961Pw4SjjxA8eobw0YcIt7fXFmaJ7r3dw1yqlonrOpIN3EQTVTSaHNREE/uLaBvKcbg9sSqZj3DvukS3UlPHUQrmb3FnNMH1Fy6Ixea526gGOduXSbQPYnX1Yd+9UdtPedvLwL1d8FzIpTB7hnA2x/QOE+IdBKwAJ7uiXAwa1XijJ7FbBWjINUQS4v7ergPL2jzlYBgzIYucYCiIfS9LQTeEEP1emq4Dy9NoAyPrVoYZaxFlowdpmfZWZjA3sSPecclBTTSx78ivQnaZO9mVrcS2k1Ta4vjGlb/rYM9OykTdqEqOpkEkQfLECfKZPOHhMyTaE6QWUpTvXD5UnVEAKBcwDAPnMMfZ3AqapjGxUmT8jRvrQuv24ix4LlbvIIZpQDyBMzm79VlpmliLpRyObmC0tdU+3w1OWaxyTa+V9BQzuI6LERDXqhEwcNoHay3UrMjhe8ZNNIxmclAT351o66v9HoptnPysCAw/JoLs5aKQY337Kt2AaCuxR98vzZg3Y7ts27W5rdm2Rx7BPPkemYw3IfHEh/jgj/4Aff1tBIIBNE2jmC/S2tlai4dthl61ogbOEnnkGSmsf6uQHOBjHz0j2q6HFdkVCvOzBE2DjuEjcr+UEl3g7DL20gJ2OouTy8gz0Q2JQfolH90jxB9+H+1nzoGmY+eLKE9hdg2SfPo5iYVvR6LRVklGqh7DX0A5uczG7fz4pO+2byReaZgSM2/iUEHd52u/0bQ4m9hf1JcMbJIvM4bOMHxqgGBwmDPHkozPZ7l04Q5tXW2szK3g5rP82I9/D8PtYX71dz0Ku7UN81GfYVmt/+s60svia69s69rMZ/MowLYdWtrizN4ch1KOUO9ALQ7nk3PAouupDxO0TJSCx852c3tyjbFyhcrmDiAHhflb/MV/zm+N3TaK1h5xfzcqdm4GhVziHVIe1Gg5TmqGS2PLPPOeIT7/ZmtNGrBib22YDbJwirRKCYoRwKk4ZMdvr39/yr0niXZ38dwzw3yl4pK+/sbWbOJiFixXLM3ZqqUbDGPGWjZuF44LSVfsrZ1btoOmiRV7mNz2TQDvPK3aJpp4YLjjbzIbDfG+9w2zmC6xulrk0SeH+cWPn+T3Xpvh8ugSf/+JQT53fZ5C+j4E2cMtGMluFu9M7yjo7izO8PK3NALBAIn2BK39PSilSI/fkfrI/jNSN+mUof8UP/ODZ0iXHGbXbL78tVHSS6laYtFb8U/sOtIqrGtYzrdXLdX0QsPjtM48xS/+1NPcTdn8xddusQwbW4HtBs/l0iujfPJnPszzI6fJX3pRzqtpMHgWFidkO3/8fkxSeTB7g+KcvtENn5omr2n86edLuK67vSvec4V8nXLt7wEL3djoSLMiFnYo1rg734rKwuGtkHFsomEcVGnJ/aBJnE0cPHSD+PlnCEVCWKEAr7x8B8cRa9AwegmbBhHLZHZyiU/88hclqcTdQzsoH4U07u3Xdv7vMoMEugewwhadPW1EIgGujk2KnJtThr6HCEbClPMdsDZHpCXG73zlNqmlNEopclcvrLfWekv/g40AZmsSxwrD9PW9tcrawzjtTJYXRld488o8y1feFGtwN4TjYp0uT8qYpq5S8T7EP/w7j/OrM9O17OhijsTDT5KemJDFSccRgh3dIpO3PCOJO6rumnRDCDGXws2lhGzrS058hGLVVmaVWvJPfo1yJoPVmti4rRXdvsn2dijlDpfIRBOHDk3ibOLgEQxz/vEhnjqeZCpV4uuvesxPzNLS3sraaoGf/M2XUUrR0ZOUmGPBJr9wn1mQ9yAKTdNo72qlpcXCdRXD54a5c7GaZLI6hwofW4+lFq6+wt0rSt77rj7YarmEYgfXVUM3oH0Qz61aY37brIPA3E3++o+rYu2+m3tzMo1uyL2wCyJUEWur9Th1Hf7gK7d59ol+Wo8dY21lar37TdqK0nHyBLnuXnRdpzBfbbq93SLAj4O2SF2vUgrXcfGmb8p9jiTEOi1m5bltPkYxgxuLoes6mq6tf7YnHBbTpokNaCYHNfHuQCgG5SIvfPUq37i2RN52mL85BpNX8DyPSDTI/KVLLLxxkfNnuvk//u4jYnHWx8Uara3UDUk82ikb0ylTTq9hWSbFohBDKBQgeWRQEn6Uh2EaJE+fk3N6bs19uJubb7+K47cbt+dCflUmjLWFgy2N8JWJ6mPDmzNQo621BLBSDpanNvx57tpNbk2vUcgWNt6X+Vss37xJaW6awsI8WryaObvdRFixhVRnblDO5qisLIr71X8WypPksVCs9lkkIWMDSM3gjF2ifPsN7MV57KWFZibtdwkOi3JQ0+Js4uBgmGj9J1HZVSiXefUbVyG9uO5WOzrSzfHBBBOXE1iJBMc6o4yliiJiXm9F5Bp0sXlu1e23y3/K3Cij38qhJ3uIJWJUyhWKKUm8sQZPYJgGlXIF2vvQOvqljdnS9O51gPejZ2uYksUbCMHMdXnfNSzuzXpyNExoHxDizKX25qY9CGRXate73XWvzqLrTxBLxEj5Y20fFJJLTUMwAspDLezQdNxPzAlYtD76PtbmlyA1jVMu1sivmBWC7x6uJfxYUcyufpzJG9A9TOeAuII/+sHjTCxkefVzXxJr3c43rcl3KKSO83A8uyZxNtE4zKBMgI32NgyEMAOmNCyeuylxprpY1djNWW7fmIGFMewlg//nP9mEIiFKEzc3Hmcvhf/3ssh8wXDdIOM4EjN0yhBJoJTCqci5rLBVd/rqMas1oCi1sdmyFZGYnJ9RHI7vmJwEQDDMmR/4Pv6XT53iOzNZfv33WyhnMnI8vx6x/tqnr0nZxdtNmj52ex7lIq//1at86Ac/wDdmT8PUVYgmhAxbe6XZ+Gbr3e9z2j4IkRZJJCoXSa+kIbe6vYhFdkXucWuPfK+CYR577wjfscucfGSYMyNJvvnyBF/99jjFfFE8EYWMfP928h5svvdNHDocEt5sEmcTe4Cm7a0sQnlUUku1zMxiZsM3v3j52xu3n7xKqX7y0jSZzBrNhtwNvmatT3bLk7XYXKKbQO+RWjxsM3zCUkrGEu+Qa3Er4DokHn4vJ071cOEP/qR6nffICHbKFIsV/sUfXaZcdhka6eX281d3jpO6zsYyn8OOpQnmFx/j6ece4aXPzsHsqNwrK7o96UYS8rkRgMzy+n1QCxNbk3RCMSkTcR35nvjP0DD5zl9dgNVZRpenGL1yBOySxED7ToHnQGFt5+9SKCbu+aWJw7NAaWIjDkh39n7QJM4mGodvMZhB6D8FK9O7W592ASo21uAIdrZb4mG7uTWd8sYVv1L7N4kZpkzQ2yXWFNZwnYEtZQzr6BhEDx7HW1siNngEK2yhaRrLo7ck8eXit7hwo2VPYx2/NgEL1RKYYPjgY3Ch2Nb7uxt88Yj7lPm7+dVv8tTf/jj/9Bf+Lv/hTy6Sf+NbtYWBH8v1J8Fqv80tAalcqhZPrS4cWh95D5lUBu/2a7XtApaQ8tyovPcTl6rXEYhFqWSzu9/jSkm+203SPNRoSu418c5DaroqOFCREo56C2k7N1csSefZcyRaI6ymclQ628lcfnXnrNB4O+imWLX+sfZLn7Vi7xynNIOYgZ3/Fax4DHtpgXDvAEPHuggEdDxPYQZOM//SokzIyqs11I4kxIW4PLn9tXou3H2z9r5+QjdMsWhD1S4jSxMPfg/i7ejdR/FyaZi/3dg+3SPiXp25fn/nzK7w8u/+IWuf+iGi8SjrdmMwLIlYhbSUu/SerD3rxTtbfXHFDMFjZylrGqRmWHv9Rbk/9c3LnbIkE9XDMCESh8wiylNbY8ebcVh1gJs4lGgSZxONw59cqlqkO8KKSIuozCKFXIGe3has3gSVikdm+USNNJL960LmdB0DK4Kmayg7v/+xpp0EC8wgJAe2NpPehEhPH72DHQQCNavUMHTiZ58kOznBUz/wLOePJfm9z13kqadHKNgOozcHSL381e3P29IpxOTH7uozUKt1hMGjZyjXuyPvB5oGoTje+KW9WVMLY5K49CBQihuvj/KBj51nafZx1Njr0DGE1daG7buzQxHIVS3h1t6tLulSjvLSXG0sTlk6orT21shSqdoCpX1QFnZGAArSS9N13Wbs8rsAIp93OEzOJnE2cf/wyciKyIRVLYSPnXmCeGuclfkWlKewbZf5mRSZybviVvMzJ1s6xH2rG2JpztxAOeW9u8uCYdlnt8kxkpC/b3bXOWVx0SUe2fUUnrexrZRlmTx0vJ1csQXz/CA//z0jnOiMYVc8bsysUSxWyK3tUtsZDMvknl4AM0jvB74X13FZHL0t9yezRLlYunci1r0UjJTaWi/p10lW7KqurLf1GK4D7j7Upk5eYXLqGK0draze0SC9IF1PKiU5593L8l2Id+wcP88ub7QI82uyfT3C8WpW9Zy8d8rivnXKqHxGvqPlYtOyfIfjkIQ4m8TZxH3AMCXGWcxJkkzXETr7OkknWtF0DTNgkl5J41QcyqUyd8fmsGcmxFXqT/SROKzO14hh+qr89Cc4qLWMule5RzAsFpzfENqKyARZ7ybdrHO6aX+7UARNxwya0smjHgpC0RDhcGD9ozNDbfzN0x38+fVlPv1wN48NtrKarzCTyjM/n2P8xZdkwt/pP31tXhYRAAGLlhaLhbk1GXOiEz2WwJsZ3V1YoaUTOgZg9tbu221eiLiOfGaYWA89iV20xQtwQPG9mYl5Pv6xc3wl/R4q6TUhak0TsrMLgCefbVcPa0Xke5Cuq2E1TB5+5mEufyckmbLLk7Xr9++3W4FEB5QLQsh+zPNeiVtNHGo0k4OaeGejXIJogtjAIANHOnEcj3wmTywRIxA0yWfyeLO3oZjBro9/+l/87Yis5zjBtnZpZr0wJlaJGdz+/H6D4nJxfWIM9gxSjlc7WsyO1rZN9gtB1xNp51GZqFempLShKsfmxNsxBk9A1XNrLy9BYY34I49v0CcYX8zyH9NFknGLkGHw+69PcXU+z8xCThJ/dusrChsTdYpZbn7xizVCW53FM4Pby8xtghkwcRrYbguqiVf27IQsTA4wKca5fZGeH3kvx072c3dMx3M9KqYFmUWxunOpWvzZ76jjL3z8XqlWtOYtqNhc/sKXpQRoc7Y0CNk6ZXHVeq5813qO31+9bROHCu+K5CBN01qB/wicQ1zUPwncBD4DHAUmgE8rpVY1CTL9GvB9QAH4caXU6wc5vibuE6EYGAGC4RCe6zFxe47S6hpoOpVyBd3Q0Q2dwNBDVGbGZMLbqazEDMpEpxuYiTYpCQnHZELcaaIzg5hHz0nN5fS19VZR5UIRKx5DeYpyOC5kaEUwO/tQ7b240zfF4tAN4oNDlAolKpXSxo4Z1ZR3TdOo2BWJuZWLrC2v0dsn+qeuq1hdK5ExNBZSBf73VIFMzsZ1FZPjSzKmvWLzvWkkJpddxinlNu7r389C+t5+LaX21uA5khAy32sGsOfy+q0lHjnRweiVCcisiMpPubi13KTrqCzK/NiuryK0GcXs9tajbqAPnsa7e0WOb0WFODNLe5cq9JWJ/NKXJt5WSNL14WDOg7Y4fw34r0qpH9Y0LQhEgH8B/LVS6lc0TfsF4BeA/xn4JHCi+nof8BvVn00cNgQjtB/pxy7a5FZWRURAeWAGqQTOkOhIYAZMIaaOIYKhIOVCsZYUFI5DKC7xKM9F6xmW3ph+HeXyDhqmPiIJIvEIdsnGPfYIzuoSaBpGMCikuboiE64VkUnU13n1J79gmOzt6+Jm9i1XELLuOgKAvbIsE3d1H2d+ktl2Ic6KXSEYCtLZFUfXNVbT0n5qYSFH9uqFty6OVp8U46OtD6O1A/fu1ZqubKPi5jvBj0lrmiTpVOzt46K7jPPya+N83+PPcOTEAP/sB5/jPT1JJjN5fuOFu3zr69fxigX5PlTKQubb9FJtCF5V0zbeQe+jj3J8JEkuX+bii9dg8sq992/phPYBWJ2j/4knWFlIURq/3iTOJjbgwIhT07QE8EHgxwGUUmWgrGnap4APVzf7XeDrCHF+Cvg9JUuKlzVNa9U0rVcpNXdQY2ziPhBJ0HXuYbq642SzNrn5hRrJVWyYu8WZD/9NvvdsFyt5h9/+4wvyt/SiNCzu7xdXXblC2XOhXMAMmOsWnrc4ubMgt1+CYOfp7mtjZSlLam5pXVzBzSzh+uLfAOG4kObCHcyBk3jdxyUjc2liY3srH54LMzcpa/rWmGF6gZWXqq7FYBisKNHYE6SW0gRDQSp2RZKfHrQes77M4n6wNoebWYT+U3T2d6LpGoujY42XoWwHMyiE6bvXo60indeoKIMV4eEnjvHhoXYe+3txfuPFu/z7L40x0Bvnlz9+ij8aSHDxTgrPe5TrV6bJF9JblYJ2g6ZB+yDBzl4S7QmOHE3yyfM9LOYqfOWVKVpbQ/z8z3yEf/8HLVtFN0DuuRWF1h6Gzp1A0zROHX+UvrYIUysdPD91p/GxNHGgOCQG54FanMeAJeA/aZr2KPAa8D8C3XVkOA/4Ct79QL1PZrr6WZM4DxMqJbJrWVzXpZQvVTNijRp5lovcGl0iHg6QLVbW6yPtcAtmPIHneXieh6ZrWP3HxC2qV0nzzhu149Qf03/fc0KK3F2HbNbGCltY0Qi2XZBJfEuXjCxMvAFK4aQWsLr78Nwold36LG5HfLohMdHUdLWPZAWj/wSdnVHuXrouSUC7JR81ilBMtH3vXtnewtE0SXLZLVYXbiFx8izpVJqlS6+BphM+coLiytT9KzBV7I372gXoPiax5XiHSOr50niVUk0PNhgmfuZx/sGnzxO3DH76d77D0vwquVQa5m4xCly59kH+7idO8o+/5xgR08D+3hHiwef42d99nbvPf2m9U82Ot+zc0/zIp87z4ZFWlgtlXA9eHFvlD742zuzkEsVsHhyHQMDgl/7+0/zCL1yt3b9gmPYnn+HjHxyhJRygOxbgaFsIU9OYy9nkyx6JsMHzbT17I/ImDgzvBq1aE3gc+Fml1Cuapv0a4pZdh1JKaZq2pzuhadpPAz8NQCC2T0NtomFUbIrXL2AffZiAFcBqTWDnO2sNkz2X5TdfZ64vwdG+Fk6fG+DO7SWK0Ti6oUtcUoGdL2BFI+vuWV3X8XzN0mirZFLOjdYmzUgCMxbDGToHnks+k6enP4nrJli8dXF7166fbeuUMdq6MAxDLFtNb0yX1I/Dea7ItQVCohnryoLgwpdfFsLeLzee8kQQf6fjhWLQeWSjeIKPYBja+tDjrSKSP/7GeiPpWOIxBj7xSW69elks6fxqLXa3V5hBiCUJWkG0E48SCAZwKg4qXv1fVFK609bZxs/9N+dwPcW//ewVlq5eqWU824X157X4wpf5ty99DZL9WJ3dmKbJw48N8fM/eIp/8sY1CEUlBl1IbxmvdfYpfuEnn+IPvzHB73/mVZRSqFJRyoucsngHqmUvjnOW9/cnsY6exs7l6Bke5KPvP8LRZIjf+osbLNwUizzQmiQUCaHrOsFQUL6f2Qa1mZs4cBwS3jxQ4pwGppVSr1TffxYhzgXfBatpWi+wWP37DDBYt/9A9bMNUEr9FvBbAHqk65DcxncZAhaBYK00w+gcwE10idWRW4FQnHPD7RxrD2PoGjMzGTRNkwlWKdTSFORS2LEkJPtAq05SQ2dR0zeg8whWJIwdfATKJYwq6eq6TrSzjfZOKS0wTZ1IJLhzW69QDKu7b/2tbugEggFK3SMce+QE4y+9unNyjKah947UuqWUCmBFsMIWSinstVWZoFu6dtdA3QvswsZ2aptRysHsTSGfii0kFm2Dlg4000RNXcdbGNu4j1IsvfYS0Q99mE/88AewAgajd1cxDJ18voxSipXFNNnFFclytfO1TFZ/MRJJQLSN2MAghmGABoFggGPD7Xz4bBc6GreW8uSKFfK2w5PHknzgSCt/9MY8n//LS1RuvFobz3YJPU4ZFsexcynsUIyXs3n+5cdO8v7vf4p4OMB/++QP8I3xNL/97/5cSL9chMGz/PMff4r/+MVRpl/4mnzmeynWS1JqRPv6N6+S+xtn+fQPnidumfy9R/v40u1F/rf/8wuSwV3dpzIPFah5Uu53gdHEvkPxLtCqVUrNa5o2pWnaQ0qpm8BzwLXq678HfqX68/PVXf4c+Meapv0RkhSUbsY3Dyk8FzudxmqVZBkzaGIGTYhFKJc6UUtTvHl7mQvXXExTJ5EIMXdtVWKLfqkBSKajUwY7T1k3RA3GDGIGg6CBFQ1DNIymaei6TsAKULErLM6tkpuZpvPkCZamF2TS3CwEoGlonYOsC41okJu6S2zwCKffe5r+njjjr0d3vkYjsJ4dbIUtqO+WUqkKufedItHbTWY1g5p4c/8n2M1lNIEQ+tCZ9ZiqbuhUSrZkAs+N7Zw1Wswy8eUvMNHWh9bWTUuyhb/3N87xvcPtFBwHTylKrsdirsxaySFve+Rsl/m1ArfGV+nrixMwdDpbQqzmy+SKFe6MrXBnbJmxW0sEQ0HOn+1moD1K0NT54itT/M7nL5OeW4KFPcRWC2l55VLM5T/Gz37oKEMtUX75r27Sn4zyK7/0aT776gxLS3nef76PfNll+vKNmnt9p/uvaeA5/OJ/ucI/+egI3dEQ/+gzb3DxL57fWWDCc5u6tYcNB9Rb836gHSSDa5p2HilHCQJ3gJ9Ammf/MTAE3EXKUVLVcpRfBz6BlKP8hFLqwm7H1yNdynro0wc2/iZ2QSSBOXhyndTs1RRmSwJd19E0jUR7QlyGiHhAaj4l1tR2cn26IW5RKyquz65hrGRywyYVu0LAClAuiZVEPrOxk0X3iFhM9RZNwJJs0GgbkcFjlEtlnLtXGfrQR3j0dBeT81ku/cmfbj9BahocPY8Vsbb+reqO1A1dsm8rtiQ07XedYDguZBgIoQ88hBW2CIaCpG/dkM+7jsLcrftLSDKD0NZHpLefcEwyi9s74nzq/YN8cChJwgqwWCjx269MMzqeknjhahpKWSkPUUqeWzgu1x9tk7/55HUvJad7wDr7FP/gR9/D1ZkMz//Zi7A2z9AHP8zHnz7CM0daiJgmP/3r32Lt1a81dsCo1HxqXUdQrgup2e1rbcPxreIZTTQM++Yf4xUWd9evvE/EB0+pJ37+/7uvfb/xz555TSn15H6N5UCJ86DRJM63ET7ZlXLyux9T0g3oOEIwHkPTNez5GdA0Ah090u/SLm+d7DuP0nXqIaIxi/EXvoXZP7JBvccn48Lk2M6uVb9DxnYWRCgm2rnpeShm6Xr6exgeaSeXK3PlL764M/GE45hDp7YoCalqFXalXEGl5g+m5ZemQUsX4aERzIBJMV/EWVmQfpWVsrgQUzMPPsG3dEptZn6tVgMa75AFR6VO8s9zt55LN0RYYPam/N7aI4lS++G2Bnmm9S5jkAxYXwB/u1rVeLuMf2GsRuKahj7yuLjdl+a2rwutbhc7/6xkis+Nbr9NE7vioInz8Z/77fva95s/9+y+EmdTOaiJvcMMyoQbiomFZxdqLbuUB/O3KK9ULbVqco3X1iVtu/LpbfViV+ZXWEyvgOei6zrKU3iuh2FKQo9uiOW4I3FuzvzsOiYkM3NDyH3mmozbMFm88G2W58Ra3tX3U8ziTFzFGzxNwJKYrvKUKBtlFuV8++3OC8WgfZBEbzearpFJZSiO34D2AYy2LoZPDaDrGpqmMTXRQ/7unQcj7nqry1czajRDWHlyD0IxiLZhtvdsFWR4EGx3HNeRz32JRbXp/mdXhOzrn6tSeIuTeIHQzqRZRblDcK8AACAASURBVG7q7kZBjCYOFQ6Lndckzib2DjOI1TMIGtgA6QXMgRPiSrSCrC6uioZooWrFlHK4s7dxK7ZYMZuRXsAtZoTghs6hlKKSScPCGE7XMGY8gbN4j16e9dAN9JakWBjhFqgs1ZpQV3XzvPRKteD+HpO8XcCbH8cbOF7r11lVE9pXWBHoHqbv+BCJRIjrL70pk7wZxBo+S/+RLqLRIHpds+2u3iTjy6tvX5NrpaREByDRgZNN766Zu1+olGRRtDa/ceGiaeL12M5FfC8JRJDreZBONE0cOA6Lh7RJnE3sHXYee24Sq3cIq7WNQHcXmqbR1duGZZkkO+K4rsfKYppKuYLruFQWpmuWTCSxkQQ9tya9tjghWY1+LeDcKM5CNVsykpCyi50SOoJhceUVs3h3r+D57kYQq8it1Ihyt1rOzcivUrl7neDRM1Ke0DEk7smdEIpJYk++KmgO4j6cv7V9AkvnUSK9/RTzRWZfe43ZUhaMAInzT9Pe2bKBMCsVj3LZxTA0VleyMFc3joAlLlO7IGN4K0jMv56JN3dWE7qXZb8bYkn5LtQvcJSSJuD1Gr26AR1D6C1JvExKnu8hmWSb2B8omsTZxDsZSol4e1FECCplmcAWZlO0d7XS0mKh6xrRaJC5mTVSY2Pi2gQhtpZOSaap/yfoGhYLcLvJ3rcqOo9ghiyc0dXtJ8VgWCbQYnbjRNvSSaDvGJVcvtaFZa/XW8pRnr5dTYLZhZA6hgh29pLsTjL/xkX5b2/p4r0feYSrb7aTv/FGbf9q39JoZzv51BpMXZXFQd9DJLuTRONhHMdjeblAMhnGdRW3r0/hzIxV24KVasRlBjGOPoy7ugj25FYN2IPGTi7rzeIRe4GmofccQ9d1nMnrG++7XZCFglOuJSoVM3ipmb3JATbxzsEBZdVqmjYI/B4ixqOA31JK/dpu+zSJs4n7QykHk5ex+09hxWPYuQJuyCK1lCYe78LzFCsrBVKLqY3WnetsL/+WW9nZ/RmOi8t39iZOay/0npQM3c0uOd2QBCE2uuW0zkFJ6EnP37/1o9SGLirbQtMg0kK5UGT+zTdr8djMEq/+55sQayd0/GHRPm0fwAyFcSqOkGm1UD928mGCVpDUxZdIrTcO95jrPyXi54vjG0kqkpB9jQBuuVyz5A8TcQRCsqjxiTNgQaz93i7m6jU4mdWNCwFNk+vuPAIr0/JMnPL9afL6ggxNvCNwQBanA/y8Uup1TdPiwGuapv2VUmrHbg1N4mzi/uE64Lq4jiuk1jVMJWCSy9lUKh7zr3yr8f6H2yWkhGLVsotjYBiihmOYhNvbKWaWtrpsW7owoxGcTbFLtTSFoxsPLna+G0Ix9IGHiLfGSV/4+ta/uw6kFyjlVoTgV+dwysUtk7ZSShYbvsvVx2aB8kgCkgMwe0OItGJvryh0ULAispi5VyKR54o7OdENpiULpHCLNDFvIDbr5dKceM85br2Qk0VBa08t03Z2VEi0ERWozeMuF6UtWbxDxtcUOXjXoqoXMFf9Patp2nVE7rVJnE3sM/yJ23Nwpm/JxD1znZJxjrkpJRbe5tVhJCGEsFOMcjOibdA+KJaW79rVdNEf3U4IfnUGPf6QEG7FrgmmN5IY8qAoF/GWZ0iv7dA/1Ifr1DJDt4FdsGHu9s7WtxUR681zJca3WdP3rYJSREbOUrjx+r0Tpar9ND/2Q0/z1eevYQZNPNejIapbmuDWi0VI9mJ2D+JkqgLwLV3Qf4pkd5LU5dcaJ0678P+3d6cxkqTpYd//b0ZkRp6VdZ99VF9z7uzN2SW5Eg2SS60ErUnLtEVCkAl7AX6xQYk2IJDyB8K2vggSREkAIWilpWnIxFISRUoLeukVyaVXJDWc2ZndOXr6vuu+86o8I+L1hyeyMuvqruzp6qruen5Aoisjr8joyHzyeY/n7fQBp3LShKyFDp4Zh92QYoyZBj4FvPmw+8UOdzfUcytoEUsmJQi2sw4vDSaG3/Rlqof3kMo8B7E2A3NXZU5du9lz4Qbce2/vDGGzQHMp6ksbOg1nX5PBJU9DO5B9xLq16VxaRobulwG1pFIQmwU5Fgd9PWOkSETfyEfav85+1KnOz8ggKC/96Puv3OM//t/fBAOTZ0ZoNg6434Ev00MWbuHGXZxsnxxrL8PQxBCf++wZOV69amzKXM+DrFuqjgkpufc4F2DYGPN21+Xndz67MSYL/Dvgb1tr91miSWjGqR5Ps0b44Mr2X+uNKjz4gMDECLz07i//ndNJ4p4UJli6/eR+9XdNUXAGxwj2K4pwGEbOSE3bRy0O3a6BuvML23GZOjNEaXFMMuq9gmcY7P1+HrUcmbWwHC2PNTApA7EmX5Lm0sepeLRz6kbfiPxIWbix95QQvykLbzc2qYwOQq3HwUuB1DkOShvS1JtIUFgp8Pv/9u7j/f9qsHzmyKjax3746sMKIBhj4kjQ/E1r7e886sk0cKrHZwzu+U/gF9cl+EG0THtwsL7N/BgmncWaJ1RopGtdxVg6i+M6BNUnsNzXfqIpEJiYvP9m82Dvu29Emh3rDZxEguDeB1tF28PQMvXSedYG+qlffmPvxyezUXN0vZMxtaf4PCp4goxwzo8xPH2K1fZUoI/ax1cryT4MncYdHMOfu7m7/zM7CPUKyzdu9TYdCMBxcV0XhkYIwyFaMzcIamUdQXuS2MMZHBSVe/0acNVa+48O8hgNnOrxjF8i0T+AiRnCXD9hbUR++feSORYWsaWVJzcw49QrJNKpraXKsEhJuGhNzkPRXpnFGFktpXtu4X4am5x7YQrPcymV6jxo1GH+BubUSyzObxD4Acl0kvpeIz5TOZi4CGGIcRxsYUWatLsHPjmuDMbpLhDf3s9ECobPEs9mMMaQ6B8gzH0Kf31JVh45SOD30jKwprDQOa6thvx4cOJSPWiPfk938jx+rbZ7oNOjuAkYOkO1UpUC99H0J+2bPHkO6WP8w8DfBD4wxrwbbfu71tpv7vcADZzq8bTqNOfvyRd5IinZRK9NZk9qHcu2zQLN0qpkPpsbkpHFk/tXk9mpnT22s6G+EXmO/YJJGHSm1vSPkzp9ntqHDx1TIOoVbl+dYWh8iJUbN6UPL2hhl+9TDM/A7BXZ770Cw9h5HMchWLiF3WsaRdyDqZdIpJK0BsawDz6EVoP8Z/4Cg8N9eJ7L7P0VKg/uUXJipHNp4ok4qekxZm7PY29/b//99tLQN0p8aFRqvu5krRznfc4D/96H+y8B9zB+ExauY3PDMDkNS3efzMLhSgHW2j8Femr20sCpHk+7f6u4dLT70a29T49bgs7EOsEq7slUhYNmYU4cN+52KgU9jLWEt95hZS4ngSSRgsCVrDGZk/skUjB2HjaiAFVckh8nJkYwd3P/fWo1YOEmzXQ/5AZlvmTcw0t5fP7jEywXa9TrLSpLadLZNIEf4LgOxfUKxhjsyLS8VjtjbDejJ1KMvv4XGBhMc/2tD2U0avvn/0Hnxn6U+ZJhAJkBwjB88j+41DNDKwep54MxcPbj0KhJpnZYX2p9I/Lc1aIEkDB48st4dfe7+c3tJfLao4TDQP7eGQTWZig/ooD4Nm5CAnOrAf1juHEXv1qVDLfVkCk0bkICaLUIo+dx84P4C3f2norTFnMkiBWXIJXDmzxDo9pg+dY9vlHc5Ac+d44wtLAxz4Ybh9UZyhMXiCfiOK5DfGScRsyR91haIX7mRYyRovIXLg6xvLwpxyX6Ahv6oR+nL5/uLAr+qEFKj8vEwG/Smr+rBQtOsGMSNzVwqo/IWvqG+hmfOs+Nt8JotZTNrVVRAPkSTqTkrH+c4uhxD3diWvrHFpvSFPuw4PEkOHEpvNAu7ZbK4aU8mo0msViMYHVeitgnczLYppdPdCoHfaN4wzumhmTS+Msb8ne08ogZnsK2GphsXorMB/7DXyudJ//yxynOzJKfnMBLedS8Gsn0CBvLG/x/33pPikmMX5QMt1XHtprEc2lisRhu3CXbf45z54c5PZrlc9N5/uz2BqWq9CvO3d/+Xtdu36U8ONz5f80MyNJiT/obLgyiPmQtVHBSaa1a9VwpXXsPaz8uASGdw0t6NO5f6/RDxZO4068Si8Vo1qPSaKWVgwfRVgO/sC4BeOolCWaPqnvabmLMDmJGz2KMIZy5erBsxXElaLpxTGoEx3HAyOLVdu4G5uzLeBNnsPY0xhga+RGob0o5vIMIA9z8oHwTdPWsOI6D397vmANuAltcBcfFrs7RDP3dPxiSWZi4JFWDjMGdukDx1g0IWhRvVTBDk7z8iXOcGs/y5282qITh1hJpoTdOK5mD8jqbxRWGP/Zxzpwd4OWzg/zI+Tx9iThV36dUbXHtxgprC2v4t7+//fWXbtNsj6iGwy02MfUSLN7WRaZPqkOqVfs4NHCqj26zQPmDN2VAS9CiMXJOmtaGz0h/4+DUVvDx0h6kJ2lY+8i1EbdZm5GsI5HaXtwcZNCK35RtiZQEk1yUBQUt+of7ASjyMuGddyWo7sxcuivwWCvNtjEHO/ki/v33ZHvci2rKAgYp8gDkR4fxUlMsH3RFjnR/1FcHTnzHYJn+UTlmgS+jVs++hpfywEIQBPjNM1JJafW+vL+JS0y/eIrZhEsqk+LMuREWF/Ks3b4Ly3ew5VWuWMvi6CCJZILXLp5nfq5IaaNEa3mO/PlL+P4Am7MPsNZyYTJPf8rlw6Uqfmh5sFblxq01lr77xpNvGu/VzIfH55tTHQnNONXzpb2QdMyBlbsyST2Zk2BUWpFAme4j7sWlydGY3srFtftOd65MEvdg6kVisRjh/C1IpBh8+VUcx2FwKMPC7DobswsyMT+RkseMnt9edMFxpdLQxrzsd7UY9T2OdwqIt/ch5uz68DbqDerVuryfAzQlehNnaDaa+DPX8L1M5zjEHNzhSfzTr8qyWY0qVAqQGgMDjuvguA6NQjR4x1pS2RSffnWMH/2B01yfL+K5DkMDKWoXhnnv7WHyQ3niCZdW02dtcY3CSgF/eYbkqXPSh2rgtU+e4U5fhrGxHPVWwM3FMt9/f4GNlQ38m987PnMlj8M+KIUGTvUkJbMwOh2VjGvJl7+1EogyA7A2S6uxCZMvkhiZoJkdgLlrj2527S7iHXPkddrZjxOHapn40DD+1CWCB1dIpjw8z6VQqEmQc+MSmNpBt2twy9Z+xxxGPvcjVMtVNj98S7bvNWI4DAg2yzhu/9Z8UWstjXrjwF/sjSt/3rmyY3SsX16V/fGbkMzi9g9v3dZqtAiX7kXrfFpwXGozd/izt/t4/ZOTzM+X8f2QH/7UJGnP5exojrgb43vXlllerDM4NkilWMEfnKB+/wZkBrChJe7GOHdukMmhDNfvbXDlO2/11hpwUMbIedDd/60jZFUPjstvJw2c6skJA1idkZJuq/e39ydW1mR7vQKNKibt4WVSNPJj20u37WQMjF+SgTjrcxI4+0bktXLDxIdGJYMFYk6MwFrm378s1WsaDfBbnapGbd2fPjeBe+oSo1MjxOMOK/fmOl/m+31K1+dors/BxCUSyaioe2H5yUzIb2fuAE6w9d5ajRbh6pxkwKPnSI9PMn5qiMLGJvl8isX1Kvm8x7nJPA0/ZCgbw6bi3F0uU6/7mJihUqyQzqU5e2Gc21dd/MV7lO43eaveoFFrQHlD3ntp+aMtPr2fiRcYvTiN58UpFzcp3Ln98P97pbax2lSrnkPNmlz26gurlaF2Xb6Qk2kaa6vRIJ/63s/luJJpxhy8vhyN9qAYx91qUnX6h7cCC0R9jpkBWJvBb/cTPkoY4K8tUUx6kkHutz/d2pnr/fdpxhxpLu6f2P/+7WbpXkeE2pBmtSbBv7gkhQViDt7wGBOnh3GcGIEfMDaaodkKCAJLNunS8kMCCyulOvOLFUqFTS69OMZbf3KNwdEBXjo3yPJikVriIgkvQTafodX08VsD5AcyrCxssLm4ILVtn9QXVdwjls6yfPmy9H/b8Oj7TNUzRUfVqpPLxGDlgQSBh30InLj0SdbKNDarkrGCBNPABzdB8OAKwcQlGdFqrWR8G/MSqA6a/RkDmX5azRaOGw3UGZyS53nUhzQMugYUPaTwgbWdfsxePvitBjz4YNdj/JZPsVBl9b5MifmT2XlIJJl++SwXhqUf992ZEq5jGBpKUa02mRrOcO6l06yvlphdqdA/mGXq9CBzM+vk80kSCXnvS0sVNpcWpcXgSX5J+U3C+5cf3Syv1H50VK06scLgYAtKt7NXkBJ07QDVnuLSzvruvdd5zMQLvPRTf43F+Q1iToz1W7dldGw7e40G1GzjJoinkjQ31uibmmLg1RdZmV+l1Uv1Ib/58CkSmX7cU5fw5+/0vph2e38dV9Y/XZshlU0Ri8XIjQ5RKXqMnxmnsFpgdDTLi0NpAgtXFjYZziYY6UsyN1fiO//5LqfPDuI4/VSrLay1FAo1cvk0tVoL141RqTRZvL8o1ZLSednXXhaJftT70KCpPiLNONXzIZmVggdP+oQ2plOofK9RnbFoEv/63FYQff2Ln+Frf+PT/O9/cIPf/ca7kMrKfTY3ZOHjhRuycku3RIqRyWEKiThjkwMUCzVaN7/X+/uJJ/eumhMFHn/mhhynx5UZYPDiBbh4gXPnh6nXfSqVBOWbH7IQBKRyGRYXy/zGW3NMDcr6mDPrVV6b6mNqqo83rj9g/eplaFZlX/tGSQ/2E0/EaTaazNSbskrL6VfwzlySBbWdOO7oafyVOWmObtfx1eLq6oTTwKk+mp3TQ56UgUlZgaNahNTA7hGuQ6fx+gdopbOyLmizxq0by/wff3iT77z5QPorF25KtmRisp/txZzjSam16jeJjU1TWC1irWVkJEOpVN8dNPdapWSnvaraxBxZRSaXpdVsYcsbUiTBcXtvtgVKayW++BMfo1RrUq/79PV5DL76Sdbf+VNSn/kCAwMp7s0WWSvWCcOQ2Zkif/BHV0l4CZi/3sn4amWIJ6k+WJNjY4z8AAEZ5JQ8LfNtU9PEnBhTr3+GTCaBMXDzyjDNqwcoZK/UITgmCacGTnUMdM1jJAwgOyjzGatVaQJ1XMlsvUznMbkBMBD34jT6RmD1Aetvfpt/c/eeZJrNuhQ4rxYk6LWfNz8o64fWShBzCCtFquuLpE9NM5BJcOZMP5XXfpDND97oZFbtxbbbzb17ZVx7Bc3+cXBdqZZUjErVnf8UX/7yp/j2d26w+e6fHvwYZQaIe3HuzRVJJl2mxrKkPZdazWc9kWL9+lXq1XP0DfZxdqqPdy8vsnx3FlbvU485u5tJ12bk/bT3deIF+TFhOoOt2lWNMpkEo8Np5hbK9A/3s3wYI26VegQLUmf5GNDAqZ6uREq+qNuZ6uCUBMSVe1K8vVmDoVP4q/Od8m0mBgPjeAMDneex0Kg1pKpOpl+mNZgYpLLEkx4mlcQMDNBYcqERTXnY3MAvbUC6r7MSyso9GD3P658/z4PFMgvzRSlJ1724diIJky9AaVW2H2S1lDCQLLm41Mkuc0N8+cuf4n/+wnneeneezWT24Bl7YYGaE+fG5SapbIr7yQRYWJ9dlH7fzQLVyhrVdD+LDybJ9edg7sreK5jEHGn+DlqdEa7rszJid2xq1wJLYWjZKNap133y/SmWR6YPXl5QqSdFBwepE6t/XGrazl7pFDNor6pSWJTb5q7Jl3qqT/o5vTRurm/b0zQ2qzB/neD8xzGuizUGTr2Ml05t/+LvG5Lg1B7Fu3irk9nmx3BGT+PGXdY2avh+yOp3/zQqq9f1CS2t4g6P4/t52ceHaZf822uB7moRP7Cs1RtMTfWxPv0KrVJR9ulR/Ya1Msx+SBBzqLSnwMSTUlC9rdWAzQ1SZy6QyqYotwvrD07Jv6VorunYeSkKEQRy7ErrcrxzQ3sEzZDVlTLJVIIwtPh+yOiLL7C8+mB75SXQAuzq0OngIHXyuAnc/mGCIMBGU02Yvy63ZQclEAyMwey1zhJiQQtncKwzVQQpBkCtLMXfl2ZkBKiJSfOsyRK0AilR5ziy6HPfiAwQan/o2l/4xSUCwD11jkqlyfL8urzuzg+nDcGCm+3DHz0n/YH7jTRt1fdfrDnw+f1vvsub7zwgm0tJ32NfnlbsJZkzuTNgt3lpmZ7Tfn6/KZd25pvp75QKNAbHdViZW5Egnh+HeIK4FyfI9RMG4db7icVlKbHY6HjXe2VX8CwXyoRhhmTKY2FmhfpGoZOROy6cfQ2Kq4dTbUipLsckbmrgVE+JMZDKYWIG++DG7iZKY3BHp/AX7knd2PLq1mjd4O772OnXOqt61OuS0Y1fxMn1E1RKMP3xaCFppGLQ2ix+VOT9oZlQcYlGMsvdzZqsvLFX5pdI4cQdjDHExyaoxZO7KyO1Wfvw5tc732f9vsvkf/VTJJMjzN1fkYpHmY9hQ0uwsdyppuO4kB8jMXZKBjtZaFRrGMfBrs51Bky1m4JTfRD3qFfrsvRZblgCbWmZVlRnd2u9zzAgdFzCsQsk8hJQ/dKGZObtvmRjoG8IL+lRLVcprZfw525v7+eNObB8X36YKHWotHKQOknSeRg7TyKZwG/526vzGANnXsNNevibVcma6hUZuJLMSqZVK22tRALg5XMEmZexYfRB8lJ4KY9GtYG/fBMGpuR59gqC7Wyw+7Z2ST4vvff+r9yjYWJ4g4NYa/H684S5j9G684EEkXT+0XM52wth1yuQzLKyUiGfT8kxqfhbGXWQ6ZfnKa3Imp0TpwHwmz5BqwWLt7GDU1Dv6metlTvZZzKLX1yCdD+xgVHCVkt+PNQrkmn2T0B2ACobcr28SrOyLre397+y3jlWWelXttbKsmepvs7tsL1EoFInhAZOdbiMgdFpvLQn9VYfXNmVAbabDJ18joZ/Rvr8mjUJZMOjJCbObhVUb/OLG9uaTBuDU5Kh1iudvrxufSMS4JIZKT6/c4CP40q5vuwQ5AYhFoO1ecnqAh+W79AIAxKDw5iYkVJ/4xdgY4GxT/8AS5c/fHhTZWZAAuf4BdxEgsJqgUqhgu9vPxau5+G353uGgcynBJla054v227e3ks7262VCddnd7dt1cvRAKzJreXKGrWGHPPhMxIs283V1QKsz9NYbsr+twcV6ahadQSk5N5R74XQwKkOl7Wwcp/GCpKZdGdlUW3ZrXqzlu2354bx+vN7Pm1iYIimiUnfYOB3mjeNkQwrO7g9M3ITeEPDNIrlvUfF5oZxRya31g3FQqM7Mw58WbTZTeDl8wR+QLIvx9/9xS/yL3//ppToi3u7s6+4J0Ennd8+KhhoNVs0G00JYFuHy3b6T6sFuZiYrC/6sCbgREp+HGQHolKAFtYe7N4fE5NqSQlp9o45MSlZmMyBtcQHR7b+P8LgFGEQEtx9X+vKqqNndXCQOknaZfKgsyZm3wgkUlz89MsUNjYpF8o0SmUJQG1dg2walaoUO2/WJLDEPQmWqRwMnYJaRaZUhAGMnJWA3R04y6s02o/dS2Udv7yKnx8je/Y8laXlvQPV2gyNjXmoV/Czg3zn+iWW59dkX0+/Avc/2J5Rtwc5DewuAt8sSlBsRN8FXtrDjbvSN1lYhKHTxLJ5wrvvPTybjXswcQkvu72pueldwo2725q525l7LBYjDENSmRS1mRIELRKnLtJsNKHZxE2n8TdWogXCW/u/tlJPkQZOdTIls8QnzxGLxTAxw8pSQQoEAIT+9nmHyQx+05cv/nblm5RkRmxEQc1vYuIJ3PQoYf8wQaNOIpUkDEP8kWmZ6gKP7otrB9R4UrJOL713H2k7mDouuYsvc+3mGvXZu/L4e+/vfkxU4N1Le9s2B340hzSe3MqyG+m8ZMrtEbTLdwnDMzIlp13ZZ6dMP87kRdyEfJRtNEG8Wa5AvYwdGtuW0QetAH91nvSpaU6fG2d9rUzm7EU2H9yR/ufoPfi9LDKu1NNyPOKmBk71lNmQWCwGRn491jelPFxreW57WT1rYfYqQdyTOrOtBjgu7umX8CslGTyTG4ZaCXv3PVoxB8YubDXtOjEHPzvQCZwH4aVx+wcp3rzaWY1lL4kUDJ/hSz/6IkuFGnPf32PA0bbnzZDNZylvSBNxGIa4cRc/3b/9PZdXMSOnse0MPQwk05x6STLa9bldr3H6h77A9Nl+3vhP16VoRLUoGWhUTN4PA4KBMWwQQGFJRisHPtVmjVt3ZdqPc+pFcBOEhZXO82vQVMeQZpzqZDFGVvfoG942T7CxtiJTGbqbVdtsCGPncbwkgeNCzCHmxPAGBrD5fpoba9J/2NiUfsul2zQSr25ld4lkgubECzJ1JJ6Ui5eW16qVdo80aNbwN6LqQA+bwjJ8houffpl6K+DqteWoHm5UlL5ekTq7mxsSpKJBNmt37kmwH5gEN05+KE/yxVcofHe5sx/WYu9f3h60Ah9W7suPhL6RXQUYZj68ydrSMP7d9zsZdXcT88Y8trC4OxC2j7cxBLe/r4FSqR5o4FSHb2Qaghamf4REMrG1OQxCyexMDEbPSZDo7suzFvwWbi6LOzRMGIQ0H1yHkbMkUkl53NJtqUbkxKNCCNJcaYzBxAzewCDNdBYbBHgZ6V8NB0do3Xp3e3+nm4Azr+K6rvTtJVL7Ty+plVle2KBYqLIyvyKBO9Un+9s3wpf/2y+wWqrzxjf+RN5PO3OMJ/Hyso9hIP2LBSe+fT/26oOt7TOgCWD+OtWHjbK1dveKML3crtQxohmnOhncBOQGpIJPq0U2nyXwAwI/oFqO5m1OvCBZZGhp1ssymGjsgkwdcToDhKTqjSXmujTXlqVeqpsg1j9CWI7LNJTCIk1jonJ+OWhs4p56AcdL0FhflzmfGWlqjZ5UmnMTUsnHxAy2f1iaiPcLSMUlSrMpTDYPm6WuQghFyI+xWqpz7cpipyhA4MslDGg2miSSia1iDvtWGVJKbWOtFkBQJ4XfhLvvErgJ8DL4wwNk+tIYHpfHggAAE0FJREFUY8jmsyz4n2RwbJDN8qZ8KJI5aVKtrOHmB3HincASLM+Al5EVUboGEYWtlqx64rgSTI2RjHFjHpJZHNeRwTjLd+DUKwAyNaXakLq445fwBga3pqEEi3f3bjpuCwNYuo1d93bf1qrz9hu3SGVSxKdfoXXjnU4zaLTySK4/x/L9BVkx5WFFEw4i7kWZ+RNYbFqpY04Dpzo5rJX+NxOjuLhCvZoDA42NDVibYf1BE9wE5tRLJEYnac7dgbyMBm0UyxBPyMjQWkmKAhTyssIJyPM++AB/+KwEy7ELMi0lNyyvmx8DgwwoijnEU0kJjkEgq53EvU7QJBrt2m4WTeelMEI00IbsoDTJllelH3GvUbrVIq0b79CafIlkfx+t9lJkAH4Te/c9VivT0nTruFJntlbe3ceYG5Lnf9TqKVq1R50kxyNuauBUT1E6D+uzNNxp4kkPktGUj6gp085eozl+UbI9Jw75fikA0KwRuglp1mzWJHDlhiVQtgujr9yTQDR4KuoTvC6ZZxjQaGxK9mktrbVlaUJt1aXowdSlraBpQ4s/cy2aCpMlcfYlaVbeLMh++k15jUc1r4YBzH5IfWmPgghhIIHdcSE7hDc2hd/yOzVqHRf6RqVu78yNQ/lvUOpZpRmnOnmKSxLk7r9PK57slKFLZiUgDkxKyblECupljDF441KrFQONtVUpDQdgYsTOvEI4d6OTIXYPLmqXh6sWt1e96Z6esj6HX6/gTL8MJioO4MSlXzYzgIkZKUiQyslrZAYkuB20WbQ7aCZSMHGJ3FA/m6VNwpmrUCuRzFwg4SVYKW1sBezs9AUCP8D3MlqxR6ljSAOnenq6pl3gNyGdx81O4RfWwcSkjF6juTUdpFFrQLUEqSxeJiVVdBIpqWSz9oCwvSpItzCQQDw4JYG6u2rRTjEHYg5+yyfmxGRtzGYNJi6BDWnWpRxe/PQLkqkmUpItHpSJAvHwGfKnTpHNZ2g2WoRByObQaRL5frykh+M6jL10ida5aaqVKvFEHC/lUduv6LxSJ5RmnOpkMwY3m5UVNxqbkJPpJjSi+rCJFDyIytcNnYbMlAwKinsPD4Yg/YIPm6IRvT6j5/CGhgn8gNb8XWm+zQ5KsF65h82PMfb511mcXcEbGadRb0hQ3m9qyE7DZ0kMjzE6NUwsZliaXaWxsY43MMj5T75IPO5QKNQorBZI59JkcilK6yXKzbKslnJMviSUOha0Vq068QIfv7COXy1In2YyKwN25q91AkYiBaPnZfBOW7scXVtuqDN4p1epnAw+amxKv+fEi6QG8tTu35S+xlSOViugMXtXstxWfe81OPdTWiZ97hzVSp1qpSrZ9Op9GmEA06Mkky59fUlK6yVK6yXWrl2VLNmJ4z+qCINSJ9HxiJsaONURaq+DCZDK4cZdWkOnJXhEJfYS/QPbKg0xfKpTkceYx5+Gkc5LH+jcFWmyTaSIJZNk+jLU0nnpz1yfZf4KEliTWVnLcu7KwTNBE6Nw/0HnfSZSEgzrFcqlOq1WyGalJiUHi4VOPVoNmErtYnUha6V2qG8SSw/ijU3SSPfB8j0Yv7i1moff9AmW7kcLMnf1lR602bTNcSUwxRxYn5d/ByZxBkYJw5ByoYw3KGUBg748juvQCHzp2/Qycv/uwNYe7Tt+Acrr21d3cVyZQtPWrEHcwx07Q2mjRCmqj2BDK9ls34jc51FTUJQ6oTRwKtWtsUkY9Est2lyWhj+1rTyfG3cJBielbuteRQP6RqQ5tXsUajvItVf6SGZhdFqKnQct2FyTAgJAMHMNklkaAxNSgH2zQqJP5op6o+PSt7pwY3e2OX6pU91o5/JbmwXJMt1EtEZmAF5mW1GHRqUqfarWSoGEVJ8GTqX2oYFTqW7rc7TqFUikcMdO4+1sojXAxuL+0zPqlU5wtFayQ5Bm1slLEAQk0ilMzBCm07TmouZTE5PH1StyKS7JNBC/STOcxqSzJJIJ3FwffjEvwbA91cVa6ZM1Rv7uLmLQ3pYdhHQ/sWSSsLgGmR0Lc1ekCIQ0BY9HhSKMDgxS6hiLHfUOKLWlWoTCIv7SDI3NKliZuxn4gTRn1stb61tKQfcosrb7OsNABu84rgSs4pI0gQJeNr3V7Nuq1TsDipq17SuOWCsBNOZAeRV7732a9SaO40ixhFSO1Kufw7n4GclW2wUcdlb+8TKyD/lR3HSacLMs25bv0Zi9S2NlSUYRZwfk/vWK9IUWFjrNv7khGL8o71Wpk8526tX2ennSNONUx09hESrrNDIDUF7Fd+IQT0hQNAb34qeljN7ynAwUygzgjp0hDEPCamWr+RWQoLa+AJPTWyunsHJ/++vttaRWs7bVJGzvvU8j1SdNsX2j/MhffIFitckb60vbV3PpVq9IFtlq4s9ekf3YGikrWXOrVtreX2ptpyB8+zms1SIISrUdk4YYDZzqePKbnUWeF292tueGmZweI5WKY188xY3378CDy9H0jZYEm+Gzkrm1lVZoDcs80GZhY/uAnYMI/K71K2O88948hdXC9tfYS2MTiov715OtrEvpwHaxhp2jaVsNrUWrVJfj0sepTbXq+OvOwqpF1pYK1Os+g4MpLn7snGSYpRVI98PpVyRT7B5t67hSUACi/smP8OGrFln5sz+gtTz36Gkj1j56jmlumMTYKQmgSql9taejaFOtUr3ym2xe/nM2F8+yNDdKvVrv9G/WK7JqSveUEIAwIAxCYk4MN5uV5cc+ylxJaztzLvfjuDKi9mHTZayFuas0C/3g7rFEmVJqm8PKOI0xvw78VWDZWvuxR91fM0717AkDyPRRr1Th3vudIFhalvU1d/ZZhgFBdRMs+CvzT7bAQMzpDFLqlsxurb/5SJuFTrO0Uuoo/AbwpYPeWTNO9eyxVgJme0oIdKZw7DXH01pYuEFjI9t7/+aj9I9L8Fx9sH37o+rpKqV6Yzm0wUHW2v9kjJk+6P01cKpn086ssm9U+jp3BrDu+x/G6NTK+qPX51RKPRHHZXCQBk71fCguHU0A2yvDVUodio8QOIeNMW93Xf+qtfarj/tkGjjV82Ov+ZjQmdepUzuUeqZ9hMC5aq397JPaDw2c6vmXzPa2HJhS6tg5TqujHOqoWmPMLxpjPjTGXDbGfN0YkzTGnDPGvGmMuWWM+dfGmER0Xy+6fiu6ffow902dIOU1bVJVSu3LGPN14A3gRWPMrDHmKw+7/6EFTmPMFPALwGejeTEO8DPA3wd+1Vp7EdgA2jv4FWAj2v6r0f2UUkqpzqjax7k86qmt/Vlr7YS1Nm6tPWWt/drD7n/Y8zhdIGWMcYE0sAD8KPDb0e3/F/BT0d8/GV0nuv3HjNlrgpxSSqmT6LhUDjq0wGmtnQP+IfAACZhF4B2gYK1tz0CfBaaiv6eAmeix7UrYQ4e1f0oppZ4tz33gNMYMIFnkOWASyNBDZYaHPO/PG2PeNsa8bX3tt1JKqZPiuQ+cwI8Dd621K9baFvA7wA8D/VHTLcApoF30cw44DRDdngd2Vci21n7VWvtZa+1njZs6xN1XSil1fByfIu+HGTgfAJ83xqSjvsofA64Afwz8dHSfnwP+Q/T3N6LrRLd/2x6XscdKKaVU5NDmcVpr3zTG/DbwPcAHvg98Ffh/gN8yxvy9aFt79NLXgH9ljLkFrCMjcJVSSqlDrVXbq0MtgGCt/RXgV3ZsvgO8vsd968B/c5j7o5RS6tlk0Vq1SimlVE80cCqllFK90MCplFJKHZC1sgbvMaCBUyml1LPhmGSch11yTymllHquaMaplFLq2aBNtUoppVQPjklTrQZOpZRSzwAdHKSUUkodnEUDp1JKKdWTY9JUq6NqlVJKqR5oxqmUUuoZoH2cSimlVG+OSVOtBk6llFLPBs04lVJKqQOyVjNOpZRSqifHJOPUUbVKKaVUDzTjVEop9WzQplqllFLqoHQ6ilJKKdUbzTiVUkqpA9JatUoppVSPQs04lVJKqQM6Pn2cOh1FKaWU6oFmnEoppZ4NxyTj1MCplFLq+LPoqFqllFLq4I5PH6cGTqWUUs8GzTiVUkqpHhyTjFNH1SqllFI90IxTKaXUM0DX41RKKaUOTkvuKaWUUj3SjFMppZQ6KJ2OopRSSvXmmGScOqpWKaXUiWaM+ZIx5rox5pYx5pcedX/NOJVSSj0bDqGp1hjjAL8GfBGYBb5rjPmGtfbKfo/RjFMppdTx165V+ziXh3sduGWtvWOtbQK/Bfzkwx6gGadSSqlnwKENDpoCZrquzwKfe9gDNHAqpZQ69mxt5Vv1d39t+DEfnjTGvN11/avW2q8+7r5o4FRKKXXsWWu/dEhPPQec7rp+Ktq2L+3jVEopdZJ9F7hkjDlnjEkAPwN842EP0IxTKaXUiWWt9Y0x/xPwLcABft1a++HDHqOBUyml1Ilmrf0m8M2D3l+bapVSSqkeaOBUSimleqCBUymllOqBBk6llFKqBxo4lVJKqR5o4FRKKaV6oIFTKaWU6oEGTqWUUqoHGjiVUkqpHmjgVEoppXqggVMppZTqgQZOpZRSqgcaOJVSSqkeGGvtUe/DYzPGlIHrR70fx9AwsHrUO3HM6DHZmx6Xvelx2e0gx+SstXbkaezMUXrWlxW7bq397FHvxHFjjHlbj8t2ekz2psdlb3pcdtNj0qFNtUoppVQPNHAqpZRSPXjWA+dXj3oHjik9LrvpMdmbHpe96XHZTY9J5JkeHKSUUko9bc96xqmUUko9Vc9s4DTGfMkYc90Yc8sY80tHvT9PizHmtDHmj40xV4wxHxpj/la0fdAY8wfGmJvRvwPRdmOM+afRcXrfGPPpo30Hh8sY4xhjvm+M+b3o+jljzJvR+//XxphEtN2Lrt+Kbp8+yv0+LMaYfmPMbxtjrhljrhpjflDPFTDG/GL0+blsjPm6MSZ5Es8VY8yvG2OWjTGXu7b1fH4YY34uuv9NY8zPHcV7eZqeycBpjHGAXwP+MvAK8LPGmFeOdq+eGh/4X6y1rwCfB/7H6L3/EvBH1tpLwB9F10GO0aXo8vPAP3v6u/xU/S3gatf1vw/8qrX2IrABfCXa/hVgI9r+q9H9nkf/BPh/rbUvAZ9Ajs2JPleMMVPALwCftdZ+DHCAn+Fkniu/AXxpx7aezg9jzCDwK8DngNeBX2kH2+eWtfaZuwA/CHyr6/ovA7981Pt1RMfiPwBfRApBTETbJpA5rgD/HPjZrvtv3e95uwCnkA/6jwK/Bxhkwra787wBvgX8YPS3G93PHPV7eMLHIw/c3fm+Tvq5AkwBM8Bg9H//e8BfOqnnCjANXH7c8wP4WeCfd23fdr/n8fJMZpx0Tvy22WjbiRI1GX0KeBMYs9YuRDctAmPR3yfpWP1j4O8AYXR9CChYa/3oevd73zou0e3F6P7Pk3PACvB/Rs3X/9IYk+GEnyvW2jngHwIPgAXk//4dTva50q3X8+NEnDfdntXAeeIZY7LAvwP+trW21H2blZ99J2q4tDHmrwLL1tp3jnpfjhEX+DTwz6y1nwI26TS7ASf2XBkAfhL5YTEJZNjdXKk4mefHQTyrgXMOON11/VS07UQwxsSRoPmb1trfiTYvGWMmotsngOVo+0k5Vj8M/JfGmHvAbyHNtf8E6DfGtEtLdr/3reMS3Z4H1p7mDj8Fs8CstfbN6PpvI4H0pJ8rPw7ctdauWGtbwO8g589JPle69Xp+nJTzZsuzGji/C1yKRsElkI79bxzxPj0VxhgDfA24aq39R103fQNoj2b7OaTvs739v4tGxH0eKHY1wzw3rLW/bK09Za2dRs6Hb1tr/wbwx8BPR3fbeVzax+uno/s/V7+srbWLwIwx5sVo048BVzjh5wrSRPt5Y0w6+jy1j8uJPVd26PX8+BbwE8aYgSib/4lo2/PrqDtZH/cC/BXgBnAb+F+Pen+e4vv+AtJ08j7wbnT5K0ifyx8BN4E/BAaj+xtkBPJt4ANkJOGRv49DPkb/BfB70d/ngbeAW8C/BbxoezK6fiu6/fxR7/chHYtPAm9H58u/Bwb0XLEA/xtwDbgM/CvAO4nnCvB1pJ+3hbRQfOVxzg/gf4iOzy3gvz/q93XYF60cpJRSSvXgWW2qVUoppY6EBk6llFKqBxo4lVJKqR5o4FRKKaV6oIFTKaWU6oEGTqWOOWPMdPfqFUqpo6WBUymllOqBBk6lDlGULV41xvyLaP3H/2iMSRljLhpj/tAY854x5nvGmAtRRZZ/EK0R+YEx5q8f9f4rpXbTwKnU4bsE/Jq19lWgAPzXwG9G2z4B/BBSveWvIZV+PoHUU/0H7ZqhSqnjQwOnUofvrrX23ejvd5BVOaastb8LYK2tW2urSDnFr1trA2vtEvAd4AeOZI+VUvvSwKnU4Wt0/R0A/Ue1I0qpj04Dp1JPXxmYNcb8FIAxxjPGpIE/Af66McYxxowAfxEpKq6UOkY0cCp1NP4m8AvGmPeB/wyMA7+LrGLyHvBt4O9YWRpMKXWM6OooSimlVA8041RKKaV6oIFTKaWU6oEGTqWUUqoHGjiVUkqpHmjgVEoppXqggVMppZTqgQZOpZRSqgcaOJVSSqke/P9Y8odtKqfBBQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot result\n", "fig, ax = plt.subplots(figsize=(8, 8))\n", "wat_occ = ax.imshow(landsat[0], cmap='Blues_r')\n", "fig.colorbar(wat_occ, ax=ax, anchor=(0, 0.3), shrink=0.7)\n", "\n", "ax.set_title('Water Occurrence')\n", "ax.set_xlabel('ncol')\n", "ax.set_ylabel('nrow')\n", "\n", "plt.show();" ] }, { "cell_type": "markdown", "metadata": { "id": "dRvSOCNaFxTu" }, "source": [ "*Export result as GeoTiff*" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "id": "1B8EyV_9Pr-1" }, "outputs": [], "source": [ "# Data dir\n", "data_dir = \"output\"\n", "\n", "# Output raster\n", "out_tif = os.path.join(data_dir, \"landsat_water_occurrence.tif\")\n", "\n", "# Copy the metadata\n", "out_meta = image_source[0].meta.copy()\n", "out_meta\n", "\n", "# Update the metadata\n", "out_meta.update({'driver': 'GTiff',\n", " 'dtype': 'uint16',\n", " 'nodata': 0.0,\n", " 'width': landsat.shape[2],\n", " 'height': landsat.shape[1],\n", " 'crs': image_source[0].crs,\n", " 'count':1,\n", " 'transform': image_source[0].transform\n", " })\n", "\n", "with rasterio.open(out_tif, \"w\", **out_meta) as dest:\n", " dest.write(landsat.astype('uint16'))" ] }, { "cell_type": "markdown", "metadata": { "id": "tToFTu4DVGrp" }, "source": [ "**Create water occurrence for Sentinel images**" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "id": "DIsW5HayVHbE" }, "outputs": [], "source": [ "image_binary_sentinel = []\n", "\n", "for i in range(6, len(image_source)):\n", " # read all image into numpy array\n", " img = image_source[i].read()\n", "\n", " # make binary image, 0=no water, 1= water\n", " img_binary = np.where(img == 10 , 0, img)\n", " img_binary = np.where(img_binary != 0, 1, img_binary)\n", "\n", " # Append to list\n", " image_binary_sentinel.append(img_binary)\n", "\n", "# Sum array to create water occurrence\n", "sentinel = sum(image_binary_sentinel)" ] }, { "cell_type": "markdown", "metadata": { "id": "YBuCVBcKGC0X" }, "source": [ "*Visualize result*" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 419 }, "id": "SAKhGLi2F_xv", "outputId": "5d6e1d20-e71b-4430-d232-ba75c4497953" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdQAAAGSCAYAAABaJjAdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd5yb1Zno8d+RNNJI03vzuPeCbVzAgOkYDAQICQmpJGxCSNvdm83uJvfu3WRT9ia7m2xLNglJSAgpJAEChGpT3HDBvbcZe3rvI41GGknn/nE0zNiePtKMNH6+fPTxzKtXr45krEfnnOc8R2mtEUIIIcT4WCa7AUIIIcRUIAFVCCGEiAAJqEIIIUQESEAVQgghIkACqhBCCBEBElCFEEKICJCAKoQQQkSABFQxJSilvqqUeuWiY2cHOfbgMNe6USlVFYU2piulfqSUqlNKdSmljiqlPhnp5xFCTA4JqGKq2AZco5SyAiilCoAEYOVFx+aGz40apZRtgGN24HVgBrAOSAP+FviOUupL0WzPCNp2yTEhxOhJQBVTxV5MAF0R/n098BZw+qJjpVrrGqXUJ5VSJ5VSnUqpc0qpzwAopZKAV4BCpZQ7fCtUSlmUUl9RSpUqpZqVUn9QSmWGHzNTKaWVUn+hlKoA3hygfR8DpgMPaK3Pa617tNavAn8JfEMplRq+VrFS6lmlVGP4eX7QewGl1Kf7tfmEUurK8HGtlJrb77xfKqW+Ff75RqVUlVLq75VSdcAvlFJfV0o9rZT6tVKqA/iEUipNKfVzpVStUqpaKfWtfl9EPqGU2qGU+jelVKtS6rxSamO/58tUSv1CKVUTvv+5fvfdrZQ6pJRqU0rtVEpdMZa/XCHigQRUMSVorf3AHuD68KHrge3AjouO9fZOG4C7gVTgk8C/K6Wu1Fp7gI1AjdY6OXyrAb4I3AfcABQCrcAPL2rGDcAi4PYBmngb8Er4+v09AyQC68IB7EWgHJgJFAFPASilHgC+Dnw83OZ7gOYRvDUA+UAmpnf8SPjYvcDTQDrwG+CXQADTg18JbAA+1e8aV2G+nGQD/wL8XCmlwvc9CbiAJUAu8O/hNq8EHgc+A2QBPwFeUEo5RthuIeKKBFQxlWylL3iuxwTU7Rcd2wqgtX5Ja12qja3ApvD9g3kU+D9a6yqttQ8T3N5/0XDp17XWHq21d4DHZwO1Fx/UWgeApvD9azHB+m/D1+nWWu8In/op4F+01nvDbS7RWpcP/Xa8KwR8TWvt69e2XVrr57TWIUyAvhP46/DzNmCCYv+55nKt9U+11kHgCaAAyAsPo28EHtVat4Z73lvDj3kE+InWeo/WOqi1fgLwAVePsN1CxBWZOxFTyTbg8+Gh2Byt9VmlVD3wRPjY0vA5hIcsvwbMx3yxdAFHh7j2DOBPSqlQv2NBIK/f75VDPL4JE4QuEA7I2eH7izGBKzDA44uB0iGuP5RGrXX3Rcf6t3UGZri8tq/TieWic+p6f9Bad4XPS8b0fFu01q0DPO8M4CGl1Bf7HbNjvjQIMeVID1VMJbswyT6fBt4G0Fp3ADXhYzVa6/PhIcdngH8D8rTW6cDLQG80GWgLpkpgo9Y6vd8tUWtd3e+cobZueh3YGJ6j7e99mF7b7vBzTB8kSagSmDPItbswXwh65V90/0Dt6n+sMtyG7H6vLVVrvWSQ57u4XZlKqfRB7vv2Re+ZS2v9uxFcV4i4IwFVTBnh4cx9wJcwQ729doSP9c6f2gEH0AgEwr3VDf3OrweylFJp/Y79GPi2UmoGgFIqRyl17yia9yRQBfwxnMSUoJS6HfgvzFBxO/AOZlj4O0qpJKVUolLq2vDjfwZ8WSm1Shlze9sCHAI+rJSyKqXuwMzljpjWuhYz5P09pVRqOAFrjlJq2OuEH/sK8D9KqYzw6+odYv8p8KhS6qpwm5OUUncppVJG0z4h4oUEVDHVbMUkxuzod2x7+Ng2AK11Jya79g+Y5KIPAy/0nqy1PgX8DjgXzk4tBP4zfM4mpVQnpkd51UgbFZ53vRXTa9sDdADfx8zL/mv4nCDwHkxiUAUmAH8wfN8fgW8DvwU6gecww60AfxV+XBvwkfB9o/VxzBeNE5j35GkGGKIexMeAHuAUJtnrr8Nt3ocZGfhB+JolwCfG0DYh4oKSDcaFEEKI8ZMeqhBCCBEBcRNQlVJ3KKVOK6VKlFJfmez2CCGEEP3FxZBveMH7Gczi+CpMVZwPaa1PTGrDhBBCiLB46aGuBUq01ufCFXGewlR6EUIIIWJCvATUIi5cZF4VPiaEEELEhClTKUkp9Qi9dUqVbZVKzJjcBgkhhHiX9jY2aa1zJrsd0RQvAbUaU3qt17TwsXdprR8DHgOwuHK1Y8EHJq51QgghhtR96IcjrT0dt+JlyHcvME8pNUuZfSUfpN9CfCGEEGKyxUUPVWsdUEp9AXgNsAKPa62PT3KzhBBCiHfFRUAF0Fq/jClgLoQQQsSceBnyFUIIIWKaBFQhhBAiAiSgCiGEEBEgAVUIIYSIAAmoQgghRARIQBVCCCEiQAKqEEIIEQESUIUQQogIkIAqhBBCREDcVEoSQgghJppSqgzoBIJAQGu9erBzJaAKIYQQQ7tJa9003Eky5CuEEEJEgARUIYQQYnAa2KSU2q+UemSoE2XIVwghRNzacPsdurlp2NHYAR04sP840N3v0GNa68cuOu06rXW1UioX2KyUOqW13jbQ9SSgCiGEiFvNTU28vWffmB7rTFDdQyUZAWitq8N/Niil/gSsBQYMqDLkK4QQIm5pQOux3YajlEpSSqX0/gxsAI4Ndr70UIUQQsQ1zQii49jkAX9SSoGJl7/VWr862MkSUIUQQsS3KMVTrfU5YPlIz5eAKoQQIq5FrX86SjKHKoQQQkSA9FCFEELEtZEkGE0ECahCCCHiWhSTkkZFAqoQQoj4FhvxVAKqEEKI+KWJmXgqAVUIIUQcG2GRhokgWb5CCCFEBEgPVQghRFyTpCQhhBAiAmTIVwghhJhCpIcqhBAibvXuNhMLJKAKIYSIYzpm5lBlyFcIIYSIAOmhCiGEiGsy5CuEEEJEQIzEUwmoQggh4lyMRFQJqEIIIeJarCQlSUAVQggRt3QM1fKVgCqEECKuxUg8lWUzQgghRCRID1UIIUR8i5EuqgRUIYQQcU2SkoQQQogIkKQkIYQQIgJiJJ5KQBVCCBG/Ymm3GcnyFUIIISJAeqhCCCHiXGx0USWgCiEuZHeCUuDrmuyWCDG8GKqUJEO+Qog+didzbttAypLVkJY32a0RYkT0GG+RJgFVCNEnNYc7riqms6EZx7RZMG0J2OyT3SohhtRbz3e0t0iblICqlCpTSh1VSh1SSu0LH8tUSm1WSp0N/5kRPq6UUv+llCpRSh1RSl05GW0WYkpzprDg3vfyzX94P28erIFgD77aSrAoEuaukN6qiGl6jP9F2mT2UG/SWq/QWq8O//4V4A2t9TzgjfDvABuBeeHbI8CPJrylQsQbqw1mryTrmlth9krILDJBMcEBFqu59Z5XMJ9Hv/IQv3x4LTtLWjl3pgbczdBaA+5Welqbmbt+HYlL10lvVYghxFJS0r3AjeGfnwC2AH8fPv4rrbUGdiul0pVSBVrr2klppRDxIDWXD39gLf/7prm0uP1owB8IUdreyfNHGjl8op70dCfrr8jnvQvzKO/0cPs3XqWj9DSLb76WMquVriNvg7cDlIWSbW+Tt3IV+bfeTtneQ9BcOdmvUIg+MZKUNFkBVQOblFIa+InW+jEgr1+QrAN6x5iKgP7/eqvCxySgCjEQpchduoyHr5zGm+ca+N5zpwgGQ7g7urjx2tlcNTudtbPSWZyTxIGaTh5+bDc1u7aDtxOAE9v3c8v7ruegK5GWvVshFMQ2Ywn177wNWdOZu24VJbsVNFVM8gsVInoJRmMxWQH1Oq11tVIqF9islDrV/06ttQ4H2xFTSj2CGRKGhOSINVSIuJOYzJceWIrdauEvv/4s1Je+e9dze7fwHICymOFegB7fhY9vLOON3/v5+GfvZU9eKqdffplATSlZK6+i+fA+SnZ1MHPNCsqOJkH1yQl7WUIM5rJeNqO1rg7/2QD8CVgL1CulCgDCfzaET68Givs9fFr42MXXfExrvVprvVrZnNFsvhCxx5li5kcB8mZz/fQsvvLSCWg4d+F5vemNoaAJpP2Dqc0OyZlmfrW1hl/95EWuvaKArDXroaud5qMHmXfzjZCURtk7B5h35QKSll9n1q0KMWnGmpI0BZKSlFJJSqmU3p+BDcAx4AXgofBpDwHPh39+Afh4ONv3aqBd5k+F6GObv4YbPnofM27egHXealZft4jmLj87X9k7uq/uSkF6QV9gbq7k8f96huycFBwL14CnjbNbtpOUmU7hiuXUVDRRUJzN4rvvxL7oKkjPN7e8OeZPISZKjCxEnYwh3zzgT0qp3uf/rdb6VaXUXuAPSqm/AMqBD4TPfxm4EygBuoBPTnyThYhRzhQ++7GrKUqzk3hlHnsrilg/K5UHv/s61J4Z3bV6fFB1/MJj7fWUnqziu391A999Ko2G7a/hObQDT2oOc9ZfQ25uEkrBuuvm0xOcSyAQwu8PcPZkLZ5QEDoaI/dahYhxEx5QtdbngOUDHG8GbhnguAY+PwFNEyIuzc9J5Cevn6ejo5vkZAe/eWoPnD8YsesHSg7w41eL+Ox9i/lOczu+E7uho5HSV16k1OEyy3EsVtOz7fZAwA8dDcNfWIgIiZEp1JhaNiOEGAmH64I50AZ3D0lJCRzbdTwySUKuNDNkW3Pa/B4Kcva11/h1aAPf/Oy1/GxzAWcOl4KnDdrqoK7EnKdU7GSHiMtHDNXylYAqxETqHwzHyu/t+znYQ0jDgqJ0jmRm4HWshO4uM9Tq7YBgYPTX7+m+dKjW76V0+06eznLxf9+7iIT3LaGy08uOkjY2vXESb8lRyJ8L1acubJ8QUWamQ2MjokpAFWIiJGfinLUQZ7KTlhPHoL1+7Nfq/TqenMm0q9exvjiTT296B29LK1itJoi60szQa9A9+utfnP3bS1l4Z+txPrZ5P9mzZ5CZlczHbp7Fp//3Bn6wYz6vPvGCBFMxOWIjnkpAFSLqEhzMvO5a3J3dNJ04Pr5garVB7mwyZxZz360L2DA3k3967RSVWzaPr9c7Eu31Zm5Ua5pqz9AE/N9XnRStv5lvfGgZHd6N7HzyaTOf6vfGzjicmPJi5f80CahCRILJWjcuDiTOVFqbO2nfu2V8z5FRSO6yZVyxJJ9Hr5nOnup2HvzHF6DqxNBDu5Gc27z4On4v1W+8xOdqm/n6Z66hpmYDZVu3mPWsnc2DX8dmN0PfoWBk2iUua7Hy3U0CqhDjlTuLWauWYbVa8PsDNNe34Wltg9ZaCPhJm7+E9ub2sV9fKTLW3sQH71yMK8HCzbMy+fGuCl58/Dlwtwz/+JyZptfYWDa2OdUR8J3cwz/8IMDnPr6O/z6UZ57H22mGnQeQfuW1uJKd1OzZaZKbhJgCJKAKMV7eDs4fNxWJElNTCAVDzL1iLikpiwkEQni9AdrPlw5zkUG40tj48P38/Y1zqO/qZuv5Nv7yyQOce33zyOcrG8tMqcFo9ga1Jlh6kNcPTOfGe9ax5fHfDxpMAdpOHWPuPbdgX7+esjc2RX+4WkxpkpQkxFTR2WxuFivdOTMh2EO1zYrFYkFrTY+/BxKTzTrNUND8nJINXe2gQ9DtvjSgKAXJmdz3qfv4wrqZ/N0Lx9m38wyh0gOjH9/SGvQEDK2Ggpzcspub/vp+Qp/4ANv+uHnw+WJ3C83NXRQWplKWUQgN56PfPjF1xUY8lYAqRMRYbeDzkLV4KV63F8/J/aaXphQ4kqBoIdSexTZtHtkF2VitFv7xwaUkJ9h45D+24jm0o+9aWcXceP8N/M362Xz6V/s59eKLQ/b4YkZrDb/58zG+/vBqEmy388ZzO/vWqfYXCnL+RBnFRVew+LorObFdQ1N57EyGibgSK//XSEAVIlJ6fOBpRWtNbmEmn374s5xq6CLFYaO9O4DDZuHwuYXsf3k7PVlp/NUHr+CXOyvJTk3k8S/dyO8PL+TZZ/dDUwUrb7uab29cxDMn6zj16qb4CKZh7U2t/M333+LOO5fzsYdv5ckfe6Hlkv0soKuDoNYsnZPF9IIbOXKygZrdb5ueuxAjpImd72ESUIWIpOwZ/Ozz1wLw092V1DV7sNut5KY7Kc508Z/vX075bXNp9vr59a4q3nnhLehq5+0rr+FfP7WGL//7A7R395Bst/GdLSX8+YmXzZBwvLDacKYk4T22ixd/fIK7PvMgqzZex/7f/P7C82x2bn3/DWzfepoTKS4+/p4l3DA/kz3LC3juqe19dYidKeDrGnz+V6ozCWQOVYgpqWDhXFw2Gx/8tzdpP7jDBAJrgulhOlP4zdJVdLR0oCtPXpBU1HLqOH/1Qz/Lr5xOXrqTY2ebOfPaayaYxBOtSXQl4s0qhqYKXnnlCB/5wBr22+x9vWylKLrhNrz+AL4zB/CFgvy4s4s7Nizh6llpzPnLO/jd5rnUbHsdS9F8QqEQlB8dOEO5fzC1O02VJwmwYpJIQBUigj546xy2VbTQfuZ4XwDoDSTeIdaittXhPtzI24d3xvfazFCQ9uZ2shcsoKm1htC5Qxw9N5vsq26g6cBu815kz+D9N83mf3695933xndiN89XlbJ39Ro+smEun71nIT9Wiuqtm00ylyMJfB7znioFCYkXZjlbrCTMXkZPzXnZ4eZyE6Wt2PpTSlmBfUC11vruwc6TgCpEpChFUZodf0CP7UM9HgNp/55nWKhkP00tRZA7G+rOcmjTTj75uXspuHcJbd4AuSkJ/MdTh+g5s//Ca3U0UvPmy/xryVIefuh6PnP3fH5mVVS88arpfTpTzbrb5EwTYJsq+h5rtaGUMkFXXHYmYEzir4CTQOpQJ0lAFSJStOZApZt7F+eYXtVUX1tps5ut25orL72voxHSrWaz8boSfvHdn0PhAlPEocc3dPnFimM8/piXBx+6hYfvmMeP/UHqtr7aN5QbCpp1tf3nTxMS8Xe642+IXERENEf5lVLTgLuAbwNfGupcS/SaIcRlJKMQMotwd/ewujjDVCea6rKnmy8OAwn4TaD1tJrA2+OD8iOmVzmSWsb1pTz10xcpbfLy6D0LmbHh7r7nCvgh2GOCqsVqbnYnNpfrwhKQ4rKhx/jfCP0H8HdAaLgTpYcqRARYswu5954VfOm62RypaTf7hEZY2pobaS8vj50iCN1uHIUz8dWXDtxF0NrMc/bufDPa5TCtNTz5m5187uHrefTOebw2I50tm45AXal57oVruPKq2aS67GQmO9Ba89Tj9RcOBYvLw9h7qNlKqX39fn9Ma/1Y7y9KqbuBBq31fqXUjcNdTAKqEOOlFGuunc8/bVjAv2wt5cmfb47KWkqllNnjNFZ0NuFrTDQ9xcEqMfm6TFBNLzB/jnY9bfkRfvJrO++5ezl3LM3lhgU30+hej9Wi8PiDbN5dQW15PQG/H0eSk1sfuIkDB6toObhLhn/FSDRprVcPcf+1wD1KqTuBRCBVKfVrrfVHBzpZAqoQo6GUCSA2u0mUsSZAZiHf3LiI10vrefK//zj0Livj0PbOW1G57pgFA6ZO8HC0Bnfz4MOxwyx3CZ7dx3M/Ocv2FWtwJNqxWi04HDbKzlQTOLP33fN8ickcyUxjydICsq/5MMfONnF2z2HTRllKM2VFM8lXa/1V4KsA4R7qlwcLpiABVYjRSc0183cWKziSsGXm8NAH15KdYuefvrf/8tw5xZVmMnAHSk7q1ZugpZTprbqbzTG7E/ucK/A31Q8dnLvaad75OjhcZr7a7zXzs/11u2l4+3UadoRg5gpWXTOf+z5yM3WtXqqrO6jcvStqX3bE5IqV70sSUIUYjQsSahoJhIJ87IoCfrC7grZDu+Nz6ct49XSbIv9DSc40vXl/l0kumrYYKo9DKIiyqMGTmy7m6wKtyb/iCtztbtwHt194f+/731zF/rdD7Pd3M/OK+eTmJpF443rO7j9pqjDFyiewiIBRJRiN/Vm03gJsGeocCahCjLN8XUjDtoPVJrBcjnp8wy8R6r/LjreDxKREulNzsGTmk+hMJBRMo8eTN7IM4IZz1HlasRXOgaT0gUcFPK0myLtbKKs5TZndCen5YHeRsGANPafeGdtrFbEpRr4fybIZcXmz2lCzV/b1kJQyH/wjXX7R0cDfv3CML945zzwu3kVr2Um32wzTBntwFhZjtVpJmj6LxStmkZ6VQk+3zwwd929H3hyzdvViWoO7hUDtObCF/96sNrMmNind/B4MmOHd3i9KAb9J6Ko9jd1hN3PgYsrQY7xFmvRQxWVPd7lND8qZYj6Uu9pNosxIqh31+Dj4TinfunuxGdKMd72VjyI9JBoKmipH7ha8zZUmACYkcqziPI78Ymir7+vhW23gSmPGisV0dnTR4vdCa82ltXx750PtThOMExIhs9D0WL0dF+4zGwpCZhGEgnjqI7+kSUwiHTsj+BJQxeUtGOjb2eTdXlTAJB0NJ6MQOhoI9ARIsFhG9piJlpKFrXA2gbMHRja/G83qTqEgpGT1BUJvJwC+UBC0xj59Af6aMkjPI29mIXWVjfjqqyEUxDZnJYGz+wb95LTmTENZFK5kF8G0JLo9WQRLLipt2FQBiSnm7zfBEVdb4on4IAFViF5a9/WChgs+eXP44hfu5GhlG7X1bnzBCCUj2exm7m+gnVXGwu8l4A/ETrJUOIheIFwEw+/vMl9S3K3U7zjxbiF8NXsl85YUc7L8+MBrSxMcJLoSCYVCdFRWmC8F7uZLg6+nzfRiu9rMLjbnDsXO+yLGzAzfxkYXVQKqELbwnNoo9h0tWrqQh1YW8T23n1sW5/C/njocmcLske41+brg/MHIXnMsehO/hnp93k7wnr7wWEo2sxcW090dGLxQg6cNT0ODuXZT+dDjf+4WyJ5OqKkapi2CimOjfy0i9sRGPJWAKi4zDpfpwehQ3wdv8RLo9vQN/Y5A9bFTHK5bxu/++ynTy/F747K3k7BwLT3lJwfuOUaKUia7d7Q78FhtJBTOorq8ge7q8qHPbamC/LlDV20CE3TdLeEvUB7InRU7pRzFmMVIPJUsX3GZCQZMMM0o7Mvs9XVhSUoZ3XW62sh2OiA1x/Rs4y2YOlxgsdJTWx79En16jNvZAT31VXQf322Skobi6wJ/t+l1Dpdt3VsQwt0SlZrLYuJpPbZbpElAFZeX3gzWluq+BJzeod6RFhcAMpesINWeEBubWbvSRr8MJCHR9Bzb62P3y0AwYALpSD/5GsuwJyUNH1C1NoE0JdtUvRJxL8q7zYyYDPkK0VJNaDQf3MD0mZk47dbhKwRNBJ9n9EHR3RKdtkymgB9/wzA92f46GkwpyeF6v0KMkPRQhYBRj/+0t/vo7gnGxgK4YCA22jEREhxD90Bba8w86kh4O03v3uGKTNvE5ImRyg4SUIUYg7qqJryBKA6VpuaYggWXE5vdVDoabPjalQbFi4e+RigIbbUXHhsqCLfXT42CHJe5GImnMuQrxFikZKTgD4bMLiujWG4zIkqRNn8J7S3tcG6US16sNjM3GE/JNjY7mWvWs/HGeSzMc/Hy4Tp2/e6FS99Xnwdqzg7/fmtt5oddaWbtaeY0rKkZBM8fuXTZzlQc+r7MaCAUIyM0ElCFGK2sYn7/xfV8fdNpCEShspDWtO/fOvKhy/6Cgahsbh5N9nkr+Z9H11HZ6WVNfgb3Ly7g2uO1l+7/GgxAcIRfXqwJkJRhAmrDOYKNlthNvhLjFhvhVAKqEKNjs7P0pqvo8Pew9ck/Rb532kvroddTDsXvjWxbosx/Zj8f/maAkN/HDbev5LcPreJ9dyzm53WNJjDWlYSDac/IK0gF/H3LY8bzXorYJ7V8hYhf7e3dnGjqvLx2LFHK9Jij0csLBgiF6+5u/VMHL64t5Ks3zeEjywoIac2Ztk62lbTx8hun6Ti4Y+RtuBw3exeTSgKqEAPp3T6sufLCgvEBP5Ul1eQmLTVrOS8TavZKsguyaXx7c3S7Ay3VfOZbr5A9LQ+L1YLDkcAVC3O4b3keH115E5/7RQoVm16M3vOLODQxG4yPhARUIQZis0NiEiRnmXlSu8sEV4COJnpCIa6+90ZKS5poPHmy774pSrc10hit8jL9hYJQfoSmcszuPVYblW9pXrI7uf7D7+GhDXP55o60uJsnFtFjkpImuxWGLJsRYiA9Pig7bNY1drWb5RW92uv51+dPc+WsTJ76y/VMX7V88to5UZorR59xPF6hoPl7CPgh4Gf/vjJe2FuNY9aiiW2HiHmxUilJAqoQwwkGLlxuEeyhdPcBfvq7d/CHQvj9kvASdc5UEuwJlJ6pY/mqWcOXFxwtm90M8V9O8+JTiNTyFSJeaQ3NlbiSXWwpa6Fu55bhH6NU7GxA7kwxu6zEk44G2k4dw332OE6HjZzV10T2+gE/uFLNUhsRd2KlhypzqEKMkdVmZV9Z68iWqdjsJC64ku6TeyO3efhYKcuFiVbxoN+ONVufehkcSZG9vs2O1WYl2B3FbezElCcBVYhxcDlG+E+ox0f3mUNjD6b5c02SVF3p+Ne+drVPXFJP78bikdTZbG6RZLESDAYhLU/2R40zkpQkxBSQnZvK+rmjGCIcT8EFZwrKajXl9GJN3hxIyRr4vqzi2KpJnJxpyjNezO81SVct1RPfJjFuMuQrRBxLWn4dn9s4h+89e2pinrDskPnnHyslYfrr6R78PnfL0HuO2uxmIwCfZ2IKMVis2OasJFB5yuw2c7GLa/2KuBAr/yykhyrEaFmsPHT/Cv7xZ3up3DxBRQYmYg3oWDhcLFi/lpR5Swa+P6t46Ixcm52k6bMgoyg67QOTDKZUv5+5rIpyTHljzPCNxj8n6aEKMRilzM4t4WSYd1ltJFgVvm6f+YCOt6LrDhf4usZ3jYxCbDlFBEoPUna2Fr9v4J5dan4uHT3dg/c+fR481VXQ2TS+9gzF4TKfnt1uaKsj0F4/to0HREzSQEgqJQkR47Qe+IM+OYskuxVfXXX8BdOMQhyF0051yp8AACAASURBVPEd333hcaXMkpGu9pG9Jk8rAaUgGMB3YndfD/AiHUffGXpIWGtoLIvu+s+Lh3alWP6UEyuDN1H7mqaUelwp1aCUOtbvWKZSarNS6mz4z4zwcaWU+i+lVIlS6ohS6sp+j3kofP5ZpdRD0WqvEAMa4F9q8sw5dAdCEC9LLKw2yCyCBAdZixbja2oc4JwEbIVzzHzmSPi9FybwDPaJNm0hZBQOfp2kdApuupO1H36fSRiKJlfawAlJQkRINMc9fgnccdGxrwBvaK3nAW+EfwfYCMwL3x4BfgQmAANfA64C1gJf6w3CQlxgooomKMU162az+WBN9LZuiySLlem33MH/+9oHSVu+jmAg2JcI1MuZAq40AjWlke8putuGXDOaOGsRRUWpvLPtRPSX8nS7J38NsIiCseb4xlHpQa31NqDlosP3Ak+Ef34CuK/f8V9pYzeQrpQqAG4HNmutW7TWrcBmLg3S4nLncA3dCxqPBMclh5ravCybkxX783AOF7nX3UpaWiJf/d4m3O1url4zg889cgtr77u177XZXaBDJiO3qWL0zzNUFahwHd7BdB/fzb7f/tHUTY728Hm8Dc+LEYuVpKSJ/kTI01rXhn+uA/LCPxcB/bfrqAofG+y4EH38XlPEPtKUujRD1e7kwGu7WFaQFDulBC+mFExfyrT1N9Ht6eboy69DxTGCJft59YkXeH7beZbPyuQjX/4kas6VZo5zJIUSlBp4rjRzGgnzV0FS+qX3pWYPP4cqgU6MQ29S0lhukTZpEwpaa62UitgrUko9ghkuhoQIF84WsS1aSSZaXxJobn34Ae5fmcfpBm9szsdZrKStWs/tN87nD3/cAxXH+gKW1uBuofrNl3my/mruvH0JusdP7qqraTh50uyoEwwMHOCcKaQtXYUr2UVXZxftx/b1DXm319OTnHFpjz0pHWW1oi/OkhYiwmIlKWmiPxHqlVIFWuva8JBuQ/h4NVDc77xp4WPVwI0XHd8y0IW11o8BjwFYXLkx8vaKqWb/gUo+va6Yx36/I3bmUC1Wlt7/XqoqmkmwJzCtOJ0/PL33wmDan9Z0H9/Ns+G6wg2dLRAKUnDdLXg9XtpOn+jbrk4pSEgkbekq/N1+2kvPgiMJx5yl+MpOmeUwAT/0+EiauxjPoR19z+NMxZXiwpM3BxrOmbbEyiefmFJi5f+qiQ6oLwAPAd8J//l8v+NfUEo9hUlAag8H3deAf+6XiLQB+OoEt1mId7XWNnCyyYO35OjIg4PVFr1kGLuTObdt4NqleRxLc/D289tpfLty+LZp3dem1hqwO6k9XUpGcRFFq1fTVNdMwB/AZreRlJJES30LnO/bD9XXaL+wAlJHI57+a1tdaZCYQigYwpGZic/dbJKTWmskMUhMWVELqEqp32F6l9lKqSpMtu53gD8opf4CKAc+ED79ZeBOoAToAj4JoLVuUUp9E9gbPu8bWuuLE52EiA6lYMYV5Bbno7Vm+bICCjNc7CxpYfZNN3HuzTeH76UqBTkzoa4kKk1csPEO8nKS+Ml/PWeGXJPSoMVy4RC4zT50ST2L1cxD15ymteY0rXYnlumLSctKIxgI0nLm1CVz1JaZywjVnTeJTGD+7P05JcusaW2pwmtR0Fxl1oIGA5A5DZrKpacqIkrHyP9PUQuoWusPDXLXLQOcq4HPD3Kdx4HHI9g0IQZms5O28hryCjJoqGsjNT2JT985n+1nmqmud3P0RD2vV9RDUwW3fOQuVn3hQ/zxF69Cc+Xg19Q6eruXKMWSuVk8++x+sybUlYYjvwjfxdWbdGh01/V7CZ07ROt5M9xLT3dfAEzNgZ5uQqUHBn98V7uZe+7N/u12m/ZM5C434rJhkpJiQwxmVQgxgZSC3NksumY5axblkZlko7q1G9eyfFZOS+LbvzlC445NlzzsjSdfYPndt/HXX7qPx35/gK4jbw/+HNHIYlWKxCVXs6wwhZdcDnwA3g58Hu+lAXS4IdaL25eUjmPmQnylx7AWLyRYfbYvEPZWUhqqR9D7fAmJWBMSCMoQr4gmfRn0UIWIC640Pv/5jVw3I42v/PYI5VveeLfO7ROJyYMP6XraOPz0n/B47uZ9dy/jaa3xHt8TneCZUQhJGVidToI9PTiSnMxZUMQtK83a2+z8LKpPhedpa8+O//mCAVNW0echWHbUbEZusZpqS6NZp+r3Eiw/Pvx5FquZZw74ZShYjEms/F8jAVXEh0gUdO9lscK0Ray4ZhFXL8rjdG0nP/yX34Kn9cIP9OHmR4MBSl56npLjK9h41woal09n33Ob+uYSxyItD2x2EnIK6Gltxp6ZzewFRTidNopykmlx+9i/u5S6mjZ+dLyClPQUvvShFfy7p5uWAzsj8x6Fgtiy8gk0V5pg6nBBYgokZ/QF1ASHuS8la/A1rClZZm7WajPveY/v0nNSc5h+9VUsW5DDSy/sh/Ij42+/uOxID1WI0cieYYqoj2eTboeL7NXXMm9+DvevKmDP+Q7e2FfF2ZdfHHvmqdZw/iCv/LqWzIVLeP9n3seBk/Wce+XPo+9t2Z3MvvYq3G4TeDyJdjx1dZzaVgV+LweDPaadfi+9oakV+FZnF5/4wCp2z8jk8POvjH85j99LoOxY3+8WqynaUH/O/J7ggGmLTdbvUF8eMov65o8HeS8Sp89lyfwcvP4gBIbYN1WIOKBiJbJHksWVqx0LPjD8iSJ+KDX24cCkdG7+6D0sK04jpOHZLeeofme3CTwD9ZrGY9ZK7rp7Ba0ePzvfOnHBUpNBOVMgdxazFs8kFNKUb9tiepqjGT7OncWG91+PPxBiyy/+EPnX1Z9SYHcO3xvuDbyh4OA9T5vdBOwEx8RsMC4mTfehH+7XWq+O9HVnL16uv/2bl8f02A9fOS2ibYrxYqRChPUGU7tzwPq6l0hKp/i2u8m65la+8a1Pkp2ayM+ePsh/f+Mxqt94yQxTRiPonD/IS796Ga01D7xvNRTMH/p8m511D76H+9+/lp6eIOU7d5olJqOdi20sY9Mvn8PrC7Dy/fdGt4qT1iaYWm0kr1wP+XMHPi8YIDk9maSM1MGL7gf8ZtTB741ePWYxxWm0Htst0mTIV8QXh8sEm6GCYWYR//srDzAzw4FVKTw9QZ59/iCB1sbo9tx6tdeza/Nhbr93DZ97dAP/889Vfb05pSB/HikF+VhtVgqKMshLd/Ls73dB1QgSeAajNXS1s+eZ19j4iXtJWX4NnQe2Reb1DMZmx5nkxD3YB1MoiPvQDrMtW3CY4Vyb3fTU28YxEhEJvfWIA/7IzdmLqIuVgVYJqCK+DFfE3e7kS19+H3/YVkbJsfNkFmRjS7ChLCo6BfQHE+xh565zFN21lOk33EzF5pcgZybTrljMJ++YR1dPiBZPgMPnmvjzi4eh5lRkntfdQktnN8tWFLNziKWiEeHrGnBJ0QUGqIc8IG8n+EomP5gmpoC7eXxz9WLCxUg8lYAqYkDvDiaR+DC1WGly91Cyaz/OwmL83X5a9myZ2B1NcmaCNYHOA9v4XU+A9961jJqZHyU5MYGOLj/f/P4rpj3eTvPhHcl1mvlzefCqIv7tmZN9mbixwGI1vdChAlW0/44sVkjPN8UqWmv62mKxQkq2+buQedy4o4FQjHRRJaCKyZdVjErLRpcfG7pE3kik51PR5AZnKt5OD3jdE789mM0OtacB8B7dyW+bmtj4vms5V9XOiRf+PPBrTHCYWrejWXKTlkfWkmU0n68w5fyCAQoXz2dtQSZfeWAxtbfP459/+CZUnZj8+rmJyZBdbAr2TwZnCusefA8rZ2ZyrsFNm8fPmdP1tOzdbqo/dTRO/nskxiw2wqkEVBELggG01zP6Enm9lHo3W/S9H76RaRlOttjtEAqZtaUTrb70wiBee4ZXfnBm6EzlkX6Yu9IgZwZ33bOKG+ZnEgxpMl2reO1UM+0ePzmpiXxt02lef/LPTL/mGr79v27lTOO1dHb3jH+edjy8HVAZxedOz4e2uoHvS0zmlofeSzCk+Z/vPGHea4erb7s5v3f8X+SEQAKqiAWtNWZB5VjY7DBtEZn52SSlOLlnaQ4/2VYOjRVmuHMytlgbrEc81LBUKDh079SVhrVoHg++fzVNHd28c7Cal373hvnC0Fv8PncWVrudYMUJ8HVRselF/s9WFyRncfdHbuO//+FOvvW7Yur37ACfZ2LnK+1OsCZEpZZv8sr1fPljq/nWT3YQOL237w6lYNYKvvb5G9l6uoktj/++L3D2H3oeTyEOMfmk9KAQEVI4H7TG0+nhgTsW8tyRBnZvO2nmwiZ6qDdaEpO5/sPv4S9vmMU/PHOcUy+/MvBcZPkRLnnFvi7wdfHif/+CV+au4tGPXU33DXOobvbw+uvHCZQcnJihzkhnzPauhU3N5aP3LedwVScLl83gWG9AnXEFjz58PTfPyuAbL5zk2LN/mjr/P4gLSHF8ISIhOZP0/Bw6WzsBKGv08NoTL0ydHofFChmFfOhTd7JuZioPfutVAmf2ja1nqTXBs/v44beOQ3o+C9at4OMfWEMwtJonfrbJDFNHg1KQXgDdneZLwHCF9TOLTK97oACsFCRncv0HN6KUIsFmoam1i9f3VpGW5uC9V08jK+NjpLoSyEx28MrOcn78H89ObHa3mBQx0kGVgCrijNUG+fOwp6by2Y+sZWamg4NVHgKhELsO15p9OHt3RIlnSpGx5gY+/p4lbD1Wx+/+48nI9PL8Xmg4z+nnz3P6BUXhTRv5xlfu4Zs/eZueswcG7q0OlIUdzoy15hYTrCkxbVPq0gL3BfOZtXwB54+WsnD1Aq5cmMO+Ew2c2fy6yXLuLzGZ1Flz6ersItDeQnJhIfZEO05XIikpDqYXpuDxBThb2kLt+WqT4NT795zg4PSJteTkZ9Dc0E7H6aNT54uVGJZk+QoxGkpB3hxW3bKKaxbmcr7Rw+GKNn72x/N4Tu4HVxqzr1nL/KuuICFhBT09Ic7sPgSJSeBuNQkrcRRkc669je984kr+6Q/HqNi+NTpFBrSm5s2X+ceyWj7/iWtx+1bxxG93Q9khExQTHCZwWqwmA7k36cfhoui6m/jKA4u5bW4eh6rbeLuinRPV7bR2dNPZ6SchwUpnp4+UFAcnNm+BzmZO2e3MmpZGe7uXWTfeRENtC56mVkiwUzirAKUUjTVN2BPtpGTMICnFSTComT0rgwP7yji167DJZh6ol9vjw31wO5MwYy4mmSZ6PVSlVCKwDXBg4uXTWuuvDXa+BFQRHwoX8ugjtzAny8F/PneKuooGAn6/qRGrNfi9nNv8GiRlkDRjLouWTeOBj9/KtXPSONvYTVmjh/1H62ipN72WQE+AlIwU3G1uguePRDfLMzHZZOe2VI/sfIeLjdfP4cs/fYfWd96K/njWuYP88LvnueHD7+GN797Px39SaMozKosJpJ1NEAyQuHQdd9++lNsWZpLhMKUEf7q3greO1lFZ0UZ2TjJL52QxO9vFzpJm3G4fJ7bs7ivs0FJLVf0c6ne+hWXGUgpn5pOWmUpjTRM1Bw9BWy1ojR/ooi9Prfat6L58IYbgA27WWruVUgnADqXUK1rr3QOdLAFVxDaLFVJzuOU9ayhv9PD0ptM0nTpleqxd7RcGmx4ftNXhaa9n31EL+2x2/piejy0jmxlzC1m5NJ8Ft82lJ6jZf64FR4KVadlJbNqeSePON6LXgx2u7N5FHHOuYNfRWlr37xhdMM0qDmcLj6FOsaeNrY//jkdq7uI7H13Op2qbsSXYmLuwgFuWF5DlsnFNUSbHmtv50RvnOLr/vBnq7Vd3uBE4OdRzuNJYOCODo1oTKtlPVcnomijEYKI15KtN+nDvwEdC+Dbok0lAFbFp2hKyphdwzZrppLrs7Dpcw7mt20ZWyUZr0MF35wsDDecpPQ2lSvGK3Wl6Xn6vGdLMnsHff3ED/1LbhK44boLDYHt3jlaCw8xJ9vhG3ju1O8kuyOLs9t2j6zVnT2fx9WtIS3NQVtZG7eHDZu2n1iNfOhQMUPrGG/wg28X3/2o9ea5EdlW18/M/HaWt5AxFK1fQXN9C97FdI29Xf51N5vtB75ysEJGgIRTFQRyllBXYD8wFfqi13jPYuRJQxcTqXTM5ELvTJBX5u/jG32wgy2XjyV3V7DlYTcOpM+NfU9q7S0ovXxdUn+SnzxfysQ9ehde/isaObnLTnLz21hk8HR4C1Wf7dkQZ7bfg/kHZajObhw8VWJUi96rrcToTTOWekbJYyV28iGAwZILpoUPgcJE6ay4+rw/fiQFHpwbW7WbPr//InmdS+0YBwolK1edzSc9Op3vkV7tEWX2nFFEQETXOOdRspdS+fr8/prV+7ILrax0EViil0oE/KaWWaq0HLPklAVXEBquNmTffykdumYPWcKTKzdO/3Qpd7eSvvJL8xQtpby7Ae/pgOJM0FLG5xZbdb/Krw/tMCboeHwVXLOP+OxZz1fRk2rqv51BVJzsP1dDW3Dn27FGtURl56I7GQUsPLr33PbjdPko2vTa64eeUbAL+AKdf23xBwYaOpopRDzcDJoAO9BrdbThnFDDWare2gtmcPV0XO2scxBShCY29+GDTSPdD1Vq3KaXeAu4AJKCKGDBQMLFYKbrxdlYvyef7v3qHlPQUGo4fe7c3V7dzC+TMJDk3m7x1N1BcnI7NZqGpqYvGujZQ4O/24y0/a9Yw9gajkX5wh7c+663iU/tWJb/YlcwvCucxY940gsEQDoeNxcuKSLt6FnsPVNJSeg4azo/qpVttVgIXD3dmFZM53+yZWnKqmu6Te0deaMFmJ2vt9VgsFhoPvnNp9aNI9wTb66k9MfYtlF0pLgI9Ui9XRF4Us3xzgJ5wMHUCtwHfHex8CahicjlTWHrX7eTnJPH0T5+Djka8GYXkLFlKY4nTFBwI9kBzJW6fh6JrV1FT08G0aWlkZ7vw+QI4nQlkZ7tomZPP3OnpdHT5CYQ0ra1eykrq8FaeM8PFPd3mX55SwwetbjecO0j5uYPmd5ud0uIl0NmMLSuPpMIiPKMJqKEggaqzFw4DFy3ii4/eyp6zTbzz/Jsmm3akPdPCBVy9YRVpSXZe+9WfJ2aXFF8XNJwb88M9HR6cyc4INkiIqCsAngjPo1qAP2itXxzsZAmoYvJkFfPAJ++gssnD6796vi8otNbQuLsJ8ueRuORqlFJYbVbSs1I5vW0PtNZwce2bMwB2J6eyp2NLSSXRlUhOfgZfffgqEm3rON3YzdtHainKT6aiuoPK8/XkFWWTn59MMKjxens4fbSCwLnDA/fsbHYTUJoqCHg7CIwlaan39aXnc/37b2VWXgo7TtRzcNMuaK8f3bUay9j9SsC0ayT7jUbKOLoCVpsVV4orsmtFpy81GSnVJ2Qo+TJltm+L0rW1PgKsHOn5ElDF5MidxQMfv5Utu8tpfHvzpR+GAT9UHae71vZub9I91G4tYBKHak4TwOS5u4F/3F0I6XnkFOUQDAaZXpjCmiV55OQk0dMTpKPDR3aWi8xUB3euWUd580qOn2vm+L6z0N7QV8zAasORnoHPk2fWlNadHf1rttpwLFzNJz+wGodV8cNf7SZwdoylBHt80SsXGCXJ6ck0VDVE9ppZGQQDQbw1FpPZLS5LUilJXNZUSgZFGYk07t4ydEDpPzQ7ln80rTWmxxsend20C5Nxq3Xf8KrFCnYnm9LzyVswl8LCVK6/7QoWFKVxpKyF6TnJNHZ04/UFcC8qoqnJgzs3l2AwiNYaX6cbR3Iyvi4vdHWYZTkWq5mTTUonrbCA1Ixk7rhmBscqWvnZb/YQaKg0mbwT9UEw2ZuNJzjIzEqh5URk90N1HxrlWl0x9ejY+V9AAqqYFLq9iU5faHLWJF48fxoKmjnTuhLq60upD5fd2+ZKh+RM9iclk5yejNfjJbcwm7lzs8hJLSIrJZGOLj9aw9xcF63eAMcqzbDulTMzKWvyMDsniauLUznW4OH1Y/Xsef6twfftjBalIHc21J6ZvPKLaXlkZjohEOGgHiufpGLSRHPId7QkoIrJ0eOjON1uyvJNZs/pYr0f0D0+M6/ZXk8AzFIRpag8rqns3dDckQQ+D7bZV+BwOgCwWCyEQiGOHqoyRd7bmvh50TQ8HR50+bHJWYOpNdScmtzg40iivt4dnZrE4rIn+6GKy16zJzAxe3FGSu8/Wq1NwA1/EQicPUDAajO9P5vd3G+xvluIwj3K5TVDGm4eeTCT/YET7KGjzSNFHcSUJgFVTCy70yQP+bs4XdsBmdPeXf8Zt0LBvqHUaH5BmL7UBOqyw9F7joFYrOMfKrZYaS0rj0x7hLiIDPmKy1NvD6XbTUNLF460FHxVk9ukuFFzZmyVj8YqNcdUTMqeDo1lY+/lOlNMb32k9YyFGIVYmkMde9kTIcaiX0+npyfE8lWzTK812pSCpHST4RuvLt68O5osVjOSEAqOL5gmJoMr3WRbCxEleoz/RZoEVDE5tObsiSpsNsvEZPrmzWHJHbfgXLwWUrKi/3zRZrFG933TIVPKEMYXxJMyTDnIyZ7DFVNXeLeZsdwiLY6/rot4l+hKJMEaXrMZbc2VHH8b08uLp0SowSSlmx5kb9as3WleW6SWxUQqAHY2SSKSiLpY+b4mAVVMHKXAmmA+YJUivygTX09wYgJcjw/qptCO1heXG0xwTOz86kgpJfufiqgyc6ixEVFlyFdMHIsVkjMv+IBtaHBPXrGBqcTTdukXk/T8yQ9m1oSJnfdNcJje+mS/bnFZkh6qmDjBwAVVgs4cOcftd63kXEahSXwRkdVeH/lgphQkJJrh5pHoajePScmKbhF/pVh6/3spzE2io6uH3W8dg/Ij0Xs+EVNiJctXAqqYPOVH2LEjhbzFC6k/YzcBIN7XpMaS4YLpWIpEaD3yYNr/MROwvZzPF+Ds+Vaa6tugWdZiXU5iZMRXAqqYRFrjcDq4clkBB4D6g62T3aLISM/HXjADf+mR0QefiZQ3x4wYdEd0Q7WBRXtYX2vOvv6GSdKSKYTLitZa5lCFAGgpr6Slsxu73WY2AJ8KfB6KZ+ez6K6NE7PGdqyaKvqWxkTKZM5dejslmF6mtB7bLdJGFFCVUqVKqd8opR5VSi2JfDPEZav2DKUlTaa49UQsnxktu9MkUo2Gt5PSV1/i5KHzrH3wPsifG522jdd4CkUM9nflTDWJQUOx2mD2SrOv7ERRyrQ5qxgyCi+8z2aXJKY4FxrjLdJGOuS7GLgKWA/8q1JqAXBEa/3eKLRJXIayslxUxVpAzZ2FPSsXf9mJ0T82GICyQxxPS+LeB29gy9vTaN+3NXYme8Zr+lKoK710uHiEc+B5xXnUt07QnHlGIbd88Fby0p1Mz0ikOxDi2a3nqdh7ABwubGkZBOorJ35bPTHljDSgBoGe8J8hoCF8E2Lcenw9OJ0Jk90MY9oSsCioOAbt9fhba8a+vZzWeA7t4M9uLwuvmElq+l1UbtkcW9vVjVXFsbH36oIB6rdvnrjhWW8HXn+QexblsCAnhe6eEKuLUvlMfRuehgYCZcdlW7k4FkvrUEcaUDuAo8D3gZ9qraOY/y4uKxYrvm4fzc1d0S3woNS7+5cO1UtMyspAa01XBRH7kA2VHuBESx0z117J7Ntup66qGb/PT6Dq7PiyX5Mzzbyzr8sMTdvs5vVNRKGM8QbDiZzr7Olm5zOvc+zwQj71vuVkuGzMTnfiae80m66LuBcj8XTEAfVDwHXA54BPKaV2Atu01m9ErWXi8hAK0t3SQoPTYerHRoNSpK2+gfs2LGL/6UaOvfDSoNm3niNvR34uV2toqaZsUx2kZOMomsW0Wflcfe9KtrxTSd2psyZBaLgSfUnp5pzeQN/t7ntMuPrUZZmUY7UN/SXClQZZhXQc3MH3j++HzCLSigqg9uzEtVFEVVytQ9VaPw88r5RaCGwE/hr4OyCGUxhFXLDawOGivbzcZGmCSRrJyINzByP2NHPm5XLTnHRunZvBQ4fOEyrZP/CJWkevhxcKQns9PpudsoozlB3NY+m6JWz81AaOV7TyzjOv9C37GCgwXpzF2j8Ah4KX57ClUqiZV6DPHRy8m+Ju6atz3O2GmtO015ye0GaK6DEZu7ERUUcUUJVSzwDLgVJgO/BxYE8U2yUuF6EgdHvMh15Y8sw5WCwWOkZ6jeRMk9wyWO9Maw68sZ9/D2qOHThPqPLkuJs9Lu5mMzxbfZJjz57h7KI1fPkTV5OQcBdJDhtlVR2c2RX+MtG/glS89D6H6zFGWGZeJs1l1oGfMzHZDIv3fllTClJzTRERMWXEVQ8V+H/AQa11nPyLFnFD6wvnsSxWktOSqTt6bGSPt1hReTPRdeeGno+sL+XwH87FxmSLr6uvNxkK4ju+m2//wxFIL8CVl4/VFq5JG/Cb3npX26W9z7FUOZooNvvEBVStad67feDnU4qUxVfSWX4emitNuwJ+E2QloIooGGlAPQx8Xil1ffj3rcCPtdYxuL2FiGsF82ltaDUfgIPJmWnuDwWxL1iNPdGOu2YE3/ViNQCBCZj1pXTVl14YLB0uSC8Ab4epoWuzg6cVx/T5+FrCW6N1tY9t+UlSupkv9nWZ60JklrGMZOjZYoXChVB7evzBd7Csaa3pbAiPBtidqOJF6PJj0CELFKaaWOmhjrRS0o+AVcD/hG9Xho8NSin1uFKqQSl1rN+xryulqpVSh8K3O/vd91WlVIlS6rRS6vZ+x+8IHytRSn1lNC9OxB9XVia+8lNDn+T3vpvA5K8tx11VOTWWovTqH/jDgZbOJhNU2+txzlrIvMXTKFqygG/8n/ex4RP3jT6RKsEBrvS+xKa0POwzFkb2dQxGKVP2MNgT/Z5s1XHz/vm96POHzWvtHf4VU4LGzKGO5RZpI+2hrtFaL+/3+5tKqcPDPOaXwA+AX110/N+11v/W/4BSajHwILAEPDSXBAAAIABJREFUKAReV0rND9/9Q+A2oArYq5R6QWs9hpX2YkIkOAYPbhar+TAd7EM0MZkER8LgtWWtNkjOunC4LtKL8V1pkVt6opRZzuL3jr93rPW7vUfvmUMcO38KlOKPmU4SEiyjWxNqtZlefsP5vsSmziaSZxbTYrFGd662N5haE6Dx/MjOV5bh25Q93YxaXPw+O1PM47vazc/d7qmxwby4QJTWB4zaiAs7KKXmaK1LAZRSszFFHgaltd6mlJo5wuvfCzyltfYB55VSJcDa8H0lWutz4ed9KnyuBNRYZLFCSja0VA98vyvNJBDVl176wacUSQtW0H6udPAPvPQCbFn5BMYz/5VZBB2NAy9RcaYw/br1NNe34Tm8Y+zP0ctiNXOgzRHuQfebgz36zJ9GnwSUmmt6vP3fA3cLLZV1LHjPPQSDmoQEC2UldXjPnwJPa2SGy5WCtDwTHGtOjSxwK4t5H4c7t7Np4Db6uvp678MtSxJxK1ayfEc65Ptl4C2l1Bal1FbgTeBvxvicX1BKHQkPCWeEjxUB/SfNqsLHBjsuYlEoOHgwBfB2YE1JN3Nn1ou+y6Xmmn8U7iFqhihFoPz4+NroTO2bL7yY3cXt62bwtw+tiUxR+2AA6kqiOxwdCg5/fZvdzMUmZ5pNx/1dA+9NWnWc0y+9RMnW7Zx8+yDOZCcf/sIDMHPFpdcbLZvdzAX3dJv/R0baCw4Fhw+ESg0+b5uc2fd4X5SLh4jJMcbC+NGIwcP2UJVSVsySmXnAgvDh0+He5Gj9CPgmZtj7m8D3gIfHcJ2B2vkI8AgACcmRuKSItGCAYOlB86FesOCChBRb/nS66mr6Phj7D/VlFJI8cw7u4/vGvx1a3dnBP1Tb63nymQOkZqSO7zkmgXXeaoLujksr/1htJqs12GPey672od/DgN8E285mWpoqOFSQzvvvX83Tv2o1S3h6i8yPRFoeM9etxeGwcfpQKbhbL8zEVsoUq2+rG3vv0WKFgvmmxztgD9XTN+wupiRNHG3fFl4q8yGttU9rfSR8G9NXbq11vdY6qLUOAT+lb1i3Gijud+q08LHBjg907ce01qu11quVTepNxASL9dLeTG8vtrvTZK1arJCcScDTdUFmb9IV15q5NgBnCu7K8oE/FJUyAaM3Y3U4w/RQuo/tomH7a3H3ARw8d+j/t3fn8XGd5aHHf88sGu2rZdmyZVvelyy2Y2cnBAgh5AIByhJSIG35FLgNt3A/9LZQblu4Lb3doLdwAzRcUgjNQoCQBAgkzuokJon3fZNl2ZK12tq3Gc3Me/94z0SLtcyMZjSLn68/8/HozJmZ8+qM5jnv9rxw/oz9QcTWzCL91f2ddiDOTMF0onCII088wesHWviLL96BZ81We5Ez02u4PXjXXs3ffu33+bs7r8DlEnvOJ47cdrkpWLwk9tV8xhKZvO80YsRvLyaUmgPR9qG+KiL/F/gJ8OYCisaYPbG8mYgsNMa0OD9+AIiMAH4SeEhEvokdlLQKeAMQYJWI1GID6Z3AXbG8p0qhqTL+gE21V1yJu2YtofZGm5A+0t9oDAMt56Cw1Nawpstq4/baL9UCp/cg1ty4Xh+UL7b9hGOSS2ScUHD8xUK0Nb7I0mVTNRuHgjRu+xXfbruBT991Da8eWcb+R382ZQDzbbiWG25czWduWMIb53q5d/9pjr68a/LfbSjIwL4JfdWxzq8NBSE0zQLpmZIMQ81KmlRQow6okU6Ur03Y/vapniAiDwM3A/NEpAn4G+BmEdmIbfJtAD4DYIw5LCKPYgcbBYF7IkkkRORzwNOAG7jfGDPLTjSVNvrOE4pM2zBhpLIGEwzYmuuFRlurmanPywkcK6++ktZzF+jf+3L07+9yQ816CIVGa3fZwJipR0qP5fbYJtf2+un3c7kZPH+e79z7a771tQ/x63l38/SPf20Hd0WIUHLVW/mXz16LMYaf7m+joa2P40eap69BjpWTZ7sD2qMY/avUGOkyKCnagPpu4PeAZWOeM20JjDEfm2TzD6bZ/+vA1yfZ/hTwVJTHqTKJMaPNh6ERTNMxmLeUa27ZyOsPPWZrjufPTB9URchZtYnCQi/9F7pifP8wdJy1TaGX4mAVT479HUz8Mpo4dSjSzznYw5/+79/wwNfex4c2fYqvPnyQcy8+DaEgntVb+N5/u4Hn67v5/rcfR8qqKKsso6/LadqPpqk5txDyM6//WqVemsTTqAPq40A3sAcYdralSRFUVjDGNjuGQ5w41gbF8/HNq8Tf2TR9sCueT1FZEYf2NkBTjDl6jZl8tOulwOWGBSuh9/zM+zpZqQA4e4hPfmWEJ/75Tn71xZt45J0reeFIO4srC+gfCfHILw+St2gp3hwvnQ1n7BSpaPV32lthuU1icSle5KiYGZN566EuNsbcltQjUZc2EwavD3dhMV319dB/Af+RppkvPXvauLCz2375an9Z9EwYAsN2Cs1EE9MPTvy9Nh/na78+ypKqIkJhw48+cRUALpfw6J+/g2OdfczLy+Fz39lBT2td9McUeZ/AUPTnMpKfV13S0iOcRh9Qd4jI5caYg0k9GnVpcrntLRQk1HTcfkmO+KNvx8mmtINzZeKiBDHa89Cj7HEWbX/+xY12wfSRIKa/B4b7WbRlC1u3LOXZXXEk8Tdh8OSQv24Lgwd3TP18EdsPrAFVpYloA+qNwB+IyGnAjx19a4wxVyTtyNSlI7I6SWjEJiDw5toBKr6CqVcFKV9kg/D5s3N7rHOlsHz8AuLpyBn8NFlWqXM7tlNe/q5JnhSFET+I4B+a4ULJmEtzDVh1kUwclKRUcojgW34Z/q5O8OaAuMhfvIzB+mn6RHs77JSXbBUOxZafN93kFVNZnh//aBFjCJ3cldhjUlnJkD6rzUQVUI0xWTSnQM2JaJOaA3hz8bc5E//ziiAcYrB5hpR6wcDFtbfCclur7TufcYkZLpKIZdRSqbeDQX/QJq3P1lYElTbSpYYabS5fpWbmcts+La/PrmYSbc7XwNBoFp2hPtuMF0+/aE4eFM/L7pprboak1QwG2L+7gevf95aL8yL78lNzTCprpUsuXw2oKjGqVpC7/mqYt9R+6Q/1Rl9LHBke/7PLHV9y+s5zcHpv7NmSMomvIPa1T1Nk6PDrNDR0seD6m8c/UFiRHmVwe2B+bWY3rSuIcy3UZNRqNaCqxOhpY8Q/YgcROcnVo1ZYPv4LtrgS36ors7umGa+etsyZHhQO0Vx3lms3VY/fPnZeayoVzWPJxg32IkWpBIh2UJJS0xvuJ1S3O752FLfXTpWI6O/E3+xKj1qMmp3+LoYCoejWNJ1r3a2cffklmxVKZayMG5SkVFRiCaYitlm3osY+r7t19LFgALqaE3982aJqBQULFmKMIRQKERgOYMIGOs6kXYL//MXLbNPa2AumWLg9TiqcJAXjoT77fzoGfBW1dBmUpAFVza3cQqhajifXh9vtJjc/1yYFKCllpOFIdEndL2UuNxUratl4RTV3XFnFi3VdvL6vhf7eAXo6Zn762NeZiwAy2HKOQHB59Pl8J/Ct2ULN8gXU/erx0Y1VK/CUlhM8vjMxB+n24F6+Mf4WFpVy6XLWNKCquVVcCW43wUAQT4GHgD/AiH8Et9eNd9l6Rtqa7ICm3EL7BTzUa7/kcgttAMj06TCzFQ5xYc9rPLcHnvMV2Bqcs52h3uhfp3yR7Y+NZjS12wOlC6G7Jfb8uh0N5OW8DalZhzkV02qPAPiP7aTu9IRRwX3nWX3DRo6c2JWYABgK2vVkNZhmrEzL5atUYrSfhnZ7d3hMMAgaY7+4wyHILcS3aBnVS+Zz5uQ5wnW77XqnsYwczmaRWvxsavPRLqkGsORy/ujj1/HM787Q9MLTMQfV+sZurtiygv31e2MPWuHQaLNsxGAPR37xmAZABSRvCkw8dJSvSp3IotiRv4bI/aE+/Edep6Wxg9z8XPvYhcbMT3aQTmJaxHuE/fUXGB4KxDVQ7PhzL3HLxoVQtSLm504p0d+gol+Fava0hqpSS8Q2Jw52j8/LagzDdQd1juB0YslGNRtNR9n96Mn4FyEY7OHM+SFWbb2Mk78+Pf54XW67qPhgjx24lKrcvOmcM1nNKF0GJellmUottxdKKm1QnWi4336JV6+B2k0aXCeqqJmbqUXhGdJARvH8nz22iw0rKy5O2FFYbgcsjQxDfunsjjNaLrfty88viT6bl0pr6ZIpSWuoKrWCAThzYPImt5w8ii+/mj/64BVUFXr58pfrL41m35w8+6UfGLL3p+or7W7NnJrV2YOMBDdC1XJ7viOK5sFAl23u726Z+XU8OXhXbmSkqS72KUIFpbgWrqBiQQXXb6lhXnEuP/vNYXp2vhjb66i0YuehpkcNVQOqSj1jwEzSbBkYovfkEX7ybAF93X0XpyhMd26PDRhj59jOZNE6PvDRGzha38mxl3cipZVTj47NlGAKEAry1FP7+cTvX8+Pv9M++jsJBkYHmnlzbW14unKFQ4z0dEefjCEnz7ZsjPihYjFl88vIy8/h4IkORkbC9DQ02P10ofKMlqx4KiI1wANAFTZ232eM+bep9teAqmYvv8TWopLRl9fbQcv2Z2KfrpEGyrbcxM03LOcX3/5x9E2mQ70889wxAv4AlevW0fHqtuie5/ZQeMV19Pf0Q/3e+A86iUz9XvYcr6Vm6xYatz8/GsQi0328vplHcYdDky+M7vXZJuOJ6+eOfb2zh7jQ18kF/8DFLR2hkfE/Vy6zQbs3lsm9KiWSlJfXEQS+aIzZIyJFwG4R2WaMOTLZztqHqmYvrzi5fXkZGEwBuurr+cUPfhlb/2PnOQb2v8LIqf10HDkc/aV3QRmr1y0kv3CalVxS3QdtDIe3baempoS7/+xuqFgMnU2jZRzojv53VVQBizeMfu6qVkb3GQz6oXLppMc2znAGtohcwsImvttMjDEtxpg9zv0+4CiwaKr9NaCq2bvQmPjmslR/+SdCR8PFNaZojfjt6jlRcK28inXveAtNjT0MHp1mUW63N/W/194Odtz/IHtOdLBk7bL42+qG+8krGzOoqOlwdOkqQ0E8uVEsutB34eL5r+qSJiLLgE3A61PtowFVpafShVmTHN+18irbLJ4k4WCQE4fO0r77tZkXZY8EsLyilAbXM/UdfOb21XaUbzxG/Awd3BF7og9PDuFQnHmFVVoygInzHzBPRHaNuX16svcQkULg58AXjDFTpiTTPlSVnjI5Of6EAS65+bkM+gpiH6HsybGBb6al8Br2E3Pv9XD/1LVDTw6UVNm5wUmqpfXu2c4TG6ptQE1kQv/cwhkzSIUvRDGaWGWUWXShnjfGbJluBxHxYoPpg8aYx6bbV2uoKrHcnsnXMS1dYOeTTlUrSnVTZDxELj7unDxyVm2y9wvLKb36bQwe3BHfBUJxJd5FCcwuNNZ030BF81h301Zci1Yn572d99/56nFWbb0svrmgpQsu3uZyM++q66BgmvmsLrftR1VZJVkLjIuIAD8AjhpjvjnT/hpQVWLl5Nm8uxMVVjB/5bLJm3HzihKblm6u5BZC9drxQdWECdQftPdDI3TXnYj/8rnzHCMnY08oP2teH+1tvYQ7mpL6NubsYT78lmW4ll0e+5MnG30bDnH+d8/bwU2TKV+Ea94izQedhZI1KAm4AfgE8HYR2efcbp9qZ23yVYk11Dd5k1vTYdqbj00+tcY/aAc2ZZqhPjuFY2zAHNuHOdQ3+ybTVIxwbj/NhY6G5GccH/Hz+GuNlFSU0FWHvTDxFdjpKjO991RTtKba7vZQsKSWgSO7bUAVSU75IheMurbqnLFZj5LzWTXGvAJE3XymAVUlTuRLKvLh9uTY/KyRoDDdl2CmfgFl4nFX1NigNV3f5RxlnjlT30ZJeTFdizfwPz73DjYvLOKnB9p4/Im9hOv3Je73awzBkSAUV7Lw8sv40NtX8Oyecxw/eMauZhSvihp7jMN99n+3185pTVVOYpVSGlBVYuTkkbdmE0PHdo/W0kIj41MKRnLPdjSk5BAVLL/9fbz3+qWIwH0/2cPggVdTdzBeH5u3LOPymlI23XUFD7zWxDf+36ssW72Ie+6+nv94LJ/+1jb7OTIGettjD1Q5eVC9mpLKMmqWzuNtd25lIBDicFMPBQU5fPwjV/PILzwEjk45E2Lca7mWrCccDsNA72iw72rOzAurLJImmQc1oKoEMWFcLtf4PtKJKQX9A4kdfCQC5YttLtjZrA2a6Xz5Uyewz8mDvGLc82u468NbGPAH+bd/fZyS2hV8+L1X8KMjO0dHJHt99gIo2X2MIvjWXcNn79pKnz/EYy+c4jtHT9gLrXCI+lN7+O6pLXzqY1ezriqPmsJ8QsbQPDDE3z98kPbtT8/8HrmFXPae29i0eh5etwuvW8hxu3joN0e58Mb2N8t89MobWbluMUdO7h0/l3rROq67dRPzS/LIy3Hj87op9Ll5ZmcTZ+tb8Qdyoevc1P21ak5pLl+VXUb8DOx/5eLtuYVQvsj2NQ73w+J1ULkEBvugtW527ykum6XJhBMXUD05ozWiTLFgJfR3jfZDe32QW8iqm9/ChpUVDAdCnG3pY9uOM7Tu3w9dzfT0ncfcuo6CDVcz0HDSNllW1toE9fEE1LF9kpGLpsl+hwWlfOizH6aqxMdjL52m8cVtF18IGEPw+E7+/e/227IUzweXmw3XX87nP7iBr+x4wQa/qfpB5y3hj+95L/Vtffz4e78cTZAxyf55hXmsry3nSORC0O1hzXvey00bq9n2+ll+t20/9F+A0irobNZUhGkqXf5cNaCq5PIPQNsp+4kvqeLKa1YTDIY5/Mq+2b92OGQz5ICdKhFJRB9vcM0vYf6Wa2lvbIeGfenzVzqTpiPgzUVWbKZ2zWK2bFhAMBymsb2fx//zOdskOTGTVWCIB37wDPfc827es/qd9PhH+Mn+Vl7d1UhfVx9DXT225j8yPJqnOdKsGQk+5YsgtwhPQT6+PB/+IT8ul4uSihJCwRDhcJj+7n6CPZ24ikqpXbOYO25Yys9fOs2ZF5+zxzTdoKtgwMn1awd2Hd1fRMW7Vtj3n1/LP/3VB3m5rotf3vufo+Urq+YbX/sI3/ttHSd/89T4ck9yPs8fO8aid69FatZhOlt5+0duId/n4fvf+y20148+J5YFDtScMiQ1l29MNKCq5DJm9EszGKClpZeudqc2VVZta0Sz/WNwe/AtW4t/0D9aS8sttK/f1RzdqFGgeN1GrtiwgOM+L41dLZnxJSqCq/ZKPv6Rq7llVRkP7Wpm2/Y6uvY5WZOm69trO8W9f3Mv95YvpnDJMq7YtIRPvHsNd6yuQgS6hgP0BoIcauvnwWdPce50C4UlhRSVFjIyEuTGrTWsqiyg/vwggWCYls5B2tv7aapvJVi3B4orKVi2kiUrN+DzeWht7uab//4iNB6Oq88x3HiUrqG3Ur75ev7qk5sIG7h+eSllf/4pHnj4NWitY8t/uYnhkTAnn346unSYFxr5/qN7+Nwf3EDn4Agnz/Xw/BM77EWgygxJWts0HpIukT2RXPnzjW/NR1J9GGoyLrf9MhWx6fhm0wclYjP69HbYAOofGF/jySuyTcKl82F4EJqPT/typVe/jRuvraV/eIQXf/iz9J6vKELuhmvZeu0KGpt6mD+/kDdeOmyb0acbuJNbaH9vU03nifxOyxay6rKlbFxTSU15Hjlu4VhrP8srC6jI97C3sY/eoQBHTpyn6VQztJ+2r2nCSf12K9lyM9++53r2t/bzjX98FFxuvv7XH2VleT5Pn+xk7fw8/uWnhy7uZ61eY0c1T9Vk68mxn0tPTnqf9ww2vO/e3TNlJYpHUc1ac9UX74/ruS/99xsSekxaQ1VzK1IzMWb2AzqMsU2SJjx5Wr+hPiiporismFCwgIER/7QjjPu7bVPx4HDw4uW80kVeEcXrN7NyTRUF+TkcPHCO7r07aIhlRZvp9jXG1sy7Wzl5ei8nf4kNsi73aEajsWuWJip4uj1Rzbnt2f0Sf/eLMk5s/92bLQhf+Z8/tBcK4ZCdxxpZDi7SZ+pyc9cnb+bxpw4xeGCKgBopjwbTjJQuFUMNqCqzzRSU20/T6/GRW1pMQfViBvwDo/2CE/4Ig60N/OqhNrtCTLotGef14V52Ode8ZQ0AO36z0zZLxtp0Gk//cqTZPlm/E5ebnNVX2QxTM02LMYbmxvPjl1br77x4Tq0nBxatfbN5+aGHXoULyc38pFInTeKpBlQ1x/JLoLjSToBvPpb8vwRn4NJwp7PUV0EZ+ZWVDPb02SbgsbW1uRjBGWsmnWVXctUNa6mZX0hxXg4PfOsnUydk8OTA/FobNHs77LzfyMAaT479vfd2JH6pvXiUVY8ubWfCBE4diHr9UZfLBfOWjF5MTVa7DY3YwVoutx0p3NNxaU+tynJaQ1WXJl8B7qJSQu2NiQmmIvbLORScfu3RSJNwbweDnU1QWUvums0Mn++w/X/Jnpjv9tjjLCjF5fUSHnGy6QQDcP7sxfu73JRfczPXXLWE7S+fZPfJg6OjbScqqwZfPuXLlrJqzXwO7Wtk4FDbaC28dAHkFlJUXU1ffSj1g618+Sy48kpaf/cS5DuJ7GNYPCAcDuMryOfNS6EFqy6+OHN7wevj2o++F6/HJhdpae2n7pXX7OckTb6A1ewZNKCqS4XXN74W2NVMKJFLs1Wvpaiqkr6zDdE/xz8I544w3F8FuYW4lm8kfPZI8vrPcvJs7XCoFy40Ena5wZvLwutuYvnyMvbubGCw7tCbQV9WbOaO921i/7EOfvPDJ2Zc3sw9r5pQKETnscO8fuoU9J0fvcAQAW8uDHTRt/dUegQS/yCtL2+zxzi2ideTYz8vM+Q/DgwHuPlta3lmyA/dbXD+jC2X22P7UufXctvtVxI2hv2H22jbvw8Gushds5m3fOgWDh9qoXP3K9P3JavMkUajfDWgquQRsZPyk5mazZtDKBiKvRkzMvjG5SY8Uh3d8yMjiWP96w2HbHKByO8gFIRQPy0vPEXLS24oq2bx9TfS2d6N2+OmuqaCp357iMDxXVH93kInd43f4PaMrslqTHqmegwFIa8Iz5K1BE/sAmPwrd7MXR/czH/8/b9PXe6CUjxeD/3DQWpWLqKxDtusW1LF6ms3ctW6+fQOBth/tJ2mN96wFxcFZWAMwyf383LDcaqvvpbaW97F6SMNcObAXJZaJYnWUFX2Myb5q8i0NzDYHIi/dhkOTX+MY7P+mLBtSow1eItMvzDAhUaaXumyTbcmzPH9O2Y3AMibO/0xFlfirlp6cSBOhPwSWxOP5gvO5ebGt67lxKJKml98GpfbxWUL8yefulJWDbmFuAuLGezoYNfvgvzJx6/m/xw9zaZbr+Oy2nJ+u72eh/9zO7SeHP/78w/YVoJwCPo7aX75OXJWbqR2/TIafT6CfU4SC82CpGZJA6rKbJNNl0mkonl2KsaFxvhXEKmosc+frolxuN+mZ0yEmQbfDPcTOnM46qkqMTFhGxDHltXltgEtNDJ++0A3L247xKZrVzF8zc1cflk1u872TT44Kb+EvPIyhk4egOF+Aj0lnGxdz2f/+G0EgmEeeWK/XcjdGDxrtrK4dgENzz1j3y9yixjxEzj6OqdP5TkJQBaCe76tzaZJTUfFJl1Omy4wrpInvyT+ZPhuj72lWt956GyaVZN16ZIa+8WdLgJD9jbbb6HJzu1Q3/jglZMHpQsoWLMR74orbbrCnLzRxxv2cfzwOXy5Pm67rJJX956b/Li6mhk6e2r0YmGwh1//xxPsO93J3rrzhEIhWzMH1l6+lH+6ayPumRYuDwzZWumZAzaFZbp8K6sY2dSD8dwSLQ2+sVTWikywj0duoU1+P8OAnKQbmzoxHm4PIpKeA2Bm26/tck//u8kthOpV0N/NQNMZm4A/t8iujuP2vDn4aLC5kcHmRg5eU8PZNyZphs4tdFoiJrRG9Hey89UT5Bbk4j916M1m7kM7DvL7h84SOu30jy5YaQNxqkc3q6Swo3xTfRSW1lBV8phZDL8b6E59ME0EcVFaXmibjbPNTBcaw/1wep+tBfa02Sbz7hbbjF5QZoOqMXbakDE8u6Ph4n5Mr4/q628azdI08RC62u3C4UXzqLj2bfY1m48z0tZk3weguxVf9TLbD6uyj0FrqEpdEkIjNJ9ph2Aa1lDngjE2sJYuwFe9DP+ZY+SWlzN85oQNqkO9b/ZRt7/cdPEF2Iif5kNHp65ND/cTDAQhGODC8eNQWAGF5bapPvKc4X78J/bY/l2VldKlhqoBVamZFFXY/rZ4BiUVlOH2zNA0eikY7LEJGZautdOc/IO2LzUyWAmm/lZsPz316waGCPX3jiaGiATTiTXddMgOpbKeBlSlZpJXbGs78QTU/BIGW50l5GJRusC+ZzZM5SgohdwiRhpPjK4OFJzFVKexQiPja57Z0E2gYqbzUJXKFNPVkOJ5bmT1lulqrU7SiYwkQt5l1zHU2gznzyALlnPl1hUc3neGkaa6xAa9UHB250dlhTSJp8kblCQiNSLygogcEZHDIvJ5Z3u5iGwTkZPO/2XOdhGRb4lInYgcEJHNY17rbmf/kyJyd7KOWc2BeKfRpDOR2KbFiCu6YJns/MLJYgxDA0P29+LJIa8gj66uYUYCIzaBwsSpM7OVqb8nlRCRXL7ZPigpCHzRGLNHRIqA3SKyDfgD4DljzD+IyJeALwF/AbwbWOXcrgG+C1wjIuXA3wBbsL+73SLypDGmK4nHrpJlpkXFI2tYzqSg1I7i7G5N7Soi5YtAXHjmLSBYtze6vtJwKPuDQP3eN+8OHt/LmTNF0Ntuz21Xc/pUKVTmS6NcvkmroRpjWowxe5z7fcBRYBFwB/AjZ7cfAe937t8BPGCs14BSEVkIvAvYZozpdILoNuC2ZB23SpJIkoaZMhtNMT3iIr5YH1VhAAASV0lEQVQC8irn2z65VNd6+85HH0yn48mxFwrZxj84foWXdPn2U1njUqihvklElgGbgNeBKmNMi/NQK1Dl3F8EjE2q2uRsm2q7yiTRfplGmwCh8xxD6bAQeCTrUCKUVc8uGYZSKqWSHlBFpBD4OfAFY0yvjKlNGGOMiCTkMkFEPg18GgBvGqV5U1YymjhTHUzBjjKNtpl6JslcFcblzv5mZnXJSpdGj6RmShIRLzaYPmiMeczZ3OY05eL83+5sPwfUjHn6YmfbVNvHMcbcZ4zZYozZIp4EDnhQajpDfVP/Nbs9dl5kXhHMr7X3UzFy1+2BpVekR25kpRIufXL5JnOUrwA/AI4aY7455qEngchI3buBJ8Zs/6Qz2vdaoMdpGn4auFVEypwRwbc625RKb1UrbBagoT47urW4EtfyjXZg1lzSqSUqi10qo3xvAD4BHBSRfc62vwT+AXhURD4FnAE+4jz2FHA7UAcMAn8IYIzpFJG/BXY6+/0vY4zO3r7UxdKEWVgOFYuhtS7+JdjiMeKHQWdE80A3DHQTLqp4c93TSHJ4wDYbl1SNH7yTSNONrFYqk6XRKN+kBVRjzCvAVMMv3zHJ/ga4Z4rXuh+4P3FHpzKeJ8dm24kmqM6rYfn6pTR4vYTrdsf2PiJ2ek4wEPvaq+EQrkWrx79n3wV7q6ixwT0cAq8PatbDueOj75nobwgRW1tO5RQjpZIkXTIl6WozKjMFhqILpp4cispLuGp9Fde9da3tz4yFMXbB63hG33Y2ER7sn7zf9MKYgevBAHS22MTuANVrba060bQPVamk0oCqsltohMBwgKYLA4RCBqqWx/Z8l9sOKIp2fuxEvvzJn1tRAzUb7OsbYxNUDHTZ/tXQSOKm4kQYo82+KmtFVoqM9ZZoesmqspsx+C90cHCflxxfDgzGWNMMh+Dswfj++oyB4QGKL7+a3iN7RvtMnZR89F24+L2G+zXwKRUDA4TD6dHkqwFVZb+2U/SfP2Nrg9EmjgAb9Irm2UWx49Vygt5gAE/NWoInd41eGredGr9fXhFlV2yla9f29Jhfq1SmSKNBSdrkqy4NoWBswRRsAE7EAKHuFsLhsG3inarpeOFK3v/OtVC5bHbvpdQlKF2mzWhAVWoqgSHovCiHSOxG/ITr90HveahePfnKNI1HePDnu8mdV5n63MRKZZh06UPVgKrUXAiH7MCjs4cmn7oy4idw7A2Ky4ptM7NSKkrpkylJ+1CVSia3J7o+0cJyKF1Ijs9j10tNFM3hq9Sc0RqqUslStQLf2i3R7TvQBf4BmrY/b7MlJYLLjW/t1tTkD1ZqjtjUg+nR5Ks1VJVeipzkBhOnlIxVVg3DfeNT90VjrmprIpRuvZkVq+bT3NxHS0etXUlmur9gY+JfbWZ+rX2+r2A0CUVvB5gwHq8Hf0mVXdRbqWxkNFOSUpPzD06f4s/tIW/xMny1G2JLMi9im1WTMeDH5R5fC/QV8JkPXsEtly+g50IPuL02+1Gyaoq9HXYAVdc5O8XHP2C3G8PA+Qtzn4xfqTmmfahKTWamDEHGEA6H8eX58PsKos+va4wNPLHILbR9oDMlWvD6bNAc7AFPDvmrr6S1L8CJ5h4Gj++1FwnJyM8bMdw/dY7ejgabiD9WFTXQ2ZQ+E/yUmk6afEw1oKrMEg7hP30Uf0GZ7WssKE1eZqFoE8kHA1BQZvd3uRk8U8ejTxqCI8HR1W3iCUyeHJi3xCaBiOf5eUW2iTsQ4/xbsE3EGkxVhkiXJl8NqCrzDPbYW7rM1wwF7XzVsmobwIZ6GTq4Y/J9Y+nHNWHIL4IFq6DlRPTHk5NncwiP+ONfrk5HBisVM+1DVZkr3RK+BwZtrbKkavJgLxJbf2YoCGcOQv80A7Qm43JDyYLRvlSlspz2oSqVbZxFxKdkDPR3xvaaoeD0I54nM9wPF87a+axGa5oquyUrOMZDa6hKFZTGvzxbolUuS8haqMXrN2teYHXJSJcaqgZUpYxJn8W3XW6oWDzrl+k9+AbkFkyeN1ipbGPivCWYBlSVvYoq7OCcmQz2xD94J8KXP7ug7PVBRQ055fMSk5A/MARt9VC6YPavpVSa0xqqUsnW3zn7QBmtytr4a4P5JRRs2AoFJQTq9sXeZzqV4X5orUvMaymlZpQm7VxKJcFcDlQ4dyT694tkbcrJh4IScouL8OZ4ob8r9jVblVJpMyhJA6pSsRCB0oUX5xKO5Q86t5DCletxuVz0dvUy3HCcYbcXetsTf7xKZbsk5vIVkfuB9wDtxpjLZtpfm3yVioW4oKwKCivif43AEP1t7Qz2D0Jfl22a7mnTzERKxSt5g5J+CNwW7WFoDVWpiGjWLg2H4PS+2b1PKAjNxwlOlvxB1y9VKiaG5M1DNcZsF5Fl0e6vNVSlwDblLlprR9vOJN7FFCcOWpr4OvklULUi9tdV6hKXLqN8tYaqFNjA1nrKJrpPlpmS7UdyFCulYjKL4DhPRHaN+fk+Y8x98b6YBlSlImZaOk4plW3OG2O2JOrFNKAqpZTKXEkc5Rsr7UNV2SVdlnRTSs2dJI3yFZGHgd8Ba0SkSUQ+Nd3+WkNV2WX+cjufc+wcUaVUVkviKN+PxbK/BlSVXfovaF+oUpeQZE6biZUGVJVd0mnBcaXUnEiXgKp9qErNFbfHJm5QSmUlraEqpZTKXEla2zQeGlCVSqbiSpsBqe+8NkcrlSTp0uSrAVWpZBGx6QZzcnXUsVJJpAFVqWxnDLSfhrwiTXivVBJpQFXqUqG1U6WSJp2mzegoX6WUUioBtIaqlFIqc+koX6WUUiox0qXJVwOqUkqpjKYBVSmllEoADahKKaXUrOkoX6WUUiqraA1VKaVU5tJRvkoppdTsGbQPVSmllEoIDahKKaVUIqRJQE3aoCQRqRGRF0TkiIgcFpHPO9u/KiLnRGSfc7t9zHO+LCJ1InJcRN41ZvttzrY6EflSso5ZKaVUhjEGTDi+W4Ils4YaBL5ojNkjIkXAbhHZ5jz2r8aYfxm7s4isB+4ENgDVwLMistp5+F7gnUATsFNEnjTGHEnisSullMoUaVJDTVpANca0AC3O/T4ROQosmuYpdwCPGGP8wGkRqQOudh6rM8bUA4jII86+GlCVUkqljTmZhyoiy4BNwOvOps+JyAERuV9Eypxti4DGMU9rcrZNtV0ppZRKmybfpAdUESkEfg58wRjTC3wXWAFsxNZgv5Gg9/m0iOwSkV0mOJSIl1RKKZUJjInvlmBJHeUrIl5sMH3QGPMYgDGmbczj3wd+5fx4DqgZ8/TFzjam2f4mY8x9wH0Arvz56dGgrpRSKslMUmqb8UjmKF8BfgAcNcZ8c8z2hWN2+wBwyLn/JHCniPhEpBZYBbwB7ARWiUitiORgBy49mazjVkoplUEMadPkm8wa6g3AJ4CDIrLP2faXwMdEZCP219AAfAbAGHNYRB7FDjYKAvcYY0IAIvI54GnADdxvjDmcxONWSimVSS6BUb6vADLJQ09N85yvA1+fZPtT0z1PKaWUSjXNlKSUUiqDpU8fqgZUpZRSmS3bm3yVUkqpOaE1VKWUUmqWkjSnNB4aUJVSSmW2NKmhzknqQaWUUirbaQ1VKaVUZtMmX6WUUmq2dNqMUkoplRhaQ1VKKaVmKZLLNw1oQFVKKZXZwlpDVUoppWYpffpQddqMUkoplQBaQ1VKKZXZ0qSGqgFVKaVU5jLoKF+llFJq9tKnD1UDqlJKqcymNVSllFIqAdKkhqqjfJVSSqkE0BqqUkqpDKbroSqllFKzp6kHlVJKqQTRGqpSSik1WzptRimllEqMNKmh6ihfpZRSagoicpuIHBeROhH50nT7ag1VKaVUZktSk6+IuIF7gXcCTcBOEXnSGHNksv21hqqUUipzRXL5xnOb2dVAnTGm3hgTAB4B7phqZ62hKqWUymBJHZS0CGgc83MTcM1UO2tAVUoplbHMUMfTw/vunRfn03NFZNeYn+8zxtwX77FoQFVKKZWxjDG3JfHlzwE1Y35e7GyblPahKqWUUpPbCawSkVoRyQHuBJ6cametoSqllFKTMMYEReRzwNOAG7jfGHN4qv01oCqllFJTMMY8BTwVzb7a5KuUUkolgAZUpZRSKgE0oCqllFIJoAFVKaWUSgANqEoppVQCaEBVSimlEkADqlJKKZUAGlCVUkqpBNCAqpRSSiWABlSllFIqATSgKqWUUgmgAVUppZRKAA2oSimlVAKIMSbVx5BwItIBDADnU30sc2AeWs5souXMLlrOUUuNMZVzcTCpkpUBFUBEdhljtqT6OJJNy5ldtJzZRct5adEmX6WUUioBNKAqpZRSCZDNAfW+VB/AHNFyZhctZ3bRcl5CsrYPVSmllJpL2VxDVUoppeZM1gVUEblNRI6LSJ2IfCnVxzNbItIgIgdFZJ+I7HK2lYvINhE56fxf5mwXEfmWU/YDIrI5tUc/NRG5X0TaReTQmG0xl0tE7nb2Pykid6eiLNOZopxfFZFzzjndJyK3j3nsy045j4vIu8ZsT+vPtYjUiMgLInJERA6LyOed7Vl1TqcpZ1adUxHJFZE3RGS/U86vOdtrReR155h/IiI5znaf83Od8/iyMa81afmzkjEma26AGzgFLAdygP3A+lQf1yzL1ADMm7Dtn4AvOfe/BPyjc/924DeAANcCr6f6+Kcp103AZuBQvOUCyoF65/8y535ZqssWRTm/CvzZJPuudz6zPqDW+Sy7M+FzDSwENjv3i4ATTnmy6pxOU86sOqfOeSl07nuB153z9Chwp7P9e8B/de7/CfA95/6dwE+mK3+qy5esW7bVUK8G6owx9caYAPAIcEeKjykZ7gB+5Nz/EfD+MdsfMNZrQKmILEzFAc7EGLMd6JywOdZyvQvYZozpNMZ0AduA25J/9NGbopxTuQN4xBjjN8acBuqwn+m0/1wbY1qMMXuc+33AUWARWXZOpynnVDLynDrnpd/50evcDPB24GfO9onnM3Kefwa8Q0SEqcuflbItoC4CGsf83MT0H/ZMYIBnRGS3iHza2VZljGlx7rcCVc79TC9/rOXK5PJ+zmnqvD/SDEqWlNNp7tuErdVk7TmdUE7IsnMqIm4R2Qe0Yy9sTgHdxpigs8vYY36zPM7jPUAFGVDORMq2gJqNbjTGbAbeDdwjIjeNfdDYdpWsG6qdreVyfBdYAWwEWoBvpPZwEkdECoGfA18wxvSOfSybzukk5cy6c2qMCRljNgKLsbXKtSk+pLSXbQH1HFAz5ufFzraMZYw55/zfDvwC+8FuizTlOv+3O7tnevljLVdGltcY0+Z8WYWB7zPaBJbR5RQRLzbIPGiMeczZnHXndLJyZus5BTDGdAMvANdhm+Y9zkNjj/nN8jiPlwAXyKByJkK2BdSdwCpnJFoOtnP8yRQfU9xEpEBEiiL3gVuBQ9gyRUY/3g084dx/EvikM4LyWqBnTHNbJoi1XE8Dt4pImdPEdquzLa1N6Nf+APacgi3nnc6IyVpgFfAGGfC5dvrLfgAcNcZ8c8xDWXVOpypntp1TEakUkVLnfh7wTmx/8QvAh5zdJp7PyHn+EPC80yIxVfmzU6pHRSX6hh09eALb3v+VVB/PLMuyHDtCbj9wOFIebN/Ec8BJ4Fmg3NkuwL1O2Q8CW1JdhmnK9jC2aWwE26/yqXjKBfwRdqBDHfCHqS5XlOX8sVOOA9gvnIVj9v+KU87jwLvHbE/rzzVwI7Y59wCwz7ndnm3ndJpyZtU5Ba4A9jrlOQT8tbN9OTYg1gE/BXzO9lzn5zrn8eUzlT8bb5opSSmllEqAbGvyVUoppVJCA6pSSimVABpQlVJKqQTQgKqUUkolgAZUpZRSKgE0oCqVwURkmYxZyUYplToaUJVSSqkE0ICqVIo4tcujIvJ9Z83JZ0QkT0RWisizzlqUe0RkhZNR6J9F5JDY9XE/murjV0qNpwFVqdRaBdxrjNkAdAO/BzzobLsSuB6baemD2MTrVwK3AP+crkvzKXWp0oCqVGqdNsbsc+7vxi7CvMgY8wsAY8ywMWYQm/LuYWMTsLcBLwFbU3LESqlJaUBVKrX8Y+6HgNJUHYhSanY0oCqVXvqAJhF5P4CzSkc+8DLwUWfR50rgJrJ51Q6lMpAGVKXSzyeAPxWRA8AOYAF2LdwD2JWHngf+3BjTmrpDVEpNpKvNKKWUUgmgNVSllFIqATSgKqWUUgmgAVUppZRKAA2oSimlVAJoQFVKKaUSQAOqUkoplQAaUJVSSqkE0ICqlFJKJcD/B/GhAawRq5wkAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Plot result\n", "fig, ax = plt.subplots(figsize=(8, 8))\n", "wat_occ = ax.imshow(sentinel[0], cmap='Blues_r')\n", "fig.colorbar(wat_occ, ax=ax, anchor=(0, 0.3), shrink=0.7)\n", "\n", "ax.set_title('Water Occurrence')\n", "ax.set_xlabel('ncol')\n", "ax.set_ylabel('nrow')\n", "\n", "plt.show();" ] }, { "cell_type": "markdown", "metadata": { "id": "KB3RnqOrGFK2" }, "source": [ "*Export result as GeoTiff*" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "id": "WUyC2efFVni2" }, "outputs": [], "source": [ "# Data dir\n", "data_dir = \"output\"\n", "\n", "# Output raster\n", "out_tif = os.path.join(data_dir, \"sentinel_water_occurrence.tif\")\n", "\n", "# Copy the metadata\n", "out_meta = image_source[8].meta.copy()\n", "out_meta\n", "\n", "# Update the metadata\n", "out_meta.update({'driver': 'GTiff',\n", " 'dtype': 'uint16',\n", " 'nodata': 0.0,\n", " 'width': sentinel.shape[2],\n", " 'height': sentinel.shape[1],\n", " 'crs': image_source[8].crs,\n", " 'count':1,\n", " 'transform': image_source[8].transform\n", " })\n", "\n", "with rasterio.open(out_tif, \"w\", **out_meta) as dest:\n", " dest.write(sentinel.astype('uint16'))" ] } ], "metadata": { "colab": { "collapsed_sections": [], "name": "Water-surface.ipynb", "provenance": [] }, "kernelspec": { "display_name": "Python 3.7.10", "language": "python", "name": "python3" }, "language_info": { "name": "python", "version": "3.7.10" }, "vscode": { "interpreter": { "hash": "54697fe90b6b48df13497a1835b00f8636d17074f5f61d6ce3a5a8b26dd66d49" } } }, "nbformat": 4, "nbformat_minor": 0 }