{ "metadata": { "name": "", "signature": "sha256:3cbddd3612af71b216e58d4d95995eb7209344bebb859b95002415549343601d" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "%run ../common.ipynb" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "img = imread('../MP.tiff')\n", "# make a circle mask\n", "r = 30\n", "size = 2*r\n", "x,y = np.meshgrid(range(size), range(size))\n", "R = np.sqrt((x-r)**2+(y-r)**2)\n", "mask = np.ones((size,size), bool)\n", "mask[R>r] = 0\n", "# discard values below 10\n", "img[img<10]=0\n", "\n", "# opening\n", "imshow(ndimage.binary_opening(img, structure=mask))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 6, "text": [ "(,\n", " ,\n", " )" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAfQAAAGcCAYAAAAvcVK9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X9sHPWd//HXOk5iJ9i7MTS546QvKJdLr+KbgFM3Ze14\nahvHjarkik5XCYnk29JTvl96HK6d+zYcRFVbkTRBoQflW3H1H1DRXnK69nSgpr+ArENsirvBmAbZ\n4kpyKndVyBEoWW+yJDj2zveP1EvseNfr9czOZ2aeD7QS7ITZ2cnMvOb9ns/ORGzbtgUAAHytwusF\nAAAA80egAwAQAAQ6AAABQKADABAABDoAAAFAoAMAEACuBnp3d7csy1JXV5ebHwMAgG+dPn1a69at\nU3V1tbLZ7JRpb731ltra2tTU1KREIlFwPq4F+tDQkDKZjPr6+jQ2NqbBwUG3PgoAAN+qq6tTb2+v\nbr311qum7du3T3v27NFzzz2n3bt3F5yPa4GeTCbV0dEhSWpvb9fAwIBbHwUAgG8tXrxYsVhsxmnD\nw8OKx+NaunSpampqdO7cubzzcS3QU6mUampqJEnRaFSpVMqtjwIAIJAmJiZy/z5blla6tRDRaFTp\ndFqSNDo6etXZRyQSceujAQAB5OWdyp3OrGK/S0XFh3V3Op3WsmXL8v/ZeS9VHvF4PHcBP5FIKB6P\nu/VRAAC4zrZtR16zfcaV1q5dq1/96lfKZDJKp9O65ppr8v6/rgV6fX29qqqqZFmWKisr1dDQ4NZH\nAQDguqxtO/Kabnx8XO3t7Tp+/Lg2bdqkY8eOqbOzU5K0c+dO7dq1Sxs3btSuXbsKLl/Eq6et0XIH\nAMyF1y33S9lxR+a1sKLSle/CjWUAAAgA1wbFAQAQJF52CIpBoAMAUISszA50Wu4AAAQAFTowBzO1\n3BjgCYRD1s7O/oc8RKADRSh07SzftFKC3ulrdJxsAM4xveVOoAPTOBWqk/MpJlTdGmzj5MkGALMR\n6MAfuBmqhQLUi5GzpXwmJwEIu5luCmMSAh2Q+6GaL9RN/xnMlebScQCCiJY74LArQ3B6uMxWDc82\nP8yOgYGAmQh0GG+uA9Kmv2dK2JRysuEXVO8IA9NP/gl0GMvJwWn5gqbcO2jQg8/UkynACWb/aI1A\nh2G8GphWbqaf6Tul0OURAM4i0OG5coXblaEelkA1CdU7/I5R7kAefvm5FtxB9Q6/IdCBPyBMkU/Q\nxxYA5UCgw3UEOYo1/bIIAQ+T8Dt0hBZBjlJcud1w3R0myRp+SCPQ4TiCHG6hNQ/kR6Bj3ghwlBsD\n6uAFBsUhcAhwmISqHeVi+pGvwusFgL8Q5jAV2ybCjgods+JACb9gZDzcxKA4+BZBDj/i+jrcYvo1\ndFrumBFhjiCwbZttGaFBhY6rcABE0DBwDk6g5Q7fIMgRdAQ75sP0IySBDoIcoUOwI4gI9JAixAFG\nxWNuaLnDKAQ5MBWj4lEs0wOdUe4hQpgDhTEqHl7p7u6WZVnq6uqa8v6vf/1rbdiwQZZl6cUXXyw4\nDwI9BDhIAcD82bYzr+mGhoaUyWTU19ensbExDQ4O5qZ97Wtf0w9/+EM9++yz2rNnT8HlI9ADjiAH\n5o6TYMwkazvzmi6ZTKqjo0OS1N7eroGBgdy0s2fP6vrrr1d1dbUymYwuXryYd/kI9ADjgATMD/sQ\nruRWhZ5KpVRTUyNJikajSqVSuWkf+chHNDIyonfeeUfDw8MaHR3Nu3wMigsoDkSAMxgJj/l6+cUX\n9XKB69/RaFTpdFqSNDo6qlgslpu2b98+/e3f/q1qamp0880367rrrss7n4jt0ZGfHcQdBDngDo5Z\n3vPy+BaJRHT892cdmdfN1y6b8l1effVV9fT06Lvf/a7uuece3XXXXWpoaJjy/7z77rvasWOHvv/9\n7+edLy33AJi83keYA+5hP4NbLff6+npVVVXJsixVVlaqoaFBnZ2dkqQnnnhCbW1t+vznP68HH3yw\n4PJRofsYBxbAWxzHysvrCv3X7zpTod9y3TJXvgvX0H2GEAfMwfX1cDH98Eug+whhDpiHUA8P7hSH\neeO6HWA29k+YgArdcBwoAMAMph+OCXSDEeYAYA7bNvvSCi13QxHmAIC5oEI3EGEO+A+PYQ0+0w/N\nBDoAOGwy3An2YLGzXi9BYbTcDUN1DgQH+zPKiQodAFxEtR4cpp+fEegG4WweCC5uQON/jHJHUQhz\nIPi4SRTcRIUOAGVGte5Ppg+KI9ANwBk7APiA4YdqWu4eI8yBcGLfh9Oo0D3EDg2EGyPg/cX0QXEE\nukcIcwCTuKbuD6ZfQ6flXmaMcgUwE44LmC8q9DJihwVQCJW64Qw/hBPoZUKYAygGoW4uWu4gzAHM\nCZfmUAoqdJexUwIoFaPgDcMo9/AizAE4gTa8GWi5hxRhDsBJHFMwGyp0F7DjAUAAGX5oJ9AdRpgD\ncAutd4/Rcg8PwhyA2zjOIB8qdAAAimF4hU6gO4SzZgDlQuvdG6Yf5mm5O4AwBwB4jUCfJ8IcgBc4\n9ngg69BrBt3d3bIsS11dXVPeHxgYUGNjo+LxuHp6egouHoE+D+xQALzELWLLzKVAHxoaUiaTUV9f\nn8bGxjQ4OJibtn//fv3whz/USy+9pO9973sFF49ABwCfI9T9LZlMqqOjQ5LU3t6ugYGB3LRrr71W\nqVRKFy9e1NKlSwvOh0AHgAAg1MvAdug1TSqVUk1NjSQpGo0qlUrlpt17773atGmTPvaxj2nbtm0F\nF49R7gAQEIx+d1mJP1t79ZV+/fqV/rzTo9Go0um0JGl0dFSxWCw3befOnUomk1q+fLk2btyoO+64\nQ1VVVTPOhwodAAAX1X+8WXf97wdyr+ni8bgSiYQkKZFIKB6P56a9//77ikajWrhwoSoqKnTp0qW8\nn0OgA0CA0Hp3TyTrzGu6+vp6VVVVybIsVVZWqqGhQZ2dnZKk++67T7fddpsaGxvV1taWa83PuHy2\nR3/7QWkLsfMAME1Qjq/TeXm8jUQi6nvpnCPzshprXPkuXEMvEsENwC8mj1dBDXbMjEAvAmEOwI8Y\nJOcww6OAQC+AIAcA5GTNzgQGxeVBmAMIAo5l4UGFDgBAEWYaoW4SAn0GnNECAK5ieKDTcgeAgKNI\nCYc5B/qbb76pFStWqLW1VZs2bZJ0+Wkwzc3N2rp1q8bHxyVJBw4cUFNTk7Zs2aJz55z57V45sOED\nCCKObfMXydqOvNxSUoXe0dGhI0eO6Be/+IXOnDmjF154Qf39/Vq7dq2eeeYZXbp0ST09Perv79e2\nbdtmfYarKdjgAQQZj1udJ5cezuKUkgL9yJEjsixLjz76qF555RW1tLRI+vCxbydPntSaNWtUUVFx\n1aPgAADeItSDac6D4q6//nqdOHFCixYt0mc/+1mdO3dOy5cvlyTV1tYqlUoplUqptrZ2ynsAAPiZ\nm+1yJ8w50BctWpT7982bN6u2tlanTp2SJKXTacVisSmPgpt8z3ScsQIIE+4iVwLDA33OLffz58/n\n/v2Xv/ylVq1apaNHj0qSDh8+rHg8rtWrV2t4eFjZbDb3HgAAcM+cK/T+/n599atf1eLFi2VZltav\nXy/LstTc3KwbbrhBO3bsUGVlpbZv367m5mbV1dXp4MGDbiy7Y6jOAQCzMb3lzuNTRaADCCeTjsPF\n8PrxqS/95D1H5tW4uc6V78KNZQAgpChmgiX0t35lgwYQZgyOK57pLfdQBzphDgCEerEihmdGaFvu\nhDkAIEhCXaEDAFC0rNmPWyPQAQC03YvANXQAgC8Q6oWZHuihvYYOAECQhLJCZ0AcAGDODK/QQxfo\nhDkA5Dd5jKT1frWIbfaguFC13AlzACgOx0v/CU2FzsYJAJgP0wfFhSbQAQBzw6j3aQz/HXooWu5U\n5wBQGo6f/kGFDgAoiEr9sojhFTqBDgCYFaHOw1kAAMAsuru7ZVmWurq6rnq/tbVVra2tqqurKziP\nwAc6138AwBlhP55GsllHXtMNDQ0pk8mor69PY2NjGhwczE175JFHdOTIET3yyCPavHlzweULfKAD\nAOCIbNaZ1zTJZFIdHR2SpPb2dg0MDFz1Z/7t3/5Nf/mXf1lw8QJ9DT3sZ5MAAO8N/ntSg79J5p2e\nSqW0cuVKSVI0GtXIyMhVf+bZZ5/VAw88UPBzAh3oAABnhXlwXKm3fv3ERz+hT3z0E7n/7vnx/5sy\nPRqNKp1OS5JGR0cVi8WmTD9x4oT+5E/+RFVVVQU/J7Atd6pzAHBHaI+vLrXc4/G4EomEJCmRSCge\nj0+Z/vTTT8/abpcCHOgAAPhBfX29qqqqZFmWKisr1dDQoM7Oztz0n/70p9qyZcus84nYHp1qudmy\nCe3ZIwCUiRdtdy+P7ZFIRMe/O+zIvG6++3+68l0Cdw2dMAcAuIHHpwIAAofiyTyBqtDZwAAAruFe\n7uVBmAMAXEXLHQAAuC0wFToAoLzCdpOZSHbC60UoiEAHAKAYtNwBAIDbqNABAChCxKblDgAIqFBd\nRzf8GnogWu78ZA0AEHZU6AAAFMPwQXG+D3SqcwDwVmja7oZfQw9Eyx0A4C2KK+/5ukJnAwIAlI3h\ng+J8HegAAJSL6T9b823LneocAMzCcdlbVOgAABTD8AqdQAcAOCbII95twwPdty13AADwIV9W6Fyn\nAQCUneEVui8DHQBgrsC23Q0PdFruAADH0UktPyp0AACKYXiF7rtA56wPAOAFRrkDAADX+a5CBwDA\nE4ZX6L4KdNrtAADPGB7otNwBAAgA31ToVOcA4C9B+z26bY97vQgFUaEDAFAEWxOOvGbS3d0ty7LU\n1dU15f2LFy/qr//6r3Xbbbfpy1/+csHlI9ABAPDQ0NCQMpmM+vr6NDY2psHBwdy0xx57THfeeacS\niYS+/e1vF5wPgQ4AQDHsCWde0ySTSXV0dEiS2tvbNTAwkJt29OhR/fjHP1Zra6sOHTpUcPF8cw0d\nAAAvlXpjmZffPKHBN0/knZ5KpbRy5UpJUjQa1cjISG7af/zHf6i7u1vf/OY31dLSos985jNasGDB\njPPxRaAzIA4A/CloA+NK8Ykb/0yfuPHPcv/d0/fzKdOj0ajS6bQkaXR0VLFYbMq0T33qU1q4cKFW\nrVqlt99+W9dff/2Mn0PLHQCAItgad+Q1XTweVyKRkCQlEgnF4/HctMbGRh0/flwTExN68803tXz5\n8rzLR6ADAFAE255w5DVdfX29qqqqZFmWKisr1dDQoM7OTknSfffdp127dmnDhg3avn27KivzN9Yj\ntkf97Lm0YGi5A4C/OdF29zILIpGIhr76sCPzWvfg/3Xlu/jiGjoAAF7L9xtyUxDoAAAUgcenAgAA\n11GhAwBQhCwtdwBA2AXh9+g8nAUAALiOCh0AgCIwyh0AgADIGj7KnUAHAKAIplfoXEMHACAAfFGh\nRyIRbv8KAPAUP1tzyOTPHQh2AIAXTA90Wu4AANf5/TfofuCbCh0A4F9BuLEMo9wBAAgARrkDAADX\nUaEDAFAE0wfF+SrQGeEOAP7k9+vnkjRheKD7puVOmAMAkJ8vKnTCHAD8KwjVuWR+y71ghX769Gmt\nW7dO1dXVymazkqT9+/erublZW7du1fj45WfDHjhwQE1NTdqyZYvOnTsnSert7VVjY6Pa2tp06tQp\nl78GAMBEQQlz6XKgO/FyS8FAr6urU29vr2699VZJ0pkzZ/TCCy+ov79fa9eu1TPPPKNLly6pp6dH\n/f392rZtm3p6eiRJu3fv1vPPP699+/Zp7969rn0BAIBZIpFI7oXyKRjoixcvViwWk3S57T04OKiW\nlhZJUnt7uwYGBnTy5EmtWbNGFRUVufcuXLig6upqLV26VOvXr9fIyEjJC0i7HQD8JajHbdMr9Dld\nQx8dHVVtba0kqba2VqlUSqlUquB7kjQxYfZ1BwAAZpNV1utFKKjoUe6RSETRaFTpdFqSlE6nFYvF\nZn1PkhYsWFDSwgX1LA8Ags627atecFfRgW7bthoaGnT06FFJ0uHDhxWPx7V69WoNDw8rm83m3luy\nZIkuXLigTCajY8eO6aabbnLtCwAA/MHvwT7h0D9uKdhyHx8f16ZNm3T8+HFt2rRJe/bskWVZam5u\n1g033KAdO3aosrJS27dvV3Nzs+rq6nTw4EFJ0q5du7Rx40ZVV1frqaeecu0LAABQDqa33CO2R6dL\nxY5+9PPZHABgZqWMgPcyDyKRiP71/r9yZF5/tfdfXfkuRt9YhjAHAJjC9Fu/Gh3oAACYwvSWu2/u\n5Q4AQFB1d3fLsix1dXVNef/rX/+6brnlFrW2tuqRRx4pOA8CHQBQdn68pOrWKPehoSFlMhn19fVp\nbGxMg4ODuWmRSETf+ta3dOTIEXV3dxdcPgIdAIAiOHOfuKvb9slkUh0dHZI+vAvrle677z5t3LhR\nx48fL7h8XEMHAMBFr//n7/Xv//X7vNNTqZRWrlwpSYpGo1Nul97Z2amvfe1rOnnypL74xS+qr68v\n73wIdAAAilDqoLiP3rBMH71hWe6/n3nx5JTpV95ddXR0NPcMFUlatuzy/7dq1apZP4eWOwAARXCr\n5R6Px5VIJCRJiURC8Xg8N23ykeTvvvtu7pHl+RDoAAB4qL6+XlVVVbIsS5WVlWpoaFBnZ6ck6Stf\n+Yo2bNigv/iLv9BDDz1UcD5G3ynOj6MgAQCz8+Od4nrub3dkXv9n7+Hw3SkOABBMtm2XFOpemohw\nY5mSUJ0DAFA8KnQAAIow04A2kxgZ6FTnAADTmB7oxrbcAQBA8Yys0AEAwea3AXGSNCGzu8dGBnok\nEqHtDgAwCi13AADgOiMrdIkqHQBgFlruAAAEAC13AADgOip0AACKQMt9HriODgAwhemBTssdAIAA\nMLpCBwDAFKZX6AQ6AABFmIiYHei03AEAZcf4KOcZX6EzMA4AYAJa7gAABIDpgU7LHQBQdn582prp\nfFGh03YHAHjN9ArdF4EuEeoAEBR+rc4nvF6AWfiq5e7XjQAAALf5pkKfVEylfmXwU9UDAJxAy90F\nk4F9ZVjnq94JdwCAEyYMjxBfBvokWvAAAFzm60AHAKBcaLkDABAAjHIHAACu822FXmiAW6EBcgyM\nAwDv+HnsExW6C2YLZUIbAOC0CYdeM+nu7pZlWerq6rpqmm3buuWWW/TEE08UXD7fBXqxYW3b9ox/\n1s9nhwCA4BkaGlImk1FfX5/GxsY0ODg4ZfqhQ4e0fPnyWfPLd4EOAIAX3KrQk8mkOjo6JEnt7e0a\nGBiYMv2f//mfdccdd8xa0Pr2GnqxbNumKgcAA/j9WFzqjWXO/Od5vfNf5/NOT6VSWrlypSQpGo1q\nZGQkN+25555TS0uLFixYoPHx8YKfQ4UOAEARSq3Ir73hGv158x/lXtNFo1Gl02lJ0ujoqGKxWG7a\nE088oS984QtFXW4m0AEA8FA8HlcikZAkJRIJxePx3LQ33nhDt99+u/7hH/5Bjz76qN5444288wl8\nyx0AACe49bO1+vp6VVVVybIs1dfXq6GhQZ2dnXrsscf06quvSpKeeuopTUxMaPXq1XnnE7E9+o1X\nqddSSlncmT6Ln7YBQPk4cf3cy+N2JBLR5vtvdmReP9l73JXvEtqWu98HZwAAcCXftdznerc3ghsA\nvBWU4/CEbfb38F2gz0VQNiIAgPdMv/WrLwN9tip9tiDn+jkAIGh8GejSh6F9ZTgXU5ET5gBQPkHq\nlGZpubsrSBsLAMBctNwNQWUOAAiywAc6QQ4AcAItd48Q5ADgraBdEjW95R7IG8sQ5gDgraCFuR8E\nrkInzAEAbqDlDgAIlaBW56YHeqBa7lTnAOCtoIa5H1ChAwBQhAmZfbJCoAMAUISs4U3gQLXcafUA\nAMKKCh0AgCKYPiiOQAcAOCLoXVLTAz1QLXcAAMIqcBX6bM9KBwCgFKZX6IELdAAA3GB6oNNyBwAg\nAKjQAQDzFvQBcZJkG16hBzLQuY4OAHAaLXcAQOBRRHkvkBU6AABOM71CD2yg03YHADjJ9GvotNwB\nAPMWhkFxpgtshS5RpQMAnEPLHQAQaGGpzmm5AwCAgrq7u2VZlrq6uqa8/9BDD6mlpUWf/OQn9bOf\n/azgPAIf6GE5cwQAuCtrRxx5TTc0NKRMJqO+vj6NjY1pcHAwN+3v/u7v9MILL6i3t1ff/OY3Cy5f\n4AMdAAAn2HbEkdd0yWRSHR0dkqT29nYNDAzkplVWXr4y/v777ysWixVcvlBcQ2dwHADAKxffOquL\nb6XyTk+lUlq5cqUkKRqNamRkZMr0v/mbv9HTTz+tf/qnfyr4OQUr9NOnT2vdunWqrq5WNpvNfVhr\na6va2tqUSl1ewAMHDqipqUlbtmzRuXPnJEm9vb1qbGxUW1ubTp06NcvXBQDAbHY2UtJr8R/VKbpu\nZe41XTQaVTqdliSNjo5eVYk//vjj+s1vfqMHHnig4PIVDPS6ujr19vbq1ltvzb23du1aHTlyRL29\nvYrFYrp06ZJ6enrU39+vbdu2qaenR5K0e/duPf/889q3b5/27t1b3NpyEdfSAcB5YTq2utVyj8fj\nSiQSkqREIqF4PJ6b9sEHH0iSqqqqcoV1PgUDffHixVedKbz++uuyLEv333+/JOnEiRNas2aNKioq\ncr3/CxcuqLq6WkuXLtX69euvah8AAIDL6uvrVVVVJcuyVFlZqYaGBnV2dkqSurq61NraqpaWFn3l\nK18pOJ85X0M/efKkYrGY7r77bh06dEjXXXedamtrJUm1tbVKpVJKpVK59yRpYmJirh/jCq6lAwBK\n5ebv0B999NEp//3YY49Jkv7xH/+x6HnMeZT7ZMV+++23a3h4eErvP51OKxaLTXlPkhYsWDDXj3FN\nmNpDAOCmsB1PS72GPv3llqID3bZtvf/++7lq+8UXX9SqVau0evVqDQ8PK5vN6vDhw4rH41qyZIku\nXLigTCajY8eO6aabbnLtCwAAgFla7uPj49q0aZOOHz+uTZs2ac+ePfrSl76ka665RitXrtSDDz6o\nSCSi7du3q7m5WXV1dTp48KAkadeuXdq4caOqq6v11FNPleXLFIvWOwDMT9iqc8n8W79GbI+SzeuN\ngUAHgNJ5cQz38rgdiUT0x//r047M6/T3n3Xlu4T2TnFen1AAAOCkUNwpDgD8wg+XBENbEBnecifQ\nAQBFC22YS7IL39fFc6FtuQOAiUyuzsMc5n4Q6grdD60tADABYW7+KPdQBzoAAEVz8aYwTqDlDgBA\nAIS+QqftDgCF0W6/jJY7APjI9PDihB85jHIHAH8gzOFnVOgAAo+WcelYd1cwfFAcgS6uowNBRiCV\njnU3lekxQaADCJzZgogT+NkR5v5DoP8BVTrgf8WEEPs5SkbL3T8IdSDY2L+LQ3WeB6Pc/YUNGfAn\n9l1nsB79iwodgK8VG0BU57MjzGfBjWUAwDmEjjtYr0UwvOVOoM+Aa+lAsLA/wxGGBzrX0AH4RilV\nJGGOsKBCz4MqHfA/9uHi0G4vkuEVOoFewORGzkEBKI2X90Znvy0OYT4Hhm9StNwBOC4SicwYFIQH\n4B4q9CJQqQPFmy20S72cxcmA81inc8Sd4oKD6+qAM+a6LxE8zmOdliBr9vGfQJ8jQh0obHL/KKZS\nL/T/FzMPAB/iGjoAV9i2nXvNxeT1d8LcPazb0kSyzrxm0t3dLcuy1NXVNeX9b3zjG2psbFRjY6N6\ne3sLLh+BXgIONsDclBru8/k8wHG2Q69phoaGlMlk1NfXp7GxMQ0ODuamff7zn9dLL72kn//85/rG\nN75RcPEIdABlRdh6i2LEPMlkUh0dHZKk9vZ2DQwM5KbdeOONkqRFixbN+ndHoM8DOwZQGrdCvZxd\nAD/imDVPWYde06RSKdXU1EiSotGoUqnUVX/m61//uu6+++6Ci8eguHlikBzgPfbB2RHmDijxTnFj\no+9qbPTdvNOj0ajS6bQkaXR0VLFYbMr0p59+WmfPntUdd9xR8HOo0B3AjgKUz5XX46nI4QeLotfp\nmv/x57nXdPF4XIlEQpKUSCQUj8dz01577TU9/vjj+s53vjPr5xDoDiHUAfcR3qXh+OSQrO3Ma5r6\n+npVVVXJsixVVlaqoaFBnZ2dkqSdO3fqzJkz+vSnP63bb7+94OJFbI/2kKBuYBxwgOLM9RjAvlWa\nIB1rvdwGIpGIVnxyiyPzejt5yJXvQoXuMH7SBsyOfaQ8WM/hwqA4AAggwtwFhneJCHSXMPodmBlB\nA9/ieejhRagDBLgXWOfhRKADQIAQ5u6J8LS1cKNKR5g5FS7sQ4UR4mVi+HbIKPcyYGcDALiNCr1M\nqNSB0gR5v5k82Z/Pd6RgKCNa7phEqCNsbNueV+CEZX8pNdgJ8/LiGjqmINQRNoW290KB5Jf95Mrv\nMNdlnvzzBDOcQKB7wIk2GxAE863gvTZ92Uvdt0s5Fvh5vfmW4cdsAt1D8zmzB4LC76HuBdaXR2i5\nA0Bhfg11L07E/bieUB4EuiG4tg6gEILcAFToKBahDnyI/eEygtwcEcO3R24sYxh2XgBAKajQDURl\ngrDJdyLLvgCj0HJHKfhpGwAYxjb7+am03AHAUFyCw1xQoRuOliO8MlOYuLUt+vVnawgZWu6YL0Id\n5VQoWN3cFtnGp+IEx0C03OEEdm6Ygm3RfaxjlIIK3UcYKAdT5NsWpwcR2yqChKetwXG04OGmuVzP\nnu3Psa1eNpf1QHVuMMNb7gS6T1Gtwy/CGuqlBDNhjvkg0H2OYIfpwrZtlnqTHMLcB2i5AwirMIV5\nMYE80wk4Qe4jtNxRDmFtawJ+RIjDDQR6gBDqcAJhA+RBhY5y4po65oMwB/KzDb+Gzo1lAooDM+aK\nbQbwTnd3tyzLUldX15T3n3zySa1cuVLbtm2bdR4EeoBFIpGrXgifmbYDtg2gBHbWmdc0Q0NDymQy\n6uvr09jYmAYHB3PTPvvZz+r5558vavEI9JDhwO2tcoeo13/fXn9+OXGZKwRcCvRkMqmOjg5JUnt7\nuwYGBnLTrr32Wi1YsKCoxSPQQ4iKzDu2bZflJ0tu/h1Pfofp3yXfnwFQWCqVUk1NjSQpGo0qlUqV\nNB8GxYUWHezKAAAJIUlEQVQYo+K9M7neJ4N3tr8HE07ACoU3gzERCtnSRrl/MHZeH1w6n3d6NBpV\nOp2WJI2OjioWi02ZXuz+T4UeciYERZhND/ZSrm9Pr5jdqIxnmx/VOMLAtrMlvRYtXKKaJctzr+ni\n8bgSiYQkKZFIKB6PT/vc4vYtAh3wsUI7eqkB6/bJQVBxcoxS1dfXq6qqSpZlqbKyUg0NDers7JQk\n/eQnP9G2bduUSCT0uc99ruB8IrZHeysbv5k4eHujlP2hlKd38ffrHo5p7vNy+41EIvrjupscmdfp\n90Zc+S5cQwcMcmULfrY/M9d5wj2EeUhwpzj4CYObvDF9fbP+/YEgh0kIdMyIYAemIrxhU6HDz7j+\nChDm+AMCHUFB1Y6wIcjhJwQ65oxgdwY39jEbYY7paLkjsAik0hAU5uLvBoUQ6Ag0rrGXjvVlFsIc\nfkegwzGEe3FYN2YgwDFnVOgII8IdJiPMUQrTW+4F7+WeTCbV1NSk5uZm7dixQ5K0f/9+NTc3a+vW\nrRofH5ckHThwQE1NTdqyZYvOnTsnSert7VVjY6Pa2tp06tQpl78GTMbjWmEStkWUqtSHs0x/uaVg\noN944406cuSI+vv7debMGfX19emFF15Qf3+/1q5dq2eeeUaXLl1ST0+P+vv7tW3bNvX09EiSdu/e\nreeff1779u3T3r17XfsC8I9inx4GuIFtD0FXMNBXrFihRYsWSZIWLlyokZERtbS0SJLa29s1MDCg\nkydPas2aNaqoqMi9d+HCBVVXV2vp0qVav369RkZGXP8i8BcOrCgHTiLhJNMr9KKuob/22mt65513\nFIvFVFFx+RygtrZWqVRKqVRKtbW1ed+TpImJCRcWHX7H79nhNIIbbvL1NXRJeu+993TvvffqySef\nVDQaVTqdliSl02nFYrFZ35OkBQsWuLT4CIIrqygOyABQmoKBPj4+rq1bt+rhhx/W8uXL1dDQoKNH\nj0qSDh8+rHg8rtWrV2t4eFjZbDb33pIlS3ThwgVlMhkdO3ZMN93kzDNkEQ4EO0rBNgO32co68nJL\nwZb7j370Iw0ODmrnzp2SpL1798qyLDU3N+uGG27Qjh07VFlZqe3bt6u5uVl1dXU6ePCgJGnXrl3a\nuHGjqqur9dRTT7n2BRBc+Q7QtOhxJYIc5WJ6yz1ie3R0ZCfEfBDqmMSxJDy83O8jkYhii691ZF6p\nD37vynfhxjLwpflU78UGACcNZiPIUW6mV+gEOgLFyYP89HkR8OVHaMMkBDoQEDxdzn0EOFA6Ah2Y\nA3477wyCG36Uldn7PYEOlIBqvTQEOfyMljsQUDxRbnYEOFA+BDrgAFrxUxHkCCIqdCBEwhDs0y83\nEN4Ii6yLd3lzAoEOuCBo19inhzYhDpiHQAdcUuzv2E2+Fk9wAx+i5Q5AUnHhONOfKVfIE95AYfxs\nDcC8uNm+J8SB4Jj1eegAvFfqI2WnP2ueZ88DpcvaWUdeM+nu7pZlWerq6pry/ltvvaW2tjY1NTUp\nkUgUXD4CHfCR2QKawAbcYzv0z3RDQ0PKZDLq6+vT2NiYBgcHc9P27dunPXv26LnnntPu3bsLLh+B\nDgCAh5LJpDo6OiRJ7e3tGhgYyE0bHh5WPB7X0qVLVVNTo3PnzuWdD4EOAEAR3Gq5p1Ip1dTUSJKi\n0ahSqVRu2sTERO7fp0+bjkAHAMBD0WhU6XRakjQ6OqpYLJabVlHxYUyn02ktW7Ys73w8G+Vu2u9t\nAQDIx83Misfj6unp0ec+9zklEgnddddduWlr167Vr371K61Zs0bpdFrXXHNN3vlQoQMA4KH6+npV\nVVXJsixVVlaqoaFBnZ2dkqSdO3dq165d2rhxo3bt2lVwPhGbUhkAAN+jQgcAIAA8C/R8P6JHad58\n802tWLFCra2t2rRpkyRp//79am5u1tatWzU+Pi5JOnDggJqamrRly5aCP3/AVKdPn9a6detUXV2t\nbPbyKNVi129vb68aGxvV1tamU6dOefYd/GCm9RyNRtXa2qq2trbcCF/W8/wkk0k1NTWpublZO3bs\nkMT2HAi2B1555RV7+/bttm3b9pe+9CX75Zdf9mIxAuW3v/2tvXXr1tx/v/322/ZnPvMZ27Zt+6GH\nHrJ/9KMf2WNjY3Zzc7M9MTFh/8u//Iu9f/9+rxbXdy5evGifPXvWbmlpsScmJua0fltbW+3z58/b\nyWTSvueee7z8Gsabvp5t27Y3bNgw5c+wnufvv//7v+0PPvjAtm3bvvPOO+2jR4+yPQeAJxV6oR/R\no3RHjhyRZVl69NFH9corr6ilpUXSh+v45MmTWrNmjSoqKljvc7R48eLcT0ls29bg4GBR6/fChQuq\nrq7W0qVLtX79eo2MjHj4Lcx35Xqe9Prrr8uyLN1///2SpBMnTrCe52nFihVatGiRJGnhwoUaGRlh\new4ATwK90I/oUZrrr79eJ06c0JEjR3T48GENDg6qtrZWklRbW6tUKqVUKnXVeyjN6OhoUev3yvek\nqTeJQHFOnjypvr4+nT17VocOHZp13Uus52K99tpreueddxSLxdieA8CTQC/0I3qUZtGiRaqurtaC\nBQu0efNm/emf/mluHafTacVisSnrffI9zF0kEplxXc72niQtWLDAk2X2s8nt9Pbbb9fw8DDr2SHv\nvfee7r33Xj355JNszwHhSaDH4/HcU2MSiYTi8bgXixEo58+fz/37L3/5S61atUpHjx6VJB0+fFjx\neFyrV6/W8PCwstls7j3MnW3bamhoKGr9LlmyRBcuXFAmk9GxY8d00003ebz0/mHbtt5///1cFfji\niy9q1apVrGcHjI+Pa+vWrXr44Ye1fPlytueg8Ori/Ze//GW7ubnZ7uzs9GoRAuVnP/uZ/fGPf9xu\nbGy0//7v/9627cuDWzZs2GDfeeed9qVLl2zbtu0f/OAHdmNjo71582Y7nU57uci+cunSJfu2226z\nly1bZre3t9vJZLLo9Xv48GE7Ho/bbW1t9u9+9zsvv4bxZlrP69atsy3Lsr/whS/Y2WzWtm3W83wd\nPHjQ/shHPmK3tLTYLS0t9sDAANtzAHBjGQAAAoAbywAAEAAEOgAAAUCgAwAQAAQ6AAABQKADABAA\nBDoAAAFAoAMAEAD/H3ensXrw3QqFAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }