scottdraves/imgdiff.ipynb

## imgdiff.ipynb
{
 "cells": [
  {
   "metadata": {
    "trusted": true
   },
   "cell_type": "code",
   "source": "import matplotlib.pyplot as plt\nimport matplotlib.image as mpimg\nimport numpy as np\n%matplotlib inline",
   "execution_count": 1,
   "outputs": []
  },
  {
   "metadata": {
    "trusted": true
   },
   "cell_type": "code",
   "source": "img=mpimg.imread('/Users/spot/src/beakerx/test/resources/img/groovy/bigInts/diffcell2_case1.png')",
   "execution_count": 2,
   "outputs": []
  },
  {
   "metadata": {
    "trusted": true
   },
   "cell_type": "code",
   "source": "lum_img = img[:,:,0]",
   "execution_count": 3,
   "outputs": []
  },
  {
   "metadata": {
    "trusted": true
   },
   "cell_type": "code",
   "source": "plt.imshow(lum_img)",
   "execution_count": 4,
   "outputs": [
    {
     "data": {
      "text/plain": "<matplotlib.image.AxesImage at 0x10d59ee48>"
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAB2CAYAAADY3GjsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJztnXl0HVl95z+3qt6iJ+lpt6zVkizJbtndbbdXmWVCBxLoZIBOOoTOOQmQDA0ZOCGZhKFZJsscSJgEQoYJCWkIgWSgG8IyEBrohqZZbXd7ba9tSbZkW5K1L0+ypCdV1Z0/qt7T2/WeNsvu+znHx0/16lb96v5+91e37q13v0JKiUKhUChuf7RbbYBCoVAoVgeV0BUKheIOQSV0hUKhuENQCV2hUCjuEFRCVygUijsEldAVCoXiDmFFCV0I8VohxCUhRJcQ4tHVMkqhUCgUuSOW+x66EEIHOoDXAL3AMeBhKeWF1TNPoVAoFNmykh76fqBLSnlFSjkPPAG8YXXMUigUCkWuGCsoWwNcj/m7FziQqYBX+KSf/BWcUqFQKF56TDE+IqWsWGq/lSR0kWJb0viNEOIR4BEAPwEOiF9cwSkVCoXipccP5FevZrPfSoZceoG6mL9rgf7EnaSUj0kp90op93rwreB0CoVCocjEShL6MaBFCNEohPACbwa+tTpmKRQKhSJXlj3kIqU0hRDvBp4CdOBzUsrzq2aZQqFQKHJiJWPoSCm/A3xnlWxRKBQKxQpQvxRVKBSKO4QNmdCNpgZuPnQA61X35V62oR6javMaWJWMHgyi3bN98dy1Ncz82gH0HdvW5fyK2w+tsBCjpnrJ/YRhINvvRRjOQ7QWCBB+YN9am6dYRbTCQvSKiqz8vVqsaMhlrTCv9FA4MIQ9M3OrTcmIFQphjBViu3+bvX0EevuwbqlVio2MPTWF0JfuR0nTRL8ZxracaLJnZvB959ham6dYReypKZieRi8Krts5N2RCFz4f9swMenkZ1ugYSIleXIQ1MYkWCKCVFGPeGATbQgsEEIaBKC9FhqaSjqNvqgDbRhYEYGgUvB6swaG4/cK/sg/fd45DtssgCAFSogUCmL19Ke2X4TB6eRmUFDnbpmcwbwxgNG5BGjp29zWkaaa007rUtYxaU6w1WmGhk4yFhjUZAttCDwbB54PiQoRlY3ZfRW/ditXZDfbird1oakDqGlbnlWgci0Ae1shoynPpJSXYF7qSY9KNPaO2BjwGzIXBMDCv96K3NCEWTOTYBFYo5Jy3rjYu/q3x8TWrH0UKpEzyd6qYQcqM+YKO7E63IRN6lLISGB0DYGFnI2a+gW9kFjkdxqipAtOEPD/mpiB6aI6Fmnq8150GoldUIOfmWKgtw3NtGOkxmDvQjDZv45tfiDZI++W7CHSMYuWypo2UWK+6DzG9QKg5H++UjTTA9GsUPnEULRDACofp/61tVD4/jRnw4Ds1zsg72ik/NY3t1dELtsHpC+jl5Ul2CsOIJnvFxmDid9opOR/C9uiIU5fQS4uxRsewt21BHwkhPQbM3oR9OxH9Y4jd29H7RzFvDKAHg4zv30zxM5eZ+bUDTFfpFHUv4Pvu8ZTn0ttama8sAAlz5R6EBQsBQdHjx0A6z4MT7bXkDS8gNYHv1BVG395O2dmb2IaGZ3YO3esBoSXFv+8Hp1RsrQN6RQVTr2gi/8nTDP6XPY6/v3cSvXUrdlEgKWZ4/mzGfJFtQl/24lzLIShKZTa/FI32cLc1Y3Vcdnopbs8kgv3yXXj7J5B5Pqzzl6LbjYZ6WDCxS4LM1RbiH7iJ1AULRX68h89jzy8g2+8mXOKl4MQ1Fpo2I35+OveLEcLtVUus4WHA6cHZU1PRp4m4xCwERnUVZl8/wjDQmrbAwDCyvjrJTl/fpOqlbzC0QMAZAhQCvbAQe2sd+vAEof215HdPIywLM+hHWBJx5AVnHqVvgJmXtWJ5NfL/40Q0FiLxkQlhGGiN9Vhd3SCl0yYWTCehx7SD2KdBa2QUze/HnptD39bsJAxDi4t/saVWxdZ6IgR6URBrYhKjtgazt4+ZXzuQMmYy5YvvLzxxQkq5d6nTbeweeiwpbjxyfBL72lSKnUEbnSAwOY15vRfh82GEw9Gxbk//OOLnV5m7fw+en7yQvF5BlghNA4+++LffD1Op7RG6DrZjgTRNhGlBYQEihZ1inSZ1Fdkj5+ejn63pm+jzJvbYOIU/nMaevonUdbRw2GmUgHWhA7F3J74nj8H+uxeH1wwDkZcHSyR0dB2Z5wWhgbQQXq+T0BOIbtedOLTn5pz/g3novcOYNwbi4t8ITa9KfSiWRhgG0rKi/rbHxkEICn/4YsqYiSubkC+yZUO+5ZItoqgQrbgo5XdWdRl2mTMZIcPhxTKGgZ2fx/SbDuI7ew1237UutqYjlZ3mjYFbaZIiS7TKCkRJMdI0o76LND5jSx16/yihhw8idQ09Jk6lmcO0uTvEks1EaiJWdZlziJj4V7G1/kT8rVU6a2uli5nV4PbpoSciBNLrYWF7HYh69Ol5zCIfRs8IAPrAOHJ+gdk37Kfw7BBWST4j9wUpPz2NCM1S8NVjSK8HbaIQWpqwOq/kbIIsCICmwf67MUamkFM3czuAbae0U3Rcc2bIFRsWKQRWzzWmfvMARfl5aDNz3LyrAi1s43/hmjPc4fUQfPwoRtVmpl++Df+Tx9AKC5E1FeibSrEudsVNnMaiFRZA2EQeuhdt3kK+2JObgW5sTSXEf8UT51RsrRNaYSHWxMSivy9dRrt7G+NtxUkx4/nBidU556oc5RbgGZpC9t7ACM3hebEXDA1h2k6in5l1xp7y88i/HMIeGEILzVL5vevovcMw4U6Izs0hr/czt6XEGaPPBSkR0zMwEcIYDiH93uhYei6kstO+ubFf11SANjEFUlJ0bgIxOIocG8c/MIP/xBVkTQVibh77irNAnjk4TP5z3eh3tTiT8R09zNYWouUH0h5fhudBSjxXBrB8OnJ2Nr0xaW4KZl9/Uvyr2Fo/rMkQYldb1N8iLw86epJj5vjqzWlsyEnRXNECAezZ2bSvHS71/YrO7fc7j0+r8Ni0lnYq1o5cYyAywaoXF2G31Md9J8IL2GcvOR0Gw0B4vSv+PYaKq1tLdEI9dluOMfMD+dU7bFI0A0sF/Fr+QCkyCbUqx9rgP6RSpCbXGIj42ZqYRAsvoIUW/S7d1w2R1qp1FFRc3VpS1f9q5o1Y7oiErlDcrthnXoy+faVQrJTbdgxdoVAoFPEsmdCFEHVCiGeFEBeFEOeFEO9xt/+5EKJPCHHa/ffA2purUCgUinRkM+RiAn8spTwphCgETgghvu9+9wkp5cfWzjyFQqFQZMuSCV1KeQO44X6eEkJcBGrW2jCFQqFQ5EZOY+hCiAZgN/Ccu+ndQogzQojPCSFKVtk2hUKhUORA1gldCFEAfA34QyllCPhHYCuwC6cH//E05R4RQhwXQhxfIJxqF4VCoVCsAlkldCGEByeZf1FK+XUAKeWglNKSUtrAZ4D9qcpKKR+TUu6VUu714FstuxUKhUKRQDZvuQjgn4GLUsq/jdleFbPbg8C51TdPoVAoFNmSzVsuLwN+GzgrhIgsHP4B4GEhxC5AAj3AO9bEQoVCoVBkRTZvufwMSLVy1XdW3xyFQqFQLBf1S9Fc0PSl91ll9B3bMn5vpBHD0AoL18Kc2wYtkH4lw6ULL+1nvaIi9Re5rtoZOeVG89cyr2MtSFvXiiRecmu5GHW1XHpPLfm9GrObJY2PHkmSt0uH0EREbwAAbVcbHW8NEujVsPxQ9+HDjiTYvLP0ac9H2rF8kqavzSKOvADA9Q8dYr5tBmva4xokueuvxphrLMXz9KLGpNG4hau/WUPNRw/H2XDlo+2Ypc6CTQUVN6n/49QLL4lAXpx6Us+X78E4XcBslUXLHzwX3X7zoQMM7dUwpgUlHTYFXzmass7M673x2xq3cOG9leRfM5ipsml5z9G4erz8sYN4x53+QsNXBqLrza92Ob11Kzf+2mBi0EmIYk6j4VsmeZcGc16UKnKdiX4G6Pz7A0ivu1GTbHvX2fgdDt4DR89EF9ZKRLtnO/aZF9OeO9FffV/fgX2yiLkqi5Z3Of7q/Nf7kAsaWAIkVNaNU/bILPONm9B+eiq1PS43HzrA4F4N76Qg2GMTfNzx88yDBxh68yzzk84LC2JWp+ZHksA3jyddR7prSIyPbGI84udYH0fIVNddf3cQ34hGw78PJUnpdX7hPnzdfhbyJVvfeyTuu5sPHWBsu54yxjs+u5eCF73MVUia3nfEuaGnWZI4W3rffyiprSViNNRz8c8rCFz0MbfJpvkDp6KiF+NvaWf4lQtgCSprx+F12Z33jlg+N2uEIPy6vfifdoRy9bZWxGzYUd1Osz9S0vv+Q9T+1eEkR9sv34X2M2daQW9pQg4MI6or4cZQVHUdWDJAzPv34BuYwrqwqAQ78IeH2PzJ55LKRXRL0/GV3iO8qbYdvXIT1uAQwjAYfOd+qp5yGo4eDGKFQgifD9qamWopXAxwV8dQmib2xCQyHCb8dAN9J6toejS+gcy+cT+Bb5906nHHNphfWEy+tTXMtlVFG6++rRl5tRd7bm5NyiUSevhgNGFlhRD0PtrOntefY7A9FOevSH0ZNdWYff1pD/Hwi/08vr06pa9nn2pk9onNlP7LkaRyif6K1G3hiX7M672OQnzt5rjYiKX/vYeo+eSJOFWigfccwnzlJLW/fj66bf61+/B+7xjgdETs0xdSLusKEH5gH77vHMt8DTFtI1pvGYiN8Vg/69uasS51IXw+NJ8vY10bDfUsVJU4nSMh0JsbsS5fBdtC7N6BNj0bF+PgLlN719bFOI+JcWtwiPAD+wgc7nBWvvT7nTXLvZ6oL3Il0t42/f3huLYW6x9wnjo6/nszzR86iQyHMWqqkQF/NKb11q2OnrJLtsvnvqSGXISuEzjWE12S1LrYmT6Zu/R+4BBVR1MvdekZWkysVlc39tQUVmc3eD3Il+1a3NG2MGrT/7jWytPiGuzII+3UfvVqfGJwH4FFdSV6SxNGTXXScTr/90He/MDb4rZpjfWUnZ+LBkok0MOvuoe56gBFT19c3FlKZHEhhMPopSV0fHo/+b9nJiVzgIKfdy/W44WOuOS6sKUC748Xe1VWx2WnoaxRuURKTo6k/S4dVUfnUiYlaZro5WXMN25C37HNSbAJdH3iII9vr0av3JT03cg72tE+WZ4ymafyF0D+0+eiPV4rFEqbzPWSEipOhZOSBRCXzAHyroxFP2uDzud0TzD+Z87E/Z3uGiJtY6lkDvExHuvnSNLSfL4l61p6DDw3xt0/HIEZoQn0tla0mXBSjBtNDdhzc/FxHhPjAAiwpx2lMXtuDmt8fNnJHBbbW6wd4Vfdk7TfwrYaWv5tfFF2MuEGJkwnZxhNDbmdfxk237ZI04TimLFK9+lEKyxMPW4qJbV/eRj92ZNJXyXeQZESPRhEa2vBGhllqt4f/ar3/YewqkpT2qRvaybQORa3LW/MRgbzQQiE4YyK6WVO+ZmtJVhlBfQ+1MDo29ujdo+9rZ3qn0ro6Ik7ltV5Bf3Zk06P3KXjs3vxPXsG35PHktTn5zflY01MYt4YoPWdzycNtUQpjdFyja1HAFvGCRpreXmL63qvdrlYv7k3vZxV7aVM6WNwk978AgtFHuYr8un84A5m37j4k4uuvzvI9r9NU0fATJVwhKITSOmvyE27sCBp/0gcxDL6n7fj/Wn828Idn95P1aeej9umV1QgZheTvjXmJMXYmIjuW1yUdINIeQ0Z2kbSMRNjPMbPmnvDFvmBjHVt1NUiFkzMnmuLJgQLkJaFdbGT6/85/mba+fk9yJDT4UqM80iMG3W1+J48Fr/mvKajt25d8prSEWlvESJtLRHtZ6fjhq+0XW3ODcmNgYXKIszaMkbbN9P/J4eyPv9LagxdLy5C+ryrc7DJKexX7MYzGGKuvhjv5DxM3ES+6CST4ONH6f9GG0V5cxR/zkYeSxhzdR/NZZ4X+/QFgOijWWGXE4hGdRV4HBdFAjnyKLz5qDMhau5tg+fPYgYEpf/yHPidG4mciZcsE14vY1+v51dqz8O9x+NVxt3H57lf3U/B6T6ykVSQeenrURoaWsw4tD0zg97SBFNTa1Ju8QASo7YGs7cviyvIHisUiia0ph+B2Hc3mqtDu+U7ZvxNL2bwvecj7TT8j9RjqKn8pRUUYE9NITJI00XQi4soPT2OHZN8r/7FIdr+4gpmoiiGFj/BKefnnWNUVcYlSICp+7cT+PqizZmuYUlSxDjE+zn6lKA5fct0dW12XkFva2X2jfspOD8Cw6NInxG9uVd/7DAHX1jga5d38cq6yxi/MYw1NZ1kUmyMW0PDaDu3Ywc8aPMWNxsKKOiYiLv5LZfxJ1tSt7XEKrpnO43/3MOZ0TnyX0v0eiJzbkVHoQi4kP4Qcby0xtA13ZnwcgNe8/ux5+bixi+XKp80ph0IIBpqYXQi7hjCMOLu/EbjFjAtp/G7vUq9tBhrZDS6j15eFvd3xD5t53bsc/GTUbFj6dar7qP/ZX6kDpZP0viBI3R94iDVP5EEvpG6MQrDQFoWwutFLy/DHh3LSUUl9voidsaOfSbWaeTa1qocOGOsWHb6p4psyTDnET13bQ1Dr6lnul4gbCg/YxLon+XKrxdQ9/Q8Nw750EywfCAsaPzCNYbvr2MhXxAuAyRp/aWXlGDPzMT1lCM3Kq2wEK0gHzQteu3C56P/3XuwPc75tHloeKKPcEMZvq4h5NQU0rLj5l60QABRszl+6Gr/3YgXOpwx3bpaut9Sn/YaNv3jkcwvEqSJ8QixfjaaGpBj40lPi7F1nXiT1ouL4vbXdm7HvtAZ57f5X96L75nT0TgXQqSMceHzoeX5oaIsOia/miTmAuHzoQUCWOPjiztpOprXk7YNZjuG/tJK6DiPbtg28uZN7Nk5tC21MB5KK/BsVG1GBguwLnVFXyG0zl8CnAkzKsthcMRpgK7T9GAQNlc4AT84gh0OI+7ain3mkjOB4/OhFTtDCHE3EiHQfD5EUz1ichpsm4WGSoyLPeDxMv6arZT+5DoyPI8wnAYjFxZS3wTchK8Hg6DrCL8PezLk9HorKpjfUYevZwS7IMDY7hLKnr4MZcXYPg/iUvdi71jTYGwSa3gYfVszIjyP2XMteoOSs7PYU9NozQ2IscnoObSd29EmprBHRhF11XFvq6xmOfPGwOKE2vTNZTdGfVszDI1ijY/H+Vlva8UsCWB09YMtoSSI9HuQl7qTJ7oqN2ENDeek3ZnoLy0QQCspBk1Dzs5ijU2gV5QhAnnR+R69rRX78tWUY+dp2X+3I2Y+EQJdg5KiaPLSy8sQeXnYE5MZJ9wTibSN2HqLtI10MR6ZiI34Web5nBcTevuyqmujtgY5PR2XzPVgELwehN8PHgOrfwAtGITiQqyubowtddgFAcTwWFKMg/tapGmCtJNuKrkSaW8yHI5ra8aR8+ibKsBjYF7pcSZnG+qR4xNgS6xQyIm7vgH63rKD2v/oR45POr7SdZ4a+AeV0NOh3bMduvvQykuxR8YyBrEwDPS6Gqz+AfTqzUivxxk7l3KxJ7VgJiUSvaQE4fchS4sQE1PIcNhJvJqOXlaa9gYiDANtawPYzqO7HBh27HN7jUZtDRg6djCA6OmPf5uGxbcyEnu65v178J26wvyuRrx9k841pOuJukMwejCIKCvBHhlD6BqipBg5PRO1XezegXb5OqK8FDkyFmeLVliIEALZVAsdPXETcKtdLpIMU71NkC1G4xbk9Ax2KJTkZ2EYaAX5yNoqZGc3WFZ6rc8cX3lL5y/tnu3QdQ37nmb0Cz3OPhUVyMqypKe1JYnETl0tSIk9GULoWnzy0nSEx8ip/iJtI1JvWkOdU2dCyxjjEONn08SeC0frLFNd6+VlWE3V8PzZpOPpwSCz7a3knehm4a56PC/2Yo+POz3hVDd6N8aNmuroGy+riV5SEtfWoh2syFO4poO00e7ehjY8gV1R7Iypu77SdrUhrt5gbt9W/Mcu89TYZ1RCT8dqqakvRWSyL5deD7g/itH1lOWEz+fYnuMxhWE4vbAcy63kmFphIXJ+PilJrFW5tSSTT9bsnEu8oprz8W7BNWQi3fWtxM5Uw1W3gmx9l+5aE2NdDbkoFArFHYJ6D12hUCheYqiErlAoFHcIKqErFArFHYJK6AqFQnGHkNUvRYUQPcAUYAGmlHKvEKIU+DLQgCNw8SYp5Xi6YygUCoVibcmlh/4qKeWumJnWR4FnpJQtwDPu3wqFQqG4RaxkyOUNwBfcz18A3rhycxQKhUKxXLJN6BJ4WghxQgjxiLutUkp5A8D9P3ntUEAI8YgQ4rgQ4vgCt/Zlf4VCobiTyXa1xZdJKfuFEJuA7wshsv7tsZTyMeAxcH5YtAwbFQqFQpEFWfXQpZT97v9DwDeA/cCgEKIKwP1/dRdDUCgUCkVOLJnQhRD5QojCyGfgl4BzwLeAt7i7vQX45loZqVAoFIqlyaaHXgn8TAjxAvA88KSU8nvAR4HXCCE6gde4f9/ZZKEGv+qn3Lk94/dG1ebU5Taaivx6slLF+hX4eal6v20U7JdThyut93XiTm4bS46hSymvAPem2D4K3LYrbSUukJ8NiWrwYvcOOn6vgMB1HcsP9X9x2FnCdX4epKT7o+1YPsnWr84hfu6ISV/7s0PYbdPMTzsKPMJjs/3Dk8w2lUZFfMFZzrXn4Rpq//JwnA2X/6Ydq3QBgPLKEJv+YD61rQkq8te/uhNxMshszaKKPMD0bxxgaK+GMSMouWRT+ESyuHKiqnvEvgvvqyRwzWC22qLl3c9FlyQFR2fTN+b0FxqeGIhb13w55SLoJSXxwgAJyPZ7Ec+dS1Ksz4bIdcb6OVt/TfxOO8asTcG/xwuKdPzDfvA5BxO6pPUdjmTc9G8cYHC/hndCUNRtE/xSfL2n8tflL+7GXtCQlgAJdXWjFP7ePOHmSvQfObJnQ9/cztzJUuaqTFrfuShFN/vG/fS/TMc7IfBOQeUnnbi6+dABxh+eZmbClYKb1an5IQT+33GQ1pL+SlfvYs8OLr09n7w+AzNf0vjokbjljSN+jvVxxrp+6jhIScen9+OZ0Gn4j9lom0rlw0wkto24a3CVglZK36OHmA864iXpMBrqufg/ywlc8DO72ablfaei9TP2tnZG/9M80hLU1Y464yJZ8JJbbdGoreG3nznM59/8APLU+cw7u4F7/YOHqPvI4aS1ru2X70L7mRNUeutW5I0hRHUl3BiKX6d8iTWyF169B/+N6ag4AMCNPz5E1SeeS1ZIWmJZzqf6T/PL1bui62sLn48b79xDzXcHsTouxymRs7OZydZCgo8fjV5vRBHdnph0guuZWi6frGPre+MDc+bBA+T/xwmkaaLt3I4Izy8m7doaZnZUOY0QRzxCXu3FnpvLuZzdc91R0Gmo57M/+RJvfej3U66HHeGRjis8tr0lN6ELIbj+gXZ+4Y0nubxvLmd/GVWbGXptY5yIcqSe0ynYhx/YF5UTFLt3IE+dRwsEkJaV1l+ploTte/QQtZ84gQyHmXnwAMHjfZjXex2hhboqrItdyAM7Gd8eiLPPqNrsCGik6NTM/ep+/E8eAykz+iuWxHoPPXyQoq87CSoSi9E6ifGzvq0Z61JX2uuL1DXXbzht6uA9cPRMnI+iZdPEaoTEthHLwB8dgv80zuY3XkxZNlv0YJDBh3dQ/k9HnLXPxybQgwVJ9axXVNDx6FaaP+DUkVFbgwz4ozrFiZrFarXFNHS/dQv/+sp9Sydzl+sfOkTVkdSvW3qGFzULrc4r2FNTWJ3d4PVgv3zX4o4RcYE0WH49LpmPvKOdui8nSGFFBITrqtBbmjBqqpOO0/n3B3jdA78Vt02rr6H8bDgaHJEbzdwv3sNsVT4l340J4BhFdL20hI7H9mG8lZQNpPBwd1R4wD5/Ka6RzzdU4PvRYtK1Oi4jCvKXVU7z+bBfsZvOd9bwXw48lDGZd37yAI+1NqX9PhNVR8JOMs+CRH+N/UIDpV+IF2aOKtg3VaLt3J6kYB+4vPiUoQ84ilP2zExGfyWil5Sw6eTiuvEFT52N9k6tUChqY7jMx6bvdceVvfTfGrFnU19v3vdfQHidXnImf0VIrHex726Cjx+N2mUNDjkiz0VOHcT6OXKdWgrBaoip60ibOuoKLse0qaViNUJi20hkpckcgMpyijudp2Z7YtJRKMvPT9ptYXstLf8WitaR2dsX99QjTEfMxmhqyOn0L7mEXvfhw0kCymmRkroPH8b44Ymkr/TWrfHq8q7Cj7azBWtklOl6f/Sr6x88hFlVkvIUelsr+bGK6IB/zEYWFYAQUbV3vawUgJnGIszyQnrf1MDIO9qj471jv9tO1Y+BF+MbnNV5BeOHJ+IU3js/vwf/M2fwf/v5pJ5DuNLpTZg3Bmh95Fj6x9fimATlBmJkbFJYi6ruAFpBATI8v6xyIj+A9tNTND56BHkzvSBJ5ycPcNffLFMcWsqUPk5For+M2hp8IQttZ4tjsyv6HFGwny/ysFARoONP25h58IBzjIoKxMxiMrXGnOQufL6M/kpk5PXb8f74bPRmrxUFk/bRK8rwf/t5zBsD0W2jb2+n6rCdskesFxchw+HF79L5y427VPWuD4wz8J5FpXqjcQvXfrM+qgsa62etoMDZlkIYO7auI21Kb2sF4tvUkrGK0z4S20aEjsf2UfV/nk/5Xa5E/AfOTZ2D92CXJPtF++mpOOFssXuHc7N0fblQVcxCXTkjL6ui732HksqnI9v30BWJTE5hv2I3nqEp5rYU45kIw8QM8oKT5INfOsrwt7ZREpil+LN2cs/SfWS0fZ7o00Lk0bGwy+n5GzXV4GqHRtTZI4rom484j83m3jZ4/qyjIv+5eBX5WITXy/Q3a3hD7Quw80S8Erk7tDT7hv0UnnAU0ZdC5nnTfmf7dIyYcWh7agq9dStMTeVcTgsuPYEV+q2DbHnSXGzQsRMdq0UKfwHYpYUEekJgWo5GZH4e4tJlR9YsFML/bSdRbH12UcGeULwavZx3bnZ69eaobmj8CZJ2AAAYqElEQVRGf+Ek3rIXJrHD4egwnCxIToqxjP1uO/t+/xTXvyUo+0xq8fDQq7eT/9XnonMVaf1lW2nr3bzeS8WZTRT9rIyLw5VsKRmn+hcOR+2M9XN0+FCL6Vumqevgl44S+l4TM/PbKP5SijaVATOwOGEb2zZ6PtxO2591Y6aTFFwBeksT4YAH42jyzTiCtquNnZ+7yMmxKYxXE71pRuYHio9AMZDts8NLbgwdViDtlWJsVQ8GEeWl2AV+xOzi+KLe0rT4ORhEBAuxhkfiekV6WyuY1uJjZ4JdUdXzpgZHWDaWiO5njDBvNgifz9F+jH26wJkQ8vSOLjmhlA16SQnWZCiurrKZrEpZrmpztHcZGYON4uqzAhn1K3MmzRh6or8g+bo0vx+tojzttYo9O5Anzscp2SddV+z+afy18Et78f747JJSa6nGpiPj1npFxWK9CYFeVooIFmJe6VlSnzXytk66ejfv34Px7EknRluaYHAkeo1L+RlS17Xe0sTclhI8PziRtk2ltDXNRLpRW+PMb6TIgZl8kg1i9w7QQJ7IbmhXb2vFutCR9vtsx9BVD32FSCkhz4cYHIubaLG7r0U/W6EQRlkJ+qYKp6FHXosbGnUEY12Ez+usaZmAHfAnbdMKCpwx+/OXMO/fQ/8rfEgNLJ+k6dEjdP7vg9T8WBL4enxPTIbDziSUYSAtC+H1opeXYb/QiZmDxqowjOjYavTG404AWuPjSRNPcnZ2ReVSYlsMvb6Z2UqB5WtBSKg4bXKjXaf538ZzF1NORxp/OQbGJwPNvcEkHcK9Vn1wAllchD25mCysUAgtEEDUVsUlMEjwl2k6iSw/gHHhBmZioi4vw56aju80lJdhj46BpkU1dCM3B+GPGbeWEqu5BnHK6RxEhkTS+WvwtVtS1/sXx7HPvIjv7DUst26sy858UGTyN9HPRlMD0h12ylTXdvc18mbmMEnRptJg3r+Ha6/0gQQzb7FtVP0cQls0bM8WLB9o89D4pT7mGsvpf6WPho9n3/tPRK+oQF68jD2Xeo5C+HxohQVx12e9eDlavyvhJddD17c1YxXlwfNno72fpLdSYjCqNiNLglgXOqLvhEcShR4MwuYKGBh2hGkj6uTFRVBZ4TT2gWHscBjampEvvOgEts+HVlyE0LS4XglCoOXlIRpqEZPTYNssNG3GON8NHi9jv7SVsh9fR4bnEV4P2DZyYSEuMCJBEVUZDwZB1xF+H/b0TWdbRQXzO+vw9YxiF+YxuruE8u9dhvISbL+BuNiNPTPj9Kw8BgyPYw0Po7e1ImbmMHuuYTRuAdNCzs5iT02jNTcgxiaxJ0PO5N7O7WgTU9hj44iazXGvLS6nnLZzO2axH+1np51x6uYGRO+N1K+epnitbsm4aGuFG8NY4+NJfk7nL724CCrKwGMgxkNI24byEqwLHehtrSyU5OHp7AdbQmkRMs+LvHjZSbj778YYDiEnQqBrUFqM1dWDXpCf1l/6syejttqXr8Ylbi0QQCspBo+BnL7pvF1RUYY1PIpRVUm4dTPe093IcBhhGAz81g4qP3/KuclUVDjnm5hMenJN56+4uE2od62wEK2sxCmX584LVG0Gnxd7aCTOzzLPh5gNL46vZ6rrygqszm705gbsq71xbSpTrEbrKKFtrAVGQz0z2zYR6HCeXGTAD6aF3XMdvaoSNM152hYCo3ELcnwSLMu5qe/cjugbpO+td1Hz7QEYd2Nb13lq4B+USHQqtEAArawUs+8GekWZM5svZTSgEhGGgV5Xg9U/gF692QnQi53Oo2REYXzBTB6KKSlxGklZsRO4c3NO4nWHCdI9qgrDQGttgsjkYORR1R0GMGprwNCxiwsQV3qTbkSRR8XEnu7Cq/fgP3GF8O4mfH2TTi8t3et5keGc4iJESTH2yBhC1xAlxcjpmajtYs8OtM7rUFke90gNzvCREAJ7ax2i61pcA8q1nB4MIqVEKy/FvNqLUVkBfp/jt5gGu3iQzK8dpsJo3IKcnsEOheL9LLTM/tq9A31oHLOmDL2z10l8bi9LGAZaQT6yvgp5qRssy7npR3xZV+u8WRSaAk3E3ZxS+UuvqEBWlWGfSf/kod2zHbquYd3bgnGxxzmmEBgN9c65fF7nBjI6ERcfwjBA11MOX2TyV/zJneuKTHLL1nq0a0PM7mkgcLYv7vXNiJ+llNg3Zxx/LdE29JIS7OZa9N5hx9ZIm4IlYxXSt43VRAsEok9CcdsjHazysmgeQNpo996FNjiGXVnqTJK6dSh270Dr6Wd2/1bynr/MU2OfUQn9VhMJ7Fx7A1ogALqespzw+RBeb87HFIaByMtb1Z7JUsfUCguR8/NJSWK55TYqmezN5MuMfl6Bv9L1QJdbr8u1RRgGWiCQ9gaw3J7yWvaw15t0MZBY59mOoauErlAoFBsc9cMihUKheImhErpCoVDcIaiErlAoFHcIKqErFArFHcKSPywSQmwDvhyzqQn4U5xfpL4diLwX9AEp5XdW3UKFQqFQZEU266FfAnYBCCF0oA9Hhu5twCeklB9bUwsVCoVCkRW5Drn8InBZSnl1LYxRKBQKxfLJNaG/GXg85u93CyHOCCE+J4RIuT6sEOIRIcRxIcTxBW6PH4ooFArF7UjWCV0I4QVeD/y7u+kfga04wzE3gI+nKielfExKuVdKuddD6kXsFQqFQrFycumhvw44KaUcBJBSDkopLSmlDXwG2L8WBioUCoUiO3JJ6A8TM9wihKiK+e5B4NxqGaVQKBSK3MkqoQshAsBrgK/HbP5rIcRZIcQZ4FXAH62BfRuLyFrN64i+Y1vG742qzSm3RxYGUyyDNfRzRBgiaXsGP6/3d5mI08pNgRbIrJqkWFuyEriQUs4AZQnbfntNLFpjjMYtXHpXFfm9GsKG6s+fw47VUMyAiJFHA0c+quOtQQK9Gpbf0SsVPp8jKSYlPR9pdwQnvjaLOPIC4IhOz7fNYE17XIMkd/3VGHONpXiePh5n59XfrKHmo4fjbLjy0XbMUmdp3YKKm9T/cWpRChHIg5gV3Hq+fA/G6QJmqyxa/mBR9OLmQwcY2qthTAtKOmwKvnI0uc5SqA0ZjVu48N5K8q8ZzFTZtLznaNw65Jc/dhDvuNNfaPjKQNx66Mspp7c0cfGPysm/ZoCE+n86jz07F+e3bFSRMhEpn4ufOz61D++4nuTnCJ1/fwDpdQ+mSba9K144IZ2fI1z+m3ZaHxtMtnUNy6UTWu781AFa/9tpJI6uwI2/NpgYdPVg5zQavmWSd2nQWe88XZ25/or4OdbHsHSb6viH/XjHdbZ8eyaprjs+u5eCF73MVUia3nckbhnlmw8dYLxFJ9hjE3w8uxhfCX3vO8RCoaThQ6nrMlUeaHyCaA64/DfteCcFUoemL/RCd8rDJPGSW22x/08OUfN3zyNNMyoDFisXF4ebaHrff4javzqctM62/fJdaD9ztP/0libkwDCiujJZMGOJ9bnN+/fgG5iKk6Aa+MNDbP7kc0nlllo69Cu9R3hTbXt0zWdhGAy+cz9VTzkNJ7ImtPD5oK2ZqZbCxSQuBEZ1FdI0sScmkeEw4acb6DtZRdOj8YE5+8b9BL590lHR2bEN5hcWk3ZtDbNtVdHg1Lc1I6/2Ys/N5VzO7rmOtnULPQ+WU/+/XL+5jU9v3eqUdX205/XnGGxfhmyYEPQ+2r5YPhc/b9+yqG0ZUy5SzxE1pnSk8zPAzK8dwDthphSwXu9yfY8eYssXr2VMeqGHD0aTZbo6sy51xfk5IocXSfZLtqn9dzv1neCj8AP7CBzuwJqYRPP7EXl54PVE1z2ff+0+vN87hrarDfv0BbRAAGlZGWN8OejBIINv3kH5Y0ectc/HJtCDBamFWBLLtrXC6AQze7cgLPB+75gjFO/18v25L6rlc7Nh5sEDFHz3heykn2IbbOvWJLkwPRhE1ldjn3sxLrh733+Iuh9MIY8ly1rp25od+a+YY9186ADBC+NYFzsRuu4kP3dh/PAD+/CNhRk4UIBnRlL2z8+DbTH2tnZ8UzaF33auJXER/9jeUcdn97LtXan1KK1X3RdVxslEpCHGVU9ErLj9XsRz56J1Fbvmc67ltKJgSnWc6TcdJPit0yuW7EpJDn6OvZ5EP+vBIDOv2IZx0+Lq63xUHbHI+3+uurwQ3Pz1/Ul+jiDb70UzbWyvjvj56UUpuBzLRcmhnPbcubjv5l+7j7zrIaRHj1Oqjz12rELUUnUW6+eIIIRRtRlZmL9kmyo5OYJ1qSuurlP2rjUdvbkBq+NyVJHJvN4bbQdLaaauGvvvRr85v6Tub6p2Ef2uvIynhv9JaYouhbarjYLuqeUlhckp7FfsxjMYYq6+GO/kPEzcRL7oOCX4+FH6v9FGUd4cxZ+zk5O5mzRknjfaSCJBVtjl9MCN6ipHVguiyjy+7xwDYPNRV1h3bxs8fxYzICj9l+fA7+iPxiqbg6MiP/b1en6l9jzcezxeRd5tkHO/up+C031ko3+eVg0ekIaGFjNsEZWzm5rKuZxWUpy0n9izg8LLy/Rbrizh55G/FczOp/azFQrhe9LxV9OPQOy7G819GtSbG9P6WQsEkHMmWngBfXQe2dKEWDAxe67lXM6+1ocMh3MqR211NN6Mpgb06yFHWm5kElFcBJvK4xOvlHGi10vVWayfo+o+mpZVmxp5soXZ+TaKH1+sa2toGG3nduyAB23e4mZDAQUdE4hZN2FrYtHU+XkA9KrK1GpXq4je0sR8vmfxCS4VMXkgHVOvaI6fvczASzah69uaYd7MqLSdCWtwCG1wCIJB8kwLO+BHWPairmhLE9UPOolaD44gamviFcpty7Hh5mJSEl6v89jpJvioKG9TQ7IBQmDeGEAvCWIBmz7ljIlGklzisIw9NUXJr0xxxFeIvi1BRV5KZPu9TjJPI8WXSKqeWuScnnM9WHb8k1+kceVaLranCO5k3py5ZI9ntVjKz4UfD1LiDlGk9HMM8thZ7D07ABBz84t6qQl+nv2FHfifPoUVc+1i9w7oyb2cXlKCFQ7nVE7s3oHodZ4IZGgK60qPa4Rz4zfy8hYvKvI041lMJUvVWTo/Z9OmIvWtB4NxdS3PvejI2lVXUnBxDAwds/tq1B6jtsY5RmEhVii05slc29WGLcTST7sp8oBzAB3N68He1UrhMxezP+8ybL2t0cvLHJXx/sFlJ/NYpGki/V7E8Jgj/upidy8GjBUKgcdA3+S+4aDpzr/R8bixe+FLfZe2A/6kbVpBgXPsCx1Yr7qP6x86xLU/O0T3X7YD0PWJg8w8eCDZ3hgVeYRA+HwYNdXopy5lnczB1aCM2OI+FRg11Y5N4+PoFXFz6MgYnc1cyuEmRr28zOnlX7+xbsk8lnR+9l+6Ef2c5OcYIteqD06kPUfEz9q8jUxIeGZx+h/lZSon/LmXM4t9oGd40yc2Tm3L0Ss1k8fk09VZop+NpoaonzOVs7uvRes7VV3L+XkYHAbbxrrodFj08jL04iLsyVC0nBYIOLG0RugVFY4Y+KnzmXdMkwe0QACjpgopJeL4hfQaril4yfXQRUE+Xb9bxdaPjaO3bgVDjyqFp8Ko2owMFjgivXc1A0QTih4MOsK5vQNxwrB6MAibK5xxxcER7HAYuygf+4xTTngMtOKipHNZo2POhE5TPWJy2gnMhkrExR70igrGX7OV0p9cR4bnEYbuCBAvLMCzJ6l71jmG5vdjA61/ej4qsBxVkZ8MOcMfFRXM76jD1zOCXRBgZHcJZU+b6E312D4P4lL34jCJpsHYpKOkvq0ZEZ53Hv3rauLU4PUd22BsMjomKitKMTwe7JFRRF314tsqyyintzQhLJtL79hE60fGnR6NoUfV5yM+YmjUSRY7tuWc9PVtzYvlc/CzeWPQqacEP+ttrZglAYyufrAllATR/B6sS91ofj9zLZX4CgNJfjZqqjG6hpEtjY6Y89Aowu9H751ALqOc2def8/n03gmscBgtEMBsqcUoLXae4q73I+prYH7BiWOfD83nw7zWlyySnqbOEv0s83wwG45OGi7VpsyuHnR3KCm2TQFowSCYJgyPRu2JiEiz/270zl7nGkuKsC47vfdMMb4cjLpaZu/aTF6nk09kng9MC/tqr3Pz8RjRm1SqPKBXVGDXb2JkR5DijlKMgQlkwA9L3BsivKQmRZcjLisMA72uBqt/wFGD93qc8UMpHRXymRnkgpkc0CUlCL8PWVqEmJhyesau2ndGFXnDQNvaALYzkCwHhh2bI0rxtTVg6NjBAKKnP+nunU7Z3Lx/D75TV5jf1Yi3b9K5hnRv37iP1nowiCgrSaukLnbvQLt8HVFeihwZi7MlqureVAsdPXGNM5dywjCW7KFEfCSnZ7BDIbSGuqiPssVo3BItn4ufxZ4d6P2jyX527dIK8pG1VcjObrAspGlGxZOp2pTsZ5ykN/3Ld1P40y5m9jcRONKFNTGB0PXcy42P53w+33ePO3XnxoF82S6Mrn5md9WTd/paNK4i7SnVBGOmOovzs2liz4UXJ6GXaFN2UzV6/yhoWlxdGzXV0WGb5JMtth17MoTQtbgbSKYYz5XIDSvV9dpTU9GXG5bKA4kokegNQOTHPbneRDKqwft8CK93WTem5arIL/eY6VTml1tuo5LJz5l8mek7cMe/x8dvebnIDSGXR/9sWE4H61aU2whsyIQuhJgC1n8AdHUoB0ZutRHL4Ha1G25f229Xu+H2tf12tRuys32LlDL1z4xjWO8x9EvZ3GU2IkKI47ej7ber3XD72n672g23r+23q92wura/5N5yUSgUijsVldAVCoXiDmG9E/pj63y+1eR2tf12tRtuX9tvV7vh9rX9drUbVtH2dZ0UVSgUCsXaoYZcFAqF4g5h3RK6EOK1QohLQoguIcSj63Xe5SCE6HHFO04LIY6720qFEN8XQnS6/6cUxV5vXIHuISHEuZhtKW0VDp90fXBGCHHfrbM8re1/LoToc+v+tBDigZjv3u/afkkI8cu3xmoQQtQJIZ4VQlwUQpwXQrzH3b6h6z2D3bdDnfuFEM8LIV5wbf8Ld3ujEOI5t86/7GofI4TwuX93ud83bDC7Py+E6I6p813u9pXFipRyzf8BOnAZaAK8wAtA23qce5n29gDlCdv+GnjU/fwo8L9utZ2uLa8E7gPOLWUr8ADwXUAAB4HnNqDtfw78SYp929y48QGNbjzpt8juKuA+93Mh0OHat6HrPYPdt0OdC6DA/ewBnnPr8ivAm93tnwZ+3/38X4FPu5/fDHx5g9n9eeChFPuvKFbWq4e+H+iSUl6RUs4DTwBvWKdzrxZvAL7gfv4C8MZbaEsUKeVPgLGEzelsfQPwr9LhKFCcoA27rqSxPR1vAJ6QUoallN1AF7dImFxKeUNKedL9PAVcBGrY4PWewe50bKQ6l1LKafdPj/tPAvcDX3W3J9Z5xBdfBX5RCLG4ju46kcHudKwoVtYrodcA12P+7iVzIN1qJPC0EOKEEOIRd1ullPIGOA0D2HTLrFuadLbeLn54t/u4+bmYoa0Nabv7KL8bp+d129R7gt1wG9S5EEIXQpwGhoDv4zwxTEgpI+sFx9oXtd39fpIEGc31ItFuKWWkzj/i1vknhBCRZTFXVOfrldBT3Rk38us1L5NS3ge8DniXEOKVt9qgVeJ28MM/AluBXcAN4OPu9g1nuxCiAPga8IdSykwLnWwo21PYfVvUuZTSklLuAmpxnhTuSrWb+/+GsT3RbiHETuD9wHZgH1AKvM/dfUV2r1dC7wXqYv6uBdILLd5ipJT97v9DwDdwgmcw8ujj/p9iWbcNQzpbN7wfpJSDbgOwgc+w+Ii/oWwXQnhwkuIXpZQRPZkNX++p7L5d6jyClHIC+BHOGHOxECKyhEmsfVHb3e+LyH54b02Isfu17vCXlFKGgX9hlep8vRL6MaDFnZH24kxSfGudzp0TQoh8IURh5DPwS8A5HHvf4u72FuCbt8bCrEhn67eA33Fn0g8Ck5Ehgo1Cwnjhgzh1D47tb3bfXmgEWoDn19s+cN5EAP4ZuCil/NuYrzZ0vaez+zap8wohRLH7OQ94Nc4cwLPAQ+5uiXUe8cVDwA+lO+u4nqSx+8WYG7/AGfePrfPlx8o6zvY+gDOrfhn44Hqddxl2NuHM7L+As6z8B93tZcAzQKf7f+mtttW163Gcx+QFnLv776WzFedx7lOuD84Cezeg7f/m2nbGDe6qmP0/6Np+CXjdLbT75TiPwWeA0+6/BzZ6vWew+3ao83uAU66N54A/dbc34dxkuoB/B3zudr/7d5f7fdMGs/uHbp2fA/4vi2/CrChW1C9FFQqF4g5B/VJUoVAo7hBUQlcoFIo7BJXQFQqF4g5BJXSFQqG4Q1AJXaFQKO4QVEJXKBSKOwSV0BUKheIOQSV0hUKhuEP4/3s6n+5Y1ONKAAAAAElFTkSuQmCC\n",
      "text/plain": "<Figure size 432x288 with 1 Axes>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ]
  },
  {
   "metadata": {
    "trusted": true
   },
   "cell_type": "code",
   "source": "plt.hist(lum_img.ravel(), bins=256, range=(0.0, 1.0))",
   "execution_count": 6,
   "outputs": [
    {
     "data": {
      "text/plain": "(array([2.7917e+04, 2.2200e+02, 2.5900e+02, 1.8700e+02, 2.1900e+02,\n        9.3000e+01, 1.1000e+02, 1.9800e+02, 1.6300e+02, 8.9000e+01,\n        9.8000e+01, 9.7000e+01, 1.0400e+02, 1.2400e+02, 1.1900e+02,\n        8.9000e+01, 1.0700e+02, 1.0900e+02, 7.4000e+01, 9.3000e+01,\n        7.0000e+01, 1.6000e+02, 7.3000e+01, 1.2200e+02, 8.2000e+01,\n        7.4000e+01, 6.7000e+01, 3.9000e+01, 6.1000e+01, 4.4000e+01,\n        3.4000e+01, 3.8000e+01, 2.9000e+01, 2.0000e+01, 4.3000e+01,\n        2.5000e+01, 3.3000e+01, 2.2000e+01, 1.3000e+01, 3.1000e+01,\n        4.2000e+01, 2.3000e+01, 1.6000e+01, 1.4000e+01, 1.7000e+01,\n        1.8000e+01, 2.1000e+01, 1.4000e+01, 5.0000e+00, 2.0000e+01,\n        1.8000e+01, 1.1000e+01, 1.8000e+01, 7.0000e+00, 1.3000e+01,\n        7.0000e+00, 1.4000e+01, 7.0000e+00, 3.0000e+00, 6.0000e+00,\n        4.0000e+00, 6.0000e+00, 9.0000e+00, 3.0000e+00, 2.0000e+00,\n        7.0000e+00, 6.0000e+00, 7.0000e+00, 3.0000e+00, 2.0000e+00,\n        7.0000e+00, 5.0000e+00, 4.0000e+00, 8.0000e+00, 8.0000e+00,\n        4.0000e+00, 3.0000e+00, 1.1000e+01, 2.0000e+00, 5.0000e+00,\n        1.9000e+01, 1.6000e+01, 1.5000e+01, 2.0000e+01, 2.5000e+01,\n        1.4000e+01, 1.6000e+01, 8.0000e+00, 7.0000e+00, 0.0000e+00,\n        4.0000e+00, 6.0000e+00, 4.0000e+00, 9.0000e+00, 9.0000e+00,\n        8.0000e+00, 1.2000e+01, 1.4000e+01, 2.0000e+00, 9.0000e+00,\n        8.0000e+00, 1.1000e+01, 9.0000e+00, 1.0000e+01, 4.0000e+00,\n        1.1000e+01, 1.0000e+00, 1.0000e+01, 1.0000e+00, 2.9000e+01,\n        9.0000e+00, 1.6000e+01, 4.0000e+00, 7.0000e+00, 1.0000e+00,\n        5.0000e+00, 0.0000e+00, 1.0000e+00, 4.0000e+00, 0.0000e+00,\n        0.0000e+00, 7.0000e+00, 0.0000e+00, 4.0000e+00, 3.0000e+00,\n        2.0000e+00, 7.0000e+00, 2.0000e+00, 3.0000e+00, 2.0000e+00,\n        1.0000e+00, 7.0000e+00, 2.0000e+00, 0.0000e+00, 9.0000e+00,\n        4.0000e+00, 3.0000e+00, 3.0000e+00, 0.0000e+00, 2.0000e+00,\n        0.0000e+00, 1.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        1.0000e+00, 0.0000e+00, 0.0000e+00, 6.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 5.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 1.0000e+00, 0.0000e+00, 0.0000e+00, 1.0000e+00,\n        3.0000e+00, 0.0000e+00, 0.0000e+00, 2.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 1.0000e+00,\n        0.0000e+00, 1.0000e+00, 0.0000e+00, 2.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        1.0000e+00, 0.0000e+00, 0.0000e+00, 3.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 2.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00,\n        0.0000e+00]),\n array([0.        , 0.00390625, 0.0078125 , 0.01171875, 0.015625  ,\n        0.01953125, 0.0234375 , 0.02734375, 0.03125   , 0.03515625,\n        0.0390625 , 0.04296875, 0.046875  , 0.05078125, 0.0546875 ,\n        0.05859375, 0.0625    , 0.06640625, 0.0703125 , 0.07421875,\n        0.078125  , 0.08203125, 0.0859375 , 0.08984375, 0.09375   ,\n        0.09765625, 0.1015625 , 0.10546875, 0.109375  , 0.11328125,\n        0.1171875 , 0.12109375, 0.125     , 0.12890625, 0.1328125 ,\n        0.13671875, 0.140625  , 0.14453125, 0.1484375 , 0.15234375,\n        0.15625   , 0.16015625, 0.1640625 , 0.16796875, 0.171875  ,\n        0.17578125, 0.1796875 , 0.18359375, 0.1875    , 0.19140625,\n        0.1953125 , 0.19921875, 0.203125  , 0.20703125, 0.2109375 ,\n        0.21484375, 0.21875   , 0.22265625, 0.2265625 , 0.23046875,\n        0.234375  , 0.23828125, 0.2421875 , 0.24609375, 0.25      ,\n        0.25390625, 0.2578125 , 0.26171875, 0.265625  , 0.26953125,\n        0.2734375 , 0.27734375, 0.28125   , 0.28515625, 0.2890625 ,\n        0.29296875, 0.296875  , 0.30078125, 0.3046875 , 0.30859375,\n        0.3125    , 0.31640625, 0.3203125 , 0.32421875, 0.328125  ,\n        0.33203125, 0.3359375 , 0.33984375, 0.34375   , 0.34765625,\n        0.3515625 , 0.35546875, 0.359375  , 0.36328125, 0.3671875 ,\n        0.37109375, 0.375     , 0.37890625, 0.3828125 , 0.38671875,\n        0.390625  , 0.39453125, 0.3984375 , 0.40234375, 0.40625   ,\n        0.41015625, 0.4140625 , 0.41796875, 0.421875  , 0.42578125,\n        0.4296875 , 0.43359375, 0.4375    , 0.44140625, 0.4453125 ,\n        0.44921875, 0.453125  , 0.45703125, 0.4609375 , 0.46484375,\n        0.46875   , 0.47265625, 0.4765625 , 0.48046875, 0.484375  ,\n        0.48828125, 0.4921875 , 0.49609375, 0.5       , 0.50390625,\n        0.5078125 , 0.51171875, 0.515625  , 0.51953125, 0.5234375 ,\n        0.52734375, 0.53125   , 0.53515625, 0.5390625 , 0.54296875,\n        0.546875  , 0.55078125, 0.5546875 , 0.55859375, 0.5625    ,\n        0.56640625, 0.5703125 , 0.57421875, 0.578125  , 0.58203125,\n        0.5859375 , 0.58984375, 0.59375   , 0.59765625, 0.6015625 ,\n        0.60546875, 0.609375  , 0.61328125, 0.6171875 , 0.62109375,\n        0.625     , 0.62890625, 0.6328125 , 0.63671875, 0.640625  ,\n        0.64453125, 0.6484375 , 0.65234375, 0.65625   , 0.66015625,\n        0.6640625 , 0.66796875, 0.671875  , 0.67578125, 0.6796875 ,\n        0.68359375, 0.6875    , 0.69140625, 0.6953125 , 0.69921875,\n        0.703125  , 0.70703125, 0.7109375 , 0.71484375, 0.71875   ,\n        0.72265625, 0.7265625 , 0.73046875, 0.734375  , 0.73828125,\n        0.7421875 , 0.74609375, 0.75      , 0.75390625, 0.7578125 ,\n        0.76171875, 0.765625  , 0.76953125, 0.7734375 , 0.77734375,\n        0.78125   , 0.78515625, 0.7890625 , 0.79296875, 0.796875  ,\n        0.80078125, 0.8046875 , 0.80859375, 0.8125    , 0.81640625,\n        0.8203125 , 0.82421875, 0.828125  , 0.83203125, 0.8359375 ,\n        0.83984375, 0.84375   , 0.84765625, 0.8515625 , 0.85546875,\n        0.859375  , 0.86328125, 0.8671875 , 0.87109375, 0.875     ,\n        0.87890625, 0.8828125 , 0.88671875, 0.890625  , 0.89453125,\n        0.8984375 , 0.90234375, 0.90625   , 0.91015625, 0.9140625 ,\n        0.91796875, 0.921875  , 0.92578125, 0.9296875 , 0.93359375,\n        0.9375    , 0.94140625, 0.9453125 , 0.94921875, 0.953125  ,\n        0.95703125, 0.9609375 , 0.96484375, 0.96875   , 0.97265625,\n        0.9765625 , 0.98046875, 0.984375  , 0.98828125, 0.9921875 ,\n        0.99609375, 1.        ]),\n <a list of 256 Patch objects>)"
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAD8CAYAAACcjGjIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAEURJREFUeJzt3H/MnWV9x/H3Z1ScmzqqLYaUujJTE6vJEBvsYrKhLFAwsZjgAolSDVmNg0U3s4juD4xIglvUhERxGBrKogLzx2i0rmsYi3MR5FEZP0d4hgwqBKpFdCHTod/9ca66s16nfU6fX+d52vcrOTn3+d7XfZ/v1edpP71/nJOqQpKkYb826QYkSUuP4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6hgOkqTOikk3MFurVq2qdevWTboNSVpWvvOd7/ywqlbPNG7ZhsO6deuYmpqadBuStKwk+c9xxnlaSZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUOSbDYd1lX5t0C5K0pB2T4SBJOjzDQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSR3DQZLUMRwkSZ0ZwyHJ2iS3JXkgyX1J3tvqH07ygyR3tce5Q9t8MMl0kgeTnD1U39xq00kuG6qfkuSOJA8luSnJ8fM9UUnS+MY5cngOeH9VvQrYBFySZENb98mqOrU9dgG0dRcArwY2A59OclyS44BPAecAG4ALh/bzsbav9cDTwMXzND9J0izMGA5V9URVfbct/xR4AFhzmE22ADdW1c+q6vvANHB6e0xX1cNV9XPgRmBLkgBvAr7Ytt8BnDfbCUmS5u6IrjkkWQe8FrijlS5NcneS7UlWttoa4LGhzfa22qHqLwV+XFXPHVQf9f7bkkwlmdq3b9+RtC5JOgJjh0OSFwJfAt5XVT8BrgFeAZwKPAF8/MDQEZvXLOp9seraqtpYVRtXr149buuSpCO0YpxBSZ7HIBg+V1VfBqiqJ4fWfxb4anu5F1g7tPnJwONteVT9h8AJSVa0o4fh8ZKkCRjnbqUA1wEPVNUnhuonDQ17K3BvW94JXJDk+UlOAdYD3wbuBNa3O5OOZ3DRemdVFXAbcH7bfitwy9ymJUmai3GOHN4AvAO4J8ldrfYhBncbncrgFNAjwLsBquq+JDcD9zO40+mSqvoFQJJLgd3AccD2qrqv7e8DwI1JPgp8j0EYSZImZMZwqKpvMvq6wK7DbHMlcOWI+q5R21XVwwzuZpIkLQF+QlqS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEmdGcMhydoktyV5IMl9Sd7b6i9JsifJQ+15ZasnydVJppPcneS0oX1tbeMfSrJ1qP66JPe0ba5OkoWYrCRpPOMcOTwHvL+qXgVsAi5JsgG4DLi1qtYDt7bXAOcA69tjG3ANDMIEuBx4PXA6cPmBQGljtg1tt3nuU5MkzdaM4VBVT1TVd9vyT4EHgDXAFmBHG7YDOK8tbwFuqIHbgROSnAScDeypqv1V9TSwB9jc1r24qr5VVQXcMLQvSdIEHNE1hyTrgNcCdwAvq6onYBAgwIlt2BrgsaHN9rba4ep7R9RHvf+2JFNJpvbt23ckrUuSjsDY4ZDkhcCXgPdV1U8ON3RErWZR74tV11bVxqrauHr16plaliTN0ljhkOR5DILhc1X15VZ+sp0Soj0/1ep7gbVDm58MPD5D/eQRdUnShIxzt1KA64AHquoTQ6t2AgfuONoK3DJUv6jdtbQJeKaddtoNnJVkZbsQfRawu637aZJN7b0uGtqXJGkCVowx5g3AO4B7ktzVah8CrgJuTnIx8CjwtrZuF3AuMA08C7wLoKr2J7kCuLON+0hV7W/L7wGuB14AfL09JEkTMmM4VNU3GX1dAODMEeMLuOQQ+9oObB9RnwJeM1MvkqTF4SekJUkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1JkxHJJsT/JUknuHah9O8oMkd7XHuUPrPphkOsmDSc4eqm9utekklw3VT0lyR5KHktyU5Pj5nKAk6ciNc+RwPbB5RP2TVXVqe+wCSLIBuAB4ddvm00mOS3Ic8CngHGADcGEbC/Cxtq/1wNPAxXOZkCRp7mYMh6r6BrB/zP1tAW6sqp9V1feBaeD09piuqoer6ufAjcCWJAHeBHyxbb8DOO8I5yBJmmdzueZwaZK722mnla22BnhsaMzeVjtU/aXAj6vquYPqkqQJmm04XAO8AjgVeAL4eKtnxNiaRX2kJNuSTCWZ2rdv35F1LEka26zCoaqerKpfVNUvgc8yOG0Eg//5rx0aejLw+GHqPwROSLLioPqh3vfaqtpYVRtXr149m9YlSWOYVTgkOWno5VuBA3cy7QQuSPL8JKcA64FvA3cC69udScczuGi9s6oKuA04v22/FbhlNj1JkubPipkGJPkCcAawKsle4HLgjCSnMjgF9AjwboCqui/JzcD9wHPAJVX1i7afS4HdwHHA9qq6r73FB4Abk3wU+B5w3bzNTpI0KzOGQ1VdOKJ8yH/Aq+pK4MoR9V3ArhH1h/m/01KSpCXAT0hLkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpYzhIkjqGgySpM2M4JNme5Kkk9w7VXpJkT5KH2vPKVk+Sq5NMJ7k7yWlD22xt4x9KsnWo/rok97Rtrk6S+Z6kJOnIjHPkcD2w+aDaZcCtVbUeuLW9BjgHWN8e24BrYBAmwOXA64HTgcsPBEobs21ou4PfS5K0yGYMh6r6BrD/oPIWYEdb3gGcN1S/oQZuB05IchJwNrCnqvZX1dPAHmBzW/fiqvpWVRVww9C+JEkTMttrDi+rqicA2vOJrb4GeGxo3N5WO1x974i6JGmC5vuC9KjrBTWL+uidJ9uSTCWZ2rdv3yxblCTNZLbh8GQ7JUR7fqrV9wJrh8adDDw+Q/3kEfWRquraqtpYVRtXr149y9YlSTOZbTjsBA7ccbQVuGWoflG7a2kT8Ew77bQbOCvJynYh+ixgd1v30ySb2l1KFw3tS5I0IStmGpDkC8AZwKokexncdXQVcHOSi4FHgbe14buAc4Fp4FngXQBVtT/JFcCdbdxHqurARe73MLgj6gXA19tDkjRBM4ZDVV14iFVnjhhbwCWH2M92YPuI+hTwmpn6kCQtHj8hLUnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpI7hIEnqGA6SpM6cwiHJI0nuSXJXkqlWe0mSPUkeas8rWz1Jrk4yneTuJKcN7WdrG/9Qkq1zm5Ikaa7m48jhjVV1alVtbK8vA26tqvXAre01wDnA+vbYBlwDgzABLgdeD5wOXH4gUCRJk7EQp5W2ADva8g7gvKH6DTVwO3BCkpOAs4E9VbW/qp4G9gCbF6AvSdKY5hoOBfxjku8k2dZqL6uqJwDa84mtvgZ4bGjbva12qLokaUJWzHH7N1TV40lOBPYk+ffDjM2IWh2m3u9gEEDbAF7+8pcfaa+SpDHN6cihqh5vz08BX2FwzeDJdrqI9vxUG74XWDu0+cnA44epj3q/a6tqY1VtXL169VxalyQdxqzDIclvJnnRgWXgLOBeYCdw4I6jrcAtbXkncFG7a2kT8Ew77bQbOCvJynYh+qxWkyRNyFxOK70M+EqSA/v5fFX9Q5I7gZuTXAw8Crytjd8FnAtMA88C7wKoqv1JrgDubOM+UlX759CXJGmOZh0OVfUw8Lsj6j8CzhxRL+CSQ+xrO7B9tr1IkuaXn5CWJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSZ8mEQ5LNSR5MMp3kskn3I0nHsiURDkmOAz4FnANsAC5MsmGyXUnSsWvFpBtoTgemq+phgCQ3AluA+xfqDddd9rWu9shVb/7VugPLknQsWirhsAZ4bOj1XuD1i93EcGCMCo/58shVb/5/+z84iA6sGxVQh1snSfNlqYRDRtSqG5RsA7a1l/+V5MFZvt8q4Iez3HbO8rHDv56pPtO6Q5jonCfEOR/9jrX5wtzn/NvjDFoq4bAXWDv0+mTg8YMHVdW1wLVzfbMkU1W1ca77WU6c87HhWJvzsTZfWLw5L4kL0sCdwPokpyQ5HrgA2DnhniTpmLUkjhyq6rkklwK7geOA7VV134TbkqRj1pIIB4Cq2gXsWqS3m/OpqWXIOR8bjrU5H2vzhUWac6q6676SpGPcUrnmIElaQo7qcJjpKzmSPD/JTW39HUnWLX6X82eM+f55kvuT3J3k1iRj3dK2lI37tStJzk9SSZb9nS3jzDnJH7Wf9X1JPr/YPc63MX63X57ktiTfa7/f506iz/mSZHuSp5Lce4j1SXJ1+/O4O8lp895EVR2VDwYXtv8D+B3geODfgA0HjfkT4DNt+QLgpkn3vcDzfSPwG235Pct5vuPOuY17EfAN4HZg46T7XoSf83rge8DK9vrESfe9CHO+FnhPW94APDLpvuc4598HTgPuPcT6c4GvM/iM2Cbgjvnu4Wg+cvjVV3JU1c+BA1/JMWwLsKMtfxE4M8moD+QtBzPOt6puq6pn28vbGXyeZDkb52cMcAXwV8B/L2ZzC2ScOf8x8Kmqehqgqp5a5B7n2zhzLuDFbfm3GPE5qeWkqr4B7D/MkC3ADTVwO3BCkpPms4ejORxGfSXHmkONqarngGeAly5Kd/NvnPkOu5jB/zyWsxnnnOS1wNqq+upiNraAxvk5vxJ4ZZJ/TXJ7ks2L1t3CGGfOHwbenmQvg7se/3RxWpuYI/37fsSWzK2sC2Ccr+QY62s7lomx55Lk7cBG4A8WtKOFd9g5J/k14JPAOxeroUUwzs95BYNTS2cwODr8lySvqaofL3BvC2WcOV8IXF9VH0/ye8Dftjn/cuHbm4gF/7fraD5yGOcrOX41JskKBoejhzuUW8rG+gqSJH8I/CXwlqr62SL1tlBmmvOLgNcA/5zkEQbnZncu84vS4/5e31JV/1NV3wceZBAWy9U4c74YuBmgqr4F/DqD7yA6Wo31930ujuZwGOcrOXYCW9vy+cA/VbvaswzNON92iuVvGATDcj8PDTPMuaqeqapVVbWuqtYxuM7ylqqamky782Kc3+u/Z3DzAUlWMTjN9PCidjm/xpnzo8CZAElexSAc9i1ql4trJ3BRu2tpE/BMVT0xn29w1J5WqkN8JUeSjwBTVbUTuI7B4ec0gyOGCybX8dyMOd+/Bl4I/F277v5oVb1lYk3P0ZhzPqqMOefdwFlJ7gd+AfxFVf1ocl3PzZhzfj/w2SR/xuD0yjuX8X/0SPIFBqcFV7XrKJcDzwOoqs8wuK5yLjANPAu8a957WMZ/fpKkBXI0n1aSJM2S4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6vwv7skWJXglvBwAAAAASUVORK5CYII=\n",
      "text/plain": "<Figure size 432x288 with 1 Axes>"
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "name": "python3",
   "display_name": "Python 3",
   "language": "python"
  },
  "language_info": {
   "name": "python",
   "version": "3.6.4",
   "mimetype": "text/x-python",
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "pygments_lexer": "ipython3",
   "nbconvert_exporter": "python",
   "file_extension": ".py"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "version_major": 2,
    "version_minor": 0,
    "state": {}
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
No results found