From d91ee32f445b7430cac14086814aefce1977ca50 Mon Sep 17 00:00:00 2001 From: davemlz Date: Sun, 21 Jan 2024 17:07:26 +0100 Subject: [PATCH] DOCS: Added stackstac tutorial --- docs/changelog.rst | 2 + docs/tutorials.rst | 3 +- docs/tutorials/stackstac.ipynb | 1383 ++++++++++++++++++++++++++++++++ 3 files changed, 1387 insertions(+), 1 deletion(-) create mode 100644 docs/tutorials/stackstac.ipynb diff --git a/docs/changelog.rst b/docs/changelog.rst index cf13f37..4ebbe9d 100644 --- a/docs/changelog.rst +++ b/docs/changelog.rst @@ -6,6 +6,8 @@ v2024.1.0 - Added support for Google Earth Engine. - Added the :code:`gee` argument to :code:`cubo.create()`. +- Added support for :code:`stackstac` keyword arguments. +- Added the :code:`stackstac_kw` argument to :code:`cubo.create()`. v2023.12.0 --------- diff --git a/docs/tutorials.rst b/docs/tutorials.rst index 54c9940..582804f 100644 --- a/docs/tutorials.rst +++ b/docs/tutorials.rst @@ -8,4 +8,5 @@ Tutorials tutorials/cube_visualization.ipynb tutorials/using_collections.ipynb tutorials/visualization_lexcube.ipynb - tutorials/using_gee.ipynb \ No newline at end of file + tutorials/using_gee.ipynb + tutorials/stackstac.ipynb \ No newline at end of file diff --git a/docs/tutorials/stackstac.ipynb b/docs/tutorials/stackstac.ipynb new file mode 100644 index 0000000..bf8388c --- /dev/null +++ b/docs/tutorials/stackstac.ipynb @@ -0,0 +1,1383 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Add `stackstac` keyword arguments" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "In this tutorial we will learn how to pass `stackstac` arguments to `cubo`." + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "import cubo" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "The following chunk doesn's pass any argument to `stackstac`, creating the cube using the default values:" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/home/dmontero/anaconda3/envs/cubo/lib/python3.9/site-packages/stackstac/prepare.py:364: UserWarning: The argument 'infer_datetime_format' is deprecated and will be removed in a future version. A strict version of it is now the default, see https://pandas.pydata.org/pdeps/0004-consistent-to-datetime-parsing.html. You can safely remove this argument.\n", + " times = pd.to_datetime(\n" + ] + }, + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray 'sentinel-2-l2a' (time: 2, band: 3, y: 64, x: 64)>\n",
+       "dask.array<fetch_raster_window, shape=(2, 3, 64, 64), dtype=float64, chunksize=(1, 1, 64, 64), chunktype=numpy.ndarray>\n",
+       "Coordinates: (12/47)\n",
+       "  * time                                     (time) datetime64[ns] 2021-06-13...\n",
+       "    id                                       (time) <U54 'S2B_MSIL2A_20210613...\n",
+       "  * band                                     (band) <U3 'B02' 'B03' 'B04'\n",
+       "  * x                                        (x) float64 6.011e+05 ... 6.023e+05\n",
+       "  * y                                        (y) float64 5.66e+06 ... 5.659e+06\n",
+       "    sat:orbit_state                          <U10 'descending'\n",
+       "    ...                                       ...\n",
+       "    title                                    (band) <U20 'Band 2 - Blue - 10m...\n",
+       "    common_name                              (band) <U5 'blue' 'green' 'red'\n",
+       "    center_wavelength                        (band) float64 0.49 0.56 0.665\n",
+       "    full_width_half_max                      (band) float64 0.098 0.045 0.038\n",
+       "    epsg                                     int64 32632\n",
+       "    cubo:distance_from_center                (y, x) float64 908.2 ... 873.7\n",
+       "Attributes:\n",
+       "    collection:           sentinel-2-l2a\n",
+       "    stac:                 https://planetarycomputer.microsoft.com/api/stac/v1\n",
+       "    epsg:                 32632\n",
+       "    resolution:           20\n",
+       "    edge_size:            64\n",
+       "    central_lat:          51.079225\n",
+       "    central_lon:          10.452173\n",
+       "    central_y:            5659638.0946523\n",
+       "    central_x:            601722.4825156148\n",
+       "    time_coverage_start:  2021-06-01\n",
+       "    time_coverage_end:    2021-07-01
" + ], + "text/plain": [ + "\n", + "dask.array\n", + "Coordinates: (12/47)\n", + " * time (time) datetime64[ns] 2021-06-13...\n", + " id (time) " + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEhCAYAAACDefxEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAAsTAAALEwEAmpwYAABX1klEQVR4nO29fbwtZ1nf/b1mZq21X89bAklIggk26hP7IFAQ1NpGoAjUSqs+fmhRgtqmUqHUVi1I+4i1WEVrW1uVJ4VosSpFRU19oAhorK2CAUU0vEbekpiQnJz3/bJeZq7+MXNgnX3/7n3mJOecfc4+1/fzOZ+z9qxZM/fcM7Nm3b/7d12XuTtBEARB0IdipxsQBEEQXDzEQyMIgiDoTTw0giAIgt7EQyMIgiDoTTw0giAIgt7EQyMIgiDoTTw0giAIgt7EQ2MOM3u8mZ0ws3Kn2xIElzpmdpeZ3ZR57yYzu/f8tiiAS/yhYWafMrNnn/zb3T/j7ivuXu9ku3KY2cjM3mhmnzaz42b2ATN73pZ1nmVmHzGzdTP7HTP7grn3ftzMPt599iNm9uItn73VzD5qZo2ZvaRHe55kZu/v9vV+M3vSlvefYmb/s3sQf9bMXrHNtr6ma+9RM/tUZl+/171/r5n9y275i7rtnzCzja7tJ/8+0a3zMjN7n5mNzeznxLazfSbWfZ2Z3WNmx7rz8P1n0idb1v1eM/uz7nx80sy+d8v7P2Rmf2pmMzN7TW473bpmZj9qZg93/37UzKx774vM7DfM7CEzO2Rm7zCzL85s591m5mZWzS27ruuX9a6fnr3lM99tZg90fXKbmY36fnY73P1L3f2OvuufDczsO7p2Hu+u2beZ2Wr33rbnq1vHzOwTZvahM9jnHV2ff9mW5b/WLb+p+3tf178PdG34mJm9cm59N7O/9EiPvTfufsn+Az4FPHun23EG7V0GXgNcR/vA/zrgOHBd9/7lwFHg/wEWgB8D3jP3+R8EvqT77NOBw8BXzr3/XcCzgPcBLzlNW4bAp4HvBkbAP+7+Hs615UHgRd37q8D/tc32vhz4VuAW4FPi/Q8BrwVK4AuB+4Gv37LOTcC94rPfAPxt4GeAn9vy3rZ9Jrb1xcBy9/pq4C7gG/r0idjW9wFPAapuu58GXjj3/s3A84DfAF5zmvPxD4GPAtd07foQ8J1zffsdwAFgAPwQ8BGxjRcB/xNwoJpb/gfATwCLwDcCR4DHdO99LfBZ4EuB/cAdwI/0+eyjvBfkuX6U2/zr3bE8ufv7QHcOVvucr7ltnAA2gaf13O8d3bn7t3PLLgMeoL2HbuqW/Szwlq6fC9p7+ZvmPuPAXzqbfSLbe653cKH+A34eaICN7iR/H+2X8edumO5k/mvg97t1/nt3Mn8BOAbcSfeF3a3/JcA7gUPdRfDN5+E4Pgh8Y/f6FuD3595b7o7vSzKfvR34Z2L5/+L0D43nAPcBNrfsM8Bzu9c/DPz8IzieZ6MfGuvAjXN//zLwqi3rbPtF0p3Ln9uy7Iz6bMtnrwb+FPi+Pn3SY3s/CfxHsfy/cvqHxu8Dt8z9/R1kHn60X4YOXDa3bC/wMeAZW+6BLwLGdF+c3bLf4/MPpF8EfnjuvWcBD/T5bI/++BTdjzrah87P0f7Q+RDwvdud60fyD/ge4NfPYP3kfAG30X4/vBX4Tz23cwfw/wL3AmW37GW0P3Lu5fMPjT8D/vY22zkvD41LVp5y92+lvaH/lreS1Osyq76Q9hfw1bS/cP+A9ol/APgw8AMAZrZM+8D4ReCx3ed+2sxuVBs1s582syOZfx/scwxmdgXtjXlXt+hLgT+ZO8Y14M+75Vs/uwg8be6zZ8qXAh/07mrt+ODcvp4BHDKz3zezB83sv5vZ4x/hvgD+PfBiMxt00spXAO96FNs7ybZ9ZmavNLPfnP9At+wE7Q29THvOT24r2ydm9vdy57aTkr6aR3c+/mTu7z9BnPeOv0b7xf7w3LIfpv2SekBs9xPufjyzbbXfK8zssh6fPRN+gPb++0La0c3N261sZh/c5v766czH3gt8rZn9oJl91bzMJrafnC8zWwK+ifah8QvAC81s2PP4/oL2Yfic7u8XA2/ass57gNea2beZ2Q09t3vWuWQfGmfAz7r7n7v7UeDtwJ+7+7vcfUb7a/fJ3XpfR/sL+Wfdfebufwz8Kq3skeDu/8jd92X+PfF0jTKzAe2F+V/c/SPd4hVaqWWeo7TS0FZeT3sDv+N0+8pwun1dQ3tjvwJ4PPBJ4Jce4b4AfpP2htwAPgK80d3vfBTbO8m2x+HuP+LuXzf/prv/SPf+U2hHrCc/f7pt/eI25/Y1tPfjzz6io0j3fRRYOTmvcRIzuwb4KeCfzi17KvBVwH/ssd2T217NvH/y9WqPz54J3wy81t0Pufs9tL/ys7j7E7e5v/5R5jO/RytlPgX4/4GHzewnTBtjXkN6vr6BdmT1W93nB8DfPINjfBPtD6MvAfa5+x9sef/ltPf8y4APmdndtmVO83wQD43T89m51xvi75Xu9RcAT5//RUOrEV95thtkZgXtl9WE9gI6yQlgz5bV99DOe8x//seAv0wrn/VKc2xzk8vdiOF0+9oAfs3d73T3Tdr5lK80s71m9v1z23p9j30fAP4H8K9o5x2upf1FKG/+M6RXn23FW/6Y9jh/8NFsy8xeRvvL8m+6+7hnu7eydd97gBPz59fMHkP7hfbT7v5L3bIC+GngFd0PodNt9+S2j2feP/n6eI/PngmPA+6Z+/vTj2Abp8Xd3+7uf4tWSXgB8BLg78+vs835uhl4S/ejcZP2R+O2I6ItvBV4Ju09/fOibRvu/sPu/ldoZfK3AL/c3R/njUv9oXE288LfA/zull80K+7+UrWymb1+yxfx/L+sRNH9cnwjcAXtXMZ07u27gC+bW3eZdjg/P4T+QdrJ1ee4+7G+B9cdy8l/n+m2+cQtv2SfOLevD3Jq/37udXfhn9zWd/bY/ROA2t3f1N2Q9wJvBp7ft/3bcNo+Ow1Vt/7JbW3XJwlm9u3AK4Fndcf1SDnlOLrX8+d9P+0D43Z3f+3cenuApwL/zcweoJ2nA7jXzE7KL0846SAS21b7/WwnfZ3us2fC/bQ/Fk6yrdRprV03d3+d9oeKuzfu/m7gt2l/YJ3crjxf3QjumcC3dO6mB2hHxs83s8v7HKC7r9OqGS9FPDS2rHuMVlJcBq7vs/2zxrmeNLmQ/9FqhPOTh9eRToT//bn3T5lIpZ20vbt7vUr76+dbaYelA9o5g6xj6BG2+fVdu1fEe4+hHf5/I+0v8h/lVPfUq4CPA1dmtj3sPve/gX/QvS62WffTtPLTiPbX0bx76pm0k5ZP6vri3wG/t81xFd3+ntdtZ2FuW3toXTd/r1vvStq5pR/eso2b0O6pqtvev6G9GRfmzvG2fSba+A9p3StG60q6H/jHffpEbO9FtHMI8hrp+m2Bds7kX3evy8y630k7x3Y17a/yu/j8ZPUe4A8RE7PdcVw59+9ptPfA1XP9/x7gx7v9/x1OdU89tzuGG4F9tF+y8+6p7T57E+2gLXdNfIrPT4T/KPC7Xd9fQ/uj5GxPhL+Adi5y/vw+BLzodOeL9t768Ja+vBL4BPDy0+z3Drrvme7c/dW59+Ynwv9ld35O3qevpr3HVrr3wz11zg++vUg+013I38OjeGh0f38xrZb5EPBwdwM96Sy29wu69m3SDv1P/nvRljZ9hFY2uYNT3V1Oq7nOf/b7t1y8vuXfTdu058nA+7t9/RGdVXHu/ZfSuokO0zrPrt1mWzeJfd8x9/4zaX8FH+1u3P8MLIltqIfGa8S2X9Ozz74feHv3uqCVyQ51ffex7n3r0ye0Xzp3zf39SWC65Xy8fu79nxPtfkn33lfTyk8n1zXgdV3bDnWvrXvv5u6za1v29XjRV9eRWm6v6/plg9YV+Owtn/mntLLtMVqNf9Tns7Q/sP73NtfEp/j8Q2OJVvM/wrlzT/014N3AQVoJ7WN0zrjTna/u+kkeDrSuzPedZr93MPc9s+W9+YfGv6B1UB3rzvEdnGqZ/9xDQ1xrr99ybd3F3PfGmfw7eVEFQRCcV8zsDcAvu/sjNWMEO0A8NIIgCILeVKdfJQiCIHg0dHE9iud5a/W9aIiRRhAEQdCbS91yGwRBEJwBl/RDw8xeY2b3WZst9gNmJn3/XXbJX7E2++WHzewr5t57ebf8LjPLpSKZ39bjzey3uu18yMyuO4uHFARBcE65ZOY0rE0v/BJ3f8mWt/6du//4aT7+H4D/4e7f1OWSWeq2+TW0tt0vc/exmT22R1PeRJsO4Z1mtkKbNDEIguCi4JIeafTBzPbS+rffCODuE3c/0r39UtpApnH33oPdZ0oz+zEzu7NLnPYPu+U30vrf39mtf8LbKNAgCIKLgnhowMu6L/bbulQLW7meNljvZ83sj83sDV2qCWgzzH61mb3XzH7XzJ7WLf8O4Ki7P402gvMfmNn13fpHzOyt3bZ+LJMMLQiC4IJk1z80ui/0DwBvAL5+bv7ia2lTQX8hbaqL+4F/KzZR0Wa9/Bl3fzJtVO0r5947QJsG/HuBt3R5h55Dm63yA7Tpli8DbujW/2ra6POn0eZUesnZPeIgCIJzx66f03D3p8O2cxp07/9n2vTbW7mXNl3Be7u/f4XPPzTuBd7qrW/5D82soa0EZ7QpBU6JdDWzZwAfcPdPdH//Ou0D542P8PCCIAjOK7t+pLEdZnbV3J9/hzavyym4+wPAPfb5msrPos19A/DrwNd02/oi2kRiB2lrVLzU2poXJ2s0L9PmTtrXpaiGNp9S71rCQRAEO82uH2mchteZ2ZNoE319ijaDKWb2OOAN7n7Sgvty4Bc659QngG/rlt8G3GZmf0Zb2+Jmd/cup851wB91ctVDtGUaazP7HuDd3fL30ybeC4IguCiIiPAgCIKgN5e0PBUEQRCcGbtaniqqwsvRFkerGliZWAY67E49ZuV6mY02YmW1rhwB5hqq6DuCzGxTLZabfJQjVRM7UvsW3eaN3rf1/Cmkdp3vj3S53r9YJj5rmevDa/X5R7ybdtVzISY8ija5n0FD+16H6pxnjlteH+r60h+XKMXGVOPlMWauBblNQSkO6FGe9NnG7KC7P0a9t6sfGuWo5LIbTy2fa+O025uRPmm2LkomL4h1VVXnZd219fpmsqxYGqX7nk7TZbX+NmzEF1Djdfp5cSEVnrkERuld5HW/b/NC3IC1elgCxUIapuIqckWct5k6P8BgUXzBi2OvRml/eq3DZky0szme9vGMdFkprq9icSj3Ux9Nj6lQq87EF0rmu7gWHWqinSYeWPKUA+Ug3eZsmp7joejj2Wbmqd6IY1cPV/GwbhbENmfpMQJUW39IAj5J11PXjN4iNNP0nUJ8xfss7aNyoPtjNlb9Ic7lXnGBbGRaWokfMOKoPvvHB7M12EOeCoIgCHqzq0caePrjxSoxbJtmfk6JXy+NeswW6ZO6cv3L2sv0V4GLUYFPxC/Jgdwk6td+tSdduTme/nJpMr/GfEP8OlYD9qHqo3SZkY6c2o2KX30zMYIQx75weTpCA5ieSPdViF/GXojRR+antZKipuKXeVGlx16Ka8HXxU9boFpNb0nVHUzEr1DxCxqg2BDnvUn30wzS9SwzfGnEsVej9CSNj6XHWS7JTVIOxL5KcUzi3vBGjJwyad2mk3Q/w2HaH+PjaduHQ90ftTjvjYn7RclGmUGBizapn/n1Wnq9Vyt6POAz0X5xD25HjDSCIAiC3sRDIwiCIOhNPDSCIAiC3uzqOQ0rjcHqqbr3bCO1Onl5BtbNtXRRKTTBei03pyEcO1OhjYvnuSuNFCiW0/0XwqXRlMremjl2ofvKQNA63c90M9WCB0sZV5Jop4l5luma0OUX9OVbir5rmrRNzThdb7Cot3nsvuPpfoRDTu27Fq63TLejTJV+SMzR7BeOuTUtjps4Ryz1c+ZQ6fNWVuk21TU/PCA+L+bLAIpKzOccFtZE5UBaTudTCteTgGpusBCWv1LMUyyM9DzaxuH0+pqNRB+L/iyGuo+HY+HeEvNwg9X0gFydc2A2TftzYWlBrpsjRhpBEARBb+KhEQRBEPRmV8tT3sBsa5CLcH6WmYKrNhK2xELIW0Mx3m0yNt6psFqKgDIT+kXTaNtqc0wcwELaJhPWYBm1irajFqJNk/W0TUq2yf4+KUUw3FRIDQupVdn2ajlmckjIJINUVvBhejzHPntMbnMggvEGK0JOGQt7rZKccsH9E3GOF4T9WmygEXIKAEJ+rQohaQg7urgM2nU3hXQyE20/JD4sbNoAtZB+axPnWLSprFMpSEbXA4Ww0tefSYNuXZzfo2I/AM1mWoCzKMR9IAI9JxlZsRSBeCib92baJmnTBirS75rJ8YwdPkOMNIIgCILe7NhDw8yuNbPfMbMPmdldZvaKbvkBM3unmX28+39/t9zM7CfN7O6uPOtTdqrtQRAElyo7OdKYAf/M3W+krV73XWZ2I21VvHe7+w3Au/l8lbzn0ZZMvQG4hbZUaxAEQXAe2bE5DXe/n7YuN+5+3Mw+DFwNvAC4qVvtvwB3AP+8W/6mrrTqe8xsn5ld1W1HYgbVFj1X5Hmjzsw/FCoJWZVqwbWyt2U8lcMD6ednqRwq9dhqqK1xNhKatzjQ2aFU6FwQWj/AsUOpxlstKxtw+tlKpBHJZHTAhR11NhaJABeFlvsXuo8HV6bH1IjUEyrJ3PCA7uNKzRGJJIajcaoPNyohXc4+LdKV1GLKzNXPvYyOXYg5kVpNDIh5CgZab2/E3E2lUuSIy8sHGZu3EOJL0e/KlmwiWaJMfAnMhH3bD4ivQnHoRSaNCMK2qqzOLmy0Bxb1NbdRpvOnk810/+LQmSl7PVBsCov7yplkz75A5jTM7DrgycB7gSvmHgQPAFd0r68G7pn72L3dsiAIguA8seMPDTNbAX4V+Cfufop1pRtVnFFieDO7xczeZ2bva9QjOAiCIHjE7Kjl1swGtA+MX3D3t3aLP3tSdjKzq4AHu+X3AdfOffyabtkpuPutwK0Ag8WBT45tkQuE/bAQ0gdAI2o4qEycQ2HhU7UFQDtxCyEvNaL+RHMk8xBcSTWAUmTOVTUHToiMsACLj0m3ORN2vUqM4V0850sRtQ4wFZlER48V8tIREdGdy/qr6kJsiDatijZlImmnh1NJY/PYiWTZ8LL0lhqMhbSVsZ1ORVR2sSKyuh5OL6RqlJGSZsI6Lqysyt2qrNsApk6n6DsTUftFRvKqhaxZiMtT3q3ivlQR3QAz0c5K9dFQXNuZ68NERLkJWdKFZffQMeE1BvYspLJVLWrFzE6kUrKyZAOY+PyozKbPluyke8qANwIfdvefmHvrduDm7vXNwG/MLX9x56J6BnB0u/mMIAiC4OyzkyONrwK+FfhTM/tAt+z7gR8B3mJm3wF8Gvjm7r23Ac8H7gbWgW87r60NgiAIMJmEbpcwWBn65U/cUuZWDDlnhzNuDmUmEYn3ahHF22QKm5gq/SkcGc0xUdhJFZACEAV1TBSlKZfS9Q7fl4mAHqYSUVGKIj1iaOuyqoz+fVIM0mH9RJRRVTWtVBlV0GVLVanacq9IJJhJXlmIEqW1kBpsSUXii2NvVI1gYCyuG3GPupByZAQxoOb2bEW080QqAapSsQCFuA9EnSs2j6XH2Wg1RkonlZCdlCxZTESSyqxlT3xe9GetwqpVLWOgaVLnmDBaam0tZ/PytO/Gx4VsLUpDL16/qDcp7gMX3zUP3fXw+939qWobOz4RHgRBEFw8xEMjCIIg6E08NIIgCILe7Oost7jD5FRd0kwVa08tawDFINUpZ5sqTFRF/GbCczfTbTYTsU0RuUmmkFFdC/ukiLodP5RqpIOBttsV4vMmwntNFeMRUa9UWsOfPZz2XbVH2E5FhKvnrNInhIYvIs+bdWHDHek+Hgt9uyhEZPKamH8Yboj2ZOZjxOTNHjE/tmHi82J+CABR5Mu3Zn8mI7dn5tFMzOc0tbAbCzt4IyLU25XFMnUtqcS3KoOzOOeATGs9E/MC1aLIjJyxsqpE11OR2cBERmuG2vY+PiHmzGaiMNQVS8myRn+lyazWqoDVdsRIIwiCIOhNPDSCIAiC3uxyeQrYMpwzYcMrR5mo16l4pi4ISWJdSC8imR3ATEhR1VQV1BGJETMuTRXl7if6yWhVpga0LBxUCMlNtL1aSiU4lfwRoFkVVkcRzS4jtY9mEq2pxUKSKPcLuU0kNgQwIZNUAxHtLLqzEVHVhTo/gInzfuJYqjU0q+lni83c7Sxs3krjEWHeKsEeQKOuORGRXgo5pVTeXGB6NL3ATek+Lq5DtZqwnQO4CmcX14yKRjBh2QddQKpQduFhamt2ZbMGBoup7FTsEfuepdtUkfQAM3FQw8z5yBEjjSAIgqA38dAIgiAIehMPjSAIgqA3u3tOgzSlRd0IjTVjf1TFlUqhjU9LkZ1TaLkAw8vFXMVxYa0T8xemiuQARTqFwKZIybAojn2z0dbgQlgLlduwfIxII6JSHZjuDxeTAEphna6Jgkn79FxUvSbsk6XS69Njnx3L2FZFeuJmJvpoJuy+4tCnIu0FAJNUn1ZtN1EESaXhAJJ5vbYB4hxV4joc6XaqtBuNuEBqlcpDFB0DMLGqq4y4o/Rry9eFiK+s8MBAzAc1G+m1VIniROP79MSiLYsCVCr3jfr4cqZgm2j/VF0fwhrciOsQoBLpcFxk2d6OGGkEQRAEvYmHRhAEQdCb3S1PGanWIUbwtciKCoCQnURtJExVVtLOPOqD6fjUZ0KQUdHGj9XP+MnBdMiKsI6OCxHlnSsqI2x4QxE9XosiTo1oe7EsdwObQiISUfe1sCqvH9RyzMLVop+ERON12m/DjBwzEcmATahjSqKR0deZyOJC1Ha3/SIb75qSQzJSo+oOFQE9FJlR1QUPNOKaN1Hf3IVFPWeFLfeoDhXXwtH0OBdFwaJNdFi0b4qMA2LXKtMCq/orc3w41dxUpDbiXpfZcIGxOE4Q95tQ5spMwdNir8hIcYYVTmOkEQRBEPQmHhpBEARBb3a3POXpMFrJLsKcAkA9FrKPSBRXiOI3tQxLBvO0y6uRcLKoYkAi4SBoV9PC5aJYlHCtKDcZ6ER1mwfTsFflihrsSd0c60dFXwKFcCAtCIeZqvdsmaJDvibcbCKxoo1SSWOykXHHiJreLopSVUKam55QyR912wcierwWrrlaZCGoltJ+B0DIkgPluFGJM5cyfXy037qliIpG1M5u9y+ueRduR6H6jNfEeSv115sqWjQbK0lV7FtFqAOjK0XtbuEiXBHfNbNCn7dySdwzQuc14XqbPaSlSnV9Zr6qssRIIwiCIOhNPDSCIAiC3sRDIwiCIOjN7p7TMKMYnnqIJgqjqIyuoDOOqiy5tYjSLHKZc0VhmMmaiDwX0agzocsDVEV6GhvhB51upMWAKqU5g5zoGV4pitKoolaij5cvF0I0UB9KbZH1IXGORHcORxm9XWQnVllqp2vivGXmt0pR5GemrNZrqZZc7hN+41l6LgBcWFxNTboJbXqYuY7HYk6kEsfjyqaZKeZTimj8WhR7mm2I48lEw7tIzVqoAlSH0vO2IDIr52zvrjI4LIj+UFkZMllu1TXPMO2PdZHRtjym5x9MFKuqJ+l1454ee/WYXCbhtO/qzNxgjhhpBEEQBL2Jh0YQBEHQm10uTzm2dRgurJOFSqoGNGK43UzSoaQ3qhhPRjoRtsShKIRUC/VitJjRTlRdazFcXhC1gF1IWyDr8cjiN8c+dTxZVomo2SoXobpHRN2vieSAIpvdNBOt7GIIviKS3G2K83utSMAIcO/DaX+qKOJiNY2qrkWtZxX5DdCsCClpOT0ZMxHhLi2rwHC/+PymiNSuhEQj5CWAej1drlQOF0V/GmEBBqhVQa+hsGTvTTt+cjDV0Rb26zQEjZCIZyIrhC+qUPrMd8XD6bLyMiFvi/t/NtEVkwbie8VEOs/B3nS92YN6m6X4XjrTkUOMNIIgCILeXHQPDTN7rpl91MzuNrNX7nR7giAILiUuqoeGmZXATwHPA24E/q6Z3bizrQqCILh0uNjmNL4cuNvdPwFgZm8GXgB8SK7twBZ5W1lmpyI1CEAltFN7UFgIl1MtuKn1/IMJ3bmZieI3IhtmnSkWPxS2RJXF1MV+ihW9zZnIzKqmPwaPEXY/kb6hyGSPZSZ0dGEnrZZEEaXcebtMrCvmohZECpSH1zJzLyNxjlXVIFEsqhZzL8VSmsIEwD2dj5kdU1ZUdR3JTTK5/4T4fHrezMT1IeypALVo50TMcyweSOd4mkwakSJtJo3Q4GtxD1X7RVqUzFzlRGRWVvMkLmzaY1H8CmDxqnT/tbgWKuE8t8x3xUxkEZkqa++xdNlIzccA64fSYx9lMvfmuKhGGsDVwD1zf9/bLQuCIAjOAxfbSOO0mNktwC0AhXBKBUEQBI+ci+2hcR9w7dzf13TLPoe73wrcCjBYHvjWyNla1IUerGayg3oqC5T7xLB+mm5TFnABKtHltahRroax9USfrkkjQldVYRXx8YkY2gKUwj85WxeR1kImqTdEVPRQ9/FMWAhN9DvroniNyMQLMPP0mCbK1izOW5VL+SlkklJlURVR4sNlsc2MlGTCx9sIu7Kvi2NfyhRhEpmMTdSvLhdE7W1ZsR0KIblVolhULSS82XGts45W04t+uikrDCWYyEwwE5kaAEwEj5cian4iLtnRAW3Jnt4vsgsIK+xYRLPbTDQIaEQBqqGS4aapP3/9uD5v1eVCOs7IvDkutp/idwI3mNn1ZjYEXgjcvsNtCoIguGS4qEYa7j4zs5cB76D9vXGbu9+1w80KgiC4ZLioHhoA7v424G29199St7hYTYeCqs41IIvFmBhuuxwu6yHfbFE4P5REoxxIM93ORhSBLvenxzk5lg5jLZP0r1b1nkWxl4FQwRqRzM5FckCASshbM+HsKUTxrPG6qCoFDEWt60bIKY2qeZ6r3a2KMG0IaW6rXQ8oxLn0jLMHIRtVKuK3EX000nKdKmDlquCSOPYmZ9kTuypE9PZMZA00mW4AZqKwVCFcQKooVSmcX02mulq1pPpTHJCI5FfOKwAX91F5Qn1XpPsuH5uJ5G+EXCcsVcOl9HqfDTLSnHBqbQhn4XZcbPJUEARBsIPEQyMIgiDoTTw0giAIgt5cdHMaZ4KZUW0pwqQy1xaZ7LG+T0SjHkw1XhM2S8tEZJYi66eJiO7x4VTfnSn7IVAJq2QtwoNVAahS6LsAzabQlxHWQLHIhY69KQpNAQz2pMeuMgQ347Q9g8VMVHUt9iXsk6qI0/So7mM1fzFYFRkDmrRN1STdz1jMbQEsTNMOVdH51qTtnIrswAD1grDXimytYzHPUYnswACVmGNS95aavzBRnAjARebfkZgwnIm5KBPtGWanjUR0v5rjKcU5zxRXq8W8QCUKqW0eS6/jqk7nJAB8cy1ZtriSznOo4mxVZt5oKuzOKkvFdsRIIwiCIOhNPDSCIAiC3uxqecpx6i2Fg2qR8Gup1DJHc0IknxO2QhdJ9xCR5wCNCEd1UXd48Jj01AxMt7PeSIeczVTIS02q0Qi3btsmMbxthB10JuzKw8eKY9yr92Mrwv4oiihVIqJb2kaBUlilm2MiaeBl6b6HmewAs01xLYhjL0RBrM1p+tlyoG+9saoxLpSGxRVRiEjILgAIaU1JooWpmtRaritEUslGZQIQFveRsAADTITUubGZXguDxXTFqaidnbNPjxZEkk1VYnwirOMqSSWZomVim6NBeg8X4r4CmIjzNh6nkpXKVjDM2I1tQRRhWslkxMgQI40gCIKgN/HQCIIgCHoTD40gCIKgN7t6TgMH35JBdmFfqinOhNUPYKqsm8oyKwrAI/Rh0MWAEJlVVfYGK3VqgHJZPPvXhKC6J11UqPkYYONQ2oDB3nSbynY6+YvNtI2L+lKrjwhrsJCi1efrTLbUmcgAWwib52xTFOQaZ/pY2JpFnSwqdc2IA2pUMR2gWBDnciP9/PGDqbadSwYxEm1HpY8RRZSKBa13NyKLs7KZz9ZFMa+MxR2ht5fCXjt+ME0f46LfhplrzlXqmonod2Fr9kwWZFX4qz4s+kjMiZKZ32I5XbRwVWrP3bwvnc/Jff8MxOLpRnq/bkeMNIIgCILexEMjCIIg6M2ulqcMSwqMTA+LIeey7gYVmWzrwoYnIjLJyA+NqEtdimJCxaqIhM3Ur54cFxKAkqKEDS/TTKqBiGZ9OB3GVmJoXQoZyzPe3tlYbFNEAY+PpkPwwT4dSVuKyPf6SPr5mZCihnty9kMhqYjTXlu6Tbf0mlkY6sjiphbZWoW0Vgvn5fAK3fbigKhfLTIeqzrsdaZGeLMuIqDVeRfR7JtH+kuAKp3u8Ko0KtrEfkTtKQBmShIVuxauZCxTDbQR0p66N2aiOFqxqeWhUviAx0qKWhfy9kqmQJn4rhhcpq/FHDHSCIIgCHoTD40gCIKgN/HQCIIgCHqzq+c0wGGLBdJEiglmacoNAC9Tra8RKSFUVTlf0s9jlVHShN5fj1Pt04RVEGAg7JON0IJdpEWZiKybAIt7RKqFKt3PYDldb/PhdJsmqv4BLIi0HeowZ2KeQ2UrBahPiDQkYt5qcZrue7PQ27QyPR8yWaunlm4T2VKLQqeOUHsvxPWxcEW6nzqXRkTMRS3tSftjQ1nPlYUYqA6IaoCi30uRVsWWM3MvYm6vHosUOSLVSi0q6g336XuwEFbcRmT9VXMs0zUxpwDg6XdFdbmY43kwPZ56Sc8p+DQ9pqJJ1y2XRTXQ/Zk+Fl7+yfGo3BcEQRCcI+KhEQRBEPRml8tTqW12pgI6M6Mzm4kIV1LpZSisueOjWiqohOwzE1HZpfALZtQpmlG6rioM5U3a9tWMbXXdhTVZRMOraPpyNd3eYCQWAo0YgtdiZF3WwlZ4QksnKqrahbw1EXJIJuCXRuyqXBR9fEJ9VhV70rJgdZU4eBHhrpaZuLYAphuppOIPCUl0lB7kQNhwAcZpUDZDkfl2cJnISHtMW25nU3EdC6tztSCkSlEAymf6600FS/tILBQ33CDTxxufEhmPB+lx1jORCVjIvgAIe6xqvJPew+NDGcndRaEtkUF6O2KkEQRBEPQmHhpBEARBb3a3PGUGW4d+ojBLIVxSAC6cBsNKOCqEdtFkwlHVMHphr3DxiFrVqggSgImkcC6G8Ii2j4XkBFCI0e1gSTid9qR9N10TMkNmBG61qistVhZtHxRa5qiFzFEPRCSucAYppxJAPRHuGpHkbizkh0Lpinsy9bzFuSxFFLKJ2uoz4VoDKGpxLRwQjpsNke1ARDAD+Eba9xtCspLVjTJFmFxEmTeL6TVnwlFV7hVfZaKIEeikgS6SQg5FHx/9ZMZt+AVplPr4nuPJsmo5lZJsrPtDyd6b4ne+cnRVq3qblSiKlTEMZomRRhAEQdCbHXlomNmPmdlHzOyDZvZrZrZv7r1XmdndZvZRM/vaueXP7ZbdbWav3Il2B0EQXOrs1EjjncBfdvcnAh8DXgVgZjcCLwS+FHgu8NNmVppZCfwU8DzgRuDvdusGQRAE55EdmdNw99+a+/M9wDd1r18AvNndx8Anzexu4Mu79+52908AmNmbu3U/tO2ODHzLnIYpF1vO8ia0xnpRZN0UmrOLrJcALiLKJwfFPEsa8Esz1u1sRHRxJaKQJyKl7SCTxVRFzq+LTLGNiMQVAfIsXKYvtenRtJ0q6n6wT/SRiEYHQEXdixNfmrDhZqJjy6W0/a4y/B4SFuIiXS/3a60UlspG+X2Pi2zJw0zRIXHem7GwSgvb6UzMOQGYCW18NW3nRGRwNX1r0Ij7bbSazp2Mj6V9PBK25tx+XMwVIBLNTjbTfY/EnADouZ/B49IqSiqa3TczmQ3U17OJ4mir4tqUW4R6IgpDZeY1c1wIcxrfDry9e301cM/ce/d2y3LLgyAIgvPIORtpmNm7gCvFW69299/o1nk1bWjdL5zF/d4C3AI6700QBEHwyDlnDw13f/Z275vZS4CvA57lnw81vg+4dm61a7plbLN8635vBW4FGKwOfOanDgeLSfogyVlZZUFdka9sKkpVF5laPtakw1u7LF1v9llhJ1XtARDRzoqFvaI++jFdZ1tFw6ro3AWR+G56IrUlNtOMlCTaXi0LqUFEGzcPZiKLhfxQi1QAviCsm56x3J4QfT9INQ0Xcp8fFUV/LssUpRKWbJUdQNWqbgotrVUr6cWoItJtJZUuVL8BVEsioaZIpudKKBnq+21UpJ83IecMD6TX8URcx1WVsdJvpm0aXZ320eSe1EMsa7gDhZBEJy6kKPFlYZkft4VIWjo4kB7T5IS4jpVGDJiI8J8evggSFprZc4HvA77e3efPzO3AC81sZGbXAzcAfwjcCdxgZtdbK6a+sFs3CIIgOI/sVHDffwJGwDu7QK73uPt3uvtdZvYW2gnuGfBd3iVLMbOXAe+grbt5m7vftTNND4IguHTZKffUX9rmvdcCrxXL3wa87Vy2KwiCINie3Z1GpAa2prQQR1xmMnnWaRYAGpG6ohKaZHE8k4F1KdU+pw+LeQWRBqTM2B/Zky5yUWx+KlJcVKu5Y0/XXXhcmiqhGYttDoQO3mQKSIlCOZsPpHMFxaoocpPRxquB0LJn6XHWG6KwksjKCjA7ks4BFCIL8kQU8zFxzstNfes1Q2GPFSJyIyyeVanP5WwqUl+I86Gcl6orAVzYlUtR6Gr9gbVk2fI16XUEMD6Y3gcjkVJmQVidbVlkpF7PZBIeihQs6yKrs8hibCobLjB+QGS5fbyYvxSZdwcirQlAPRXzr5WwjosUKoXItgwgEl1T7j2zWYoLwXIbBEEQXCTEQyMIgiDoza6WpwzYGhhdC3utrWcsq8IyVy6KKHEhcxSLepvKeVqKIWsjCqO4yO4JMDssMpaK2kqlKMBSH89cAkL5MYRdUEQ7U4rhcqY/GhGhOnqcaLyQQyx39arstSpK/DJhET2h7YeFkDBtkFo/q2OpLmALIjtw5pJrVCYBsWgo2q6uQ4DBZakcNJmmslEjMvyayBgMOgP09GjqRx+q2u6iVj3IRMbUwsK8KdLpTkTRs4WM7DwRGW1NSKoqynucid4u96TtHBwW2a/3pl8AY2GZBSiV7V1c9M1M3JdCKgQoxPfXcC2KMAVBEATniHhoBEEQBL3Z1fKUA/WWRHWFitgVjgSAmajDXCnnhqinTSYatTwuIkc3RN1eUVJ7OtMR0MNVEendiOH6Yrre+JDeZrEiZLgj6bouEgFWImp2lkkEaAg3iGjndD2VUzwjK1arqbxVHhBusmMimd5yxk0m3DGspG2yJXHehSQ6ySSJUwWoKuGOmYjEiCauTYBGSG4DESU+EFkMJhkpSRbaHoksmyKqeXxYyzGLjxXRykoOEokqlctKmNsAKJdExgBxfVqR7qdcySQ3FUW6JuLaLjaECy/jLPSpaKes8S2ujzWRgREoxe0+6ZlR4iQx0giCIAh6Ew+NIAiCoDfx0AiCIAh6s6vnNKw0ypVTdUlvRPF6HThKuU9ofWL+oRylOvZUWDwBTOiH5QERMSz00NH+tKgLQC3mVAoRTbp5JNU5m5nWrItG2ElH6eVSC911ILKVrq+lNkmAQvQdayJr51KqL6uCVgAudOPJhiqElPZ7XWvNenBtuv/muIgiFvp0uZjOsdSbuj9MWIux9LddIQp/qfkhgKGYI1q7N0134Csi6n6P1rvHD4piPkVquS1EZuXhsk4BPV4TN6IqYHVQWKr3i3Zmsle7CvovxbpqWkBPATIT100lCnehrtmVzNfwhrDcijkNE/seLurMBjYRc4ui4Np2xEgjCIIg6E08NIIgCILe7Gp5CidJzFZaOjT2JrVOAjARQzwVUS4SxYkS3e26IzG+FQ7ETVHs6dTSI5+nWhCnsRLDUCFZLR3Q8oOsZS6SEw7EMHi8lvZnISQWgFIkHWxEfwwmwsa7J1dkR0TICgtyKaQPq/SJq1UUsUoQuF9EFh9K+6PMyGCFSB43O5jqJOXlaZS3N9pmWd8vZKMlcR+Ii3b8Wa3dCoerDHM3kQjQRMYA0GrScDFdOBWR2hsiQedAJPIEEPW8sLEoaiXky0IkEgUYig5phFxYi9uyOJSRWcUyE+eoHqTfKaXre2OqEqmOMkXoMsRIIwiCIOhNPDSCIAiC3sRDIwiCIOjNrp7T8AaatVM1wFr56FSmVgCRybNxYX8UOuVslik6VKZzANNNoUkuif1k5gUQxZlqVeRnkh573Wj74/pnxDyPsnmKw5wJ6+TosStyP7UoDDUSNs36qMg4rARiYLKWasS+LDRn0fjSM9fCQNiaG6FvC23chKXSa30uZ4fEfkTSX1eZUUWKHIBapCFpRPqYobCYz8RnAZk2Q1o3xTW7eSRjvxa26lpkLC5FxmSVNbhZ0/NGMjuyizkRkfm2ztzXiAzBtZi+LBfFuRhnUvmsinW3FpVDz9E0Az1PMhBzMtPMIeWIkUYQBEHQm9M+NMzs5Wa2/3w0JgiCILiw6SNPXQHcaWZ/BNwGvMM9k6LzgsNptjRVFVspR3oY67KOs8hSuynqNWe2ibB0lqpkspQvdLdPN1OZZraeSlHDFRGZrBOOMtiTSkTlirAWK8lKLVN1wwHbFMN6UVCnrEQ97kzRIQqRCVRkYHUhjc2EXRdgpjIRi7rUps6RkMEyVwcmZDRrhI1XVGbKZaR1lRlVZXoVKkm9kclOrORT0e/qq6K6Skui/lB6THUtbLzCm1uJ7MRTkVkAoBZSVClOyFTUe1fyMkAtOm8gIsKVHCua0657LG1/JYq21yKTsLREA7U4pmGljynHaUca7v4vgBuANwIvAT5uZj9sZl94RnsKgiAILnp6zWl0I4sHun8zYD/wK2b2unPYtiAIguAC47TylJm9AngxcBB4A/C97j41swL4OPB957aJjwKDcnTq8NjH6XA5Nzx0FfGrwpWVfJDZ6Oy4CPVW0cGNSsqmt7kwSof7k6EY1ott5pTG8oBI0Ccim2uRcFBdVNNDOrK4ELpALUbLjaqoI+QQgGq/aLuQBRrhsipXdSTtUPSntI4JR1etDl0Ur4KMbFWKNi0JSUJIfQALQrqZTUV0/7Xpfpp7M0q0OnRR770QCTGLdf21Y8IVVarrWGQclDKYkNAAikoUbBJFmIaror75euaaW1Y124UrUhTKKjLFs0xISfVIfSeJ+1q4+EB/pyn35nb0mdM4AHyDu3/6lJ27N2b2dWe0tyAIguCips+cxg9sfWDMvffhR7NzM/tnZuZmdnn3t5nZT5rZ3Wb2QTN7yty6N5vZx7t/Nz+a/QZBEASPjB0L7jOza4HnAJ+ZW/w82kn3G4CnAz8DPN3MDgA/ADyVdmD8fjO73d0Pn99WB0EQXNrsZET4v6OdD/mNuWUvAN7UTby/x8z2mdlVwE3AO939EICZvRN4LvBL2+3AMGyrFj4SRUzU3AXApihaJCyAg6GI8s4UHbKFdF0T9lwTGq1lHKbjdaHXiyy1Q2UBFDo0wMZDJ5JlC5ellt1KZSE9JjaYM2kLfblYEZMaIrNpmZnTEDWpZEGdcqQimHVDXdg8m03RdnHeCml51XMns2PpuvK0H0n3MxAaPMD04XQerTyQ+rzH96WTL5kuxlbSvmtEplkTEfbTJqOhT8VclLgHrUqPs9ojoqcnuj9sWcwX1qLfN8Q8hSpiBAxFpLWruUoVNZ/5+pmK+aBhk37eF/tl6AVwYR1HzZNsw45EhJvZC4D73P1Ptrx1NXDP3N/3dstyy4MgCILzyDkbaZjZu4ArxVuvBr6fVpo6F/u9BbgFtPMiCIIgeOScs4eGuz9bLTez/xu4HvgTa8MWrwH+yMy+HLgPuHZu9Wu6ZffRSlTzy+/I7PdW4FaA4crAGZ46nHRRTKd2bQctV4VMIqKVZ8JTWS2qMG8oV9Oh8fiBdLheVCrJXSY6V4yCyzQAmvGJVKYYqAJOwPK16QZmY1GY6XDaJhND9UJllAMQEa6zI2k7Va7GcjEj8Ygkd6o2fCFkjkYl4stgA5WpMl1UT0XjRV8C+EAkn5yJyGBhJ1WR0gC+R4U7C3lLSIUqehq09XQk+nMq2jRaFhkYAVfCh6hBX5g470KiJXNf+yFhez2QtqkRUpQd19HTKpNAKeq9K+m0EZZZgFLI3oU49PE4tdwbwtoPmLIwZxI75jjv8pS7/6m7P9bdr3P362ilpqe4+wPA7cCLOxfVM4Cj7n4/8A7gOWa2v8uD9ZxuWRAEQXAeudBSo78NeD5wN7AOfBuAux8ysx8C7uzW+1cnJ8WDIAiC88eOPzS60cbJ1w58V2a922gTJgZBEAQ7xI4/NM4lTho2X4oi6rmQ+0IUtZm50MZFGhEb6fmHqZgTqVRRGTUlouZYgEZkB53WqVVR2TxN6NgAk/tTPdiWhB4q5iSUFlwsZLRx0v4YrAo750RkuRWpHwCGop21yNqL+HxRZuaNhF6/PEy3aeLYT6h5iswUjyrs5GvC+l2l+y6UnRMYyEkvka1ZpNgRsny7LzFHNVVpRIQ11zJpM1gXBaiEHb6eCctuI/Ytzlm7zfQ4Zw8Ku7HKXJMR9ItFYaUX806zh0Rht0zxrIGo7rZ5RBzTKL0+ygWdSVjNn/hSxledIYowBUEQBL2Jh0YQBEHQm10tT4ElGWSnx/vZWwGa9X7DYISNrRZRwACIUaOLIk4zsczX9TPeRI3wok6PSaoUI33s08V0aFyJYXS9LIpKzdKh+iyTSbMYCalAZGuVSWEz0slEyRxCunERNVu63qgLCdHE+Viy9NhHjYiKHosQdeCIyMZbLAifpfBuWqbWfSmuj8kJYWUV0sUMXcioEEWxClGsWmXebYTEC2CeHuemsImXIhewTUVEuIjyBrDLxfUlaqG7KG5WZoowyYwBB1XW3/SzjZBoAVxIgL6g2iQi1w9nou5F8a3qDOPZYqQRBEEQ9CYeGkEQBEFvdrc85Y5viXwtl0SU5d6Ms+eIqEWsCsWIAixNJhGgN0JKEvWzZyKKeJCps12rmtjCzTEVNbU3RfQ1wEAWlRGSmUicV4gkhlXm90kjVJqpclQJ81OzppOy1UJmKU+IgjiWXv6zsZZOEFLDcXGcjXDijcr0s4cXtHzQqOtLXEoqmZ5KlgiZGuNC3qouT7c5GGbky4fT9tcnhKQpLi9f1dv0DZFZYSVdVxn+GlHMy3JZCFTGANHHpbhXx8f0/SITJrqQ8Ibifsm0s56JWuYiuWBZCVlRRc0DqMSMKo3BNsRIIwiCIOhNPDSCIAiC3sRDIwiCIOjNrp7TsMIot0QX15upptisay24FJrkTERPToUOXonIXoCmSbXkmdKchaVxJrKdtuuK4koi66aQ8BnszUSOCum2WEz7SWUmVfMUjbDWgi5QZGJeoJ6Jwjna/Ui5tfAWUAg7qgsrKlWmWpRIszudpuf9qIiAZkOci4xfuBDZZxtVeEdkVnZRxAhgJo5J2Sxnh0S0suojYCAi/JvL0z4uRVpWmc0WQMwhzDbSuZNJKfajipaJomMAjbDSmji/jbiOqpWMJfuEsuKni1SE/SyXWVn4c8s1Ye8XEe7K2guAmBfN1cTKESONIAiCoDfx0AiCIAh6s6vlKW+cZkvN5eEBEa28ocdytaoVMxHD9XE6vKwzcoyLobUqrqRqhGeUAvy4GBqLnwMDsc2cpFEU6Zi1Eln2xiqCeZ+wg2aKDm0eSyWRhctXkmW2rmyFmah7cUi1kDlGe9J2bmRsvLJYlZAA1bWgQvEbtR6AkE+VrKis202mSFcjsgOUouB8Jc4bEy2J1iLSWxUoq4fC2ivaAzCTKQtE20V0/0zYwQfaHYuJyHUXUeZCnZa16tudCdu9qMftwi88XNa/3SdHxHVjaX8uP2Y1WbZ+WGccUIlQbZKx52aIkUYQBEHQm3hoBEEQBL2Jh0YQBEHQm109p4GRWswOi9QPtdZtlWvWRqlOORymttVaaPCg7Y/FSFgqhe20rPT8g62I+Q9RDKgRWURLYQcFYDm9NAqhxw4W0vweflzM8aD7Y3ElFVmbzbSPS6HF1ke1V3BpJT0fjbB+zoq03wYqXwnI7KAIzbwRltnBitCmhQbeklppG9HvxUzYrDOTXo2YTzJxmCbml3JVmIplsXBNVS1KF3nG1jxdE/MSI2URFfNT+9JzrtLmANTCcusTcb+JOc0q85Wp7rehSC2iMr1MMxmxi0WRjVdY8TcfOJZ+ttLXl9VqnJCZp8kQI40gCIKgN/HQCIIgCHqzu+WpxvEtVsuxiK4txTKARg2jh+m66w+mOoVl6v6WInrbF/vZ/WZHUukCZCAtpRraisy7ORuvKo6kbKuzE0Ju2yushpmo00ZIEi7qdNeTVJorM1H3Jx5aT9cdpZ93kQ2XTJZbSzdJI6LUh3uF9VJs00UEMoByzVZCkqhE9PXmkUyhq0LIhWJVUxl2RQZVAEQWZtVOEzbaWUa6LStx8LIolpDm1lRxokz2alHouxF13EtxjkTi2bZFKlO2kCpdZapey0jERSq5FaqYmZDHC1EcDXRGjDqTHTlHjDSCIAiC3sRDIwiCIOjN7panzBI5qRBS1ExENX/u81sohJuk3C+G5cIBBGDCBVSIIfhYROJWBzLFa1SxFuG4cSFTmHIFAc1ArLsh6mTvEwWPROGeRtTjBrC9IkndwXS9SrQd4QYDWNgvpKhB2p8+Ee1c0dJJISSNgUhoty5qM5eqSJaQdwBKkciwURHuR9NrdqSi1oFaSCcmnEXjh4VMkvmGKPcLiUjVsFf12lf1eStUUax1ITupgkciYWAhCk0B2EBcxyJKXKhDNBsZKUdEyE9VpgcRUe5r+h6sLhP3oHArKvlxksnAUAp5rFSR+NsQI40gCIKgNzv20DCzl5vZR8zsLjN73dzyV5nZ3Wb2UTP72rnlz+2W3W1mr9yZVgdBEFza7Ig8ZWZfA7wA+DJ3H5vZY7vlNwIvBL4UeBzwLjP7ou5jPwX8DeBe4E4zu93dP3T+Wx8EQXDpslNzGi8FfsTdxwDu/mC3/AXAm7vlnzSzu4Ev7967290/AWBmb+7W3f6h4cDk1MFUIwq7N6LYOsBgb6ofmioGJCKLbVVHb7uIwG5GQrddFjq0yDIL4Kpi0kiE/M6ElisKCYEuTGXeTxtnSexH2VuB5qiIoF4VobjiHNWZSP5CabTiMIvMPItCWVRLUTxrJI4dFZ07y8x5iY+PDwptXET2Dla1cOCHhC26SndUmFhPzOUAFMr2OkyPqRHWc8v4vJU9d7iaZnB10gyuJmzrNHruxIXltxj2y3xrmUB+Fb09OC6s42LOysQcDUDhIrOzqHBWb4h2Zizdk7GwFovvxO3YKXnqi4CvNrP3mtnvmtnTuuVXA/fMrXdvtyy3PAiCIDiPnLORhpm9C7hSvPXqbr8HgGcATwPeYmZPOEv7vQW4BaAQefeDIAiCR845e2i4+7Nz75nZS4G3ursDf2hmDXA5cB9w7dyq13TL2Gb51v3eCtwKMFgeuG1ViVS9omFGShLqRy2iPAsVOZ4pslMra58I6bYmPTUmCucAFAupnDPbTKPHS1EPPGf3Q1gLy30iOeFUdZKICM9cauWisBWWYggvCi4Vg0wNaCH91CqyuBT7dr3NStgvlR0U0XYVsl8sZq45kfhuJBIwKpmkVucCqFWNcRGFPJulGpzXOoGjLQopStXJ3pP2R7Wg5RCh/FKL69gWVBEmceyZgknlkojaV+vuE0k7M1b6WsmNi+l+TFhhi0xB741NIaOJ9QoRye+qWBxQqsJlozN7DOzUT/FfB74GoJvoHgIHgduBF5rZyMyuB24A/hC4E7jBzK43syHtZPntO9HwIAiCS5mdmgi/DbjNzP4MmAA3d6OOu8zsLbQT3DPgu9zbvMFm9jLgHUAJ3Obud+1M04MgCC5dduSh4e4T4Fsy770WeK1Y/jbgbee4aUEQBME27O40IoXBwqm6opo/cFHUBRCGOzBhr0WkkyhXM11biyJQyra6ININZDKwNsJKOxiIQkSi+EwzzRSAEdl8fSO1+5XCkjlV/ZHpDtXHLvRhtU3LWWaFxqtSSjTSwqz7WKUckbbGoZgnEWlepFUZmAkdu9rTr3BOORBWZaC6RiycinQnJ8Q8S5OxBoupjkpYP2fH08/P5PFAIXT4RmWFFdmahyvpuZyllysA1T6RqkVYzP2IOJfi/gVgSRyTOMcz9f2TmXspyvR8FCIljcrmOzV9HRdDMV+YyaKUI+xFQRAEQW/ioREEQRD0ZlfLU+aObclgW8hISR1pXYtAa0T2WEpR2OREJsPmvtSG56pGr5CSCmHXBUDYQacPp40v9wvLbZnJhimG2yakExWUPRDy0lgMqwEKISEUI2FVFDKJktAAGlFfXUlRjThvReZaUHKOF2LdE+kiFxG3KqsqgC2n26xFv7ulHW8iMzFAcTyVOWohnQz3pp+tdd0v6nvF+RB1w8uV9FyqyHOAiciIq6KVp6I/BvuF7TRT7n0qJDObqXMk+jNTsM1EcSMXUtBI2mPlJqWNuFwRX9lrQrIqdOh6KTIWqPrm2xEjjSAIgqA38dAIgiAIerOr5Sl3aLaMvEQOMBkZDGCiqE2h3BwieVyjJCcAlQhQSACFiL5ujmmtoBTunPEkbVPlYhgr2g7gx9Ih60zVgK6EK2kiot7JaAWVkH2EfGFCphgf1OP6UkiApUhIp7xbKmEgaCmqFEkYG5WYUUWJZyKLETXkEVG8KmNAuaL7uD4mnFKiUNb0oMgCIGrFAzRCyrKRcDCdSPXHietrbqAchyKJoQlZ0ZWEl/l287HoDyGJqp/UvqK32Yhb04V2W4xTiXgqZDmAclk49kTCwUZInbmySo2QzV1E2G9HjDSCIAiC3sRDIwiCIOhNPDSCIAiC3uzqOQ0cEmeiCHrNZblVEyAuLLuNyHZqKkIUKJU1T8wVoAq4iCJIAIg5jYHI0Kl110w2TJGhExGRPhMFbRb2pf053czotmqKSEwBmLI6Luj+UFHqsuiROL/1ZmabQttX0b2DvcJSORUFfrI6stiPyDhcDIQV9fC63KLS61VRrGokovszvytVhL+K0FfzHJXoNwAXNvNqQdybotiTiXk0G+lz2Yj5j1JYzyciA3RuDlBN2amiZc16Op+inNsAiIzaUxGRXtViP6ooFW2SjK3YGT4FYqQRBEEQ9CYeGkEQBEFvdrU8ZQaDLQWSxIgvW3t76kJ2EtUAG/HxgRgyAmyKIfhAyAIudJtS1FsGXZSmkW5UUdQlo3jVKgJbSRJ7UguhKlRFRlqbieF+ISyIqrbyQNpooZmpPhbJ9MRvJllQC3BRxEkVtPGxkKdcyDYZU2SzIeQLpZ6WojBToyVAVRu+EVZpJVmpcuugkzA24povhTzl69pu7OIcbx4VSTJFVoZiWZzLTCGiRlxz4zVVN1xY7pe1ld6n6brK8usDkelB1JoHaETFuEocey0yQvjWWIOOUuiKRUbmzREjjSAIgqA38dAIgiAIehMPjSAIgqA3u3pOwx0mW+19QrOezjLpKESmSFNpDYQ+7aqAOzBcFesq36myJWayurqw/BYiBYFPU8+tN9puLAs+id2bsK3aRCwbZHyFq2Kj62m/VysifcKG1m1NpO3YEMdTCjtnUevzVgs7qYv9mLD2utCci6x1U1idxWE2onBYuVdnNvUT6XkvR8KSLWzeg1XReGA2SedERsLaOxb7bjLzOWruZiAmBmTaDJGl1jOW22JJzReK/SjbvMigDOBiPsdEgbNGzG8pyzyAi0nERhT5KpW1N2OlH4v7YCGXOydDjDSCIAiC3sRDIwiCIOjNrpanzKDYKiEIe2uuG0oRheyFsjUK+2LmedyI7LOFGHIqlaTI1BJ2ZWc9JjKwqqyumZrYJjKRFsLi2ogIdykfiMykAAMRvT0T1tFaFDeyzHnzWhyTkgCELbGptMSjimIhsga7sEWXotBUM9P2WOHyxkTGAGVlrYVtFKAciIy2IpJf+a9nY51ZuRCRyetCUh2Je6hZlZukEYWhlG1VlM7W5zdXW0jImi5qjNtRUeNbl2GX2RKqxbShJuzojSwMB2Ulri9xnLU4F4W4PgAWV9I2TdZUtbk8MdIIgiAIehMPjSAIgqA38dAIgiAIerOr5zRwYIu+7bXQrEWFLYBGZbkVqQ6aKk1xacpaBzSqUp0oDM+isNvJLQKF0ElF28tBeuy1yv4KmEi96RORWkTosbUJy23OAvigSJshbLg2E7ptZp5EpQ2tlsW6KhvuRqY/xLyTi/QPpTjORqaJ0JqzqjBYq32Lc1GKFDcALi2qok3ColqIfm/XVXZ0kZ5DWZ1n+mvH5HWssv6K60v1sbJzA8WqmA8Sc3PVspg8EZXvAFhRWaXFnISYX2qEFR7ABun3im+K+0WcdxdzpwC1pxNsnrER59iRkYaZPcnM3mNmHzCz95nZl3fLzcx+0szuNrMPmtlT5j5zs5l9vPt38060OwiC4FJnp0YarwN+0N3fbmbP7/6+CXgecEP37+nAzwBPN7MDwA8AT6UdP7zfzG5398M70fggCIJLlZ16aDiwp3u9F/iL7vULgDd5mxb0PWa2z8yuon2gvNPdDwGY2TuB5wK/tO1ODJotFshCRG7OMvbHYlUMT0XmzOEV6VB9er/wiAKVkAp8Of38TEUWiwhkgGpFRa6KAvIi8twyMomyhDbCOjpYTpfZhtiPkMsAEJHrSnZSFuA6J9iNVZZdUbjHhFyXsTUr6URtczZMI6grIbFkEitTK4lIJCxQmZknIhIfoFwRy0Wg90BkdR0fyxTPUteSiIC2gVgvI1V6KdopzrvlhdpTKI5lsumq62spvRbcxbGLYksAhWjSTNwvlWh7dZWOuvcjIsu2yLZMI+6hYUa6FQXn/Jho6Dbs1EPjnwDvMLMfp5XIvrJbfjVwz9x693bLcsuDIAiC88g5e2iY2buAK8VbrwaeBXy3u/+qmX0z8Ebg2Wdpv7cAtwAUmVxNQRAEwSPjnD003D37EDCzNwGv6P78ZeAN3ev7gGvnVr2mW3YfrUQ1v/yOzH5vBW4FGKwMnS2STuPpWL9a1VHAxSx96EyEG8TuF4ViVG1jYCYK0BSLYigpHA3KwQMwFRJCKSSvRkQMe6Zec1EJB5KoOyyMUiCi5nPSmiz8I/wZjUgYKOt+Ayr/pHKDVQtKBpOblMN6WcBqIqKNVcR+xqVVinaqiOGpiGavhNQHwFicdyEFTYUkmbuOlUJkIlq6UFHeot8BpsdEEsZFId0I1UjJYJk6RJg4caX4gTmbikwPmUSTLlxmgz3p+TBxgdWHMpKoaFMjzqVN035rxP0LgKyPfhG4p2jnMP569/qZwMe717cDL+5cVM8Ajrr7/cA7gOeY2X4z2w88p1sWBEEQnEd2ak7jHwD/wcwqYJNOTgLeBjwfuBtYB74NwN0PmdkPAXd26/2rk5PiQRAEwfljRx4a7v6/gL8iljvwXZnP3Abcdo6bFgRBEGzD7o4IN/AtWW0HnnrmapVaFJiJ6MtKaKe1cKy5KJIDUB5INeJGzEkMLhMa6wmtfRaLYk5mKnTbBZH5dpyZa1CbPCT6Q0TCuihok6m7QyWKYqksuS48jZbRbUfKwiwKVSl11uQkDbjIClujorJFf4h5m3JJzxU0aqpDzN2U4jqcrenruFpVBbnEfI6KZs8UpSpHQpsXfVRUwuIp7iuAgSqKJbIGN8p1Ks5bKYotASBs9/VxNVEisiJkLNkIK6xymTfK8l/q8zYUX8+iO6ktXa+UE27aym+ZaascYS8KgiAIehMPjSAIgqA3u1qeMncGW4bX9VgUtNmTq1/dM7JZ1BL29czzWBUIEkPrqUiQYhlrnImEi6rYkxrXlyqUFaiFbFUN0/Bva1KpQSUsRMghAKWQwUpRPacWw3qUlbVdW6wq5Jg9ySLKw/paqEfCOropCveIIl++lu57oq4DoBQR4a6yAwh5qhTFmkBbkFXxK9VHUvoEGQ1fLgopScg5M2FlB6SEqaLMq4HQU4Q9Nhd1jyhaJOp+UYqkobly2koNcmETV7Fjo1p/DSt7LQNxzSmJN2elXxD7Wrg4LLdBEATBRUg8NIIgCILexEMjCIIg6M2untPALLFvKvWuWc/kjhCpFgpVL0UVyclY3miElXYj3ehwT6rbzjJ2P99MdediUei24jhzRZiUyNsUwton5ngMUaBHZMNt969SNQgNfyQyeap5DpApWFTW3vpQqg/PZrogDsreq4pAiQy/jciEMRLHA5k0Myo9iCiiNN3IZGsW6V8aMcdUCA2+yWSKLVZEcaSjKq2rmEfL5GopRAEqNcezeSTdz2BRaP1i3gfAN4TVWs0RqSJdmbarO9PU3Fol5qwy8w/Ks6tSo5iLuZdlPfmi5lnKDb37HDHSCIIgCHoTD40gCIKgN7tbnnIotvju6hUh5QiLJwDC7udL6WrNUTEMXtJRryZsmsMFoV8Ie6yp9J5AMUiHnC5qHqsspI1oD4CJYbSqKz0VikQppLHZcd0f8hIUGXFl8Hcm9b3KnFtNRZEdYTf2TJ1tVftbSUTNgpCXROZaVecawJVt9oSoKy+ivItZpkiX2pUo0lNtpJLoRGQsBqhPiOtrr5CihJ10uqm3ORPy6VAIyibtrSlFLg2BiGZ3EYpfq0s28zO7EJmIVX3zQtjEbZYpSlUJmVdKWeJ4VIQ7YKJgWy1qmW9HjDSCIAiC3sRDIwiCIOjN7pangK3RwYti+L+p5CWgGIokZGLVZiEdxw4qXfe3XhJawZpyVIhiPnIQDiZGl7KdwhikJCcAE64VMQLHhiIBWpN+tspcakoOKlaV3CakoMyo2qcioly0XbXJJxm3kIjGr4WEV4pjb4QLxoVEAmDCyeIiaZ+Joj91JtK63C/aJPajEuep3JMALjJaFsItNBW16ktxPLl9NYWIPFf12oVs4yrMG/AyleHUXaDu/0ZE9wNsCgfkSCUnFb/Ti5zsLJIgKpXUKvVdobMQ2glRIO3MAsJjpBEEQRD0Jx4aQRAEQW/ioREEQRD0ZnfPaRg0W3TJmcpyW2jLrainw2Qz9dyaigxuvlA3afPedNVJqoea0Pqr5UxhlYmwPwobX6Eyq2YspoUUmIXeX6RzNyqSdpwrSjUQ2vyJVMstRNS9TbQWrBTqUsxJNGIyyJZ1fzTCWlwsiHVFv5fC3komCldlCFbRyi4s2UNRZAtgfFjMCwi7sso43IjofgAX11ItIvERttMmk4WgXBC2aHEuCzH3YgvCbixs1qALjM1cWGaFFVVc7gCMVtLibipDL+I6tKHwwgO1SE8s59bEcRajjOVffOWf4ZRGjDSCIAiC/sRDIwiCIOjNLpenSqpi76nLNvan6y2sy49XovDOrFhLlg0naTc+4Skn5Db//DMvT5aVf/H+ZNl09s5kWSOivAGUgiCciig1R+VkA2hE2K2PheQlwo0nk1Q+UMnoAMoNcUwqueFU2B839XmbHE3XHayKKPFVUb86F7iugq03005uRHR+JYr55H6ulWsiC4GwqNbTVNJspvp2liXbV9JlsyPi/C5o6bYRhYwQNu1CXGB1Jhq+Fp7wAanso1zRA1GwKJfAsRLJ/EwlcBQFqGZrGXvsRMhOqmCSsNI3yg8OmEhuqr4BygWR9LMUqSuApkp10VLpktsQI40gCIKgN/HQCIIgCHoTD40gCIKgN7t6TqOsFli5/ItOWXbguX89WW+6/iL5+cn0ULLMf++/psvs9mTZFQePy20+89evTZa96TXfnixr3v4VybKZPyy3WUrLr8g4KiyinsmwibBF2l5h1xOStwt9WVknATZOpPtZeGy6Xv1wul61L9W7Qc/TmLBKTsdp4wdLepsU6WSHlelG/UiqGc+EdbPJpIRRevv0sEhTI4o4qfkUgFqkoyhFOgxbFKk4hJ0bwMRESaG6Tuj1lbIgAyyl52Pj/nSeY3FFzG+J8ysLpgGIVC+yOJIoBlYtZb4y1T0o7Mau5hXV/BDQiLQsygrfbCiLuvCIA8Uonb/IFmLLECONIAiCoDfx0AiCIAh6Yy6GULsFM3sI+PROtwO4HDi40404y+y2Y4rjufDZbcd0IR/PF7j7Y9Qbu/qhcaFgZu9z96fudDvOJrvtmOJ4Lnx22zFdrMcT8lQQBEHQm3hoBEEQBL2Jh8b54dadbsA5YLcdUxzPhc9uO6aL8nhiTiMIgiDoTYw0giAIgt7EQyMIgiDoTTw0ToOZ7TOzXzGzj5jZh83sK8zsgJm908w+3v2/v1v3S8zsD8xsbGbfM7eNa83sd8zsQ2Z2l5m9Yu693LbMzH7SzO42sw+a2VMupGOa21ZpZn9sZr85t+x6M3tv1/b/ZmbDbvmo+/vu7v3rLqTjUdvplp/Xc3QWj+e7u+vtz8zsl8xsoVt+Xs/PIzimF3X9+adm9vtm9mVz23mumX20a+Mr55ZfyNecPB67wL4XeuPu8W+bf8B/Af5+93oI7ANeB7yyW/ZK4Ee7148Fnga8FvieuW1cBTyle70KfAy4sfs7t63nA2+nrcb4DOC9F9IxzW3rnwK/CPzm3LK3AC/sXr8eeGn3+h8Br+9evxD4bxfS8ajt7MQ5OkvX3NXAJ4HFuXPykp04P4/gmL4S2N+9ft7JfqWtavLnwBO6bfwJn7+PLuRrLnc8F9T3Qu9jP987vJj+AXu7G8+2LP8ocNXcif/olvdfg/iCnXv/N4C/sd22gP8P+LtqnxfKMQHXAO8Gnkn30Ogu5oNA1f39FcA7utfvAL6ie11169mFcDy57Zzvc3QWj+dq4B7gQNfXvwk853yfn0dzTN3y/cB9W9va/f2q7t9Fcc1tPR7x3o59L5zJv5Cntud64CHgZ62VYN5gZsvAFe5+f7fOA8AVfTfYDY+fDLy3W5Tb1smb/iT3dsseLWfzmP498H3AfJrOy4Aj7n4yZel8uz93TN37R7v1Hw1n63hy22GbbZ2Lc3RWjsfd7wN+HPgMcD9w1N1/i/N/fh7tMX0H7S/rU9q3pe0X0zU3fzyf4wL4XuhNPDS2pwKeAvyMuz8ZWKMdKn4Obx/33mdjZrYC/CrwT9z92Nb3z2Rbj4Kzckxm9nXAg+6e1qo9v5ytc3Ta7ZzBth4NZ+v87AdeQPsF9zhg2cy+5Zy0+PQ8omMys6+h/ZL95+epnX05q8dzgXwv9CYeGttzL3Cvu598+v8K7cXyWTO7CqD7/8HTbcjMBrQXxi+4+1vn3spt6z5gvvjGNd2yR8vZOqavAr7ezD4FvBl4ppn9V+BhYJ+ZnUz8P9/uzx1T9/7ebv0L4Xhy22GbbZ2Lc3S2jufZwCfd/SF3nwJvpdXWz/f5eUTHZGZPBN4AvMD9c4Vkcv19wV9zmeO5kL4XehMPjW1w9weAe8zsi7tFzwI+BNwO3Nwtu5lWi8xiZga8Efiwu//Elrdz27odeHHnlngGrbxwP4+Ss3VM7v4qd7/G3a+jnWD8bXf/lu5X0e8A35Q5ppP7+KbuM4/qF9RZPJ7cdra2+5yeo7N1PLSy1DPMbKm7/p5Fe/2d1/MDZ35MZvZ42ofct7r7x+Y2dSdwQ+eUGtJed7df6Ndc7ngupO+FM+J8TqBcjP+AJwHvAz4I/DrtRNZltBPAHwfeBRzo1r2S9lfIMeBI93oP8Fdph5cfBD7Q/Xt+95nctgz4KVq3yJ8CT72QjmnL9m7iVPfUE4A/BO4GfhkYdcsXur/v7t5/woV0PGo7O3GOzuLx/CDwEeDPgJ+fOw/n9fw8gmN6A3CYz98r75vbzvNpXUZ/Drz6Irnm5PFwgX0v9P0XaUSCIAiC3oQ8FQRBEPQmHhpBEARBb+KhEQRBEPQmHhpBEARBb+KhEQRBEPQmHhpBEARBb+KhEQRBEPQmHhpBcB4xs6d1dRAWzGy5q6Pwl3e6XUHQlwjuC4LzjJn9a9pI5UXaHEb/ZoebFAS9iYdGEJxnurxJdwKbwFe6e73DTQqC3oQ8FQTnn8uAFdpqbQs73JYgOCNipBEE5xkzu502nfz1tFXXXrbDTQqC3lSnXyUIgrOFmb0YmLr7L5pZCfy+mT3T3X97p9sWBH2IkUYQBEHQm5jTCIIgCHoTD40gCIKgN/HQCIIgCHoTD40gCIKgN/HQCIIgCHoTD40gCIKgN/HQCIIgCHrzfwCe/AXWRI1HIwAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "(da.sel(band=[\"B04\",\"B03\",\"B02\"])[1]/2000).clip(0,1).plot.imshow()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Now, let's modify these default values by passing the required modifications as a dictionary." + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/home/dmontero/anaconda3/envs/cubo/lib/python3.9/site-packages/stackstac/prepare.py:364: UserWarning: The argument 'infer_datetime_format' is deprecated and will be removed in a future version. A strict version of it is now the default, see https://pandas.pydata.org/pdeps/0004-consistent-to-datetime-parsing.html. You can safely remove this argument.\n", + " times = pd.to_datetime(\n" + ] + }, + { + "data": { + "text/html": [ + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
<xarray.DataArray 'sentinel-2-l2a' (time: 2, band: 3, y: 64, x: 64)>\n",
+       "dask.array<fetch_raster_window, shape=(2, 3, 64, 64), dtype=float64, chunksize=(1, 1, 64, 64), chunktype=numpy.ndarray>\n",
+       "Coordinates: (12/47)\n",
+       "  * time                                     (time) datetime64[ns] 2021-06-13...\n",
+       "    id                                       (time) <U54 'S2B_MSIL2A_20210613...\n",
+       "  * band                                     (band) <U3 'B02' 'B03' 'B04'\n",
+       "  * x                                        (x) float64 6.011e+05 ... 6.024e+05\n",
+       "  * y                                        (y) float64 5.66e+06 ... 5.659e+06\n",
+       "    sat:orbit_state                          <U10 'descending'\n",
+       "    ...                                       ...\n",
+       "    title                                    (band) <U20 'Band 2 - Blue - 10m...\n",
+       "    common_name                              (band) <U5 'blue' 'green' 'red'\n",
+       "    center_wavelength                        (band) float64 0.49 0.56 0.665\n",
+       "    full_width_half_max                      (band) float64 0.098 0.045 0.038\n",
+       "    epsg                                     int64 32632\n",
+       "    cubo:distance_from_center                (y, x) float64 894.1 ... 887.9\n",
+       "Attributes:\n",
+       "    collection:           sentinel-2-l2a\n",
+       "    stac:                 https://planetarycomputer.microsoft.com/api/stac/v1\n",
+       "    epsg:                 32632\n",
+       "    resolution:           20\n",
+       "    edge_size:            64\n",
+       "    central_lat:          51.079225\n",
+       "    central_lon:          10.452173\n",
+       "    central_y:            5659638.0946523\n",
+       "    central_x:            601722.4825156148\n",
+       "    time_coverage_start:  2021-06-01\n",
+       "    time_coverage_end:    2021-07-01
" + ], + "text/plain": [ + "\n", + "dask.array\n", + "Coordinates: (12/47)\n", + " * time (time) datetime64[ns] 2021-06-13...\n", + " id (time) " + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEhCAYAAAC3AD1YAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAAsTAAALEwEAmpwYAABSKklEQVR4nO29e7xlV1Xn+x37nFOVpJIQAiQqDwEl0IIEkgBq0xAIHUFt8Yr4QVCwoW8aG2yufcFGm3QjvtDuK4KIITd0NLagSCcQXwkRDbHlIaSFkJAQXiEkEorwyrNSVWeP/mOtSnbt/Z2n1q7apx454/v5nE/tWnutueaaa6491/yNMceIzKQoiqIoFsHoQFegKIqiuPdQg0pRFEWxMGpQKYqiKBZGDSpFURTFwqhBpSiKolgYNagURVEUC6MGlaIoimJh1KAykIh4SETcFhFLB7ouRbHRiYirIuLUxnenRsQN+7dGxS5qUGkQEddFxDN2/T8zr8/MIzNz9UDWq0VEbI6It0XEFyLi1oj4WEQ8a2qf0yLimoi4IyL+NiK+feK7/xYRn+6PvSYiXjh17NkR8amIGEfETw+oz+Mi4vL+XJdHxOOmvj8pIi7rB+ovR8Qr1ijraX19vxkR1zXO9Xf99zdExJn99hf05d8WEXf2dd/1/9v6fV4eER+NiLsi4vel7Gabyb6/GRFfjIhb+vvwi/O0ydS+r4qIK/v78fmIeNXU978cEZ+IiJ0R8dpWOf2+ERG/ERFf7f9+IyKi/+6EiHhPRHwlIr4WERdHxCMb5bwvIjIilie2PbRvlzv6dnrG1DE/FxE39W3y3yNi89Bj1yIzH52Zlw7dfxFExEv6et7a99m/jIij+u/WvF/9PhERn4uIT85xzkv7Nj9xavsF/fZT+/8f07fvTX0dro2IV0/snxHxnXt77XORmfUnf8B1wDMOdD3mqO8W4LXAQ+leFn4IuBV4aP/9/YFvAs8FDgP+K/ChieN/CXhUf+yTgK8D3zfx/cuA04CPAj+9h7psAr4A/BywGfj3/f83TdRlK/CC/vujgH+2RnlPBH4KOAO4Tr7/JPCrwBLwHcCXgB+e2udU4AY59keBHwF+D/j9qe/WbDMp65HAlv7zA4GrgB8d0iZS1s8DJwHLfblfAJ438f2LgGcB7wFeu4f78W+BTwEP6uv1SeClE237EuBYYAX4ZeAaKeMFwGVAAssT2z8I/BZwOPAc4BvAA/rvvh/4MvBo4L7ApcDrhxy7j8+C3ut9LPOp/bU8vv//sf09OGrI/Zoo4zZgG/CEgee9tL93/9/EtvsBN9E9Q6f2284F3tm384juWf6xiWMS+M5FtkmzzvvjJIfaH/CHwBi4s+8EP0/3Y333A9Xf7F8BPtDv82f9zf4j4BbgI/Q/6P3+jwIuAb7Wd5If3w/XcQXwnP7zGcAHJr7b0l/foxrHXgj8v7L9f7HnQeV04EYgJrZdDzyz//xrwB/uxfU8Ax9U7gC+a+L/fwr8wtQ+a/7Q9Pfy96e2zdVmU8c+EPgE8PND2mRAeW8Cfke2/w/2PKh8ADhj4v8voTE40v1YJnC/iW33Aa4FvmfqGTgBuIv+h7Xf9nfcM2C9Hfi1ie9OA24acuyA9riO/qWPblD6fboXoU8Cr1rrXu/NH/BK4N1z7D9zv4D/Tvf7cD7w5oHlXAr8Z+AGYKnf9nK6l6AbuGdQuRL4kTXK2W+DSslfQmb+FN0D/6+yk7x+s7Hr8+jeoB9I94b8Qbo3hmOBq4H/AhARW+gGlLcDx/XHvSUivssKjYi3RMQ3Gn9XDLmGiDie7sG9qt/0aODjE9d4O/DZfvv0sYcDT5g4dl4eDVyRfW/uuWLiXN8DfC0iPhARWyPizyLiIXt5LoDfBl4YESu9dPO9wF/vQ3m7WLPNIuLVEfHnkwf0226je+C30N3zXWU12yQint+6t71U9S/Yt/vx8Yn/fxy57z1Pofvh/+rEtl+j+xG7Scr9XGbe2ijbznt8RNxvwLHz8F/onr/voJsdvWitnSPiijWer7c0Dvsw8P0R8UsR8c8nZTwpf+Z+RcQRwI/RDSp/BDwvIjYNvL5/ohssT+///0LgvKl9PgT8akT864h4xMBy14UaVPaNczPzs5n5TeCvgM9m5l9n5k66t+XH9/v9EN0b9rmZuTMz/xH4n3SyygyZ+e8y85jG32P3VKmIWKHruH+Qmdf0m4+kk3Im+Sad9DTNWXQP+MV7OleDPZ3rQXQP/iuAhwCfB96xl+cC+HO6B/ZO4BrgbZn5kX0obxdrXkdmvj4zf2jyy8x8ff/9SXQz3l3H76mst69xb19L96yeu1dXMXvubwJH7rKr7CIiHgT8LvAfJradAvxz4HcGlLur7KMa3+/6fNSAY+fhx4FfzcyvZeYX6WYJTTLzsWs8X/+ucczf0UmlJwF/AXw1In4r3HHntczerx+lm5m9tz9+BfjBOa7xPLoXp0cBx2TmB6e+/1m6Z/7lwCcj4jMxZVPdX9Sgsm98eeLznfL/I/vP3w48afKNiE6j/pZFVygiRnQ/ZtvpOtgubgOOntr9aDq7y+Tx/xV4DJ08NyiEdUwYv/sZx57OdSdwQWZ+JDO30dlzvi8i7hMRvzhR1lkDzn0scBHwOjq7x4Pp3ij1x2FOBrXZNNnxj3TX+Uv7UlZEvJzuzfQHM/OugfWeZvrcRwO3Td7fiHgA3Q/eWzLzHf22EfAW4BX9i9Keyt1V9q2N73d9vnXAsfPwbcAXJ/7/hb0oY49k5l9l5r+iUyKeDfw08G8m91njfr0IeGf/UrmN7qVyzRnVFOcDT6d7pv9Q6nZnZv5aZp5MJ8O/E/jT/vnYr9Sg0maROQG+CLx/6o3oyMz8Gds5Is6a+qGe/GtKIP2b59uA4+lsKTsmvr4KOHFi3y10csHkFP2X6Iy/p2fmLUMvrr+WXX/X92U+dupN+LET57qC3dv37s/9g7GrrJcOOP3DgdXMPK9/YG8A/hj4gaH1X4M9ttkeWO7331XWWm0yQ0S8GHg1cFp/XXvLbtfRf5687/elG1AuzMxfndjvaOAU4E8i4iY6OyHADRGxS955+C4PKCnbzvvlXlrb07Hz8CW6l4ldrCmlRueO3Hq+9vgik5njzHwf8Dd0L2C7ytX71c8Anw78ZO+ddRPdzPoHIuL+Qy4wM++gU0N+BhlUpva9hU6y3AI8bEj5C2V/GG4OxT86jXLSuPlQZg31/2bi+90MvXRG5c/0n4+ie3v6Kbpp7wqdzaLp8bSXdT6rr/eR8t0D6OSF59C90f8Gu3t//QLwaeBbGmVv6o/7e+D/7j+P1tj3C3Ty1ma6t6tJ76+n0xlVH9e3xRuAv1vjukb9+Z7Vl3PYRFlH03kNPb/f71vobFu/NlXGqbj313Jf3q/TPayHTdzjNdtM6vhv6bxvgs6r6kvAvx/SJlLeC+hsGNpH+nY7jM5m8yv956XGvi+ls/E9kO6t/iruMaYfDfwDYjjur+NbJv6eQPcMPHCi/T8E/Lf+/P8Xu3t/PbO/hu8CjqH7EZ70/lrr2FPpJn2tPnEd9xjqfwN4f9/2D6J7aVm0of7ZdLbQyfv7FeAFe7pfdM/W1VNt+S3A54Cf3cN5L6X/nenv3ZMnvps01J/Z359dz+l/onvGjuy/L++vA/3Xd6Lr+47+SvZhUOn//0g6LfUrwFf7B+xxC6zvt/f120YnLez6e8FUna6hk2UuZXfvtKTTfCeP/cWJ7y/t95n8O3WN+jweuLw/1/+md8Wc+P5n6Lyhvk7nOffgNco6Vc596cT3T6d7i/5m/2D//8ARUoYNKq+Vsl87sM1+Efir/vOITob7Wt921/bfx5A2oftRumri/58Hdkzdj7Mmvv99qfdP99/9Czp5a9e+AfxmX7ev9Z+j/+5F/bG3T53rIdJWD2XWpfihfbvcSefV+IypY/4DnSx8C52NYfOQY+lewP5+jT5xHfcMKkfQ2Ry+wfp5fz0FeB9wM51Edy29Z9+e7lfff2YGDzqv0o/u4byXMvE7M/Xd5KDyGjoPsFv6e3wpuy8JuHtQkb521lTfuoqJ3415/3Z1rKIoioOGiDgH+NPM3FtnkeIAUYNKURRFsTCW97xLURRFsV7065qMZ2XnynxIUTOVoiiKYmGUS3FRFEWxMDb8oNJH9twaEVcO3P/HI+KTva/72/d8RFEUxcZhw8tfEfEUOve/8zLzMXvY9xF0K1Wfnplfj4jjMnPr/qhnURTFocCGn6lk5mV0ft13ExHfEREXRZfz4u/6eDvQLfr73cz8en9sDShFURQTbPhBpcHZdIuVTqZb+LgrcukJwAkR8fcR8aGIeOYBq2FRFMVBSLkUTxERRwLfRxeMbdfmXWGul4FH0K3OfhBwWUR8d2Z+Yz9XsyiK4qCkBpVZRsA3MvNx8t0NwIezC9T4+Yi4lm6QWUSY9aIoikOekr+myC7C5+cj4rlwd17pXZFW3003S6GPLnoCXVC4oiiKghpUiIh30EW1fWRE3BARL6ELuPaSiPg4XXC1Z/e7X0yXnOeTwN8Cr8rdM+QVRVFsaDa8S3FRFEWxODb8TKUoiqJYHBvaUD9aHuXS5qkU0+sxcYs977Ku5zesTvOce55rGkJrxhwDTySHt2bhQ4u0Hee6lXZ622iFLjXOZGWOZduh0ueMA/0MDmVf66nHr8NDuA7XvvPOnTdn5gPsuw09qCxtXuJ+/2z3bJ4hD2i2bvTOgT8Q1sojv9O5ar8Qst949tzR6D1a+6XZTWE/eGMvM+RHT9sprJ6z5GqjjZdlbytAMqiv3rWqRY5Whg0Wo5XZiXyMhk/ux9tn72XulG1y8qX7rHihO2Tb9oH9sIX1ebsdI9voJxpLX9B2t3PLvWzVSfusPVty21o/3/bSYc+W9dlWmfa82vF6PU2kona4tPs8/SPkvm/9x5u/0Nq/5K+iKIpiYWzomQow8yY9ljeF1qCe8rZvrwr25hXyBt/E3jTszafxnmSzCpupsGpvPl6mvk3aK4rNaHQy5jO0SCnUmiNmZyXLW/QGkTKBybGc307dfL21E8m1yzadEe1ozFhlppSbZJscr/0AtC/ZBenbts5e0H6TO2fLHNvbupdILM9e58j6l80grDn1ugE5jz8v+zjTsOdlDrVAW0ruh8q9TXl1X3XxmqkURVEUC6QGlaIoimJh1KBSFEVRLIyNbVOJIJZ2H1fN02Hc0LfDPEpE/zQdPho6pUmqucN0VimzdTdFDx6JDj++aw7N3M5v3nCyn9qtWh7FS8PKVO+vHe79FXLjzP14LF5V0bD9rIqnl16nebOZraLhARXL0h7mOWc2DbOZ0fA4tNdNsedYP4LGfZNrStP2G+akkH5sHnbmlWlee9HwLQ+ze8l+q2LTWF7y9hhLnXZYg4jn2mjJbYP6bJmjl9qIGl574jE5/Ru5J2qmUhRFUSyMGlSKoiiKhbGx5a9McnpaqrPCxjTZ3C7Nd1GG7rG67jnmyjmSqXvTpdjkFZHUQupu7q4AYzne6mlShNbTZKEWJg2Ju2ps8sPN5XRJ23OWVXM9Bsars7LBaDpaA/h6NdX+GlKVLZ4UlU9lNlcDXXKxBYTm2t7ox7aoUd2HrX+0vF2l/qsiTev91cW+LRdrqfx22U3u+c6GVLW6fbbTpkmEQxcVg0vQsmbWpHZ7foHZ30P8d2Et1nWmEhHXRcQnIuJjEfHRxj6n9t9fFRHvn9h+TES8KyKuiYirI+J7++3HRsQlEfHp/t/79ttfEBFX9Of7wES4+qIoimI/sT/kr6dl5uMy85TpLyLiGLpUvT+cmY8Gnjvx9RuBizLzUcCJwNX99lcD78vMRwDv6/8P8HngqZn53cAv06UELoqiKPYjB9qm8nzg/My8HiAztwJExH2ApwBv67dvn0jZ+2zgD/rPfwD8SL/PBzLz6/32D9Gl+y2Koij2I+ttU0ngvdHFJHlrZk7PHk4AViLiUuAo4I2ZeR7wMOArwLm9jHU58IrMvB04PjO/1B9/E3C8nPclwF9ZhSLiDOAM6Nwhp71LNWhda+jV6LDDgtm1QkSY22SYm7LJwQ2JeKgbq7mmLpmbMB6scTwwdMzoMHHvbLgtasgNC7MiPTm3uRa8JOFbQm6S6cstfXvpiFkxe+kwuZni5mz2NQu+2O08u8m0eQ1m3Opz9itgNiYN6OhtnBbk0rzDLUDmSiu6tLg0S1+yoJv2DLUYi01lpwVw1XBBXmZa/zZbmhSw1HD/Xdo0u327dJCxPC9Nm4q51s+Zc2u9ZypPzsyTgGcBL4uIp0x9vwycDPwg8P3AmRFxQr/9JOD3MvPxwO3cI3PdTXYtsNsVR8TT6AaV/2gVysyzM/OUzDxlZP7bRVEUxV6zrr+qmXlj/+9W4ALgiVO73ABcnJm3Z+bNwGV09pMbgBsy88P9fu+iG2QAvhwR3wrQ/7t1V2ER8VjgHODZlea3KIpi/7Nu8ldEbAFGmXlr//l04HVTu70HeHNELAObgCcBb8jMmyLiixHxyMz8FHAa8Mn+mAuBFwGv7/99T3++hwDnAz+VmdcOqmTOriA1+ckiy3bbZaOtaDVJq+VCa8P8wMCh0ViFrTKdTYlFarprm+kT+OrqzbPdyaL/hrirWhsBjMWFdrRJb9JsmQ3d0tppbFKESDtWd3CZbsdt4odq6T+WZ68nG+dRWUr2dddQL1IjEstemsOnNdu358By6xwusmND+tO+ZH3eNCh7XhqyUlrf1h1l2zx5SqSeKiE3cg0tS//eKdtWTaJsdIaQtl9qRE1osZ42leOBC/pQCMvA2zPzooh4KUBmnpWZV0fERcAVdD9952Tmlf3xPwv8UURsAj4H/Ot+++uBd0bES4AvAD/eb//PwP2At/Tn3GkeZ0VRFMX6Ea2UqxuBlSNW8thH3He3bfPMVHSVlsbekWNbg/8+zFQ0E2XreHn7GYsRuWXstplKyAxC3y4tA2Dj9Wa8UxYVjuQ81o+3eXuMxMCZanW1haxzGNDtDXHgTKW16NRnKnIa6wutxY/q9GGGenG6aM1U1BBsF7+vMxVbqLhvM5VxK97dgDJbM5WxxFfTBaIyU1lpzLgPk36zTWYgNlPRhcqg0wybqXz5H7Ze3npp39gr6pGOZfevJSuZZ5NJGdZ7WiXacycPmXrXtALxidSmTmoiX+001xEaq6PtoTeZb45gh/rjJj+Ymjq48dzkjmEBJW3waclfIT8aek2WUEt+CMaNPmNJy9IGTzu8qVQNG6g09a8FdARdGb4kA4iN23mHP3D6bmfSskmp1udaY4dJh9KPddhsFjqsL8n7Uuux5g5Z0W9y8ar0zaWj/KdfvQ4bgTdblPtTURRFsTBqUCmKoigWRg0qRVEUxcLY2DaVkOQ9EuVThd8GZkw0g20zodZ4mGZu0WpbLrQq/lo01LtkNW7D+B9m7DYrpdle1MzSWOGrphIpUwzbrTbWVfoi2GuVGu6dmlxNl5BbZFrZ1vDkTtHRV8SwYC7S2XDb1jrpknxbud9KdCVFqmHb7H2Nfiz300xc3pXMjtcyQsrRdt/VJtKq+zC7RFjfbFhVbOtYzFEW2WGeiA1N42SDmqkURVEUC6MGlaIoimJhbGz5K5mVM2xW2PLF14CSsskknJafuEk75kspPp/tFdOzjG0NhSS6GokbaH+22SoN1TxMNmhIHhoE0VZmq5vy8AaxREa2WnzUWNuw8zbRHUTaseRq6vLdiEyo7r+WuEsli0Ybq3v4wHM3lWEp0yIp2AOz4jqdJcUaGtVCNdd5PGXNBd8ScrUkRmt7aTyVLVvLz+SZWdos5xlJ/7AlCbhL8tKcMRJrplIURVEsjBpUiqIoioVRg0pRFEWxMDa2TYVZvdLsH2Eun3IsaCgjUkOVNDRzSzp0l+yo9oI54lKZ3i/ukK0IDerWa+aPw01klkJ3uDhvNh213Zjbd0vv1zhQjX2ni2zFV7NEW2J/UdORbDNtGwAJfKxxpSyybMOOZ31WQ7KYLavRQWy7efCGn0jLzLvE3iC/YHG4ubbLfq2MWva8Dky+l41+bC7z2o/NvtUIDeTx4ew8EjpKS/RQOkOfjbvrNd/uRVEURdGmBpWiKIpiYWxs+SuYcQH0VdR+uLkF68pq268Z7dZcY6VMO7i5gnx221iiy84VjFSkA5N2NBGRJdSaw23bXGhXt4u7aSO69EjcLi2Eu7kkq3s3Dddr6wuW997SELSkKovGLyumVUJqyYGmzpqcKC60LWlGGejl3IquEJvl/CJNWxuP5ETNCAOWmkE691h+GMYSlaKr0+y2kH4YIkG33Ng1n7y4XavbduNx83NVlOKiKIriAFGDSlEURbEwNrb8BbMtMIdsYKuwlTkkghRPL52+2rS/sbzZkvRYAiqVXOZYMW0y3+ptsprXgu61pvgi0y2Jh9ySecLINUIjOZLJMLKyu+XdM5YkYbpS3pI1SZ/Te97VYNgmy7xowVJxycUCm/qq9IZnpGbClMPFS23U9LYcdu3aco0AnXoeeV7sGWLT7KbRYd7IFgjUJDG7dE10hzs8ajBOjTw5PEf9vNmBa6ZSFEVRLIwaVIqiKIqFUYNKURRFsTA2vE1l2p1zJmkXjVXQQK4OjHxq4uccdpaxHW/aZyOZmNlPlsResLpNdN9GIiO1F0hDjWQ/VWhbNhV771EXb7GT6Ipjr8BINPedd0p7NGwII3mSNLi0tae6keppPMKARQ+WRho1fGgtcjLLYh/TKA5zREKQ/XSleStogbrLiquv2AvMVd+OBfxV2+opK/xbWBNbPbU1tzUMm+ZJrkaZ2U1NN2Wr0xzRvqFmKkVRFMUCqUGlKIqiWBgbXv6acS8dmBRKj6Wx2ltW+GoyHRqykuoGdnBDUjOXUUuVblPihg4TluNeJIIdt8mqY5lOWx5tgKXNch6RqubwmiTE3XZZ2s5kpaVGG5sSKmqNBufTW9kIKGkKlrmBmo/0qBEZQr1QLcKByaONRl6VTmuKb0iAzKYHq3kUW5BLOZElpRs1lgTYPdYoEJZ0rBFxAXEf1vNrZIeG/KVu0lLPldmf+Wb0DIkIYCaBtVjXmUpEXBcRn4iIj0XERxv7nNp/f1VEvH9i+zER8a6IuCYiro6I7+23HxsRl0TEp/t/79tvj4h4U0R8JiKuiIiT1vPaiqIoiln2h/z1tMx8XGaeMv1FRBwDvAX44cx8NPDcia/fCFyUmY8CTgSu7re/GnhfZj4CeF//f4BnAY/o/84Afm8drqUoiqJYgwNtU3k+cH5mXg+QmVsBIuI+wFOAt/Xbt2fmN/pjng38Qf/5D4Afmdh+XnZ8CDgmIr51f1xEURRF0bHeNpUE3hsRCbw1M8+e+v4EYCUiLgWOAt6YmecBDwO+ApwbEScClwOvyMzbgeMz80v98TcBx/efHwh8caLsG/ptX5rYRkScQTeTYbRpSdxgh7s4hrmRmruq2TRaya8s9IJJqmKniYYLreq8FkbDXFNbwYPN0GN2AIswo9Fm53DvVN1Y9G0Jo9EdLgmkzM6i973hiim6t3rwWiQcObfZrKDVF2xHO3XDHiSu05q/SlzWm0ncZF+ruun1FkUaYLxD7A1mqJGHyNqz5bY9llA8drzZVFZbyebk2RzJw6W/KcteUYvGvLpttjMsSUdu2X400lPLONlgvWcqT87Mk+ikqZdFxFOmvl8GTgZ+EPh+4MyIOKHffhLwe5n5eOB27pG57iY7y/BcV5yZZ2fmKZl5ioU7L4qiKPaedf1Vzcwb+3+3AhcAT5za5Qbg4sy8PTNvBi6js5/cANyQmR/u93sX3SAD8OVdslb/79Z++43AgyfKflC/rSiKothPrJv8FRFbgFFm3tp/Ph143dRu7wHeHBHLdDE/nwS8ITNviogvRsQjM/NTwGnAJ/tjLgReBLy+//c9E9tfHhF/3JfzzQmZrEHOTiHNdbAlRUjraTIgW6XamibLlNhcdd2ltxHNVM5lrs8xmiPv/UAXXFspb+ceWU51UKlKozHvnCPSr0h3oqyonGf3B9yF1pJvaXIzk0ZautLQ10BLzGY+zsB4+6xkYnVaMhlG+0wjwZlInOZCrxEkQH2fzR06B7Zx6zym9nigYHGbbrjGhyXqs9shUbnHDVdw0xO121jUZTkPQMg91gjga7CeNpXjgQuiu8pl4O2ZeVFEvBQgM8/KzKsj4iLgCromOiczr+yP/1ngjyJiE/A54F/3218PvDMiXgJ8AfjxfvtfAj8AfAa4Y2L/oiiKYj8R88bKvzexsmUl7/eo++++cZ6ZihiCLUfCXDMVTewhddLcFq20o/trpjKL5qGwxaBHNN5vbCGb5C6x/SzuWFepYSlo7a1x1AgEtyoWdF2YZxZwse3NNVNRo7pcTyOG03rMVGyBq89Uhs3Mu0KH5eHRBZm6wrMxU9kuThvWl8xHoOUsYzMVc7aR9pxnpjKWjSPJP2R9G4bPVP7pw1+63JaJQK2on+0X5uXRODal89kPiSYnanieWEA4+yHTYHiNlbfqUWYBKe1gdQPygWpsucHtIbMf0UYrj23ws0RX0ka5059GO14DbFqZjR9m/cGU/ex6RuZJ10ioZffSpMOUlerNgAvLsz8D9oPnSeCGJ0LTZFP2I9oI1DiyAcx+69XrTgb9lsombZ/2YqleUY1BtnWyKXbKtY8O9za2AKHLci9jk8hfzaidQwOJtin3p6IoimJh1KBSFEVRLIwaVIqiKIqFscFtKjFj4FUXxZb+aBGFVTq1430814CgYoMYS8TYnXfY0mrXty1xl7k42rHQMOg1ojnPnEeucccdYgQANUaODrNua4bphhFB2tN0Y9XWWy60ZnSV/mHRmM1m1gp2q5Ftrc9pPq1G3eWaNHmWXaMZz2kYhweuzG4mkBK7xJI9MLbNzGsSAbs7Xs4tUQ/0d6HRjS2hVogtbCQ/ICOx94Eb/805ZbRJ7CyNCBart0t0he0tTwGnZipFURTFwqhBpSiKolgYG1r+ikhiyoXX1o+stAK6iT/japiGJMc23H/HIiekOapLlZaPbvihCubLPzI30NZCd3XXNUnNpvOyX2vVrkmMJslZUMZGbDdd4qOJrmY3tYIy2joZvcfSbKtyf1tro4auol6ye9GIWrBqbtuaVF3cf1vSsNwPO95czlca64usx3mSsGGBLzUYJe6ebnLxeOfQRUOQO6T2GuhRDm6sabO1QKuSz360zYJptuYTFpj1IErSVRRFUWwsalApiqIoFkYNKkVRFMXC2NA2FUvGYjlWxi0XR9OYZVtDufVKmX5q7ohLluiqUU9L/CPaq1WpFT4ld8hGDaMh2r6F4bDYW6Bhc8ZmpxE30miEzFDXafHGHsvxYw0t2wiVIvdNw86Y++64EcfH6i6Xae63q3c14j1ZgrOBCeiayeZs14GvsM3wKXYycw+/a7Zzmn0rGq66S0cMS2rHkoUQatjx7JfW3OXlXtg28JAuS1vsOZBn8A7vxxbnq/X716JmKkVRFMXCqEGlKIqiWBgbWv6KnJUObEW8ySCAuwSuDotWq27CjZ3N1ddWpWernnYqufOac7shRVh4cNNHrJ5LJhW1cmaLZGPurmOTIhpukx7V1xIZzTacuZb2Z5utk0V3tVs0R9RmjfhgrqXSbiaDAoyOMJ1OKioRtKMRbtskPZN7UqIe7Gis4Nb7OdQd2xqu0R66gnxYgGRYaURYthD98gxZ385GiP6RRXwQ+dyiPrf616q4PreeoxY1UymKoigWRg0qRVEUxcKoQaUoiqJYGBvapgLMpFh1t8mWa6ron+YSbOFPWmn47I7Iviplt0KqWPZdy1QoOnrLDXVsqXIta6ZssxTD5moL7u6q+vQ8962hUU+zsiIZEVtPjNgbluSGuI4+u18jMhA7B6a7DvFhtayiLdyOZ7YwP17dpC2Ns9wLzRDZ2D7eMesL3g5BMn1y36xhgDRFsRTZ6MeIy7pF+ta+2WhkjyBuYVYs1ba30Vhsi2MLubMGNVMpiqIoFkYNKkVRFMXC2ODyV8D0VFmn7Q1pZjwssu3IVicPU2CaZVrSoGZCrYFykSVhsmi3AHG4rLzVREZyrEy9zc0YGhGFLcKAuCR7BFtIkcVGmyRCgR3b8rGWm2RJvvQqpZ4qo4JGbTYsedXOht4zvktcp+X4oZInQJprrOxnK8Cbq/Q1isOwFeQetcDPo8q0RrWw8NCN1e+3iP41MJKzPWvgUltKkrCR/AC1ZEt1vR4oF999vrn2LoqiKIo1qEGlKIqiWBgbW/4azSag0Tzcrfm4TAtNlrKAbHYe8BzTOsU3z6ZGmRrE0M6jziTDA1+aHDjaPKsbjCXw5KhxHl2ZbUKKtHGrzKG+LCb3mAzSYiTtkbKs3JJktT0Ope3E22lVoj00JTVLHCbHW0KuUSNKpHoMiTSj/biFRmwQjzA5NEzetIgHtBSgYf2r6cUoHlh26RaFISwAK+4FuSoysEWGsMgh4O20JNvWYl0HlYi4DriVrtvuzMxTZJ9Tgd8GVoCbM/Opax0bEScCZwFHAtcBL8jMWyJiBTgHOInuus7LzF9ft4sriqIoZtgfM5WnZebN9kVEHAO8BXhmZl4fEccNOPYc4JWZ+f6IeDHwKuBM4LnA5sz87og4AvhkRLwjM69b5MUURVEUbQ60TeX5wPmZeT1AZm4dcMwJwGX950uA5/SfE9gS3aqvw+km27cstrpFURTFWqz3TCWB90YnnL41M8+e+v4EYCUiLgWOAt6Ymeft4dirgGcD76abnTy43/6ufvuXgCOAn8vMr01XKCLOAM4AWNq8NGNHSPWBbVydeUPq6lWLYNsoVLebvq3hlJ2BLpLqYtlYhW1RTrVKFvFV2mPUcJsc3zl7vLkKW8IluxfgNoixRGf1hdUNW4e1h7hOm0uyRXKOUSP6r7p3iq3BojCIbQ48oZdd5UgSODWjOFiftdXzOwdGR6Bh8jtsdtuq2HNsBfmo9bxIl7Wo3Gr/FM9haCTvMjdlW37QSmBnrtPmxm73qLGiXt31G675LdZ7UHlyZt7Yy1qXRMQ1mXnZxPfLwMnAaXSziw9GxIcy89o1jn0x8KaIOBO4kHvMf0+k6w7fBtwX+LuI+OvM/NxkhfrB6WyATUdtmq+1iqIoijVZV/krM2/s/90KXED3wz/JDcDFmXl7bzu5DDhxrWMz85rMPD0zTwbeAXy2L+v5wEWZuaM/5u+BGceAoiiKYv1Yt5lKRGwBRpl5a//5dOB1U7u9B3hzbwfZBDwJeMNax0bEcZm5NboofK+h8wQDuB54OvCH/THfQ+dV1iZFjjANp5FAytDV6zbNbqx+9yXksp9Ns1uez6KLjURzGVmAu5amJhLBdskNHofJecRnU91qgbGsom4Fw5uhIUWMV2crb4EzTe6ZaxW2uZxbnSzh2aaGa6rJWhrYULa1dCWTcUzykOM1ERloX1w+wrQdObRxe82133a1IKRWT3MTBgiTz6wv2H1rlGn1XLX+of7QWiSrFqRys+xortit51r0r5bk22I95a/jgQv6h2AZeHtmXhQRLwXIzLMy8+qIuAi4gs4icE5mXhkRD7dj+3J/IiJe1n8+Hzi3//y7wLkRcRXdbTg3M69Yx+sriqIopojm28sGYNORm/L+J95/t20auro5UxnYdhYHqbXwauBMpbXIytCZymaZQSzJTKV1HnlLu+uWgTMVM5q2QnHfIXWXN2tbzGnh28Fjfw2fqbRmELJAVWatY0v3aovTGk4GOlOxCYClZWg969b05gxhWYcbz4CmW9Ap1cD60Jip2ONiz4s9b62ZipUp/TOs8q0JgGxblRWI1m7N0PdjqdNm+60x55DGtZvDi1zTTf9w0+W27hAOvEtxURRFcS9iQ4dpycyZsCqpoUYab7wSdsL0R9OIm17KmhxpYOiYRqFL9hZtfoYWLaT1ti9vrStHz/ozTofBgYabcCP6r7oayyzJXD7NvRJwV055Q7NZxagRTVm1cDne3KHt2FHL5dxmJRah2bpHa2KtthLZTQpozX40UrC4cttMWF11u0Jn2Sb3XUIDWd6unebu3jqPXLupBc0Zt9QzpX+mzG5bYVLMbmYzXIvurBHWac3SdNcmNVMpiqIoFkYNKkVRFMXC2NDyFzBrkNQkXa2VyAMTZWk++eGykhnvzG2y6UJr8ohFLlbJoqWpzW4y6U4lF5EnWlKVGaZXRS/S8zSkKq2nRQPQ3O9apK4gt5zuI6uoOR60pCqTz+TclnSs5Q+d4koe22Y7nXrGt/qxuSRL/vOUbYjDCMBYJFJ1bJHoCOZfYY4HXaWGJRizSxzf2XoIZzctyXWODpvDrV91OtPaTV9tlGjRvptivVMzlaIoimJh1KBSFEVRLIyNLX/l7LRYp9Pz0FphPE0rgZSthLZ89HLnLBkPoMHjzCPNkv5EYzaviqDIExo4087T8EYxOXIka19s/ch4p7dHLEvjmdeOSTOtAIrmiTNQdrAAl62V+0g+ec2LbktCGt1D16SIrLVs0SJa3l/D0q+rDGz9Hbx/2rWnBUHVCBStd2qLimFrSoaHA1gSD1KT80biUdpKDGfnT/OGM6m8JYWKdJjzZKajZipFURTFAqlBpSiKolgYNagURVEUC2Nj21Ri1jZhWnIzZpJGpjUxWTTzljYvq7B98bysvBV3xFaVzEtxVc69ersbVULiYi2JnUh1cGnj1YYrpsX0CnOBNXfoRu/WWEh2etOsG+6VI9GtLTqD3XdbMd2K/qsJysytd46Yfmr/MJdzXabvZa6KPWpkURzM7btlb7RkYraqXXyFl2ylect8as15l+03vN3tGdTV77Zfoz3Clj/I8SNzd29cvNk2LaLHWtRMpSiKolgYNagURVEUC2Njy1/MTi0txLW6loJP/aVFzctw3JJRNtlWC2YnU9qGK6bNiS24ny0bXj7au4jmqLcV0xYo0qSdRkBJJNCj6TW66rclo5istWNY1rOR1QdA5MCxFmkRG2z1ecsdWjaatmKSWOsV0vqS7Ka9qxVNX84/FvlqaSRSZqOeKgkOTCth3cvSMgDaTuqObWWa+y6tQKTS7paPXgJCApooz9zYx9uHR4ZY3SayZyvAZ4OaqRRFURQLowaVoiiKYmHUoFIURVEsjI1tUwlmhtWxRGfViKvdN7JJXH01vETDZdQ2muZtempL+xwYYtXcmUM0b4Cd4gJsYWKWD5eoqxIKwhIrgUdNXTI3Y7PnNIRjTUBlidDUvbyRTlh0fHUFNVuHue82ov9aaI+UdtcoxY3uYRF8xxYayApomJiGpt/VqNwNe5JGBbfnyPqs2kkaz7U9b/K8aITm1mu6/VTY74ptaqUeN7um2isHhm4BluS+WRTstaiZSlEURbEw9jioRMTPRsR990dliqIoikObIfLX8cBHIuJ/A/8duDibS8wPMcYwnl4pa9PcRrInTehlUoK5sDZWtNoKdFXflodNncFlLV25O0cUWkv6vSQumuYCa1N0XXkP7n4sUtVIpKbWqnTzLx2JTKd+l9u9zFWJHjwySU9WLJtM1pK/bAU5sjJbc6q3XFPtcJO1rMtZBAjwyNi6en5Y7naAkNDY1pc0z7om1JpDZtPHTfZrRNs2V2GToEyGtQjJ/c6z22RXX6XvRVoBC0/SlZmvAR4BvA34aeDTEfFrEfEdc52pKIqiuNczyKbSz0xu6v92AvcF3hURv7mOdSuKoigOMfYof0XEK4AXAjcD5wCvyswdETECPg38/PpWcR0J8R4xKaKVMNzmkKJVjW063aiSBfJT/cuknc1e5shWe981cBV2U/6SbSbzmdwj3lu6krj7ZmaLyokqGzSKHKbCqKdVa+W/SYfadpZMTIOQNnqIRWzQjGtSnUZkCPNiMsk3jpDAl3dokY0l+cNWkGvgSlDPSuvbFuxUA476WXS7tZ0nuptjlb5ITXbuME8+AEvoZRKjeQw2ErZZwNR5bR1DbCrHAj+amV/Y7USZ44j4obUOjIjrgFvpAgrszMxTZJ9Tgd8GVoCbM/Opax0bEScCZwFHAtcBL8jMW/rvHgu8FTiaThx8QmZuG3CNRVEUxQLY46CSmf9lje+uHnCOp2XmzfZFRBwDvAV4ZmZeHxHHDTj2HOCVmfn+iHgx8CrgzIhYBv4H8FOZ+fGIuB/tqDlFURTFOnCg16k8Hzg/M68HyMytA445Abis/3wJ8Jz+8+nAFZn58b6sr2Y2s3IXRVEU68B6r6hP4L3RGSXemplnT31/ArASEZcCRwFvzMzz9nDsVcCzgXcDzwUePFFWRsTFwAOAP87MGUeCiDgDOAM63XjGddLE9VbYVLF/qPueRWwVF1TAXTHt/GIw0BXP/TcziM7acmM1bPW8ucba6ndLWqa2JBorpkUzt/ejMLdrGu6hdn5p49Gm4Y9M2oksEoLZxxoyOqvD3L6tNS1HFuD92OwStoC84QpudkSNrC0u2s1eaDYMuyZdQS5u6NqP5nDBtagUzdXvtnHoiRotIs+G/VRYRPJWlBCNND6fR/G6DypPzswbe1nrkoi4JjMvm/h+GTgZOA04HPhgRHwoM69d49gXA2+KiDOBC7nHXLUMPBl4AnAH8L6IuDwz3zdZoX5wOhtg5UjpaUVRFMVes67yV2be2P+7FbgAeOLULjfQLaa8vbedXAacuNaxmXlNZp6emScD7wA+O1HWZZl5c2beAfwlcNJ6Xl9RFEWxO+s2U4mILcAoM2/tP58OvG5qt/cAb+6N7JuAJwFvWOvYiDguM7f2Ls2vofMEA7gY+PmIOIJu9vJU4A1rVjKZDR5ni6B3zmGaMRXF4ts18snb6ujx7bP+BiMLmudp3tV10WbUdp0WBBBg2RIcmXxmEkEzQKedaOB5TL6aJ7jf0KCOrdXNgrqia2BBO9jLtICWGmTS3IQbQQRNplMVRhI4tZZmW79ZNtd2a5CGpKar/OU6Q54DyV3VlFztfsRm6wtybCPigvV5jaQgndakxH7nGTR5nuhX+puAS2Uxz/PK+spfxwMX9A23DLw9My+KiJcCZOZZmXl1RFwEXEF3i87JzCsj4uF2bF/uT0TEy/rP5wPn9uV9PSJ+C/gI3aP7l5n5F+t4fUVRFMUUcW8J47U3rGxZyfv/s2N332ihuFsxpAx5S0p589Kw6LBvM5Vm5PthBkWNZdSYqWja5aEzFUt73MrYLG+IFkMqdKbhZdoljS2Glb0Ftxa3GbqiUqpji9gab+tqlN/XmYrMvizlgc7cWimbzZHE4rNZX1iHmYoZ9M1xAOaZqQxzPGjtqzHK9nGmsiSd3hY0WmpnaKRwkLp/+aqbL7d1h3DgXYqLoiiKexEbOklXhETRVc/Sxtv6iiSgGpokx2wSoAmCzKVYw4I03kTt7WPn2KZUcmzDD9VWAKXZZMyd0arZilJsurO6csqxjVAnNvuKzbOPgrr6NnR4q/+y3LcVuZ4dEmJmRyM0kE5+1H94eGhaT6g18M289Vo6NBGaTbgbrr6q7cs1mb1gbCFqWgnGrJ2sf82hYFj4Fo1SbOdp2ErtHmkCOpuVtBL6mVIzp0txzVSKoiiKhVGDSlEURbEwNrT8RTIjkTQT4gihRlObv5qU0HBTFhnFV/4OW7EMLRdLOY8ZPRsJkywplUZ3Nb3GrqeR7ClNpbNdVc1rGCPt/EMjSbcUD9l5Wd7ZtuyUSL9S3F2NCEPbpH/tMDlPNIuWimG9ayzuw54TqhH5WPu87WjSWwuLTGGRpAeu3G85BIhR3m67OS60/G80+q89rxatcI5X/7E5XQyV8/B73GqnFjVTKYqiKBZGDSpFURTFwtjQ8lem+cnLNLnlKWH+47Z8Q7w8WlNK97QYltxIXYMa51oS2cHWO5gMAo166jTb1iHMEczSkp7ZOgJrz4Z3Tm6f1cpyp+ULlzJtXQXumbRDrnOb3DeTyVZH3u6rFhhRvMzGd0obW3BOUMlVZUtbb9VSiweuBbKoA5q0rNt59jTyDOpaDXuG53mlbid13/08DTlQn1dLYDfUkw+X1Mx7ayQSdmvtnXm6zruSsWYqRVEUxcKoQaUoiqJYGDWoFEVRFAtjQ9tUGEEcvvu4GpJAqqX3p9ggTL9ULbiR6Ng0YjsPtkLXiyRNSxc9V20VjdXN6uYsmvlYIwxIdRqxiNTF2+JamZGnIfibFr1ksczmSZ4lbI/ZC90p9rUlcS9fbenoZhvYMcy1VXV9Gu6ytp+u7G72OtlV7EGbJcaYPINdicMiKaxa1GUpb9SylUrbm6suZodrxbCTbWY70nvRWrkv/Vtj9dmSgpZLsdmU5+jzUDOVoiiKYoHUoFIURVEsjA0tfwWz01IN6d5ycTRZTJakmqvuPHm4m2Hyp8/TSjok03TTDdRLuBHqXatk7tRWH2vjRiK02Cn7agBEOU8roKS5gpvEKPe35e5qQT/N5XxV5AVVNxpBBHW1uG1TaaXRj1ftLomkJr8WLWlGA44ODTLZco03hdOqbs+l3d/VVioAi1Cguw47D5CqR0ob6yr9lvw+26AmN4+OkMC3c/gJZzOHhFMzlaIoimJh1KBSFEVRLIwaVIqiKIqFsaFtKsCsD5/ZJVq6sUUKFtdDdVdtZR1VN1YNDztDNH3/hmnuWqdWGA45frRs2q0lMpMIuA3bjaa7Fo3ZXLlXLTETc0QfNvtWI2ORhnSxTZawTVxbWzq6pt81W4eU2XQplu0DzYWezAuIwyRpmdjH1H7ReDg0gLh2eXP1N/dfP8/qNgt5baFsxD286WEtdkBLQCf1HDdtYWaTkf2kzzRtRC036zmomUpRFEWxMGpQKYqiKBbGhpa/cjw71dWVpg0JSGfPm8RlVBIJZcP918q0/NZ26tZ0XuUmW6Qv5xlLMi6AkUl6tqJept4pK8jNXbXb16IRiDxh+7VkS4taMDAne6svxIot7R4WxdpkjGxFsR7oCjo6XGQUWXkP+H0zWUpkMpX9gJFdqCmZ0u6rLZlOpbJhCbVc1pkj45rtZfeyIeNaX1KJ03SpVuJAc+E3N3o7fLm1Sn920zzux1AzlaIoimKB1KBSFEVRLIwNLX8Rs8mddOHrHNO/0WaRqiz/esOLaNkCNZos1cjpbmiAPvM8sWNbUoBdksgWFk1AHMKagkNYgqFtch5NWtYKGDjbHhZcUGWyVnQFaw9NrmbHireSBFrsCrCV7sNWazfzoG2SoI5WpMm4LWnGtCGr513SRq3bJrnjrX/Z/V213O2NB1s952yluklVc/xWjGVFvNLYbSRdxGRTbc7Ws6GRMg6iHPURcV1EfCIiPhYRH23sc2r//VUR8f49HRsRJ0bEB/vv/iwijp4q7yERcVtEvHL9rqwoiqIw9sdM5WmZebN9ERHHAG8BnpmZ10fEcQOOPQd4ZWa+PyJeDLwKOHPi+98C/moxVS+Koijm4UDbVJ4PnJ+Z1wNk5tYBx5wAXNZ/vgR4zq4vIuJHgM8DVy22mkVRFMUQ1numksB7ows7+tbMPHvq+xOAlYi4FDgKeGNmnreHY68Cng28G3gu8GCAiDgS+I/AvwSa0ldEnAGcATDaNGJaRByLON5yTbUEVqo+2n6tltdT2SpsqY+JrAA7Z4V8TfYk2u2oIZmntZPsZ9Lt6LDh7zIpLslL1p624rmVMWnTMPdSlcxbrr62Kl2We2uyJ0mi1rbdWCI0qc8dsqJdr7uB2CWsf7QS2OmDIN1zVTpYMzK2tP1oNPsghNhe1K1e7DnQsCPafbeEfA37h0YfNu9hW2rQsH9oPeV3IeXZaEU093s8n0/xeg8qT87MG3tZ65KIuCYzL5v4fhk4GTgNOBz4YER8KDOvXePYFwNviogzgQuB7X1ZrwXekJm3rRUqvh+czgZY2bIypwd2URRFsRbrOqhk5o39v1sj4gLgidwjXQHcAHw1M28Hbo+Iy4ATgWtbx2bmNcDpABFxAvCDfVlPAn4sIn4TOAYYR8S2zHzzel5jURRFcQ/rNqhExBZglJm39p9PB143tdt7gDdHxDKwiW5geMNax0bEcf1AMwJeA5wFkJn/YuLcrwVuGzSgTM1VdIVu09/VptkDTrJGoSq1aSIi2daSTMS9M2xfkx0akefG20UuEslEFw1LfVpu25awySIMoO6Zw1c3m0Sgt6iVpEvlIgu6aSvARRdqyGxDE4yN5PhmAEWTQkxTE4lPk9rRkLCk7uZ5vCRBFaHxbA1V9Gy/VsQFle4seKM+hIPPb781YYnDGsqL3TeLYDGyZQrNFQnDoya0WM+ZyvHABb0UtQy8PTMvioiXAmTmWZl5dURcBFxBd5nnZOaVEfFwO7Yv9yci4mX95/OBc9fxGoqiKIo5CA0tvkFY2bKS93vUsbtvlPbQcPQ0DPVD46qvx0ylZUvS8PF2noEzGmB1u8Tvsjc3ezu1BaKtbmgGcJ0VzNan2bPF8Dh4ptKeUs1uszD1GndMmGOmYkbooemNoTVTkR3nmKksHWEWY9kk/XBJFmNCSwWQt31xSLBFmnm7W9U19a85CWyyGUBrpmKLEsWArukaWs4hMlORe6wzlea6S0utMHv+mz928+WZeYqVcKBdiouiKIp7ERs7TAsx64Yr7p1Lrei/muhq2LZR421sbNFD7e3Udmvp/YO3yttcy2N0DQ+73fcz99+B1QFPwmQv1vqC15q5yaahM/ZWaI+BbqhD3a5b4XHs7XRosqdWuI2V5dmfAZv9WMTqVtKy2CQbxQ6n7tAWeoWGvdOeDQuFY9Vs/PpZrjt1CbbIw2YToZHgbOhksuEar7MStYua6/HwcEPNhF4NaqZSFEVRLIwaVIqiKIqFsaHlrwiZ6pqbcOt4WxFrhnaL+Now/g82ytvK10ZNVaEYKBFoIiKGR3JVSW67JaVqtIfqRSbTDVsRDzTcKYcZd1sqnV2n1UndWC26QSuisNw3XdVuLtYWrRrvNyapLR1mS+K1SHKbNLJc+9JmWRHfCAyxKkZsU3bGd1rUZ5FhG67LKhNagrKB/i8Ao4HJxDTaRMv91xwnhq4KaOnaco9bzjotaqZSFEVRLIwaVIqiKIqFsaHlr8yclS3m8kyyAHcSRFCTPTXKVKesYYHjdD8gRHYwmU5llIa3k64ZMNlC5t4WJLI1wbZZ+sDlMIy3m8bn62lMHrGTmxcPoNKOLnMZKFv6mgzcc22gF2IriKmtwtbj/Yq0TPPK0sRfst4pW/1Y+qfeD/HA0jZqBQcVmc1iKo6WZjt8HqZFqmypib/Ee6vpDWdSl55ctrX6scl8c049aqZSFEVRLIwaVIqiKIqFUYNKURRFsTA2tE0FmInwmhppt6Hxig3BtNvQ6KqNMtV9b5i+bauTweNaDb3zujoZ17fH5lYr20ZmuxkeMom0BGWybbx9dhu03aTlTFJoa1dzKR4Wl8qEcG0j3MZkNpGcRwjXldlSd7MLtCIuSFysEFuJuqG3wqtZX5KbORa7hPajZqReObf0hbG6bbdW1A9MqLXDbEyNOolLtD3rFum7GW1iAVGKa6ZSFEVRLIwaVIqiKIqFsaHlr2CYFNLMwy3TUndxHOZeCTC2VbJD01Y36qmuh1b3OVYIrxFYfurcA111G26TNvUPza8t7t2tXOe2UaQVUyJVSqDh7jr47MNdl3OHNYjvO3NsQ0fRMPl6/LBoEYDLhLavSDgqF4O7u5pcZHKxyHGjxkNo7vY7Ney/RR1oZvSb3WQSo6Y2aLn1D0yTIcE9LVpEd/6BwUnXoGYqRVEUxcKoQaUoiqJYGDWoFEVRFAtjQ9tUEpH3h6Z7BXUftNSdGnaioWlqmAaLEaF6f8Od0Y0Ds+c2F8dWLh8L7aExROxgq05DtxUXXHW7tm0tLdhSRou9wVyCNZVx43hNKjWwPs1OZ6FFzM6jIXDncGO37jmwH3X7SqgT69sWyLkVFVcj8Fr/mN3P+3HjPLJUYLQsG+eILm11UhuR1qn1EJqdxqok7v9yboC0a2/YO1vUTKUoiqJYGDWoFEVRFAtjQ8tfEbMur7pCuJXLZ/DK8DlkA5un2zR7qGTRYGg002xMk4cmWzdZS11bG683S+Lumior2UrkhlOvJk2T/cwVs5FASiPjDoxIrPVpNbvIjnq8RrtulGk7D+xezQRO5u4qMoolr2ol6Rroxa7uw/O9Pg+L2ozci3HrF9XuW0NKndmvkVxt6DWNxe26hUdJriRdRVEUxQGiBpWiKIpiYdSgUhRFUSyMdbWpRMR1wK3AKrAzM0+RfU4FfhtYAW7OzKeudWxEnAicBRwJXAe8IDNviYh/CbyeLijBduBVmfk3a1YwzRtzeHa7ED9B9Q4dqFm3Tj++azY8rEbWaEYetdMPiyis2RxphfbYB9fWltekZOFr1Gjgfo1AKfYkqPtuq1QLmTHM5qZZHhtGFbvFIfXMwXFS0Hs8tnAwFoqmVU+LziPH75SQO6Nlf9fVU5k7tbWnRdtu/PqpOcpcpK2erUjhcu1DIydnKw2oRVQfuNSgFXrKXJ99+UCb/WGof1pm3mxfRMQxwFuAZ2bm9RFx3IBjzwFemZnvj4gXA68CzgRuBv5VZv5TRDwGuBh44CIvpCiKolibAy1/PR84PzOvB8jMrQOOOQG4rP98CfCc/th/zMx/6rdfBRweEZsXXN+iKIpiDdZ7ppLAe6ML7frWzDx76vsTgJWIuBQ4CnhjZp63h2OvAp4NvBt4LvBgOe9zgP+dmXdNfxERZwBnAIw2jRhPyRYtN1S9OFk9HzJOjw6TRELbPPuVBh+WabbJDtmKoLvZ3AQHrvJvRcuVGXluH5aQS5WzRruPTTqUCAOWrKmV3UiTiUnbxaqU2VLj9MZJlWw3c7tuuXEOlORG0m6rO7xIVqRIix5s7qYaMRo02VPLT3p6v5YHrLpeD0vypbJQIwEdS9IX5FnHnkuTDfHnwBPgSd8+vCEHagRv21PuZUMKtejY8/wmwvoPKk/OzBt7WeuSiLgmMy+b+H4ZOBk4DTgc+GBEfCgzr13j2BcDb4qIM4EL6ewndxMRjwZ+AzjdKtQPTmcDrGxZma+1iqIoijVZV/krM2/s/90KXAA8cWqXG4CLM/P23nZyGXDiWsdm5jWZeXpmngy8A/jsrsIi4kH9vi/MzM9SFEVR7FfWbaYSEVuAUWbe2n8+HXjd1G7vAd4cEct0XltPAt6w1rERcVxmbo3O/eg1dJ5gu4z+fwG8OjP/fmAtySnXlxiJR1fraE02NWwlsSX4gcaUWKfOw+tpgevUM8mSWrVWpdvKcJMDB84Fm7nOVdsZ5j3WklE0mJ7JI+b5ZvITjYCYdhqV/rRIP8/QVdR27c2EWpaX3NyAhiYi8/awYKsmZY5b0sxO01wH5mS3KAyN/mEr3U2+siCirQ6v0RXUlW+4Z6QGqRzogZmNBHb6vA+ULXexnjOV44H/FREfB/4B+IvMvCgiXhoRLwXIzKuBi4Ar+n3OycwrW8f25f5ERFwLXAP8E3Buv/3lwHcC/zkiPtb/TXuTFUVRFOtIDDWe3RtZ2bIp7/eYB+y2TdeeNI63MOiD30jEGAgt491sncby1tZ827dZ0dD0ty0jnb7RyLk19P1wI67F+fIQ6HIae7MFxmZctsmcGUhbqaUHLm9SQ73GXGuke23McKex2VjDb8FnP1KpGM1O3eaJr6brOqSe+z5TsYOHzwBG5qSwWWYqukam0R4yy9NZuNW9dc+te+tzKcdbKH/QWcmqOBV99cqvXW7rDuHAuxQXRVEU9yI2dJRigpm3gBDRu/VuqPYTzTpkxzbexlSzl7cci8TaeGu0WYBdpy/3bgm6Uk9rD03sZNv8/cb08ZA3Xp2VtCIfHzHb7TUac2tWIugsQNxLre5aXqvXmVuwRleYQ4HQlfJ2f+V6LFEUHqFgLP1jZPe3kcDOzyN1t8PtGWjZNW2juMtrxOo5Vp8PddG2yNTQiK5gfVYvs6UMzGKRi9eiZipFURTFwqhBpSiKolgYG1r+CmA0NeHTmWJj6FUpw4zQIlmYOzK4G6sZllWVailVtkrW5KKB6c+7U4mBVZwMPKf7cEO91l19cC0JUsv9V462aooU6W6c7opuco+6Kc8hmai9WY63+2PRHqCxMtvazpxIGi7W9stiq71N6WoGMBzowTvYjb0ZVFH6trWxyGfe3yFMAjdp2HxDGhKjO4cMk8+tf4AHJ7XIAWtRM5WiKIpiYdSgUhRFUSyMGlSKoiiKhbGhbSqQM4Jyjm1lXWvxkenWs4wsiVJLNpYyNfKp2T9k0VZXplVqUJFky71TkwHZfsN2VC0XVCD3hYFSZqM9VHM3w4K40Nr9gZY5SxbMqdu1uOpaSBNotKdtkoWKjQi6tmDX3UjlXrQiCpt3+nZZsGumvcbCYHOjt3usLrjWbI1fv9xmJxf7iRj81GbWwqIh2yLJhu0npX8OXTxpoWig8QwP96wHaqZSFEVRLJAaVIqiKIqFscHlr5iJy5Uyf9XYXaC5oy3Zk63Wbrn0eZL7YXJRa5aqU3JLbmS7NXMwydTdkk3ZSnWTvxrtoeqbuKaqlNG6b+b2bREO7F42ha6hGoEcrzGxGkcPjTEmN64luaoSoqvfJZJBK4rDXcMiFKg03JBcPa7e7Darky40b0UpNkltYPQMj92HR3iWMpfkWW09g7rdLtRUsu2NQode+xrUTKUoiqJYGDWoFEVRFAtjg8tfzMyflzYND5E9NslEtISxLCVuBrNTaWhYgL1xSzYYmktIZTI/l0uCmix9UKEqk7XOs1m8brTdhstffqHWF1pSwDCJQCMpqMtgowBb2T1wZXbLw04i2nuEA9nWcCJibPfDlBmVVhqFav508Zwb6CFnEiFArIgELvu5JNbwwNwhWpu0+6p5mTWjTQz0KjWZrunROkdQ2QY1UymKoigWRg0qRVEUxcKoQaUoiqJYGBvbphIBy9NNIFpjw/1uyZJFpYVntVMf4VUSkXp16daZbbkqdhrRgvtKzW4zjdciBLQWutuKa3VTtkiustK8lYjI9H4VyGXTqvuMquRu+rSmLW40iNmtrCtoamkpr7VS3dppqBtq42aqXcM2WT1b/cPa0+xmmlray7Q0v+qmbBGazW27lXxK3ILNndn6YSPXHLHZlirIjlKnZqRwjZCgIadlW8v+aoav+ZbU10ylKIqiWBg1qBRFURQLY4PLXyOWlo/efZMElMzlFT18eUVW1EtGLsvDfdQDjtQyt6189+y2r9wxs218x0dmz73tTi1TV4ab56Ad21gxPTRLmE3dLRBeU2az86vrs8gtDdlyVfKNxyY7t2xr+NBqki9joKtvS0YZ75T+aY2n+eQbrum2WYI3eu73lhwofc6u3dqzFbRzmOe0rmo3F+loBdgUN3iTezWCxByrz1UKtWfVli4AaZE29ESyqRGwVJ/XVpSABjVTKYqiKBZGDSpFURTFwqhBpSiKolgY62pTiYjrgFvpHCR3ZuYpss+pwG8DK8DNmfnUtY6NiBOBs4AjgeuAF2TmLf13vwC8pD/m32fmxWvVb2nlcI487jG7bTvsiAfMXscxT9Ljx9tnx+Q7/ulTM9t23vwXM9u+dckyAcG3/chjZrZ9+uvPmNl2/Tt/cWbb6l1XaJkhYWIs4uxOzaykRbp4q5GTxX14+2x9msmeLDSI6MGW3CgbyY1GYj8ZmkysZfwJC0djmrvp+JaYTc+Ct/FAF+lmuCH5FTC3bwsZPbawHsBIXHhHA+00jdtGSv8ci5uyRRk2u5Um82owuHs0EsOpnch+fdV9txUZ2xpv4Dyh9bxZ/3KTcpP9Yah/WmbebF9ExDHAW4BnZub1EXHcgGPPAV6Zme+PiBcDrwLOjIjvAp4HPBr4NuCvI+KEtPj0RVEUxbpwoOWv5wPnZ+b1AJm5dcAxJwCX9Z8vAZ7Tf3428MeZeVdmfh74DPDEBde3KIqiWIP1nqkk8N6ISOCtmXn21PcnACsRcSlwFPDGzDxvD8deRTeAvBt4LvDgfvsDgQ9NlH1Dv203IuIM4Iz+v7dd/Y4/ntWrZvidPe8yJ5/r/rk/sPtM7E9fufBz7Udmr+fQ5t52PXDvu6a6ngPDt7e+WO9B5cmZeWMva10SEddk5mUT3y8DJwOnAYcDH4yID2XmtWsc+2LgTRFxJnAhsH2eCvWD0/TgdkCIiI+anelQpa7n4Ofedk11PQcf6yp/ZeaN/b9bgQuYlaNuAC7OzNt728llwIlrHZuZ12Tm6Zl5MvAO4LN9WTdyz6wF4EH9tqIoimI/sW6DSkRsiYijdn0GTgeunNrtPcCTI2I5ugiLTwKuXuvYXcb8iBgBr6HzBINu1vK8iNgcEQ8DHgH8w3pdX1EURTHLespfxwMX9KEIloG3Z+ZFEfFSgMw8KzOvjoiLgCvogkWck5lXRsTD7di+3J+IiJf1n88Hzu3Luyoi3gl8ki7O6MsOAc+vg0KGWyB1PQc/97Zrqus5yAhNG1oURVEUe8GBdikuiqIo7kXUoFIURVEsjBpU9pKIOCYi3hUR10TE1RHxvRFxbERcEhGf7v+9b7/voyLigxFxV0S8cqKMB0fE30bEJyPiqoh4xcR3rbIiIt4UEZ+JiCsi4qSD5XomylqKiH+MiD+f2PawiPhwX+8/ieiCzfeOFX/Sb/9wRDz0YLoeK6ffvl/vz4Kv6ef6/nZlRLwjIg7rtx/M9+gFfXt+IiI+EF24pl3lPDMiPtXX79UT2w+564mD6Ddhr8nM+tuLP+APgH/Tf94EHAP8JvDqfturgd/oPx8HPAH4VboQM7vK+FbgpP7zUcC1wHf1/2+V9QPAX9FF3/oe4MMHy/VMlPUfgLcDfz6x7Z3A8/rPZwE/03/+d8BZ/efnAX9yMF2PlXMg7s8C+9wDgc8Dh0/cl58+BO7R9wH37T8/a1e70mXX+Szw8L6Mj3PPM3QoXs9B85uw1+1wIE9+qP4B9+kfzJja/ingWyc6x6emvn8t8iM88f17gH+5VlnAW4GfsHMeDNdDtz7ofcDT6QeVvrPfDCz3//9euvVJABcD39t/Xu73i4Phelrl7O/7s+BreiDwReDYvr3/nM5l/5C4R/32+wI3Ttez//8v9H+H5PXIdwfkN2Ff/kr+2jseBnwFODc6meec6NbTHJ+ZX+r3uYnOrXoQ/RT88cCH+02tsnb9KOxCw9HMySKv57eBn2f3fIL3A76ReXcOuck63309/fff7PffFxZ1Pa1yWKOs9bg/a9VlrmvKblHxfwOuB74EfDMz38uhdY9eQvdmvlvdpup9qF7P3Rzg34S9pgaVvWMZOAn4vcx8PHA73XT0brJ7ZcghhUXEkcD/BP6f7MP4721Ze8lCricifgjYmpmXr1dFB7Ko+7PHcuYoa19Z1D26L13svIfRRfPeEhE/uS41Xpu9up6IeBrdj/B/3E/1HMpCr+cg+E3Ya2pQ2TtuAG7IzF1vEO+i61BfjohvBej/3WPU5YhYoes8f5SZ50981SprPcLRLOp6/jnww9Hlwvlj4OkR8T+ArwLHRMSuxbaTdb77evrv79PvfzBcT6sc1ihrvcIFLeqangF8PjO/kpk76BYQfx+HwD2KiMfSpb54dmbuOn+rvQ/V6zlYfhP2mhpU9oLMvAn4YkQ8st90Gt1K/guBF/XbXkSnhzaJiADeBlydmb819XWrrAuBF/YeH99DJ198iX1gUdeTmb+QmQ/KzIfSGUD/JjN/sn+r+lvgxxrXs+scP9Yfs09vYAu8nlY50/Ve1/uzh7rMdU10stf3RMQRff87ja7/HdT3KCIeQjcA/lR2AWd38RHgEb2n1ya6fnfhoXo9B8tvwj5xoIw5h/of8Djgo3QhZt5NZ2y7H52R+tPAXwPH9vt+C92bzC3AN/rPRwNPppvCXgF8rP/7gf6YVlkB/C6dx8sngFMOluuZKu9Udvf+ejhdLLbPAH8KbO63H9b//zP99w8/mK7HyjkQ92fB1/RLwDV08fT+cOJeHMz36Bzg69zznHx0opwfoPOS+izwnw6RPqfXw0H0m7C3fxWmpSiKolgYJX8VRVEUC6MGlaIoimJh1KBSFEVRLIwaVIqiKIqFUYNKURRFsTBqUCmKoigWRg0qRVEUxcKoQaUoDiIi4gl9TozDImJLn1PjMQe6XkUxlFr8WBQHGRHxK3Srvg+niyf16we4SkUxmBpUiuIgo49h9RFgG/B9mbl6gKtUFIMp+asoDj7uBxxJl/nvsANcl6KYi5qpFMVBRkRcSJc64GF0GfxefoCrVBSDWd7zLkVR7C8i4oXAjsx8e0QsAR+IiKdn5t8c6LoVxRBqplIURVEsjLKpFEVRFAujBpWiKIpiYdSgUhRFUSyMGlSKoiiKhVGDSlEURbEwalApiqIoFkYNKkVRFMXC+D8WY4UJpfOO2wAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "(da.sel(band=[\"B04\",\"B03\",\"B02\"])[1]/2000).clip(0,1).plot.imshow()" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "cubo", + "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.9.16" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +}