Commit d708d573 authored by Ines Filipa Fernandes Ramos's avatar Ines Filipa Fernandes Ramos
Browse files

change of location of notebooks

parent cd8c412e
This diff is collapsed.
%% Cell type:markdown id: tags:
# Simple RGCs simulation
%% Cell type:code id: tags:
``` python
import numpy as np, array
import matplotlib.pyplot as plt
import matplotlib.gridspec as gridspec
from scipy import signal
import matplotlib.image as mpimg
import matplotlib.cm as cm
import math
from mpl_toolkits.mplot3d import axes3d
import torch
from collections import OrderedDict
import neuralpredictors as neur
from neuralpredictors.data.datasets import StaticImageSet, FileTreeDataset
import random
from numpy import save
```
%% Cell type:markdown id: tags:
##### Simple RGCs simulation:
%% Cell type:code id: tags:
``` python
_default_2Dgaussian_p = (1,1,1,0,0,0,0)
def gaussian_2D(xz, sigma_x, sigma_z, amp, theta, x0, z0, y0):
"""Two dimensional Gaussian function
params:
- xz: meshgrid of x and z coordinates at which to evaluate the points
- sigma_x: width of the gaussian
- sigma_z: height of the gaussian
- amp: amplitude of the gaussian
- theta: angle of the gaussian (in radian)
- x0: shift in x of the gaussian
- z0: shift in z of the gaussian
- y0: shift in y of the gaussian
"""
(x,z) = xz
x0, z0 = float(x0), float(z0)
a = (np.cos(theta)**2)/(2*sigma_x**2) + (np.sin(theta)**2)/(2*sigma_z**2)
b = -(np.sin(2*theta)) /(4*sigma_x**2) + (np.sin(2*theta)) /(4*sigma_z**2)
c = (np.sin(theta)**2)/(2*sigma_x**2) + (np.cos(theta)**2)/(2*sigma_z**2)
g = amp * np.exp( -(a*((x-x0)**2) + 2*b*(x-x0)*(z-z0) + c*((z-z0)**2))) + y0
return g.ravel()
def mexicanHat(xz, sigma_x_1, sigma_z_1, amp_1, theta_1, x0_1, z0_1,
sigma_x_2, sigma_z_2, amp_2, theta_2, x0_2, z0_2, y0):
"""Sum of two 2D Gaussian function. For the params, see `gaussian_2D`.
However, both share the y0 parameter."""
return (gaussian_2D(xz, sigma_x_1, sigma_z_1, amp_1, theta_1, x0_1, z0_1, 0)
+ gaussian_2D(xz, sigma_x_2, sigma_z_2, amp_2, theta_2, x0_2, z0_2, 0) + y0)
def ELU(r):
if r>0:
return r+1
else:
return np.exp(r) + 1
def RF_circular_CS(vis_field_width, vis_field_height, x_rf_center, z_rf_center, polarity, plot=False):
x,y = np.meshgrid(np.linspace(0,vis_field_width,vis_field_width),np.linspace(0,vis_field_height,vis_field_height))
if polarity==1:
sigma_x_1, sigma_z_1, amp_1, theta_1, x0_1, z0_1 = 2, 2, 1, 0, x_rf_center, z_rf_center
sigma_x_2, sigma_z_2, amp_2, theta_2, x0_2, z0_2, y0 = 3, 3, -0.5, 0, x_rf_center, z_rf_center, 0
else:
sigma_x_1, sigma_z_1, amp_1, theta_1, x0_1, z0_1 = 2, 2, -1, 0, x_rf_center, z_rf_center
sigma_x_2, sigma_z_2, amp_2, theta_2, x0_2, z0_2, y0 = 3, 3, 0.5, 0, x_rf_center, z_rf_center, 0
z = mexicanHat((x,y), sigma_x_1, sigma_z_1, amp_1, theta_1, x0_1, z0_1,
sigma_x_2, sigma_z_2, amp_2, theta_2, x0_2, z0_2, y0).reshape(vis_field_height,vis_field_width)
if plot==True:
fig = plt.figure(figsize=(5,4))
ax = fig.add_subplot(111, projection='3d')
ax.plot_wireframe(x, y, z, rstride=3, cstride=3, label=f"x_rf_center={x_rf_center} z_rf_center={z_rf_center} \n amp_center={amp_1} amp_surround={amp_2}")
_ = ax.legend()
return z
def RGC_response(rf, image, plot=False, seed=None):
Img_barHat = image * rf
if plot==True:
fig, ax = plt.subplots(3, figsize=(8,8))
ax[0].imshow(image)
ax[0].set_title("Image")
ax[1].imshow(rf)
ax[1].set_title("RGC RF")
ax[2].imshow(Img_barHat)
ax[2].set_title("RGC Response")
plt.tight_layout()
if seed is not None:
np.random.seed(seed)
g = ELU(sum(Img_barHat.ravel()))
spikes = np.random.poisson(lam=g, size=None)
return spikes
```
%% Cell type:code id: tags:
``` python
#Generate the receptive field of one RGC
rf = RF_circular_CS(vis_field_width = 64, vis_field_height = 36, x_rf_center = 30, z_rf_center = 30, polarity =1, plot=True)
```
%%%% Output: display_data
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAP4AAADvCAYAAADSHGlvAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAACt2UlEQVR4nOyddXhcZfbHP3cs7t54kzQu9VJaCpTisshiRRcWyuKyLPyQRXZZbHFdtIsVh0Jx2lIK9Tbaxt2TycTH5/39Mbm3SRpPKgv9Pk8fmMmdd+69c8/7nvec7/keSQjBYRzGYfy+oDrYJ3AYh3EYBx6HDf8wDuN3iMOGfxiH8TvEYcM/jMP4HeKw4R/GYfwOoZnEZw+nAw7jMPYvpP018OEV/zAO43eIw4Z/GIfxO8Rhwz+Mw/gd4rDhH8Zh/A4xYnBvx44dwRqN5lUgjUGThE6n25/ndRiH8buBq6srERERaLXaA/adIxq+RqN5NTQ0NDkoKMigUqkGR/Gj9+N5HcZh/C4ghECv11NbW0tsbOwB+97RXP20oKCgziGM/jAO4zCmAJIkERAQgMlkOqDfO5rhqw4b/WEcxv6FJO23dP2wOBzcO4zD+B3isOEfxmH8DvGbNvxnnnmG5ORkli9fvl+/Z/369fz666/7ZewTTzyRzMxMUlNTWbFiBXa7HYC2tjaWLVtGQkICy5Ytw2Aw7JfvB/jwww9JTk7mmGOO2W/fAZCdnc1XX321X8a+4ooryMzMJCMjg3POOYfu7m4AzGYz5513HvHx8cyfP5/Kysr98v2HGn6Thm+z2QB44YUX+P7773nnnXf26/dNxPDlcxwNH3zwATk5OeTn59PS0sKHH34IwMMPP8zSpUspKSlh6dKlPPzww+M+79EghMDhcPDaa6/xyiuvsG7duin/jv6YiOGP9T4++eST5OTkkJubS1RUFM899xwAr732Gn5+fpSWlnLzzTfzt7/9bdzn/T8JIcSw/7KzsyuFENuH+bffsXXrVpGeni6MRqPo7u4WKSkpIi8vb8hj161bJxYtWiROO+00kZCQIK6++mqh1WpFWlqaeOKJJ4b8TFdXl7jssstEWlqaSE9PFx999JEQQohvv/1WLFiwQMycOVOcc845oqurSwghRHR0tLj33nvFzJkzRVpamtizZ4+oqKgQISEhYtq0aSIzM1Ns2LBBNDc3i7POOkvMmTNHzJkzR2zcuFEIIcTf//53cdFFF4mFCxeK888/f1z3wmKxiFNPPVWsWrVKCCHEjBkzRH19vRBCiPr6ejFjxowRP3/PPfeIzMxMkZmZKaZNmyYuu+yyIY+rqKgQM2bMEBdffLFISUkR9913n/Dw8BAzZswQt91225Cfsdls4tZbbxWpqakiPT1dPPPMM0IIIbZv3y6OOuooMWvWLHH88ccr57tkyRJx++23i7lz54qEhASxYcMGYTabRWRkpAgMDBSZmZli1apVoru7W1x++eVi7ty5IisrS3z22WdCCCHeeOMNcdppp4ljjjlGHHXUUeO6jw6HQ6xYsUI8/PDDQgghjj/+ePHrr78KIYSwWq0iICBAOByOcY05Fdi9e/dQb49on5P5d0gbvhBC3HXXXeLWW28Vf/nLX8RDDz007HHr1q0T7u7uory8XHkvOjpatLS0DPuZ22+/Xdx4443K67a2NtHS0iIWL14suru7hRBCPPzww+L+++9XxpMf6ueff15cccUVQginQT/22GPKOBdccIH4+eefhRBCVFVViaSkJOW4WbNmid7eXiGEEIWFhYoxDv5nMBiU8Y4//njh6+srLrjgAmGz2YQQQvj4+Ch/dzgcA16PBIPBINLS0sT27UP/hBUVFUKSJLFp0yblvSVLloht27YNO+YLL7wgzj77bGG1WoUQQuj1emGxWMQRRxwhmpubhRBCrFq1Slx++eXKeLfccosQQog1a9aIpUuXCiGcBn3ttdcq4955553irbfeUs47ISFBdHd3izfeeEOEh4cLvV4vhBCis7Nz2PtYUFCgjHfZZZeJ4OBgcfTRR4uenh4hhBCpqamipqZGOWb69OkjPjP7Cwfa8CdTlntAcO+99zJ37lxcXV155plnRjx23rx54yJB/PDDD6xatUp57efnx5dffsnu3bs58sgjAbBYLBxxxBHKMWeddRYAs2fP5pNPPhl23N27dyuvOzs7lT3l6aefjpubGwCJiYlkZ2ePep7ffvstJpOJ5cuXs3btWpYtWzbg75IkjSklJITgoosu4pZbbmH27NnDHhcdHc2CBQtGHU/GDz/8wIoVK9BonI+Tv78/+fn55OfnK+dqt9sJCwtTPtP/Pg63r/7uu+9YvXo1jz/+OAAmk4nq6moAli1bhr+/PwBeXl5juo9vvPEGdrud66+/nvfff5/LL798zNf4W8Mhb/h6vZ7u7m6sVismkwkPD49hjx3pb2OFEIJly5bx3nvvDfl3FxcXANRq9bD7S4fDwebNm3F1dR3xHIuKijjvvPOGHGP9+vX4+voqr11dXTnjjDP4/PPPWbZsGSEhITQ0NBAWFkZDQwPBwcGjXtt9991HRETEqA/8VN3H1NRUNm3aNOTfx3IfhRB8/PHHJCYmDnh/y5YtA86xq6uLxYsXDznGu+++S0pKivJarVZz/vnn8+ijj3L55ZcTHh5OTU0NERER2Gw2Ojo6CAgIGNe1/i/ikA/uXX311Tz44IMsX758ygMvy5Yt4/nnn1deGwwGFixYwC+//EJpaSkAPT09FBcXjziOl5cXXV1dyuvjjz+eZ599Vnk93Gokr/jZ2dns2rWLXbt2sWPHDrZv346XlxednZ3U19cjhMBms7FmzRqSkpIAp+ewcuVKAFauXMkZZ5wx4jl+8cUX/PDDD6N6TRPBsmXLePnllxUDbmtrIzExkZaWFsXwrVYrBQUFI44z+D6ecMIJPPvss849KbBr165hPyffx8H/UlJSEEIov6cQgtWrVw95Hz/66COOPfbYg0KoOdA4pA3/v//9L1qtlgsvvJA77riDbdu2sXbt2ikb/+6778ZgMJCWlkZmZibr1q0jKCiIN998kwsuuICMjAyOOOIICgsLRxzntNNO49NPPyUrK4uff/6ZZ555hu3bt5ORkUFKSgovvfTSiJ/vv/fq/15XVxenn346GRkZZGVlERQUxNVXX40QgjvuuIPvv/+ehIQEfvjhB+64444Rv+OJJ56grq6OefPmkZWVxb333jv2GzUKrrzySqKiosjIyCAzM5N3330XnU7HRx99xN/+9jcyMzPJysoaNfNxzDHHsHv3brKysnj//fe55557sFqtZGRkkJqayj333DOh8xNCcOmll5Kenk56ejoNDQ3K9V9xxRXo9Xri4+N54okn9kt25FCE1P9hG4ycnJzKzMzM1mH+PPwm8TDGjKEMfrjjbDYbQgi0Wq2yr5dXp9/DKvVbxp49e0hOTh789n77UQ/5Pf5vFYMNfjT0N2xJkvb5/OGJ4DDGg3EZ/hnPbVSiLK5a9aS//P2rjxj9oEHIy8vj4osvHvCei4sLW7ZsGfYzb7zxBk8//fSA94488sgB+/sDifEa/VAYbNjjnQj0ej1Lly7dZ9wff/xx2ODWt99+u0+cJTY2lk8//XRC13AYBw/jcvUHGb7nZL98Iob/v47RjH6kv1mtVoQQYxJBGTzOYY/g0MYh7ep/ft2ion4vD+/xx4GpWOXlccaCyXoEh/HbxiEd1f8t4KGHHhq30Q93bH/jLCws5Mgjj8Td3Z1///vfo47Z39DlGIHD4cButyv/HA7HlE1Q/yvw9Jy447py5UoSEhJISEhQUoKDcd999xEeHk5WVhZZWVn7rQhpvDgc3NuPEELw0EMPjZpqGwy73a6w4IaDv78/Tz31FJ9//vmEzu1Q9QhsNtuAax/8+lBBW1sb999/P9u3b0eSJGbPns3pp5+On5/fPsfefPPN3HbbbQfhLIfH/9SKf8011zBnzhxSU1P5+9//rrwfExPDnXfeSVZWFnPmzGHnzp2ccMIJxMXFKTn09evXc9RRR3HKKaeQmJjIihUrcDgcw37XN998w6xZs8jMzFSCYD09PfzpT39i3rx5zJw5UzG6N998k7POOosTTzyRhIQEbr/9diXXbjQamTVrFhdddBEAb7/9NgsWLGDWrFkDymy9vb257bbbmDVr1rBst/4IDg5m7ty54xJofPDBB1mwYAGZmZmsWLFCMfRjjz2WW2+9lQULFpCens727ds555xzSExM5K677sJut1NWVkZSUhLLly8nOTmZc845h97e3mG/64477iAlJYWMjAzlob/sssv46KOPlGPk1Xb9+vUsXryY008/nZSUlH1em0wmLr/8ctLT05k5c6ZSJfjmm29y3XXXKeOdeuqprF+/Xhn7rrvuIjMzkwULFtDU1ARARUUFRxxxBOnp6dx9991jvneD8e233yq0YT8/P5YtW8Y333wz4fEOOEYi8h8KRTr9IRdl2Gw2sWTJEpGTkyOEcBbPvPDCC0IIIW666SaRnp4uOjs7RXNzswgODhZCOIt4XFxcRFlZmbDZbOK4444TH3744ZDf09zcLCIiIpSCH/l7RyoaiY2NFe3t7cJoNIqoqChRUVEh7Ha78PDwEHa7XdjtdpGfny9OOeUUYTKZhN1uFytWrBBvvPGGsNvtAhDvvfeesNlswmq1ihtuuGGfgpOMjAxx//33C5vNpvy75557xCOPPDLgveH+NTc3K/+/fPly8emnnwqbzSaOOuoocdtttwmbzSaeeOIJERYWJmpqakRPT48IDw8XTU1NoqSkRABi7dq1wmq1issuu0w8+uijwuFw7FPN1traKmbMmKG8LxccXXrppQPuuYeHh/Lb9C+wGvz68ccfVwp89uzZIyIjI4XRaNynqOeUU04R69atE0IIAYjVq1cLIYT461//Kh588EEhhBCnnXaaWLlypRBCiOeee045h7EW+sh47LHHlDGFEOKBBx4YUKgl4+9//7uIjo4W6enp4vLLLxdtbW1DPnOHi3RGwAcffMB//vMfbDYbDQ0N7N69m4yMDMBJvQRIT0+nu7sbLy8vvLy8cHFxob29HXAW8UyfPh2ACy64gI0bN3LOOefs8z2bN2/mqKOOUgp+5GKQkYpGli5dire3N0IIkpKSKCsrIyoqasC4a9euZefOncyfPx8Ao9GocOzVajVnn322cuwTTzyxz3nZbDbFQ5gI1q9fz2OPPYbRaKStrY2UlBROO+00AOW/aWlppKSkKAU1sbGx1NTU4OvrS0REBAsXLkSSJC688EKee+45br755gHbAUmS8Pb2xtXVlSuuuIJTTz2VU089ddRzG1xg1f/1xo0buf766wFISkoiOjp6VBq1TqdTvnf27Nl8//33APzyyy98/PHHAFx88cVKenKshT7jxTXXXMM999yDJEncc8893Hrrrbz++utT/j3jxf+M4VdUVPD444+zbds2/Pz8uOyyywYok8pFHyqVSvl/+bXMIR+8Tx3vvlWMUDSi0+kU13m4whMhBJdccgkPPfTQPn9zdXVFrVYrY9xyyy389NNP+3z+7LPP5q677hrXeYNzkrruuuvYsmULkZGR3H///VNy/+T35JVE/symTZtYu3YtH3/8Mc899xxr165Fo9Eo2yuHw4HFYlHGGlwYNJZCof7jydcoQ2Y3wr6/x1C/+2iFPl1dXVx99dUAPPDAA4SHhyvbCoDa2lqOPvrofT4bEhKi/P+f//znMU2CBwL/M3v8zs5OPDw88PHxoampia+//nrcY2zdupWKigocDgfvv/8+ixYtGvK4BQsWsGHDBioqKgBnIAeGLxqRX/eH/J5Wq8VqtQLOvfTHH39Mc3OzMm5VVdWAz1ksFiwWC48++ijbt29nx44dyr+tW7dOOEgkG0VgYCDd3d3DlhSPhJqaGjZv3gzAqlWrlNJlGBgM7O7upr29nRNPPJHHH3+cnJwc7HY70dHRbN++HYDVq1cr92U0LF68WFFRKi4uprq6msTERGJiYsjOzsbhcFBTU8PWrVtHHevII49USrH7KzONVugzf/585fXpp5/OCSecwHfffYfBYMBgMPDdd99xwgkn7PN9DQ0Nyv9/+umnpKWljema9zf+Z1b8zMxMZs6cSVJSEpGRkQMeurFi7ty5XHfddZSWlnLMMcdw5plnDnlcUFAQ//nPfzjrrLNwOBwEBwfz/fffc88993DTTTeRkZGBw+EgNjaW1atXj5j++vOf/0xWVhYzZ87k7bff5oEHHuDEE0/E4XCg1Wp59tlniY529ibp7e1FpVIpq6zZbEaSJNRq9T6R7cbGRubPn09nZycqlYpnnnmGvLw8vL29hzwPX19fRXcuJCSEOXPmjPv+JSQk8PLLL7NixQqSk5NZsWLFkMd1dXVx5plnYjabEULw2GOPIYTgT3/6E2eddRaZmZmccMIJeHh4jCl1+Je//IVrrrmG9PR0NBoNb775Ji4uLhx55JHExsaSkpJCcnIys2bNGnWsp59+mgsvvJBHHnlk1IrGkeDv788999zD3LlzAaduhLwlvPLKK1mxYgVz5szh9ttvJzs7G0mSiImJ4eWXX57wd04lfjdFOuvXr+fxxx/nyy+/nJLx+ru2g2G327FarUPW4w8Fm82GyWQa4O7L7mj/XLvsrup0OtRqNWr15GnTY0VlZSWnnXYa2dnZk/7ewfdtcNrw90gmOqSZe4ex96EdbaUay0oGTpVXm82Gu7s7KpVqn8/JHoBWq8Vms2Gz2ZAkCYvFgsPhQKVSKR6BSvW/sXMbK4dg8P8fxtThd2P4Rx999JDBl/nz52M2mwe899Zbb5Genr7PsSOt8uOFEAKj0YhKpRqz4o1sAFqtVsnfy96A2WxWJgKNRoNarZ7SiSAmJoZt27bt8/7ZZ5+9j3TWQw89NOR+dzgcqmSi3zJGM3yHw+GQfstttEaq6uuPqTT6/q79ZFlp/V1+0Y+GazKZEEIof9doNPvFYOTU2FRipInAarWi0+l+UxPBVD1X48FoT11+S0tLSlBQUMdv2fhHwlhd+/4Y6UEc7NpPJeRAoFqtVtKLskdgNBqViUD2CP5XDKb/eVqtVrRa7W/GIxDC2S13rPGgqcKIwb0dO3YEazSaV4E0BqX+dDrd76JN9kRmYyGEsjL1f89isSj79Yl8n7yaT6aPusPhUP7B3hjCWCYhq9W63zyHscJisYypLPl/aRJwdXUlIiJiqN91v538iIY/Cn7THoBMMJEj7ON5gHp7eykqKmLmzJmAU/SisLCQxMREAgMDR/ysnBEYyhBl8U1ZKHKysFqttLe309bWRmdnJxqNBj8/P/z8/PDy8trnHPbs2UNERAReXl5T8v0TwbZt25QU2lCQtwUOh0P5zTQajfLvf8nT4XBU/8BB9GnbydHzibjjKpVKKXEtLS2lvb2d2bNnH3B3bjRotVqCgoIICgoCnNsQg8FAfX09XV1duLq6KhPBWHPu+xujnYM8Scu/m7zdsVqtAyYCrVarZEL+hyaCKcNhw+8HIQSdnZ37BI/GC0mSsNvtbN++HT8/P+bMmTPmsUY6Tq6j319wcXEhNDSU0NBQhBCYTCYMBgOVlZX09PRgt9sVroGbm9v/hMEMNxHIwqVyOvT3NhEcNvw+2Gw2rFYrO3bsUApRJgqDwUBHRwezZs2a0uYM+9vwB3+Xm5sbbm5uTJs2DSEEubm5OBwOSktLMZlMeHp64ufnh7+//wB+//4+r8l+vv8Y/ScC+e/9twa/1Yngd2/4g137yaz0slG0t7fj6en5m+rIIhtEWFgYHh4eOBwOuru7MRgM7N69G6vVio+Pj7I1mEwA8kBiqIlAXgTkv/8WJ4LfteE7HA6sVqsSCJrMD2oymcjNzSUgIIBZs2YpxShTiQO54o8GlUqFt7c33t7eREdH43A46OjowGAwUFtbi8PhUCYCX1/fKVHRORDXPtREYLVaB0wEcrpUJkn9L04Ev0vD7x/wASadT29paaG4uJjk5GT8/f2VwN5vDXKGo9di54q3c7nrxHjSpjkj/CqVSlntAaUPnRwjkCRJ+bu3t/eE+P79I/UHCrKhyxBC0NbWRlNTEwkJCUiSpMQH/pcmgt+d4cszuN1un/Qq73A4KCkpoauri7lz5yr55cmszL29vRgMBvz9/fdxlw+VFb9S30t+QxfrS/SK4Q+GRqMhICBA2e7IqcPm5mZKS0tHTR0OBTkYdzDRPzMgF1RZLBaF9i3zNOStwaFaa/C7MvzJ5OYHw2g0kpubS1BQELNnzx4w1kTHbWxspKysDD8/P2praxFCKMEzHx+fCZ/rVEG+b/oep6dU1NQ95s+ON3U41D3sX7V4MGG325UJaCiPQJ4I5GdMrq2QOQSHwjX8Lgx/vLn50R6w5uZmSkpKSElJGVJVdbxwOBwUFRVhMpmYM2eOsrLJq2RLSwulpaVKnX5XVxeenp4H7QFq7XEq5xQ29Ux4jNFShx4eHspEIKcO5SKkgw2HwzHsVqX/RCB7Z7K4CjjFWzo7OyelBTAV+M0bvjwDjzWAJ7vTQx3ncDgoLi6mp6dngGs/GZhMJnJycggODiYpKUk5X9h3lWxvb6ekpITq6mq6u7vx8PBQVF7d3NwmfS6jQX6Q9X2G39hppsNoxcdtchH8oVKHPT09GAyGAalDT0/PEZWRDxT6r/gjob/0FzjvX2FhId3dY/eU9hd+04Yvr/Ljce1l1t3gH1Z27YODg0lMTJyS1ba1tZWioqIBnsNIe3hXV1dcXFxITU0dYBzFxcWYzWa8vb2ViWB/pdP6u/oARU09zIvxnfLvkA09MjISIZwtw5uamuju7mbr1q0HNXVot9snFJyUJImenp5JNfGYKvwmDX8ytFvZ8PujqamJ0tLSKXPthRCUlZVhMBiYM2fOPuSXsUwqg43D4XDQ2dlJW1sbNTU1CCHw9fVV4gNTpdbjNHwLXq4aukw2ipq7p9zwh/pOORNgsVhITk7e76nDkeBwOCb8Hb29vYSGhk7xGY0fvznDn2xuvr/hy3tvo9E4Za69xWIhNzcXHx+fcVF5YeSovkqlwtfXF19fX8Dp7bS3t9Pa2kpZWZkSRff398fLy2tCHov83W09VuID3ak2GCe1zx8vZE/sQKQOR8JEV3xAiV8cbPxmDH9wbn6i0VPZ8Ht7e8nNzSU0NJSkpKQpce3b29spKChgxowZyr59f0Gj0RAYGKhUA8pR9NraWrq6unB3dx8QHxhPLUFrt4WYAGdMIa+uc79dw2AMF3vZH6nDkTCZIONhV38KMZW5eZVKRXNzM3V1daSmpior6GTPr6qqisbGRmbNmjXhQNxk8viDo+gyX0AOnnl5eSkTwXCejfzdLd0WGjpM9FodyvsHIsMwVoObitThSJjMit/d3X3Y8KcCDoeDuro6JEkiMDBw0oSczs5OzGYz8+bNm5Kgkc1mIz8/H51Ox7x58w6JdJQkSXh4eODh4UFERAQOh4Ouri7a2tqoq6vDbrcr8QFfX98BD7nV7qDT5CxoCfDQou+x0mW24e26/wNsE51gJpI6HAkjpfNGQ29v72HDnwz6u/ayfvtkjF527bVaLTNmzJgSo+/q6iIvL4/Y2FilJdVoOBhluSqVCh8fH3x8fIiNjVX2zG1tbZSXl6NWq/H393e60Ma9HWmi/d3Q91hp7rIcEMOfijz+WFOHI1UdjjWdNxQOu/qTwODcvFqt3kcpdzxobGykvLyc1NRU6uvrp8S4LBYLeXl5ZGRkTNkPfaAIO4P3zBaLhba2NsxmM7/uyFOOC3R3rnqt3Rbig/Z/wGp/bCmGSx0aDAb27NmDxWLZJ3U4mRX/sOFPEPIq3z83P1QKbqxjFRYWYrFYlJbTjY2NkyKJ2O129uzZg91uZ86cOVOSCeiPg8HV1+l0hIaGUl9fT1BQBGwtBMDF7myTnVNUQYybGT8/v/1al38gmHty6nCkqkOLxUJ7ezuBgYHjTuv19vYeVOkyGf8zhj9Sbl6tVo+7i2xPTw+5ubmEh4cTGRmprCQTnUTA+aPm5OQQHh5Od3f3flHRPdhFOm29TlffU6cmKSacL0rLkdy8MZlM7N69G5vNhq+v737JqR8Mrv5QqcPt27fT2dlJbW3tuFOHFovlgImWjIT/CcMfLTc/XsNvaGigoqKCtLS0fXrNTdTwZZJPWloaPj4+NDY2TthID5VilMEQQtBmdKZL/Ty0+PVRdU1CTUxMDDExMdjtdiU+UFlZOcBwvL29JzUZHgpcfbkqTy7JHW/q8GBP3DIOacMfa25+rMYqu+E2m4158+YNuRqN1/D78/f7ZwImOoFYrVZaW1vx9fXdZ5twKKz4+h4rKgn83XUEeDrPr713L4VXDgTKDSRlt7ixsZHi4mJcXFyUtOF4U2mH0oQon8d4UoeTSOOeCDwNqIFXhRAPD3Pc2cBHwFwhxIhKMIes4Q8liTUcxrLid3d3k5eXR0REBBEREcOOJ1eBjQWy6k5gYOA+/P2JGGlnZyd5eXn4+voq+8lDrSy3rceKWpLwc9cS1Gf4cnpvKOh0OoKDgwkODgacNQ+yNyAHuuSJYjQX+FBY8UdD/9QhOK/XYDBQXl7OFVdcgdFo5D//+Q9Lly4lPj5+1PEkSVIDzwPLgFpgmyRJq4UQuwcd5wXcCIypNdQhafjjpd2OtrrW19dTWVk5pGs/3rFkyFr5SUlJQ2rrDdUAcyTU1dVRVVVFZmYmWq0WSZKw2WwYDAalDFin02E2m+np6cHd3f2grH6tPRYE4OeuwdfN+fh0moc3/MFwc3MjPDyc8PBwhBB0d3fT1tbGnj179tHtG+yRHUor/ljRP3X4yy+/sGTJEiwWCy+99BKPP/74WIaYB5QKIcoBJElaBZwB7B503IPAI8BfxzLoIWX4gzXQxzq7D7fi2+12du/ejcPhGNa1HwyVSjWi9yCEoLy8HL1eP6JW/lg9B4fDoWQW5s2bpxSigHM/2d+N7O3tZdeuXVRUVChEEHm1nOrswXDQ91ixOwS+7lrcdc5AVo9pfIFVGZIk4eXlhZeXF9HR0djtdqXQqKqqakDgzMfH539ixR8JRqMRPz8/rr/++vF8LByo6fe6Fpjf/wBJkmYBkUKINZIk/W8Z/njr5vtDrVbvY2Syax8ZGUl4ePiYx5MFMIaCnJv39PRkzpw5Iz6EY3H15Vr8kJAQkpOTR/2Mu7s7Op2OtLQ0Jd/c1tZGfn4+DodjANtufxiIEAJ9j9W54rtpcdX2Gb5lYoY/GGq1ekAE3Wq1DvB47HY7Pj4+CuvwYKz+k4mxyBoKUwlJklTAE8Bl4/ncIWH4k5XEGrxK19bWUl1dTXp6+rhzpsO5+nKBTUJCgrJfHQmjGbHs3soCneNF/3xzTEzMgGq80tJSdDqd4g1MlZE4hKCjj7nn665FrZKQJDBap8bwB0Or1Q6IDxQXFwNQVVWlcN7lGMiB6lI0Gd2/CdJ164DIfq8j+t6T4YWzt+X6vt84FFgtSdLpIwX4DqrhT0W7KthrrDabjd27nVufsbr2w43V/xyrq6tpaGhg5syZuLu7T2ic/uNVVVXR1NQ0pW21BlfjmUymAUE0uQhnMtsCs31vw0R/d2f2QqOSMFkPjCqOWq3G19eXgIAAhWrb1tZGUVERZrP5gIhzTIauO8EVfxuQIElSLE6DPx+4UP6jEKIDUBoySpK0HrjtkI3qT8a1HwxJcras2rp1K1FRUUREREx4rP4GKxfYaLVa5s6dOy6a5lArfv+Cnblz5+7X/aqrqyvTpk1T+OjDbQvGI9LRj6aPbz/DN9sOjOH33+P3p9pGRUUNYNjJQiT94wNTVZM/GbruRAxfCGGTJOk64Fuc6bzXhRAFkiQ9AGwXQqyeyLkcFMO32WyUlZURFham9KmbKIQQ1NXVYTQaWbhw4aR50LLhd3d3k5ubS0xMDNOmTRv3OIMNX2YKRkVFER4ePqlznMi5DLctKCsrQ6vVjmlbYLTtvR6ZvKPTqOidoj3+aBgpqj8Uw85gMAwQIpGvcTJCpZNZ8SdamSeE+Ar4atB79w5z7NFjGfOAGn5/1769vZ2goKBJ0RdtNhsFBQWoVCo8PDympPhBpVLR1dVFbm7uhGIE/ceRPQeZ1Zeenj5qOvFAYKzbgsHc+/6GL6/4LhoVncaxp/Mmg/HsrwdnRMxmsyJL1tXVpZTi+vv7j4tYM1n1nUOhQAcOoOEPzs1rNJpx8+v7Qy55lVfkX3/9ddLnaLfbqaiowGg0cuSRR06KZy6n84qLi/dpuHGoYbhtgcy9lw3EaHUavloCLxfnw++qUSEAs82Bi2b/ptom00nHxcWFsLAwwsLCBpTilpSUYDKZFKHSoRiTg89hMiW5h4LsFhwgw+/fhFC+aWq1WulQOh4IIaitraW2tnafktfJEDzkenx/f39lYpoM5C47QUFBzJo163+GeDJ4W2C32zEYDOj1etp7zYAKD51KiarLufwukw0Xz/07sU1VHn84odL+FXjDxUAmq75zKDAw4QAZ/lABvIlU1MnBMY1Go5BdZIykhz8a5DxxamoqOp2OoqKicY/RHx0dHTQ0NBAZGUlCQsKkxjrYUKvVyrZA5LcCVrxd1Yq2Pzbnb6jv6iVwPxv+/mLu9RcqlYVI2tvb0ev1+wiVTnaPP5nA81TigBi+3GNs8HvjMfzOzk7y8/OHDbbJ443nR5FX5e7ubsUVN5lMk6rHr62tpaamhrCwsEOi7noq0du3xw/2dlO0/QPKckDfSV5hKZ21YljJrqnAgaLsjiRUajAY0Gg0Ck9iPEKl+4PAM1EctHTeWF19IQQ1NTXU1dWRmZk57I0bir03EsxmMzk5OUpb68nW4zscjgH04KqqqoNeSTfVkGtx5By+JEn4ebgCnYRFTWdWjA/t7e2KZNdURdJlHCzKbv/Cm/r6ekwmE8AAoVLZIxgpPnCo6O3BQTT8sQT3rFYrBQUFaLXafVz7wRiNY98fMmtuqAKbiRi+0WgkJyeHsLAwoqKiJqUKZDabqaiowNvbe8hClYMJecX399j7cPu4Os/PYLSiVqsHSHbJkXR5WzCeSryhcChw9R0OhxIM7S9UajAYyM/PH1Go9Hdn+EPN9P2LUYZCR0cHBQUFYxaqHMvWQQhBRUUFra2tw7LmxmuwcpXe4C47EynL7ejoID8/n/DwcDo7O6mqqkKlUuHv709AQMBBf2jkrJ2f+15WnNw3T+6n1x+DI+lyJV5/pZ7xbAsOheq8wV10+guVysHQ9vZ2DAYDFRUVCr/AYrEozU7Hg9Fq8SVJWgFcC9iBbuCqwSW7Q+GguvpDGWp/DfqsrKxJU2RlWK1W8vLycHd3H7HAZqwG238SGa4N1ngMX44NZGVlDaCbDl41zWYzjY2NB7QiT0a3xXk9vv2aZPq4Ox+hjlFy+UNV4o13W3AoGP5oUf3BXo/FYsFgMPCf//yHX375hZtvvpnTTz+dyy+/fNT9/hhr8d8VQrzUd/zpOAt2ThztOg6qqz94j2+1WsnPz8fFxWXcGvQjrfjyShofH09ISMiI44zlwbLZbOTl5eHq6jrsJDLWrUf/slyZFty/OnDwqrllyxZMJpNCve0v1LG/3eCePlff23XvYyNPAl3jJPFMZFtwKLj6403n6XQ6QkJC+Ne//kVOTg7//Oc/2bFjx1jHGLUWXwjRv5WRB3vLKUbEIbPiy9VvcXFxE2oqONSK3z8wOJ4Cm5EwVirvWFZ8OcAYGBg4prJcWUpc1reTaalNTU2UlJTg6upKQEDAuNloY4WssNXf8D1dnP/fNUna7lDbAoPBMGBbYLPZJsWVnwpMZvIxGo1kZGQwd+7csX5k1Fp8AEmSrgVuAXTAsWMZ+KDu8e12+wDXfjLGOXgikSv1VCrVqIHBsULW3x8LlXc0I5a9kKH66I3Vpe1PSxVCYDQa0ev1SttseQ/t5+c3JdcvU3a9+hm+u85pBN3jUOEZDf23BVFRUYqAZ2NjIzt37pySBqATxWQIPFardb9sz4QQzwPPS5J0IXA3cOlonzmorr7FYmHXrl24ublNur1Uf8OXRTimqiBGzvf39PQo+vujYaSYQ319PVVVVVPmhYDTWNzd3XF3dycyMnKA2m1FRYWyhw4ICJhwfb65b17tv+K794lx7M9CHVnA09XVlblz5+7TANTDw0PZFuzvuvwDvN0YrRZ/MFYBL45l4INm+N3d3bS3t5ORkTHqvnsskA1Nls6eTIFNf5jNZoXKO3PmzDEbzFArvtx222QyMXfu3P2aqhusdms2m9Hr9QMKceRtwVhr1y1DGb5u/xv+YAzuhSfX5RcWFiq6fXK2YKrv8URXfPlZGOeEO2Itft94CUKIkr6XpwAljAEH3NUXQlBZWUlTUxPu7u5TYvTyd9TV1SnNKafiB59MW+vBhm+xWMjJycHf33/K2m6PBy4uLvsU4uj1emprawEU13k47Xur3YFcnNff1ffoM3zjATT8/hhclz+Urr88AU7FtmAyMYbxfvcYa/GvkyTpOMAKGBiDmw8HcMWXJAmz2UxeXh4eHh7MmzePzZs3T8nYRqORyspKPD09ycrKmvSPK6vu1NfXT9gd7+/qy7LZY5lA5JqG/Zm66l+IExsbq2jbNTQ0UFxcjJubm7ItkF1neUVXSeCm3TsxyCu+8QCo8IwlPTqUrn9bW9uUbQsmytW3Wq0T9RRGrMUXQtw47kE5gIYvEzfGqlk3VrS0tFBcXKywqCZrLHa7HZPJRHt7+7hVd/pDNl5Z2jsrK+uQ4WkPRn9tOyEEvb29CjHJarXi6+uLWeMknrhp1QPusVuf4R8oFZ7xQu77N9S2wGKxDCARjcVLnOge/1AqyYUDZPhCCBoaGpg1a9aUpZmEEJSWltLR0cHcuXPp7OxEr9dPasz+rbKTk5MnHQlva2vDaDRO2dbjQECSJEXFVnad29vb2VHeBIBWclBdXU1AQADu7u6oJAm1SsJi3/+GP9lJfaq2BRM5j0NJhAMO4B5fruaaCsgBNz8/P2bPnj0pbryM1tZWioqKSE1NpaSkZFLnarFYKC4uRpKkKdl6HEzIRBv/HjXQiq+HCxqNhvLycnp7e/H29kYjgXU/r/j7o+BpLNuCiaj0DIXf5YoPQ0e5ZZWa8bhOMqkjMTFRKZuEidX3w94GGW1tbQr1djKTiKwMFBYWRk9Pz/+00fdHT18uz9tVqwQJ5QIVjaoVo02wY8cOxZC8vb2n9NoPBF138Lagt7eXtrY2hRfh4+OjSMeN14P7Xa74w2E8NfRyNqC5uXnILcNEjFXm73t4eDB79mzlPCZq+DLBJyMjQxHs3B84GJz1LrOsp79vgYqbTkOP1UpaejrtfQ0jCwsLlUBaQEDApFtDH+hr7r/lkVV6Ojo6aGlpITs7Wym+GeskdyjV4sMhYvij5ZFlDr9M4Bhqohjvii+vzNOnT9+HIjxewxdCKB1zZYJPd3f3hNxTvV5PQUGBIvQgV+XJD9bB8iDkFd/Xfd/fykWrBqzYURMSEkJISMiAQNpg7b6JyF0fbJ6+rNLj4uLCnDlzlOIbeZJzd3dXvJ2htgW/2xV/JNruSJCVd4Yy0MFjjdVYZZLPYM0+GeMxfKvVSm5uLt7e3gMIPmPtndcf/Rt3AEoPuZ6eHry9vZVGEgdlxe9T4Qhw29fw5fRer8WupPeGCqTJctelpaVKu2x/f/8xNQA9FCrz+k8+cvGNPMkNtS2Q6dIajWZCe/wxlOTeAlwJ2IAW4E9CiKqxjH1QV/yhKvRkyHr5NTU1IyrvyBhLNVx/5txIkfaxdrqVvYa4uLh9iEjjKct1OBzs2bMHh8PB3LlzsdlsCCEGFK3IWQu5caZc2TYVyjZjQbvRWaHj5zGU4cu5/OHvf3/tPtjbLrusrExRuQ0ICBhWfORgr/gjncNw2wJ54s7OzmbDhg1ERUWNOT7Q9yyPVpK7C5gjhOiVJOka4FHgvLFcyyHh6g+G3OVWCDHmApvRvIf+lXCjMefGsuL3389PxmuwWCxkZ2cTHBxMdHT0kOclSZIi9mAwGEhJSaGjo0MpYZ0I/Xa8aOsrzfMZYsWfCG23f7tsWeVWr9cPEB/pn1Y7FFb8sdJ1Bzf3iIqKIj8/n507dzJr1iwefvhhTj755BHH2Lp1K4xekruu30c2AxeN9VoOOVdf7jgz3i63I62wciZguF72gzGS0cr8gc7OzhELdsay4suNO8ZDCZYkCa1WOyD6LBuNTL/dH95Ae5/he7vue72eLpPj6/dXuYV902qenp54eHgcdA3DiXodwcHBhIeHs2zZMs4777wxxaLq6upgDCW5/XAF8PVYz+mQWvHlVTQtLW3cHWeGU2yR983jIQ8NZ/jyft7Ly2tUrfzRDL+pqYmysjIyMzP38Rjkz64r1lPe0oPBaGNxvB+zIvfVZO/vDcBeo5lqb0CO6nu57rviefXV5PdapqY0d3Barbu7m4aGBjo7O9m+fbvSOPNAiI/0x2S76MhVkVNN5pIk6SJgDrBkrJ85JPb4g6vWpsJdtdvtSnutiTS8HGz4sgDHaEFGGSN1yy0vL8dgMAx7rUIIXv21hmfXVxIb6E6NwcjKLbW469RM94bT7Y0cPSOQMJ99ueaDjWZwMY6cKRhvwUqXaW8efzDkFb/dOPWFOnJtvhACu91OQkICBoNB6YXg6uo6IEi4PzGZAp3xRvX7yslHLcntK9C5C1gihDCPdfyD7uqbTCa2bdtGcHDwlFWt9fb2kpOTQ2Rk5IQaGAw22ubmZqX33VhLfYda8eWGIK6urgrjcDDMVjv/93khX+Y3c1JqEA+cMgOHgM0VBt7cXMuu2k7yvyvnoe/KiQ9yZ1GcP4vj/JkZ6Y1WPXD1G6oYR54Eurq6lB55AQEBo062shvfvyRXhrzvl7cD+wNy37zBPfHkaLosde3j40NAQMB+K8k9UA0z+1R6RivJnQm8DJwohGgez/kc1BW/t7eXuro6srKyBijUTgZy0U5aWtqE2xXJhj+4HmA8nshgwzcajWRnZ48oDtLSZeba93LIqeskLtCN+dE+dJps+Llr+aXcwK7aTmaFaLj5+CRyGnr5ubSN/26p5c3NtbhrVSyI9WNxvD+L4vwJ9d6XMDM4NiB7A7m5ucDQ3sA3u5t57IdyWrqcKrpDGb78Xodx/xn+cPtrWXxELtLq6OhQdAdUKpWyzZmKeMdk22eNJ53XN2mNVpL7GOAJfNh3bdVCiNPHNP74Tn9qIISgrKyMlpYWQkNDp8TohRCYzWYqKysn3aBSpVJhtVrZtWuXwuob70PT3/Dl4OJgCe7+2N3QyTXv5tButHJErB9lrT3c91UpUIq7TkWvxcEpqUGcHWkkJdSTjEg/OnqtbKlsJyHIndRpXmwqN7C22FmoFBPgxrEzAlkU50dWxNi8gcEBtYCAANYXtWG02BUFRw+XfR98WW67awrltwZjLFH9wdH0wQKesoczUYXiyTbMHK8wzBhKco+b0MlwEFx9i8VCbm4uPj4+JCUl0dw8Lg9lSMhBN4CZM2dO2sWzWq3U1NSQmJg4Jk3/oSBfr9zgczgdf4Bvdzdx+ycF+LppWXXFXOIDXbHb7fxU2sa9XxbTYbQhAWsKWlhfLLGgvJjGLgsFDd2cnRXKXSfGo1U7vZSn11fy+qZaukw2/rulltc31eDpomZ2pA/HJgZw5HR/QobxBvoTUrq7u52luXVtRHgI9vTtHocyPV+3PsFN0/4z/IkY3WABT7kLcH+F4oCAgGHFRwZjssG9Q6ml2gFd8eWVT05fdXR0TKpVNuxl9sXFxVFdXT3plE9zczM1NTWEhIRM2OjB+aCaTCZaW1uHDS4KIXj+pwqeXVdOVoQPz52fQZCXC1arcyW/64tiNCqJNy/JJC7Qnc2V7XywqYT1pQbsDud1FjZ288ov1SyM9eOL/GY+2NnA0sQAHjo9CYcQbKlo593tdfxU2sZPpW0AJAZ7KFuCzAhvNKqB5iwH1Dw9PWk213NsvB972loA2LJly4BMgU6nw2OKlHZHwmTz+P09nP4KxY2NjQPER0aqxJvMit/b2/v75OqbzWZKSkoGpNUmWlEnQxatlJl9tbW1Y+L+D4X+0fb4+HjM5jEHSPeB1WolJycHSZLIzMwc8oE1We3c+WkBXxU0kxrmxf2nJSndZj/YUc8/vykhNtCdZ89NJcLXeb8ifV0pMdjRqSWuOzaWXoudjWUGXvy5mhd/rgYgPsid45ICsdodeLtqKNf3srWqg5kR3txwTCw5tZ1sLGvjjU01vPqr0xtYON1PCRL273jbbrTRZbJhtO+dTONSZ+KKBb1er6yc3WrnStZ7kF398WCwQnF/yq3FYhlAuZUn7Yk+W/JnDyVNhgN2Ji4uLsydO3fAjzcSZXckDG5CId/QiU4kcoMMNzc3Zs2aRWtrK0ajcdzjwN60X1xcHGVlZUM+rE2dZv7yXg4F9Z0EeGgpaOjijBe3EOHriqeLhsKmbo6I9eXxM5PwdnMa4tqiVu74vBB3Nbx6fgoZUc4a8tPSQ7hmVR5VbUZSw7yoMRi58/MiJJx777ZeK4vi/Hjq7BRctGrmRPkQ6u1CTl0n/m5a5kR5s726k+/2tAIQG+DG8clBLIn3x9rnVch/A6jrMJMV4Y2XlxcxMTFYrVYKa5qBJpranN6XTCCaSinp/UnZHUy5lQU69Hq9olAcEBCAyWSaUF3+wSYeDYUDuscfbAQTMVSTyUROTg4hISFKEwoZEymnlZmC0dHRSoOMiZbltrS0UFJSoqT9ysrK9jkmr66Tv7yXQ7fZxvMXZLI0KYhag5Fvdzfx6sYqatudnVh3Vndwx+dFLIzxprnbyptb6kmb5sU1aSoSgpwPX05dJzd8WIDVLnj5gnQWxPphdwi2VBq4/6sS6jucXsvGMgMnvbCNhbG+mO2Cb3a3MCfKh6fOScHHTYsQgtW5Tfzjm1IaO838Z2M1L2+sHpKsU9duIitiL7lKq9USFRYMlKF2cSM6OnqAN9A/UzAZwz2QXP3hFIoNBgN6vR69Xq/8fawewFDP/8HEAfU9Bqe4xmv4cpfb5ORk5Ufpj/GOJ6f+0tPTBzAFJ1KWW1lZqfTRG26l+yq/iTs+LSDQU8d7V8whKdTpIjuE4OOd9XSabNx1QjyR/m5sKG3jp+JWfi4zAE4124xpntR3dWG22Fhf3sQ9X5YQ7KXjuXPTmB7oJK/Ud5h4+LsymrssPHR6IkfE+rGxrI0NpW2s2d2CzS6Q+r7zk+xGFsf7U9bSwwNflxDm48qL56fh5aLh57I2XtxQqRB3ZNQYeve5Lpmrb7I6lEYY8j66ra2N+vp6Ojs78fDwUGID463PP5hcfVmhuLu7m+DgYCRJUjIggBIkHGlyO9RW/YO66RhrFZzcbaepqWnE6PhY+9XJDS/1ev2Qqb/xGL7MENRoNAPEPPrD4RA8s66cFzdUkBrmxSsXzSSgby+9rdLAdatycQjBfy5MZ36sP5IkkRHuQ0lzN3UdZo6I9cPucLBqRwN2AY9uzcZiF0T5ufL02UnE+DvvR3ZtJzd+WIDN4RxrbrQvAMcmBvJlfjM2u+CMjBACPLT8UmbgibUVPLG2AoBATx1/WzadCF9XJEliaWIg//q2FIAIX1fFE/lgZyN/yAwl1Hvvb6BVq5CkfavzNBrNABHPnp4eRW9gvFF1mcBzMCFH9b28vPDx8RmgUFxfX09XV9eAunz5ObVYLAe8weloOHSiDcNAZrvpdLphRThkjKUmXx7PxcVlWEMdq+HL246wsDCioqKGPMZosfO3Twv4drczbVnQ0MWFr29nSUIAapXEfzfXEOnnynPnpRET4ORyV+p7WfFeLg0dJh7+QzKnpYdQVlZGXSK8Vgjbqztw0aioNpg489VsUkM9CPV2YUOpgVBvF547N4XpQU6WWH2Hib+syqeqzcg/TpvBGRlOuvENRwvu+6qYz3KaCPDQ0tFr4S/vFxDp58r5s6bxY3ErXX3iG/NjfKnNbgSgtdvCea/v4qmzU5jZr3ZAq5IwjSCx3b8+Pzo6ekBUvaioaFS1nsHtqQ8GhkrnjaZQ7OXlpfBBxoNvvvmGk046qYjha/GPAp4CMoDzhRAfjWf8g+rqjwZ5/z3WVlijufoylXe08cZi+HLvu+G2HQBtJgfLX9/O7sYu/nZ8AsclB7GhRM+6ohbe2lKDQzh16qP83dhZ04Wnq5ayll5u+jAfjUrijYuzyAz3Ii8vD4uk5dVC2FHdwfVHx3L1oigKm3r4qbiVD3fVU9DYAzhJNK/+WsOR033x99Bx5+piLDYHL12QxvwYJ7HFaLVzx2eFrC3Wc/G8cG47bjpGi50filpZtb2Bx34sB5w5e51GNaDqztdNjdFi58p3cnntokxlv69Vq8YlsT04qj5YrUdeNeVCnEO5Hl/GUArFNTU1fP/99+Tn53PSSSdx+eWXc+655474PXa7nWuvvRbgJIavxa8GLgNum8i1HLIrvsyPH0+l3kgGK6vojoXKO5rhj6X3XW5tBw9uNmPFyksXZnH0DKcAxR+ydPxSpsch4Kh4f4K9dGwsM/BTSaHyWT93LfedmkhysCs7d+7E5h7A/euaaeww89hZKZyS5hT9iA/24J1ttTR2WliWHMjRCQH8WmZgQ1kbX+Q78+46tcQfZ4bi6+KcFDtMdq7/oIC8+i7uOD6O5XOdE6CHi4bjk4NYnduEJIGXi5pOkx1fdy36Hotyboa+QpwADy3Xvp/Pfy/JJC7IA1etih6LHavdsQ9LcDQMVusZ3AXYzc0NIcRBd5fHS+CROxv/9a9/RaVS8dhjj42JsLZ161bi4+MpKysbqRa/su9vE1KFPeiGP1hpVwhBSUkJXV1dIwbKhsJQK74ceGtpaVFUdEfDSJV1xcXF9Pb2jtj77pNd9fz9iz346OCty+cwI8QZxKtrN3LNuzmUNvdw14nxXDg3AkmSsNkd3P1FIatzm/ByUdPea+WGD/Lx1kFiiCf5jbW4alW8cUmW4l53mqzc+EE+WyrbuWZxNNcdHYskSZyZNY2Vm2t49LtSAj11eLuqeWd7A+9sbyDEU0uP1YHZ5uDxMxM5LmmvBoDJaueGDwvYVtXBP09PRKuCv35WREuXGZ16YFAtLcyTR89M4uKVOdz08W5WXT4Tb1cN+h4rHUbbAC7ARDBUjr24uJja2lrq6uqU2MCBLsudbDONqKioYbeE/VFXV0dkZP/CvFFr8ceNg274/ZV2+9N5R6t3HwoqlWoAL8But5Ofn49Wq2XOnDlj/tGGMnybzUZOTg7e3t7DauU7HIKn15Xx0oZKADy1EmsLWxBI9JhtXLcqB4vNwYsXpHFkXACSJNFjsXH7J7tZV6znwrnh3HFCPJV1zXyxvZRNLTq21TiVei02O0+vLeeohAASQzz41zelVBuM/OuMZM7IdO7b7Q7Bw9+V8M7WOo5LCuSRM1Nw06pp7DTx1pZa3tpSi60vN//8hmo6TTZOSQ3C7hDc+mkRWyraeeDUGZyWHsLT65xBP29XDXV9gT03DRhtcPbMUCL93HnszGSufCeXe9cUE+iho0JvpMNonbTh94fsPss8e29v7wFluUO1+9pfmIzhH0pCm3AQ9viDIRt+b28v+fn5k2qxpVarsVicbulkSnMHG35PTw85OTnExsYOS+PtMdv426cFfL+nhZNSg4n2d+ebnGqeXFvOk2ude2ZPFzW3LJ3O7ChfJEmiocPEX1blUdrcwz0nzeCCueHU1NTQ1tiI3SOY/Lw6FsT6cumCSHZUdfBTqZ7Hf3ByA1SSM1of4KHFbLNjd8BfPylgXbGeSxdEcttxcaj7qLjZtZ28s7WOcF9Xnjg7ld2N3by9tZb7vy7js9xmJAlyaru458Q4Tk72x263U9rijBncdWI8f/20EJW0t022h8752MyN9uXGY2J5cm0FqWHOB7u120Jc0NgCWfUdJtYX66nQ92JzCNx1aqL83EgK9SQ5xBOdZq+Ryem8obyB/sG0/e0NTCSl2N3dPe5a/Jqa/uI7o7bHHjcOiRW/rq6O5ubmSfeLlycR+UFITU1V5JzGg/6GL481ONffH/XtJq55L5vipm7+78QZXLIgEkmSOMJLz3ctXryzrR5fNw0mq50Hvirh0e/KSAnzpLi5ByHgxQvSWTjdj8LCQnpMVt6pcOX7wjr+OGsad5+UgFatYklCIMlhXtzx2W48XTTEB3nwc2kbPxS24qqR0KrVdJtt3HBMLCsWxwBOY3ljUw2P/1DGrEgfnj0vDT93HclhXpyVFcr7O+r5x9fFOATceXwcf5wdjt1ux+FwKKv8kjhfVBI4+sVk5b8BXLYggp3VHWwoc9YBVBmMzI8dudqytt3ICxuq+DKvGQGoJfB109JjsWPqCxBq1RKzo3w4ItaPI6f7IYaohR8qmDZYpENmEe5vb2AkjFdhd+7cuZSUlDBSLf5kcVANX+7EYrPZxiyqORIkScJgMGAwGMa8nx9uHIfDQVVVFY2NjSOOtaumnWvfy8Vss/Py8iyOSnAG8Xotdp7baWRbYxdnZ4Vy98kzANhW2c5/t9Sysc9QAP79QxlrfGxkhnvzSZGZgoYu/rosjsv6JhAhBK/8Us1Ta8sHGLDJaueT7Eae+LGMbrMNATyzroLv97SwJCGAitZevt3TwgkpQTz8h2RcNHvvb0OnmXe31aGSJDQq+HZPK+fPjVBiKvoeKzq1hEYlKUbvqYVuKxQ17W0UopIk/u+EOH5+oQ0B7GkcvomI6CMNPfK902vJivBmV20nD5w6g9MzQrHZHdz1RRFfFbSQHOpJa7eFJ9dW8OTaCrx0EounWzk2JZQjYv2G1AXor+Tbn3/fv/mnLNJxIGMD43X1NRoNzz33HKeccsqwtfiSJM0FPgX8gNMkSbpfCJE65u8Y5zVMCv3dJKPRSG5uLjqdjtjY2Ekbvd1up7KyEovFwsKFCyf1wwohMJlMiqDmcGN9ltPA/322G1eNipuXxjE7yheApk4TK97NYU+jjavnBbLimOm4aNQIIciv72JjWRszI725ZWkc28r1fJtXy+oSB5+XOCPx82N8CfNxpdtsx0Wr4oE1xXyS3cApaSH84/RExYB3Vnfw1NoyPHRqVl6ShU6j5qeSVn4sbOWln53y6l4uGuIC3ek02gjycn5uV00HN3yQh9nm4D/LM9H3WPjrJ7t5/ddqru7zFrrMNnzctDhUeympssBOVZsRq9WKSqVCkiRy67uVyeHrghb+dEQkkX4DOe2GXis3fVjAztpO5kT5cOXCSK59P5+TUoI4ItaPHoudf35TwlcFLfx5YSTXHx2DJEk0d5l5Zn0ln+c28WNJO18VGlBLkBnhzaLpfhyVEMCMYI99XPCh+Pf9df1lb2Aq+uKNhp6enn3k10fDySefjBBiRv/3BtXib8O5BZgQDsqKL7vPKSkptLS0TKrZJTgnkZycHPz9/bFYLJMyelmGW6VSkZaWNuSezu4QPPljKa9srCLC15V2o5V/fF3MI9+VkBLqRVlrDw4H/PsPCSR4OMezCYl3igXrKno4PSOEB05NwtjThcWlmZijYrn7q0o0aomMcG9y6zq55aMC1JIzzdZpsnHBnHDuPilBOZ+PdzVw/5oipge688IFGUzr09/zcdPwdUEzaglOSQuhpdvCCxuq+M/Gai4/IhI/dy1P/FjONF9XVl6azvRApwv6ZV4T/91SyyULItGoJKx2QbCnjvJ6p7CHp05Nt8WOSoLmbotClnI4HPywpxlfNzXtRjsmq51zX9vJDUfHcsyMAMw2Bz8WtfLm5lrajVZctSqePy+V5W9mo5Ik1hW38vXuFuXexvi74eWqobCph8QQD74qaObz3CYyg7X8++xU6noEG0vb+LGolWd+quKZn6oI9tSxKM6fRfF+LIjxw2sUbwBQYgNFRUVj8gYmQ7k91BpmwgE2fJkq2z+1ZjAYJlWaK/P3U1NT0Wg0QxbGjBWdnZ3k5eWRmJiodLsdjG6zjds+zmddUSvnzwnn7pMTAecq+tovlfxUrFfUal78pZ4lMwKYFRnLSxuqyK3v5twkN04I6SA/Nxuj0UixCOPfq8tIDPHkhfMzCPF2weZw8P3uFh742inCAfDe9jp+LW9jSbw/+l4ra/KbOXK6H0/+MQ3Pvnr44qZuVryXS6fJxvPnZ3BUglNKvKqtl+fWV/DKL87S3YXT/Xj87FR8+2nkX3lkFBe/uYtPdjUwJ8qZMgzx1PCfdU5+wawoH9QqifXFejqMNqU4xWSxsbG8nWVJAXye28w0Hxd83bU89G0pD/VRfgGOiPVlV00np6YFc88XxZS29OKmVXFGRgg7azopbu4hJdSDXqujj0ZcgatWhcnqICPci6tTJbzdtIT4u+HpouGTnEZ8XDVcMHcaZS29fF/Ywic5jagkmBnhzeL4ABbF+Q3pDcBeyS7ZG2hvbx/Q5UeODcjewGTVd37XUf3W1lZMJtOA1JparZ5Qaa4snd3Y2Kjw941G44S9B1naOysrCw8PD4qLi/c5ptZg5Jp3sylr7eWekxNZPi9C2YNvrTSwvljPzAhvbjtuOjl1XfxUoueNX6t5tW8mmBPlw4LUcNylDno62vikSs3XpVXMDtVy3wlh+Lk6H9CC+i7++U0JQsDKS7KY5uvK+mI964paeXtbHQ4BGpWEu4uaHwtbWRzvT2FjNzd+mI+ni4a3L5tFUujeB63bZCe3rguV5HSBO41WtIPEN2ZH+TIr0ofXN1Wjkpy5Zg97Nz/18U1mRngzK8qXdcV6bA5Bj9mGh4uG7dUd9FjsHJ8SwlcFLRitDr64JJOn1lXw5pZ6rj8qkkVx/lQZTGyqaMfHRcNHuxrxc9PywRUz+ee3pRQ39/C3ZXFcNM9JJmrqNHHXF8VsqWxHp5bIrevi2jpIz99DYogXXxU04+Gi5vWLMpXipF6LjWvfL2B7dQeNnWaeWlfBU+sqCPbSsSDGl2MTA1kQ46uIhvSH3Apc7rkwuB2Wn5/fpLIEv3vDDw4O3kdzbiKluYOls/t3uR3vWLL+32iCmtur2rn6nWxA8J/lWSyKdz4kJqudOz/bzVf5TZyeHsL9pybiolUzM8qPuCAP8uu7UKkkssK9ya3v5JaPd6MCvN00tBttnJkZyu3HRNCmbyUnJ4etDTZeyTMT7KnjpQszmd6XGjslLYSvC5pxCDgtIwSdSmJDaRvf96uVD/TU8dhZyYrRCyF4b3sdj3xXSoCHjv9eOpMOo40bPsjnke9KeeC0pAHXeMmCCG76sIBv85yZI1efALoqmgCI9HdnTrQv/u4a2nptbKls59jEQL4qaMZdp2ZBrB9uOjU9Fjs6nY6Cxl7ig9y58shohBC8tLEaP3cNb2ypRQKePTeZ+74q4ZdyA3edEM/5c5wl0Va7c8XfUtnO5UdEcMOSGAoau/nol93sarXx4a4GAHRqFW9sqmFxvD+zIry5Z00x26s7uO246Vw6P4KmTjO/lLfx8a5GVuc1szqvGY1KYlakt6I+FBc4dM++/gKesjfQ3NxMZ2cnOTk5SmxgrBmo8SrsHggc9HSeRqNRcu9jgclkIjs7m2nTpu3DghrvJCILcLi7u49IGPp4Vz33rt6D3SEQwJM/lpFd20FWhDdP/VhGfn0XNx8byxULo5RJ6J2ttfzr2xISgj154fx0AtxU7MrOoczswfNb2jD0Rco+zWkkp66To+L96bV48UF2A2khbtwwy5Wm0lzMrb6Ydd7c+XUNDR1mHj8rhZP7KLt2h4P7vizm4+wGPHRqWrstXLIym8QQT05JDSavvovvC50R/n+dkax0ur1oXjhvba3lonkRzAjZ+0AeFe+Ph1aiuMUpQvJjsZ4ADy36HivR/m5IksTxKcGs2l7PhzvryQj3Zk1+E3+cNQ1XrRovFw31JhP6Hgs7qtu5enEMWq2WTpOVjeXteOrUOAQsS/Tnga9KKWnp5c5l0zl3lpOAZLLaue2TPfxU2saNx8Rw5ULn75sZ7s0ObzVflZuICXDjvFlhZNd28WOxns9ym5TzX5oYwBGxfgghCPF2wU2rpqChi7RpXlx9ZCQ7azvZWGbg3z9W8O8fKwjx0rE43qk8tCDWTykvHvxMyelAh8PB9OnT0ev1lJSUYDabB8QGhgtQ/+73+MMReMbq6sv7+eHUasdTTjsWuWub3cG/fyjl9V+rWTjdn2uPjmVbZTvrilt4dt3eQpaFcX4khnhidQjUOHj421Le3VbHMTMCefSsZLCa2bkzG6vXNF78uQabXfDq8kwi/d2cLnxxCys31yJwuvARAZ50uAWxMNaX3Mpmbv+wFCEc3LnAiwxfK2azGUmt5e4vCvkyr4mzZ4Zx78kzaOm28GNhC5/sauCJPtLQ3GhfHjkzeUAjjKuPiuHTnEYe/6GM/yzPBJx72LLiIuZHuLKuwmn47b02UsI80fdYifJ3rm6npIWwans9v5S18dLPlVjtglPSQnh2XTlWuwOHgM9zGp0G3kcJ/n5PC1a7wNAXr/i+aG8q8/G1FawpaCYr3IvNVR2UNPdy1wlxnD9n72+yJr+Zp3caSQ715KXz0/F113LRPND3WPjT27lUtPYS5u3Cj0V6fizSE+rtQrivCzurO8mK8ObF89PwcNFw9IxALptv4fK3c6lqMxLp58pXBS18tKsRrVoiJdSTZUlBLI73JzbAbcDzKu/x3dzciIiIULwBWamnvLx8QGvz/t6ArPB7KOGgr/hjWaWFENTU1NDQ0DBqPf5Yoq+y6OdIBB+THf7yXg4/lehZPi+CO0+cgVatYk60HwnBHtz6UT46jYrEEA92VnfwS5kBV40KDxc1+h4r582ext0nzaDd0EZJSQnNumnc93kFgZ463rgki/g+F/6MzBDWFbcigOOTA3HTqvm5tI1vdrcoira+7loeOTONrBAX9Ho9W3bl8eS2HooMDlYsnMZ1xyagUqmY5uNKQrAnjV1mPF3UJId6sa2qnT++soMnz0klJcz58Pm6ablqUTSP/1BGXl0nySHu5OXl4efnx4ULQ1lb4VQsvmR+BBv6BDrlbjmyXr/NIXh3Wx2Bnjr+/E4OJqtdSek9/kMZkX5uJIY4r/GT7AblvurUEnYB1y2JIcLPjcLGbjaX6/nvNucxrhqJrZXtaFSCo+IDWFvcxkPfljHDT8UrF6Th1ReQbO22cPV7edQYjDz9x1SOmRGguPfvbKtjR3UnAHn1ndz08W4WxfmTGurJP74ppa7dxPPnpnJknD9Wu4OdNR28vLGabVUd5NR18fiP5UzzcWFxnHNLMC/Gd8gCncFKPUajEb1erzT38PX1Ra/XYzQax234YyjLdQH+C8wG9MB5cuHOWHDIG77D4VA6586ZM2fS+f6amhrq6upGnEBq2nr5x2YTjT293HdqEhfMdaZLhRC8srGKf/9QSvo0L549N41gb1fMNjtr8pp59PtS9D1OF/6DHfXkVLeR4mPHxSeA93aUkRXhzbPnpRPg4STJ1BqMrHgvl+o2Iw+dkcQfMp10YLvDwUPflPDe9nrctCoMvVaueieXKH835kT5sKnCQkuX4P+ODSfL18KWLVvw8fFhU7OKpzc2Eh3gxvPnZRDl78aO6nZu+3g3F76+k5cvzFBYdX+cNY3n1lfw4Y5aTg3tJjIykrCwMCS9U2FHAi49IpK3ttbiqlEpq1+Ql9Pw50b7sKO6A7UExyUFsmJxDJ/nNPDyxmokCc7KClVoybtqnEaolkCjVvHCuWksnO40llmRPqwrbsVFo+KSBRHou638XKrn+yI94PRaovxcOX26DVet06Nr7rbw53fyaOoy8/x5aRzRd00h3i4Yem0UN/dydII/F8yZxqaKdjaWtfH4D3s9tKNn+GN3CIxWe992oJttVR0sSwrkpmNi2FzZzsYyA5/nNvL+zgY0Ksia5kF6oJozA3qJ8Xcb0nvt7w04HA7a29tZuXIlxcXFnHvuuZx88slcd911owYJx1iWewVgEELES5J0PvAIY2yRDYeAqz+S4KYsdBEaGkpUVNSkpJdkgU6r1TpiL72tlQb+8m4OPRYHp6cHMzPSByEEVrvgntV7+CyngZNTg3nwtETc+jjr+fVd/PtHZxrx9Ysz8XbV8PGmYrbU9PBJqQNoxE2rIjHEk8LGbuZG+7KnsYtrV+VhcwheuyiTuX218la7g398XcyHOxs4JS2Yf5yehL7HyoYSPatzG/mkTxDDTatmZ6MVH59AjsxI4Pm1pby3q4H0QBXXz9Ig9bRidAtkdpQvH/15Nsvf3MmK93K5fVk8xyQGEOLlwnGJ/nyZ18RV89IICw2ixmDksrd2ASCAJ34sxyEYUHSjU6vwd9cSG+jByktnKe83dpqUSW96gBvPra9ECPgyv1E5xstVwyvLs0id5lz9KvW9XPF2Nl0mG68sz2ROn2KQw+Hgni+L+DS7EV83DdUGE4/tgFf3bCVzmhfZdV1Y7YIXz09TSFNCCJ7fUMXLG6s5KTWIf56WiFatYuF0f86bHcZl/82h3WgjJdSDzRXtrCtuQ6eWCPV2odpgYkmCP4+emYxGJRHl787iOH8uf7sHfY+FpYkB5Nd28EZtD29kbyfIU8eqP80k2Gt4Zqjc6vuf//wnP/30Ey+99BKbNm0aU2ZgLGW5fa/v6/v/j4DnJEmSxBgJB4fsii+74yMJXYwVFotFicYOFujsjw931HHfl4UEeerw1Tn4LLeZz3KbCfFywSEELd0WrlkczbVLYpQfcHVuI/d8Ucg0H1devCCDCB8d+fn5nJroRVm3mvKOdpYkBCDh3Pu+v6MenUaFze7A103LSxekkxXpCzg5Ajd/mM8v5QauWhTNDcfEopIkpvmoCfFyoaipm2k+Lly+MJLCxh7WF+v5Zvfe+u450T48cmYKvlpBa2srBbt3s7HaxPc1gup2p1E++HUxD34NUb4uhLpYMdnh7Zx2XHd3sXJzDY4+f91dq+LLvCbctCoiBrHwQrxdaOjYy9fvMds46+XttPe10PJx05EZoeWZ9RXKMRqVxDEzAmnoNBHh50p9h4mr3slBCHjzkpnKNsQhBP/8ppRPsxv546wwlidqaG7vpk0Xwpr8Zn4qdVKDVTjpyYvi/Fg03Zcv8lt5a1s9Z2WGcu/JCUqBUqW+lyvfycVodbDy4kxSp3lhtjnYXtXOcz9Vkd/QBcBPJW2c+fJ2FsX7kxziyfMbKuky2Xh9eQbp4d40NjZS09bDf/ONrC9pY3dD14iGL0O2w+joaGJiYkY9HsZclhtOXxttIYRNkqQOIABoZQw44IY/FsFN2R0fT2vr4SDLXcfHxw9b9WezO3jkuxL+u7mGxfEBPHFOGmWF+fiHx/J5vp7XfqlSikf+u6WWcn0vSxICKGrqYeXmGuZG+/L0uWm4SnZ27tyJyjuEe39soK7dKZ11ep/cldFi48GvS/gspxGNSqKt18ryN3aRFeHN7Chfvi9soa7dxIOnJXL2zGnK+b27rZaHvikhJcyL589PJ8jT+cC1dpv501vZlLb0EuChZXtVB8ue3sS8GGcPvbVFEjuqLUT7uXB5phc2s5G3Cm3M8NeiwcrOFqFcE8DxyUE0dJjIq+/ikgWRvPRzFQ4BsQED01ZxQR5sr2pXXn+S3Ui70crTf0zjlo8K2FnTMeB4leRs4vF1QTOf5jQiAZLkFOn85+nJJPelH20OB3d9XsgXeU1cfkQkp0XasFktHDknk4KGLnLrOwn01HHbsjjKWnrZUKLnmZ+qeeYnJzEpLtCNI2O9MZqtuLtoKNcb+fM7eTiE4PWLMkjsy2Do1BJbKtvJb+jirMxQLj8igl/KDfxc2sb72+uQFcRmR3qzp7HbqY/ocLC+opf1JQZOTQtmcV86dyw41BR24RBb8R0OB3v27MFut4+7tXV/yCWcY+ly22m0cvNH+WwsdZaz/u2EGahVEiqVit0NXbyxqRovVw0vnZVMt9nOTyV61ha18m0fzTTQQ8vieH9qmw101Zdh8org/76oRELi9YuzFFfUanfw0LelfJbTyMmpwfzj9ETKWox9lNUmXv3V+fAGeWopaurh1/I2ZkX58MzaCt7cXMMxMwJ47KxUJeVUqe/l6ndzaOmy8My5aRyXFERpSw9fFzTzwY56Nlc41XkXTvfjvlMTifB1qtiYP9nJBwWd3DhTx63zXPi8UmJ1YSf/vTSLOdF+HPHoz+jUEiuOiiHAQ6c09uiPxBBPvsxrot1oxctFw9tba5kZ4c2y5CAi/VyobDMNOP7NS2YyJ9oXi93BW1tqeHptBSoJus12bvwwnyBPHUfGOdV+8+q7uOHoWBYH9CBJGpKTk9lZ08GKd3Pxddfy2kVZRPk7F4Prjo7hrx/v5ts9LUT7u1HfYebWz4rRqSWSQjwoaenFVaPitYsySAjey2149Pty3t5Wx3mzwvi/E+NRSRIxAbJ7n0OP2c6iOH/y6zt58Bsn+9DHVUWHycFp6cE8eGqi4lGM5VkcL913jGW5dTjbaNdKkqQBfHAG+caEg274ciS+v15+dHT0hGdIOT1YU1ODXq8fUcWnUt/L5St30tRl5sHTkjm3L4UkhGBNqZH/5haRHOrJM+emMa2vm036NG92N3TT1mNlfowvbb0WnujTqAv00NLWW0qIlwuvXZRFTJ/BdJqcBSqbKwys6FPLUUkSqdO8aOk288amGoK9dJw9M4zdDd18uLOet7fWolaB3QFHxPrxwKlJitHvrG7n2vfzUEsSb1ySRWaEk2Ib1/d9+h4LsQHuRPu7sbGsjZOe3cIpacGcFKPmlBiJrQ2ufFIpcf4xKXj7t/BVcScr1+ahnRdMh8lGjL8bOrWKhGBnVH7wip/Ut3IWNXZT32GkxmAkIdiDC1/fMcDotWqJzGleyt79hz0tPL22ghkhHrx8YSZCwMYyPT8WtfJFbiP2Pg3CdQW19EZ5cNaCGH4tN3D9+3mE+bjy2sWZirqvxe7g9k92892eFq4/OpYVi6Ox2gU7qtv5eFcD3+x2kp2MVgc3f7SbI6f7sjDWl3WlBj7ObuKieeHcftx05Tmr0Pdy5du5WB2ClZdkkhji6aSY6408/F0pmyraWZbgPS6jB2ftx3hLgsdYlrsauBTYBJwDrB3r/h4OAVcfnESaHTt2kJSUpNAmJzP+aCq6AJvK27jxg1y6TDYcAlbtqKW528yShADe317Hh7ndLIrx5Ilzs/Dsy4HLAbkOo41nz0vnmBkBVFVVsaemhY8qNWyqaEfCWfK6/M2dLEkIICPcm7e21FBjMA2I3MNeFz4p1JMX+/rmgbO+/8/v5FCh78Vdp2ZThYElT/7CzEgfpvm48s3uZqb5uPLyhZnK6me22bnr80K+KmjmD5mh/P2UGbho1NR3mHjj12re317H6jxnVFutkrA5BO/saOIvS2I5MdUp2jmvzXmvQl3M5Obmkt3qfDwGr/gJQc7X17yXq2yB1ha1Dmio6a5z8uxvODYOgPd31PHAmmJmR/nw/PkZSiHNMTMCeX97PQBXLoykvrGZPL2DV3e08+qOrQD4u2u54ZgY/PuyIWabnZs/LGB9iZ7bl8Vx2RFOoo9OI6HTqPipRM80H1f+cXoSxc2yIGkj7+5wBhqj/FyI8tFSb+glzNeNyjYTV7ydiwPBa8szlAkP4KuCZjZVtLM01p2/HR0+LqOHiZF3xlKWC7wGvCVJUinQhnNyGDOkSVQdTeiDVqt1AMmmtraWPXv2sGjRoknv500mE7/88guxsbFMnz592OPe21bLP74qIibAnTtOSGB3QxfrilvZ1W9vmhzszlkJGqJczfj7+lDS48qDP9bi46bl+fPTSQz2YM+ePVgc8N8iwbe7WzhnZhg3HzudzZUG1hU7twSyQm36NC/OzArj6BkBBHm68Pj3pazcUruPC1/V1svV7+bS1GnmkTOTOS4piPz6LtYXt/LxrgZaup0sxwhfV5YmBXHMjACmB7o7S15rOrhl6XSuWLg3A1La3MV17+6iutNOsJeO9l4rlr5eeBLw0oVpeLroWP7GTuXaHz8rmcXRHjz0dTFfl3Ty8jIPgvtUb8oMNq56N5cusx0J50MwI9iDQE8tv5a3D7jPwV46Tk8PpcNk5cOdDSxJCODJc1Jx1TqvtbXbwp/fyaa8tZdH/5BEoKmO0NBQwsPDeXNTNY/9UIanTo3J6sDqEHjo1MyP8aOqrZey1l7uPXnGAKLP5goD167KJcTbhdcvzlK8A5vDwd8+3cPXBc1OQlK3haYu532M9nelqcuCTq3itQvTmBHiqXihz/5UySu/1HBmZggXJaoI7KPqjgeVlZXcddddrF69elyf68N+CwwcNFe/f3rN3d190kbf3t5OQUEBXl5ew9Y+2+wOHvqmmHe21nL0jED+fXYanq4aFicEclxyMH9+exdNnWYSgtwp1xv55y8OPHUqwrw7KWltJtpL4vYFrnhYDOzYUYjWy59HN7eTW9fJrcfF8acjnMIZJ6WGoJIkvt/Touxfd1R38MBXxTzwlVPBtsts5+TUYB7+QzKaPlXa7JoOrn0/D4DXL94rrJk2zYuv8pto6bY4y0/j/NlY1sa722pZubkGeRG6ZH4E588JV4z+q7wG7lpdiKtWzTPnprE0MZBeq50vcpt4Zl057UYbV7+bj1qCEC8dGpVEXYeZ/Ppu1hfr+bXaiK+7jqiEFISxk6e+KeDd3c48v7+7Bo1KornbytWLornj8z1IwH+WZ7CzuoPCpm66TTZe+7VaYSS661T8UNjKonh/esx2rnw7m+YuM8/+MQW3jirCIyMJDQ1l1fY6Hvu+jHkxvjx3fjoSElsqDHxf2MKa/CasfROXPBEuSQigvdfCjR8WEOnnxmsXZyoBUKvdwR2fOY1eVicSQlDa0stHO+t4b3s9NofAZHVw5bv5HBHrw5GxvhQ29/LO9gbOznJmCUqKiycUczoU6bpwEFZ8m82m6OEFBQURExPD5s2bmT9//oSrn/p3zS0rKyM6OnofmawOo5WbPsjj1/I2/rQwituW7U35bCzVc9MHeeg0Ek//MZWZkb6YbA5+LtHz1NpyKtucFFa1BDPDPYlz7WGap5q3iux0WeCBk+M4NWtvpd7rv1bz7x/LyYrw5rnz0vH30CGEYEd1B7d/upvGzr2deEO9XTh6RgA+rlre3FxDiLcLL1+YQXQfRdZktXPHZ3v4bk8LF8+P4PZl8cp5b600cN37TkENF7WKbosdjUpiXrQPahX8XNZOaogbz10wk5A+xl2vxc5j35fy/o56NH0uv5tWhXGEZhg6tUSUvxulLb1oVBIPnBjF1vJWPivsIdBNossKZpvg6kVR3Njn2tsdgvvXFPHRrgYWxfkT4KFlY1kb+h4rEqBRS6gkiXtOiGOatZa4uDiCgoJ4/ddqHv+hbB/voL3XylXv5lDY2M1Nx07HanewoURPdm2n8iD6uGm47bg4TkgJxtNFg8Xu4K8f7+b7whZuPS6OKxbure3Y3dDFFW9n46ZV8+x5aVS3mdhQ0sqGUj2GXiev5KzMYP7v+OmoJImSkhIiIiLGLPUuY+vWrbz33nu8+uqr4/pcH347K357ezu5ubkkJiYqogj9lXbHg6HkrodKD1a09nDFW7uo7zBxwdxwrj16umI8b2+p4aGvi4kLcue5c9MI93Oyssw2B+9tr6OyzchVi6JZHO/PN7l1rC1qYXuPAGyoVRInJHhj62xm06YafPz8+W+BmS92t3FSajAPnZGkqOWUt/Zy5+d7aO+18ux5acyM8OnLELTw0c4GbA6BSoJYf3d2Vnfg5aJBkuC6VXlk13Zyx/HxXLJgb253fXErt35cQICHjpcvzCTS35Xsmk5+2NPCx9n19FichmyXVHy0q55jZgTi46bhuvfzKW7q5vIjIvnzoiiOe3ozvRY7fm4aDEbnNV08P4LpAW7c+2UxYT4uNHSYKW1xrvQ6tcT6il5UKneghzaTQJIgxF1innsrZWXg4xfAP36s5bs9LaxYHM31fdLfDiH4IreR+78qxmoXWO0O7l5TQpiXjqPb2ugwNvNVQTMnpgTz8JnJ6Po8odZuC1e8nU2V3siz56axpK9HwdWLY/hoZz33rSnCy1WDwyG454si7l9TzMxIb9p6rJS19u5z7wrqnUbv4aLmjYtnEuXvRmqYNyemBPHvH8p4fVMNf8gM5b5TZ4AQmM1murq6UKvVA5SHxvK8HooluXAQDN9oNO4jqikb63h6j1utVnJzc/eRux5cqPNLmXM1dwjnHvG9bXV8tLOeeTF+WGwOtlW1c3RCAI+emawE8Sr1vVzzXq6zAWVfHr6+vp5lwT1EB8Tw8A8VBHjoCPFy4ZviDr4qgkBPLZJooaXHxkkxai6d4aBd30pAQAA7arq48cN8tGoVKy+dSXq4c9U4PSOUgoYubA7BrEgfYgPd+blUz0+lzqyMi0aFzSH42wnxXDxvr8rShzvruX9NEcmhXrx4QYbCrEsO9eS5daX0WByckxlMVIAH64r1vPBTJc//VIlKArUkce2SGP60MApXrZr3r5jFGS9tw2C0oVE5d+3nzAzD1udO2/upbHq6qDkqIYAdVe3KHtkhAAGXL4pj1qwQ6ptauf6DfHKarVya4ckFaV44HA7UajU5tZ089E0pvm5anjt7BpWlxbTqQthc3c0HO+qxC6dXZbbZ+TynkSUJATiE4E9vZdPUaebFC9I5YvrePfbnOY3ct6aIrAgfXrowA1etiuyaTtYWtvDBrgYlvvLOtjpq200sSQjAVaPiL6vy8HLVsPLSLML7sjVCCB7/oYw3NtUMUDuyWq0UFBQQHx+Ph4cHDodDWVjsdrsyAQw3CXR3dx929cHp6g9ekXNycoiLixvzzCjLXU+fPp3Q0NABfyspKcHHx4egoCDe2VrLQ98UExfozosXZhHm48qumnbW5DXxya56JSKdGubJcUlBHJsYhL7HzE0fFqBWSTx7bhozI30oKyujq7ubrxvdeXtrHUsSAnj8rBQ8XDS09Vj4NLuBl36uoqfvQXPXqZkf5UVWkIoWQxfvFVkJ99HxwnmpTA/xBZwu918/2c264lYuWxDJbcviUPVtFT7PaeTBr4sx2xxK4UuUvxtHJwTQYbTyeW4Ti+L8efKPqYrUdUu3mSv/u5MyvYm7T4jj/Hl73VqZXahRSQicXW1dNCrmxfjSbbKxq7YTlQQuahCoSAr1JC7QnY/76MGuGhV3nZTAQ9+UEurtgqeLmrz6LvzcNbhoVDR1WhCAn7sGCQlDr5X/OzGekxM8aW1tpa2tjaIOiSe39xLi5cKzZyXQWu1sJ+7m7qFsCY6Z4aQS/1Sqp6GvxbdWLSEBd500g7NnhqHqm+A/3FnPfV8WMT/Wj+fOS1eCo70WO9euymVrZTu3LI3Dw0XNTyV6NlcYlBZfrloV1xwVw+npoYR4u/Tl9stYubmGC+eGc9eJe41+165dTJ8+XfFOZcgTgBBiwEIz2BtYtWoVra2t3HnnnWN6tgdhv7n6B9zw7Xb7Ptz8/Px8IiMjR21tBXtbYQ0nd11WVoaLqzuv7Gxn1fY6jk0M5LGz98pTVep7WfFONjUGI9csjkalklhb1EpefZcyho+bhntOmsHSxACK9uxGaFx4KcfMTyV6LpkfwV/77bNzaju47v08rHbB42elYHcI1hY7XXiZu+7louasFG/SfawEugrUHn78a2Mbhc29/N+JCVw4d+9qvraolds+drbSfvnCTFy1zvTUD4UtbKowIITTGI5PDmJpUhCL4vxp7bbwp5U7MBhtPHF2Ksck7WUofrCjnge+KiI1zIvnz8/A21XD9qp2vi9s5vOcJmXyiw1wpcfioLlroDZCfJA7j52VSmKIJ7+U6rn9sz0Yeq34u2v54cYjcNWqaTda+Tq/mafWldNlcv62WrXEvGhflswIRKOS+Ne3JUR4a7kxU43a2sO0adMICg7l4fX1fFXQzNV9FGU5TrKuqJW/fbYHk82Bo08HIcBDy+L4AFSSky14VLw/T5+bpmynesw2rnkvl501HTx0xl7GJDjTt39ZlYdWLeGuUyvXmRjigaeLhh3VHWM2+sFwOBwIIRRpchmSJPHmm2+i0Wi4/vrrRxxjGPx29vhDYayluXKr7JHkrnuscMd3FeQ29vLnRdHcsjQeVZ+Rbq5o44ZVuUgSvLo8g7kxfkiSxJVHRvOPr4t5f4dT/77bbOO2T3bjqXUG88rbe2jsMClNL2R8u7uZOz7bQ7CXjhcvyFCEK4+Y7o/JaueLvCZmBHvgEIKVu5xMukhfV9p667HYHFyToWGWVzd6vR4/Pz8+2NnAP74uJiXMixfO3+vCn5YewtqiVoSAE1KCcFGrnP3u852imuA0tMfP2mv0Qgie/6mSFzZUclS8P0+ck6asikmhnjy9rhyL3cF1S2Jw1arZUKKnqh8NV6OSOCsrlFuPi8PLVYvdIfiusAVDr5VT0oK59+REJfDWY7axcotTZ+DlC9Nx1ahZX6JnfbGeh74pAZzbltmRXjT2dHHGkfPo7O7hlk92s73BwoWpHpyf6oHdbkej0VDc3MPf1xThqlXzzuWzCPZy4ecyPT8V6/m6oBmzzYEEmG0O3t/u1DX0d9dx9bs55NV18dhZKZyUujezs63SwHXv5zPNx5U3LskiyFNHaUsP64tbeWdrHUVNPVw8P4I7jo8ft9EDA2TkYK834HA4+Pzzz5kzZ86oYxxoHBIrfklJCb6+vgQFBQ35Gbk0FyAlJWXY/VRZSw9XrtxOfZeVSD9XzpkVznFJQcQFefDBjjoeWFNEtL8bz52XRpS/e78WVntYV9yqtLBqaevko40F7GrXsaW6G4GzDfSieH+OTQxkSUIAn+U08sSgyD04swc3fJDPtqp2hVEmSRI1BiMrN9Xw/s56Zd8c4KHliGgv0v0c7Kjp4rtqO/MjPXjq3Ax8PJw56JZuM9e8l0tRYw/3nTpD4fDbHYLXf63i6XUVSEDflpzEEE+WJPhT2tLL2qJWzswK5b5TEpVGlnXtRv78Tg4NHWb+fXYqxybufbA/2FHHA18V4+miwSEE3WZnlmBOtA/tvTYKm7oHrMzKPX87B6PVzssXZigsQnCqEP3zmxKi/d3wdZHIb+zFLsDbVY1O41QMuu24OM5O9aGlpYW2tjaquuCxbUY8XDS8fvFMhTwkhOCFDc5YxfwYX5JDPPm5vI0yJeiowuZwcM1RMVy1KFq5Xjm3H6YYvYsy3j+/KeHdbXVcMj+Cv03Q6IeDEIIHHniA5uZmXn311YnSz387rr7D4cBqtQ54r7y8HDc3tyHbU5nNZrKzs0ctzd1Q0srNH+ahkeDIaHdquiVy65x14F6uGrpMNjLCvXjx/HT8PJw/fmOns4VVcVM3d56QwPJ5EbS0tFBWVobeNYJ7vy7Dz13LisUxFDV3s7aoVdl7gtPIHvlDsiJfVWMwsuLdXGrbjfzj9CROS9/ran6Z18Rdq/cQ6efG42elUNLSw7oiPRtKWuntS6WFees4J9mDBA8zfm4aTDof/r62GUOvlSfOSWNJwl5W4/vba3nw6xLi/F147dI5dJttrCtu5cc9LeysdV63m1bFKWnBHJsYxIJYP6rajFz1Tg5mm4MXzk9nVl8dAcBbW2r417elzIv25dnz0nHTOQNl3+1p5qNdDZj6zjEh2IOjEwKcE4YE17ybh0Yt8eryTOU+CCF4eWMVz6yr4NjEQP660J/mxnriktJYX9rOv38sU6THNCpnx5yjZzjbgd2/phgvncQd893wVjvbZQcGBvLGzjZe+9UZbX/wtCRlq5Vf38n17+fT3GVG1Udx9nLRcGScP9N8XXh7ax3R/m68fnGWooMghOAfX5fw3vY6LlsQyV+XxU250T/88MNUVFSwcuXKyWhI/LYNv6qqCrVavU+Pu/5y18P9EEIIVm6u4ZFvi5kR4skDx03DU2UlLi6OspZubng/j9KWHoVlFuCh45gZAcQFuvPaphpMVjv/PjuVRXH+1NTU0NzczC5jAE+tqyQ93Lmayy53p9GZS86t61K06ACmB7qTFubF+hI9kgTPnpuu8NOFEPxno3NlnhvtyzPnpuHTpyLTabJy3ao8tld3kBnuTWOniaYuCxLOSrMagwmtCm6ZrWPO9GCCgoLw8vLiuXXlvPRLDfMiPXh++SwlwNdudI63q6aD0zJCMFkd/FLW5hTAVDs74rjr1Dx7bjpzY/ae3zPrKnh5YxXHJQXy2Fkpyp65vdfK1e/msLuhmxuOiUWrdspr76juwN733Ljr1Nx63HROzwjFQ6dBCKGkxE7PCOGaWV7oW1vIyMig0+zgz+/kUNLcw6NnphDspWN9sZ71JXqlV59WLfGHzFBOTQshfZonne0GHv2hgq9Ke1kW68rtS6MJDgpCq3W27/7TW9lUtxl55tw0ZkX5sKncqXb8Q2ELnX3xhvQwT45NCuLoGYHEB7nzz29KWLW9nsuPcPYYnGqjf/LJJ8nNzeXdd98dtqPyGPHbNvza2lrsdjvR0dHKe7LcdWZm5rDpEIvNwQNrCvlwZz3HJQXx6Fmp9HS00dHRgWtgBCveyaZS38v/nRDPyWkh/FzWxtqiVtYWtSr7xCOm+3F6eghhtKFTCT6qUPNxdiMnpjjz8PI+tq7dqZZTpTdy36mJnJUVRl27iXXFrXy4s56SZueD6+um5bikQI5NDGR2lA+PfV/GR7ucXXD+eXqS0giyrt3EivdyqNIbeeiMZE5ND0EIwZ7Gbl75pYrvdrcoNzjG3435ke4kedv4sayTjfWCpdM9eOzcmbjqnJNIQ4eJq97NobrNyCNnpnBiSrByj/6zsYqXfq5EkiTsDoEEpId7syTBn+KmHr7d46Qb//2UvQUoTZ1mrnwnm5o2E/8+J4WliXu3YV/mNfJ/nxfiqnVeS7fZrgTzjFanlNUFc8JZnqyls6OD9PR0WntsXPF2NnXtJp7+Y5qi+Q9OTsKNH+bj564l2s+NXbWd2BwCLxc1vu5aagwmzp0dxi2Lp9Ha2kprayudFnhsh4XmHjsvnD8wzfdTcSs3fJhPhK8rx8wIZGtVO/l9wVt3nYpei4MrFkZxy9LpU270zz//PJs2beL9998fV4v3YfDbCe4Np8JjNjtdaCEEpaWlSvuq4XL7bT0Wrnx7FwX1XaxYHMONx8ahUkmY1GpyG3p57JOt2B2Cly9IZ8F0fyRJ4uTUYGoMRr4uaHau0tO8+LWsjTvKDahwbgk6TDaWzw3nzhMTlNRRbl0n167Kw2p38MryTEW+apqPCxabg9LmHtLCvDhrZhhbK9v5usDpHsvNJo+ZEcD/nRCvGP3uhi5ngYvVwasXZTKvT31HkiR21XTw3e4WMsK9ufeUGeys7mBdcSsf5eqVffx0Px1zQrVs3roNPy8PejTe/N939fRa7LyyfO94AJ/nNvLSz5WkTfPmxfPTaeqysL64lR+KWnh2fSXgzDq4aFRsrTQwJ9qXunYTV76dQ4fRyit9QVAZq3MbuevzQlLCPHn5wkw8XNTsrO5gbWELH+c0KrnzTSVNmLu0/HFhEnUdZq58O4e2Xiv/uXDgeF8XNPO3T3eTHOrJy8sznc0zzTY2lrXx5I9lVBuc1X4f7WigvKWXJTMCSQtL4PEvC2nttfG3+Z7QVEShxY+goCB2Ndu55ePdJIZ48upFmYp31dRp4q+f7GZ7dQdXHhnFzcdOvdG/8sorbNiwgY8//ngqjH6/4oCv+EKIfeS0W1paMBgMTJ8+XZG7njFjxrD7+ZLmbla8m0N9uxGHgJgAd5YlB7EsOZjcymYe/rGKCF83nj8/jZgAZycVi83BPV84RR5OSw/hgdMScVgtZOfkUCP8eXJjk9K1BiA1zIuliYG46VQ8tbaCIE8dL124N3Jv69PFW7W9nhNSgvjXGcmKd1BjMHLFW87VzbOPl6+WnHvZmAA3vshrwsdNy8sX7q0TdwjBkz+W89qv1RwzI5DHz3b2tgeny335yu0UtZhICfGgtsNMp8mGTi2RGOxOcXMvLmr461xXZseFEhQUhLu7Oy//XMUz6ytYHO/Pk/2i+t1mG9e/n8eWynZOSAnCaLGzuaIdi92Bu1aNzeFAq1bx/PnpAyaRt7c6qwnnx/jy3HnpSnMKk9XOLR85q+UuWxCBo7eDbfUmivRWHH2ltmqVxA3HxHLBnAjlPD7LaeDu1YXMjPThxQsylJRrf3799UfHsHC6v7IlkBt2SsDxKUH8cdY0Zkd6093ZwRe7qvn3pnZifdU8fmosMeEh6HQ6HEJw35dOrsBVi6K5sS84OZVG/+abb/LFF1/w2WefTWVn3t+Oqz+U4be1tVFXV0d3dzfR0dFMmzZtmE873cJbPsrDXavmwdOSqe8w8UNhC1sq2pQV0c8Fbl4SyWmzY3DROkk2N3yQz86aDiXS3tHRwZ49e7D6RnHnmgocQvDMH9MI8NSxtqiV7/e0KLl9nVrFWVmhnJwWwsxIH0xWO7d+XMCG0jauWBjFzUunK95B/zZWT52TysI4fwrqu/ixqJXPshto7quumx7gxvEpwSxNDCI+yJ27vyhiTX4T58+Zxl0nzlBc7lqDkT+t3EFTt5VH/pDEiWlhWO0OdlR3sHJzDT+V7NVeSAvzZO40FxI9LXxX3svaGjsnJPrxyNnp6Pr27foeC1e/m0NxUw//OD1JyXX3Wuys3FzDixsqcQiBQzhZfnOifTg6IYD6DjNvba1laaJzUuqfO792VR7bqtq5+6QEMtw70Ol0xMfHs62qvc9TEqhVztp4F42KBbF+eOjUfFXQzBGxfjzbj4BjsTm49eMCfixq3YdfX91m5JKVO+k02kgK9WR3QzcWuwMPnZq4IHfy67tICfXimbPiMXU5W2LZ7HbeLZX4oayHqxZFceMxU7vSA7z11lt88MEHfPHFF5Nq8z4EfjuGDyhuvYyamhqKi4uZPXv2sHLXcvHLY987G1O+cEEmoX2NInvMNm78II+fS/VE+LrS0mXGbBe4aSApQEtFpwOjVfDQGcmclBpMQ0MDNTU11KrDuP+bMkUvL6ZPcKK/4GVqmJeT9FLdjtUu8HXVICToMtm444R4Lpq3lwP+a3kbN32Yj7tOzYsXZJAc6qWc+3M/VfLihkpmRfpwVEIAG0v17KzpwCGc/HeLXXBmZij39tXRA+yu7+TPb+/CYhe8cEEGc2P27mNlYk7aNC9uXRrH9uoO1hcPJCLNCHTl/GQXInQm/Hx9sLn4cMc3tTR1mnnyjwOzBGuLWrnlowIi/V15+cIMmjstrCt2xkPKWp0pMy8XNX+cPY3jEoNID/emy2RTgn//PD2RGKkFLy8vYmNjya7p4Op3c52tri52tgHbXtXOT32iobJ3lRTqwbEzgjhmRiDTA9246aMCfi5t4/9OTOCifjTlSn0vl/83G5PNzmsXZZES5kWvxc7mCgNvbalhS2W7cmxWhDdHzwhkcZw/KzdXszqvmXOS3DgxwoGvry/+/v5UVVURFxc3aaN///33WblyJWvWrNkf1NzfruHX1NRQU1OjCGcMBYvNwd+/2MMn2Q2ckBLMI2em4ta3QtS3m1jxbjYlzd3ccXw8y+dFYLY52FRu4J2ttfzaJ0GlliAjSMOcMC2ZgSryzQG8tLFa0cuTG0h2mWzc/FE+v5Yb+PORUdx4rHM17zbbeH97Hc//VKmw3dy0KhbFBbA0KZBOk41HvytleqA7L16QQVjfpGS1O/j7l0V8ltO4T069qKmLa97Lo6nTjEbt7FDroXPy4SN8XHlrSw2eLhKvXTJ7gHSUnM8eTMzpMdu4+l0ncy02wI3adhNWu8DHTUNqiDs5dV04hOD2+Z4cmRhGUFAQLi4ufJ7TyN2r9+7b5Y47NoeDe79wnvucKB9UEuyscQbe/Nw02IXTU/jXGYlMszUSEBBAVFRUX+48jyAvHa9dlEW4r6ty7i9uqOS5nyqdra8ifdhQoienrhNHHyPRahdcODec246LU7ZOpS09/OmtbBwOwesXZw3o/rM6t5H/+3wPs6N8uf7oWDZXGPipRE9Bw94J8NolMVy7JBaHw0Frayt79uxBrVbj4eFBUFAQgYGBE3LPP/nkE15++WXWrFkz7qq9MeK3ZfgWiwW73a7U48fHx1NYWDis4d/6UR5f5jVx8fwI/u/ERIWJl13TwV/ey8Zsc/Dvs1JYFB+gxAVkdZv4IA+uXhTNzpp2vslvQG/ce9ozAl341+mJJEc4V766diPXvJdHpb5XidzL+LlUz80fFeDlouGZ89LoNNr4saiFHwtbaOl2Zim8XTVcsTCKk9NCCPd1pdts46YPnZPItUti+MtRMcr5lTR3c/W7TgWgp/+4NxW1triVb3Y302N2BslmR/lwUmowx8wIJNjLhQe+KubDnfX8ITOU+0/dO4noeyyseDeXwsZuHjw9kT9khtFjtvFruYGPdtbzc5mzMYZGBXOjfJgVoiXB3ciWBivvFVqYF+XNcxdk4OniNHqzzc5tH+/mx6JWrlsSwzV9595htLI6t5Gn15Yr/AOtCmaGe3BSRjg6tYr71xQTHeDGq8szFVWh/mm+wbn4WkMvV72TS2WbUfF+XDUqjpjuR1KoJ+9tq0OjVvH6xXsbkQB8mu2MEcyL8eWFCzKUmIjdIbj14wK+29OiGL383GVnZyvufW9vL62trbS0tGCz2QgICCAoKAhvb+9Rpd++/PJLnnrqKdasWTNkV6cpwm/L8Lu7u8nOziYgIIDY2Fjsdqc67bx584Y8/sWfynmqryVUUqgny5KD0alVPLuunBBvHc+dl0Z8kCeSJGFz9G9hFcCjZ6aglRzk5uai8w7g7+taKGrqGZCHj/KSSApxZ3OtCSHgqT+mDUgPvb+jjn98VUJ8sAcvXZCh1LZb7A7u/aKQ1blOaq7V7qBC76zdjw9yx9Brpd1o4/5TZnBWP9VcmULqqlXx0oUDtwSvbqzkyXWVxAe4sCAuiI1leir7xpQFPM7OCuP+U2coDMZag5ONN5QL/1OJnps/zCfE24Ubjoklr66LtcWtVPdpDABE++q4MsOFMJ2FgIAAPHz8ufe7WjZXtu/jcvdn6j19dhJ7ikspM7qzpbaX2nZnBN5Nq2L5vAhOSg0mKcQTATzUx5K7YE44d520N2PSYbRy1Ts57Gns5tGzUjh2RiDbqtpZV9zKd3taaO2LicwI9uCElGCOnhFAUognH+1q4O9fFrFwurNIx7Wf0d/5+R6+zGvihqNjWXFUjPO3GmT0g2Gz2dDr9c5UYWcn3t7eBAYGEhAQsE8u/ttvv+WRRx5hzZo1k5aKGwW/HcMXQvDrr78SFRWlyF0LIdi8eTNHHHHEsJ+raevl+z0tfF/YzK6aDoRw8r/PzAzh9IwwMiK86TE7g24by9q4bEEktx4Xh7G3h/z8fLQBkdz5dTUt3RYe/kMyJ6QEU9XWy4+FrXy0s14R2/BzgcXRHpyYHsbCxGk8u76S136tZnG8P0+cnapEsjtNVm78IJ8tlQOpuZX6XlZtr+O97XWKUkyEryvHJgayNDGQxi4zd68uJNLPjZcvzFTcYIcQ/HNNIe/tbOTYeB+eODdLSf/l1nZw68e7qeunZR/eN2ZcoAfPra/AYnfw4gUZimoP7E29JYY6xS1l5prN7uD2T/fwze5m/D20tPVNgNN8XJg9zZ0dtZ00dtu5ZpYnZ8+OJCAgAK1WS0F9F39+JweNWuL5PybRW1+q7JM/3lXPvV8UMc3HFX8PLfn1XQicyj4eLhrKW3sHUGOhLyX7Tg5lLT08eU7aAPpwbl0nV72Tg4tGxclpwWTXdJJb5xTd8HbV0GmykRrmxWsXZyp9AW0OB3d+Vsia/CZuOnY6Vy1y8kJGM/qhntGOjg5aW1vR6/VotVoCAwMxmUzU1tZy//3389VXXw1LMZ9C/HYMH/YN7gH8+uuvLFy4cPQvFYKmTiPf727mx+JWtlV2YHMIAj20OIRzBbnrJKcWm9zVtNcrkju/LEencaaoMvrq4YUQvLGphn//UEZKmBcnpwWzudzA5goDVodAowKbAxZEefL0eZl4uTkNp67dxDXv5VKp7x0QGQcnN/yGD/Jw06p5+A/J1LabWFvUyi9levpCAwR4aLl9WTzLkoNw1aqdbvWHefxYYuD8mcHcfWqKsiI2dpq46p1cqtp6efgPKcyN9mV9cSs/FrXyS1kbtj5CzlEJAZyRGcriOH88XDT8d3MND39XyvwYJwVXTpVZ7A7u6DP6K/syEvoeKz+VtPJ1QTObyg0InKW4i6f7kBmkItbVSEMvPLHdhI+bhhfPTaGtuojExET8/Px4a4uzM/CiOGe1nJvWycNfW9TCCxsqlUo4d62KRfFOum9qmBc3fZRPrcHEM+emDdCp31ndztXv5uLvoeX1i2cqk2Nrt4VHvy/ly7wmhSPhplWzcLofi+MD2Fim54fCVm4+djp/nqDRDwWTyURVVRVXXXUVpaWlLF++nEsvvZS5c+dOaLxx4Ldl+IMFN2Fshi/XPjscDqVJQafJ6ixb3dPKTyWtWOwCb1cNc8PdSPG24h4YxuM/VDI9yJ0Xzs9QHiJn5L6ED3fW78PSq+nbc1a1GdGqJKwOgYsaZoboSAn35vPdHVhsgqfPTWNB7EByy92rC4kJcOelCzOY5rN3Nf/HV8Ws2lFPqLcL3SYb3RY7rn018ZWt3VS3W7hxSRRXL4lTxitt6eGqd3LoMjlVfft/1w+FLdz2cQF+Hjoyw73ZUmmgw2hDo4IwH1dqDCaWxPvz9LnpiufQY7Fx0wfOTj1DpcqufDubtl4rVyyMpK7dzE8lreh7rIqmn7eLmuuytISoewkLCyM8PJz3cgw8s76C45ICefysVOW7zDY7t3zk1Bu48ZhYkkI8WVusZ31xqyIYqpLgnJnTuOyISCWjsq3SwIr38gj2cjYXlQUzAWUyOzYxkIf/kMTO6k7Wl7SyrqhVEQa5del0rjhy6oxexq+//spf//pXPvjgAwoLC+ns7GT58uWTGnMMOGz4Qxn9YJisdn4pa+OTrWVsqzPRbXWuhv4eWq4/OpYTUoLxcdMOiNz3b1UFUNHay4r3cmjusvDImckcnRDIpgoDPxa28M3uJrrNzvOe4afmpGR/TsmKIjzAi1d+qeapteXMj3FmCWT302yz87dP92rm/e34eOwOwbaqdr7Ma1L05CVgTrQvSxMDWZoUSFOnmb+sykOnUfFyvzgA7E3lpU/z5sULMvB112JzONhe2c6/vitV6MMAGeHeHJsYyLxoXx75rpS8+k7uP3Vgp56ipm7+/E5OX1lthqIQZHcIXvmliufWV6BRSYo6b3ygKwsiPahu7WZDjZmjoly57+R4ggMDUKlUGK12bng/j1/KDfuUMtcYern4zV0Yeq2EertSY3BusWIC3EgM8WRdUSsRfm68cXGWEhgEFC2+ZUlBPHZ2iiLLZe3T1/92Tws3HRvLVYtigKk1+q1bt3LTTTexevVqoqKiRv/A1OG3ZfhDqfCMZPiyyIHcIWe4iKvNZiMvLw9vb28io2PYUd3BmvwmNpa20dxtUSrBylt7aeu1cv+gyP32qnauez8PjUri+fPTB5SYvr+jjge/KibK3435MX78Uqantt25ZfHWQacFjorz5ZnzMpVVTy6a2VnTwd+Oj+fSfrpvJc3dXPnfnXRZ7Ny6NJ7WHitri1sVo5UAbzcN/zoj2dl7r0+g4qWfq3h2fQVHxfvz5B/TlEi22eZU9PmhsJWrF0VzUmpQnyBIq8JTByd9+E8Lo8iK8EGtksiu6WDFe7m4alW8etHAqPnHuxr4+5eFZEX48K+To9iRX0yzJoQN5e3srNlbAXj8DF8yAiBc24PWxZ2ndpooaDLy4GlJnNnv/lbqe/nTW9kKtTg93Ju6dhPri1v5JLuBPY1OVp63q4ZjZgRy9IwAFsX58862Op5aW86JKcE8cmayksmw2h38ta+pxl+XxXF5n77+VBr9zp07ufbaa/nss8+IjY2d1FgTwG/f8Ddt2jSk0q4QQqnfH0nc0GQykZubS1RU1D5yXA4hyK/v4oc9LfxQ1KJEybMivFmaFMRxSYHk1HZy9+pCovzdePGCDCL7GkX2p9IuSQjg8bNTlCq03LpO7vhsD1X9IuThnhJHRnsyLy6YF35toMZg5OE/JA8QhtheZeCad3PQqiVeuSiL1Gl7J5iXf67kmXUVuOnUCu99mo8rxyQG0Nhh5seiVs7ICOWB0/am8rrNNq5blcfWqnbuPCGei+fvnWCq2pzEl7YeCzOCPSls6sbmEPi7a0kO82RbZTvBXk4d+v7NMWW3+sjpfvz9uHDqqpwFUxqdjntWF/F5biOL4/3RqSV+LW/vazmtQquS6LbYuSzVhaOjdAT16fHX9zi18+wOwWsXZZEUujcXv7aolZs/ymd6oDsXz4tgS2U7G0r1dBhtyl4+LcyLJ85JVc7RandwW5+Cbv9JdSqNPjc3l6uuuoqPP/6YhISESY01Qfz2DX/r1q3MnDlzQFGOrGIyWtPBjo4OpbPucMw/5aT7NNV/LHROArsbnKuMSnJKXf/rjBTmRPv0Ke3a+b/PC/m6oJnzZk/jrpMS0PRNPi3dZv7yXh57Gru4+6QZLI4PYG1RC9/vaWFHdQcC5682P9yFC+eGc1RKBDqNmm8KGvnbp3sI9tTwxqWzifD7//bOPK7KMu//78Mui+zngLih7AqoaDmOmnuZKIummE2WNrbZL6enbJke0xonbSYbZ5qmZvLpMdcEydzSUlMzU59yAVSQHUE4K0eWA5zt/v1xODdgliAIBuf9evHHMbjPdTr357qv63t9v59vk9FE8+SW9+YMpc5gspSYXlZwvLFLrJODhAcipUyJ9GfMIB90ehNPbrWUuq5OaOkBcLmimiVbLmAyW/zuh/bpTU2DkW/zNGw5c1V8ajvbSxgz2BJ0mxDqw/Yfr/HPY0VMi/TnxbF+XLtawrBhw8DOgZfSLWJrfpLRYDTx1SUlbx/MRduYkWcvkTCsrwdxgU54CDo+Oq/D0d6O92eHMWxQgPh9fnVZwYs7LUU6/17QVFRjMJlYsSeHLzLkuDnZi36GEQHujA/xJfNaFd8XVLaY6DpS9JcuXeLxxx9nx44dREZGtuta7aB7Cf9mLjxnz54lKioKFxeXVu3nrcjlcoqKioiJibmtphxl2joOZ6vYnVHB5QqL204/bxfGDvbl7FUtOfJa/mvyIBY1605ToKrlya0ZaGr1vDt7CBPCmm6yk/kantth2ZuHS904X1qF3iTg7ggB7g7kVRoJ93Pmk8dGiRlyJrMgWn/d7GluLahJiAnALAgcy1WLRTrWdlh/TohkxtCmVcWPJVqe2ZaJm7M9Hz8SKxYXAaSfL2fFnmyG9unN4jH9OF2k5UiOqoXf/5BAD14Y44tjnZrY2FiM2PH8jixO5Gt+Ylctr2pg8ebzXNPW8/d5lvjGkRwV3zTbujjYSZge7sVvAsDfXoeXpycZWif+fKSUmKDefPhwrNhWq7lHQPKwQFbFh1Os0fHNFTWHs5WcbzQaaZ5j0JGiz8nJ4dFHH2Xr1q1ER0e361rtpPsL3+q029zC2OpYetM3FwSKioqorKwkOjq6TdbcP4e6Vs+RHJXF2LJAg9FssZR+cKiMKeF+3BPsTWZZlRgH+CClKRAGFrvn/96TzSA/Vz56OBZZb2dq9UZO5Kl5/5tC8hu3GI4SGCZzZFK4L5OH9uWdw0UcylaJx2vWz6yq0YtP8+bHhgaTmV0XKnj7YC4NRjNCYwXc8H6eTA73w93ZgdUHcunj6cJ/HokVTxcAi2lJ4xJ+/dym4hijycyy1IscuaLC29VRdMgZ5GfpInu6WEtORQ1v3tDCu0xbx6JN51HXGvhwfoxoQAKWY7klWzJwdJAQ7NOLzGs1mAQBXzdH+nk5caGslhAvCa+M8aR/oBQ/Pz+cnJxY17i1emiExSPAGnjVG838IS2Lb66o79jyPj8/nwULFvDpp59aVjldS/cXflZWFn379sXd3f2WQTxrO207OzvCw8NvuwPPL1HTYOTYFRWHclQcz9VQZzDh5mRPvdGMj6sj/3kklrBm+fNWl517B3rx97nR4tPLYDLz2udZ7LukZmaUL/GxQXyTq2qR6gvw2wFuvD49nAFSy36/RGPJxlPVNPDenJbGFdanuauTPf9+OAaDWeBwjuWaVxqfsM4OdsyP60N8TIDYf/4fRwv58Nti7o/yZ21ilBiE1JvMvPz5JQ5eUvL0uAHMGAAFci3ldv58dUnBj41bAg9nex4YYrXy8qL8eoMYrPv3glgxPwLgdGElz2zPRNbbSexjd73OwPE8NZ+eKhVz6V0cJNw7wJPhMgcGu+jYV2DgQJGB2bH+rJwZhX3jd6s3mlmWmsXRXHWLk4KOFH1xcTEpKSls2LDhbjHI7F7Cv5kLz+XLl/H29sbX1/cXhazX68nMzMTf359+/frddjvtttBgNPF9gaV328FLSnR6E84OdowZ5MPEMF/OXb3O5xcqLC47CRHiUVNtg5Gl285zuqSaJWOCeH5yqDje8uv1PPbpOcq09fi6OqK02nh52hEpc+VEST1I4MP5LQ0sv8lR8cLOiwR6OvOfBU0FMICYSNPfuxd+7k6cL7VU/8k8nPDs5cgVRS1JsQG82SxPvs5gEpfwL00dzFi/BhoaGoiKikKjM/LEZouL0fxRQVRUNfBtnkb8/IIg4Ghvxwfzoxk1oCnH4ESemud2ZP2kjx3Atv8r460vr/DbQd48PKovJ/I1fHPlhm2G1IXHhzjh42jA29ub3t6+rDpUxvE8TYtGmR0p+tLSUubOncuHH37I6NGj23WtDqR7C99kMqHRaCgoKEAQBPz9/ZHJZD/Zs9fW1pKZmSn2WesKjGYzZ0uucyjbsiWw3rCBns4s+k0/pkRIkfV2trR9+vTHxgYXg0i5p8lWLF9Zy+9vSMzJV9ZyJEfF5+eviT3m/VzgvmB3Hojuw+iwQPZkyvnv3TlEBLjz0cMxorNv87Lf5r55lTo9h7JVfHCsCHm1ZZwezvZMCLPkCsQG9ebF9EucLbnOGzPCiHatQiKREB4eTnlVA4s3nUdRref9eU21C3qjmc9+LOPdQ/kYzS3r9ieF++HiYMefDuQy2N+NDY/E4u3a5ERjPSmYGObLe3OGiisOk9nMi+mWFYePqyOaxm3GQJ9ejApyJau8issqA0uGuzP/nv74+vpiNps7TPTl5eXMmTOH9evXM378+HZdq4PpnsK/WRBPr9ejUChQKBQYjUb8/f2RSqU0NDRw5coVhgwZgoeHx63fpBMQBIGL5dXsz5JzPE9DQWPdelSAO1c1OuqNZt6bM4RJzRpcnLt6nae3ZeBob0nMiQps+ixfXVbwUvol+ni6kBAj43SRlh9KtJjM0Mse6kwQ6ufCJ48Ox8e9KSvQWgCTPCyQlfFh4smD1TDkmytqnr1vIGFSN47kqDmaq+J6nVG8q+YMD2RaoB6ZlyuDBw+mWFPH4s3nqak38eHDLfP/z5ZoeWpbBr1dHPn4kViq6o0cyVFxOEcpWl27ONjx8KggHhwqIzLAUjy14bti3j1cwLRIf95JbkrAMQsCq/ZZKg4X/cZSX1Fe1cCxxtOMU0VaBGBVfBjTBrujVCpRqVTU1dUREBBAcHBwuxxv5HI5s2fP5q9//SuTJk267evcIbqf8PV6/S0j9waDAaVSSXFxMXV1dQQFBdGnTx/c3d07ZYnfVvKVtRzKVrLvQil5GstTa7C/K1Mj/JkS4U9FVT0v7ryEtLcz/1kQK+YKAGz/wZIgNKxvb/45P0b0B9Dq9Lz6RTbHctVN+ekOEBfgxIQwP85c0/NVtqqFYyxYYhTPbs/kh2It/31DH/kybR0LN1pagnv2ckCjs0wCw/t5Eh3kwe4MOQD/WRDbYmL6vkDD0s8yxXP/wGZBQ2tdfKCnC/7uTmQ01tgH9HZG6uFERlk1D0RZRG+dmMyCwBt7c9h5rpzf/7Y/yyY1BTbrDSae+yyTkwWVrHgwjHk3LO+DgoIwmUwolUpMJhO+vr5IpdI23RsqlYrk5GRWr17N/fff36bvupPoXsJXqVTY2dnh5OT0i0E8q/GmTqcjIiKCyspKFAoFOp1O/KJbUzvdWVgbfzg5OeEu7cfhHIvN848lWsuS2E6CVy8HVsWHc1+oH/Z2khYdbyaE+vLunCFiNp7RbGbl3iukny9nXlwfXpwymFNFlRzJtgTyqhpr9vu42zF/mB8PDh9AoLc7lTo9S7ZkkCOv4e3Elsd8JY1Pc63OwD/mDsGxsohqe08uXrdnT6ZcLNft792L6UOkTAr3Y2gfD47nqnk+9aKlzv6Rlvt2ax+7ewZ68c+UGFyd7NHUWhx8NpwsEROm3BtNRiZH+DFmkA9rDubxRUYFT48fyNL7mrwK6gwmlm7P5FRhJW/NihCzK39uT28wGMS6+traWry9Lcab3t7ePxsv0mg0JCcns2LFCuLj4zvqFuhoupfwd+7cycqVKxkyZAiJiYlMmTLlJ15lJpOJrKws3NzcGDx4cAtxm0wm1Go1CoWC6upqfHx8kEqleHl5ddkkYDQaycjIwNfXt4VNONAoAjV7Myo4e/U6hsbMuYlhvihq9HybpyEx1nJ+33yZ/mL6JY7kqFq0mgZLBeLT2zK4UFbFPQO8KFbrqGj04w/2lKBpkKAzCKx/qKmlNECuopYnNp+3lPDOG4JRUUBQUBCBgYH8UKzl6W0ZeLg4kBgbyNkSreif7+niQHWDkX7evfh04fAWOfRWA85xIT6sf2ioWOgkCIJ4LJcQG8DkML/GtloqNDqD2OfgvlBf3pgRJhbj1BlMPLs9k9OFlaxOiCAx9pdFfyNms5nKykrRwLW5y471yFer1TJ79myWL19OUlJSu7735tTX1zN+/HgaGhowGo3MmTOHVatWUVhYSEpKCmq1mri4ODZt2tRaF97uJXywfEFnzpwhLS2Nr7/+mpCQEJKSkpg2bRpKpZKCggIiIyN/0XjTeh2NRoNcLqeqqgovLy+kUukvzvYdTUNDAxcuXLhpuvCN1DYYOZ6n4XC2kmO5amr1JhztJEyJ8GNKpJTxoT6YzbD0M4uB5Y1GGMrqBn6/5QIFKh1/SY7i/iip6MefdvYaaefKMTa26OrvIWHsQA+mxwTh4OzCU9sycbK341/zhlBTdoWBAwcilUo5kafm/+3Ioo+XCxseGSYajWh1Bt47UkDq2WuiUD2cHbgvzJLlV6iq5R9Hi5gU7se62U2VeYIgsOarPDadLiVlZB9enx4mnsU3GE08vTWDU0VavHo5iJl+QwI9GBfiw6nCSks77YRIEmIt/y9vN3ovCAI1NTViXODq1atkZmZy/PhxXnjhBebOndvqa7X2/Wpra3F3d8dgMDB27FjWr1/PunXrSE5OJiUlhaeeeorY2Fiefvrp1lyy+wm/OdYIbWpqKunp6VRXV/Pkk0+yZMmSVnXQbX4drVaLQqGgsrKS3r17I5VKb3lE2B6sJw1hYWH4+Pjc+g+aoTea+b5Aw6HGJh+VOgNO9nZ4uNhTqTOw4sEw5sY1r2yzlM6qagz8Y95QxjRzCbKm5woCrIoPp0hTx+FsJRdKq8T04V6OEl4cH0h/Ow3hYWH4+vpyKFvJC2kXCZG68fGCWPGkACxVgKv2WVpRv5scxbnSKg7nKDmSreJ6Y5caf3cnnh4/kMkRfvi7O1tKkL+8wvYfrv3EeKN5Uc1/TRnMot/0o0ClE12Ns8qrkUjg7WadbjvyyC4vL4/nn38ejUaDk5MTK1asYObMme265s+h0+kYO3Ys//rXv5gxYwYVFRU4ODjw/fffs3LlSg4ePNiay3Rv4Vs5c+YMzzzzDKtWreL06dPs378fqVTKrFmziI+Pb5OwrC4qcrkcjUaDu7s7UqklO6wdvcxaYK0RGDp0aLtPGkxmgbMlWg7nqNh/UY6qxlIHP7K/F5Mj/An27cVru7MxmMx89HDLZJkfSyzLdHdnBzY8MkxsNAmwO6Oc13fn4Ghv6S1gEizVhKP6uODb242dmWqig3rz4cMxYikxNGX43Rfqy98eGiI6/1qW8PlsOHmVQX6uNBhMlF1vQALEBHlgEiDrWrUYobeKXm8y819pFtvsGysVdXoTT2/L4McSLatnRZDQxuV9a9DpdMydO5ff/e53PP7441RXV6PT6ZDJZLf+4zZgMpmIi4sjLy+PZ599lpdeeonRo0eTl5cHWMxlp0+fTlZWVmsu1zOEX19fT11dnWheKAgC2dnZpKWlsXfvXjw9PZk1axYzZ87Ez8+v1ft5QRCorq5GLpejVqvp1asXUqmlF93t9jazbkdut0bgVuO9XFHDoWwlX2c3HZM52EmYP7IP80YGibn3x3PVLEvNIsDTmY8fGdYiPXd/lpxXdl0mIsCddQmDybl8Ga1rEN8VVXM4W4XBLGAHxAU4Mi3KnweH9cfbvRcfNlYI3nj0JggCa7/K49PTpTw0og9vzAhDAlxRWE4zNp8pFW2zg317MSXCn8kR/oRJ3cSGGzduXWr1Rp7amsG5q9dZmxQlBiI7UvT19fXMnz+f5ORklixZ0ilxIK1WS1JSEm+99RaPPfaYTfi3/WaCQH5+PmlpaezevRtnZ2dmzZpFQkICMpmsTZNAbW0tcrkclUqFk5MTMpkM/8ZGjK2htLSUiooKYmNjO6RG4FYUqnR8fVnBwUsKLsstKbmD/FwZ7OfKkRwVYTJLZZvvTZbpcf09WfPgQEoKrhAdHY2bmxsbTpbw7qF8ooM8CPFz41iuCo3OiB3g20uCsk5gSrgP6x6KbnH0Zl3CN+8lD5Yl/Ku7LrP/ooLHftOPPp4uHM5W8kOxJThodc59eFQQy6eFtMhsfHJrBhdKq3gnual0uSNF39DQwCOPPMIDDzzA0qVLOzX4++abb9KrVy/Wrl1rW+p3BIIgUFxcTHp6Op9//jkAM2fOJDExkaCgoDZ9ubW1tSgUCpRKJQ4ODuJKwNnZ+Se/KwgCBQUF1NTUMHTo0A7bMrSFiqp6DjdmDf5fkRYzIPVwYlqklCkRfozo78mnp0r566F8xoX4sHJKX4oLcomNjcXFxUXM158+RMqaRIuphckscPaqljUH80QzDIBQL4vv3vSYvmw5p2TXhYqfFBLpGx1wvrqsbGF7BZZ+dUu2ZpCrqBU9892d7Rkf4stvB3uz42w5WWXV/GV2syafHSh6vV7PY489xrhx43jhhRfuuOiVSiWOjo54eXlRV1fHtGnTePnll9m4cSOzZ88Wg3sxMTE888wzrbmkTfg/hyAIlJeXs3PnTj7//HPq6+uJj48nISGBgQMHtunLrqurE7MGJRIJUqkUqVSKi4sLZrOZ7OxsJBIJERERd0XugFZn4GiuikPZFuPNBqMZdyd7avQmRvTz5E/T+iAvLWbYsGE4OjmxtjHSPnt4ICubdcY1mQVW7bP0lls4ui+JMQEcbgy45TSz8Rri78SLk4MZFRKAnZ1di2q5G0t1mx/LvTUrghlDpXxfUMnhHBVHspVU1hmxt5Pw7uwopkV2vOiNRiOLFi0iLi6OV155pVO+r4yMDBYuXCj6SMydO5cVK1ZQUFBASkoKGo2G4cOHs3nz5ps+WG6CTfitQRAEFAoFn3/+Oenp6Wi1Wh588EESExMJDQ1t05ff0NAgTgLWakI/P782X6ezsJT/ajjUmNxTbzTjYg/jQ/2YHOnPiTwNezLlN420v/aFxZL6xnwBvdHMs9sz+K6gkiAvF8qv12MWwMcFRvXpRUkNXFbU/cRXr1Zv5NltmfxQouXPzSL0YOlU9PstF7h4rYp1c4YyNdJSc9HRon/yyScJDw/njTfeuCu/r1ZiE/7toFar2bVrF+np6cjlcu6//36SkpKIjIxs9c2g1+s5d+4c7u7u6PV69Hq9WD/g7u5+6wt0AUUlpXxzsYwCQ2+O5mpQ11ocaPt5u/DEmP5MivDH182pRYPK5pbUYEkgWpaaxfE8Da/eH8rv7u0rJiJ9dUnBiXyLI5CzPfwmyJkpkVKmxfRDYmfPU9sswbo1iVHERzdFzavqDfx+s6V5xro5Q5gS0fGiN5lMLF26lKCgIFavXv1rFj3YhN9+tFotu3fvJj09neLiYqZOnUpSUhLR0dE/e8ZfV1fHhQsXCAkJEW9Ia3qoQqGgrq4OPz8/pFIpHh4ed8VNVlJSglqtJiYmBnt7e0xmgfOl1/nyooLjuWpKtfUW046+vdHWGclX6URhW9HpTSz9zLJMf2NGOHPjmpKoahuMPL3N0p9v/qggrusMHL2ipkZvwtHOYuNVaxB448EQ5o5sWvo3F/17Dw1hcnjHi95sNrNs2TI8PT35y1/+0mkJXHcQm/A7kqqqKvbt20d6ejpXrlxh8uTJJCQkEBcXJ94s1dXVZGVlERUV9bNJRCaTCZVKhVwup7a2Vqwf8PT07PRJQBAECgsLxcDjzW56QRDIkVuO3g41M+2ICnRnSrg/UyL9kXk48fS2TM6XXmf1rKbsOUDsjptZVs3apEgetB69mcwcv6Jm1f4csS2ZvQQifOyZEOLFfVFBvHmwkOyKGv72UFPHnI4W/fLly7G3t2f9+vXdQfRgE/6dQ6fT8eWXX7Jz504yMzOZMGECQUFBFBcX89Zbb7W637nVU0Aul1NdXY23t7eYOnynJwFBEMjNzcVgMBAVFdXq9yvW6MSYgNXHzsPFnpoGE8/dF8yScQNu2uPur80Ccjf+t3VzhuDv7sThbBVfZyso1tQjwTIRvDren/gRA8SU1o4U/euvv059fT0ffPBBdxE92ITfOdTX1/P666+zdetWAgICGDVqFElJSYwZM6ZNiT7WQhG5XM7169fx9PREJpPdkfoBQRDEts9hYWG3Pckoqhs4kqPiwEUFPzRWE8o8nJkc4ce9A7341/Fi8lW1/G3OUCY263Gn1RlYvPk8ecpa1j80tIXxqLbOwCOfnKVYU8ffZkcS5WUWK+j0ej0DBgxgwIAB7ZoYBUFg1apVKJVKPv744y45Yr2D2ITfGeh0Op577jnWr1+Pk5MTR44cIS0tjVOnTnHvvfeSmJjIuHHjWltZBVhuTGs5cWVlJR4eHshkMnx8fNp9k5rNZi5evIirqyuDBg3qsJXF9TrLvv1QtpLv8jXUG81IgNHB3jx8TxC/HeSDi6Ol9Hbx5vMUqur4x7yW/e+sE0K+Usff5zZ5BlqDpX5+fjQ0NIgTo7+/f5trKgRBYM2aNRQWFrJx48buJnqwCb9rMRqNHDt2jNTUVE6cOMGIESNITExk4sSJrT2PBZrqBxQKBWq1Gjc3N2Qy2W3VD5hMJjIzM/H29v5JGXBHotNb2pIduCjnRH4l1Q1GejnaMTrYm8sVNVTqDPxzXjRjBjfVUWh1BhZtPk+BUtdiQrjZnl4QBLGwSqPR4ObmJhZW/VJWpCAIrFu3jqysLLZs2XLbqdd3OTbh3y2YTCZOnDjBzp07+eabbxg6dKjoKdCWnH1r/YBCoUClUuHi4iJOArdKAzYajVy4cAGZTEbfvn1/8Xc7Er3J0p/vUI6SQ5eVqBobao4Z5MOUCH8mhfthbweLNp2nSF3H+/Oi+W3jhNCaQJ61jNY6Mf5cJqUgCLz//vucOnWKHTt2dGja9NWrV3n00UeRy+VIJBKWLFkiVvTNmzePoqIiBg4cyI4dO8SakjuITfh3I2azmdOnT4ueAmFhYSQmJjJt2rQ2n/Fbb3iVSoWDg4NYP3DjtsIaFOvXr98ta//vJGZBIKO0iq+zlRzKVnG1sg4JEOTlgrJGzz9TosWy4duN3lszKZVKJYIg4Orqitls5uTJkxw+fJidO3e2acXVGsrLyykvL2fEiBFUV1cTFxfHrl27+N///V98fHx45ZVXWLNmDZWVlaxdu7ZD3/sm2IR/t2M2mzl37hypqakcOHCAgQMHMmvWLKZPn94mTwGwxBqsN7ydnZ2YOgyWxiPBwcFd5jJ8MwRBECv0PJwdiAxwZ9RAy9Owo47s9Ho9P/74I8uXL6ewsJAlS5awYMGCO97pJiEhgaVLl7J06VKOHj1KYGAg5eXlTJgwgZycnDv63nQn4R84cIDnn38ek8nEE088wSuvvHK773/XYjabycrKIjU1lS+//BKpVEpCQgIzZsxos1lHfX09CoWCiooKampqCAwMZODAgR1eCnwn6MhzeoBNmzaRmprK1q1bOX78OABz5sxp93V/jqKiIsaPH09WVhb9+/dHq9UClonO29tbfH0H6R7CN5lMhIWF8fXXX9O3b19GjRrFtm3biIqKut0x3PVYj9vS0tLYt28fnp6eJCQkEB8f3+qnttXlJyQkRJwImluPu7m53foinUxHi/6zzz5j48aN7Nu3r1M+b01NDffddx9//OMfSU5OxsvLq4XQvb29qaysvNPD6B7Cv7EW+e233wbg1Vdfvd0x/KqwuganpaWxZ88eXFxcmDlz5i96ClgzCG90+bFaj8vlcvR6PX5+fshkMtzc3Lo8dbijRZ+ens5HH33E/v37O6WngsFgID4+nvvvv58XXngBgPDw8G611O/UFKeysjL69WvK3+7bty9lZWWdOYQuRSKREBoayquvvsp3333H//zP/2A2m3n00UeZPn0677//PqWlpVgnY61Wy8WLF4mJifnJDe/o6EifPn0YPnw4cXFxuLq6kp+fz+nTp8nNzeX69eu0Y1K/bTpa9Hv27OGDDz5g9+7dnSJ6QRBYvHgxkZGRougBZs2axcaNGwHYuHEjCQkJd3wsd5JOfeKnpaVx4MABPv74Y8CyZzt9+jTvv//+7Y6hWyAIAteuXRM9BRoaGhg2bBjZ2dmkpqa2aT9vrR9QKBTU1NTg4+ODTCbrlPqBjhb9wYMHWbt2Lfv27cPX1/fWf9ABnDhxgnHjxrUo3vrzn//Mvffey9y5cykpKWHAgAHs2LGjzfGa28C21O8pCILA5s2bWbFiBWFhYVRWVjJjxgwSEhLa7AVgNpvF/gNW63GZTIaXl1eHpw53tOgPHz7Mm2++yf79+++qE4xO5o4Jv1PTnUaNGkVubi6FhYUEBQWxfft2tm7d2plDuOuRSCSUlZVx9uxZvL29UalU7Nq1i9deew2FQsH06dNJSEholaeAnZ0d/v7++Pv7t6gfyMnJoXfv3mLqcHsngY4W/fHjx1m5ciX79u3ryaK/o3T6cd7+/ftZtmwZJpOJRYsW8cc//vF237/HUVlZKXoKlJSUMG3aNBITE3/RU+Bm3Jgm6+HhIabJtjV1uKNF/91337F8+XL27t1LUFDQrf+ge9M9lvo2Oo6qqir27t1Leno6ubm5TJ48mcTEREaMGNHmSaCqqkpMk3V1dRX7D9wq/72jRX/mzBmWLVvG7t276d+/f7uv1w2wCd/Gz6PT6di/fz9paWlcvHiRCRMmkJiYyD333NOmJ3jzXHmVSoWzs7OYK39jPnxHi/7s2bM8++yz7Nq1i+Dg4HZfr5tgE76N1lFfX89XX31Famoq586dY+zYsSQmJrbZUwBubj1uTR3uSNFnZGSwZMkSdu7cSWhoaLuv142wCb893GUVV52GXq/n8OHDpKamcubMGUaPHk1CQgLjx49vc0VbXV0dcrlctBkLDAwkODgYFxeXW//xL3Dp0iUWLVrEjh07iIiIaNe1uiE24beHu6ziqkswGAwcO3aMtLQ0Tpw4QVxcHAkJCW3yFLAu7/v27YvRaESpVGI2m/H390cmk7W5fiAnJ4eFCxeydetWhg4dejsfq7tjE35H0sUVV12OyWTi22+/FT0FoqOjb+kp8HN7er1eL/YfsPYesKYO/xL5+fksWLCATz/9lGHDhnXkx2PRokXs3bsXqVQq9qj7la7ubMLvKO6Ciqu7CpPJxKlTp9i5cyeHDh0iLCyMpKQkpk6dKnoKtDaQZ60fUCgU1NfXi5OAu7t7i5yD4uJiUlJS2LBhAyNHjuzwz3T8+HHc3d159NFHReEvX77817i6swm/I7hLKq7uWsxmM2fPnhU9BYKDg5k0aRJ79uzh3//+d5taShuNRjF1WKfT4evrS11dHV5eXqSkpPDhhx8yevToO/ZZioqKiI+PF4XfRUU27aV7FOl0JQaDgdmzZ7NgwQKSk5MBkMlklJeXA5Y4gDVi3VOxs7Nj5MiRrF27lnPnzrFs2TLeeecdampqWLp0KZs2bWr1xOjg4EBAQAAxMTGMGjUKT09P1q5dy9ixYwkNDcVgMGA2m+/wJ2pCLpcTGBgIQEBAAHK5vNPe+26kRwi/p1RcdSR2dnacPHmSDRs2cPLkSd555x0qKipISkoiKSmJTz75BKVS2apr2dvbix2Ot2zZwuLFi/nss88wmUx3+FPcHIlE0uWly12OIAi3+/Or4dtvvxUAITo6WoiNjRViY2OFffv2CSqVSpg0aZIQEhIiTJ48WVCr1V091Lses9ks5OTkCKtXrxZGjx4tTJw4UXjvvfeE/Px8oaamRqitrf3JT1FRkRAXFyccOHCg08ZZWFgoDBkyRHwdFhYmXLt2TRAEQbh27ZoQFhbWaWNpB+3R5y/+9Kg9vo2ORRAEioqKxHJie3t7Zs6cSWJiIn369EEikaDRaEhOTmbFihXEx8d32thu3OO/9NJL+Pr6isE9jUbDO++802njuU1swT0bdzeCIFBWViZOAnq9nqlTp7J3715ef/11kpKSOm0s8+fP5+jRo6hUKmQyGatWrSIxMbEr6unbi034Nn49CIKAXC5n/fr12Nvb86c//amrh/RrxSb87oLJZGLkyJEEBQWxd+9eCgsLSUlJQa1WExcXx6ZNm9rUostGt8Z2nNddWL9+PZGRkeLrl19+mT/84Q/k5eXh7e3Nhg0bunB0NnoKNuF3IqWlpezbt48nnngCsCyJjxw5InrDL1y4kF27dnXhCG30FGzC70SsCTFWowy1Wo2Xl5dYLtvTXIdtdB024XcS1qKRuLi4rh6KDRuda7bZk/nuu+/YvXs3+/fvp76+nqqqKp5//nm0Wi1GoxEHBwdKS0ttPnM2OgXbE7+TePvttyktLaWoqIjt27czadIktmzZwsSJE0lLSwNsacM2Og+b8LuYtWvXsm7dOkJCQlCr1SxevLirh2SjB2A7x7dh4+7Fdo5vo2dy4MABwsPDCQkJYc2aNV09nG6D7Ylv466lJ7ZVvwHbE99Gz+PMmTOEhIQwaNAgnJycSElJ4YsvvujqYXULbMLvgWi1WubMmUNERASRkZF8//33aDQapk6dSmhoKFOnTr0rLMh6elv1O0l7lvo2fqVIJJKNwLeCIHwskUicAFfgNUAjCMIaiUTyCuAtCMLLXTzOOcADgiA80fj6d8C9giAs7cpxdQdsT/wehkQi8QTGAxsABEHQC4KgBRKAjY2/thFI7Irx3UAZ0K/Z676N/2ajndiE3/MIBpTAJxKJ5JxEIvlYIpG4ATJBEMobf6cCaL2l7p3j/4BQiUQS3LgySQF2d/GYugU24fc8HIARwL8EQRgO1AKvNP8FwbL/6/I9oCAIRmApcBC4DOwQBOFi146qe2Db4/cwJBJJAHBKEISBja/HYRF+CDBBEIRyiUQSCBwVBCG860Zq405ie+L3MARBqACuSiQSq6gnA5ewLKEXNv7bQsB2btaNsT3xeyASiWQY8DHgBBQAj2N5COwA+gPFwFxBEDRdNUYbdxab8G3Y6IHYlvo2bPRAbMK3YaMHYhO+DRs9EJvwbdjogdiEb8NGD8QmfBs2eiA24duw0QP5/yeyH7TCXrezAAAAAElFTkSuQmCC)
%% Cell type:code id: tags:
``` python
#Generate the response of one RGC
x = np.load('D://inception_loop/RGC_sim_data/data/static27012021_original_images/data/images/'+str(555)+'.npy')[0]
RGC_response(rf=rf, image=x, plot=True)
```
%%%% Output: execute_result
2
1
%%%% Output: display_data
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAASMAAAI4CAYAAAAh2Y/oAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABMUUlEQVR4nO29eZidd3Xn+f3WXWqVVNotS0ayscE44AWEsYEmYHDGIQS789AZHEicjrs9mZAZ6DBDTBYSJsnTzjIGupOG8QRip4dgiFnsdghgFJO0aZCR8SoJW/IuWVJpK6n2usuZP+4r655zq+6tUqmqfqr6fp6nnrrn/b3Lue9769T7fu8550czgxBCzDdt8+2AEEIACkZCiERQMBJCJIGCkRAiCRSMhBBJoGAkhEgCBSMhRBIoGIkpQfI5ku+abz/EwkXBSAiRBApGYlqQ/FWS3yf5KZL9JJ8h+eZs+Ysk+0jeULf+z5F8mOTxbPwPw/5+heTzJA+T/P36OzCSbSRvJvl0Nv4Vkivm+C2LOULBSJwKbwLwGICVAP4OwJ0A3gjgfAAfBPCXJHuydYcA/AqAXgA/B+B/JXkdAJC8CMB/AfABAOsALAOwvu44/xuA6wD8NICzARwF8Fez9q7EvELVpompQPI5AP8OwAYAv2tmF2TLX4daYDrLzA5kyw4DeKeZPTLBfj4NwMzsP5D8BIDXmNn12VgXgH4A7zaz75LcCeA3zWxLNr4OwAsAOs2sPItvV8wD+fl2QJyRHKh7PQIAJwJR3bIeACD5JgC3AHgtgCKAdgB/n613NoAXT2xkZsNZIDvBRgBfJ1mtW1YBsBbA3tPyTkQy6DFNzDZ/B+AeAOeY2TIAnwPAbGwfandaAACSnag9+p3gRQA/a2a9dT8dZqZAtABRMBKzzRIAR8xslOTlAH6pbuwuAD+fCeBFAH+Ik4EKqAWuPyG5EQBIriZ57Rz5LeYYBSMx2/wGgP+L5ACATwD4yokBM9uOmkh9J2p3SYMA+gCMZat8BrW7qu9k2/8QNfFcLEAkYItkyL6B6wdwgZk9O8/uiDlGd0ZiXiH58yS7SHYD+AsAjwN4bn69EvOBgpGYb64F8FL2cwGA95tu1xclekwTQiTBjO6MSF5D8kmSu0nefLqcEkIsPk75zohkDsBTAK4GsAfAjwBcb2Y7JtumyHbrYHf9XsIa0Zc4HkbzPmfTOgreDqG2UvT7s1zYX8jprbYH71qF7mqwi2GBTfB+qs3fY8M5CftoK7XYOuw+P+bttpLffzXPYKOp3UB8O9O9xC0+jvGaIRfPjzfbRv0BCoP+mlSL/qKWepyJXLv/UFSrfn2r+P3nChVndxfG/Tgb32AlXKSxij/JpYp/0xbWjz403mI0/wy1pFWIaHVNyycXlI8eQWVoaEIHZpKBfTmA3Wb2DACQvBO15/9Jg1EHu3FF/n86uYDhrFn4443jgdyaVc4eu+AsZ5eW+Lc3sN5f1NJSf07aj/qzevw8f7xyT4w2wZ8h7291/aizbaLAM+gDaMPnJu8XcNzvo3tP+KCGQ8Tg0fu0fw9dB/wfy8jqordXMox7m/5vD9VitIP/Yf34wW0rhT+0Nr/9+DJv25LwHyT8YS7Z6c/vugcGnD1wbrezD1zpd7fsvKPOHhnzb3D0SIezl6877uw3nvWCs3sLI4j0lzqd/cyA/1wfGPARcnTE+1Aa8DY7/ElmCIATfg7rxyvx77Lp6kAprB+Ol+8/+SHc85lPTbqbmTymrUddKj9qd0fr40okbyK5jeS2ko3FYSGEADAH36aZ2W1mttnMNhfY3noDIcSiZCaPaXsBnFNnb8BMixdbPJY1kA+PKLnw7L3U76/S2fwRI9q9T3r78CVBT+n0jzzlZWEHJe8fcxM85oVN4mNZy+fxKEuFx6S4fnxs6z/f/4MoDPoDtoWnoKhRNWg44Ra9LTxWxse2crd/A/Gxq63gx5ct8Y85hbw/gUOj/gQMnbPE2X1v9DYr3p/enc7E2P6Vzl71VNCEXhx29t7f8SekM+dPWLna+BmvBjHyyLB/bFvdM+TsvvAsXhoOWulg0FLjZyrqbJFWOp7Fv6PmWmz9Z6jZA+JM7ox+BOACkudmdUXvRy11Xwghps0p3xmZWZnkbwL4NoAcgC9ktUZCCDFtZtTPyMy+CeCbp8kXIcQiZo6bq7G5LhS+2rdq+Fo4Fx5GK379+PxfLYRn26B/tB/363f0ez2g3OF9zY0F/SPq8R0xNSF8pTrUeLpbZny0eH4f7/ULciPN9ziy2r+n/LDfvtQdtg/Ha+/3C8Z641fxfv1Sb/iaudenEnR2ersYNKC2Nn9OOwoh7yfqF8H9wnqvt/R3eD2mQQ9p98fLha/J97wqF8a9/Y6zn0Ezjpc7G5YdHW9cVs9wqdB0vK3ofYzf3K/7tv/cFQf9+sOr/PjI2qC9rgifkaiNxnMY/wzqV28iV6k2TQiRBApGQogkUDASQiTB3Dfkb6t7Hg2aUNST2jr9szK7uvz6uaDpjHo9odLuRZ3cWHONyNqCJhQe1Xue8/bAeV4vKMfQHvI5WJ5Az4l5OnEfDbVModZqoHktWdRwxn2aDdp9tQPGlodaqxF//NEVfnx0TfM8oVyHt4vFMB5yr6JGVAjjUSMaL4dcrnC+ujt91r8F/3u6fMnO6Li/6Bt6j/n9Ffz+2sLxon99o/6Ej1Ya/+TiPrqLPjfp0IAvWSmN+30UQv1cKepmoUyp69ljwQ4OtahXtWLQmDb4cpVjG/05LNf92TaUA9WhOyMhRBIoGAkhkkDBSAiRBHOqGTGfR27VyVqfytm+7qdt2OecxFozjPln46jxcMw/kOZDt4bOQ0FTCnlEsZdPzrvToKcUg95iubg/Px7rtACgsmq8YZnfSdgm9OephFynWEvWqvvDwMaYUxI0oM6QJ9Qe8oBC3VM+5AEVQn+fQs7budAiJGo+0Y5ETamj4DWgSqgFi/4saQ/9hoI/5SC6xTqyle2Dzj5W8i1FxoNGNBXNKNrjY16D6Qi5WT0dXsfqe8b/XbUf9ueEI6F7BmOuWPNctbYxr2l17/R/aF27YqOwk/vbc2zyiYB1ZySESAIFIyFEEigYCSGSYG7zjPI5VFctO2lX/fN+pTs24/FYtxdI2spB3wgaU+9P/PN8eYnff9t42D7vY/Pxjf705P2jN8aXerv9kN++3BOe/ZdPkGQxHv4fxF4zYbzY7+2Ow3712MO52h18WBk0oK6Q95MP5yS4mw+1Y53tXj/IT9SzqY5iPhyvpSYU6q5iD/Cw/Wg51GGF7VtpUt1Fr8fEurB4PIyGlrBBE6oEjWmifkaRYyNed4o+LuvyGs2BI/6DyJBodOwCn5+3an+/P2AlnKOoU0ZNqdAibIT1xzac/Juv7o2JdSfRnZEQIgkUjIQQSaBgJIRIgjnVjIxsqGtpSsh3aNCIYuudsH7UgKIAkgv7CxIWend7/eDoBV5z6vKtchrykophjq72o42aUX7Ab8RS7KsdEolCU6aBDX68/6fCMUKeUOwpHYkKTqwdywdNKc75Fa9uPmg2XQWvMbVhenlFrai0hXrFUOu2tNMLfyNBE+oImlZ78D9qSHFOs3i8OD5eadRMSmFqoJGQV7Skx2tEUWerHA9TFYWP/dDZ/hr5iZDQ2ASqEusNveZk4Rwz1LId2rzC2Qf/1clrPr4rHvwkujMSQiSBgpEQIgkUjIQQSTDH/YzMCzPhWbVBLmgxRXisBWsg1rYFqi22j3rO2v/6E79CzM8ohjypWOOTn/7pbu/3eSzPXe8n7R0+L2hOxdAXPNZ+tTWvgyoUm88tX4pzwQWBIs4D3xM0l6gRtQeNJh8a3jTUaYWCv7GQVxTzkjrzYaK3QEPeUKCxx3bz8xftnqKvAzsyEnpyobGHUiTmbh0e9PuIfbLi3HYNPa5CnhDHQk3oUt8/yaKmFIj5fYMbvN3Ze1KnY35yzVJ3RkKIJFAwEkIkgYKRECIJ5r4Hdv0j4+RlKgAaNaIGDSauH/OMYr+jkA/B0IM7NxzqpvYdcXa1VW/gst+eUSOqTPC8XAx6QdChDr31bGcPnxvydGJ/oVDblgu1ZNWQ05IP/X3ivGWAt2PeTDXUMZVDT+rjw77OamnQUKarEcXarmroCR41qbj/PP01GK96f/MhT2hlh08mi/2J4vbFtpBXlfd6TEeuUcM60Ob7ZI8HHawSzvFwv59nLUo6baXmGtLw+T4PaP8V/nilpUH3et6f8+59/hx17/G5W9VLBpz9qtWHTh4rr35GQojEUTASQiSBgpEQIgnmXjNqkftTD4PG0hbruIb9s6o12MN+h0GPsRYaUFRPmAsiV7BbjU+YZzQeHujP8pVDhy8JPkbZKZYVlZvn/cRkrjge66KixhTzbizoGVGDiv/vYp5NnIesFVGTaUXUgFoR85KixlQM/sbzEfsZRU2pbI1CaayHi327B4eCRjTY/M829lrPhz+DF38pzGXXHnpkh3nZhlb77fNLvY62tNvbFxdD46/6fTfpyq47IyFEEigYCSGSoGUwIvkFkn0kn6hbtoLkfSR3Zb+Xz66bQoiFzlQ0o9sB/CWAv61bdjOALWZ2C8mbM/u3W+2IVUPb0Mlnbh497sZt2PdtqY6FedIb5ndqEUtjg6LYh6X51o21c1EDir2A28K4heOXJ8ixWO1zPp653mtGlZ6gkYQ8oupY0K1iv6HQQ7stjMe8mVZ5SpWgJ7R3eI2lI8wTHzWpctCgoobSk2s+l330NxI1orj9WNi+I2hEPXl//CVhHrZy0IDGQh7UsiZ6CQC8NNTRsCzmbkXNaGwk9OEOHyNWwnx/QYZbud0vOL459Irv8aLS+p5j/njhHB4b8xrWcKl57/q+wZP1lfH8u+M03QsAM/sXAEfC4msB3JG9vgPAda32I4QQzThVzWitme3LXu8HsPY0+SOEWKTMWMC22vfjk35fR/ImkttIbhuvDE+2mhBikXOqeUYHSK4zs30k1wHom2xFM7sNwG0AsLRtpVWfeeHlsQYNJmg8bCi6CbEzakKRuH4rDSjWvsVmwq2IGlGshevxfWIAYN9Va5w9vrzFe4rzrHUEjSQf84hij2m/ebG9ueYT+/d0hByT2I+5FPSOoTGvJ4wHzWmk7PWQqNl0h9quaEcNIub5HB9vPgdZb7vXKaNGFIn6ScwzGqn499MZatEm0kwGwzkaCTaPBDvUnhWPeXvdD/17aBv1ut/qLT7Xa+x93qenj3jdshB0w5gXdWTI72/oW2d5/4ZOnjMemYFmNAn3ALghe30DgLtPcT9CCAFgal/tfwnADwC8muQekjcCuAXA1SR3AXhXZgshxCnT8jHNzK6fZOidp9kXIcQiZk5r04gJdKB6Zpg31FITaqURRVr0/p0u1WU9DctGzoq6WFgh9rQO/ZDjPGhxnrNK6C+0LGg+sX9RzHFZ0u41nEKoDYtzyR855ud9Hx4I874NeE3lpeDvYJfXR9Z2DTo79gcaLvv1Y3+jqNEsb/dfoqwq+vPRFmrRYh7UQMipaXX8OD5abvyTixrR6FGvcxVGQx7RiLfP2hpys8aa1+8t2+XPwbG7fK7b4df59bv2++O1P+v3v+FZf41YOjrpsfMjk/umchAhRBIoGAkhkkDBSAiRBHPfz6ieVhrQdDWiWCs23TyhaszdDHZler10Yu3awSsb64krHeEYYQ6s6EKrec/KYV6zrm6vJ0xXI+oueI1mKGgmz+1f6eyO7b5uyVuNvXX45DJnV/u95vRSzudhDa0LtW6Xer3i7BW+rirWxnXkfI5M1Igih8Z8bthgyWtgsS5rrNK8sft4uXF8dNDvo3A41AuG/kRrHwx90Ev+PTT0fo+f62BGDan3qZCbVok9tabXI4r1fzdNeojpzkgIkQQKRkKIJFAwEkIkwZxqRgZfK9U05who7D8Ue0i3yhOKPa+DzQ7//M8ur3BYpx+3rmCHvKVyrx8f2ODtOAc5AFh4gLeQNxSbLkUNJPa8bth/XD9oRKu6fJ5NnNfr+WNe5xp8xGtEa3aE/kelcM6bt+BGW8kvKB73mk65y5/jwpD3/6WNPidnqMfXmp3d43tmxdq2yKExnwu2d9BrWseGw2ckvJ+uUOuXz/nz0d/fWJ/Y/mKoPQt9xdc8FObzG4taathhNX7Omvd6j73m4zVqWL/cvAazwZ9Six1OspkQQswLCkZCiCRQMBJCJMH85hm1qkWLPa9jnk855Nws73V2ZZ2vuRkPmk616I9fGPDP5vnjoS/MYa8/IO/1jFyfr8mxtnXOPn5uqNMC0DCNViG8p0Kc682fk2rIKyp0es2iM/QnihpK7Cf08LMbnb1iqx9ffSjUxoVrEHNcGv7dheFqIdSSLffHi3JHqSvmovnz011srgnFfkOxH9GewV5nHzrmNaRY+7e0y39GcqEH9+HjXiMqPtvYA7sYPlbLn/Kfw8KgtxvyhgLWMDdh0F4rzfOEqvkWf5d5r3G1jQf/SrFv+9TueXRnJIRIAgUjIUQSKBgJIZIgrdq02G8o5i+E9auXvsrZg+tCDkgItflhv33Xbj8Dk+3d7+1xr7dUi15viP2SqiNePygc8K3Bl2x4IyKH3hBzRprPixZh6HfU0ek1k5Uhj+hY6An90g/PdvaGh6OeMM06pKhnRN0vrB+vUfz3GEvHOvqDRvR86Km9KWhcnb52rRwOODzuPzOHB30/59gzuxpyeI4O+PXL4+EzEa5f7wE0sHJHmOvemutwDbVnQbdrmM++xf4QNKKG/cfatEiD9hs0oxYa1cu7mdJaQggxyygYCSGSQMFICJEE86oZsRjm6A4aDJf6HI+jV6x3dv8FPpaWloTn+3Zvr/1h6PUT52nrbMwBqSdqSBYnPQ/zplmDftK4z1iHhBHvo3WH3Kd2/zzeFvJeRoZ9LtNTx/xkv+1P+/e4ekfY33hznxs0HsTxoG/E99yiVq0acmRyLXJqwlT3KIZasJhHVA61eceCZjR8NPRjei7MaXa2vx5da70m1x5q0wZCP6SlL4bPDCbIC4pts+IpaJ5G1Lj/0NerQZNqUSNabffbt4037+vVtHatyeXUnZEQIgkUjIQQSaBgJIRIgrmdN62tDW11uoyV/PNz9UJfF/Xce0NdUJg/yoL3y197yNkH9/Y6e/9bfOzd9/bVzm7v83OEr3o01CE9cdjZlaee9g409HXx+k/3Pq8nAMDwWV6TKIem0ZVB73N5eejh3BXO4ag/Kblj3oexVf55/9Dr/PjSZ/176DwcauNa9ZBqRQv9o1Uvndj/qL0/jIcdRDvf1jznpbjf5ymt2Bn6EZX8+a2u8W9gaMBrcmv+OeSmVRo1o2re76OhnK/aIk8orh97YocaTAsnva2h17w3W+U5teqJzXL9OVQPbCFE4igYCSGSQMFICJEECkZCiCSY46RHcxMlDv7Ma93osfO8mDq+1gu+HA3N0Nb45uvDYyGJMibsdVaajo95PRt7f84Lcy+9wzejf82f+wkPyy/scXZumZ+QcHhNaBwGIOdrJBsKQ4vHQ0P9434fpaVBsA76YG4kbN/lVygt8wc8+GZ/jop9fv9nbQ3nMNZkNtT1RvEzrN9KD28haBcG/QrDpcZzPB3i+T6+qXmhcnXHEmevDhMgFgfCFwANjc/Q2HAuJCHmRluc8ziZaThGTHKstDd/T638s1BY25DkeIrozkgIkQQKRkKIJGgZjEieQ/J+kjtIbif54Wz5CpL3kdyV/W6cSF4IIabIVDSjMoCPmtmPSS4B8BDJ+wD8KoAtZnYLyZsB3Azgt5vuqLcTh977Uy/bA5v8w+jo2lAU2h2ay6/2Gk1Hwa9fCUWnG175krN7273GlA8CzdKCF3AePewbjx044jWgZz94jrM3fsonRQ6+7QJnl3wOZ42Y9NeYE+fID3u7eMzvID8SioVDg/8le0KSZEi4O3yR11zGVgbd7F95vaHzgN9+2TOh8LYc/YkCiTdj1mNuvEWS4qAfPx6SEmPSYySOj6z2dseheH799qseDxcsSkKxZjTfqBnFJMKoEcXxhnMYYRTmghkvQUiKzA+FiTS7/TktDvm/Q8Tma+ONyb1ToeWdkZntM7MfZ68HAOwEsB7AtQDuyFa7A8B1p+SBEEJgmpoRyU0ALgOwFcBaM9uXDe0HsHaSbW4iuY3ktvLo0ESrCCHE1IMRyR4AXwXwETNzMz2ZmWGSL2HN7DYz22xmm/MdjfOMCyEEMMU8I5IF1ALRF83sa9niAyTXmdk+kusA9E2+hxqFVWPY8Gu7X7bP6hxw48dLvsjwyJhvdv7Uj1/h7LbnfSytXuUnUXzhiNfUX4C3uzt88/rDh72o07PMCwQXb9jr7O05P0njsfde7OzRXv9wPr6s8Vm/LTxeN+TpBM2hEvq/xf5uUU9oC3MaDqz3l7ww5P+HdO2Pmo0/xwPn+wMOBp1vcKPXnNb+IOgRo82bt8U8q4Zm80GfKIT1hwb9CRoNk0J25f0J6cyHCxAmRGgrewd7d/n1Kx2h8Vjwt2GCgYma28feamGfjOlxIa8nnrP4GYi5TXH9tlBYWwnN1I5t8p+ZNVu9cMk4ccVy/3c0uubkNakemjzkTOXbNAL4PICdZnZr3dA9AG7IXt8A4O5W+xJCiMmYyp3RWwD8MoDHST6SLfsdALcA+ArJGwE8D+AXZ8VDIcSioGUwMrMHMGH3ZgDAO0+vO0KIxcqc1qZVqm04MnpSxL5oqZ80cXUxakg+j2fprqCHhOfzoYe9JrTiSr//IwNeQB8a9bVsl577orMf3nGusx97epmzO89zOj4GXuGferv3+mfxiTSjaiilqoTyuoZ/A0FyiHlHDeMjsS7Jj4/1ep9zoSF/+9GgIT3uPzL9YRJK9nr94OBm/4bO+oE/fq6VhhQ0loFz/P6ihlYdbJ6DsyTkko1Xm9dpLX/Ka2KlJS0mGkXUZ6Lm1fRwE9LQfC0kCsXas4aJMxvq4eJnwn8GjrzGfyiHX++10yV7/d9RpRgmznhVqDF93ckP6fh2NVcTQiSOgpEQIgkUjIQQSTCnmtF4OYe9h3pftr8x6PNyLl7ra8me/Bev2eSXhzqh0Auo4CUn9D3qk8IvefMuZz/+kq89e/LgGme39YRJG0e94DLyjK9V4wr/PPyKu7xmlb+sMUl9eE3I4zm3YRVHsb95zkjUJKpxnsyQd1QKeaiVmKfUFvOOvN2z0x9g8CJ/gGqHX3/wbK8n9D4d8nyC/4df68/5WDjHsT8Sx/yCqAm1h8Ss5UWvh/Q8FyZA6Jhc4wAa9Zg4cWe8PhM102/QlVrU68X6tkrUhFrojA0aUhgffoMXIs8/66Cz9/zbXmeXdvq/g+7X+RrND2x65OXXn+2cvApDd0ZCiCRQMBJCJIGCkRAiCea4B7Z/PI3pEbv7Vzk76h3Dq/3z/vVv+mHTY93z7OucPTDuk1JescrXsu1+Lmg6pVAjFJ/FQ8pJg15zwD9rL/t+aIYDwP7VJmcPhlyltvHQ6ya0kin78r0G3azaIsck1oqVO8L6ofYtajTtx0I/oMM+R8VCP6WYVxVzaPpf6T+So2uCXhLfTrDj+To07EWxVy31JZRbfvJqZ5/9UuwlFPxt1T46yjFRU5pg+1if2CDixH3GntcWz1GLPKRgHn+FP+exBnO0Eia27Paa0vhl/kP5+tW+F/z+8ZOaUmmiE5ChOyMhRBIoGAkhkkDBSAiRBHOrGVWIyrGTQtBllzznhr//5CudXdzkcxJ+9cJtzn5tp382fX27z1P67dVbnf3lAb//vpLPj3jmJa9Z2bB/vs2NNq97snzIFyn401s676yGbToPecGgcNzrWtVYS7Y8aD7DoT9ynDouaD4xJyXW97WVmvesjnlGUX9oP+rXD6e4wZ+Dl/lz1JDX0yJnpmHOsWKoVwz1h0Nlf0JX/ZO3aT7JpxrygirtLT4DQRKJetBELbkbehxZc52pIbesIZepqYsodfsd5t9zyNmxL3il6ndYCf6VKn58x1H/OS/VbT9cnnxeO90ZCSGSQMFICJEECkZCiCSYU82o2A9s+sbJB9wfDr7GjefCs+r//qbvOPtDvb7f0NOlQWc/Pu5ryy4o+GfhN3Q85+yHsMnZ69f0O3vP0VCrFvWXoM8Uj4TYXvArtI00Top2+LU+Uajc0zwvJ2om5c7mekOs34uaTcxJaQv6RdSgogbS0H8puFMIc9fHvKFWGlBD7VmLOcCsw+cJLe3yJ2BF0euQ+bHm/hQH/f4Y9JOYRxRr++K8dRNPWxF7XjevB2ysNfN2rHWLPh66xNtv6PV/JwdHfX+i4ZL/EEYNKZbbFXL+nK3oOHnOd7ZNPg+e7oyEEEmgYCSESAIFIyFEEsxtD+wO4uirT4oMXfv8+IiXaHDHn77H2Z+5xtfEPPXTdzh7b8WPPzq23tkXFn1/obEgyKzo8NsfGGheJ4bzvWa16oHQkHmF75k9cL5/Fgca83jyod3L+NKQuxR30DDnVsMh/P6DhtSg0cQ5uMInJPbQbnAnbF9a0qL/UJQQoqYS+ilV4ic2ikgV/4aWtfs3PBLqrBpq3QLlzjgvmh/PBU2pEtYf7w560MSi0bSIuWHxPcR51ypFv8Jr3vqMs8vhopRjXlGwY16RBQfGyv4ideends+jOyMhRBIoGAkhkkDBSAiRBHOrGRWBwVecFAkqveEBPITG42M+iaX7v/veNJcv/TfO/oNX/zdn9+a8BvT146939nf3+V42KzvDHOKhP9HoRt/fmQc7nd3zw2ed3f/285wd+10DjfOaRV2qMBDqgJYFvSDWooUcFYbn+ZjpFDWkqGG1rHuKkk3QrCx8wqqhfq/a3iIPJxzPCmGetk5/kWLP7iMjPo9rqMuLXlF/iTk6jRO5eTPqM/H8dxz1/lWLE4hU4ZC5cf8eYy5Y1IjiPuM1OnyxHz+v4PsPDZb8ORmv+MSl4XGvs7Xn/Xs6Z9kRZ79z1U+cncPJ9/NovrGn1wl0ZySESAIFIyFEEigYCSGSYG77GTFoCKHHdKHXP8sWLvYaTtQ7jj2zwtn3nfVaZ8c5sl7V4fOMYi3c4aAvjG1oaE7s93e7Twoaf80GZ4+u8M/qxYHGHJOG2q5ArDuqBI0l5r1EQSHWRjXkSsX2QS165cT1S0uDphVr69rDDqb776+h73jQzMoh5yVc0yPHvM54sKcx18sdroVGFvWaWMvX0B9qzO9wIskk9gGPulO1EJOzgg+l5rlLpV7vw57BXr//8Kbi38WGZcecHfuIryr4fLuz8n79Ql1BZCEKsXXozkgIkQQKRkKIJGgZjEh2kHyQ5KMkt5P8ZLb8XJJbSe4m+WWSLR44hBBicqaiGY0BuMrMBkkWADxA8h8B/BaAT5nZnSQ/B+BGAJ9tuicDWKp7Pu3xz7K5vLdX9QRNJuQ/rLrQj3/n2Qv9+BI/nlszeS8VAOgf8nlDHPTHe+Xfh4nqw3xU48v86WzvD3VVE/RPHl/avDgq5uE09POJ01DFObLKzXNQ4vYxzaZhXrVOv7+R1XEe+OBvqBWzFrVZrAbNKxdFnJj3E4vpQn+hkdBju0U/olZzljVoaAGWW7y/2O8agLU3z1WK17xtPJ6T6GN8j378lUt9/6J8m9dxesLkfFF7jSwPBZXr834+wvG6CQZzTU5gyzsjq3FCoSpkPwbgKgB3ZcvvAHBdq30JIcRkTEkzIpkj+QiAPgD3AXgaQL+ZnQiZewCsn2Tbm0huI7mtMjQ00SpCCDG1YGRmFTO7FMAGAJcDuLD5Fm7b28xss5ltznV3t95ACLEomVaekZn1k7wfwJUAeknms7ujDQD2Nt8ayI0CvU+dtA+/PvRJGffuHBzwOSFrl/qJ5JcVfdJG34/9zdmRip+0667lfj6nFU8EvWOj9+cVj4Q5zX7i52mrbvL7Kx73z9ajy33Nz9jyRn0o5p2UQhpM8VjQUGJ/oTDPWH4sFi55s0maB4BGfaHcEXrVrIgCRtSImu+vVQMhiz2j47/LtuaaTKvxJ/eudfaaFpJRW+g/3WSq+Np4PopMQTOLDaMwBY1ozJ/Uxr7lfp9tI6F+r+y/W3r6uJ8f8PBQ6MNe9m+yVArzBwYd75OX+prQx8bOcfbR8smbkKGKz1GqZyrfpq0m2Zu97gRwNYCdAO4H8L5stRsA3N1qX0IIMRlTuTNaB+AOkjnUgtdXzOxekjsA3EnyjwE8DODzs+inEGKB0zIYmdljAC6bYPkzqOlHQggxY+a0Nq0wWMaaB07mOCx5sdeNH7zE5/kMbfR9VPY95HtKH/IlMVj1TPN8iIZ+xqGmZ9PXDjp7/KwlzuYSL8CPrfCa0PDa0F85PATnhxv1jHJXrEvy4zHFo9gf+hMFjSk/FGrFljTXUBrm3ArHK3s5AeOhzinmETXUdhVb9CsK2zfkzMQ8pfiJjdvHuq325v2Ooj7T0M+o4fQFDS/2JwqfqYacoAkks/g5jJ+batEviP2Oyh1B0wn1cLGGsu1FX/PZtcJ/bvve6v/Ouvf7c9j95GFn/+EH3+/sZbudiZ69J/Pzjr34GCZD5SBCiCRQMBJCJIGCkRAiCea2n1GlAh49/rLZNewbMK8bXOnskd0hTyfMIRbrugbP9s/OxeN+/e79vrasuN/nLY2v9RrR+FJ/egpLvKYV65pyYd72ahiP+lBtJW8y1o61sNtKoRYrzGsWNaBqGM8dR1OGz2muwyHWtsU8n2J4g2NxHjGPxe3jClEjCsSe2Az7q4btiwMxMSrsMLYSCjlB8ZpHvaZB7xltTPQqd/qTmA/rTJCa5F0MmlNu1F8zlvz+Shf6fLzi814DWvVIEGNDvSNLfv/nfcnXulW7fF5T28jJfL22kvoZCSESR8FICJEECkZCiCSYW80ol4MtO6nLMGhGxT1+/qXi3pBj0uUn5ap0+2fTsdV+vHgkzHMWaniqHX77I6/xgsqSF/3z7YEre509uNGZqHSGObP6vP/5CZoWFELuUdR4SktiXosfz8V5z4ImFPOS2kJLppgXNO7L+Ro1n0CcZx2hJ1VbmOcMIe+nOh5Ep6jZtJqavrmEBKvG2rCgsQVNpzDgL0DUfBoIb6/SHkU+bzL2mwKQH2leMBi3iT2zo07FcrgGx/wHr3jc5xlFcvv832F5g9dyR8/1tW0dT+7z24/4fkjlNSfzliw3+fnUnZEQIgkUjIQQSaBgJIRIgjnOM6qCgyefV23MP1si5/UDdvm8nvKyoBl1ePejRtQW8i2qYf2xVWF/oc4ozoHVMEfYKt/vqK3TH298TXjWfzIUeqFxrvtckA9irVp7f9hBLKUabzEe84IKzcdjT+mGfkFhHrOY19PQHyisn+v057AS50GLx4/uxP2FeeArlbC/YyEHJs5zFudBGw4aV9BrqlEjCkQ9J2pMwATnLFDNh5rKmOs04s9h7lBIHgv9j/o3+z5cfW/0+9/03/yHMn/MN906/hpfo1kt+rylrsd9a7NmOlE9ujMSQiSBgpEQIgkUjIQQSUCbIO9h1g5GHgTwPIBVAA61WH0+kX8zI3X/gPR9XKj+bTSz1RMNzGkwevmg5DYz2zznB54i8m9mpO4fkL6Pi9E/PaYJIZJAwUgIkQTzFYxum6fjThX5NzNS9w9I38dF59+8aEZCCBHRY5oQIgkUjIQQSTCnwYjkNSSfJLmb5M1zeezJIPkFkn0kn6hbtoLkfSR3Zb+Xz6N/55C8n+QOkttJfjglH0l2kHyQ5KOZf5/Mlp9Lcmt2rb9MsthqX7PsZ47kwyTvTc0/ks+RfJzkIyS3ZcuSuL6ZL70k7yL5E5I7SV45G/7NWTDKpsf+KwA/C+AiANeTvGiujt+E2wFcE5bdDGCLmV0AYEtmzxdlAB81s4sAXAHgQ9l5S8XHMQBXmdklAC4FcA3JKwD8KYBPmdn5AI4CuHGe/DvBhwHsrLNT8+8dZnZpXe5OKtcXAD4D4FtmdiGAS1A7j6ffPzObkx8AVwL4dp39cQAfn6vjt/BtE4An6uwnAazLXq8D8OR8+1jn290Ark7RRwBdAH4M4E2oZefmJ7r28+DXhuwP5ioA96JW+5+Sf88BWBWWJXF9ASwD8CyyL7tm07+5fExbD+DFOntPtixF1prZiV6a+wGsnU9nTkByE4DLAGxFQj5mj0CPAOgDcB+ApwH0m9mJnirzfa0/DeBjONnAYyXS8s8AfIfkQyRvypalcn3PBXAQwN9kj7l/TbJ7NvyTgN0Cq4X+ec9/INkD4KsAPmJmrmHNfPtoZhUzuxS1O5DLAVw4X75ESL4HQJ+ZPTTfvjThrWb2etQkjA+RfFv94Dxf3zyA1wP4rJldBmAI4ZHsdPk3l8FoL4Bz6uwN2bIUOUByHQBkv/vm0xmSBdQC0RfN7GvZ4qR8BAAz6wdwP2qPPb0kT3Szm89r/RYA7yX5HIA7UXtU+wzS8Q9mtjf73Qfg66gF9FSu7x4Ae8xsa2bfhVpwOu3+zWUw+hGAC7JvMYoA3g/gnjk8/nS4B8AN2esbUNNp5gWSBPB5ADvN7Na6oSR8JLmaZG/2uhM1PWsnakHpffPtn5l93Mw2mNkm1D5z/2RmH0jFP5LdJJeceA3gZwA8gUSur5ntB/AiyVdni94JYAdmw785FsPeDeAp1DSF350PQW4Cn74EYB+AEmr/BW5ETVPYAmAXgO8CWDGP/r0VtVvgxwA8kv28OxUfAVwM4OHMvycAfCJbfh6ABwHsBvD3ANoTuNZvB3BvSv5lfjya/Ww/8XeRyvXNfLkUwLbsGn8DwPLZ8E/lIEKIJJCALYRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBITks1YMUJykOR+krdn3Sbr19lM8l6SR0n2ZzOY/En9TBEk15H8PMl9JAeyGSY+mfXuicfcRNKyYw5mPtwc1qn368TP2bN3JsRcoWAkmvHzZtaDWj+by1CbRAEAQPLNAL4H4PsALjSzXtRmWSmjNoMESK4A8AMAnQCuNLMlqDVf6wXwyibH7c2O+z4Av0/y6on8qvt5aYbvUyRAvvUqYrFjZvtJfhu1oHSCPwPwN2b2H+vWewHAH9St81sABgB80Myq2TovojZt0FSOu43k9uy4983kPYj00Z2RaAnJDag1i9+d2d2o9bn+aotN3wXgaycC0Skc9woArz1xXLGwUTASzfgGyQHUppjqw8m7nuWofXb2n1iR5J9lutEQyd/LFq9EraXvdDlEcgS1R7z/glqr0+hXf/YTx8QZioKRaMZ1mc7zdtSmH1qVLT+K2hxk606saGYfy3Sjr+Pk4//h+nWmwSoAPQA+mh27MIFfvdnPdaewf5EgCkaiJWb2z6hNA/4XmT2E2kSSv9Bi0+8C+Nckp/05s9pcbLcCGAXwG9PdXpx5KBiJqfJpAFeTvCSzPwbg10jeTHIN8LK2dG7dNrcCWArgDpIbs3XWk7yV5MVTPO4tAD5GsuN0vAmRLgpGYkqY2UEAfwvgE5n9AGoTIr4NwFMk+wF8C7Wv+/9zts4RAG9GbRqorZn+tAXAMUxdlP4H1B4L//1peisiUTRVkRAiCXRnJIRIAgUjIUQSKBgJIZJgRsGI5DUknyS5OxY0CiHEdDhlAZtkDsBTqBU+7gHwIwDXm9mOybYp5ruso733lI4nhDjzGR3rx3h5mBONzaRQ9nIAu83sGQAgeSeAawFMGow62ntxxUU3zeCQQogzmR/uuG3SsZk8pq1HrWbpBHuyZQ6SN5HcRnJbqTw8g8MJIRYysy5gm9ltZrbZzDYX8l2zfTghxBnKTILRXgDn1NkbsmVCCDFtZhKMfgTgApLnkiwCeD+Ae06PW0KIxcYpC9hmVib5mwC+DSAH4Atmtv20eSaEWFTMqO2smX0TwDdPky9CiEWMMrCFEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBJCJIGCkRAiCRSMhBBJoGAkhEgCBSMhRBIoGAkhkkDBSAiRBApGQogkUDASQiSBgpEQIgkUjIQQSaBgJIRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBJCJEHLYETyCyT7SD5Rt2wFyftI7sp+L59dN4UQC52p3BndDuCasOxmAFvM7AIAWzJbCCFOmZbByMz+BcCRsPhaAHdkr+8AcN3pdUsIsdg4Vc1orZnty17vB7B2shVJ3kRyG8ltpfLwKR5OCLHQmbGAbWYGwJqM32Zmm81scyHfNdPDCSEWKKcajA6QXAcA2e++0+eSEGIxcqrB6B4AN2SvbwBw9+lxRwixWJnKV/tfAvADAK8muYfkjQBuAXA1yV0A3pXZQghxyuRbrWBm108y9M7T7IsQYhGjDGwhRBIoGAkhkkDBSAiRBApGQogkUDASQiSBgpEQIgkUjIQQSaBgJIRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBJCJIGCkRAiCRSMhBBJoGAkhEgCBSMhRBIoGAkhkkDBSAiRBApGQogkyM+3AyIxqvN8fP17XLTo0gshkkDBSAiRBC2DEclzSN5PcgfJ7SQ/nC1fQfI+kruy38tn310hxEJlKndGZQAfNbOLAFwB4EMkLwJwM4AtZnYBgC2ZLVKn2uKnBTSb0c9s+yfOXFoGIzPbZ2Y/zl4PANgJYD2AawHcka12B4DrZslHIcQiYFqaEclNAC4DsBXAWjPblw3tB7B2km1uIrmN5LZSeXgmvgohFjBTDkYkewB8FcBHzOx4/ZiZGYAJ78HN7DYz22xmmwv5rhk5K4RYuEwpz4hkAbVA9EUz+1q2+ADJdWa2j+Q6AH2z5aSYAS10lgYdp5WuE8db6Tjh3x3JFg75cYvrx+Pp++AFw1S+TSOAzwPYaWa31g3dA+CG7PUNAO4+/e4JIRYLU7kzeguAXwbwOMlHsmW/A+AWAF8heSOA5wH84qx4KIRYFLQMRmb2AIDJ7q3feXrdEUIsVlSbttCYrkZU9Ruw0nwclWC30pii5pMLykCbty3n12ccl4a0YNGlE0IkgYKRECIJFIyEEEkgzehMZ6YaUanix0vlMO5tjJecaVFTiscPmg+KBb99wX8EGWzLhf0FTcni/qUhnbHoUgkhkkDBSAiRBApGQogkkGZ0pnGaNSKOjvv1R0adaaPRHvN2zDuK/oS8IhaL3u7q9Bt0dng7akzhI0v696s8pDMXXRohRBIoGAkhkkDBSAiRBNKMznBa9SNiOWhGQSOy4RFvDww6uzrixxtq0Vr0J7JSWD9oUG3j3h8GDYo9oSFfW8wzCu60qmUTyaI7IyFEEigYCSGSQMFICJEE0oxSZ7pzhcV+RLG2LOYJRY1o2M/gwrz/iORWrfTbL1/a1B32D3j3Dh3xx4saUtB8Yp4S8qFYLRfsVv9elXeULLoUQogkUDASQiSBgpEQIgmkGZ3pxLyiWItW9rVosbasGmy0eQ0m94oNzu5/g584uP/8oNkElj3rNabeh3zeUPXZF7wd8prY3u7tDm8j9j+ykFekPKMzBt0ZCSGSQMFICJEECkZCiCSQZnSm06I2DVEzCrVgqPrxXO8yZw+9ZrWzX7rK7//ayx90doF+f3c9tNnZ+dFVzu450u/symGfhxT9jRpYwzxu+ZCnZC36HYlk0J2RECIJFIyEEEmgYCSESAJpRmcYDf2LItUWGlKL7dnhe1APr/YfkY2v3OfsPz9rq7ML9HlH289f5+y+1RudvaQz9MCOVIJG1Or9izMW3RkJIZJAwUgIkQQtgxHJDpIPknyU5HaSn8yWn0tyK8ndJL9MsthqX0IIMRlT0YzGAFxlZoMkCwAeIPmPAH4LwKfM7E6SnwNwI4DPzqKv4nTQqmd12fc/Kgz5PJ6XDvs8pH8Y9nYHS85+/shyZy8dCnk/Jb9+A22xybXyhBYqLe+MrMaJDlyF7McAXAXgrmz5HQCumw0HhRCLgylpRiRzJB8B0AfgPgBPA+g3sxP/RvcAWD/JtjeR3EZyW6k8PNEqQggxtWBkZhUzuxTABgCXA7hwqgcws9vMbLOZbS7ku1pvIIRYlEwrz8jM+kneD+BKAL0k89nd0QYAe2fDQeGJtVWMGkqLntEM/X+iBlMNPbGX7vI9rIcf6HX2f+j/pbA/by7d6Y+37KnjzrZj3o7EHtwN76+tef8i1aKdOUzl27TVJHuz150ArgawE8D9AN6XrXYDgLtnyUchxCJgKndG6wDcQTKHWvD6ipndS3IHgDtJ/jGAhwF8fhb9FEIscFoGIzN7DMBlEyx/BjX9SAghZoxq0850oiYSNZM4z1jsKV30uao25ntic9fzzl435vOCVu5Y4rcPx2/f3+/HX/S1bXHeNBZC7mzwt3GeNBURLBR0JYUQSaBgJIRIAgUjIUQSSDNKnfjvIs4VH7BcyEMKeUVx3rG2nm6/+9APqTo05A/w5DPOLDxbaOpPZTzUnoWe2zGPKPrDLt9fyeI8ablpzpOmf7/JoksjhEgCBSMhRBIoGAkhkkCa0RlOq1o1i3k57T6Ph0Ejiv+dbMTnAVVHfR5SNeQlNUC/x7bQY5vdvniaXd624G/UjFSLtnDQnZEQIgkUjIQQSaBgJIRIAmlGZxqt8o4a+huF8ZinEzWn2P8o1K6xK+QNxbnuI6H/EAshLylqQsUwHjUi5RUtWHSphBBJoGAkhEgCBSMhRBLQ5nDucpIHATwPYBWAQ3N24Okj/2ZG6v4B6fu4UP3baGarJxqY02D08kHJbWa2ec4PPEXk38xI3T8gfR8Xo396TBNCJIGCkRAiCeYrGN02T8edKvJvZqTuH5C+j4vOv3nRjIQQIqLHNCFEEigYCSGSYE6DEclrSD5JcjfJm+fy2JNB8gsk+0g+UbdsBcn7SO7Kfi+fR//OIXk/yR0kt5P8cEo+kuwg+SDJRzP/PpktP5fk1uxaf5lksdW+ZtnPHMmHSd6bmn8knyP5OMlHSG7LliVxfTNfekneRfInJHeSvHI2/JuzYJRNj/1XAH4WwEUArid50Vwdvwm3A7gmLLsZwBYzuwDAlsyeL8oAPmpmFwG4AsCHsvOWio9jAK4ys0sAXArgGpJXAPhTAJ8ys/MBHAVw4zz5d4IPA9hZZ6fm3zvM7NK63J1Uri8AfAbAt8zsQgCXoHYeT79/ZjYnPwCuBPDtOvvjAD4+V8dv4dsmAE/U2U8CWJe9Xgfgyfn2sc63uwFcnaKPALoA/BjAm1DLzs1PdO3nwa8N2R/MVQDuBcDE/HsOwKqwLInrC2AZgGeRfdk1m/7N5WPaegAv1tl7smUpstbMTszDvB/A2vl05gQkNwG4DMBWJORj9gj0CIA+APcBeBpAv5mVs1Xm+1p/GsDHcLLhykqk5Z8B+A7Jh0jelC1L5fqeC+AggL/JHnP/mmT3bPgnAbsFVgv9857/QLIHwFcBfMTMjtePzbePZlYxs0tRuwO5HMCF8+VLhOR7APSZ2UPz7UsT3mpmr0dNwvgQybfVD87z9c0DeD2Az5rZZQCGEB7JTpd/cxmM9gI4p87ekC1LkQMk1wFA9rtvPp0hWUAtEH3RzL6WLU7KRwAws34A96P22NNL8kRntPm81m8B8F6SzwG4E7VHtc8gHf9gZnuz330Avo5aQE/l+u4BsMfMtmb2XagFp9Pu31wGox8BuCD7FqMI4P0A7pnD40+HewDckL2+ATWdZl5gbbqPzwPYaWa31g0l4SPJ1SR7s9edqOlZO1ELSu+bb//M7ONmtsHMNqH2mfsnM/tAKv6R7Ca55MRrAD8D4Akkcn3NbD+AF0m+Olv0TgA7MBv+zbEY9m4AT6GmKfzufAhyE/j0JQD7AJRQ+y9wI2qawhYAuwB8F8CKefTvrajdAj8G4JHs592p+AjgYgAPZ/49AeAT2fLzADwIYDeAvwfQnsC1fjuAe1PyL/Pj0exn+4m/i1Sub+bLpQC2Zdf4GwCWz4Z/KgcRQiSBBGwhRBIoGAkhkkDBSAiRBApGQogkUDASQiSBgpEQIgkUjIQQSaBgJIRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYLRAiCbXWKE5CDJ/SRvzzpD1q+zmeS9JI+S7M9mG/mT+lkdSK4j+XmS+0gOZLNBfDLrsxOPuYmkZccczHxIYsYXcWaiYLRw+Hkz60Gt98xlqE14AAAg+WYA3wPwfQAXmlkvajOilFGb7QEkVwD4AYBOAFea2RLUGqX1Anhlk+P2Zsd9H4DfJ3n16XxTYvGgYLTAsFpnvm+jFpRO8GcA/sbM/qOZHcjWe8HM/sDMvpet81sABgB80Myey9Z50cw+bGaPTeG421BrDvbycUn+WjbP1lGS3ya5MVtOkp9ibb6649mcYa/Nxm4n+blsLq4Bkv98Yrts/M0kf0TyWPb7zXVj3yP5RyS/n237HZKrsrEOkv8fycPZneGPSK7NxpbV3RHuJfnH2dRaYg5RMFpgkNyAWmP33ZndjVpP6q+22PRdAL5mZtUW60123CsAvLbuuNcC+B0AvwBgNYD/jlpXTaDWWvVtAF6F2lQ4vwjgcN3uPgDgjwCsQq2z5Rezfa4A8A8A/hNqnQZvBfAPJFfWbftLAP4tgDUAigD+j2z5Ddmxzsm2/XUAI9nY7ajdJZ6P2l3lzwD4d6dyHsQMmM82oPo5bW1BnwMwiNqdjaHWDrQ3G9uQLbuwbv0/A9CP2kwPv5ct2wXg16dxzE3ZfvtR+6M2AH8BvNw99B8B3Fi3fhuAYQAbUWuK/xRqk1K2hf3eDuDOOrsHQAW1IPLLAB4M6/8AwK9mr7934v1k9m+gNvkgAPwagP8B4OKw/VrUJqLsrFt2PYD75/u6LrYf3RktHK6zms7zdtSmClqVLT+K2nxh606saGYfs5pu9HXUpqIBancmL68zDVahFjA+mh27kC3fCOAz2SNRP4AjqE2euN7M/gnAX6I2w3AfydtILq3b58vz65nZYLbt2dnP8+H4z8PPeba/7vVw5hsA/FfUHl/vJPkSyT/LZl3ZmPm8r87X/we1OysxhygYLTDM7J9Ru7v4i8weQm3Sx19osel3AfxrktP+TFht3rRbAYyidjcC1ALK/2JmvXU/nWb2P7Jt/pOZvQG1qc5fBeD/rNvly1NaZd8KrgDwUvazEZ5XYArTDJlZycw+abVpwt8M4D0AfiXzcwy1GV1P+LnUzH5quudBzAwFo4XJpwFcTfKSzP4YgF8jeTPJNcDL2tK5ddvcCmApgDvqhOb1JG8lefEUj3sLgI+R7ADwOQAfJ/lT2b6Wkfw32es3knxTdmcyhFoQq9eq3k3yrdmUVn8E4Idm9iKAbwJ4FclfIpkn+T+jFszubeUYyXeQfF0mTB9HbTaYqtVmRf0OgP+b5FKSbSRfSfKnp/iexWlCwWgBYmYHAfwtgE9k9gOo6TRvA/BU9ijyLdQ0lv+crXMEtTuGEoCtJAdQ056OIROlp8A/oPZY+O/N7OsA/hS1x6LjqE1j9LPZeksB/L/Zus+j9oj453X7+TsAf4Da49kbAHww8/Ewanc0H822+RiA95jZoSn4dhZqExAeR21et39G7dENqN0hFVGbD+xott6pPLKKGaCpikRSkLwdtRlMf2++fRFzi+6MhBBJoGAkhEgCPaYJIZJgRndGJK8h+STJ3SqSFELMhFO+M8q+In0KtWLKPQB+BOB6M9sx2TbFQrd1tPee0vGEEGc+o2P9GC8NcaKx/EQLp8jlAHab2TMAQPJOANei9vXohHS09+JNF//6DA4phDiT2frY5yYdm8lj2nrUpe2jdne0Pq5E8iaS20huK5WGZnA4IcRCZta/TTOz28xss5ltLhQaenQJIQSAmQWjvairIUKtOrxljZAQQkzETILRjwBcQPLcrIbo/QDuOT1uCSEWG6csYJtZmeRvotaWIQfgC2a2/bR5JoRYVMzk2zSY2TdRq6QWQogZoXIQIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBJCJIGCkRAiCRSMhBBJoGAkhEgCBSMhRBIoGAkhkkDBSAiRBApGQogkUDASQiSBgpEQIgkUjIQQSaBgJIRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEnQMhiR/ALJPpJP1C1bQfI+kruy38tn100hxEJnKndGtwO4Jiy7GcAWM7sAwJbMFkKIU6ZlMDKzfwFwJCy+FsAd2es7AFx3et0SQiw2TlUzWmtm+7LX+wGsnWxFkjeR3EZyW6k0dIqHE0IsdGYsYJuZAbAm47eZ2WYz21wodM/0cEKIBcqpBqMDJNcBQPa77/S5JIRYjJxqMLoHwA3Z6xsA3H163BFCLFam8tX+lwD8AMCrSe4heSOAWwBcTXIXgHdlthBCnDL5ViuY2fWTDL3zNPsihFjEKANbCJEECkZCiCRQMBJCJIGCkRAiCRSMhBBJoGAkhEgCBSMhRBIoGAkhkkDBSAiRBApGQogkUDASQiSBgpEQIgkUjIQQSaBgJIRIAgUjIUQSKBgJIZJAwUgIkQQKRkKIJFAwEkIkgYKRECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCRQMBJCJIGCkRAiCRSMhBBJkJ9vB0RiVK35eBvnxg+x6NCdkRAiCVoGI5LnkLyf5A6S20l+OFu+guR9JHdlv5fPvrtCiIXKVO6MygA+amYXAbgCwIdIXgTgZgBbzOwCAFsyWwghTomWmpGZ7QOwL3s9QHIngPUArgXw9my1OwB8D8Bvz4qXIhlY8ZqSzbaEJI1q0TAtzYjkJgCXAdgKYG0WqABgP4C1p9c1IcRiYsrBiGQPgK8C+IiZHa8fMzMDMOHXMCRvIrmN5LZSaWhGzgohFi5TCkYkC6gFoi+a2deyxQdIrsvG1wHom2hbM7vNzDab2eZCoft0+CyEWIC01IxIEsDnAew0s1vrhu4BcAOAW7Lfd8+Kh2JmtMobCjCubtbUZrXFDnNB86G3W2pO0X9pSAuWqSQ9vgXALwN4nOQj2bLfQS0IfYXkjQCeB/CLs+KhEGJRMJVv0x4AMNm/o3eeXneEEIsVZWALIZJAtWkLjRYa0XQ1IQRNiFW/gOPl5v60+f931Q7/kSOnqSmpdm7BojsjIUQSKBgJIZJAwUgIkQTSjM50ZqgRsRxEoagJxfFKGC8FzWhk1Nt5/xFrq3Z6d9oL3g55SQ0KUNCgWmpK0pDOGHRnJIRIAgUjIUQSKBgJIZJAmtECY7oaEUsVv/54KYwHTagc1g/7rx4f8ONB42mLeUz0xdNR4bF80Ii8xARKQ1ow6M5ICJEECkZCiCRQMBJCJIE0ozONVrVZUZMJPasbNKKxcT8+6m0L41FzibVlLPiPlI2OBTvkIQWbHR1+vCvYsXatEPOSptkvSSSD7oyEEEmgYCSESAIFIyFEEkgzOsNpmVdUaZFHFDWikh9vyMtZtsSZ1a72pv61HR/2+z90xNvlkMcUa99y4f9lPue3b4vjLUQi5R0li+6MhBBJoGAkhEgCBSMhRBJIM0qdaeYVsUVeEYNmFPOIokZjK3udPXjBMm+v8xpOudNrMN37ljp7+WNBY3rhpeCPz0vCaChGK3qbhaAhmfKOzlR0ZySESAIFIyFEEigYCSGSQJrRQiP0sG7oPxQ0I8Q8n64eZw5v9JrPwUv8R6bymkFnv/bsfc5+eOcmZ+dKvc5eeuio39/Bw3790EObnaFWrRx6aIf+R7GWTaSL7oyEEEmgYCSESAIFIyFEEkgzOsNoVYvWsH6o9bKoKYXaLAt5PGPLfB7P6DqvMf32xVuc/eu9e539wfzbnb398YucvSxqQAELtXXx/YiFg+6MhBBJoGAkhEiClsGIZAfJB0k+SnI7yU9my88luZXkbpJfJlmcfXeFEAuVqWhGYwCuMrNBkgUAD5D8RwC/BeBTZnYnyc8BuBHAZ2fRV4HG2qpWWTTWql9PNc6r5jWawkjoLzTu93d2wecJfX/Ur7/94Fl+fwNB4wrzsjH6G/sVqf/QgqXlnZHVOJHZVsh+DMBVAO7Klt8B4LrZcFAIsTiYkmZEMkfyEQB9AO4D8DSAfjM78W9tD4D1k2x7E8ltJLeVSkOnwWUhxEJkSsHIzCpmdimADQAuB3DhVA9gZreZ2WYz21wodLfeQAixKJlWnpGZ9ZO8H8CVAHpJ5rO7ow0A9jbfWpwSUSOJ/Y1i7VW0g+bCQqjlCrVpNjzi7M6XvL3sJ74H9oc7PhCO780lP/HHW/2UvzuuHh/wxw/vry3UpsX3Yy1qz9S/6MxhKt+mrSbZm73uBHA1gJ0A7gfwvmy1GwDcPUs+CiEWAVO5M1oH4A6SOdSC11fM7F6SOwDcSfKPATwM4POz6KcQYoHTMhiZ2WMALptg+TOo6UdCCDFjVJt2phPnno8aUZhnDIVwycfCeOh3lH+hz9lnjfrx5bt9/6Oo0XTs93lIfGG/s6sjXpNq6/A9shlsi/7nWmhm4oxB5SBCiCRQMBJCJIGCkRAiCaQZLTTCv5eGntBBc2GX7ydkw6POLh/wmhEP+R7V7U+HedAC1RG/P5ivXWtrDxpRV6dfP/Q7iv2Wpt3zWrVtyaI7IyFEEigYCSGSQMFICJEE0ozOcBr6G8XarTBVfYOmkgt5SWH1tk6v4VRHx/z+h4eDQy16ckeNqKe7qW0doWdfC41ItWhnLrozEkIkgYKRECIJFIyEEEkgzehMo0V/o5YaUsMVjwt8Xk/byhV+f+PjfvXQD6mh31LUpEI/JcTas85Yi+Zr55RXtHDRnZEQIgkUjIQQSaBgJIRIAlqLvJDTejDyIIDnAawCcGjODjx95N/MSN0/IH0fF6p/G81s9UQDcxqMXj4ouc3MNs/5gaeI/JsZqfsHpO/jYvRPj2lCiCRQMBJCJMF8BaPb5um4U0X+zYzU/QPS93HR+TcvmpEQQkT0mCaESAIFIyFEEsxpMCJ5DcknSe4mefNcHnsySH6BZB/JJ+qWrSB5H8ld2e/l8+jfOSTvJ7mD5HaSH07JR5IdJB8k+Wjm3yez5eeS3Jpd6y+TLLba1yz7mSP5MMl7U/OP5HMkHyf5CMlt2bIkrm/mSy/Ju0j+hOROklfOhn9zFoyy6bH/CsDPArgIwPUkL5qr4zfhdgDXhGU3A9hiZhcA2JLZ80UZwEfN7CIAVwD4UHbeUvFxDMBVZnYJgEsBXEPyCgB/CuBTZnY+gKMAbpwn/07wYQA76+zU/HuHmV1al7uTyvUFgM8A+JaZXQjgEtTO4+n3z8zm5AfAlQC+XWd/HMDH5+r4LXzbBOCJOvtJAOuy1+sAPDnfPtb5djeAq1P0EUAXgB8DeBNq2bn5ia79PPi1IfuDuQrAvag1tEzJv+cArArLkri+AJYBeBbZl12z6d9cPqatB/Binb0nW5Yia81sX/Z6P4C18+nMCUhuAnAZgK1IyMfsEegRAH0A7gPwNIB+MzvRX2S+r/WnAXwMwIl5klYiLf8MwHdIPkTypmxZKtf3XAAHAfxN9pj71yS7Z8M/CdgtsFron/f8B5I9AL4K4CNmdrx+bL59NLOKmV2K2h3I5QAunC9fIiTfA6DPzB6ab1+a8FYzez1qEsaHSL6tfnCer28ewOsBfNbMLgMwhPBIdrr8m8tgtBfAOXX2hmxZihwguQ4Ast99LdafVUgWUAtEXzSzr2WLk/IRAMysH8D9qD329JI80bltPq/1WwC8l+RzAO5E7VHtM0jHP5jZ3ux3H4CvoxbQU7m+ewDsMbOtmX0XasHptPs3l8HoRwAuyL7FKAJ4P4B75vD40+EeADdkr29ATaeZF0gSwOcB7DSzW+uGkvCR5GqSvdnrTtT0rJ2oBaX3zbd/ZvZxM9tgZptQ+8z9k5l9IBX/SHaTXHLiNYCfAfAEErm+ZrYfwIskX50teieAHZgN/+ZYDHs3gKdQ0xR+dz4EuQl8+hKAfQBKqP0XuBE1TWELgF0AvgtgxTz691bUboEfA/BI9vPuVHwEcDGAhzP/ngDwiWz5eQAeBLAbwN8DaE/gWr8dwL0p+Zf58Wj2s/3E30Uq1zfz5VIA27Jr/A0Ay2fDP5WDCCGSQAK2ECIJFIyEEEmgYCSESAIFIyFEEigYCSGSQMFICJEECkZCiCT4/wHQAwkva/EFYwAAAABJRU5ErkJggg==)
%% Cell type:markdown id: tags:
# Generate simulated RGCs responses to small dataset (lurz - training examples)
%% Cell type:markdown id: tags:
##### RGCs response generation to smaller dataset:
%% Cell type:code id: tags:
``` python
#Get indexes of 2000 images to remove from the training data of Lurz dataset
tiers = np.load('D://inception_loop/RGC_sim_data/data/static27012021_original_images/meta/trials/tiers.npy')
random.seed(1)
n_remove = 2000
train_ids = [i for i, tiers in enumerate(tiers) if tiers == "train"]
remove_ids = random.sample(train_ids, k=n_remove)
tiers_new = np.delete(tiers, remove_ids)
```
%% Cell type:code id: tags:
``` python
#Retrieve image sets from evaluation data set of lurz2020 #5993 and remove 2000 training examples
paths = 'D://inception_loop/RGC_sim_data/data/static27012021_original_images'
paths_new = 'D://inception_loop/RGC_sim_data/data/static23032021_small_set_original'
images = []
i = 0
for n in range(5993):
if n not in remove_ids:
x = np.load(paths+'/data/images/'+str(n)+'.npy')
np.save(paths_new+'/data/images/'+str(n)+'.npy', [x])
np.save(paths_new+'/data/images/'+str(i)+'.npy', [x])
images.append(x-np.mean(x))
i+=1
images = np.vstack(images)
```
%% Cell type:code id: tags:
``` python
plt.imshow(np.load(paths_new+'/data/images/'+str(n)+'.npy')[0][0])
```
%%%% Output: execute_result
<matplotlib.image.AxesImage at 0x207ddcd3670>
%%%% Output: display_data
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAXAAAADdCAYAAACi7Cm7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA0v0lEQVR4nO2deZBdZ3nmn/fuSy+3d7XUkmXJkhfJlgzCxqwG45RNQiDLJCEzFFPjiklVUgU1VAhkaiZLJVOZVAKTmsmQMYHgpAiQYXUYAniME3aDvEi2JVmytbdaavVye737/eaPvk7Uep6LWlvbh3p/VSpJr84533K+8/XRed7FQghwHMdxokfspe6A4ziOc2n4Bu44jhNRfAN3HMeJKL6BO47jRBTfwB3HcSKKb+CO4zgR5bI2cDO7x8yeM7PnzeyDV6pTjuM4zoWxS/UDN7M4gIMA7gZwEsCPALwzhLCv3TmpRC5kU4XlRtF+iOufK1ZvsFHYQiYlzm3qTsWMzzdhi7MNwrR0sLCp0xt8oDUv4n7U62yLqbkT40noOW6m2B5fqPH5ybhoW0+InPumsIm+N7KiHQBNYQ5iSIlFMcfiem1nXRys1kIzycfFqvqqMbGMmwm+ploLjXS7ORY2cX58oUq2kBadBxDa3E9qRwyznuFzk/Ni4IDeA8T6UmOP1drMcUW0JdppphO6T4JmktuPl3kd1zp5IcZ42pe6pJa3GFJp4uRECGHgfPvKe8/cBuD5EMJhADCzzwB4O4C2G3g2VcCrt963zGYVXnmNnpw8Pz4xx8bJIpnqN24gW2JiXl4z5NJsE4un2s0/FNptWBAPjnrok7O8McbnKvqaYnOKjU9zOx1i7sTGWO/vkM0sDmfI1v34GJ+/pkC2RkYvp+TEItmsIjaSPLc9c0OXvGa5h+ezITaNvqd5PmPiB6f6oQ3oh7bazeNcWMNz3HVM7KoAUrNsL/XzJppc4M2huFlvttlJPja5yLbOx46Trbp5jbxmvYPbUj8UYlVuZ/p6fq4Gv8vrdemafH51kNfnzCa+Zu6s/qGQPzzDxio/b6VNvaJDei0sDvB9LzzH+8rJuzrJ1nVUv0BWu8QPbjGkJz/2/mPq/Mv5hLIOwIlz/n6yZXMcx3FWgct5A18RZnY/gPsBIJPUb1OO4zjOxXM5b+CjANaf8/eRlm0ZIYQHQgi7Qgi7Uon8ZTTnOI7jnMvlvIH/CMAWM7sWSxv3rwD41R93Qogb6l3Lv3MmnzlCx8Wy+kuMFDFzWT6/Ko5r940zxd+7lYiaWOTvlvF5/b26meFvh7FF/v4mf3y2+ZFqo+NkCwm+faFDzMcEfw9MHi/LdrpeEIJlN3/TS4zP8nFruuU1mzmej0RR6BlCZMpM8HgAoJZlTUIJWkqUbYql0BDHAUC5wPbuI3zfCz+c4JOloAw0c/ytP77Ac6QEtv5n9DffxSGej1qW2y/dPMJti7UNAAtruP1qJ09edoK/7yZY9sCRfyO+NwPoOcDnF/bxmu3by7pJrCSeKwDNrPh+r7StLvX8670iPctzXyvwd3kTn7tnNuu1kBTSXNfxNmKv4JI38BBC3cx+E8DXAcQBfCKE8OylXs9xHMe5OC7rG3gI4asAvnqF+uI4juNcBB6J6TiOE1F8A3ccx4koV92NcBlNIH5+hNQQBRcBBw7L08MGFjfrwz1kS54ukq06ogWU5NPcVtgwTDYZTdjQzvnxo6f4/C4WAhvD3KdGWkceht6NZEsWtRBJqIhNIf4CgDWEgLLI7YQy25JtonprIugndHOgRmVIBxcpUgvcVrrI40zOsvClgrLyZxZkOyHGnlOLa/j85DQfVxnQAWnlPhbYUnM873Pr+fHsfkGH9OVP8v0oD7DAVu4V12wz9vQMn1/PiGhZESEZhA6YPSubkQJfdpLnM/PDQ2QrvXqrvGZmjNXBWi+v+cwkr5m5DSJoD0BCPG5qfaVm+f6W+7QwWhNLfvYavQco/A3ccRwnovgG7jiOE1F8A3ccx4kovoE7juNEFN/AHcdxIsqqeqFYvY746fNSSiZFF0bYCwQAgghRT8yUREOs+Bav4/BlAEgN30i23CnhdZEUIdltQm5rN/WRLS5SbqamWcFOnuEQdQCoiTD1pkjfGpsXXgrK40SkcwWA0CMSjtWEF0uerxmUBwva5B4XaRFSE+wNUVmjPVMSJZ7PyRvZe6BjjOeoluP7Vs+38f4RXa9l+fzSWu5nerpdamBxfp8Im98j5qOXPUMAIFZnrxyV4jYuuqTSNABA1xw/W+kN7DllwvuonuH5nNyuvTvUHE9u42P7sIVsmZMiJQMgX02TZzm+P4hnKFHWKXs7DhbJVlrPz0vfM2L/EPneAWD09bwvqdzu7fA3cMdxnIjiG7jjOE5E8Q3ccRwnovgG7jiOE1FWVcQMqQSqG5eHzsfL/MVeheECQLWHBZysEjFE7bv+J0WNPAC1HhYRklNC7BC5nStCuAKA5ByPqdbJU10e5PE01mmxNTvBY0ocOsP9FCHuYZ2oeZjQol2tl8O/k6d57mxeJHxuk3M9vsiCqTq/2cFpEWY3aEGpkRIFZtsUESbEYep6gBbYVKHjSoHn05patMuMcph3sijygWd5zag6qkuNcaf69/Acy9qwKtUCgGYHr0VVLDgxxWKrumvxLYOynY6TLEhXunk801t5Poem9fqIzbIA2+xi4b3Sx2OstykcXbxFpL4Q66Z3r9g/2hQRX/dttl1MYXN/A3ccx4kovoE7juNEFN/AHcdxIsplfQM3s6MA5gA0ANRDCLuuRKccx3GcC3MlRMw3hRBERVdBAGKV5YKJFBteOCpPz6RFAdGeAh8oIgJV8WIASB/g3N3NQRbT6p1CQN17Ql6zev1assWnWMhLTPLYa/06h3RqTBQRFsWGLSnEMCGgWBvBL66KxKrivCKCNojozHYsbuc5qnXwPUrN634qoWd6i4iWFRFwSpgsJdsUsi1yOyYCTss9fFGVNxwAbCuLYT0HRVTuHN+L4la9PvKnWYhMlNjW7Obz52++QV6z0s1jqnUI0e6AyG8uolAHdwsBFYA1xD1uCmGzn9dXeVCL/hhge6Wb19f8CI9R3XMAiIm1uDggnBtEP2tdeqtNT4l7vEmM6Z/k6f4JxXEcJ6pc7gYeAHzDzB43s/uvRIccx3GclXG5n1BeF0IYNbNBAA+b2YEQwrfOPaC1sd8PAJkUJ2RyHMdxLo3LegMPIYy2fh8H8EUAt4ljHggh7Aoh7Eomuc6d4ziOc2lc8hu4meUBxEIIc60//xSAP/hx54S4oVZYLgY2ckIAGWdxDgAaW0bYuMgiQExE+dW62hQqzbA4GZvgyMNkiYWJxrp+ec3prStLETl4Yppsqf1T8ppKnKxu5si2xIwYpyo2nNIRbCqir9HLP3jjkyxImSh+DAD1QU6vmyizEqhSbtaUgAogNSvS0c7y2GN3T/LJX2URceF2kZYYQGO/iEwVGX/jFZ7jhK4VjGoXj9PELYrPshDYeVw/slrU5XWYKPFxqigxoAvxBuELUBOpeIPxc9UupWrmND+vi+t5zeUP88Q3c/q5rvRx+93P8w3pPshiaTOl57jay22ZmJCZTXxcQi8vlEQ/VaRwOy7nE8oQgC/a0gOfAPB3IYSvXcb1HMdxnIvgkjfwEMJhADuuYF8cx3Gci8DdCB3HcSKKb+CO4zgRxTdwx3GciLKq+cBjpSqyT59cZgsdIjR4gL0EACAxIXJ/Cw+LRh97sWT3j8lr1tcLT5If7CVT7BYON5YhwACyk6xs50+IHMEzrKpbRocGN3t5TMkDo3x+QhTxvWaAbO1SC6TG2bskKFVeFaMu60LJ6WPCs6bEHivxEe4n+vV8JBbYrafwPHs5zN7N53b9PK+FRE175Uxu5WvmH12Zd4fyNgGA9Awfq/Ke91TY82lug/a66NnH961vL3toldbyNXPPidzuACZ28ppbHF5ZaoLEInsJzW3QBZlDjD1OVAFhlat+aptO31Du5X4O/4D7VOvgdRyv8PMLAJkz7EqibA1RKLnWqdfX4hAfWxpsUxBB4G/gjuM4EcU3cMdxnIjiG7jjOE5E8Q3ccRwnoqyqiIkAhPPy/FpNFAAe5nzcAJA8cpovKYr4xhZF3GpOix3xA8fY2CPaL7MgpPoOAPkTLBDGJ1mADWuEaDdZlNe0khAIVTHarBDYijwfNqNzM4ecEA2FiKmKPJvI4QwAzU6ee1MiqDg/M6HD81WB60SJRar4F1kQr/0yF4Pe2nNWtgNhP/aF6/k4oTvVGlqMquXZHhM5xpOnec10pLX4rIrz5p/jvqunYGqbTl3RMcrrqy5SX6hc6KkZthUO6twCZSFUVwe5YLjKid+/R6+5mc0seFZFTu7MWU5XoIqNA0BqksVnlftbpUUoDeittpbjtdB1RCyGNvgbuOM4TkTxDdxxHCei+AbuOI4TUXwDdxzHiSirK2KaUV7rMMdiWqLI+bgBIIgow+ZmzhFudZHjVxT2BYC4iORUfcIMC0rNNZznGgBiMyLqcopzfyMu8ij36yjU0sYC2XINIeAIYVUJlmFBC0qm8m+r3OFq3pQA2oaQ5ntpFZEjXOtJSImC0CqP8+QreI62Z/n+3Nv3tGxnXYLv22/13Eg2VVR46JsslgLAwo0iMlYUVW528nymT2vxuZnhe9ToZSFQkZ7VQmCqyMJ54QU+rpEWxY+FYFga1NGIMaHFN1M8H/EyC5Onb9djVNGuZ3dwn0YeZeeEeFlHWI++iSuKqQLIceFv0HVMi/HJs/wcKgG3Hf4G7jiOE1F8A3ccx4kovoE7juNEFN/AHcdxIsoFRUwz+wSAnwEwHkLY3rL1AvgsgI0AjgL4pRCCUOnOby2OZv95QkDoosNik6JqLAA0WOSSwldWiCVtogRDgoVE6+Y+NTuEQKeKBQOAuGZ513Vkyz7HIlco6rFnj4i2VPtinM1JTuca69cCrIruDKLQMeKiMK8odAwAoUtEqwmhOQgxLFZikQkA6t18zRd+iW07d7Dq9tMDLFgW4lrUTYpq1MlFnveFNUKUbeii1/UMj7PzEAv3IS0KWfdz6lUAyOzn1MKhJtLJvvJasqnUr4Au4pvfNy4O5Haq1w1xOzMc9QgA9W4RVVtk0a/eycdlpvUzuChSsuZP8bHlAR5jLasjaBNCh6yLSMpKgY/rOKnflVWq5rRIE92OlbyBfxLAPefZPgjgkRDCFgCPtP7uOI7jrCIX3MBDCN8CcP4r3NsBPNj684MA3nFlu+U4juNciEv9Bj4UQnixrMlpAPz/pRZmdr+Z7Taz3dW6rvzhOI7jXDyXLWKGEAKANh+DgRDCAyGEXSGEXamEKJ/mOI7jXBKXuoGfMbNhAGj9LpQNx3Ec52pyqaH0DwF4N4A/bv3+5RWdFQKstlzxtqLIk61C2aG9Q0x4YtisyH+9IHKEA6huHiRb6iAXvZU5xpV3BrT6n8qzgl7dyCHVqaM6L3WYZCef0MOhvVBFomdFMeg2NIc4lH/8Nm4n/wucm/36gvbkOFvmeRqd42tOTnDf139ee10kFthzIneK30fmbmHvoY0pnuNa0I/CliTfd5XDeXEN26pdugBxRuSVjs3y58XyJl4fyTldOLo5UODzhzkkOzPGz9bELp1/X6UHUB4nyvNJ5TKvDvPzCwApEU6unlfL8XzmxnVO/oZI1aDym6vc3bJKM4AQ44P7fsCeZOVreD5rnTqPuwVen81+Tg8CXZP9wm/gZvZpAN8HcL2ZnTSz+7C0cd9tZocAvKX1d8dxHGcVueAbeAjhnW3+6a4r3BfHcRznIvBITMdxnIjiG7jjOE5EWdV84I1MHPNblotXHd/kr/NBhMwDOqe2zYt8ujs5XDg5ufKh1tezeGQi93ZsXocGB1Gwt97BYdGpce57Y1yLmLFOFqSszIJWyHM4ef0VHMZ/+B060faG7Xw/XttziGxTVRYXj87rXOav6T9Mtl9b989k++kcxyrfd8Pr5DWf+YubyVYXFXtH8kWyFWJa0Fb83ew2siVKLGZ1iVQHi2v0+1FTZHqY27GGbLEaXzNzVAvSlQ0895kzLIzGplnEzJ/WRY2nt3BH80+v7J3P5rltq+s814sbWdxMzrOgPb+O12zxOt2fymZeS3PH+fy8iFoffFw7UWRFTYHaGhbjSwM8b/Gq9rSudvNeMTci9qqH5en+Bu44jhNVfAN3HMeJKL6BO47jRBTfwB3HcSLKqoqY8UoT+aPLBYL6dhYcleAHAEEEPmZPcP7s9IFTZGsKsRMAEiainkSEpYpG7D6qc1XX8hx1lT3NokpD5DeOF0R0JQBLcRRas4dFoWM/y2LW1rs5J/b9vSwsAkAtcN8/8e03ku3VOw+SbSirBbbuOIuGexavIdsHPs6hBb/4Kyx2AsDuXyiSrVHmdbOj8wTZvr24lWyvz/F4ACBjfI9VAWIloA48qQvZqiLAUzfx/TWh5Yc4Rw4DQC0volBv4yi/dd/k47LH9X2zJoubp+9dT7bOkxzhmB1lITA5pZPZpUbZGeDwu9aRrXwtH3fnjfq+xUR6ppMjBW7ncR5PotRGbBXRthu+Mkm2whz3c3H9ygsVtyuqrPA3cMdxnIjiG7jjOE5E8Q3ccRwnovgG7jiOE1FWVcREo4nYzHIx0UosPJUHdUSfEgJrvRyxFetgcTC2tyivqYrrqpSZPQdZEFoc0ulCF4b55+LCEPeznmdRJHPDZnnN+RE+tud1nNJ1e8dzZPuZ/j1kG0xo4Wq8zsJVfI7Hsy5bJNuZio7oO1rmAsqnSizWDuzhef+HO7bLa37gxm/w+XEWtG/PsO07ZU732YBODfy3x28nW9cRXocqUrc8qKNdGylR1Pg4K5Yzm1hQbiZ0PzsPs0ifmRDCqFjvC9fpNK/xEh878ASvm3ont9MQqV/bTDHqa3ndlDdwlHH/ALe9e4xFSAAwkSc2EePxZLZwMWl7RjsSjHyT21f7Tz3P22q5oNPJxkW0bceYTpGr8Ddwx3GciOIbuOM4TkTxDdxxHCei+AbuOI4TUVZSUu0TZjZuZs+cY/s9Mxs1s6dav956dbvpOI7jnM9KvFA+CeB/Avib8+wfCSH86cU01kwnUNrcv8yWOcY5vlUeZABIHOUCorKwcJwV35AXxX4BWEnk9BZFWlOjRe7Pgi64WxYh7sXtwtulg9XmzjU8HwBwZx+nB6g0+Pb1ptgbIS4U+dGaLmS7Jlkk26ZXcTj6QoM9LDZkdd/7k6zeb8py3vP/JcKn33nNM2QDgGKD7+cbs8fIlgR7Q8w1OO5dhfYDQO0zQ2QLCV4zIcFrLjWjvQnO7hSFbIWTQkycnhvTucxjVRF3H2vj9nEeDZEaAADqojAwwH3PPSty+uf4uMnbdRqAyXt5TAPdvI43dPH6mizrZ7AzxfdossRrZnaK18K6Z9lzCQBiZXFDRFH1ye089tCuALrwKlocXPmHkQseGUL4FoCpFV/RcRzHWRUu5xv4b5rZ3tYnFv065ziO41w1LnUD/yiAzQB2AhgD8GftDjSz+81st5ntrtV0RkDHcRzn4rmkDTyEcCaE0AghNAF8DMBtP+bYB0IIu0IIu5JJ/b3KcRzHuXguKZTezIZDCC8qFz8HQCtN5xGrNZA5tVzQUiHI+X3j8vwgBAN0sDDR6Bfh4BM6dByqMLAQYKZexYWOJ2/RwsTgDhZbt3dPkE2Jft+f4PzoAJASilZCJIxuiqTp65Ocs3iuyWMEgHKTRb971/DtHasWyFZp6uXUDPye0BHncPTt61ioPVla+de5lBCKJpp8f08JAXff/LC8ZscpDu9fHOQ5yo+KfO85PR9VEbneTPDajlfE+mpqgX/8dg7/7t/D/+ONzXJO7lKfDh3vfY6FwGoXj6m2oZ9szTSrsmdvE0I+gFSC1/GuQRbOHzu9gWzruzgUvh2vGjhOti+fuJVso3dph4ehH/E9nrqBxfzsJN+j6a0rE5QBYPBJXWdAccEN3Mw+DeBOAP1mdhLA7wK408x2AggAjgJ4z4pbdBzHca4IF9zAQwjvFOaPX4W+OI7jOBeBR2I6juNEFN/AHcdxIsqq5gMPyTiqQ8ujFJPTLAzEJnUkFAo6bzGdX2IRIOS1aHfi51m86ngTi5DbeljI29bBohsAHK9wPvOFOosdqoBwX0a7Wqrzr8myOJmLsWi3p8RRhrmYiEBtwxOzLB4Vqyz03N2/X54/IyIfn5rjYtI7ukfJ1inEznaUhcj96RkWqcaqLNo98Y83yWuOiGLY+RKLbovDfH8SbYrTZs+wvTzAIlcjzcedfaUujps/wwJhaYjXfLzA/ewcFVGcANJHeX0lO/leloe5TyffxHn+EwN6bd+x4SjZTi3yPerK8JpdqOuc/Fu6ONI3DiGiqluktVbMree5y43zwdlxfgZTMzof+OIgz1PuyMqFWX8DdxzHiSi+gTuO40QU38Adx3Eiim/gjuM4EWV1ixo3A+ILywXGegeLEEkhkAG6gGjytPjgb3zNA+/R4s+bb91LtlKDhYXOJItpZ2paVJ0S/d/WwSk3x6scMbq+TUrWfzh4M9luGWFBaTAzTzY1nrf07JPtLDRZqNnZeZJsO0Tq1j85eq+85sQiz8dwJ0fGqrS3Ozo5Ig/Q0Z2Haix8HVzgdLB39zxLtr2P7ZTtLIg5zp/g9Kfdz/J9W9hckNdMcjAk0i+wGCYCbSF0bwBAZpwFvrhIf1rpF+M5rEWz2VvXcPsioHBhmDuV3srXHOjktQkAvUkWN7NxdkTYUeB1uKfIYjgATFQ4bUd3kvvZcUiIrW1084bQS1WB6vS0sE3oNMDZI7xuTBRVb4e/gTuO40QU38Adx3Eiim/gjuM4EcU3cMdxnIjiG7jjOE5EWV0vlBBg5+UzjldYam906bB3paqffQMr5cW3sOJ768gReU3lXVJQbgKC8Qp7kQBAb4rP70+w14XykBjOaI+AnetZgd99mEPk/8PO74l2uJhsX1x7BCi+dnYb2f7q6/eQrb5dh0qv7eMxZYSXwaTwHFjbq71yXpXmsPuvL1xPtq15Tovwu0+9jWzrqzp+OneK10d8QRU15nehdt4dM2/uY6PIZd51jNf7vPD4AICuOve/meJjk7Mc5r14jfamSk/xPZrZxF5KmZ/mOW5W2GVjfYe+lyWRg/6ubvaSOl1nL6MmdJ7tRRFif3SGU1z07eMxFjexZwoAVHq4LXX++XscAMSn2jxvJeHyktX7n8LfwB3HcSKKb+CO4zgRxTdwx3GciHLBDdzM1pvZo2a2z8yeNbP3tuy9ZvawmR1q/b7y4oWO4zjOZbMSEbMO4P0hhCfMrBPA42b2MIB/D+CREMIfm9kHAXwQwG//2CuZkdgTE8JkTAgtANDoESH2Ip/vhsEpsm3p5PzAADBWZgFnc44LEE/Xue2DRS50DOjiqTMNPr8mwsEPzLGwCQDPTbAQ+Z5bv02269IsKJWbLMr85wNvl+3UvsEFalXO49JP8T26RoiVANArcpy/pucw2VSh40JMC8o/qqwj26JIA7Ajy/citu9usi0MaxEzVucFVrqZw9E7TrEY3/m0Ls7dfZTX/MTN/ChWuliELBzWYdaqaLeVWGytr2MBNV7WYz+7Q+QTv4tzhKtNpCfHjgQdCf1cv7rzebJ1xvj8vz17B9kWajofeE+a183it/h5Taf5vhWe13Nc6+T7Eavx+ohV+ZoLN/HzCwCZcSGSj/H+1Y4LvoGHEMZCCE+0/jwHYD+AdQDeDuDB1mEPAnjHilt1HMdxLpuL+gZuZhsB3ArgMQBDIYQXMzSdBqBfHR3HcZyrwoo3cDPrAPB5AO8LISyreRZCCNDFiWBm95vZbjPbXatpP2HHcRzn4lnRBm5mSSxt3p8KIXyhZT5jZsOtfx8GID/4hRAeCCHsCiHsSiY5UMNxHMe5NC4oYpqZAfg4gP0hhA+f808PAXg3gD9u/f7lC12rmYphYd1yAahrD+/7IcNiFABU+thez3J01Lp8kWz/PHbdhbr3L0yUOXf40UmO4mo8pyMx029lgW6uwYLQtTkWhI6XuB0A+OXNT5Dt2XkuyPy/H30z2Ya+z3M0+G2dZ/uFX2MRs/EmLjK9XohUm7tY/AWADVkWZWLGwtkNaS4S3SfELAA4XeaoPMX7nvglvuY+FpkykywsAkC5jx+Rch/PZ6wmBK7rtcitIvVSRT6uIR4DFbkMAPPbuK14RbQzxaKZ6g8AiCWLjhQLfKf3s0C345UvkG0krSMxpxr8vD29uJ5s1abIOx7X901FYnac5DWXnhaRlDUt6s6v47WQO8P3QxWTrmfaRIz289gLSfFerR/XFXmhvBbAuwA8bWZPtWy/g6WN++/N7D4AxwDwk+I4juNcNS64gYcQvgO0STgA3HVlu+M4juOsFI/EdBzHiSi+gTuO40SUVU8nG2ssF0zqAxwJqVJzAkAzyV9yslMsOBRVUeHe0/Kar+nmKLDH5zeS7cw8iw2l6zn6DQByMY44y8RYLFEi5Hcev1Fe88mDIirvEF9zqMBz1LOHRcSxt22Q7dQ2sWj474SAqqI7b8yyCAkADVEJd0OS+5Q3nrcT9YK8ZkpU/H1ilsfU+TW+b9mzHKE4t14L513HWfTLCzEsxHmM9Q6dllRJhovrhOA4w9dsKIELwMIaFSUojtvFwnvXES3avernnibb7i9wce1N3+c1c2ILZ9Z4x+CTsp28eF7UvUzEuJ8x7b2M4zMFsvUf0YL4+VQLOrozPcttzVzLgmXHSR7P7Aa9vtJFvmaIr/y92t/AHcdxIopv4I7jOBHFN3DHcZyI4hu44zhORPEN3HEcJ6KsqheKNYDk7HLvgWaG1fPELHsJAEB+ju2ldZxfpfZeDkd/5pZr5TV/sO4WvuYa4WWQFgp4Sf/8++yTd5Kt8AKfn5lgN4HNNZ0zubiJ1W6Vnzh/ms/f/z4OOx8a0fnR81X2nFBpAO7t2ku2zpgo0ArgUJUTVV6X5PD875Q4fLozrj0HHineRLaTf7iFbIVZPr/aw14G5V4dq1bp4dzfuTN8LzuPcf7paudFeFOdYVtV1BpWRYkBYOAJts9cx89GJ6dHx/gd2pNjpMaeE/NbRCj9Hfxc3nftD8mmPIcA4GxdeMYk+JrpGJ8/XeX7AwCzh9gLptDBYfchxvNuDT0fuTPcJ+WFMr+O15fKAQ/odAfVgvZeUvgbuOM4TkTxDdxxHCei+AbuOI4TUXwDdxzHiSirKmLGqnWkjy8Poa4PslJjdR3a28ixONBMtEuUuJzevbrgbt8PWJgobyhw21kWDBcHdNuFQyycxUss/ixsYJFprlv/TO17WhStFfN06AM8R/du2Ue2Yhvx56bOMbJtzbBttC5EoriuuKRyf/9N8ZVk6xRFjX/3qbfJa3Y8wnPX1RRCkQpxF7mZlZgEAA1x7Nx6vkexOs9nVojUADC/lu9R50kW6Kod3E673N3VAguOyRLP++RGfuTvuX2PvKZaI7/+6n8i20SN0xWodBLHqpxrvt35TZEA9UyZxc7ZikhaDuCa/8tzn5ri57I0zOtI5YAHgHqWxUWRUQKVHpEvvq6f6+paPrYi0mHgIXm6v4E7juNEFd/AHcdxIopv4I7jOBHlghu4ma03s0fNbJ+ZPWtm723Zf8/MRs3sqdavt1797jqO4zgvshIRsw7g/SGEJ8ysE8DjZvZw698+EkL40xW3FgBrLBdWVM7kRkZ3KyaKjWbGWYSs9LP4kp7QEX3NLj42Oc8CSLwi2j6jr1ntEcWXcyyC5k5y9F7H81r4ik2xiLn/QyNk2zI0SraEiIC7p49zPQPAYILb+doM54B+W4FzO8fb5Gb++iyfP5jidv726G1k2/A/2kQzplgkS4+yUG0VPm5uZC3ZGmktSKt8zdUuPra4RUTFjuq+ZydYbG2m+Jqdx3ltx46yoAwAtTdw0e65Ee5Tfec82X50RueGH8jzsWNVjup9Rf6oPP98DixslvbpCufv/9XhH5DtS9VXkK2Q4mcIAH7wFs61nxvj51IMp21+9NQcr4V4lW29+/kZnrpR5wNviPuendDPkWIlNTHHAIy1/jxnZvsBrFtxC47jOM5V4aK+gZvZRgC3AnisZfpNM9trZp8wM/YrcxzHca4aK97AzawDwOcBvC+EMAvgowA2A9iJpTf0P2tz3v1mttvMdleb+r87juM4zsWzog3czJJY2rw/FUL4AgCEEM6EEBohhCaAjwHgD5hLxz0QQtgVQtiVivG3LsdxHOfSuOA3cDMzAB8HsD+E8OFz7MOt7+MA8HMAnrnQtUIyjtpwYZktvsiCTmJep8wsD/EPgNxBLo4b6+OIrXbERyfIFjq4nViWo+fqXToKLDXD/VeRclPbOQKt45ROO5lO8a3KrmGRSTFV5WizzanxFZ0LAO/q/R7ZHpq9lWyTVR4PACw0eO6+cfgGsl3zYZFS9SJSazZ6+L4F47F3v8D/E0wu6sjUeFWI10URaVdj4amRahN91y1SKIuoyaY4v3rzNfKa+eMcBTv6s7zm+jtYeE/GdZrXRpPbj4P72ZfgdVgLvF5nqvp5OVsSUbUiNfH1HWfI9uiZrfKaaXGPmqJWce4037dEuU20q0gPXJZRlzz2/j3668PEDl6zSiRvx0q8UF4L4F0Anjazp1q23wHwTjPbiaUi20cBvGfFrTqO4ziXzUq8UL4DiMQEwFevfHccx3GcleKRmI7jOBHFN3DHcZyI4hu44zhORFnVfOAwLiLaEEWNUye1d0WiU4ejnk+lj9Xuahe3AwCpNayAq6KzqVkOj6116OmLC4+CxAJ7l/Q9zUp7bFF74MxdXyBbCOx50Jnia97efYRsZxuiYi6AjPE4n6isIdupMvdnWweH8QPAn3/vbrJt+WseZ3ySw+vTNe0hUbyBPY0aab7vuVPsdRE7coptAzrMW9XhTZbYGGuI8OucXh+VLrYnRFYGVXC3HVPbeD5Ghk+Trdrg5yCT0J5PW7q48PX6DHt9Ha5w0epHJtjLKBXX7dy7lvPVr4lz0evPHd5Jts7P6nXcLdbN4qB4XxVTXMvp91qR8l2uD5V2o929VPe9rNOmS/wN3HEcJ6L4Bu44jhNRfAN3HMeJKL6BO47jRJRVFTFDzFDPL28yM8aCZbNL50yJVVgECSkOtc49wyJVLi3iaAFM38YCnTU4lDZ3XAiOJT19tS7RlvhRWc9z361NLvSyKna8hwWc6Tey+NMQSs14XYs/FVGldVHEIBeSHBr8sU/pmh43fvx5soUZ7mfYupH7M6DXgsqVnTzD15x+5QDZ0v2b5DUViwN8P/JnRIHqARbYM5M6t3tmhpWvihDZGyLy3ESYNgBMvJ7bGhKCpQqPzye1cH5djkPXVbHiz5/iPN2VBvdzbV4XFt+VY5H9T07dQ7aBv2AVMVXUDg+q4Hd+lPtUGuJJThf1fUvM8JoLCZ7PZlI7TCimtosc49xMW/wN3HEcJ6L4Bu44jhNRfAN3HMeJKL6BO47jRJRVFTGtGZAoLxdwQnrl+Z7lNRdZXIQQNps5nYu48LX9fM1OjmprdnOu69hiG7FDCBtlIXLFyyy0KAEWADrTXKS1+ygLINO38TjHqgWyDSVZ8AOA0Qof+4XdryTbNQ/xuQPQ81G7gUuoJop9ZDMRPZc5ydGZAGBlVnrOzzUPAD1PcL73RrcIqWsTKZc/yGJts1MIX2dESJ4QwwGgvL0g7eejiuhWuttEZ4pA0JTI873Q4LU5WWojFIuc3N+b4eLJ8zUWuX91w4/IdrzC9xwA/nL0TrLt+y4LzZtPT5MtNqNFzMZQgWxqfXUc4MjS2prLqycQ6+VKyc02zglrv803LjnH/WQ3gFZbP757juM4zssV38Adx3Eiim/gjuM4EeWCG7iZZczsh2a2x8yeNbPfb9mvNbPHzOx5M/usmelIGcdxHOeqsBIRswLgzSGE+VZ1+u+Y2T8C+I8APhJC+IyZ/SWA+wB89MddqJmKYWF4uZgXr/C+rwrJtqORGySbBRHdtKAFNhtgYaXRxSJXrZuFq1hVpzptpjgSq9TLtpBgW+6Q/jmYOV4km9W5/b7/ygLKI9vuIFtMZ/ZE9xEWrq5r8sG1Tl46SZFyFwAaIq1q6OE5VsWsraw7avOiMPG4iDwUhY6r3aJAdV5Hz6WzIvXrLAuosRlO7RvaRP8W9rJwdvoNvA7rgQXLzJR+NmJpXgsNcf7r1x4m20PfZZEaAP7o+beTrW8jC4l3rOFISiWcL9R1OujRB4Vg+XiRbCHLzgnVgs69mpzmPK2l9Rx9nJ0XEb2ndMSolXl9NvsKbMtxP9s5PKRmRVH3OR0Zq7jgG3hY4kWpN9n6FQC8GcDnWvYHAbxjxa06juM4l82KvoGbWbxVkX4cwMMAXgBQDCG8+OPjJAD2FXMcx3GuGivawEMIjRDCTgAjAG4DwOU22mBm95vZbjPbXSvzfzMdx3GcS+OivFBCCEUAjwK4A0DBzF78QDgCQNbTCiE8EELYFULYlcxw+TLHcRzn0liJF8qAmRVaf84CuBvAfixt5L/YOuzdAL58lfroOI7jCFbihTIM4EEzi2Npw//7EMJXzGwfgM+Y2R8CeBLAxy90oVgtIHdmucLajLNS3kzrnytWF6HFvaz4pqdFAeIurYCrQP7YHHtihB72QglJ3c/kDJ8/8N0i2ZrdOoRZ0RSeMVYVCvYstz34PVbky2t1uHAzLvIbZ/gepadYvW+0CRee3cCznFzke5meEV45z7PXAwDAuE9W4/mITwtvFRE237aAsIiGVyHZU69ZS7Z0UXsp5V7gMa351iTZFq9lj6Kx16w8+0W5xsduznCh4sHruG0AmNrDudR7szyfExVOMzFX5+flnx+/SbazriiKAIuc2rEif4JduJn7COjc8NkTIgd9lj2Fqr0i1QKA1CSPPaR4juc28ZeG0OZVueOE8GgS+087LrgaQgh7Adwq7Iex9D3ccRzHeQnwSEzHcZyI4hu44zhORPEN3HEcJ6Ksaj5wNANi5+UDTwghbva6NgKb6K0Sw0r9LJrFqzo3c+ZQkWy19RzW3JSCpQ5rbnSwMNIUYkfyOAtKSpwDdEh5bLLIB4rw7TDHOZOzs9onX+VRbqaFUHMNC7DtiviamKZSP89nXYilzWSvvGb+CN/jyZ0FstU6VxaOHqvp9RESfH5pkNtRYlRihsVjQOcjj5XF3IkupbcV5TVDSRSezrIYNtPgtm/s5eLFALDwehZbu5N8zbESh6gf+9K1ZNv6qM5Br9Z8bIYFQ0wVydT1hE610OxmIbEm0jcoSgO6RkEjy2KtCZ16YYjXdnZC7xUzm1jsHRhtI9wL/A3ccRwnovgG7jiOE1F8A3ccx4kovoE7juNElNUtatxocqRglcWbWJ3FAgDITLJgUenhITRSLIpUO/XPqrCLkyhWO/jYtCg02kzrHNLzBZHjvMKKVEdiiGyp58fkNetCgKmtZ7G31sl9Ss5z33Pf12VSY2VRvFnkTE6IIs0La3X+63qO70e8zPORELaOg0V5TUV6VhSJPss2VUy63KcfhexpFiIrvTz2Sp8QztvkMo+/wGmDrJPnvdxT4HYqWmDryLO4eNfgc9y2UJTPlLTTwLZuXotPTY+Qbewb68nWe1A8q/1aRMw9K9a8yOmPbhZLm11topmb6h7z3Cnhvfs5LbYGETU+s1XMndhqSoN6/0ks8DhrQxyBixfk6f4G7jiOE1V8A3ccx4kovoE7juNEFN/AHcdxIsrqRmKGACstF8RUGs92UZP1DhbouvdMkK3ZydFNqigxACRFgdqO+srSW85s0QUqVHRWepaNMZGWtB3J/ce5na0sKCUX+WdyQhR4bW7i9KcAEBPCW0iKwr4LF9H3BZH6VQSmxUS64EaXvm+J00WydRxm8UlFkSpbuwhYJaZlD3NRYpWWtE2gLiwjxiTaV8KX0vYA4K6Rg2SrBV6zj57dSraetIh6BLB3mgX+qc/ymht6np+hWhe3Pb+23XYzTJbsbi6+XLuWC5irdNIAkBvlMXUdKJLNpriAcXNIR/+WB0WBbBHpm57im1RvEwTaSIv7vkaveYW/gTuO40QU38Adx3Eiim/gjuM4EWUlNTEzZvZDM9tjZs+a2e+37J80syNm9lTr186r3lvHcRznX1iJiFkB8OYQwryZJQF8x8z+sfVvvxVC+NzV657jOI7TjpXUxAwAXkwonWz9aqOFX+hi4BDXhChKPMGqNgBURYh6ZaRANhXC3Ezp/2yUhllZVqH42TPcJ+VJAQC5KW5fhfIn50UofraNAp1nGTsuPGhUjvDGiCj82sadYWZbD9lypzhMO7Eg8pM39DUbIuVAYp5DmBfX8tjbFUqOz3M+c0sJj4Qs25T3j8pFDgBnX8Gh0rmz7H1kYj5jIn0CACQLPM7k+BzZ+p7m+zvxRp1zvSKS5U9VuZ+dIp/3kSLnvweA8CW2D+zh3PK1Tn4uVWHxhEjp0I7aDRyeL3Pyt3EeauT5vsvC1Z3s2VIe0M9g7gh7OVW7CmRbWCs8wUTIPACki5e2lb7Iir6Bm1nczJ4CMA7g4RDCY61/+iMz22tmHzEzXfbdcRzHuSqsaAMPITRCCDsBjAC4zcy2A/gQgBsAvApAL4DfVuea2f1mttvMdlcb2t/UcRzHuXguygslhFAE8CiAe0IIY2GJCoC/BnBbm3MeCCHsCiHsSsXbZA5zHMdxLpqVeKEMmFmh9ecsgLsBHDCz4ZbNALwDwDNXr5uO4zjO+VhoF5v74gFmtwB4EEAcSxv+34cQ/sDMvglgAEsywlMAfj2EwArH8mudBXCs9dd+ABwHH118PC9/ftLG5ON5eXMlx3NNCIG8ES64gV8tzGx3CGHXS9L4VcDH8/LnJ21MPp6XN6sxHo/EdBzHiSi+gTuO40SUl3IDf+AlbPtq4ON5+fOTNiYfz8ubqz6el+wbuOM4jnN5+CcUx3GciLLqG7iZ3WNmz5nZ82b2wdVu/0pgZp8ws3Eze+YcW6+ZPWxmh1q/c1KRlylmtt7MHjWzfa2Mk+9t2SM5ph+TQfNaM3ustfY+a2aijM7Ll1ZKiyfN7Cutv0d9PEfN7OlWNtPdLVsk1xwAmFnBzD5nZgfMbL+Z3XG1x7OqG7iZxQH8BYB7AdwE4J1mdtNq9uEK8UkA95xn+yCAR0IIWwA80vp7VKgDeH8I4SYArwbwG637EtUxvZhBcweAnQDuMbNXA/hvAD4SQrgOwDSA+166Ll4S7wWw/5y/R308APCmEMLOc9ztorrmAODPAXwthHADgB1YuldXdzwhhFX7BeAOAF8/5+8fAvCh1ezDFRzLRgDPnPP35wAMt/48DOC5l7qPlzG2L2Mp4jbyYwKQA/AEgNuxFFSRaNmXrcWX+y8s5SF6BMCbAXwFSwF0kR1Pq89HAfSfZ4vkmgPQDeAIWrriao1ntT+hrANw4py/n2zZfhIYCiGMtf58GsDQS9mZS8XMNgK4FcBjiPCYzs+gCeAFAMUQwot5cKO29v47gA/gX0sl9yHa4wGWEkx/w8weN7P7W7aorrlrAZwF8Netz1x/ZWZ5XOXxuIh5FQhLP24j595jZh0APg/gfSGEZcmPozamcF4GTSxlzowkZvYzAMZDCI+/1H25wrwuhPAKLH1S/Q0ze8O5/xixNZcA8AoAHw0h3ApgAed9Lrka41ntDXwUwLmZ2kdatp8EzpyT4GsYS29+kaFVbenzAD4VQvhCyxzpMQHLMmjeAaBgZi9WPojS2nstgJ81s6MAPoOlzyh/juiOBwAQQhht/T4O4ItY+kEb1TV3EsDJ8K+1Ej6HpQ39qo5ntTfwHwHY0lLPUwB+BcBDq9yHq8VDAN7d+vO7sfQdORK0Mkp+HMD+EMKHz/mnSI6pTQbN/VjayH+xdVhkxhNC+FAIYSSEsBFLz8w3Qwj/FhEdDwCYWd7MOl/8M4CfwlJG00iuuRDCaQAnzOz6lukuAPtwtcfzEnzsfyuAg1j6JvmfXmrx4RLH8GkAYwBqWPrJex+Wvkk+AuAQgP8HoPel7udFjOd1WPqv3V4sZZZ8qnWfIjkmALcAeLI1nmcA/JeWfROAHwJ4HsD/AZB+qft6CWO7E8BXoj6eVt/3tH49++JeENU11+r7TgC7W+vuSwB6rvZ4PBLTcRwnoriI6TiOE1F8A3ccx4kovoE7juNEFN/AHcdxIopv4I7jOBHFN3DHcZyI4hu44zhORPEN3HEcJ6L8f1n07dbM4zZZAAAAAElFTkSuQmCC)
%% Cell type:code id: tags:
``` python
#Generate receptive fields of several RGCs #2304 RGCs - haf ON/half OFFù
import random
rf_ON = []
rf_ON_center_coord = []
rf_OFF = []
rf_OFF_center_coord = []
i = 0
image = images[0]
for width_center in range(image.shape[1]):
for height_center in range(image.shape[0]):
if (i % 2) == 0:
#theta = random.uniform(0, 6.2)
rf = RF_circular_CS(image.shape[1], image.shape[0], width_center, height_center, 1, plot=False)
rf_ON.append(rf)
rf_ON_center_coord = np.array([width_center, height_center])
save('D://inception_loop/RGC_sim_data/data/static23032021_small_set_original/RFs/'+str(i)+'.npy', rf)
save('D://inception_loop/RGC_sim_data/data/static23032021_small_set_original/RFs/center_coord_'+str(i)+'.npy',rf_ON_center_coord)
else:
#theta = random.uniform(0, 6.2)
rf = RF_circular_CS(image.shape[1], image.shape[0], width_center, height_center, -1, plot=False)
rf_OFF.append(rf)
rf_OFF_center_coord = np.array([width_center, height_center])
save('D://inception_loop/RGC_sim_data/data/static23032021_small_set_original/RFs/'+str(i)+'.npy', rf)
save('D://inception_loop/RGC_sim_data/data/static23032021_small_set_original/RFs/center_coord_'+str(i)+'.npy',rf_OFF_center_coord)
i+=1
```
%% Cell type:code id: tags:
``` python
#Generate responses of simulated RGCs to the image set from evaluation data set of lurz2020
import time
start_time = time.time()
i=0
for image in images:
responses = []
for rfon, rfoff in zip(rf_ON, rf_OFF):
rgc_on_response = RGC_response(rf=rfon, image=image, plot=False)
responses.append(rgc_on_response)
rgc_off_response = RGC_response(rf=rfoff, image=image, plot=False)
responses.append(rgc_off_response)
# save numpy array as npy file
save('D://inception_loop/RGC_sim_data/data/static23032021_small_set_original/data/responses/'+str(i)+'.npy', responses)
if (i % 1000) == 0:
print(i)
i+=1
print("--- %s seconds ---" % (time.time() - start_time))
```
%%%% Output: stream
0
1000
2000
3000
--- 3275.217703819275 seconds ---
%% Cell type:code id: tags:
``` python
#Retrieve responses of simulated RGCs to the image set
import time
start_time = time.time()
responses_all = []
for k in range(len(images)):
responses_all.append(np.load(paths_new+'/data/responses/'+str(k)+'.npy'))
if (k % 1000) == 0:
print(k)
print("--- %s seconds ---" % (time.time() - start_time))
```
%%%% Output: stream
0
1000
2000
3000
--- 1.4416766166687012 seconds ---
%% Cell type:code id: tags:
``` python
#Generate data - pupil_center npy files
#Array with list of two values- coordinates of pupil center - [759.87785056, 472.71767702]
for j in range(len(images)):
pupil_center = np.array([0.0,0.0])
# save numpy array as npy file
save(paths_new+'/data/pupil_center/'+str(j)+'.npy', pupil_center)
```
%% Cell type:code id: tags:
``` python
#Generate metadata - neurons - unit_ids npy array
unit_ids = np.array(range(1, len(responses)+1))
#Generate metadata - neurons - animal_ids npy array
animal_ids = np.repeat(1, len(responses))
#Generate metadata - neurons - area npy array
area = ['retina']*len(responses)
#Generate metadata - neurons - layer npy array
layer = ['RGC']*len(responses)
#Generate metadata - neurons - scan_idx npy array
scan_idx = np.repeat(14, len(responses))
#Generate metadata - neurons - sessions npy array
sessions = np.repeat(6, len(responses))
# save numpy arrays as npy arrays
save(paths_new+'/meta/neurons/unit_ids.npy', unit_ids)
save(paths_new+'/meta/neurons/animal_ids.npy', animal_ids)
save(paths_new+'/meta/neurons/area.npy', area)
save(paths_new+'/meta/neurons/layer.npy', layer)
save(paths_new+'/meta/neurons/scan_idx.npy', scan_idx)
save(paths_new+'/meta/neurons/sessions.npy', sessions)
```
%% Cell type:code id: tags:
``` python
#Generate metadata - statistics - responses - all
responses_max_all = np.max(responses_all, axis=0)
responses_mean_all = np.mean(responses_all, axis=0)
responses_median_all = np.median(responses_all, axis=0)
responses_min_all = np.min(responses_all, axis=0)
responses_std_all = np.std(responses_all, axis=0)
#Generate metadata - statistics - responses - stimulus_frame
# save numpy arrays as npy arrays - all
save(paths_new+'/meta/statistics/responses/all/max.npy', responses_max_all)
save(paths_new+'/meta/statistics/responses/all/mean.npy', responses_mean_all)
save(paths_new+'/meta/statistics/responses/all/median.npy', responses_median_all)
save(paths_new+'/meta/statistics/responses/all/min.npy', responses_min_all)
save(paths_new+'/meta/statistics/responses/all/std.npy', responses_std_all)
# save numpy arrays as npy arrays - stimulus_frame
save(paths_new+'/meta/statistics/responses/stimulus_frame/max.npy', responses_max_all)
save(paths_new+'/meta/statistics/responses/stimulus_frame/mean.npy', responses_mean_all)
save(paths_new+'/meta/statistics/responses/stimulus_frame/median.npy', responses_median_all)
save(paths_new+'/meta/statistics/responses/stimulus_frame/min.npy', responses_min_all)
save(paths_new+'/meta/statistics/responses/stimulus_frame/std.npy', responses_std_all)
```
%% Cell type:code id: tags:
``` python
#Generate metadata - statistics - responses - all
images_max_all = np.max(images)
images_mean_all = np.mean(images)
images_median_all = np.median(images)
images_min_all = np.min(images)
images_std_all = np.std(images)
#Generate metadata - statistics - responses - stimulus_frame
# save numpy arrays as npy arrays - all
save(paths_new+'/meta/statistics/images/all/mean.npy', images_mean_all)
save(paths_new+'/meta/statistics/images/all/median.npy', images_median_all)
save(paths_new+'/meta/statistics/images/all/min.npy', images_min_all)
save(paths_new+'/meta/statistics/images/all/std.npy', images_std_all)
# save numpy arrays as npy arrays - stimulus_frame
save(paths_new+'/meta/statistics/images/stimulus_frame/max.npy', images_max_all)
save(paths_new+'/meta/statistics/images/stimulus_frame/mean.npy', images_mean_all)
save(paths_new+'/meta/statistics/images/stimulus_frame/median.npy', images_median_all)
save(paths_new+'/meta/statistics/images/stimulus_frame/min.npy', images_min_all)
save(paths_new+'/meta/statistics/images/stimulus_frame/std.npy', images_std_all)
```
%% Cell type:code id: tags:
``` python
#Generate metadata - trials - animal_id npy array
animal_id = np.repeat(1, len(responses_all))
#Generate metadata - trials - condition_hash npy array
condition_hash = np.repeat(" ", len(responses_all))
#Generate metadata - trials - frame_image_class npy array
frame_image_class = np.repeat("imagenet", len(responses_all))
#Generate metadata - trials - frame_image_id npy array
frame_image_id = np.load(paths+'/meta/trials/frame_image_id.npy')
frame_image_id_new = np.delete(frame_image_id, remove_ids)
#Generate metadata - trials - frame_last_flip npy array
frame_last_flip = np.random.randint(11000, 30000, size=(len(responses_all)))
#Generate metadata - trials - frame_pre_blank_period npy array
frame_pre_blank_period = np.random.uniform(0.3, 0.5, size=(len(responses_all)))
#Generate metadata - trials - frame_presentation_time npy array
frame_presentation_time = np.repeat(0.5, len(responses_all))
#Generate metadata - trials - frame_trial_ts npy array
frame_trial_ts = np.repeat("Timestamp('2021-05-13 13:14:43')", len(responses_all))
#Generate metadata - trials - scan_idx npy array
scan_idx = np.repeat(14, len(responses_all))
#Generate metadata - trials - tiers npy array
#Generate metadata - trials - session npy array
session = np.repeat(6, len(responses_all))
#Generate metadata - trials - trial_idx npy array
trial_idx = np.repeat(0, len(responses_all))
# save numpy arrays as npy arrays
save(paths_new+'/meta/trials/animal_id.npy', animal_id)
save(paths_new+'/meta/trials/condition_hash.npy', condition_hash)
save(paths_new+'/meta/trials/frame_image_class.npy', frame_image_class)
save(paths_new+'/meta/trials/frame_image_id.npy', frame_image_id_new)
save(paths_new+'/meta/trials/frame_last_flip.npy', frame_last_flip)
save(paths_new+'/meta/trials/frame_pre_blank_period.npy', frame_pre_blank_period)
save(paths_new+'/meta/trials/frame_presentation_time.npy', frame_presentation_time)
save(paths_new+'/meta/trials/frame_trial_ts.npy', frame_trial_ts)
save(paths_new+'/meta/trials/scan_idx.npy', scan_idx)
save(paths_new+'/meta/trials/tiers.npy', tiers_new)
save(paths_new+'/meta/trials/session.npy', session)
save(paths_new+'/meta/trials/trial_idx.npy', trial_idx)
```
%% Cell type:code id: tags:
``` python
```
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment