{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Downloading with ASF" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import terragon\n", "import geopandas as gpd\n", "from pathlib import Path\n", "from utils import visualize_sat_images, list_of_dicts_to_tree" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# we will use the sample data to download a minicube\n", "gdf = gpd.read_file(Path(\"data/TUM_OTN.geojson\"))\n", "gdf = gdf.to_crs(\"EPSG:32632\")\n", "arguments = dict(\n", " shp=gdf,\n", " collection=\"SENTINEL-1\",\n", " start_date=\"2021-01-01\",\n", " end_date=\"2021-01-02\",\n", " bands=[\"VH\"],\n", " resolution=10,\n", " download_folder=\"tests/download/\",\n", " filter={\"processingLevel\": \"GRD_HD\"},\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# credentials for an NASA earth data account can be created here: https://urs.earthdata.nasa.gov/users/new\n", "credentials = {\"asf_username\": \"\", \"asf_password\": \"\"}\n", "tg = terragon.init(\"asf\", credentials=credentials)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Finding collections" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "['AIRSAR', 'ALOS', 'ERS', 'ERS-1', 'ERS-2', 'JERS-1', 'NISAR', 'RADARSAT-1', 'SEASAT 1', 'SENTINEL-1', 'Sentinel-1A', 'Sentinel-1B', 'Sentinel-1C', 'SIR-C', 'SMAP', 'UAVSAR']\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/mnt/datastorage/home/adrianhohl/terragon/terragon/alaska_satellite_facility.py:75: UserWarning: Currently Terragon only supports the following collections: ['SENTINEL-1', 'ALOS PALSAR', 'ALOS AVNIR-2']\n", " warnings.warn(\n" ] } ], "source": [ "# explore the ASF data on their website: https://docs.asf.alaska.edu/datasets/using_ASF_data/\n", "# collections can be viewed using the following command\n", "collections = tg.retrieve_collections()\n", "print(collections)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "WARNING:asf_search:WKT REPAIR/VALIDATION: The following repairs were performedon the provided AOI:\n", "[\"'type': 'GEOMETRY_SIMPLIFICATION': 'report': 'Shape Collapsed to Point: shape of 4 simplified to 1 with proximity threshold of 0.004'\"]\n", "WARNING:asf_search:Query Time Elapsed 1.7983815670013428\n", "WARNING:asf_search:Product Subclassing Time 0.001031637191772461\n", "WARNING:asf_search:Page Processing Time 1.1920928955078125e-05\n", "WARNING:asf_search:Page Time Elapsed 1.8896667957305908\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset> Size: 2kB\n",
       "Dimensions:      (x: 28, y: 18, time: 2)\n",
       "Coordinates:\n",
       "  * x            (x) float64 224B 6.976e+05 6.976e+05 ... 6.979e+05 6.979e+05\n",
       "  * y            (y) float64 144B 5.326e+06 5.326e+06 ... 5.326e+06 5.326e+06\n",
       "  * time         (time) datetime64[ns] 16B 2021-01-01T05:26:52 2021-01-02T05:...\n",
       "    spatial_ref  int64 8B 0\n",
       "Data variables:\n",
       "    vh           (time, y, x) uint16 2kB 65535 65535 65535 ... 65535 65535 65535\n",
       "Attributes:\n",
       "    crs:          EPSG:32632\n",
       "    data_source:  ASF\n",
       "    collection:   SENTINEL-1
" ], "text/plain": [ " Size: 2kB\n", "Dimensions: (x: 28, y: 18, time: 2)\n", "Coordinates:\n", " * x (x) float64 224B 6.976e+05 6.976e+05 ... 6.979e+05 6.979e+05\n", " * y (y) float64 144B 5.326e+06 5.326e+06 ... 5.326e+06 5.326e+06\n", " * time (time) datetime64[ns] 16B 2021-01-01T05:26:52 2021-01-02T05:...\n", " spatial_ref int64 8B 0\n", "Data variables:\n", " vh (time, y, x) uint16 2kB 65535 65535 65535 ... 65535 65535 65535\n", "Attributes:\n", " crs: EPSG:32632\n", " data_source: ASF\n", " collection: SENTINEL-1" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# we can now download the data, in the background we search the folder of the corresponding tile for the correct files\n", "# the download may take a while depending on the size of the data and the number of files because ASF is not providing a STAC API\n", "# therefore it is possible to download data with up to 4 concurrent workers\n", "# make sure to have the right packages installed (remotezip)\n", "ds = tg.create(**arguments, num_workers=4)\n", "# or\n", "# items = tg.search(**arguments, num_workers=4)\n", "# ds = tg.download(items)\n", "ds" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA9IAAAHqCAYAAADoJXY3AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAACiYklEQVR4nOzdeXxU9b3/8deZNZOV7AlbWEUxoAjIprKETcUNK95b22qv7fV3a/2VW7zVaqu016X1V5cWW27tpWJFi61KW6vssoogBNCACoLsAgGyJ5NZz++PhGgEJDPJZGaS9/PxmIdk5nzP+UxM8p3P+X7O5ximaZqIiIiIiIiISItYoh2AiIiIiIiISDxRIi0iIiIiIiISAiXSIiIiIiIiIiFQIi0iIiIiIiISAiXSIiIiIiIiIiFQIi0iIiIiIiISAiXSIiIiIiIiIiFQIi0iIiIiIiISAiXSIiISEWvWrKGioiLaYYiIiMhXOHz4MFu3bo12GHFHibSIiLS5VatWMW7cOB566KFohyIiIiLnYJomt912G2PGjOHw4cPRDieuGKZpmtEOQkREOo5gMMjw4cMZPHgwCxcu5IMPPqB///7RDktERES+5B//+Aff+c53GDVqFBkZGTz//PPRDiluKJEWEZE29dJLL/HjH/+YXbt28YMf/ICysjJeffXVaIclIiIiX+D3+xk0aBDf//73ueaaaxg4cCAbN27kkksuiXZocUGJtIiItJn6+nouvPBCfv7zn/Otb32Lo0eP0r9/f5YvX86oUaOiHZ6IiIg0+v3vf8+TTz7Jzp07sdvtzJo1i507d7JkyZJohxYXlEiLiEib+dWvfsWCBQsoLi7GarUC8PDDD7NixQrWr1+PYRhRjlBERESqq6vp378/v/vd75g+fToAZWVl9O3bl7/85S9MmjQpyhHGPiXSIiLSJs41AdfU1NCvX79mk7WIiIhEz8MPP8zKlStZt25ds5PcZzshLmenRFpERNrED3/4Qz788MOzloR9uXxMREREouOzzz7jggsuOOtlV6cv0frZz37G7bffHqUI44MSaRERabV9+/YxcOBANm3axODBg8943e/3M3jwYO6++27uvvvuKEQoIiIiAP/+7/9OeXk5f/3rX8/6+ssvv8x9993H7t27cblc7Rxd/FAiLSIirfav//qvuFwu/vjHP55zmzfeeIM777yTPXv2kJqa2o7RiYiICMDOnTsZNmwYJSUl9OvX76zbBINBLr/8cr72ta9x//33t3OE8UOJtIiItMrmzZsZN24cu3btonv37ufczjRNxo0bxxVXXMGjjz7ajhGKiIgIwLRp0+jTpw+/+c1vvnK7VatWceONN7Jnzx6ys7PbKbr4okRaRETCFmpy3NKkW0RERNrW6eR47969ZGVlnXf7libdnZUSaRERCVs45dotKQMXERGRthNOufaHH37I0KFDv7IMvDNTIi0iImHx+/0MGjSI73//+yE1EDtfYzIRERFpWy+//DL3338/u3btCqmB2Pkak3VmSqRFRCQsrbml1axZs9i5c+dZb5UlIiIibef0La1+/vOf861vfSuksUePHqV///5nvVVWZ6dEWkREQlZdXU3//v353e9+x/Tp00MeX1ZWRt++ffnLX/7CpEmTIhChiIiIAPzqV79iwYIFFBcXY7VaQx7/8MMPs2LFCtavX49hGBGIMD4pkRYRkZDNnj2bFStWsG7durAn1dMT+9atW7FYLG0coYiIiLTFieuamhr69evH3Llzuemmm9o4wvilRFpERELSVmVerSk1ExERkfNrq0upWnM5V0elRFpERELSlo1Hwm1+IiIiIl+tLZt7+v1+Bg8ezN133x1Sg9GOTIm0iIi02M6dOxk2bFib3QojnNtxiIiIyPm19e0mw7nlZUemRFpERFrsuuuuo3fv3vzmN79ps32uWrWKG2+8kb1795KVldVm+xUREemsNm/ezLhx49i1axfdu3dvk32apsm4ceO48soreeSRR9pkn/FMibSIiLTI6tWrueGGGyKS8E6bNo2+ffvy61//uk33KyIi0tmYpsn48eMZM2YMjz76aJvu+3SCvnv3brp169am+443SqRFRKRFvvWtb5GXl8cTTzzR5vveuHEjEydOpKamps33LSIi0pl8+umnXHDBBZSVlUWkBHvq1KlMnDiRe++9t833HU90vxEREWmRQCBAZmZmRPadmZlJIBCIyL5FREQ6k0AggN1uj9h1zJqzG9iiHYCIiMQP0zSJRCGTiqNERETalubWyNKKtIiIiIiIiEgIlEiLiEhITq9Kt/WjNR5//HEMw2DmzJlNz91xxx0YhtHsMXLkyGbjPB4P99xzD1lZWSQlJXH99ddz+PDhVsUiIiISCyI1X2ulu4ESaRERabFYnJQ3b97Mc889x+DBg894berUqRw9erTp8dZbbzV7febMmSxatIiFCxeyfv16ampqmDZtmq79EhGRuBdr83VHo0RaRETiVk1NDbfddht/+MMfSE9PP+N1p9NJXl5e0yMjI6PptcrKSubNm8eTTz7JxIkTGTJkCAsWLKCkpIQVK1a059sQERGROKNEWkREQhJLZ7jvvvturr32WiZOnHjW11evXk1OTg4XXHAB3/3udyktLW16rbi4GJ/Px+TJk5ue69q1K4WFhWzYsCGseERERGJFLFaRdSTq2i0iIjGjqqqq2ddOpxOn03nWbRcuXMjWrVvZvHnzWV+/+uqrueWWWygoKGDfvn389Kc/ZcKECRQXF+N0Ojl27BgOh+OMlezc3FyOHTvWNm9IREREOiQl0iIi0mKROhNtmiZ+v5+0tLRmzz/88MPMnj37jO0PHTrED37wA5YtW0ZCQsJZ93nrrbc2/buwsJBhw4ZRUFDAm2++yfTp078yFsMwwnsjIiIiMUIrx5GlRFpEREISqUTaZrNx6tSpZs+fazW6uLiY0tJShg4d2vRcIBBg7dq1PPvss3g8HqxWa7Mx+fn5FBQU8MknnwCQl5eH1+ulvLy82ap0aWkpo0ePbqu3JiIiEhVKpCNL10iLiEjMSE1NbfY4VyJdVFRESUkJ27dvb3oMGzaM2267je3bt5+RRAOcOnWKQ4cOkZ+fD8DQoUOx2+0sX768aZujR4+yY8cOJdIiIiLylbQiLSIiLRYrTUZSUlIoLCxs9lxSUhKZmZkUFhZSU1PD7Nmzufnmm8nPz2f//v088MADZGVlcdNNNwGQlpbGnXfeyaxZs8jMzCQjI4N7772XQYMGnbN5mYiISLyI1HwdC58DYoESaRER6XCsVislJSX86U9/oqKigvz8fMaPH88rr7xCSkpK03ZPP/00NpuNGTNm4Ha7KSoqYv78+Wdd0RYRERE5TYm0iIi0WCSbjbXW6tWrm/7tcrlYunTpecckJCQwZ84c5syZ0+rji4iIxBKtHEeWEmkREQlJrCbSIiIi8jnNrZGlZmMiIiIiIiIiIdCKtIiItFgsl3aLiIjI59RsLLK0Ii0iIiIiIiISAq1Ii4hISLQiLSIiEvs0t0aWEmkREWkxlXaLiIjEB82tkaXSbhEREREREZEQaEVaRERCohVpERGR2KdmY5GlFWkRERERERGREGhFWkREWixS10iLiIhI29J8HVlKpEVEJCQq7RYREYl9mlsjS6XdIiIiIiIiIiHQirSIiLSYbn8lIiISH9RsLLK0Ii0iIiIiIiISAiXSIiISktOr0m39EBERkbYTi/P1448/jmEYzJw5s+m5O+64A8Mwmj1GjhzZbJzH4+Gee+4hKyuLpKQkrr/+eg4fPtyqWFpLibSIiLRYpCZlJdIiIiJtK9bm682bN/Pcc88xePDgM16bOnUqR48ebXq89dZbzV6fOXMmixYtYuHChaxfv56amhqmTZtGIBAIO57WUiItIiIiIiIiEVNTU8Ntt93GH/7wB9LT08943el0kpeX1/TIyMhoeq2yspJ58+bx5JNPMnHiRIYMGcKCBQsoKSlhxYoV7fk2mlEiLSIiIYm1M9wiIiJypkhWkfn9fqqqqpo9PB7POWO5++67ufbaa5k4ceJZX1+9ejU5OTlccMEFfPe736W0tLTpteLiYnw+H5MnT256rmvXrhQWFrJhw4a2+4aFSIm0iIiIiIiItNiqVatIS0tr9nj88cfPuu3ChQvZunXrOV+/+uqreemll3j77bd58skn2bx5MxMmTGhKzI8dO4bD4ThjJTs3N5djx4617RsLgW5/JSIiLRap1WOtSIuIiLStSM6t48eP5/XXX2/2nNPpPGO7Q4cO8YMf/IBly5aRkJBw1n3deuutTf8uLCxk2LBhFBQU8OabbzJ9+vRzxmCaJoZhhPkOWk+JtIiItJjKsEVEROJDpOZr0zSx2Wykpqaed9vi4mJKS0sZOnRo03OBQIC1a9fy7LPP4vF4sFqtzcbk5+dTUFDAJ598AkBeXh5er5fy8vJmq9KlpaWMHj26jd5V6FTaLSIiIiIiIm2uqKiIkpIStm/f3vQYNmwYt912G9u3bz8jiQY4deoUhw4dIj8/H4ChQ4dit9tZvnx50zZHjx5lx44dUU2ktSItIiIhUWm3iIhI7IuFuTUlJYXCwsJmzyUlJZGZmUlhYSE1NTXMnj2bm2++mfz8fPbv388DDzxAVlYWN910EwBpaWnceeedzJo1i8zMTDIyMrj33nsZNGjQOZuXtQcl0iIiIiIiItLurFYrJSUl/OlPf6KiooL8/HzGjx/PK6+8QkpKStN2Tz/9NDabjRkzZuB2uykqKmL+/PlnXdFuL0qkRUSkxdRsTEREJD7E6ty6evXqpn+7XC6WLl163jEJCQnMmTOHOXPmRDCy0CiRFhGRkCiRFhERiX2RbDYmajYmIiIiIiIiEhKtSIuISIuptFtERCQ+aG6NLK1Ii4iIiIiIiIRAK9IiIhISrUiLiIjEPs2tkaVEWkREWkyl3SIiIrEvUvP16X2LSrtFREREREREQqIVaRERCYnORIuIiMQ+zdeRpRVpERERERERkRBoRVpERFpM10iLiIjEB82tkaVEWkREQqJEWkREJPap2VhkqbRbREREREREJARakRYRkRZTabeIiEh80NwaWVqRFhEREREREQmBVqRFRCQkWpEWERGJfZpbI0uJtIiItJhKu0VEROKDmo1Flkq7RUREREREREKgFWkREQmJVqRFRERin+bWyNKKtIiIxL3HH38cwzCYOXNm03OmaTJ79my6du2Ky+Vi3Lhx7Ny5s9k4j8fDPffcQ1ZWFklJSVx//fUcPny4naMXERGReKNEWkREWuz0NdJt/WiNzZs389xzzzF48OBmzz/xxBM89dRTPPvss2zevJm8vDwmTZpEdXV10zYzZ85k0aJFLFy4kPXr11NTU8O0adMIBAKtiklERCTaYnHO7kiUSIuISItFalIOd2Kuqanhtttu4w9/+APp6enN4nzmmWd48MEHmT59OoWFhbzwwgvU1dXx8ssvA1BZWcm8efN48sknmThxIkOGDGHBggWUlJSwYsWKNvl+iYiIREuszdkdjRJpERGJW3fffTfXXnstEydObPb8vn37OHbsGJMnT256zul0MnbsWDZs2ABAcXExPp+v2TZdu3alsLCwaRsRERGRs1GzMRERCUkkm41VVVU1e97pdOJ0Os86ZuHChWzdupXNmzef8dqxY8cAyM3NbfZ8bm4uBw4caNrG4XA0W8k+vc3p8SIiIvFKK8eRpRVpERGJCX6/n7S0tGaPxx9//KzbHjp0iB/84AcsWLCAhISEc+7TMIxmX5umecZzX9aSbURERKRz04q0iIi0WKSujTJNE5vNxqlTp5o9f67V6OLiYkpLSxk6dGjTc4FAgLVr1/Lss8+ya9cuoGHVOT8/v2mb0tLSplXqvLw8vF4v5eXlzValS0tLGT16dJu9NxERkWiI1Iq0VrobaEVaRERCEsnGJampqc0e50qki4qKKCkpYfv27U2PYcOGcdttt7F9+3b69OlDXl4ey5cvbxrj9XpZs2ZNU5I8dOhQ7HZ7s22OHj3Kjh07lEiLiEjcU6OxyNKKtIiIxJ2UlBQKCwubPZeUlERmZmbT8zNnzuSxxx6jf//+9O/fn8cee4zExES+/vWvA5CWlsadd97JrFmzyMzMJCMjg3vvvZdBgwad0bxMRERE5IuUSIuISItFsrS7rf3oRz/C7Xbzve99j/LyckaMGMGyZctISUlp2ubpp5/GZrMxY8YM3G43RUVFzJ8/H6vV2ubxiIiItCetHkeWEmkREekQVq9e3exrwzCYPXs2s2fPPueYhIQE5syZw5w5cyIbnIiIiHQoSqRFRCQk8bIiLSIi0pmp2VhkKZEWEZEWU6MRERGR+KD5OrLUtVtEREREREQkBFqRFhGRkKi0W0REJPZpbo0srUiLiIiIiIiIhEAr0iIi0mLxdPsrERGRzkzNxiJLibSIiIREibSIiEjs09waWSrtFhEREREREQmBEmlptbfffpt/+7d/48ILLyQpKYlu3bpxww03UFxcfMa2W7duZeLEiSQnJ9OlSxemT5/Op59+2myb3bt3c++99zJ06FC6dOlCRkYGY8aM4dVXXz1jf4cPH2bmzJmMHTuWLl26YBgG8+fPD/k9tCQugGeeeYbp06fTu3dvDMNg3LhxIR+rpqaGmTNn0rVrVxISErj00ktZuHDhGdutX7+e73znOwwdOhSn04lhGOzfvz/k44m0pdOl3ZF4iEhkab4OTUvm60AgwFNPPcXUqVPp3r07iYmJXHTRRdx///1UVFSEfEyRtqT5OrKUSEurzZ07l/379/ODH/yAt956i1//+teUlpYycuRI3n777abtPv74Y8aNG4fX6+Uvf/kLf/zjH9m9ezdXXnklJ06caNpu2bJlvPnmm9x888389a9/5aWXXqJ///7ccsst/PznP2927D179vDSSy/hcDi45pprwoq/pXEB/M///A8HDhxgwoQJZGdnh3W86dOn88ILL/Dwww+zePFihg8fzr/+67/y8ssvN9tu5cqVrFixgp49ezJ69OiwjiUiInKa5uvQtGS+drvdzJ49m4KCAp555hneeustvvvd7/Lcc88xZswY3G53WMcW6agef/xxDMNg5syZTc+Zpsns2bPp2rUrLpeLcePGsXPnzmbjPB4P99xzD1lZWSQlJXH99ddz+PDhdo6+OcPUaQVppdLSUnJycpo9V1NTQ79+/SgsLGTFihUAzJgxg1WrVrF3715SU1MBOHDgAP379+c///M/+eUvfwnAyZMnyczMxDCMZvucNm0aq1atoqysDKfTCUAwGMRiaTgftGXLFoYPH87zzz/PHXfc0eL4WxrXl49XWFhIVlYWq1evbvGx3nrrLa699lpefvll/vVf/7Xp+cmTJ7Nz504OHjyI1Wo941i/+tWv+K//+i/27dtHr169Wnw8kbZ02223kZ2dzbe//e023/eBAwe49dZb9aFTJII0X69u8bFaOl8HAgEqKirIzMxsNv7VV1/llltu4cUXX+Qb3/hGi48r0hY++eQTBg0axKZNmyKy/x//+MeMHTuW++67L6RxmzdvZsaMGaSmpjJ+/HieeeYZAH75y1/y6KOPMn/+fC644AIeeeQR1q5dy65du0hJSQHgP/7jP3jjjTeYP38+mZmZzJo1i7KyMoqLi5s+O7c3rUhLq315UgZITk5m4MCBHDp0CAC/388///lPbr755qbJD6CgoIDx48ezaNGipueysrLOmJQBLr/8curq6igrK2t67vQkGa5Q4mqL4y1atIjk5GRuueWWZs9/+9vf5rPPPmv2B6+1xxKJBJV2i8Qvzdct19L52mq1npFEQ8P3AGj6vopEQyzN1zU1Ndx222384Q9/ID09vVmMzzzzDA8++CDTp0+nsLCQF154gbq6uqbqj8rKSubNm8eTTz7JxIkTGTJkCAsWLKCkpKTpBGA06JO6RERlZSVbt27l4osvBmDv3r243W4GDx58xraDBw9mz5491NfXf+U+V61aRXZ29lk/CISrLeIKxY4dO7jooouw2Zo3zD99/B07drTZsURERM5H8/XZtXa+Pl0qf/r7KtLZ3X333Vx77bVMnDix2fP79u3j2LFjTJ48uek5p9PJ2LFj2bBhAwDFxcX4fL5m23Tt2pXCwsKmbaJBt7+SiLj77rupra3lwQcfBODUqVMAZGRknLFtRkYGpmlSXl5Ofn7+Wff3v//7v6xevZpf//rXbVq+0dq4wjlenz59znqsL8YjEsu0eizScWi+Pvfxwp2vjxw5wv3338+wYcOYNm1am8QjEo5Iztd+v5+qqqpmzzmdzqbLOb5o4cKFbN26lc2bN5/x2rFjxwDIzc1t9nxubi4HDhxo2sbhcDRbyT69zenx0aAVaWlzP/3pT3nppZd4+umnGTp0aLPXzlYCdr7XFi9ezN13383XvvY17rnnnrBiCgaD+P3+pkcgEGh1XOdimmazY/n9/ogdS0REJFyar9t+vi4rK+Oaa67BNE1eeeUVXaYlHdaqVatIS0tr9nj88cfP2O7QoUP84Ac/YMGCBSQkJJxzf1/+nTJN87y/0y3ZJpL02y1t6mc/+xmPPPIIjz76KN///vebnj99/dDZzuCWlZVhGAZdunQ547WlS5cyffp0Jk2axEsvvRT2L8u//du/Ybfbmx5FRUWtiuurrFmzptmx7HZ7022rMjMzz3ksOPuZdpFYomukRToGzddtP1+Xl5czadIkjhw5wvLly8+6oi3SniI5Z48fP57Kyspmjx//+MdnxFBcXExpaSlDhw7FZrNhs9lYs2YNv/nNb7DZbE0r0V9eWS4tLW16LS8vD6/XS3l5+Tm3iQaVdkub+dnPfsbs2bOZPXs2DzzwQLPX+vbti8vloqSk5IxxJSUl9OvX74yzVEuXLuXGG29k7NixvPbaazgcjrBjmz17drMPCqc7AIYT1/kMHTr0jNKVrl27AjBo0CD+/Oc/4/f7m113dfr4hYWFIR1LpL1FKulVIi3SfjRfN2jL+bq8vJyJEyeyb98+Vq5cedZruUXaWyTnVpvN1qzx37kUFRWd8Xv77W9/mwsvvJD77ruPPn36kJeXx/LlyxkyZAgAXq+XNWvWNHXiHzp0KHa7neXLlzNjxgwAjh49yo4dO3jiiSfa+J21nFakpU3893//N7Nnz+YnP/kJDz/88Bmv22w2rrvuOl5//XWqq6ubnj948CCrVq1i+vTpzbZftmwZN954I1dccQV/+9vfznq9RSh69erFsGHDmh4DBgwIK66WSElJaXasYcOGNX2ouOmmm6ipqeG1115rNuaFF16ga9eujBgxohXvUkRE5Ktpvv5cW83Xp5PoTz/9lGXLljUlAyLS8HtWWFjY7JGUlERmZiaFhYVN95R+7LHHWLRoETt27OCOO+4gMTGRr3/96wCkpaVx5513MmvWLFauXMm2bdv4xje+waBBg85oXtaetCItrfbkk0/y0EMPMXXqVK699lo2btzY7PWRI0cCDWfAhw8fzrRp07j//vupr6/noYceIisri1mzZjVtv379em688Uby8vJ44IEH2L59e7P9DRw4sNkZsFdffRWATz/9FGi4P2VycjIAX/va184bf0vjOr3v02VfVVVVmKbZdPzhw4dTUFDwlce6+uqrmTRpEv/xH/9BVVUV/fr1489//jNLlixhwYIFzRqznDhxgjVr1gCfnwFfvHgx2dnZZGdnM3bs2PO+N5FI0Iq0SHzSfN3287Xb7WbKlCls27aNZ555Br/f3+z7mp2dTd++fc/73kQiIVJza1vv90c/+hFut5vvfe97lJeXM2LECJYtW9ZUkQLw9NNPY7PZmDFjBm63m6KiIubPnx+1e0gDGKY+vUgrjRs3rinhO5sv/ogVFxdz33338e6772Kz2ZgwYQK/+tWvmk0ys2fP5mc/+9k597dq1SrGjRvX9PVXXYfV0h/vlsQFcMcdd/DCCy+cdR/PP/88d9xxx3mPVVNTw4MPPshf/vIXysrKuPDCC/nxj3/Mv/zLvzTbbvXq1YwfP/6s+xg7diyrV69u0XsTaSu33XYbGRkZ3H777W2+74MHD3LbbbfhdrvbfN8i0kDzdYO2nK/3799P7969z7mP22+/nfnz57fovYm0lU8++YRBgwaxfv36iOz/Jz/5CePHj+e+++6LyP7jhRJpERFpkdtuu4309PSIJdLf+MY3lEiLiIi00ulEet26dRHZ/09/+lMl0qi0W0REQqTSbhERkdinuTWy1GxMREREREREJARakRYRkRbT7a9ERETiQ7w0G4tXWpEWERERERERCYFWpEVEJCQ6Ey0iIhLbIlVBJp9TIi0iIi2m0m4REZH4oLk1slTaLSIiIiIiIhICrUiLSFgMwwh5TEpKSgQiOVN9fX3IY/o6HBT5/RQFAozz+0kHNlssFCUlnXNMVVVVK6KMX1qRFhGJL6mpqe1ynHDm+bKyspC2t5smE10uxnu9jPN6GeT3A/AvaWmscjrPOe748eMhxxbv1GwsspRIi0inlGCaXGWaTDJNJgeDXOTzEQCKLRb+x+HAAtzn9dItGOSIRcU7IiIi0dInGGRS42NsMEiyx8MJw2C1w8HcxET+q7aW6zyer0ykRdqaEmkR6RxMk4HQMBE3JtEJwCFguWHwuNPJGpuN8saV9lTT5D+9Xq73+5nrcEQz8piia6RFRCTSkk2TccEgExuT5z6miQ941zD4pc3GhpQUdtpsmI1z9kV+P193u/kv0yQQRsVcR6W5NbKUSItIh5VumhSZJpOCQSaaJj0AN7DOMPipxcIyi4WPAAyDBLu92dgqw2CV1apEWkREJMIM0+QS02RiMMjkQICRpokd2GsYLLNYWGGxsMZioeb0ye4vzdn/dDq5p66OUT4f6zVnSztRIi0iHYbVNBkaDDIuEGCSaTLMNLECHwKvWSwsNwzWGQb1LTxb/Xe7nWfr68kJBilVeXcTrUiLiEhrZQUCjPV4GOX1UhQMkgPUAKstFu61WllusbCvhXPvdpuNQxYL13o8SqS/QHNrZCmRFpG41i0YPKNJWDnwtmEwz2plhWFwOMwyr7dsNoLANL+fP2piBlTaLSIi4bGbJsO8XsZ6PIyvr6ewsUnYNsPghcbEeaPFgi+cOdsweMvp5EaPhweSk5tKvjs7NRuLLCXSIhJXEkyTMYEAExuT5wuDwWZNwlbabLzj9bbJNVJlhsE6q5UblEiLiIiErJffz1iPh3H19Yz2ekk2TU5aLKxxOvl9cjJrnU52V1a2ybH+6XRyl9vNcJ+P9zRnSztQIi0isc00uTAYpKgxeR4dCOACjhgGK2w2HnM4mjUJA9q00cjfbTZ+5fGQYZqU6Qw3oDPRIiJydknBIKO9XsbX1zPW46FXIIAP2OxwMCc5mdVOJzvt9oisGG+22zlmsTDN41Ei3UjzdWQpkRaRmNPFNBnn9zMxEGCC309306QeeMdq5edOJyutVj62WKAdEtt/2mw85fFwjd/Pgi81NxEREenMDNNk8BcusRpZXY0d2G+1strpZHVCAhscDmrboc+I2VjefY3Hw0PJye3yGUE6NyXSIhJ1p5uEFfn9FPn9DA0GsQIfWSz8zWZrKNe2WlvcJKwtlVosvGu1coPPp0QaXSMtItLZZTdWiRX5/YwPBMgxTWqAtVYrD6elscbpZL8tOinGm04n/+Z2c6nfz3bN2ZpbI0yJtIhExbmahK222fiT3c5Km40jMdIp+x82G//t8ZBqmlTpDLcSaRGRTsRumowIBJqS50uDQQC2Wyy8aLez0mplk9WKzzBISUqKaqzv2u2cMgymeTxKpFGzsUhTIi0i7SIBmlaciwIBLjpLk7Bii6VNr29uK/+w2filx8PVfj+vaGIWEZEOrncw2NDU0+/nykCAFOCEYfC21crvHA5WWq2ciJGT3V8UMAyWOJ1M83h4JClJ5d0SUUqkRWLYRRddFPIYn88X8pju3buHPGbo0KFfvYFp0qe+nlFVVYysqmJIdTUJbjfHbDbWp6Tw+8RE3k1KotJqbRqSG3IUZ2cPI9mtqak552v1wHafj69ZLKzIzGxFZPFPpd0iImd3xRVXhDwmEAiEPGbHjh1tfpxk02RsMMjEQICJwSB9TRMf8K7Fwv+z2VhhtfK+YXzeJCwYbHh8SW1tbcixpaWlhTzmq+aMfzqd3FZfz0C/nw+/8HkgnM86hw8fDnlMLNHcGllKpEWkzaT4/VxeXc2oykpGVlWR5/PhMQy2pqQwt1s3Vtps7HE44vIM8eKkJP6zooLEYJC6GDwLLyIi0lKGaXKJaTYlziODQRzAp4bBCouFH1utrLFYqI7D+Xqdw0FlY3n3h6oikwhSIi0iYbOYJhfX1jatOl9cW4sV+DQhgbfT03k3NZWtKSl4GhPPr1r1jXWLExP5cXk549xu3oryNWDRpBVpEZH4lG2aDbeSbGwWlgvUAGssFu6z21lusfBpBzhR7DMMljqdXOvx8ERycrTDiSrNrZGlRFpEQpLt8TC8rIwhJ05weVUVaYEAVVYr76Wk8PeCAjampnK8A96/8ZDdzk6Hg6m1tZ06kQY1GxMRiQenm4RNPN0krPHv7PuGwYs2GyssFt61WPDF4arz+bzpdDKjvp5+fj97otRBPBao2Vhkxf9pJxGJKEcgwPCyMr63Zw/Pb97MXzduZNbu3eR5vbySk8O3Bwxg0iWXcH/fvvw9K6tDJtGnLUlMZLzbjfMs14VJ+5o7dy6DBw8mNTWV1NRURo0axeLFi5tev+OOOzAMo9lj5MiRzfbh8Xi45557yMrKIikpieuvvz7ur4cTkc6tdzDId3w+FrrdHKitZXF9Pd/y+/nIYuFOu53eCQmMSkjgIbudtY2dtjui1Q4HtY3l3SKR0nlP0YjI2ZkmBXV1DC8v5/KyMi6prMQZDHLC4eC9jAxeKCigOD2do51wclqcmMisigquqq9neWJitMOJilgp7e7evTu/+MUv6NevHwAvvPACN9xwA9u2bePiiy8GYOrUqTz//PNNYxxfOskzc+ZM3njjDRYuXEhmZiazZs1i2rRpFBcXY/1CEzwRkViVbJpc9YVbU51uErbJYuEph6OhSZjFgr+TnQCuNwxWOBxM83h4phNXkWnlOLKUSIsIyT4fQysquLysjOHl5eR4PHgNg/e7dGFer15szshgX2Ji8yZhnTCR3utwsNtuZ2ptbadNpGPFdddd1+zrRx99lLlz57Jx48amRNrpdJKXl3fW8ZWVlcybN48XX3yRiRMnArBgwQJ69OjBihUrmDJlSmTfgIhIGAzTZPDpcu1AgJGBQFOTsJVWKw/abKy1WuOySVhbe8Pp5H+rqigIBDigk6MSAUqkRTohi2ky2OPhyro6rqqrY/DevViB/YmJrMnKYnNGBu+npeHRxHOGJYmJ3FFVhb0Tn+WN5BnuqqqqZl87nU6cTudXjgkEAvz1r3+ltraWUaNGNT2/evVqcnJy6NKlC2PHjuXRRx8lJycHgOLiYnw+H5MnT27avmvXrhQWFrJhwwYl0iISM9J9Pi6vqGBEZSUjKivJ8PmoAdZZrdzvcLDSZusQTcLa2kqHAzdwbX09v+ukq9K6RjqylEiLdBJ5fj9X1tVxZV0do91uugSDVFksbHC5eKqggM3p6ZQmJEQ7zJi3JCmJ/1tZySi3O9qhREUkS7v9fv8Z9xN9+OGHmT179lnHlJSUMGrUKOrr60lOTmbRokUMHDgQgKuvvppbbrmFgoIC9u3bx09/+lMmTJhAcXExTqeTY8eO4XA4SE9Pb7bP3Nxcjh071ubvT0SkpWzBIINqahjRmDxf2Hhv5t2JifwzO5uFZWVsslrxatX5K9VZLKxqLO9WIi2RoERapINyBoMMr69vWnXu7/MRBD5wOnkxLY21iYl84HQSMAy65+dHO9y48ZHdzn6bjavr6qIdSodjs9k4depUs+e+ajV6wIABbN++nYqKCl577TVuv/121qxZw8CBA7n11lubtissLGTYsGEUFBTw5ptvMn369HPu0zRNDH04FZF21q2+vilxHlpZSVIwSLnNxntpafwlL49NaWmUNfZ52PGlyh05t38mJPC7qiq6BgKc0D2lpY0pkRbpKEyTvl4vo6qquLKujsvr60kwTY5ZraxLTGRORgbvuFxUqly7dQyDJYmJzKipAb8fOuFtNSLZbCw1NbXFYxwOR1OzsWHDhrF582Z+/etf8/vf//6MbfPz8ykoKOCTTz4BIC8vD6/XS3l5ebNV6dLSUkaPHt2atyIicl6uQIDLKiu5vKyMERUVdPd48BsGJcnJvNitGxvT0tidlISpE3utstzhwAtc6/EwvxNW3WlFOrI63ydAkQ4kNRBgZG0tV9TWMqa2lq5+Px7DYHNCAk9nZLDW5eITh6N5kzBptcVJSfyfqipYuxYmTIh2ONLINE0852iCd+rUKQ4dOkR+Y/XF0KFDsdvtLF++nBkzZgBw9OhRduzYwRNPPNFuMYtI52CYJv1qaxlRUcHl5eUMqq7GbpoccTrZlJbGb7p0oTg1lbpOeHI2kqosFtY4HA2JdLSDkQ5Hv60iccRimgyqr2dMY/I82O3GSkM36RUpKaxLSuJdu516NR2JqA8cDo5YrXR77bVOl0jHyu2vHnjgAa6++mp69OhBdXU1CxcuZPXq1SxZsoSamhpmz57NzTffTH5+Pvv37+eBBx4gKyuLm266CYC0tDTuvPNOZs2aRWZmJhkZGdx7770MGjSoqYu3iEhrdPF6GV5RwYiKCoZXVJDp81FnsbAtLY05vXvzXpcuHLDbdbI7wt50OnmquprsQIATnawqT83GIkuJtHQYLpcr5DGhlJGeVlpaGvIYgOTk5JDHJCUlkR8IMM7rZbzHw5VeL+mmSaVhsM7hYEFKCqudTo6cnhh8Pmxh/HGrbWxkEopwriO97LLLQh4TjoMHD4Y8xuv1hrT98pQU7li0CObMgU524iIWEunjx4/zzW9+k6NHj5KWlsbgwYNZsmQJkyZNwu12U1JSwp/+9CcqKirIz89n/PjxvPLKK6SkpDTt4+mnn8ZmszFjxgzcbjdFRUXMnz9f95AWaQeXXHJJyGO+3NW/JXJzc0MeEy5PdTWXut2Mrq5mTE0NA+vrAfg4IYG/p6WxPjmZ7YmJ+E7PGV4vvjDm34yMjJDHBAKBkMeEw+fzhTzmXJVEXyWUz2+r7HaC1dUUVVd3uvJuJbyRpURaJMY4TZMxjfeInOR2c2EgQBDYbrPxx8REVjkcbLPbCegMdlQtS0nhjkOH4N13YcyYaIfT6cybN++cr7lcLpYuXXrefSQkJDBnzhzmzJnTlqGJSCfStb6ekZWVjKyoYGhVFUnBIGVWK+8mJ7MgM5MNycmcUpOrqKqwWtnodDLN6+10ibRElhJpkWgzTQaYJhP9fooCAa4IBHABnxkGa5xOnk5OZq3DQXknW/WMdVsTEiAvD157rVMl0rFS2i0iEg2uQIChVVWMaEyee3yhSdi8rCzeSUnho4QENQmLMYsTE/lZeTnpwWCn+jyluTWylEiLREEX02RsINCUPPcwTeqBDVYr/+1wsMJq5SOLhaQwysGlfZiGATfd1JBIP/mkrnETEemADNOkX11d06rz4JqaZk3C5qSlsaWxSVg4l0lJ+1iamMh/l5dztdfLy1qVljaiRFqkHViAYY3l2hP9foYFg1iBXYbBP2w2VlqtrLdacSsZiy9f+xrMnQtbtsDw4dGOpt1oRVpEOrJ0n4/LKysZ0fjI9PlwWywUp6by65492dSlC4ecTp1AjSMnrFY22Wxc18kSaTUbiywl0iIR0hWY0viYCGS63VQAa6xWZjqdrLRaOdSJyos6pKuugszMhlXpTpJIR6q0W0QkWqzBIINraprKtS+sqwNgd2Iib2VlsTEtjQ9SUj5vEiZx6Q2Hg9l1daQEg1R3kv+Xmq8jq3P8FIm0AycwCfgVUAIcAf4X6AX8FpjkctErKYlvuFzMt9uVRHcENhvceGNDIq3JSkQkbnStr2f68eP8cvdulhUXM/ejj7ihtJQDLhc/69OHa4cM4VuDBvHbnj0pTktTEt0B/NPhwAFMCaOzuLTO3LlzGTx4MKmpqaSmpjJq1CgWL17c9Podd9yBYRjNHiNHjmy2D4/Hwz333ENWVhZJSUlcf/31HD58uL3fSjNakRZphQv5fNV5LJBIQwK9FPhvYAVQ1rhtsm6p0zHdfDPMmwclJTB4cLSjaRcq7RaRePNVTcJe7NqVjWlp7EpKUpOwDuwzq5Vim41pXi+vOp3RDqddxMrc2r17d37xi1/Qr18/AF544QVuuOEGtm3bxsUXXwzA1KlTef7555vGOByOZvuYOXMmb7zxBgsXLiQzM5NZs2Yxbdo0iouLo3bbSiXSIiFIDQa5ks+T556AB1gLPERDAr0jeuFJNBQVQVoavPpqp0mkRURinmlykdfLtZ99xojKSi6prsZumnzmdLLxS03CpPP4p8PBj+rqSDJNanXSpN1cd911zb5+9NFHmTt3Lhs3bmxKpJ1OJ3l5eWcdX1lZybx583jxxReZOHEiAAsWLKBHjx6sWLGCKVOmRPYNnIP+eoh8BYtpcqnfzzivl/FeL0N9PqzAR8DrNCTOawB3VKOUqHI44PrrG8q7f/7zaEcTcbr9lYjEqoxAgCvdbq5yu7nS7SY7GGxqEvabnj3ZmJbGoYQENQnrxN5wOHi4ro4ir5d/dIJV6Ug2G/P7/VRVVTV73ul04jzP9zUQCPDXv/6V2tpaRo0a1fT86tWrycnJoUuXLowdO5ZHH32UnJwcAIqLi/H5fEyePLlp+65du1JYWMiGDRuUSIvEirxAgPGNifOVXi8ZpkmlYbDO4eC+lBT+Wl3NwWgHKbHl5pvhxRfh44/hwgujHU3EKZEWkVhgM00u83iaEufBXi8AH9rtvJaczFqXi0M9e+r6Zmmy32rlA6uV65RIt9qqVav4yU9+0uy5hx9+mNmzZ591+5KSEkaNGkV9fT3JycksWrSIgQMHAnD11Vdzyy23UFBQwL59+/jpT3/KhAkTKC4uxul0cuzYMRwOB+np6c32mZuby7FjxyLy/lpCibR0ek7TZITPx3iPh/FeLxcFAgSB7TYbz7tcrHI42Gq3E2g8g11aXR3dgCX2TJ4MSUkNq9IPPhjtaEREOqwePh9XNa46j6qvJ8U0OWWxsN7l4oXUVNYlJHDiC+XauUqi5Uv+6XDwf91unKaJR9UJYRs/fjyvv/56s+e+ajV6wIABbN++nYqKCl577TVuv/121qxZw8CBA7n11lubtissLGTYsGEUFBTw5ptvMn369HPu0zRNjCj+P1QiLZ2PaTLANJkUDDK1ooJRXi+JwDGLhVUOB88kJbHW4aBMk6+0lMsF117bKRJplXaLSHtyBQIMr61ldHU1Iysr6e334we2Op38Pi2NtS4XOxwONQmTFnvD4eABt5vxPh9LvtTQqqOJ5Nxqs9lITU1t8fYOh6Op2diwYcPYvHkzv/71r/n9739/xrb5+fkUFBTwySefAJCXl4fX66W8vLzZqnRpaSmjR49u5TsJnxJpCUmvXr1CHhPOmaK6xns4tpU002R8MMjEYJBJgUBTk7CNdjtPJCez2uHgI6s1otdN1dTUhDzG5XKFPKa2tjbkMQcPhl6sbrfbQx7j9/tDHhOOgoKCkMds2bIl5DF33HFH07+HHz/O3du28V8338yJlJRzjpk/f37IxxERCUfv3r1DHpOZmRnymDZdETJNBtTXM7q6mtHV1Qypq8Numhyx23nLNFlht7PaYqHKMKC+vuHxFTweT1hhhNMF+NSpUyGPCQaDIY+xhHGiP5z/r+EI52fhfNfUnk11GNWBpz9bbgM+Mgym1tbyauPlANL+TNM85+/nqVOnOHToEPn5+QAMHToUu93O8uXLmTFjBgBHjx5lx44dPPHEE+0W85cpkZYOyWKaDDXNpsR5uGliA3YZBm9YrSy3WFhnsWBLS4t2qNJBfNCtG16rlWEHDrC4sDDa4USMVqRFpK2l+/2MrK5mdE0No6uryfL7cRsGm5OTeTI/nw0pKRxwONjduDol0lp/s1j4j0AAu2ni68DVDJFsNhaKBx54gKuvvpoePXpQXV3NwoULWb16NUuWLKGmpobZs2dz8803k5+fz/79+3nggQfIysripptuAiAtLY0777yTWbNmkZmZSUZGBvfeey+DBg1q6uIdDUqkpcPoapoUBQJMDgYZHwySCVQCqywWfmC1ssJq5eCX/li2vCBF5Kt57HZKunbt8Ik0qNmYiLSONRjkspqapsT5IrcbC7ArIYE30tPZkJzMtqQkNQmTiHndYuHHgQDjTJPlSqQj7vjx43zzm9/k6NGjpKWlMXjwYJYsWcKkSZNwu92UlJTwpz/9iYqKCvLz8xk/fjyvvPIKKV+o8Hv66aex2WzMmDEDt9tNUVER8+fPj9o9pEGJtMSxczUJ22oYPGe1stxq5T3DaGoSJhJpWwoKuGv9etJraylPSop2OCIiMaOrx8PIykpGVlUxvKqK5GCQcquVd1NS+HNmJu+mpHAyjEuGRMJRYhjsBW4KBFiuEzYRN2/evHO+5nK5WLp06Xn3kZCQwJw5c5gzZ05bhtYqSqQlfpgm/b5wa6ovNwl73GLhbYuFU0qcJUre79EDv8XC0IMHWXHRRdEOJyIiVdotIh2LKxBgaHU1I6uqGFlZSYHHgx8oSU7mxbw81iYm8pHLpSZhEh2GwSKrldsDAe4xzQ676BIrpd0dlRJpiWmpwSBXNibO47xeegSDeIBNdjv/LymJVU5nU5OwL98UXqS91Tkc7MzPZ9iBAx02kRYROSvTpL/bzcjKSkZVVXFpTU1DkzCHg41paTybmsrmlBRqG29N5VWTJ4myv1ks3BsIcIVpsqaDJtISWUqkJaZYTJNBXi8ja2uZ4PVymc+HDfjEamWJ08kqh4N3HQ7q9AdPYlRxQQF3bNhAmttNZRhd1+OBrpEWEYAuPh8jqqsZVVXFyKoqsnw+3BYLxSkpPNO9OxvT0jjodEb0jhgi4dpiGByiobx7TQcs71YFWeQpkZaoy/H7ucrt5iq3myvq60kPBqkyDNY5HNyfksJqh4NDUWwkIBKKrT16cLthcNnBg6waMCDa4bQ5de0W6byspsmgmhpGNybOF9bVYQF2u1y8lZHBxrQ0ticnq0mYxIfG8u6vBQL8p2l2yMsMNLdGlhJpaXeOYJDhHg9j3W6udLu50OcjCJQ4HCxISWGty8WGQAB/B/yDJh1fTUICH+flMfTAgQ6ZSItI59LV42lacT7dJKzCamVTaiqv5OSwMTWVMocj2mGKhOVvFgv/NxBgpGnyrj53SoiUSEvkmSZ9fD6uqq/nKrebkfX1uEyT41Yra10ufpuWxjsuF+VfWHX219VFMWCR1tnSsyffeO89kurrqU1IiHY4bU4r0iIdV0IgwNCaGkY3dtj+cpOwd1NT+TgxkeAXkg6lHxKv3jUMjgI3BoO82wErKdRsLLKUSEtEpASDjK6vZ6zbzVi3m26BAB5gc0ICz3TpwlqXi4/tdl03JR3S1oICvrlpE0MOHWJ9//7RDkdE5NxMk4t8Pq5qrBIbfuAADtPkM4eDjampzOnWjc2pqdTqEivpgEzD4B8WCzcFAtzX2LxWpKWUSEubON0k7Cq3m6vq6xni8WAD9tpsLEtMZK3LxaaEBNwd8GyfyJdVulx8kpPD8AMHOlwirWukReJfRiDAFY29Sa6srycnEKDOMNiYkMCvu3dnY2oqB9QkTDqJv1mt3BUMMtQ0Ke5gP/OaWyNLibSELcfvZ2xjufYXm4RtSEjgoYwM1rpcHLbZMDrYHyWRliguKGBGcTEurxd3B7t+UIm0SHyxmSZDPJ6mVedBXi8W4GO7nUVJSax1udiSkIDXMMjMzIx2uCLtaq1hcJKG8u7iDrbgo7k1spRIR1hiYmK0QzinhBCv3XSaJkPr65kYCDDR76fQNAkC2ywWnrNaWeFwsNliaWgS5vE0PIDa2tqQY8vNzQ15jNPpDHmMq4Penkiib0vPnnx982YuPXyYd/v0iXY4ItICF198cbscJz09PeQxPp8vpO3zPR5uPnmSCV4vV/h8pJompwyDNQ4H/5uSwmq7neOny7WDQWjsTVJZWRlybD169Ah5TFVVVchjRFpj3Lhx53xty44dfL28nJVXXNGsEmP8+PEhH2fVqlXhhCdxSIm0nJtp0j8YZILPxwSfjzF+P4nAMcNghdXKk1Yrb1utnNKKs8gZypKT2ZuVxdADBzpUIq3SbpHYlBAIMLS6mhGVlYysrKSgvh4/sMVm47cuF6scDj6w2Zo1CRORButzc5l65Ai9a2rYl5IS7XDajJqNRZYSaWkmxTQZ25g4T/D76RkM4gHetdn4pcvFEmCnYei6KZEW2FJQwI3bt+Pw+fDa7dEOR0Q6EtOkv9vdlDhfUl3d1CRsU1oav+venVfLy6nuYKWqIpGwPTOTGpuNMcePd6hEWiJLiXQnZzFNLgkEKGpMnocFAtiATywWFtvtvG23847NRl1j4uz1eqMbsEgc2dKzJ7cWFzP4yBG29OoV7XDajM5Ei0RHms/HiKoqRlRWMqKykmyfj3qLheKUFOb06MHGtDQOJiQ0neyuDqNMW6Qz8lssbMrO5orjx1nQr1+0w2kzmq8jS4l0J5RvmhQFAkypqWGc30+maVINrLbb+S+nk7dtNg7pNhcirXYiNZWD6ekMP3CgwyTSKu0WaT/WYJCLq6sZUVHB8PJyLqytxQJ84nKxNDOTjWlpvJ+SglerziKttj43l6KjR+leU8Ph5ORoh9MmNLdGlhLpTsBhmowJBpkYDDIxEGBQY5Ow7VYrzzcmzltstoYmYSLSpjYXFHDNjh3YAwF8OkElIueRV1/PiIoKLi8vZ2hlJcmBABU2G++lpvJqTg6b0tI42cHuBCASC4qzsnBbrVxx/DgLO0giLZGlRLojMk36myYTg0EmBQJcGQySBBwDllut/Mpi4W2rlVp1rBaJuC0FBdy8fTsXHznC9p49ox1Om9CKtEjbSQgEGFJZyeUVFYwoL6dnY5OwD1NS+HO3brzXpQu7kpPx+P3RDlWkQ/NarbzXWN69sG/faIfTJtRsLLKUSHcQqabJ+MYV54nBIAWmiQfYYLHwmM3GCquVki81CQvt5lciEo6jXbpwJC2NYQcPdphEWkRawTTpW1vL5eXljKioYHBVFQ7T5JjTyaYuXfifXr3YmpZGjU0f0UTa2/rcXB58/33y6uo4FsO3sJXYoL/SccowTS4zzabEeUQwiA3YbRi8abGwwmplrcXS1CRMRKKnuKCAoo8/xhoIRDuUVtM10iKh6+L3M7q2ljE1NYyurSXH76feYmFbWhpze/ViU5cuHHS5dEcMkSjbnJWFx2JhzPHjvNa7d7TDaTXNrZGlRDqO5H0hcZ4QCJAFVAGrLRZ+aLez3GLhgBqOiMScLQUFXP/BBww8dizaobQJJdIiX81qmlzidjOmpoYxNTVcXF+PBdjldPLPtDQ+yMvjg9RUNQkTiTH1NhvFWVlcoURaWkCJdAxzmCajv9AkbHBjk7BthsE8m43lFgvvWSxqEiYS4w6mp1OanMywAweiHYqIREhXr5cxjavOI2trSQkGqbBa2ZCUxMKMDDYkJVHaeD/59C5dohusiJzTO7m5/FdJCVn19RzRZ2z5CkqkY8lXNAlbYbXyVGOTsJP6pRaJL4bB5l69uPKTT8Dvhzi+9lGl3SINXMEgw2prGVNbyxU1NfT2evEDH7hczM/M5J3kZHYmJBDUnC0SVzZmZ+MzDMYcP85f8vKiHU6rqNlYZMXvp7kOItU0Gde44jyphU3CRCT+bOnZk2t37IB162D8+GiHIyKhMk0u8HgayrVraxlaV4fDNPnMbuedpCSeyclhU1IS1brNnUhcq7Xb2ZaZ2SESaYmsmEukhw0bFvKYHTt2hDzGGsZEl5SUFPKYYDDY7OvTTcImBYNMCgYZYZrYgE8MgzcNg+WNTcJqv5g4m2bD4yukp6e3OraWCOcMVG5ubshjCgsLQx7Ts506Ii9ZsiSscSdOnAh5jM/nC3mMI4z7i+aFMVEcC+N634SE9ukVn52dHfKY733veyGPCeX/6Z70dE4mJpL12mtxnUhrRVq+yOl0hjzG6/VGIJLIyAAmAVOByUBXoA5YDfwXsBTY5fNBRUXDo4VycnJCjiWcOdvj8YQ8xhLGdduHDx8OeUxtbW3IY8L57BYIs8ljOMcK5/POkSNHQh4Tzlya3E73RbaFUXEVzphu3bqFPMYfwi3k3snN5Qc7d5Lu9VIex/dt19waWTGXSHdEeY3l2pODQSYEg01NwlZZLMy0WllhsbBfK84iHZthsKVnT6a+/jr85jcQx02GNDFLR2UFRgJTGh/DAAtQArxMQ+K8Dgg9PRWRePJuTg7/98MPueLkSd7o2jXa4YRN83VkKZGOAIdpMsLnY7zXy1VeL4Mbf4iLDYN5VivLLBY2GYaahIl0Mpt79mTqsmXw7rswZky0wxERoCcNSfNUoAhIA04By4G5wDLgs6hFJyLRUOVw8EF6OledOBHXibRElhLptmCa9AkEmOD1Mt7nY7TXSxJQahgsMwyetFpZabGoSZhIJ/dJVhbk5cFrr8VtIq3Sbol3icBYPl91vhAIABuBJ4ElQDEQeiG1iHQk63Nz+d5HH5Hq81HV2HE/3qjZWGQpkQ5TSjDIlY2rzuO8XgqCQbzAJrudp5KSeNvhYKfVSlV1dbRDFZEYYVoscNNN8Prr8OSTaiIo0k4KaVhxngJcCTiBAzSUaj8IrAQqoxadiMSiDbm53P3RR4w+eZIl+fnRDkdikBLpFjJMk0t8voZVZ6+XYX4/NmCv1cpyh4NVDgcbHI7mTcJERL7s5pth7lwoLoYwmivGAq1IS6w73SRsCg1NwrrR0CRsDfAjGpuERS06EYkH5U4nO9LSuOrEibhNpDW3Rlb8drtpB3mmyW1+P897POxzu1leUcH33W7KLRZ+nJzMsIwMRmVk8EBKCsudTiXRInJ+Y8dCZmZDeXccOl3aHYlHKObOncvgwYNJTU0lNTWVUaNGsXjx4mZxzp49m65du+JyuRg3bhw7d+5stg+Px8M999xDVlYWSUlJXH/99WF1H5boswJjgJ8Dm4ATwEIamoX9mYakOgO4BvgNSqJFpGXWZmUxtLycpBA6fseSWJmzOyol0l/gME3GBQI84vWy0e1mr9vNc14v/UyT5202rk9LY0BmJrenpfGCy8VB3StSREJls8ENN8Crr573tnZybt27d+cXv/gFW7ZsYcuWLUyYMIEbbrihKVl+4okneOqpp3j22WfZvHkzeXl5TJo0ieovXG4zc+ZMFi1axMKFC1m/fj01NTVMmzYt7NvlSPvqCXwXeBU4CawH7gY+Be6kYRV6MA23qVqBOm2LSOjWZWfjME1GnToV7VAkBnXu0m7TpJ9pMjEQYGIgwFXBIEnAcWCF1cpTVitvW61NTcKS4vg+ciISQ26+Gf74RygpgcGDox1NyGKhtPu6665r9vWjjz7K3Llz2bhxIwMHDuSZZ57hwQcfZPr06QC88MIL5Obm8vLLL3PXXXdRWVnJvHnzePHFF5k4cSIACxYsoEePHqxYsYIpU6a0zRuTNuOioUnY6Wudv9gk7CkayrW3oCZhItJ2ShMS+CglhStPnGBFGPcJjzatHEdWp0ukU0yT8cEgkxuT516miRd412LhF3Y7K6xWSgwDU2XaIhIpRUWQltZQ3h2HiXQkVVVVNfva6XTidDq/ckwgEOCvf/0rtbW1jBo1in379nHs2DEmT57cbD9jx45lw4YN3HXXXRQXF+Pz+Zpt07VrVwoLC9mwYYMS6RhRyOfdta8EEoCDNHTWVpMwEWkPa7OzuX3/fhICAepVjSpf0OETacM0GWKaTAoGmRQMMtI0sQGfGAaLrVZWWK2ss1h0fbOItB+nE667riGR/tnPoh1NSCJ5+yu/309aWlqz5x9++GFmz5591jElJSWMGjWK+vp6kpOTWbRoEQMHDmTDhg0A5H5p9SA3N5cDBw4AcOzYMRwOB+np6Wdsc+zYsTZ6VxKqr2oSdj8Nq84fRy06EemM1mZnc9enn3L5qVOszcmJdjgh0Yp0ZHXIRDrPNJnYmDhPCAbJBqqA1RYL/2m1ssJi4aCtQ751EYkXN98MCxbAxx/DhRdGO5qQRCqRttlsnPrSdWhftRo9YMAAtm/fTkVFBa+99hq33347a9asaXrd+NIJUtM0z3jubHGcbxtpO1ZgBJ+vOg+noXlLCQ1NwpYC69D1zSISPZ+5XOxJSuKqkyeVSEd4v/GmQ2STDtNk9BeS50sa/+cWGwZ/tFpZbrGw0TDwf+HDkQozRCSqpkyBpKSGVekHH4x2NDEjNTW1xds6HA769esHwLBhw9i8eTO//vWvue+++4CGVef8L9yypLS0tGmVOi8vD6/XS3l5ebNV6dLSUkaPHt0Wb0XOoSefJ85FQBegDFgO/B5YBhyJVnAiImexNjubfzl0CHsggE/l3dKoxYl0e52hHzp06Pk3Mk16ejyMqqpiZFUVl3m9JAHHgJUWC09bray0WDjxFTE/8sgjIcf23HPPhTwmv53uO9e3b9+Qx5xuwhOKb3zjG+1ynF69eoU8ZuDAgSGPCYfb7Q5r3Pr160MeE0734KSkpJDHZGRkhDymW7duIY+prGyfqxlrampCHnM6IQvFl0uHW2LChAlN/34oMZFuv/gFd61c+ZVj3n777ZCPEymxfNsL0zTxeDz07t2bvLw8li9fzpAhQwDwer2sWbOGX/7yl0DDXGO321m+fDkzZswA4OjRo+zYsYMnnngiau+hrYQzZ5/vWvRwuUyTqxovsZocDDKA0JuEpaSkhHxcWxiVZ3V1dSGPCYfP5wt5TDi/d4mJiSGPcblcIY8JR3Z2druMufjii0MeA7Bt27aQxwSDobe6S0hICHlMOD+n/na6fVM483w4v9/h6NGjR8hjTv8dOZSTQ+L+/dyUlMQHYXxGjZZYna87irhZkU4KBBhWXc2oykpGVlXR3evFZxhsT07mscZVZzUJE5F4siYnh4d37iTf7eZoO3147SgeeOABrr76anr06EF1dTULFy5k9erVLFmyBMMwmDlzJo899hj9+/enf//+PPbYYyQmJvL1r38dgLS0NO68805mzZpFZmYmGRkZ3HvvvQwaNKipi7eEyTS52DSZ3Jg8jzHNpiZhyy0WHgwGWQlURDdKEZEWO5aRwdH0dIZ8+mlcJdISWTGbSBumyYV1dU2rzoNrarABB51ONqSl8W5qKltSUnBbrezYsSPa4YqIhGxTRgYei4UrT5zgLz17RjucFouF218dP36cb37zmxw9epS0tDQGDx7MkiVLmDRpEgA/+tGPcLvdfO9736O8vJwRI0awbNmyZisfTz/9NDabjRkzZuB2uykqKmL+/PlYVbYXsgzTZEIwyKTG5Lkb4AbWGgY/aTzZ/TGAYeD1eqMbrIhIGLb17s24HTuwBgIE4mSe0Ip0ZMVUIp0LXHvqFCMbk+d0v59ai4XNKSn8v5492ZiaypEIlZ6JiLS3epuNzRkZcZVIR7JrdyjmzZv3la8bhsHs2bPP2fEbGkoq58yZw5w5c0I6toDVNLncNJkYCDDJNBlmmliAnYbBXy0WllssrDcM6lUlJiIdxLY+fbhm61Yu+OwzPgqjTDwaYqXZ2Ny5c5k7dy779+8HGi65eOihh7j66qub9vezn/2M5557runk929/+9tml2Z4PB7uvfde/vznPzed/P7d735H9+7d2+x9hSqqibQDGMPnTUcuBdi/n48SE1mUlcW7qal8kJxMQBOxiHRQ67Kz+fFHH5FVX8/JMK6VE2kvPRpXmyc23hGjCw1NwlZaLPyvxcIKi4Ujmq9FpIM6nJnJidRUhnz6adwk0rGie/fu/OIXv2jqR/PCCy9www03sG3bNi6++GKeeOIJnnrqKebPn88FF1zAI488wqRJk9i1a1dTJdnMmTN54403WLhwIZmZmcyaNYtp06ZRXFwctUqydk+k+/F54jweSKahSdgy4P8BpYMHU263t3dYIiJRsSEzE59hcOXJkyyK4lnVUMTCirREnss0ufIL5doXmiYB4D3DYE5jufZmIKjkWUQ6A8NgW+/ejNq1i4VXXknQYol2ROcVK3Prdddd1+zrRx99lLlz57Jx40YGDhzIM888w4MPPtjUoPiFF14gNzeXl19+mbvuuovKykrmzZvHiy++2NTHZMGCBfTo0YMVK1YwZcqUdn9P0HC7xohKAW4AfgfsBT6hoUtnMvAIMAToCtwOvAxKokWkU6m129mans5VpaXRDkU6O9NkYDDID/x+3vR6Oer18g+/nxsDATYYBv9qs9HN4WC8w8GjNhvvWSxKokWkU9napw8p9fX0O3o02qFEnd/vp6qqqtnD4/Gcd1wgEGDhwoXU1tYyatQo9u3bx7Fjx5g8eXLTNk6nk7Fjx7JhwwYAiouL8fl8zbbp2rUrhYWFTdtEQ5sn0gZwGfBjYA1wCvgbMBF4C7gOyAQmAL8EtgOxca5ERCQ61mZnM6iykvQ4aMJ0+hrpSDyk/WUAM4B5wB6vl60+H7MDAQLAT61WLrXb6e9w8D27nUVWKxVKnEWkEzuQk0NZcjJDPv002qG0SCTn61WrVpGWltbs8fjjj58zlpKSEpKTk3E6nfyf//N/WLRoEQMHDuTYsWPAmbcUzc3NbXrt2LFjOBwO0tPTz7lNNLRJaXcuMJmGcu1JQA5QDbwN/ICGe0TGx4+biEj7eycrix/u3s2YEyf4Zxj3525PsdJsTMJjBS4HptIwZw+n4Yz6DuBVq5UVahImInJuhsH23r25bO9e/nLFFTF/291INhsbP348r7/+erPnnV/RFHrAgAFs376diooKXnvtNW6//XbWrFnT9Lrxpe+laZpnPHe2OM63TSS1OpGeT0NZNkAx8L80JM7vAr7W7lxEpBOocjh4Py2Nq+IgkZb41Q043PjvMmAF8BwNc/YRwGmLqRt5iIjEpG19+jChpITex4/zaV5etMOJGpvNRmpqaou3dzgcTc3Ghg0bxubNm/n1r3/NfffdBzSsOufn5zdtX1pa2rRKnZeXh9frpby8vNmqdGlpKaNHj26LtxOWVpV2p9OQRC+jYVV6GPAgsBYl0SIioVibk8OQigpSfbH/11Ol3fHph43/HQlkA7cCf6QhiRYRkZbZm5tLpcsVF+XdsTxfm6aJx+Ohd+/e5OXlsXz58qbXvF4va9asaUqShw4dit1ub7bN0aNH2bFjR/wm0uXABqAGUJscEZHwrc/KwjBNRp88Ge1QpINa0vjfOiAYzUBEROKYabHwfu/eDYm0TgK3yAMPPMC6devYv38/JSUlPPjgg6xevZrbbrsNwzCYOXMmjz32GIsWLWLHjh3ccccdJCYm8vWvfx2AtLQ07rzzTmbNmsXKlSvZtm0b3/jGNxg0aFBTF+9oaHUd19+Bh4AEoL7V4YiIdE5lTic70tK48sQJlnyhtCnWaPU4fq0GKoEbgZKoRiIiEt+29unDVR9+SM8TJziYkxPtcM4pVubr48eP881vfpOjR4+SlpbG4MGDWbJkCZMmTQLgRz/6EW63m+9973uUl5czYsQIli1b1nQPaYCnn34am83GjBkzcLvdFBUVMX/+/KjdQxraIJH+Gw3dtycC/2ztzkREOrF12dl8d+9eEv1+6mL4elU1G4tPPhrunnED8N9RjkVEJJ590rUrNQkJDNm374xE+t/+7d9C3t8f//jHtgqtmUg2GwvFvHnzvvJ1wzCYPXs2s2fPPuc2CQkJzJkzhzlz5oR07Ehq9e2vdgMf0zAxi4hI+NZlZ+MwTUaeOhXtUKSD+jswFOge7UBEROJY0GLh/V69VN7dybXJfaT/DlzfVjsTEemkShMS+DglhatOnIh2KOek+0jHt8WAF538FhFprW29e5NbWUnXsrJoh3JWmq8jr01y37/RcO/okW2xMxGRTmxtdjaXnzpFQiAQ7VCkA6oCVqFEWkSktXZ1747b4eCyOOjeLZHRJon0JuAYmphFRFprbXY2CcEgl8dwebfOcMe3vwPjgLQoxyEiEs/8VisfFBTE9G2wtCodWW2SSJvAGzR0AhURkfB9lpjInuTkmC3v1qQc//4B2IFroh2IiEic29anD13Ly8mtqIh2KGel+Tqy2uyy5r8BFwAXttUORUQ6qbXZ2Yw8dQq7yrslAo4Am1EVmYhIa33Yowf1NltMr0pL5LRZIr0SqEETs4hIa63LziYxEGBYeXm0QzkrneGOf38HrgYc0Q5ERCSO+Ww2dvbsGbOJtObryGqzRNoDLEHl3SIirXUgKYkDiYkxW94t8e9vQCowPspxiIjEu219+tDz5EmyqqqiHYq0sza9Y9XfaejcndeWOxUR6YTWZWcz+uRJbMFgtENpRtdIdww7gb2oikxEpLV29OyJz2rl0hhcldacHVltmki/CfiB69pypyIindCa7GxS/H4ujcHybk3IHcPfgesBI9qBiIjEMY/Dwc4ePRiyb1+0QzmDkujIatNEuhxYi8q7RURaa29yMp8lJKi8WyLmb0A3YFiU4xARiXfb+vShz/HjdKmpiXYo0o7aNJGGhjPcRUByW+9YRKQzMQzWZmcz5uRJ8PujHU0TlYl1HBuAk6i8W0SktUoKCvBbLFwaY6vSmq8jy9bSDVNTU1u03dvBIM6aGqa7XPT+0Y9CDmjw4MEhjwnHgQMH2uU40vEEwrwlUWJiYshjWvp790UpKSkhjwmGcR3u3r17Qx5z8ODBkMeE49ChQyGPeeedd0Iec8UVV4Q85tMQrqFa6PPxLz4frFsH49UWSlrOarW2aLs3g0FuNE0etlrD+jvg8/lCHtOrV6+Qx5w6dSrkMdLx3HXXXSGP+eCDD8I61tixY0Mec/LkyZDHlJSUhDymuro65DGZmZkhjwlHXl7onZLCeT/Z2dkhj/noo49CHnPRRRe1aLtqq5WPunbl0r172Th8eMjHkfjU5ivSBy0WPrBYuCaGVlBEROLR+w4Hn1mt8Npr0Q6lGZ3h7jj+bhhcDPTV919EpFW29upF/2PHSK6tjXYoTVRFFlltnkgDvGWzMcXnwxLmyp2IiIBpGCxJTITXX4cY6d6tSbljWQG4gev1/RcRaZXtvXphGgaDwqjYixTN15EVmUTabqcLUKDyaRGRVlmcmAhHj8LGjdEORTqgOsNgBXCDPhiJiLRKTUICu/PzuWTPnmiHIu0kIon0+xYLhwyDC3ftisTuRUQ6jWKnE3JzY6q8W2e4O5a/GwajgCz9PxARaZWtvXrR7/BhEuvrox0KoPk60iKSSGMYvGWzMWDXLtA3W0QkbEHDgJtuakik9fdUIuBNw8AArtHPl4hIq2zr1QsjGOTiGCrvlsiJTCJNY3l3ZSV5x45F6hAiIp3DzTfDgQNQXBztSHSNdAd0wjB4F7ghRq7DFxGJV5WJiezv2jVmyrs1Z0dWxBLp9VYr9U6nyrtFRFpr7FjIyIiJ8m5NyB3T3w2DItPEpf8XIiKt8kG/fgw4eBCnxxPtUDRnR1jEEmm/YbD7ggsayrtFRCR8djvceKPKuyVi/mEYJAIT9fMlItIq7/frhy0Q4OJ9+6IdikSYLZI7/3jAAAaXlNClooKKLl0ieSgRkY7t5pvhj3+EHTtg0KCohhKJs9Hxeoa7traWX/ziF6xcuZLS0lKCXyqP/vTTT6MUWWj2GAYfAtcHg7xhidg5dhGRDq8iNZUDubkM3rOHrRdeGNVY4nVujYRIzNcRTaT39OuH32plwK5dbBoxIpKHEhHp2IqKIDW1YVU6yom0fO473/kOa9as4Zvf/Cb5+fkYhhHtkML2hsXCncEgVtMkEMfvQ0Qk2j7o148pmzbh8Pnw2u3RDkeIzHwd0UTa63Syr3dvLvz4YyXSIiKt4XTCddfBq6/C7NlRCyNS10fF61nzxYsX8+abbzJmzJhoh9Jq/zAM7gNGmSbrlUiLiITtg/79ue6dd7ho/37e798/anFEam6Nxzk7EvN1xOu3Ph4wgIIDB3C53ZE+lIhIx3bzzbBzJ0S594Sal3wuPT2djIyMaIfRJrYYBp8B18Xp/wsRkVhxsksXjmRlMfiTT6Iah+brz0Vivm7xinRlZWXIO7///vv5oKCA60yTgg8/ZPvFF593TE5OTsjHSUxMDHnMddddF/KY66+/PuQxffv2DXlMOA4fPhzymNTU1AhEcqbXX3895DGPPPJIyGP2798f8pjk5OSQx8S6urq6kMckJCREIJKO77333gt5TDj/fwY1lnI7g0HWWiz84aqr+N/z/K0sKSkJ+TgSuv/+7//moYce4oUXXghrLooUv98f8pjc3FyWV1Vxo8/HL7p0gRasStfW1oZ8nHD+3lxwwQUhjykqKgp5TCAQCHlMOOrr60Me4/V6Qx7z5WsAWyKcckebLfQCx4suuijkMSKtceTIkZDHfPzxxyGPyc3NBaCPYfBve/aweskSvOo9EXWRmK8jWtoNUJ2SwsH8fC7evbtFibSIiJydx2JhXXIyk6qqzptIR4pKu2HIkCHNko09e/aQm5tLr169sH/pWritW7e2d3itssTp5Pb6ei4MBPg4jORIREQaLE9N5Z4TJxhdW8vqlJSoxBBPc2skRHq+bpdZ8sP+/Rn/7rvY/H78mphFRMK2PC2NXx06RDevlyMOR7TD6ZRuvPHGaIcQMesdDmoMg6kejxJpEZFW2Ot0stfhYFJVVdQS6c4u0vN1uyXSU9eupe+BA+xqp1JnEZGOaF1yMh7DoKiqij9lZUUlhs6+Iv3www83/fvb3/42t912G0VFRXHdsfs0r2Gw0uFgqsfDM0lJ0Q5HRCSurUhN5V/KyrCbJr4ozBGdvdlYpOfrdinYL83M5GR6OgOjfMG9iEi8q7Na2ZCczOQw+la0hUg1LomXSfnLTp06xbRp0+jevTv33nsv27dvj3ZIrbbU6WSI309eO10vLCLSUS1LSSEtGOTyMPpJtAXN15+LxHzdPle+GwY7+/fnoj17MOL0my8iEiuWp6ZyidtNrs8X7VA6vX/84x8cO3aMhx9+mC1btjB06FAGDhzIY489FlYDxFiwwuHAB0zxeKIdiohIXPs4IYFDdjuTqqqiHUqnF4n5ut1ayH3Yvz+ptbV0/+yz9jqkiEiHtCY1FV9jeXc06Ox2c126dOHf//3fWb16NQcOHODb3/42L774Iv369Yt2aGGptFh4125nahhdokVE5AsMg+WpqUyorsYShblOVWTNtfV83W6J9MGuXalJTORilXeLiLRKldXKpqQkJkapvFvOzufzsWXLFjZt2sT+/fubboESj5Y4nVzh9ZISxu2TRETkc8tSUsgMBBgaxu0vJTLaar5ut0TatFj4qG9fXSctItIGlqemclldHZlh3C+4NWLl7Pbjjz/O8OHDSUlJIScnhxtvvJFdu3Y12+aOO+7AMIxmj5EjRzbbxuPxcM8995CVlUVSUhLXX389hw8fDimWVatW8d3vfpfc3Fxuv/12UlJSeOONNzh06FBI+4klS5xOHMAErUqLiLTKDpeLozZbVMq7Y2G+jiVtPV+3693BP+zfn5yyMrJOnWrPw4qIdDirUlMBmNBJJ+Y1a9Zw9913s3HjRpYvX47f72fy5MnUfqmhy9SpUzl69GjT46233mr2+syZM1m0aBELFy5k/fr11NTUMG3aNAItbLTVvXt3rrnmGk6cOMHvf/97jh8/zvPPP8/EiROxWNp1im1TR6xWPrDZmKrrpEVEWsU0DFakpjKpurrde0XFwnwdKyIxX7frTSI/6dULr93OwD17WJuZ2Z6HFhHpUMptNrY0lnf/NSMj2uG0uyVLljT7+vnnnycnJ4fi4mKuuuqqpuedTid5eXln3UdlZSXz5s3jxRdfZOLEiQAsWLCAHj16sGLFCqZMmXLeOB566CFuueUW0tPTW/FuYtNip5P/qKuL2m1bREQ6iuUpKXyzrIxL3G62JyZGO5xOKRLzdbueLvfb7XzSq5fKu0VE2sDy1FSG19aS2o7l3ZEu7a6qqmr28LRwRbSy8XrxjC+dVFi9ejU5OTlccMEFfPe736W0tLTpteLiYnw+H5MnT256rmvXrhQWFrJhw4YWHfff//3fO2QSDbDE4SDVNBml7vAiIq2yLTGRk1Zru5d3x8rlWLEgEvN1u9edfdi/Pz2PHCE5SvdTExHpKFampmIFxldXRzuUNuH3+0lLS2v2ePzxx887zjRNfvjDH3LFFVdQWFjY9PzVV1/NSy+9xNtvv82TTz7J5s2bmTBhQlNyfuzYMRwOxxkTa25uLseOHWvbNxeHPrTZOGixcLXKu0VEWiX4hfJu4jAJlbNr90T6o759wTC4aM+e9j60iEiHctJuZ3tiIhM7yBlum81GZWVls8ePf/zj88bz/e9/nw8++IA///nPzZ6/9dZbufbaayksLOS6665j8eLF7N69mzfffPO8789QKTMYBkuczob7SeuDn4hIqyxPSaGbz8fA+vp2O2asrEbHUoPQthTRa6S/3PQFoBbYm5/PBR9/zOqz3LNr9+7dIR9n9OjRIY85ePBgyGM2b94c8hin0xnymHC0121W9u7dG/KYV199NeQxdrs95DHJyckhjxkyZEjIY4YPHx7yGIBHH3005DG+MEoq3W53yGP8YZQGd+3aNeQxFRUVIY8JRzhNI5KSkkIe09Ky40jakJzMv584QVIgQK3VGvHjRaqk6/Q+UxubqLXUPffcwz/+8Q/Wrl1L9+7dv3Lb/Px8CgoK+KTx8qK8vDy8Xi/l5eXNVqVLS0vDmlfinfcsHbrfsFj492CQi9xu3red+ZGhS5cuIR+npqYm5DEJCQkhjwnn7+dTTz0V8phw/PGPfwx5TDCMW5GFU1kRzvfNGsbfniuvvDLkMeF8D7Zu3RryGGj4exGqxDCuff1yMtESLW2G+EXhvJ/2Es6cUlJSEvKYCy64IOQxxcXFIY85l+2JidQZBpOqqvjQ5Wqz/X6VWCnBPt0gdPjw4fj9fh588EEmT57Mhx9+2Ozz19SpU3n++eebvnY4HM32M3PmTN544w0WLlxIZmYms2bNYtq0aRQXF4f1d6i12rXZ2GklvXtz3caNOLxevF/6BomIyLkZpkmh28246mrGVlczoL6eesMgvZ0S6Vhhmib33HMPixYtYvXq1fTu3fu8Y06dOsWhQ4eaPlAOHToUu93O8uXLmTFjBgBHjx5lx44dPPHEExGNP168a7NRbhhc6/OdNZEWEZFzS/P7uaqmhnE1NYypqSHRNEkL44RQvIuVBqFtLTqJdJ8+TH/nHS46dIj3+/aNRggiInHDFQwyqqaGsVVVXFVTQ5bfT4XVyrrkZP6Qnc07ycnUtGMSHckV6Za6++67efnll/n73/9OSkpK08pbWloaLpeLmpoaZs+ezc0330x+fj779+/ngQceICsri5tuuqlp2zvvvJNZs2aRmZlJRkYG9957L4MGDWqapDu7gGGw1GbjGp+Px9ppBUVEJG6ZJn08HsZVVzOuuppL3W6sQElCAvMzM1mdksLH7VSt2hBOZFakW7vf8zUI7dKlC2PHjuXRRx8lJycHOH+D0E6TSJ9MS+OzjAwGffqpEmkRkbPIDwSY4vUy2evlypMncZomnzqd/KNLF9akpPB+YiKBTnwd79y5cwEYN25cs+eff/557rjjDqxWKyUlJfzpT3+ioqKC/Px8xo8fzyuvvEJKSkrT9k8//TQ2m40ZM2bgdrspKipi/vz5USkRi1Vv2e38i89HQSDAAX1fRESasZkmo/z+pjm7T1kZbsPg3aQkfpafz5rkZE6GcclirPP7/VR9qUeL0+k872WtX9Ug9JZbbqGgoIB9+/bx05/+lAkTJlBcXIzT6YzJBqFRq9P6oE8friwpwRIMEgzzJtgiIh2FYZpc8oWJeJDfjx94127n17m5rE5J4VA7nsU+l1i57cX5YnC5XCxduvS8+0lISGDOnDnMmTOnrULrcN6226kHrvH5mKtEWkSELsEgE30+pni9TPD5SDNNPrNYWGa38/9yctiUlIQnBvKbSM7Xq1at4ic/+Umz5x5++GFmz579leNONwhdv359s+dvvfXWpn8XFhYybNgwCgoKePPNN5k+ffo59xfNBqHRS6R792bqli30+ewz9pynQYyISEfkMk2uakycJ3m95AWDVBgGKx0O5rhcvO1wUGWxkJuVFe1Qm0S62ZjEnlrDYE1jeffcMJp+iYh0BP0aK8WmeL1c7vdjA7ZbrfxPQgJLHQ4+sFrBMMj9QtVTtEVybh0/fjyvv/56s+fOtxrd0RqERi2RPpSTQ3lSEoP37VMiLSKdRl4gwOTTJdteLy5gr9XKIqeTJQ4H79ntnbpkW2LTW3Y7T7ndpAeDlMfAKouISKRZg0HGNK46T/F66RsM4gbW2u38KCmJZQ4Hxzrx30ObzdbiO2101Aah0WvBaRjs6N2bwZ9+yutXXAH64CgiHZFpcrHXS1FtLRNqayn0evED79nt/DIpiaUOB3vjrBuyVqQ7nyV2O0+73Uzx+1mou22ISAeV7PMx7ORJRp44wbCTJ0nx+zlmGCxzOHjI4WCt3Y47jnKWWGk21lEbhEb109sHffpw5Y4ddD11is9iqHRRRKQ1nMEgo91uJtTVMb62lrxAgCqLhbWJifzW6WSlw0FFJz6LLfGn1GJhi9XKNT6fEmkR6VC61tYy8sQJRp44QWFFBVbT5JOUFP7esyfzT53ifasVM46S51jUURuERjWR/qRbN9wOB4P27VMiLSJxLdvvZ0JdHRNqaxntduMyTQ7YbCxOTmZlYiJbXC78hnFGh8t4o2ukO6+37Hb+q76eBNOkXh8qRSROWYJBBlZUNCXPPerq8FosbMvI4LcXXsim7GxONvaD2N54m6Z4FCvNQaHjNgiNaiIdsFr5sKCAwZ9+ytLhw6MZiohIaEyTgR4PE2prmVBXx2CPhwCwNSGB36Sn83ZSEnvt9g552YoS6c7pLbud2fX1jPP7WdIBb+UiIh1XUiDA2GPHGHHiBMNPniTV56PM4WBTdjbzLriArRkZeOLsMquW0NwaWVH/ifmgd2++vWwZXaqrqYihLnciIl/mCAa5vLaW8TU1jK2uJt/vp8YwWJuYyJ/S0liTmEi5bg8kHdQnViu7LRau9vmUSItIzOvm8XBVRQVXVlZyWXU1NmBvSgpv9OjBpuxsdqemqmRbWiXqifSHBQX4LRYG79vH2sGDox2OiEgzmX4/V1VXM66mhtE1NSSaJofsdlakprLc6WSzy4WvE03EKu3u3Bbb7Xzd6+U/TZNgJ/q5F5HYZzFNBtXWNiXPferr8RoGm1NS+FWPHuzo1YsTLle0w2xXsdJsrKOKaCJ96tSpFm33YU4OF+7axaJu3SIZTjN//etfQx6zdu3akMc888wzIY8Jx6hRo0Ie05LW818WzvWdmzZtCnnMV914/VymTJkS8pjTXQND8eCDD4Y8BmDnzp1hjQtVYmJiyGPuuuuukMfs3r075DGn7wUYaQcOHAh5TJ8+fRr+YZoM8HopqqtjQl0dl3g8AGxzOnk2PZ23ExP5pLFke+/evW0ZtkjMe9Nu5wceD8MDATa1YxlkQhj3r7700ktDHlNcXBzymHDU1dWFPCaczyAHDx4MeUx9fX3IY2688caQx2zZsiXkMeH8bT958mTIYwCeeuqpsMaF6uqrrw55zMyZM0MeE+73oT0cOnQo5DGn/yYkBQJcWVfXUClWW0t6MMgpq5XVSUk8lZ7OhqQkahube+4K43dI5KtEfUUaYEv37nyruJhErzfaoYhIJ+QArmxMnIvq6ujWWLK9LjGRl1NTWZ2YSJlKtptoRbrz2mK1ctwwuMbna9dEWkTktB5+P1PLyxlfW8vwujocwC6Hg1e6dGFVUhIfJCSoYqaR5tbIiolZcGu3bvzbli1c+tln0Q5FRDqJbOAa4DpgMpBy7BiHbTZWJiayMjGR91wuvJqIzxBLXUCl/ZmGwRK7nWt8Ph4OY5VYRCRUFtNkiNdLkdvNRLebAT4fXuC9xER+mZ3NquRkjqhvw1lpvo6smEikyxIT2ZuRwdAjR/go2sGISId1MQ2J83XAyMbnNgGPA9u6d2d3B+2yLdKW3rLbud3r5YJgkCPRDkZEOqSkYJCr6uspcrsZ73aTFQxyymJhlcvFM2lpbEpLo1aVYhJlMZFIAxR368a0jz5il89HUGeVRKQN2IGxfJ489wZqgGXAncBbQGnjtn0cjmiEGJdU2t25rbHZqAGu8fn4Q7SDEZEOo5vfz0S3myK3m5H19TiB3XY7f0lOZqXLxVaHo6lkO0FJdIuo2VhkxU4i3b07M0pKyP/4Y44MGhTtcEQkTqUHAkygIXGeAqQCB4E3Gh+rAU/UohOJfx7D4O3G8m4l0iISLsM0uai6mlsrKihyu7nI58MHbHI6ebxLF1YkJnJIvRgkhsXMT+ehtDSOJyfTc/t2JdIi0nKmSX+/n0n19Uysr2eo14uVhpLtJ2hInj+IboQdim5/JdBQ3v0/dXXkBAKUamVIRFrIFQgwrLycMWVljCovJ8Pno9xiYVVCAnNSU1nrclHd2GVbWk9za2TFTCKNYVDcrRtjt2/n3dtuA/0Sicg52E2Tyz2epuS5VyBAnWGw1unkvi5dWFBRwfFoB9mBKZGWpTYbfmBSfT0vJSVFOxwRiWE5Hg+jy8oYU1bGkIoKnKbJfpeLJTk5vJORweKKCgLqTxIRKu2OrNhJpGm4DdY1u3aRtX8/J0/f01VEBOgSCDDB42FifT3j6utJNU2OWiwsd7l4KCGBDU4n9Y0T8fGKiugGK9LBVVgsvGuzMcXtViItIs0YpsmFNTVNyXP/2lr8hsH7qak816sX72RkcMTlato+UFkZxWhFwhdTifTurCzqk5Mp2L5dibRIZ2ea9P1CyfbwxpLt7XY7zyUnszwhgZ3qst3uVNotp71pt/Mzt5ukYJBaVZGJdGoJgQDDKioaSrbLysj0+aiy2diYns5L3buzKT2dGl3v3O40t0ZWTP1EBy0WDg0eTM9t2yiePj3a4YhIO7OZJpd7vUx0u5lUX0/vQAC3YbDO6eTHXbqwMiGB47oeUyQmvGWz8QtgnMfDm19YXRKRziH7CyXbl1VW4gwGOehysSwnhw0ZGZSkpqpkWzq0mEqkAQ4OGUL/DRtIOX6c6tzcaIcjIhGWFgwyvnHVeVx9PV1Mk2MWCysTEpidkMA7Tif1Wu2KKVqRFoBDVisfNpZ3K5EW6QRMkwtqaxlz6hRjysoYUFuLHyhJTeUPPXuyITOTQ/pbEFM0t0ZWzCXSRy6+GL/dTsH27eyYMiXa4YhIBPT2+ZhUX8+kxpJtG1Bit/N8Y8l2id2OqbPYMUml3fJFS10uvl1Tg8008et3VqTDcQQCDKusZHRZGaPLysj2eqm2WtmUns7Cbt3YlJ5Otd0e7TDlHNRsLLIimki//PLLIY956623ODVkCP127uTk7be3aIzf7w/5OCLh+vjjj8Mal5ycHPKYQCAQ8pgBAwaEPGbv3r0hj3nttddavK3VNBkZCHCf389Uv5/+wSD1wBqbjf9KSGCJzcZnp1edPZ6GR6NKNSERaRfl5eUhj7l98GC6lJRwXVoaW9LSWjTGFsZ1ku01z+/cubNdjiPtxxrjlwP17ds35DFpLfxd+6L9+/e3eNtsv59xNTWMr6lhVF0dLtNkv93Om8nJrMrNpdjl+vzE2Zeae4b7GUkkHsXcijTA8VGjGDRnDo7KSrxh/LEQkehLM00m+v1M9fmY5PeTDhw3DJbabDzkdLLaZqNOK1hxSWei5bTdiYkcczi4sry8xYm0iMQY0+Qij4fxNTWMr61lUH09AWCry8WcrCxWJSWxz+FQc884pPk6smIykS69/HIwTXLee4/DkyZFOxwRaaE+wSBTfT6u9vsZFQhgBz6wWPiDw8Fiu51tFotKtkU6EsNgXXo6V5aV8XRBgT5oi8QJRzDIyLo6xtfUMK62lny/n2qLhXVJSfypSxfWJSdTEeOr+SLRFpOJtLdLF8ovuojcd99VIi0Sw6ymyeWBAFc3lmwPCAbxAGutVu5vLNk+rEZhHYqukZYvW5eRwS3Hj3NBXR27dU9pkZiV5fcztnHVeXRtLYmmyUG7nWUpKaxKSqI4MRGfToZ1KJpbIysmE2mA4yNHcsFLL2GpryeYkBDtcESkkcvr5eLDhylyu5ns95NhmpQ2lmz/3Olklc1GrSbiDkuJtHzZ1pQUqq1WriwrUyItEktMk97V1Yw8cYJLDx/mksaS7e0uF7/LzGRVcjJ7VbLdoanZWGTFdCJ90R//SNa2bZSOGhXtcEQ6tazqagYfOMAlBw9ywbFj2IJBdlgs/NFuZ7HNRrHVqpJtkU4qYLGwoUsXriovZ16PHtEOR6RTsweDDD51ipEnTnD5iRPk1tdTa7WyzuXi5bw81iQlURFGwz8ROVPM/ibVdetGdY8e5G3cqERapJ0ZwSB9Tpxg8MGDXHLwIN3Ky/FbLOzKz+cvI0bwQc+ezF28ONphSpRoRVq+bF1GBlM++YS8+nqOqYpMpF2leTxcfvIkI0tLuezUKVyBAMcTEtiYk8PG7Gx2ZGTwycGD0Q5TokBza2TFbCINDd27eyxZAoEAqOGBSEQ5vV4uPnKESw4eZNChQ6TU11OdkEBJjx7847LL2NmtGx6HI9phikgMejctDa9hcGV5OX/Nz492OCIdm2nSq6aGESdOMLK0lAGNt4nclZbGwj592JSdzf7kZJVsi0RYbCfSI0bQ7y9/If2jjygvLIx2OCIdTo7bzYjSUsbt2MEFR49iDwY5kp7OugED+KBnTz7NzsZUszD5Al0jLWdTZ7NRnJrKVUqkRSLD48G+YQOOpUuZv24deW43bquVrZmZPF1YyHtZWVQ6ndGOUmKM5tbIiulEuvKCC6jPyCB340Yl0iJtwDBNBlRWNp3F7l1Tg88w+CQ/n1cvv5wPevbkZGpqtMOUGKdEWs5mXUYGP9y3j1S/nypdgynSasbJkzhWrMCxdCn2Vauw1NYS6NaNzVlZbMzO5oOMDHyq2JSvoGZjkRXbM53FQumIEeRu3MjHd96pEhWRMCT4/Vx26hQjGhuPpHu9VNrtbM7O5uW+fSnOyiKzd+9ohykicW5dejo/2rePUeXlLM3OjnY4IvHHNLHu2oVj6VIcS5di27IFwzTxXXYZ7nvuwTtlCoGLL+a33/xmtCMVEWI9kaahe3fPxYtJPnCAml69oh2OSFzIdrsZceIEI06c4JJTp3CYJgeSkljetSubcnL4qEsXgl84MZUZxVglvqi0W87lpMPBzqQkrlIiLdJyXi/2d99tSJ6XLcN64ABmYiLesWOpeeopvJMmYebmRjtKiVOaWyMr5hPpU5dcgt/lInfjRiXSIudgmCaDvV6K3G6u3bCBvtXV+A2DkvR0/njBBWzKyeFoYmK0wxSRDm5dRgbfOnIERzCIV/0VRM7u1Cl46y144w0ylizBUl1NoGtXvJMn4508Gd+VV4K634vEvJhLpJOTk894rnzECPI3b+b4d75z1jF79uwJ+TjHjx8PecyAAQNCHvO3v/0t5DHhCOc4gUAg5DGfffZZyGPCUVZWFvKYZ599NuQxkyZNCnlMrHAFg1xRX0+R280Et5ucQIAKi4Wtubm80rs3xVlZ1NrtLdrXH/7wh5CPH84Y6Rh0hltOc7lczb7elJ/P/zl0iNH19WzKPHutiy2M66fDGTN//vyQx7jd7pDHhCM9PT3kMevWrQt5TE5OTshjwvHCCy+EPObEiRMhj8nKygp5zPjx40MeA1BTUxPymLN9hgXANOHjj+GNNxoeGzZAMAjDhuH+3vfwTp5MYNCgFl/CuGTJkpBjk85J83VkxVwifTanrriCCx95BMeJE3hVLiadWJ7fT5HbTZHbzRi3Gyew12ZjUVISK10uip1O+l90UbTDlA5Mpd3yVQ4kJnLY5WLMiRPnTKRFOgWfD9at+zx53rsXXC6YOBH+53/g2muha1fcJ09GO1LpwNRsLLLiIpEuGzmSoNVK5jvvcPTGG6Mdjki7MUyTQq+XiY3Jc6HXix/Y7HTyRHo6K10u9rdw1VlEJOIMg3eysph4/DiGaWKqSah0JmVlsHhxQ+K8ZAlUVkLXrjBtGjzzDBQVNSTTItIhxEUi7U9JofLSS5VIS6eQEAwy2u2mqK6O8XV15AYCVBkGq10unktNZXVCAlW63YVEkVak5au8k5XFrYcOcWFVFR+lpUU7HJGIMnbvxrZ4MbbFi2HjRggE4LLLYOZMuO66hn/rhJJEiebWyIqLRBrg1Jgx9Pnd77DW1BA41zUoInEqx+9nQl0dE+rqGON2k2Ca7LfZeCMxkZWJiWx2OvFrIhaROPBhWhrldjtXnDypRFo6Hr8f67vvYn3rLWyLF2PZuxczIYHA2LHw2982rD536xbtKEWkHcRPIn3FFfT7zW/IeO89TkyYEO1wRFrHNLnY62VCXR1FtbUM8noJAMUJCTyTns7KxEQ+tdsJBIPRjlSkGV0jLecTNAzezcpi9MmT/KFv32iHI9J6FRXYli/HungxtuXLMSoqCObmEpg6Fc+jjxIYNw6Sks7dbEwkSnSNdGTFTSLtyc2lun9/MtevVyItcckZDDKqvp6i2lrG19WRHwhQbRisTUzk+bQ01iQmUqGSbYkDSqTlfDZkZXHN0aP0qKvjkG69J3HI2LMH2+LFWBcvxrphA0YgQGDwYHx33YX/6qsJDhkCusWbxDjNrZEVN4k0NKxKd//LXzB8Pkw1WJI4kHW6y3ZdHWPq60k0TQ7YbCxJSuLtpCQ2JyTgU8m2iHQwxenp1FssjD5xglcKCqIdjsh5WYJB+p44waWHDnHpoUMkzZ+P6XQSGDsWz69+RWDqVMzu3aMdpojEkLhLpHs9/zxp27dTMXx4tMMROZNpctHpLtt1dVzSWLK91elkTno6bycmssduV+MRiVsq7ZaW8FitbMnIYMzJk0qkJWa5vF4Kjxzh0kOHGHzkCMkeD5UJCbzfowfpc+cSGD8eVK4tcUxza2TFVSJd27cv9Xl5ZL7zjhJpiRmO0yXbbjcT6uro1liyvc7l4k+pqaxyuSi3WrGoBEykzTz++OO8/vrrfPzxx7hcLkaPHs0vf/lLBgwY0LSNaZr87Gc/47nnnqO8vJwRI0bw29/+losvvrhpG4/Hw7333suf//xn3G43RUVF/O53v6O7Vp5a7Z2sLP7r449J93opdziiHY4IANnV1Vx66BCXHDrEgGPHsJkmB9PTWTVgANt79GBfVhamYXDZdddFO1QRiXHx9cneMDg1ZgxZ69eDzrBIFGUFAtxSXc3/lJay9dAh5peWMs7tZnliIt/MzWVoz57cnZPD68nJlOu6Z+lgTq9Kt/UjFGvWrOHuu+9m48aNLF++HL/fz+TJk6mtrW3a5oknnuCpp57i2WefZfPmzeTl5TFp0iSqq6ubtpk5cyaLFi1i4cKFrF+/npqaGqZNm0YgEGiz71dn9W5mJiYw6uTJaIcinZgRDNLv+HG+tmULj/ztbzzx2mvcsmULQcPgz5dfzr1f+xoP33ADr192GZ9mZ+ve59KhRGq+DnXOfvzxxxk+fDgpKSnk5ORw4403smvXrjNinT17Nl27dsXlcjFu3Dh27tzZbBuPx8M999xDVlYWSUlJXH/99Rw+fLjV36dwxdWKNMDJK66g22uvkbx7NzVfWHkQiSjTJO/ECS7as4dbP/uMS7xeALY5nfw2LY2ViYnsVsm2dAKRKu0O1ZIlS5p9/fzzz5OTk0NxcTFXXXUVpmnyzDPP8OCDDzJ9+nQAXnjhBXJzc3n55Ze56667qKysZN68ebz44otMnDgRgAULFtCjRw9WrFjBlClT2v19dSRVDgc70tIYc/Ikb3XtGu1wpBNJDga50u1motvN5FdeIcXjoSohgfe7d2fRkCHs7NqVevXakU4gFuZr+Pzk9/Dhw/H7/Tz44INMnjyZDz/8kKSkJODzk9/z58/nggsu4JFHHmHSpEns2rWLlJQUoOHk9xtvvMHChQvJzMxk1qxZTJs2jeLiYqxRWLiKu0S6cvBgfCkpZK5fr0RaIsrq99Pn4EEu2ruXi/bsIb2qCo/Dwds2Gy+lpLA6MZFTWm0WaVNVVVXNvnY6nTidzvOOq6ysBCAjIwOAffv2cezYMSZPntxsX2PHjmXDhg3cddddFBcX4/P5mm3TtWtXCgsL2bBhgxLpNvBOVhbf2bePBL+felvcfeSQONLd76eoro4it5uR9fU4gI/tdtZceCHbe/Tg06wsTF1iJRIVHfXkd8zNag8++OB5t/lJQgL9Xn2VBz/9FIBvfetbIR8nJycn5DFvvvlmyGPaS01NTchj7GGcjW2veyTW1dWFPGb//v0hj1m6dGmzr9P9fq6orGRsVRUjq6tJCgY5arezNC2NtX36sCU5mffef5+rQz6SSMcRqWZjfr+ftLS0Zs8//PDDzJ49+7xjf/jDH3LFFVdQWFgIwLFjxwDIzc1ttm1ubi4HDhxo2sbhcJCenn7GNqfHy1c7cuTIV77+eiDA94JB+uzZw/LGFYU+ffqEfJzEMG6hFU55fmlpachjwvHxxx+HPMYWxomI0yeYQhHO77cRRjXWl3/vWuL48eNN/7aYJsNNk2sCAa4JBik0TbzAOouF+2023rJYOGCx4P7gA6aFfKT2k5WVFe0QpAOL5Iq03+/v9Ce/Yy6Rbol1mZlMOXGC/Pp6jiYkRDsciWemSd/6esZWVXFVZSWDGhP4HYmJPJ+Tw9q0ND5JSFDJtkg7sNlsnDp1qtlzLZmQv//97/PBBx+wfv36M1778gd80zTP+6G/JdtIyxxyONjlcFBUW9uUSIuEK8k0mRgMck0wyNRAgBzgJLDUYuFxm43lFgvV+t0VaRerVq3iJz/5SbPnOtvJ77hMpN9LT8djGFx56hR/6dYt2uFInLGbJmP8fqb6fFz70Ud093pxWyy8m5LCz3v0YF1qKmW6dkrkrCJ9+6vU1NSQxt1zzz384x//YO3atc06befl5QENE29+fn7T86WlpU0TdV5eHl6vl/Ly8mYTc2lpKaNHjw77vUhzbycn8y8VFVhNk4CSHAlRV7+fSfX1THS7GeXx4AQ+MgxetFp5y2plo2EQ1M+VyFlFakXaNE3Gjx/P66+/3uz5znbyOy4TabfVSnGXLlyhRFpaKCMYZFJj8jzB5yMVOGIYvJORwdq0NDYnJ+PVtVMi5xUr95E2TZN77rmHRYsWsXr1anr37t3s9d69e5OXl8fy5csZMmQIAF6vlzVr1vDLX/4SgKFDh2K321m+fDkzZswA4OjRo+zYsYMnnniiDd6VAKxITuY/ysoY6nbzXhgl2tK5GKbJJV5vQ/JcX89Anw8fsMnp5Cc2G29aLOzTfC3SIpEs7bbZbJ3+5HdcJtIA6zMzmbVnD6k+X7RDkVhkmgwIBpnq8zHV52N4IIAVKLZamZOQwBKbjR1WK7169Ih2pCIShrvvvpuXX36Zv//976SkpDSVdaWlpeFyuTAMg5kzZ/LYY4/Rv39/+vfvz2OPPUZiYiJf//rXm7a98847mTVrFpmZmWRkZHDvvfcyaNCgpkYm0no7nE6O2WwU1dQokZazcgWDXOXxMNHtZkJ9PTnBIBWGwdsuF8+mpLA6IYFqi4WTupWaSFzqqCe/4zaRficjg3uB0WVl0Q5FYoTNNBnduOo81e+ndzBIHbDaZuM/XS6W2e0c11lskVaLhRXpuXPnAjBu3Lhmzz///PPccccdAPzoRz/C7Xbzve99j/LyckaMGMGyZcuabqMB8PTTT2Oz2ZgxYwZut5uioiLmz58fldtodFiGwdtJSUysqeHx7OxoRyMxIt/vZ2LjqvPo+noSgD02G68nJrLC5WKLw6FLAURaIVZuVwkd9+R33CbSZQ4HO1NSuKKsjN3RDkaiJtXv54qqKq6qqmJUZSWpwGeGwVK7nfvtdtbabNRrIhbpcFry4cAwDGbPnv2VjU8SEhKYM2cOc+bMacPo5MtWJCfz9cpKBni9qI6sczJMk0sDAaZ4vUz2ehkcCOAH3nM6eSItjeUJCexXfxKRDqmjnvyO20QaGsq77zh4kL1eLwGHI9rhSDspaOyyPbaykktqa7ECH7pc/LaxZPsDq1VdtkUiJFaukZb48l5iItUWC0U1NSw5/+bSQbhMk7E+X1PynGuaVBgGK+12nktLY3VCApWqFBOJmEg2G2vr7ePx5Hd8J9IZGfzH/v10++gjDl5ySbTDkQixmSaX1NQ03KKqqooCj4d6w2BTSgqPde/OutRUTjgcYd1HWkRCp0RaQuUzDNYmJSmR7gTyAgGmNCbPV/p8uIC9FguvOZ0sdTjYZLPhNwxcLle0QxXp8DS3RlZcJ9IHExM54HJR8P77SqQ7mBS/nzHV1VxVWcmY6mpSAwFO2GysS03l6a5d2ZSSQr3OYouIxI0VSUk8XV1NltvNSSVRHYdpckljyfYUr5dLAgECwCabjV8kJrLU4WCPeg6ISAcU14k0NKxKT3//fdYHg5hKrOJan0CAqxubhY3asQMb8LHLxcKsLNakpfGRy4Wpkm2RqIql5iUSX9YmJeEFRpWW8kZBQbTDkVZIME2u8vmaVp7zTZNKw+Btu525Lhcr7HYq/n97dxocV33uefzXe2tfWla3hGTFYAMGLxg7tllu8G5rgp3gDFAhuSE1VIopqrjlAooXoZJ4pmpMFTMBbuxJqDvjGxPAgVvFUhhCsA3CxBi8BRmbzXa8YBtJrd2W1K3ezrxA9CCwcR+pj3rR91PVBTp9zumnu6g8efp5/v/m/5MBGUe+tlbuF9I+n35y5oyqjx1T2+TJmQ4HJjgMQ/Pi8eQu25cnEgpLetvp1CN1dXq7tFRtrH0HgLzQ53BoT2Gh5lNI5yR/IqFlkYiWR6O6KRpVoaRjdrtecru11ePRu0Mj2wAwXuR8If1RSYkGSkvV0NxMIZ0DSg1Di4cK56WxmCoNQ202m7Y6nfpvXq+anE4N2GyaVFWV6VABXABrpDFSbxQX66H2dhVFo+pnh+bsZhiaPrTeeUUkollDI9t7nE49UlCg191uHbbbJZtNNgpoICtly2Zj+SrrCum//e1vpq/55J/+Sd9pbtaeH/0o5d2aDx48aPp1ysvLTV9z8uRJ09eMRCgUMn1NJBIxfU1lZaXpay41DC3q79fC/n59NxSSS9LHbreeLS/Xm0VF+sDjSY5sFw093n//fdOvA8B67NqNr/r000/NXXD6tFRfr++2t+ut2lprghpy9uxZ09fU19dbEMk3ffXnXFJl+rOWTP8kjMcwdFMioWXhsBaHw6qNx3XOZtNbXq82eb1q8nrV/ZV7Fg7989SpU6ZjA2A9cqu1sq6QHomT11yjK3fuVHlrq3pqajIdzrhnNwxdOzioxQMDWhwKaUo0qoik9woL9T+qqtRUVKTP6UQAwPhTV6fDpaWaHwxaXkgjNRMSCS0ZHNSywUHdFI2qyDB00uHQawUF2u71arfHoygdZwD4hrwopM9ceaWiHo8ampsppDOkOJHQ90IhLR4Y0IJQSJWJhDrsdjUVFupxn0/vFBaqn41HgLxARxqj8a7fr1uPHZMrkVCUvDD2DENXxeNaNlQ8z4rFJEn7XC49VliopuJiHXY6U57wA5C9yK3WyotCOu526/RVV6mhuVkHGhszHc64UReNaslQ8Tw3HJZb0icul/5cUqI3Cgp0wONRwmaT05kX/5kBANLgvepq3XnkiGZ0dmr/hAmZDmdccBuGro9EtCwS0bLBQdUnEuqz2dTkdmtTQYHe8HjUOfSlhpuJMQBISd5UOCeuuUYL//hHFfT0KDSCtcy4OLth6NpwWIuHiucrhka2d3u9WldZqe0FBTpDAgbyGmukMVoniovVUlCg64JBCmkLVSUSWhEKaVkkogWRiIoNQ5/Z7drq8eh1t1vvut2K0HUG8hqbjVkrbwrpUzNmKGG3q+GDD/TJ976X6XDyRrFhaFEsphWxmJbFYppw7py67HY1FRToX8vL9beCAvUxmgeMKxTSGBWbTe9WV+t7ra3631ddldxsEqNkGJqaSGh5JKIV0ajmxOOSpL87nfpdYaG2ejz62OFgZBsYR8it1sqbQnqwqEitkyerobmZQnqU6hMJrYjF1BiN6sZ4XB5JH9vtetrl0juVlXp/aGQbAICReM/v1+qTJzWlt1eHmSIbMZdh6IZY7IufqIpG1ZBIqE9Sk8ulf/F49FZhodr5shsALJE3hbQknZg1S/Oef16ucFhRrzfT4eQMm2FodjyuxqHO87REQlFJ7zgc+rXHo7+6XDoxlIgr+VyBcY3RbqTDh+Xl6nW5dF0wSCFtUmUioaVDhfOiaFQlkk7bbHrd7dZfXS7tdDo1OPRlt4siGhjXGO22Vl4V0idnztT1zz2nug8/1PHZszMdTlYrMgwtiUa1IhbT8lhM1YahLptNW51O/S+PR284nTpL1xkAYIGE3a49EyZofjCoJy+/PNPhZDfD0HdCIS0Ph7UiEtF343E5JO13OLTe69VfXS4dYmQbAMZcXhXSfVVV6qyrU0NzM4X0eVySSOg/xeNqjMf1vXhcXkmf2u36s8ul15xO7XE4FCcRA7gIvolGOrzn92vp55+rpr9fLUVFmQ4nqzgTCV1z7pxu7O7WDd3dqhsc1ICkt1wu3efxaKvLpTa6zQAugnxtrbwqpKUvdu+e9uabssViMsb5zy7ZDEPXfqV4npFIKCbpHbtdv3G7tdXt1jESMQATGO1Guuz3+TRot+u6YFAvTJqU6XAyrjQa1XU9Pbqxu1vze3tVHI8r6HbrnfJyPVZRoT+3tSnMl90ATCC3WivvKs2TM2dq9iuvqObIEX0+dWqmwxlzhYahhUOF84p4XAHDUJekbQ6HHvV4tM3hUO9QInZQRAMAMmTQ6dT7Pt/4LaQNQw3hcLLrPOPcOTkkfVxUpD/X1GhnRYUOFxYmR7bDwWBm4wUADJN3hXTnxInqq6hQw4ED46aQDsRiWh2NqjEe103xuAokHbbZ9JzDodecTr1rtzOyDSAt6Egjnd71+/Uvhw6pLBJRr9ud6XAs50gkdHV3t645dUo39vSoPhzWoM2mvWVl+p+TJumdigp1jIPPAcDYYLMxa+VFIb1ly5Zhf9eUlWnuu+9qy1e+yf264uLisQjNEjbD0PRoVEtCIS0JhzU9GlVM0rt2u/67263XHA4dTaHb3Nvba32wAAB8hcPhSP77vkBAtkOHNL+jQ9vr6y94TU9PzxhEZo3iaFRz2ts1LxjU7I4OlcRiane59E5FhX7X0KC9paUa/MpnciFdXV1jEC0AIFV5UUh/3e5AQDefOKFLz57VsbKyTIeTFt5EQjcODmpJOKzFoZACiYR6bTa95fXq30pK9FI4rG66zgDGAB1ppEuvx6OPKyo0v63tWwvpXHNJf7/mBYOaFwzq6p4eOQxDR0tL9XJDg3ZXV2tfIiGDnA3AYuRWa+VlIX3I51Of06l5ra05XUhXx+NaHAppaTisGwcHVWAYOuFw6JXCQm3zerXH41FsKBF3Dw5mOFoA4wGj3Ui33X6/7jh8WJ54PKXObDayJxK6uqcnWTzXDQwoYrer2efT76dO1Z7qanV4vcnzjRzusAPIHeRWa+VlIR2327XP79e81lb9+YorMh1O6gxDV0ejWhIOa2kopJnRqOKS9rvderS0VNu9Xh11OvmtSABA3ngvENB/+eQTXdPert2BQKbDSVlRNKo5HR2aFwxqTkeHSqJRdbnd2lNdrY1XXKH3fT4NjvNfDwGAfJa3/wu/OxDQgjNnVD0woGBhYabDuSCPYej6cFhLw2EtDodVG4/r3NDI9r8XF+tNr1c9OfoNPYD8REca6fR5UZE+Ky7W/La2rC+ka/r7NW9ovfO07m45DUP/KCnRlokTtXvCBB0pK2NkG0DWYLMxa+VtIf33CRMUtdk0t7VVr1x6aabDGWZCPK5FQ13nfxocVKFh6KTDodcKCrTd69Vuj0dREjEAYJzY7fdr2alTshuGElmU/+yJhKb29CSL54n9/YrabGr2+fTE1KnaM2GC2gsKMh0mACAD8raQDrlc+qCqSvPa2jJfSBuGpkajWhoOa0kopFnRqBL6YmT7X0tKtL2gQIcZ2QaQA1gjDSu86/fr1n/8Q1O7u/VhZWVGYymMRjW7o0Pz2tv13fZ2lUaj6na7tXfCBD15+eX6u8+nMCPbAHIAudVaeZ0JdgcCuvvQIRVHIuob499ldMXjmt7ZqVu7u7UkHNYl8bj6bDbt8Hr1p6GR7S5GtgHkIBIz0u1Iebk6PR7Nb23NSCFdGw7rxu5uLThyRNO7uuQ0DB0vLtZf6uu1u7panzKyDSAHka+tldeF9J5AQPccPKg5waDeqquz/PXKBgf13WBQc9vaNKujQwXxuE45HNrq9Wp7QYHe83gUIREDADCMYbNpj9+veW1t2jh1quUTWnbD0LS+Pt3Q3a0bu7t1aSikqM2mDyor9W9XXqndEyZk9f4qAIDMy+tCusvr1afl5ZrX2mpNIW0Yajh3TnODQc1ra9PlQz9ncbi8XP8xebL2VFfr5WPHGNkGkDcY7YZV3vP71fjZZ5rY16fPSkrSfv/CeFzzenp0Y3e3ruvpUUUspm6nU7vKy/V/6uq0p7xcbp8v7a8LAJnCZmPWyutCWvqiK/2fjxyRKx5XNA2j1M54XNO7ujS3rU1zg0H5QyGFHA69X1Wl382Yob3V1er1eP7/BRTRAABc1AGfTwMOh+a3tqatkA6Ew7pxqHi+9uxZuQxD/ygo0Jbqau2sqNCHxcXDNjcb20VgAIBclveF9HuBgP75k080o6ND+/3+Ed2jJBzWNZ9/rhknT2pWe7sK43EFvV7t8fu1p7paB32+tBTpAJAL6EjDCjGHQ/urqzW/rU3/MWXKiO5hNwxdcfasrhvaZfuyoZHt90tLtX7iRL1TUaHPvd40Rw4A2Yncaq28L6RPFRfr86IizW9tTb2QNgxdcvasrj1zRteeOaMpHR2yS/q0vFzPX3aZ9vj9Ol5SQrcZwLjDaDestNvv1wPNzfKFw+pMseD1xmKa09Wl+R0dmtfRoYpoVL0ul3aVlenf6+q0u6xM/eyyDWAcIrdaKy8yS91F1j9/NHmy5h85opcvuSS56+Znn3027BxnIqGru7o0r61Nc9vaFAiFFB4a2V4/fbr2Vlfr6e3bdYWkf7bqjQAAMI7tra5W3GbT3LY2vdbQkDyeSCSGnVcdDuu6jg5d19mpa7q75TYMnSgs1F9ravSuz6ePysq0valJjWP9BgAA40ZeFNIXc2DSJC05cEANbW06EQgkj5dEIpoTDGpuMKhr29tVFIupw+vVnupq7fH79YHPpwgj2wAwDB1pWKXf5dLBykrN/1ohbRsa2b6+s1PXdXTosv5+xWw2HSwr0/+97DLt8vn0ObtsA8AwbDZmrXFRSB/z+3XO69XMEycU9ng0/cQJXf7pp7qyu1sOSYfLyvTipZdqT3W1jpWWMrINAECGvOf3666PP1ZFOKwreno0t61Nc4JBVUYiOut0ao/Pp80NDdpTWal+lyvT4QIAxqlxUUgbdrsOfuc7WtrcrOXvv6+I06n3Kyv1+6GR7S42HgGAlLBGGlbb7ffrv370kTa98YYckk4XFWm7369dPp8OlZUpYbdnOkQAyAnkVmuNi0Jakt6YOVMht1ufXnKJPr3kEv2jpSXTIQFATiIxw0rthYX645VXKmGzaXd1tT4vLtbAwECmwwKAnMNot7XGTSHdUlmp52+4IdNhAACAi3j+sssyHQIAAN9q3BTSAIDRY7QbAIDcQG61FguNAAAAAACWePvtt7Vy5UrV1tbKZrPppZdeGvb8z3/+c9lstmGP+fPnDztncHBQ9957r6qqqlRUVKRVq1bp9OnTY/guvolCGgBgypdd6XQ/AABA+mRLvu7v79fMmTO1YcOGC56zYsUKtbS0JB9/+ctfhj2/Zs0avfjii3r22We1c+dO9fX16eabb1Y8HjcdT7ow2g0ASBmj3QAA5IZs2WyssbFRjY2N33qOx+NRIBA473O9vb3auHGjnnrqKS1ZskSS9PTTT6u+vl7bt2/X8uXLTcWTLnSkAQAAAAAZ89Zbb6m6ulqXX365fvGLXygYDCaf279/v6LRqJYtW5Y8Vltbq2nTpmnXrl2ZCFcShTQAwASrxsRG8q15vq65AgAgHazM17FYTGfPnh32GBwcHFGcjY2NeuaZZ/Tmm2/qt7/9rfbu3atFixYl79fa2iq3262Kioph1/n9frW2to7uQxoFCmkAQE7K1zVXAABku6amJpWVlQ17PPzwwyO61+23367vf//7mjZtmlauXKnXXntNhw8f1quvvvqt1xmGIZvNNqLXTIe8WCP9+9//PtMhAMC4kS1rpPN1zVW+e+WVVzIdAgDkPas38ly4cKFeeOGFYcc8Hk9a7l1TU6OGhgYdOXJEkhQIBBSJRNTd3T2sKx0MBnX99den5TVHgo40ACBrpGtM7Eu5uOYKAIBs53Q6VVpaOuyRrkK6s7NTp06dUk1NjSRp9uzZcrlc2rZtW/KclpYWHTp0KKOFdF50pAEAY8PKXbtjsZjKysqGHf/Nb36jtWvXjuiejY2NuvXWW9XQ0KDjx4/rV7/6lRYtWqT9+/fL4/Fk7ZorAADSIVt27e7r69PRo0eTfx8/flzNzc2qrKxUZWWl1q5dqx/96EeqqanRiRMn9Mtf/lJVVVW65ZZbJEllZWW66667dP/998vn86myslIPPPCApk+fnpwoywQKaQCAKVYlZqfTqc7OzmHHRvPt9u23357892nTpmnOnDlqaGjQq6++qtWrV1/wukyvuQIAIB2y5acl9+3bp4ULFyb/vu+++yRJd955p/7whz/o4MGD+tOf/qSenh7V1NRo4cKFeu6551RSUpK85rHHHpPT6dRtt92mUCikxYsXa9OmTXI4HGP+fr5EIQ0AyBqlpaWW3TtX1lwBAJBPFixY8K1F/euvv37Re3i9Xq1fv17r169PZ2ijwhppAEDKsunnr8zKlTVXAACkQ67m61xBRxoAkJPydc0VAADIfhTSAABTsuXnr/J1zRUAAOmQLZuN5SsKaQBAyqzctdusfF1zBQBAOlDwWos10gAAAAAAmEBHGgBgSrZ0pAEAwIWRW61FRxoAAAAAABPoSAMAUpZNa6QBAMCFsdmYtSikAQCmUEgDAJD9yK3WYrQbAAAAAAAT6EgDAFJm1Wg3AABIL/K1tehIAwAAAABgAh1pAIAprJEGACD7sdmYtSikAQApY9duAAByA7nVWox2AwAAAABgAh1pAIApdKQBAMh+jHZbi440AAAAAAAm0JEGAKSMNdIAAOQGcqu1KKQBACmjkAYAIDeQW63FaDcAAAAAACbQkQYAmEJHGgCA7MdmY9aiIw0AAAAAgAl0pAEAKbNqjTQAAEgv8rW1KKQBAKYw2g0AQPYjt1qL0W4AAAAAAEygIw0ASBk/fwUAQG5gszFr0ZEGAAAAAMAEOtIAAFPoSAMAkP3IrdaikAYApIzRbgAAcgO51VqMdgMAAAAAYAIdaQCAKXSkAQDIfmw2Zi060gAAAAAAmEBHGgCQMtZIAwCQG8it1qKQBgCYQmIGACD7ka+txWg3AAAAAAAm0JEGAKSM0W4AAHIDm41Zi440AAAAAAAm0JEGAJhCRxoAgOxHbrUWHWkAQMq+HO224gEAANInW3L222+/rZUrV6q2tlY2m00vvfTSN+Jcu3atamtrVVBQoAULFujDDz8cds7g4KDuvfdeVVVVqaioSKtWrdLp06dH+xGNCoU0AAAAAMAS/f39mjlzpjZs2HDe5x955BE9+uij2rBhg/bu3atAIKClS5fq3LlzyXPWrFmjF198Uc8++6x27typvr4+3XzzzYrH42P1Nr6B0W4AgCmMdgMAkP2yJbc2NjaqsbHxvM8ZhqHHH39cDz30kFavXi1JevLJJ+X3+7V582bdfffd6u3t1caNG/XUU09pyZIlkqSnn35a9fX12r59u5YvXz5m7+Wr6EgDAAAAAMbc8ePH1draqmXLliWPeTwe3XTTTdq1a5ckaf/+/YpGo8POqa2t1bRp05LnZAIdaQBAyvj5KwAAcoOVuTUWi+ns2bPDjnk8Hnk8HlP3aW1tlST5/f5hx/1+v06ePJk8x+12q6Ki4hvnfHl9JtCRBgCkLFs2LgEAABdmZb42DENNTU0qKysb9nj44YdHHK/NZvtG/F8/dr73eLFzrERHGgAAAACQsoULF+qFF14YdsxsN1qSAoGApC+6zjU1NcnjwWAw2aUOBAKKRCLq7u4e1pUOBoO6/vrrRxJ+WtCRBgCYQjcaAIDsZ+UEmdPpVGlp6bDHSArpSZMmKRAIaNu2bcljkUhEO3bsSBbJs2fPlsvlGnZOS0uLDh06lNFCmo40AAAAAMASfX19Onr0aPLv48ePq7m5WZWVlZo4caLWrFmjdevWacqUKZoyZYrWrVunwsJC3XHHHZKksrIy3XXXXbr//vvl8/lUWVmpBx54QNOnT0/u4p0JdKQBACnLpjXSb7/9tlauXKna2lrZbDa99NJL34h17dq1qq2tVUFBgRYsWKAPP/xw2DmDg4O69957VVVVpaKiIq1atUqnT58ezUcEAEBWyJZ8vW/fPs2aNUuzZs2SJN13332aNWuWfv3rX0uSHnzwQa1Zs0b33HOP5syZozNnzmjr1q0qKSlJ3uOxxx7TD3/4Q91222264YYbVFhYqC1btsjhcKTnwxoBCmkAgCnZkpj7+/s1c+ZMbdiw4bzPP/LII3r00Ue1YcMG7d27V4FAQEuXLtW5c+eS56xZs0Yvvviinn32We3cuVN9fX26+eabFY/HR/z5AACQDbLly+8FCxac9x6bNm2S9MVGY2vXrlVLS4vC4bB27NihadOmDbuH1+vV+vXr1dnZqYGBAW3ZskX19fXp+qhGhNFuAEBOamxsVGNj43mfMwxDjz/+uB566CGtXr1akvTkk0/K7/dr8+bNuvvuu9Xb26uNGzfqqaeeSo6GPf3006qvr9f27du1fPnyMXsvAAAgt9CRBgCkLFu+3b6Y48ePq7W1VcuWLUse83g8uummm7Rr1y5J0v79+xWNRoedU1tbq2nTpiXPAQAgV+VCvs5ldKQBAFnj7Nmzw/72eDwj2gW0tbVVkpI/nfElv9+vkydPJs9xu93Dfkrjy3O+vB4AAOB86EgDAEyx6hvuWCymsrKyYY+HH354VLHabLZvxP71Y+d7fxc7BwCAbEdH2lp0pAEAKbMqiRqGIafTqc7OzmHHR9KNlqRAICDpi65zTU1N8ngwGEx2qQOBgCKRiLq7u4d1pYPBYEZ/lxIAgHSwquilmP4CHWkAQNYoLS0d9hhpIT1p0iQFAgFt27YteSwSiWjHjh3JInn27NlyuVzDzmlpadGhQ4copAEAwLeiIw0AMMWqjrRZfX19Onr0aPLv48ePq7m5WZWVlZo4caLWrFmjdevWacqUKZoyZYrWrVunwsJC3XHHHZKksrIy3XXXXbr//vvl8/lUWVmpBx54QNOnT0/u4g0AQK6ic2wtCmkAQE7at2+fFi5cmPz7vvvukyTdeeed2rRpkx588EGFQiHdc8896u7u1rx587R161aVlJQkr3nsscfkdDp12223KRQKafHixdq0aZMcDseYvx8AAJA7KKQBACnLpo1GFixY8K2x2Gw2rV27VmvXrr3gOV6vV+vXr9f69estiBAAgMzJlnydryikAQCmZMtoNwAAuDA2G7MWm40BAAAAAGACHWkAQMqs/PkrAACQPuRWa9GRBgAAAADABDrSAABT6EgDAJD9yK3WopAGAKSM0W4AAHIDm41Zi9FuAAAAAABMoCMNADCFjjQAANmP3GotOtIAAAAAAJhARxoAkDLWSAMAkBtYI20tCmkAQMqsKqQBAEB6ka+txWg3AAAAAAAm0JEGAJjCaDcAANmP3GotOtIAAAAAAJhARxoAkDI2GwMAIDew2Zi1KKQBAKZQSAMAkP3IrdZitBsAAAAAABPoSAMAUmblaLfNZkv7fQEAGK/oSFuLjjQAAAAAACbQkQYAmEJHGgCA7MdmY9aikAYApIxduwEAyA3kVmsx2g0AAAAAgAl0pAEApvANNwAA2Y98bS060gAAAAAAmEBHGgCQMtZIAwCQG9hszFoU0gAAUyikAQDIfuRWazHaDQAAAACACXSkAQApY7QbAIDcQG61Fh1pAAAAAABMoCMNADCFjjQAANmPzcasRSENAEgZSRkAgNxAbrUWo90AAAAAAJhARxoAYArfcAMAkP3I19aiIw0ASMnmzZstvX84HLb0/gAAjAcOhyMvXiPbUUgDAFLS1NSU6RAAAMBFTJo0yfLX+PGPf2z5a2Q7m0HPHwCQIpvNZtm929vbVVVVZdn9AQAYL/bu3au5c+dadn9KSAppAIAJH330ka6++mpL7k06AgAgfaz68ru3t1elpaWW3DuXMNoNAEjZVVddZcl9BwYGLLkvAADj1bFjx9J+zw0bNlBED6EjDQAwpaWlRbW1tWm735NPPqmf/exnabsfAAD4Qrq70pFIRC6XK633zFUU0gAA09KZmOPxuOx2BqQAAEi3rq4u+Xy+tNzrhRde0C233JKWe+UDCmkAgGl9fX0qKSkZ9X22bt2qpUuXpiEiAABwPun68juRSFi66WiuoQUAADCtuLhYTzzxxKjvQxENAIC1QqHQqO+xa9cuiuivoSMNABiRWCw2qnVSBw4c0IwZM9IYEQAAOJ/NmzfrJz/5yYivp2T8JgppAMCIbdmyRatWrRrRtaQfAADGRiKRkMPhGNG1R44c0eTJk9McUe6jkAYAjJhhGCPaKOzUqVOqq6uzICIAAHA+TU1NWrRokenrKBfPj0IaADAqe/fu1dy5c01dQ+oBAGDsmV3n3N7erqqqKouiyW0U0gCAUTOTmHt7e1VaWmphNAAA4Hw++ugjXX311SmfT6l4YezaDQAYtWPHjqV03oYNGyiiAQDIkKuuuirlcwcGBiyMJPfRkQYApEUqXelIJDKqnb4BAMDotLS0qLa29lvP2bRpk+68884xiig3UUgDANKiq6tLPp/vgs8///zzWr169RhGBAAAzudiX37HYrER7/I9XlBIAwDS5tsScyKRML3JCQAASL++vj6VlJSc97mtW7dq6dKlYxxR7mGNNAAgbUKh0HmP79q1iyIaAIAsUVxcrCeeeOK8z1FEp4aONAAgrZ555hn99Kc/HXaMVAMAQHaJxWLf2LfkwIEDmjFjRoYiyi0U0gCAtEokEsPWVR0+fFhTpkzJYEQAAOB8Xn75Zf3gBz9I/k1pmDoKaQBA2jU1NWnRokWSSMoAAGQrwzBkt3+x2vfUqVOqq6vLcES5g0IaAGCJHTt2aObMmSovL890KAAA4AJOnz6tYDCoa6+9NtOh5BQKaQAAAAAATGDXbgAAAAAATKCQBgAAAADABAppAAAAAABMoJAGAAAAAMAECmkAAAAAAEygkAYAAAAAwAQKaQAAAAAATKCQBgAAAADAhP8Hnaso3MpB5CMAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ds = ds.where(ds != 65535) # remove no data values\n", "arguments[\"bands\"] = [band.lower() for band in arguments[\"bands\"]]\n", "visualize_sat_images(ds, gdf, list(reversed(arguments[\"bands\"])))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### special arguments for asf" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# downloaded tile tifs are removed by default when creating a minicube, this behavior can be change with rm_tmp_files=False\n", "ds = tg.create(**arguments, rm_tmp_files=False)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# resampling: using a different resampling method in the reprojection process - default: rasterio.enums.Resampling = rasterio.enums.Resampling.nearest\n", "ds = tg.create(**arguments, resampling=rasterio.enums.Resampling.bilinear)" ] } ], "metadata": { "kernelspec": { "display_name": "tg", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.13.0" } }, "nbformat": 4, "nbformat_minor": 2 }