{ "cells": [ { "cell_type": "markdown", "id": "fc865825-fe1d-4da8-83e6-8996422cc66d", "metadata": {}, "source": [ "# Deconvolution on CUDA-compatible GPUs\n", "We will now apply [pycudadecon](https://github.com/tlambert03/pycudadecon) for 3D image deconvolution." ] }, { "cell_type": "code", "execution_count": 1, "id": "f9023ea3-8166-4484-818d-f350ae2400c8", "metadata": {}, "outputs": [], "source": [ "from skimage.io import imread\n", "from pyclesperanto_prototype import imshow\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "from pycudadecon import decon" ] }, { "cell_type": "markdown", "id": "3541b287-b03c-47bc-944d-7a856b906035", "metadata": {}, "source": [ "We will load an image showing fluorescent intensity along lines. This 3D image was taken with a confocal microscope." ] }, { "cell_type": "code", "execution_count": 10, "id": "8473a58b-6977-4d72-a107-b3c5dfb56460", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(88, 320, 256)" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "image = imread('C:/structure/code/pycudadecon/test/test_data/im_raw.tif').astype(float)\n", "image.shape" ] }, { "cell_type": "code", "execution_count": 11, "id": "b12e8550-c366-471f-b060-989331db129a", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARIAAAD8CAYAAAC7DitlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAD2XUlEQVR4nOz9aaxlS5bfh/1WxB7OcKfMl2+quZssUuqmQFqiSEmGBQq0RX6w3TIECaRhSzZo0zA0GIY+iPQXGgYaoGGZhmxBgtuwTBIwSbUECyRkWZxMQTYsNdnU5J6KrK7qrnpVb8zMO51hT7H8YUXEjnPuzfde1bvlyq7OeMh37z3DPvvsHbFirf/6r/8SVeXVeDVejVfjswz3wz6BV+PVeDV+449XhuTVeDVejc88XhmSV+PVeDU+83hlSF6NV+PV+MzjlSF5NV6NV+Mzj1eG5NV4NV6Nzzx+YIZERP6giHxNRL4uIn/sB/U5r8ar8Wr88If8IHgkIuKBvwP8N4B3gL8F/GFV/aUH/7BX49V4NX7o4wflkfwe4Ouq+g1V7YG/APzUD+izXo1X49X4IY/qB3TczwPfLv5+B/i9L3pxI60uZA3JOZL0Py0eE/gk7ym+RUQ4fKWmg8ZD2bE0/V4eXwBx9h7nIUzgPUwhPicQNJ5LPL8XnVc6tnf2Hufm8wzBfp8CiHDXM0znXJ57/D2d530f+2muUzqH79cZLd+br9sLPvcFnyPOoekalNce7s6D78Vr/pjvf/CZgPj4d/ld0j2977yPjn3D849U9fVPf3J3xx/4x9b69Nn0qV77t//L7i+r6h/8LJ/3gxo/KEMi9zx2cFtE5I8CfxRgwYrf6/5x0DSxnP0uDnGCTvFCO/sbQMN8OHGS/3aLlrDvEO/jp9oxpW3RfkDqCp0mpKpgmuyYVbwMzoyErFf296KFYTRD0PWoKnq7MYMQwryA4nEIAeraJui+Q05PYejte6yW6H4Pw2jfxzl7X3y/qtpne5eNipQGJn1fDWhQ+37x9zsXP10PDcdPgAbccknY7eJise8gVY2OQzSi83U7MBJS3NbyDst87MM7Ht9zNBvE+3gNjl5z3+/Hn1U+fnxOxfmk+6/jmO+vjqPd43QtQrz+rvi9+Dyp/Dz3iuuXxl8L//av3/3w72189Gzi5/7yFz7Va+u3f/XJZ/28H9T4QRmSd4AvFn9/Afhu+QJV/RngZwDO5LGawZifF+/RoIc3UgMa3IHBcU1N6IdsYGyBuINFJk4Iu302IgS1xTuFg4UiVQVO0H0HVQVPn9vneA9NbUYEYBjsZzQGqmqejffovkOdmBG4uUnfFZeMRTRiyZipqk1wbIeE6CXVtf0+DOgU8nM6hmxEzNDatdBpsr99MfnlaIFoQKr64Bql53Qc8jVPx87PiwPC4bHtix0ce755h+892AzU7mk+lhQLm/IY9ngycNnQOX/wWeJ9NBA+fz+dJjt2vKd2resDo2PfdcrzQ/XIOB15Nem46b022XiAoUzHBvg34PhBGZK/BXxVRH4M+A7wh4D/7se9QYPOk+Ueg5FGNjhxgoauyzdZnMw3Ob4WAO+R9HDQ2aA4mb2RtPv3oxmO3c4WcF3BONq/OJHy4yHMnkUyFNGI5O+l5j2ErsOdntjkPj1B9p15PDe3iHfIyYlN4CrupLdb6DpoGpvUpXHQYOcY1L5bWuzOvJvkkYgTlNkzs0U33Ln2efFFI/IiY3RgAOI1K/92y6Udp+/LG3uwQRy/B1WkbvJ5pefyXEiHSc/H+59eNx8nHMyTfE/TOSdDme5hYZR0HAuDFq+lTtnz02k6OpeRez2h72MoEL7vGPPlGT8QsFVVR+CfB/4y8MvAz6rqL37i++JNl8p247tudphfkxZLmnhpNz0e6Rghhg6LFkRwZ2c5XJCmMcMAc6jiPdLU83mlieMcUldI2+ZwJoUgItGIOAtPzJOYz48p2OQdJ/RkhdYVsljYOVae8PgUxriTNrV5QapmsOI/8S4ubs1ueg5/olcC2PeJHgrMxiJf22IXTMY7PabFcdCA1M3BPcqX9shDCbsd2nX58/IxygWuR16myOEiLRb6fYYvP1+eRzo/DejQczzye47nRzSaOH94nsmbivMtP6f64nDqM4zwKf97mccPyiNBVf994N//9G9IF+po8qQdOGEdpYvpZrf6YNdL70sjGgaCWtjiBLZb80ySp+H9vICDEnZ73HoJVYVMwX6mkCJ5PTm+Duk7m7FyttClqc1DiGEPTpBFiy5b8A5ta3TZWNhTV2jlYBjQfYCus+MNQzZSOsVrFLGc0pgcnFdwB0YlPy/ujpEor/GxYcjhZVzU9pYZe9GgM96S31ScT5jM2yhCggMjcnyPj0bpAR17MeViPpwvMs+VpkGahrDZmld1z2I8Dp3vG/7xI6anz3KY9JBDUYYfgdDm5WG2FvG9eG9/V/WMBWTwbLhzw/OuG0OVPNkLA0QMZZL7jyo6jPNB4q6uw2gLt6mhbS30qSrzZCrzIKSuDYSFIpYvbHLyHI4nnciMs0wB2Xa457cGCO573LY3T2i1iMYiZOORwyZfGIFkRIIeGg0nB3+XoU4K/9I1LseBEUm4RmE07ix4DTOwXXgf/vwsho7+4H5lgxK9IklhW8ZkEqA8e4qzhzLfK6kqM/pFGHvf0L4n3N4WHsbHhBDZ23B3jjs9e34QJj3kUGBCP9W/l3n8wDyS73lEXKT0OMoswsFEzj8LbCS5+VMBRsKMHSS3Nh4nYyRpZ0tpwYQ51BX0A6rmjYRnz22ipQxAmnjxd5kmpGlw6xU6jmjXH3gTDAM6DJbVeX5t37froarQ3Q76AW5uCBFYTUZOFq2Bu4XnA7PhUo2YT98jMmcnEpBLxD6kMC7pmuXrU3ojJVCdnkthkXdoWs/l+0qPQAPT1fULvQxgDsn6vsAkZuA3A8CFoZe2RYcxGzipGwtjZMbFDsOQezJWhUdyBwdKz4UJVTnAaKRp7FyOQ50HGq8wkh/A0DFmYApc5OBncQPLG20v0byL5oxN8gyyJxK9haN4OeEMOppHosMMrhpe0iDrtR3He1vIIhG70Ay6hptbdLuzRZJSy87NHJIQLJszBbTv0a7L5y8nJ+YJVXPIJd7FsCyGLvE4qppDqRxaTQk4dObJJKynfsF+kTy2ZFjL6xx/z4ZdAyHhH4Whzkak2O3z86UHUHgc2ve45XJ+X3lvkyeaQ9r5PTMw6uImU4Q30cu0w0i+pnmEF4QkpUcE2XPJWIuqGekU+h7RCj7rUGBS/VT/XubxUhmS5HKn3VCqGte28cljRN7l14d9N++6xbFgNhA5syEFT6NpbPFNwSZeWviJ5wFoP8BuH72aCVku0d3Odv2mybyRDIgWKc3kgpc8GNQyLdp1EViuoKnNSHhnBuvsdD6f5dLA3bTwSq9EBLdc2Lmk55MHVl6zAhcpvZHjkCQvkgTY3tm1yce4A6QeL2qYQUzIHkcKazKuUgKX0bDN4O89wOY9QGeZRcr3uAxby+fyR8U5UAC1abjV6tA70ogTpbD5gUf4lP9e5vHSGBK3WqHjMN+oOKESRiB1dRDC2Ev0YGEc3OQyRSpujrGL9zNNuLMT3HqJeId7dBEX7cImzThmAFN3e/u32Uam62QhSWkkAOraFnbiN4QYaqWsSzRU0tQ57qaPKe/bDXhP+PDpnGqOnJX8We7wloXdPvNXDr5fwlCKcAWYQxxxBzybvCsfG4x0LMgYQnmc8n4cvF/cvZ7Qfenng5GOf8T7yJ+fzqf0Dsr3lIboePMpPJ85FD4GpfUQPAZwHte2M2Z0dJ0+y9BPiY+87BjJS2NItB8O4u78eOJFREZoZnWWo5gYBxM7HLnNap6HSAQtI91d+wHaFtrG0rHRiKSsi6asTsy+SNPMGATgLs7tuUQiSwai9B6axhZ832fvQRYtslzM5weEq+uM20hTG38lqHlPOocxOeXsHe5kPX9OpuX7/DOHR+naRNBZnBRp7yK8Ka7rwZ91c7gAS8wgXfsyrdzfTcWWx70XuEyhVmFwZoNReFrO3zF0d4zUvTgJh+ddhHcwY0/ZaBV8pdKQHcyzzzBUYfiU/17m8XKArQfEoZB3CdfY7p49jQQEHs33Mj15MJzAARnKzQBo8ipi3E/XEW5ubSIlL6QfssGwxeiQNi7ozryLlEIW73O2R6dpLhtJvzsBdchqZceKnomerpFdh64WyPUtvPYInl3GkOjuojxIBbtIotvtM7U+P18Arsd1PAm0tEsvqPeIRHp+ut4HuMFRWPQCMPPeVKoq/uyM6fo6hwnJUT/wIEuvKc6DFNocG7gDpux9owiJDti6wTJJuYTgiPWbj+8PjZSOcxiUjNXDhTjCdG9FyW+s8XIYkgh46QSubTOoF/oXuKjFsEkRQTp/NNnvAxBx5kF0XQ55JBXieQ91ZR7KNB2kdHUYDbOADKDqcGOPR4MjIvZeyIsmYzXx/Qw92qllhUKk5XuH7Ds0hJmWnzI10SjYIefam1yTExRk9lKA6IX4/P1zmJafm4+VjuGWC8J2e8doHyxEsMyaE1LG7BgMP6jZibjIFEsF7lLsiwzcMUEuXWe4YzQ+kVl69Fw2TJlkVryuLJEoGbYRyzlkwJbA8sOBrT8A2OX/7+PlMCSQJ1ToD2s+ctybWJlaMFrHwYwI82STuEAP08RuTvU6h15dZ2BVRGC5sEUNBqjGsEfOTtGuRy5O7dj7Hj1dobVnvFjgujdQ76jev5oJbSHY64Yhejx9NDC9/T2QMRYAruYanDljMB7gC6VBE+dyFkEAKokZpRV6fZNBW/rBjFWRfUpYjZbeTrw+0+3mHgzqyMvIqWGHfZG7ozQiUtWW/XAeOGaJhjvHPZ4Pd+p+8nNHNTMpY3NsQDKVIOIg0SMpPZ8Ubh54PgkgL19fejkPDLq+8kh+ECNN1PRnMflKRmtJDpoR+MoWYSKoeTfvPtNkdSsFbpFDlq6DurFK3TRB2hbd742/UBm9fXrj3IzIqmZcefRRTfBC2z6mebplWjXIFPCXW2QXuSO7PepcTt1mbCN6CQd1K4nWHtm1OUyJ2IicrKHr7JzSYxFj0WFAHp0XoGsCly31bMV/8fqJoI6DiuKymrhkDJfeYPb+7ln0B3ySI5wj8VxK72M2AHezMxmoTgu3KOgsPdY8P1Tnor2yniYZg1yLVPCLorFwTWMwZullFOzYO+eneuChfdZhhLRXhuRhRklrrub6lhIcK7MFLyqdh0MAL2ddIliaK38Tu7PcwbrOdvC2KNKLHBS5NjaqCwE9WcK6pr4dcUNAKzvGdNJSPd8aY/XGsi86Rs8iZlUyKJqyN/2AVB6WCzQVATY1YVHhbrvoIYUYegUzCvUj1AvSjTbBQ0CGEW0b2OyMgr+PYVtfcCFySFRcB7tw8zXwHsI4G+3ja+49TP3BdS4XbbmI5+JKyedhz5MNVvZc2jafT87cRaKZNM1BtXLmsjTN/P3S/Mk8FDdnbo6qhdMcy9yYdOxsRNwdg5QeP2b7PsRQYNCXJufxfY+Xw5CokjQtDuNrd7Ab2Q30dxDznEkoQMR8jJRxiQbErZfGGymxh7ZF1bAO9p3t4N7BBDr1hmFUFRIUNnsaVaSPO/44If2ArhYwjMho1a+ZEJUMm3PoGBfZojVj1jaEkwXTusXfdHRvrXGj0j2uaa4WDKeeahuoNiOuG5H9CJUjNB5WDSpC9WzDdG68B1k2yK6382lqRI0sl8ltiZND8tbC/HvCSprars99WYlj71DcQa3NjD24A08y37sSC/m47E5Rbay73b01OuV7Mh+mwDEOvNXSY4qckMODFV5HxoJKJvURtvKAQxGmlyd5+n2Pl8OQHI8iLs4l44kPksrjS2Gj6JYfxPfHdRFJr2QTOQK+CJ92u3nBZ1p5iC75gDs/s12+6030SASubhDnrEL3ZMW0bvA3e6aLNa6uUBF7nXfIZmdaPVfXlgnKbFXF3eyQbcf0+IT6dkBFaJ9B8nZV4r/KwaLCXW3xt8r4xhmh8UyPVhAUf7Wz7yRi5y5C6Id8HXKGaRjxjy7Q7Q76fs70pPqVlE3qh8PFm0LOEncqd/QjRqzdmPuwkLJkv9j9U5ZpKkKRtKiTd1SGKcXv5iUcHvcgK1N6XWU2J4ZFJcg6p31jSNZ1Gbs6+PmAI+jDHu+HMV4uQ3LsUhZpx+ydFGJHJV15psTPsXJ6jlBM/vRRJeM1aYrEx1PIo84mabi6LmLkAPu9TajaGKmy2VHtOrTyyBhMEmDR4LoBrU0USQZjn8pymY0VbYM6+3wJSqg9w2nNuLK/Aca1R71QbSfqyz26XiDdQPXRrYU0q4WFN/1gn98PZgRiAaM0dax4LrgbyRsrvjequEcX6H5vi8cJonLXs6AEMe+7h7OY0HEIkGntR5maZDzuTR+XuEpa8MlLiN9Hh2mWxSzxjfI1yVDkDWoGaHUcLFRKRiPNqbIWiESMbB40vHmFkTz0SEajjJ0LTyS9xjUF1+Hg7Uf8hmQ4slfiD/Jss9xf/NxkWNKOnDIwcScHIikuwFCorWU8xUNV4foBvbk1Z1UEt1gYpX6ajOwWZQS0tgXdP1kxLRzj0jEuhMXziWobqK97tHZUl3umVYMbJvNKJkVPFsgwwUaRXSRKjZN5JN7D64+R61u7RrE2KJ0/mGeUsjdp6DASnj6bHzv2+NLrpsmuZXHPDo30/fKPB+M4E3PP76XI1XyOc5HeQWFn4obo/Jxr2wxmp5DmwMjE+5PnwovIc6TQ2d891oMMYfoRwEhenm/wovRf+ZJxMHf9HtS8fG1Js08ciZwSLitPj6UcIRfticisfDbN5fzZc0k7/25vk38Y0dsN4YOP5tBgtzdvZhhnozZOhJXV0UyrhlA7ujPDQkIlEKA79/jbjvqDW9vhp0CoPbu3Vkzrmqn1OXSaHp+YEQFks0OfX5kRSfqkpfc2TXM2pK4PwWZAp2C1TUXaNYPfwbIVx97IQYan9AQjLpEB0MRtSc+VdVGl7EGZ8SkA4uRJ2XlO+NPTo3OfDrCZ0HWHnuvHLPyDWp38oObPv5MCf8DQxnJF7lP9e5nHy+ORwMHOdpwVyKnHOJK3kuL+/P742lnkRw0XcQJhLsqTJlLTm9pAVpjDACg8GYwclibzcmGiQ5GEhptBWyOcxfqgovQ/pxf3ezNs7wyIM+Zu/Y6yXJhhWf9iB0FZV94KBesali0uAF6oF5762ZbQ1rhth/QD/uloHk/Xz1oldT2nhb1Dnl8ZUS5yT9ybr6NO0Pc/wrWNyRlMk3klu71dx5jV0indkyKlWtynZFzSc8dh6By+pBtn+EaWgCzvfXr+KDS5jwE73dzcCZGOX3NwjkeA7bGUI+keHUkZZHGn6Uj+8aFqbVTo9WE1Tn4Y4+UxJMXkEyc5e5NxkZhRKIWepdAoLfGUREqTSLEnJHFmmUOZGNbIOFraNxmsRZs5HKm4Tp3DnZ3O2Y9Yj4M4EyHad4hzVtCX3p94I5UBdsnD0SkY8FpVBtaKwPOr6LHETMo42vG1t2RW10M/0N4YUOz3/YxtdL0JLUVZAvfm66aj0tTgHePFkqptcM+uCNc3lj2qPOF0gbtaHN6CnL2J2Ek0lAByAAj6gwV/oD6fF3ENeleVXuoKphkPyWnchK2kRZ/FqWzPtkV8ZDQKTkhO85aGqCCUZdA+emTHtTxmrMpjH9HzU7YnPi7H5RefYYQfAYzk5fGXEhmqrg5xkfL5mJZMhoU8AeJLklEpwo/U4gGYGaT2R3xPsNAjkcZSzYxIVkcDDIgLagu2622xgy1k5/LrKHe+uHsn7VeqKlPik1ejm20+lywJ6CSLTeu+QzdbdJzstdsd+uy5VQp3vWma7MzjkqaB7Y7w+JTpfGkhlAiy2dlnLxe48zOmixVu2yMn1nZDFosDzCTpwuoUSIr79oQ7+JnDnXuKKHP9U7nYI7el1A7JurnHIGbmFVUHf5efcaeY775x531Fxiid633g6T3Znnv1XT/jMLDVfap/nzRE5N8UkQ9E5BeKxx6LyF8Vkb8bfz4qnvvjsaXu10TkDxSP/wMi8v+Nz/3v5TgGvme8VB4JUJCMCrc0gl05bveYYUngYVDEzelfm9zxuM7cc0lckJRqTHyJUuRoMO8kPaYJSBUx4wFk/dUIZOoQMxF9j5yeoJtNZrCyjLKMgJ6ucjOs8fEav+nYv33C4v0t/eMl1dYyLW43ggP/7NbOqR+sH85kBDkdR6gbZL00T6hpCK9fINuO7suP8dsR14+ENmrNJoO5aGBpLF3XGyFNmxpZL9Hd3oxQxIRSTREwY0NFtkPu88QLcPy4+NLuw7GxifcqejNlyJDTsqWHcpwlUiumSwxh8fewbgujcafQr8zw3EfRj59RhjsHodonySF86vGgYOufBv414M8Wj/0x4K+r6p+MPbj/GPAvi8hPYN0dfhL4HPDXROS3qeoE/BtYz6n/BNNd/oPA/+PjPvil8kjyLnS0w5UpXwtJCtc41oeUXkwuZy+8EyOEzZR0WSxmjyNpszox0LTkUKSCLomU93Gc/93cmkFJwOvWQg9pGzMiiXeQtEdWLeGkZThvCIuaaeGY1g3T0jEtK26/tCSsaqSb0JtbtPLmVVVV4eYzn9dqCYuWsGqYHp8wnHj2b7bc/PgJm8+1jCc146oinKyMKBfxIOlHZLtHrm4MLE7XpK5mtbnEho14SV5wRztxee2Td/jCRXaUgs/3uljEM6FwZpGW4cWLPIFMh4dYk3XoKc0AfOKJzCJL9xmRGSieqfGzYXm4ZfOQYKuq/kfAs6OHfwr4M/H3PwP8E8Xjf0FVO1X9JvB14PeIyNvAmar+x2ou6p8t3vPC8VJ5JGnXuCNQFCdwCa4e9GIpXwuzx5CMSyGSnKp9w+XVTJ9PGYzlYmaAJmm+KAZNgXkAhoEkRbS2tQXQdcj5mS2CZUto7LXD4yUyBHZvtvhO0QrG04bmesTtB9qnAf98S/20Qt5/Rm45enkD42gkuBh66TBYWHRLxgCqcYJh4OT5mvGNM1w/4a626KqFMeCubo1rsmjQVYu72cE4mmzCojWPpG0zMS/JSKZK5iQ1QPIOj5ir4o5Svsds1iPwM78meTgF3V2axs7jvvTwPenbTGCLfwM5JXy/MFNBWjzqdZOfP2bmps8WZ933Utj2QEz56QdLSHtTVd8FUNV3ReSN+PjnMY8jjXfiY0P8/fjxjx0vhyEpXccjDklJdTby0XiwiyXl+BCFkUz/I06CZCgm5vauqQI2gp44l9XKpK5twVZW96L9gKybXHwny0VmpkoATtaGjyT5RLBU7BQpRiJoUxFqR7WfcKPi9xP712pcp9TbkfGsxe9GpvMl1QdXprC22c7UfkBvb+2cxsnCnBDJWf0AobOK5dMTAKr3r9CTJdrUuKfX1vrCCeAIiwYqB89uZh3ZCCgDJiPZdbOyW1Pb31H7ZA4t3cwcLu5DWZUdH5xDiiPQ9cBI5PeNBxm4g8K+dP+Sx3QP7+NYRyQ/XkoC5AdfEM6k44wz4/aApDaGA+/rsw5FGPRTL8MnIvLzxd8/EztWfj/jPuulH/P4x46Xw5Akt7FI/SZSUYi6IQcvL/kBWkgPlDFw8kqiR5KzNsmINHVO15K0RDZbZGXMU90a5Vz3kY+RetWkEEgUTtfWj6apGM8WuG40WYGrHaGt6d5Y0jzvaN/fIPuBkysLfdr3Em4A7mZnjNRILAvXN1lgSMCU2yBXJmfPKYK4qaWDbnb2vmEwoPhmg56fItt9VrD3HxjrlnG0ymbv7Pi7vXHMNFg5QD+ngw/S4PFa5nR3wfmYBZ+lWIDcu7DzsQpmq3hvXlFkl+ZWnTHEyIJCJYZ2zCV6Adv2jgGJ53NcSVx6NPfiMvmAD+SKMIOtn3J8pKq/+3v8iPdF5O3ojbwNfBAff1Fb3Xfi78ePf+x4OQyJcNiGIhqAZERKTkJK7SbOQin9d0BsKpmZiQ+RaONJEjFhJLHtRNjtcRGDsYMKYbs1vQ/vLI17ehLDH0sLE9R+qtI9WeD3AfUrplVFfT0wXLS03x2sRWdQC1UiXiPOmRtfN2a0ut7Aw77PUgLJOFA3hr2IoCdL5HaHLhqm8xXuamvne7u173azse/77AoN09yDZ4EZy+XSUtD7/Zx98g76yT4rmPxk9g6czNczZWpSylyLHbrcqY/SsfeKGhU/dSy4I+JmHZMyzODY6yjo+MeGJaWEy89TxfLAMs+vF7BUy888NjgfWyLwPQ5FftChzV8C/lngT8aff7F4/M+JyJ/CwNavAn9TVScRuRGRfwj4OeCfAf4Pn/Qhn8mQiMivATdYRn1U1d8tIo+Bfwv4CvBrwD+tqs8/9kDKgYs8f0AhKZBGAbLml5VaGrXMzMzolZSsVIW5luaoX0yelBPR7Q+z+JFELsl+byHAMCD9wPTaKdOyRoKyfOeG6aRlWla4IeA3HX43ICGg6yV88NQyN1F4SJN40dCjA3N2KbnvYwf7vaVn69rU2yrDasJrZ0azHwPhYo0ME92XHuH6gIwBbRz+tkeGCffsxhZM21go5gyotZ45I1pFRfyqQochy0+65cJ6K3sPEtm9yZhPRb+cdD9i6GknH4qw9KjDX3w+e6DJKCWqeoGHzAbG5Xuc7rmWc+O+9O19noMI7uQkli1Q4CszjpPDM3HmVRXNuey+RE/pwXgkDwPeisifB34fFgK9A/wJzID8rIj8EeBbwD8FoKq/KCI/C/wSMAL/XMzYAPxPsQzQEsvWfGzGBh7GI/nHVPWj4u97000fewQpMJBpysQzmCt4D1K7ZV1NMcEkhTPpsFU1Sx+mFHIaRZ+Z/HeZFUlGxvsslqxdH7krIVfLEqC63iNjILRmUPx+wm96ZD8wPV7jdz1UHjlZmyF4fmUkODFRaxOGdogE85BSGNPUZD7H1FmKOS5+mhrpnWEwAsNrS8aVRxrDY8aFt1SyiBmgm1tj1iZCXQyRUtVtykgRNWuthCB6C30/h3SJrOZjY/bSK4z4VaarJy/S+xeX4IcJ7eMC5h726LE3o+FQiyTd6/jYvYzVo/RvuLnJBuQAh0lzaTz0lpB61jg5MmifdajyYOlfVf3DL3jq97/g9T8N/PQ9j/888Du+l8/+QYQ2P4VZRbB003/IJxmSNHQGuPLkPOq2V2YN7vNK7I+Y4eiLCZVdbplZqyn+P+piNxuQCCz2AzBkvECCGvUc8Je32bNwQOrAJ5s94+tnDOcN+rilvh6oam8AaC5Hj96HhugFiRmRdJ5TAILpioB9p+sBWS2REJienCHdhO9HhouWxYd73LYnrBqady7tPAZj9vLGa7P84rVxXXSzNbykH+zz9p2lgEWQ6M4nMeyDnV/cIRhcLqqjquGZNVq8l3nh59J8ot6H3M3CHWRVYG6xEbkgZehxb+pZj7REyszNMbekfD7R+Y8M20MOA1tfUeQV+CsiosD/MSLIL0o3ffxREq3aCRq191Of2sxkJRmTcPj7kXdy0LumSP3mWpRkOJLqmFhzb9rWQo6um4FYafLzmvrcajA8wXtkt0eGCj0/MQZp75HRSF3Vh9f4mxa8GM3+epsNCF5IeibmmVhaF0DWa/Tmdj7fyOJFBJoKPVvDsyv8e8/jYzWLbzy19w4j7oPnlv3Zd6gGuDgjZZZSgR9g39mvoDPvQ87P0JuNnYePvVyiCtvh/Qq41Yqw3eb7lI38vSzYu2FNTuGn6t17Sv3vZE2Oa29KrONoJE/jfjJbQakvx5Ge6x1tlPvo+Q8wXgkbwX9VVb8bjcVfFZFf+bRvFJE/irHnWLCyFJ2Tueetv+dGHz1eAq1u0VpD6zSZ082PbMfcViJWzWaBouS1eG9AZeSJZMPRWvo30+hd5FeIGF09nddmi2qwsCEej37A3W7jOQbTHkk79epkxkg08mNO10Ycu92YmHMIBsgOQ2TFGrbCh89ANWY4HNzGhZmuQQplAFmt4PkVUtcGrlaVnUtqKdrUFnJtOvueo4VNSbSaaZrlHuPxQtcRdhErSmpwsZp6xhnIIcOdJlquilR5Zm8zKZsVolVWN9W80Ms4+FkeP4Ueere6O06+O3NKqupQC7YItWcV+jnj81A8EkVeCRup6nfjzw9E5N8Ffg8vTjcdv/dngJ8BOJPHeqdI7+NSb94bVyTtXE7m/qzpNWmCpbjeuUiB97YjDePcBCvWk6SMBU1t/1IPnNVq1vBIPXjBetxE6rwR26p8DGkbc4uLnrNSNbCOHoCIGY26Qhc1Kua1hGVFdbkg1N7kHG+2iCrKOH9OMkhB0anAClK4FnSuRI64SObB7PZzYy617Exi4MpqFTv6RdB3sTBlesieWzbGTTXXwaTFFsbD1hdwCMDai2G86+HkX9MmcEDJP2yJccwrSff+Tq1OMn4plXz0syzMu5eCHww3Oxa1zozYBxq/qT0SEVkDTlVv4u//OPC/4sXppo87WDFxOHRtS1wktp8QUsbgHtVzivfH3RcMMMzhgepMl++J5KvGjIWPhiZVDUcgUs7PshgzlUe3exNuLtKMut8bflFbhzypi3oeILx+YTF9xFG0tgK+sGoY1zWhEtygTG+d4ncjXgRkjVzdWgp7GOZKZXHgAuKqOfRJYV00jHKytlAthXDDmI0HISCrFW5dm/dSFtKBeWZVhXvtcfbMdL+P4VEdPTTj5Vg17TiHYsx4VvYwj7I65b0t5SHueBjH+EZ6by7M7O8PX4p5kDkox5qtia2aRZ11xmxelPUpSGkPMRQIPwLCRp/FI3kT+Hejy1sBf05V/wMR+Vvck276uCEidyeUFOxJOADqskcCdzM4ZfVvYmQCIne/amK55naYqSF422acIFf/9j2cn5qUYV1l+p80dfYwaGrDPeI5WvHeCJH/Id1I//YZzS5iEsOIVh7/bIO/qdDKoc4hsXDOXW8tlHEO3e5sMQ8F+JzS0yur7dGbW7g4Q9//yM5Ldab9Jzc9tRqdJnQbZQ9S8eNuZ98jSjQigp6tjRbuHZxZV0DAPJXUMyel2UkdCWePMRl701Ad53tbhC+HN2UGNF9kIO6IUU3TgTL8/MSLcZY7x7tHi+Te1x8zdD/z+E3eaU9VvwH8znsef8oL0k0vPFYImVmZUX0n5G5ucRwIFpV0eXtyPmCZBk6krkL1rBQ7pjZshn0/N5cKalhH8mDqivGLT6g+vDY+SMJLVM2IpPNbLq3K1nvzWuoZjZfNHq0rmg9uCWdL3OXGUsH9YFW6O0vvutr6DbNoDRiNim00tbFtoxdD9AA0poT1+sbwkKBm/GLfYnHOQq7JvJKEjaTz136wMCyGe3J2mm6k/fCCqFpV9L6fjc5mOyvlF+GWTgnbmLkmL7jp8+9lNqRsq/kisLQY6fn7cJSUzZG6yTKN+bvlxlczEa5UxM/qbFl/Vu/8/jAYCa+yNg828vbuDj2MY1A16IyNQJG+lVzWn1ivSWc1SykEtUzJoo11MLEkfxjAt7nJN94d6py2NeF0RfVR5B50PbpoSalVPTshrFr81YbhrXNcP9FftITWEWoxJTKBaru22pqFZ/G+FdRJaiO53c+cjYTT3NwawFpknebFmngsIadxzcu4yaFUAls1bGbODMwZqySqNAwo4NYrtLKMkNYe2XbWKydgcgO7ziqFU8i3Xs2i0gmMTYYlne9RQd+BcFW4p7Tfbpidd+KiJHD6Hn5J/p17XgORi8QsOp1fH+4eB+7yUI48k7Ii/OF4JPKbPrR5uHFAHTjCPmJa+KBReOIXHGuVTES32s+EruiBQAx1IsjqknSAd3Nf3t3OPJQQrF3n3pTNnIjxMcB27Uwdt0UUFhWubwmNZzypmVpHtZ3Yvt7Q3AZUYFo6tKqpthPdkyUAzVVvgOrCBIkSSGvGxIGLO2dmj0ZKfboGIaZm6wpB8t+amLpNnTNMqdF5uavazhy5Ml1vnsmuQ6ulPa6KXN2aFzQauCsna3S7N88hFSsO41zIJ3JIqWc2IAfaIiUeBnPWrqilKbMl99PrlReGGBn74NBgJKq8zMeZvRoLgQ4wmGjQcrfA4js9FLP1R0H8+eUwJPCCXesoBXf0uLgwa4tCDmlS0Vua6AIzvwFzeYFMBTeWpkOJNTjeE26tu55eXSPX5lkgzhabczlNKoPi9mZYXDdR3XSEpkK90F4FmuuRUDtcb135+rOK5ft7ExdKi3kKFp5sZ22QbCgGaw2RwE/qBqlCVq4/7N1j4GvyWDI2UKS9D/rHgIU3U4DJ+DOy3yOX1/bcNBEiXjMzgGVubRq5OG69NK3XNAojctAmhCN8o5SOGMNhyFEU6uXrUYSviXdSHu/g9WFC9Qg/S9KLyYikrM0RL+WO2HOYINzVyHmIYVS838QYyUOOBFAeNrC+W8chkfUodXrdDPKVk8yqe5sZgYfZvS8YmUlmIHWdy2FQ2yKJKZmIXBpm72a3s8W4N+anv4xq7+sGGQMujEgIrMaAdJOxUNcNKCzeubaQpfLm3VzfWn+brYkuaz/NWZl0jhpM4nGzmSuXg0JJxguahanFx6yVCJys54xMmUKN6mIsF1EFzsISaVvzYjRE/dpgwHNy7Xd7uwbDkK+r7rv5nP1hBq3M1mQRqgJEv9PRzt40exPpWEc1MMeYyAtB1eI8VOXAgJSvOWh/UeiSpPO6tzvfg4wfjXYUL4UhyWBrWZLOoTEp04iZ+er9IYZQAKupaXVeDBC1Nxp7T9vC7cYWZduaCHRMm5JaOWTMRmbuQ8zmpAI2DZofq75rXBSJxqoqeBeS9GIT7pGwgAKwTaJLZem+nJ/ae1SjUnxJ6TYeidSxDmYc5+/sPbpoLAzznrBaICEQVg1uNxj7FtBlw3je4rcj/XlDfTvgbqMLP4wmPJ34J5GPkhuTRxKXeOvRm7AXgYNiu4NWqumxYwbsPBnuegXx8Y9NAycWKhErSdXD+b33xCGFWtrBY6WBOTZQDzws/fvKI3mYIRR4yNyW84DFGuUDANyyzgJDGhePPSEZD7EdzM1K8P1g4cMUMp6QFOPziECr1LUZhAjA6u0tsp4bgNtnudiSwmcvYpY6qM2QLWLrCoickoncKjSedybSxfBKFgvjo6zXB+FEOD9Btp0BvPvOOCJdb4V4Pno3CztvjYppUleE8zXqveE4YyDUjlCbAU7zt7uoWW0GtHJWObyorAHXPhLZnIuaJX5efCRwciagJeGpdP3TuUsKMeCQKQqH3kbymAoOUakEnx6348TPkCOcJG1GZfiWnorXOauxFazZewFWKLRXykyPPphheVVr88AjG49y57qXIu9mbANXYAEu645IXUE1E4t0HyUQ+wFZtXC6zsVs4bUzJLbWDKvYQlMEtx0YHi+oL/fIxSkSqe50PTpOSEUGRqWu5/TyYmmNxZvGKn7dcu6LQ/S+oro9THPIFetPqCNzNmVxmhrpB2SYkKBMj89w247+9RP8ZjBAdLO3UOn5taWA43fQpiYsakJb4bc9oakIjae+3DOtaqrNwHjSsPigw390zfJ6h9xszANxUX3+DhGuAESDzgm30oiUXJ7sPRb38wVkw/u8E7dYRPJgfOgIE8ltKtJ5xJDoDjckvbdIKd8RMLqvKK9MGcfn3Gpl3/XFzfm+p/GyN7/6NOPlMCRxk7vThvMoS5NZhYkAlVplRiJaNiKQwxgdR9u9xxG5ODMy2Xaf9TymZU3/hRPUw7By1NvA1Dp8HzGDqcUtAm5rPBBJqVfijhxFkli2lsZNO3aRPqY1sSEpe+021ouG1FcnKO70xCj3TsxYRSOiTW1yAauG/RsrmsuK/lGDW1fUtyM+febj85xV0o1xTvxth+sn9m+uqG9HXDcxnjTghdB4/M4Ib9Nrp0g3wcUpYdGgtcNFASWJfXSSEcxCT0nT1XvzgMq2HwUucgCSFtXb8z0vMKESHJVCS+a+cVysmY5bYi7HgtmqB/PsoLK4ED06mHfpvXGE3W7m0HzGYTICr0KbH8jIxLPExBxD5iWk1KKFNcVOJGJAXyo0S6FEmhjTaGSqZcv4ucdo5ehea5FR2bztCR6mhRAqhyiE2tE+V7aveU7eG9HqEfWHmwN8g/NTS5fuO8MPhgH2asS0GHJI6oMTjQNg4c0UjEGajjdZX2FZLXNTryy4NE6wH9DTBatfv2I6aTn5hfcNA9l11kj8g6dRKDqGYBenhJMlosq0bizt/LjBd4Hu3FPtlGZSXA84we0mdFlDL7jNHq5u7Fy2O+OHeA/7LmfJ8B5pXZZlTNXGYF5jkqaEIpzJC/hI2S4BnaVBSKDqdCgXUUyS/P78d2FUMrnsDjNV7vd0Y+uLGRe5Xwk/f4+PIcl9r+MVRvKAI6tkwQFAd4yB5H4yqjkzY67vXNMiC8s85MZUGpBH50bzVsVvzM13fc3NF2t8B8MjQQJUW9g/gcVTGJdCtVOm1uE6B49XVM6ZbGLbZJZn/ty6tkZW+/2s43F7G7MZ7nB3THRu782opOzJ5HKtzsHiHUb8YMaw2nXGZO0XFuZdxcI67zNTVXYdTgS8w4eAes9y2xMWNfV1j9tZ3Yl0MU1+dYNbtDNTt6rMu4nZqXT9S2mDA3GhYcidDVMD9sx/ocC6jjglGeMoMzWllwGU7T7vFOeV3kPxnsxQTcc9BkzLeq4S3/oYI2LnMt0hvn2WYdW/r0KbhxlJKawovstA60HtRQONn8WJE0ehrlBsV0wVuzmzAohXI3ypwmsXxn243bHsR2Q6we8nTlrPtHC4LtB9VOP7wLBytFfRBfeC6ydjfYbaGnanDEZd56xKGpnOLs7CGDCP4fTEziPW0GTa+unajNPTS5J6mQkR1UZCSwY14Q/eWzMucWYkn18VEgaR5t8PhtdMk3lDTiwaj2Q93e9n99/5XPyXW4z6opnYOBohtKrsXDSGkdl4tZnoF7rO6pZCBCY9kbRm99afnTBd3x5Q6Q96E91Dn8+6NGWrzao+zMjcU6Mzq9jLPY9z572HxLfjNLDLnydNAwV15vsdRpF/ZUgebNypDi05EjHMSWDZwYgVqKmORk5P0Zub4sCStVp13+GqyuplVGGcaD/aId3A8No6SyQ2lx2h9qy3Vp+htUeGCZXo9osQzla4saiYjYsrT754znnhOdM70bpCdh3h9Qvc02s7/yePkO3eoKLlgvHNc2srsVog+95CpatraNzsqSSMYRzhZmMTvB+sA58zVmxScSOln5OhJWqZTMZWpTZ5R02G3PuIh8QUamSsaqr4TcS8dOxUn3NziyyXuO02YypSVTNZLS7GcLuZF3XyTo41S0oCWgpxqjp7J66p7wga5aZWAEctSQ5qsvKxizquQrW+pMLn8y5/JnnIBxmvPJKHG4n0xQywWmwcbBHFhkkSwxgib2GmZLvZyNzcZJr7rMuqIBYK6XZngGEKR7Ymqtz8+keW+egH9GyN3/fWma7ySGc0cynEjaQbZzwmVsym5t8pNNFhmDU7YnUwlzeWDn73qWVGQjDpxSZ6OZsdVSygY7vLKVBVteK+Oi5g52GItPahz8Q1vY1eSizsU4lNnVRhu82LQxYLC7tiM7AkUmQp7VgYOIyoWi0OIYUBARE91ByJwDDemypbrMFJHQnLbEfYbilV7u6fD2WoE8HTkoQm96uiWchzGMKk8OuFSvTxeKmo7xCsdYebwzGR7YHGjwKz9aUxhdr36DDm+Nc1tbmyMT7XWKwGHMTnIbXXDDoblbggmKYc6mjfz5M/VYJGcR+63hbObm91K7cWBlmKOCDdkBtya8ywSNfbe5Py+iJmearZNmceSqLCx/CCcbRz6Dp0t8M9v4YPnlqoFCb06sbOv1iE8vYbhXdT0NRTGNcPc3gV4rVSjQr1w5xyTrt2Ama7zgxfxJykjT2A+yjqVGqMTNOclu56dLfPhXx6c4teXWchJZLAk0i+JtPtZi5viEDzx+qIlAxfsHOM4czBa9OvqVq7/PtFI5QhjDUcu5OtyS1AYwFo+ZoHGilr82n+vczj5fBIipTbzCsYD3AS8UR8IL0n4gBRrDhTs8NcZ5LTtNNk5DOY03bLhfXWXS8NQ6g9XJxY+LLtLNRYL61orfL4p3O4JF2fO+qZdocVu7FeHTxuhXCxUG67OwAIc5XuMEKkvjOO5in0+9j8KpLTFgu4vrVdM3kkwXrzSmwpmhqfS/Qk8LE9adsYmDtFI9Z1UWMl9ulZLedzPVnZ9VsbLiM+Gk3Rua4yyRem3jypxaaLUg0p5NIZFE8GzFUyM4LHWfn/WCqi7OOch4aZZHa04O+t47nn7zxv4KCm5z6AtcwkZU4KkJuPP+B4Fdo81EjaDwXYmvgFmfEKVhFbMiYjAS0bkTRc0RtG9UBWEOctq+Md4XyFDBP9kxXqhNBYZkWmBW48w+9HnBfCosY/u0U2O/NS6gpdVnOYMAXDPm63ZALXaJ5V9gbsC9rPsr4nFRomT2bozRgkEeqg6LTJtUCpQ54ZidauzTDFhR45LtHopsyRqYLpXCEcm57jPSoButhhMLXZmCYToN5HbMMnfZAws20TLpXadpTl/2DnmsoCEsZ1D4A+ywRUM8U/hT0lyS29tjQwZQpY3IuzKUfVvzrqAfiqE3feeyCvaF9oliN4QK/klWbrQ45sP2aDcpAOtgdnSnwipKWNKbJKM7s0BFtoF+dWfBdDEHd2ii4a1Dm6N0/y+0Pj6M49119xLD9QpqVw8p2J7rxl/V5r4OvpCh80M1912SDbjvDoBLfp0CbVu0yGSez20BdexhBDrdpqfXTfmbFL37EsvU+LKmmQiIMwGp+jtzAsSSRQV7l7oO67GBKEWULhtp8JZCmLFUv9dRxzq4tc5wOWPZpCpuvnhux2sQ/xJzvg7J1kTGJO/5bGJC/+BLDGx5N2yEFD8lBQ78fCezimvmc1+urOIr/jbYgj6ZGUGZssn3hk6O4cq0w9P8BQYPwR8Ehemm+Q6PE5rk1kplyG7jLGkfQv0muywrn3BrQCLuqVgoUziYGqtSmX+W4iNI7tW7UZkR+zS7H5gjAu4PK3eqYGxpVjOKmtpcSyNW+kNSOiKxM40tozrWoYJ8JZ1PJIoOxiYbhF6Q4fNULP2E0WU45q9MuFtddsG9zjR9DUuEcX8+uaaDx8lHVcLZHTU2QR33d6Er2xkMHq2aPQ2XBEw6VdPxu2MUpEOm/eipPcMjSLQCVWr8YQpTfleZGZb3KAGaUmWym9Ly6n8XNlcKlDYxNg9jjuwTzK12rf3/9ee2F+7OgAh39GDCY3ZktCzyVr9p73fZYR1H2qfy/zeDk8Ei29kZkCn7CRY09EPHNxXhquENaB2T0fR5MthBxChSalAmFshakVVu8qflC6c0eo4exbE8HD8v3OuuZtO2ORLlvc3oBWudkR1sYerZ5ukNst/ibqmHS9ZWnq2mpstOBIuGIB+9gbZ7GIjcujtzUFA2Pr6G1styYqdLux96matxO9mnQ1zFhM9vldZ8YmZZsod9UkHRgQovfjFR2ZJQKem/eTOvDpbm9K85vNfNkvzgmbrXFz2jZzLHQccbEuxwxIXAjDfB/TOEj9q6m3z0bD5+eyN3DELTnAnrLnVLzuWF0tPn/wuCoweyZZU7YMcZIgUvr7IYa+Cm0ebsghYJZ0R2bNzmIXrSvzSESsWVRKB9cVDJJdb50mawp+EptNrVe5j4x6R2gcfh9oNsqgMC2g2tu5vPYLHaLQfLjBPb9Fl62lZrc72O1QcVHsRpEoiZhK95kmS+WmdOE45n41DDE1mZTXB2sFoXUVm3LFsKk3/VbAiF4xo6JX15amHSM/RSQCqfGfBjQyUlVjm4kinEnarzIW6u4iJHnKLIDkqtwMLKvDV5Vd49R9r23tOzvB1TW6WqB1hbu6NQOnzVz7lLJsJX7CocdZCkOncfBc9mKOms1ThDa5RifeH+HAkzhuwnUY8hxme+4zPAdNtR7IIzHz9cqQPOxIbMmCfVgyORMpSmpD7yVJKB41eSYKKIuPKvBBbbGL1axU3cj2tz6ifdrRNI71d0a2b7fUNxPr93rqZzvCskYmtbTws8uoIlbsiPmc0zlEuYIU/pyuoXJWWdxUpmcyTsg2VuoC7EbjkKRz7FOtiqWMU28cXGXGo47qbAUYiWpue0GYImNWEB+NgRMLb8CulzgzAikVHLM5UkdJhSgy7U5PzKuIbFnd741zstsXKeyIdzw6I5wsIASTOxhGA6Zre51lgArtkFRJexxWpPqb5I26oq/NxGwUivli1z7ceaysnYkHv2McXtSFL9f9hAkozvO4pecDjVceyYMNyZma1GaizHDkWDoBrokRyXSYHUgjCeyI2MKaJjQyZFPos/raR4RHJ6y+/hSc4+I7z83z2O7RZUv1wXMyUS7u9qm1Q+6hO02WSUoU+D4yI521mhhPT3GVY1rXVM93yIixVWOVsKV6Y6vOlM5NBXHTNIsg7fbI4wur6E2jH+xznZ8FozXYY/0AYcyGNHfNSyHVMOVWG9I29pplC9e3SF0jJ85S3/ue8OQc6QYLzxbNDKj2sen4OMKHI1K/jnt+ayDu1TVaeDl2P2UGYIvwo0zxH/R5Lu7nfa0rpKqzjurx48cSjOmzjr2SA6Gkkmx2DLaWxuk4c/QZh/LKkDzgmLM1h6LPR6m+VHxVgG4HGYpiITKZmG/CJwRisyhFVgu43eJvtxlv0LM18uzKZA9rb95FpLMnkFFTW0onFqa0raVj0+JoG+v38vgCVKnfeQp1hdvGdg/TZOzUGNpkNbb0FZMAU8JRpglUzGh1/TyhRWDRFpW3DTDBoNloWDgUG3UlOnvMXkltLTkZR5QKaZ0Z0FgjYypvbWb1yu3OMkVBjaq/XMxtLUJA2hZ5/5lR3+vqQHA6XS85CB1coRmSpAg41Cih8Ey8B5XsYcxp2ynzju6tlSnmzoE27H31PNw1HFlhLT6WPLSkCvcQ4s+KMIaHM0w/rPFyGJIEtqYdqqytAeayc58BQmBe1In1miZgIrOltGQCChetLZ7bDbJcxv4wMYV6bW0bdLPD7Y20JWMXuScyL3gnEByqI2xjb5cIIGrEKmSzQ/d7WxSLFun6WI5fW1o4YQ853Rpp6idrkkgzquYhxTYQOk5ReW3IZfpS13Pm5GSF3m5NMjLqmABmZBKZK1L2E1Et0df1+ta+73pthmi9soZcInB5ZfcmnrOVIdzG628AZGp9kUSkKHRXsvRlJM3pFAwf0dgT2Bcp8BjOyHJpAG6aHslgxNeUlbwp5DkgpR1naCKAmz8jhUpHWib3/lSNSmrhTn/ghxqvMJKHGmmOpFRb6aqWLSopsjWFVmvBu7T3eD9zHVLok/QxfFyAyXNJFbspXQv23pipQMNMhJvG2UVOrR5gFika+hhbR7nFot4mMWu1bY10luprUuUwGBjsnDXXiriKnq/tGpQ7/eXNXA0cIvAaAWWCAbysV+Z1nJ/A1W0EPuvZGFSVpZjHEU7mbnuyXpoxTJmW/d6+XwJdC75OxqZCyPdFUr1KIfcARPW3am7rkGjwSS4izG1ItDAiB6ni3L+mOszOMAPz6f4f6JuUc8PJ/Fz2VLjzujtENI68nYcKb/RHI7T5xKshIv+miHwgIr9QPPZYRP6qiPzd+PNR8dwfF5Gvi8jXROQPfNoTKXkViMtsyixopEdeSKmmVupkHBmReE72XErhge2cIaCvXdgiXRgJS1cL408sWvNYvJ9rWxIWEmnoubw+1fVEwRvd7Qy3UAN5VdUe23eZaKaTpWi1N29FvDNpQzCgsh9xmx3q3Ix7iMDzK3vtcgkXp8hqZTyZpUk8ptStDlGXpDOtWkmp6OXC2no+OrPU7uk6AoyjhTax6E63O0s1p/sh7pBlWlY5JzyKewx9vGeSwPJP4GKUXJFcLhGv10HIUzzvIph83FSrNA5SV3Nbz/hcDqW8n41V9mKO6AXp3NKG9mB6JGZIPs2/l3l8GrP6p4E/ePTYHwP+uqp+Ffjr8W9E5CeAPwT8ZHzPvy4in07Z9phEVFr/FAenrEmsKUnEpwy4AhqLzebjWpm8O1lHnCUSp2LPXrm8AZEcCsi+N89gHzkciVyWdnrIBW25oDDWleSRS9bjYxFMJX2vWL2cSXWxz24qFNTt3vCaccLte1QEPVla69Df/kX07ITpzQu0bZjeOEdbyxaNX3xCeOMRw1feJHzhdcbPPSa8dsb05gXj51/LzFtEDPBN55jObRwtJNtba4wDGYTJtEMy9hLDmvw76XAFFtI0c/0PzN5QqY5Wz07xCw1G2lSSx5EwtPh32BYezMGcKuZBUKtVsgNkzosVOE4HjFWp6ruZmVRjVKaYH2g8pCERkf+5iPyiiPyCiPx5EVn8IDb+4/GJV0NV/yPg2dHDPwX8mfj7nwH+ieLxv6Cqnap+E/g68Hs+8SwKQhoUu1LyUtLkO7qBSRkt97EJaWIXx4r4iHkIwRpnOx8Xi85ZGCyjkhfzOMJHz6yyN9G+U31KpL0DyGKBOzu1BVcWqQ3jHB61raVcz8/s9wSmRvyCKcTy+73hDbFGR/sePrrEXd6YdskUqJ5tmB6vkd60Wd1umIWpgfGRdfEbT1umdc142uI2nYG8U7BszDCasUzeW+Wtbqesbi28pnT9c6vQtItnD80bk7aqoK5xywVuuZhDutScPRpQKHCPxGyNxuJYVlHqag5fUt/hacItl3cV6/PECHnTyA9N00FGJgOmR7UzL1RIS4YlYycPQ0hThCm4T/Xvk4aIfB74F4Hfraq/A6u6+kP8IDb+o/H9YiRvquq7AKr6roi8ER//PPCfFK97Jz728SMS0tLkKidZytYchzIZRI1ZHBNirmajE0Iu4ZeqyulSc2MdGiYLYcCU4U9W2XuZtSdSgaAHN1nR2zRZeBCm2LR7gmoBvfEz3HKBXt2YwUjnnRp1D4NV2LKy86k89D6nHnWcyO0XpjETt3S3M8P1zMKaqrNCQN3vzQCIIF2PfwZVTF9LP8TevZaK9rdRZS3t0m0zYzzjNHsVMUxMMgQHTcOiQBNgeq7L1sSwV6souFSSwebUu1QV6j0uSSqUeMaiNTnJCAiHfsgGJWmjHjRLi0BpbimRXmd/5J93JAqOxYngMAR6QX1NPtYRJ8XC7k+c2Z9qPDDYWgFLERmAFfBd4I8Dvy8+/2eA/xD4lyk2fuCbIpI2/v/4+/nQhxz3XZG7wSYgIn8U+KMAC1Yzig+Hu1ICO5MrW0yaMoTJ8Xf8PZOtpgmWJqgsKXXsY2uIRQvvfwhvPLEFd7s1bGS3N8wBDPgsFmReiH3sH+N8FBOyrJLe3EbwdUKkQlHk7CRnh7ImSeVJUovAvJBFLBSqq7nQLygadjEF6WZgF4xl23Wwi3iCBkvRltW4SXdlMgxFw6y1AswCTEkgKV3PxD1JdTnnZ2Z0L04ttl9UhKYCATcEQuMJjaN553IOn/Z7mAJutSTc3M4ZNR+/8zAzUnUcTYcmGZtj1bTjfOtxNiaxaMeBg56/L8zI3LUEsxLblOfksUG6l8j2fQ59QLBVVb8jIv8K8C1gB/wVVf0rIvKwG/894/s1JO+LyNvxpN4GPihO5IvF676AWcQ7Q1V/BvgZgDN5rAc3GQr3VGa85D6gtSzkSo/76FGkNHJuDq7o0M0u77MrMxg3m1yPo5fXtrPuO8vwVJWJPSdKeuUtzeojOS7VpSSaOszd+FIZfufs8+NxtIkd8Hb7rL4uTQ06zMJFYFmW5B0lTyBmTqiW82c0xguRhTFhJV2XYTTvaYociPWKzIdJ17Sq5hAlEu+SKp00UQM2ZZ9CC01N//oaUWVceMa1p302EBofK6PJmadcz8OY07lzUWbI3qJA5qTkTNCBJzrNWZkjElvpkaRaHde2Vk1czKe5e+Oh51ECril7cyAncE+m5qGMSP6MT29InojIzxd//0xcRwBE7OOngB8DLoF/W0T+ex9zvE+98X/S+H4NyV8C/lngT8aff7F4/M+JyJ8CPgd8Ffibn3g0mQloB5yBMvVaFruVxiSlGVNNSV3lWps8aVNWIXoiQFZ610R73mzn3XixMGIZoJPJGep1bM+w7wygO+CWeEOb0vmVJfUpXZqKB6MXEi5OkH5hO3dlsbpkVmrESI57uoTJDFms3dFFA6uF9Q/2URM2fY8Qw7RoHFxKB6eU8yKmfKtZVcwyO5V5WHWTDWIaertB1iuad56jqxbvPfq2ablUW2viFZrKQOqunz2u1PMmFfHFa5NreQpv8nDzsIyZ8/XMeQmK1LFtayrwS+njaAySJMFMZgsHRMdSNDrT9osWo2VGJ89DwC3a2OLiBWHQ9zW+p4zMR6r6uz/m+f868E1V/RBARP5vwD/CA2z8nzQ+0ZCIyJ/H4qsnIvIO8CcwA/KzIvJHMDfqnwJQ1V8UkZ8FfgkYgX9O9bhj631D8wTJ7R0/ySM5NibpSEli4PgTEhNTxLRLIWYl2nnBVJUt1imGITBX5dbV7GGkcCKoKbhvtubZhDBnYKbJOvr1gxXlgT03zB5HaKPkwWDeji4b3PvPZuKYKrJY2WcmYFSE8clJzrZUlxaOpV6+cnZiHkjkx8iqNYMzjHCysu8xTga4xrogpmBZq9qEsfVkBV4M0L28IbUnJUzGu9nucUFht2d1szWsKUpJusgG1rZGz9YWyiXMZrOLnpctao1yCrmeJl2bUiEtaOxX5A6wM0v9Hsku6ryZzLooaXMZ5/Rw5JHc109nPpSad1c3oBKLQHcHn/NQ43vwSD5pfAv4h0RkhYU2vx/4eWDDQ27894xPNCSq+odf8NTvf8Hrfxr46e/pLJQDzsjsVYTDSZU8FMhpxQMiVOYwOGuFAJnncF+hlUkVWpPsfMzVyjI7EaDNKUstDFgJ2EWXXfvedDtSmwmwxRmZqjKMaCzWS2LS7nozq8JH8edUg5MBSRFbgKuFhUDxXP2mhzEgu85wi+GQgJUZwMMAbW1C1j4WEXpnxq1tmM5aXDeifsV40uCmwLjwVPsJtxvhfEn1ztP4vcUMVFWZ8awq68I3jKbM72N4smwJTWUSlk2Lf3YdJS1DDpPyEMmqbPduEnooy5i9i8jv0eGw1eud/sOR9JafKzede/gm5ftIntJ9uMhDVf8qTOHBMJKfE5F/B/hPsY38P8PggxMedOO/O14OZisctlgoGy8d19ukkepu4uNJ9+VOhXDcyTNfIeMo0WUveBJ4j263RkPHMhMA+uYTuLzJBLQctozj3CrCjbPRSZ+RuA/JyF0aZ8RaSPis8yFtY827tgM8Pp9bZCayW5hmQHffUf/6h3ac2+0BNwKn8/dI3f3CBB89s9BntTLDMgXkZIXcXoJ/ZEZt3yNDYDprTELhW89gnGbR7IQrpfArXf+qQkY7L1mvTAe3qXD9iFYOf7UxI5aA6qY2iYSqsgrkhO/UsYwh37cCTJdDzd68qNOUT/jHMAIF1yMZVTc35crjmA9SZnUOqrvnkBvI7OuHIqTBw2ZtVPVPYFFDOToecuO/Z7wchiSFAsl9zbvNaKXo+SbODcMzsSumfgluxkpi+jcrgiVAc7lEVpHHsN0Zpdw5mBykRVZVZhyGITNk5XZrWELbGJ7hvXk8y0XuhHcgVqSRR5J+F2fueSJ+lU2mItAp+xh6PbeMi0aavtS14TEJwI0eAUlMOgGNOmdhiPqu+HoGncE8rZTNud1CmHDvm5HR9RJ/vcHdNtTPr6wgMZH0VhYS0Q9WvdwPaOUP+CtSVdZeY7tnulghLhCWFUzx/BYVzjlriB4bg8NgtUeQtVnyfePQg0jiVmUri4ylTRN3eEjiOGg8ft+GdATGHmMiB8aoaUxc6b5m559hKA8a2vzQxsthSOI4cDUPVLKYi+LS8+JiWnScsRTIlaZJDEjSbtg05nqDGYvHF3ao9RKtPf3jJe23n5vh2O4NWI2hAGCL+2Yzp0NXCyvrT2Bt0hRJbnlS6ItaIRLFf2S7j9R7wz5SGwttKsKiwt32FhKFGLZUMVVN9Lb2+8g1sXTzrGvqsiteyhse4EsJn0harM7PavtXN+jZiYUfy2XWfwWyromul7ZoVwsrcly08zUYRzOu3lO9f4m2De4mGH4yTrjt/oD8l+6pwCwDUajc5bafMNc1eU8SvSqV86SyOqSyYfkhQ/Zu+jcVgc78kHtEi0q2deKtHIVHn328/PT3TzNeDkNSqMhrERZklXg44I8cx9LJDc4ur3PWUW+achZDmiaSwCpLbZ6fgSpaOWTbUVeO8ckpboi07cohux5diIGHm+2slerEKmZTAWA65nJpbvrKUq7ExZ1U0EzwaEVYNSaapIqMgfFsYZ87TEznC6pnMTvj/axdkto+VNWsyLZcGFYRFeB1b/1+k7SC5jS0VeSmMoCwu7R2F3Vt+MNqAe/dxAZi0Wspa2lSId12P8s3ioPb3j4rikXruMvZI7Z7896Ssn7Z0yd5CBEXkXYBOzmQyszeJxwIV5XZuPR46O+GGQdzKc0byByRQ53YGOKIWv9jOMjOlJKOGYStqgcjpN1T1vMbbrwchqSsz0i7SlnVe0Bhdjl8MRc2zDokJQmrPHziHbSteSPex1i+tipbIDSe8aSmvgyE8wXVs01kngb0/MR2/NGMFN6hF2d5YWpTw9maUDlLgbY1WpkR6R631LcjWgkyKuPS40ZrTF5fD4TW47qJ0HimZUV93RMWDW4bW2omL6rrZ43VNoZMqeTfO6ga0wvZd1bBO07Z85AmUvG9B+dxqyU6BcLG+DC896EZudTqNKgZp5ILExQdIx6z3ZmBjCxhE4puD7kpFea5dH2usBbnUImC2GmDmKa5OPDOtBAT3U6KahEfATJekr2UDMIyL3gNsabJE6JXVc6zA2MSQdWw29l7EimtxEcS5yQJHD3QeBXaPNjQwxRccikTAFukgVVntqWGWU1eHXMrzsQWTYalj/1zt1tzwVdL8zIWNe56x/DWKSrCuHKoWxAqyduEv45tNdsGWsw72O5tUTQ1YWU1LQDjusINdi6imNcBTK2nP6/wXUC90DdCtbPfp9ahXkBhahy+dlA7RNU8IrDwoPJWUPjWOfSj/R61ZCV22SMWMbLvzAuImRw5i/ICasJHqX5GRIwvMo5ZRiAXP0Yjo8NofYCSDkeqTyGFD1UEqXdZ25UqkuGG0Yxv8lKcIFGrxLCmGJpM/lAuM+mZuKKS18/qdln2MvFlClA16cUkbE2Hca7RiYbmoBwjZYCYcZSyyrnEUTLnJLFmH2BY1ubhMJcf1nhJDAl3YlhgrtNIaeAis5Nvbqy8Pe62lw+b2JWFII32A9LUuA8vYdHSfONDdNlSPzc+RP9kjQRFptg57oOnM9axWljYUVdoW6G1QxRcN9I/aoDAtLSd1+8DblLGlccNtkCToRmXjnHZ0txMqBdCJbg+MC0rUCz1Wle5N7C0TaYchmWN1DE7EoFdbRtEFfWCNgZ8+ps9w5MT/H5kaj3V5d50W/vBFnmsXmY1c08kXf+orZKLISOvJV3bpIRmKea4UAEmk2HEe6PHQxZkyoYkhUxgIWDlEA13CXgyN7E66FuURha2mufDAbfoiFZf0t3FydyU6x4uCiSDczQ/AbdeHQgvfdbxKrR5wHHQWS+OgywNzMbkiGtwMLmgEDV2BzF3HkNvaVPvbfENg+lVPL1GVwsWv/bUshK7LhYDRiNSxc9OJfRjQBRCJfTrBe2z3rySTvHdxLRwjCuHxA1uOPFUu8DUWAsM32n0SKDaBsYTz+KDjtB6wrJCWk94tMRve9NZPVmaUUmezqoiNB5/swcRQlMjw4RWDrft0bYm1A5/G6iGuCIKbotGwhzRuxLvDAvRMGeckkcYVdeUGF6ke1GMsogyk/2G3oBuJ5b1Sfc33cNU3Ne2CP3MzxaZWbFF4y2cs2sv1WF6uC7CjRIDKbgn+dCxE4FGjd34prup33ly2vvisSwkfBipRXgV2jzgKKo9ocjdR4S9wEUO3pWk+xwZdLX3z5NcEhV8ucj8B6o2Kowl91+MAt+2cGVMTmlt8huY6o2SvmzR9QK53RHWCws/ugEfAn4jTMsa1wfC0jOceWSCYeVAQR00t2ZEQi3IZOEP0RANJ55hJQyrJc11QCIpbDipqCuhe1RTbc2AVVt73A2B/klFs67wO9ODFcEIZ61laJpnOwuhmpqw8KhbU930uH7CXd7agh1GWC0I64UxVr2z4sNpmkMfcbBc4JrGUsFJJe3qOgLAyVAFskxm9Ej0+mYON4lGPaa9dbuzax2BYRHNinYGUtfze6ai5iaEbIzEu5nJWhqQOJekNfZypsmHApi/z3gkGn16bfRkXkhO+wxDkVeG5EFHIh2VVb4lLpJSdAUnIHkrIoI6Zte3VEgbxplEpQFplmZolvYz66KmfriA9p254CHYrr3r0EVDOFmitUNfP0OmwNRW+G2PTIoK4MB1E65x0MG0cCw/Grn9XE21V3aveZZPJ0ShP3H4wR4TtZCn3ipuUAurFMaFRz1s3m6odkqoxDCWM0+1UzQIUyOoE8Z1hd+bRyeKEXHbyryTzrIj9fM907ImLO1x2UQ5xRbG8wXqhLoz/Ejqek6NxmyLJNW4wSQDwnqJ64dYpRyraZPBiNcv++2RFJjDl90+Z9hm0aDCw/Q+NyfPYHzJWM6CU4XMYwpVctgUwdJSZ7UET/PmNeGa+qBGR6fJWohEhbZ7+SYPNH4EIpuXyJCkUCWCaqrTYeyb03cpVtfsrQC2cxRqW8k11mlCJp/d9awEH1mpstnlnjRME6yWMzs01sLoqXE+hscLUKhveqZFxbjy+L1jOKupn++RPlhDcQHfB8aVY2qdLXoHEmD7hjdseYKpFkINboT9aw7fKQTwPaiHequMreAmZWqFcSm0VxP1zcTUOsaVp72aCJWw+MDqQKZVjQyB0FZ2LkB1uUW9tSqtPrw2RikwPVqzf2OFGwJ+iO577CKICHJ+ahmV5KXsdoarOGvbKVc3pDYfxKrh3Ngr3YOS+5MIeQc4x0z4ArJAFF2HnKxz5XIKTVPDrqxfm46ftE5ixk+8gxQSxVCr9CSSdk0yJqHr7ngmiYSWzz8yYQ/6E3/WoaAPRJH/YY6XxJDojL5PGEcBDrI1Bzf56CamsOegsz3kSZtFjxI3Qq12JrWTlMb6uGSG7WJJaGvcvqf73BnVTW8p1mhExhNbrL4L9OcNbrT+vzhBBgs9JCjNjWVlfARYTbcDhhPDSCTAsBbcBFNDzBbBACyeKfsLodrDsBTqrVLtFRmVUAvN9cCgNeqhfd6z+dKK5fs9obHrFBpH88zqasKqYVrW+N1AeHSC7AbCylLUvg+4PtC888z4JLt9FDUKBsym1Gt0v+VkPXt8rid5ibJe2/OrlfW/Sa1Lvc6s3CkqzE3G1RH83B+obQ91ZKICHW0EbSV5RQsTz44p2LwBED3SqMSmKRN1vBEVbNjME6nnsOmgcK/M3mQv+OGo8flzfgRCm5ck7ySHsWr5ewL1CqoyKc3nkkL8IUBLWgDB4tnjOhtZr2els7ZBL06Rmw3aVoRIGMPB+GjFcOrpH7Xs3zCi0nDaoCKE1oNAc9UjqgxnDaF2dK+ZiNC08FS3A26wkCR7IJVQbxS/V9Sb5zGcmFcyrjBjEaA/E8aVEEycjVDZe6els1BmGfU+JugvatrnIzgY1qYREryw/cKK4SwSp0QYzheMJzXD6+ZhSVDqq47qxnoap3ICWbRZNjFrurpITLvdzOzW2NEwL9xYdSsJk0pA93oNzjgn0ja58DBr8wYzNtr1UYS6n9PM0zQbBO/NKKUNJoHBMGf2IFYq6zxP4pzKvJAUqkSeyQH28olT1c3/Hmhke/0J/17m8XJ4JOkqHaugpceKit+5BcIhs7UUxsnGJB0i1p7I+dpc9GULq4UptANae/T1C9yzG/NMAG0r64ypLdPC+gTffKlh/d6ICjTXA9u3Wnhcs3viCJVQbRUEqn2FOqi3FSqGf4wLwfeRmzIoboSpsfe4a/t9+b5GfMOMSbW18AaBfi1UnRmU5nJkXNukby4HhpOKUDtC7XCD4ncBNwS4UcaTmumsYVzaufTnnpNftzAo1BZGVjvjoch2H3v+RlwihDkTk3oTkxZtUZ0dxvneVJUBtanIb+jnRllBDyjy4n0Gc60YMnojGmIpgGZCnv09k9iydwl38ZKSFV3KI0YM7kD1rO8LLgkHKmuZSl9IGPwg8JEfBY/k5TAkwqEHUmqR4Gc8JKbtDjyQRI8PBcM1pX5T3A22a95u5gbfsQCNKFlAwHbZqNMhWxMHXr7rrP7EORbvmbs+ni0Yl572+Uh/UXH26yNTK3Rn3ngjrbC4DIxL8yL2jyWDpb5XxsmwERmhvbKwZ2osxEFhcalUu8CwtmtS7QNutLDH98q49qYs780bkVHx+4nN5wyUHVceUYdET0ZF2D82Lkt9awtruGjx3YQbzSgkbkdqDM5kcgBc3sSewzIXEU7BQhfvDCh0xnOhrq126KPnyKK2xV83WRaSujKi2mqZNV6SIpxbr2YWLZDFopJi3KR5s8iFmnC3PqcwLuaVFp5DJKtZ/6JDfgnEjE/qnZML/yZyapjZY3mw9K+Sw8bfyOPlMCRlGfWBETkcyYjkIq0XaZCk3Sm5zrVpcTCavKFiIZDcbm2ynZ/aBB3jTpc0RbzDX+2sU15tKUutPc13ntN/7oLQmMwgAt1FQ7UPFppMjnEpOUSRYEbDoYxLoblWfGe4iBvNO0ECvlPcBL4LDGtb+OahWHamvbSZG3KIFJDR8KHucY06ywaBGRw7trJ7o6G5joatFmQKdOceN3lOrjprAA6ERyv81R5dVIwnDdu3GhYfPaL9YIN7dmPGpqoQr1k4WzdbE8UWa5kxXqyorm/RzSaGAAUnBWJV95ClIFkuzJtQnXESMK5G1FBBBYm1OQoHbNd8XO/Nsykyd8bcjQBrBHotBVzNXkYZNqfyjFJSID9eH6Z+H1Dc6GUPWz7NeEkMCRlUnYu1krZnEfYA4DNyn5ttwwwAlsYkGhHxzlibTW3kjZTqTce43drOFmUARaOqWVXnkvm0m1I5hrcvGM5rI5MJdOdGKjPsQmIWBqoN+A7qWwttJFjIEjz4wbyFqRXGhTA10J8I7bUSvKNfCxC9FKDeBgNuuxDDKEuthtoR2nRsoeqMJDesLFWsXmgvLRzTk4rJweZLhpGESth8+QTfqeE5/QTO6o4uv9qCQnPlmNZtNiQ6DIZx7PYZrM4LbtdRjdE7iJuBdv1R/x43c05iB78McnpnNTvEFHCIRJuY8cnpYsgyjfm1hVdi1dbmvUp6bawSz80WSu3Xo6GFVzzT6VP9zvDAKWB5lbV5uKG5dWPoB1MSTySodENT+i+mfjW5vXEctOlMw8/GaCarhdnd7ieoZw6DUblj+jJJF6Z6kX2Pnq0zruJ3gf2Fp+rMy5gWBpxOLdQ3gEJ7aWnbJLrkRph83KT3iu/NsPh+Yn/us3EJkeqvHvqlcPLuxO6Jo9oqt2/XFr4sxD67Ne9EgmEv41IYW0ErWDwLbN/wSPD0J46T7w5cf6mmvTbMJdQQqor1BxPNZcD1E9pU+O3A8qPA8oOO2y8sqPZ+FnNuG8NAEkenrqBqzcuItT2aFmeSvawbKx5MWElq+TGFuXWoiHUj3BwW8Bl2UgORDpD6DSe+UBqRbyLTNHshaYQjMmOq0yl/prmWsjjJmCRMJeuieI6p9595vPJIHm4ksMst65lAlMDXybrIu9XKnivlAxKPIFUBp+M1dQYFs5SiWFGZiJgKexMFjxIzs6qsUjWpf11FSnjbGLay3ePGCdePhNfXnP/qltsvLlm/N7F93RFqYfWeGYjm2lK2Uwxh/KCgMK7NwOweC4+/1rF7vaE7swxQvTHma3cu+L0ZHpmU2897xgX0p/Z4f2ahRL0xg7J+b6I7t4U0xq/q98aq9Z0yrBz1Vrn6Sh2zPwbgTq3hLMPKPAYVsVAJqHaBmy8vqLeBUBsJLRUJ2iKb5tRtvokyX2enZAJhWYTpyU3XUytULZTzs+hTwlMS5pXlJQqi27GEZnospGb0cxh0EAaXgH5BNrMpp7maeP5eZjgylpK0Th4II3kFtj7YKFJ0Mb13LLUHGGmocCkPUr+ucJm96XZK7ucStUHCZC0WkucyTqaSdockVbi2cSc27RIfK35rXD8xrs076M8czY1a2ASMrZjH0dk/UWVsHeNCOP/GRHfqWFwqt59raK8Nr5gaYVh7y/bcKsOJsHwauP2cfffxxDyd/RsgE4xrRb4rDGu4/rLH9+TQyRX0B8siKbefdzRXyrSw54e1sPwoUHVKfT3i9gMEmE5bxtOa+mZg95qnvpnwuxE9WVlWZ7Bewhop6llNDsyjSwLJIgaWqhiwKmJSiolQmDghzkOYYt2LO6ynSvc5VvtKkipI5DOYK4JjyUMOZ0XmHjopm1NuPseGoJBxvFOFHr2SnDJ+6OzNj4BH8pLwSIodoKytgTks0QieRW6A1OZZpAbQB93oxZiXWUSn0MnIP4NNTpJehneR/RjZsuLmfi/VjNZr623RqXkLobFMTHMbGNYWlvjeKPPVdsLvQ8ZCJBiO0t4EQgXqhauvWLVvfwbTArZvK5vPCW5QpkbwqXp/D6GCYa05VOoewXCmhBq2b9lPdbB7Itx+QeguBJVIaLuxczr5zsSwtCxStVfq24nQOMKiZjpf0F80dOee/qxm+dGI30/0jxr2X75AT1fI6YnhSomBmhZvP5gRicVwc//j2BvI+0xGy2Po57amYO9JLTqnaTYWsb+QRjA1sZ8NjD3ifpQ4WSF0leU7sxSnHkhkZhJjke7Nxwim91q2zDhuLfrZhnzKfy/veGk8kpL6nhiKuR1nLN7LtTW+qJeIu41rW3t/3B2NuxCVxFprBk5VWTNugGEkXKytd+6us2rfxFtQNYwkMWNTxqbQOZmWVjRX304MJyZWVO0sCyMBxqVw+/kaN1gGxXfKsBb8deD2bR+NinkX+yfCtFT2byoyQrdU3ORMLzniLhLMmDRXNqGqW2FaKKExPsXyQ6E/tddqDeNS8Z0gW3vMjVBv7NxO3rUq3KkV00DpA9e/ZY0blXobqLeB6y9XjCuhvayZGlg+DYyrC07+zqUB0Gdro9P3Q/YUNGq9sN/PbNYSx4qyl3J2amniy2vLkPXxuajMNm8iIWIwlmWRoYc66uZC9oqMsh+B21JXJWmuxHt2qKZ3WAhqlHs3kxYTTgJ3jcaLqoS/3/FwCaAf2nhJDEkxedIjCWzNsXHklByPqPuZm0I7N0+ipNdaVVE5LArujBPhbIX0I2FRIZUznoIvJtL1rTWRGier/BWBZWNyiOvW3P3G4TtrljWcOGTSSD5zRj5TwyHqnT0eakd3bt+hvxDcAPWN4p9CPwiLD4WptUxPaKDaqhXS3Zr34XtjuQ5rM0BukFzDk97XnZqu3LQMjCvPuBYWH4HvzGNxk3lF1e2AjGbQQuNYfjiwe6PG7wJaCd1jodrA9k2hubL3Ld7vYgo4XudlazIMmyhEPQVUwqyWNo6GL0GWL5DXHpH6MucCu9QuNHYzzOGKi9jNdjd7DF2Hjn4uxKyiodnts6ea23DE+ZDqrhJxUSioBLHgM2UMYQb0c1uUY50ceLj07yseyQOOIkbMnghYOFMIG6XMS9YkLdpIHO48Ib8OQCJ4q5uthTEnK9z11ghpImy/dGrl/3XkfyjI9BrVLhhpqw/WRyYE1Hv8tmc8bZlqZxW3Y0CdZ/HcNEjqXWD3mtXmdI8EcGzfFKotczr3RmmuLcsSvAGyMhk+0l8Yy3VcGRjqemXxkUEwU2OAbndhqeH2uRHVqr1y/ePGVxEB2TgWz2yOjksIjRmG7swxfsXju5r2ysKbahfwQTj7xpbucYsEWH9HqToTd3Ij1Dcj1fUehpHx7Ue42x53vYkd8xzaFcVtKQI4OckAqD6+MIxl0Vh18dUNnJpnItu9MY23kZAWW1Qkg3MQdkAMSeciwLxxxHBLmjpyV2JPm2REmpWFXkc1OsA8z9J3gMO2KPHx3MzrATGSVzyShxoi+eaF3m6ULFp0382AaqqjUNspEpCmMT433ZFIv3Zu7lRXeBl4zP293cLjc5LAz+LDPf15g3rLrty+7Vl9ELj9XMXJd2Fcehb7gbAyboUbJvCCHwJaGTW9vpnoY/bFjYZHaGXGY2qgubFQZjixEEVT2lahuVaG0wiCnhqW0j0SQmMAqzpB1Jix08KMglbxuRj6aGUZnYSphPS4B4JhMMYbmWAS6luLu9XBsPKs3w2M65pQmyzB6Ts9oTaeTHPV4zc9YVmjpy1uN6CtN48k7tqS9F0TdX0KUHmmR6dICPRPVmbwbnv8ONn3eXSKDBPsuixDoFGbNhfwlXwRX80p53Q/wVK+SRYyKesXILw7O40EOZlD3uWSsI0CRcMItbdzSOzXYzkLiJsZD06T/1EAWz/RkIjIvwn8N4EPVPV3xMf+l8D/GPgwvux/oar/fnzujwN/BNuX/kVV/cuf6kziTZM63bRwkLkpf8+U6ASwenIv2SwfkISL0nsjOckmymQ74mqJbEwFfflsgy4bwxveq5iWFav3Am47mJziMFE/39B//hHBCf1ZTX/qGFuh2SibNx2L5wHfYwSzVlh9EOjP7Pf6RukeGfi5e8MK9mSCaaG43qQC9q+D6wxAbZ86hlMzMH4n+D4air0Zh2aruB76c6F9rnSPhcWHMJxaWDQtFBAWT9WEp2M2p/3I47qUWjaDtXrf5CGHyshxAOPSMlL1zchwWpvUwLvPmN68IKnD6bJBrqPqfF0hU5RdGAbzFKqoQOfMOxNVbn58zaP/9zNT02884bylEsFtdsg4T8fslZZarkMEW/dJYiASGGNBIFVswpXkCuI80s3GPJNpLgXQrivmUAHQ5jmSvN+oS1tm8lKvpYeikvwmCW3+NPCvAX/26PH/nar+K+UDIvITwB8CfhLrJfrXROS3fao2gIkUlKo5j6UCHAd4SUmFTqFOqrHJEywbnkgyKgVvlkujdzuX1c9kY6617yrcrbemTrc702htKqbHJ2glEMWHqn3SGVEjn0UOiDoDNkMtVmmv9jo3WIUvDtQrzaXVw7hRGM4C9bXDDfa8BGguJWIYtuizFEEfQyY13GT/xDyb/tzCn/YZTAsjnI0rialhm6/rdwwrqWKWtnk/RBElZyGOmOdSdXbthrOK+nq0qmf/GqH1+O3IuDbjgiqyHzIl3dpOePBqIUxnPYVxgAgrxdp3Xt7ghxEf+94QQtE4bAZDD9p5AqR2IRETMe5PZ4B4YsM2NUxWEiH4LFgFzIzXxCFKc8J70D5jI8pU4CSHNIQ7shafcchvBo9EVf8jEfnKpzzeTwF/QVU74Jsi8nXg9wD/8ad6d0FAKxtd5crRoqpzlgaY9SSsPiYi+tEjSVKLOgXLyrRt5oykiS9gWZo0obsBFo3V7jS11bqsGtx+pLrumJY1BFtwzW1gWDnGJVk53sV525/GupF0lcUWcHMtDGthOFfap8K4gsUHjuFM8XvBbwXEvIr20sIYC5XMGzFGKjE7ZMZERtAWxsbYtKFSfB+/Z8wQqYfQW/ijLvJNxHgvMmFAsNp5t9fBaoQmcAuPVjC1DjcEwmlDtR2YFhXh8QrXTaYr65zJMYxFBa8I9B6/24Nz+Ljgte+tPemus9BiGHJ2jj7GZ4nfk0JbF+8nkCqOw2ZrcwBmY3PENypFpDXNjYJun8HfhMOlbBDGK0lGhfwZRxXqn2WowG9yivw/LyL/DNbt/F9S1efA54H/pHjNO/GxTxiRXZhEjLJKWnnzmPkiaVfxc4bGte0dcEySUWpbU0fPJeZVTldaWnlCUuvNKPAs+x51jrBqcd1gpLLzFjcEplUFAqv3O/qzmuCF5TPji+wfCVOqnTlXQgOhVfRigNsKXU24qwqtFbdzdI8tTTuuwO9tQddd8kLMuxnXpg5fX8/qac5HUlrkpmgV5R7FLqf6aFyc8VNcTzZw6qB7LLRPFRcLAI3NCbvXzYiNvdUNTS20V2b43AgXXx8IHvqLJSis3t1ZL59ljf/md9CSRex8LEWw2qUUmmRg/L0PLRuWiF7rlRmOhd1LV9c5TMoEsQSIpjkh3mj7SYypbQ0/c3Zf6brZqzhiuQrM5MQh6vkWG1mu74E5a5Pmp384j+RHASP5fq/GvwH8FuB3Ae8C/9v4+H2m9d7LJCJ/VER+XkR+ftCYVUmKVEGzhkj2RmCOY0vdzmR0+t4YjZFsJk2dyUUzwSxqX0zRJU/CwqtF3EU1tlIw8WRd1BaGNDYRZQyM64ru3DMuHf1pTXfmIg4g7C+E/tzYpuPKPJTQKvq4N2T+bMBdVqhXtFLCwmQZu0eB/tFEqJXusRow6s0YjEtwnYVIEg1E8Pa7G4X2uYUu9bXgRgNRZYLVe4ktbLom08JSyOMqVhI/V+pdUQskds6LZ0p9Y4+nyuRhZXR8gP7Ucf2liusvejZvOfavL+gfLUzo6a3XizJ+w6Ly7/HeZiJb0uPtIimtHyyjppaGl0U7hyreW2fA5TIaiyYKH8XNp66tKRhkwSrtrdSBpjbF+FQFnoSlY6o5D1eELWmOHgD1RYHfNN1tnfFZhn7Kfy/x+L48ElV9P/0uIv8n4N+Lf74DfLF46ReA777gGD8D/AzAmXuswIFC1QFxKGEdqgdeyUGlb+QfhC7WzkQCm6zarH2hwbq8mZcS+7U4h9aVZYmWbTq5nBkYz9qcGtZU46OgXhhObKdGYFgRMx7QvRaorxyhgua5Y3/moHcwCjIJ1c7wCxUYH4+cvXHL4/WWygV2Q817T8+ZrmpEBemEJh5rXMQwxcPqPWXztnFRxgVoPeMgdcwUWagF178lII87lr+wxHfRyG2s2tgAXzOA00JZfBSNzNPJZAlOnYG01xOuU4ZTT3upfPB7lNW7DokEvOrpjiQslfOZqUYm1djEIrhc6QtIhRmS5cL4PrcbZLEgXN/gTtbmbSwaM0Qh9i1WtZarSbf16tpwkZQ5mkJsUzo3gi/ZrtlYlCzbcm4m6c6SuFaG3fBwPBJ46Y3Epxnfl0ciIm8Xf/53gF+Iv/8l4A+JSCsiPwZ8Ffibn+KImUw0yygenlre6ZLOxJERSXwCk2JMpeZuBuWSOpdzdzICstmZh9IPcHVjQkaqyDDh9hOuGy18cDAtPX0UHEoYiNW+WLpWAtTXjlAr9Y0wtYrcVLi9A4GwDEythRFhEZCd5/q9U9758BHf+ugR33n3EVPvEBXaD0xhflpEWcYIwvrOUrvtc8NLfG94h4thjwG4ZlS2bwe46FE1zyV4e/24EqodTEvJ2M3Jt6G5Utyk5luKGZr2asq9eaptoN4FLn5FWH6ghFZwUwAvyFg0MI+7dtaESeLQdZP5QbkOqrZ2FLrZZgKaO1nHfsvWYZBycdc1bHfofm9iVQkvS3T7rjucPym7V3q0qnPonOZYJLCVnkimH8jdpvYPMhIh7dP8e4nHp0n//nng9wFPROQd4E8Av09Efhd2GX4N+J8AqOovisjPAr8EjMA/96kyNugsOpNJZ7abKZonEHAIkKWdxvus5Zpcziy5mFJ1QFZC9z73o81Vv6pGUGsbxkcr3H5keGQtGvq2MW0PB1dfqag3yiCWKk3gp0zmGUxLW/QAu7cCejbgrmrDLZzit47QKn4rtB/GiufBEZqKaaFWrS9Q3wr7L/c079YMZ4Hq1hmHZC+MJ4H1O47ukSKThT1GgTcQ1g0G+LpeEHUM2xa/E0JlhDcJppUyrI3XMpxYCpoguMYA2am2mqHFpXUC7M7surtJqXbK6bdH6s1If1qA0Te7omBuQvyUFdXAUvTpemeZxqRlcnKS++UgYhvAm0+MY+KdGYzIT8nV4YnjMZFZsXlzSKzmssYGZq4RhQccCm8jAfaFBmzWBy4U6V/kzXw/4zdL1uYP3/Pw//ljXv/TwE9/b6chc1FWeX+SO1mOaECkaQoGrM43GQomY20TMmV1JrJrrQ7b2VIWoa6gHwinKyQ29R7WFcPa0V6O7B95Nm872udKf26YQX8Oi49sMXevGbg5tbbj969NrN/YoAoXbz/n3Q8uqJuRAWtpMZ0Extcn3GVFWAWW366sbibiIMMJMAn9mwPLbzZMKzMO3esjrovG1plBcAMMZ+aRTAsDXxdbBVHCTvA7YVwTwxorBNQWlh8FhlOh2lioNZzGaGAv+ErxA2xf92w+Z+9xgxmrJFA9rCurLZoUN0wWIg6jhRWpD7AIUvkZNG3bjGOxjr2FmtqEpJoK2fWE0wUyGPuYAieT9RK9vZ29hWkyAlqi20OWjki6rO70xNTWEnifaPNJikCOvN9UNRzpBBmrq2O6OIk2PeT4ETAkL0n1L3Pqr1CIl0J6T4oaDIi7W4jd5qsqIuk+TwqdwlzL4bxN4EWbf1o/4H5O+SZEf5hiPU1t1bt9VCGrhPV7gam1UCDUFrbsnwibLyjDSWD3uZHpdKJ/ewCvbG9avvjoknc/uGB9uudk1UEdTEypnXj8+jX/lX/w61RnPbsvDebdNIrfCdPZSP1RRf1hTfdawHUxg/KRp33m6E+hfSbUG9i/rkytMq5g+8UJ9bB9W9g9EZobyxzNRsbwk2kBN19y3P6Wke2XRvZfGOh/Ysvmt/Xs3zDq/rAU/r7/0S/wF//I/wb3Dz/n5u8dTDKyV9wQaJ73LL57Q/1si//wykLE2EvXKO6KjpMB2eKQ01PDsV6/gOWC4e1HhNM14XyF1t5YriHgrrbG6RExDyKVPHT9HC5BbtN5UBjoXW4pKt7beUBmybqL83netO2hUUnHyNyVuW9SBm9dMUcfaIh+un+f6lgiFyLy74jIr4jIL4vIPywij0Xkr4rI340/HxWv/+Mi8nUR+ZqI/IHv9zu8NIakvDkHGhMxs0JhJPLfiRuQwNckLdDUuJO17YgpA5ArORNPIPbyrSvzRmA+vreG3gg01yNaQXtteqzLZ4H1u4H2MrB6X1h8pDSXQrVxnP9yRft+hX9a4249rlK+9s230Um4fR5bQOzMa1qcd3RDTVDh9GQHXlGnaKPGSp2E8SwwngXqWOnrRsvgdI8D00qZWgtllh8IzaXQnyuEaDBiKnhcWUiz/UKgP1O2nwv0F8rUwP5JwG0dNIHqZKCuJ1wzWer6Qnn6D0z83rNvcCrK586uQeI5jIobY98bEauhaUz4OeuURP6H7vdIZVkXVE2jdbNHvcNvOvNgnFhfoHFC2yYXBbLdmXRiEoXOWSA384tyRXaxgTS1GZuYacn1WyFYAWCUKkibUQ5nSo+2aWbdk+NaroceD4uR/KvAf6Cqfw/wO4FfBv4Y8NdV9avAX49/HxNI/yDwr4uIv/eonzBejlqbNGLMm1Bzq4OIMW1JIIJDtzSxWwsuiu52Rj5LtRWN9fCl662XzaJBK8d0usi7kRtCbL1pVHa/nxhOa5qrke2bjckaNoYfIDE1Oxi/wt3A9vNKqJWwUKQXeK9FForWAbccef7dc2gD0jv6b68J5yPfbh6xagaeq3FBFu9VlklpLcOjXhnXFvL4nX3+ya85xpXhIqG29G5zZSlo/4Fn95ZSbYWpUYYzcvow/PgOgtA/Friq8XvHdD7iLysmr4ydh52n6ex7ffG3fMigFT/9/u+ncgG/nJBgpQGLd/dM6wZ1DvVRaT82Gc86H11noGlqxD6OEMMyWS9twW73uEWNu9oY9X21NGwjeQgh3ufYRO3AK83ZO2uFIRLV7FXN62yiYesHa3WRmLNSerA+q9zLemX4zW4XAdpqlqUoqPcHHKfPOh4wtSsiZ8A/CvwPAFS1B3oR+SkM5wT4M8B/CPzLfFYCaTFeLkMCmV1oNy4YeJeMSaqjmabc2Ch1edPtNooTkTVFGEczJpWf6ywWLWG1QFtvTb9HK1ZzfeD2cwvqm8lacXYhq7VbGjRmbVrT/ZBALK6DUJu04tQo1UaQMVLb6+gFDQ73bGEdHR6P+FvL4GgVuLpZMg4eqgCDZ1wYQc13wrhUhicjBEF6YXo0Un9Yc/PVwPI7Hr83vMINRpl3o2Er9Y0wrpVQKVpBWE1IExg3NdJOuDownY6MrUM6B0HQrae69TTPrao4eHjn/Ud86/XHfP3mCZULTJ0nVLB4NjKdtMgUrIQAYLuzlPp+yqLO1sYiegpJRW0dCWuqJroNBtJ2vamwTQFta9OImYJlYCBvMLk+pq6Qup3rb0JAnbPM23Ix0/VTmjZVCveDGZvEmk2tXHd7ZClzKBQUiJtXJNGl/tIPxmpN4+EO9+NY/dv/RUR+J/C3gf8Z8Kaqvgugqu+KyBvx9d8ngfTueKkMSamMlgCuO+h4bNqUmzdF11e8j+ldAYxVKZW3yZ2kFkXQZYvb7mETwbtFlVtynnzzlmlZ0zzvmFY16s2QQAQhPUgN7RV0F2TW6tQovQOtlf71YJmUnRAaxe0sszO+NoJT6ByhsRCmake8V6bJitoUcBP0jybqa2dxcaWsL7Ysm4GP3juDr2yo/+6a7rEyLa0ex6QGLHNUbQ2srVSotsaYbZ87useJql/TP7GyAF1OuGtHWAVOv15lLEUdhBXo3vN//8ZP8qXHz/mVv/UVms6KBJ/99pb1B9aDeHkdM16Pzwltjb79iGlRUT/bMry2orod8B9dE06XuOe3BsSOo3kCiZcTLMwQ1cxD0dRwC2L7kb5I60eKuifr0eQsXfJIC/HwXLOTWLUJU6mjVspgxw7XN7YBxVSxHOnf5DYoD8lqxTakTzmeiMjPF3//TORjpVEBfz/wL6jqz4nIv0oMY1700fc89n2ZtZfHkCRl+JjqdacnlvKLcWqmu8N8s2GWXyw70Set1ejmooq2LeF0AQHcOKFNjaji3r9kenJuxxpNST003sIcFaamwveBxeXEVAt9MM0Sn9j4KizfE/rzKOz8fmXA5/lEfWk7uD7qEadU9cRADb2junLwWOj3nrCvkN7hBvMkqm0ks10J1fmO2k9cXq9gdLivnTCuzANyvWWIhjUgZkgQ8Dus/sxHsaPHQLDK4WmpVDe2EMYJpscD1dOacQ3NpR1jWlqq2m08u7Di79y2SKW4jaM/N2Pz6O9Yt7+wNhKf23SERYXbj3SfW1JtakLtmFYVcr5mPG2pt0ZXVyfWvL3v57YgdWPeQMzUpX7NOoVc2V2GExpZyEluIItEG/Mtvsi8DQnDAbs2Ud9F3NwXOnkoiSldVJAffF4qJvzhUOQ/UtXf/THPvwO8o6o/F//+dzBD8r6IvB29kbeBD4rXfyoC6SeNl8eQwFwDUVXGRoWZTJRGwWJNSuOpRFyi7mqqANXNzh6LhV7STdY8e9WilYN+ZPzSE8NFvBDaygr7RJB+RMWDwPbNmv1jh9+bnqodzBZUqE2ISEZonhl3vbkS2suKqYF2B5tVBaMwnIxUH9aMZxMyOqbLNoZBgp6MqDjqa8/weKT5sGI8Uaq/fc7mVAmtsrhy7L8wUD232+YTBlnD/o2AeqivHO2zWUB6OInp5POY+fAW+nRf7GF0SB0YHw/IWCNjDGsq46CsviOMJw6oWL+juWFXcxOYFp76drKG5CeGN1XPNsgwcvq10f6+nTk/fjfk36UbDbuYJmMaJxX5xAkRwV2c2+5fM4PwKcVbVPMCJglQVbFrosbq3yln80wxbTS9WRerglXnHj2qSBXramKIlM71oAlXCq9Si9gHGN9LRuaThqq+JyLfFpHfrqpfA34/xun6JeCfBf5k/PkX41v+EvDnRORPYdX6n5JAene8VIYEYpp3tTIhmrKSFzKlPT+WwLaJ2a0dRljGidTUNhnGKWdrqu2ecLYi1B5qE3KeVg3BO9w00V/UVJuJ/vGCajcSasF3VhA3rE2oaFxFAtpJ1BWJBXHWr1foHkF7aUQv7QRGwe8FujoaIIffCfVtRfcooG2AziGDfTd/VTEuNaaZwXUgX+p4/Scvudm3XLoTmueNNR5fG8M1hS9uAK1ifcyt1cwMjbJ83wDa7o3JOCqdx58OTJ1n9WjHrlKmy5rmmTOvZSvc/vhIdWPK9rvXjTrf3Nh39ptAfblnvFggk8Z6JODyGgkBcQ53tTF5hn1Ul7vdWMuLmF6V9dqA783W+D5ppBAnyi9K1c742DBYe8/tzhb2EHVnUhFeDG1kUbRj9Q45M9FqdZK9Calrw9ZSPU8xr+xGzCS1rEYf3Yc76mmfZTwsa/VfAP6vItIA3wD+h1h29mdF5I8A3wL+KfgsBNK746UxJKWrqDc39rNsOeCcea1RIDin6hL6PowRaPUzkWgRGzeFgEzRmLQNstlT7WwihFWLG0wuEYH6erRK25uBaeGtB++Jtx4yq0gn39qCnZZmUIYTZfVebDhVQ39hQO14GvB7j+uF8cmA7D3N61vGm5b+XFl+vaXaCv6p3QYXCW3VYPiL35mEgBvAfW3JO88bMzqDUO0xtXpnOMjyA83sWzeY3uv+NWE4NQBYPYwrhQnCKtiiF8UvRrYfrmke7VlcbHn0EzueLG/5YHvKdz644OyLlzx79xy9rpkWQvVhwPdK+/6G8WyBGyb8TZQCCMFIfl0/p95vE7dEo3p/QU1vG9jurN5pubDUbGpgHo2B7vez55+alKe+wZBDXlks5pqYcTTjVBAadR+lBYYejd38mKYZrE3pYphTwiU+Fw6BDKkqOHSMvv/xgNitqv7nwH3hz+9/weu/DwLp3fFS8EiEaDQSP8C5jNBnVmG6qakBUnouTgagEC4Kc7/ZA81PZ3FuqvQFqEz9DMB1o2ludBPj0lNfdSbys5lorwN+p7TPTRgoNFiB21Nh8dRqUoKH7pEdd1oZejqeKPrWnnrdo17pNg3uumLxay3DuS3yameLH43hipgRETWjFSrDQpbf8Vz8FzXNc49MsH1LGE7tfbdfsGySevNSrn6bcUGm1tLH/ZnxU3QZQBRZTZmf8NaXnzJ0FTfvnPHG6oZ//MkvManw1c9/wO1/+Rr1s8q8kNirp74dmdYN1fWe6v0rkwLwzroSJkZrypSBhRrLBVLXSNvMKdztLkopThl/yFgX0TtIxMHUSaC7h6ZeJelNq/KmbUntWK1vjlr2J/VMWi3h7GT2YlMT9YST1PVsTOCQ/RoxvAdL//KwhLQf1ngpPJJ0jbLGSJK/gwOkPnfPSyS1NLyPqd+QMwG2e7UG0p2tbWeMNHhSf9rKm4IXQFPhbjsaMGDSC9OiornsGVcVvlOWz4P1m2m9FcU580CsObcZg+ZaCK1ayvZzHbpbsFybjMAQU6/hbGS/dEZCk4p6I7goNjScW0sKCbHe5jU1cx+IwkbCuFJ2b0K1EXafiyzP3ngn+0WwdG4bkM6hdWC6CDA4JAhMgts7msd7Pvfoim/8+hs8vVqjm4r6xvHe5oy/NP5OPrw84WbR0j+ekEk4+1XHcCLsH3tCLUDN+p1AOFuZUe6L6lpP3OGLexgi0zjppk7BZl+81+pqW+iJ4ZyOlZToxUGYcFGzJLeiSBolTdQuSS1FokCS6mDzpg+5P5Jud9Y4Pg5pmtzHOMt4lrjcMYcJHi4FrN9T1ualHS+FIUkjSwck4V+N6HqiSCdNi7LOZlJrJ5H6yjaLzEXJ3suus342QWHRMDxeEVpPfW0Cx+OJyQZO50tCE7NBk4UKSaFr+d6O7eeWyKRWdzLZ7tyfWjFdfWWgaxBoLq2pePeswQv0v3IGP7aFSVie9/S+Ynp/STgfmNYT05MJ/1HDdDFCFXBVIHSeN7/4lGVlhu7vvXiPv3P9Bo/bLb/44VtcXy8ZVNDBUsd+W6MCbm/p3MRVCZWgo2P95obu62eg8Nt/76/xweaEX3vvNVAYni+orzx+J3z4/3mbb39hoHm/YrtQFteWrh5WliVavzcgo+L3I9OiwqvidgO5khdMwjIp0dWV/Z08i6aeX9f1OfzMzcbLcCJV86oibY32WDiUvNHYnzkRyHJ4kjMvgohJLbqzUzMgiVfy7NI2sNqKPQVilfI9BqIgPSaVPvwD4hovubfxacbLY0hi+jdnaFI9RdIfqeZTtbae5m5KjMWtIbi5xbJaobtdRu1NV7RHm5rxYonfj4wntvCGU6up6S9qXGeu7dQ6qu1oJfytZ1h7ZGoMhL2orFl3K0y1KaE1l3buobJy/FAbMaz90MfyfyV8e8Xp+8Kwbhg+N0JjXoI7HQi9Z7oYkTpwcrFlu1nwD/z2X+OtxQ1vtVf8juU7fGd4xH/r4j/jb9z8BO8sLxgmz+bpCgI0Fx3D3sPJiG49rjNDVl87BhdorzybeokXC5F+8eufx91WVBvTik0VxFXs5rf4Th3Fo4kNviRq0FqNUbMZLU3eRryq63PTdaYAj87R61vLiOz2RWrVGcCZ7u8UDMuIpQ1WaxXm+9+2BrjWTZYiIGVsym6JqrNw0qJGtzsre0i1Wk4iWW6PNEsLpWDmk2S+STRMRQVxOteym4Hq4Xz8zOOVIXnAoQE4bOlok0tMRvGesu3kqjJhky7YhM2VwJFGnzrlZVo9UF/1TGsLkULjqG9H+osm4hPBPBMBvxtpValuBzZfWNE+NzB2/6jC93Nxm5sM8OzPTTG+2ghTFO0aF8Y6vf3KRPPMc/KrFeMSW8TPK8LZZExXL2xuFvgq8K3rR/y2kw+4mRZcTws8gb9y/ffx7v6MZ9slX3n8jK+PjuG7a9x3TnCnyuoLW26mNcGZZ9K/NnH6q57NFywrFJYBNzmW37TvObVKfemodlEgaQ3NFYQh4jKD1e9obdwWNygyWmg3LTxT41je7C2sAEup7jtbsIsFenObvQTdWsmCXt9YQ/a6xhjcfm5y1Q8xvPUmdLTbRw91RF18z76z1+67OS28WhnnqKqsLieB77FSV+q5UZrd8GjYsuZIoVSfwHximlms/atGHVjzmgtZggcYLzv+8WnGS2NIcoMiVwj1qiJRsBfn5kKwfoitJ+JOkoC2eJxcWBV3QVK3NxHcdmA8N9R+/9jYqyff2nLz5RXN9URohP60pt4Eo8l7+8zQeBYfdsikTMuK9Xsdw0kFOLZvOCRAH8Pu4cRYral5VZ0A0+94psbSx9NCWX7g2L2u1FeeUCsTEEZBO4ee7vgb732V6+2CX33jCR9s7eD95Ln57inh/Ipp9CzeN+/Df2nDP/jWt/nW6SN+9Zc/hzqlunb0ZyCDIK+PVO+0aARuU6e+FJ9PbSSyiZHY3GgNtUINvjM+SnMp+C7QPapj9iqBlPNPqWtIjdkhGwhb3HNoo2+cwYfPYTJNV5I0ZsIielOfFyCLUi3auYF5VRl2kjgjZyfm3WTWq889cnQIsYrYslzJe0lhdAqlLeU8zPo23sB5jWFzqjZ+8NDmR2C8JIZkJvwAc7+alHJLeiNVlQlLMNOixZMnXmbBRv6ApfwiQ7HydG8scaNy/aWGqlO6c2H7xgm+U65+vKa+tYzM8gPBjY7l05H944rwekW9S53nlOrpQHhUMy6E5iqmhp3pr7rBqnUlmNSA7yQrjEkgF+ENJxpbaNolqC494+sD1ZM9m32DqrC/XPCffve3cPLFa7abBeFZA7Xyy7/yBcsKfWnE7R26r7iot/zez32DP/X0gtWi5/l3z9Frb8JHv96CA78RxiXUt/F8JmhuNEtGhhrGiGXvX7f7Ut0K+88P+H3Nh79rYb1yPEyN0J+esXg6Ut32uGWTa1MQQVYL5PImZ1UyX6KpLdQE6+XrjR8ii4WVNkTjr/vOUvYpRGpNcjF1BiCyZKkapB+Q9co0UcbJCgDXq9m7KLoyHuiMpCLR2KsYb83ZcjFfXc3zKQG8PmZ2dg87/X8jj5fEkMyMwQP6cSrQS/U1YHhHFPSVxmT7DOWPaWDnZ6ZjFwWhdx3aWmhT31iM3dxWTI2gIvi9Mi2Ek3esNN9fmtdQ38LVV+q8g48tSIg9ZE5XVHtjk6Y2FL6zzEl9I+ZhRCPSvTbhOiOlJbUzo8Jbu06JGRnXA4NjGh1VPSGi4BRtAzffPeXi89dcXtW4nTMB6UZx64HmyUi3rfn17WO6UHO66miqkeeisSQgFhTuJHbOMwynuTbJgXFpoVeuEk4NtgJUO+Oi4JX9E+XkW7H3cLD7hgrVzuN3Jpg9tR7XT0hQ3JWFmLrbWaiTFvM4WogD9rtrY0haUNCn4p7v9lFC0WpitKlhae+Z1o1VH8diOq0c1eUO1ku09qgX/LNbq+cpUsZ5jsWKYk3ynapWwNe2ZmCK8zig6D8UIe1V1ubhRwa0IKcJ825QEIQ0lpknjYksQCPOdpa6mlOBk6mOp0ZN1QfX4Bxn1x1UNpmGi0UuzlMvdI8qmg/ipBFrdvX0d3h8J7EVg2VuQu0thHlqNTJugNW7wuZLAa0Ud2ParTJYRXD/xkhz3tFXC1gE6i/sCMERgiDfXDOeB07evOX2ozW7GxN/9jsxD6dVbn/5EXWwc3KDeRZy7enWRjD721/7CovzjqGvTPd1jCxgZ7R9iHT61w3DAXLjrXFlgkfV3gxadRtbhz4OvP6TH/KPvPFN/svLz/OdzRdZPFWaWwOYfaeIKsN5Y8rzvbJ8t49tJlLWKzapOj+xa3p5bZ6C8zMnRBUI6L7PXBIrtvMWyux2triX7VyA6QUJimhAhskMB5jSWjfgrrvMbAYsBEshc2TXHvBBiopz4ycN9nz8mSuAU03Xg038hzvUD2u8FIS0cmSth2Q0Ug6/lBKAXKFJmhQprx+BvPQ7Tiw9mKjTQdFVm2PcIVK8p4Un1A7fG1dkWFvbhsWzie7MeuBWtxaWTGvLdAyP7BxDbYuwvrbfqxuhunaMK6V/bULe7JiWajok1y3+fGB9sTMBsN7U7YeLCa0Dt++eQDQ8qRvfeBbMq1imxSu4LirAt4aFyCQwOPaXC6ZtBb3D7axdhRuM/TquLfPiRvMyEhiMxvqawej4SS6y2oK+2dGPnr/xna/y1bMP6R8F9q8LblBW7w3Ut4GpdXTnPjZdD6gIWjvzIEIgCRJJN8DVrYUPdZ2LLWW5jGS11lTUos6rVAbEysqep22MmbzvLVO0H3BXW9z1DtnskW6gutzhLm+Rm62FPV1vnRLHKZPUsqeTqnyPRYsKinxZFHhgdB6IRyK8IqQ92Eixau6il4DTAieR5XIW/YUZFEuyeDBncCqy/oTFvAuLt5cLdNEgw8R4vkBrh99aZaiEChxc/VhLvVXqreEY/aln+WwCsUre9rmw/o5jWMP5Lxnjs9oooZEopmx3fHg8IcvRkgOj8OS3PuW3XnzEL3z4Fl+6uORyv6RygdWj5+zGmt15zfvffmRpYcEakD+r6M+V+rljeBQyPiSxQ161s5DE740CjzozQFVquhUYTk1TFsg6Kfs3JvRkYnhNOPm7NbdfmbD+PYHmw4r+ycji3Zr9j49wU/NsPOXkYkc3RRffWRp4PPF0pz5zam4+X3P7ZeHJf1GxendndTbDOBv2KMNIcOCLToppk0ieZ9sgTQRMU73M6dq++868DOljm9VpTglLP6DbPawWZjSubs2IXF7FFHLIAkh57sVmaxk/id6IxqpxYCakRcJkFq9+qPGSG4lPM14KQ5IBseLmqJpORUbTd7sDIG9WuIqTImhO72miyEOuEpbKz5KK0dWtLvdMqyZmICa6RzWL59bac1w4XARWZVSWTye2r3vW7waGtem39ieWNt2/bmFGqBU3CNXfd8W0a/jv/46f4+eefYWLZsdXVk95Z3/Bf/srv8Dfevpl/tE3v84vXb9NHzybvmFRjchigpsav3VMq8B0EvBbI5al/jbVzryH/ZNgKvHBvJKpjbPRaxRgUuobSwG3H1RZaMn3Ysf34DfWJjTLbtSB4SKAGpbSfmDGc1Thdlzz/3z/J5ComL95y7F63zAWOjuP7VsGMle7QH/R0g4TbtvPmZcwd8LLGZNEb68qw7ZUwbc5DauL2sKWpiIsKsO0dgPSjzPDdZBsUKSpoz6CFW1qpMUzzOLg2csF6CMKnkLnacoU+awNK3LY2TGp1z/I3H/5vY1PM14KQ6Iwu4pFg/AMdtX3nGZiveLnLns5do1eTMoW7HQuMx9GaGrzRFTxm47x0RL1QrULyBhiFzs1F72C/szjB8MFfBdAHdXemoi7UfGD5+bLYsBlEwidne8Xmmecvr7nt7bv8bXubf7Ji5/n37v+XVz3Le/sL3h9ccv/69d/nO62xT+t0bORxUfesjHngfpZhevNw1i96xhOle6RGYkf/53f4RvffULobXGcvbZht6+p64nt9QKc4p4vkMFCovHUqO69V6obT/3cxd6+Foa5XgitY3y7QwA3Vvi96b02157h1FHdmti0680DG5dJXd6Mx+m3oN4qoRKrx1nWhLaCixV+0+Fu98Y2TdmOtFCXC9vtl8u8UYTzNVp7xnVNaBz7R6ZYPy6EqjMlu2o/IYMZfrcfcZe3ZoD6iG10UURJIn42jFlr5KAIL+FxqaZmmjIXKXvJySjBTOF/qPEKbH24UdLjNTVcKiUX4zjQgyhlF0vvJCqilcK9Oo7I7RRTyAF6h64XhvJfdQwXC6pxZFhVaGX4SKqmVQfdwsSG3NJT7W1C+15xvdLcBNpnnubKuu9N2xX8tg3/2t/5fez7mn/yq/85z4Y1/+vLrxBUeHp5ws1uQd97xn2N3PosMJTqbKqrFEaYt3H79wzIxlNtHOFi5HK35OJiw49dPOOj3QnX+5blYmC3rxGvLFY9+98KIoprJqpfXyNf2pi3cbOKWAnUN8RwiBwy6bbKxX8yQPPc6n7qWwv3ls8mNm9YY/FQwfKjial1+EhkS+zg6byiuRpx3URYGbXdXW3Qk5lnIuNcdKmrBbLv0CyVGcM8Z0ar2tm1mWpB1t4yZQH8tjOvZbVAbrazt6Gp9WsinhVeUFJci0BvLstI4XDbHmZ4itRxIk8+2Nx/5ZE8/CjbdGayUCE0IyVDNQTkZG2sybIlQarnSCFTwlHa1gDYEKA1weFwvkZFcP3EuKrYvFXRn5mR2HzBsjBXf+/E6tuOYa3UG5vUzTW0z2M9jirLpwF1JhrU3Ab276xB1lRnwl/8m/+1zGRdvu+oWuhXFoLUoxXfqfOZXi8ThEYJbWDxbsVwpribChmE8OUdtlaVz5/e8NbihoUf+dCf8M0PH+O9MgK7qwVMggZher9lOplY/uIJ00JprxJmglHpNyaE1D0Wql9e0F4q3YVdyypVI0fOTbVXLn+8sszRCM11IDRJvkDpTxyrD0d2r1c0t4Htmw2LpwPDWcXiA4dvDPOQMSC7Hq2A3d7A1qubHLLK9Rb1jmbbMbx+wtmvT7g+ZBEqjVkbf9sTGk8VX58lC5I+a+ISpd43Lm4qy8VMWEvhzFB0Yiy1g2H+Gx42YwOvMJKHHGXsmX5PXfZS0V7GSIoS76ykVrwmAasUmIuklF4IxkMQsc56agSjcVXhJstq4Iyd6gZjq0rsKaM+MlRjinT7ptDcWoijXnGd0gqWMt4GxoVDJmvG1VxDeGbiQvWNhRQSoLk2g5Iahvu98TjCIlCfd/g3N+hHK+u+d9Iz9p5mOfCTr73H3/zOl/iHvvpNNlPDl06e83e/8wbBK3U70vcN7tZS1lOrUCur95Ttm2LckmFudL57bGpq9Q1UozFbmysrTNy/JjTD7K3cfgGmZWD5nmP/mvXg2b5VsXo/UG8DzSYwrpy1B63E2MK10etToV91meplPHK9P0j3EslhUpsUJm1D/X5A2wa33RNOlyQFO6aAJJ2R2NSqFHzOXJFCdkDj3NFt7FWc5lUE7xPlQDlKDeeJFL2ThwpvInfnN/r4REMiIl8E/izwFhbN/Yyq/qsi8hj4t4CvYG07/2lVfR7f88eBP4LlGf5FVf3LH/sZkHP0wB06cq4ItoPPIU3himYKPfG+pNgYDnQmtI5aGTEdOK0WuL1JA441PPnPb7n87Scg1mTb94qbHH5PDnWqrXWa05217QRHcz1Z9e2oNM87+kct40VsUFURuR/K9k1HtVHAUrKps97Uxn+LCIqeDjgf2F4vcKuRcF3HPKB9pV95/gany45/6xt/P1Piojxt0Ccd0zsrpFYzRreeauuYNoZ3yAQu6rJWW9h8MVDfzIvCd8r+sRk5fw3L9zWq0wvNtdJcwf41x/k3J66/7BnXkUinMDWOahdwfaDaWqo9NKbvorUz4LUfTQOmi5T1RWP6rVGIinFEQ7zHXWe6IbsO2VrDLHetxiWJRkRbe79JKE4zQdEZVpIK/zTNgbLCN1f1TlnbJpVn3PFIyrafD+yR/GYJbUbgX1LV/1REToG/LSJ/Feud8ddV9U/+/9r711jbti2/D/r13sdrPtZzP86559y3y7Gr5BgHkoBwFEUJCk6IqBAlKJEIQUQUH0CKhREuOx9AyhdLgIEPAamQEUY4mIgqCRMliCSyFRmiSqoq5djl63rd53ns59rrMR/j1XvjQ+u9jzHX3ufec+5Zt84+9+wuLa215pqvNeYYrbf2b//2/xtjfhEVmf2zt0x33gH+A2PM3/ejJNwO0PAoNpPxkjl1fta6y76/M8p8lh9IVGrr1KA6SvEZkYmgNEDxzCOrhvrZnuGkRpzl5Nt73E1LaErEWcZVoR44tcou2iFoALDgusBwVGjabaE/LeiPliqI1KZsw+JLoElli9pIaDtG6347gGlg//Wexbcr+l2BlwKcUL7XKGD6eMnXfu5D3n9+ohoi31thBMYvdYp9rDyMluZKiXB+oUFMrBLQ2vuabYRC2eriYPmhZf8wegVH5bVyK9RXCpqKATcIbKC59GoqfmNong8UO80+ijbg9gHXafZRPd0q7nC9RRqdczLJoCy52o0zYaJIp9dOjcvAuDEGbrY6FNjGNv04ahbSqEueGX1WS5MiOgdE98Qgghl6Doy0YkAxUWAa53SeK6nxMW1q87Imyw8Ye3daJPnkv9un+yzWjwwk0Q8jeWLcGGO+hXpf/Dx3ZLqTuzaJY5ACx5zpGu0mUtaiQ3tmmsOZCc/kbk0QGJIvSvSfHUeM1LoTxfkdu9lpJ+exEE7W2O2esFpo92arLUmxhtJom9GMAd8UjKtCfViswTfaBbGjMCws3aniKNsvWS1fnGYA40LLnO7UUOwlM2qROCx3VbD/2oDpLYwGaTzjSjOge3/ogg9fHNM0AzcfHikesgyYqxJpAm41woe1Yhq92lFYD+W10N7TwNefqDNgqMlykbbXoDIca3vY19DeNzRPRZmuowZMIPJmAvsHJdVNUGbrKPQnBcv3e8pNpxlGVGbPwlFp5H8OaqZ2/dFShacKixwtCU2FbXulwY8eE3VLcFYFvasS2ew0sCSLiCoS1vadTgCbSRuWpEvSdVkUS8ZxKpXhEGyN2W1qHSttHz0/I7HxLjs3XziKvDHm68A/APwqd226Yw0QJzeTNUUqcYYxfvCzyeD5YN9cvSrVvTFVzXTotDOVcdTcukyUkiQODdgrnWYz/YC92iDrBe75Bn9vrfaSpSPUBa4dcZ2myc2zlnFZsn9Y0p6a2D7WUoiIfyhTVrOR/kQzkeFId/xxrUFm+UgYVhbZRoOtAspHOuvT3heef+cMd69j9/0VThTDKTcFzRMYlw7rS9r7kjMcMbB/GNj87MDx2Y6q8PSj4/rxmub9ku5ceS92MFl/dPM1PasXjyzdmWHxe57j371hXFeE2tGdFtSXA+OqUlLaQlm1rhOlq+/7SWw7AZwian610elcWS3g4goTOzxm15JNzQDb9vRvH1E92erns6gx3aCt42jnmRz1JLZ4JXhoWyQppM3xDWOnbDUFs5TZRleCXPLMMpfcuckBR8/P5Kd0J+uLgpGkZYxZA78M/GkRuT4oMW7d9RW3vXSojDG/APwCQMMyZySJsAREibzhUDntVv//1pNOHR1rNT32fnJ+I9a+S2XJSqcyforszwDdcVSuQ+F0SjUi/2aMMx3OgVNpAd84fGNB6SU6/IaWB+NCVeddp5kI6O5jxggqV4pXVFd60SdsotibTFl3vdpxIobNEfiLGueV+g4KghatBo+ija1cx/RlwRaB3a6mPtlGoR5D99BjRsNQe8orlztOZoxDfo2+r/bcEn72mOWTAd9YQmm4+KMNxR7GWqivAt2pxfWR/7FqsCFoJpKyDqefBUVBuHeMvd5N0g4SS7xEIut6ZLWg/s6zl8uRKNjNMGhmkdinUXMka5GkOZ0EtqayJD7OOKuzP/NN6LaUYgoeM41gAKpZtnsHK2Lzn/v1sQKJMaZEg8hfEZFfiTd/KtOd6BD2SwDH9l7ML+3BoJ4k57zbZU36UJOafPqeRHyTcpr3k7pWUWCO1sii1l1rvdSL2VrCssJ9eHG4K4WA3Owx6xVms8PFnUuaEndxrZyFwWN3FrsqMV7ozhY5ALTn2q1pLgL7B9rFqK6FMeIkrlePmDT9qeQuYdFpJlNsVSTJ+GgEXoKsRu79f6vIbdFyS5wGL1+hbetWtCx5LvTHSufvdg2hgsvvazQrjL6uX4i2l4+EsAiMa2H1/YLg1E6jutFp6M2XLZgS48H4CBQ32jYWB/VloH6hQKp98kLxi4SFgM7FpHPp0XMkCKauYvvVTW36aGPBzVY/ozzZHVfwU8DoerL9iHMzvVZUjiDJOkbQNJHeEEG2W802ikL1TW5uYsfIT9lGmgSOWFua9jWRt3S3dhR391Sf1fqRhZ7Rq/YvAd8Skb84+9NfQ8124GXTnX/RGFMbY77BJzXdCfrhzUFWJZ+ZQ6A1sRPThx6C8g+yTUVkM44jpolesPsWc3mjqbQI4+lSadWD15mMVMenHXC1RAqngjrbfZzxGJGyUH/aXYftR8oXe0LjWDwZePuvP6W5CjQXkSHrNaMoN8Kw0gwhYSmLZyFnMNWN4BcmBw1faxbTnRHp7WB2Bdt3lfQ2LhXnkEKp6d2Zlk79SZzaPdWsYliD7bUMKrYGBFW9l+jB85bHH4+Y0UTyF0gJrtX361rly4y1wQ1C83xk/b6nuRTq60DzpKO6GZAitmQh4yEHczRdrwBsFCfCWQVYZ/Mspqr04rbqlGjqegLfYZpxSSrxRVS+S0N4ValYWJrp6Qf9Hg3ScFbLKheHAeNQpzk/02BTlQetYohEyb4nGa9l57+Pzsg/8fqiDO39SeBfBv62MeY3421/HnXtujPTHROHs3KLLcnoGaO8kMhmtXVN2O2mdDCVMbctFvftlK5uNhOwN45wrOPsxbMbxvtH2H5keOcMt+0ZVxVuN6jBeD9ir3aE9RJ7syUPjRmDuYzKQFE0qf5ui5QF4XjB6r09pvMcl5ZQK1O2Oy0p9nqxLp4O3HylYn/fZm+c9r6opWYjFBslXBU7DSD9sZZHq+85iq0GGCMQUhdSNKCYDvoTBXOHFVnpvnvLY/caRIb7I8OpwW0c/nikPmvpP1xRXhvqi4LyRsWWir2aa/ka5cXshRBLpfpFz+LDAYJ2wcaqxl33jGdLxnePca2nfKGBl25QvZBhRFYLPX5W1At40xLun6kXc8QiZLefVO6K2JFLBLNEVhwG7caNI5yfZnastG20bdUgY8pSg0CSZxTB0CP3z7XzE+dzpN9qcNq3YGN53XW5K5hsLPKKuq13tl7zIPFx1sfp2vxNPrqM+yc+4jGf0HRHsuhMSi0TfySdPIAGmv2e5EgPsZ0LE0PRe93xEpsRdKcCTZujWPHwpTPKD1/QvlXj2oru1FG/KPELS/3MMK5UHNm2Ff5Eqds4i9QVUhZ68qc2po3gXQi4x5eEe8eYfiS4Srs8HhaPPf1JxfZLBeNXqqhropO4ZjRaWgyG+rmjexAov2fBxt0oTJlMd6pZTkjZQwHN08hFieVGd6b6I8Ve77/6ng7f2RFCVVC90Pkb3hkZ+iLP3NQvREuuqHlS7gLBWfpTw/r9EeOF6kXLeFTjm4byumesHb52hLMGX6v2rfHCeNLEGZgS0w2EusIMowbu3YCUDn++xq9K5MGS+ttPCasFFm2XmxkzGZiwj8xqnlHs5+BpAlqDIGP0zfEBup0CtADPXyjOlcsno8Eiyngqg7oi7PUzzqVMlPbUc3QCqD/Vki9g1+Ynt+JIeULXzWzM3IdJXT6EQ+p8fnjSHImZSRQQzlquyU8lCDw8hxCw/cj48ASA6rKjPV/SnRXahl0VhMoyuBK/LFRr4+wBdtRdutgHio2my8VVO42x79R20j6/BmdxNztc4ZBlw3jSUO5GitZR7IX21LL+gTAcWWwHfQ2yHukEpFZWbHWtBLFQTkBsudNyZTiKpDgDPg7P7d8SxlWgfuomVq5X/MT1CswuP1QgNwQDj+qM9tkBNl8xLJ5CfSn4UjtKGFg8CRR7j9uP2Os91S52wKJaXXGtH4PUGrB97dRb2UAoHa7QzAxRvdsCkNJhu5HiqlNg1hjs5Y2SBZN4s3NQGs0oqnKSQ4yfuUki0BCDehR4ToZpqVUbOzmS6PPGavaRvIIjDifpuaIhWyJF5swk6GMPCJJ3sb4IGckf1JqkFmE+AZzZhWZW6iRd1jnoClN9C/lESLeZxUJr26sNslrgnl0zfumM5nGHb4oo9iwU2xHbjbkkGRtHczGyeUc1XssbrySsbRxd1zel76GIY+ptp25uiRz17AXOnBOWJev3oTstWT1SNXr+riOUsB0Lyk3BsAbfWKxHW8Rr1YTVuRahO1NR5vqFShckvddQQHllKDY6TGfC1LnxjYoYFVulwe/fDaq81mupNa5Clh5IbWex2g1SjGfU1u6gHTB11ovBoO0ZTleUL/aISDZiD84QSktYG8xRMQ3etUIoaortiOlHtUqFqeyBjGElRTNAdUYglivDFBiSTWcqPdLMVVVO1hVlYkXbfNuBV3CiF6RMdxgPJs4PyJJzesIdrdcd//g46zUJJDJRkOfYSOz95+ndFESSHSdM499z39bU4vMh12SyVxo2QWC/Vwziu33GWOpSp4JloR0Hc7TUB57U2CFw8ns7huNSKfBXoxpsPW01GJTKc+jePaXYDrQPGlwX8JVlXFmqyzF2Ykbcpmf9KO68Iiwe14iB09/V3XxcOvpjx7DQzs6wirT8PREoVUNzO2h7ltiB6Y91OjfEIDKuDOUmljOF+tZIocN5p7+t1H4jQID+1DEuheFIuP6m4eg7htUTr63lnWc4ctTPFXMSo9wYM3iMCKGpcPuBsJisPfT/dthB6E4s1UaPxbAwYART6bEYj1WdrrjSealMYrNGhzFsIhwOU7BIGEh03lPsQ8sgUzaKnzj0/EiyEZHdnANBZE5LP0ziSnEzMMkpcL4SXT7ppoyjEifv8PT/vK/XJJBwONuQlLqTEM0slcwly+x3YCIOJYWrEHVMjNUTLynO19XEMUiK4v2gt68W9A9XFJcd3VtLXOvpzktsL+zvOeqbwLA01JeOYu8ZTxtcpzM2sih1KG2pk7HdWaFSA51kFijG4NeVOtONHqkL7E7biLb3+GUJS82AOC8UKD0z0S9HOSJJiLrc6a5YdBLnfYi7vbaSi61ottEYVh8G2jNDeaX3ac806JQ7wXqh2Fsu/riezeWNwYiWNuOxxXWWchvwi5LiJrJ8Y3gOVaGfkxeGsxoThGFdRGEogynVRGxYWIpOqK8l+hgL40pPvWIXJQa8TC1EH3Sb9mEin6WhvgSaFgX4uFnUtQabJBCdytmEt0VsxTYLNSX3s6G+MLV69UYlPYauw56eIJdXcRhwnES2YOKm3MF6k5Hc8TpQ6Z5rk4COjkf0PrvyJYJaxEay/GIitUVwjDEon6Trcp0sEnJ5ZBqdsQmVXgTdW0vGpVXs5FSlF10vdMeW5eORYu/R0Xqv6X7pMGNg//Yi7y52FFYfjJRXLVJqfW13g3YvmlIH1/b9tHOKpXixw113+KOaJgjtvZLlY73gt6IZxrCGUCsrttgpmNr8/Zesm44PvncPu3Usnmj5U7TC6rFnWCp9f1wCGNYfevoj5bkQOS3VpWVYC8ffDZG6D4tnI8PaaYl1VLD7Uk259YhVzVY7BFyrv4NyTMqbkXHpqK4DdhTqK6JXMMqCDUqpH44crhXEBcZ1hWtHDf5AKHVsIWEuJg33JVzGGA0Ui0bB1iQPkC1BI38obh75nDoIIm6ax7plvpY4IuHiEiQotvaTWjEr/Lyv1ySQHOIfcxm7VKJgrVoKyESbpyomJqvIZC7uLIZINtrvFR9p4+CeiLYGEy06PaYsdCeunbJMTyzbhyvqK8UN2nsqrxgqQ/AWOwTstgMRQlMwrkq9MI0K7yi+EBQzGAO2G/X9NyXjaUNx1SEr5VGItdhtq7KCY8CIUF73SGEpd4ZhaWmeKc8ENCMZjnROxzd6Lj69XFNcFoRKGBcKsPLC8PyPOeoL2L0rNE+Ua+IbRygMyyeeUGiJU96o9mux19d3A7GjoBoj/ZFmQ/vzAiOweK7+PcNaTcWbZxF7cIbqSslhflFge0/7QEFdO2oL2S+sSlh6zZpsH0tOL1kXVapSj8dejwMStV8S7lQulYqPgrxmSHiJEgqTiRaDdltMJJ8x+iy7mHG2uWLa3CoWT65zsryieynwfNoz/01Gcpcr6ZDASzVpav9mkBUm1Sp7yIYFMuAmrfb/ZdRa2p4cIZttbgNnXVhQKnzpdHe12mKtnwY27zpOvj1SNi7aWEJwBrcLelIag931VLue3cNTLWEq5V9s3qlZPh7wC0t5NSj79V5NKA3FosDXhup6xG0HhodH2i7tRsZlqRhDH6KDn17soVQQdDhSJXs7WEIT2Lx/TPVcj1F5rW3jUMDuHb1Q2weqbA962+KJ2ooCGTQsNzqg51o9Hn5h8QujwarW1xer4G+5Fdozy7DUEqm5DLT3ShaPO3xlcfuR0DiGI0coVCLRjkJ/pDKSxfUYWbmpRLKUm07pkclqYqGKajrj5Al1QfEiZinO6aCltZmjklvAIhpEUms3KeIl+cRIKks2JbJPIG6Ruzbz8yLbenqPiQZf2RXyrtabQHLHK05V5iwk6oxgonnSXLt1JlRka23FiveHgkYy+y6i/rBAMm3Ku1FUTTODZ/H9K8Kqprqw9KcVi2dKQ188G0ki0bYPFC92egKLYCL34Pg7Fe66pX9rTbEZ6E8qlSDcB4bjkvbMMUb2qh0s4gzNhaVoFYcpdh6/iF2P0uD2ghTa2gXYvqM8EfvcMPTaObK9oz8Jalx+4uPFKDAY3D7Nmhi6hyPFtcPtlS8yLlRrRax6wzSXZB1UsQbrhe3bJT5KMjYXgXIbaC61POlXFhOJa9Zrl6Y/rRgXBnF1LHM8tlctkmIzUD9XsNMvFW9JwQFA53Qm06v2YY0d4hUmUOw9+6+dUmwG7BiQ2ikhznvELrBXccBv2ejn4cNEm4+1g7RdJj4a5wg3m2kINAka5XaxOdDIMWJ0U0oGX3dY7hj5/EeSO4SeP82SaTR7LiCDBgdJbbxsmxgUD5npkEgf8YZ+OKA0K04iec4CYzWIzJz7ZKsnc2o7msEjhaV5smP1Ya9dmlIFjW3ERaQqkMIpQS1mTO5iCz5Q3Kj8n/VCqA2bd5UI5aNgcn+so/vjQsWTuiNLdTNQXexxuxHbehaPWgVCt9G1bhRW7wurDwL1C6F5bqIeCoRl5Np47XbYrcOMBr8MVJd6wbutxb/d0Z959g80OPZrVcoXZyh2KqLcnaofcnfiqK+8yhDshWIXtPszROmDSIhzgyruA9QXHa4XXK+fnx0COIPtvJaMzhIWBcVWP8/iRcoqHGFVEwrLeNwwHFeasRjD7n5Bf+wYl45hbenPKvqTivG4VlX5avoccVaDSBfnYNLF7qziKZk3YrMEAZB5STLcwlT8YbkjqaPEDOT/tEs+wddrvF6TjMRkHoDuCpN1AHCAmeQMZRyxkSeSW8JRz9XEWYqDV1hVUxm0XCFViRm9iuG08cTzIatxVU8vwRjcC33OGpTElOY2IkNWjlaEk6UyNRcFLl4kZvAUvcdtB8qjSqUSd5Ok4f6hob4UqhuheeEJhWX/7ipfwN2xY/F0INSW5sIzriyrx4H9uZLNmmdC/4cMQyMU107bvN4gVaD+stL3y8JzE46xraH++g1fO3/Bdxb36LdrQmWorg3bE0u5Fa6/XtA8U93W8lqzpVBqyWO8sPlywf6BUSnIEk5/3yt5rxN8pYFl946quwenGY2v9XlsYRmOHMPCsnw80L+9wHWB0lncttfMDjDrSrtYFoZVpTamXtjfs2zfsnmievkkArlxAtom1bxExR9HSJPAyee3iJyioO6Mea4qyR1YC9JHlbSZMlpS4AuCsbNy5jXt2hhjHPBrwPsi8s/cpZLhD1uvSSCJy5pMRZ6bZc0p88nrJgnaAJqZzMoes1zoCbBcwG4/XfzDCKfHGuCbMgsIm7pCmgLTj/ijBtur0lfxfKvvqXTYTQwio0eOFOhjtdDMpHSEWgNf+3BBdTXQnZcUW93JqssOvyjpj/XEL7eqA1LdaPepPXe43hIKw/6eBVE5w6LTSdx0oo0Lq9qxhU7f1s/BBM1wTNDZG/NwIARDUQSuL1aU5y3+0YLd0xXf9pa6GtmLArbdicnareWNkteKjT4/QQf36qvA7oFlXBqGtcSsBPqVWnJcf0VLMTsYyr3gzwvKfaCvncpUDjpxVmw89TPl7Sw+aBVoFYkMWMWl3IsdYVkTnKW68fTHLsoqSCwHoe60o1R02vmpLvYa1Pdt1pTJpLRUKpvJ3nVOFzBVpaVKEsfyHknEttnmlZjSgE4u35WnTTpf77Zr868B3wKO4++/yB0qGX7Uem0CyQEuEuSwe5PavTObTmAyV4ru9DotqmbWJILZ2/do315R7MdocB0Yly63JE1SXVxZ6iuvwKKHxeM9w4M1bh+VxTstjcZ3z7G9hzFg2x6/1E6Lrx3FZqAUQATXBg1SXvCrkt3DitUjz9gY2nPtWnQnNgOo/alejMNK2L9lGI4C+4eW6ko7Lev3NTMot1C0nmGlz+FLg28Mbq/2Ed11hS8C/c7BIuAfLWieWPWGYYF5VFDto/RAHP5ze7RbE7SFm8DRYusZjpXhm4JfsVex6+PvtHTnVTRSVzq9r7Vk6o4czaW+R1XcD1A4deY7dtRXJXYQqqtB2cGxNY6z2G6g/lDFkepljZSO7v6CYj9GTEQV6YyAbQe16kwDgZvdbMRCuSSSaPXJMKvtcgZyIAUQZ7PyKEYC7INk7MRYc+dBRF/kbp7GGPNl4L+Bzrn9T+LNP88dKRn+sPWaBJJZ4HgFiHVAk080+LIk+dYQp0WlHzCrhfIOAFnGOtrBsFKi2PXXasbG4Bc63m/7iQ7enqnAjw7KKetyOC4pdp7qagdlgbtpGc6XGC/svrqift5jrFE8wBiK7aAude1If1rjF47urECc0Z290fH/9p5QXcJwBN27PW454q8rzGio39nQP1pR/v0bNu8dYUZDubUUO227ulYlCopt0F1blLi2/BDGpVXV9t4QjEoFaMdHWLxfxItRsZX6RbTEiDyZxA/Rgw79SZE1VYIziNFZIzsAooEmlNotac+VyWq9/m1/z+nzRpNxtw+EytJ4LYlCpRR6Y41ySJoCc73PXRsAe7MnHC1YfOeFYhrdQDheYnYdsqwx+17LUhEVRJp3UjKxMXb5UqclgfQw6ftmOUZzOIEOLwWOO89I5BOVNveNMb82+/2Xoq5PWv9b4H8GHM1uu1slw49Yr0kgMZGNKhMgBhk9z8zWpKdZKXiZ69gkgGTQLCTNvBSWUDtsFxjXepsvozL8oPKHcqLaH67Vk6O+FG1vnjv6tWH9yLM7rvDNOQCu9RRX2qpctnFeZPDYvdMTfNVA5D3YUehrw+6Bo9wK+weHthbdfWF8p8O8qJDGU521fPneJe1Y8PWfu+BFu2BTr5CV0J1VgKEbLOPbGlS6I0u5V8Ja0apmSHVl8ZV6/7qn2tlI6u+IqLyjg82XhZPfQ6d6g2RN1lAYHV40FtcpwNpcKOA7Li2uVaU4cdrmLa89flGyapWkNy6cdqNeeM3IRAlovrYUrce1mvWoRUXA7gZsNyBFFHIeY6lRqJiVvWECzscR2/VI22KeKrgqiXQYZ24ya9nPzo8kv5nwjiT+7L3ykSDTDfRDdhhTTARINKBkM7b0XHe1Pn4geSYi/+Cr/mCM+WeAJyLy68aYf+xjPNerouGPnRu9Jl0b9MObT2NCFLWxL90vZS8Tl8QoySxJJhoziQ7Hk9a1eowWzwNm1B0RFFcQA/2xZKWy9sxGwFADz9gYzKgX1bhyDGdqYWGv96qwtqrVm/Z4oWZbN5GbIMLyww436M4fKtUKGVcBKXVmRvYOaTzGCv2u5Nvfe8j1vuH3n97nw0dnlMcdx3+rwvbkgbqkM1JtAnYQ2tPUTYHmqZp4lVuy9cTiItBcCK5XbVixsHxsOHpvxDeW8lr1UENl2b1VqkyBg+7Y0Z5buhNHH5modghYL4xLh18W+KaIpDqj9HivrWs7Kj4SnIlDiKrpOq400GDQ4chGFfqzT1FyEEj4VwoicwazMdO8DeRuHzCTnJCJVzRfyd9mHKdzJ3UKM0fJTUEkBap4++SPczf1SCKk3YGw0Z8E/pvGmO8CfxX4x40x/1eikiHAj6Nk+HHXa5GRpNCY2YRRl0RvNNOQXpqviTKKMtthstiMiNbMhdLWi02vJ/C6VJHiXqhv9KnFKVhptugcSydZ/6O+DogDRG0YhmO1nxyOHMZburfXSKSOL9/fRU7EHn9UKwdh8Noq9sL6/ZHdw4Kj72pWUj+3jCtVlrdbBw879XIZLOVRT9eWrJYdQ18gwbL5RlCXvdJgR5VRdB7GxuaBvt0DHdv3jdLjl0/0/VfXQS/u2mCCodwaVo8la5yE0nD9zYV69pZKPhNr6I/0WLtO2L5jKTfaPVk99rg24Hr1r0lBRCyYUQOFrwzBWcZGZ2zGUi08i73NpZMYDRq21Qlg00Zf4G6YSlzvkWHA2ChmBIdMZj9OczVz4liax0oaNbfL4nS/lJkM40FZnXRtXlKKj0El4SbcUVJi7iAoicifA/4cQMxI/qci8t8xxvwvUQXDv8DLSob/ljHmL6Jg6ydTMry1XotAIqCBIIyTNkSaiYBpt0pC0CndLIqoORFUq7MotLMS51dCU2CSOdMQlFG6C3RnBUWrNPDOk4V8QgHEXXxYKavT9YKMRqd5G0t37HC1MDYF5V47D92DBjEmA4H2Wqdhzai/X32j1GzkSlg+CezesoTK0J147L0Ovym1A954huuK4mhgVfccNR11MRIeGr79nbfAKD+kvDHYXvGIUKDm5gPs76u2SdEqltG88EgB9QuVQuzOlGDmq5jBrB2+JJLIlB5vBw0iwzGU14q9FPtYEg7a6jXeUl31dOc15c0kBWC9UFwM9KeFOusNas1R7vTYh8JgRsVMMks0TMEI55B1of698XM3ib6ezc5GHRS0Tj+4FETyZG+8X/StMVU1abCmdm/s3uSAEc+/PGEcN6g8gU7CRkIuv+9Mj+QnzxG5UyXDj1qvRSDJa+5rk02xwixjsQdgK8TW3KKZHOiHAbNYYNpep8lLp3MtziBGVbyqG4+vlKFZFtp2HZba9myeB7oTQ3WjbUdx0B0bTFBBIhOmrKVfqaHUEA2tt+82rL+306AyBkJT4BeFaqce6YVpB8VkhrVgRsNq2bE43XDa7Pn+xRnN6YAPll1f8s2z5wA826/5ylefcdPW3GwW9D8zcPODNeKEYmepXliaC+H4e1HF7EoDx/5+RblRMpgR0Q5VZRgWBhfNr4aViRgJ7N5y2F4YjpUs1z5QFz7Xq3Sj6wyutYgVhtUyEthUBkGDKsixlkBuHxDjYidMM5HUDg/OUtz0MfAq98MkD6JOW8QyH3dIGUg6R8aIhWT5iIh1pA0nTeqKRDLiMBsCjddKzEhkrmcSg5UEweAJs45Nut02NaHtDrlNn3LdtUKaiPwNtDuDiDznzpQMP3q9JoEk4SLuEGi15oBHYppaA0fcZfJJMY5QR1HhIdkg6HCXGdVJz3jB9l4H6Qple/o6tS2jlUNQnkSotMYv2kB3ogSwNP4eCr3w1O9Xd3YjOqMi1uCXJW7TEZxyS8Qa1u+HLJW4fxCB0jOUrl6OPL9cc1ksOF3v+dnzx3z35pzeO6wRWl/yJ+69x68++Rr/xbfe4/vrMx5dH1G+u6WuRra7msEvGJcG1yqbFmD3UGdcxGrdX260xBnj3EwSigbFbfojLe1U+kCzp1Ba+mOdOK6u4+RxqeWTa6G9b7NcQX9kWD1Rk+9QGXzklvvGEmKJNa4c9eWA3as/kMpQxixhiK3YMo4rQO7IqR6JnwSq2g7juEUVsBM1ILd7q4yT5Bmt+HtS3oOU2doJeI2cJIMHCQft39wyvkODrNedtfpx1msSSMzh7jBz25sviTT30HUv4SYM/TTODmqBEDkloXK66+4H+nsLXKcZSSg03V48DwwLkwfkmueBolUgc/VoUOB1UL5E/bwllDpabwf1bzGDV3e4biAsK/xxQygdYgxuP1LsHf2RArih0Auy3Bq6e56LyzVlNfLO2RUn1R4vhp85fsbC9XSh4B89+R1KM/IPrb/Dqdvynyz+ELuziv/P9/8oi2rAGuGmWmCtqqetPlTW5/q9XrsulWX7dglHjmFhKHeB7amL1hKA6P+sDoCGUCmvJJTq9etrE71/VRO2O1MyXNFGbdjC0J0LzVOl+i+fjUqZFzAIY62eNyEeb4LyP2gHTOTmmC4ymeMAXWadDgPZCWC1VHHnNL09H+Sc843SeARMOBro+TX3SEqUAmsm7om+ODIctnhzEImb3J0T0n4KAslr07VJTnoZgU/AV+QD5N5+BMjyyZKU0oZxuq0skWFQwHXXKnFp0Eld4wVfqnLZ2Gjdv32obdRhrbu1G/RCGlaq9uX6QLEb41BdiR0DxVWr2cxRw3i+wq/KHEBC7fALh+1V8NmO2qIVFzkcLYplbC1+W9BuKi73Db9/cZ/Wl3z75h6Xw4I/efy7fK+7zyAFD9w1gxR8qbzkYXXNH7n/hFXVM3gHVo+B+uaoiZUJMKwLfGWprwKui4SzACffHWleBJrnmk0oe1QDRfIENgHGRt8rQHcuDEeB3dvxeMfZHjXSmszS9/eUe+Ibqzan+8D6e1vq5x2r9/cUN1FOsRsx/aDyiv2gpemsy5I3imiTqUOWE3M1ixwFidKIpZLPgORxZB/en/R9b1Hfs4BWei2YOj9JojO9F+cmbCQGlJwNf9oVCYwf6+s1Xq9JRhIPUga+Zi28WXmTDLNMmo1ILMQURPoBGPRDLwr1ga1LDSJWBYjdfkRiO7i88YTKItbGC0NT9+7E6ACcQCi0m1GU+p6aR1t9TWNwFxu1kjxa4o8b7T54wT3e4u8d0T7Qkmr7Vkm10QvZdUL7IHZMCsEuR6wV+rHgZNHy+y/ucVRr+vzXL3+Wo6Lly9Vz/rP917kYVzgCu1DxZHdEEMPuuqHaaibRnscRf3FUhTJ0xcHmS04v9AS0jhKd8VSusT/Vx5sxclvWQTOKUcu4UAsEkEWAAG3lMOcdYbSUH1Y4Y7j5KlQ3UYm+M5g+qKL8GGnwgN0P6n2T8IwQtOtCoVlHWebvB5YTGQTtkU6B9CRSldTywnVsxSWBK2sJj55k+c6sdRPPM1MUCHGTSqV1EnfO3JJY0rjII0mi0sSs5M66NnfzPJ/lek0CCVO/fp6RzAWOjIHSTdO/MYhoUDFQVfr4+LsMgxpmtb16qzQVxb4jHC20PHKW6mLP/p1VHoVfPvH0xwoMul4YF5ZyE7s9+zj5u+sIRw3u0QuSBoa52VFcbwnnR5jBs//mPcrNwLjQ7ADRSdv2nlpi9icBKQP2aECC4Q+/+5hNX/NHTp/w5eYFv715i+t+wVuLawCOXcuvPP4mP3f8Id/d3eOiW9KOBRfXS6Sz2E41WaurGRB8ZLEjlNtAfa2YDjAxWL3QLq0O3Amq2XoshJOBs3sbukEtF8rCc3O9oFn2OBfY72pCHaiakfb5guHeyBAM9VPdwceF0ujtRhgXjmLnteu17bUEXDe4Z9dIGr5cqHVI5o2UhZqWWYugQYBB0HZaBOHjBZ3mZ9JzQSSWDWTeiWFOTJOc9ebb5p2afCreKmMSj2R2+11hJIlH8nlfr0kgmU/3ahCZrAB0x5i3fQEF6UBPIi+YwsZUNWqwzob6zHqpJ+rMfb681kG6+llLbQ3DuqR63rJ4bBlXJUaE5lk8mXzA7nrEOUzX43bt9FyJDVlXWV7AjhJ1RcAOwnhmWH/gGRcuD+7ZziJS8fAbzyls4B95+Ps8rK75D57+LP/V828ziOO7+3s8qDb8v57/CZZFz2+8+ApPtyu6oWT7YoHZO+rnLos8F512ZuqrqHQf267VdVBleVFMqF/re+iP9aLvTw39vYBYwWwL7n11xweXx7x1csPgHd84u+B3nj7gbLknBENzPPDWesP37Bn799eYoCVOsdPuzuqRktnWH3SEwmIHr4OPN3ts+kxid8Zs9vi3TnFPu6k7E8LEG5kJex8IgCemM0ytXSJmEn1oElfkIOuASWEvBZH5JG8OEH7ijJSFlk3WZYYrclft39e/bPk46zUJJMT0s8ipqKT0c27YnEqdBLBGhD+77IWAcZUi+oApkids0LKnKnNtjgg2WkZKVVBvVCENEYqrDhulBUw3qLp5CGqKlbKhcWo3UlfaIdq14CtKwK8qqivDcORoXgibd5SoVW6heWJp3/FKi3eeb66f8V57yvf35/zzb/86lsDf2n6Vry4u+JVv/xf4yuklPlie7ZZcPDuiXvVUj0rcXrEQ66G+CtSXSkIrdp5QKcajk7MRmDRx47b6s+tUF2X1vlBfOHZfEkIpfOc338UvA9/b1tSLgat9w/6y4YO+oKpGrm+WiBi8txx99Zrr5yvq71b4Cqqb+NyDYDuPu+6xmzih6z2mD5p1jCpKRVFgv/cYySJUEUBPVhT2kEiGNdCP0ch9poqXzoU8O2MnlbPsb5Oy3Y+oJYzFlAWh7TQTMRZkiK3lGMzmHKa7MMjiC5KRGGO+AvxfgLdRqalfEpH/nTHmfwH8D4Cn8a5/XkT+3fiYT6hzIBMzdUYImqeb82lg8VEhPtbCgAaRaGegJ19kQxqrwkVxpzP9oJOi6cI3JksmSlOoDmsUGaaYpo3NzXYiOc10XhnGqFw+gXN+WekQn4zqhucM6w8C3bHyOcaFpX7sGFeW98ozPnh2yh//yntUzvOb26+y9yX/5Olv8cFwxj/51b/HD/ZnfP/6jP/y29/nb/pv0HUl40KncRmU82F7A2LjgF6pOiErDSbKYFV8JFSaNRQ7wEB3rqzU9ssDxfMC1xr6c09x7eDK0TcVcmUpS6HYVey/OuCuHZt3Be8t3XUNVkcMqkuT7USra0+oHa6LzFJrNZjXlXZn0sBc36vnLmhXJgURN5UbZrFANpv4OYuawe/2k9h3VU0KZ7OJ3rmYeM405p40t8sTCRkHkZjdZh6JsbFbo8Br6O8oisAXpv07An9GRH7DGHME/Lox5t+Pf/vfiMj/an7nH0/nYKpDdcfQk0jSxZnKnUhIk3Gcpa0RMynLzCkwZa0nJTNEvogAnbGqKWKNsiYXjQaJRYN7fDmRl2I2JEGUixKfDxeznKKYJPLqKk+lhlVDcbkDY7AiuK2jfmYZj2olex05tZc4NtjBsG8q/MLzdz74EufHO3734j7vHl/zy8N/iaftmqOy5XtXZ5Qu8O/92h+nvHLUzw12Cc0zSOJCialrohDQuDRKgtvAzTc9Uge++c3HPLlZ09jAsu4JYqid57Te84PrE+yXoBsd/bMVvrE0Ty32Uin3RnTg8Pi3SjU0/3BFdans2gRS20FoLj31807f136YZp4Sx2P0U2lSFLBQCQjZ7iP/J/rUxKE7E9AgEmdnxFrMvp3o77NsdC7NmVTxUok8tzoxxhyUOgQ1EM9dn3RWOjfxSIJM+MjckOsO1hciI4kjyGkM+cYY8y1++Ljxz/MJdQ5ytTnPRpJad9JwZUY0uq38nRiNqbWXgogxU1obn980tYpBB9FPMFlDzhmQ1k674zgiIWZKyeoA9PadZC0Ue7NHmgq7U2X58cFR1i0xg6e46bLYsRinvI0rrTPGpaG3wgsDJ2s1i3q8P8Ia4e89eYtuX7I+asFK3ryaZ2BHbVm7DroTqz41A7heSXbjQiAYFu9s2D1dsRtKmnLkfLHjuq/5Y+eP6IKjtp7LbkHnNXDfWM1wpIDmw9gWH6C+EHyturG+MVoixUC2eO6jvIEO9olBMSPnoHRaNh4tcc+usp1E+3Nfpr1fsv7uFtsuMT94rNPbPmAK/dzzfMt8klcE+l4/y1QCzT1nUiaTg8hMixUmway0UTiX8bS5nOIrW7wxiNzZrI0A/vMfST4RRmKM+TrwDwC/ik4b/o+NMf9dVNrtz0QJt0+sc5APY2qpJXR9zlxEwTa7XBJ2uyn4zAMKKF0+zk2Y0xNVjW+a+N9qEFGBmxaKKAINGd03VaXMyapSn5XokWPK+JizE+U+lIWm19FfhX5Qo/GF4i5uqzoZ+y+vdb4knizdmcoiVhuh2AdWj7XTMbxXMqwrrk9WXCzugxX8kcfdONw7e/Z/95RKwO1M5KLoBV5dKz+l3CVBIlU8U80RDQb7JyvK044nz445Pd1ysV9ijPD/e+/rDH3BatlxdblE9gXVM0fTG+pLKLbC+sPo4wOqvdqoxcRwVFLsdShxOCoorwf6s4rFhzonY/fqLjieLXCbnv5MTcfDmarC775xysUfLdm9E9ifr7n3Wy3V8wXSdZhUe8QLPJHPzHqlxzwEJaj1Q/wsD/kn0vckS8404DnZSxCnf1NAIGfDc62R+XRvLpkBCRGcvSuwlS9IRpKWMWYN/DLwp0Xk2hjzfwD+DTQO/BvA/xr473PQgsnrpUNljPkF4BcAGqI95hxQy5kJB2Br2O0iCu/zSTLnlDCM6lHTrDS1bfTCpiqRRkFR4pyGHK8UJO16rcuDZiWcrPX33V7xlO1egdrke1sWhJMV1hhtJxtDaKJz3A+e6WMGz3i6oLrosYPHr0qCs6pGH4T+pGDxrGdYq/bHsFR3PNca2vsqORBKS3lj4bdXSj2vBd8Iqw+iqdVGg8ew1CnbUChvRIxS8JePhKu/TzDLEWOFr7z1AoCHyxt+/btfxVrB947rRwvlfgR12qtfCM0LVY0vrwcl6W17BZy7kdCULH5wQ//WCjuEHGiax3uksNgIXBvvKS5bBZ6ftxjvMfuesF6w+tYTrr/2Dmd/13DynU5nbdKkbhCwAdlroMBYTBUz1ASoRkV4mZc5IShG4hw2uuVlyU6Y2NMJf0t0gQjCGvwUTJKBGuQOjylKkIi73aGK/Bema2OMKdEg8ldE5FcAROTx7O//R+Dfib9+LJ2DqOz0SwDH9nwW/m0mBs39bbIb/Pw5knVnxElMU+vjjNGdrVQPFOpq+rCqkuAMZmeUKFU6wskS9/yG8PYZ7kJBvXC8hNUCCossKkJd4FclxWX0X0mboAh+VTEute17/Y9+lbO/dUH/cIXtlVa/v79k8cEWUxW4bYeUjvqJ1wvyxZ7QlBR7nSDePXCc/T0NELuHSpzrznVex/WTAFH9PLZWx8CwVrDSLxzdibrqSQlXPwOhEk5Od+y7ktqN7MeS71zeY7Hs2e8quC6onylIW+zg+PtqQdpcDJhROyjupsW0KoVIgOJqhzQ11eMtOEMxxGG8RYndDDE7ix2vuqS83KgJuNWSxm12yOkRD//TG0w7qkeNDxMOFTsmJjJXRSa91WzDanUC2BSFOuxZmzlF4qORWuIaVeXkHwwT4TGV0rNz7naGcnBezuZu3mAkh+tHHg2jvda/BHxLRP7i7PYvze723wL+Tvz5rwH/ojGmNsZ8g4+rc5CmL2eTv0nxKn9osfWbh6rShGdZKi0+niymKBTpH+OY+G6vmYqP+qJ1iawXjGcLVTjrVb/EPb7M1hL2aqtBxKn487iuMENQ9bPYft794ftIXTKuS7AGt9PdezhfYlvP/mGt9glrZdVKZRnPFmy/tmb/7hEvfnbNcE8V6G0X9U8MVDeqLuY6baPWL1QkyAQlky3e2+CuWxVd3vRqTt567KB/942SzIwHc9Zz+XTN2dGORTFwsV3y7PExu+sGntW4nUUKWH4olDdCd6QiRGPjcO2ofr/7nrCupzOmLPBHNVI7/Kqif7jGr2tM5/U4J7OqUoWhGf1EMkzzME8vcO8/y0Ekg6UA0adoLuidPv88j+W9DvrFDSLhYXmUIgpcJUa0qeKmEgHUuSiSpOwkfkm2MZmJRaOYieSA47iTJZ/g6zVeHycj+ZPAvwz8bWPMb8bb/jzwLxlj/gT6L34X+B8CP77OQdot0j3ngSWtJAIdZicLaAsx7RTt1MJTnENbu+z2GGcj18PDoqF6fq2yA8tFdHCLcn8hwOix17vovGep94OWNU4zmbDQQzeuShXt2Xl2X14iFvZv1Sye9OzPLdVWMYzdlxrqyxE7BsqNYgvnv3VD+3BBOQrlZqB5/5rlUm0zFu8NjMcNu3cajn9/R6gdxXbAPblE6gozespHPeIsRT9A4XBtRbEp6O5VEArsCNu+wbWGR/6cZy8eKp4yavlz729rcMII5S5QXY3abdlNlhr6NarqmzHYm52yhJ9v4u/KrTEpGIAGhqqEfTsN4DkHweuxjtq6crNR0HSxiINzRoPITB1tYqDKBJZGQ7J51iHOoVaAg5YyVTlRA+BABPrgfIMDXGQ+kDeXEJiDrGqodTftXwMZP/s8r4/TtfmbvBr3+Hd/yGM+vc5Brl/TvI1T4aN5WximDk76mcMTQPoBu1oizmldXRSZpGaKQsue3R6ub/R559T8tJMlDsSihj7g76/1Ioqkq/Fhxf7cUV+rSHLdakZhRz1JjDH4yuArLUXKS9Wz8IuScjsihUGC1QniTZd32XLwnH5wSTheUlxsFRuI5DhVgPNKrjNqJmW3yg6tn/dU1yPtecnqibB9y2GCehM3z3R2qLkIlPsQzcCF8lJLruLZZuLchNTVigE2KYmJKODsrAa1zU4/q5RBFE7xpViSmtUyHj/Fk4yzMUtxkdIepR/m4s3OasfLTbMvmQxW1xoYEi4i0dM5smGNs5mjksvjBOKn86auCfs2097nCvF5MO9WaZMYsdoWNtyV+fdPg9Pe68FsFabs47ZSFSinIKlWuSqDX5JUrPo+a0rMmaeJFyC7vaa9vfI/8t/j2Dp1HYfD5KBVrASliWxm2g7xgXLfIccr7LIiVI5iFyj2MCwti4uR6qJn+5UF6/d7tU1IkoSlVWanNQzHlYoiNxYM2NISTkoW7wf8usZ2I3hBTha46zZf3Po+er2IKy2pKAst1wqrWitBEIH193eE0tI8huFESXK2j4bmPmjQShydtptKkpiRYa3OwkS2qQF1F2zKzNXQgJJa9XqbLBtM16s9iDFqJOaDXuiJuZrc8FwF6XOKvA8AYlflgNuReCNdN20EItjEgIYMvKcgggRkDJqdDhPWNn/elzISY3mptyshey7JONx9afM5X69HIDFkEtJLK5+kCRCLF1O0F5BXoDyy32diUmY9YrJqlomBSmSMu5cGH2k7zHIJwefuTxL7xTmki9yFtsOEgH024sqSYf0lFj+4pv3yMbYPuJuO1fs6qCZOdy6/KhlWBdWogsjDSlmu1dVI9eiG/q0jXCeERUmoLOOq0Yu+9/h1rart206NupwluQJSOBWgbgr6kwrjhf60oNgF2vtlpqvXz3qksBTXrfIWnEFqDUB21xNOVhDQDCcOzkld6e9VCcOIP1lp52XXRU5NQBa6+2cFeGfh8gY5WWP2nR6zeF9TFJFwFlnDyV61rLTsSRhKCujOTR2ZWOZMtPdpbioRFFXnNRxQ49PYRaa9z86xPIQ378jgDvgjLw3qScgUhLtZb2Zt7nBNnJED4+8EiM3+JmOYAkJiK853h9QKjmluokxLUiGPdo2mqTM1O++4VRXbjRUSdsodSaPsseSRMQ5z9YN2iURY/p0PwFqav3Wp5ZUx2E2FPzvCvdiy/5n7iFXfX9da+tNCg0Zp6I8L9g/uUV8MjOuC6lK9chEoNuohbIaAe7HVnTxrrhRqOxozieLDF4g5x0b91FBbumMbGa9Cd69S8erzBWIN1bMdw/kS146EqqB/sKC4GbCtA2ewo5+sNLd78AE3jBnEZPQaaLoeqas8kUsfW6X7Duk0kDD0E/g97DVYp+MelcxkmNqyxtkJW0mBx7k8qIm1mdmcMBTjXN5k5jYS6RKdlyVzoSL9/nKWYZyWXaEf8n3SursgEl/r8x9HXpdAkhiHfrKkyOCrzztBCiyJD3DgvneL2Zixk5kEn7rLz7ROElNVBNM02n6sSqRU+nt6PrNcEK5vJmJbP8CiQXa7eBEMmLLU529bzOmJBqaIZRTbUfkiwP5hqTKPR4aiFW7eLaivA/uHJcNSsQ4jet9QWcwo2G5Uuch+hD7Sv9uYjawWhGWFdZbi2Q1UJfU+mnTvlpGfIfRnVZSRtBR7j19WlM92qrhfF9SP0iwLSGQV5ywjAaEwTdoao0EMJnOqxP8R0XJSQuTxhEz8w1jk5kYDSwr2seTJ071zE/n0mabfY1aSXRajCrzgJ7Hm1J2Zt2lT2WKNCk/PW70pQ51hInNAdR5ETFHq7Xcqtfj5jySvSSBJszR6UmTk/FabLi2dp4nueolLMlMDz94mCduoKt0dZyri0vfaMSg13Zb9XjOQpsJcb3N6n7AVU6veSU7Vuy4GOG0HKmAYMKvV9F/1I2HVMKwLilb9hM0u0J6qVOP+nlp2dsfact18Wf9PXytjdVja6F5X4XZRlhA0E/GaGZh9h9t3+PvHOvBiDHjBblvKF230kDGsfn93gL2YEDBXMXjIQksNH7SrldzuUqliYuaWcCQTM7MkGJQu9mCnjkvUW+XhOTx6OnVcJAaLSB7U7CFOWqcMY15aJPwjDv7l14tZYnYeSFIBc25IauEmfOSW7kgGb7OAkculziuXBAVtRV7dfvhxlnxBujZ/MCuWDTGIZN+Q2yvPzbhpaA+mdt8YMY+UhSQP1wQEQiQoVXpbWeiOGWc/AHh6oYEl6sLKTlmtzLKTKYiEzGswdY3cbJCVCieZXgfWTOlYvHdDaArsotAZlKCTuMkrd/uWUGzVCrO9F+OBBxMMy6dCd1qy2I/IosS0lWYRi3rq2DQVdtshVaGK7Ns9UpXYS23REhXkiut4TAuH7NuJWJXUyeYj/BH0lLaDkyN9ztVCGb7X22lid9dO07oSMEdrqCvk+kaf7/GzaZ6JeP1J0NGFssBgJkV4a2BIreKQs5R8lsz1V9M0eMRVsnd0vnPMWvy8vJ11bRzYqszzNplPklq8CTtJP6fHWYPcIT1e39jdPt1nsV6TQKIBIrffEmkIDsuaRBSKK5kg4UwuXYAptY1Bw9TVpG8xswRN8zOmjDu5tXrfJAtQlvDusfqspLbrssFc3miw2e2mtHgcMWcnSmgbPbLZwtsPVOZx9FgR6md7Qqmt2FBZXG/ZP7Csf6DzMb6C8ibKAqRN3kG9CdrlEehPz1Sw+sMtUqpK/e6dhmKvx6vcjJRdjzTlFHCieI40lQK0YwzCaXwgHTcbO2c+5MzCGAM7nUvCWsVLolhQXtGkKpeHfUz9xWugXiwm/RbvI/gdP8MhmsDPs5CmQW5uJh5I0lydvc+sBm9ngSU9PgkbyWF2YatSafNzLshM9Hne4j0od+Yckrukxqf3+6a0ueM1JwzFn1PKecArAZJBdAbb0kxF6rQkOnVs95qy1Mwi4RvbnQaKYdCTP2VAxmBWyzihqpmMHC0Zz1eUH7xguL+m9EFNrHc7HeIbPfLiUp/v6QXZ/e/Zi9ip0A6F23eYowW2V/8Zd3+B9WW03bSEChWk7ifJwmGp6uubd616z7wrlFcG/w+d4FpYPNHhPTs6yo2wsAa3XWmG4qxaidZHWsr0I/50pdqy2zghve/0Ik//fzfAcqHHSUSPH7GcSgElMUtFtMMVVczScCNNrcG0i8zTqpxwgNQ2bpqZhshwiHm0LSTMKXXNQANI7OLIbOPIZmpxnCJ1aOYzMuIV69HuTDgEUGcl9IFy2iuWKqR92hP91noTSO5qycusQ88EiCWNkrRSRpFYjQn9T1OiyYEvCRR1HVJqZiKJxMWU8uYL3wdti15vlEQFWcnLbfTCLJ9r92Q4bSjKd9UoewzYsiBUJfbZqBfibq8XoffQh5wlmXbAeEFKR/WinY6ALXQu5xykMJQbUdvMBp2iXWlJVGzUkziIjvO7Acy1UHSqM1vs9MIO60bH+EsNEP1xAwLlTewEmYXu6lU5qbs5q4ONEEf4jQLP+8jzOFpqx6R0mG2L2etQHDGLk2FQrOj8RGn/jy+mEYU4IQ3kyWqMRe1NFecwdaVG4ScnyHanow7GIK12gBLlPWefMTNJFIBpnCJxjqbAMHVpwuHmNDPGeuUcTaLJJ26KsTkY3ZmMwB0R2z7L9ZoEEmaMwcM6F5g+sDTCnTKTFETSSoLCCc+o60m3NcSTsSw0G/E+8yMA/dlFsNBFeQFjtHvTD9gg+regJUL1dIsUVk2nThvsLrYniwIpHByt9AKJICaF02nktJv2A3YYaTbKEK2el+AM1VVDqCztqaO+1GDSn5DFg5rnasxVbJQf0p4ZVo8C9ZWnP3YK3vWjiiytFxACoSkpL1rG05r+pKK66glNoZwSmISXRZSlWld6IR8fwYsrLYHKaBZuDMNZQ1E6wsNj3FWLvbzRLKTvNXg+ea4XZddN2cJuH19iVoKkUiVmiub4KJdGeWq7cNEsS7kzWJNlA3JA6PtZKzd2Z+aaJDNV+Pn8TAZh0+mWSp15WR2/z8cw7jIjMcib0uYuVxLVNdZMSP2r7zgRxNIAV5YbCBkbmbeDTRwgy1iAxI5DXSuptiy1C4GCqdL1ehEFlRbQLk2vF1jcsRlGzACIUD6Pj921+vcQFJ8oi0hqC7m8SUSqRDE3mx2yqHFtT1g2VFe92nxWJjJlA24wdGdKbBuXOpCH0WzE9aocP6wt9cVAedNj9j1UJXbXaYv2wQnuaqtB0RjsriOcr/T9lIW2khslvZko+mRKnZUxpbbIE7YiZ2uqRzeEdYPdDdibrR7Sy2stRVPgTp4zEejOfJ+ZG17OJBOzNEljDoN2p1Jgyp/7gPT+QMdXwpjPH5hwkNyxSdlG4pHE0lfC1C1KMorakQlTMJGAOzkmbOL/OG+D3+X6KA3Zz9F6PQLJrYD8UhBJDMd0guRU1GYEPQ9wgYJ3QTSdTkj8skF2LaZwUNfIbq8nvQ+I7yb3Nufg3qm+JWcj21JBVkSUnh53kKRFYi83sTQLGMYcOMxmp3V/3NEl8lZMUBo5RtN64+OQ4LMXmG6FKwuq90b8vTXdWc3iSaA7VxKbWLh5V392PZReRZNcLzTvX2spJoK52mhQsxa76xneOgHALwvkXkP1omW8fwRAYYx2YhLjVA+oBs7NVoWPNztln0Y9GJtbruRu2hy3OBi4iyBrEurWgzfjikQ+ikT9EED5PjMSYbLUNE2dqfFzfdUUSFIQyRhIkgRIsoszTdeDcy5lN/H50nP4q2s9pyJ/5Ee2iD/pelPa3P16lRWieI8Bwrz0iR2cbFORTuAUcPZ7BetSvXt+glxcYVaL6UJrdFQ9abdKiBJ/afgsnsDD2RJ33WnrNeprMAYVit4PSGkZ3zpV4R9nsyC0iSVS1tnoYkcozsiYYUSaCikLZFHpc7u1iizFC9pd7Vk+3yDLmvJKpQhs7yn2jRLVfGBYFiw+1B3T3OwymCnDAG2rGrWFy0N/xi+x7YA4pxO8KOAq64VmItebiVWa2qnRl9es19rlqqOTYeLwpHZo8pi5vcMmTZA0bwMTMA75WM8HMFMWY25lAMmSM2vxRvzjQC8kDdxJQILN2e7ct/eA2ZpYsAcENo/gZhwTi4zEz/TVwejHXW9Km7taPyRbTEBXMnQGpu+pvEHbgPT9pCxuTOZIcHmtJ8Z2rxlJvNBMVUVN1siqXCqmkC5IWTYUl3vCosyj6lI7VQArLH4VNUqswa8qpZCnMsaHqU0Z0UDpeu0cLOPYfJUEiUpkrXMy4cGS4rLDXW4U6BxGzJMdWIuLLdjmSinapu2p4qQr0SdGjlY6MOdqzcCqEm62mcTlnoA5Osr6tJxoCWc2+6nFW5XKD4HYMo/HabvLs0YQ27S7HYidwE8zXZg5myjL7IaYDOCRfgo4dtoM5pgFkEulzPNJePyc/RpB1PRzbtfGbk3WF5keoGVM4o7AYRsYNAub3y8TDyea/Z2tN4HkjpZMqWUyJJrzCuYi0Nk1fqY6f3vXgljq5AnPMrIoo8pX7BAlm0fjXGarSmRkmuVCrSsWde5+GO9hQEV9nMFtB3AGd9NNMygxUIXjpRLSRq8ZknPqsxNbp2bRKD8jBArAtB3+/Jji6bW+v7ZTnGDfTjqy6WJte72YY1al/6KKNTF6DUA+YI5Waln64nrKLoLM3k8B20hdD0Ed7axVrkzqiG1mowrO5ZF/8V4Ddw4WBcZFU7LUup17/8TPhJgVJJnMjJMkXGzmGJBL3NSaTUTFRGt3DjyHUgAzLCR1WeKDp9IlnVox4NiqVHsJCdio0Sve406OtbRJk7/6T9wxIe3N0N7dLcOEb9wKInNkXcHV2SRmXAkfyYNckagm+ke4ihentXqimzE9MBLaEmclnmxBkMtrzHqF2e71Oasysj09bvCYMai04KWWFQlIDMuFjsAD/nwFZ8ssFMTglX6PBgjOT+DiCi6vVdj6eqMndxAtf1JmJQKnR/RvHWvrOCq32c5j2xGpHcXjKwVERSdgib64md9xeqwBpuuRVNY5S1hWasNROCWstb0G3cRMTcclyjUAU3BIXJKUSQwzCv1ssA7I2UfGsuLz5s8wSSHmiz0C5sOo2YZzGBPZppFwpubyU0YhQTD4lwb0Eq72SmxDIvs5BpUUUICX8JGDU9a5u8E2BL5wKvI/uWUmHsk8EzGzenuut5k4IrOVg0iiy8O04wWdUclC0XCQhmut7hRYBcUEEtaSJoOLAppa8RSnWYFrI+i328PJEf7eWsuepmBcV3qhDz63YF0/IscrBTZXSx23T6l7UcDQ6YDgdqfZENrCDudHKh8AGB/ozyt8ZXG9ww4l48Ky2itd3vRj7gqZLnogL5f6Hk+PkdWC4d4Kt1OB5lBY7HoR/+8xt4CTyFN+f4ngZx2y3WqgngcUOABTpe8nzCN9LrPxhURxn8/X3O7WyTDpg8xnsGRINhJMWWxqzQbRgAMTC3WuhvaKssTUtZZot+47795Mt5mX2safdr3BSO5y3a450yzLeKunnzgCSRw6pbAJL5n72EQTacysRQw5qxHvtebfRywhtYhdwjh04M2Upb5GAv980Is8ErBMo2WKu7qJ7WNHfbFRIDfqrBbXinMkoprE3S45A6aLM1zfkIzQ8QHZ7THbLaxWVO9dIMuGyllsPyopbjswLhq9qJ0BCuzVVpXs05Bf7LSYqxtkGKj2UcQoBJUC2LWKmYzjlMVFi0pTV5HXEbOYEAcTJWh72Dl9zFxzNWUiKfOIQUNEsmlVNgKXCeDU9uo0EjGfzj009g4HAWHy49XSJMxFi5w7wELmQSefY/2sGzM/BZ1DxvBS2zcHtbuCSe4okPwQV8xz4P8OfB2VRf1vR+sYPrkr5qvXaxNI5u7vc/AsnQT59wh+AZMMX5rDAbKua1L7SsEiDriZqpowAqdiRnl2pI4eKHH25kDLJPEhfFSrj+phxkfZw2GE4xW0fT7pzHavTNpZZ0l35qlk0LF8yES6qIOh5YRk2UfZ7bX1aS2FDzpL44VxWeLaQPvOkmLr8Y1jseuy2DUxI0M0KOGccj5As6CbjR77hB2li2botUtzfaNBIKmzGQNVoZ0oO+nBiPeTXgixhJm1gOfyABm7uEVFP/jcbyu5z/gd2omZMoUEiBrnDmZppnImGZLPJ8lnJdErgNODzCUFkRnecmekNOGl7PpTrI9yxfzvAf+hiPwFY8wvAr8I/NkfzxXz1evVAwV/0EsOgwgcnkD5dmMPMotMSEsr6U2k3TCe6DKMmEUzAXVlqWl2VepQX1PrfE1VaclRV4QHp3D/XKdZj1aY9RpOj3WHBuU8dDE4jV4lCdF2sjSqaSLHa8zRSinkKXAMvV7Qkf5tmjqWDNPEMjDNlPgQ54WsBsDdXieLty317z4ilCaf6K4dqZ/uFPBNrdWEPS1i+ZIlFqxOMQ+j/s/GIpvtxK3phymIxI5LlkRMO356v7P5KBGZgmYKYkEmPOR252227KKZPttXlCXxJNGp3blqWZ7EvtXGnTFUTdw4NMud4R0z9fjcjk6vn9+rHLSX8+veyYpg68f5+lHPJPKhiPxG/PkGSK6YPw/85Xi3vwz8s/Hnnye6YorId4DkivmJ1+uRkRgzpaezHn/q4MwZiXgyRT4T0nKNbqcUek5q837SZ03sSmN0l69KfVxTazAYRsJ6gZQWf0/nbfYPKopdoDuxNJeeUKjbXX3RUVzuJzJX1CE1+04vyKQTG8sXUps5dSySxEFifUYhJukHDXypnNi32uVJ5dtun4/b4je+pzNEizoTv0wQHdhLk82AJM/ddCFECUXpetjulHwXX8NYGyUQCxVsOiqVmAaK36QWr8ghoOrjENDMFiLZZmZR7ZRRzgJK5nN03UvlxXSOTOdFAkSnLsxs/P9VokMxiwF/WMLcClK3uzz5dd3h/3vn6+OXNveNMb82+/2Xoj/US+uWK+Zb0XoXEfnQGPMw3u0Tu2J+1Ho9AkkSpDGH4rsmBwx3OIszO1lNWcS691abEQ4JTmmFGepvzCQGvW+1PDEGjheYIIxLReaNFzbvOPYPDf1zlUxcPgkUu4LwcI09XVA+ulK9jstrDRTjGOdD/BRQEmEudYtmJ6ctClVHM0aBzqThMetKmNJNlPGEzyyXmqlsd3D/bCq1NjtYrzRIA7JeYv1KA86MK4MPGohEZ4kSazTjGGWJ7PdTx8Z7xWNiMLGLhhCD9IFGCMTPyB9ObCeRb/9yOXv754OgMuN0pDLnVQN2L/0cv09zMlO28lI3RgTEY8pqwuYkINEoTEvgOO18l8xW/7Gf65mI/IM/6k6vcMX8yLt+xDv6xOs1CSTkEuYl4lDiB9j44eWBPZnS2bkxUtI0ySPmE/U5kaPyQFdV6VxNXakivfd6wV3tkFXD4r0NUhW0DxrsYKmuQIyh2AnD0nD1jYpQwvJpoG4K6vcu9UJMIkEzXkR6fSGBeOPkGtd1uVuRd744jJYc5GSnmihZBCgqtGcJQ4DrjT7+3hn+7Xtc/ewR5S5w9J99SDhZcvlH1tTXgepFj7x9RHmxx0ZKvDy7wNQVYbubui/eq/dMkqxMXa8ZXpL5LUwBPIOlMZCLn3Z68R4j5jBIxI3DrVcxKM1A1Lkr3u3uy20gfnbevJRZzH9O2Mp4Czthuv0ljCQG/7vPSO6Qbg+8yhUTeGyM+VLMRr4EPIm3fyxXzI+zXhOMhMP6FnKr7yXafCoT8mTwjHNSRtxhJnSTLT0hlzWJKJUp3pD1TZRGrtodJrZui53n7Hf2HL3nWTwPWB9NwFuddyl2Osouyxo5XmsAiSWN9APJSgMbhXNSeZPS/fTe5p2ltIJ2jzLPxhjM8RGmcJkDY1ZLeOehliGFsmHd5YaT393qe9upnUVz6Vl8uMN1niIN9+1b5PmL6DUzTi34uVp7wkES6S6WS6aq8s+KOzm9LWEWs3mY221UM/uf0ibiN9uDC/WVZU4Wkhqm39OfZhmKmZ8j6baqwq7XJNbrwes7p7cV5cFjpjdsDoPInWEk3BlG8lGumKj75b8Sf/5XgP/n7PZP7or5ivXaZCSpTZf1ImZeIwf3gemiT9yD1M5NNW4CXGecBZN22ZiiJu0L8QGT2o6Fg+eXUFdZQd2aNWXnkUXJ0e9cMp4usH00wNr0+rzDGAf0ZmVMTN9zKVOVOYuy985VavD4SDOio/UUELP+qcymmqf/l0pLDbXNCHCyRkZPqAp4+wyMIRQWKSxu19N85wLT1NirHatHL6b3loJG00xaIRLy+8x8juRc56asxJQFgk7Z2kWDUM5maGaBII3vv2qGKnVTXtVJmXVobF0rJpLvw+EOnpipbQvWkaZ4Dy76OY3glb/LS9jKRzrppbKU6e1+qnW3XZuPcsX8C8C/bYz5V4HvA/8CwI/tivmK9SMDiTGmAf4joI73/3+IyP/8TnvThsxitItGgdHbKPnt9C9nMBFxT7jH/KSdT5gm0lNKt2MrlGFAQmwFx8dnKcCokmaifkg4apDCEHBIYXA3EVcY/QTiZmwkRO3XfpJ6BOWchBDFpCU7+Unb6YxLmg+qC50P6gfs8VG0nLSTXkpykdt32oYOAQJ0D5Yq6Vg5zais1UHDFERnMz7AdBwSfjAf949dmAyYpouz7SY2cdtNGRZktqnqsR5mF7ZRMDt77c7LndvJWGKaziUBbpPJMtdonHW5pk0pl76ZWCYvE8mMJVNUX0FGS8LhUwvY3H15c0c8kh/iignwT3zEYz69KyYfLyPpgH9cRDax/vqbxph/D/jnuKvetEy7l3TRyMgaZBDm9bLeIRy0/HK9G8HJObCU2pAHmp5p0rSqIj9idgbPBr6EUZ3iLq60e+IshIbygyu1htjFwbV9p5yKRGWHSbowmUKljMca7XrMOBuqhRKxlLIgPHuOPT6exKYjF0YB0FKxk5RNjCOyAW42GuzKksXjCyiKqWZNSmfbTfQxNorJRA6LcS7zatI8TcZCYsaXiHivFCWCA8nLHPznhK+IU4V2RhRL7NDZxO70x1cxSu1hUJhhJhmonQWaeUaRtUasywEFYzJ57QB0nQetUX2kM+M2DYLeIaahb/YLwGwVPUOibwFl/BK0B/2Pxdv/MvA3gD/LrDcNfMcYk3rT//FHvkjm+0zlTKI5z0fEc/swnUgJeIWpvRoV6ecAYAom+oTx9nSxp/JnHDU76D1YJYlNVp+9aq6mzsa+0+G1fsgaI9khbvYv5dR+0ehF3ycCVQwwfSwbVsusA+Ie3NfAlIBlMVlTVrooVRC0u5A7IMHqC87bsE5nZzBqWWGWjQ70xWNgTJgCWMSGKOOUcxo3SIBxxHLmFHZJrNnEAD44adL/mYSqyoNd3Nhks1pqGTIL4Pnx859vYR3A1D2Zv16wU8CZv51XCRKluRpulTG3gdjcLTIHj727ro28zOr+HK6PhZEYYxzw68DPAP+miPyqMebuetNCNsbKLmjzWjsBsLeL0jkf4LYMozMvp8LxQpkHkFReiFfPGgC8BrKMbXSdvrc0VJd23CCwjztgCiKpnJplKBKp5KYoJsm+LupqWKM/771iJc5q5iGiJlKbzYSVxIs+WUdkvY5VoyXDyVolD0A1U9tR5SBXDVJYEMHuh8whMf0AN9s8Y5So5QfZx5wxHMWZ8wBeKlEytCEHx/wgi0ydnBm5MD9e75wfI56XB+/yyZCCTCLulYfAa2Y9z/VXw+HfDT80ELzUFjb2sO07GyG4k/VTkJF8LOhZRLyI/Am0PfQPG2P+2A+5+8fqTRtjfsEY82vGmF8biM7yCT2vqqnli54UOTOZlzUQRYoTW3NGVEuYRNpJUvtyNn2aHiPRijJnLvMTd15TZ3bqmHdiGZPJk8+BSV/P5E6MqSpYNJi6wiyXmKM1pmk0SDmXhYPkZqOkNms1qCTfHe/1xL29qwZN0WWv6u7mekuoCuz1HjFGBZgKC2NQvAQIK20bm+1euSbeays5Zhy5OxSZt+l/MW6m4D6XJJiD4al0ij8f6IDcKnNelWXYSjOXjHfxMmg6DxTAyxf8PBty81LmR2Q4s99fMgmXcBg47pJHAtxV1+azXJ+oayMil8aYvwH8KT5lbzoy8n4J4Nicy5TeTuSy6SS9dfI5R1aVL0sYE5di+rDnbns52xHB3j+ftEjHkfDwDLPrVO4w6rKafasC0aT0vIjdEm2jSsRipB80AC0aCHofuk6DU6szH+ZorUxVY3R6GNQrFw74GtkQfLvTYOe9qtHXlRpUdb1S19sOCX6SCVgu9Lljml5+97F+36vQtRFVrDf7XkWd0kqBtq4z5d0erfV/HgYNltZimBTN8txOpM3nAD9ft0qYxAWKN2iAeYWfDKA+uymwp4AxL4tmmc2rWsO3uSUHGcmrLvzbAO6MpKZl3+ycO6DLfyTB68dYcpddm89s/ciMxBjzwBhzGn9eAP814O/xE+xNi59Ovnl6fJvpiIQ8+ZtPvKbW4JK6FN7riVkWuqNe32i7ExSUvIrU76MFsl7qBdR2eXguIf2SaOmJizLjhuT3EzEPikLfx2x2RAHPvQary2vMZodp+8lwKglM9wMmclzk6lrfy8UVcnmlQeTmRjMYib66+1bvO0xK8EmciKgRIoVVHRI3BRwKlzVtTVlCWelx92Gado6g6wFjNTGL0/98+7NLdPeP4oDoh5ov4tuZQtZfnV/kM3bqR7WAX7ptHgBexfkw9mXs5gB0vYWJ5Ne/NcD3aZeASPhYX6/z+jgZyZeAvxxxEgv82yLy7xhj/mPuqjedatbczpvV036mnJbSTe/zCTefKs2q5bP2pWlq3Z3aDrNeZbZm8vCVfavdmXHUYb2un5ijc5BuxkvJu1Ka2Wm7A6JZFqqOSvKStEWqEtlqphA22ykIXXhtFUfMI1xdZzzlQCw5qo+FfRsp9vECE9HyqakJRwtMOxBWDdYYwrrG3rSqehbLG+NDlnFMnjQqp6D6KyFpc0DMPOw04o9+VjZpsM4u/oMgInG+Zdaat3UdPyuPXS6n17l1Ub6q86KBwb0yeB20/2+fWrc7QgnfuM0luf168+dNnjrzoOPc3ckIfHyK/Gu7Pk7X5j9Hh39u3/6cu+pNC3rBvGr3iN2ZZHZ0wHg1dhq9n1PiZ0FE2k4H4MpC255xnkQv1OGg3pXdnixgnNuZoqVLwgZSDZ84F7fmeWQcp6E6OMBOpI+/pywqZRGJMl9H6cQo9ZfXXEx5mNXvw5CNvGT0Ghi+8z44i1utwFps4SZBaYhi07XiKWdrvc+LG5IFRYievPk9JQmEdKFGHCp75kIOGvPPJm8dxuZW8FwnJCQtk1cArjAFgJfEmW+vtPnMpBDnbdpXZhy3SpqXlOHTZx+lJNzpCWHOur3L4T2Rw8/3c7rukOf76dcc4DoIGPPBL72RzGAtJvwDmKXdoozT9Wp6DOiFEHU5TVFM8y6pmzOfy/ERqZ+Bhlm8OD6XaepIwxclmYWgLdqYYRwYW89kAgCllS9UCNqsVupX2/f6vpo62mqEqVSbd0tSi7fV4Biur5HrmymgjhFAvtroid9OvBeSCPXjC1WvL6OtaF1hV8uDY30gypzwi3jB2aaeqOXxYrwttpwlMOdgeXrudL95AEkbxDzAzD+/H7U+Cr+YsVczi/o2oJr/Zg7AVX91fSsgFW/A1lvrtaHIwwSuHjimpQ7LrfoaUK5HerAHMXYCar1XVfhUnjgLYrCLhYKekVEK6M6e1NFTezcpoSU9z6LAWBWRNoUCq7mrcbRWPKKuVPfDqbZHEl4yZZEvZGDibQxRzyRajJqmVnB22ai2qwjm5FhnZcYxD+jl1027mUxtV+mHyIcZdJo4lla53AoSOS2xBbzZI7udHntrMlM1E8oASUc5lpkqFKUBL+NQItiqOTSvGoecmbxSrX1Wvqbb8qRtvO8B0PrSCSOxLD68yA6G6zKBbaIDfORrBI9IGtBz3KbRp0znpanjT7nkpyAjeT0CybyvPytlMlDKDHSN60BRPBGz0tPFnTKTfWpVFst6ooWDaB7F6JHFMaFy2HbQQTYRpb2LaFmQyhfvkUgFN02tgWYYcstYrMGMI7Kb+adEKcccPFKmk9qX6SLoBxVovtmoO14qLapSwdX5vE0alJtlKBqMbAaIE49FqffVNDvjHFxda8ALgn/ydHqsmWV/iY8zl67UA58DWn7vZRltN14Bwsafk2PdS2SxVF7kEYfZZ/qKzePlqdxw2F2ZPeft9/KS4hm3y7AZDf42XnKL/HZ3Gcnrn218nPV6BJIZRpJYjxiLTRjHLD1O4KspD0uOvAta9KKsY3Apimw7KadHMHqGB2swhnHpMAK+tpQ3I8Npjeu8SgVsB8SA7b0qoY3awTCRtBbWSygsZttqNlIUyssAfb26Ul+ZWOYAuTVMzhxmZZJEfdaY0dD1CqDuW/WV6ZRrY1ZLvXCGHpoGUzjC1bUev+hVYyBnCXmGKK0UfLoO2alIUcqOZE51TxfdLShgnjWG7UwweZZp5M8r677YrMh+AKLPZAImjGN2oeaAeAu7eOn8+REX9ZxWf+sxL2Edr2KtzmZ1DgLdXSzhp6L9+3oEEiD38GdpY9jvc9Zx0FKcsSkPZkQg0uR9HrrLplRNoxqqVUn1wSXheIlfLAnOIBZCZWNGIWAhlLF93EecxBhM2yFlQThdaSekqKAqCXWlCmnRKzhlHvLug0hWUy6H7QZ4fqk2F8lrOA3HtVEcaLebyopixhdJQSLhI0Fg6PHXqYMzKeZLmk+KnYYDUl5KoyMBL5PpIHoTm8OM4Nb3l0SIZl21eOOEfdz6220TbgkyY4y+gp/xit9fuvBFohDRTEPk9pBfymbSew8+P+ZV52F63gzapt8P/n43GYnwimD2OVyvTyCZU6TnBKaimMbI45p7vR4QoIxJpnaq6pVallUJpoe2VXe4YcBcXbN4vlQDKR8m3kU0uDJJab0fJqJYWWKC4B5fIkdLQlOolacIw7rEdYFhXdA8axFn8IsCE4RhWYAF2wWquqR7uKBoPaGwlFcddtvBfj91hLY7REaMBG1DJ4GjCHqaRZNPdHu0VhnEVF64l60uEwtX4YRpAM9E3IBoV5q5Iq/aweEgyM9bqh95IZhyIp/NSWTzobx8sf4QXsjsos66q5nAOBHcUucrjSHo60qWYMydtvl7nmcgs2BmVyvCdovy9c1LQe3ufG3kzoLSZ7len0AikzL33DE+9MNLBKc5sj8PKDKjvotzmLrOqX3iZMh2l8f6ZbeDvtCyZBjzrI0pnAaf0WehZkX9fWayml2LGz1h1WC8xy4KQmkpbwb8sqA9r1g8bhmXBdVlz3hUYgSGs4bqqkeMwQZwFxttx65XyicolaZuQLORts2ixLkDlAKrjXyXWNbk45NOepcmjIucicxFmcUyk1YwEyaiDz4I0orDlLOSJGYZM1LXAXBq7AGDdWrXv4Jl+irK+kE2FN9TiEOHryotUlCZU+tljIFmev6kv6rnzW1ZgilYhN0uBy+3XuFvbg5f7g6ziI/Uqf0crdcjkBgOTozMF5mdMAc4SRAOivcUQFIQCWqeJLEsOPA6CZLH5/FoNyVdqNHtXva9YiJdpwFl6BWnCB4GHwlKFpzFXqlYcv1ig9Ql4WyNGTyNF0wQquctUjvK6x7bjmpgNXeY2+5guzvslHivokeREp+HDeOOKamMATXSiqVO7ljNAFw9HgraJsBVvWimlrrMR6FSMJllJnkYb5ZdpOMv4/Bqj91E5kryEEnIKH/mP4SNevt3Mxv2ewUQm9/brXW7e/NRw6CHD7IvvaeDIPIR2dqnWj8FGYmRjwKw/iDfhDFPgS3w7LN+L6/Rus+b4zFfr+Px+JqIPPg0T2CM+X+j/9vHWc9E5E99mtf7Sa3XIpAAGGN+7eMoZH9R1pvjcbjeHI/Xe71WzNY36816sz6f600gebPerDfrU6/XKZC80jHsC7zeHI/D9eZ4vMbrtcFI3qw36836/K7XKSN5s96sN+tzuj7zQGKM+VPGmN82xvxetLX4qV/GmP+TMeaJMebvzG47N8b8+8aY343fz2Z/+3Px+Py2Mea//tm865/cMsZ8xRjz140x3zLG/JYx5l+Lt39hj8nnbX2mgSSqrv2bwD8F/BzwL0VfnJ/29X9GdW/n6xdRn6A/DPyH8Xdu+QT9KeB/H4/bT9MagT8jIj8L/FeA/1H8v7/Ix+RztT7rjOQfBn5PRL4tIj3wV1FfnJ/qJSL/EXBx6+afR/2BiN//2dntf1VEOhH5DpB8gn5qloh8KCK/EX++Ab6FWph8YY/J52191oHkXeAHs99/tAfOT+868AkC5j5BX5hjZIz5Oirt+au8OSafm/VZB5KP5YHzBV9fmGNkjFkDvwz8aRG5/mF3fcVtP5XH5POyPutA8rE8cL4g63H0B+LH8Qn6vK/oK/3LwF8RkV+JN3+hj8nnaX3WgeQ/Bf6wMeYbxpgKBdD+2mf8nj6r9df4CfkEve7LqO7BXwK+JSJ/cfanL+wx+dwtiR6yn9UX8E8DvwP8PvCvf9bv5w/of/6/AR8CA7q7/qvAPbQz8bvx+/ns/v96PD6/DfxTn/X7/wkcj38ELU3+c+A349c//UU+Jp+3rzfM1jfrzXqzPvX6rEubN+vNerN+CtabQPJmvVlv1qdebwLJm/VmvVmfer0JJG/Wm/Vmfer1JpC8WW/Wm/Wp15tA8ma9WW/Wp15vAsmb9Wa9WZ96vQkkb9ab9WZ96vX/B2UTgeo9PSntAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "imshow(image, colorbar=True)" ] }, { "cell_type": "markdown", "id": "1088dbd9-e74f-4953-a9ba-97632b080603", "metadata": {}, "source": [ "The following PSF image was extracted from images taken with the same microscope using the procedure [explained before](deconvolution:extract_psf)." ] }, { "cell_type": "code", "execution_count": 12, "id": "fc5768b9-896f-49d8-951e-e946976519e6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "(61, 128, 128)\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAT4AAAD7CAYAAADzaviDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA3TklEQVR4nO2dbawlR3mgn/feO/eOBzIC44AmM1bsrEZ8WQvElmNCFBEcgkMQw/5g1+wSrF0k7yJvAtlICbP5EeWHV0gbRQEpeHcEie0Ni9chsFgoBrxOrCiSPxgSB2OM4yFOYPiwMQnYYGbux3n3R1efW6dPdXVVd3WfPvfWIx3de6q7q6o/TvVb71eJqpLJZDL7iZVFdyCTyWSGJg98mUxm35EHvkwms+/IA18mk9l35IEvk8nsO/LAl8lk9h29DXwico2IPCoiZ0TkvX21k8lkMrFIH358IrIK/B3weuAs8Dngbar6peSNZTKZTCRrPdV7JXBGVf8eQERuA04AzoFvXTb0IM9J17qYv337Zg/VTiaTgGf456dU9Ue71PGGn3uOfuefdhr3+/wXzn9GVa/p0laf9DXwHQW+Zn0/C/yUvYOIXA9cD3CQQ/zUys9DIulT1orT0u3tJPW1akck2fksHDEjfNP5lPuVpDj/LtcxtN9tadM3V5+G6Cfw/yZ/8o9dq3rqn3a4/zPHGvc7cOQrF3Vtq0/6GvjEUTZzV1X1FHAK4LBcWGxLMViI7A5EKR8oR13egdXXZuh5+vqfamBNeY36+OF2qbPNsSurxd9Js1TTqn7XMdWyIe57a5QdnSyw/TT0NfCdBS62vh8DvtF4lKyABjxwPuyHIuUDMsQgU93u2E8OrBebtja798fXl7p9YgfKrgOr64ee+sdvD3blgBczAKbGd24124aa5Sgw2QO6nb4Gvs8Bx0XkUuDrwLXAv+2prUwmMyATssTnRFW3ReQ/A58BVoE/VNWHGw+c7NRLCAsX8SP7sbI6Ly00HRdQbzJJry2x96AP/ZxVFi0Bu+p1SXUhkp7vWbC3xUr6Lehb0pu2g7KVp7r1qOqfAX/WV/2ZTGZ4FNjJU90eqHvr1ZXH6mJi9w/QuzlZhG4olL6tiCUuqTeUwL5FS8CL0PvGSvpNz2j1/g08G8o6vkwms69QYGcMKqeOLOfAZ7/xYiWK6v5Nb8tQt5Ox6CSbXCHqtvVBzL2pk3K6SI1tqOtHjM6uaVsTdttOv78i0lTWir4Orfddfg3fCAe+ILN8CteIlO4KIb5ZTaQaKEN+jEMPgFVcU7W6+zCwymBl/UDR7Hnz8264RrKxUey2uVm/fxe/PNcL1VwTta/NQO43imYdXyaT2V+owtbyj3vjG/h6ibqYVq679Y7N+NAU6RGyX4q22hB7r6r7pXaG7vDsTM6dc2+wnx07escn6YUQey51+7um5oCsrqI7KZ91YccZmLVcjG7gy2Qy40WBSZb4eqQv/dOAeq3WYURtXWgWxVD96+I+kmIG0ZdUF3NsqC7Y7KPb2/PJIzqSJb4+6epv1wHZ2EDPn29Vvx1J0NqbfuwDXZ/EZngJve99Wd1DIje6UP0dLPjZKByY88CXyWT2EQps6fKvWDHegS82VrKtv52DOWkvAq9P1SIzfowRl7Q2VCxwQ8aZ4Bjg8p5qxf0lVTzxAmY8PhRhZw8s1TPegS+TyYySieapbnr6SAQ5pGLb1VbIOTkiFKKNI4uIFunSZsocfV2QUoKZ7BoFQqUz3z01z0J0ZIWjP424dJiu8o7sFR3f+GRWnexOG5LVqe4HQ6QYcMr/Qyj3rdvf1Zad4LIOx3bd3o4zkKgiGxvTaIKZ/rpIYe1rE6ESWlaWV88hRZsu2gwS1b7Z313PQt1xNuXzEmLY8V2jXl6Cwo6uNH4aaxF5sYg8aH2eFpH3iMiFInKXiDxm/j7fOuakWbXxURF5g1V+uYg8ZLZ9QKT5ho9v4MtkMqNFgQkrjZ/GelQfVdVXquorgcuBZ4FPAO8F7lbV48Dd5jsi8jKKhMYvB64BPmhWcwS4iWL9nuPm07jI0XgGPlspHPOmCpHWXPuU7ZRv1pgpRW9v0+7o+fOzxhlHX+XAeqFwVy0k3nIqnoomiTi0rO0zUe2HKrK2NlUdzLfVYZZR7Vv5vemZDDgn2dgIn11U64uZxUSgKmzqauMnkquBr6jqP1KsxniLKb8FeIv5/wRwm6qeV9XHgTPAlSJyBDisqvdqsVburdYxtYxPx5fJZEbNJEzHd5GInLa+nzILjLm4Fvio+f9FqvpNAFX9poi80JQfBe6zjjlryrbM/9VyL+MZ+NoqrCO82Ofaq2wbasGWRmINJPa2AKX2TOxmH641PccCy9pa2D2yY2r7yPYT2PaUyGe81q0qJBFpTzOSwrgRNFF8SlWvaNpJRNaBNwMnm3at6U5duZfxDHyZTGYJkCDjRQS/CPy1qj5hvj8hIkeMtHcEeNKU163ceNb8Xy33Mh4dX1csnUatPqfO6mWOta2oXp1Q37j0Pz49lG01DtGHhViZQ+lBjzRHeU5GH7lwidymyWqeOrPOENfbQyrjhsXb2J3mAtwBXGf+vw74pFV+rYhsmNUbjwMPmGnxMyJylbHmvsM6ppa9I/H5pjUt1oKd1jGCaAvn1G4sq84l7oNX3dDkM1fXlz6zOPvO37fNmVk5cgob+GyuHDpU/PN9727B7CRyYBaRQ8Drgf9oFb8PuF1E3gl8FXgrgKo+LCK3A18CtoEbVKeLcL8LuBm4ALjTfLzsnYEvk8n0jiJsaZphQ1WfBV5QKfsOhZXXtf+NwI2O8tPAZTFt752BzzeliM3cYX/vur5qyLGuPtjdcUk/xkVjbnusdLvALDguqucqa2u7xhifdOTrz5Ap2rtMbWOfW9c5OO7P5Nln4/rhIcK4MWr2zsCXyWR6R5FkU91FMr6BbxqCY94qTStOlXSRQLpKMSligFtIjV4lf+hqa7GSj93fum0JqT3H1veq5zXC2sbWDt1+ByKNF6Ok9RmIyMUi8hci8oiIPCwi7zbltbF2QdgRFUMZFMYQiZGqfTuypM9zqrbTNWogZaSBKxqlLAu9LjHW2lTYceDVT9NxA6FKkljdRdOlh9vAr6vqS4GrgBtMPJ0z1i6TySw/hXFjtfEzdlpPdY3/TBla8oyIPEIRKnICeK3Z7RbgHuA3O/Vyt9FWh80YAQaaonn7c2B9PlWRa6rblyGhr/NtK3mkdjcp63JN6ZtUAHXGkrpt1fp8xpMmdUbb53t1tX51wh6kwWzcMIjIJcCrgPupj7WrHnM9RUYFDnIoRTcymUzPKJITkQKIyHOBPwXeo6pPB6TCAsAELJ8COCwXtnzVhRkEZmJTI9+qZW67Luno5/rjSkyZ2mBTZQiH57b1B8YYJ+lPk8RV14+646rlPsm1Lxcge71HV78TS337XuITkQMUg95HVPXjprgu1i600uLP2gGgZpAoqXuQYny8fMdRM+C1/ZH27WO3CHxTdBehA4iv/hTU1Rtq+Ajdt486qs/FgM+HApMlMF400cWqK8CHgUdU9fesTXWxdplMZukRdgI+Y6eLxPca4JeBh0TkQVP2X6mJtQvGvAGdkl7oW7JOsZvSZaTar5C6Y332mt7kMVJDjNQbE2PaRmrqo9+xtJHWY2cQfdQBcRJe4mmuwlJYbZvoYtX9K9y5sKAm1i6TySw3qrInprrji9zw0TX2McKdoA+jRm2/hqhnZZW5tV/b1Jmi74twFl9EdE5bXXNTnTH19GHcyANfJpPZTyjBqedHzd4c+GJdEhx4Jb2WusPWqe1TWDbHbBnuC9v6GWqt7dudpi0h+fsGIXkG5oUw3oGvi8tI2wchdO0KX/0eF5TWmYMXmYBhmekr1dYink3fczhgUtrCnSVLfJlMZh9RxuouO+Md+IZUQpdMAiM8fG/8DolLezGoLFLSSyGJpJZmXPeuq9Es1XE+FYrd7z6MJhHshbRU4x34MpnM6CjSUi3/VHe5hu5qnrXUmSdC86yF5nTzHe8q3txENz0henWErDCXiphcdCESTllX5d5OzyOVNFO25bh3wdcsRR6+8jxddZV9cz1fdlnX568jE5XGTwgi8jwR+ZiIfNnk9Xy1L5+niJwUkTMi8qiIvMEqv1xEHjLbPiABCQPGM/BJQFfK5KT2QzxXj+Oc+xrQXAkvPT/m3rD6bS+ROdevkDLw97vtj85Vp11XJfFs7XnU0XR/Pf0ObqvrgCPCdJnQkLpi7tlAFNlZVho/gbwf+LSqvgR4BfAINfk8Ta7Pa4GXA9cAHxSR8oG6iSLT03Hzuaap4fEMfJlMZvQUIWsrjZ8mROQw8LMU8f6o6qaqfpcin+ctZrdbgLeY/08At6nqeVV9HDgDXGkSoRxW1XtVVYFbrWNqGY+OL2YthNjIg5SxozZNhoyQ5Jc2fU5fYq5LHz5/oWunhOJT8LuScVbjt11Ggtg2Y0lhRFm4W1JwyNpFInLa+n7KpKIr+Qng28AficgrgM8D76Y+n+dR4D7r+LOmbMv8Xy33Mp6BL5PJLAWBkRtPqeoVnu1rwE8Cv6Kq94vI+/EvU+FqVD3lXvb2VLeq26vTj6RUFtt1eRTqnameVxvdT00dcmAdObC+W+arv40Os7weKYwFvvrt765tsfdlkUaF0GvlumeyEqZDD6C06jZ9AjgLnFXV+833j1EMhE+Y6SuVfJ5ngYut448B3zDlxxzlXsYz8CW6MTP4fgDeviT6Mfp+KCl/9G3bqRkQdGuzSAvmszKWdJkWhw4kPgtoS+TAeutjZysKeDmE1tGF6ipyPQ7SKYwbqvot4Gsi8mJTdDXwJerzed4BXCsiGyJyKYUR4wEzLX5GRK4y1tx3EJADNE91M5lMMInX3PgV4CMisg78PfDvKYSxuXyeqvqwiNxOMThuAzeoavnWfRdwM3ABcKf5eBnPwBdi3OhLGb1gT/hWhPYxxbkMEe87VJp9cy61iW6HiuIIrSNUEqxeN/u4hIuoK7CdKEmBqj4IuPSAznyeqnojcKOj/DRwWUzb4xn4MpnMUpATkQ5NX29jnyN0X1JOtd5YaWPAjBzJpUaXu0nX1clSpP/vy+2py7OUwnUmpQEpIjJjzCzXwJfJZBaKkhORJkfW1upXhE/SQISU1Lc05dIr9r0IT8sEqknwWdf70pPF5lSMWJqgc99CGeESo1niS4WArK7Oxkv28aNMVGfrTMpNfUk51U2RfikFHV9gMy9Du84QVcHYjFqOPjrPz8Y14FUHw5yINJpxDHyZTGYpUITtSTZupEFBd8YjysOug6vL7SFU0uucWLRt+vuQ7W2JleBC9vOsADdzrX3Tvr4NByur3aebDoNDq+e+2o+BJfms48tkMvsLzVPd9AzpotHESuTNdUhCnVPI+1wvesSrwwyV4CA8I0uoJNVF4ursipLACdiSbGW1uEYzEl8fDvqJY6Gzjs9gkgGeBr6uqm8SkQuB/wNcAvwD8K9V9Z+DKuvy405sBZ4btHz1j2XAbtMPxyBVa1mfiQbwtNM0FXVdy6H9JuuoDhRlDHnTynsh/Xdd4zZ9DNnf5TeZiL0w8KXQUr6bInNqiTODaiaTWX4UYWey0vgZO516KCLHgF8CPmQVn8CdQbU7vjRTqbNRVNvy1e8oD173IiQzR9cIFF/9lZTvM8e4stukuM6uOnz1xqbwT5HtpOyPfX2anoExSP0lfWZnQRo/Y6frVPf3gd8AfsQqq8ugOoOIXE+RJ5+DHOrYjUwmMwS6R4wbrSU+EXkT8KSqfr7N8ap6SlWvUNUrDrARlo+vmlCxzzds+cZsKT0EOzeHvJmrK8vN6aEi+9i0fyppoa6N2Bx2Vak0REJ2OTgPkZS2zfPi6lvk/ZxbKa7Hha5UpfEzdrpIfK8B3iwibwQOAodF5I8xGVSNtGdnUM1kMkvP3khS0FriU9WTqnpMVS+hWPbtz1X17dRnUG2oMMBdwKWPqpL6LRcqAfSVRh1mz9mnH3O179rm2z8lTfow+1xirl9bZ+WQ8L2u97GNxFinSw3tl+r8Epkhv5WW7HeJr4734cigGkxIUHaTT1hXt4O+6BJwHpJgQLW+jRbRH67olSRxyi4WaRgw11bW1+N8LxeRGmzBrlOqsDMZ/8DWRJKBT1XvAe4x/3+HmgyqmUxm+VkGq20T44vcCJGGWsapyoF1d8rxFIS8hbsuzBOCaaNznDB0ilNeKsy1jb5Wi4iHbuOgnjj1/DJMZZsY18CXyWRGzj43bvRCSP60Dkv51S4wE9Jm37RxgfCg5883SzCuOoY41xBSuIXY5aHHx7bXA7K+jqxHLn9Z577Sg4GjaptyfUIQkX8QkYdE5EEROW3KLhSRu0TkMfP3+db+J0XkjIg8KiJvsMovN/WcEZEPmGUmvYxr4HMw45/ku7qhPxRXRIaPalt9DQyy0r9l01HHnP9X2yD+kP3KfV3nWf3hprCO2uWxETIhfnB9rAVN4Euris+Km3BBcUhu1f05VX2lqparrTlDXkXkZRTeIy8HrgE+aPIEANxEEQxx3HyuaWp09ANfJpMZD4VVt9dY3bqQ1xPAbap6XlUfB84AVxpf4cOqeq+qKnArAWGyo9fxNSrTfVMbn3K4rVtLaJRF6BSj7MeC1lTobKwIUU9U96tzQ6o5LtiFpu7at3EBCTKypTMa1BLyPDUlSU3cz8BLeVE5fTWcUtVT1aqAz4qIAv/TbK8LeT0K3Gcde9aUbZn/q+VeRj/wZTKZcRE4lX3Kmr7W8RpV/YYZ3O4SkS979nU1qp5yL+Md+EIcdu3yEL2d/WYM18DG7R8ruYVGhbTV6bWRQFO6abStq81xPh2X9iBRt+ljiuepWodO/Pc5pX6PdJEZqvoN8/dJEfkEcCX1Ia9ngYutw48B3zDlxxzlXsav4/MouYNTP5U0/fhdCvAxpBvqkCwBncRNdVwK/iEs23VtrKzOh2OFUva/LzVCGyuwy1gWm6TA9Uz6jBuxz0ADGvBpQkSeIyI/Uv4P/ALwRepDXu8ArhWRDRG5lMKI8YCZFj8jIlcZa+47CAiTHa/El8lkxoeCpglZexHwCeN5sgb8b1X9tIh8DkfIq6o+LCK3A18CtoEbVKdi/LuAm4ELgDvNx8t4Bz6flGXehtHr8DYqgefrqEqVM212ifuNPTbEKOOapjZdR9/2SEnJe62aqOtHF2mtb4NRG8NYdUoaGp3RJZojsdtNiqmuqv498ApHeW3Iq6reCNzoKD8NXBbT/ngHvkwmM0oWrflJwXgGvpir6XIoDsleUaco9rQ9t/BOXT9cpIy59B3Xl2tOJK1dY1z3IFayHQtNxqRYKdR3HbwGDeu4HKs7x3gGvkwmM34UyANfQtq+xWN0WiHHd92v6zFt6OLqEuIukTKHoX2fyzx4q6vz0qKvrbFKezCsI3qovjq5ji9pdQthPANfHX0nDu0ybVpkUtPYfs+4S1hrxbro+jLx4Yjg0O3tuGvZ5ty7nkPq6XWXpLSxJI3ckFRW3YUy/oEvk8mMiyzxDUCXbCEhx9btUyeB2G/+LpJiV+mhy9S7dH/q07nXR1Xaib0esQaPRJlskhJ63X3P3yLQbNzIZDL7kSzxjZRQZ2Yw4TwR7iaLlh5ckmiIfqwpI0rbrDIpYmpTXVOXJOmrv6Xk6VqEqTUppN3ByRLf+Aj1rUs9xRvK0BE64MUOZCHWXUc/ZG0N3QlI/NB0fVJYkmMGVNc0P7CtTgNe5bycFu3IOpr3TxySP0A2rr7ZewNfJpPpj+zHt6Sk8kVLXWebel3SUcpknLHx0m3rcm2vZomBXbeMFNe7i8TfRbqvHBPtytOm3cUkIh01+2/gy2Qy3cgD3z6lLxcDn/7LU7+sHSgO39oMd0q26/Ztj9nfIW1ODQHbW3Ht2PiksyEz5JSExIW3qS+ERbuzwJ6Y6nbSeorI80TkYyLyZRF5RERe7VseLjmplwN0JYUMYWiDRqV/urVZr3Bvk9SyqS+ugV/EmQxz2rfYdrr0p8uxMUlBU6o4Qp+9Ln6tiRBt/oydruae9wOfVtWXUOTWeoSa5eEymcweQAUmAZ+R03rgE5HDwM8CHwZQ1U1V/S71y8OlxyUJ2WuhVt+gIQr2OgmiS/p1zxqtcmB9Oh0MaiuFhBaLqz92maeNuXV7Q9saIN39ysGDrBw8OFtYPZfye50UmKKf9nPnU0m0WnqgDyk74DNyukh8PwF8G/gjEfkbEfmQyZ0/szwc8ELXwSJyvYicFpHTW0QunpzJZBbHPh/41oCfBG5S1VcBPyBiWquqp1T1ClW94gAbHboxrbD42DqmEGknVO/XRXLyLATj1M/FvvljJaQmCW6uk7tSz1SCs/tYJ6WJ7C4UFNq/EOkn0QJIk3PnmJw7F7azL5KnaySO5/oFtdOzZDzHPh/4zgJnVfV+8/1jFAPhE2ZZOCrLw2UymWWndGBu+gQiIqtmxvgp873WOCoiJ0XkjIg8KiJvsMovF5GHzLYPmNXWvLQe+FT1W8DXROTFpuhqihWQ6paHWxwe6UBW3bq3KX1ZI33E6viaJKSYOpoO3dnZDVEr+9q2TRcefeiUyU5rqTha5+gixUyi3G11tV7KrdM1u/ozIImtuu+mMIqWOI2jIvIy4Frg5cA1wAdFpHxQbgKup1hy8rjZ7qWrVfdXgI+IyBeAVwL/DXgf8HoReQx4vfkeTtspjP2Dqf54PA9q63Ui+qDrtNq+buU1CBlImvozo8xfMZ8eXC+AJGvAVq/fyu7gMjP19uE4P9nYQDYC1TKlQaShnejnL3Qw7JNEU10ROQb8EvAhq7jOOHoCuE1Vz6vq48AZ4EozqzysqveqqgK3EmBQ7fTqU9UHgSscm5zLw2UymeUnUKK7SEROW99Pqeqpyj6/D/wG8CNW2YxxVERK4+hR4D5rv7OmbMv8Xy33Mr7IjS4GBNf/FnJgPU0qoVB8MaauqIGuUxb7+C6redVJEFaWD1tFoBMt/5mtq6mPde1D2tTsMRKkJ5pDz0d6H/gifEKjRlzx29VjA2KqkxKmw3tKVV1CEQAi8ibgSVX9vIi8NqA+V6PqKfcyvoEvk8mMl3RW29cAbxaRNwIHgcMi8scY46iR9mzj6FngYuv4Y8A3TPkxR7mXxIm6FkCIrsns05u0V9d+6cZS0TnNuYOENFEq5UN1dhHXpdaFwujzZHV17jPV9QGyIsiKzJQFt2vak7UD5mPO0+MCFI3rOtfpyvo0Zrkco0OOqRo+FmFwm+lTwKepCtWTqnpMVS+hMFr8uaq+nXrj6B3AtSKyISKXUhgxHjDT4mdE5CpjzX0HAQbVcUl89lQglNCHp09U/Vl5K1OTNgYV7zGuaWF1SmWXNVBOY6dT2HKAW3UNaivoTjmV3JnZ31Qyf4hrGhnyUoqdHvbxLKXCp1qw710K9Ufi6a70m4j0fcDtIvJO4KvAWwFU9WERuZ3Cc2QbuEG1XDyGdwE3AxcAd5qPl3ENfJlMZvwkfj+o6j3APeb/71BjHFXVG4EbHeWngcti2hzXwNfXG9f1dm2b5aKmjqTT6DoJoI7QhaWD2p43YEwlvXKaC2hZ787OVLtc7jeVAHWCTlam/xc7dfBDS7BfkvUy2j5DdfiiQkLbGmid3mXJvtLEuAa+TCYzfvZAPr79MfClcBnpUkeIu0FKt5ZQHPofWREoXVVKHV8p+R3cQNYLiUlKfeXWFuxUdHtGH6mbm0ip98NIijsOPSRES1Fl9EXQQkcW04SoXejQb2cdKfZzLTUQc3wMWeLLZDL7jTzV7ZM+31hjaj+0/pQ6HNe52ZZYI+mVOruplHfoENtHipjxZ48eAuDAM9tc8JWniuPOFQ6+YrLtTJ6e7AoHtlRY1fe1uMZtQw2nOsrQ4+0++qznrkNLqdRuK/acQ+97jFNzF7R3q+4gjHfgC/Fwd5n7Q903Gh7AaWRC+UOJ9doPJfSH0CUSo4rT1aV8mi0fwRUzGK4Xa3pMXnCYf/ylIrrox3/mqwA8/uQLuOiTRwB4/gPfKvY3g5wc3ECfdkwtfRENVez1bxO8jFoNeCWR98DZVvSUOGCUabp+qckSXyaT2XfkgW8AfM6coR75vnrrNpdv6/JviqlvB0NGtBtGiFRlGzfMVFcniqwYKWNS/i2O23zBBehLvg/AZ176KQC+8C/O8W8e/S8APO+hYoorT/9g5rj5k4lw1wmcViZnESqWyKn0zD51z+dkJ70D8x4Y+JY/ZC2TyWQiGb/EV5LChaCOkPqa8qD5tidQOE8lvb7cJ6YGh5XdUDXjiCxGAjnw3XNsf73Q8Z363o8B8MUfHGX9n82hzxZp3NUYObBcV2YyuEjFqRn8mWxCSflc+IwKXdrxHRuqQ3RJzHW66rW1WReiFOwBiW9cA19obOL0h9NC4d9leuw7LuJhnNk/ZMCsa9NH6wFyfmDSzcJAsfrU01zyqWJVsj/46lsAWH9aOfK33yv2e/qZ4ritXT8+54+uqrC3DRghpL6PLlIPeNM6Wj67NhHt9zHoZatuJpPZf2SJLzFtvNR9b+HqdKVPhbXLRSSEWINHqNQRJDm7+yorRfl0erpZTLP16e9z8O+K1/2xx9em7ej3jTHDSBelYWg3a0tDH2Oln9j7uLLabQrtazfkftj7hEi2tsEpZP9YibkDwt4wboxr4MtkMuMnD3wDEysdRb4FS5cRdBLm6Bran7ZuCi6ajovRQ9XoKLUae1tKcj/84e7uBzfm6iglPC11fDNxuZb0F+PAnII+paHY6xxDaL+7Ls4Ug2aJb3hCojm6VB+bqsjXvj39GGgaAsRNcUMD+8sp7/Y2mAiWmRdDOTBWB7zQH2TTNLHtfY6N+klRr2u/klgf07o+DhWeVkc2bmQymf1GlvhS0/eUpy9c/bY95sdwTnYffT6RrvhdK4FBKdWVLiszxE652rr0tEinP4d13MrBwkVncu5c9LHR+8Ukm6i7TzHnbCdXSMUIHueujGvgy2Qy40bJA19yEsXZDk5qHVIsLikiVA8U4oZTOjJPVigXFCpdXmbSTJUZbXYsXWmsXraLFBgSveDYZyrpxUqSvntct4ZuSKaZ1MlpExs/9sJUt1Osroj8mog8LCJfFJGPishBEblQRO4SkcfM3+en6mwmkxkBGvBpwIwVD4jI35ox5HdMee34ISInReSMiDwqIm+wyi8XkYfMtg+YZSa9tB74ROQo8KvAFap6GUUit2uB9wJ3q+px4G7zff9QrhWbWtrzrZPrWoO21AW51s6t62+T/sghOdguK7q5iW7WWMZddYeWhW6r2V935tc29tbv2G/q6tTUprVt5eDBQn+Y4jx999CHCMHrHYdWOWn+BHAeeJ2qvgJ4JXCNiFxFzfghIi+jGF9eDlwDfFBESuXlTcD1FGvtHjfbvXS9GmvABSKyBhyiWMH8BHCL2X4L8JaoGtve4NSELtxdQayVyJK2E6vUdg1o1Wtrb7fLAtrUiU7dXOz/Z7qwdmC+P7G4jomtJ8FLqM2qbJNz55icO7e7SHqnDlj3IuZaqhYvrFTT3RBpL8i1UVVVv2++HjAfpX78OAHcpqrnVfVx4AxwpYgcAQ6r6r1aLP13KwFjTuuBT1W/DvwuxaK/3wS+p6qfBV5kVjfH/H2h63gRuV5ETovI6S16ym6cyWSSIoEf4KLy920+18/VJbIqIg8CTwJ3qer91I8fR4GvWYefNWVHzf/Vci+tX0Nm7n0CuBT4LvAnIvL20ONV9RRwCuCwXKjWhrZdaiZGed0yW4g3G0ZTOvOu7i+203Td1K/6PTrSoyI52MaN6u5bDiNHAyuHirU8Js8+W3/cIoxbHdQXaiezTeE0vWjjXljzT6nqFd5qVHeAV4rI84BPiIhvUXCXiKueci9dpro/Dzyuqt9W1S3g48BPA08Y8RPz98kObWQymZFRLiru+8Sgqt8F7qHQzdWNH2eBi63DjlGo1s6a/6vlXroMfF8FrhKRQ8aKcjXwCHAHcJ3Z5zrgkx3a2MWh04jWm4Qo8NtQdTtoo5y3t7fVh9XlkKvW5dIr+owcXQwTLr2i7xSefXZX2usD1/UIKavqQdvcn04uKQ3PzpC68TRW3R81kh4icgGFIPVl6sePO4BrRWRDRC6lMGI8YKbDz4jIVWYcegcBY07rqa6q3i8iHwP+GtgG/oZi6vpc4HYReSfF4PjWtm1UGpwvarnEYC1jibRI2b6rri6xw1PrYJnmaTKrdK+22XV63difyLpCp85NL6hU9Jkyqw802GrbxBHgFmOZXQFuV9VPici9OMYPVX1YRG4HvkQx3txgpsoA7wJuBi4A7jQfL51MTar628BvV4rPU0h/mUxmL5JgjFXVLwCvcpR/h5rxQ1VvBG50lJ8GfPrBOcYVuQFxb/DUCRjbRAvUbQvZDuGrp6VcUNyFr4/2wtjV9mOnWKl9HEPr6msNja6kWiDepWKAXs5pL0RujG/gy2Qy4yYPfD0Q84bqS/qpy2Rib2+zzUGwY2z1XF3SbhtpyiUpVMpmdKm+GOAQya+pf31JKimflbHogksGdvnJEt8AzEyzhsL1Y/Yp7OuOjXn4bCV35XhZW5s/f9cPOfZhd1ktrTpko8iyrOfPz+/nosnqWN0ndTB+W7q+LGB8g2FfOFw5l5HRD3yZTGY8CFniS0tdTOiQkp6zA7tSF0T0J3DaGVJvp2vgk0QapK8ZSY8aybNNWyEsiwRVZzRxSdOLSLTbh2/fyG9JCOMZ+DKZzFIgY38ZBTCega98G7Z90/csIURLXaHGjZQSrUuicPTD6UJT7udxEWrsawpdWdu6bEL0iW2kr6ZYa5sYqbePZ9euM6XUFxiZMXbGM/BlMpmlIOv4+sD31vM5oXbJerGsonv1esS60NQtkpSSqrThs4i6tndwUp/Rn/ok4RDn8L6ek9RO+CV9urNkq+7AdH1Aymy0dl1LluJopo6mH6rNwL5etW0ETsentLnnpr5gNULsimcpCV2NL3I6HBwR1IYllRNslmvgy2Qyi0XzVLcf+oxJ1ZF4X4ZKDy6ltC+ipE0bsfRlRFqk+8qiXWdC2o3s24w6IzV54MtkMvuJ7MDcF6myVcCgeq1edCpjlKraHtskyfeRXSS0jkUbxmLPNXa5gMRSnzgWllo2xjfwpSCF4cBlXYR6H7c+lMhlm3ayzy4scvrrum6hVnpH/SsHDxaHlouBu/bre9BsijtO4TvY1G61ren3ihEvFdmPL5PJ7EeyO8t+wiOxyOpqe4nPJyl0eVv3nbjUJkSKiYl6CKh/RtJr6keXdF2hUrct1VX74pEeZWNjLiZ6uu3AevhzNeci1ON9zxJfJpPZb+wF40aXVdbGSXVFr5AVtaqEJMssHYgnO7Nv5cBU7LK2VkQWqM4roMsyu84q1VXSqpi+Ofvd2DnPfl1WGLOlolQK99D72fbcfdfR9azZ51n+X72fdtc26yW6XvTGXVHmz8v1aUBELhaRvxCRR0TkYRF5tym/UETuEpHHzN/nW8ecFJEzIvKoiLzBKr9cRB4y2z5gVlvzMp6Br1TGdv1R2A+56yYE3pjGNup+TIH16/b2bGSB78fpqq/TSmmOH261vT4syr770pbQurruZ5dN77/185GV2e9NhDznTfdpQcik+RPANvDrqvpS4CrgBhF5GfBe4G5VPQ7cbb5jtl0LvJxi/d0PmhXaAG4CrqdYcvK42e5lPANfJpMZPaUfX9cFxVX1m6r61+b/ZyjW5D4KnABuMbvdArzF/H8CuE1Vz6vq48AZ4Eqz6PhhVb1XVRW41TqmlvHo+FqsLSpra+hOJea2Ieg7KKFokzI81I8qJF42ZntM2776m1Ycq9sv9LpU22vqb6jxYdERFlXs6xcggc8kcvW5pKRwx+mLcGn9IhE5bX0/paqnXDuKyCUUS03eD7zILBKOqn5TRF5odjsK3GcddtaUbZn/q+VexjPwZTKZpSDQuPGUql7RWJfIc4E/Bd6jqk971HOuDeop97LUA5/u7My/fRrevEEZO1I5vvbxxk5RZ+jasqUUniICoiGaJkgSD3QR6YU6qTQis8p0dlKHq44h3ZJCSXTJReQAxaD3EVX9uCl+QkSOGGnvCPCkKT8LXGwdfgz4hik/5ij30qjjE5E/FJEnReSLVlm05SWTyewNUuj4jOX1w8Ajqvp71qY7gOvM/9cBn7TKrxWRDRG5lMKI8YCZFj8jIleZOt9hHVNLiHHjZuatJG0sL2GsrDa7apRUJYAY61eTy0ZTu3Vv+L4kvbZ0sQo6XBSmbjihbZX3s2oJr+ynOzvNEpF9XMh1TmkR7XK/fW4eTf1zudP46NsKrMCONn+aeQ3wy8DrRORB83kj8D7g9SLyGPB68x1VfRi4HfgS8GngBtWpl/a7gA9RGDy+AtzZ1HjjVFdV/9IoH21OAK81/98C3AP8JpblBXhcRM4AVwL3NrUzpa1IHznNWrnggqK5Z5+Nqwu6Tz9SJSKNqSN0/wbjkGsqWjtNdSVL9Q0g1WOr5bFT3UW8gKr+mE37hRrIRvQyTeHArKp/hVs/B3B1zTE3Ajc6yk8Dl8W039adZcbyAtiWl69Z+9VaWETkehE5LSKnt3CH7GQymRGSwIF50aQ2bgRbWIxp+xTAYblwfp+UbzqHpOCU9EJwSTFd+lPXRnW/tq4OsW2Fxnhax9UaJOqmdqGOxL62h3Z/sdv0SXVt7osri0tMXamyvgSyn0PWnjAWFwItL5lMZi+ggZ+R03bgi7K8tGohVHndR72upuzY2r7pK9QutK2+j+tKjFSVMjSuajjoK+zOfJ+J5W46vkpPRg4BZEcbP2OncaorIh+lMGRcJCJngd+msLTcLiLvBL4KvBUKy4uIlJaXbWYtLw0NrYRPs0oGFO+DV+xKOb3q24erzRQpNhqletwSKvOnOKa6U8NOmZU4NMplWkH9fo3PnOv5GOC6yRjvTSQhVt231WyKsrxkMpk9wJJMZZsYT+SGBqR0GFJqaFtHagmo3Ley30zcZ9t2Yvtat3ljo6iuJqFmbVtNLit1+6WKpQ6hoY9190BWV/uPErIlvcEiPJbDatvEeAa+TCazFOwFq+64Bj7Xm9x+k0VKR0Hxn3WkfKuFSjYRfQk+pwGMCV5JL6TeJgku1mWki97S11aEJKnb22FuSU39CpBsXdLl9NlvioZpQ5b4MpnMvkJZCqttE+Ma+Fxvkiadhcuxt8yG0UbSG4KGsLApKfU1TW2W0sP6+q4E11ViCaVLXSl0W6H6whT63hgn5Yb2Voxu1bXw0syzn9qtZfnHvZENfG3oorQfimp/6n6kAT+GlUOHiirqIk86GmW809ZFR0y48BnFugyKfbvVJKjXudKci5iU+CHVje331YLlH/gymcyw5IEvk5SAB2rywx/WbwyUyIKNPm3jUFPSxZDR1/S3pIsDeKj6oC6ON7SPEOYqFooCeUHxTCaznxA0T3UzC6CtMtySfvoy+kz1j+eMnjA2BDGUPuJlq3WH1JvIWV1WC12k876MUYc9WX6RbzwDX9uUPAF08udz1OWtx1aox0xhlgGXddDq/5zBpctUsFTIu6ardXWmuKa+uttGizT4bwY/T77+DUWe6mYymf1InuqmpHyj9nBRp2/ULitlVeuqIzRTxhgenliJrGHfasyurK+jm5tBx861oZFZTmLaaEOTMaLLtNf3/NnZXgLaWTl4sDgs1NWlDWN4djsynoEvk8ksAXsjSUFaz8aU9JFIsc6RsyppulZ5C139LWaVuEVTlTqq19teES3gXujm5q6Eh5H8fFJ8SL11x7quc7U++3vo8+Tbr9xmn1Psc+q6HiHXKHA2NNncYrK5Fd6fWJRUq6yRaulaEblcRB4y2z4gnlXJS8Yz8DU9DHWDUQwh4W/lmhpz2XYnoBPkwHpzG2Na/DkU3w+y6UcX8+O071kX1YbrOlfrs7+HDsCy0viCnFlisyf1TLXNYAZ4/kS18RPIzaRZuvYm4HqKjO/HHXXOMZ6BL5PJLAfVF6LrE1SN/iXwT5XiExRL1mL+vsUqv01Vz6vq4xRr6F5p1vw5rKr3qqoCt1rH1DIeHV+TEtvp1uCwq3dZccpznKwdKHbZ2qzZYWTuKa441bZp4+3jQ/zKZqabFbeUDnHKnfYvaUrXXu2fo52+/CCTuF31HmMMTILqvkhETlvfT5mVFZuYWbpWROyla++z9iuXrt0y/1fLvYxn4MtkMktAsET3lKpekbBhl95OPeVexjvwtU3bozp37NTEXyrbQ6jsVyvpVfcvJQqd9Cv9tcla0rU/da5ALt3g9P8aCWpuv0jn4JBzcUl3sfqvlPewwZ3KKen5fgchiw31sNJaz7OaJ0TkiJH2QpauPWv+r5Z7yTq+TCYTjgI7k+ZPe6KWrjXT4mdE5CpjzX2HdUwt45X4Er5VkjtzVt+09ps8hUUtVNoNCZMaW145u98lXfSyPlLeCxhOYm46piTk/FQTu1epezbRgoRL176LwkJ8AXCn+XgZ78DXhY4PaGM8rsuNIqjiwEGoNAh4gvzlwPru9Dt0AHRNg2KnmD5C4mxd9TXV3/J+JjEW9G2sanwWUsQfJw6uTXRNUi1dq6qngcti2t6bA18mk+mHcKvuqBnPwFenlPW5mNhST0noehZ2O1VDxvb2/CpVZp/g9WxdhL4pA/o/c94+qa5EVpim1bCder2NBEyX7Tanfx0uNH1JTp76k7ud9LR2bekU7zSgjcU9ymaMfYqk0bhRE1by30XkyyLyBRH5hIg8z9rmDCvJZDJ7hEQOzIskxKp7M/MhIHcBl6nqvwT+DjgJjWElflzOtQ0X0PmGTOSuoNvbTomh15XbQuI+PfuUTtZOynWJQx7Kaoyu77g6KbDL6mRVyrhc17kP9UMToQxbTI1ubTa7S40FVdjZaf6MnMaBzxVWoqqfVdVyBLiPXT8aZ1hJcG8ifY6m8ZKO8rptUX1JnCrL7pezjyFtlfs4rlWyH0/M29tlMWw4Lvr+lPGnZb1trJQhyQcA2diYptcCdgfduusRmchhrum1Nf+x1WQMffjlxbJPJL4m/gO75uOjwNesbbXhIyJyvYicFpHTW3iWNMxkMuNiDwx8ncQiEfktCp+aj5RFjt2cV8HE7Z0COCwXqimMal+3t50K56hU3rWV1xkJ2vt12f3qPGVe1MNVlThiF3wn8Nzt46p1uNqsc9cpjwuMfJhbV7hpauu7DwHPWlRi25r2ZG1tzgjXn2So+9uqKyLXAW8CrjZZEaA+rCSTyewFFLQHPefQtJrqisg1wG8Cb1ZVe4UZZ1hJ92568OUfc+lOyv2b9Cp1dQ0pzjfppaqk8NAPlRRCr0Hba2UfF6r7jG07JHddl3s9UG5G3d6e72efz2m/IWuD0Cjx1YSVnAQ2gLtMstP7VPU/NYSVtKNLyFBlf1lbg3Ipv+qUxiZmZa86UoRh+Y53/tC7P3CyuuqffsVGeoSoFmzfy7Gl9wql2u9UYXix169vVPfH8pI1YSUf9uzvDCvJZDJ7hGV7KTkYT+RGHQkvsu7sQKhRwfUGj+mPvV/1WNcbOvCtPRN/2sOaq8njWkMkkdCV6ZqIuUepEyNU6+pSt31fY6/fABKz7geJL5PJZHZZDneVJvbXwOeTwnz7+vaLbRfcb/FAHc2MRBajF+sr9dOQxN4zmx6k495oo6/rKY54DmV/u7NkMpn9h8Kuz+ASs5wDX+TbbeXQoWL3H/5w1kVi7ITqa0LdSlJIhikWBWqrh2ra31fvCJb8XDl0iMmzhfdXklyBNkNadfeAH9/yDHx1WY4DfkTlwxbclCvd1SKoSwAARZopj4HE+cOyjSvQ7LbjW+vCE1kx03adf10fxNSbYuof+pIw2M9h9IAX+rIfxLixBEJDA8sz8GUymXGwByQ+0RFM+UTk28APgKcW3RfgInI/bHI/Zlnmfvy4qv5ol0ZF5NOm7SaeUtVqOrvRMIqBD0BETidehzP3I/cj9yPjJC8vmclk9h154MtkMvuOMQ18pxbdAUPuxyy5H7PkfuwBRqPjy2QymaEYk8SXyWQyg5AHvkwms+8YxcAnIteYdXjPiMh7B2z3YhH5CxF5REQeFpF3m/ILReQuEXnM/H3+AH1ZFZG/EZFPLbAPzxORj5k1kx8RkVcvqB+/Zu7HF0XkoyJycKh+1KwjXdt2X+tI5/Ws+2XhA59Zd/cPgF8EXga8zazPOwTbwK+r6kuBq4AbTNvvBe5W1ePA3eZ737wbeMT6vog+vB/4tKq+BHiF6c+g/RCRo8CvAleo6mXAKsVazUP142bm15F2tt1pHel2/Ui/nvV+RVUX+gFeDXzG+n4SOLmgvnwSeD3wKHDElB0BHu253WMUP6jXAZ8yZUP34TDwOMbgZZUP3Y9yidILKUIqPwX8wpD9AC4Bvth0DarPKvAZ4NV99aOy7V8BHxmiH3vxs3CJj4i1ePtERC4BXgXcD7xIVb8JYP6+sOfmfx/4DcAOghy6Dz8BfBv4IzPl/pCIPGfofqjq14HfBb4KfBP4nqp+duh+VKhre5HPbqv1rDMFYxj4gtfi7a0DIs8F/hR4j6o+PXDbbwKeVNXPD9mugzXgJ4GbVPVVFLHTg+lbS4z+7ARwKfBjwHNE5O1D9yOQhTy7XdazzhSMYeBb6Fq8InKAYtD7iKp+3BQ/ISJHzPYjwJM9duE1wJtF5B+A24DXicgfD9wHKO7DWVW933z/GMVAOHQ/fh54XFW/rapbwMeBn15AP2zq2h782bXWs/53aua1i+jHsjOGge9zwHERuVRE1imUtHcM0bAUa2N+GHhEVX/P2nQHcJ35/zoK3V8vqOpJVT2mqpdQnPufq+rbh+yD6ce3gK+JyItN0dUUy4QO2g+KKe5VInLI3J+rKYwsQ/fDpq7tQdeRHtV61svOopWM5qX1Rgor1VeA3xqw3Z+hmBJ8AXjQfN4IvIDC2PCY+XvhQP15LbvGjcH7ALwSOG2ux/8Fnr+gfvwO8GXgi8D/oljDeZB+AB+l0C1uUUhS7/S1DfyWeW4fBX6x536codDllc/q/+i7H3v1k0PWMpnMvmMMU91MJpMZlDzwZTKZfUce+DKZzL4jD3yZTGbfkQe+TCaz78gDXyaT2XfkgS+Tyew7/j/Rjy2+of0tWwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "psf = imread('C:/structure/code/pycudadecon/test/test_data/psf.tif').astype(float)\n", "\n", "print(psf.shape)\n", "\n", "imshow(psf, colorbar=True)" ] }, { "cell_type": "markdown", "id": "6ee88be1-4a3c-46cc-8756-5441be595f4e", "metadata": {}, "source": [ "We can now deconvolve the image using pycudadecon's `decon()` function. It accepts file names as well as images as input." ] }, { "cell_type": "code", "execution_count": 13, "id": "c05090e0-9b25-4511-a4e5-a528705dfd95", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAANYAAAD8CAYAAAAL1Fp+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAKhElEQVR4nO3dXaycBZ3H8e/PLi++YKQKhFUirOkmiyZU04CJm4277K6VG/CCDVxsekECF5Bo4k3VC/WCxN34cqeJRrK9cGUb1EA2ZLU0GmNilJcgFstLYRErTSuKEb1Awb8X85wwlDl0ONN/Z+bM95OczMwzz5x5zpN8O89MJ88/VYWkk+s1894AaTMyLKmBYUkNDEtqYFhSA8OSGrSFlWRnkoeTHEqyu+t5pEWUjv/HSrIFeAT4F+AwcDdwbVX97KQ/mbSAul6xLgUOVdXjVfVH4FbgyqbnkhbOXzX93rcCvxi7fRi4bL2VT88ZdSavb9oUqcezPPN0VZ0z6b6usDJh2UuOOZNcD1wPcCav47Jc3rQpUo+76rafr3df16HgYeCCsdtvA54aX6GqvlxVO6pqx2mc0bQZ0nx0hXU3sC3JRUlOB64B7mh6LmnhtBwKVtXzSW4Cvg1sAW6pqgc7nktaRF3vsaiqO4E7u36/tMj85oXUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUYKYz4SZ5AngWeAF4vqp2JNkK/A9wIfAE8G9V9cxsmyktl5PxivWPVbW9qnYMt3cD+6tqG7B/uC2tlI5DwSuBPcP1PcBVDc8hLbRZwyrgO0nuHQbJAZxXVUcAhstzZ3wOaenMOm3kfVX1VJJzgX1JHpr2gcdPdJQ2k5lesarqqeHyGPAtRkO9jyY5H2C4PLbOY53oqE1rw2EleX2Ss9auA/8KHGA0uXHXsNou4PZZN1JaNrMcCp4HfCvJ2u/576r6vyR3A3uTXAc8CVw9+2ZKy2XDYVXV48AlE5b/Grh8lo2Slp3fvJAaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IanDCsJLckOZbkwNiyrUn2JXl0uDx77L6PJTmU5OEkH+jacGmRTfOK9V/AzuOWTZzamORi4BrgncNjvphky0nbWmlJnDCsqvo+8JvjFq83tfFK4Naqeq6q/h84xGi0j7RSNvoea72pjW8FfjG23uFhmbRSZp3oeLxMWFYTV3Siozaxjb5irTe18TBwwdh6bwOemvQLnOiozWyjYa03tfEO4JokZyS5CNgG/Hi2TZSWzwkPBZN8HXg/8JYkh4FPAp9hwtTGqnowyV7gZ8DzwI1V9ULTtksL64RhVdW169w1cWpjVd0M3DzLRknLzm9eSA0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ02OtHxU0l+meT+4eeKsfuc6KiVt9GJjgBfqKrtw8+d4ERHac1GJzqux4mOErO9x7opyQPDoeLacG8nOkpsPKwvAe8AtgNHgM8Ny1/VRMck9yS55088t8HNkBbThsKqqqNV9UJV/Rn4Ci8e7jnRUWKDYa2NSR18CFj7xNCJjhIbn+j4/iTbGR3mPQHcAE50lNakauJboFPqjdlal2XigEhpYd1Vt91bVTsm3ec3L6QGhiU1MCypgWFJDQxLamBYUgPDkhoYltTAsKQGhiU1MCypgWFJDQxLamBYUgPDkhoYltTAsKQGhiU1MCypgWFJDQxLamBYUgPDkhoYltTAsKQG00x0vCDJd5McTPJgkg8Py7cm2Zfk0eHy7LHHONVRK22aV6zngY9W1d8B7wVuHCY37gb2V9U2YP9w26mOEtNNdDxSVfcN158FDjIaJnclsGdYbQ9w1XDdqY5aea/qPVaSC4F3Az8CzquqIzCKDzh3WM2pjlp5U4eV5A3AN4CPVNXvXmnVCcteNtLEiY7azKYKK8lpjKL6WlV9c1h8dG0A3XB5bFg+1VRHJzpqM5vmU8EAXwUOVtXnx+66A9g1XN8F3D623KmOWmknnOgIvA/4d+CnSe4fln0c+AywN8l1wJPA1eBURwmmCKuqfsDk900AE8cwVtXNwM0zbJe01PzmhdTAsKQGhiU1MCypgWFJDQxLamBYUgPDkhoYltTAsKQGhiU1MCypgWFJDQxLamBYUgPDkhoYltTAsKQGhiU1MCypgWFJDQxLamBYUgPDkhoYltRglomOn0ryyyT3Dz9XjD3GiY5aadOcu31touN9Sc4C7k2yb7jvC1X12fGVj5vo+NfAXUn+1vO3a5XMMtFxPU501MqbZaIjwE1JHkhyy9hwbyc6auXNMtHxS8A7gO3AEeBza6tOeLgTHbVSNjzRsaqOVtULVfVn4Cu8eLjnREetvA1PdFwbkzr4EHBguO5ER628WSY6XptkO6PDvCeAG8CJjhLMNtHxzld4jBMdtdL85oXUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUwLCkBoYlNTAsqYFhSQ0MS2pgWFIDw5IaGJbUYJpzt5+Z5MdJfjJMdPz0sHxrkn1JHh0uzx57jBMdtdKmecV6DvinqrqE0cienUneC+wG9lfVNmD/cPv4iY47gS8m2dKw7dLCmmaiY1XV74ebpw0/xWhy455h+R7gquG6Ex218qadj7VlmDRyDNhXVT8CzquqIzAapwqcO6zuREetvKnCGgbMbWc0RO7SJO96hdWd6KiV96o+Fayq3wLfY/Te6eja8Lnh8tiwmhMdtfKm+VTwnCRvGq6/Fvhn4CFGkxt3DavtAm4frjvRUStvmomO5wN7hk/2XgPsrar/TfJDYG+S64AngavBiY4SQKpe9vbnlHtjttZluXzemyG9KnfVbfdW1Y5J9/nNC6mBYUkNDEtqYFhSA8OSGhiW1MCwpAaGJTVYiP8gTvIr4A/A0/PelgXyFtwf4xZxf7y9qs6ZdMdChAWQ5J71/hd7Fbk/XmrZ9oeHglIDw5IaLFJYX573BiwY98dLLdX+WJj3WNJmskivWNKmMfewkuwczj94KMnueW/PqZDkliTHkhwYW7ay52lMckGS7yY5OJy78sPD8uXdJ1U1tx9gC/AY8DfA6cBPgIvnuU2n6O/+B+A9wIGxZf8J7B6u7wb+Y7h+8bBfzgAuGvbXlnn/DSd5f5wPvGe4fhbwyPB3L+0+mfcr1qXAoap6vKr+CNzK6LyEm1pVfR/4zXGLV/Y8jVV1pKruG64/CxxkdMq8pd0n8w7LcxC+yPM0AkkuBN4NLPW5K+cd1lTnIFxxK7OPkrwB+Abwkar63SutOmHZQu2TeYc11TkIV8RM52lcdklOYxTV16rqm8Pipd0n8w7rbmBbkouSnM5omMIdc96meVnZ8zQmCfBV4GBVfX7sruXdJ/P+9AS4gtGnQI8Bn5j39pyiv/nrwBHgT4z+9b0OeDOjqS2PDpdbx9b/xLB/HgY+OO/tb9gff8/oUO4B4P7h54pl3id+80JqMO9DQWlTMiypgWFJDQxLamBYUgPDkhoYltTAsKQGfwEcmcbOuwEn8AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "deconvolved = decon(image, psf)\n", "\n", "imshow(deconvolved)" ] }, { "cell_type": "code", "execution_count": 7, "id": "cefed646-2978-4db8-ac05-92eb7bcc4b0f", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.0" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "deconvolved.min()" ] }, { "cell_type": "code", "execution_count": 8, "id": "e7e10f95-cb93-4a29-b110-ce6371780d2e", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.0" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "deconvolved.max()" ] }, { "cell_type": "code", "execution_count": 9, "id": "d7dc240c-5a45-4881-860a-4a825b87f0e0", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\u001b[1;31mSignature:\u001b[0m\n", "\u001b[0mdecon\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m\n", "\u001b[0m \u001b[0mimages\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mUnion\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mSequence\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mUnion\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\n", "\u001b[0m \u001b[0mpsf\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mUnion\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\n", "\u001b[0m \u001b[0mfpattern\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mstr\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;34m'*.tif'\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\n", "\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;33m\n", "\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m->\u001b[0m \u001b[0mUnion\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mList\u001b[0m\u001b[1;33m[\u001b[0m\u001b[0mnumpy\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", "\u001b[1;31mDocstring:\u001b[0m\n", "Deconvolve an image or images with a PSF or OTF file\n", "\n", "If `images` is a directory, use the `fpattern` argument to select files\n", "by filename pattern.\n", "\n", "Parameters\n", "----------\n", "images : str, np.ndarray, or sequence of either\n", " The array, filepath, directory, or list/tuple thereof to deconvolve\n", "psf : str or np.ndarray\n", " a filepath of a PSF or OTF file, or a 3D numpy PSF array. Function will\n", " auto-detect whether the file is a 3D PSF or a filepath representing a 2D\n", " complex OTF.\n", "fpattern : str, optional\n", " Filepattern to use when a directory is provided in the `images` argument,\n", " by default \"*.tif\"\n", "**kwargs\n", " All other kwargs must be valid for either :func:`rl_init` or :func:`rl_decon`.\n", "\n", "Returns\n", "-------\n", "np.ndarray or list of array\n", " The deconvolved image(s)\n", "\n", "Raises\n", "------\n", "ValueError\n", " If save_deskewed is True and deskew is unset or 0\n", "IOError\n", " If a directory is provided as input and ``fpattern`` yields no files\n", "NotImplementedError\n", " If ``psf`` is provided as a complex, 2D numpy array (OTFs can only be\n", " provided as filenames created with :func:`pycudadecon.make_otf`)\n", "\n", "Examples\n", "--------\n", "deconvolve a 3D TIF volume with a 3D PSF volume (e.g. a single bead stack)\n", "\n", ">>> result = decon('/path/to/image.tif', '/path/to/psf.tif')\n", "\n", "deconvolve all TIF files in a specific directory that match a certain\n", "`filename pattern `_,\n", "(in this example, all TIFs with the string '560nm' in their name)\n", "\n", ">>> result = decon(\n", "... '/directory/with/images', '/path/to/psf.tif', fpattern='*560nm*.tif'\n", "... )\n", "\n", "deconvolve a list of images, provided either as np.ndarrays, filepaths,\n", "or directories\n", "\n", ">>> imarray = tifffile.imread('some_other_image.tif')\n", ">>> inputs = ['/directory/with/images', '/path/to/image.tif', imarray]\n", ">>> result = decon(inputs, '/path/to/psf.tif', fpattern='*560nm*.tif')\n", "\u001b[1;31mFile:\u001b[0m c:\\users\\rober\\miniconda3\\envs\\bio_39\\lib\\site-packages\\pycudadecon\\deconvolution.py\n", "\u001b[1;31mType:\u001b[0m function\n" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "decon?" ] }, { "cell_type": "code", "execution_count": null, "id": "bdd0bf19-bf60-49e8-ae4d-3ba1c1d3e367", "metadata": {}, "outputs": [], "source": [] } ], "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.9.9" } }, "nbformat": 4, "nbformat_minor": 5 }