generate_hrap_index_raster.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import rasterio as rs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "!cp ./data/hrap_grid_stere.tiff ./data/hrap_grid_stere_idx.tiff"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "fl = './data/hrap_grid_stere_idx.tiff'\n",
    "with rs.open(fl,'r') as srs:\n",
    "    rast = srs.read(1)\n",
    "    profile = srs.profile"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# generate the index raster\n",
    "k,l = rast.shape\n",
    "idx = np.arange(1,(k*l)+1,dtype=np.float64)\n",
    "idx.shape =(k,l)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# update the profile\n",
    "profile['nodata'] = -9999"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/tbarnhart/miniconda3/envs/py36/lib/python3.6/site-packages/rasterio/__init__.py:160: FutureWarning: GDAL-style transforms are deprecated and will not be supported in Rasterio 1.0.\n",
      "  transform = guard_transform(transform)\n"
     ]
    }
   ],
   "source": [
    "with rs.open(fl,'r+',**profile) as dst:\n",
    "    dst.write(idx,1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "np.save('./data/hrap_grid_index.npy',idx)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x7f06f7fbf4a8>"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHk9JREFUeJzt3X2sXdV95vHvY0x4a8AYJ9TYpgbhpqFokoAFJBmNGN7G\nMFHcPyC8VImhVB5VUGjaUYM7o3EmTSpHikKcJqJcAQ1ESQyhZLAoxKEEq4oEDjZE4S0MDlB8wQVs\nDKVDErD9mz/2uvhgn1fuOWevc9bzkY7uPevss/e++977nLXXXnstRQRmZjb+ZtS9A2ZmNhwOfDOz\nQjjwzcwK4cA3MyuEA9/MrBAOfDOzQjjwzcwK4cA3MyuEA9/MrBAz694BgDlzjoiFRx9d926Y2QjY\n9PDPtkXE+6azjiVnnRnbtm/vdnvrImLJdLaXiywCf+HRR7PxJ+vr3g0zGwE6ZNa/THcd27Zv7zpz\ndMisOdPdXi6yCHwAYndNG1ZN2zWzWhU4jlgmgR81Hvzyfuk2YHIlIn9RYyWzPpkEPkV+2tqY8t/y\naCjw95RP4LumbWZDVV7m5BH4QZGftmZWk0IzJ4/AhyIPvpnVqMDMySTw21y09fUvM+u7wE06dWoV\n+OX9TsxsGHaXFy75BL6T3cyGqrzMySPwC72AYmY1KTRz8gh8KPLgm1mNCsycTAK/zLvesua7RW2s\n+aJtraLAT9us+fdhfabMKhGxu7xKZjaBX+KnrVlJ8qvU5bY/g5dP4Gf3x2BmY63AzOkq8CV9Fvhj\nqo/ER4BLgbnAGmA28BDw6Yh4U9IBwM3AScB24IKIeLbtBqLO0TLNrDiFZk7HKQ4lzQOuBBZHxAnA\nfsCFwJeBayJiEbADuCy95TJgR0QcB1yTlutC+OGHH3508eiXYW+vft026cwEDpL0FnAwsBU4Hbg4\nvX4T8HngWmBp+h7gNuAbkhSdGvDcS6dHeV0AMxs5BdbwOwZ+RDwv6SvAc8CvgB8Bm4BXI2JnWmwS\nmJe+nwdsSe/dKek14AhgW+N6JS0HlgMcPe+oIg/+9Ph42QBk1pNmoAqsZHYMfEmHU9XajwFeBb4P\nnNNk0akEavYXs086RcQEMAGw+EMnhAPfLAPF/B9GQT/rHt006ZwJPBMRLwNIuh34GDBL0sxUy58P\nvJCWnwQWAJOSZgKHAa903kx5B9/M6lRe5nQT+M8Bp0o6mKpJ5wxgI3AfcB5VT51lwB1p+bXp+f3p\n9R93bL+HEo+9mdWpwMzppg1/g6TbqLpe7gQepmqK+UdgjaQvprIb0ltuAL4taTNVzf7CjnsRFHl6\nZWY1KTRzuuqlExErgZV7FT8NnNxk2V8D5/e8J+N2AaWki19mo2jcMqcLmdxpO4YXUMbt57H6uRLR\nR2OYOV3IJPChyAY1s14UGFCDVd7xzCfw/cdsZsNUYOZkFPjltaeZWY0KzJw8Aj8Cdu2qey9sHMzo\nODyUGaUOnpZH4EORB98GwBUH61aBfyuZBH4UeXplZnUpM3PyCPxCb4Iws5oUmjl5BD4U+WlrZjUq\nMHOyCfz85rs0s3FWYubkEfgRsGtn5+Vs9Mm9aCwHAbtdw69PgZ+2RYryekZYhgLYXd7fYiaBX+YV\nczOrS/8yR9JngT+uVsojwKXAXKqh42dTjTT86Yh4U9IBwM3AScB24IKIeDatZwXVnOC7gCsjYl0q\nXwKspppP/PqIWJXKj2m2jXb7mkng4xq+mQ1XHzJH0jzgSuD4iPiVpFuphoQ/F7gmItZI+juqIL82\nfd0REcdJuhD4MnCBpOPT+34fOAr4J0m/mzbzTeAsqsmlHpS0NiIeT+9tto2WMgp81/DNbIj6lzkz\ngYMkvQUcDGwFTgcuTq/fBHyeKoyXpu8BbgO+IUmpfE1E/AZ4Js0nMjX8/OaIeBpA0hpgqaQn2myj\n7Y62JekDwC0NRccC/4vqtOQWYCHwLPCpiNiRdn411SfcG8AlEfFQ241EwG7X8M1sSPqUORHxvKSv\nUM0M+CvgR8Am4NU0/StUNfN56ft5wJb03p2SXgOOSOUPNKy68T1b9io/Jb2n1TZa6mbGqyeBDwNI\n2g94HvgBcDVwb0SsknR1ev45qgnOF6XHKVSfOKd02k6JF1Cy5DHXrRTdZ84cSRsbnk9ExASApMOp\naufHAK8C36fKwL1Nfbo0+weLNuXNurW1W76tXpt0zgB+GRH/ImkpcFoqvwlYTxX4S4Gb0zy2D0ia\nJWluRGxtvVpftM2GT7RsELLrjttT5myLiMUtXjsTeCYiXgaQdDvwMWCWpJmpBj4feCEtPwksACYl\nzQQOo5oKdqp8SuN7mpVva7ONlnoN/AuB76Xvj5wK8YjYKun9qfztU5Zk6lSjdeAXepuzWTFy647b\nv8x5DjhV0sFUTTpnABuB+4DzqHrRLAPuSMuvTc/vT6//OCJC0lrgu5K+SnXRdhHwU6qa/KLUI+d5\nqgy+OL2n1TZa6jrwJb0H+CSwotOiTcr2ObKSlgPLAY7+7fe5hm9mw9WHzImIDZJuo+oWuRN4GJgA\n/hFYI+mLqeyG9JYbgG+ni7KvUAU4EfFY6uHzeFrP5RHVp6SkK4B1VN0yb4yIx9K6PtdiGy31UsM/\nB3goIl5Mz1+caqqRNBd4KZW3OzV5W2oDmwBY/MFF4Yu2ZjZUfcqciFgJrNyr+Gn29LJpXPbXwPkt\n1vMl4EtNyu8C7mpS3nQb7fQS+BexpzkH9pyarGLfU5YrUvehU4DX2rffg9vwzWy4ysycrgI/tU+d\nBfy3huJVwK2SLqNqx5r61LqLqkvmZqpumZd23EBAuJeOmQ1LoZnTVeBHxBtU/T4by7ZTXaDYe9kA\nLu9tN8r8tDWzupSZORndaes2fDMbogIzJ6PAL+/T1sxqVGDm5BH4Ef2bUDi7GzzMLDuFDueSR+BD\n/06vcrvBw8aDh5wYP67h16XMCyg2QsqrDI65KHL8rjwC30MrmNkwFZo5eQQ+FDm/pJnVqMDMySTw\no8hPWzOrS5mZk0ng4zZ8MxuuAjMnj8DvZ7dMK9MMd8e1HhSaOXkEPhR5emV9VOA/r01TgZmTUeCX\nd3plZjUqMHPyCPygyLvezKwmhWZOHoHvG6/MbKjKzJxMAp8i29PMrEYFZk4egV/oFfOizPBYNJaR\nQjOn2xmvZgHXAydQtX79EfAkcAuwEHgW+FRE7JAkYDXVrFdvAJdExEMdN1Lg6VVRyvvfstwVmDnd\n1vBXAz+MiPMkvQc4GPgr4N6IWCXpauBqqlnUzwEWpccpwLXpa3sFnl6ZWY0KzJyOgS/pUOA/AZcA\nRMSbwJuSlgKnpcVuAtZTBf5S4OY01eEDkmZJmtt2IvNCr5ibWU0KzZxuavjHAi8Dfy/pQ8Am4Crg\nyKkQj4itkt6flp8HbGl4/2Qqax34hQ5Vug+PuW42JO6l026ZE4E/jYgNklZTNd+00iy19vkolbQc\nWA5w9BGzijz4+yivwmGD5hngWnOTTlOTwGREbEjPb6MK/BenmmokzQVealh+QcP75wMv7L3SiJgA\nJgAWHzMvSjy9Mhs8nzk35V46zUXEv0raIukDEfEkcAbweHosA1alr3ekt6wFrpC0hupi7Wtt2+/f\n3pBr+GY2RAVmTre9dP4U+E7qofM0cCkwA7hV0mXAc8D5adm7qLpkbqbqlnlpV1so8PTKzGpUYOZ0\nFfgR8TNgcZOXzmiybACX97QXEUXOPmNmNSk0c/K40xby76Xji19m48U1/BrlfvAj8w8kG03uiluf\n3CuZA5BH4Bd6emVmNSk0c/IIfMi/hm9m46XAzMko8Mv7tDWzGhWYOZkEfhQ5roWZ1aXMzMkj8IPp\nX0DxxS8z61bgGn6tptueVmB7nA2BKxLjq08Xbfs1X4ikZcD/TKv9YkTclMpPAr4FHER1Y+tVERGS\nZjfbRrt9zSTwy7xibmZ16WvmTHu+kBTeK6lucA1gk6S1KcCvpRpo8gGqwF8C3J3W2WwbLeUR+IFr\n6GY2PH3KnH7NF5KWvSciXknrvQdYImk9cGhE3J/Kbwb+gCrwW22jpTwCHxz4ZjZc/cmcfs0X0q58\nskk5bbbRkgPfzMrUfebMkbSx4flEGt4d+jdfSK/l70oegV/o2NRm1sKgL5b3ljnbIqLZ4JHQv/lC\nJtnTPDNVvj6Vz2+yPG220VIegQ+u4ZvZHsPIgz5so1/zhUhaB/yNpMPTcmcDKyLiFUmvSzoV2AB8\nBvjbhnU120ZLmQS+e+mY2TD1NXOmPV9ICva/Bh5My31h6gIu8Cfs6ZZ5d3pAFfTNttFSJoGPa/hm\nNlx9ypx+zRcSETcCNzYp30jVx3/v8u3NttFOV4Ev6VngdaoJMndGxOJWnf7b3VjQkrtlmtkwFZo5\nvdTw/3NEbGt43qrTf9MbC9qv2hdtB2qGJ28x24cDvyc93VjQcSLzAg/+0PjD1Ppt5IecKLOS2W3g\nB/AjSQFcl/qg9npjwTsCX9JyqtuFOfqwQxz4ZqNk1P9f3aTT1scj4oUU6vdI+kWbZbu6USB9aEwA\nLD5qTpR48M2sRgVmTleBHxEvpK8vSfoBcDK931jQev0E4W6ZZjYkpWZOx8CXdAgwIyJeT9+fDXyB\n1p3+m95Y0HFPCvy0NbMaFZg53dTwjwR+UPW2ZCbw3Yj4oaQH6eHGgraCcm+8GvmLX2YjyG34zUXE\n08CHmpQ37fTf7saCNlspN/DNBsFdcTtwL516FfhpazYwBYZZzwrMnDwCv9DTKzOrSaGZk0fgQ5EH\n38xqVGDmZBL4bsM3s2EqM3PyCPyce+m4F43Z+Mk5cwYoj8CHfE+vct0vG22uSNSvwP/tTAK/zNMr\nM6tLmZmTSeBT5KetmdWowMzJI/AL7SJlZjUpNHPyCHwCdpV3emVtzHAbtw2YA78+UeDBtzZ2+e/B\nBimIAiuZeQR+ALv9D25mQ1Jo5uQR+FDk6ZWZ1ajAzMkm8KPAT1szq0+JmZNH4EcU+WlrZjUpNHO6\nDnxJ+wEbgecj4hOSjgHWALOBh4BPR8Sbkg4AbgZOArYDF0TEs53WX+IFlFEm3ylqo841/LauAp4A\nDk3PvwxcExFrJP0dcBlwbfq6IyKOk3RhWu6Ctmsu9ALKKIt956U3m75hdceNMiuZXQW+pPnAfwW+\nBPy5qurd6cDFaZGbgM9TBf7S9D3AbcA3JCk69Lt0t0wzG2Z33BIzp9sa/teAvwTem54fAbwaETvT\n80lgXvp+HrAFICJ2SnotLb+t9erLbE8zs7qUmTkdA1/SJ4CXImKTpNOmipssGl281rje5cBygAUH\nH+AmHTMbrgIzp5sa/seBT0o6FziQqg3/a8AsSTNTLX8+8EJafhJYAExKmgkcBryy90ojYgKYADhx\n9ns7tfiYmfVN1UmnvMzpGPgRsQJYAZBq+P89Iv5Q0veB86h66iwD7khvWZue359e/3FXaV7gBZSB\nci8as/YKzJzp9MP/HLBG0heBh4EbUvkNwLclbaaq2V/YcU0RRd4EMVg+ntZ/GpdB7QrNnJ4CPyLW\nA+vT908DJzdZ5tfA+T3vSYGnV2ajJsZpULsCMyePO22BKO/sysxqVGLm5BH4QT3taeNyempmvfEE\nKPWq5Yr5OJ2eWhY85MQI6WMbfj+GnpG0gmqkgl3AlRGxLpUvAVYD+wHXR8SqVN50G+32M4/Ajyiy\nT6yNHw85MSKi7xOgTGvoGUnHU3Vw+X3gKOCfJP1uWtc3gbOourw/KGltRDzeZhst5RH4UOTplZnV\nqE+Z04+hZ1L5moj4DfBM6uU41Slmc+okg6Q1wFJJT7TZRkvZBH6JXaTMrD59zJx+DD0zD3igYZ2N\n79myV/kpHbbRUjaB7xq+mQ1V95kzR9LGhucTaaSAfg4906p8Ro/Lt5VF4Eer9jRfADOzAYjoqYa/\nLSIWt3itX0PPTJVPaXxPs/JtbbbRUhaBD7S4aOtavw2Au+Nan0bL7NfQM5LWAt+V9FWqi7aLgJ9S\n1eQXpR45z1Nd2L04vee+FttoKY/AL3QgI6uJu+Pa4CdA6WnomYh4TNKtwOPATuDyiNgFIOkKYB1V\nt8wbI+KxDttoKY/AB7fhm9lw9Tlzpjv0TER8iaqnz97ldwF3NSlvuo128gl899Ixs2EqMHMyCfxw\nk46ZDVGZmZNH4EeZAxmZWU0KzZwsAj+A3QWeXtm+3BPXhqHUzMki8MG9dKziPwMblhIzp5tJzA8E\n/hk4IC1/W0SsfDejwbVU5gTyZlaXQjOnmxr+b4DTI+LfJe0P/ETS3cCf08NocJ02UuKnrZnVp8TM\n6WYS8wD+PT3dPz2CHkeDazeReaFzEZhZTUrNnK7a8NPg/puA46jGZv4lvY8Gt63dNnaXePTNrDYl\nZk5XgZ9u8f2wpFnAD4APNlssfe1qFDdJy4HlAPP2n8lu3+4+EuRxaGwcRBSZOT310omIVyWtB06l\n99Hg9l7XBDAB8B8OPrDE5rSRFAX+k9j4cZNOC5LeB7yVwv4g4EyqC7GtRmprOhpcp+2UePDNrD4l\nZk43Nfy5wE2pHX8GcGtE3CnpcXoYDa4TV/HNbJhKzJxueun8HPhIk/KeR4NrvQ3YNUZ3vc3w7aJm\neXM//HqN08HfNU4/jGXBdYj+CsarktmtjAK/vINv1i3/e/RfiZmTUeDXvQdmVpISMyePwPcUh2Y2\nTIVmThaBH5Q5NrWZ1aPUzMki8CGGOja1L4CZlc4zXtUnhjyuRXm/ZxsCd8cdIWPWFbxbWQR+qbc5\n23hxd9zRUWrmZBH4UOYFFDOrT4mZk1Hg170HZlaSEjMno8Av8OibWW1KzJwsAj8Cdjd0kfK1LzMb\nKI+lU69ouGJe4O/BBswTt1ijAE+AUqcCe0jZMBX4z23tlZg5eQR+lHkThJnVpNDMySLwAzfjmNnw\nlJo53UxxuAC4GfhtYDcwERGrJc0GbgEWAs8Cn4qIHZIErAbOBd4ALomIhzptp8QZ5HMjXy23gpSY\nON3U8HcCfxERD0l6L7BJ0j3AJcC9EbFK0tXA1cDngHOARelxCnBt+tpWiQc/NyWe4trg5VqNKLGS\n2c0Uh1uBren71yU9AcwDlgKnpcVuAtZTBf5S4OY0cfkDkmZJmpvW03wbVKcOZmbDUGrm9NSGL2kh\n1fy2G4Ajp0I8IrZKen9abB6wpeFtk6nsHYEvaTmwHODIGfsV2SfWzOpTYuZ0HfiSfgv4B+DPIuLf\n2rT3Nnthn0MbERPABMDv7f+eCDfqmNkQlZg5M7pZSNL+VGH/nYi4PRW/KGluen0u8FIqnwQWNLx9\nPvBCp22EH3744UcXj37px/YkLZB0n6QnJD0m6apUPlvSPZKeSl8PT+WS9HVJmyX9XNKJDetalpZ/\nStKyhvKTJD2S3vP11DGm5Tba6Rj4aeU3AE9ExFcbXloLTO3UMuCOhvLPpB/sVOC1du33UB3UXeFH\np0fd/2h++JHDo1/6tL2pTi0fBE4FLpd0PFUnlnsjYhFwb3oO7+zUspyqUwup1+NKqg4uJwMrGwL8\n2rTs1PuWpPJW22ipmyadjwOfBh6R9LNU9lfAKuBWSZcBzwHnp9fuouqSuZmqW+alXWyDEk+veuWb\nRa3flG0fmsEK+jN/Qb86taRl74mIVwBST8glktYDh0bE/an8ZuAPgLvbbKOlbnrp/ITWPavOaLJ8\nAJd3Wu8+7+v1DWY2bSVXtPr9k0+zU0u78skm5bTZRktZ3GkLFHnF3Mzq00PmzJG0seH5ROp08rY+\ndGrptfxdySLw+902Z2bWTo+Zsy0iFrd6sV2nllTz7qZTyyR7mmemyten8vlNlm+3jZayCHyA3Y58\nMxuifmROF51aVrFvp5YrJK2hukD7WgrsdcDfNFyoPRtYERGvSHo9dYDZAHwG+NsO22gpj8CP8Riq\n1EPRmI2OPkVOXzq1pGD/a+DBtNwXpi7gAn8CfAs4iOpi7d2pvNU2Wsoi8MelScfXIWwQXI/ov+hT\nJbOfnVoi4kbgxiblG4ETmpRvb7aNdrII/KqvgNPSrBn/ZwxCmZmTSeCPR5OOmY2OEjMnm8Av8Nib\nWY1KzJwsAn9c2vDNbDSUmjlZBD6UORmBmdWnxMzJJ/Dr3gEzK0qJmZNF4Jd6emVm9Sg1c7IIfMDz\nqZrZUJWYOfkEft07YGZFKTFzsgn8EtvTSuE7RS1HDvwmJN0IfAJ4KSJOSGWzgVuAhcCzwKciYkca\nSGg11VgRbwCXRMRDnbYRlHnF3Mzq0a8JUEZNN3Pafos9U2pN6Wn6rm7UPW2aH374MRqPfhn29nLQ\nzYxX/5xmcmnU0/Rdnea0BTfpmNlwlZg577YNv9fpu9pPYh5ljmthZvUoNXP6fdG26+m4JC2navbh\nCM3wBChmNlQlZs67Dfxep+/aR5oTcgJg4X77x67yjn023IvGShO4ht+Lnqbv6maFJR58s3E1YwRq\nESVWMrvplvk9qgu0cyRNAivpcfquToIo8vTKbFzlXoErNXO66aVzUYuXepq+q5Pc/0DMbLyUmDlZ\n3GkblNlFyszqUWrmZBH4UOanrZnVp8TMySLwI2Bngbc5d2sULoCZjZJSMyeLwIcyT6+6VWJNxAav\nm3FVxlmJmZNF4JfaJ9asTiUG3pRSMyeLwIcy73ozs/qUmDlZBH4AO6dx7N3GbWa9KvEMJ4vAh+kd\n/BJPzWzwSm/jHmdu0qlRdfALPPqWtRJrgKWoWhXKy5wsAh/8z2Vmw1Vi5uQT+OV92JpZjUrMnCwC\nv9TbnM2sHqVmThaBT0AU2J5mZjUpNHOyCPzpdsu0Mrk7rk2Ha/g1KvHg2/SU2AZr/VFqJTOLwA/C\n3TLNbGhKzZyBBL6kJcBqYD/g+ohY1ek9ruGb2TCVmDl9D3xJ+wHfBM6imtT8QUlrI+Lxdu/z6bmZ\nDVOJmTOIGv7JwOaIeBogTWi+FGgZ+KV2kTKzepSaOYMI/HnAlobnk8Ap7d5Q6m3O485j0ViuSs2c\nQQR+s85y+xxZScuB5enpb67j9UcHsC/DMgfYVvdOTMMo7/8o7zt4/9+N35nuCraxe911vD6n68XH\nxCACfxJY0PB8PvDC3gtFxAQwASBpY0QsHsC+DIX3vz6jvO/g/a9LRCypex/qMIiz7geBRZKOkfQe\n4EJg7QC2Y2ZmPeh7DT8idkq6AlhH1S3zxoh4rN/bMTOz3gykH35E3AXc1cNbJgaxH0Pk/a/PKO87\neP9tiFTiAEJmZiVyzzkzs0LUHviSlkh6UtJmSVfXvT97k7RA0n2SnpD0mKSrUvlsSfdIeip9PTyV\nS9LX08/zc0kn1vsTVCTtJ+lhSXem58dI2pD2/5Z0gR1JB6Tnm9PrC+vc77RPsyTdJukX6ffw0VE5\n/pI+m/5uHpX0PUkH5n7sJd0o6SVJjzaU9Xy8JS1Lyz8laVkdP4u9U62B3zAMwznA8cBFko6vc5+a\n2An8RUR8EDgVuDzt49XAvRGxCLg3PYfqZ1mUHsuBa4e/y01dBTzR8PzLwDVp/3cAl6Xyy4AdEXEc\ncE1arm6rgR9GxO8BH6L6ObI//pLmAVcCiyPiBKpODBeS/7H/FrB3t8Wejrek2cBKqpsuTwZWTn1I\nWI0iorYH8FFgXcPzFcCKOvepi32+g2qcoCeBualsLvBk+v464KKG5d9ersZ9nk/1T3o6cCfVzXHb\ngJl7/x6oeld9NH0/My2nGvf9UOCZvfdhFI4/e+46n52O5Z3AfxmFYw8sBB59t8cbuAi4rqH8Hcv5\nUc+j7iadZsMwzKtpXzpKp9gfATYAR0bEVoD09f1psRx/pq8Bf8me4UOOAF6NiJ3peeM+vr3/6fXX\n0vJ1ORZ4Gfj71CR1vaRDGIHjHxHPA18BngO2Uh3LTYzOsW/U6/HO5vdge9Qd+F0Nw5ADSb8F/APw\nZxHxb+0WbVJW288k6RPASxGxqbG4yaLRxWt1mAmcCFwbER8B/h97mhOayWb/UxPGUuAY4CjgEKom\nkL3leuy70WqfR/FnGXt1B35XwzDUTdL+VGH/nYi4PRW/KGluen0u8FIqz+1n+jjwSUnPAmuomnW+\nBsySNHUfRuM+vr3/6fXDgFeGucN7mQQmI2JDen4b1QfAKBz/M4FnIuLliHgLuB34GKNz7Bv1erxz\n+j1YUnfgZz8MgyQBNwBPRMRXG15aC0z1PFhG1bY/Vf6Z1HvhVOC1qVPhOkTEioiYHxELqY7vjyPi\nD4H7gPPSYnvv/9TPdV5avraaWUT8K7BF0gdS0RlUQ22PwvF/DjhV0sHp72hq30fi2O+l1+O9Djhb\n0uHpTOfsVGZ1qvsiAnAu8H+BXwL/o+79abJ//5HqVPTnwM/S41yqttV7gafS19lpeVH1PPol8AhV\nD43af460b6cBd6bvjwV+CmwGvg8ckMoPTM83p9ePzWC/PwxsTL+D/wMcPirHH/jfwC+AR4FvAwfk\nfuyB71Fdc3iLqqZ+2bs53sAfpZ9lM3Bp3X9HfoTvtDUzK0XdTTpmZjYkDnwzs0I48M3MCuHANzMr\nhAPfzKwQDnwzs0I48M3MCuHANzMrxP8HgYKs5/xWnywAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f06f6f64828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.pcolormesh(idx,cmap='Reds_r')\n",
    "plt.colorbar()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x7fea04f0ba20>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW8AAAD8CAYAAAC4uSVNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGYxJREFUeJzt3X+MXWWdx/H3Z1qhUoQWCqS2ZFtC42pMhG4DRTbGpYiA\nRvijbiBm6WKT2T9YxR+JgPtH465/SGKskmwaJxRF4/LDigthCWxTIMY/rLa0iwVkO6LCSKUU2upK\nqrbz3T/Oc6e3t3dmzp25d+555nxe5uSe85znnPPcU/zOc5/znOdRRGBmZnkZ6HcBzMyscw7eZmYZ\ncvA2M8uQg7eZWYYcvM3MMuTgbWaWIQdvM7MMOXibmWXIwdvMLENz+10AgHlSnJH+jqgpXWlLxxNa\n9nPSMSekq11elcrX7nzjp6WzTnCeMuc9/j1byqg2+Zpuyrh52pWp9dxt87RLTNdq+0XG+eLN12q9\n2Dj5xi9Uu2Mn3tZ4+8dWp7PdesMn2HfSf7gaf98J283ln+B64+WbsMwTnXui49pdo+3GOMe2u/5k\n6Sdm2rlr94GIOGecE5dy9YeujANvvFEq785du5+IiKunc71eqETwfgcDrGU+AxIDwIBInydua5L9\nA+kf+fj+Iq3dttrsVzqn0rrSNcXxnyhSERTE8XOMbTfO27Ktpms2ttV0TZq+29i6QANFvmLRCZ8D\nzdsD7fc3yjswcHy9cc6Blvxqc76iEDpxveWzccNPzMOJeVr3t0trl3es0I19Ax1ta7z9jfPOmXPi\nNSY830Aq60DLvqZtgIE5TWUfaDm2zXnKnLd5u3HelnNoLE+b84933bH0dHxrmdtev+n+jK3TPn3s\nk+NlbNzvTs7Tci7NX/Cb8WJJWQfeeIMdP366VF7NX7BoutfrhUoEbzOzGZf5uE4O3mZWQwEx2u9C\nTIuDt5nVk2veZmY5cvA2M8tL4Jq3mVmWHLzNzHITuNnEzCxHow7eZmYZcvA2M8uLH1iamWXKwdvM\nLDd+YGlmlqUY9evxZmYZcs3bzCw/mbd5l5pJR9JnJT0naY+k+yTNk7Rc0nZJeyU9IOmUlPfUtD2c\n9i/r5RcwM+tYRPmloiYN3pKWAJ8GVkXEe4E5wA3AncDGiFgBHATWp0PWAwcj4kJgY8pnZlYxUXKp\nprJzWM4F3i5pLnAasA+4AtiS9t8LXJ/Wr0vbpP1rNDaFhplZRcz2mndE/Bb4KvAyRdA+DOwEDkXE\n0ZRtBFiS1pcAr6Rjj6b8Z7eeV9KgpB2Sdhyp8F83M5ulYrTcMglJ90jaL2lPU9pZkramZuWtkham\ndEm6KzUrPytpZdMx61L+vZLWTXbdMs0mCylq08uBdwLzgWva3YrGIRPsO54QMRQRqyJi1by2h5iZ\n9UpX27y/DbROUHw7sC01K29L21DEzhVpGQQ2QRHsgQ3ApcAlwIZGwB9PmWaTK4FfRcTrEfEX4CHg\n/cCC1IwCsBR4Na2PAOenAs0FzgTeLHEdM7MZ1J0274j4ESfHuObm49Zm5e9E4ScUcXQx8GFga0S8\nGREHga2c/AfhBGWC98vAakmnpbbrNcDzwFPA2pRnHfBwWn8kbZP2PxlR4YYjM6un3j6vPC8i9gGk\nz3NT+lizctJoch4vfVyT9vOOiO2StgDPAEeBXcAQ8F/A/ZK+nNI2p0M2A9+VNEzx1+iGya5hZjaj\nOhuYapGkHU3bQxExNMUrj9esXKq5uVmpl3QiYgNFe0yzlyjaZlrzHgE+Xua8ZmZ9U372+AMRsarD\ns78maXFE7EvNIvtT+lizctJoch4BPtiS/vREFyjbVdDMbBbp+Us6zc3Hrc3KN6VeJ6uBw6lZ5Qng\nKkkL04PKq1LauPx6vJnVVHcexUm6j6LWvEjSCEUrxVeAByWtp3hu2GiNeAy4FhgG3gJuBoiINyX9\nG/CzlO9fI2LCjh4O3mZWT13qRxERN46za02bvAHcMs557gHuKXtdB28zq6fybd6V5OBtZvVT8Vff\ny3DwNrN6Onas3yWYFgdvM6uhcLOJmVl2PHu8mVmmXPM2M8tP7kMuOXibWQ0FePZ4M7PMBDDq3iZm\nZplxbxMzszy5zdvMLEOZ17zLzGH5Lkm7m5bfS/rMVCbYNDOrhAgYLblUVJnZ41+MiIsi4iLgbyiG\nMfwhHU6waWZWKaPHyi0V1elkDGuAX0bEb+h8gk0zs4pIDyzLLBXVafC+AbgvrXc6waaZWTU0Xo/v\n3Uw6PVc6eEs6BfgY8P3JsrZJO+kOSBqUtEPSjiNdmtHCzKy0GtW8rwGeiYjX0vZrjeaQkhNsniAi\nhiJiVUSsmtc23puZ9dBsf2DZ5EaON5lA5xNsmplVRP5t3qX6eUs6DfgQ8E9NyR1NsGlmVhkBUeGe\nJGWUCt4R8RZwdkvaG3Q4waaZWTX49XgzszxVuCdJGQ7eZlZPrnmbmWUmqt2TpAwHbzOrJ9e8zcxy\nE5Uet6QMB28zqx/PHm9mlinPYWlmlptqDzpVhoO3mdVT5g8sOx0S1swsfxFw7Fi5ZRKSPivpOUl7\nJN0naZ6k5ZK2p5nGHkijsiLp1LQ9nPYvm+pXcPA2s3rqwnjekpYAnwZWRcR7gTkU8x7cCWxMM40d\nBNanQ9YDByPiQmBjyjclDt5mVk/dG1VwLvB2SXOB04B9wBXAlrS/daaxxgxkW4A1kqY0JraDt5nV\nT9CV8bwj4rfAVylGVt0HHAZ2Aoci4mjK1jyb2NhMY2n/YVoG/SvLwdvMaqij8bwXNWb9Sstg4yyS\nFlLUppcD7wTmU0xc0+aCxSET7OuIe5uYWT2V7yp4ICJWjbPvSuBXEfE6gKSHgPdTTLw+N9Wum2cT\na8w0NpKaWc4E3pxK8V3zNrP66V5vk5eB1ZJOS23Xa4DngaeAtSlP60xjjRnI1gJPpjkQOlYqeEta\nIGmLpF9IekHSZZLOkrQ1dYXZmn4+kKY/uyt1hXlW0sqpFMzMrKe68MAyIrZTPHh8Bvg5RUwdAm4D\nPidpmKJNe3M6ZDNwdkr/HHD7VItfttnkG8DjEbE29Vc8DfgisC0iviLp9lSI2yjae1ak5VJgU/o0\nM6uOLr1hGREbgA0tyS8Bl7TJe4TjU0ZOy6Q1b0lnAB8g/eWIiD9HxCFO7PLS2hXmO1H4CUXbz+Ju\nFNbMrCu61Nukn8o0m1wAvA58S9IuSXdLmg+c15gVPn2em/KPdYVJmrvJmJlVQP6zx5cJ3nOBlcCm\niLgY+CMTt9OU6gojabDR9ebI1HrKmJlNXRfesOynMsF7BBhJDfNQNM6vBF5rNIekz/1N+c9vOr65\nm8yYiBiKiFURsWpe23hvZtYjXRzbpF8mDd4R8TvgFUnvSkmNrjDNXV5au8LclHqdrAYON5pXzMwq\nI/Nmk7K9TT4FfC/1NHkJuJki8D8oaT1FX8fGE9THgGuBYeCtlNfMrFoq3CRSRqngHRG7gXZvGK1p\nkzeAW6ZZLjOz3onwTDpmZlmqQ83bzGzW8ezxZmaZcbOJmVmm3GxiZpahCncDLMPB28xqqNrjlpTh\n4G1m9RO45m1mliU/sDQzy417m5iZ5SdwbxMzsyw5eJuZZcjB28wsM43xvDPm4G1m9eSat5lZbtzb\nxMwsT5nXvMvMYYmkX0v6uaTdknaktLMkbZW0N30uTOmSdJekYUnPSlrZyy9gZtaxRlfBWT4BccPf\nRcRFEdGYUed2YFtErAC2cXxG+WuAFWkZBDZ1q7BmZl1To+Dd6jrg3rR+L3B9U/p3ovATYEFjlnkz\ns2qowezxSQD/LWmnpMGUdl5jVvj0eW5KXwK80nTsSEo7gaRBSTsk7ThCdf+6mdksNAuaTco+sLw8\nIl6VdC6wVdIvJsirNmkn3YGIGAKGAM7RnOreITObnSocmMsoO3v8q+lzv6QfApcAr0laHBH7UrPI\n/pR9BDi/6fClwKtdLLOZ2bQEQWTeVXDSZhNJ8yW9o7EOXAXsAR4B1qVs64CH0/ojwE2p18lq4HCj\necXMrDJq0GxyHvBDSY38/xERj0v6GfCgpPXAy8DHU/7HgGuBYeAt4Oaul9rMbDq6OKqgpAXA3cB7\n05k/CbwIPAAsA34N/H1EHFQRSL9BESPfAv4xIp6ZynUnDd4R8RLwvjbpbwBr2qQHcMtUCmNmNjO6\nOrbJN4DHI2KtpFOA04AvUnSl/oqk2ym6Ut/GiV2pL6XoSn3pVC46na6CZmb56kKziaQzgA8Am4tT\nxp8j4hAz0JXawdvM6qd7XQUvAF4HviVpl6S707PBaXWlLsPB28zqqXzwXtR4JyUtg01nmQusBDZF\nxMXAHzn+tnk7pbpSl+GBqcyshjoaVfBA07AgrUaAkYjYnra3UATvnnelds3bzOonKIJ3mWWi00T8\nDnhF0rtS0hrgeWagK7Vr3mZWT93rw/0p4Hupp8lLFN2jB+hxV2oHbzOroe5NxhARu4F2zSo97Urt\n4G1m9VThtyfLcPA2s/rp4huW/eLgbWb15OBtZpabII7lPaqgg7eZ1U8Ao655m5nlx80mZmb5Cde8\nzcwyU/GJFsoo/Xq8pDlp1KxH0/ZySdsl7ZX0QHq7CEmnpu3htH9Zb4puZjYNo1FuqahOxja5FXih\naftOYGNErAAOAutT+nrgYERcCGxM+czMqiMgjo2WWqqqVPCWtBT4CMVUP6SpfK6gGEELTh5svDEI\n+RZgTcpvZlYZEVFqqaqyNe+vA18AGn+GzgYORcTRtN08oPjYYONp/+GU38ysIkqO5Z1z8Jb0UWB/\nROxsTm6TNUrsaz7vYGNw8yNTG4vczGzqMm/zLtPb5HLgY5KuBeYBZ1DUxBdImptq180DijcGGx+R\nNBc4E3iz9aQRMQQMAZyjOdW9Q2Y26xSV6rzDzqQ174i4IyKWRsQy4AbgyYj4BPAUsDZlax1svDEI\n+dqUP++7ZGazz7HRcktFTWcmnduAz0kapmjT3pzSNwNnp/TPMfF8bmZmMy+CGC23VFVHL+lExNPA\n02n9JeCSNnmOcHzWCDOzasq8QcBvWJpZLUV1W0RKcfA2s/rxZAxmZpmqcHt2GQ7eZlY/4ckYzMzy\n5GYTM7P8VLkbYBkO3mZWT655m5nlJcI1bzOzDFV7xMAyHLzNrH7SZAw5c/A2s3pyzdvMLENu8zYz\ny021pzgrYzpDwpqZ5SmKganKLGVImiNpl6RH0/ZySdsl7ZX0gKRTUvqpaXs47V821a/g4G1mtRPA\n6GiUWkq6FXihaftOYGNErAAOAutT+nrgYERcCGxM+abEwdvMaqlbs8dLWgp8BLg7bQu4AtiSstwL\nXJ/Wr0vbpP1rUv6OlZmAeJ6kn0r6H0nPSfpSSu/5zwIzs57obPL4RY3J0tMy2HK2rwNfABqNLGcD\nh9L8vlDM67skrS8BXgFI+w+n/B0rU/P+E3BFRLwPuAi4WtJqZuBngZlZr3RQ8z4QEaualqHGOSR9\nFNgfETubTt2uJh0l9nWkzATEERH/lzbflpZgBn4WmJn1QmMuhpI174lcDnxM0q+B+yni4teBBZIa\nvfmWAq+m9RHgfIC0/0zgzal8h1Jt3ulJ6m5gP7AV+CUz8LPAzKxXRiNKLROJiDsiYmlELANuAJ6M\niE8ATwFrU7Z1wMNp/ZG0Tdr/ZEyxz2Kpft4RcQy4SNIC4IfAu9tlS5+lfhakdqNBgNPbHmJm1iMR\njB7raT/v24D7JX0Z2AVsTumbge9KGqaocd8w1Qt0Onv8IUlPA6tJPwtS7brdz4KRiX4WpHajIYBz\nNCfv3vJmlpVeTGEZEU8DT6f1l4BL2uQ5Any8G9cr09vknFTjRtLbgSsp+jP2/GeBmVmvdKnNu2/K\n1LwXA/dKmkMR7B+MiEclPU+PfxaYmfVK7nXKSYN3RDwLXNwmvec/C8zMeqLiteoyPDCVmdVOAMc8\nqqCZWX5mfbOJmdlslHnsdvA2sxoK17zNzLITlB+ru6ocvM2shvKfScfB28zqJ9zbxMwsO714PX6m\nOXibWS252cTMLEOZx24HbzOrJ9e8zcxy47FNzMzyE9DryRh6zsHbzGop856CDt5mVkPhl3TMzLIT\ntJlYNzNlpkE7X9JTkl6Q9JykW1P6WZK2StqbPhemdEm6S9KwpGclrez1lzAz61SUXKpq0uANHAU+\nHxHvpph4+BZJ7wFuB7ZFxApgW9oGuAZYkZZBYFPXS21mNk2jEaWWqpo0eEfEvoh4Jq3/gWLy4SXA\ndcC9Kdu9wPVp/TrgO1H4CcUs84u7XnIzsykKYLTkUlVlat5jJC2jmM9yO3BeROyDIsAD56ZsS4BX\nmg4bSWmt5xqUtEPSjiOV/nFiZrNRHWaPB0DS6cAPgM9ExO8ljZu1TdpJtyAihoAhgHM0p8K3yMxm\no8i80liq5i3pbRSB+3sR8VBKfq3RHJI+96f0EeD8psOXAq92p7hmZt0x6x9YqqhibwZeiIivNe16\nBFiX1tcBDzel35R6nawGDjeaV8zMqiL34F2m2eRy4B+An0vandK+CHwFeFDSeuBl4ONp32PAtcAw\n8BZwc1dLbGY2TQEcq3KDdgmTBu+I+DHt27EB1rTJH8At0yyXmVlP5R26/YalmdVU5hVvB28zq5+q\nt2eX0VE/bzOz2WKUKLVMpJ/Dhzh4m1ktdam3Sd+GD3HwNrPaiSjG8y6zTHye/g0f4uBtZjUUpf8H\nLGoM5ZGWwXZn7ObwIWX4gaWZ1VIHM+kciIhVE2Xo9vAhZbjmbWa11K03LPs1fIiDt5nVTtnAPVnw\n7ufwIW42MbNa6tJEC30bPsTB28xqqRsTLfRz+BAHbzOrndnwhqWDt5nVUmQ+uImDt5nVUt6h28Hb\nzGoq9+BdZiadeyTtl7SnKa3ng66YmfVKYzKGMktVlenn/W3g6pa0ng+6YmbWS7lPgzZp8I6IHwFv\ntiT3fNAVM7NeGi25VNVU37Ds+aArZma90q1RBfup2w8sSw+6kkbmGgQ4fdw+7mZmvTHZRAtVN9Wa\n97QHXYmIoYhYFRGr5jl4m9kMCvKveU81ePd80BUzs146FuWWqpq02UTSfcAHKQYkHwE2MAODrpiZ\n9UqUmJ+y6iYN3hFx4zi7ejroiplZL1W5SaQMv2FpZrUTVLsbYBkO3mZWS655m5llJgKOVvjV9zIc\nvM2sltxsYmaWmUY/75w5eJtZLc36roJmZrORm03MzDLjZhMzswwF7m1iZpYlN5uYmWXIzSZmZpnx\n6/FmZjkKCLd5m5nlxzVvM7PMFL1N+l2K6XHwNrPaCYLRzJtNpjoN2oQkXS3pRUnDkm7vxTXMzKZj\ntORSVV0P3pLmAP8OXAO8B7hR0nu6fR0zs+mo6wTEE7kEGI6IlyLiz8D9wHU9uI6Z2ZQ0ugrmXPPu\nRZv3EuCVpu0R4NIeXMfMbEr8enx7apN20l2SNAgMps0/beIPe8Zy5XdPFwEH+l2Iaci5/DmXHVz+\nqfir6Z7gAKNPfJM/LCqdvYJ6EbxHgPObtpcCr7ZmioghYAhA0o6IWNWDsswIl79/ci47uPz9EhFX\n97sM09WLNu+fASskLZd0CnAD8EgPrmNmVltdr3lHxFFJ/ww8AcwB7omI57p9HTOzOuvJSzoR8Rjw\nWAeHDPWiHDPI5e+fnMsOLr9NkXIfnMXMrI568oalmZn1Vt+Dd9VfpZd0vqSnJL0g6TlJt6b0syRt\nlbQ3fS5M6ZJ0V/o+z0pa2d9vUJA0R9IuSY+m7eWStqfyP5AeLiPp1LQ9nPYv62e5U5kWSNoi6Rfp\n3+GyXO6/pM+m/272SLpP0ryq33tJ90jaL2lPU1rH91vSupR/r6R1/fgus1lfg3cmr9IfBT4fEe8G\nVgO3pDLeDmyLiBXAtrQNxXdZkZZBYNPMF7mtW4EXmrbvBDam8h8E1qf09cDBiLgQ2Jjy9ds3gMcj\n4q+B91F8j8rff0lLgE8DqyLivRQP8G+g+vf+20BrV7qO7reks4ANFC/oXQJsaAR865KI6NsCXAY8\n0bR9B3BHP8tUoswPAx8CXgQWp7TFwItp/ZvAjU35x/L1scxLKf4PdwXwKMWLVAeAua3/DhS9hC5L\n63NTPvWx7GcAv2otQw73n+NvG5+V7uWjwIdzuPfAMmDPVO83cCPwzab0E/J5mf7S72aTdq/SL+lT\nWSaVfsZeDGwHzouIfQDp89yUrYrf6evAFzg+VMPZwKGIOJq2m8s4Vv60/3DK3y8XAK8D30rNPndL\nmk8G9z8ifgt8FXgZ2EdxL3eSz71v1un9rsy/w2zV7+Bd6lX6KpB0OvAD4DMR8fuJsrZJ69t3kvRR\nYH9E7GxObpM1Suzrh7nASmBTRFwM/JHjP9nbqUz5UzPBdcBy4J3AfIpmhlZVvfdljFfmHL9LVvod\nvEu9St9vkt5GEbi/FxEPpeTXJC1O+xcD+1N61b7T5cDHJP2aYoTHKyhq4gskNfr5N5dxrPxp/5nA\nmzNZ4BYjwEhEbE/bWyiCeQ73/0rgVxHxekT8BXgIeD/53Ptmnd7vKv07zEr9Dt6Vf5VekoDNwAsR\n8bWmXY8AjSfo6yjawhvpN6Wn8KuBw42fm/0QEXdExNKIWEZxf5+MiE8ATwFrU7bW8je+19qUv281\npoj4HfCKpHelpDXA8+Rx/18GVks6Lf131Ch7Fve+Raf3+wngKkkL0y+Qq1KadUu/G92Ba4H/BX4J\n/Eu/y9OmfH9L8XPvWWB3Wq6laIvcBuxNn2el/KLoQfNL4OcUPQ36/j1S2T4IPJrWLwB+CgwD3wdO\nTenz0vZw2n9BBcp9EbAj/Rv8J7Awl/sPfAn4BbAH+C5watXvPXAfRRv9Xyhq0Ouncr+BT6bvMgzc\n3O//jmbb4jcszcwy1O9mEzMzmwIHbzOzDDl4m5llyMHbzCxDDt5mZhly8DYzy5CDt5lZhhy8zcwy\n9P/38nM+pjquIAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fea103974e0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.pcolormesh(rast,cmap='Reds_r')\n",
    "plt.colorbar()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# iterate through the cells\n",
    "for i in range(k):\n",
    "    for o in range(l):\n",
    "        "
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}