mwcraig_test.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "from astropy.wcs import WCS\n",
    "f = fits.open('ey-uma-042R.fit')\n",
    "data = f[0].data\n",
    "hdr = f[0].header\n",
    "w = WCS(hdr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x10bde1550>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgwAAAILCAYAAACEg8S3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X9sVHW+//HX/CylpZ22WzFSdO0SKfS6f1z81Y3SgLDq\n6v4DK3ETk01gi79yQ9wssuo/sslu1GyW1Yhkl7gbf2QRd73ExOSCsEGWDbgoau4aWrOQFS5chdKh\nZabTzq8z94/9tl+7LXymvGd6pvJ8JCZyZl79fObMmemrZ+acE9i7d29BAAAAFxH0ewIAAKDyURgA\nAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2egCT19/frgw8+0JVXXqloNOr3dAAAuCxk\nMhl98cUXuuGGGxSLxS5634ooDB988IF+9rOf+T0NAAAuS08++aSWLVt20ftURGG48sor/Z6C7wKB\nQFl/fqFQKOsY1p9d7sdvHT+fzysUCl3w9mCQT/cuJpvNKhKJlHUMz/PK+vPLrVC48El3XdufKz8V\nrONb3wOm+/Pvt2J+D1dEYeBjCApDpRcG130oDBcXCAQqfh35/Qv3YgKBgO+vkXIr53tIud//SsHv\n7a+Y38OV/QoGAAAVgcIAAACcKAwAAMCJwnCZqPTP7ypdpX/+XulcX9jDxfH6RSXgXfAywRuODYXB\nhsJgw/Znw/tfabAVAgAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAp4q4+FQpcJyt\njd8Xj/L7OHPreQL83v78Ht/vC+dI9nXg99UW8/m8KW+dv9+P3+/xrfye/1S8BtnDAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAicIAAACcwn5PAP/k97XcreMHg7buGQqFfB3f+viteevj93v7sSoUCuafkc/nfZ2D\n53mmvHUb8Hv+pXgO/WSd/3R9DU7mcbOHAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOF\nAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2eQKWwXsv8\ncs+Hw7ZNKRi0dddQKGTKW+dvHd+at64/az6Xy5ny+XzelJckz/N8zVsfg9/zt74HFAoFU946f7/H\nt64/P8cvdt2xhwEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAA\nOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOIX9nsC/sl5T3C9+zzsUCpnywaCtO1rH\nD4dtm2IkEjHlq6qqfB3f+viteatcLmfKZ7NZ8xzy+byvc7COb12HnueZ8tb5W/PW91C/x7fyc/xC\noVDU/djDAAAAnIr+s6Snp0e7du3SRx99pNOnT6u+vl4LFizQmjVr1NLSMnq/p59+Wu+88864/Ny5\nc/Xyyy+XZtYAAGBKFV0Ytm3bpiNHjqizs1Otra2Kx+PasWOH1q5dq82bN+vaa68dvW8kEtH69evH\n5Gtqako3awAAMKWKLgyrVq1SW1vbmM+qlyxZotWrV2vbtm164okn/v8PDYe1bNmy0s4UAAD4pujv\nMLS3t4/7YtucOXN0zTXX6MSJE2OWFwoFeZ6nwcHB0swSAAD4yvTV6kKhoHPnzqm1tXXM8nQ6rbvv\nvlvpdFqzZs3S0qVLtXbtWlVXV5smCwAA/GEqDHv27FFfX5/WrFkzuqypqUn33XefrrvuOnmep0OH\nDumtt97SsWPHtGnTJvPhdwAAYOpdcmE4ceKEnnvuObW3t+uOO+4YXd7V1TXmfkuWLFFLS4teeukl\n7du3T0uXLr302QIAAF9c0nkY4vG4Hn/8cdXW1mrjxo3OE07ce++9CgQC+vDDD50/u1AojPsPAADY\neJ6nfD4/5r/JnPBr0nsYksmkNmzYoMHBQT3//PNqbGx0ZqLRqOrq6pRIJJz39ftsWwAAfBVNdEbf\nkYMUijGpwpDJZPTkk0/q1KlT+sUvfqGrr766qFwqldLAwIBisdhkhgMAABWi6I8k8vm8Nm7cqO7u\nbj311FNauHDhuPtkMhmlUqlxy1955RVJ0o033miYKgAA8EvRexi2bNmigwcPqqOjQwMDA9q9e/eY\n25cvX654PK6uri7dfvvtmjt3riTp/fff16FDh3TTTTfp1ltvLe3sAQDAlCi6MBw7dkyBQEAHDx7U\nwYMHx9wWCAS0fPly1dbWqqOjQ4cPH9auXbvkeZ5aWlrU1dWlVatWlXzyAABgahRdGDZt2uS8T21t\n7ZhTRAMAgK8G04mbvkqsR2dM9/xE356djHDYtilVVVWZ8jNmzPA1b7242syZM0156wnRrM//RN9d\nmox0Om3Kl+Jn5HI5X8fPZrOmvHX+1vGtJnN4XzlYXwP5fL5EM7k0U3EKAtsaAgAAlwUKAwAAcKIw\nAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAA\nACcKAwAAcKIwAAAAp7DfEygV67XMA4FAiWZyaazzD4VCpnwkEjHlq6qqTPkZM2aY8rNmzTLl6+rq\nTPmGhgZTvqamxpS3bj/WfC6XM+WHhoZMeUlKpVKmfCKRMOWHh4dN+WQyacpnMhlT3voeaM3n83lT\n3roNep5nyvv9O2QqsIcBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADiF/Z7AlwUCgWl7TfFQKGTKB4O2\n7hYO257KSCRiys+YMcOUj8VipnxjY6Mp39TUZMo3NDSY8rW1taa8dfux5rPZrCmfTqdNeUk6f/68\nKV9XV2fKx+NxU976HCSTSVPeb4VCwZS3vgfmcjlT3iqfz/s6fjHYwwAAAJwoDAAAwInCAAAAnCgM\nAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAA\nwInCAAAAnGwXEK8ggUBgWueDQVt3s+YjkYgpP3PmTFM+FouZ8rNnzzblv/a1r5nyV1xxhSlvffzR\naNSUt24/nueZ8gMDA6a8JPX395vyvb29prx1HVoVCgVfx7duA9b553I5U946f7/l8/lLyk1mvbOH\nAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEA\nADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2ewJcFAgEFAgFfxg6FQqa8dd7hsO2piEQipnw0GjXl\nZ86cacrHYjFTvqmpyZSfM2eOKT937lxTvrm52ZSfMWOGKW/dfguFgimfSCRMeUk6deqUKW99DVnX\nQS6XM+Wz2awpn8/nfc1b15/neaa89T3Yuv6nA/YwAAAAp6IrVU9Pj3bt2qWPPvpIp0+fVn19vRYs\nWKA1a9aopaVlzH2PHz+uzZs365NPPlEkEtHNN9+sRx55RPX19SV/AAAAoPyKLgzbtm3TkSNH1NnZ\nqdbWVsXjce3YsUNr167V5s2bde2110qSent7tW7dOs2aNUtdXV1KpVJ644039I9//ENbtmwx7/YB\nAABTr+jf3qtWrVJbW9uYz/qXLFmi1atXa9u2bXriiSckSa+99prS6bS2bt06+rlsW1ub1q9fr507\nd+qee+4p8UMAAADlVvR3GNrb28d9MXDOnDm65pprdOLEidFl+/fvV0dHx5gvcS1atEgtLS169913\n7TMGAABTzvSlx0KhoHPnzo1+N6G3t1f9/f2aP3/+uPu2tbXp73//u2U4AADgE1Nh2LNnj/r6+rRk\nyRJJUjwelyQ1NjaOu29jY6MSiYT50CEAADD1LrkwnDhxQs8995za29t1xx13SJLS6bSkiY/pH1k2\nch8AADB9XNIhC/F4XI8//rhqa2u1cePG0ZO+VFVVSZIymcy4zMiykftMZKITbwQCAQWDnC4CAACL\nfD5vOsHVpAtDMpnUhg0bNDg4qOeff37Mxw8j/z/y0cSXxeNx1dXVXfSwymAw6NuZHgEA+CoLhULj\nDl7wPK/orwpMqjBkMhk9+eSTOnXqlH7xi1/o6quvHnN7c3OzYrGYenp6xmV7eno0b968yQwHAAAq\nRNH7+vP5vDZu3Kju7m499dRTWrhw4YT3u+222/Tee++pt7d3dNnhw4d18uRJdXZ22mcMAACmXNF7\nGLZs2aKDBw+qo6NDAwMD2r1795jbly9fLkm6//77tW/fPj366KNauXKlUqmUtm/frtbWVt11112l\nnT0AAJgSRReGY8eOKRAI6ODBgzp48OCY2wKBwGhhaG5u1q9+9Su9+OKL+s1vfqNoNKpbbrlFDz/8\nMKeFBgBgmir6N/imTZuK/qFf//rX9eyzz17ShAAAQOXhT/4SsR7dYc1bDz292OGuxaitrfU1H4vF\nTPkvn8r8UsydO9eUb2hoMOVnzpxpyv/rN6cnK5/Pm/KpVMqUl+yPYaLDwSdjaGjIlM9ms6a89aR4\n1ufQ7/EthwuWIu/3e/ilCgQCRT93nOAAAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UB\nAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADiF/Z5ApbBeC916\nLXO/r6UeDts2haqqKlM+EomY8jNmzDDlY7GYKV9TU2PK19fX+zq+df0XCgVTfmhoyJSXpHw+7+sc\nrPl0Om3KZ7NZX8cfHh72dXzre6DfeevvoKnAHgYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEA\nAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4BT2ewKV\nwu9roVtZxw+H/d0UQqGQKV8oFEx5v69Fb338kUjElK+qqjLlrayPX5IaGxtN+cHBQVP+7Nmzpnwi\nkTDlrfOvrq425a3bUDQaNeWz2awpP91/B1zqe+BkcuxhAAAAThQGAADgRGEAAABOFAYAAOBEYQAA\nAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABO\nYb8nUCqBQOCyHt/zPFM+l8v5On4mkzHl0+m0KZ9MJn0d3/r4J3NN+3IIh21vJaWYf01NjSlfX19v\nyl9xxRWmfH9/vylvffwzZsww5aurq035wcFBUz4YtP396/d7uN/jF4M9DAAAwInCAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAyXYRe4zyPM+ULxQKvo6fy+VM+VQqZcqn02lTPplMmvLnz5835QcGBkz5WCxmylvX\nX1VVlSlv3X6DQfvfLuGw7e3Mug6sz2FdXZ0pH41GTflIJGLKW9e/NR8IBEz56e5SH/9kcuxhAAAA\nTpOqdENDQ3r99dfV3d2tnp4eJZNJPfbYY7rzzjvH3O/pp5/WO++8My4/d+5cvfzyy7YZAwCAKTep\nwjAwMKBXX31Vs2fP1rx58/Txxx9fcHdGJBLR+vXrxyyrqam59JkCAADfTKowNDU16c0331RDQ4M+\n/fRTPfTQQxf+weGwli1bZp4gAADw36S+wxCJRNTQ0FDUfQuFgjzP0+Dg4CVNDAAAVI6yHSWRTqd1\n9913K51Oa9asWVq6dKnWrl2r6urqcg0JAADKpCyFoampSffdd5+uu+46eZ6nQ4cO6a233tKxY8e0\nadMmhUKhcgwLAADKpCyFoaura8y/lyxZopaWFr300kvat2+fli5dWo5hAQBAmUzZeRjuvfdeBQIB\nffjhhxe8j+d5yufzY/6znpAIAAD88wR96XR6zH/ZbLbo/JSd6TEajaqurk6JROKC9wkGg5f92boA\nACiHcDg87oyanucVfabYKdvDkEqlNDAwYD59KgAAmHolLwyZTGbC6wq88sorkqQbb7yx1EMCAIAy\nm/RHEjt27FAymdTZs2clSQcOHNCZM2ckSStWrFAikVBXV5duv/12zZ07V5L0/vvv69ChQ7rpppt0\n6623lnD6AABgKky6MLzxxhs6ffq0pH9e5eovf/mL9u/fr0AgoG9/+9uqra1VR0eHDh8+rF27dsnz\nPLW0tKirq0urVq0q+QMAAADlN+nCsG3bNud9nnjiiUuaDAAAqExc3hoAADhN2WGV5VYoFEx56/ke\ngkFb9/J7/plMxpQv9rCcCxkaGjLlrdcs6e/vN+XPnTtnyluPHrJeCdZ6ynbr2Vut279kfw1GIhFT\nPhqNmvLWdWidv3X86X5IvHX+1m3Yr/xkcuxhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOYb8n8GWW\n64FP92ux5/N5X/PZbNaUHx4eNuXPnz9vytfU1JjyiUTClI/H46b8FVdcYcqnUilTvra21pSPRqOm\nfClev9afEYlETHnL+1cpxg+FQqa8df0Fg7a/P/1+D7c+f5cD9jAAAAAnCgMAAHCiMAAAACcKAwAA\ncKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCi\nMAAAAKew3xP4qrBeS92a9zzPlM/n86Z8JpMx5YeHh30dP5lMmvLnz5835QcGBkz5pqYmU966/qzb\nTyQSMeUlKRQKmfLBoO3vJ+v4gUDAlLfye3y/3wPhxh4GAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOAU\n9nsCXxYIBHy7Jrv1Wux+563Xgs/n86Z8JpMx5YeHh0358+fPm/IzZ8405QcHB0156/yHhoZMeev2\nY33dhsP2tyLrY4hGo77mQ6GQKW99DVufw2w2a8pb529lHX+6/g6YTI49DAAAwInCAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAyX4R+hKyXg/cwnot+Eu9FvmIUChkylvXnfVa8Nls1pRPpVKmfCKRMOVnzZplyieT\nSV/zmUzGlLc+f9bt37r9VQLrOrAKBm1//1m3oVwu52ve+h5off78fv6nwqQKw9DQkF5//XV1d3er\np6dHyWRSjz32mO68885x9z1+/Lg2b96sTz75RJFIRDfffLMeeeQR1dfXl2zyAABgakyqkg4MDOjV\nV1/V//zP/2jevHmSJv7LvLe3V+vWrdPnn3+urq4urVq1Sn/961/14x//2NwiAQDA1JvUHoampia9\n+eabamho0KeffqqHHnpowvu99tprSqfT2rp1q5qbmyVJbW1tWr9+vXbu3Kl77rnHPnMAADBlJrWH\nIRKJqKGhwXm//fv3q6OjY7QsSNKiRYvU0tKid999d9KTBAAA/ir5URK9vb3q7+/X/Pnzx93W1tam\nv//976UeEgAAlFnJC0M8HpckNTY2jrutsbFRiUSC7zEAADDNlLwwpNNpSVI0Gh1328iykfsAAIDp\noeTnYaiqqpI08TG9I8tG7vOvCoXCuGNpA4GA+RwJAABc7vL5vOl8ESUvDCMfRYx8NPFl8XhcdXV1\nCocnHpZyAABAeYRCoXEnCfQ8r+ivCZT8I4nm5mbFYjH19PSMu62np2f0/A0AAGD6KMu1JG677Ta9\n99576u3tHV12+PBhnTx5Up2dneUYEgAAlNGkP5LYsWOHksmkzp49K0k6cOCAzpw5I0lasWKFampq\ndP/992vfvn169NFHtXLlSqVSKW3fvl2tra266667SvsIAABA2U26MLzxxhs6ffq0pH9+5+Avf/mL\n9u/fr0AgoG9/+9uqqalRc3OzfvWrX+nFF1/Ub37zG0WjUd1yyy16+OGHL/j9BQAAULkm/dt727Zt\nRd3v61//up599tlJTwgAAFSesnyHAQAAfLXw+cD/Y72WupX1Wur5fN6UDwZt3dE6fjabNeUHBwdN\n+fPnz/uaP3funCnf399vyg8NDZnyw8PDpnwpDqf2exuynpAumUz6mrfO33oG34nO3TMZfr8HWll/\nB03F7zD2MAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoD\nAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAp7DfE8A/Wa9l7nmer+PncjlTPpPJmPLpdNqU\nTyaTpvzAwIAp39fXZ8qfPHnSlP/a175mygeDtr89qqurTXlJyufzpnx/f78pf+bMGVO+t7fXlE8k\nEqb80NCQKZ9KpUx56/NnzWezWV/Ht+anAnsYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAA\nOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIBT2O8JVIpC\noeDr+J7nmfLBoK375XI5U94qFAqZ8plMxpRPJpOm/ODgoCkfj8dN+f/93/815Wtra0156/qzji9J\n+XzelD937pwpf+LECVP+7NmzpnwikTDlrY8/nU6b8tls1te89T3Y798hlzr+ZHLsYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBE\nYQAAAE4UBgAA4ERhAAAATmG/J/Blfl5PPBAI+Da2ZH/suVzOlA+HbZtCPp835a3Xsg8Gbd3Xmj97\n9qwpb13/1rzneab8mTNnTPmamhpTXrI/hoGBAVPeug34nT9//rwpPzg4aMpb3wOsz781b30P9ys/\nmRx7GAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4\nURgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAU9jvCZSK9Vrilzvr+svn86Z8MGjrrplMxpQPBAKm\nfH9/vylv5XmeKT88PGzKnzlzxpSvrq425SX7Njg4OGjKDw0NmfLxeNyU7+vrM+WTyaQpb30NZrNZ\nUz6Xy5ny1teQNT8dfoeVpTB8/PHH+tGPfjThbZs3b9aCBQvKMSwAACiTsu5hWLlypebPnz9m2VVX\nXVXOIQEAQBmUtTBcf/31Wrx4cTmHAAAAU6CsX3osFApKpVLmzxYBAIC/yrqH4dlnn9XQ0JCCwaC+\n+c1v6oEHHhj3EQUAAKh8ZSkMkUhEixcv1i233KL6+np99tln2r59u9atW6cXXnhB8+bNK8ewAACg\nTMpSGNrb29Xe3j76746ODi1evFg//OEPtXXrVj3zzDPlGBYAAJTJlJ2HYc6cOfrWt76l/fv3q1Ao\nXPC494mORbUeIw8AwOXO8zzT+R6m9MRNzc3NyuVyGh4evuCJWigHAACU3kQnyCsUCkWfdGpKTw39\n+eefq6qqqiRndQMAAFOnLIVhotPkHj16VAcOHNANN9xQjiEBAEAZleUjiZ/+9KeqqqpSe3u7YrGY\njh8/rrffflvV1dXq6uoqx5AAAKCMylIYbr31Vu3Zs0d/+MMflEqlFIvFtHjxYv3gBz/g1NAAAExD\nZSkMK1as0IoVK8rxowEAgA+m9EuPAABgeprSwyqLcanHiFoPx5wO1yIvJ+u15EOhkK/jW/m9/Vjz\n1vWXTCZN+aqqKlO+Eg6nHh4eNuWHhoZM+VQq5WveOv9MJuNrvthDAy/Ees0jv98DLjU/mRx7GAAA\ngBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAU9jvCZSK9VrigUDA1/GnO+u15K2s69/zPF/z1vlb138qlTLl\nI5GIKR8M2v92sa6DXC5nymez2Wmdz2Qyprx1/VlfQ9bn3/oa9Ds/FdjDAAAAnCgMAADAicIAAACc\nKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgM\nAADAicIAAACcwn5PoFQCgYApb70WuXV8q+lwLfWLsV7L3m/W9W/NZ7NZUz4SiZjy1u3f79ePJHme\n52s+l8uZ8tZtwLoNWl/D1rzfr0Hr8+/Xe/hkxmUPAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCi\nMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwCvs9\ngVLx61riIwKBgCnv9/z9Nt0fv3X+1nwwaOv+uVzO1/ErgfU5yOfzvo7ved5lnff7Nei3qZj/9H+V\nAwCAsqMwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMA\nAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwCvs9ga8K67XcA4FAiWbiD+u12K2P37r+rfMPBm3d2zq+\n39tfJWy/1jlYnwO/8/l83tfxL/e8ld/jF6MshSGTyeh3v/uddu/erWQyqdbWVq1Zs0aLFi0qx3AA\nAKDMyvKRxDPPPKM//vGPWr58uf7jP/5DwWBQP/nJT/S3v/2tHMMBAIAyK3lh6O7u1t69e9XV1aUH\nHnhAd999t375y19q9uzZ+vWvf13q4QAAwBQoeWHYt2+fQqGQvvvd744ui0aj+s53vqMjR46ot7e3\n1EMCAID62ujdAAAJcElEQVQyK3lhOHr0qFpaWlRdXT1m+fz580dvBwAA00vJC0NfX5+amprGLR9Z\n1tfXV+ohAQBAmZW8MGQyGUUikXHLo9GoJCmdTpd6SAAAUGYlLwzRaFTZbHbc8kwmI0mqqqoq9ZAA\nAKDMSl4YmpqadPbs2XHLRz6KmOjjCgAAUNlKXhjmzZunkydPKpVKjVne3d09ejum3nQ4i1gls55J\n8XKXy+X8nsK0xvqz4f2vNEpeGDo7O+V5nt5+++3RZZlMRjt37tTChQvV3Nxc6iGBsuMNx4ZfeDbW\n0z5f7nj9lkbJTw29YMECdXZ2auvWrTp37pyuuuoq7dq1S2fOnNGGDRtKPRwAAJgCZbmWxOOPP67f\n/va32r17txKJhL7xjW/o5z//ua6//vpyDAcAAMqsLIUhGo3qwQcf1IMPPliOHw8AAKZYRVzeeuSQ\ny8vZVHzGxud4l65QKFx0/fn9pchKvzx1oVAo++fwfl/e2so1vmsb83sb9Pvy0q689favumJ+D1dE\nYfjiiy/8nsJlz+8Xi9/jF4MvntkMDw/7PYVpjZPe2UyH9xg/ffHFF/q3f/u3i94nsHfvXt/XYn9/\nvz744ANdeeWVo2eEBAAA5ZXJZPTFF1/ohhtuUCwWu+h9K6IwAACAylby8zAAAICvHgoDAABwojAA\nAAAnCgMAAHCiMAAAAKeKOA/DhWQyGf3ud7/T7t27lUwm1draqjVr1mjRokV+T63iffzxx/rRj340\n4W2bN2/WggULpnhGlWtoaEivv/66uru71dPTo2Qyqccee0x33nnnuPseP35cmzdv1ieffKJIJKKb\nb75ZjzzyiOrr632YeWUodv09/fTTeuedd8bl586dq5dffnmqpltxenp6tGvXLn300Uc6ffq06uvr\ntWDBAq1Zs0YtLS1j7sv2N16x64/tz66iC8MzzzyjP//5z/re976nlpYW/dd//Zd+8pOf6Je//CXX\npSjSypUrNX/+/DHLrrrqKp9mU5kGBgb06quvavbs2Zo3b54+/vjjCc8a2Nvbq3Xr1mnWrFnq6upS\nKpXSG2+8oX/84x/asmWLwuGKfjmVTbHrT5IikYjWr18/ZllNTc1UTLNibdu2TUeOHFFnZ6daW1sV\nj8e1Y8cOrV27Vps3b9a1114rie3vQopdfxLbn1XFbmHd3d3au3evHnzwQa1atUqStHz5cq1evVq/\n/vWv9cILL/g8w+nh+uuv1+LFi/2eRkVramrSm2++qYaGBn366ad66KGHJrzfa6+9pnQ6ra1bt45e\npr2trU3r16/Xzp07dc8990zltCtGsetPksLhsJYtWzaFs6t8q1atUltbm0Kh0OiyJUuWaPXq1dq2\nbZueeOIJSWx/F1Ls+pPY/qwq9jsM+/btUygU0ne/+93RZdFoVN/5znd05MgR9fb2+ji76aNQKCiV\nSnFa44uIRCJqaGhw3m///v3q6OgYfbOWpEWLFqmlpUXvvvtuGWdY2Ypdf9I/t0fP8zQ4OFjmWU0f\n7e3tY37ZSdKcOXN0zTXX6MSJE6PL2P4mVuz6k9j+rCp2D8PRo0fV0tKi6urqMctHdq8fPXp0zAsH\nE3v22Wc1NDSkYDCob37zm3rggQfGfUQBt97eXvX390+47tra2vTXv/7Vh1lNP+l0WnfffbfS6bRm\nzZqlpUuXau3ateNe55e7QqGgc+fOqbW1VRLb32T96/obwfZnU7GFoa+vT01NTeOWjyzr6+ub6ilN\nK5FIRIsXL9Ytt9yi+vp6ffbZZ9q+fbvWrVunF154QfPmzfN7itNKPB6XJDU2No67rbGxUYlEQrlc\n7rL9HLkYTU1Nuu+++3TdddfJ8zwdOnRIb731lo4dO6ZNmzaN+yvxcrZnzx719fVpzZo1ktj+Jutf\n15/E9lcKFbt1ZTIZRSKRcctHLk7Fldsurr29Xe3t7aP/7ujo0OLFi/XDH/5QW7du1TPPPOPj7Kaf\nke1tooujfXmb5A37wrq6usb8e8mSJWppadFLL72kffv2aenSpT7NrLKcOHFCzz33nNrb23XHHXdI\nYvubjInWn8T2VwoV+x2GaDSqbDY7bvnINburqqqmekrT3pw5c/Stb31LH330EZd6naSR7W2ia8az\nTV66e++9V4FAQB9++KHfU6kI8Xhcjz/+uGpra7Vx48bRo03Y/opzofV3IWx/k1OxdbSpqUlnz54d\nt3zko4iJPq6AW3Nzs3K5nIaHh/ncbhJGdgWP7Br+sng8rrq6Ov66uwTRaFR1dXVKJBJ+T8V3yWRS\nGzZs0ODgoJ5//vkxHz+w/bldbP1dCNvf5FTsHoZ58+bp5MmTSqVSY5Z3d3eP3o7J+/zzz1VVVUVZ\nmKTm5mbFYjH19PSMu62np4ft8RKlUikNDAwoFov5PRVfZTIZPfnkkzp16pR+/vOf6+qrrx5zO9vf\nxbnW34Ww/U1OxRaGzs5OeZ6nt99+e3RZJpPRzp07tXDhQo6QcOjv7x+37OjRozpw4IBuuOEGH2Y0\n/d1222167733xhzSe/jwYZ08eVKdnZ0+zqzyZTKZceVfkl555RVJ0o033jjVU6oY+XxeGzduVHd3\nt5566iktXLhwwvux/U2smPXH9lcaFbsPa8GCBers7NTWrVt17tw5XXXVVdq1a5fOnDmjDRs2+D29\nivfTn/5UVVVVam9vVywW0/Hjx/X222+rurp63Jd/IO3YsUPJZHL0Y7ADBw7ozJkzkqQVK1aopqZG\n999/v/bt26dHH31UK1euVCqV0vbt29Xa2qq77rrLz+n7zrX+EomEurq6dPvtt2vu3LmSpPfff1+H\nDh3STTfdpFtvvdW3uftty5YtOnjwoDo6OjQwMKDdu3ePuX358uWSxPZ3AcWsv3g8zvZXAoG9e/dW\n7LffMpmMfvvb32rPnj1KJBL6xje+odWrV/MXchH+8z//U3v27NGpU6eUSqUUi8X07//+7/rBD37A\nqaEn8P3vf1+nT5+WpNEvShUKBQUCAf3+97/X7NmzJUmfffaZXnzxRf3tb39TNBrVzTffrIcffviy\n36XpWn81NTV6/vnn1d3drbNnz8rzPLW0tGjZsmVatWrVZX1I26OPPqr//u//nvCLyIFAQH/6059G\n/832N14x6y+ZTLL9lUBFFwYAAFAZKvY7DAAAoHJQGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIDT/wGKnqNabRh3HQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10b868bd0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "x, y = 325, 720\n",
    "from astropy.nddata.utils import Cutout2D\n",
    "cutout = Cutout2D(data, (x, y), (30, 30))\n",
    "plt.imshow(cutout.data, origin='lower', cmap='Greys_r')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<SkyCoord (FK5: equinox=2000.0): (ra, dec) in deg\n",
       "    (135.40390644, 49.73659214)>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from astropy.wcs.utils import pixel_to_skycoord\n",
    "coord = pixel_to_skycoord(x, y, w, origin=0)\n",
    "coord"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x10c1e7210>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAgwAAAILCAYAAACEg8S3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3X9sVHW+//HX/CylpZ22WzFSdO0SKfS6f1z81Y3SgLDq\n6v4DK3ETk01gi79yQ9wssuo/sslu1GyW1Yhkl7gbf2QRd73ExOSCsEGWDbgoau4aWrOQFS5chdKh\nZabTzq8z94/9tl+7LXymvGd6pvJ8JCZyZl79fObMmemrZ+acE9i7d29BAAAAFxH0ewIAAKDyURgA\nAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2egCT19/frgw8+0JVXXqloNOr3dAAAuCxk\nMhl98cUXuuGGGxSLxS5634ooDB988IF+9rOf+T0NAAAuS08++aSWLVt20ftURGG48sor/Z6C7wKB\nQFl/fqFQKOsY1p9d7sdvHT+fzysUCl3w9mCQT/cuJpvNKhKJlHUMz/PK+vPLrVC48El3XdufKz8V\nrONb3wOm+/Pvt2J+D1dEYeBjCApDpRcG130oDBcXCAQqfh35/Qv3YgKBgO+vkXIr53tIud//SsHv\n7a+Y38OV/QoGAAAVgcIAAACcKAwAAMCJwnCZqPTP7ypdpX/+XulcX9jDxfH6RSXgXfAywRuODYXB\nhsJgw/Znw/tfabAVAgAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAp4q4+FQpcJyt\njd8Xj/L7OHPreQL83v78Ht/vC+dI9nXg99UW8/m8KW+dv9+P3+/xrfye/1S8BtnDAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAicIAAACcwn5PAP/k97XcreMHg7buGQqFfB3f+viteevj93v7sSoUCuafkc/nfZ2D\n53mmvHUb8Hv+pXgO/WSd/3R9DU7mcbOHAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOF\nAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2eQKWwXsv8\ncs+Hw7ZNKRi0dddQKGTKW+dvHd+at64/az6Xy5ny+XzelJckz/N8zVsfg9/zt74HFAoFU946f7/H\nt64/P8cvdt2xhwEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAA\nOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOIX9nsC/sl5T3C9+zzsUCpnywaCtO1rH\nD4dtm2IkEjHlq6qqfB3f+viteatcLmfKZ7NZ8xzy+byvc7COb12HnueZ8tb5W/PW91C/x7fyc/xC\noVDU/djDAAAAnIr+s6Snp0e7du3SRx99pNOnT6u+vl4LFizQmjVr1NLSMnq/p59+Wu+88864/Ny5\nc/Xyyy+XZtYAAGBKFV0Ytm3bpiNHjqizs1Otra2Kx+PasWOH1q5dq82bN+vaa68dvW8kEtH69evH\n5Gtqako3awAAMKWKLgyrVq1SW1vbmM+qlyxZotWrV2vbtm164okn/v8PDYe1bNmy0s4UAAD4pujv\nMLS3t4/7YtucOXN0zTXX6MSJE2OWFwoFeZ6nwcHB0swSAAD4yvTV6kKhoHPnzqm1tXXM8nQ6rbvv\nvlvpdFqzZs3S0qVLtXbtWlVXV5smCwAA/GEqDHv27FFfX5/WrFkzuqypqUn33XefrrvuOnmep0OH\nDumtt97SsWPHtGnTJvPhdwAAYOpdcmE4ceKEnnvuObW3t+uOO+4YXd7V1TXmfkuWLFFLS4teeukl\n7du3T0uXLr302QIAAF9c0nkY4vG4Hn/8cdXW1mrjxo3OE07ce++9CgQC+vDDD50/u1AojPsPAADY\neJ6nfD4/5r/JnPBr0nsYksmkNmzYoMHBQT3//PNqbGx0ZqLRqOrq6pRIJJz39ftsWwAAfBVNdEbf\nkYMUijGpwpDJZPTkk0/q1KlT+sUvfqGrr766qFwqldLAwIBisdhkhgMAABWi6I8k8vm8Nm7cqO7u\nbj311FNauHDhuPtkMhmlUqlxy1955RVJ0o033miYKgAA8EvRexi2bNmigwcPqqOjQwMDA9q9e/eY\n25cvX654PK6uri7dfvvtmjt3riTp/fff16FDh3TTTTfp1ltvLe3sAQDAlCi6MBw7dkyBQEAHDx7U\nwYMHx9wWCAS0fPly1dbWqqOjQ4cPH9auXbvkeZ5aWlrU1dWlVatWlXzyAABgahRdGDZt2uS8T21t\n7ZhTRAMAgK8G04mbvkqsR2dM9/xE356djHDYtilVVVWZ8jNmzPA1b7242syZM0156wnRrM//RN9d\nmox0Om3Kl+Jn5HI5X8fPZrOmvHX+1vGtJnN4XzlYXwP5fL5EM7k0U3EKAtsaAgAAlwUKAwAAcKIw\nAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAA\nACcKAwAAcKIwAAAAp7DfEygV67XMA4FAiWZyaazzD4VCpnwkEjHlq6qqTPkZM2aY8rNmzTLl6+rq\nTPmGhgZTvqamxpS3bj/WfC6XM+WHhoZMeUlKpVKmfCKRMOWHh4dN+WQyacpnMhlT3voeaM3n83lT\n3roNep5nyvv9O2QqsIcBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADiF/Z7AlwUCgWl7TfFQKGTKB4O2\n7hYO257KSCRiys+YMcOUj8VipnxjY6Mp39TUZMo3NDSY8rW1taa8dfux5rPZrCmfTqdNeUk6f/68\nKV9XV2fKx+NxU976HCSTSVPeb4VCwZS3vgfmcjlT3iqfz/s6fjHYwwAAAJwoDAAAwInCAAAAnCgM\nAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAA\nwInCAAAAnGwXEK8ggUBgWueDQVt3s+YjkYgpP3PmTFM+FouZ8rNnzzblv/a1r5nyV1xxhSlvffzR\naNSUt24/nueZ8gMDA6a8JPX395vyvb29prx1HVoVCgVfx7duA9b553I5U946f7/l8/lLyk1mvbOH\nAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEA\nADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4hf2ewJcFAgEFAgFfxg6FQqa8dd7hsO2piEQipnw0GjXl\nZ86cacrHYjFTvqmpyZSfM2eOKT937lxTvrm52ZSfMWOGKW/dfguFgimfSCRMeUk6deqUKW99DVnX\nQS6XM+Wz2awpn8/nfc1b15/neaa89T3Yuv6nA/YwAAAAp6IrVU9Pj3bt2qWPPvpIp0+fVn19vRYs\nWKA1a9aopaVlzH2PHz+uzZs365NPPlEkEtHNN9+sRx55RPX19SV/AAAAoPyKLgzbtm3TkSNH1NnZ\nqdbWVsXjce3YsUNr167V5s2bde2110qSent7tW7dOs2aNUtdXV1KpVJ644039I9//ENbtmwx7/YB\nAABTr+jf3qtWrVJbW9uYz/qXLFmi1atXa9u2bXriiSckSa+99prS6bS2bt06+rlsW1ub1q9fr507\nd+qee+4p8UMAAADlVvR3GNrb28d9MXDOnDm65pprdOLEidFl+/fvV0dHx5gvcS1atEgtLS169913\n7TMGAABTzvSlx0KhoHPnzo1+N6G3t1f9/f2aP3/+uPu2tbXp73//u2U4AADgE1Nh2LNnj/r6+rRk\nyRJJUjwelyQ1NjaOu29jY6MSiYT50CEAADD1LrkwnDhxQs8995za29t1xx13SJLS6bSkiY/pH1k2\nch8AADB9XNIhC/F4XI8//rhqa2u1cePG0ZO+VFVVSZIymcy4zMiykftMZKITbwQCAQWDnC4CAACL\nfD5vOsHVpAtDMpnUhg0bNDg4qOeff37Mxw8j/z/y0cSXxeNx1dXVXfSwymAw6NuZHgEA+CoLhULj\nDl7wPK/orwpMqjBkMhk9+eSTOnXqlH7xi1/o6quvHnN7c3OzYrGYenp6xmV7eno0b968yQwHAAAq\nRNH7+vP5vDZu3Kju7m499dRTWrhw4YT3u+222/Tee++pt7d3dNnhw4d18uRJdXZ22mcMAACmXNF7\nGLZs2aKDBw+qo6NDAwMD2r1795jbly9fLkm6//77tW/fPj366KNauXKlUqmUtm/frtbWVt11112l\nnT0AAJgSRReGY8eOKRAI6ODBgzp48OCY2wKBwGhhaG5u1q9+9Su9+OKL+s1vfqNoNKpbbrlFDz/8\nMKeFBgBgmir6N/imTZuK/qFf//rX9eyzz17ShAAAQOXhT/4SsR7dYc1bDz292OGuxaitrfU1H4vF\nTPkvn8r8UsydO9eUb2hoMOVnzpxpyv/rN6cnK5/Pm/KpVMqUl+yPYaLDwSdjaGjIlM9ms6a89aR4\n1ufQ7/EthwuWIu/3e/ilCgQCRT93nOAAAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UB\nAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADiF/Z5ApbBeC916\nLXO/r6UeDts2haqqKlM+EomY8jNmzDDlY7GYKV9TU2PK19fX+zq+df0XCgVTfmhoyJSXpHw+7+sc\nrPl0Om3KZ7NZX8cfHh72dXzre6DfeevvoKnAHgYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEA\nAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4BT2ewKV\nwu9roVtZxw+H/d0UQqGQKV8oFEx5v69Fb338kUjElK+qqjLlrayPX5IaGxtN+cHBQVP+7Nmzpnwi\nkTDlrfOvrq425a3bUDQaNeWz2awpP91/B1zqe+BkcuxhAAAAThQGAADgRGEAAABOFAYAAOBEYQAA\nAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABO\nYb8nUCqBQOCyHt/zPFM+l8v5On4mkzHl0+m0KZ9MJn0d3/r4J3NN+3IIh21vJaWYf01NjSlfX19v\nyl9xxRWmfH9/vylvffwzZsww5aurq035wcFBUz4YtP396/d7uN/jF4M9DAAAwInCAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAyXYRe4zyPM+ULxQKvo6fy+VM+VQqZcqn02lTPplMmvLnz5835QcGBkz5WCxmylvX\nX1VVlSlv3X6DQfvfLuGw7e3Mug6sz2FdXZ0pH41GTflIJGLKW9e/NR8IBEz56e5SH/9kcuxhAAAA\nTpOqdENDQ3r99dfV3d2tnp4eJZNJPfbYY7rzzjvH3O/pp5/WO++8My4/d+5cvfzyy7YZAwCAKTep\nwjAwMKBXX31Vs2fP1rx58/Txxx9fcHdGJBLR+vXrxyyrqam59JkCAADfTKowNDU16c0331RDQ4M+\n/fRTPfTQQxf+weGwli1bZp4gAADw36S+wxCJRNTQ0FDUfQuFgjzP0+Dg4CVNDAAAVI6yHSWRTqd1\n9913K51Oa9asWVq6dKnWrl2r6urqcg0JAADKpCyFoampSffdd5+uu+46eZ6nQ4cO6a233tKxY8e0\nadMmhUKhcgwLAADKpCyFoaura8y/lyxZopaWFr300kvat2+fli5dWo5hAQBAmUzZeRjuvfdeBQIB\nffjhhxe8j+d5yufzY/6znpAIAAD88wR96XR6zH/ZbLbo/JSd6TEajaqurk6JROKC9wkGg5f92boA\nACiHcDg87oyanucVfabYKdvDkEqlNDAwYD59KgAAmHolLwyZTGbC6wq88sorkqQbb7yx1EMCAIAy\nm/RHEjt27FAymdTZs2clSQcOHNCZM2ckSStWrFAikVBXV5duv/12zZ07V5L0/vvv69ChQ7rpppt0\n6623lnD6AABgKky6MLzxxhs6ffq0pH9e5eovf/mL9u/fr0AgoG9/+9uqra1VR0eHDh8+rF27dsnz\nPLW0tKirq0urVq0q+QMAAADlN+nCsG3bNud9nnjiiUuaDAAAqExc3hoAADhN2WGV5VYoFEx56/ke\ngkFb9/J7/plMxpQv9rCcCxkaGjLlrdcs6e/vN+XPnTtnyluPHrJeCdZ6ynbr2Vut279kfw1GIhFT\nPhqNmvLWdWidv3X86X5IvHX+1m3Yr/xkcuxhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOYb8n8GWW\n64FP92ux5/N5X/PZbNaUHx4eNuXPnz9vytfU1JjyiUTClI/H46b8FVdcYcqnUilTvra21pSPRqOm\nfClev9afEYlETHnL+1cpxg+FQqa8df0Fg7a/P/1+D7c+f5cD9jAAAAAnCgMAAHCiMAAAACcKAwAA\ncKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCi\nMAAAAKew3xP4qrBeS92a9zzPlM/n86Z8JpMx5YeHh30dP5lMmvLnz5835QcGBkz5pqYmU966/qzb\nTyQSMeUlKRQKmfLBoO3vJ+v4gUDAlLfye3y/3wPhxh4GAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOAU\n9nsCXxYIBHy7Jrv1Wux+563Xgs/n86Z8JpMx5YeHh0358+fPm/IzZ8405QcHB0156/yHhoZMeev2\nY33dhsP2tyLrY4hGo77mQ6GQKW99DVufw2w2a8pb529lHX+6/g6YTI49DAAAwInCAAAAnCgMAADA\nicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInC\nAAAAnCgMAADAyX4R+hKyXg/cwnot+Eu9FvmIUChkylvXnfVa8Nls1pRPpVKmfCKRMOVnzZplyieT\nSV/zmUzGlLc+f9bt37r9VQLrOrAKBm1//1m3oVwu52ve+h5off78fv6nwqQKw9DQkF5//XV1d3er\np6dHyWRSjz32mO68885x9z1+/Lg2b96sTz75RJFIRDfffLMeeeQR1dfXl2zyAABgakyqkg4MDOjV\nV1/V//zP/2jevHmSJv7LvLe3V+vWrdPnn3+urq4urVq1Sn/961/14x//2NwiAQDA1JvUHoampia9\n+eabamho0KeffqqHHnpowvu99tprSqfT2rp1q5qbmyVJbW1tWr9+vXbu3Kl77rnHPnMAADBlJrWH\nIRKJqKGhwXm//fv3q6OjY7QsSNKiRYvU0tKid999d9KTBAAA/ir5URK9vb3q7+/X/Pnzx93W1tam\nv//976UeEgAAlFnJC0M8HpckNTY2jrutsbFRiUSC7zEAADDNlLwwpNNpSVI0Gh1328iykfsAAIDp\noeTnYaiqqpI08TG9I8tG7vOvCoXCuGNpA4GA+RwJAABc7vL5vOl8ESUvDCMfRYx8NPFl8XhcdXV1\nCocnHpZyAABAeYRCoXEnCfQ8r+ivCZT8I4nm5mbFYjH19PSMu62np2f0/A0AAGD6KMu1JG677Ta9\n99576u3tHV12+PBhnTx5Up2dneUYEgAAlNGkP5LYsWOHksmkzp49K0k6cOCAzpw5I0lasWKFampq\ndP/992vfvn169NFHtXLlSqVSKW3fvl2tra266667SvsIAABA2U26MLzxxhs6ffq0pH9+5+Avf/mL\n9u/fr0AgoG9/+9uqqalRc3OzfvWrX+nFF1/Ub37zG0WjUd1yyy16+OGHL/j9BQAAULkm/dt727Zt\nRd3v61//up599tlJTwgAAFSesnyHAQAAfLXw+cD/Y72WupX1Wur5fN6UDwZt3dE6fjabNeUHBwdN\n+fPnz/uaP3funCnf399vyg8NDZnyw8PDpnwpDqf2exuynpAumUz6mrfO33oG34nO3TMZfr8HWll/\nB03F7zD2MAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoD\nAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAp7DfE8A/Wa9l7nmer+PncjlTPpPJmPLpdNqU\nTyaTpvzAwIAp39fXZ8qfPHnSlP/a175mygeDtr89qqurTXlJyufzpnx/f78pf+bMGVO+t7fXlE8k\nEqb80NCQKZ9KpUx56/NnzWezWV/Ht+anAnsYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAA\nOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIBT2O8JVIpC\noeDr+J7nmfLBoK375XI5U94qFAqZ8plMxpRPJpOm/ODgoCkfj8dN+f/93/815Wtra0156/qzji9J\n+XzelD937pwpf+LECVP+7NmzpnwikTDlrY8/nU6b8tls1te89T3Y798hlzr+ZHLsYQAAAE4UBgAA\n4ERhAAAAThQGAADgRGEAAABOFAYAAOBEYQAAAE4UBgAA4ERhAAAAThQGAADgRGEAAABOFAYAAOBE\nYQAAAE4UBgAA4ERhAAAATmG/J/Blfl5PPBAI+Da2ZH/suVzOlA+HbZtCPp835a3Xsg8Gbd3Xmj97\n9qwpb13/1rzneab8mTNnTPmamhpTXrI/hoGBAVPeug34nT9//rwpPzg4aMpb3wOsz781b30P9ys/\nmRx7GAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4\nURgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAU9jvCZSK9Vrilzvr+svn86Z8MGjrrplMxpQPBAKm\nfH9/vylv5XmeKT88PGzKnzlzxpSvrq425SX7Njg4OGjKDw0NmfLxeNyU7+vrM+WTyaQpb30NZrNZ\nUz6Xy5ny1teQNT8dfoeVpTB8/PHH+tGPfjThbZs3b9aCBQvKMSwAACiTsu5hWLlypebPnz9m2VVX\nXVXOIQEAQBmUtTBcf/31Wrx4cTmHAAAAU6CsX3osFApKpVLmzxYBAIC/yrqH4dlnn9XQ0JCCwaC+\n+c1v6oEHHhj3EQUAAKh8ZSkMkUhEixcv1i233KL6+np99tln2r59u9atW6cXXnhB8+bNK8ewAACg\nTMpSGNrb29Xe3j76746ODi1evFg//OEPtXXrVj3zzDPlGBYAAJTJlJ2HYc6cOfrWt76l/fv3q1Ao\nXPC494mORbUeIw8AwOXO8zzT+R6m9MRNzc3NyuVyGh4evuCJWigHAACU3kQnyCsUCkWfdGpKTw39\n+eefq6qqqiRndQMAAFOnLIVhotPkHj16VAcOHNANN9xQjiEBAEAZleUjiZ/+9KeqqqpSe3u7YrGY\njh8/rrffflvV1dXq6uoqx5AAAKCMylIYbr31Vu3Zs0d/+MMflEqlFIvFtHjxYv3gBz/g1NAAAExD\nZSkMK1as0IoVK8rxowEAgA+m9EuPAABgeprSwyqLcanHiFoPx5wO1yIvJ+u15EOhkK/jW/m9/Vjz\n1vWXTCZN+aqqKlO+Eg6nHh4eNuWHhoZM+VQq5WveOv9MJuNrvthDAy/Ees0jv98DLjU/mRx7GAAA\ngBOFAQAAOFEYAACAE4UBAAA4URgAAIAThQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAU9jvCZSK9VrigUDA1/GnO+u15K2s69/zPF/z1vlb138qlTLl\nI5GIKR8M2v92sa6DXC5nymez2Wmdz2Qyprx1/VlfQ9bn3/oa9Ds/FdjDAAAAnCgMAADAicIAAACc\nKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgMAADAicIAAACcKAwAAMCJwgAAAJwoDAAAwInCAAAAnCgM\nAADAicIAAACcwn5PoFQCgYApb70WuXV8q+lwLfWLsV7L3m/W9W/NZ7NZUz4SiZjy1u3f79ePJHme\n52s+l8uZ8tZtwLoNWl/D1rzfr0Hr8+/Xe/hkxmUPAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCi\nMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwCvs9\ngVLx61riIwKBgCnv9/z9Nt0fv3X+1nwwaOv+uVzO1/ErgfU5yOfzvo7ved5lnff7Nei3qZj/9H+V\nAwCAsqMwAAAAJwoDAABwojAAAAAnCgMAAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwojAAAAAnCgMA\nAHCiMAAAACcKAwAAcKIwAAAAJwoDAABwCvs9ga8K67XcA4FAiWbiD+u12K2P37r+rfMPBm3d2zq+\n39tfJWy/1jlYnwO/8/l83tfxL/e8ld/jF6MshSGTyeh3v/uddu/erWQyqdbWVq1Zs0aLFi0qx3AA\nAKDMyvKRxDPPPKM//vGPWr58uf7jP/5DwWBQP/nJT/S3v/2tHMMBAIAyK3lh6O7u1t69e9XV1aUH\nHnhAd999t375y19q9uzZ+vWvf13q4QAAwBQoeWHYt2+fQqGQvvvd744ui0aj+s53vqMjR46ot7e3\n1EMCAID62ujdAAAJcElEQVQyK3lhOHr0qFpaWlRdXT1m+fz580dvBwAA00vJC0NfX5+amprGLR9Z\n1tfXV+ohAQBAmZW8MGQyGUUikXHLo9GoJCmdTpd6SAAAUGYlLwzRaFTZbHbc8kwmI0mqqqoq9ZAA\nAKDMSl4YmpqadPbs2XHLRz6KmOjjCgAAUNlKXhjmzZunkydPKpVKjVne3d09ejum3nQ4i1gls55J\n8XKXy+X8nsK0xvqz4f2vNEpeGDo7O+V5nt5+++3RZZlMRjt37tTChQvV3Nxc6iGBsuMNx4ZfeDbW\n0z5f7nj9lkbJTw29YMECdXZ2auvWrTp37pyuuuoq7dq1S2fOnNGGDRtKPRwAAJgCZbmWxOOPP67f\n/va32r17txKJhL7xjW/o5z//ua6//vpyDAcAAMqsLIUhGo3qwQcf1IMPPliOHw8AAKZYRVzeeuSQ\ny8vZVHzGxud4l65QKFx0/fn9pchKvzx1oVAo++fwfl/e2so1vmsb83sb9Pvy0q689favumJ+D1dE\nYfjiiy/8nsJlz+8Xi9/jF4MvntkMDw/7PYVpjZPe2UyH9xg/ffHFF/q3f/u3i94nsHfvXt/XYn9/\nvz744ANdeeWVo2eEBAAA5ZXJZPTFF1/ohhtuUCwWu+h9K6IwAACAylby8zAAAICvHgoDAABwojAA\nAAAnCgMAAHCiMAAAAKeKOA/DhWQyGf3ud7/T7t27lUwm1draqjVr1mjRokV+T63iffzxx/rRj340\n4W2bN2/WggULpnhGlWtoaEivv/66uru71dPTo2Qyqccee0x33nnnuPseP35cmzdv1ieffKJIJKKb\nb75ZjzzyiOrr632YeWUodv09/fTTeuedd8bl586dq5dffnmqpltxenp6tGvXLn300Uc6ffq06uvr\ntWDBAq1Zs0YtLS1j7sv2N16x64/tz66iC8MzzzyjP//5z/re976nlpYW/dd//Zd+8pOf6Je//CXX\npSjSypUrNX/+/DHLrrrqKp9mU5kGBgb06quvavbs2Zo3b54+/vjjCc8a2Nvbq3Xr1mnWrFnq6upS\nKpXSG2+8oX/84x/asmWLwuGKfjmVTbHrT5IikYjWr18/ZllNTc1UTLNibdu2TUeOHFFnZ6daW1sV\nj8e1Y8cOrV27Vps3b9a1114rie3vQopdfxLbn1XFbmHd3d3au3evHnzwQa1atUqStHz5cq1evVq/\n/vWv9cILL/g8w+nh+uuv1+LFi/2eRkVramrSm2++qYaGBn366ad66KGHJrzfa6+9pnQ6ra1bt45e\npr2trU3r16/Xzp07dc8990zltCtGsetPksLhsJYtWzaFs6t8q1atUltbm0Kh0OiyJUuWaPXq1dq2\nbZueeOIJSWx/F1Ls+pPY/qwq9jsM+/btUygU0ne/+93RZdFoVN/5znd05MgR9fb2+ji76aNQKCiV\nSnFa44uIRCJqaGhw3m///v3q6OgYfbOWpEWLFqmlpUXvvvtuGWdY2Ypdf9I/t0fP8zQ4OFjmWU0f\n7e3tY37ZSdKcOXN0zTXX6MSJE6PL2P4mVuz6k9j+rCp2D8PRo0fV0tKi6urqMctHdq8fPXp0zAsH\nE3v22Wc1NDSkYDCob37zm3rggQfGfUQBt97eXvX390+47tra2vTXv/7Vh1lNP+l0WnfffbfS6bRm\nzZqlpUuXau3ateNe55e7QqGgc+fOqbW1VRLb32T96/obwfZnU7GFoa+vT01NTeOWjyzr6+ub6ilN\nK5FIRIsXL9Ytt9yi+vp6ffbZZ9q+fbvWrVunF154QfPmzfN7itNKPB6XJDU2No67rbGxUYlEQrlc\n7rL9HLkYTU1Nuu+++3TdddfJ8zwdOnRIb731lo4dO6ZNmzaN+yvxcrZnzx719fVpzZo1ktj+Jutf\n15/E9lcKFbt1ZTIZRSKRcctHLk7Fldsurr29Xe3t7aP/7ujo0OLFi/XDH/5QW7du1TPPPOPj7Kaf\nke1tooujfXmb5A37wrq6usb8e8mSJWppadFLL72kffv2aenSpT7NrLKcOHFCzz33nNrb23XHHXdI\nYvubjInWn8T2VwoV+x2GaDSqbDY7bvnINburqqqmekrT3pw5c/Stb31LH330EZd6naSR7W2ia8az\nTV66e++9V4FAQB9++KHfU6kI8Xhcjz/+uGpra7Vx48bRo03Y/opzofV3IWx/k1OxdbSpqUlnz54d\nt3zko4iJPq6AW3Nzs3K5nIaHh/ncbhJGdgWP7Br+sng8rrq6Ov66uwTRaFR1dXVKJBJ+T8V3yWRS\nGzZs0ODgoJ5//vkxHz+w/bldbP1dCNvf5FTsHoZ58+bp5MmTSqVSY5Z3d3eP3o7J+/zzz1VVVUVZ\nmKTm5mbFYjH19PSMu62np4ft8RKlUikNDAwoFov5PRVfZTIZPfnkkzp16pR+/vOf6+qrrx5zO9vf\nxbnW34Ww/U1OxRaGzs5OeZ6nt99+e3RZJpPRzp07tXDhQo6QcOjv7x+37OjRozpw4IBuuOEGH2Y0\n/d1222167733xhzSe/jwYZ08eVKdnZ0+zqzyZTKZceVfkl555RVJ0o033jjVU6oY+XxeGzduVHd3\nt5566iktXLhwwvux/U2smPXH9lcaFbsPa8GCBers7NTWrVt17tw5XXXVVdq1a5fOnDmjDRs2+D29\nivfTn/5UVVVVam9vVywW0/Hjx/X222+rurp63Jd/IO3YsUPJZHL0Y7ADBw7ozJkzkqQVK1aopqZG\n999/v/bt26dHH31UK1euVCqV0vbt29Xa2qq77rrLz+n7zrX+EomEurq6dPvtt2vu3LmSpPfff1+H\nDh3STTfdpFtvvdW3uftty5YtOnjwoDo6OjQwMKDdu3ePuX358uWSxPZ3AcWsv3g8zvZXAoG9e/dW\n7LffMpmMfvvb32rPnj1KJBL6xje+odWrV/MXchH+8z//U3v27NGpU6eUSqUUi8X07//+7/rBD37A\nqaEn8P3vf1+nT5+WpNEvShUKBQUCAf3+97/X7NmzJUmfffaZXnzxRf3tb39TNBrVzTffrIcffviy\n36XpWn81NTV6/vnn1d3drbNnz8rzPLW0tGjZsmVatWrVZX1I26OPPqr//u//nvCLyIFAQH/6059G\n/832N14x6y+ZTLL9lUBFFwYAAFAZKvY7DAAAoHJQGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIAT\nhQEAADhRGAAAgBOFAQAAOFEYAACAE4UBAAA4URgAAIDT/wGKnqNabRh3HQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10c05fd90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "cutout2 = Cutout2D(data, coord, (30, 30), wcs=w)\n",
    "plt.imshow(cutout2.data, cmap='Greys_r')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArAAAAHHCAYAAABHi2tSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJztnXuwVdWd57/nvgBFLnJFMhHzEBCQcarHR5Pb1QHtiRPT\nSIHQobQqme6gNDFW2iHVJubRVWSquida1d0yBqNj2Vai1Ua67YxUmDaRDl41GtMxZDAloPjAR4ni\nBa5c7uU+z/xBzvHcc89jP9bj91vr+6lKRc7dZ+/fe/32PmuvVdi1a1cRhBBCCCGEKKHFtwCEEEII\nIYSkgQ0sIYQQQghRBRtYQgghhBCiCjawhBBCCCFEFWxgCSGEEEKIKtjAEkIIIYQQVbCBJYQQQggh\nqmADSwghhBBCVNHmWwCbHDt2DL/61a/woQ99CB0dHb7FIYQQQgghVQwPD+PQoUO45JJLMHPmzETf\nCbqB/dWvfoW//uu/9i0GIYQQQghpwje/+U186lOfSnRs0FMIPvShD/kWgRBCCCGEJCBN35b4Cey+\nffvwk5/8BLt378Y777yDzs5OLF68GNdddx3mzp074diDBw9i69at+O1vf4v29nYsXboUN954Izo7\nOyedd8eOHdi2bRsOHTqEs88+G2vWrMHVV1896bj+/n7cddddeOqppzA0NITFixfjhhtuwIIFC+rK\nnGXaQKFQQLFYNH5szEiykyRZbGBKvxDtlFenEG3iC9qSEDdoy7U0fVviJ7APPvggnnrqKVxyySX4\n8pe/jKuuugp79uzBn//5n+PVV18tH3f48GHcdNNNePvtt7FhwwasW7cOzz77LP7yL/8So6OjE865\nfft2/O3f/i0+/vGP4y/+4i9wwQUX4I477sCDDz444bjx8XHccsst+NnPfoY1a9Zg48aNOHr0KDZt\n2oS33norsbJJSONoTUHhE0l2yitLoVAwJIkdTNk663na29uNXN8GeW1T+r70GNCApJqgmXqx6CpG\n017H1rsozMn6VOdaM1tpsmXiJ7Dr1q3DokWL0NraWv7s8ssvx/r16/Hggw/iG9/4BgDggQcewNDQ\nEO655x7Mnj0bALBo0SLcfPPNePTRR3HVVVcBAIaGhnDvvfeiu7sbmzdvBgCsWLECxWIR999/P1au\nXInp06cDAHp6evDCCy9g8+bNWLZsWfnan//853HffffhW9/6Vn5LkFRou6szRTOd89glBJuOjIz4\nFmESpu2qyUeVg5FEuUOIeVvksY0rm6a9zvDwsAg5YqaZrTTZMvET2CVLlkxoXgHgnHPOwUc/+lG8\n/vrr5c+efPJJdHd3l5tXALj44osxd+5cPP744+XPdu/ejePHj2PVqlUTzrl69WqcPHkSzzzzTPmz\nnp4ezJo1q9y8AkBnZycuu+wy/PznP5/0ZDcNhUIh0x2HprsUGyQJ8hhtlCf5XRQO+iQfhUIBLS26\nXx1oFAOuddM0WLomiW1oP13EWH9tkqtaFYtFHD16tDy39fDhwzh27BgWLlw46dhFixbhpZdeKv/7\nwIEDADDp2AULFqBQKODll1+ecGytua4LFy7E0NAQ3njjjVw6ZCkCLguH1qBncZWH6WYuNorFIsbH\nx32LkYtisVh3qodL3fLEj8bY0yhz7Jj2maYxUeo0lEpyNbA7d+5Eb28vLr/8cgDAkSNHAACzZs2a\ndOysWbNw/Pjx8tPS3t5etLS0THqxq729HZ2dnXjvvffKn/X29qKrq2vSOUuf9fb25lFDPGyWiUuq\nf2mph+tinHXuluR5uUnJmpclH1X7SsJUj0bx00xfTY1AibQyS4vbGMcGjXGWlEb+dDm1J891Mjew\nr7/+OrZs2YIlS5bg05/+NIBT81qB2hO1S5+VjhkaGqqboO3t7RPmygwPD9c8tvqcrpCYyKZkyjN/\ns5paNzJEPmNjY75FqEnWuVsSmrW85HmBTOMg7OOFOWkvt2SN22aNSVY0xhGpTyN/avF1pgb2yJEj\n+PrXv47p06fj29/+djkppkyZAqD2RO3SZ6VjpkyZUjdBh4eHJzTB9Y6tPmdekia3ROfWksllwa11\n/dIT+VBIak/tcyTTIPFmLmR81h4fvnapbygvt4TQmJD0xFiLU4+0/f39+NrXvoYTJ07gtttum/CU\nrfTftRqXI0eOYMaMGWhrO7XwQVdXF8bHx9HX1zfhuJGREbz//vs466yzyp91dXVNmFJQojR1oNb0\ngiyEltw+9DGVRBKTMak9tc+RTENoOUPqI9XXJmpFadqMpLrj+umzJN1JeqTmp01SNbDDw8P45je/\nibfeegt/8zd/g4985CMT/j579mzMnDkT+/btm/Tdffv2Yf78+eV/l/67+tj9+/ejWCxi3rx55c/m\nzZuHl156aZKD9u7di6lTp+Lcc89No4Z4kj7By1twbBQs3+uQ2kCinaURg45EJiZqRWnajIu6U3qI\n0wxJNZCEgaQpOSZI3MCOjY3h29/+Nvbu3YvNmzfjggsuqHncJz/5SfziF7/A4cOHy58999xzePPN\nN7F8+fLyZxdddBHOOOMMbN++fcL3t2/fjqlTp6K7u7v82fLly3H06FE88cQT5c/6+vrQ09OD7u7u\nxAVBC0mf4JlamJ00hnZuTgw6knCwNbgmOW+eZR9JPmK/0ZY0JccEiTu/733ve3jmmWfQ3d2Nvr4+\nPPbYYxP+fsUVVwAAPve5z6GnpwebNm3C2rVrMTAwgIceegjnnXcePvOZz5SP7+jowPr167FlyxZs\n3rwZl156Kfbs2YOdO3fi+uuvL29iAJxqYB9++GHcdtttOHjwIGbMmIFHHnkExWIRX/jCFxIr63vR\nbG2L3Pu2l1ZoNz20tbU1bCik+LKeHFLky4ov+W1ds955Nfip2aoQ0uVPQgg6SEBKPBR27dqVSIpN\nmzZhz549dYvov/3bv5X//dprr+HOO+/E888/j46ODixduhRf+tKXMHPmzEnf3bFjB7Zt24a3334b\nc+bMwerVq7F27dpJx/X39+Ouu+7CU089haGhISxevBhf/OIXcf7559eV+cUXX8TGjRuTqEcE0Cwp\nsiSNiUSTkKy+ZfB9/Xo0kkuqzFkISRdCtMH8c8fdd9/dsK+rJHEDqxGpDayEZGj25ImEiYTY80Gs\nepP4YKwTzaRpYONZ78cxjebaSCgubF7TEcrcKQmx54PQ9A4lHol5Qov1WjD+CcAGtiEmF31mwvmj\n+iW/WBaDzwPjVTaxxWNeGM/uMf1ydaUPGf/ZCG1FHTawDTCZJI3m6aXB5C4r0oLRFtVPm6UVP5Nb\nRvrekS1WQsmlUPSohvHsniS/8qWJN/owP0ltWM8v0nzABtYzaQPC5C4r0oLRNtVJaWqwLi2Cnuba\nlZjc6tT3jmyxIiWX8vpaih5aYG7lg/EmEy1+YQObg3rFi0VNJtVJaSpJS4ugp7m2S7QUIylozt9i\nsaha/mZI0425RUgybOQuG9gc1CteNoqatMJtG5P6hjbvp1Ke0HSTjBZbhdxUhawbyY+WHDVNM70l\nbPZkI3fZwCohtsLtYv6xq++bLqqV8vjSLel2xxpxeWNKiEakNoqactSkDZutnx7qqkPhjkKekJrY\npjH54lHoaCqqSUm63TEhoZCktsdS/0Osaa6ptKHNuAnZV2xgDeMjWHwUTZMvHqUli761vhPLYFNJ\njDqbJhYbxqJnUpLUdi3NAn0rC1NxI/WXMZPT3iqRqS1JRZLg11awTG8EUes7WbalNYFPX8Ryg2UT\nLU1KXmLRMwnTpk3zLYJRmv3knCdnbeR7aDXEFNV2kfrLmMlpb5VE18CaXEdVE7YHI9N3fibktTn3\nVMJ5pNAspzTqayJ2Qqsn2vQxKe/g4KCxc6XFtd2LxWKu9UJt5LvGGuKC2O0SXQNrch1VLbgogBLv\n/EL1pzSy5JT0ZshE7IS2Gom2fNImbz0k6yFZNuIGnzUpugY2RlhkiCli3ulLQtOt0W6S4ZP2uIjR\nV7Z1zlOT8srGBtYQpjc1iDHRiHzYQJFqNNcqqU/aiR00+yrEm628/mADW0VLS0umN+ZMrx0pNdE4\nQZ+ESrM4lJqTeTCReyHaJSZYf3Vgopfw9eKyLaJqYJM4Y3x83Nobcxqptlk9e+TZ6cOUjaUnmwRo\no/q4znUJvoi5vrm0v9TljYC4Y4A0RvoKR3KzygJM1PQktZntnT6SJEko/uWi1nEgcUkzk8snSW7a\nALf2l/iSq08k3LwRM/gcU2RXGI9IS7C0u8BIkz8vMTVemnQNLc5CI+kvKEn/nua7vps2xqZcfNY4\nzXGhWXYbsIGtg7QmIu0uMNLkNw0TWQahx1kaJK4xHbN/YtCddTA9muNCs+w2YANLVBJCInPwSYfp\nn6RdbnQRQrwSeYQeV+3t7b5F8ArHiMZE1cAyGIgkQh98TFP6STrmtWg1w/obDq58OTIy4uQ6jeDm\nIXKJqoFlMBCiH+axTug3fdRr3mLypUZdY7lZjKqBJclw8faw9DeUTRFLISHuaG1t9S2COJhndsjb\nvNEvfsjqN23+iqOLIHWpFbBJ3x7OE+y+31B2hca79xLaipk0bNlvbGzMynl9ktdWmvNMI0lvougX\nXRSLRVV1nw1sDTQ5MC1pl9VpRNrvZnlypNkXmmUHOPjkJQT7Jd2gxFcDqj3HtCL1JorxkB9NdYsN\nbAWl4K/lwHqJ4Tph8v70bjI40+qepehJS6Y0OkuT3TUcTOSR1idJNyjxFeux51hIVMdmlvrBeIiL\nKBrYpE8RsiyD43qOkKSf3mMsFll0jrWRS2qrWO3jgxhzVhqM94mUHspUxyZjlTQjigbW9janeWCS\nhg993BitL4qwESFZCKEemIx9SQ9liC6iaGDTwoGpNrSLOdJuDRwLmn42tHndGH1PTuF6hRZNOUdI\nJWxgaxBTcpaKV6MidtpppwGIyy62Sbs1cC1CbHIYY6egHeLF9RNJzbEWYg2MAVN+i6qBTToXVjKm\nHF/9wlqjRnZgYMDINW0SYyHTPPAQQsIhbf3lbnpxY8pvUTWwkufCJsXUU7nq85Tu+lkQdJF1IIix\n4TcFbUfIRNKOGxxnkpG31oReq6JqYGPAZ2FIMh3BFtVPkpshJbF9LYumYQCR4qNqNNjOBqbmZvr2\na+X1pe8I6MJWnI8vl7y1JvRaJTt7SSKkFJdSsvhMmqTXlpLYGt/AzRpv7e3tqY5P6yNtA7G2KU2m\nYtV37lVeX3r+ubCVifn4aZCUg0Q3bGAt4ipRfQ8ILgit6GnWJ2u8jYyMGJZkImkGYgn29zWlSYLu\nJF5iGK+IG9jAGqDegKApUbNs8+oSTbZMQmj6aEOS/V03lJJ0901ozXxo+oQMfZV/Cg8bWAOEMCDk\n3duayUhINkKoHxopFArB2T40fWwgZeMT+ir/FB42sMQIaZMxzZ1XTM2xTV19PGWX/mQ/dmLKrWok\nTSexQT29XOkr1a4hbnziG1++ZgNLmmIjONPceVUnvutkcXk9m0Uu71N2X9dkE2yPkAfVpGSxgfTV\nC4D6ernyOWMrHnz5Wn4WEu9IK0Su5ZGmvxTSNPZ5Fjr30XgTOUh8Yih99QKXSH3SmpdQ9bKJa5ux\nga2CQUtIMtI09qEsdO6zPvi4toR6aDMWKp/uS4056SSxm4Q4SktM8aB1ZzQ2sFXYWHvSBPwZlRD/\nmCjQWWuGjwFV6yCe1MZanu5rbAAr0RpHsaDVP1E1sDaKgCvHaym00oj1qZUEGUht1A4WCuZ9lqhl\nY805kTdmNOuujUpb0+7NyWMjPRXJAFoHDpIdUz5Pk2QS4kyCDCQstM/7jDknfOgeQvOWRYeSrVtb\nW6OOuRLNbJjHRlE1sIRkhYWIlJA0MJuURZJesZHX9s2+78O3IdTMPDrwV9NT2IyDqBtYFuxwoW9r\nQ7vkR9LAbHKJuax6MaayUWm3Stt3dnamPlcz30mK2WZIjCeJMpHIG1hNSd0MJthENPvWli817zzE\n+E6GC/9yRyEz1LNbX1+fY0lkITGeXMvEepeMqBvYkJCY9FqQVixs+VJzjNiUXZr/pSMtjkz5T1oc\nJJVHmtxp0S5/FmzOC40JNrAkelgs4iZW/4fSOJjyn7Q4SCqP9i1xpS5daRNpsaYVNrANSJIoISQT\nkYGmZYqIfmIaRGOo07633HaFhrgN1fbS4IjZgCQFQUMykey4LETalykKgVr+1jgYudDDl62yPFiI\nsU7HqLMUaHs3sIFNQUhb5mmR0zcsRDKojlcb8VvvJTeNMeBCD9e2Kvk8yTXyyMHaSGLHRb01ARtY\nw7h+2cTGtpShPIXSjm2ba/Jp9Tw/G3mmpVGV6jfbcrnyj5Y4IMQW1TkgNSfYwCrC1hOPJD+3SQ3g\nkLFt83rnl9ogAYxDQK4NJDz1THseybFOCGkMG1ih2FwLtPplIakDYlI4CJmlWTxIs3dra6uV80rT\nUzOulg1Kex7ttY+QmGEDKxSba4GaeFmo1hvzp512Wu7zZkHaIBR642Pa3nmnrNjastFnXIUWQ9Jy\nlJAYCH0aGhtYkolaTfDAwIAHSeQRy3I2puCUlcnErj8hvgmhbruYhubTTmxgiXWqA7z0k28IBSIJ\nbEYISU/I9SFk3UKBdTsZPu3EBjYFWYuOlmJlS87qAC/95MsCQUJDS65rIE198GV3G6vAmDh/Fhi7\nRBvRN7BpkjZrw6WlUdMiZ+hwINELc8gPvuxuextXl3ppil3WSAKwgVWVtLGhpUiZlpMxSYguNKzc\nIUEGU/iee0lkEHUDK2XbQ9dokclUI5dnmaUktmLDSYgMJNY2QEaNsLFmuE9q6SNJPg3Us5cWO0bd\nwLooKhIKVzW+ZbLZrNYizzJLvm1FwkTKACFFjmpszy0NDVd+lG7ftDtMxk49e0n3c4moG1hiB1eL\nlhM5cHBIh8kcyGN7qbloe25paPjwY5r1myX4UWqsa6bk12b+teV/NrAWqLXbVZ5zaaNYLBrTXwMx\n6VoPDg7+CNn2UueWaqzLpkmzfrPEGKUP81Pya7MbTlv+58hrAVO7XZXOpRFT+msgi66SiqckWUg4\n5I0ryfPPtdZlrdioUfSheSpt6mJcYQObEw7+stDiD0nFM8/LEHlekCNhk/fpqdSnr8Q9kuqlBiTk\nhgufsYHNiZbEkhDQLtDiD1P4Xn+y2QtybHDNY3vKSmVM2awbeXM1tlzPSyxjAAlvxYl6sIE1gG9H\nS/6pjdhFgl8bxV+eFSB84zuv62F7ek5lTEmIL1NI9acrQvIlMUOjnNAQL2xgDeDb0b6vHwuxDYBJ\nn/SFGn+h6hU6rl8kIUQr2nMi2gY2tmYkKVrs4kNODclu8udlKS/iJfG1lrjNSh79JNrGpEzV59KQ\np4Q0Ikt+SMxz20TbwJoqcqEFjevizwXLzSKl6TRJEl+HHg959JNomzwysWEloZMlpmPMg2gb2GYk\nbaw0BY3EZluT/SQg0YckO/RnelgzCMDcIWxga1IoFJzsBOM6AVn49VPPhyzmeqj0lcac9BlrjHNS\nQmPuVJJmJzNSGzawNXD1Fq72BPQNk/0DbMeSRls3ktmnPlryvt4SaD7l12K7RtiIPY35mRUXG2S4\nIM1OZqQ2bGArkBLYJWzKE8L2p7ElO5uudDSSOc/mDbGgeQk0ydjIJY35mRXX6wezLshFfxdjkOrA\n9h24NopSoVBAoVAI8mWf0HExSPmOeZ/E1ATkJeY4aQTtEh6sC/bImy9sYBvgI3BtPxktFotMSFIX\n7bHBBsINWeNE2i8/puMlb/6klYfxTjSTN19kVRPFmCoksTwZ1Vh4NcocG9ob8DxoiE9p9U3acopp\n5fEd7xpijoQLG1hD+C4k2tBoL2lTTAipRGpOhfLSTSNs2V667lJjLnSkx4Ur2MA6hoEXDsVi0ehP\noqbONXXqVCPnITpoa2uzcl5fTxVNf18zMetO6qMpLqy+jG7tzIqxaXBNgVcPU/bx1cybvK7Jn0RN\nnevkyZNGzkPy4Sq+R0dHrZw3hFpFCPGLzTrCBrYGLNyNMWUfV3YOZetJPr13Q6Wd89hca5z5pp7N\nY43/WPUmpBlsYA3CQiOHersdaR4c2RC5wdVGJqQ29Wweqy9i1bsartBAqommgXURzDa3+WQyfkAS\nW6QdBKUOEvQ7kYjLuLR5rXq7jRF5aFuhoRGs62aIpoG1tSmAq2unOUfoyRHTagCSi3BaQvaTBpLa\nP88Nog1sXou7jREXhDKNTRqpXl8dHBzED3/4Q+zduxf79u1Df38/vvrVr+LKK6+ccNx3vvMd/PSn\nP530/XPPPRff//73J32+Y8cObNu2DYcOHcLZZ5+NNWvW4Oqrr550XH9/P+666y489dRTGBoawuLF\ni3HDDTdgwYIFadQwRktLi9MCWCgUEgV+o2OSnqPWcUm/6xqJMpHJ0E9+aWb/Un7H6KeWlhZxa9RK\nQGrN10bsNrQVR6ka2L6+Ptx///2YM2cO5s+fj9/85jd179bb29tx8803T/js9NNPn3Tc9u3bcfvt\nt2PZsmVYt24d9uzZgzvuuAMnT57EtddeWz5ufHwct9xyC1555RVcc801mDFjBh555BFs2rQJd999\nN84555w0qhjB9d27yye5tY6LPQkJSUJHRweGh4d9i5Ea7fmdZ5CMqXlNYyftMUFOkcbnNppNW3GU\nqoHt6urCww8/jDPPPBP79+/HDTfcUP/EbW341Kc+1fB8Q0NDuPfee9Hd3Y3NmzcDAFasWIFisYj7\n778fK1euxPTp0wEAPT09eOGFF7B582YsW7YMAHD55Zfj85//PO677z5861vfSqMKsQDv1vMRiv1s\n69Hs/L7tqLF5DYEQcscFtFN8pPG5y/jIW6tTzYFtb2/HmWeemejYYrGI8fFxnDhxou4xu3fvxvHj\nx7Fq1aoJn69evRonT57EM888U/6sp6cHs2bNKjevANDZ2YnLLrsMP//5z62thUiSk+RtfxuEMrcy\nlIHFth7Nzp90ikwocUOIC7Tmi1a5YyDvWGHtJa6hoSGsWLECK1euxKpVq7BlyxYMDg5OOObAgQMA\ngIULF074fMGCBSgUCnj55ZcnHFtrruvChQsxNDSEN954w4IWJCsamzEbhY7FUyaxzvWUgskd7Igb\n8uaLr1qYVW4JtduVDO3t7U6uYxorexB2dXXhmmuuwfnnn4/x8XH88pe/xCOPPIKXX34Zf//3f19e\nuqS3txctLS3o7Oyc8P329nZ0dnbivffeK3/W29uL3/u936t5rdLfP/7xj9tQhwhH8sYKaVePCKWp\n4ksxcZA1Zhkb8aGttkmQ19WL1CMjI5m/63PcstLAbtiwYcK/L7/8csydOxf33nsvenp68Ed/9EcA\nTj2lrdf5t7e3T5hLNjw8XPPYjo6O8rnyElIDQfQRUuyxQYkDXzEbY62WoLMEGWLGlu3z+NVnPDj7\nHeezn/0sCoUCfv3rX5c/mzJlSt3Of3h4uNycNjq21OROmTIlt4y2HSHhJwnJ2P5ZUZP9TcmqSWfX\nuLTNzJkznV1LAtUbBJjeMMDXfHufmB6fsthNQ/Ma2vx2n5swScdZA9vR0YEZM2bg+PHj5c+6urow\nPj6Ovr6+CceOjIzg/fffx1lnnTXh2MopBSV6e3vLf7eJiQKcNkgkJWEjWUzJafupnaYklTwtIhRc\n2ubYsWPOriWB6iUGbS45yBjPRqh2C21+u0Rd2tpq/3jvumdx1sAODAygr69vwpOI+fPnAwD27ds3\n4dj9+/ejWCxi3rx55c/mzZuHl156aZIz9+7di6lTp+Lcc8+1KL2fHVskBW4jWSTJaQM+DSWuiClG\nYtKVkJCot+qT617AeAM7PDyMgYGBSZ//4Ac/AABceuml5c8uuuginHHGGdi+ffuEY7dv346pU6ei\nu7u7/Nny5ctx9OhRPPHEE+XP+vr60NPTg+7u7rp3BCVC+1mBuMHknC8Xye3iSXleOTTiSh8TMaLF\n9qHf+PpESwwQ/fiMtdQvcf3oRz9Cf39/+ef8p59+Gu+++y4AYM2aNTh+/Dg2bNiA//Jf/kv5qei/\n//u/45e//CV+//d/H3/4h39YPldHRwfWr1+PLVu2YPPmzbj00kuxZ88e7Ny5E9dff315EwPgVAP7\n8MMP47bbbsPBgwfLO3EVi0V84QtfaChzoVBAa2sr973OgYlGTuMLACHJG8r+9T7Qsq0yEJ7tQ8fn\nzkd5thbXjhadpMvpU7bUDey2bdvwzjvvADhl2KeeegpPPvkkCoUC/ut//a+YPn06uru78dxzz+En\nP/kJxsfHMXfuXGzYsAHr1q2bdL5Vq1ahra0N27Ztw9NPP405c+bgxhtvxNq1aycc19LSgu985zu4\n66678C//8i8YGhrC4sWL8fWvfx1z585NJLv0QKiHBLldbmNLZOM6HiXEfzXS5CHZkBBbPq8f4pay\njXxa+TfpOpVklS6nTwq7du0K1jovvvgiNm7cCOBUA8xgkIGEQaMSm/Lwyb8OpMUkSUdbW5vx3RgZ\nE8QEjKN03H333Tj//PMTHRvNdijj4+PqgijN3JJ6x0qcCyXNDzblYfNqHxMx7nM9U9KYJDaysZW4\ntDpF6iM5jzSvPiSdaBpYjaQJ/HrHsghng0VED2nm+0mD+dmcUGwkMf5KSJYtCaHECEkHG1gFSCgu\nse1drqEgZo2LkHyZxgYafErCRWL8ldY3lyhbrKTxhckNQiT0GWkJZyRzgC8HSygu3BpUHll/mgrJ\nl65yQ2px5/KAJA82pjgxHt1h0n8S+oy0sIFNgQ8HsxicIkQ7uG4+NBYoKdi0XZ4YCOXF1BDzO1ZC\niMdqKuNTQqy6kkGCro1gAyucEItBFqTZQXpiEz1Ii20fNLIBc434pjI+JeSrKxlqrX8tCTawCZHm\nOHIKH37xuSyKhDiUIENI83gBoL293bcIYpHQMCRFQm5IJKtdaE9ZSMvFsEYBiyTdrSRWYpofbHJr\n2bTnklBATMuQJXZczeN1FdcjIyNOrqMZDfVVQn5KpNHGAo0+z2NPDfFiCk26mpSVDSzMGVRa8eL8\nSqIBybEjWba8aBr0gPRNkE8kytQMHzLbXP4x5NytRpOuJmVlAwtdzk9DqHqFitRBz6Rctc4lVe/Q\nCeWlVInTWrecAAAgAElEQVR1TqJMzdAiM+uFf0wu35UHNrApYfIQW0gdQEzKVetckvSOLb9d6yvJ\n10QnGmPIZp5lOXdeeaTsMMkGFsDpp5+e+FiNyUPqE1vDQhoTW37Hpm8tpC2RRMLD5nsDWc4dSt6z\ngQVw4sSJmp/7uLMhbgklkUtoiD8NMsZOTD6StkRSHlz7jVOC3FBtU+1xago2sA2I+c6G6CzExWIR\nbW1tvsVoSJYcKRQK6OjosCBN3NSzacx1TGPelzC9nm6z75ieEqTZ9jbRlo+u/MgGtgElJyR1BpMv\nLLQVjRKjo6NOruN6lYvh4WFn14sF2lTn060suafhgYwG2+chlh7BlR/ZwDag5ISkzgg9+QiphPFO\nKrE5l9TmwG8qjrlsIWkG/WYWNrCesF3sYrnT80GzxbcJkYjt+LQ5l1TDwJ9HRunTfgiRCBtYT9gu\nyBoKfhI0NYXabK7JtpVwsM+GtvishdaYbYaraT+E2MBXXrKBDYBQizogc9A1/aKELyTaNgkc7OOl\nWCyqyjFiBvrcHiZs6+vFPTawAaCxEQm1IGnyRag+IGGjKcfS4CMftdSAUHwu0d6+bZvn+mxgPSEx\nkCtxOV/ONdJt7wrfhSsEGEvh2CCJHhpeJpN+zZgp2TuUnPENG9gqWlrsm6RQKIgvHNLly0PIuhG3\nMJbCsUESPaqPYSMiF8m+MZkzkvW0DRvYKsbHx61fI0nw2gpKqcEuVS7iBvr/FLRDNnzZrbKWe3uR\nxcFDFx/XyksoN1bNiEXPWuiJxsiwFZRSg12qXMQNEhqBanzIwTw4RVrbS7CbLxlcPHTxcS2bSKkx\nWpBqLzawGZHqUJfQBsQGEpoRQI4cUgltPmglJnRjfbRH3t0xfceXNqTaK8oG1veyEaGgwQbVvtb0\nE1geOHg2hzbKh4b8z4oJ3UK2j2+4OyYBIm1gpQU1B1J7VPva109g1T7WsspDyLGZ1kaSbMEnhHHC\nTTyIRLiRQcRIaqhrBSIHuvxU+1iSzxthe+H4vOd2+US9ls8kvDzk8xyhI632mdjEQ5pOptCql+T6\nmhRftYQNrAW0JhJQOxClDHQm7arZR66x6f+85/b9Uomv3GD8usGFf137Uko9b4atF/mk5Y7k+iod\nNrAW0Bo09RJbyg4xWu1KJhLKz6A+B8KsudDe3m5UDknNQOlpfEmmadOmeZMljV1iqmsS7GL6vJJy\nQCPcSpYYoV5ih7hDTKPzsyDZxcTPoK5pbW2d9JnGxmNkZMTo+WxPMUlD6Wl8yS+Dg4PeZOGT+dpo\nzJlmhKiTS7iVrACkF45QcGFn29v9nX766VbOS+wxNjbmWwSxcADPj6law61KSV40xU50DWxe59R6\nEgOEucWgRB1cDpa2rnXixAkr53WFxLjwgQQ7SJCB5Md0reFNhWwk562m2Imugc3rnKRPYmw3tC6f\nRMZCEptKLjyuMBUXp512mpHzlIjxZRhJP+GbJESdCCkhoXaEQHQNrKZlbxrJygQw78skNqXdzTEw\nMGD0fLH6xrbejfLMVj2N1ZeEkORE18BqKoyaZPWBdvtoXxbM51MyPqFzR6M8056DUijFc/X/x0SM\nOoeED/9F18A2I/Ykil1/l5gc/ENcKULqtdMSQk6FoIMPktqtFM/V/5/nnNrQlNO+kOx7H/5jA1tF\n7ElkeqoDIRLwOQc9hJqiXQdNU8eSnNOlPqzvtfFhF+15aBo2sGQCWZLSRFKxSMpHs4/SxGgSPbUM\nJJp9ZhIt/kpKCKuxaKfejQRzLh3cyIDkojKAfBUrFkn5xOKjkPQMSRdbhNZwhKaPZCRvvS4R0zts\nsoEl0SSclMIuRQ7t0I7EBKHVv9D0yQvrhBxMxyYb2ECQnqQS5JNS2KXIoR3akRB9NBsLWlrMtiVa\n60R7e7vX65v2gw3kSygcKUshSU9S6fI1QkLzTQghIdBsLBgfH3ckiWxGRka8Xj+pH3yOj2xgc6J9\nKSTTaLhrS0szv7DBrQ9t4xap9uaawcQVpvzNuEmGz74lvG7DAwz0D6i+a5s6darxa0izdwg3HraI\ncT3LalzqKTUWuWYw0QbjRj5sYHNSKBQY6A04efKk8XPS3mEQix9j0ZMQCWjNt1hu6E3CBjYnWpOF\nEEKkon0w9yW/drvFjOsNKkKADWwdGEiEEFOkqSesPfofDHA97bgwlbOS/NdIpzz6mqxvbGDr4DKQ\nQnzxiRByiiTTjCRsJkL8EsqNS4zjWYg520inPPqatFV8kVaFhKLhetkQCToTP9TzvfSYaG1t9S1C\nZpIUbFsDoM83sqXHlDRCaYK4DFY6mCfZia6BrQ4Wm0VDamCGUihJeur5XnpMjI2N+RZBJab8WjpP\nmppmM6ba2tqsnZvoQuo4CzSXrdavM5L1MU1eXaNrYF0O1NKbgtBJ81OWxKKRpPiReJDgbyk1bXR0\n1LcIVnHpawlxlQfJ63Q3k63W36XkWFoq7ZzU5nl1ja6BJfGQ5qcsiUUjS/Ej+bH18kJeQve39kbK\nJCZ8nXS6UOhxJX16TihU2tlVTLGBJbnwkcQ2rqmxGOWVWaPOacmio42XF0K1tUm9Qm+kXON6ulCo\nMV4itsZfA2xgSS58JHGWazYrrhqLUVKZ6+kuRWebA58UHaXIYZokeoXe2EhAgo1DiXHp9bIREuKg\nhAtZ2MAqQFJQakVD8bGF5Dli2nayYy6mR7p/Q/glQ7qNNeHClrZiJu9DDZO4sCMbWAXEuFICcQf3\nqU+ORHlDzGGXy73l9amNX4RCIRY90+K7jvi+vimibmCZXOEEMqDDnxpkJLoIKYdL2Jq/KSX/NPtM\nylJqWpAScyESdQPL5EqOhiSU+oINd1kittCQl5LwmX/afVWSnzUsHaHZq1AoiInlqBtYkpzQkrAS\n27ppsJ3UheGlFEqpaIittJjwucS40e4r7fL7QGIc5qVYLIqJBTawJHqkFBmfTaTUheF9FUqpDX0M\nmPC5pJdZiAx8+Drvr4KMz8awgSXRI+VuUmoTqQWTxZ6+iAMpuU/SkzbfNfm6JKsGmX022WxgCbGI\nqeTmnXhzNBR7F4T6EzyJl1rxGEu+S89Fn35gA0ty4Su5tOzGZSq5YynWvpE+WCTB5U/wRAea41rb\nWtGmiVn3ZrCBJbnwlVyl6/raylLKgCBFjlAIZbBwGReMQXdktbXmuNYgu/QckC5fVtjACkBacEmT\npxG+G2jfSJGDyMJlXDAGm2OqpkqztYuxQsN4JH1LZWlxYwo2sGQSoQZ7PTQUyCyEqhch2gi1prrQ\nKxTbhaKHJNjACsBkYMfQtJjWMdTCEqpehFRSWQ9cbkpCiE/a29t9i5AImznDBjYwYmhaYtCxFr7W\nJuWgTRrR2tpq7dxJYq+yHnBTklMwZ8NnZGTEyXXyxpLNnGED65k8TQmLVFz4WpvU9qDNONbN2NiY\ntXNraRilQbsRU0iOJTawBskyEOdpSnwHFhsPYgLfcawFafk2bdo03yKoRtNub5pkLSEtXzSgzWZs\nYA1iayCWGlTN9DUld0tLeGFq2qcu5wFqJo1tpNmxOt98yzc4OOj1+trRtNtbFll9xydvjNNjymau\nfB9eZxAgJoJKwz7Q9Zqw8fFxYzJJwXRxdTkPUDNpbrqk21G6fMQcruu3iesxPmXhMoZc+Z4NrHKS\nBqXUYlKvYZAqrw2yFBbfTzdCJaa4s4mm+NTQHLqOS+ZBPiT+amjTp1lzKG/uybMySYX2QuNKfskD\nahYbaPc7mYzkGE2L6fi0aZuYmsOQYswmee3k+1dD137OGtN5c4ENLMmNhqLIhq8+GvwXA41iNHYf\n+W76QrE/62AytL/PEoufo21gmz3iD6VgZSGt7i6SJWZ/5CHtOppEJvV8pCEvfP2c2sw2aaZfMUfS\noyE2bVKtf2trK+MIZutBtA1ss0f8oQZakqcJEhvSUP1hG+12axQnsQ+QgA7/+vo5tZltNNhOM6Ze\n4s3yfdukeTBQOjbLeskh1jiT9SDaBjZGqguEz+Tg4BEmLrf5lRZDLp80ah3YtMptEqk28C1X3pd4\na33Hlk5p5MtTp6TVuDS4iCc2sJ5xtl5aSwtaWlrQ1taG1tbW8nV9F608aJY9VCQWXFdx4vJJo0Q7\nJ0Gr3CZJaoPSDVFM8yZd3gBrptpOEsdCF7ZnAwu/znexTWfpf2eccQY+9rGPoaurCx0dHc4LpGlC\nLU4kH9XxzDghNrBdN0s3RNriN49dtOmal6xTpKrtlNduaX0mpWdgAws5SWMyKKrP1dLSgg9/+MO4\n8sorceGFF2Lq1KnlY6ToHxtSikBoMJ7NoyVW88gZ+rx8V1NctNnFJ1KmSKW9lhQfR9vASizISYMi\n65vlx48fx+DgIM4880y0tramlo+YRUoRIOaRWF/yoCVWY51vmATfa5PmIbR8SoomvX3IGm0D64q8\nTq31/TSNbunY8fFxHDt2DEeOHMHg4CCGh4eDL9gx4/slPR9I0Ne1DG1tbU6v5xIJ/kxKR0eHbxGM\nUigURO0mFetYpUlvH7KGW/2a4MrYzRYnt7nUS+V3x8fHMTg4iJ/97GcYGxvD4OCg6jvyPCSxu3Zs\nLsRt8twmzyfBpyUZXMkyOjrq5Do+kODPpAwPD/sWwSghrX0ruWbFMBbZJNoG1hR5AtDGdou1zlks\nFlEoFDA2Nob3339/wrVjTB4XL86FalfTeoVqJ9ukibGQ45GYIeQYybIerau5qaHa3BVyfiNQiqQA\nTJJ0lXfWNmWPeY6tpJggZkjzc7aLn75drUMZAjHXoqRIjxGX00mk24J8QKonsIODg/jhD3+IvXv3\nYt++fejv78dXv/pVXHnllZOOPXjwILZu3Yrf/va3aG9vx9KlS3HjjTeis7Nz0rE7duzAtm3bcOjQ\nIZx99tlYs2YNrr766knH9ff346677sJTTz2FoaEhLF68GDfccAMWLFiQRo1ocZmYWXYdsUWepwua\nnkwkkTWtPpr0twkbRv9kjUVJtSgpzLuJhG4L+jsbqZ7A9vX14f7778cbb7yB+fPnA6h9Z3T48GHc\ndNNNePvtt7FhwwasW7cOzz77LP7yL/9y0pyt7du342//9m/x8Y9/HH/xF3+BCy64AHfccQcefPDB\nCceNj4/jlltuwc9+9jOsWbMGGzduxNGjR7Fp0ya89dZbafUmSjBx5635zeQ0+ieR1dRyKZpesNEG\nbVubyliU9IKRDXzXHeKW2Pxt6iXjVE9gu7q68PDDD+PMM8/E/v37ccMNN9Q87oEHHsDQ0BDuuece\nzJ49GwCwaNEi3HzzzXj00Udx1VVXAQCGhoZw7733oru7G5s3bwYArFixAsViEffffz9WrlyJ6dOn\nAwB6enrwwgsvYPPmzVi2bBkA4PLLL8fnP/953HffffjWt76VyQCxUu+Or7W1VdQTi9ASO+2dtlT9\nY3ii7Qva5xSNYiXWF1BjQkKtkCCDaSToZOr6qW5j29vbceaZZzY97sknn0R3d3e5eQWAiy++GHPn\nzsXjjz9e/mz37t04fvw4Vq1aNeH7q1evxsmTJ/HMM8+UP+vp6cGsWbPKzSsAdHZ24rLLLsPPf/7z\noN/GtUG9AJLUvIaI78Lhm0r9Ne3+ou2paAi20pQroT8R9oEE/0uQwTQh6WQ86w4fPoxjx45h4cKF\nk/62aNEivPTSS+V/HzhwAAAmHbtgwQIUCgW8/PLLE46tNdd14cKFGBoawhtvvNFUNm2DUKzY8tN5\n551n5bwkGzZ2f7EVO6WVPLTgc5AKaYBsRGU88IkwIe4x3sAeOXIEADBr1qxJf5s1axaOHz9eflra\n29uLlpaWSS92tbe3o7OzE++99175s97eXnR1dU06Z+mz3t7eprJpLqyaBs+82PLTK6+8YuW8MSI1\nHqtjx6ScmusHMQ/jITtS6wfRhfEGdmhoCEDtnUlKn5WOGRoaQnt7e83ztLe3T1gcenh4uOax1efU\nTr3EZrEkvsizG5xvtD05JcQFWaY88GbQHKxJZjDewE6ZMgVA7Z1JSp+VjpkyZQpGRkZqnmd4eHhC\nE1zv2Opzaif2xCby0B6T2uUnejDdmNhqdKqnPCTZkrjZrpIhkEQP3yvjkA8w3sCWpg6UphJUcuTI\nEcyYMaOcLF1dXRgfH0dfX9+E40ZGRvD+++/jrLPOKn/W1dU1YUpBidLUgVrTC4gsQilyxC++48jG\n9X3rZJKQdEmLxJ3qkvgj70vQoTRkNpYiTEPMuZMF4w3s7NmzMXPmTOzbt2/S3/bt21dePxZA+b+r\nj92/fz+KxSLmzZtX/mzevHl46aWXJgXP3r17MXXqVJx77rkm1TAK31A9BdcU9UsodvY9WLraAlor\nIekSAj78EUqtcQ1zJx1WOqtPfvKT+MUvfoHDhw+XP3vuuefw5ptvYvny5eXPLrroIpxxxhnYvn37\nhO9v374dU6dORXd3d/mz5cuX4+jRo3jiiSfKn/X19aGnpwfd3d2JfgLxhck3VGsVBknFIossTFo3\nhGTnkG4KQ/ILIUCYMR1SzQmF1F3fj370I/T395d/zn/66afx7rvvAgDWrFmD008/HZ/73OfQ09OD\nTZs2Ye3atRgYGMBDDz2E8847D5/5zGfK5+ro6MD69euxZcsWbN68GZdeein27NmDnTt34vrrry9v\nYgCcamAffvhh3HbbbTh48CBmzJiBRx55BMViEV/4whfy2kENtQqD62LR6IlRiIUrRLQ/9eOyRXGh\nPV6Jfh+y5kympaXFq10Ku3btShVR1157Ld55551TX/7d07bSm77/+I//iDlz5gAAXnvtNdx55514\n/vnn0dHRgaVLl+JLX/oSZs6cOemcO3bswLZt2/D2229jzpw5WL16NdauXTvpuP7+ftx111146qmn\nMDQ0hMWLF+OLX/wizj///Jqyvvjii9i4cWMa9YgltBcvQkg6Qs35UPXSiG1fuPJ1aDGVR5+77767\nbk836TppG1hNsIElhEhE44ClUWYSPozLsEjTwHJSByGEOEbjgKtRZpIdSe9WNCKUuNRib0mwgbVM\nKEvuSEsuafKkQbPsRAcaY0yqzFLlso2GxjCrbyS+kOXL3q2trdav0dLSYiWP5HkxMGwEpY9Ar7ym\nhAbad3HNYwPfshOzSGxwNMaYVJmlyqUZUzmT1TdpXzySmOOmGBsbs36N8fFxK3nEBrYOaQI2T3Bn\nXf6rUCh4SyrfDXReuJMKMQljQR4hNxwhoC1ntMkrASdPdq1fQSlpAjZPcGfdAaVYLDKpMuLTbhxY\nCbEPayMhfnHxZDeqBtZX86ChadEgYwhwYCVENqyF9jBpW/qJRNXA+moeSteVnHBsrEgjJMcuISZJ\nUgtN5IOpnGp2Hp/TzaoxOc5wzGpMDDcLUTWwvmHChYPUhLZFabOSJMRmG+0knauW1q8hx4GJWm5q\nPGh2ntJ0s5D9QSaTJb7qrc4gtXdhA0uc4qOI2phMniahpQ8cSe2TVGepxY7UJulctbR+zRoH0vNF\nIxIX+6ef5aFtu1w2sBEQe6FwMZm8EbUGDkk+8W0fQiqR1miFgESb+l6akeiHDaxjfCSq7ychleeR\nWEh9IMkOHDwICQttOe27Hrqylza/SCeaBlZK4NS765T0UkAJV3O0iF/oH0LCgjl9iqRjoit70S9m\niaaBlRg4lTJJeimAEO1IuWElhPiDY6JMTNXnaBpYQkg8cOCSDW8wTuHbDjau3+icvvWNCVfLaGW5\njqn6HGUDW2upCK2JpVXuakLRg/iDMWQWm/asN4C52sJbCr5vtGxcv9E5fesbE67W3PXp0ygb2FpL\nRZhygu2iWn3+RnJrKvAsbBMxsfSXJv8D+eVlDJnFhz1dbeFNCNFPlA2sTWwXVVcFXlvzExomlrbS\nNsBrk7eEhlxpdEOkQX5CmjF9+vRUx+eJ+3oL/udBex4mkd+0jmxgf4et4NEalFqbiSxo9VFSQtfP\nN1Jypa2tre7fGt0QSZFfE8wpefT396c6Pk/c21jw32Ue2pgfm0R+0zqygf0dtoKHg4N8Qt9m0cR8\nQyKf0dFR3yJkpt4TLakxyrpOXGEjB2zMj/WRq2xgDSO14CZFk/wmZY1xQHKlc14/SYtJn/JIs4Up\n6j3RijEvQ8PV2/DaqadbdQ74skGz6/rIVTawhtFecJPKXygUygHtK6FsJnbIhTILeeyRNyek5ZRP\neaTZInRYB/Lj6m147STdWteXDSTang1sFT6bIB+ToPNQamKlBDYLpT18zM+SFOuhIM2m0uSpxvc2\n3CROOP4kgw1sFT6bIB+ToBsh8ScDIhsTA3cprmzGV6x7n0vLWWnymCJUvbJSLw+k5Yc0eUhj2MCS\nTMRcoFnk6qMlLmLd+5yxS3xQLw+k5Yc0eUhj2MAKx/SAUzl3tdk1miVzsViMLuHzTJmI9akfSU6l\n76S/fQzoizVt8kqk0XJtpqCfSBLYwHqgMjmbLYhscsApNV/V58xyDS0/CZkm6UT7Zt+1SUg3FaHH\nUzWVvtPgR2kyJpn21OyY2GKuFo1s4GK5NmlxRWQSTQMrqShVJqeNBZGTXNfWuWIqPCGvHytFr7Tx\nZGOHHKIHE+8RxFTD6mHDBqHlppQaGTNhRVQDWJQaw2ScTJK34U3FlTT728oX23qOj4+Ls2VsFAoF\ntLS00A9kAi4f1riguka2t7d7kiQfmvM0mgbWBJod3Qw2+JPJ8za81OkFNmltbW16jAs9Q7ClZkrz\n7Bs1sTZq6bRp04yfk5CkjIyM+BYhE5rrJRvYCmq94FSJZkeHhIYbiRhjZWxszPk1NcRCTCR9ydFG\nfgwODho/JzEPc9YMLuwo3VdsYCtI8la9KYfaDIyk507yxEwiMTaHpDa1YkF60ZWArdyv/NVCep7a\niJPQ5nkmIcZfmyTAX7PYwKbGlENtBkbSc/t4YhYDbKCy08h2Se3aLP5NLAPkYikhm9jM/ZaWFoyP\nj6duYl3nTdoanES+0OZ5JkF6k+MTlzEd47jDBpYQmGmcSmgo6L6LXZa1h5PatZluJpYBcrGUkCYq\nbV5qjm3sROgTV/L5zs0YSWpzyWOB9PyxARtYoUyfPt23CHUJscCaaJw04VsnCb9A2CTEHGmEBJtL\nx9QvCMQ8tLlO2MAKpb+/37cIdZGS7LE1CYBZnWO0nyuk5Ege0sQHY6k5pZiQaCuJMkkkhLwOCTaw\nyom58JgoJtrsxwJKXJHmhVbfcakpj33bqhY2ZNLkE1fQJmaJsoE1EURSAjFr4ZEiv2+S2M/1Czuu\nfBPqJgwkPVl8KKkRkySLRHy8TESfTIY2MUtUDWxpiRMTQaQ9EDXIL2WZL9cv7GjwTSWm5WVD7B5t\nMUc+IEm+8GUiYhsfdTuqBjbGJU40E8oyX3nXhoytoTM5ANqwXQz+qNQxBn01w4aRSMBHHEbVwBLi\ng7w3ThygsmPDdjH4o1LHGPQNCd4w26G9vd23CKLhE1hCCCGiibnBkah7tUy8YTZHpW1HRkY8SmIO\nWzHMJ7CEEO/UK3AatlGOGVf+ibnBkai7RJlCIUTb5tFJWu1mA0sS4zN4XVxbyktjvqlX4EytYxni\noCABDdtc18JVXcn703olPmuhST1qIa1JIbXx4SdptZsNLJiwSfEZvC6uHcpLY7aRVsRq0Wi+GvNd\nFq7iyeRLvD5zoKRHSD8FS0dizXDhJ4l6V8IGNgAaBZn0AEyDb118Xz8GTNm40Xy1NIXfx/qZhCSB\njaY7tNr6tNNOy/V96XqzgUVtJyUdTNIOOjYGqXryl65V+m/tA6TvZPJ9/RiQNj/L5Rq3jC8SGtrH\nHO0MDAxk+p4Wv7ndYkgRSQeTtIOOi0FqGoDFxSIKhcKp+VKFAorj4ygWiygC2AcgW1gTkp1CoWA1\n/jU0gLZktG1bUh/avj6x2kV7TNiU3aRt2MAKIq1j6x2/GMBzAFAsAjXmdV4M4NeZpUxOM32S6qu9\nGJBT+PRh6DEUsm6NkOBX39evhQS7xAxtXx+TtmEDKwjTT3OvmzIFM7u7ce655+LZZ5/FaQcP4t6h\noTwiJiZJAbX1lDuLLCQ8Kv0euv9bW1ujfAkxdL9mxbZdWFOJBDgHNmD2Aug5fhxHPvYxHF+wAAdy\n7CSSdk6MpOImSRbijnp+1zK/Kw0xNq+1CHl5K0mwpurBRU74yrt4Ms4zaRxsKhhGx8bw5ptv4tCh\nQ3j//fcxOjqa+RosWEQalS8ppoGxHC6ufFsr5kwu0+WSUF7yJbVxkRO+aiqnEDgijYPzBkPp2+Nj\nYzh69CgefvhhDA8PY/HvlhYqGrgGIb6JZXqAdGL8OVmzvtX+0qyLLXzEtOtrhpC3bGADZlGxCIyM\noHDsGIrFIs5XGqyaEk2TrHmIRc+kxGyPWPV2henYor+a48NGrt/89xUHXIWAJOIB4NRKBMoKVktL\ni9qf42IZHGLRMym0xwfE3MznRVKjQcJBUgxxFQLSkH0ALmryd8lUN6+Sko/ohs2VffJuRmFi6T2t\nhKybdipjL/Q41AIbWMuYDPSkS+UMAtht5IpxwuIULvSrbEwtvZcG5jtJAucN2yFP/nEVAsuYDPRQ\nl8qR9vYrixMh7vG13E+M+S6t5tomNn2rcam/y1Vh2MCSSZgO9mbn0zCAaCiAGmQkpB7N6oCJ+DZd\na6TmXAg11ySx6VtNWv3zrGns0tbRNrBSC48ETAeg1OKRJgak6lCJdBmZcxOxYY+QbSwxviXKBLiV\nK+SYi5UkL1G3trYmOlfpOBtxEm0DG8udOKmP1MEnVGLaGSsJNuKPMU1KuMorxlx8FAqFxFMaS8fZ\niJNoG1jTxHrHm2UHF0nyS8G3TXxenwNgbXzHBNFN7HlVK3+YU2aQEltsYBUiJXhK2Ji0rbHQtLVl\nX9TDt099X59Mhj4hWTBROzXW32pq5Q9zSi5Z5t2ygSVNaVTMKtfFM4nGQjM6Ourt2iEMOEQGocdS\n6A2eidpZeQ7JukqHtktOls2L2MCSpiRZm1FjwxkStH84+B70ssaSb7mTYrrBC52YdDUNbWcXNrBC\n0GJPNj4AACAASURBVFL8a8EkJcQcJvLJRz3RUgdCfwJL0kFfnkKjHdjACkFL8Sf50FgkpJN0OZeY\nYD2pTxbbVOetVvtKb9554+UPjXZgA+uQEJsXjTrxjfuwCHWHujT4WFNWUu7bliWUOaHSm/dYV/OR\nTBo7ubYpG1iSC1svcdkkxiZSk39IeipjOs8uOvXOmeXvLnEpiyS9XRCqvqHqZZo0dnJtUzawDgk5\nYULWrRbaGsLY/BMzSd7m1Ra/JDmmbmBMwVhzQ4x2lhXphOSEu8+QNJiOFy2DiK/41WIfzWRZjsgm\nEmplDHEnwc6N4FayhDRBehKTfEhfb5jx1xjaRz55NmSRiu+444um3EqWREoMd88kGb4HIuKOennP\nemCXWhuy0Ob5cPWiaWx+YgNLvLJixYqmx7BpIY2IdRqAZEzYsF7e+6gHscdEUptrtpNm2UtIHStt\n2ZYNbARITswdO3b4FoEoh9MA5BGaDbXo47vWu7CTLR21+FgjtmzLBjYlvgtEFpiYH6DRfyGTxB/0\nmTtoa91obiCTwvGMlGADW4d6SVoredIsW8IB2y+hLEgeCkkGIw5Y7shqa+ZSPDSKEcZBvJR87zIG\n2MDWIU0hT7NsCQdsOWiyc2VR4CDhFtp7MtOmTZvwb0255JPQYymmOGBNnkjJ9y5jgA0sIZBfjCqL\nQkyDhGZCXjpncHDQtwjeyeJf5q4/bC7BR79OxsU4ygaWELAYkXTThpLgaumcPEi8WdOCBv+SD7BZ\n15lHk3ExjrKBzYn0wJUsn2TZiHykb2rgg7Q2MalziPkcok5poQ2aw53t/MAGtoIsL1glCVyfQSZ5\nUJYsmwZ8TJqXRL34idUegN+cCjGfQ9QJSJcjodogBGL3TVQNbLM5S7ZesAo1yGJuFCTgY9J8FlzH\niXR7EOIb5ggJgagaWM5ZMguLIEmC6zjxdWMV0w0ddZWBZNmIeUL1d1a9ompgCSHhk7dhzlpMY7qh\no672SbLKQUx+MIH2BjBUf2fViw0sUUtbW5uV86bZmIKER6iDhA+yNAzamwxTZPnFkLZrjKvcph/c\nwJGaqGV0dNTKedNsTGEbaYVQkjySZNGGK9vxnQG3lGzH3PBLaDEsNZ6sPML6zW9+g6985Ss1/7Z1\n61YsXry4/O+DBw9i69at+O1vf4v29nYsXboUN954Izo7Oyd9d8eOHdi2bRsOHTqEs88+G2vWrMHV\nV19tQ4UJFAoFLwHp67ouCFm3JCTV34SNTNpaks8kyaKNRis4SLdrLRk1yO0SKbbw6RfpMSFdvkqk\nymnnN9jfsXbtWixcuHDCZx/+8IfL/3348GHcdNNNOOOMM7BhwwYMDAxg27ZtePXVV/G9731vwk/E\n27dvx+23345ly5Zh3bp12LNnD+644w6cPHkS1157bWKZWlpaUj9h8+U8qUFjgpB1S4JL/WO3tTZ8\nDmwaYqWWjBrkNkGz2PARO21tbXV/DdOwrJuvfDN9TU0NsSmsNrAXXnghli1bVvfvDzzwAIaGhnDP\nPfdg9uzZAIBFixbh5ptvxqOPPoqrrroKADA0NIR7770X3d3d2Lx5MwBgxYoVKBaLuP/++7Fy5UpM\nnz49kUyNmteQAiCrLlkafCKLkOI4LSHorln+EOwvmWa29WF7W1O5GhHqr0p5CEWPNFidA1ssFjEw\nMFB3MvqTTz6J7u7ucvMKABdffDHmzp2Lxx9/vPzZ7t27cfz4caxatWrC91evXo2TJ0/imWeeyS1r\naIU3qy5sXmsjdQ5QLUKK47TErLsEKu2vKWeIH0yt+MFY04Upf1l9AnvbbbdhcHAQLS0t+E//6T9h\n48aN5SkFhw8fxrFjxyZNMQBOPYV99tlny/8+cOAAAEw6dsGCBSgUCnj55ZdxxRVXJJarVrPKgY80\ngvFhjmnTpmFwcNC3GCQFWW7wpf3UTeRhYsm7YrHIWFKGKX9ZaWDb29uxbNkyfOITn0BnZydee+01\nPPTQQ7jpppvw3e9+F/Pnz8eRI0cAALNmzZr0/VmzZuH48eMYHR1FW1sbent70dLSMunFrvb2dnR2\nduK9995LJV9a47HY2oO2jQ9pzStjsDmmm1HaOx+aY5Y//8eD7Ti10sAuWbIES5YsKf+7u7sby5Yt\nw/XXX4977rkHt956K4aGhgAAHR0dk75f+mxoaAhtbW0YGhpCe3t7zWu1t7djeHjYghYfYOouURIm\nZcrzYkEaGSTaMQsubU+aI8l+ef3pIx4k2S8WNNtcs+wkHbZ97Wwd2HPOOQd/8Ad/gN27d6NYLGLK\nlCkAULP5LH1WOmbKlCkYGRmped7h4eGaTbAkJCasSZlqnatyjgvvtj/AdIOR5k1bkp40djNh47yx\nEUKOuII5EQf0c7g43chg9uzZGB0dxcmTJ8tTB0pTCSo5cuQIZsyYUV5Gq6urC+Pj4+jr65tw3MjI\nCN5//32cddZZ9oUnqeBAWptGdrG1s1iz6+alcstLk4OFhIEnjd2a/QohDYkyuYQ1Kh1a40WLn03Y\nV6uPsuK0gX377bcxZcoUTJs2DbNnz8bMmTOxb9++Scft27cP8+fPL/+79N/Vx+7fvx/FYhHz5s2z\nK7hBsgaY5MCULFsapk2b5vX6PpajMUHlKiM+ni5rQKIuEmUicgk9XvKOY3m/b8K+ofuoGisN7LFj\nxyZ9duDAATz99NO45JJLyp998pOfxC9+8QscPny4/Nlzzz2HN998E8uXLy9/dtFFF+GMM87A9u3b\nJ5xz+/btmDp1Krq7uy1oYYesASY5MKtlq3wipwlpLxdJo6WFO08TQhrjetqNKTh9Rx9WfrP8H//j\nf2DKlClYsmQJZs6ciYMHD+LHP/4xpk2bhg0bNpSP+9znPoeenh5s2rQJa9euxcDAAB566CGcd955\n+MxnPlM+rqOjA+vXr8eWLVuwefNmXHrppdizZw927tyJ66+/PvEmBibgSzPNqbfub1a4uYIMNPmA\neUp8E2sMlnROon9W+9h82dGW32KNB5tYaWD/8A//EDt37sQ//dM/YWBgADNnzsSyZcvwp3/6pxO2\nkp09ezZuv/123Hnnnfjf//t/o6OjA5/4xCfwpS99adJ8wFWrVqGtrQ3btm3D008/jTlz5uDGG2/E\n2rVrbahQFwagezQ1TqGirfhqkjUPvv3i+/ppcXkzrMkuNrCpv8mnpdUxbEvu2OOhGhO1o7Br165g\nrfriiy9i48aNmb+vrTgTf4QYK7518n19X7jSO1b7xgh9TbRw99134/zzz090LCe1NcBkwkua62Ob\nemv2hkyIg4NvnXxfPwsSltKqpJE8Gu3bjJjqbBqk+brSTxJ9JlEmre+W2IQNbBW2AldaAbFJvTV7\niR0kFttYkZbn0uSxTWz6asXFT/Z5kCiT6XdLQoAN7O8oNQEmd0iKnWY2KBQKtFNKatlLYrElhJxC\nU43TJGsttMtP0sEG9ndUNwES1nSzSRL90tqgepmlJG+gSreTNNJuvUuILUKIrxB0MIn2eqxdfpIO\nNrB1CD0RkuiX1gZcLUAW2naGIroIoUbmfgs6QR6FYCeTsPYQU7CBJcQy0gq21jeSpdmREI15VMJX\nPkm1mZaNWrL4TWLtNCGTDo9ZRKJj8yL9Dc/YkFawpcmTFFdLS7mGOaoXzb7TWgdsoeUXxCx+k+hr\nEzJF38CGCH86JiQ5lTlhoqimzTGJgwupTbVv6TsiiSQvTodE9A1s6AWoWr/Q9SV2CK3wVWI6J5hj\n4dLMtyHnSWiE4KssL07nRZLdom1gszhByxwZQkzDpkwvXADdHdrzRFJzYpvq7WQ14mLag+RfHaLt\nyLI4QcscGUKIXTQNeDEugJ7GP5p8aZuk42JoNpPSlLW1tfkWYRJSbFOLaBtYILwk9AFtKAe+vOeO\nrEVdil/qySFBPtfb8UoeoLUhIX40Mzo66lsEVUTdwBp5Cy7yaQVSij8Lp/ztGdMixacm5ZDil3py\nSJBPggyAnPhzRSN9k/6MbMt3Ln0Rm9/rocEOcXdfBkg7rUBDUGjEZOGMxUfS9ZTSyEiRg7glNr83\n0te3LWxeX/IcT58ksYPvMYQNLD5wggtnmEoO34HjAls6NjtvLAUsBj1jyBNCSHZiqIO28G27qBvY\n0uBWcoL0Ceym16uUji0d05zXha/ZZNkj1DwJPWZC188ktBWJNQaibmCzDm4uB8XQmlZtiebC5iH4\nNXZcx7WpmJG6xBZzIjkx2Urb+OGKmGKgkqgbWA2EFpih6eMCjUVbo8x50BbXJf/YXmIrTxzYjqEk\n548tjqXDjSRIJWxga8AkIKbJE1PamiNAp8yVhF4DXPknz3Vsvy+Q5Pza4zg2fPor9JphAtM2YgNb\nAxat+tQKwJATN83SMo1gTOmC/goH+pK4IMQ4Mz22m7YRG1iSiloBGGLiJiGENQ+JfejP7NB2cZPF\n/4wZc9Qa4yTZlw2sIyQ5PSmxb9IAZGtS8/raRmOsMf6ayewyPmObAuKLENfkPP300zN/N/YanGXF\nGJv1U2MdNY2knIw7OxwiyelJSbtJAzmFRF9LlKkZxWKx4YDhMj412k8jIdm5FLsnTpyo+7dmaK/B\nLhs+m7FTb6lNNrR+YQNLCBFLSA0NiQvJO1u5woWePpvIWPwoFTawxBu8eyWVMB5kQr/Eg0Zf21yt\nQqM9YoINLPEG714/wFSh1FxwGQ8yoV/iQZuvTda7WF5QljhGZJWJDSzxguskkpi0lZgqlKYLrnS7\nEdm0t7f7FkE9UnOwra3NtwhBNpi2qbSZlNjK6kc2sDVwvXSHlCAKGS55lQ3ufCMDrW+jj4yM+BZB\nPVKbtNHRUd8iiERTTTQZWz701lkVLZPFqWm/U+lsqQXKJqHoHIoeWYldfxvUGgi0v41OSCzEWhN9\n6M0G1hOxBrkUNN0lk7hgbSBS65NUuchkYvAVG1gSJWwSSIzEMKhpo9InNhfjN4FUuchkYvAVG1gS\nFBygiWR8xmehUEg0qDGH3FLpkxiaDhIevmoGG9iAcBVEkgc4EwNAPf0k601OId1HPhuUZjubVR5H\nksEtRgmQ3P+hxomvmsEGNiBcBZHLYM2b8Fm+X08/DuzyoY8mUh3/MdrHZtNQb4tREhd5bg5DbWpd\nEG0Dq3VZmtjIOzD4GFiqCxILVHpoMzPkjf8Q/JDVBq2trYYlIY2Q+jQ7qTxZ40zSzY+2vkiXtAbR\ntixNoySSlvCxU12QJBUoLdBm9UmS76aOCd0PjWwwNjbmUBIi9Wm2NHlsoq0viraB1UajJJKQYNru\n3AjRSpJ8N3UMISFiYhcxnw+ONI+3Ju2m1wo54VNLs2i7cyPEB7brDutaOtI08dpsq01el5jYRczn\nDaDm8dak3aJpYOu9zOAyyVlQSCxofkLQiLxTeWwPenyqmpy09VjiHvKN0BoLGmybF1M6xmCrRoQ5\nytTA9JvlWQZorQWlHqWfYZolkc0ka2triz6JJaL5CUEjQsvhmKn2ZZo6wjiwRwy2NaWjy5sqieNs\nNA2saUIdoNMwOjqaaHF0mwVpdHQ0ioJHZJN0kwAil9j9J7FBIcmR8OuO61/e2MCSXGgt+tVLtsRa\nvKu3sYzVDnnRmgfELLXyR8si94zh9OTxmW9/28D1gz02sMQI2uY8Vi/ZEmvxrt7G0oYdQizUhNSi\nVv4kzamki+ETOeTZmjnWMcckuroOIhZOqSD1qFeoYx6sY9bdJdrsLKGp0WYzQLbMEnwaKtE2sJID\nnvjF50tpsSB9zmgWH/MloIlIyBMpdpZgi6S4tpkJ20jxM3FLtA0sA57Uw8RLaZoGLB+Yzj/T9s4i\nH2vKRGiPD6At6hO7bbSPFV43dPB2ZUICJuSiLLHghmxv10j0LyGhYqJ2+cxZn7WXDWxCWNTDgb7M\nB5vFsKF/SUyEMB7EmrNsYBMiLUBCSDpfSPOlC/LGi414YwzbgXb9gCS2sGUv+kEHjcYD+lA2bGCV\nUUqoGJsw7VSvueqSvPFiI95C3ofeJ6wNH5DEFrbsFZofYszBpDtdxWgbCbCBVYaGoshkrk31mqtp\nidmuGuKexEkseRl7DjbSX7JtWltbfYtgDTawxDiSk1kzIds1yy5gsTQORDYh5yXRz9jYmG8RatLR\n0ZH7HGxgCSF1cdUkZtkFLE/jEPsWwlmgrQjJh6Ycsr0F8vDwcKbvVcIGlhBSF+nrtWbF1BbCrvUJ\nfbkcKfERErSpHEzcdLsizRbIvmADa5CsAcYCQ2IhabHTMm/LdfEuFotB1wv+HJ+denERu0015EsS\nGWP3Yy3YwBok6w5NvgJTUmJLkcX2FqJ5ziXFRi4wPW8rJNuF+lSc5IMNTm002EWDjBJhA+sYSYGa\nVhabT5il2MX3FqJa33SVDm1XH9rGDSZuFHizIQv6wy9sYFMSc8BmHeg4QBJCKpFYR23JVJoOY6IO\nspbKgv7wCxvYlPgIWInFPi82dGppYTgT/4SYr7WopWdS3SUO/LZk8rWMUSxxSOwi+d0ejviGseE0\nicU+LzZ0Gh8fN37OEGBj75ZYNqmopWeItUor9AUxgeRfXjmyGYZFQzcaG4lmsLGXD+sGkcJ//s//\n2bcIqbnwwgt9i6CKUB5qhKGFYGw0RCE2WWmxZQM2EqQaLUt6acJW/rI25mf37t0T/p3Epr7t/vzz\nz3u9vjZCeajBBtYyNhoiH02WtEGcjWZ8+HpqIHUrxkp8NxBpsZW/zc6rzU5ZMK1jEl+xHseB72Um\nq2EDa4AYiqKGQTwNXG9VH42eGsTuMzYQyYjBTtrnYEuSRRq+beN7mclq2MAaIIaimBffiVeNqfVW\npekVK7Zy0JR/GSdh0NbWBiA8f5rMn7y24XhaH9pmImxgiRNC3UedBSVsTPmXcVIfTc3g6OgoAPqz\nEbSNPzTlkgnYwCohtsDMAgsnIW4xUZeYtx+Q154aXriKGdu2jy2X2MAmIM+C3SauBXwQmDEUnxh0\nDA3fPvNxfd86SyC2AdMkpekIleS1p40XrkKL85I+vn+x821XX9c3eV02sAmwvWB3pUObnTeGASMG\nHaWTtsj49pmP6/vW2RQ+B1Lfg7hPStMRpBNKnJco6VMsFr3Gny+7lnRudH2bdjGpNxtYAUguEDEP\nMDEjOSZdoSn288jq09eMM91oypFaxBh/IS2LxgaWNERLIBNiGk2x7/OnSU7fiBdNORILMeUGG9jA\niSmY00C7kFBx3VT4nr7BXDaLzzmiJD8x3VSwgfWEq+KgfVFrW0hK8hjsTfwTapxJyuUQkDJHlOjD\ndbywgfVELEW31pu2SYmleMYSC8QvjDOSlqQxE0utzkJMtnFdY9jA/o7qIKvXeLkIRt8BXx2EteRp\nbW1NdK48b9pywCWx4DvniTtaWsIbdqXVakn5JM02IRFeJmWkOshqNV6FQsFJMEoL+FryjI2NeZCE\nkDCRlvPEHuPj475FCB7X+SSpYdZIVvuxgU2BzaSodiATIhxcb3pBCPEPV2eIh9huQE3HWVb7sYEV\nQrUDazmUxUkWSf1hq7iFXDQZ60Q7vldnIB8Qw9JyLihNf5ESZ2xgFZEnaEJNKJ/YTuKYfea7QIZq\ne5fbYjdDk41dyqrJLlqIYWk5F0ib/sIGViGlApfmZQBNCcUCfgpNPguNUG1ve1vsNGiysUtZbV2L\ndZWEBhtYQaT9SVra3VAaGumqaWDTiOaBTLPsxD8xx0/MdbXS7zZiIOa48gkbWEFo310mjcwxF1Pf\naLa9ZtkBnXkN6JW7GonxE4ptJVPpdxsxkPac9LkZ2MAKRWKhbYYJmZnYJGQ05jWgV24NSLQt63B+\nTP7KmNYfsfiPDSwRhban0BpkrIUtuTs6OtTapIR2+YleSrHnOwZNN9UdHR1Gz5eFPDbN8l2TNky7\nI5rEmyIbsIEViO/iJQVJSVjPJ5JkTIrNDTmGh4frnltLXGv0aezkjS3XzU09SrEXWgwODw/7FiGX\nTbX4Q4ucpmADK5DYglADIfmEb51/gKumWkvzrpW8sRVDc0PcUy/vWQ/M0OZbgDQMDw/jvvvuw2OP\nPYb+/n6cd955uO6663DxxRenPperbWHTIlUuQkLERa4xpwmJk3p5z3pgBlVPYG+99Vb88z//M664\n4gp8+ctfRktLC2655RY8//zzqc8lNYCkykV03TVrkrUZ2nVhTpOsaI99k8Rsi5h1b4SaBnbv3r3Y\ntWsXNmzYgI0bN2LFihX4u7/7O8yZMwd33323b/FIBGhqRDTJ2oyQdGmEhkGqUCiokDMUpG0p3tnZ\n6e3asdSBWsSseyPUNLA9PT1obW3FypUry591dHTgj//4j/HCCy/g8OHDxq8pqVBLkiUEpNmztbXV\ntwiZkWZLrdgcpEz5qFgsqhxMQ4pRn/bv6+vzdm2TPpSy2oNr0ujb7FgJtlPTwB44cABz587FtGnT\nJny+cOHC8t9NI6lQS5KlhIQAzoo0e46NjfkWITPSbOkTqTmhoTm2CWM0OVL9aWNZKldx4cOm1ddM\nOxe/2bESNm9Q08D29vaiq6tr0uelz3p7e12LFD2mk19q4SQkKTE2SjHq3AzNtUyKPzXbsBofNq2+\npm+/2ri+mgZ2eHgY7e3tkz4vLZA8NDTkWiTShLQFyHeCEfm0tEwsWT7X/9RMrHq7grUsG5VxSRuS\nZqhpYDs6OjAyMjLp89ICyVOmTKn7N+IHFiBimvHx8Qn/9rn+p2Zi1ZvIhnFJ0vRtahrYrq4uvPfe\ne5M+L00dqDW94NChQ9blIoQQQggh+UnTt6nZyGD+/Pn4zW9+g4GBAZx22mnlz/fu3Vv+ezWXXHIJ\nvvnNb+JDH/qQiL2YCSGEEELIRIaHh3Ho0CFccsklib+jpoFdvnw5tm3bhh//+MdYt24dgFMKP/ro\no7jgggswe/bsSd+ZOXMmPvWpT7kWlRBCCCGEpOA//sf/mOp4NQ3s4sWLsXz5ctxzzz04evQoPvzh\nD+MnP/kJ3n33XXzta1/zLR4hhBBCCHFEYdeuXWpmTQ8PD+Mf/uEfsHPnThw/fhzz5s3D+vXrUz1y\nJoQQQgghulHVwBJCCCGEEKJmFQJCCCGEEEIARXNg0zA8PIz77rsPjz32GPr7+3Heeefhuuuuw8UX\nX+xbtCD4zW9+g6985Ss1/7Z161YsXry4/O+DBw9i69at+O1vf4v29nYsXboUN954Izo7Oyd9d8eO\nHdi2bRsOHTqEs88+G2vWrMHVV19tTQ+tDA4O4oc//CH27t2Lffv2ob+/H1/96ldx5ZVXTjrWhv37\n+/tx11134amnnsLQ0BAWL16MG264AQsWLLCir3SS+uM73/kOfvrTn076/rnnnovvf//7kz6nP7Kz\nb98+/OQnP8Hu3bvxzjvvoLOzE4sXL8Z1112HuXPnTjiWOWKfpP5gjrjj1Vdfxfe//3289NJLOHLk\nCDo6OnDuuedi1apVuOKKKyYcKzVHgmxgb731VjzxxBP4kz/5E8ydOxf/+q//iltuuQV/93d/hwsv\nvNC3eMGwdu1aLFy4cMJnH/7wh8v/ffjwYdx0000444wzsGHDBgwMDGDbtm149dVX8b3vfQ9tbR+E\n3/bt23H77bdj2bJlWLduHfbs2YM77rgDJ0+exLXXXutMJw309fXh/vvvx5w5c8rLy9XaWcmG/cfH\nx3HLLbfglVdewTXXXIMZM2bgkUcewaZNm3D33XfjnHPOcWIDSST1BwC0t7fj5ptvnvDZ6aefPuk4\n+iMfDz74IF544QUsX74c5513Ho4cOYIf/ehH+PM//3Ns3boVH//4xwEwR1yR1B8Ac8QV7777LgYH\nB/HpT38aXV1dGBoaQk9PD/7n//yfeOedd/C5z30OgOwcCa6B3bt3L3bt2oUvfvGL5eW2rrjiCqxf\nvx533303vvvd73qWMBwuvPBCLFu2rO7fH3jgAQwNDeGee+4pL3O2aNEi3HzzzXj00Udx1VVXATi1\nDfC9996L7u5ubN68GQCwYsUKFItF3H///Vi5ciWmT59uXR8tdHV14eGHH8aZZ56J/fv344Ybbqh5\nnA379/T04IUXXsDmzZvLvr/88svx+c9/Hvfddx++9a1vWdZeHkn9AQBtbW1Nl/ajP/Kzbt06LFq0\nCK2treXPLr/8cqxfvx4PPvggvvGNbwBgjrgiqT8A5ogrli5diqVLl074bPXq1di4cSN+/OMflxtY\nyTkS3BzYnp4etLa2YuXKleXPOjo68Md//Md44YUXcPjwYY/ShUWxWMTAwADGxsZq/v3JJ59Ed3f3\nhDV6L774YsydOxePP/54+bPdu3fj+PHjWLVq1YTvr169GidPnsQzzzxjRX6ttLe348wzz2x6nA37\n9/T0YNasWRNuXDo7O3HZZZfh5z//OUZHR3NoppOk/gBO5cz4+DhOnDhR9xj6Iz9LliyZ0CwBwDnn\nnIOPfvSjeP3118ufMUfckNQfAHPEJy0tLZg9e/YEX0nOkeAa2AMHDmDu3LmYNm3ahM9LP3UfOHDA\nh1hBctttt+Gqq67Cpz/9aXzlK1/B/v37y387fPgwjh07NmmKAXDq7u2ll14q/7vkk+pjFyxYgEKh\ngJdfftmSBuFiy/4HDhyoOUdp4cKFGBoawhtvvGFKhSAZGhrCihUrsHLlSqxatQpbtmzB4ODghGPo\nDzsUi0UcPXq0PG+POeKXan+UYI645eTJk+jr68Nbb72Ff/qnf8K///u/45prrgEgP0eCm0LQ29uL\nrq6uSZ+XPuvt7XUtUnC0t7dj2bJl+MQnPoHOzk689tpreOihh3DTTTfhu9/9LubPn48jR44AAGbN\nmjXp+7NmzcLx48cxOjqKtrY29Pb2oqWlZVIha29vR2dnJ9577z0neoWELfv39vbi937v9yadszK/\nKuezkQ/o6urCNddcg/PPPx/j4+P45S9/iUceeQQvv/wy/v7v/7781IP+sMPOnTvR29uL6667DgBz\nxDfV/gCYIz6488478eMf/xgA0Nraii9/+cvlX7Cl50hwDezw8DDa29snfd7R0QHg1N0dyceSnYR2\neAAABh5JREFUJUuwZMmS8r+7u7uxbNkyXH/99bjnnntw6623lu1csnsllb5oa2vD0NBQTZ8Bp4J/\neHjYghZhY8v+zK/sbNiwYcK/L7/8csydOxf33nsvenp68Ed/9EcAQH9Y4PXXX8eWLVuwZMkSfPrT\nnwbAHPFJLX8AzBEf/Mmf/Akuu+wy9Pb2YufOndiyZQs6Ojpw5ZVXis+R4KYQdHR0YGRkZNLnJeNN\nmTLFtUhRcM455+AP/uAPsHv3bhSLxbKdazWf1b6YMmVKTZ+Vjq2VPKQxtuxf71jmVzY++9nPolAo\n4Ne//nX5M/rDLEeOHMHXv/51TJ8+Hd/+9rfLK0QwR/xQzx/1YI7Y5SMf+QguuugiXHHFFbj11ltx\n0UUXYevWrRgeHhafI8E1sF1dXTV/ci5NHag1vYCYYfbs2RgdHcXJkyfLPzmUfoKo5MiRI5gxY0Z5\n+Y2uri6Mj4+jr69vwnEjIyN4//33cdZZZ9kXPjBs2Z/5ZZaOjg7MmDEDx48fL39Gf5ijv78fX/va\n13DixAncdtttE34KZY64p5E/6sEcccuyZctw4sQJvP766+JzJLgGdv78+XjzzTcxMDAw4fO9e/eW\n/07s8Pbbb2PKlCmYNm0aZs+ejZkzZ2Lfvn2Tjtu3b98EP5T+u/rY/fv3o1gsYt68eXYFDxBb9p83\nbx5eeuklFIsTd6Deu3cvpk6dinPPPdekGsEzMDCAvr4+zJw5s/wZ/WGG4eFhfPOb38Rbb72Fv/mb\nv8FHPvKRCX9njrilmT/qwRxxS+nn+0KhID5Hgmtgly9fjvHx8fKkZOBU4jz66KO44IILJiwFQbJx\n7NixSZ8dOHAATz/9NC655JLyZ5/85Cfxi1/8YsLSZc899xzefPNNLF++vPzZRRddhDPOOAPbt2+f\ncM7t27dj6tSp6O7utqBF+Niw//Lly3H06FE88cQT5c/6+vrQ09OD7u7uCYtakw8YHh6edFMNAD/4\nwQ8AAJdeemn5M/ojP2NjY/j2t7+NvXv3YvPmzbjgggtqHscccUMSfzBH3FJrHB8dHcVPf/pTzJgx\nAx/72McAyM6R1j/7sz/bnOhIJcyePRuvvfYa/s//+T8YHBzE22+/jTvvvBOvv/46vvGNb2DOnDm+\nRVTPN77xDTz++OM4dOgQXn/9dTz22GPYsmULpkyZgr/6q78qv4U4f/58/N//+3+xa9eu8hymrVu3\nYu7cudi0aRNaWk7dP7W2tuK0007DP//zP+PVV1/FiRMn8PDDD2Pnzp34sz/7M24BXIMf/ehH+NWv\nfoX/9//+H1588UW0tLTgzTffxJ49ezB//nx0dHRYsf9HP/pR/OpXv8KOHTswNjaGV199FVu2bMGJ\nEyfwV3/1V5gxY4Yvk3ilmT/6+vrw3/7bf8O7776Lt956C3v37sUPfvAD/PSnP8Xv//7vY/369eVz\n0R/5ufPOO/HYY49h6dKl+A//4T/glVdemfC/0tMg5ogbkvjj8OHDzBGH/PVf/zV27NiBQ4cO4a23\n3sKzzz6L22+/HQcPHsR//+//vbzMleQcKezatavY/DBdDA8P4x/+4R+wc+dOHD9+HPPmzcP69esn\nPB0k2fmXf/kX7Ny5E2+99RYGBgYwc+ZMXHTRRfjTP/3TCVvJAsBrr72GO++8E88//zw6OjqwdOlS\nfOlLX5rwc1CJ0h7Kb7/9NubMmYPVq1dj7dq1rtRSxbXXXot33nkHAMovQRSLRRQKBfzjP/5j+UbN\nhv1r7WH9xS9+Eeeff75FjWXTzB+nn346/tf/+l/Yu3cv3nvvPYyPj2Pu3Ln41Kc+hXXr1k1a5B2g\nP/KwadMm7NmzZ9JPlMAp//zbv/1b+d/MEfsk8Ud/fz9zxCE/+9nP8K//+q945ZVX8P777+O0007D\n4sWL8dnPfnbSQyOpORJkA0sIIYQQQsIluDmwhBBCCCEkbNjAEkIIIYQQVbCBJYQQQgghqmADSwgh\nhBBCVMEGlhBCCCGEqIINLCGEEEIIUQUbWEIIIYQQogo2sIQQQgghRBVsYAkhhBBCiCrYwBJCCCGE\nEFWwgSWEEEIIIapgA0sIIYQQQlTBBpYQQgghhKji/wP1Hx/M3gBRDwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10c21bd90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(data, origin='lower', cmap='Greys_r')\n",
    "cutout.plot_on_original(color='red')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}