Commit 2aee9566 authored by Bryce Hepner's avatar Bryce Hepner

Creating 2d hist, trying to find good scale

parent f52387ba
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from WorkingPyDemo import *\n",
"import matplotlib.colors as colors"
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {},
"outputs": [],
"source": [
"scenes = file_extractor(\"images\")\n",
"images = image_extractor(scenes)\n",
"newnamesforlater = []\n",
"# list_dic, bins = make_dictionary(images, 4, False)\n",
"file_sizes_new = []\n",
"file_sizes_old = []\n",
"list_dic = np.load(\"first_dict.npy\", allow_pickle=\"TRUE\")\n",
"bins = [21,32,48]\n",
"# np.save(\"first_dict.npy\", list_dic)\n",
"errors = np.array([])\n",
"diffs = np.array([])\n",
"for i in range(len(images)):\n",
" image, diff, new_error = predict_pix(images[i], False)\n",
" errors = np.concatenate((errors, abs(new_error)))\n",
" diffs = np.concatenate((diffs, abs(diff)))"
]
},
{
"cell_type": "code",
"execution_count": 68,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAD4CAYAAADmWv3KAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAANXElEQVR4nO3db2xk51XH8e/pLk1EEW5CVlXYDXgrrxatKkErKyQUpKpAu2lqglBfZFuJCFaxhIgofyTYiBdR1TdBQpQCURWrpEEVbAglgjhZsYJQlEiNQrwCld1uQ9w0JY5a1iXB/HuRBA4v5jqduLZ3xjP2nTnz/UhWPM/Mes7jJ/r5+tzH90ZmIkmq501tFyBJ2h0GvCQVZcBLUlEGvCQVZcBLUlH72y4A4Jprrsnp6em2y5CksXLu3LlvZuaBrZ4fiYCfnp5maWmp7TIkaaxExNe2e94WjSQVZcBLUlGtBnxEzEXEwtraWptlSFJJrQZ8Zi5m5vzU1FSbZUhSSbZoJKkoA16SijLgJakoA16SihqJP3QaxPSpRzcdf/7um/e4EkkaLR7BS1JRBrwkFWXAS1JRBrwkFWXAS1JRBrwkFWXAS1JRXk1SkoryapKSVJQtGkkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKL2D/sLRsSPAR9pvvaxzPyRYb+HJOnyegr4iLgP+CBwKTPf0TV+HPgksA/4dGbenZlPAE9ExE8DTw+/5N5Mn3p00/Hn7755jyuRpHb02qK5HzjePRAR+4B7gJuAY8CJiDjW9ZIPA38yhBolSTvQU8Bn5uPASxuGrweWM/O5zHwFeAC4BSAivg9Yy8z/HGaxkqTeDXKS9SDwQtfjlWYM4CTwme3+cUTMR8RSRCytrq4OUIYkaTO7sosmM+/KzC9c5jULmTmbmbMHDhzYjTIkaaINEvAvAtd1PT7UjEmSRsAgAf80cCQiDkfEm4FbgYf7+QIRMRcRC2trawOUIUnaTE8BHxGngSeBoxGxEhEnM/M14A7gLHAReDAzL/Tz5pm5mJnzU1NT/dYtSbqMnvbBZ+aJLcbPAGeGWpEkaSi8VIEkFWXAS1JRrQa8J1klafe0GvCeZJWk3WOLRpKKMuAlqSgDXpKK8iSrJBXlSVZJKsoWjSQVZcBLUlEGvCQV5UlWSSqqp6tJ7pbMXAQWZ2dnb9+r95w+9eiWzz1/9817VYYk7TpbNJJUlAEvSUW12qIZNVu1b2zdSBpHHsFLUlEGvCQV1WqLJiLmgLmZmZk2y7gsWzeSxpHXopGkomzRSFJRBrwkFWXAS1JRBrwkFWXAS1JRBrwkFeXlgiWpKPfBS1JRtmgkqSgDXpKKMuAlqSgDXpKKMuAlqSgDXpKKMuAlqSgDXpKK8i9ZJamoVm/Zl5mLwOLs7OztbdaxU97KT9Ios0UjSUW1egRflUf2kkaBR/CSVJQBL0lF2aLZQ7ZuJO0lj+AlqSgDXpKKskUzAmzdSNoNHsFLUlEGvCQVZYtmhNm6kTQIj+AlqSgDXpKKarVFExFzwNzMzEybZYwdWzeSetHqEXxmLmbm/NTUVJtlSFJJtmgkqSgDXpKKMuAlqSj3wU8AT8pKk8kjeEkqyoCXpKIMeEkqyoCXpKI8yapvs9VJWfDErDROPIKXpKIMeEkqyoCXpKIMeEkqyoCXpKIMeEkqyoCXpKLcB69d5YXOpPZ4BC9JRRnwklSUAS9JRQ29Bx8RbwI+Dnw3sJSZfzTs95AkXV5PR/ARcV9EXIqI8xvGj0fEMxGxHBGnmuFbgEPAq8DKcMuVJPWq1xbN/cDx7oGI2AfcA9wEHANORMQx4Cjwhcz8VeAXhleqJKkfPQV8Zj4OvLRh+HpgOTOfy8xXgAfoHL2vAC83r/nfYRUqSerPICdZDwIvdD1eacYeAt4fEb8PPL7VP46I+YhYioil1dXVAcqQJG1m6CdZM/N/gJM9vG4BWACYnZ3NYdehyeAfUklbG+QI/kXguq7Hh5oxSdIIGOQI/mngSEQcphPstwIf7ucLRMQcMDczMzNAGVq33a32JE2eXrdJngaeBI5GxEpEnMzM14A7gLPAReDBzLzQz5tn5mJmzk9NTfVbtyTpMno6gs/ME1uMnwHODLUiSdJQeDXJCWZLR6rNa9FIUlGtBnxEzEXEwtraWptlSFJJrQa8J1klaffYopGkogx4SSrKgJekolrdJulfso4ft1ZK48OTrJJUlC0aSSrKgJekogx4SSrKv2SVpKI8ySpJRXk1SbWi3+2W3oJP6p89eEkqyoCXpKIMeEkqyh68xoKXSJD65zZJSSrKbZKSVJQ9eEkqyh68StqqZ+9+ek0Sj+AlqSgDXpKKMuAlqSh78JooXgNHk8R98JJUlPvgJakoWzTSNrZr6di+0agz4KUdcq+9Rp27aCSpKANekoqyRSMNma0bjQoDXtojBr/2mi0aSSrKgJekolpt0UTEHDA3MzPTZhlSq7x8gnaLf8kqSUXZopGkogx4SSrKbZLSmHG7pXplwEtFGPzayBaNJBVlwEtSUbZopOLcZz+5PIKXpKIMeEkqyoCXpKLswUt6A7db1mHAS+qJwT9+bNFIUlGtBnxEzEXEwtraWptlSFJJrbZoMnMRWJydnb29zTok7Vy/++zBts5esUUjSUUZ8JJUlAEvSUW5TVLSyHAr5nAZ8JJGnsG/M7ZoJKkoA16SirJFI2ls2brZngEvqRyDv8MWjSQVZcBLUlG2aCRNjElr3RjwkrSFcf+BYItGkooy4CWpKANekoqyBy9JfRqXm5x4BC9JRRnwklSUAS9JRQ29Bx8R7wE+DlwAHsjMvxv2e0jSuGljT31PR/ARcV9EXIqI8xvGj0fEMxGxHBGnmuEE/gu4ElgZbrmSpF712qK5HzjePRAR+4B7gJuAY8CJiDgGPJGZNwG/AXxseKVKkvrRU8Bn5uPASxuGrweWM/O5zHwFeAC4JTP/r3n+ZeCKrb5mRMxHxFJELK2uru6gdEnSdgY5yXoQeKHr8QpwMCJ+JiLuBT4L/MFW/zgzFzJzNjNnDxw4MEAZkqTNDP0ka2Y+BDw07K8rSerPIEfwLwLXdT0+1IxJkkbAIAH/NHAkIg5HxJuBW4GH+/kCETEXEQtra2sDlCFJ2kyv2yRPA08CRyNiJSJOZuZrwB3AWeAi8GBmXujnzTNzMTPnp6am+q1bknQZPfXgM/PEFuNngDNDrUiSNBSRmW3XQESsAl/b4T+/BvjmEMsZF5M470mcMzjvSdLvnL8/M7fchjgSAT+IiFjKzNm269hrkzjvSZwzOO+269hLw56zFxuTpKIMeEkqqkLAL7RdQEsmcd6TOGdw3pNkqHMe+x68JGlzFY7gJUmbMOAlqaixDvgtbjgy9iLiuoj4fER8KSIuRMRHm/GrI+KvI+LZ5r9XNeMREb/XfB++GBHvancGg4mIfRHxDxHxSPP4cEQ81czvT5tLYxARVzSPl5vnp1stfIci4q0R8bmI+HJEXIyIGydhrSPiV5r/v89HxOmIuLLiWm92w6SdrG9E3Na8/tmIuK2X9x7bgN/mhiMVvAb8WmYeA24AfrGZ2yngscw8AjzWPIbO9+BI8zEPfGrvSx6qj9K5/MW63wI+kZkzdO4zcLIZPwm83Ix/onndOPok8FeZ+QPAD9KZe+m1joiDwC8Bs5n5DmAfnetZVVzr+9lwwyT6XN+IuBq4C/hhOvfiuGv9h8K2MnMsP4AbgbNdj+8E7my7rl2a618CPwk8A1zbjF0LPNN8fi9wouv1r79u3D7oXJX0MeC9wCNA0PnLvv0b153OdZBubD7f37wu2p5Dn/OdAr66se7qa8237idxdbN2jwDvr7rWwDRwfqfrC5wA7u0af8PrtvoY2yN4trjhSEu17JrmV9F3Ak8Bb8vMrzdPfQN4W/N5pe/F7wK/DqzfGex7gH/PzsXt4I1ze33ezfNrzevHyWFgFfhM05b6dES8heJrnZkvAr8N/AvwdTprd47aa92t3/Xd0bqPc8CXFxHfBfw58MuZ+R/dz2Xnx3ipPa4R8UHgUmaea7uWPbQfeBfwqcx8J/DffOvXdaDsWl8F3ELnB9z3Am/h29sYE2E313ecA770DUci4jvohPsfZ+cuWQD/GhHXNs9fC1xqxqt8L94N/FREPE/nHr/vpdOffmtErF/5tHtur8+7eX4K+Le9LHgIVoCVzHyqefw5OoFffa1/AvhqZq5m5qt07gL3bmqvdbd+13dH6z7OAT/wDUdGVUQE8IfAxcz8na6nHgbWz57fRqc3vz7+s80Z+BuAta5f/8ZGZt6ZmYcyc5rOev5tZn4E+DzwoeZlG+e9/v34UPP6sTrSzcxvAC9ExNFm6MeBL1F8rem0Zm6IiO9s/n9fn3fZtd6g3/U9C7wvIq5qfvt5XzO2vbZPPgx44uIDwD8DXwF+s+16hjivH6XzK9sXgX9sPj5Ap+f4GPAs8DfA1c3rg86Ooq8A/0RnZ0Lr8xjwe/Ae4JHm87cDfw8sA38GXNGMX9k8Xm6ef3vbde9wrj8ELDXr/RfAVZOw1sDHgC8D54HPAldUXGvgNJ3zDK/S+Y3t5E7WF/j5Zv7LwM/18t5eqkCSihrnFo0kaRsGvCQVZcBLUlEGvCQVZcBLUlEGvCQVZcBLUlH/D+BOEmtiI4pzAAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAERCAYAAABrWly6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABYk0lEQVR4nO29e7g1WV3f+fmuqr3P9b322/c72KBAuLYExTEEYrjEAJkxCYxRTJjpGYdETZhxROfRJM9kok8Sb49Gp0eI6MOgBo32GKMhBDXm0dYGAbmIdMRAQzdN3/vt93LOrvrNH1W1d+3adVl1Tu2zzzmsz/vs5z27au31W1W11q9WrVq/75KZEQgEAoGjj1t1AQKBQCAwDMGhBwKBwDEhOPRAIBA4JgSHHggEAseE4NADgUDgmBAceiAQCBwTgkMPBAKBEpLeIelBSR/1TP83JH1c0sck/b/LLl9rWcI89EAgEJgh6WuB88DPmNlzOtLeBvwC8HIze1TSVWb24EGUs47QQw8EAoESZvbbwCPlbZKeLunXJX1A0n+S9OX5rv8R+HEzezT/7cqcOQSHHggEAj7cCfw9M3sR8L8C/zLf/gzgGZL+s6Tfk/SqlZUQiFdpPBAIBA47kraBrwb+taRi81r+fwzcBrwMuAH4bUl/zsweO+BiTgsTCAQCgWYc8JiZPb9m333A3Wa2C3xa0p+QOfg/OMDyTQlDLoFAINCCmT1B5qz/OoAynpfv/mWy3jmSzpENwfzpCooJBIceCAQCc0h6N/C7wDMl3SfpzcA3Am+W9GHgY8Dr8uS/ATws6ePA+4H/zcweXkW5IUxbDAQCgWND6KEHAoHAMeFIvxSN3Mik0YA5Gj5PLIYBfk82Zqm3bX9W9VTVx656pO+TNtgeBl/bx+pcPmRmV+6nRK985Yvt4Ycf90r7gQ/8yW+Y2YFOYzzSDj2O19lev60zna9TTbXLZLLTnc4SzHZJ0+58k/QCPpXSLPFMZ0h+Nx7J89jNcM51Ho/1aKw++RVEUUSSJN55D5nfKm33OUdDI8mzDvml65vWh+Hr0OS/7rdMDz/0OHff/RNeaePRK87t115fjrRDDwQCgYPFYEU3YR+CQw8EAoE+BIceCAQCxwADDvHMwODQA4FAwBuDyWTVhWjkSDv0wefQp/75rXL+vq/tPmUc+niC7YO33Ydl1KGhOZTn0jjUQy5Heh56mqY4p9Y0mZiOKInq1OKcAyIvu8JhVvymLU8hIrIpVV1ldKW/u+m27fCZNVPYs3ymS2va/Dh8bPdpYNl19KuKPtexz8wRH9uFzaFt+5zzYn9XtSiXbdg61E2fOtTH9rLq0P7IX4r6fFbAEXfouyTpTmtlz5yQo+1QnRNmKYbhNEKotnLMKm6KU4zkGhtPttmQotypN9kuHG9S2dZ0PDY9rmbbfvN3nXNz44F9GmRj+WrKVLetaqfLto+Tqvt923XsSlelzcHUOSBf2+3XsbDdXC7f4+6y4bO9y85hrkODUPTQg0NfDucvPshu0jTXW5gJJ4dwU8c635vJ5ldP0kn+3ZEFK8079aLBpumEzFELS4VTvbM2S0mSIs8IMZrmM7OdOd40nZR+ZwtlnP4mn1MuRJImLT7bSNPFObnV48Hm57O32i7KjLCWoSlJtb3UunM5V+Ia2/MOzeb+b6Jqu9rI6xp85/z7Gpt15aw7pv3Yzupb+7mss+1zHQvbTefTx3bx915tt7GMOjQMhiz1+qyCI+/QAS5efpSLO49iuaPNyJx4Gcmh/LVBkS61lEkyqaQTWDTtgRcBE2XHO8vD4TSfp9lkIaAnsz3GrFzJjCTZXTiecmWeBmtU86sM48zyTFuDiabpGhpsre26fKhvNHV5lhtak/MrflvXKNsY0nZdnj62m5zfzHaz069u87GdBZe1H7fvdSwz1Dnfr+22c7SMOtSbQ9xDP9IvRctMkotcvJSyuX6OzJk3P0qaRaTsklpKUze3cOpZBGda2+OdpXXIYtJ0B2h3pmYxxg5mSWtFz5x4HpvZ5qARhuVO3LP3akbaZZv6BlNvu7vBlhu5l7zCgC/tlmHbt3yFbam7V7oX2511yDNP33TlY+h6ulim7SGvUS/MYDJcZPHQLK2H3rZytqS3SrJcP7jQF/5RSfdK+oikF+7FZsokH1rpfglZjHF7HAmSbw+ju8FKyl/k+lW2qOOlLxS9Zd/8/C+5z3hk17num9803z30Fg+aZRyPL1Hk9wK/b9rDbrvPOV/KePqX8Bj6TwMLwjSSbgT+MvCZ0uZXk63ycRtwB+AnlhAIBAIHjaV+nxWwNIdet3J2zg8B38l8l/J1wM9Yxu8BpyVdu6yyBQKBwN443NMWD3QMXdLrgM+Z2Ycrj6DXA58tfb8v33Z/TR53kPXiOSbvdAOBwFHB6BWAeNAcmEOXtAl8N9lwy54xszuBOwGcRof3zAYCgWNICP0veDpwK1D0zm8APijpxcDngBtLaW/ItwUCgcDhwUAh9B/M7I/M7Cozu8XMbiEbVnmhmT0A3AV8cz7b5SXA42a2MNwSCAQCqyWbTuz1aUHSuqTfl/RhSR+T9I9q0qxJ+vl89t/dkm7pKt0ypy3WrZzdxK8BfwrcC/w/wP+yF5ux1vA9JBE3RnkWmBlODp8HmSyoKaZLtyULQhHqsJ0V0pF4dwb8jjtJ025hkByvedM9VjHqq3MyZLpl0FczZkj6rLI05IpMq7a9ynNeyniIl6KXgZeb2fOA5wOvyjuzZd4MPGpmX0Y2meQHujJd2pCLmb2xY/8tpb8NeMt+7K2PT7MWnyJzqGqNlgRwihCOyKXsJpfrypdrvLhcOiBiklyibr63kQJpHlnqSNPd+nRmQIpzUZavi0lqbAN5mlyaIBVm7eN22XzwCEgbnWw56lXOkTY0tnJYdWvgSg9nXo7wg+bGVghc+UYYlvNvSl+ej2zWHAi1V5ttolzFuSzSDnUDKud7kLbn6pCn7aEc67LrkBcDqS3mPu98/nWUf6qFfR3wD/O/3wP8mCRZy0Edg2kiYnvjGtZHp2EuHD1zrk3MwoMjxtH6wr5ZhLyRpoaciKP5J4AsnVG9Ds6NcG5+8epM/CvBsFJ4ssO5tcUjUqbQWNjOYscXe//1IkYO55p7/3Oh0a0CZM0OYF5iYf63TcEcRcMuGlqX2NJeGmJbIEnZdp/ftgUElcvYLuZmnba7zkeVrhD3LlmCrnIfhO391KGm3++3DnVj2SwXnw+ck3RP6XNH5XgiSR8CHgTea2Z3V4xNZ/9Z1qN7HLiirXRHO/Rf4uTm9URulDm+hd3CuWzxWLl6USlLsx7WyK2zm14uhdDP9yyy32ZOPbVdkmQXOcsXMKmz7YiiMZPJZeSofWIoJD+lNcx2crtuIUx8mr+iLNjeklphpOw7ODmc04L2TDnt1MFIc+N9bT0qSaSWNkqbFj2yup5ROc8kSRYaY/G9znaRZ3G+6m42vrab8i/yqKPo8TXZnl3HYW2Xe/9NTwHlclXzbQqVLz9Z1FGcy7qebm0dqjmmtjpU2N5PHaqzXfyubpgny6+2SP0w+sxyecjMbm/MKtMVeb6k08C/kfQcM1uIrO/Dke6hO41wqnfmUFS+XMu6JZ80d9aOmYZ4291dRBgJ0NwDySpm4XSbeymz38cUl6PpeACsOJ4WbZA0Tach+b49zMXytKdvOj9tPdG92q460Tbbe6XNdtv38u9905Zpc36+mi1dAllt+bbVIR8OWx1qGrMfrrc+zEvRStkeA97PYmT9dPafshd0p4CH2/I60g498hDBLypXm5MsaBmhmcPy8fKuPNM0xWiWKK3m6mc8u/kkHW9L+44hejWgvIw+0qd9NTe6ylns95G6Hdr2rA51O7m++i6+MrLLGD7wsT10Hepju8+5jKJ+C2LsmWIMfZ8vRSVdmffMkbQBfB3wx5VkdwFvyv/+BuA/to2fwzEYcllVfpWRiuZ0eLvqwW33sb6UF0ieDC1c1Se/Vdte5Tlf5Qyh4TlAkbZhIkWvBd6p7IWZA37BzH5V0j8G7jGzu4C3Az8r6V4yGZU3dGV6tB16IBAIHCg21CyXjwAvqNn+vaW/LwF/vU++waEHAoGALwNNW1wWwaEHAoGANwYDB0sNSXDogUAg4IsR1BYDgUDg2BCGXJbFEsJ6B+bw3ssDyyFc8WPNIe+hH+l56F1zsYF8/U6/+cFpWh/Svphn5DUnuchryLUNVdLI6MLMcy6vp96Gyyfq+51Lv+AQ8BNw6nMu6yJR92O7sOln26/35qtxspQ6VNJi6aLPvH7fQKry/20MXYeGwQ71ikVH2qGnNsmCfBoWUs6CRqA4TNeQrqgzQohuFUQBTjGiWXci225ZOtVrdczj1xjK2bTaVlHSet2UvpT1V9oaeXm7j91yGh+tkzbbPnk1pW87ni7Nlr62q1o2PufSJ18f+tjuez6HtN30m6Z86/5uynNfGNlLUZ/PCjjSDh1Szl+8nzTdXXDWMweQ9SwtFXWHWzj91HaRch0UxXN5zFm0CZNkJ3fSmWBWfcNLSdNJ1htKVSvVO9NjWVRnrK2AlpImu/lv6sWRnHMYRai88pD0xaym+deEKTdV/kJSoNW2R1h5+XfVXl2TE6tGTNbZrpa1i2oYepPtrijV6u+quil1dufzGs62z77quaymq+oE1cka7MV2+amkTx2q10ryr0PDYn3EuQ6cI+7QwUh48uL9TNKLC8FiaWrzmiaW9VgLMnnclNR2537n5Ig0ytPMeipJukOmp0MpT1fTAzeSZDKfjtmNosjXLM2dec1x5Q5mmq8lc7az41ts5IWa4wyVPjM1O+vQm2h1UvmNQiz2jNoeu51zrY2trpHXKfZVG/lB267mu1/bhYRE1blVbdU7rhbdn2odqsljL1ouc3Woha6bbPX7kNex2sEbVBv9EDv0I/5StMA4f/FB1kanWYtPApnzrnv4cooysScmpKmRNuiMS8JZnA3rWEqS7jSmw7JY/Ezhslm7RXI4RqS2m91I0npnDkXFzHva1vz4NlXBw6aOtp75xmItFdxn7LLcUy/y9dFE6Qo532vD89F42W8eTfm25V2nUuhj26e8y7Ltey67rmNfga0h69DQkg5TDnlg0ZHvoZfZmTwJQL0rn5Fd7GZnPpfODPB5aefydO0VUhJOtDrzAktTGob969N3jsMrH37yG47wGofuoaGxLP2QIZzfXvA9Ht9z2QefHnJhe+ghiFXaXuU5nxJ66AeHr4NRH9Ws1WmALUF/bHUOuK9yn++siaMgcLWMc7mMtEfB9qrOOUD+wm3YPAfk2Dn0QCAQWCoh9D8QCASOAV+qgUWS3iHpQUkfLW37Z5L+WNJHJP2bQuA93/c2SfdK+qSkVy6rXIFAILB3vnSnLf40i0sqvRd4jpk9F/gT4G0Akp5FJt7+7Pw3/zIXfg8EAoHDxZdipKiZ/TbZKhvlbf/ebDq15PeAG/K/Xwf8nJldNrNPA/cCL15W2QKBQGBPFEMuX4I99C7+DvDv8r+vBz5b2ndfvm0BSXdIukfSPdXFl7Owfb+36obzT2u+6212n04z63Wth64XacuC1VW8pof1EKPqM8/bd3bCKpdRW8bxLMP2oEE1K7bd5zwObXuWcXDoc0j6HmACvKvvb83sTjO73cxuV2lV58itcWL9OqQI59pHa5wTjpjYraGWU2CkIIfTiDhaaysTWfh2TBSNW9PJGWYgjei8ochhJlw0ak+W6120TxGzPDLWsrmQAwgzGf5zffcyH7nrN7557kWTpCtd3+PpMyfbx3Yfx+Zj21c0a5m2u+hrW1qCDIAZNkm9PqvgwB26pG8Bvh74xtIK1p8DbiwluyHf5sU4PsF27sxzGzQNwWsqWpX9H7k1IjfvMDOnZ3O9T+ciRvHGQn5mKVYKKMrCosdUnbXlAUrlpwopxrlFZ71YCYVcx4QkK8TF6hpH4XhnxyOBixbPUV3jqtUjKR1vF9WGWISO+7Dfecx1dnxst9n1PUfVsPtl2fY9xv2eS1/bvuenT7nqtG66zmVR5QYPMCpkM7o+K+BAHbqkVwHfCbzWzC6Udt0FvEHSmqRbgduA3/fJc2v9SjbHV871tGfaGHVDMIvh8U4xo2gdyB2vs5LA1SxPyRFHG1CEvcuyXnyJwrZzI4oniKlejOoeA0VU6oHXVdJpKLPqnXo1FFtzCovZ04MqtotDK9trsr1gD2vdX9cjroZ4V9PVOb+uR+Y6LZA6tceq7ewBpV31r6/taj51tovf7Md2WeCqoNoLbq1DDdSF81fT19mu2tuLbd86NF9/h6lDvflSHUOX9G7gd4FnSrpP0puBHwNOAO+V9CFJPwlgZh8DfgH4OPDrwFusqkRVQ6Qx4yjXbmksx/xYed0jWzYK4XBuRJr3tusqQrEtcmMkozqGv2g7c8CFxECtc7Rs7F1Thcf6dOm0gmQ3KR+tldxCo+00TadDL23CSGW9mHJaX5psVx/z22z7UE3b5IDKmuV76b35lLPJdnmbTySjj4Paq+2hrmO1rHXH1CTE1UY5zUHVIS8OsUNfWmCRmb2xZvPbW9L/E+Cf9LGhqd55Z2k6U6RpIa5Vdp5NuRlpPm7eReKh2ZKpB2YV1m+BBOs87qwSd+dVbihekrNY52HX9caHsO1DWbypLc+6nvNB2/bVJPG1XeCz2MOybPuey668V1mHWjnkof9HWpzLadjiy/sFijUulrGQ5xL0gfxf7Pmfn2WIRy0jbbA9DEO9HF4mh/VchpeigUAgcBwYaAxd0o2S3i/p45I+Junba9K8TNLj+fD0hyR9b1fxgpZLIBAI9GGYIZcJ8FYz+6CkE8AHJL3XzD5eSfefzOzrfTMNDj0QCAR8GWgM3czuB+7P/35S0ifIgimrDr0XYcglEAgE+uA/D/1cEdWef+6oy07SLcALgLtrdn+VpA9L+neSnt1VtNBDDwQCAU8M6KGe8ZCZ3d6WQNI28IvAd5jZE5XdHwRuNrPzkl4D/DJZjE4joYceCAQCvhgwSf0+HSjT//hF4F1m9ksLpsyeMLPz+d+/BowknWvL80g79CQdduWQNPFbXFZyPebb+tnuE+jhL1zlV0afAJe++J4fSV7zpvuQJIn3uVyl7WXMmx6+DnmuLbuE6Y2rrENtWOr3aUPZCXs78Akz+8GGNNfk6ZD0YjJ//XBbvkd6yCULJjWaIkWLNZ5FhJHgnJqDhpRFgJqlSM2BIc6JJBVOEWl3MCtOY1LbyW3X55lds+LT3NgyTZQUS30XRHVI1hg5OLM9O4NtttM07aWuWNykms+lmzqMNgdT3e97AyrK3Lavr+22PKth6F22++Jjuy1isshjLwtmD2W773H3qUNttrPr2Mt0PcW0xf3zUuCbgD+S9KF823cDNwGY2U8C3wB8q6QJcBF4g3VcuCPt0FNLOH/5QbbWrlxQTXQuW92+0DUxE2Zpo2NNbYJZgsyBq4+ynEUB7iBFOERqk5pGn9040jTb5xgDk8YKbZZQKB009XSywIlcbkDK/bkt5FmUpXC8hbhiW2OyNG1102VxLaFp3k22Z8dl0wZZF55epGmjSWGvyzn42O6yX2e70E3pCsFvuo5t27tuHG22i/1tPeoi+KZ6TNU8q9+LPJtsFw6377mss9WnDpXTtNFXpbGTAR6szOx3aNMsydL8GJlcijdHesgFYJJc4MmLn5/2wCFzqNnteKY+KAmZmxOuKv5PbGfOoWbpFpUIjYRJcnn6Xbm0bmbTTf/PnP6klE5gMVJU05uZUJataa54RloeYsp79WVhpmmDnXPPymUN6rI0bA+9rrKiY9V2lTbn09brarJd/l7sb7PdpO3RZbvJ+fWx3YSP7TZNkuq1qZbTN1S+nL7L9l6Ox1dnZa91qG1fVx3aM2ZY6vdZBUfeoQOktssTF+5jkhbOVvlQy6JanKWa9uYNI7HL1A1fRC4imqobGkm6S5Is6rJIQsTMhktSkmRSmw6LIL9RZL2p3Vrbi1Sc+SzT3FGXRIxqB+9qHjd7jKHW9spqVCv32oNqS9/UYOuU94Zg2bbbzpGPU+7qCbfRdnP1KeNex/zbbqC+tpdRh/ZM6vlZAcfCoQMYKecvfiFzbemiMy/IHLADJiTp5do0s7TF6ZnQdoWKG0WaTmqdeRmnmOwBolu0KxtSyf5vMT4d2+5+Sewy/RtPYaauF11C056IT6Pp+8LOV8BpGQ6963h8xv7L+By7r23oN4zgo3NS2PS17YPvC9U+tvvUoWW87MfAEvP6rIJj49ALmhx5TUL/ZB5p+7zlHzztKm33YBm2VykeddzO5TLyG1o0a5XHPeUQ99CP9EvRQCAQOGh6BBYdOMGhBwKBgC/1E+AODcGhBwKBgCc9Q/8PnODQA4FAwBcDj3jClREceiAQCPQg9NAPJcO+Ae8zParPTKrBZ131Kuewxpdhe/BpaT04bucy2PbJlEM9hr60aYuS3iHpQUkfLW07K+m9kj6V/38m3y5JPyrpXkkfkfTCPVjk5Ma1WfSmx1SpiDXG0SZtjt3MsnnjrC9ICyykcyJyI6Jo1Go3tQRDOLfWWUYpwmym8dJgHKcsIjZyi9GtlcRZ4JGn2FLnPHBm84x9zvmQc7bL6YacGlcXhdmE7zzrIm0XXRGTZcor3rfRpoOyH9s++NpeVh1a1s1+CHGuZbHMeeg/Dbyqsu27gPeZ2W3A+/LvAK8m0/m9DbgD+Ik+hpxizmzeQhxtUIhcuajeuTmn3J8JKWLsNmuddVYZZvolo2iDuMZZz9JlUgNOjjga16aTZhNUnSucelMAlCs5lwjnah6mrFisuqQBUrswtOXHXAod7woaqghNLeZYhGx365f0pexUmxp5dfsQtgvNkCYb5e19HL+v7aqNJttNv+lro2m7Tzrfm0lB+7nstt0Xnzq0V8z8PqtgaQ7dzH4beKSy+XXAO/O/3wm8vrT9Zyzj94DTkq71sRO7DU5v3EzkRtOoRTNDNDf64mRn6cXIbeDy0afM8UKtHIDGjOON8jFiliCVo90MyRG58Vw6uRQjneZbqD46N0aa3XwyR16nI6L53n/phlOOtJNUCa5adPrTfNUs2VvtWZXPZblnnpaiG+vyKju/Kl22q9ol1d91RVY2NeS2Rj6U7bZjrrNd7s1Or2OHk/fp/dYJcfnYrlL9nY8QWVsdKv8uTfdeh/Z6HfeMAan8PivgoCNFr7ZsLT2AB4Cr87+vBz5bSndfvm0BSXcoX9IJ4OTGdUiLYdBp7qylco9mVnkquRK7dWKt4Vybqlv2NDCON6fOPLM1X3GLipnJ8RpoQqbH0tT4Y5wb41zUWPEsFxtzLs6f5+odSvE9c+rFzanlUVVCeYPpapjSTGmxzfYs625hpDa7RZ51TrVLR6Ytz6ZGvliHmpUF22jrETbZruZZ1fcuzlPbDaS4WbTXoW7b+7mOPnVoSNtlNcYmfIenfDBEmvh9VsHKXoqamUnq/WBiZncCdwKM4xPWNbZdplELvdhvCUa7VnQhjmQdc5emOiNMED69qfkeRZPtfmQ3gU6Z2dLfbbatrpffgY/EbVdeXY1+L7Z9G7mPk6/7zRDXsZpH2yIO+1Fb9MHnOva1vco6tGfscM9yOege+heKoZT8/wfz7Z8DbiyluyHf1oqLBtaJiPxepGS66n62hV/FzcbdhhHNmtqueXLZb56++IpHmRlRw/uOvRJF0aG3Df1e5g4966fPi2ffczk0q6xDrfa+RF+K1nEX8Kb87zcBv1La/s35bJeXAI+XhmYCgUDgUGBAtlhO92cVLG3IRdK7gZcB5yTdB3wf8P3AL0h6M/Bfgb+RJ/814DXAvcAF4G8vq1yBQCCwZ4x8CcjDydIcupm9sWHXK2rSGvCWZZUlEAgEhmKFsWydfAlHigYCgUB/0uTwLiMRHHogEAh4ssqgIR+CQw8EAgFvVvfC04cj7dAHX1m7R37e08j2WpYhbPeYO7VKEaWhF/Htk98qbR8FYbGjUEZYwkLQrbYOr0M/vINBHqSWUitdUiKTLek+zEzjZZT/pkvnJAK6RcAk4dQu1pXZdhSXwmtuckdkXGHbR1GysJcJjHXbFt22fecPl4NRvOf197DdlraPbd852762C4a2Xff3fvLua3tZ19Ennz6294Wnjsux03I5CJJ0h53JU7T2gy0iYjTVaqmjCGlPbYLTGNGstwHGJL2McJ03iiKaVLntep2KInBjkv+mvmLOfpvkAaA+jW3m1OvynDaaDhGlqdAReeh/y+n20QHpstNGV/j5Mmz3cSx7sd0HX50Vn2PsG0zWWoc60vWxtew6tB8MSFPn9WlD0o2S3i/p45I+Junba9JIPVVoj7RDB3jy0gNc2H14Liy9uKiOCFf0fJU59jpnbZawm1wu/TbOHXbZ0WXObDe5BIU4lbk5Ya2y7STdJUkn+TaHGOX5VE+5kaY7pbIs9vAKVUOzSWEk7wE0aHYUMrkUabTQyBdCtRt0QhYElWh3Cmb1mjWttql/ZO5yfnX7y/k0n8tF2z5OqW1/VWelybavXkw1n66y+R53l+0yVduddaghz1q1zhZxr7o6NKTt/ZKavD4dTIC3mtmzgJcAb5H0rEqa3iq0R96hA1zafZwnL30es3Ta2XS5Uy4jORzFsEq2LbFdJrZbSSes4qxTS9mdXFyw7RThFEM+FGFmJOkO1W5sdgOIYe5GkZAku1QpO5hC+THNbw6VTCmc+rSC13afZzeAckNYaADSVD1xLs8ap1J27F29orqGVret2si7esZVyYK69L62y/iGu5dt7+UYu/Kso85ZN9muph3adtPNu7xtaNtNdurS+V7HXpiw1O/Tmo3Z/Wb2wfzvJ4FPsChI+Dp6qtAe6ZeiZSbpJZ68/HlOrt+II1roSRZIgjQm4SKJJTSNH2QVLMJsgmmXZFLjUKdpHTJI0su0iXZJynQepr3tLoeVkolrNdtGAqPUI29M6P/iSIJcYKp1iAP/BlPXg2xjSP2Qrl7cXvLdi8DWQdsu5znU+VyG7VXWob5kof/eyc8pV4XNudMyccE5JN0CvAC4u7KrSYW2URbl2Dh0yNQSHa7RmRdMewPeV8avt2Zp9+qxmdQnJImvgBO0iOwVmXpOpxFR5EiSlhtE2XbkSJOO3mwxru6Tn3O9Vrvx7SWvaiZGn+Ppk9aHKIpa1Rf3mvaw217lOS/wGE4peMjMbm9LIGkb+EXgO8zsif2W7Vg59EAgEFg2Q81DVzat7heBd5nZL9Uk6a1CeyzG0AOBQOAgMIMkdV6fNpQNE7wd+ISZ/WBDst4qtKGHHggEAj0YaBDnpcA3AX8k6UP5tu8GbgIws59kDyq0waEHAoGAN8OE/pvZ79AR+bcXFdrg0AOBQMATo9dL0QMnOPRAIBDowWEW5zpWL0W34nOs6zSO9vUFHRGnuZHt6IrOPNdHJ9mMrlwIUlrIUxEb47Me2i2WyRG4cadtKSZNfddKzIKbukiS1EvbJosc9Z8O2DVVFPqtvn4UxKN8j6caRToESZJ42x5y2mAf20XaIelTh5YxZTF7KSqvzyo4Jg5dnNv4Mk6v3YxTxJo7Qay12pQxa2zpHJFGbOosZ+IbqHeEYj3aZuzWiaM1Tq5dR9TgrEfRmDhay536KUbxZkM5jUw2wOFcTBxtNB6RczHORVnQktYaylgIj7l8fnuMc+03lCydw7nmG4VzRX7y0owpnLlPuPw03440fZDqtXeq+7vy3ovt4ndttrvK15XHXvMtjqevAJiPcJbvdRxCL2avtpdFirw+q+DIO3RHzLVbz2VrdK7UixYjbbIWbc2lXXcn2NTpWW/SxFibnItumRPvEhHrbhuniDQ1LDWEY3t8DSM376xH8TpSeeRKjOMtNtZOl7ZZ6VOETGeOPXIbLDrriCKycxZuPV7QjcmccvW3IoqanfosXFvUXf5yqHZZ56XaQOoaTZOoWUGapq09rKYw8Trb1WOqS1dnu4k623XOsE53pM120Ttvs90kMtV2Lsth9/s57iYbQ13HtieoJse7jDo0pKZLUFtcEsJx3dbzGUebtdoJkdbYiE4CYis6y5gtqg7QUog04oroFkZsEGnEmtskc6gVx4LYGl3B5vgMIEbRBk7RYqU1iN2YzbWzWSmdq+3oFhUzcuvMLkU03bdwvIqJorX878yZ1wkoZdK+TWJgM9sZnqJNEi6K8j/rG02apnPyuk2NrsnBtA1NFMfRdjxN9srlbBKFqrNdpPWx3UTbsZb/b7I9FZprcH59bZdvFE3XcZW2i+s4dB0aamjO8BPmWtWL05U4dEl/X5lk5EclvVvSuqRbJd2tTCry5yV1DjKP3AZOI5pkTCw1nGK2dSUxa43h8WbZzWHTnWGkuh7zXOkZaZNRtN7qgNLUiFyMU6GUWG982niIoWPsP7vzZz1rqblSz3rgwk8ToF3gqthebuittvfYdpry7GU7x+eR20dnpcnRtuXlg6/tLvGxvrarveb9Hnef8vhex/2Mu7fVoaEIQy4lJF0PfBtwu5k9h8yLvQH4AeCHzOzLgEeBN3fl5fMYlaaW6Y34vOBzfi9RUku9BPjT1DASz5cz/gJF8nhhWSgR+r7Y9HKA+XH4HI/vAgWQ6X102s7z8mnsUVTz1LQP22nP4/bFrw752/alqBu+apZ99Hd8bIOfdO/QdWgowpDLIjGwoWzweZNMPezlwHvy/e8EXt+Zi+eLD+975eGdjbR0hn6J5DPr5TiyzJdxgdVjQGLO67MKWq1KDVNF9oGZfQ7458BnyBz548AHgMdsuoLDVCYyEAgEDg+WLT3s81kFXbeR3wWQ9LNDGZR0hky4/VbgOmALeFWP398h6R5J90ySS0MVKxAIBDopIkUP60vRrkjRsaT/HvhqSf9tdWeD5GMXfwn4tJl9EUDSL5EJ1ZyWFOe99EaZyFwg/k6AzbVzq4sqCQQCX4IIO8TDiV0O/X8GvhE4DfzVyj4D9uLQPwO8RNImcBF4BXAP8H7gG4CfA94E/Moe8g4EAoGlsqrhFB+6HPq1Zvatkv6wbumkvWBmd0t6D/BBsoVS/5Csx/1vgZ+T9H/m294+hL1AIBAYkqPcQ38b8K/JeuqDOHQAM/s+4Psqm/8UeHHPjPyS4TmB5RDfeZfN0JoovsvSLcN2HwY/7hUeS2D5GDBZkU6LD10O/WFJ/x64VdJd1Z1m9trlFMuPLHy+3a875yDxW19Qaey1RqZTRJoYzrXP83bO4RRjmnisNO/yRaHbKeYPdx2Pc9naoT7zjcHPEck5LE0714ks5i772vaZ5+x73NAvMMWnnIVNH9t95osPbduXrM341aE+17HPGrBdtvvWoaFFwNo4zLfsLof+V4AXAj8L/IvlF6cfKQlmKc5FtY7VyTG2MSfSM1zSRS7rUm3lKBzzmrZxFnOBR2m7bE6OjegMl9MncK5eK8O5rCKOo5MkXGDXLtbaLips5NZxljJJL04bXB1Z+hhobuSzgAzl6esr+6yxFD2O5oZWhGMb7c5gLhikZ83vtF2KMqy/jq5TP2Qvtqvfm2z3cbjlfNrOZd+8usqZ76WIXva5jm151e3zvY5edagnbbaHeHgyO8J66Ga2A/yepK8uZqUcLoyL6eOs62SuqTLbI8SmbbGVnMz+TraIXMQFzs/5mczxZqH/Qoy1QWQxT/EwuLIA0uxJwCzrpa+70+xyHrRDWuldGzbtYcRsgYvYSc5T9nJZxQWnEZmqomOkLSbphfwmU69JkfVeHM6BNO/AMgdkmCWlSuyoLpwVRa5yE5xpa9Q1CjObixRtEpSalpFmR9jU6NocTDlisk3QqcvZNh1boTNS3VeOlGzTmWlyQOU8y+eg/H/V0TUd935sl+myXdVZ6Wu7tQ6V8uyqQ+X/i317rUPDhv4fXloduqQfNrPvAN4haeGMrHrIpeBS8gRrbhPHmMyVi1PJGdaZKSMKsZ5uEBHxpHsi91+5LIDNRzZGGrFtV3KRR8HtAs1iQyPbRu4iO8kFZs7aSJP50xVH60gRO8kT03SSA4sqIlGO2G2Rcomsvlvt08e8s571wLNGs+iUzDJNGefalPA018iLdHVdG18NkPIQVpvTLbYv2K5J3/T7Nh2PLttNDsZHobEtX+g+7qab1H5sl8+lr/5N+bc+x+NjuylNk+227WXbfevQkBzmBS66hlyKgKJ/vuyC7JfL6QViUtZ1gtPpWUbUa3uNGHMqPcPj7lFApBVnXuAUsZlewVM8zITLjXYlMWITuYhLyeO0jTNEbsS6TnNp8lg2Zp5GtQ1ZEs7WSewSqe12HHnRuCattosbQPewQLnx9xOlaHwsL1Qhrftl6V7Hidt+5zsM09d20w1nL3n2dTx9bA/t1HyOaa/Xsa2sy7qOfTBgcogH0buGXD6Q//9bkq7M/z6EQy85lnDOrupcXSgiYqQxk6nSQD1ZD3wdc7skHZUkW/zC56VQxHi0zWSy0zr1JutlxETOSJKucvqNOWbDCnFnfpCNL7nIkQ70skkIOU+RNLrGgPunWwZ9xs2HfKkJdL6Y3mvaw257lec843AHFnW+dZH0DyU9BHwS+BNJX5T0vcsvWiAQCBwujCOs5SLpH5CF5X+lmZ01szPAnwdeKunvH0QBA4FA4DBheS+967MKunro3wS80cw+XWwwsz8F/hbwzcssWCAQCBxGhuqhS3qHpAclfbRh/8skPS7pQ/mnc2Sk66XoyMweqm40sy9KncvbBwKBwLHChh1O+Wngx4CfaUnzn8zs630z7HLoO3vcFwgEAseSZKBpi2b225JuGSSznC6H/jxJT9RsF7A+ZEECgUDgsGP0Ciw6J+me0vc79yBy+FWSPgx8HvhfzexjbYm7pi0e3EJ9g+B4Sk+yaSdwLa8HDOMylxDCdSzMfDm9wK5dwjFqXV4stYSUfNX7lhciZsYkuZQHQHQtCp00hu1X8/TBzEjTBD+5MiNNBlzHEvOesgj9jmlV9JkS56PT04e+ejWrsj20xkpfrZxl0COw6CEzu30fpj4I3Gxm5yW9Bvhl4La2H6xqTdHBiRiz7k5y2V3mcfcICfVzrRMSHokfJNEuu7pMovrAHTPjPA9ykUeYcJlElxudxyS9zIXk4VzrxBqDZ8xSdtMnSG2X1HZzh10faZfaDqllo1ptmh6Zo9jtXMsyi5Y0Mmdu+XKsTc6wSKPOfKF7HU3LNWKKdD4aJT62fctX2CvC+1dhG/xt++TZ97osw7ZvGYaoQ1XbQ13HvhQ9dJ/Pvm2ZPWFm5/O/fw0YSTrX9ptj4dDX3BabOjWNSEyV8Jh7hJ1KhOcOOzwcfYFEk6kjTZiQuN05J5xawuO6jwv2WOaezUiZsMuFuZ6WmbGTPsWF5NFpOsh1XCqXNE0nXJo8SmqTUroJcumcUzfLnh+SdGf6HerXw85677tAewRd9jSQVNIV4fBWTohTFjox+53/Su5ZOefTG7M8yunqGmQ1pLyv7TZtkHLapujcujybqAY0NR/PfESnz43Xh3J+PudyCNt9ruOcrk/P67hf2+U8h3bocHDz0CVdo/wAJL2YzF8/3PabrjH0Q89WdJrIxnN9zeyaG0+6x9lim7V0g0u6wBMuc9BlX2sYqSaYS4nSMYnt8qh9BkjmHtkKvY3L9hRjNhAROzzJTnqRxp5urpqYpDvsTJ7MrM35z0xECyWYZZGmiV3M0y3aLjsH56xxOKas35HlkzAf0l9U9kJ+OAvLV64dU240zcp17VofaZpOb5JVLZFCN6UqtlSl65G5TuOlGh3YZLvqGNr0YuqOt0klsCzElaWbP5dNthtlEzpsl7VLqo67eh2rzq0tkrJJA2avtqv0tV13HaMo2ncd2gtDDfJJejfwMrKx9vvI1ogYAZjZT5Kt4PatkiZkq7u9wTrufEfaoUfERKw1C/VgPKXzPBU9yYRJ41BImhrOpVzgMc6nD4Lqx3sLO5fTp0hsB5TSNjZqGLvJBZKk2ekXjTxJL5HapFVlsfg/u7405lk0tJkzr+8BTRt5fgxi/6HSU9u0N+iqg5HUqbE+O672467aqbNd99su2mzXiXsNabuN6o2iSYyrj526c9lah3LahMB8y+J7HZMk8Uo3JMags1ze2LH/x8imNXpzpB16JnDVXknNjEnDOHmZNDV2OI+Rdt6CjRQj8bPd4sxnttO8t12vrFi1bpbS1Taziu3fsHwbYvk3bbZ9NFbKtrtenvmIUfWxXeBz3Hux7XuD6PsU0ma77m/f3+xlf8FhqUPLduaZwdWF9ftwpMfQ5YYdH+sz3uabdhljeMuwPXQ5nfyrlu9CDsvIb5W2l1E3hrZ9FMoIw1/HJg7ypeheONI99EAgEDhoVjhTtpPg0AOBQMAbkR5l+dxlIOm0pPdI+mNJn5D0VZLOSnqvpE/l/59ZRdkCgUCgieylqN9nFaxqDP1HgF83sy8Hngd8Avgu4H1mdhvwvvx7IBAIHCqOrB76MpB0Cvha4O0AZrZjZo8BrwPemSd7J/D6gy5bIBAIdGGen1Wwih76rcAXgX8l6Q8l/ZSkLeBqM7s/T/MAcHXdjyXdIekeSfdkUwIDgUDgYDBCD71KDLwQ+AkzewHwFJXhlTwaqvaUmNmdZna7md3uWPOzaPWh8wsk8ppu12utQvOcfuWteeG3fij4hz4vYw5vki4GfTTZHlrAqS7gZAjbPnn2sd3nnA89zXD4OrT/AKkqdcFa9baHr0ONWDbLxeezClbh0O8D7jOzu/Pv7yFz8F+QdC1A/v+DXRklLdGfkDleh+Mqu4rI4tbKIcSV7iZO6gpciwqic46R1jkRX9Oqqlg0grXRaZyizooZuRHOjVrznGlkjOZsNNmGGMm1H7cEcsi1Kz9WBaF8G1qnbQ/2arst/679TbaGsL0Xca1l2e6uQ922s7T9RcX6XP922/vLoy+HeR76gTt0M3sA+KykZ+abXgF8HLgLeFO+7U3Ar3TmRcIT6QMY6UKQkXOO2GKutmvYZJvr0htYs7WFCytlzvx0cpbt9BTXpLdxFTehmlMjibG22NSVbESnOTO6GafFG0Xx3RERuRFr8VkiN27p/QtwiJimhaBmGhngXIy0Rp2CX1ndT3KkaYTU7EDmexT15ZvT5yj93ebsxN5CwJscSVXAbPE69g8Br6apO56q3bp0dbbbtFl8bZfxFdjysV3dX7UVRf1s+0TG+tr2kU+os19nqynP/RBmudTz94B3SfoI8Hzg/wK+H/g6SZ8C/lL+vZMJOzyafp7EdqZO3cmxwQZXplcR5XrnDsdVyTWcTE9Oe8HOOZxFXJFcxRobQFYZznAdN/IVOKJS5RDrOs26Tk+3xW6ds6NbiLWOy3u4xT7ZvMzoyJ1kFG9ByXZGRFkX3bmYyNU5a5GmVsrTAeO5HvjM+TG1kz2OOqBsY+b053pWym4sVdvT4YHiXNSIMdUJXAktyCOUG1uTE+hq1HW26+jSC6nbVpbarUvX5YjrqB5zk87KMmxX1Qd9zmV2g1+03aV02Ga7rvw+dajud0PUof1wmMfQVxJYZGYfAuqE31+xp/xIeDT5PKeiq4lZ56SdZCvdXhi+EOJ0epaYEQ+5LzKyEaeSK2oXw9jSaW615/MZPsqEXbaic8RaHLN3ijkV3chT6YNcSB/FyWHpYi9BErE2IY64PHksd7pRba9Hiojc2lQPveycq3mmaYyUkq0IWPSKF9NlDl3Abt4IGh6Tpek4YUFa6pmX01WFmRp7hszEuny1Oaa2Gxpik3PYD+Xx4rZe57zt7t6pT9mqY9VNx133xLBf29XruF9xrbo0befSpw71se1bh/ZC9nLv8AYWHatI0QvJY9ym5zFi3JpumxPENsYSWsesx1rnam7jqdFjtK1SIoktdyWXeAJMrS9g42iNlE12J5c7xngdWEwUpa0VsuysM2Gv9jylEWk6aU2X9cBT1NBjmkuX/9nVaIQo5IS76OucB+199XxR2S2mNkvrI0gFfsfv+2K+6Pn7ioD55LlK277Ca13KnPshiHMdIFHHknLldG3OvMDJtb4kLdNnJoLvSzF/gaI+syCGm3XTh2WIhR0V8aijcC6DbT/CkEsgEAgcAxrnUx8SgkMPBAIBXwySQyy3GBx6IBAIeGIc7jH04NADgUCgB4fYnweHHggEAn0IPfQDQvk/H1LLVqXvSp/SHHlXZja32yetVxG90/Zd/NczpV85l2K7n9bIqljOeT9etofWCFrlcc/yXUq2g3Bspi3GjLjNPZsRjrhlrVHD2NEOO+ySKGnVgrnMJZ5yT2SBNS2klnCRh4ncmDhqFgwzM1LL5oDH0Tpt0wfN8puNRa3TJrO5vlmj6Y7eS0lTHxGjWRBQ69SvmujG5hz9Fk4u59V1PL62++Bru+8cZ99Fo31t93GUvnPLfWz75lfgI+7V5zr2qUN1EhH7xYCJ+X1WwbFw6Jts82y9iA1tUvTT45rKYRiX3WV2mGQbZEw0Ia2R0nnKPcmj0RfnHH7dotSJ7XA+/QIJu1kaiajGqZsZpl2SafRnFmRUpxlTFZuUYqJoMVgqSzfBLC2lrTluMyTDTzJocWJWbaOoNNRGnQ9sX4t5NzXIqjRBE3X7fMWeXEO5fW13/bbNdlO+e3VQbUJpdRHITbb3evP0sd0lKrZXBnXqFtQWl8oVuopn6rnEjCj8mqWGBHGpB5ySckEXSUpDKMVJT5SQ5r11w3gifoQn9Vht771c6XbSCzyZPEhqybTXkqZZz7rQY8nspCRcJrVkoYcRRWPi3FlnjheqDjXrrTsiNy5tSzHbAaxku0kjJNN/qxMxmlX2ZtvV46a0rdU2zT2kNgdTjhxsC65qsw3tGi9N++bzrJNlmNdZKZ/TqiJkX6dctl0vCVGv8VJnu0qXwFbXuWyyXb2OTR2Kpn1ddQjabwY+dajJ9l4oukVBbXEJjFnnJm7D1USHFs49RiQkPMVFDGsM1U6VkkS7PBJ9gUu6UOvMixuFmXHJHudi+gh1Y+ZFZYrcGMPYTS+RDXc0NZYo6623VLqskTucxpglmO3O2arazr9B/jTS5Owz+8p7ox7aJaX8vcL9O2wX9puOvUnnus3xlPNrEuEqp2/CV0emzn6Tbd/hirp0XUMtXZHFTTeK9jrUbnu+DjW7k72cS1/b0H4DHX4cf5geuqR3SHpQ0kcb9kvSj0q6V9JHJL2wK88j7dBjRrXCWlUucAloP8lm8CRPsKtdkpYKkKZGwi6Xkydbx98LdpOL0HIjyWxbrq/SrtuSKd6lmE067Wb5JtPfteVJPr7dqU7YZ8w4/9dtu5R/m21fuzY7Fh/bfV+cdZXTR5dkrw6m7Xd9j3tI2+VzOdR1LNO2eEVxk1qW7Vqbnh8Pfhp4Vcv+VwO35Z87gJ/oyvBIO3QXDawT4XwvvCHPcTnfJ72sQnqlHFznxOfFVV+ijgUzygz94iqK/G33Seubn6/zGPq4+7CMOjQ0vudnGS8/G22VOj9dn868zH4beKQlyeuAn7GM3wNOK18EqIljNW0xEAgElk2PxSvOSbqn9P1OM7uzh6nrgc+Wvt+Xb7u/Pnlw6IFAIOCN0Suw6CEzq1v3YWkEhx4IBAK+2IFGin4OuLH0/YZ8WyNHegw9EAgEDhrz/DcAdwHfnM92eQnwuJk1DrdA6KEHAoGANz2HXFqR9G7gZWRj7fcB3weMAMzsJ4FfA14D3AtcAP52V57BoQcCgUAPhtJDN7M3duw34C198lzZkIukSNIfSvrV/Putku7OJ9H/vKT2hUFhYUX5/WI9puia53xe72vvPWvQX0NklUJGiZdmTMbQgR9D65wcFdt9OAriZ4f1XIbQ/3q+HfhE6fsPAD9kZl8GPAq8uSuDxFJQ+5mT4Kp4o/NAnROnOMlYI6IOvY2YNTbd6c61RoVYj0/T5a2dc0RulEWCds6nzWQFuuYHZ5F78fTvNttF+Tpt95i77NQt9lTeN9Q84q4Q+Lp9Q83B9z2ewl6/NWP9bHdFbPpqp/RJu6zrWD5PB2m7jRD6X4OkG4C/AvxU/l3Ay4H35EneCby+K59ddni0QXMFZf7n2vUxt2xu88zN00TTEPeF8rClMdfpDF/B09hio9apO+dwRJyyc1xlT+cs19eKaxXE2mAzPsvJ8fU4RbUVblZZM0ftFDUKSmWV3OHcOk7rzSeGoofgkEbT3y/kWbJjli0e3Slm1TPQpCnoo7ytrRfY1JjbGrmP7TotFl8bbdexGo7e9lvfABRf22X7PqJkXbar+4e07UPfc9nE4IqLAwUWLYNV9dB/GPhOZjeyK4DHbBbTXkygX0DSHZLukXRPYpd5iIe4X/fPKSa6KFNbvGlzja1R1os+EY/4c1tnWFM056wFnNEGZ7SZKSUq4ul2E1dyZk5WwEmMbMzJ5CwRMZI4qau5RpmWTFGJncvUHtfcCWKXqS7Gbo2T4xuJNMYtRFDOQrWzPGIiN5pLUVRIs7LE6pgo2gJmPbLZzarcu3Fk71lqem65RkahI5+mBjbfGMvOr/y7atlqG3ZpeKiukZdt1w2f+crnlimXtc52neOta3xttutuFNWQ+0ymYVHIqsleef8QtqvM6lC97S753OJ3dfu7jru2DtXYrqtD5b+b6lDd3762e5NPW/T5rIIDd+iSvh540Mw+sJffm9mdZna7md0eKXOYT+kpPqvPkjBBgnXnuGlrjXE0f3hjF/HszdOcjkcIcIirohNsu3m5W0lcx9XcyDW43OGs2xabyamFHvmGTnK9nkXMOE8ZseZO4jT/vtkpYju+jrVom/IQzKKvFFI8p6wIuYpjpUI7xcTRNsVlzPJa7GVnWhcRVETMLJ1f4CP7W1SrRW2D8B1+oUFYKV1cXKRusZGmxujbSKsOpkvrZEjbexXnGkKTpE5wq0sLpiv/LoGsPuXpS53t/V7HvZANuZjXZxWsYpbLS4HXSnoNsA6cBH6ETKcgznvpnRPoq+xohwf0eZ4/fhqnRnHz47LE09ZOEKUxscXEar6nXeFOk6bicS4Sp83vaEda51r7Cu7XnxAxbh0m2IjOsZtcZJLu1KaZltPFWUVU2qq1Ijkcm6Q8OdeDr7Od9dIj0jRpXalJKFvLSR2NpqSC1zp0gcCKfNtXiSrSQHeDXeaLMB/bmfKmX36FgFQbPk51Ztv/5bhzrvN4+tj2ya+gSyGyyNPHdtEO+jjowZUWscFmuSyDA++hm9nbzOwGM7sFeAPwH83sG4H3A9+QJ3sT8Ct9806VcmJcPwZdRhInorVWZ16wxhprHePVkPXAY4+XlQCR87uP+r0kXZSMbS2nc57L9M1ebHYY98irr+3hxaOWgXzOzzTtsOXsMy48tO0++Q09fr1K2wVhlosf/zvwDyTdSzam/vYVlycQCAQWCEMuDZjZbwK/mf/9p8CLV1meQCAQaCOLFD28Qy4hUjQQCAR6MJBOy1IIDj0QCAR6sLr43m6CQw8EAgFPslkuh9elHz+HbvTQRfHMb2B6DcENbL+P7aGj3ZahLbNKrZFVauUE2wdvu2BVLzx9OEyzXPbNrWtXkolXtaczg5Ox48RauxaLmREr5rRtEbWcKsOYaIdtzjGifYpjYruAGEcbnbazaVczrZXmtAnSeCGYqS7PNAV1TJs0DCf5Rbu1RPHV2e9Kl9nuVy2HnJbXJ68maYG6PId2LsVcdB/byxAgW4btoeoQLCFCtLDP4Z7lciwceoTjKzefzg3xFYBIySLY665natmagJJYi8SZ9bjWXaZmXLIUJ0ckx1ltM655oDGMnegSu9rJIkrZZsOdqC3nbnKBi5NsTVjnYsbRJnXOOqu0s2CgTAqgKV0mBSSEFBFFawvpirRFd19yqMH5W7EI9bTRNNxQahpWl3Prmitv5A7Sw6dWw/l9ZAL6SgkMNRfeJ0ag0OppKste8u0Vm1ASwjpo223f95LvsmMT/Nz5aoZljrxD39CYl24+g7PRFq6kQ2IGuJlTN6tRQUshEpxZj4lKdWBixqU0zX6T2jTK8ZQ22dRa6ecpl9xTJExKvQExYp0tnaHwTGbGrp1nx56kcKpZL8cxjjbm5ATKjrfcu4ncaE7jpXDmZrPQdjMQjqgkZVDcHMpjN2maZg5bMWXvacx622XbWSMvPc3kPdNys2nT2yinaeqxWbV8c7bbG/28Fk49aZq29habdFbqbJft+PRWq/IDddIMvrar+frYbsujarvpOla/d/W863rnXcez3zrUZNv3JuWHrzsPPfTeRIr46o3b2HDjhVfPZmAIE6QppMq+L5xny07C6bWYkWCCsdNQYYTY1jqnoy0SJlzUeVItOgpLIdKIbc4iEzv2GLvphVp9DRCjaAMRZYWhOQRbcsTRGLMUmACLlbto5E5rmDU/ek63KUKKGh1qcS4L2YByz7zOdpmq86tDUoft+Ube5Py6nEed3aYbw7JtV9PUOb/qTaopz+rvmgSuytt8bPs8LbTdpJqcblWIy8d23x5503UcAgMmJF6fVXCkX4puaUysqPHFoaVG4cNlXXoj2TDNbtLxqGTgMC7pqcxsQ/JMfMpxcfIoqdvt7lXgOgMWit52IUrZdDzpdPDbT5/Dp7pXG1pjT7v2ptVi26OxlW23jct2KSjupZy+tn3yatrflKfvOV+W7f3k25Snz/Esy/YwGKbDO8vlSPfQXeRR/NKQSxepul+oAkws9X6ES9XtzDP8ehG+L6Qgk9P1cwJ+x1Os0uTTgHwFnAwjitpfTpdt+uQZRZFXGc1Wa9vrnPew7Yvvi0Xo92LT96Vm+f82+oiA+VzHITjsL0WPdA89EAgEDpp0RS88fQgOPRAIBLyxlc1g8SE49EAgEPDEyGS6DyvBoQcCgYA3xoRJd7IVERx6IBAIeJJNLD68PfQjPcslEAgEDhYjJfH6dCHpVZI+KeleSd9Vs/9bJH1R0ofyz//QleeR7qEPrtXQS/jHN+HeijKE7eMmxBVsD8cqbftyWM/lED10SRHw48DXAfcBfyDpLjP7eCXpz5vZ3/XN90j30NM8cKgN50Qk6JryK8G6/Fa8HBHjTESd2iBiTae8xKac/ObRZutY+mmXYN1RdnmufnOTO6IWy7b7LuTrP7e+23afOdvFosc+Nn0iFvvY9pmL7mt7L/heR1/Tfea3r7IO7QfDsuhwj08HLwbuNbM/NbMd4OeA1+23fEfaoV9OUy6laaNTd4IrxsbzT6dcu5F9r0MuW+H+mq2Ip58YEynb1sQaMU/X9ayz1ujUJSEcZ6KbOOmum9NrmStjHtocuXXW4pN0LQwtiVG8hXNx4yLFs3DxEW0PYTNn4abpumyDgxaBJN+Aq2n0Zeniddtup84J1P2uuq3NsVZFwNps+2zz/W2Xverv+wqlddmYlxBot+2zrS7foerQXn67H3oMuZyTdE/pc0cpm+uBz5a+35dvq/LfSfqIpPdIurGrbEfaoQN8+skLPDGZLIxsOODGjZRbNlMiB9evpzx9O11w6hLEDk7GEElsjRzPOLXGmhNRSbFLuQjgWI7YOSJF3GRXc5oTuJo7inBYKpwc2+4qroiehqgq6uV/W5Trr8SsRScR0YLeRpErKL9ZrBFFY5ruZkVP0SnGaY1qr77q/CRhFi2kW9AGkUgTy3RxqsfcoM9R2Kv7OzuqmbxsUyPviphs0lmp2otqoovbbJfFoPrart4o6pxNV48+rYnOrdNZGdp209NGWx0ql6O1DrEoWOZTtq7j8RH32i/ZAhe7Xh/gITO7vfS5s6e5/w+4xcyeC7wXeGfXDw7coUu6UdL7JX1c0sckfXu+/ayk90r6VP7/Gd88P/fUJb5w6fJUCyUSPONEylUVafIzY3jWyZSRm6nDrsdiy80/Ko6duO3EiBOxw+W9dSHGcrhKD+NqznINV8z3NImQublt6+4kV8dfQcR42qMSEZa6Sq/FMXLbRG6tIiClhd6N05g4WtRfX+ylusypm+Z6c3W92TTNe+G0NAxlN5Zq9elSwmu1TeZAinPWJzS9zXa5kWcOqN7plx1MX9vF77ttN+vF1NHm8MoObT+291KmtutY/e1enOsqbXdjQ8nnfg4o97hvyLfNLJk9bGaX868/BbyoK9NV9NAnwFvN7FnAS4C3SHoW8F3A+8zsNuB9+XdvHt3Z5QuXLnEyNp59MuXEqD7dZgx/7lTK2MFWDOuu4dFN4qbNmFMjR2QwsuZhhtPuBDdzNVm/PG4cXom1xlXRlyOLwYQ15CmJWBs4xpgZWce46fF0xCjeqj/YSp4wwtJ2x5v11DOHndlued7WolNvopCvbbVd6qn7jp/6CjMVdocQuKqm6+ppFraHenHn89SwbNtdN+/yk8VB2+5bh3qXg8Tr08EfALdJulXSGHgDcFc5gaRrS19fC3yiK9MDd+hmdr+ZfTD/+0myQl5P9kKgeKR4J/D6vnnvpilPP2F0LERE7ODMmhhH7XdwSZwYRYzjqPNuv641ImK6Xqs6RaxH28jjJWjkRkRR7PEyzveFqohiXxEj5yd41KMX1Gd8cy+95IOmz/EMPbbbR4xqaOGqVdpe5TnPMFLPf625ZJKpfxf4DTIf+Atm9jFJ/1jSa/Nk35aPYnwY+DbgW7pKt9Jpi5JuAV4A3A1cbWb357seAK5u+M0dwB0Aazp5AKUMBAKBjCEDi8zs14Bfq2z73tLfbwPe1ifPlTl0SdvALwLfYWZPVF6cmKTa56X8xcKdACfjaw7/ZNpAIHCMsOKF56FkJQ5d0ojMmb/LzH4p3/wFSdea2f352NGDqyhbIBAINGM+4+MrYxWzXAS8HfiEmf1gadddwJvyv98E/MpBly0QCATaMLIFYXw+q2AVPfSXAt8E/JGkD+Xbvhv4fuAXJL0Z+K/A31hB2QKBQKAFCwtclDGz36E5YP8VB1mWQCAQ6IWBWRhyORCeed0Of/ElD7Cx3q5XHLuUv3j9F7nlxIWOHI1rN4xnn+yUjCF24tnrV7GujnukwRV2E6d0ZUeOsOa22XLX0G0dRvEJv+mLFuHUMEm/gr8siedc9DzS1AtvAbLVvRfvqxkzJEni71T6pD3stld5zjMGCyxaCkdabbHMX3nuBe742vOMIuPqc5/jN3/3ah56ZDGKcmu0y4uufZhRlHLTiUt86OFdfv+Bk1SdpjCu3TS2osy3nFkTv/egsVNzndYicWocI404Ha3zsUsP8lhyaSGdTGyzjXOOrfQUG+4UDyT31h7PZnyGmHUkMY42eOTSn5HWCutnkgEom7e+s3ue1HZq83QuQsrmrTtiJpOLtekKzZbif2t5qy9UCjBqHjt0zk19tJwjbWjoksBm/1uHZ88CodrTlCMlfQS0fPIsp2tLPxO48gjU8kS5lk4hLNZ0PMu0XeTbHqBmOKfG6Ny+ZPPKZ4F2XbZ9r2MfDCM9xLNcjnwPPXbGW1/5JHd87XnWYsMJxqOUl7/0AZ5+8xNzac9tXOTF132RcZTgMGJnPP/ck3z90x4idrPKEcu4eTtlO7ZpgPupMbzsWnGy0rk9uRZxai2adjwjOZ6zcTU3jk/NpYssYttO4HBYmsWUntbV3BI9D8esZy0c26MrGbkNsmhNiLTGFetPZ6SNuTydizJnXvr1eLTNeLQYOSq5vBFkDcIpIo43F9JlodXZ30WjyCYlVW94qgmicjhX00fI7iJzTkUdok7TdB1PJ9UQ90XTszDxttD+rpDyNtvVstfZblMD3EsATNmZ7cf2XrZ3ySSUz2Wa9neobecjq7/ttguWFil6iF+KHmmHHgl++I2P8heecYm1eP7ixZHxouc+wotf8DBSytPOPslzr36U2M3LSsVKuX7rMn/ztgc4MZqw7oybtxLGjnmZOYM1B19ztbhuI3NvZzdiNuI6tyZuHp/mWZtXIRNrNmbTtvLQ9lLCVGzoBE9zL2LMBo6Y7fgqIo0p1wczQ0ScXruZjehMfnwjRN0QS6bcuD4+nX2TKFQZy04gTQ0nRxxtUtVuKTfCWeOJp0M6kWuOnM3UIkflDTW2cwdTo7xX1wiFphLETY29rZFXbS+WucV26UbRJRxW3V/tRXYJTrXZrrNX5Fn0moe03bS/9jpW7FbT1QmL1SlEDmV72cNwh9mhH+khl5uvmHDruQlRfQwScWTceuN5Tl2+RHJRRA2P7xEpJ8fGK65/hI89fDZz0A2VIhI854xI5cAJa+iBOMS5eJMborM8njSP6VsKscZc576cR9wDQHOewrEdX8VOmo39N1UaM3AuJoo2SNPLtWmgcNwicmOS9FK+ra0iOsSi+t98nrm+uLLH42bb83a6hkJ8NVuqQlR1v6sXJWsWhPKRF6hzML5juD6224YPkiTZs+020a6Z7WYJ3To7vtexsw510LcODYGFWS7LI4rU6MynaVwKl6y2L1tGGI9eXvda4GInzaRvk67HyRTOJ0n3Cz6Dy7qEz8vPJH/D3i27amRyEd2YpX7jjR0StjPbKU6aql924bMoRTGW7mW7xyIXPmnrJGxry5irAfra9nFAfWz3xcd2n3HoPmmHtt13YZP9sKretw9H2qF7ed9e+fln6JvUN8c+h9LWY9pLuiztsC+Qipea3mkHpE9+q7a9qlk6q7TtS9+byYFgdqinLR5thx4IBAIHipF6PvmuguDQA4FAwJMi9P+wEhx6IBAI9GBVQUM+BIceCAQC3ljooQcCgcBxITj0ZbHKl/SHe4LAygmnJ3AsMf/pwKvgSEeKek1pknAjIOo+1I14l8jjFrcW4z3HeiSHc91TqpzHeqQAcRR52zarjzSs4huQoY6IyTJpmviFtEteAk5FpKiP7bpAm3rTfrZn0Y3D2u4zb3rI9TH7TAcs5tb7pvWxXf6/Dd8AIxheBKyJYgm6wyrOdaQd+lNPxexOhDU5wiiCU1uc/aevJL7pNIxbvHUsnv86x7P+G4haxAidg2s3d/iGWx5mLUppqm9ywgledHqbU6OYqMMZbNlJTnK21alLwmnEmdEtOEWNld05hxDro7NEbs0jxHsNp0Uhs4U8JYqHuqY8Zw1VwGJIenP6dttt3+vyqguHr6atpu+ynUXg7r/JVG23Xce27015L8P2UNfRx3ZTSP9+bQ+DHerQ/yPt0Hd3Hf/ht67j4sWIxCqHMorh5mvQ//R6oluu5Ow//cusf9VNWfe6jBNad1zxV8+y9YxNnv91KV/zN9PMqaua1Lhp+ym+4tRj3Li1yzc//YucHk+IK2GociIWnBnDRhzxopPb3LCxRtRSUYXYSk9xNr06k72qqaSSg9Qxjja4Yvx0Yq3h3LzxouJGbkzkIsbRKUbR1sKNYpa/A0QUrRFFW8Birz5Lmwl7yTmMuDZd8T2L6ixEnBYOYw6rSOo2OZJMz2bWs6yen+rvygJSTXarPeQ22+X/u2z7yhSU03c57zZFxSbbTY7Ox3adXkyV8rbqMTedy2rk616Pu4yv7eFIPD8Hz5F26ADnnxrxG++/jocfWSNJ88OJI/jzz0Lf+HVoLetuaxRx8u+9hO1vej6McycYO+IzMVd+wxWMr5x1y296Frz6W1PWt8DFAmXO/CtOP86NW09N/c/JccrfetoXuWX7MrErGjusR3AqNlzp8fIZW5t8+YkNio569vi+2PDX2OTK9AYii4hKFVM4ZLMejlPM6fgWNqKTU/GtrAfvcIzmHm1H0TZro9MUd6hZhY+mvy3yjKNtqs66kD8t55maA0XTPKcNdq5notKnWWxqeuKKcrQ0auU3ivINanojaXGkzrnWXlqds64bFqneKKrOr8l2m/NrulHU2a7m62O7zlaT7SaBq66bX5XqjaLuBlqnxdJWxiJNd2+72/beCT30pTOZOH7rP1/Ff/ncKSyK0F/7WtzLX1T76Lb16mdy5v/4i2jsWL91jXOvP0u0uaj0cuYaeO13pJw+k7DmEl5w7hHOri/qjI8cvPaGR3jxufMAbEWwHdX3Cq5bX+f20ycRhZxr/fHEjDiX3kCUjoFsfN3VLF4hie3oWk7EV2fpFIPVKyHG0Tob4yuAQp+6Xu1OcjhtlZQcm54WhJkoV6EkbeqVKLfb0bBqHGiT+mJV5MlHG8RXN6XpexM+OitdafZmuz1fXw308k3K1/l15d3Hdtv3unx9dIyWS+r5OXiO9iyXOcSn7r+KL/uh16DtcWvK8bOv5qpvvQV3sX3ForVN+JpXXOTh37mIS5ovkAS3n3uKex8/gXU83p0axZxeG/Ho5XaRfIdjy05BbCQdFXgjOsOF5JHOR8vIjRiPttjZvdg5vpzaiCjqbjxyEWk68Wi8xcIIHRW9x8tCIazDqfWlT16+LzZ9F1voY9tXjKqvYNgqbfdZVKSLvuJe/hxuLZdj0UOfI/JYhg1Q5DeuJoHH5IZpWr90nrY9Xir2zrPjZWHfdH1YhmjWgYky7dP2UTiXwbYvoYceCAQCx4NDrFIZHHogEAh4Y53r3K4SHXZN5DYkPQl8ctXlKHEOeGjVhShxmMpzmMoCoTxtHKaywHDludnMrtxPBpJ+PS+PDw+Z2av2Y68vR92h32Nmt6+6HAWhPM0cprJAKE8bh6kscPjKc5g5fi9FA4FA4EuU4NADgUDgmHDUHfqdqy5AhVCeZg5TWSCUp43DVBY4fOU5tBzpMfRAIBAIzDjqPfRAIBAI5ASHHggEAseEI+vQJb1K0icl3Svpuw7Y9o2S3i/p45I+Junb8+1nJb1X0qfy/88ccLkiSX8o6Vfz77dKujs/Rz8vqV3kZtiynJb0Hkl/LOkTkr5qVedH0t/Pr9NHJb1b0vpBnhtJ75D0oKSPlrbVngtl/Ghero9IeuEBleef5dfqI5L+jaTTpX1vy8vzSUmvPIjylPa9VZJJOpd/X/r5OcocSYcuKQJ+HHg18CzgjZKedYBFmABvNbNnAS8B3pLb/y7gfWZ2G/C+/PtB8u3AJ0rffwD4ITP7MuBR4M0HWJYfAX7dzL4ceF5ergM/P5KuB74NuN3MngNEwBs42HPz00A1wKTpXLwauC3/3AH8xAGV573Ac8zsucCfAG8DyOv1G4Bn57/5l3n7W3Z5kHQj8JeBz5Q2H8T5OboUMpdH6QN8FfAbpe9vA962wvL8CvB1ZFGr1+bbrgU+eYBluIHMMbwc+FUyzdqHgLjunC25LKeAT5O/dC9tP/DzA1wPfBY4SyZ18avAKw/63AC3AB/tOhfA/w28sS7dMstT2ffXgHflf8+1LeA3gK86iPIA7yHrDPwZcO4gz89R/RzJHjqzRlpwX77twJF0C/AC4G7gajO7P9/1AHD1ARblh4HvZCbzdgXwmM1WtD3Ic3Qr8EXgX+VDQD8laYsVnB8z+xzwz8l6efcDjwMfYHXnpqDpXByGuv13gH+3yvJIeh3wOTP7cGXXYTg/h5aj6tAPBZK2gV8EvsPMnijvs6z7cCBzQiV9PfCgmX3gIOx5EAMvBH7CzF4APEVleOWgzk8+Nv06spvMdcAWNY/3q+Qg60oXkr6HbEjxXSsswybw3cD3rqoMR5Wj6tA/B9xY+n5Dvu3AkDQic+bvMrNfyjd/QdK1+f5rgQcPqDgvBV4r6c+AnyMbdvkR4LSkQlHzIM/RfcB9ZnZ3/v09ZA5+FefnLwGfNrMvmtku8Etk52tV56ag6VysrG5L+hbg64FvzG8yqyrP08luwB/O6/QNwAclXbOi8hwZjqpD/wPgtnymwpjspc1dB2VckoC3A58wsx8s7boLeFP+95vIxtaXjpm9zcxuMLNbyM7FfzSzbwTeD3zDCsrzAPBZSc/MN70C+DirOT+fAV4iaTO/bkVZVnJuSjSdi7uAb85nc7wEeLw0NLM0JL2KbMjutWZWXsrrLuANktYk3Ur2MvL3l1kWM/sjM7vKzG7J6/R9wAvzerWS83NkWPUg/l4/wGvI3sb/F+B7Dtj215A9In8E+FD+eQ3ZuPX7gE8B/wE4u4Lz8jLgV/O/n0bW+O4F/jWwdoDleD5wT36Ofhk4s6rzA/wj4I+BjwI/C6wd5LkB3k02fr9L5pze3HQuyF5m/3her/+IbHbOQZTnXrKx6aI+/2Qp/ffk5fkk8OqDKE9l/58xeym69PNzlD8h9D8QCASOCUd1yCUQCAQCFYJDDwQCgWNCcOiBQCBwTAgOPRAIBI4JwaEHAoHAMSHuThIIrB5JCdk0tYKfM7PvX1V5AoHDSJi2GDgSSDpvZtsdaSIzS5q++/4uEDiqhCGXwJFG0p9J+gFJHwT+es33N0r6o1wL/QdKvzsv6V9I+jCZ2mIgcOQJDj1wVNiQ9KHS52+W9j1sZi80s58rfwd+m0z3/OVkkatfKen1eZot4G4ze56Z/c4BHUMgsFTCGHrgqHDRzJ7fsO/nG75/JfCbZvZFAEnvAr6WTIogIRNXCwSODaGHHjgOPNXxvY5LYdw8cNwIDj1wnPl94C9IOpcvm/ZG4LdWXKZAYGmEIZfAUWFD0odK33/dzFrXJDWz+5UtIP5+MpW+f2tmBy2TGwgcGGHaYiAQCBwTwpBLIBAIHBOCQw8EAoFjQnDogUAgcEwIDj0QCASOCcGhBwKBwDEhOPRAIBA4JgSHHggEAseE/x/9IY18krI9ogAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"yfilter = diffs < 1000\n",
"betterdiffs = []\n",
"# plt.scatter(errors[yfilter],diffs[yfilter],s=1)\n",
"# plt.rcParams.update({'font.size': 20})\n",
"# plt.xlabel(\"Error to the true value\" )\n",
"# plt.ylabel(\"Least Squares Difference\")\n",
"# plt.show()\n",
"\n",
"bins = [0,21,32,48,np.inf]\n",
"# for g in range(3,-1,-1):\n",
"# i = np.where((diffs[yfilter] < bins[g+1])) and np.where(diffs[yfilter] >= bins[g])\n",
"# print(g)\n",
"# plt.scatter(errors[yfilter][i],diffs[yfilter][i],s=1,label=g,color = ['red','yellow','green','blue','purple'][int(g)])\n",
"# plt.show()\n",
"# plt.hexbin(errors[yfilter],diffs[yfilter] + 1, mincnt=5)\n",
"# plt.xlabel(\"Error\")\n",
"# plt.ylabel(\"Certainty\")\n",
"# plt.colorbar()\n",
"# plt.show()\n",
"plt.hist(diffs[yfilter], bins = 50, log = True)\n",
"plt.show()\n",
"plt.hexbin(errors[yfilter],diffs[yfilter]+1, mincnt=10,cmap='magma',bins = \"log\")\n",
"plt.xlim(0,150)\n",
"plt.ylim(1,150)\n",
"plt.xlabel(\"Error\")\n",
"plt.ylabel(\"Diff\")\n",
"plt.colorbar()\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.8.10 64-bit",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.10"
},
"orig_nbformat": 4,
"vscode": {
"interpreter": {
"hash": "916dbcbb3f70747c44a77c7bcd40155683ae19c65e1c03b4aa3499c5328201f1"
}
}
},
"nbformat": 4,
"nbformat_minor": 2
}
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