notebooks/Grid2D.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from grid.grid2d import Grid2D\n",
    "from mpl_toolkits.basemap import Basemap"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "xmin = 118.5\n",
    "xmax = 120.5\n",
    "ymin = 32.0\n",
    "ymax = 34.0\n",
    "xdim = 0.25\n",
    "ydim = 0.25\n",
    "ncols = len(np.arange(xmin,xmax+xdim,xdim))\n",
    "nrows = len(np.arange(ymin,ymax+ydim,ydim))\n",
    "data = np.arange(0,nrows*ncols)\n",
    "data.shape = (nrows,ncols)\n",
    "geodict = {'xmin':xmin,\n",
    "           'xmax':xmax,\n",
    "           'ymin':ymin,\n",
    "           'ymax':ymax,\n",
    "           'xdim':xdim,\n",
    "           'ydim':ydim,\n",
    "           'nrows':nrows,\n",
    "           'ncols':ncols}\n",
    "grid = Grid2D(data,geodict)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASIAAAEACAYAAAAX2nuLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFWJJREFUeJzt3XvMHXWdx/H3h164FagoQZBuigZcFuUmFBRYClaDhAUT\nNgqRgnjZROUi8QLFxYXdDSi73iKarEpJFZCwBVmMsFCVByEqtOUOxSpLkQKFSiPIpbRP+90/zlSe\nPjxnzmXmnN+cM59XMuk5c+b85ts+7be/329mfl9FBGZmKW2ROgAzMyciM0vOicjMknMiMrPknIjM\nLDknIjNLzonIzAqRNE/SQ5IekHSVpC0l7ShpkaTlkm6RND2vDSciM+uapJnAJ4EDIuKdwCTgROBc\nYFFE7An8InvflBORmRXxArAe2EbSZGAb4CngOGBBdswC4IN5jTgRmVnXImIN8DXgjzQS0J8jYhGw\nc0Q8kx32DLBzXjtORGbWNUlvAz4LzAR2BaZJOnnsMdF4jiz3WbLJJQTih9XMEooIdfvdbv79jjvf\ngcCvI+K5rL3rgHcDqyS9OSJWSdoFeDavzcKJqOFfOjx+BJjdwfHbd9h+N7br8PgbaAyDO9GP30en\n57gCOLnlUZvbscPju9Hpz+M7wGc6/E7Vfh4XA/M6bH+HDo9/vX/v4Nh/fv2uR4DzJW0NrAXmAHcB\nLwGnAl/Nfr0+r92SEpGZDaopBb4bEfdJ+iGwBNgI3A18j8b/JNdI+jiwAvhQXjtORGY1VzQJRMQl\nwCXjdq+h0TvqSwxdmpnmtKV6e+oASrJP6gBKclDqAEpwWJKzbp3krJtzIuqaE1G1zEodQAkOT3LW\nIkOzsnhoZlZzVUgCVYjBzBKqQo+o5Q2Nko6W9Iik30s6px9BmVn/TO5g62UMTUmaBFxKY/b7SWCx\npBsiYlkPYzKzPqpCj6hVkpsF/CEiVgBIuho4HnAiMhsSg5CI3gI8Meb9SuDg3oVjZv02CJfv23wO\nZWTM65kMx+V5syq6Hbij1BarcMWqVQxPAjPGvJ9Bo1c0zuzSAjKzPIez+f1GXync4iAMzZYAe2Sr\nsD0FfBg4qccxmVkfVb5HFBGjkk4HbqaxBORlvmJmNlwGoUdERNwE3NSHWMwsgcr3iMxs+A1Ej8jM\nhluRy/eS3g5cPWbXW4EvA7sBxwLrgEeB0yLi+WbteM1qs5qb0sE2XkT8LiL2j4j9gXcBLwPXAbcA\ne0fEvsByWiw96R6RWc2VmATmAI9GxBNsfiP0ncAJfYrBzAbRlE6ywGjupycCV02w/2PAj/O+6ERk\nVnOTc7LA7Rvgjo2t25A0FfgH4Jxx+78ErIuIiRLUazG0EaeZDbEpk5p/dtQkOGrM+6+81PTQDwBL\nI2L1ph2SPgocA7y3VQyJygl1quuyTe2b1vtT+BwdmN6HcwzDn9XVrQ9pJa9H1IGTGDP8knQ08AXg\niIhY2zKGUkIws4E1Zcti35e0LY2J6k+O2f1tYCqwSBLAbyLi083acCIyq7uCWSAiXgLeNG7fHn0M\nwcwGXgWyQAVCMLOkKpAFKhCCmSWVc9WsX5yIzOquAlmgAiGYWVIFr5qVwYnIrO4qkAUqEIKZJVWB\nLFCBEMwsKU9Wm1lyFcgCFQjBzJKqQBaoQAhmllQFskAFQjCzpHz53sySq0AW8OL5ZnU3qYNtApKm\nS1ooaZmkhyUdMuazz0naKGnHvBAqkAvNLKniWeBbwI0R8Y+SJgPbAkiaAbwPeLxVA+4RmdXd5A62\ncSTtABweEfOhUaZ+TP2yrwNfbDcEM6uzYjc07g6slnQ5sC+wFDiLRk9oZUTcn63QmMuJyKzucrLA\nyFMw8nTLbx8AnB4RiyV9E7gQOBx4/5jjcrORIqLNaJs0IAW0UW+k2Fl63D7DsZD6MJ3Di+e352oR\nEV3/A5EUcXoHx1/KZueT9GYa61Hvnr0/DLgAeAfwSnbYbsCTwKyIeHaidt0jMqu7AkOziFgl6QlJ\ne0bEchqL6C+NiDmbjpH0GPCuiFjTrB0nIrO6K54FzgCuzIosPgqcNu7zlsMuJyKzuitexeM+4KCc\nz9/a4xAyn+jxHI7nJHyOXhiGn3kJBRa9DIiZpVeBLFCBEMwsqQpkgQqEYGZJ+el7M0uuAlmgAiGY\nWVIVyAIVCMHMkvJVMzNLrgJZoOUyIJJmSLpV0kOSHpR0Zj8CM7M+KbAMSJkhtLIeODsi7pU0DVgq\naVFELOtdWGbWN4MwNIuIVcCq7PWLkpYBuwJORGbDYKvUAXTY2ZI0E9gfuLMXwZhZAhWYI2o7hGxY\nthA4KyJe7F1IZtZXgzA0A5A0BbgWuCIirn/dAUsveO31LrNh19llxGZm4y0ZgaUj5bZZsEckaQXw\nArABWB8Rs7L9ZwCfzvb/LCLOadpGqxUa1VhwdgHwXEScPcHnwSeKrfLYkp/E9jl6YRh+5geWsELj\njR0cfwyvO99EC59JOhI4DzgmItZL2ikiVjdrt50qHocCJwNHSron245uP3Qzq7SCdc0y45Php4CL\nI2I9QF4Sgvaumt2Byw6ZDa/iV80C+LmkDcB/RcT3gT2Av5d0EbAW+HxELGnWQAXmy80sqeKT1YdG\nxNOSdgIWSXqERm55Q0QcIukg4Bqg6UqNTkRmdZdXTmgpjNyd//WIeDr7dbWknwCzgJXAddn+xVnZ\n6TdGxHMdhmBmtZCTBWYf3Ng2ufCyzT+XtA0wKSL+ImlbGrXMLgReBI4CbpO0JzC1WRJqEYKZ1UKx\nLLAz8JOsmutk4MqIuCW75We+pAeAdcApvQvBzAZfsbpmjwH7TbB/PTC33XaciMzqrgJZoAIhmFlS\nXrPazJKrQBYoJ4TPl9JKc9PW9vgEsNW0l3t+ju126P2zwtvQ+99HP86xdV9+H6/0/Bzb8Zeetn9T\nGY0MTSIys8FVgSxQgRDMLKUYlGVAzGx4bahAFqhACGaWkhORmSX36pZTOzh6XU9icCIyq7kNk9JP\nEjkRmdXchgosWu1EZFZzo05EZpbahgqkAS8Ba1ZzG5jU9taMpEnZevY/zd7PknRXtm9xtkpjU+lT\noZklVdIc0VnAw8B22ftLgPMj4mZJH8jeH9nsy+4RmdXcq0xte5uIpN2AY4Af8Fo1j6eBHbLX04En\n82Jwj8is5kqYI/oG8AVg+zH7zgXukPSfNDo8785rwD0is5orMkck6Vjg2Yi4h81rm10GnBkRfwOc\nDczPi8E9IrOay5sjWjLyEktGcpdkeQ9wnKRjaFRI217Sj4BZETEnO2YhjWFbU05EZjWXdx/RfrO3\nZ7/Zr424vnfhnzb7PCLOo1FaGklH0CikOFfS3ZKOiIjbaFTzWJ4XgxORWc2VfB9RZL/+E/AdSVsC\nr2Tvm3IiMqu5sh7xyHo/t2WvlwAH53/jNU5EZjW3rsll+X5yIjKrOT9rZmbJVeFZs/QRmFlSXgbE\nzJJzIjKz5IZmjuiEt19RRjNNbd2HQnjDUphwWo8L+kF/ChMOy8+j1393yyiwuK4CNafdIzKrOQ/N\nzCy5oRmamdng8uV7M0vOQzMzS86JyMyScyIys+RercDl+7aWih1fKsTMhkePygntKGmRpOWSbpE0\nPS+Gdtes3lQqJFodaGaDpYxExOtzxLnAoojYE/hF9r6plomoSakQMxsSo0xqe5tIkxxxHLAge70A\n+GBeDO3MEU1UKsTMhkSPygntHBHPZK+fAXbOayA3grGlQiTNbnbcwxdc99fXO83ei51m75Uftpl1\nZdnIapaNrC61zbwh14qRx3l85PGmn7eTIyIiJOVO6yii+eeSLgLmAqNkpUKAayPilDHHxAnxo7xz\nFOaHXtvnh16rdY5e/92dq2uJiK6nTCTFeXF+28dfpH/b7HxNcsR1wEHA7IhYJWkX4NaI+Ntm7ebO\nEUXEeRExIyJ2B04Efjk2CZnZ4CtScrpJjpgL3ACcmh12KnB9XgydDg591cxsyPSonNBXgGskfRxY\nAXwo70ttRzC2VIiZDY8elRNaA8zJ/8ZrfGe1Wc35EQ8zS87rEZlZcl6PyMyS89DMzJJzyWkzS85z\nRGaW3NDMEc2lt4949OdW/H48UjAcj0Zs14fHSPry83hpbc/PsdVLvW1/bglteI7IzJJzIjKz5DxH\nZGbJDc0ckZkNLl++N7PkPDQzs+SKDM0kbUXjifstganA/0TEPEn/ARwLrAMeBU6LiOebtdNuFQ8z\nG1JFqnhExFrgyIjYD9gHOFLSYcAtwN4RsS+wHJiXF4N7RGY1V/TyfURsuulrKjAJWBMRD4855E7g\nhLw23CMyq7midc0kbSHpXhrVOm4dl4QAPgbcmBeDe0RmNVe05HREbAT2k7QDcLOk2RExAiDpS8C6\niLgqrw0nIrOayxuavTyymJdHlrTVTkQ8L+lnwIHAiKSP0ii8+N5W380tJ9QOSXF9vL9QG634WbNO\nzuFnzdo1DM+aaWcKlxN6WzzY9vGP6h3jywm9CRiNiD9L2hq4GbgQmAJ8DTgiIv7Uql33iMxqruB9\nRLsACyRtQWPO+UcR8QtJv6cxeb1IEsBvIuLTzRpxIjKruSL3EUXEA8ABE+zfo5N2nIjMas5P35tZ\nck5EZpbcq+v80KuZJbZhNH0aSB+BmSW1YdRDMzNLzInIzJIbXe9EZGaJbdyQPg2kj8DM0vLQzMyS\nW5s+DZQSwfGP3VJGM801XWCyRD1+ONHn6JB/5v0zmjoA94jMzInIzJKrQCLyUrFmdbe+g20cSTMk\n3SrpIUkPSjpz3Oefk7RR0o55IbhHZFZ3Gwp9ez1wdkTcK2kasFTSoohYJmkG8D7g8VaNuEdkVnej\nHWzjRMSqiLg3e/0isAzYNfv468AX2wnBPSKzuitpxVxJM4H9gTslHQ+sjIj7sxUaczkRmdVdCZPV\n2bBsIXAWsBE4j8aw7K+H5H2/ZSKSNB34AbA3EMDHIuK33QZsZhWTl4geGIEHR3K/LmkKcC1wRURc\nL+mdwEzgvqw3tBuNuaNZEfHshG20quIhaQFwW0TMlzQZ2HZsDWtJEf+X20Rxvrmtfufwz7wt+tfi\nVTy4toNKPidofBUPAQuA5yLi7CbneAx4V0SsadZsbo8oK5h2eEScChARo/Tnr4iZ9csEl+U7cChw\nMnC/pHuyfedFxE1jjmmZ6VoNzXYHVku6HNgXWAqcNabWtZkNugKX7yPiDlpcfY+It7Zqp9Xl+8k0\nSoV8NyIOoNHRPLfdIM1sABS4fF+WVj2ilTQuwS3O3i9kgkR0wTdfez37kMZmZuUbWdHYStX7grct\n5SaiiFgl6QlJe0bEcmAO8ND44y74bK/CM7OxZs9sbJtc+KsSGq3As2bt3Ed0BnClpKnAo8BpvQ3J\nzPpqEBJRRNwHHNSHWMwshUFIRGY25Ipdvi+FE5FZ3RV7+r4UTkRmdVf1q2ZmVgOeIzKz5DxHZGbJ\neY7IzJLz0MzMkhuaRPSdUlpp7sUetw/wQh/O0Y/fxxCssdO3c/RhQZsYhAKLFZgj8uL5ZnX3agfb\nOJLmS3pG0gPj9p8haVlWYuirrULw0Mys7ooNzS4Hvg38cNMOSUcCxwH7RMR6STu1asSJyKzuCgzN\nIuL2rHrHWJ8CLo6I9dkxq1u146GZWd1t6GBrzx7A30v6raQRSQe2+oJ7RGZ1lzc0+9MIPDfSaYuT\ngTdExCGSDgKuAXKXi3UiMqu7vEQ0fXZj22T5he20uBK4DiAiFkvaKOmNEfFcsy94aGZWd+s72Npz\nPXAUgKQ9gal5SQjcIzKzCS7Lt0vSj4EjgDdKegL4MjAfmJ9d0l8HnNKqHScis7orcPk+Ik5q8tHc\nTtpxIjKruwrcWe1EZFZ3fvrezJIbmodezWxwORGZWXKeIzKz5Apcvi+LE5FZ3XloZmbJeWhmZsn5\n8r2ZJeehmZkl50RkZsl5jsjMkqtAj8jrEZlZIZLmSXpI0gOSrpK0ZadtOBGZWdeyhfM/CRwQEe8E\nJgEndtpOKUOzFy4to5Wc9vtw52c/6iuu6cM5/tKHc/Tjz6of5xiWP6vEXqAxy7SNpA3ANsCTnTbi\nOSKz2ut+tjoi1kj6GvBH4BXg5oj4eaftOBGZ1V7ebPWvsm1ikt4GfBaYSaOI939L+khEXNlJBE5E\nZrWX1yN6d7ZtctH4Aw4Efr1pcXxJ1wHvAZyIzKwTrxT58iPA+ZK2BtYCc4C7Om3Eicis9grNEd0n\n6YfAEmAjcDfwvU7bcSIyq71idzRGxCXAJUXacCIyq730z3i0vKGxjLsmzazKRjvYeiM3EZV116SZ\nVVn5Nac71WpoVspdk2ZWZYWumpUit0cUEWuATXdNPgX8uZu7Js2sytIPzXJ7RO3eNXnxmPgO2wIO\n96O0Zj2xItvKlX6yutXQrK27Juf52ptZX8zMtk1uK6XV9AsSteq7PAIcImlrSaJx1+TDvQ/LzPqn\n4pPVZd01aWZVlr5H1HJQVcZdk2ZWZdWfIzKzoZf+8r0TkVntuUdkZsmlnyPyHT9mtVfsqpmkoyU9\nIun3ks7pJoIkiej2jSnOWq6OV36qqPtTB1CS36UOoAQrkp25+zurJU0CLgWOBv4OOEnSXp1GkCQR\n3TEEiWhx6gBK4kRUHSuSnblQj2gW8IeIWBER64GrgeM7jcBzRGa1V2iO6C3AE2PerwQO7rQRJyKz\n2it0+T7KiEARxdqRVEogZtadiFC33+3m3+/Y80k6BLggIo7O3s8DNkbEVzuKo2giMrP6kjSZxhTd\ne2ksFXQXcFJELOukHQ/NzKxrETEq6XTgZhoruF7WaRIC94jMrAL6evm+jBufUpM0Q9KtWUGBByWd\nmTqmbkmaJOkeST9NHUu3JE2XtFDSMkkPZ3MWA6fuRSr6lojKuvGpAtYDZ0fE3sAhwGcG9PcBcBaN\n9aUGuVv8LeDGiNgL2AfoeFiQmotU9LdHVMqNT6lFxKqIuDd7/SKNv/i7po2qc5J2A44BfgB0fdUl\nJUk7AIdHxHxozFdExPOJw+rG2CIVk6lhkYp+JqKJbnx6Sx/PX7rsf7L9gTvTRtKVbwBfoLHg3aDa\nHVgt6XJJd0v6vqRtUgfVKRep6G8iGuTu/+tImgYsBM7KekYDQ9KxwLMRcQ8D2hvKTAYOAL4bEQcA\nLwHnpg2pc+OKVOwKTJP0kaRB9Vk/E9GTwIwx72fQ6BUNHElTgGuBKyLi+tTxdOE9wHGSHgN+DByV\nLQk8aFYCKyNi06N/C2kkpkHz1yIVETEKbCpSURv9TERLgD0kzZQ0FfgwcEMfz1+KrIjAZcDDEfHN\n1PF0IyLOi4gZEbE7jUnRX0bEKanj6lRErAKekLRntmsO8FDCkLpV+yIVfbuhsawbnyrgUOBk4H5J\n92T75kXE/yaMqahBHjafAVyZ/ef2KHBa4ng65iIVvqHRzCrAKzSaWXJORGaWnBORmSXnRGRmyTkR\nmVlyTkRmlpwTkZkl50RkZsn9P9nWNkAvwPqQAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x109e92450>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(grid.getData(),interpolation='nearest');\n",
    "plt.colorbar();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The coordinates at the center of the grid are 33.000,119.500\n",
      "The data value at the center of the grid is 40\n"
     ]
    }
   ],
   "source": [
    "lat,lon = grid.getLatLon(4,4)\n",
    "print 'The coordinates at the center of the grid are %.3f,%.3f' % (lat,lon)\n",
    "value = grid.getValue(lat,lon)\n",
    "print 'The data value at the center of the grid is %i' % (value)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "sdict = {'ymin':32.7,'ymax':33.3,'xmin':119.2,'xmax':119.8,'xdim':0.33,'ydim':0.33}\n",
    "grid.interpolateToGrid(sdict)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAT0AAAEACAYAAADBf6HuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGl5JREFUeJzt3X+sHeV95/H3B/+iEIpl4bXBvl2nG6M1VRuTdB0Uw+ay\nSZBjUki6EY6rFEKjFKVyN8qqXUISCSztKmx/ZFtKSr3EqUhDIFGDqZPYwpe012uqxIkbG2hsFxxy\ns7YxF0vEJo0X8DXf/ePMvRwfz5w7PjP3nhnm85JGnjnnmfN8OcDXz5x55vkqIjAza4pz+h2Amdl0\nctIzs0Zx0jOzRnHSM7NGcdIzs0Zx0jOzRpnZ64mS5gFfBf4tMALcEBHHUtqNAC8Cp4CTEbGi1z7N\nzIoqMtL7JDAUEZcC306O0wQwGBGXO+GZWb8VSXrXAfcl+/cB7+vSVgX6MTMrTZGktyAiRpP9UWBB\nRrsAHpW0S9JHC/RnZlZY19/0JA0BC1Pe+nT7QUSEpKzn2VZGxBFJ84EhSfsjYkdv4ZqZFdM16UXE\nu7PekzQqaWFEPCfpYuD5jM84kvx5VNImYAVwRtLrkjTNbIpFRKGfoM72/9+i/RXR891bYDNwE/A/\nkz8f7mwg6TxgRkT8TNL5wDXA+uyPvL1AOFNhGBjscwydngY+1O8gOnyZasZ0T7+DSPFZ4LZ+B9Hh\nwlI+5b/nbPeZUnrrXZHf9O4E3i3pKeA/JcdIukTSt5I2C4EdkvYAO4FvRsS2IgGbWTXNyrn1W88j\nvYh4AXhXyuvPAtcm+88Ay3uOzsxqo8hl43SqS5x9sqTfAaT4tX4HkMIx5XdlvwOYMr/Q7wByctLr\nakm/A0hRxf+ZHVN+V/U7gClThUvXPJz0zKwUdUkmdYnTzCrOIz0za5S6JJO6xGlmFeeRnpk1ipOe\nmTWKp6yYWaPUJZnUJU4zq7i6XN66RoaZlWJmzi2NpBFJT0jaLel7yWvzJA1JekrSNklzM85dJWm/\npKcl3TpZnE56ZlaKggsOpJWVmLQkhaQZwN3AKuAyYK2kZd3idNIzs1IUGeklOtfYy1OSYgVwICJG\nIuIk8CBwfbdOnPTMrBQljPQ6y0rkKUmxCDjYdnwoeS2Tb2SYWSmypqx8N9kmcUZZifY3u5SkOOsV\n1530zKwUWaO4qzh9bZk/T2mTUVYiT0mKw8BA2/EArdFeJl/emlkpev1NT9J5ki5I9sfLSjzJayUp\nIKMkBbALWCppiaTZwJrkvK5xmpkVNitvNhk745UFwCZJ0MpJ90fENkm7gK9J+ggwAtwArZIUwL0R\ncW1EjElaBzwCzAA2RsS+bt0rohpFyFrX61UrDFRFb+t3ADWyst8B1MSFpVRDO35uzt5eqm81NDOz\nCbNm9DuCfJz0zKwUM2uSTWoSpplV3aw5/Y4gHyc9MytHTbJJTcI0s8qrSTapSZhmVnk1ySY1CdPM\nKs93b82sUWqSTWoSpplVXk3u3hZ+9jbPqqWS7kref1zS5UX7NLMKKmFBvelQKOnlWbVU0mrgTRGx\nFPhd4J4ifZpZRTUh6ZFv1dKJ1U8jYicwV1LaYoBmVmczcm59VjTp5Vm1NK3N4oL9mlnV1GSkVzSE\nvEu0dK6oUI2lXcysPBVIaHkUDTPPqqWdbRYnr6UYbttfkmxmVq4dwGPlf2zBbJLcI9gFHIqI35D0\nVeDS5O25wLGIOONGqKQR4EXgFHCyrZraVIT52qqlwLO0Vi1d29FmM7AOeFDSFUngo6QaLBiOmU2u\ncwH3O8v52OJTVj4O7AUuAIiINeNvSPoT4FjGeePlI1/I00mhpJe1aqmkW5L3N0TEFkmrJR0Afg7c\nXKRPM6uoAtlE0mJgNfA/gP/a8Z5orZp8dbePyNtX4avwiNgKbO14bUPH8bqi/ZhZxRW7M/u/gD8E\nfjHlvauA0Yj4Uca54+UjTwEbIuLebh3V5KdHM6u8jGwyPArDaXXMEpLeCzwfEbslDaY0WQt8pUvP\nZ5SPjIgdZxmmmdlZysgmg4ta27j1/3xGk7cD1yUPMpwL/KKkL0XEjZJmAu8H3pLVbUb5yMyk5xKQ\nZlaOHicnR8SnImIgIt4IfBD4+4i4MXn7XcC+iHg2rcsu5SMzOemZWTnKm5zcPo93DfBA+5uSLpH0\nreRwIbBD0h5gJ/DNiNg2WZhmZsXlLAHZTURsB7a3HZ8x2yMZ9V2b7D8DLD+bPpz0zKwcFXiuNg8n\nPTMrR02ySU3CNLPKq0k2qUmYZlZ5vrw1s0apSTapSZhmVnk1ySY1CdPMKq8mhYGc9MysHDXJJjUJ\n08wqrybZpCZhmlnl+e6tmTVKTbJJTcI0s8qrSTapSZhmVnm+vDWzRilhlZXp4PX0zKwcBdfTkzRD\n0m5J30iO75B0KHltt6RVGeetkrRf0tOSbs0TpplZccUvb08rAUlrMdHPRcTnsk5IauXeTWuF5cPA\n9yVtjoh9Wed4pGdm5Sgw0msrAfkFXivnKCYv7bgCOBARIxFxEngQuL7bCU56ZlaOYpe34yUgX217\nLYDfl/S4pI2S5qactwg42HZ8KHktk5OemZWjx8JA7SUgOX1kdw/wRlrLwR8B/jSl10h5rSv/pmdm\n5ci4ezu8G4b3dD0zswTkeANJXwC+kXLuYWCg7XiA1mgvkyLOOlFOCUkBt/c7jBp4W78DqJGV/Q6g\nJi4kIib77awrSRGP5Wx7JZn9SXoH8AcR8RuSLh6vaSvpE8B/iIjf6mg/E/gX4J3As8D3gLXdbmR4\npGdm5Sgnm4jXLln/SNKbk+MfA7dAqwQkcG9EXBsRY5LWAY/Qunje2C3hgUd6NeSRXn4e6eVT0khv\nV862v5490psOHumZWTlqkk0K372dbDa0pEFJx9tmVX+maJ9mVkE93r2dboVy81nMht4eEdcV6cvM\nKq4mI72iYU7MhgaQND4bujPp9e363cymSU1qZBS9vM0zGzqAtyezqrdIuqxgn2ZWRQUXHJguRUPI\nc+v3B8BARJyQ9B7gYeDStIbLbn9iYn/+4DLmDy4rGN7rz6V8p98h1MZlbOx3CJW0b/go+4aPThxv\nWl/SB1cgoeVRNMxJZ0NHxM/a9rdK+ktJ8yLihc4Pu+yO3ywYjplNZtngfJYNzp843rS+67S2/GqS\n9Ipe3u4ClkpaImk2sAbY3N5A0gJJSvZX0JobeEbCM7N6ixn5tn4rlJuzZkNLuiV5fwPwAeBjksaA\nE8AHC8ZsZhV0qiYjvcJhRsRWYGvHaxva9j8PfL5oP2ZWbY1JemZmAC/PmZ2z5StTGsdknPTMrBSn\nZlTgB7scnPTMrBSnqvCMWQ5OemZWijEnPTNrklM1SSf1iNLMKq/o5W2ygMku4FCycvIfA++ldefj\nR8DNEXE85bwR4EXgFHAyIlZ068eFgcysFKeYkWvrYrzu7fjjrduAX4mINwNPAbdlnBfAYERcPlnC\nAyc9MyvJy8zOtaVJq3sbEUMRMV4SciewuEv3uVdyctIzs1KcYmauLUNa3dt2vwNsyXgvgEcl7ZL0\n0cni9G96ZlaKrEvX7w+fYNfwiczz2uveShpMef/TwCsR8ZWMj1gZEUckzQeGJO2PiB1Z/TnpmVkp\nspLeWwYv4C2DF0wc/9X6M9Ybyax7K+nDtC5735nV73iZyIg4KmkTrcWNM5OeL2/NrBRjzMi1dYqI\nT0XEQES8kdaCJH+fJLxVtC55r4+Il9L6lHSepAuS/fOBa4Anu8XpkZ6ZlaKkeXrtdW//AphN65IV\n4DsR8XvtdW+BhcBDyfszgfsjYlu3Dpz0zKwUZTyGFhHDwHCyvzSjzbPAtcn+M8Dys+nDSc/MSvFK\nxnSUqnHSM7NS+NlbM2sUP3trZo3ipaXMrFGc9MysUfybnpk1yivM6XcIuTjpmVkpfHlrZo3iy1sz\naxRPWTGzRvHlrZk1ipOemTWKk56ZNcrLNZmyUmgRUUlflDQqKXPRPkl3SXpa0uOSLi/Sn5lVV9Fq\naJJmSNot6RvJ8TxJQ5KekrRN0tyM81ZJ2p/kmVsni7Poysl/DazKejNZ/vlNybpYvwvcU7A/M6uo\nKSgB+UlgKCIuBb6dHJ8mqZV7N608dBmwVtKybp0USnpJ8Y2fdmlyHXBf0nYnMFfSgiJ9mlk19bpc\nPKSXgKQtfyR/vi/l1BXAgYgYiYiTwIPA9d3inOoaGYuAg23Hh+heu9LMamoKSkAuiIjRZH8USBsw\npeWYRd3inI4bGZ1FeCO1FbD3jocm9ucPLmP+YNdRqpn1YN/wUfYNHy39c7MuXUeGf8JPhn+Sed5k\nJSABIiIkpeWOzHySZaqT3mFgoO14cfJaqsvu+M0pDsfMlg3OZ9ng/InjTev3lfK5WUlvYPCXGRj8\n5Ynj/7P+sc4maSUg/wYYlbQwIp6TdDHwfMrHd+aYAVqjvUxTfXm7GbgRQNIVwLG24aqZvY68zOxc\nW6eMEpC/TSt/3JQ0uwl4OKXbXcBSSUskzQbWJOdlKjTSk/QA8A7gIkkHgduBWck/yIaI2CJptaQD\nwM+Bm4v0Z2bVVeKzt+OXrHcCX5P0EWAEuAGgvQRkRIxJWgc8AswANkZE16FroSgjYm2ONuuK9GFm\n9VBSCcjtwPZk/wXgXSltJkpAJsdbga15+/ATGWZWCj+GZmaN4vX0zKxRvJ6emTWKL2/NrFFeSZmO\nUkVOemZWCv+mZ2aN4t/0zKxR/JuemTWKk56ZNYp/0zOzRvFvembWKJ6yYmaN4stbM2sUX96aWaP0\nevdW0rm0lpOaA8wG/i4ibpP0VeDSpNlcWosQn1FGVtII8CJwCjgZESu69eekZ2al6DXpRcRLkq6O\niBOSZgKPSboyItaMt5H0J8CxrI8ABpP19yblpGdmpSgyTy8iTiS7s2mtgDyRwCSJ1qrJV3f5iM4C\nZJmc9MysFC8zp+dzJZ0D/AD4d8A9EbG37e2rgNGI+FHG6QE8KukUsCEi7u3Wl5OemZUia6R3Yvj7\nnBje1fXciHgVWC7pQuARSYMRMZy8vRb4SpfTV0bEEUnzgSFJ+yNiR1ZjJz0zK0VW0pszeAVzBq+Y\nOH5h/V9lfkZEHJf0LeDXgeHkN773A2/pcs6R5M+jkjYBK4DMpDfVJSDNrCHGmJFr6yTpIklzk/1f\nAN4N7E7efhewLykGlHbueZIuSPbPB64BnuwWp0d6ZlaKAvP0LgbuS37XOwf4m4j4dvLeGuCB9sbt\nJSCBhcBDrXsdzATuj4ht3Tpz0jOzUhSYsvIkGZevEXFGrez2EpAR8Qyw/Gz6c9Izs1J4aSkza5SX\nX/GCA2bWIKfG6pFO6hGlmVXeqTFf3ppZg9Ql6RWepyfpi5JGJaXOjZE0KOm4pN3J9pmifZpZ9Yyd\nnJFr67cyRnp/DfwF8KUubbZHxHUl9GVmFfXqqXpcOBaOMiJ2SFoySbPcKyCYWU3V5PJ2OlJzAG+X\n9DhwGPiDjhUUzOz14KWGjPRy+AEwkCwQ+B7gYV5bDfU0S2f+9sT+lefAVX4y+Az/+HK/I6iPnf0O\noKJGkq10Y1PxoeWb8qQXET9r298q6S8lzUtb5fS2evxFYVZrS5Jt3PayPthJr0XSAuD5iAhJKwDl\nXdbZzGqkKUlP0gPAO4CLJB0EbgdmAUTEBuADwMckjQEngA8W7dPMKuhkvwPIp4y7t2snef/zwOeL\n9mNmFXeq3wHk41sFZlaOsZxbB0nnStopaY+kvZI+m7x+h6RDbQ82rErrVtIqSfslPS3p1snC9K0D\nMyvHS72dllUCktZ0t89FxOeyzpU0A7ib1grLh4HvS9ocEfuyzvFIz8zK0eNID1JLQP40OZ7swYYV\nwIGIGImIk8CDwPXdTnDSM7NyFEh6ks6RtAcYBf4hIn6YvPX7kh6XtHG8jkaHRcDBtuNDyWuZnPTM\nrBzFRnqvRsRyYDHwHyUNAvcAb6S1HPwR4E/TTj3bMP2bnpmVI2vKyg+HYe9wro9oLwHZVvcWSV8A\nvpFyymFgoO14gNZoL5OTnpmVI2vKyr8fbG3jvr7+tLclXQSMRcSxthKQ6yUtjIjnkmbvJ7204y5g\nabLoybO0qqd1nUbnpGdm5ej9iYzUEpCSviRpOa1L2B8Dt8DpJSAjYkzSOuARWjdANna7cwtOemZW\nlt6nrKSWgIyIGzPaT5SATI63Alvz9uekZ2blaMqzt2ZmgJOemTWMk56ZNUpTVlkxMwNqs8qKk56Z\nlaPHu7fTzUnPzMrh3/TMrFH8m56ZNYp/0zOzRvHlrZk1ipOemTWKf9Mzs0Z5ud8B5OOkZ2bl8OWt\nmTWKL2/NrFF6nLIi6VxgOzCHVjW0v4uI2yT9MfBe4BXgR8DNEXE85fwR4MUkgpMRsaJbfy4MZGbl\n6LEwUES8BFydFAb6NeDqpO7tNuBXIuLNwFPAbRk9BzAYEZdPlvDAIz0zK0uB3/RS6t6+EBF725rs\nBP5zl4+YrD7uBI/0zKwcJ3NuKVLq3u7taPI7wJaMngN4VNIuSR+dLMxCIz1JA8CXgH+TdPy/I+Ku\nlHZ3Ae8BTgAfjojdRfo1swrKmrJybBiOD3c9NSJeBZZLuhB4RNLgeAlISZ8GXomIr2ScvjIijkia\nDwxJ2h8RO7L6Knp5exL4RETskfQG4J8kDbVXI5K0GnhTRCyV9DZaBXyvKNivmVVN1uXtGwZb27j/\nuz6j4el1b4FhSR8GVgPv7HLOkeTPo5I2ASuAzKRX6PI2Ip6LiD3J/r8C+4BLOppdB9yXtNkJzJW0\noEi/ZlZBPV7eSrpI0txkf7zu7W5Jq4A/BK5PbnaknXuepAuS/fOBa0ivjzuhtBsZSbHdy2n94Nhu\nEXCw7fgQsJjWtbuZvV70vspKVt3bp2nd2BiSBPCdiPi99rq3wELgoeT9mcD9EbGtW2elJL3k0vZv\ngY8nI74zmnQcRxn9mlmF9Hj3tkvd26UZ7Sfq3kbEM8Dys+mvcNKTNAv4OvDliHg4pclhYKDteHHy\n2hk+2/alXXkOXOV7y2alG0m20jXhMTS1xpQbgb0R8WcZzTYD64AHJV0BHIuI1Evb2zxr0GzKLUm2\ncdvL+uCGPIa2EvgQ8ISk8WkonwJ+CSAiNkTEFkmrJR0Afg7cXLBPM6uiJqyyEhGPkeMOcESsK9KP\nmdVAEy5vzcwmNOTy1sysxYWBzKxRfHlrZo3ipGdmjeLf9MysUWoy0vMzD2bWKE56ZtYoTnpm1ij+\nTc/MSlKPOxke6ZlZSXorhybpXEk7Je2RtFfSZ5PX50kakvSUpG3jC42mnL9K0n5JT0u6dbIonfTM\nrCS9LZ3cpQTkJ4GhiLgU+HZyfBpJM4C7gVXAZcBaScu6RemkZ2Yl+X85tzOllID8KW2lJpI/35dy\n6grgQESMRMRJ4EHg+m5ROumZWUl6rwGZUgLyh8CCtrU3R4G02jpp5SgWdYvSNzLMrCRZs5O/m2zZ\nUkpAXt3xfkhKKzNx1qUnnPTMrCRZd2/fmmzj/jzzE9pKQL4VGJW0MCKek3Qx8HzKKZ3lKAZojfYy\n+fLWzErS893b1BKQtEpN3JQ0uwlIq8GzC1gqaYmk2cCa5LxMHumZWUl6nqeXVQJyN/A1SR+hVcvo\nBoD2EpARMSZpHfAIrRsgGyNiX7fOFFGNaoyS4vicfkdRff9YkzoEVdBZgNnSrQciorNM61lp/d7W\n/Xe711xRuL8iPNIzs5LUY5kVJz0zK0k9HkNz0jOzknikZ2aN4pGemTWKR3pm1ige6ZlZo6QvJlA1\nTnpmVhKP9MysUerxm16hZ28lDUj6B0k/lPTPkv5LSptBSccl7U62zxTp08yqqvelpaZT0ZHeSeAT\nEbFH0huAf5I0lPLs2/aIuK5gX9Nux6twVcWWZHiC1tKyVVLVmKpoBFjS5ximTgNGehHxXETsSfb/\nFdgHXJLStG/P2RXx2Kv9juBMVfyf2THlN9LvAKZUPUZ6pY1jJC0BLufM57wDeLukxyVtkXRZWX2a\nWZX0trTUdCvlRkZyafu3wMeTEV+7HwADEXFC0ntorYl1aRn9mlmV1GPKSuGlpSTNAr4JbI2IP8vR\n/sfAWyPihY7Xq7HGlVkDlbO01PT1V0ShkZ4kARuBvVkJT9IC4PlkjfsVtBLtC53t+vklmFkxdfr/\nt+jl7UrgQ8ATySqnAJ8CfgkgIjYAHwA+JmkMOAF8sGCfZmY9q8zKyWZm06Evs9AkzZM0JOkpSdvG\ni4KktBuR9EQyqfl7UxTLKkn7JT0t6daMNncl7z8u6fKpiONs45ruSd+SvihpVNKTXdpM6/c0WUz9\nmBifZ8J+0m66vys/SDAuIqZ9A/4I+G/J/q3AnRntfgzMm8I4ZgAHaM0XnQXsAZZ1tFkNbEn23wZ8\ndxq+nzxxDQKbp/Hf2VW0piQ9mfF+P76nyWKa1u8o6XMhsDzZfwPwLxX5bypPXNP+ffVj69fzBtcB\n9yX79wHv69J2Kn8gXQEciIiRiDgJPAhc39FmItaI2AnMTW7OTKU8ccE0TvqOiB3AT7s0mfbvKUdM\nMM0T4yPfhP1+fFev6wcJzka/kt6CiBhN9keBrH/hATwqaZekj05BHIuAg23Hh5LXJmuzeApiOdu4\nqjbpux/f02T6+h11mbDf1++q6Q8STNkqK5KGaA2pO326/SAiosscn5URcUTSfGBI0v7kb/ey5L2L\n0/m331Tf/cnz+VWc9D3d39Nk+vYdTTJhH/r0XflBgikc6UXEuyPiV1O2zcCopIUAki4Gns/4jCPJ\nn0eBTbQu+8p0GBhoOx6g9bdutzaLk9em0qRxRcTPIuJEsr8VmCVp3hTH1U0/vqeu+vUdJRP2vw58\nOSIeTmnSl+9qsrgq+N/UlOjX5e1m4KZk/yZaf6OcRtJ5ki5I9s8HrgEy7xz2aBewVNISSbOBNUls\nnbHemMRxBXCs7dJ8qkwal6QFyeRwuk36nkb9+J666sd3lGfCPn34rvI+SFCx/6amRL8WEb0T+Jqk\nj9BaeOIGAEmXAPdGxLW0Lo0fSv4dzATuj4htZQYREWOS1gGP0LpjujEi9km6JXl/Q0RskbRa0gHg\n58DNZcbQa1xM86RvSQ8A7wAuknQQuJ3WneW+fU+TxUR/JsZPOmG/H99VnrhoyIMEnpxsZo1SsSUy\nzcymlpOemTWKk56ZNYqTnpk1ipOemTWKk56ZNYqTnpk1ipOemTXK/wfpQ4ijfqaVCQAAAABJRU5E\nrkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10a4a3890>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(grid.getData(),interpolation='nearest');\n",
    "plt.colorbar();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}