Skip to content

Instantly share code, notes, and snippets.

@taoning
Last active January 5, 2023 01:00
Show Gist options
  • Select an option

  • Save taoning/c8658872d95deae690b1bed1adb7557f to your computer and use it in GitHub Desktop.

Select an option

Save taoning/c8658872d95deae690b1bed1adb7557f to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "eb258acd",
"metadata": {},
"outputs": [],
"source": [
"import datetime\n",
"import os\n",
"import tempfile\n",
"\n",
"import matplotlib.pyplot as plt\n",
"from matplotlib.colors import LogNorm\n",
"import numpy as np\n",
"import pyradiance as pr"
]
},
{
"cell_type": "markdown",
"id": "5ad6fc91",
"metadata": {},
"source": [
"## Call oconv to the generate an octree, which returns the actual octree in bytes. Then we write the octree into a file"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "01475085",
"metadata": {},
"outputs": [],
"source": [
"octree = pr.oconv(\"Objects/materials.mat\", \"Objects/ceiling.rad\", \"Objects/floor.rad\", \"Objects/sky.rad\")\n",
"with open(\"test.oct\", \"wb\") as wtr:\n",
" wtr.write(octree)"
]
},
{
"cell_type": "markdown",
"id": "59250113",
"metadata": {},
"source": [
"## Call rtrace with the generated octree"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "50db0455",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"b'3.638473e+02\\t3.638473e+02\\t3.638473e+02\\t\\n'"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"pr.rtrace(b\"1 1 1 0 0 1\", \"test.oct\", header=False, params=[\"-I\", \"-ab\", \"1\"])"
]
},
{
"cell_type": "markdown",
"id": "4e868e1c",
"metadata": {},
"source": [
"## We can also build a scene using the Scene class"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "e22adc11",
"metadata": {},
"outputs": [],
"source": [
"scene = pr.Scene(\"test\")\n",
"scene.add_material(\"Objects/materials.mat\")\n",
"scene.add_surface(\"Objects/ceiling.rad\")\n",
"scene.add_surface(\"Objects/floor.rad\")\n",
"scene.add_source(\"Objects/sky.rad\")\n",
"scene.add_view(pr.load_views(\"Views/v1a.vf\")[0])"
]
},
{
"cell_type": "markdown",
"id": "33dedbd6",
"metadata": {},
"source": [
"## The scene instance can be then used in the render function, which calls rad, rtpict, vwrays, rtrace to do the actual work"
]
},
{
"cell_type": "code",
"execution_count": 25,
"id": "33df8227",
"metadata": {},
"outputs": [],
"source": [
"image = pr.render(scene, quality=\"M\", variability='H', detail=\"L\",\n",
" nproc=5, ambbounce=2, ambcache=False)"
]
},
{
"cell_type": "markdown",
"id": "a18e4c73",
"metadata": {},
"source": [
"## Here we use pvalue to convert the image into a numpy array of dimension (512, 512, 3)"
]
},
{
"cell_type": "code",
"execution_count": 26,
"id": "ce920c71",
"metadata": {},
"outputs": [],
"source": [
"xres, yres = pr.get_image_dimensions(image)\n",
"pixels = pr.pvalue(image, header=False, outform=\"f\", resstr=False)\n",
"imgarray = np.frombuffer(pixels, dtype=np.single).reshape(xres, yres, 3)"
]
},
{
"cell_type": "markdown",
"id": "0a83bc40",
"metadata": {},
"source": [
"## We can then use this numpy to do many other things, for example display as falsecolor"
]
},
{
"cell_type": "markdown",
"id": "3e95a9cf",
"metadata": {},
"source": [
"## we need to first get the luminance value from RGB"
]
},
{
"cell_type": "code",
"execution_count": 30,
"id": "a0cfbc05",
"metadata": {},
"outputs": [],
"source": [
"lum = imgarray[:,:,0]*47.4 + imgarray[:,:,1]*119.9 + imgarray[:,:,2]*11.6"
]
},
{
"cell_type": "markdown",
"id": "f8710ef0",
"metadata": {},
"source": [
"## We can then display the luminance using matplotlib"
]
},
{
"cell_type": "code",
"execution_count": 39,
"id": "387cb049",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x12a7bd840>"
]
},
"execution_count": 39,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUkAAAD+CAYAAAC3MvfRAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACDOElEQVR4nO29abBl2VUe+K1z75sys3KuIWuQSjWoNJSEJIQQxoMG00hMUtDMbpBBptoO6MBhu21oOxo7wh2Bo6PBpk3YXWERCIeNjAcagTGYZmja3WYQSEgqFUg1qrJKNVcOlZXDe++s/rH32nuttfc599z37st338u7ql7ec/bZwzrD/s631tp7H2JmLGQhC1nIQurS7LYCC1nIQhYyz7IAyYUsZCEL6ZEFSC5kIQtZSI8sQHIhC1nIQnpkAZILWchCFtIjC5BcyEIWspAeWYDkQhaykIX0yAIkF7KQhSykR3YEJInofUT0p0T0EBH98E60sZCFLOTaEiI6SESfIKJvuJrtzhwkiWgE4KcAvB/AGwB8JxG9YdbtLGQhC9nbQkQ/TUTPEtFnXXoXyfo7AH7+6mq5M0zyHQAeYuZHmPkKgI8B+MAOtLOQhSxkb8vPAHifTugiWUT0NQA+B+DZq63keAfqvAXAE2r/NICv7FXi8AFeuuHo1lojgHozDJub3llHT+X5EMe9bbTVfxIxy9bm2U+supqB+w/Psn3XXlcdcpWnrLizvq2L12SH7ssMhEFYbxswCC9//pnnmfn6rdb1te8+yC+8uDko7x9++vIDAC6ppPuZ+X6jG/PvENHtrmgiWQBAREKyDgE4iACcF4noV5i53dKJTCk7AZKDhIjuA3AfAIyvP4Lb/9f7YrrOxWbfgJICRyKAGSAqH1YpT7EudvXk46E8M9XTVTmK7VudbPtUScvptv6+8yzTOZ1XDTRzfg66m4Ns9XX1FzrG1kKecAWkzdp5l0AW2zPnUAff2r0rdbLXvnaMlY65bjnWAbaVthiTrq/OX8nX2U5MVfdikj5DRM5d67PejvDcpUNYb0cAgN9+748/vp02nn9xE7/3a7cOyrt06uFLzPz2LTRTJVnM/IMAQER/GcDzVwsggZ0ByScB3Kb2b41pRuJb5X4AWLvrZs7pcACiyqS0gHYsHSv9UGcnZAQkLUCH804ASNs+xzwaYDm278GzpiszlUBZKaDPW6lUqdWdE3wnZXfesd4SW/QZVtthl1MDZz5gQZN8GX8tffNR/9p1ctmydjGv1CmAVgM2f01r9XnxWpT3K78wuqA9PxfqbqZrxf7IliVfAzL1bLQjPH/pEDbaZssWSCmMzauHTXUNmH/mare5Ez7JPwBwNxG9hoiWAXwHgI/3FfC3kDn+gZBWcktp+a8ozxnoUjrrfJTqkEN6W7cPc4xMvZ3nMSBPPR+ZdkNK+WBzR3eyna2GwF1Mrd55hIP6lwwrrRgo2JCvxV9TrZ9+QfUBpBZhkQUQqj8DhmSvaW+9bl+n6Wcmt2ifj5oUTbv7XjD9KUSXJfdEb7QjPH/5IDa4CSya6i/maYUBtOBBf9uQQSTrasrMQZKZNwD8IIBfA/AggJ9n5gcml+t6ZAJAScc1D70HRFWmJulhYp2PEphm0CwhqnrbBz4LrMC7BsY1nfX51jq6BVSlDped1r8MPNPs1LsjPd0p9hy2UtIgeOZ9RYWVXQ8E5rZxTi+aUY1nV0E3IGVGbvflxVDS4HREP0FGf0J8MQlAEaeXStfLeahYBstGh8AgDyYTe9bSDvxvGzI1ydpp2RGfJDP/CoBfmaqQsX2lnniIXMZISwIoKBObcjmKtl721OSuUkKg7QWMaBl509wJg0DMYFIdTPQhDQe5dDbhw1Y2EbuM9m4WwAw0FZOfY4IxgSOgFS4CzmZgte1URyUPZYDS19jfrlyjXHvv4621bOtLLy9xmZBvy95Duf9ybJjJWdaogS2fhT6XEkNTLfISd3UULokpxN4P/eIlbLQNXrh8IDDI7bG5qjAY68PN7SNEdD+AX2LmX6plIKKfA/AuACeJ6DSAH2XmjxCRkKwRgJ8eQrJ2UnYtcFMIA0wCKdYBr/11Mat7uORhUYAYkYukPGAqqdZpKlYdUh3XgKiadkARwLPWAww4UM6vqirLONDzZY1pGI95tlgDIoIGan1Klqb6DqfBv8aiuqQW9AHKe9FVOjND3ZZ9BQIw7BEd51jTqzCblTshtzvQylB1WH23K56Nhru43jZ48fIBbLSjGbVTaxnYHA6+Z5n5vt76mL+zI316krWDMl/TEuOdZxaaaEEtCxmbJR1y5l/X05zMKWPKRgPK12mEymMVP6SYWTYvF77OWv36PHWnrQFLp4YVk7uWr8sPqOPiyTcs18dmNFBEuusO7ancrYfU27ddA17m8tr1qaPPLZ2pesEC+np3G8oE1xZxukZcLTFcpN7afdhoG5y5cgCb3ORnb9JJb1Gugk9y7mQ+mCTTIDZhI8DxkYkUwZt8SPvWa+XN6Eq/V3WSoqEqu6pDTGcCwGJWpSi5fli80YZk6nYznPLcaxFtw7rRcx2p0n7eKzSlvrpifZVXBPQLoSb+FtSG7WitEgBT3rdmeIUpJvbdzyIn6ihtmnvZ9XLRZjBDuzK2i1f5WlhmusENzlxZi1FsrYPKNCNhAJvX4Dex5ohJlkEKjmxSd0XPKENnd1UVdYS8NktGueItT/qHEgiG+ipsEugFu5SlMBVLsxgd++nlwAr2p3hea1m5+LcEFR9x97DaZ2x3ASy5na5xhl3XUwNBF2kybpKOegb57bjUJ7O6Hn0JaZiSzjupxb4Xptd3kwlnr6zGKDbnv556tivtwL/9JPPBJJOUzIYrT3liVQwQRQ9menEr8BPgUvSQIyoG5hdhMCJB8DdSAqKI0cnsIt1LfMCDMzGTB9oPbdGBKN3R9bi/8LJg4wMTZiZDOWrMq+bD1ZdOX1nT4QsmTSafPmIYOpUA39fBy5eK/HSzMslG0WzVdcm4QFb/mvJDLJMKKy+uG9m8+kzMta0wzdr1nwReNduD7JUHgbHBDc5eWcMmNyl4pwNFuq1ZcT8GT+OT3DcyN0zSPGIVSuXfj9XhMxVWKPtDxzCW2rg2Yl3sdND+zLRfqcMGi7wPsmgqAUKnphVVO8FLM2Lk61I/29jRSDT1QGvHQaqrk/a0D7ZohStpqq78l9Exg4ZtT59vmghgrIFuseMfrM55K/9NfIrUc7YdkLLXV4/cZLQgvLyxgpZthJ2I0VDptZwVs2QG1gf+IUa3iegbZ9D0rsr8MMkIcGnoDutOmo+DuvxXHQOSFaOQt23K5+mPsFH1iOqIdhllVvlsNQqYUAwRAtgwlCL4wLHruo7u1a2xQG8ukytXY4jdYrmMMbcjk64zSDYDubs6aJcvspKxeHEI3y2vB6Oh/rMKOjGqwah03G4VIwi8egm9yxq3AlCe5crrreUG5zdWsNk2sa04DE2XIV2D3IdZwCRhc3g9E6Pbe0XmiElmMy9I5BHsHni2rFD7MTkGgBKjZKlDoaQCJ2GDXaxPfoxf0ytjfsgfKX5N3bJHZR4/PKlSKp3rEJkEhb2lPOsju1GAQgHeZUtEdlaLz5+vG5kXUx8ryj7s6SChHPJUaqR1MroIKNL2GVsGxBpHZ7QALm6OA4MkRhPzBl8kFIt0Txth0IymScIAWh72t59krpikEJfMKMvDed+xkNrrHYoFxmRiYZ02L8cnXd7KJlBCtoNkVcvxkOUUu9x6V+fV7DLVzf4aeF/lBJpWa6eSXSL4NqikWFOhR24a6toYds3RPC9MXs3G+lmkaTLdnzqDNFmrOlXy1donr6sfg0rqXygzY1aoYG0Uqw3hyuYYGzxKx2RiQJrUAKSLHnziuaZZyRRMct/I/IAkLMBk0zZTOTOQOyFZBiCoDm1AscIOJOgjLIFiJmu2K6SU5kSfnmdFzqMOSl0gGI5laq+hvSyXTGvBSgNyeUt73LL55YB84nNfNzW75oP7F4wBFn1qlXbJbNdfBF3qkhrfNemU9ABycoGOGkDaMh7Eu+/3NNLFqxmEy5tjbHCDBmGGl1hBwl51zCZfZwXqM0BKxgIk50LMGz4+dRLJ1NMRM2Bkf1kQOy4yVxw5iPEZZU4YgBLpKTN+RKejxFVlL+T3HU01PWEBBylJbPWrsckaG6uatzUGCMCO3UsXuGDvum7Nn8J7IlM7q1HW2+qTc5kRAjX1IB2f65WrPN01DBcdOc/PVkfDtXbd1MztAKW6SilNALKVoBtl0BPrpxXmDlIrY2X/dGE5bVEYwDrPjYfuqsncnDGrjRD1jaBmnlVSs3FUOb3PSGZGOEQpX8ir6q6aInnmTXFMYbceL+l1QGxHJyd/qTkXvRc6bLniTNG8KVPmoc48dd4XMvlO1GVmhg6a3yC2zgyE6c9f5UQOa1c+3zXNOkm1XTuL3FauZ4jIfaZUuZ2n03u9By7rNp0kT3n682Mgm3SNcyS7IVaMPfoqUxmgQS63Xe020Qz6wyK6vQOifTvpJz7EqtPk6DfCW1XnhX+0bU9OvipSzBEAk4oKszIjEy3oZhZJ/cI8Y1XeaNPJNqqr/UB3/iokQfsUa3l8W5qxhL5FnZ3el8r1dbfTO0ZxEN3y15vcdfBi2Vx/zWy2rd80K+dntpgWIpOWRZqzb9CewbTin6E2Pf/AiBhtZI8NkIZuUYxwM4sZHrVNLyPA8YptSTt8KN2+iW7PD0ii4l9TAEgqjxkWJAzCoKjNJ36nzFaz6Q3IM68G/qTFMaxOCVilXjecB0x5aiKUfxXZZLWBntrA8XiesqcArBZs6HXLDwbpGmeKbTrAo0oemw5TTt+7BHId/UwDcMls+6T+EvNjIEndiXK4VPd1TNaAe6mwfqFO1LFfeyRdSW3nVZ7Cik9I29kvmYeLJbAE0kVmZYpvVxgLn+SuSWJEca9gIp0UIbCo8EaNWdluIx7zYCctpmorbVhQUh2RkR3lqiuyypWAW5HZ+uPVzdC6IsDFC6B2zVJmDbIdLKw472Gmq9ZP+BeQgaMKPprCd4i/dnpH61O7516/WqroSCaNVAn/Iuqoufry2JokHzCypZPAkBlMMmQtvIjz8Tje04AlsttGgeX2hbB5Dfok5wIkAah+UBg3YTsiYXgA5EHP+S04qoglSSfSPUoqzSaLLucZYxbPTSS1dh6ZOdYGuhvoHbgy9xAxwKosr/RCKDo22WviJb0MNHh7s9Vq0DsMh5yOZVMdSvS8ZirT8OrwKGfMrkYNkF6feG3kZcNyDlsDRs15a9dwRG0y41uK44Tj1FtmoFUVCFNPpjbHfeSXgIDlrAI37fyEMa6azA1IcvwndGQJarjOm264B9LS3BLTSADVmH9p7ARnoJRybjhHrTP4sYUy+0SDTh6fGQEhda76x7Ryx6kYibotxxynMbdtSzZTDfwnVVmCEhvWrvXV0js+snfZtNowm26wJZfDDynKrLFujKYjjv1OAuM+0S/HsJ9HESzTRgJIPf+hVY01kJlCpXUkeWUssAXL7aMkM+EK78yK5/MscwOSIuZ5VGBphqcI4LEqI3nVohZibvhhMobhAdC+RKNLZKAsO1CMVLWhI6tR8RIo3UIXlmlNMu+sSWl8tBUQKv1t6jqpAXUC1ZQOWqaT2gWSz9YDkN/WevpzG7RuZOfbybI9/SEwr4PSRp1hV9tlySK3fjbUfdsO7OQ6Qt0HR1cg97OV9hnBb84U2KSAH3FgljGvvghNyJ4bUWA5C2m3ddZ7U+aDOzPCnS2Gzags7h7n797A9rCON2bVJWPQIJfT39vxC1CYemI2SZe3efVxtL3OHbTnXhsupE81jCGstNLZD7q5W8mSHQCSNJ+Zu5xzX7OapeV8GmZrqnOBUBlI/JcJrbe4XmcxAEm90OwZ6PtW4HO8H+ysjmlFA2uaNQPg0PgyxrSJETFGaNHEoT4NtWjAIV2G/MiwHyANCZI0Qgj05PIqzxZ11sLAYgjQbkp6+NTCCfmAAJaLeirWIWwy5FPmLwuTg/Fl5sdGo2AZkQ71ITFS2Q/MVTvwKdOu0IrpnqyOi/nTPysnnFiVDToFu3x8+vT6F8PwTDvo3xj9LdxpNKuxYc3h9Dssz4pRDNPrbBGqI19uexJwmfuQmK69/wUzJjasu3aOQ6TGYUWfEbU4OLqM5WYjsMLIGMPntigyyFC+5RDdJubkq2yA7ItUL+vwkMUXS6pzFkxyqsDNYgjQzCXe6ASUQALLHEQJwKSIVQK7NFxH1WFnwdQ6ChzYIOdP4JyBUQOl5JT1KUO7NtINrU/UNQ8B8g/thId4AijmgUcV3lghrhkcypdCpcmCWUliDZw0i0xA3AF8E9VU1zSlm2mjdZakSxk/MPWcV6pRE/WtgYsf85A1kQyM68aXsDa6EoI0IBAILRoQM4jCx7xaNAkgJXDTQBa2jW6jamAn9yEh6NsVxiJws8sSnoL0+BMsq5QbbZgj0jOc/H4ebFXkOtSr2SAqdcfuJA8X2Y5YDrlRQArEOeCWnXp/Z4KPyKz8vmawechSCQgmil8BzZymjhkamoHEL9SRIsBJBwfMHfRPn49l6bUyVvqCPFI8WAQyza4fwASYysE+XLma0lZ+AcBbLsNOo9BCv6QZwSS+fvk8VpoNtNwkYAuPWhuYH7cgIlD8bbnJ/sl4HdoIkPrtLhDWRhYJuDHC25TNGQSA9prMD0i6py+blQrokB85GKCKKZyBUmUzK8j4RTRSIKZjLGFtOEsZGNGAIA8nQ4YsDXmsnMpFvfU67DXQgZ0kwsCR2baPiutB1rqtkgFKdVlDw76LM3LXs4Iwdre7K2uY6xvOYnmnsMdKUM4BMqU0xVTT9exmnpPFnnRDjJNLL+PQ6HJkj4gckmJQhBJbs1ejVcwYyWoSU7sVVonwUm4kn+SdAZdkENZ5fiDjasncnHHRfzhbDSYim1hfjSdklpHNck8XO8zVBB6SkBlgmqtsfJpQKwxpYKaUR85LjqWOqO05hzK9YwzVuE7NnOWcasxHvW1sPQY4/DjNUH9jhr14ACv3e8fiOWJZyjCuw2z1KmupsdHu9rRa+QVQe2FNK/J8ZifMiFrcsHwex8avoAXFrxu20fpoEpvM7DGky1TZlhsFjJFVpsfbscp4T9tY5yy+O8OABGWuKZkbkAQiwAFFp/asMmEMAbKKt4BQqEdMF8uGNF4IyGizOc1KoNyenrpo7HsjapiQwtNibKOonEpl1kcutTj3tJdzdeFKbcaS2ddEU1wUJs+EWTy2CnPdQ525Dn0OtWEo6Vil/gLs4kkPA654dSewQfeqtOczqJ1aq6G0LjuiFjevnMHxpQtx/nMT88dtZVbLM0QIgNkSxWe/TYyzZUJDeQyl+CyFVQrVlGmNo2HvoF5h0MLc3lUR8yFuU6CQmSklMpmBMnf+UNAuMyaMTvnIOHIkkiKqQQDFR8Viw0XndoCZOxNBxkfKdi1AhO6kpLuf+SPdV1hUd9lKoo+iR4oubIRyYof48yegAnhy2v5FUBv0LnmoqEA2rc4awCexTj1A21fdO/jeAeoQOCi1Ke9MQ4zbVl7E9csvBwYZrYlNCs8zUZgTTVpLpnjSlrmxqj6NE1YAqVllfhbq92srsgjc7KpQvusk5mxCKgOiOjjiI9WaOSZxz23fLDzpnp5tZHaoB5e7KYCS17FTP7yGzIboFM6p0MsDq2FrkiFzFz8S0V+vXG/NT2fBq94GKuDDlfJsi/a+IDyAdnVoy/Rs/tCQNrd9s7V6tQ5b4Uh1/p9TGmLcvvo8blw+F5lYAyLCZspF+aJyE9hiDOYE10Ib0xGJQH6+xITODBJVVplWLd+mMGMxd3v3RBtC8npUhEXbzpD0RC2rkWYRccZzQkZr+urIbQYgzQhDW9k/CTt7Bt40lc4qZyTDligO7bDl0hVQM0hyO0DJ2tQLQs5NXUFzOYtrrMGjwshMRV3dvqtBNntFwU5hyOypnC0HmqxyDsTjtp4HpNXvA8ji7LaAkDUtdVsjYty59ixuXn4JLZoAMIQEhqF4GNBDcWAPGT8l4vds2sg4KUbDGY3yY3LMl1YtV4tgyPGuSRbTSAjcDJ6WeISI7gfwS8z8S9tufBdlTkBSk8j4+CpmJuZhGvahGaOY3wrgzAIYJi9yZxB/o6YCkg8RgGLPFXZo52MLbykxw46HtOlighfLqNXAUEDNm/xlgwVTlpOhdC0diKkyfrmvrJnXP71nKkc0w1SMz13feo2cgmZhbGD34HhzC92xxGrJlpgID1s2ReVp0O0HGRHjnrUv4bbVF9Fygw0Aev6PgCAxYzOCnYyGQALAFkQNNiM7DHlaSMAG3AIktgPltR7VG6AhOxd8uzJF4GYxmHymwigALB+IQtL5I2CQKpKG3CCZxAWjVA8lm4cxA5Fu048NlKg0E+fFYFIz8tByYjC5+2fualoofKb2cmjgNMwynQ/yS8Md90CW3QRc0c2de6WNXj+eucZxg+xLpIZSmQnaVYP0sZrU0i0jT2di8utr79O2gh+5Pg3o4YxG1OKNB57C7avPY5MJm8RABENmwghttFaQwDA8m5FhpuHiciwHdOQEWN3E+AiCYqkwqFxlqrlxtiAyK+hak/kASRFF7IDIKhV7NLYUK75DioFKomeUKR3JBPFAZFSRjuunKjpGaViTbgN52mJ2G2SWaWcXaYCvXhbXycuXQNXvCduJfdUFMJM/ml8Q8Yzcr6vPMbhcUxcrzFdRWCS7/PpFoKvX16Nv2I/n5jYTYxqK5a96jWs3xLhz9Vncs/alaESPimeDwI6RZaBENLflN0W/YX2TMvMm6SGnwoyWGjRpJk54rrbKlb0shgDtpigfY0oCUvfUQBkIWx4zmMGPYjkBx7hf+AKRWJ80lEALZafzw2MsGKtTgO+IGQTsAQ1Ammn6dA3S5TFdV+fMGH0e0MwYqWX5DIbV3YPOsG4mdyDp28Ei07a77r5dYz047XyQqs60S8A0bHIKU9vfBXV3AQBjavG2g4/hjtXnQMTRB5ke2NBeJneGGRIaENoQ0OEGrAI22Z/eQkxxQvy+dTreIE2djcOGWOWZhYT6FiC5e6KBzgySDpJMaEkX0JPhMulhg1SkBp4by0MxQZWfkYbEQJWpTT2ENoMF4jjraHXJoA6tpxnQrQG5HNhdPya8y7VDfluzXQcIqZP52SeOgWupvsR8gvKFdnZSzj5jBWYl6E3DgfK1aSpujPyeq8HxZClPJ48mGFGLrzz0CF6/9hQ20cQPeEkBGSjO2Ihqbso7JD6TG5FFBnMbilFSel6IIkBxixZNGA7GlKYxtkRoYmQ8fAcH6f7OJr5Ni8837KaYh08HMRR7DIwiJpoCGUCNGasAJdTFKqviFx7QFLjosWY5CKI6TKwrzxuP+URvlTeZ1qosquyvfn0aD1AOBXpnvGjdYrHye9OSNY8I0D5Jn6+LEeo33ESmZtheySL9wPSwJyeuXgCUa/H6aXD0Mt2q8NnxoPeJGH/m0EN484HTcQwkh9EKmlCTu1+OSWbglLUDkFhjfm7bsA8YoDT+SQnsAMn8ZgqLZGxXGJgmur1vZCJ3JqKfJqJnieizKu04Ef06EX0h/h6L6UREP0lEDxHRp4nobVNpYwgUKboI339MOqcngjTZzFlUArPaT0+jBjGdL9SpTfU0gLfwZVHZFlQHrT6kquuz3tY/uZ2Wc/vyuVDRS59bzU1fNF9QwO78zqj1UGa2uZLaJfYaWnYXRgB0tSg5uxivBlSVwiVbHRrZJrOV78yIWrzj4GN484EnsUSbGKPFmNr4m/cbtBhRi5H86m20ce3INq0dKXkak4eLPGHdSE6/pH7DNkLdMxhMzhyHIA34208y5Gx+BsD7XNoPA/gNZr4bwG/EfQB4P4C74999AP7ZYE1YgYwHy7TtANFzF67k4wwguZyqUtpTjNEDnQFWDcQK2GwZMttmYQJYYCi3wqg5NoBdoYgc8hi9AMP86vrUWi/rJ3SAFDvgcS+RXlNd6hUGVmGoud18U+rvx3y+tv56JLdkyxXA7MhL5gHJgaYxtfiL130e773u87iuuYIxOAFkBsFNjLGJsQLFMTarQCmgWCy2axbiZZNHA6MHzAyWmAlIAmEw+ZC//SQTz4aZfwfAiy75AwA+Grc/CuCDKv1nOcjvAjhKRKcGa8MagHx6ZpUFmGrm5dmgSxP2qAGGK4BTZYqpHVJ6aKBUjE7ySD6tZwowOSPTsWB/GUSPnF4HWS6gQh8lk0Q+XZeq9KvyE6pqjKqvqQAwrV+1SdNu+SqUdE7DwKRG+2v1S68n8sdK6XYp2CNjavG+w5/HVx96HMvEGEVQ00CZAZMLMBzBAmPBLCFA2SagtIySO4GSyKe1ZgHlrQojvMCH/O0n2apP8kZm/lLcfhrAjXH7FgBPqHynY9qX4ISI7kNgmxidOGrtLM4gYDqcD+oQ8tRFGRKU0jJySppdhDceZs1GZHZMrpJEC51WoSU6Op6TXZRctW2GGEV6VptW2CcpUp9hxx5n7yvM7E10yTqV4FdvM4hezMsuIhEj9ZXrpM9fR56L4TtDWI+urOOKZX0G1gn/Oil5t7R079pz+POHnkALYCNdzHgfGQC1xlcov2n2C4WJiczAKD4HoYzMlmGAWzAaNIi+yOiTZPFNMmJ/YIDy4rwEHfm2X1jcniw+KbslYWamLXxliJnvB3A/AKzcfmtJifRDFdOLaXqsdzOAyoNoZuhIGoC0yAXs05Mi1GkcoxoEzkiBjzzUQthMBPYYbdSrC+VxkXnoUaqbuQSSFDl3g6zNtg42uG5dqVOfX5qFg3wtBUw0P62aomR+zPVT/N1s16Tm1bTgWQ9kmbTipdTFuZEi3UNFuwF8uTG1+LrDX8RXHTyNZSJsqJe0nlAABlqKQJdAEPk+8gijFMlGjFA34ZOyqQ9kYMy/Ur9Op5QuAZoQ+VYL9c6A3DEwzWDya35a4jNEdIqZvxTN6Wdj+pMAblP5bo1pk0WDY9qXXplBMTFKVTAPIRFKYzuTrj+xSgWUAYw4FdVDU0zeGlCitnRXhGKFoZLJMrd8wAzniXXbC9MtBfjBAkoN9DQAECYzLXK//pz97QtpXie74VlkV6BmqFQIfmcEvyxrrkYur+7vmFp8y9HH8O5DT4GxiY00jlGAUABM5lBrYGsNaxzF5SkYIS0P9QnDgUZJAQuUAYhlULlmlHpcZZwNriLfNdfJtDLl3O19My1xq9z54wA+FLc/BOAXVfr3xCj3OxEuVGFqV0VFp6vHODOqIk9MSwEWzg85a3+mK2Oj3qVvUeow3da1bfyZLG36pqI+Efh9BFrYYa67hIkQWa9xNCrRB7Z+g0ru/ZFAr3bdVX659p4DZvC0Jr98wa+mGVSZnF+O1RlkfUYNpz8N4ulPjYfswwgyJfWXE22pt6y9iPce+hLGBIzQYEQNRmgwpgZjEEZEGCOAXvBNBh9lE32JY7TBd5l8kWz8leMYkNH+SvFB2t8W9WCN9VmK71b8lLMQ8ZZO+ttPMpFJEtHPAXgXgJNEdBrAjwL4MQA/T0QfBvA4gG+L2X8FwNcBeAjAKwC+d7AmCiDSjqE74XXIKr1mfie2psckkjInWRiGAK+aEaPq0gOstdtTT1U0/klFGQX0QhP1+dLCVDUls749mJ0C5rWpXLDUCCyukJtkZIGli2rpzkUaSiq6phtTVlNjs+WclXz9auokZko2v2e5cOc+iUVC3aOQP7PPwPQYf+HQ0/imw19MUeKGCCM1Iyat3xifC47sjZONwMHXmNggICuJs2KJTXyBhry5TPBFZqYqS6rJghfCKFsCGsUwE2eYBZPkxTduqsLM39lx6L2VvAzgB7asTUIS3ZWheiXlbVLmbEKZTHV8UCdUrUAQUoTSohUgYZFs522zakMeeecXZKV7zXTP9nrYDtkNqgOiS6xflNJzvfNsnv6nnpWN1bXwhb4O6VK74wJD5edla+3roT25Z2rgtiqoIf3s9KxqocHP1yVZp0MD96QV9S1Ri28+8jj+3KFncKDZMMcatVRZAMQ2aslhoQkwWgqrgnM0t1NeAsbxpso9HQlwQpZRi4pFkz1PVWzB3KQ+kIETDihjpJnjKkIzkMUCF7sp/tlOT2+FWaokhSUZXAXYLG3JQCkJqg1ZNq2chij0ToAzHJfvzIT0nEcULcAugqPxFTrA62JQ+hLlgI/2pToeJLrCtafq7wIi054HZMNqbX4NpJNmsdRYZa5Fa9QFXXV9u4Cz67r6NA+YS9Tivz36OP7ioaeq4wwJQIMGIzAYTWR9lPyTTOGLhgKELTFGiR0G0Gs4+ykZnH2VCjwzuDbIwUmO7LP0SRorJdY1mxk3tO8Gig+R+QDJOilBemwN0KhjiZy5RSTEFE3YWmGQ8elLQ3dcGcQypPppin57MzeCqz4NHfDRjFQ+FC9Aw6mwBIcywCbwho/2qqi7nLLSoQYItamPNm9XLwr6JjJclJPt/CLrjUyTAKTNU//+TbpKjkH2B4S8jnWALAFYpzCA162exXsOfal3IDYR0HCYHjjiJgVsAoPkBJgjkk/AtpExQrFACtwzASIlMAWHsZKeWXIR9QagPiomHwBrgDQ0aLvCANYXILmL4k1sLZolSh7NKqWUkEPJ74cEqfozq8x+zuI5MvWrSLa8phUAJgA150D17UQsO/xxnE831VhhZ+mcHfsdItJ+n/jZO10mvum/DkEptiPWYB0gc15djmt5UeePfnzpEEyQK17L+7rVs/jLx7+AsQzT6RHxTwqTbAXkom5t/B3Js4bgqWzjbxgoHu7xSA0bGqWXcJMAU1Y05+rwoLpPUsB3+7JgkvMh3mQW0fhD7ld3KW02OwaqgdGAJyGwyng0t+2G/2igVOauXhCi0Fe0i1Qs+RdJ1a/Mo14zOebzDMlfEs80qwPgCzX7TWTdlnQTD251Bgt1cu58Uh0lUCWANC8vju3Ib07vG7hengebPDrvmFrcu3oG3338IRwdXanUUJcGTYhWAxH8rNmdfwEmxkiWM+MQ4WYO0WkxrRsZVM7i52zjtW8itGa4JSCmxntDAEWgbDmyyRkwSQD7bjbNEJkfkPTWjwZLf7wa2NFdVZmHHigrgBFYTiwn2RMA5KhlyFEBsMQMJVX5BVXEPHfIoIgZsB4rSIPBlSlmv2njwdMy0AI5a9dZX+Oqr60+IN1WWTI8OOYpdQUA7mJ3FrBsuXygPpwI6b6U59Cd5tmozvuuQ0/j244+GhjdFEKExCYzq0T8Hg2Sud0Q0DBHttkm4GxACv6UbxISwRZ/pI14ezYpvsdg7qvP1M5AFtHt3RZWv+ZpVr1eP921Y6l8BCEooFT5WLFHqcuCF0XAcRF0Md8FLmPdTDX/oAY8ynUJKxVWCSCtQyk6KOAAMvjZoE+N+UlZDaQ2Mi1pw4NFftiN1swyM3XBjU5qBKoFwo6B3n78pPDifp/qMPZo07KMqcW7Dz2Nrz18emqAzLoGs1s+1NVERtgEqMMoatFGEOTIC2WOdg7UUIyGy4DwON5RgSqr2TfsZ99Am96BuTZbO6VCFub2rglZk8zbQQr4JHvBiHqAEtBgqcEh0hsVgUwgpfL7cZKhfCqQ54VHMK6voo1selZYqTHlZV+AOvNNkyfVPzCa7Ldhtm3bEkjqZHAV8TqU4JSHPklUvs48pT7RrPYJipLtTgLIPpJ97+oZfOsWGKRvi0AYRU7YYhMjUIpyC4sU/2NDHPfF3NZmd5OCPPLwcPosA8eAX8jXqGc4PmRgZDNbTPvtSohuL5jkrkl4NuWRdmCYc8Vf9yB3Mk5KlWcWpxih4KEC3fRMJvYWgbJgh8jA7UEbmk1G8HCmbSomKqYyuirNCMuOnXV06RXw8ellfQr0i5IdEkE8nUBFSr9lOdYR6phWp0hLxyYvAlJjj67qJPesnsV3H39oWwCZ6iZCwxR/m7yaD4RFBvDK7DLAaQOK7DL/Mjg+uxzNeFa1INSkxk766Yra9J4Fk2QgTMW8xmQ+zrjAPMrp+rlIVlt4ZysrTv2RtThVGpvyOQ+7+vQ+m3xBu6Q2B+DRvyF7KM9KB2GRrOpPzNLhjF5KzT/bzEqvClRoHXxZP4yoS8pyOZWKA/Za1kQDpE4bAsY+/j8MIOUBcGUrcs/qWXz/iT+dKkjT37awyfDVmoaiv5HCl7XDfjmlME9TtL8EWeosrzWZpyS2yVyXqYfmGHGcysgzW09yikV3jxDR/UT0jTNpeBdlbpiktUmRZsuwMEOdL4kuoOvIbDGnC7jk0Xe6QjH+hIkWs2UUa0pWfzKDZEUgG6027FBYn/dPkgBBQjc0RftIAaA+gNAsuCtQos1VNmmsWF89cKPb0GxQY2XZTq4Tplz9uOTR7oWhzLYwPPqQG8BrV8/ivz/xpzgyI4BMehCBOMzlbrkJ6z6y9UvmPzG/tbnto91NYpNIv3nxDJmJIyMx8ne6WwAUPvMwixPjqcztfbPAxXyBpIixnJ0JLsdZ7SRk8OkOKGN6dZgQFPdQbCx9ZAyA9ghSfLB1pxefZpqvTWYvDx4WldQUyATG6VxtWdFLDwEyl07Ga2qTmktQSiCn/rYi1j+LyqKubOq3vK4GXvkFlf2P+p5zvVillolpBNyzchb37QBASnuBR3KMWiOzSuWbbCjP1RY22MbfEOGWoT4tmGTigfgpGwWogIyfbEERVKEIQjuTRSeCJlt9YvauzCdI6h4lfUcxwSSGDDpWCVhGabI4NucGhwNQ86/VXGxkJphBjeP/dlB5HtaDlJbR10XEIxvNs250OTke9E7AKcyW7WwKy/L8JbYfKbPDf7J+mcGVqFRjmSVrzNc4tywHSmaYoufpJ/8rr41UZeW8usCX9IY69NqV2ZrYNWkg4Bf+hBE3yC+nBvJFbR3UyVMWc3RbRijkiDhDGCRipDs/0w04+yHTdZ08KH6ILAI3uygyjS/soMIA4yEBS+mALfJTp21c6LRsyhqgFM6mgVIOmfy53lBdBDmNlcjgovALMOCJwuzOhmVuMquQpy368ZHZx5SvWwCiEvgyu3TGaI8JL3r28Y9uU9iySKN4Tx8z5r4UcC+Aemu5TLUJDZDRxN5JgASQTG4J5DREaBhpumJLYmJbNsmJSVJmlWl4UIhuIw0yD1AsUxpN4Abh07IUX6Kz+lriAiR3WZIf0qJMCZYCVImZ+XyZ4SkKkwE09aYIIsl0dMwwVp6rj/k1UEJHsBVbowzOmY0KmkbWGnXP4CljJzOt0mxT9/7qNEV17jKTJaSwYbN1rMp1kdnKemhduoYcWSabQVSv7Sg6+XJZXzndepQ+l5mu518/voQPH//CjgOkCIEiU6S4zXEfcWZOySZb5YscxeciG+Cspi2KP5KD+Q2GRLjZEQWZx71dYRA22vmI9V5NmR+QNKYBkLuA6gjkd2WYhE2vssrUhmOM7IGDYHyZqvKA3Tq/m5pYGeKTcTGDVzInk4nUPTA8m9YwZcylmzhO0i0IQWoKW1cJynegACMyP2XbPYyxSNYvAqVnX/CoBqiT5LWrZ/F9x7+AE6PLA3LPRgQaiRgjphTgE2bXRDZJiU3m6YUaFPO6knG4kArkNMx5vnb0ajZpaBABaPT7etuy8EnulrD68+aVOKy9CKB6X2UVMG29pa/SMUrl35FKxG9W9n81CDstdqGXUssgogdpg2HHSGYKDWqEgToGEM01ATdtLlu/p/0tAzkZ6MVHNgluMiC6dNX7SG14UK51LT9e0uShfKTeLbO+5FK93Lz0Cr7t6GO4fnypWtNOiYATRRskQlYCQ80m5fEfxeeHidGyDAHSPkjJ3aorLGY2518gbQciMQMGyAtze1clBTkkgfMvEWUToqMvS7BE2XPZ9PZlFYNMaWrIUSBt6vUuwMIZ7MzK5NH0T2PLSZhjPDFheqrHhxk1IUEv+iv6pCh2B0r45OQXZXXKqTpKzMwPq6nxV5ugWLTKQ6gwSGj2KXr1Mdxc16Q2pW4T7R8gr5vxOMhppTS5A7wJgDaUfZTEAqTxBeaCOMIwg58ylIMMByoGlguLbCOwbT9ww1iA5O5KBxFMjBGIrDJmroBHHlsJX0utZpumQFWv52iUQ2ag4rOsDcnx56UZI6AYXY2ZMtLQIYEFra1mrKpIQZg1WFu2l4ciNXKOxJ1sT/FcBa41cKxvJxau2rfn45iofqlUtRluN96y9MquAiQgLFICN4Q2vhAbSGBGsUnDIqE8kTKXW5ZMiyWFbVKYeshp2BESqWghYyZncz4LkJwHkTn6UapMx7NGfTCBQ97PxzLIZTRRQKjSDLuTIFHcNh04VWdnVXvmKNt6/UkDbZ5x5pNNjXCiqiUwpwHkmu5Vr6C7ZC7yL4wzAaJhv/UXggBnaYqX7TMoReZrJjlrnSotDQXJe1bP4luPPrarAAnEexJnWMm/EpCSl15aGUilN/F5CqCa49gcWWTDMpYyR7JHcRxlyw1GACguwtuKOb5NYRA2F4GbXRSZnudAQvdhLTn27MoUx7sMSnZJ2ZzWKgno1hbRYAWaQfWscAZp5Z+UNmtMSfygJk3aBpJvluJQEef/FDCtDSD3olljAr5O9NNAlveLxSwUcOegi2WnXQDXt4iGZq69p5XYOuO1K+dw3y4zSC3ZJykR7gCAMrNKRrFp32TJJpX5DSA82W1c9CKb2xzN7ZabBMDMLWhGM5AXgZtdFFKgkxPkYAUsK9lSumKRha8ylc2sLFVsbFZ3HDlfimQrIATU0CBwMXg88c0iiJPzcczrF8QwxNKq0v3I9h6s2F81cBV2arSomdY1c7oDEB14yqWuY/QAcEzlQ955A0gAjkGG+0uc/ZMypzt/FwcpTYb5hPe3LJXWqOXP1PREACPzFcU4D4qAWQwmZ16Y27sn8fVJLqkTwFhtizkTH7KadA4VKhI0TUSaRM2qHuOTlG0S5prr0lMUAUCmf1eDONbizu1IPQa8MYEtZj3s9RRQLst3DrfRvkydv2JW1xlgrQ0NkPHFYU8v54e7XYVuZJ6FE6PLu+6DrEoaWI4UnEl/MU18k4H9ySybwBclUCMzb2SqYv5EAydQ3ATSUmvCKGdpIPMCJHdHjFWXGKLiKIouFkAqicIaiyh2zKJn6vQBZRdSV8YCGqCUmgS0IbN5xAxWw4gUSKWhO2E+WWZmDBXASTzU6qrOsQggkQUev4CEnBPpfcUCKZYK4MapDg9o5Mo0zmVgF7awAGkHvCe1q68tI6RucdTtQLOBbzn6GI42cwaQ0KCYRk7GaDZDTPA8FEjPtkEExzx2EkAK4gh4NmnIT3huAlDKQylA2fm6mepMFkxyFyVaFkFULwydQVNKVcAlhWSywNpDFvN2ZohFpoLNOtA00x4VOPa+vsUFUBlKpExtUvn7Fq3w5qq8J+qBlFL0sCFhhQwFgBYVoc1ghh5rWbfyNUAS4M7NXhWptzhOZT6Ru1fO4YOHv4jXrpybycySnRAq/ssAKdfQDgMC2OURsztHwgOQZusmD0wH2vjyjb7K+utmalkwyd0UwSOLNwBU/6Bsvnr8SlJLVPQnMVQfRa7MZLHRcAK3nIbuKLUzc5T6dUAlMk7ldHUd2QZ2zCFtYpvTEzCOBr0CuXiZ1ElbqQVUEkC6S5ZYLWUWqY8DbnGNTiCrLFZR1NXDdBQw61wNGPesnMVfOfH5+TOxnSQ2mYBQRbhjWg7cZPOa03MFNDFQw0xqhSAgm95xNCY1II5LsMZRGrNY3oIZ2JzVF8X2kMwVSALIDLJmFadBX7pAiYsJSKKpUmSKYNZrelepkVrdp5Y3nYtB5WoWk12zw8pUxZKdeUCpDY0Rcz/mMACWf0vwkxV7XHAlbvhhO9JW7Xs1ecGKMs3qUrlGlXz+OuwVgAyiY9yRHcZxvwwY87oVNpki3IjDv8QcD+b5KM3AkVemuqeKUbZE2QLapiyi27sthemM1DsSAEbvtjw02qwmKJM9VUKVRNR7XU0XXb5uI1rLPRYJT6uesog04yVNTexTQdWbjqfMniVqU9yeq3SaoqEiH+e0gpXmweT5MnAu56tWNL8GuiafLqsuYjksKC85AgA3ji/i248+ukcAEin4lmGyHCspbpIU6UYeuJMH/lgrKE9VzIEcGZ4e8jW5/W0K49o0t+djZCgj+WIUOcl/lfS8T2W66v/EgTXCpSOaPKYdjvWV+FNuw+cty7HWDWrbAX+tfm9eJmxUqtqgCJKZTa4OqaA620Wl63JFHWTLyWcIQru6XqlzSoC0B1VbcsZiHzCua9bxPccfwm3LF2ql51byEKB8NkSaX+ZVgdIf5fnbRPKphtaky/YILUbUYkQhT7hPbdqexRm0POwPi8837IAoNEjuO50ORTRUv01jK/0MF1cn4BilvNm1nzOVU699v6/MdFI6atPWm896gV09ODwXrF2QPDUxmc6xGHUVAdLA9qaSwQBXTNB6aUAzwEe6nFyXDJA1AExvX+ow571iVfZoMxCAVy1fwPcf/zxuuMqLVcxCBAhlq0n3ldJCuXaIUBx4rgI2+TGW5dLCNqjBJuty4Vdm2xDPZtHdjlF2NVl8vmHmohgWUOKTEWOW5YwcF6QwRQy+OQqngbi2Lfv6U3PCCrumzqnynkUlc7gT4XxFlnFlcONKfssEAQVuwhbJlVflukxpDax6KFBDlm3m6Yb5mL+UvTNuFGBzR/4bly7i6w8/gRuXLlbrmXshhPGS4Mwqkc8zzcRBNFIYQLy22pDK5p9+2ANQhscrRLZbxG9xx2mLs5Br0dyeG5AsRuB0gSWhHtQxmRTLLNqJTNAUr7BJLQVT7clLFUYZy6QFJ8CKsRUKdtddmMWZyQXwVdHSUi0HuMrPGPP7L+p5YNR1eBGTL5wqZV08cyS96U1u4VX5XOXgrUsX8F1HH8VdK+eq7e8FCdfDzroJLzGK72FO4BjuB8yQn/R+V/3D3g3tl8z7FKPe25UQ3Z4PD93VlLkBSUCWDrMuPAOMGiglX2GuZrZIsA9VAk62WXtZpKVCruDEUyqTyA8jsowx5MnlU+CEMjOUerSihDhekbioSwdlEpg5lkmRBWrTWpvSeuC3XqBC0hvFVDWI++tRM6n9cKFcT0i/c/k8Pnjki3saILVksMzzuIWPy364G3m8ZBM7Aftn12/r8ZEIrLJN0722L1OY2/tG5gcktZlaY5W+xwlLA/JiEumYM4fVw5TAdhJAAvXeXmGSAfTKjq4Zklmhx1cZ9cnYGYdwCBiRYnQOxBrNJBX7EmaYhoMAkd0JEMIApAYswyA1I3QgraOzGiCrEW8T3MnnmFSWe6nKA4QDzQa+/ehjePXyy9gP0hBFawIhaMMwzDLde8UatdktA34QjxmgjIf08KBEO2ax6C6uTXN74pUjotuI6LeI6HNE9AAR/VBMP05Ev05EX4i/x2I6EdFPEtFDRPRpInrbcHUyTJhodvylFjl6XfnL/VB1SAE/D3QVDE3Hin0u0qqPiqSnX2VSSz2qPmGV4Vf/1cGo0ccNy1NgSDbiSeq3ccxRfIvCQBtXt9ShGWY4PcVcFciSOmbZsBrq4vMLs5Z21LW+c/kc/tb1n8WrlvYHQIqEryf63/xlxZyeF7oI2/l+NdD3L3x+tlH7DRgjamN0W5ZU254wAsAP+dtPMuT1sgHgbzLzGwC8E8APENEbAPwwgN9g5rsB/EbcB4D3A7g7/t0H4J8N1sZjkQfKSj5zPxJQkj2uy1fvH1sg9aCqy2lzWY9bcsAo+aiRYwxq4kPchL8EirCAkQIjlPObvArE0lAQ8syRHdhp5ujLOmBNbbDJ5/eh9VWdlFy9CUhJXQvNfDUjjXJidBnvu+5J3LL0ykzG+M2XqDg3qQAO5SmeDVHeRh0oPTCOIMOAFFCC0/4spIejmL/9JBPNbWb+EoAvxe3zRPQggFsAfADAu2K2jwL4bQB/J6b/LIepLr9LREeJ6FSsp6chvUMoAhjaHNZma0xP2OiYn8fIVJf+AHLRCZU5o6PRDecO2yjfIKnflJbzStTXAKEH0wQwMGAoJncGH8X6DDN0dRQgafPqcXcCXGbbM0dXZ76+ltHm25DdA4V7QqXU3kVrzQa+/8Tncfs+Y5AiYRoiQb53k32TOVhj52FLQagLx5YBmAfdrRCEMB5z28IAL6Yl9gsR3Q7grQB+D8CNCvieBnBj3L4FwBOq2OmY1g+SgAVB71d0wNjL6NUx9vsOfFk/ZIYdum0Ig1SsUQGlHEv+Qw1eykyGAJABRAdCMb2BzdOkOhxrNKa1Zo6WWcKDJjKwCRuR7XReKMExAaMB0niNTDqlNO/yLcEx5Lhl6QK+8fBp3L708j5kkFnMrJvom8yrAcmgqtgH4iMqIFqsraef4ZgxDwMiiWPORPabKT1EBoMkER0C8O8B/HVmPqffTMzM1PfFp3p99yGY4xgfOQbTkxgJINMjoHqZB76sSEe6lNMgXP3rAEgHgl2/feCowagbIC04St5UzoCizlvJo4DNA2HarzDHBHaOOfpgTAbNbMrD5fGMPgAol/cp7q/RBr7p8BN489pLlZu4v0SzyZAAhM+ExPvKslwapUUqkvGTgDKX09eeKMz/ru1vVxbR7Q4hoiUEgPxXzPwfYvIzYkYT0SkAz8b0JwHcporfGtOMMPP9AO4HgNVbbnNOqXAnDKjpI4YZwgBdYoqkjqV8bPPJ8SbWLB5asXGIQ1DQm9Tq14CkY4ImCCJ+Rc8OG+tnnGRS11lj3ZQ2bNLXncxhd8wAo1xKBf7pmbAgCZfXUBdymyWNxCpt4r879gjetLr/AVIkMEhK1zoEb+Qg4hjJuHRFBMoWmlFKXk7v9VYuPsl+YKSzMLcZCyZZFQpX9yMAHmTmH1eHPg7gQwB+LP7+okr/QSL6GICvRJieNNHUTh/egjOlFVBav2P5K88MNypdA2KT28qv5Qp7FN9jLzg6tthUGCQU2BHHqYIZGEt/ozLF1XYy0R2zs6ApIKa3remszekMfDmvAT6y7FGPibSXPoOi6YfmhSfD9cmkye8abeC7jj6Kt66+sK9NbC/hkdRvEEoPMQMYRQYofsWWOL+/4z95bDFBm+aEuD4lZsciA0peQzcoyhAm+dUAvhvAZ4joUzHtf0IAx58nog8DeBzAt8VjvwLg6wA8BOAVAN87SBOqXH/PGHWay8/m6UEGyEaDoGeQss2JTZJKF/ADGNTUTWjjQ0xgaJllE9vUPsamAoo10Kv5F7tAsVEAqgHQD/72/sg6Q7TAmCFNAaO+RbHemscjD+vRvTWA5hpt4C8dexRvXX3xmgJIkTxWEkiLXzAwArAJjuMqAURG2CogBAIQIjJGII+3bMUaAxKrnIUszO2KMPN/Qfnci7y3kp8B/MBUWlBkfzUzTfsiHQCmtPiUadPbm9NmPzHJkjmKeW2AsADIDGhNE5ieGVbjTGkDngUwWpaoZ82UjLAEOQOgQAK1xgFcYB4KRGNeE8GWy5vA0zHFVH+uV/Khsl25zbFkkDXawHcde+SaBUiRAJRs/JNhHckm+CUVEIqPUke9BVgJIfCjL7R0n9kEpWkR3d4tYQJYNKm9qVwvDG/Hiu/Rg6Nnko45hjGM1r9I0XcogGkYooBe05agSHoMYQmMeqB3Od6xtcwTeVsDlw60GJB05jH8MZWWL6eY/2HbXuIMoum4Y441xjiNrDWb+I4jC4AE5D3eoA1LUgSfYy2TM73lxaWBUcBSZvFwZJUzu8QLJrlLQkC7lJ0xkQBVpWZ2J3CM+1VgVGlEHNijmNVqsHbTWMbYNNaMHqkAzKjJjLFxgJeG1DjgTCBJspxYN2Ms2aALtsglqJjNAqY6vcYC5bi+pqTSGfZSo7I9RLQXcoU28e1HHsZb115Ec60jZJQwJlIBpVwsAMmvS+JfjFeT8srlgGKMKo8s0zcTvyQvAje7J8RoV+LC8HEYRArjmYdF8pfbdV9j3hbWKD5HMaPzbBYYUAy/rQFHMwuGwkwGPYaxmMVCNlDSF3mWfOGUuABC+zkFZQ4rZ5MwxBogWrALugW+YYFU9hpq3WWe3Mt0fkYYo0fguAhr2D/QrOPbjzyCt629gGZG84n3iwSglLuI+CVEIC9zodOELQKtXhgGUO4Rjt2IMDun5Gyq2UsyHyA5AujAZnwy4tsqbtuVx+PN9t7oZF7ngIL1M2ZwbNJ++BuNFPA1Yb4rCWMEMGriPgXQ1ODoAVEHTqpDeIAi4ALkY11AiHgMEXTG1EaTLOqbygMHm8tYa9axGj+t6lcPb2Lew6OLuNgu4cjoYjgntDg8uohX2hVcasc4PLqEQ81lMxeYwGm17DRFDnrKnEyhI2zwCC9sHMSR0RVs8BjnNldxrl3B29ZewI3jS2l+8kKsyLJp6ROziiT4NPkuDogg34dPXYGiGa5M9BlpOKN69o7MBUiujNdxx8ln8Ozlw7iyOQpskgFu41LwbYQMpt7oWtcAb2tOh/RR0xpTekRtZIttZo1gjJo4F1YBYtpXfkI9zU8DYTB7g/ixiXo4DRFjiTZBYCw3mwCAk0vnMaYW6zzCWgQ9AuPg6ApaJoyoxWqzjiXaTCA2kgUN4m9iuhroSINeBP+4fXx8Ic/DVnnr84YdQBLiAg3ACIRbll7BCHHRBmrUAg4NZjJNbp9KXl9SfftO8YHkg0xgygDFiLYyv0Nd2kSfgcxmgfM9JXMBkjesnMcPvPa38czlw7jSjgAEcHz0wgk8d/k6MBO+dOkIzl5ZC47oCJYbPIIdfxckmarKZ6jNZAHHUdMq07lVZnUGxgQMChTz6jjWlEZ6s2eNhPWNqMWR0UVc4TFWm3UsNxs40FzBSrOBJn6bZIk20RBjmTYSmFmAtQCcF49A3pZ08YPKt07idWnQmnMwC1PofQCaRdr9MIbPsscMkAEU9Yo2+Tcv57CQPslTFQMu5Vk40YsERGOK0wByGSKk/ZQeZLcl2rKbYyGir0EYkvhTzPwpIrovTl7ZkswFSBIYS80mbl19SaAADMLtB18AEG70xY1lXG7HuNKOcaldwno7wumLxwAAL60fwJMXj4JBWOcR1ttRMDDEz6iCKsISBRjFlJbtwMLEtA76NYmVlf7AJq66Mqbg/Tk4uoxDo8sAAkAeX7qQwGul2cAmU8yb69C+RB2k0T7HDILaZylAaBmjzzeKV1UzQ513ZNLg8ljzehSva2CL2bwOEXkKaR4gI4ukxCJn/QTtT6E4uDzcYxlMHi6egKUAZXbThM9oeKCc1TXfI+Mkvw/AXwPw94joOIC3bKeyuQBJAGG9OxKA5BRF48hOrhtfxkFchg4h3HnwOTCATTTY5AYtE85urOH8xio2eITnrxzCmY0DwWQnoGXC+c21arDFMMTEzAAkwACWmw2sNesgMK5fPo8l2sQSbWC52cRKsw5AlUugGraBsL1kgFb5HU2gxgJdDuhYFolYV5pJo8zpDIK5DsM2K77GxJSRgbfmexRwDC+evCZiYJCyHUb9jRYAuWWh+G94Tig/JwVTJOOjbKEDcpXhRNuRvQGS55n5DIC/RUQ/BuArtlPZXIAkAWlRUDFkmXK8NUVcgbiMVIQSCnlHaONHvoC10RXctHwu1HMQps4NHuHcxmpq9ekrR3C5HYOI8dLGQWxwk9jb8fGF4KdEixuXzyVzeDWCodWfUztmrKLJZEFQH88RbhfkUWCpjyXfpgZIYpPf+EgVYHqzWpvqMj5TPhebF3cNAGnBMZvXBMKIKgvHUjavFwC5NZGFMMzohBCnUc+XPC/Z+SR5krk+K4X2gLkN4D/KBjP/MBH9D9upbC5AEkAckoI8ThL6K3GRXQIAyei9+GiwBMKVsUqs9pHSx9RidWk91Xn90nlIEP1KO46+nQDMK80mckulSP2doFg9JoEbDYyWWRqT25niGhxNlFwBHBnGCAOimjlmd4ECx5SWQVJMa/1HFf9jDSC1yb0AyO2JuBVro+G4siXLsLHKNxM99gCTZOZfdPv/+3bqmxOQDJ3VjvgJr0INLskUVzDKREUeAU4NUwmSyD5okj4eXbEa6Temc3ybARUmm32CKOVltV0HSBQMMmtnA0MVH6UeetTJJK2fs1H5NDjmvOXQHmtea4C0eQREE4tcIOS2hfUWl+lsc2RSUZTfjhI0q/mNV02I6O0A/i6AVyPgHQFgZn7z0DrmAiSlszLyCAOKAMicoUeDR+ZQ+uGIUOR8gcJENXjmo768FHGP1cBnowaUhlGSbVMft+Y1DDD6fGZMpjpmI9uOaZr0DIhNkV4CpOgjx3WUWn9BOgztkfSFbFeYGfY/pK3WcUv9lUUJsMicjJkRwD3AJJ38KwD/I4DPYIsDmOYCJLtEgii1qVDS+bOJnKAQgKyrJyJBGA+PmVZqdjj0OdA+obRP9WO2jGWS/vsjNWbozXATIIICS2VWAwr4KqAq6Vpfcn9W7+lgT66pYd4LGSTyCYdWAWLLbSQSOY1jwKZVJINZbcffmQ1v3Hsg+Rwzf3w7FcwNSIabqaHFw0mtjO3K094/AWFAWdScjvSUy2X7c1rJPkWREiB1/VWwMgzX+zK7wWgSUHlA0+zDcpXM5fWnBmQ2iAw9adO1obgC9yJwM0nk+fPgyMwKCC04yrVP4MkWGANoXtMg+aNE9C8QPlZ4WRLV4uETZS5AkkHYRFruwbtdCqC0gRqRyXCV6mG3D5M8QN+sZx4Y7WqrAGgbc2eokcHBtqw+N+ujJLQR11ty4+BivpZDkIRhI5vEeRHW0ELUhcW8Tl9VCfrFeyB1ERC/Bc0xLQYG4mDm4AIIgZoWQBNHCrQpqBPG+OUvAy4QU4OiNqG5AEHFFKMd3UZwbDmDpPxpYGzjUzgTkNRBg70j3wvgdQCWkC8DA9hbIAkgLoIQxEa1c1o+rsVBkQFYe6yPlebjNU7lDfE81Md6F5VWnMsV/FgFmzzIBjekgFhkYsokTmyUGZvIACqrBzXxGEEWYg1Deggc08V7KL7LAGDhMwLim6QcyOHsuxRfZRinp5kux6l0eShKEz9E1bDsZ4BsSD5bYLyZ6dqUV3NvS3pG2PoUjX8xmcgKKIGULoxRs0XNGEPebJFlsCRDKrYruxXdJqLXA/ghACcRPmU99FPVX8HM92yn7bkASUYYEK6N0QIc2fzAAqAGUy11dunr1Hrk7jlpCG4G0WSqS6sxgchrI0Wzv9Afr0W7wZyDJ+oVAuRxjxBGmLpPDshQAi83/EeNHmigAJg4rTRjouJQ4yZ1OkEBpAJQztHwBL5AAEiCAU77Syp/diRY1wOlaz3JlbCTwm5DAyDSLxRDVACogVIxxwRw0CAJky5MsQaKsvIP699ZMcAZgiQR/TSAbwDwLDPfq9LfB+CfIAzN/RfM/GPM/CCAv0pEDYCfBTAUJP8/InoDM39uq3rOCUgSNlmZ2ynd79fNYxuAqTNOCy1evOk+QV+mAgBNbQaxIxvshHeltTLRiVU6NFCpsnH2Rec4y+qvBUYLhGV+HSBqInjmSLowVKhB6ZbxesbpPy2kQZQgnzIQkzyz0HytEPflUwf5GLl82hkie/nqJYdBylWOKvTWgt3XYChglhl8Pp72OadpX2INDDPo6X2yx9kCYTqu0nOZuWSSPwPgnyKAXqifaATgpwB8DcLnqP+AiD7OzJ8jom9CmG74L6do450APkVEjwJxyh724BAgANgwXCmIjz6XAEhum1ED1a75ppOCQrZL5fZUQupmIpM8o1njzELT53MziewE1calk9JMl/PHADZcnZjdgPSchyILLtgny3hWMjN3RHezVmYsqxkpyIOmYrmUQTScZwZSMcNNOQBgUk9Nviri7tTg2sU2NazmuV15nIQAmgZSefexKuOBU3zBAoC5rsqfA0JAM0lSZcnk4zSP27FGUDS99Z2dkQxnpCeJ6BNq/36/yAQz/w4R3e7KvQPAQ8z8CADEDwp+AMDnYpT640T0HwH864F6vG+owl0yFyDJkBV9LDvM3T2npQLQN14/SHnf5qkwUe7OW+oo3YnNfk1yd7JzaDvzx+hJ8kFWNPFXQz/6TUUbMcH1ikUeLKGATephhGXQJKAT/I0VNirgR1aHYHZnXfVCwsJCAyg6cFeLNJTDloKuGnwTYyYFmlJeXT5Sv54n+u0yp34eYdIkF1fyaNBDykPVNL1ieAbIfHdbtYaBB0LPGvNvfoIC25wNizQnMFmeZ+a3b6GVWwA8ofZPA/hKInoXgG8GsILwscFeIaKvAvC7zPz4FnQwMicgSZAVETUAlua1fyuSeTBRHM/58pu9GxSZ7TGq5NF12pzlWXWVSPvEqU3PicuypcMgMU8qr4wO8GjgLVcdgvFZCljrPOY4hPXFEapcMl2zcpHyv5pv9rA17+W8NYMVYNSmOyMEjzjVb88l1VOBBvmkr+XzuW39XPmy+TnregHLr31hCwc1AKkmSHB8NkWTBHRMRp9W5dd9I/ggPdtU+7MESX2iV1mY+bcB/PYURb4HwE8R0ecB/CqAX2Xmp7fS9lyAJABssOVDmjV6eLCQUD603kldcq+yLl+HqXfKB8Ob6qbK2EdTf5Ukk70Eyy4z3gCKrlQf03UoUCtWRgcnE1cApCFbXppo4mtHgDWtdxmPG+YpLJk1e9THVbow1HQNFLAK+HI+5usz14bcteDM7ruZZT6mrQfvUul6JoHyOdRMz1tLUMekjD47C7r29SaLvWg99OSKnQBJ2vlFd58EcJvavzWmTSXM/NcAgIheB+D9AH6GiI4A+C0E0Px/mXlzSF1zAZLMYbn/ms+RXb58zD40+kF2fAusEvt9NH3HunkeJ4DK+6ZWRTrFjM1lLBMd4tPUXEjjIrmGy+/f6Hq4YJm6jGZ8Wi9v5gqw2jzOpO4Ax9Qe2/oL9wBFPRMrzT219iKh0GihlzkvsLvycHnYpNSAsUtqz5h5tnnC8QSQXvuYbwBxMMCZPic6A9l5JvkHAO4motcggON3APiurVbGzH8C4E8A/AQRrQF4N4BvBfDjAAa5A+YDJCEgCdSYoeQxnIbtMZvfP0C2XqCbHXYBraRpI1uXEcCRfw0LccxWGI2uRMpovlGBY6dRKbmePvbZzTY1YAEWaE27iS2WTC/V68BRN5JAln33tiClhyPVAVEDo9ovgDuzY7tv2y5ONr7Mkva1F6K11V1tutVSiudaVVf3wZd11Vhm1q//xT+NEMNZPL1yhIjuB/BLzPxL1fqIfg7AuxCCPKcB/Cgzf4SIfhDAryEMAfppZn5gal2J/kbP4UeZefDyaXMCkhLdVje048Gosc20zwPymJbD4+gBqwsgu/WfnTnTJ308pnaGadaMzkc2X11z6zss2iILNu5mAVBReCqvYANhjxVWnNrILDSzP6t3wZQ5M8T6dM/ynCkt0UfGZ4r4cpN69SvQv0jYfv/VaOtfrfqVkF6F7IqoNoHu56vr2a6+XmfFAIdHt88y8329VTF/Z0f6r2BAcGaCXBd/70FYdFfmb38jgN+fpqK5AckwTtKb08hpEwDQPyASyNlom/SpBiB0ihXaSN+RkcV0D9AVLNMGgODXObN5AATgMo9xsV0Gg3CZx6mNdsDDwqaj6K7SVTZ3KOmEFuRKRtslFgwVM53QWYrgkNHNMUmq0Cg4MDJsOSTI2qF1UzcDLvk02Y8VE1fKqXZ8/VVWnN7IZF4I+qy6fJK1a+WtDX91hOXr1GR1sK2321/uxfNkFNbLzGRWYLvDwsz/AACI6HcAvI2Zz8f9vw+1KO8QmQuQBGgySKa07JO5eHkZBGCjDaPlnnz+OFaWNnBlY4zzr6xhgxs88swNOHLwFbzq8PNYW7qCEwdexuuvP40jyxexOl6HfhR1576uuZQ6lXxwTH4v8xLObB5IZS7yMi63S7gUQXSDR2kueqq5hgrFA2cNSm3ueZkIlJnIVfJ2AWHStgAB0hkZ0bfazUZ7TX3u0J1s2VodmTHm86hdH3EFGBcG119PhhXG/PWp5f0IMRSWPOh6YK3VOxybOu74jMBtLyy66+RGAHqx2CsxbbDMBUgy4rREGRzLhMtXljEabeKls4cAAE8+cxJnzx8MHZRCnudePAoixsXLywACWMr3hpkJgrtnzxzCF5+8ASCgGbVYGm3i1mPP4+Sh8/iKWx/G2tIVnDr8ElZH63aIDMv0OAbAWKINEIAl2sQhBaJ6rBoBeLldwWVeyufHhGc2DhvgvNguYx2j1IGrzx67321KbsVz066GSZnLFT0M2PUDb828J59QCSSVulfa8AEr+SXLr/pAzL8Yam6GPkmfVKCskh49oIOLvq3aC0m3Ph1I1mUmjxDjakS3Zy0/C+D3iegX4v4HAXx0mgrmAiQvXlrBpx64E5cuLeOZ544BDFy6vIzxeBMvnTkUnxJhkHnXCMUHNXbq9LCqKRlMAJoGV9DgkUun8DCdwu8/fDfQADcdeQkr43W84aYncdN1L+F11z+F5fEGRtRi3OjZsLpzWtNYfq9rLuE6XAppEUhvWDpr1H2lXU4D6IFwXmc3D+BCu2LyXeYxWiZc4TE2eFS9fjl80h+BzcEknVYTmpShUqK7JsNryLG5whwvG52GqQKKWaqoyiSAnCT+mumAFcMxLFddOfhna1K7b3XemHOR+ncmsseYJDP/L0T0nwD8uZj0vcz8yWnqmAuQvPDyKv7w98NCHWbeM8XZEzUzVQCRFRBCgSeF9MRWBCSJwH6/YTz93HGgAR579kYQMQ4feAVH117BiYPncfuJZ/HmU49j1LQ4tnYhLRyRO4t9FA1bck59ObZG6wCtG9A4PLpoTpMQfLXCUje5KTq0+EovtcsAGFd4jLZjcQ4G8MLGIYyIcbFdwphaLNEGLrVL1fxaXym/FU/XiFqMaTN8c5zWcW5zzbkjfJsTobs8plGqUryXiQ5sX5djeDOfy7aVy2M70uX7NLpss43BMhwkJ0a3r4YQ0UcB/BAz/1HcP0ZEP83M3ze0jrkASTDQrGd/k6SZJ0B+1RMiYEebSE+JHhJmwBAANRFYHcvkhlLdDAANcO7cQZw7fxBfpOvxycfvwC+P347ReBN3n/wSxqPw+duDy5fw5bc9gmNrF3DiwHkAMNPwRDyIapEv4eWTljIZhKXOMZVjXwnAKtZxZHQxlegKvjCA25efD9vxepJs76Ck6x11efTK9Th95Tg2e/lP/d3YKRMCFdOAJHzUvLPNbvY7ua2BzNuVLlwUXYrtkEzhk5wY3b5K8ub4eVkAADO/RERvnaaCuQBJYqDZQH7zpu9liqNHGUzqARHz2lI39SsgKvmin4jjcWojaMY5iwKg3MZFEhTAbmyMsU5jfOaV2xULBf6fh96Aw6uv4OiBV9A0Lb781oexuryO42sv49XHngfAGDUtlkebkLPoHqdou1y2GutsNL034vTGDEa2U4nnSwcxzBTB4tLuLGresfwsmIGn1o/hCtcewW7fZJLyceiQ7Pfz+TtrnyIyXLt3vfkop5R5JrUVxSieXUD5mblqvHIvSENEx5j5JQAgouOYEvfmAiQBoNmoPbLaYVVxXkVzuUhvYJknHKDq8jpvUx7P5ShvIwJTzHv+5QM4d+EAQMBjz4UA0dJ4AweWLwMEnDx4DneeeAYAMBq1eOvNjwbQJMaYWhxbe9mMPSxP03ILf6XSgGHVD70XTC6TrI5e41Y1aNXk3Xa9rQMpgXHnyrO4eekMPnfpljRSwGvS60msWLeD2vZ0vqOJoX7KoWNk/cDw/nq7zWvZr+9x9XxmCpl7zCcJ4H8D8F+J6N/G/W8F8L9MU8F8gCTDRs1Uj0ygZKhVPpZ8UZFtEqA/uWjq0YGdBI663sgoSbWrxrgbdpoyyqYCTQDY2Bzj7JUxQMDZlw/ioWdPpbK//idfBpkkvTJex50nnwkdg4Bjay/jhkPncOq6l0AAxs0mjqy+guNrF9KJZ5NeXQyVoi9j2lZs1HOfFKU3pXPt9rJL211QO0wIjAPNZdy79gQ+ffFVFaCsl7HaDWtJS+Fa6FC/D/z6pjTqZsv706VVf9sDVUbsBC6FZ4drvp/uAWHmn41Ltr0nJn3ztAvwTgRJIloF8DsISxSNAfw7Zv7ROLfyYwBOAPhDAN/NzFeIaAUh7P7lAF4A8O3M/FhvG8gX31iWZAkku0Iyo0QOUlz9NRWPT3IK4OSPwQCtAkxWvsr4nKUvo26iZKYVRspEhtTGdcGyrnJuBGxsjmIZYH19jD8+fXuuJyke9kdNi0Mrl3DdysWU9rZTj5p9ALjx4Dk01OKGA+ewPNow1xZASqNo/uvL3GXODwFPDxbTguYqreMta4/j9PpxPHL5BveVSyu1RSSGStZ1etDq0qPfR6q5eFcNFT06Fdma3jOfDbb3mCQiKO7oyuSXAbyHmV8moiUA/yWG1P8GgJ9g5o8R0T8H8GGEJdU/DOAlZr6LiL4DwD8C8O0TW/EuG1ZvfXLHcnYDPrmgK8Aw9ZuysZcXYMm5XvFdGmtGo0MNNFv5ZgtswRoQw6VFhZiAzXaEs5sHcPaVAynfE2dOeGTC8mgdRAEMidgdZ7zm6HNgEEbNJo6tXsCY2rKHUWCu77/90zgwvmLMd5m+V94D7ZUTX6LNMUmWaQOvWX4Oa3QFD1+5Aa+4YVCzkPRC6ACNPk2HTQedRnqGaXU2ZV9kpv1K/6j7KrcnBMx07vZekYkgycwM4OW4uxT/GIG+yuocHwXw9xFA8gNxGwD+HYB/SkQU6+kXfUcZQFOyRclHDgv1QgMEBqtXcqqjcfmBNAg4PWhsmknNWpNT5dUPqOyzYqUZ+SCzTRKzFaaryxhXgjREKS0zZnvdrrRLAAGX15fMMY5P9acuvdqelAHnXN+915/GuNk0mVhdr6BO/VbqMZoO/1OOLiEwTi2dwcHmMj7xymuwgfqY0N2QPjY2FDOmiq731FJtn9W2y6unt85Ehlc0L9Htbcukr10BAIhoRESfAvAsgF8H8DCAM8wsdt1phBWFAbWycDx+FsEk93XeR0SfIKJPbFy8EBIFNIAEhEkccpkApMOMQMTYMEhrX1pGmABLgFPl83/Uhj+TVx9nVb5V+21UmsnWIXXW6ml1WZ2PwvdhWyryQ/KnMvHLhHEwfvH9UUY6PqYW73v1Z7Dc1JbZC7A3fN56CBKFP2lKD8Guy3WjS3jrgcdxw/jcxDb2ksi567+2588ulzabdmdwEul5nvS3n2RQ4CYuTvkWIjoK4BcQvmO7LYnfu7gfAA5ef1u+rBG0ZKB4YdLWbL4um0LRv2SW99tVJY2EY7MVSpmSdJnWZrNNO2apWKxmlokZE5I/tTx9VbN6sfg6AIDEX9tB6l9//Cncc/xL1WNW8ve/Q72TegWpWFpQyt4yfWUYx0YXcHjtIj5/6SY8sV68X68JCVfEX50stSOzMKknyh4L3MxCBjFJkTgo87cAfBWAo0QkIKtXD04rC8fjRxACOMNE2WncKOApbbe0z4SCWZrfShu+XtL5PYu1auV9xzg16zOSMcylR7CssEjD+tQn1U0biQXLRVBt+Wug60503R7/qpsfQkNt52WzYlec6Stj3zuZJeljXkZo8drVp/Ga5eewVBlEf62LZuZ9zHQWbFTLtcgkJ4IkEV0fGSTiyr5fA+BBBLD8lpjtQwB+MW5/PO4jHv/NQf5IwJrSdaeWBUSPXIlFxQPa7+ZB0fVSo2AN5DSb9MfJJStg8y4DA8a6sGuL/KEuVNEPJStFuDzH1CJHLqfyX3/gPF534qnU6DRAWSPxNbF4rU1A2Sa1HYDy7pVncMfys2kZu4VMJxpAZ1ThsL99JEPM7VMAPhq/h9sA+Hlm/mUi+hyAjxHRPwTwSQAfifk/AuBfEtFDAF5EWH69VxK58fdR9T4/NCilyb4zpXuZTYWdVqxoq0NFtzRnPOZJ5rLUrdJTWqlqSGMKQRbvItB6cmXbXRMvKbikG9R6MEAN4y/c9iCuW7lUNDNZQk4xvyeZ3n0MsoV+a2elX738Ag41l/GZS7d2zNBZyFWRfQiAQ2RIdPvTAN5aSX8E4Ru5Pv0Swqj2qSQtVuEAyXTYms2Lyr77NaAFBxyNA6IedPCgpaPbkqH2DLEGtpjPgpdDT1WnwTcPsgrAky/S7ti2jFI57fUnnsK7b/8cZOCxxn7V7ETheLKTfZRaLBttYwrHY2KYnxi/jLcfeBRfuHwjnts4PEX9C5mlLIYA7bIkgNRAQRYk+gKredXnXJ5lIDiQBoV3fhdJte3ZZi/4ZQXqrNPpZFhgasfR4D5kYhTfXalemwro6vyS/MaTT5oB6C2o+vmDfhEzPaPvELB03ogEkPlRoFT7oeYS7l19Ep+7xHhu4/DsTMiFDJfhj8W+GQI0PyDpWV+F2RVmMoC09Jnkicdr7JCBPHvG59Og3MUqu4C1SyaBqGOhNTO8S2r5ivIV/6X3c548+DLeeetDlfoFnKYHS7vu9rSS3yQeQAHCEm3iTWun8fT6ETx0+UZc4u5l3hYye9lr0xJnIVNFt3dSUkCGUGglxwhwYOa6MCGDHJV16Hw8KV8NNFV6jdEOTas4JKs+wz4HXhWEUn4qMtZBi/HnX/0gDq+8MlHF6STclyHjKWvt+sCOj4bLwPOvOPAIjoxeWfDJqyVDgzb7zG85HyCpgUkArEHwFxJK8NNmpAMzA3yN23ft+PZ7GaRntH1m9YDzZFT8O12g6EGzp93iUFcbAK4/eB5/5tYvdNbVF2QZJrTlD1JpsBQNLHACq8063rb2OG4an9kC413ItEJT/O0nmQ+QBKoMrTCZKXYUojztMKXbenSar8fsS54ag2xcvV133wN4n1mu9aqBr8ubT6BnvyivAEOVse+FwCKvW7nYaxwPWsyhKgrceHtgKReVXRqDMKZNvH71KdwwPrcAyqsh1yCTnBufJOuVdnQ65Q8nGRapfItFkKXGCOWYbLPL79uFytNjYhvQzbhQ9216fUjvcPEpCq97se+uVzWKrUQ/uweWruBtpx7rzoziMm1Bsn9yO/WUZfVXYwgjanHv2mmc2TiIBy7dsvBT7qDst4HiQ2R+mCTgACmzRSay41KUSe7N6SowaYAs2unQYQJzLNpzx6YjTlzOGooKd9Xv9bDlysZ10rjZxIfe8js4vvoy7PyXonmD+1sTOzNnq6JZpA/oSDT+2Phl3LXyDA7QlbKChcxGFkxy9yQwRumYbAHAs7NJjE799jIrB7CAu789ICh5a+yt07fpj3nHpAPAIiil2q3KEAwi4NR1Z/CGG06jaYY9zVuHNl1DjFZPPY6yWzSIC1+9aeksDjWX8ODlm3F2wEK+C5lCGIvo9m6LgBY3sVs6xlYDrDRt2a8gDpdfsUMT4BGTXbU/0a+o6u4M4njA86AHzgDp6+gByD7JTLkbhAiMd9zysPmoWP+AHZ4ZQUhegm34KIu6zG9IPTS6jDetPoGTo/PbamMhFRnOJI8Q0f1E9I27o+jsZG6YZGFGcoUZOudULQBjAMZHwKGBJP72zLgp8up9naZcAJ1SM5OhyjgEIn9sAmhPErl8Jw+8jK+89SHry+yBPzNAf1swqefRbF9qmujLtNps4E1rp3F6/Rgev3JyMZ1xRrIHv5a4bZkPJukYnGGG2n+oVgXyTNCwtQojTOwxDi3SbRkd4Oqv6FqLpPedW/lccWqDVL6ObPV29LmjH6A1i/2ymx5Pn34Qf16fdGD4NmQ2Pso+kfpH1OLVyy/gjatPLlYSmpVcgz7J+QBJ5I5cW7TCm9oGQOO2AUfkOmoDwYEef18PcFZ9ot58d8Dsqi/PW2cQ4AZKM57U81epcAjcrIw38OW3PGLKTCqXcXoWT75Dduw8UDLCvO97V5/A0VE5aH4h08liqbTdFs+YhFFKWg2I4IDLMVKg9FcK0Ez0J9ZYKmz+SX08A6jyP6IOgkX7rm7z8FX1dk+oOj5uNvEN9/wRXnXk+fTN7b7IdvVcMAuSYFuWb4bPQjwEZ28L48T4At689kWcHC/8lFuWYHoM+9tHMj8g6c3tmpnqgLFgfLoslelV8OkyWXVZve90qe57BiuV+1esoE5JsArdCv0qVdWmI0rSoeXLeOupx9BEUMpMawhAWefuLIlC/mb49oGySy9ZfHaJNvHG1Sdx98rTGO23nnwVRN711xqTnBtvtmBI1TxW4GhgvWLaFmWROw9rMPLHa320wmy9zrX0zrTYWlFOAj9d5WssdoidLLhMwHvueAAnDpw3nywdugyFDd7MSoJiMth8VsOCaq0A+T00pk28avkFLNEmvnD5Jqzz/HxwbE/IPgPAITIXIGnM366IdB9z81IxrU1vUcBafEdH1VEb3N3LJqEYa1Ev18EOPp9raxIAmzrrpvbJA+fxztu+ENRigEg3OfSp52Jve8CpNchvx50GS5FTS2dwZHQRD1y8BefatR1pcz8KDfzIwH6S+TG3RSoMsAAUskGervGNRaTb16XqMMDaYcprcDWstEfP/MGajnxej+JAf1rfI8sUAPHPv0bmaNtS4RSGs0mdd5bM0l/KWYyjHNLmgeYy3rT2BE6Oz0/lm71mZWhke59dyvkCSQFAxQBr5ncaOlNje1KHDtaodCnP3myX+uDKueh675ChzpOy51TTyTxYDtSHmdoliyQGThw4H8ZFVnQS6BsintRuvx9kr6heRIN3mFECFpRXm3W8efUJvGb5uYWfcoBciz7J+QBJxd46o740AKw8C9Tl4crp4xqYK8cm6a7BN+uiXqs94FZly1DXwx9350r+uJM7jj+TxkVSZJZS6zRL49YIwvb6QuWErhKf827phhi3Lz+P160+tWCUE0S+GT/pD4sZNzsgtf7qACGlKX9dFURcL6gyOPn1dfiyNbbofjvrhzreJR6cyfzUyyqmaIYXOVkZr+Ndd3zOpW4NBHbW+BWxDoGdbpPVLxHjxvE5jNdafP7STbjIyzvc+h6V4Y/PvplxMz8gqUVYEncDU2cPkt7lrM+JJmuiFa5sTYcJulefpEp7VeYq0nfuCRz7VXnTTU/gVUeet4msEDZVunX2tBNgxkD4eiSs2c0zXByjJkSMk+PzWFlbx6NXbsBzG9ftWFt7UvahKT1E5sPcBqD9iFWTWkSBhwnMwObVdXQNGpfgT9HJfVmvg2KghWmvsMe0Xam/k31qU7pmopu83Szy3Xc8gMaATMhvgWaa6DanrZ2RfGGSs0IFcXYCIGsAf93oEt64eho3L7008/b2vFyDgZu5Y5JVs1pLzcT1INqoPHKcTBZbV6WePr9mtX1hkFq/mtSOd7TFOiPZLKF914oy1f/8a/4Erzr6vMXtNN5p+qeYkb+gOOHSzFzCddgZFtlV45havHblaSzRJp5aP7YYT4n4mO8zABwicweSvb69LrNXgJGRzWXP4GJaV6CkMPE9SHexwUlSYYhJH3cOqIF7TSbYuCcPnsN77/wMRtS6ZixAThOkCCVz+Z0Bx7o+hJ1hkZM0GVGLu1aewYnRy3j0yvV4afPgVdVhHoXaaw8l58bcrg2h0RHvqsns02S7Nvyntu32CyvUAWSnNaHGQhbR9eKkenTvM1VqLLd60YDXHH8Wh2NEu1J8S9aQhtadY49dbxDs+LjJLiEAx8YX8Ka10zgxenlXdJgbGWpq7zMcnRuQBFBdwiz55TwQdpnGXvr8i+rXgHEtT0c5jaxdc7/NOfWgVa8PtDgv50WPeY4dvIBvfuPvVwqEBqlzb7LUPJI74yXMfzoCDVxdsNQtLdEG3hiB8loeJjTFEKB9I/MDko6tSdok36A51sESC3BSoNU1btJL8g92mN9TzZjpYsW17aHCwSQ9vvYyDi5frre5jep1aT0Efecgy1ISwta+4z07DcICGfeuPYF7Vp5Gc60OPL8GmeR8+CQde+zME6VvMQst1fo6zPoiT/RNJi+cATHHpvrALmY1JTzCKNbcu1BH2nYMMoL3yQMv4/u+/LcwUq/yxMSLarf3JGvDeGegy3Leeel3Y2pxy/KLWG2u4MFLt+DyNbbi+bUYuJkbJtkJaA5Aqyatr4NUIKcjb2dbsOUKgIRLmIQQlMG22pZn0F5n0aOooNT1DTc8gWNrFyqq8sze7t783Vlu5w18SsOCoH53Q06MX8adK89cW/O+w8Uf9rePZG5AsiZ6PGQxqLsLLEl1Kh/AUfV21qXTunyBqlPwhPzVRTh8gMi3XTnW+wanMLXu7bc+UhwKz2zZjf2CFZPEn6IvuTNWVvdN2k2AFDm1dAZvXn0C119DQLnwSe6mONDzC1AMmfWSzNVBJqskuHzkkjU6uCBN53Ae3WbPseo88w4mmRhhxwoC95x8Cqeue8k2Q/XtrUgXBGhw3BnYKt8o8wRHRIw3rD6JVy8/v++BUrrAtbbAxXw5VDxIwIFe5Zgp0tNLZQzkpHypzeJGl3e+yObqrQJyBeyrdfi6JFFotZKTB87jQ1/+f+PA0pVCRwDpcw0Wf6d/kkX9qwcG0mJu2d8HPxtnp6cu1mRELe5Yfg5LtLm/v8y4D03pITI/TDKK/iIigNwvPMD0mNtdvso+5icr7kwyw81KRF35UAHISr4B5FjV1fFwEvCO2x4qxkUC6nmWedD62BZ4X7580SdYpM9aaleIivb0d3KuNkAmrYjxquUXcNfKMzjUXNoVHa6GTMEk980qQINBkohGRPRJIvrluP8aIvo9InqIiP4NES3H9JW4/1A8fvvQNgZFtyNoCkB1lSnGPep6OoCQ/EZibbUQtS3rjxcf7fK62dPprzellQCwPNrAG288XSmAuDRaeGprVW5lnKQYvBp0Nc/bOckho8K/GtnjPHCcU0tn8Ja1L+7fLzPywL+4ChAz/9LuKDo7mYZJ/hCAB9X+PwLwE8x8F4CXAHw4pn8YwEsx/SdivokyDZh1jXs0gFfzjThi4v2OaT/1+g4fpK6rou+kWTfacKy6CpTeZnaN04kIeM+dD+BVR91KP6Y6im1ZI3naLyXW1NxNYfcLzI9uK8063rL2OE4tnUnz3feLXIs+yUEgSUS3Avh6AP8i7hOA9wD4dzHLRwF8MG5/IO4jHn9vzD+gIaCgJArANAMzx2t+TPXbzcx6gGqoeADvaMfo59ucAPhdsjTawNtvfdiMi7R6BCD0ZvFWwVF7B68+JHXfzERg5iDiLTKiFq9beQpvXD2NJdrcbXVmIwxgk4f97SMZyiT/MYC/jfxF3RMAzjDzRtw/DeCWuH0LgCcAIB4/G/MbIaL7iOgTRPSJjVculEABy8aqj3+l31THU3ZV0AtGXNWjq91YYjJ2TGKjqq5Jff7NN30Rp64703k8D/0hA4zTDv8RfXT53RUqb8EcjJ300hDjhqVzeMPqk/vGT7lgkhUhom8A8Cwz/+EsG2bm+5n57cz89vGBcnUVHblOzKvCsrhjexIj42LDlYvHemfteLO47vhLZQZPsZQsRkmr6OrSOt5z52cHsELyqm7Z1PZ1DM07W8kcudZOinLvWPvTy8nxebx17XEc3w8LZCwGk1flqwF8ExE9BuBjCGb2PwFwlIhkrMOtAJ6M208CuA0A4vEjAF6Y2EqHWVocc/lkkHnv+Ecpp4CMJI/xfWZIrvogK0AoJaqBGge4E/XzWUgd1KyagK+49WG86mj3ZU3Nsm3MM8qtiobt3eduaiEMWaR3jhglACw3G3jj6pN7HigXTLIizPwjzHwrM98O4DsA/CYz/yUAvwXgW2K2DwH4xbj98biPePw3mSe8WhyoVFmj/7phVz/wdljFZ1hdeq3vzk4y1b2vlFGWUcy46od09ZriTrfV8Tq+9rV/jKZnasNOPKfk9rrwv8y7U1poT6ljunPWUwkBKO9dO427Vp7ZmwEdnuJvH8l2Rr3+HQAfI6J/COCTAD4S0z8C4F8S0UMAXkQA1l4prqvvBxXK0ul7VAClTXYDUL6Aa6+4x55xquRBaCD6VHQsRJ9rR0d/441P4MjqsCEmFP+tmcfTmtxc2dIpV5e7mbsKqPYDeY5nPidgKVos0SZevfw8Vmgdj125HhfalV3VaxohALTPgjJDZCqQZObfBvDbcfsRAO+o5LkE4Fu3pI0lBjYdPRaURkKVpncH39qImdXIeGxDY2bhEvD61/L3nIdlpOIbyKnL4w28967PViPaRXXRF8GMbU9LrMnVB0Yt5Q2qjZTUgZx5AUwAuGnpLA6PLuLBS7fg7ObaHATDhgntM3/jEJmfGTeCIB2sDZjg6+jwafaPbbQVdg5A9+Un5XM6GIZTMcOlb1dPzy2s+6abvlgdF9mlMCMD5OymE9YdtLvbffYKzGQ50FzBW9cew6uXX5jhvdlBWZjbcyA1AFHSte5jkX9SUMQsKVSpu2buK0bY5XO0pjIKS77z/Cj/WLdD6Yt8z50PbKFDbS2a3V1bANwJt+uqS9U3qqYs7sa87knSEOOO5WcBAE9cOY7NOeItpTD2W+R6iMzPHXGAoX2OndML9S8Q7qEq291ryQBQ70o8XB4qgK8LMJUvsvpoccduR4jw3hufmIJFxqrS4haErYyPLOpTW17Dq08kbGvZ/6pyOHN7nsZRilAEyrcceByrtL7b6vTKIrq926Ke3wrZy3nIbce/tLit+PZqjDIN84nsoquNWntSt/dTejfBpIekw52Q6/TOycA4Xn/DkxN9kdb3yOpSsfndqnT5IbuI+M6Cpn0YpJ1am/MKkCJEjKOjV3Dv2hPzPe97MU5y92WqGTOVLAkcJwERsf3Y2IA2Jq2eDmlbl+mr17FlV8rIzYdfwpederxfQeTn0zc5C49d3y3p6haTbsNspE79u3Sa5y58ZHQRr199CidH53dblVIYoE0e9LefZO5AEhhoWtfKSZ4arYkImsz3WjkvKq1gj7VtlXfighj6UNV3mVOuW7mEb3z9H+LAUuUDX13VV2yeWU1NtKdv69gdrlYzLaxc7Y+IbUcONJfTeMq5k0XgZhfFW5hdZmvHs27WivSsTPabDjN+CPh25a/URaKP16VyThZIFYorkDu8ehF3nXy6Uz8ia+GQ2bLAuF2zO9fIVZDtMsevrpRsMsHoHnGYjajFrUsvAgAevnzD3MTur8UhQPMDkrp39TFH1wvZp6eyFTT0SR60KsA2hGWastwBkEOkMi6yIca773gAa+P6quNABSBTeftAzyLCXbp5d7PT9F1k/4KQjRzh9j7KeQPQEbW4belFHGwu4/ErJ3Fm88Buq7Tv/I1DZG7M7UHWkARoUGH1GhxR2fZtdAHZVl7Y3tqblKemD+phwbfc/Bi+4raHB6vCkM5PcX/2DIRSvTzjwUVb0aQrXb81ybqN4/xuDYrzBpAiDbU4OT6PN60+gWOjC5ML7KQwwjpgQ/72kcwNSBqpRqXtIdLpxl9IVbAZxAbZblZZZBe75Q7QHioVkB+PWrznzs9iqdmol3HF/awagcms5mwAU0ELyg9D5Dy7Czva/9K9ajkRzy1AalluNvCmtSdw18ozu/ZaIjCIh/1hH32+YX7M7Y7ASBrWo/P1+irrD9BQAEtdPmJt1U86oXwnkKpjOrLOkuBm19x+7DncdN3Z7sYqehOFh9n74sI2G7/kdqSJbTD6Te7d91Hm1r0XlZ3pPe9guUSbuHXpRazzCE9cOYF2N65sO5gmnmXm+3ZSlasl8wOSWoQduqivHiRegp57wDWYelNXW2Js8+pi1eh05bk0QD7J3HbbEgQpdAPwdff8UW9EuybMiJSy3uFnAZAahMtj9jcD6W5J6XTW+nB8Oc07QIqMqMVdy89gldbx+JWTuMRLV69xMbevMZlfkIxPssOwDtFoqMzLaQCrA8SqTcACYy8ATFa+2vYt172EVx8bPrsmzc8miV7vpIiRXZ7Y7jPHPrHgLgPM9wKLNELArcsv4vj4Ah6+fAOe2zh81QzwazG6PT8+Se9f9ABZ+AbZ+SM7opl9SKai0VX22KWbbLKppl7Ot+1ZcKVMQ4x33dkf0faiJzrkzzZIlcJXtw9fVOz1c1PvMr56ot9qpRY5wBWu914SOZMDzWW8cfU0XrX8/NXzU16DM27mk0liAoMkKCrHqoQr5P2cnrHFtML32KeQ16NLOvyXBsblPJy89ebH8I7bHpqkUYeeFFllUGCWfkjAx4y5elk8r98dkKw5uf23fqKOcz4UyIvWTi+Q8cUrJ3ZkNINpeZ8B4BCZD5B0vcmwuq5eVjsmLK3me9Tl3LaDWVt/jxi26nWZiA71BhowXn/DaSw1039hL0e3WV2C2Xca+w4SP58eqK50mnnrQ0XfADKpNZl3YOyThhh3Lj+LY6MLePDSzbi8U35Kxr77EuIQmStzO5mtfYEQPVDaj4Os1DlN+51g3He82nBHtoJylUB56vAZvOXmx4dVqNUjbemQArLZPtQ1fsZD3ihXXSY7pLXWtcUv5nlBDC9EjBPjl3Hv6mkcbKYL9k3VzvAhQPtG5gYkEwOMT27386ke7Zi/dym1WlHx3flyE+rQLHUqQCaUQ5kq5VbGG3j/PZ+cOqLtNSLHJMXknoXY86aUtnfgJDPgrq4s4LgX2eXR8St4+4FHceP47M40sPBJ7rLUzGMR6nHwDQVIZLCauPhET89PE1r6XAGVukpgZVPHiQPnB63049uiiN4UfbT5MubKZ8kofQikzyep03YfSIMW/kqk/Y43816Lfo9pE69ffQqrV9bx5JVj2MBoNhUzgHbvXIdZydwwSQDdfkDla/MZOhmaLkt2d2LbHUBt1Ot7VrrAO1WkqGzM0xDj3puemDrSGsZFqmAEZ1CUKYOzWGxXJHPHMtVu1e7WbkvpozSparxk+jQt9iajHFGLu1aewetXn8J4ZoMbB7LIBZPcYdEsDagySJmcktae1GWlZ/axwOoBVUeXXprp1uMCneWq+K86382HX8LXvvaPp2J8afUfzt+yEZBtIfOVZ/vASmgmXwJ7Zrq1+QjgAN1RtaC76NyyXdV8L83G8SJ35YalcwCABy/dMpuK9xkADpH5AckEOGz3dRoP8O31AJeQLq4d74s9KFO9KhpgXX+stuXaaYjxrjs+N3hcZHoXqBdJ7XRq4yNnCZp759vRQ/wxQVKUPg4032vgKKK1vmHp3GzGgjKAzWtvys38gGQSx0wIPohdAk+XT9Edo2LD7U8ATrM7pI5en2WuOES0H6tkrEv5uAfN+gntbDq7vgbS4p740t9gUa4QlCxyrwLnyfEsVjpngBcgubviwANQ1jZXmJxnfxXGpkG2czZNF5hxZZNUQhcgdupb71x/4Y4HthTR1oySyMawOeWYLYip25Hq9pdi9wM0XaK19xrnY12el700JGjHZGFu77bUb0CnedxV3IPXJJkAnhMBsqfvkK/DybG1C7jr+NaW6U+nKy8XJjBxisZ1LWO2XdEACZSw470PwLwAZ43u2zebZcqyodwVe5BFzkwY12R0e75A0tnWsjuRCXaJDuxM20u7SIevqyNPdaiRy9cQ4+te90ncdN2ZqVTTn2swS6ORBaZKkzORDIYS6ODqJdF551vKEFSVb17LACmyYJK7KNUPVwWZaOV0AZU+NrU+A9vuatc3Wzm/mw+/hLfd/OjUqunnNC+Nllsksz37h9oDYJdvcmd47Cykbm7LFatNr5RS8vkH/z3va0YWIDk/oteOTNIXXXbmMUlSHxOc5n4PYZWV7EVDii2//oYnsbY0fKUfwcLyOa21tPNBlS637PyC42Txa2Wa86i8Ma85gNycfk2BvS7zBZIKWQoC2GFqmw6p9nsZoEbRvvq8HkMAu2Ly2pTwe3D5Mr7qVZ/vUbKitqpUTG6i0lWaT892dy7Stie+DfHozb+53X3zbNgmbNfOZc+uRbldWTDJXRYd5JhkJsfj5BFCyk/B+nT7ia1NKqfq96u2sc+kpx9SaOODb/yDqX2RXvQ0RK1mCZAh3ywBy1/2vFczx+cRLH2oqfRI+kdQA6Pev6ZkAZK7K72LVAwJoHQd83VMMMGn+qqiYGAEv5JBlps3H34Rb7tluC/Sf1c7KZpM6rB3NcHIt54X0lCujquky9aky+dSZ5CMckZOcAfP91nOVngR3d5NmbjghLUji7ITn9UBgMr6WM1mrPV8cQ902ZcVxb7i1kemXnXcVKkYjABzF2ub5YK7RgfVgpjxXo+9IZPhvHDByPVP00H39uycwcIALwaT757oudidDq0J9tt2LJ+ig9fa6HO0keIfBqXIpB1YvoI33fTFrSsKFAsvlB6HPKx8pwI3dcO0Djjzb5D6t6A9j8IQdxHulgkN1Vdp33eymJa4e1IsP1azgjokzceeZJJXojIcy1d9obUoTpfu1UrVbyz+9a/7oy37ItPl6SA/+fLtDHvU7ei9PvN674GGBcwaXIZnjWzqHvs87ZaEeZpPyu4bGbRUGhE9RkSfIaJPEdEnYtpxIvp1IvpC/D0W04mIfpKIHiKiTxPR26bSaAhAquNp4dyKGe7zDuqx04Zlu0z9Cq09dfglfOVtX9gygGkfGCMzyhLTr14HtZe3u11Wf3tDxIkAaK6+d/TfIbkGl0qbZj3JdzPzW5j57XH/hwH8BjPfDeA34j4AvB/A3fHvPgD/bFDt3qbRL/RK3hSo0Agx5b1J2bsAsXZ1Km1YfyrDTBOSZALuuf4pHJhiXGRdSmX9pdrZj0FlyberHzw6XMlzLKVn12tvHrfKep77VbhtB/3tJ9nOorsfAPDRuP1RAB9U6T/LQX4XwFEiOrXlVri+W5h4XWyytk/FZnfbut4OpxP5MhU5vvYyvv51f9TX2iCRFzWRmN3laMhZfq6hqkOxL55Qvd9fZv6ky+QoQ1HF+3vC1xb3D3AOZJHXKJNkAP+ZiP6QiO6LaTcy85fi9tMAbozbtwB4QpU9HdOMENF9RPQJIvrE5oULnQDU6QeEMrO78noyoLe7nltf5wA7sTjkEhpivPNVX8Da0np3JT1CgPoSojapS8Yajux8p6x5JfJ7pJzLPf/S56vpZu8mt1pezeTdL/5JRhgCNORvH8nQwM2fZeYniegGAL9ORH+iDzIz05RPAjPfD+B+AFi95bb6qrEdLFIwLk05rIGe35fXgQ7S1PJ5qRxPi1cMzH/z4Zfwnjs/uy1fJBh58LhfQbuz3M4FcOT8vR+05F1lub0l3lnQvZSaDtzsR2EAvEvTEonogwC+HsBhAB9h5v98tdoexCSZ+cn4+yyAXwDwDgDPiBkdf5+N2Z8EcJsqfmtMGy4dwFcQQ233TOqdfVIro01sZ+PrwHXRHSod5LXXPzXVHO2aKoBYMtQLNHroD/Uwoe2Kv2QetLtI/N6QmnPcppvzrXwXZ18KM8DtsL8BQkQ/TUTPEtFnXfr7iOhPY/D3h0PT/H8y8/cD+KsAvn3m59YjE0GSiA4S0XWyDeC/AfBZAB8H8KGY7UMAfjFufxzA98Qo9zsBnFVm+XBRz1qyeLV/0PsgNa2psT/5JdsFOnuwB8eYlvSI6axrr5DpU9edwfvv+VRHI8MkqUEc2KQbH6mlnEuta5itiCeytvqPZv2Sd++If6hKMedY+WCYNr333kuiW7jlQX8D5WcAvE8nENEIwE8hBIDfAOA7iegNKsvfi8evmgwxt28E8AsUbL0xgH/NzL9KRH8A4OeJ6MMAHgfwbTH/rwD4OgAPAXgFwPdObEEzQuc7NOMnh9TRdbjCTA0brdXVwWbLruPopZK7Tz6Ng1v6jnZFOA9JIYXwWk39nW39jZudG1RO0HC8v8Q/eKWTw6wGX5nTve/GTs5wxg0z/w4R3e6S3wHgIWZ+BACI6GMAPkBEDwL4MQD/iZm3HwGdQojnIBJFROcB/Olu6zFQTgJ4freVGCB7RU9g7+i6V/QE6rq+mpmv32qFRPSrsd4hsgrgktq/P8YhfJ23A/hlZr437n8LgPcx81+J+98N4CsBfB7BYv0DAJ9i5n++1fOYVuZlxs2fqvGXcy1E9Im9oOte0RPYO7ruFT2BndGVmd83OdfOCDP/JICf3I22tzNOciELWchCZi3bD/zOWBYguZCFLGSe5A8A3E1EryGiZQDfgRAM3jWZF5AsfBVzLHtF172iJ7B3dN0regJ7QFci+jkA/xXAPUR0mog+zMwbAH4QwK8BeBDAzzPzA7uq5zwEbhaykIUsZF5lXpjkQhaykIXMpSxAciELWchCemTXQbI2BWkXdSmmSe3Yupnb1/U2IvotIvocET1ARD80j/oS0SoR/T4R/XHU8x/E9NcQ0e9Fff5NdNKDiFbi/kPx+O1XQ0+l74iIPklEvzznel69NV6vcdlVkBwwBelqy8/ATZPCrNfNnJ1sAPibzPwGAO8E8APx2s2bvpcBvIeZvwzAWwC8L05X/UcAfoKZ7wLwEoAPx/wfBvBSTP+JmO9qyg8hBAxE5lVPYKfXeF1IEGbetT8AXwXg19T+jwD4kV3W6XYAn1X7fwrgVNw+hTDwHQD+DwDfWcu3S3r/IoCvmWd9ARwA8EcIMyieBzD2zwFCVPOr4vY45qOrpN+tCODyHgC/jDD3cO70jG0+BuCkS5vbe7+X/3bb3B609uQuy7bWzbwaEk29twL4PcyhvtGE/RTCSlG/DuBhAGc4DPfwuiQ94/GzAE5cDT0B/GMAfxuATFA+Mad6Apj9Gq8Lqcu8TEvcE8I8/bqZOy1EdAjAvwfw15n5HKnVeedFX2beBPAWIjqKsNTe63ZXo1KI6BsAPMvMf0hE79pldYbIzNd4XUhddptJzt0UpIrs3LqZ2xQiWkIAyH/FzP8hJs+tvsx8BsBvIZitR4lIXtJal6RnPH4EwAtXQb2vBvBNRPQYgI8hmNz/ZA71BLALa7xew7LbIDl3U5AqsrPrZm5RKFDGjwB4kJl/fF71JaLrI4MEEa0h+E0fRADLb+nQU/T/FgC/ydGRtpPCzD/CzLcy8+0Iz+FvMvNfmjc9gV1c4/Vald12iiKsPfl5BD/V391lXX4OwJcArCP4bT6M4Gf6DQBfAPB/ATge8xJCZP5hAJ8B8ParrOufRfBLfRrAp+Lf182bvgDeDOCTUc/PAvifY/odAH4fYd3RfwtgJaavxv2H4vE7duE5eBfC8l1zqWfU6Y/j3wPSb+bt3u+Xv8W0xIUsZCEL6ZHdNrcXspCFLGSuZQGSC1nIQhbSIwuQXMhCFrKQHlmA5EIWspCF9MgCJBeykIUspEcWILmQhSxkIT2yAMmFLGQhC+mR/x+LOznCb+2BkAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Using viridis color map\n",
"cmap = plt.cm.viridis\n",
"# Set up a logrithm normalizer with min and max\n",
"norm = LogNorm(vmin=1e2, vmax=1e4)\n",
"# show the color mapped image\n",
"fimage = cmap(norm(lum))\n",
"plt.imshow(fimage)\n",
"# add a color scale\n",
"plt.colorbar(plt.cm.ScalarMappable(norm=norm, cmap=cmap), label=r\"$\\mathrm{cd/m^2}$\")"
]
},
{
"cell_type": "markdown",
"id": "6802890f",
"metadata": {},
"source": [
"## We can also use the external library to denoise the image.\n",
"## Here are use a library call scikit-image to denoise our image using a wavelet denoiser\n",
"The code below are taken from https://scikit-image.org/docs/stable/auto_examples/filters/plot_j_invariant_tutorial.html#sphx-glr-auto-examples-filters-plot-j-invariant-tutorial-py"
]
},
{
"cell_type": "code",
"execution_count": 40,
"id": "ae6c4ed5",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Clipping input data to the valid range for imshow with RGB data ([0..1] for floats or [0..255] for integers).\n",
"Clipping input data to the valid range for imshow with RGB data ([0..1] for floats or [0..255] for integers).\n",
"Clipping input data to the valid range for imshow with RGB data ([0..1] for floats or [0..255] for integers).\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA1YAAAEYCAYAAABfkWdeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACcfklEQVR4nO39e7BtSXoXBv5yrbXf5/2+91bVfde9VV3V1Uy3JJAausF0S91qDYORZJvx2Mgv7AiMscFgxg9ADiaYsDEPg8MTY0PIyB4LhIfBPMYeghEhZuwRBiTT6BF6gdQgdXVX1a1b933OPmv+2Cf3zf2dfHyZa6291977+0WcOHutfH2ZK/OX3/dl5lqqLEsIBAKBQCAQCAQCgSAd2aIFEAgEAoFAIBAIBIJlhxhWAoFAIBAIBAKBQFARYlgJBAKBQCAQCAQCQUWIYSUQCAQCgUAgEAgEFSGGlUAgEAgEAoFAIBBUhBhWAoFAIBAIBAKBQFARYlitAJRSf1Up9c8vWg6BQLDeUEr975VS/2MD+ZZKqTue8N+qlPqjNZanlFJ/Win1gVLqRyrmdeNC/sIR3lNK/aRS6rBKOQLBOmMVuEcp9Vml1FeM67+vlPrsxe/fr5T6/jrKqQtm2yil/rBS6l9btExtgBhWLYBS6h8opd5VSo2Me/+SUuqHOOnLsvxCWZbfV7NMXjIRCATtxAWfPFVKfaSUeqCU+v8qpf5VpVTjfF+W5X9dluXnmy7HhFKqC+DfA/AfXVxrQ+bRxd9XlVJ/SSn1uYhsPw3gcwBeKcvyG2uW94eUUv+Svi7L8jmAPwXg36mzHIFg3lh37tH3Loygn1ZKPb5okz+llLoRm39Zlh8ry/KH6pP4JZRSv0Up9TdrzPI/BvB/vGiTtYYYVu1BDuDfWLQQAoFgJfAdZVluArgO4A8B+D0A/svFitQYfgOAnyzL8h+R+ztlWW4AeAfA/wvA/10p9VuYeV4H8A/Ksnxcn5he/DcA/nmlVG9O5QkETWHduecHAfxvAfxmANuY8M/fBvBPzEso1+p4kyjL8pcA/CQmdV9riGHVHvxHAH6XUmrHFqiU+mal1N9SSn148f+bjbCpB1QpdUcp9Tcu4n1dKfUDF/f/pFLqD5M8/6JS6t8MCXbhfflzSqnvv/BE/T2l1OtKqd97sdL2i0qpzxvxv0cp9RMXcX9OKfVbSX6/Wyn1S0qpf3yxMmcuJ/eUUv+xUuoXLjzN/7lSasBvRoFAoFGW5YdlWf5FAP8UJor7W4B/nOntKEqp33kxvn9JKfU9Ok+l1LZS6r9SSn1NKfUPlVL/nvZIm15QNcEfucjj4QVvBMu/CP+3DY74FwLV/AKAv+Fpg18uy/KPAfj9AP7PhqxXlVJ//qIeP6+U+u0X9/9FAP8FgF91seL1B5RSuxerXl9Tk+2Bf0kp9Yoh7z9QSv1649q6bUcp9QcB/GoAf+Ii7z9xIeNXAHwA4FcG6ioQLAXWkXsuOOBzAH5DWZZ/qyzLs4t2+JNlWf6XF3G8+pEJyisA+kqpH7hI+3eUUu+QuL9HKfW/AnislCqUUv+OUupnL+L/uFLqN17EfQPAf46XHPegprb5IQDfHmizlYcYVu3B/4JJp/xdNEAptQfgLwP44wD2AfwnAP6yUmrfks9/COB/BLAL4BUA/+nF/e8D8M8YJHQA4Ndj4inl4DsA/JmLfP8ugP8Bk/5zDcD3Avi/GHHfBfAlAFsAvgfAH1FK/W8uyv02AP/WRdl3AHyWlPOHALwO4BMX4dcA/AdMGQUCgQVlWf4IgK9gotQD4XF2gom39RqAfxHAn1RK7V6E/acXYbcAfAbAP4fJOKf4PIBfc1HONoDvBvBeqPwLjvhdmCgodzHhCh/eBvBTgTgA8N8BOAJw74IH/3sAP3ZR9j8B4Hcopb71QgH6VwH8T2VZbpRl+fsw4bo/jYkX/jUATwH8CUaZMyjL8t8F8MMAfttF3r/NCP4JTLzbAsHKYM2459cD+JGyLH/Rk8apHzHwGwD8OQB7mOhuf0Ep1THC/xlMDJudsizPAPwsJu2+DeAPAPh+pdSVsix/ArMct3ORvmrbCIdBDKu24T8A8K+ry4eYvx3AT5dl+WcuPCD/N0yWXL/DkscpJpP/1bIsn5Vl+TeBKbl9iJfL0f80gB8qy/KrTNl+uCzL/+FisP45AIcA/lBZlqcA/lsAN9TFaltZln+5LMufLSf4G5gYeppUvxvAny7L8u+XZfkEEy8ygImXCcC/AuDfLMvy/bIsPwLwf7qQVSAQVMM/BrDHHGenAL63LMvTsiz/CoBHmBgk+UW831uW5UdlWf4DAH8YwP/BUt4pgE0A9wGosix/oizLX2KUrzniyxdb8X5/oF47AD5i1h+YKCXfAOCwLMvvLcvyRVmWPwfg/woH15Rl+V5Zln++LMsnF/L+QUwUuzrxESZ1EQhWDevCPfsAfsmXIKAfhfC3y7L8wQu96z8B0MfsKvcfL8vyF8uyfHpR1p8ry/Ifl2V5XpblDwD4aQDWM6M1tY1wGIC578MUuFGW5ZeVUn8Jk0PMP2EEXQXwD0n0f4iJN4Hid2OyavUjSqkPAPzhsiz/1EXY9wH4ZzE5b/DPAvhjEeKZBthTAF8vy3JsXAPABoAHSqkvAPh9mHg+MgBDAH/PqMv/YuRlenYOL+L+7ckYBwAoTM6fCQSCargG4H3wxtl7F04UjSeYjO8DAB3M8pGVi8qy/OsXW93+JIDrSqn/DhOPZz9Q/lVMziSY+fvwASZKVAhaxvcx8TRf1VtgLpBjspp0CUqpIYA/AuDbMFm1B4BNpVRu8GBVbAJ4EIokECwh1oV73sNE73EioB+FMNWXyrI8V5M3CF61hV+U9c9hskPoxsUt3Y42hJ4Np22EwyArVm3E7wPwL2OWLP4xJqtQJl4DQA9r6/ME/3JZllcB/FYA/5l6+Xa/7wfwGy725b4B4C/ULDvU5PD1n8fkDTHHF0vMfwWTAQpMvDmvGEleNX5/HRMj7WNlWe5c/G1fHEAXCASJUEp9Ayac8jdRbZx9HS9XxTWsXAQAZVn+8bIsPwngTUwUiX+bUf4vYZYXXgvI9L8ioMxc4Ddisg3npzBRQH7eKH+nLMvNsiy/6Ej7OwHcA/BNZVluYbLNCHjJa48xUUo0TjxylI77b2CyNVEgWBmsGff8NQDfqIzzlyYY+lEIU9kutjO/gpcr8YDBLUqp65iswv82APsXZX3ZKIvyUB1tIxwGMaxah7IsfwbADwD47cbtvwLgdaXUb744kPhPYUIWf4mmV0p9lzGoP8Bk8Jxf5P0VAH8Lk7NSf14vF9eMLoAegK8BOLvwzpivQP2zAL5HKfXGhRf439cBZVmeY0IEf0QpdXRRn2tKqW9tQE6BYOWhlNpSSn0Jk+2631+W5d+rMs4uVmf+LIA/qJTavJi8/y1MnDa07G9QSn3TxRmAxwCeAThnlP9nAfwWpdSbFxzx+wJi/RV4tuUppY6VUr/tIp/fe1H+jwD46OKw90AplSul3rpQAm3YxETpeHBx5pXK9KMA/mmlVEcp9SkA3+mR96uYnBExZbyGyRbF/9mTTiBYGqwj95Rl+dfw8g2kn7zQ1zbV5JXz/wLC+lEIn1RK/ZNq8ta/3wHgOdycMcJE//vaRT2/B8BbRvhXAbyiLl6PXlPbfAbAX42oz0pCDKt24nsxGRQAJvv7MTns+DsxWWr+3QC+VJbl1y1pvwHA/08p9QjAXwTwb1ycH9D4Pky2wfyZJgS/2Jf72zEZhB9g8srRv2iE/1VMXsLx/wbwM3hJCs8v/v8efV8p9RATD9C9JmQVCFYY/71S6iNMVmb+XUz245uHvKuMs38dE2Xl5zDxQv83mHyHiWILk4n6A0y2jbyHl997cZZ/wRF/FMBfv4jz10N1BXBfKXWV3H+glHqMyTabLwL4Lr0t+kJJ+xImh7R/HhNv7X+BySFvG/4ogMFFvP8ZwP+ThP/7AG5f1PUPwP9SoD8G4DvV5O2Cf/zi3m8G8H3l5JtWAsEyY9255zsxMbh+AJNz7V8G8CkAfy2kHzHw/8DkLYsfYHK27J+8OG91CWVZ/jgmZ9D+J0yMqLcB/H+MKH8dwN8H8MtKKa1LJreNUuoKJg7/vxBRn5WEKkvXrgTBKkIp9Wsw8fBcL1vw8NXktZ9fBtAj+6oFAoGABaXUvwLgzbIsf8eiZYnFxfagHwPwa8qyfHfR8ggEAj6WmXvqhJp8zudny7L8zxYty6IhhtUa4WJZ/L8F8GNlWX7vAuX4jZh4dIaYrKCdl2X5v1uUPAKBQCAQCAQCQVXIVsA1wcXK0AMAVzBZzl0kfismh8h/FsAYwL+2WHEEAoFAIBAIBIJqkBUrgUAgEAgEAoFAIKgIWbESCAQCgUAgEAgEgooQw0ogEAgEAoFAIBAIKqKIidzv98uNjXq+1Wp82XkhWPfybVBKoSzLVsrWBNalnjaUZYnz83MAwEcffYSnT58udWMIN61O+TYIN60PVo2bgHr5iWKRfWVdy/ahrXLVgVWumwv0uJS+LsvSyU9RhtXGxga+4zu+wxpma3A9GfrCuXmlpOF0AjNOKH5s2b66+9KlyGCGuRQQbrgvf7OT+eqly3C1gU2WkOJEy+fInCJbHWgzAZ2dneHJkyfT+v/gD/7ggiWqDuEm4SbOeBZuEm5aBHz8BKTxCjdtHeN2GctOTbvIvOsIbyptHekXlXcVaCcPMJkHnj17NuWnH/iBH7CmaXQroM3SM++lvDiD5uHLyxVXgz7IUHxbuDnx0riuMJ3OVVZIBh84aUPKVqjNuITpK0+3pe0ZhOR35cmBWcc6B3KMEjxvmM+LKi7rCpcXynXNzXPZuMlV9jJzUwjCTe2BcJNAIGgz9Ar6+fn5DEf7eCpqxcosyDbp+Lx7tonMNcHZJinfZE9lMMunygUtxzWx0rhUHpusMfm47leZVKpMyJxJPTZ/rudF5xtjrHHyj8krNQ9fvm1RYnT7roPiUoWbbHkB9XKTiTZwkwuryk0xXBMbv05uqopl4iZAjCof6ja0BQIBH6enpxiPx9Pr8Xgc5KnoFasqhkRMXqG4Nu8sNap0PO7kSPOKkTnGeOKU7yvHZy1XmZjqWoIOWfOufDlp6LNJKceVVyyo98Lsf22bCE9PT1decQkp3kAzCmUMN9F4wk288DqfWxXOCMWti5tiy6ZYJm4So0ogEKwSkrcCxmx7iQmPlYHKEztp+Lbp1CFrSj51rhrZ0qeE6XBuXULGElXCUus1rzSufGxKdNsUhLOzMzx9+rR1cjUB39gVbrLnFZumKW4KlVtHHIB/XquJFe0m0rgg3CQQCATVYM65XJ5K2gpoK9jmqYvZjuBKk7IKljI50XS+PHzbdOqQ07byZouTkn9V77Vv+5K+5myJMdPHbuOp8ozNfFIQs7q2SOh2Go/HU8WlbZ7qRUC4qTk5OfL58q9rZY3LTZz0sdxUFXVxUxuNKI115abYFeUmy5xH2aHy1+GZC3hYdF+wOUJtYVxEG1YuI8pGjr57dXlK65j4XJ7hmIed0vi+OqXIYKIuZcSVJ6ctOfKn1DF2e44pl+mNTjHSOPVrgiRix4xWXNZpi42LmwC7USTcVE0m4aZqMPua2feEm9YLHOMntZ04fbxqHr68q5RfdWfOopX1VUPV9mzr86CGFD1TZb60AuAdIWn8A8G25TNzMrGF2e67yMWl0LgUEhdc3s1UpcSWFyeuju97eDFe1xiECJybp5abowjFrpTFwte+qflS5T21n8QgNn/tDQbC/WldsShucsnignBTPQqiziekdMWOF+Em4SaBQCDQiDas6pxEYvMLKSp1TiChvGIUIVNZ45ZddUKsY8tRap6xBlhToF7susqzeS9cKyXzAK3XeDye+dYCjUv/VgXLwk1V+8e6clMdebbFk20aT6vMTRTryk0CgWB9kPRWQJ/CYIZX9V768rbJ5FIafPKH7seiTg9h3RPLvCYq7van2LZK2UJUR51jPMhNKi6u7U5mmabism7KiXCTHyn51GF0ctKsIzfVgaa3hKbI4ZJpnbmpChZtDAsEq4yyfPmdKr3tr6oDMfnlFZwCbcqEz7PbpsmVA1d9XPn7ts+42qWq1zwmv9A2mbq9vaG8ONt2QjDTp+ZlSxf7jOuAy2jQ7XR+fo7nz58HZVr1iVq4qT5uis3fhVXjpth4IdTBTaG+Kty0WNTppIgpUwxYgcCNsixxdnY2vaZnqoD4MVr5jJVrj77P++rKx4aqHtaQDFxPclXodqqal8/TF8qf0znmQcLcrSl1bzdy9dVQHuazC235WgS04vLixYsZ5crVV9Zlu41wEw/CTfYymhrXMQZvKI+6nl1TEG4SCATrhijDyjex28KoJ81MY9u2w50cQhO0TTlxyXB+fh4sL9YzG7PdyEzD3Z7kUwA5nuHUvDmYxypUFaRsneF40xcFpdTUG3x+fm7tA+ugqKwiN6WsvIXixhpn68RNupxlGCc2bmqj3MJNAoFg3ZC0FTBme4xPOeWSqa082xYv/d/mdYzxElJZYyZbV319ctS13Yjm45LBhzq886EtO5zyqqIO7zutxyKMQrNMW53Oz8+ny9hmHJfSPg+P/CKxCG6yrSy57rWRm3xlrhM3ccurinXhprIshZsYCNWVs9PC99xTnD0x8qXIpPOtIrfAjipt01S7Lvp50b5kLqjQxRWlFLLs5ZpTSj9M+o5VKMxFnK40MdtAbAYVR9FwGWf0Xsira1OYbPebRsh4acILWKfHlxrDdcGluMWCKgymV3ieCozPUNeKi0kMdOz5lJU2e7pTsChussXnGlw22Vz3hJvcEG6aLzeZ5Qs3hRFaca1iOOn8U9LGOjhi+bBKek7ebe0fbZWrSSxLnctydqeabWymGvwajXzHiuvhMInZvOe6D1xWNLjgTua2sjnlcR9MjPeHQ7gxCNWfU55vcuC0MZ0w6x6MvgksJS9dp5BXbx4w+2ZZljg9PZ0SBG1727NoQqFdNtTFTba4y8RNvnpwINwUj1XiJloP4SaBoDpkDKwGKhlWnImOo5i4QD1crvxcJG+TJUahsOXrm7xdZXLLc8VLmYxj5DTT1K0w2crgImZ1gaZLSetbTbChrnZyKeRUIdTPRnuDaTz625aH696qoWluouUsGzf5ZOFg3tykw9eVm1zlN40QN5nXwk0CQTU0yW+C+SH5desaPm9ZaDuILR/bpGOLG/IOxywV2+LYJnHfxG6T31V2HZMG1/tqItQmIU+tr06+MI5sTcWN3bpgPuOUNk6BbWXWJo/+TxUXM67tty1u3XVoI1z1jO3Hy8BNvvzmzU0huNoktT1C6VeBm3T8eSteIW6iYcJN8eD0zyrtso5tuqyQ59MMXDymETpTFTuGkj4QzAnzeWGrEIMrvOrqECedz0sbu2LlMtpiPcYxcHkCOR7COr2IHNlTlA6dd+rqQ53ypOTJUcL127X0PZsx7PMSu+6tAmJWF1eNm8yyKObNTSnwGXrCTfXKk5KncFN7UKX/yIpIuyHPZ/EwV80pV9nCbEjaCuibHDgKRMrkEhM/tIpGy0/ZzmErs4oH0ydDHd5OXUbo+VRRNrntWHd+NH7q5ByzWlUHaD+l1+ZfWZYzH66zKS603i5isIWvCpaZm2zlCzf56193/22Km8z8hZvWk5sEghCkz68GKp+x4k6GttUejhIU453VJG92TpviFJI1VA4Xvskhxgip04sRksc3sOeh5NDJNWYFz5zom0aVMmzKhk2RMccAVVp8iktISVkHpWUZuYkja6gcLpaRm3yYx5ivwk2xaRYF4SaBQCCohspnrLjgKOy2iYcqIpSEbfnY8nd5cmkcW1yuF5gaJz5jJdYzalPKXPn6Vp1ivddmmlA6X962SdqVBye/1Pw5+bja3FZeCkIrEq5+Y/ZPW3/1/aZ1q1qHVUIsN1HFUccTbvLXKcRNKVgnbqKyNDF+Xf2Ooglusl2vEupwxqXkVSUtTV+lXJo+5Gjg9AFX+lXsP0C4XlXqvcxtFtMutvHgW3Axz19xkfwdK9tg892ra7XDR8q2OD4ZXTKHFAXqsfY9NF89YuKH8uIoFnXl5UJIQY3NN1ZxiUnDyWMeROOatEJE4LtHFSKXUd1273ks5sVNPmNIuCkeXGN/3bmpjrxSy3X1vzq5yXYtmMDncI5Ju8zKs0BQF6rMBxxuqv07VjbLz2cNxuTlK4ObNjWeDdRTWsX7G1OOr4zU+rg84jR/jgy+8DYRu27XeRobtB30+DDb3vzvuhcT1yzbd73qEG66HKcqYvqQcBMP5hgXbhIIBILlQvJWQJ/nNeTVMuO6PFS2+7FEG/Ke+e6l5GWTlWP92jyDNI6ZZ+i+SxFZhDeYE94mNCkr9fxSmPfyPHfGc4GzBYfKsYoQbqqfm2Lj+Awh4SY/qhiRqRBuWi5UWcESCFYJIcfYIvik0hkr7pYRXzzuxMxVHFKUFZ9HNBWcbUG2NLGIaVNX+lC5Ifnr2DpkyhGTT0qZuqx5b5GwtTNVMpRSyLLM2SepgmJTWGhc3+9VhXCTGyncxM2vzrjCTen5xEK4qXm4+rKrn8f0O25/EYNMsIrI89w7B9D52BwHZTl5q6kGzccMo+GuMmvbCkgHbMjrFbpvI2Jbg6TIFhseGy8WnHzr8GDGPBOzDP0XilMFdRmVTZUVi5BHXreZOXEqpZDn+ZQkaNvTe7Y/sywXAdDnteqTq3BTOtrMTboc4aY4zJObzDCzLOGmMNqycmfjN4FgWcDhEJ+em8JBtWwFNLcQhDy1MZ4VTiVpGRwPj82LFso3JBenDq48uWW78q/Ly5eaTx0ea/r8mkad7WYDHR/ms3LJ0e12k8pI8RCbZa8alpmbbPkLN9nBybeO8oWb0rmJXq87NwkEgvpBjaJFc0fyWwH177omNzphhxTFKmXa6lDFqNDgWsYccJUm2/0UpdCVjtM+HCNW/w7Vv0kPJZ2sm/aG+lY7zPbMsgydTmf6Wk/OKkbMuKOKi+33KmAVuMkWT7jJHT+Ubww3hSDcJNzUBOqua2w/iolfhxOVi1XuA4us1zK2aciR6dMtqePINVeZeVU1zGr/jpXPm0rvu+6ZFaMVrMNbS/OKzbvtAz7GCOIaTSGFyBWnTgU3FbqOi/BiUE+sra0Gg0GyVzy0aiWKy0usIjfZ0rX5uTbBTaHyXHGFm4Sb5olFe9FdsPGeLy63HrEKakjhDcXlpIvJY1Fpm3Q+t8GBlJrOTFsUBbIsm3KFPvsJTJ69eS7Kpr/7nLC079AzVpy+VdmwSvGO1PHgfF5kbrnmpOZLx80vxfvNTedC1bY0J9WUvOoahLHe8hjCnocn2KU4mM+byqyUQr/fZ79lK6R8+FY4bF6btk70dWFduInDQ4vgpqoQbqoHlJNo2cJN7cc8+skyYVWN70VjWdpUy2n7eK9tLo3Nuyr/VPpAMNcgiVVqXORL8wp5WkJeXlua2Ek0Nq1rEkuF6xm0hXjq9pBUNUK5nv9Y2JRAsz9TuZVSGA6H6HQ6Xo81VyH29YGQh3hVMC9uoumEm+LQlr7XFDelGoHz5CYzzMVNRVHMpLflK9wkEAjmjSzLpk4foH2rwZVWrGyKRArZ2ojfpvxwjC2f0sSZuFyKhU1mrjwcJSvFwxwCNw+fhe7Lw6fAhvpDSPkNxU/xQqSkiy3DZ2CZ8UajEbrdbpAQbF7mUP+1KSwuxaVthFQXUp0NHG4y78dyk89r75KbxnHl67oOyWPLP5RXVb5qAzeFykzlmFjMm5tsZZvXJjf5+CGWm2i8deMmV30WbUTKapig7dBb/XT/zLJsZrXK90p0rs7kuvbNQS5Ueisgt2Cul1ff93lxUwY+R/GwISQbRxmqk6hiDZGqefjy4YZzyuR6audlaIbg6+8hpUwphc3NTZZRFUKKl9dUZrQ8qwabouBTzs14rnv6ftVxReVI5aZQ+fPmphi0hZtijNmm23EduCnGKaTlWQeEDF6OA8Q3pnzxQnm7HDShcgWCOlEUBTqdzrS/6TNWGmdnZ8E53ry23XOFmwadzeiyysupFC0w5A11ESjX4+pKl+pJdOXJ8fKa3kRX/U2E2ibF+q0DPgUzpDCa6WPT+dL6wl15xhjErlWFqn2IKgDm/ZBc29vbUYpLrLeFyuf7z5F5WbDO3MSVdRHcxGkTX7nrxk1VIdwkEAhWDXqFShs5bTfok1asOJOIb9ILeV845Zt5+ZQl10TmytdljXK9aKF4OrwOxciXzqWcNekRpt7GkHeUY6hy8qJowuNJJ38a5pJDx8+yDDs7O+j3+9N6m390Kdsssw5Fx6W4LANJxWBZuYkjc8hT1hQ3ccZTTB/iGI4htImbOPJoNMFNNoTqyOUmGy+Z+dsMuBRDdx24SSAQpCHPc3S73Rl+aLPjpZYPBJv3gDiPL41v5hOTl1l2jOc1dpJPqQ8nPMY4SoUtfZ35cZQ/M25M+aly2ibtlLxsSifnOWrFZW9vb6q42JBlGc7Pz6eyuhQaVxkhz7ArjzaTUyqWkZvqGP9NcZPtvnBTdVRtMxt8dUzhJv1KYy43UdiMpXXkJqoIcp0gOm4orautfPzFlYMzVjhl0TguB06sbE2Mo2VGW9siZu7RL6fQc2a3251eA5MzVWZ/0TylQfulyVn0WuteHLm4vFTp5RW2Ccg36E2CtV270sfIEEPGrrjcCcCXb1MTd1UFxKVwpCjnrvw5ctXZnj7URTK+NnMpAkpNvq+wv7+PwWAQ9Aab5ECJIoSY57gOExJ3jHG5qYoxLtyUln5VuamKg8eXnylPG7nJF7ZO3ATEr4b7DJDYuLF5cQ0eTrxY45LL37FGmaAZpD5TbUBpdDod9Pv9qYFVFMXMWwDH4/ElQ8o3Z/reIGi79o2LUDiQ+Lp1l2ci5FmxKRohr5qtfJp3SB4bkYTu0bCqoHLUOfC5kxCXEEPxQ+G+58vJl8ZNeQZVFDyaxjaQzLqZ/6kX5fDwEMPh0KqwVOkDMcoKN2zZEeIm27UvHpcjfHGEm+bPTZywRXJTHTDLFm4SCATLgBAvaOitf+bK1TKh0hkrl4fRvAbCygjNK8ZQ4xBxClH76haqN+eeS4lKNWY4cVzPI7ZMn0LBUdC4z62Oyb3qJJ3qCVNKYWdnB5ubmzOKi7m1Rv+2eYi5soX6aKgdFqUYNgXhpsVxEwfz4qaQ0bQK3JQqk3CTQCBYBLgOLW1Y5XmOPM9xdnaWxEOLQq1nrPR9IM6z6yLQGGUoVlngTGo2zx8nLuc+Nx+OomTLgyNzKkIeyVDZXKUiVfEwld7UNkj1Ris12WJz5coVjEYjAJhRVjQ50H298/DIcBWaZYdwU7j8FKRwUVOypOQv3CTctEhwnJGctNz53xU3Ni8bd3HSr+IzFFSDqz/keT6z9a/f76PX60Gpl1sEx+PxTBqzH9IzV/S3K0yn9fVT29bnUL+utBWQO3Col9fnXbUJbVMifMqKq4xQ+dw6uGS15dfEBM31QtswD6KrQw6Od9kMj+34Ibg8rKG4WZbh2rVr2Nramg520wMMwKrAmN7iVDl9cq/DBBdadXGl0fHmwU2+e/PmJk6+TXAT12vZBFaRm7hxhZvaCdszDDlXbfFceS/KkcFNr7Hu/WCdkWUZ+v0+iqJAt9tFp9NBp9OZhtuMI/Pa9/IKG3eZ38Ci8WPg6rOVX7dOlc5Y5SIkKK2wSxFyxQ+FhZSxkFJQxUPqmnRDbclRShZBUqme1DrKbGKS1vlRxdkMM+9nWYbXXnsNW1tb03agCov2Gts8xTqPqkveHKdDyFhcVvg4oS3cFGMwtZGbfGUIN82WaZbdJDfZwkwZhJsEAkFbURQFNjY2ph8CBtKNnTag0lZACpvXV9830/k8wdwJKNajW5dCwikLCO/hN+PGpI+VpW4lxpcfZyDUKU/Tk3CMApllGQ4ODrC3tzezLK0VFPM39Qo3qcDYlP1V9xYLN6XJxeWmtqJt3BRTdtX8hZsEAsGyIssy9Hq96dmq8XiM09PTRYuVjEqvW3eBM4ED/gPfGq44tniUqM30sURtW5Xj1i/kDTdlS/Hmco2vqhNTrNc9ZOBwJszY52VbTUqFz1Dz5Z/nOW7duoXd3d2p0kGVFu0N5iowsX2Vm6ZORX0Z0RZuqvIMqnCTLy963RQ31dHnYlcdfdzE5ZzUuUS4SbgJiDOwfTxD49n0iNRyYxEaV6GyfQ4Ijj7RRjQpV1Wdb9Gy6Th5nmNjY2PmPNVwOJxykS0vc/ue2TfK8uVnIvS1CTNMX4e+E2rGt53fCvXr5DNWNu8mVSxs3k+OwhIK5670xEy8NH2VtD642iY0aacoTDaZbN5NW94xhiM3jBNuk9EXLyZfrlwxbZ3nOe7cuYP9/f1Lg9U2wIHL3mFbn6t7282qKy3AcnGTKy9XmJm+aW5ypa2Lm1xp6uAmH3xtFFJOXTL64vnKjEEVbrp79+50pcomH70WbmoOIcOem9Y1Bl3hPkPcJpMvfZW8UzjCl3YV+4mL49uCWLlc+rh5popuAQRmnT86vgnKUbQcyk++OSTE/THzrUbyGStXAb6Jmd6z5ROacF3l2dK7FC2qQNjypmXaJqJYQ8emuPjkNhG6tikltnqkyM6FK9+6y62an/kcfPKGkOc57t27h4ODA6uiob3BSr30CJu/Xf81YhWYkDJHx8MqTkyx3ATM9id9PQ9uomkoR7gmpdBzWyVuqqt/CjfNQrhJIBA0Ddc4LooC+/v7KIoCo9EIZVni7OxsARI2g8pbATkE6DNwbPn4JlfuSgyNn1KHkMc6Jr/QJBmjsNkQUupC90KyhsJCil5Mvtz+lDrx+ozkmDwPDg5wcnLCUjBs22/M8vT/ul917BpnqX17mZCqnM2Lm5oK82EZuYnLDcJNLyHcJBAIFg3XXKOUQrfbRVEU6PV6ODs7q2xYUaflIpG0FVDD5X2lnl+fB4rjGbUtKfsmxJBnmV6HvMi0Lr48uTL5vNm+esSAm5aj+ITy88lap+exCS+m75naUBQF3nzzTRwcHCDPc3Y5pldYX9s8wmb5KdtuQv10VbGK3ETrJtw0f27i1rVpbuKAw01lWV7aWiPctDhwFEHb+OSGu8aQLSyUPpRWsN5w9Yksy3B4eIiiKLCzs4NOp4Pt7e0p5+R5PrMV0OQiDfM7VpTDyvLyGSvzv9lvbc5R85qeqTLz4qLSihXHyxia7H0TZwrpUm+XzyMWqgenDj45aNpUcDzh3PIWRYZNGF4p/cPWN0IKtomiKPDOO+9MvcGhQ5B60FKlRV+bHuI6DovTsm3XMfVdVgg3+eWgaVMh3BSfpwtNcJNLIRiPx1HcRDnKlDkW68pNTRolizB4mnC4CJYXvueulJp+n2o0GiHPcxRFMRPuc9rQ8FBZNE5Mn+TMUyFDK8qwivGEmNc+jwfXe8IZwDFWpc1qdYX5yKLKpByqU0x6V9ulkFyddbV18ph8fUhRXPR/W7/j5Pf666/j1VdfnRn0Nu+JzeNi8wqbSoxLlrqUO5cCswoQbgrnY6IqN/nA4THhJrcsdXOTiRA3UaMqtGol3CTgQoyu9YJpDBVFgePjY/R6vem5z6dPny5YwuZQ+4pV1XghsjXzcylTNs+y/m0Ls01kXFls15z6uNolxuvsS+eamJpUejhyUnDLqTrRVq2PfsPW3bt3vdv/qJFlu29TYnyHxWPHjS/MFncVJ7tl4iabLG3nppg+I9zkT5tarkad3AS4X8Eu3LScmLdBw+VUM24ovo/HbGmXyYhb1T6vX5uulEJRFNMVq6IoZlbDgTRHlKvvxDpmqjjOXGmTz1jZJmHzvi+eSwmwlWNTWmK9liFvdgrxcyeM2Mk7pOzEdILUARurUFRRTpokQDq4UhXDTqeD+/fv480334RSCuPxeGqkZ1l26b8J00NMDSrAf1jcJautP/vGBle5XnYsKzf5JhjhJn86X7xl4qZUCDctB0ynrr7W8CmJrvgp5dvSc5wZ3HxDsPV5F1eHwmP6nS2+K40PVfitSl5tyjuU1uSFq1evotvt4uDgAJ1OB0dHR8iybGpwdbvdaVrNW+Y1lcXGX+Zvl1OTXlOjzvYpCt/nKehYtiF5xcrleXVNuDGeU52vL04VJcKnEFVVGEII5RdTfmy7VKkLh5hi8uLKEpqYm07f7Xbxzd/8zbh+/TrOz89nFBcNWzuYSot5z7dq5VJaYl51zPHg2fr/KkG4KQ2x+Qk31cNNqfXudrv4lm/5Frz22mvCTUuM0LiOHfdVjLB5G7RVDcaYchZlcDSZ17wRmjNN9Ho9dLtdDAYDdDod5HnuNOJNo4wTbjPSbem44aHxFTvvJn8g2Ly2xTGFMS1J86GY/6nwnHJDyohPsUrxovkQo2RwwzgyxJCuTwGhE2ZVcJ9l00ReZXJWSuHTn/40bt++PfWgaMXF7NPmb+od1vWjSorrni7X/E9/265t9eYqMKuCZeKmmHi+sjioi5vodUofdEG4KQ5ZluFX/+pfjVu3bk3fliXcJBAImgSXi5VSuH79OjY3N3F0dAQA+PDDD9fGWVLpA8E2BUTf1/F8kzaNw1UufHLQsJActrxC8bj140zgVTxPNqvdpwi5ZOLE4aTlgquc1aFQVUl748YN3L59e+bcQp7nl/qW/qNbb7RyYsal93zeYddZBq7iYevTq67ALAs3xTyDtnCTTRl2yWXjJp8Mwk1xuH79Om7dujXlJrPPaqRwk3lfuGlxCM39ofB5Yd2fk+AliqJAlmUoigJ5nk9XrPTWP8oZlGtjnDW++Jw+GeMgDM1lNtTy8grbfZcS4vJW6TiUMGIVGl9am5ITI1NoUg/V1SZDLDH50oY6IyeNC3UpLjYZmkDV/G/duoVv/dZvRVEU0y02WrHQyow+yE09w/S/T2GxKS51vyUwpNivImzGEb1eBm7y1a0ubnIZgVzjx1WecJMdVfO/ffs2Pv/5z89wk+ahqtyk0wo31Q9ufUKGFI0T004uHnNdx8rm44tYPhEsF/b29tDv9/Haa6+h1+tNz1gBk2fd6XSmfYV+qubs7GzmO1UhaM7Tv02Ht8l5+tr3plTbvdC1hqv/1vbyitAA4SgDrsmfowCZ6Tl1oOQdo9zY5NHXPiuaI48NXPKJsdRdBFqV6HzEXEVRii07lbS73a5VcbENXuoFtt2zbbvxGVVUWalDebG1xypOajHGjStdG7nJlK0JbuLKY8OqcBM3/0Vz0507d/C5z31ubtxkPhPhpvmgzfVvs2yCxUEpheFwiI2NDWxubqLf73vfUMrJr4m43PyqOnVqeXkF59qWLlYR8nlZNULlheQLye6S1aZYcb3JPjnrROrkFSOLK9+mjCraXlXb7du+7dvw8Y9/HOPxGGdnZzOKyPn5+VSB0QaRVlrMa9e2G5cCo+Wuorxwle51mBTbyE0hnmk7N4UMx6poOzel1LlubvrCF76At99+28lNNr6pyk0up4+rXsJNqw15RgIN2hfu3buHk5MTvPrqq+h2u/ja17526W1/64JKWwGBakqMa7LmhvkGechzHMo/VrY6Joa606d4AmmcqoZOTNmxRmWdJN/r9fDGG2/grbfemhpP+r9LKbH92RQXU0mx3Yt5C5erbW1v4wopObSvrxraxE0crmgzN1VFCh9QLJKbYuKZZdUBzU0f+9jHkjipLm4y+7ZtS6CvbWlbrCM3ufqDjwtccTjtqeNRjvCVH+rjMXKF6sUZh674beC0NqLpNjHz73a7yPMco9EInU4Hm5ubGAwG0/Fu8ofJO7a89GvYNWzfuTLDqb5D49J28F3T+KF5htPGyVsBaQGuAcyd3ENKAs2LouokTfMxy4wdxL62qKrwcBU2br5cIrV1rlAbhu7ZwrjGlylbSFYOPv/5z+NX/spfibIsMR6PZz6GSRUS2xkp259JJLqt6LcYXIqLy7Byta1vUgr1mVWBcFN83pw+wkUVbkoxeHzcxJUxFN8Mr8pNqUrPt37rt+KbvumbZrjJZVhxucnkIi43hZw+MW27btwEuOtj6xt1GiwuBZQjp08ZpnFtiCnbHD8c+Iy8kJHGkScGVfKKVdzrLCsmPjVAdnd3sbW1hbfeegsHBwc4ODhAt9tFWZY4OztDnuczDh96xurs7Gwm76JwmyP0eVLDqyiKGe7Vf/rad8bKpdvS8mMQZViFFFkbWVIlQMO8thGpLdwmj0+ZDHmIOEqG6XGzpXNNmC6rl6s8VFV4bOlDebny5d73Kf0+uWMMIlsf4So+Lty7dw/vvPPONL1WXEzjiv43By/943qEbX+6HlzlhdMfqih1y4JFc5PP+PKVZyvblZ9N/jq5KRTPJ4svPid93dzEicdV5GLrZeZdlZvu37+Pj3/847VzE3croHCTQCCgUEphY2MDe3t72NjYwHA4vOSYWWdU/o6VvgeEPcA2bxSdgFxlmggpMDaZbGlp3jb5OUqHqeC40rjiuOriq2MdE3aVtC4lpUo5sbJw+g0H9+7dw3d913dhOBzO5G0aRK6tek38meWb7cLpFz5wDI1lR8h40td1c1OMceXLX7jpJYSbJkbVd37ndwo3CQSChUNzjMbVq1dx69YtXL16FaPRCOPx2Plx8HVDpbcC+iZVF1lyDSGuh5ZLxjEKhuteSlod3wyrqrC4UHUispURuq4DMXmGnmFMPq+//volo0qHUUWFHhJ3vTJde4bpvSaUlzqeBcdjv4yYJzeFDCXhJuGm2DLv3bt3yajSYcJNywfd9hqh9gk5k11hNNzHGbZrmxxVx35qWYLFgj6LbreLfr+Pra0tbGxs4OjoCFtbW9MtfJRnbOegzPnK3GZKz1jpb1+Z12ZedHWMck/MXOGbAznXNkQbVjHKiDnR0LQ+hcCVpy1OTHqbfL5wW/mcckJKkCt+VWKh6WPbIFSXOkmPPn9u/DrkUGpiVNkUFzOOT7EIeYdN0ohRZMz4+rfrGfj6E60LDVvlCczHTdQI0nD1fRom3JSGWG6iCMVfJW5yGVVmnLZxE/f5rjs3aYT6lq1duOEhHqC8FzLUKHzpY7kulLdgvvC1+WAwwP7+Pu7evYtXXnkFBwcH2NjYAIAZ7jCNJQ3zjcoa9APl1LAKXZugZwHp+XZaFu3zvmtbfA1Xe9W+YqXDuIPTN4nZBrCrPJshp+OYZYQUqNAE4ZMhBBeBucjI5x3iKDwx3qWQsuhKb2tzV34h+Vxy1em5VOqlUTUajZzlUUXCVFhs96k3OEZhcSkvZjpbe7kmTM5zW0XvIJebOHmYcYWbwvLWzU2h/rmq3HTv3j38pt/0m6xGVZu5ydb2JtadmwSCZYQ5Dnd2dnDr1i0cHh5ie3sbnU5ngZK1G5W/Y2WbWGyvSnQpCCGFwrymE5mPfG2ead9vW1kuWX3pfPm5ZOZOJC6litOevvRcGTj5mtexypEtfJ5GFY1re2Y2JYOGcRQT1x9VdMzyfc+chlPEGtXLDB832QwI4aZ2cxMHq8BN2qhqmpvotXCTQCCgoBxycHCAe/fuYXNzM8hR647aXl5BidTlXXWFc8pzeTW5io0rL5tSYubl8xr75HbJYZOnClLyqKpYcPKpq251KDAxRpWZhqt0pCgqnLT6t/k/tHLhAu2TdSqGbQCXmwB7Wwo32eWpgmXjphh56+QmrlFlpqnKM8JN80VMnbhOmrrLrTPvlPFh41Bf2TS+q+y6eDEksyt9yphIRV35dzodjEYjbG9v48qVK3jttdewubnJXqlSavbD5LQNQlv98jyfXuvPTJjhJureCkjPc4XiUyStWLkUFNc9n5EVGnwupcJFwD7FyrwOxeOA6401wznl2OJwCYZbB1/+IWKogpR86ij7+vXrUUaVr3yXUmNe284xVFFeXM+4LgV4VRDLTSZchpELi+SmGKWjaW7ixm+Km+rCorjpxo0b3jNVMeVzuMl3T7ipOaQaODYjIdVIiy3LFu7KO4YvOeEufjXlqNMArRNNyRLKtyqHmb+LosDGxgZeeeUVvP322+yVKhsn2Haxmdecl1e4DCttOFF7w+wjPsPKZp/Q+LR+oXZO3gpIhbJ5finJ6ns+ZcSVP/3NCXOVwR2MtnSxSk9VpSZG6aGgnSXFEx/KP5a0Y5DqQaJl3r17F5/97GejjapQ3VxEkaq06D+TGGzKdmx/dNWrTRNRE3Bxk0850Gg7N4XK4Nyri5tSUCc32dItCze9/vrr+MxnPhNtVFXhpip/dXCT71mvCzcJBG2Ab5zt7OzgnXfewd7eHnZ3d+VMVQSSX14BYIZgdZhLgdD/bRO7eT9GWeQSOs0npADFeINd931lpXhxQ3XzyaDhqncKbPLXNRmmGpK2PK5du4bv/u7vrqy4hAynUBz628yX5mXW33zDl5bL13d8CuWqKytcbrIZL6vOTdyyhJvceVZxctnyeuWVV6yfewhh2bnJxlGrzk0CQRvh47LBYIBXX30Vo9FIzlRFotKKlU2BC02YNB69Z8I38VPlKHTfV4ZLgTHDXfKE6h+jCFWdsKmMrrxCZbgm2CZRd5llWeLevXv4db/u10UpLlxPMG1n27VNcXEpKba0Wh6bUkR/m3FDdauzv7UVwk3hviHcxIPNIKmK+/fv49f+2l+7UG6yxRFuag6pDq8QT9A8aDiHu3xxff0ulFddz3AV+8Oi4GrHTqeDvb09bGxs4O7du9jZ2cH+/v7MWacQfH2Pw0vmtbkdz3Zt8k7oTBW9jjlDVZZzet26WaHQhG1rWO6E7jNebPm5BjZXZjNPW11sdeOU2zZS8CmYqfDlsYi637t3L8kbrOFTqjU4Sksonk2Jofds/aiuSaxtfbMKhJv83MThrkVj3txUVxkx0Nw0GAyS0tfFTfSejZeEm5pHyFCicXy84wr3xfWlnadxtC7Pe5HwtW+WZRgMBtjd3cXNmzcxHA6TOQoIO4J8stl0KHodG9+UK2TUcfNyIfnlFXpAupQT817I0+rz1NJ4tIF8eVBZfXWh98w8bGlscpqTjy89LZejAPrCfXLawqqQV2zappUa28RwdHRU2agy/7hy0P8cg8qlxFCFx/xfVXEJ9e9lRoibdBx9b924ydYevnKFm9Jh46bj4+PKRhXlphBH2Z6FcJNAIDDH2XA4xMc//nGMRiMcHh5eevOegI8ow0qTOVUqzHBOHhquSd+M61JCQmlSlQOu541D/CEZUhCqi4mU8kJpYpS0phQXm8Kice/ePXz+859PUlyoosLxDIfk1P/pn+s+DQsZ7DH9YZXB5aaQQaEh3BQP4SY/N92/fx+f+9zn0O/3o/O1cVMqP4U4yHWfGlfCTfOBtNllSJvUB92ORVFga2sLo9EI3W63de3bNnl8qOWtgBo2EjV/Uy9ZyGPrU0I4ioHP28zNo2q8kCfb1SYu2CawUHpOfrZ4tmfD6dycZ5oCKivNR3+nKnWlSoN6hbnKjOvZuBRLW3qqoPv6k09p9Y1LU/5lIqtY2NrAxjnrxk2ue01zky9NKGzZuUl/p6pObtLXZpjtN5XPvBfLTba+IdwURlnOnuswEeIHn16lw23j2AVfu7p4wpfWVXbseLLxsi+uDz7HEye/UFxf/FBedXFi1fRZluHg4ACDwQBvvfUWhsMhbty4gTzPk8be2dnZJV1JX9v6vk9Psr1e3ewPVEYabpZHxx7n+dC+F9uvK33HykUIroJ9io5NYBdh28r2kXrsBMBVQjgKk6vcUHkmYjp5TP6xSlKMbKFJNWXg0r5l5vH666/XcqbKpbjYnquvjWyDNfRnKix0THCVVt/EyO2vywzhJh43+coNlWeiqvK7Dtx07969Sg4fUzmpyk0upbBOborpS+vETS5UrSflplD/DRksMXnR/FzPl5N3jDFkjgOunE2h7rKr5BerK3a7XQwGA+zt7WE4HKLf71cq36VHmWEuWX3XIR6z9RtX37OlpdchQyo0B9W6YkXvuQaSa0CFjBIzra2iLiXHRt42BSakRHBlo+G+PF2Kli+OS34bUurDkauK8sGB65nb8rh7924tiov2bJikwFVkTLnNfs6pr0uBsfUf2/PhTjBVFMdlxDpxk0spEW6Kzz+EGG6quopuU1RiuIn28aa4yWZg2coP1XNduEkgWDQ6nQ4+9rGPYTQa4caNGyt/pso1jzeBSitWrvs2BcKWzuedcE2kscaPTUZfXlSRoUpQFQMsVv6Q0pKahwshJSomPb1fpUP7lMB+v49Pf/rTGA6HSX2DY1TZFBj62wdTObH9mfWzKS0hJTulzquqwKwrN7lkrVomJ10dZYfQFDfVBVv+g8EAn/70pzEYDBbGTaEy6+ImrtHLqfOqctMyQp5HGtraZlquLMtQFAUGgwGGwyE6nU5Q5pC+My+jpQq4umhVZ12lDwS7FD3uhB/zIFxeXF0OvUfL5+TlUxRcdfSVYwNHpljY0sXm5ZoM68ibkyalnNFohC996Uu4ceNG0ENPQfux/js/P58qMXpfLlVgOHAZgrY4VG6bEe/yCHPqaoMvv2WFcNNl2YWb4mTjxgul3djYwJe+9CVcv359bbkptd+sIjcB8YqnTwl0jXdOub5nZItLucYVZis3VFadzzeV99qEOuR25aGUwt7eHjqdDu7cuYPBYID79++jKIrph75jVq00H2mMx+OZZ0CvKXxhNh3H9R0rwH7GynQ40bQhWWwczJUdqGhYhQYWTRMiiZBHNuRd5SoYIZljyMJVrvngbWEhZS5l0uekscmU2ha+9Jy+YZs4uGR7584dfPGLX8T29vZU0ciyzDoBuMrS8WzeYHqeQceNVWRMBYX+0edtPhuXwmKb5FxEsayTSypWgZtCEG6alakKN4XSpXLT3bt38YUvfKEV3BRSXqpyk6vfx/ZTwWVwldJ5tKVtrmniOXKdD4I0dLtd9Ho9jEYjDIdDdLvdS8ZUlb7l6rOxTgVO2S4+5aY349E5hzM/+FDpjBVVMmihHMPGpQi4Bm1IeXANTJ/S5VKWaH60bK4SwyUL2iYupcBXpi8+Z6KjcWKVIzNeilycNrp58ya+/du/HaPRCGVZYjweTz0WVEHQsCkbLmXFVE44hlQqaVBFhvYVl1Jv5uGSh2NYrLKCI9wk3OQrsyluunXrFr74xS+2gptSeEnXg8tNLkNKuEkgaBeyLMNbb72FjY2N6UoVx4hwwZYuZOysC5JWrEJei5gJ1Wcdu0jYRvQpRGxTYFwThU++0D1fHj65OHVKmXw4ikOsohKTLqUcjRs3buCLX/wier3edBnanOSp4uJSgPV/+hc6yxBjcLnqTsmMKjAuJTY07mz3XYbFKmJVuYnKINzUTm66efMmvvCFL7SGm2LRJDf5HANmWYLmYXPkmNep+aXkxXFwtQEpXBKTX1PQK+WdTgdFUaDX66Hf70/v0zHo4w1bGHUKhXinLoPL5ejRZfj6IkcG3zzAMUZred06hzSpUCFvFY3rK9M22fsUi5DSEpKzyqCvSyHhKEgxSkLTMsXEcaXb3t7GZz/7WfR6vamiocPKspySRRXlxafA0PQxslOlhd6jv13/zTxtdQr1i7qIra1YRW7ylW2LI9zUXBxXuu3tbXzmM58RbrLUi9uvVp2bXN+xAsKKu89g0edjXHF93/AJwfdMQrxG++YiDaUYw4gq6pz4ofSx4amgeY9GI+R5jjfeeAODwQDXr19Ht9vFs2fPAEw+CmymzfPcK4fPkDo7O7t05sqMS69tspu/fd+xCoVTR1NT37FyPavKWwFtCCkmVHCfwkHvcxQfOrh9eXIUGo4CY05AHDl9oHnVCdcEGFN/Tkf0xU9Bt9vFl770JWxtbVm9wVmWTScYqri4JimbEmMqLrbvxlAy8dWVKizmfarE0Gfu+u8yFKqS8aqhjdzkS+uSQ7ip/dzU6/XWjpvMZ9IUN1XltXWCS9GTNlxvdLtd9Pt9DAYDjEajmZWq0Nxim0doeCqqOlFCvN40QmOq8nesbAh5KWIepu2+LU6swWULcym3VOGy5RNbHx9cE5qvLF8Zofiu35xrwF3Hujr/9evX8alPfQo7OzszioWpsOjyOMqLlln/D/3ReLZr2hauutPnavvtU6yrTJwuxWedsEhu4qT3yS3cJNw0D24yf7u4ydbf6+YmgUBQDa+99hp2dnZw//599Pt9b9x10wWaRPIZK5uy4CNuSpwcIg0pQWaZPtlCik+sIkTrwfUAxICj4PnqwZnAQ/ApNdw8Y8ql7QhMvMHf8A3fgMPDw5nlXFNxoR5hmzJgQ0iBoWGpoDK56hwypmifqKuvrQpWkZtCdbTBVWZdqJub6G8O2sJN3/iN34iDg4PWcFMKT9Fx4TKkhJvmB9s8wU0n7V5d32ki/3kgz3NkWYZut4s8zzEcDqff0bNtGTW30Cl1eduoa66kv23X5j36v27QetF7ofmIG1dfh+pR6YyVbRBroo+Z2KiSQhvGpRyldu6QwuEzRqoOKJfBwzGEXGGu65Di5fvty99VZswzsT1nDX1/d3cXn/3sZ7G3tzf1AmuC0OlciotLebORQUiJsZEDV5kJGVM2BcaU26bEmPkIZiHcFM7bhUVxU5XfvnqY101x0+7u7sK5ycdHIW6i5Qo3NQfXs4htq1SjS6f1xV82wyxWIeamdd3jxKWKOjffOtDv99Hr9XDt2jWMRiPcvHkTw+EQ4/EYT58+RafTccrh0lU0QoaU7TtWZlyOIeYCbVPNrxqmgagNRrMceh2ShXJo7NhtZCugFsZGwDavWCidyxsWQwK2SYA7GcQMEq4yYj7sWAMophxunjFtktpOLqXUlsfW1hY+85nPYG9vz2rQmNtt9KCyvdbY1vc0XAYTx7iy5UfraQsz6+5SYGz92zWGUrBsk2fdSOWmUF7mtXDT6nKTflEFl5s0muImmsaWH62nLSyGk4SbBILFwTUm+/0+NjY2sLm5OT1TJZg/ogwrGymHiNRF8Jz0PoJ1EXsIVZQlzqTLKctlTHLu2+L50oSUFfo71guZqiz5cO3aNXzyk5+cnlugHmFNFqYio5SaWeo228HmkTVBDSeX0mIzrkJ1p8q6LZwqKjHPTjDBPLgpxdgSbnKnieUmn5xt5yaz3i5usiGGm+j9UN1t3BQyooSbBILFw+Z4AYCDgwMcHR1NV6zMt/wJ5oekM1aAfQKm97leLRrHp4zY8vQpFfq+z1MXuudSGFyTjWvS5ygJMQpSqPzQ5OdTTkKKi28CDT1DH7rdLt58801sbW1ZlQr9G8DMOQYXbIooDaf/uV5h27UNZnuEFJeQR1jnVydWRRlaV27ypV1WbuL+5lzTsFRu6vV6ePPNN7G9vd0IN9nC6X8XN/kcCz7UyU2p7RqSbx1gc7q5wptsk9DYjh1robRVw31xaV8NxffVy5eWE+7LOzbcRJZlyPN8ugVwa2sLGxsb1jeRauePeT0ej1nc74LL+WwLs/13wfU8XbKFrk0epjsJbPL46sWRP3krYAqJutKYxG4qGiHSDimeVFmypXVdp8TjKA8uZcs1GG3pOPddedji1aG0xBCuCycnJ/jUpz6Fra0tAHZvLXB5YJhxzO2AZrhLHp/yQq9tsrjgUlJomL425fQZUlzFfp2xbtxkG980nnBTOJ4PV65cwSc/+cmpwwe4zBnmPY15cZOZpi5uqtPhI9zkhq1tfIZWTB/igvYBbhrb2PKN3xhlmSq/NL4tnMZ1cWMorc0h4ivbJ6cvn6rhnU4H3W4Xx8fH2Nvbw5UrV6b609nZmVduDi/G9Ct6xopykflijNh5jz4T+t0qesYqz/NL3Mutk1KzL/Gg8nLGSrRhZWaqC7SRgot8Q94ul4fBLJ+G+RQcW5oq/+vKK1U2lxyhtgkpK9y2tuURC9qHgMmhy3v37mFjY2M6IOkkTyd/DdPQsm0H9CFWSQp5ZLh1dykz+trMl1sfSi7rpsgINwk3ue5xQfuQUgq9Xg+vv/46i5soDywzN9lkdT2rUH1i4goEAh76/T62traws7ODnZ0ddLvdRYtUC5adH6INK/rKRk4DmF5fLslT+BQkMzxWOeDG8ykuvjx8MtYRTuWIiWOrQ2y4LSxkzduMqs9+9rPY3Ny8pKyY+bmUCuByv7SVaetD9LfNGxzjBbbd8/2nsqV6hCmoR2cd4PPy2sL1PeGmaoYcty7Lyk2f+cxnhJs8v3V+Pph8tI7cJBA0geFwiP39fRwcHGBvb29lXlax7PyQtGKl4ZsQbF4uV2OFJjzuROy7Z7tvu45RJFKUEW7eoftU/ir3XHFs8WgaWx8ITbJmeL/fxzvvvIPRaDS951ulqqIAu+DzDIfS2MBVZHQ+XGWFY1hxvMFUvhiDbVngawfhJuEmXzqNfr+PT3ziE8JNnt9mHhx5hJvqR0wbhXjFxyXca/M33YLFTQvEb9eLDY+Za2P6O827TphnqjY2NnBwcID9/X30er3p681txyQ0fKvmNrl915xV87qclTHPS6nLLzKLed26zTj1xbehka2AOp7535cPh6hdZXAUCtckHquAVFFKQsqJT0ZbOtc9XxxfPpxrTh4x6Ha7+JZv+RZsbm5O79G+YPOmumRMHcQuMrB5hX1k4VNazN8+YjPl8F276rHO3mDhJuEmXx4x4HATvb8u3ETbIMRN5kH6deWmRYGOU33PdU3jUyUzdB1KT8NCXODrg7HXIbmrXPvkrhNKKRRFgY2NDRweHmJ/fx97e3sAMPOGUo2QQWAzvMyyKHzj12ZYcdKZ5dmcnvo3Z97gXvteXuE6w++SxYaktwLSBqDgTHq2BvQpQfSaO4H7FABOfHqfo6i4rrnpbHnEKjG2vF1luMJ896qiKAp8/OMfx2g08ior9B6Vh3o3bQPY7HM+hDwsMUoBVVp0etdvVzxaBx84RhVto6YmgEVgHbnJ/N6UcFM96HQ6eOedd7CxseHlJt8YNu/p+8JN8VsWBYsBx5EnmD96vR729/exvb2Ng4MDDAaDRYtUK1alz0UZVnoy5HpxfROLbzKnZB5SOGxxuMpFqAwaHlrapnn7lDTuf+49rlJSl8LieqahyVErLleuXJlpI99/V5k+mWxxXYoSlTt1GdulgPgUF1e8kNyCl1gmbvL9Nq+FmxbHTScnJ9P4wk2X4/nu1QXhvMVB2r6dKIoCm5ub079VOVOlsSoGfdLr1n0ToUaVSdv2akSfsmGGh+K67rkUEdMjHJI1RQ5bPW152cLoPU4crvIZiu9LZ5t8NfRK1cnJibMOIfi25Jj5UeXHp/jELluH4FO66G+zXFe96iAbV/1Xgcg0hJuEm3zpfNzU6XSEmyK5qS7uWDduCt33jbVQOP3tc7jYtuOZ1/QDszSvkDMnNpyW5ePzqmeofM+D80x8qHNcFEWBTqeDnZ0dbG5u4uTkxPn2PzqGYq45jhRfvWK5yfa8XHMXfQEQXQGnr1c3+y11RNF62eroaydOP6j08gpOGGdStsXzKRm2PGxKSEh54Co1tvtcmX114bQTNzx0HcrLDAt1vFR0Oh28/fbbU2+wrewQfAqJjRRC21o45dk8yz7PNa2PTT6uktK0R3iVtt4INwk3pWKZuMlUIELc5IJw0/wRy0/cuD7+UurySxxijKOQ4eUzfpRSlwwzblrgslFHwznnvVzXMfxku+aGpcD2/Lrd7sxKlYmQsUTzdTlzYq9tSB2ztvYO9WtXei2nyxiyxfeduUrZIVCLYUU9wDaPsCutSwGh91y/OWE2z65L2QnJEXrYIcXF9ZvTaWKIIhTflcZ3PwW6L+R5PqO4+BSolPzNa52vTYEx47jyC8UJ5W+TxXfPVg8NffC7KlZRUaEQbpovN8VMdK6wmDi++ynQfaEoiqXiphB8nGSTxXfPVg8N4ab1A5djbNec9L7yuHLNG006GDqdDk5OTjAYDHB8fIxOp9NIOYJ6UYthxZ2oQ3G4yoRLMXEpMBwlxxduyulTUHxt4ZqsY9rMdR0TjzPpcqDTcxQBbVQdHx975bL9riprk6QXAldB8clY17NaB8VFuGn5uakOcPp7m7hpERwl3LQaWMTc5hvLMVyxKmiyXnmeYzAYTP9W7UzVqiL65RXmXkYN1wTBnah9E71NOdFlmUvJNmWFq8hwlKQYpYSjzITaSP/m7nMNDW7fc/Gl4ez7990rigJvvfXWJcXF1rY22WJIi26tqUsJos+Bs5pllmu7H5JL329KcTHvr8KEJ9y0GtxUB2K46e2338bR0VGwbW1y1sVNVSDctBxQSiWvNtC2pm1Ct9RR5bsoXqp5NK8sy2bi07zyPHfmbXKg7Vqnp3Xh1svXxygfuuL4wlxjAbB/yygGvrxD6TqdDvI8x+7uLgaDAa5duzZ9TrY247RpG8YRHde0jcw2t52pMkG35/m+U2W7DrWZb3ulr89qRK9Y+d7vboNJ/KFJ3XXPpniYnYzzZ8vPV5ZPJrMOrsk2pJS42ojCN8C5gyV1UHGery+ONqr0YXBf+1I5YyduOjDrgmvyN8NtslQJN8utQogh2VcNwk3rw00ccLjp+PjY2abCTW6jSrgpHqnt5euLQPyHdGOMI5/RZksbM25C9Qpdh8ZUlW8R0TDOnM1N78tDt2lRFOj3++j3++j1et50IXnq4tiqzmnzP7c9zHDOPGWL64rPCQuV5UL0ipXp/eCeR6F5mL9pg7u8wDblJRQvRoFxPbhYMqBpffXnponJqygKnJ+fI8/zGVLsdDooyxJPnz5Fp9OZDlaXUuaTy7z3/PlznJ2dod/v48WLFyjLEo8ePcLW1hZu3rw5s1JVpU1DqGOyjy0PuEy4ofNWZloX6lJaXONyFZUa4abl56Znz56hKIrGuOn8/ByPHz9eCDfNE01yU9N1WUVuAtz85Lq2pbf9Dl3b+i81nEKGlH4Dm42LQi+vCBl9dYGWRdvbVm9fO9HVRVrnmFXAmHoqNVmx6nQ6ODw8lK1/S4oow0q/nYQqLa7/Nrgmw5CSQZWWkFfYjOPLjzuRhsjLV9dQGj1ItdJhDmo90PI8n3nFpqsOrvbRS/L6q9JFUUQpfuafmbc+vFwUxfS5j8dj9Ho9dDoddlvbSNv1vFxtXidp0wmeKiY2ZcT05vpk86EOo0o/k3krdIuEcNPletjCXPcXwU36eh24ydf2sVgFblo3mHMwcPmtY7Rd6DP0jVd6HeKM2Gv6PxTuQuycHVtn85pjMPqufYaTUurS685N49NmWMWMm263i06ng263u1ZzOMU8626Wxd3t4kOUYdXtdnHjxg189NFHePLkyYyyMh6Pp+RAScMnmGvgpiotvrS+SdBGDBwFRedvTjpZlk1JVCsi2hN7dnY29dia6fWkR+tnaxtfXVwKgf7ThFFFYaH52cJTnpPZ1j4SdpF8HdDPIXTPlCN2WdsEPUtQtT6mZ3EVPb8uCDfZYdsGNG9u8v1eJm5y8ZXrWa0SN9WxG2BduQkA+v0+PvaxjwGYtMPz58+nbXB6eopHjx5N456dneH09HR6bTtP4hv/pgFnjit9TZV+M37ojJWvLC2rqyxqcOR5fqlsE91ul2UcubZz+4yb0BbH0Nky2sa0XjHn1ui1dvrodLSN6LWvjW08xB3Hde8AonnRM1W2c1S2365wF6fYVodNR8Z4PL7Ud3x5cebeKMOq1+vh1q1bGI/HGI/HM0K///77ePHiBQDg6dOnePjw4bQCSk08VS5vFXfi9U2qAO8tXLYyXfJomA9cf6xtPB5PJ+5er4eyLGcGiEsB0Wlccvj+u37Te5x28CkRMYpHrJLCkclsG04nDoWlgg4wqmC4ttDY0tB41KCi11XhI8V5bpmcF4SbhJua5ib6HNaJm+qsw7pxEzDpt3qLq81Q0vwEXFa0KT+Fzjb6zkHpe2ZetB/7jAJfWbpuZngVA4TyETXiQi/GiDGsQuGhs2Ux2ytDhhW9roJFjiXKL74wjpw2Dgbitvtz83XxpG0usCHpdetFUcxsrwAw82HF8/PzqXLz9OlT9Ho9PH78GE+fPgUw8ch88MEHM15lTRzmxGtWhDvp0rS2//Q3Tdftdme8PL1ebzqIdXnmBOB64JzfMf9D93wKAVd5iPlzKS6u+y7ZzDaxhbnqGWpjWztyYCMEF0nYBrdZFlUU6jagfOCWtUjyrRPryE39fv+SsSLcJNxk5ivctHjot7zZ6l2WJc7Ozmbim8/s6dOnM4aXPjtohlMPvAmTL4C4DwTTa2oo9Xo9ZzjdIkzr5TPaaF5mmQJB25H0VkC6RYD+Nj0RertJt9vF7u4ugMngunbtGsqyxOnpKc7OzvDs2TMAE5Iwl8mVmniU6VY730RO75uy63jmNwH0OQHfBEnz9d2PmUx9ilVVxSV07foz20r/pu0Xq9Bwy3W1Cfc5hMKahE2R0fddSoyZzpaWU14Irsk8Jo9lQJPcVJaTlysINwk3+dq6SljTMMcBvefjphSPsHDTLHT/s6Esy0vGj4nxeDzzDOi2KXOrs7424VuxsvXv2NUWM0/fNjWK0PMNjSMXf9nCfbwZm5aTd2z6VLnXBZSH6Pxu/td8VoVHqrZ50oqVS2mxVcYVppUgTSibm5szS+RmIz579mxmz/GTJ09wdnY2zdPcI6nU7P7Yfr8/VaB6vd4lD7CrjqE2CMXnKDMhRcb3P0aJiblX5Y8qlSEFxRfHlo+vLdsOGzFoxC6J29L5yqXxXWUvO9rATY8fP57yUYibBoPBtAzhpvlwU4h3QuXa6roM3ET5xxdG+SGlLsJNs1Dq5XesfOdBbOEbGxsYjUbOcJoXDdNOa53OXO0yV/CByeqWaRx1u91Lq1C+ZxQaAzHcFmMQ2sJjVuJoWbZVPjO+b4tj6NyTyzh1yR3adhjbLlUQMyarjF8fX4XiUz62cZivr9BjAZw5lSL6deuciT2WwF0Kj4YmFX1/e3t7+vvs7AzPnz+fyXMwGFwqIxZa1hgi4HZol1LiC/NN7JzfPgXGdu4h5c+sAzduDLjp6FkDzrM004Z+c9PPE6Gx5Yq3KmiCmzjtORqNZu7ZuMnsf8JNfG7ihMVyEyevFAg38WDWc124iSKlntyxblt9NJVGapTRPm97oUToZQhcOTngcHgdsOXjytum31QtS7DaSFqxivFiVSERH/HqON1u99I+XgqbZ9pWHud+SGGh93yKSSi8DgXGpZjY7tWlxKQQia1tKHxLu1RJoQqMrTxXv7R5AKkiY65e2FYybOlDnsYqWAeFJIS6uckWbutT9J5wk3ATDYvhJh+WkZtoeQI/UnQms0/RPmh+YBa4/OFc34pXyngRCNYd0YYV4FYk9D3fdQpp+MqOSUN/x6YN3XcpNj6Fh6OwuP677oWUmViFpory4ivDhM97G7NNhKsMh/qPLa7tv++3T2kJxfX9ccBRpFYRwk3CTT5OctUvlZvMcPOatqmZh629Y/qMcNPqwcdL+h4nvm08+zgmxDmx4VXys+kNdYVz4leRlYbFcrkPPkM31Fd8RjPnOlXO0LWPd1KvfXWoIpsNofAkw0rDNinQycM24XBJpM7OmQLuwApdc5SaqkqLLcyl1HDu2cJD4BhAZv/w5e3rUy4FxxXmIyZXPHodo8C4FBPbvRTFpUpcWr9VVWKEm3jXy8xNMXA9S8pHsdwUW45wk3ATcPlMD+0HIR6iaX1t1SRXNWmYhXgsdAYr9jXzvvDQK+xjt0tWGes03DwTxO07tvmR1sOVhis3vWfKSetRluWlcPOafnJAvzTKdW2eG4zNm8pGXx7D4aZKhpULVIHx3ecoN64yUhAiKRdCgzwUJ6TAcBUZXz6uvHxy+OBTMqhSYsajbcydIHU6U2bzHi2X1t0nk+vaJ6NLYdG/XYpMivICvCSeOpUZrhzrgnXgJtu9VeQm31gWborjptB94SaBoP3gziGuObBq3r7xGgqjHMZN20ZUMqxclY29HwqrE/Myqjj5pcjCUSjoPbMsM44tjMrpey76A6s2GehfKJwqIi6lzKfYmf/pb9s9s/4+5cX87fvvUg5s1z5FJRSeqiS57q0ihJvS8l1lbjLLEG4SbhIIVhHm2KH84QsDLq9++fKODaNyxFzT/G1x28QZ0YZViOS5YSkKTqxcHMRY965JX99zGRmu33TydHk+bYPBp1i4FJyQPLZrszyfolHnvdBv1z2zjXzP1fYcYwYx/V9VkamifKQoKzaZVwGrxE1cw4aOdeEmPueEDCbhJuGmumDWyxWeElY1LTdOE3DxBxehuD4jITRX+GQJjWGaF91WaH6Cw1YuLcsMp1veqCz6LY6aF+mrw2mdzLRcbnLBx1Nadtd1WV7erkfrHbNdj8b1bQ0M5W3bduiqo0aUYeUihxhFhDtwUgZY7MDkxqeTnVluVQXG/G1TZlxp6CROFSruhM9RILjKBVdhob/p93toOp/s5n8zru03hev5uwjCfCZUoeEqML6wUNwq6VzhtvouI1aNm7jw9eEQN7ni18lNrv/CTcvNTT5+iimDG74K8HFRqJ6+8Cppm4avbDqOaFzf2HHl58qbhlM+9OltNtm0ck15Qt/znbmiq+k0rmkE2PLy5e07G2YLD7VxqM1NuLjJvKbcZP6ONW5irmPOXNGyQ0abDZVXrGIVFM4gTyWKWAKxxacKgy2eK4xOhqbS4Qvn/Pbd8ykp9DpVobHF5So8gFtBCaXn1C/UTua16/mYYbY+TtP6FBjz2mfIcIweuu3G9TtFebLVdZmxDtwUihfDTRzu4vxeZm5KyWte3GSbN5aRm1y/14mbBPNFyEgLGXwuQ4qG6/FnhmdZdunaxa0chHiacr55TY24ULv4Xl4Rkj00//pWz6hxQ+OHjJ/QqlKobN91CpINqzqVllRlpQn46kcNJRrGUYZ813UpM/QeR9nxxeXGT8knVU6znrbBzr3nAh1ooTDfPY4SY4bFpEtRiFz/lx3rwk02CDcJN7nCqnJTKrcINwmaRMjgcMXl5EuNIfrtLxM244bG93ErDaf1oEZdaMXKROiaw1eu+BSheTRkvPjCqeFkCzfDfEaaLb4rL5fsIUQbVrRQV0FVFBmXIpCSX13gyFRlEnUNPM79WOWmzv9VlJ3YOJy6+uK74Hu2ob7tUwRc8UJKBMf4qpLeJe+yQ7jJH6eue8JN9XCT7ZpiWbgphm/WkZuA2e1IJmxKpi1OShgnvElUUdRDYTHhtnHnG4sxBojtmho/Jnzb8UKGEw23ye2rty+t7ZobphHTT32GEo1rM6xM0LFFdQEaTvPybe/zle2qb6UzVlylJZYUYvPw5R2rCAEvO5CZvsoEaObpu5d6naLIhOLGKBmccFdY7H/fb1+9Xc/QfNbmvTr6echL41J4XNccRSlURkjJWlYIN/nTCDeFw2nYvLjJhbq4yXVdhZtijKJ15yZgUo9lNaxiuSMGVRT5kHEV+iaTL22e515ZfHnrF0j4wn08kee5U7ZYw4qGV2mzqgj14ZjtejQ8ZHilGk62cB3HN69WXrFyZZxKBBwSiCGKFPmAy5Map8yYyZKTLjRAXPG56WIVGd+90HVqPhz5bHH0tU/xpEqqhu2awjawOAqM7X7IG5IaPybtKoDLTb4w4SZeOuEm4aYq8WPSClYXXK5LSZ9qcCp1+W16NJympWMyJJdv3FNHCg33ca3NsKLhLoSMrqpGdKhNfIZTqE1txo+Z1vc8Q4aVq2yfPK0+Y8Ul1yZIuGqeqYqMLyw0WdN7nMk95b4rTorSw83bF0c/K65SUwUpZB1SNnxhscpRKA9bnFVRYoSbeBBuEm6y3RNuElSFzyioml8MbMZPbJku+cuyvGSgVKk3jR9jWFGEVsNCBmPIQVYFVYxRn9yAe5utztuXnhMe25eitwJyzjE0paikDJR5EzP1JnPi+zwYZp4p9+h9rkITE5eThpsv57ctb1sbxnqDzfCmVh5MhcJVToqyE7qOGaPLCBs31WVgCTcJN3F/2/IWblpvbgImdTo7O/OGp4Tp8KpGRFuRqtzHOjBiuCeUX8hACV3HnoOiZcXwJkWVrYCxz6oOu0EjNPeH8k7lI1dYLS+v8BVQlTBCE0ldxGCW5SuTeiGpdyNWHo7yYsaNDaui+FQhG67iEYrr8/768va1KZcAqJIT+2zrGhOcgV4lXdsnVy7mzU0hCDf5w4SbLkO4afbeqnAT4Pe6N8VF8+SpRaHqWLLFr+JMqmqkhbbrha5jZInNPyZtDKo6QasYQ3XmpVHpO1apAzrFi8bJN1YOV1yu7Ny28CFW6YmdeGPSc5Uaej9FGUrN13Y/Rm4dxlUYafwqk1DVSS5WGaqyArOMmDc3pcRJjS/ctD7cFAPhJsG6w9ZnOOOBjjVu37PpDmbaqhwQK3cK33Hy4sSPQSht1fkz5vnFlJ3CSUmGVZ0e4Do8wylxbWmpwl0XyVcZRGYescZA1Xh1KAWxSpRtoNN6p3iCfWXG5lMFthWFKvHMuE2FLwuEm+Ih3MQvd9W5CeCvgMXEazJcIEhFrAGTmja2D9PvVoVAjboYo4/KxnFIc5FquHJQ94JLimy+NI2cseIUXIeFmOc5yrJEp9OZSfvixQtkWeY90GYri+MVTkEdCkqMbDEDoC5lqE7PbEpYXZ6WJpUWLmgfaHIlZFWwLNx0enoKpZRwE7PMOuIJN9UH4aY0lOXsGavQirLZxvStZVR5pme3dJj+3+l0oJSactH5+TmUUtMXHZivFj8/P5/h0bOzs5my6TV1ZoWUaXObm+9Nbi6YDscqaDo9J9zFpzGGUSht0yvmdaVtGk2t5rvCGtkK6LtPJ1ab5fjee++h1+uh3+9jNBoBmCgqRVGg2+0CwJQYdJiJs7MzZFmG8/NzjMdjPH/+HMDLN4eMx2O8ePGCLfe8wFVQqi6pxmwfifGQVyETTh4cbwrX8xuLNpEG99nF5LcKaJqbzs/P8f7771fmJl2OcJM9vXBTHISblgNmXUKH7UOGlYnxeDxjcJiGldn3siyb/tZ/eZ7PcBQ12qih5atTKNzGrzHPl6at0u9DDoIqPBbq/2Z4yKgOhceEVU0bE6fq86kbVXkkVp+odSugWZDty8fvv/8+iqLAs2fPAABPnz7FgwcP0O/38fjxY5yengIAnjx5gk6ng6Io0Ov1oJTCK6+8gu3tbRwdHaHf78+8VpI+xKKYVCvLMhRFgX6/PyOLNrhMmU9PT/HixYup/M+fP7d6lVeB6F2D2ueN9iFGwXGVkxqvKW9vk6RgEislIXpfwMc8uanT6aDb7SZzk/4t3DSLdeSmJuPGwsVNPgVZEA+z/V68eIGPPvpoen12djbjXHn27NnM9enp6QyHPX/+3DpelFIoigJ7e3vI8xzdbhe9Xg9HR0dTQ6vb7c6sqtPXdXe7Xe+zpqtwoWv6oVbpR5fR1OoKJ21Vg9KMEytLmwyxqojeCjgej6GUmnpRXrx4MR0sX//61wEADx8+xLvvvnsp7fPnz6HUy20wZkPS5WtdzuPHj6GUwgcffAAAGA6H6HQ6uHHjBobDIUajEQ4PDwHYv1htmwQ04ZjldbvdqQdal0/TnZ6e4unTpzP3Xrx4Md3eo+tp1jkGtgnMhqbIKNZr40oXSl+nJ9pXjituVW/KKhHAqsCctE1u0pN33dykV5qWgZtsZQk3tYub6uIU4aZ2wtSdqMFBDanT09MpvwATJ495fXp6OuNYefbsmdewGg6HM0cj9Kp5lmWXjkxQwyr06m9dN/2fGk42zqMypsKXNnZlfZ7Oglj9I8apUXX8V8m7jlXEVFQxCEPpU+SKMqyeP3+OH/3RH8VgMMAv/uIvziwTm8qJKZBNkfUpLebytXmtfz9+/BgA8KM/+qNT4uj3+8jzHLdu3UKv18Orr76KPM+nxEHLNMt2NZrtC9ZFUWAwGMzco0vlOs3jx48veZW19xmYkGPKMnDTqLP8tihDZvoUb4prVakpiHIUD81N/X5fuOkCwk1utJGbUiDctBx49OgRfviHf3jq6NEr44DdIInZfkefiWkMvXjxAl/5ylem8bIsw0//9E9Pr7vd7tTwyvMcOzs72NramuZ5eHg4DQdeOq10ON3qTLcVhmT3gZ7nopxF24yG0ZV/Kod5T587M69jx7qrLFt6Xxy9ZVMb376PEcfKVgVV826aOxbBTa5nEWVYPX36FD/3cz8XXbjP2+lSZKhSYwtXSuHFixd48eIFlFL4O3/n7yDLMvzYj/0Yer0etra2sL+/jytXrmBjYwNKqZllb5+8AH+ric0bvbW1ZY3v2o5kQp8Re/r0KfI8x3g8Dn55+smTJ8jzfKocFUVh9Ww3Ca1Mnp6eotvtXvKsVUWqMuRKb+ZjKsuudFW8/KkeqjpQdSJYBjx9+hQ/+7M/G51uXtz0d//u34VS6hI3Xb16FaPRaLptmSMvwOcmW55VuUkphWfPnq0UN9Xtba2Lm0J5CjctB87Pz6eOF5th5YOvjbSxZIKuTOtxp+PT1bKyLKdnrfr9/nSLMzAxzIqimDpzqKOG6j70tylbrLOGGj+uOvrSu659efvK0KhTibetnnHloEjlkTryXuRYDhmr80YtbwW0wVQA6GCjHjZbRzIVlbIsZw5e0nC6EvHkyRM8efIEH3zwAX7hF34BX/7yl5HnOYbDIa5duzZNm2UZbty4MXNewiWTq36mzJy4Lk+PCR3GMQI1tre3Z+RYVEczn8/p6SnefffdmZWDNoISbl2kSfM1YVstaQouw2GVkMJNOp15P5abzLQhblJKOblpNBrh6tWrtXNTzGQYw036RR0cLAM3mfK1FcJNy4vxeDw1rIDF1tVsc+380YbVgwcPALyUbzAYzBhWGxsbUw7IsgyvvPLK9LooCuzv7yfLQkENxqIovI4lynu+vhVyhMQaEbE7YFxp9e4KbRCPx+MZI1ywHKj88gqf58vlVXClcd3TnhKazjYBUKVHe46BiVf7vffemwn78R//cRRFgddee216wPPGjRtTbws9XE7rTO+72iyEqkTr86Lb4s1j4ux2uzg6OsK7775baeWqSrumlDUPBYuOA0F1pHCTjWtiuckWJ4Wbvv71rws3QbjJV5Zw0/JCG/GaH3zgroBqzvM5U3xlmTLptOa9LMumO2iUmqz6moaV3lWj1GT16smTJzN1NFfMbfxnc06Z177dBD6eCLUxDbPFTVnRrQpzEUEblvo/Z5WNW0YKYvtsnXmHsMiybajlrYCcZVnAPnG6PMZ0wLtkcXmZOWWUZTn12PzkT/7kdDB++ctfBgD0+31cvXp1mu7mzZvo9XrT6+FwGPVQYo1QmtZW91j4PJW2cqpAqcke7itXruCjjz7Cw4cPkwZAHWTCrVfTy9XzKstW3jzKXCS437FaBW4aDAa4cuXKNN2qc1Pdk59wk18G4ab6YfITNeZDRkCIa3xpfS+g0Hm7nAv6hRsaDx8+nAn/R//oH83kZb54p9vtYn9/f2qU9Xq9KS8ppTAajWbOhQ6HwxlDjG4zNGU0nVfaEKHbFGM5w+c0m5dxpZSaaZPxeDz9dJB+URx90+KqrPouUva6y07+jpXNq+WbUF2eZBqHKi5VJjubZ4QqN2Zc/VuTyenpKX7qp35qmv5nfuZnpvlkWYZr167NHCS/evXqdNsLMFm61q+M5xijPi8MR0Gso3O42jtVedJku7OzgzzP8cEHH1w6KDoP1D1phzxm85JhnuUtC9aFm/QrmjncdO3atZmzVcvGTSHlkcYVbhJuWiZUUY7rfNYhA5uWRY0dumqvXxqkjR/g5ct3tEGn8xuPx5c+R2EahVQ2vTXQ5qTSL+NYdiilZlYEQwj1hZhnP28sUra6y44+Y2XzKtosfZty40pn8+rSvMzrkFJjG2gcz7Gp7FClRntEzLr+/M///KVrU4bBYIButzuj0Ozu7qLf72N3d3d6bzQazXwvInUlids5UjzKPmXT1p42bG5uAgA+/PDDS1+Lr4pFKUNVFYQqci+y7LZhHbnJTC/cdPmecJNwk2B+KMty+nIb4OVnKUxe1DxlGlca9Hp7e3tmBWxra2vmw+u7u7vTNPqbXBr6bYcatjepLgOyLMPGxgbG4zG63S6eP39u/YC8oF2IXrGywaaMcOLRMJsi4cvDppzYjCKfXD5lSefhm/jNa/oNnMePH+Px48fT79wAwC/8wi/MxFFKYXd399KrkpVSODg4ADBReFwvstjc3Jx6erQH1gafV7qOCdSn1JrXm5ubGA6HePfdd2e+zVEVyzqRp8pdh7x1KWBtxjpwE8d4EG6avSfcFIZwU3sQGhMu/qIOp9CzsXEITW8bT674piFlvsLctjJMdS36IeRnz55Nw/VbVE3DqtPpTOU3P/au86arX2ZZPsNLqcm2RlpmU7C1v5bffHusNlTpFvjUMZPi1JoXQnVqUm5f2a6w5K2AIQXEV6hLEfB5kmmZ1FNMy6PeZJqHTcFxxfN1OKXUDEnYyvURIQC899571np+5StfsaY17+3u7uLs7AydTmdKAC4opXDv3j0cHR05ZanaQX2kq6/zPMfR0REePXqEBw8etGLyrCLDIshokUpPW7HO3BSqj6vceXGTVnxcEG5yQ7hp+UENHBM+Q8YWh+PcMdPZnDyuND4uscnqO9fkMpx87WCmMZ0+APD+++875QIuv0XQl3ev15tx9PR6vZn0pqFVFAWuXr2KTqeDw8PDS6tjtrKqwDaX6DL12xvH4zHOzs5weno6s4LFfTOuC4vkuzqcaG1BtGHFVS5CebhgUwRcRGDzIrsUEJ9yYlM4XAoJTUffVuhL62qLWGIzDTlNNiElR6nJ9h+95cWFOhUZH/I8x/b2Ns7Pz/Hhhx82Vs48EDuobf12XmgbAdUJ4SZljedrkzq5icr2/vvvex1S5u/hcMjmpqbHjOam8Xh86bD+siGVmxaBVeamdQR3fkudB33cFHIY0Q8In56ezqQxzzNp+ZRSC99SqNTLs1ZlWda+bVlQD5K2AtJJMbSi47v2xbd5FTnpQx5j89qmmFC48rF5f33l+K5d6V1hoXxtct+8edP6imYX5qH4b21tIc9zfPTRRzMfMlxl+DznPqyyR7cuCDddDnflVTc32fJxKTtVuGleBpY+47Gu3BQD4SYBRR3OrKZwdnY2Y1jpb0dpmG8k1Ktb3W436tuiTSDLMvR6PYzH4+nKlclNvpVRwfxQ64pVCnxeYJviw10GN721NoXI59FweZt9Rowrjim375qmt9UpJKMrHJi8nvnWrVvOclxwtWNdyPMcW1tb07MNcjDTDc64W2cFR7hpOblpOBzi5s2bznJcmBc3DQYDfO1rXxNu8kC4qRp8bdNku4XGjmsMc2QLpY2FL31s3iEDRO8K6vf70z/uR9Gr8lGojXX+enu13v5XlpNPdFCH/7Kgirxt45bkM1Y+cD3E9F7sRO2Ty+UddaXzeVNDCok5uXO8qTbCcZXhkpOWF6pz7GqVT94mUBTF9IOdosCkIzQ+20ZAdUK46WX+wk31QZ+tEOOqGtaZm5Sq94UHrrFVNS3HgWKCjuuQXDQ8Jj09P1TFqOMYlFmWTR0ru7u701fGhxx4nPCqMF9kMRgMpi8G0d/j06txMYbwMqNtRqT7tJ8H1GsbGoguUgndC5GRrTGpte6z3qkX2SUTJy03XSgspoPQQe5S/vI8n77Fq0pZTXfeoihwcnKC0WjUaDnrDP0czb9Vwqpwky1dndzEKY973xWXw01FUWB/f79yWU33406nI9zUMFadmwTLAaUUOp0Out3u9K+t0CtX+i/P8+l/+oFlwfyQZFjZCM9UZDie29gHHuvxoVtEOJ4Z875LEXDd87WJbbKwyWMrMyQHR87Dw0McHh5eihfyrNjQ9ISXZRn29/dnvqUjEHCxKtzkkrEObqJhbeWmWMyLm3Z2dhorQyAQLB79fh+DwQCDwQC9Xq+1BopSk5dZFEUxNabMv7bKveqI2gpoUyB8W0xc1648qEfYBFVEXJMv3bZiKiO+/M1ybHDJbOZLf/vK8oWnKhY2dDodvP7669avdldRQuqUkSLLMmxvbwOYfLCz6itEBauPtnCTK455b9HcZOa1aG66e/eul5tSDawmuWlnZwdKKeEmwdwQ67gRVEOv10O328X29vbUsGpLu7occUpNPrdjznuPHz/G2dmZ8NScUcsHgrlbaUL3XQZTzMqKyxscSu9SrEzYDCfOYKs60cemN+PfunULx8fHyWXbYCpu+roJbG9vT19qsS5v5RLUi3lzk4+rhJua5yazHEC4SdB+xHBUTB6rhirO6hgoNXlpRa/Xm34AnRpWbdqmaspC31yoz4SenZ3J9to5InoroO5gscqJLy8Nn+eX43F2bZMJbZUJyeWTjatUxSCUD3cbTpZlOD4+Dj6T2AE3z8HZ6XRwdHSEXq83tzIFywnhJjfaxk15njfCTfOEcJNAsJqgWwGXBeYcqFewtJHY5i2Nq4Zow8rmFeAqFpzJ33Yda7y4yvfFMcFReDh1d9WZc/aAszUohOPjYxwcHMxFMWmyDH1wfGdnx7ptSCAAhJt8edJ8m+QmDl8JNwkEgjZCGyTaGNHb65YRnU5n5kUcy1qPZUPSd6xc16Yy4/OiciZ9n6Lgm9hTO46Wl+ORNeP4vNpmnimymNe0rBCyLHOeraob9MxIE1BKYWdnB6PRSLbfCKwQbmqem2xlhsqiyPPcebaqbgg3CdqCLHP7sWMdO76wGN5JcSrVFTcWdedtcqk+g7S1tYVut4srV66g3++jLMvpq8ybKNeG2Hr6nne3250648qynH4Y2bU1sK27A5YNSWesbMpJ7OD2KTiuh2s7n1AHbEoBZ79uigx1bsnxQXuEY1CH8tH02YZOp4Pj42M8ePAAjx49aqQMwfJCuMl/n1uWC9xthz6sMjcdHR3hwYMHePz4cSNlCFYPnPEUY1jFhs8TVY2GpqANK711bmtra/ptvVh9rUknTgi03KJ4qeKXZYmiKKCUwng8hlLq0kstxLCqB8krVpytMTSMm7cNqRM6Z9uK7/yELjskH7eOVTyooQPq+l6e57h3797UI1xH+8egaVIpigI7Ozsoy1IUGMEU68pNdSlWVQ5/U26i6U1uMlfSV42bOp3O9FMRwk0CwXJgNBohyzJsbm6i1+t5VxeXFUopDIdDnJ+foygKnJ+fC0c1hMpvBeRsv3GlCRkqIaUi9p6vrFA4VcpcigM3P/rb5/HmtoW+Pjo6mvEIc40401iry+vSlPemKArs7u6iKAo8evRo+qVxgUBjnbjJV7cmuSk2X7paJdwkEAgWDf0dKPpWvVWD/oDweDzG+fl5bTuoBLNI3gpouxeaqGInZc6k5zJCOHnFKD0pipYLNu+uTRnyeYBd+Q4GA5ycnFw6vxC7DaaKwlHF8x0DrcAMh0P88i//shCEIJmbfGm5cUNxVpGbbOXbytLcdHx8bOWmOo1NH+bNTYPBAF/96leFmwTJO15i+2mVscQZy3WU2waY8m5sbKDb7eL4+Bi9Xm9lX0SjtwbqFSv9d3p6irIs8fz58wVLuBqoZb3T583UEzPXQ2yLa8uDk6frOhUcbzaNW8dBxRj5b926hXv37lXOZ5nQ7XZxdHQ0s59YIAD43OQzZISb4uK4jMfbt287uWlV0ev1hJsEjSOkZ/nixqT15bfMyPMcnU5n5V9LnmUZsixDnufI8xxFUUxX61a1zotA0nesbP+r5AG4vxOjPa6pB7ZtSoRPQfLFMc9H+erikonbZqlbVLIsw6uvvhqdbh6wtV1dUEphMBjg6OgIg8Gg9vwFywHhpvZyk1JKuEm4SSBoJUajETY2NrC5uYmNjY21MDKUmpy5Go1GGI1Gwk81orETelwvML1vwty+lqIscb3INmiFyVcHU0Zunhy4lKZQfkdHR9ja2vLmyZHBfD1nVXAVzzrQ7Xaxv7+P4XDYWBmC5YdwE59jUtPR+8fHx9jc3GSV4UNdvGTLax7cJMqLoA6EVtsFPOj206s3eiVnHdpUKTWzgpVlmfSrmlBpf0JoYjfvufbl2+67zgPFbKOxxeUoIyn3XHWLkdcWJ0bmLMtw//5969tsXG3NlaMKzLLn8RrSoihweHiIhw8f4qOPPsLZ2Vmj5QnaCeEmf9046X1xhJviURQFjo6O8PDhQzx8+FBearFmSB17ISWXk77JvFPTV+HQ2LIoynLyfSr93apXXnkFW1tbyPN8+lKHRRkWMQ6xqnnrena73emr2HXbnJ2d4cmTJwBw6bXsITnXHZU/EMy9b/43JzJOfqZCQxFzKJ07cGOVGY7ikoKYV7MfHBzg6OjIm0+sbHUoG/RtXvOAUgrb29vo9/v42te+JsbVmkC4yZ/vorjp8PAQh4eH3nyEmwSrjKoGTNWy25r3oowXXbZeoep0OkvxNsCmHEC6fxZFgbIsL701cJ4cuQpI3gqoH4RJGNxJ3EYyrrQhxYMT1yVDzH1uPJtHICaP1M57fHzs/PZCqvelbkXMVFzngV6vh8PDw6UgTEF9EG6qHk+4qVkINwkEi0Ov18NoNJqeMVrF71bFIssydLtddLtd9Hq96SqW/pM24iOqpUxlxTUBmeE2hSNmmdN2Hev5sZXJ8fqailDMJE7jhrzftvLo/VDanZ0d3Lx50xmeojDUrWC4nn3Tikyv18OVK1dqOd8haC8WzU1m/jF5zJObKObFTTdu3PDKINwk3CQQzBP6bJH+btUiV87aBOqUNM9fiWHFR/IZK64RZEvHnchTZdH3bNtMXIoVBZ1Uy7JElmXetLa8uXVJXeIdDAb4Vb/qVzXywoY6D3bTbTepW4BSkGUZ9vb2AACPHz+27hcWrA4WwU3c+9zVMB+Em4SbBKuDOvShJvvqKhodg8EA/X5/erZqMBisZD1jYRpUWZbNrKg/fvx4+sFz2zk02Sr4ErV+x8q8TlFAYr29rvx8XmmaxizXdd+Xpy0dBzYvKdd7bd4fDofodrtJ5XHSuORIxTyVFhNKKezv71s/UCpYXbSNm7h50HIXyU0xK2vCTfHQ3HR0dCTcJJgr1tWYyLIMRVFMz1aFHGDrCLp6ZRpc+g2CAjsqG1ZVlYdQHJsSAVT7uKXZUXyweUVdW0ZS2sGnLPnSaeR5jo9//OOsV/imeBOaOrRoKoLz9nLoD3ZyFD7BcqOt3BQqg8MFNrma4ibbtSudhnBTPPr9Pg4PD4WbBIKGMRqNsL29jZ2dHWxvb4tDwwFzPup2uxgMBhgOhxgMBtJmHiQZVnTyN/+nWLGudLZ7MW/sct1PnTC59XMpR3UrekdHR843AbryTKl73UqG+UYxem8e0GcbNjY25lamYD5YF25yvTY3BOEmPxbNTf1+X7hpRWFbAaD8lBJedeVgEfnWsdqRkp9+A2Cv18NgMJhuoRbYQdtYf+tL/ze/fSV4ieTXrVMvKVf5cIXFvE3PFidGmeHGjZXfbBOf4sWVxwadR57nuH//flLaGHAPuFfFvAemUpPtN0opOduwIlg0N7n4KKYs4SY+hJsEy4SYMcsNqxrOSVslfSh+nXn70p6fn6MsS2xsbGA4HOLVV1+dOn5snzxo6rXmdSDmJUtVQJ99r9dDr9ebXj969AjPnz/H6enp9Jt81CFV51nYZUL0ihVXaalCIty0VJEK5RHK1/c2sSp5+8p0pQu9pSrWI6zLi8U8BsOiPEZagZGzDauFRXGTLR8fN8UYX23iJhtMReT4+Fi4qSJMbpK3cQkE1aDPVOm/dVLym4K5clUUye/BW0kkbwVMSeOb+G1hNqWBxnN5X20IeSE4Rpttq4itLqH8Q/KGlMLr16835u1ZBEzvxiLONpycnMx4YwTLh9TVixRucsWZBze57jXBTVyY6V977bWl4x8fFs1NV65cEW4SCCpAf7dK/4khUB3dbnemTcU5/RKVPhBs++0Kr3KuyYSe5DhKCGfZ2aVshA4wp3iGbW8B9L290IXhcIiTk5OospYJi1DKOp0Odnd30e/35162oF7Mg5tsvDEvbtJh8+Cm2HxGo1GQm5YZi+QmMa4EKaCOH+oE8oXX2d8X6WxRavJGu36/L4ZVIujz0+erzJVA/SFhvcreRD9aBlT6jpXv7AHdjmNTZmK31biufTKZ8bll+uKnenxt+aS2we3bt5Mm2ZTXCc9rj6z57Ba1v7nf7+P4+Bhf//rX8fjx47mXL6gHi+KmWA4w7wk3tZebaJmL4qaTkxPhpiVGlX6Twi2LxqJlMrfQ6q2AOzs7ODo6EgdqIsw2Nc/UApOzaqenpzg9PZ065hb1EqBFI3nFih4cpgoGZ/uMD7bJ3pYnLd+3/YariPjip35vJRQWajNd7nA4xK1bt6JJq4rnYJ4eB9/h+nlAKYWDgwM527DEWBQ3ueJxxnUd3BQL4SZ+GYBwk0CwrCiKAv1+f/oCBhk/9aMoCnS73enfog3rRSLp5RWhySz1rSWuJWifQcPxIPqUnhQZY8Ex0EKTs87j9u3byd6WKq8zblp5cL1NxlQk56nADIdDHB4eCgEvEebFTfR+W7gpBRwDjctNd+7cmet2NeEm4SaBgAN6Hki2AtYP/Qr7wWCAfr+/1oZVLb0rdQuKL27K69dpfJ9ik3K2KbasmLJ99zVOTk5w9+7dJDlTt7CEzpo1AaowUdnntR1nMBjg6OgIH374IZ4+fdp4eYL60QQ3VclLx58nN8XcTx1bV65cwZ07d4SbhJsECahzBb2ucjnxm5StCsqyxHg8nq6e7O7u4vj4GIPBgPX5gnU2CnwIbbFXSk2/bzUcDqevuT8/P8fz58/XZmtgo2Y75yC3K05I+dBxfPdivcE+hYbrtQ2Vkzpg8zzHvXv30O12k9JXwbzPPoW2O82T9Pr9PrrdLh4+fIiHDx/KN2VWBOvETdxyYsvXEG4SbhKE0ZTxVNW4mecKeZW8Y4w6rczneY7hcIj9/X288sor07BVNZwWWTddrl4N3NjYQFmWODs7w/n5+czZKy1rytnaZUClvQS+bTehLTmciShW+eDAZnH7thCFHnyKguSK78ur2+1if3/fmW/TSFXiqsA8ALlI70aWZdjZ2cHBwYFsv1kSCDfNh5uUUuj1etjb23Pm2zSEm4SbBAITeZ6j3++j3+9jMBjI1r8FQL8lsNPpTL8dppSavk1wlZH8HatUL0cobSg9PZjOTecqOzQx+raaxHSO1DqXZYk8z/GpT31qIR5hE5xn10SZJhapyAyHQ1FgWg7hJl55sbLZYHJTp9Nhl9cEFslNpmEn3CQQLB5ZlqHX601fWCHfWJo/9IeDzY8y+xyFq4RKZrxtYg81mGtLC2dJ0KdIhOT0le1L44oT2zF89fOFKTXZs7pIj7DGIpZtaZmLHpDD4RBHR0d4+PAhnjx5slBZBLOgY3be3JSCZeemLMuwu7u78HG5SG7SWHQbaG6Sc1eCdUee5+h2u9jY2MD+/j6Gw+HCx+eqw+YY1MiyDJubm9N74/EYjx49AoDpFuZVOnMVbViZjcedzDgd2jVxU4UjZtuLKw7dOhKzL5Vj3NW5zzXLMrzzzjuVVqvqkiXVsK0LbdmPq882vPfee3j69KmcbWghFsFNdZQh3LScaCM3PXnyZKWUlXXFovvUMkKpyYdre70eNjc3F77bZx1hcqJSavo2a33uSvPTovXKJpC8FdD8bwtzhfvi0/uuPf+2P27+Ol8z79jVEKr82OTjyGLGd6U5Pj7GjRs3Ki1l19FpF32WoG3IsgyHh4cYjUaLFkVgYFm4yRa2SG7yybgM3CR4Cc1NGxsbixZFIFgIOp0Otre3sbW1ha2trbl+BkJgh55LsixDnufTb4r1er2FbyWvG0mGlen1dCkrLoUgVVEwy/Yp+Wb+sd5q3xkJW3xal1A6l8yuLSVKKezv79dy8HLZlQ/qEW7DwXEAM69xFSwebeYmVzpu3incFCozJIuPmw4ODoSbINwkELQNVHGXl1csHtThmOf59G/VzoYmbwW0bSnhTv4hrzJHMdETl02GkOxNLK1X2XLkqkOWZbh+/Xp14Yx866j7Ire92PrYIl8xmmUZBoMBsizDeDzGixcvFiKHoP3cxJG97n5cZWXOFmZy02uvvVZNOJKvcFO9EG5qD2J31nDDmg5vc9km9PYyfa7q8PAQ169fR7fbxXg8jl6trwqOHrqIsn1oeu4xjSel1MyZqxcvXuD8/Bzn5+cYj8cLd0pVRbKZaPOIchSPmLxdg8GcwDieaNc9G1LestXEIFJK4fbt2xgOh0npbfnVNXAWpSjoss3zJ4tUXEz0ej2cnJyId7gFEG4Kp4kp15buzp07reSmRUO4SeAC9djbeGIRMq1C2eYKcafTQa/Xw2g0mjlbFbOSvOyKfQqarrPZ37Msm76KXb+O3Vy5agNvVkGUYcXZOpO65cSM43vbUhXPjy+Nj+RccjWJk5MTfOITn6j8mlCbBzU1H9/1PFCXctwUzLMNy04Mywbhpvly0zvvvFMLN5lYdm5qM4SbBKuMoiiwu7s7/RuNRtPvJXHmAMFioT/kPBqNLhnEy4hKK1Yuj5zLk5sC+mYRmn/IK8j1VvvCbdtLQp7rEHyKklIK9+7dq+XbC3Xt9zfTL+r8AJWhjSSZZRkODg6wubm5cvuGlwU+brKh7dzElcHMl8oVgxA33b9/vzZuMv9Xzcd1PQ9QbmojsizD/v6+cJNg5UBXQPI8tzqk2jo21x1KzZ65aqNuF4OkM1Yhxdqc7G3KRUjxCXkYuG/NMr25IQMw9kGaSoGvfr56us6CjEYjbG1tRckTgtlOdRglbej4ZhsD9dWtDuzt7WFjYwNf/epXMR6PFy3OWoDDTTTuMnBTLISbFj/+AVzipjZAKSXcJFgJUC7Wqx5Xr169tO21LeNPYEeWZej3+yjLEoPBAM+ePUOn08HZ2RnG4/H0b1kQ7bayeWl9cVNWr2x52xQm7rYbzqDyvf3K51UO5R27NSjLMty+fbvR13gvq+Li6hc0Tlu8Ut1uFycnJ9PvNwiaxbpwE5VlntxU59mqWJnmkb5KuSFuct1bBISbBKsEpV6+xrvb7dayoi6YHzR/Zlk2fY7m/2VD0jsoqWKQ4mX0xbdNPinKSIziEKqTzata5WC4S76DgwO88cYb3rTrCp8CaR4Wb5N3qtPp4PDwEB9++CE++uij1ihWqwrhJt5qXYp8h4eHuH//fjD/dYRwk0CwOBRFgf39fWxvb2Nvb0/68pKj2+3OrFh99NFHePr06aLFYiP5desanG0lvvSx4bpMX9yqk1doiww3TqyMSimcnJwshYXeli13JkLbqxaFPM+xt7eHbreL9957T0i/IVTlJlse3DCzTF/ceXCTr5x14KY2ou3c1Ol08P777ws3LQht6hPLBnOlqtPpIMsynJ+fS19eYujxkGUZyrKcrmQB9Z3NbRK1fDUt1qgyvXcpe9CrKC3zMP5Sw0ejEW7evOkXbsFo01Y7DbP/tU02ExsbG8iyDB9++CGeP3++aHHWAm3hJk4+beemGzdueNMuGqkrlE1iWbhpc3MTeZ4LNzWI1H65qDHf5rLLssSLFy9QFAW2t7exu7uL119/HUopvHjxYrqtLBVt4pB1gW2nRpZl6HQ6AIB+v4/xeIzz83M8ffoU5+fn07jUkI7Zmt8EKrkftQLi8nCaW1RoGMCbaGj+sdt0YuLUsVWjSvpbt25hMBi0egKmimfbZLW1f5vkHA6HODo6krMNDaNt3MRBm7np9u3bS8NN5nWbINwkENQH8y1y+k9W1FcT5nkr/Rr9NiOpF7omCB1mU058q0y+MF8eVDnyraaE4nAUqpByU8XjvL+/P12tamunabPnlSoobT44nuc5Dg8PRYFpAMJN7vAQfNykV6uWhZvaMtYBHje1BZqber3eokURCJzI8xzb29vY3t7Gzs5Ooy/7EiwenU4Hg8Fg+tf2l5NU+o7VTEaGp8B3WNf0JIe2lZkeSJ/CFFKSXEqP7Tqk3KTC5z0HgDfeeGP6tq22Try2rVFtU7RC/a4t0ArM7u5uq+RaBfjGt3DTZYS46c033xRuqgG+ftcm5HmOo6Mj4aY1wDI/X6Umr1jv9XrodrveXQqC5YY51+rn3ul0UBQFiqJo3XOv5YwVYFcgbPvd6Rv16vCyhvbV+8qistjy4pQRgit9WZbY2NjAwcGBU8Y60MTZgzZ1ZBO2lxa07ewF8NLrps9dnZ2dLVqklccycROVc1HctL+/75WxKoSbhJvWHaHn37b+0QbkeY5+v49XX30V3W5XtgGuOJRSKIqJybK9vY2yLKdnqz744IPpGwQBLLwvRJfu8tLavLy+VYOQlzTkMfbJYiuTE8clM1W4fLKkhN++fbvxbWF1EvMybG1p45YgFzY3N3HlyhXZflMRrj6+rNxEw1KNqirhd+7cWWpuaiOWjZtOTk6EmwStglIKnU4HnU5Hvlu1hqC7S/TZK/PtgYtcxUo260LnBXxpOGcNfHlxG8z0CMbIaSJ2Kw8N88kAvHwT4DJ6pFLemtYUlt3jJ9tvmgenXZeVm2LB5aYbN2402h+bXAVrC5adm4qiEG4StApZlmFrawubm5vY3Nxs9KPlgvZBv6xEbwXsdrvo9/sz20IXicofCDaNJe7WGZvnlZMHJ67Ni1ulkUNl+GTmGB56tarJ7SBNbLOhRvKiO3JIkXLJ2aZtOHr7jVIKH3zwQauUw2VBG7jJF75M3KRXq9o0RkJYJW5qE4Sb2g2uM4fDDW2HPl+zvb0tL4BaY5j9t9Pp4Pz8fPrdK71F8Pz8nLXDpG4kGVahg9saIY+oK08OuIpSSrmm3DFKWSyGwyFu3bpVOR8XqJJRVxn0ebaJqEODiK5AtkFmiq2tLSil8PDhQ5yeni5anKVCG7jJh2XhJvObek2PkSa5qe78m0Rod0QbINyUjipjtS4HjI0bmyq7ar70HCIwWakaDAYYjUa4ffv29GPAdck8D1Q5UtI0qu6CmBdoPx6NRjO67vn5Oc7Pz63f5JuHkVVpKyBXQJeym+q5TdmGaMvDplibsroU75iJmrNaNY/tMDZPfhW4FL5Fg2s4tUVeFzY3N3F8fIxut7toUZYOdXATjRNTdh15rDo3ceWoml9bnSeAfztpmyHcJJg3lJq8uEBv/5IzVQIfiqKYnsFbRF9J/o5VaMLiKN4xB8DN8nxbYKrKzpVJK28+Jc7XPkVR4ODgYCGTaBNltvFAtu3Z2J59m2Q2oc82yDc64lAHN3HKMMsSbmov2shNgH3lf5m46fDwULhJMDfoF1XoP4HABqXU9LyV/j9vRG8FpF5TGhaTB2cpmIu6VmHMPcj6nq8cn9eRLmWbcY+Pj3F4eFhZ5hBonZrcDmj+byOoglt3mzSBoiiwt7cHpRQePXq0aHFaDf08fVvoYvKqIw5Qr3LcBDdRnJycTD//0CSEm2Zh22HQZpk7nY5wk6Ax0JXnLMumTp95raYLlgOmU6osSxRFMZ1TxuMxXrx4gbIsp69jb9phVfnlFSbo5OhScOrarhITxyeTLazucs24+/v7uH///tzetV/3JG0zQNsOl1G1DMjzfDqZPHjwQL4p48G8uIkbry5uMvMK5VllnB8cHMyNm2L5dlXh46a2t0ue59jf3xduEtQOm2HV6XRw5coV65vf2j5WBM3Btf270+ng7OwMz58/x3g8nr7Mgp7Lq1sXTJ49bRXhdvS6FBKqKHG2vVRdJePK5fOaj0aj6Uc3m0YTxsMyvlmItgPX698mbGxsYHd3d/qRPIEd68pNnHgcbtrb22OV01Ysm7HG5aa2Qikl3CRoHP1+f+ZPIOAgy7LptsB+v49Op9N4mcksaJJ/U2/Nsim9sQpRHXlw09O8bF7yoiimb7NZlonThmUyqkyEVqraXqfRaASlFD788MPp8rZgFuvKTdx4wk3txDKtotugz1s9fPjQ+jYugSAVSin0er2pgiwGvIAL3Xf0atXp6Wnj/JTcO7me1hjlxtwjadu2YyszpIhwvX91buvxvdXrlVdemcvZqnnB9mzbrti4FJi2y60xHA4xGAzw8OFDPHjwYKmVsSYg3OSGj5teffXVuZytmheWZTxr+LZXL0tdRqMRhsMhPvzwQ3z44YfCTQZSn1+d4z82XdPhPujtWnqV4fr16+j3+9NvFc3rKIVg+UD7XZ7nyLJsek4PwPTM1Xg8nvl0BN1hktKHk15e0VSaqi/FMOPH7t+va+Ly5TEYDPCJT3xioYTQxCFxmmebFQDXthug3XJTKKWwvb0NYOIh1ocy1xnCTeGyXRgOh3jnnXdWnpvaDpOPlpmbdnZ2oJQSbpoDmjRumoZLNlO51W8DPDw8XMgb3gTLDaXUdF7TBlZZTr51pQ0q82xoHc6gxmdRav1VOcNk87S6yjQVl7reGMiBr6xXXnll+jabeZEdrX/VcttM0imgqw/Lhu3tbZycnIj3LgGL4iYzjXCTcJMLwk2CdYXevjUYDKZnY1ZtfAsWgyzLZs7qNbGtNJrxQudTTIPG5pV1eXxjBg3HQxPjsaxbmbAdVh8MBnjttdcubQNqErRt6yhzGSd4HziKbtvr3Ol0cHh4OJdDmW3GsnBTbPymuWk4HAo3tRDCTYJ1RlEU0w+9ypkqQV1QavZj0004fmrdCmju3/e9zpj+5uQdA64XlFtelTMYGlevXsXBwcFC3/hUZ5m257ds226A2a1ZrvMzy4DBYICjoyO89957ePbs2aLFWQjqOoug+0LdPALwFeF5c9P+/v7KcJOZZ9sNDx+EmwTrDqUUrl69Oj1nJaufghRQrtT8WRQFBoPBtF+VZYlnz57h9PQU5+fnyfNHrb2UTuC2ycD87Zv4qCc51mtc96RT5U1yWZbh+vXrjcnmQ1OKhU0JsymkbYe5sqn/07otiwLT6XRwfHyMnZ0d8fA5wFnJDimtvlWuENrITa+99lpjsvkwb54QblochJsEKVBKYTAYYDQaTd9WOm+eEqweTB4tigLdbhedTgedTgd5nlfuY7UyXMwKVCgs5o1ZsSsnXE8mV4EK5XVycoKDg4OV8QTrfF2KyqqR3jJ5i5WaHBwfDod499135YOdF0jhJtequnBTPWi6zHXkpjZDuEkQg83NTfR6Pezu7soZK0Fj0B+e1vqs/q0/KHx2dhbtkKtlxcp1fsEX35ZWI+ZAOd2u4yufyhmCbmiuomMrezgc4u2330ae53PzmIbq1sR5hmXyBJsIKaim4qKf3zIQfLfbxdHR0dp7h6twk+16lbhpNBqtJDfNI895IIWblgHCTQIbaB/WKwj6rNUyzLuC5YR+Fbv+U0rN/I5FrczmUxpiFFLXlkJ6hsvMW8PnTU5Ris2VGVvakIf54OAAe3t7TpmahE/mOlFla9Si4VKUV+H8mFZgHjx4gCdPnixanIWiLm6iaYSb0iDcFIZw0+qgytaiKs821pE0z7KB2b6sz7QMBgMURYG7d+9ia2sLnU5n+nHXumRrK5ZxLGtU2ZI+b/jmHv0tPl2fBw8e4PT0FGdnZ9Nvq4VQecUqdvBw4NpHThUXU4aQV5pu5fKB5uHL03dWLMsy3Llzx1tWE7BtXVvEeYZVg+2sQ9uhv/+xtbW1aFHmjia4ieYt3BQH4aZmINy0vjBXpOnfssFcKej1euj3+4sWSVAD2sxJdJxkWTZ9Y6B+e2DMbhIgccXKXFHS/30F0omU412lSgEt05YfVw7udhRuPFe5+vzCIrwQVLYmOzb12gPt81CE4Ftp0PdscZcBSins7u6iKAo8fPhwpc82NMFNFKvCTfpNgIviJqD5sSTc1G6sEzcJ/BgOh8jzHMPhcLoFUN4CKFg0er0eiqLA6ekpyrLEkydPgitXtZ2xct3nLt/aziO48nDl45s4OZ5gjpzc8jY2NvArfsWvmHmN47wnPNqOWo66y2gi33nDtdVGKyrL5Am2QSmFra2ttTvbUAc3udIJN6WDGgNNlrHM22uA9eGmw8PDteImgX21Tc5UCdoEvYKa5/nMmyl9iGIx2zYUGg7YV0u4yoDvzACNY5PBLN+MH0Lqyo7Ncw1M9mnqLQ6Lmtht226alINuP1pmYqT9TN+z/ebUtS3toc82PHz4EI8ePVq0OLVBuMktD5VlY2Nj7bjJLKMtYzEVq8pNvV5vJblJYAflbL0KcOXKFWxtbaHX612KLxA0Adq3TG7V/VDfy7IM4/EYz549c/bJqBWr0MTuInufImL+1gQfGkCmd4MrMydPbvk22c3rPM9x7do1a9x5wLVK1YQcPm/qMsK2AmdTymPr2qb26Ha72N/fX6mzDcvKTU3EFW7CTN70d5vGYgyEmwSrDKXU9GxLp9NpVb8UrC/MOUsp3tsCo1esQuEx3q+U7QspXscYA4yTV8hzfXx8jNdff70V3uB5yhDyqrfFMxoLqqC52rZJJbFuKDU527AqWFZuqhNcbrp79+5Cucn8vQxjpc0QbhKsCjY3NzEYDLC3t4ft7W05XyVoJQaDAcqynBpYNkQZVlxlg6tgpHgtmzCqbJ5AF0LbcvI8xxtvvHEpvuu6KcR4uOsAfS62ei7DpO56rralYlf9lkVh1ApMnueLFqUyhJuWh5s02sRNywDhptVFk89kHs875FQNpQUmW64Gg8HaPHNBe+GbF/Vv822BFLWuWNE4ofipXrSmtjfEpPF5hA8PD6dxKNnMU5lYFJZl8g7BppC5wnzp2gyl1EocGBduCsc9OTlZS25axVUy4abVwaK3a1bJO8Spvrz16nqn00Ge57hx4wZOTk6glML5+Xnyx1kFgjpB+6DmJN1vbai01qoLjFkZMePGpjPh847MazBSGYqiwP379+e6UsQB3XffVBmrBtu5Ba7nt03Pfx0h3HSZm+7du7eW3LSKEG4SLDv0m9b0X9u4SSBIRfLLK2KUCV8+nO0NprLjKl/fo0t4PsR4ijgyjEYjHB4etuKgNJ1wY9qlKmK2Ly0LbHVyKYWrVO9lxLpwU4wMGxsbwk1YTaNCuEmwjCiKAv1+f/ohYP2Kdd/ZFYFgGRDVe20TYMrWAuoZtoFOACkTQshrx1GiuJNTURR44403WjFx0+c0ry0/NN95KkxNIeQFtvWjNvSBdca6cBMX5kr6otEWbloFCDcJlh1lWWIwGGBrawvdbnfKv9JPBcuMShuYU84h1H24OzZv12oKtyxXGUdHR7h+/br17MI8YZY/75Uj257/kJLadnAVXFFa2gXhppc4Pj4WbvJw07JCuEmwzNCOlZOTExwfH2Nzc1NeXCFYCSQZVq4JOjRxpygjOl3Mtj1unhzQOtnSmecXYvKuG1RxWcRbvtqw1ahu+JRA1xtjBIuBcNMshJsul6+vVwG+fi3cJGgz8jxHt9tFp9NBr9eT7X+ClUF0T7Z5HTVsE1fKdjDOdpyq4MpiHhJ2YTAY4ODgICrfOmHzxJpyz2tC1WVReZb9gLpLfl1XHW7GC3mQBfVDuOkyBoMB9vf3o/JtAsJNzcB3DlC4SdBmdLtdDIfD6Z+sVglWBckrVoB/Yqyydz5lS07VCcFXZug8xK1bt7zvtG8aVGnUdaF1moci4zJAVgmmwmJDqC+uWnu0CTHcNK/tgk0qqyFuun379kK5icoj3NQshJvai9S2rfJMmnyeqav8+uUU29vbODw8xHA4XIuxue5YJ+6pdMaKu71lHg0WUjBsHlNuep8CduXKFdy5c2cmXii/JuHzSs57642Ww7eSsGqg/Wtd6t02xPBBm2SJNcJ89bh69Spu3749Ey8kT5MQbloshJvaiyaNp0UaZnSsn5+foygKdDodHBwc4ObNm1Dq5Xer6ipXIFgkorcCpmyZCcWj8TlbbWK9zDRNrNfZ5n3W5xc6nQ5b7nWBbcVyVdqFs5WmynYzQRpi2le4abFYpBzCTcJNgvkjyzL0ej10u130ej0URTHzLSvpj4JVQdIZK3rtOndAPYI+xHjPqioJqYoL/X18fIzDw8MkGZpAinLWlBx0u8+qESVHsdZbctqm1K4qbNxk+w+sPjednJxMz322DSkrc02UL9wk3CSYH/TLKvRfURTI83y6NbAt+otAUBWVtgIC9m+6uBQcV/rYidZWRmjLIZ1MfeDGOTo6mi5ft2GCMttyngfDXXLo38DqEWXIM+zygi/yuawb6NazdeKmw8PDVnGTLn/R3GRCuEkgaA42Hjw/P8dwOMT+/j42NjZmwmxpBIJlRO3vt0zdsx47oFJejsEtg75FyYbhcIgbN25cSjNPuBR3V9g84Fq9XCVw6uXqQyle8lXzqC8Kwk3zg4+bFoV1MCiqcJMZzsWin6mg3aA7ms7Pz9Hr9bC7u4ter3dplYrGFwiWEY18OKBJpSW1jJizYSHkeY5v/MZvtBIDML/JhrOffp6g20za4pmuGxzlPDQ5zGOMCC5DuGnx3LQIUG5aVQg3CdqKbreLnZ0dbG9vY2dnB/1+fyZ8lcelYL1Qi2FFiTo0ecV6M21ejarbc3yyhFZcTk5OpmerbFsr2jDZLEIG87mbisyiz1Q0hVCd6Bu4bNtmBc3Cxk2h+LH507IWyU1XrlwRbrLAxk369yqOQ+EmQduQ5zkGgwEGgwH6/f6l71a1gZsEgjpQ6+vWQ1ttzD32GhwPG1U+OCshVLF3pUuZUPb392deDdqWlRmqJCxCLt8zbks7zQuu8WBTeNetbZqGj5tc8Zedm5RSwk0MGYSbhJvahDq3CrcJeqxlWYaiKLC7u4vbt2+j3+/j/PxcPggsWFkkG1Y24o0h4Nj4sQpMFW+tL97Vq1dx7949dvx5YtHnq0wZTDnWaWIO9UlbnHVqn3nAtSIu3LR4tI2b1g2uegs3NY91bUelFLIsQ7fbxebmZtABLxAsO5IMK6p4hLan2NLGxAcuT4rcrQqxRh8tjyo/BwcHU08L1/O8KCx6kjTbsE3t0iRcXmAf1qFd5gVXX5snN3FRNzeZq1VVjcu6QceAcNNiYHMC+LAu7SKoH9qYyvMcnU5n+u0qbWhJ3xKsKpLPWKUMCvOFBjGGkbldI6Z8Wzpumbb8r1y5gvv371+6n+J5bhKhM27zQIpiuyqwbTnS9wH/81n0c1sFLJqbuGlT08Vwk4k2jEHhpsXCxU0abXg+gtWANqD0R4D1Kr02qtZp3AnWC0krVikevpQ92zFKjqs87pmH0HkGpRTu3bs3c36Bm/+84Wq3RZ1pMNHG9qoCTt8J3af5ubBqbdcEUrkpZfWoTdx0//594aaKcrSxvaoghpuqYtXaTlAdZTn5btX29jZu3bqFjY0NFMVLlVP6i2BVkbwVMDa+ufXCld62FSdl8Lm29LjI37W1xsTGxga2trYupWsrOSx6HzPdHmnKtUqIUVSoQp26FUzgRhWjype+aW5ygcNNm5ub2NzcvJSurf1FuGk+EG4StAH6bYDdblf6iWAt0Mh3rDSoMhC7zcZ37SrHjBuz5SO0ne9Tn/oUBoPBpTht2jZB23jRy+1UaVkXUrX1O5viYl5rtKk/LRti+pdtvC+Sm3zwcZNSCp/85CeFmyKxrtxEIdwkaApFUcx8t8rkKIFglRFtWHG2rthgnmHglBFjFNm8jrETJT3zYOa5LKtVKedE5gFfn1hF2PoGPbtgU1jM/6vaNm3EsnGTiWXjJv27LVg3bgLsxr5wk6BuKPXyjJW8rEKwTlCR2wW+BuAfNieOQCBYAK6XZXm4aCGqQLhJIFhJLD03AcJPAsGKwspPUYaVQCAQCAQCgUAgEAguo9EzVgKBQCAQCAQCgUCwDhDDSiAQCAQCgUAgEAgqQgwrgUAgEAgEAoFAIKgIMawEAoFAIBAIBAKBoCLEsBIIBAKBQCAQCASCihDDSiAQCAQCgUAgEAgqQgwrgUAgEAgEAoFAIKgIMawEAoFAIBAIBAKBoCLEsBIIBAKBQCAQCASCivj/A4LBkvNRTAVyAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 1080x360 with 3 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from matplotlib import gridspec\n",
"\n",
"from skimage.data import chelsea, hubble_deep_field\n",
"from skimage.metrics import mean_squared_error as mse\n",
"from skimage.metrics import peak_signal_noise_ratio as psnr\n",
"from skimage.restoration import (calibrate_denoiser,\n",
" denoise_wavelet,\n",
" denoise_tv_chambolle, denoise_nl_means,\n",
" estimate_sigma)\n",
"from skimage.util import img_as_float, random_noise\n",
"from skimage.color import rgb2gray\n",
"from functools import partial\n",
"from skimage.restoration.j_invariant import _invariant_denoise\n",
"\n",
"_denoise_wavelet = partial(denoise_wavelet, rescale_sigma=True)\n",
"\n",
"noisy = img_as_float(cimgarray**(1/2.2))\n",
"\n",
"# Parameters to test when calibrating the denoising algorithm\n",
"parameter_ranges = {'sigma': np.arange(0.1, 0.4, 0.02),\n",
" 'wavelet': ['db1', 'db2'],\n",
" 'convert2ycbcr': [True, False],\n",
" 'channel_axis': [-1]}\n",
"\n",
"# Denoised image using default parameters of `denoise_wavelet`\n",
"default_output = denoise_wavelet(noisy, channel_axis=-1, rescale_sigma=True)\n",
"\n",
"# Calibrate denoiser\n",
"calibrated_denoiser = calibrate_denoiser(noisy,\n",
" _denoise_wavelet,\n",
" denoise_parameters=parameter_ranges\n",
" )\n",
"\n",
"# Denoised image using calibrated denoiser\n",
"calibrated_output = calibrated_denoiser(noisy)\n",
"\n",
"fig, axes = plt.subplots(1, 3, sharex=True, sharey=True, figsize=(15, 5))\n",
"\n",
"for ax, img, title in zip(axes,\n",
" [noisy, default_output, calibrated_output],\n",
" ['Noisy Image', 'Denoised (Default)', 'Denoised (Calibrated)']):\n",
" ax.imshow(img)\n",
" ax.set_title(title)\n",
" ax.set_yticks([])\n",
" ax.set_xticks([])"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.4"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment