{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "bf6256ae",
   "metadata": {},
   "source": [
    "## 12/08/2022 Coverslip4\n",
    "## Registration with QC metrics\n",
    "### starts after tif conversion"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "c1337f9f",
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.patches as patches\n",
    "from pathlib import Path\n",
    "import numpy as np\n",
    "import os\n",
    "import sys\n",
    "import glob\n",
    "import pandas as pd\n",
    "import xml.etree.ElementTree as et\n",
    "import datetime\n",
    "from skimage.io import imread, imsave\n",
    "from imageio import volread as imread\n",
    "\n",
    "import tifffile\n",
    "import pystackreg\n",
    "from pystackreg import StackReg\n",
    "from skimage.filters import threshold_otsu\n",
    "\n",
    "from pystackreg.util import to_uint16   # make sure version 0.2.5 (not anything below)\n",
    "#from ims_to_tiff import convert_to_tif   "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "4097d2c8",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/mnt/disks/store/102222_D10_Coverslip4_Reimage_Processed\n"
     ]
    }
   ],
   "source": [
    "print(os.getcwd())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "9adc6d6d",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define input params, directories\n",
    "allfiles = sorted(glob.glob('tif/*'))\n",
    "allCycles = []\n",
    "for x in allfiles:\n",
    "    if x.startswith('tif/Cycle_'):\n",
    "        allCycles.append(x.split('_')[-1])\n",
    "\n",
    "_allFOVs = sorted(glob.glob('tif/Cycle_0/*'))\n",
    "allFOVs = []\n",
    "for x in _allFOVs:\n",
    "    allFOVs.append(x.split('F')[-1][0:3])\n",
    "    \n",
    "    \n",
    "for i in allCycles[1:]:\n",
    "    if i == '0':\n",
    "        continue\n",
    "    try:\n",
    "        os.makedirs(f'tmat_Cyc_{i}')\n",
    "        os.makedirs(f'reg_bin_Cyc_{i}')\n",
    "        os.makedirs(f'reg_Cyc_{i}')   \n",
    "    except: \n",
    "        continue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2c700364",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "210"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(allFOVs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "8f573a3c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['0', '8', '9']"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "allCycles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f539d4a9",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "cycle 8   FOV 000 \n",
      "Saving Registered Images... reg_Cyc_8/Cycle_8_F000_reg.tif\n",
      "cycle 8   FOV 001 \n"
     ]
    }
   ],
   "source": [
    "# # loop through registration\n",
    "for c in allCycles[1:]:    \n",
    "    for fov in allFOVs:\n",
    "        # read images\n",
    "        ref_name = f'tif/Cycle_0/Cycle_F{fov}.tif' \n",
    "        mov_name = f'tif/Cycle_{c}/Cycle_{c}_F{fov}.tif'\n",
    "        ref = imread(ref_name).astype(np.uint16)\n",
    "        mov = imread(mov_name).astype(np.uint16)\n",
    "        \n",
    "        print(f'cycle {c}   FOV {fov} ')\n",
    "        \n",
    "        ref_max = ref.max(0)\n",
    "        ref_binary = ref_max[0] > threshold_otsu(ref_max[0]) # nuclei\n",
    "        mov_max = mov.max(0)\n",
    "        mov_binary = mov_max[0] > threshold_otsu(mov_max[0]) # nuclei\n",
    "        \n",
    "        sr = StackReg(StackReg.RIGID_BODY) # create instance of stackreg class\n",
    "        tmat = sr.register(ref_binary, mov_binary) # calculate tmat\n",
    "        out = sr.transform(mov_binary)\n",
    "        out = pystackreg.util.to_uint16(out) # make sure to convert later into bool. Save first then convert.\n",
    "        \n",
    "        # save binary images\n",
    "        fname_to_save = f'reg_bin_Cyc_{c}/Cycle_{c}_F{fov}_bin_reg.tif'\n",
    "        tifffile.imwrite(fname_to_save, out, imagej=True, photometric = 'minisblack', metadata={'axes':'YX'})\n",
    "        \n",
    "        # save tmat\n",
    "        np.save(f'tmat_Cyc_{c}/Cycle_{c}_F{fov}_tmat.npy', tmat)\n",
    "        \n",
    "        # REGISTRATION\n",
    "        reg = np.zeros(mov.shape, dtype=np.uint16)\n",
    "        for Z in range(mov.shape[0]):\n",
    "            for ch in range(mov.shape[1]): \n",
    "                reg[Z,ch,...] = sr.transform(mov[Z,ch,...], tmat=tmat)\n",
    "                reg[Z,ch,...] = to_uint16(reg[Z,ch,...])\n",
    "\n",
    "        fname_to_save = f'reg_Cyc_{c}' + f'/Cycle_{c}_F{fov}_reg.tif'\n",
    "        print('Saving Registered Images...', fname_to_save)\n",
    "        \n",
    "        # save registered images\n",
    "        tifffile.imwrite(fname_to_save, reg, imagej=True,\n",
    "                         photometric = 'minisblack',metadata={'axes':'ZCYX'})\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "65938805",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e99d07b8",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "219a37ed",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJIAAAI4CAYAAAAxusOtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABn9klEQVR4nO3dd3xkVf3/8ddJtvddlmULyy5dijSl2hBEkSKgICgKKorY/epX5WvvYgW7YqNYEFEQhB+CFBGkg4CC9LIsuyzbKMvW5Pz+ODebSTaT3CRT7kxeTx6Xmblz596TTDbzzueee06IMSJJkiRJkiT1paXeDZAkSZIkSVJjsJAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSVovhPCeEMKTIYTnQggb1bs9A5W1f4s6HftlIYR763FsSZLUOEIILwkh3J/llsPr3Z7+qHfeCSF8MoTw83odXxrqQoyx3m2QVCEhhEeATYA2YC3wT+CkGOO8HK8dDjwD7BVjvKOa7cyOdzWwF7AOWAVcA7wvxrig2sfOjr8v8OsY46aD3E8Eto4xPlCJdkmSpOaSZZ6dgekxxtUl668ALowxfjd7XPFMke3zeSACTwO/Bz4WY2yr1DH6OP7nga1ijG8ZxD72pQKZTVLl2CNJaj6HxhjHATOAJ4Hv53zdJsAo4D/9PWBIBvL75P1ZW7cCxgHfGsA+qiaEMKzebZAkSY0rhDAXeBmpkPO6bk/PYQC5q8xxesssO2d56xXA0cA7KnHMShhEhpRUR/6jlZpUjHEVcB6wfce6EMLIEMK3QgiPZZew/SSEMDqEsA3Q0T15eQjhymz7fUIIN4cQns5u9ynZ19UhhK+EEK4jnenaIoTwghDC5SGEpSGEe0MIb8zZ1uXABcAuJfsvu68QwkYhhItCCM9k7fpyCOHakudjCGGr7P5BIYS7QwjPhhDmhxD+N4QwFvh/wMysO/lzIYSZIYTPhxDOCyH8OoTwDPC2EMIeIYTrQwjLQwgLQgg/CCGMyPZ9TXbIO7J9HB1C2DeE8HhJW7bLvlfLQwj/CSG8ruS5M0IIPwwhXJy178YQwpZ5vmeSJKkhHAfcAJwBHN+xMoTwILAFcFGWIa7PnlqfKbLtDgkh/CvLEf8MIexUso9HQgifCCHcCazo6wRY1tPpOrrmrd72v1sI4fYso/whhPD7EMKXs+e6551PZDnr2Sy37R9COBD4JHB09jXdkW3bU4Z8ewjhnuz1D4UQ3p1t21tm+3XJ8V+X5azl2f636/Z9+t8Qwp1Zpv19CGFUn++cpLIsJElNKoQwhnTW6YaS1acA25ACxFbALOCzMcb7gB2ybSbFGPcLIUwBLga+B2wEfAe4OHQdO+mtwInAeOAp4HLgt8A04BjgRyGE7elDts/XAw9kj8f2sa8fAiuA6aRQdnz3fZb4BfDuGON4YEfgyhjjCuC1wBMxxnHZ8kS2/WGkAtwk4DekywT/B5gK7A3sD7wXIMb48uw1O2f7+H23r2s4cBFwWfZ1fAD4TQhh25LNjgG+AEzOvv6v9PX9kiRJDeM4Up74DfCaEMImADHGLYHHyHqSxxj3zrZfnylCCLsCvwTeTcpiPwUuDCGMLNn/m4CDSfltXW8NCSG8gNQ7qiNvld1/dtLsfFIBbArwO+CIMvvdFng/sHuWt14DPBJjvBT4KvD77GvaueRlpRnyUWARcAgwAXg7cGoIYbc+MlvH8bfJ2vdhYGPgElKBbkTJZm8EDgQ2B3YC3tbb90pS7ywkSc3nghDCctJ18AcA34TUdZj0gf0/McalMcZnSR/ux5TZz8HA/THGs2OM62KMvwP+Cxxass0ZMcb/ZMHlQFJo+FW2/e3AH4Gjemnr90IITwOLSYWaD2TrDym3rxBCK/AG4HMxxudjjHcDZ/ZyjLXA9iGECTHGZTHG23rZFuD6GOMFMcb2GOPKGOOtMcYbsnY8QgpZr+hjHx32Il2yd0qMcU2M8UrgL6TQ1+H8GONN2ffwN5ScJZQkSY0rhPBS0uVr58YYbwUeBN7cj12cCPw0xnhjjLEtxngmsJqULzp8L8Y4L8a4spf93BZCWAHcA1wN/CjH/vcChmX7Xxtj/BNwU5n9twEjSXlreIzxkRjjg318beszZLb/i2OMD8bk76STcC/rYx8djgYujjFeHmNcSxoqYTSwT8k234sxPhFjXEo6ybdLzn1L6oGFJKn5HB5jnEQa7+j9wN9DCNNJZ2jGALdm3X6XA5dm63syk3SGqNSjpF5MHUoH8Z4D7Nmx72z/x5J6DZXzwRjjRNKZoclAxyCKve1rY1KwKT12b4OJvwE4CHg0hPD3EMLevWy7wb5CCNuEEP4SQlgY0uVuXyUVvfKYCcyLMbaXrOv+PVxYcv95UuFJkiQ1vuOBy2KMi7PHv6X3XtTdzQE+2i0PzSbliw59TqgC7EbKF0cDewJjc+x/JjA/dp2ZqcdjZZfMfRj4PLAohHBOCGFmT9uW21cI4bUhhBtCGtJgOSm79Sdvrc+sWe6ah3lLqhoLSVKTys4s/Yl0luilpF4/K4EdYoyTsmViNvhiT54gBYxSmwHzSw9Tcn8e8PeSfU/Kuh+/J0db7wK+DPww6znV276eIs30Vjpzx+xe9n1zjPEw0qVlFwDn9tD2Li/p9vjHpJ5YW8cYJ5Cu9Q99fU2ZJ4DZoesgkt2/h5IkqcmEEEaTLqd6RXYyaiHpUvmdQwg79/7q9eYBX+mWh8ZkvcQ75JqCO+vpcy5wPfDZHPtfAMzKclmH3vLWb2OMHT2wIvD1Ptq3fn12qd4fST2JNslOiF5CZ97q62vsklmzNs/GvCVVjYUkqUmF5DBST597srMzPyNdcz4t22ZWCOE1ZXZxCbBNCOHNIYRhIQ36uD3p0qye/CXb/q0hhOHZsnvpYId9OJM0c9zrettXTNPV/gn4fAhhTHa9/3FlvgcjQgjHhhAmZl2dnwE6egc9CWwUQpjYR7vGZ697LjtW98LYk6TBMntyI+ms18ezr2Ff0qWB5/RxTEmS1NgOJ53M2550GdUuwHbAPyiTW9gwU/wMOCmEsGeW68aGEA4OIYwfRLtOAd6V9Vbvbf/XZ+1/f5YDDwP26GmHIYRtQwj7ZQWhVaQTl6V5a27ofWa2EaRL454C1oUQXgu8uuT5vjLbucDBIQ3wPRz4KOkSvX/29c2QNDAWkqTmc1EI4TlS8eMrwPExxo6pZT9BGmDxhuwyrb8B2/a0kxjjEtJYRR8FlgAfBw4p6Z7dfftnSR/6x5DODC0knY0a2dP2Pbx+DfBd4DM59vV+YGK2/mzSAIury+z6rcAj2dd7EukSOWKM/81e91DWnbtcF+z/JY1n8CwpcP2+2/OfB87M9tFllrrsazqUNEjkYtKYBMdlx5YkSc3reOBXMcbHYowLOxbgB8CxoecZ1j5PSaaIMd4CvCt7zTJShnvbYBqV9QK/BvhYb/vPMszrgROA5cBbSCf6espbI0kFqsWkbDYN+L/suT9kt0tCCD2OU5nlvg+SCkLLSLnrwpLne81sMcZ7s/Z9P2vDoaRBzNfk+JZIGoDQ9bJXSWo8IYSvA9NjjP0Zd0CSJEk5hRBuBH4SY/xVvdsiqb7skSSp4YQQXhBC2Cnrhr0H6WzZ+fVulyRJUrMIIbwihDA9u7TteNLkKJfWu12S6q+nLpWSVHTjSV2cZ5Kum/828Oe6tkiSJKm5bEu63Gws8BBwZIxxQX2bJKkICnVpWwjhQNIYKa3Az2OMp9S5SZIkSU3PDCZJkvIqTCEphNAK3AccADwO3Ay8KcZ4d10bJkmS1MTMYJIkqT+KNEbSHsADMcaHshH2zwEOq3ObJEmSmp0ZTJIk5VakMZJmAfNKHj8O7Nl9oxDCicCJ2cMX1aBdkiSpjmKMod5taHJmMEmStIFyGaxIhaRcYoynA6cDhBCKcV2eJElSkzODSZIkKNalbfOB2SWPN83WSZIkqXrMYJIkKbciFZJuBrYOIWweQhgBHANcWOc2SZIkNTszmCRJyq0wl7bFGNeFEN4P/JU09ewvY4z/qXOzJEmSmpoZTJIk9UeIsXEvcff6fEmSmp+DbRePGUySpOZXLoMV6dI2SZIkSZIkFZiFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLlUrJIUQfhlCWBRC+HfJuikhhMtDCPdnt5Oz9SGE8L0QwgMhhDtDCLtVq12SJEnNzAwmSZKqqZo9ks4ADuy27mTgihjj1sAV2WOA1wJbZ8uJwI+r2C5JkqRmdgZmMEmSVCVVKyTFGK8BlnZbfRhwZnb/TODwkvVnxeQGYFIIYUa12iZJktSszGCSJKmaaj1G0iYxxgXZ/YXAJtn9WcC8ku0ez9ZtIIRwYgjhlhDCLdVrpiRJUlMxg0mSpIoYVq8DxxhjCCEO4HWnA6cDDOT1kiRJQ5kZTJIkDUateyQ92dFdOrtdlK2fD8wu2W7TbJ0kSZIGzwwmSZIqotaFpAuB47P7xwN/Lll/XDZzyF7A0yXdryVJkjQ4ZjBJklQRIcbq9EwOIfwO2BeYCjwJfA64ADgX2Ax4FHhjjHFpCCEAPyDNMPI88PYYY5/X39utWpKk5hdjDPVuQyMxg0mSpEool8GqVkiqBUOMJEnNz0JS8ZjBJElqfuUyWK0vbZMkSZIkSVKDspAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJyqVohKYQwO4RwVQjh7hDCf0IIH8rWTwkhXB5CuD+7nZytDyGE74UQHggh3BlC2K1abZMkSWpWZjBJklRN1eyRtA74aIxxe2Av4H0hhO2Bk4ErYoxbA1dkjwFeC2ydLScCP65i2yRJkpqVGUySJFVN1QpJMcYFMcbbsvvPAvcAs4DDgDOzzc4EDs/uHwacFZMbgEkhhBnVap8kSVIzMoNJkqRqqskYSSGEucCuwI3AJjHGBdlTC4FNsvuzgHklL3s8W9d9XyeGEG4JIdxSvRZLkiQ1PjOYJEmqtKoXkkII44A/Ah+OMT5T+lyMMQKxP/uLMZ4eY3xxjPHFFWymJElSUzGDSZKkaqhqISmEMJwUYH4TY/xTtvrJju7S2e2ibP18YHbJyzfN1kmSJKkfzGCSJKlaqjlrWwB+AdwTY/xOyVMXAsdn948H/lyy/rhs5pC9gKdLul9LkiQpBzOYJEmqppB6NldhxyG8FPgHcBfQnq3+JOka/XOBzYBHgTfGGJdmoecHwIHA88DbY4y9XoMfQqhO4yVJUmHEGEO929BIzGCSJKkSymWwqhWSasEQI0lS87OQVDxmMEmSml+5DFaTWdskSZIkSZLU+CwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXHIXkkIIU0IIU6rZGEmSJHVlBpMkSUXSayEphLBZCOGcEMJTwI3ATSGERdm6uTVpoSRJ0hBjBpMkSUXVV4+k3wPnA9NjjFvHGLcCZgAXAOdUuW2SJElDlRlMkiQVUogxln8yhPtjjFv397laCSGUb7wkSWoKMcZQ7zbUmhlMkiTVW7kMNqyP190aQvgRcCYwL1s3GzgeuL1yzZMkSVIJM5gkSSqkvnokjQBOAA4DZmWr5wMXAr+IMa6uegt74dkwSZKa3xDtkWQGkyRJdVUug/VaSCo6Q4wkSc1vKBaSis4MJklS8xvQpW0hhGGks2GH0/Vs2J9JZ8PWVrCNkiRJwgwmSZKKq69L234HLCddn/94tnpT0vX5U2KMR1e7gb3xbJgkSc1vKPZIMoNJkqR6G9ClbSGE+2KM2/T3uVoxxEiS1PyGaCHJDCZJkuqqXAZr6eN1S0MIR4UQ1m8XQmgJIRwNLKtkAyVJkrSeGUySJBVSX4WkY4AjgSdDCPeFEO4HFgKvz56TJElS5ZnBJElSIeWetS2EsBFAjHFJVVvUD3arliSp+Q3FS9tKmcEkSVI9DPTSttIdLAEmhBBeH0J4QcVaJkmSpLLMYJIkqUh6LSSFEC4ouX8YcCVwKHBhCOFtVW2ZJEnSEGUGkyRJRTWsj+fnlNz/BLBfjPHhEMJU4ArgjGo1TJIkaQgzg0mSpELq69K20uvfh8UYHwaIMS4G2qvWKkmSpKHNDCZJkgqprx5JO4cQngECMDKEMCPGuCCEMAJorX7zJEmShiQzmCRJKqReC0kxxnJBZQzw7so3R5IkSWYwSZJUVCHGxp291alnJUlqfuWmnlX9mMEkSWp+5TJYX2MklRVCuGvgzZEkSdJAmMEkSVI99XppWwjh9eWeAqZXvjmSJEkyg0mSpKLqa7Dt3wO/oevMIR1GVb45kiRJwgwmSZIKqtcxkkIItwLHxxj/3cNz82KMs6vZuL54fb4kSc1vKI6RZAaTJEn1NtAxkj4MPFPmuSMG0yBJkiSV9WHMYJIkqYCctU2SJBXaUOyRVHRmMEmSmt+gZ20LIRzS22NJkiRVnhlMkiQVSe5CEvDibo93r2RDJEmS1CMzmCRJKow+L20LIbwAOAzYlDRzyHzgwhjjPdVvXu/sVi1JUvMbqpe2mcEkSVI9DejSthDCJ4BzgADcCNyU3f9dCOHkSjdSkiRJZjBJklRcvfZICiHcB+wQY1zbbf0I4D8xxq2r3L5eeTZMkqTmNxR7JJnBJElSvQ10sO12YGYP62dkz0mSJKnyzGCSJKmQhvXx/IeBK0II9wPzsnWbAVsB769iuyRJkoayD2MGkyRJBZRnsO0WYA9gVrZqPnBzjLGtym3rk92qJUlqfkPx0jYwg0mSpPoql8H6LCQVmSFGkqTmN1QLSUVmBpMkqfkNdIwkSZIawuXA6Ho3QpIkaSi4mvUD5Vx+OYz+BzAOuBiYXLdWqUbskSRJaniXAt8GPgKMLFl/FPCDbP2COrRLlWGPpOIxg0nSELcGuA4ujfDtl8NHWmHkP4C9gJvgqMPgBz+Aj2wCC44DHq9vczUwXtomSWoKl2a3BwEXAa3Aq4CrgH3pOovE1cCewC3A89m69wIP1aCdqhwLScVjBpOkoefSS4HXwkEXw0WvgdaWXjLY1bDnnnDLaHj+WuAt8N5H4aHvAT8B7q516zUQFpIkSU3hYOBC4BJSMam/12jvAtxR4TapuiwkFY8ZTJKGkIvSzcEHw4UXwyUHw0GhnxnsatjlObjjT8A+wHeBf1e6oao0C0mSpIb2x5L7RwADrSzsArwZ+CVwLHAecOegWqZqs5BUPGYwSRo6/ngEKYiFCmSwq+CXO8Gx/4HzPgB3GsIKrVwGG9bTSkmSiuR3wOs7HvyWVAGKwK9J17YdB/wMOAlY1fu+TgF2B3YDdgKuqUJ7JUmSmsHvfleSwQbpFGD3V2YZ7Gq45skK7Vg1ZyFJklRoZwLHlK44BlhHKiS9idSvug04GhgBvBP4JmnmkA8By7vu78Ds9lXZ7ceBp/ByN0mSJCB1OTojy2DH9LFtP3TJYHvCx+fCU9+BOwxhDcdL2yRJhRSA00l1IX4GvCu7PYHe+1SfBbwRGEXqvfQBYCnwJeDHwBMbvuSvwMnAvyrTdFWYl7YVjxlMkppXAE5/J7zzZ9U/1l//CifPg3/9lDQ7igrFMZIkSQ2lJUDbD7IH7yEVgd5D/y/MPwt4DrgfmA18B5i/4WZHABcAxwP3AdcPoM2qDgtJxWMGk6Tm1QK0BaC9Nsc74qdwwTfh+AfNYEXjGEmSpMZzP3Bqdv+9A9zHcdntGcChwCbA/wHzum72DuAVwGI8ISZJkoawQDrxVgunwzuWwCvGm8EaiYUkSVLxtABfZeDTgvTkbdntU8DaDZ8+NLv9C+lSN0mSpCFrQY2OMwUO3RfYAf7yBfjr7TU6rgbFQpIkqThagc9BaIEvfKJKxxhD6t30MzbolQRwCGkc7y8BN1epCZIkSUUUgC8AjK/RAY/MbreBQ/4A8XYzWCNoqXcDJEnqohX4VCrmfLka+z8R+Ey2fBrYfMNNDgVeVI1jS5IkFVwM8OVP1+HAh8OhO5vBGoGFJElScbQBn09FpM9mS9W8C9ieVLiSJElSZwaLwDdqfPAjsYrUICwkSZKKoxX4QI2O9XvgDmBNz0/vB+xSo6ZIkiQVTq3GSSr1KtjvhWaworOQJEkqjgBs3Xn3JOBH1TrWPNJUbR8FNtvw6WnApGodW5IkqaBChJN+DD86te9tK24TmDbRDFZ0FpIkScXRDpTM1rEr8L5qHet/gW2ADwI7bvj0OcDV1Tq2JElSUUXY9f1VzGC9eQzOWWwGKzoLSZKk4hgO/DTdjaRxsWtiP2B211V7Ai+o1fElSZIKIgInRuCsOhx8G9hzhhms6CwkSZKKpZ00flGJc7IlVuuYH2WDwR3fBry6WseTJEkquuNrkMG6uxPe9pgZrOiG1bsBkiR10Q6cBxzduepN2e0I4AjS+EkVtxdwMzAfbgUeBR6sxnEkSZIaxJv+lG5HHAFHhCplsBK37g6PHgkPnodBrMDskSRJKo524K/AH3p++g3AX6jSWbH9gc3T3WuB/wEursZxJEmSGsUvgTdUOYOVuPZH8D/PwMXbV/lAGhQLSZKk4lhL6nLUi9eRak1/JdWdKuZi4P5090PAQZXctyRJUiP6QHZ7GbwuVimDlfgQcNCPgYuqdABVhIUkSVKxRPqcquO1wIHAlVQwyLwe2K5SO5MkSWoCB3a9rUoG63A3sKDSO1U1WEiSJBXLOuBV+TY9ALiGCgWZ3wJ3pbv3Y46RJEkqp6IZrMN5cP+tZrBGYCFJklRMNwI39b3ZK4HrgLbBHu/twJHARnAK8OfB7k+SJKmJVSyD3QcsAY6FU/Y0gzUCZ22TJBVPG3AScAe5TnW9nBRkRmWPd2IAH3CfAPaCR/aEJZf098WSJElDT8Uy2DHwyEWwxLGRGoKFJElS8bQCZwC7ki43e2HfL3lJdrsd8E9St+itgOF9vO5hYCqwFHjmZPgkaVYSSZIkZXbMbv9NGs9yRyCkVQPOYA/D1Odg6SR45hn45NNmsEYRYqz2BH7VE0Jo3MZLkvLrx2/7m4DxwNuA7wBTsvVbACOy+48Aq7L7JwDHzYM/vB2uuGLwTVXlxRhDvdugrsxgkjSEbAvcA9xLml03kgpK3bql5M5gjwIr4YQT4Lh/wh++D1dcDFxazS9CA1Eug1WtkBRCGEUaf2sk6UfsvBjj50IImwPnABsBtwJvjTGuCSGMBM4CXkS6QvLoGOMjfRzDECNJQ8EAftvfCBwKLMoeX07qfj0C2B24pXTjE0ifTM8Poo2qGgtJ/WMGk/qpBdiM9BeupK62AB4g9T7akjQjSR8jLfeawR6HWw4DbqtOc1VZ5TJYNQfbXg3sF2PcGdgFODCEsBfwdeDUGONWwDJSfCe7XZatPzXbThpyArBpyePZZe5L6t2edAYY6Jxd5DFgTcfKJ4GVwNeAl9aydVJVmcGk/hgHXGwGkzawGfAg6y9h40FyVRDWZ7D5wDo44Am45jF47DFYczAWkZpATS5tCyGMAa4F3gNcDEyPMa4LIewNfD7G+JoQwl+z+9eHEIYBC4GNYy8N9GyYGlkrsEkP60eTeopuSTqN/BApvCwAHqdrwJGGjGr9tt8X+BRwGuAA24Vlj6SBM4NJJYYB07plsHbSSYUZMLod/r0QtpyZZbAIsxeYwTTEzMxuA+kHf6AWkvq5nge8hfRHjRpOuQxW1cG2QwitpK7TWwE/JNUwl8cY12WbPA7Myu7PAuZljV0XQnia1PV6cbd9ngicWM12S9XSSvqhhhRGbu1l2/kl9x8HZmCA0dAyjexs1jSYFlnfvWjRNJgW0sONgacGsvOlwFpgLLCC1H9DaiJmMKmr1lbYaA/gum4ZbAmwN2n68eeB7WD+o9lz6+Dx4TBjGmy6aINdSs1nFvB4SQYj3e+wiM7nymawp0gF2l1IRdp9qtRW1VVVC0kxxjZglxDCJOB84AUV2OfpwOng2TA1llbSuCyXAJP7+doWUlG/hTRo3dLKNk0qnFbgCWDaVAgL4cklwPbpuamL09nh6dk2Y0g1oVyWA23Aq4B/AX8mjQb594o1XSoEM5jUqRXYfTu45LoeMthGpCISpA+UR0ueC9AyFRYuhJYWM5ia3Ejg8ZIMRuqU9GTJJlPJkcFegP9QhoBqjpG0XoxxOXAVqd4/Kes2DemEQEfHi/lklx9nz08knSOQGl4rqRj/C+DVg9jPJNLp4YkVaJNUZJOAZ6fAkqeyLhFbke4sgcUhfXgtIp0NmUT6N9GxlPUM6VNoKqmIBHAYFpHU1MxgGupaW2GfifCLMfDqZ/r7YtZ3uZg00QymJpf9cE8CniV9CCzutslicmSwjgdelN7UqlZICiFsnJ0FI4QwmjTO6T2kMHNkttnxpPPBABdmj8mev7K3a/OlRrIL8P9IHSpuLn0ikC6tySGQRkYNpE4VUjMaBwwHFk6CuR1/xnb80JcJJIuypzuWiQDPkbpVP0dKQ8+SrtP/bzVaLRWLGUzKtMIuL4P/txy2/xXc/MqB7SYEWLbMDKYmNQ6YBuOezDIYMDfnS3vMYA9lD2ZUtJUqmKoNth1C2Ak4k1TLbwHOjTF+MYSwBWnq2SnA7cBbYoyrs6lqzwZ2JXWGOybG2OuQXHarViNoAUYBOwP/7P7kDOA60rSa/TQWZypXcxk9Jv3ynz4cli8f+H4mPQ9PzyYNgLE3KRGpoTnYdv+YwSSgBVpeDqMugZ1H95DB+iuy/hS8GUxNYyKwPE32s5R0ydryQexuEvB0x4NZpOvf1NDKZbCazNpWLYYYNYKXANcGYEQPT64mXY/c/X7H4z60UL3JrKSa2ggmLx7kJfWr0s306bD8acfPbiYWkorHDKbCezm85Atw7WeAf1Rgf5H01zbAajOYmsQ04Mk0dlglhjXqKEStBgtJTaJcBqvJGEnSUBRI3UOHBVKviFXdlntJ/UZXkebOmVby3LJsJ8N7P0ZPtSmp4UysQBEJ0oCpo2Hh02k4pI0H3zJJUgMKAYb/E4Z9msoUkSAFu46cNsIMpiYwhYoWkSB1AjeDDQ0WkqQqOQBYA1y9D+nyte7mAA9n96fSdUoESNMgrOj9GFmWkRraxGEVCjDDWD8X6WzKTEkrSWp6B7wK1lwKV7dW6QCrYdVwM5gaXISJ6yo/wdr6DFatf38qBAtJUpVEoK2FdAarrZ8vHk1nEamPf6Wr8fe0GlQLMBVauk8JMlBPk+agnZjqstMqtFtJUmOJl0PbV6nurJxrYHWLGUwN7GloqUJYWp/BHiObD1TNyEKSVCWXvwYOaANOIQ2UNBDDyT1WktRQZgBtMOWpyp8J6+CgOpI0NF3+ajjgitocywymhjQNprTB0iqFMDNY8/N3n1QNh8Cr/x9cCamIdEN1D7eGPodTkoachcDMejdCklRzr45wZaQmo2GbwdRwZlD1GW3XZzArSk3LQpJUaUcA7wYOrHdDpILaFHgiDcRYqavaypmPvaolaci5nPRXzsuqfJw2vLZNjWcBbLxJjTLYo8AWVT6Q6mJYvRsgNZ3zSaen/lqh/Q0jnVGLWPpV84jwVEjjzC+pd1skSc3l1cCl2BtCKuMpapTBatArUPXhn6VSJR0N/LFK+w5Ae5X2LdXS48CsKvZIWkaa0laSNDRdBuxX70ZIBTUDNn6y+j2SAJgLPFSLA6nWLCRJlfR74NfARVXafy/FpNU4Da0awFzSz/D89LAqJ4u77dRe1ZKkqmgF2sxgakChNh32Hn0EtjCENSULSVI11KErtb23VXhbAw+Sflhr+AM7F0+GSdKQczXVHyMJYA0ER9tWI1kAbFKbQ80NZrBmZSFJqqQ3Auelu5cB+9fw0KNIQzNJhXU/qZhUC08B09JdL8+XpKFlfQarxQfAKBi11gymBjKdVEyqATNY87KQJFXSH0h9mw9LD6sypFEA1m64+jnsVq0GUKtEUfLp9hCweY0OK0kqhpoNK/k8PLcWRhjC1CgWArNqc6iHNofN7ZLUlCwkSZV0JPA81Rtwu0MP8y1OIM28MKrKh5YGbCvg3hodazqwKN3dFni4RoeVJA0hq4BJMGEMLFljBlMDaavBMebAtlfAw57Na0oWkqRK+iPwVtYXev4OvLJGh14GzCRlGqmQAlCrcSQeJU0LR48d+CRJTa4mGWwisAqWrTWDqcE8Bcyo8jHWwdrt8Gxek7KQJFVSO13+ao3ANcDLK3mMCIzZcPXGwLOVPI5UaQ8A29XoWCNxBHpJGsIicM318PKKhrBuFgOtZjA1oKnAraQe3NXyAOvHq1TzsZAkVdqFpEvcMu3AdcBLKnmMlcAzdPlD2bNgKrwtgFvq3QhJ0lDR3g7XXQcvqWgIK5Gd2DODqeEsJk2AUs0f3p2A+VXcv+rKQpJUaW1s8Eu5HbgB2LuSx9mU9QMXT8HLd9QAHqbLP4LFwNxqHWtzUrdtSdKQ1r4r3PBX2Ptu4MUV3vlkmNJmBlPjWQzMfb7KB3kep21rYhaSpBppB24G9qzQ/qY8DO1LYGpLGh9JKrx20vSCmQg8RiqEzq7kcTYHrgA2gu3xZJgkDWn/gvbN4Oa9Yc/n+tw6lylToH0KTF1uBlNjWp/BnoHZFQ1hmW1h+wVmsGZmIUmqoTbS5ci7D2IfMcDGi+C/L4Dp28KSms1vK1XAo8AunQ8jsA64vZLHWAbsASyBp6nhFNCSpOJpA5ZB2zNw6wOw+8bAbgPbVYyw8cbw3//C9HthiWPxqYFFYF2E2x8HNqvwzpfD09EM1swsJEk11kb6o3kGsOtAdhBh8U6w81Pw1OKKNk2qvlnAGcDOnaueA/aq5DHuBdrhhREWVnK/kqSG1tYGty+GGXfArnlCWBtpOjZIf3XPhMV3ws47w1M74WU7anjrM9g8YE6Fdro9vHCxGazZWUiSquEK4NjyT7cBrcAl2eOlwA697K6dNCTSpsDsANzsL2c1qCeA19Jl/KIIPEQa87Ei9gOehkV4JkyS1FUb0NoOl9wJ7Notgz0PbFWycQu03wCbbgqbzobZC4HpsHAhBjE1hS4Z7DHS8AADsQ2pKrUDcC8sajeDNbsQY+OW0kMIjdt4Nb9RwOHA73p+upVUGIL0i3Y+G44TMxy4P3u+tfsO7E6tRjaDVFQqEejsWT2MNGvsgEyCXZ+GOzHENIsYo7/xCsYMpka2PoO1QPvskgwWST0zNuuWwR4tefEc0mXaUhPpksHmlGSwQJospS9jganAPNg1msGaSbkMZiFJqpaXAOcATwJfA84b2G46emnc37EiAtuWrpAaUCtp9pwbym+yLXA38ALgvhy73Bm4bld4xZ3wL8+ENRULScVjBtNQsEEGk4aILhmstMv4vdnKu0lZbgdgLez8R7hub3jFCvgXZrBmYiFJqrWxpNL+SmAacOMg9xfp7Ht9zyD3JRVBK2lQ7H+W32Q70o/7dsB/6L0j3j3AW8bDv54zwDQbC0nFYwaTpObWYwYrXQnr/ya55wXwlvs8kdeMLCRJ9TSadFprK+CPA3h9O2mGkTOzx7tUpllS3bUCewLX9b3pTt0e30aqQ63rWLEX/PtmaDfBNB0LScVjBpOkoaPHDHYLrNuL9UHs31hEakYWkqQiGAUcTN+Xua0G9gWuJ40KuRfwI+A92fO3Vql9Uj20AvsA1/TvZbsBt+0JXAu8gtTrzwTTlCwkFY8ZTJKGrt2A23YjTUXtp0FTs5AkFcUo4BDgDyXrlgNvBC7LHreTemh8nPTL+SZSQcl5FtWsWkjdizq0kgpEfRkOvIj0b8RPhKZlIal4zGCSJDU/C0lSkYwCdi95vA64Gdi7ZF2k6x/SLwP+Uf2mSYXxsuw2AFeTeh39PXu8L6ng+gXg1ZRc36ZmZCGpeMxgkiQ1PwtJkqTGtT9wRXZLdh9gP+Aq7I3U5CwkFY8ZTJKk5mchSZLUHC4mXR7qJ8CQYSGpeMxgkiQ1v3IZzBFXJEmN5adYRJIkSZLqxB5JkiSp0OyRVDxmMEmSmp89kiRJkiRJkjQoFpIkSZIkSZKUi4UkSZIkSZIk5WIhSZIkSZIkSblYSJIkSZIkSVIuFpIkSZIkSZKUi4UkSZIkSZIk5WIhSZIkSZIkSblYSJIkSZIkSVIuFpIkSZIkSZKUi4UkSZIkSZIk5WIhSZIkSZIkSblYSJIkSZIkSVIuFpKkAdgOeGG9GyFJkjTEdMlgWwJHArvXrTmSNCQNq3cDpEa0IzAauAuYC+wEzANur2ObJEmSmtoLYMdtsgz2EMx9Gez0I5h3Ftz+SeBFwJPAjfVtpiQ1uxBjrHcbBiyE0LiNV8PZDpjTw/pDgPcBVwDfKlm/GLilBu2SpGYXYwz1boO6MoOpZnaA7WZnGew44E3Z+gvgkPnwvvfBFVfAt24GTgb+AosPNYNJUiWUy2AWkqQ+bAfMBN4LvL4fr7sFeDdwWzUaJUlDiIWk4jGDqRa2A2b+CN77npwZbBFwOtzyGTOYJFVCuQzmGElSHw4FzqJ/RSSAF9O1h5IkSZJy2A6YBoduA2fN7EcGuxY4FV68oxlMkqrJQpJUxtbAXsA1wA0D3MeEbB87VqpRkiRJzWwb2Pr7sNd74ZpvwQ2H9eO1U4C3Az+HCXvBXjuYwSSpGry0TerBlsDHSBOBbFRuo02A4cDjfe/vRlJBSZLUf17aVjxmMFXLln+Cjx2RZbBHgDHAtAHu7F9w4zthr+eBeyrUQEkaQry0TeqHbwIPANeXrpxJmp5tWracDHyUVGnaCZhdfn9jSD2cJEmSVF6XDPYr4OpB7GwXGHMLbP3rwbdLktTJHklSN7NJuWV/SMWjidkT/we8Ffhe9viDpMGTHgU+k90/vvx+7yHN8PZQ5ZssSU3NHknFYwZTxW0Gs8fCr74P++9f2V3fcw8ccgg8ZAiTpH5x1jYpp9NJM8uOmw78BMh7bX4fhSSAu0idlyRJ+VlIKh4zmCruMjj9gCyDVWH3d90FOxnCJKlfvLRNyunTwHUAp5G/iPQssBLYuPfNRgCzBtwySZKk5rU+g1XBiBEwyxAmSRVhIUkqMXUq/HYWvGYWaWCjvP5Cunbt56RKUZkRurcF/jbINkqSJDWbqcBvgddUaf/bbgF/c6wkSaoIL22TOkyBy/4EB7wCWAaMAkYPcF+XA8dl99cASzufegDYG1g84IZK0tDipW3FYwZTpV12GRxwQBUPcBc8sCvsPRkWG8IkKRcvbZP68kfgFdn9DwN/HcS+DgAWkAbiPrPrU1vRbTY4SZIkVdcw2Go3uP6SejdEkhrfsHo3QCqSZ0kdiEac2deWOd0IvLlC+5IkSWpS6zNYtQ6wHWlilO2qdQBJGjrskSSNBSYAw+ANpKvSKuZlgGe+JEmSelXxDNaTFqozJZwkDTEWkqRLSGMYvbRK+2+ly1hL7cDzVTqUJElSQ1oJrKvuIdq3guerNS2cJA0hFpKkVwM3V3H/ewCXdj68P1slSZKkzGEMbnzKHO6/D/bYubrHkKShwEKS9A9gryru/3pgX2B4ergt8O8qHk6SJEkbMoNJUmVYSJK6aSNdflYxAdgTuLaSO5UkSWouFc9gkqSqsJAkdXMYcD4VDDIvAU4jFZMkSZLUo4pnsFKRVKmSJA2ahSSpB0cCf6jkDkO2SJIkqayKZ7AO/wZ2rMaOJWnosZAkxWzpQXv5p/pnd9YP6N3L4SRJkoaukpBUsQzW8+4lSYNgIUnakzQgdjfHAK3ArypxjJuAF6e79wLbVGKfkiRJzeRQ4IIKZzBI1aMd4d57zGCSVAkWkqRenAW8ncqdvfIsmCRJUt8qmsFuh9gCbFeJnUmSLCRJvTiO9I/kJxXa378xw0iSJPXFDCZJxWUhSaqR24Cd6t0ISZKkIns9cG5ld2kGk6TKspAk5fBe4LuD3MeuwJ0VaIskSVJTO5r1M96+9yeDzGA3wa7vhDsNYZJUMcPq3QCpUXwYWAv8bz9fd9118LKXVb49kiRJTS9WIINFYOdKN0yShi57JEkALyNND3J5+U1OBT7Sj11eBuwPxH0gXueUs5IkSf32Pjj1e/3MYJfB/vtDjNlCupUkVYaFJAmgnc7KTxkfBUYCw4HP5dhlBNqvBkaQClWSJEnqnwgfjf3MYK+C9k8CnwSuq2rrJGlIspAkdWgl/Ys4DPjzhk+3A+uy5cuk+lDHclLJdhdk6w4B/h7hleuAtuo1W5IkqZm1fxTWfQfWfR6+PAJGfLKPDNYCf2+BV14LvKL27ZWkZhdiA/fzDCE0buNVPCOAS+i8zK01/0tb6RxwrJ10HT8AVwH7VaqBkjQ0xRhDvdugrsxgqrlhpO7ebaTg9T/AN8tksL8BB1ESyCRJA1Eug9kjSeqwBjgYmAxMAP6Y/6VtwOpsWZ9Z/gZ8jV7HXZIkSVIOpT2820iB63RoGwur3wWrz4O1Y4GxpG7hFpEkqWosJKlfvgo80235Rl1bVGGrgTNJvYjWDXJfbcCVwOsG2yhJkjTUNX0G648PA5sBdwILgX2At5LGulwInFO3lknSkOClbcrtG8AHgFHd1q8GVgK/Ad5f60ZVw1jSFzUC+CVwdD9e+0fSBfrHAkcCKyrdOEkaery0rXjMYLU1ZDJYb74APAt8C/g48EXSSbvzgKuBl5IGTDoCeAdwYF1aKUlNxUvbNCA/BRaTTvKMZsMAA2kWjUnACcBPatayKlpB6o30PP3vFn0I6Zr8nwNnVLZZkiRp6BiSGaw33b8JI0k9j95Lym3HAguAn9W+aZI01FhIUq8+Rhrq5wekkzu9GQWMq3qLaux9wG+z+58Gflhmu45TgSPhoJFw9kXAu2rQPkmS1JSGfAYr9SVSVe07Jet+SrrEbQXwB2AzOOh/4ewJtW+eJA01w/reREPZM6Qu03k/k48gfa6/u2otqrFnSMWkjwFPk6YGWUsKLsA7gS2ATx4ObzgoXZr/EeCGNaTBuyVJkgZgyGewDl8jnaz7Aqm3+IeAT8E7W2CLpfDJu+AN34T9L4SPfARuaAfeXM8GS1Lzs5CkXn2PNNRPXmNIk541leVp+dzn4B0LSNfnA/wYxi+CYZ+DYy+CSz8O1wP3roQr6PvsoSRJUjk9ZrBPkC7l+iEbjLTdlBkM4BvwuWnwjk+SBoqakJbxwLD3w7Fr4NKN4Prj4N6D4YpR8I7b69tkSWp2DratXv0WeFM/X/Mc6UqvkyrfnLr5IvDBjWDiOtIFoZOBNwJTSIMSvBae/R94/mLY5EOpl/XC7LW3Al8B/pQ9fgh4dU1bL0mNzcG2i8cMVn09ZrCvkcYK+gLwemAWKaRkmjKDfRE++EGYOLH8Ns+SOittQpbBngcWwK23wleOhj9tD3wLHjrIDCZJ/VEug9kjSRU3Dphe70ZUyNeBG4GpwMQlwP+SRrSEVET6Laky9GsY/zcY/1LgezD2g7BlttkMYA9gbva4v+N3S5Ik8QnS+IsReA0ph/y+6yZNk8G+A19/Ldz4KZg6tfciEqTeSeOz+2OBLccAW8KMGbDH72HuKcB+sPYe4DHS90+SNGAOtq2qeCXwo3o3ogJmAt/9ArzxTtJsIGcDfwdeAPwFeAK4E/gHaaDHPUjT0pYYQ2cRCWBz4GJSgelv1Wy8JElqHtOAjYBzgaNIAzN+bcPNmiKDnQozj4Lv7gtvfOPAdzNmDMw9hPQ9eww2Pwou/gbM+KsZTJIGwx5JKusbDPyEzQS6Fk8a0SnAQcCUXwAXAI+TBm88jDS17OdJ3Yv+BhwOfBb4D+lMVy9GkrLfZTTpWAaSJGlQymaw75GKRx3Xzx8NbEW6hj7T8Bns23DKvlkGm0Mqng3GGNL3aCWMPBP2fwou+4gZTJIGwx5JKut04OZ6N6JOvk6akW3K50lnsTYHniKNRTCdFOCeIFWEfkxKbZ/PXrwv8O3e9z8S2BHYGLikkg2XJEkNr2wGOwI4JLv/ZuBFpEGRmsjXN4d37gZTdmPwRaRSo4HdYOQrYMevmMEkaTDskaSyHgCWlq74H1Lx5Pc9bt5UtibLLlsCVwJ7k77+udkGbyP1HZ8JbAHcBdxGuuTtTuDJfMcZAexaoTZLkqTmsEEGA/g+6fK2z5BOBf+FNF3sI7VsWZV9A7Z+eWXrRxsYBbwGRvwDdl1K6mkuSeoXC0nq1ZdJYzm+BlLBZAj4KvDS0hVHkEZwnFWybna2dDgIeAnpm3Vetu5m4HfAO4APlj/eFFKt6lnMMpIkKemSwSAVjD4HbAY8Sjq512S+eja8dG+6BbEqGAVsC1MON4NJ0kBYSFKv7iZNVNb2ITjodcBppMLIGuDX5V/3T+DT1W9eVfzhi/DyS2HjfYH9SL2O+jKzZLtNsts1pIGWbgMuBR4E3rfhS0cAu5HqVZIkSVCSwUjnq4A0U+xDvb+uoTPYUfDyrdNlZ1W3Bkb80wwmSQNR9TGSQgitIYTbQwh/yR5vHkK4MYTwQAjh9yGEEdn6kdnjB7Ln51a7bcrn38DJ28NfNgNOAP6Pzrnty3iKVD9pRLfvBoumATuRr4hUzm7AWaTBluYCNwE/6HnT0cCJgziUJEmlzF/N4d/AyaSr2PJq+Ay2Sd/bVcQU4CwzmCQNRC0G2/4QcE/J468Dp8YYtwKWkUoTZLfLsvWnZtupIO4CPv0zeP1n4E8fp+w4SdcAryd1xGlInyX1HtqXdKnaYGxMGhDzJdn9lwK3kmZc6WYEA58hT5KkHpi/msRdpB5Grwf+1Mt2DZ/Bam00cBSM+JUZTJL6q6qXtoUQNgUOJk1K+pEQQiBdLPTmbJMzSXNd/Zh0afLns/XnAT8IIYQYY6xmG5XfHXfAHeenK7Ta3wBH7gx/+X3XK9weBW6oU/sq4ZMvhReOhN23ADat4I6nAG8A2oELKrhfSZK6MX81nzuy5UFSlDiS1EupqTIY8EJg91oedBSds+BJknKr9hhJpwEfJw1VDGkShuUxxnXZ48fpHMJ4FjAPIMa4LoTwdLb94tIdhhBOxB6o9fFW4AG4cwJ8aXu48EKYAwyn1+GSGsf/wU07wFEzK1tDWm8ycBQptOwHfKzr02OAnwAnVePYkqSh5DQqnL/ADFYEdwJfAi6kyTIYaQSAo6hSBivnaeB9ZjBJ6q+qXdoWQjgEWBRjvLWS+40xnh5jfHGM8cWV3K/6cDawAtgeeBjuvBDOvj11rnmwrg2roFfD32bC/GoeYzJwHHD4hk+NBN5UzWNLkppetfIXmMGK4k5SLLuAJspgwN/+B+a/sMYHHQ7sbQaTpP6qZo+klwCvCyEcROqDMQH4LjAphDAsOyu2KZ1/t88nTaj+eAhhGDARWFLF9qk/riV1fH+G1Lc68+/6tKaxLaTzIgJJkirL/DVENF0Gey2wWY2POQZ4N7CODXqKS5LKq1qPpBjj/8UYN40xzgWOAa6MMR4LXEW6tBvgeODP2f0Ls8dkz1/p9fkF8w+6FJGaygeBrVKX5u2rfawRpFncejAKB8mUJA2c+UuN6qQfw/Z31/igz5BmI37QDCZJ/VHtMZJ68gngnBDCl4HbgV9k638BnB1CeABYSgo/Um0sAdakH7zV1T7WFOCjQBsbJJZ24MlqH1+SNBSZv1RoS8+H1S2kXts71uigz5P662EGk6T+CI180imE0LiNV7G8kzS37pyUoD9DDXomLSfNp/OtzlXPkK4pkCR1ijGGerdBXZnBVA3HAJ/5E2x/RI0OuBCYke6awSRpQ+UyWNUubZMkSZKkfvkDUOtL3CRJ/WIhSQJYSbrUjHSN/EXAf6t9zKfp0htpFWm+ZkmSpKFoFHDR7+C/d9X2uGYwSeofC0kSwHuBLdLdM0iT1K2s9jHHAx8C3pWO9R3gc9U+piRJUkGdQZbBLgHuqc0xzWCS1H/1GGxbKryZwPWkWWG3rdZBhsNz28BZz8BzwKeqdRxJkqQGMRO4/iwYE2DbTwDbVeEgz8JzZ5vBJGmgHGxbAjgZ2BJ4LTCrc/UJwMeobDFpOXAZ8MZHYf5c2LSC+5akZuRg28VjBlO1nQB87PgKFJMWARemu8vJMtgSmH+yGUyS+lIug9kjSQI4Jbu9jC6FpI65kV8O7E7/c8zCbJelFpC6UK8ClvVzf5IkSUPBLwDOzDLYflkG2wXYCfg36Zq03Ute8DDwD1J1aL+SDPYw8K60iRlMkirDHklSqcuAA3p+6njgQOCFwA6kS/dXkoY6WgZsDNzY7TV3A18qd6xHgbmDa64kDQX2SCoeM5hqaX0GOwp2+Azccwms/DuM/yIsmwwb3ww3zgfOB54HPl6SwR4APlO3pktSQyuXwSwkSaU+C7wHmF5+k2OAo4ALSAWkGcBDwPbA9/tzLAtJkpSLhaTiMYOpHo4BjjoSLhgNy86GGa+Bh7aA7X8M338JsA/wzTo3UpKaiIUkKa+/AAfX4DgWkiQpFwtJxWMGUyHMASYCd9a7IZLUnMplsJZaN0QqvNuAJ6u4/xXAX4FrqngMSZKkZvdCYL96N0KShh57JEk9+SLwbmDaIPezjDTI41bAvaSRIv8IvG2Q+5WkIcQeScVjBpMkqfk5a5vUH58FIl3Pcm0BrCb9q5ndbfuFwDPANt3W/wv4EPAB0jX7n8YikiRJkiSpYdkjScrrNCAAy4F3kKaXXUga6+h60qwgb+32mluB99WshZLUlOyRVDxmMEm92o00dtW6ejdE0mA4RpI0WI8DLwXmAV8gjaX0LWAv4O/AJtn90sUiUiHsCAyvdyMkSZKGgp1IJ2AnmMGkZmWPJGmgppIuc1tY74aoJ9sB9wFbAz8FjgYmAI8Aa+rXLEkDYI+k4jGDSVpvB+Bu0rAQ28F2N8B9B8HWT2cZ7F6YsDk8MgLW3E8aKkJSQ7BHkpratvU46GE43lFBbFNyu2223AC8GLiWNDTVROBvwNw6tE+SpKbQQ+DaFjYcI1LNLwtc22wL3AbbvBC23Ra2/Q3cMBZefC1cexd8+i6YOB3+dgnMvQt4FTCyvk2XNHj2SFJTuI/UA6Wt3g1RzWxOGrIK0oR4OwJ30Hc2eRVwDbC2ek2TVGH2SCoeM9gQ9QBpJtoS9wHb3QttO+KH6xCxORDuB7aEe0O3DLYT6SzehJIXvAz4DbxqsyyDHUI6u2fPJKnwymUwC0mSGtIjwGZ0FpP641Wk4PskMJk0frpZRiouC0nFYwaThq5HZmcZ7GGgtX+vfdUCuG8jePINMPkOWL4IVhvCpMLy0jY1vRn1boBqZibpErWyf8VMp9cK09+Ax4D9gYtI46JLklR10+j3H96NwAw2dMycCXMfhfgYA/pZ/tub4bH7Yf+L4KLHYK+D8VI3qQFZSFLTuBlnhRgqHidNktejjYHrgPF97+cSYPdKNUqSpN5sBFxK6srRYTJp4o5SE2m4P6zNYEPH44/DJoPpI3oVsENJBvsj6dI3SQ3FQtIQNYnOkwgT2TDDNJoppEuyvTS/uU3JFkiT5bVQsqLj/iPAK0k/1FNIPZMmZ/fL/MabkD1tCJYkVcVk4O/A22DS09CafaBNvAqG7UfnB9wU4GfA64ERwFgKX1SaMqVbBpvSw0aT8a+OBjeF9F5DSQargPUZbDwGMamB+Ct9CJhA+lu643Y8cAvwItIPwJ+BnYFxNG5B6T7gQdLnz8Q6t0XVsyRbwkTSGz0RWEQ6yzuRNPXsaODRko03Bu7K7r+QHn/rXZg9fXi2G3OMJKkSJgBhAky4B8IbYfyFcMvD8KIHoWUR/Hln2PmvMG4JDOv43PoTcChwNHBadltg990HD7ZmGWwiMI8UKktdR7ouXQ1rCbBkMYQKj1h3IbBkHBx+Hkw82AwmNQoLSU1qHJ2f4Q8Dc0jFlm1If1NvCdwI7AEcQpr16nzSIMTjSL2VRtG1sDSSOv9y7x5KOlaPg6mkk11rgftJtYT1xjKwEZlVXI8By0ijZLcCi7P7OwIru207js7fdP8iVU3L/OY7N9vNm2jcoqokqb7Gjcsiyzh4eBzMeRDu2wS2+Q/cNQe2fDncOBv2uKdMBvsdjLoYhh1Kmv3qrAJksDLGAVOnwuTRWQa7H0a3lTzZYXvSdelqTOPINWTAgJ0P574alo+GNw0zg0mNwFnbmtSK7HYL4B5SkaUvBwCfAvYFjgKOBS4A/pA9/x1S2Dm1gu3MrQV4mh4/xJ59Np0Ba2/vXDeaknrCY6TuV09VuY2quggwJnvwHP0vEK4kpfE9SP8wVgHtG262GngbcM6AWimp0py1rXjMYGWMgRWLgTmwxQNwz4RBZrATgd9kGWwNnLquOs0eqGdJPXnbl5N6AK8tyWDzSZ+380teMIr0IetPT+MIwLMQx5L+wBhDVU/Qrn4rvO13cE5b39tKqr5yGcxCUhMaSeqsMRt4iNStuhLWAicD3wWK/rt9CakX1nPZ45E4vXujGwWsHAU8z8ADzHakazm3yR6/FPgnGwTatwDnAeso/s+6NBRYSCoeM1jPRq6CZSMrnMHWwtoAJ38SvvsdaCvSB1MPAatLBnscRu4KqztO5t0NHEbqPq6GMWoUrHyeziphFX8jvwU47x2w7qyC/axLQ5SFpCFiOOnsUI/jMg5nUKNRf4jUw2kR8E1SR45G+f0+j3RSbEG9G6IBCaTOQ8Q09mhF7QNcv+HqdcA7gTMrfTxJ/WYhqXjMYBsaDjy7CkZWenDs98CHdoIt3gOLPgXf/GbqhV31P7Lz5MaVpCEEWku27fY6M1hjCwFWrQJGVCGDlbEOeOcJcOYva3RASWVZSBoiVlPml3wraTarvejsYtzKhpWgVlKFqI/vbBvwdeAz9HhlUP0MI336qPkMg2FrYQ0VPhG2D3ATG/zcHw/8lvSz7i8aqb4sJBWPGaybYbB6HYx4lqqPzdgGfP2r8JlPVTGDjQSWkr6WPNnqedI1bmtJgw1uSme38J7yphrOsFiFDFbG8W3w23dC2xlmMKneymUwB9tuIi2U+eXeQioe7Uk6HdSSLXfROfBwx3Ip8Jrsdb18UnyZdJn01yvR8MFo6XZ3BRt2x/KnvPEFaF2b8mnFA0wrcAWwf9fVZ5KO9+5KH0+S1HRanoMwklRM6T7pQ4V9GVjRAl+v9AdiaV5aTWdPoxW9bNexagydvZAm0VlEgtQlaUaF2qi6aAXWttVu7pozj4e1Z5jBpCLzT+wmspoeZvQIpEGmdwJuIX3IryWdWdqONItVW8lyCvBB0rSz5wLH9HyszwFfqWjrB6jk1MjKAMO7nZoLAAuBaZ2PPa3dYELqVl21jmb/AF5B2R8Mf14kSX1ZPRqGP0PKUmP62npwPgd85RPANyq40ymkGVC7a2PDE3Q9XO62kg0z6PrPz5nAgvRZ7mdq4wlkGWw4tese1NJ5bEnFZCGpiUR6+P2+jBQOniSdDZpDOq3wQJmd/A34HmlwmHPoc9qquv+CH8b6L3rkalg7gS6DPi4FpkwjDexEurqvjc6xllVwHd3h2yDGKueXy4CDNlz9Q1JtVZKkcuIaiLUaQAbgaxA+VsH9LSXlxTxaN1w1kg3rS913+cgj0LaFGawRRSC2UbvgfxbwDjOYVGQWkprEOtJECmtXkz7NV5B+60/stuG8bP3Wvezs/wGvy3fcj5LqTnXVBrTAmhEwvCPFPAeM33DTCGwF3FezxmlQ2oBhENuhpSX9wqpqMeli4PCuq95LAX7GJUmFtW4djB42qPlMBqQQGQxgLawZ1kOv+G5iLMlgHcMrqPAi2QgYLek9rCUzmFRcFpKaxDBg1VoYMYI0vVUlulWfR+cAxO+qwP6qJRsgfAQlIW4caVyorJC2EJj5OMydnabjVYNoJVVJswpSOzU4GXY+8MbOhz/Gs2GSpDLaYNgIWLWudjNaAfBJUj4rwl/Zw2HE87A20uO0wQtJV7fNnQsPGcIaVk0yWDdmMKm4LCSpd8eSRh3+KfC+njd5P/CT2rWoU19TlUyAyU+nrtVAmkFkXnWbpArLeiR1qHqPJEmS+msdXT6raubb8P774Cd5Q9hOwH/6eYyOEzp9GUEat7LD0zC5HZaWu1zuhcAd/WyL6qqlHWItg9jPgZNqdCxJ/WYhqUlULcP8BjiOdAri+6TizSe6bvIDKvx7fhxdZ/sop+Ond13J/VV0ng17mvU9kqYDT1SsgaqXmpwNO4I00LwkSX1phXXt9akjEfuZwe4Edui2bgZpHM1So+mcea7bCZ1ejSSNzTmGlL+eTqvNYM2hvaUkg7XQ9wndwQrAj7BLklRQFpKawToY3gLrup8NqoTSac7KTHn2flLX04rpGN9oJKkw1JfhdH6YjaZzsO1JwMPAFFiAM882rHbWD+zZir2SJEnFMnx4Giep5j4K7/9eBTLYxqQZfjusBCbTcwZr62NfY4Hns/tTgKUlGewRYPPBNFT1tD6DlQ47UE0nZActwuWbkjZgIakZjIC4CkauhbV9jXRYCV8FPtv5MADvBn4x0P1NJ6WMUpFUECod62n9/KM9bFvu/lRg6fqZZ9WIWlhfIO2YmbBksr7K+yNwVNdVp5IGNZUkqbsYYeRIWFvr0bYDhDDIDLaA9KE6rdv6cuNt9jmi9ob312ewMick1RjWZ7AIsfQkbpUPagaTislCUjNYBQyD9laItfiA/jTwpa6rBpUNFpLGL+pJ6YdUmUEce5W9Ps7HLkmNqh0YteGqqmWYFrr8ML+PNPTDd6pxLElSU2hvh1jaK7qGwrsg/HwQO+iYWAXS8AJPk/JWT8MMDOCD117EzaUdGN4O7SPou4faYJwO7zvRDCYVlYWkRreadCYp+8N3LDXIMKWBI/MzBjmxW94Pou7brSZfBWsOnV2S/gtsnfN4KoYefj6qmV1KtWfHMghLknoztq0OdaQfws9Gwbve3c/XzSU19lHSpW2rSZe3PQdslD2e0M99PkfvJ/y2xKlzm0Tb+v9VwZtJY7S+B9p/bgaTispCUiNbSfqQL/ntWpNL9L9M6gX1xc5VHX9sD9TUbJfdx3vs03hgBet7ZQGwnA27Y5d+Y3YCHuh3E1UvLaT3eNSGT40m/dyNI/+QWrn8Gt7+prTP0yu1T0lSU6vHMEm0Q/taaMsTwl5I+qC8i3QSbgQwDKY+BavGw5PbZ88/l90u7WEfK3tYB/AsPV8KV6ou3yBVy2igbQyMW1vhDHYmvP0YGNkGp9fi8jlJA2IhqZGNJH3Id+uRM4n0WZ5n4rMBGQZ8ixQmvppWnQD0OvvsHsC1PT81HbiVNK7jHFKvqsWQvr7lfbRlDenMWemMIiOzHYzr5TWe2mgc7XSeHY2kH+7s/VtDZy1xEemtn0rKs4Px1uGw7xmw7Dh4xyD3JUkaOiZNgjFj4LmqhbAN9ZnBOvybFLaOAS4BxsD0LbMMtgbmrIGxI2FxGymc9dS7aEoP654lhbdJ1OXSPtXHGmD8GlgxARaNgZFjYOrzFcpg74JlZ5rBpCKzkNTIxpNOB3QrJK0i1XimA89U69gdA0pmgy6eBXyot+1vAQ7otm4OcE9q/mzgXmDVBFixID01cjUsnZWjLStJ08w+lbVnKqmbdg1DnKqs4wzoeFLFqIenZmVPLyH97C8HNttw81x+dhxcNQWm/RbOGMDrJUlD06pVsHJllTNYN31msA6R9KF5N/CidD+syjLYZFh1N6yYnmWw5bB0E2AZ6a+Fji+m9EP1GVKIm076AK7VF6zCWAmwCmathPErYck0mP4sLJ8Fmy3JmcEOBK7sfPiz4+Cqs2FamxlMKjILSY1sBensT5lunyuoYsebk1k/c9tPSQMS93oSqp0Nu0O3ANsA96YcMgbgGRi7GzyQjWc0dkXJtst62f8KUvJZm91fUfLcQ6SqwvTeGqiGsIJUMSrzVMfb/jypiPQ4nf8GdgAe7mP3R5FqkpeshnUrYMU6e+JLkvpvxUyIE+k9uwzGt4FP5cxg3UXSB+UsOjNYAGbA2Ps6r/4fu4L0obiYVG1aTtfP4I6e391zF6QP4P6OsaSGtf5HYAU8Pws2ewIe3zz7N7CgWwbbk/RzdRvwauBvwKFw1ESYOBEu+R2sW5f2ZwaTCizG2LALnTNRDu1lMpF1xJ7+mxSJkyPxqSq9Dd//JnH0ANu9WQvxzsnEGROJiyBuDnEJxKUQ2ycRF0Cc1P3rLHm8JLtdDLGlt2NNIrKIyMwCvFculVkmE2mnx5/57v9NzP4NtJTc71huzX6OD8keD89eM/ZZ4oijCvB1uri4RCDWO2+4bLjU+2eiEZZJECcvqUIGO4UYRxO/P5I4evTA2rYZxDsDcca2xEX3EjefQlwSiUufJbbPIC5YQpy0kPR525G/JhNZku4vWUJkSi8Z7BkiE+v/HrjUd5kIcfJEYss84sRdiZMnEye3ECdDvHU8MbYSDyE9Hp69ZizEEQVou4uLS1rK5oB6BxFDTIWWqUTWEsv9NyUSF1bhbfjmCiJfHVibWyBOyW63gfjUJGJ8gBjHEp9/lDgT4vySbZ/s9vr27HZqdvsEnR9CGywLsZDUbEvOQlJP/10TiTvEVGidGjsLSLyJyMVETiAysgBfo4uLS4TyIcalfku9fyYaZtmIOGUqceFUYpxPjLsR44OD+NZ/kRjHECPEbw6iXS0Qp8wktiwmbvMk8alAjFOJcTPi848RZ84izm8nsojYsrAkg01N69qnpsdlM9hGREIBvv8uxVimEGlJ96+BuAOp0DqVXrK7i4tLIZZyH0ghCwMNKYTQuI2vhrV0HXS6myeAg4E/k3ooD9RXSOMvjvgefPpr8OwK+jey3tbA34HHgNcDt8DGS+DenWDPjeG2p2CrabD8yXS5fcdYj9OBhSW76Xg8H5j9OGw8B54sN2vJNFLXbGd/aB7TSW/+AC7QnUrqob9Bl+mlpGssjwfOHUzjJFVSjDH0vZVqyQzWP08AB28Mf14Ks+8hZaH++hx85dsw8iQYMSLLYINpVAswDTZuh3vXwJ43wm0vgq2eyTLYRjByc+BRmD6sJIO1w/SFsHBmlsFIQ1P2e+ZdDUllM5ikQiqbwep9RsuzYRVceumRRCTOiMTWSJweibMi8b6c3+oPR+Kvs/ufj8Tx2TLhK/1o225E5hH5K5HtiDxIZDqRViKziC07EmfMI4Z5xK3nEe+7L71uLsSHIZ3FeLSH/T5GnDWP+MhmxNZ6f/9dar/MKllW0+vPf+7/TiAywEsFXFxcqrPUO2+4bLjU+2ei0ZYZpJwyHeKsTVLOia8ixjuJ8VBinEWM1/aQwT5M/PWs9PznxxHHQxw/njhhQuXa1gJxRiCGmcStx6V8yCri3FlZBptNynCPZa+ZReoVDHEWxEcwg7m4uLg061I2B9Q7iBhiKrjMyZZVxL7+uyoSXx6JcyLxrkg8PLs/JxKvzr7FJ2WPfxyJz0TiZyJxQule+lNIGpm1bT8iV5FCSMdzG5OKRMvSNuO2J65Y0fn8HFLIeXBOD/ttJ7J52qbu33+X+i4dP/vbEXmWPv8NEIm8jsgtJY/fTmRMAb4WFxeXLku984bLhku9fyYaebkK4stnEucMJ941g3j4iJRj5mxCvHoOMc4hnjSHOGcO8cdjic9A/AzECTVo27hAXDGH9Jk6jDjnYWLLw8QHI+uLR5A9D5EH0za01v/76uLi4uJS+aVcDvDStma0is7rwcqYS+qOvJbUJfmp7GUAM4GxwALgC8D9wBWkKTyfLt3JV4FP9bNtI4AZwKMl62YAd5Fmbjgdwu6wxQ7w4O+AXYE7oWVHWHs3tL4A7gW2vRfYltQ1/P5+tkHNayvgQWALUpf9G4HJJc/vD8wDLgb+F/grqT/+6Oz5+aSZbCQVSvTStsIxgw3cXCqQwaokkD5CHxwB3A4cnj5O1wKtwL33p/i13tbAPcAovFZJkppQuQxmIakZvYCUBG4hjfcyCLOA5+gWXr5DKlQ9Tf8LSR3mAKcDryEVkh4DbiWNTXMv6Q/754B/A8cCv4Pt3gj33APb/QfuOZIUXHpyG7APnalMQ9c2wD+BN5BS+AOkNLwFKcX3a75kSfViIal4zGDV12MGq5VAqhj9Nz3cjhS7Om676HGlJKkZWEgaip4FxlV4n98CvkE6PRXoOvp1f4wAXgZ8DHg7aRTKJaTTdNcCa4B3kjrU3QXsAvyK1EOpnd4HWH4hqQDlT4cAtgfuwzOlUgOzkFQ8ZjBJkpqfhaShaFc6Cy6XA8cAZwGbDHB/3wC+TerZMVizSZcXTQK2BHYi9eX+MqmX09WknkXrgL2z17wY+DFwEqn3kiRpSLCQVDxmsCEmkC4X36PeDZEk1ZKFpKHuxcAdwEPApgN4/SnAacDRwDLg7EG2ZxSpl9Ea0nhJZwFvJc0dexupsPQaUgGsY7ynALQBN5AuXctcR+rc1D7IJkmSislCUvGYwZrMBOA80niV5exFymAlzGCS1NwsJCnZEzgf+CjwNeDzpMGJf0K6BKjUZ4C/Z/fvJ13GtgVpjJl5FWzTCGA30uVHjwIHAcNJPZO+RBptsh14JSmtAPyj8+Uv6/pQktRkLCQVjxmsyQwj9Tb6Z/9eZgaTpOZmIUmd9iHNxLEzafyhFaQeSxO7bXcnaSqRWhkGvAS4iXTZ2ytIl7i9gnSG7Mpu21+e3b4ax0OSpCZmIal4zGCSJDU/C0lqHOOAxcDh2eNImqa9uwOz20tr0CZJUt1YSCoeM5gkSc3PQpIaRyupSHRxvRsiSSoCC0nFYwaTJKn5WUiSJEkNyUJS8ZjBJElqfuUyWEtPK6VCaSHN6iZJkiRJkurKQpLqbyLw3V6ej8AVNWqLJEmSJEkqy0KS6m8tcEsvz0fgzBq1RZIkSZIkleUYSZIkqdAcI6l4zGCSJDU/x0iSJEmSJEnSoFhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpTLsHo3YCgaB+wK/KPeDZEkSWpWU4EtgPuAvdKqLhlsFXB1XVomSVJDs5BUZVOAucCDwIuzdXOBzwLvyB6vBm4iZZxrats8SZKkpjRlDsx9Nzx4J7z4tLRuLiUZ7ElYfSTcdBPs9bIsg60D/l6X5kqS1DBCjLHebRiwEELhG/9q4BTg58APy2zzFPCubJsZwB7AP2vSOkmSii/GGOrdBnVV+Aw2BV49E06ZCj+/qpcM9hS8613w8wuyDLYC/vlaoB24rmatlSSpkMplMAtJVbARsHV2f0/gtJyvWw68ETgf2I/US0mSpKHOQlLxFD6D7QN7vgtOu5nyVaQSyynJYG1w0z+AV1atmZIkNQQLSTV0APAjYKsBvr4NuAH4EBCB2yrULkmSGpGFpOIpbAabDD/aIstgewE/6N/r24AbVsKHXpZlMIOYJGkIs5BUQxNJ195/Z4Cvfw54FamYtBoYVaF2SZLUiCwkFU9hM9jR8I5zKpjB1gIjKtU6SZIaS7kM1lLrhjS7icA7GXiAgTSjyA3Z/QBsN9hGSZIkNbmqZDCDmCRJG3DWtgp7HfAtgLHABNJF97N62HAN8Fjf+xsB3EnKMJE0+5skSZK6et14+Nb0yu1vBHBnK2x3PsSD4UFDmCRJgIWkihoHTO148FLgeOD/AWf1sPGDwCuAhcCm2bp2YN6Gmw4D7gMeBrasYHslSZIa3ngYNwWmvp7BdUfqwbB1cN9B8PDlsOV+wKOV3b8kSY3IMZIq6BjSxCBTRgNjSIM8/qWXFzySvaijD/VSYCdS16Mnum7aBszJ7s+vUHslSWoEjpFUPIXKYMfBMV+BH46BKVMqv/u2NpizOXAdzN+s8vuXJKmoHGy7Rt4KnHUYcBRwKXB2P3fQTuqttE3PT0cc2EqSNLRYSCqewmWw0fDWE+Cs71dp/+0QZ2UZbGGVjiFJUsGUy2Be2lZJo0jXt10HTKL/RSSAZ4B9SNfIRWDJhptMBRYPsImSJElN581AtYpIkCpIC2BqhMWe0ZMkDXEWkirpWODng9zHJOCp7P5K0iVyJQKwCNgIWDbIQ0mSJCmf9RlsMiwzhEmShjDPqRTd+J5XP1zbVkiSJBXTcKAVeL4GxwrwcA+9xSVJGkosJFXKMGAkrANWVWqfo0gDco/d8KlJlTqGJElSIzse2AnWnVTBDNaLSTU4hiRJRWYhabCGAaOB9wE/hHOBEyu171XAFqQZ3EZ3fWp0D5tLkiQNOT8H3l/hDNaL0ev/J0nS0GQhabA+SBqs6Jvp4dHALyq179HAk8BMuoyuHYAVwMhKHUeSJKlRtQDD4Og2+MXa6h4qACsCjHSMJEnSEBZiLNbsrf1RmKlnW4D3sn62kKOBcyq5/7V0Xt7WLSA5H7IkqdmVm3pW9VOYDAbwbmBX4CQ4+mg4p6IhrAcRaDGDSZKaX7kM5qxtlfAB4LQq7n84sAZoJw0mKUmSpOSn2W2gdn3tW4G2Gh1LkqSC8dK2Koikmk9V+I5JkiRt6FiIZ0N7G1UMYtm+V1dx/5IkFZxliUoIpOpR5lzgyGocp4U0LZwkSZK6+jWcOwyOHAYcVqVjRFJP8eFV2r8kSQ3AQtJgfQzYjDRrWzfVHjygOIMTSNLQ8HHg1Ho3QlLvjoB4UZX2HYB2iF7WJknV9w3g0+nu+gz2HuBndWuRMhaSBuubwEc2XH0+cGgVD5uN8yhJqqLTSL9vO5avAx/O7pthpGKqegaL0GIIk6Tq+DkQ4bQI8WMQv5R+7349wocjxB3hZ++qdyPV6LO2PQWsABbXuy3qYiq+J0Xje1JMvi/F43tSPHNijBvXuxHqKoTwLHBvvduhDfg7rHh8T4rJ96V4fE+Kp2wGa+hCEkAI4ZYY44vr3Q518j0pHt+TYvJ9KR7fEykf/60Uk+9L8fieFJPvS/H4njQWO+ZKkiRJkiQpFwtJkiRJkiRJyqUZCkmn17sB2oDvSfH4nhST70vx+J5I+fhvpZh8X4rH96SYfF+Kx/ekgTT8GEmSJEmSJEmqjWbokSRJkiRJkqQaaNhCUgjhwBDCvSGEB0IIJ9e7PUNJCOGXIYRFIYR/l6ybEkK4PIRwf3Y7OVsfQgjfy96nO0MIu9Wv5c0rhDA7hHBVCOHuEMJ/Qggfytb7vtRJCGFUCOGmEMId2XvyhWz95iGEG7Pv/e9DCCOy9SOzxw9kz8+t6xfQxEIIrSGE20MIf8ke+55I/WAGqw/zVzGZwYrHDFZcZrDm0ZCFpBBCK/BD4LXA9sCbQgjb17dVQ8oZwIHd1p0MXBFj3Bq4InsM6T3aOltOBH5cozYONeuAj8YYtwf2At6X/Zvwfamf1cB+McadgV2AA0MIewFfB06NMW4FLANOyLY/AViWrT81207V8SHgnpLHvidSTmawujoD81cRmcGKxwxWXGawJtGQhSRgD+CBGONDMcY1wDnAYXVu05ARY7wGWNpt9WHAmdn9M4HDS9afFZMbgEkhhBk1aegQEmNcEGO8Lbv/LOkX9Cx8X+om+94+lz0cni0R2A84L1vf/T3peK/OA/YPIYTatHboCCFsChwM/Dx7HPA9kfrDDFYn5q9iMoMVjxmsmMxgzaVRC0mzgHkljx/P1ql+NokxLsjuLwQ2ye77XtVY1vVzV+BGfF/qKuu++y9gEXA58CCwPMa4Ltuk9Pu+/j3Jnn8a2KimDR4aTgM+DrRnjzfC90TqDz8/isXP+QIxgxWHGayQTsMM1jQatZCkAotpKkCnA6yDEMI44I/Ah2OMz5Q+5/tSezHGthjjLsCmpLP4L6hvi4a2EMIhwKIY4631boskVZqf8/VlBisWM1ixmMGaT6MWkuYDs0seb5qtU/082dEtN7tdlK33vaqREMJwUoD5TYzxT9lq35cCiDEuB64C9iZ1YR+WPVX6fV//nmTPTwSW1LalTe8lwOtCCI+QLsfZD/guvidSf/j5USx+zheAGay4zGCFYQZrMo1aSLoZ2Dob5X0EcAxwYZ3bNNRdCByf3T8e+HPJ+uOyGSr2Ap4u6earCsmuGf4FcE+M8TslT/m+1EkIYeMQwqTs/mjgANK4CVcBR2abdX9POt6rI4ErszOYqpAY4//FGDeNMc4lfW5cGWM8Ft8TqT/MYMXi53ydmcGKxwxWPGaw5hMa9f0IIRxEus6yFfhljPEr9W3R0BFC+B2wLzAVeBL4HHABcC6wGfAo8MYY49Lsw/UHpFlGngfeHmO8pQ7NbmohhJcC/wDuovO640+SrtH3famDEMJOpEECW0lF+3NjjF8MIWxBOhMzBbgdeEuMcXUIYRRwNmlshaXAMTHGh+rT+uYXQtgX+N8Y4yG+J1L/mMHqw/xVTGaw4jGDFZsZrDk0bCFJkiRJkiRJtdWol7ZJkiRJkiSpxiwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSaqpEMKBIYR7QwgPhBBO7uH5l4cQbgshrAshHNntubYQwr+y5cKS9ZuHEG7M9vn7EMKIbP3I7PED2fNzq/4FSpIkFZAZTFKlWEiSVDMhhFbgh8Brge2BN4UQtu+22WPA24Df9rCLlTHGXbLldSXrvw6cGmPcClgGnJCtPwFYlq0/NdtOkiRpSDGDSaokC0mSamkP4IEY40MxxjXAOcBhpRvEGB+JMd4JtOfZYQghAPsB52WrzgQOz+4flj0me37/bHtJkqShxAwmqWIsJEmqpVnAvJLHjwOzQghfDCG8rsxrSo0KIdwSQrghhHB4tm4jYHmMcV3pPrsfL3v+6Wx7SZKkocQMJqlihtW7AZIUY/xszk3nxBjnhxC2AK4MIdxFCiaSJEnqJzOYpIGwR5KkWpoPzC55vGm2LpcY4/zs9iHgamBXYAkwKYTQURgv3ef642XPT8y2lyRJGkrMYJIqxkKSpFq6Gdg6m+FjBHAMcGEfrwEghDA5hDAyuz8VeAlwd4wxAlcBHbOLHA/8Obt/YfaY7Pkrs+0lSZKGEjOYpIoJ/nuWVEshhIOA04BW4Jcxxq+EEL4I3BJjvDCEsDtwPjAZWAUsjDHuEELYB/gpaQDIFuC0GOMvsn1uQRo0cgpwO/CWGOPqEMIo4GzSWbOlwDHZmTRJkqQhxQwmqVIsJEmSJEmSJCkXL22TJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5WEiSJEmSJElSLhaSJEmSJEmSlIuFJEmSJEmSJOViIUmSJEmSJEm5/H9jL0Vj8tYLzAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1440x4320 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "c='8'; fov = '201'\n",
    "ref = imread(f'tif/Cycle_0/Cycle_F{fov}.tif')\n",
    "mov = imread(f'tif/Cycle_{c}/Cycle_{c}_F{fov}.tif')\n",
    "\n",
    "ref_max = ref.max(0)   # max proj of z for each channel in reference\n",
    "ref_binary = ref_max[0] > threshold_otsu(ref_max[0])  # binarizing for channels - T/F\n",
    "\n",
    "mov_max = mov.max(0)  # max projection by each cycle in for loop (ref above)\n",
    "mov_binary = mov_max[0] > threshold_otsu(mov_max[0]) # binary of moved image\n",
    "\n",
    "out = imread(f'reg_bin_Cyc_{c}/Cycle_{c}_F{fov}_bin_reg.tif')\n",
    "\n",
    "fig, axs = plt.subplots(1,2,figsize = (20,60))\n",
    "axs = axs.ravel()\n",
    "\n",
    "im_reg = np.zeros((2048, 2048,3)) # empty color image\n",
    "im_reg[...,0] = ref_binary\n",
    "im_reg[...,1] = out.astype('bool')\n",
    "    \n",
    "im_orig = np.zeros((2048, 2048,3))\n",
    "im_orig[...,0] = ref_binary\n",
    "im_orig[...,1] = mov_binary\n",
    "    \n",
    "axs[0].imshow(im_orig[:500,:500])   #before reg\n",
    "    \n",
    "axs[1].imshow(im_reg[:500,:500])   #after reg\n",
    "\n",
    "axs[0].title.set_text('Before Registration')\n",
    "axs[1].title.set_text('After Registration')\n",
    "for ax in axs.flat:\n",
    "    ax.set(xlabel='0:1500', ylabel='0:1500')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "cac9fd61",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJIAAAI4CAYAAAAxusOtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAACAh0lEQVR4nOzdd5xcVd3H8c/JpgdC6CX03qR3RCkqIFUsoKKoSFERsIIoYkWwoYKioAIKUh4VKSrSQXqv0kInBEIgCSGk5zx/nDu7s5vZ3bu7d+p+3rwuO3PnzpmzbfLd3z3n3BBjRJIkSZIkSerNkHp3QJIkSZIkSc3BQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkqV0I4XMhhFdDCG+FEJaud3/6K+v/mnV67Z1CCE/U47UlSVLzCCHsGEJ4Ksst+9e7P31R77wTQjghhPD7er2+NNiFGGO9+yCpICGE54DlgQXAPOA24MgY44s5njsMeBPYLsb4YDX7mb3ejcB2wHxgNnAz8IUY46Rqv3b2+jsD58cYVx5gOxFYJ8Y4oYh+SZKk1pJlnk2BFWKMc8r2XwdcHmP8ZXa/8EyRtfk2EIHpwMXA12KMC4p6jV5e/zvA2jHGgwfQxs4UkNkkFccRSVLr2SfGuBiwIvAqcHrO5y0PjAQe7esLhqQ/7ydHZX1dG1gM+Gk/2qiaEMLQevdBkiQ1rxDC6sBOpELOvl0eXo1+5K5uXqenzLJplrfeDRwIfKaI1yzCADKkpDryl1ZqUTHG2cBfgQ1L+0III0IIPw0hvJBNYfttCGFUCGFdoDQ8eVoI4frs+B1CCHeHEKZnH3coa+vGEMIPQwi3ks50rRlCWD+EcE0I4Y0QwhMhhI/k7Os04B/AZmXtd9tWCGHpEMIVIYQ3s379IIRwS9njMYSwdnb7/SGE/4UQZoQQJoYQvhpCGAP8G1gpG07+VghhpRDCd0IIfw0hnB9CeBP4VAhhmxDC7SGEaSGESSGEM0IIw7O2b85e8sGsjQNDCDuHEF4q68sG2ddqWgjh0RDCvmWPnRtC+HUI4Z9Z/+4MIayV52smSZKawieBO4BzgUNKO0MITwNrAldkGeL27KH2TJEdt3cI4YEsR9wWQtikrI3nQgjHhRAeAmb2dgIsG+l0K53zVk/tbxFCuD/LKP8XQrg4hPCD7LGueee4LGfNyHLbbiGEPYATgAOzz+nB7NhKGfLTIYTHsuc/E0I4Iju2p8x2ftnr75vlrGlZ+xt0+Tp9NYTwUJZpLw4hjOz1OyepWxaSpBYVQhhNOut0R9nuU4B1SQFibWA88O0Y45PARtkx42KMu4YQlgL+CfwKWBr4OfDP0HntpE8AhwOLA68B1wB/AZYDDgJ+E0LYkF5kbR4ATMjuj+mlrV8DM4EVSKHskK5tlvkDcESMcXFgY+D6GONMYE/g5RjjYtn2cnb8fqQC3DjgAtI0wS8BywDbA7sBnweIMb4re86mWRsXd/m8hgFXAFdnn8cXgQtCCOuVHXYQ8F1gyezz/2FvXy9JktQ0PknKExcAu4cQlgeIMa4FvEA2kjzGuH12fHumCCFsDvwROIKUxX4HXB5CGFHW/keBvUj5bX5PHQkhrE8aHVXKW922n500u5RUAFsKuBD4QDftrgccBWyd5a3dgedijFcBJwMXZ5/TpmVPK8+QzwOTgb2BscCngdNCCFv0ktlKr79u1r9jgWWBf5EKdMPLDvsIsAewBrAJ8KmevlaSemYhSWo9/wghTCPNg38v8BNIQ4dJ/2B/Kcb4RoxxBukf94O6aWcv4KkY459jjPNjjBcCjwP7lB1zbozx0Sy47EEKDedkx98P/A34cA99/VUIYTowhVSo+WK2f+/u2gohtAEfBE6KMb4dY/wfcF4PrzEP2DCEMDbGODXGeF8PxwLcHmP8R4xxYYxxVozx3hjjHVk/niOFrHf30kbJdqQpe6fEGOfGGK8HriSFvpJLY4x3ZV/DCyg7SyhJkppXCOGdpOlrl8QY7wWeBj7WhyYOB34XY7wzxrggxngeMIeUL0p+FWN8McY4q4d27gshzAQeA24EfpOj/e2AoVn782KMfwfu6qb9BcAIUt4aFmN8Lsb4dC+fW3uGzNr/Z4zx6ZjcRDoJt1MvbZQcCPwzxnhNjHEeaamEUcAOZcf8Ksb4cozxDdJJvs1yti2pAgtJUuvZP8Y4jrTe0VHATSGEFUhnaEYD92bDfqcBV2X7K1mJdIao3POkUUwl5Yt4rwZsW2o7a//jpFFD3Tk6xrgE6czQkkBpEcWe2lqWFGzKX7unxcQ/CLwfeD6EcFMIYfsejl2krRDCuiGEK0MIr4Q03e1kUtErj5WAF2OMC8v2df0avlJ2+21S4UmSJDW/Q4CrY4xTsvt/oedR1F2tBnylSx5ahZQvSnq9oAqwBSlfHAhsC4zJ0f5KwMTY+cpMFV8rmzJ3LPAdYHII4aIQwkqVju2urRDCniGEO0Ja0mAaKbv1JW+1Z9Ysd72IeUuqGgtJUovKziz9nXSW6J2kUT+zgI1ijOOybYls8cVKXiYFjHKrAhPLX6bs9ovATWVtj8uGH38uR18fBn4A/DobOdVTW6+RrvRWfuWOVXpo++4Y436kqWX/AC6p0PdOT+ly/0zSSKx1YoxjSXP9Q2+fU+ZlYJXQeRHJrl9DSZLUYkIIo0jTqd6dnYx6hTRVftMQwqY9P7vdi8APu+Sh0dko8ZJcl+DORvpcAtwOfDtH+5OA8VkuK+kpb/0lxlgagRWBU3vpX/v+bKre30gjiZbPToj+i4681dvn2CmzZn1eBfOWVDUWkqQWFZL9SCN9HsvOzpxNmnO+XHbM+BDC7t008S9g3RDCx0IIQ0Na9HFD0tSsSq7Mjv9ECGFYtm1dvthhL84jXTlu357aiulytX8HvhNCGJ3N9/9kN1+D4SGEj4cQlsiGOr8JlEYHvQosHUJYopd+LZ49763stboWxl4lLZZZyZ2ks15fzz6HnUlTAy/q5TUlSVJz2590Mm9D0jSqzYANgP/STW5h0UxxNnBkCGHbLNeNCSHsFUJYfAD9OgU4LBut3lP7t2f9PyrLgfsB21RqMISwXghh16wgNJt04rI8b60eer4y23DS1LjXgPkhhD2B95U93ltmuwTYK6QFvocBXyFN0butty+GpP6xkCS1nitCCG+Rih8/BA6JMZYuLXscaYHFO7JpWtcC61VqJMb4Ommtoq8ArwNfB/YuG57d9fgZpH/0DyKdGXqFdDZqRKXjKzx/LvBL4MQcbR0FLJHt/zNpgcU53TT9CeC57PM9kjRFjhjj49nznsmGc3c3BPurpPUMZpAC18VdHv8OcF7WRqer1GWf0z6kRSKnkNYk+GT22pIkqXUdApwTY3whxvhKaQPOAD4eKl9h7TuUZYoY4z3AYdlzppIy3KcG0qlsFPjNwNd6aj/LMAcAhwLTgINJJ/oq5a0RpALVFFI2Ww74RvbY/2UfXw8hVFynMst9R5MKQlNJuevyssd7zGwxxiey/p2e9WEf0iLmc3N8SST1Q+g87VWSmk8I4VRghRhjX9YdkCRJUk4hhDuB38YYz6l3XyTVlyOSJDWdEML6IYRNsmHY25DOll1a735JkiS1ihDCu0MIK2RT2w4hXRzlqnr3S1L9VRpSKUmNbnHSEOeVSPPmfwZcVtceSZIktZb1SNPNxgDPAB+KMU6qb5ckNYKGmtoWQtiDtEZKG/D7GOMpde6SJElSyzODSZKkvBqmkBRCaAOeBN4LvATcDXw0xvi/unZMkiSphZnBJElSXzTSGknbABNijM9kK+xfBOxX5z5JkiS1OjOYJEnKrZHWSBoPvFh2/yVg264HhRAOBw7P7m5Zg35JkhrUFkAo37E5AztF8gw8OxXeGFCvWsfKwPLlO5bPdlZTBO6De8t3xRi6O1yFMINJkvqkUwbrb/56EpjRcfdZzGAl7RlsRdKqsNWW5S/Il8EaqZCUS4zxLOAsgBBCY8zLkyTVxX3AHGB4acdNpKXY++EDH4CP3guXAxcU0bkW8CXgy+U7DgZ+WuUXjbBgAQwdVuXXUZ+ZwSRJJZ0y2P3Am3SfwdYHnui+rQ8AH8UMVu5LwJdPBo6ny1nTKpkLjIAF5CsSNVIhaSKwStn9lbN9kiR1ayTwdvaRJUinssb18IQtgAc679ob+GdMl/7zr+NuHA+cXIPXCTAkwHwaK6S0ODOYJKnPFslgAK8BS2e3VwdeoMdwtTfwT8xgFZ1AKvCcVLuXHEK+DNZIayTdDawTQlgjhDAcOIhUlJQkqVuRdF3imaU7S5OuO9UGvJyOWQd4YoNs3/3puJ0itGXbP2NHW+pGoDZnxGr7UkrMYJKkPlskg0VgOTpy2POwTlx0MNJOZYf8s6wtVVDjL0zeDNYwJ/tijPNDCEcB/yH9TP0xxvhonbslSWoCC0mDkN4AFl9Y9sCqQEiPbzS/8z+M82vXvdZwKumL9uN6d0RFM4NJkvqrUwYr7ejy+EaYwfrlh8Bx9e5EZQ1TSAKIMf4L+Fe9+yFJaj7z6RhJ3W4BTOi4qYH4KinQVNt8WDAaRtXgpdTBDCZJ6q+KGQwz2IB9m/TFO7F2L7mAfBmsoQpJkiQNxLwK+9YlTS9X35wGHFW+o42qp4a5wFJtMH0qzFusuq8lSZKKYwYrTnsGWwB8l/TF/V4VX3A2zB0HSwHTqfy97KqR1kiSJKlwc3DefX98nbRu5s9KO34CfKP6rztzAYxdrvqvI0mSqssM1j+dMtgC4EdUL4O9TRpONietdTU259MsJEmSpEXMI2WLE8gu1jaflGi+VvxrvQ0sA+m05ZLw9tvFv4YkSVIzqHkGe7vz/TwsJEmSpG4dD3y5dGce8CvgS8W1/xawFh3rKPBWcW1LkiQ1q4oZbBzwhQIanw5vjYO1VirLYH0QYmzewWYhhObtvKR+mUI2ckFSTYwmrfU4HJhKNkV/BOl6vx8Cftf/tt8E1gDemAVLrZz2vfHGosfFGPNciVY1ZAaTBh8zmFRbFTMYpBx2CP3PYFPhzbVgjanpantLZbsrRLBuM5iFJElNZSEdQyknkS4lOh6vBiFV05js45HA+sDjwE8hBZmDgd/3vc2pwHrAa28BawKvdX+shaTGYwaTBh8zmFR73WawjwEHAocBuwIX9tLQRGCLdHPqQlhvSo/Rq113GcypbZKaysplt7fINgOMVF0zs+23wFeAGdn+v8+Bwy8APg1cT/qL4oDKbbwAbF26MwUWjofXxgPrkC/JSJLqygwm1V63GexSOPxTwGTgUlIGGw/sU/bkJ4F3wgsvwNZbZsdOhoU5i0g9sZAkaUA+RrqYE8D7geeBc6v4ei+X3Z7U5f6TwKgqvrY02M0kTUf7C7Aq8C/gh7OBi+GGg+BjL2c7V4WH9oY9yp47/wW4fztY9VXYfDPSL+/LwCs1/RQkqWWYwaTBY5EMNgt+ODU9dsOcLIO9DFwND60Ke6wKfBY4FeZvD/e/mp63eUH9sZAkqU+OJlXDSxYjXTES0jzeJ4B/A3/qY7uLAw/14fgJpCHV5XYHZvfxdSX13VvAi8DFwHbAWrPgd69lQ63npAfXuwZ+vxY8tRbsthbwUVjwC3hxG3hoImxZt95LUnMyg0laJIORlkr6aemAubDei/D7F+GpO2C3jwIvp9GDL5J+14vIYK6RJKlPlibN0Z1Hmmr7U9IZqEmkIDIWuIJ0lYEX+9DuENIb4VM5j18ItAGlN4GHgK1IVw+XVHtj6DztYWPgB6Sz5C8BqwyDZ1YFnu57266R1HjMYFLtmcEkVdJrBgOe6Wfb3WWwof1sT9Ig9TopQCwLvA2skO1fArgauBLYk94DzPKkNeF2ze4vJH+AAdiEFGAeIAWgT5KClaT6mEk6G17yPPAg8Gx2/5l59KuIJElKzGCSKuk1g1XhNZ3aJvXi26QrXAMcQ5pqWslewI9q0qPa+hbwkS77/gicVuHYrYFjgR2AM7tpb1XS2bI3SJewvBe4oYfXD8A9FfY/kn3cGPgM6c3S0+NS45hNmv4gSf1lBjODSeq7WmQwC0lSL/4M3Ax8CTiOjrM/5fYFfkm6+uIpwI7AHcCva9THavgeaTjk+cCNXR57lc4LLJYsSbqK9xLAR0lfg9+TLlV5fnbMCGBD0pmrJ0hX/HhH9tgo4KYubUbgc9ntO7Kt5HbS1/pXpCHWkiSpdZjBzGCSGpNT26RePEta2PDLwN+BYaQzYm/T8Y/r8qQhwbcDmwLjSGdx7gD+CzwGHF7LThfgfODzpHBwBXAUcGDZ4xcBfwPuJ4WdE7o8fwlgW2AD0oJu65C+FqOAlbLbw4GpwEHAVVk7x1foy93Zx6902f8V4E7gq3j5WUmSWo0ZzAwmqTG52LbUg28A7wXWBcaThgiWqq8LSAsTlrxCmru+HGnRw5VJ81M3J12q8SLSmbKTgE9Vv+uFWB+YDnyYdDWQ8jNRL5KWO9kb+DEprADcSgo9p+Ro/3VgP1LIe5p0Nq0319Mxp7/kGtKaAPNzPF9S83Gx7cZjBlO1mcHMYJLqr7sMZiFJquDbpDnm7yCduSnCm8DDwGbALcBk0uKEzeB00iJt00lz9S8CziUNJ/8RaZh0yWRSwMlzWcmXSeFwKLAHcCjwgV6e8z7SgpLl3gtci/PzpVZlIanxmMFULd1msM/SsWBSuauAX2S3dwROzG4/SJoPhxmsEjOYpDwsJEk5fI909morKs/Db3cgcHCF/TcDP8lub541+D/ag0y5t0lndqbSuGHmWNKZqpdI/T2ctFDjY6QAcj1paPNR/Wx/Nil8zCB9Dd5DyoOSVM5CUuMxg6loFTPY0aRKBaRhN2st+jxeIF1/HtJcuKmkatQU0jCeJ0lhZTXS3K3PDcIMtjrwNeALHbvMYJLysJAk5XA1HXllEZ8hXRYE0njjDSsc8xJwV3Z7EmlFxGPoWL3wBdKKkWVeAVbsZ3+rbUvS0Of3kD6VZ0jDyT8ALEUKMe8GxtCHYtKKpLOFn+/Y9TopHP0ZmAN8LGdTfyUN+faNQGptFpIajxlMReuUwb5MulzYJ0lDZfKaSJrTtkPZvjdIK1aPJc1526DjoUGTwcYC25AqR12YwST1pLsM5mLbEvBd0lDhio4gTR7fBDigl4ZWzjZIhaTHgKXLnjeNjtNsU4AvpkUhz6IxF4K8N/t4H/AW6eoe84B/kK788U5S/+8gDRW/h3QGrUczgH913J1Oqq39ifRlersP/bsEA4wkSc2sYgbbhlSlWK6PjY3PtnJLkQLGC6RhT7+gPayMo0Uz2Makz/lC0uf8PBVX0jaDSeqvIfXugNQI3gOcBmxU6cGtge+TJob3xYosuiLhONLlMXYhDcP+RTr7U+EEUUM4DNiZ9HVZg3Sp2I1JWWT9suP+C1xOCji/6NrIWsB3gFWA80ijka7seHgEaaDXB0nDrD/bh/5d0odjJUlS4+mUwb5MygrbAtuRrmdflMVJw3xu7NjVshmsNK1vKumk5h6VX8MMJqm/LCRJmd3psqjjEcDZpFM+29FpKPSAjSJdo/ZOmEX9/jH+FimcdOcJ0sCqCaTR4muQFndchjTE+nrSUOjlSVfsWB24oWsjM0jXp52ZNbR/54fnkq748UVSCLqwn5+LJElqTu0ZbFfSdLbVq/AiSwIfJY12yrRsBnuFNFdwCvAb0rCvCsxgkvrLqW0a9L4O/J30D/O65Q/sQlpUuxrGkq7ZelS6eSoV1+OuusdJ9Z1KPk1aJ/yJsn0fJw2tvhc4pGz/+sBipMxyBPAocB1pjUwmA5dlB55BmiaYmU4ahn0uMJx0NmwYab3yY/vzCUmSpKbRKYN9kTTkpprGkabMnZDutnwGg1RMunTR1zCDSRoIC0ka9D5EOrMzurTjMNJC2ptX+YWXhGknwXe+m/7xroe/drP/U6SZeXdUeOwqFr2qx5OksPM06ezeKNJVVxYxFfhLx902Oi7AMhf4PelNaULvXZckSU2uUwbbh3RltWqZRKoaTU93p5Fm3g+aDNaFGUzSQFhIkugy8OgN0oT9dSsfW4jXYeqP4Zsz4MwqvkxvvkhatPHFLvunAneSptWX7EWaqz+fylc4WZG0rvg1pGD2wRyvvxhpOYT5pGx3AumKvWfk674kSWpyB0Ja7Xod4BzS2pRFjkx6gTQiZ1fgl2nXVOCbmMHMYJL6yzWSNOj9lrT+87OlHaOoXon1VVKImQpv/hjOPDNdtrXChTRq4m1gYYX9l9E5wEA6WzWXjiWj3irb1iENC98PWJs0arzkBboPajOAX5HCVMzakiRJg0N7BpsF/AE4GXiq4BdpI41CKquQvEnKJmYwM5ik/nFEkga9P5KGEXfyf9nH/YF3FPhik2HKD+Ci59II7nr7Qw+P7U8abv5Q2b7rqXyp1weBkN1+krRY40dzvP4M4AekBR6PBn6c4zmSJKk1tGewPwMfI13Zdr2CX2Q8aV3KLdNyQRdhBgMzmKSBcUSSBr1PkM7gtJV2nA/cRzo1M29gbb8AXFB2/7Vl4az9YOIv0lk4SAstnjKwlynUAaRZfaNIiy5Cmul3ImkdgR+Q5u3vSgo67yRdKaQ01/9x0sVBLs7ut1G2/lQXi5EWlPwR6aygJEkaPDplsFGkYHAFaa5Vabu/j41OSM974UK44KWsjfPgNeAs0hXQzGBmMEkD44gkDXrLkBYYnAR8jWydx6VIE85vzjaAd9PnBbiffhHO/BvstjxctBtMugjOXTat5/3DYrpfmAOAVYDlSFcAmQk8kD22H2n49QjSGay/ZY8PzY7dkJTb/pkd/wQp9EwC5pCWnTofOLjs9WaQgs6JpNHm36nGJyVJkhpWpwz2h27W2h7LovnrHuAWYBNSVaXcHOBOePr3cOaesNvpaRTSJNIVysxgZjBJA2chSYPeadnHX5OG6K0P7PFHWLPrgR8irUTYNcw8Av97DW7cJZ1Ve1+2+wXgslkw6Wb4+X1wzkdht1PSHPfnqvGJDNDKpEUXjwKWJ32a5WexLiMFnBOBu0hnudYHNgJuIp3NmlR2/FPAl4A1gI+Q1iB4M3tsLKku91VSwJEkSYNPxQxGlwx2Mx1zurYhXZLsOlKw2An+F+DGlWDt6+B9a6ZGXhgPl10Kk06Hn5PW8N4NMxiYwSQVw0KSVOb07OPn6agX7UQ2Xb80bnj37OMmpEBzM9x2IXzhEDhwTVh3Tbj2Wnh0PPxyU2Bd+MmlsMopsAHpH/yNavLZ9M2vSHPsXyCNKt+ZVBTbn3SJ3EeAxUlnD0th7yPA74AvANsCm1Vo91nSqHKy48ja+Bpp5uDDxX4akiSpCXWbwS4iDSmCNDxnb+Du7P5/4bZvwxfWhQP/COtuDte+Ex49o/0CbfyENNrHDJaYwSQVIcRYadm25hBCaN7OqynsSMorK2X3t6Rj2aSh74f7DoT1roOH/5SGSm+7LWyzDZx+OunStXvTWJPve/Be4FHg5ez+50gB4xOkxR5/TQp2C0iXoT2ZNBf/VNJQ8TyWIIWjy7L7K5KGbi8/wL5Lam0xxtD7UaolM5iqrccMRlrOcj1SVjmMVEzZho6CVDMxg0lqVN1lMEckST3YDPgPcGN2/2PArOz2yH/Bhf9K//ivke278860Aen00SM16mg/vQ+4F3idFMAm0xFiSpeL3YB09mpdUmBbj3Qm7HHg7/TtU1yMFAoXAFeSvpb/GMgnIEmSWtJm9JDBSFcn65TBsq1ZmMEkNTNHJEk92Ip0lY/n692RKvkJ6QomT/Vy3GeBV0nDrU8iLcr4U9KlaFclDa3Oa0XSHP/P97GvkgYvRyQ1HjOYqs0MlpjBJNWTI5KkfngXaUhxq4aYr/Xy+KbAVNIVVSp5jDRXfy/S5XQfyPGakzDASJKknpnBzGCSGpeFJKkHP693B+rsPcD/SIs/lnuGdOWTA4FPkubdn06+ECNJktQbM5gZTFLjcmqbpD45hFSBfhNYgXRFkA2BH9ezU5JamlPbGo8ZTKo9M5ikWusug1lIkiRJDc1CUuMxg0mS1Pq6y2BDat0RSZIkSZIkNScLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXIbWuwNqLCsC04CxwBjgVWBmPTskSZI0CJjBJEnNwhFJ6uTPwIeyj08D+wBLkQKNJEmSqsMMJklqFo5IGqSWA94A5gPLAMOz/SOAP5UdtyRwJnA38Bvg7Rr2UZIkqdWYwSRJzS7EGOvdh34LITRv5+vsduBY4FngGmCTHM/5BvBH4C0MM5Kk2okxhnr3QZ2ZwfrPDCZJahbdZTALSYPY/4AN+vG87wGnALOK7Y4kSRVZSGo8ZrCBMYNJkppBdxnMqW2D0OKkxbH6u0DWt0nDsE8B5gCzC+qXJElSK+s2gy0GtHXzpFnA3HTTDCZJagSOSBpkxgCPAqv1dNAoKoeZOcC8zrt+CnytmK5JklSRI5Iajxms7ypmsNGkatCNwE7dPPFw4OxFd5vBJEnV1l0G86ptg8go4Al6KCKNIP1E/B2YUWE7iopj2IbSsVCkJEmSOquYwUYCt5AKSD0l8uGkgtNoFslhZjBJUj1YSBpE/gcsDVQ8hTgcuA54L93/VPwc+BKp4JSNWBoKfBH4VaE9lSRJah2LZLDhpMuxbU4ajbRjD08+A5iZbV/p2G0GkyTVi4WkQWRdYCzwfNcHhgI3k0LMVcD7emjkx6Qh2CcBbemqIz8n/SB1N7VfkiRpMOuUwYYBDwIbD6zNYzGDSZLqw8W2B4E2IAAvACtUevBOYIs+Nnpi9txvpruHkX6YPtP/bkqSJLWUihnsUWCdATRaOg28MH0wg0mSas1C0iDwCrAk3Qw/uw/YZACNZ0tvLYzteUaSJEn0ksH6Ywjp0m3DaD+ZtxAzmCSptpzaNggsS8obrxTdcAC+DPwSfk86IyapcXiZK0mqr04ZLDDwN+aTScsLlDGDSY3HDKZWZyFpkHgVWLHoRr9BuvYs6cq0fyi6fUkD8ioVprNKkmqqPYM9DyxgYNPaKjCDSY3HDKZW59S2QeB10rDqTp4E1qawcnnFK8FJqrkZwGL17oQkCeiSwVYFngVWL/Y1zGBSYzCDaTBxRNIgsDQwpevOdYHHi3uNs3CRR6nW5pP+gCjfugaYSWWPrV3T3kmSKmawgpnBpNozg2mws5AkSU2mtLBqXy/3/CSwfvHdkSTl8QKwWr07IWkgus1gM7MHxlZ+nhlMrcZC0mBV8LvZ4cAfi2tOUhcLyz72ul7rtOzAZTrvLnggoiSpL4pYbPv4rI1vduwyg0nVtUgGW0BHNWlOdntU9uA0OjLYy8D49FwzmFqNhaTBbKBh5kekn6CjvTKBVC2BtDZr+cdO5pIuCQTwVnbQ2OzAV7P744Gn4cm10913VL/bkqRq+BHw7c67zGBSdfSYwQIdGay8SFy6/SqwJTAx7X4SM5hai4WkFvcGFRbaLspxwKnValyN6GQ68uuXSfPDz6xfd1peGzCPjjfqim/YAZhF+maMzg4KZU8YQppO8R4IE3zTl6RaqUoGOwH4ftGNqhmYwWqrUwabD0Pmp4/Mp3PO6q6SOwR4qeM5YS0zmFqLP88tbllgarUa/zGpmERa6PGz1XodNYxvkv69nAf8hPSP7OHZ/XnARfXrWssYScfXczYV5uDP67K1lW09hZknO55z/8aweeE9lySVWySDrQ48N4AGv0o6gdflMm1msMHBDFZ9i2Sw0he8jXx5q6shiz7nfsxgag1D690BVdeC7ONKwIvACqUHNgQepv/rJH0fTv0OnJiFmdLCc2pdPwaOoePfw5LSoBeAD5PeVD5U2661nG7fmOf09GC+Rt8BPBFSLpIkVc8iGWxBDwf35HPAH7IGs9x1KnBi9rAZrPWZwWpnaOl/M8t3DNA6wIQsg2EGU2twRNIgMIm0dMry5TtL16zsq68DI4DvwlcWptHV8+gIS2pdbcBwer5S2BBgf+AyYBfgmup3q2UsDrxe6YFAOi02m/QNGKB7N4X1Hhl4O5Kk3i2SwdYj5ai8q+5+ljTkaB6dqkVfwQw2mJjBqqtiBhtOIbkLgEeA2XDveuktQGoFIcb+VBMaQwiheTtfQyNJZ8K2Bm6h/eIBKcjcA2ycs6EvAWfAifNTm8NJw2ytqre+nwNH0bGmc28WkALMfqR1CNWzJUgzzzYk/a6OKn8wUOyp5tkweyvY/lF4oMBmpWqKMbqecIMxg+XT7wz2MeAfpH9EyypFJ2IGG2zMYNXVKYONhFGvk7LXqB6f1j8bwOzHYXvMYGoe3WUwp7YNArOB1Ugz2VYsf2AOsA1wKz1P1j0CuCAdf/z89A8adMzTVusbRv4AA+mM2VAMMHkF0loaz1GWW4ZSlQXOthwJT9wFs3bEFCNJVdZjBis/SfA+4LYuT+wy1Oh4zGCDkRmsutoz2GgY9TywMmml/KJtCls+kaa2zapC81KtObVtkHib9Adq6Ru+OemNjFnATsDdZQcfQCrPl7ZzSPOE56d/lEqL0PkP1ODwS+CwfjxvZ+BfxXalpQVgsdKd4cArwKrlO4txI7Dae2DhQ8W2K0mqrNsMtgMdWet6UtYqbRXmq5nBBh8zWG0EYLG3ScOSnqvSi7wNN8ZUWHZNM7UCC0mDyAbAUtn2R+Dg7PaDM4Hd4P1Lwc1LAVcAb5Zt8+DY7Njf1KHfqq9vAH/ux/NuBQ4suC+tZlnSnPxnuj4wl/QLO6H413wv8MRbmGIkqYYqZrCZ8GCWtd6/AG7u5rnHYgYbrMxg1dMpgy2W3XkCGFu913xv9hJSK3Bq2yAyvez2bqQa0QLSWYuhM9I/VD8mrcP9R9I/QPuQ5mX/gTTKWoPP94GP9+N52wPnk+boq7IppK9Tp1AxgrSgBqS/HAo2HRdmlaRa6zGDYQZTZWaw6umUwQJVyVzttgSeMYOptVhIGqTKl16Zln38OOnKIreRhmH/k7Sg43QMMIPZd4FxwGf6+Lzh2fNU2QrAfXQZFjqKdGps2Sq+8LuAx6rYviSpR2Yw5WUGq46KGayaXseR4Go5FpLUrrSuXGne/Zx6dUQN5U3SZYbnA4fXuS+tZAhdFl6F9NfCO6nKlLbtSQOdXn2N9M2UJDUMM5gqMYNVR3sGWxJ4iDQiqVq2hu1fyjJYFV9GqjULSZJ6dDLwUdK/taqiMcAjVC3MTAIm7gw8XZ32JUlSscxgVTYd2At4sErtbwPcD5MWwMQqvYRULy62LalHvwBuIV1URsVYCfhv151DgNVJl/Mo2LuAibuQ5kx4vWhJkprCLzCDFa1TBlsIvFDFF3sR3mURSS3KQpKkHk0Gvgac18fn3Q58svjutITXgL1JBZ52M4GtqvN6z+wK82/FIpIkSU3EDFa8ThlsHOmLVUXP4IoCak0WkiT16hXSJWj/2IfnzASer053mt480nrXtwPvLu0cBZxTpRc8g+ou4C1JkqrCDFasThlsBnBYFV5kJ2BT0qXhpBZlIUlSLpOAbwNn1bsjLWQ+8Gjpzizg88W1vSvpykC77w6TP4lhRpKkJmUGK9584NEFpPUpC7brozD1Idh9fhpVJrUiF9uWlNtE4PvAAuBzPRx3O3B0TXrU3JYC/l26s5B0LdqCnAp8BLjpIZj3SnHtSpKk2jODFWsp4N9LAFcX2OiuwCw49c0sg+GqAmpdIcZY7z70WwiheTsvNbFVSOtCfxA4BvgXcC9wImlRyM9RlRM8LWUZ4HJg+9KO0cB/gHcW0/7ee8PVb8K8OzDFqOnFGKt5cWb1gxlMqg8z2MC1Z7ClKXbE9mLAzLQG09UYv9QaustgFpIk9dvKwHrAy8AbwMakRQwfqmenmsRKdLmKx+LAm8W0vS/w76Vg/tRi2pPqzUJS4zGDSfVlBuu/9gxWhULSvjPTaHMX2Far6C6DObVNUr+9lG0lr9arI+rkJgwwkiS1MjNYYzKDabBwsW1JqrFl6dvVVyRJkjRwZjCpGE1dSFqz3h2QpH4YBexe705I0gCYwSQ1IzOYVIymntq2JHB+l32fJF38SJIa1RTgKOCMendEkvrJDCapGXXKYG8BhwNn1bNHUnNq6kISq8HHn++861MYYiQ1tlHATvXuhCQNhBlMUhPqlMHmAP+HhSSpH5p6ahvLAL/NtsxvgLY6dUeS8pgD3Nd152zg6Nr3RZL6pTyD/RYIZjBJja9iBpPUZ809IgngiOzj3PThsGNh7kI4FlfMl9SYhpEu2dvJPOAPwK9q3h1J6p8jym7Pg8OOMYNJamyLZLBZwNeAn9SlO1LTau4RSeW+mG0/gS+0eUZMUuOaT+dL9kpS0zsKGAJfwAwmqXEtksHmkE7inVCX7khNq3UKSSVfBn4I3xsKPwSG17s/ktTFDNKJr5O7PjAP+HbNuyNJhfoeZjBJjaliBpsL/JSBZbCTaJ8hIw0GrVdIAjgOvt6WCsuGGEmN6C0qzGKbB/wI+H7/2z0ZmPMj0lBtSaqlE4GTzGCSGltVMthP4OR5aYCTNBg0/xpJ3fk6aeziz7A6LKl5zAe+S5rEf3zfn34yMOcU0uLdklRLpbP5pwAL6tkRSeqH/mawHwPzsgxWjX5JDag1RyRBGld9Mhw9HEbUuy+SVMFM0kjqM7o+sAD4JvDzvrX3S+Bzv4LhphhJDeBozGCSGlMhGeznWSMnwC/nw+dwJKYGjxBjrHcf+m2rrUK8555eDlocxr0F02vSI0nqu3HA1EoPBNLY66NyNPIbGL0AvvdVOHGuA5LUWmKMod59UGc9ZrCRtJ+WH4cZTFLjGkcPGeyX2e2jSBWn0seSY4DsT+nRpHEMJ2IGU2vpLoNVbWpbCOGPwN7A5Bjjxtm+pYCLgdWB54CPxBinhhBKv6rvB94GPhVjvK9afZOkphBJp/SHle07HDgr+3h2OuZsYOExMH9+uoKtpMHNDCZJA1TKYJCGGR1d9rHM2cBC0qw4M5gGk2pObTsX2KPLvuOB62KM6wDX0TH7dE9gnWw7HDiziv2SpIYyFzgH+DMpiJxX/mAEjizbzs0+ntOx7wtHwpHz0zqRkoQZTJJy6TGDlRzZ5WOZL2S7zWAabKo2IinGeHMIYfUuu/cDds5unwfcCByX7f9TTPPs7gghjAshrBhjnFSt/klSo3gb+AxpLZEFwGF0vDkH4GPlB38mfbjg0I5dC6veQ0nNpK4Z7EJ8U5LUNPqUwTIXlN327U6DVa2v2rZ8WTB5BVg+uz0eeLHsuJeyfYuEmBDC4aQzZqy6ai+v9ne41LP0kprEHODT2e2Ds4+BzjPb6PK4JOVUmwx2CO3B61LMYJKagxlM6ptaF5LaxRhjCKHPK33HGM8irRDCVlv1/PwrPwEfm+2CZ5KaVwQ+Uu9OSGopVctgl9O+8OyVpDP5ZjBJzcoMJnWv1oWkV0vDpUMIKwKTs/0TgVXKjls529dv//kPnLXAM2GSJElUM4P9O/t4ALAA/kOqNpnBJElqTdVcbLuSy0mDnsk+Xla2/5Mh2Q6Y3u+5+dcB18J++8CFc2DUQHssSZLU/KqXwfYiXfNtQbq7H2mpJDOYJEmtqWojkkIIF5IWdVwmhPAScBJwCnBJCOFQ4Hk6Rgv+ixRBJpDWPPv0Ig3mtQf8d34aingr7ZlGkiRpUKhJBnsLuHnR3f/FDCZJUqsL6SIdzWmrrUK85x7gdrgd2A4IO8Fwp7RJktQyYoyh3n1QZ1uFEO/JbrdnMGA4ZjBJklpFdxmsbottF+Jt4B5gB9gRuIM0V695S2OSJEnNxQwmSdLg0tSFpFmPwYNbd9zftn5dkSRJGjRmAQ+W3TeDSZI0eDT11Lb+XLpWkiQ1F6e2NR4zmCRJra+7DFbrq7ZJkiRJkiSpSVlIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkqQB+gRwYb07IUmSNMiYwaT6sJAkSf30OWAhcF52/wPZ/auAnYH/1qdbkiRJLc0MJtWXhSRJyuFU4JvZ7a8BC4BfAyHbDgT+mt2WJElSMcxgUuMJMcZ696HfQgjN23lJTaUUTk4AvkfPVfhIOis2JPsI8AiwWbU6J7W4GKN/HzQYM5ikWuk2g30GOAu4GPh42mUGk4rVXQZzRJIk5fAjYA7wXbp54/wIMBf4Rwo8bXR8bMuesy7weA36KkmS1Co6ZbCvkvLWXFIRqQ04qGNfmAttpY+3msGkarGQJEm9+BnwVWAYKZC0O4KUbOYAF2QH7FO2r7TdCRsBD3d9viRJkrq1SAYbkt0pD2Xl+4Yt+rgZTCqehSRJ6sGvgBeyj4sYAgzPtqEV9pW2rWDIPenmmsATVe6zJElSs+uUwb4NzAJ+0IcGtk7PGTILhj9gBpOKNLT3QyRp8BoK/JguCzh+Cfg++d9BhwCbA2/BkGdg7U3gSdIwa0mSJC2qUwYbCozsYwNDyp7zDhjykBlMKoojkiSpF8NJI6TbDQPGACP60MiQ7DkbwZBHYC1gQmE9lCRJaj3DgWHfB74+wIaGAKPSBzOYNHAWkiSpgt8CU0kXBGl3QrbzOwNoeAiwAQx5FNbAICNJklRukQw2gr6dvOvOGsBDKYqZwaSBsZAkSRV8FfgnXXLLacCZwKgBNj4EWMIgI0mS1FWnDPZD4KiCGm4DxqabZjBpYCwkSVIFbwGzu+48GvhiQS+wAkx+OF1JZKmCmpQkSWp2nTLYGAZ+Aq/cysADMBkzmDQQFpIkKa8xwGIFtdUGSy8D1wNLAI8X1KwkSVJL+BFwaMFttgHLwNKYwaSBsJAkSXmcQBqRVKC2ZWHF+9Mb8YrFNi1JktTcfgxcUoV2l4e2e1L2MoNJ/WMhSZK6cTxwQenOb4E/F/wCbcDKMB3YtOCmJUmSmtXxwAVTgRlVaHwoMD7dNINJ/WMhSZK6cSzw/tKdN0iXECnaOFh4JzxXhaYlSZKa0bFkGexk4A/Ve52FmMGk/rCQJEndOBO4utovMhTGrgO3V/t1JEmSmkR7BpsMfBP4XXVeZyxmMKk/LCRJUjcmkgYitTsT+E3xr9MGbFh8s5IkSU2pUwZ7lVRQKspk4L3pphlM6h8LSZLUg1OBi0t3JgE/BM4orv1p0+A97ymuPUmSpFbQKYP9juJGJc0DHoFpgBFM6h8LSZLUg+eBV7LbZwNnvwy8WFz7Y+bDcffB7sU1KUmS1PQ6ZbCJcPZLxbY/BjgOM5jUH0Pr3QFJanRnAMsDu5V2XAisCnxh4G0PA3YGvgXsPfDmJEmSWsYiGWygXgE+nG6awaT+s5AkSb2YQBqEdFBpx4vAUwNsdArwcWAuzAF+PsDmJEmSWk2nDPYX4K7sgb2BL/ajwTl0Wl3bDCb1T4gx1rsP/RZCaN7OS2oqGwBrkoLMwQBrkFZnfB9wdD8anAX8Fd78ZDoxVvWrw0lNLMYY6t0HdWYGk1Qri2QwgNWAL5M/g70EHEnKX9enXW9iBpN6010Gc0SSJOXwWLZNABYCn3wWeBZ4JDugL8WkyfD6kXDYHJiLAUaSJKk7i2QwSAsoPQTcCNxHKip153ngs8C16e7rwGGYwaSBsJAkSX3wBPAjIAKHQAonp2U7ticFmq/30MArwGEw6sq0uOOR1eysJElSi1gkg91AOqG3OmnK242kcPUjYD3gu6STfp+nvYgEMAozmDRQFpIkqY8eJ12SFtLI6ieeg01/AXddCkcvBuwJXAl8EJ7+C5z7Hfg+wMukBbqvhNl0yjSSJEnqRacM9gw88Qxs+irc9XU4egJpiNF1wLLw9LNw7hvw/as6t2EGkwbONZIkqZ82AJYmDTI6ILs9chU4eiPgQWBLmHwnXLgvPAysMB2O/iucALxNuvibpN65RlLjMYNJqqeKGYzOKw1MJmWth4EVssfMYFLfdJfBLCRJUgHWBtYiDfO8EngS+EX22Khs+zWwF3BBHfonNTMLSY3HDCapUZjBpOqxkCRJNbAGsC/p4iB/y/YtS5qLf369OiU1OQtJjccMJqnRmMGk4llIkiRJTclCUuMxg0mS1Pq6y2BDat0RSZIkSZIkNScLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJyqVohKYSwSgjhhhDC/0IIj4YQjsn2LxVCuCaE8FT2cclsfwgh/CqEMCGE8FAIYYtq9U2SJKlVmcEkSVI1VXNE0nzgKzHGDYHtgC+EEDYEjgeuizGuA1yX3QfYE1gn2w4Hzqxi3yRJklqVGUySJFVN1QpJMcZJMcb7stszgMeA8cB+wHnZYecB+2e39wP+FJM7gHEhhBWr1T9JkqRWZAaTJEnVVJM1kkIIqwObA3cCy8cYJ2UPvQIsn90eD7xY9rSXsn1d2zo8hHBPCOGe6vVYkiSp+ZnBJElS0apeSAohLAb8DTg2xvhm+WMxxgjEvrQXYzwrxrhVjHGrArspSZLUUsxgkiSpGqpaSAohDCMFmAtijH/Pdr9aGi6dfZyc7Z8IrFL29JWzfZIkSeoDM5gkSaqWal61LQB/AB6LMf687KHLgUOy24cAl5Xt/2R25ZDtgOllw68lSZKUgxlMkiRVU0gjm6vQcAjvBP4LPAwszHafQJqjfwmwKvA88JEY4xtZ6DkD2AN4G/h0jLHHOfghhOp0XpIkNYwYY6h3H5qJGUySJBWhuwxWtUJSLRhiJElqfRaSGo8ZTJKk1tddBqvJVdskSZIkSZLU/CwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXCwkSZIkSZIkKRcLSZIkSZIkScrFQpIkSZIkSZJysZAkSZIkSZKkXHIXkkIIS4UQlqpmZyRJktSZGUySJDWSHgtJIYRVQwgXhRBeA+4E7gohTM72rV6THkqSJA0yZjBJktSoehuRdDFwKbBCjHGdGOPawIrAP4CLqtw3SZKkwcoMJkmSGlKIMXb/YAhPxRjX6etjtRJC6L7zkiSpJcQYQ737UGtmMEmSVG/dZbChvTzv3hDCb4DzgBezfasAhwD3F9c9SZIklTGDSZKkhtTbiKThwKHAfsD4bPdE4HLgDzHGOVXvYQ88GyZJUusbpCOSzGCSJKmuustgPRaSGp0hRpKk1jcYC0mNzgwmSVLr69fUthDCUNLZsP3pfDbsMtLZsHkF9lGSJEmYwSRJUuPqbWrbhcA00vz8l7LdK5Pm5y8VYzyw2h3siWfDJElqfYNxRJIZTJIk1Vu/praFEJ6MMa7b18dqxRAjSVLrG6SFJDOYJEmqq+4y2JBenvdGCOHDIYT240IIQ0IIBwJTi+ygJEmS2pnBJElSQ+qtkHQQ8CHg1RDCkyGEp4BXgAOyxyRJklQ8M5gkSWpIua/aFkJYGiDG+HpVe9QHDquWJKn1DcapbeXMYJIkqR76O7WtvIHXgbEhhANCCOsX1jNJkiR1ywwmSZIaSY+FpBDCP8pu7wdcD+wDXB5C+FRVeyZJkjRImcEkSVKjGtrL46uV3T4O2DXG+GwIYRngOuDcanVMkiRpEDODSZKkhtTb1Lby+e9DY4zPAsQYpwALq9YrSZKkwc0MJkmSGlJvI5I2DSG8CQRgRAhhxRjjpBDCcKCt+t2TJEkalMxgkiSpIfVYSIoxdhdURgNHFN8dSZIkmcEkSVKjCjE279VbvfSsJEmtr7tLz6p+zGCSJLW+7jJYb2skdSuE8HD/uyNJkqT+MINJkqR66nFqWwjhgO4eAlYovjvN6whgM+BzwIeBDwEHAu8BvgbsXreeSZKkZmMGy88MJklSbfU4tS2EMA+4gM5XDin5UIxx8Wp1LI96Dqv+NvCNsvttpGQ3nzTMqw2Yl92+Fdi11h2UJKlFDMapbWaw7pnBJEmqje4yWG9XbXsI+GmM8ZGuD4QQ3lNEx5rRD0hnuIZXeKz8C+olVSRJUj+ZwSowg0mSVH+9rZF0LPBmN499oNiuNL7TSF+M6dntvHbKnvcm8N8q9EuSJLWcYzGDtes2g32ejpD1O+CD2e1L0sNmMEmSiudV23I6HTgMGAHMIg2hHtn1oE8Av8puX0EakH44KdRk5gN3ATtWs7OSJLWQwTi1rdE1RAY7GjgFGJUdOJc0p21M9vEKzGCSJA3AgK/aFkLYu6f7rew3wGdJAQZSXlmkiPTp7MBx2fYR0tmwvYEpwD/TYUOBbbJddwMrA156RZIkdccM1iWDHU8KUifTUUSCNN9tTHZ7GB0ZbArwHzOYJElFyV1IArbqcn/rIjvSyI4H/t7dg18AXgF+ASxWtn8EMJYUapYmXTrklbQNvTHtGge8DOxcfJclSVLrMIMBfJ+Upb5FClJjunlSSSmDLQ3sAlyTiklmMEmSBqbXqW0hhPWB/UgnbiIwEbg8xvhY9bvXs1oMqz4T2IcUODrlla8BxwCLkwpGffEwsEkaYv0qKRN1TYiSJCkZrFPbzGBlGexHwFF0PmnXV3OBW4DdzGCSJOXRr6ltIYTjgItI09HvJE0tD8CFIYTji+5kI/oW6RPvVEQ6AfgmMJ6+F5EA1gOuT2fFIrDvAPsoSZJaixmsSwZbgoEVkSCNUFo23TSDSZLUf0N7efxQYKMY47zynSGEnwOPkpY4bGknAe/tunMcKdD013DSSo/Xw4Jd09BqSZKkMmYw4L0/A26vTvsLMINJktQfva2RtBBYqcL+FbPHWt5ypNlrQEo0nyio4eHAKumLe2NBTUqSpJZhBgMW/ynwbuCgghpdl/YLoJjBJEnqn95GJB0LXBdCeAp4Mdu3KrA2aab64HI26SvyqYLaWwWG/Qd22B1uIK0DKUmShBks+SLwcWDJgtobAayZbg4DdsAMJklSX/VYSIoxXhVCWJd0tdTx2e6JwN0xxgXV7lzDOSTbli2ovRHAeinIrF9Qk5IkqfmZwTK/AzYA9i+wzdWBy4D9zGCSJPVHbyOSiDEuBO6oQV8a31+ATYED690RSZLU6sxgwPPA6wW3ORLYsOA2JUkaRHpbI0nlniddJ7ZIywOXwlLAlQU3LUmS1PROAf5RYHvP02nNJTOYJEl9YyGpF98DrinfcQbpYrxFmQb8OK29vW2BzUqSJDWz9gw2AZhUYMPLAF/quGsGkySpbywk9eJ/wKvlOybQseRlEZZgMC2ZKUmSlEunDHY6xY1KegP4bUFtSZI0CFlI6o9zKG5U0gzgjwW1JUmS1IoeA74F/L2AtmYCtxTQjiRJg5SFpBx+CVxXvuMx4KmCGp/VtXFJkiRBlwz2KHAS8H91644kScJCUi73AN+gy1pJl2RbQaYBhxXXnCRJUtNbJIM9Ajxe7GtMwwwmSVJfWEjK6W7S8kjtHgF+AFxcTPuzKfaCJJIkSa1gkQw2EC8Bx3XeZQaTJKlvLCT1wfnAZ4GrSjseBu4bQIOTgeNgKvBd4LQB9U6SJKk1dcpgl2V3Pkvf16ycCtwFnNxx1wwmSVLfDK13B5rJbdn2PyACewJcTVowG2AH4OA+NPgmcA0M+yFs8E14oLiuSpIktYxOGexe2PPe7IHFgINyNPAccCGwN2kI0qNp9zBgA8xgkiT1hYWkfrgd+D6pmPT+B+hIH3PJV0iaBJwLfBjemAffexbeJl0MTpIkSZV1ymAANwPHlh2wOXBIhSe+Shp29D/SokgXwBvA9zCDSZLUVyHGWO8+9FsIoa6d3550JdrFgSnABzYB9sge3Aj4ZDdPfBjYGfggPH82rF7dbkqS1NRijKHefVBnDZfBSg8cmD3wZ2AdYBfgLNJJvD93buN5zGCSJPWkuwzmiKQBuJ203vb7gG0BHso2gPVIK1BtB5wHrEqaC/c7YDK89gb86myYXutOS5IkNblFMljJw6RFj/4KrAncRFpgqcxrwK8wg0mS1F+OSCrA9qQa0eakqfft1iUVkv5EKiTtAZyV1tj+ISnESJKknjkiqfE0fAbrhhlMkqT8ustgFpIKtDXwQdKijft2c8yrwKl4dRBJkvKykNR4zGCSJLU+p7bVwN3ZtjkQgH1I8+/vBj4EvAL8FAOMJElSkcxgkiTVjoWkKrgfOBF4EXgM+A9pKPVzwM/q1y1JkqSWZgaTJKn6nNomSZIamlPbGo8ZTJKk1tddBhtS645IkiRJkiSpOVlIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpRL1QpJIYSRIYS7QggPhhAeDSF8N9u/RgjhzhDChBDCxSGE4dn+Edn9Cdnjq1erb5IkSa3KDCZJkqqpmiOS5gC7xhg3BTYD9gghbAecCpwWY1wbmAocmh1/KDA1239adpwkSZL6xgwmSZKqpmqFpJi8ld0dlm0R2BX4a7b/PGD/7PZ+2X2yx3cLIYRq9U+SJKkVmcEkSVI1VXWNpBBCWwjhAWAycA3wNDAtxjg/O+QlYHx2ezzwIkD2+HRg6QptHh5CuCeEcE81+y5JktSszGCSJKlaqlpIijEuiDFuBqwMbAOsX0CbZ8UYt4oxbjXQtiRJklqRGUySJFVLTa7aFmOcBtwAbA+MCyEMzR5aGZiY3Z4IrAKQPb4E8Hot+idJktSKzGCSJKlo1bxq27IhhHHZ7VHAe4HHSGHmQ9lhhwCXZbcvz+6TPX59jDFWq3+SJEmtyAwmSZKqKVQrJ4QQNiEt3NhGKlhdEmP8XghhTeAiYCngfuDgGOOcEMJI4M/A5sAbwEExxmd6eQ1DjiRJLS7G6MLPfWAGkyRJRegug1WtkFQLhhhJklqfhaTGYwaTJKn1dZfBarJGkiRJkiRJkpqfhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIGkReBhaWbevXtzuSJEmDghlMktRKLCQNAq+TQssKQCjb/kfnUPPOenVQkiSpBZnBJEmtyEJSi5sKLElHcCkXumySJEkqxiIZ7AVgAbC2GUyS1NwsJLWoGcB8YAm6CSiPs8i46puy57yvul2TJElqWYtksFeyHSuTkvcTpAx2f7Z/Ptz0TjOYJKl5WEhqUUOANioUkZ4F5gHrAo9ktzcHboEh86BtF8+MSZIk9VenDDYZWI7OoWwIKYNtmu1vgyE3msEkSc3DQlILmgGMqvTAC8BqwFBSUmnLbt8JbJ/dvhpPh0mSJPVDpwz2GrAMlatDXc/2VTz7J0lSY7KQ1IKGUSGLvEgaUl0ppAyj4ydhKFxxJey1V/X6J0mS1IqGAWEqMBtYmj4Xh64AjGCSpEZnIanFTAeGl+94BZgJjCd3mBk2DP42BPYvuG+SJEmtqj2DDQdG0PcRRv+CYe+Fv2EGkyQ1NgtJLWQasDhdcstIYHTXnb0bcTG07VdQxyRJklrYNLIMNoVu1hfIYQRwGYzYM810kySpUVlIaiFj6FIvmgiM7V9bB42Cy4YOvE+SJEmtrj2DLcbA1joaBQcNhcsK6ZUkSdVhIamFLEe6CFu7sfQ7zJwN7HkucMAAOyVJktTiFslgA3A2sGdBbUmSVA0WklrIVCAW1NbiwPDF6LLgkiRJkrqaCsSJFJKbFi+mGUmSqsZCUqt6hrQ2Uj8dAlxZVF8kSZJa3RbA3IE3YwaTJDU6C0mt6p3ArP4/fSowp6i+SJIktbpXi2nGDCZJanQWklrVyxQ3z02SJEm9Ww+YPYDnHwJcW1BfJEmqEgtJLWYjYF0GlmEOBa4H+CLwzwI6JUmS1OI2AtZ9HmYv7H8bh74K1w9gRLkkSbXgBd5bzITs40KArYB7SJeizekzwCXATIBJwIwieydJktSaOmWwfvjMZ+CSW7IMJklSA7OQ1KK2BYY8AbcuyFdHOhh4mLRG90yAo3FotSRJUh9tuy0MGQK3AovdDCzR8/HtGex5mGkVSZLUBCwktahHgJuBUTsDbcA1wJKLHncQ6Qza43Q5A/YcML2qXZQkSWo5jzySZTCAnUnrBVTKYAfBhJPg8VNh5iPAkzXspCRJAxBibN4VmUMIzdv5GtiGjkWw/rUVLHkVsHTH4x8FLgfeBvgscBTwB9J0uCdIlw2RJKnOYoyh3n1QZ2awnvWawT4Kl18Ob68FPAu8VfMuSpLUq+4ymIWkQWI7YNh2wGXAcmnfPcAsgMOAi4E1geeBaXXooCRJ3bCQ1HjMYPm1Z7BhHfvuuQdmuai2JKnBWUhSsiMwssu+u3BRbUlSw7KQ1HjMYJIktT4LSZIkqSlZSGo8ZjBJklpfdxlsSKWdkiRJkiRJUlcWkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlUvVCUgihLYRwfwjhyuz+GiGEO0MIE0IIF4cQhmf7R2T3J2SPr17tvkmSJLUi85ckSaqWWoxIOgZ4rOz+qcBpMca1ganAodn+Q4Gp2f7TsuMkSZLUd+YvSZJUFVUtJIUQVgb2An6f3Q/ArsBfs0POA/bPbu+X3Sd7fLfseEmSJOVk/pIkSdVU7RFJvwC+DizM7i8NTIsxzs/uvwSMz26PB14EyB6fnh3fSQjh8BDCPSGEe6rYb0mSpGb1CwrOX2AGkyRJSdUKSSGEvYHJMcZ7i2w3xnhWjHGrGONWRbYrSZLU7KqVv8AMJkmSkqFVbHtHYN8QwvuBkcBY4JfAuBDC0Oys18rAxOz4icAqwEshhKHAEsDrVeyfJElSqzF/SZKkqqraiKQY4zdijCvHGFcHDgKujzF+HLgB+FB22CHAZdnty7P7ZI9fH2OM1eqfJElSqzF/SZKkaqvFVdu6Og74cghhAmkO/h+y/X8Als72fxk4vg59kyRJakXmL0mSVIjQzCedQgjN23lJkpRLjNGriDUYM5gkSa2vuwxWjxFJkiRJkiRJakIWkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi4WkiRJkiRJkpSLhSRJkiRJkiTlYiFJkiRJkiRJuVhIkiRJkiRJUi5D690BVdeupG/yDcAOwAjgZmB2PTslSZLU4sxgkqRWZSGpBewM3Ai8GwjZvpuBHYHLgftI3+g/AUsDHwUuA2bVuJ+SJEmtZGfMYJKkwSfEGOvdh34LITRv5wt0M3ACKci0Zft2A64ERgF7AhcC47LH7gG+D1yDQUaS1PhijKH3o1RLZrCkxwy2Lew5FC68E8bNT4+ZwSRJzaS7DGYhqQVsDfwK2JaOs2GL2Bx4BJgH+wOnA+8jnSV7Bni76r2UJKl/LCQ1HjNY0mMGmwYsASwLTEm79t8UTn8c3jfHDCZJanwWklrQ5tnHe4GtgLvoOBu2iMnAx4CbgHlp1wEbwKeehT/OTmfGDDKSpEZkIanxmMGSThlsMzpXk24GFiMtljSNdELvSeC9cMAE+BTwR8xgkqTGZSGpBS2khxFIG9L5mnw3AUsB7wFeJQWZB+CAA+EHD8MPSHP2DTKSpEZjIanxmMG6ZLCNgPuBYT086T3AecDhwHVwwJyUv8xgkqRGZSGpBXVbSFoXeBAY2cOTdwd+A3wAnns4XUHkYNKZNUmSGomFpMZjBssy2HrZjkfpYVh4BRsAj8NzmMEkSY2ruwzmVdua1JqVdq5BSjUPka4x25P/ZB9XhqOfgEfnwssF9k+SJKkVrVl+4390HgHeFyvD0ZPNYJKk5uOIpCa0KvB8152rkFZs7E9pcEtgIuw0DW6ZM7C+SZJUNEckNZ5BncFWye48Tw9rDPRiV+DXpKFI98FOwC0D754kSYXqLoP19xyK6qhTEWlFYCXS2Oj+ji9bFvgX6dIjkiRJquh5gBeybSDlzetJ09skSWpCTm1rMsuX31kOmMjAggzAVenDkqQ1IucNsDlJkqRWs3zvh/SbGUyS1EwckdQklgKWAV4p3/loce1PB84BtiuuSUmSpKbXnsGWApYuvn0zmCSp2TgiqQkMBZ4GxpXvXAJYG5hazGt8ELgHmFlMc5IkSU2vPYMF4PXqvIYZTJLUbByR1ARGAquX71gcmJZtBS0/ei3piiHvoW9Xr5UkSWpVi2Swoo2Ga4eYwSRJzcVCUoMbTprONq1856TqvNbuI+GLQ2CX6jQvSZLUNDplsDFVepH/AlvD7sAXMYNJkpqDU9sa3FzSlLZOiy8uA8yqwotdA3vtRrqSiCRJ0iDWnsFGA69V6UXmAAvTzb2q9BKSJBXNEUlNYHjXO9OLf415QJxHe5iRJEka7IYPBabQZaHKAu0A8+6GWKXmJUmqBkckNbjhlC2+GIAZpOHVBV8fdmfgtl2BW4ptV5IkqRkNB2bOJ61NWY2R4Jmdgduq17wkSYVzRFITiMCC0o0RwOziX+NW4J3FNytJktS0IrBgATC6eq9hBpMkNRsLSU1gLql+1G5YnToiSZI0iMwFRrSR1jIq2kKc0yZJakoWkhrcXNKlZzupxjpGBhlJkqR27RlsAWkxiKKz0mbA/QW3KUlSDbhGUhMYSgoz7UpBJhT3GnF74M7i2pMkSWp27Rkskk6/VuHEm+fyJEnNxhFJTWA+aTZbe9AYQuFhZntcZ1uSJKncfGBYG8QF1XsNM5gkqdlYSGoSC4C20p1qTG2TJEnSIhYsgLa23o/rsweBLarQriRJVWYhqYlE0my2UMB3bR3gsez2ljirTZIkqTvtGSxkNwYwKrw9gwUzmCSpOVlIakLto6vbSN/BfoxQehJYH9gU+B2w7e147VlJkqQeFDHD7Ulg/Q1h0yHwu/tg2wLalCSplkKMzbvEXwiheTs/AKXq33yy9bYX0K+S4EbA/4CwNcT7cMqcJKkhxRgLvLyEijDoM9jCbHRSf6wHGz2ZZTBcbFuS1Li6y2COSGpCC7Ot0wLcfbQ+KcCwOcR7sIgkqa4WA6bXuQ9/BT5U5z5IamztGWwYxKGky7oNJZ3d683q6dj1syISWESSVH9mMPWHI5Ka3DBg9rCyiuDMbGcv5m0A77gUnvgI8HDVuidJi1gZeDq7PRVYF5hC+lvsLWBsjfvzL2C37PVLfySW+0ZZn75Tu26pjCOSGo8ZLMtgZBlsKPA2PWewlYGJMA94B/BEtTsoSV2YwdRX3WUwC0ktYHj2cSYwdDjwJjCi5+dstA489jzE+Xg6TFJNrA08kt0uvUVFYG6X+7UKMjeQLrs9jJ6H55YPNPgxMIm0WO4x1euaurCQ1HjMYMkiGazca6S/0qZm9+emDxuRFtv2CyipVho+g10C7Js9eARwXrppBqs/p7a1sLnZtjgwby6MGwezRgOlrRRg1kj31x8Njz0NcR6mGElVtzHpj6yHSWGlvM4dKtxfLDt+JinQFO32rO2dstde5B/CP2QHfDbdLZ+98jXgZ9ltSVokg82FWaWdywOvlh1EWlrAIpKkWmnYDHYdjJgJQ0ovdjFwZdahszo6MbS0nWgGazSOSGoxY0i/d2OAV0hvBoyG1QNMeRuIaeS1XzhJ1bYNcB2pUDO6n20sJF2gsgj3kc5ijera5mXArmX3j8ruH0iag1LuNODbafdZwK3A+4BPF9RHVeaIpMZjBltUxQxGWhppSnbbDCapFho2g90KbZuwaBibTcdZu0rKCvLzzoazvmwGqxWntg0yrwCbkQILwAwMLpJq492k2kwbHX9I9VckzdYdN4A2HiGtCbAYXQLRf0jX3R5D5+Dydna/6zQVgDmksAPM/RPMPzoNyy69114PHDCAvqoyC0mNxwzWPTOYpHrpMYOtBjxIGim5IzAh2z8TGL9oW4VmsHugbW0qhLF+mAtzZ8H8i2DYkWawarOQNMiMA6bVuQ+SBpc9gfNJhZXFezpwJdIY6ynAeqTJ+M8Cs0hpo4tImqG7dD/79TywavmOm0hjvccysPHRc0jh61Lap8HNA/4N7DeAZrUoC0mNxwzWvXGYwSTVVq4MtgbwDGmoUXmFqBS0IAWZFTqeMuAM9iCsujKwBMUNbyop5bArgUPMYNViIUmSVFXDgQ8AF/V24BBgA+Ah4HXSpPxlSGmlNP9jPqnglIn0fVG/R4FlgaXIssvtwFrAkhQ7wX42acjB1cDBaQHJXXt+hvrIQlLjMYNJUuPoNoNtCNyY3W4jhaKeRNKwyiIy2EOw1IbQVnQBqas5wF8xg1WJi22rV5cCL2fbe4Gzgf2BX2T7nH8qqTt7A88Bv+npoDVIbyYvkcYfDyFVepbJHg/Z/WVJZ8Ne6vz0l/vYp12ABZQVkbbK2i56lcaRWbsfBC5Mu7YG/lnwy0hqXWYwSf1VMYNtSXrzuImObNVbEQlSFisqgy1VgyISpAW6l+i4awarDUckqd2ypCCzI2lQwGjSTJORpCVEpgLHkcKNJJX7EPB/3T24AWm0zlA6DZfuVQReJM3pz7wErJLz6fcD7yArJE0gjUaqttkw+3J488A0W2+7GrzkYOCIpMZjBivWssClF8GOO8KUT8Poa81gkvLplMF2JA1LGg4sN4BGi8hgL0FbhbWXqmIW8HeYfXCatWcGK44jktStK0i/bHcBW2T7ljkbRj8LS++TAgyk2SCnZMceU/tuSmo2m5PeMK4hrX3UlyISpLNiq5Dm82cqLKHUrRUpfjp+r0bCyGVTdtuMNG1fkrrTnsEugS32BVaGZc7PMtizMObZdMCSnzCDScphJCksDaSIBCmDdQldfcpgd0Pb8gPsQ1+MAg6AkX8xg9WKhSTxedKikKuTfgf5PXBQtuN3wBMd21JPwOqHw4nAl+vQV0mN5QDS28OvKj04gvQ+MpCzUSFr44m+Pe1OyhaGvJEuq21X2bbAhfA46f1VkrrzeWDaJbD63jBqVLZzedL7Xvk2Ns1KWR0zmKRkkQz2LuDcevWmw513wtKbUfxSAr0ZRVpp2wxWExaSBrlLgauAdUs7/gAcSMf1IlfMHizflkx/oH0D+FotOyup4VxHqj2v2PWBLYG/FPQiAVinb09Ziyy/3ADsQLqMSa2MBvaFdS+En9fwZSU1l0svhasehXX3IjuT14Nvklav/awZTFLSKYPtAvyZvg0b6k0AHun709ZaC4bWuohUYgarGQtJg9TlpHVv1yQt5j/6z8ADwIfpKCJ155h07DJfHPioSUnNbTowsevObUnz89eoeXfavZt0QTjWobZFpJLRMGFf+Eavl7CTNNi0Z7A1YcMNYfToHE9akRTYsqGWy2AGkwa7ThlscYoffR1I7zt91J7B6mU0TFgjFdxVPRaSBqkTSMMg1y7tWBvYlPQm1JsVs2NXgIOBr1ajg5Ka12KUvbkU6O58h11PGtk9tgpd6Is1R8NJ+1DhWrySBrMTgF9dAmv3533yKOCwdNMMJqnRNEQGewDW/BScVM8+DAIWkgapbwO7kkb/cQ6wUf/aWYE0T//rBfVLkioKwFb5Dt08O7QeA5HKjXkU9t8N/nxynTsiqaF8G9h145wjkbpaGVgp3TSDSWo0mwNb7QHD3qhjJ96CMY/D/qTZfqoOC0mD1GZkRSRI18fOMxKpq4Phos+lJZUuLahfkprLnqSFX6sukhaRzGEv0nDvulsFJh4Kv3io3h2R1Eg2uwBGr9b7cd36NHBYGuxoBpMGr5plsD7YC5h+DzCvTh24nzRykzTt7xd16sZgYCFpsDuLPi9i225V2H5N2AR4qsAuSWoeywPr1+rFvp3vsNuAfYH3AFMOAl6tXpd68vjL8Imz4d76vLykRrU1ZWfz+mE1YHXYHjOYNJhVPYMtBN7bt6fcBuz7N3jPkjClGn3qzTTgwXTVtk9gBqsmC0mD0Dm0j4pOi+IOYBLrah+E9Y4ceJ8k1d7P6fes1nbXkM6I1WT21s/g/TkPvZl0NZM5R1C3ifrTp8Ndd9XntSU1pk4ZbIBWA9YrqC1JtVV4BruD3Cfccguk+bOZ3BnsN3Dd2zCn4O70xXTACFZd9bown+poZ2BMUY2tAXttCI8BZxbVpqSqOJPsErGZnUgnxru7ssavSSGlJxOzrTR764SBdLA3R8K/r+rjc/4M7E7vl9Yu0MPAt0gnxSSp3M5nwZgVez8ur70wg0nNoOoZbDKckPOiJH3yu46b/877nOuAg4ALaL/SZE3cDw9/ywxWKxaSBrPTKOQykU8Adw68GUlVcBYwLru9B4suh/bOHp47nnRxoPOAJUjT3f+vm2NXJc3UOBk44WHSWbHv9avLFUXgI+f3/XlHXA3nzIJli+tKjx4EjgFuepwqV9UkNaX3kq5sWRAzmNS4ap7BHoATvkcxI5MWkhZhuzzLYH19/n/giEPgnHNg2RqFsAcnwzG3wU21eblBz0LSIHQs8Ftghd3peHcbgDVJi3ffN/CmJBXoTNLlofs7GGerbFsTuAz4Rw/Hvkxa8HU0wGTgxn6+aCUROBj++te+P/WfwKyjSWluuQL7VMH9wFfIAswU0jVwJanMscfCb38LK6xQTHtmMKkx5cpgJwJ/A/YmLXj2SyAbVdSvDPYK8BvSlLSBrMK9gLTAUFnu6kcE45//hFmHUd0MdjtwRpbBJllEqiULSYPQjgxsjceu1gO2Af5YYJuSBu6D9KGI9CXgBmAH0sqEZae4tyTVch4iTduq5IVsG0+6oNB6T8LxPyCNL+6nhcBnI/zhM/CZv/S/na9eBmcMheXOIF0vuwruBY57AG74ZbZjcnVeR1Jz2/EyGP0TBv5etDfwEKx3sRlMakSLZLAfkVbHLvd+Uu7aGFiZNPft+eyxc4FVYMvXIF7dhwz2Kqz3azj+GWB14KQ+dnw+8FlYeCF8FvgD8Jk+NlGuahnsZtKic8/CvTfBcaQYq9qxkDQIfZhs7dnvAD8jvXFJaik/p4fZEyeS0ka59wPvA9YFngUmZPsvAm6ErfaHb80HrkxnxrozkZR9VngVwulw0Ivwp5XhxH6cGYsL4ZzPwbBzU5v99X9A299g7FBgCfjWt2CVVQbQYJk7gIfvhqV/D9e9APR1DSdJg8qHfwhjizgzvwlpuMLFBbQlqVCdMtiPgCVJawYtUeHgPcpu71p2ey3gauCh9KteOi+XO4OdCwctA396GU5cEfgcacrbksApXZ44mzQvH9JZvPPSCcRzgGEUkME2g7Ej0/1vAQOKYNcBlwCPwx03p+La0tlu1ZaFpEHoO8BPgGXXBkYW0OC7YMcD4WMXwwAGDUgqyI9JeWGRX+8TScOLP0T6I6Sr0r/sawK7ZbfXIw2RfhW2WpACQKDnIdYArwC/mAz3nQXXLQ2TJ6dRzScuTQorPwN+UOGJ82D+l9MAqRhTf8/q5bXyuAjSH1zfgDdHwjKk11izH23dAPwd4DZ4/EJ4/ilY+z8FdFJSy/vOx+AnSxSwbtvVtL8R7wh8DDOY1AgWyWAH07+T9u8iFXWyUUxb/Ru+9e8+ZDDgvilw3VkweUlY7mk48XxShWtmlyfMoz1szSfLYNlDhWSwp2ifavcmsMxX4Ut/hzWfyQ44ksqXsPsbiyyVcMOD8Pf/ptuPkwZwrV1AH9V3FpIGofFk3/iDSX9NDdQMmPF691cdkFRbh9NNjfgTwBp0vmxIb3bOPt6TPmy1NHwzu2zHP3p56iukk0a8DmeckZZkmz6WNJn/L8DsNOz7BMqm8s+HBWfAGX3oYp+8DBedmm5+9ChYc3X43e/gXU/BnZ+CRzZOj32ENF2k5Kqr4Npr0+17yXLNw6Q/5oCnqtVfSS2lPYMN1D3ArenmDMxgUqPolMG+TeVRSHntTEcO2zKNTOpzBgPOmArjzofpAG/RHrIWyWCk5ZEKz2B/6rh5EcBU+OgNsObL6aJw75oId64Jj2THtGew64H702DvLIJ1ZLAyZrD6sJA0CE0jvUkU5n546FrwhLxUfyfRTRHpq6Sxv5/uR6M70xFkdugIMpGeh1h3NQ342ZvA2dmOn6W+vk0aBl4T53Xc/O1kuHJluPBCWPk5uOw5+Me66bENSPWuu7Jjb74Zbr21Vp2U1KqmUXAGI62dYgaT6q9TBjsB+BrFXaVxB+BbA8xgXfbVPIOVXJAu/HQlcCGw8qWdFxRfJIPRXjdXA7GQNAgNI1Wav0Axl8S+k/RGIKn+vgqM6LrzS8A3KeQqjaUgswBY8O+BNzebOgSYzJ//3HH7L8Djf+u4fzkpxNxT4z5Jam1mMKl1tWewr5MKSWMKfoEdgF1T/iqiIF3XDFZ2+y+kaWolZrDmYCFpEDqWNOU09nJcXndhiJEa1hdJp8gGMrS6qx1gwbdgfqRlFpfuurbI5XXphaRWdyxmMKnlHUfxRaQyC0hrGbUKM1hzGlLvDqj2ziZdxrGIi4ZIanDfoNgiUmaHHWD//YtvV5Ja2dlnw2emDiCD3Qz8FLiluD5JKtBRFHMxo+5sBTvsAvtX8SWkPCwkDUInA1MAzi/d6J977oHTb4GbiumWpAE6kjRtot2hVPWM2MZ0XNxNktS7k0+GKQPIXkwDnod7psHpmMGkRtGewX4AjK7iC+0MnAQb72IGU31ZSBqEPko2QOFk0gTUfnrtNXjkapjk6mdSQ/gZXdZH+hEwtnqvt+MG8PHDgV2r9xqS1Eo++lFYYiCjRPcFTofX9k5XOJpUUL8kDcwiGaxaHgSGwo4Hwsdr8XpSNywkDUInkw2p/gCwZP/bGb8nbLcDrP5cId2SNEB/BubV8gWXh3U+Ae/6TC1fVJKa18mbw3L/JI0s6q/7Yfx9sB2weiG9kjRQNctgk4Dn0s11gHfV4jWlCiwkDWbbM6BLUj4LnL4a/GW9ojokaSBuoPjLSvfoJVjx0TTFTZKUw9dJC1X2dyjR/cAP4Nm/paltXReplVQfNctgewBbAXfBipjBVD8WkgahfwEzgKu+DtP+CEztXzvrPgYb/wj4U3F9k+plJ2AFUn31Q9m2OrAFsHb9utUn+1DbS3E+vRv8bHf4zZ01fFFJamLtGew/MG1aPxq4GPg7rIt/QKp1mMH66Gp4+o9pOt1vavWaUhchxqIuQFp7IYTm7XydPQn8AvgasPo9wJZ9e/7jwA9PhvO/WXTPpNrbETgDuAJ4L2m6AKSzvWsDE4FTgKdJ4X0o8L/adzOXGZQNNJwMLFud13mGdOGgM88hnV2XqijGGOrdB3VmBuu/9gz2IKy+Sc4nPUQaCv4nePzv8EPSNVOkZtcpg20N260I3AenLwNrrwoTH4FTnmmiDDaNqlwtt9wzp8NPj4Yzq/syEtB9BqvlyWs1mF+XbtwGrAWMy/nEJ+CiZ+D8J6vRq9rakbRY5fR6d0R1sx3pbM4mwGblD2wCX1wJeBRYDMJs+OuzabH6ScDxte5oTleTzooN6+3AAfonBhhJ6q9fA9wCrEa+PzrPoj24XURrFJHMYGrPYFvDZksDJwI7AD+BL24KvA/4LYTL4K9kGexZOP6JunW5R1cD+1wDw/ajqkHMDKZGYCFpkLoNeJ70Br7Y0cDWdAzD6MWEc+CZU6vWtZrYnnRlzt8CpwFPAPcCawAvAK/Xr2uqsZ+RikjtNiGN4jkeeA/wS2AVOHQKHHpJOuRiYLcX4KWn0s9OI/kgaQ3XKp8MkyT1U3sG+wIsFoCP0fOb9mPAi+nmBNKI0GZmBlPJz7aCTZYATqXz7Iivld0+Eg49Eg7N7l58Mex2YgNnsA/DEq8DS1XpRZ6n8T5xDUoWkgapT5HWabuIbBrMA8AG9P7X59Pw+xfTlQma0TakS3OeD6ya7SuNzDoa+ATpjMcfSFe2e500O0itZUtSiAUYW9r5DtKovB+SJuuXHFN2+/D04UDgwP+DS0+A0yakn5FG+Td9R8re2G8nLco4vEovNonG+cQlqUl8irIM9nkgkIZalGew/wFTsts/Aa5MN39Pk2ewreH8kWUZ7H7gLTPYYNIpg/2SNAKpDw48EA4cApeeDqe9ApMbqKDUnsFuI42mKjp/vUAqujkcSY0gxti0GxDd+r/9BeJbEB+E+CbEeCYx3kaMU7t8qZ/K9t9GfPYg4icaoO993baEuB3ESRBjaduUGLcr28Z2PHYcxGshfgfisg3Qf7dit/NJP/uxfLu6H29Df00/O/9ai7hOA3xe20Cc2fXzern4t9+Jkfill4l8vf6fs9vg2OqdN9wW3er9M9Hs2yIZ7OayL+9jxLgznd/LIT5LE2ewLYiTtiPGl7r8KH2ajhy2RPo8zWCtvbVnsI2J8cEBvhVdRvzXdsR11qn/57VIBptc8NvuC8R4JHEixC81wPfRbfBs3eaAegcRQ0x9twshrg/xHIh3Z9u0f5d9mZ8mPrNXx2PHQVyhAfrd1+0Vyt7YNybGLYlxQpcfqcOz/Ut0Dm6faYD+uxW7bQjx8fKfiXWJ8Y4BvB1dSfzPlsS11qrv53ULxPl0/vmNVxHjnGLedidG4t2ReMwkIl+t//fRbfBs9c4bbotu9f6ZaIWtUwY7m3j33cRpdxPjTh3v4c/Q5BlsU+IrW5Lvj+rPYwYbBFt7BruruLek//ynATNYwYWkiSel94FjGuB76Da4tm5zQL2DiCGmcba1IG4C8ZZfE595gPjAA8QH9iDukz0+HuIpED/dAH3Nu22UfU5TIMYNiHETYnyxlx+tL2THZWHm2xCXboDPxa247WE6B9V4QzFvS9deS9xkE+Iaa9Tm89iky/2NqTAiiRw/8zm3kyLpv5/U/3voNri2eucNt0W3ev9MtNrWnsFIxaMHsq2pM9hGxE2mEqf05Ufr/cRIx2YGa73tYYhxHWJ8uNi3pWuvJW6yOnGNGn0e3WawjYixjWILSS8TTzqy/t87t8G5dfeD6RpJTWiDCvseK9v/OOm73lcnka1z94V0Ne8bgfHAW9njnyJN2f9TP9quh/VIU5THQrpe6I2kSfe9OSP7eCBwCXyXNEf/t8V3UXWwFjCyfMfqwJhi2t5tN3jwQbj1P/DxPdJ6iNV0AbAp6Wcd4G66fG7rkq6XW6QpuGiFCjGEjp/dSFrjYv3s/mN16ZHUu6pnMFokg60Ht90GY8f2fmwnq5K+mBOBN81graY9g10CbFxs27vtBg+eDrfuAx+nThlsPeBc4PNAWwEvMgmYSlqV/vcFtCdlCslg9T6j5dmwvm3r0PlMTWlbl47hlBtCDAW93s8hfqABPu/+bG+Wvj5rEONr/fgRO5oY1ybGxYknQlyiAT4nt4Fv/4E4m7Lfn1ur8PZ0M/H68dX9PNYu2+aWfz4Q45rEOIQY7yTGrYlx0sA+n8mR+ORrxC8+SeTL9f8eujXnNhTi6tntNoibQ5xA+pl9G+IGEJ+COI/OP9/Q/dkwt/pt9f55qsdmBsu5rUl8c8YAf8SOJsbF09fUDNY6W3sGu79Kb03XEeNKxOur/Hl0ymBrkf5eCMQ4gxg3I8ZpBXwurxLjx4mTIT4J8YsN8P1za96txwwWiBuslWWwtYlrr5Uvg9U9iBhi8m2rkYZqxpzbetlzigozzbi9CTGuSoxvDPBH7TPpa/odiGMb4PNyG/jWaWpbNQpJkXjb9el3cMUqfQ4Lyj+H0rYqKcg8RFoLbLUsiAzg85gSiYe+TsQh1W4D3NYgLRK6GsQtIE6F+A46fn5nk0JLXC3btxpxIcTh0G2IcavfVu+fp1puq2EGy72tSmQm8c0ifsw+3vE1/Q5msFbZHobqFZIiMd5EvG35GmSw1bJtbva66xHjWwV+HoelpTkObYDvmVvzb50y2BDi1NWI7yj9DK9PnD2buPba2c/e/PS+u3A14vDhxO5+SJ3a1gSGAs/18TmPZx83IA1VWwF4FVhYXLca1vjsY1gJeIjOl9PtjyWBMXDMzDTC9HcDbE6Dw+Yj4K5l4dnXYLsC2x3f087Hgc2APYFHgHH9f503gLeBbwHn/QDnFWhAhpL+HVqJzv+ePRSyncAI4KlhwFOk+Q/PQFgV3jGxlj2VOjOD9c34J4CREApu9xjMYMrpXbD53+Gu/aucwZ6m8/S1xxm414FZWQabmWWwAprV4DZ0KKywfFkGWwZ4IP2ZXDICeOqpsh2rQHgO3rFVD+0W3VEVb4UBPPcx0nzHW0lTkV8ppEeNbQLZ/OvHgcULaPCnwEw46bcGGOV37w5w6pXwzfeTgkEBVgBepEJAfwYYnt1+YmCvMRWYDRwzDf5vVrZz5sDa1ODWBmxNWrOOIcDSpN+J5Ujv0U9WeFJpcYsX6fjZluogVwZbno435rKgNagzWMFOwgymnObAvevBqX+Hb+5UXLPtGaz0plB0tRSYegTM/lsqnP5f8c1rEGoDtl4DbquUtXp60gu9H2YhqcENA+7ryxOWyT5O6dhVKpAvRVqwsL9nxMYAo8ruv00aKzeE/v+dOYQ0aOKNfj6/lhYjff6zejtQAnYELt8c7rwE+CAwbeBtTuq6o/T7XmCYOQz4G8CxeBpMAzKU9P6+HnBLaedKpBVJdyHfao6BtJqpVAc9ZrBlym6/TAo0kE7rzu14aFBlsGV6P6S/zGAtZiowj/RLVrTbYMefwuWfgzsLbLY9g71MsUWkaaSvBXDYnCyDSQPUnsHa4JYlq/MaQ3o/RPU0D9iot4OWLNteBl4DVux8yDTgvwxsdNMpWdOl7QfAV4AfDaDNFUh/YATSD3vpYyW9/Q50ery7Rgbgh8AXim9WdfAmML8WL3QvDP0+jL20Su2/QvplLCCIzQDemAFz3yD9VTFn4G1qcNuYNDqhvYhU+qt1Bbwsm5pCjxnsCVJl6DU6p+lSHiszjQIz2Ch4bUn4wZLwlSXhR6N6e2b3Cs9gE6nOcCTMYK3kTWD+rsADVXqBXYDjgY+lP6b7euHAblX43S7EB2HGcvDGcjD3yiq0r0Fp4zaYsCTcsjXFVlTLWEhqcMuQ40rbU0h/+L1Bxx+UL5OKSdm75zrAalQY0dCLUXT/N+qXhsP3x8LwsenyriP7GB4C6bK22wFrkKYajycNGx3bZRtC+jR7+sfgKdKMCQBeophpbWpJ7yNdRhlIP4TVqiptB5ueAlfuMrBmSr8H3e8YuA/PhKU/AlcsTfpFuqjY9jX4PAS8H9Ib+FjSX+QP17NHUt90m8HGUhY4uniFtAjF6I5dhWSwkdnr/hB4A770Bnz/DRh+WpbBxtJ5yFIvCs9gY7v/kkjlOmWwatkJuBo2HQMDrc2MJf2Otf+tVYUpbR8m/f5cUXzTGqQe2gze/wZwe/Vew0JSAxsJPNvbQWNIY5pjhcdezrZx6bCnSLMKRpPeA0fR8QNQfrvcr4GDs5dpDzPDsx3fAKbDEdNh+nQ46ywYMwZGjOj9cwvAhunpTCcFmKVIAWaxsv3TSUtpjMn6Ny27Xcky2ac7cgxVeZNX67gBWBNYALA7ZUMmivdAG+zedU5CHwQ6fheA9AswLdvRVvk5/THy49B2TXHtSZuT1oZhPdLP60M9Hi41lB4z2ES6DyOQhvq83XG3kAz2E9Lv0Zc6H3PEEVkGmw5jzulHBgvw9BhYavQAM9jz8PKwqg1IUgvplMGqaWt44MoU8/qrPYNN7+3IfpoFLEi/NwVGOqkjg1WRhaQGNpscZ3dmkAJLd5PuVwKmpZNjK5DCyv9IJ4bvA7YgvUleD+xAqhGV/1AcAexKOmt1RGnnd7Md3+n8Up/4BLz1FvzlLzC8wuKopeA0ijTl4ZFK/a2QQIaShsGW2phB57/JR5VtQDptVnYmsChzqdF0KFXdNrRfsKDqttgKrnsTuLmgBmdQlULpPy6G9w8kbUlSC+kxgy1On96Ln6ObDDYCwii4PsAOw2H4kAoZbBi8NQqO6GVV008Abx3Yjwy2AinTvVh2QCmLjYKhI7MMNqqXDLYk7eu8VIMZrHW0Z7DZVLeadDtssQtcV8WXGLB9gP/CP8hG8EpNxEJSs5tDCgA5y9hPkYZXP0y6ksjdpDf0IcNgWICrgPcPhSHZT8Z5pLNhQKroDO/9tfY7AC74K4ss5b4F6QTd23RzYnop0pDw4V02Ot8Ow1MbpZNub5W1OxzS16TSCK0B+hrwi+KbVZ3cTKrWV9u9wM4PwPCt089nX5YzquXFquZ9CBZeVcMXVEsLeLE1tbCefrh7yCCdMtgwuPtu2OZtGLIjDLserto5/TFZCufnAQefRQo4R+br2n4HwAWXLLo/VwYrHfA6aT7P23RUjt4Chqff7YoZbHjBv/Pz6HSS1AzWWm4GNn8XafGw/q5A35shcO8w2JmOPyP6lcFyjPLrt2uBnRf5cZcGJITKJxSKZiGpCcylh5M8Y+m+aJLzzNAdQ2Gb/8L174ZdhsEVV8K+e6U3204n204jBaSv9dzeP4CP7QNtV9BedArkvETgEtlrlLYZpDNjc0jDP4eRThNmHyr+g7AMnYaUD9h8fHdvQduS1p5bCFX9Hm9N5x/px6j8u9B131DSNIJOP+PVuLpJ5gNXwNV7OStUxdgGuKl0p4o/t1K1VcxgM7KdlfLXWNpzSidDSb8Lpe0h4B1wB1kG2xF2aYMrhsG+w2DYsPTHQF/9A/hYgLayk36dMtiwLlulf5BG0zGfeigpUw2hU/CaPSy7Wfr9zgpNhTkYuLDA9tRQ2jPYLhQ3YrurHWDrW2DO0AFksEDl3+eCfQC4GjOYirHNNnDTTb0fN1AWkhrcXFIhfLHuDnib7t91lqRjTli5ti7bnaR39BuyF9wdLh0Cc9vgoNIxfXhn+2DWzBV7kIY4kabN3dFbPyqNdBpOx7Veh9BxSd3s2LnAkLbO+5hBz2sX9NURwFkFtqeGsQNp/nB8L+lf8GoVDAPt77Zrka5kVf5jvzRpQdfyX4FIunpOp4unzaVqKeNK0q/a/tVpXoPM3cD2AO+gh+unS42t2ww2ItsWkE5EzCe9aXc396qNNHVsbtm2foXjrk6PXToX5s6Fgz5Gn5P6B4G5+8IVZav2tmewTbr0YS7wQh8aD6R/lNrSx7nAkLfK+lj6Okg5tGewar7INsBdtAeutdpgQlsqtLaRPi7dBpOzfaXj2jNYrrPgA2cGU5HuvhO23776r2MhqQmUTgYtoMJU4mFlDyyg49240pzj0h+zE+gIPvNJ4527+keXYw6jX3/ADgnw7lBhLeN1urQ/nxyXp8uE7PghdBSXhsCCmRBLZ8bKvxYFWFhsc2og7yY7GbYn8O8qvcjmwIO0/8yuNgTml7aQpjusMQSmD0nTSkdmj3UqIlX53XohKcBcWt2X0SCwE2UnmB8mFZOkJtVrBitt80hTwOZ2PYh0ubYV+vHi5wOf7MfzMkNI/8bdUr5joEoZrHTiZUS6vWAhxNGkL9ZAA1OX0GUGa13vBm6OVPcbvDmd/t5YbT7Mfx7mrwjz58NT82GN+TB9PgyZDyNnpkPnDKPy73PRghlMxemUwarMQlITmE8KMqUliqDLwIlhZQfMzR5cnHQ1t3IvZY2t3o9O/A74Qt+esiewYDe48cYKDxYxqqJLqhs1EubMg7gQFo6FhW8V8O9ShBjh88DpA21LDekW4F3UIKRuTOeib2mbDEsuCdMWwJgFsGABzFpAGkpd+j0ZQtUvb/J+Bn6JXAnSkhfvzG77x5+aXa8ZrGQElUck1Wmuyp6kfzZuLO3YCri/wBfoMmW1PYMtBgtnwsLYz9//CBwA/F/HXTNY67oFeNeuEK+htv9gjCdd6pk0gWMaaTLDArKJELUqIgFcmy52YgZTETplsCr/TllIaiKlk0CRNPKy4s9G6fqRsyo8Np50ydp62xB4osD2hgIL0/v9WGD2GFj8bWiLMH0gv0AR+AR89pxUR1NrisAupPUO23fU0jLAGxX2l0JMDYpIUjXcC2yGxSS1hlwZrJIpwLJV6lRe25PmmxZpIZ2KZO0ZDFh8cWibkWWwvrwBRGBv4LKOXZ/FDNbK2jPY7sC/aIx/MEZQuyKSVCX33gubbVbdYpKFpCYS6ZjN1S8vAysX1p3+2QR4tLovMZo0Sy4uAeOm93Z0Dw4CLiimT2pc76RsUeC9SdM6G8VQaltEupy0wIZUkIdIlzqXml2/Mtg00tXQ6mkn4LYqtd2lmARdMtgQKp8o6c77SMUEDRoNncGkJvfQQ7BRFUNYiNUe81RFIYTXSBO4ptS7L+pkGfyeNBq/J43J70vj8XvSeFaLMdZ7TIe6CCHMoNjxxSqG72GNx+9JY/L70nj8njSebjNYUxeSAEII98QYt6p3P9TB70nj8XvSmPy+NB6/J1I+/q40Jr8vjcfvSWPy+9J4/J40F6e2SZIkSZIkKRcLSZIkSZIkScqlFQpJZ9W7A1qE35PG4/ekMfl9aTx+T6R8/F1pTH5fGo/fk8bk96Xx+D1pIk2/RpIkSZIkSZJqoxVGJEmSJEmSJKkGmraQFELYI4TwRAhhQgjh+Hr3ZzAJIfwxhDA5hPBI2b6lQgjXhBCeyj4ume0PIYRfZd+nh0IIW9Sv560rhLBKCOGGEML/QgiPhhCOyfb7famTEMLIEMJdIYQHs+/Jd7P9a4QQ7sy+9heHEIZn+0dk9ydkj69e10+ghYUQ2kII94cQrszu+z2R+sAMVh/mr8ZkBms8ZrDGZQZrHU1ZSAohtAG/BvYENgQ+GkLYsL69GlTOBfbosu944LoY4zrAddl9SN+jdbLtcODMGvVxsJkPfCXGuCGwHfCF7HfC70v9zAF2jTFuCmwG7BFC2A44FTgtxrg2MBU4NDv+UGBqtv+07DhVxzHAY2X3/Z5IOZnB6upczF+NyAzWeMxgjcsM1iKaspAEbANMiDE+E2OcC1wE7FfnPg0aMcabgTe67N4POC+7fR6wf9n+P8XkDmBcCGHFmnR0EIkxToox3pfdnkF6gx6P35e6yb62b2V3h2VbBHYF/prt7/o9+f/27ufFrruM4/j708Q0gpCmiWSRiabBgFSQZKGoFQkJgk1D0kWQiKWp9B9wIUVdKAS6cGMj1IVgCxK0NURrZltICt0YGm0woAgx9EeCbaD50YVaSfu4ON+Jl1DxTrm958y97xcMc74/7pnvzMNwHp7zaylWJ4A9STKd1c6PJAvAA8DPWzsYE2k5zMF6Yv41TOZgw2MONkzmYLNlpRaSNgOvj7QvtT71Z1NV/b1tvwFsatvGasrapZ87gTMYl161y3fPAVeA54G/Ader6mabMvp3vxWTNn4D2DDVBc+Ho8BjwHutvQFjIi2Hx49h8Tg/IOZgw2EONkhHMQebGSu1kKQBq+5VgL4OsAdJPgb8Bvh2Vb09OmZcpq+q3q2qHcAC3Vn8T/e7ovmWZB9wpar+0PdaJGnSPM73yxxsWMzBhsUcbPas1ELSZWDLSHuh9ak/by5dltu+X2n9xmpKknyELoH5ZVX9tnUblwGoquvAaeCLdJewr25Do3/3WzFp4+uAt6a70pl3H7A/ySt0t+PsBn6CMZGWw+PHsHicHwBzsOEyBxsMc7AZs1ILSS8B29tT3tcAh4DFntc07xaBw237MHBypP/h9oaKLwA3Ri7z1YS0e4afAv5SVT8eGTIuPUny8SR3te2PAl+le27CaeBgm3Z7TJZidRA41c5gakKq6ntVtVBVW+mOG6eq6psYE2k5zMGGxeN8z8zBhsccbHjMwWZPVmo8kuylu89yFfB0VT3e74rmR5JngF3ARuBN4IfA74DjwCeAV4GvV9XVdnB9ku4tI/8AvlVVZ3tY9kxL8mXgReA8/73v+Pt09+gblx4k+SzdQwJX0RXtj1fVkSTb6M7E3A28DDxUVe8kWQsco3u2wlXgUFVd7Gf1sy/JLuA7VbXPmEjLYw7WD/OvYTIHGx5zsGEzB5sNK7aQJEmSJEmSpOlaqbe2SZIkSZIkacosJEmSJEmSJGksFpIkSZIkSZI0FgtJkiRJkiRJGouFJEmSJEmSJI3FQpKkqUrytSR/TXIhyXffZ/wrSf6Y5GaSg7eNvZvkXPtaHOm/J8mZts9fJ1nT+u9s7QttfOuH/gtKkiQNkDmYpEmxkCRpapKsAn4K3A/cC3wjyb23TXsNeAT41fvs4p9VtaN97R/p/xHwRFV9CrgGPNr6HwWutf4n2jxJkqS5Yg4maZIsJEmaps8DF6rqYlX9G3gWODA6oapeqao/Ae+Ns8MkAXYDJ1rXL4AH2/aB1qaN72nzJUmS5ok5mKSJsZAkaZo2A6+PtC8Bm5McSbL/f3xm1NokZ5P8PsmDrW8DcL2qbo7u8/af18ZvtPmSJEnzxBxM0sSs7nsBklRVPxhz6ier6nKSbcCpJOfpEhNJkiQtkzmYpA/CK5IkTdNlYMtIe6H1jaWqLrfvF4EXgJ3AW8BdSZYK46P7vPXz2vi6Nl+SJGmemINJmhgLSZKm6SVge3vDxxrgELD4fz4DQJL1Se5s2xuB+4A/V1UBp4Glt4scBk627cXWpo2favMlSZLmiTmYpImJ/8+SpinJXuAosAp4uqoeT3IEOFtVi0k+BzwHrAf+BbxRVZ9J8iXgZ3QPgLwDOFpVT7V9bqN7aOTdwMvAQ1X1TpK1wDG6s2ZXgUPtTJokSdJcMQeTNCkWkiRJkiRJkjQWb22TJEmSJEnSWCwkSZIkSZIkaSwWkiRJkiRJkjQWC0mSJEmSJEkai4UkSZIkSZIkjcVCkiRJkiRJksZiIUmSJEmSJEljsZAkSZIkSZKksfwHMP3SuU1LmUAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1440x4320 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "c='9'; fov = '202'\n",
    "ref = imread(f'tif/Cycle_0/Cycle_F{fov}.tif')\n",
    "mov = imread(f'tif/Cycle_{c}/Cycle_{c}_F{fov}.tif')\n",
    "\n",
    "ref_max = ref.max(0)   # max proj of z for each channel in reference\n",
    "ref_binary = ref_max[0] > threshold_otsu(ref_max[0])  # binarizing for channels - T/F\n",
    "\n",
    "mov_max = mov.max(0)  # max projection by each cycle in for loop (ref above)\n",
    "mov_binary = mov_max[0] > threshold_otsu(mov_max[0]) # binary of moved image\n",
    "\n",
    "out = imread(f'reg_bin_Cyc_{c}/Cycle_{c}_F{fov}_bin_reg.tif')\n",
    "\n",
    "fig, axs = plt.subplots(1,2,figsize = (20,60))\n",
    "axs = axs.ravel()\n",
    "\n",
    "im_reg = np.zeros((2048, 2048,3)) # empty color image\n",
    "im_reg[...,0] = ref_binary\n",
    "im_reg[...,1] = out.astype('bool')\n",
    "    \n",
    "im_orig = np.zeros((2048, 2048,3))\n",
    "im_orig[...,0] = ref_binary\n",
    "im_orig[...,1] = mov_binary\n",
    "    \n",
    "axs[0].imshow(im_orig[:500,:500])   #before reg\n",
    "    \n",
    "axs[1].imshow(im_reg[:500,:500])   #after reg\n",
    "\n",
    "axs[0].title.set_text('Before Registration')\n",
    "axs[1].title.set_text('After Registration')\n",
    "for ax in axs.flat:\n",
    "    ax.set(xlabel='0:1500', ylabel='0:1500')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bbe6c7fc",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}