Skip to content

Instantly share code, notes, and snippets.

@votti
Created June 3, 2019 10:55
Show Gist options
  • Save votti/da664fcf1bfa189adfd08a383b6fb410 to your computer and use it in GitHub Desktop.
Save votti/da664fcf1bfa189adfd08a383b6fb410 to your computer and use it in GitHub Desktop.
IMC ome cropping example
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from imctools.io import ometiffparser\n",
"from imctools.io import tiffwriter\n",
"\n",
"import matplotlib.pyplot as plt\n",
"\n",
"import numpy as np\n",
"import copy\n",
"import os\n",
"import re"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## An example how to crop an imc ome.tiff"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"# a path to an ome\n",
"fn_ome = '/home/vitoz/Data/Analysis/201805_cp3_processing/ometiff/20170906_FluidigmONfinal_SE/20170906_FluidigmONfinal_SE_s0_p1_r1_a1_ac.ome.tiff'\n",
"# the output folder\n",
"fol_out = '/home/vitoz/temp'"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Load the file"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"ome = ometiffparser.OmetiffParser(fn_ome)\n",
"imcac = ome.get_imc_acquisition()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"def crop_imc(imcac, x, y, w, h):\n",
" \"\"\"\n",
" Crops an imc acquisition object.\n",
" Args:\n",
" imcac: an imc acquisition object\n",
" x: x-coordinate\n",
" y: y-coordinate\n",
" w: crop width\n",
" h: crop height\n",
" Returns:\n",
" a cropped imc acquisition\n",
" \"\"\"\n",
" out_ac = copy.copy(imcac)\n",
" out_ac.image_ID = get_crop_fn(imcac.image_ID, x, y, w, h)\n",
" out_ac._data = imcac._data[:, x:(x+w),y:(y+h)]\n",
" # out of some legacy reason the axes need to be swapped to xyc\n",
" return out_ac\n",
"\n",
"def get_crop_fn(fn, x, y, w, h):\n",
" \"\"\"\n",
" Creates an output filename containing the crop coordinates\n",
" \"\"\"\n",
" re_fn= '(?P<base>.*)(?P<meta>_s[0-9]+_p[0-9]+_r[0-9]+_a[0-9]+_ac.*)'\n",
" base, meta = re.match(re_fn, os.path.basename(fn)).groups()\n",
" new_base = '{base}_x{x}_y{y}_w{w}_h{h}'.format(base=base, x=x, y=y, w=w, h=h)\n",
" return new_base+meta\n",
"\n",
" "
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"imc_cropped = crop_imc(imcac, x=0, y=0, w=40, h=40)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEYCAYAAADPvfYMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXecJNlV5/s9cSNNeV/V5dpNm+keq5nRzGgkIQmBECA9SYhlMSvg7YLwb4E1sDxYzOKWXWDfLh6E25VgWSFYJEBCCwjkxpue6ZnpaVfdXV1d3maZzIx77/vj3IzMqu6e7vHTNfH7fOpTWZGRN0xGxYlz7u/3O+K9J0OGDBkyZHitInqldyBDhgwZMmR4JZEFwgwZMmTI8JpGFggzZMiQIcNrGlkgzJAhQ4YMr2lkgTBDhgwZMrymkQXCDBkyZMjwmkYWCDOkEJGfEJH/cZXr/oiI/M6zvD8mIl92Neu+GiCK3xORBRG5X0TeLCLHXqSxPyMi3/ZijJUhwwuFiBwVkbeG15uu+7Dsu0RkSkRKItLziu7sy4TXTCAUkYKIfEhEzojIiog8IiJf2fD+20XkaRFZE5G/F5FdDe99nYh8Ibz3mS3jvjlcMI0/XkTe37DOD4jIpIgsicjvikih4b17wo13RUSOiMibtoz/jWGfV0Xkz0Wke8v7Xy8iT4X3T4rIm69wHt4qIm7L/n78uZ5P7/3Peu+v6ub+XNZ9PhCRd4VzuCoicyLyYREZaXj/W8N38m+2fG68dkMA3gR8OTDivb/Te/9Z7/3Bl2qfL3EMeRH5xbBPJRE5LSK/3PD+mIisb/nefuUFbG93uM7XwnX/ZS/CMXyviDwoImUR+f0XME5eRD4ajtk3fEcvdP9+S0SOhev/W1+MMV/tEJHfF5Gfblzmvb/Be/+Z8Oem615EcsAvAe/w3rd67+de3j1+ZfCaCYRADJwD3gJ0AD8G/Em4IfQCHwvLuoEHgf/Z8Nl54L8AP7910HDDbK39AO8CSsAnAUTkK4AfBt4O7Ab2Aj8Z3usG/gL4T0An8AvAx0WkK7x/A/CbwAeAAWAN+LXatkXky4H/CPzfQBvwJcCpqzgXE4377L1/91V85lUJEfla4CPA/wf0AjcAZeBztfMYMA/8kIi0X2aoXcCY9371pdzfZ8G/A+4A7kS/y7cBj2xZ591bvrfvfQHb+6Mwfg/w/wIfFZG+FzAewATw08DvPt8BRCQOLz8H/DNg8gXuUyMeA74bePhFHPNax9brfgAoAkdfuV16BeC9f83+AEeA9wMfBL7QsLwFWAeu37L+twGfucKYvwf8XsPfHwF+tuHvtwOT4fW7gKNbPv8M8C/C658FPtLw3nVABWgLf3+htu5zOOa3AuOXee8ngP9xufWAMeDLtq4b/v4AcAaYQ2+sz7buNzes+2OXWPd/Af8DWAEeBw6ggWIafZh5R1hXwjj/dst+RsATwE+Fv78VvbF+HPjxhvXGw3H+C2ADsOhDzE9uPf6wj/86XDNL6INSMbzXBXwCmAEWwuuRhs9+Bvi2K3wvnwC+/1neT8/Rc/iu9wH/EPZ3FvifYfkB9GGhrWHdzwLfeYXxfgL4aDj2FTSg3HKJ9X4a+P3ncj0CP4QGvf++5f1x4K1XOdZXo8F9OVwnP3GZ9T4HfOtVjnkn8EVgEbgA/AqQb3j/BuDT6IPWFPAjlxlnCH3onQdOAN++5bxe1TUf1u8APhT253w43+YS2/wgUEXvGSXg443X0iWu+z8CVgEf/v6753K9Xcs/r6WMcBNEZAC92I6iF/Njtfe8Ph2dDMufy5jNwNcCf9CweNPY4fVAqL1L+Nk0DHDjpT7rvT+JXtQHRMSgGUSfiJwIJbVfEZGm57LPLwZE5DDw62gwHEKzjJFnWffXgG8CBtF/6uEtq70b+O9ogHkE+BQa3IaBn0KzZICDwE70JpLCe++AP0VLPo34MeAHtpaXvfcfAr4T+KLXTOvHL3OoXwe8E9gD3IwGWMK+/R76dL0TfYh6rmXLe4EfFJHvFpGbRGTrdfF88B+Av0HP4wjw38LyG4BT3vuVhnUf4+qu9/eg57sbfcj781BOeyHYEcbbhd68ny9W0YesTjQofpeIvPcF7psFfgCtNrwBfZD9bgARaQP+D1r9GUIfPP72MuP8ERrUh9B7xM+KyNsb3r/aax70/pKE7b0OeAf6kL4J3vvfAj4M/IK/ROXnEtf9N1C/Bjq991/6bCdmO+E1GQjDP+6HgT/w3j8NtKJPzY1YQktUzwXvR5+8/6Fh2daxa6/b0IxuSES+QURyIvItaNbXfJnPNu7XAJBD/6neDNyK/lP86FXs55CILDb8fN3VHuBl8LXAJ7z3/+i9L6MBxz3Luh/33n/Oe18B/j36BNqIz3rvP+W9T9Cbbh/w8977KvDHwG4R6URvTqBPxltxoeF9ALz3j6KB4Yee8xEq/qv3fsJ7P49ml7eGcee893/qvV8LweVn0BL8c8HPoWXub0JL8+fD9dCIP9/yvX37FcasosFlyHu/4b3/XFj+Qq73h7z3Hw3fxS+hZbS7r+JzzwaHZupl7/368x3Ee/8Z7/3j3nvnvT+CBp/n+j1sHfMh7/293vvEez+GBqTamO9Cqzu/GM7vivf+vq1jiMgoOhf3Q2G9R4HfQR8ca7iqaz48wH8lWj1Y9d5PA78MfP0LOc7XOl5zgVBEIvTJqwLU5lhKwNa5o3a0TPFc8C3AH/pQf7jM2LXXK14not8D/CBaVnkn+oQ5fhX7Vbth/Dfv/QXv/Sx6Y/qqq9jPCe99Z8PPn1zd4V0WQ2j5Bkgz6stNsm9dd+0S6041vF4HZr33tuFv0Jv5bHg9eIntDDa834h/j2YKOy6zf8+GxvmqtbAPiEiziPxmIDUtA/8IdIas/argvbfe+1/13r8RzWh+BvhdETnUsNp7t3xvv32FYf8tWmG4X5Qp+M/D8hdyvTd+d456lvNCMOO933iBYyAidwUC0IyILKHZTu+VPneFMQ+IyCdEyW7L6HRFbcxRtHJ0JQwB81sy8DNsroRc7TW/C30AvlB7IEKDc/9zOa4Mm/GaCoSh3PQhNJt6f3jaAi2P3tKwXguamV31hHF46nsr8Idb3to0dng9FYIg3vt/8N6/3nvfjT4hHgTuv8x+7QUKwDPe+wX0JvRStQ9ZpZ6ZEm7qlyNTXEBvCrV1m9Hy6OXWbWR0Nj3LulfCMfQc/JPGheFh5/1cokwVKgAfA37keW7zUvhX6Pd2l/e+HSUtwcVl76uC937de/+r6Hzj4ee7U977Se/9t3vvh4DvAH5NRPah19XeUNqr4Rau7npv/J4j9LuceL77WNvVF/j5Gj6CzsONeu87gN/geX4HDfh14Glgf/huf6RhzHPofeJKmAC6t5zvnej83nPFOXR+t7fhgajde3+5snbWXugq8JoKhOhFfQhl3zWWYP4MuFFE3i8iRTRrOBJumoiICctjIBKR4iXmRT6AEm62PiH+IfAvRORwYDH+KPD7tTdF5HWhLNoO/GeUoPGp8PaHgXeLSjRa0LmCjzU8Wf4e8H0i0h/G/n6UdPFi4BmgKCJfHY71R9EgfCl8FHiXiLxJRPJhPy93bX00HNM9Yd2f5PkHDI8SWH5UVGbSFDK930EznF++zEd/EmXadj6f7V4CbehT+2KYf7zcHONlISLfLyptaRKROJRF27iYOfpcxvwnUpeRLKA3Reu9fwZ4FPjxcC2/D53z/NOrGPZ2EfmawO78fvSmfG/YXhz+Twxgwtjxs4z1bPteCGMB5MNYV7pO2tDMa0NE7gS+ccuY+TCmALkw5pXugW0o+aYkItcD39Xw3ieAHeG7K4hIm4jctXUA7/05dBrk58I2b0aJKh++wrYvgvf+Alre/0URaReRSESuE5HLlYCnUKZ6hmfBayYQiuoCvwOd15mUuhbrm7z3M2gG8TPoDeMuNtfcP4De6H4dnY9bB7aWpb6ZzSQZALz3n0RlEX+PlkPOsPlG+W/REt45tJz3vobPHkXLOx9G2WNthIn6gP8APIAGrafQm+bPXNUJuQK890thW7+DPrmuUi/Zbl33KPA96BP5BfQcPtu634fOe1xAy3HT6A31+ezn/0S/nx9Az+OTQBPwRn8ZDZT3/jRaHm95Ptu8BP5L2OYsGhQ++TzGWAd+ES2/zqLn8/3e+0Y5zMdls47wz64w5uuB+0SkhGZK/zIcO+j1fQf6Xf088LXh/+BK+N/APw2f+wDwNQ2VlR8Nx/HDqPRhnaubs74UjoXPD6PEkXW0LPhs+G7gp0RkBX2Y3Vry/5swzj3Ab4XXX8Kz41+jAXUF/Z9PZVXhgfTLUaLLJHAclb0gIt8kIo0Z9jeg8qkJ9MH7x733n77Cti+Hbwby6LW+gD5cDobtvjl83zV8CDgcyqh//jy3t+0hm6ezMmR4eSEirSg1fX/DTTrDqxAi8hPAPu/9P3ul9yVDhhcTr5mMMMOrByLy7kAwaUHLwY+j2qYMGTJkeNmRBcJtCFFvz622byUR+etXet8C3oOWiCaA/cDXb2HaIiLvFLXDOiEiP/xK7OSLDRH5jct8L7/xahlTRP76MuM9b3LRi3k9irJfLzXWN72A/XvRjznDtYWsNJrhVQdRhuoz6PzLODoP+g3e+ydf0R3LkCHDtkSWEWZ4NeJO4IT3/lQQ3f8xmkVmyJAhw4uO50VtzpDhJcYwDcJtNCu8iJYuIh8kWHJFJn97i21BjMG1FpClNV0nl8NXq5s/ZwzeWiSXA+/BRPhyBTEGcjGUK/jmIqxencmJ5PP4agVBuKoKiwhiIvDgnUPyOXy5AiLgPdJUxG+U0/3xAqxtbPo83iOFPK4QI8trtfOBxyP5AngH1ulx5nO4phzRehVfqSJRpMddG0sEYqPnIArPxgLeOiSO8UlS33ZzMd0XEWHdl6j48othB5chwyuGLBBmeDXiUjfWiyJM8FL8LYB26fZ3ydsRkyca2Y1dOa4rWUA8RAa8Q247jH/kaYgcpq0TTIRbKeGjCmJy+EqVqK2VqKOdZG2L3jky4Cxb4V9/C5WuPM1jyzAxjV1YSIMVgOTyRE1F7MqKLhPB7N2DrG3gFhZJbjuAufcJkAg5uJdoqYRfWUGKRezcgo7RXsCvrxN1dSHNRZIz56AC4vN4qWzex4rF9PXhFhaIrt+LHxtHRkexR48R3XoIGZ9CmppIxs8j+QK+UoHbboL7HgfvidracKUSmAiJDGa0j2RiEokEv55g+vuwc/OY9la+sPCx5/7tZsjwKkNWGs3wasQ4DQ4mXIV7iZhwKXsH1QTTrR2Y5DY1ZpFIA5N/6Cg4q5nj2pre0Pv7MNftpvKWmzQQdHbgl1eI9+7WgFaDdxAZ4h0Dm7Ydz6zQfP8Y9ugxiML6DZmhT6q4tTXN8ABz3W6YXcBOTRP19RB9/jG8tfhqBXf0GHZ6Bre+gVtcIto9gpgIVyrhrQXv8At1m1BfrbAJXi1e7cwMZscArjmP2yjjT58DEdxjT2Hn5vHFPOaGg+A8pqMd7g+SNxGkWMS0tem2qhWS8fNE+VyaGa69fjd4hyutQvPL7vGeIcOLjiwQZng14gFgv4jsCe4zX48Kwi8Lbx0b77qTaM9O3Omz2Ll5oqYm/ENHkUIwxBEham7GDPQjh/eB84gxJOPn8VOzFB8/hzm0n+T8BdzaGsmpsU0BrYZalpb+fVw176a9naf+4950W5LLQ2Qw/X145/HWYtrbsafOUrllD1FrC8nZcQ3Q99ySZpG+XMZXKriNDezxU7iNDZCIqKkJOzuHXV4mHh0hHlS7VH/PLcR7duHvuQV/zy1Ezc1w983YqRlKu1uImopIcxNiDNHN1wPgxs5BNdGSblsbprUF09GO6e/DLS7h1tbS4yAySGsL8egIpreHwl8/GN6LkMqWQJwhwzWILBBmeNUhOPB/L+oo8hTwJ8GR5lnR/HdPYI+fSjOX2s3cl8v4JMF0tOPW1lh4+178kyeImopE+/cQ3XRQy4PtreBqc2eX+NeQSMuDSfWit+zMDHZ5mdFP6OeipiZMbzemvRXJ5TA93eA9trRK3N9L4fgk9PfqciC674m0bOre/DpK/6Q+JRrvGABniTratfQJGrzbW4luOUT00NP4xWXM46fIXVgk6uuB+4/iraXlYw/i1jeQOCbq6SZaKKXzfrJcwvT24Eur2NIqdmlZz5e1RF1dmP4+zWKd1TnMJMGvb2B6ezVgJ9VLn6cMGa4xZPKJDNsC7dLt7zLv2DSHF+/dTTJ2DtPagl3Wm7zp6daSHuArFaJCAVepYvp6cHPzABoorNtUdowHd4AxuNm5Te/FI8PYmVl8ucEhTqROyIlzaUkz3jGAnZklam3BVxOqdx7E/P3Dm+YTTXs75HP49Q18pQqR4KsJZv8e7LETOr+4fy++KQ+nxpGdQ7gTY/hyWTNBESQXYxfr5VPT2YFdXMIc3KcLnMOdPou3lnhkmOT8BeJdIySnz2jpd3AAt7wC1uoxr6xs2kdEMG1t2JUV7pe/ZymZycgyGa5pZI9zGbYNJKfcL9OrzSyS02fAWSWpBLjSKjiv5UJjkD2jmK4O7MycBizAbWxsnnsTwVcq+FJJmZQmQmLdlp2a2RwEIWRLif62lnhwANPZSTI5hRnohzhGBvsxn3kkXR+AyOB3DmFn5zSLrWqg9m+4CXdyjJnvfIOyRStVotklJI4p72hVtugbbyXq78WtroIxmL4+Jek0N2OXlnFvuhXmFrDHT2NPnkEKBSSfB+c0yzURUigQD/SRTFxA4hi3vo5bWSG65dCmcyH5PBij+3YJ8lCGDNcaskCYYVtAcjktbwJ2di4tIaaZTGSI9+7WUt+tB/HWYd9wE/ap47idg7o8BDApFIhaWlKijOTzuNIqfqOM6e9FWute3Wa0oRVfTYoQXpu+PkxXB76iZJm1992JnZrGzs7h2prTUmgN8UAfcn4Kuf0GLeXu24MtrSJfeIxoz04GPzUBkcHNLeDm5rFLyxQePoUrl4lnS9jzkzqfVyxiZ2bwSRVfTUAizL1P4Hf0EeVzRDcdwK2u4isVkgtTOn85foGouZnkwiSmsxO7sKAPCoUCnDhL1Nqq53BkmKhQwC4u4splbPeL5VueIcMrhywQZtgWSLWCtcDiHRLnkFi7ZSVvuxW/oiVRnCPq7CD67COaNZ2brLM9RfDlsmZWDeNJPk/U18v5r9mtcgvnMV1dSqgBolsPp3NnUYsST9ziEogoS7NcpvVTTyhh5sB1SFV/N86xJVMz2IUFzNQi8d7duDPnMa0tWr49fRY7foGopVmlDbkcprUF8jnikWGYntM5O2dJJi7Uy7NhGRIRrW1otns0SEu81/e8U2LOgpKAfKWCFAp45xER3NoarlQiamnGTs9ocI5zRK2tmLnVl+YLzZDhZUQWCDNsS5i2NpUe3KDzYvHfPYzkNSjaphx+ZQUz0I+vVvDD/Tq/FiBxjOnqIt6zSwPI3p3K8Dw3zuCHHsVXKpi+HohjzTxFiOaWkdhoRtbcnM4j+tIqUUuYuxsZ1O0/cxL35HHs8VPE/Q0N1J1F7riR5PyElnW9Q3q6NHglyvB0Je2w40ol3IGdkCQk4+exi4v1ecbDBzBtbVrC9F7LoLkYv7JazzYPXJdmzfHIMNUvu71+/LkY09WpAbJSRfJ54l2jSHMTpq9Xg3Nnh657xRaBGTK8+pEFwgzbBlGhkN74faWCPz9J0poHVE/olldwa2vkT05il5axe3ZgOjtwjz2lZJMWLfN5a7FLy/hYRfj26ROs3LVT583W1jTbm19UDWJPt5YSL0ySTE4BYGdnU4KJ29jAr68rY/SZk5vkGNFNWqLVHdSA4h98QufkvCd5080kZ8aJikWVYtTKvN5j2towU4u4pWX9u7sLRFh/7524Z05jV1ZwG2VMV5eei0oFKRbwDz9JPDKMaykiuZh4dAR7YZLc3z6c7odb36ByUEu+prVFS6hnzmFn5jTorqzoMZZKEJuX46vNkOElRRYIM2wbuEpVs7O2VnyS4DbK5M/MYrq6iMZncKtrSkhZXyfes4uoVIFCQcuh1iLFgpYqJSLu79XAFdD6+dNp1ug3ypjebqJ8Djpa8bsG8U4DVLx7VMeI6gEitSjznqhY1HGcRS7MQUerCtsbAqSMqaON+YwGJ1cuq1BeIiSn2RkECUWSYN96G6e/93rigX6a/vx+zJAK/svvvE0Dsvdg9LhNdxduYRGeOI4vl0nOjafEHkR0DnB0iPyxCc0oB3pVXF/btzgm3r1TXxtTl5tkyHANIwuEGbYNoqaiZl5z86y++3YNZpPTSuwIlmA4C4UC1R0duCeexk5N477kVha/9nXY+QWimw5oFjg7F8gv+i9i5+axy8tExaIGn5lZZHQI2ajgjhwLgUpITp8hHtqR6gbLX/V6JZyIYAb6ceWyuszk8kpoGb+APzuxqTRrl5YxB/dh9u3RgJnPU/6y14GzmB39JGfHkd5uJcF0dZG790n2fGQSmorEoyP4YoF4oJ/8YkXnJLu66udmYUmJMkk1zUKjYhHT1aVZ6KBmlm51TUX9Y+Pqv+o9UbGgx7JRVllIkqRM2wwZrmVkgTDDtoBEEW51lXhkGNPWRttnniGZnoVbDlB+5x1450kO7QbATk0jX3hMy41A9A+P0PGR+7SsOj6lmRCA91TecRty62FMd6dmk8GLM9q7Sx1lRFJbM7ntMEjE+a/ZpXOQPd00n5hX5ibg5hc1eLW3E3V2YHp7kJZm3OpaKseobdc+cxI3dk4ZrLtGKI6vqCykWiXeOQLlClE+h3S2K3vz+CmS02dIzo1jnzlJMjlF7vQUJ//TXbh9I7jSKpLLIcZgensw1+/T4Hbz9Xjv1ce0pQV5+Cnc6bO4UkkJMXtGVU8YqSWdNDfj19fxSTV1tsmQ4VpHJqjPsC3QLt3+Dbd8D+7I00Q3Xo87ekzfaLy+axle0L7VzLClvQ2/uoa0NIO12Nn5tBSJdyqpyOW1i8PqKkSGqFio25CBBsfuTjXDjgzx7lGSU2NBnK/ept4pSzNqa8Ovr6ssI8wfIoI5tB/75DP1/Y2MShgOK8vUPnUc96ZbyR+fIJmaTkX7jebeNf1jdOth3KNPaqm4UMBtqLtMbR9q843x6Ahubp6ov5dk7Gx6juKdw1BNSM6rxWu8Y4BkcirYximzNmpr496Nv2KpPJkxZjJc08gywgzbAlIsIOeUrBKtbdSzOmiQVPg0CMajI0SdHUhHO8n5CeziIsnEJG5pWfV3gaUp+Tzm0H7MyCAyNJA6t7j19dSHszbv6HOxBgrv0iAogcDjrU23XTPQ9lWdm4v37CJqbsafu5CyUHVFNeJO2ov4M+eJh4eIHzqmpJwG0b45cF04PjUFj1pacI8+mZZbox39qbg+uvkg1XfcoWSbzo50/tC1NytJJ9jI2QtTqQzDdHaQTM9qIK1WdJ5wxwBuZSXVbmbIcC0ja8OUYXugmtSt0xaW0hu06euD3k7s0yc2ZYPEBjszsyVLdOnfUUGlFm6jjBw/rQ4xA/2ASgaijnYl2jgPnW24sXHsyTHVLUoE3mL6elXfuLqqAau9HVtaxfT1YKdV8C5xTHL6DO5Nt5I7ekZ1gqtriAks0iQh+uwjOIJ36qUqOEbJOVEIbNXDuzD3HtVSZhyTnFXyjV1ZgSPHaOrvJSmtauBcLqkF3eP1DNo7MF1d2JkZopYW7NIyks+z8fabaXnsPL6pgD03oR6oF16arzNDhpcTWUaYYVvAF3JEHW0qYF9c1Lmwzk6l+Z84g3/DzbqiiAaHsbP1rhQEj09Qf09nceUy3muLImlq0jLm4pL+XSxo0F3fgCRBKtW0RGmGBpSUI0IyNYNbWkEKBUxfHzYEH19axXR2akNfazH79hA/cpzSW/YjA72YDp1DbCTQ6M41NNMlkFza28E6tZKbX8TOL2DuPZoyR91dN6ZzmEiE6enWuVOU9WlaW/B7R1R+Ebpg1Ig8UXMzEkg2vlym+aEzShQ6OYavJloqjbJbSIZrH9lVnGF7YG1DmZ7eExUKKpZP1F5M8jniJ04jxjD77Xfr+uHmLnFOWaK1bDKpah9CwN1xSF1mSiVlUAZPUV9NVJ6xtgbWkpwdx999I5LLUx3qDl6mTUTFghpsl8u4YFmGRNoHMbi44D3u3ATSVKT1sycAkM52/PqGZnBbBeuhlVS8Z5eSZEqr2M5momKReHAg7QqRnDmnq3/+UUx3F6ZLBfB2dlYD7eH9mNFhPY4jx7Bz81rynFvVeUYR3Po6fm0dUANzOtpw84uq17zzBuLhIbxzL9EXmiHDy4csEGbYPhBRlmVPN769BSkWkdsOEXV3YUurRHt30vvb94YAGbo4JMrCNPt2B9JKRDJ2TrOnIydxgdSSElK6u3Crq5j2dhXeL6ugXb54ROcpv/iYBpf1daS5SYOlCFFrC3iHObQPyefrLFEROLwPOzePW1rm/LsH8aU12DNKvHNE5yj7+jSr7OpC4pxans3OI3fciERCtFFF9oyqh+mbbk1lFaavj3hwB3Z+ATu/oPKLQ/t1/554Gjt+AUIgi0dHMMODaYcL1UTu1LlQ1MA86W3VJsMbG3DvETXnzjLCDNsA2VWcYVuh/NabcAuL2KdPYGdn8Q8dJTk3TlQsqAzA+zTj8xNTajR9bgJ36kw6RtzfGwynlXDjqwnm0H7V2QVvUbe+sWm+bpPoPBBkah6jABJca+xTJ9TlJY41c+3uwkzO6XaShJGPalNh9+RxLnzViGaiCwv4aoJdWFD9n/e4UomV3c0gEf7YaeyxU9pt/t4nNMNdWmb5LXu140VHu5ZBu7rw5y7oA0IcK2N2dIh4dEiPpRxaS+0cUSLNRuiqEQKjfPGIjtPTHUrMOTLWeYbtgEw+kWFboF26/T0d74NCATs7q2XI4HxierpwKyXN7PI5ne+bmdU+hK0tuHKZ6htvJP67h+o2ZjWzbWPUm/PgPuzx0yp/aGnRnn+hkW26foOMAUhLjPGOAcjnsJPTwY3FQRThNurtmyQXI0a1eqavDzs7Wx8bUhNtaWpShxzvdA50aVl9RKuJZnyHD6jIf2KSaEc/bno2zUpT39F8vt5jsEYUSk26EzUNSEvHsa5z834nqsj5AAAgAElEQVTk6EllmAbSjunp5ouLH2OpmvUjzHBtI8sIM2wLSD6XklGQCDO4A6xFTIRvayG57YDOzU3PYC9M4splTHurljatJf7MozpQaNmkzE+PNDVpYDgznpJO3PqGur90d4WNb44Dx37lllRmgPfYmVmSM+fw1UQlFUUl6ZjQzsm0t+KrCVFPtwrXFxZS39RG6Ye3oUu9s8ElZkHJN8GCLSoWkZU1klNjuI0NkrGzm7SOEJoRBxNwvEcilYDUjkXiWM9NVyempzst6fLYMcjlVD/Y1KS+rM7XvVIzZLiGkQXCDNsCNbanEmYc5X39mJEhbS90/BTxoyfAWkxHO2Z4ECTCLpcwB65LvUBNVxdRS4taiRn1Co16u9PuDxLniAd3ILcfVpH6bu0moYSbemXlwHc8oOL2YD+Weo06nVO0i0tKREkSzMG9SFtbCEhCPDqkAvy2tlRQT2RUQzjQjw2idnP4AGb/XuKhwXACPORyUNtWTY8YmWDYHem+h6wPqPueguonrU0JRrWejqkkxVoNont34asJbnVV2bkv3VeaIcPLhqw0mmFboF26/Z1v/EHki0d0Qa0Z79AO7ORU6qgiQfsnJoLr9+Iee+qiBrlRayuuVMK0tWnGWHOS2SjX7dTiHFFri5YYjUEOXYeMncfvGiKaXsDOzte73Eemrl+E1G3GdHdpAK9qqyNbWiXK57SFk7XaOT5JiEeGScbPp2XW5MIkUbGohxkCl69WIDLYt9yCWa0SHQk9B2vHiHbnqDFfgTTY+yRRu7dcTucU29ognwPrcKVVoqaingfQ+cEdfchSCbewyL2VT2bOMhmueWSBMMO2QLt0+7vk7XVLs9sOwxMnlOjS1YEf7Eem5nCLS8x/4+10f+QhxERqPRbaNkXNzUhzcyqQr83xmUP78fkYGTuPXVzSbC+pprZpkotDF3vtFbipk0TYn2ezegPNzly5jOnvw07PBFNrbTZc8wet2arVypdpN4u+XpZfP0zbZ0+pG836hnabqFQwO9QEwC+XdE7T2YsCM2gPQ/vUcUxnJ8SxEnRCRiv5vJZfI21LJXEuDdT3Vj/Fsp/PAmGGaxpZaTTDtoJPEtyX3ApHjuvc3+CAsjCfeFq7PVhL758dxVurMgDAdHUSD/TjVldxCwvY46dYee/rArOyCdmowPEz+A0lj9TE6rXOFL5chmpVy49DgyqXCMxMgnawsS1To9VbDbV9sQ0eovHwkGaN1kJBJRxRT7cSWoJWUvJ53PwCrZ9+Mu1U4a1Nu0fYiSkoV/DDA3WPUdi8PwDVhKi1FSnkkdgQdbSrDVtbW535arRUGzUVMfv3EB28DsnlXuRvMEOGlx9ZIMzwkkJERkXk70XkKRE5KiL/Miz/CRE5LyKPhp+vavjMvxOREyJyTES+4qo3Fmm2FP3jo/WyZH1QZUZGQuW2fUoYCXCLSyoDEEnnz1r/131pUEnGzuHWN7Rbexxjz1/QcqS12qleREkj3nHnX6nNmiuXteS4oz8tp5rOjosF8o0IWaNPEqpfeqt6oC4sIcZgJ6a0+0S5nAZxOzWNXV7WAL6ygl1YIpmcSskzOlaVZHIK9+Tx+nZqQbhhX+zJM8hgv2aDS8u40ip+z3Aq6o+Hh0LrpkiXTc+pQXiUJYMZrn1kpdEMLylEZBAY9N4/LCJtwEPAe4GvA0re+/+8Zf3DwB8BdwJDwP8BDnjvn7XxXbt0+3v6vg5frqibS1sL9pmTmsEN9GNn5zD9ffhmZYBGO/rx84tISzPeOez0TFquNAP92Bkl3eC9ljcbO1AEJqbk86pNhHr5syG4mM5O3Noa0e5R3MkxLcEmSSqT2GQGfjlcoox5VRBJt+8rlbQDxdZ1Lrf9TSXdhr6MkosxQzuw4xfw1nKf+xuWXVYazXBtI8sIM7yk8N5f8N4/HF6vAE8Bw8/ykfcAf+y9L3vvTwMn0KB4RbiVElIspkEQEbzz2Nm5IJwfh0U1kHbTs0hPF8nkVAh6HrN3p0ollpbVsaWgjWijPTvDwahUofzOOzRABDmDt5aouVl9P4PFm+TzuL1Dav59YTq4zWhXDF+ppIzQePdOlSIERLce1hchg91UzoTNfQupaxYR4cK/uoeorU1lD01N6RyhxDnI5TSjFCEeHYHIsP6e1zcMJJu2k3atT09uaCVVqeAuTKk8JBLInqMzbANkgTDDywYR2Q28DrgvLPpeETkiIr8rIkGUxzBwruFj41wmcIrIB0XkQRF5sBpVERHc8jJ+/EKq4ZPQO88uLgLg5uahqRiCjFOCTCjvuXMTuEpVSSvDg6EVU6IONIGliUQUP/0IoFmgXVgiam2lcvf12JWVdF7QVyr4h59MvUUlVsG8DPRiRoeJmrTT/fr+PpV1DPQT7xpFxs4Tj47ourffgOTyujwYhDfaw0kcQxS8R4cGGfzl+9QNx3ncRlmP+/yFVPpQC2xucQm8o+nP78d0tKfzhWb/Xj3OBulFo1GAr1Yw3V3IyGD9O8hnc4QZrn1kgTDDywIRaQX+FPh+7/0y8OvAdcCtaDOfX6yteomPXzLv8N7/lvf+Du/9HTmfVzeUSkWDgLXEOwZUEN/cjOns1GzPedzikorrJ6bwtc4RaKmzZoxdM63GO6L2NogiJr/rTny1okzRmuOKs2pW/fePqtC8WZv7SpxLMy7var0DHW7sHFSqqZdn7v88ovOJpVXsxBR2uaQtmpyHx7RJbzJ2Vvcv9DeUfJ7876wi1+9Ls7TkwpQG/dKqskO90zGAqKWZqK0VO7+gjM/1dTUdCH0Mo6Yipqcbe2KMaGhHSorBu1RikWaMXR24M+eR198UeipWn/c1kSHDqwVZIMzwkkNEcmgQ/LD3/mMA3vsp77313jvgt6mXP8eB0YaPjwATV9xGJPhKhXign6ipGFiga4gIMrwDu7iYBiNQHaAZ6FNSTch41t98vRJgDu6tC+K913ZKa2sM/Op9aWnSO5+K3e3ionqZOqflz7078UmV5k88rDt35w3EOwaU8ek8yYUp3EZZ+wQ6q/N36w3zd9Zqe6SkWneosVYDT7kMzrP+tlk4MQbOhfnLIOEIxgJijFqvOQ/VKr5cUYJPILxILsaeGMMuLhEN9Gn3iUjwK6sqjQgZde0hoTZHaE+c1jGOPINpa0Mku4VkuPaRXcUZXlKIiAAfAp7y3v9Sw/LBhtXeBzwRXv8F8PUiUhCRPcB+4P4rbcdbzV6SyZDlVavqyWkMLK5oEHD1xNJXKypSryEyNI8tqQ7vzHn18uzuqs+biRC1NOsY3qnwPanq6+bmtDWTtxYqVSTOkbz5ZmWn3n9UmalLy0T5nBJOWluI8jllsnZoL0TJ59JMzicJUVNTff9qmVlk1IGmv1dZrEXtSkHt2MKcohnoV5JMPqf+oKtraQlYtYCx/i1CcvqMzilai1vW+dGavjKdJ0x7GipxpvzWmyCSrA1Thm2BLBBmeKnxRuADwJdukUr8gog8LiJHgLcBPwDgvT8K/AnwJPBJ4HuuxBgFZTNGHe1KkKlW8Osb2nWho00NrL0jHh5U4gik82ISx0gcU3nHbbiTZ4haWvCH9gBolpSSVyLVJe4eVdZnpRqE7r11FmmhoJ0fQtYUf+EouIbyIgSrtiCIjyJ1r5me0YwSdX8xrS1IPqf+o8ViKEEq2zQKy+3sXLBsK6nUIaniNspEN16vpd3zE2kgdaurDdvOBV2jyj/sW1+ny0JJWPbvqZ9Ut6Ui7X1oSeXI/81DMNh/Nd9/hgyvemTyiQzbAu1Rj7/n4AdhakbJIuWyGldHgl1c0ozupoPIWhl76mwqWsdZTGcHbnU9lBdz2KVl7FteR7UtpvnTRzTTC/IJINUo1mQVjU4xUVubCu+LBaSgXe0xhqitVS3L1tbSUqgUC1CtEvV0k1yYUvPtSlWzy7ANM9CPL63iK1UtdVYqmi02N19kqP2cIKJBDW0z1SgVkaYmfK0zxo378Y8+qaXWfJ6ovQ23vKIyktIq99lPs2RnM/lEhmsa8ZVXyZDh1Q/J5XCnzqa2aFGhgF1eVlZoLo+vVnBPHNcsrrUlsCgd5AoaZEwEIthltUkrHJvATM3gQ8CUQiGdp6tlkW5jQ70/p6brAvbQyNetrBBVi0ihoMSceWWt4izOWi3VVipK7jk/oeOtrmsmaK1qFk2kLFdIy7reaunzBQVB0Pm/Btu0ml6wsUWT2beHtR3NFLyWas3wIG5qJm08HO8YgNnsFpLh2kdWGs2wPRCIMLWszW1sYPr6KL/pcN0ztF07rNvl5bSjhC+XcevrRAN9mq211LsxSCRE7a2qC0yzv2CrVuvyUMgjsfY4jFpaVC6xvq5BpakptIIym9xcpKlJHVmsTbWK3vkgc6izMH2S6HLnMd2ddWLPi1nFaRDrm55upKlYt3Y7P0nhUw/XT/GZc6ljTu0c+UrloiEzZLjWkAXCDNsCPjAcfbWiHRIiA0lC7tPabDceGkSKRRXOt7drxhXHqU4uOTeh1milknaA39Ab/sQ3Xs/GW27U7Ck0r41aW4gOqOYuOTehBJdqokHMO6LWVqKbr1eNX2dHKnuo9UqsZYKz3/p6fK0cGWmgXnnfbQBIsaB+n+2tarEWMsOX6OQB4JZWsPML2OlZqCrxp2byneoKGz7j1tdfun3KkOFlRBYIM2wPeK+tkVpaUsG4XVjQUmh3F8n5CZLJKZU7rKwQ7d+t5b62tvoQibZDSt1aIsPQHx2j+IVjGhRuv0G9Rw/vxp85r0zVSFSz2N9blzCsrmnrI+e1Y0RbW53s0lQE54laW+n50L1qB1cuY/p6cavrdHzySQ1+IvjVNSXCNEg8GrHJZabB6FuF9iZlj2o/QtFgXNMDbjXdJsx9SoS9+wYtvVqrrNTeHiXxdHWk84i1c/6s3qkZMlwjyAJhhm2BGvnFra4qexSId41i+nqhMWDUurs/fRIpFFJT6docX9TczPpX3KpjRoLbPYhfX1d5wYNPqEbvi0cov+H6IB/wIWNbUB1hIJnEw0M6T9jcXC/bFgqpUN1XKnUGq0TYmdkwZxlr8POeyX9+W70E23istUAHaVAzvb1atrRWO917h+lox87MBZmHT4MxEtXlELA5KDpL9Pkjelz5vIrvA6PVLZdSv1TT26OZd0a2y7ANkAXCDNsDkd7cTVcXVJPUHcYtLCrhpEEugYjKLZqb0zm3WlBypVWaPvUoUVHn7eTYGQ2wNXJKZIhaW2k6OZvOj3lrNShapyzUjnZ8uYwc3ldnnNYMugsFFr7h9ak9ml1YSjMxb63OH0oE1Sr9v/aFSx6qtxZ/142ps4wYEzK3Xs1oa7ZoubzqCFFpR63Rb9TSXO9sD/V5wsikJt+mtSWdF/Xzizr+QJ9KKiKDX1vHzi+8BF9khgwvP7JAmGFbwId+gHZpGUZ2aMYTGXV5CTIJRDCjw2qM3VRM/Ue1jZJNnVRqgS/K53ArK8GFpl+F83ffiORzJGfGMW1txP29Wl4NrM65b3uDyi/m5vFPnUpNtk13l75eWaH7Tx9LfUxNV4c21+3uVPuz1TW1YtvaKQIN4qanW4k5Dz6FO3pMM7RqBTc3j52dxZfLJNOzunxtTS3T4pyK68PcqSuVcF3taYkzLXU6mwZFW1pNO2T4JMFvlPHtLVra9Q727axntBkyXOPIAmGGbQE1m7aY66+DCzPpch9HDe4oHnvuPFFLC9LZoWzPONbsaEuXByIDxtSlF3PzGkzuP4pbLmG6O7UUWa7g1jdCx/qI3t9/AMnVS7Gyow9fTXDLK5APfQorFXAe09ONWy7hb7gON7+I5GJtm3SJcigiKl9YWsYM7UjLnVF/bzoHaHp7techqCxjbQ03pedCjEGam5AoEIfOTmAG+jF9fUTNzUQNc6U6gNMya2ur7pP3uONjqT7TPXGcaP/uzHQ7w7ZAFggzbAv4srY2YmKqnukB7sjTyuBsgFtdxZ6f1OwvSZQlCXUT7ZyyON3aGqanC7Nvj5JaXIPdmDHqY9rSHLSFKqGQQkGXx2ql5s6e1zJruYxfWw8EE/X6dAsLGmxOTegYI4Mqbt+KWoAODE6/osxW09tDMqbaSbe2pszSmlQjlHx9bd6wuRnbqy2j/EoJBnq18e/sLNLbrQ8Sm06oft6trKi7TTVRBqwIrlLl9M/eiR8bh+R59ErMkOFVhiwQZtg2iEeHtDQqkQabYG3mjjydrlMzkfZJlcnvu0szvkQ7z8eDAxAZFr7utrS7xOKbd+POnE+lERCs1CpVosEB7NQM8ehQGDNBiqEBbvDxxFrcRjkI8MuavUWipJrAIvWVCvHOYe2huLUJb0Pz3GRMdXw1cb6dncO0tRHdfD2mu4uoWMDOL2i3ja4u7aXY0U7c36ukoEeP4d54i1qrnTyTGmtTqWLn5rWf4pa+hIDuZz4HnW1IawumvZU9P/xFDaoZWSbDNkAWCDNsH5RVZmB6ukNW5FRW0NyczoNFoa8fEjH0u4+n0gTJ57H9XUguputPHsZtbBDlc3R+4ZzOGQbNXK3DPNbiS2vaxWFyWofM56Cro95o11okn9cSaKWKRErSmf+G29QztFIJLFJfb/u0FaKlXdPXV3eAcTaVeNjSKu7I0+qOMzSg/qojg6pVLJchjnXOMK8BP3d8Are+genrIerUfU3Gz+uDQ09Xur1GuNVVyOWYu3tAM+ANdaRJ9z9DhmscWSDMsC0gUYS3DnNwHwR6f9pl3Zi0rOnKZcZ+5g1qnL26lmZAbqOMPDMGh/cpuaWnG6KIZLRXJQVNTelN366s6JzZ8jJQ9+p0pVWYnlN5QVcn5NQuzQ/2pF0qJI7p+eNHlL1ZKOC9Z+1Lb7zEAel+mZ5uomIRt2uAeKCfjXffqf6jwS0HZ1U7WSphT5zGLi2zfFNvSnSx0zP1LLNmmu0sdmqa0huvq3ulekdy+ozOs7a3s1VnmNxyHT1/dUznMltb1CRcIi7TKjJDhmsKWSDMsD0QCX5lRf1GA4mlBhHZ1In9ul86pizNxmzGO7U8m5hJ+wy69XXkwSeZ/IF7dI4skErEGKK+Ht1sX4+WN4MVmlvfUHLMqgrSo9YW5Pw0YgzT//SG+vaMoXr3YVypRPET91+sx/Oh32GSaHZ6bhqamyj+5UO4+UUluQSBvFtdDfOOmiW2fOx+Vt9/J9fdXyAqFIhHhkN/wjzJ+HkN8kDTXzygkg60KbHp7Um1lZsa8gLR5x5Vqce+nbiVEm59XT1bbdaGKcO1jywQZtgW8IlF9mg/X19N6t0TgOSGPSTX79SsrKkJv76h/QGtav5AS6ZuYwM7M8eJn78DiY061TQ30/VMlai1BVta1ZKj0/k6nCc5O55qDDUw2CCez4Mx2IUl1R62tdH33x9RI/DBAZVCfObhywvSRYh2j+JKq5h9e1QYH6nFm+npws7NE/X2aNm0xhQNQc10d9Hyp/dz8i7tV+hm51Kto+ntxa+uaVBvblaCT5zTuczFJS2nQio/ady/KJ+Dx4/pZ/J57aBhsltIhmsf2VWcYdvAPn2CaNcw0YE99dIhYNYqyBcf05WC+4t2erdKrvFeBfmozKDvYbDzC1TffBO+UqHwlw+onjAXp6zOWv/AeKBfg1F3J7a0iuloJ969E7e0op0oIoFqBT86gPeeeO9u7MQkM99488WklIYmwBLnkNIaprcbe3JMxz91Fp9USSaniEeH8B2tmN6etGOG2btTg1SxiOnvQyLBLi0rwWVoBziPnZkJQWwdt76B2yhjhgb0oWB0ODUM8M6H/o11rWHNqBxrlUFqIsgywgzbAFk/wgzbAu3S7e858O0wu4BdWmbsJ+9k9499Ud+MTGrBZt9yC/Fng4VYzSd0xwDJ5JSWPHeP4sbOgUREHW0qSQgdJ1JRfkc7vlJl+gO30HYuofBXD6TtjCTOEdUytI0NdXlpKkIQpUs+j124CkeW4A4T7xjAzs4R9XTj5hdVZhHVDQBSMbwx9Wyu5iwTfksur0Se/bvUA7XWcinOqR5RIuKBPg2wQ4PYqWmdU/UO0911keG3xNoE2S0tc2/1Uyz7+cxwNMM1jSwQZtgW6Gga9PcMfRN2/EI9GwzXtsRxmuHYyWl8tVJvqxRQCyg+tEZylSq4epBEBHP9PmRhGbe4pO97p4Gjv4vFw210fvQRcB4xStwhEuT6vbjHnlIR/6WCYMN+poEreJJGA31Mv22Inj98QN9uasKvr6dNgqOWJuxyCXNwL358UvsIAubwAWSjoqbdc/P4JCEeGVZ2aOM2GzLS2pxg4wOCZsoN/QobzlWtDHuf/1uW3VwWCDNc08gCYYZtgXbp9nfnvqKBBbk5CKaoOcjUAoFESnYpl9P1ap6cmz7bGAwatH2bENo0pQG2McgB0S2HNCO7HML6aff5hvFMTzfS2kJ1qJv4qTHsgZ3wwBPpcdTmOOM9uyCx9aDXuI+R6hYlF4Pz+D3DRMtrJGNnN60XjwyTnJ/Qkm9XlwbvhnMWtbbiSiUkn+fejb/KMsIM1zyyOcIM2wKqkwsByHsq73y9at2SRHV3xmAO7gWUTFLrxo6zRN1dqR9pVCxQ+dJb6/KBmv9mY0C5XHPc0Ml903oNeNYgCJjr9ylbdW2NuW9/A2b/XrVBO3CdGomXK8i9R3Re84EnqLzjdm27lM/jNjZwb34d9tx5qFZZ+Oa71UotlyceHWHj3XfqRqJI5wYDwSc5M56WWOPdO5V52rDftT6Mm45xfV3LxdZq6TdDhmscWUaYYVugXbr9XfL2i7Iw09uTMkRry/0bb0U+/6i+39eHW17WnoUXJlVDl4uR9jbV1TVACoX6PNxzQWQudozZipCdSiSYHQP4liZkaQW7sIjp6sQuLKpJd2sLdmVFg+Ps/EXzd/HoCMlwN9x7JN3fqLlZmbS2bqodjwzjl1ewpdXwAFDENRybNGgvTXcndnZu83b27sYXcnzxqd9kKSuNZrjGkWWEGV4WiMiYiDwuIo+KyINhWbeIfFpEjoffXWG5iMh/FZETInJERG676g01lvj27MIGgknU1FTfl88/iuTU8aXGokymZtSPs7QKO/ouCoJE5uIg2GhH1pgZXU2W1PjZSHsElr/iNm38OzHJma/px87O4SsVksmpdNt2ZYV4aBB/ZlxJQA29Ftfedxf2wiTy4JN6KkI2J2k51GL6+oh3DOCbi9jDu1n/v25Hcnk1077lEPHOEcyh/WqrVutDuLh00bmwZ8c59h09Kds2Q4ZrGdlVnOHlxNu897d67+8If/8w8Lfe+/3A34a/Ab4S2B9+Pgj8+vPZWHL6jNqaxTkW37tZrhA1FbHzC0gcs/zOw9pyaX1dm+wePXbxYFszulp/w5rwfGvptBG1JrgNQdMc2t/A7HTE/b0UPvVw+pGRn/tCYG4qeWbi39xD1NJMdNNBkgtT+CRR0+4kgcgQ7xql5eMP6eZCebYWJO38QloKtbOzGmDPTyLW0/S/H1B26HW7kXOT2PMXkHnVE7r1ddbec8em460djxkdZt8PPpCSZjJkuJaRBcIMryTeA/xBeP0HwHsblv+hV9wLdIrI4KUG2IRws067TYRyo69WaP/jB9LVJJen+rrrADAD/bT+r/vg4B42vvr1mzu3p2OEnzC+GehXRumu0ZScU3N1uSRCYIyam4lHhlV+MTa+KTD79lYti4bu8rVyquS0n+Hobx9FBvvhxFkkEqK2NrVpK2p/QDsxxcr7bsf01dsypRZz4Tji4aHUcg7AP/hEnVS0uo5dWMIMD5JMz6bnrvnP7lMm6SYXHk9y+owSbzKv0QzbAFkgzPBywQN/IyIPicgHw7IB7/0FgPC7PywfBhpdqMfDsstDSE223ZGn06Dlk6q+HcytTbu2Ior+UecIk/MTSlB5/BlaHzyzWVJRC241ckwIUC7MyyVjZ9UhprcXM7wjdJ03mz8fgqjk8rjVVZLx8/hqgrS1hpUiZXyePqctoeYXMfv3bgowUaHA/LsOIeUKUVcnQMrktLcdTIN960fv1zJqw3xovGeXbibOUd2tp1ezPW38GxWLaeNdnKU60kM80Ifp7NQsuKaPbHSQSUu6Fxt0Z8hwLSILhBleLrzRe38bWvb8HhH5kmdZ91KTbBfdcUXkgyLyoIg8WPXleueEcPOOR4aJ9+xS0khoteTWN9KsL96zSyUFxYKaTk9O6cCRIbrl0KZtRcWiWpGFhrWIEN16GEDNt4OlW9zfq22aIhP0ijninSP4akWzt8jgK1X80jKmt1fdarzHJ1XtVN/SjH3mpDq9tLeHlk0bdP/FkyTnJtKOFvGOAaJCgWp7nnj3aBqcKu+4PRhiawC2584H39II+cJj9ay2u1OPJ4p0f02kJdb5VdzSsjJDQecV9+/d5Cka7xjQwF4q4V3mLJPh2kcWCDO8LPDeT4Tf08CfAXcCU7WSZ/g9HVYfB0YbPj4CTFxizN/y3t/hvb8jR6Ghy4I6vNjpGezZ8dQWzTsNOKa3F4Dk1Jh2jVhfvyiz8XG0yb3FW0flbTerJCF4kMq5SSDMxYkQdXZg5xa0PVEuVmaltSRntQzqrVOhfGuLZmErK8roLGsQ90lC1N6WdpWoZae1gCiRqJm4cyRT07hKleanJnETk9T6L/pIiDraqbzjdvCe5a+9A7N/rxoASJRmtXZunqipGBr5etziElFTEXf8tJqGb2xsOrZaZo0IbnklZIkGycgyGbYBsqs4w0sOEWkRkbbaa+AdwBPAXwDfElb7FuB/h9d/AXxzYI/eDSzVSqhXgunq0m3mYu06MTKkr4NmEF/329QVG8p7qVDeEW0kmI72tLeg6emi+PBpJaqEzvb0awcKyeexs/P40ipmaECty6oJhOwzHhoMGsMqPqliFxfVEi1JtHVTZLQpsDEkExdCpqZ6wqhQgEiwgbnpdw1ukjn4haX6Pnd10fQPR/GJpTC9jrnhIJ0PTeNOn9V1ag8KNV2kMbiVFXy1ks4nemvVmScymCPR72gAACAASURBVIF+pKlJGwYHYo/p7UUG+zE7RzR4Zxlhhm2ALBBmeDkwAHxORB4D7gf+0nv/SeDngS8XkePAl4e/Af4KOAWcAH4b+O6r2ooIdv8Ik//PG4gG+nDlMm5mDtzmbK+R2BI1FfV3czMA5sB1mlmdPINdXNTO9LccxC5oV3jT1RFE9h5ZXMGMDKWMTCkWcDNzbHz5LZp5FdUqLZm4QDw6osHXe50XbG7WkiQQtTSTnJ9Aav6f1Uo63yk9XdjFJc0ckwT36JNqyH3HjZi+Hs1aczkA3O4dSD6HX13FTC8g80vYE6c16/NORfO7RtNzYBeUNSuhnZOvVNLzaK7bpR0vGrNliVS3uLBEMnZOyTKZoD7DNkAmqM+wLdBhev0dX/KDSoJpNJwuFLSbfChz1gKQW9N+gdLUhDQ3YWfmkFxM1NmBm5uv25yBZj9r67iBbrUkO32GeHREPUdLJcy+PbixcS1fGoO0tWKnQpU3sC9Nlxpxu6XlzWSWkWFtk+S9BsJiQbtJeEc80J8uTyantB3TyTFtJWWtdpSoJtipmdRfVfJ57N2HkcQjn3+UqK0t9SAlMul+2JBJRvlcug1XLmPa2rDLyzqvWguMYV9NZ0c9KFuLL5fVazSzWMtwjSPLCDNsC3jvMJ9/fFNHBtPejogQ7dmpWVShgFtfxy4uEgX2qFtZwS0uUf7K2zRz7Ggj6upC2lp1Hs9a7LkJkskp3GNP4VdKmN4e/NpaSiiRQCTxlYqOPzWdit1NRzt4R+nN+3RerquLtJuFCHZmVhvvdnZoBltaxbS3apmyrQW/tKxzcqjEAYm0J2KlQnL6DLZG8AGipiZk/x6izx5RYoyIskHvVg1lvGsEKRbx6xvEA334u27Elcu635UqUWsr0tGu2bG72EbOLpfqJdv29jSjzZDhWkcWCDNsO9Tmu9zaGq5cZm1fD3JIJQmms1MDUW9n2pTXVxMKf/0wks9hnzmJnZ1ViUSkzXklF6cSCL+mwc8tl5CmJm10O1u3OTOdnZjODp2XjGPs4hLm0H6a//JR7RSfzxEPqyTS9PelZBQ7NU28a1QzMeu0ddJ6GVepKvkmjvFVJaxE+/fo59vbMcOD2krq/2/vveMsuaoD/++peqFznJ6eHCSNciRIGBAmGH6AwcASTDAIVracbWy8wHptY3udwGuwveuAbMAiGTA5GmMMIkkgCWWNRpNnenpmunumc35V9/fHuRW6+3VPT5B6uvt8P5/+9Ht177t1q169OnXOPWHrZuKJSeKHHtV4xIs1TtJdd4mmWwtDKge7iHr6YMdWok0dEPh6jc1NhK3NuPFxTdadd44BPfYkg00SZ9nWrOWmTBc0VgAmCI2VgVdewo4Owh0XpGnLJAwpf+0u4od2Z9rg5TuIHt2DNDSot2fJr7GNjhLU16tnaHMTYUM9BKI5Pr2nKWGINDUiYUDQ1Ijz5tOkrh8uJhoYVOF1wRbCy3YQ7dwNgVDpPkal+yjxiZMUNm0kOt5D0NpCoXOthnb0ntCq7yWtExgdO44UC4TrOgk3bVBnnIZ6ODlIUC4TDQ1ROXhY1zQHh7Qo8Pp1qunuO0TY1orc/xiI0HeTJgsImhogckw3ltL4SBfFRP2DBBdfkJ6HVBtM+sSOcMd2wkb1ao127VHN11ZWjBWArREaK4ImaXNPK70odQpJShIdfsVG1r/3B17zCtI4wKCuDtnQyfS6Zgo/fixzeEmK1frg+3hsTFOd+SwvuDhNRh3UqKk1SLS4MNCwg6TCfI2usyVCK1zTjpTLufVDDYQP6uvT+MBoZBS55lKCyWmiRx7LEnYHoaaBSzxGk4wwxRJy2QUEQ2NEXd3psSOi65yjo2n9wCQOkjCEq3YQ1RcpPXZUM8Uc71kwOD5c05625xN92xqhsRIwjdBYGdRrUm0XuzQGrnLgEOvf+wP1gvSaUhKbF4+OEh86QnjHgwRNjRTWdWrqs8o0hc61moB7eFidaYoFn4FlGhc7gpqyBsJPTYGo8Atqa7KQjCRUYnxcqz5s30h40XaQgOktHRrPGDs1uQZhWkLK+ew14cAIMjymoSBJ/URAGjSUI2hoSA87qK8lGBylctBnpnnypWkR4HhsjHDHBQTbvKfoJRf4slSOYHiCwW01UFszrxDMrwFGJwfUwSawlGrGysMEobEyGB3HTU9R2LIxFw/oNTefJUYCIb7xOs3n6YveBnV1RL19VI4eI06cUdpbNA7QOXUumZwknpyksHFDJgBjrT0YNjcRdnSolrWuI11PkzD03qEx0nWc+MBhze5y5wMEV+7QBNqtLWpK9cmz45GRNI9nfOIkeJOtlLRSBq3NGnrR0pyGe0RDI0xvaCOoqyOorycYm4YdW3X9MQyJ9uwnPnhEHWmOHCdsb1Oz6+79tH74Tir7DlDYsmlmrtQgTD1Qk9fJ2mDU27u46hqGsYwwQWisKCoHNHhcymWCmhpdLwSinl4qz7yawn17oFLRKusTk0TDw2rqFNGUaIEgYxN+zS8rk1ToXIubnCJoaiJobCRs1+K+0cCABuiXS3CiP11PS3GO6MRJXbdbqx6j8UO7Cda0UTnSTdjSkjrMJAT19QTr1qZhHG5qiri/X2MC40hLLdXWMPHS6wlbmwkf2oeUSkipRHD4GMHxk0THejRTTU57JIqYvmyLmoh9sWKCULXJaQ2VkGJJnW1aWtSsWypmSQfSzD3OhKGxojBBaKxI3OQk8cSEakQP7CLYuonwO/cTDw8T+VyaYVuL3tSvvyIr2hurRiZhqBpVXx8SiFZkcLGWMervp9LTp16ol15E2NRE5dhxxGtpgMb5taumhoiaXXfuo7B1E4UN63AF1bSiwSEVWOWyruGJIJvWEx3u1rXM9jZdr/RrdGFLswrVoRHqvvGAllgaGcFNTmqAfLnM1IXrCNtb1ct1fBxXmSYeH9ewh+/fn6ZZCzdtmHHOomc/iaC2hsqNV2vGmyjS1Gygx3DxhRQ2b6KwrpPC1s0UNm00gWisCEwQGiuGkVffoNrbti2aQSUx74GmCcuttwU7tqvTRxDCnQ9k63DXXUrhgm1p7s8kKL+wdk1aaSKoqUF86IGMTWgibxEqR47qvrdtJiiXqVy8CdmswiYeGdXE3gcOUTnSTbRnP8FVlxDfeHVWKWN6SjXIx/aqt2hbK1QqhJ1rqRzvyeL7vJMNoOuZDQ1IbY3Oofsowffuw01MaHiEBEihqKbVqy5JvUCDmjLxsZ60+C4iFH/0KJUrtiPTcZqkHBert2m5jFQijUsMQ+KWBqLjvTOrUhjGMsWuYmPF0PSVBwFwwyNUDnerWXG6kpUTgjSkIt69X8MNNq4nvGyH5vQUQR47RNTWAIFohQfveekmJoj6ThD19uIi73na2oobm0gTUBc2byA6elzX+CYnkR/cT7R7f5rRxsUuS7MGxA88SnD7vRQ2rsdFsVbKKJYI21pxU1NeUAfQptlg3NSUOvCEoZowa2uhUCBobsKNT2ixX7RSxvRVFxBs1jVNXKxjPagFh8OLtquGODGRessCHPmlaynsPEB4107k2suRBnW6kUJRHXqGR6Gtmbh/AOk6nnrXGsZyxwShsWKIx8fp/p0biE6cRIoF4rEx1dxEiJ51jWp0O7YTdqxRgVhbS+XIUVzXMU0zVtLSQrLrIGHHGjV/bt1IYV0n0dBIGlQupaKmOZuYIOrrY/RnnqweqSOj6khzw1XeS1QrRszwtIwiCus6AU0Q7p5+DW5sjKCtBdc/oELV5xUN6ut17bLrGKEvxAsas+imp4j6+4kHhzjx7M2aAWfnbk0kMDVNeOfDSOyIx8fTlG7J56M9+zMB5jPgFDZtZP377iAaGtF6ibu0T7i+UwVebY2uU+7eryElY+MaU1haoCCxYSwTTBAaK4awsZENf3kHAHLpBT4mUM2bpZ1dAIxe3EbfT20nqK0h6u2DWNOsuclJLYfknObmdA4JhGjnbqIT/YQXbtXg940bIAhUC5vQ/vVfuldNmif7Na/orsPI2IQKko41fk1O0hJMUf8Ahe1bifr7Ce7aSTw0QryhQ/N47thOZWObrgVWKsT9A7quOTyc5k8FUqEcDw/T/PG70lJTulErSFT2HahaXSP5bGH7Vq12X1dH5XCXn297VpqqtZmoo1kz4tT5/KZNDRqrODamFTR8gV/DWM6YIDRWBiJZMmsJCAZHAV9VwleVd5UKDfceoeWjdyL1mjgaMsGQHys63qNCJNAwiGjvQeLONm2PY+KhYaSk+UvDTg2bCMpaEzEaGCDuHyC84hIVjs5pdprODtW4JieJG7RIb9BQr9rgfY9o2MXAMMGDe4kGh+ZUzUhzlPpySSlxlJVWmk3i6ZkI0SD0MYwlKgcOaRo4L2TD9jainl4K27aoYN93CHf3Q1QOdxF1HdXqE0Mjmi8VsmTehrHMMUForAikVNSSSZA6pSRZVIKasq7p1dRouaMwJPL5QcMrLiFobkrNl5MvfiqFrZu9JhlT2LqJoLaGoFREDh6FQqjVLCYn03qBcd+J1CSZVJsgDIl27lGz4/p1GgwvgrvhStXEDh1VR5WC1k3EOYKLLyA62Y9s2+RrEkrumHLhC8l7Ea16D1m83+yAd1+pvrBxA+ElFyHXXZrGOEoYkmaWkgDaVcBFh48g111BYX0nlec+ibClmaC5kbC9TbXkvj4VqJUKhrESMEForAyiCHfdJQQ1NT4UQVOqOR82UNiySdOTSZAVn3WOeNdeXScrqsCq7R4lOnpcQw8qFXV8GRtDaspISxNx30lcFBOuaVcNanhYA+ubtMpEeOE2AOLRMW9inCA62a+1/WpKyA8fwg0NE49PaK7SDWuINnVQ2LyJ6JHHcJOTyIBWgHeTk6l3aJ6gpsZ7dDo10yZerN4RKGxtJbxsB1Iu0/MrP6Gxiyf71XP2/l3qBVoqqmlzWkM3wuYmoj0H0hJLPLiLysZ2ah7ugmIJNzlF5eKNaupFkxOEHR3qjGMYyxwThMaKwFUiCn3DWspocjItoSSFAkF9fbqeV9i8IY3XA40bnLpCU5BFvX3ENaqhBQ31KlQbG1Vb2riOyoFDBD6ODxGivhMAyJUXI/V1KnAOd/sJaSaWsK0FEhPsxJSWPxod03W5bZ2wrwt3z8NUDncRNjVRWNdJpfvovNpW2NSkVR/CMC2oe/zN1xE0N+kaZFMD8dgY0c7dBLU1rP/EowRXXazre42NhJs3ahq6qenUo3b8uVdpXKWPWQxqNV2d3P8YUzs2MPKM7biJSYoHehi7ehNBQwPh5o3EQ0Mz1yYNY5liSbeNFUFzqdNdX3kWoNlRwrVrqBzp1rWvNWuIh4YI162FySkNfi8U0mTaaWLrxERZU6PhD0+5kmDfkTTJdHDlpcSP7KawdZMG3ftk1mnB37papKmR6FCXrlNu3Uh8QHOAhh0dalJ0TlOyNTcQ7T2gWlyxoI46yTzmw6eFC9e0E4+OeeEea2hDFCGXXYgcPJrGLKYfK5WIJyYJL9wKpSLxrr3sf9dT2f4nP4bYEa5pIx4ZTbVY0NJM8fiEaoylIrJlI9Fje6k898mU79+fjvuDox9nMD5hUfXGssY0QmNFkNTq09dTTG9bS2HbFp96rKhmxsEhKJd8vb2KalCNjRQ2rEs/G7Y0a+WKrZtxdz2INDbo2li5DHsOUFi7hujwkZnOOUGo6dr6ThAf79V9hgHxwSOqiXWu1Rydfl0v6u0l2rNfU5y5ONUYqwnB1JHHC3QplYj6TtD7s1dqou/aWoKmRnWs2XsYwiBNDu6S4Pso0gD6A13EtUWkUODC9z6qBX6np4iHhlUIxhHxM67RtcOpaX1wiCOktpZ430GkWKJ012PEg0OaNm5qCnuQNlYCJgiNFYEU1ElEiiWNF7xvN1FXN25yUh1kCgXi0XFcfW3qfRk0NRCPjqnmCOopOTSC+4lriNobNSF3VzfS1KhB7s5BMQvJkDDMBFUcETY1IQ31BO1txBMTWaaYpOyS34dOVNJx5jODJmuBYWurFhROtDwJWPvpRzTRdxRBna+8MTEJU9Nq+k12F0Wp96mLItw9DyMN9UQDA+k8ZOtGwou2aW7WOx/SnKVNjbhigfjG64j6+wnXryNoqMddspXg4gv0uDvXELfWn+1XZxhLjglCY0XgKmraDNpaNP4uyZgiouWMQLW0vQdVoIgQD43MCDuQYonoJ69B7nwAqcT0vk6D8Cv7D1I5dhwgreAO+ArtgZoOEzPp5CTR0WOLmPDCmpQUCuppWixqPtLJSeLBYRV+QDyqRYalVCQ63I34ShWRD2kIm5vU27VWy0fF3oyKc9DRRlBbS2HLJtV2Tw4S7z9M0LGGcON6pLGBeHgEiR3FY4O6xjo+TjQ4RNjVS9RUoyEljzxG0D96el+UYZyHmCA0Vg7OEfX0El5xSZorVApF9fqsrUW2bCS4aBuDr7+B+BnXED3tCsKmJnVAaWwkaGmmePduHerRfbR/6EeaZsynPZt+5pVqzpSAwuZNWd3BpMbh2JhqlLPj/yAtlJtki0m2IaJarNdkg/p6Jl/8VIL2Nq12saaN8IItEGv+TzcxqTUImxp0LsChdz4F2bSesL1Vs71AKjyJIvWiLRbStnjfIYK2VqIjR9VzdGICF0VUuo5oxftCQdPK1ZaRiSnGnnMF0tgAccTEFZsIHt5/SkFuGMsJc5YxHldE5BLgk7lNFwB/ALQAvwD0+u2/65z7qv/M/wRuBiLgN5xzXz/Vfpqkzd0QvkD36UMJkiK7lWPHCTvXIkXVmirdxyhs3cTYxR3U7j1BpbOZwsP7mb7qAkpH+qnsPzjDhFnYuGGGJpg41AQda6h0HdH8oZdfhOzcC5AWBq5yMtSxZXKSwqaNVLozzbGwYR3x4BCyZQPxrr1qlh0ZnVEAOGxuIhoYICiXcVGM1JRxE5Pe/Kmm2WhkND1+8FrwxISGg0xPZQ5BvrxT2KHhG8QxYf+odwIqam3H7VtVcz7Wg9TVarLtgUEmnns19TuP4wohP9j9z1ah3lj2mCA0njBEJASOADcAbwFGnHP/Z1afy4F/Ba4HNgD/CVzsnFvAnRKaCx3uaTUvVpNooCWUXGWawsYNuPpaZGRM6/t1dhB1H9c4vSjCPf0aRjeUCacctV+6h6C+Li2QK+Uy4cb1xHU1xA/vSvcVNjZqiMbkJIVtW9IMLYnwnU9bkkKB0Z95MnWfvzut6OAiXfcLmhuRcpm4f0DXD30ISFBT4zXbAgSBOqh455qgXNaEAR3txL0aykEUqZD0AjCoq1ONuFxWgQrZGmmpSDw1TWF9J65SYfSp26j5yl3aVltL5cmX0HttLe2PTFL41n0zPGzDzrVEx3v4ofumCUJj2WPRsMYTyfOAvc65gzJ/HbuXAZ9wzk0C+0VkDyoU71hwZBGktobJn7yCmm/c6wvrBkTHjmt+z5FR1WiO90IghK3tmlT7+/eRK12bpg0LW1s1FrHvJG5sbIZGFg0N6S6LJaIudbRRD9IFYuqCkMnnXUvTPd24pgZNbF0qIe0tyMQUbmgY55yGOTTUazLrMETq6xAR3PgEUltDPD6emkRdFGvO05MDGohfLuMqFY689Xo2vf9BmJggaGvVxNtTU5oFZ2ICnBdoxSJhqUR0vAcXO+puHydGBb00NxF89142nLwEunuIIF0HjZ7zJOKKI8g7ARnGMsbWCI0nktei2l7Cr4nIAyLyQRFp9ds2Aodzfbr8tjmIyC0icreI3D1VGSU6cZLy1+7K8nDGEVx3GdLSrELCF+sFdP2tpkxw9aW63jdLMMfDw7iJSfpecUWqQcVjYzOqs7vKNEF7W1o0dz4K69cRlIrUfPcR1UbL5bSIb7TvkIZUjIwS9Q8StjZDUjnDO/REJwe0msbwiFaclyCrFRhFGu8XuzTP6pZ/PQDT0yqoj/emRYYJ/M/da8xEkT4ggGbB8ZpwNDSkpmCASa1yIdddOiO/afD9B3Rd0wrzGisAE4TGE4KIlICfAf7Nb/oH4ELgWuAo8FdJ1yofryplnHO3Ouee4px7SqlYn4YyBGUtNiuFAu6eR6gcOKRCLHFOKRSQ5kb1tNx/BKanCVtaUu/SsKMDqa3FVaZp+5gK1rC5CXymmkLnWu3rHPHAoFZg8M4wVakpQ7GIbN6gJsuTA7jRUc1VWirqOqPX8piuaEB8EOoaYBQhxQLRwCDgU7cFokLNC8SgoT4tO+Wmp7zzS0zY3pqZUZubdE3w8oshjgga6rXChIt9hQxHuGaNxkQ+81qkVCJsbSXas5/Cuk7cvY+qEAbCb9+r53/a4giNlYEJQuOJ4kXAj51zxwGcc8edc5FzLgb+CTV/gmqAm3Of2wR0n2pwN11JvTdl++Y0li5sqPfCooHBN9ygWlMU4QaHcdMVFSwiahJtb9FqDL29xCMjmqh60wYKWzdBqahrilGEGxtXb8zGRoKmJqS2BikUCdpa5yS9VvPpUdzUlGaZ8Z6mXLwNNzXF1DOu0Li9thYVZpOTBI0NBDVlpLGRsKMdNzWFlMuq2QWiptLpCkFLM5WevrQwsK4jhsSDarqtHDtOkFTZ8MH10aOaWzUaHCJoaCAol6kcPaYPCc3qiRp8/35EhKi/H4Cpi9Yz9YInzahGH164VY+1ofa0LwTDON8wQWg8UbyOnFlURNbn2l4BPORffxF4rYiURWQ7sAP40WJ2kIQmRDt3q8CqVDQpduyIR0Zo/vhdqaCK+vuJRzUVWdTb551dJBUmUy94Mm50jMqBQ1QOHCI63pOaGOMdm3WdcHoaWhpxI6MQCNGx43OywwT1tQQXbEkrViRrdbKvi6C2luLt92vl+Z5ewtYWiJ1mbYki4t6+LEA+1vVDF2uibQJBasqEHe3Eo+MaJxiGmm3GZ9ORYkmD7CUgHhlFasqabNtX1dAJBpo1xzmivQdUS/aJAwCkXObElbWcuFzjGRPza9RarwJ4bPL0rwTDOM8wQWg87ohIHfB84LO5ze8RkQdF5AHgOcBvATjnHgY+BTwC/Dvwq6fyGE1QwaYpz4JSUTOzlEqE7W0UOtcixQJBfR1So8mqCxdsI7xwW1pfL9qzn6ClmbC5ibpHj6cxhBo7KKmZ0d3zsGpek5PI8ChBexvEPkPM7DWzMOTEDR2aFLuocxl/zlVq3hTR+MZCUTWznj6kWGD/v16l+UGTUk/1tT5LTaxFd6NItdmT/SpgK9N6XJddqOfAnwtXmfZOPDHBxdthelpDKcJQzcXDw+pIMzmp4SANDelrNzVF2LkWNzVFxz/ewaYP71JP2mKJwrbNBDsP+OOzW4ix/LHwCWNF0Fxe536i5b9pYdz+wdRcqE4k3nRYmzPjxbF6hfoYvegnryG8/f40bk8KqjmF69dROXJUc24msXj4nKTbN8GDuzS7TEszUV+fD9vIvEeD+npk2yY4coxoaESFcVMTUi6pR2tLM25ySk2RI6Nqmqyt0fW+XJhEPDmZBvOnBYWLBU2avXaNOgINDGq4xPgErjKdrjsmSbclUGEupZIed/LbD0LChnr1Wh0Z0bXHwKeASxKBB2F6HvVYNVzjR/ItBqM+85gxljUmCI0VgYgMA7tO2XFlsAboW+pJeLY65zqWehKGcTZYHKGxUtjlnHvKUk/iiUBE7l4tx2oYTwRm4DcMwzBWNSYIDcMwjFWNCUJjpXDrUk/gCWQ1HathPO6Ys4xhGIaxqjGN0DAMw1jVmCA0DMMwVjUmCI1ljYi8UER2icgeEXnnUs/nXCAiB3zWnftE5G6/rU1EviEiu/3/Vr9dRORv/fE/ICJPWtrZG8bywwShsWzxhX7/Dk3ofTnwOl/YdyXwHOfctbl4wXcC33TO7QC+6d+DHvsO/3cLWtXDMIzTwAShsZy5HtjjnNvnnJsCPoEW9l2JvAy4zb++DXh5bvuHnXIn0DIroblhGKfABKGxnFl0Ed9lhgP+Q0TuEZFb/LZO59xRAP9/rd++Us+BYTxhWIo1Yzmz6CK+y4xnOOe6RWQt8A0ReXSBviv1HBjGE4ZphMZy5oyK+J7vOOe6/f8e4HOoCfh4YvL0/3t89xV5DgzjicQEobGcuQvYISLbRaQEvBYt7LtsEZF6EWlMXgMvQIsWfxG4yXe7CfiCf/1F4E3ee/RpwGBiQjUMY3GYadRYtjjnKiLya8DXgRD4oC/su5zpBD4nWuC3AHzcOffvInIX8CkRuRk4BLza9/8q8GJgDzAGvOWJn7JhLG8sxZphGIaxqjHTqGEYhrGqMUFoGIZhrGpMEBqGYRirGhOEhmEYxqrGBKFhGIaxqjFBaBiGYaxqTBAahmEYqxoThIZhGMaqxgShYRiGsaoxQWgYhmGsakwQGoZhGKsaE4SGYRjGqsYE4TlARP5QRD66yL6/KyL/vED7ARH5qcX0PR/w5X8+JCL9IvIjEblRRHado7G/LSI/fy7GMs4dIvJsEela6nkYqxsR2SIiIyIS+vedIvIdERkWkb+afW9aaKxTCkIRKYvIB0TkoN/BvSLyolz780TkUREZE5FvicjWXNtrROQHvu3bs8a90R9E/s+JyCtzfX5LRI6JyKCIfFBEyrm2p/sb77CIPCAiz5w1/uv9nEdF5PMi0jar/bUistO37xWRG09xHp4tIvGs+X7pVOdvNs65P3POLermfjp9zwQReYk/h6MickJEPiYim3Ltb/bfyf+Y9bkuEXm2f/tM4PnAJufc9c657zrnLnm85lzlGEr+ou/y38l+EXlfrv2AiIzP+t7+31nuc95r/gzHExF5t/8OTojIe8TXYTqXiMi7/Pf5U7ltZf/bGvK/td8+h/sr+fPUNWv7tSJyjz9/94jItWe5n3P6ffgxbxWRXf43/+azGOdKEfm6iPSJyDkp9SMi60XkiyLS7b/Pbedi3PMdySkJAM65Q865Budc5DfdAvQBTc65tzHr3rTQ2IvRCAvAYeAngWbg99G6aNtEZA3wWb+tDbgb+GTuNK02swAAIABJREFUsyeBvwb+Yvag/obZkPwBLwFGgH/3B/3/Ae8EngdsAy4A/si3taEFSf8SaAHeA3xJRFp9+xXA+4E3ovXdxoC/T/YtIs8H3o3WbmsEngXsW8S56M7P2Tn30kV85rxERF4FfBz4G2ANcAUwCXwvOY+ek8A7RKRpnqG2Agecc6OP53wX4H8CT0GruDcCzwHundXnpbO+t187050t4po/E24BXg5cA1yN/hZ+8SzHnIGIXAi8CphdtPcPgR3o9/gc4O0i8sJztNv/AfTMmkcJLSr8UaAVuA34gt9+2jxO3wfA/cCvAD8+0wFEpABMA58Cbj4Hc0qI0fvkK0/VcZWxFXjEZbUFF39vcs6d9h/wAPol3AL8ILe9HhgHLp3V/+eBb59izA8BH8q9/zjwZ7n3zwOO+dcvAR6e9fnHgJv96z9DC5ombRcCU0Cjf/+DpO9pHPOzga552v4Q+Oh8/YADwE/N7uvfvxE4CJwA/tcp+r4p1/f3q/T9N/QGMww8CFyMCooe9GHmBb6v+HHePmueAVoN/Y/9+zcD3wO+BLwr16/LH+fNwAQQoQ8xfzT7+P0cf8dfM4PoTarGt7UCXwZ6gX7/elPus98Gfv4U38uXgbcu0J6eo9P4rv8B+HTu/buBb/rztqhrftZ4rwbumbXtbcDnc9fjLbm2m4E7TzHH24C3+dcbAQf8in9/EfoAI7n+X0ML+M44H8CR5Lrw7/838InF/Bb8MfSgwvUts/psB3YCL5p1PbzA7zM/t0PAC0+xzxcDj6DX9hHgd/z20/4+ctfWnwM/8tflF4C2Kv2+B7x5kdfNHwKfRn+DQ/lr138n7jSuwXcCe/3xPgK8okqfgv/ety1yzLf472QYffD/xVntLwPu83PfO993AjwduMuft7uAp886r3/ir+kR9N7RDnzMj3tXfr7ApcA3/PW6C3jNPPv8CPoAMO7HfTuqIDl/Hv4FfeiY8u2/yKx700Ln5rTXCEWkE73BPoxqEfcnbU4l716//XTGrEOfVm/LbZ4xtn/dKSLt6A1ptulIgCurfdY5txc9QReL2pOfAnSIyB5vUvt/IlJ7OnM+F4jI5ehN943ABvSC2bRA378H3gCsR7XzjbO6vRS9YFpRrejrqHDbCPwxqiUDXAJsQQVninMuBj6DmhPy/D7wW7PNy865DwC/BNzhVNN61zyH+hrghejN8WpUwOLn9iH0yW0LepGfrtnyTuC3ReRXROSqc2RSfBtwtTcN34gKppuc/iLP5Jr/IrBdRC7Lbfs59Lti9pj+9al+Q7ejAgnUWrPP/we1cHzXzxcReTUw5Zz7an4Ar/lvOIN9A6wjuwZvBv5uliXh/wK/i36nea4AHkjm5nlgEfv8AHrjbkR/5/+VG+9M70FvAv47eg4qwN8u4jOn4mWoMGxBb/5nyl7gRvQc/xHwURFZf5Zz60GViCZUKL5PRJ4EICLXAx9GtfgW9Bo6MHsAfw/4Cnqu2oH3Al/x9+WE16L3tI2oEnIH+jtvQwXxu/xY9agQ/DiwFngd8PfeojcD59wb0QemxLrznlntb0bP93t8+/tZ3L0JOE1nGREp+p3d5px7FGhAnwryDKImqtPhlaht9/bcttljJ68b0aeNDSLyOhEpishN6Amvm+ez+Xl1AkVU8N4IXAtcB/zeIua5QUQGcn+vWewBzsOrgC87577jnJtEBU68QN8vOee+55ybAv4AfRrK813n3NedcxVUyHUAf+GcmwY+AWwTkRbUFApzzWTJtjX5Dc65+4D/AN5x2keo/K1zrts5dxJ9QrzWj3vCOfcZ59yYc24Y+FOym/li+XNUY3sDahY74q+HPJ+f9b39wkIDOufGUEH1XvTp/tedc8k612lf8/67/aQfMzHdb0O12WpjDgINpxDqtwM3ikiA3rTeAzzDt/2kb0dEGlALyVurjNGQ29+ijiXHNGo5mPYCdgR9wEJEXgEUnHOfm2efZ3LPmAYuF5Em51y/cy4xWZ7NPegjzrmHvPD8feA1/kH5bLjDOfd551zsnJv9ELBonHP/5n8zsXPuk8Bu1Px/xjjnvuKc2+uU29HfdOIbcTPwQefcN/w+j/h7/Gx+GtjtnPuIc67inPtX4FH0ITzhQ34/g6glYq9z7j9z96XrfL+XoKbLD/mxfow+iL/qbI7zTFi0IPQ/uI+gmlWyxjKCPl3kaUJV79PhJuDDs54SZ4+dvB52zp1An7x+GziOahv/iZprTjWv5OL8v865o865PvSG9+JFzLPbOdeS+/vU4g5vXjagJksgfZo9sci+Y1X6Hs+9Hgf6XLaQnBx3A/rQAapZzmZ9rj3PHwC/LCLr5pnfQhzLvR7zc0BE6kTk/aJOTUPAd4CW07kZOeci59zfOeeegT7J/inwwVna18tnfW//tIhxf4RqWYKu8SSc6TV/G/B6L9zeCHzKC8hqYzYBI7N+D7Pnt9d/7lr0ZvZloFtELiEnCFFt4iPOuf1VhhnJ7e90jgXghL+xJYyhwrseFcq/Ps/nzvT8vRL9jR4UkdtF5CfOcjzI/Z7QpYIisx4Cz4DDp+5yakTkTSJyX/LwhmrBZzU3EXmRiNwpIif9mC/OjbkZ1UJPxQb0XOU5yEzr1Oz70Oz3yQPYVuCG/EMq+kB7JveYs2JRgtD/eD+AalOv9BoGqHn0mly/elQze3ixExCRzaiJ58OzmmaM7V8f90IQ59ztzrmnOufa0BvLJai9v9q8LgDKwGPOuX5UYJ4TD64qjJJppvibesc8fY+iF2DStw41N8zXN+/RWbtA31OxCz0Hr85v9A87r0TXw2bgnw4/i5q7zhVvQ7+3G5xzTahmA3PN3ovCOTfunPs7dL3x8rOZmIj8KnrNdKPrEQlndM075+5EHyJvBF5PZhadM6Z/vZjf0O3o03PJOXfEv38Tahq/z/d5HvAboh6hx9Dr7VMi8g7/Wzh6hvuejx2otvtdv7/PAuv9/rf5sa+epe1efap9Oufucs69DDWhfZ7s4eRs7kGbc6+3oFpntYfA0+Gs7yuiXq//hCoc7c65FnTt/ozN/qIe958B/g/Q6cf8am7Mw+h5OxXdqADLswVdtz1dDgO3z3pIbXDO/fI8/R+ve/aiNcJ/AC5D7bN5df9zwJUi8koRqUG1hgcSlVpEQr+9AAQiUuPNq3neiC52z34a+TBws4hc7tcefg9dEMWPfZ03izahX26Xc+7rvvljwEtFQzTq0fWxz3rzG6i9+tdFZK0f+61kZqqz5TGgRkR+2h/r76E31Gp8GniJiDzTe839MfN/J5/2x/R03/ePOHOB4VAHlt8TDTOp9ZreP6NP0++b56N/hK4ttJzJfqvQiD4hDvi1hwXt+NUQkbeKhrbUikjBm0Ubmes5ejpjXowu+P8cen2+XTIX/wWv+VPwYXQNtOKc+96s7b8tIhtFZAP6gPAvixjvdvRm+R3//tuoJva9nCXgeag2ca3/60YdCf4ut+/fE5FWEbkU+IVF7ns+HkIFTLK/n0c1gmvRG9+3UQeG3xAN3UisS/81dyhFNAzjDSLS7B/Ch/wYcHbfx8/5+0sd+tv7dHLe/D5r0N9Y0d+7zsSnQvw4Jf++RnJhYPNQj970e/1n3kLm/5CMW0N2Xyn79wtR8v17gYpoCNwLcu0fAN4iGooS+Gvx0irjfBX1tXi9/739LPrQeSb3zy/7sd7o7+VFEXnqLGtOnuNo9MA5ZzFxhFvRH861wDHJYrHe4JzrRTWIP0Wfwm9AF0oT3oje6P4BfRIeR5908ryJmU4yADjn/h01sXwLVb0PMvNG+Xb06e0was57Re6zD6MLpR9DF4gbUVfohP+Nei89hi7e3uuP4azxdvFfQYXKEVRDrBp87Of5q+hi8VH0HC7U99fRtb6jqOmnBw15OJN5fhL9fn4LPY+PALXAMxKtu8pn9qOaTP2Z7LMKf+332Yc6vfz7GYwxDvwVan7tQ8/nK51z+XCYL8nMOMJqa1dA6vL+UeDdzrn7nXO7US34IyJSXsQ1vxAfQW9oH5m1/f3o2umDqCD5Cplj00Lcjl7biSD8HmqNSN4n67DHkj9UgPQ75xKz6LtQk9hBP95f+t/eGeHXevL7OwnE/n3kdH375ejvfgB1Vnm5374QbwQOiJrQfwm/3noOvo9/Qa+dGuA3cm3/gV5bTwdu9a+fxemz1X820VDHUYvMvDjnHkGv6TvQm/9VwPdndUu8J0HX6BZcj/RKwG+gmnQ/apX4Yq79R3gHGnSN9XY/d0TkH0XkH32/E+ja3tvQpZm3Ay/xS0ynhZ/TC9Dvqxv9Ht6NF/CiCUW+lvvIn6MPbQMi8junu7+FkAWWIYzzGFEniAFgxzzrP8Z5hjdn9wBP8gLWWCJEE3x81Dl3XmduMp4YLMXaMkJEXirqYFKPmoMfpIqLs3He8svAXSYEDeP8wgRhDq+Kz077NjJLPV9KXoaaELpRp4TXLuRZuFLwpplq38s/nmdjVhtvxK9VHwB+EzUpLXa8N8wz3tk4syx230/4b0FEHp5nn284w/Hm/T7OYo5fm2fM03IikyxPZrW/LWcxv3N+zKsBM40a5yWiab7+BgiBf3bOzUnTZxiGcS4wQWicd4iGnDyGZrjpQh2bXuedCAzDMM4phaWegGFU4XpgT+L5KSKfQM3C8wrCkpRdzTzOrFKjXubTDdnlHnjn+3Ak56wY6UYX++Q+uVC39FWyLchFrsRu5udy/ZJeMx43F/Pwmd+3f60Z8GYNluzHzydqyjIFFtaqQ/F0nOUnkF59nRy/TE6TNfoxavU8xaXcHPyuw/HsGCWKGZ8eZKoyds4rZRjGE4kJQuN8ZCMzM3R0oW7xMxCRW9Cky9RQxw3yPAjmJqUJL7oIgJ6fyFKl1gzoDb35zixaJR4c0v/DGm4qxawggoR+OT3U8RPhCuAmVODE4xO5/uGMz+UtL246n5AFiCNmI+Vs/GQsN6VC28VuTltQq2Fkgy/I8gi0/+oBAHpGG9JtpffrOSgN6RzK+3qznfqxhq9eC8DQ5uz2UBzVfbY9NJJ1H53ijj0fmDN3w1humCA0zkeqaRhz1Cjn3K1ojBdNQbuTcplwTZZsp3JUMzvFNXqZ157MazM63Mh1WWao+n0qMMKekwBEJ05m+/LaIrEXdjmhlwiyoCaLaY4ntF1C3eYms3DPVMB6DU9qM002EaZSKMzZlhDU12VvplWji0a00szI+sz/rf8bGntc05cTwj4ZWTih+3ZDmWCLx8YAqDvcDUDtdE5bTh4wckJb1nWmWrRhLGfMa9Q4H+liZvqrTainrGEYxjnHBKFxPnIXsENEtoumk3stuSwYhmEY5xIzjRrnHc65imgOyq+j4RMf9CnmFvoQbnKSqCfL9CRP1vWywR1qehxdlz33tT+iJkWXN8J2azH1ZI0wT+KMkvirJGtzAEG7rrvlTamJKTGeyjmjJFOdnplNLL/mF/i1x3g0K6qdmFKT9cZ4dCxrS9YIS5rCt/lgtv7o/Jz7L87muv4/ume0JeZQndestcs8VdYx3dg4xPNVDTOM5YMJQuO8xGmNu6+esqNhGMZZYoLQWBFIEBDU1hFPZE4pcVk1oel61X6aD2RajfPKoQszldBt1fKMYZ86o1S6szKKrjJTW8prdVFvL7NJtLjZ2t8MqjigJOSdftzGtf6F1xwfyjK0JeOHG3RJta4r0/AmOjSUonwy04RPPk1LvTU/pk4y4fTatK3inWQS7TKuMvW8JuymptLQEcNYztgaoWEYhrGqMY3QWBE453BRROXZ16bbykcG9cWlqhkNbcm0mfZHVAvLKzRB74COVVlgrSztnAtST+L8ctpfEm6RhFTMWCv0GmC6vtfanI3VqCEccW6tM+jpB2D8Kq3LXDuSOdS6fj3GNAyiIxvr+PWq2XX+KNt3aVDnOLZJtd7GvqG0bfKFTwKgOKrHXzqQzSHRFtNFUsBFj2OlVMN4AjGN0DAMw1jVmCA0DMMwVjVmGjVWBFIuEVy4leJgloXl5FPWADDdoA4xI9sys15tn176rQ8OptucD1mIR32x77wTizeFpmEUOfOpq+Lsknw2nqjSlox1haZ+k4FcdpfeE/6AMieeyUs36MemdP5Rcy6n6pGj2t1nogkOHE2btn1O57jn51rTbevu1Gffhn0aIhKtaUrbikNqQp1qVUefaEfmSFOb5CH12Xp0ioLMjQ4xjGWHaYSGYRjGqsY0QmNlIIIrhkS1xXRTEjbR0K2a1Np7stCK8U6v9dRnibULSRUJ7+iST7qdsGA4xCJJc4Xu1bzirpTNOQmRCJoa002lHtVUo0YNtg9zGmTFB8EnuUzzT7aVZnUSWndHpgnXH9axZES13srGlrRNKrrvhoc1sYArZ/NKNNV8zlQKBZy5yxgrANMIDcMwjFWNCULDMAxjVWOmUWNlMD2NHOkh3Lwu3dRwNCl3pP8GLsrKJJV8fb3w/j3ptijJ7+mdWVyuxFA+o8ppkTjGFOf+1CpPUmeZwmCuRNO4rzmYq1945AWaZabip7/lb/dm/f24btIXFc7FK4Z3anrWhty+k7yoA0/T8lONezIza9Sg56uyVh1oCgdyjjEdOod8SSgXOwskNFYEphEahmEYqxrTCI0VgYti4sFhokuzrCuJ80fdXq0K0fBwFvIQdWmYQVzN+SXJ/JKvEp+GM1TJIbpAztAk3CKfwka2aDhEkrnF1dembZU2DY2Ybs4cdVp267zrd/uCwbnKFGGLzyTj5xpfmB1/eFKzxsQnB7IJFXSuLXd06Vjr29KmsXW6z8YDmq/URbnKEskYVUJKDGO5YxqhYRiGsaoxjdBYGfhco8EPHkw31TWodpWsa51u6MOM+nxVAulTEi0pryH5nJxJ/8LWTFNjyIdDrNNA96icq2045bXRSi7k4ZCvjzigGl5QV5eN5TXBwWdtB2CsIxtr0ucwbd6/Md3Weo9qoVGntsXFrH847UM39h3RDYXs9pCcQ8ltc1Z5wlghmEZoLBki8kER6RGRh3Lb2kTkGyKy2/9vXWgMwzCMs8UEobGU/Avwwlnb3gl80zm3A/imf28YhvG4YaZRY8lwzn1HRLbN2vwy4Nn+9W3At4F3LGrAWc4q0dDQzPa86TKukj1mljlzxly9GTBvGkzbfP98iIUU1WQ5fcOl+n4oC5EIe3zpJO9HU+wfT9sqTRojURjLOfY8quESYZvPAnPxtrRtqlX7d78omnFcAIUTmhlmrDN73m31Jl6Z1n7FPUeyeQ3ovKLE5FnNMSa3LWxqQkbsWdpY/thVbJxvdDrnjgL4/2vn6ygit4jI3SJy9zST83UzDMNYENMIjWWLc+5W4FaAJmnznh5ztb553zNPGES1/kmQfRUHkbT47kQu2NwH45f3++K2QfbM6aZ0n+GgaoIymc0hLGg/2Z9paskspFbDLI48O8sPOt6p8/lv1/wQgM/ffn22n1DbhnZkx7HxX3Q+4rVYN50F4M8+tqA+q3IR+5CNfEhJNDSEyxXqNYzlimmExvnGcRFZD+D/9yzxfAzDWOGYRmicb3wRuAn4C///C6fz4aA2S6OWpgOrVi/Qk9dwksoKyTrgjLXCBcaIvUKXD2tItL7p9er0WjzUm32gxdcA9MHtw5dmaeHqv6EOtNFELu2aX3uMWxt0rJFMcwv9IX7tc08DoJyVMWTt3artffsD/5Ru++n3/IyO5esK5qtJJOuliTYb5wL303Myu/90boeGsUwxjdBYMkTkX4E7gEtEpEtEbkYF4PNFZDfwfP/eMAzjccM0QmPJcM69bp6m5z2hEzEMY1VjgtBYGYia7/LmvNkkTi2QVXeYUWg2cYjxJtG8qTMNm/BmymgwF5qRhGJs6Mz692soQrj/mHZZ157tO/RZanxWl9ov/CgbatZcAAKfIUdOaoaZ9oeykI9gQud69FnqQHPrb/1N2vbFVz4JgOf/7FvSbcWTvnJFNDcbTrKfGceWtPlzEY1k59dFUVpI2DCWM2YaNQzDMFY1phEaKwRRbS1f5SEJHvcOMfmA93h2sD3McYjJ5xqVOl8hInGgqVZpYjzTLseffAEAUUmfNcsncgH14+rEMnSRalktjzanbZEPas+Pn2howUbVOCfbMgef6QYdY3ytHvct7/vNtC3wuUPXPZbVL0w14SQH6qYsD2ml+9icY0rnkJyvnAZZWLsG6bNbiLH8MY3QMAzDWNWYIDQMwzBWNWbXMFYEghbPdTmTYpL0JM0TGp6ikKw3+4VJ+aaxsbQp6jtRta921H3mC+wmRYFLYxpP6Eq5zDITarJt+eZu/36uw86MnJ7NGnco/WqerOs9me3HFwxu+YH+lOOT/dm0/LjTN16dbfOZa4qD2jbRlJlZa7y5NB5SpxwpFbM5T03PGBMgOtGPq8wfX2kYywXTCA3DMIxVjWmExorAOUc8NV21LXGaoUpezBmFZn1IQeIYMqMAbtKnirNMkpFFRrMqEqU+1aYm12k2mHA8p+ENjs4Yq9q8Z8xrXMeNk9CF3HFIyWeD8Y49wdbM+aXgtz36skzr2/4l1VBPXtkIQNv9Oachn5VHkpypeQ061H0WNmRZcOL+AaRimWWM5Y9phIZhGMaqxjRCY+UQRzOC5pNqENXqCyYhFdVybeI1yCRfqDae+pnR5XNzNuk6Y83+k3P6xcc0j3h+DTJlVlA/5PJ8RnNrDgbtbdqnWTW8fW/Igvqb9+o65eb/zMYqHdHwjDX39PrPNWXzT9b/Ojt0N03ZmmcwqHN1J7I1SKmrgwl7ljaWP3YVG4ZhGKsaE4SGYRjGqsZMo8bKIQhnFMed05xzfsmHASQkpkdx3jSaN0+mZtMqIRiJ80oxywEqY94s67PCREMjWf94Vp7PfJaa2W0w95hybcNP3w5AfZeaLhsOZZl1Gg/NPcaJLZqTtDbJgpM3G/sixYnTjxzsytq2b9ZtpewYpy5ch7s3C7EwjOWKaYSGYRjGqsY0QmPJEJHNwIeBdWjhhVudc38jIm3AJ4FtwAHgNc65/vnGSZmV/7OwzjuOJIVt+wfm9q0WGD83ygI3PTW3f9LmtSqX09yi3t4Z/YNccHo8Ec3YlryHXNhEzjlH6uv8+KrFhRuzEIZwSifb+yQN06jrySYfTOm4QxdkTi9tX3tM5+dDMYJcYWJpbZ6xH/GJBQA4qU42zjvlAJT2HkMmq4esGMZywjRCYympAG9zzl0GPA34VRG5HHgn8E3n3A7gm/69YRjG44JphMaS4Zw7Chz1r4dFZCewEXgZ8Gzf7Tbg28A7FhpLamsILr6U6fZsHbB/i2o7zQdUUys+kAW8k6wRVqkikVJFW5z3PRAPD887VLW1y7TuYWOmZcWjSUhFLlWc3xZetA2Aqc4s5GFsjf6E2x/yQfe5VG6916pG17I3Fwbi082ltQd93USAsFiY0WdGaMnmDfriZNY/Ghi0FGvGisA0QuO8QES2AdcBPwQ6vZBMhOXapZuZYRgrHROExpIjIg3AZ4C3OueqFAqc93O3iMjdInL3VGX+yvSGYRgLYaZRY0kRkSIqBD/mnPus33xcRNY7546KyHqgp9pnnXO3ArcClLdtco/+eiPFE9klfcFnNGTh5OXqSNKxJ3MaKTTptrg3qyqRZHpJwixmZH6pFuqwCKo5xKTzT0yQo2NV+mem1CR0Iwl9mG7MjrFmUMctHpmbwaadNTpWmMsHmuRd9aEihbVrsvlMehNqqM/HUp85y4jPW1o5MWs/zmEYyx3TCI0lQ7SG0AeAnc659+aavgjc5F/fBHzhiZ6bYRirB9MIjaXkGcAbgQdF5D6/7XeBvwA+JSI3A4eAV59qoPraSW64fC+7P3xJui0JKRjepu/bPt4353NpWESOqjlAE6qET4h3Msk7l6RjVQvwT0Iqasq+z+SC/cVraDV37wFg6FWXp22d/3lE9+01yXwFjOl6nVfDPYfSbYlGl+RalVzliyDJO5popcd70zaXnJPZVTdiqz5hLH9MEBpLhnPue2hN3Wo874mci2EYqxczjRqGYRirGtMIjRXB2FAN9/7XJbQPZs4bw5v1Oe+iD3sTX21WoinJrLJYFjJ/ukmfKSbnXBKPq4kzbFMHl6jvRO4Dmv2lWr7TJLNMPi9q5QrNJ9p/qW5r2ZN9LjGJEuuYUxdkZZjqutRZKCnsqwOHM47D5U294uMgg8DPL+ewk2a8yRR4Nz1lzjLGisA0QsMwDGNVYxqhsSKQGIqjwui6TGNZd4d38DimTjLRUBaimBa7rZJXdCGSwr8uyuX0rK/142cVJpIwCDfmKznkK1OEczWu9HNNmu8zKY4LUOxWB5eGZh1jqiX72ZYPqLNL1K7ZaUpHspSscYPONc450Mw5jpxGFyeFhRfIp2oYKxHTCA3DMIxVjWmExsrAgVRgw7dyRSp2HwQgGp+rESXB7NXyiabrdDu2Z20+oJxhX7WhJqvaQBKCkNMIxbdHA5qbc0YtxCQUoVoohtfU4u7j6bbosm0ATLT6sItKpsWNXKWVKOof9aEhY9mxypBf8wuz590kXMRFfq0wqpIg4DSTBhjGcsc0QsMwDGNVY4LQMAzDWNWYadRYERRHY9bfMYZ7dF+6LS06W8XFPw2HmM6cQNJsK76QrytlP4/KGs1SU/ROMFMdWahE4Yc7AQh9YVvdODMMQnKFeZmYaRKVIHPwqfSoiTPsaM+G6lUnn2bvoDO5JsuZWvvD3frC5wx1jbliuj6PatUSUFUy6qTz8Y49kjOpVs2QYxgrBNMIDcMwjFWNaYTGikAmpins6SbKaXixd4hJNJy8FpRqgvncmeLzddapxtVzfUvaFpVVa2t9zDusTGfhE7J1k744miuS4fedOsZUSV+azCvwlTAAnA/Ed7mA/8qF6wEIR72jS06DlFadYxoG0Zc5C7nEiSfnlJNqu4mTTBXHmNShprJwHlEpl2HSco0ayx/TCA3DMIxVjWmExorAVSpEvSfStT8g1cp5CLk1AAAIUElEQVSqskCIgFunwexTTZm2M9GhY41uUK1x7T05jbCiwezFQqZ5RTv3zDt+qgn69Gvxhix4PhhW1dHVZ+ngXFGfVweu0DXIlocGsjYfIpHUNMyv5QWNOq/8sbpZx50P9J+zbijZc7Ikx5bTRt3kpKVYM1YEphEahmEYqxoThIZhGMaqxkyjxpIhIjXAd4Ayei1+2jn3LhHZDnwCaAN+DLzROTe/v39CHKWVIGbspzD3Mk9ybeZNicnr0DuZFIczs1+tL2DRc6M62Zy4KguHKA6pc40rZGENW4bU2SU6pg40+QwuiUmUNjV1hoOZY8zk1jb9XDlnZvWOOrV9um85mqtk4Z1+0uPIZ6uZzoruZjv37d5EWjWMYlYfmGtSNYyVhGmExlIyCTzXOXcNcC3wQhF5GvBu4H3OuR1AP3DzEs7RMIwVjmmExpLh1Oc/SdBZ9H8OeC7wer/9NuAPgX9YzJjVnD+qVU6oqkElWo+v7dd0MNOoJtr1p1K/VzXBcn+mLRZH9HX797vTbZXDXTPGD5ub0rapyzTcIgmDKPVkGmFhQOsElg5kY7FOg+XpUycZyQXNV06cnHNsCVUrRiyk2SW1Bhfok9eurSKFsVIwjdBYUkQkFJH7gB7gG8BeYMA5l9xlu4CN83z2FhG5W0TunmZukVvDMIzFYILQWFKcc5Fz7lpgE3A9cFm1bvN89lbn3FOcc08pUq7WxTAM45SYadQ4L3DODYjIt4GnAS0iUvBa4Sage8EP58dZwPkjn9MzLcNUxQwY7TsEQG1uW2FMHVwkVoFb25uZTUu9atqMj2WZZZIco/GEaqpSn5VhKozqZ/sv1Ywy9YVsXiPr1fS6ZjhLRSNjasZNstREvb1zjzFhhoPL3OY0y061zDKLiAk0c6ixEjGN0FgyRKRDRFr861rgp4CdwLeAV/luNwFfWJoZGoaxGjCN0FhK1gO3iUiIPpR9yjn3ZRF5BPiEiPwJcC/wgVOOJOrIkVSQAIhHvRNKEipA5hgT+GoQ8USVcAufj5Mwl6OzoupVuV81ovI9WeaYKFeQN9u5jltYr4Vzp7dm2WMm1qhW1rJLP1c4muUHLQ606ov+wXRbJXl9DkIY5mjM1ZyFqoRPJOfVTdparLHyMEFoLBnOuQeA66ps34euFxqGYTzumCA0VgZO1/3caBaKkLj6p+taOQ0n0QSrhQMkWlO0K9P6xGtJNTWqGcW5KhfJuHltNOxQTRCveSaVIwAaHtSaifElW3U/61qzse5+ROdQm+UaXZQmWEWLO+P+1SpSmCZorGBsjdAwDMNY1ZggNAzDMFY1Zho1VgQiQlAq4qJ8EVp9zktCJfLhE2mfRWZKSco7pYV2c04m1fKWVrqPAVDYsG7uWD6UIth3RN/X5Eou+SK90UDmLLMos+fpOtJY7lDDSDGN0DAMw1jVmEZorAxEIAxxOa3MzSq+UC3AfIYWOCvwfkabD0BPAtITbRMyDW/0p69JtzX/4CAAU9vXAlDaneUEmLxUM8YVB3SucjzLF+omqjileO2tmuaZTj0J4B8bm9O2YPHdhagWWmEYKxDTCA3DMIxVjWmExorAxXEWQJ+Q12jmIVn7A1KtL9EEg8bGrM3X9ktTpa1pS5sGr9HqEI37ssB6V68J2hJNMF6bhUiUerRfXO9DMfoH0rZ4qkoNwaStiiaYHGO1tcskMUC+FuJimBN2YhgrHNMIDcMwjFWNCULDMAxjVWOmUWPlEs/MHpNWnMhRNWNKYm4cHk43pU4yienxRJYftOGAr1MhWXjG9HqtVlF8TM210p2rGCE+rOMRvy0XwjEn32c18o4rs51YZlSfWMRYhmGYRmgYhmGsbkwjNFYGIkixNCM8IA03SBxQ8tUUUi2Ruf29U0rYnjnE4HOLxj7QPdyYBcpPtuvnarszZ5nigGqCbnBoxpiQhTokmlreKaWq9jpL66uaBKCa1hfPrTm4GEcYc5IxVhumERqGYRirGhOEhmEYxqrGTKPGysC5OVlTEnNkNTNoNZNiYkJN4gfd6Nic/olZM2rNYgwrdfo8KVP5LDX+GdPH8oVNTWlTNDKzYPCMw1iEWbJqNpwkZvAUmWPSArsL5VW1OEJjlWEaoWEYhrGqEefmupQbxnJDRHqBUaBvqedyhqxhec59q3OuY6knYRhngwlCY8UgInc7556y1PM4E5bz3A1juWOmUcMwDGNVY4LQMAzDWNWYIDRWErcu9QTOguU8d8NY1tgaoWEYhrGqMY3QMAzDWNWYIDQMwzBWNSYIjWWPiLxQRHaJyB4ReedSz+dUiMhmEfmWiOwUkYdF5Df99jYR+YaI7Pb/W081lmEYZ4+tERrLGhEJgceA5wNdwF3A65xzjyzpxBZARNYD651zPxaRRuAe4OXAm4GTzrm/8AK91Tn3jiWcqmGsCkwjNJY71wN7nHP7nHNTwCeAly3xnBbEOXfUOfdj/3oY2AlsROd9m+92GyocDcN4nDFBaCx3NgKHc++7/LZlgYhsA64Dfgh0OueOggpLYO3SzcwwVg8mCI3ljlTZtizs/SLSAHwGeKtzbmip52MYqxUThMZypwvYnHu/CeheorksGhEpokLwY865z/rNx/36YbKO2LNU8zOM1YQJQmO5cxewQ0S2i0gJeC3wxSWe04KIiAAfAHY6596ba/oicJN/fRPwhSd6boaxGjGvUWPZIyIvBv4aCIEPOuf+dImntCAi8kzgu8CDQFIu+HfRdcJPAVuAQ8CrnXMnl2SShrGKMEFoGIZhrGrMNGoYhmGsakwQGoZhGKsaE4SGYRjGqsYEoWEYhrGqMUFoGIZhrGpMEBqGYRirGhOEhmEYxqrm/wc1fxhCozEkWwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, axs = plt.subplots(nrows=2)\n",
"axs[0].imshow(imcac.get_img_by_metal('Ir191'))\n",
"axs[0].set_title(imcac.image_ID)\n",
"axs[1].imshow(imc_cropped.get_img_by_metal('Ir191'))\n",
"axs[1].set_title(imc_cropped.image_ID)\n",
"plt.tight_layout()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Save the cropped ome.tiff"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"fn_out = os.path.join(fol_out, imc_cropped.image_ID)\n",
"imwriter = imc_cropped.get_image_writer(fn_out)\n",
"imwriter.save_image(mode='ome')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Load and show the cropped image"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"ome = ometiffparser.OmetiffParser(fn_out)\n",
"imcac_cropped_loaded = ome.get_imc_acquisition()"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7f2a73e945c0>"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJztnXuMXOd53p93Zmcv3CV3ucvLLimK1M2m5ItoR5XVpglUOTYUJ4ic1mnjNoVaGHAKxIXTBkXU/JMLGsApkqgtWqSVYVUykMQ27DhSC7uNoshNU8uyJVuWKFEWL6IkXsRdknvn3mbO2z92aHP3fT5xzs7OkMvz/ABidz9+c853Lu+cmed7v+c1d4cQoniUrvQAhBBXBgW/EAVFwS9EQVHwC1FQFPxCFBQFvxAFRcEvREFR8AtRUJoKfjO718x+YGZHzOyB9RqUEKL12Foz/MysDOBVAB8CcALAdwB83N1fTr2m07q8G71r2h8AWHcXbV/q6whtpRrfRnlmMTbWeGfPMjIIi018V7QvSoneWbwOdP85x0CvbquyOtm42DkA4E6OLTUstt3Eeaxt6QltHTsWQttSVua7GovtqXvJFpZIY2JcPfEezTp5XyOnpjzH7wWrrWyfW5rEYvVC8pa8lDiixrkTwBF3PwYAZvYFAPcBSAZ/N3rxAfvgysYSvwiM8s030/bRvz0Y9zXBT1b/t06EtmxyivbNpqdDm1U6Y1s58QGqHI8t9Qbm8/EGzebmaV9j202Mgb25+1KV9qVkiTufjasrHhsbKwD4YnwTdvIGmNpGqaeb9p388G2hbehXjoe20dk++vrO/xbvpc4pfr66jo3FxsTxTr93R2ib2sPDrzIbz8PgwRnatzy78jw+feRztB+jmY/9uwG8ecnfJ+ptQogNQDNPfvbRIrxlmdknAXwSALqxqYndCSHWk2ae/CcA7Lnk7+sAnFrdyd0fcvc73P2OCvhHXiFE+2nmyf8dALeY2Q0ATgL4RQD/+G1fYRa+F5a3DdGu1dNnQlvWzYfbcz5+v7ca//448774zaT3GP/+Vx49H9pq52KbJwRDZOQ7e+J7PPtuXerm32uz+bgNK/O+vhC1BKZbLHcm57EnCrRJLaKDiFqp4yWUehOfDJeisFabmaVdZ0bi82z8iRtDW/dZfn/4lthWnuf6kU/F7+HZhQu076Y3w3MRPUtEfAa4DpbQXmx458qG1L1IWHPwu3vVzD4F4H8DKAN42N1fWuv2hBDtpZknP9z9awC+tk5jEUK0EWX4CVFQFPxCFBQFvxAFpanv/LlxD+pzbfQs7Wo/FjO1Jm/hqcGzw/E9bOhlknoJwFl2wqlR2pdl+DFSqaYsgzWV8VYaipllbGZhuXPcRrbIj5eOK6Uys74k666UyFLMZqMCn5pZYBmJ2SxXymmGX2eF9u1/PWbjObk+4+/g12HkL6Iqz14PcGU/V/ZkihxZlX5hbtVrEynhBD35hSgoCn4hCoqCX4iCouAXoqC0VfCzUgmlnpUpnBlZygoAWVcUZJZ6ufDSfzwKJJ54W/Ny3IbvHaF9y2djumn11Fvx9dXGRZ6U2FYbI8tDEzARLY+IlyRHWimDpWr77riUdfk/SHrtwcO8Kzm28q49pCew6UQU4ea3xzX+Xef5DXL+ruHQ1v9qYjntUjy2KknjBbhAmeW4ZA0vjU4si6Zjanz3QohrCQW/EAVFwS9EQVHwC1FQFPxCFJS2qv3uHowvqncfoH27Tk7Gxv1RtQWAqeujEjr0MlepmRhaGpugffOo+A2TMCxlam5KwWfmIUnjD5b2mzKGYGm0W/tjv83c/CQjqdql0XHad+4914W2nhmu4Pt4vBeYkQYAYHsc75k7o9K+89s8HbpzMp7zC9dxk5HNZ6Px68K976d9K7PxXuo8zlPb6YwByxUH4KsuZR5fZj35hSgoCn4hCoqCX4iCouAXoqA0JfiZ2XEA0wBqAKrufsfb9u/qROmmvSvaKpPc3fX8HdtC21IfT++d2RfFkJ6z/NC2vkjEI7IOHQCy2TnSSMSylIhH1oGnRETPkUbLxpDN53h9arzvihWRbII41I6d49slpaoW9u/iQ1iM16zWnyjldvJ03BVxCgaA0vHYd99X4zk/8ktb6euHvxWfh33HuK9DbVu0+q1McSFxcWtMya7dwlOfe8h5rBE3ayCWQ7PGbR3WRe3/e+7OZUshxFWLPvYLUVCaDX4H8Bdm9ly9LFfAzD5pZs+a2bOLNW7TJIRoP81+7P9xdz9lZjsAPGFmr7j7X1/awd0fAvAQAPT3jLSoNrQQIi9NPfnd/VT95yiAr2K5bLcQYgOw5ie/mfUCKLn7dP33DwP4ncu8CF5ZqTTXergLKzPu6DvFUxx3PBcNQeZ2ctfYWm9s70gYILA02mSdO/b69TDYaAHJmnhH3wxNzlxymREHgNKWzaGtc5TPpNQ2RwfgMplZAIAqccRlNQgB/jSr9se08OGn+b3U+yZxIJ4hsz4AqrsHYt8qPzd9L0WHaO/i9z6bTUkdL1bNeniOBN9mPvbvBPDV+lRDB4A/cff/1cT2hBBtpJlCnccA3L6OYxFCtBFN9QlRUBT8QhSU9pbrWlqCnVwpfJT3RLdUAOg7TYS1hJYxcXNcy945yzuXv38ktNUS6b0sDZaKgAln1baSStmtNH6Jq++P6b0dk1FosrmEzwARAk9+ODr6AkCV2A9c/5+O0r7sGHyBpzM78S8of+ul0NaXOC+sbNrEXbtp381HokBZ6+OCcHVHTAXuOJ5I2d0ez1k2x9PgQzm1HJPpevILUVAU/EIUFAW/EAVFwS9EQVHwC1FQ2uveW8uQTa40Rqjt546tLE1y09HztG/fSzH9s3YimjoAQJYn5ZaYZlhXTEtdbajwI3LU1GuyTh4zDlneRjyPdj032GBust4bU2Org9x0Y6k/Hu/AYW5e0ns4XsvUrEt5IDryglwHAMhuivdT+Xx02c3Oc8dmdMTrMPD0Cdq1NhJnBi4Mc7V/8/G4otVrPMUYbGypeyExy9MIevILUVAU/EIUFAW/EAVFwS9EQWlvei8p11X65ou066a+KColUxzbuG7eydpyJMS2XOW+crgCs9JNqX117CWC6lRijf1wdLStdZESXouJcl/VOK7eN7jzLSaiCFfalPAZIOLe5E/eQLte2B7Hu0BKjvW/xlN2tz4XRc/aTiI4AsgqcV/lpYTXwbGTsTHhQMzu85RbcUjvzYGe/EIUFAW/EAVFwS9EQVHwC1FQLiv4mdnDAH4WwKi7v7veNgjgiwD2ATgO4B+6Oy/EvpoGs9ZqU1EQSpIjO44acKZqnzco2KVEl5RI0+i+Uj4BVokC2NIH9vO+U3E9fnk0liwDACe6ZWU8mldWt5DF+AA6LpBMy1f4Gv3yYDS/xDv20b6LW+P+Tv104j4i173jXDTKvLCTP/e2EvHWlvi+KkeiiFee4Oe2xu6RPFl7ib7lLSt9Amym8ed5Iz0fAXDvqrYHADzp7rcAeLL+txBiA3HZ4K8X4VidiH0fgEfrvz8K4KPrPC4hRItZ63f+ne5+GgDqP3m5Uaws17WEhPe4EKLttFzwc/eH3P0Od7+jAr4SSwjRftYa/GfMbAQA6j9jORIhxFXNWtN7HwdwP4DP1H8+tuYRpFJYc6xlz9OXpgLnWROdIw03T+plqTsq2tl8Ip2ZOAh3vRbTUpc3HN/ffZGnQ5cno7JvC7FvuYM/M+y1qH6nroz1RJ+Ak3eTGQAAczvjefz7tz9D+/75/4nlIr0cXz91Cx/Z7kfiebTErI8vRafgXNe8l/siZMTXgPlIAHFWzBMzV3T/l+tgZn8K4GkA7zSzE2b2CSwH/YfM7DCAD9X/FkJsIC775Hf3jyf+64PrPBYhRBtRhp8QBUXBL0RBae96fkKph6eK0rX7eUTABEw4SdU+Z+m5NOV3HcaVMR0ysb6dCXZLI3EtPgBU3hiLjQOxdBQAal45vT+WU+t94iB9eW2elPZKpChnW/tCW2UmUWKN3Apf/+pdtG8XSVHe8WwU5r7xuc/S1//Mv/+50Jad5mW12H1D08fBRVom7AGJ+y51j67e31LKTDaiJ78QBUXBL0RBUfALUVAU/EIUFAW/EAWlvWq/RSUzpXjmgaXGunPlmKqmqfRcouwzBT5p2kGU7tpkwqSElQbbtZN2tfFoGFF+7S2+2eGh0OblhNswcaPteezbcZv01aDnsURcmAHAzkdX36GDXCkvzcfrcPoneSrwQ//qP4a2x//B+0Pbh/7RP6evr5wn5iNEqV8eWOPHm7zubLPkHqvN8DgJswiJ+57up+GeQohrCgW/EAVFwS9EQVHwC1FQ2pzea1EESznfMhfVxJpmlkKa5XH/zeUHQFx2N8W16QCAJlOBbY6ndM792I2hrdbJ38e7zhH33rmY7goAUzdHoWnglViqqpZwqGXHlhK6SrujmLkwyK/vUl8c19wOft988sFPx32RElrDr3JXYSYUJ0uhXRdLflVPceE1D9S5OiFKd+zYtuJvO9t4SOvJL0RBUfALUVAU/EIUFAW/EAWlEQ+/h81s1MwOXtL2W2Z20syer//7SGuHKYRYbxqRBh8B8J8BfH5V+4Pu/vt5dmYAzFaq+J5Qv5kJabImXsIwomESSmqZpGpmFy6EttrZc03viynl3stnEawaz0PnBe7I62QWwOd5eu/Ak4dJ38bTodkxlPu5cYiNR0V709jqwlD1vhbHO/BNfutm52PJyIwcw9JPvJe/njgTVyb5rMv8ljg70Z2YGcimYjqzdcYaggDgi3E2hh0DANTOrTxerzY+m7TWcl1CiA1OM9/5P2VmL9S/FnAPKSHEVctag/+PANwE4ACA0wD+INXx0lp9i6rVJ8RVw5qC393PuHvNl8uDfBZALJPyo74/rNXXqVp9Qlw1rCm918xGLlbpBfDzALid6yrcHRkRMxreL0n5Xd5w4yWKqDNqYr02S7NMOerSYeVI72WurzYby2cBQOfZKBQtDEc3XAAozxERbjKxNpyMN8/1oud2jh9DxtanJ66jdcZzw9KsAaC0N6bcdpC+r9zHH0Q3/I8onJ5/92bad/D7JA034UZtrPRaSqgux/PQsSu6KANANj6xcj/Vxt17Lxv89XJddwPYZmYnAPwmgLvN7AAAB3AcwC83vEchxFXBWst1fa4FYxFCtBFl+AlRUBT8QhQUBb8QBaX9tfpWqd3MeRcAMqKOpkwVGCnjjzz11UBmF1idPFhr3kM9VcttS0w77n6t8STM7K1R3k5Slyk53I5TzsZ0hiUxE1IaGozb7ecK/LF/Ek1C+o/GdOg9f8nvpc6T0ahk23Ok3iEAI6nLNB0aAHZuD03ZFp6+XZqM18HPxbRlALDVs0/zjd+LevILUVAU/EIUFAW/EAVFwS9EQWm/4LdKLGLCXu5NkpTb1PpnRiq915wIfkzUSgqGOXwGWGprYrt2gYihCZfc2tRMbEw5CLPxsr45Xp/r+ibO1/TfuSG09Z7g4mTfG1Hc2/xG4/fC/PWxDFhPwkWZujMvcV8Flqptr5/g271hT+xLUpwBYPGmlWm//j3uEcDQk1+IgqLgF6KgKPiFKCgKfiEKioJfiIJyxdN7U3QMxzTNlPnBakODt91Po4o2GvcI8YTCm0ftZ7MInlDKa2Mk3TSxrxJxiM3mE2m0DfZNpeyyNGfr5eYnLA22vJsbVpQX44UYez83L9k0GvuWFuMxTN3IU2sHv/5qaKsx4xEAJZJCbltjbUOAH68Rd2gAwPmYYuyJdObOoytrA9pC4+YrevILUVAU/EIUFAW/EAWlkXJde8zsKTM7ZGYvmdmn6+2DZvaEmR2u/5R3vxAbiEYEvyqAX3P375rZZgDPmdkTAP4ZgCfd/TNm9gCABwD8+tttyHq6UXrH/hVtS0NcEBq/Poop/ce5AFZ5gTjEptJ7GxQck+QQDJvdfzYdSzzlJU96LRMdS5uj0JTNptb9k5Jjib7lm/eFtsWdvLTXhW3xNh06mHAFJuXJxg5EYW3gaEKkJSXhSglhrjYehblyJRFSZLvMWwIAbM+u2EhEQACoTaxsX+9yXafd/bv136cBHAKwG8B9AB6td3sUwEcb3qsQ4oqT6zu/me0D8D4AzwDYedG7v/5zx3oPTgjROhoOfjPrA/AVAL/q7nwJGX/dj8p1Vfl8qRCi/TQU/GZWwXLg/7G7/1m9+YyZjdT/fwQANYZbUa6rI5HUIIRoO41U7DEsF+k45O5/eMl/PQ7gfgCfqf987HLbmt9Rwiv/cqWAVDnHh3DjV+I69PO38ayu7UditlbHFt43GzsX2xLGldQngPXNUa++VbDsPCCdzcdwJkoRwS69ryguprwO2Lr5pc38XuiejMdQOdm4YekQtoW2rJwoa8VKwiX8Hjp2xO36QkJILLPsR/4wNFJerHquweP1eA1TNKL2/ziAfwrgRTN7vt72G1gO+i+Z2ScAvAHgFxreqxDiitNIua6/AZCq/vfB9R2OEKJdKMNPiIKi4BeioCj4hSgobV3P39uzgA/cdnRF2+HPv5P2Zeu1p/fx7Q7+ydmGx5Bce09ouHxVihzr+Y2khabSPxnr4ZJb6o4p1cwFOc++jKjcAND97JHQNvWx22jfnX95MrR5YsaBueQu9cZz2/fcG/T1TFVPlX6zxbh2vkRKeAEAmFfCGV4GzNl9l5g5CrMpWUqei+jJL0RBUfALUVAU/EIUFAW/EAWlrYLfhalufO+vVgp8Q5M8HXF6T3xfuvnzXCBBT3doSpkutgIm1gH5BDtfiIJOKZH+mc1Fwa08GNNlAaB2NqYzp5xJGy1xljLwZOnQ1XfFUlsAML4/9h04wvdPxb2MH8PijdH4ddOJmCruc9wPgImhqevoTDi1hAdDKd7PKeGUnl/jQl4QsHOk9+rJL0RBUfALUVAU/EIUFAW/EAVFwS9EQWmr2m8ZUJldqVrODnMVc/hpkuL4Fk/jrU1FV7GUIt1oCa71oNQdZyG8xgdQ6o2GJLWpqFID3EzDL3D1mplppFJuG03bLW3hJamwc3toqpziJhR9/XFciwP8mnUdj2m0taFE+aqT46Et64vXISNpwCnYdQQAJ8p6NpuYZWqydFsr0JNfiIKi4BeioCj4hSgozZTr+i0zO2lmz9f/faT1wxVCrBfNlOsCgAfd/fcb3psDtkrL2PVUFGgAAIdfD021VEom2xVxogWQq9wWEw1Lt5B0VeK2CgCYjuIPWzMPACBrw5EQ/IxsY3XZph/ur1EHYqBhUcoSAlh26kwc1637aN/5rXFfpSq/ZjPvGQ5tva8kPByI8GlTJOU2IXoyvwev8fPiCVdfShudnBulEQPP0wAuVuaZNrOL5bqEEBuYZsp1AcCnzOwFM3tYVXqF2Fg0U67rjwDcBOAAlj8Z/EHidT8s11WbU7kuIa4W1lyuy93PuHvN3TMAnwVwJ3vtpeW6yj0q1yXE1UIjaj8t13WxTl+dnwdwcP2HJ4RoFc2U6/q4mR0A4ACOA/jly22oMpthZFXarr9yjPYtMcfUHEYFSYMNoswn3VnLxNihM263uo3XBayQNNzF7fzTT8czh0JbeWsijZbNQhBVHwCMGWHM53AVJrXrqqNcaS9vHwptHWO8oHM/SXNe2BZTnAGg55nDsZHUyQMA30zOL6vPmMOBOI/jc4pWpFk3SzPlur62/sMRQrQLZfgJUVAU/EIUFAW/EAWlvev555fQceTUirZaIjU2I+m5TDQBEimZqZRbkma5PJMZsU1RgBq9M7rk1rq4J8HWV0kK6xJfz297r4uNp0dpX5Bzk0zZzVFxjJ3f0pYoZjpxDwYAJ47J1ZtGSE+gPEuuGREXAcC2xnPO1tIDAM7GdHFnqdOJVGYq8qbSeHOk7NJ7tNp4aa0UQaxeULkuIcRlUPALUVAU/EIUFAW/EAVFwS9EQWmr2u/VKmqrUi1TabhM0c7FOpgn+HB0o13cEtXU+e18rLO74izCjucSan81utFWOrgiXTt0hLY3SmrWpETq/WW74jkoTfMpBO8lbsUV/nyZeFdMXR44OMG3S8w4fJaPgaXGljYTp9/E/eE57ht2HnOlAhs/N8aue2ImJNQRVK0+IcTlUPALUVAU/EIUFAW/EAWlrYIfgCC0+EIOgSVRgouRKrHEBKHU+ukySQutTEdBpWeMj2H0J2KK8bn38FTiylRMJfYOvr79+qmYRlt7i6cCs9RUJuwBAAajCFeejPta2DtIX17rikJVKvW55yzxVTgd190DoO7IyTXvLG13iaT3psjh7pxL3Muz3TY5/erJL0RBUfALUVAU/EIUlEYMPLvN7Ntm9v16ua7frrffYGbPmNlhM/uimfHMESHEVUkjCtoCgHvcfaZu4f03ZvZ1AP8ay+W6vmBm/xXAJ7Ds5Z+LXGv0c9QtzyUIpQSWLGbjbXk9ikfzQ/w09h6N4l7XOM/AqszE9qH/d4r0BKpvnoiNifXp5f4toW3xVuIdAL6evnM0Cn4dEwuhDQA6j5PxDnOjTZyN2XzGzDcBVM+d59tokFz17psV2yyxnr7J7abE7lzHtorLPvl9mYtF4yr1fw7gHgBfrrc/CuCjax6FEKLtNFq0o1y37R4F8ASAowAm3P3i284JqH6fEBuKhoK/XpnnAIDrsFyZ51bWjb320nJdS+AfF4UQ7SeX2u/uEwC+AeAuAANmdvGLyHUA6BfUS8t1VZAoTy2EaDuNqP3bzWyg/nsPgJ8CcAjAUwA+Vu92P4DHWjVIIcT604jaPwLgUTMrY/nN4kvu/j/N7GUAXzCzfwfge1iu55ebplMkwUtKecoPIIfqWjv2RmhjCbcdF3i6rGXxk07PGE817RyLqnqWSNllpbmyef6Vynpj345ZPobx/dGpt7cjntuZEZ6ivI2s87cLfNaFuQ3XxhJ50nkg19e5hQKFrtHP496bYz19HppR9VM0Uq7rBQDvI+3HkKjMK4S4+lGGnxAFRcEvREFR8AtRUNq7nt9immIoN1Qnm40CWHL9M0hZLFaXHkA2n8M/gJRuAmmzKleUusajSNP1HDffrE3N0HYKOQ8dI8O069LeaMA5v42nVA/8II6h43Qsf1WZ2MrHNT4ZmqqkDcC6GKy2glat0c/Tl8VEMOpcB/TkF6KgKPiFKCgKfiEKioJfiIKi4BeioLRX7feYdutM1Qc3L0imOBLVNKXq59kuU35rP4hqvSXSjru7o2qbESfa5f+I403NhJS3E2U/MbtRno3H0PfiMT6Ed+4NbbVhouw/+zJ9vfcQx+RWqfqJc96S/a3HvnL0bYWyz9CTX4iCouAXoqAo+IUoKAp+IQpKWwU/Mwtpt15LrNEvx/el1Bp9tp4/OYY8QmKj26zw08jWrKfEI1ZeLOVAXD31Vmjr2MXTexlsjT8AlI6djH3JuHxLXPcPALUJksrbKmGunenBV2kqcrPoyS9EQVHwC1FQFPxCFBQFvxAFpZlafY+Y2Wtm9nz934HWD1cIsV40U6sPAP6Nu3/5bV67ErNghuEJRdu5wSzvm8OdlSr7eVyB2esT7q7MCZbNYgBcgZ/9mdtp3/5vvh7aFm/YQft2Ho7lFBb28+JKlYl4LexMrJPnCadgSkIpzzO7QV9PHIyBxAwLIU+NyJbRzhRlQiPuvQ6A1eoTQmxg1lSrz92fqf/X75rZC2b2oJnRVSiXluta9Mbf2YUQrWVNtfrM7N0A/i2A/QD+FoBBAL+eeO0Py3V1Gln1JYS4Iqy1Vt+97n66Xr57AcB/hwp4CLGhuOx3fjPbDmDJ3ScuqdX3e2Y24u6nzcwAfBTAwctty7OMu/IyUmJIk9BU3IRgx8S90ubNseMSVydpGu22Qdp38vZtoW3zMe7o672xaBgT9gAg2xHX43eO8u1mvcR/YHwiti3mUGMT5BH32L2QFPZIX+bCnCzB1SJakVbeLM3U6vur+huDAXgewL9o4TiFEOtMM7X67mnJiIQQbUEZfkIUFAW/EAVFwS9EQWmve28emJstUUyBtMkH7ZvHGZWpzNPTcVyJVFHq6nsu1r4DgL7jUcGHcZOSpZGB0FZ5NeGCfGqMDYz29ZdJX3bO89Sjy0Nqu0265Dp7fYtmkzYSevILUVAU/EIUFAW/EAVFwS9EQWmv4GcWxLHU+mm63juVVppLHGxsXwBPQS0PkfTcRAmujLjZlndzl92FoTiGnlM8DbcyEcU9n5ziYyDHkFoLz0QwloKaS3jNIdYlt5vDg4HCxpAYV6vScK90Ki9DT34hCoqCX4iCouAXoqAo+IUoKAp+IQpKe9V+94bdUZlKnUfBT6qrLGU3MYvAjDt8NppIpPbFVPXaVmIGAqC6Kb4P22LqGMh7NjGsAIDyli1xDDMJQ5UGlflWKdd5rhkz6ACad9+1rmhosh7HezWaeejJL0RBUfALUVAU/EIUFAW/EAXFlgvytGlnZmMALtaa2gbgbNt23j50XBuPa+nY9rr79kY6tjX4V+zY7Fl3v+OK7LyF6Lg2Htfysb0d+tgvREFR8AtRUK5k8D90BffdSnRcG49r+diSXLHv/EKIK4s+9gtRUNoe/GZ2r5n9wMyOmNkD7d7/emJmD5vZqJkdvKRt0MyeMLPD9Z+xUuZVjpntMbOnzOyQmb1kZp+ut2/oYzOzbjP7tpl9v35cv11vv8HMnqkf1xfNjHuxX2O0NfjrxT7/C4CfBnAbgI+b2W3tHMM68wiAe1e1PQDgSXe/BcCT9b83GlUAv+butwK4C8Cv1K/TRj+2BQD3uPvtAA4AuNfM7gLwewAerB/XOIBPXMExto12P/nvBHDE3Y+5+yKALwC4r81jWDfc/a8BnF/VfB+AR+u/P4rl8uUbCnc/7e7frf8+DeAQgN3Y4Mfmy1w0RqzU/zmAewB8ud6+4Y5rrbQ7+HcDePOSv0/U264ldrr7aWA5iADsuMLjaQoz24flKs3P4Bo4NjMrm9nzAEYBPAHgKIAJd7+4vvZavCcp7Q5+Vn9K0w1XKWbWB+ArAH7V3bk98AbD3WvufgDAdVj+JHor69beUV0Z2h38JwDsueTv6wCcavMYWs0ZMxsBgPrP0Ss8njVhZhUsB/4fu/uf1ZuviWMDAHefAPANLGsaA2Z20W3jWrwnKe0O/u/dOJWiAAAA4klEQVQAuKWurnYC+EUAj7d5DK3mcQD313+/H8BjV3Asa8LMDMDnABxy9z+85L829LGZ2XYzG6j/3gPgp7CsZzwF4GP1bhvuuNZK25N8zOwjAP4DgDKAh939d9s6gHXEzP4UwN1YXhV2BsBvAvhzAF8CcD2ANwD8gruvFgWvaszs7wL4vwBeBHDRJO03sPy9f8Mem5m9F8uCXhnLD74vufvvmNmNWBafBwF8D8AvuXuOcs4bE2X4CVFQlOEnREFR8AtRUBT8QhQUBb8QBUXBL0RBUfALUVAU/EIUFAW/EAXl/wNBmHlXcnYqggAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure()\n",
"plt.imshow(imcac_cropped_loaded.get_img_by_metal('Ir191'))"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "spherenv",
"language": "python",
"name": "spherenv"
},
"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.0"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment