rsignell-usgs/wms_test.ipynb

## wms_test.ipynb
{
 "metadata": {
  "name": "WMS test with Basemap"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "from IPython.core.display import Image\nfrom owslib.wms import WebMapService\nfrom mpl_toolkits.basemap import Basemap",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 63
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# Mercator Basemap for Baja California\nlon_min = -118.8\nlon_max = -108.6\nlat_min = 22.15\nlat_max = 32.34\nm = Basemap(llcrnrlon=lon_min, urcrnrlat=lat_max, urcrnrlon=lon_max, llcrnrlat=lat_min,\n        resolution='l', projection='merc')",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 69
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "urlbase='http://motherlode.ucar.edu:8080/thredds/wms/fmrc/NCEP/NAM/CONUS_12km/NCEP-NAM-CONUS_12km-noaaport_best.ncd?'\nwms = WebMapService(urlbase)\nimg = wms.getmap(layers=['Temperature_height_above_ground'],styles=['boxfill/rainbow'],time='2012-10-16T00:00:00.000Z',transparent=True,\n            srs='EPSG:4326',bbox=(lon_min,lat_min,lon_max,lat_max),size=(256,256),format='image/png',elevation='2')",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 65
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# add on the additional custom WMS arguments supported by ncWMS to define color range, colorbands and scale type \nurl = img.url + '&COLORSCALERANGE=271.2%2C308&NUMCOLORBANDS=20&LOGSCALE=false'",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 66
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "Image(url=url,format=u'png')",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "html": "<img src=\"http://motherlode.ucar.edu:8080/thredds/wms/fmrc/NCEP/NAM/CONUS_12km/NCEP-NAM-CONUS_12km-noaaport_best.ncd?layers=Temperature_height_above_ground&styles=boxfill%2Frainbow&elevation=2&srs=EPSG%3A4326&format=image%2Fpng&request=GetMap&bgcolor=0xFFFFFF&height=256&width=256&version=1.1.1&bbox=-118.8%2C22.15%2C-108.6%2C32.34&time=2012-10-16T00%3A00%3A00.000Z&exceptions=application%2Fvnd.ogc.se_xml&transparent=TRUE&COLORSCALERANGE=271.2%2C308&NUMCOLORBANDS=20&LOGSCALE=false\"/>",
       "output_type": "pyout",
       "prompt_number": 67,
       "text": "<IPython.core.display.Image at 0x3848190>"
      }
     ],
     "prompt_number": 67
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "# What's wrong with my syntax on warpimage?\nm.warpimage(url)\nm.drawcoastlines(linewidth=0.25)",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "pyout",
       "prompt_number": 70,
       "text": "<matplotlib.collections.LineCollection at 0x36e6250>"
      },
      {
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAANQAAADtCAYAAADHq1t9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGe9JREFUeJztnWlTVFmax/83M9k3QTYBRVBkU0QBV8SltKuruqqruiN6\ntn4x0dEfYd7Oi/4E8wl6JnomJmJ6pqM7qrSWrnItBTdcQXEF90IRV2RPyHnxN5tkyf2ee+7NfH4R\nhJhknvOg+c9zznOexfD5fD4IgmAKLt0GCEIiIYISBBMRQQmCiYigBMFERFCCYCIiKEEwEU+oHxqG\nYZUdguAogt02hRQUAPjOmG5L1Dx8CvT2A5trgRWFuq0Rkh1je/CfOWLLt6oU+NlO4HyfbksEITSO\nEJSf2lXAk+e6rRCE4DhKUOtWAbcf6rZCEILjKEG5XIBEHgp2xlGCEgS74zhBLc+TbZ9gXxwnqI01\ngHcGOHxetn+C/XCcoACgoQrYsQE41Ak8f6XbGkGYw5GCAoCsDODnu4A7j4CLN3VbQ0ZGgesDuq0Q\ndBI2UsLu7GgCjl/UbQVJTWFEx9MX/HtJAdBYDUgEV/LgaEF5vcCl28DElG5LSFoqULQM+KCNf3/6\nAjh2AUjxAB2b9NomWINjt3yzs8C//RGoqwQ+ChFbpQO/s6R0OcXVWA188QMwaRPhC+pwrKBcLuC3\nnwKnrtjL29e0Fui5O/+x5XnAp+3A0QsSOpXoOFZQAN+oGWlcrexCZjowNrH4cbcb+HgH8HgIuHQr\n+Osv3wIOngIePVNno6AORwsK4Orkduu2Yo6LN4GWuuA/39oIFOQC351demV9NULv5cu3wLdngIlJ\ndbYK5uN4QU1O67ZgPtNeevtCsXoF0L4R+OIk8HZ06edsrAF+sgXo7AHOXjPfTkENjvbyATz4330M\nrK3QbUl0Z7msDODzDp6rcrOAvCw+7glYbd1uYH8bMPwa+KqTq1tRvrk2C+ZihKocaxiGLTJ2w9HV\nA1SWAhXFeu24eR/IyQLKi6J73cu33Nr5AJTkA54gH3PHLwJ7W+K1UogXY3vwFHjHb/kAYGcTcPMB\n8OKNXjuePI9eTADPVGVFfG0wMQFyQewEEkJQALdGZ3qB0XHdlgjJTMIICmDdie/O0TFgNe/GgOxM\ntXPMzvJDQy6I7UtCnKEC8XqBL08Bv9jNy1+r+OEyI+BTFLt5RseBy7fnRJWXDTTXAI+fA1PTLBMg\nqCXUGSrhBAUA45PAiUvWhiQd7Z6L4bOSN++AK3foJcxIA/qfAAe2qBd2MhNKUAn5z56Rxi8rKSsC\nfnzOP60kLxvYHRB4u6Yc+NMx4B9/Yq0dAkmoM1QgVsf31VXS06gbj4fFQO0U35hMJKygrMZOLm2P\nm2UCBOtJWEGlpljvDTMMYMYGb+T0VIkB1EXCCqowz/qL3uYa4Ord8M9TTXqqfZIuk42EFVRRPjBk\ncQGX/Fzg2Uv9hWPS0+wXNJwsJKyg8nOA56+BvnvWzvvRduD6PRaP0YWsUPpIWEEZBu9jUlOAb05b\nGz2xZzPvg3QhZyh9JKyg/KytoLCOdDPNwyqy0vWEQAGyQukk4QUFMGrgo+0MzTl33Zo5a1bq2/aJ\noPSRFILy01DFsCQrKCmgg0IHHo/cQ+kiqQRlJbovenXPn6wknaCSJSRnStzmWkg6QRlGcogqPTV4\nARhBHUknqMI8Fj2xghSPvpWifSPQeVXP3MlM0gmqvMi66q06PX2GAVSXAf0WXhUISSioZTnA16eB\n+4Pq5yrWEP4USN1qppQkwxbXLiSdoAyDF66Xb6u/ePX59Hvbtm+QQplWknSCAoCZWWb0qk4TH3rF\nVUonBbnA6IQUdrGKpBTUL3czCU9158OHz4CVmotvAsDezcDxS7qtSA6SUlCbalnHb8VyljhWlTc1\nOs5Ksrpxu1mRVlrpqCcpBeWnrAj4pB24eke3JerZVAtcua3bisQnqQXlx8r6fTppXsf+U4I6kuSt\npAdd6RvBKC8ChjXXf090RFAKaagCrvXrtmI+Hhs1p0tERFAKqShe7AjwLfgSEgsRlGLWVjD8KJiA\nrBaW7ovmREcEpZjqisjj+WTFcj4iKKiLIvALpLWebWgiea6IytmIoMCzzl/PqktZL8oHZn36Oyz6\nkTAkdYigADRWAz/dBvSa5JFb6ly0swk4HWaVCvV6s9i1kR8eghpEUAGodilvawRO96idIxxuN7vJ\nd2m2I1ERQVlI0fvI886retPTS5czRf6BBTlhyYYIKgCPO74M20i2ajuagNY64NYD4NiF0FVtVW79\nWuqAawNSzMVsErIlaDxcH6BzYm9L9Hc2sbz5vV7gYCcLcQbruqjq6mhmhtnLP9+laIIEJVRLUFmh\nFtBYzU/vg6fYvzZSYl1JPB7g8w7g+3PB51O1SrndwJaG8C59IXJEUEuQl81P7Uu3uGKpxuXifKd7\ngSGLq82WLmfm8sOn1s6bqIiggmAY3PalpnD1UF3oxDC47bs2ADx6pnauhbTWAz135TxlBiKoMNSs\nZI27g6dC1/Mz65yzr5Wp84HOESvC7366Te6nzEAEFQGZ6cBnHSzJZUXW684mnqesTP3wePh7CvEh\ngoqC9o2s6/fXszzIL9wiGTBvNWmt55+XFBeSCcQt74a4UVxIK/FYvYJfo+PsNeW/Q9qz2fxU+g1r\ngNsPGV2xo8ncsQU1yD2UCYxPskPip+1zj8Xrwwhc6R4+ZffFfa1xDhqGI92ct7gAaKxKnlob0SL3\nUIrJSAN2bACOds89Fuv2b6nXrSrl/dixC7HbGAn724AP2lic87++lRLOsSCCMonleczOVVX2uKTA\nmjf4tJcexoYqqeMXCyIoE6lcAeRmzU8DiWaVCvfcimLg8VAsloXnwg1u+U73AptrgTcjwIDGTvZO\nRQRlMg1VwMTk/DejseDPQEL9bCHrVtFJYTY+H6Pf97cBuzfRfe5ysQa8EB0iKAW0NfByNjCMKJRg\nIl3FVBVYGRwGygoXzzUrgooaEZQi9mwGum8A78bmP24E+V4ndx4zIiSQ9dXW9NBKNERQCvl4B/D9\neaZoBBLPBfCKQuBHk50Fs7OMPA+kuCC6aHuBiKAUYhi8mzrUZZ6Hrq6SIVBm4fUGT3BMT+N2UIgc\nEZRiUjzAgTbgW5MuyP3nqLEJ3ksd6Q79/HCcuExHxFLUVQKPhjjHzEx88yQLIigLyM5k0uIPl80Z\nLz+HuVodzYArjoPYyCjgnQHSUhf/bGyCF9ZbGoD2Jm5drcgNczoiKIsoKWA3wws34h9rUy0DdT0e\nICtjseMjUrIzgaJlwOHzLBwTeNa7dIv3UQC3fh9t51xfdwGv3sb/OyQqEstnMb39LAZTv9qc8cYn\n2YWxKJ/ntKoyBu9Gy9gEozxmZoHKUl4gLxU76PMB/34I2NcCVJfHb78TCRXLJ4LSwNlrvPdZVWr+\n2Ee7GY8XDw8GgXfjjB9cyP98D+xqZtRGshJKUJK+oYFt63nQz8pgDKDdqAyxwrXU6a0paHfkDKWJ\n/W2s3jo+ae64ednA6xFzxwxk3SpuB822O1EQQWnkk50sdGmmS3rTOuCK4ibc+9vid9cnKiIojbhc\nFNWhTvMuft1u9TF4Lhc9gGZ4LBMNEZRm0lJZruz7c+aOq1pU5UXc9r0UF/o8xMtnE548ZzDqThNq\nR4xP0pPo/59tWgsULot/3KX48iSLdCZTq1FxmzsEM1zeC5mdZRHL3n7g1x+aXyfi7Shwvo/fV5cl\nx92UuM0dwOMhoKzI/HFdLqB5HfBqRE3RldwsOikAoP8xPxR8AJpr1K2KdkbOUDbh5gPzoid0saaC\nK2x7E3Djvm5r9CCCEkznzaj6bpB2RQRlAyYm2dBaZVWj9FTrOmyUFNB7ee9Ha+azEyIoG+CP5v7T\nMXVzbN/AYplWdaLfXMv5ki3rVwRlE+48irxLfKzsbaHHzyr2twHHLiZXcqIIyiYMDrP8sUqsvisy\nDODj7eZlKzsBEZQN+MPXzGf62Q71c1ldXjktlWXVOq9aO68uRFCaudYPTEwBQ69YdiwRKSkACnKT\nw5UugtLI5BTw4zCQkwksywZqV6mf0+3Wc6ZpqKJD5Nn74p8v3iSmF1AEpZE/n2CRlJExZsHWrVY/\n5+pS4IGmBtXtGxkCdbSbkSG3HwYvYeZUJPTIYvxHmPFJhhp1NLMakoqjzVI+iJUlwMkr+mLu9gfE\nKtav5u/+5h23hIH4fHysZiUzm52CCCoOHj4FzlxjPyWAl6cNVcyaDcbgMHDxJuPqPlbshAgUqV9c\nLpd9+j6lpgDZGUBhHuMNA5mc4gru87HKk1MQQcVBRTHLcO1t4d/HJ4G+e4zA9r9pfT7mDlVXACcv\ns0fvJ+3AiUuLx/vjYeBX+xaXRTYDH+xTSz2QbeuXfvzbM8Df7WOpNCfhMHPthcs1/wyQkcYiJoHM\n+rgqdb3vk5uRxscnp9idozigkVpeNhtiN1SxHJjZBC5Mw6/tHQ2el+3MHCtxSsTJ+jXAwVPBi5YY\nBs9KezbPiQkAPtzGA/qT58CqEuD4RYrsZzuBq3fUFkHZ08LIjK4edXPEy+Q08I0DL4QlwdAEZmdZ\nfbWydM5TF+kx5eRlbh1nZufax0xNU1iqef2Wompr4F2RXRgZZRjWh9t0W7I0krFrETfv0yV9YAtg\nRLH2n73GS8/ffKLKsqXx76jO93FF7Gi2xzZr6CUvutev0W3J0oigLGR8EjjczWjr8igycGdn1WTU\nRoIBOlJOXqZHLRq7VeDz0Wnjd/bYjVCCkjOUyWSksSfU4HB08Wu6xOQnN4vex2v94Z+rGjuskrEi\ngjIRH+bOTq317K/05Un7ly4O/KytKDa/Q2IsTE2z/5XTEEGZxFIbgMJlLLF16RbPV06gsRq4fk+3\nFexEn53pvA6KIijFGAZd5k9s8KkfCpsET/yNXc1s9nbplm5LokMEZRFZGcDouG4rIqOu0j4r6uZa\nhmo5BRGURbTUARcc8sZYWQI8tsGKOjPD+6inL3RbEjkiKJMI55hK8cxvuWk3FtpvQH199HC43cDa\nCmb92iWgNxwiKBMJJyqXS/+bdCmWsrulDvjyFL1tOtlYA2xfD3zVZf9zKCAXu0oI9g86Mgqc62MR\nSJ+Pd1aNVUBOlqXmLSLYB8G0l5e9PjCKIjXFSqsWc+EGc6f8zQ903VdJpITFRLo7GZtguseDp3Sv\np2iI/Y/kPRkorL2b1aSXRMr4JOv9Db+e/7jfm2oFIihNRCosrxc41AV83mHtp260U3VepStb90q1\nFP2P+W9nRSayhB7ZHI8H2N8KfGdy07VQxKLbqWl7iglgo4IBGxR9kQRDhfjftJGsVDlZQNMatgfN\nyeRjRcsYuRAN/3sE+Pv9kdmViPh8emMBRVAWYCAyUZUVze8RdfUOzwtrKyKbx+sF+p8wUTEtNbgt\nsfJunOe+zPQ4BlFIy/tL4NZ6fTbIls9Con0zb6xh0uHrkcief64P+OePly4oacQw/0I+bee2dGyC\nb1x/50K7kJ8LvNbcnEAEZRGxvpk/aGVu0MLilF4v3ciBZ+OxCeD6AMs6B85p1g7IMOg4+eY0I+hL\n8lljz06kpbA9kC5EUBZjRPnlel9w/7uzc48NPGbX+MoSPt57F3gzAvzhKzoOHg6aL6a/2W8A+Tn8\nvnIFy38d6rTPhfXWRuDsdX3zi6AcQFoqg0S7eiig8UnWnCguYF+pglxWS/rv33FbprLAix9/fb/l\necCBNkZV6FwZ/KSmUPQLg3tHx4H/O6q+vroIyiGsKARKC4COTYtrLVQUA/9wYC7rt2YlyxyrYnSC\nHd/9zQ3S04DPdnGlskPM3e5NvIr4umsuufPbM8DOJp79/PXVVSCCchBrKuaXIgvGqlLg4TN1drTW\nsclBbuZcvpLLxRXTLunraytYmff6AM+ghcsY7ZGdAbgVvuvFbZ6gzMyoK/xSVsRA1exMYHqGXRGb\n1po/T7wYBluhLuTweRYTVVGMRlaoBGVXM89VfYrS2dsagCu3maLufe+BTPWweKcdtn2hOLCFoUoq\n2ulILF+CM/CEotrZxHsa1Tx9AVwb4Kd//Wr188XDuevv+3JVRvc6ieVLYqrL6RHsu88qQqrd26XL\n2bLGCVm2Wxu5ir8bM29MEVQSYBhcobatB/583P5bMisZmwCOdJv3byKCSiIy04GfbgNOXdFtiX0o\nKWAe1fGL5owngkoycrLYmOCN4pg3w7BP9EQolucBp65yq3rrQfzjiaCSkD2b6ag42s0vFbUaKhXf\nhZlFz11Glvh8jOyPFxFUEuK/n/mgDdjXymiCo93m9ouqLJ1rz2NXevuZlPjL3cCjocguzcMhgkpy\nDIPu7Q/aGPT6wCQR2KmX71L4fMCjZ8DGtQxTaqwy5wNABCX8jYYq4IYJ5wgncLqX+Wb+Hr8rS4Bf\nfxj/uCIoYR7rVvJwrrsen0qmvXSXlxfNb5odLMs5GkRQwjyqy7n1+csJc5wVdtz2nbrCiHQVSHCs\nsAh/b9vvzsYXQFpWyFi/Ms0dERficqmr3iQrlBCU/Fzg5dvYX7+m3BxXtJMQQQlBaasHuuMoxOLx\n8BI5mRBBCUExDN7NjE3otsQ5iKCEkOxsootZiAwRlBASt5vdQmLF5WLhTbtxf1BNkLAISlBKexPT\nI+zGzfcX2D9cNico1o8ISlCK280IjJ67ui2ZIy+LRVzaN9Lp0tUj+VCCg6gqA4Ze2ce5samWgjIM\n4LMORkucMemcKIISLGFfC3D0gm4rFlOzkitoQS6TDA91xjeeREoIluByAZvW0RGwY4PeLohLUbea\n3UW8ATXkv+oEsjIY17i1EViWE34cEZQQlpkZdi9cW8HM1lipKOY2q6uHOViftJtnoxlkpAF3HvF3\nTU2hmPa2vP/9eyKL/xNBCWEZfMFUh3jE5Ke8iF/HLuhvjraQxmpGdizLZr9jf8lmt3tx95NgiKCE\nsFQUAxvWhH9eNNSvZuH+hipzx42XwAq4Kwrnvp/ysuFAVkbo14tTQgiLihSMFYXOqN3n50AbO0qG\nu1OTFUoQIsDtBnY0hX+erFBCSCan2BVQBakp9gxLigcRlBCSh8/YqVAFm2vn2uEkCiIoISRpKcB/\nHFJzjspMt6bbopWIoISQXLkD/Otv1Lm33W5gZFTN2DoQQQkhyc1Se1fU0QycuKxufKsRQQlhUVmj\n3DCAllrgwg11c1iJCEoIybQX+P1BtXOUFTGOTnUDAysQQQlBmZnhylFRrH6u3ZuY7Od0RFBCUN6M\n8q4oM139XIYBbGkAzsdRZckOiKCEoLx6C6x+fwf16Bnw5Um1lWBLlwMTk5zXqYighKDkZQMfbuX3\n//kNPX6dV9XOuasZOHtd7RwqEUEJQSlcBmRnsnpsZjoT7VSXFBsctubMpgrDF6w/PADDMOA7Y6U5\ngh0ZHGaKeIqH5bdmZpk6roLvzwE/2apmbLMwtgPBZCMrlBCWFYVs9eJysTvHvR/VzOP1Am6HvyMd\nbr6gg+V5wPNX5o97vm9+vyYnIoISosYfJT42we4anVfN6bKxegW9iU5GBCVEjWEA6WnA735PR8WW\nBnP605YVsX6Fk5GMXSEmOpr55Q+cLc5njYhVJeHrLiQyskIJMWEY86PQN6wBvjnN8KGhl/rs0o2s\nUIIpGAbwL//E709cYrBrdXn043jcDMhNceg7U1YowXT2bOZlcG9/9K9tq2fNvqUYeAJ88QNw+Lw9\nm2EDIihBEa31XG2iLcKfngasW8WSXQu5Pwh8vptlkb84yeqzdkMEJSijfjXDiA6fj+51VWXcMi68\n6yovAh4PMabw8w7eW928b5a15iChR4JyXr0FTl4BPm1ntEWknO4BJqeBkTG2w6lZyfPV/ra55/Td\nY8HMvS3WlXUOFXokghIsYWIS+OYM8PF2buui4Uwvz0w3HgC/2scVKpChlyx31lpvnr2hkFg+QTvp\nadymfX+eDoto2NoIjE4Av/10sZgAoLgAePHGHDvjRVYowXKOdAPrVgKrSs0b89lLRsXnZAJ997m9\nVIWsUIKt2N8G/DjMyAqzKCkAXo2wYVpOpnnjRosIStDCtvV883ebWENibwtQW2neeLEgghK0sWEN\ns4KPXzR33Gkv8HrE3DEjRQQlaKWqDGisYj9bs6If9rfRTT86bs540SCCErRTXMDt2l9OMB0kXgyD\nTonvzlnfLkcEJdiCrAzgs13A16fNqSBrGIwpvPkg/rGiQQQl2AaPh3dVp3uBByYkLKa46fiwEhGU\nYCsMA/hoO6vWft0VXwCsx8OwJSuRi13Btvh8rFcxOQ3s2USBRMulW4wDbK2jR9EMJJZPcDSTU0xa\nXJbDMKRY6O4Dht+w/sXyvPjsEUEJCcHAE2BiCmioiu31Ph9w8BTwWUd8dkjokZAQVJfHV2bMMDjG\nM4U1L0RQgqMoyo+vCMz6auBaDKn5kSKCEhzFpnXA5duxv151EqIISnAUhsE66/FEQESTNRz12OqG\nFoT4Werov7MJ6OqJcTyf2opJIijBtvgC/gzUQIqHEeWxCOP+IANyVSGCEmzHQgEt9XhbAxtqR8v1\ngbk2pyoQQQm2IpJFxwc2gHsVZc7TkW7eYal0TIigBFsQbFUK9fyqMqA/gjY6015WnG2pja08dDSI\noATHUrMyfF+ql2/ZxOCTnUB+rnqbHFqSXUgU4nG4+QDkZjPdfVnO4p/3P6YTIt5Qo2iQFUpwNFsb\ngXNLFHq5eJMpIB+0Lf6ZSmSFEhyNYbDRtdfLZMK3o8CVO9wOqnSPB7VHos0FHZh5tzo5xYvezDQg\nLxuoLFXbRVHSNwRbobK1kxX9AiR9QxAsQgQlJBS6GxuKoARLUf2Gt6hFVFBEUIJgIiIoQTAREZSQ\nUMgZSkgqdJ9xVCOCEixF9wqiGhGUYCmyQgmCiYjbXBCEiBFBCYKJiKCEhEK300MEJViK7jOOasKm\nbwiCsJhgsgmZsRtCa4IgLIFs+QTBRERQgmAiIihBMBERlCCYiAhKEEzk/wFtqhWi1+PJ2gAAAABJ\nRU5ErkJggg==\n",
       "text": "<matplotlib.figure.Figure at 0x36e1b10>"
      }
     ],
     "prompt_number": 70
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 68
    }
   ],
   "metadata": {}
  }
 ]
}
	{
	"metadata": {
	"name": "WMS test with Basemap"
	},
	"nbformat": 3,
	"nbformat_minor": 0,
	"worksheets": [
	{
	"cells": [
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "from IPython.core.display import Image\nfrom owslib.wms import WebMapService\nfrom mpl_toolkits.basemap import Basemap",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 63
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "# Mercator Basemap for Baja California\nlon_min = -118.8\nlon_max = -108.6\nlat_min = 22.15\nlat_max = 32.34\nm = Basemap(llcrnrlon=lon_min, urcrnrlat=lat_max, urcrnrlon=lon_max, llcrnrlat=lat_min,\n resolution='l', projection='merc')",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 69
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "urlbase='http://motherlode.ucar.edu:8080/thredds/wms/fmrc/NCEP/NAM/CONUS_12km/NCEP-NAM-CONUS_12km-noaaport_best.ncd?'\nwms = WebMapService(urlbase)\nimg = wms.getmap(layers=['Temperature_height_above_ground'],styles=['boxfill/rainbow'],time='2012-10-16T00:00:00.000Z',transparent=True,\n srs='EPSG:4326',bbox=(lon_min,lat_min,lon_max,lat_max),size=(256,256),format='image/png',elevation='2')",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 65
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "# add on the additional custom WMS arguments supported by ncWMS to define color range, colorbands and scale type \nurl = img.url + '&COLORSCALERANGE=271.2%2C308&NUMCOLORBANDS=20&LOGSCALE=false'",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 66
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "Image(url=url,format=u'png')",
	"language": "python",
	"metadata": {},
	"outputs": [
	{
	"html": "<img src=\"http://motherlode.ucar.edu:8080/thredds/wms/fmrc/NCEP/NAM/CONUS_12km/NCEP-NAM-CONUS_12km-noaaport_best.ncd?layers=Temperature_height_above_ground&styles=boxfill%2Frainbow&elevation=2&srs=EPSG%3A4326&format=image%2Fpng&request=GetMap&bgcolor=0xFFFFFF&height=256&width=256&version=1.1.1&bbox=-118.8%2C22.15%2C-108.6%2C32.34&time=2012-10-16T00%3A00%3A00.000Z&exceptions=application%2Fvnd.ogc.se_xml&transparent=TRUE&COLORSCALERANGE=271.2%2C308&NUMCOLORBANDS=20&LOGSCALE=false\"/>",
	"output_type": "pyout",
	"prompt_number": 67,
	"text": "<IPython.core.display.Image at 0x3848190>"
	}
	],
	"prompt_number": 67
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "# What's wrong with my syntax on warpimage?\nm.warpimage(url)\nm.drawcoastlines(linewidth=0.25)",
	"language": "python",
	"metadata": {},
	"outputs": [
	{
	"output_type": "pyout",
	"prompt_number": 70,
	"text": "<matplotlib.collections.LineCollection at 0x36e6250>"
	},
	{
	"output_type": "display_data",
	"png": "iVBORw0KGgoAAAANSUhEUgAAANQAAADtCAYAAADHq1t9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGe9JREFUeJztnWlTVFmax/83M9k3QTYBRVBkU0QBV8SltKuruqqruiN6\ntn4x0dEfYd7Oi/4E8wl6JnomJmJ6pqM7qrSWrnItBTdcQXEF90IRV2RPyHnxN5tkyf2ee+7NfH4R\nhJhknvOg+c9zznOexfD5fD4IgmAKLt0GCEIiIYISBBMRQQmCiYigBMFERFCCYCIiKEEwEU+oHxqG\nYZUdguAogt02hRQUAPjOmG5L1Dx8CvT2A5trgRWFuq0Rkh1je/CfOWLLt6oU+NlO4HyfbksEITSO\nEJSf2lXAk+e6rRCE4DhKUOtWAbcf6rZCEILjKEG5XIBEHgp2xlGCEgS74zhBLc+TbZ9gXxwnqI01\ngHcGOHxetn+C/XCcoACgoQrYsQE41Ak8f6XbGkGYw5GCAoCsDODnu4A7j4CLN3VbQ0ZGgesDuq0Q\ndBI2UsLu7GgCjl/UbQVJTWFEx9MX/HtJAdBYDUgEV/LgaEF5vcCl28DElG5LSFoqULQM+KCNf3/6\nAjh2AUjxAB2b9NomWINjt3yzs8C//RGoqwQ+ChFbpQO/s6R0OcXVWA188QMwaRPhC+pwrKBcLuC3\nnwKnrtjL29e0Fui5O/+x5XnAp+3A0QsSOpXoOFZQAN+oGWlcrexCZjowNrH4cbcb+HgH8HgIuHQr\n+Osv3wIOngIePVNno6AORwsK4Orkduu2Yo6LN4GWuuA/39oIFOQC351demV9NULv5cu3wLdngIlJ\ndbYK5uN4QU1O67ZgPtNeevtCsXoF0L4R+OIk8HZ06edsrAF+sgXo7AHOXjPfTkENjvbyATz4330M\nrK3QbUl0Z7msDODzDp6rcrOAvCw+7glYbd1uYH8bMPwa+KqTq1tRvrk2C+ZihKocaxiGLTJ2w9HV\nA1SWAhXFeu24eR/IyQLKi6J73cu33Nr5AJTkA54gH3PHLwJ7W+K1UogXY3vwFHjHb/kAYGcTcPMB\n8OKNXjuePI9eTADPVGVFfG0wMQFyQewEEkJQALdGZ3qB0XHdlgjJTMIICmDdie/O0TFgNe/GgOxM\ntXPMzvJDQy6I7UtCnKEC8XqBL08Bv9jNy1+r+OEyI+BTFLt5RseBy7fnRJWXDTTXAI+fA1PTLBMg\nqCXUGSrhBAUA45PAiUvWhiQd7Z6L4bOSN++AK3foJcxIA/qfAAe2qBd2MhNKUAn5z56Rxi8rKSsC\nfnzOP60kLxvYHRB4u6Yc+NMx4B9/Yq0dAkmoM1QgVsf31VXS06gbj4fFQO0U35hMJKygrMZOLm2P\nm2UCBOtJWEGlpljvDTMMYMYGb+T0VIkB1EXCCqowz/qL3uYa4Ord8M9TTXqqfZIuk42EFVRRPjBk\ncQGX/Fzg2Uv9hWPS0+wXNJwsJKyg8nOA56+BvnvWzvvRduD6PRaP0YWsUPpIWEEZBu9jUlOAb05b\nGz2xZzPvg3QhZyh9JKyg/KytoLCOdDPNwyqy0vWEQAGyQukk4QUFMGrgo+0MzTl33Zo5a1bq2/aJ\noPSRFILy01DFsCQrKCmgg0IHHo/cQ+kiqQRlJbovenXPn6wknaCSJSRnStzmWkg6QRlGcogqPTV4\nARhBHUknqMI8Fj2xghSPvpWifSPQeVXP3MlM0gmqvMi66q06PX2GAVSXAf0WXhUISSioZTnA16eB\n+4Pq5yrWEP4USN1qppQkwxbXLiSdoAyDF66Xb6u/ePX59Hvbtm+QQplWknSCAoCZWWb0qk4TH3rF\nVUonBbnA6IQUdrGKpBTUL3czCU9158OHz4CVmotvAsDezcDxS7qtSA6SUlCbalnHb8VyljhWlTc1\nOs5Ksrpxu1mRVlrpqCcpBeWnrAj4pB24eke3JerZVAtcua3bisQnqQXlx8r6fTppXsf+U4I6kuSt\npAdd6RvBKC8ChjXXf090RFAKaagCrvXrtmI+Hhs1p0tERFAKqShe7AjwLfgSEgsRlGLWVjD8KJiA\nrBaW7ovmREcEpZjqisjj+WTFcj4iKKiLIvALpLWebWgiea6IytmIoMCzzl/PqktZL8oHZn36Oyz6\nkTAkdYigADRWAz/dBvSa5JFb6ly0swk4HWaVCvV6s9i1kR8eghpEUAGodilvawRO96idIxxuN7vJ\nd2m2I1ERQVlI0fvI886retPTS5czRf6BBTlhyYYIKgCPO74M20i2ajuagNY64NYD4NiF0FVtVW79\nWuqAawNSzMVsErIlaDxcH6BzYm9L9Hc2sbz5vV7gYCcLcQbruqjq6mhmhtnLP9+laIIEJVRLUFmh\nFtBYzU/vg6fYvzZSYl1JPB7g8w7g+3PB51O1SrndwJaG8C59IXJEUEuQl81P7Uu3uGKpxuXifKd7\ngSGLq82WLmfm8sOn1s6bqIiggmAY3PalpnD1UF3oxDC47bs2ADx6pnauhbTWAz135TxlBiKoMNSs\nZI27g6dC1/Mz65yzr5Wp84HOESvC7366Te6nzEAEFQGZ6cBnHSzJZUXW684mnqesTP3wePh7CvEh\ngoqC9o2s6/fXszzIL9wiGTBvNWmt55+XFBeSCcQt74a4UVxIK/FYvYJfo+PsNeW/Q9qz2fxU+g1r\ngNsPGV2xo8ncsQU1yD2UCYxPskPip+1zj8Xrwwhc6R4+ZffFfa1xDhqGI92ct7gAaKxKnlob0SL3\nUIrJSAN2bACOds89Fuv2b6nXrSrl/dixC7HbGAn724AP2lic87++lRLOsSCCMonleczOVVX2uKTA\nmjf4tJcexoYqqeMXCyIoE6lcAeRmzU8DiWaVCvfcimLg8VAsloXnwg1u+U73AptrgTcjwIDGTvZO\nRQRlMg1VwMTk/DejseDPQEL9bCHrVtFJYTY+H6Pf97cBuzfRfe5ysQa8EB0iKAW0NfByNjCMKJRg\nIl3FVBVYGRwGygoXzzUrgooaEZQi9mwGum8A78bmP24E+V4ndx4zIiSQ9dXW9NBKNERQCvl4B/D9\neaZoBBLPBfCKQuBHk50Fs7OMPA+kuCC6aHuBiKAUYhi8mzrUZZ6Hrq6SIVBm4fUGT3BMT+N2UIgc\nEZRiUjzAgTbgW5MuyP3nqLEJ3ksd6Q79/HCcuExHxFLUVQKPhjjHzEx88yQLIigLyM5k0uIPl80Z\nLz+HuVodzYArjoPYyCjgnQHSUhf/bGyCF9ZbGoD2Jm5drcgNczoiKIsoKWA3wws34h9rUy0DdT0e\nICtjseMjUrIzgaJlwOHzLBwTeNa7dIv3UQC3fh9t51xfdwGv3sb/OyQqEstnMb39LAZTv9qc8cYn\n2YWxKJ/ntKoyBu9Gy9gEozxmZoHKUl4gLxU76PMB/34I2NcCVJfHb78TCRXLJ4LSwNlrvPdZVWr+\n2Ee7GY8XDw8GgXfjjB9cyP98D+xqZtRGshJKUJK+oYFt63nQz8pgDKDdqAyxwrXU6a0paHfkDKWJ\n/W2s3jo+ae64ednA6xFzxwxk3SpuB822O1EQQWnkk50sdGmmS3rTOuCK4ibc+9vid9cnKiIojbhc\nFNWhTvMuft1u9TF4Lhc9gGZ4LBMNEZRm0lJZruz7c+aOq1pU5UXc9r0UF/o8xMtnE548ZzDqThNq\nR4xP0pPo/59tWgsULot/3KX48iSLdCZTq1FxmzsEM1zeC5mdZRHL3n7g1x+aXyfi7Shwvo/fV5cl\nx92UuM0dwOMhoKzI/HFdLqB5HfBqRE3RldwsOikAoP8xPxR8AJpr1K2KdkbOUDbh5gPzoid0saaC\nK2x7E3Djvm5r9CCCEkznzaj6bpB2RQRlAyYm2dBaZVWj9FTrOmyUFNB7ee9Ha+azEyIoG+CP5v7T\nMXVzbN/AYplWdaLfXMv5ki3rVwRlE+48irxLfKzsbaHHzyr2twHHLiZXcqIIyiYMDrP8sUqsvisy\nDODj7eZlKzsBEZQN+MPXzGf62Q71c1ldXjktlWXVOq9aO68uRFCaudYPTEwBQ69YdiwRKSkACnKT\nw5UugtLI5BTw4zCQkwksywZqV6mf0+3Wc6ZpqKJD5Nn74p8v3iSmF1AEpZE/n2CRlJExZsHWrVY/\n5+pS4IGmBtXtGxkCdbSbkSG3HwYvYeZUJPTIYvxHmPFJhhp1NLMakoqjzVI+iJUlwMkr+mLu9gfE\nKtav5u/+5h23hIH4fHysZiUzm52CCCoOHj4FzlxjPyWAl6cNVcyaDcbgMHDxJuPqPlbshAgUqV9c\nLpd9+j6lpgDZGUBhHuMNA5mc4gru87HKk1MQQcVBRTHLcO1t4d/HJ4G+e4zA9r9pfT7mDlVXACcv\ns0fvJ+3AiUuLx/vjYeBX+xaXRTYDH+xTSz2QbeuXfvzbM8Df7WOpNCfhMHPthcs1/wyQkcYiJoHM\n+rgqdb3vk5uRxscnp9idozigkVpeNhtiN1SxHJjZBC5Mw6/tHQ2el+3MHCtxSsTJ+jXAwVPBi5YY\nBs9KezbPiQkAPtzGA/qT58CqEuD4RYrsZzuBq3fUFkHZ08LIjK4edXPEy+Q08I0DL4QlwdAEZmdZ\nfbWydM5TF+kx5eRlbh1nZufax0xNU1iqef2Wompr4F2RXRgZZRjWh9t0W7I0krFrETfv0yV9YAtg\nRLH2n73GS8/ffKLKsqXx76jO93FF7Gi2xzZr6CUvutev0W3J0oigLGR8EjjczWjr8igycGdn1WTU\nRoIBOlJOXqZHLRq7VeDz0Wnjd/bYjVCCkjOUyWSksSfU4HB08Wu6xOQnN4vex2v94Z+rGjuskrEi\ngjIRH+bOTq317K/05Un7ly4O/KytKDa/Q2IsTE2z/5XTEEGZxFIbgMJlLLF16RbPV06gsRq4fk+3\nFexEn53pvA6KIijFGAZd5k9s8KkfCpsET/yNXc1s9nbplm5LokMEZRFZGcDouG4rIqOu0j4r6uZa\nhmo5BRGURbTUARcc8sZYWQI8tsGKOjPD+6inL3RbEjkiKJMI55hK8cxvuWk3FtpvQH199HC43cDa\nCmb92iWgNxwiKBMJJyqXS/+bdCmWsrulDvjyFL1tOtlYA2xfD3zVZf9zKCAXu0oI9g86Mgqc62MR\nSJ+Pd1aNVUBOlqXmLSLYB8G0l5e9PjCKIjXFSqsWc+EGc6f8zQ903VdJpITFRLo7GZtguseDp3Sv\np2iI/Y/kPRkorL2b1aSXRMr4JOv9Db+e/7jfm2oFIihNRCosrxc41AV83mHtp260U3VepStb90q1\nFP2P+W9nRSayhB7ZHI8H2N8KfGdy07VQxKLbqWl7iglgo4IBGxR9kQRDhfjftJGsVDlZQNMatgfN\nyeRjRcsYuRAN/3sE+Pv9kdmViPh8emMBRVAWYCAyUZUVze8RdfUOzwtrKyKbx+sF+p8wUTEtNbgt\nsfJunOe+zPQ4BlFIy/tL4NZ6fTbIls9Con0zb6xh0uHrkcief64P+OePly4oacQw/0I+bee2dGyC\nb1x/50K7kJ8LvNbcnEAEZRGxvpk/aGVu0MLilF4v3ciBZ+OxCeD6AMs6B85p1g7IMOg4+eY0I+hL\n8lljz06kpbA9kC5EUBZjRPnlel9w/7uzc48NPGbX+MoSPt57F3gzAvzhKzoOHg6aL6a/2W8A+Tn8\nvnIFy38d6rTPhfXWRuDsdX3zi6AcQFoqg0S7eiig8UnWnCguYF+pglxWS/rv33FbprLAix9/fb/l\necCBNkZV6FwZ/KSmUPQLg3tHx4H/O6q+vroIyiGsKARKC4COTYtrLVQUA/9wYC7rt2YlyxyrYnSC\nHd/9zQ3S04DPdnGlskPM3e5NvIr4umsuufPbM8DOJp79/PXVVSCCchBrKuaXIgvGqlLg4TN1drTW\nsclBbuZcvpLLxRXTLunraytYmff6AM+ghcsY7ZGdAbgVvuvFbZ6gzMyoK/xSVsRA1exMYHqGXRGb\n1po/T7wYBluhLuTweRYTVVGMRlaoBGVXM89VfYrS2dsagCu3maLufe+BTPWweKcdtn2hOLCFoUoq\n2ulILF+CM/CEotrZxHsa1Tx9AVwb4Kd//Wr188XDuevv+3JVRvc6ieVLYqrL6RHsu88qQqrd26XL\n2bLGCVm2Wxu5ir8bM29MEVQSYBhcobatB/583P5bMisZmwCOdJv3byKCSiIy04GfbgNOXdFtiX0o\nKWAe1fGL5owngkoycrLYmOCN4pg3w7BP9EQolucBp65yq3rrQfzjiaCSkD2b6ag42s0vFbUaKhXf\nhZlFz11Glvh8jOyPFxFUEuK/n/mgDdjXymiCo93m9ouqLJ1rz2NXevuZlPjL3cCjocguzcMhgkpy\nDIPu7Q/aGPT6wCQR2KmX71L4fMCjZ8DGtQxTaqwy5wNABCX8jYYq4IYJ5wgncLqX+Wb+Hr8rS4Bf\nfxj/uCIoYR7rVvJwrrsen0qmvXSXlxfNb5odLMs5GkRQwjyqy7n1+csJc5wVdtz2nbrCiHQVSHCs\nsAh/b9vvzsYXQFpWyFi/Ms0dERficqmr3iQrlBCU/Fzg5dvYX7+m3BxXtJMQQQlBaasHuuMoxOLx\n8BI5mRBBCUExDN7NjE3otsQ5iKCEkOxsootZiAwRlBASt5vdQmLF5WLhTbtxf1BNkLAISlBKexPT\nI+zGzfcX2D9cNico1o8ISlCK280IjJ67ui2ZIy+LRVzaN9Lp0tUj+VCCg6gqA4Ze2ce5samWgjIM\n4LMORkucMemcKIISLGFfC3D0gm4rFlOzkitoQS6TDA91xjeeREoIluByAZvW0RGwY4PeLohLUbea\n3UW8ATXkv+oEsjIY17i1EViWE34cEZQQlpkZdi9cW8HM1lipKOY2q6uHOViftJtnoxlkpAF3HvF3\nTU2hmPa2vP/9eyKL/xNBCWEZfMFUh3jE5Ke8iF/HLuhvjraQxmpGdizLZr9jf8lmt3tx95NgiKCE\nsFQUAxvWhH9eNNSvZuH+hipzx42XwAq4Kwrnvp/ysuFAVkbo14tTQgiLihSMFYXOqN3n50AbO0qG\nu1OTFUoQIsDtBnY0hX+erFBCSCan2BVQBakp9gxLigcRlBCSh8/YqVAFm2vn2uEkCiIoISRpKcB/\nHFJzjspMt6bbopWIoISQXLkD/Otv1Lm33W5gZFTN2DoQQQkhyc1Se1fU0QycuKxufKsRQQlhUVmj\n3DCAllrgwg11c1iJCEoIybQX+P1BtXOUFTGOTnUDAysQQQlBmZnhylFRrH6u3ZuY7Od0RFBCUN6M\n8q4oM139XIYBbGkAzsdRZckOiKCEoLx6C6x+fwf16Bnw5Um1lWBLlwMTk5zXqYighKDkZQMfbuX3\n//kNPX6dV9XOuasZOHtd7RwqEUEJQSlcBmRnsnpsZjoT7VSXFBsctubMpgrDF6w/PADDMOA7Y6U5\ngh0ZHGaKeIqH5bdmZpk6roLvzwE/2apmbLMwtgPBZCMrlBCWFYVs9eJysTvHvR/VzOP1Am6HvyMd\nbr6gg+V5wPNX5o97vm9+vyYnIoISosYfJT42we4anVfN6bKxegW9iU5GBCVEjWEA6WnA735PR8WW\nBnP605YVsX6Fk5GMXSEmOpr55Q+cLc5njYhVJeHrLiQyskIJMWEY86PQN6wBvjnN8KGhl/rs0o2s\nUIIpGAbwL//E709cYrBrdXn043jcDMhNceg7U1YowXT2bOZlcG9/9K9tq2fNvqUYeAJ88QNw+Lw9\nm2EDIihBEa31XG2iLcKfngasW8WSXQu5Pwh8vptlkb84yeqzdkMEJSijfjXDiA6fj+51VWXcMi68\n6yovAh4PMabw8w7eW928b5a15iChR4JyXr0FTl4BPm1ntEWknO4BJqeBkTG2w6lZyfPV/ra55/Td\nY8HMvS3WlXUOFXokghIsYWIS+OYM8PF2buui4Uwvz0w3HgC/2scVKpChlyx31lpvnr2hkFg+QTvp\nadymfX+eDoto2NoIjE4Av/10sZgAoLgAePHGHDvjRVYowXKOdAPrVgKrSs0b89lLRsXnZAJ997m9\nVIWsUIKt2N8G/DjMyAqzKCkAXo2wYVpOpnnjRosIStDCtvV883ebWENibwtQW2neeLEgghK0sWEN\ns4KPXzR33Gkv8HrE3DEjRQQlaKWqDGisYj9bs6If9rfRTT86bs540SCCErRTXMDt2l9OMB0kXgyD\nTonvzlnfLkcEJdiCrAzgs13A16fNqSBrGIwpvPkg/rGiQQQl2AaPh3dVp3uBByYkLKa46fiwEhGU\nYCsMA/hoO6vWft0VXwCsx8OwJSuRi13Btvh8rFcxOQ3s2USBRMulW4wDbK2jR9EMJJZPcDSTU0xa\nXJbDMKRY6O4Dht+w/sXyvPjsEUEJCcHAE2BiCmioiu31Ph9w8BTwWUd8dkjokZAQVJfHV2bMMDjG\nM4U1L0RQgqMoyo+vCMz6auBaDKn5kSKCEhzFpnXA5duxv151EqIISnAUhsE66/FEQESTNRz12OqG\nFoT4Werov7MJ6OqJcTyf2opJIijBtvgC/gzUQIqHEeWxCOP+IANyVSGCEmzHQgEt9XhbAxtqR8v1\ngbk2pyoQQQm2IpJFxwc2gHsVZc7TkW7eYal0TIigBFsQbFUK9fyqMqA/gjY6015WnG2pja08dDSI\noATHUrMyfF+ql2/ZxOCTnUB+rnqbHFqSXUgU4nG4+QDkZjPdfVnO4p/3P6YTIt5Qo2iQFUpwNFsb\ngXNLFHq5eJMpIB+0Lf6ZSmSFEhyNYbDRtdfLZMK3o8CVO9wOqnSPB7VHos0FHZh5tzo5xYvezDQg\nLxuoLFXbRVHSNwRbobK1kxX9AiR9QxAsQgQlJBS6GxuKoARLUf2Gt6hFVFBEUIJgIiIoQTAREZSQ\nUMgZSkgqdJ9xVCOCEixF9wqiGhGUYCmyQgmCiYjbXBCEiBFBCYKJiKCEhEK300MEJViK7jOOasKm\nbwiCsJhgsgmZsRtCa4IgLIFs+QTBRERQgmAiIihBMBERlCCYiAhKEEzk/wFtqhWi1+PJ2gAAAABJ\nRU5ErkJggg==\n",
	"text": "<matplotlib.figure.Figure at 0x36e1b10>"
	}
	],
	"prompt_number": 70
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 68
	}
	],
	"metadata": {}
	}
	]
	}