water_equivalent_diameter.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "water_equivalent_diameter.ipynb",
      "provenance": [],
      "mount_file_id": "1o5kwNA95p-kPL41d99srmmhy3PKylRvN",
      "authorship_tag": "ABX9TyMAKmbtS/XkEYnMEmY+z/j6",
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/siddharthbharthulwar/Synthetic-X-Ray/blob/master/water_equivalent_diameter.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Lc9wkCDyqIKf",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "import cv2 as cv\n",
        "\n",
        "scan = np.load('0002.npy')"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "6b4Ns5mVqzG3",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 268
        },
        "outputId": "bf14bd97-c317-4aa7-df4e-8010f400751d"
      },
      "source": [
        "plt.imshow(scan[64], cmap = 'inferno')\n",
        "plt.show()"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD7CAYAAABqkiE2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOy9a6xlyXXf96uq/Tiv++rbz+l5czgkRXIkmgytlyVLsgDZ0QuwYTgODAdQoC8JECAfYief8iEGnHxw7CCAAwFKkABO5ASxLcVR4kdkyrIs0RQtiqI4JDUcznC6e/p9+z7POXvvqsqHeuza++x7u6ene553NS76nH32s3bVqrX+679WCWstp3Iqp/LhFflu38CpnMqpvLtyqgRO5VQ+5HKqBE7lVD7kcqoETuVUPuRyqgRO5VQ+5HKqBE7lVD7k8tiUgBDip4QQ3xRCvCKE+OuP6zqnciqn8vZEPA6egBBCAd8CfhK4AnwJ+PestV9/5Bc7lVM5lbcl2WM67+eBV6y1rwIIIX4F+DlgUAkIISyIx3Qrp/LW5KT3kE4YYmD70LGnZLT3jtjb1tpz/a2PSwlcBt5Ivl8B/mS6gxDiF4FfbLfkj+lWTgVAJAPUHjMwBQJEhkD5/fTAXsb/L++zLVyswWLj9dPPJ93LqTwOqV4f2vq4lMB9xVr7S8AvAQghT3vCY5QwuON32/gP/UHrvlt0VAROThr4DPx2/PUFprOvwIA1/rqn3eDdkMelBK4CTyXfn/TbTuUxiBtkMg4mSGZcIRkcoCvbDSfO6Ok2awYUSF9Mcp3eNkCgsNSdfUR6Xts73ppTJfGY5HFFB74EfFQI8ZwQogD+EvBrj+lapyIyQLr/w4DyCkAMuVmDA7g7WLvWwKPoJn3LoX9tmSgm/zk8V/ztVB6HPBZLwFrbCCH+Y+CfAAr4H621f/Q4rvVhExd4OUZskwyeYGIP+fVBUfjZ2JpkwLXm+vCxtFbHsQMzNfmHrQt37mF3ZPg8EtHBIk9diEcljw0TsNb+OvDrj+v8H27pD7D+53SfxyQPPTN3MYGHE+kCEdYgEKeK4G3KuwYMnsr9RYgcEuCsHXjJQFqZkWXX32bgmGgl+N9Fuy1eu3OO9twrn4P1sSJD5r9/lsEo5P2sgoFt3n0Q8XregrE1j0bZfDjkVAm8x0QMjZAHAuHkMYMryP3OEVyB+12n9/044LFzzQR0FDL5resqOMWjebABPAQ4Gn8eATHMeWol3E9OlcB7SFo0P8hJaP0Q+t735bvHtfH/49yIdmB28YDVAXuiwvAYQzuoUxke3C0IeZIiau+hfZa+tUISDjWI06jCfeVUCbzL0hn4dsj0h8GZseMGrA721W5vkm0nKJVorg+Z5w+IOxzDP1i51gPfW/dcXYuhf0ywOlxkxIrahR69nCqEVTl1mt5tSQdy+PygMfgHltRcfqvHnixdUtFbldSVcPcmUPc955ACOP6YpE1Pw4yDcmoJvMPSsufCrJ/6wwxaAKszHytWQGp6r+zbkZNYfyTMvr41YU60PvrhvHAfw7N2f+C2oUorUvfB9FyL1XsWIvdEoj4QGlwFiRBlez9p+3pK84ddTlXjuyki7bB939mc8NuqrFJ9B9yHwfPL3t8J13yImXQVW+j+Fv6OPb84Dufw92tN7xqr1k74fcVaEBKR/Puwyqkl8A5Ih9a7MpMO+dnJtsjz92SeDu+/H6brm8TtNRyNmPi9+3sq6W/32+94kG/VGnlQtH/VYjh+X3oWQIKpYEgBQktQAn0CUu7Po8EOgIwfAjlVAo9JhFA9kC8N4Q0NoH48PzW/ZU8ZmO6+fp+Tsv4c0/BkZH51Rj0pEnA8+NedtR8Ug3grCiPdvxcdEccDmIGlmD5v+swpG9PaBw1Vvv/lg/+E75rIgVn/UTW37M54K/H4Bxm0TlYBtQcg6SQS3JB2lj1u0D+AMlhxWd76/sPZjwO/iyEF13eLPhzD49QSeESyGuOHriKAlsySD8e3wzErOQCpaRsumH4fMuGPI+LQ2b46+x8nyX315LiZ312zz0noAnPdpCHoKrTjQMywf/t7eL4IBq48b8qSBFfGpg82JucTeXvcB5xr8OFQde+EiAefQYYINJ0Z7BgA7sHCcavHrs6OD/naRatQhuXBLZD0fA+0733att+mwySlVZfFPct98AcPIH5Q5dQSeFviQ3gdf/t+HdlgbX2Cj368L97v2K1v237vZ+etKIBerkEXD1idybvXPe64Ybl/ZaKTztcFJo8jQJ0MbrbHnnRf6e/HZV26UOTymOzH97ecKoGHlJSf7uRBgK0UpAvf79eJTzpfuJd2AA0Nzg5SL4ZM5f5136olM3TPqQwh/90oQkDvuwPyfte8H36xGmYN11htn5MtLYv2bdfiMR8UF+GDpdLeIYn+f5+F1gHrjgPpkkIfDwSE9V9RwBXu7xp09rFDA3yIotwdoA/CURhMLDrm2fv3dtJzdPMJzDH7ngRGprjH8fyB4WsO3UPAIjIehjPxXpVTS+CBxXfCPkDnf+sCVccRc9zn1RReiRAZUhRuD1OtzE59k7wT2ko6an9uigMgpt2uzoKrIJk8xqrouyon+fQnkZOGSUTDz3TcMSdJCsB2+QonYQdD5w9A4xCg6VyE5n3vIrx/7/wdFxlN6SDHl9960DCbswpaJDqc2HWw1ZnvrYXQ+oqjP7Omg2t4ADzooHt83SiY728vRyHIw9/n0PXdpHASqPj+kPf33b8j0vXj+6E9geqEk9oOK3t/eDegaxIrNUGKAmsbtJljTIUUBVIUCG92Hp8LMGwK9xWV+0uYcYNoeDC5uwrpeNrvAyi+jrsTrnkSlyG9bn0CgDeAv3gylXsX7lpdJbgKED6IdNstbAzKs/bXfJA06PeunLoD95Xjw0PRZLftvveTYVDORLPTorGxQ0uHjoncV8sZlpNchbcmb31OOB5k5Bi/eTgCASe7CCcP4OBa0VE8xwGAJ1/3waIFXTGeIfr+LHV2qgTuI8MkIG9Ke99wNUTXP8cq2ORmukSSPH5jl97CGCGERGDcsImo9IPTch+8gz+cAmifZygufxJglwB9Pq9i+N5O5jX0w6KWumPNDN3bA1kEvRLunWP9b12lkrsMyPdh/sGpEjhG+rO/FGOfhipjOqo2i9gZhgf68TNyNDN9OrGUU6xtOp3Y2ArhrydFBsKBht1zy955T47fP8i99s9xXJz9uMElRIaUBZkcx7ayfkBJkWNsHZmEImlP98xN3DeeN5KUpFeKqnNvFo02lTuvbTCm6t1rV5H026D/vAEsHZL0t276dnANXF1Il3vw/pBTJTAgq7O/jApAinSWkVgBWD04uPoyvE/AHDI3MGzSWW0NInex6dCJj5mhug8g4SE7YXqPxz3PcEpugnUI107uzytMaoeB+OhKqkCU6OMQBouJ52yVgEL6a1laxWCRWOHaxABCmFbJWCIH4X5tdz/F3Vo2J4c13b2/f1yDUyWQyNCKPUqOAeLsYljEgSLlyIeI3C82MPLAd8pVHz3W+feiTTIY5BSl1qj1fjSPra3RtkaKEiEylBxFEHEohGfhWAVwvxm+vcdjzOdEMTrwMidTI1RUVKBtw7K5R6MPqZvd3g0cY9Y/0IpG3f3C9ax3yaTIfPuU5GraHmKDkm0tCGMbtFm6UGwM8aWSph677/6pu7ez4pbJ9rt349w9vLetglMl0Bff0drwDyvmKdAD8PKkzDU9vkCX/tqREzr+ajzbeAQ8i9eUvsaAsZV3S1blrSLiQ7O8i1SstkU0w0XjZ2Tjfzd+Nu8pk+QG+zyEaGL7wVxka4zUBiOxzoxNSjuitCUCQebbsVuZUKBQLujq79VYz8hIVi1ZWs2Cirv5Hfa4yVLvsWz24/u0EXfxz2lTDsCqJdA+w/F5E4L3dgLSqRLw0sUADEKUSFlEv7/dsTW1jV0iRUmu1mj0IcYuuyj1IKBIsn31dxuLjaY8f4m1tZ+dDFJkZHKKkm5wVs0+loZj+1nvPob83UFl4WfYIltDefO+MUuMran1IcZUaHs0aDk8UE5BXIOw/U3JMbmaciF7kSfMJZ4sRjw7g81cc3a0JJeGXGi0Da6CO1hbAKeUVdzm7AUlDAaBtoK9quBuNeO7hxu8dvQk1+VdbpSvUZs5tTmKSqAxc4cvJDUhrF0mAG7PQrJB2a9yGqwA8R4uZfbQSkAI8RTwvwAXcN3vl6y1f0cIcQb4+8CzwGvAX7TW7rz9W318MhQCbGcCE8E6i14xtS0mAbp6ygJIzcS00wRzeq24xERucU5fZEtMWVMZk0yghCAfMBSM70fWwtJYKmO5KysOqdiRdznkHkt7wLzZSWbmcIybnR0jMSQyrVKYw2xcZptIIanNHG2qeLy1xpvR3mqyw2DboNshJEqMULLkTPEcZ+xFnpfnWM8VazlI4caYElAqy1hZJspQKEOlFY2RLITFWoGxAukHvLHtO1QymOECISzSXz884SwzPD2F7TKnMhepzUX/LuHOwnKzqvhm/k0O9E2Wet+3SWsJBQAyumQJJXsIYBUopwjeo4ChsPbhtJMQ4hJwyVr7b4UQa8CXgZ8H/gPgrrX2bwoh/jqwZa39ayefS1qGFs58B+S4xT6EKJGiQJtDjs/skxEbMGbhZgp/fLqP8+clyqPlUmTkckIux1ywz7JlN3h6NGK7FGzmhrOjmlIaCmk6joTBdWxtBbWR7DeKhZbcXEj2a7i1rNmxc/bkHnvcRuPwBHDKStsaazVLve98YltFlyK4NkD0raf5OQSS/epNjF1ikn06VZMGy6SH9g2U6Dw+e6FmlHLGM+ZjPJHN+N4tuDCq2CqWaCsxwLzJMP5Zhwxta49P7XUDvysyMTeCwgjbSqUd9CosbxyNuXqk+L39XW6p6xyYO2hbY5J2dEpxSaP3B0KG7Rvrh0DdJPJuWgTVl621n+tvfWglsHIiIX4V+O/935+21r7pFcUXrLUfO/nYd14JRJZXz9QX5G5Q22qAoGMGOrxcxQQS5ZCrKc/mn+WCOcNLGyM2csNW2fi9nMnq7ifMbqCk9WbscPIstIMg/W6sm0kBMtk+V2MkCy24vsi4V8G35vsciEPuiusRma/skbsPJJoaYxvmzQ7aVGgzj+fq+74BC8Aax3D0++RqjVKt8Yx8ibN2g09vlGzkhu2yWXkGJSxK2M5ABaISEMn2/uAXx8Xykn1TpTCEzoTrCtz7qK3AWkFj4bBRvLwr+M3m99irr9HoQ8rsDBZNrQ+9y9AkHI62qMngPYWsyXfFIhhWAo8EExBCPAt8BvgicMFa+6b/6TrOXRg65heBX3wU13846ZFZjuWA9zPPbGI9yGgOCo+Ql/kWhZyxqS4zthOmZsplucbWSHJh1DDLNRt5q1yWRsbZvfH/D1/ZXz/5XQjbEpSlQeA6tBKWXJqoWGojGStJYwVjJanNjIWesq/PoKxACsHSaCyWBsuCikN5yN38TZbmgNrMI6oObbxeiRwlSgo5QZEzZRNFhkR6IC/ncjZlPW+ffbOooilfGYnxg86QzNAnPH94pgeVVElIYTtuQ2znxDKQwlIKG5XoSBkujUte3P0UO9lTHBS77Jtb1OYII0uMcW/AUJGWkH8/ydu2BIQQM+A3gb9hrf0HQoh71trN5Pcda+3Wyed4Zy2BtvY/DJv3hSet1HSq+/oyU6EgiJJj7yMuydQahVrnM/JHeXo04vPbNdOsYZw1HDYZjZHUxnWQdNYL2wbv089gUlgaPzv1f4fU912d6YI5ba0gkwYJUUHkwlIojRIWa6EyintVwc1Fzu2l4MbCcKAb7nDAUiyYiyOklUgkEztjakdsyRGXxhnbpeWZ6ZKx0kzzGm0E2kqOvB+fDr7w7E2yLbUKgvSVQOoShWdPW+8k9kQmbJzpj4nVHIvvK2FR0tAYyX6d8T/ceo09c90pAm/mL+udjksVpZdqHjEDv+876xo8BktAOFTp/wT+nrX2H/jNN4QQlxJ34ObbucajlrSijzPZkjCVCEy9pmf2O78+osJ+2yQ/y3p2ief1RziTFazlkq0CZlkywJsMnQz+1ftZ3Wb7M2LPZegPANM7h7Wis01CB4I3VoAVWGGjuR2sESEs67kmk5K1XFCbgsqcQVuHwDtrwwF3SkApLbPcMFKuvSojMXUen6GvAI6bK4faYYh43H+WB5UAIkphIVgevXsZul7YJ7zDkTL81OQF9usX2KkaXuEa1/W3MNka2izR+qh7grQfpSSuoBxsSFV+9+TtRAcE8MvAy9bav5X89GvAXwX+pv//V9/WHT5yGa7nJ2hZby2BxMT9++BWJsuoAH54u+TCuGKzqFj4mW9pJJWR1Lbl0inh9H6TbLsf+/8k9v1J+6b+cHodl4PgfquM7FgTYf9ppplmmnNle99BlLcmtDfpwy/GuzONFSy0iufqi/UDMPX/U0tm6Dk7ttoAwPegYqyIVtiQZWKS95Jes0mUci4NP3Bun3mTcW1esrx9kbvqmttPZC6kHJ/LU4lDBGGQP/Luuw9vJzrww8BvAX9I22b/BQ4X+N+Bp4HXcSHCuyef6/G7A8EFUHKKkg6td6Gu9nUHBNj91rQuASBFySg/y1b+FD+af5LLE8tabrxJ7f5XfmZ9K5Lu3x+wqfRBrrTjGlhxFSKtZcBsTvdJu2X/HOHYk8C0eK4eUBnucWigDoX1+tctlI4KI5wvgqe4mVn72d3YVRO/337BHQj7p3cfju+7FtYrjf5z5dJELKPSkoWR/NZNw3V2ed3+IbU5otELGnO4QlPuRxDcxneqMMkjdgestf+KlhrTl5942PM+DhHC59OLNllFeq6+tlWbwBJj8J5qigSRI0XB5uhZzvAE5/U5ntq0XB5XjLMGm8x8dTq7JNcfNGnpZtWvhLQGjkv3C1dS3rxF2O65jrmXIQkD436gWxc9SbafcK2+sggijtknDFSJc3Nkz6IR4Xn9fo5TMGBF9BVIcp1UOVmG2zkV0/scFHcuLNOyQgJPTaaUiy1q+zHuyjfZtVdQjNroAXT4BCug9COK0j2MfOAZg0LkjPLz8Xut92mMG/jGNDR63xf2yF2sP2SDkSNlQaE2OF+8yH949gXW84aNvOawWVIbyWGTteZtvJ49cXYF17ltb+AoP4iDPy9xsxf4UJnfN5izAeSTwnYAxyA64ROEwRJCX0PWQZbMeMcBcuFzOF+/2/YHc9gmhEX7+0jvMR3gSth4D0AHCDWANYEK7AFE//yhPZRvh/Dc/TKwAQzVPYBVcIzi8u9hSEGkSm/pJ4DvP3uEtoI7y/P89q3zfEFcZ1Jsd3kWOHr1KsVbeuvz3UlD/kArASlGCJHRJHFucOZ+o+fR97fWxa6DSSZQ5NkG4+wMn+H7eSIvmWYVuUfp3TnagdwHl6AdNH2fm2S79T5oxw3wqH0Y6AGwy+kChO753AALHTw9l+4NeHdTEvy5w/61kSfi03Fg21XSTmqNrNY29Jf0IGT8TKusCmniAFbCECj+1kKW2ApB0fbDo+nzhnsNs3SY7cP9hzbKPG5RadU5LjaRxzruZx2k14AWI5lmhhfXC8T+T/FVvsVuc9W7mZJMTtFm7rMMh2jagen4znIIPtBKIFNrnlvfQhJKTr0S2I/bnPkfYuCO5DPLL3BWPM1PX5RcGO+xWxVoKzhsXPqKBDRtZ7wvgDcwq/SjBcEMzqXp+MUpbrAcCCkWUiMEZKJ7BYMgS3gIwdoIs3kmDTQZVXLO6DN7F8O1SbfDpwBbfDxhYySg7/MHyyFsC21WKE0uDSPVhtYMIlrGomdEHWcxGyuojUInLoNN/P5ARgqxf2FkRwnEdI94D+k1VwHWlWcCjrRCAiOl+dTmId+3Zfnua9vc1N/A2gYpC3I1jZTjjkRwuuyGEd8heWSMwbd1E48BGBQIhE8DdnTelr4Lkky5GL82S1Jjt8jWOJe/wOfki3x6y7BdNOTSRE2fSp/N5s6wKkMgXxgIIQad95REJmwc3NDOVtqETpnu6zv3gP9dGx8C7Q0sa92Aa4yMJCXbe8Z0MFvaMNtx/nxkP+IUWS6Nd1vC4Ozem/QWgGNIun2U7LaW8PSsqJgSMZGvKGLbOKsh7N/uF54vRGaUNGjjKMqpQgjP2b1Ocs+0733I7UuV441Fwe2l5B8efIW95jqL+jZSupyRkC7eF2eZPi634DEyBt9rEiIBsTEDo88ThKTIyeU4csIDEJirKbPsPM/oZ3h2HZ6bHXFY54OEnj719ST66sr9JQi/8rO+kiYaiEFJCBGsAxvPr1S7T+d8PSReSUcAitaB6N6jthJj3IwuhAXvFqxwEJLPwSJIXaH0nFlybKE0uTDRDB+SCNTRzuAC2zHlhQBtHDdhxZIKioNu2+ikMEsgLWmtOhaOEhakAY9VBBmMZvTaAbpcjODGBbclXOfZ6ZxzZcaX9p7naiZZNvdimDkUpLG2XnEN3umCJB8wS8B1EykmSDlKkn9AyTXG+RlyOQHgqLlDo+dosx/zBZ4dfT+fEE/x55+euzi/ViuD/f7P0mXxpccaoJSGUjlzL8xIYZ/MK4QwcBsrV5RAe6521s5laz662TctYiLiwAqDw3HyJY2W8VxKmDiLhustdBbTccP9pwCoewaLkm5GH2UNwQivTEalFdqIaHGEXICgIDbKRc+lWLVY3DN321hbuYIhaCvjc4ZnDPtqK53l45+70oojrdBe8U2zJirDeZPFaE945lQkRO5HOuv3CVGpksyF5bXDkn9zd8lXxe9xb/EqSo7c/QVeQd8N6BQ7eVShww+FJeAby4cCAwNQ+qw4IFI9jWkI3brMt5nm53hJPc2TE7kSGz5OAfRDav0Zrz9bjjyIl4l2UOd+wGvkip/tOrdAI6JVkHa2MDjC4A6mcRj4AGlEPHx2SUYGlDuH8ECZkpZZvoz3MMpCfr6J10+VUTC5tXUDPPjc2ghGqmac1dyeT1hoFZVJiAQo2Z19pbCDPr+1JyuE46Rtx5DY49oQYyiU21pBx80L4caR0m3EBQYJXuE5gvswhCMEKZTmTKH56GzE4cFLXBnPuFt9B22WSJG5EGK/JFxM0X78k/QHSgkEN0BGs9/5X6NsA3DZbfP6tivWGfcpuZR/kufsk/z05Tkj1bBTlTEKcJICCP/3ufv9sGEA+2Z5TYrFC9HO4h3SEiHBRsUOlvuQ2NA9aB1ChO05cqk7CiAF7MCiBBRC0xiJNpIb8wnTrOHZrdseVBOMywWZ0khpyLKGvKgRYaA3GbpRLJYl+0dTjqoSibM25ibnzPiQrbU9blx5hoPaKeOR0syyhpFyOAu0CsR5/kmHTz46y4Q4m8fnoZ2Fwxm0bRWgs1As0up4XIOiEJpxVjNvchZacVjnWFoFMM4al8rs273SiiaxDgLQOFKauVYdVmif92GtYKkVm0XNnzy75Fy5xtWj7+P/aK767MxQM7GNGMQoVVwE9fEqgg+MEpBiFD8b2yCsifn7aWGIWD7KNhT5Flv5M/xA8Twvrmu0bThsjscAwgxwHA/AQJzlA/qeiZbR5mZb0/GVpXAx9MZKNyATN0IIy15VcGeZcWMhqQyMlOPw1wYKCaWC9dwwzQznR0tKa8isgcwphaBkUgsinVmVMBghKJUmkwZjJVWTUemMeVVEJaekQXmmXPiujaTSGcsmozYqDkqJe6aqKvw1EqIPLRYSIEYpbGKlrPrnLszotmnfTn3fPaYDC/dMrl3b38O1hLCMsoatyQF3DteQyzKGShdasdCK2tOpM2HJhUEqd1zlLZr0GZSR3TBvYv3Z9P4BoxUbRUMmFT+j/xyvVQd8ufq/kbJAoDoMw7YWQfbYi5F8cJSA9Ov4+eQfi4lMQBcBcBIweCkLJtk2T5rn+ei65uMb+9yrCtfBwjlJTPzw8pNrDoFjYZYI6byFNy0DcUfhO6r3/0On0lZEHCI+k7DcrTKuziXfOpxzIOas2wk1DYdiwZqdMBIZF8uCM6VkpEwsRgLOGkjxgvS8QQIeUHo0XxtJoxVVo9A274QohYBGSwyCkWqcr++/xwQdLKVnUjY+bBYiCmkOhfBWQ18BpDN+8OcD6g8OJ6m89dInMAW+wShrOsqgL7lqWF87YF4Xsb2FaBXAEocTCO+2CNsSjXTPIhCJi5T2hRBCju3sJ4JZppllmonKmO6v8Qd6HHNWhjgE7RqRp0rgBHFztBTO5xciIxSMDOZWrqY0Zh4zvPJsg58ofobnZhnft3UELLlXFREIVAnSK3vhIG2Fm8X999D5CtlFwkMe/2GTMfLputOsIZeaUdYwySsypdk5mlJ74OygzrlXK75ws+aWuMeu9FVt0CzYwxpNJkdum6l90pMiW5ZkVcl0f4uSEWMz5oc2trg0bnhxfY9R1lCo1iTuYxlSWMZZHZVApjQTaah1FkOSYVCqzMSBVhvFkc3JRItHKOEU36ioGJVLRqqJNRNy2bWApPDlx5JJvTGSRZNx9WjKjXnB1SPBjWXDvqnYEwdUoqKmohZOsY/tNA6iC+Ys5/KSz25rzhQ12+UiKtv2mpZcaYpyybRc0GgVLZtSao50xrzJWGrlgU0ZrbeJ0ow8d0MEhefdtEViJYCLFgyFj4ObOc40L67Bf1T8LP/fvR3+cPHrHscquhWcwGMDoWjNo88v+AAogXaZMGfuDy2s6RKFLJpJcZkz+TM8N824PGnYLJbs1gWNVwApqu/OuTqTCDyq78NLme/coaBlmvgSOmAYOIVs/ExlabTizaMp+3XGvVqxV7syYVfkG+zYaxxWtxC4Gv6hjmFjl/5ZWlt3iQM/l/KATJbkcsRrB+ue1rzOdllzfnzkow/DHSn35j5AprSP3/tYunEzcjoZGVrrwCkQ7VwCYSkzF5pdVkXcL4TOtHEKDwLlWbFXFSz8oFv4hJzrc8XtheXNasF1eZMjtc+R2UGbOhY7BSjUrG0HtWBen2djf52tQnFnmVNIyyTTnCmXTLKa0re9bhTG30ehGrSV1AJyY9DSRAJVbYV7zzL0NqDXP2KIltY67PejVAKTc5xpLgrLU3KDa+OPsVtd6Vit3SXVBMN5lm9f3vchwpAclKk1LLqt+4ZjDAIdX+vHxn+Zj60VfH77ACEc37zy4E6KYAcZykpzM+xX8E4AACAASURBVLqJfPhAiBHC1cYL/HVXVKSO1sLm+IgiaxgVFbf2Nrh5NOXvvzbldbPDK/qL8b5VTG5qSSNhZZ4QY5YiJxb9xKBEyTTbxvhttTnCYpAi5/Pih/mJC4JLkznTrKZQOprhaYhLSTeAt2b7zKaHVJVTjvPFiFpnLOucRZPRGOWxC2cBXN68y8b6Hvd2HQA7Gc+5s7fB7cMZ95YjausGfuZJUYH3n0nDzcWIf32z4NVqn6vyNY7MDsbW6GQ2DM+ZFvsMC5ukbZQuYhLa6Ix8iifMJf7shZzn1w65vH6PQjXkWcPu0ZR5nXNmeog2kjuHs6i0Fp5JeVDnkR8RLLoUH9JWoH3KePj/ftyRFEiWwnJ7mXNnqfiV/S9xd/Ft2hLy3aKllvptYgMfwBChQMTswNhRfOFHgYwluS2aUbbJTJ3n+UnB5XETwb86eXErNF7/f1AAmTf7M/8CVULrrYyKnSOXhlFQDsKy4VF2gN35hKu7m3zx9havHcA37RvsilveBHQJJrlnOtZmTrv0Vgh7hlh/3X63IIWkEBPnKghNjQuFSgHfFTf4lzcvsp7PKKVgo3BhQhfjh1zApzb3GOcNSz/ID+YTlnVObRS7yxHX52NePyjYrQWVgVD1WApYv3aWaWaZ68Cag8MGjhqojEXbEK4UyMT0HynBfm34tr7FjrrFkd6JFo/2RVtjDUMgF+PWGkJjLBRyhiv+eUTjZ1GnIFwbHXCHaxJ+69bTfH13g2em6zwzm/PZS1c82GmZVwXWCgrVkKlgLY2RITvUTxRLIyPWEyRYfbl3M5RXBmkFoyGrIPYtK5gogxzBJ3df4uroAq8tv+Rc2QBHHbN0/KOS960SiKsFBfZVsnBEWI0mk6XzFxWcyZ5xTMBpw6Xxkrluw2/xRfpzp/x6k5h3ShpGKvix1lsDglpn0Sdcy2sK7/eHmXI2mlPkDTsHM+7Ox3z7YMY/2bnDHze/Ezt4ulJOJktirQNBogjaEJKxDUoUhGWvpMjJ8HUSMMxjO0lumG/zpv0GwiiUyLmwfD6WCSttyYySJydlpNfu1wVKWCrf+W8tC17ezfjt5avcNW9QmYP4HtIZOszEQXnJyNNoC7bohAQzUhsYq1nqvfbZRenfZ4jvm3gNZ/0YhKijBVTKc2hqz//w1pIK2R2w0Lss2OU2r5ItRzy5/DjfvzzL911oXZJ5nSOFy9Uo85pcNVQ6Q2LJpeawLiJGIBMl0HIlBLnvFwqofSQhpTof5x4YYJJpJmi+Z33E+tGTXFF/iFvYxRBXojpVAsdJuzpOSAhy69WVZGoUO+Mk2+YlXuBHLjnQKiiAIKFgZ2ryO1Zblxcf/P5AdNmriqhEpnlN7kExJS2FbNicHLExPSDPG3Z21/nfvv0krx/WvCJe5x7XIn3Z2Jql3o+DqS/G+hRUf8tK5BTZjMoc0Jg5Qkhqc8Rd+xrr2SXW7Bb74ibWLr3FoJAip5QzFDm73GYi1tkyZ/hIucYTY5hkh+TSUIbYO65Ows1FyT+8scdddYd9exOLIRNlHMyh4Kj0hVadQqp7999NUA6DdckBSmRMsm03k9tlHOjuvbgqTkEJZLLE2IaaNspzpO/Ec4ZFW937au8hrIGobc0V8Q3+6fwi3/7S0zw1LnhibPm5j3ybIm+4unMGbSW5Ujxx5jYAr926wEg1kU1YG8lhncd8jxAqTDGgmm4V6ba3tq2Q9rUAFj4/q1nLM7618xluqte4t3iVsJJTiBI4y+/Rli1/XyuB1YUzDaDijLI0FVJmbIsnOTdSnB/7XIBjCECpBLM+freBZ++vlFgRwseTRx4DCDTaQjVIabh9b5Mru1t882DBd9UVblQvx5ldCtW1BgQRHVaywBUybVYUhBQqHhMGSm3nGBsW/JCtlUBOJkuk/27Q1GKJRqOEoFAJ9VUE7jvsVc5XvSJfodJHxEVCBcRCmYm1IoWKbgJw4uwVFvGQXok0QoENM3/fF/bUb/JoGQTRsYSXz8RL2inuJwAyrNUszBG1OWJf3uTg6NMcNjOu720yK5Yc1gWCJZNCUxbuvKMQ7rQy9ocQopRJKvRbqSh1HFYwVpr1XPCUPcdSLNgVr3UYkq0lOFDQ9G3I+1oJhFJNxidhWL9NyYJCzajNnM3sKX5u8yk28yYy10ImGByvqcNvMfc9xn+tiyYYycgTbEqpUT5j7qmzt1BKczQfc7AY8+beJv/NyzNeFr/P3OxgPYMlYADzxvnAsbSZVRxUN8jUiM3sKQ7NHeb13TjDATRmGf3fgHuE2W5h96iEC4VmsqQ2R8zUeS7Yp/kuX2ehd9nILlPZOVfkK8jFi+zVM86PCqZZE4k4tRX8mzsZV6oDFnYXIRSldEi8Qcclu1IAU4qJ++4Xb42UWNrBGUz9FN9o7NIj/g2VPui4AM4fNg738INA+6IwQsj4WYrVrtyuIgW1nbe8EVsz1zt8XfwOLy8lv/HNCZd5kZ/bXuf7z9/i6ctXuXdvA60znjhzm7v76+zvblEqTak0e3VOY1wl5QktN4SkX6UFTYP0sxHT/QEqK5hmmj//tOY3bzzHm/UfUbM/sNT6o5X3pRJwGYF5TMKAlixEmAGtYSt/hifMM2wVOlbDPS5LLP0cF8PAdrjr2ggqq8iEJVM6+v5rpR+QwjAaLTiaj/nilWe4tcy5vVS8ob7JvN5pzXqIndYkWIbwVkzAA/oVCMK2LvbhEXM01lNQFTlGNBgLuZxg0OyJexij/VoBznQvGCMQ1NZXG/Zpxdb/1cZSJ4OwMm0lXSVyDGFJLoPBUIoZWtQsuOfvaShVNg17+UHdm82lyMizscvxSMz6sKKSW8kn9IU2YhDPe8JCr4VXZNouMX6JOa0bbmVX+MrOFofNeV7b32CaNUyzmqc277KoC0+ldpZg6S2CxuNASz8hSE9dds/y9mS7FHxy8SO8LH6HeXW9fTYUVhiEfXSLnL7vlIBAoOSYTE0p5JTKHGJMg1IOUApgoLY1z9kXeLIYc6ZwmVpDdOCQWhpAQCUsE59VBj4ebolZeEdasebDbEpaZsWS8xs7fhEKGE/n3Lq7za+8Aa+rV7hVv+KvYzomayZKhFAs9V4c9EIkSgBDzSJGAcLg0WZJWPlXKrdQaLAKjGdJFmLi/WvIxQxNzT17LYYSM0pKO2JqphRkGGtprE+USSsQ+3OCm61FAurlckIoemZ8SvaUTRrRsMuV9pl7AzINcyqvCBuzjPgCOAUzU+c5Mjssmnl48d3Ub+FusFWCrUsVKkfH6tG0bT+RW1gMS3vgIwo1SuYc6Jv8hv1nfGFHIu/l/Fj2ozw/E5HxaQB8pmKpNC7L0kUMaq8EhL+nACqnpcz6LsCqemz747zJuDhq+DG1zrV7T3SUgGuKHARYuzjmLG9N3ndKACAUA6mRflA0GF1H8/Cl/M/wifEGF8eWaWYGB3+n1p2ntRZ+TbpAABGizfILWXLr0lAqx8A7N9uLob97h1N25hP+7h98nCtHDX8sv8bSHAwiuwLFWfU8YzvhVb6EEJJcTjuzvLWGeXMPi0aKjO3ieRQZN6pvdJDzOimd1pg5Sy0xsqEJKaoSjNUxjCZQVPKIqZ2xLWY8NcnZLi25rDuLgShh+d4t+Ei9ycW9H+eKvcdr5is01oVig+LRdkkuxyhR0oiGhmU0zwWKXI5deXbOM2eP3eaqcxcEXFQvshRzbtWvdP18ag70zdZC8AO9MgcYa3yFnqZzzOrCq10FFKywiiP/XXtLwGBie2q03/779o/59u4Wf7i7xROjnBfWNRdHS2Z5TeMZkCOlYzhwob2F6Ilj96uCHO9MtLZRt+y7YZwZnrn7PIvxHruL1wjLvYV7fVTy/lQCuFVx3Uq5gUQDgSJ8OVvjxXXNet6QCdupGw9d7nxqtoWqte4vFLAMSkCReRZdoTSF0oxHC4yR7B9NuXm4xtWjCf/P/GXuNq/Hl3TcmnRjO2FiJ0nYqw2Ltb6rSzXN5Zg1u0VmM25GBL71t2O7JIMzdGwRFtK0OqL4E9aZ2SmbRcb5keX8qI6VjdKKRxdHFVUhWZoCDja5m11ibnapzVGMDoSZV3nyknvmMKO7MG3BhJldx0jj3AjhBuTMrKFk1hnA4TlqO++a+Egas+xwI/qiZEEmymgdpm0ZzhvaLKDu7nM3k9NiuF1/h7vidV5H8sL8s8zyLaZZFku+gXMLFp4n0hiJ9UogJB8pabC9kmudKsgD24IUnsF5Pptw1zzLnriCtRVhJSzXrx5NlOB9xRjM1DogMWYR+QGpXJp8hr+w9knOjxo286YL0tAtjxVSR1OLILUEHM9dx/JcQRlYK3hi4x6jvOJgMeYrty7wP109ZE/cY84ee8312KnCbKV8lpi2yxjzH6kNMlFSihlzu8ve8uqK6axkwbn8BZ7UT3FX3uNI7HN9+fXQZnG/1HrI1IhMlDSRbKNinP4l8UNcLka8tKUZKUd86odCoa1UHHLtXfKToDKC1w4zbi4M/0r/rg9RtgSd59WfoBEN36m/5K8tKdQMgVvePJMlpZxReTbjRflRFuKQa4s/8ASfdkFTY5oOeBjCj8EKTNeOcJbUmJfED/Ox6YRfm/8ue8urSJk590eW0e3aLJ5GU3PU3InYSGUO275Iu86icysbT8Ue82nzGZ4el/ypC4dMVIOSNtKdKyMRwHpRRfZoKOEWipRknqGa0otDG6dk7pSnsjSSO8uMv33zn1LpPQTK4yQmJsuZB3YLPgCMwTDjp6vfSpG5mgCjT/G8fpZzZeMYWEnYJgz6oAzSLt+PDoSIgGR1lRr3sgS3D9YwCL61u84393Kuia+y9ANCComxdAAtd03d8WVrc4QWNZkvduJiwXJFEVR2zp444EDssrQHnd9Sn9c9i0x+c4N/orZZY9txAsZjzo8tZ0dVjHGnyVFppeJg5morfHhUMFJwrlQoIfnE/mfYlfvcVtc41LepzZx9uYem7mABoR20Xbo6JgkPYJ+7cTXktp1ClmfbNZUokUJS6cNBnMFdx7AnD7ixKDv7SI+/1J6K3Q8xxusmXAQlyshchDYacy27hZ2f5cV5yVah2B4tWvcxvC+tXPao1A6JsC09OEwuxgrHdUj6V1rYFbzbYAUTpalzwSw7x75fCTmIFBmGZtAqeivyvrIE0jwBB0ZVKDliml/gF7Z+hO1Ss5EnlWttoLG2gz+4BpnopvymOe8j1RBKbaXpqNqDQH+0u8aVQ8k/mv8LKn0QmXEApZphrGahPULe47z3ZZydwdiaRbOLlG1BlLB/6LQhEjC0LYTJlM+kdIrIMFLrfMJ+ho9Nx3x6s2KW15RJdd9QSqzbxm0WZKhxsO+LbgSrIMi1ec6Xd2peFt/g1vIbsXpT6tYMPXOwjKrAFEQhZRZn+0CdDv76ONtEotitrrjn9gqiHxEIAGumXNSo1oeMszNsZJc5NHc8JtJVnMY6inJoz7X8Iuuc47Z5ncocJECmiu/3R9QP8vxM8LntPVcUBWKtgaWRTLOGNb8oa0goC22b5piEnIoQkRmqXD1SmnmT8feuLrgiX+XG0VcSWrwjmjW+r91fHpMlINwKn78HXLXW/rQQ4jngV4Bt4MvAX7HWVied4/7XyFFy7POta7RnySk54tP5T/Jcvsl2qZko05m9Y2OmMzptRlso6x0klNEKacAGgZTW88AFe1XBvSrn9+8avs21mKQTefwoRmIdK4wzPwd8/UCSCZ0qdDSAkdpkLDda8IyWNedMP+I25VdONqYFyKxwzDKB5LPyRziXlTw7g3Nlw2ZRUUgX0WgLYAy3d1xPQQTSk46ZfzpRrGcKzac3Ci4uXmKH7+H3+Vq0CqTIUSKLNR3neqclPnlOQyR79RB+ILoz1ho+x+c4Uyi+YV7gjrrJzepbnYSqsB84CyL4/a6op4qKRYqc2oc5uzn7qxbYkMUQEpa+Zr5Lvf8Uz81GjLOGMplEwC1Iov0sroRBk1CHpW2px0a2i48NEIgkbbTmT22u8bXdT3JDfDX+HtwCV4Ho4fGBR5Gb+J8ALyff/2vgv7XWvgDsAL/wdi/gaMDT+OKMdYh4Jsd8fLTJnzijWct0h+EHxHh3IK32kdnW9PefpSvy2VinABpf8cfgavDtNxlvzjO+Ll7jO8svdkGmMJCZkfkZmXidrv8e0poBGr2IA76QE6ZsdvZzz59H5RG2qWRbem4pJLmc8PHZiO87Y/jExhGXJ3PGPvMxLUGWFsQwdGetNFU473Vyg1MQa3nDs9OKT282fG5bcs4+yVSdJZOl/xsxkVtM5Jbzsb2PL4WM4dC0fWJ4LypL175PTzI+sWH4SLHJeX2po0zTZ49VfK3pKAgdlbR0irNTxHP1PRnPWuzzHNy5NdeaP+I1c4s7VRZXIAp1C0KVokOfTRrauS8hiuDeZZLNSQCmnYTaht+zecSz09xHXULIOEwuLonuYeXtLk3+JPDvAn8D+E/9SsU/Dvxlv8v/DPyXwN99+KuESMDSI6Nuxi2zM1zKPs7FEZwdVXGmAjoNahNTK/1NQPT9SejAtZHsNzkjzw6zFirjKs58Y7fgHx99lX19HZmg2kJIxtkWpZhxyE5kvylZRIDJYlg2+97vHMdjS7WGtg21PmSvvsaBuOl48IlrENmEnkwUTNpMjrE4kpQDzCr+HfWTfG6j5JnZkrWsYS2vfFu0s1HbssJbPzZaBf2CI5kwSOW2LbQi5Fk0VpAJyyyv2SycMvmFco2bi5f4wu0FV9U1btbfjGE+nSgwByYuO4M/tke2jhSKebMTn/MfHX2JzaML/ED5NDPK6BJIWrcgTbEu1AxrNY1dUpsj9jwPoVPH0TakmPyxloAIVkYTLT6Aq/ab/PKdHX5m+gk+v924jFHZWk2hXFltJNOs9qFAEcHC0NYjpWMa8hB3IJRxdwxVuDz+LLfqV1hUtwjLmik5QmtHfHqYef3tugN/G/jPgDX/fRu4Z9slVq4Al4cOFEL8IvCL979E6Cgh5OP872m2zfP2MrPczeYNvTp68TrDRUHSwp0pI7C2XWVSW/dCr85H3FgYdptr3pxfNSENOrLaOrN/D+xK02M7x5vGvVjRZg0GgKpTHzGJbQcAUMmScb7BE3nJpXHDmq8mFOoCksxKfUDq2HUBBszL9IkDnhJqJhYjl1Tz/HiCmF+myo44NHccGSiN4/cGXxjUIf+gL3v1NZbqgLvVZY6oVtowYCOhLdx5fS6CNQ48o4vPpNhA4DQoP8sa0Sr3VCJ9WeTUZs5d8zrX5x/j2nzEE+NFxJissJgEmA7YS2jPwarKA+0frAKBTznOLM/rZznMdlhUtwjLnUuRY0SGsOqhXIKHVgJCiJ8GblprvyyE+NNv9Xhr7S8Bv+TOdZx3ivfpNMZUkRpsreGj5pP8xWcrKq1iZVjodtJ0gcsAagkfvw3hwIDgAuzVhSdqtIxBieVenfO/3rrBPW76gh9uVgl5/xbDUh+w5MCnAetoltbedRFIRmoTbZcdgk9lDjv+rEB6/oMB2c6SUgYQcOk6d9JRS7XNZV7kJ9bPc3HUsFW6FZMz4QqfhOfoV0gmPr+r3Rf2C+0l/ewVZjRow6jak6cCXVYIS4bh3GjOzz4559v7a3xr/yX+3/lXuau/E5OcQlu4azTJM6xTyhlzvROthxAqlDKj0gf8hvlncf/UdchEyVSdZW52qcxBQkPuJWeFPuWZpUE5CeEUwFhtochpWLoIgTBxHcuAy1ja965tw+/qr/HKjSf4Kxc3OD9eUPp6E0pYjny7BeypkCZaAq4ys3DJSFZ0iFr96SEA2pfGFT//JOxf+QR3xbdiW7j/M7filpm/ZUXwdiyBHwJ+Vgjx54ARsA78HWBTCJF5a+BJ4OrbuAbREvA+XKamPJ1/hiezMXA0eES3vLbTpLG+nPfF8sQ3DuZZWwosRAMEO1XJ7UXOPW6ytAfHglfhZdSme08plyHw1UMmXp8cFKTMnGG1bPY75xFCIuSqhdGYBUJKNgvNNNOR2x6eK1UA7n7b7UH6BTlTbkUwVUNlZCva8Gk7u3X93+2y4gXgB6pPc808z8t8mdrM22IhQpKJMaE2QuAbOApxF5QLM3XgWaTPbq2b5yt75Nu3Wdmnk2JsKgo1ZSO7zL6+yVLvRXdraQ4Yyw0yJp226GMXqSztAbvyDneqbXJZsp7XsZ+VCc+i8kVXQ4QgdQlS4Dq0ZYgUNB3XzHE7RhTOBTALp0gTRfUwPMKHRhOstf+5tfZJa+2zwF8CfsNa++8D/wL4C363vwr86sNeQyAQIu8guYWc8pnsKZ6edTO2UjAFXOdtbLuwRForvkXJE1fARwvc4iCep28lNxYF1xeCQ32bpT7w96ViaK7RCxozj7Naoxc0ehFNeJkAYrWZ+/TZvB3UBF82i513rLYo1bp3BUIRE0dgaUG3Mu4fSnGtZYapr+sPTokFbsNg+3Y4Et3P6fc+3yKlYQtvLTQ2yTmwgq1ywYvre/zguYYf3JoykVtxBg1Sqhm5nDg8wy5Z6j20qWLEIwVBlbckwmwcfjfGMUcDccmYJvIBYkq2yMlEGQHbTI64oJ9gLDciY1PbpaMloylsEd2S8F5ClCYlNIFT+gfmDtfnglvLLKapF9KtXxAWMqmNjIvJhghU6Jfh7aT9Ny1GkkZkCqWZiIxCrbv3bw3aHPpnfbjh/DjIQn8N+BUhxH8F/D7wyw97IouFBAGe5Oe5pD7O57YbCmVYJCWtU2nXuHMdN5iwKqkGnC4pHeoCzvLa4QsJLvDbtxu+Lb/rzys937yteedmhyya//2U2c6Mf0woSoqMkdpgoXdp9IKD5qbb3yuRYEq7slr9sFXG98gf5Nl8DSXq9tmFjbOfy4bsFlHpf3bLpHfbsfYFRmtfSzEtpGGsiOh2oFingGyQ86MFI5XzY/OXeGWxz1fsP48hu6U+cGQeH++2uBh/8M9HaoOJ2GJHv+Fi/EnbxcEowmD0s2FCMnKDNkeJMlphpVpDiow9T74K1Gz3fhQLs0cjlj78myznlvAFguIw/n00dsm/XLzKc/Mn+eia8YO/4ajJaYBx1nSsgbDALcnnUMpcQlzjMHVxY9i6zvnYek558JP8c/vrLJsdhHU4mXnISPwjUQLW2i8AX/CfXwU+/yjOK0g7bs5W9hSX9Fm2ygpr3TLhQ/FVd2zLEgwKoF9DMFJj8S4EYTEQVzPwqMm4Jm5zp3n1PveZzozdz6lfOmROCiSFnHFWPM1N+R1qfdhyBFBIz78fMpPdORUX1ZTtslsjcQjUexgJbkBwk0LJcFeTMTxX91op2OjWJ4BnpiX36glCtyE8Qw3kHZwlgHuZHMVsR6Bj5rfP7nIIUpqx9MzA2A7CgY3ahu9uwByJfbQJ7MbW0mzMEiPqCLxCAmRa4r7OemtZkXfM60zFjNrMmAlDmbk6ltZCLgEjaSAyMFOXNUxkrYvVKoAufgNYwXapXS5LU7KEbnhQyLecZvyepQ2nboCUI4psjZ+efoyPrFUsvYJWvVk9XTa7XZGmOzACpzvUzzvSimnWMMlrKqMQnub59d0pL+8a7uFm5eCzCh+Hd5zzg5VB6cqFGSq9xyjbZKy2WJoDXAXdJK7tB/SsuMDHzEv8/IWcX33zHF/k/+rMcpNsmw1xkVu8QpVQRoNksuQTG5Ltsk5m4nZWD52o4xJ4hZd2sHR1oc7zCMOkXMTzpjkUfYkr7XiwMbgKI6X5ns19dut1RnrTZwMGjoWrgwAB8XcFRAWSI7PDgb1JbebRqhqqtBzei3Mrqo4P7xiBdbTgtF2i9ZIlewNRHO3ckeCyJSxECESkmsroaK2ESI1As5QL3jg6yzhreHK8gzZuoOdGshQqhg5rIx1JzXMATNJ2MFyQRCTWweXxkosjQb4/iRGCNk19jBUGa4bxsiF5FGShxyKhajBAka2xkV1mLbeMfGc1iXYMjeXWtGslFgdJBkeQkMiRC9tZIDT4YDuV4HWzw9IedGYFd506Uk2DtH5oWxNO25rKV8ENcXKbuBPuet7vlZYclwUpfbZfIacAzHH02hgKS/4AcmnJpVOAbyXc132mLkkInO+qpI1u1GDIcJDp1gUL3T0atkvDp/g8a9nFjmsUw50i1CuUNHaJtksaD/al0i3F7ejS0/wcuZom9+4TjtARW7FeIYS/EGGI503CrisW3UBeRmpBWAxHYp+XdwW7VcFkPGd7esDW+MiXWHegnhRujg5Lzrn2YoXolkr/nbaL2Y6RcgSJMnwYec9aAtHEEZJZdp6nzUcYK78gRrIUVGAEWlqkNZhWodHTfQMxaKlVLCASlgTXfmasjeTmwvCK/mInNh/q9vVN827eepsDr03F3Nxd3S95p41dUtFQ6YJcWEZqw6H9QjGV2yzsHvfqN5AeFIzptJ5D4AZBl3WWos9hRnbXPbneTaiei3DnSJcwC4trxOekiyGkLkBctKQXnXhqsuDHz404uvU8B+JWp83COgJKlD7kute7N9NaSAkeI1CUasYF8Ty3su9S68Pe+2g6FkOaYtxRAD5KEVyKdrvp7BsYj9JXaHL1EN2595o3+Q2zy0ePfoC1tQPW1g6oq5y9N57xg7yhtjnG5wvkHujLpUH5tre00az0fYXwdgjPaisYiRm5mrI0C6cIwu4D+RonyXtXCViDEBmZmvJx82n+1NmCtbzpMAPBdfA+vSRQLUMRyGgp4Aa/W8XWo61JRz87OmJnOeZLd2bcbA5iCa3w0kdynTW2uWW+w1LvodL1D3sSfNPjst6CLJp7XCle5Xduv8Tr4hqNXrR+q1+Iw81qruBFmkRUqBkzdTbSVdvnl8e6Q+m2lSw2AoZiW5IRAc0m/u72b5OqwjnT9NhQfDX8bq1gmtU8PbWcvT3jSrbNYXMnzuouIQxfBKUb0w8DsVNwNeg1XLGUm/I1KnPQAQbDewjVl8Lx/XRhpzIz8mwS6cLBonPkIAdShrUNAiGsCUQtkaFtgxQuWvD6m2Y//AAAIABJREFUoeAr3/4o07yizGpeuHiN/cMpb9w7E8O3YQ3EtlCt7aRvp+nFKeBqICYdPaefQRSK7zb/2g187xrgsZIHlfesOwCOAFGqNS4UJR9dO4quQCqy9xfE4DpmXOnFbw+LjcQFLKX7y6VmbeRYX9891NyTexH8CUUzxmKDLXMmhppcss5w0ZD7P1solLnkQN/kFX2THd6M3AN3r3O0T7jpi8WV/h7ZaXf7fWb7VIZCh2HWTtcHHGK4pTkH/XOmeQepEsqlYbNYspZljMV6wu5r4/o6DYvSdXvcvazWX2zMkrneibhN5z49mt9+by0I97n9Py2dnkoIMfbvoZ8u7vCHhlvLmm/e2+DG4YyjqmBr+y5r08MIrnbyMPyENfRbf71ImRxjreBCUXJBX+gC0w8xpN+7loDwVYLU97Cet2hqSGmFNMe/K+mAJ4kKhKFkk3CgsYIXtm+yubGLMYob++v8Pt9gbncBCMUo62bOpNiiIItotBWBC5DFjhli4WmhinhfSXFRcC+szDYRSK7ab1LruYt5c7BybDh/S5GVTOQWm3rLLX9Or9It3dBgmiw0RBYK+4Q21cYtQOok5dyLlXOENQwRFmHbQqWNXuUoSAuzTLA9P889cQ3jk55CHQLXLgopaGdqkWNsy94LDMNMlCz0XjTlW3wlLPU9R7PE6DpaAYF/UZt5nOndA0E6pMO5CjXFWMN+fT1eO6x6ZGzdyRMJWMEfiK/y+s0L/K2P53zsmdexVsYl6wqfWVhpVx174ZfAc66gW+lZ13lCOWbFjQsT3Ke2DOP9Nb5cSYQYEdasFLhl3x80QvCeUwJhZSEpCsbZFs9xns3CRjOpL0PeTzCjwMe/Q4mwhKXlwl1t5mFd58wXI/brgsr6GvtChVUNscLNUnviAGN1J4EonZnSGWoolx6IbgQQlUZKDLqfCxH2yW1JiVt0s0xiz29V+vkE7vwPll9w3HboWgmpwlkv4OxixuuUmGShkoC5xOOtQQ1wA9w9uTyNtMRZysnolAwTrbWRFgtpr6s6YK07Pv3sXIRQ3MS5eXoFMA6yNAfsSIO2zyKE5WBvjaoqkrR1t5JRMP/xk1Xey1tJ+3EqYdt63jDNVNzSfeYHDxW+55QAvlJQka1xQTzPn7lkyIVm4VcN7kgSYombkg6brg4c0mSVxxAyYZjmFVuTQ/bmE27sb3BnMeaV/UmHwaekm41yMWJh93lV/1uUyCjktLNCcOhAIYzXL3gR0GolylgwA6D0s0pljlCYaF0MNg0tIKhNRcmIDVmwVS4ppWGvXi3Mkvr2qVXQL5oSKuPeT5HEAis9ayxt93D+FINIachPTmoEOV/bn9EIV7XHIfIqErKM97Ehi0BsqzCNn8kP23YJfrCgYxWkbQfOQosgn5+5c+nWcKx0d3k195xtVKfMXJr4QXMzKuqUvuw+58CSysBhXbC3t8aVO2dR0rI+mrO3GLvlzn271b6Pak8kCu3VdwBTDkFo9828ZrPIvQWwROtFBNQFOVbU8AALmL73lIAXJXJKO2LitWdKVe0Pkb4ikLThwdwvOmmMQ/1T0tCkqHji0nXmRxP2D6b8+pXz/PF+TU27ko+2daSQBjM+0HRPYgKaZJWY4DIEMzKw1wBqMWcit3hCfIQ35XfYM9f8OduMOCCayymN2GBobEvcmWbumgutYmRkyHpKJVWU95NOPkaPoxEiAulvfQn7buQNi1LC/moZNpO0qcsLaP381OcP0Zr1/AkkCk3NQu9GN6z/PgLNOq29GGRNngNgX9yiMQtXrp2UeeiWi2vMgoZFO+BtaItu5MA9R81v3tjg5mLED152jNNl45R0Y70S8FGBcKdpBMB418o964DFmwCHmRpRNTUPW1PgPakEhA8VlTZvudcMDPYewOKObSmsKqwdKIwrT24FhWgXlBzlFesXbzM6HKNUwx/sHfHH8o+o9VEEixxq7RJ3gmJYKSF+DBOw+10xlhtU9oi5X71HCFepRpFxUayzxxb74npUAAHQclmJKpq7AS139xLuAUqlIwmqSBbNHMoehG48Pw3nPYikFoA6JsadhhFTZTTL66jUAybQN63TwiIhv2JlNSMhWWMbRcZSLGjksjvThzayJgJ+IlYdaglKM7OBRKJlw5xdtG5IoYxQIWmh9yI/JLUCUrAxPJOl5ncP7rJXb/FnP7GH1opqb4OwxJuSBmPkSqarBIQ0HRygr8atFZFaDLjCsmLx7tKGH6Vkao1xdoZP8XmeGY1ZmqOVggtBA6aNk5pJ6XZtJfPGE4N8BmGhNCNVc7AY87Uvfy/r0wOMkcwJ5BRfoMG2HSgnEIJa3zEF+KCNZYeZv+0UbtGMA+1zAvyABje77XKdP1KGQ3Mn3ncA/kJ6rJvJXDZcQNHHdsJUKYSookLsl0xLi6kYG1KpLdZaaqNiNEBrGRsybAsJ3ukgH8RlBralx6dpyyHnYKTayj2FnFGbo1iR2Ng2tJqJMloAhScDOazA5QW82XwD6cE6bWuUKGPewcfNp7kmr3Ol/oO45kJ7z4ZL2cfZNtsccMhcHFH7qr1uUZQMIZRffLVmofdc1WHGsXScu/cpuZy4/XwV4/COd+QtvrsY87XvPkvh1yiotEIIB5BmwlL59pV4AlsCrKbqOIVowzsNVYy3sqe4i6FZ7rfugMg8pfr+hUbec0pAoMjkiEv5iHMjYtXbTnUgEgQ8+P/x+FYCOSh04VBJSGJZGy2cT6YVN+6dYX9ZMpdzjNEOrbYtLdXF6GsXQx4Ag44DAMNsbq1PUImZae2qOxZNbeYcijsxDTkcl4kSJfLeOQMCrLyl1LU5XKUkF9ZrkOTi/j6hwwLcWVQEl7qpwavPdn9X4yRxFXplx+UJUQ8lZHzG1NUSyVAIbRlWZQ7t4qy0breO5ru/XSVySrXOzK4zoWBH7LAUcx+edOcOi6KG9rfS+FWjJJWRiStwfIi4snMWVNxbulqEod3C83dm+9QiSM4xFBIPvwf3a2Y22E1CmK7NMuwDZha+55QAOBT4s2c0Z0fLuIjo/UT1TN6UVqyNjBaAwSmTJ5+4RpZprBH8d7/1Q/zje29y077mz+UaNPjhxjaxenC4PwBJSwhKqwGDU2a5nFCKGeucJbMZtagw/t9d3vC57C50tdBtmGuUb7jjzBZWGmp5lCiQsNKwYV2MOFOKxAowICXCOlKUEC3rL+1M2oi4v7YSbVyVXAmxjf5/6t4s1rLsvO/7rbX2cKZ77lR1a+yunkk2qSZFkRJpUYpFSYBtwXYCK87woiQK9BDACfJm5CmPDhAg8FMAIUCiBEYsR1Gg2FCcRLZMSyQjURSnJptDT9U11607nXumffZea+VhDXvtc29VN9UQUF5Aoe4995x99rDWt77h//9/xopIqT4zET2OQAmLTmZs8Bie6Bkkf+oxYiZyFvqITJQUchQTpQ58E2DRjt8fDCd0dRusNZ1GplY4Zt/X+VcY7eHB0rnrKz1hUD7LJ+2n0dawQnNgb7PSLTw8hCcKyCgZyG1K23MhBxVWBT6CKxHWtF6G8uVMi6ayU07lKafNxUhVdyVCG70t/PzUdPE9KTT+vLDXWidbLgRcZZcjuc1E3IxGU4ocKz6YHPlTYwSEUCi5Edtxt6+niLcEG+BxA9qj/86jEzfWM7bi51s4bZZpZNZgtWLeCE541O4EyBhDal8+6ohVrJXwHveztjUr5kzEIxf/U3PFPs/L+TbfaCT7vHMGHuvO3VUmJvKIpT3tJKEav0NJkbPPKb3ZJp/eObvoglFcJRBrJbuCFsY6ZZvau6kpfj0TZ3MKXQj2ma9Ee2VcF254Axw9jDbECN/ftwN6cuyP5zD+ob16REiiohOfQrLXcyPteTkjgG35HO5+uN3/Y/LzbNo+l8qc95ZL9uUjtDmbnNzMrpKLktq2hmdl57GHYStY4vpMtM8uR6DRngQlkVFpqH2P89a0d9LOW6chr/W4ak3kD0jLbpExqP19DHNFhPe1Xs3jxlNjBKToMyquxDqutp52SZsDWI/30/ZiLZy1WyJskd2tCyaFReU1Qlq0Viw1zMxBvFnRLZU5mKULIBL8/zpqreuydmmpDRVzexCNy0fUx/nchYb9u5c5lQexxTfxHCU5jj57qh861puto1vqGpzklLLHPfEOC3uJxuyeSehlXspqrt0jzoST/0pbaActhcr31gOiB3AmibhmZM/HbDj2oMRghAAb0G2g/G1RXmpbAiM7ZCY30dKBdxpTofU0Gl1/F5PvbGG/KSQ73P8OKSs5vZBoHokL/MLOmFzAXAveWjY8su9FDoYUbXViz16nbwpuyrcx3pgs7YRFc5TMkXXYsIy5gSgq4ysBYfMKz0nStrhLRzjtoHJlvHy5PGczDPyDvZ5gtHJGIHAIQoNXpzXwZLzAU2MErG2o9ISXs8/zvNohk63W/XpUozxyEJ9k6nCuk5LLeVoDSri6+YO7l7l9eIF/ff8i35qdRAPgavbTx0J1A299/fVwj88AUYTslKfeVm/xe3dfBOAF8xHekNNYMlS+Tp525XEAlZyB3EYgmQlnrGK5ShkOV5fJpGHTKwunrMlyjQEYmmeGLkMuSdqGDJM6Z6A0/ax+IhDovJyASkBZwVPIRLsgUxxCoTQ/sznm5mzIH5n7riVbNmLmE2uBQQnwkvgp+pR8R3wlYi0ck7PplBgFssMdsNaQyT69bMy2uEppevz+0QMa0VCzYiIfUut5+1mRM1LbDNniIbdphIMjS5EzFyWNWRL6KYTkrE52WWNJ8jrdMDbwA9ZHoXRsXBIQse4egvXzP8z3FB1r4n00XOo17E5GlNk2jaddK1n6UO39KwZPjxHwXO4dO+LZoSQT+omJpxACRGZV8rdWJPMcTIE3GEfTDd6ZbPKVwyX31W1s09JKUznq9Yf5+PM/z6rLeAwlcnpiE4C78j7XzWU2fC/CygC0+nl6rXYeuAsKR4bREONRIRTHK8VA5WzkdVuSI/AjXEMVaF34NPvsdOt8nG0lyzrH+Kab54UB6/fyvGcUcg7rXlk8L5x7fLXfMG0y9Kx2nY3JY1I2NmpFckWNGSnJG3Uv4jYsLmaPmoUpYCh5JkpkZKIk99Lst8zrkZ6cyz65HLiGMbh25ZkoKU2PI+6y9PqDQtQdoZGgZOyuq8sQDed85l7RxVaEvErY1Wsropca7td5uZg0JxCQr8PM0JMqeicdd/kDkImeGiOAbWj0KVc3Sj65PYtsvzA6yishCQjRW4gsyuTnIIpZrDXdWGnFGwcX+c5xyRviK6wap093pjzlfwc6bmg85SQcWO+IA5CrPoUcIYRkh6v86vYznDaCO3ODltBYyzX9AhN5wpG9G7v9ZtK30TLzmGQ6MotOPT1l3/3R4YxX+iOuDVq+QgviWYtFgbxw3tC6JHYuHMFnZRSHq5KtYkUu9ZmJl/4c4v+QqTY43kHIP6yPVJn4Qq/i4dLFrCs9pRaOG6B8pj9c56ubkgtlw1ce7jLR95MyXFfSLbIr5TAaEWs1c31AJaZIoaIBMKbhav4qN+wlPrermNSS/+P0B8zsESfmjku+CldLCiQkYyuXtzBpCbbbaq3dSGqgj0TGEGy932D3vpw/0qrX+icDKK7IGkBR1UfxfLRvTX9WNO7seHqMgJBIWZAL19SiopsLCCMQKtYzph32FcSEIWufr30jEfee1n0LjUEhxZm3Cy7gywMFFbpqQ+dZf2NdllpYyVyd8u4M5o1hv1kicASPQ3VAQ0XJCC1qL6bRegIh8WTRNKbuTPjwHQ/UQy6s+h2sfrjuNBHozvVslt4dj857gusadAHWFYvXocAAlQk6eUkZ7Akjl4bt0vCa+Fkeyn0OuYvxenk6xumKnrT0lEGyrgHQ9b7C7yO1hxSKSXMvvhawGiooVkkf+lnNzWnOaWNY2qkXgDlLDQ+fT3f6dE6kAK8QHiqc9xYk7FIsy5NG6gWkc7jTV6MDFgplcnPGUnwQsZGnxggIcopsk0x2yyPnxkNAKh4CbXkqxJ2hQ05wmYJLvEy04MGJeuRy4DnnS6QAiyMHpa3EXbmqLVm53b9VEZKim40GLzGu3cI9Mgv+aXMnGpZRdhFFzsPF6xTZJs9kn2Ql51R6wkpPCTJm4HAFroTVoERLPgI3Oe/W32FLbVMbFa8rGABoqyJntAVEu2Arq+LungtDkXkFJ1w7tkwaStmc8TCC8Qj9CQCGWdNBI8ZzDajEAEISlmv9Jf/e9Yw/fvgM32igkgtqW7Fk4sqGYkChnFFy99npLnZIW0ki0FrDVXMDgeDQvksqAwZE6HZtJHNxyl0J31i8FRO0qWCLu4+twY/JStHCuoOBDuClUo5c5yM7pxADcpM5OTFPEsqSex7mrThnU0s9gBD2nscmDN6Au3iTcAe81/oBBEaeGiNQ5Nu8pH6GQeb0/wJ/uk4SfeuWNNyg1F0NP4eyjLtJrYEIS7mOSRenDtNYjfY96tKscxjpz413tVKXfP19IZNt4z7a3ckX+ohM9rje/2mklVRiCRaCfp7AucmpUZEevmps3aHeAszkjLdOr3C5X7FTVARef+oJrDMFw99SbcGgw9gxIqpt1R5wBcI4A6BT5pv/f2UkVrik13lqxiGRW0odM+ef3pXcqJ7hzlxwsGr4c/ENMkoK20f63HZImA6yXRbNMdpWLvGnxuyJ55iLU1Z2Ecuoo2yPykyjURUon/R1O/ZUP2IujiKBqZ9tR+SiO8+zQCZ3Lc7oGOs8NYFrXALEZjOZ7DE221xWQ1676FpvHC8GMSEbkqc2md/rXbRNutF5w5D6mwa3BkJebJ074JiVEmvP6iyk4+kxAnLIM3KbnrKdnTrgpNMSYfBm0zpq1Gn3xws3JwhknClxxePImOk9yzwLAJ0Ex4/pgFbOkxiL2PIYTiRYc38atV1gMVy0e2g0D5UjDrk8RNKBJyENQYtgS7v0ACzFjDtzxSjL2S2XcddIMQEtYi65znOQgcEARLmwZFcPOgPGT+C0TBuGNhLl+xeeJ28er0UahHclnh3OeWYAG/mIW7OM1xcDFBm5Lah9RWPFHGO1S6YKibbQU2M2xWWu24scmxETOaFmhRGGUrgQa0Vb7Qkl1ty3EgveXiZLj8Gvzk3uOfRi19Nzxtp5AIUYYDGxmpDJkpEZslUo9rYP0T7HdVr14qLuHv/xI3jC6eiqE1ukEOeGpEJInxx8PFbgqTECm9k1fnoXBl6i+nFgTJv830HBJSXC+F7vgoXXA3XWSC/0aNqFBm1eIJBaBJLN7Bql7dO3AwyGRjQcywfUduHyCcIZkhBSLPQRmrqtZScua9xJfPIqkz1ui7ddGy27SUGfWlWumw41lZ4kLqc7VqW7giNhAk/NAX8wf5ueeoFnh0m5dC08EMKiraQVAzFIZ1LPvG8dOhwSgQrtKgDCdVxYacXKyM4iN9qp65ZKR9xBfC6PBcDA86MZw6zHj6YvUft796WHDTUGIxq0rThp7sTs/r87+itcLA39zPD68QavVyseibs0uFzMSO6ylz/n5ggNB/Z2FDF1i1j68GDB1O7H18IzUzJo+pvopQUUKUBtDUJUHCfdkcJ4qbfBixuag5Mt9nYOeO0nv80PX/8Y7x3vxOsNIxXMCeXCgAWA1sgauqVWK7o5nzBUbJEXmp//G2AESttnq/D46g/4mXRf+XFx7I11/9IRHmImg4Ck23EbkbNiFRf0kG20HDETR7gWWssWtELoNdBKZLXosqzjHRirWdgTCm9AAshkhytIK5ll06hnN9UPffORbmY6fEdjHRHp4fJ57i4G7BQrsnPYfY9bgOclq2J7tuROxxKkCPH945OAAShkhCAXZ6XhzvtuKZyi9EDllFLQV5IfLifsq4c0zZKgNbCTP88Fe5VrA00uLPeXGcerhqWYx3uWiRKDQdNEA65wUuLG5sgE+hso3usjfX4pjXs9Aac910Elwie7JVwsa+arEq0V5dYpuWrOJGsfdy/c/T5fOKc9v26OIZxft2rx5LzAU2MEcgp2y4pFk8WcgEncyZQ5GGms/rOPs5i56MI1wVtfCwstWWobd/12J3TJndouaMySk+auf71VmP2I/RQbouQeR5zKCUf2ltsZhIeNChWrDo/LzgY0oRCuPdlCHBHYcT9bvsq4EGgLS21Zasufi7c50rdYNscRAadEFjkF1moW5ohv1feY3N3jV67lXOgt4/elJcG0bJiKi4Tf0//XR3geGcaVA62ItNi0SYYED+RqP5sKZp5HaY6JQ2EZKcmLG/Dq1oQ33tLcX323syB/sXyNz+xWbOYLfjgZ8j8dfymCrvrZFpkokSgW9oQTfSfqCGxm1+iLMT2GzOSJo3b7dvKpUvG67HjKHs1l35G0knBtXf9BovjIeMkr24fsz4bsLvqo4RKlTEzyWdvezxQenN75IIKTer7roZsQIfztIlud+MzyfZODH8oICCG2gP8B+IQ/z/8E+AHw28BzwLvA37XWHn2Q43UyyeckQvBfEtpinVs7TQxC99jd75k1kpU5W2IK+PUgLtGlnzom4e3sFgO7wY7ZwhjDodWERufnyVe1n2+58WGE42tbs5FdZsg2b9QHbK1GfHTU40ofhhlcrl5iVr/EvWXNPY54S3+N8PhCqVIJ1zj1LQzz5lKs169n9GN/wnNOMw0D0gz/4+XFvKcgTUy25rQt32KZ8hzjEhKS6a6ojet2/JFN6xqqYPmp0Ta7s7/BN8TX2eEqr8rrvLShGWYNPzod8PYUVqYbJilyNthxz1QZRnZMaXMuyyHDTLBVCB4srnFQV3xDfBWo6Oc7jsQkRmywQ2YzJJIjuc/B6u34XKOojF/8jxObPa0zTquSXBoarZjfuUi1WmOFnvvJ8/E9601JOuxREV5psStS5CChsTVpO7/18WE9gX8I/HNr7a8KIQpgAPxXwL+w1v4DIcTfB/4+rj/hE4e0MtZTgc4uD113P5QI1TmGYN3ltLaFsMZsuRXMG0FlTMd6u+81HrRztkWVyx0YHq3eJJd9roqfo7b9aDwC+63FFSRJpGCNhQsL4nVHTAJscpGL5gLftn9MX23zqvgsF0vNs8MZSrqE6c3pgG8eXeRt7eS4BDJJbo2Y6ofMeMRcX6G2TmcuflfICfh6foOMizR2wvXJ1if1FegYCeENrGijzlDZCbTtrtyb6BiXdbdYW0kpDa+Mp7FS8antJdcGOTcfXeVlcY1fvrJis1ghsbwzFdysZp3yniInFyVjPcYtY8nVos84F1zuG8Z5w+X+knemA+7Me3x75hZmqcZscZlts8OeGlIqQW0swkge8WYyD5IFJThjBIIS1WkjOa0LtoqKxigmD3dZ1G2JN3il4jEbWjqXI1S+8z32CcZZIvHt7PUcy1+CERBCbAI/D/xHANbaFbASQvxt4K/6t/0WrkfhE41AKTd5QV6kNlXUVF+3kCk2IPwtvfxgAwNlUwpL5kVEwm6oaNltlvNuettjwFiDoaWzBhJPqg/wI/UmxmoKG0Qx5mSixCJbIFEaUyb5gDA0NZvZNT5qPsaVomCUCb4/G7HQR3xl+oj7ix3uLzf52OaccV7z7HCBEn16J3+N7ywPuSveiiWtdHJ+ed9wZ77Bz1+auAlk4bQpMLjka6iYBI3BxxGC2vt7tu4fXte0+oQGkFaANAjhJN1CSBAMemXUGe8iKkL7sCT1WjaLFcOs4dfVFXqqYbNYUXqo8+UejLIRr+i/SWUsK2OpjKU2hhOWHMpjDrjN200PpXOGC6ci5MYpmobazJEiR6I45j6H4jbv2T656bFpdpnIo9heLoxOs9JYKgbl1aIEkvdmkoHq85Hdh/TKisW8T5MAqh5XuoUEH0A7V9efUStA4r4bTEwIrpK29koNyNhg1dw58z3w4TyB54F94H8UQnwS+DrwXwCXrLX3/HvuA5fO+7AQ4jeA3wDoiRE7herERanL87iIJs1+nzk+DiuQYgvMOTuPe2+bSQh95XpqjBI5KzOPi8sZifZzc3MUPx9YbU+s9fiRlqCsdaWsZ/slu6WlpyyjxS4ntuZI7nOvKZGzIeO8x7JUbBUr+pnmcl/xzqLvhUcyT6YJtFnJPY4oF7v+O9wOm8biKQ5g/Z7+uCPkBtafkyvvArb1JJQIIqdeXDO8RrdVfBghS55Lw7PDduHVRrIyisq43EmhWuJS7a9rISrm4pRFfczcL9SpHBEkxwvlhENDLqYUI+bmiEVzRCOX1LLnuj7ZypeI256G3esMJWDXVzF4gZPacFIrRoM5RbHC+qpM+gzg8chKkfzfkSA/Z8TSeECUeqKbEBkqaWN/3vgwRiADPg38PWvtnwgh/iHO9W9PzForHnOF1trfBH4T4HJxyd4YmSgrLkkuOpS4/OdS5KDE0YmDivD6UKKla2rTNQCFhFwIjDVksu1EE5RsXxM/y2sbfb50ep99ecspzPq+dmEiNPpsNjnEi0HpxiRMt+A2rlNht8wOn9mt6HlwzU+evMwdrnCfO7wlv8sbqylfejRgU1zm39l8nlK6iTGWOdt2jykPOxNTW7hpvo3hE2g7RON25GHWnCnXhXu6PtYz2IEmDG3eIACN0hJhmxB0/k6t3a6/8tiPzEuL5VjuLfsoYdkuV/SUpRS6g2vo0me75/POdMjdRcYfTG8z4RFzcxTBVSnKz+o2oec0HR0sO/R3aITjA/TUmJf0C7yjbjG1D2JpuPbErkyUrHw1Ie04Fdqfx43CmUQMklvNCaP5NhevPqAYLjB1xuHJJoeLQfx8B0fh/w+vRKIRZ0Fx4MIvlZbEUx4FvquSGp5JWq6PD2MEbgO3rbV/4n//HZwReCCEuGKtvSeEuAK+re8TRiYsG5nuZPjF2uJ/0ui8N8Sf8W/elfUuRdiDh5mhpyRStzcttfL78oAfnV7iRB7QGActtuegCsOInWt9ODGQ28z0I5bNSXxP6DcQsrgOGtxnLqZ8/XCHZwaGrUJ4vAs6AAAgAElEQVRz2Kw4kccs7EnsQmTNlIl6yHeOn6GUEiUEE1MTGmMKq87sUJVYcmd+gY1Ms5GfHxOeJyO2Lhl+nvcUDEATkZdt89fQ+s1aB91GJJ2ghGXWZMyajD89cBiInXLIx8Y1z43a5F7aQAWcQu+jZY/DVcadueLh0nDcLDmWD1iaSUeXQSS9DNIQzGE12vg9wMEtzoW/Lw9Y2BOkyKjslIaK2i4xtokufi779NU2mppFcxw3hvR44fsWcsG02eTwwQWGoxlSGWoPrQ73z11jN/cVQoB16bzwXiO6Wg8dY23CJlR44+QYseZ8ZwP4EEbAWntfCHFLCPERa+0PgF8Evuf//RrwD/z/v/d+x8qlYVzUzJosLuTsfXardShxGMontprE5XI67hYpRSS5jHPNMMuRdYsJSMet5lvcRmG1S/KNs8tuYuiTdmIlySiXIDOxxHjBXkWrmkVzGN+XkpTCMUo15pgH/LPlXT5ffZaXNzJuy1scmlsx1nfXb1g0x3xV/CHKlvTlJmN5gdL2HOtOtNx2R2opaUTDm6cZN4aC7dJ5LcY7q+4cWu8ofc2Bgrqtytchv6HkGCsC0pD57rsGh4hrrNMzUMLQ93kZa+H+osfdRcbvL3+PRs8YNlf4D+Qv8Nyo+wy6ElyCW/OCH0wsf7D6523J7hwMTNiplSjjcwn/pMo6zy3gOmqz4D3zbcCBxkKfxFRRqlBO/uxZ8zJLseSO+iGVnqLtKnYgSlmEU3HCQbPNW/evsjucMuzPqZrcl1Tb3ECYozbxftISeLwfqaFYA1811rrEtFkgREaRbRDYku83Pmx14O8B/8hXBt4G/mPcuvwnQohfB24Cf/f9DiJwNf3c7yAhCQhtyW/d5U9DgtRtTNs9K2HOZKJ7qqGUmtoKeirvqNCmFr09N/e3uXG15DOKQihy2XduJs5drMyUB/ImlZmi/MKXImOk9lxdWh9FoRGHCnSiot8R7/DmZMCJvU/aTttiGGYXEUIyrR+4HcosmIlH9NQmf734AqWCeWO5Wc24pW4yN0fM7RHfmU0o5JgXN9o8QLrQOafBq7WCOuEEhFyKFK082MrjAgxEotaiyZjUOd8+ynm00jzSczIUPVHy8XHB1X7DjdGMN04yvrU8RJuKzfI5vph/nmeHDalCceA+GC80e7Iq+KPDGXfVHSw6kriCPFuhRjEcCEnZ9XJtUC+2wtAYZyy1TUBAHuzjpObXdAoIJKKa99SPkCh6jNGy7lYLaBO/W2aXscy5NxsxqwvG8wFLnbcegHDVlXCK56Fe14fyG2Qwwhq8MpRNwswVK++BfpBuVh/KCFhrvwl85pw//eJf5HgqqYPGR5ckRJ7U+jO4UOfRj9P4dZTXLr5C0FcDQqeb92ssmgpbumO2N1cKhbGBfWbAN8Jw8anbeZTIKUQ/En+CGOTKOIFLpXKO9C0O10hDwY11oiJ5NApusjv+weW+ZZhZprXguO7xwPaZc0RjK47kIYtmfCYDnXLbU+8gjPOSroF/YKwjwQRMf5hicy15VGV8f3HKfXWXR+YdpJD05TYXFp9hlLl7vF813BY/xNqGTJRsF4JStnDmMAJP4bTO2a9y3pXfZ9Y88s/beVarhAhm8CrFiQZhZ4740Cm91PNCuw6OQ6TG0HVOntp9ctn3HaFzUuGZ9DsHtk8vk8yazOVIbCve+rjxOAOQzuv1BHqTSPK3ua3EA7A80RA8FYjBsHhDHBQSftB2ZelSZN3Nqs+ZqOlImXTgsuQvXb3LxvYxepUzWX2a/LjvJbwUlde3a8+r1bqTXqEmUz0qPcXQil6EElMpe1EIpAVtKAwuplzaaUSuudxCSjuuKeWIQgw41Q87JUYEkZhSqCEpZRXgj49cLH0gD5iwz0K7qkWpxlwyewxzdzbrC7s2KqoAGaFaYM9aadB5YxJrYWlUbK75cNHjy/uWO/aQI7mP8f0BK6ZYbSKENqNktycolGZS59TedQXJwfJH/M/1bf7m6m/xBdlnt6yi4OmkLtivCv7Zgznvyu9SmSlCKDJKAt6/kKN4n0LlRuE1nRNPIO1FGH4PCbMgNBJSzkFQhOjae1RgR0S24lQ/ZCC32cgvcqzvnvEIrmRDrg8k/awhi4tXkvmGONYKll4DUgkbeS7a6110PIZkPaTDWsFcKyrjNo9cbUQjaGyNMSuU7COEpHkMfeApMQLizE71QUZIRLljuJFWEdJdRWKR0pLlNXm5QmWaYVZ3rHjKD+9ARAlAoe5ouecuI/w4ikZoRFKZabtTeTXdlpbq8O41yzOTCZxQRp8xV8xlChQlypFqrOVCUTDXhmN7gkF72TGHVZiJJZUeeJUlJwCaLnQRwD5nzrlNyFkbym+Ce/M+p41k0Qj2K8Fd+4gDcZ+ZdtqHsdW3N17Wahoq9peWQip2C9/iVQSy1opVs+KkbpjUGVvFCmnc4ri7KHnzVHBX3WZa78ecSgrgalF87YJPE7zt1bYLOn3eqSxZKiaTdiiSgo6WYNpqrGaJsF0vIiSZh5lkMzdxTqaJ1lyYWFZNId3nPQP3me7cDkNb4cqkuNzFGSasn2vnVIPjeEqMwAeLh8LfQ+04cAlIPr8OfAk3Xvm+A1IahLTk5ZJe1sQd1XhdO3AJtuCGj9QzKJEzaRz0oTHLOMECrTcoCGHx7qHsqAOFCb9oDp1HIUtqs4iipeHhpUpF6RBInjMvcaMc8ktX5vSU2y1rHy9r03BvUbJ6eImlnHFq72GNcwlv5z/kueqnqa1AWYkxbRu2eL9CvB+Xi4vxnQCmZOl7GGTC8s0jxQ+rY97im9hwz2zrsSiRU4qRM0ZeXGOpT/iaeIv56fNc6TvOYiEGpFXvE7vk3mKDa32JlHBvUfL1Q8sf1v9Xx10PakNtr8i2r2ML2mml3tz7225QafzfeT4izKPGe4DtYl9/f0o/r/SEWiycAcSFG6Fn5G4Jl/sr/0zdfW+s7MxP6z3axsgzJcLHoQhTJmFtHAS+9qrM2jdjUbKfJKwbntSi7KkwAmGcZwDOSwam8mIpSCWECconvVLgSaE0g3zF6ekIYyTbVx+SSUMhXIJJCwcgCaKemVf5XRifYHksTOP8EVzNdGcIbbTCxES2jU6cOo3vdeDlq5UoeZmf5NX+mGeHDu7aU5o3T4d86WHjG2Zbv+82FELxGi/xmf4rrIwhl4ILpeTqQPssfZdQFUIBICYCg8x7UCgOk/et0wFfPzS8xR2O5QOsTndGFe8bEFl8sfmnUFR2ynv2gC892OFSL+OG+hi37MscmQVvye/yXDbixrAml4ZZk/GdY8UDPYkiKsa7t8LH/fEa1jAY4f52Y3nn6nf7BdYJsOfxLeGjXBgtYCyln0fmaOIduAamA0p1lsB2Xk4gJcZ1GpIk8zxUywS4vJlwxKIKwWktqag7595K45/t0Lw+niojEEYaA3VKJsFI+N0+KAqHbHdwmZQ05KJbelHCUOY1i2UPawXb/tgZbheQVlHjFWWQsb/ctNmPC9SdR8rSaneg8xKLYWfQpvHH6HsPo5UkCy5uELtI3ddSjngx3+ALe3PGXvRTG8n9heJLy/+l3dVEzkb5DD/DF3hmqLg+aJg1EiXgUm9JkciOp/dT21YDIJVlCwmoQmqCMvD9peSP9b+M7dizBIUWDSfuWlbM2+v0BqJmyQNxkwPu82/nr/HKxorLvYz7yw0Ws5e41Bdc7i9RwrDUineqKcfyiJw+tYGA1gtbwblNTMMzSIx/KA0igvGVcSHjOQ8h72OtIeDtzzxH/1zO+95gBAMdORcD+nKTXHheQPBan5C/kmvvs0kCUKy9z4JbH/54swZqkYaQAavQljif+sQglrMNFpLJGsZ5SDcJSF/TDuiq0JIc3ORW0pApQyY1VZO7akGdIbBsmDHH8oCaZefYC33kKwfOlQzcAadX150IKSIrNQyN1QirIk899DNI3cwQAgQxi0yWGF0zzq7wb+Uf55WxZpTXsZ1YOuKk96HEkVxwayaZ1Rlfrt7llAP6YpOXzQ0+Oh7wykbFVrGKBiHIkUthHeNQQEELPglQ40WTsdROjzETJUK20mfaK/BqkbNiTiDPDOQ2Q7HJMUGAxRnYFXN+cLqiNgU/tbPg2aHlY5tjxvmMXBiOVyXHtSJDociQqERNyZVzM1n6UAQguPw52lbJ/TRob1gDSMsdwz3HXG7F11yFpu0nkXYX1rZ9Zut6E00sSS5beTHmfNR+is9uDLjUr8hF23os5Uto3/kp/NNWMMyaiLeofC+CFCGbnbM+rIVpbZ08HcQkYLgOgQsNpMxo9KMzcwieEiOw1JK7i4Kr/TZuCYs/tblhh+qECP59nTZPa8dvdfeTWMy7wrkNO3zTAQ1Z6xqNSNXeosAdOLsTtBOk5RakXPQsotVkjFPde1q3svUkcjlgw27zzNCyVTQxzxEefE/BsLzKsjmO0tKlGnFZDqmt5fZyxSN7k0pPqNSUOVe9AbPRAIT70uLXk+uxzhcJvPelVqx0Cm+VkT2XNm0JhsG57N5I+NAgvM9QM2HJ8SqPoJfaCKZ1TqUVdxcFxyvJQOT07YBTUWJ8U9Wa0Ca8j0Rh0G0PwsRFb89ToWRB4fkCbbxvUGFORCShXzhCdjy7rp7gWvnwHO8vlwM2RcFer6GUJuZa1keoXKUjELoyYWiEXPtrG+qGn0M+pzIWk6QOQ2/M9HrMEyKCp8IIPGwe8buH9/nPru7G16IrJFrCSWgtPkogsJVXDw5hQZdb4HID2ggaLamanCJryLMG3bhLLwlGQHekpoWQZKoXY0bnMjYY27rroRtNqj+IdXFkIVt8eMhcG6v9YswxvgKx9GCXXPZj2euSfJEb7PHa1iSy6bR1mnoSy16v4efUF/kz9Q0OFj9AiIxrvMLfeXbB//5eny/V/zdCSAo1YlNd5kbe51PbU4ZZTabc/Vin8aYJ1VAJyKTLYB+uMua6O4tKOXL3ABX78/mbgEUzMwfMaNushbBH25oTdcyDVcH9ZcGklnz1+JRd2Wcjy3mjPiCzip8cbpMvxszZohIlKxbUYk6hRlyyzzGwAySC1+2X0bZyJbc1+W8jnMDImD2mHLAyc7/r11Q6gRYn4K8wUoWgljNCfFaubXnLBWlsRSEGXOEl9noZO8WCzIdhwdgGPIar9bf3MhjiQurIuxDeg1jvv2GSOe5yYDDTmkY2vuuQ15bICoditJK6OcHYrqebjqfCCIB14pC0rMGUD2CEJccl9/LEujoDQUSrCeGkqZSHqqYddLSVrHRGpjRKafrbE0ZlhWLs9OFFTpaVBLmws0iysy4/dOv94fXQHksksWLYJYPrHF4DN+GUyGOSa8oxU3PRxezWnMFCODSkQFgX613sfZTSlHztYMRDfUou+xRqxA5X+ULvWW4MdewylCYD0xEIVk1CBIpPx7YeWQQqdeTYM0o1Zii2qW3FSXMXiYkJss65x0avhjdOBHOtOZJH1KbhZFWiRUOGYtZYxrniJ+0zfLO5w4oF19Un0DRM5BFTTjAYrGmrBO1urmKZsDYLFnLipdncMhT+Wa6LywbjLlAxqRlG2m8iNfYhHwAw5gJfHF/mSr+hSMLRc6XbkhxXioMJQwlLKY0DAiUbWzoHQvWg8mXY0H4McCpJXh1pnVOyPp4SI9CONNMfRuZj/Vwacg8FDqKXBRplLbPaubu9hCvfOa4VrBrFoIBMaYrdE4blEokgs1ksbWlq5iIsWH8uPlPvflbR1Q27X2O6yr/utQWy4yWkOPSWExDkyISQDhGIYWYOOBGXqcyQtMlHB/fgG34KkXHdvEBBxlcnRxyoh+RmwJg9rplLfPHyxBnOgA1I8BNp3ToYgKVWrmOTenwOwlpN7TP1meyRiz6FGLCr96hExZG9CbjuvLkcRMxC6BmQ2YwGzevNPSqWzDhirk5RZAzsBtJKThvDtb7ixrDhRw+GnMpDXhCXONQL3uAdV2JNsAEh/pci65QKG1Mx5yhWFFpsgFd0oo51OOMNujMijV9URfQsggGrzSKGfVJkMTTZMpv80pVHaCtYNhlLnbWAn7WQIBjWsPG1MnrOW8iEQSrr+kYK29n92zngPlOFkAvddq+qT305esj7jafOCIS4v2P9vDXNhenUuK31CELRNtXU1iGyWmIG8XVtJc/fuMnoyiPqozGny37nu2fmoD0PkccYMLjp7vWuNyBQFGoU5asDuKiVBPfttUUeY9OQYGvLT8RFEjLRj8Rdfvvmx/nUdsandyZYCxrXgPXhUvGl5l8zr12iZ9fvtqdelz8TJX+1/wLPDjW5XLQS42sl2JS2GwxFZUJc6bLPuTBcG1S8Pe1D3VY82oWjqVnQ2IrbcopBx9BG24aeXzgrO213a+nSee8138ZY14fQqk0KMeChfRcEPGKPXvUiL25YBvSQQvGGvclCTqibRWd3y2Wfnfx5ZuaAZXPMbv4iUigO63f9tbQKwq7Nm0NcLvQR2qwinkDJgtB0xj3rFl8QjtMaGANUKJsjRM4WF3muGPHZn/gzHj24yLv7lzArkST33L3OfPQ+11n0BHpK+82rnbPWWvBemaD1kFPpPOdBQB9fVUrDWZyWQPBs/lKUhf4yxjrmP74u2pbkZzoQJy5/C6SgW1cJrwGqaFBlzemdPSbLPjpATNExgRUsfufzSbIoRaS5HcLtFJYsKZudpfaGRppuB1feqyCGCeAEQaRQ1Cz5kbzJ7uxFrg96bGQNShpqLTmpBafVLQAyNURb56RXYon2oJHLPcuVfhUNgBQWa89CgqO4JxYjWlWmMISAftbQU3SFTSOsyKC9QTzDr0Bj0EhaL8liMMLQWMd9CHV17RuqBC6FxbDkeQA2VcG2uchDccv3a+jeV4EkF6UHamXuZ5vTU5vud3rRlQ8NTJRfOOtksQ457BzpufURXrtoLnChJyhHc/ony05F4MxnbKvEJEIyUJk4R51SUDdfs95fI/aEEJZCuE0mlkPTe/M+5CF4yozA+ggCFaHk5yxkW77SxuHRoYVe1lphva+mEq8heATvvX2D/NZVbj7a4xuH2xyLfWYcRTZaQKRpLyCxnizKpOMOBGCJsa4XQC77lGocOeZdgpH7P3ATgtBIqKOHob0e/kjtIZBMzH3+oH7El+72+LXtT3Ktv+K9ecF78xpjK4psh0IN+Yb4OtYa12nHu+jbRcNOEXa0FjWpcOWptLFIyENLYen5yZsJh24DKKRhI4ed7DmO9C1WSe+DFIATkmkhwapEHu+r8cnBXPZZiDlaNozyS9RmzsrMWOgjVmIaPaylPmFuG1Ym469csPxEc5HfOZozEfvRwIYdvLEV+/Wb9NQmW/kzzOwRfbHJT/M5NjLFOBf0lMvuP1xYHtU1P5RvMVQXKNSIWb3fMcQh6dtKzdURxyGTZ+papGusMPzShQGvjE+59aPnqVYF/XzFvM6BLAKEUgMcEJhZ3Pndc2m81mZL4fb0eF+pSUV2ArP2QpkxWDnOgDauwhaqHh0i0WPGU2METKfI0cUJFNI44Ip/LXDYhei2wRIeReXIRZKefz3Fxh/NRwhhOap61Eawa7eYyCMWeKkwbwDWNQYCW8zF+U6Iois0Yrz4SMjyzxlkuwzFNif6PrVZJN6Du+1pecolEd1iCpgFZ4g0RjTcmwuMLdhfCqa2Qsk+LxQ/wyWzwxt8l5qFn5S1h0CfxVWsS3e5e3OWkGKxa4Ab2C40nxQv8E2leWTe6ezGYbFsqasYDCeJll0q1tKTY66ZF5BIGqE5FjnNmpBri6GQ1BhOa8nFsmGYtR7beVgN6z25WiydQWTOfTvhZFXQW+X0pUIIwaFesBArclvSULX4ANuebzieywk5kFIQgFGe+dleuyM0fXLnmBvbB9R1Tq0zVkk+YF2hKSQFQxXAXbv7PfOrIMiKauE9BroguvAZCYwyQVn14nn/uOOpMQLQxqEqKVcFZmGRxPy1dVJVuTCRbgotWlAb6XDvSnbEMQyCe7MRtZFMaucePtPrsb/aZCLuu9Kd1R0xj5AXCPx1RUsdHchtapbMm4OoSzDOLvvyn2aX61w3e3xHTliZ1ntwYqQmQmEzWZKLNoseat/BJRVI3los2K9yFkYzEVOKbJOfG1zl1c0V797fpLEVpexTmWmcyOsYgCi1Js4StsLkNIBARC8gjGuDJZt5zqM7z3Cq9juuv6sC9Lmhb6AxfF8exZKoc/dd8nSLS3y8v81pbZhpzb0EcRhGaxAdHv5oJbjUc95ILSqsMWSyF5vHhs82Vke+vxIl2tb8QH/ZnYdpKLMNl5swM0o15rJ4mSm1axyTiMK4eegRhQlxSKAYyG0KMeDAG4GQwynEgJ/52PfYvP6QN772GlWdM18VsTLQyoK5+5wKsoaqgLaSXGqEtCjbztmVWVcy7q4NJS0XSstw2ovn/uOOp8oIQBsLPY5GpIRxrl2oECh34U7lVcbEVgc85GOqYDkzX0rcLRsKqXj74Zh7tKSfUo08UchEmnH4WwgT3HfksfYfxtwcRXfygNucikNHgfUTDej0Mwjx8EIfdbjpId4G53Jrq1Gi4Gd2C7S9xFL/Ci+OVvSzhuv6Gg9lwQn7DOQ2PUax+1DKptRreY6ULCQ9arM2kqV/BgOlI6qw58uzv7w34FPLL/DdyZIH8pDb+nV3TabiHfUOANaYDoEqGLNH9j2+WjXUoqIWFUszwVrt+y8abzBceJHJkgfiAX8yq1jpC2QSFvYkIv4CHLsLJ8M/oxY1GAxv6BDk7knNhEdUZnq2FBjLnm1VJ1KOfRnO8Rk0xmoymSG9wRRKszU+pWpyloshQXId4UqwtW3RgaHH5lKr2Dw39IkI6NBQOjQ+HEhh3eHZKmHZ69VsiBIpyugJKNFzpDQ9e2JSEJ5CIxDGeSjrNgmYtuJy73YqNDYSYPDgjND22aHfiNlWJQ09oCwNPQ8YCjF6KO1p6zL24W/um/QZZFo6Amw1l30qPWFm9mMsGeW/bLefAeBiOf90U0y+FE5Hv8EgheBafxUXZCYNK60YyYKpHXEqDinEgKEZOXaaUSjRQqrPG+shg2s62sJYlXSJq9z/7aWNGdf6GQs9IJ9f5EBt+95+NRN9/8y5hyclrGKpJ9xj0sFftGjLFmsh0EgUU/uIUx6yu9igL1VSwkvd97NJwoAMdNeXNoZtW8tVdurzMqbznjTud+fTlhUbW2GEpkV8ugRjSY+mybBG0B/MKU7GaCOQsm0pbmjva9BjcOfSLmohBBiHAUmfz+OencFtaGXW0FcOGqxNCDlzD+lu+EvtQPSXOR6nbLMOGw4jyGYttXIyZR4pKKWICcJcarQvzQTqZiENmdecB2dOXCtrFd3KNIGXGoNJc49AKQ2j3XmcMcikO4ZOPQZBB6ji0H3DjiusfKfbQvQpbZ+pnHLQ5PQzp9abtvxamgaJ4Dn7EofymCO5z3dP9jipFa9uTgiqNuu1aud1tWo+SlqGWR1xFtfGx5RZjRDW7W6rgpVW1EbxxcsNR1XJa7PP8MdHU77Pn7XPQih2xHUumgu8LX/gvZzszL0KKLyIdffAlsZqps3DuHN/S30N6SsvSuRe5q1L59a2FV11bbtaMVcl83i/w5IL4UxsQosrDY6zy+zYy9zjTacOZQ2Qo4SMnwneYSZKXjAf4UY55N07Gt0o9l6+yWw6RB0Huff2vgcVpDCHc2k6DVvP48so4chEftpEabEwBM5L6ytBT20x1fcwdoVJeme+33iqjECaJJHJhXe63VpxboIrllqE9YhBG0EwIRQIDyTWXr2HcLHMuF59lLvyTZ/ASxlmftKuJY4cqShQfvNO2TCeE630FRBdzHgtfuKF8CF8ZqQuoJJHU4kFBkMpelR6EDX9Vk3GaZMxY8lCLnCnbVBk3F1olFB8dPOsdFgYcdL5fEDIvYyKFYOiosxqpDegrozlrq9Qms3egq2qhxRj7i1GLKqf4L54JyIJG9GwECsHxoGILUgJVEFqLRX+7NJ6neafy7c4GDTW4UvTUp5d8yLaZKR7jsaHBeeJi3Q1AzxegLS+7p5rF0NJzA9dL4b8xJZ2uSidgZURvo5ow1b3nW0VJiQFUwMQm7QkzyulFz+egwi5FBQemIVNypkfIEfw1BgB1yyqO8JCTQFE4OKr6E4ZicGVCgtcf7r1fgPaBDJF+C4XDoT8wGtbmivVLr8zOeCERP4riRfDBA2KwjLGmRmFGrHSU8caTEpmwtdvQ5yZ7vzWuuTgvPEeg08WAtzQz6OQ3JcPHYTYdyRusor96ieojWSU1zxcFuxXGffVW8zNEUrmjNhl0+zyunyTw9Nr/OxFRaGIYQH+XiIsIkmYNkaihGFYNlzf3Wfv+j3uvXeNxbIXJ2GRNZRZTaY0V6/dZTEfMLx3hULu8uJsh989aDgU99HUTMx916jVtHx/IRQSE5OsQQat+8xbOnAoKYZErQuLWqZgOsIiD8pN6WtnDIY38EEVKhieTHrZMrHyXkeo4ug2F+HDnJCw/dyFhr/x0e+ivZT48njEYtmjMZLC95FQ0iJ1i2dJx7ruRTpPwT+35IXziHXg+miM2OWY99zC/wD4gDCeGiNwt/ku//j2Hp/dGvDiqIo7U8hYR6qwCAvbLX4lHJQ4jEJqv/vLaAhCAqwlcRhKKSIZY5y7Wmpuy2TXTqnC7sE3ZklDFWPEoBjUYgxUAmXNMdZ1Hk7r5wFw4zLNfZ5TH6cWNQd+8QDsi0coMkrbYy5DDT6nslP+aF/zTL/go2PJd08U7y0XVHLqQ46amZTU0uHkD+Qh/+TmNZ4bKV7aqNgtK3rKXWvQE1h5yurF/ox+XjPuO5f3ZH+HTGkG/QUmlmRbOmy17GGNYGdjwgtGMi6G3Jpf5tZih++IbwEudAqus/tOHROhwZ0PRvNa9iqXzR4X85LTpuHL+l9Etz8FJgXp9pS9CK1bH95nEibd+giGIE28hlHZKYdCO4CRylnoI4JsS+cY1mCFYT4vCiYAACAASURBVLes2Ltxh9nBFstFj3d+9EIsQ6f6DUK4BW8TrzS8vj7CJqeES0uGey8SDzktFTZGUirYNRd4oIaYJBfgvKYn+RBPkRGYV7f4ivhHXJ79p7zo9efDwg8IKyHdjQ03RhtBntno+od2VZVRPmniD77G2DIICqljNWHg0XiBuuomR8jeeyFIckdpXXPfQ0srN9mDRqGJXkP42Z3vCuPlxTJRMjbbvFxscdpoJuIocgeO7QMUOZd4tpNoa8ySb8g/ZT7/FHu9IW8tp7wpvoM2ddwFl8ZQiwWlGHHCPv+y/h6fOvl5NrI+w8zJrSPa3gGVlw7f7C8Y9paMhlOqquRkMmY4mJNlmqZRvhJj0Nq9v6kdYWs0mpEpzeZyxsemQ5QoeH2RAqxa1qTwFZjQBchdl3P5L5s9PjLq8fyo4aDK+Oqh4xu0sGuJCTu2KF0oQYuQk8ku7Ty0FuILZ0tnabycLvBKT6nFnM3sGgV9lmJC4OWfF2OPi4rBlQNMnbGqCt49uhCz9iJ18f2ibtZ28iA4GrAE1ral8nRI6AiNhA3S4nAzhYRt0adUI1Z6ArFJzvvjBp4aIwABm207rn9IAjZWojy3O5eu2VPnJvvpYhCxWWVwY7OYnXUlmGAIjHVWVPocwi9u7vHu9AL/qvmjMw+8MhN/PjnPyE/Qt31+xJ/7vEBGS0zp7hrWGkb5HoqcqX5IQPSN7CYD26e2lsoaViycmjFl9AiO5D6Vbev+4bjHcsKbpwMWoqIUI2b2IOLjrXWCpxVTpMgYZ5d5T9zkHx+VfPb0ea4NBry2PYuNQja8YMn2xoQ8r6kqt6MWxYqlB6Ao1fiGGYKicCKtZW+JzDUqa+hXBcNlycdOttF2G7XIvceUx8pBcLWDVFiKMsxkj7viPovpDrfmA5YmdFkuI2LPYsg9/MsIHVmXyxgaBPkv2aWErysKi9RAla7XYOxT4Fme1jBp7uGIREHRqIVNG6sZqW2e0y8wqxsO3niejcuP2NAK+baNuQC3gcHcC9nk0jBQruISOyv59mzBGCAgy5oIK5Y+b2Di9bg9LXJr/Ps2MsP1Qc73ltss5FHUmfgg4ykyAmetbFcKS8SkoBLGNamQ6adbDPx6pjWUaeJxkxRRSlne6zWuC9JMRTx8ELRMy4JG+MTTGqEk9C9Id6BgFOI/X4J06DeL8H8xyWeMT24trSa0v2qTX5pKLDlpGs8VOFsDDqo3QdegslNOzUPeW+3R2CHPDguGWcMwqxnmKwb5ijx3ScAUk641GO9tnTeCMIvKNEWvYqs/Z2u6QeyRgOMVtBDWbk+FkM3PREklFhzIR1R2TC1XoNtEHbifDZpCDBjYDQC0aHgo5n6B45GcZxWf1o1z+F2thQsBCWq9cKywj+9FUdg+e3mJoWF2OqIYLFjNexFzkSazAz07VAQyFUBJbaI7Td6e17L9SUNbQaEso8zJu7eG7t+46kBam+0qqCx9w8tcGhQ6QQi22e0QR5VZgw4y2UnspYRBBTWu5LWNzFD7jO7V/pJZMyCft0zAFhjkOv9qW3Gz+UabwCJNSrXxZUwCIlmaSaf8B3DKIVYaNvIxFkVhBmhb0+BANEGCPJeDyFJ0EOWKKQfcFT0OuM2iOY47Zhp7L5pD97OqY4LyDft1buptNo9f41o/45Wx5srmEXt7+9S1261c91yBtRJtGqRsJbGsFdR1jtYZq6pESu2k2zdPKUdzLu0ccjAfsWV2mUjJggnrtGvntod6fUPuk5kLO2FmD3hk3j7j+SgfolVmygVxjU/3LtNTrvHq/7m4T9NUMcuf1vBTj8x1HnLlxUyWlGoc6cit2lPAi6S6hC5cWb+/W2abj++4EvODwx3evHuVxiiUdEStMLSnZ/eUZlisSFvBB1qcY4jKuDmttMN35DLIknskhWh7bqhkYzNWMMw0ed+SnWYfCBuQjg+eQjxnCCH+SyHEd4UQrwsh/lchRE8I8bwQ4k+EEG8KIX7btyj7sUeQFQ+w4XYitlyAdPdvwRcy/j39XPi9e/5416zNuo5zzU/xSa7wEpkoKdSIQo0IKkLrhCKHXpMxCx1i4M73EGSy65hDCGO/0hzWK1Z27tqLsTxTA0/5BeDyE7ktIrYhDOM9h5Ahl97VDshGi6GyU344W3B34cVDpEVIn3cxApW5hS2lJvcqTO4cXEswYyRaS4wRNE1GvSpoqgJTZyjVMC6WvJpf5LK+6iHW3XNcj1Ebs2TKAZWd0iowtSW+oGEYEoQTecTbiznzxmfExQV6ajN+xoVnedR/iN9LKz+mk2cgfXk3CMEEnkAmS4KmZGB8hnucyZI+JeNcR3Sf02MMJWj3nSF5PcpryhT2btSZ3b7t7dDO53R+nzfaEq9r49dTJpY3f5zqwF/YCAghrgH/OfAZa+0ncJm0fx/4b4D/zlr7EnAE/Ppf9DvSsa4rmBqDbnPGVoIplXVy72tdrWBklOhazM284fMXBc/Ji2SipC826YvNZLfNo5vpFljZJo5i+bBrBALGPPQnSFl4d8wJ9+QjartgZeadzkOp+nB4LQBm+rbbc76Nt+uYpwgJyPQYja14nT/l1ryOk8sagfHdmUNYIJWh7FUUxSre88a3GQ//jJFUq5xVVVBXBSrTjPtzXt0yXMmGsaT2pNHYikl91zMzu2jMsFhT9aZT/ZA3xJ8zbZxO4NhsMZDbtBJgtcutJM1fAykseHWBcQltT4igQgQtWczlDHqRohxCw4IBA5EzLpxEesrDSDEBocfARr5ikHnYuZe6034OpyPS3WUbLoiECPa4Bj2BTdjPmnMK7e8/Pmw4kAF94XpBD4B7wBeB/9D//beA/xr47z/oAb/U/H+8feej/J29TXaKOtpKB2cVgIpTOjbOsJDJNsuaKueEjjthRBhnUnkolCbzUk7DrKGfNXz/ZAO1ypmZA1J+QMqpr82cXA7YyW84QQt90il9hcnU+CRNkC6X0iXOphwjhGQpZlTNNIHQdhmGFhMnK7jFMRcLB2tJymUuvu65/n5cAeCRfW8N0ahd8kkIhnnNoL+gGCyZnIwd8WrVOm5lr0Iqg/T3VltJtWoNQXtMQTYbUJYVWaa51FsxzPqYVdsFCE/2SSW6tA5EnJZVGYFEMQwr4z0PORhtK74u3+CN6QYnYt91BjauM7ASOfOmFYcRQiKsig1GYi+EBCEaWo5JmfQUTL4r5hKEoifH/MrgNfZ6hsZoJquSVg+z9URTWTEdNiSc5xXnIS0LNtUQCKXwVO8hmMbzenOEtaCEjTwGa6sPHBJ8mNbkd4QQ/y3wHk4v6/8Bvg4c21bi5DZw7bzPCyF+A/iN9dePlm+yyI9Y6r/ezQ3gbth5mntBkmk9AVj4mviKLLpa0N7cFNUlhUV7gzHKG7YK2K72uGOPcEKW7QQJPzvWmo4CFevvUV7oorbztYy1JysJqMSSlV34XapVMwoYeJVo2wU0mEGzEHMK0UdIyVQ/7BxbICltL4YonWRZEMcUgkI1KNUgletnUNc5JqjZCEthRDf5anxIYAUm5GWEYbUqaKQhz8PCtchzNi3pcxNBZLU9564Me/p6OGfh6z/BrT9p7jDxIVLgAITrD9JgUmaeL6LafIwnb7UKUbrjraSch3UsQiFHDNnm+kCzmbtCdW0kyyZz3a18KJBKiglhKVVD7iXvtc8/1Vri2u8lbr846/bHkuEHGJ07aM374gPC+AsbASHENvC3geeBY+B/A/7aB/28tfY3gd90x5LxbI2Zs/Qb1zpHIOgIDIq0A40HXwjZWk8flN3Ye0Ce10xON5gsBhwvB/6Y9ty+9tq6EuJmueBzF0652r/Ab+2vuKe/HxNbYVcWSGoxR9uKY30nupohgRR093JbcihuBd/FT9LagY9sxaF+1x1zzX0FF1/11TZDtjm2d6NoyNJMWDLhp8UX2Ctz/t/V16IXsrJTlnbCxN4FXMvukGQMSjvGNpRSMCzdTh9c+0YrGt3G8FJYpDTUTYbxcWzIclstYhXGWIHRimVVMlv22K8KZo2J4B5FHqscLmzyOgnBQxEArYAHcIZuHcFCSEo1aDkd/nNStt5ElAkzKzLVi8hDcD0GglaADipDiVEv5Kj1uoIoqsc2PG9e5dl8g828dnF+Uk0Jtf10tzcIeqrmEy+9yWDrlN6FY6Z3LzI9HvP2netUTRbbw6c9IkOoqxFUOqMxadjn10liGHKRNOz1lY8PagDgw4UDvwS8Y63dBxBC/C7ws8CWECLz3sB14M4TjnHOeLwL47yBs1YxdZ/ia1jyvEZljmyzNZyxOz5h/3ibWV3QIGNVIeQOJNDLGnY2JgzLJdfHJb//4AIHvl9hmHSlGKHIqcQ0xqFhwgoPYpEic58RhlKO0JQxDpUij4ak8QIiUuRx5wkLoPZ6hEsxjfFohL0CIVYKJcV0tCXFts4dRErcz1A3GbpR6DqR1w4LXRoarRBGYozCGBH/DkSdhlSlyBjJoi64v8g4aVZoUeNUlJIQiZyaBSm8N1tL4oVyq7svbXgURoo4DBJtqR5goG1Hw2A1z5mPMaDke9k3PQipPV4wEE5AVMX7FAy+k4XLuZaNeG5EVBIO4zzUH+DKr0WFyhqaVc7s/i5ZUTPePULeueavNYEEJ4nCtlfG+XM+HZEqT8uT+HHGh6kOvAd8TggxEEII4BeB7wF/CPyqf8+vAb/34x02xISiY8sCeCIVbQTnLWTrEuO4SSx9PXZZF1y48IhXfuFP2R2fxPdqI5jrLB5XCcMwX7F35QGvfu4bfP5v/QE3en0GcjvG95nsMTbb7OhdCjkgk71oAEJyqTZzVnrKTD/iVD9kyDZbXAZ83C5L+mLMkO2YSAzHD3XzzDMYF/qIk+ZuZCVGYgxOanqhnaKRjiAc6UtaLkschDbAiZVUeuJLdobJqudw7pXLA6QyZFnWUK0KFsseqzqL0GF3DsIlD4Wl8R6CkoamyZiuSr4/qbkrHvlE5zJ6PorcNyJ15UGBu1anSKxiYjBg9d05hex920E6LOJgbFW8f5KUCxCSto2t+MzGJr98KecizzJSF+jJseuSJHfZ4SpbXGYgt9098zJvwYMC12fhlbHgU9vTRODmbEIwHZv9GdsbDmQ2Odjije99FJFpNl64G7tmdXMrbXJQ4nQFtJU0IRGbYl2S7xP+XJZa0YjmzIbwfuPD5AT+RAjxO/8/e28aa0mS3ff9InK7y9vr1V7VXb1Nz0zPcHp2kiI5tIbyUDRhWYRhiF9sWjZoCRRgyYZsEobNL/og2YZhG4JkEaJESaC4WIYoiYJEc5eH5Kychb1M793Vtb5Xb7t7bhH+EBGZkXnve1W9VHV19ztAoe7LmzczcokTZ/mf/wH+GCiAb2DM+38N/LIQ4m/YbT//Zo7fbp5pzmkyzEUpq0pBt91/CO7z9a1TrC0N+dj3f4V0d4WrX/kIGxt7LC+PeOH1B5jkEZkK6piA9dVG+yuUWUS0lfLnLh7wyN7j/PzuU2RigtKmSMYRVTg8fCSN6eiaiYSyU61YB/oGYF58iSIrS2RgJoWjHVNePMHg1etqxTbDcSASloJNNkTMiUQiC9cJd55k0v+946RXOmekCq6P+1xIE4TUlGVg2pdjsgCu5l1Ih683Zm5ZGoSblBpp02Muo7C9v84gi9mIQ/bSJbuqmq7BWk0oRU4hUtOQxAZJfQyBf81Oat6/wHO5aoXkp1udAqhWwmo1VTw3TNmexexyvWo861KSpTTH9VmZ/DjBBT7IJzqbbMQFqZIkniXgiFlzFRh0X5hXE/v06S16KyMOtjfoLU34+GNfY3DtJDdevkhauDjPPBlpM8XoAovNgLb7hctyHeQhe1nIWBxUDUjuVN5SdkBr/bPAz7Y2vwx85q0cF63IlMG1+7Xz4GiXJULrZkagpYiVFuxO+iRhTuexPYpvdLm1s8GDD79GPxmwvH3S9t9rwhhKLRlP+pRlQKcI+cj5yyxFp/kne73qpZmqAw9zaF6YUCQERMwYIG2OXukSRV4pBrPCmYlT6JSS2pXwX/xSuxx2WB3fEWVU9OW6SycUxLLFWUizOaqLrDsYLkBW5qTkDPKu4cIrTf5fu2ApArQJrFbWlXeD/echpaqsgmkRkZUBy5EgSev0oPPBlU5RonZP6uudJ8N0/AMNpKZX/ltTvteITd0KOJqxGuV4Q+4wzPvMMGxGjnTDxSoAWwNSZ3+c1XZar/LEak4nUNWEbDM15x4dWyhLAqnprYxIVsbk104TxDmdR7a4/txDvLZ1mlLPG+FugvvHrupmvAB5m4dTCM20lOymgpkeohfcz6PkPkIM1qL0jL9z6zd4JPgU/9WFpWrFd0GpzK5SUhRz/lLbLLs1WuZrv/BnuHjhKt/1Q3/I9lOPMNpb4cnPfYndFy/y/z310Yo9B2mq6q7trfPJ73qKE9//IvnlPp1vf4DTz3yImwEM9Zbnh9aTLVWjOlWHJKKDEgW5UtWL5fDuzsysglpCVp16lC4JMYShhUrtfJT1Cy5gVg7IxYTv6DXW0z4n5AXiqMdO9hIIh26s70Eku3TFKmO1U6UrJSbzsjNeIrh+plIAPqV1oQJQAVIoVCHJirCxTxvwUipJNyx4ZKlga5ZQFnmVNXExD7+gym/e6pdZu7Jj4w/W+zhIdqHyar+V8CylzhmVWxaxWcO2TSbCnG+HqwyDHmucI9NTBsUN+6yCavInwVJlvUkR0ZPrfFI8waXlgOVoSly5QCYO4mIDrl7FVQmeP7HD+UuvMx0sMRkssX5yh+HuKs/+4/+AWR5VVoMLqhqCXA2yVgxFKZmUIYXNkDnV5gPbzDWaTkU7acDTownD4jplOTx0bi2S+1IJAEyyq9zqnaPQj1evmW/+N/yjQyKhjmJrd9xnZW+VpVurpLOEoggNwi3OOd0fMc4NY477TaZCijRGjzXFuEMQFvzQZp9v7j3El8UtCwoK6hXYW32SwPToM6t82TBd3f41As6rF2g0Q22mLpRuBnsc+86+3EUpRW5ZZNr03y6+4HD8LtglMOwNUsC0iBimXUMgIupW7i4GA6CEpcZScuEKViqTGz+3ukdaRBZv0UcX86Z1u3gITJASMV/p5/bRzfUPRK1Acj2rr8nb7s5lzh9UZCfr+iJD6WdAFofFunKVVTY53ws4kRgciX9/HMeFE2FjU8vJjK7FV8ymHbIsRpUB40mPgd/sxllYQrMo919YPkJf4R4lsxJ25T6lSt9QZgDuYyWgLY7eX+mlDRamVju6yLS0n4GKR8D5+FoLCh3w/I3zPH/jPBdWd+l1p1x/9mE63RlPfvzbvPDMB7hysF4pk6KU3LxxCvElzd7OOkmS8j/9tb/Pv/0nP8ZXn7eRZ1Pxj0aRYlJ3Qkg2xQNkImM7f9HCihPTklwbDoFS5DZqbsxi95Kn5aCaxD4/HtBQFj65yV75OgNxg6L0SVFtxF0a0NAap5mIIcNyq0IjBhWdl2aYR5Sqz0Z3ShLm9OKUvAzJPaXoT3x3nyNZIqUiCkomWYJC8NEvfBGtBEu/9xm+vfchiknaAvo0q/Ei2TUTX03QSOpaOXBNSRxaD5y5btrFTdUBpU4ZFjeo25V7ilIrFEVVKZirCUg4yzqdssM2LzZMf63LilxWIDmrHuSBcJUn16d0LDmIU46u6q9drNaNch554DJCaCb7y2zvbTCcdSm3TZ2AY8JWCPDeXXefXVCwxJC8ZGXQABxBU2n4RXKDTHNdvYB6A9WDTu5bJQBmBUmkMTidbnM3o9Sm0MJ0bvF9KNDM12M7GUx7TPOYQCr60ymF7U68kqQczDoV4nB/vIS6Knnw0VeIuikHz1yiG+X8WP/T/MFoi+viZcNbT9lYxUfs4/rC1d6dedFcLMEPcLXZbtsknG3asrp3QVRZFL5P7bj0pIgoydkV11Cq2XZdiICYgI044/HNLU6s7fH6zTN28huzNpQlhTIBQschoGwWwOTCjek7zeMKEDN+7QxIRRJnnEhKzoYfZE9fM/Rgtq9C4OosLEW4D2Zq+vqmLNpVQ5pna8qkp/qgxhi4+1KlDA2ng3tpHJzY3ZuX5VVmcoxWljKs9ZosB2dYZZOPd9c511WVAnCitLDdhpu/cxmD7e1NHFnrOO1UDUWr31cIQBPHqkqEhakoLZVgVobk1uWtmIRacTGAjlSMi5BnBx2uZJNGb4s3IvexEvAisNJRhdVEjYUWpMqQXSDr1Er1a+3dbOr67VGWQJYghGaWxRS2m3E/mTFIkwpxuDfrMsoSPnzuqwRLU17/7c/SiTK+cH6Hl55b5yZBFSj0I8ozPagmomE9qIlJ/dW+jvgX9prqR+HTkYMhKXX0ZI6Pr1YQNVrOKIaiCoYVOiUrRg03wUkiA9aTIQ9cvMLaQ1e5/jsnybKArAjpxhlJnENmdy5tObcSRFJ76UFJWoR0w5woKNnfOkEQFsRJxnpccFGdZSh3yBhVAU0pwoosxAVFfeou/7lX33lB2Dqe0kQZuhiCgwFXigBqi0uXXC++U1lSClf+XVd4rrLJRX2SD6/mnO5OTTTGdgXykYDt4rRKCQxXTX1Aq1+AexfBK4ATTf5Hx4kxLcLqLmg87oCWGAIdwVMHKVeC1ymKaeP+3ancx0pAsj35E/7HK1M+H/0pPn+mNnOcaTYrg6q8uCIhkdYq0GbyuwIh7SkFAKmNP5wNV7i4uU2vP2F/0mOSx+RKVjjwr/3W99FPZqytHnDq3E36Z3Z4Zu/HYOtJviW+SWEhxS7g51Yd18PenDutJqaZsEUFuw28Ahb3Yq+GBkgyKK6b62203TKf02JIGHQsLHlqfWIJJGhRGjCSkJUZ7lsUWpd0pOT08gHjYZ/8+QftfYX9WZf9WbcBvso9jMC0cPRrtSnqJIoMM/Fg3ONjp6/z4HKPv/HM4zzNlyqrxW/C6lODm/MHhJ553q41qPbzFEAoEpaCUzyhP8RleZMb6gV7HtPwxFV/OjGfQ7SoFYvpnnSe7w4f4WIfLvQy1uKsivqXSpJrYajehaomvMPrO8yJ0rKa/ELoRkVrfY3gmLDaMrHxlDb8HeqIiAISmx7fSRMujyO+rn+L6WyXwqtdeSPylkqJ77YoPWF38i1283wu6u/yp7mSFKUXfGvZd+1iF19Mp6Kwqo5z2t1/QFujZW4OVpGBIoxzgk7KRzf2+MyJgK5Y4TAeO586vEayNVmHFjWLFEhi3SWkRtG5hhzu5dcLmI3bk6QqQPI4DOrjRXQCQa87JctiBgcrBnchFbkKSMuQaR5aenGJ6+js4gO5Cqp/zk1QCGbTDuksQWvB+vKAh89c50LcZyO8VI3RKQPz9xt/Yev7IStwlen4bOolHPrPzzxUz0P4AUHzWykiNuWDXCjP80AfzndzTiazirdSa9fk1txFx0xl6MIL4qAgDswipLkzMpB2ibAfiPX7ELbBQQ7VKmy24No05PpUMS12KcohWs8TzNyJ3MeWANiECAGCRCpS3MQNKsDEIDP+/cnOdA5c5HDYrglHHJSV+eUHu67sbiJ2DU98JBW5qqmgp0VUuSNbr59l8OxjfOE//ef8cKT46l//L3hRvM4BtoWZrQtwEssesejZgFOJT03uzF9He+U3MA0JG8hAMESaAtmqgy/wG076cNhpYdJonWDVTAOL0AtEyIa8yIlEsLw84ubWSfanfdZ7I4TQHMw6SExhlWPAQbnIuCXBxJTJOmsgkqZ45luvPkwvynj84mW6/Slxb8qfvzjg0t4j/PzeDdJyBNSVgm0WJj942ORtsPgG28kplAmJNNDtTE/I9ISXuUlIyFn9MFfl85S2lRsaClI60tRNpOXAU86mhdiPbzzAqU5OP5xUjW/dhHQTMbR1Aq5MfSnOOH9im1makBURO+MlcmsFLMKtgCtld9dUv4dFaTIBWRmQ21iXe7+F3be0Acl+WJApyTAP+WeD59ktXrMKLXqvKgEAycviKr994xJPrBasREUdUcXGf7RgUkREsjbVnFR0ZAuQWU5KZRCI/TilUAFZGVSWh7BosOtbp1jqTjh18hbp1Q2E1Pzlxwd8efsh/tHeoGoi6l5iE9RKK1x/ICKSYIVCpfZllpUf60ptDQQ2JBMZrtNRk75M3nb1dFkK87lururcjThY4tPhQ1zsFQyHS2gESZgzswG+86v7HEx7DFPfErEBKp+/QRs/eK0z4eTqPr3+hNHLj5CWIbd2N1jNhyxrwUMbOwRS852D7+FVdnhVfdM7rinhrcftWouBqx1w9Ri++K6CowLbF1sEFrWJcqnFOshqfudSiRGR6PKo+hAXoi6bnYzlKPeYqusV2SUnHUI1kiXrnSkr3Qn95RFZHjPNYhzHhRMl5vkCHMdltY8NBObK9Ncsj7BaI6uESi14ftDhtTHsl9dIiyFKp29aAcC7QgnAy5Pf4/VslZ9K/qOKHhxqhJXSglEe2fZcZRW5rXgF5GILoRK7/1Iyo1ABoyzB1XdHtpPRs7dO8eSZq5z+0Mtcf+pRsizmR//6r/GRX/kYv/zba6AhY+Kx+JRVK3K3eq1wiqkcNKCuVY9CjSXD6BiWHfIqgGVM++xQ96FSPK00oblHNUJPCklXrPK50ymnuhNuHazRjTOWkhk3h6v0ooyHH7xMeO0Moyxu3CMHDHIKFaAT5JzZ2OHsY6/SObPLS1cucmvS55W9Tc7kpiHnAxevcObUFjuzJ/nyrVNcntagoTblmrkGp/Tqc7smL24yV63PcanDCRO1U8VMItkllB2ErlOTrlOzqf9I6IglPr3W45MbQ1y/xUiWZGVApoJGMC6QilgqIlkSByUn1/bo9yf0V4fs3DrBMO1UJr5bUNriKN4DTH9HlxJ058uVxQVAAxbsbkNsqcZyJfnGfsZX1O+SFvsonaHV9A1jA3x5VygBrVPy4oCZEuRaEInmJbsIaqok+2mHbliQBEXFHehM10CohiYGYym4vO3OeJlIlpxaGlTgmJ1Jn5KARCqu9nZTDgAAIABJREFU7m8w+d3v5dGHXmHzxDWmf7xOb2nCzz9xll986VF+b/qyrRMoccVGbmWWBMwYkTNrTFJ/ohfadN1FGvCQH0xssxi7t8RNfGM5uEBb3VzV3D/zUv5A9HkeX5Z85ORlkrBm51Va8Kf/w98gOjdGXzpP/CtTrn19nawMKNrBLSmRouTS+g4PPnCZkz9+C4ZDOMh49OJllrdO8tT2GfamPaOglSSKCj7/+LM8vHqWh7Z+gH+9NeZ5vo7PAxDJHlKEVaFUbCs3XV1BqWs4tKY0CE2pSMQSYZDgeAVKnZtOPMhKGRsrSgEBiVziUfUhPn+yw7nutPLnAWa2bLewQcBOUHCyP7INQjTnTt+k058w2F1nNOpzY/ukwQF4q/dc7EoLpkVUQY1LIUBKpnmIo7wvbI9CZ/77x3MpymER8soo4Yv7B9wIrhDokLrb0JtXAPBuUQJoIOcgg70s5GTi97PzGFyURAtNpCShS8FRd9dRIkCKYmHpp0KQFiGEsBIUKGmivf0oIw7cCiTZGi3zKCDjnPHWBkGc8+Qnv8nXtk/x7PgUeZCSatPtVooQSb3SlV46D+aRgVorSnLEEYAPlyr0LQBwJq9rLOJcEqM84mCJRC7x2HLAh9fGbNjKtumsQyAL4jgj+lhA8eCnyE8+wvKD/4rTzx+wNVpBFVGFi5fUeevl3pilM7vMPvo5kle+gSgu018bsDTqm4q2ImSYdkgmOUu9CacvXkMIw970/OAkg+wx9rlp0phq5F2HtQ6EROgmeKrxvGyZtUNEBh4XgRnrfIpOINngLOfiLh9cGVTU384ydApAa8FKnNKPMtaXB7aqsmT5xD5hd8b+rQ2msw5bo5W5c7SJQbRdvatmIlqDUmQqoNCiUgCFXdygrg/wgUK3ZhGvjzXPqj+yCiBYiLB8MyJ023F5B8SQikS33S8OT7KZPMZ/e+ZDxLLGl/lGsrvZJoiSVyv/rAyYlUHVgCMOjP9XlHU9fOSZcUtxSi9JefjRl0lWRshQkY+6TAdLvPzaA8zymI984Dni/oywO+OFbz7Bc7dO88WtFS6PC76sv1LjCGxDU5/w061aPmbeiY8daAcNF3H2OdPaJyYpVIoQkm6wzmfFZ/j8Gc2nTl9nc3UfKRWjcZ8rBxt8z5PfZP1jLyEChVgWqMc/SHb6EdTyOUZ/9Yu89OqDvLK/AdQBqkgoPv7gK5y8cIPVz7wKm2voXp/tX1jlxvXTPHPzXLWarSczVjpTHnvoVbI05uBglTjKKMqAX3vugzx7IPm36b+rMAPdYB2ATE0aSk5rU6btMwW7lKO0vQBqfsUmTZsjVNkQF/grZ8+xFmdV3wUhTKNQlwZ0tOBf+MxX6G3uMbh6mt7GAf2Hr/HSb32Wy9unbdcmYyG1eSzchAeYFWHVjrz91DIPROQsWV86QUk3LLg57XJtGvIP9/+AWXnQ2Gea3cA0GbnTOZx9XWv9qfbWd4Ul4CQvdxmrHYR9UK4fQVtcVDdTAYGtNpRQaVo/7+3EB3+4NFhWhLasOCTsZGSTLuPhEsNZh0kec2t7k+XZiJVTO/R7E84tDTg/6iNFyHTwSW7IHba5XL2wyg+C2cBWu2+Ar92dMvAVRLunoQn+1VwD5lokm/GjLOlVTqsNPrIe8MTGTc6fukmnN+Pq1bMMZj2mecjuzU2i72R0z+wQpFNk+AJR3CFdPlelrmKp6qi2NOjANIsNGclqD3YO0K8M2Nu5xGDaM1YXdmIpySSLubV9gn5vwulz19FlQJbGfObkNpvJGvmNH+BaPuam3GLMXtVdyNyPuiN0o8egRUqG0jQkadOoVTEHTHHW4/pJHoj6rCejxvW4VuGOsdcFAAc7axRpzP7uGlkao4uAwWiJaRHNBfygjjOZdCrV++f3DmxUw+omX4bfgDcOSgol2UkTnhtEXJkUjPNtlM4JZbcqOjPX+dYX8XeVEtBW20e2gksKzczmsn1xrZ7GRYhjYY2EohMVzMqQTJsKLdcRpmImqkw5Zfxh1eHpVx8iEIrV7pQ0j5gWEXuzLqmSfPvqA5ztD/lQZ0a3N+WsuMUHJkuc70Z8bF3wh1vn+P1swlBsUZBWL6YjwdAolMc8XJcV2356qkBKU3/g19UjTEuzUHQxjVHrMmJzHPiEeJwLfckHV1I+snmFD3/oOZbObaPykD946qOM84hSC/7wlUfpXH6IH/r0V+msjFCvBiwNv0F3sMu1g7NkZchaMiUrDW6gF2VEgWJ3uMKp4ibpk59D/MIf8tJXP8pz188zK0ND8oIGizvIs4Bb1y7yyYuvcukHv87gKw+Q5yF/6t/7Ip8Z9vneFx/id155hC/deog/KEZVo9JqUtvqQNcMxJB+FLa+IDToP8+ykqJmYg5lh0h0+ZGTXZ5Y262qIAtqX7xKv0W1Qv76K49Wn6XQBC/W8aS2O+kWHb+3Y26PDQ4cRIMmrH7mtTsbSkUkNCtxymujZb4ziPg306e5MfsTnJNQqow4WppjY3or8q5yBwDCYIWPxz/KA9Eyj6/A6W7Gss2dOu3qR3ZDqzAiq+HbFOMuJWRWvLKlDOr0YljRQJt0pHvAnaBgtTPj3MYt+ssjXnrtQQ7SDtuzLjtpzE4acGMq2MuaLoKTuqinbnAxLWpSESlCjxZbNbAAjlWnUCmR7LISnuXh8hEudBIu9eFEUvD42h4rnSmr/RFFEZIVES/unKxefHP/NeeXD+hGBi780Ee+Q/8zO2TP9CgmpvFoPjaukKvCzPKIXnfK+plbXH3lItf3TnCQdgy2XigKm/IKhGnFNS4iTnUnPLi5xd5omUIFPPHY8wAMD1aZzhKmaYfnbp3i5izhmf2Ia7OcK2KbfbbI9KQKEvrtwZJgCaVLwx1oA6mx7NMPNvmofoLNJORUR/DE6oT1ZFYFikvrj7sMkGv24b8XgGW4bpYL+6a7cw3cL93kb6/65vnNm/2BMO3gI6nYSRP2s5DnB4LtzLAz3dQvk5YjVsNzLLPBSXWC5+VT7KWvUBR7b9ASeA+4AwBFOeBrs/+bG8H3sjT5EOuxJIgNWSNUQXMAq/ExPHlYcy9QFX9AroIKqOGCM5JmbCCURlOnZUgoDKouFAopTI53XETsHHRY741YXT+gn8wotclSXOxNudRXXJ/2uDYNeXpgOu24/nx+0YsQAbHsGWIScYBSRaO2vrIO9Nj+Jqzq9Lthj45Y4VR5jseXOnxoNaMXmDbt670xpZLcOljj+niZqbWOTJFQbRa/sr9BJBWryYwT10/SvbZN/PCAODbn03sFvZurTLfXycYdppMe43GPV77xMcZ5zKwIcX3zDHkmGPJQXa221yd9rl9+CDAprwt7a8RJRhjlXDx/g2TjgAevnGb/1gYbLz/G88Mu8cFpLquQoRwwFbYAixLX2j0UCUrYIiFBFVfY0Gf44ErE+V7OpaUxEl2lOR1vn5uQjinYF5/Yo/QUGsorArImv2+JVsS3+ICgxbn/dhBxkAfcmEr+qHiOETvMioFJZwYrbOpzbOglLnY7vJwmb0IBHC7vOksAQCA41f8U3xd8guVI0AngZKJYiUpOdWembZltXQY1MaS0VoEQmtAGD+sVngZOXGIe4nKUEQaqekl8jIFDgKVlSGTZZC6s7hKGJUoJrh+sc2W0XJmD0yJkNwu5PpV8dbzLFfF81dbcNRoBg4xzAUF7f6i5CExjkUj2eEh8jHNihY9vSPqhIpGaFcuECxAKVbVly8qAaRFWL2wgdN2m3KoYIcxvwkARy4InHnmJTn/KaG+FLIvJ8oiZBRFdvHiV0WCZZ648UN0bn39/rg2cVKRlSFrWfvKJzrRi4Vnrjul1p6ysDog7Kb31AaoMKNOIyy9dYme0wuXhCuMiZD8L2UkDxoUBFQUCliNNR5ouPLHURFKzFJnmIM51BCreCGXfi9CzDN0q7zgTzf02VX1jW216GM7eD1I7d8BJxcgM1XvXDw1V2V4a8/RBxAuTCbtyn6mYMNUHnNYPclGuM1MlM13wgnzaoktLRtlVCvvevDF5j1gCYIIh03KPa2pMr4jpyIBeYFfTLK4eVE0/VjOyuHys1gKkIvB0oAS00EgtKnx2oSVCGQuijeZy9OYOZEIJ0zwmVqZ9VxyULEV5NTl61uyDmEuzNZLiCVKZk5MxFRNmjMiYGP+WAC1KItkzpj6niIgJdUhETKRDHo6W2ewIzndndIKyAWpxlFTOXM1t+guoLAHn2wZ4DVwRTPOQMRG3tjfpDFKmMxN8E0KTxDlSlkzHXSbTrjGXpXGbRAuE5awA81ujYEpZB3NnZUigpEnVzXpkZUgYlGgl6a6OkEGJ7Besr+9XxUnTPGKQJWzECeMiqK6nEygLFlNVUc8cjNx7foFntZQtX73Z0aruuARNpnq/me3tRApdJS1LLdhOE6aF5OYsZGtWsC/MpI51TGZJYDKlmOmCKSmz8oBMjSnKMaUHNns75F1pCdS/C4jDTZai0/xI53uJJAxzzSPLcLGX0rd99MZFWOV/I8tP4LIBriBDYvx7NzEyi+V2mYjlKKvOu4j11a0grgKxExasJgZe+srOSca2+s791lFEB0JzkMW8NIp4YZhzWdxkn5sGBETEaf0AF+Uq37UuOJkUVXrL+bGO5NInmminrBwoxb3Ia3FaKYF2paALkgJM8xAh4ER3wnpvxObmDicuXUMEJb//W59jnBt+Rh+E5Wo1HOFGW3H6BTZuH/+8sTSB2kC6iawMbbddpd1zK6xyK22/BFP9GVgYrplu/iqfWzIPrQ0fQOKsoAWmulu5nQWVW+vQSUXCat8b9zwb1+llroySMim/vTRhPw/54hZsqwlbwXVW1DpLusfjvSW0hq+k19jlOgf5FarmqBiOxrIcvgU34D1kCThx2YJZOeBbs1vEOkIgWJ6t0A8jbs5iIqk5mWTGJPcCNKUWCC3QFqKpMSuT9FYQh9eWGJMfHA5hvguSEzeBilKyP+syTDsUWlpadF2lkaJqfw1xxoN9WIlCPpBfYFZeqOoiliNYDhWnOjndsCC25nO7x4Lfc7HOL9SsQI6w1fUUMPePqulFNVntGN31A4yymPMnppz+6AvceuZhdnZOMKv4BkFK00nHh2tXrDkekAuoILPVtWOtMvu5sBwRynaaKpAVW1STcNak9oqK9syl+pr7uMloysotF4J1ffzGtf4Y6+i+6Sjsp/caJB/u/fEWAifa3j+HO5gUIdenCVenIfuZptAFK6LDUnmJITN25T5PTUoykbMtLpOWA0MG49U/6DfQVeiNyLtaCQBoXZCXY75T/j5h0GE5PMPJ6cdYCkN2UuiHcKZTP6jKJ/Q1tzbUHzlUmQDnu7lsQmqtAmMxmE6+TlwEuWKhtS/RQR4zLULWk7QCKJlqMRrw5Ugq+mHBpaXm5Pb7Khba1c/7ABrr3rjVVDeZa9w4YguBdVJTWjsSS1EdX3iWoQNUjbKYJEmJnoy4/jun+ebNcyw72LEtnw28YzsF4KfS3BgDqSnVYs48Q+BhJhfScg+CYUT2Jrcfb5hZ5ezu/dw12ixAEpR10JKmAvAVgrtvLuCXKVlZjpVrQ31sDQhvcXELitLCnlMTy5KdMuGlUciL4xl7YsAZ1lkOQjYSybPjnH1ucrn8BoWaLkYCLo4tvi3yrnYHwAQJEZIwWDUEEjKmF56gL9Y5UZ5iUyzx2RMBqRLMSsHjKzNWorzCbPsIL/ewzZjmfUpXTeabgKFUJJZ/3xXYQA0X9UuWUyUJrZJxq+y08F9iIy5a7V7ICpm2IMXpWyRuYrRr2t3k8Cen1jR8YefGRC1LwwVC46CkG+bszbpkZdBQYs4y8rvmuPHU97OevP6429aDbyW4a3F8EU4p+cfx729un6nv3rXF5/WrMAM2+u8z+vi0YM46bLsOpX991AuIswD20oRBHvD0gWGnjqRgkCtGZcHz8iUyPaEkZ1rsk5fjqo6kKY5Va8wb6S+4WN6D7gDYqjOtqxxxqTKG2XXGYhsVKyjPcHW6Tqkg15qLRUASmECeponUAu8lpQ4CVawuNmiYVy+g+d+hF82+ZmUWzq/VqkqRuf1DWmYlNaCkAYFGG2CQ97KVWjYm9SKZ68KMJcU4ZH+F14rNU4L1hIRhlnBj0qcXFOZ6PSulehYOhefStYecrx7b/HjrY4mq1NZZKYFWhyo6/3pdnKf6W+jaVdCCUptzuuIolxXKPKVYeO8BNHP9Tvxqv4oCTwu0tTbHRcCwCLgyS5HAchCRKkWuFcNyi1m5T1HaFKxWiEMIahxCkrk35O2Rd70l4CQMVsCrpXdiqvm6bEQPcq68wAe7y6zEkJZwsqP4wMq4go26Vdn8zpqv3qpRrQje8Z3P5178OCirunO3yszxH9qiJletZtqm6bkVu/EbG6ga54bkpNcw7w9/hn6Kru7wXANhfLPcKS7Xxs0wB9Ul1UlYWOspYGbvVdiyTAotq7Lubnh0jbtTJq7O3rkD7lrdyq68azD3tAZ4dez97niNUtpUcsrb3wVInZ/f4A70VvuK3p4afOZ+575PpMlIuDqEWRlyedJhJw3YnmlWIkGh4WuTW+zKm+xkL1dAr6I0fID+yi9lbMdf4BOKmeaqNbL0zcubtASEEP8A+FFgS2v9EbttA/gV4BLwKvCfaK33bE/C/wP4EWAC/ITW+o/fhtHfVpTKrCYNW8U4OYWCYbnFtQCYXmBl1mE5NFRde2lCaifDSlRU5r576BXVNPOrtxQ2bONF3o2prBqoMeevViufroucfCWDF633TelF4sx3AGfE+pH26vr9e9Tyq30FUFgLQ+rSdHW2vQT8vHk3yphkMSWymvztiSdtJNwF8trf+7/xJ7wu57MJgdDgAbcMpYqsjlfa+1AIjSL0pg3VPu5etie9u0O+ie+7g7S+l4C0mSUfTHaQR2Slw3cIdtOAG1PNlWxMlIbkFGzLK0zKvZa/P1841hYhDEyYt6la8DC5E3fgF4C/Dfxjb9tPA7+ttf6bQoiftn//98CfBR6z/z4L/F37/10XpWeWvGMVaFJXaa2Y5LcY5zfZEt8hCZb5QP5pSr1MP4yYluZBf3w9q15gZyKW9mXzU3sNsRFiJ6UWNXrRmqZxUIKgojYTAjqyBuoo7Zc7C0Jcw1Q5ZxU4hqRMB1U+XNoIv+9nS+9ctZnuVtTaInC+eKnqKLoQGi2oUm1g/P4kzEmLCKG0pXr37r/bD00/Kk3kXsvqWoDKKnJZF5faLFsTv+J+EMY+9H34Uusa86AN2SxAJvTcxPZBOpVF1LLk3BR0vIG+ZeCuywHHHLy3FxQUWjItQm7NIq5Nbf2GhlRpLmcjXhDfYlLsmGYgWlV3yF/5awVwuCLQthv23ZQ7cgeEEJeAX/csgeeAH9RaXxdCnAV+T2v9uBDi79nPv9Te7+jjv3V3oD5WgBRdpO1Rb6T2pQz6LqIfnaQjVlhlkxW1zJJI+O4TZjW5MZNEEkIBJ5KSXqBYTzJCO+lckUjbT/RXkgYLr2eKu/xyZNlstBfRbwbjvKh/hciTVVFUHJQV4m2R++DG4AJn7d55hfaDaE1yVYdB8MUF/lzhUWSxAZELinrX69+DdvqtDcRpb2sHDM1+9fW4eMC0CEmVSeG1G9L6E9spgbaV4PZrdBGifpY+4MhlivazmHERsJs5HAK8NhZcm6Vck9vMxJhcp0zUnvX3p2hdWAu1qQSaFkDNemyuM0frAqWnb7MCeHsDg6e9iX0DOG0/nwde9/a7YrfNKQEhxE8CP/kmz3+IKBsgjNC6fWnmVdBaUeqUYXaNETfZla+wHj3Ihj7DjekZNPD6JKcXBPRDgRQBZSRYjiRIRf04mz4k2BSR5z44ccAgoEIjIs2qaVZf8zs3Eatqs9Yk8VOCPuJt4Z2wE6ZUJujnx9GqCLsdjxJNRaO0qAKc7rypMjBsN+l0UBIyryzMeO21tvLxLuZQp+OMOMukraiq+3cH+TFfoSwKny1SznPHaP+mdZxJKdnPA16z7nknENxIM7bEHtfKZ8jKcaMjsNaFy/FU2/xS8KOvp6DmDby7pOBvOTugtdbidmHgxb/7OUwrc2sJvB3i/MWUUuUEctlqXLe9sP8rm3aTqDJnR73EnnyNq+M1enKdc/pBVsQKJxLBC8OSmSqRhMRS0gskl5ZgPS4535sioconO8BKmyOuMoe96LPDokd2FXcQ57Y4i8Of7CYoVS5UAM4VMBNRN6re/IkSBBBTzq28/qqbK2PxDDKLDJSqCtS5Y45EVKXFHGTZT28uCnZWcOW2/02dDq2vx8BtTfpVU5SipkK3Kdfa2KaK5xgyMYegtMfCcPYrbTI8Ci/mYa2rRCpiWZKpgHER8sKwQ67Mc5wUMC40XymfZ6z3yNIRuZqibFmzb/abu99keHKjm7cCHCIwB1SVDrzbk9/Jm1UCN4UQZz13YMtuvwpc9Pa7YLfdU9GY1BzQ4KXzH0LdzcfcfFXmTLWiDFJuhT2SPKIcd5ipEqU1gTD+6Lgo2U1DShUAXYsCFPRDAwZqN0uF2hRu58i1jSf4frF7IatodGsCu89VIdAhloATpxSU714s0Nltk9vlzsvWmJ3L4x+jotCyx3HWhBMHZfaBTf4x20FQh2Q0f4gqi+HGVJXr0gz6OSsAbB1I+xq9z6FV1A4EpLRgXISMqe/1rBRsz2BcKCZlyUyXTMk44AaZGlGUs2riHjVh60l/hwG+u4QMPEzerBL4l8B/BvxN+/+/8Lb/FSHEL2MCgge3iwfcLdFohJDVhHcluTUgw0feGQuhVDBVKVfyXa6JCFmGfCb4AptRzIlEMsw1O1nB5YntMjQQBELSCSRPrsODSWpWeM8dEFDhCkJvIkVCVwFFTa0o/EfvwEXOPzU7mEkwyGOWopyExSnFtvhgHhe5ryLluk5lVpTbqq6R73g07oFVdHW6sXZTJmUA1D56uy7DSeUueErAxTa8B1iNK1M1L7/bw9WC+GNwaidaoFwa7MHW8kpsOXWmTOr19UnMTirYSRUnEhPse3k24npwgxvF0yhVLJj0zfzLYio4OffdIqVh3lXuqQKAOwgMCiF+CfhBYBO4Cfws8GvArwIPAK9hUoS7NkX4t4EfxqQI/3Ot9dduO4i3MTDoixQdhAiRMm5ZBH6Udj5W7PP2byaP0WGJru7R0z2WSDidxHQCwTDXZEoxUSUn44iVSHCmq+iHJSeSjLa4MlN/hV2UCmy4ErSKlLx9/Ny/40FwxTy+Ge7nzR0C0QUf2xaHsmP0J3JwB0rGmdj+sXzgkf/Kt6vv2qhHqDEBUE/4duTetwT84/l1/K5iMrBjUBgI+G4WMcgl40IQS3duuDJRvJ6PSEVGKmbsc5OZGjDJtxvupLmOJtXb4bLou6YDYyzTwi5Sb76HwNHyJgODWusfP+Srzy/YVwM/9cYHd3fEpQ0lMW3G2voB1j7Z/O9zbqUvVA/4RPwID6vH6IeC5cj0hJ8pGOmUQTojTkOSoNfwQ/01Qzr/tVVrfpi4131RXlvaLIX/KsVBSSQkgVZV9sFnCC7s6u5cCX9izcGPMT61ErpaWaGOmvvxDYd3KMvAoCrtMaoJqwVtAK+7NpctcRPejc2lEduy6H75Pj/UCqRyZagzNIWSjJRkL5W8OtZs5TNWZEQnkKzHkr08Zyu4znb+IkXpGsqUNlCncPRv5v7diXlf36fDgoI1tXx6z60AeA8hBo88PgIhu0TBMq6P3+J8rf+bYG57KLtEslsxAClKYtFjXZ3knFjjRBxypmui37uZYCctuVVO6RARColCsxqGbCaSD69mnEjSqkzVFTYtimK301jtiVNfhxdR90qD2xaEXx3ZXq0dyYqf+agJNOuJnVffCVvdqOYVSQNH0Ixp+FPBv47Uqwpsix+X8K0AZ624DlSuj2CpJZMyYJRHvDaOGBVYDL8pLHt9rNnOU14OXmKmh2ityJQN9nkVfO08f9Okv72ZX7VWX0CdrrVCqRlKH04z//bJe7R24PaizASwHYNdY887rcqq2n5L0xijLFNm5b7ZJkI6wSqRSBipPlEugIBCabaylCEzBvKAdbVOpENSclSeABE3ZhF5BZetV1O3ivk5a18adezeBDZjrS+qbP2mkbPXAqVkAwXp9mvcOTvJnHtQ+C575UpApwgbSsC/tT782jE6OXFmvuOHVNQBOxdkdeNoj09ZU99NOY3lgPAsFAHMlGRWCnZS2MkKUlWSyIC1KGA3zxiICcNyi6wcUViyDn9St1fuhQvGbUA/R6cFHUHsvcsGtOV9YQn4EoUbBDKpOfyERNn+d/PxAvMABQFShtV37vvAYb0rn07hqMBCmXAifJjz5Tke7fZIAthJNYOi4CYH5CJDasljwSnO9STftTatVlS38s4sKKmNqBPWRK/746m5ldeBhHxMgQMIbc06DAvJThowKaBQWDdBM8hLUq1IdUlp74c0qAEiJF0ZEAhBIIx1UWjTUKPUmsJ1ExKSSEoiIYikIJKwHgt6oWYlMjGTQGimRcCklOxlklkJuYKHlkrW4oKLfZOMd4And32upqAuCTaB0nERsjWL2E4FN6Zm3Ge6AZGEWMI392dck1tcKf6EUmVzE33RszeyyAJoSx1LckjVRZZm3UkqsJDglMIuKPdG3reWQFOUfQGciQa1pq6q4hdYCe6hNjoILYglGDBSRqky9sVViiBlOj1HRMhETEnFjIkY0tPLJCyR2M6bu1mETs2Ycmt2O9qsRNYU6w5OPFWSYREwKwWBcCZ7PfBcm2yHzx84K00Fnct7g+HnSwKQ9vtxIciVYiSmTMWEkoJCGCUZ6pCAEKlldYDSefuCar+ImFCZfbtFQkTIMI/ph5KNJGAlknQCdz8hljaFCYwLSaEiCjXf3UdpM04zbnsvtCArJbtZwM2p4GZacFXvkomMq6M+ESEJEVfkNQZsUaqszum/gfVvsQKY3+ab/H6AWQiJkJJOsEYoE3I1JS0H91iL8HckAAAYVUlEQVQJLJb3nxLQU4SOkEHT8mg8ZJeNw+t/5/W7d9t8a6K9qiiVMUxfZySuclM8Xb0IzkroBR+kq2O6oUBpuDULGBaCSQGzUiMFrMUmWHUyKegENaBoqiS7WcjlsWQvc0ATTe5cF2CiCgo0OQUm66wZyiElBWtqnWU6nIpjVmPBUmiURCAEw1warkQxYZ+b5HpKWowaK6WzfNzfUoY4/n93H1zr706wSix67KkTrKRLDIsuJ5OAfhjQD811JtKY+rkSjArBvhLcmDkyD/PPfY6EJpQQWTcq14YnYmumuTYzq/0N9TxpOaRUKVJExoorrYJuIfqaz35RkLg90ZsQ9GZbuaa1UE1+2xEplAmrwRl6epl9uYXWJfciEnA7ed+5A+Z8Easd01wiV9PKfPMn8lz7rwV+XWXeLXgZgIVmoXMtItklkj16cp2AkEBH5CKl0CmJWKKvlnhYnuREEnCyo1kKDdtRoQVpKRkWgmEOmdJMitokj6UkkYJOIOxkgUzBpNDcSFMGTNkNdqrOR8bUNyZsScFMDylUSqnTqseBc5eEkDSaoFCvfP4L74sUIYKAQISEskMoEmLRIyAk0oa52DWUUx4E2e/M5Jp4+ecUBEgtyYWhb8+YkKtpNXbXkKT9TBbLopKio6V9nW5sQkiW47NscI6L+gxLQUggYFSU7OoJqUhJxYxX0i9RlpN7FBB0cuwOeGIeeiKXWA3OVC+hY3rJyhGlLnDtvcACP8ThL9PCnLFN9NfRZWmtBEjVkEyMmYq96hhuhe2Ea6RynbVyBdIErQOGYUAgDOWWm9S5wqAYy9KEl7TCVDdIQikQGhIPDqyAQpRM9B6FTknLEa7rkYtpKJ0v9G3N9XlKrVHgY6wD6SnyihzTugs5IMSwtVI6gFSddmtaVKqxr1FCRaWUzG/Lxv11Y/elCQ5rpuz8sd4Z6q+13Sq/UHYrRbfCKdbVOue6MWsx9EPNfhbRnS7zeqHYk9sU5RCt0yPPea/kfakEtC7Znz7LR3s/xl86v1L1JfjGXo+bU83zepuB3DP48HJkVkOXWqR+OasJ4tDyniLwMxGul17dhtxWlmnQdvXzX+y0HJCVI74hriJUQFCE1TkXtRx3vekExpSnBJ03lVU1fmvlAHNBUX8y+cf3xa90s0du3ds6yNYeox9UdeesvnPWgMob52n2XZQt1B4GDGZNfj+t55/TIUXbx/PH2HbpFimT6riybn3m2sQ9Lj/LuuhzOok4yBX7pKxEcKZT8LETu0zziO1Zl797Jefq6IsWD/DOKwB4nyoB9wLdFK/xpVuf5qElxUZccKFXsJlIzmenmJanmJUGN54qzX6ZMhQTDuQOI7VDriYUorYUXBxBCNlYFRuuQmUZ1CvIUeARbJGTM2/dd+1VyXTeba6sR5m/c25Ny23RWh2aQm3muuePc+g5j/rOc58O26Z1idB+m3K58BhHHbepUOoJ3qb1cuXm7jOYe5sEKyRyiVP6Ij3dZVV06AYBsRSc7wl6gWYlKtlOA27NepQabs5Cvry1ybCQbM/guvjmodfwTsn7MiZQiwI039f9i3xitcv3n9pjozthzaantBYcTPoMsg4vHqxwdRrx0rDkZbXNgdxhqg4orO/clnYgzX12/7uUowtU+S9l5fPaFccXtyI2z1Uf9yh3pdq/BXtdVE9xtNkMTfO6/tyeYIvJM+fNa3d+YWMIzd82x+WDbyoLw6fmXqBM/PO2gWAN681+H8i4EeQMRMJZ+Rgn1DqfWutyvpvz4fVdlpKUju1JoW0689ZomavjJb622+HatOCP+QaD/BqzzNXSvVOTf3FM4FgJAMudR1kNz3OpfJglGbEWGS6BXgibllRkLa7N8Mwy2xRKkirBuDC57mEOr09nDJhyQ77eCFYtUgrubzgcUNKegIt8alVBWss5M9oplDblmvl9MLdtkSyaMG0f+zCf+ShkZh13aHLqud8JES6MtSwy0xupOS9IdxjYRwhJIBICERIHS3TEMpvqLGfECieTkPM9TT9ULIdlhc500ONYKnJtiE120pD9XLI908xKV0tSMCZjO9hiovfYnb1kUYGTxj2793IcGFwg5mEMZy8y5GVuBN8ikF36+iQrxRk2yhN8sLPGZhLQD0vW4pT1zpQkLCqwSlaGTLKYnVmXvSymE3TYSWNkLpmKKVM5IQumlDpHUVKSo3RBoWZVrAGYizfA7SLazYCaL4tW7sVuRL16+r896rzNmovm+drjPgyOvdg1aMYiHBmHIKpck8N89PbxnX8f2G7O0n6WtrGrEJKAiFj0iHTCkl5hWfW4mHQ534OLvYyHlgcsJzOWO1Pb5CQgK0IKFbAz7TFLE14dx7wyUrxe7PGaeJZZOWhRh7dJQu8P878t73NLYF6ECBAiIZR9wqBDJHtVjjfCpLiW9Cp91ed8uMxaJDnf05ztZqwnKctRRhyYPoRgwDqm83HAIEs4yGK204jXxwF7meJ6PmEsZuzKm2R6QqFTcjUxSkPXwUKXalu0YrcRi+2Vu70at4N3/r6H3xe5cAL7QTTf2vHdAl+xNHAWXoCvFmuCy2Th+BzWwl2ztNaCwyVIEdGT63RY4rQ6xUbQZTMOeKCvWE8Kzvcm9KOM5WRmW7GZ978oA6Z5zCiLGeUxN6YdBnnAyyPYzQqus8dQDswzwqQl3bNSumjgD5TKUDoFXbwjBUGHy7ElcEdiON1SA6XVKbkYVy+4FBFCSAbBCpHsMlYPsDZdYVj02M0S1uKY9bggsR1xO0HZYABKZMlKlFd8AuMi4Ey+xKxcYpyfINeaXGlSW+9goLuKAvu/KJmIEQUFmZiS6xklucE66JJSF/P+sFjsUszFFVpR9PZ3Tty98NGV1WRvbfNN7xqV2Y7ANwNwgYwJREQke0gRIgmIRY+QhFjHJLpDR3cIkUgkMQFSCFM9KYTlLpDEUlTpudWoYLNj2rU7XoNRmpCpkFKZLkZZGTDII4ZFwCCX3JiaRiGvlfvsy132y6uk2aBCHLoxL4IZG/fm3hKDvBU5tgRuI4FcqnxTmAesQBNEY36TkAQrnJaPcE6d4lwnYSUSnO0aXPzJZMZGd0o3yuhEGWFQ0knMSlJ1v1WSWZowyRLGWcLWpM/Q4uP3M8mtmWY3LxjplH05IGXGjBEFKZmaVPl/X8pDgD7mvGpu/0VQ6UXgqHlff3ElXTsuIkVYBd3M/xFdsUKkE9bVBn1iejLkZBKyFMH5XsFqVLCRpKwlU3pxRi9OCcOSOMpwTV/KMjAmfBFSlAF5GXIw7TLKE14ZLrGdBrw6gsv5kO1gi/3yKrmakttGIO3n2RYXw5iPlai32DD0bstxYPBNiRCBKZ2RnYUmcTuS7aMCO8EaXbnKMhskusOS7pEQ0g8ieoG0xS2iKm3tBZo4UCyHhvGmHxZVMMqvty9t5Z37nCufYcexA9WEJTVdulEyrhrQMBphm3mYvwsLqVfa7b84ZyBacF5Totzc5moWpDA1CvVvHesQFR+jXxNg6v9rRiXHAOQ3agFzfanlG5yUQVWQZOofBNMSstJ0npoUmlxrDsqMkZiyJ7crHEimxg2EYdu1MWP2J7xz07xgq5pxdwlB3g45dgfelOjKb+20tltwUMsiAKMQijJnVE4ZcZVbwlUg1i+P29d1SOqGa6xxhhW1yhm5zHIYc7aXsBaXrMemJXknKFnvTA01VpgTBSWhdTekVISh6cLjT/wGZ38ZopSgcIzBCLIirAJfrq23o/XOyrrlec0YbP53dGCBrLs0J2GB61sQSFX9E2jCwETZpdRIO+aq36BPFaYClDKlzmZM0oxRS7IiIFMh4zxilEdMyoDr04hhLtiaKW7lGXtiwJ7cZqaHjPNtr2CoBkj5lshRAc5Fz/vw710G6P4M/h0lx5bAHYojJgllH6jdgTbxhBCh9/cimV9XHTjFBblCkSBEQCS7SAICIkKRIJEEOrLesCTSsanU0wldIk7GSbXC9kJBKGAp1MQBdG01YiA0PVvK67ovu3LjeuU2n13HX9fia5FygZqctE1fZq7QxPlnZUBaBkyKgEwZKrJxIUlLs2LPSshKzUFRkmvFTBdkFKQio6RACWUqGikoRY6ymRQ/LuLqBtrIQR9rIMQ89qJ6xgsmeTN9Ob9/WU54pxiB3rgcWwJvSUzj07wClDS+816y26X1muAa9xsMYQkpoJhnJ2z/3kggE6QwxUhxsMSJ4oJRD1rSTbskRPRFRGKp0uPAuB7dILIlxJpAOFYeY5aHlpBUAql1M3qBMmZ+tXLX5nuh6zhG7rkdQMOKmJTCliqboqdcwdBO+InOmFpOv6HYI2dGqkemiKlMqyyJD2xa5I/PYwL817u54t/+Oc3/ri0VJ+C71AJwcmwJvAmRsocUNaLsja0sTXz+0eLv10ToOTBN+zxti8Qh6haNp42VP0z8QJ9qld+WKp3L3zuyFajRjK6+Yg6e3JJFKcjDRVWmvahcLKdk5cLncximwf/uqFSo+VzYNOD9UAj8RuTYEnjbROu8Avk4NuPm94fj/M0q5G85qoLtzqraFuPmjRnrV0Iai/V2Cujo72vf2vRiWghCUvOw3UUK507gxPWYfJGNbfPZiDtbledqKLw4zyLFqVRmFdkifMO7V46VwJsQrXM0DvOfEIp+Zaq6qkEXafZXfj9vXn93mMxPRvOSMre93h97ntA7L9X2phtyOMy3eQ63b73qmuv2J8micR5mBdVSm+W3t4yax1StSdpUCm1ZdJ5FAKpFY/T/Nr0B89uO9d0mx0rgTYt98XRBVuwiRWJWRwdxXWBazhFxeBNoXiH4L1pTibRX4MU+sjfGhS9++xxuu+u2eIiLU/3mcLjyImnWHLTHUytHrdvk5L7irK+/Xfm3kOrriHt6VEFVLa6wqOBedAd+p+S9pdLuuUg02lgGDiVmi3lgAZLubZBFKy+4F/n2j3MRL+KdSXMC+f8OPdcCNF1ze1PZOIjz4mrFxXK4+7D4d3ceEHQkK4V9vvdX+e/bKceWwNskSs9QtlmFEAFRsAF4VXyHBJsqwlMxb6LWv6lXyLqgJvT2pfE7d2xfKuXhgEJVALG9ss9bETXa77CCIBedp3WMOkLfXOHbk31R4K7tjvjXVh9pEXBnfv9FJv6iMdbjKdX4Pbvyt+W2ak0I8Q+EEFtCiKe8bf+LEOI7QohvCyH+uRBizfvuZ4QQLwohnhNCfOFuDfz+FQVaUaoppZraYFJ5aNR7ITutBzqqxUziRZ2U3HdvedyHFPO0x+BvX0yQ8cbGdCdWxZuVozID/j5m1c/svwLmYirvXbmTu/4LmN6Cvvwm8BGt9XcBzwM/AyCE+DDwF4An7G/+jnDcWu8bMS5CqUbVP0dzvjiKrxrWQhueaopuFpN4+qhD//vmeWRjdb29KV9PXgNgihoxDL92oFGQtOBfdX0tV2W+4vEopXFU7MEdd/G1LLKsmpaGH5MpUGpS/buf6L/uttxJL8J/J4S41Nr2/3p/fgn4j+3nPwf8sjYMiq8IIV4EPgP80dsy2nepKD2FqhlINIcxcNV+hwWyFqWsGspioXVA8/iNYy6qBpiXoyjBGvstdD0WiRv77c99VObADyg2Fd/h6dZFBVJKZ+hDuBLfT/J2qLq/CPwb+/k88Lr33RW7bU6EED8phPiaEOJrlaP6HhXD0ptX/+Z9Yl8WWwvtz+1+CEeL4s7M/fnz+qmzRQG7Ouh5Z5PIX/nv3AU4GmA0f/7D7ms7EJl7/94f/v8ieUuBQSHE/wAUwC++0d9qrX8O+DlzHPne1gKeaJ1Tlq5rjySQy43VfRGA05myUFsFVb5+YWrwsJX+9um9RWWyzfH7Ef6qY0DjvM7891NutQvRPP78ar8onXm7egz/mtr3oHl/jFXm4jHvm9fuSHnTSkAI8RPAjwKf1zX2+Cpw0dvtgt12LJ64l094k8LI4jy2+c3RK9U8OMiXxbGE+ne3O64vt/v7djDkw3ARb05ufxwP5KTNv+PJ35Q3pQSEED8M/HfA57Su2BMB/iXwT4UQ/xtwDngM+MpbHuV7VFwA0YkUHQtBrgNzftHMwmPcdvU3shBHb7szHxZI9I99uN9fxxgcUlFTIKqUobMOmmnP20GGfaXWruRrxgsOu2bz971r+/3uldsqASHELwE/CGwKIa4AP4vJBiTAbwoT5fmS1vovaa2fFkL8KvAMxk34Kf1+drbeoGhySmX5CYREyHrCtCfI4ZVwbeafw0FEh6YtF6zkPjx3cQCzNvul7YLUdG+sUjiii9M8lHfx+OcDnU1xaVl33mM5Wo6rCO9LMWnBIFittjTBQ0YWm+qL/W2X6jvKVF9U59Dedhis143xsONq7Zq3Hl7PX1+DX/fvczT45z4MEOQovt7d5b13R46rCN9F4rAGQ/OnVmBX4sCSmjT2nmsNZsRfiQ8zwY+uWaiPUSuhwydwO43ZxAKEjW2L9jfb2uM6GglYlEMcQKs6xvsox/92yLESuI+l4Unp0vrZxfwkksmhq+vRE+6Nm8qHuyFvz++OrrCs4wvV3/cdrfe7T46VwLtI2oFEAIFAy279t12xfVpsV94812zkkH6D0JyEVX1D6xxHWRcVctBnZ/YUkAtKtqV2eZp+v9b3O5Pvu1eOlcC7XDSaQo0r5h5sY9IKR6CNH96eyEdJW2kctd9RiuROzuPEJwPVLfMeURyn9u6iHCuB94AoNTn0O0eQukgJLOo/4NcpKDsPfQuiXRZ9u76Ah+0317dBZWg1XTzRj+f+XZVjJfAeF40GNSVvUGItDpoJERGHJiOhdd05GW4X0W/HGuoAXlYccEdc/Me+/Tsmx0rgfSAaDY2JeFhD0JYr4SuLOzD7q/1F/bfB5accR+vvX7lPcAJiGxgDt97psQCbHI/Dl+NxNOXdPI4HtdYn2xvvCyUAIIT42iIgw/E4jsdxPI67O45jG+1YjuV9LsdK4FiO5X0u95MS+Ll3egBWjsfRlONxNOU9N477JiZwLMdyLO+M3E+WwLEcy7G8A3KsBI7lWN7ncl8oASHED9s+BS8KIX76Hp3zohDid4UQzwghnhZC/Nd2+4YQ4jeFEC/Y/9fv0XgCIcQ3hBC/bv9+SAjxZXtPfkUIEd/uGG/DGNaEEP/M9pR4VgjxPe/E/RBC/DX7TJ4SQvySEKJzr+7HIX02Ft4DYeT/tGP6thDiE3d5HHen34fW+h39BwTAS8DDQAx8C/jwPTjvWeAT9vMypn/Ch4H/Gfhpu/2ngb91j+7DfwP8U+DX7d+/CvwF+/n/Av7yPRjDPwL+S/s5Btbu9f3AsFO/AnS9+/AT9+p+AD8AfAJ4ytu28B4AP4Jh2hbAdwNfvsvj+PeB0H7+W944PmznTQI8ZOdTcMfnutsv1h1c7PcAv+H9/TPAz7wD4/gXwJ8BngPO2m1ngefuwbkvAL8N/Gng1+1Ldct74I17dJfGsGonn2htv6f3g5q2fgMDa/914Av38n4Al1qTb+E9AP4e8OOL9rsb42h99+eBX7SfG3MG+A3ge+70PPeDO3DHvQrultjmKh8Hvgyc1lpft1/dAE7fgyH87xjiVlelcwLY1zWf1r24Jw8B28A/tG7J3xdC9LnH90NrfRX4X4HLwHXgAPg69/5++HLYPXgn39031e9jkdwPSuAdFSHEEvD/AH9Vaz3wv9NGrd7VHKoQ4kf///bOphWiMIrjv1Neig12ygIlW8pCsVBsWFgrxcKnkJUvoCzExkqiSJq1l7VQQiFEsfCysrY4FueZXC8TE/eZqXt+NTVz79Rz5n+f/j3n3LnPAZ5U9TDNcX5BBbb8nFfVTuxZjg/1mUh61GOdrFqwHatr+doGr2TE0OAn/tLv4zvKwQRK1qtARCoxA1hW1Y1w+FFEGsP5RuAp5TB6gGERuQVWsZRgFqiT/A6bcTS5B+5VdS98XsdMIbYeA8CNqj6rPYO8gWkUW48khTSIPncT/T5GgyH9OY5yMIF9oC1Uf6uwhqa5tAcV2yt9EThT1ZnEqRwwHt6PY7WC1FDVSVVtUtVm7LfvqOoosMt7j8cYcTwAdyLSHg71Y1vHR9UDSwO6RaQmXKN8HFH1+EQhDXLAWLhL0A28JNKGfyfR72NYv/b7GBGRahFpodh+H2kWeYoogAxh1flrYCrSmL3Ysu4YOAqvISwf3wYugS2gIaIOfbzfHWgNF/IKWAOqI4zfARwETTaB+lLoAUwD58ApsIRVvaPoAaxgtYhXbHU0UUgDrIA7F+btCdCVchxXWO6fn68Lie9PhTgugMFixvK/DTtOximHdMBxnBLiJuA4GcdNwHEyjpuA42QcNwHHyThuAo6TcdwEHCfjvAFae/jQxD9E2AAAAABJRU5ErkJggg==\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": [],
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "n94j1vMGrD8u",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 281
        },
        "outputId": "bc1f2c14-170e-4068-9572-f8459ecfb872"
      },
      "source": [
        "plt.hist(np.ravel(scan), 300)\n",
        "plt.title(\"Distribution of Hounsfield Units (HU)\")\n",
        "plt.show()"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEICAYAAACqMQjAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAYcklEQVR4nO3dedRkdX3n8fcHENxl6yEIaKMyRjQTNYySMYkGFBp00s4cFzxGWmTkGDUTR53YLhn3uExcwkRNUDgCzoiEZMZWMYgLMZmMQGsUBCW0LLK00NIIrij4nT/ur6G6qN+zNE8/W79f59TpW79772+pW3U/dZenOlWFJEmT7LTQHZAkLV6GhCSpy5CQJHUZEpKkLkNCktRlSEiSugyJZS7JXyb5kzmq6yFJfpRk5/b8vCT/aS7qbvV9NsmauapvFu2+Lcn3k3xvvtu+J0b7Pb5tplnvhUn+cYr5c7pdx+r+7SSXzWF970jyirmqr9X5niR/MJd1LmWGxBKW5KokP03ywyQ/SPJPSV6S5M7tWlUvqaq3zrCup061TFV9t6ruX1V3zEHf35TkY2P1H1VVp97TumfZj4cArwIOrqpfmTD/KUmunVC+3XakMzHe77ncNtO0+9EkbxsrW5mkkuwy3fpV9Q9V9ciRdad9303RlxXAscBftefTbqteQI7148+A1yXZdVv6tdwYEkvfv6+qBwAPBd4JvAY4ea4bmckOYIl6CHBTVd240B2ZpaXa77n0QuDsqvrpXFZaVRuBbwO/N5f1LlWGxDJRVbdU1TrgucCaJI+Brb/5Jdk7yafbUcfmJP+QZKckpzPsdD7VTln88ci3w+OTfBf4Yucb48OTXJDk1iSfTLJna+tu3+q2fFtLsgp4HfDc1t432vzRb3w7JXlDkquT3JjktCQPavO29GNNku+2Uy6v7702SR7U1t/U6ntDq/+pwLnAg1s/Prqtr3+SFyfZ0F7XdUkePNbXXUaWvds32yR/luTmJFcmOWpk2RcmuaIdLV6Z5PmT+j3eThvzyUk2Jrkuw6mpiaeikjwtybeT3JLkL4Bs6+vQ6rsqyauTXNTq/ESSe7d5d74vOu+7eyf5WJKb2vv0wiT7dJo6Cvj7e9LXKZwHPH071b2kGBLLTFVdAFwL/PaE2a9q81YA+zDsqKuqXgB8l+Go5P5V9e6RdZ4MPAo4stPkscCLgH2B24ETZ9DHvwP+FPhEa+/XJyz2wvb4XeBhwP2Bvxhb5reARwKHA/8tyaM6Tf4P4EGtnie3Ph9XVZ9n2NFc3/rxwun6PkmSw4B3AM9heB2uBs6YRRVPBC4D9gbeDZycwf0YXs+j2tHivwO+PsN+f5RhezwCeBxwBHC302NJ9gb+FnhDa/87wJNm0fee5wCrgAOBf8OwLbfSed+tYdhWBwB7AS8BekcKv8bwum0P3wImvS93OIbE8nQ9sOeE8l8w7MQeWlW/aOeHp/vxrjdV1Y+nOKQ/vaq+WVU/Bv4EeE7vG+ssPR94b1VdUVU/Al4LHDN2FPPmqvppVX0D+AYTPtStL8cAr62qH1bVVcB7gBfMoi8Pbt9q73wwBNRoX0+pqq9V1W2tr7+ZZOUM67+6qj7criecyrCNtnx7/iXwmCT3qaqNVXXJdJW1b95HA69o2+5G4H0Mr8O4o4FLquqsqvoF8H5gLi7gn1hV11fVZuBTwGNnuN4vGMLhEVV1R1V9tapu7Sy7O/DDsbLpttVM/bDVv8MzJJan/YDNE8r/O7AB+Fw7hbF2BnVdM4v5VwP3YvhGek89uNU3Wvcu3LXzhK13Zj9hONoYt3fr03hd+82iL9dX1e6jD2D04udWfW2hdtMs2rhzHFX1kzZ5/xa8z2X4Nr0xyWeS/OoM6nsow5g3juwo/wr4VxOWfTAj27B9aZhqm9/e6h51L4Yw++WkMdHfNpOcDpwDnJHk+iTvTjLe3hY3Aw8YK5tuW03q/5Yx/GLk+QOAH8ywz8uaIbHMJPm3DDunu93B0b5Jv6qqHsZwUe6VSQ7fMrtT5XRHGgeMTD+E4YP2feDHwH1H+rUzw2mumdZ7PcPObrTu24Ebpllv3Pdbn8brum6W9Uxlq76200R7tTZ+3IrvO7L83e6i6qmqc6rqaQxHF98GPjyD1a4BbgP2HtlZPrCqHj1h2Y2MbMMkYettOu67wMqxsgOBa6rql3dffFpbvQ/aEe6bq+pghtNrz2A4PTjJRcC/nmV73wUe0sYJQJL7MgTo6BeJRzEcne7wDIllIskDkzyD4Vz4x6rq4gnLPCPJI9oH5BbgDu769ncDwzn72fr9JAe3D9pbgLPaaZN/Ae6d5Ontm+AbgN1G1rsBWJmR23XHfBz4L0kOTHJ/7rqGcftsOtf6cibw9iQPSPJQ4JXAx6Zec1Y+DhyX5LFJdmt9Pb+qrqqqTQxh8ftJdk7yIuDhM6k0yT5JVrfQuQ34EVt/W5+o3Z3zOeA97X2xU5KHJ3nyhMU/Azw6yX9sp/L+M1OH2N8AT09yRBvPgxm27WyuwYza6n2X5HeT/Fr7UnErQ8D3xnw2wzWm2Tgf+Bmwtl0kvx/DXYHr2Tokngx8dpZ1L0uGxNL3qSQ/ZPj2+HrgvcBxnWUPAj7PsLP5f8AHq+pLbd47gDe00xOvnkX7pzNcJP0ecG+GnQxVdQvwUuAj3PWNevRup79u/96U5GsT6j2l1f1l4EqGD/YfzqJfo/6wtX8FwxHW/2r1z4l2IflPGHagGxlCYPT8/4uB/8pwCurRwD/NsOqdGALteobTh08GZvpHXscCuwKXMpyWOYvhaGS8798Hns2wo7yJ4T3yf3uVtmsiz2N4v2xmeB+dD7x5hv0aN/6++5XW11sZLh7/PcP7YJLTgKOT3GemjbVrRk8HnsLwfryC4ZTbc7Zcn0uyL3Aw8H+2ZUDLTfxPhyQtVUn+FLixqt4/h3W+B/hOVX1wrupcygwJSVKXp5skSV2GhCSpy5CQJHUtux9t23vvvWvlypUL3Q1JWlK++tWvfr+qVoyXL7uQWLlyJevXr1/obkjSkpLk6knlnm6SJHUZEpKkLkNCktRlSEiSugwJSVKXISFJ6jIkJEldhoQkqcuQkCR1GRKztHLtZxa6C5I0bwwJSVKXISFJ6jIkJEldhoQkqcuQkCR1GRKSpC5DQpLUZUhIkroMCUlSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHXNOCSS7Jzkn5N8uj0/MMn5STYk+USSXVv5bu35hjZ/5Ugdr23llyU5cqR8VSvbkGTtSPnENiRJ82M2RxJ/BHxr5Pm7gPdV1SOAm4HjW/nxwM2t/H1tOZIcDBwDPBpYBXywBc/OwAeAo4CDgee1ZadqQ5I0D2YUEkn2B54OfKQ9D3AYcFZb5FTgmW16dXtOm394W341cEZV3VZVVwIbgCe0x4aquqKqfg6cAayepg1J0jyY6ZHE+4E/Bn7Znu8F/KCqbm/PrwX2a9P7AdcAtPm3tOXvLB9bp1c+VRuSpHkwbUgkeQZwY1V9dR76s02SnJBkfZL1mzZtWujuSNKyMZMjiScBv5fkKoZTQYcBfw7snmSXtsz+wHVt+jrgAIA2/0HATaPlY+v0ym+aoo2tVNVJVXVIVR2yYsWKGQxJkjQT04ZEVb22qvavqpUMF56/WFXPB74EPKsttgb4ZJte157T5n+xqqqVH9PufjoQOAi4ALgQOKjdybRra2NdW6fXhiRpHtyTv5N4DfDKJBsYrh+c3MpPBvZq5a8E1gJU1SXAmcClwN8BL6uqO9o1h5cD5zDcPXVmW3aqNiRJ82CX6Re5S1WdB5zXpq9guDNpfJmfAc/urP924O0Tys8Gzp5QPrENSdL88C+uJUldhoQkqcuQkCR1GRKSpC5DQpLUZUhIkroMCUlSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHUZEpKkLkNCktRlSEiSugwJSVKXISFJ6jIkJEldhoQkqcuQkCR1GRKSpC5DQpLUZUhIkroMCUlSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1DVtSCS5d5ILknwjySVJ3tzKD0xyfpINST6RZNdWvlt7vqHNXzlS12tb+WVJjhwpX9XKNiRZO1I+sQ1J0vyYyZHEbcBhVfXrwGOBVUkOBd4FvK+qHgHcDBzflj8euLmVv68tR5KDgWOARwOrgA8m2TnJzsAHgKOAg4HntWWZog1J0jyYNiRq8KP29F7tUcBhwFmt/FTgmW16dXtOm394krTyM6rqtqq6EtgAPKE9NlTVFVX1c+AMYHVbp9eGJGkezOiaRPvG/3XgRuBc4DvAD6rq9rbItcB+bXo/4BqANv8WYK/R8rF1euV7TdHGeP9OSLI+yfpNmzbNZEiSpBmYUUhU1R1V9Vhgf4Zv/r+6XXs1S1V1UlUdUlWHrFixYqG7I0nLxqzubqqqHwBfAn4T2D3JLm3W/sB1bfo64ACANv9BwE2j5WPr9MpvmqINSdI8mMndTSuS7N6m7wM8DfgWQ1g8qy22Bvhkm17XntPmf7GqqpUf0+5+OhA4CLgAuBA4qN3JtCvDxe11bZ1eG5KkebDL9IuwL3BquwtpJ+DMqvp0kkuBM5K8Dfhn4OS2/MnA6Uk2AJsZdvpU1SVJzgQuBW4HXlZVdwAkeTlwDrAzcEpVXdLqek2nDUnSPJg2JKrqIuBxE8qvYLg+MV7+M+DZnbreDrx9QvnZwNkzbUOSND/8i2tJUpchIUnqMiQkSV2GhCSpy5CQJHUZEpKkLkNCktRlSEiSugwJSVKXISFJ6jIkJEldhoQkqcuQkCR1GRKSpC5DQpLUZUhIkroMCUlSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHUZEpKkLkNCktRlSEiSugwJSVKXISFJ6jIkJEldhoQkqWvakEhyQJIvJbk0ySVJ/qiV75nk3CSXt3/3aOVJcmKSDUkuSvL4kbrWtOUvT7JmpPw3klzc1jkxSaZqQ5I0P2ZyJHE78KqqOhg4FHhZkoOBtcAXquog4AvtOcBRwEHtcQLwIRh2+MAbgScCTwDeOLLT/xDw4pH1VrXyXhuSpHkwbUhU1caq+lqb/iHwLWA/YDVwalvsVOCZbXo1cFoNvgLsnmRf4Ejg3KraXFU3A+cCq9q8B1bVV6qqgNPG6prUhiRpHszqmkSSlcDjgPOBfapqY5v1PWCfNr0fcM3Iate2sqnKr51QzhRtjPfrhCTrk6zftGnTbIYkSZrCjEMiyf2BvwFeUVW3js5rRwA1x33bylRtVNVJVXVIVR2yYsWK7dkNSdqhzCgkktyLISD+Z1X9bSu+oZ0qov17Yyu/DjhgZPX9W9lU5ftPKJ+qDUnSPJjJ3U0BTga+VVXvHZm1Dthyh9Ia4JMj5ce2u5wOBW5pp4zOAY5Iske7YH0EcE6bd2uSQ1tbx47VNakNSdI82GUGyzwJeAFwcZKvt7LXAe8EzkxyPHA18Jw272zgaGAD8BPgOICq2pzkrcCFbbm3VNXmNv1S4KPAfYDPtgdTtCFJmgfThkRV/SOQzuzDJyxfwMs6dZ0CnDKhfD3wmAnlN01qQ5I0P/yLa0lSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHUZEpKkLkNCktRlSEiSugwJSVKXISFJ6jIkJEldhoQkqcuQkCR1GRKSpC5DQpLUZUhIkroMCUlSlyEhSeoyJCRJXYaEJKnLkJAkdRkSkqQuQ0KS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpa9qQSHJKkhuTfHOkbM8k5ya5vP27RytPkhOTbEhyUZLHj6yzpi1/eZI1I+W/keTits6JSTJVG5Kk+TOTI4mPAqvGytYCX6iqg4AvtOcARwEHtccJwIdg2OEDbwSeCDwBeOPITv9DwItH1ls1TRuSpHkybUhU1ZeBzWPFq4FT2/SpwDNHyk+rwVeA3ZPsCxwJnFtVm6vqZuBcYFWb98Cq+kpVFXDaWF2T2pAkzZNtvSaxT1VtbNPfA/Zp0/sB14wsd20rm6r82gnlU7VxN0lOSLI+yfpNmzZtw3C2zcq1n5m3tiRpIdzjC9ftCKDmoC/b3EZVnVRVh1TVIStWrNieXZGkHcq2hsQN7VQR7d8bW/l1wAEjy+3fyqYq339C+VRtSJLmybaGxDpgyx1Ka4BPjpQf2+5yOhS4pZ0yOgc4Iske7YL1EcA5bd6tSQ5tdzUdO1bXpDYkSfNkl+kWSPJx4CnA3kmuZbhL6Z3AmUmOB64GntMWPxs4GtgA/AQ4DqCqNid5K3BhW+4tVbXlYvhLGe6gug/w2fZgijYkSfNk2pCoqud1Zh0+YdkCXtap5xTglAnl64HHTCi/aVIbkqT5419cS5K6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHUZEpKkLkNCktRlSEiSugyJbeTPhEvaERgS95BhIWk5MyQkSV2GhCSpy5CQJHUZEnPMaxSSlhNDYjswKJYWt5fUZ0hIkroMCUlSlyGhHZqnmqSpGRKSpC5DYg75rXT+zdVr7raTJjMktMMyGKTpGRLbiTsgScuBIbEdjQaFoSFpKTIk5oEBsX1ty+s7aR23k3R3hsR25o5H0lJmSMwjA2Nu+XpK258hMQfcWS08t4G0fRgS0gjDRtqaITFHZrpzcSe08NwG0swZEgvAndT24esqzT1DYhu4M1q8pts2M9l2bl/pLrssdAekbTHbHbk7fmnbeCSxQNxpbT9bXltfY+meMyQWkDux7Wc0KHydpW1nSCwwd2CLk9tFGhgSs7C9dhzukGbH10uaP4v+wnWSVcCfAzsDH6mqdy5wl7aL8R3fVe98+gL1RFts2SZuC+3IFnVIJNkZ+ADwNOBa4MIk66rq0oXt2fZnaEy2EEcRK9d+xtdfO6xFHRLAE4ANVXUFQJIzgNXAsg+JcZ5iWViL4fXfElQe4Wg+LfaQ2A+4ZuT5tcATxxdKcgJwQnv6oySXzUPfAPYGvj9Pbc2n5TouWMJjy7umfs4SHts0luu4YHGN7aGTChd7SMxIVZ0EnDTf7SZZX1WHzHe729tyHRc4tqVouY4LlsbYFvvdTdcBB4w837+VSZLmwWIPiQuBg5IcmGRX4Bhg3QL3SZJ2GIv6dFNV3Z7k5cA5DLfAnlJVlyxwt0bN+ymuebJcxwWObSlaruOCJTC2VNVC90GStEgt9tNNkqQFZEhIkroMiY4kz05ySZJfJjlkbN5rk2xIclmSI0fKV7WyDUnWjpQfmOT8Vv6JdhF+UUjypiTXJfl6exw9Mm9W41zMlmKfxyW5KsnFbTutb2V7Jjk3yeXt3z1aeZKc2MZ7UZLHL2zvt5bklCQ3JvnmSNmsx5JkTVv+8iRrFmIsozrjWtqfsaryMeEBPAp4JHAecMhI+cHAN4DdgAOB7zBcVN+5TT8M2LUtc3Bb50zgmDb9l8AfLPT4RsbzJuDVE8pnPc7F+liKfe6M4ypg77GydwNr2/Ra4F1t+mjgs0CAQ4HzF7r/Y/3+HeDxwDe3dSzAnsAV7d892vQei3BcS/oz5pFER1V9q6om/eX2auCMqrqtqq4ENjD8fMidPyFSVT8HzgBWJwlwGHBWW/9U4JnbfwT32KzGuYD9nIml2OeZWs3wnoKt31urgdNq8BVg9yT7LkQHJ6mqLwObx4pnO5YjgXOranNV3QycC6za/r3v64yrZ0l8xgyJ2Zv0UyH7TVG+F/CDqrp9rHwxeXk7jD9lyyE+sx/nYrYU+zxJAZ9L8tX2UzQA+1TVxjb9PWCfNr0UxzzbsSylMS7Zz9gOHRJJPp/kmxMey+VbJjDtOD8EPBx4LLAReM+CdlZT+a2qejxwFPCyJL8zOrOGcxjL4p725TQWlvhnbFH/Md32VlVP3YbVpvqpkEnlNzEcHu/Sjibm/adFZjrOJB8GPt2eznaci9my+HmXqrqu/Xtjkv/NcFrihiT7VtXGdgrmxrb4UhzzbMdyHfCUsfLz5qGfs1JVN2yZXoqfsR36SGIbrQOOSbJbkgOBg4AL6PyESPtG9CXgWW39NcAnF6DfE42dp/4PwJa7MmY1zvns8zZYin3eSpL7JXnAlmngCIZttY7hPQVbv7fWAce2O4MOBW4ZOZWzWM12LOcARyTZo53COaKVLSpL/jO2kHcCLOYHw8a8FrgNuAE4Z2Te6xnuPrgMOGqk/GjgX9q814+UP4xh428A/hrYbaHHN9K304GLgYsY3oj7bus4F/NjKfZ5rP8PY7jL5RvAJVvGwHDN6wvA5cDngT1beRj+w67vtO17yEL1vTOejzOcevlF+5wdvy1jAV7UPlcbgOMW6biW9GfMn+WQJHV5ukmS1GVISJK6DAlJUpchIUnqMiQkSV2GhCSpy5CQJHX9f4jNmWH8dSzWAAAAAElFTkSuQmCC\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": [],
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "7U0SbwsYrPR_",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 136
        },
        "outputId": "d7d6e14d-81f5-4a5b-a231-cefe7e2e815d"
      },
      "source": [
        "def Aw(ct_scan):\n",
        "\n",
        "    a_map =  ((ct_scan / 1000) + 1)\n",
        "    aw = np.sum(a_map) * 0.2\n",
        "    print(aw)\n",
        "\n",
        "def Dw(ct_scan):\n",
        "\n",
        "    \n",
        "\n",
        "Aw(scan)"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "error",
          "ename": "IndentationError",
          "evalue": "ignored",
          "traceback": [
            "\u001b[0;36m  File \u001b[0;32m\"<ipython-input-4-6a99dcd5c2ba>\"\u001b[0;36m, line \u001b[0;32m11\u001b[0m\n\u001b[0;31m    Aw(scan)\u001b[0m\n\u001b[0m     ^\u001b[0m\n\u001b[0;31mIndentationError\u001b[0m\u001b[0;31m:\u001b[0m expected an indented block\n"
          ]
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "kJZf0TOYxzE_",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 55
        },
        "outputId": "6d6586bd-614e-4978-eb85-d28c6b3f99e7"
      },
      "source": [
        "!pip install pydicom"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "Requirement already satisfied: pydicom in /usr/local/lib/python3.6/dist-packages (2.0.0)\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Z6c9E16w1THv",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 143
        },
        "outputId": "15784750-4429-4e3b-8b31-b073674f12df"
      },
      "source": [
        "from water import DICOMwaterequivalent\n",
        "DICOMwaterequivalent(scan_path, 0, window = False)"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "(28665.53275358887,\n",
              " 191.0447326485878,\n",
              " 31635.448450780437,\n",
              " 200.6975435399357,\n",
              " 64525.59102245011,\n",
              " 286.62961141025136,\n",
              " False)"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 16
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "sN1PuOnbx2pp",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 1000
        },
        "outputId": "7e53912e-a40a-44f2-8cdb-99cc7c469750"
      },
      "source": [
        "from water import DICOMwaterequivalent\n",
        "\n",
        "tuples = []\n",
        "for i in range(0, 233):\n",
        "\n",
        "  print(i)\n",
        "\n",
        "  if (len(str(i)) == 1):\n",
        "\n",
        "    scan_path = '/content/drive/My Drive/volume2/00000{}.dcm'.format(str(i))\n",
        "\n",
        "  elif (len(str(i)) == 2):\n",
        "\n",
        "    scan_path = '/content/drive/My Drive/volume2/0000{}.dcm'.format(str(i))\n",
        "\n",
        "  elif (len(str(i)) == 3):\n",
        "\n",
        "    scan_path = '/content/drive/My Drive/volume2/000{}.dcm'.format(str(i))\n",
        "\n",
        "  tuples.append(DICOMwaterequivalent(scan_path, 0, window = False))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "0\n",
            "1\n",
            "2\n",
            "3\n",
            "4\n",
            "5\n",
            "6\n",
            "7\n",
            "8\n",
            "9\n",
            "10\n",
            "11\n",
            "12\n",
            "13\n",
            "14\n",
            "15\n",
            "16\n",
            "17\n",
            "18\n",
            "19\n",
            "20\n",
            "21\n",
            "22\n",
            "23\n",
            "24\n",
            "25\n",
            "26\n",
            "27\n",
            "28\n",
            "29\n",
            "30\n",
            "31\n",
            "32\n",
            "33\n",
            "34\n",
            "35\n",
            "36\n",
            "37\n",
            "38\n",
            "39\n",
            "40\n",
            "41\n",
            "42\n",
            "43\n",
            "44\n",
            "45\n",
            "46\n",
            "47\n",
            "48\n",
            "49\n",
            "50\n",
            "51\n",
            "52\n",
            "53\n",
            "54\n",
            "55\n",
            "56\n",
            "57\n",
            "58\n",
            "59\n",
            "60\n",
            "61\n",
            "62\n",
            "63\n",
            "64\n",
            "65\n",
            "66\n",
            "67\n",
            "68\n",
            "69\n",
            "70\n",
            "71\n",
            "72\n",
            "73\n",
            "74\n",
            "75\n",
            "76\n",
            "77\n",
            "78\n",
            "79\n",
            "80\n",
            "81\n",
            "82\n",
            "83\n",
            "84\n",
            "85\n",
            "86\n",
            "87\n",
            "88\n",
            "89\n",
            "90\n",
            "91\n",
            "92\n",
            "93\n",
            "94\n",
            "95\n",
            "96\n",
            "97\n",
            "98\n",
            "99\n",
            "100\n",
            "101\n",
            "102\n",
            "103\n",
            "104\n",
            "105\n",
            "106\n",
            "107\n",
            "108\n",
            "109\n",
            "110\n",
            "111\n",
            "112\n",
            "113\n",
            "114\n",
            "115\n",
            "116\n",
            "117\n",
            "118\n",
            "119\n",
            "120\n",
            "121\n",
            "122\n",
            "123\n",
            "124\n",
            "125\n",
            "126\n",
            "127\n",
            "128\n",
            "129\n",
            "130\n",
            "131\n",
            "132\n",
            "133\n",
            "134\n",
            "135\n",
            "136\n",
            "137\n",
            "138\n",
            "139\n",
            "140\n",
            "141\n",
            "142\n",
            "143\n",
            "144\n",
            "145\n",
            "146\n",
            "147\n",
            "148\n",
            "149\n",
            "150\n",
            "151\n",
            "152\n",
            "153\n",
            "154\n",
            "155\n",
            "156\n",
            "157\n",
            "158\n",
            "159\n",
            "160\n",
            "161\n",
            "162\n",
            "163\n",
            "164\n",
            "165\n",
            "166\n",
            "167\n",
            "168\n",
            "169\n",
            "170\n",
            "171\n",
            "172\n",
            "173\n",
            "174\n",
            "175\n",
            "176\n",
            "177\n",
            "178\n",
            "179\n",
            "180\n",
            "181\n",
            "182\n",
            "183\n",
            "184\n",
            "185\n",
            "186\n",
            "187\n",
            "188\n",
            "189\n",
            "190\n",
            "191\n",
            "192\n",
            "193\n",
            "194\n",
            "195\n",
            "196\n",
            "197\n",
            "198\n",
            "199\n",
            "200\n",
            "201\n",
            "202\n",
            "203\n",
            "204\n",
            "205\n",
            "206\n",
            "207\n",
            "208\n",
            "209\n",
            "210\n",
            "211\n",
            "212\n",
            "213\n",
            "214\n",
            "215\n",
            "216\n",
            "217\n",
            "218\n",
            "219\n",
            "220\n",
            "221\n",
            "222\n",
            "223\n",
            "224\n",
            "225\n",
            "226\n",
            "227\n",
            "228\n",
            "229\n",
            "230\n",
            "231\n",
            "232\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "3BMff4-q0W-W",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "water_equivalent_areas = []\n",
        "water_equivalent_diameters = []\n",
        "roi_area_mms = []\n",
        "roi_area_diameters = []\n",
        "roi_hull_areas = []\n",
        "roi_hull_diameters = []\n",
        "imgs = []\n",
        "\n",
        "for item in tuples:\n",
        "\n",
        "  water_equivalent_areas.append(item[0])\n",
        "  water_equivalent_diameters.append(item[1])\n",
        "  roi_area_mms.append(item[2])\n",
        "  roi_area_diameters.append(item[3])\n",
        "  roi_hull_areas.append(item[4])\n",
        "  roi_hull_diameters.append(item[5])\n",
        "  imgs.append(item[6])"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "uq4-GZdK06UK",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 545
        },
        "outputId": "1bed5642-cbef-4f56-9ab1-732d777ba61e"
      },
      "source": [
        "plt.hist(np.array(water_equivalent_areas), bins = 50)\n",
        "plt.title(\"Water Equivalent Areas Distribution\")\n",
        "plt.show()\n",
        "\n",
        "plt.hist(np.array(water_equivalent_diameters), bins = 50)\n",
        "plt.title(\"Water Equivalent Diameters Distribution\")\n",
        "plt.show()"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAEICAYAAABGaK+TAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAWN0lEQVR4nO3debxkZX3n8c9XVgXCIj3YgNAYHSI6IzIdhYhKEpXFBbOMAkaJo0HjMnESx0GdRByTvJhMMkYnjkrcMAouEaIjGGDQAC5BG4IRRGRrZKcRWUSjAr/8cZ4L1cW9fW/frnv7Ps3n/XrV6556zqlzflV16ltPPadO3VQVkqT+PGxjFyBJmh8DXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4Ji7JW5J8YALrOS7JxyZR06YsyR5Jfphkswmt731J/rBNH5Tkukmst63v6Ukum9T6HuoM8AWS5M1JvjDWdvkMbUfMsq6JvohG1nlfe+GPXg7Y0HVX1Z9W1SsnUeckJFmRpJJsPodlf7st++LFqG02rZ57R56fq5N8OMm/nVqmqr5XVdtW1b1zWNeXZ9tmVb26qt4xoforyWNH1n1eVe09iXXLAF9I5wK/NNUrSrIc2AJ48ljbY9uyC2YdwXVDe+GPXr62kLV04GjgNuBl61poLm8GE/S1qtoW2B54FvBj4IIkT5z0hibVi9fiMMAXzjcYAnvfdv3pwJeAy8barqyqG5K8PMmlSe5KclWSVwEk2Qb4ArDrSC9s1yQPS3JskiuTfD/Jp5Ls1G4z1eN8RZLvAV9c3+KT7JXknFbPWUn+amo4Y7pPBElWJ3lWmz5uZNkvJHnd2LLfTPLrbfpdSa5NcmeSC5I8fR017Z/kq0lub+s4aGTePyR5R5KvtJrPTLJzmz31Bnn7uj5lJNkTeCZwDHBwkkeNzDsoyXVJ/luSm4APr+s5aLf5dJKbktyR5NwkTxiZd1iSb7dar0/yxpnu95Squreqrqyq1wDnAMe1da31CaP1tK9q6746yUuSPB54H3BAewxub8t+JMl7k5ye5G7gl1vbH489Nm9Jcmt7nl8y9ri/cuT6/b38JFOP+zfbNl88vu8keXxbx+1JLknygpF5H0nyniSntftyfpKfn+1xeigxwBdIVf0UOB94Rmt6BnAe8OWxtqmd/BbgecDPAS8H3plkv6q6GziUtXvLNwCvB17IEDi7Aj8A3jNWxjOBxwMHz+MunARcAOwMvIOhZzofJwNHTl1Jsg+wJ3Baa/oGwxvaTm2bn06y9fhKkuzWbvPHbdk3Ap9JsmxksaMYHrt/A2zZloEHHu8dZvmU8TJgVVV9BrgUeMnY/Ee1be/JEPKzPQdfAB7X6rkQ+PjIvA8Cr6qq7YAnsv5vsqcwdADW0t7w3w0c2tb9S8BFVXUp8Gpab76qdhi52VHAnwDbMeyf4x7FsB/sxrAfnJBk1mGQqpp63J/UtvnJsVq3AP4fcCbDY/R64ONj6z4CeDuwI3BFq1ONAb6wzuGB8Hg6Q4CfN9Z2DkBVndZ6V1VV5zDs1DP2RhlejG+tquuq6icMvbHfzNof7Y+rqrur6sczrGPX1vMZvWyTZA/gF4E/rKqfVNW5DC+0+TgV2Lf1bmEIxVNazVTVx6rq+1V1T1X9BbAVMF04/BZwelWdXlX3VdVZwCrgsJFlPlxV323391M88Elnrl7G8CZC+zs+jHIf8Lb2mPyYWZ6DqvpQVd01Mu9JSbZv6/oZsE+Sn6uqH1TVhetZ6w0MbybTuQ94YpKHV9WNVXXJLOv6bFV9pT2u/zLDMlP7wjkMb6QvWs96p7M/sC1wfFX9tKq+CHyekTd84NSq+npV3cPwBri+z+kmzQBfWOcCB7aP1cuq6nLgqwxj4zsx9LzOBUhyaJJ/THJb+3h7GEOvZyZ7AqdOBS9Dj/FeYJeRZa6dpb4bqmqHscvdtN5km55yzdzv9gOq6i6GF/zUgdojGemJJnljhqGjO9r92J7p7/eewH8cfbMBDgSWjyxz08j0jxjCYU6SPA3YC/hEazoJ+HdJRgNjzVjAzfgcJNksyfFteOVOYHW7zdR9+w2G5/iaDENV63vweDeGsfq1tOfsxQxvLje24YdfmGVds+0n0+0Lu65PsTPYFbi2qu4bW/duI9fn/Zw+FBjgC+trDIH0O8BXAKrqTobe0+8wBOjVSbYCPgP8ObBL+3h7OpC2nul+MvJaho/Jo+G7dVVdP7LMfH9q8kZgx/ZxfMoeI9N3A4+YupLhwNfoUMa4k4EjW0htzXAsgAzj3W9i6M3t2O73HTxwv0ddC/zN2P3dpqqOn8P9mcvjcHTb7kVtjPv8kfaZ1rOu5+Ao4HCGg47bAyvabQJQVd+oqsMZhg7+juETw/r4NYZPcw9SVWdU1bMZ3ty+A/z1DPXPdL/GTbcv3NCm19oXGIZb5uoG4NFJRnNoD+D6GZbXGAN8AbWP2auA32ftF9uXW9vU+PeWDEMHa4B7khwKPGdk+ZuBR458/IbhgNSfTA1NJFmW5PAJ1X1Nq/vtSbZMciDw/JFFvgtsneS5bRzzv7f6Z3I6Q2/1fwCfHOlxbQfcw3C/N0/yRwzHAKbzMeD5SQ5uvdut2wGx3edwl9YwDCs8ZrqZbcz9RQzj2vuOXF4PHJWZv3GyrudgO+AnwPcZAu5PR7a3ZTuwuH1V/Qy4s9W3Tu1+75Xk/wAHMYwNjy+zS5LDW+D+BPjhyLpvBnZPsuVs25rG1L7wdIZjNZ9u7RcBv57kERm+LviKsdvdzAyPO8Ob5I+ANyXZIsNB6efzwKcgzcIAX3jnMPSyRg8OndfazoX7hxn+M0Mv7AcMvbfPTS1cVd9h6MVe1T6u7wq8qy1zZpK7gH8EnrqetY1+s2Xq8htt3lFtfbcBbwM+OlLPHcBrgA8w9JbuBmb8nnobAz6FoTd60sisM4C/Z3hDuAb4F2b4OF9V1zL0aN/CEMjXAv+VOezDVfUjhoNfX2mP3/5ji7yQ4at5H62qm6YuwIeAzYFDZlj1up6Dj7b7dD3w7TZv1EuB1W145dU8+IDpqAOS/JAh6P+B4U3uF6vqW9Ms+zCGzsENDM/dM4HfbfO+CFwC3JTk1nVsb9xNDPvlDQzDX69u+yTAO4GfMgT1iax9oBaGsf8T2+O+1rh5O9D/fIaD9LcC/xd42ci6NYuU/9BBc5DkOOCxVfVbG7sWSQN74JLUKQNckjrlEIokdcoeuCR1ajF/kIedd965VqxYsZiblKTuXXDBBbdW1YPOtVjUAF+xYgWrVq1azE1KUveSTHsmtEMoktQpA1ySOmWAS1KnDHBJ6pQBLkmdMsAlqVMGuCR1ygCXpE4Z4JLUqUU9E1PS0rHi2NOmbV99/HMXuRLNlz1wSeqUAS5JnTLAJalTBrgkdcoAl6ROGeCS1CkDXJI6ZYBLUqc8kUeaJ0+E0cZmD1ySOmWAS1KnDHBJ6pQBLkmdMsAlqVMGuCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTBrgkdWrWAE/y6CRfSvLtJJck+b3WvlOSs5Jc3v7uuPDlSpKmzKUHfg/wB1W1D7A/8Nok+wDHAmdX1eOAs9t1SdIimTXAq+rGqrqwTd8FXArsBhwOnNgWOxF44UIVKUl6sPUaA0+yAngycD6wS1Xd2GbdBOwyw22OSbIqyao1a9ZsQKmSpFFzDvAk2wKfAd5QVXeOzquqAmq621XVCVW1sqpWLlu2bIOKlSQ9YE4BnmQLhvD+eFWd0ppvTrK8zV8O3LIwJUqSpjOXb6EE+CBwaVX975FZnwOObtNHA5+dfHmSpJnM5b/SPw14KfCtJBe1trcAxwOfSvIK4BrgRQtToiRpOrMGeFV9GcgMs391suVIkubKMzElqVMGuCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTBrgkdcoAl6ROGeCS1CkDXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4JHXKAJekThngktQpA1ySOjXrf6XXwltx7GnTtq8+/rmLXImkntgDl6ROGeCS1CkDXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4JHXKAJekThngktQpA1ySOmWAS1KnDHBJ6pQBLkmdmjXAk3woyS1JLh5pOy7J9UkuapfDFrZMSdK4ufTAPwIcMk37O6tq33Y5fbJlSZJmM2uAV9W5wG2LUIskaT1syH/keV2SlwGrgD+oqh9Mt1CSY4BjAPbYY48N2Jym+B98JMH8D2K+F/h5YF/gRuAvZlqwqk6oqpVVtXLZsmXz3Jwkady8Aryqbq6qe6vqPuCvgadMtixJ0mzmFeBJlo9c/TXg4pmWlSQtjFnHwJOcDBwE7JzkOuBtwEFJ9gUKWA28agFrlCRNY9YAr6ojp2n+4ALUIklaD56JKUmdMsAlqVMGuCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTBrgkdcoAl6ROGeCS1CkDXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4JHXKAJekThngktQpA1ySOmWAS1KnDHBJ6pQBLkmdMsAlqVMGuCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTBrgkdcoAl6ROzRrgST6U5JYkF4+07ZTkrCSXt787LmyZkqRxc+mBfwQ4ZKztWODsqnoccHa7LklaRLMGeFWdC9w21nw4cGKbPhF44YTrkiTNYr5j4LtU1Y1t+iZgl5kWTHJMklVJVq1Zs2aem5Mkjdvgg5hVVUCtY/4JVbWyqlYuW7ZsQzcnSWrmG+A3J1kO0P7eMrmSJElzMd8A/xxwdJs+GvjsZMqRJM3VXL5GeDLwNWDvJNcleQVwPPDsJJcDz2rXJUmLaPPZFqiqI2eY9asTrkWStB48E1OSOmWAS1KnDHBJ6tSsY+CS+rbi2NM2dgkPGet6rFcf/9yJb88euCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTnsizhHkCRp9met7W90SOSa1Hmy574JLUKQNckjplgEtSpwxwSeqUAS5JnTLAJalTBrgkdcoAl6ROeSKPxMY9acYTtjRf9sAlqVMGuCR1ygCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5Jner+RJ6e/muJJ2yoZz291h4q7IFLUqcMcEnqlAEuSZ0ywCWpUwa4JHVqg76FkmQ1cBdwL3BPVa2cRFGSpNlN4muEv1xVt05gPZKk9eAQiiR1akN74AWcmaSA91fVCeMLJDkGOAZgjz322MDN6aFqUziJZFIncnlC2OT1un9taA/8wKraDzgUeG2SZ4wvUFUnVNXKqlq5bNmyDdycJGnKBgV4VV3f/t4CnAo8ZRJFSZJmN+8AT7JNku2mpoHnABdPqjBJ0rptyBj4LsCpSabWc1JV/f1EqpIkzWreAV5VVwFPmmAtkqT14NcIJalTBrgkdcoAl6ROdf8feaTpTOrEjIfiSTOTvM/ru66lfuLMlKWyX9gDl6ROGeCS1CkDXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4JHXKE3k2Ib3+V5HFtFROwNDGsak9//bAJalTBrgkdcoAl6ROGeCS1CkDXJI6ZYBLUqcMcEnqlAEuSZ3q5kSe9f0C/qb2hf2lYqmdLOTzrOk8VPYLe+CS1CkDXJI6ZYBLUqcMcEnqlAEuSZ0ywCWpUwa4JHXKAJekTnVzIo/mb6mdfKNNy8Y8aWaht73UTwiyBy5JnTLAJalTBrgkdcoAl6ROGeCS1KkNCvAkhyS5LMkVSY6dVFGSpNnNO8CTbAa8BzgU2Ac4Msk+kypMkrRuG9IDfwpwRVVdVVU/BT4BHD6ZsiRJs0lVze+GyW8Ch1TVK9v1lwJPrarXjS13DHBMu7o3cNn8y11QOwO3buwi5sG6F5d1L64e616ImvesqmXjjQt+JmZVnQCcsNDb2VBJVlXVyo1dx/qy7sVl3Yurx7oXs+YNGUK5Hnj0yPXdW5skaRFsSIB/A3hckr2SbAkcAXxuMmVJkmYz7yGUqronyeuAM4DNgA9V1SUTq2zxLflhnhlY9+Ky7sXVY92LVvO8D2JKkjYuz8SUpE4Z4JLUqU0qwJM8OsmXknw7ySVJfq+175TkrCSXt787tvYkeXf7KYB/TrLfyLqObstfnuTokfb/kORb7TbvTpIJ1L11kq8n+War++2tfa8k57dtfbIdLCbJVu36FW3+ipF1vbm1X5bk4JH2BfvZgySbJfmnJJ/vpe4kq9vzeFGSVa1tSe8nbb07JPnbJN9JcmmSA5Z63Un2bo/z1OXOJG/ooO7/0l6PFyc5OcPrdGnt21W1yVyA5cB+bXo74LsMp/n/GXBsaz8W+J9t+jDgC0CA/YHzW/tOwFXt745tesc27+tt2bTbHjqBugNs26a3AM5v2/gUcERrfx/wu236NcD72vQRwCfb9D7AN4GtgL2AKxkOMG/Wph8DbNmW2WeCj/vvAycBn2/Xl3zdwGpg57G2Jb2ftPWeCLyyTW8J7NBD3SP1bwbcBOy5lOsGdgOuBh4+sk//9lLbtycaoEvtAnwWeDbD2Z/LW9ty4LI2/X7gyJHlL2vzjwTeP9L+/ta2HPjOSPtay02o5kcAFwJPZTiba/PWfgBwRps+AzigTW/elgvwZuDNI+s6o93u/tu29rWW28B6dwfOBn4F+Hyro4e6V/PgAF/S+wmwPUOopKe6x2p9DvCVpV43Q4Bfy/BmsXnbtw9eavv2JjWEMqp9hHkyQ292l6q6sc26CdilTU89SVOua23rar9umvZJ1LtZkouAW4CzGN6db6+qe6bZ1v31tfl3AI+cx/2ZhL8E3gTc164/spO6CzgzyQUZfu4Blv5+shewBvhwhiGrDyTZpoO6Rx0BnNyml2zdVXU98OfA94AbGfbVC1hi+/YmGeBJtgU+A7yhqu4cnVfD292S++5kVd1bVfsy9GifAvzCRi5pVkmeB9xSVRds7Frm4cCq2o/h1zRfm+QZozOX6H6yObAf8N6qejJwN8PQw/2WaN0AtPHiFwCfHp+31Opu4/GHM7xp7gpsAxyyUYuaxiYX4Em2YAjvj1fVKa355iTL2/zlDL1cmPnnANbVvvs07RNTVbcDX2L4iLVDkqmTrUa3dX99bf72wPdnqXshfvbgacALkqxm+DXKXwHe1UHdUz0squoW4FSGN82lvp9cB1xXVee363/LEOhLve4phwIXVtXN7fpSrvtZwNVVtaaqfgacwrC/L619e5LjWxv7wjDm9FHgL8fa/xdrHyz5szb9XNY+WPL11r4Tw1jjju1yNbBTmzd+sOSwCdS9DNihTT8cOA94HkNPZfSAyWva9GtZ+4DJp9r0E1j7gMlVDAdLNm/Te/HAAZMnTPixP4gHDmIu6boZelPbjUx/laF3taT3k7be84C92/RxreYlX3db9yeAl/fwumQ4BnUJwzGpMBw8fv1S27cnGqAb+wIcyPAx7J+Bi9rlMIaxqLOBy4H/P/Kkh+GfUlwJfAtYObKu/wRc0S6jO91K4OJ2m79i7IDSPOv+98A/tbovBv6otT+m7ZhXtB1nq9a+dbt+RZv/mJF1vbXVdhkjR+Lb4/DdNu+tC/DYH8QDAb6k6271fbNdLpla71LfT9p69wVWtX3l7xiCrIe6t2HokW4/0rak6wbeDnynrfdvGEJ4Se3bnkovSZ3a5MbAJemhwgCXpE4Z4JLUKQNckjplgEtSpwxwSeqUAS5JnfpXgndguyc8aW4AAAAASUVORK5CYII=\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": [],
            "needs_background": "light"
          }
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAEICAYAAACpqsStAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAWyUlEQVR4nO3debhsVX3m8e8rICTiwHCDQISLQ+wQ06J9Q9QWg61RcAia2ApxQKNN7IR0ErUNzhCTNDFRn6RNYmMkICDiRIsBFWIIoFH0YkBRRAZBJi8XCINDqxd+/cdeJcXxzKfuOWddvp/nOc+pWnufXb9ateuttVfVPpWqQpLUp/usdAGSpMUzxCWpY4a4JHXMEJekjhniktQxQ1ySOmaIa0ZJXp/k7yewnSOTnDiJmiYhyX5JLl3pOlaDJO9O8qYJbWuPJN9JslW7/i9JXjGJbbftfSLJoZPa3pbCEJ9Dktcl+cSUtstmaDt4jm3tn+TaCde3f5K72pNn/OfxS912Vf1ZVU3sSbhUSdYmqSRbz7LOkUl+lOSO9vONJO9Ksutonao6r6oeuTxV36O2iT/+c9zeVUm+3/rh1iT/muSVSX78vK+qV1bVW+e5rafOtk5Vfauqtq+qOydQ+0+88FfVgVV1/FK3vaUxxOd2LvCEsdHFrsA2wGOmtD28rbvZzBJe17cnz/jP5zZnLavcKVV1f2BH4LnAg4ELxoO8R7O9eM3i2a0v9gSOBv4IeO9EC2PRtWkCDPG5fZEhtPdp1/cDzgYundJ2RVVdn+RlSS5po58rk/w2QJL7AZ8AdhsbLe+W5D5JjkhyRZKbk3wwyY7tb0Yjz5cn+RbwzwstPsleSc5p9ZzVRqUntmU/MTIcH3GNj4baoezhU9a9KMmvt8t/leSaJLcnuSDJfrPU9Lg2Kry1bWP/sWX/kuStST7baj4zyc5t8ehF8tb5HG1U1Y+q6qvAC4CNwKunu99j/X9Hkq8lee7Yspe2Wt7Z6r0yyRNa+zVJbhw/xE+ybZK/TPKtJBsyTFf81KQe/yTbJTmxrXtrki8m2WW2fmh9cVtVndb64tAkj2q3cVySP2mXd07yj227tyQ5r9V3ArAH8PFW92tnqG26I6WHJflC2y8+NnbfZtz3khwAvB54Qbu9i8b2jVe0y/dJ8sYkV7fH4H1JHjil3w5tj8NNSd4wVx/1yhCfQ1X9EDgfeFJrehJwHvCZKW2jgLkReBbwAOBlwDuTPLaqvgscyD1HzdcDvwc8B/gVYDfg34G/mVLGrwA/Dzx9EXfh/cAFwM7AW4HFzimeDBwyupJkb4bR3emt6YsML2o7ttv8UJLtpm4kye7tb/6krfsa4CNJ1oyt9psMffczwH3bOnB3fz9oIUcb7fD+YwwvttO5oi17IHAUcGLuOWr/ZeDLwE7tvn0A+CWGo68XAe9Ksn1b92jg5xj64uHA7sCbJ/j4H9rqfEir55XA9+fTD60vvgBcO0NfvLotWwPswhCkVVUvBr7FMKrfvqreNkNt03kJ8FvArsAm4K/nUeMngT9jOKLavqoePc1qL20/TwYeCmwPvGvKOk8EHgk8BXhzkp+f67Z7ZIjPzzncHSD7MYT4eVPazgGoqtOr6ooanAOcyczhAcOT8A1VdW1V/QA4EnjelNHMkVX13aqa6cm6Wxs9jf/cL8keDGHzpqr6QVWdC3x8wfd+cCqwT5I92/UXAh9tNVNVJ1bVzVW1qareDmzL8ASa6kXAGVV1RlXdVVVnAeuBZ4yt8w9V9Y12fz/I3Uc8S3E9w4vGT6iqD1XV9a2eU4DLgH3HVvlmVf1DezE4hSFA/7j16ZnAD4GHJwlwGPCHVXVLVd3BEEazvVey0Mf/Rwzh/fCqurOqLqiq2yfUFz9iCNs921HMeTX3P1eaa988oaoubi9ibwKenzYNuUQvBN5RVVdW1XeA1wEHT+m3o6rq+1V1EXARMN2LQfcM8fk5F3hiOxRcU1WXAf/KMFe+I/Cotg5JDkzy+XY4eitDOO0804YZRrOnjsIXuAS4k2EkNHLNHPVdX1UPmvLzXdrIrl0euXr+d/tuLZBO5+5AOgQ4abQ8yWsyTCPd1u7HA5n+fu8J/NfxFxyGEdP4yPfbY5e/xzDKWqrdgVumW5DkJUkuHKvnUVNq3zB2+fsAVTW1bXuGEexPM8y/j7b1ydY+k4U+/icAnwI+kOT6JG9Lss0s25/OTH3xF8DlwJlt2uiIeWxrrn1zfPnVDFOTsz0f5ms37rkvXw1szT37bXPsR6uOIT4/n2MIpf8GfBagjX6ub23XV9U3k2wLfAT4S2CXqnoQcAaQtp3pRjXXAAdOCeDtquq6sXUW+68mbwB2aPOxI3uMXf4uQ+gA0EZIswXOycAhGeait2N4b4AM89+vBZ4P7NDu923cfb/HXcMwOhu/v/erqqPncX8W1Q8ZPo3xbIajp6nL9gTeAxwO7NRqv3iG2udyE0Og/8LYfXtgVY3CY8mPfxshH1VVewNPYJi6e8l8C0zySwwh/pmpy6rqjqp6dVU9FPg14FVJnjJL7bO1jzxk7PIeDKP9m5h735tru9czvACOb3sT93zBvVcwxOehHSquB17FPYPgM61tNB9+X4ZphI3ApiQHAk8bW38DsNPoDZjm3cCfjqYpkqxJctCE6r661X1UkvsmeSJDmI18A9guyTPbaO6Nrf6ZnMHwxPljhvnKu1r7/RmeQBuBrZO8meE9gemcCDw7ydOTbNXeqNs/yc/O4y5tBO5imAOdU5Kt2zzoyQyfUHnHNKvdjyEwNra/eRnDSHzBWn+8h+F9kJ9p29s9yWi+eMmPf5InJ/nFFnq3M4TiXTOtP/Z3D0jyLIb5/BOr6ivTrPOsJKNpodsYjghG297APPt9ihcl2TvJTzPsNx9u01Jz7XsbgLUZ+zjkFCcDf5jhjfvtuXsOfdMiauyaIT5/5zC80TY+gjmvtZ0LP55y+B8M87j/zvAG3Wmjlavq6ww735Xt8Hk34K/aOmcmuQP4PMMbaQsx/omH0c9vtGW/2bZ3C/AW4H1j9dwG/A7w98B1DKOjGT/H3OZsPwo8leENvpFPMUwbfIPhsPb/McNhdlVdAxzE8KbZxrbe/2Qe+2JVfQ/4U+Czrf8eN8OqL0jyHYYgOg24GfhP7Y3Eqdv8GvB2hqOtDcAv0o62FumPGKYkPp/kduCfaO8NTOjxfzDwYYYAv4RhvzxhlvU/3rZ7DfAGhheyl82w7iNavd9h6I+/raqz27L/Bbyx1f2aGf5+OicAxzFMbWzH8PyYz773ofb75iRfmma7x7Ztnwt8k2Gf+70F1LXFyNzvW2hLkuRIhjfFXrTStUhaOkfiktQxQ1ySOuZ0iiR1zJG4JHVsWf9pzc4771xr165dzpuUpO5dcMEFN1XVtOdwLGuIr127lvXr1y/nTUpS95LMeKa10ymS1DFDXJI6ZohLUscMcUnqmCEuSR0zxCWpY4a4JHXMEJekjhniktSxZT1jU9KWZ+0Rp8+47Kqjn7mMldw7ORKXpI4Z4pLUMUNckjpmiEtSxwxxSeqYIS5JHTPEJaljhrgkdcwQl6SOGeKS1DFDXJI6ZohLUscMcUnqmCEuSR2bM8STPCTJ2Um+luSrSX6/te+Y5Kwkl7XfO2z+ciVJ4+YzEt8EvLqq9gYeB/xukr2BI4BPV9UjgE+365KkZTRniFfVDVX1pXb5DuASYHfgIOD4ttrxwHM2V5GSpOktaE48yVrgMcD5wC5VdUNb9G1gl4lWJkma07xDPMn2wEeAP6iq28eXVVUBNcPfHZZkfZL1GzduXFKxkqR7mleIJ9mGIcBPqqqPtuYNSXZty3cFbpzub6vqmKpaV1Xr1qxZM4maJUnNfD6dEuC9wCVV9Y6xRacBh7bLhwIfm3x5kqTZzOfb7v8z8GLgK0kubG2vB44GPpjk5cDVwPM3T4mSpJnMGeJV9RkgMyx+ymTLkSQthGdsSlLHDHFJ6pghLkkdM8QlqWOGuCR1zBCXpI4Z4pLUMUNckjpmiEtSxwxxSeqYIS5JHTPEJaljhrgkdcwQl6SOGeKS1DFDXJI6ZohLUscMcUnqmCEuSR0zxCWpY4a4JHXMEJekjhniktSxrVe6AEkLs/aI06dtv+roZy5zJVoNHIlLUscMcUnqmCEuSR0zxCWpY4a4JHXMEJekjhniktQxQ1ySOmaIS1LHDHFJ6pghLkkdM8QlqWOGuCR1zBCXpI7NGeJJjk1yY5KLx9qOTHJdkgvbzzM2b5mSpOnMZyR+HHDANO3vrKp92s8Zky1LkjQfc4Z4VZ0L3LIMtUiSFmgpc+KHJ/lym27ZYWIVSZLmbbEh/nfAw4B9gBuAt8+0YpLDkqxPsn7jxo2LvDlJ0nQWFeJVtaGq7qyqu4D3APvOsu4xVbWuqtatWbNmsXVKkqaxqBBPsuvY1ecCF8+0riRp85nz2+6TnAzsD+yc5FrgLcD+SfYBCrgK+O3NWKMkaQZzhnhVHTJN83s3Qy2SpAXyjE1J6pghLkkdM8QlqWOGuCR1zBCXpI4Z4pLUMUNckjpmiEtSxwxxSeqYIS5JHTPEJaljhrgkdcwQl6SOzflfDCWtjLVHnL7SJagDjsQlqWOGuCR1zBCXpI4Z4pLUMUNckjpmiEtSxwxxSeqYIS5JHTPEJaljnrGpnzDbmYJXHf3MZaxEvZtpX3I/mhxH4pLUMUNckjpmiEtSxwxxSeqYIS5JHTPEJaljhrgkdcwQl6SOebLPFmShX+flCRf3bp6IM3krcaKcI3FJ6pghLkkdM8QlqWOGuCR1zBCXpI4Z4pLUsTlDPMmxSW5McvFY245JzkpyWfu9w+YtU5I0nfmMxI8DDpjSdgTw6ap6BPDpdl2StMzmDPGqOhe4ZUrzQcDx7fLxwHMmXJckaR4WOye+S1Xd0C5/G9hlphWTHJZkfZL1GzduXOTNSZKms+Q3NquqgJpl+TFVta6q1q1Zs2apNydJGrPYEN+QZFeA9vvGyZUkSZqvxYb4acCh7fKhwMcmU44kaSHm8xHDk4HPAY9Mcm2SlwNHA7+a5DLgqe26JGmZzfmvaKvqkBkWPWXCtUiSFsgzNiWpY4a4JHXMEJekjvn1bPdiC/06N0mrjyNxSeqYIS5JHTPEJaljhrgkdcwQl6SOGeKS1DFDXJI6ZohLUsc82UfSPcx0EthVRz9zmSvRfDgSl6SOGeKS1DFDXJI6ZohLUscMcUnqmCEuSR0zxCWpY4a4JHXMEJekjhniktQxQ1ySOmaIS1LHDHFJ6pghLkkdM8QlqWOGuCR1zBCXpI4Z4pLUMb+eTVphM30d2qS2sxq/Vq2nWlc7R+KS1DFDXJI6ZohLUscMcUnqmCEuSR0zxCWpY0v6iGGSq4A7gDuBTVW1bhJFSZLmZxKfE39yVd00ge1IkhbI6RRJ6thSR+IFnJmkgP9TVcdMXSHJYcBhAHvssccSb04r7d52pt1sZ1Mu9D5P6szMhdrcZ4RuyXrY35c6En9iVT0WOBD43SRPmrpCVR1TVeuqat2aNWuWeHOSpHFLCvGquq79vhE4Fdh3EkVJkuZn0SGe5H5J7j+6DDwNuHhShUmS5raUOfFdgFOTjLbz/qr65ESqkiTNy6JDvKquBB49wVokSQvkRwwlqWOGuCR1zBCXpI5tsV/P1sOH9CXNz+Z+Pvd8IpMjcUnqmCEuSR0zxCWpY4a4JHXMEJekjhniktQxQ1ySOmaIS1LHDHFJ6tgWe8amtLn1fJbfatVLn66mOh2JS1LHDHFJ6pghLkkdM8QlqWOGuCR1zBCXpI4Z4pLUMUNckjrW/ck+C/3Q/WI+pL/Qr4Ba6FdJLbSm1fgVcyv19VmrsS+0fDb3c60HjsQlqWOGuCR1zBCXpI4Z4pLUMUNckjpmiEtSxwxxSeqYIS5JHevmZJ+V/JD+pG57SzzRQNLKciQuSR0zxCWpY4a4JHXMEJekjhniktQxQ1ySOrakEE9yQJJLk1ye5IhJFSVJmp9Fh3iSrYC/AQ4E9gYOSbL3pAqTJM1tKSPxfYHLq+rKqvoh8AHgoMmUJUmaj1TV4v4weR5wQFW9ol1/MfDLVXX4lPUOAw5rVx8JXLr4ciduZ+CmlS5iDj3UCH3UaY2T00OdPdQI86tzz6paM92CzX7afVUdAxyzuW9nMZKsr6p1K13HbHqoEfqo0xonp4c6e6gRll7nUqZTrgMeMnb9Z1ubJGmZLCXEvwg8IsleSe4LHAycNpmyJEnzsejplKralORw4FPAVsCxVfXViVW2PFblNM8UPdQIfdRpjZPTQ5091AhLrHPRb2xKklaeZ2xKUscMcUnq2BYb4kmOTXJjkovH2o5Mcl2SC9vPM8aWva79+4BLkzx9hes8ZazGq5Jc2NrXJvn+2LJ3L1OND0lydpKvJflqkt9v7TsmOSvJZe33Dq09Sf669eeXkzx2BWv8iyRfb3WcmuRBrX219eWq2TdnqXG17ZfbJflCkotanUe19r2SnN/67JT2wQuSbNuuX96Wr13BGk9qj+fFLQO2ae37J7ltrC/fPOeNVNUW+QM8CXgscPFY25HAa6ZZd2/gImBbYC/gCmCrlapzyvK3A29ul9fOtN5mrnFX4LHt8v2Bb7Q+extwRGs/AvjzdvkZwCeAAI8Dzl/BGp8GbN3a/3ysxtXWl6tm35ypxlW4XwbYvl3eBji/7W8fBA5u7e8G/nu7/DvAu9vlg4FTVrDGZ7RlAU4eq3F/4B8Xchtb7Ei8qs4Fbpnn6gcBH6iqH1TVN4HLGf6twGY3W51JAjyf4UFeMVV1Q1V9qV2+A7gE2J2h345vqx0PPKddPgh4Xw0+Dzwoya4rUWNVnVlVm9pqn2c4n2HFzNKXM1n2fXOuGlfRfllV9Z12dZv2U8B/AT7c2qful6P99cPAU9p9WfYaq+qMtqyAL7CE/XKLDfFZHN4OrY8dHf4z7KDXjK1zLbM/sZbLfsCGqrpsrG2vJP+W5Jwk+y13Qe0Q9DEMI4pdquqGtujbwC7t8or255Qax/0WwxHCyGrqS1iF++YMfblq9sskW7VpnRuBsxiOVG4de+Ee768f92Vbfhuw03LXWFXnjy3bBngx8MmxP3l8m375RJJfmGv797YQ/zvgYcA+wA0Mh4Sr2SHcc7RzA7BHVT0GeBXw/iQPWK5ikmwPfAT4g6q6fXxZG1Gs+OdVZ6oxyRuATcBJrWm19eWq2zdnebxXzX5ZVXdW1T4MI9l9gf+wHLe7EFNrTPKoscV/C5xbVee1619i+D8pjwb+N/B/59r+vSrEq2pD69C7gPdw92HpqvsXAkm2Bn4dOGXU1g6pb26XL2AYdfzcMtWzDcMT+qSq+mhr3jCaJmm/b2ztK9KfM9RIkpcCzwJe2F5sVl1frrZ9c5a+XFX75VgNtwJnA49nmL4bncg43l8/7su2/IHAzStQ4wGthrcAaxhe+Ebr3D6afqmqM4Btkuw823bvVSE+ZV72ucDoEyGnAQe3d6/3Ah7BME+1kp4KfL2qrh01JFmT4f+4k+ShDHVeubkLafOG7wUuqap3jC06DTi0XT4U+NhY+0syeBxw29i0y7LWmOQA4LXAr1XV98baV1VfrqZ9c5bHG1bXfrkmd3/a6KeAX2WYvz8beF5bbep+Odpfnwf88+hFfZlr/HqSVwBPBw5pL9yj9R88mqdPsi9DRs/+QrOYd1x7+GE43LsB+BHDvNjLgROArwBfZnhAdx1b/w0MI4hLgQNXss7Wfhzwyinr/gbwVeBChsOuZy9TjU9kmCr5crvtCxneXd8J+DRwGfBPwI5t/TB8YcgVrb/XrWCNlzPMg47aRp9OWG19uWr2zZlqXIX75X8E/q3VeTF3f1rmoQwvdJcDHwK2be3bteuXt+UPXcEaN7XHdNS/o/bDW19exPBG/BPmug1Pu5ekjt2rplMkaUtjiEtSxwxxSeqYIS5JHTPEJaljhrgkdcwQl6SO/X8vD1/ZEOUaPAAAAABJRU5ErkJggg==\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": [],
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "UCwWn6Dq1JTa",
        "colab_type": "code",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 35
        },
        "outputId": "145df736-017c-4322-9ac5-3497036088ad"
      },
      "source": [
        "print(np.mean(water_equivalent_diameters))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "245.78901400991444\n"
          ],
          "name": "stdout"
        }
      ]
    }
  ]
}