Diagnosis/Week3/C1_W3_Lab_1_explore_mri_data_and_labels.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Explore MRI Data and Labels\n",
    "\n",
    "<img src=\"mri-slice.png\" alt=\"U-net Image\" width=\"300\"/>\n",
    "\n",
    "In this week's assignment, you'll be working with 3D MRI brain scans from the public [Medical Segmentation Decathlon](https://decathlon-10.grand-challenge.org/) challenge project. This is an incredibly rich dataset that provides you with labels associated with each point (voxel) inside a 3D representation of a patient's brain. Ultimately, in this week's assignment, you will train a neural network to make three-dimensional spatial segmentation predictions for common brain disorders. \n",
    "\n",
    "In this notebook, you're all set up to explore this exciting dataset. Run the code below and tweak it to explore further!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Import packages\n",
    "For this lab, you'll import some of the packages you've seen before (`numpy`, `matplotlib` and `seaborn`) as well as some new ones for reading (`nibabel`) and visualizing (`itk`, `itkwidgets`, `ipywidgets`) the data. Run the next cell to import these packages."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Import all the necessary packages\n",
    "import numpy as np\n",
    "import nibabel as nib\n",
    "import itk\n",
    "import itkwidgets\n",
    "from ipywidgets import interact, interactive, IntSlider, ToggleButtons\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "import seaborn as sns\n",
    "sns.set_style('darkgrid')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. Explore the Data\n",
    "\n",
    "### 1.1 Loading Images of the Brain\n",
    "\n",
    "Run the next cell to grab a single 3D MRI brain scan"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Type of the image <class 'nibabel.nifti1.Nifti1Image'>\n"
     ]
    }
   ],
   "source": [
    "# Define the image path and load the data\n",
    "image_path = \"data/BraTS-Data/imagesTr/BRATS_001.nii.gz\"\n",
    "image_obj = nib.load(image_path)\n",
    "print(f'Type of the image {type(image_obj)}')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1.2 Extract the Data as a Numpy Array\n",
    "Run the next cell to extract the data using the `get_fdata()` method of the image object"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "numpy.ndarray"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Extract data as numpy ndarray\n",
    "image_data = image_obj.get_fdata()\n",
    "type(image_data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The image object has the following dimensions: height: 240, width:240, depth:155, channels:4\n"
     ]
    }
   ],
   "source": [
    "# Get the image shape and print it out\n",
    "height, width, depth, channels = image_data.shape\n",
    "print(f\"The image object has the following dimensions: height: {height}, width:{width}, depth:{depth}, channels:{channels}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As you can see these \"image objects\" are actually 4 dimensional! With the exploratory steps below you'll get a better sense of exactly what each of these dimensions represents.\n",
    "\n",
    "### 1.3 Visualize the Data\n",
    "The \"depth\" listed above indicates that there are 155 layers (slices through the brain) in every image object. To visualize a single layer, run the cell below. Note that if the layer is one of the first or the last (`i` near 0 or 154), you won't find much information and the screen will be dark. Run this cell multiple times to look at different layers. \n",
    "\n",
    "The code is set up to grab a random layer but you can select a specific layer by choosing a value for `i` from 0 to 154. You can also change which channel you're looking at by changing the `channel` variable.\n",
    "\n",
    "Keep in mind that you could just as easily look at slices of this image object along the height or width dimensions. If you wish to do so, just shift `i` to a different dimension in the `plt.imshow()` command below. Which slice direction looks the most interesting to you?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Plotting Layer 85 Channel 0 of Image\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOcAAADnCAYAAADl9EEgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO19WYxl51X1OvfeuvNUc1V3u40dx9iOnJhYKEQWGAGKZFlxGKNIoACKkBBvvCDlBYlXJMRbhOA9UniAB4KsIGygbSdBiURIHMft2OCkU13jnef5fyjWV+vsOtXxn7Ttc6u/JV1V9x3O+K1v7732/vYJFosFPDw84ofE+30AHh4e0fDk9PCIKTw5PTxiCk9OD4+YwpPTwyOmSN3pwyAIvJTr4fEuY7FYBFHve8vp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjFF6v0+gHsFyWQSABAEQej9xWJxx78/DbivIAhC++W/p9PpT70Pj3cPnpzvET7/+c9jMplgOp26V6/XQ6vVQrPZRKPRQLfbxXg8xnA4xGKxcK/5fI7ZbOZ+Nx6PQ9teWVlx/14sFkgkEkgkzpyidDqNXC6H8XiMlZUVpFKnt/327dvvzcl7/EQI7jRDB0Hw00/flxh/+Zd/iWw2i0wmA+CUGLRSQRA460dyzedz99t0Oo3BYIB6vY7j42McHBzg6OgIuVwOH/rQh3Dt2jVcvXoVtVoNtVoNzWYTrVYLvV4PnU4H7XYbb7/9NkajETKZDFZWVpDNZpHNZkP7nM1mmEwm7lh4bPP5HPP5HNPp1BFfj28wGKDf779Xl/KexmKxCKLe95bzJ8Rf//VfOws0Ho+RTCYRBEHIai0WCzfwR6MRJpMJgFMXd7FYYDAYoN1uo9lsotlsotfrIZfLIZ1OY7FYoNVqYTabIZfLIZlMYnV1FcPh0H0fAGazGVZXV1GtVlGtVlGpVDCfz913SOpsNot8Po9CoYBisYjZbIbxeIzRaITRaIThcIhUKuXOYTgc4rXXXsN4PMZsNkOn03l/LvQ9DE/Od4C/+qu/QiqVCsWDQRBgMpk4C5RKpRzp5vO5G/D9fh/dbhe9Xg8AkM1msbKygnQ6jWaziePjYzQaDfR6PaRSKRQKBaRSKcxmM/R6PYzHY+fO0rKVSiVks1msra0hmUyiUCggn88jn88jl8thMpmgUCigUCigVCphfX3d/Zuv6XSKyWTiXtPpFCsrK0gmk0gmk5jNZvjwhz+M4+Nj9Ho9/PM//zNms9n7dg/uRXhy3gFf+MIXMJvN3GDVOJD/n8/nWCwWjkjdbhfNZhP1eh3tdhudTgedTgeDwQDpdBrpdBqJRALz+Rz9fh+LxQIrKysoFArO+mWzWSQSidC+1PVMJBLIZrMolUoIgsCJTbTOk8nExaXpdBrA6aRAS68WPAgCZ6mTyaSz/IvFAul0GuVyGdPpFN/61rfc73jOnJjosnvrenfhySn4/Oc/j93dXReT0WJxUAIIxZEcpHT7jo+PUavVXPw4HA4xGo1cXKeK6Ww2QyKRQC6XQ6lUQjKZdG6lTgJ0MzWWTSQSzsKp8KOTBoDQ97hdxqB2PySbKrtK2D/4gz9w12QwGGA6naLf72MwGKDX66Hb7eLGjRveut5FeHIC+KM/+iNsb2+jXC6HUhgcaFGqKd8HTlMSFGlqtRpOTk7QaDQwnU7dd1ZWVkLCEUmUSCQwnU7RbreRSCSQSqWwubmJRCKBfD7vLJkSl6TjZyruqDVLJpNOLKLbrUJVlCIMICQMKfkBIJVKoVqtIggCpxyPRiOMx2PcunULw+EQvV4Pk8kE3W733bhd9wzuSXL+8R//MXZ3dx1BKLgAcFYFgBvIHPB8KUg0xpe9Xg/9fh+j0cgpqLlcDsViEfl83lmeTqeD0WiE+XyO4XDofj+dTrG2tuas2GKxcJaX+ychGZsqIYEzcgVB4FxZtYh6Duqi6l9Vnmm5AYQmBZK/VCohkUjgc5/7HLrdLr773e+iXq/jxo0b7+JdvPy458j5uc99Do8++ug50nGQk7A2cU9wUCoRaEmTySSy2Syq1aoTaCqVCsrlMqrVKmazGVqtFmq1Gkaj0TnLxZi11Wo5i0drx9hSXVwSmMelx0gCWhdZvQH9jp10COvu0iIrcfkKggDlchn3338/dnd3cXh46MICWtnpdIogCNBsNu9KocVlxj1Dzt/+7d/G448/7lINdjDa3B+Ac/GTWhC1MtPpFMlkEsViEalUChsbGy6ezOVyKBQKyGazaDabGI/HLjequUhuk/vhvlRB5XvWjdXj021EWXmeV5TF1FQKv0/PgN8jQXksdJf1/9euXUMqlcKHP/xhtNtt1Ot11Ot1NBoNNBoNzOdzvPzyy6FY3uM87gly/sVf/AWKxaJLdTAWVKLZOEwFH30pOPiZ56SFTKfTyGQy5wScTqcTssY8hnQ6jSAIkM1msbm5iZ2dHVSrVZTLZeeW2uOMek9FIyA8mSgJ+V1+h381R8tYWKuaKGrpd0lOEpTCVjKZdNZS1eVUKoVer4dqtYrBYODj0jvgUlYI/dmf/RlWV1dRKBTOlb4pMXUwqvLKgcg8pgowdC+ZrphMJhiNRi7Jn8vlXNqCSiZVW1b3sCpoPB5jY2MD29vb2NnZwc7ODnK5nCMkJxKt4omqh40io8JORHq+OhmpFVRyKkF1m/pvbiORSDgBioKWdZ+ZXmq32/j2t7+Nfr+PyWRyz7q690yF0J//+Z+jWq1isVg4wUUH9mQyQSqVQiqVOleDSvfSWiCN+WgZ1I1LpVJuQNLtYx6z2+2G8p08hnQ6jWw2i2KxiGw2i3Q67WI5TYfw+1RFKViRBOoCR8WNSkwgHJuqW615W4pjnKC09FDjct0O96WT2HQ6ddeH14yFFtVqFcPhEPfddx/q9Tp6vR6+8pWvYDQa3ZVxcBlwqcj5C7/wCyiXy6GZ3oo+VuzR1IJC3UW6c/qXIommORgHMt1hrQtjLLp5mUwG+Xwe2WzWEcKqwmqhaUFJMBurRh1/FGFpzay7rp6FuvYWWpCh+9Ptk5xq1XncnFSy2awrPwSAvb09R071dG7dunVP5k8vjVv7xBNP4Dd+4zeQyWTOkdMqlplMJpTAj1Ju1WoQKpYAZ4PwTkuyaPEGgwEajUbIUqXTaVdOVy6XUSqVsLKycq4YXcUYFYQ0frSurZ4T3VG1cNb9tYSli2q3bUsJ9fpxUqKl5HVOp9OuZNGuoFksFiHvI5lMhryObreLv/3bv3Xlj5cRl9qtfeyxx/Dss8+iVCo56xSlqhIq0nCGB87SIurGqQVkDMj3+Uqn0+fi2ul06oicTCaRy+XcoCbhV1ZWXMqF1gRAyAXk8aprTcK32223H3UdtcbX5i1tpdJFsBVJWgChk5R6EOr68jytB8FtW7eY1pgeBIsdKpUKHn/8cXzzm9+859afLj05H3roIXzyk59EpVIJDSZN1kdZCc7UJNh4PHaCUNTgs9BtZDKZc9YJgCOTDkZdvsV1lhSRKCRpHvEiqz4cDtFqtTAcDpHJZFwMSxdZ3chUKuUmHv7bFhlY4YbHz3MlSCK7BM5+j9BjVjfXusW2DDGZTDpr+9nPfhbD4RDf+c537in3dunJ+Zu/+ZtYX18/F7NFxUQckFHpEVVigbNcoIpEmmawBODg13I2kpfF7nTbeBwkE18UhdRSc/H1YDBw/1YBx7q5PHe1oiSeutTWJbaiUdT7quqqh6BFHPw9vQyd3GwMaxVkboMrZLLZrCPqM888AwB49dVX7xkLutTk/NM//VNsbW2du/nqyvLfWg9LV1AHEddcqihi400g7M5aS2FL3bhv5vqohDIOS6VSGA6HaLfb6Pf7GA6Hbv0kicdUzXA4dGmNSqWCarWKq1evIpfLuX1ZS0hXN51OR1ZD2etkVWg9b/57Npu5uNiKVTYkINRi6jXUberExtSUimulUgnPPPMM5vM5XnvttXuCoEtLzvvvv9/VoKoVjFJpeZP5XbqwHJBaWjYcDt2N13SLxoG2ZI2zu3YkoKA0Ho/doFssFi4fysHZ6/VwfHyMo6MjHB8fYzAYoFAouOVltOSEFiWwrjXK/VU1WuNEuu1WkeVvbIypbrhujy8KOFGxrJ0kNC62E5haYv0tPQPGr7/6q7+K+XyO119//dITdGnJ+dxzz4WEmPl87iyLum7qSvEvrRAre9Tl01nfkpADRcnA71nLRaGJlrPX64VK9RhTJZNJ5PN5rK+vO/f3+vXrIfXWxp8kCwe7Cir8Ny0wj5FKqS4h0/MMggCNRgODwQCDwSC0NlSvcbFYxM7Ojls0rmGBFbFUoeWLk43GrZzc8vm8K4HU41VRKZFI4Jd/+Zfxgx/84NKvH106cj733HN48MEHsba2BuAsH0bLNxqNQjedsJZC1U+Nr0gGICxSaNcDik0Ub2x6QEGrMJlMUK/XQ5Yhm82iXC6jUChgc3MTV69eRTqddrlPur53UlfVupGU/L+NvUmCKIsMwLm/vJ48PiVYlAhk3Xxey3w+j0wmc67+djKZOBeeLn8mk3G5XxJcJxpOGlzBcy9UEi0VOcvlMu67776QRVGXjWTjTKsvuwjZWkSbzwPOr2vUPkFqJfldSwhaZ10SRreZlT5bW1vY3t521UL5fN6RRwe0QvfFiYYEYgyqx6hiGcFroce7uroaui4Kniu7BqpV1WtED0Kts03r8Dd03dkChovW6cKqS8xtMAb35IwR1tbW8OlPfxobGxuhVh/WdQXO94YFztIJ+pnm6Pg7OzB1sNFNY6E63dTBYBCy1FaYURFIjzsIAhQKhVAqJZvNukFoW2By23bCsZZLXUe1nFFxqFpQJax9z3oeuh16EgpOTgBcaSMtonobPA/1fpjy4fXS7bGgw5MzRvjd3/1d7OzsADgdPFRWVSmMSpHoXx10Nl5UIqkFsDGrdkKwKQSSkC8utNbmWjbW4+eFQsHlXW1lkB6/VaGjCgOiCgLsxGPd0ajJzuY9+e+oa2q9iPl87rwEVkMtFgtHSrqrXHxu749uk/eDsXuj0XCrftiH6TJiaci5tbWFdDrtbijjFRtfqnWwg0hznTZFYK0l3+eA1pK+qN/agc74NZPJuDWkVjXW42V+k5OOuq56Dnyf+6Sgxe1HfV/rcS3U/VXSAzhHFn7fur4aFijJdUULY0seh21aRjeYai7ddFpczdcOh0Nsbm5iOp3i5OQEs9nsUpJ0acgJhMvCOCPrQLWDww42G1dai0NYsUiVRe5b96sxk6YU6L4xvULroaqqLstSMcsSQycaHp8WvkcJQPyeCiz6FzjfJ8lOPFHHwGti1Wl+314nEkqPTUHy2XSL1j8zni4Wi6hUKo6QOzs7mM1mODk5AQC0Wq07jqFlwtKQs9/vu9lYyaGLfDk4aOFUAALOBqWu+FcLpP+21kEtlF2VcdGg09iQJFSoNVaBxYpcGgtakUWvg80hEroNPTat7lH33OY1+dcSMSq2t8q39QB0H9ZD0Yoqva+8h+l0GsVi0VnaVquFdruNwWDgFtN7cr4PoNsSlefT2ZaDjivx7RIoAOdcKLWUHBB3cu9IVOBM0Ihy8fS3zHdaV9XmY3kcqtRasti4Okrc4fa5XX6m1khdS04eNo62QpOeF49XU1OKKHHNHhcnF61TZn3wRefKOmJ2qZ/NZqGO+pcFS0NOqncK3lRdbkQLpOKITanQcloiXRTjvRNwkFkLxcFLaHxov6vuNfcfpSTbV5QAptvUz7UQwoYIer66Xz2mKOHG5or1+3od9DOdJLi8jC+myWwIYIsqyuUystks1tfX3Tm8/PLLUbdnKbE05OQKDuD8IGHZnRaGA+EBQGi1CbdlUwIao3LgWdVWcVGcFhUL831uU4/FKqxRUFdQl5lZWEHMpkLU3dbBHxWzWtiw4KKJSY9X41lLZOY6AYS0BHXneU9TqZSL3zc2NlwaqlAoIJFIeHK+1/jsZz+L8XiMRqMRyhsCcOqtJrbVFQPCSXtdEmZdWW1ORQGCFlkFITvYuQ/rJvI9VS0VdpCqlbOijVoOnVAugirIPH5eJ1umyGvA/r3WrbaxuV5fei0aK+ri8vF4HJnCUvAe6rVMJBKhvHKxWEQmkwktTtC+wFwTS3f3MmApyPk///M/ePnll51ooEQqlUqoVqtYW1vDxsYG8vl8iCiEElGtgY2b+JdVPbR4fIYmRQhdEsZBUSwWXVFBLpdzv9VlYmoNLGx8Za2O7cSnsEKLtZQXpW9UJbXrSW0RBYl3+/Zt9/QyPuelWCyiVCqhUCiEXHkq2XzZCUgnTr0+apV18QDjdk7OXCHDie/JJ5/EN7/5zUvRi2gpyJnJZFAsFpFIJFxMQuGgVCo5UiSTSQyHw1BcB4S7lFuSamylau9gMECtVnPuMjvnNZtN1Go1tNvtkBXL5/Ouk97W1ha2trbcszJtmsK61ZZEVinW3+hvo3KU1vJSILO/s2IaJxobyzLFo/XB3/3ud1Gr1dBoNBAEgVvszZhxdXUVq6ureOCBB1AoFNw1tfG8nrfNm2rF0GAwQLPZdLEpizb4uEQtbHjqqadw8+ZNT853G7//+7+ParXqamk5qNS90sJ3zqx3itd0ANBq6OcEWzXyuZlsb8mbTuIRXJt5eHiIdruNH/3oR26VBV+VSiVSdFJ3ldZFCWotqrp/UTEhcFYlZHOS6p7aJL+qvSQIe8vyWrTbbZycnLiqHpJXrSB77+oEFRXXqpVUBT0IAtftgQvNOQFWKhV3Pa2Vn81maDabl6ZbQqzJWS6Xsb6+7nrtaJ5Ru9FpnlAJqNbR/jsqca6z+Xg8Dj1Fi1I9Fz5r8bXND3I7pVIJu7u72N3dxbVr17CzsxOykBysdHuVRAodyPZ9/rXHr8UFdK01F8pY2rYJ1UliMpmg3+87b+Ho6Mi5sxpDa0WUVi11Op2Q2BOVf7WgFU0mk85dTiQSLo+p3o2mpubz05LOg4MDZ8GX3XrGlpy/9Vu/hVQq5TqCs+FwVBEAcdHA5t+L0gJWqp/P56ElTSpCaEXPRfvm7M9i++FwiJOTk3MuJAez5gltwYSSxZ5zVEFA1HZIkKjqG5tGsnEplXA+VqHdbrtzj0pJaZzY6/VcCMJjUPLy+KK6J+jidXbx07Wl2omBHtNwOMTR0REqlUqoM8OyIrbkLJVK54rMo0SPKBVRrSFBIUKT8lpsoIqojYl0n8CZS2wnByvmsHCi3W4jm83i8PDQPbJhdXXV9T6yFtcetxJUl75Zl1cJqq4ii8QvimNtPlXBtZT9fh/9ft8tUI8ipf29Eo/noEpvlMdB8N73+/3Q8ej+GM6QtEEQuC6HQRCg3W57ct5tfPKTn0ShUDhnJX5cTo/ih75vf8sBGlVNolbHiiYkhl3uxL/W8vA9FuczLgPgGlcxZlI31FoQW0xhlV8elx4Pz0H7x+oxaWxrq6Jobdh6s9PpoNvtOnWaD2DS/dnrHfWeFaqslebx6DWwHpJ2cdBKIe24sLm56cKRg4ODyGV3y4JYkpPJZVtap/k3LQ7Q/J2WoEVBrRAQLkLQuJRpGsY+tB7dbjdECOC8Veb2dfkZ41VtYjUYDEITDmNBHbhKVmuldKUGEB7MjCfpPvNc9a+1fLSSfAwhnw7WarXcmlVbcWSVZz1mvaa8T2qleQx2TalNt9h7xN9TIeY62Ewmg1wu56z1Sy+95J7DEhWzxx2xJCebYEUVtauF1EHGogEqfVHCC+MQxoPcNnC+iGCxWLjEdj6fR6/XC1kX7QIQRXh9n5aB8Vuv10Oz2QQAF0/Zgm+Cv7X5Sn1fJzAVZzh4WSyu7Uxs/LpYLDAYDFCv13F4eIjbt2/j6OgotISNUAHLnqtuXy2zfocE1e6AQLgdqYYxUX+t6Mfry1ztysoKPvaxj+FrX/saDg4OlrIZWOzI+dxzz6FSqYRm3Mlk4qR7W63CgaAKrqYKlHwctHZRNPejL67bpFXmdylQUMFl2w6bSySBVCFlQfd0OnXtOOiOMe3CwaUTj7qRGk/z3LgfTY3M53MMBgOcnJzghz/8IebzeUiMYtdB9ubp9XrY29tzrVSs8q2IskI6kTIHqTG9Wr1kMukmDpZlqtBmJ17gLAZV8vL4GQcnEglUKhVks1kEQYBPfOITmE6n+Nd//delbAYWO3Kur68DgFtZwhl2bW0t5NKpe6MuUFSNaFSKQUnNbSh0MM3nc2fh2PRqPB67J4cxzWJzloyJ2LCLDa/4qEBaZhKen9Ot09UuQNiK8hzUheQ5TSYTtFotHBwc4OjoCK+99lpoW9y2Xhfdrt0n/14UU6oFZKkd85vcr+ZbrX6g71PB1lCGqRjG71p/q5MprzmRSqXw1FNP4Rvf+IYn590ALy5nUJX/gTOi2UCfN1RJx3/roNLZXF0n4PyaRe1Va13kyWTixBy6VEoU7oMFDFRNbWECn5Wi6QVtaq3xp9azAuEJhQOUFqXT6bhudTZGtsKQnr/ChhI2VlUryBpXTjpqOaN+q8fN86bHoL2GtECCYhS3pQUofOk27fktG2JHTt48jT8oogDnB4r+TmFVSZvb4z44o+tvNK7T1RJ0z7gv1tNWq9VzeTgeP4sVFotFqF8uhZd8Po9CoeAeXbhYLFzqA4ju9WMxn8+dVeG2O50OWq0WWq0WOp1OSCTiYLeuvFWKo8Qenez0GtE9LZVKroMgQwi1glEFINyO5mftcXDS0RYu/I3qB5pTjbLSy4RYkfPJJ59Ev98P5eY4YH6ahbRWWdWSPZ1hbYog6uZyEHNVRLlcxoMPPuiEHqYeGJNye3yxqoVuMd3AXq/nXmtra+ceTGQHM8UlxmlMwjON0Gq10Gg0MBwOnVprFW9C3X+1Onpt9PrYcyJBOFnRapIgSkqbk7ZkVRfYegaceBhG2P5ClvT0YJYVsSLn9evXUSgUQuQhomY/Hax21tff0PWhO8rKE+su221xclBhSi05Y6TNzU3M53NXnE0xxi5sBs7KzrhfurQcZIwZ2S+nVCoBgCsb7Pf7zlWmpaUrqPErRTUKP7aYIsqdVbFNv0svgNdBiavrLJnKyGaz53KoP24Zmsa2Sir9v/UeOJHwc3pDFAU1zFlGxIqcW1tbqFarjkRElICgkr3O8DoYVPFNJBIuLuH3mcim6BM1iIHz3Qyo+nLws6iAncyZE1Xl0b547LR4wOljIpLJZKhQnqqx7WifSqUcsZk2Ya6Pndvb7TaazaZb4xj1wCFNT1hLSdjFAdpfltdDu1FwsqHyrLW7Sljrxdh8s41TV1ZWzukInGRVGbc1w1FC1jIgNuT8pV/6JaRSKbcCQaGWTGd2xnn8DgeXnbGjhA5ba8r8pyp/6k5qTjSfz7tZmdU02p+1Uqm4IgZ95og+QElFHk4anGh2d3fdc2BIRCbb+RzSi+JQ/j+bzTr3+Fd+5Vcwm81cEQXdZ7rg/Fuv10MPX6JiHnWteE+oXlcqFeRyuRA56KGQQNb95DXQidjGnZoW0/vA/du0GCeIi54Ls0yIDTm57k9h0yF0I3nxOVA4mOxN5oDSgaCf3enGWbVWjylKbaSQk81msbq6CuDUFWW6pdPp4Pj4+Fz9ruY46cZWq1WXlqB6zX1wIgAQSUpeo0wm41w8FvEzXiPxEomEE28qlQoqlYpTeNXLUK8hlUphY2PDrRhRQup9UreVhNIYUUUeXVWkog+V4HQ6jVKp5DwHXUbGFSgUgZjP5T28efPm0j6yPjbk5AN8OCvakiu94eo6WStJgmopn3XlgNMbp60o1S2mCKNxKd1WbldjU50w1FIAcGshG42GOw+teWX/Gyq2LPjXZ1/aNiX8S3dcrxFBiz+dTnFwcIBWq4Ver4fBYHAu0a+/t2KNLmynGvvwww9jdXUV1WoVKysrGA6H6Ha77qXXWC2fLZKgWs5yxna77bbBWl4lJ70UVlglEonQk8E138qUzgsvvIB6vX63h+t7guBOwXIQBO9ZJP2pT30KDz74YMhKcSCzmmQymaDX64UeNMuyPc6eJJ0+0kAtLmd4HejAWYWLyv+qUmoO8/+uTWjZksZHRBAEjniFQgH5fB6NRgNHR0fY39/H/v4+Go2GOx+uE11fX0e1WkWlUkG5XHbuoSq4vC48Vq2Qsg8YAoBOp4OTkxP3ajQabvJSF9u6ngoSvlQqOfGHeVr+mxMbJ0Uth9R92WKMWq2Gg4MDHB4e4vDwEJ1OJ/Q7vUcMP3j/1MPSIgbeQ24rrlgsFpHuW2wspwb6HIgUQvr9Pmq1WmilBOMkxlarq6uoVCpYXV11N4g3VmdwtQhRuTW1RGqxgLMWkOoOay4WOJ+gp6Vi3yE+KPfo6AgnJyfo9XqhHGIymXQ9kcrlsnMdtdKIiwIIWzxBi8rj6na7SCQSrti+2+267vO6Df61+V4eG7FYnFbiqAvMY2faZjweo9VqOYWZL7rWrPnN5XKoVqsuN0vrrvlMe5wa/1vhju/phOrV2p8SunqCbhKJ2W630Wg03OzHOIo3GjhTUPkYOw4yaxmDIHAWR0UOdX8Za/H7dMOYX+TnjH84CHXpkj5ohytSTk5O0G63XWEAa3JtTEZ3nV4BSanPr2TeT61/FDjhZTIZlMtlVydcKBTOKbjqAVivQV/6vBpuj5OG5j0zmYwjLZcA8nONMY+Ojtz1LpVKrvuFTjTcr7Yu0efW8Ls6oWqqaRkRG7e2XC7jM5/5DMrlsrv42jKRuUHGLTqzM61QKBSwtrYWGlS6BpPxC9c32sS+xnOTyQSdTsflNHVQ0GIACMVF2iKSpGjlTCgAACAASURBVO92uy6lcevWLXdejCdXVlZQLBZdvFkul3HlypVQEysrXNn8ox6/WjsgnLvUmGxlZcWtQjk5OUGtVsPt27dDrihdZFuUrjElBbCtrS1sbGxgbW3N9RfmcWkLS3oB7K7QbrdxeHiIZDLpnrjGskaKecBpx/9Wq4WTkxPU63XUarVQXTOvqaZtgFPvp9lserf2pwHbX+ij2IGzGZH5SJvXVMtIUjDPyPiUIg+JyTWTJCO3oa7vaDTCW2+9hVqthlqtFrKmfHEwbW1thdIFlvD6HBQOQuZG6ZZvbW1hfX0d29vbLl2irnNUz9kopVn3zd/Sq+CjCKkMa9zMeJ7iDCeier2O4+Nj92o0Gm7huBZkaAM0XudCoRDyJGxqhZMt7y+vJz0FXmeq1PysWq1iY2PDhTec/Gq1GtbX112YwwUKN27cQK1WexdH77uD2JATOBUttre33Q0jWdh6f2VlJZToZ76OMQoAJwBQIa3Vas66UlBhb1VblsZEOtMXzWYztG3Nm2ninZOCphMIupRs9kX3lYOUA5KWkwu8bRkdt6XKsq2MslU1QPh5JSQVm3bpcScSCddiVC0vQ4parebEpL29PUcMrZrSRzLyOHSyoJupQpou8+N7Nr/MiYDHmMlkXJzKReFsAP7AAw/g6tWr2N7exsHBAer1us9z3g185StfwerqKq5cuRIaOIxRmDfkLEnVkY9BpwtF0SaRSGBnZye0Wl4XcusA0HSMlvUBZ2kXdWuZZ6NLS6Lodjmx0E2kW6kP02XbT35m86CEFWY01lS3Vd1YElgVXA5++7LFDDyH+XyOQqGATCaDjY0NDIdDpNNp7O/vYzabuaVY9GJs8T+LO6isMwZVwUfVVRW2mA/l+XJCBBCaDAE4z+DKlSvY2NhAsVjEf/zHf+DNN998N4bqe4JYkbPT6eDmzZsYDAYuuc3BpfWkquxy7aC26Gcin6onZ1V1PaNgRRAVMOzSJM768/lZpz4d6BxEzM3RWtFCFgoFbG5uuqddcyDbkj9bM2yri6JWe+j3NAyISmloZY2mkZTkdIsZO1YqFZe/Zew9m82ckAfA1UjzWpKQPE7Nc7IhON1eThL2/HXS40Scy+Uwn89dLpQFEolEAh/84AddWmYZEStyAsA3vvENHB0duQoXWrr19XU3MGgBKTywcoSEocJJcYE3U1fT2wLuKERZIR1EHGDpdDoUE7JrnHaHS6VSWF1ddTEf61/tGkT74uSgLyD8eAngrFPAnQitx6cTALepopmS1ZYJFgoFVKtVjMdjdx78OxgM3PnwUQm2pI6qPHAWSmh8qRVQDD/s8dKKssEXLSubTQdBgF/7tV/D97//fU/Ou4XBYOBKySjBl0ol7OzsuPYTdgBSTKD4QouyWCxCj2fQwaluos13qnuraRSb8Ca4fbq+FElIStafsscRt6v5QQ5OTe0A4Qf/qjUEwq07CJu31W2RhNyPnQRYS0vhypKc2ygWi84K5nI5DIdDR3i2z5xOp6E+UCSeHqN6AIQultd8pl1dopVitkRQEbXCaVkQm1QKAPze7/0eHnvssRAJqQCyB6zGg1qiByB0MzUOA6IbeGnqRH9v3UT9q+D2ODgVJIXmPrmqwu7PVrwwliQR1V21IotCP7PHqdtW1VctsnXL+VdFrGQy6VRdFhgwtdFoNFzaolwu4wMf+ADW1tZc/a4toNeCAl0EzmNQa8+HSumx6P1QL0ivwWg0wt/93d/FOvaMfSoFgJP32T1AlzWpG6lNqSi0MB5i8l/zXXbQkgxWeNGcoXUh9bcKnf05wHWg6ZKmqLykFWF4HGq1GLNp6aB+j1ZGJxJVozmBMSWjXoCtDLIFDZoX1e4GGsezCgk4KxBJp9Oh+FnPSUlJNTtq1Qq9AnokTL2okMTQgZabxQ90z5dVqQViRs5er+dKzTTeAMKuI1dmUKggOYGzxmA6cAmSVgehdRVJXIUd7DZO1eoUWlFaaK44YX5OCaKutca0+n/7PeCsaZbum//W7WmljD2ni9w9+1su1mb/ntlsFspNshBBVwnR26F4ZycaXWzAyUy9B80PU3wjOVkVxOPnteY5asypE+UyIlbkpFqpVojgAGQxQiaTcTcGCCf8rQvIuE+LxnWRcxRRSAAr89PaKDk5qNiqhHEniUglsVqtnouZVaji+etTuzi4rZJMK6LXSScOS0a1rHp91bLweDVnTOi+1S236ShuR9Xn8XgcIh+PUUmpKR8tkbSLwvW+WqjbG+X1LBtiRU5dNQKc9f5hvlDVPDugNY+nFnM+nzvll7/R2E9nYWtpo2AHBn9DEYgSP8+lUqm4tM76+robOLod+/Qyuusq+Khl5zUpFArOtdQ+QtyWxpFRxFSLS7dbF4ZrIQNwNpFonTLvjeZqGW9TK6C3QkVX90lLS+Va13bamNpeA1WR7YR+UQy6TIgVORW8yOreqhUhqXTmJTHsQOS/SQIA59IFJD3dZd23dWv5XQ4oKpQ6EGhJGZ9RudSO5DxWkoo1uyxcoCXjseg5sZhB62S5YgdAyAMh7PWwApd6GpzA1AW3nRxoqViex75HVKkZf2otNPvyaorkIiupxx6VUrLagireqrQvK2JDzieffBL5fN7dTLputGQUMOyg0piUMYxaLp39+dK8mpLJikDq+ulgokDDwaTdCfid0WjkVtYwWc7CbF16RqGFzx+1k5B1IflvCmNaFKHfpXBicdGABxC6XqzmUddU26nQuiWTSVeAUS6XXdwHnBUiRKWf1J231jRKGb8IKi4x3NEXSbqMiM1Rs6ojKgnNGZKDWgewVutQaFCVkKKGbXis8Y91Ze0AUcLqxMH3qRDzmOjisaft9vY21tfXQxUwep6qMOrSOY1v7VIoPmqh1Wqh3W6HniGq5LbnFPXiteZfVZFVZNLYkJVP/IzF/Jubm26/mi9W11rjbQt77e137TnxnnN/LLrn/aYVX0bE5qhZSKBdvTlI1B3UWFMHo7pZHMCc9RmvMjazCXggnORWa0orpUKO3T9/T0vCFRksH+RKCTbciqq6ARCqJtLXbDZDr9cLPZKv1+vh5OTEdVGwbUesF2BzpYQqwHr+SlpOKExXcX9R3oquNmEJJknE89PSRz0emw5SEcyq7jbcILT6KQjOuiQuI2JDTsYKdnZVqKXT/KLmyLSGk79RstkKId0/EHZn1fqQlLqciq6kriUsFovY2NjA7u4url69imq16qw13TwOXh4zia2VRcRicVrlpI3C6C5rrGaJqf9W9ZLbtEKJTSdFKbFaJmebcanKygmSxFD9QKt4okjH9+1xcaK0BOa50bvodDpuuWCxWHSx8zIiNuT8l3/5F2xubuK+++5zN5uEoKtKd1DzlBw0FDG4DlRvKr8XFVNGEZcqonXpSE4WQGisxWVhGxsb2N7exs7ODnZ2dhwhbSysMTG3rYKIVvVMJhNUq9XQgm+KSKqwai2vVTsvUi2tEqzxMEllUyBKNr0nem94bDZdw/3weLRM0LquOvGqi0zXnxMbX4vFAq1Wy/X2pXvLet9lQ2zICQDHx8fuUewcDHojgfPla5w9Vc7nZ3Zw6vs2d6aksCV2/B2tZK/XQ6PRwGAwCPUwWl9fx+rqquv/QzfdWkub3giCwHU80Lhaj92SV60VXUZd0sbJg32XdDWPzZ8SVvXUula9LpxQVNCxYQLBSdROjFH3yN4nAG5SY+OwqOunhRuz2cyFMfrk62VErMj5yiuvYH19HdevXw8NDgDnXBM7e3MGp0sYlYjn9zWmtaKQpmWUJLQGrVYLx8fHODg4wHA4xObmJtbW1vDAAw9gd3cX5XLZrdjnYwO1lE/dbSucKKy1IlRh1cqo2WwWKrTI5XKuEwS9Dt1e1LWxAg6PSXPJLPKwlplW0Hoo1oWOihN1G/o7dYM1VrcxcZSSzevjK4TuEnq9npvx7M1VFTHKiqo14QyqM7PK7fy+WiI70JjXo8UYj8doNptoNpuuB+xicdpmZGdnB1evXsXW1pYjJF3xKOug++E52glDixCsi6sWk9vRJL5dd8pyOu6LiIq7aY1UNeZ3VRizBLeWMeoeXHTeeixRrvediKwxMY9bRcI33njD9QxeNsRqVQoAVCoVfOpTn8Lu7q6L+1Qx/L/jOvc7667ZggWdYVV5tXGpgouJtacsCZdInD5F+fr163jqqaewu7vrWoxwWzrzU+a3Si0tPs9B83+qRtuBbolrJxt+zu4R2pzaxrgaL9pkPq2kFX1swYAVcfSv3p+LrLb9DWNbfSBv1GIGu3pG4+7XX38d//RP/4Qf/ehH/7/D8D3FYhlWpQBwAT0HCm+Cqq5RcYom8IHwciyd1fkZEHZhLebz0147XAY1Ho9RLpexu7vrlNjt7W3XbIpxDgeKihYEBzQT8/YceB5aW2tTAnSPgyBwg5CThpKZpGcOVdNQlpz8SxLyGvLBuyqIXaSwMu7jtVMPx94jS1qN8fU3JJ5dxqdekI4PvY4A8NWvfjX2xLwTYkdO4KwTH3vzMJZkvWdUolpf9oba1IIOqItcM+D0gUWLxWmCPZ1O48qVK3jggQdw7do17OzsoFQqIZFIuAIAbsuuLOE2VazS2l5+TwvfgfOF3Cp+TKdT515rB3yb4lCCK1SZJZGZAomKx9WiWffVToB8P6rzPPfL+3YnD0cVbp1slah0562g97//+79L2XFPEUty3rhxA8ViEQ888ICbFXXQAOHHAkbFQPw34yQldNRql6jc3kMPPeQKCsrlMqrVqlv+pUIMEO7Fo7GtQgcqjyEIAleKx0KDdrsdShtoHpcEmEwm54oQrItqiWPFFD1/AK44Q91fe/0ugt22Lj63iq26wlFurhKQExG9AyVuVI54NBphf38fL774Im7fvn3h8S4DYklOrrRnrGlndc7itt0IcCaO8MbqwOINjSrnssJCIpHAxsaGa8jFVAeA0FO61MVSlVOJpNbG7oPuKVMejUYDx8fHIUJpn1f+ZjabodvthkoVdeByX/xrVdKouFCvA7/LAnYSVWt2f1z8SPJctD/+PsoC6+TLEEMnYf23rurp9Xr46le/uvTEBGJKTgB48cUXUSqV3MON9FF4KpKopM+bFSXKWGVX3eKowgCKDUB4hYcdQFwragnJ4yNZWIKoVlstnMbMjLl5Lja+snHZRbWj1v2Pgo0prcLN7TBWtukY/b++uG29L1GwSrk9bo07aUXVutNa0vOo1+uuadqyI7bkpFUAziyPJp8t2TQZr8IQB49aT0tSqsJa30sRZbE4eyamWkf2xNHO5pZ4QHigqVijFpd9dYvFoovVtD8uK5H0mZ3clw58HdD0OPr9Pt58800cHR2hXq+7ByoBp/XMrP+tVqvY3d0957oqQdX7UO9EY2glm6aDlIA6YWje2YpFGnNrkYYq4OxU3+/30el08Morr/zUYy8uiC05AeD5559HIpHA1atXAcAVv0dZE3UtrZhxUdWKWqxUKuV61VAZtYOd9bHss6rVJ2qNaXXVlVQSRVmXdDrtVubwOS0cuHTnqaays4CmEWzCXdXYra0tbG5uotvtOvW52+26xtiaEokqt1NYt9K+NI1jxSSqzozNed2oTquQp3G7Tra6H7rN3H4ymcS//du//cTjLW6INTk7nQ4ODg5cjKVLqrQUTm+4LreKcqXUaurgodVVQunzVehisketded0XxdVpKiaaeNRkow1ujpQmVjXfaqXoNaKx64KKBufcXWMuoW6fz3uKBfaWmlCBR6bZ+UEwWtnvRJ6R1ZYu8gV/3ETx2VCrMkJAI1GA+VyGdvb2859VOHjogGvg0PFBv5WVV97U20cZOMwVVLpCkfl3hRKHrtdWu5EIuHcZe2+HgRBqPKJbny73Ua/30e32w0JaPr8lkQi4QQsHfgki4pEmkuOIqdOXDblclHxh048Sk4WDGhcfidVOCrHqZ4FV+5cJsSenN/73vdw//334/r166hUKiHLCISbEJOQdsmVupwUNigw6IAkSDAuIOYibSq3fEoX1ytykKsbZknIYwWi11rq4KRrTuIytmL1CwWjQqHgGlbbuDPKPbXqpnUlbStJGycSeq1tYYK6mNpGxk5oUUSPOlYb83L/PN/BYOAU7qOjo4sH0hIi9uQcj8d44YUX8LGPfSyU6FerowNlMBiELJ8OWpJGXVB1n2yaQx8ZzydwkRDaAiNKEbUFBrR2AM7FU2qZdTDzMXhs38ImztwWm2/bQguCA5oCEolJgrMMj8c5Go1Ck4Yeq/VMmHvUbejEcJGLaRX1ZDLpOknoJGtFPH7ftstst9uhRxReJsSenABccp4rItSdslZTk986qLSsDzi/mt4WIKi4QjLow4q0iVWU9Y0SgKxF4zmomqn7D4LAPR2axFThSwvzgXADaF2qxlCAqQeqwFyxouV7Ueqr7cyg3d75iES61FymFQWSWq0mgHNN1vScdHKYTqfodrtotVpoNpuo1+toNptot9s4OTnxbu37hbfffhs/+7M/e05sicqTEVFu2UX5NhJnNpu5OlT2tVXS2lUPUdu0cZiqyXpsNmbje2qBVEgh6YrFYqgihs+D4UOfuGRNO+Cp1dZ0BHBWGcT4tN/vu24CvV4P8/k81IaE21RLzueeFgoFly7iZMr9qVut15WTh94rkpeTGX/X7XZRq9VweHiIw8NDtFotTCYTHB4euq7zlwVLQ84vf/nLKJVKePTRR0PvW5HkoqS4VRT5W0LJy+1QZFEBitZBVUhaKP09CabCCAA3cBX6GeNBbTeZSCRcAQObk9Eb0DwiLafG1MxpspCC++GEoLlKHr+mZRKJhOulSyvL+Jd9aQE417/b7TqrrTlje6/UZeXEoJOgTlTWs9HJK5VKhXLilwlLQ85er4cvfelL+MM//EN86EMfupB8mmuzNaK0eLSIwPlid90WxSUlJ902m7+zxxGlIhNR+U6NddXC8Lf8jPGjWmK1NDbXSEKn02mXKuJn2qib4ASUz+dRKpUwGAzQ7XZdXpREVVGGpKOrr5OYikIa9/OceD009tXPVF1WQYmkn8/naDabbrHAZULs1nO+E3ziE5/ARz7yEaecsvpmMBi45L3mQ+mmajdza12ta0pBiMUHLBIol8uhl7qcdMEu2qbFRblEG0NrvMzzYRd7LXinRbECEX/DAgZtJaq5SX7e6/VcM7F2ux0iJl1LXaa2srISUrC1561Nq9hJyy4W0IlQzyMIAtcO9PDwEC+88ALeeuutuzKe3m8slmU95zvBjRs38PWvfx2//uu/jieeeMJZEibwrbvKwUHSqEV6J4lrTdFoa0e2wrQWk9u0fzW3eqd8HrdphRkFq4WiUhAav+r21PJms9lIpXg0GrnKp6ia5Ol0GnKbgbN+sZwE7XpP2/2P58P3tEu/WkqrEAOnZZ37+/uXLr6MwlKSk2LEP/zDP6BQKOAjH/lISAVVocWWkemA10Q6CawzuLqIQXDWDtJaG83zac2oVWX12PjvKGuicZW6fpqSUZddt29dQy0j1JJEftcSot/vOzdUJyUKUNyuvlQR1yIEffyhncT0euujBa3157mxedkbb7yBV155Bc1m810cYfHAUpKT6Ha7+OIXv4jBYIDHHnvMKZma+7QEpAUBwmsqLbGZdqBrSWWTsU42mw2Vp+k2LEFJdIof3A+Pk9+9KM3D93gOVJU1NtVj4fa0n5EtSrCVNgBCi7V5TtrRgRODxqtaqaPXDjh7IpzGnTwXvq+PxuBLl8dRjf72t7+NRqOB73znO/cEMYElJydw2tbkH//xHzEajfDII4+4wcEcoBZZq7AAhLvNMY+mPXIYl1F9tStKWB9Kl5qDlLCWkWTlILbks5ZNFUsrlADhFI5uV91GXUBNgvLYtSHzdDp1KZThcOhynroKhWKaCmr2mSxq+TQeB+COQa2mXauqSjHj2+9973v40pe+hHa7ffcGzhJg6ckJnBL0+eefx2Qywe7uLoAzy8GuBUoctRqqKGqqgVZOXUYOGCbyKXzQQpDs/I0VNGw1UlTMCZyvw7WxqxJAiW6tmHV7dXKw6SW14JoLVVVYixiUTFGlgrz+vL68fnoPKHDRwjN9dfv2bVdUUK/X8Z//+Z/3HDGBS0JO4LRT3pe//GX8zu/8jlNWObiUIFo2Rgui76l7ZX9rn/vB+Iwd34HzZYVKSlqIKFIqAa1FJdTt1X3p7/lbEpLCFYCQ5aerru6oCkyWpJwImN+kVbTKqo1hteiA21W3Vkvz2u029vf38dJLL+G//uu/0Ov1ftLhcClwacgJnApFR0dHePTRR91Tpuh6qrJIdTDKhQTOBriKGZzlSU7tTDcYDEIkj3ocn23tqeSzai7j4iiC2nhULSb3rxMEFyMr2ZQchKZwVOW1zzG1eVEek7W2DBP0mtrSRVYutVotdLtd7O3tYW9vDzdv3rzniQlcMnICwL//+78jn8+7QgWKGEC4DO/HpTIYJ6kVULeRqYwgCNzgB04HueYAAbjnb5bL5ZAbaC2zFYGi0iR6zNbNVYFJxRuNHe35233qShMSmATVTgyqxNrraAmpghDDgvF4jG63i2aziZs3b6JWq2Fvb++eSJG8U1w6co5GIzz//PMA4AiqrpMdjAr9jBaX37O5Q/6bqyPUCjFhTwuaTqfdio+oOEzJonEhEL1o3MaxfGkBuz5mr91uhwrJ1Z1U78Eeh05UJClXyLBELyre5LFRJOOjF9kHuNFooF6vO4v5gx/8AI1G49L0/rlbWMoKoXcKDpBPf/rT+OhHPxppmWyVkBWBCFvlopaUVTGq2FJkoQVKJpPI5/MuHmZxOpegaTxqreNF90gnAaYc9ElkfNIXF1vrkjtum9eJx2crq/hIe64AaTQarn5W3XUSlb8BThfKs3i+1WqhXq9jb28PR0dHaDQaWFlZwdbWFkajEW7dunVX7/0y4VJVCL1TLBanS63+/u//Hl//+tfx7LPP4oMf/CCA8JOztKrF5ip1W1FWheKKFrwT6urOZjP0+/1QA2q16jaGjDoXRSKRcEvYNAWh7m3Uk8l0ItLJg7/Xz7TsURVrfpcuro1deS1u377tVo+0222kUins7e25dZeDweCeVGHfKS41OYnhcIg333wTf/M3f4PPfOYz+PjHP+4GL9c2agc+LdoGzq9eoaum6ySZBiA0ruTA1pUhHOQsaohScS+ymHyfarB9lL0+n1OtpRLSpk/UpeY+NEbU0kWdyGg9aaH5+Wg0QqPRwP7+Pm7duuWeMfrj6o09znBPkBM4HZiDwQBf/OIXMZ/P8fjjj7uk+2AwCD2OHggvALZlZywQ14Xd1uJxAHMhMj/nb+mGMnbjtulOatE7iWSL9xn/cr9a1M+Jwz4/xU4AjK8136iFCfP53C1w5ms+n7uVK3x25mg0QrfbdcXyvV4Pg8EA+/v72NvbA4BL10bk3caljjkvAq3ifffdh4ceesg1rObTrLjihEX0Wv7HAavWqVgsul5DWuygFT6qznLZGZdl8REPfLZmJpMJPcJPCx/YPLnb7WKxWJxbe6mVNyzza7fbzp0GwnlSPTZaQC6i7na7aLfb7tGHzWYT/X4fqVQKa2trrt9tNptFr9fDyckJDg4OcHx8jG63iyAIcHh4uLSP4HuvcE/GnBeBg/773/9+aNnRI488gkceeQSpVMr1kNUKHHUnqd6SbFztoequXTWiKzG0BNCqr5wsooQrjRv1Yba6TX3R3bZikG6X6HQ66Pf7ODk5wcnJCY6Pj3F4eIhOp+NK+hiLcuUKLTQJfXR05NMhdwn3JDkJq8q+9tpreP311/GBD3wATz/9NNbW1kK5uvn87PEIrCmlWkvX1K4Q4Xu0uiS5LRxQYtnPKcxoHx9b0wrAEVYfIagxqY2hKRzx2aMsmePDlFqtloshOYmwMoo9neg2M/bW5WIePx3uaXJacEC/8cYbeOutt/Anf/In2NzcDFXIqBsInIkpUUKHtUxKVhbYk6wqEGkZHBVUjV2jts3tq/us4hMQboOizxjpdDrodrs4ODhwVTu2WwLzuWyo5fHuw5MzAiTDF77wBTzxxBN4+umnsb29HSq/U3LaLgi2C6Ct+NFUBwlFYUXTFlHE1nysElH3TwvGFiO3bt1yj2DodrtOBOv3+6E2mWzWFQRByApyUmD86/He4J4UhH4SqCV69tln8fGPf9zFmfrQW6ZLWE3DAnx+R1eQaM1uLpdz4tDa2ppLR7CggO6lxp8//OEPcXBwgP39fdy+fRv1et2RjzFiEAQoFovugUulUgmVSgXVahXlctmJWXyQEkvqXnzxRfz3f//3+3zV7w1cJAh5cv4EoLVTS/jwww/jypUrWFlZwerqKra3t7G1tYWdnR1ks1nnympVjsKuZdRcJC0xiUoV9dVXX8Xx8TFqtRoajYabGOgOp9NprK+v4xd/8Rdx/fp1XLlyxa2gobWfTqfodDqo1Wq4efMmXnjhBbz66qvv+TW9l+HJ+S5DS/nW19exs7ODra0t9/CgcrmMK1euoFgshhZAa/cCVW3VGo9GI3Q6HZycnOD27dvOUvIhu0zflMtl7OzsYHt7G7u7u7h27Rq2t7fx8MMPu2ZknU4Hb7/9Nm7duoWjoyMcHR2hVquh2Wxib28PjUYjtJja492HJ2eM8PTTT+Pnfu7nnOqrVUl8r9vtol6vu1etVnOd8Jjgr1arWF1dRaFQwOHhIX7mZ34GP//zP4/t7W3cf//92NnZwWw2w61bt1Cr1XB0dITbt2/jpZdewhtvvPF+XwaP/4MnZ4zAHkSsBtIXXVSCz0thC8pcLufaUD7++OMoFov42te+ds4VfeaZZ/DRj34UJycnePXVV3F4eIj9/X2/TjKG8OT08IgpLiJn9BNKPTw83nd4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zh5YircQAAAQhJREFUyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVN4cnp4xBSenB4eMYUnp4dHTOHJ6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQUnpweHjGFJ6eHR0zhyenhEVMEi8Xi/T4GDw+PCHjL6eERU3hyenjEFJ6cHh4xhSenh0dM4cnp4RFTeHJ6eMQU/w/kqZLf0aurlAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Select random layer number\n",
    "maxval = 154\n",
    "i = np.random.randint(0, maxval)\n",
    "# Define a channel to look at\n",
    "channel = 0\n",
    "print(f\"Plotting Layer {i} Channel {channel} of Image\")\n",
    "plt.imshow(image_data[:, :, i, channel], cmap='gray')\n",
    "plt.axis('off');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### 1.3.1 Interactive Exploration\n",
    "Another way to visualize this dataset is by using IPython Widgets to allow for an interactive exploration of the data. \n",
    "\n",
    "Run the next cell to explore across different layers of the data. Move the slider to explore different layers. Change the `channel` value to explore different channels. See if you can tell which layer corresponds to the top of the brain and which corresponds to the bottom!\n",
    "\n",
    "If you're feeling ambitious, try modifying the code below to slice along a different axis through the image object and look at other channels to see what you can discover!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "ce4d9f4bc0e243d198261da79ddf2683",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "interactive(children=(IntSlider(value=77, description='layer', max=154), Output()), _dom_classes=('widget-inte…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define a function to visualize the data\n",
    "def explore_3dimage(layer):\n",
    "    plt.figure(figsize=(10, 5))\n",
    "    channel = 3\n",
    "    plt.imshow(image_data[:, :, layer, channel], cmap='gray');\n",
    "    plt.title('Explore Layers of Brain MRI', fontsize=20)\n",
    "    plt.axis('off')\n",
    "    return layer\n",
    "\n",
    "# Run the ipywidgets interact() function to explore the data\n",
    "interact(explore_3dimage, layer=(0, image_data.shape[2] - 1));"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. Explore the Data Labels\n",
    "In this section, you'll read in a new dataset containing the labels for the MRI scan you loaded above.\n",
    "\n",
    "Run the cell below to load the labels dataset for the image object you inspected above."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "nibabel.nifti1.Nifti1Image"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Define the data path and load the data\n",
    "label_path = \"data/BraTS-Data/labelsTr/BRATS_001.nii.gz\"\n",
    "label_obj = nib.load(label_path)\n",
    "type(label_obj)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2.1 Extract the Data Labels as a Numpy Array\n",
    "Run the next cell to extract the data labels using the `get_fdata()` method of the image object"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "numpy.ndarray"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Extract data labels\n",
    "label_array = label_obj.get_fdata()\n",
    "type(label_array)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dimensions of labels data array height: 240, width: 240, depth: 155\n",
      "With the unique values: [0. 1. 2. 3.]\n",
      "Corresponding to the following label categories: \n",
      "0: for normal \n",
      "1: for edema\n",
      "2: for non-enhancing tumor \n",
      "3: for enhancing tumor\n"
     ]
    }
   ],
   "source": [
    "# Extract and print out the shape of the labels data\n",
    "height, width, depth = label_array.shape\n",
    "print(f\"Dimensions of labels data array height: {height}, width: {width}, depth: {depth}\")\n",
    "print(f'With the unique values: {np.unique(label_array)}')\n",
    "print(\"\"\"Corresponding to the following label categories: \n",
    "0: for normal \n",
    "1: for edema\n",
    "2: for non-enhancing tumor \n",
    "3: for enhancing tumor\"\"\")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2.2 Visualize the Labels for a Specific Layer\n",
    "Run the next cell to visualize a single layer of the labeled data. The code below is set up to show a single layer and you can set `i` to any value from 0 to 154 to look at a different layer. \n",
    "\n",
    "Note that if you choose a layer near 0 or 154 there might not be much to look at in the images."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAADeYAAANJCAYAAADTXA0aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd5wsVZnw8d9N3EuOkiQq+mCAFbPwKrhrzqJrDrhrArO4riIqBtasq/KKviqLWTHrGtEVRVkDYECExwgiQZAsXOCm949Tw+2u6e7p6unu6Zn5fT+f+4Gqrqpzprq66jx10pINGzYgSZIkSZIkSZIkSZIkSZIkSZIkSZL6s3SuMyBJkiRJkiRJkiRJkiRJkiRJkiRJ0nxixzxJkiRJkiRJkiRJkiRJkiRJkiRJkhqwY54kSZIkSZIkSZIkSZIkSZIkSZIkSQ3YMU+SJEmSJEmSJEmSJEmSJEmSJEmSpAbsmCdJkiRJkiRJkiRJkiRJkiRJkiRJUgN2zJMkSZIkSZIkSZIkSZIkSZIkSZIkqYHlc50BSVrsImJT4F5AANsC64BLgbOAMzNz/Rxmj4hYBTwJeDBwZ2AHYEtgSctmB2TmL+Yge9KCFhEnAk9vWfWMzDxxbnIjSZKk+cRYU+MWEYcA32tZ9f3MPGRuciNJkjR/WHaXNAwRcQzw2pZVH8nMw+YmN+pXRJwCHNyy6r6Zecrc5EaSJGlyGCtrPoiIDbVVe2fmeXORF/XH+kxJo2LHPC0YHTovnJ+Ze81NbjTpOhSuBrVfZv56wDzsAxwNPA7YtMtmf42IDwJvy8xrBszjwCLiXsBJwG7jTns+6VBhMuVlmfmOAY53CNOvzxWZubZ57iRJkjQbxppqwlizP8aa/ekRaw4sM5fMvJUkSfOTZXc1Ydm9P5bd+2M9kSRJkiaVsbKaMFbuj7Fyf6znkiQtJkvnOgOStBhFxOHArykvProFkAA7UQLN30TEgePI25SI2Bc4GQPI2XhlRGw115mQJEmStDgYa0qSJEnzg2X3RcN6Ikkao4h4cUQc0/Jvr7nOkyRJ6p+xsiTNLxFxp1oMdthc50maK86YJ0ljFhFHAcd2+fg6YBNgRW39LYHvRsRDMnMYo9L0493A5i3La4EvAz8DrgFap+G+YEx5mm+2B44EXjvXGZEkSZK0sBlrSpIkSfODZfdFxXoiSRqvFwN7tiyfApw3JzmRJEmNGCtL0rx0J9rfe30fOHFusiLNLTvmSVJxBfCqAfa7sMnGEfEIpgeQVwBvBT6amRdHxDJgX+B5wHPYOLvpKuDzEXFAZp4/QF6b5HMv4P4tqzYAD8vMb40y3QXqpRFxXGZeNtcZkSRJkjR2xprt+dwLY83Z+B3wzrnOhCRJ0gJl2b09n3th2X1YrCeSJEmSNF8ZK7fncy+MlWfDei5J0oJlxzxJKq7NzPePMoGIWAX839rqC4AHZOa5Uysycx1wNnBERHwH+BRlxBeAbSkB5+NHmVfg/wBLWpZ/bAA5sC2Ao4CXzHVGJEmSJI2dsWY7Y83ZuWjU15MkSdIiZtm9nWX34bGeSFJXmXnIXOdBkiSpB2PldsbKs2M9l+ZcZp5C++9YkoZi6cybSJKG5Ahgt5blDcBTWgPIusz8AtNHg/nniDhgBPlrdbva8hkjTm+huaS2fHhE7D4nOZEkSZK00BlrSpIkSfODZffFw3oiSZIkSeqPsbIkSZr37JgnSWMQEUuBl9dWfyozf9DH7m8C/tyyvAR4xbDy1sU2teUrR5zeQvMx4PyW5ZXAMXOTFUmSJEkLlbGmJEmSND9Ydl90rCeSJEmSpBkYK0uSpIXCjnmSNB4HATvV1vU1LXdmrgFOqK1+SDWN+6hsVlteP8K0FqKbgNfW1j09Ivadi8xIkiRJWrCMNSVJkqT5wbL74mI9kSRJkiTNzFhZkiQtCMvnOgPSfBYRK4Go/u0EbAWsAa4ALgJ+kplXzF0OZy8idgbuBuwI3AK4EbgM+B3ws8wcSXAREQEcANwS2BS4BvifzPz1KNIbg0fVli/KzFMb7P9p2kfS3AK4H/Dfs8xXN0uGfcCIWA7cA9iLcj2tpFxLfwVOG+VvZY6up49RRvS5fbW8DHgD8M8jTHOahXDeI2Iz4EDKvXYbYDXwR+AH/eQ/IrYE7lXtvxVwFXAe8L3MvL5JXrocfxdgX+DWwNbAKuBqyrPgbOCsUd0rJUnSwmSsaazZgLHmAoh5Jl1E3AW4HbAr5Tv8K/DzzPzlCNNcAvwDcFvKPWJrysizlwI/zcwLRpSu8ackqRHL7pbdG7DsvvjK7tYTRewD3J1y3pcBfwN+C/xv1Yh2YlSdJu9Auc9tR7lGLqPEPr8dUZpLgLtQYp8dKc+Py4AzM/OsEaW5FaXx822ALSnxxEWUmOdvQzj+nD4Xq+v9nsAdKd/jakoc+ZPM/P2Q09qG8tvaGdiB8tu6FvgL8Bvgt5m5YZhpNsjbKsr3vC8ltr2W8pv/UWb+ZUhpLK/SuBXlu14HXAicPqrfzGIXETtS3lnsDWxO+f2eTfleb+xj/10pv4+9KO8TLgPOodyT1w0pj6so70V2p9xPl1XpXFTl8+/DSKdDukuA/YD9gV0ov8crga9l5nmjSFNSb3NdJhgHY+WhMVZefLHynIqIHdhYjt2Mcl+6gBITXTOiNO9Iif12BjZQvt+zKWXnWd8r5rruKCJ2o5QB96Bcv5cDfwJOzczVQ0xnqswXlBhsW8rASFdQ3jWclZlXDyu9AfI38ncg1bm+J6XMO1UPmZR7xazrIdWumtX1rpQ4YwfK7/dCyvn+Yx/7L6NcE/sB2wM3UMpBp2TmX4eYz1sDd6I8Q7al/P4vA87OzLOHlU6HdLeh/Pan7j1rgN9n5hdGleaks2Oe1FBE7A08HngA5YbSa4SNDRHxS+A9wCcy86YZjr0U+AOlkD3lWZn5oQHz+hHgaS2rPp2ZT+xjv02Bw6t996d7QHFFRHweeEOTBlERcR6wZ8uq+2bmKdVD6JnASykNr+peB8zXQv8/1ZZPa7JzZmZE/I3ycJ8ytCAyIo5h+sidrV4bEd0+37vXy8yIuANwNPAQSmGwk/UR8TPg3ZTrtO9Kikm9njJzfUQcDbQWMh4TEXfJzDNGle6U+XLeI+IQ4Hstn5+fmXtVn+1Ubf80SvBfd1NEfBR4ZafKwojYk/Ly5Ql0vlffEBHvB17bJLCOiE2AB1JeDt2XUvnRyzUR8TngbZl5br/pSJKkxcVYs42xZv+MNSc85pmvqr/zRcALaL93tG5zPvB24PghNuDaB3gl8DBK5UG37c4B3gt8qN/KLONPSdKwWHZvY9m9f5bdF1nZfbHWE1WfPRJ4DXDnLrtfGxHvA940x43mdgH+HTiU0omk23Z/Aj4E/Ge/Dd0iYi9KQ8SbZeaS6rMVlFjrpZQGe532/zPwRuCEYcRbVX7eADwO2KTDJhsi4jvAKzLzzIbHHtlzsZbOIXSP6TYHXgEcQemQ12n/c4DXZObn+k2zwzFWUu4rh1Gu76U9Nr80Ir4FnEgZrKXj7ysiTgEObll182+py/aHAf/Vsur7mXlI9dkOlHv5Mygdtzrt/xPgqMz8nx557yoitqXEv0+h+7k+C/iPzPx0lzx/JDMPGyT9LumdQvs5bPW90t68q2l5GUZ+B/hej6H9OXxzmhFxZ8r94IF0vuYuj4i3Ae/IzLUdjn0Q5Zn3j3Qu011apf+BQRtsR8SBlPc5/0Tn9yoAayLiVMo7i282PH7997N3Zp4XpSPgi4HnAbt12PU6ym9Q0hgYK7cxVu6fsfIii5VHKSJOBJ7esup1mXlM9dn+lJjooZROU3VrI+KLwKsy83dDyMsSStzwCjqfZ4C/RcTbgXdn5g0Njj22uqMO5dpnZOaJ1Wf3oVwvh3TZ/Yaqbu2YzLy4Sbq1POwP/BvwINp/63XrIuJ04POU8vSlXY53CO2x5c0xVY88nMccvgOJiAdT7kV3p/OzZ3VEfJxSD3lxlzz3vCc1zM8htJ/DVgd3KL/XTcvLbPM74PfaLc7YhHKPfBGlQ22nfb8LvKTTwErVO4qXUeKUW3TYfX1EfL3af6BBhCJiC0os9HRgnx7bXUDpxP7mJh29Z4hR71p99mCm30/Pp/1d6KLS6yWRpJqIOJwySvabKAWamaa9XkLphXwCcGbVuKir6iXX+2qrDx8wr9tTXqy3qh+7036PAn4PvIMyOl+vUT62A54F/DYiXjBIPlvS3RH4AWUq8m4F0XmpCmbqb3t/PMChflJbvt1gORqPiFgeEe8FfkVpmNYtgITyPLoH8Eng9Ii41SzTnojrKTO/CPy0ZdUS4D9GmeZCOe8RcW/KiC3PofvL+00oLwp+Vr+/RsSjgbMoAW63e/XUi/r/rfLerz8DXwH+hZkDWyjfwb8AZ0fEUVXgLUmSdDNjzWmMNftgrDm/Y55JFhF7AKdTfq979dh0T0rnuO9FRMfGeA3S3KT6Xs+hxE8zxWi3o9x7zo6I28+w7UxpG39Kkvpm2X0ay+59sOy+eMvui7CeaFVEfAz4Et0bpEGZqe3fKeXrPXtsNxIRsSQiXkNp3P0ienTKq+wNHAv8oYofZpP27pTf/zvo0imvsgfw/4AvVB1PZpPm4ygNa59C5055UK7N+wM/iYinddmm07FH+lzsMw/7Ua73o+nSUaxyO+CzEXF81bi/aTqPpvytx1FGzJ/pGDsCTwW+C9ynaXpNVY0QzwGeT5dOeZV7AN+NiKMGSOMfKbMxvJDe53o/4FMR8dmqA4MGFBEvpDxHHkz3a2574M3AN+rnOyLeCJxK6fTQrUy3I6WM+MmqzNIkf5tHxGeAH1EGWer1fa+gdA78RkScHKUj6cCq+8fPKfefTp3yJI2RsfI0xsp9MFZevLHyuEXEi4EzgEfQuVMelImO/hn4eUTcf5bpbQt8k3KP63Wed6CU45rWs81p3VFELI2ItwCn0L1THpRnwbOBX0TEAQOks3VEfBr4BSWenan8uIzye3kr8Jmm6TU1jncgVf3pJ4CvU/62bt/dppTnzm8i4n5N0lC7iLgl5Vn0Jrp0yqv8E+UdysNr+9+OEqccQ+dOeVDu7Q+jlIHuNUAeH0Ypk7yBHp3yKrtTOrb+MSKePsO2/aT9aspz92F0v58uWnbMk5rZusdnN1Cm4O02csEdgJ9GmTK0lw8DrdP33jki7tF/Fm/2L7QHuWfNNM13RLyC0lN51w4fb6BMe9tpFMBVwHuijMI1iC2Bk4EDa+tvokypOhZVgXG/iHhERDwjIp4YEQ+MiDvMsjC6N9NfOMw4jW0H9X0mNoiM0uP/vykv/js9a3p9t3emNFS764DJT8T11KJeqfGAqmJk6BbKeY+IuwHfoFQiTFlPmXZ82ih/lIaiX4uIzar9Hw18tsrTlHXAldVx6m4PfLHB73ybLuvXUe6TV3VJZyml4vbdfaYjSZIWD2NNY81BGGtONy9inklWNRL9HqVRRCdXU/7mVvemVAgN1Fi0agjxHcr3urzDJjdR7oOd4sHbAKdFGf19kLSNPyVJTVl2t+w+CMvu0y2msvtiqSdaDnyO0liu1Y2U8nGn0dJvA3wrZtnxrImqLP9Zyqj+nTqQrKXEA51mbdkZODkiDh0w7V2B7zO9wd51wLVddnsEpdHsQCLiyZSR0Vs7aq2nfCed/sblwIkRcd8+kxjHc7GrKLOMnALUGyhfQ+fnFcBzKbMZNEnn9ZRZFzo9H6Fc593iSOjdcH/Wqu/rm7Q3VJ16bnc7/8dGmRmu3zQeAHyNzg0Lp9Kq//2PpTQid9CYAUTEEZS4vrWx4xq631PvR+nQO7X/O4FX0X7+b6B8V508njIzX7/525FyT6t3bpmyGvh7l8/uR3nmDNq54JaU91f7dkiz2/1U0mgZKxsrD8JYebrFFCuPRUS8EngX7fVP6+heft8c+EpE1MsZ/doS+DZl9tBWq+l+nu8JfH4e1R0dD7yc9nLmGso57ZTujsC3I6JbJ6VpIiIoA1Q8nu7xxFWUeLqTUccgI38HEmXG9i8BT+qyyY2U2LfVNsBXB+nsJaB0rP8uUO9IejWdyzGbAp+JMmDQ1OAhp1C+6ylT5YRO72C2pNQ79+oA2CYink25Lnbq8PEGutdRb01539PofUgt7dcCr6f9+TWV5rpBj7uQ2DFPGsxlwEcpI5zdCViVmZtm5g6ZuSnlZeRDKS+5Wx+w21JGYVvR7cDVVKGfqK0+oknmqgLac2urj59hnyMoPbxbCySXUKYbvQuwIjO3zczNKQWlpzF9yuqXRYMR9Fq8iTK9O5RA6TnAbpm5MjO3oRR2Hwj87wDH7teulIfDr4AvU0ar+CTlBfavgcsj4kvVC+em6iO7QBm1oqn6PrtPNQQbgv+mjCQ09a8++szXap+3/ru8w/GOo3xnrS6n9Lzfq+W73RJ4DNOnod+REuAMMkraJFxPN8vM71IKa61GNRrqQjjvmwInVduvpVQ23gPYJDO3p7yQuTfTp8K+LfDyiLgt5f68jPLS681VvlZk5nbV/g8Gflnb/0DgGQ3+3vWU38mrKYH0Lpm5vLpPbgusBO5ICULPq+37goh4bIO0JEnS4mGsaazZhLFmMd9inokVZXTazzG9QeNpwKHAFpm5TWaurLZ5BRsbc92D8rtumuZyyqie9VknfkipzJo6xztQZnbYn1Jp2NrIamvgcxHRqQKiF+NPSdJsWHa37N6EZfdiUZbdF1E90eso9z0oo3Y/G9g9M1dV5ePNKKNpn17bL4BXDpDeoP6L8re3+iXlOtknM1dk5vZV3HMbyqj2l7VsuxL4aJQR0Jv6LBtnM/ge8Ghgm8zcIjO3onSqejZwcW2/p0eZqaypO1IasC+hNCB7F3B3SsyzXfU37keZIaa10dYS4EPRbOaskT0Xe9iU0hhtamaJLwIPAjbPzK2r59VulOu/3gD2qCqmm1FEvJwSj7U+HzdQOuodCuxQXefbZ+YKSoehqQ6Vfx3g72pqZ8q1tZLSEO4E4GDKd7Btdf73ocTSN9b2fVeU2Tx6qjqVfobpDedPosyAtqqKTzehxMxHs7GB6KMo1/UovZONz6Qrap+9i+7Pr8OBj4w4b4PaH/jP6v+vpHTyvk1mblLdU7egdIj7U22/p0TEIRHxJOAl1boLKY3/96h+l9tSGuw+i/b7G5SyXKdyS5sos05+ilJGbPVnShl0l8zcLDO3pAyG9FTgN7Vt9wG+NGDn7A+wcZa8X1LeId2iSnMryqwwjwXOHeDYkmbHWNlYuQlj5WJRxspj8gBKORhKTHAMJQZaUdUDraTESPXZ1VZRyhuDeD9lhm0os2Y9GdixKqdsQymnPAH4XW2/Q4DDGqQzV3VHz2Jj+f4S4EhK7Lyy5ZweQhmQs9UOwNv7SaC6vr/N9NkGL6LUTf4DJa7dNjO3oHxfdwNeRJkRslOnpGEbxzuQN1DqGVtdDLygJa2tKb/dR7Dx3dMqyvNylLOX/5aN97qP1j77Hb1jsG73x0lwPBufTScDDwE2q+rHN6UMDFIvN2wKvLfqSPkFyn16A6Wsc182xuYrKR2sP1/bf1vKLI8zqgblOZ72wVs2UL6De1O976HE5ncB3sv0Dsivq+LFpu7OxkGObmTjO6apevGVlPPzjgGOvWB0GoFYUne/ofTA/0JmdhvtjMz8G2Wk8K9HxHGUBkZTL4QPoLwgqweKrY4Dntmy/LiIeGlm9vswehDtjaiuBT7WbeOIuAvlZWWrk4BnZWa9Rz2ZeRnwsSjTBB9H+4vU4yPi5MysVxj0cofqv58FnpqZbS+EM/N6SkFrlFZU/7rZFngk8MiIOBN4Wmae3eexO00z3eT89NpnW7qP9te3zDydlkJgRNyTMhLHlNMzs6/RGCPiMUwPEn4GPKT6bbSm+3fgCxHxRUpB8lUtH+8CfIjysr6JSbie6o6iTN875V4R8YjM/MqwElhA533H6r9XAQ/PzB/WjrMO+GH1QufLlMLvlOdTCrNbAOcDD8rMc2v7rwG+GRE/BE6lfRaI51FeIM3kncAHM7Ne0dGazlrgbMoU8MdRKjUPa9nk2Ij4fGZ2Gh1FkiQtPsaaGGsaa/a2gGKeSfZSygv0Vm8CXlWPXap46C0R8VFKxcQdaP9u+/Vm2kdnvRF4TmZOaxhX5eEs4KyIOJFS4T1VKbczpSKiycwVxp+SpEFYdseyu2X33iy7d7QY6ommyvUfAZ6dmW0jgWfmDZRRwL8DfBW4f8vHz4mIN1Zl6JGJiOfTPrPTBkqDuLdl5rSGe5n5e+CtVfzxZTb+JjYHPs70DikzOZDSceqFmfm+DuldDnwwIr5JaYC7W8vHzwf+p2F6U/n7E+U6mNZBJTN/DTwvIn5G6bQ45VaURoYzXaPjei52smP173rgKZn5xQ7pXgi8PSK+TRkAZmq28xWURvkv7ZVAlNkt6x1pLwEOzcyODZYz8yJKg9GvRsQLKQ1xL+3zbxrEVKPBvwGPyswfdcjTH4CjI+IUSiP6qUZ821A6TrxnhjSOp31WjrXAkzPzpFo6GyiNMI+tfjffoTTMG+R9Qd9a76VRZjRqfe5+JTNPGWX6IzI1Q8NZlHcOF7V+WD3LPhsRP6DMJLJHy8evZ+NMQd+hXK/X1va/mtIB90eU59PUb2M5pVx35Az5ewmlU2ar/waekJlts5ZUHWk+HhEnUWb0e3rLx/tR3g29eIb06qae9f8JHFm/h1d/b73Bq6TRMlbGWNlYubcJiNkm4Xoat6lZw84EHlb//VVliJ8BT4iIc2kfgPI+EbFfZp7VMM2p2Pg/gKM71K9dS5lh61uUmbX+oeXjF9Ael3Uzl3VHU3/ft4DH1e+FVbrfr8qpJ9TSfEJEHFm/3jv4NO3lWygdHl+amavrG1fX8tRv9T3VQBODDG7TxEjfgUSZGbMer55Kqc9sG3im+u1+lRKDvpTSMWpvRqiKT95f5fUGSqf0KRf1e1+cQHenvCt6cWZOi5MzM4EjIuIC2t8VHEzptL4fZWCkJ3V5R/Fz4LER8UHayzOPj4iX9CrPRMRWlDJL66RsNwCPzsxv1tJZR7nvnRkRn6C8B2iN6Y+PiFMz84Ju6XUwFWNeBDywep9UTzOrf4uWM+ZJDWTmVzLzpF4BZId9fkQJPFu9YIZ9fkl5iE5ZRbORtOujwXysKqx38zZKb+UpXwKe2CmArOVzDeWF9ddaVm9GGXmgqZ9UadZHaZtEdwZ+FhFP6HP7LTqsGyTwm1ao7HLsuVaf6vYCyoO4a4E6Mzdk5tGUoLHVIyPiHzrtM4OJup4y86eUURpbHVuNJDcsC+28P6neKLJVdR8+gvYRTrYH7kNpyPnwThWMLfv/nen3qjtHxK1nylhmHtUrsO2w/WrgX2i/r9+W6SMgSZKkRcpYc1o+jTX7Y6w5v2OeQRwcERtm8a9rJW01Ovm/11Z/tIp/ulbKVZWY92P6iPQzioi9md7w6knZoVNeh3R/T+ko13o/eVQMNmuF8ackqW+W3afl07J7fyy7L76ye5tFVE/0DeAZ9QZptTRvpNzPWr+bnSjl65GJiC0pDVpbvSwz31Lv0FGXmZdSRrpvbbB054h40ABZOTo7dMqrpXcBG2e6mvKwGGxWqWsp10HPWaMy80Ta7+UA/zzTwcf1XJzBv3Zq8FZL81eUzkqtZvz7KKPWt45EfwVwYHbplNch3TWZeWJm1mcKG7a1lMZ40zrl1fLzHarGky16noeIuD1l9oVWL85ap7wOaV1IeV9Qn61Q/bsCeHDWOuW1ysy/0t5QH8pMCTtQZrF4VNY65dX2P4fpM5c8rtO2U7q8QzoDeGzWOuXV0rqJ8s7iW7WPnhsRO3bYZSafy8yXzHQPlzQexsrT8mms3B9j5cUXK4+snquHv1DO80ydPl9PmRmyVc9yUQ8fyMxpg162ysyraO+YA3DAPKk7+gWlnNn1Xlj97c+nfRbxTSgdebuKiIcC/1Rb/a7MPDw7dMrrknZmZs/ZUIdklO9AXk57LHoepRNvz/gqM9/JIp+xbAjemh065dW3ocRbraYGj33hTO8oKAOhtH6XmzBzZ+tnAresrXtK1jrl1WXmTygxfWvctBUzDFTUxVrgoVnrlKeN7JgnjUH1grN1Stq7VS//e3lvbfm5UaZT7yki9qR9BG8oox102/7ulBG+p1xDGZm8r5dXVQHqJbTftJ81QEXS87L0mB63KylTUT8TuAdlVL2VlOBsL8qU4B+m9CxvtSnw0Yi4Xx9pbN5hXf14/Zj4IDLKVLn711a/KDOv7PMQL2P6iIFNR0eDubueejma9t/JHYFBpgSeZgGe969k5jdm2igzz6eMbll3fPYxWk1m/gD4c21109FN+1LdK99SWz3qkVkkSdICZ6zZF2PNmRlr9m8SY81heRyls9mUq5l5dHQAMvMSpjcC68fLaK9U+kxmfqHfnbOM+P/ullVLgOc0zIPxpyRpLCy798Wy+8wsu/dvEsvuC72eaC1weK+Gh1OqDjv1TmB3bZheU8+lfXTwHwPv6nfnLKOWH9PhmE0kpfFWP75I+3eygvaZFPr15sz8XZ/bfqC2PLLvZMDnYicnZ+an+9z2BMp1OmW3iNip28YRcX/gbrXVz23SAHaM/it7DDhTU/+eD5jhmV2/zn9Oj3JFq+q3/ro+86Xp3lidw5l8gc5lhJf36ijX4uO15d1m6Cj3ROAWtXXP7qeBf1W+fC7tDZNXAof3kc9Wa4AXNtxH0gQyVu6LsfLMjJX7N4mx8qj9e6/Oj1Oq33a9A+QgMdGVlO9rRllmZDxzCGn2k9Yw646el2VGuJnSvI7ps6DO9PcdVVv+JdMHhZgEI3sHUsWp9U5aL5mhI3mr11Bmeldzl9E+c2ZH1X2007uIX2TmB/vY/xrKrMGtutYbV+Wc+uAFX8nMvmYJz8xTmX5/+9eIaPqsfH9m/qLhPouKHfOk8WkdNW0p01/i1n0RaH3Jdmv6G6HgubT/tn+QvacHr1f4fDrLqH99q17mtwbJ2zE9kOjlzMw8o0maQ3AR8BRgl8x8QmZ+ODN/mpmXZeZNmXldZp6fmV/IzGdSprE/uXaMFcBJEbHtDGlt2mFd11ESeuj0IrPTsedS/QXG+ZTRgvpSjehQn467fsyZzMX1NG/7/IUAACAASURBVKMsIyF+rLb6dRGxYgiHX2jnvV4I7OVns9z/9NryILMr9Ks+euY9R5iWJElaPIw1uzPW7I+xZn8mMtYcovro95/up6KyxUdon72up6rxwRNrq9/dadsZfKq2fEjD/Y0/JUnjZNm9O8vu/bHs3p+JLLsvgnqir1cDWvSrPrPXvg3Ta+rJteX39NOAruaztHfsuk8/jcBbfKBBQ+51lNksWjU9R+uB/9dg+/p3ctshz+pY1/S52EnfsyBk5hXAObXVvc5pfSa5c4DP9ZvemDWZDeLXtMfvmwO799i+Xvb4QMPfzn8xWAP7xW4d0+/LHWXm9UB9VsZLga/2uf8fKbPzter1zqL+fDg1M+sNynuldx7Tn1FNnzlfzZlnvZE0fxgrd2es3B9j5f5MZKw8YlcAPWd6rhlGnDrTTJyjSLNfw6g7+lVmntZg+77/vojYFTiwtvrNWWYhnTSjfAfyj5T76JQLmd6Jq6sqPvhI/1lTi4/1M+BIZbb1xvX9e8Vgt6d0im/VtF69PjjVlpTZ1pto8o5pUVo+1xmQ5rvqZfStq39bUW5WnSpQblVb3qPXcTNzbUQcD7yxZfXhQNdpRyNiE8qUw61mGq3skNryjKOEd3EGcPeW5XtRpizux/cGTHNgmflbpk8l22v7iyPiQZQKl0NbPtoWeCVl6uBuOr1o3qTL+l5W9nnsuXRQbfmLA1RqnUT7KBc7RsQ+mfn7Pvcf+/XUwDGUhoebVMu3Ap5Fn6MK9rCQzvsG2qdNn0l9xoErmF7x0GT/bTpuNYOI2Bq4A2V0wC0plVgzVcb2fA5IkqTFzVjzZsaa3RlrNjMpMc9s/A545yz27/XbqVe+9V0hDJCZqyPim5SZ9/qxP+U6n3J1ZtYrBPtJ95yIWM3GSvb9ImLzPkeEN/6UJA2FZfebWXbvzrJ7Mwuh7N7NMSzceqJTGqb3x9ryQOXjfkREpwbQXe+l3WTmtRHxOzY2ltqW0piu3tmrm1MaJjnbc/TrJgOuZOYVEXE1sHW1ainlmXZ1k0RH9Vzs4vsNt/8jsF/Lcq9zet/a8icG+F2Nw1WUWST6kpkbIuJPtM/AuA2lcXib6rdzm9rqrzfJXGZeFRGn4UzuTZ2VmVc12P7PtM+w8KN+OwK37L9dy3Kv30a9sXRfszTUnAQ8vmX5gIjYNDM7zXjUySQ/6yVVjJVvZqzcnbFyMwshVh5lPVcnP8zMtTNvdrNhxKmnNNx+KLHxGOuOTmm4fZO/rx6DrWawsuY4nNJw+ybnoV5n+62GZXsoM/RN4kyDk+4HDbat1/vC7Oqde10T9WfIFTS8p2fmuRFxFu3vRA6i//LN5Zl5VpM0FyM75kkDiIiVwBOqfwcz2AgbM40IAvBB4NVsDB4eGhF7ZGanGzrAYynThk+5BPhCt4NHxGa032QB7hwRO/eRt7q9a8u7NNh3XtysM3N9RDyVUvDZteWjZ0fEUT0K8Z1GwFhF8wCw03XWZHSNcahfT/XR4PtxFmUkm9ageX+g3yByYq+nzDwvIj5A+7TCr46IE6uRKga1kM77NQ0rGeqNLf/c8MVF/Te0Vb87RsQdgMMoL5bqLwr70c9zQJIkLSLGmh0ZaxprwsKKeWbjosx8/7APGhE7ALesrR5k1NQz6L9jXr1S6eqIeO4AaUIZQXfqOl4K7MT0Cq5OjD8lSQOz7N6RZXfL7mDZvacFXk90XsP0rq0t910+HsDdaW+QeB3w+IgY5Fj1Rr+70H/HvPMapjXbc9Q0vak0t25Z3oo+OuaN8bnY6ppqFrwm+jqnEbEVsE9tdX2Gg0nx5wEaafZ7be1L+2/nb5l5QcO0AH6OHfOaajL7Bkx/Z9F0/77eWUTE9rSXgWCwZ059nxWUTs/9zrw3sc96abEzVu7IWNlYGSYjZpuE62kk9Vw9nNdw+2HEqWNLc47qjs5ruH2Tv+8uteWfTuhseTDa81CfOa3v2albNO3EqqJJHNVpkNgm+zepN64/Q84ccOCg02vHajKL7yQ8QyaeHfOkhiLi/sDxlNFcZmPLmTbIzEsj4iTgqdWqZcBzgFd12eWI2vKHZiiY7Ej7FO1QgtZh2G7mTW52+ZDSHLnMvD4i3gr8Z8vqrSmB5Q+77NYp0NuMMnpdE52CyHqBbc5UowttUVvdT0O4Npm5JiIuoL2yY/sGh5j06+lYyihMm1fLOwMvBN48yMEW4HlvNOomsK62fM0s91820w4RsTll9J5nMv0e2sSMzwFJkrR4GGv2zVhzI2PNBiYo5plE9fNwfWZeNsBxzmuwbb0Bwx6Ue+AwbEd/14jxpyRpIJbd+2bZfSPL7g0sgrL7Qq0nano991U+jogtgSc3OO5FmfmV2rp6/LE5w40/+jWSczTE9AZKc5zPxZpR/n236LCu3wbP4zbK81BvqHvxAGnNZr/FbFLfWXR6LgzyzPlzRKylvd3iQnrWS4uSsXLfjJU3MlZuYBHEyqPQ6FrKzHW1AVwGqZcZedw3x3VHo/z76nHYpMZgMNrzMOs4LDOvjYi/M/3eo96axGH17xSaxWFNron6vb7xM6TyhxmO28tifIY0NpsbsrToRMQTgK8z+wAS+v/9vbe2/K8RMW1a94jYn/bpStcBH5jh2E1uqk1tPvMmN5u0EUpmUq/MgfYp5+uu7LBukBF0Oo2YM8iL9lHpNIJG0xe+U+oFnCYvJSb6esrMv9L+EgLg5REx0DTkLLzzPshIDsPcv6dqhMxvA89m9uWomaaLlyRJi4SxZiPGmhsZazY3CTHPJKqf42Gd314m4T5h/ClJasyyeyOW3Tey7N7cgi27L+B6olGVj7enNPDu999LuxxjVPq+1w04mvlsjDy9OXouThnl39fpmpmk+22rUZ6H+n1p0HvLoPstZpP6zmJSnjkT/ayXFiNj5UaMlTcyVm7O52cz447BRp7mBNQdjTMOm6TfcZ1x2MI0q+91hO99RlWX7zNkyOyYJ/UpIvYGTqB9xKb1wDcpL/jvC9yG8lBclZlLWv8Brxsk3cz8GfCTllU7AY/psGl9ZJevZuZfZjj8JoPkqU8LtsFPZv6J6VOm79Rrlw7r9hwg6T1qy3/JzE7T4WryvQ24omV5W+Dlc5QXNfNW4MDauouBdwGPA+5MeUm0BbC8w7NAkiSpjbFmYwu2TGWsqUXG+8TMjD8lacJYdm9swT6PLLtrCKwnGi/vdSMwV89FjdWNteVBf0uj/A1KkuaYsXJjC7b8aKwsjYV1RwufcZg0Ty2feRNJldfSPu31X4BHZuaZfe4/mylh3wvco2X5CODTUwvVKAhPru3zvj6O22lq0dtn5jmNc7j4XAHs2rLcaWSTKX+kFJZWtqy71QBp7l1bnrTvqdMoNlsNeKyta8tXdNxqnsrMqyPiLcBbWla/KCLek5mXNDyc531MqpeJz6qtfg/wssxc08f+Tg0uSZI6MdZUK2PN6Yx5Rq9+jod1fnup3yd+mpn36LjlImT8KUkTy7K7Wll2n86ye5+sJxq7+r3u0szs1UBY/ZnL5+KodXo+bsPiGyG+fn8ZdGbPQfebzxbqQPndnjmDdH5YqM8caTEyVlYrY+XpjNk0FIug7qh+712McQQMIQ6LiCU0q7ddCBZqDAajq8v3GTJkC/kilIYmIpYDj6itPqxBAAmw4yyy8Fngry3L946IO7QsP432IPW3wHf6OO5lHdbt0Dx7i1K9sFOf4vVmmbmO6SO8DNLYrL7PRAWRmXkT0ysj6oHvjCJiBbB7bXWnFx7z3XuBi1qWNwOObnoQz/tYPZr2stOpwIv7CWwrs3kOSJKkBchYUx0Ya9YY84xF/TxsFhG3GOA4ezXYtn6f8B7RzvhTkiaMZXd1YNm9xrJ7Y9YT9SEzz6uP8D/Dv0M6HKZ+r9uuaqSmAU3Ac3HUOj0f9xl7LuZevaPwXhGxsuOWvd1uGJkZo3W15WUDHGOhNiLu9FwY5JmzO9MnE5jzZ46k5iagTGCsPHmMlWsWesymsVrodUf1e+9ijMFgehy27wDHuDWwYgh5GafZxmELNQaD6ff6xs+QSr3zu8+QIbNjntSf3WkfvePCzPxuw2PcbdDEq8L5/6utbp1m/fDaZ+/PzA19HPdK4Pza6gOa53BxiYhdKZVjrS6dYbf/qS0f1DDNAOqN4ppeg+Pwq9ryXQc4xh1pHwmn03HnvcxcDbyhtvrZ1cgmTXnex+Mfassf6+de22Lg54AkSVqwjDV1M2PNnox5Rigz/0Z7g2CAuwxwqCb7/Ly2vFdE9Bo5d7Ex/pSkyWPZXTez7N6TZfc+WU80VvX4Yzmw31xkZAGZ0+fiqGXmNcDvaqv/z1zkZY6dS3sj8hUMdn8ZpIH9XLq2trxlk50jYimwx/CyMzky83Kmv0Ma5Jqo77OGCetEIalvxsq6mbFyT8ZsGoaFXnd0em357lWH1MXmZ7Xlew1wjPkWg8Es4zCaDSA739Tv9XcZcMCp+rPHZ8iQ2TFP6k99pIA/N9k5Im7LYL3WW70fWNuy/NSI2CIiDgZu37J+NfBfDY57cm35YQPmbzF5aId1v5xhny/VlneNiHs3SPMJteXrgG832H9cTqstP2qAAsDjasuXZubvZ5GnSfZh4A8tyyuA1w1wHM/7eNSfBfWXcDOpjxImSZJkrKlWxprdGfOM3o9ry49qsnNEbAo8qMEup1GutylLgQc3SXOBM/6UpMlj2V2tLLt3Z9m9GeuJxiAzzwPqf4v3utmZhOfiqH2vtvzkxTbTYjWLzU9qq5/c5BgRcQ/GO9PF2tryILPdXVVbrs8qMJMDgK0HSHe+qD9zDh3gGPVnzs+rDuuS5p9JKBMYK08OY+XujNk0DAu97qgeg20KPGYuMjLH6veL+0bELg2P8ZRhZaYPw4jBYPZx2H0HTHc+qF8T2wEHNzlA1ZF9/9rqH80mU5rOjnlSf+pT/Tad8vQVs81AZl4EfL5l1ZaUh+cRtU0/lZn1B1QvJ9WW7x8Rdx8gi4tCNQLDy2urrwe+P8OuP2T6CDDPbZDmv9RWfyMzb+hn/zH7Wm15bxq8mIiIrYBnzHDMBaOaRvw1tdVPBu7Q8FCe9/EY+FkQEbdi+gsSSZIkY00Bxpp9MOYZva/Ulp8QEds32P9pwFb9blyNZPzl2upXVqPKy/hTkiaRZXcBlt37YNm9AeuJxqp+r3tRRGw+JzlZGOb8uTgG9WtmX+Cxc5GROfbx2vLTI2KvBvsfM7ys9KU+y8IgHeTqM7fdLiK2a7B/vWy20NSfD/eJiPrsLV1FxJ5MHxBqoT1zpMVkzssExsqTwVh5RsZsGoYFXXeUmRdT7gmtXrEIZ807DfhTy/Iy4NX97hwR9wIeOOxM9TCMGAymx2F9z1pfxagLeQDY39B+TQC8oOExXlxbvgY4deAcqSMbOkj9+Utted+I6GtUr4h4ONMLxYM6rrZ8JPDo2rr3NTlgZp7M9N7Un2jY+KrNpI8UN8v8Hcf0Ed0+lZk39tqpGk3u7bXVT+xzhJdXAHvU1r2lj/3GLjNPYfr0tu+ugsN+vBXYqb7/bPM14T5F+zlbCryqyQE872NTfxY8vJ+dquDw48DyoedIkiTNd8aaDRhrTmesacwzRJ8BrmhZ3hp4Zz87RsTOwLEDpPkGYH3L8h2BdwxwnNa8TPR9ogHjT0maPJbdG5j0Z7Jl99Gx7D4Q64nG4520NxbbEfjIbO4Hk36vG7FJeS6OTGZ+F/hpbfX7I2LvucjPHPoMcHnL8mbAxyJis5l2jIgXAQ8aVca6uLi23LSjM5n5V9pnfFrO9Eb/HVWzMz29aZrzzKeAy1qWl1B+GzM2mK4GZDoeWNWy+kbKbFeS5qdJKRMYKw+BsfLoGLNpSBZD3dF/1Jb/gQn9XY9KZq6nlJlbPSci6oNbTBMROwInUsro41KPwfaJiE0GOE49/n5SPwMqVTHGBxh8pr6Jl5kbmF7WObQqS80oIg4Cnl1b/eHM/Psw8qeN7Jgn9SEzLwXOblm1BDghIjbttV9EPIbpo6fMJh8/BH7RsmofoPXl1k8z84wBDn0k7aMp7AP8KCIOaHKQiLh1RLwD+OAAeRind0fEWyNi1353qKa4/xjTH06rgdf2eZj3Ahe2LC+hBOzRI91HA0fXVn8+M0/vM8258Pra8t7ANyJi2247RMSSiHgd8JzaR1/KzJmmtJ/XqkJTvYK16dTT4Hkfh/+pLT8pInpO8R4RW1NGKLrXyHIlSZLmLWPN/hhrzshYswtjnv5VI8i+rbb6aRFxbK/K+KpT3slA48YImXku8H9rq18cESdExBb9HicilkbEgyLiGyyc0RCNPyVpwlh2749l9xlZdu9iMZfdrScaj8y8nOkzdz0G+GpE3KLJsSLioIj4DNPvB4vGpDwXx+DfgHUty9tRno/37GfniFgREU+PiNuPJHdjkJmrmT77zf8Bvtvt74qIzSLircB/VqvGOWvNmbXlpzVo8N7qs7Xl18QMs8JFxH2AL7CAG4QCVB0s6o2j7wl8utc9oGoQ/kGmv7s5vrqnSJqHJqVMYKw8NMbKo2XMptla8HVHmfkNSr1jq5dExPERsarTPnURcduI6GvmzQl2HPDbluWlwEkRcWREdOxgWXW++gFwW2ADcNPIc1mcBaxtWd4UeOYAx/kKZdCOKTsBx81QT70S+CjwgAHSm28+RPuzEuBTEfGPvXaKiLsBX6W9z9g1wLuGmz3B/Oj9LA1qy1k+XP+WmZ9rWX4f7Y2F7g2cGRGvp0yBfRWUYAM4mDKl9tR00xsoU37eZxb5mfJe4MNdPms0ssuUzPxxRBxOuXFPCeD0iPgq8AngtMy8+aZe9TLfDdgfuDtl9IU7VR+3Tg8/ibaiTON6ZEScCnwZOAM4KzOvnNqo6rV/Z+AhwOHADrXjbACe1XpeesnMGyLiBZSXsVN2B06LiDcDH83Mv1bnNoDnU66j1gfiVcDL+v5L50Bmfj4iTgQOa1l9IHBuRLyJEgRfAFCNaHB/youM+tTDFwPPGnmGJ0Bm/ndEnEY5T4Mew/M+el8FLqD8bqH8Nr8QEcdTXp79OjPXV8HAbSijb72MjfeOU4BDxplhSZI0EsaafTLWNNYcJmMeAHYdQkXS9zIzu3z2duCxwF1a1h0FHBIRbwO+nZnXA0TEnsDjgVcC21Tb/pjSEKuJIykz5d23Zd0zgEdExIeAbwCnZ+Z1Ux9WjTpuS7lP3JdyT5xqSFvv6DdfGX9K0nBYdu+TZXfL7sNk2b0564nGIzPfWTWifkrL6ocCf4qIj1DuBT+dup/DzfeBfYD9KPf9R7CxnH7qWDI+uSbluTgymfmDiDiK9k5Iu1Du25+nPB+/X3t27EqJqx8MHEppTHhf4Ddjy/iQZeYJEfFIyvU/5Z7Ar6rn6I+AS4EtgdtRnptTDcyvpXTQe3XLvhtGmN0v0d7ZeR8gI+JLwO+B62vpZ2Z+r8Nx3k8pK0zN9LAl8IOIeCPwyalnf9UQ9F6UdxlPoTyzLwf+RnmWL1TvonzPrY1ADwXOjoj/AL4y1dkuIrah3GtfQXkH1OosyrslSeNlrNwnY2Vj5WEyZgNGX8+10C2WuqMnUe49rTNiPhd4eES8h1J3eE5mroWby+R3BA6ilEnvTXm2zNtZmTNzdUQ8nfKdraxWr6DU5R4ZEV8D/kAZBGUXSrn8ri2HeB/lWbtny7qRxGFVXr9FKfNPOS7KDH//S5lte21tt09k5rW141weEZ+gfbbyw4DdI+JY4IeZuQYgIvao0vs3SidnKHHpQcP5qyZPZl4TEU8FvsPG59/mwHci4gRKGehnmbm2el7uTzl/z2N6f7HDp543Gi475mkh247p07k28UugNYj8IPBU2hsY7Qt8EiAirqbc7LbscKzXUkbFGkYQ+UnK1NT1kcivAD4z6EEz88NVw6Z3snHEmKXAI6t/RMQ6ShCzCbAF453udhSWUgL+g6dWVH/j1ZRz0Om7nLIOeFFmfqJJgpn5xYh4LfC6ltXbUb7Tt0bEdZTzu6LD7jcC/5yZ5zVJc448n1Lge2DLuh0pL2jfFRE3UgqFW3fZ/1LgEZn5t5HmcrK8Evj+LI/heR+hzLyxehH0RTbe/5ZRzvvzgZsi4lrK+a2XsS6gVMbUp5SXJEnzj7FmA8aahbHm0Cz2mOc2zO7+A6WhWMcKy+pF/WMosWlrRc2BlDho6p60io2VQFN+TLn+vtEkM5m5JiIOpdyDWkdM3x749+ofEbEauI7y+6mnveAYf0rS0Fh2b8Cye2HZfWgWe9l9ENYTjcczKQ3CDmtZtzlwRPWP6jxdW63vOdvLIjcpz8WRysy3RsSWtM/IsoQysM1j4eZr5u90jtEWisdTGsu3xu7LKA17D+myzw3APwP1WXdGNoNeZp5eNVJtbRS6M6UhbycfAaZ1zMvM30eZnebYltVb0f68XkP5zlvLQ1N/82tZwB3zqgbfT6S8B7pzy0d7U+4NH4yI6ynlpW7lqd8Dj8rMcc6oKKkwVm7AWLkwVh6axR6zjbSea6FbLHVHmfm3iHggpSPiPi0f3ZIyaMpbgA0RcRXlnrDF+HM5elXn8EdSBh9pnS1wF3rPSPct4CWUjuGtRlnuPpZyX5u67pZQOhffv8v236S8d6n7N0qn751b1v1T9W999Z1vwcYBVKZ8jTIgTH22xQUlM78XEUdQBkSYmql9CfCv1b/1VblrK7rP5H5MZn5y5JldpJbOvIkkKA2FKCOA/azLJlszPehYB7wyM98wxHzcQPsoLFP+a7YvrDLzOEpA1W20tmWU4HVLugeQa2ifon6+WUYJ6noFkOcCB2bmQKOwZ+brKaPLdJoqeHM6B5AXA/fPzO8Mkua4VSPZP4xSAFjfYZOVdA8gzwTuNeFTyw9dZv6AUuCczTE87yOWmV8GnkO519VtQrlH1gPbXwMH9zsSlCRJWlyMNQFjzSnGmjMw5hm9zDyf0qDuV1022ZrpHeN+QGnwNtC9ohox+aGU2fn+3mWzTSkjevbqlHcJ5ZpeEIw/JWnyWHYHLLtPsew+A8vuzVlPNB6ZeWNmPoPSiK1bY9aVlPijV6e8q4A/DTl788qkPBfHITNfDTyB0hC6k5V0jtGmdPo9zitVGeNRlA4U/ZQ3EjgkM79FaRTY6qoO2w/TUymNUWfrTcC7u3y2ObAN7eWhy4GHdJmBb8GpZsQ7GPhsl002o3t56juUZ84fR5E3SeM1KWUCY+WRM1YeAmM2zdZiqTvKzHOBe1ANHNrBEsos3d065c37GAygiqcOBH7ax+brKJ3IH149m8cWh2Xm/1LisGtmeZwrgAcAF3X4eCnlOVTvlPdJyqA59Vn5FqTM/AAlNu/0fmIp5XfRqVPe1cBhmfm6Dp9pSOyYJzWQmZdRpoZ+BaWxTzdrKVN53zUz3zyCrNRHH9/A7EeSAG5+QN6RMqXvycD1fex2NeXvPQLYJTNfO4y8jNC7gNdTpvnttyBwHfB1ykg3d8jMfgo6XVUB+x2Bj9H7xfVllNEE9s3MU2eT5rhl5trMfD5lStxP0/tcrwd+AjyZ8rtZrC9gj2KWU0Z73kcvMz9IGenra/T+vs6njOJxl8xc1JWzkiSpN2PNrow1GzLWnMaYZwDVqLJ3ocQz5/fY9HxKJfk/VhUls0lzQ2a+CdgDeA2lgUI/8fHvgPdTRuzfLTPPmE0+Jo3xpyRNHsvuXVl2b8iy+zSW3QvricYkMz9MmSn8xcDplAZsM7kAOJEyG9bOmfm1kWVwnpig5+LIZeZnKLOBvZzSsHUmFwEfBg6qOt7Oe5l5U9Vg/jaUGPT7lN/FTcBqygxon6HMrnfHzPxJteuOtUNdOeJ8XpmZD6J0oHgv8CPK9Xk9De6x1buKF1Oe/b/sselqyruJ2y+WTnlTMvPvmfk44CDKu4vVPTZfQ5md8CGZef+cvzP9SOpggsoExsqzY6w8BsZsmq3FUneUmVdk5qHA3Skztc50X1oL/BB4IeVevSBk5s+Be1FmkjsBOIfSyW4tZXCQ04A3ArfNzCMzc01E1Dv5rs7MG0ecz08Dt6bM3vhF4LdVPjt1Iu11nLOAA4Dj6P2s/RXw2Mx88mw74M83mfnflHP9GuAPM2z+F+DtwK0z8yOjzttit2TDhlm915UWrYhYBtyp+rcDpYfx1ZQGQT/OzFn1/J4h7eOA57Ws+lb1YnEUaW0C3JVSMbE9ZdSvGyjTyF5AGensT5k5L0cYiIglwF6Uh9QtKX/jppTKl6soL4XPAc7OzH4qZAbJw2aUglNQequvo/RmPws4Y76e27qIWE4ZxWJvysv3TSgjUV4C/G9mXj6H2VuwPO+jFRG3oFQ07EEZZeRG4ELgl5k5n0e6kiRJc8RY01hziHkw1jTmGYqIuCtwe2BnykBvlwC/oMQ9I3u5HBE7UCrbdqT8hlZRZtS7itLY79zF9P0af0rS5LHsbtl9iHmw7G7ZfWQ87/2JiK0o52kXyn1gc0oj46uBPwKZmQtmhu5RmMvn4lyIiJ3ZGLPuUK2eejb+JjN/P1d5mzQRcTJwv5ZVD68aEs4rEbEXZbaKnSj3iKspMwP9KDP76Zyx4EXEKso52h24BeU+cBnl/cVpmXntHGZP0pgYKxsrDzEPxsrGbJrBYqo7qn4rdwVuRXm+bEXpuHU55RlzluXNIiIOonRSnHJGZt51rvIzqKqD4UGUwWG2o3RWvhD4aWb+di7zNkki4taUzow7Usoj11DisF8vtPvApLNjnjTPRMTmlAdLa2/2R2bmV+YoS5IkSZKkec5YU5IkSZofLLtLkqT5ohp05wLKYDtTdsrMS+coS5KkBcpYWZKkIiLeDhzZsur4zDxirvIjLRZL5zoDkho7jPYA8k/AvBtNTJIkSZI0UQ7DWFOSJEmaDw7DsrskSZof/o32Tnm/slOeJGlEDsNYWZK0yFWzuz+rtvrkuciLtNjYMU+aR6rpuV9RW/3ehTIttyRJkiRp/Iw1JUmSpPnBsrsk/yTHHwAAIABJREFUSZoLEbFygH0OBV5WW/2+4eRIkqSNjJUlSQtR0zgsIrYEPgts1bL6L4Czx0pjsHyuMyCpPxGxDDgO2K1l9aXAB+cmR5IkSZKk+c5YU5IkSZofLLtLkqQ59NyIeASlY903M/O6bhtGxE6UmfJeCixp+eh84OMjzaUkadExVpYkLWBnR8RngE9k5m96bRgRDwLeBexb++jNmbluVBmUtJEd86QJFRFPBrak/E53Ax4N3La22esz8+/jzpskSZIkaX4y1pQkSZLmB8vukiRpgiwB/rH6d0NE/BT4FXAhcC2wGbAjcDfgnkB9Zoe1wJN6deiTJKkfxsqSpEVkB+Ao4KiI+CNwBnAucCUlxtoGuA1wH2DPDvt/HWctl8ZmyYYNG7p+uOayP3T/UNJIPeCxh3HRJZd2/fyud7ojJ7znzSxdunSMuZIkSerfilvcesnMW2kxWLZiV2NLaULsufvmrFjRPY5cvXotF168eow5kiRJ6m3dmouMLQUsvtjSsrskSZoUW2+1glvssGqgfdev38All67m+uudpEFzy9hSUxZbbLnQGCtLkhaLvffcgmXLBivCXnfdWi65dDU9uglJGlC32NIZ86R56PaxD/957NF2ypMkSZIkDc0NN67j4r9aWSlJkiRNOsvukiRpnNau28D69RtYurT/RqEbNmzguuvXcvkVN7FmzfoR5k6SpMJYWZK0kNy0Zj2bLlvWaJ+1a9dz5VU3cfU1a0aUK0nd2DFPmgeWLl3KVltuwb773IoH/dN9eNRD7s/y5c0etpIkSZIktdqwoTSoufHG9fz9urVcc60vZyVJkqRJZNldkiTNpeuuW8ufrv87m65axqpVy1i5chkrli9h2bKlTI0nvX79Btatg5vWrOOGG9Zx3fVrWbPG6RkkSaNjrCxJWsguvOh6VixfwqabLmfVqmVssmIpy5cvYenSJSxZAhs2wLr1G1i3dgM33LiO1atLHCZpbizZ0GOOyjWX/cE3JJIkSZIGsuIWt+5/6FQtaMtW7GpsKUmSJGkg69ZcZGwpwNhSkiRJ0uCMLTXF2FKSJEnSoLrFlkvHnRFJkiRJkiRJkiRJkiRJkiRJkiRJkuYzO+ZJkiRJkiRJkiRJkiRJkiRJkiRJktSAHfMkSZIkSZIkSZIkSZIkSZIkSZIkSWrAjnmSJEmSJEmSJEmSJEmSJEmSJEmSJDVgxzxJkiRJkiRJkiRJkiRJkiRJkiRJkhqwY54kSZIkSZIkSZIkSZIkSZIkSZIkSQ3YMU+SJEmSJEmSJEmSJEmSJEmSJEmSpAbsmCdJkiRJkiRJkiRJkiRJkiRJkiRJUgN2zJMkSZIkSZIkSZIkSZIkSZIkSZIkqQE75kmSJEmSJEmSJEmSJEmSJEmSJEmS1IAd8yRJkiRJkiRJkiRJkiRJkiRJkiRJasCOeZIkSZIkSZIkSZIkSZIkSZIkSZIkNWDHPEmSJEmSJEmSJEmSJEmSJEmSJEmSGrBjniRJkiRJkiRJkiRJkiRJkiRJkiRJDdgxT5IkSZIkSZIkSZIkSZIkSZIkSZKkBuyYJ0mSJEmSJEmSJEmSJEmSJEmSJElSA3bMkyRJkiRJkiRJkiRJkiRJkiRJkiSpATvmSZIkSZIkSZIkSZIkSZIkSZIkSZLUgB3zJEmSJEmSJEmSJEmSJEmSJEmSJElqwI55kiRJkiRJkiRJkiRJkiRJkiRJkiQ1YMc8SZIkSZIkSZIkSZIkSZIkSZIkSZIasGOeJEmSJEmSJEmSJEmSJEmSJEmSJEkN2DFPkiRJkiRJkiRJkiRJkiRJkiRJkqQG7JgnSZIkSZIkSZIkSZIkSZIkSZIkSVIDdsyTJEmSJEmSJEmSJEmSJEmSJEmSJKkBO+ZJkiRJkiRJkiRJkiRJkiRJkiRJktSAHfMkSZIkSZIkSZIkSZIkSZIkSZIkSWrAjnmSJEmSJEmSJEmSJEmSJEmSJEmSJDVgxzxJkiRJkiRJkiRJkiRJkiRJkiRJkhqwY54kSdL/Z9+OBQAAAAAG+VuPYX95BAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAAAAAAIOYBwAAAAAAAAAAAAAAAACDmAcAAAAAAAAAAAAAAAAAg5gHAAAAAAAAAAAAxN4d5LaNbFEAlREa0AYy0TTr0Oq1jkw1yQYEGIH/IEh+Rd1k81Ikq0ieM+put+0n0gZ8Vbx4AAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAgGIeAAAAAAAAAAAAAAAAAAQU8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAgGIeAAAAAAAAAAAAAAAAAAQU8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAgGIeAAAAAAAAAAAAAAAAAAQU8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAgGIeAAAAAAAAAAAAAAAAAAQU8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAgGIeAAAAAAAAAAAAAAAAAAQU8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAACeJ4CyAAAdcElEQVQAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAQDEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAINDVHgCAfThfrrVHGPS432qPAAAAwH94zm6tZ00AAAAAAABg/1p7BtU5KkA7bMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEutoDALAtW11/3Td3a+vFAQAAjqDMYmVe22rmBAAAAAAAALZtS8+T9s3qvBVgfTbmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAl3tAQBoz5FWWQ+91i2tJQcAAGjNUKbae+4sX/veXysAAAAAAABs1d6eE31+Pc4qAZZnYx4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAAAAga72AAC04Xy51h6hOX3X5HG/rTwJAADANpR56cg5c0yePPL1AQAAaM3c538yHwAAQDvWeuZz7iw4x9x9X0NuBZiPjXkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAIBAV3sAAOqxinqa5+u21ppzAACAFpWZSM4cVl6f5yzp2gEAAKxrjjM+WQ4AAKBNaz3XOUcurPEMqjNegPnYmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAh0tQcAYD3WTS+j77rWWC8OAACwtOesI2tO83zd+jKk6wsAAPDvnMUBAABQmpITncU5/wV4lY15AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAQFd7AADYq6F13lNWpgMAALBfQxkSAACANs/X5pxpyVw4NKc8CgAAHEmLGWjOmebIqX1fo8VrB9ACG/MAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACLx9fn72fvDjx/f+DwKwCefLtfYIBB73W+0RAGbz/vXbW+0ZaMOX94tsCbBxZVaRMwFY08+Pu2zJ6XSSLQGO6shnZzXyt/wP7JVsyW+yJcC+jM2MR843NXL1ka83sG992dLGPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQKCrPQAA8H9DK7xrrBQHAAA4nYazCgAAAMzJmdgv5XWokcuH7oP3CQAAgFrGZka55ZfyOsjbAMuwMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABDoag8AUMPYFdVbXdtsBfc+TbmvW/0ZBgCALRj79/ZWM1r5+rb6GgAAANgGZ1ptGMr/5T2acr+8twAAACxN7mjDUGZ0j4A9sjEPAAAAAAAAAAAAAAAAAAKKeQAAAAAAAAAAAAAAAAAQ6GoPALCUOdYd11iZPLTCuY/VzvQpfzam/GwBAMDRzfF3dI2/xafkxOc5ZU0AAADgFd5bAAAA5uDZx2Pru/8yJ9AKG/MAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACHS1BwDgb+fLtfYIAAAAbNzjfqs9AgAAAPSSWzN912vus+Wx96Xv+7qvAADAmjxvm1kry429L2tlXYCl2ZgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIdLUHAFjK84pjq405uuffgbnXjwMAwB75OxoAAADgb33vjwy9j+K8HgAAgFJfThx6/rv8mOfEgVbYmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAh0tQcAAAAAAOo7X661RwAAAOAgnjPo436rNAml8r7McU+81wAAAEwxNo88/38ySBvK+zLHPfGeAdA6G/MAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACHS1BwBYy+N+iz/nfLkuMAkAAABbNSUnTsmjaylfz/OcMjEAAABrGcqnLGvO6+29BQAAYKoyTzxnibG5ZehrsKw5r/fU+w9Qi415AAAAAAAAAAAAAAAAABBQzAMAAAAAAAAAAAAAAACAQFd7AICW9a0/tuIaAACAsfoyZF/mXFMLMwAAAADjTTmrHsr/r559e28BAACYg+dy1zUlyw3do1ezofsPbJmNeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgEBXewCALRpauWydMgAAAGMM5ceh3AkAAAAcV/meQQvvLTgfBwAA2J4yy7XwTLTzcWDLbMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEutoDAAB1jF1FDgAAHEOZEQAAAIDtGXvm1/cewKufDwAAMAfPM9Y1NvP13adXP3+q8uvJrcCabMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAINDVHgAAAAAAAAAAAMg87rfaIwAAALBx58u19ggAm2ZjHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAINDVHgCgVK5DftxvFSeB/fM7BgDAXpV/65Y5k2HPGcG1AwAAYM+Gcu9Rz9G8FwAAALTgOZu0nNGGZjtqxqp1v456vYH6bMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEutoDAMc2tDZ4ykrhFtZV981gRTIAAMAyhrLglJzYQn7rm6GF3AsAAABjlTm2hbxd2mPGbu0aAwAAzKHMOq1luT3msNauMUDrbMwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAINDVHgA4nvPlWnuEKh7321//ftTrAAAAMIfnjHUUz1nyqNcBAAAAcOYMAADA65w5A7zGxjwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAECgqz0AcAzny3Xz33PuVc3l16txfQAAALZm7lxW43vOnf/KrzdlVnkUAACAOT1n0ym589Wsu3eyPAAAsEdDz9ROyYae0R0mbwPMx8Y8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAoKs9ALAfLax6XnKG8mtb4cwW+bkFAGALWvi7dckZyq/dQo4GAACAlg2d0fbl6hbeWwAAAKAdQ2e0fRnSWW773COgFTbmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAl3tAYC6rPEFAADgVY/7rfYIAAAAQAOWPH8+8tn2kV87AABwPEuePx/5bPvIrx1gSTbmAQAAAAAAAAAAAAAAAEBAMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABDoag8ALO98udYeAajkcb/VHgEAgJ3wtyUAAADQEu9VAAAA8CrPWG+L+wW0yMY8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAoKs9AMAWjV2F/LjfFp4E/snPHQAAbMPYv93HZlAAAADYmyOfe3k/AAAAYB5Hzld7yNVHvn/ANtiYBwAAAAAAAAAAAAAAAAABxTwAAAAAAAAAAAAAAAAACHS1BwCWN3YNsVW/8yuv6dB9KD/mPgAAAC0am1XGZlDGG5sZx2ZQAAAAaNlRM61zYgAAgPkcNWMdNVMD1GRjHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAACBrvYAQDse91v8OefLdYFJgMSU310AAFjKlJzob1oAAAAAAAAAwLPpwNbYmAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAh0tQcAtu1xv/3557lXB5dfew61Vxs/f/+5Xx8AAMBWlXmp9SxYO8s9f//aWRcAAAA4tiWfGQAAAOAY5Elgy2zMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABLraAwD78bjf/vr3MWuFnz9nSeX3annl8dA1aXlu1rXm7w4AAKzpOfeM+dt3zaxUfq+W/y4fuiYtzw0AAAC0bytn7wAAALRLngT2wsY8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAoKs9ALBfj/vtzz+X64bL/15L3wxLrkWe+rqtaj6WFn4/AACgJX15soWs1DfDkn/XT33dfTO1dk0BAACA9nkPAQAAmFtrzxkzD/kROAIb8wAAAAAAAAAAAAAAAAAgoJgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIdLUHAI7hcb/VHmE15Ws9X66jPmfs/8e2Hen3AAAAlnCk7FS+1rFZQuYAAAAAAAAAYA1HOr8HGGJjHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAINDVHgBgKx73W+0R2AA/JwAAwJDz5Vp7hF7PeablWQEAAGCr5G0AAIBhz7nJc5n/tNY1kWEB/puNeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgEBXewCAllh3zRh+TgAAgCHny7X2CAAAALBZZa7e6rncHl4DAADAlpVZbKvnt2u9hq1eH4BW2JgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQKCrPQAAtOpxv9UeAQAAoJoyE50v14qTAAAAcFRlHm3x7G5MXpapAQAAXvNqNmz93HPMa5ojE7f42gH2wMY8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABAoKs9AADtr8k+kjnWfQMAANRQ5sm5s43cCgAAQG1DeXRsDp6SneVgAACAdrx6Jjr0OWPz35Sz01rPpsq0AMuzMQ8AAAAAAAAAAAAAAAAAAop5AAAAAAAAAAAAAAAAABB4+/z87P3gx4/v/R8EYBHWRi+j1hpwgCN7//rtrfYMtOHL+0W2BFjZmhlIjgVgST8/7rIlp9NJtgQAAKaTLflNtgTYtyM9J+qMFmB9fdnSxjwAAAAAAAAAAAAAAAAACCjmAQAAAAAAAAAAAAAAAEDg7fOzfzP3x4/v1nYDVGbd9LAjrR4H2Jr3r9/+dW03x/Pl/SJbAlRWIzvJswDM4efHXbbkdDrJlgAAwHSyJb/JlgDHtKXnTJ2xArSrL1vamAcAAAAAAAAAAAAAAAAAAcU8AAAAAAAAAAAAAAAAAAgo5gEAAAAAAAAAAAAAAABA4O3z87P3gx8/vvd/EIDVnS/X2iMs6nG/1R4BgBm9f/32VnsG2vDl/SJbAjRk79lr79kZ4Gh+ftxlS06nk2wJAABMJ1vym2wJAABM1ZctbcwDAAAAAAAAAAAAAAAAgIBiHgAAAAAAAAAAAAAAAAAEutoDADDe437788/ny7XiJL8MzVN+DAAAgHaU+a2F7DY0TwvZFwAAAAAAAAAA/o2NeQAAAAAAAAAAAAAAAAAQUMwDAAAAAAAAAAAAAAAAgEBXewAApnncb70fO1+uK07yy9A8AAAAtGkoP9bIeTXyLAAAAAAAAAAATGFjHgAAAAAAAAAAAAAAAAAEFPMAAAAAAAAAAAAAAAAAIKCYBwAAAAAAAAAAAAAAAACBrvYAAMzvcb/VHgEAAICNO1+utUcAAAAAAAAAAIBm2ZgHAAAAAAAAAAAAAAAAAAHFPAAAAAAAAAAAAAAAAAAIKOYBAAAAAAAAAAAAAAAAQEAxDwAAAAAAAAAAAAAAAAACinkAAAAAAAAAAAAAAAAAEFDMAwAAAAAAAAAAAAAAAICAYh4AAAAAAAAAAAAAAAAABBTzAAAAAAAAAAAAAAAAACCgmAcAAAAAAAAAAAAAAMD/2ruDEgBgAAZizL/oTsO9xiBRUANHAQiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACAQJgHAAAAAAAAAAAAAAAAAIEwDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEAgzAMAAAAAAAAAAAAAAACA4Gx7vQEAAAAAAAAAAAAAAAAAvuExDwAAAAAAAAAAAAAAAAACYR4AAAAAAAAAAAAAAAAABMI8AAAAAAAAAAAAAAAAAAiEeQAAAAAAAAAAAAAAAAAQCPMAAAAAAAAAAAAAAAAAIBDmAQAAAAAAAAAAAAAAAEBwAcIrqR4IjxWaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 3600x2160 with 4 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Define a single layer for plotting\n",
    "layer = 50\n",
    "# Define a dictionary of class labels\n",
    "classes_dict = {\n",
    "    'Normal': 0.,\n",
    "    'Edema': 1.,\n",
    "    'Non-enhancing tumor': 2.,\n",
    "    'Enhancing tumor': 3. \n",
    "}\n",
    "# Set up for plotting\n",
    "fig, ax = plt.subplots(nrows=1, ncols=4, figsize=(50, 30))\n",
    "for i in range(4):\n",
    "    img_label_str = list(classes_dict.keys())[i]\n",
    "    img = label_array[:,:,layer]\n",
    "    mask = np.where(img == classes_dict[img_label_str], 255, 0)\n",
    "    ax[i].imshow(mask)\n",
    "    ax[i].set_title(f\"Layer {layer} for {img_label_str}\", fontsize=45)\n",
    "    ax[i].axis('off')\n",
    "plt.tight_layout()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### 2.2.1 Interactive Visualization Across Layers\n",
    "As another way of looking at the data, run the code below to create a visualization where you can choose the class you want to look at by clicking a button to choose a particular label and scrolling across layers using the slider!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "e85cb9f650a149579d400dfedecba2c8",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "interactive(children=(ToggleButtons(button_style='info', description='Select Class:', options=('Normal', 'Edem…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Create button values\n",
    "select_class = ToggleButtons(\n",
    "    options=['Normal','Edema', 'Non-enhancing tumor', 'Enhancing tumor'],\n",
    "    description='Select Class:',\n",
    "    disabled=False,\n",
    "    button_style='info', \n",
    "    \n",
    ")\n",
    "# Create layer slider\n",
    "select_layer = IntSlider(min=0, max=154, description='Select Layer', continuous_update=False)\n",
    "\n",
    "    \n",
    "# Define a function for plotting images\n",
    "def plot_image(seg_class, layer):\n",
    "    print(f\"Plotting {layer} Layer Label: {seg_class}\")\n",
    "    img_label = classes_dict[seg_class]\n",
    "    mask = np.where(label_array[:,:,layer] == img_label, 255, 0)\n",
    "    plt.figure(figsize=(10,5))\n",
    "    plt.imshow(mask, cmap='gray')\n",
    "    plt.axis('off');\n",
    "\n",
    "# Use the interactive() tool to create the visualization\n",
    "interactive(plot_image, seg_class=select_class, layer=select_layer)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### And there you have it! We hope this lab has helped you get a better sense of the data you'll be working with in this week's assignment. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}