{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " ___ ____ __ __ \n", "/ __)( _ \\( \\/ ) \n", "\\__ \\ )___/ ) ( Statistical Parametric Mapping \n", "(___/(__) (_/\\/\\_) SPM12 - https://www.fil.ion.ucl.ac.uk/spm/\n", "\n", "SPM12, version 7487\n", "Octave, version 4.2.2\n" ] } ], "source": [ "%-Welcome\n", "spm('asciiwelcome');\n", "[vspm,rspm] = spm('Ver');\n", "tlkt = ver(spm_check_version);\n", "fprintf('%s, version %s\\n',vspm,rspm);\n", "fprintf('%s, version %s\\n',tlkt.Name,version);" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%-Working directory\n", "wd = fullfile(getenv('HOME'),'data');\n", "cd(wd);" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "%-Input data\n", "ref = spm_vol(fullfile(wd,'T1w.nii'));\n", "src = spm_vol(fullfile(wd,'bold.nii'));" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAJMmlDQ1BkZWZhdWx0X3JnYi5pY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpTNDAsAAAAJcEhZcwAACxMAAAsTAQCanBgAAAAddEVYdFNvZnR3YXJlAEdQTCBHaG9zdHNjcmlwdCA5LjI2WJButwAAIABJREFUeJztvX+IXFWa///M51NCBdqkAp2lB1IhWRJN6BYcMGCiAR1HSO+uCwac6ICzCUzCOjAjzk5mYXYyLIYZWVtGMrOMYMM6jH9oK6ggTARdFNTuQFwykArRieyEVEOatcHqSUH3koJ8/3h/6/m865x7b/3ojrlV9X79UVTduvfcc2+dW+d5n+c5z/na9evXTQghhLjZ/J+bXQEhhBDCTB2SEEKInKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRCwo3uwI3ittvv/1mV0H0AZ999tnNrsLaoAYvOiHnDX5gOyQz+9Of/nSzqyByzW233Xazq7CWqMGLbPLf4DVkJ4QQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5YIg6pEKhUCjckDD3G1dyt+SnJkII0S1D1CEJIYTIM0NkTTcajfyUXCgU1rY+EEY37hqFEOJGI4UkhBAiFwyRQsoVay5l2hYoCSWEyDlSSEIIIXKBFFI/sRqVI20khMg5UkhCCCFygRRSPzlX5CgSQgww6pAG6u97kK5FCDFsaMhOCCFELpBC+n9w0h1JDSGE+IqRQhJCCJELBl8hdZ6kp2dVtPpQAgUjCCGEFJIQQohcMPgKaa2ES4aIWb2yuaHaSPJLCNEXSCEJIYTIBYOvkNaK1Sut1ZfT9hRtBVxXaknSSgjxVSKFJIQQIhcMhUJKs/TbKoDEr+Kjgi3xDl9BAN6NjiQUQogbjRSSEEKIXDAUCmkNZUHirKZgSyenC8pJlEE33YUjOSWE+CoZig4JrMlIWsYhXn4nHUkPfVjn3PSeTAghekBDdkIIIXJBn3VIlUqlWq3yx6WlpexDPOS60Whw3IGzVnXj8rPrg5Py2RNrklZg22rjwLW9utVzI+75wNNDgxeif+mbv4alpaXDhw+Pj49Xq9Xx8fFjx44dOnSoXC6fP3/+xz/+8d69e292BYVYS9TgxRDSNx3Sa6+9tmfPnmPHjpnZAw88MDExUS6XT5w4Ua1Wjx8/nvF8tnXqrLnHJTtawd/3fMab6xzq+XbJp9UVPTd4IfqXvumQ9u/fjzcYsqhUKhMTE2ZWLpfn5uZuZs2EuAGowYshpG98SOVyuVwuz87OHj58+OjRo9iCr/bs2ZN9bJozpkF0WI1sr0+wJ2K7O3HnxHVYvaOlq+vi83ZbVXEjWE2DF6JP6ZsOycympqZmZmZOnjx58OBBM3NnrwxGMZCowYtho2+G7GZmZubn50+ePImPExMTlUrFzODy7aSE7AynPlM120eS6AEqFotmtrKykrinv+/Q+/KVzSKKF9ewLqsqbhyrb/BC9B190yEh/vXQoUP4+Lvf/W56enpqampubu7IkSM3tWpCrD1q8GII+dr169dvdh16Z3Z2FkPt8Ve33377n/70J97S4Wp7PSzQADqUPm13vkHqJFEPxVuGShXddtttn3322c2uRRd01eCFCMh/g+8bhZSIgl/FUKEGLwab/u6QuiI7i10n2U7TyvTMC2mHdCV6bpBGyb78oRJGQoh8MkQdUjaJ0QprtRaR/u57RuEVQgwP/RT2LYQQYoAZCoXUdo3XTg7x7dkrzybuIHpGC3MIMTxIIQkhhMgFg6+QghmvoOeoa7mLuiJviiQ/NRFCxEghCSGEyAWDr5Aylp9I3N6VER3nIS0Wi0EOoRvE5s2bzWzbtm1+XjMbGRkJalWr1fAV4K/q9bqZlUolM5ufn8eFI2HawsLCmlRSikQI0TlSSEIIIXLB4Cukzme8piXUaTvXtasEQqvhwQcf9NNB9zQaDUgclmXYoVarmVmpVILcgaLC+1hIoRAz2759uxeOnefn583s3LlzN+66hBDCpJCEEELkhMFXSIkuogzvUVoeIN8S+42ys/L0zM6dO81s165dXixED0TMjh07zOzChQvQRrHjauvWrdj+0EMPmdnFixfNDO+vXLliZpcuXfKzLCwsuKKypnsJZUJa3XnnnWY2Pz///vvvr9UFCiEEM/gdUiIZ3UZ2wHcQRH4jBuvQkezatQsFLi4umtl9993nO2AMDR3DhQsXOACB+1FfpSnuXOORxlKphHE89H849t577/WdP/nkE2x/5JFHzOzzzz83s7Nnz67JJcdofrEQQ4iG7IQQQuSCIVVIiUN2Xa3oitc4lGA1dj1iFqBO6vU6huYgXABkyocffshHeXSDv2LYzd9zeDfexwfWajVIMQzc4aT8HiN7hULh9OnT1hzBW15etubQX2Kwe8/3pMOoEyHEICGFJIQQIhcMkULqSgClRYoHq1TEAQ4dVsNLuOOOO6ypP6BIoFQmJiagZvD69ttv+3s+6ebNm1HON77xDTN79913rVXZLC4ujo2Nmdljjz1mTTWDLdgZwRH1eh1RDzg7ToT32BkCaPv27T6R1k96yy23mNmf//xna1VgttY6RqpIiMFGCkkIIUQuGCKFlJY6KPHbWBW1FUPu++nkFMViEeHUExMT1hQuiK+DIllYWMDGCxcu+LHYAeoEu/30pz998803vVioIo4RLxQKCMNDXBwq+dFHH1lz3ite77jjDkglCKmXX37ZrwU7oLR6vY4S4N/64x//6JeAU589e/YrXo42b/lbhRA9I4UkhBAiFwyRQooJHEXZ6YLarljh2ihbXUFhjI6OQpFAedx9993WlD7Y4qF0o6Oj1gzAe/jhh83s29/+tpnt2bPHzF544QWkQ8WJsDNqgveNRgNvkPuH3UsAzqEzZ86gYvAkec5WazqH4E86c+YMnF7QW9BJrNi2bt2Kr7IlS7fTjDJkkLSREAODFJIQQohcMPgKKdEYj4PorNUMz16xwpKkT+BkikvgFSL27NkDsYL4Oo5tgzZaWVlBPNv999/vx77yyivWDMxD3N3Y2BjUDIse7AxJ5MthYAfoGE5HhFN/+eWX69atM7OrV6+a2ZkzZ6wpwvCK0hYXFyGh+Fo4Vesdd9wBnxNPSwoCFOMb29YPJBkkxDAghSSEECIXDL5CamuMJyapa7uEeYfLT/BuUBKTk5P4CtoIDhjsMDc354eUSiWoEGyEjvGpRWa2e/duMxsdHYWQ4h2gZiCD/KtYqwGXVgjhw9khyxBrh2hAFPLUU089//zz1hRAp06dMrN9+/b5DqOjo5icBI3FE6fa5llQIgYhhpzB75A6JJjxygRdVFdxxvxfjJAE9B+1Wg1R3chPim4JXyEGAaEB1uwneEQOoQTYbevWrRgiwynQDSBcAqV5R4vYh1tvvdWa3Ql6xEqlYmb1ep1H9vAeeVTRI2Ksb2pqimvC/Sg62pWVFR+TtObcW06zlBg933bFqYxvFfYtxMCgITshhBC5oP8UUrVaXb9+/YYNG8ysUqmUy2W8TyMxqCEjfiFtAqzvlrZsUoaRziu9Yrdt27Y9++yz1tQNGF778ssvrTnMdccdd2BgDRt9LQkz+/rXv25NyVIoFHzKqjW1ESQUdjt37hxLE06xChEGndRoNHwirbUOJ6JMBH/v27ePl6PlSbsIsnj44Yc5shxDf9BtbYM+gi3Zoqe3xLj9SLcNXoj+pc8U0tLS0oEDB/BAHjp0aGZm5vDhw7Ozsze7XkLcENTgxVDRZwrpueeeW79+/dLS0uzsbLlcPnHiRLVaPX78+N69e9MOSVwWtisSY8RjVeRfsZqBNnI14zv/x3/8B2QKdsYSruVy2cvcunUrjoKmwcITANvxOj8/z54naCMEcEOXWFMPIQWqz5b1GsKrVCwWUT2Od8CBLumwkX1ICGfAiVwnIaiBBRn8ScFSs2uVYagrZ17fSageGnwiHWYBzpjGAILMWEKsLf2kkKanp7ds2YIhi0qlghRq5XIZTnUhBgw1eDFs9I1CqlQqlUrl5MmTH3/8MbZAT1jTAO+QOLY4CN3OdjhlxIP5brAiEXWG/RGiDYngK0TwTFUciy0Is/aJsdAx2A2lYYdNmzYZrVCOu3Ht2jXfDWKoVqtB9OA1LbR6+/btn376qW9kvxQviL5161a+G/BCQRIhXHBxcRFB5w899JA1k8biQhCz56lgOwmyT6SHKEf3/PURa9Xgu6Ltj8IC+gbdUn6a+LnYuHGjkUBPm229vLzM9Y8naPsWvpbYJYzdMJlBfGX0TYc0PT29efPmqakpDFl87WtfQxyzNcOO08joaYTILT03eCH6l77pkI4cOfKXv/zFzObm5iYnJ69cufLf//3fZlatVsfHxzMO9H6IXTsdzobpyhh3JYHUPjzD9PXXX7emQnIjFLtBx3B6VrzfvHkzB7nBPIzt040bNyKajh1FCMxDyf6GsxNBzcAthFxBvhwG74bTQZx5aTypiA/B5Vy4cAEfccmw5bEzVBeE46lTp3pYA76HHfrXCum5wa8VHa66ckPPDlivB4m+WP3DdYqRA0vSRuzXvOWWW/BxeXnZkqShDzOgnIyZ8h3eorbphldz+OorkBP6pkPCALqZrV+/Hh7dAwcOTE1Nzc3NHTly5KZWTYi1Rw1eDCFfu379+s2uQ+8g9MjH1pnbb7/9T3/6k2XmBEpM/5MYOGcp9gVMs8AZA1kAxYB0qFAkrG+2bduGFAy8YASv47Bz506YZrg6NtNg4vlJWfMh4A1bXN+gWMwZwoQnJHfwCqO2kErwJLF1idNBDK1bt45vGoID8R5lzs/PY8YSl4AUFXyX6vU69CJXvsN8DV3tnL3Dbbfd9tlnn6XtnEM6afCJdCV01sTk7/Z0LjXiZyoD9oyyWrJm04oluB+C4QRWSIB9SI1GI1shBYUHV8T1YeGVNhjTOV39TTl//dd/nfMG3zcKKZFOgl87nGjpZIcwJIZ98267d+/GUBjyvKEr8vgCo7mriU+RNTsnjD/4R57Qyls8qAFwEm6c7tKlS/iIxw9RBuh7eKxvZWUFo4LoDrEbSkblER2+vLzMcRMoGaEN6JBGR0dRAo9PIpAB/TT6P2udb8ujMW4fZAySJG7vhL4Ytcig22jv3kj7s0uj58iU4D0/Cx12A2x7+dORPVE9uKJ4nns81hcflQ2fmkvLPuRmjY7min4K+xZCCDHADHWfnOH9TlPBgSHDY2iI7W40GtAfiH6GhoAFB02AbxuNBgcL+OxULpPlDkcZsAHoyx1xTVhwjIyMQMFwTDlfI3ZeXl5G4ZBN2Bkqja93YWEBu8F1jLgvaCOPX0d9MGcWazvhGn2dJ+yM7ERffPGFVz7RluzK7l5NQIRIJOO56Ipg5Dx+9TdxIw8OjIfp8Phgt2vXrvFAQhzE1JvoaUtXMQuJrbGTUbiMaKw4cj1xh5wjhSSEECIX9EGfuXo6cfolxj5k2CM80g1bzFOawi+C6GcEEWALfEueQRV+GsgLzC3l1EF+dl9nyCvAI+au2LCQK9cH5Zw5cwZvcHYcCyHFK01s3LiRBQr8SagqdJLHL/BVsCcJJW/evBm+K642+6IQdnH//ffz2rWYWNOV8yD41doibdQta6spvaHyXFfGm18cueMzHMysWCzyCsWsk7zloAViKIKfBX4tlUrx081tr60vs1s6DD3I1jEZB/ZcZq6QQhJCCJEL+qbnXA3ZQZaJY8ppQZy+ncemfRFYM7ty5QqbeIhJ88yk1nQ1FYtFBODBZwMJhUhxFy4cdIQyYe4FS+TxVxyA5zNzEfYNycIJUVgSeXLVeNSeT7dr1644VSt2cF8X7FO+5zgdG7O+HiCqCrXEKw0GJWS49EzcJNqmZQo0VhxWGvgOsRFKCA2Do0+9scVJsAJVzRNp2ZPk7lsz27VrFz8mOAQtEIMNTtsQ7bRMRdmFdM5QuT+lkIQQQuSCoVBIaSZGW1s71kZBlB3ew8AH9XodggC5RJFBB54S7IY4t4WFBczLgfskjo4rFos8Ms6mIm9x+YJptgBeH9iDi4uL7t3xo2AJ4nS+IAVLFnwFWxWXDyVnreYhBwdCM83Pz6MczMOFV4w1Fm7OF198gaPgeEP8IWoVGNQ9GIZDZVH2wM2KuYoVEiv1QqGA1oImwW7OIAaVhU4w+RRbcCxeQTwwsH37dk7HhY0YqEDkpwv9tOaUqJPSPNCuI1cTpphYjQFDCkkIIUQuGAqFFNM2IoXnbAe2CVs6UDlsi1lr9h3OoIpXzOApl8vQEJivA9sNWUw8Mi1e/SGegeHKBsXCuEPJ2D4yMoISOEcR3qPa8G/59A4YhvFdgq3qW1BVqCK8YuTdBSKnVeXsD7CCq9UqVCOrPdQq0EnxPc+IVorNz0Q3xsCbmdl0m0Yo7XYl3vw0JVEoFDg6jp8sNI8dO3ZAxOMjVFH8i7uTJm4GvnPaUezm2bRpUxythyeCBwaC2L8Og9xYsXVLH80ZuhFIIQkhhMgFQ9cPx3ZcogCKvwrmk8O8uuuuu3xnnwMEOwvOIcTUsTZyHQBp4jrGWrMk+MJ6vEOGYsMr56ADtVoNRyHWjo3QxBgn9vfEqb18IydxCGKcPFudNT1bnMoB6xPu3r2bg6CQ2RrvOUQquNi2yibNTehm/pBro95Iizj1LWn3lsX9yMhIbPXjSYEHcWJign1F3PgTk8tBl3OtgjSmgEvg9deXl5cbzcwO1myiaHs8ec5zRaIEDAx4nkkvylpbb6KgV9vrhKHrkECa3rfWdpMR9YCeBo8QvKw+2IU98dePf2d+otD0S6USP6486IdySqUSd0XsBPYwBOyPY9HfcAW818Gjgsw92Jn7Ld/ZwxasNV8RwEzYTZs2caeIVzzS6HuCBWqxMyrszmS8IqERD9PxQKIPRabF4nc4XhTsrCG7zgluWs+DSIGhw4O0GLZFEqlyucz/5mwtBf0N4Czd/IOuW7eOnyx+9DiuZ9OmTXwKfsoCsy8YrPad/YbgMUF7ZuNsNR3SKgfuEqey9AUashNCCJEL+qz/7I0gbtu3Z1jcGb5Zlj6wlZBWB+8XFhagYyA4WCLApsPAnYeB+ioS1moJerg2f8XaKAh4hUCBVoNXFiwvL3MMN2RTYj4VrglfMod9e25/gJ1ZJo6OjmJAAxfLV4T3KOfTTz/F/riNMI2DZEKdBOtnb+x5t0HixtnIiT9QrGi9GfO8AjQDtGc0GJbjGZUvFotp4xYYxPMnIq4JnlA8CB7sA/yptGiFJJ4/y1N3XVTFAwm8xWeCt3UWxPSsk/pOGDlSSEIIIXJBv3akXZFtawczXrNnwBUKhfvvv9+a8dxxVtBLly4h9hoj4xxICjxmIY6YYNerO5kAW5esjVxIAfiHefarF4L4C1ZsOBFiK8bGxnhwPx4E9y1xNC22wD4tFAocEY5TwGUFIIb+/Oc/QwnhI58CAnR0dJS9Yvx78Y8yhIpnNazGM5F2qzPmhwazuVECngUOrgFXr15liZM2w9RamyjvEKx1yS4ibqhonx6twCdibYQaeugEB/vENbSmUxm74W8BbTsxJkjtNg0pJCGEELlgKBRSNkFkMBObadu2bWNDLJ6pevXqVbaJ8B4L2bEVVi6XOZKNV6zw2YLQLuxq4uSqHh3E6grSB+V4MiGe63r69Gkzu++++6xp04HFxUWWTXGUHV+OtVqpqLwvlcY+JD4112p5eRmON9ihHJh+55134j5D4fFvFJPxqwUHyjK1jrVRsFssTDPuMOB8u56wh6c/43cPFtZDm4mj4xKjzOP6BD8xx6/GM1VXVlY4si4InLPW3EV8lLWKMwfeUJ5djvnpuPxG0grlab7txN06pN/n1UohCSGEyAX92pF2TqIRHZNo6MWulPvuuw+WO1YojyN5vvvd7/7617+2pmXEoWvsLykUCgg0wimwM6Yy+JwkjF/DlkyLLCo0lwxnbcQzJ7744gt8hP5gc4znq7rLit1UjCd+ZVGIm8DX5fOZuPLwIXH6123btmH1DdQBX/FJx8bG+IrS7L5Ez99qwvAGjx5M5tVY67GLEZRKJTQM/NYcqOayI57E0/bU7L904cJ6iLVaMMMpngjI2Y1dIfFoAbuFvJ5oqJzCGG5m7OCDBNnz6tIuML7e7N36VxsBKSQhhBC5oL+7057JmJOUNlHJR5kxWAzBwUPPlUrFzC5cuMCpR2BYsTrx6UpIM8rZTnnxb7e8OOUPm37uOuI8CGzoYWckOfWv/vjHP3odgjlJKBB14FXScaCfBRsx34Jvo4cIxj4AnChYITCOIQQfffSRmd11112oOVINpXkUEn+vjKH5IfQkfQUmcxxZx00C7zdt2oQHB80+du00mutKpOG/mvucrKmN8B6vnqkhntvH6t+XiIwFme/AH9OcZ4XWzEm88gWibVHOxYsXefIfYD80e1vjSx4ehqJDShvVaSuZsQV/35DhIyMjvFgqN33Eelpzaie6EByLnT/44ANrdjClUglNEDvwLFE00Pn5eQQg8Axc/LOj2vh/t9YnIb4uf+YxbsYDgHy9o6OjPBzB1eaxx0ajgccMG7nn8xOhr+JUY9gZF+6dHPZHjnB2IPv0WAxgJobYOsEIRsbAZqIjWrQlI4gAuCnDMQLclnxONy90xD+rlxbbN7zFX+MfkTMrBvELcfARSMytF5/OZ4KzmZVoEvFuuBUISvIxZ0Q2caBQkK8osfCAgW/AGrITQgiRCwZfISUO2sRfBVM+ecwKu2EE6ezZs2wPsrr3sGnOGYzAVqgKTE2FUikWi1BFvNwqs337dg8TsCj81MiaQ2TBgw8+aK2SxWMZ4jTeiXnEUW1Ule1KHhK8cuUKwhA4zStwUQV/NeskyDIOT/cc5GwI87jc22+/zQb122+/bUkDd4nRtLyDvx/Cwbq1IsMPb/TgxMN0aBIYo5uYmOAonrj8YrGINxh5Y4JDeLAucfZFWobToJx4FWauTyDIOC4pODyIlfAd0HQxXjI6OooBD57MABDsHpw9mE5rHQznDAZSSEIIIXJBnykkjMPC5jKzSqVSLpc3bNiQfVS2dZxoevBQLxZ1xUTORx55BNY624AsNazVvcTDxJBEbgHB2oKeQOHwsrijKE6Rwq/YbWFhAWf/zW9+Y2ZPPPGEtVqpLGKCElgnnTt3joOzYdz5Qhh+Q5CD0lpdBZwAxq+FlzuKzdjgt8BuQTnf+MY3zOzNN9+0SO4YWb58RdlegX6ktwa/hgQaNJYdQUNl0x6SCDOdd+zYEcyTtdaYoEKhwFOqOU4h+HE9SZW/xrHdQYXTBFCwQ6IXOfY+Jo6y8P8GV9inVWAhaZ78zq+J2VcT6d/G3An91CEdP37czKrV6uTk5MGDBw8dOlQul8+fP//jH/947969N7t2QqwxavBi2OibDml2dtbMTpw4sbS09POf/3z9+vXlcvnEiRPVavX48eMZz2dGvnrGLUG2U3hZ2EceecTMnn/++XjFPMSD+aoTGCaGCwTpQ7AzDEYY/kF4NE++84mxiNLBRxQeZz5tNBoYpP7kk0/MbHp62swef/xxvlI2D9lKRQmQZS+99BISs8JVxhKKPWQLCwtB8Lq1Tie0pmyCPciL8vFSsPV6nYPaUQ7by6VSCVGL2AiRijUG0+LFE39cv0t9Z1f23OA7xG9I27jwtOhTj43kdsJOQahttCifwQ1YSXCrtkgPBbi3lZVNrEus9UHmuvmsjDQnU0YQXcZ03TQp5n8XvIwFduO8WV9++WUswoIn3SKP3eDRNz6kjz/+eMuWLTMzM++8887JkycrlQrWvS6Xyx5vLcTAoAYvhpC+UUhm9uqrrx49evTy5cuHDh0aHx/3gXWf+5mI2yxtPUmxEQRjCm6eIBaILUH8QcAC2rp1K8LeME8WYCkKBJ7Bk1Sv11kQsMkJarVakJLHWnUbVE6tVsPHxx57zMzef/99a/qTsMWvy5e98Mojec/rr79uZg899BAEGfuNANuYxWIRcYMoAX+RmGzrMgj7x7km2Sq0Vn0GtxOUqBuk+AqnwA+Bq4uNx4AMCVXoLI9UTuitwWeQeO3ZisRJ888Vi0X4frhw+BrR7H1adDxQESwUmSgI4ppwlF3sASoUCnGm1MTL4RNxxF0s9Sz93yMxyDBRdeH5wmCAh9rGp4jfB7ciYxgg8Ur7qMFbf3VI+/fvP3jwoJkdOHDAmv5eoxmpifTXf5C40fRRe+itwQvBtDU1ckXf1HXLli2XL1/G+6WlpYmJCUiQarU6Pj6efWy2zZU4mA4ja3Jy0prGPlLa+MrHPEAMLw6SDjQaDfiKcBQsJmgI944YZWoAsJUgrXziEaYEYSA+noXu+onNTOyMqv7qV78yswcffJCXmYCBhn80KCT4xjx1As+34Og4PwtqhfrgRNgCdu7cyV4xVlocROeztfjOB+PvUEWwtZFACO43VDv++RLp3wH31TT4NNqKj4yj4i3jgv/kAAAgAElEQVT4iSGPvFgMJ3CWIP9xWcfw5BtsuXbtWlyxxKryrwmdxGtLumRh1xRPwgukhge5Wati831i4RLUJ00VBYmROBQwzgG2adMmjBPEMjHt2rM39il90yEdPHjwwIEDU1NT58+fP3r06OTk5PT09NTU1Nzc3JEjR2527YRYY9TgxRDSNx2Smb3xxhuzs7OPPvooLK/gYwZpUXaBvRNH3aBkmOqYfuQxRQBjwRhLQa6Es2fPIusdSvB52tYMxsN7L4eHnqEAfMUKrJIOOPEdZ4r0pdABR+shJu3dd9/FVcCDheni2AHKKVhknSdzcMmun7ADLhwOHp60ND8/H5uKkIzYgnv1+eef8wJ9nP0hmIYFPQT3Em8JQrPSSIyY6hd6bvDZdDuGk/bIuG8VTatB897wCHCAaCN9oZAgPi1eTA94fB3K4eQm8CoFUiye6xNor9jXmBhlxyfNuDlxxCArfvc98ysnRN68eTOPNGRnkUir5ADQTx2SmQXRrqsMfs34UfFnjZ8f8dD+OHHbxZ8+8EQpnD6Sx6yC+AWeBshjDkH6UX7O44RGpVKJZxriFOjJ8C+/sLCAFV1RE/QHCEzAdYF6vc5LN3EO1iCPOIrlGFaMqnEXbq19J2cJwuvWrVu5G8OJeH2mRqPBXz300ENm9swzz1jrAGlAPNrTyVd5Zm0bfG/EFg8ngPe02Z4iyFrTTfn/Kf/dp9mI/pH/uIPEWnGeraB8fhZ4sC7ogfg5SkvVU2g38zqxwhyMnhGzEPTQHC7UNtFq/45FZ9M3Yd9CCCEGmz5TSD3QoZXh9guMenhlwSuvvOI7lEqlOH0kAjoRpnz33XdjKANGIq+PwmWurKywlQQ1AKvfzSvOJITTcSCAz1Tl3KmsKmCgXb16FafATFUfmrPWqFO3K3G6YOEAL82T9nPUA6dEGhkZYenD5h60oC/CxAYjdBt28Ehxdk1j5Q6eZgvZF0d1W9Iv3keRdWvIKsOrAtOeQ3IwBut5s2Da43fBs5AYCZ2tjdrWNki6ygmEgh3ihV/jACJvOfHU2oBYAAWVT5NQceova31CA7JXWkmkbexDf4XYmRSSEEKInNBn/WcPZLgBYz+kNQUBDD1ewhJD5NVqFQXCqw97B5HTiO2u1+vwqcCEhyJhf1JQHzhOcCIODPWh+XgRWLbFNm/e7L4iaw3U9gPZJIQKgRuGq1EoFKBdIMsQPs6pY10R8kqA+ArBGu4tixdDwxacOohfZy8d7phfCF8LjF/PU2lk5/LdiD0BfI3DRrfh3U5wCP9kaNtBm8QEWPyIrNRZJVuzkafVKvGkHlNupAZi5cE5WK31seVbAfzwWL6DNB9VYoWDYA12BgcOVNQKHl+EgQRT4/k+x0ENHQ4G8Lf9iBSSEEKIXNCvHWnnBMOp8bC1j3HDKkFeFmz8zne+Y2Y/+tGPrJlF37OCIr4OQgplQluMjo6ygYZQTpQGXeJOGuRdjVeoc28Nu3PiavuwOE9EZbPRLUEWOmz9cV6fzZs34w20EWctCkrjZdbY2vVFzXl9QoArglj0wXoUiNNBY+EGoswHH3zQ12i3pgOMY8QRSXjmzBke/eeo98QffVBjkxLp2UwOvCa4q6yNsAN+oM2bN8d6GgRNN9G356+BLomDv71W8XUF7iXsyZm6EpUWjuIHkJdCD5xG/DGQaNyogmRIXBmMqWC2A9+fRMUWJC6KGdQ2LIUkhBAiFwy+QrIUSz/+FlYJZiAh7gtmPue8ccMfpuK3vvUta2bQ2b9/v5mtrKzA3udTYJCdI9M8BRE+YvEIzOGA6hobG+PRZKgHdg75+hexpwdXhGpcu3aNv4rH333KCC9+we4ujlkqFAqQKayuOOVJohcHG9mF4KeDlQ3j0ReJR8kcH4XTIZPTCy+8YE3BOjc3x1YkVzVxkuOg2pVrSyARfHF6a03sC/Xva5Nza4nv87Vr12J7H4okcCjGE3oSFxmJJTjq4+I7OHvw3muIndnXCwJ5l+gG5h3iuUfBsEe8VgvLoIwAv8Dr1kMb7qOBASkkIYQQuWAoFBLDVgknPigWi/DxYCNynzz99NPWtMVgqvtqNJzLh9dNGBkZgR6C6ImHnj0ajTMOYBSesyQUi0UOWOLUCbGDx1pnpEMbubkH64+FFIc/ud7y5TOsdd4GcGnFEylgk7LTqFwuxxOM2HgEly5dwv44HZ8C1vf8/Dz0InbgZEiY74IK3HHHHUGuVUuSv9ZqZvaFqbh6Vhll56/QMeyk5Ido27Zt/IDwHUboGg5nL461PhGBQgJpGUsDOHOEq/Y4yo5viJ+LZzL5qn3Wqo0ChRSHwAXFJu6AWsUZVYIJfDyhMDHXg5HLKvvqOv8qb0ghCSGEyAWDr5CCWfoZLgf4ijAbBl/B6MN2DJRXq1UsE8c5SWFquQDCKhK8eB37b4BbOtBJcKJgC6LvduzY4UkN/Cu8h2hw+5eNKbyiAnh/7do1z+lgrcYmwpB8Tgm+grOK89rxdArPEQk7DvcEE/VRq2KxyDYgbiBuAofAjYyMsE8C+g+hdC7RUAde3RzKj8u58847oZDSBNAw+5BWPwMJ8HQZbo0+8SjDfWhRCFxiDJtRmGtaMF4Q1RaPOvgEuLgm7DfygLo4VXEcKWdJjqVgtYv4urgceFgvXbrkLi4vjVcCdPA0JTq5jSbexQyG+h+KDilW4nHrKZVKvMgpRuHQhtBEMP5QKpUQdIDWjN4LW9CmR0ZG4JxHvwLfe2IQJ3ch/ND6GAK3e07SikPwh764uBjHLHD+U8+ailPwkAJ2w4TfUqnEvmIML+A93zEfGOSNPIbmY328/7333mv0RBmlIOKRHNQTN/PcuXM81ZfTGqFbwrioteaGSbM8hjPsuyvifsL/8blD4jkGHu0d+975zxoPTmJkdjB+zuXEj0wQPs7tDTt7CFIcEMETD3z6OTeDeGe/LXEnHViWfAPZkuPSarUadzzxlQa3nROt8uUE9Dwqm1s0ZCeEECIXDE7XmoanMQWxjexbYDrBMOeMojCsYLN46iBWEhAZrtNhQGFxIx764+VhSqUSHP6QCFwf38KWWjz71cvEcB8vU5QYLYqJjRhAwLVgiAxh7i6k4iE7ZmFhgUfM+NWD49nFzStcxAkurVXtBUNtceAJX5EHlUCPnjp1yjpgCLVRbERnzDBlIx0/9MaNG+M1fnj5CQ+/5rAabvbeniGSAqlkmTECGYOB8eRT1ygYAeb0RfzM+mBDnP+Xm1naskz+VVA3DkmIdRKnJgouwR/keG572p9V4l0aDKSQhBBC5ILBV0iWZE0kjrrybmxY4RUhAIuLi1A2cJzAo47BdE9SwstMxC5TVwOQLF6stS5xtrKyAs3hC256mYGTE7tBG6E+bMb6R/b0wBLEqT2snP00OAWLRS8T9YE1B7uP7cpr166xhznwGxkN0LNHlydIeg3hK8LdYHMVi/XBRVev12EO852Jf03BxAEF1ur7gRuPo06s2U7QZnDPg0juDrOCsk4KWnWshBJ/xPjhDaQGZ7fitodHAFvm5+ezT+cXnhaA7ncsvsxYXRWa84v5KUA5+D/xeRdxqrAgAH2wkUISQgiRC4ZCIcUk6iQO6OKYGU6QevbsWUgE7AAHDFIHYed9+/bF0Ts8vO7nYkGAY2F+grGxMV5ygkvjupVKJbiCcCw8QIhq82yn7Bzi5f6wBTHixWKRF3eIHWku7HjpAU6ICYF46623xslj+Ma6y4FvvvvVjNaYgNpDuCNfOOoGx16tVoOsDCKDTXRJHEvtHhfOecNt0t17GZHcaaeLf6kgXC3YyOVzADrn8vGdIdzjlSCCoQWOXmO9nqgdQbAUS9plxjczUJOs3vzU8T0J1hJMq1Xb29hHSCEJIYTIBX3Zi3ZFIX356gzLHbIjMcyMlQ1C6bAUAlIKrayssCpibcSl1et1dg7hKxj7Hj/GFlY8IO67YWOlUgm+wvY777wTb3iiLrKUvvnmm77z1atXWfrwnEHONpQ46wLayEfS4WDA/eQVJYIoLM6aCpMWegvlfP3rX8fYOk+MxekQKOi3hSO+sDPXMPjph2oeUvZlJk4Zxitm17mJHUegcdsOfD8sDli+eDpdwNray4klQqwDGo1GWgSag6bFs/FYpgROWVSeZxwmaov4r4AfT2sVYbwlCO0D8QopV69e5WN57lcQnsp/BfHlB5VPu5m5RQpJCCFELuiDPnNtCcaR/T1i3hwE5CCUC1/BhgoSQcZqaWxsLM2CYxvqwoULnq3Vv4K2CA5ku4bjcGBDzc/P41hUFfFRwFfTWKHlAd977z3fGcGB0IK1Wg2xSSxc4rnlIyMj7EWAv4cP8cJxXbgiDkPyyKLgKN8B4YKlUglpZ9lq5gA/qNI333yTs9ymaYK+sA3XnLZXHZvzIJ7i4294tg0IpE9cWuD7YTUcOJ/a2vtxfRIdRaiPR68l1sp9Nuw94kUogtlacaydb0lb0Qbb8Vjdeeed7tO1Vq3v2yHgeE4e69HANRtH/aXdqP5i6B7UOJzBJ82hhUHp428REh5dEUafCoUCT/bkSG5PHRT3efwkoGFt27aNd+A1N4PsO7G7lS/Bx77Qa/KJ0Bn4yrO4orvvvtuazwC2IxqiUCjEk2F5iw/lxYOQSI+EwcAvvviCqx2P83j3j/uMnOg8aoodxsbGcJ9x21Ft7sP8nxGn4+Vl46e00bqs58DTQwecFtbcaK4TBniSOKjX62xq8B9u0C3FgSf4Cw7SJGaPMgVh1jzl3Fe25aM4WT4aGDoqryo6yHidMK9M2zaTGAvur6jbjh07cOvwuPFz4QPsqAk/PnyXeBzeOuiQ+tQI05CdEEKIXNCXvWhXtNWz7Hu0pp2F9ZB++9vfWqvpt7Kygj1hc8WBAD58wcufBNmILUWQQbK4koiTCvM6rS47ELfNRiIOxCX4RxSOV8gO1LDYXLIWx3LqWL5wXO/i4iLXh3UbRN7o6GiBEs+wjYyTuprEECKv+wkTEnr09OnTnG2FbwXXx+cbcjxIhok9VDqpLWmTxLmxraysoJ3g58B95oZqSWmLE2V9YlB1Wq3i+iT+dqgVhwB4+WhOTDAsGY9McqhFMLk1lkFBOdzs+THkulnrHHlecHlxcZE38kkzrqItaT9xbpFCEkIIkQv6pucElUqlXC5v2LAh8WM2cQhpMDMutpjYkevTA9n2R4AsD6n7hFa2krAbe1/m5+eRPQUDxIicxg48HdVaTUIeoPcx5dOnT1tTWEB2QBvBtePXBR5++GG/Fizs5CEbMIRRAocPBOuB8p1BTeA98lSbbAhzxAHbhta0stlvxNpr69ateAP9hxJwr3DfPH0nLhyCKW312IzQ/5yzmgbfljS/Efv2V1ZW3MtorULcC4kDIuJAcGttyfgpMyz3tN8rcKLEJ71y5UrsWGIhHgQjoP3weivc2ALPVjxcYa1O1rhu7Ouy1ruNEvDv4TvEJwqUJe8wePSNQlpaWjpw4MDMzAxezezQoUMzMzOHDx+enZ292bUTYo1RgxdDSN8opHfeeWfPnj3Hjh1bWlp68skn169fXy6XT5w4Ua1Wjx8/vnfv3oxjYxswNs1GR0dhdGACKUdkwpZ304zT2vMWGFZXrlzhMLzYteNWj69nYU2ZEoy8x/YRm0Wo28WLF9n6gyPHtZGRpcm66vHHH7emTsLrpk2bIHQQvB6Ht3ohnOaHbwJUYLlc5ovlWDsujZMkeSVxE3wsHmHf7J8IUisZxfjiJmSkWuk7VtPgeyZuZj7RG78Oy1zfOS2uMnG9Cfhp4kX/EgOamcSJELx6iy9PDNEDEHWNRs67+fPFzTWe5xBUI566G3yMQ8n9NV5dAsQzahMvPwiO79At1EfeI9A31R0fH5+ZmZmdnT1//vz4+HilUkHPUS6XkSJBiEFCDV4MIX3TIZXL5fXr1586der8+fMHDx68fPkyAn7MDEu0ZcCjrrE544YS7BpIBISE8Vg5zHOf0fL6669bcwE6no20fft2mPb79u2z1mUsmFqthqPg/MDpeIUIrye7lNgUha36ySef8Pp4OGngs4lnleJWIIHQyy+/zLvBlvTV9oLK12o1fMUrsAHXNzAJcRN4GXK+5z5/iO8wLx546dIlDpmLZeKJEyeMXE0x/a6WVt/guzKi+TVwYLC/kM38tkk/A9hvlCaDMirvW/irxNUueFIaz8DlzFK+kZc+iZd6WV5e5hJifRNIFn7WYj+0td46tHBXllwOK9HEoZ3OpU9/eU/7pkN68cUXJycnDx48aGYPPPDA/v37fYJbW4Ox73SrEGrwYk3or8bQT3VlJiYmkE60Wq2Oj49n75w48mutP9XIyAjsGvhCOPFBHBhjzUg2XwPCKAUIrybOrhS2fXzFYp5Vw9LHWqP7gklOZvbhhx8aTcR55JFHrOlDwoG+tHkcaATTDILsqaeeMrNnnnkG1cNa4NjImhKFLC4usujh+Sg4fH5+HlISFw4Fia+w3WeHuOL0+sTZFvxErAJ52lOtVkuM+7IkKdzXdNXgVwnfwJWVFf6t8VwECilNPTA+xSf+dRLD8GKdFDhR4gpn7MBpDvxBSHMq89XV63VOzBqXnOjZCpxVRiMHvDNq4kMX7EvjZp+hkGJx1u/0TYf06KOPHjp06PLly+fPn9+/f//k5OT09PTU1NTc3NyRI0dudu2EWGPU4MUQ0jcdUrlc/s///M/Z2dlHH30Ug+lvvPEGf0wj0f7iLdjh888/R75OfMTwSDx/28O6YDAirwG7dnx9CrZ3YCvx4ni1Wg1GPacf3blzpzX9Sf4mTqiKtb0RFGdmP/jBD7wcwDX0j5wojxUSKn/8+PGpqSk/ETKuIoNcsJAg+414xNzT+EOZwRXH6bziOepGkVF+Ok/5BamKynNoHzh8+LCZ/fM//3Oc7pMJfv0+siV7bvC9ESsSd8nAq4fWEk/oabTmCeTF8QJiIZXoSYon4gTaKy48KI0To6DCcSrIQNlwK+IdVprr+CFWMGP4Ky0s0B8ffoRRQx4XqdfrCA7k+mRUmxVbRkRi39E3HRIIol07DH5NHB8IvvW/SE4ynbjeD7cPjDthYApTSjdt2oTpmVyCD2f59tHRUQSkcjtDHTAY1Wiu9AO4W0JniTIfe+wxjkqPE6Ra65KvPPWPo9XN7OTJk0El44yTQcnot3hF3V27duHvkgM9eHawj15iB3RF9913n0VPF3omHrHkCF2cGv+SXiynVgroo66I6a3Bd0JiV80OfDAyMoJGxenwg2bAH/FHGQ9weZh1hzXhLYlh1mmlNRoNVJX/+jlAJrFDQqfLQdu+G++PvwJ+MD3qPYb7xXq9jmPxROAQPCner6PPw0/A7ZntXe/+4xuYPWTaL/TNxFghhBCDTZ8ppB4IhhTaDt3AYPn9739vZt///vetdWDBD+FBLU79WSqVdu/e7XtyMEKwZg8rBqgizoITnJdfMbSIJEAeQ8GVZ0Hm8eWQcRhP4yhzT5rCko5n9vFdSozFCGb2YSMfBfsURihO9/nnn2MjbiAbjIGE8uDj4LdDUMlLL72EW5fmKvdD+t1y7IoOw77TAr6BR9CwPOXbGATWx18xvtRp2nBFIFkyfq+058IbSXwVGamM+F5BlPNggEe9Z6i3+O8lfnx85DMeioSe82h1nk7L072DcAaebtxfcXTZSCEJIYTIBYPTtWaQEaxpNGmOpcaZM2esVb64TxIGC09iLdAcT/f9xBZcYEPB7Q94JQgf9Y49zDx47SEPeMMh0ZzExcxeeeUVay7kCkdXHErQaDQ4SpUjuXn7ysoKZ67kdc/8QnjOLHva8B4SrV6vI4IDqgh3GAdi50uXLvEqZ3zrYGm6CRmsapj4E4tEWBXFXv14LY/gkEB9xq9BKtJYISU6ddJI9CQlHpUovOLt8Vd4hWThRptYQmJNYm3kJbBLiYt1rxLcpXiceagjdqCmXfgAIIUkhBAiFwyFQkoblWajwxN3slkXW9+I83ZgTMG0QfDbLbfcgnAaTsyK0qBOfCJn7KGB4PCZp7HjhLd4ABsLMkSEcy79er2OOsCHxPHTcQyPtQaqceCsh+SxbIp128WLF9mjxosQ4tVjzTncCLB16cFdbKfz/GIPZcTHOOkL6NPgulWS5lFrewj/ah7AjY9x6qmMtEzx6Tw8NR4zCKz+tkF08Sn4EA/nC5yIaXWLv+L27POCM2qS7RsLfD/sIuLBBmt9EuEi4sVZMDCwceNGn1+RWKt+RwpJCCFELhgKhZQWWcfB/h6uFqsitlyKxSJPU4CdwsntfRkLfEQwHtxFkC/ufeFwtXgukX+MJQvLhVKpFA8385XWajWIEux21113Wau+4Vvkb+JJrMHYd2xLeigdBsGRbw0BgagAgg8DpcUXzu8bjcYnn3xirfFavKCZ50PClBeOUIp/+uGkq+ArzjEaTMTBDmnCxRVJvEJ5ojyNdRII4tnYC9X2GoMQwbjw2HFl6YKML8RnX/Epgp2zK5m4vB6H4LoPj2PqMGDD7lt29QXXMkiNXApJCCFELhh8hRREAcW4DmB7BIksEQMGawW2yYULFyB3EKLGNqC7YeLheJSMuDJQq9UQWsb5hAJFwsqMs4zwSRcXF9l7xMtG4NTz8/O8wjp2u/fee61VBfoQP0+ZikOqrNV0Zf8N+6WsqY2gIHEb2SPl4NggDSVe4YrjcXa2TN3Tho1YukL0TFo4qLUG4MXLIrjYjZ2LaWIoPkVck7TI2CAzQuKcpDiejReecEmUeLHWqpk8ODAtlVHiPKTYPZx4jUEUIh4TjKmgPceNv9BMpgzphucrdl/1r2Ya/A4pGGWKh+/8PWfmRlPgOZ7+WKJBYBiKp5Fiu6eYQ/fDg0goE6mJRkdHsd4aZ/gOhgLwkRdMivO2rayssOeTR/m82pwiCK+eLs+o7+EoA44vCLKY8E3jB8YP4fVqcRM4Sbn303jw0NnHA0F//OMf4+wp6FlRT4SP+79PPF4E+vfJvNF43HNa1DIoFov8/84DyL5bWjrBIHIhHtpKjDtINIPiHfjUwbguV5ifCL9wLrntDNM4/oJLKzSTW3K0QtxNej8ar6ccJE3n3VC3YHmqIOTEOuj/+ggN2QkhhMgFg6+QLEnGJoY5QCZjI3zy//Zv/2bNKaWwWXbs2BFHKnPeyZ07d/L4AA9xQJKz7vGvOLjZp8einDRvqlebx7u4NGi1kZERTlTMxmxgx7ENyOoEBPPy4o1ebR4JZO3I8/6+/PJLfIROwkgjtvjAJqLkEd7NQ5FQSOfOneNqswjraiBo8OgwliFws6ORsPXtuyV68o3UQGJ0TMaWtKryR25age4PRiz81dtALJd5wi83laCSfOHBV/HlB0N/QWBFsCVYMTYeIPVy0NTxELGzwCsf69pBatVSSEIIIXLB4CukQtKi92yhuP0FCx3ahU1y4KPnseXFHqBgJB3pGjk1ju8QDzSzpbmysoKa8ITEeOHLer2OY9nJ5CoE14uPSGuPEuCA8aRBFmV24XvFp/bd4imxPn0yXsoFAdwIUgC1Wo1XnuX1mTD12EuGPI1lomu4DLOXGSQrck0IvPqsRTjs+9Zbb413Cx4i1uVeeOJ7azXq08IK4lP49mKxGGsj7IAWNTo6Gi/CEqcM9onn8Umhk3gVtKB62AGB8n4b40xXQaONvVD8X+HKD+0fNzPW/VevXsV4DMcTDZLrVApJCCFELhh8hZRIHL3j8/LY2oK+AZ5yETtwYhsYMnARjYyMYAYo4EBt9icFI/gsevxAnpHKhhXnTi0UCh999JG1Wm2wvFD5q1evesC6NZUfdmDtNTIywlWNg4XYFPWzQ2n5ioJGwVR8XQhO5XtVLBYRZMjGHdukKysrL7/8sh/F0zYRke8OpzRjMBDE/Wsz3iASnR/xtM3NmzdjJUle9jeei+3E0c9+ltgDGrxyI+dHhnWbZzfmCqORwK24b98+fjY9E5W1tmcPe+NWl6YCE3dAm0y8Fr/DXEKaby9wL/GfQBxlt3HjRs5B3JXHri+QQhJCCJELBl8hucnQ1pqAbQLdgFgvtpKgMEZHRyEyWF4ExhEmn8JlwgO+wBWPZ1m1pvXHM0xLpVIsWXgxcpTzwQcfsKGHQEHkUfXt2BPzeVErXj09cMbw6HZiChnExbHhyVvGxsbiiUoYFvf0tSgZBiZsQJ456xOMcJ/ZYITSQmkQWMVikWcypYV79bvZuLYkmuqBy8Sas77GxsZYVYMg500sfThqNGP+UNDAoPI9KZQ1HbpBdldut5zXGM/d7t27UXOeUMhT8bieQa1ix5VPKseD4zLFSNnEvqjEO5w4p4pPF+8QryhfKBR4vCGt/ET64imQQhJCCJELBl8hOfEoOfAxZZ4Hc/r0aTP78MMPfWfYiRs3boQJjyylXBrUSbVaZSOR5xMgm46vwsAL63FNvMzYnROk5wE4L0f6wWDE+3Xr1qHCANoIpiK2w/TbvXu3r2cRnJQ9CouLi3yxMZ45Asfy2oacO7XQzLDg6/VZc013lP/888/jMn/605/6JWMH/jXXrVv3q1/9KvEGimwCL2Yc8+nPBSty1jRB/FjaKYJpTBlRdhAEHD7Ha4p7dCWeKQgmjkPzLFxoJ9gBTxm3Rn984kzEcSbfxLsU531IvPDE1phWgis/fvTivywPxuuKPnoipJCEEELkgqFQSGnTHYIJFvASwYziOUDQNPBzXLx4EW94UhFKwG44xJpODnzFxhrKX1xc5KA7HAvTz0fPY2MqTl0VTMSBy4rdVx68xBGDqA90m3u5UAe+IsDp7+r1Oo4Kpkz5IfV6HR4snnAeLESG3ZDTjxc85LU8fGf4FfDVvn37rLm6PG7LwsICD6bHlqA0UxruUInVcJDLjhU/K6TAd5I4ZyjYYkmPYeDO4aPQktmBuri4iBaCdghRBU8P5rHV63XEYfKxsfSxpCeLcUWI6uEPga/R14TMFkx+UXH8YXAr2FXGTiCxxD4AACAASURBVLhg5zhzXdvm3UdPwVB0SGkEPyqnHEX7QyCpp/kxs4sXL+JJRrA1R77idefOnfyfjvccOOtBBLyaKp4udHW+QmvsGo2n3U1MTPBDgn95PtHZs2c59ho1wbPKXZSXEI8Y8H/Q6OgoHhgO1sDD40kqUQceA0GZPOC2sLCAy+Tux8dkcCCqio2IyOBpyKj8m2++mf2Y9cVD+BUTxGGnRWMHzvY4jXcwthwUbtH/YHyiIDMQPzhsyqCRBKlIPera8S3cu3BsQuL/cvaIpROPngWpjOIw8bg7yfjK/xC444c1FhgBaTEUg9HUNWQnhBAiFwyFQoqNoNin6m9gjrH/E1lKEeywuLiIwSi29xHg4IXAUkMJfFI2bbwcDiflsFS3TOOpiGxejYyMcLJRHgHDzmNjYxwgiyvibLCuS9g+jacW+oJMfF0c7oHdLl68iPogKSpuGk9l9YhhTtHEJjCCdxcWFu6//35rmr0Y/5yenrbWQHBcdXAtiT+uCPCGilEvNvM5qrvQXICHJyFkDHbFiU0TfwLeGISPx5OyeSpoMB01HkW4evVq/LwDjoBIDEZPHJPgZzOYk841SQzT4KLS6hOMvceqiJ/35eXlnht2YjB63pBCEkIIkQty3WcuLS395S9/gWcFVCqVcrm8YcOGxI/ZxDZFYNRwAALMcHiVDh8+bE0T6ZVXXoFggqUWL6y3a9euOHcOYOvSl2mITaEgRQoTG0du9HlmI2sdMXcJxWEauC5sd3OYgw5YzQR2a7A8oLV6oXfs2IHCoY3eeecdM7vvvvvMDBmV2AD0YyEToY38ZqJ6vG4s/2oIbUgcmuf3/aWT1rbBp+GNDW4/KCReeSFoeLEHPmO9iURxYJnmeeDTiksL/JrxKXi1Og9D4JkSnNjXy2HpEz9rXggrG1bn8dRUa5VBQQqi+DYGmom/iqM5+BoHmFwrpNdee+3VV1/1j4cOHZqZmTl8+PDs7Gz8UYh+Rw1eDDn5VUgHDhw4f/789773PXw8depUuVw+ceJEtVo9fvz40tISf9y7d29aOYWk5Sf8K4usOZjecFpAVbzyyitm9oMf/MDMdu/e/e6771ozPQ/mh8K0h3OlWq3ybL54UB4axWOKOA4NuP8mXv6OA43c1wU3DGaYQrdxqlafqYpa4ZI/+OADI61mUTZ+HrsP/AGx8Ruk2sRNQwz9vffea60SCrFzCPmzplbjxSxwXWNjY9gTx2IBC9wQaME333wz/u1iH1JAngXTWjX4DgnWFI/FbuK9BcEdzjbbvczsm++lsdOUF59M3J8llLs8uT48wTZI9JW9riAf7jvw05oYPZgWshjEdvPfQqCT4qkdPKiwsrLScxvOZ7MPyK9CeuONN/zhNLNKpYJOolwuz83NBR9vWi2FWCPU4IXIr0KK8bH1PXv2xB/TyLALEm1AtrbiWaKlUglWErw1iGpjj8vIyAhkE89yhSBA+TjEh7DZ8mJbKcjqz1/hdAjzq9frd955pzX1mWcmNVq8iyuMrx555BEv2ROuBLMrgrsUmHj8nn1jbndDx/DOvEqFewJwr6CTkKQVW3bv3o2bhhlIUKXQXrhwz4eUEdcUb+kLI9FZZYPPjqoKxC7fQNbZjeY0o9jXwi7PRIIKxL8RtxxvqBzJxh7WRFcT3qPN+EzVWM3E02ODWsVqKTgd79Z2NciMuUd8ycEt5VPEfuJA2HXekvurzedXIcUgRZuZwUIMPgoxYPTW4AtNbnT1RF/QX+2hP2ppZhMTE5jOUq1Wx8fHg4+rL99dTTBYYK1/61vfsqbWef/9981s3759HN/laQWM8iNwtB4vsge1BOfK7t27YenwzCEWZwsLC7zYBODM/Hj/5Zdf4j4gb0rsdioWi6gJTsRGmaeEMMot64LJSM0YJWXABCOuMCeYGB0dZQ8WT1TinBTr1q3zJdi9NGyBtbt161ZoIGSkxdlhuSNyL8N7lPb7Zu+QN3po8F396bj0iY+C08IDvbg5JUZ48tlj3Z929uAQa9UfrGy4fC+Tj+XYUWttsVxJRBW6MzXRhWOtjd99NhmpmeMVBQGn/4mPstYBhsTLX427qF96oIC+qfTk5OT09PTU1NTc3NyRI0eCjze7dkKsMWrwYgjJdYd07Ngx/vjGG2/Mzs4++uijMMmDj12RYXrANoEOwAQaWO4Ipdu3bx/i606dOmVmv/zlL83sxRdf9B1GR0choTAZKJ7ljsAzz4zA5idv+fLLLyEpAOsJlAND789//jM7qOD39kSoOJZDlXiNCdTTM58iqg0bcUs5iR+o1+tseELEeIWNFh5k1QiwAwL8JiYmeAELFp2esA7H4odgVwFvyTA/gy19oY1uXINnghk28RrewE11XoyRdwtETLyyeKCZ4t8IBEF6cT6FWKkHJ4IER3PdsWMHO114BlJi2BuvlhLfJa9bhuDgy2RRyI9zo3U1mTjtxcrKCoZk+anPaNWxpIurlLFDbsl1hxQTRLv2HPza1vuNVsITSH2M7oknnrBmh/Tyyy9bs29AS7p48SJGmUC8xis6qlKphAJ5RA4EvlmeEMd/Cp55CL0CvsKkXQzfgWq1Go+BIOMOdvZFOVF/7OCrhVrr6EGpVMJXuGQMA2IQ0scbeR4fL02LQ9DH+0K3+Oqhhx6yZng9Ai4ee+wxhHpjSSq+ura0/XH7iNU3+IwgeHzFy5LGN8rXH+Jo7CCuOh7UAsGQXVpQg3/LO/C/fPCPH18RpzCemJjg0XI+NYbsfHYB+irOjByPOra1eIKRz7hz4vnv1jTLcCwq7LM40Lw5JxlfadtIigz6aLy6n4IahBBCDDB9ppB6wIdr2hoI8TAChpIQZYtv33zzTcyQhVH/9ttvm9nk5KQ1ox527NjBK89CZLCF4qHbLLzwysvLVioVxEnzIhE8vof3mzZtwrEwrxACHsSw8qAE7CwoP16QyZdx4lw++Aq7eRA5RgXZUkb5GIu77777WGnhdAjd5iGL7du3I2wE9+f111+35lxjlPbKK6/wDgBqKdHY7yMb8MaRKIbSxpoajQYb3fHMa7CwsOBhJtbaXF0HpAmL7GGluIY4L3RMLDUSwxk4Tt0DduJJFBy5jvZ87tw5NGmclNset+2VlZU4Lhz4iFwQvB68enaxOBgdWzBS9/nnn2PdliAUIviBgmoMXoOXQhJCCJELBl8htY0PTnS6ArhnoJBgENVqNYQdI+MqHOyYF/LUU0+Z2TvvvBNnYeETgdHRUVijrIrYNLt48SI0VjxAj/BfOGnWrVsH85AjXznxqx+FUyACAif1FTa5bjCE4ShilYND8K3Xx9MgeWmffvopzEzOUbtCKwS6K45NcogqiE7Izddffx3aCCXAV8eX4xWWNgroPOQ3LV0Q38wrV664U9M3BvIlHgAI6oM3sZMpcTf4tDjMh91Xlv4IexqqOA479ktt2rSJmzfP30iMfeATBVvikRV2ngV3CU8WJ7TFn8znn3/Oii041t/7t101+D56OqSQhBBC5ILBV0huv2QHRyZGBmML8ngi+qtUKmGoF94aGPWw3+EImZycRAAeiKeOgnq9jq84PSuP6fuE2TjgG2YjBEelUsFX0BOxc8jHuOFk4tBq1igjIyO88h6HBUL/ubHGoXfYDYPg2LlSqXDWJTZm2RQ9c+YMLhDxftCC2AI78dKlSzgRYu04JVLwM/WR9Xej6W0uZNpR+Inn5+fRqHiKN7cfj/yM/U8drpUQaBEWBHEoeWKFUQfEp129ehVpq+KlYbieviwLVzuOOM9wXCWSFrPutwIn5exfeO/Lf6T9HF6BwW7wUkhCCCFyweArpAyDum3kD9tuCPq66667YCpCNsG9hIlHCMnzCT0QPfCscGJTD2BDYA+ALIAmQBhbvV7nhcY5LIfD84rF4kcffWTN6UQwFeFk8pIh5ng2CUsWbK/VaigB/id8xYPaPsjOdeBJSz4PEVFSuBZ2NWFn6MjDhw9DR+IrCE3U84UXXsDleLidJVmmg20nfjWw8Z7oHcE++I0wMICAyWDWTlo5GdMz4wC8xAk9ceiatU44jbd4OWkJhLwcTANKu3A/PM67Gnukgo1cDifo8t3idcoLrbllM27jYCOFJIQQIhcMft+bOOqaMRYchwzBRIIrZevWrbyOHL76xS9+YWZ/8zd/Y2avv/461uqGZ4XDhHgiwsLCAkeyseyARimXyzzTG64m1Aclw76z1rwMOBYVwLwoT/aK0uCtQa04GO/zzz+HpIM+YzEHL46bb3A1BUF3fPdwXYDjBrGEBEzs119/HZYgYuqQnA2uOB/6f++994L7k/gzdRhlp2C8NNoqG7Q9yG40LV6WxXMLsQSPNUqjuW5CW2dM4twjfs/6Psj8aySGMIEvDngLTp0WXhiQ5mAutK7/Gcs7n0fIlYwTCCVWMvBjxadLq1v/MvgdUqIvNC0RiCVJb6ZSqSAZXaE5T9aa/8v333+/mX344YfouvDowmnPyce878HQFife5oG7QqGA9spdEf+/+3Z0jRhw4+cT9UEFfCNqghIwz9Q7FR4JROGIg8C/D7q9QqHAzm0ObcDlbNu2DWfEV5zsDkOUOHxxcREfH3/8cWv2efx8Xrp0CSH1bY2JDp/GQXpou6JtNDZnpYpnX66srKBZoplxND/KXFxc5EiZ+J/d+620ITvggdEZOcIT6+941Ez23NLE+Ol4TC8xqD3x1PEfC99MWHUbN25EzxQHenj58f1hQ9Z/x8E2rTRkJ4QQIhcMvkIqJKUOytDp2bP8Ll68CAXDo1Xww//2t781s7//+7+HRICBv2/fPmsKBQ5O9XWDMPIGmw5hERhVs9aI8DhTpBuSMFph4vEoHCRLrVZjsxcxBQiygFpym46NRBYrd911lzVHHi5evIiNHPvAKV9LpRKit1nGYZgOlwDf+MjIyN/93d9ZU9ghcoHLeffdd9sOJTGDbTauhoybFo81xYltrFU0Q/3jl8Ij0Ha+apyP1ZKeMrf9E+MdjJp9mtLycTBuSPEwXfDUx6XFo4WW1LSC2vLHWAbdeuuteB45kAH4IF72dSW27cFr8FJIQgghcsHgKyRLMsfa+pDi92BlZQXBAj5P1ppWDxKtHj9+/Lvf/a61pmdFXDhMIeiSYGIsr0zh8Hg94stjy6vRaHDoAa4FgQk+ps9mL+QdRA8H79ZqNVwXFBt2Qxw2gLIpFAoQQBzMCi3oKz9hxi4uEydCxDyEIw554oknHnvsMTP7l3/5Fy8NNxOhDY3W/DdtXbiDZyp+xbDwZR+MDzCgNaKRo3n4wlesgWKVkxgeHUeKBz6ktDGMDKUVeFni3F1Mh1Hm7tqJ/zECbRQUbq1/OH5Xea0yvCJGKTE8h3cLCh9UpJCEEELkgsFXSG0NisQA0Lb7w+0BMx8GI8Kan3rqqV//+tdm9qMf/ciaZhGMIIgGRCstLCzwymBwMnHCnkajgf25Vp7g1WhqKsamET8N59Dp06eNcqJgT8gmvOd1bCFfPOwN0XQcF843ZNu2bSiB5Q4H1J09exZiDuXgLnE8LhbvePzxx3/zm99Y04eEa8QM32AKYTadx/TLyZRG7FkJJAIrAPyI8GtCTPtEiLRoNPeatP2l0jRxIEECH2pQ7ZGREQ5yS6NQKPAIAZ8okOPswuFsSYEgY+nGXih8u7y8jAcEDw50EsYSgmtMu64M19cgIYUkhBAiFwy+QgpIsywSbe0MYLZgSB2Zb2AnPvPMMxABWI3i+eeft6YKgTbatWsXDoethPxA8cBxtVoNUvJYc3IiJ2kdGRlBCWxGYQaSu3Y44yqAbYvAuWAonDcWaB4urMLl5WU4k1DJYCEMM3vwwQddKlmrNmKv23vvvcc5WIFPmbKODcAO/Und/rhDBbtz8FOyAnAHBlvrUEhQtxMTE2zUs57g94nzbBKD+rjtxSl64/0dDFQ0Go34qEBCWUpQX3xDPPgtDntzARTvwG3bZ5TjwcGtw0gGO4fc1xukQYovf7BbshSSEEKIXDB0ConJmI3UIUHCUDPbunUrpiXh4y9/+Ut/hUEEHTA2NobQO56XgPfIOvqNb3yD52kjiA42KV5hzC4vL6dNnPJ11Xi+CE6BuVDQSXj16R1smvEYN7bUajWoIpQALxf7gSqVCrQRm8BIx/Dzn//cmsppenqaI7XgT+qKDJ9Q26n1IoabhweC8g58V2Hse/IO3o1bMrsYE1MMBFviCUDcJr2eiTFsRu2Ql/jjiMHEiT5ps38SFVt8SGJN+NZ5ahJ+cPhAv0YOtU3L6DrwSCEJIYTIBUOnkNJUUbduhtgsQl67hx56iNdNePDBB83sD3/4g5n90z/9kzUlwqVLl6BgIBGQy4DDcs6ePcuzNDxFo7UmiBsZGeEgHzir8B5hddYalYdqc2ZYVGBkZIRnk/AgOLuLfPUz6CrsgAvBsPinn36K+uB08BvhJmCeE7xumzdvxrEvvfRS4s1s+3NkBNFlHDhUsXZptn8i3NiCWLv4dmELfJMfffTR3XffbUl58IItaT6bxBPFisS1e1q0Hlqmz2dK0zR+dbG7i/2sXnIQc+ivnN01OB0/OLhLZ86cwWMb18d1IftT0xj4pjtEHVL2n1FvvzQfhf/ll19+GbHgaNwYzcP02N///vdm9pOf/MTMXnnlFTwDOAq9FOcn3rhxI0fTYrcLFy74SRGlHYSlYhwM+YrQAbizFCUE+Xu88rfeeis/bDxdF7XyIHXUByMP2I0jFwqFAkZyvv/973t9MAjJ2aDHxsbQf3NKpK5+iG49vUPVFYG2nVAM7g8v3xXklOLbiGiXubk5/OL4rXm34DWOww4CuOM0r/EwsrXabTwRwuc58IQK7jaCMPT4XgVjaHy6jDsW91jcFeFZu3jxIp836PJR23iWxRCiITshhBC5YIgU0ipNj8RJfImjGRzjAOMRQ2SYKotps2NjY8iRw8YU5wrywiEpeHUJ7ID40TNnziCOHBuhtKBIXACxyYnEP1BXwVKwLJjwnkf5YIR+8cUXrLTi+/P4448jhAGDh1jhggcAPT8Q7klw67i0robsuA4ZseCDrZPikc9ud2MDH03Fmmo4dtrjZz137hwGADDZAEfxIUF75lYXCAX+deIsWaCQlCs5Hm+0pMjpxPJjMZcRsxAHfwdfcX4gaCOM1HmqsLTlJzZt2oSvMMEjQyAONlJIQgghckGuFdLS0tJf/vIXXt2uUqls2LDBt1QqlXK5vGHDhoxCEp2lGfZyhvTJ+MiHw6KZnp42s4cfftiaWUqx3t13vvMdMzt58iTiAhBajbhn9iqtrKzwGhAQFnDbAthQe/bsQeEYgIbfiM2rYrGIcuDZgr0GCw6mq+cr8vhUrxUOhEfBh+k5nhungHqD0+iuu+6CPkMQhEdM+HtEMdRqtTTXdOLt7TDIu60VmXMzc00a/OphU/2WW27hDMKA73mtVnvnnXesGa+c5gEKwrXbnp2bGTcMz+XKpSWKmFiRcJCFh1knuru4nOxYjOAjniwMVGAMAI9P4uUHp0ubDJvzRru25LpDeu2112q12rFjx8xsaWnp8OHD4+Pj1Wp1fHz82LFjhw4dKpfL58+f//GPf7x3796bXVkhVosavBhy8tshHThw4Pz589/73vfw8bXXXtuzZw+e1QceeGBiYqJcLp84caJarR4/fjzj+cxQNt36ITohiPtCLBnyjULKQGo8+eSTMDzhavrhD39oZs8++6xRJDSPs/MC0oi1g323efNmnIIdTtjZLwSWFzbyAmuw3cDi4mJgP/ohHApVKBSgh37605/6RjgPEOb+wgsvcLwfx4WzCky0LjPueZqcCoKA+9pFtFYNvi0ZTia+dR4VyVM1Y//TysoKGiTUOVogz5kN2iHXIbEmQWIea9VbhWZSVK4wN9pG0roS8W6JJXAFXJ3EYyp8oMfgQRvBeca5U720+PJxIb7oDLIkZ/9BBT7sNPr0EbA8d0hvvPHG1NSUf9y/fz/eLC0tmVmlUsE8m3K5jD9ZIfoaNXgh8tshBWAYfXZ29rnnnjt69Ojly5d9YB36Y/XwDNOeD7ckYwr/IEgahHUiJiYmYIEi+yr0xJEjR8zs6aefNrP3338fXyGEyVM0BuVDHvlHzEBi221xcRFeIsxIxVcwzRCn51t8Goe1WpeImkO19+zZw3OPsOIGogrdXcS2LdxmuITVa9MOPUmDQc8NvkMj2pLsfd6Cn7hYLHK+3TShYE3HCVoR0g3v3LnTaB2/7CC6YIYpx3byidz3w86q4Epjx1Ka/8aJFVKihstwVmFWFrQRu10TvaGoPFzCniEMGY/4XyjDz5rtgrW+fSL6pkMys6mpqfn5+ZMnT5bL5ampKbQAa/7dCzFgdNvge5gJK0Su6JsWPDMzg4cTHycmJhDHBZfvmpxildZEhrOK5yhAT+zateu+++6z1vnkzzzzjJHWQWYHDNBDCcH+wv8RDvGUP7Cz8J8VpLWHC8eT81tzkhC24NtisQgDHOU88cQT1jQ/kQQWHoKXX34ZcgdwNkyUViwWkQcWM5DaziZJvHXBnpZk4badENa2tJzTQ4NvO8GoE3s/2O5r0kMn8boJ7uBh3wycKGgtiNND+ykWi5zyIB6NcAkV5xXl90EChYx0O7HPJmPaEx/Cjisvn1NqsSyr1+s8K4tzRgT3HEfh+eKIRJTMDt3ECw/INj761zTpm3pXKpVqtXro0CF8/N3vfjc9PT01NTU3N4eRrq6I22KhdV2THv6zMg6M/w7OnTuHngbDL3AP4BU9xJkzZzztmzVHzDAGgnEwTyyEETnOxIPmjnl5jUYDJSDKHL0gYsTvuusur9W2bdvwSGAcD+EM+CdC8nJQKpV4zMFHY8zsgw8+QLU7iZvv/PZmhDN0W1S3O99cemjwGX9DiYNs2U290Zwki93QIXGAQ/BXy3/Z6JD4f3n37t28IDLPcMBucR45SzdHHC7Hx/e4R+Ea+hXFl8nv+aT1ep2tLk7A7x0kHhzO6Bj3qd4fc1fE9fETcb/bR811rch1h4QQI3DixIng2zfeeGN2dvbRRx/leRtC9C9q8GLIyXWH1Jaeg18z7PfVB3x3uAM+YmgLo3AIlt29e7eZbd26lc0xzDaFDOK42Hq9jqhrZANCCcCtMMgddrfCBsTAHbacPn2a1zHiNJcQWB5Wjh1QJlQRRmmCpUXXyriLrfjBG53rnK4afMYIXsb2GB6RQ6uDNx4jcoXWXD58Umh0NONt27bFualYLY2MjPB4F6/6yvhEb05BhEAbT6eN0B4eEOMbkkjaH4KHfXO1OeBiZWUF2ghPBC8BFRAPjQbrxgbB62lVHWyUOkgIIUQu6G+FlGdiP2rbeGX4b/BaLBZh6GG9V85z6oPOZlYqlaBsgC854SUvLi5yNn4eXg+WPmJBhtNhC3trP/30UwSvt40yaHtzOrQBb5xgHXhWr43MrNFo4E7Cocg+eYCViy1JAfDyxBcuXPCp3Jbi4Mn+ybzRslbjGdxQ/NZ0uLIgiwMcLD0uPAj7hgDix4dfr1y5ggcW8UQZUTy4aZCVOBG2+EOdFhA/VEghCSGEyAVDrZBuqJuhw5i9tJmGKysrmHbKU2Jhn8JdBEuzVCpxPDeDMnfu3JkWRwTgdvLkj3A1sdHH8/sSPW1cWofBipIv/YL/oGgG+FkRpekqx51JlvS7Q7589NFHwTxQixw8LN/jZfFc8WA3LgeN1l2ePJDQYUtLWzyw0Zo6CGVCM0H5Xbp0CY8JNsbPhT8dyAzEl4zARb8Va+t87VOkkIQQQuSCoVZIN8UYSZQOiUE+8TxBWKmYmgpfTlBCvLBeo7mEeXwitsiCGLl4mD6x2hkhcKuc1BXQoUMurVbDTJp07u1wzuWDmW3Xrl3LNu2x/cqVKxzGGWdftRRPj7U2UZ8XFc9k8qmssSuIy/THoZMW4o8Pdobaw4VgFOHLL7+M1yaPT1dori4BvxHUEi/y4gppTejfxi+FJIQQIhcMtUJiOo/76rxAS7KYOj9L9jwS/5YL5Cwm3dY2rliHczjSCuy2Dhl0NQ+pf83DG0HPP0SGlOfZPyMjIzw7J21CT71eh6bnGUiQ8sGUo3hWU/DKM+HgyGFBH0TrxXFxgSDjmFWugO/sCs+aqggX4lOO4tOlTaKy1ii7GzfxaM3/zb4ypJCEEELkAimk/5/eFEBXc7+7Ki2YAM+7cfmJTp3ELTzynuj7SbuioMy0mL22XrHefD/ZR/WvJZh/Em8sm/butuT5NBk6CbFwyMQINbNr1y6jqDkOnOPfnU/nuUKQChKSBXOhvKFywmJ2LwXXFa9uDjxrKnbGbqg8kkFw3GmhUAh8Rdb6yHjls5Vf52QPugTPe989HeqQeiSjWXSSXbTtVxl/BwHZCSg7jzLIdkoHp+OvOmz0N2Kore+et1wRmy89/JS1Wg0L+SDDXtwkghaIsS+khMc/OwIcduzY4UtqWVJD9SEy/MVz34Ae0buEuKknPhe8Mc4ajrNcvXoVPRMvAguC03FQO2rF3aTPe+3qbyGR+LrWdnj85qIhOyGEELlgcLrWm0jb0aq1osPBsQxWP/aVJg3bpkfqobY97C9WQw+aqV6v+0oo1hQKHPsQlMP5hKAhfGwNceSY6B1MXTBKYcWpVDmcwYfdeAEwwOsY+ZWmtV7esrKywllTQbwIhTWjFRDPjYnDPEgeDNll0OFQ2yCpIkYKSQghRC4YzG62Q9oa9Xkzz3urTyfX0q3hFs83XE3hnddTrCEZcjntY6JziMMHOJEVLywbHMtJUd1Jw0sPIzMW/FIQUr7KAwd8g2DuKvw9WH8SqisOaggkC3t9ghVgWRvFa+gl3h++/EAhZTCooqcrpJCEEELkgqHuk7sKe+tf1uRavgI/2SDd874mMYyYPwY6iefJ8lfwtXj8dFoCU1cS2AHL1kEGQRVxlqAvvviC1/Hj7Fa+nArcOZyVGDqJnTrLy8t89kAV+en8KH4N3C5xJwAAB8NJREFUwr4RoYdIP2gspFu9ue7SPn2apJCEGF56+9vq+c+Owxw65+233+7tdFhhuVuwqFgP8MpkndPzzezTLieboVZIQohsGq1reCf+CfJKj4AnBgHPLhqX5pFy+AhpgtKglszs7NmzLoMS/Vi8BSWgX3nvvfesOdUJkXteDmsjTDByt9DFixcT8wzxq697iZg6VJUlXYBcpJ0ghSSEECIXSCEJMaRg8lmHs8TafsU6iZMX4NV9Pxl5TDhFUDCTKV6dISMmjfeE9IFrh6cHLS8v8+mCsUT/6CKPRQ+0ETxk1pyBxJ6trpxwbelNWvWpIPva9evXb3Ydbgi33377za6C6AM+++yzm12FtUENXnRCzhv8wHZIQggh+gv5kIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDsnMbGlpCdk+QKVSWVpaantUpVLhozo/cHZ2lnfr8Cgzq1arvR3YQyWr1WpvV9fb6eLdejuwt0oOG2rwiedSg7/p/N9//dd/vdl1uPm8/PLLn3zyyT333GNmhw4dqlar09PT5XIZyYZjlpaWvvOd7/zP//zPW2+9dfHixW4PvH79+i9+8Yu/+qu/2rFjRydH+bF/+7d/+4Mf/AAfOzxw9+7ds7OzXs8Ojzp+/Ph//dd/vfXWW//7v/87MTHR4elmZmaee+65t95666233vrFL36xd+/en/zkJ53fzOeee65YLHZ+uqWlpX/8x39cWlr693//9/Xr13d1M4ccNfgANfi8cH3oefjhh2+77bZnn332+vXrf/jDH372s59dv3798uXL//AP/5B2yIsvvoj9r1+//s1vfrOrA1988UXfrcOjwM9+9rNvfvObtVqt89MF33Z41Mcff4zdarXaD3/4w84PdM6dO/fss892eNSrr76Km1mr1bBbhwc+++yzr776Kg58+OGHu63k0KIGH6AGnx+UqcHeeOONqakpvK9UKjBYyuXy3Nxc2iH79+/HGxfLHR545MgR379cLnd4lJlNT09v2bKlWq1u2LCh89Nh/+PHj69fv/7o0aMdHvXxxx9v2bJlZmbGzE6ePNn56Zznnnvu5MmTL774YidHjY+Pz8zMzM7Onj9/fnx8vKvTrV+/3sw2bNhw/vz5bis5tKjBB6jB5wf5kEJc+e7Zsydjn3K5PDs7e/jw4aNHj3Z+IJiamnruuee8MbU9qlKpVCoVf7Y7P93S0tLmzZsnJydLpdKTTz7ZeSVfffVVM7t8+fKhQ4c6Px2YmZm555578D/S4c1cv379qVOnTp06tWXLls5Pd88990xPT586derJJ5/Eg9p5JYWjBm9q8LlBCinEHZvZRsfU1NT8/PzJkye9TXR4oJkdO3bs6NGj3/rWt7797W93ctT09PTmzZunpqaq1erx48fxj9DJgZOTk5OTk2a2d+/eU6dOdV7J/fv3Hzx40MwOHDjQ7dXNzMy89NJLnR/14osvTk5O4nQPPPAA3nRy4N69e59++ulTp04dPHhwdna2q0oKRw3e1OBzgxRSCxMTE5cvXzazarUKNZ3IzMxM8HB2eODx48fxnMCe6vCoI0eO3HPPPbDCJicn169f3+GB09PTaLhdVdKtNmsO0XR4oJlVKpXx8fGuri6mwwNPnTp1/vz5Y8eOlcvlvXv39ny6YUYN3tTg84Ry2ZmZYUj92LFjZnbgwIE9e/bMzc0dOXIE1lbM8ePHz58/D9VsZr/73e86PLBSqTz55JP79+8/f/78Pffcc+TIkU6Ocg4dOoRz9XA6GGUdng67+VEdns7MpqamtmzZgkM6PKparR46dAiVHB8f7/xXWFpaOnDgwP79++fm5p5++umJiYmubuYwowYfoAafF252VEUe+fjjjy9fvnyDDqzVasFuN/R08W69HdVzPW/o6XAzEYi1mkoOOWrwibupwX/1SCEJIYTIBfIhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQuUIckhBAiF6hDEkIIkQvUIQkhhMgF6pCEEELkAnVIQgghcoE6JCGEELlAHZIQQohcoA5JCCFELlCHJIQQIheoQxJCCJEL1CEJIYTIBeqQhBBC5AJ1SEIIIXKBOiQhhBC5QB2SEEKIXKAOSQghRC5QhySEECIXqEMSQgiRC9QhCSGEyAXqkIQQQuQCdUhCCCFygTokIYQQuUAdkhBCiFygDkkIIUQu+P8ASKsLpRJ9e4oAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%-Display\n", "A = spm_slice_vol(ref,spm_matrix([0 0 70]),ref.dim(1:2),1);\n", "F = spm_slice_vol(src,src.mat\\ref.mat*spm_matrix([0 0 70]),ref.dim(1:2),1);\n", "subplot(1,2,1), imagesc(rot90(A)), axis image\n", "subplot(1,2,2), imagesc(rot90(F)), axis image\n", "colormap gray" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "SPM12: spm_coreg (v7320) 10:19:01 - 15/03/2019\n", "========================================================================\n", "Completed : 10:19:19 - 15/03/2019\n" ] } ], "source": [ "%-Coregister\n", "x = spm_coreg(ref, src);\n", "x = spm_matrix(x);" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGkCAIAAACgjIjwAAAJMmlDQ1BkZWZhdWx0X3JnYi5pY2MAAEiJlZVnUJNZF8fv8zzphUASQodQQ5EqJYCUEFoo0quoQOidUEVsiLgCK4qINEWQRQEXXJUia0UUC4uCAhZ0gywCyrpxFVFBWXDfGZ33HT+8/5l7z2/+c+bec8/5cAEgiINlwct7YlK6wNvJjhkYFMwE3yiMn5bC8fR0A9/VuxEArcR7ut/P+a4IEZFp/OW4uLxy+SmCdACg7GXWzEpPWeGjy0wPj//CZ1dYsFzgMt9Y4eh/eexLzr8s+pLj681dfhUKABwp+hsO/4b/c++KVDiC9NioyGymT3JUelaYIJKZttIJHpfL9BQkR8UmRH5T8P+V/B2lR2anr0RucsomQWx0TDrzfw41MjA0BF9n8cbrS48hRv9/z2dFX73kegDYcwAg+7564ZUAdO4CQPrRV09tua+UfAA67vAzBJn/eqiVDQ0IgALoQAYoAlWgCXSBETADlsAWOAAX4AF8QRDYAPggBiQCAcgCuWAHKABFYB84CKpALWgATaAVnAad4Dy4Aq6D2+AuGAaPgRBMgpdABN6BBQiCsBAZokEykBKkDulARhAbsoYcIDfIGwqCQqFoKAnKgHKhnVARVApVQXVQE/QLdA66At2EBqGH0Dg0A/0NfYQRmATTYQVYA9aH2TAHdoV94fVwNJwK58D58F64Aq6HT8Id8BX4NjwMC+GX8BwCECLCQJQRXYSNcBEPJBiJQgTIVqQQKUfqkVakG+lD7iFCZBb5gMKgaCgmShdliXJG+aH4qFTUVlQxqgp1AtWB6kXdQ42jRKjPaDJaHq2DtkDz0IHoaHQWugBdjm5Et6OvoYfRk+h3GAyGgWFhzDDOmCBMHGYzphhzGNOGuYwZxExg5rBYrAxWB2uF9cCGYdOxBdhK7EnsJewQdhL7HkfEKeGMcI64YFwSLg9XjmvGXcQN4aZwC3hxvDreAu+Bj8BvwpfgG/Dd+Dv4SfwCQYLAIlgRfAlxhB2ECkIr4RphjPCGSCSqEM2JXsRY4nZiBfEU8QZxnPiBRCVpk7ikEFIGaS/pOOky6SHpDZlM1iDbkoPJ6eS95CbyVfJT8nsxmpieGE8sQmybWLVYh9iQ2CsKnqJO4VA2UHIo5ZQzlDuUWXG8uIY4VzxMfKt4tfg58VHxOQmahKGEh0SiRLFEs8RNiWkqlqpBdaBGUPOpx6hXqRM0hKZK49L4tJ20Bto12iQdQ2fRefQ4ehH9Z/oAXSRJlTSW9JfMlqyWvCApZCAMDQaPkcAoYZxmjDA+SilIcaQipfZItUoNSc1Ly0nbSkdKF0q3SQ9Lf5RhyjjIxMvsl+mUeSKLktWW9ZLNkj0ie012Vo4uZynHlyuUOy33SB6W15b3lt8sf0y+X35OQVHBSSFFoVLhqsKsIkPRVjFOsUzxouKMEk3JWilWqUzpktILpiSTw0xgVjB7mSJleWVn5QzlOuUB5QUVloqfSp5Km8oTVYIqWzVKtUy1R1WkpqTmrpar1qL2SB2vzlaPUT+k3qc+r8HSCNDYrdGpMc2SZvFYOawW1pgmWdNGM1WzXvO+FkaLrRWvdVjrrjasbaIdo12tfUcH1jHVidU5rDO4Cr3KfFXSqvpVo7okXY5upm6L7rgeQ89NL0+vU++Vvpp+sP5+/T79zwYmBgkGDQaPDamGLoZ5ht2GfxtpG/GNqo3uryavdly9bXXX6tfGOsaRxkeMH5jQTNxNdpv0mHwyNTMVmLaazpipmYWa1ZiNsulsT3Yx+4Y52tzOfJv5efMPFqYW6RanLf6y1LWMt2y2nF7DWhO5pmHNhJWKVZhVnZXQmmkdan3UWmijbBNmU2/zzFbVNsK20XaKo8WJ45zkvLIzsBPYtdvNcy24W7iX7RF7J/tC+wEHqoOfQ5XDU0cVx2jHFkeRk4nTZqfLzmhnV+f9zqM8BR6f18QTuZi5bHHpdSW5+rhWuT5z03YTuHW7w+4u7gfcx9aqr01a2+kBPHgeBzyeeLI8Uz1/9cJ4eXpVez33NvTO9e7zofls9Gn2eedr51vi+9hP0y/Dr8ef4h/i3+Q/H2AfUBogDNQP3BJ4O0g2KDaoKxgb7B/cGDy3zmHdwXWTISYhBSEj61nrs9ff3CC7IWHDhY2UjWEbz4SiQwNCm0MXwzzC6sPmwnnhNeEiPpd/iP8ywjaiLGIm0iqyNHIqyiqqNGo62ir6QPRMjE1MecxsLDe2KvZ1nHNcbdx8vEf88filhICEtkRcYmjiuSRqUnxSb7JicnbyYIpOSkGKMNUi9WCqSOAqaEyD0tandaXTlz/F/gzNjF0Z45nWmdWZ77P8s85kS2QnZfdv0t60Z9NUjmPOT5tRm/mbe3KVc3fkjm/hbKnbCm0N39qzTXVb/rbJ7U7bT+wg7Ijf8VueQV5p3tudATu78xXyt+dP7HLa1VIgViAoGN1tubv2B9QPsT8M7Fm9p3LP58KIwltFBkXlRYvF/OJbPxr+WPHj0t6ovQMlpiVH9mH2Je0b2W+z/0SpRGlO6cQB9wMdZcyywrK3BzcevFluXF57iHAo45Cwwq2iq1Ktcl/lYlVM1XC1XXVbjXzNnpr5wxGHh47YHmmtVagtqv14NPbogzqnuo56jfryY5hjmceeN/g39P3E/qmpUbaxqPHT8aTjwhPeJ3qbzJqamuWbS1rgloyWmZMhJ+/+bP9zV6tua10bo63oFDiVcerFL6G/jJx2Pd1zhn2m9az62Zp2WnthB9SxqUPUGdMp7ArqGjzncq6n27K7/Ve9X4+fVz5ffUHyQslFwsX8i0uXci7NXU65PHsl+spEz8aex1cDr97v9eoduOZ67cZ1x+tX+zh9l25Y3Th/0+LmuVvsW523TW939Jv0t/9m8lv7gOlAxx2zO113ze92D64ZvDhkM3Tlnv296/d5928Prx0eHPEbeTAaMip8EPFg+mHCw9ePMh8tPN4+hh4rfCL+pPyp/NP637V+bxOaCi+M24/3P/N59niCP/Hyj7Q/Fifzn5Ofl08pTTVNG02fn3Gcufti3YvJlykvF2YL/pT4s+aV5quzf9n+1S8KFE2+Frxe+rv4jcyb42+N3/bMec49fZf4bmG+8L3M+xMf2B/6PgZ8nFrIWsQuVnzS+tT92fXz2FLi0tI/QiyQvpTNDAsAAAAJcEhZcwAACxMAAAsTAQCanBgAAAAddEVYdFNvZnR3YXJlAEdQTCBHaG9zdHNjcmlwdCA5LjI2WJButwAAIABJREFUeJztvXGIXFWaxv3O91WgAm1Sgc7SQpckoI6hW1AwYKIBHUdI764LBpzogLMJTGQdmBFnJ7MwOxkWw4ysLSOZWUawYZXxD22FKAiTAWdRUDtCXCJYIc5E1pBqsNk0pHpSkF5SkO+Ph3q/p86599bt6o59q+r5/VFU3br33HNvnVvnfc77nvd849q1ayaEEEKsN//PeldACCGEMFOHJIQQoiCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFILSelfgevHNb35zvasg+oA///nP612FtUENXuSh4A1+YDskM/vLX/6y3lUQhebWW29d7yqsJWrwIpviN3gN2QkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIRiiDqlUKpVK1yXM/fqVvFKKUxMhhFgpQ9QhCSGEKDJDZE23Wq3ilFwqlda2PhBG1+8ahRDieiOFJIQQohAMkUIqFGsuZboWKAklhCg4UkhCCCEKgRRSP7EalSNtJIQoOFJIQgghCoEUUj85V+QoEkIMMOqQBurve5CuRQgxbGjITgghRCGQQvr/4aQ7khpCCPE1I4UkhBCiEAy+QsqfpKdnVbT6UAIFIwghhBSSEEKIQjD4CmmthEuGiFm9srmu2kjySwjRF0ghCSGEKASDr5DWitUrrdWX0/UUXQXcitSSpJUQ4utECkkIIUQhGAqFlGbpd1UAiV/FRwVb4h2+hgC86x1JKIQQ1xspJCGEEIVgKBTSGsqCxFlNwZY8pwvKSZRB6+7CkZwSQnydDEWHBNZkJC3jEC8/T0fSQx+Wn3XvyYQQogc0ZCeEEKIQ9FmHVKvV6vU6f1xaWso+xEOuW60Wxx04a1U3Lj+7Pjgpnz2xJmkFdq02Dlzbq1s91+OeDzw9NHgh+pe++WtYWlo6ePDgxMREvV6fmJg4fPjwgQMHqtXqmTNnfvKTn+zevXu9KyjEWqIGL4aQvumQ3njjjV27dh0+fNjMHnjggcnJyWq1evTo0Xq9fuTIkYzns6tTZ809LtnRCv6+5zOur3Oo59sln9aK6LnBC9G/9E2HtHfvXrzBkEWtVpucnDSzarV68uTJ9ayZENcBNXgxhPSND6larVar1bm5uYMHDz7xxBPYgq927dqVfWyaM6ZF5KxGttcn2BOx3XncOXEdVu9oWdF18XlXWlVxPVhNgxeiT+mbDsnMpqenZ2dnjx07tn//fjNzZ68MRjGQqMGLYaNvhuxmZ2fn5+ePHTuGj5OTk7Vazczg8s1TQnaGU5+pmu0jSfQAlctlM1teXk7c09/n9L58bbOI4sU1bIVVFdeP1Td4IfqOvumQEP964MABfHzllVdmZmamp6dPnjx56NChda2aEGuPGrwYQr5x7dq19a5D78zNzWGoPf7qm9/85l/+8hfeknO1vR4WaAA5pU/Xna+TOknUQ/GWoVJFt95665///Of1rsUKWFGDFyKg+A2+bxRSIgp+FUOFGrwYbPq7Q1oR2Vns8mQ7TSvTMy+kHbIi0XOdNEr25Q+VMBJCFJMh6pCySYxWWKu1iPR33zMKrxBieOinsG8hhBADzFAopK5rvOY5xLdnrzybuIPoGS3MIcTwIIUkhBCiEAy+QgpmvIKeo67lLloRRVMkxamJECJGCkkIIUQhGHyFlLH8ROL2FRnRcR7Scrkc5BC6ToyPj5vZ9u3b/bxmNjIyEtSq0WjgK8BfNZtNM6tUKmY2Pz+PC0fCtIWFhTWppBSJECI/UkhCCCEKweArpPwzXtMS6nSd67qiBEKr4cEHH/TTQfe0Wi1IHJZl2KHRaJhZpVKB3IGiwvtYSKEQM7v55pu9cOw8Pz9vZp999tn1uy4hhDApJCGEEAVh8BVSoosow3uUlgfIt8R+o+ysPD1z2223mdmOHTu8WIgeiJhbbrnFzM6ePQttFDuutm3bhu0PPfSQmZ07d87M8P6rr74ys/Pnz/tZFhYWXFFZ272EMiGt7rjjDjObn59/77331uoCxXVllWs85pyrl7EPV8AXfmRnJw5MbMD4Kn4kA7ADyvRVYOLH9sYbb7R2S8aDsLi4mOdRDa6xZ/dwUPl4ZMVyD7Ss6B+m75y4g98hJdLD7x3Me40fkrX67dGR7NixAwUuLi6a2X333ec7YAwNHcPZs2c5AIH/RPz5jP9Z4uehUqngPwL9H4699957fedPPvkE2x955BEz++KLL8zs9OnTa3LJMZpfvO6sVcas4G8X/+lxR1Iul7mhwiTibik4KRfLZSYueXzx4kUzu3z5spfsrGh9ZB8nTyshrnBcYOL7xNkpieSZSpFoN/cFGrITQghRCPqyF109iVbGilZ0xWscSrAaux4xC7DCms0mhuYgXABkygcffMBHsdWGVwy7+XsO72bzkA9sNBqQYhi4w0n5PUb2SqXSxx9/bO0RvCtXrlh76C9xNKPne5Iz6kR0JafdvSJ6+1kzxIGRssErN/sVDYkH5fBzgSbqzT4t533X0/WQeCz/vUqrVW+/YNFmpndFCkkIIUQhGCKFtCIBlBYpHqxS0YPZErhJb7/9dmvrDygSKJXJyUmoGby+8847/p5POj4+jnLuvPNOM3v33XetU9ksLi6OjY2Z2WOPPWZtNYMt2BnBEc1m0529fiK8x86wLm+++WafSOsn3bBhg5l9+eWXFg3Qr61p1keG3lDBDiHekqFo+Zlybc3xDjfccENQQs5p7O6LSqsnTnflyhU8JjlVEX+VEQyVdo2JdP336Oo8SztR/y4ALYUkhBCiEAyRQkpLHZT4bayKupoz7vvJc4pyuYwg1MnJSWsLF8TXQZEsLCxg49mzZ/1Y7AB1gt1+9rOfvfXWW14sR9PCACyVSgjDQ1wcKvnhhx9ae94rXm+//XZIJQipV1991a8FO6C0ZrOJEuDf+vTTT/0ScOrTp0+nDYJfJ/pulPxrZq28R2mGubMi50fG7xWEzAWvK40fu3r1alCaR5nHcmet2lLarPm0jWknTfz3+Jqfr68ZKSQhhBCFYIgUUkzGPIAVmSE4xLVRtrqCwhgdHYUigfK4++67rS19sMVD6UZHR60dgPfwww+b2Xe+8x0z27Vrl5m9+OKLSIeKE2Fn1ATvW60W3iD3D7uXAAzGU6dOoWLwJHnOVms7h+BPOnXqFJxe0FvQSazYtm3bhq+yrbaVTjPKMF0HyTzsU3oQFrGGCNweaaX5bl1VSPw+2Cd+QoOkXNkldGVFtyXjiRiq4FIpJCGEEIVg8BVSoumRGJES21wZbqfYNEtMBBJP84Y62bVrF8QK4us4tg3aaHl5GfFs999/vx/72muvWTswD3F3Y2NjUDMserAzJJEvh4EdoGM4HRFOfenSpY0bN1p7NvupU6esLcLwitIWFxchofhaOFXr7bffDp8TT0vqGsjU1ZYceMOwyHTNZZDzkDgMj7d0LTY4JO3ZDKRPWskZantFwW+Jx64msm5NMjX0L1JIQgghCsHgK6SuxribNnlsnGA2wIrmW0BJTE1N4StoIzhgsMPJkyf9kEqlAmMQG6FjfGqRme3cudPMRkdHIaR4B6gZyCD/KtZqwKUVQvhwdsgyxNohGhCFPP300y+88IK1BdCJEyfMbM+ePb7D6OgoJidBY/HEqa5D4UM1Vj6opI0ZJKaYC0jT0Pw+IyFIYgqGxLGQ7FPzzsHznlOdxH6pnIdkfLV6n1ZfMPgdUk4yfu/EppnfV+mvCElA/9FoNBDVjfyk6JbwFWIQEBpg7X6CR+QQSoDdtm3bhiEynALdAMIlUJp3tIh9wHxDPNXoEWu1mpk1m00e2cN75FFFj4ixvunp6TinC3ZAR7u8vOxjktaee8vhtvmDgLPvcxAEPNhP6TriHUnXOxx3ReghMqIVgjKzrcDEx5B34HjxxDoEs3HjmnSNhug6UJ/TkM0g57DhoKIhOyGEEIWg//rher2+adOmzZs3m1mtVqtWq3ifRmJQQ0b8QpptErhMY5srw4TklV6x2/bt25977jlrm3UYXrt06ZK1h7luv/12GHfY6GtJWHtxF0iWUqnkU1atrY0gobDbZ599xtKEU6xChEEntVotn0hrncOJKBPB33v27MEO2Jkn7SLI4uGHH+bIcgz9Qbd1DfoItmSb5CsNi+hfVtrg15yuEfzxx8RB2rTFhHJqL39NG9plSeQf0Q4RsMNT1y9fvhzMYbfOB9lLi08XHBiLMDwCiYot4xpX1HoHdXy7zxTS0tLSvn378EAeOHBgdnb24MGDc3Nz610vIa4LavBiqOgzhfT8889v2rRpaWlpbm6uWq0ePXq0Xq8fOXJk9+7daYfkHKfOINEvGqsi/4rVDLSRqxnf+T//8z9hZ2FnLOFarVa9zG3btuEoaBosPAGwHa/z8/PseYI2gj0IXWJtgw4pUH22rNcQXiW3BDnegRdJ86my7ENCOANO5DoJQQ0syOBPCpaazY7Fz8+KnHl9Z0v20ODXikTpE5MoVeMIb5dHaQ6YxHDteL7qyMgIuyTj9V+2bNnCvzUaPF7RdKH7g+rxKTgSZ+vWrXwtOB3nEW61Wthhy5Yt1rnyy4qCv1udKZEymmt8bJ+27UT6SSHNzMzcdNNNGLKo1WpIoVatVuFUF2LAUIMXw0bfKKRarVar1Y4dO/bRRx9hC/SEtQ3wnMRjr0FYS7bDKcNm8d1gOiHqDPsjRBsSwVeI4JmqPN6NMGufGAujDLuhNOywdetWoxXKcTc4lSSswkajwQtFpw0933zzzZ9//rlvZL9UiyYkbtu2je8GvFCQRAgXXFxcRND5Qw89ZO2ksbgQxOx5Ktg8QfaJ9BDlGIQC9wVr1eBXxIqCuxLjwdKer5GREbRYNEv2ZaKxXbx4kdtb7IVKrGHimhe80WNNLXLqeMWs0y+LavjjCemDEjBt/P333zfSSSgWX7mu8vcIgk2rfxpdYxEHSRUxfdMhzczMjI+PT09PY8jiG9/4BuKYrR12nEZGTyNEYem5wQvRv/RNh3To0KG//vWvZnby5Mmpqamvvvrqf/7nf8ysXq9PTExkHOj9ELt2cs6GWZEZ4koCqX1KNMP0zTfftLZCcgMHu8FgxM48eD0+Ps5BbjDTYntwy5YtMO7YUYTAPJTsbzg7EdQM3EKw7Hw5DN4Np4PZ6KWxscmH4HLOnj2Lj7hk2PLYGYYwhOOJEye4nJx3u4cd+tcK6bnBrxcZ2sjMbrzxRl5dBRri4sWL1tYQ5XIZre7KlStBCUH5aOR45bRVwBsqjsXpAMp3zYRjPdmxv+IQPHGjo6NxJmK8opyzZ8+yzsNXrJCCVK2J3iNLibLLaL3927Cz6ZsOCQPoZrZp0yZ4dPft2zc9PX3y5MlDhw6ta9WEWHvU4MUQ0jcdkvPKK6/gzfHjx+fm5h599FEfW0/ETY/YcE6cgbQinRSYP2Y2Pj6OSTxQDEiHCkXC+mb79u1slPFq5XgdHx9nx1Jsmrkgw0YUjoA33jI2NoY3UELQNBzO5wv0oWLQN3iN79XGjRv9jZ+II/e2bt2K+vuqFtZOUcF36f777+egu64xRbG3rzcLsR9H3lfa4FdDWuO3zpsffBWXwOVALkxOTrLySPRoxk2dZ/z448M+UTwLeHXFw4ocjZNHDvwsnJqEy8SQQ0YWVzhNcSGtVguz8djtyvcH5TebzUSPdTb92FxXSf91SEye4NecEy2d7BCGxLBv3m3nzp0YCkOeN3RFHl9gNHeVFX0cbI1/fOscZ+Aey0MeuHFzEm6c7vz58/iIZxJRBuiceKxveXkZDxj+NbAbSkblER1+5coVjptAyQhtQMc2OjqKEnh8EoEMGL5DYLp1zrflscqMEQy+873R78/21xDtnU1itxQ/MjxchpE6/NyJ5aAlW5JZFqcDL5fL3BWhBXK4hAewYCOvocyBCR5mzbvxIxmElXNN8B623d13342aICaI23MQTJTd9lYzUj1I9FPYtxBCiAGmvxXSKsnwfqeFRwfz1/AVTCTEdrdaLZhgiH7mAQEYiW6gcbCAz07lMlnusM3FRqgPUHBNWHCMjIxAwXBMOV+je5J5vIITrvD1LiwsYDfYlYj7gjby+HXUB3NmsbYTrtHXecLOGNiEW7vVOXcy7TfqiszMmAyVn03OEBIPs44nCWBoEQppfHwc7SQeVEAz27BhA0sKHr4Lqs2rcOHBiUPJjaIbrDOUKZjAztPP+UK8nHjpLx45v/POOzHcjekceCJYXflJ4yietWUwkglJIQkhhCgEQ6GQ8tgOpaT1kDLctjzQDGPKU5piHBlxAciaii3wLXkGVfhpYDbCTcqxBn52X2fIKxAYpNgNC7lyfVDOqVOn8AZnj8ffcciWLVvYdoM/CVWFAYvkRq1Wi6+CPUk+sA7fFVebfVEIc7j//vt57VpMrOFbmt/ht6Kg/CFkRXNdmYwYE35GvNmwRoeAhjZy5z/H70AZA3g3S6US5E62E9cbPIdFcE0WFxfTNFZwCTxCwA81X2MwApFYK1ZaeMUjEByyoobag9wZjBYuhSSEEKIQDIVCyh5GT5yzFr8PLBceIPZFYM3sq6++YlMRMWmemdTarqZyuYwAPPhsIKEwGO2WHQ+po0xOauKV56+CqX9mdv/992OCLSSLW5HWaTx6ctU4OJVPt2PHjjhVKw+vN5tNNg8BTseGrSdcQVVhV/JKg0EJGS49E187rCH8p+GQOTgIeW3J5eVlVvloJJgDCzyMDYocDw7jLZMnHgRfxdWLp8MHDYyfAlb8XlSaNvLScMkY3oDixxMHFYgrtRyiJ1EOxrsNNlJIQgghCsFQKKS0eKGutnasjVpJWeKD9CQQBMgligw6sJuwG+LcFhYWMC8H7pM4Os7nW/AsV8Bb3JqDNQpgPCKcz4fUcV68hxbB6Xz2BksWfOWj/9ZWctZp6HFwIDTT/Pw8ysFsQZiNbHLi5ly8eBFHwfGG+EOOwrL0X60rQxtQd52IG3wgnbndop1AIeE9fohGo8GRbGh17pg0s8uXL6MJ8UQ3lMm5gq5evYo37O8JpsomuoKs0z0Z+Ja4tQS7pWW38geQ41eRXwNXh+hWPAj+/PK08Yz/n6FtvVJIQgghCsFQKKSYPLMrrHPo2Q9h2w0qx6eaA8zZxrGcQRWvCECqVqswnTBfB7oBQ+oemRbH/8RWoSsbzlaJkj0zCkrgGCeetwH/ls8m8SHv4C7xjHevKpuuCPNzgchpVTn7AzRTvV6HamS1h1oFOim+5xlj8XFEWfDbSTz1DLc6NAYI6CtXrrDrhefVsRa5cuUKuxLxFaQ8Z1P1N2lhrhs2bIjHLQIxFAfgdf3FY4Xk22MdEygwfk5x4bgubMdzEaScyFmrrhXuOYSysEghCSGEKASD1sF2JbaREwVQ/FUQTQR78K677vKdfQ4QlAecQ4ipY23kOgDShDM8cpYEX1iPd8hQbHjlHHSg0WjgKET+YKSbl/sLsnUFIUYWjbZzjBPv4PXxbHXW9mxxKgdMaN+5cyeHBWLkndP4+zz5nu1c3lJKSrAr0ki0u7l5sO/Q2ncV2hfNO87dYO1fNtY0UNjLy8tpcZWJdYuFcqKqDuoQ75BzC3/Ft6JcLkP5YQd+suAh8weKhxPixCvW+T+TERg8qNoIDOZVdSVuzYnRxhleR/Q0eDIRPuCDXdgTf/34d+b/dPw7VyoVlv886IdyKpVKnI2YZ5t6fXAs+huugD8baP3I3IOdud/ynT1swSiG27dgJuzWrVvjCFp4odH3BAvUcjpwvMcNufnmm+Hy5WE6Hkj0oci0WPzEwbq0/kZDdisi6Bj4tgf/xUZ5rDmRfBxKEBQe/wRxHlUnCAHnOsSTZBuNBo8NciRF4uMc5+cO6hkPFPMV+eXHQ5q4FRi0xzC+1wejlEGOouw+ZkgarYbshBBCFIKhUEhpWjjD4k6TzKVSiaUPTCHW5gsLCzDQIDhYIsAswsCdO0t9FQnrNLI8XJu/Ym3EAsvaAgVaDd5mcOXKFY7hhmzicQx/5ZrwJXPYt9uDADuzTBwdHWVPdWzAopzPP/8c++M2QhQGyYTyBOtnb+x5NxHDQ7VBhAu+4tHgxKSofPNZ9PBYVkA8Pfbq1ausjTJURVyH4HlPG/8IyozbIQ8GWueTyK94WnFbXCHhK55QkZgra2ilvBSSEEKIQjAUCinb1g5mvGa4E/H+/vvvt/bQcJwV9Pz584i9xtxAyA6WLB6zEEdMcGpUdzKBxLVijYQUwAxTnv3qhSD+ghUbToTYirGxMY7AjmNqfUs8mM6u6VKpxBHhOAVcVgBi6Msvv4QSYkc3uyJGR0fZK8a/F/8oQ2hCfs2USiX2GHFsN36gUqnkMxksPfqm1WpxoiBuMwHxxNgAFhaJE2N553hpWtBoNOJ4h9i/VUpKHRS0vXhKBosz3BB3CXNIVKCQFHdjUkhCCCEKwlAopGyCyGAmFjHbt2+PB4vZTLt8+TLbR3iPZbtQDuzEarXKkWy8YkWpnY4I2oVdTZxc1WPw2CiD6YpyPJkQz3X9+OOPzey+++6zdqAgWFxcZNkUR9nx5VinkELlseWGG25gHxKfmmt15coVON6ghzgw/Y477sB9hsLj3ygm41cLDpSiWimB3yVNB/jSjrECCGaY8u8Sa6P8ccyJDiojfypXlQ9hCRLkBIrxy4+D6PyKjFpUWq3gYfX56Xjl++N3T43TpJCEEEIUhMFXSIlGdExgSvNGNhXvu+8+WO5YoTyeJfO9733vN7/5jbUFAYeuBRMDOfUkdsbUHB8EhxUJDZERERRrIx6hvnjxIj76iL9/xfNV3WXFbirGE2iyjYybwNfl85m48vAhcfrX7du3Y/UN1AFf8UnHxsb4itI8Comev9WE4Q0ePc+jDCbixClNfaHuOFqMp1cD9z5miN20mrBXyd9gI0/QjhdnCWAfbaz+g5oEM3DjtWD8qnm9lTQh5W+4khwuKIAUkhBCiEIwpJ1zxpyktBkJbs5g6gwEB9s4tVrNzM6ePcvrdsMgYnXi05UQlcTZTnnx71KpxLFMvHgE183XPYsXj8DOSHLqX3366adeh2BOkk9092N5uT8/CzYicQPfRjc/4/n5OFGwQmAcQwg+/PBDM7vrrrtQc6QaildRS/wY+/w0vaMHAzzQ/Sx6uFXDTVgul3n5lcTsVmmFB9kZAuGVp6ocuQfcWxMHcHJDdWJtHWTPSovTSxxTATwigidleXkZbwAvoQknE3tbV0rPIriADMI1dCVtVCf4CdPiwtFAEe09MjLCi6VyE0cIuLWndqIL4TDZ999/39odTKVSQRPEDjxLFA/A/Pw8AhB4Bi6efFQb/+9eBx5qS3zmMW7GA4B8vaOjo6gP5yviQTMvH0Ht2Mg9n58IDx57rbEzLtw7OeyPXMgonK9xx44dGMCMMzczwXhsxsCmYmpXSjBmxX/THDC9detWjLhyd9LV4ItDJFqdSXrSfimfKss7cLfUbDbT7JI4tMHf8AzcoIbxhTNBt8S9aWBFcS+FQHk0e7wixie+XfGJBhsN2QkhhCgEg6+QEgdt4q+CKZ9xumKMIJ0+fTqe2sbDbs1mE0YQLwwDVYGpqVAq5XIZqoiXW2VuvvlmDxOwFC8xKoDIggcffNA6JYvHMsRpvBPziKPaqGqcAQU7f/XVVwhD4DSvwEUVZgGzToIsYxPSc5CzQc0W5TvvvANpCPvxnXfesaSBuyCYmAnU0hAO1vVAHKRTLpd5IxoJpjF4nAsPQLEWSRyRCwq3SBKxQGFFgi1Xr15NG8nAA+UFpo2L+JY4fVHiUrOJs2X5ctIWqA3mhKTVyh9PjrmIyTkK2u9IIQkhhCgEfda1wjRDOICZ1Wq1arW6efPm7KNWFG/KogdgUVcM8j7yyCOw1uMJbi5B2L3EaVdgwbkagLEJPYHC4WVxR1EcY8qv2G1hYQFn/+1vf2tmTz75pHWOlbOICUpgnfTZZ59xcDZsXl8Iw28Ixr79WE6K6uBaOLg20WHAh2C3oJw777zTzN566y1LsnODeYtpft1+F0a9NfgeiOWLv8YzTFmqlstlDAMEusFSwgfi3fykiT7C4P2WLVtYD8XC5cqVKyx34uAIjxdPm+EQXH5c7cSRlYzrMrOxsTEMKrCnlv8iNm7cmFZCzrvUv42c6acO6ciRI2ZWr9enpqb2799/4MCBarV65syZn/zkJ7t3717v2gmxxqjBi2Gjbzqkubk5Mzt69OjS0tIvfvGLTZs2VavVo0eP1uv1I0eOZDyfPgoM0uyIxCAfXhb2kUceMbMXXnghXjEP8WC+6gS8LHCBnD171nfGjFEY/kF4NKfq8YmxSB1AUzA7AAAgAElEQVSEjyg8znzaarUQ1PfJJ5+Y2czMjJk9/vjjfKVxqkfOUQRZ9vLLLyMxK1xl8Ti+r6QXBK8b+RJwCtiAyLfPi/LxUrDNZpOD2lEOBxZXKhVELWIjRCrWGEyLF0/8cf0u9Z392HODXz2BO8RXdLWktlSpVPiJ4CbKhwTRaD0Ho4+OjnL751hqT++bGN5tnU6jSqUC1yZL+YxYc3bEBjNh47jwWLKMjY2hqXNQH4v7kZERHm7JaK6DpIdi+saH9NFHH910002zs7N//OMfjx07VqvVsO51tVr1eGshBgY1eDGE9I1CMrPXX3/9iSeeuHDhwoEDByYmJnxg3ed+JuJWT1dPUuyHgBUPN08QLcPhYfiDgJTZtm0bwt4wTxZgKQpYcxj7bjabLAh4LhFoNBpBSh7r1G0+JI2Pjz32mJm999571vYnYUswsZGHsJG858033zSzhx56CIKM/UaAnWruMEAJ+IvEZFuXQdif51TxSWGWWqc+g9sJStTNRnyFU+CHwNXFs0wCMiRU4hB8Yemtwa+UtCkv+O02btyY5st0YRSvPs5+Ry8zntadUY20UDojZWbtpsXiO5Bi8cxZPM4+Ey4jh5DRwAlXPjF6MHui25YtW+JEQRlRvtlVStwyGMqpnzqkvXv37t+/38z27dtnbX+v0YzURPrrP0hcb/qoPfTW4IXoX/qmQ7rpppsuXLiA90tLS5OTk5Ag9Xp9YmIi+9hsmzowqQDMmampKWsb+0hp44lJ2JUCLw6SDrRaLfiKcBQsdGgI944YZWoAsLkgrXziEaYE8TR4Ns1cP7FhiJ1R1V//+tdm9uCDD/IyE/DW4B8NCgm+MTcYeT45ewX8LKgV6oMTYQu47bbb2CvGSouD6Hy2Ft/5IIE/VBHi+pBACO43VDv++RLpX7NxNQ0+J7HbI3EfbnXshsFPMzo6yj9fHPbmAiuWCDzMkGEoJE4pY8dVV43CcXddU5oGbSY+RXAsX3jauumVSgXtmacDBiGvnCg2W28l0o+NPKZvOqT9+/fv27dvenr6zJkzTzzxxNTU1MzMzPT09MmTJw8dOrTetRNijVGDF0NI33RIZnb8+PG5ublHH30Ug+nBxwzSouyCgJ94wgFKhmmD6UcjIyNsMGLuDsZSkCvh9OnTyHqHEiCeUAKC8dxQiuORoAB8xQqskg448R0v1exLoQOO1kNM2rvvvourgAcL+fGwA5RTsMg6Z6Dgkl0/YQdcOBw8PJo/Pz8fBytCMmIL7tUXX3zBC/Rx9odgGhb0ENxLvIXvXgbBDv1lRfbc4JkMjZj2K7OyKZfLsUMIO2AM4MYbb4yfr3h18MQZY8EyDdkBZn44R8TxOg5oOb7yHpfDCYK9AnF8Xfzn4MQRgxlLjyeWg4uFrPSMq0Z/OPHNz8gpPKj0U4dkZkG06yqDXzN+Y/xZo3EgHtofbG6U+NMHeCqq1SpaP3cbIIhf4JVVOdVKkH60RfNb44RGlUqFs2jjFOjJ8C+/sLCAFV1RE/QHCEzAdYFms8lLN3EO1iCPOIrlbDF4xrgLt86+M844uW3bNu7GcCJen6nVavFXDz30kJk9++yz1jlAGpDmMc7+qsisvsGv6GJ5NAyNxGdN4IfD2BfMFx9bjuMdeCZA4ogc/+0GofnxP3uwJRjpCqodRAfww4jtQV7tjE6It6TVKvhDSPMIlEolTljMs3r9QnhEMV51OqhSH7XhFdE3Yd9CCCEGmz5TSD2Qc4pZoKwxHAFee+0136FSqcRpPxAwjTDlu+++G/YjTCHOMM9lLi8vs/0Iqw1Wv+dY5ExCOB0bdz5TNW3KHqzUy5cv4xSYqepDcxbZpzivj3v4V540FmfhCscpkUZGRlj6sCqCFvRFmNiWhG7z9XWMhiLxESt38DRbyL6cxmMfRdatF7EW8eaBmxxny0WzL5fL3CyxM5oENJYTjARaJDXi+sSHW5JY4YE7rww/m54aOC4/1jRdZ1sHO6dNpw1uJstKtN5gdAT3hPMw9VCffkcKSQghRCEYfIWUaFjxlsA9DkEAw8oT2lvbO1Kv11EgvPowxxA5jdjuZrMJnwqMIJhF7E8K6gPHCU7Ehl6z2UxbBJYl3fj4uPuKrNP89APZkQOjDG4YrkapVILtBlmG8PE4l4kLIPYrIFjDvWVxJhVswakDbzZ76XDH/EL4WhBtjBqyV2B5eZnvRh5bVSSSEQbNCgbtGc3D8/ayXI4PSTxFxi+Spo1c/vJHF+5G2gJfwUfLU7N5/CDxRLEfyEcy0khcnyIu06OBcOswXBF4ktgTxl/1qfuzN6SQhBBCFILBV0jBYHEcVOP+EhhTyMuCjd/97nfN7Mc//rGZ3XHHHUZZQRFfByGFMqEtRkdHWVdh0ihKgy5xJw3yrsYr1LnRx+6cuNo+Vs7xr3EAUqs9IZGjpAAbj+Pj43gDC46zFgWlIYY7DuH1Rc15fUKAK4Jx7SGCKBCng8bi7JMPPvigr9FubYuSY8QRSXjq1CmOOeSo98QffbAjlFYKS400XetbuKFyHL//lGkKKcM5lDM3Lr/3KLs4EVGcwcjIk2opWYLiWb3BchtBrGBi3TKq7Vv4LnEAut8E/gruLn6+hqTRSiEJIYQoBIOvkCzF0o+/hamCGUiI+4K1wjlv3PCHJ+Pb3/62tTPo7N2718yWl5dh7/MpkFyVh9p99gY+YvEIzDOF6hobG2OTCuqBnUO+/kXs6cEVoRpXr16NrT82XX2Fcl78gt1dwZwSyBRWV2zNJXpxOPLKd0OxkEFwyPki8SiZo/VwOmRyevHFF60tWE+ePMnWKFc1Z+zW0OK+Fm4SHOgVRNnxrYNLz2V0rLFYdgTrVsSCw5tEmkJK/NVYE6clN3L4UfL6pLXYxDYTy69AsQVzoYKTOjzPD80Vgyi+ikfiDK3gpGn3ZACQQhJCCFEIhkIhMTwUzokPyuUyfDzYiPQ/zzzzjLU1BOxBX42Gc/nwugkjIyPQQxA98TC9R6NxxgG4Uth6KpfLPK+CUyfEDh6vCYA2ctMMzhsWUjxJyPWWL59h0Wx2IzuaQ5Vg3LHTqFqtxhOM2DcGzp8/j/1xOj4FNNP8/Dz0InbgZEiYh4QK3H777UGuVcu0KIdqOH5FsDmfOEEtdqL49tihGLhhjJREGokhebGQSvxBg7lNrGk4KDTOx5OftIC3wAPEQwtBbmIuASG78Ib6MxIrqkQXWg+V7yOkkIQQQhSCwVdIpaQkWiBwOfAkanwF2x/bEVlUr9exTBznJOXk9pVKBatI8OJ18TC0TzuAToITBVsQfXfLLbd4UgP/Cu8hGtzgYv2BV1QA769eveo5HazTzsLce5+uhK/grOK8dmxXltvrE0Jv4Z7A5YBalctlNqhxA3ETOARuZGSEs/FD/yGUziUa6sCrm0P5cTl33HEHFFKaAJIPKY1gOg48l4krROAHxew6tEnO99HqTOaWqI2MXDvcfoLFGuK1GxhXS7F4YkXicXFpUX9+SEYoIB8YnyjIx8qpUnigIvCQxcMkHBkbXAXuMD/agctqUBmKDimO7IyHAiqVCi9yilE4DEnhDxfzOiuVCj+66L2wxRsZnPPoV+B75/9oh7sQ9sSW2gl7uD/gJK04BH/oi4uL8d8B5z/1rKk4BQ8LYDeMHlQqFY5xwLhZkB3Z6HnjjTyG5mN9vP+9995rnUMcnoIIsxcB6omb+dlnn/FUX05rhG4J46J+VDDnMbjbCvtOhH8jTv3Jt67RaMBwiYeSgrG4OFyI3wdphOK+J6fpEMRh884cI+Bn5PngPETvE8/jEbzEQTM+BX8V2LtpfWoQ+xAPkCbWgW9aPIo+kGjITgghRCEYfIUUJP9Ii6RstSeQwjDnjKKQz7BWPHUQKwmIDLehYN1gcSMe+oMtj0MqlQqcmZAIXB/fwuZYPPvVy8RwHy9TFA9QWHvIBSHauBYMkSHM3YVUPGTHLCws8IgZv3pwPMeA8AoX8Qq51qn2AuM6DjzhK/KgEujREydOWA6kjRwfOQAcpOMzASwleWgwGAX9ATLGSxM1a378108LqoaS87S8HBeOQXWMMPtjlaanE7fEOwfiKQ4iD7wD8UAiHx4kwQIc1J44qzcuJ76Q/kIKSQghRCHo7+40J7HVlmhH8G4cWo1X2F+Li4tQNnCcwKMOPzwmfm7cuJGXmYhjT10N8NqRrLQQDbG8vAzNga8SJR3vBm2E+nA8t39kTw+sLbcozWx8fJz9NDgFi0UvE/WB7QbTmOXL1atXed5u4DcyMvTYM8x5ML2G8BXhbrBdicX64KJrNpsIiOA7E/+agvG7xLeLI1yC/Ln4jSAveIdAOXHGLHYNBnBMeTzlNk+1UQKHGvEleJAFt7p4PnipncsnQ9LxxthxFVcsKCFxbRR+z/ELiU4mLtN/lLSb1u/aCEghCSGEKASD0Kn2QKJO4oAufAWrhBOknj59GhIBO8ABg9RB2HnPnj2xlcTBoH4uNtBwLEw8MDY2xktOcGlct0qlAlcQjoUHCFFtnu2UnUO83B+2IEa8XC7z4g6xI82FHS+HwalafU02FkzxjfXcS3zz3a9mtMYE1B7CHfnCUTfYlY1GA7KSVaO0URqBkc7GeJoE9+3B5FMjWx4f0zxJwanTZrxmHxVv5GBrFg2tzpm8aKLxFbmTkk+dKIPiCifWJ658IMj4Iw+cuGuZS2BvcTCK0NvE3n5BCkkIIUQhGHyFVEpfvjrDcofsSAwzYyMRoXRYCgEphXzMPZ4ox6U1m012DuErGPse7xQbsFyy74aNtVot+Arb77jjDrzhibrIUvrWW2/5zpcvX2bpwyPvnG0omNCDQ6CNfOAeTh3cT15Rgm+Fr1JRprX7YMy6xwK+Op4Yi9MhUNBvCzby5F+uYWKg19BKqAwHamLSIIsmtLIPMng6eN4PjmJ14j4brklArGNiEeOF4ERoJBx36quWYE+MHPDyE96S05L0xFVKrENilF3aHQ4UEk+rx/t6vR4vpIJr4Qaf8W82GEghCSGEKASDr5ACghA1f4+YNwcxaQjlwlfwJLlVFcsXqKWxsTG2kmJXSqmdH8iztfpX0BbBgWxtcXQTbKj5+Xkci6oizQHw1TSWaXnAP/3pT74zr6bcaDQQsMTChU1gN9zYswV/T7CCNQrHdeGKOOoJ5dTr9Xjda+yAcMFKpYK0s2xXcoAfVOlbb73FE2jSjMfBCEBaW/yepLk5/Wby5KTE1Dj8G3HOEbwPTpRYhzTivK7B9Ck0icR1LngHfqb8IYqbN9cqYx5SouzmuEHGl92LFRKLe+ts6oEIs7YPKb8DqU+F1NA9qLGm9omu+LEx7IO/RYxCoCvC6FOpVOLJnhzJ7amD4j6P/3bxh759+3beAX0DD7h5xXh4hEv2b7E/ek0+EToDX3kWV3T33Xdb+ynFdoxplEqlOJaXt/hQXjwIifRIGAy8ePEiVzsehfDuH/cZOdF51BQ7jI2N4T7jtqPa3Ie5Txun4+VlY1PAPdh9+pReD4IRpxgP0Od8g5zFzg/n9DwxPB4bEP/tcrHW7tJ4VoB1RgSgeQQWTzy6yOPn3nPEPXH8lAX1SYx9yBhXdCqVSjxdlzPaWeejxyEbGXeJGYy4Hg3ZCSGEKASDr5AyJrsB19T4CKsE6yH97ne/s05Ta3l5GXvCYIwDAdyu5MEEniuaIcggWVxJsLriEUKOLxgZGUHcNltbOBCX4B9ROF5hV6KG5faStTiWU8fyheN6FxcX44mNeA+RNzo6io/xUp44qatJDCH6bbe29Q09+vHHH7PFzbeC61MqlXAijgeJh1+Cj/1rP64hrc481lAhfLtcdqRNIA0UCYYTMgz5rlHd8bhF1zjsYA6vURg6EzzjgGN/0gITPIggvi5v/GkR88G5EkPPrf2k3HDDDWkxHXxjfcQyLcq8a00KjhSSEEKIQtBnCqlWq1Wr1c2bNyd+zCZtHpz7SGLjhUeEIWWazSbb/lh+grPE+4RWFivYjUer5+fneXwcYansPU5MAYmNvODC8vLyxx9/bG1hAdkBbQTXjl8XePjhh/1asLCTh2xAMKEEDh/gZKmJ0avwHmE3V364txxxwMastbUa+41Ye23btg1voP9QAu4V7puH8OLCIZjSVo/t32DZ1TT4rqRF3wR+CNx2bvaJS69mzHi1qOXEBC6ZtNWMLEnZcJMLasXH8sJC7jRlVR1LtLSNQfmJN4G/8lrFl+8eMr6W+Fb4dv44eFE8faOQlpaW9u3bNzs7i1czO3DgwOzs7MGDB+fm5ta7dkKsMWrwYgjpm470j3/8465duw4fPry0tPTUU09t2rSpWq0ePXq0Xq8fOXJk9+7dGcfGwWmxBTE6OgrTAxNIMRmNo9HcQoFJxYspsA746quvOAwvdu24HefrWVhbpgQBdbGfhs1G1O3cuXM8Mg5HjmsjIxuK7cHHH3/c2joJr1u3boXQQfB6WiSrJ1fl0gBUYLVa5YvlWDsujZMkeSVxE9wFhbBvtgSD1ErYmTMbJVrufcpqGnxO+CmA7mQ94T48bgwcrhYIl7TA6IxoNKbVmcuHDwncqPFM7TRdYinaCJeAqR3I4RuvvBd8TIsL96j3NAcYCELp4ltaLpdxRVxV1qNxRq6BpG86pImJidnZ2bm5uTNnzkxMTNRqNfQc1WoVKRKEGCTU4MUQ0jcdUrVa3bRp04kTJ86cObN///4LFy4gsMrMsERbBuyiiO04t1Ngg0AiICSMvT6e3B77v/nmm9ZegI5nI918880w7ffs2WOdy1gwjUYDR8H5gdPxChFeT3Yp8RQfaJRPPvmE18fDSQOfTSuaVYpbgQRCr776Ku8GweSr7QWVbzQa+IqzsADXN7ABcRN4GXK+5z5/iO8wLx54/vx5NoFjmXj06FEjV1NMv6ul1TT4nGRHanmUHYsMfpRcCqf5n4KwsXg+TWIAW1wax8hVKhXUxHMHW/TIMHED8EELtDRuZnE648RiA29ZWguMQ3kTPwYrDSaukp5RmQGjbzqkl156aWpqav/+/Wb2wAMP7N27F4Nd1k4il0H/uvjE0LKaBi9En9Kv/9STk5NIJ1qv1ycmJrJ3TjOggkFe2FzwhXDiAzb03JBBJFuQQN4oDA9Ti9iVwjMnKpUKz1uCTmLpY53RfcEkJzP74IMPjCbiPPLII9b2IeFAX9qcg+UArEIIsqefftrMnn32WVQPa4FjI2tKz/bIoocnTuHw+fl5SElcOBQkvsJ2n+rvitOitJj80/AAOi8eiCttNBrBdPfg9wIDYFeuqMFnEzwOaVFkgUJKnKNj5NSJEw0Ep0v7FTIiITmFh+cuShN2/runzfgJ5v2gOXFq4HiZ8FbnbK1Y1hs9sJauyH25GS6HfcZ+xxLdb0ZCKp5sN0j0TYf06KOPHjhw4MKFC2fOnNm7d+/U1NTMzMz09PTJkycPHTq03rUTYo1RgxdDSN90SNVq9b/+67/m5uYeffRRDKYfP36cP6bhVka2YfXFF18gXyc+clp7TsPlYV0wr5DXgO0dX5+Cw8Ng1PDieI1Gg1OEYefbbrvN2v4kfxMnVEVcEILizOyHP/yhlwO4hv6RE+WxQkLljxw5Mj097SdCxlVkkAsWEmS/EXuncODFixehzOCKixci4xvib3AtOB0k1KVLlyBV2WHApujBgwfN7F/+5V/YnxET/Pp9NH295wafk8D25xBHlqqJpn0wBs6hd9zs4+hWi/RZTJpiczdP4uIXccnxsSyD4LC0JM0XNJK0uUfBehDxdbHL011WrPiDrLVxUCuXELjiut7GPqVvOiQQRLvmDH5NU8H8rf9FcpLpxPV+uL1i3AkDU5hSunXrVkzP5BJ8OMu3j46OIrs2NyzUAYNRrfZKP4C7JXSWKPOxxx7jqPQ4Qap1LvnKOZs5Wt3Mjh07FlQybvpByei3eEXdHTt24O+SAz14drCPXmIHdEX33XefRf9c6Jl4xJL/8nBq9KnW6Y7mLtPpo66I6a3B5yT+U0v8lUH8Zx3czDhxOx+YmNEn8RTZZbZarTjdcDDSGM8njTukYIXWuAIZQeTB6GXakCZbYzfccEO8xhjwbglf4d8mxkuOO/5Bom8mxgohhBhs+kwh9UCikzPDXobp8fvf/97MfvCDH1iKkOdBLXZLViqVnTt3+p4cjBCs2cOKAaqIs+BY0jACXjG0iCRAHkPBlWdB5vHlkHEYT+Moc4zP+HotnC42tsISYzGCcF5s5KMgWTBkgdN98cUX2IgbiGMhrQIJxeMYvAVBJS+//DJuXdrsxQyH+dASqArcWwxhQZUGEySCESovx0Ou05R0IDjitY8zgvKDGCKjjEHZ6yolDidydDgOX1xcjFN88am9tLT/Co4UDzZyoA3GQj0WA1fBi8DiuahWqxiKx59A2imQZ8hShgEGACkkIYQQhWDwFZJ1mwAIC2hhYYGlxqlTp6xTvnhYAewansTKwaDu+0lzk7rVA7c/4JUgPCtPHIPLXiUPecAbDonmhD1m9tprr1l7IVc4uuJQglarxXN4OZKbty8vL+MjagIxxAHlbnhycCpHeEOiNZtNRHBAFeEO40DsfP78eZ4wyLcOFqUHBAerGib+xAIEfv7sGZ2uUTizLbfkIGFonNsmGI2IA08C8REfC4KS44c3uJAgKj2osL+PBwASHcxpDcnvEj/vAE8E51xutVeT4WsJXE1pM/fT6hYwGPNnpZCEEEIUgqFQSGmjwDx87Ik7eaw8tr4R5+1ggBiRYAh+27BhA+JkODErSoM68fHl2ECD4PCZp7HjhLd4ABsbaBiGdreQmTWbTdQBPiSOn+YrdTnFgWqcssVD8lg2xbrt3Llz7FHjRQjx6rHmKDb2SbgSjVcH4PnFHsrIkcppcVP9azBeP9z2jwUQcPs9nuUauIiyZyJnOIqC8YM0FyCXmeHUAeVyuauYYNKi9TLaTDAuwnqIQ+m40TYaDU5fy/fcy4z9aon6jxeuHLzmLYUkhBCiEAyFQkqLrOPJKx6uFqsi1hDlcplHgRESA9vfhQvbRwjGg7sI8sW9Lxz/E88lss7gPZYsbHlVKhU28dj341E9ECXY7a677rJOfcO3yN/Ek1iDAW6cgrWRh9LBp4V8awgIRAUQfBgoLb5wft9qtT755BPrjLPiqCSfvQEzkyfnxj+9AIk+CWyErA8SfRpFx3HkGEtn3yFuTkHm09j2D4j1WWLwJD/OideYKCyyT5oRpxfLpmA5jODs1inrIWhc/8WVD6L+uJzAV8flJF544vb+QgpJCCFEIRh8hVTqtnayGzUcToZElogBg3sJptDZs2chdxCixhacGzjxRHGUjLgy0Gg0EFrG+YQCRcLKjLOv8kkXFxfZe8TLRuDU8/PzHO2D3e69917rVIE+v50nZ8RxSpZkJMZ+KWtrIyhI3EY2qx0cy3feX2Gzc1Af25XuacNGzwQj8hP/vrGE2rhxo6/7YJ0Nw0cUeFAhLsFdI7HTNFGx5ZQsaf5CdzJlX3V8VNoOaadOlJj4r4ijaoPaxvrv8uXLcXbaOCeFP18r8pP1EYN5VUwwyhQP3/l7jsvE08VzPD0nFR4/DEPxNFJs9xRz6H64kaFMpCYaHR3Femuc4TsQ8vjICybFeduWl5fxf8Hqnv9iPD0wD0h6ujyjvxiOMuAngQc2gwET7i38EF6vFjeBxyK8n8ZgJj/AzKeffhqc19qDdagnwsdLpRJnS0oM2zXRSfDHnTZQFsyEha0Tz5sOkt2lnS4xS/eKuoHE2Ad+Nr0avH/aOGH+/i+tPomx7/EOPgbOlUR6TB7xTjwpP4D+FzHYTVpDdkIIIQrB4Cskywwh5e0Y9sFG+OT//d//3dpTSqEkbrnlFrZlYLwg+BvcdtttGE0CHB2AkSvWPf4VD3e4DxnlJJqrXG0e7+LSoNVGRkYQb8rl8MCLv7LFFy+dGWdYCTZ6tXkkkLUjLgfi7NKlS/gInQTrG1t8YBNR8gjv5qFIKKTPPvuMqx07hAMGL0B2lQRKIs3l7ookHj/wLfHiraxOgsHexJoEbzLeB8S6JEP6pJ03g677BNEKIM6PvnXr1njJrnhcztJHUDNON0hIIQkhhCgEg6+QArcHv2c7rlwuw0KHdmGTHHj4NRtNsHrYAxTkR0G2Sk6N4zsEAsU69cfy8jJqAvOT/ca8jGaz2cSx7GRyFYLrxUcMW6MEOGA8aZDR6gDxveJT+27xlFiPXuVoDpSDAG7Oq99oNHjlWV6fCVOPvWTI01gmuoaLFVIig2pR9kyrncsnTj/Kr1evXuXAnDj/jVHYghee6NXnphWnycmpbIKo7vhxDgRZV53UtWHE15JYT/6IZxPPvt/kIMQjrdqYxhD8BP5+eXk5bX2KxLvUd0ghCSGEKAR92YuunsTgVDapYL/DxgFwFF26dAk7cGIbGIxwEY2MjGAGKOBAbfYnBTGjLHr8QJ6Ryn4azp1aKpU+/PBDa8+/g3kFowyVv3z5sgesW1v5YQfWXiMjI1zV2IJjt5OfHUrLVxTEgXE6SwR/870ql8sIMmQ7jlfmXV5efvXVV/0oXBdKRkS+O5yyPRN+t6WQnEBD4Laz+uRfn38USzLqgzCzNI9UqT1/NqiDpYiM2F3qpbHG4ufCK58dPZhYh6AmacQ7BIqNg+h4Z/97iZWWa6a0+iC4NBCaK61nvyCFJIQQohD0a0ean4wgn2ALBAp0A2K92MSDwhgdHYXIYHkRWD2YfIqhXugqNljcsuOZhr6ssrWFVKVSiSULL0aOct5//322ZBEoiDyqvh17Yj4vasWrpwfOGA7Di/0K1o6L8xSowZaxsbF4ohLcQp6+FiVD9GB8nGfO+gQj3GdWjVBaKA0Cq1wu80ymtGXWhlkeJWoRI1sbPwG7LrjhBalxEsuJNyYekuaMiUP7Mk7ky2HwaEGsmaxTFSWelDdmTKWKnV5dIzZxS4OQVI45TMxXGz+AAH8I+HPw+bOD2qSlkOSIyMsAACAASURBVIQQQhSCwVdITnaKFI+ygzb6+OOPzeyDDz7wnWGwb9myBSY8spRyaVAn9Xq9RZOB8AorHtl0fBUGXliPa+JlxgP6QXoegPNypB8WwsD7jRs3osIA2gj2F7ZDM+3cudPXswhOyrF2i4uLfLExnjkCx/Lahpw7tdTOsODr9Vl7TXeU/8ILL+Ayf/azn/klYwf+NTdu3PjrX/868QYK0DXAjD1JHNDliyMEy4JYp5nvkoVdSnHoY54fJW2fRJGXNuDh2aS4McRJY91HG4uw4I7FPtRgB44/5MsPEmKl/RB4AFvtleDjZQwB/kOazWb8DxajKDshhBBiVfRlL7pSsmcSuGUHLxGMep4DBE0DP8e5c+fwhicVoQTshkOs7eTAV9gISYTyFxcXOegOx8KkgsIYHx+PB/FjKyyYiAOXFbuvWu0cfRwxiPrA5nIvF+rAVwQ4/V2z2cRRwZQpI6MYHizoIZaJPGmp2Wwipx8veMhrefjOiB7EV3v27LH26vK4LQsLC5wYIi1OSZopkdhjFOOL1rP/ksW057Ljm8yZPhIt+pzTg+JIucQSMkpLi9lL3D/tLyLjpBs2bGB3adcFUGL1xskvrNMTxv8/QUKHxMLTqt1H9Gu914TAyckpR9EhIT+Np/kxs3PnzuEPGsHWLOfxetttt/F/Ot5joIwTU/qCsGht+JdHV+crtMatKvasTk5OchPHvzyf6PTp0xx7jZqgx+IuykuIUwtzaMPo6CgeEg7WwLOE0R6vAw+VoEwecFtYWMBlcvdT6pzxiqpiIyIyeBoyKv/WW291HbvI+HZ4SByMiv9/OcjbPfNsE2BAGzYHhm2bzWacaCojojr+T88IO4o7Eo/q5p05JMEfnLSghoz/9GAhqHi3OIZiw4YNeBYQnI0nC++DvpDtNs655fOO+dHDM4V/JDw4wVjloDZsDdkJIYQoBEOhkNg8ZIsp0UCDzY6BNRhByFIK23BxcRGDUWzvI8DBC4E3GCXwSdk083JgUrHSAqX2XM54eQU29EZGRjjZKI+AYeexsTGYdZyqh7PBui5hYcfTfnEiX5CJr4vDPbDbuXPnUB8kRcVN46msHrPLKZp4UAJB8wsLC/fff7+1XesY/5yZmbHOQHBcdXAtOWNzBeCgBpYFQewywk+Q8RbNLJjRyQogFmStzgmkacNlAfEIWKlzebOMsam0rwKtli3LLElXJboAsBHCPSMVsg8kWDv4CPFTLuywJ+42VCmXMPDtWQpJCCFEISi0QlpaWvrrX/8Kzwqo1WrVanXz5s2JH7NJixN1w4QDEGCGYwz34MGD1jbtX3vtNQgmWO7xwno7duyIc+cAHnr2ZRriyaceLhHPSYxtOo+s9cxG1jmE7RKKwzRwXdjukoWDDljNBL6lYHlA67QHb7nlFhQObfTHP/7RzO677z4zQ0YlOCS82jgWMhHayG8mqsfrxvKvhtCG4DcNnIKJd6zIrG2DBxlKglsdr1GCV9/CNnucLMo6m3qGmz1NN2SsvJAz6oEPCWqVNi6S6NNKhIPX4539PR7AOGjba4WWz34jzCjH3JKrV6+yyxbiiZ+vxJUAB49CK6Q33njj9ddf948HDhyYnZ09ePDg3Nxc/FGIfkcNXgw5xVVI+/btO3PmzPe//318PHHiRLVaPXr0aL1eP3LkyNLSEn/cvXt3WjmlpOUn/CuLLA6Y3nBaQFW89tprZvbDH/7QzHbu3Pnuu+9aZ6ARTHuYjfV6HePIbK9xshNolCADCgsp99/Ey99x6Kf7uuCGgekK3capWn2mKmqFS37//feNtJpR1B/fmTgg2HUkw/6tkZER3DTE0N97773WKaEQO+eD45BTvJgFrmtsbAx74lgsYIEbAlP0rbfein+72BAOKLJgWqsGH5AzAhgalBOtBncJTR0hlFgoxOd08xQIPl3XJDc8ZtD1p0z87RJ34I0ZkiU+NrH9xHoo2CFt/i9LIp+ui9ELjur24F4eIOHTxR7oRLpqyr6guArp+PHj/nCaWa1WQydRrVZPnjwZfFy3WgqxRqjBC9FP3amPre/atSv+mEaGgZZolPFyZPEs0UqlAoMO3hqEG7HHZWRkBLYkz3KFIED5OMQtJp6vwwZRYNnxVzgdwvyazeYdd9xhbX3mmUmtbX/52D0nIX3kkUe8ZNi5bjDG5mowNB9PiWU715efgI7hndnx4PN5ca9gJGJIHVt27tyJm4YZSFCl0F64cM+HlBjyFFQ+7WPB6a3BrxR25KCRBCmD2e2B247gRk/yy/Pw2EJnQZ/o2kkk/kE5p06g0TN+8Xi3YF557BBKjLjLCJkD/BTw6nlBBXgJFU8zZu0/h5GRESgqLjPxXmWr/H7XScVVSDH4/cwMFmLwUYgBQw1eDBt905FOTk5iOku9Xp+YmAg+rr58dzXBhIG1/u1vf9vaWue9994zsz179nB8F89R9/wIHK3Hi+zBooRzZefOnTCCeOYQi7OFhQVO1gJgJHJ03KVLl3AfMKwfu53K5TJqghOx48pTQhjllnXBZEnqrdFoYIIRV5gTTIyOjrIHq0UTlTgnxcaNG30Jdi8NW2A8btu2DcY4MtLi7LAxEbmX4XJI+32zdygaq2/w2cayS4Q4HJQPD1aWgzxFK4JQGxsbixMc8KuvJ5m4mB4fyK039vX6gfEy6olOYn6aYvdM2sX6ifzUcUZj3tmf9zgiDgd6Bi8uAW3b04zxrYiXUIkneGXQv9oI9E3tp6amZmZmpqenT548eejQoeDjetdOiDVGDV4MIYXukA4fPswfjx8/Pjc39+ijj8IkDz6uiAx7GfYIdAAm0MByRyjdnj17EF934sQJM/vVr35lZi+99JLvMDo6CgmFcBpe6Q4nReCZZ0Zgw4q3XLp0CZICsJ7gOeFffvklO6jg9/ZEqDiW0+BzWkzU0zOfIqoNG3FLOYkfaDabbPbC0GP3gy88yKoRYAcE+E1OTvICFiw6PWEdjsUPgZPiEN6S6DBIDKfsC220tg1+RcYyi4x47TvrFBloYPj1R0dH0XLw68Q/ROCqiX8Ib/xxUFms4UZGRthZFYueIHwuPmnQPOLJT4Hvh8/O/l1/jaPs4nK82lBFaMP4x/Dapq38MiTaCPTZNQTRrvmDXwO6er/RTHkCqT+QTz75pLU7pFdffdXafQP6j3PnzmGUCcRrvKKjqlQqKJBH5ACarD8JHDzKkQ6eeQj/C/gKk3YxfAfq9Tr7lvEKpzR2Rq+zvLyM+mMHVDJYwQjVxle4ZAwDYhDSR0hwOlSbl6bFIejjfaFbfPXQQw9ZO7weARePPfYYQr0xbZCvritdf9w+Yq0afExGWLyvhGQ0JAXwHrYFXFmVSoUzUWX8e6ZZgUGHxBvjcbnl5WUeuGaCKAwQT10PJlrEeU4zYhbwynNgW60W7kac7Zur7f0WHha0ZDgFvUPKGFxN3J6xc/92Tv0U1CCEEGKA6deOND+lzmyM2Xv6ex7UgvMW37711luYIQuj/p133jGzqakpa0c93HLLLbzyLERGPJ/O11/guvHysrVaDXHSvEhE7DjdunUrjoXNhRDwwEzjJVx5yIUXZHL/KufywVfYzYPIMSrIFiXKx1jcfffdx0oLp0PoNg/F3HzzzQgbwf158803rT3XGKW99tprvAOAWuq7Ga99BItyn2TNEgH4yiZGyoaz7ib64YPYBOsc+/KZALwzj855UEC8MhCf1Jfs4uRecQiSny4tIsOJFRLuTzBkF69bxu8XFxchKPGPgdiQIJAkLqGH9ty/2ghIIQkhhCgE/d2d5qFrfHCQcYSBewYKCQZRo9FA2DEyrsJIhO3z9NNPm9kf//jHOLyVTwRGR0d5pUhehtXXcYDGitOrIPzXzVhYgpyGhG1D6xz3RwQETgpXU+C7hp3LCU64nvjW6+NpkLy0zz//HBKHc9TyCoFuq7KTAKIKohNy880334Q2Qgnw1fHlBJldpI3WBP7FE3PnAFbVn376KcfLcBgLH+65qTjaBc8UmmupVEqbVcrqJFgimQNh/BmJJ+Tyk+UPJiskFitBhHc8+zuoG+cE4ueC79vZs2cxgoIp3omLWcQRJbH3Ljgqgz59IqSQhBBCFILBV0gZ8aaBfEmLzkIeT0R/VSoV5JeEtwZGPex3OEKmpqYQgAfiqaPA5wlyelaOH/UJs3HAN+xBCI5arYavoCdi59Bye8lnOJlii9IzY7IPgMMCof/cfOPQO+yGYCHsXKvVOOsSyxefIGlmp06dwgUi3g9aEFtgXZ4/fx4nQqwdp0QKfqY+tQSLSawnfGZobLN72BjENxo5XIzxqpKuJPAVdubdvPEnxl57OTfeeCOvSIIGj6aCeNFKpcIeUNY0rLRc+sSh1UGmoq56KxYr/HxhJOO9997DA8gPMr96AG3G7F2QFs6euFvfPR1SSEIIIQrB4CukDIM6p32BVwR93XXXXTDNIJvgXsLEI4Tk+YQeiB54VjixqZtvCCoDkAUw9GBjNptNnt7BBhqH55XL5Q8//NDa5iHsLDiZvGSIOZ8Aa52SBdsbjQZKwIA+vkKt2ID1EXx2DmH0HNsvX76Mhd1wLTykjp2hIw8ePAgdia8gNFHPF198EZfj4XaWZO71nenXF7ACcKXLSppxU52D7vCrcXCpv7LHkfO34vCLFy9yI/dYOOvUBFeuXImneCcG48X5SQMPUBz7yktZeuReHGQYhOPGU2sBRg4gjD7++GNPRRiUwEVZkhLNOKTfY+pipJCEEEIUgkHrYGMSnUMZA6yx3QTLBa6Ubdu28Tpy+OqXv/ylmf3t3/6tmb355ptYqxsGkSfCsU7Tb2FhgSPZWHZAo1SrVewJQxWuJtQHJSMJgnXmZcCxqACienymPUqDt4ajm1C3L774ApIO+ozFHLw4bibD1RQE3fHdw3UBHm1HfBHmJL355puwBBFTh+RscMX51Po//elPwf1J/JlyjpX36ZD6usA6yQVHRngq9kT74dhItEwX0Ow+YT0Bb+j8/HycTyFx0T8U6O3fIjnFs/riRSb5EIdPxJmZFhcX05w6HinH8/w8RNbaDw5GAs6dOxd7j/g2Li8vJ67W0TN92tQHv0MKFHqw0aInLTv3Rq1WQzI6nydr7f/l+++/38w++OADdF3oFeC05ybufQ+GtjjxNg/clUolNFDuivj/3beja8Q4Aw9xoD6ogHU+bCgB80y9U+GRQBSO0QY8q/hzKZVKnL+cQxtwOdu3b8cZefkcvMcQJQ5fXFzEx8cff9zajy4/sefPn8cQUFdjIueD16fP57rjf7scgIAf0cG9xWRPgB327NljbRto69at3IaxA5q6B0xzxA1eeVY48OAI/kE5pXej0fA2Zp1RD0GwDz8s3KWh8QfJvPmkQQp/fIWN6IrQCX388cfWmVgrKIcr4FFOQ95QNWQnhBCiEAy+QiolpQ6Kxw0Sv4p10rlz52Bt8WgV/PC/+93vzOwf/uEfIBFg4MNI5PynMCR93SCMPMC8QlgELErrjAhnmzRIeQJLEPKCR+EgWRqNBqsZxBRwnK6bpTw0wWLlrrvusrateu7cuRLl3/SQcT+wUqkgeptlHIbp2IQcGRn5+7//e2sLO0QucDnvvvtutjM5QCNy+UkcM0jbzSUCJ3GP54da+xeHIMDvjpYJpX733XezIIhjDbwy+Ij2xrEPgcjgY3nkwNotzZORW+ekXU9+H8RbW6d6c73CFxvnTrX2U4CnHtoIrxi6SLzb8WhNkMQ2ja5xWP2OFJIQQohCMAidaldio6yrDyl+D5aXlxEs4PNkrW0iIW3ikSNHvve971lnelbEhcPygi4JJsbyyhQOZ99BfHmQjRE15NADXAsCE9y6ZCsS8g6ihx25jUYD1wXFht0Qhw2gbEqlEgQQrzbrSTZxXfBR4zJxIkTMw4TEIU8++eRjjz1mZv/6r/9qnQYsQhs8WoF/o7QfJW2jSGSlapITBeE34l+/1WrFjwyLBm/naMPQxJxox39cbpBo1TzF1Z/fWCEBH0XggAg+RdCiWBvxPAe+EEuKd+doo4WFBUgxeGTxGvjYLJIvnGCp60zYxHIGNYpHCkkIIUQhGHyF1NVGSJwo13V/uD1g5kNPIKz56aef/s1vfmNmP/7xj61tKyFQG6IBI8sLCwscTcTj3T5mjf25Vp7g1WhiIHxFiJ+GcwjhPR5Fij0hm/Ce17GFfPGwN0TTcVw435Dt27ejBJY7HFB3+vRpGMIoB3eJA3CxeMfjjz/+29/+1jonIWKGr1uXeYy7/DH9fWowXj/yuxxY9LDfJbDWY78sew3dR4LRAm6N7q2JTxT7b/yksYDmJVmtM38r18fhqDbOaxXLqQBE7uFB/vLLL/EEwTfMGiu4hFjSBTmTVsRgeIxipJCEEEIUgsHsZjNIs5ETbe0MYKHAb4TMN7DInn32WYgArEbxwgsvWFuFwKTasWMHDofIQH4gPjUspnq9HqTksXZOIE7SOjIyghLYPkVck7t2OOMqgKMIgUyByckbeR4uLM0rV664wWvRQhhm9uCDD7pUsk5txF63P/3pT5yDFfiUKcstZXL6k1b644oYVhKAHajW2WZYfLiMRhtG80Cjgp5wvRWv24LWiHLQzAI3KrciqO0vvviChxO4PXt8ndFE3ThVTzw/yc+LRxiOIlzIV199hauIZSJLtMXFRY4e5PywvSmkmOwJlP2CFJIQQohC8I1r166tdx2uC9/85jf/8pe/ZO+zGr9CbIlgyT5rDzHjI0yqX/3qV9ZpCo2NjWEwnQ00vIf1dOedd/LEcgTRwQZE8Fui64trhTXFK5UKj9RzrB10UrCANFug8EgFI/tQRRgx50UEggwLbGA+/PDDZvaLX/zC2sppZmaGc8zAn8SsJogo/8966623/vnPf+66W18QN/jsGXWrAT+rry7PTYtve7AFdUBrQduDJwlT+rZv3+45HazT94N26PPteC4UwG541l5++WVOA4HRAk726oekTXTjqNRms4nzYmCA3UXBSuosEPlKMc5x+fJlFIuBBAjEtDyqvZH4QweFF7/BSyEJIYQoBP094NgDaTkaVupmiO0R5LV76KGHeN2EBx980Mz+8Ic/mNk///M/W1sinD9/nid2IJcBD16fPn2aLThO6sUJ4kZGRnjUHhYc3sP8tM6oPLYlsZuPdLOdy0Pz7C66fPkyDFhWV7gQmLGff/456oPTwW+Em4B5TvC6jY+P49iXX3458WZ2/Tkygui6zmYfQn/S2uokFs3lcpkdgfwTJK7LwJGZ7FX64osvuKHyXJ9gob94ahFqAhGDtu1nxEbMjYNYgVpyH1Ja2BsOWVhYQPOGKsLgR9CK4nlRfJ8RBLtlyxZWe9cjbV2/e4/AEA3ZfT1/RqVSCbHgPoRl7UGJZ555xsx++tOfmtlrr73GbRcPIXodf0J4KAxPBR4SXAKitL0j4b8A5CtCB+ATY9mv6/l7vNo33HAD/3HwdF1O4uLTEuNMmj4giYDvH/zgB9Z++FETdimPj4+j/0YP6rfOcv9AK/01E/cv/ghGftIaPIhngq8eD7YOkpZaZ//n7SqY3xqU40mn4h6CYxa2b9/OxeL/nbOpBhMq4pzfGPEul8txHiDshi7NnzgeWMsYIo67Sd6yYcMGDmT42vKo9l2D15CdEEKIQjAIKi8nq7RHEnV6XGar1fJFUa0tF2BzYaosps2OjY0hRw7sLwwFcK4gLxzmJ68ugR2QdPXUqVOII8dGjAcit6MLIM6hAiMU6ipYCpYFE97zKB8su4sXL6KSaclRHn/8cawoAenDmVR4aahXX33Vh1b41nFpKxqy4zpkxIIP7ZDd9RjPcVXB6jwxSY+ZLS8vp/0ErXZqongMjQO40WB82aTERKhmNjY2xvN2ET7AzZUXj/CjOKY8iNHImChiNGaQ+JVvuXTp0toudzSoSCEJIYQoBIVWSEtLS3/96195dbtarbZ582bfUqvVqtXq5s2bMwoJbO04FjZnjplgS4bdBCNoZmbG2hHPGLbGenff/e53zezYsWOIC4BbCHHPOBASxFeQ5BBSBBQAmHi7du1C4QgnhbeGncnlchnlwLMFnxYEGWfmL5fLPBCPWuFAHqZvNptxTheoNziN7rrrLugzBEHw3EC8RxRDo9FInL6adntzBnl3NTwLbpmuSYPPSaszfW1X0n4djw7g1RrRUHnagCuJtJ+g1Wqxl5F35synX375JTuruFZok5VKhWUTx+8AV0vxFeWJnA5qGCgtgOgJqC5UPrHB5+R6+P8KS6Ev74033mg0GocPHzazpaWlgwcPTkxM1Ov1iYmJw4cPHzhwoFqtnjlz5ic/+cnu3bvXu7JCrBY1eDHkFLdD2rdv35kzZ77//e/j4xtvvLFr1y48qw888MDk5GS1Wj169Gi9Xj9y5EjG85mhbFY5rTL7dHiDWDLkG4WUgYH21FNPQTHA1fSjH/3IzJ577jmjSGgO7+HIH6yBBtPMw9XY4YSd/ULYVMTy6pAsMGCBpzbhAXReYcw9BNBDP/vZz3wj5h4izP3FF1/0oXzrjAtnFRgYel3veZqVGgzx97WLaK0a/HUi513l2d9xvJwlCaDAL8uni5/WUtJ6m3yKL7/8kgtPi/rLT+wQir1cvhGvCKgLZs72zMCrIqa4l3r8+PHp6Wn/uHfvXrxZWloys1qthnk21WoVf7JC9DVq8EIUt0MKwDD63Nzc888//8QTT1y4cMEH1qE/Vg9H7/R8uCWNROMfBMlDsU7E5OQkpBKyr0JPHDp0yNrTld577z18hSX+YHmxZkL5kEf+ETOQWNMsLi5iWB8zUjlsCXF6voUNOrYKETWHau/atQtTlFB5rLiBqEJ3F3G0FdxmuITVa9OcnqTB4Ho3+Ay7O9GB15VYkXC0WxCNluizyXAixtvjWnE4aEAc9ZdILIZWWgLXYTUBdUPlN2L6Kcpuenp6dnb22LFj+/fvt3ZAp7X/7oUYMNTgxbDRN33v7Ozs/Pz8sWPH8HFychJxXHD5rskpVmlfZzirOEsC9MSOHTvuu+8+60zP8+yzzxppnd///vfW9hhBCSE6Dv9HOMQzNUCgcNJGrwZcOJ4xxdqThLAF35bLZc6U+uSTT1pbYyEJLOL3Xn31VcgdwOPmKK1cLiMPLGYgpcnNIPQxpypKPOS65ndYR1bf4BNjxvKwGpOc7y3PK3K/S+IUJf7Ir3GWhK4hcF1/3K4KKXjP4w2x98hPx3po9ZONVpRlo+cfuoD0zTXUarV6vX7gwAF8fOWVV2ZmZqanp0+ePImRrhUR/zH5X2TP/1k5Z2WCzz77DD0Nhl/gHsAreohTp0552jdrj5ghHTLGwTyxEEbk0GOhz0PUKfIDtVotlIAoc/SCiBG/6667vFbbt29HjAPG8RDOgODvF1980atdqVQ8YYzXAU8pkoufOnUqT9x8/tubEc6w0qJWuvP6svoGH89w+NqIu6Vms8kZffj/3Sdop4UPxIEJ1tmZZfRYvLPHhRsF+8RRNsHfAh8b18cn/OaceNADXf9YEqu95tX4eih0vRFiBI4ePRp8e/z48bm5uUcffZTnbQjRv6jBiyFniJKrFhlYNMgztHPnTjMbHR2F3OFBAH71MAdEXcPuQwnAV5TxnCvWDvWGhEIJ2LJhwwZex4iTWvLKNK6QUCavjRQsfLlWpuL1EKxO8XNN5ic7ueq6EPwELIDQqKDasQNCpS1J7gRSJm1Mj88SwNoI3fny8jLP+E6sMJ+O6xbEcxdEdifehKBuxW/w/RTUIIQQYoAp9JBdXxOP4HeNV4b/Bq/lchkx3Fi6gvOc+qC8mVUqFSgb4EtOeMmLi4ssdHhAPFj6iBcRwOmwBTvg9fPPP0fwetcog643J6ddufoZymK9yAjzQetFDE4Q7xA3Ud4SKCQuOTFoO5jWbSSq0nxOQcmxm6qYTWvdBfGaIIUkhBCiEAxCp9oz1zUIOGfMXmymudMI0055SiwkCxaPQOxcpVLheG4GZd52222J+Vd8C9xOHiwEVxNUGqfuD2zDtFH7nMGKxbQxxdqSEajNr4nOIQ6wTlRIGzZsCHZz8FWa9ElsmTw/IZhpEK/jF++2JgxGmNwqkUISQghRCIa0HwbrYqcnGlaJs3aCaYPWFiuYmgpfTlBCvLBeq72EeXwitlI9Ro6N1ngcP6hnvDGQUGtlReZ0yKXVSqwLXc38xPlDcbMPSgsEU3CiDRs2xBNs+XQo+fLly7xAZVcp38MM3JWyIlXUdXpZn7Z/KSQhhBCFYKgVEpM/7it/gba6QedsI8i/jYOXEvNLdq1tXLG00fOcBa60DhmsKMtDn9qGQ06GFglUTizcg91iX1Si4yoOn0s8aTGb06A6maSQhBBCFILB7GZ7oDcF0DV8rufSSkkLkcUD4olOncQtPDSf6PtJu6LEKCl+n9Mr1pvvJ/uoNde1Q0jXX7znndeKxKaV83dnhZRYGsjwhvYpfepJVYfUIxljcWlNYUW9V+LOXZ+ljHK6RgTkDNHu2immcT2G2vrueSssOfuV+Nf/+kmbKRGQFhwRPKdFo2u0wmCjITshhBCFYEj74bWl62jVWpFzcCyD1Y99pUnDrumReqhtD/uLlVJkY7yroF9N3HNaCUUY7BpanSSFJIQQohAMXQ/M5Mx5Whx6q0+ea+lacqKfLH7fW+H56yl6YzUWd7aSsOv5w13XEP804X5dW+DX08g9RqnvkEISQghRCPqyF10r1ipou+CsybV8DX6yQbrnxaQ3nbRWCriYfM2VX9FowRAihSTE8NLb33HPf+I63bqfruAMaT8sxFDRg8Vd8ECv1VevmH/oa3LDC/urdUUKSQghRCHo145UCHFd6V8rW/Qvg9zmbr311vWughBfH2rwot/5xrVr19a7DkIIIYR8SEIIIYqBOiQhhBCFQB2SEEKIQqAOSQghRCFQh2RmtrS0VK/X/WOtVltaWup6VK1W46PyOMZRrAAABntJREFUHzg3N8e75TzKzOr1em8H9lDJer3e29X1drp4t94O7K2Sw4YafOK51ODXnf/33/7t39a7DuvPq6+++sknn9xzzz1mduDAgXq9PjMzU61Wq9Vq4v5LS0vf/e53//d///ftt98+d+7cSg+8du3aL3/5y7/5m7+55ZZb8hzlx/7d3/3dD3/4Q3zMeeDOnTvn5ua8njmPOnLkyH//93+//fbb//d//zc5OZnzdLOzs88///zbb7/99ttv//KXv9y9e/dPf/rT/Dfz+eefL5fL+U+3tLT0T//0T0tLS//xH/+xadOmFd3MIUcNPkANvihcG3oefvjhW2+99bnnnrt27dof/vCHn//859euXbtw4cI//uM/ph3y0ksvYf9r165961vfWtGBL730ku+W8yjw85///Fvf+laj0ch/uuDbnEd99NFH2K3RaPzoRz/Kf6Dz2WefPffcczmPev3113EzG40Gdst54HPPPff666/jwIcffnillRxa1OAD1OCLwyBPjM3J8ePHp6en8b5Wq8FgqVarJ0+eTDtk7969eONiOeeBhw4d8v2r1WrOo8xsZmbmpptuqtfrmzdvzn867H/kyJFNmzY98cQTOY/66KOPbrrpptnZWTM7duxY/tM5zz///LFjx1566aU8R01MTMzOzs7NzZ05c2ZiYmJFp9u0aZOZbd68+cyZMyut5NCiBh+gBl8c5EMKceW7a9eujH2q1erc3NzBgwefeOKJ/AeC6enp559/3htT16NqtVqtVvNnO//plpaWxsfHp6amKpXKU089lb+Sr7/+uplduHDhwIED+U8HZmdn77nnHvyP5LyZmzZtOnHixIkTJ2666ab8p7vnnntmZmZOnDjx1FNP4UHNX0nhqMGbGnxhkEIKccdmttExPT09Pz9/7NgxbxM5DzSzw4cPP/HEE9/+9re/853v5DlqZmZmfHx8enq6Xq8fOXIE/wh5DpyampqamjKz3bt3nzhxIn8l9+7du3//fjPbt2/fSq9udnb25Zdfzn/USy+9NDU1hdM98MADeJPnwN27dz/zzDMnTpzYv3//3NzciiopHDV4U4MvDFJIHUxOTl64cMHM6vU61HQis7OzwcOZ88AjR47gOYE9lfOoQ4cO3XPPPbDCpqamNm3alPPAmZkZNNwVVdKtNmsP0eQ80MxqtdrExMSKri4m54EnTpw4c+bM4cOHq9Xq7t27ez7dMKMGb2rwRUK57MzMMKR++PBhM9u3b9+uXbtOnjx56NAhWFsxR44cOXPmDFSzmb3yyis5D6zVak899dTevXvPnDlzzz33HDp0KM9RzoEDB3CuHk4Hoyzn6bCbH5XzdGY2PT1900034ZCcR9Xr9QMHDqCSExMT+X+FpaWlffv27d279+TJk88888zk5OSKbuYwowYfoAZfFNY7qqKIfPTRRxcuXLhOBzYajWC363q6eLfejuq5ntf1dLiZCMRaTSWHHDX4xN3U4L9+pJCEEEIUAvmQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQqEMSQghRCNQhCSGEKATqkIQQQhQCdUhCCCEKgTokIYQQhUAdkhBCiEKgDkkIIUQhUIckhBCiEKhDEkIIUQjUIQkhhCgE6pCEEEIUAnVIQgghCoE6JCGEEIVAHZIQQohCoA5JCCFEIVCHJIQQohCoQxJCCFEI1CEJIYQoBOqQhBBCFAJ1SEIIIQqBOiQhhBCFQB2SEEKIQqAOSQghRCFQhySEEKIQ/H9nG4uVfZh8HgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%-Display\n", "A = spm_slice_vol(ref,spm_matrix([0 0 70]),ref.dim(1:2),1);\n", "F = spm_slice_vol(src,x\\src.mat\\ref.mat*spm_matrix([0 0 70]),ref.dim(1:2),1);\n", "subplot(1,2,1), imagesc(rot90(A)), axis image\n", "subplot(1,2,2), imagesc(rot90(F)), axis image\n", "colormap gray" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "%-Save transformation in NIfTI header and reslice\n", "% spm_get_space(src.fname, x \\ spm_get_space(src.fname));\n", "% spm_reslice([ref,src], struct('mean',false,'which',1));" ] } ], "metadata": { "kernelspec": { "display_name": "Octave", "language": "octave", "name": "octave" }, "language_info": { "file_extension": ".m", "help_links": [ { "text": "GNU Octave", "url": "https://www.gnu.org/software/octave/support.html" }, { "text": "Octave Kernel", "url": "https://github.com/Calysto/octave_kernel" }, { "text": "MetaKernel Magics", "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" } ], "mimetype": "text/x-octave", "name": "octave", "version": "4.2.2" } }, "nbformat": 4, "nbformat_minor": 2 }