{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Area of Interest Definition and Product Downloads\n",
"\n",
"This notebook defines three forest fires in northern Germany in 2018 and 2019 for based on newspaper coverage during the fire and information published later by governmental bodies.\n",
"\n",
"After introducing relevant articles, it uses the OpenStreetMap API to query products before and after the fire and thereby builds the foundation for further analysis.\n",
"\n",
"## Setup"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from datetime import date, timedelta\n",
"import os\n",
"from pathlib import Path\n",
"from sentinel_helpers import search_osm, plot_product_extent, geodataframe_on_map\n",
"from sentinelsat import SentinelAPI\n",
"from tqdm.notebook import tqdm\n",
"\n",
"api = SentinelAPI(os.getenv('SCIHUB_USERNAME'), os.getenv('SCIHUB_PASSWORD'))\n",
"\n",
"out_dir = Path('resources/forest_fires')\n",
"out_dir.mkdir(exist_ok=True, parents=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Treuenbrietzen 2018\n",
"\n",
"_Deutsche Welle_, a German public broadcaster, [published an article](https://web.archive.org/web/20210222122952/https://www.dw.com/en/forest-fires-near-berlin-could-burn-for-days/a-45203042) covering fires \"some 50 kilometers outside the German capital\" on August 24 2018. The article was published while the fires were burning and noted that they \"began on Thursday\" (August 23 2018). It contains a map designating the position of three towns and villages which had to be evacuated and takes note that train traffic was disturbed as a direct consequence of the ongoing fire.\n",
"\n",
"### Area of Interest Definition\n",
"\n",
"Using information from the article, the evacuated places positions can be queried from OpenStreetMap:"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
Make this Notebook Trusted to load map: File -> Trust Notebook
"
],
"text/plain": [
""
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"import geopandas as gpd\n",
"import pandas as pd\n",
"\n",
"frohnsdorf = search_osm('Frohnsdorf, Brandenburg').head(n=1)\n",
"klausdorf = search_osm('Klausdorf, Brandenburg').head(n=1)\n",
"tiefenbrunnen = search_osm('Tiefenbrunnen, Brandenburg').head(n=1)\n",
"\n",
"evacuated_places = gpd.GeoDataFrame(pd.concat([frohnsdorf, klausdorf, tiefenbrunnen]))\n",
"geodataframe_on_map(evacuated_places)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Comparing these search results with the map given in the article, they appear to be correct.\n",
"They are saved for later use."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"evacuated_places.to_file(out_dir / 'treuenbrietzen_2018.json', driver='GeoJSON')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Product Subselection\n",
"\n",
"To calculate the dNBR, the NBR has to be calculated before and after the actual fire.\n",
"It estimates damaged plan material, so it is best calculated before processes such as regrowth and plant activity have taken place. Similarly, the pre-fire NBR should be as close to the start of the fire as possible.\n",
"\n",
"Because the exact location of the fire is not mentioned in the article, a bounding box containing the three evacuated places is constructed and a geometric buffer around it calculated. These buffer calculations take place in EPSG:25833, which is a metric coordinate reference system, as opposed to EPSG:4326, which uses angular units."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/opt/conda/lib/python3.8/site-packages/pyproj/crs/crs.py:53: FutureWarning: '+init=:' syntax is deprecated. ':' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6\n",
" return _prepare_from_string(\" \".join(pjargs))\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAGTCAYAAAAxyUkTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAA3L0lEQVR4nO2de7wcVZXvv0dAB0UNjwARQRlfkTsohhnkjuMUMb4yowhqLxXBBMEgjKCRDz6iueAgGRQxihrknUhEWXMVxQeKaKYY7wwqREYQgswIyiMCjonjAx/IuX9UBU53Tp/q6n261t6n1vfz6c/p03V+vdZv716r61RX7xobHx/HcRzHSY9HWCfgOI7jDIc3cMdxnETxBu44jpMo3sAdx3ESxRu44zhOongDdxzHSZSoGnie50tcG39s95yG1jJ2ilrL2MNqo2rgwPWuTSK2ldYydopay9gpai1jD6WNrYHv7dokYrvnNLSWsVPUWsYeShtbA7/ftUnEds9paC1jp6i1jD2UNrYGvtm1ScS20lrGTlFrGTtFrWXsobSxNfC5rk0itntOQ2sZO0WtZeyhtLE18HWuTSK2e05Daxk7Ra1l7KG0sTXww1ybRGz3nIbWMnaKWsvYQ2lja+BnujaJ2O45Da1l7BS1lrGH0sbWwM9wbRKx3XMaWsvYKWotYw+lHfMLOtTjnJU3TTlgxyzdZ6ypXBzHaTdR7YHneb62Ce05K28an+o2bA51CfEbqk9Raxk7Ra1l7BS1lrGH1bZyDzxkL9r3wB3HiYVW7oHHQtv2FkK1lrFT1FrGTlFrGXtke+AisifwKWB34EHgXFX9qIicArwJuK/802Wq+tVhkmga3wN3HGcmMMge+APAiar6TOBA4B9EZJ9y20pV3a+8BTfvPM9XWmitCM3Zarws5ynFvH280tBaxh5Wu23VH6jqRmBjef9XInIzsMcwwQZghZHWitCcrcbLcp5SzNvHKw2tZeyhtLU+xBSRJwNXA38BvB1YDPwPcC3FXvqmSTRLgCUAc+fOvXnBggXnAbOA7YHbKPbqLwOOA34JzM2ybHGe5xdnWXZEnuergaOBU4FVwKHANRTLL95PsQjMXOApwG8pTog/I8uyY/M8X5tl2eETfq4EVmxYP/veqXzOnXffrsCyLMuW9j5H1SGUufPue26Pp+XA+X08XQksmsLTOopvaE3qCfgBMIdi8o8ErgDmUb7hltvWAwuBiyZ6KrXPKp/rbOAk4ETgEmA+sMF6nvp4Ogp4/2Seep5jMk/PKZ9vMk9d87Rh/ezDK+b5FQnM03LgmizL/nKIeVoHfAJ4ZczzNImnB4F7pvDUd57Kn/9aeq41T6V2fZZl80YxT1mW3cEkDNzARWQHIAdOU9XPi8huwM+B8TLIHFV940BP1oc8z/fNsuyGUWtjOQYe4jdUn6K2ydixzHMq4zUTtJaxh9UOdBaKiGwHfA74tKp+HkBV71HVP6nqg8B5wAF1g0/CPCOtFaE5W42X5Tyl+Bpp43ilqLWMPZS2soGLyBhwAXCzqn54wuNzJvzZocCNwyTQw8bqPxmJ1orQnK3Gy3KeUnyNtHG8UtRaxh5KW/khJvA84AjgBhG5vnxsGfA6EdmP4hDK7cAxwyTgOI7jDMcgZ6F8G5jseN8ozvmeU/0nI9FaEZqz1XhZzlOKr5E2jleKWsvYQ2kH2QNvkvVGWitCc25kvLb+QG82G9bftHrLbzW/vJSE52nGMmer2ClqLWMPpY3qq/QUp+dYaK0IzTnF8XLPzWktY6eotYw9lDa2Bn6RkdaK0JxTHC/33JzWMnaKWsvYQ2lja+DLjLRWhOac4ni55+a0lrFT1FrGHkrry8lOgi9m1Y173pqZ6NlJj6j2wH052eb0VuPlnpvTWsZOUWsZ2y/oUAPfA6+He96amejZSQ/fAzck1b2FEJr0PBMunZfqayRFrWVs3wOvge+B1yNVzz7Pzkwntj3wsy20VoTmnOJ4WXpuiunc80/1NZKi1jL2sNrYvol5kpHWitCcB9b3+Tblm7f81uAeZWOeZwiW45WEtvu1Hfy6TsLzFqLaA6dYAN1Ca0Vozu555mM5XilqQ0nKc2wN/BIjrRWhObvnmY/leKWoDSUpz7E18PlGWitCc3bPMx/L8UpRG0pSnmM7Br7BSGtFaM7ueeZjOV4pamvR57Ohc7b8VvMYeuOeY9sDn2WktWKWsd6CWcb61JhlqE9Ra8msprWxNfDtjbRWhObsnmc+luOVotaSxj3H1sBvM9JaEZqze575WI5XilpLGvccWwM/0EhrRWjO7nnmYzleKWotadxzbA38MiOtFaE5u+eZj+V4pai1pHHPsTXw44y0VoTm7J5nPpbjlaLWksY9x9bAlxtprQjN2T3PfCzHK0WtJY17jq2Bn2+ktSI0Z/c887EcrxS1ljTuOaoGnmXZYgutFaE5u+eZj+V4pai1xMJzVA08z/OLLbRWhObsnmc+luOVotYSC89RNfAsy46w0FoRmrN7nvlYjleKWkssPEfVwPM8X22htSI0Z/c887EcrxS1llh4jqqBA0cbaa0Izdk9z3wsxytFrSWNe46tgZ9qpLUiNGf3PPOxHK8UtZY07jm2Br7KSGtFaM7ueeZjOV4pai1p3HNsDfxQI60VoTm755mP5XilqLWkcc+VF3QQkT2BTwG7Aw8C56rqR0VkJ+BS4MnA7YCo6qZhkpjANUZaK0Jzds8zH8vxSlFrSeOeB9kDfwA4UVWfSbFi1j+IyD7Au4BvqurTgG+Wv4eyt5HWitCc3fPMx3K8UtRa0rjnygauqhtVdX15/1fAzcAewCuANeWfrQEOGSaBHu430loRmrN7nvlYjleKWksa91zrmpgi8mTgOcB3gN1UdSMUTV5Edu2jWQIsAZg7d+7NCxYsOI/i8kHbUyxifiDFUorHAVfmeb46y7LFeZ5fnGXZEeX5kUdTfEq7iuJY0TUU71j3A5uBucCv8zw/GTgTOCPLsmPzPF+bZdnhE36uBFbA7Cl95nk+G1iWZdnS3ueoGqM8zw/o8bQcOL+PJ8nz/PtTeFoHHNbPE3Duw544ErgCmAdsLNOZA6wHFg7geQfgRIqrY8+nuEbfLGD7Km2p3+JpGudpUk/PzfP8tsITF9FnnvI8P3sAzyfQZ54G0B5Mn3kaQLsvfeep2xPFHE/wxEn0naet6mk5cCxw+RD1tA44JM/zbzHyedrK0yPL8e3nqauehp2nLMuOnXKiCv2iQeap7A1vKcd62ucpy7I7JstvbHx8fLLHt0JEdgBy4DRV/byIbFbVWRO2b1LVHQd6sj7keb4ky7JzR63d+kKm3Ux1IdMQbS8hfuvq3fPWjMpzquOVqnaq8a4a61jmaljtQGehiMh2wOeAT6vq58uH7xGROeX2OcC9dYNPwjojrRWhObvnmY/leKWotaRxz5UNXETGgAuAm1X1wxM2XQ4sKu8vAr44TAI9HGaktSI0Z/c887EcrxS1ljTueZBj4M8DjgBuEJHry8eWAacDKiJHAT8FOsMk0MOZRlorQnN2zzOQ7n/rZ7Nh/U2nbPmtzr/0JVY11bZaBgPPlQ1cVb8N9HvRLBgm6BScQXEwv2mtFaE5u2enCquaalstg4HnqL6JOcinwqPQWhGas3t2qrCqqbbVMth4jqqBl6dNNa61IjRn9+xUYVVTbatlsPEcVQMf5DzrUWitCM3ZPTtVWNVU22oZbDxH1cDb9q7dxr3RNnq2JMW96FTn2PfAW/au3ca90TZ6tiTFvehU59j3wIuv5jautSI0Z/fsVGFVU22rZbDxHFUDp1hXwUJrRWjO7tmpwqqm2lbLYOA5tgZ+pJHWitCc3bNThVVNta2WwcBzbA38CiOtFaE5u2enCquaalstg4Hn2Br4PCOtFaE5u2enCquaalstg4Hn2Br4xuo/GYnWitCc3bNThVVNta2WwcBzbA3ccRzHGZDYGvgcI60VoTm7Z6cKq5pqWy2DgefYGvh6I60VoTm7Z6cKq5pqWy2DgefYGvhCI60VoTm7Z6cKq5pqWy2DgefYGvhFRlorQnN2z04VVjXVtloGA8+xNfBlRlorQnN2z04VVjXVtloGA89RNfAsy5ZaaK0Izdk9O1VY1VTbahlsPEfVwNu2BGUbl1Zto2dLUlwSNtU59uVkW7YEZRuXVm2jZ0tSXBI21Tn25WRb9q7dxr3RNnq2JMW96FTn2PfAW/au3ca90TZ6tiTFvehU59j3wPP8bAutFaE5u2enCquaalstg43nqBo4cJKR1orQnN2zU4VVTbWtlsHAc2wN/EQjrRWhObtnpwqrmmpbLYOB59ga+CVGWitCc3bPThVWNdW2WgYDz7E18PlGWitCc3bPThVWNdW2WgYDz7E18A1GWitCc3bPThVWNdW2WgYDz7E18FlGWitmGestmGWsbxuzWqa1ZFbT2tga+PZGWitCc3bPThVWNdW2WgYDz7E18NuMtFaE5uyenSqsaqpttQwGnret+gMRuRB4GXCvqv5F+dgpwJuA+8o/W6aqXx0mgR4OBL5roLUiNGf37FRhVVNtq2Uw8FzZwIHVwMeBT/U8vlJVP1Q3YAWXGWmtCM3ZPTtVWNVU22oZDDxXHkJR1auBXwzz5ENwnJHWitCc3bNThVVNta2WwcDzIHvg/XiLiLwBuBY4UVU3TfZHIrIEWAIwd+7cmxcsWHAexSeu21Mc9zmQ4t3nOOCUPM9XZ1m2OM/zi7MsOyLP89XA0cCpwCrgUOAaYG/gfmAzMBdYk+f5ycCZwBlZlh2b5/naLMsOn/BzJbACZk9pLM/z2cCyLMuW9j5H1aDkeX5Aj6flwPl9PI3leb7nFJ7WAYf18wQsftgTRwJXAPOAjWU6cygulrpwAM87UHwb7BKKc1I38NA8Ta0t9Vs8TeM8TerpR3me71t44iL6zFOxtkSl5xPoM08DaA+mzzwNoN2XvvPERVPpp56nreppOfCEUle3ntYBfyjjjXieOKnH05Xl+Pbz1FVPFWPdd56yLDt2Km2pX0Tfeer2BOxVaibzFDRPWZbdMVl+Y+Pj41UeEJEnA1+ecAx8N+DnwHgZYI6qvrHyiSrY0rxHrT1n5U1Tmj5m6T5jo9D2EuK3rt49b82oPIeO11T6OmMNzdWUpTZkvGJ5bQ+rHWoPXFXv2XJfRM4DvjzM8/QSUtghWitCc3bPThVWNdW2WgYbz0OdRigicyb8eihw4zDP00vxr2vzWitCc3bPThVWNdW2WgYbz4OcRvgZ4CBgFxG5EzgZOEhE9qM4hHI7cMwwwXupOp41Kq0VoTm7Z6cKq5pqWy2DjedBzkJ5narOUdXtVPWJqnqBqh6hqvuq6rNU9WBV3Vj1PINQHrhvXGtFaM7u2anCqqbaVstg4zm2b2IebaS1IjRn9+xUYVVTbatlMPAcWwM/1UhrRWjO7tmpwqqm2lbLYOA5tga+ykhrRWjO7tmpwqqm2lbLYOA5tgZ+qJHWitCc3bNThVVNta2WwcBzbA38GiOtFaE5u2enCquaalstg4Hn2Br43kZaK0Jzds9OFVY11bZaBgPPsTXw+420VoTm7J6dKqxqqm21DAaeY2vgm420Vmw21luw2VjfNja3TGvJ5qa1sTXwuUZaK0Jzds9OFVY11bZaBgPPsTXwdUZaK0Jzds9OFVY11bZaBgPPsTXww4y0VoTm7J6dKqxqqm21DAaeY2vgZxpprQjN2T07VVjVVNtqGQw8x9bAzzDSWhGas3t2qrCqqbbVMhh4jqqBD3KJo1ForQjN2T07VVjVVNtqGWw8R9XAy+s8Nq61IjRn9+xUYVVTbatlsPEcVQMf5KLBo9BaEZqze3aqsKqpttUy2HiOqoG37V27jXujbfRsSYp70anOse+Bt+xdu417o230bEmKe9GpzrHvgef5SgutFaE5u2enCquaalstg43nqBo4sMJIa0Vozu7ZqcKqptpWy2DgObYGfqSR1orQnN2zU4VVTbWtlsHAc2wN/AojrRWhObtnpwqrmmpbLYOB59ga+DwjrRWhObtnpwqrmmpbLYOB59ga+EYjrRWhObtnpwqrmmpbLYOB59gauOM4jjMgsTXwOUZaK0Jzds9OFVY11bZaBgPPsTXw9UZaK0Jzds9OFVY11bZaBgPPsTXwhUZaK0Jzds9OFVY11bZaBgPPsTXwi4y0VoTm7J6dKqxqqm21DAaeY2vgy4y0VoTm7J6dKqxqqm21DAaeo2rgWZYttdBaEZqze3aqsKqpttUy2HiubOAicqGI3CsiN054bCcR+YaI3Fr+3HGY4L20bQnKNi6t2kbPlqS4JGyqcxzrcrKrgZf2PPYu4Juq+jTgm+XvwbRtCco2Lq3aRs+WpLgkbKpzHOVysqp6NfCLnodfAawp768BDhkmeC9te9du495oGz1bkuJedKpzbOF52yHj7aaqGwFUdaOI7NrvD0VkCbAEYO7cuTcvWLDgPGAWsD1wG3AgcBlwHLA4z/PVWZYtzvP84izLjsjzfDVwNHAqsAo4FLgG2Bu4H9gMzAXel+f5ycCZwBlZlh2b5/naLMsOn/BzJbACZk9pLs/z2cCyLMuW9j5H1cDkeX5Aj6flwPl9PN2R5/meU3haBxw2mCeOpFgQZx4Pfy13DsX5pQsH8LwDcCJwCTAf2MBD8zS1ttRv8TSN8zSpp2/keb5v4YmL6DNPeZ6fPYDnE+gzTwNoD6bPPA2g3Ze+88RFU+mnnqet6mk58ECpq1tP64Bby3gjnidO6vF0bjm+/Tx11VPFWPedp0EuJJzn+SL6zlO3pwmayTwFzVOWZXdMlt/Y+Ph4lQdE5MnAl1X1L8rfN6vqrAnbN6lq8HHwPM/PHvbqzHW056y8aUrTxyzdZ2wU2l5C/NbVu+etGZXn0PGaSl9nrKG5mrLUhoxXLK/tYbXDnoVyj4jMASh/3jvk8/RykpHWitCc3bNThVVNta2WwcDzsA38cmBReX8R8MUhn6eXE420VoTm7J6dKqxqqm21DAaeK4+Bi8hngIOAXUTkTuBk4HRAReQo4KdAZ5jgk3CJkdaK0Jzds1OFVU21rZbBwHNlA1fV1/XZtGCYgBXMB2410FoRmrN7dqqwqqm21TIYeI7qm5gUn9RaaK0Izdk9O1VY1VTbahkMPMfWwGcZaa2YZay3YJaxvm3MapnWkllNa2Nr4Nsbaa0Izdk9O1VY1VTbahkMPMfWwG8z0loRmrN7dqqwqqm21TIYeI6tgR9opLUiNGf37FRhVVNtq2Uw8BxbA7/MSGtFaM7u2anCqqbaVstg4Dm2Bn6ckdaK0Jzds1OFVU21rZbBwHNsDXy5kdaK0Jzds1OFVU21rZbBwHNsDfx8I60VoTm7Z6cKq5pqWy2DgeeoGniWZYsttFaE5uyenSqsaqpttQw2nqNq4MU6zM1rrQjN2T07VVjVVNtqGWw8R9XAqxZnH5XWitCc3bNThVVNta2WwcZzVA28vApF41orQnN2z04VVjXVtloGG89RNXCKSwhZaK0Izdk9O1VY1VTbahkMPMfWwE810loRmrN7dqqwqqm21TIYeI6tga8y0loRmrN7dqqwqqm21TIYeI6tgR9qpLUiNGf37FRhVVNtq2Uw8BxbA7/GSGtFaM7u2anCqqbaVstg4Dm2Br63kdaK0Jzds1OFVU21rZbBwHNsDfx+I60VoTm7Z6cKq5pqWy2DgefYGvhmI60Vm431Fmw21reNzS3TWrK5aW1sDXyukdaK0Jzds1OFVU21rZbBwHNsDXydkdaK0Jzds1OFVU21rZbBwHNsDfwwI60VoTm7Z6cKq5pqWy2DgefYGviZRlorQnN2z04VVjXVtloGA8+xNfAzjLRWhObsnp0qrGqqbbUMBp6jauBZlh1robUiNGf37FRhVVNtq2Ww8RxVA8/zfK2F1orQnN2zU4VVTbWtlsHGc1QNPMuywy20VoTm7J6dKqxqqm21DDaeo2rgbXvXbuPeaBs9W5LiXnSqc+x74C17127j3mgbPVuS4l50qnOc3B64iNwuIjeIyPUicm3IcwHkeb7SQmtFaM7u2anCqqbaVstg43nbYQNOYL6q/nwangdghZHWitCc3bNThVVNta2WwcBzVIdQgCONtFaE5uyenSqsaqpttQwGnkP3wMeBK0VkHDhHVc/t/QMRWQIsAZg7d+7NCxYsOA+YBWwP3AYcCFwGHAdcmuf56izLFud5fnGWZUeUV2s+muKacasorlxxDcX6ufdTrOI1F/hhnucnU3yj6Ywsy47N83xtlmWHT/i5ElgBs6c0lef5bGBZlmVLe5+jakDyPD+gx9Ny4Pw+np6e5/meU3haR/EV20k9AR942BNHAlcA84CNZTpzgPXAwgE87wCcCFwCzAc28NA8Ta0t9Vs8TeM8TeppuzzP9y08cRF95inP87MH8HwCfeZpAO3B9JmnAbT70neeuGgq/dTztFU9LQf+utTVrad1wJPKeCOeJ07q8XRvOb79PHXVU8VY952nQc69zvN8EX3nqdsT8LfAB/t4CpqnLMvumCy/sfHx8SoPfRGRJ6jq3SKyK/AN4HhVvXrY58vzfFGWZWtGrT1n5U1Tmj5m6T5jo9D2EuK3rt49b82oPIeO11T6OmMNzdWUpTZkvGJ5bQ+rDdoDV9W7y5/3ishlwAHA0A2ch9/pmtY2RvcLZjYb1t+0estvdYuTRDz3EJpzip4tsaqpGV/Lk9C456GPgYvIY0TksVvuAy8Gbhz2+RzHcZx6hOyB7wZcJiJbnucSVf1aYD5zjLSpkqLn0JxT9GyJVU21sZYb9zx0A1fVHwPPHlbfh/VG2iTY+njdw4dghjj8YkXoPM34eZ5mrGqqjbXcuOfYTiNcaKR1miN0nnye62FVU22s5cY9x9bALzLSOs0ROk8+z/Wwqqk21nLjnmNr4MuMtE5zhM6Tz3M9rGqqjbXcuOeoGniWZUsttE5zhM6Tz3M9rGqqjbVs4TmqBt7GJSjbhi8n2ywpLgmb6hzXyfuclTeNT3Ub9HmiauBtXIKybfhyss2S4pKwqc6xRd5RNfA2vmu3Dd8Db5YU96JTnWOLvKNq4G18124bvgfeLCnuRac6x74HXqzi1bjWaY7QefJ5rodVTbWxli3yjqqBUyzBaKF1miN0nnye62FVU22s5cbzjq2Bn2ikdZojdJ58nuthVVNtrOXG846tgV9ipHWaI3SefJ7rYVVTbazlxvOOrYHPN9I6zRE6Tz7P9bCqqTbWcuN5x9bANxhpneYInSef53pY1VQba7nxvGNr4LOMtE5zzDLWt41ZLdNaMqvpgKEXNZ5utjfSOs0ROk8+z/WwqqkZX8t91ud/6Lcm1uiPbQ/8NiOt0xyh8+TzXA+rmvJaboDYGviBRlqnOULnyee5HlY15bXcALE18MuMtE5zhM6Tz3M9rGrKa7kBYmvgxxlpneYInSef53pY1ZTXcgPE1sCXG2md5gidJ5/neljVlNdyA8TWwM830jrNETpPPs/1sKopr+UGiKqBZ1m22ELrNEfoPPk818OqpryWmyGqBp7n+cUWWqc5QufJ57keVjXltdwMUX2RJ8uyIyy0TnOEzpPPcz3qjFefL6Y8dJGCOl9M8Vpuhtj2wFdbaJ3mCJ0nn+d6WI2X13IzRNXAgaONtE5zhM6Tz3M9rMbLa7kBojqEApwKvNtA6zRH6Dz5PE/BVOtzjHJtjmleF8TneEBi2wNfZaR1miN0nnyeZz4+xwMS2x74ocBZBlqnOULnyed5hhHDqn6pYtrA+0zcR6F60qbSDqJ3zLimzh+HvEYcZ6YT1MBF5KXAR4FtgPNV9fRpycqZUXQ3Yd+7cpzpYuhj4CKyDfAJYCGwD/A6EdlnuhJzHMdxpibkQ8wDgP9U1R+r6h+AzwKvmJ60HMdxnErGx8eHunU6nVd3Op3zJ/x+RKfT+fgkf7ek0+lcW96WVDznlNtdG0ds95yGNtW8fbwGv4XsgU927HK89wFVPVdV/7K8nVvxnEsC8mmb1jK2e05Daxk7Ra1l7KG0IQ38TmDPCb8/Ebg74Pkcx3GcGoSchfI94GkisjdwF/Ba4LBpycpxHMepZOg9cFV9AHgL8HXg5uIh/WFgPlWHWFwbR2z3nIbWMnaKWsvYQ2nHxse3OmztOI7jJEBsa6E4juM4A+IN3HEcJ1G8gTuO4ySKyWJWIrItcBTFynJPoDh//G7gi8AFqvpH19rHds9paFPN28erft69WK1GeDGwGTiF4nxyKM4jXwSsBV7j2ihiu+c0tKnm7eNVP+9uQr52GvC10Vum2PYj18YR2z2noU01bx+v+nn33qz2wDeJSAf4nKo+CCAijwA6wCbXRhPbPaehTTVvH6/6eXdh1cBfC3wAWCUimyjWVXk8sK7c5to4YrvnNLSp5u3jVT/vLsy/yCMiOwNjqvpz18Yb2z2nobWMnaLWMnZo3mDYwEXk8cBLgT14+JPYr6vqZtfGE9s9p6FNNW8fr/p5T8TkPHAReQOwHjgIeDTwGGA+cF25zbURxHbPaWhTzdvHq37evVgdA38PsH/vO46I7Ah8B/iUa6OI7Z7T0Kaat49X/by7sPom5hiTXPwBeJDJLxTRZq1lbPechtYydopay9iheXdhtQd+GrBeRK4E7igf2wt4EXCqa6OJ7Z7T0Kaat49X/by7sPwQc0fgJRQH8scovpX0dVXd5Np4YrvnNLSp5u3jVT/vicRwGuFOwPgwybdNaxnbPaehtYydotYydmjeYNTARWQv4IPAC4BfUrwLPQ74FvAuVb3dtfax3XMa2lTz9vGqn3cvVh9iXgpcBsxR1aep6lOBOcAXgM+6NprY7jkNbap5+3jVz7ubuounTMet0+ncOsy2NmpTzds9+3jFqk05796b1Vko14nIKmAND38SuyfFkorfd200sd1zGtpU8/bxqp93F1YN/A0Ui5q/j+5PYi8HLnBtNLHdcxraVPP28aqfdxfmZ6E4juM4w2H1IeZWiMi3XDv9ehHZpef3w0XkLBFZIiJjo9KG5DwK/UzXzpR5TqGmYtJanUb4g948gKcDtwCo6rNcO22x16vqvPL+e4HnA5cALwPuVNWlI9Jaem6VttQnN88J15RZ3r1YHQO/Hfgf4P3A/RQm/hV4uWunXT9xD+qVwPNV9TcicgnFqmij0obkHKpvmxbSnGcrbcp5d1P3tJXpunU6nUM7nc7VnU7n4PL3H7t2+vWdTmdDp9N5TqfT2b/T6fxHz7brR6W19NxSbZLznGJNWec98WZ2DFxVLwMWAgeJyOXAI107Ev1G4MPAh4BfiMgcACmuBvLACLWWnlunJdF5TrSmTPPuYtjOP523Tqfz7E6n82bXNqMvn2ObTqfz6Ka0lp7bpk11nlOtKcu8g14cobdOp7PdJI/t4tq4YrvnNLSp5u3jVT/vLTeTDzFFZD5wMfAoEfk+sGTCIi5XAvNcax/bPaehTTVvH6/6efdidQz8g8BLVHU2cC7wDRE5sNxWdd5p27Sp5u2efbxi1aacdxdWDfyRqvpDAFX9v8AhwBoROZTJLzfUZm2qebtnH69YtSnn3YVVA/+jiOy+5ZfS0ALgZOBpro0mtntOQ5tq3j5e9fPuwqqBvwvYbeIDqnoncBBwumujie2e09Baxk5Raxk7NO8ufDErx3GcRLE6C+XxwLspjv/MLh++F/gicLqqbnatfWz3nIY21bx9vOrn3YvVIRQFNgEHqerOqrozML987J9dG01s95yGNtW8fbzq591N3RPHp+PW6XRuGWZbG7Wp5u2efbxi1aacd+/NajXCn4jIO4A1qnoPgIjsBizm4csMudY+tntOQ5tq3j5e9fPuwqqBv4bi09hcRHYtH7uH4rJC4tpoYrvnNLSp5u3jVT/vLvwsFMdxnESJ5pJqWxCRWmsBtFlrGds9p6G1jJ2i1jL2MNroGjhwrGuTiO2e09Baxk5Raxm7ttYPoTiO4ySK1YeYfRGRuaq6wbXTqxeRvwT2pLjCyq11Ylppp0M/yfONfK4sx6vPcyb32vaaGozoGjjFmrh7uXZ69CKSAWcCm4H9gf8H7CgifwSOUNW+py5ZaadDPwUjmyvL8QrJ27X19DHNs9VX6c/qs2kMmOXaadV/BHixqt4nInsDH1bV54nIi4ALgBdHqA3SG87V0DkHapN8bXtNDVUXXVh9iHkkcCNwXc/tWuAPrp1W/Taqel95/6fAkwBU9RvAHpFqQ/VWc2U5Xim+tr2m6s9zF1aHUL4H3Kiq/9a7QUROce206q8VkQuAbwKvAP6l1D0a2CZSbajeaq4sxyvF17bXVP157sKqgb8a+N1kG1R1b9dOq/4Y4E3AXwNXAReWj48DL4lUG6q3mivL8Urxte01VX+eu/DTCB3HcRLF6kPMx1GsiftE4ApVvWTCtlWqepxrp08/xfNeoaoLU9IOorecq2FzDtWm+NoO1O5OcRmyB4H/AxwPvAq4GXirqm7sp50O/RTP22hdWB1CuQi4Ffgc8EYReRVwmKr+HjhwSmX7tEF66f/13DFgvxi106A3mSvL8SLN13aIdjXwFeAxwDrg08DfUxxX/mT5cyR643nuwqqBP0VVX1Xe/4KIvAf4logc7Npp138PyCleIL3MilQbqreaK8vxSvG1HaLdTVU/BiAix6nqB8rHPyYiR41YbznPXVg18EeJyCNU9UEAVT1NRO4ErgZ2cO206m8GjlHVW3s3iEjVlwastKF6q7myHK8UX9sh2omnQH9qim2j0FvOcxdW54F/CXjBxAdUdQ1wItXnYLZNG6o/hf7zfHyk2lC91Vydgt14pfjaDtF+UUR2KDXv3fKgiDwV+FGFNlR/Cnbz3IWfheI4jpMo5svJ9h7Un+Igf+u1lrHdcxpay9gpai1jh+YNETRwtl4Dt86auG3TWsZ2z2loLWOnqLWMHZq3H0JxHMdJFbPlZEVkDDiAYgGXceBu4LuqWvmO0jZtqnm7Zx+vWLUp5z0Rkz1wEXkxsIriJP67yoefCDwVOE5Vr3StfWz3nIY21bx9vOrn3YvVHvhHgReq6u0TH5RifdyvAs90bRSx3XMa2lTz9vGqn3cXVh9ibgvcOcnjdwHbuTaa2O45Da1l7BS1lrFD897qySy4EPieiHwW2PLtoz2B11JclcK1ccR2z2loU83bx6t+3l2YnYUiIvsAB1McyB+jeFe6XFVvmqHaZ1IskFNbG6pvqedWaY3zbuPry8RzL34aoeM4TqKYHAMXkceJyD+JyMUi8rqebatmmrbiea8YVjuIvo2eXdtc7Da+vmLyHMt64K9m+DWEo9dKXGtjz3jPbdMax27d6wu7HrQVvh54/OtEh+rb6LltWsvYbXx9WV4boAtfDzz+daJD9W303DatZew2vr4srw3Qha8HHv860aH6Nnpum9YydhtfX5bXBujCz0JxHMdJFPPlZHs/TJjiw4XWay1ju+c0tJaxU9Raxg7NGyJo4KS5Fq+vXeyeY9Vaxk5Raxnb1wN3HMdpK74eeALaVPN2zz5esWpTznsivh545NpU83bPPl6xalPOuxdfDzx+bap5u+fmtKnm7eNVP+8ufD3w+LWWsd1zGlrL2ClqLWP7euAt06aat3tuTptq3j5e9fPuwnI98OTWXLbSppq3e/bxilWbct4T8dMIHcdxEiWGL/J0IQmuuWyltYztntPQWsZOUWsZexityTFwSXDNZSutZWz3nIbWMnaKWsvYoXn3YvUhZoprLrdxrWf3nIbWMnaKWsvYoXl3YdXAU1xzuY1rPbvnNLSWsVPUWsYOzbsLq2Pgp0wRO9Y1l620lrGttJaxU9Raxk5Raxk7RLsVfhaK4zhOokR3ForjOI4zGN7AHcdxEsUbuOM4TqJYrge+A/BSinUAHqBYXvFKLa/WPJO0IrIXcK+q/k6KtYAXA/OAm4DzVPWBUeoneb4VqrqsjmYYrYg8Dpitqv/V8/izVPUHo9KLyMEU8/K7QfKcLm2pH9pzoNbM8xTP+yJV/caotCLyt8A9qnqLiPwNcCBws6p+ZcAYQ+lD6nG6a9nqizwCnAT8BzAf+DeKwfugiLxeVW+YSVqKZSIPKO+fDjwF+ALF1an/CnjjFNogvYic1fPQGHBE+WaEqp4wIq0AHwHuFZHtgMWq+r1y82qKF21fAvWXAr8pv9n2GeDrqvqnqeJNhzYk59DxCsk7UDsVFwB7jUIrIh+hqIltReTrwALgCmCpiBykqidN9eSB+pB6Du0FXVgdQnkvkKnq0cBzgV1V9fXA4cA5M1D7CFX9bXn/hYCo6lpVfSOwf4U2VP9KYCfgWuC68ucfy/vXjVC7DNhfVfcDjgQuFpFXltsm+xLDdOo3AE8DrgZOBO4WkU+KSDZA3BBtSM6h42XiWUQu73P7ErDzqLTAi4DnAQcBb6FYY/tUYGF5qyJEH1KPob2g+8nqCqaJMeD+8v5vgF0Byn8THzcDtXeIyAvK+7dTHIJBRKpepNOhfybwc4pDP1ep6hrgV6q6prw/Ku02qroRQFW/S/Ffy3tE5ASKy0hVEaIfV9VNqnqeqi4Ank3xL+rpA3xZIkQbknPoeFl5fj7FDsyZk9x+PULtuBaXINty+HLLGD3IYH0tRB9Sj6G9oAurY+BfBb4mIjnFu90/A4jITlTvbaSoPRr4lIicAvwSuF5Evg/sCLy9QhukV9VfAW8Tkf2BtSLyFQZ84w7RAr8SkadsOZ6rqhtF5CCKfxf/14j1XfOhqj8DzgLOEpEnjVAbknPoeFl5vgb4rarmvRtE5JYRar8iIv8K/BlwPqAicg2QUfwnUUWIPqSeQ3tBFyZ74Kr6TopLC/0B+EdVXVFu2kzFsb5EtXeo6nzgGIoXy9HAuyn+Zf7mVNrp0JfPcR3Fcbb7gW8PognUHkvP66t8Q3gpgx3nC9Ev7bdBVX8yQu2xbN0MB805RAtGnlV1oaqu67Ptb0eofSfwDuD48nj124DfU9THcVNpQ/Uh9TgdtTyR1n4Ts9xzHlfVTaloLWOH5u00h89ze4iugYvIDaq67yi0UpzC80GKvclfUuzxPA74FvAu7bnQ6Ai0Cyj22AfWhuoNtXsCZ1BcdeQK4AxV/WO57QuqekiF5yD9FM87steXpdZwnucCKymOHZ8ALAcOAX4ELFLVm0ekNXt9WWknw+o0wlf22TQG7D4qLcXpUh8BXq/laVIisg3QAT5LcVpgbNpU874Q+BzFcc6jgFxEXq6q/w1UHVcN0lu9vgxf12A3z+dSNKQdKBr+OynOonkZ8HGKN4VRaM1eX4barbD6EPNS4NNM/un6n41Qu4uqXjrxgfIF+1kROTVSbap5z1bVT5b3jxeRw4GrpfjSyCD/9oXorV5fVlqwm+fHquqXAETkVFX9bPn4l0TkfSPUWr6+rLRbYdXAfwB8SFVv7N0gIi8cofY6EVkFrKH7itCLgO9Hqk017+1E5M+0/Hafqq4VkZ8BXwceU6EN1Vu9vqy0YDfP20y4/+GebY8codby9WWl3Qqr88DfBvxPn22HjlD7BuAG4H0UA3Zlef9G4IhItanmfT7FF54eQlWvovi3fKsmNc36t2Hz+rLSgt08f0Ie/mbuqi0PishTgatGqLV8fVlptyK6DzEdx3GcwYhmNUIRWe/a+GO75zS0lrFT1FrGDtFG08AZbL0H19rHds9paC1jp6i1jD20NqYGPtASkK41j+2e09Baxk5Raxl7eO34+LjprdPp7NjpdB7r2rhjh+ZdPscuVvq2aa1j+62Zm8mHmCLyBIq1cF9BcRL/XeWmC4HTtnwzybXp5i0iC4FVpeZ4YC3F+cyPoviW3ZTrPoTo26Y1zvsXwOcp1hH/lhYr/A2ElTblvHuxOoSyFrhQVR9PcfrM5yiWLt0W+IRro4kdov0n4O8oLoJxFXCUqj6FYh3mMyq0ofq2aS1j3wdcD/wjcKeIfFREqr5ZbK1NOe9uLHb7O53Of/T8ft2E+xtcG0fsQO36Cffv6Nl2/QCeh9a3TRtR3nt1Op13dDqd9Z1O58edTmdFjNqU8+69WX0T877yK6TfAl5FsbA5Ulwjruq/grZpU817s4gcQ7Eo0iYRWQooxVVIqhbrD9W3TWsZ+6EzKFT1pxSLYn1QRJ4BvDZSbcp5d2HVwN8IfAh4F8W/E28pH9+JYm1c18YRO0S7iOIydOPAi4HXUXzL7yfAmyq0ofq2aS1j91vP+xaKb3PGqLWMHZp3F/5NTMdxnESx2gNHRF5CsfbvHhTv/HcDX1TVr7k2ntgj0H5BVb9epQ3Vt00bYd5RvzZTznsiVqcRfgR4OvAp4M7y4SdSLKpzq6q+1bX2sd1zGtpU8/bxqp93L1Z74H+nqk/vfVBELqW4GsdUJtqmTTVv99ycNtW8fbzq592F1XngvxORAyZ5/K+A37k2mtjuOQ2tZewUtZaxQ/PuwmoPfDFwtog8lof/jdiTYj3kxa6NJraVNtW8rbSp5m2lTTnvLkzPQhGR3SkO5I8Bd6rqz1wbX2z3nIY21bx9vOrn/RB1v/kzqlun0znFtfHHds9paFPN28er3i2m5WQPdm0Ssd1zGlrL2ClqLWMPrY2pgSe3mLqh1jK2e05Daxk7Ra1l7BlxQYf9XZtEbPechtYydopay9hDa62+yLMTxboadwMXAMuA/w3cDKxQ1U2utY/tntPQppq3j1f9vHuxXA/8MRTvPOuA3YEPAPcDq10bTWz3nIY21bx9vOrn3U3Ip7YBn7peX/4c63Q6d022zbX2sd1zGtpU8/bxqp93781qD/wRIrIjxQnsO4jIkwFEZGfgka6NJrZ7TkObat4+XvXz7sLqm5j/BGwo778ROF9ExoF9qF4Tt23aVPN2zz5esWpTzrsLs29iisg2wJiqPiAi2wL7AXep6kbXxhPbPaehTTVvH6/6eXdR95jLqG+dTmeua+OP7Z7T0Kaat4/XYLeYzgPfwpWuTSK2e05Daxk7Ra1l7Npak2PgInJWn01jwCzXxhHbPaehtYydotYydmjevVh9iHkkcCLw+0m2vc610cR2z2loLWOnqLWMHZp3F1YN/HvAjar6b70bROQU10YT2z2nobWMnaLWMnZo3l1YNfBX0+fqE6q6t2ujie2e09Baxk5Raxk7NO8uTC/o4DiO4wyP1YeYjwfeDRwCzC4fvhf4InC6qm52rX1s95yGNtW8fbzq592L1WmECmwCDlLVnVV1Z2B++dg/uzaa2O45DW2qeft41c+7m5AT3gNOWL9lmG1t1Kaat3v28YpVm3LevTerDzF/IiLvANao6j0AIrIbxVWZ73BtNLHdcxraVPP28aqfdxdWDfw1wLuAXER2LR+7B7gcENdGE9s9p6FNNW8fr/p5d+FnoTiO4ySK2VooIjJXRBaIyGN6Hn+pa+OJ7Z7T0Kaat49X/bwnYtLAReQEitNmjgd+KCKvmLB5hWvjiO2e09CmmrePV/28e7HaA38TsL+qHgIcBCwXkbeW28ZcG01s95yGNtW8fbzq592FVQPfRlV/DaCqt1MYWSgiH6baRNu0qebtnn28YtWmnHcXVg38ZyKy35ZfSkMvA3YB9nVtNLHdcxraVPP28aqfdxdWDfwNwM8mPqCqD6jqG4C/dW00sd1zGtpU8/bxqp93F34aoeM4TqLEeEk1x3EcZwC8gTuO4ySKN3CnMUTkdhF54TQ/5/NF5JbpfM4hcvi1iPz5FNt/KCIHNZiS0xL8GLjTGCJyO3C0ql5lncuoEJHVwJ2q+l7rXJyZj++BO47jJIrvgTuNUe6BnwMcAcwBvgAcq6q/E5GXAe8HngzcBLxZVX9Q6uYBFwBPBb4GPAjcqqrvLQ9NrFXVJ06I8XGK07WeVP79ojLGQcBaYBXwduDXwHtU9dOl9vHAx4CFwG+B84AVqvqgiDy1zGE/4I/AN1X1NaVuHHga8ALgE8A48Adgnaq+fOJ/HiLyKOADPLzynALvVNXfT8hvJfBO4E/AMlW9KGDYnRmM74E7TfN64CXAU4CnA+8tG/SFwDHAzhRN/nIReZSIPBK4DFgN7AR8Bji0IoYALwX2Bp5FsdbyFnan+NLEHsAi4FwReUa57WPA44E/BzKKN4Ejy22nAlcCOwJPLP+2C1U9F/g08EFV3UFVXz5Jbu8BDqR4I3g2cAAw8XDL7mUOewBHAZ8QkR0r/DotxWo9cKe9fFxV7wAQkdMoGuHOwDmq+p3yb9aIyDKKRjdO8To9S1XHgc+LyHcrYpylqneXMb5E0SwnslxVf0+xJvNXij+TFRRrNT9HVX8F/EpEzqT4b+ECir3uJwFPUNU7gW8P6f/1wPGqem+Z3/so3rCWl9v/CPyjqj4AfFVEfg08A7hmyHjODMYbuNM0E6868hPgCRSNcZGIHD9h2yPLbePAXWXznuw5JmPiN91+Wz7PFjap6m8myWGXMuZPerbtUd5/B8Ve+HdFZBNwpqpeWJHHZDxhkhgT8/vvsnlPzH+HIeI4LcAbuNM0e064vxdwN0VDPk1VT+v9YxHJgD1EZGxCE98T+K8h4+8oIo+Z0MT3Am4Efs7De9k3Tdh2F4Cq/oxiJTlE5G+Aq0TkalX9z57nr/pQ6e4yxg8nxLh7SC9Oy/EG7jTNP4jIlyn2LJcBlwKfBy4TkauA7wKPplil7Wrg3yk+zHuLiJwN/D3FceN/CcjhfeUhmudSLCR0sqr+SUQUOE1E3kBxvP3twIcARKQD/Ht5+GQTRaP+0yTPfQ/FMfR+fIbiuP/3yuf4PxQfXDpObfxDTKdpLqH4MPDH5e39qnotxd7txyma439SfvCoqn8AXknxgd5m4HDgy8Dvh4z/szLG3RQfOL5ZVTeU244HflPm9e0y1y2HSf4K+E55TPpy4K2qetskz38BsI+IbBaRL0yy/f3AtcAPgBuA9eVjjlMbP43QSQ4R+Q7wybqn1/Wecug4qeOHUJzoKY+D30JxnPr1FKcGfs00KceJAG/gTgo8g+ILLztQfHj5alXdaJuS49jjh1Acx3ESxT/EdBzHSRRv4I7jOIniDdxxHCdRvIE7juMkijdwx3GcRPEG7jiOkyj/H61ELIB6hmu7AAAAAElFTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"from shapely.geometry import box\n",
"\n",
"# construct footprint using bounding box containing all 3 points\n",
"footprint_2018 = box(*evacuated_places.to_crs('epsg:25833').buffer(2500).to_crs('epsg:4326').unary_union.bounds)\n",
"\n",
"# query the api using information from the article and bounding box\n",
"fire_date = date(2018, 8, 23)\n",
"pre_fire_treuenbrietzen = fire_date - timedelta(weeks=4)\n",
"post_fire_treuenbrietzen = fire_date + timedelta(weeks=4)\n",
"results = api.query(footprint_2018.wkt,\n",
" platformname='Sentinel-2',\n",
" processinglevel='Level-2A',\n",
" date=(pre_fire_treuenbrietzen, post_fire_treuenbrietzen),\n",
" cloudcoverpercentage=(0, 30))\n",
"\n",
"# plot cloud coverage to discard unusable products\n",
"gdf = SentinelAPI.to_geodataframe(results)\n",
"gdf.set_index('beginposition')['cloudcoverpercentage'].plot(kind='bar')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"There are a set of products with very low cloud coverage immediately before the outbreak of the fire on August 22, another set of product with low cloud coverage on September 9 and a third set another ten days later with close to no cloud coverage is available.\n",
"\n",
"The results are deduplicated to avoid reprojecting an merging. For this purpose the position of the evacuated places in the downloaded product's tile geometry is plotted:"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAiQAAAIHCAYAAACvy0foAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAACN80lEQVR4nO39e7hsV1nni3/GrFmXdd1r7Us2uWASkCTsgEkwAoKHIqTlBLQjiE6UTgwoEqG17TRN2+KhxUM3P89B3M3vtOCWdJOYiDJbGuTIvQWLpmlu5iJkEyKBRHMx7Ozs29pr1X2eP8acqVq1a61VlzFv9Y7P89RTt1ljvt+at3e+4x3vUEEQYLFYLBaLxZImTtoGWCwWi8VisViHxGKxWCwWS+pYh8RisVgsFkvqWIfEYrFYLBZL6liHxGKxWCwWS+pYh8RisVgsFkvqWIfEYulDKfVipdQ3lVItpdRfp23PrKOUertS6jtp2wGglLpFKfXf+96PbZtS6q+VUjePs548EB4XgVLqvLRtscwu1iGxxIZS6ilKqbpS6h+VUsW07RmR9wF3AE8DfnrYAhNeqK5TSmWq6I9S6jtKqbenbceoKKV+RinVUUr9t5hW8evAz8bUtnFCByF6rCml7lZK/VLado1KFo8JS7pYh8QSJ78IfBw4CvzUTgsrpUqxW7QzzwA+GwTBPwRB8ETaxgwjI/9TGrwB+L+AlymlnmK68SAITgRBcMx0uzHzq8DZwBXAJ4GblVJDnSrB+40lJ1iHxBILSikH+GXg1vDxhiHLPKCU+vdKqfcqpY4C/zP8/IeVUp8J7/qOKKX+m1Lq/L7fXRh+9ohSal0p9Q2l1PUj2HSxUurjYbtrSqn/Vyn1g+F3Lw7v1grAH4d3na8dUevbw2jDTyml7lVKnVZKfV4p9fSobeC28HV0R3tL3+9/LfxdXSn1d0qp31JKuQb+p/OUUh9WSj2ulNpQSn1XKfWW8Lu/Bp4O/HafTRdsoe85SqlPKqW+H67ra0qpawaWeUAp9X8qpd6jlHpCKfWYUur3lFKFvmXKSqn3KaVOKKWOKaXeB5RH/I+fBrwI+H3g82hnN/ruT5RSnxnym08qpf4sfL3jPqN26EoZY79zlFK/G/7vJ5VSNyul5nbQ93NKqbvCfeABpdTvK6UWtv1TNCeCIPjHIAj+LgiCfwt8hzCyp3T30X9WSr1DKfUo8HD4+fOVUl8I94ljSqkPKqXOGrDn15RSD4U6Pw38wMD3r1VKtQc+Oy/cj17c99nTlVL/Ndwn1pVSf6uU+sntjgml1I8ppf6nUupU+LhbKfW/j/BfWHKOdUgscfFSYAH4BPrE8+LwojLIvwC+D/wocINS6gBQA/4XcCXwEqADfFYpVQl/swj8FXAN8Gzgj4APKKWu2sqY8ILwGaACVMPHIvAppe8cv4S+04TeXeeHxtB7NvBG4J8BLwBWgP8SfvelsM1oubPR3QMo3WXyr4HfBJ4Zfn4j8NsD7U/yP70X2AX8k7DtXwIeCr/7aeAB4N19Nv3DFtqWgT8DXgw8B/g08DGl1EUDy/0a8CjwvNDefwn8Qt/3vwu8KvzsR4HTwD/fYp2DvAH4RBAEjwO3AK9XSqnwuz8GrlZKnRstrJTaD/w42hmGCfaZIYzaxs8Ae4D/Db0/XIuO7AxFacf3fehtcQD9//wT4A/HsC1iA+jvHvWAfcDVwEuUjix9Br0fPBf4p8CzgA/32fNTwEG083c54APvGteQcF1fAlbR/8GzgbcBXbY4JkIH9mPAV9D72nOAtwPr467fkkOCILAP+zD+AD4CHOx7/wngnQPLPAD81cBntwB/NvBZGX1CesU26/sL4P3bfP9LYRt7+z7bjz6B/0LfZwFw3Q7a3g58Z+B9G9jX99nPoU+8lfD9dfpw29TOfGjTNQOf/wJwfNr/CbgbePs2Or6z3fc7/Ad3A781YOPHBpb5FPCn4esFoA788sAyX+//L7dYVxF4DPipPp1PAC8N3zvou//f6PvNv0I7R4VR95nwP/3vW23nEdv46/C/KPR99gagASxssZ4HgF8ZaPdF4b64us26n9xXARd4ffjZr/TZch/g9P3mHWhnpNT32WXh714Uvv8i8CcD6/q9cJnzwvevBdoDy5wXLvPivnX9Y6R7iP3DjonV/jbsQ9bDRkgsxlFKnQ38JL27U9An4depvq6IkK8OvP8R4JWq162yhs5BqaDzO1BKzYch8XvCUPAa8HLgfLbmUuBwoO+wAQiC4DHg2+F30/JIEARH+t4/DCjgrC2Wj2yaAz48oPcQsEspta9v2bH/J+A/Am9VSn1FKfV/KaVeNIkwpdQ+pbuL7lVKHQ/XdSln/t93Dbx/GO30ge4eKqPvjPv54ggmvBLtdHwCIAiCBjpi84bwfRf4E6C/++R69EW1E2qYZJ/ZxBhtfDVab8j/BEro/2CwzX3h739/YFt+MlzkB3cw6+Zw+To6qvG76P0n4m/C/yfiUuDLQRA0ow+CILgbOEHvODjAZNtpkB8GvhQEwelRfxDoHJ6bgU8r3eX2b5VSF0+wbksOGbw4WCwm+CX0vvX1XlQd0PkZ1wL9oyQGT1YOuovnd4e0ezR8fhc6SfbNwL1hG+9Gd09sx7CMfrXF5+PSHHgftbmd0x9997PoO9lB+pNqx/6fgiD4gFLqU+guhquATyqlPhIEwXXb2DSMW9A5BP8G+B46qvRn6ItsP8P+g0ij6vtsXN4A7AU2+vYnBXSUUvtDx/JW4C1KqR9GRyMuB27oa2PSfaafSdtQ23wX/T+/js6NGeShIZ/181voKM1p4B+DIBj8f4c5A1ttg2CEZSK6Qz4bNpJu7O0dBMEvK6Xeg+72/XHgHUqpXw2C4NAOP7XkHOuQWIyidDLr64F3An868PVvoC8u2w3b/DrwQ8D9Q06uES9C3/1+qG+dF6HD+ltxD/ArSqm9UZQkzDO4CB2OjptmuM5C393zPeg726cFQfCJMdsb5X8iCIJHgQ+gcx0+AfypUupNQRCcDG0qbPXbPl4E/JsgCD4WalhAD4v+5hj2fidc3wuBw32fv2C7HymddPwSdJRkcKj1nwOvA343CIJ7lFJ3oLu7GsBdQRD87YCGcfeZQUZt40cGtvOPorXfP9hgEASPKaX+Abg4CIL3j2FLxGNBEIwzBP0edKSyFEVJlFKXoZ2qe8JlDqO303v7fvfCgXa+DxT6HELQ+R79/A3wy0qphS2iJMOOCQCCIPgmev/6faXUH6LPG9YhmXFsl43FNNeg76YPBUHwzf4H+sL442qL0Rwh70QnYN6ulHqu0iMbrlJ69EaUFPtt4KfC7w+gkwvP2cGuDwJHgA8pPWrkh9F3+Q8zXvLqpHwvfL427AJZDIJgDa33nUqpX1V6FNClSo+42DIJMmTH/0kp9Z+UUi9XeqTDpehE1n8ATvXZ9EKl1A8opfaGF9hhfBv4Z0qpZyulLkc7mqM4Mk8SXpD+EPj3SqlrQ63/N3DJDj99A/DdIAg+OmR/8tmc3Hor8PPoRNI/HqJh3H1mkFHb2AP8gVLqmUqpn0DnUrx/m66L3wL+hVLq/1BKPSv8b16hlIrjAvyf0EnKt4Tr+jF0pO2LQRD8j3CZdwOvVkr9ulLqGUqp17G5Owx0F+Ip4HfDZa4B/t3AMu9FX2P+Qin1wnAf/Uml1MvC7884JpRSPxh2L/6YUup8pdSPopODD2OZeaxDYjHNjcBXgiD4+yHf1dBOweu3+nEQBN9C3zUvokdzHAbej861OB4udhPwIDrE/Vdop+LPtzMqCIINdAi4AXwhtOU0OqF0sKvBOEEQfA14D/qi/Bj6wkAQBO9A63k9OlH0i+H7B3Zob5T/SaHzSL6J1rwAvKwvovLb6Dvjb6O3y6ahnX28Dn2u+CrwUXSy6tdG0T3Avw1/f1vY1grwB1strPTop9cC/3WLRT6Ezsu4Onz/wbDNs8LX/Yy9zwxh1Db+HH2x/iLa6f0EurtrKEEQ3IYeDfMT6P/la+iE2ofHtG9HwmjGS9EJqF8D/hK9f7yqb5mPoLul/g3wt2gH7zcG2nkC7fw9P1zmbQxoDKNzP4b+Lz6BjsD8B8IurC2OidPoHKg/Q3djfpjNI3IsM4zaJtprsVgsFovFkgg2QmKxWCwWiyV1rENisVgsFosldaxDYrFYLBaLJXWsQ2KxWCwWiyV1rENisVgsFosldaxDkjC1Wu3laduQBla3LKTqBrnarW5ZxKHbDvu1WCwWi8WSOjZCkjC1Wu32tG1IA6tbFlJ1g1ztVrcs4tBtIyQWi8VisVhSx0ZIEqZWqx1M24Y0sLplIVU3yNVudcsiDt3WIUmed6ZtQEpY3bKQqhvkare6ZWFct3VIkud1aRuQEla3LKTqBrnarW5ZGNdtHZLk+WTaBqSE1S0LqbpBrnarWxbGdVuHJHmek7YBKWF1y0KqbpCr3eqWhXHd1iFJnkfTNiAlrG5ZSNUNcrVb3bIwrts6JBaLxWKxWFLHOiTJc3baBqSE1S0LqbpBrnarWxbGdVuHJHnuSNuAlLC6ZSFVN8jVbnXLwrhu65Akz8vSNiAlrG5ZSNUNcrVb3bIwrts6JMnzgbQNSAmrWxZSdYNc7Va3LIzrtg5J8rw1bQNSwuqWhVTdIFe71S0L47rt5HoWi8VisVhSx0ZIEsZOVS0Lq1seUrVb3bKIQ7eNkFgsFovFYkkdGyFJGOtNy8LqlodU7Va3LGyExGKxWCwWy5YcOnh4P3A+unDZPmAvsAqsAMUbbzrwi+lZtz3WIUmYWq32vmq1+sa07Ugaq1sWUnWDXO2Sdd97x743ARcC5wFPAc4C9qCdgF3AEjAfPhaBClAGSkAxfEQ9Fg6gwtf9z9HrYfR/HwDd8DkAOuH7DrAbuPbGmw58elK9EXFsb9dkY5aReEvaBqSE1S0LqbpBrvbEdR86eHgR7QicA+xHRwSiaMAu9MW/3wmIHIEy+vrnAoW+Jodd9Ic5CH3sVsD19JyAbt/rDpsdghbQBurAifB5Azg9pvRJuQr9v5jA+Pa2DknyvBn4nbSNSAGrWxZSdYMA7cO6BSrzCy+/947D3wKW0Q7AXN9zCe0M9EcDFJsv9oMM+26YY9B/8e+PDrTC52b4WSP8bC38zIgj4Jbaz2g3C383bTsJ0UVvBxMY38+tQ5I8H0zbgJSwumUhVTfEpP3QwcOK8boFSmhnoBi+7o8GbOUIRN0G/Rf+oO91f7dA/8W/U18vKuCi8LN236MFrAPH0U5Am+QiArHTbjoPp23DGARoh9EExvdz65Akz1VAXrxpk1jdspgZ3WG3wDloZ2AVOBd98d+LvvivorsAdgElWDn33jsOr6GdAUd/tlPYf9toQP/rwfyAqCsg6hbodwKa6G6BJtoBiNUJKLjdH+i0+fs415FFCm53b450B2hn1QTGj3HrkCTPvWkbkBJWtyym1n3o4OEl4JnoRLzz0JGAVfSFP8oPWAgfg/kB/U5A9Bxw5kV/VEdgWH5A5AD0OwJNCNaAY8BjJJsfkCrdrlpL24Y0yKFuUw6J8XObdUiSZyVtA1JiJW0DUmIlbQPGIewWOBvtAJyL7hLYh9axwvD8gAqbRwsUYLV07x2H24zWLbDV64AzIwH9iYIdNkcE6sBJUu4WcFz2d9s8lsa600Q5FINu2lYkT850d9EOvAlWDLXzJNYhSZ65tA1ICat7Gw4dPOwAFwBPRY8W2I+ODKww2bBB2DkasFVkYLthg8O6BU4y0C1QcDvndNruI6NonzUUQWHnpWYPqzsXdNHnChMYP6dbhyR5vpe2ASmRSd2Ghg1G3QJwxkV+1b33jsN/dObnI3cLDMsPaIaPYxgcLWCSbletp21DWkjVbnXngi76fGYC4+d065Akz/OBr6ZtRAps0n3o4OEF9AU/6gJYQocSIwdgnl63wDzaEZgLX5fpRQQGkwUHawpE78thWyvhegJ6eQb9F/9BR2CwW6CBjgjUw/fbOgJusXNBu+U+sOO/M2MUCsFKu8vxtO1IA6nare5cEGAusmH8WmYdkuT5SForDvMDooTApb7HApudgHl6zkFUO2ABfQEvh81FjkDkDEQX/sH34XK7S/fecfhn+34XXeSjC35UMKhLrysgejTQtQMep1dLYNy7kgvDtp+JdiIK6EhIPbQhFtqtgrhcApCrG+Rqt7pzQYC5OiTGr2XWIRmR8I7+L9AXsFPh4yjwBHp8fZRVfwT4xxtvOnB8i6beBPzmhDbsAz6MvqBuFw2AMx0GB70zDjoCXfRFP/osSg5so7sCjqMv/q2+57FZ3lP/qZNHF/5ikt8awqWXc/E4uoYDxOiMALilzvntpitupI1U3SBXu9WdCzroc6AJJr6WbYV1SEbnQvTF+OP0cgsW0Xfc0UiDImF+waGDh102h/7Rz3uCe+84/Bl6IwY20F0AJ9EX/6N9j8eBR4DHb7zpQBO4JFzu4zFrNc7Jo3MfS3H1RbTjFYUqo+Gfsc923W4Wvx33OrKIVN0gV7vVnQs69KLc0/I2Q+08iXVIRmeVXnTgRPgYv5F969cfO7J4W99HRXrdKGehR1pU0N0m0QiKUtjdUkY7Rs+gNwICdDQjoh4+1tHdHOvhI4rqRJ8lyuq+9esGdCfJCrqrZzF8jpJR29v8xgilcvOyZqNyV9zryRpSdYNc7VZ3Lui/MZuWm4HXGmoLsLP9jsyhg4d/Engd8N9TNGMXege4Y5tlolEg0YiQ/iGhJXS3Tn9XThSpiUZ3ED5HeR4baAcmGskRPdboOTgTdeMkyNPRmi5G21xB/5fr9DRbLBbLrPNM4PiNNx14ftqGDMNGSEZnNzrXYipW963dcOzI4q0T/nydzXkiw4giJBNFcAaInJl5dPLrnvC9y+aZMqMqmP21K6JukQ6AU2id0+0U7++z7yTJRW8ceqNzjqCdkX4HLDZK5foVzUblzrjXkzWk6ga52q3uXNDBUB2SWq12W7Vavd5EWxHWIRmdZQw4JFM4I6AjEVvNjBkHLSYb0XIG3U7xb9kcvTmLzcN3+6M3ER16zk6UhzNu9KYQLhsN9WXIemIjRycqo0jVDXK1W925oIWhpFbTzghYh2QcFtF39lOxum9tMIdEBOW5+nMbG5WvYj56U2H76M18+HohfN9l834fOSlRJdI2Os8kGnJ8ml4UZ+ziY6Vy/fIc9S8bQ6pukKvd6s4FxpJaa7XaLdVq9bUm2oqwDsno7ELfiU/FsSPzt0/ZRC5zHhobpa8ZbnLU6M0PoiMw56AdiihKMyzaFU3THhVk6x86PTgiZ6vkq/58nFazUWoBz6NXUTWK3phwyjJLs1G6O20b0kKqdqs7F7QxVzr+9YbaeRLrkIzOIvpiMhXLezaunbIeRy6zkIvl1rNajfI3Ulh1gV6Z93b4vNVw38jJMcUctCpQrnDmzLPRY3B7BgPPsDl6ExWKa6KdsdNMGL2JE7fUurjdLOelNoNRpGq3unNBm53zEEflHdg6JKmxCHx/2kbWT5VrUzaRS4ek1Sjcn8Jql9AX7xV6kaVhTkBcbEChgZmoVhS9qaC7nwajN4Mz5NL3Hnq6n4ze9D2abC74N3X0pt0sPDhtG3lFqnarOxdEZQ9M8F5D7TyJdUhGZxEDyZ1z883LTjXdaZySJJNajVEsdc5pNd3vJLzaqHZMGX0gRnPXxFqddTOdCrgmRg7FEL3BRTs4w6I3Edt1TcGZ0ZswgtMpgfswGYzexI1b7OyXOH+R1Z0L2pgrCPlK4P9vqC3AOiTjMIcBh6TVLDwwZRO5jJB0OuqJFFZbRocnC+httxB+HntBtB5q6pFZMTF1PlQfBXp1bxYABaqAHkm1Xe7NsO6q2KM3cdPpqONp25AGVnduMHXd/7Khdp7EOiSjY8SrLLjdPcA0Ib5cJrU6TrDQ7ZCkUxJtrwr6gtkJn2MvF7+ZIJGKsCnT4YwoSFBh8lFpO0Vv+rvdBp/7vxuM3kTzNMWae+M4wXwnP7O/GsPqzgWnMXcOvBA7229qGOkqCVDTht1zGSEJUAl2kwB6VFQTPRw4cggUiTt0Kpfba3qm0h1z9GaTk9PPMJv7HZ4uvZmpo+6pM6I3KezrmcDqzg2muv1NHqeAdUjGwchG7HaV8Y2YB4Lu9EXlxmQX+gJSQV84KvQKqiVJLiNaBsiK7iHRm6mYQzs50bDw6PFk7k237fQnEMNwR6dLr3sqsehNnATdzE8hEQtSdYP5qJB1SEbHyFApt9DdD0yT3JnLO27HCZa7HR5PcJX9c/hEE+ulECEJXLI/108MzKzuEW4oOnPgjnLjMWn0ZquRU202F/ZLNPfGcYLFTjfRbtlMkEPdpob9XgJ8wVBbgHVIEqdRd+9L24Y0aLedxxJcnYs+WZfRJ/YuCZaL34yT1aTWmJGqG8bQHkf0xmGb6E3IdonFMGH0ptN2krzhyAxSdQOfN92gdUhGx8jFbG6heWWrUfzkFE1kJRQ+FsVS9/xWg8MJrW4FffLcR2+YWwrREYBuhZyF3s0gVTekqD3l3Bs1bEj9dtGbyLk5ZtDuxHFL3XPbTf4ubTvGwFSU/TXA7xhqC7AOyTgYCXOdOlb53JRNRKV/cxUObzXcbye4ul3o/8dFnwCjuRtSGO1ipAZJDpGqG2ZE+wTRG7Vd0cFotFsq0Zs4aTfdNIo+ZoF3m27QOiQJs2vPxiuOP774oSmaaJBDh6RUaV3WrBfuSGh10ZBRFziJjpg4pPKftRahcCr59aaNVN0gV/u2uqeemLSPaXNv+t9vF71Zh53z3orl1qWtRiGNaTHS5l3AG002aB2S0TEydntKZwR6lUdzdRfWrFeSckZK9EbXRAzehSVIReCFCeTqBrnaE9NtOvcmmutqouhNq1EBOI/h0ZsG2sE5CfyjQZunwci1rFqtGnVGIPEiUbnGyAVt9761G6ZsIiqBnivKc/XnJbSqVbSzNoc+IUQnlpRqBdSX01lv2kjVDXK151Z3A+3gnEQPZX0CHRk5gp6/LHoc6Xs83nvUm+HrJ9BdRqBviFaApwBPA34EeHYSYpKiVqtNO3P9GdgIyegYcUieOLJ465RN1NE7eq5obFS+ktCqltB3KUV6jkk09DcFKifTWW/aSNUNcrVb3WwdvdlDzqLaO1GtVq8z3aaNkIyOkf/KQIRkjRw6kglGSArAPPo/aqIdkxTJ7V3jlEjVDXK1W93boOhFT2aCOCIk1iEZgUMHD5cwFPI3ECExOX10YiQUIYm6aebpDfFNeXZke9coD6nare5tcNDdPlnAyLBfGyFJj70YGjK6a8/pV03ZxAbmKu0lRqlcvzyB1eymN8y3iXbchtVGSJD6YnrrThOpukGudqt7G1K+MdqEket+rVY7aKKdfnJ3p50SuzE0ZHTtePlTUzaxQQ4dyWbDTaIo2lz4XER3bUUjbVKcbdcVWhxMqm6Qq93q3oYsTflhyjl6p6F2niR3F7aUWMWQQ7Kw3HzBlE2sk0NHsljqPC2B1UT1CaKaI1GEJEU68+muPy2k6ga52q3ubciSQ2KK15lu0DokoxFNZT819Q33nimbWCeHXTatlvNozKtYojeJXn/+SMonAsdkQagcIVU3yNVudW9DlrpsTDHNFChDsQ7JaESlyKemWO6cN2UTa+Rwu7ludyXmVayincYyvVotKeePAHRTHuWTFlJ1g1ztVvc2zGKE5DmmG8zdhS0lljEUIel21bSZ6FGxr1zRDVTcdUDK6ChJAb2tMpA/AqByORni9EjVDXK1W905wVS0xnjUO3cXtpRYxOxcDBazRPtxVAStjXZM7P5tsViywCxGSIxjT9ijsYwhh8RxAhPFg3K3czsqKO+81MREOT5z9LpoMpA/AhAIPcak6ga52q3ubZjFHJKzTTcodAcam2UMlf1tNQoPGWgmAxfa8Wi3neMxNr8L7YgU6XXXKDIxI7KTARvSQKpukKvd6t6GvHXrjILxCVOtQzIaCxhySCpz7UsNNJM7h6RY7Br3pvubR3eruehIVjT5YAZOAt3KzsvMIlJ1g1ztVvcWFEg9uX4Tpq4fLzPUzpNYh2Q0FtCjW6bm9MnSl0y0kzdazcJ3Y2raRUdD5sLnLpnKHynM1IRaoyNVN8jVbnVvQYVMRGuN8wHTDWbkpJ155jEUIVlcaVxjoJncRUhK5faBmJpeQW+bBfRdiEPPMckA7YW0LUgHqbpBrnarewsKpDbbeKy81XSD1iEZDWOFyE4cXfiwgWZy55A0G5W7Ymp6F9oBiQ76cvg+5eG+ERUjkbX8IVU3yNVudW+BQ7ZGaRq57ler1ZtMtNOPdUhGw9j/tHvf2g0GmsldxnZ5rv68mJp20fVH+vNHMhQhsVOyy0Oqdqt7CwoYirBniVqtdrvpNq1DMhrG/qcnjizeaqCZjFxsR6exUflKDM2W0P/FXN9nUZdNRrBTsstDqnarewsUcCIJS5KkWq1eZ7pN65CMhrEuG0MRktx12cQUIVlF33nMobtoHDJRLr4fe9coD6nare4tcIBTSViSJDZCMgMYipDkziGJKUKyhB7y66LzR6KRNhnKaLd3jfKQqt3q3gIFnE7CklE5dPDw1FFkGyFJD2NdACt7115toJncOSSlSt34REzoyFVUf6SBdk4yRn0pbQvSQapukKvd6t6CrJ2vA+CsaRup1WrvM2DLJtxRFvI87wF0yKkDtH3fv9LzvHcB/xRdGfN+4HW+7x8f+N1TgT8GnoLu6/8j3/ffY8z65DA4ymbuowaaaaMvvhmKBGxPs16823CT0dj+ZXoHfIZyRyKKQkceSNUNcrVb3VuQRYfExBDttxhoYxPjREiu8n3/ct/3rwzffxZ4lu/7PwTcB/zmkN+0gTf7vv9M4PnAP/c8L656FLlgabX+EgPNZDQasDXFcvtiw03uQf8PkWPikrn8EYD2fNoWpINU3SBXu9W9BVm7UQqAvQbaebOBNjYxUoRkGL7vf6bv7ZeBnxmyzKOEUxT7vn/K87xvAecChyddb0oY69raOF36uoFmWuTMIWk1nQcNNzmHjrq56MTWaPK+jNQfiXCyVH8gQaTqBrnare4csWqgjQ8aaGMTo15oA+Aznuf9jed5bxjy/S8Cn9yuAc/zLgCuAOJIbowbYx5uudK+yEAzDXrzteQC1+3uN9xkAV1Bt4DuNiySyZyobq62kzmk6ga52q3urRZIxIzRCdAFJaflKgNtbGLUCMkLfd9/xPO8s4DPep53r+/7XwDwPO+30Helf7LVjz3PWwQ+DPxL3/eHZiSHjs4bAC655JJvXX311e9HlwWfA76H7vL5CPAm4G3AzdVq9bW1Wu22arV6fa1WuwV4PfAO4L3AK9GRmwuBDeA4cAnweeA1wLuBd1Wr1TfWarXbq9XqdX3PB4F3Aq8DPukUVs8qz7cvAXCcYLnVKDxUmWtfevpk6UuLK41rThxd+PDufWs3PHFk8dboeWXv2qtPHJ376NJq/SUbp0tfL1faF7U7zmOOEywuLG88p9N2jhZLnQs21kt3zy81qiePzn1sdd/6dceOLN62um/thmNHFm9d3bd2/bEj87cv79m4dv1UuTY337ys1Sw80KirlaBL2SkQOE6w3G47jxVL3fNbDffbpUrrsma9ckd5rv68xkblK9FzqVy/vNlwDxdLnae1Ws6jrttd6QaqAeCooNxuO8eLxe7ZrWbhu6Vy+0CzUbnrjDYq9ec068W7i+X2xa2m86Drdvd3u+qkcigpgkK3q04XCsHuVrPwSLHceXqrUfxmea75I42NylcLpeCp7RbfLc/Vn9vYKH2tWG49q9Uo3F8sdc7pdNQTjhMsBKhO0KW5k6Ziuf7CVqPyMDQvBKcNnVUouBAUIQgP/sDRs3B2K3quifaCrqhYX9ZZ8U8+L+k+4Pa8vtPplkCFQ4gDpV8HJSjUoTOvl20OtrEMpZPQWtTr6lRANSEIu5GaRf3aaWp73HW9bOXUEHsWwT2t1+XUoVsElTFNKgC622tyNmZP06jbiQK0CrOlaZTt5GzMnqZRthNAY34bTejzb+lut9S6uN0sPOgWO/s7HXXccYL58LzXcpxgsdN2HndL3XPbTff+Yrl1aatR+UapXL+i2ajcGT0XS/VLW033PrfYOb/dch4ruN3lbqCa4bm81Gk7J91id3+7VXiwWGpf1GpW7hloI1jec/LGWq32CXS3ywfRzsW941xzgVfWarX/Muo1t1qt/s6wa38/KgjGy7fxPO/twJrv+7/ned4NwK8AV/u+P7QSned5ReAvgU/7vv/7Y60sIxw6ePh/AkbGXM8v15+9frLyjSmbeTHwVOCBqQ1KiEKxeU6nVXrEUHNPRTvBF6FzRpro5NYumeuyaZahNIvzWOyAVN0gV7vVvQX7gE8kZc0IvAh45403HXj3NI3UarVrq9XqxwzZBIwQ4vY8b8HzvKXoNfBS4Jue510D/AZw7TbOiAL+M/CtHDsjCoNZ0orARO5HnZzlkCgCYyOV0Imsi/Tmr4mm/86YMwK9OyhpSNUNcrVb3VuQtS4b0OfPaZnbeZHxGKXPfT/wRc/z7ga+Cnzc9/1PAf8JXZjqs57n3eV53h8CeJ53jud5kTf4QuB64CXhMnd5nvdy0yJiZi8GL3SdtnPUQDN1pkhIToNuV5kuDLSA3n/b9EbYZBCVQScpCaTqBrnare4tyJpD0sXMDe33DLSxiR0var7vfxe4bMjnP7jF8o8ALw9ff5HsDXkalz0YrPdRLHUuqK8z7YiTDTJ7AR5OoRDs7nZ4wkBTu9BdNHvpDfE1GsUyS1Aik5GbuJGqG+Rqt7qHME/26kV1MVOH5PnoIIUxcnVRS4lV9AXQCBvrJRMFwjYwWKwtCVrNgqn8kRX0AV5Cb5cymSsX308hj0MCDSBVN8jVbnUP+xLdrZwlumhHaVo+YqCNTViHZGd2YfBiN7/UqBpoZp2cddkUy52nG2qqRG8Omzq9+iNZC4uGdIQWi5KqG+Rqt7qHkEWHJMCMQ/ImA21swjokO7OCwTDkyaNzJrKS18lZV1irUfymgWYctO4omaqLPuAzvB/bctrykKrd6h6Cgz5fZ4mA3kCAaXibgTY2keETeWYw2mWzum/dxAyJ6+Rs25Xnmj9ioJlVdHfVPJvzRzIaHQFdj0AiUnWDXO1W9xAc9DxwWaKDmQj7zQba2ESuLmopsYjuGjDCsSOLtxlqKlcRksZGxUTyU1RrpExvuG800iaj2CnZ5SFVu9U9BAUcS8qSEelioNJ3tVp97fSmbMY6JDuzjL4rN8LqvrUbDDWV4ajAmZTn6s8z0EwRnT/ioB2SjOePgK7UKBGpukGudqt7CAowXfJgWjr0zp0TU6vVTN1cP4l1SHZmGTDWN3rsyOKthprK6DDX4TQ2KtPOYeSyedrsLr2ckgxj7xrlIVW71Z0TAgwUNQvLxxvFOiQ7s4DBpKTVfWumNmKuHJLyXP25UzaxG70d5tAevhM+Otv9KH3sXaM8pGq3unNCGwNdNuH8cUaxDsnOLGA0QjJvZE4cMh8Z2Exjo/S1KZtYRu+vLr3umgzXH4ko5e3uyRBSdYNc7VZ3TuhgZrb41xtoYxPWIdkZo5X2lvdsXGuoqVxFSIrl1rOmbKKAzh9x2Vx/JOO0TMwZkUOk6ga52q3unNDCzCibdxhoYxPWIdmZAgYdkvVT5ZqhpjKcyHkmrUbh/il+XkGHGefpOWI5iRAVslaDICGk6ga52q3uIWTxPN3CzFw27zXQxiasQ7IzRv+jufnmGfMCTUi+IiSlzjlT/HyVXv5Ii14xtIznjwB0TBQgyiFSdYNc7Vb3ELJ449TGzHXtlQba2IR1SHbG6A7VahYeMNleXuh01DQT6y2g+zyj/JEKertkuP5IhDJWVC9fSNUNcrVb3TkhmiF9Wr5soI1NWIdkZ4w6JAW3u8dQU7mKkDhOMM3skv35I010uDGLdx5DCHI155A5pOoGudqt7mFfJmbG6JjKIbnQQBubsA7JzhidVTdAZXxUSDwEqEm7V6Kk4nl6/bGKbB7oQ1A5sdM0UnWDXO1W9xCy2K3cxsx1zVjB0AjrkCRMt6uMb8Q8EHQnng9oFV3psILurnHJTf4IkM2ktiSQqhvkare6N1Mmm93KpzHjkBw30MYmrEOyM0bDkG6hu99QU7mKtDhOMGnxoDl0Dkl//ghk80Afgg1jy0Oqdqt7gCIGJ2Y1jIku70sMtLEJ65AkTKPu3meoqSZmppBOhHbbeWzCnxbQ3TXRJHpFcrXfOlk9IcWMVN0gV7vVPUABfROVRUxESD5voI1N5OjEnhpmh/0uNK801FSL3BQHg2Kpe/4EP1tGO16L5DJ/BKCbG6fRLFJ1g1ztVvcADgZnis8grzHdoHVIdsboaI5TxyqfM9RUHTPFbRKh1XC/PcHPovojUf5ICb3P5qi7yhVaLEqqbpCr3eoewMHgtCMZ5N2mG7QOyc4YdUh27dl4haGmclQ+HUqV1iQF4Uro4b4Fel1UAblKnstdWWlDSNUNcrVb3QM4wKkkLRkDE9e1dxloYxPWIdkZo//R8ccXP2SoqdPkKELSrFfuGPMn0f8+T68IWoHc7bOVrJ6QYkaqbpCr3eoeQAGPJ2lJklSr1TeabjNnJ/dkOXTwsPF8hd371m4w1FQdwyOA4qQ8V3/emD/Zhe6mWaA3xFeRq+gI5HBqckNI1Q1ytVvdkqjVaqZmrn8S65Bszy4MDy994sjirYaaqmO4aFucNDYqXxnzJ5FDUqbXXRONtMkRlbxNTW4IqbpBrnarO0dM3WVTrVavM2FIP9Yh2Z69GE6gNBgh2SBHDskEEZIS2ilx0c5XKfzcRkhygVTdIFe71Z0jpo782whJ8uzDcGEbgxGSXDkkY0ZIXPQBs0AvibVAbuav6SeXd08GkKob5Gq3ugfIcnkCGyHJIbswHCHZtef0qww1tU6OckhK5frlYyy+ina4ovwRB30A5aVcfB91oSMPpOoGudqt7gGy7JBMTa1WO2i6TeuQbM8KhnMW1o6XP2WoqXVyFCFpNtzDYyy+jI6KFOnlkShyVX8kwj2dtgXpIFU3yNVudQ+Q5YiuCdveaaCNTViHZHuWMNxls7DcfIGhptbJ0fYrljpPG2NxF+2UFOjlj2T54N6GznzaFqSDVN0gV7vVLYzXmW4wNxe0lFjCcOnf+oZ7j8HmchMSbLWcR0dctISOSkX5I5Dr/dSZ5dLR2yBVN8jVbnUPkOXzs4nBAZ800MYmcnyiT4Rd6FwGYxTLnfNMtpcXXLe7MuKiUbn4ObRj4qAjJTnMHwHo5qZ4nVmk6ga52q3uwS8SNWM8TFz7n2OgjU1Yh2R7ljHskHS7ymQmepY98E10AzXqrJdRqfgSOn9kLvw8Z/VHIlSWT0oxIlU3yNVudfdRJLc3USMzatR7ZKxDsj0L6Lv1rDKLJ4ACvfyRBvrAzmn+iMViEUoJw/mHhsnkOdU6JNuziGGHxHGCPBbRmRpHBaNMBDiHHkmzQM/ZMl6+P1kCoceYVN0gV7vV3Ud0QzXLnG26QaE70MjMYdghaTUKDxlsLjcX6nbbOT7CYrvRXWRz6LsLF72P5jj06eRwqLIJpOoGudqt7j6iEYKzzLgTpu6IdUi2p4jh2heVufalBpvLjUNSLHZH8abn0P+5i3ZIoqhKTvNHALqVtC1IB6m6Qa52q7sPB8P5hxnkZaYbtA7J9hj/f06fLH3JYHOZ7AccRqtZ+O4IixXQSa0O2iEpkvt9tJDlHKQYkaob5Gq3uvtQwPGEDUmaD5huMOcn+9gxfsFfXGlcY7C53CS1lsrtAzsssoTuc+3PH3HIURRoOO2FtC1IB6m6Qa52q7sPB/h+0paMgYlr21sNtLEJ65Bsj3GH5MTRhQ8bbC43F+tmo3LXDov0z1/THx3Jeb90ZS1tC9JBqm6Qq93q7iPr0WsTk+vdZMKQfqxDsj3G54rZvW/tBoPN5cYhKc/Vn7fTIugcEoWOlFTozfSbY3I5NbkBpOoGudqt7j5yc26elFqtdrvpNq1DkjBPHFm81WBzudnpGxuVr4yw2DJ6n2zTG2GTc+yU7PKQqt3q7iM35+ZJqVar15lucwZO+LFiPOxmOEKSm+jBDhGSXehumnl6Q3wVOdK3NfauUR5StVvdOWJqh8lGSJLHeJeN4QhJC921kXl2iJCsoMfsV+gN940iJTnH3jXKQ6p2qztHmMghsRGSvLOyd+3VBptr0KvVkWlKlfp2EzGV0MmsLtoxiTTNQoRkKW0L0kGqbpCr3eqWRK1We5/pNl3TDc4YxiMkJ47OfdRgc3X0aJTM06wX797iKwftrS+G77vo/z3rWeojUhQ68iB53RsbJ89Nep1zc8sPn/mp3eayGKp75nNIgLeYbtBGSLbH+EVxabX+EoPN5cYhKZbbF2/xVf9w3yh/JOfl4vtpz6dtQTpI1Q1ytVvdeeLQwcPTXt/ebMSQPqxDsj3G/5+N06WvG2yuju7uyDytpvPgFl9FCa1z9Ib7KnJffyTCmfX5LLZAqm6Qq93q7iPrEZIAOGvKNj5owpB+rEOyPcYjJOVK+yKDzW2QkwiJ63b3b/UV2ilx0A5JLhys0enOmJ5Rkaob5Gq3uvvIeoTXhENylQlD+rEOyRYcOnh4hRju0tsd5zGDzdXJyTbsdtWwTHQXfWAs0CuCZjxvJ13UDIwUmgSpukGudqs7ZIHsR3gDdHf5NNxrwpB+cnExS4m9xDDs1HGCOYPN1clJYrJyhkY+dgPr9PJHnPAxA6NrnkTqMSZVN8jVbnVrCuhob9aZ1iFZMWFEP1J3oFFYQec2GEURmOxiWScnEQVFMMzOJbTzUaI3hHmG8kcAghkZLTQuUnWDXO1Wd0geHJIu05eMMHlzDViHZDt2E0OEpNN2jhpsbp2cREi6XXV6yMf9+SNR/ZEZO6nZMLY8pGq3uqMPgGHnu6wxbYXZ7xmxog/rkGzNMjF4ucVS5wKDza2Rk21YKAS7Bz6qoB2+eXoZ6YrsZ6ePSSA00U+qbpCr3eoOcYBTaVgyBh30zeA0PN+EIf3k4mKWEsvE0HWwsV7aqkDYJOSmy6bVLDwy8NEq2v559P/soLXM2F1WQehQSKm6Qa52qzvEAU6kYckY9BejnJSPmDCkH+uQbM0iMURI5pcaVcNN5iKiUCx3nj7w0QJ6/4vyR+bQB0nWh8uNSSeXRZOmR6pukKvd6g7JQ5dNF30zOA1vMmFIP9Yh2ZqogqhRTh6d+5jpNvNAq1H85sBHBXRSawGdPJyLeirjY8tpy0Oqdqs7R0TlFqbhbSYM6cc6JFuzSAxe7uq+ddMzJOYiQlKea/5I39t5tBMSjbIBvS/mQst4NPM4NbkBpOoGudqt7pA8lC3oMH2E5GYThvRjHZKtWUTnOBjl2JHF2ww3mYuLeGOj8tW+t1H+yBw6f8Rlpuav6SeXU5MbQKpukKvd6s4RAVMO261Wq681Y0oP65BsTSwOyeq+tRtMt5kHynP15/W9nUMP8S3SG+4bMHMJrQB1oXeNUnWDXO1Wd47oMmXJiFqtZvrm2jok2xBXhORWw03mJULylb63BfT/GxVBKzKz+2Iu754MIFU3yNVudYfk4ZzcQZdemJhqtXq9IVueZEYvAkYoEsOw39V9a6Y3Yh76KynP1Z8bvlxCj6pZpHfgzlq5+D5yefdkAKm6Qa52qztHBEw5kWmtVrvFjCk9rEOyNbH8N8eOzN9uuMk8eOM0NkpfC1+uopOF59HdNSVmsv5IREnoXaNU3SBXu9Udkodq0x2mn1n99SYM6cc6JFsTy061vGfj2jjazTrFcutZ4csyerhZNNy3Qm+m3xmkNW3xoZwiVTfI1W51h+ThXNZh+rlo3mHCkH6sQ7I1sTgk66fKNcNN5iJC0moU7u97G+WPtNGOSR7uKCakYDwPKR9I1Q1ytVvdIXlwSKJz7zS814Qh/ViHZGtiKck+N9+8zHCTuXBIiqXOOfRmUF6gN8S3QD4O4AnpTJU4ll+k6ga52q1uelNhZJ0203fZvNKEIf1YhyRhWs3CA4abzIVD0umoJ9CTOUX1R5r0pr+e0fwRANVM24J0kKob5Gq3utE3WHlwSKYe9gt82YQh/ViHZGti6UYouN09hpvMRXTBcYIFtEe+hD4QovojkBMNkxFMe9DnFKm6Qa52qxvtkBifciQGonIL03ChCUP6sQ7J1sTSZROgTHvPJnas2Ol0VcDmGSa7zHz+CIDKRQTLPFJ1g1ztVjf6mpqHWY/bTH/dMO54WYckYbpdZXoj9nd9ZJagSwm9A/fnj8xoufhNzHD0Z1uk6ga52q1ufU47kZYhY9Bm+i6b4wbs2IR1SLYmlv/GLXT3G24yquWRaZRSe9G2zqELo1XQ/3Ee+lunwIax5SFVu9WNjvgeS8uSMdhg+mvcJSYM6UfoDjQSsXQlNOrufaabZPrx5LETdJ11dFKrgz4YdqVrUVI4QhP9ktc9N7f8cNLrHI7d5rLYpFsRwyzxMTFtWsLnjVjRh42QbE0sDsncQvNKw01ukP0Iiauc7m70kLj+cvEC+py7QodCStUNcrVb3TniNNNf/19jwpB+RoqQeJ73AHAK3d/f9n3/Ss/z3gX8U3QOw/3A63zfPz7kt9cA70F7Yzf7vv+7ZkyPnVictVPHKp8z3GSD7Ce17g66zveAH0DvQw4zXS6+H1dosSipukGudqubfN1kTXuNe7cRK/oYx6CrfN+/3Pf96A7/s8CzfN//IeA+4DcHf+B5XgH4A+BlwAHg5z3POzClzbFz6ODh/sRLo+zas/EKw03WyX6kawmn9QPo5NtouK9i9hNaseW0JSJVu9VNvkYNTnvdeJcRK/qYOIfE9/3P9L39MvAzQxZ7LvAd3/e/C+B53p8BPwUcnnS9CbGXmO7ejz+++CHDTW6Q/Vwgl27lCDpC0kBXbBVC5VTaFqSDVN0gV7vVLYlqtfpG022OeiELgM94nhcAh3zf/6OB738RGHahPRf4h773DwHPG7YCz/PeALwB4JJLLvnW1Vdf/X70hWsO+B7wfOAjwJuAtwE3V6vV19Zqtduq1er14VTIr0dP+PNedFnbL6OLt2yghyhdgk7EeQ063PSuarX6xlqtdnu1Wr0uej7nwuP//pHvLanl1caP1zfce4rlznndrjoJ4DjBcqtReKgy17709MnSlxZXGtecOLrw4d371m544sjirdHzyt61V584OvfRpdX6SzZOl75errQvanecx+YXmz/RWHf/R6ftHC2WOhdsrJfunl9qVE8enfvY6r71644dWbxtdd/aDceOLN66um/t+mNH5m9f3rNx7fqpcm1uvnlZq1l4oOB29wSoVrerNoJ29wcb9eJcsdw50Gq43y5VWpc165U7ynP15zU2Kl+Jnkvl+uXNhnu4WOo8rdVyHnXd7ko3UA0ARwXldts5Xix2z241C98tldsHmo3KXWe0Uak/p1kv3l0sty9uNZ0HXbe7v9tVJ5VDSREUul11ulAIdreahUeK5c7TW43iN4vl5vNajcoRpRo/FATlk1DfD6UWtCs6GaxT0ZUOAzccz9/Vr52m7pt11/UdSOWUnua7crLveRHc09CZB6cO3SKocAhe4IDT0m0U1qG9AJW1IW0sQXEN2lEbJVBtwIFA6ddBCQp1vZ7iGjQH21jWM362FvW6+jUF8+Ccmi1No2yniFnSNOp26i5B4fhsaRplO+kFZ0vTKNspmAPV0JpQ4KhSuXlZs1G5q1SuX9FsVO7U59/S3W6pdXG7WXjQLXb2dzrquOME8wGqE3RpOU6w2Gk7j7ul7rntpnt/sdy6tNWofKOvjSuajcqdxVL90lbTvc8tds5vt5zHCm53uRvoarGOCkqdtnPSLXb3t1uFB4ul9kWtZuWeM9oo15/dahRVrVb7beCDwFXAvWNec+8HLh71mlutVn9n2LW/HxUEO3d5eZ53ju/7j3iedxa6q+bXfN//QvjdbwFXAj/t+34w8LufBf533/dfH76/Hniu7/u/tuNKU+TQwcMvBX4d+Mu0bRmBs4CfA+5M25AtODt8fjo6b+Q0sEqYj5SWURaLxRIDe4FPpm3EiLz4xpsOrKZtRD8j9SH5vv9I+Px9tMf0XADP824AfhL4Z4POSMhDwFP73p8HPDKNwQmxjE7WNc7ufWs3GG5yjWznkCwACtp70PkjFXoz/QqgP1ogCam6Qa52q5t85ZBMRa1Wu910mzt22XietwA4vu+fCl+/FPg/w9EzvwFUfd/fKrv6a8AzPM+7EHgYfSdvfKhQDMTmkDxxZPFWw022iKnMvSEKwHKYid5ED/0VROVk2hakg1TdIFe71U2+RtlMRbVavc50m6PcWe8Hvuh53t3AV4GP+77/KeA/oSdK+6zneXd5nveHoLt3PM/7BIDv+23gV4FPA9/SH/n3mBYRA4vENEHSrj2nX2W4yRbZPQiiqbgXdc4IIKb+SERd6MgDqbpBrnarm3yVz58qmlOr1Q6aMiRipBwSaRw6ePi3gacAf2u67UKhvdDpuKYr+f0rdDQqa5yL7pq5GNoqjJKsoiMleTpwp6CtwBV4kEnVDXK1i9ddRhepNF1rKi6uuvGmAyuT/rhWq+2rVqtHDNqT6dyDNIktQrKw3HxBDM1m9SQwhy7aVoJOEZ0/0kWMMwI6u18iUnWDXO3idbvE1NWfUV5nukHrkAxnEYil6mB9w42jyyqrDkkB3a3n6KGvuIjb55w8TEUeA1J1g1zt4nVLc0iMjyYSdnEYmdgckmK5c14MzWbRIVlCF0FbBDrQLaP3N0HREdB1DyQiVTfI1S5edzRxaF6Y9rrxHCNW9GEdkuHMo4fTGicqsCaAVXTNkQWgDkohZv6afpQwByxCqm6Qq128boeYrhsxMe31/1EjVvRhHZLhlNGjQ/JCFiMkZXQOiQs0wCmh7RR60rJYLDOOg56ENi9krmaKdUiGE9v/4jhBHMWDsuiQgO62AejoUs/ZOwDiJxB6jEnVDXK1i9etgMfTtCRhzt55kfEQugPtSGwXzlaj8FAMzWbNIdmFTu5a5MmIiAoQMbvvIE6eIm0Gkaob5Gq3uoVxh+kGrUMynNgckspc+9IYms2aQ7KCTgqeRye2liEokq9uMEN0KzsvM4tI1Q1ytVvdwniZ6QatQzKc2Eqxnz5Z+lIMzWbNISmhk1ld9Pw1ZT0bpkQKsYzWyj5SdYNc7Va3MD5gukHrkCTM4krjmhiazZJD4qAjTMv0klgL0J1L1arUaC+kbUE6SNUNcrWL152l8/AoTNsT8FYjVvRhHZLhxNZlc+LowodjaDZLB8Iqeiz+PL2ckQKU8jQcziAVq1scUrWL152l83DsVKvVm0y3aR2S4cTWZbN739oNMTTbQpdozwLL6G6a+fA5jIzUhUZI7JTs8pCqXbxuUaMIa7Xa7abbtA5JwjxxZPHWGJoNE0czQRE93LeAdkhKgLJTk0tDqm6Qq93qzhlTRXSq1ep1pgyJsA7JcGL7X2KKkEQX/rRx0Tt533BfHP2Z+LsnYUjVDXK1i9edty6bqa5zNkKSHLGF3mKMkGShy2Y3erhv5JA4PFkuPrd3EVNidctDqnbxuvPmkEyFjZAkR2wOycretVfH0OwG2YiQLKHnqplD21RB/5cdqC9t98PZxeqWh1TtonWLm6erVqu9z3Sb1iEZTmz/y4mjcx+NodmsREhctFPioG3qc5KKQjPwrW55SNUuWncFXZ1aEm8x3aB1SAY4dPDwAjGWOF9arb8khmbraGcgTUroO4RFeqFLp/e6PZ+KValjdctDqnbRugvomzBJvNl0g9YhOZO9xOjpbpwufT2OZolxqPKI7EHnjyygHZMCev8KnTunnpZh6WJ1y0OqdtG6o1GFkvig6QatQ3Imu4lxzpVypX1RDM1ukP62XEAnslbYnD8S9qt2s5DjkgJWtzykahet20Gf9yRxlekG076IZZEVYnRI2h3nsRiaXSf9HJICepZfBx1hGrBHiUr46mF1y0OqdtG6FZC7UUaHDh5enOLn9xozJMQ6JGeyQowOieMEcVQsXSfdbTmH/s8G6490+5aRuq9Z3fKQql2ybgf4x7QNGZMAOGuK368YsuNJpO5A27FEjDkkiiCOSMY66Sa17gZOo8vFN0NbXDYNgwtElVXuYXXLQ6p2qztnBOicyUkxfnNtHZIzWSZGh6TTdo7G0GyLdIvyzKG7bMroTPMKOjrSFyERHc4ViFTdIFe71Z0zAvRkqJPyPVOGRFiH5EwWiTFbuljqXBBX2ylSQDtyDto5cjmjuFwgNOHN6paHVO2ideexSmuAvnmclOebMiTCOiRnskKM2dIb66W7Y2o6rQNiER0V6c8fcdmcPwIUpA2JC7G65SFVu9WdM7rogQiT8hFThkRYh+RMlonRIZlfalRjajoth2QVnT+ygO7qCmf3HSyj3BFaNMnqlodU7aJ15zFCotDXu0l5kylDIqxDciYL6CTRWDh5dO5jMTWd1gFRQeeOlNBdXZXQloEIieiy0gKRqhvkare6c0YHHdmelLeZMiTCOiRnskiMDsnqvnXjMySmjIMemRRFRQoMnZywKXRqcqtbHlK1W905ZBqH5GZjVoRYh+RMysTokBw7snhbTE2nESFZRkdFFunN/1Ng6FxA4qcmF4ZU3SBXu9WdM7pMUVCzWq2+1pwpGuuQnEms/8nqvrUbYmo6DYdkFVhDd3PV0c5cNNJmgHqe7yKmwOqWh1TtonXnMYekiz53T0StVjN+c20dkjOJtcjNsSOLt8bUdBoHRAm9QxfpOSTd4Yvm9i5iSqxueUjVLlr3Fue9TBOgi1lORLVavd6gLYB1SIYRq0Oyum/N+EZMiWjfWaKXxLrNjMOi754EIlU3yNUuWndsxTRjJGCKaqu1Wu0Wc6ZorENyJttcVKfn2JH522NqOukIyS70QRjN8gtb5o8AlITePVnd8pCqXazuAF2LKW9Es7NPyutNGRJhHZKEWd6zcW1MTafhkKzTq2xbQe9PW5RRbk2TzZ1jrG55SNUuVneJfDokHXS3+6S8w5QhEdYhOZNYu2zWT5VrMTWdtENSREdHXPTBuEP55EJsI5eyjdUtD6naxeqGGItpxkgHnfc3Ke81ZUiEdUjOJNYum7n55mUxNd1kiiFcY+KgHbeoXHyUP7KNQ9KZJjSYY6xueUjVLlV3twCcSNuKCdgh729HXmnKkAjrkCRMq1l4IKamGyTnkKyi7wii/BEHHSnZIn8EQOUx6csAVrc8pGoXq1uRT4ekzRRJrcCXTRkSYR2SM4n1Pym43T0xNd1guvDbOEQF0RbQjkm03m2m4Q7c2K3KJFa3PKRql6pbgZ7PK29MVRgNuNCUIRHWITmTWHNIAtSQomFGiCa2S4IierivQ6/+yA6oPBYOMoDVLQ+p2qXqjveaESNtprtmGM+bsQ7JmcS6c3W7Kq7kpzrJOCQuOldkkd5/5bBzUm0eCweZwOqWh1TtVne+aKHP55Ny3JAdT2IdkjOJ9T9xC939MTW9wXQ716isoof7LqB3aAedGLVNdw3IDeda3fKQql2q7lxHSKZJar3ElCER1iHp49DBwyW2TcycnkbdvS+upknGIVlG78jzaCcoSora4X9zhCa8Wd3ykKpdqu5czmMD+jw+TQ7J500ZEmEdks3sZcc7/emYW2heGVPTG8Q8ZDnEpZc/MsZQ467UIYFWtzikapeqO7fX0S7T3cS+xpQhEXn9I+NiN0NnqjXHqWOVz8XUdBIOSQm9E8/TuysYJX8EcIUWTbK65SFVu1TdKtab2BjZYDof4N2mDImwDslmVonZIdm1Z+MVMTW9Rvx1SFbRw9uW6CVEFRjpPxNbVtrqFodU7VJ1B0lEpuNimgjJu4xZEWIdks1EE8bFxvHHFz8UU9PrxL89l9C5InP06o9EM/3uQOVUnIZlF6tbHlK1S9Xt5jVCcpoprhnVavWNBm0BrEMyyC5ijpDs3rd2Q0xNrxN/l00B/R856P+pyMgZ5qKnJheIVN0gV7tU3a08R0gmHiFUq9WMz1xvHZLNLBNzhOSJI4u3xtR0rI4UejbfNr1y8aAdlBHH4FeETk1udctDqnaRuuegmOfcmWkiJNeZNASsQzLIIjrSEBsxRkjiZhWdpxLljxTR+8+I4Uqpd09Wtzykahep24VGUhWyM4WNkMTPMjFPIx1jhATiHQ+/gI6IlNBO2xj5IyD07gmrWyJStYvUXYDyE2kbMQUTd9nYCEn8rBCzQ7Jrz+lXxdh8nA5JAe2wRVGRMfJHAOpCM/CtbnlI1S5StwP1pCY1zRS1Wu2g6TatQ7KZ2Lts1o6XPxVj83E5JHPobpp5ehVZC4xV1dbN42yYBrC65SFVu0jdCtyH0zYiJd5pukHrkGxmnpgdkoXl5gvibD8mdqOHiC2iS9SX6I20GZHOfByGZR+rWx5StYvU7UBHZA4J8DrTDVqHZDNzxOyQ1Dfce2JsPs4ISRGdN7KBHnEz5gyXTt24VbnA6paHVO0idStwHknbiJT4pOkGrUOymdgnpyuWO+fF2HxcDkkBPboGtCMywbj7btxVZDOK1S0Pqdpl6i643TyPLppmpuLnGLMixDokm4l9GuluV8WZiR6HQ7KA7qYZrD8y5qzIasyIyqxgdctDqnaRuoNuoKTOcvyo6QatQ7KZ2B2SHBLljyyhHZMyY9UfsVgsFotlZ6xDspnYSwA7ThBneC+OCEmZXg5J//w1YxII3desbnlI1S5Tt6OCPCe1TnMTfrYxK0JE7kBp0moUHoqx+TgckgJ6dE1UBK0w2XqcuEvbZxSrWx5StYvUHXTajsSCcAB3mG7QOiSbib3LpjLXvjTG5k07JMtAHZ0/0l9/ZILumm7FmFW5wuqWh1TtMnW7xe7+tG2YgmmueS8zZkWIdUg2E3uXzemTpS/F2HxUQdUUq+hh0Ev0hvsqxk5oBSjkeQKqKbC65SFVu0zd7VbhwbRtSIkPmG7QOiQJs7jSuCbG5huYdUhK6OiIi46UTFEiub1gxqS8YXXLQ6p2mbqLpfZFaduQEm813aB1SDYT+/9x4ujCh2NsvsVUTsMmov8iyh+BifNHACpr05uUR6xueUjVLlJ3t9WsxFnsMrNUq9WbTLdpHZLNxJ5Dsnvf2g0xNh+VdTfBMtCkN3+Nw0T1RyJETk2O1S0RqdrF6S4D7VK5fkXahkzBxNe8Wq12u0lDwDokg8TukDxxZPHWGJuvY84h2YXOH4kmHIwiLxPWHxE5NTlWt0Skahen2wUazUblzrQNSYNqtXqd6TatQ7KZ2P+PmCMk65grf19CJ7MW6BVEmwJxd08hVrc8pGoXp9sFmjmPkExMHBGS2OduyQuHDh4uEd9cME8Sc4RkAzPb1EFHi+YHPpvi/xF39xRidctDqnZxuh2g3mxU/jZtQ9LARkjiZS8JlENf2bv26hib38DM0OVdYVuL6P+kED6mKHxUX9p5mVnE6paHVO3idDvAqWK5/uy0DUmDWq32PtNt2ghJj91MdcEdjRNH5z4aY/MbmHEyV4A1tENyCl1/BHqT601AUWIGPla3RKRqF6fbAU61GsU81yGZJm/yLcasCLERkh6r6FElsbK0Wn9JjM2byiEpoqMkDvo/MVDbpD2/8zKziNUtD6naxelWwONuqf30tA1JiTebbnCki5fneQ+g75Q7QNv3/Ss9z/tZ4O3AM4Hn+r7/9S1+exPwenT+wTeA1/m+X5/edOPsQidvxsrG6dLQ/8kQ60zfZeOit9U8vZwRh6miIwBOFrd5Aljd8pCqXabudtN5OG0bUuKDphscJ0Jyle/7l/u+f2X4/pvATwNf2OoHnuedC/wL4Erf95+Fvlj+3KTGxswuEuiyKVdireq3xvRRr1V0188S+v9ww8eU+TXdPM+IOQVWtzykapepu+B296ZtQ0pcZbrBicP7vu9/C8DzvFHWMed5Xgt91/3IpOuMmagQWKy0O85jMTbfYnqHZBnt2MwBJ9DDfbtMHSFRsScMZxOrWx5StYvTHQB0u0pa7kzEvaYbHPXiFQCf8TzvbzzPe8Oojfu+/zDwe8DfA48CJ3zf/8z4ZibCIrqwWKw4TjAX9zqmxKWXPxJN1meiYJzUfCWrWx5StYvUrRyj84fliRXTDY4aIXmh7/uPeJ53FvBZz/Pu9X1/y66aCM/zVoGfAi4EjgP/1fO863zfP6OgSujovAHgkksu+dbVV1/9frTgOeB7wPOBjwBvAt4G3FytVl9bq9Vuq1ar19dqtVvQuSrvAN4LvBL4crjujXD9lwCfB14DvBt4V7VafaMu8LLve/MLp1+4fnrhjl17Tr9q7Xj5UwvLzRfUN9x7iuXOed2uOgngOMFyq1F4qDLXvvT0ydKXFlca15w4uvDh3fvWbnjiyOKt0fPK3rVXnzg699Gl1fpLNk6Xvl6utC9qd5zHisXO0xaWN4qdtnO0WOpcsLFeunt+qVE9eXTuY6v71q87dmTxttV9azccO7J46+q+teuPHZm/fXnPxrXrp8q1ufnmZa1m4YGC290ToFrdrtpwC939jbp739xC88pTxyqfQwV7g65Lea7+vMZG5SvRc6lcv7zZcA8XS52ntVrOo67bXekGqgHgqKDcbjvHC4XuuZ12IYDuuVDqQn0/VDrQnINSXQ/rK67p5DWnrkO0qg04ECj9OihBoQ6deb1sc1nXJ+guAA1dPKl0ElqLenbQTgVUEwIXVAB09Wunqaczd9f1spVT+reVk33Pi+Ce1uty6tAtggojOYEDTku3UVjXE39V1oa0MYWm/ra20oQLzRnTNMp2Us3Z0zTqduqWoNWeLU2jbCfVnD1N220nRxXczjlKMe8UG3PtVuExt9Q5v90sfrtUbl7WbFTuKpXrVzQblTv1+bd0t1tqXdxuFh50i539nY467jjBfIDqBF1ajhMsdtrO426pe2676d5fLLcubTUq3+hr44pmo3JnsVS/tNV073OLnfPbLeexgttd7gaqGZ7LS522c9Itdve3W4UHi6X2Ra1m5Z4z2ijXn91qFO+BwP3Lj9SesbSbq9ARj5GvucC/qtVqnxj1mlutVn9niJuwCRUE49W68jzv7cCa7/u/F77/a+BfD0tqDRNfr/F9/5fC978APN/3/TeNtdIEOHTw8HvQEZIH4lxPZb5xfn29HOcwsX8J/M2Ev92PzvN5Wvi8gd5Bm0zdZdN0oSQtpIvVLRGp2sXp3gt8UjnNlaBbOp62MRPyT4Czb7zpwNjdTrVa7bnVavWrJo3ZMcTmed6C53lL0WvgpeiE1lH4e+D5nufNe56ngKuBb01qbMwsoEepxEqx1Lkg5lVMU202mtl3Dv1fVML3UzojoO+KJGJ1y0Oqdpm6C4VgJW0bpiAAzprwt883aQiM1ue3H/ii53l3A18FPu77/qc8z3ul53kPAT8KfNzzvE8DeJ53jud5nwDwff8rwJ8Dd6CH/DrAH5kWYYgFdDJnrGysl+6Oex1TEOWPKHT+iIuxCQcLIocEWt0Skapdpu52qxDnQIUkmHSU0EeMWsEIOSS+738XuGzI5x8ZZpDv+48AL+97/9vAb09nZiLMk0CEROeLuH8R4yomjZBU0E7IMr2ISAFde8YAnXlwBWajW93ykKpdnO4AQOeNuMZHnCREgC71MAlvAn7ToC0ys6K3wMQcMDty8ujcx2JexaQOySraIVtCF4groaMjhmqziCsrHWJ1y0OqdlG6C4S1mdrN4rdTtmUaAnRUfBLeZtIQsA5JP4n8F6v71o3PkDjApA7JAvo/KLM5f8QQTWlTk4dY3fKQql2U7gph3apSuXlGD0KOCNBR8Um42aQhYB2SfgzlSmzPsSOLtyWxngkooKMjCt1lU8DofyJuavIQq1seUrWL0v3k7OfNRuWudE2Zii4TRkiq1eprzZpiHZJ+Epn5eHXf2g0xr2KSqMYc+uBaoJczYjB/BPSYfolY3fKQql2U7gJhzmGpXL8iZVumIWDCyVNrtZrxm2vrkCTMsSOLt8a8ikkckt3AaXSEpInutnEwOrePqLunPqxueUjVLkq3g67TRLNRuTNlW6YhQJd7GJuwOJpRrEPSI5Eum9V9a8Y34gCTRkhK4WOd3vw1BhF199SH1S0PqdpF6VbASYBSuX55uqZMRYCOjI9NWB3dKNYh6ZHIKJtjR+bPKJtvmGj+mXGI8keglz9imJKku6c+rG55SNUuSrcDnAJoNjJdW2onuuiSF5PwepOGgHVIEmd5z8a1Ma+iwXgOyTy6a2aO2PJHQM81IRGrWx5StYvTfRrALbUuTtuQKZjGIXmHSUPAOiT9JPJfrJ8q12JeRYvxHJJVevkjdfRwtqhSq0EKsRedyyZWtzykapepu90sxDk3WdwETD6g470mDQHrkPSTSA7J3HzsY9ajomajMod2Qlx0/khM81F0KvG0m3WsbnlI1S5K95O5em6xsz9NQ6akgz7/T8IrTRoC1iHpJxGHpNUsPBDzKuqM51T01x+J3hssiBahp8eWh9UtD6naZerudNTxtG2YgoDJb0K/bNIQsA5JP4n8FwW3uyfmVawzepdNVCa+gvaUHXSkxHD+CECQSJ2X7GF1y0OqdlG6n7xpc5xg0hyMLNBBj6qchAtNGgLWIQHg0MHDiliiAmcSoAzW9hhKndH7BFfRY+mj/JFoxzScPwKgEvl/s4fVLQ+p2mXqDlAx3MAlRhfdbT8JGyYNAeuQROwllovwmXS7yvhGHKDO6MN2y+gx6C5654opfwQwXtckN1jd8pCqXaTuoGuygGTitJn8vH/coB2AdUgi9mC0KunWuIVu3AlQG4weIXHYPOQrpvwREBbO7cPqlodU7TJ1O06Q5+HOHSZ3SC4xaQhYhyRilXDmxrhp1N37Yl7FBqNFSJbR0ZR5evkjMdQfiXBEJrxZ3RKRql2U7idv3Dpt5/E0DZmSFpP7AZ83aQhYhyRiFwlFSOYWmlfGvIp1RouQrKIdkkV0HZJo6FdMXVddSUMC+7C65SFVuyjdT3ZPuaXuuWkaMiUtJq9D8hqThkBCM9zmgBUSyiE5dazyuZhXsc5oQ5hL6IhIAR0dijlT3BVZNMnqlohU7WJ0R9WtAWg33ftTtGVaukw42y/wbpOGgI2QRESz3MbOrj0br4h5FevsvF2j7+cHPosxS15cWekQq1seUrWL0V1AR5cBKJZbl6Zoy7S0mTww8S6ThoB1SCKiYa+xc/zxxQ8lsZ4dWEY7YEvonJEoUhJjt1XlVHxtZxmrWx5StYvRXUDXbwKg1ah8I0VbpqXBhA5JtVp9o2FbrEMSsosYxlQPY/e+tRsSWM1OkY5d9Lpp1ujlj8Q4bE/U1OR9WN3ykKpdjG6HvutFqVy/IkVbpqXNhFXKa7Wa8ZnrrUOiWSEhh+SJI4u3JrCanRySEjqZNYqKTNqHOAYVSVOT92F1y0OqdjG6HeDJaFCzUbkzRVum5TSTR0iuM2yLdUhCFtC5F7GTgQiJg/aI5/qWc4i9qJGYu6cBrG55SNUuRrcCjkVvch4hgQn9ABshiY/EHJKEIiTbEXVPLdJLaHKJfZSRmLunAaxueUjVLka3QkcWgNxHSGD0yt6bsBGS+JgnIYdk157Tr0pgNdtFSHahnY8of6SMjo7EHSGRkoE/gNUtD6naZeoulup5HmUzMbVa7aDpNq1Dool5hEmPtePlTyWwmu0ckhJ6lI2DdkyKTJjUNB7u6Z2XmUWsbnlI1S5G96bza6sZe/XtuJn0/P9Oo1ZgHZKIxP6HheXmCxJYzVYOSZS81J8/4pLIpFidPE/RPQVWtzykapep2y12zk/bhpR4nekGrUOiSSBCoKlvuPckta4hRPkj8+hhv1F0JIEqtU4idV6yh9UtD6naZeput5zH0rYhJT5pukHrkGgSc0iK5c55CaxmqwjJLnQ0ZJ7e/DUBiURIugkMLc4iVrc8pGoXo3vT+bXgdqWMLhrkOaYbtHPZaCbKMp6EblclkYm+XZdNCe2AdcL3CTljKgGnJ4tY3fKQql2m7m6g8j7L8aTXgEeNWoGNkMwqwxwSN/x8ru+zAtoxsVgsFoslVaxDokmsy8ZxgiTCe8OcjFX00OYFdP5IFClJZHQRBEL3NatbHlK1y9TtqKCUtg0pcbbpBkXuQENIrMum1Sg8lMBqmvTmp4mIhvrO0Zu/JsbZfQdxEnJ8sobVLQ+p2sXo3nTD12k7UgrCDXKH6QatQ5Iwlbl2EkV0WuiCZ/246Nl9oyTWAglGhqA76CAJweqWh1TtInTPMRBVdovd/SnZkjYvM92gdUg0iUVITp8sfSmB1dTZPGFeCe2ElOk5IQnnjxQSqYSbPaxueUjVLkL3GUU0263CgynZYopJb0w/YNQKrEMSkVikYHGlcU0Cqxl0SKL8kaXwuzJ62ycYYm0vJLeuLGF1y0OqdhG6XaDR/0Gx1L4oJVvS5q2mG7QOiSax/+HE0YUPJ7CayOmIWER79mW0YxLNX5MglbVk15cVrG55SNUuQreDPr8+SatZSbPYpQkmuiGvVqs3mTZEvENy6OBhRYLJnbv3rd2QwGpOs7nGTJQ/oujljySMmKnJB7C65SFVuwjdCjjR/0GpXL8iJVtSpVar3W66TfEOCbp6aWK5FE8cWbw1gdXU6TkdFXRp+DI9xyuF+iNipiYfwOqWh1TtInQ7DDgkzUblzpRsSZVqtXqd6TatQwJ7SWQuF01CEZINeg7JKtpBWQo/r5DY/DX9iLh7GoLVLQ+p2kXothGSEBshiYcVdN2OREgoQtLvkCygE1yL6K6clIr4iLh7GoLVLQ+p2mXqnoEIyaQ5JDZCEgO7STBasLJ37dUJrGadXg5JAZ3U6vS9T7AgWkR9Kfl1ZgGrWx5StcvUXSzXn522DWlQq9XeZ7pN65DoCqaNHZcyxImjcx9NYDXraMejQq9IWge9vQsk3l0DUJSQgT8Eq1seUrWL0H3GzVyrUcz7KJtJeYvpBq1DonMsEuuyWVqtvySB1ayjt+0q2tlaojfcN5rpN2Ha88mvMwtY3fKQql2mbrfUfnraNqTEm003aB2SXrGwRNg4Xfp6QqsK0PkjZXT+SJ3U8kcAnMT+42xhdctDqnYRus+IkLSbzsNpGGKQSQuDftCoFViHBPSw342kVlauJFrVrwDM09vhHFLJHwHoCp0R0+qWh1TtMnUX3O7etG1IiatMN2gdEp1DkphD0u44jyW0qmjOhQq9/BGXVLprAFQKeStZwOqWh1TtMnV3u0pC7sww7jXdoHVIdLdGYpNCOU4wl9CqymhHK8ofiWbiTOukIXVfs7rlIVW7SN3K2TRvmCRWTDcocgcaYJEEHRJFkNTO66C1ZSB/BCBIbALDbGF1y0OqdhG6z5gDTBGkMBVHJjB+c20dEv2nJuaQdNrO0YRWpdDaopyRFPNHQGo41+qWiFTtM6+7yJAIc7erErt+ZIzvmW7QOiR6J2sltrJS54IEVrOAzhWZD58L6G2dUv4IQCAy4c3qlohU7TOvu8SQa0WhEKwkb0omeL7pBq1DkvB/sLFeujuB1ayGz3PocvEpzV/TT0HCkMAhWN3ykKp95nW7DCmi2W4VkhqoEBvhrPfj8hHTdliHZPIx2BMxv9SoJrCaCtqbL6KLvmUg6aojsmiS1S0RqdpnXrfDkBGZbqlzfgq2mCQAzprgd28ybYh1SBJ2SE4enftYzKtQfY+M5I+AkLLSQ7C65SFV+8zrdtAR5020m8Vvp2CLSSZ1SN5m2hDrkPRmxU2E1X3rxmdIHGABHVaM5qxx0ds5sTyZ4TQlTE0+BKtbHlK1z7xuBzgx+GGp3LwsBVtMEtDr5h+Hm00bYh2ShDl2ZPG2mFexik5eLQJr6HokAUOGqyWLzKnJrW6JSNU+87oV8Pjgh81G5a7kTTHO2EN4q9Xqa00bYR2ShLtsVvet3RDzKsroYmgOOkJSJBPbuT7rd09bYHXLQ6p2mbpL5foVadswJV30FCpjUavVjN9cZ+BClTqJdtkcO7J4a4zNK/Q2jcrFE75POToCAu6etsDqlodU7TJ1NxuVO9O2wQBjO5PVavV600ZYhyRhVvetGd+IfSyhs8B3oR2SJXq5JCkj8+7J6paIVO0zr3vowIBSuX55wnaYZtIIyS2mDbEOScL/wbEj87fH2Pwu9EGzgi4XXyET+SMAJZF3T1a3RKRql6m72UiktlTcLE7wm9ebNsI6JAnnkCzv2bg2xuaj/BGFHmlTJGF9W9OaZIefAaxueUjVPvO6h0ZI3FLr4qQNMUwXXdV7XN5h2hDrkCR8wV4/Va7F1LSD1lKmN8y3QCaiIwAFofM9WN3ykKpdpu52s/Bg2jZMSYAuFzEu7zVtiHVIEv4P5uZjG7Me5Y+soof7BvRG2mSATiVtC9LB6paHVO0ydbvFzv60bZiSqEzEuLzStCGiHZJDBw+vkPCEc61m4YGYml5Bb89lesV7FJmJkKhm2hakg9UtD6naZerudNTxtG2Ykkm7bL5s2hDRDgmwl4QrmBbc7p6Ymi6iE5OiSfSaZGr7Bm7aFqSD1S0PqdpnXvfQHBLHCWZhDp9JolsXmjYiQxesVFhBX7gTI0DF4QD1549EOTEpz10ziMqYPUlhdctDqvaZ1l1gi2h6gEo0yh4DHSZzSM6YaHBapDsku0k4x6LbVcY3Inq4bx2t5zS6DHCDbDklGek6ShyrWx5Stc+y7gpb3LwG3bTnCZuaLnp2+HE5btgO8Q7JMvrCnRhuoRtHAtQu9LZcAo6hM6Y3yNT2nflw7hZY3fKQqn2mdRfYwiFxnCDvw507TFax/BLThszyDjQKyyScQ9Kou/fF0KxLL+TWRe9c9RjWMwWOyIQ3q1siUrXPtG6HLc6pnbZzxoR7OaPDBJPrAZ83bUiG7qBTYZGEIyRzC80rDTcZOZVR/kiUS5KR4b4RXZFDAq1uiUjVPtO6C+ju8DNwS91zE7bFNJN22bzGtCHSHZJVYkjM2Y5TxyqfM9xk1O20GziF9nQLJOxo7YwrsmiS1S0RqdpnWrcChpbGbzfd+xO2xTQdJnNI3m3aEOkOySIJOyS79my8wnSTaAdkAZ1ktIDewbpkKql15stKb4HVLQ+p2mdat4O+4TuDYrl1acK2mKbFZOkb7zJtyEhGeJ73AHpjdIC27/tXep73s8DbgWcCz/V9/+tb/HYFuBl4FvoC+Yu+7/+vqS03wyJwNMkVHn988UOGm+zPHwnQzkkGHc3K0IN59rG65SFV+8zrHtpl02pUvpG0IYZpMUFhtGq1+kbThoxz4brK9/3Lfd+PciC+Cfw08IUdfvce4FO+718CXAZ8a3wzY2MRSDTMuHvf2g0Gm3PRkZDIIXHQobcMDkOb+anJt8DqlodU7TJ1l8r1K9K2YUraTFA6vlarGZ+5fuJRNr7vfwvA87wtl/E8bxl4EfDa8DdNEi5EtgPz6HlfEuOJI4u3GmxuBf1/DuaPZLAvtyJyanKrWyJStc+07i27v5uNyp1JGhIDbSbwBarV6nWmDRnViAD4jOd5AXDI9/0/GvF3TwOOAB/wPO8y4G+AX/d9/4zQl+d5bwDeAHDJJZd86+qrr34/+oI7B3wPeD7wEeBNwNuAm6vV6mtrtdpt1Wr1+lqtdgvwevSUyO9FT/zzZXR52w10fsUl6KFKrwHeXSovP6fZKP/33fvWbnjiyOKt0fOuPadftXa8/KmF5eYL6hvuPcVy57xuV50EcJxgudUoPFSZa196+mTpS4srjWtOHF348GAbK3vXXn3i6NxHl1brL9k4Xfp6udK+qN1xHptfbP5EY939H522c7RY6lywsV66e36pUT15dO5jq/vWrzt2ZPG21X1rNxw7snjr6r61648dmb99ec/GteunyrW5+eZlrWbhgYLb3ROgWq16YUU53d3NurtYLLfLrYb7/YLb/sFOu3xHqVy/otmooO9aKiehvgjuaejMg1OHbhFUWMgocMBp6Sz5wjq0F6Cy1vfb6HkJimvQjtoogWoDDgRKvw5KUKjr9RTXoBn+dmMPzB3VbZVO6v7mwrqekEs1dQ0DFQBd/dppanvcdb1s5dQQe1LW1N/WVpqCeXBOzZamUbZTxCxpGnU7dZegcHy2NI2ynfSCs6Xpye2EW2o/o910Hi643b3drlpTDkVFUCgUg6d22+qxdqvwmFvqnN9uFr9dKjcvazYqd4Xn4TtL5frlzUbpbrfUurjdLDzoFjv7Ox113HGC+QDVCbq0HCdY7LSdx91S99x2072/WG5d2mpUvtHXxhXNRuXOYql+aavp3ucWO+e3W85jBbe73A30PEKOCkqdtnPSLXb3t1uFB4ul9kWtZuWeM9oo15/dahTvcUvtp7ebznFQC7Va7UVjXnPvBy4e9ZpbrVZ/Z/C6P4gKgp3zHj3PO8f3/Uc8zzsL+Czwa77vfyH87q+Bfz0sh8TzvCtDA1/o+/5XPM97D3DS9/237bjSBDh08PAngP83bTum4Bno7ppnAA8CZ6O7bKKoz8uAvI+Rt1gslrTZB3wibSNiYgH4oRtvOvCUtA0ZKYfE9/1Hwufvoz2m547Y/kPAQ77vfyV8/+fAc8Y1MkbUzouYZdee068y1FSR3iyNbXRXTUbzR0DfqUjE6paHVO0zrXvLO/diqZ73UTanmWAgRK1WO2jakB2N8DxvwfO8peg18FJ0QuuO+L7/j8A/eJ53cfjR1cDhCW2Ng8QdkrXj5U8ZamoFnT+ygh4fP4/enlnK0enDHZqhPvtY3fKQql2m7lYzlurbSTPJtfCdpo0YxSvaD3zR87y7ga8CH/d9/1Oe573S87yHgB8FPu553qdBd+94ntcf2vo14E88z/tb4HJiEDEFk9Tvn4qF5eYLDDW1jI6IzKMTWitkrvZIP51ZmKJ7AqxueUjVLlO3W+ycn7YNBpjkWvg600bsmNTq+/530cN1Bz//CLr7ZvDzR4CX972/CzBdLj231Dfceww15aCTjyLP1iWFiM/oOBmbWycprG55SNU+07q3vNFrt5zHkjQkQ3zSdIMZLKCVKIlfwIvlznkGmimji9Qt0MsfKXNmd02GoiXdsce5zwZWtzykap9p3d2tvii43VmovzLJtdB4Pqh0hyTxLpto+PCUrKCTV1fplYtXnDl/TYYcErXlAT3bWN3ykKp9ZnXPs81ggWjIrUAeNd2gdIckryyhc0ZK6AzpOYZ78LN6grBYLJakyOBkpbOJdIckcf2OE5gI7zn0RtXA1rlAbbSzkgECofua1S0PqdpnVncB2DI/xlHBJDPlzgJnm25wVnegUUk8h6TVKDw0ZRMVdPhwEZ0z4oaPYWHDDpnZxk5G66PEjdUtD6naZ1a3wzazwnfaziyUzJ/kWniHaSMycrFKjcQdkspce9oiOqtoR2MXm/NHhjkkLVLIkxlOt7LzMrOI1S0PqdpnVreDLq0wFLfY3Z+gLVniZaYblO6QJK7/9MnSl6ZsYgHdXVNCe+0Vtk5eDeeQyAKFDE74lwRWtzykap9Z3Qo4ttWX7VbhwQRtyRIfMN1gRi5WyXPo4OEFdKQhURZXGtdM2YSDHuK7U/4I6KhJRmqTtBfStiAdrG55SNU+s7oVevDAUIql9kUJ2pIl3mq6QbEOCbCXFMqsnzi68OEpfh4NP1tGZ30X2Tp/BDLlkFTWdl5mFrG65SFVu0zdrWbFVLHLNBn7OlGtVm8ybYRkh2Q3KUxEt3vf2g1T/HwFHdVZQYcQF9k6fwR0l05GtnF9FooHTYDVLQ+p2mdW97b1nErl+hVJGZIlarXa7abbzMjFKhVWSMEheeLI4q1T/HwenUPioiMkJbY/WDJUyrkyC5noE2B1y0Oqdpm6m43KnWnbYIBJIiTXmTbCOiQJM2WERKGdkP75a7bjFJnZxjN797QDVrc8pGqXqdtGSMyRkYtVKiyhR6EkyhQRkgW0A7WC7oopo3NItqsgmKHpwGXePVndEpGqfWZ1b9tlMyMRkrGxERKzrJBCl8bK3rVXT/jTVfSBsQudP7IQvs9JMaL6UtoWpIPVLQ+p2mXqLpbrz07bhjSo1WrvM93mTiH/WWaVFBySE0fnPjrhT6O5axy0E1IkMyNoRqEoMgPf6paIVO0ydbcaxVkYZTMJbzHdoOQIySLblAOOi6XV+ksm+JnizPyRIqNNnpeRGX/b82lbkA5Wtzykapep2y21n562DQaY5Ob2zaaNsA5JwmycLn19gp8toIf2rgLr6PyR7eqPZBAnQyN+ksTqlodU7TOre9sCmu2m83BShmSMD5puULpDsuX8BHFRrkxU1W8V7cEu06s/kqP8EYCu0BkxrW55SNU+k7rL7DD4oeB29yZkS9a4ynSDknNI5oHE+zzbHeexCX5WDh8K7a275Cp/BEAlPqIpG1jd8pCqfSZ1F9khEt3tKpG5M8C9phuUHCEpk0KEwXGCuTF/otDbKUpoJXw9Sv5IlpC6r1nd8pCqfRZ1F9i+tALKoZiQLXEyyQ3uimkjZnEHGpVUtCuCcXfeJfQBsRsd0akwwkGSPYKcRXRMYXXLQ6r2mdRdYIfRmIqgkJAtWWPcm+sdkeyQpHLwdNrO0TF/sgtt6xKb80dGDY9mZJTNTIZzR8DqlodU7TOpW7FD1363q9YTsiVrfM90g9YhSZhiqXPBmD8po50Q0N00BXKXPwIQzGLC2whY3fKQqn0mdTvAE9stUCgEK8mYkjmeb7pByQ5JKmG2jfXS3WMs7qCdj3L4PqpFkrf8EaAwq0MCd8DqlodU7TOpWwEntlug3SpMMlBhFviI6QYlOySpML/UqI6xeJQ/sooOG84xQp9mNumILJpkdUtEqnaZut1S5/y0bUiJN5luULJDkkq3x8mjcx8bY/EV9BDfReA4vflrti3UM0BGckhklpW2uiUiVbtM3e1m8dtp22CASa6HbzNthGSHJJUum9V96+PMkFhEOyGgHYuoC2ccMuKQNEVOTW51S0Sq9pnUveP5s1RuXpaEITETHDp4eHHnxTZxs2kjJDskqXDsyOJtIy66Vf7IONERyEy+ycxOTb4DVrc8pGqXqbvZqNyVtg0GCICzxvlBtVp9rWkjJDskqWhf3bd2w4iL7kJXCNwNnERXlp2k/kibGMaLj099Fu+eRsDqlodU7TOpe4QISf2KJAyJmQAYqwR+rVYb9eZ6ZCQ7JKnkkBw7snjriIvuQjsg8+gs73l0tGPcCEmHTGxnmXdPVrdEpGqXqbvZqNyZtg0GCNBFN0emWq1eb9qIDFyoUiMVh2R139qoG9FF549EdhaYbHu1JvydYWby7mkErG55SNUuU3epXL88bRsMsTrOwrVa7RbTBmTgQpUaqWg/dmT+9hEWiyY9nENHRaK5bCaphNgmE9u5JPLuyeqWiFTtMnU3G2PVlsoqXXRUfhxeb9qIDFyokufQwcOTJIcaYXnPxrWjLIZ2JPrzRxwmm7+mSSa2c2vcDO4ZweqWh1TtM6l7x0EBbql1cRKGJMC4Ea53mDYgAxeqVNjLZNGGqVk/Va6NsNgu9JDfCtohiSIlk4yYaZKJUvMFofM9WN3ykKp95nQXGeE60W4WHkzAlrjp0BvROSrvNW2EVIdkNynNljs3P9KYdRfthEybPwK6qmsGtnNnrISp2cHqlodU7TOnu4TOwdsWt9jZn4AtSTBuhOuVpg3IwIUqFfYxwo4WB61m4YEdFnHRkZB5tNdaQHvqk9q7MeHvDKOaaVuQDla3PKRqnzndLiPcuHY66nj8psROgL7mjMOXTRsh1SFZISWHpOB29+ywyAraEVmhN9x3kvojEafIRJdN4O68zCxidctDqvaZ0+0wwrxhjhPMwhw+Ab2q4KNyoWkjpDokS+jcisQJUDs5QsvoUGEFPaFeBR0xmbQE/Gky4ZCojJSwTxqrWx5Stc+cbgd9Dt6WAJXKAAnDdBg/QmI8+i7VIVkmJYek21U7bcQCOn8k2jYumXAopiYjJewTx+qWh1Tts6ZboSPM2xJ004m2GyZg/Irex00bIdUhWWSEUFwcuIXudglQRfRBvYDO7o7yR1IZEWSWmQvnjojVLQ+p2mdOtwN8f8eFnGAWhjt3GbNSK3CJaSOkOiS7SCnZs1F379vm6xV6BWqOoR2TSeuP9JOBUKozawlvI2J1y0Oqdpm6O23n8bRtMEAHnSowDp83bYRUh2SRlBySuYXmldt8vYweC14G1tEhtIBMOBTT0p21IYEjYnXLQ6p2mbrdUvfctG0wQIdehfBReY1pIyQ7JKkU8Tl1rPK5bb520M5ItF0K8VuUFO6sFU0aEatbHlK1z5zukW4E2033/rgNSYCA8SMk7zZthHVIEmbXno1XbPFVGe2lLqGHJLvhYxYSppjRstIjYHXLQ6p2mbqL5daladtggA7jJ7W+y7QRUh2SOVJySI4/vvihLb5aoZc/8gS9/JEZ6Zet7JitPptY3fKQqn3mdI8UIWk1Kt+I25AEaKMHUIxMtVp9o2kjpDokqenevW/thi2+WkI7SkX0CKAKM5M/AlKnJre6JSJVu0zdpXL9irRtMECbMbtsarXaKDPXj4VUhyS1uh5PHFm8dYuvHPQOEeWNmMwfyYBTUxE5NbnVLRGp2mXqbjYqd6ZtgwFajOkPVKvV60wbYR2ShNkiQlJB7xDL6NE/RaabvyaDyLx7srolIlW7TN0zEiHpMmaXjY2QmCO10StbREhWw+cVevVHIKUZieNB5t2T1S0RqdpnTvdIkeUZiZBMkkNiIyR5Z9ee068a8vECOn+kgE5iLSdqVCLURWbgW90Skap95nSPNEdNsVSfhVE2bcasQ1Kr1Q6aNkKqQ5Jal83a8fKnhnzsoL3T/vlrZgz3dNoWpIPVLQ+p2mdK9wIjdpm3mttW384LG4zvD7zTtBFSHZLUumwWlpsvGPhoHr3jr6B3ihLaITE53DcDSa2dWZiiewKsbnlI1T5TuguM2GXuFjvnx2xLUox7I/w60wZIdUhSo77h3jPw0Ur4vIyePXERHcEx6ZBkYBZOJ5XJDNPH6paHVO0zpdthxAlY2y3nsZhtSYLTjO8PfNK0EVIdktR0F8ud8wY+mg8fUf5IEfMRjTap56V0x0qYmh2sbnlI1T5Tuh1GLJ5ZcLuzMrpo3FSG55g2QKpDkloOSberBjPRFTpUFuf8NWMnLJlHZSBKkwZWtzykap8p3Q5wYpQFu4GakWraY/sDj6ZtwKyQmkMyQJQ4tRsdMitjPn8EtEMidVtbLBbLuDjArJXC34nUr4tSL1Kp6XacoD+8t4reCaL8kQXM549AJhySQOi+ZnXLQ6r2mdM90qghRwXjzpI7K5xtusEZHF66PYcOHi4x4vjyOGg1Cg/1va2goyKKCQrTjLNaUvd+nRmqOjsOVrc8pGqXqbvTdmatINyo3GG6wVnzaEdhL/rinwqVuXZUREeFjzLx1x9pkvq27lbSXX9aWN3ykKp9pnSPPLDALXb3x2lIhnmZ6QYlOiS7SXGOmNMnS18KXy6GdqwCa2jHZOSx72NSJ/UISWGkjPXZw+qWh1TtMnW3W4UH07YhJT5gukGJDskqKTokiyuNa8KXK2gnYQk9f01UdjkO2zZI3SFpL+y8zCxidctDqnaZuoul9kVp22CIca8RbzVtgLgcEmAX5pNGR+bE0YUPhy/L6BwS0IXL+of+Gl8tqTsklbV0158WVrc8pGqfKd0jd9m0mpXBYpciqFarN5luU2KEZBcpRkh271u7gV7+SImeo1Akvoqqp0ndIZE5NbnVLRGp2mXqLpXrV6RtQxrUarXbTbcp0SFZJsUIyRNHFm9Fd9M06eWPVOhVap1RZm5q8hGxuuUhVbtM3c1G5c60bUiDarV6nek2JToki4w4R0EchBGSFfR/v0Sv/khAipGb+JF592R1S0Sq9pnSPXKXzQxFSMaKotsIiRmWSdEhCSMkJXpJrHHnj0SkPOOvzLsnq1siUrXPlO6Ru89thMQcUh2S1Ian7dqz9nP06o8Qvo4zfyQj1JfStiAdrG55SNU+M7rnGSNaXSzXnx2jLUkyboTkfaYNkOiQLJCiQ3Li6Nzn0LkiK/TyRxziqT+SIYqzlIE/Bla3PKRqnxndY9WDajWKIkfZAG8x3aBUhyS1A2duqXk1eoeP6o9E89ekVj02GdrzaVuQDla3PKRqnxndYzkkbqn99BhtyTJvNt2gRIdknhQjJI314nfQTkhEgdSH5CaBk1reTrpY3fKQqn1mdDuMcY1oN52HY7Qly3zQdIMSHZJCmusulTrn0SuIJiR/BKArdEZMq1seUrXPjG4HODXqwgW3uzdGW5Jk3Bvjq0wbMFKlVs/zHkBvoA7Q9n3/Ss/zfhZ4O/BM4Lm+7399m98XgK8DD/u+/5PTGj0laTphy+2OOoKeT+cEMId2kDYSWHfKo2zUjHdJbYXVLQ+p2mdGt0J3p49Et6tmJXdmXO413eA4peOv8n3/8b733wR+Gjg0wm9/HfgWeoRL2qTZPbKCYg7dbfT3wD60o9BJ0aakkBiNA6tbIlK1z4puha5uPdrCDsVAQIx7CCumG5x4B/J9/1u+7397p+U8zzsP+Ang5knXZZg05+8pqDPzR2blIN6BQECezDCsbnlI1S5TtyJIMw0gTeZMNzjqxTAAPuN53t94nveGMdfxH4F/g4g8iW1xAVot1UL/Fyr8TEJ0hBkK546J1S0PqdpnRvdY3dvdrkptkETKfM90g6NGC17o+/4jnuedBXzW87x7fd//wk4/8jzvJ4Hv+77/N57nvXiHZd8AvAHgkksu+dbVV1/9fnRIaA4t/PnAR4A3AW8Dbq5Wq6+t1Wq3VavV62u12i3A64F3AO8FXgl8GbgQnaNxHLgEtXthebX+slPHKp/btWfjFccfX/zQ7n1rNzxxZPHW6HnXntOvWjte/tTCcvMF9Q33nmK5c163q04COE6w3GoUHqrMtS89fbL0pcWVxjUnji58eLCNlb1rrz5xdO6jS6v1l2ycLn29UOxc0W46p5UKDnQ7jUe7XXWeUvxAt+M85JY6l7SbxW+Xys3Lmo3KXaVy/Ypmo3JnqVy/vNko3e2WWhe3m4UH3WJnf6ejjjtOMB+gOkGXluMEi52287hb6p7bbrr3F8utS1uNyjf62rhCVxJsh6N5OvM6G75bBBU6iYEDTgu6FSis62nEK2u6FHTlZN/zkq410I7aKIUnIUffHak2BCUo1PV6imvQDH/b2QUc1W2VTkJrUa+rUwHVhMAFFQBd/dppanvcdb1s5dQQexbBPZ2epv62ttIUFKFZny1No2wn1Yb6/GxpGnU7dSvQOj1bmkbZTqoN9bn8a2o74BKdO4vl+rNbjeI9bqn99HbTebjgdvd2u2pNORQVQcFxgl1B0FhptwqPuaXO+XGfy4ul+qWtpnufW+yc3245jxXc7nI3UE0ARwWlTts56Ra7+9utwoPFUvuiVrNyzxltDNHUaRcKtVrt2lGvucD/U6vVXrjjNRc+D7ymWq3+znY+AIAKgvFyHT3Pezuw5vv+74Xv/xr418OSWj3P+/8B16NrbFTQOST/zfd94yVnR+XQwcOfB/wUVn0h4EDnR6Dwd+j8kTnG6KuckmuAowmtawhtB1yBUTKrWx5Stc+M7n3AJ0ZfvD0HbhIDE+Lmx2+86cDizotparXaU6vV6j+YNGDHLhvP8xY8z1uKXgMvRSe07ojv+7/p+/55vu9fAPwc8Lk0nZGUcYGFQrFzXvg+6fyRlE8UnVkpmjQmVrc8pGqfGd1j3aW7pc75cRmScd5kusFRLoj7gS96nnc38FXg477vf8rzvFd6nvcQ8KPAxz3P+zSA53nneJ43hneZOGkkkUa1RuY7Lfe+0IYiyeaPdOjNn5MCM1NWekysbnlI1S5Td7tZ3HFwx4zyNtMNjt1lk3cOHTz8BeBPE17tPvRQ3wOlufqB5kblb9C1SDZIzil5cficVBfRAFGfrjSsbnlI1T4zuvcCnxx1YZ0fUrkrPnMSY9wum1uq1eprTRogZMjpJtIYmraMjk5UmhuVv0bn0yRdfyRMLEuLmThRTYDVLQ+p2mXqnhFnZGxMOyMg0yFJQ7ODTmBVlYX6i0in/kgrhXX2Uc9CUbwUsLrlIVX7zOge60axVK5fEZchWaZWq91muk1RXTaHDh4uAZ8D/iTB1ZaBc4GnAHuAx4Bz0BGSJDOzr0CHIo8nuE6LxWLJE2Uguk5I48eBp9x404HUcoGkRUh2obsukl5nB11T5Vix0ngR+n8feXprQzSxEZIUsLrlIVX7TOh20efKkSmV65fHY0riBMBZoy4c1v4yijSHZA+66yJJltE5IyVgo1Uv3oPe8EkPw62T6jw+JZH9y1a3RKRqnwndLmPeLDYbpbtjsiVpAnQUfVReb9oAaQ7JKmN6vwZQaIfEASgUW88knf99g1QdktbI2duzhdUtD6naZ0K3g755Gxm31Lo4JluSJkBfI0flHaYNkOaQJB0hqaC7a5bQjlCh0yo8TPLdRqCH+6bokBSEzvdgdctDqvaZ0O0AY+VQtJuFB2OyJWkCdIrBqLzXtAHSHJKkc0hW0A7JLuAYsOAUumeRfP4IwAlSdUg6lfTWnSZWtzykap8J3Q5wapwfuMXO/phsSZoAnWIwKq80bYA0h2SZZLtsFtEF0YroMOBcOEnfLMz3MCYq6a6yjGB1y0Oq9pnQrYDHx/lBp6OOx2NK4nTRUf1R+bJpA6Q5JIuM2T84JQ46mbUQvXcc5hJcf4YIRp1ZesawuuUhVbtM3Y4TzNIcPuPkAV1o2gBpDskyydX+mEfnqyyjnSAXKAQEaeSPRKRYdEbJKXizCatbHlK1y9QdoJKsuB034zgkxq+lEh2SpIq+rKBDYE/mjwCFoJvWXDKpI7CbCrC6JSJV+yzoHtupCrqJl5KIiy76OjUqx00bIC3EtgA8kdC6otyRaFz7ChAUCsFCp83RhGzIEIFL8jVgMoDVLQ+p2mXqdpxgsdM1el1ZCB8l9DWkgr6OFOlNO1IIHyp87/S9hp5jNczBGlYHKwjb+9QYdl4CfGGM5XdEokOSVIREMZA/AqhO2xkrYWp2cGYh4W0CrG55SNU+M7ovRJ+7S+hrpIs+nxfpXfQVoUPQabsFIKpFErC1EwDDHYH+z7top66FHoDRQI/6WQPWw8cp9KjJ48BR4PvAPwIPAd9PsPT75003KM0hifI64mYBPbx4Bd3PFkVK2m6p+9R2k79LwIaM0a2AxO4qq1seUrVPrLuAvvj33+X33+33fxYwPNVgp66W7b4P+h4O+ny9ge5qjxyBNXpOwDHgCHpeskeeevETP/Pyl1d/e4f1zyKvAX7HZIPSHJICyTgkq2hPdwntwc6jD6ZGu+nen8D6M4g7C0WTJsDqlkcutBfo3f3D5rD/4IOB18MIwqYGR5xs5Qj0Rwa66HpNnb7X7b7n/mhBE+30xOXwVW+86cCl4/ygVqt9LyZbss67TTcozSFJKok3mrsmmqhpldC7L5Zbz2o1Ct9IyI5BUsyCby1CYayCQ7OB1S2PbbWXGd73v5MTED1v1SWg2D5voP/zbvg6cgCii3+/I9AMn9fDzzbC5y0dgWK59ewUz22mmCQx913AG00bkgOM65bmkCRRqTQ6iRTpOUBRHgmtRiXvB+yEVIRenKzuDFDkzP7/6NgcfN4qEhD0fd//fsj3FdCOR+Qk9D/6owDRRb/d92ihywQ06TkCuej+mYFzW9TVPhbValWiMxKLbuuQmGcBfTI5m17+yJNdRaVy/Ypmo3JnAnZkjPoyVGZhNtAxEa27xeR5Af0X8+3YLhoQPfq7BPq7AqJHk81dAlM5AVKP8RnQvcQEgx5qtdrt1Wr1uhjsyTRx6JbmkBR2XmRqou6ZXejEpwX0yXYDIOcH7BRIvChDCrrLTJ8X0H+RHxYRGKFLoNJFH28B2hnvdwiSzgtIFKnH+AzojmpGjYVEZwTi0S3NIUmCMr1+4lb4GsKT9QzcRUzITEcKBvMC+p4bLpS77JwXwJDXwwgGfjdJXkDUJbBjXsCkyN3P5WqfAd3zwNgJqjZCYg5pDkncXTbRhchlc/7Ik+tN+YBNManVuDMyaV7ATk7AKN+PkRdQjqIBLXoOQG7yAiYl5xemqZCqfQZ0l9BDecdCojMCNkJigri7bJbQF6Rz2Fx/5MmCQcVS/dJWs3JPzHZsxbSlnaeoF9BW4e4WXcT7L/zD7viHfd8lpbyASQm3931JrzdtUt7PU0Wq9hnQXQYeHPdHtVrtYLVavSkGezJNHLolOiQ/T++OuX84WyN8vY6+cEWFcMapKbCr7/kx9ERFij6HpNV007w4rQP7Ge6YDDoGDLyOnIE2k+UFFBFYVjrl7Z0aUnWDXO0zoLsA3DXB795p2I68YFy3KIfkxpsOvCh6fejgYYUeCbMf2AfsRieknoWusHoO8Exgjt5df/R/RWWD+0sIB+gE1mLYztPC1wFwEl3u96hb7JzfbrnfiVHmdnwxpfWSsu7UsLrlIVX7DOguAF+e4HevA/5vw7bkAeO6RTkk/dx404EAeCR8TM2hg4dLwFOBHwAOABf0vd8Vvr6o3XJLwLlszmnoH/IYRW+ibon+7ogmOiExmtvguAnbk6Ddcsbum50FrG55SNU+A7rVjTcdmETDJ41bkg+M6xbrkJjmxpsONIH7w8eWkw7VarUbqtXqrTu1d+jg4f3A5cAz6Dk2u8PHIvAUdAWmKEIT0f+63+mJki/b9CZwirqpou6p2HIsCm53udMmS8WyEsHqlodU7TOge9Icu+cAeS8KNwnGdVuHJHkeHWWh0FP/dPiYmkMHDz8NeDbwg+jZLPeju5aix3n0ZifeLnoTPUdJpZFz00ZHb6IIzin6HJxuoGZlJtCxsLrlIVW7VN2MeE6fQYzrtg6JEG686cB3ge+aaCvMv3kWOsfmmeh8myj3Zjfa2TmfvhE5Qccpop2WYd1TsHnobH9ybJRkPDOFsywWi8VyJtYhSZ6z0zZgWsL8m28wRrhuu66qQwcPL6Lzbn4InXtzDr1E4xV0zs0cen8drC462EXVn38T5eA06c0Rchrt5Dw+qu3T4Kig1EliRRlDqm6Qq30GdE9qfu7P6RNiXLd1SJLnjrQNSIktdd9404E14KvhY2pCB+eF9PJvzkE7N/vQtWLOAi6mN1IqYrv8m8i5iR4Neg7OKXQOzhl02s6sVqfdFqm6Qa52qbqx53RjWIckeV6GzASoxHSHDk5c+TfnoROK96GjN7vQdwr95eMjHFAldHRmMP+mfwRVlGBcJ4EE4yRwi9397VauExwnRqr2GdA9aSVre043hHVIkucDaRuQErnVPU3+Ta1W21etVo9E7wfyb56GdnDOQScW7wb2hJ+V6dW5ieifKI++5/4E46h7qkkv/2ZTgnEStFuFsStezgpStUvVTY7PbVNiXLd1SJLnrYC4MsNY3cBk+TfbETo4PwJcgh49dT6wF+3YLKOjOU+jl38Dw2f/7Xd8+ifoi+bhiRKMTzFC/k2x1L6o1XTzXEZ8YqRqnwHdk1aStuc2Q6ggSHG+NYvFkinC/JvnABehu6jOpufgLKGrEc+xuVIxbJ7HiL7nnQr8bZl/Y7EkyFnAOTfedODCtA2RjHVIEsZOVS0L6br7Cvxdik4wfgqbC/wtcmaBv8EJG/ujN6kW+BuFUrl+xQzMfDs2Odf9dCC48aYDl437Q+nHuMk2rUNisVhyQ5hg/MPo4eHns7nAXxTB6S/wN9glNejgjFXgzzKz/BDwvRtvOvBP0jZEMjaHJGGsNy0Lq9ssMRT4+2H08PCLOLPA31PQjk+RzdGbYU5ORAdaZSiuIazAX84jJHNMWHnUHuPmsBESi8ViMUCYf/NsdPfUBZxZ4C/KvykyPGqzVf5Nf4G/Nrp7KtECfwL4MeDQjTcd+D/SNkQyNkKSMLVa7X3VavWNaduRNFa3LCTqDuvf/K9arfYLJrQfOnh4GfhRzizwdxY69yYq8DdYwXi7/JvBBOOoe2rbAn+jUCzXn91qVPJaj6MA3DfJDyXu6xCPbuuQJM9b0jYgJaxuWUjVDYa033jTgZMkU+Bvlc0F/qLrwij5N5GD0201ii10LkYeC/y5wBcn/K3Ufd24buuQJM+bgd9J24gUsLplIVU3ZFR7jBNsPg04zy21X9xuFo5wZoG/Yfk3sH2CcdIF/pzw/5mETG7vBDCu2zokyfPBtA1ICatbFlJ1gwDtwwr81Wq1Z1Sr1b+bpL0tCvxF+Tf9Bf5Gzb+B8Qr8TZNMOfPbewuM67YOSfJcBUx00OYcq1sWUnWDXO0T6w4dHNMTbEYF/p5BL/8mKvB3AZvzb9anWJ3d3oawDkny3Ju2ASlhdctCqm6Qqz0zusME4y+Ej7jJjO6EMa7b2XkRi2FW0jYgJVbSNiAlVtI2ICVW0jYgRVbSNiAlVtI2ICVW0jYgJVZMN2gdkuSZS9uAlLC6ZSFVN8jVbnXLwrhu65Akz/fSNiAlrG5ZSNUNcrVb3bIwrts6JMnz/LQNSAmrWxZSdYNc7Va3LIzrtg5J8nwkbQNSwuqWhVTdIFe71S0L47qtQ5I8b0rbgJSwumUhVTfI1W51y8K4buuQJM/b0jYgJaxuWUjVDXK1W92yMK7bOiTJc3PaBqSE1S0LqbpBrnarWxbGdasgmKZirsVisVgsFsv02AhJwtRqtdvStiENrG5ZSNUNcrVb3bKIQ7eNkFgsFovFYkkdGyFJmFqtdkvaNqSB1S0LqbpBrnarWxZx6LYOSfK8Pm0DUsLqloVU3SBXu9UtC+O6rUOSPO9I24CUsLplIVU3yNVudcvCuG7rkCTPe9M2ICWsbllI1Q1ytVvdsjCu2zokyfPKtA1ICatbFlJ1g1ztVrcsjOu2DknyfDltA1LC6paFVN0gV7vVLQvjuq1DkjwXpm1ASljdspCqG+Rqt7plYVy3dUiSZyNtA1LC6paFVN0gV7vVLQvjuq1DkjzH0zYgJY6nbUBKHE/bgJQ4nrYBKXI8bQNS4njaBqTE8bQNSInjphu0DknC/MEf/MEladuQBla3LKTqBrnarW5ZxKHbOiTJ84a0DUgJq1sWUnWDXO1WtyyM67YOicVisVgsltSxDonFYrFYLJbUsQ5J8vxR2gakhNUtC6m6Qa52q1sWxnWrIAhMt2mxWCwWi8UyFjZCYrFYLBaLJXXctA3IK57n/RfgJ4Hv+77/rPCzdwH/FGgC9wOv833/+JDfXgO8BygAN/u+/7vh57uBDwEXAA8Anu/7x+LWMi6Tavc876nAHwNPAbrAH/m+/57wu7cDvwwcCRd/q+/7n4hdzBhMuc0fAE4BHaDt+/6V4eeZ3+ZTbO+L0doingb8O9/3/2OOt/c7gJ9C77/fB17r+/4jQ347i8f4jtpn9BgfdZs/wGwd46Nsb6PHuI2QTM4twDUDn30WeJbv+z8E3Af85uCPPM8rAH8AvAw4APy853kHwq//LfBXvu8/A/ir8H0WuYUJtANt4M2+7z8TeD7wz/u0Axz0ff/y8JGpE1XILUymO+KqUNuVfZ/lYZvfwgS6fd//drQ9gR8G1oGP9C2Sx+39Lt/3fyjU9JfAvxv80Qwf4ztqZzaP8VF0R8zSMb6jbtPHuHVIJsT3/S8ATwx89hnf99vh2y8D5w356XOB7/i+/13f95vAn6G9UMLnW8PXtwKvMG23CSbV7vv+o77v3xG+PgV8Czg3ZnONMcU2347Mb3NDuq8G7vd9/8EYTIyFLXSf7Hu7AAxLwpvVY3xH7TN6jI+yzbcj89vckO6pj3HbZRMfv8jmUFbEucA/9L1/CHhe+Hq/7/uPgj6wPc87K14TY2Mr7U/ied4FwBXAV/o+/lXP834B+Dr6LitTYc0R2E53AHzG87wAOOT7fpShPgvbfMftDfwc8KcDn+Vye3ue9x+AXwBOAFcNWWRmj/ERtPcvewEzcoyPqHvmjvFxtjcGjnEbIYkBz/N+Cx26/JMhX6shn83MUKcdtEfLLAIfBv5lnxf+PuDpwOXAo8C747XULCPofqHv+89Bh/H/ued5L0rMuBgZcXuXgGuB/9r3cW63t+/7v+X7/lPRmn91yCIze4yPoB2YvWN8RN0zd4yPsb2NHOPWITGM53k3oJOD/pnv+8NOQg8BT+17fx4QJQo95nne2WE7Z6MTiXLDCNrxPK+IPlH9ie/7/y363Pf9x3zf7/i+3wXejw5754JRdEfJYL7vfx/dxxrpy+02H0V3yMuAO3zffyz6IM/bu48PAq8a8vnMHuN9bKV9Jo/xPrbUPYvHeB9b6g4xcoxbh8QgYWb9bwDX+r6/vsViXwOe4XnehaFX+XPAx8LvPgbcEL6+AfiLOO01ySjaPc9TwH8GvuX7/u8PfHd239tXAt+My1aTjKh7wfO8peg18FJ6+nK5zUfc1yN+noFQbo639zP63l4L3DtksVk9xnfUPqPH+Ci6Z/EYH2VfjzByjNvCaBPied6fAi8G9gKPAb+NHmlQBo6Gi33Z9/1f8TzvHPTQv5eHv3058B/RQwL/i+/7/yH8fA/gAz8A/D3ws77vb0o0ygKTavc878eA/wF8Az2UDMKhYJ7n3YYO7QXooXE3Rv2uWWEK3U+jl3nuAh/M0zafcl+fR+dTPM33/RN9beZ1e78cuBi9/z4I/Irv+w8LOcZ31D6jx/goumfxGB91Xzd2jFuHxGKxWCwWS+rYLhuLxWKxWCypYx0Si8VisVgsqWMdEovFYrFYLKljHRKLxWKxWCypYx0Si8VisVgsqWMdEovFYrFYLKljHRKLxWKxWCypYx0Si8VisVgsqfP/ATT1O445HkOrAAAAAElFTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"gdf = gdf[gdf['tile'] == 'T33UUT']\n",
"gdf.set_index('beginposition')['cloudcoverpercentage'].plot(kind='bar')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The products at August 22, September 6 and, because the exact duration of the fire is unknown, the September 19 are selected for download."
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"timestamps = [\n",
" '2018-08-22 10:10:19.024000',\n",
" '2018-09-06 10:10:21.024000',\n",
" '2018-09-19 10:20:21.024000',\n",
"]\n",
"product_uuids_treuenbrietzen = gdf[gdf['beginposition'].isin(timestamps)]['uuid'].values\n",
"\n",
"out_dir.mkdir(exist_ok=True, parents=True)\n",
"downloaded, initialized, failed = api.download_all(product_uuids_treuenbrietzen, out_dir)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(3, 0, 0)"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(downloaded), len(initialized), len(failed)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The files were downloaded to the following paths, which can be copied into the notebook [](02e-spectral-indices) to calculate their NBR."
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[PosixPath('resources/forest_fires/S2A_MSIL2A_20180906T101021_N0208_R022_T33UUT_20180906T131549.zip'),\n",
" PosixPath('resources/forest_fires/S2A_MSIL2A_20180919T102021_N0208_R065_T33UUT_20180919T132226.zip'),\n",
" PosixPath('resources/forest_fires/S2B_MSIL2A_20180822T101019_N0208_R022_T33UUT_20180822T161243.zip')]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"list(out_dir.glob('*2018*T33UUT*.zip'))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Jueterbog 2019\n",
"\n",
"On June 10 2019, _rbb24_, another German public broadcaster, [published an article](https://web.archive.org/web/20190622070809/https://www.rbb24.de/panorama/beitrag/2019/06/waldbrand-jueterbog-sonntag-lage-unter-kontrolle.html) noting that \"one of the largest fires of the past decades\" was successfully extinguished. The fire fighters were active for a week, starting on Monday, June 3. The fire was burning close to the former military site in Jüterbog, Brandenburg.\n",
"\n",
"### Area of Interest Definition"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"gdf['tile'] = gdf['title'].apply(lambda t: t.split('_')[5])\n",
"gdf = gdf[gdf['tile'] == 'T33UUT']\n",
"gdf.set_index('beginposition')['cloudcoverpercentage'].plot(kind='bar')"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"timestamps = [\n",
" '2019-06-03 10:10:31.024000', # closely before the fire?\n",
" '2019-06-13 10:10:31.024000', # after the fire, low clouds\n",
" '2019-06-26 10:20:31.024000', # another post-fire product within cloud coverage < 30%,\n",
"]\n",
"product_uuids_jueterbog = gdf[gdf['beginposition'].isin(timestamps)]['uuid'].values\n",
"downloaded, initialized, failed = api.download_all(product_uuids_jueterbog, out_dir)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(3, 0, 0)"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"len(downloaded), len(initialized), len(failed)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The downloaded products can be found at the following paths, which, like previously, can be used to calculate the NBR in the product NBR using the code in [](02e-spectral-indices)."
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[PosixPath('resources/forest_fires/S2A_MSIL2A_20190613T101031_N0212_R022_T33UUT_20190614T125329.zip'),\n",
" PosixPath('resources/forest_fires/S2A_MSIL2A_20190603T101031_N0212_R022_T33UUT_20190603T114652.zip'),\n",
" PosixPath('resources/forest_fires/S2A_MSIL2A_20190626T102031_N0212_R065_T33UUT_20190626T125319.zip')]"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"list(out_dir.glob('*2019*T33UUT*.zip'))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Lübtheen 2019\n",
"\n",
"The wildfire in Lübtheen 2019 is interesting because it was burning during a Sentinel-2 flyover, which ESA took as an opportunity to [publish an animation](https://web.archive.org/web/20210223154159/https://earth.esa.int/eogateway/news/german-wildfire) that shows the actively burning fire both as a true-color image and as a faux-color composite that highlights the burn scar using short-wave infrared measurements. The article the animation was published in mentions that the fire started on June 30 2019. The flyover happened on July 1.\n",
"\n",
"The animation provides an unusually clear depiction of the active fire area: There is no such information in the other two articles. This information can therefore be related to own NBR calculations and information published by official sources in follow-up reports.\n",
"\n",
"### Area of Interest Definition"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"