ponynet.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Here is the notebook where we developed intuition about the problem\n",
    "# This is, we believe, a good example of visualization driven developement for deep learning\n",
    "\n",
    "%pylab inline\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def pony(I, r):\n",
    "    assert len(I.shape) == 3\n",
    "    assert r>0\n",
    "    r = int(r)\n",
    "    O = np.zeros((I.shape[0]*r, I.shape[1]*r, I.shape[2]/(r*2)))\n",
    "    for x in range(O.shape[0]):\n",
    "        for y in range(O.shape[1]):\n",
    "            for c in range(O.shape[2]):\n",
    "                c += 1\n",
    "                a = np.floor(x/r).astype(\"int\")\n",
    "                b = np.floor(y/r).astype(\"int\")\n",
    "                d = c*r*(y%r) + c*(x%r)\n",
    "                print a, b, d\n",
    "                O[x, y, c-1] = I[a, b, d]\n",
    "    return O"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "I = np.concatenate([np.ones((8, 8, 1))*0, 20*np.ones((8, 8, 1)), 30*np.ones((8, 8, 1)), 40*np.ones((8, 8, 1))], axis=-1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f9532c97b50>"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADM1JREFUeJzt3X+s3XV9x/Hnq9RFlJ/+gQQaKMyBC4myJkMncesEY4MJ\n+M+S4hIF/1UguhgQkrHyx2JMjCMw/yAqQQMzsZrAH252TVOSbdGJpfws0AQDBUaNcc6YJYyW9/44\np3ip95b2nO/ne+j9PB/JTe75cb+vz5fyOud7vufc+05VIakvaxa9AEnjs/hShyy+1CGLL3XI4ksd\nsvhSh9a2Dkji+4XSglRVlru+efEn/naF63cCG5e95bU7bpspac11K2WZN1beat63ReTVnbPl5XNH\nWMtMWzy04WRTkqeSPJPkxnm2JWk8Mxc/yRrgTuBjwEXA1UneO9TCJLUzzzP+JcDeqnquql4Fvgtc\ndWybWD9H/CzMO37zxsxa/XnzFP9sYN+Syy9MrzsG6+eIn4V5x2/emFmrP8+386QOzXNW/0XgnCWX\n102vW8bOJd+vZ/xHU2n12/kM7Nx7dPedp/g/Bd6T5Fzgv4DNwNXL33XjHDGSjsbGCyZfh2z555Xv\nO3Pxq+pgks8B25i8ZPhmVe2ZdXuSxjPXB3iq6l+ACwdai6SReHJP6pDFlzpk8aUOjfJLOrcc/L9j\n/pktM65slizzhs1bzfu2iLy/O2G2vCPxGV/qkMWXOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDFl/q\nkMWXOmTxpQ5ZfKlDo/ySzt+v/YNj/pkDt8+WtXaGLPOGzVvN+7aIvIN3zJa35bqVb/MZX+qQxZc6\nNM8IrXVJdiR5IsljSa4fcmGS2pnnNf4B4AtVtTvJScDPkmyrqqcGWpukRmZ+xq+ql6tq9/T73wJ7\nOOYRWpIWYZDX+EnWAxcDPxlie5LamvvtvOlh/lbghukz/++p2rnk0nomjxOShrRz7zgjtEiylknp\nv1NV9698v43zxEg6Chv/aPJ1yG1HGKE176H+t4Anq2rGjzRIWoR53s67FPhr4CNJHk6yK8mm4ZYm\nqZV5hmb+O9DgL35Las1P7kkdsvhSh0b57bybDxz76KDbZlzZLFnmDZu3mvdtEXlbHKElaQgWX+qQ\nxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUPOzjNv8LzVvG+LyHN2nqRB\nWHypQ3MXP8ma6R/afGCIBUlqb4hn/BuAJwfYjqSRzFX8JOuAK4BvDLMcSWOY9xn/a8AXgRpgLZJG\nMvPbeUk+DuyfjsneCGSl+zo7T2pvrNl5lwJXJrkCOBE4Ocm3q+pTh9/R2XlSe6PMzquqm6vqnKo6\nH9gM7Fiu9JLeenwfX+rQIB/ZraoHgQeH2Jak9nzGlzpk8aUOOTvPvMHzVvO+LSLP2XmSBmHxpQ5Z\nfKlDFl/qkMWXOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDFl/qkMWXOuTsPPMGz1vN+7aIPGfnSRqE\nxZc6NO8knVOTfC/JniRPJPnAUAuT1M68r/FvB35YVX+VZC3wjgHWJKmxeSbpnAJ8uKquAaiqA8Bv\nBlqXpIbmOdQ/D/hlkrunY7LvSnLiUAuT1M48h/prgQ3AZ6vqoST/ANwE3Hr4HZ2dJ7U31uy8F4B9\nVfXQ9PJW4Mbl7ujsPKm9sWbn7Qf2JblgetVlwJOzbk/SeOY9q389cG+StwHPAtfOvyRJrc1V/Kp6\nBPjTgdYiaSR+ck/qkMWXOuTsPPMGz1vN+7aIPGfnSRqExZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8\nqUMWX+qQxZc6ZPGlDll8qUMWX+qQs/PMGzxvNe/bIvKcnSdpEBZf6tC8s/M+n+TxJI8muTfJbMcy\nkkY1c/GTnAVcB2yoqvcxOV+weaiFSWpn3pN7JwDvTPIak4GZL82/JEmtzTNQ4yXgq8DzwIvAr6tq\n+1ALk9TOPNNyTwOuAs4F/gfYmuSTVXXf4fd1dp7U3liz8y4Hnq2qXwEk+QHwIeD3iu/sPKm9UWbn\nMTnE/2CStycJk9l5e+bYnqSRzPMa/z+ZTMh9GHgECHDXQOuS1NC8s/O2AFsGWoukkfjJPalDFl/q\nkLPzzBs8bzXv2yLynJ0naRAWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGl\nDll8qUPOzjNv8LzVvG+LyHN2nqRBWHypQ29a/CTfTLI/yaNLrjs9ybYkTyf5UZJT2y5T0pCO5hn/\nbuBjh113E7C9qi4EdgBfGnphktp50+JX1b8B/33Y1VcB90y/vwf4xMDrktTQrK/xz6iq/QBV9TJw\nxnBLktTaUG/n1RFvdHae1NwYs/P2J3l3Ve1PcibwiyPd2dl5UnstZudl+nXIA8A10+8/Ddx/DOuT\ntGBH83befcB/ABckeT7JtcCXgY8meZrJsMwvt12mpCG96aF+VX1yhZsuH3gtkkbiJ/ekDll8qUPO\nzjNv8LzVvG+LyHN2nqRBWHypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk\n8aUOOTvPvMHzVvO+LSLP2XmSBmHxpQ7NOjvvK0n2JNmd5PtJTmm7TElDmnV23jbgoqq6GNiLs/Ok\n48pMs/OqantVvTa9+GNgXYO1SWpkiNf4nwGOMLND0lvNXG/nJbkFeLWq7jvS/ZydJ7U3xuw8klwD\nXAF85M3vu3HWGElH6Vhm5x1t8d8wOy/JJuCLwJ9X1SuzLFLS4sw6O+8O4CTgX5PsSvL1xuuUNKBZ\nZ+fd3WAtkkbiJ/ekDll8qUPOzjNv8LzVvG+LyHN2nqRBWHypQxZf6pDFlzpk8aUOWXypQxZf6pDF\nlzpk8aUOWXypQxZf6pDFlzpk8aUOOTvPvMHzVvO+LSLP2XmSBjHTCK0lt/1NkteSvKvN8iS1MOsI\nLZKsAz4KPDf0oiS1NdMIramvMfkT25KOMzO9xk9yJbCvqh4beD2SRnDMZ/WTnAjczOQw//WrB1uR\npOZmeTvvD4H1wCNJwmRS7s+SXFJVv1juB5ydJ7XXYnbe6yO0qupx4MzXb0h+DmyoquXOA0zvs/Eo\nYyTN6lhm5806QmupwkN96bgy6witpbefP9xyJI3BT+5JHbL4UocsvtQhZ+eZN3jeat63ReQ5O0/S\nICy+1CGLL3XI4ksdsvhShyy+1CGLL3XI4ksdsvhShxZa/Od2jvt3Os07fvNW874tIm+xxX/wefPM\ne8tl9ZDnob7UoVF+Secszlr2+pM5ecXbTtiwYaasgytsz7zx8lbzvh1PeezateJNqarZNnqUkrQN\nkLSiqlr2z+I1L76ktx5f40sdsvhShxZS/CSbkjyV5JkkNzbOWpdkR5InkjyW5PqWeUty1yTZleSB\nEbJOTfK9JHum+/mBxnmfT/J4kkeT3JtktsHvK2//9yY0Jzk9ybYkTyf5UZJTG+d9Zfrfc3eS7yc5\npWXekttGmUA9evGTrAHuZDKB9yLg6iTvbRh5APhCVV0E/Bnw2cZ5h9wAPDlCDsDtwA+r6o+B9wN7\nWgUlOQu4jskQlfcxeWdo88Axy01ovgnYXlUXAjuALzXO2wZcVFUXA3tHyBt1AvUinvEvAfZW1XNV\n9SrwXeCqVmFV9XJV7Z5+/1smpTi7VR68/g94BfCNljnTrFOAD1fV3QBVdaCqftM49gTgnUnWAu8A\nXhpy4ytMaL4KuGf6/T3AJ1rmVdX2qnptevHHTEbFNcubGm0C9SKKfzawb8nlF2hcxEMyGdp3MfCT\nxlGH/gHHeMvkPOCXSe6evrS4azrYtImqegn4KvA88CLw66ra3ipviTOqav90DS8DZ4yQechngCMM\npJrf2BOouzm5l+QkYCtww/SZv1XOx4H906OM12cONrQW2AD8Y1VtAP6XyWFxE0lOY/Lsey5wFnBS\nkiNOW2pklPehk9wCvFpV9zXMODSB+talV7fKg8UU/0XgnCWX102va2Z6SLoV+E5V3d8yC7gUuDLJ\ns8A/AX+Z5NsN815g8kzx0PTyViYPBK1cDjxbVb+qqoPAD4APNcw7ZH+SdwMkORNYdjLzkJJcw+Ql\nW+sHtqUTqH/O7yZQNzuqWUTxfwq8J8m507PBm4HWZ76/BTxZVbc3zqGqbq6qc6YzBTcDO6rqUw3z\n9gP7klwwveoy2p5UfB74YJK3T8ekX0abk4mHHy09AFwz/f7TwNAP4G/IS7KJycu1K6vqlYGz3pBX\nVY9X1ZlVdX5VncfkwfxPVho7P4iqGv0L2AQ8zeRs6U2Nsy4FDgK7gYeBXcCmkfbzL4AHRsh5P5MH\n1N1MnoFPbZx3K5OyP8rkRNvbBt7+fUxOGL7C5IHmWuB0YPv0/5ttwGmN8/YyObu+a/r19ZZ5h93+\nLPCulv+GfmRX6lA3J/ck/Y7Flzpk8aUOWXypQxZf6pDFlzpk8aUOWXypQ/8PzlW1q8n+n1cAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532dac610>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imshow(O[:, :, 0], interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f9532bdd850>"
      ]
     },
     "execution_count": 49,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADGJJREFUeJzt3V2sZXV5x/HvD0YjiLyYiDhMYEMt2JAonaRoJbangnGC\nCeNNk8EmCt5aINoYEJsiXhETY0nVC6ISNFATRxO4wDqdTMakbbTiMLwOOAkODEwZY6w1pgnl5enF\n3oMbOGc4s9fLYeb//SQnWXutddbzX3Pmt9fL3llPqgpJbTlurQcgaXwGX2qQwZcaZPClBhl8qUEG\nX2rQuqELJPHzQmmNVFWWmz948Kf+YYX5O4GlZZf8/Yv/t1ClL37t5hWXfeEe+MJlKyy7eqFyfHHF\nfYNW9+9Y3jcYf/9uWnD/bjrMsk6n+kk2JXk0yS+SXNdlW5LGs3DwkxwHfBX4MHABcEWSd/U1MEnD\n6XLEvwjYW1VPVNVzwHeBzUe2iUmH8kdu6Y9HLYf7159jed9g/P3rEvwzgf1zr5+azTsCkw7lj5z/\nefo2Ga3SsbxvcHQFX9JRqstd/aeBs+Zeb5jNW8bOuekJY7+bSi3YN/tZjS7B/xnwziRnA/8FbAGu\nWH7VpQ5lJK3GhJcfUn98mHUXDn5VvZDkb4FtTC8ZvllVexbdnqTxdPoCT1X9C3B+T2ORNBJv7kkN\nMvhSgwy+1CCDLzXI4EsNMvhSgwy+1CCDLzXI4EsNMvhSgwy+1CCDLzXI4EsNMvhSgwy+1CCDLzXI\n4EsNMvhSgwy+1KAuLbQ2JNmR5OEkDya5ps+BSRpOl4dtPg98pqp2JzkJ+HmSbVX1aE9jkzSQhY/4\nVfVMVe2eTf8e2MMRt9CStBZ6ucZPMgEuBH7ax/YkDavTc/UBZqf5W4FrZ0f+Zeycm55gCy2pf/sY\np4UWSdYxDf13ququlddc6lJG0ipMWH0Lra6n+t8CHqmqWzpuR9KIunycdzHwN8AHk9yXZFeSTf0N\nTdJQujTN/Hfg+B7HImkkfnNPapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUG\nGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGdQ5+kuNmD9q8u48BSRpeH0f8a4FHetiOpJF0\nCn6SDcBlwDf6GY6kMXQ94n8F+CxQPYxF0kgWfq5+ko8AB2dtspeArLz2zrnpCfbOk/q3j3F6510M\nXJ7kMuAE4C1Jvl1VH3/1qksdykhajQkj9M6rqhuq6qyqOhfYAuxYPvSSXm/8HF9qUKc22YdU1Y85\n/JmFpNcRj/hSgwy+1CCDLzXI4EsNMvhSgwy+1CCDLzXI4EsNMvhSgwy+1CCDLzXI4EsNMvhSgwy+\n1CCDLzXI4EsNMvhSgwy+1CCDLzWoayedU5J8L8meJA8neW9fA5M0nK4P27wFuKeq/jrJOuDEHsYk\naWBdOumcDHygqq4EqKrngd/1NC5JA+pyqn8O8Oskt83aZN+a5IS+BiZpOF1O9dcBG4FPVdW9Sf4R\nuB648dWr7pybnmDvPKl/+xind95TwP6qunf2eitw3fKrLnUoI2k1JozTO+8gsD/JebNZlwCPLLo9\nSePpelf/GuCOJG8AHgeu6j4kSUPrFPyquh/4s57GImkkfnNPapDBlxpk8KUGGXypQQZfapDBlxpk\n8KUGGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGde2d9+kkDyV5\nIMkdSd7Y18AkDWfh4CdZD1wNbKyqdzN9cOeWvgYmaThdH699PPDmJC8ybZh5oPuQJA2tS0ONA8CX\ngSeBp4HfVtX2vgYmaThduuWeCmwGzgb+B9ia5GNVdeer1945Nz3B3nlS//YxTu+8S4HHq+o3AEl+\nALwfWCb4Sx3KSFqNCSP0zmN6iv++JG9KEqa98/Z02J6kkXS5xv9Pph1y7wPuBwLc2tO4JA2oa++8\nm4CbehqLpJH4zT2pQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXyp\nQQZfapDBlxpk8KUGGXypQQZfapDBlxr0msFP8s0kB5M8MDfvtCTbkjyW5EdJThl2mJL6tJoj/m3A\nh18x73pge1WdD+wAPtf3wCQN5zWDX1X/Bvz3K2ZvBm6fTd8OfLTncUka0KLX+KdX1UGAqnoGOL2/\nIUkaWtduuYfU4RfvnJueYO88qX/7GL533sEkb6+qg0nOAH51+NWXFiwjabUm9N87L7OfQ+4GrpxN\nfwK4a5XbkfQ6sJqP8+4E/gM4L8mTSa4CbgY+lOQxps0ybx52mJL69Jqn+lX1sRUWXdrzWCSNxG/u\nSQ0y+FKDDL7UIIMvNcjgSw0y+FKDDL7UIIMvNcjgSw0y+FKDDL7UIIMvNcjgSw0y+FKDDL7UIIMv\nNcjgSw0y+FKDDL7UoEV7530pyZ4ku5N8P8nJww5TUp8W7Z23Dbigqi4E9mLvPOmoslDvvKraXlUv\nzl7+BNgwwNgkDaSPa/xPAj/sYTuSRtKpd16SzwPPVdWdh19z59z0BHvnSf3bx/C980hyJXAZ8MHX\nXntp0TKSVmnC6nvnrTb4L+udl2QT8FngL6rq2SMcn6Q1tmjvvH8CTgL+NcmuJF8feJySerRo77zb\nBhiLpJH4zT2pQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXypQQZfapDBlxpk8KUGGXypQQZf\napDBlxpk8KUGGXypQQZfatBCLbTmlv1dkheTvHWY4UkawqIttEiyAfgQ8ETfg5I0rIVaaM18hekj\ntiUdZRa6xk9yObC/qh7seTySRnDEnXSSnADcwPQ0/6XZvY1I0uAWaaH1R0w79dyfJEw75f48yUVV\n9avlf2Xn3PQEe+dJ/dtH/73zXmqhVVUPAWe8tCD5JbCxqpa7DzCztMoykhY1YfW98xZtoTWv8FRf\nOqos2kJrfvm5/Q1H0hj85p7UIIMvNcjgSw0y+FKDDL7UIIMvNcjgSw0y+FKDDL7UIIMvNcjgSw0y\n+FKDDL7UoDUO/r5Rq+3cO2o53L/+HMv7BuPvn8Ef1L5Rqx3L+3cs7xs0F3xJa2GRZ+4dsY0b37Hs\n/AMH3sL69csvW89zixV728aVl514AN62ftlF7zjMrx3ORpYfP7S7f8fyvsHRs3/s2rXiolTVgltd\nnSTDFpC0oqpa9rF4gwdf0uuP1/hSgwy+1KA1CX6STUkeTfKLJNcNXGtDkh1JHk7yYJJrhqw3V/e4\nJLuS3D1CrVOSfC/Jntl+vnfgep9O8lCSB5LckeSNPW//VR2ak5yWZFuSx5L8KMkpA9f70uzfc3eS\n7yc5ech6c8tG6UA9evCTHAd8lWkH3guAK5K8a8CSzwOfqaoLgD8HPjVwvUOuBR4ZoQ7ALcA9VfUn\nwHuAPUMVSrIeuJppE5V3M/1kaEvPZZbr0Hw9sL2qzgd2AJ8buN424IKquhDYO0K9UTtQr8UR/yJg\nb1U9UVXPAd8FNg9VrKqeqards+nfMw3FmUPVg5f+gJcB3xiyzqzWycAHquo2gKp6vqp+N3DZ44E3\nJ1kHnAgc6HPjK3Ro3gzcPpu+HfjokPWqantVvTh7+ROmreIGqzczWgfqtQj+mcD+uddPMXAQD0ky\nAS4EfjpwqUN/wDE+MjkH+HWS22aXFrfOGpsOoqoOAF8GngSeBn5bVduHqjfn9Ko6OBvDM8DpI9Q8\n5JPAD4csMHYH6mZu7iU5CdgKXDs78g9V5yPAwdlZxks9Bwe0DtgIfK2qNgL/y/S0eBBJTmV69D0b\nWA+clOSw3ZYGMsrn0Ek+DzxXVXcOWONQB+ob52cPVQ/WJvhPA2fNvd4wmzeY2SnpVuA7VXXXkLWA\ni4HLkzwO/DPwV0m+PWC9p5geKe6dvd7K9I1gKJcCj1fVb6rqBeAHwPsHrHfIwSRvB0hyBrBCZ+b+\nJLmS6SXb0G9s8x2of8kfOlAPdlazFsH/GfDOJGfP7gZvAYa+8/0t4JGqumXgOlTVDVV11qyn4BZg\nR1V9fMB6B4H9Sc6bzbqEYW8qPgm8L8mbZm3SL2GYm4mvPFu6G7hyNv0JoO838JfVS7KJ6eXa5VX1\nbM+1Xlavqh6qqjOq6tyqOofpm/mfrtx2vgdVNfoPsAl4jOnd0usHrnUx8AKwG7gP2AVsGmk//xK4\ne4Q672H6hrqb6RH4lIHr3cg07A8wvdH2hp63fyfTG4bPMn2juQo4Ddg++3+zDTh14Hp7md5d3zX7\n+fqQ9V6x/HHgrUP+Df3KrtSgZm7uSfoDgy81yOBLDTL4UoMMvtQggy81yOBLDTL4UoP+H96/a4zL\nqiRhAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532c3c1d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imshow(I.reshape(16, 16), interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(128, 2)\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f953298a650>"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAADEAAAD/CAYAAABYfj7YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAACIVJREFUeJztnWuMXVUVx3//23GkPFplKjPQSouiIE0c4ENFUayl4ICJ\n+Mngu62vRJCJJEpLzVTmg1SNIW2wH4xSHgEbKlHGqFiaBhpNiLxqK53WYkM7HbSK2CpVyzyWH/ae\n3DOPO4+77+Pc3f1PbnLPvvvsddZZ/3Xvvud/1j4yMxodhXofQCWQnMgLkhN5QXJiMkjqkLRP0p8k\n3VYtOwCYWcVfuJPzIrAQeAOwC7i4GrbMrGqRWAIcMLNDZjYAbAFuqJKtqjkxH+jLbB/xbVVBU7UG\nngqSxs13zEzljFWtSPQD52e2F/i2UejqgKuA61tbg4xVy4mngQslLZTUDNwI9FTJVnWcMLMh4GZg\nG/ACsMXMesf2637sKqzrSuzLlwXZq1pOmNljwEWT91rEwqUDnE0Lv+4u31bdEhtgaOP93LEM2trb\ng8aJYtpR10g0dXaxdvh1hjkXCp3lj1PBY5oxBjd0c0cBColOkdApRSIUKScyiCInEp3ygijolCIR\nipQTGUSRE4lOeUEUdEqRCEWlciLRKS+Igk5RRKJsJyQtkLRD0guS9ki6xbe/WdI2Sfsl/UbS3FJj\nNHV2MTy8muHnV5V7GEBYJAaBW81sMfBe4CZJFwOrge1mdhGwA1gTdITTQQV1up8Dy4F9QKtvawP2\nlehvQxuxLrBN7e3mDqWOmp2kRcClwFPegaP+BP0VOKfUfpWiU/C3k6QzgZ8CnWb22gQyVsnbeMzu\n5f6lRjPNQccQFAlJTTgHHjCzR33zUUmt/vM24G+l9h/aeJgVO/u49fiZIYcRTKd7gL1mtiHT1gOs\n8O8/Bzw6dqdKQ+XetCXpSmAnsAdHGQNuB34PPAy8FTgEfNzMjk2wv0ldrB16nTbO5eZCZ9nqadk5\nYWa/A2aV+Hh5ueOWgyimHWkWmxdEQacUiWDjKSeKiCInEp3ygijolCIRbDzlRBFR5ESiU14QBZ1S\nJIKNp5woIoqcSHTKC6KgUxSRCHZCUkHSc5J6/HZDyV0j6AT2ZrYbS+7Cldk8DiwFenxbw8lddwFf\nZ7Qa1Dhyl6SPAEfNbJekpZN0rbrcFUKlbwOHgYPAX4DXgAeAXkbTqTe3dDKz283sfDN7G656a4eZ\nfQb4BTWWu6rxO7EeuEbSfuBqvz0h6p4TWZjZk8CT/v2rJLlr5kiz2LwgCjqlSAQbTzlRRBQ5keiU\nF0RBpxSJYOMpJ4qIIicSnfKCKOiUIhFsPOVEEVHkRKJTXhAFnaKIRGgRyFxJWyX1+iqv99RD7gpV\niu4FVvr3TcBc4DvAN3zbbcD63F7alzQH+ICZbfYnY9DMjuMWLbzPd7sP+Fi5NqaLEDpdALwiabNX\nT38o6XQaSe7y+14O3GRmz0i6C6ecNlR11xGgz8ye8duP4JxqnOouT5k+Se/0TVfj1vxrnOouAEnt\nwI9wy6keBFbiiqVqWt0V9BUb+PVc/6/YPCHNYvOCNIvNC6LIiUSnvCAKOqVIhCLlRAZR5ESiU14Q\nBZ1SJEKRciKDKHIi0SkviIJOKRKhSDmRQajc9TVJf5S0W9KDkpobqrpL0nnAV4HLzezdOGp+ghlU\ndw1u6KZQWE/hsnvKPQwgnE6zgDP8KnSzcQ9Uahy5y8xeBr6PqyvqB46b2XYaSe6S9CbcWV8IHAe2\nSvoUDSZ3LQcOmtmr5p4Q9TPgfTSS3IWj0RWSTpMknNy1lwaUu9bhytMGgOeBLwBnkeSuJHfVwXia\nxRaRZrF5QRQ5keiUF0RBpxSJUKScyCCKnEh0yguioFOKRChSTmQQRU4kOuUFUdApRSIUqVgwgynp\nJOnHko5K2p1pKylpSVoj6YCv+Lq27CObCaZxzfT9uKdB7c60TVjBBVyCu7DcBCwCXsRftK7rtVgz\n+y3wzzHNpSStjwJbzFV6vQQcAJbM+MzOEOV+O51TQtKaD/Rl+vX7tglRd7lrDMoSOeotd5WStPpx\n4soIFvi2CVFruUv+NYJSklYPcKMX5S8ALsQ9BmtC1IxOkh7CrULaIukwsA63QOFWSavwkhaAme2V\n9DBOuxsAvmIhetp0keSunODUmHY0AtIsNhSJThlEQacUiWDjKSeKiCInEp3ygijolCIRbDzlRBFR\n5ESiU14QBZ1SJIKNp5wo4tTIiRJK0Xe9ErRL0iN+JbqRz6atFNWyMGoz8OExbduAxWZ2KU5IWQMg\n6RLcddl3AdcBm3xtRVVRllJkZtvNbNhvPoW7hA8zVIryUt0FsAr4lX8/I6WoUghKbElrgQEz+0lZ\nxut9g4qkFcD1wLJM84yUopo+uwsn5+7JbHfg1vxrGdNvRAJuxi07WRMJeDoOPAS8DJzEVXStxCXs\nIeA5/9qU6b/GH3wvcO1kIovUZd8cXm13D28IcmJKOpnZJydo3jxJ/zuBO6cat5I4NX6xGwFpFpsX\nRJETdXVi1i2f5dM7BjhGCyz7Q9nj1NUJeIlDTwzxj3F3Is0MdXWiq2MnT3SDWluDxqlvYp/Rxllt\nbZw+b17QMEH12EGGpXGGrcx67Lo5UUlE8TuRnMgLauKEpA5JJ9ziCjJJJ8euceMv/Qxn+gxIOjKd\n8avuhKQCcDfujs1f4v5czQK+l+lzHe5f4THgf7g/VM+a2YJxA06AWkRiCcW7NdcB/wL+y+g1bm4A\nTgB/xznbBJw9csfnVKiFE/NxB9gEvIK7wRFgzpg+zcB5uIV93gG8HfjWdAxUdNoh6XEge/aEW4xk\n3OIjjL8h+ATu6uF+3BMUFgOfl/Sgv4BXEhV1wsyuGdsm6Qrcs+UHgXm4hwsA/DvTrR84DZhjZgcl\nzcadgP8AHwQmdaIWdHoaeAswBHwRR6PZjF7jpgdoAb4k6UO4k2vAG/3+k6NG98B24M6qZV4ngD8D\nzwK7cQ8JH858fhL4wXTGT3OnvCA5kRckJ/KC5ERe8H8hdLjMWj5wOgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532ace410>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "A = I.reshape(-1, 2, 2)\n",
    "A = np.concatenate([i for i in A], axis=0)\n",
    "print A.shape\n",
    "imshow(A, interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "B = I.reshape(8, 8, 2, 2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(8, 16, 2)\n"
     ]
    }
   ],
   "source": [
    "B = np.concatenate([B[i] for i in range(8)], axis=1)\n",
    "print B.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(16, 2, 8)\n"
     ]
    }
   ],
   "source": [
    "B = B.transpose(1, 2, 0)\n",
    "print B.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(16, 16)\n"
     ]
    }
   ],
   "source": [
    "B = np.concatenate([B[:, :, i] for i in range(8)], axis=1)\n",
    "print B.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f95327eb610>"
      ]
     },
     "execution_count": 71,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADN9JREFUeJzt3W/MnXV9x/H3p5RFhPHHRGTQQGEGXEiUNRk6iVsnGBtM\nwCdLiksUfOqE6GJAWEYgW2ZMjGNzPiAqQQMzsZrAA/90TVMyXXRiKX8LNsFBgVFjnDNmCaPw3YNz\nWm/KfZf2nOt3HXr/3q/kTu7z5z6f30X5nOs61zn3/U1VIakvaxa9AEnjs/hShyy+1CGLL3XI4ksd\nsvhSh9a2Dkji+4XSglRVlru+efEn/maF63cAG5e95aaXX5wp6e/WHH+YW80bI281b9si8v66Zsv7\n2/z9irfNdaifZFOSx5P8NMn18zyWpPHMXPwka4AvAO8HLgSuSvK2oRYmqZ159vgXA3uq6qmqehH4\nOnDl0T3E+jniZ2HesZs3Ztbqz5un+GcBe5dcfmZ63VFYP0f8LMw7dvPGzFr9eb6dJ3VonrP6zwJn\nL7m8bnrdMnYs+X494z+bSqvff+54iqd2PH1E952n+D8G3prkHOC/gM3AVcvfdeMcMZKOxPqN57B+\n4zkHL//bLd9f8b4zF7+qXkryl8BWJi8ZvlxVu2d9PEnjmesDPFX1XeCCgdYiaSSe3JM6ZPGlDll8\nqUOj/JLOS/9461H/zK3HzZo128+ZN1zeat62ReTdsuzv183HPb7UIYsvdcjiSx2y+FKHLL7UIYsv\ndcjiSx2y+FKHLL7UIYsvdcjiSx2y+FKHRvklnbXXrTRJZ2U37v+/2bLW/s5MP2fecHmredsWkXfT\ny7PlseYzK9802yNKOpZZfKlD84zQWpdke5JHkzyc5NohFyapnXle4+8HPllVu5KcBPwkydaqenyg\ntUlqZOY9flU9X1W7pt//BtjNUY/QkrQIg7zGT7IeuAj40RCPJ6mtud/Omx7mbwGum+75X6Vqx5JL\n65k8T0ga0lgjtEiylknpv1ZV96x8v43zxEg6Aq8aoXXrD1a877yH+l8BHquq2+Z8HEkjmuftvEuA\nvwDem+SBJDuTbBpuaZJamWdo5g+AGf/CuKRF8pN7UocsvtShUX47b/9tM4zQmnFl+2c8zWjecHmr\nedsWkXdLg92ze3ypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk8aUOWXyp\nQ87OM2/wvNW8bYvIc3aepEFYfKlDcxc/yZrpH9q8d4gFSWpviD3+dcBjAzyOpJHMVfwk64DLgS8N\nsxxJY5h3j/954FNADbAWSSOZ+e28JB8A9k3HZG8EstJ9nZ0ntTfW7LxLgCuSXA6cAPxukq9W1YcP\nvaOz86T2RpmdV1U3VtXZVXUesBnYvlzpJb3++D6+1KFBPrJbVfcB9w3xWJLac48vdcjiSx1ydp55\ng+et5m1bRJ6z8yQNwuJLHbL4UocsvtQhiy91yOJLHbL4UocsvtQhiy91yOJLHbL4UocsvtQhiy91\nyNl55g2et5q3bRF5zs6TNAiLL3Vo3kk6pyT5RpLdSR5N8s6hFiapnXlf498GfLuq/jzJWuCNA6xJ\nUmPzTNI5GXhPVV0NUFX7gV8PtC5JDc1zqH8u8Iskd0zHZN+e5IShFiapnXkO9dcCG4CPVdX9Sf4B\nuAG4+dA7OjtPam+s2XnPAHur6v7p5S3A9cvd0dl5Untjzc7bB+xNcv70qkuBx2Z9PEnjmfes/rXA\nXUmOB54Erpl/SZJam6v4VfUg8EcDrUXSSPzkntQhiy91yNl55g2et5q3bRF5zs6TNAiLL3XI4ksd\nsvhShyy+1CGLL3XI4ksdsvhShyy+1CGLL3XI4ksdsvhShyy+1CFn55k3eN5q3rZF5Dk7T9IgLL7U\noXln530iySNJHkpyV5LZjmUkjWrm4ic5E/g4sKGq3s7kfMHmoRYmqZ15T+4dB5yY5GUmAzOfm39J\nklqbZ6DGc8DngKeBZ4FfVdW2oRYmqZ15puWeClwJnAP8D7AlyYeq6u5D7+vsPKm9sWbnXQY8WVW/\nBEjyLeDdwKuK7+w8qb1RZucxOcR/V5I3JAmT2Xm753g8SSOZ5zX+fzCZkPsA8CAQ4PaB1iWpoXln\n590C3DLQWiSNxE/uSR2y+FKHnJ1n3uB5q3nbFpHn7DxJg7D4UocsvtQhiy91yOJLHbL4UocsvtQh\niy91yOJLHbL4UocsvtQhiy91yOJLHXJ2nnmD563mbVtEnrPzJA3C4ksdes3iJ/lykn1JHlpy3WlJ\ntiZ5Isn3kpzSdpmShnQke/w7gPcfct0NwLaqugDYDnx66IVJauc1i19V3wf++5CrrwTunH5/J/DB\ngdclqaFZX+OfXlX7AKrqeeD04ZYkqbWh3s6rw97o7DypuTFm5+1L8paq2pfkDODnh7uzs/Ok9lrM\nzsv064B7gaun338EuOeoVihpoY7k7by7gX8Hzk/ydJJrgM8A70vyBJNhmSt/REjS685rHupX1YdW\nuOmygdciaSR+ck/qkMWXOuTsPPMGz1vN27aIPGfnSRqExZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8\nqUMWX+qQxZc6ZPGlDll8qUMWX+qQs/PMGzxvNW/bIvKcnSdpEBZf6tCss/M+m2R3kl1Jvpnk5LbL\nlDSkWWfnbQUurKqLgD04O086psw0O6+qtlXVy9OLPwTWNVibpEaGeI3/UeA7AzyOpJHM9XZekpuA\nF6vq7sPdz9l5UntjzM4jydXA5cB7X/u+G2eNkXSEjmZ23pEW/xWz85JsAj4F/ElVvTDbMiUtyqyz\n8/4JOAn41yQ7k3yx8TolDWjW2Xl3NFiLpJH4yT2pQxZf6pCz88wbPG81b9si8pydJ2kQFl/qkMWX\nOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDzs4zb/C81bxti8hzdp6kQcw0\nQmvJbX+V5OUkb2qzPEktzDpCiyTrgPcBTw29KEltzTRCa+rzTP7EtqRjzEyv8ZNcAeytqocHXo+k\nERz1Wf0kJwA3MjnMP3j1YCuS1Nwsb+f9PrAeeDBJmEzK/UmSi6vq58v9gLPzpPZazM47OEKrqh4B\nzjh4Q/IzYENVLXceYHqfjUcYI2lWRzM7b9YRWksVHupLx5RZR2gtvf284ZYjaQx+ck/qkMWXOmTx\npQ45O8+8wfNW87YtIs/ZeZIGYfGlDll8qUMWX+qQxZc6ZPGlDll8qUMWX+qQxZc6tNDi79hjnnmv\nv6we8hZa/PtG3ljzjt281bxti8jzUF/q0Ci/pMObNyx//YnPwZvPXPam31vhR1476zC3mTdK3mre\ntmMqb+fOFW9KVc34qEcmSdsASSuqqmX/LF7z4kt6/fE1vtQhiy91aCHFT7IpyeNJfprk+sZZ65Js\nT/JokoeTXNsyb0numiQ7k9w7QtYpSb6RZPd0O9/ZOO8TSR5J8lCSu5LMNvh95cd/1YTmJKcl2Zrk\niSTfS3JK47zPTv977kryzSQnt8xbctsoE6hHL36SNcAXmEzgvRC4KsnbGkbuBz5ZVRcCfwx8rHHe\nAdcBj42QA3Ab8O2q+gPgHcDuVkFJzgQ+zmSIytuZvDO0eeCY5SY03wBsq6oLgO3ApxvnbQUurKqL\ngD0j5I06gXoRe/yLgT1V9VRVvQh8HbiyVVhVPV9Vu6bf/4ZJKc5qlQcH/wEvB77UMmeadTLwnqq6\nA6Cq9lfVrxvHHgecmGQt8EbguSEffIUJzVcCd06/vxP4YMu8qtpWVS9PL/6Qyai4ZnlTo02gXkTx\nzwL2Lrn8DI2LeEAmQ/suAn7UOOrAP+AYb5mcC/wiyR3Tlxa3TwebNlFVzwGfA54GngV+VVXbWuUt\ncXpV7Zuu4Xng9BEyD/go8J2WAWNPoO7m5F6Sk4AtwHXTPX+rnA8A+6ZHGQdnDja0FtgA/HNVbQD+\nl8lhcRNJTmWy9z0HOBM4Kclhpy01Msr70EluAl6sqrsbZhyYQH3z0qtb5cFiiv8scPaSy+um1zUz\nPSTdAnytqu5pmQVcAlyR5EngX4A/S/LVhnnPMNlT3D+9vIXJE0ErlwFPVtUvq+ol4FvAuxvmHbAv\nyVsAkpwBLDuZeUhJrmbykq31E9vSCdQ/47cTqJsd1Syi+D8G3prknOnZ4M1A6zPfXwEeq6oZ/7L5\nkauqG6vq7OlMwc3A9qr6cMO8fcDeJOdPr7qUticVnwbeleQN0zHpl9LmZOKhR0v3AldPv/8IMPQT\n+Cvykmxi8nLtiqp6YeCsV+RV1SNVdUZVnVdV5zJ5Mv/DlcbOD6KqRv8CNgFPMDlbekPjrEuAl4Bd\nwAPATmDTSNv5p8C9I+S8g8kT6i4me+BTGufdzKTsDzE50Xb8wI9/N5MThi8weaK5BjgN2Db9/2Yr\ncGrjvD1Mzq7vnH59sWXeIbc/Cbyp5b+hH9mVOtTNyT1Jv2XxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlD\nFl/q0P8DB9C6WTZf9sYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532ace490>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imshow(B, interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def pony2(I, r):\n",
    "    B = I.reshape(8, 8, 2, 2).transpose(0, 1, 3, 2)\n",
    "    B = np.concatenate([B[i] for i in range(8)], axis=1)\n",
    "    B = B.transpose(1, 2, 0)\n",
    "    B = np.concatenate([B[:, :, i] for i in range(8)], axis=1)\n",
    "    return B"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f953276ee90>"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADM1JREFUeJzt3X+s3XV9x/Hnq9RFlJ/+gQQaKMyBC4myJkMncesEY4MJ\n+M+S4hIF/1UguhgQkrHyx2JMjCMw/yAqQQMzsZrAH252TVOSbdGJpfws0AQDBUaNcc6YJYyW9/44\np3ip95b2nO/ne+j9PB/JTe75cb+vz5fyOud7vufc+05VIakvaxa9AEnjs/hShyy+1CGLL3XI4ksd\nsvhSh9a2Dkji+4XSglRVlru+efEn/naF63cCG5e95bU7bpspac11K2WZN1beat63ReTVnbPl5XNH\nWMtMWzy04WRTkqeSPJPkxnm2JWk8Mxc/yRrgTuBjwEXA1UneO9TCJLUzzzP+JcDeqnquql4Fvgtc\ndWybWD9H/CzMO37zxsxa/XnzFP9sYN+Syy9MrzsG6+eIn4V5x2/emFmrP8+386QOzXNW/0XgnCWX\n102vW8bOJd+vZ/xHU2n12/kM7Nx7dPedp/g/Bd6T5Fzgv4DNwNXL33XjHDGSjsbGCyZfh2z555Xv\nO3Pxq+pgks8B25i8ZPhmVe2ZdXuSxjPXB3iq6l+ACwdai6SReHJP6pDFlzpk8aUOjfJLOrcc/L9j\n/pktM65slizzhs1bzfu2iLy/O2G2vCPxGV/qkMWXOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDFl/q\nkMWXOmTxpQ5ZfKlDo/ySzt+v/YNj/pkDt8+WtXaGLPOGzVvN+7aIvIN3zJa35bqVb/MZX+qQxZc6\nNM8IrXVJdiR5IsljSa4fcmGS2pnnNf4B4AtVtTvJScDPkmyrqqcGWpukRmZ+xq+ql6tq9/T73wJ7\nOOYRWpIWYZDX+EnWAxcDPxlie5LamvvtvOlh/lbghukz/++p2rnk0nomjxOShrRz7zgjtEiylknp\nv1NV9698v43zxEg6Chv/aPJ1yG1HGKE176H+t4Anq2rGjzRIWoR53s67FPhr4CNJHk6yK8mm4ZYm\nqZV5hmb+O9DgL35Las1P7kkdsvhSh0b57bybDxz76KDbZlzZLFnmDZu3mvdtEXlbHKElaQgWX+qQ\nxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUPOzjNv8LzVvG+LyHN2nqRB\nWHypQ3MXP8ma6R/afGCIBUlqb4hn/BuAJwfYjqSRzFX8JOuAK4BvDLMcSWOY9xn/a8AXgRpgLZJG\nMvPbeUk+DuyfjsneCGSl+zo7T2pvrNl5lwJXJrkCOBE4Ocm3q+pTh9/R2XlSe6PMzquqm6vqnKo6\nH9gM7Fiu9JLeenwfX+rQIB/ZraoHgQeH2Jak9nzGlzpk8aUOOTvPvMHzVvO+LSLP2XmSBmHxpQ5Z\nfKlDFl/qkMWXOmTxpQ5ZfKlDFl/qkMWXOmTxpQ5ZfKlDFl/qkMWXOuTsPPMGz1vN+7aIPGfnSRqE\nxZc6NO8knVOTfC/JniRPJPnAUAuT1M68r/FvB35YVX+VZC3wjgHWJKmxeSbpnAJ8uKquAaiqA8Bv\nBlqXpIbmOdQ/D/hlkrunY7LvSnLiUAuT1M48h/prgQ3AZ6vqoST/ANwE3Hr4HZ2dJ7U31uy8F4B9\nVfXQ9PJW4Mbl7ujsPKm9sWbn7Qf2JblgetVlwJOzbk/SeOY9q389cG+StwHPAtfOvyRJrc1V/Kp6\nBPjTgdYiaSR+ck/qkMWXOuTsPPMGz1vN+7aIPGfnSRqExZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8\nqUMWX+qQxZc6ZPGlDll8qUMWX+qQs/PMGzxvNe/bIvKcnSdpEBZf6tC8s/M+n+TxJI8muTfJbMcy\nkkY1c/GTnAVcB2yoqvcxOV+weaiFSWpn3pN7JwDvTPIak4GZL82/JEmtzTNQ4yXgq8DzwIvAr6tq\n+1ALk9TOPNNyTwOuAs4F/gfYmuSTVXXf4fd1dp7U3liz8y4Hnq2qXwEk+QHwIeD3iu/sPKm9UWbn\nMTnE/2CStycJk9l5e+bYnqSRzPMa/z+ZTMh9GHgECHDXQOuS1NC8s/O2AFsGWoukkfjJPalDFl/q\nkLPzzBs8bzXv2yLynJ0naRAWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGlDll8qUMWX+qQxZc6ZPGl\nDll8qUPOzjNv8LzVvG+LyHN2nqRBWHypQ29a/CTfTLI/yaNLrjs9ybYkTyf5UZJT2y5T0pCO5hn/\nbuBjh113E7C9qi4EdgBfGnphktp50+JX1b8B/33Y1VcB90y/vwf4xMDrktTQrK/xz6iq/QBV9TJw\nxnBLktTaUG/n1RFvdHae1NwYs/P2J3l3Ve1PcibwiyPd2dl5UnstZudl+nXIA8A10+8/Ddx/DOuT\ntGBH83befcB/ABckeT7JtcCXgY8meZrJsMwvt12mpCG96aF+VX1yhZsuH3gtkkbiJ/ekDll8qUPO\nzjNv8LzVvG+LyHN2nqRBWHypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk8aUOWXypQxZf6pDFlzpk\n8aUOOTvPvMHzVvO+LSLP2XmSBmHxpQ7NOjvvK0n2JNmd5PtJTmm7TElDmnV23jbgoqq6GNiLs/Ok\n48pMs/OqantVvTa9+GNgXYO1SWpkiNf4nwGOMLND0lvNXG/nJbkFeLWq7jvS/ZydJ7U3xuw8klwD\nXAF85M3vu3HWGElH6Vhm5x1t8d8wOy/JJuCLwJ9X1SuzLFLS4sw6O+8O4CTgX5PsSvL1xuuUNKBZ\nZ+fd3WAtkkbiJ/ekDll8qUPOzjNv8LzVvG+LyHN2nqRBWHypQxZf6pDFlzpk8aUOWXypQxZf6pDF\nlzpk8aUOWXypQxZf6pDFlzpk8aUOOTvPvMHzVvO+LSLP2XmSBjHTCK0lt/1NkteSvKvN8iS1MOsI\nLZKsAz4KPDf0oiS1NdMIramvMfkT25KOMzO9xk9yJbCvqh4beD2SRnDMZ/WTnAjczOQw//WrB1uR\npOZmeTvvD4H1wCNJwmRS7s+SXFJVv1juB5ydJ7XXYnbe6yO0qupx4MzXb0h+DmyoquXOA0zvs/Eo\nYyTN6lhm5806QmupwkN96bgy6witpbefP9xyJI3BT+5JHbL4UocsvtQhZ+eZN3jeat63ReQ5O0/S\nICy+1CGLL3XI4ksdsvhShyy+1CGLL3XI4ksdsvhShxZa/Od2jvt3Os07fvNW874tIm+xxX/wefPM\ne8tl9ZDnob7UoVF+Secszlr2+pM5ecXbTtiwYaasgytsz7zx8lbzvh1PeezateJNqarZNnqUkrQN\nkLSiqlr2z+I1L76ktx5f40sdsvhShxZS/CSbkjyV5JkkNzbOWpdkR5InkjyW5PqWeUty1yTZleSB\nEbJOTfK9JHum+/mBxnmfT/J4kkeT3JtktsHvK2//9yY0Jzk9ybYkTyf5UZJTG+d9Zfrfc3eS7yc5\npWXekttGmUA9evGTrAHuZDKB9yLg6iTvbRh5APhCVV0E/Bnw2cZ5h9wAPDlCDsDtwA+r6o+B9wN7\nWgUlOQu4jskQlfcxeWdo88Axy01ovgnYXlUXAjuALzXO2wZcVFUXA3tHyBt1AvUinvEvAfZW1XNV\n9SrwXeCqVmFV9XJV7Z5+/1smpTi7VR68/g94BfCNljnTrFOAD1fV3QBVdaCqftM49gTgnUnWAu8A\nXhpy4ytMaL4KuGf6/T3AJ1rmVdX2qnptevHHTEbFNcubGm0C9SKKfzawb8nlF2hcxEMyGdp3MfCT\nxlGH/gHHeMvkPOCXSe6evrS4azrYtImqegn4KvA88CLw66ra3ipviTOqav90DS8DZ4yQechngCMM\npJrf2BOouzm5l+QkYCtww/SZv1XOx4H906OM12cONrQW2AD8Y1VtAP6XyWFxE0lOY/Lsey5wFnBS\nkiNOW2pklPehk9wCvFpV9zXMODSB+talV7fKg8UU/0XgnCWX102va2Z6SLoV+E5V3d8yC7gUuDLJ\ns8A/AX+Z5NsN815g8kzx0PTyViYPBK1cDjxbVb+qqoPAD4APNcw7ZH+SdwMkORNYdjLzkJJcw+Ql\nW+sHtqUTqH/O7yZQNzuqWUTxfwq8J8m507PBm4HWZ76/BTxZVbc3zqGqbq6qc6YzBTcDO6rqUw3z\n9gP7klwwveoy2p5UfB74YJK3T8ekX0abk4mHHy09AFwz/f7TwNAP4G/IS7KJycu1K6vqlYGz3pBX\nVY9X1ZlVdX5VncfkwfxPVho7P4iqGv0L2AQ8zeRs6U2Nsy4FDgK7gYeBXcCmkfbzL4AHRsh5P5MH\n1N1MnoFPbZx3K5OyP8rkRNvbBt7+fUxOGL7C5IHmWuB0YPv0/5ttwGmN8/YyObu+a/r19ZZ5h93+\nLPCulv+GfmRX6lA3J/ck/Y7Flzpk8aUOWXypQxZf6pDFlzpk8aUOWXypQ/8PzlW1q8n+n1cAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532802410>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "B = pony2(I)\n",
    "imshow(B, interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0 0 0\n",
      "0 0 2\n",
      "0 1 0\n",
      "0 1 2\n",
      "0 2 0\n",
      "0 2 2\n",
      "0 3 0\n",
      "0 3 2\n",
      "0 4 0\n",
      "0 4 2\n",
      "0 5 0\n",
      "0 5 2\n",
      "0 6 0\n",
      "0 6 2\n",
      "0 7 0\n",
      "0 7 2\n",
      "0 0 1\n",
      "0 0 3\n",
      "0 1 1\n",
      "0 1 3\n",
      "0 2 1\n",
      "0 2 3\n",
      "0 3 1\n",
      "0 3 3\n",
      "0 4 1\n",
      "0 4 3\n",
      "0 5 1\n",
      "0 5 3\n",
      "0 6 1\n",
      "0 6 3\n",
      "0 7 1\n",
      "0 7 3\n",
      "1 0 0\n",
      "1 0 2\n",
      "1 1 0\n",
      "1 1 2\n",
      "1 2 0\n",
      "1 2 2\n",
      "1 3 0\n",
      "1 3 2\n",
      "1 4 0\n",
      "1 4 2\n",
      "1 5 0\n",
      "1 5 2\n",
      "1 6 0\n",
      "1 6 2\n",
      "1 7 0\n",
      "1 7 2\n",
      "1 0 1\n",
      "1 0 3\n",
      "1 1 1\n",
      "1 1 3\n",
      "1 2 1\n",
      "1 2 3\n",
      "1 3 1\n",
      "1 3 3\n",
      "1 4 1\n",
      "1 4 3\n",
      "1 5 1\n",
      "1 5 3\n",
      "1 6 1\n",
      "1 6 3\n",
      "1 7 1\n",
      "1 7 3\n",
      "2 0 0\n",
      "2 0 2\n",
      "2 1 0\n",
      "2 1 2\n",
      "2 2 0\n",
      "2 2 2\n",
      "2 3 0\n",
      "2 3 2\n",
      "2 4 0\n",
      "2 4 2\n",
      "2 5 0\n",
      "2 5 2\n",
      "2 6 0\n",
      "2 6 2\n",
      "2 7 0\n",
      "2 7 2\n",
      "2 0 1\n",
      "2 0 3\n",
      "2 1 1\n",
      "2 1 3\n",
      "2 2 1\n",
      "2 2 3\n",
      "2 3 1\n",
      "2 3 3\n",
      "2 4 1\n",
      "2 4 3\n",
      "2 5 1\n",
      "2 5 3\n",
      "2 6 1\n",
      "2 6 3\n",
      "2 7 1\n",
      "2 7 3\n",
      "3 0 0\n",
      "3 0 2\n",
      "3 1 0\n",
      "3 1 2\n",
      "3 2 0\n",
      "3 2 2\n",
      "3 3 0\n",
      "3 3 2\n",
      "3 4 0\n",
      "3 4 2\n",
      "3 5 0\n",
      "3 5 2\n",
      "3 6 0\n",
      "3 6 2\n",
      "3 7 0\n",
      "3 7 2\n",
      "3 0 1\n",
      "3 0 3\n",
      "3 1 1\n",
      "3 1 3\n",
      "3 2 1\n",
      "3 2 3\n",
      "3 3 1\n",
      "3 3 3\n",
      "3 4 1\n",
      "3 4 3\n",
      "3 5 1\n",
      "3 5 3\n",
      "3 6 1\n",
      "3 6 3\n",
      "3 7 1\n",
      "3 7 3\n",
      "4 0 0\n",
      "4 0 2\n",
      "4 1 0\n",
      "4 1 2\n",
      "4 2 0\n",
      "4 2 2\n",
      "4 3 0\n",
      "4 3 2\n",
      "4 4 0\n",
      "4 4 2\n",
      "4 5 0\n",
      "4 5 2\n",
      "4 6 0\n",
      "4 6 2\n",
      "4 7 0\n",
      "4 7 2\n",
      "4 0 1\n",
      "4 0 3\n",
      "4 1 1\n",
      "4 1 3\n",
      "4 2 1\n",
      "4 2 3\n",
      "4 3 1\n",
      "4 3 3\n",
      "4 4 1\n",
      "4 4 3\n",
      "4 5 1\n",
      "4 5 3\n",
      "4 6 1\n",
      "4 6 3\n",
      "4 7 1\n",
      "4 7 3\n",
      "5 0 0\n",
      "5 0 2\n",
      "5 1 0\n",
      "5 1 2\n",
      "5 2 0\n",
      "5 2 2\n",
      "5 3 0\n",
      "5 3 2\n",
      "5 4 0\n",
      "5 4 2\n",
      "5 5 0\n",
      "5 5 2\n",
      "5 6 0\n",
      "5 6 2\n",
      "5 7 0\n",
      "5 7 2\n",
      "5 0 1\n",
      "5 0 3\n",
      "5 1 1\n",
      "5 1 3\n",
      "5 2 1\n",
      "5 2 3\n",
      "5 3 1\n",
      "5 3 3\n",
      "5 4 1\n",
      "5 4 3\n",
      "5 5 1\n",
      "5 5 3\n",
      "5 6 1\n",
      "5 6 3\n",
      "5 7 1\n",
      "5 7 3\n",
      "6 0 0\n",
      "6 0 2\n",
      "6 1 0\n",
      "6 1 2\n",
      "6 2 0\n",
      "6 2 2\n",
      "6 3 0\n",
      "6 3 2\n",
      "6 4 0\n",
      "6 4 2\n",
      "6 5 0\n",
      "6 5 2\n",
      "6 6 0\n",
      "6 6 2\n",
      "6 7 0\n",
      "6 7 2\n",
      "6 0 1\n",
      "6 0 3\n",
      "6 1 1\n",
      "6 1 3\n",
      "6 2 1\n",
      "6 2 3\n",
      "6 3 1\n",
      "6 3 3\n",
      "6 4 1\n",
      "6 4 3\n",
      "6 5 1\n",
      "6 5 3\n",
      "6 6 1\n",
      "6 6 3\n",
      "6 7 1\n",
      "6 7 3\n",
      "7 0 0\n",
      "7 0 2\n",
      "7 1 0\n",
      "7 1 2\n",
      "7 2 0\n",
      "7 2 2\n",
      "7 3 0\n",
      "7 3 2\n",
      "7 4 0\n",
      "7 4 2\n",
      "7 5 0\n",
      "7 5 2\n",
      "7 6 0\n",
      "7 6 2\n",
      "7 7 0\n",
      "7 7 2\n",
      "7 0 1\n",
      "7 0 3\n",
      "7 1 1\n",
      "7 1 3\n",
      "7 2 1\n",
      "7 2 3\n",
      "7 3 1\n",
      "7 3 3\n",
      "7 4 1\n",
      "7 4 3\n",
      "7 5 1\n",
      "7 5 3\n",
      "7 6 1\n",
      "7 6 3\n",
      "7 7 1\n",
      "7 7 3\n"
     ]
    }
   ],
   "source": [
    "I = np.arange(16*16).reshape(8, 8, 4)\n",
    "A = pony(I, 2)\n",
    "B = pony2(I)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f95324132d0>"
      ]
     },
     "execution_count": 79,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD3VJREFUeJzt3WuMnNV9x/Hvz7vrC7jchICCZRtKIAiVUFflUpJ2GxLh\nJgqOqlaCNE0gEq+SgJIqgsALNm+qJFVCaUle0CQWbaGR4iDFaUPjWs6GtARyMeZmY1yZYsMWV1Gb\noJQEfPn3xcza683MsHvmOWd29vw+0kizM3Pmfzzj3zwzz5l5/ooIzKwuSwY9ATMrz8E3q5CDb1Yh\nB9+sQg6+WYUcfLMKjeYuIMnrhWYDEhHqdHn24Lfc2eXySWC8y3UrEmv1GvcQ8IcF6z0I/FHD9U7o\ncd0/AO9vuF6vcX8L3FSo1l8DNyeMS633GeDWgvUm2qdOOmZ3DrqP6+utvqT1kp6V9Jykbo+SmS0w\nycGXtAS4B7gGuBi4XtKbm5qYmeXTzxb/MmBPRLwQEQeBrwIb5ncXa/son+L8wvUuKlzvksL11hWs\ndXnBWgBXFa43XrRaP8E/B9g/4+8X25fNw9o+yqd4U+F6iz34v12wVungv7VwvfGi1bycZ1ahfvbq\nvwSsnvH3qvZlHUzOOL+W8lt6sxpMcnzWulPqz3IljQC7gauB/wJ+AFwfEbtm3S66L+f1kmPJZDGM\n67Wcl6NeyrhheSyHZVz6cl7j6/gRcVjSR4AttD4yfHl26M1sYerrCzwR8S/AhQ3NxcwK8c49swo5\n+GYVcvDNKlToRzpjhcbUMC71KSs5rvQch+W5Wzi8xTerkINvViEH36xCDr5ZhRx8swo5+GYVcvDN\nKuTgm1XIwTerkINvViEH36xCDr5ZhZIPvTXnAlLAkfkPXJ5YcFniuNR6HjfYWkM1LvEQWqn1PtP9\n0Fve4ptVyME3q1A/LbRWSdom6RlJT0nq1tHQzBaYfg7EcQj4eETskLQS+LGkLRHxbENzM7NMkrf4\nEfFyROxon/85sIt5t9Ays0Fo5DO+pLXApcBjTdyfmeXV9zH32m/zNwG3tLf8HUzMOD9O6QaBZlXY\nN9k6zUFf6/iSRoF/Ah6KiLu73Mbr+LWNG4Y5DmTc4lnH/wqws1vozWxh6mc57yrgT4G3S3pc0nZJ\n65ubmpnl0k/TzH8HRhqci5kV4m/umVXIwTerUJkWWqMJezMXezelxTxuGOY4TOMypNRbfLMKOfhm\nFXLwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q1CZ3nl3JtQYmuOo\nedxAawEsT/w/vPxg0rAliePGlr+WNG7Z8teTxr2y4tfdO8/MjnHwzSrUd/AlLWkfaHNzExMys/ya\n2OLfAuxs4H7MrJC+gi9pFfAu4EvNTMfMSuh3i38X8Akg79KAmTUq+TB+kt4NHGi3yR4Huh9Rc3Li\n2Pm1462TmTXq0MOPcOjhR+Z02+R1fEl/AbwfOASsAH4NeDAiPjDrdl7Hr22c1/E7WhTr+BFxe0Ss\njojzgOuAbbNDb2YLk9fxzSrUyKH6I+K7wHebuC8zy89bfLMKOfhmFSrTO8/91+oal1wrce988rjD\nScNGxg6llUscNzKaNs9evMU3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q5CDb1Yh\nB9+sQg6+WYUcfLMKFfl13um375v3mGWkHZ9sKWnHJ1vs9ZYm1luWUK9krVa94XgOStf7yx7XeYtv\nViEH36xC/XbSOVnS1yTtkvSMpMubmpiZ5dPvZ/y7gW9FxJ9IGgVOaGBOZpZZP510TgLeFhE3AETE\nIeCVhuZlZhn181b/XOAnkja222TfK2lFUxMzs3z6eas/CqwDPhwRP5L0V8BtwJ2zb/h/E3cdPT82\nfgVLx6/so6yZdbJv8nn2T/7nnG7bT++8M4Hvt1toIemtwK0R8Z5Zt4vTj7ww7/tf7GusXsdvplar\n3nA8B8XX8ZdMZOmddwDYL+mC9kVXAztT78/Myul3r/7NwP2SxoC9wI39T8nMcusr+BHxBPA7Dc3F\nzArxN/fMKuTgm1WoyK/zRnVw/mNI6zM2yvxrDVe9xL5tBeuVrNUat7if8xHcO8/MGuDgm1XIwTer\nkINvViEH36xCDr5ZhRx8swo5+GYVcvDNKuTgm1XIwTerkINvViEH36xCycfcm3MBKe6IT857XPrx\n10of783Hlxtkrda4xf6cp9U7Xy83f8w9MxteDr5ZhfrtnfcxSU9LelLS/ZKWNjUxM8snOfiSzgY+\nCqyLiEtoHc3nuqYmZmb59HvorRHgRElHaDXMnOp/SmaWWz8NNaaAzwH7gJeAn0bE1qYmZmb59NMt\n9xRgA7AG+BmwSdL7IuKB2bd9eOJ7R8+vGV/NmvE1qWXNrItHJ1/jscm5LRn20zvvj4FrIuKm9t9/\nBlweER+ZdTuv4zdaz+v4vzpusT/nC2sdfx9whaTlkkSrd96uPu7PzArp5zP+D4BNwOPAE4CAexua\nl5ll1G/vvE8Bn2poLmZWiL+5Z1YhB9+sQkV656X0UhuGHnGDqDcMfeIWc1/AQdRz7zwza4SDb1Yh\nB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKFemdtz9Om/e4pZF4\nbLnDicd7ey3tV2/L0sqhX6aNSzz8Wvq4lHmWrDVM9QrPc8lNuHeemR3j4JtV6A2DL+nLkg5IenLG\nZadK2iJpt6RvSzo57zTNrElz2eJvBK6ZddltwNaIuBDYBsz/wPlmNjBvGPyI+Dfgf2ddvAG4r33+\nPuC9Dc/LzDJK/Yx/RkQcAIiIl4EzmpuSmeXW1ME2e64Jfn7i1aPnrxwf48rxsYbKmtm0yd2t01zM\naR1f0hrgmxFxSfvvXcB4RByQdBbwnYi4qMtYr+N34HX8hmoNU70hXMdX+zRtM3BD+/wHgW+kTc3M\nBmEuy3kPAI8AF0jaJ+lG4NPAOyXtptUs89N5p2lmTXrDz/gR8b4uV72j4bmYWSH+5p5ZhRx8swoV\n6Z2X0jMsuUfcobQ+Y6Np5VDquLQ2aiQ+LGXHpdYq/ZgMy3PQfOs8b/HNauTgm1XIwTerkINvViEH\n36xCDr5ZhRx8swo5+GYVcvDNKuTgm1XIwTerkINvViEH36xCRX6dd+YXXpn/oCE5HlrxY+cNw+Pi\nx6Sz0vPswVt8swo5+GYVSu2d91lJuyTtkPR1SSflnaaZNSm1d94W4OKIuBTYg3vnmQ2VpN55EbE1\nIo60/3wUWJVhbmaWSROf8T8EPNTA/ZhZIX0t50m6AzgYEQ/0ut3EPx87P/4mGL+gn6pm1snkqzD5\ni7ndNql3XvuyG4CbgLdHRNeVRklx5J65TeY4Q7I263X8DvyYdFZ4ntrTvXfeXLf4x/XOk7Qe+ATw\ne71Cb2YLU2rvvL8BVgL/Kmm7pC9mnqeZNSi1d97GDHMxs0L8zT2zCjn4ZhUq8uu8pP5yw9LXzP3e\nBlsLhuMx6Wdc6r+vB2/xzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3\nq5CDb1YhB9+sQkV+ncfn5z8k9YBeBxPH/TLxeGi/SKz3atow5ngsxYGOG4Y51jCuF2/xzSqU1EJr\nxnV/LumIpNPyTM/MckhtoYWkVcA7gReanpSZ5ZXUQqvtLlqH2DazIZP0GV/StcD+iHiq4fmYWQHz\n3qsvaQVwO623+UcvbmxGZpZdynLebwBrgSckiVan3B9Luiwi/rvTgImfHTs/vgzGlydUNbOe9gLP\nz/G2826hFRFPA2cdvUJ6HlgXEZ32AwAwcfIcq5hZsvPap2nf6XHb1BZaMwV+q282VFJbaM28/rxe\n15vZwuNv7plVyME3q5CDb1ahIr/Oi4TeX5HYZ+xg4rjXE/uTvZ42bGjatpVsnZf6WC72camPZy/e\n4ptVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVaKDBn0w8Jn2qhw+Xrff9suX4YeF6\nOwrWeqZgLYDdhevtLVxvsMFP/Q5jou8t8uD/qHC9xRz85wrXqyr4ZjYYZVpo/ea6zpf/xxScf3bn\n6xLfDSzpMU77pliyunO9kcRfQoz1eBcxMjXF2Nmd6x1JK9d7LlNTrOhSL/UVfqTHdUunpljZod5Y\nYq2lPa5bPjXFKV3+bSsS6/X6pHnC1BSnd6mX+gm1V5e2E6emOLNLveQ3xtu3d71KEZF6t3MiKW8B\nM+sqIjoeFi978M1s4fFnfLMKOfhmFRpI8CWtl/SspOck3Zq51ipJ2yQ9I+kpSTfnrDej7hJJ2yVt\nLlDrZElfk7Sr/e+8PHO9j0l6WtKTku6X1Gu/XMr9/0qHZkmnStoiabekb0tqrFtDl3qfbT+eOyR9\nXdJJOevNuK5IB+riwZe0BLiHVgfei4HrJb05Y8lDwMcj4mLgSuDDmetNuwXYWaAOwN3AtyLiIuAt\nwK5chSSdDXyUVhOVS2itDF3XcJlOHZpvA7ZGxIXANuCTmettAS6OiEuBPQXqFe1APYgt/mXAnoh4\nISIOAl8FNuQqFhEvR8SO9vmf0wrFObnqwdEn8F3Al3LWadc6CXhbRGwEiIhDEfFK5rIjwImSRoET\ngKkm77xLh+YNwH3t8/cB781ZLyK2RsT0quujtFrFZavXVqwD9SCCfw6wf8bfL5I5iNMkrQUuBR7L\nXGr6CSyxZHIu8BNJG9sfLe5tNzbNIiKmgM8B+4CXgJ9GxNZc9WY4IyIOtOfwMnBGgZrTPgQ8lLNA\n6Q7U1ezck7QS2ATc0t7y56rzbuBA+13G0Z6DGY0C64AvRMQ64FVab4uzkHQKra3vGuBsYKWknt2W\nMimyDi3pDuBgRDyQscZ0B+o7Z16cqx4MJvgvAatn/L2qfVk27bekm4C/j4hv5KwFXAVcK2kv8I/A\nH0j6u4z1XqS1pZj+qv4mWi8EubwD2BsR/xMRh4EHgd/NWG/aAUlnAkg6C+jYmblJkm6g9ZEt9wvb\nzA7Uz3OsA3W2dzWDCP4PgfMlrWnvDb4OyL3n+yvAzoi4O3MdIuL2iFjd7il4HbAtIj6Qsd4BYL+k\nC9oXXU3enYr7gCskLW+3Sb+aPDsTZ79b2gzc0D7/QaDpF/Dj6klaT+vj2rURkeN3pMd1oI6IsyLi\nvIg4l9aL+W91azvfiIgofgLW0/rl4x7gtsy1rgIO0/ox2ePAdmB9oX/n7wObC9R5C60X1B20tsAn\nZ653J62wP0lrR9tYw/f/AK0dhq/ReqG5ETgV2Nr+f7MFOCVzvT209q5vb5++mLPerOv3AqflfA79\nlV2zClWzc8/MjnHwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMK/T90LNcW4MvK8QAAAABJRU5E\nrkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f9532634250>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imshow(A[:, :, 0], interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7f95324e0790>"
      ]
     },
     "execution_count": 80,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAP4AAAD7CAYAAABKWyniAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD3VJREFUeJzt3WuMnNV9x/Hvz7vrC7jchICCZRtKIAiVUFflUpJ2GxLh\nJgqOqlaCNE0gEq+SgJIqgsALNm+qJFVCaUle0CQWbaGR4iDFaUPjWs6GtARyMeZmY1yZYsMWV1Gb\noJQEfPn3xcza683MsHvmOWd29vw+0kizM3Pmfzzj3zwzz5l5/ooIzKwuSwY9ATMrz8E3q5CDb1Yh\nB9+sQg6+WYUcfLMKjeYuIMnrhWYDEhHqdHn24Lfc2eXySWC8y3UrEmv1GvcQ8IcF6z0I/FHD9U7o\ncd0/AO9vuF6vcX8L3FSo1l8DNyeMS633GeDWgvUm2qdOOmZ3DrqP6+utvqT1kp6V9Jykbo+SmS0w\nycGXtAS4B7gGuBi4XtKbm5qYmeXTzxb/MmBPRLwQEQeBrwIb5ncXa/son+L8wvUuKlzvksL11hWs\ndXnBWgBXFa43XrRaP8E/B9g/4+8X25fNw9o+yqd4U+F6iz34v12wVungv7VwvfGi1bycZ1ahfvbq\nvwSsnvH3qvZlHUzOOL+W8lt6sxpMcnzWulPqz3IljQC7gauB/wJ+AFwfEbtm3S66L+f1kmPJZDGM\n67Wcl6NeyrhheSyHZVz6cl7j6/gRcVjSR4AttD4yfHl26M1sYerrCzwR8S/AhQ3NxcwK8c49swo5\n+GYVcvDNKlToRzpjhcbUMC71KSs5rvQch+W5Wzi8xTerkINvViEH36xCDr5ZhRx8swo5+GYVcvDN\nKuTgm1XIwTerkINvViEH36xCDr5ZhZIPvTXnAlLAkfkPXJ5YcFniuNR6HjfYWkM1LvEQWqn1PtP9\n0Fve4ptVyME3q1A/LbRWSdom6RlJT0nq1tHQzBaYfg7EcQj4eETskLQS+LGkLRHxbENzM7NMkrf4\nEfFyROxon/85sIt5t9Ays0Fo5DO+pLXApcBjTdyfmeXV9zH32m/zNwG3tLf8HUzMOD9O6QaBZlXY\nN9k6zUFf6/iSRoF/Ah6KiLu73Mbr+LWNG4Y5DmTc4lnH/wqws1vozWxh6mc57yrgT4G3S3pc0nZJ\n65ubmpnl0k/TzH8HRhqci5kV4m/umVXIwTerUJkWWqMJezMXezelxTxuGOY4TOMypNRbfLMKOfhm\nFXLwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q1CZ3nl3JtQYmuOo\nedxAawEsT/w/vPxg0rAliePGlr+WNG7Z8teTxr2y4tfdO8/MjnHwzSrUd/AlLWkfaHNzExMys/ya\n2OLfAuxs4H7MrJC+gi9pFfAu4EvNTMfMSuh3i38X8Akg79KAmTUq+TB+kt4NHGi3yR4Huh9Rc3Li\n2Pm1462TmTXq0MOPcOjhR+Z02+R1fEl/AbwfOASsAH4NeDAiPjDrdl7Hr22c1/E7WhTr+BFxe0Ss\njojzgOuAbbNDb2YLk9fxzSrUyKH6I+K7wHebuC8zy89bfLMKOfhmFSrTO8/91+oal1wrce988rjD\nScNGxg6llUscNzKaNs9evMU3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q5CDb1Yh\nB9+sQg6+WYUcfLMKFfl13um375v3mGWkHZ9sKWnHJ1vs9ZYm1luWUK9krVa94XgOStf7yx7XeYtv\nViEH36xC/XbSOVnS1yTtkvSMpMubmpiZ5dPvZ/y7gW9FxJ9IGgVOaGBOZpZZP510TgLeFhE3AETE\nIeCVhuZlZhn181b/XOAnkja222TfK2lFUxMzs3z6eas/CqwDPhwRP5L0V8BtwJ2zb/h/E3cdPT82\nfgVLx6/so6yZdbJv8nn2T/7nnG7bT++8M4Hvt1toIemtwK0R8Z5Zt4vTj7ww7/tf7GusXsdvplar\n3nA8B8XX8ZdMZOmddwDYL+mC9kVXAztT78/Myul3r/7NwP2SxoC9wI39T8nMcusr+BHxBPA7Dc3F\nzArxN/fMKuTgm1WoyK/zRnVw/mNI6zM2yvxrDVe9xL5tBeuVrNUat7if8xHcO8/MGuDgm1XIwTer\nkINvViEH36xCDr5ZhRx8swo5+GYVcvDNKuTgm1XIwTerkINvViEH36xCycfcm3MBKe6IT857XPrx\n10of783Hlxtkrda4xf6cp9U7Xy83f8w9MxteDr5ZhfrtnfcxSU9LelLS/ZKWNjUxM8snOfiSzgY+\nCqyLiEtoHc3nuqYmZmb59HvorRHgRElHaDXMnOp/SmaWWz8NNaaAzwH7gJeAn0bE1qYmZmb59NMt\n9xRgA7AG+BmwSdL7IuKB2bd9eOJ7R8+vGV/NmvE1qWXNrItHJ1/jscm5LRn20zvvj4FrIuKm9t9/\nBlweER+ZdTuv4zdaz+v4vzpusT/nC2sdfx9whaTlkkSrd96uPu7PzArp5zP+D4BNwOPAE4CAexua\nl5ll1G/vvE8Bn2poLmZWiL+5Z1YhB9+sQkV656X0UhuGHnGDqDcMfeIWc1/AQdRz7zwza4SDb1Yh\nB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKFemdtz9Om/e4pZF4\nbLnDicd7ey3tV2/L0sqhX6aNSzz8Wvq4lHmWrDVM9QrPc8lNuHeemR3j4JtV6A2DL+nLkg5IenLG\nZadK2iJpt6RvSzo57zTNrElz2eJvBK6ZddltwNaIuBDYBsz/wPlmNjBvGPyI+Dfgf2ddvAG4r33+\nPuC9Dc/LzDJK/Yx/RkQcAIiIl4EzmpuSmeXW1ME2e64Jfn7i1aPnrxwf48rxsYbKmtm0yd2t01zM\naR1f0hrgmxFxSfvvXcB4RByQdBbwnYi4qMtYr+N34HX8hmoNU70hXMdX+zRtM3BD+/wHgW+kTc3M\nBmEuy3kPAI8AF0jaJ+lG4NPAOyXtptUs89N5p2lmTXrDz/gR8b4uV72j4bmYWSH+5p5ZhRx8swoV\n6Z2X0jMsuUfcobQ+Y6Np5VDquLQ2aiQ+LGXHpdYq/ZgMy3PQfOs8b/HNauTgm1XIwTerkINvViEH\n36xCDr5ZhRx8swo5+GYVcvDNKuTgm1XIwTerkINvViEH36xCRX6dd+YXXpn/oCE5HlrxY+cNw+Pi\nx6Sz0vPswVt8swo5+GYVSu2d91lJuyTtkPR1SSflnaaZNSm1d94W4OKIuBTYg3vnmQ2VpN55EbE1\nIo60/3wUWJVhbmaWSROf8T8EPNTA/ZhZIX0t50m6AzgYEQ/0ut3EPx87P/4mGL+gn6pm1snkqzD5\ni7ndNql3XvuyG4CbgLdHRNeVRklx5J65TeY4Q7I263X8DvyYdFZ4ntrTvXfeXLf4x/XOk7Qe+ATw\ne71Cb2YLU2rvvL8BVgL/Kmm7pC9mnqeZNSi1d97GDHMxs0L8zT2zCjn4ZhUq8uu8pP5yw9LXzP3e\nBlsLhuMx6Wdc6r+vB2/xzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVyME3\nq5CDb1YhB9+sQkV+ncfn5z8k9YBeBxPH/TLxeGi/SKz3atow5ngsxYGOG4Y51jCuF2/xzSqU1EJr\nxnV/LumIpNPyTM/MckhtoYWkVcA7gReanpSZ5ZXUQqvtLlqH2DazIZP0GV/StcD+iHiq4fmYWQHz\n3qsvaQVwO623+UcvbmxGZpZdynLebwBrgSckiVan3B9Luiwi/rvTgImfHTs/vgzGlydUNbOe9gLP\nz/G2826hFRFPA2cdvUJ6HlgXEZ32AwAwcfIcq5hZsvPap2nf6XHb1BZaMwV+q282VFJbaM28/rxe\n15vZwuNv7plVyME3q5CDb1ahIr/Oi4TeX5HYZ+xg4rjXE/uTvZ42bGjatpVsnZf6WC72camPZy/e\n4ptVyME3q5CDb1YhB9+sQg6+WYUcfLMKOfhmFXLwzSrk4JtVaKDBn0w8Jn2qhw+Xrff9suX4YeF6\nOwrWeqZgLYDdhevtLVxvsMFP/Q5jou8t8uD/qHC9xRz85wrXqyr4ZjYYZVpo/ea6zpf/xxScf3bn\n6xLfDSzpMU77pliyunO9kcRfQoz1eBcxMjXF2Nmd6x1JK9d7LlNTrOhSL/UVfqTHdUunpljZod5Y\nYq2lPa5bPjXFKV3+bSsS6/X6pHnC1BSnd6mX+gm1V5e2E6emOLNLveQ3xtu3d71KEZF6t3MiKW8B\nM+sqIjoeFi978M1s4fFnfLMKOfhmFRpI8CWtl/SspOck3Zq51ipJ2yQ9I+kpSTfnrDej7hJJ2yVt\nLlDrZElfk7Sr/e+8PHO9j0l6WtKTku6X1Gu/XMr9/0qHZkmnStoiabekb0tqrFtDl3qfbT+eOyR9\nXdJJOevNuK5IB+riwZe0BLiHVgfei4HrJb05Y8lDwMcj4mLgSuDDmetNuwXYWaAOwN3AtyLiIuAt\nwK5chSSdDXyUVhOVS2itDF3XcJlOHZpvA7ZGxIXANuCTmettAS6OiEuBPQXqFe1APYgt/mXAnoh4\nISIOAl8FNuQqFhEvR8SO9vmf0wrFObnqwdEn8F3Al3LWadc6CXhbRGwEiIhDEfFK5rIjwImSRoET\ngKkm77xLh+YNwH3t8/cB781ZLyK2RsT0quujtFrFZavXVqwD9SCCfw6wf8bfL5I5iNMkrQUuBR7L\nXGr6CSyxZHIu8BNJG9sfLe5tNzbNIiKmgM8B+4CXgJ9GxNZc9WY4IyIOtOfwMnBGgZrTPgQ8lLNA\n6Q7U1ezck7QS2ATc0t7y56rzbuBA+13G0Z6DGY0C64AvRMQ64FVab4uzkHQKra3vGuBsYKWknt2W\nMimyDi3pDuBgRDyQscZ0B+o7Z16cqx4MJvgvAatn/L2qfVk27bekm4C/j4hv5KwFXAVcK2kv8I/A\nH0j6u4z1XqS1pZj+qv4mWi8EubwD2BsR/xMRh4EHgd/NWG/aAUlnAkg6C+jYmblJkm6g9ZEt9wvb\nzA7Uz3OsA3W2dzWDCP4PgfMlrWnvDb4OyL3n+yvAzoi4O3MdIuL2iFjd7il4HbAtIj6Qsd4BYL+k\nC9oXXU3enYr7gCskLW+3Sb+aPDsTZ79b2gzc0D7/QaDpF/Dj6klaT+vj2rURkeN3pMd1oI6IsyLi\nvIg4l9aL+W91azvfiIgofgLW0/rl4x7gtsy1rgIO0/ox2ePAdmB9oX/n7wObC9R5C60X1B20tsAn\nZ653J62wP0lrR9tYw/f/AK0dhq/ReqG5ETgV2Nr+f7MFOCVzvT209q5vb5++mLPerOv3AqflfA79\nlV2zClWzc8/MjnHwzSrk4JtVyME3q5CDb1YhB9+sQg6+WYUcfLMK/T90LNcW4MvK8QAAAABJRU5E\nrkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f95328e53d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "imshow(B, interpolation=\"none\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[[   0.],\n",
       "        [   2.],\n",
       "        [   4.],\n",
       "        [   6.],\n",
       "        [   8.],\n",
       "        [  10.],\n",
       "        [  12.],\n",
       "        [  14.],\n",
       "        [  16.],\n",
       "        [  18.],\n",
       "        [  20.],\n",
       "        [  22.],\n",
       "        [  24.],\n",
       "        [  26.],\n",
       "        [  28.],\n",
       "        [  30.]],\n",
       "\n",
       "       [[   1.],\n",
       "        [   3.],\n",
       "        [   5.],\n",
       "        [   7.],\n",
       "        [   9.],\n",
       "        [  11.],\n",
       "        [  13.],\n",
       "        [  15.],\n",
       "        [  17.],\n",
       "        [  19.],\n",
       "        [  21.],\n",
       "        [  23.],\n",
       "        [  25.],\n",
       "        [  27.],\n",
       "        [  29.],\n",
       "        [  31.]],\n",
       "\n",
       "       [[  32.],\n",
       "        [  34.],\n",
       "        [  36.],\n",
       "        [  38.],\n",
       "        [  40.],\n",
       "        [  42.],\n",
       "        [  44.],\n",
       "        [  46.],\n",
       "        [  48.],\n",
       "        [  50.],\n",
       "        [  52.],\n",
       "        [  54.],\n",
       "        [  56.],\n",
       "        [  58.],\n",
       "        [  60.],\n",
       "        [  62.]],\n",
       "\n",
       "       [[  33.],\n",
       "        [  35.],\n",
       "        [  37.],\n",
       "        [  39.],\n",
       "        [  41.],\n",
       "        [  43.],\n",
       "        [  45.],\n",
       "        [  47.],\n",
       "        [  49.],\n",
       "        [  51.],\n",
       "        [  53.],\n",
       "        [  55.],\n",
       "        [  57.],\n",
       "        [  59.],\n",
       "        [  61.],\n",
       "        [  63.]],\n",
       "\n",
       "       [[  64.],\n",
       "        [  66.],\n",
       "        [  68.],\n",
       "        [  70.],\n",
       "        [  72.],\n",
       "        [  74.],\n",
       "        [  76.],\n",
       "        [  78.],\n",
       "        [  80.],\n",
       "        [  82.],\n",
       "        [  84.],\n",
       "        [  86.],\n",
       "        [  88.],\n",
       "        [  90.],\n",
       "        [  92.],\n",
       "        [  94.]],\n",
       "\n",
       "       [[  65.],\n",
       "        [  67.],\n",
       "        [  69.],\n",
       "        [  71.],\n",
       "        [  73.],\n",
       "        [  75.],\n",
       "        [  77.],\n",
       "        [  79.],\n",
       "        [  81.],\n",
       "        [  83.],\n",
       "        [  85.],\n",
       "        [  87.],\n",
       "        [  89.],\n",
       "        [  91.],\n",
       "        [  93.],\n",
       "        [  95.]],\n",
       "\n",
       "       [[  96.],\n",
       "        [  98.],\n",
       "        [ 100.],\n",
       "        [ 102.],\n",
       "        [ 104.],\n",
       "        [ 106.],\n",
       "        [ 108.],\n",
       "        [ 110.],\n",
       "        [ 112.],\n",
       "        [ 114.],\n",
       "        [ 116.],\n",
       "        [ 118.],\n",
       "        [ 120.],\n",
       "        [ 122.],\n",
       "        [ 124.],\n",
       "        [ 126.]],\n",
       "\n",
       "       [[  97.],\n",
       "        [  99.],\n",
       "        [ 101.],\n",
       "        [ 103.],\n",
       "        [ 105.],\n",
       "        [ 107.],\n",
       "        [ 109.],\n",
       "        [ 111.],\n",
       "        [ 113.],\n",
       "        [ 115.],\n",
       "        [ 117.],\n",
       "        [ 119.],\n",
       "        [ 121.],\n",
       "        [ 123.],\n",
       "        [ 125.],\n",
       "        [ 127.]],\n",
       "\n",
       "       [[ 128.],\n",
       "        [ 130.],\n",
       "        [ 132.],\n",
       "        [ 134.],\n",
       "        [ 136.],\n",
       "        [ 138.],\n",
       "        [ 140.],\n",
       "        [ 142.],\n",
       "        [ 144.],\n",
       "        [ 146.],\n",
       "        [ 148.],\n",
       "        [ 150.],\n",
       "        [ 152.],\n",
       "        [ 154.],\n",
       "        [ 156.],\n",
       "        [ 158.]],\n",
       "\n",
       "       [[ 129.],\n",
       "        [ 131.],\n",
       "        [ 133.],\n",
       "        [ 135.],\n",
       "        [ 137.],\n",
       "        [ 139.],\n",
       "        [ 141.],\n",
       "        [ 143.],\n",
       "        [ 145.],\n",
       "        [ 147.],\n",
       "        [ 149.],\n",
       "        [ 151.],\n",
       "        [ 153.],\n",
       "        [ 155.],\n",
       "        [ 157.],\n",
       "        [ 159.]],\n",
       "\n",
       "       [[ 160.],\n",
       "        [ 162.],\n",
       "        [ 164.],\n",
       "        [ 166.],\n",
       "        [ 168.],\n",
       "        [ 170.],\n",
       "        [ 172.],\n",
       "        [ 174.],\n",
       "        [ 176.],\n",
       "        [ 178.],\n",
       "        [ 180.],\n",
       "        [ 182.],\n",
       "        [ 184.],\n",
       "        [ 186.],\n",
       "        [ 188.],\n",
       "        [ 190.]],\n",
       "\n",
       "       [[ 161.],\n",
       "        [ 163.],\n",
       "        [ 165.],\n",
       "        [ 167.],\n",
       "        [ 169.],\n",
       "        [ 171.],\n",
       "        [ 173.],\n",
       "        [ 175.],\n",
       "        [ 177.],\n",
       "        [ 179.],\n",
       "        [ 181.],\n",
       "        [ 183.],\n",
       "        [ 185.],\n",
       "        [ 187.],\n",
       "        [ 189.],\n",
       "        [ 191.]],\n",
       "\n",
       "       [[ 192.],\n",
       "        [ 194.],\n",
       "        [ 196.],\n",
       "        [ 198.],\n",
       "        [ 200.],\n",
       "        [ 202.],\n",
       "        [ 204.],\n",
       "        [ 206.],\n",
       "        [ 208.],\n",
       "        [ 210.],\n",
       "        [ 212.],\n",
       "        [ 214.],\n",
       "        [ 216.],\n",
       "        [ 218.],\n",
       "        [ 220.],\n",
       "        [ 222.]],\n",
       "\n",
       "       [[ 193.],\n",
       "        [ 195.],\n",
       "        [ 197.],\n",
       "        [ 199.],\n",
       "        [ 201.],\n",
       "        [ 203.],\n",
       "        [ 205.],\n",
       "        [ 207.],\n",
       "        [ 209.],\n",
       "        [ 211.],\n",
       "        [ 213.],\n",
       "        [ 215.],\n",
       "        [ 217.],\n",
       "        [ 219.],\n",
       "        [ 221.],\n",
       "        [ 223.]],\n",
       "\n",
       "       [[ 224.],\n",
       "        [ 226.],\n",
       "        [ 228.],\n",
       "        [ 230.],\n",
       "        [ 232.],\n",
       "        [ 234.],\n",
       "        [ 236.],\n",
       "        [ 238.],\n",
       "        [ 240.],\n",
       "        [ 242.],\n",
       "        [ 244.],\n",
       "        [ 246.],\n",
       "        [ 248.],\n",
       "        [ 250.],\n",
       "        [ 252.],\n",
       "        [ 254.]],\n",
       "\n",
       "       [[ 225.],\n",
       "        [ 227.],\n",
       "        [ 229.],\n",
       "        [ 231.],\n",
       "        [ 233.],\n",
       "        [ 235.],\n",
       "        [ 237.],\n",
       "        [ 239.],\n",
       "        [ 241.],\n",
       "        [ 243.],\n",
       "        [ 245.],\n",
       "        [ 247.],\n",
       "        [ 249.],\n",
       "        [ 251.],\n",
       "        [ 253.],\n",
       "        [ 255.]]])"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "A"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True],\n",
       "       [ True,  True,  True,  True,  True,  True,  True,  True,  True,\n",
       "         True,  True,  True,  True,  True,  True,  True]], dtype=bool)"
      ]
     },
     "execution_count": 83,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "B == A[:, :, 0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def pony2(I, r):\n",
    "    a, b, c = I.shape\n",
    "    B = I.reshape(a/r, b/r, r, r).transpose(0, 1, 3, 2)\n",
    "    B = np.concatenate([B[i] for i in range(8)], axis=1)\n",
    "    B = B.transpose(1, 2, 0)\n",
    "    B = np.concatenate([B[:, :, i] for i in range(8)], axis=1)\n",
    "    return B"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [Root]",
   "language": "python",
   "name": "Python [Root]"
  },
  "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.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}