GoogleMapsWidget-SciPyDemo-HelloWorld.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:6cae7b06dd254107ccbeea1a02dd8a5fb27fda402145bd4c2443b7f981685eb3"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%matplotlib inline\n",
      "\n",
      "from __future__ import print_function # For py 2.7 compat\n",
      "\n",
      "import datetime\n",
      "import ee\n",
      "import numpy as np\n",
      "import os\n",
      "import matplotlib.pyplot as plt\n",
      "from pylab import figure\n",
      "import matplotlib.dates as mdates\n",
      "\n",
      "from IPython.html import widgets\n",
      "from IPython.display import display\n",
      "from IPython.utils import traitlets\n",
      "from IPython.core.display import Javascript"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Create a Google Maps widget."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "class GoogleMapsWidget(widgets.DOMWidget):\n",
      "    _view_name = traitlets.Unicode('GoogleMapsView', sync=True)\n",
      "    value = traitlets.Unicode(sync=True)\n",
      "    description = traitlets.Unicode(sync=True)\n",
      "    lat = traitlets.CFloat(0, help=\"Center latitude, -90 to 90\", sync=True)\n",
      "    lng = traitlets.CFloat(0, help=\"Center longitude, -180 to 180\", sync=True)\n",
      "    zoom = traitlets.CInt(0, help=\"Zoom level, 0 to ~25\", sync=True)\n",
      "    bounds = traitlets.List([], help=\"Visible bounds, [W, S, E, N]\", sync=True)\n",
      "    \n",
      "    def __init__(self, lng=0.0, lat=0.0, zoom=2):\n",
      "        self.lng = lng\n",
      "        self.lat = lat\n",
      "        self.zoom = zoom\n",
      "        \n",
      "    def addLayer(self, image, vis_params, name=None, visible=True):\n",
      "        mapid = image.getMapId(vis_params)\n",
      "        self.send({'command':'addLayer', 'mapid': mapid['mapid'], 'token': mapid['token'], 'name': name, 'visible': visible})\n",
      "        \n",
      "    def center(self, lng, lat, zoom=None):\n",
      "        self.send({'command': 'center', 'lng': lng, 'lat': lat, 'zoom': zoom})"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 3,
     "metadata": {},
     "source": [
      "Add the Javascript code for the widget."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%javascript\n",
      "\n",
      "require([\"widgets/js/widget\"], function(WidgetManager){\n",
      "    var maps = [];\n",
      "    \n",
      "    // Define the GoogleMapsView\n",
      "    var GoogleMapsView = IPython.DOMWidgetView.extend({\n",
      "        \n",
      "        render: function() {\n",
      "            // Resize widget element to be 100% wide\n",
      "            this.$el.css('width', '100%');\n",
      "\n",
      "            // iframe source;  just enough to load Google Maps and let us poll whether initialization is complete\n",
      "            var src='<html style=\"height:100%\"><head>' +\n",
      "                '<scr'+'ipt src=\"http://maps.googleapis.com/maps/api/js?sensor=false\"></scr'+'ipt>' +\n",
      "                '<scr'+'ipt>google.maps.event.addDomListener(window,\"load\",function(){ready=true});</scr'+'ipt>' +\n",
      "                '</head>' +\n",
      "                '<body style=\"height:100%; margin:0px; padding:0px\"></body></html>';\n",
      "            \n",
      "            // Create the Google Maps container element.\n",
      "            this.$iframe = $('<iframe />')\n",
      "                .css('width', '100%')\n",
      "                .css('height', '300px')\n",
      "                .attr('srcdoc', src)\n",
      "                .appendTo(this.$el);\n",
      "                        \n",
      "            var self = this; // hold onto this for initMapWhenReady\n",
      "\n",
      "            // Wait until maps library has finished loading in iframe, then create map\n",
      "            function initMapWhenReady() {\n",
      "                // Iframe document and window\n",
      "                var doc = self.$iframe[0].contentDocument;\n",
      "                var win = self.$iframe[0].contentWindow;\n",
      "                if (!win || !win.ready) {\n",
      "                    // Maps library not yet loaded;  try again soon\n",
      "                    setTimeout(initMapWhenReady, 20);\n",
      "                    return;\n",
      "                }\n",
      "\n",
      "                // Maps library finished loading.  Build map now.\n",
      "                var mapOptions = {\n",
      "                    center: new win.google.maps.LatLng(self.model.get('lat'), self.model.get('lng')),\n",
      "                    zoom: self.model.get('zoom')\n",
      "                };\n",
      "                var mapDiv = $('<div />')\n",
      "                    .css('width', '100%')\n",
      "                    .css('height', '100%')\n",
      "                    .appendTo($(doc.body));\n",
      "                self.map = new win.google.maps.Map(mapDiv[0], mapOptions);\n",
      "                \n",
      "                \n",
      "                // Add an event listeners for user panning, zooming, and resizing map\n",
      "                // TODO(rsargent): Bind self across all methods, and save some plumbing here\n",
      "                win.google.maps.event.addListener(self.map, 'bounds_changed', function () {\n",
      "                    self.handleBoundsChanged();\n",
      "                });\n",
      "                \n",
      "                self.initializeLayersControl();\n",
      "            }\n",
      "            initMapWhenReady();\n",
      "        },\n",
      "        \n",
      "        LayersControl: function(widget, controlDiv, map) {\n",
      "            var win = widget.$iframe[0].contentWindow;\n",
      "            var chicago = new win.google.maps.LatLng(41.850033, -87.6500523);\n",
      "\n",
      "            // Set CSS styles for the DIV containing the control\n",
      "            // Setting padding to 5 px will offset the control\n",
      "            // from the edge of the map.\n",
      "            controlDiv.style.padding = '5px';\n",
      "\n",
      "            // Set CSS for the control border.\n",
      "            var $controlUI = $('<div />')\n",
      "                .css('backgroundColor', 'white')\n",
      "                .css('borderStyle', 'solid')\n",
      "                .css('borderWidth', '1px')\n",
      "                .css('cursor', 'pointer')\n",
      "                .css('textAlign', 'center')\n",
      "                .appendTo($(controlDiv));\n",
      "            \n",
      "            // Set CSS for the control interior.\n",
      "            var $controlContents = $('<div />')\n",
      "                .css('fontFamily', 'Arial,sans-serif')\n",
      "                .css('fontSize', '12px')\n",
      "                .css('paddingLeft', '4px')\n",
      "                .css('paddingRight', '4px')\n",
      "                .css('paddingTop', '0px')\n",
      "                .css('paddingBottom', '0px')\n",
      "                .appendTo($controlUI);\n",
      "            \n",
      "            this.$controlTable = $('<table />')\n",
      "                .append($('<tr><td colspan=2>Layers</td></tr>'))\n",
      "                .appendTo($controlContents);\n",
      "        },\n",
      "\n",
      "        initializeLayersControl: function() {\n",
      "            var doc = this.$iframe[0].contentDocument;\n",
      "            var win = this.$iframe[0].contentWindow;\n",
      "\n",
      "            // Create the DIV to hold the control and call the LayersControl() constructor\n",
      "            // passing in this DIV.\n",
      "    \n",
      "            var layersControlDiv = document.createElement('div');\n",
      "            this.layersControl = new this.LayersControl(this, layersControlDiv, this.map);\n",
      "\n",
      "            layersControlDiv.index = 1;\n",
      "            this.map.controls[win.google.maps.ControlPosition.TOP_RIGHT].push(layersControlDiv);\n",
      "        },\n",
      "        \n",
      "        // Map geometry changed (pan, zoom, resize)\n",
      "        handleBoundsChanged: function() {\n",
      "            this.model.set('lng', this.map.getCenter().lng());\n",
      "            this.model.set('lat', this.map.getCenter().lat());\n",
      "            this.model.set('zoom', this.map.getZoom());\n",
      "            var bounds = this.map.getBounds();\n",
      "            var playgroundCompatible = [bounds.getSouthWest().lng(), bounds.getSouthWest().lat(),\n",
      "                                        bounds.getNorthEast().lng(), bounds.getNorthEast().lat()];\n",
      "            this.model.set('bounds', playgroundCompatible);\n",
      "            this.touch();\n",
      "        },\n",
      "        \n",
      "        // Receive custom messages from Python backend\n",
      "        on_msg: function(msg) {\n",
      "            var win = this.$iframe[0].contentWindow;\n",
      "            if (msg.command == 'addLayer') {\n",
      "                this.addLayer(msg.mapid, msg.token, msg.name, msg.visible);\n",
      "            } else if (msg.command == 'center') {\n",
      "                this.map.setCenter(new win.google.maps.LatLng(msg.lat, msg.lng));\n",
      "                if (msg.zoom !== null) {\n",
      "                    this.map.setZoom(msg.zoom);\n",
      "                }\n",
      "            }\n",
      "        },\n",
      "        \n",
      "        // Add an Earth Engine layer\n",
      "        addLayer: function(mapid, token, name, visible) {\n",
      "            var win = this.$iframe[0].contentWindow;\n",
      "            var eeMapOptions = {\n",
      "                getTileUrl: function(tile, zoom) {\n",
      "                    var url = ['https://earthengine.googleapis.com/map',\n",
      "                               mapid, zoom, tile.x, tile.y].join(\"/\");\n",
      "                    url += '?token=' + token\n",
      "                    return url;\n",
      "                },\n",
      "                tileSize: new win.google.maps.Size(256, 256),\n",
      "                opacity: visible ? 1.0 : 0.0,\n",
      "            };\n",
      "            \n",
      "            // Create the overlay map type\n",
      "            var mapType = new win.google.maps.ImageMapType(eeMapOptions);\n",
      "                \n",
      "            // Overlay the Earth Engine generated layer\n",
      "            this.map.overlayMapTypes.push(mapType);\n",
      "\n",
      "            // Update layer visibility control\n",
      "            var maxSlider = 100;\n",
      "            \n",
      "            function updateOpacity() {\n",
      "                mapType.setOpacity($checkbox[0].checked ? $slider[0].value / 100.0 : 0);\n",
      "            }\n",
      "            \n",
      "            var $checkbox = $('<input type=\"checkbox\">')\n",
      "                .prop('checked', visible)\n",
      "                .change(updateOpacity);\n",
      "            \n",
      "            var $slider = $('<input type=\"range\" />')\n",
      "                .prop('min', 0)\n",
      "                .prop('max', maxSlider)\n",
      "                .prop('value', maxSlider)\n",
      "                .css('width', '60px')\n",
      "                .on('input', updateOpacity);\n",
      "\n",
      "            // If user doesn't specify a layer name, create a default\n",
      "            if (name === null) {\n",
      "                name = 'Layer ' + this.map.overlayMapTypes.length;\n",
      "            }\n",
      "            \n",
      "            var $row = $('<tr />');\n",
      "            $('<td align=\"left\" />').append($checkbox).append(name).appendTo($row);\n",
      "            $('<td />').append($slider).appendTo($row);\n",
      "\n",
      "            this.layersControl.$controlTable.append($row);\n",
      "        }\n",
      "    });\n",
      "    \n",
      "    // Register the GoogleMapsView with the widget manager.\n",
      "    WidgetManager.register_widget_view('GoogleMapsView', GoogleMapsView);\n",
      "});\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "javascript": [
        "\n",
        "require([\"widgets/js/widget\"], function(WidgetManager){\n",
        "    var maps = [];\n",
        "    \n",
        "    // Define the GoogleMapsView\n",
        "    var GoogleMapsView = IPython.DOMWidgetView.extend({\n",
        "        \n",
        "        render: function() {\n",
        "            // Resize widget element to be 100% wide\n",
        "            this.$el.css('width', '100%');\n",
        "\n",
        "            // iframe source;  just enough to load Google Maps and let us poll whether initialization is complete\n",
        "            var src='<html style=\"height:100%\"><head>' +\n",
        "                '<scr'+'ipt src=\"http://maps.googleapis.com/maps/api/js?sensor=false\"></scr'+'ipt>' +\n",
        "                '<scr'+'ipt>google.maps.event.addDomListener(window,\"load\",function(){ready=true});</scr'+'ipt>' +\n",
        "                '</head>' +\n",
        "                '<body style=\"height:100%; margin:0px; padding:0px\"></body></html>';\n",
        "            \n",
        "            // Create the Google Maps container element.\n",
        "            this.$iframe = $('<iframe />')\n",
        "                .css('width', '100%')\n",
        "                .css('height', '300px')\n",
        "                .attr('srcdoc', src)\n",
        "                .appendTo(this.$el);\n",
        "                        \n",
        "            var self = this; // hold onto this for initMapWhenReady\n",
        "\n",
        "            // Wait until maps library has finished loading in iframe, then create map\n",
        "            function initMapWhenReady() {\n",
        "                // Iframe document and window\n",
        "                var doc = self.$iframe[0].contentDocument;\n",
        "                var win = self.$iframe[0].contentWindow;\n",
        "                if (!win || !win.ready) {\n",
        "                    // Maps library not yet loaded;  try again soon\n",
        "                    setTimeout(initMapWhenReady, 20);\n",
        "                    return;\n",
        "                }\n",
        "\n",
        "                // Maps library finished loading.  Build map now.\n",
        "                var mapOptions = {\n",
        "                    center: new win.google.maps.LatLng(self.model.get('lat'), self.model.get('lng')),\n",
        "                    zoom: self.model.get('zoom')\n",
        "                };\n",
        "                var mapDiv = $('<div />')\n",
        "                    .css('width', '100%')\n",
        "                    .css('height', '100%')\n",
        "                    .appendTo($(doc.body));\n",
        "                self.map = new win.google.maps.Map(mapDiv[0], mapOptions);\n",
        "                \n",
        "                \n",
        "                // Add an event listeners for user panning, zooming, and resizing map\n",
        "                // TODO(rsargent): Bind self across all methods, and save some plumbing here\n",
        "                win.google.maps.event.addListener(self.map, 'bounds_changed', function () {\n",
        "                    self.handleBoundsChanged();\n",
        "                });\n",
        "                \n",
        "                self.initializeLayersControl();\n",
        "            }\n",
        "            initMapWhenReady();\n",
        "        },\n",
        "        \n",
        "        LayersControl: function(widget, controlDiv, map) {\n",
        "            var win = widget.$iframe[0].contentWindow;\n",
        "            var chicago = new win.google.maps.LatLng(41.850033, -87.6500523);\n",
        "\n",
        "            // Set CSS styles for the DIV containing the control\n",
        "            // Setting padding to 5 px will offset the control\n",
        "            // from the edge of the map.\n",
        "            controlDiv.style.padding = '5px';\n",
        "\n",
        "            // Set CSS for the control border.\n",
        "            var $controlUI = $('<div />')\n",
        "                .css('backgroundColor', 'white')\n",
        "                .css('borderStyle', 'solid')\n",
        "                .css('borderWidth', '1px')\n",
        "                .css('cursor', 'pointer')\n",
        "                .css('textAlign', 'center')\n",
        "                .appendTo($(controlDiv));\n",
        "            \n",
        "            // Set CSS for the control interior.\n",
        "            var $controlContents = $('<div />')\n",
        "                .css('fontFamily', 'Arial,sans-serif')\n",
        "                .css('fontSize', '12px')\n",
        "                .css('paddingLeft', '4px')\n",
        "                .css('paddingRight', '4px')\n",
        "                .css('paddingTop', '0px')\n",
        "                .css('paddingBottom', '0px')\n",
        "                .appendTo($controlUI);\n",
        "            \n",
        "            this.$controlTable = $('<table />')\n",
        "                .append($('<tr><td colspan=2>Layers</td></tr>'))\n",
        "                .appendTo($controlContents);\n",
        "        },\n",
        "\n",
        "        initializeLayersControl: function() {\n",
        "            var doc = this.$iframe[0].contentDocument;\n",
        "            var win = this.$iframe[0].contentWindow;\n",
        "\n",
        "            // Create the DIV to hold the control and call the LayersControl() constructor\n",
        "            // passing in this DIV.\n",
        "    \n",
        "            var layersControlDiv = document.createElement('div');\n",
        "            this.layersControl = new this.LayersControl(this, layersControlDiv, this.map);\n",
        "\n",
        "            layersControlDiv.index = 1;\n",
        "            this.map.controls[win.google.maps.ControlPosition.TOP_RIGHT].push(layersControlDiv);\n",
        "        },\n",
        "        \n",
        "        // Map geometry changed (pan, zoom, resize)\n",
        "        handleBoundsChanged: function() {\n",
        "            this.model.set('lng', this.map.getCenter().lng());\n",
        "            this.model.set('lat', this.map.getCenter().lat());\n",
        "            this.model.set('zoom', this.map.getZoom());\n",
        "            var bounds = this.map.getBounds();\n",
        "            var playgroundCompatible = [bounds.getSouthWest().lng(), bounds.getSouthWest().lat(),\n",
        "                                        bounds.getNorthEast().lng(), bounds.getNorthEast().lat()];\n",
        "            this.model.set('bounds', playgroundCompatible);\n",
        "            this.touch();\n",
        "        },\n",
        "        \n",
        "        // Receive custom messages from Python backend\n",
        "        on_msg: function(msg) {\n",
        "            var win = this.$iframe[0].contentWindow;\n",
        "            if (msg.command == 'addLayer') {\n",
        "                this.addLayer(msg.mapid, msg.token, msg.name, msg.visible);\n",
        "            } else if (msg.command == 'center') {\n",
        "                this.map.setCenter(new win.google.maps.LatLng(msg.lat, msg.lng));\n",
        "                if (msg.zoom !== null) {\n",
        "                    this.map.setZoom(msg.zoom);\n",
        "                }\n",
        "            }\n",
        "        },\n",
        "        \n",
        "        // Add an Earth Engine layer\n",
        "        addLayer: function(mapid, token, name, visible) {\n",
        "            var win = this.$iframe[0].contentWindow;\n",
        "            var eeMapOptions = {\n",
        "                getTileUrl: function(tile, zoom) {\n",
        "                    var url = ['https://earthengine.googleapis.com/map',\n",
        "                               mapid, zoom, tile.x, tile.y].join(\"/\");\n",
        "                    url += '?token=' + token\n",
        "                    return url;\n",
        "                },\n",
        "                tileSize: new win.google.maps.Size(256, 256),\n",
        "                opacity: visible ? 1.0 : 0.0,\n",
        "            };\n",
        "            \n",
        "            // Create the overlay map type\n",
        "            var mapType = new win.google.maps.ImageMapType(eeMapOptions);\n",
        "                \n",
        "            // Overlay the Earth Engine generated layer\n",
        "            this.map.overlayMapTypes.push(mapType);\n",
        "\n",
        "            // Update layer visibility control\n",
        "            var maxSlider = 100;\n",
        "            \n",
        "            function updateOpacity() {\n",
        "                mapType.setOpacity($checkbox[0].checked ? $slider[0].value / 100.0 : 0);\n",
        "            }\n",
        "            \n",
        "            var $checkbox = $('<input type=\"checkbox\">')\n",
        "                .prop('checked', visible)\n",
        "                .change(updateOpacity);\n",
        "            \n",
        "            var $slider = $('<input type=\"range\" />')\n",
        "                .prop('min', 0)\n",
        "                .prop('max', maxSlider)\n",
        "                .prop('value', maxSlider)\n",
        "                .css('width', '60px')\n",
        "                .on('input', updateOpacity);\n",
        "\n",
        "            // If user doesn't specify a layer name, create a default\n",
        "            if (name === null) {\n",
        "                name = 'Layer ' + this.map.overlayMapTypes.length;\n",
        "            }\n",
        "            \n",
        "            var $row = $('<tr />');\n",
        "            $('<td align=\"left\" />').append($checkbox).append(name).appendTo($row);\n",
        "            $('<td />').append($slider).appendTo($row);\n",
        "\n",
        "            this.layersControl.$controlTable.append($row);\n",
        "        }\n",
        "    });\n",
        "    \n",
        "    // Register the GoogleMapsView with the widget manager.\n",
        "    WidgetManager.register_widget_view('GoogleMapsView', GoogleMapsView);\n",
        "});"
       ],
       "metadata": {},
       "output_type": "display_data",
       "text": [
        "<IPython.core.display.Javascript at 0x105095050>"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Initialize Earth Engine."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# This script assumes your authentication credentials are stored as operatoring system\n",
      "# environment variables.\n",
      "MY_SERVICE_ACCOUNT = os.environ.get('MY_SERVICE_ACCOUNT')\n",
      "MY_PRIVATE_KEY_FILE = os.environ.get('MY_PRIVATE_KEY_FILE')\n",
      "\n",
      "# Initialize the Earth Engine object, using your authentication credentials.\n",
      "ee.Initialize(ee.ServiceAccountCredentials(MY_SERVICE_ACCOUNT, MY_PRIVATE_KEY_FILE))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Google Map Widget - Hello World example"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "map = GoogleMapsWidget(lat=70, lng=-140.0, zoom=6)\n",
      "display(map)\n",
      "\n",
      "image1 = ee.Image('LC81412332013146LGN00')\n",
      "vis_params1 = {'bands':'B4,B3,B2','min':5000,'max':10000}\n",
      "map.addLayer(image=image1, vis_params=vis_params1, name='Night')\n",
      "\n",
      "image2 = ee.Image('LC80660112013132LGN01')\n",
      "vis_params2 = {'bands':'B4,B3,B2','min':0,'max':40000}\n",
      "map.addLayer(image=image2, vis_params=vis_params2, name='Day', visible=False)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "(map.lat, map.lng, map.zoom, map.bounds)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 6,
       "text": [
        "(70.0, -140.0, 6, [])"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Print out the image metadata for a single image.\n",
      "print('Image Metadata') \n",
      "metadata = image1.getInfo()\n",
      "for key, value in metadata['properties'].items():\n",
      "    print('{0} = {1}'.format(key, value))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Image Metadata\n",
        "THERMAL_LINES = 8271"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "EARTH_SUN_DISTANCE = 1.013023\n",
        "SCENE_CENTER_TIME = 06:06:59.8894716Z\n",
        "CORNER_LL_PROJECTION_Y_PRODUCT = 7597800.0\n",
        "REQUEST_ID = 0501305266288_00005\n",
        "RADIANCE_MULT_BAND_8 = 0.011201\n",
        "GEOMETRIC_RMSE_MODEL = 17.422\n",
        "RADIANCE_MULT_BAND_7 = 0.00049452\n",
        "CORNER_LL_LAT_PRODUCT = 68.46586\n",
        "RADIANCE_MULT_BAND_5 = 0.0060338\n",
        "RADIANCE_MULT_BAND_4 = 0.0099434\n",
        "RADIANCE_MULT_BAND_3 = 0.01174\n",
        "RADIANCE_MULT_BAND_2 = 0.012822\n",
        "QUANTIZE_CAL_MIN_BAND_9 = 1\n",
        "RADIANCE_MAXIMUM_BAND_3 = 710.70575\n",
        "QUANTIZE_CAL_MIN_BAND_7 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_6 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_5 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_4 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_3 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_2 = 1\n",
        "QUANTIZE_CAL_MIN_BAND_1 = 1\n",
        "ELEVATION_SOURCE = GLS2000\n",
        "REFLECTANCE_MAXIMUM_BAND_8 = 1.2107\n",
        "CORNER_LR_LON_PRODUCT = -137.82565\n",
        "REFLECTANCE_MAXIMUM_BAND_6 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_7 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_4 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_5 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_2 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_3 = 1.2107\n",
        "REFLECTANCE_MAXIMUM_BAND_1 = 1.2107\n",
        "REFLECTANCE_ADD_BAND_9 = -0.1\n",
        "ELLIPSOID = WGS84\n",
        "system:time_end = 1369548419889\n",
        "catalog_cloud_cover = 0\n",
        "REFLECTANCE_ADD_BAND_8 = -0.1\n",
        "RADIANCE_ADD_BAND_1 = -62.86801\n",
        "RADIANCE_MAXIMUM_BAND_10 = 22.0018\n",
        "RADIANCE_MAXIMUM_BAND_11 = 22.0018\n",
        "BPF_NAME_OLI = LO8BPF20130526060428_20130526074321.01\n",
        "WRS_ROW = 233\n",
        "RADIANCE_ADD_BAND_3 = -58.70205\n",
        "QUANTIZE_CAL_MAX_BAND_11 = 65535\n",
        "GROUND_CONTROL_POINTS_VERIFY = 145\n",
        "RADIANCE_ADD_BAND_2 = -64.10882\n",
        "RADIANCE_ADD_BAND_5 = -30.16886\n",
        "RADIANCE_ADD_BAND_4 = -49.71705\n",
        "RADIANCE_ADD_BAND_7 = -2.47261\n",
        "CORNER_UL_LAT_PRODUCT = 70.68735\n",
        "WRS_PATH = 141\n",
        "DATA_TYPE = L1T\n",
        "FILE_NAME_BAND_QUALITY = LC81412332013146LGN00_BQA.TIF\n",
        "METADATA_FILE_NAME = LC81412332013146LGN00_MTL.txt\n",
        "QUANTIZE_CAL_MAX_BAND_4 = 65535\n",
        "QUANTIZE_CAL_MAX_BAND_7 = 65535\n",
        "QUANTIZE_CAL_MIN_BAND_10 = 1\n",
        "QUANTIZE_CAL_MAX_BAND_1 = 65535\n",
        "QUANTIZE_CAL_MAX_BAND_3 = 65535\n",
        "QUANTIZE_CAL_MAX_BAND_2 = 65535\n",
        "QUANTIZE_CAL_MAX_BAND_9 = 65535\n",
        "QUANTIZE_CAL_MAX_BAND_8 = 65535\n",
        "REFLECTIVE_LINES = 8271\n",
        "IMAGE_QUALITY_TIRS = 9\n",
        "SPACECRAFT_ID = LANDSAT_8\n",
        "STATION_ID = LGN\n",
        "CPF_NAME = L8CPF20130401_20130630.04\n",
        "CORNER_UR_LON_PRODUCT = -137.47694\n",
        "google:cloud_storage_path = gs://earthengine-public/landsat/L8/141/233/LC81412332013146LGN00.tar.bz\n",
        "REFLECTANCE_MULT_BAND_9 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_8 = 2e-05\n",
        "QUANTIZE_CAL_MIN_BAND_11 = 1\n",
        "SENSOR_ID = OLI_TIRS\n",
        "RADIANCE_MULT_BAND_6 = 0.0015202\n",
        "REFLECTANCE_MULT_BAND_1 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_3 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_2 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_5 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_4 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_7 = 2e-05\n",
        "REFLECTANCE_MULT_BAND_6 = 2e-05\n",
        "REFLECTANCE_MINIMUM_BAND_4 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_5 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_6 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_7 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_1 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_2 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_3 = -0.09998\n",
        "TARGET_WRS_PATH = 141\n",
        "GEOMETRIC_RMSE_MODEL_X = 14.286\n",
        "GEOMETRIC_RMSE_MODEL_Y = 9.972\n",
        "THERMAL_SAMPLES = 8431\n",
        "RADIANCE_MAXIMUM_BAND_8 = 678.02966\n",
        "RADIANCE_MAXIMUM_BAND_9 = 150.09816\n",
        "RADIANCE_MULT_BAND_1 = 0.012574\n",
        "RADIANCE_MAXIMUM_BAND_1 = 761.14294\n",
        "RADIANCE_MAXIMUM_BAND_2 = 776.1655\n",
        "GRID_CELL_SIZE_THERMAL = 30.0\n",
        "RADIANCE_MAXIMUM_BAND_4 = 601.92426\n",
        "RADIANCE_MAXIMUM_BAND_5 = 365.25436\n",
        "RADIANCE_MAXIMUM_BAND_6 = 92.02586\n",
        "RADIANCE_MAXIMUM_BAND_7 = 29.93584\n",
        "K2_CONSTANT_BAND_11 = 1201.14\n",
        "K2_CONSTANT_BAND_10 = 1321.08\n",
        "LANDSAT_SCENE_ID = LC81412332013146LGN00\n",
        "CORNER_LR_PROJECTION_Y_PRODUCT = 7597800.0\n",
        "SUN_AZIMUTH = -44.640205\n",
        "CORNER_LL_PROJECTION_X_PRODUCT = 377100.0\n",
        "CORNER_UR_PROJECTION_Y_PRODUCT = 7845900.0\n",
        "system:time_start = 1369548419889\n",
        "RADIANCE_MULT_BAND_9 = 0.0024795\n",
        "BPF_NAME_TIRS = LT8BPF20130526060034_20130526083119.01\n",
        "QUANTIZE_CAL_MIN_BAND_8 = 1\n",
        "CORNER_LR_PROJECTION_X_PRODUCT = 630000.0\n",
        "GEOMETRIC_RMSE_VERIFY = 229.931\n",
        "GRID_CELL_SIZE_REFLECTIVE = 30.0\n",
        "ORIENTATION = NORTH_UP\n",
        "REFLECTANCE_MINIMUM_BAND_8 = -0.09998\n",
        "REFLECTANCE_MINIMUM_BAND_9 = -0.09998\n",
        "RADIANCE_MULT_BAND_11 = 0.0003342\n",
        "RADIANCE_MULT_BAND_10 = 0.0003342\n",
        "RADIANCE_MINIMUM_BAND_10 = 0.10033\n",
        "RADIANCE_MINIMUM_BAND_11 = 0.10033\n",
        "PANCHROMATIC_SAMPLES = 16861\n",
        "ORIGIN = Image courtesy of the U.S. Geological Survey\n",
        "TARGET_WRS_ROW = 233\n",
        "UTM_ZONE = 7\n",
        "PROCESSING_SOFTWARE_VERSION = \"LPGS_2.2.2\"\n",
        "REFLECTANCE_MAXIMUM_BAND_9 = 1.2107\n",
        "CORNER_UL_PROJECTION_Y_PRODUCT = 7845900.0\n",
        "SUN_ELEVATION = 6.951576\n",
        "MAP_PROJECTION = UTM\n",
        "IMAGE_QUALITY_OLI = 9\n",
        "QUANTIZE_CAL_MAX_BAND_5 = 65535\n",
        "FILE_DATE = 1369597247000\n",
        "RESAMPLING_OPTION = CUBIC_CONVOLUTION\n",
        "QUANTIZE_CAL_MAX_BAND_6 = 65535\n",
        "CORNER_UR_LAT_PRODUCT = 70.68376\n",
        "DATE_ACQUIRED = 2013-05-26\n",
        "PANCHROMATIC_LINES = 16541\n",
        "NADIR_OFFNADIR = NADIR\n",
        "CORNER_LR_LAT_PRODUCT = 68.46267\n",
        "REFLECTIVE_SAMPLES = 8431\n",
        "GRID_CELL_SIZE_PANCHROMATIC = 15.0\n",
        "QUANTIZE_CAL_MAX_BAND_10 = 65535\n",
        "system:footprint = {u'type': u'LinearRing', u'coordinates': [[-139.6057110296895, 70.71121962280009], [-144.20239509257885, 69.95059426786496], [-142.06816219962155, 68.49335821912977], [-137.7280385272379, 69.2112752477388], [-139.6057110296895, 70.71121962280009]]}\n",
        "GROUND_CONTROL_POINTS_MODEL = 45\n",
        "FILE_NAME_BAND_3 = LC81412332013146LGN00_B3.TIF\n",
        "FILE_NAME_BAND_2 = LC81412332013146LGN00_B2.TIF\n",
        "FILE_NAME_BAND_1 = LC81412332013146LGN00_B1.TIF\n",
        "RADIANCE_MINIMUM_BAND_5 = -30.16282\n",
        "FILE_NAME_BAND_7 = LC81412332013146LGN00_B7.TIF\n",
        "FILE_NAME_BAND_6 = LC81412332013146LGN00_B6.TIF\n",
        "FILE_NAME_BAND_5 = LC81412332013146LGN00_B5.TIF\n",
        "FILE_NAME_BAND_4 = LC81412332013146LGN00_B4.TIF\n",
        "CLOUD_COVER = 8.75\n",
        "FILE_NAME_BAND_9 = LC81412332013146LGN00_B9.TIF\n",
        "FILE_NAME_BAND_8 = LC81412332013146LGN00_B8.TIF\n",
        "RADIANCE_MINIMUM_BAND_8 = -55.99191\n",
        "RADIANCE_MINIMUM_BAND_9 = -12.39515\n",
        "RADIANCE_ADD_BAND_9 = -12.39763\n",
        "RADIANCE_ADD_BAND_8 = -56.00311\n",
        "CORNER_UL_PROJECTION_X_PRODUCT = 377100.0\n",
        "FILE_NAME_BAND_11 = LC81412332013146LGN00_B11.TIF\n",
        "FILE_NAME_BAND_10 = LC81412332013146LGN00_B10.TIF\n",
        "REFLECTANCE_ADD_BAND_7 = -0.1\n",
        "REFLECTANCE_ADD_BAND_6 = -0.1\n",
        "REFLECTANCE_ADD_BAND_5 = -0.1\n",
        "REFLECTANCE_ADD_BAND_4 = -0.1\n",
        "REFLECTANCE_ADD_BAND_3 = -0.1\n",
        "REFLECTANCE_ADD_BAND_2 = -0.1\n",
        "REFLECTANCE_ADD_BAND_1 = -0.1\n",
        "RADIANCE_ADD_BAND_6 = -7.60105\n",
        "RADIANCE_MINIMUM_BAND_6 = -7.59953\n",
        "RADIANCE_MINIMUM_BAND_7 = -2.47211\n",
        "OUTPUT_FORMAT = GEOTIFF\n",
        "RADIANCE_MINIMUM_BAND_4 = -49.7071\n",
        "K1_CONSTANT_BAND_10 = 774.89\n",
        "K1_CONSTANT_BAND_11 = 480.89\n",
        "RADIANCE_MINIMUM_BAND_2 = -64.096\n",
        "CORNER_UL_LON_PRODUCT = -144.33105\n",
        "RADIANCE_MINIMUM_BAND_3 = -58.69032\n",
        "RADIANCE_ADD_BAND_11 = 0.1\n",
        "RADIANCE_ADD_BAND_10 = 0.1\n",
        "CORNER_UR_PROJECTION_X_PRODUCT = 630000.0\n",
        "RADIANCE_MINIMUM_BAND_1 = -62.85543\n",
        "DATUM = WGS84\n",
        "catalog_image_quality1 = 9\n",
        "RLUT_FILE_NAME = L8RLUT20130211_20431231v06.h5\n",
        "ROLL_ANGLE = -0.001\n",
        "CORNER_LL_LON_PRODUCT = -144.00128\n"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Image Collection Example"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "lng, lat = (-120.6, 40.0)\n",
      "map2 = GoogleMapsWidget(lat=lat, lng=lng, zoom=6)\n",
      "display(map2)\n",
      "\n",
      "# Add an Image Collection of Landsat 8 data to the map.\n",
      "collection = (ee.ImageCollection('LC8')\n",
      "                .filterDate('2014-01-01', '2014-06-01')\n",
      "                .filterMetadata('SUN_ELEVATION', 'greater_than', 0)\n",
      "             )\n",
      "vis_params3 = {'bands': 'B4,B3,B2', 'min':5000, 'max':30000, 'gamma': 1.6}\n",
      "map2.addLayer(image=collection.mosaic(), vis_params=vis_params3, name='Landsat 8 mosaic')\n",
      "map2.addLayer(image=collection.median(), vis_params=vis_params3, name='Landsat 8 median', visible=False)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Extract information from the image collection"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Get the coordinates of the center of the map.\n",
      "(map2.lng, map2.lat)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 9,
       "text": [
        "(-120.6, 40.0)"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "point = ee.Geometry.Point(map2.lng, map2.lat);\n",
      "info = collection.getRegion(point,500).getInfo()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stderr",
       "text": [
        "WARNING:oauth2client.util:new_request() takes at most 1 positional argument (4 given)\n"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "xBand = 'time'\n",
      "yBandList = ['B4','B5']\n",
      "# extract the header column names\n",
      "header = info[0]\n",
      "# create a Numpy array of the data\n",
      "data = np.array(info[1:])\n",
      "# extract the time information\n",
      "iTime = header.index(xBand)\n",
      "# convert to Python datetime objects\n",
      "time = [datetime.datetime.fromtimestamp(i/1000) for i in (data[0:,iTime].astype(int))]\n",
      "# Extract the data columns we want to plot.\n",
      "iBands = [header.index(b) for b in yBandList]\n",
      "yData = data[0:,iBands].astype(np.float)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 11
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "red = yData[:,0]\n",
      "nir = yData[:,1]\n",
      "ndvi = (nir - red) / (nir + red)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 12
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "fig, ax = plt.subplots(2, sharex=True, figsize=(10,10))\n",
      "ax[0].plot(time, red, 'o', color=\"red\", label=\"RED band\")\n",
      "ax[0].plot(time, nir, 'o', color=\"magenta\",  label=\"NIR band\")\n",
      "ax[0].set_ylabel('Band Values')\n",
      "ax[0].legend(loc='best')\n",
      "ax[1].plot(time, ndvi, 'o', color=\"black\", label=\"NDVI\")\n",
      "ax[1].set_ylabel('NDVI')\n",
      "ax[1].legend(loc='best')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 13,
       "text": [
        "<matplotlib.legend.Legend at 0x108872510>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAm4AAAJPCAYAAADFWFiAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xt8XHWd+P9XmkBC5dLWQgulXHZEBX8CggL7c62D2kw0\nQkF9QHHBcvnuz4fIJLBekKaFlLYrKKjN+JCfLtd2FcVFwd2BNEWNcX+7wupCQbkIkQotpWDpFhAT\nScjvj3OSTtJJOm3mdmZez8djHnPmPeec+WTG3b55f24gSZIkSZIkSZIkSZIkSZIkSZIkSZIkSZIk\nqYrUAg8C/xa+ngGsA34PdAHTMs69AngSeBxozIifCDwSvrcqI14P/CCM/wo4PP/NlyRJqh7/CHwX\n+En4+ivAF8Pjy4FrwuNjgIeAvYAjgKeAmvC9B4CTwuN7gKbw+GLgW+Hx2cD38956SZKkKnEocB9w\nKjsqbo8Ds8Lj2eFrCKptl2dc2wmcAhwMPJYRXwj8vxnnnBwe1wEv5rHtkiRJZWVKge//deALwBsZ\nsVnAlvB4CzuSuEOAjRnnbQTmZIlvCuOEz8+GxwPAdoKuWEmSpIpTyMTto8ALBOPbasY5Zyh8SJIk\naRfqCnjv/xs4HfgI0ADsD6whqLLNBp4n6AZ9ITx/EzA34/pDCSptm8LjsfHhaw4DniP4Ww4AXhrb\nkOOOO25o/fr1+fibJEmSCm09cHy2NwpZcVtMkIgdSTAu7WfAeQSTFBaF5ywC7gqPfxKet3d4zVEE\nkxKeB14mGMtWE97j7oxrhu/1CeCn2Rqyfv16hoaGuOqqqxgaGvIxiYffod9juTz8Dv0Oy+Xh9+h3\nmO8HcNx4yVUhK25jDXeJXgPcAVwEbADOCuOPhvFHCcarXZxxzcXArcA+BLNKO8P4TQRVvCeBrQSJ\nnyRJUkUqVuL2i/ABQVfmh8Y575/Cx1i/Ad6ZJd7PjsRPkiSpohV6VmlZicfjpW5C5Pkd5off4+T5\nHU6e32F++D1Ont9h7sab7VlphsI+Y0mSpLJWU1MD4+RoxRzjVnZmzJjBtm3bSt0MZZg+fTovvbTT\nxGBJkkSVV9xqamqwElde/E0kSdVuoopbVY1xkyRJijITN0mSpIgwcZMkSYoIEzeNmDJlCn/4wx8q\n7rMkSaoUJm5l6ogjjmDq1Knst99+zJ49m/POO4+XX3555P3zzz+f+vp69ttvv5HHu971LgA2bNjA\nlClTRuKzZ8/mtNNO47777ivVnyNJkvLAxC2LnnSaJYkE7fE4SxIJetLpot+jpqaGf//3f+eVV15h\n/fr1PPLII6xYsWLU+5dffjmvvPLKyOPBBx8cdY/t27fzyiuv8PDDDzN//nzOPPNMbrvttt3+WyRJ\nUnkwcRujJ51mbWsrK7q6aP/FL1jR1cXa1tbdSrzycY9Ms2bNorGxkd/97nd7dP1BBx1ES0sL7e3t\nXH755ROem06nicViHHjggXzxi18cWZqjt7eXD3zgA8ycOZMDDzyQc889l+3bt49cd8QRR3D99ddz\n3HHHMW3aNBYuXEh/f//I+1/96lc55JBDOPTQQ7n55pv36O+QJKnambiN0dXRwcre3lGxlb29rEul\ninoPYCRp2rhxI52dnZx88slZ38/VmWeeyQsvvMATTzwx7jl33XUXv/nNb/if//kf7r777lFJVltb\nG5s3b+axxx7j2Wefpb29feS9mpoafvjDH7J27VqefvppHn74YW699VYAOjs7uf7667nvvvv4/e9/\nb5etJEl7yMRtjLqMKlGm2r6+ot5jaGiIM844g/3335/DDjuMWCzGkiVLRr1/3XXXMX369JHHBRdc\nMOE9DznkEIAJdya4/PLLmTZtGnPnzuXSSy/l9ttvByAWi/HBD36Qvfbai5kzZ3LZZZfxi1/8YtS1\nLS0tzJ49m+nTp3Paaafx0EMPAXDHHXdw4YUXcswxxzB16lSWLVuW8/cgSZJ2MHEbY6C+Pmt8sKGh\nqPeoqanh7rvv5uWXX6a7u5uf/exn/PrXvx71/he+8AW2bds28rjlllsmvOemTZuAYKuv8cydO3fk\n+LDDDuO5554DYMuWLSxcuJBDDz2UAw44gPPOO4+tW7eOunb27Nkjx/vssw9//vOfAdi8efNO95Uk\nSbvPxG2MxpYW2mKxUbHFsRjzk8mi3iPTvHnzSCaTuxyftis//vGPmTVrFm9729vGPeeZZ54ZdTxn\nzhwAFi9eTG1tLb/97W/Zvn07a9as4Y033sjpcw8++OCd7itJknZfVW8yn8285mYAlqZS1Pb1MdjQ\nQFMyORIv1j3GuvTSS/n617/O/fffz8knn8zQ0NAux7gNv79lyxZ++MMfcvXVV9PR0THhNddddx0n\nn3wyr7zyCh0dHXzuc58D4NVXX+WAAw5g//33Z9OmTXz1q1/dZZuHP/+ss87iggsu4FOf+hSHH364\nXaWSJO0hE7cs5jU3TyrJytc9Ms2cOZNFixZx7bXX8qMf/Yiamhq+8pWv8I1vfGPknH322YcXXnhh\n5PW0adMYGhriTW96E+95z3v413/9VxobGyf8nAULFnDiiSeyfft2LrjgAi688EIArrrqKj71qU9x\nwAEHcNRRR3HuueeO+uyxampqhjfJpampiUsvvZQPfOAD1NbWsnz58pGxc5IkKXdZd56vQEPZqlM1\nNTW7PTNTheVvIkmqdmHhI2uO5hg3SZKkiDBxkyRJiggTN0mSpIgwcZMkSYoIZ5VKkqSq1ZPuoauj\ni7r+OgbqB2hsaWRe87xSN2tcJm6SJKkq9aR7WNu6lpW9K0dibb1tAGWbvNlVKkmSqlJXR9eopA1g\nZe9K1qXWlahFu2biJkmSqlJdf/aOx9q+2iK3JHcmblXq1ltv5X3ve1/FfZYkSbkaqB/IGh9sGCxy\nS3JXyMStAbgfeAh4FPhyGG8HNgIPho8PZ1xzBfAk8DiQuTfTicAj4XurMuL1wA/C+K+Aw/P8N5TM\nEUccwaxZs3jttddGYjfeeCOnnnrqyOspU6bwhz/8AYD29nb22msv9ttvP6ZNm8Ypp5zCL3/5y6K3\nW5KkqGhsaaQt1jYqtji2mPnJ+SVq0a4VcnJCH3Aq8Fr4Of8B/B0wBHwtfGQ6Bjg7fJ4D3AccFZ5/\nA3AR8ABwD9AEdIaxreF5ZwPXAgsn2/B8zDDJxz3eeOMNVq1axRVXXLHLc2tqajjnnHNYvXo1g4OD\nXHXVVXziE59gy5Ytu/WZkiRVi+F/l5emllLbV8tgwyBNyaaynZgAhe8qHS4X7Q3UAtvC19n231oA\n3A68DmwAngJOBg4G9iNI2gBWA2eEx6cDt4XHdwIfnGyDh2eYrOhaQfsv2lnRtYK1rWvpSfcU9R41\nNTV8/vOf57rrrmP79u27PH9oaGhkj8/a2lo++clP8uKLL/KnP/1pwmuSySTTpk3j6KOP5mc/+9nI\ne7fccgvHHHMM+++/P7FYjO985zsj73V3d3PooYfyta99jVmzZnHIIYdw6623jry/detWTj/9dA44\n4ABOPvlkent7c/67JUkqpnnN81jeuZz27naWdy4v66QNCp+4TSHoKt0C/Bz4XRhPAuuBm4BpYewQ\ngi7UYRsJKm9j45vCOOHzs+HxALAdmDGZBudjhkm+Zqm8+93vJh6Pc9111+3WdX/9619ZvXo1sViM\nmTNnjnve/fffz1ve8ha2bt3KsmXL+NjHPsa2bUFuPWvWLNLpNC+//DK33HILl112GQ8++ODItVu2\nbOHll1/mueee46abbuKzn/3sSIL52c9+lqlTp/L8889z8803c8sttwxvmCtJkiah0InbG8DxwKHA\nPCBO0O15ZBjfDFxf4DbslnzMMMnXLJWamhquvvpqUqnUhJWzYXfccQfTp09n6tSp3Hjjjdxzzz0T\nnn/QQQfR2tpKbW0tZ511Fm9729tIp9MAfOQjH+HII48EYN68eTQ2No4aM7fXXntx5ZVXUltby4c/\n/GH23XdfnnjiCQYHB/nRj37E1VdfzT777MM73vEOFi1aNFINlCRJe65YC/BuB9LAu4HujPiNwL+F\nx5uAuRnvHUpQadsUHo+ND19zGPAcwd9yAPBStga0t7ePHMfjceLxeNaG5mOGST5nqbzjHe/gox/9\nKNdccw1HH330hOeeffbZrF69mq1bt/Lxj3+cVCpFR0fHuOfPmTNn1OvDDz+czZs3A3DvvfeybNky\nnnzySd544w1ee+01jj322JFz3/zmNzNlyo68f+rUqbz66qu8+OKLDAwMMHfujp/ysMMO262/WZKk\natLd3U13d3dO5xYycZtJ0H35v8A+wHxgGTAbeD4850yC2aIAPwG+RzBpYQ7BhIMHCCYnvEww3u0B\n4DygI+OaRQQzSj8B/HS8xmQmbhNpbGmkrbdtVFfn4thimpJNOV2fr3tkWrZsGSeccAKf+9znJjxv\nuKr15je/me985zsce+yxXHrppfzN3/xN1vM3bdo06vUf//hHFixYQH9/Px//+Mf5l3/5FxYsWEBt\nbS1nnnlmTlWzAw88kLq6Op555hne9ra3AfDMM8/k8mdKklSVxhaUli1bNu65hUzcDiaYODAlfKwh\nSKxWE3STDgFPA58Oz38UuCN8HgAuDs8hPL6VIAG8h2BGKQRj5NYQLAeylTzMKM3HDJN8z1KJxWKc\nffbZrFq1alTVayJvfetbOe200/jqV7/KDTfckPWcF154gY6ODj7zmc9w11138fjjj/ORj3yEv/71\nr/z1r39l5syZTJkyhXvvvZeuri7e+c537vJza2tr+djHPkZ7ezs333wzTz/9NLfddtu4yaMkScpd\nIRO3R4ATssQ/NcE1/xQ+xvoNkC1r6AfO2v2mTWxe87xJzyrJxz0yXXnllaxZs2bUIP+xx2MnAHzh\nC1/g1FNPZdmyZRx00EGj3qupqeGUU07hySef5MADD2T27NnceeedTJ8+HYCOjg7OOuss+vv7Oe20\n01iwYMFO14/nm9/8JhdccAGzZ8/m6KOP5sILL8y5BCxJksZXLVP9hrJ189XU1Dhovsz4m0iSql1Y\nHMmao7nllSRJUkSYuEmSJEWEiZskSVJEmLhJkiRFhImbJElSRJi4SZIkRUSxtrwqS9OnT3fz8zIz\nvI6cJEnaWbVkLVnXcZMkSSo3ruMmSZJUAUzcJEmSIsLETZIkKSJM3CRJkiLCxE2SJCkiTNwkSZIi\nwsRNkiQpIkzcJEmSIsLETZIkKSJM3CRJkiKiqvcqlSRJ1a0nnaaro4O6/n4G6utpbGlhXnNzqZs1\nLhM3SZJUlXrSada2trKyt3ck1hYel2vyZlepJEmqSl0dHaOSNoCVvb2sS6VK1KJdM3GTJElVqa6/\nP2u8tq+vyC3JnYmbJEmqSgP19Vnjgw0NRW5J7kzcJElSVWpsaaEtFhsVWxyLMT+ZLFGLdq2m1A0o\nkqGhoaFSt0GSJJWZnnSadakUtX19DDY0MD+ZLPnEhJqaGhgnRzNxkyRJKiMTJW52lUqSJEVEIRO3\nBuB+4CHgUeDLYXwGsA74PdAFTMu45grgSeBxoDEjfiLwSPjeqox4PfCDMP4r4PB8/xGSJEnlopCJ\nWx9wKnA8cGx4/HfAlwgSt7cCPw1fAxwDnB0+NwHfYkeZ8AbgIuCo8NEUxi8CtoaxrwPXFvDvkSRJ\nKqlCd5W+Fj7vDdQC24DTgdvC+G3AGeHxAuB24HVgA/AUcDJwMLAf8EB43uqMazLvdSfwwQL8DZIk\nSWWh0InbFIKu0i3Az4HfAbPC14TPs8LjQ4CNGdduBOZkiW8K44TPz4bHA8B2gq5YSZKkilPovUrf\nIOgqPQBYS9BdmmkofBRce3v7yHE8HicejxfjYyVJkibU3d1Nd3d3TucWczmQpcBfgP8DxIHnCbpB\nfw68nR1j3a4JnzuBq4A/huccHcbPAeYBnwnPaSeYmFAHbAYOzPLZLgciSZIioVTLgcxkx4zRfYD5\nwIPAT4BFYXwRcFd4/BNgIcF4uCMJJhw8QJDgvUww3q0GOA+4O+Oa4Xt9gmCygyRJUkUqZFfpwQQT\nB6aEjzUEidWDwB0EM0I3AGeF5z8axh8lGK92MTu6US8GbiVIAO8hqLQB3BTe90mC2aULC/fnSJIk\nlZY7J0iSJJURd06QJEmqACZukiRJEWHiJkmSFBEmbpIkSRFh4iZJkhQRJm6SJEkRYeImSZIUEYXe\nq1RSBelJ99DV0UVdfx0D9QM0tjQyr3leqZslSVXDxE1STnrSPaxtXcvK3pUjsbbeNgCTN0kqErtK\nJeWkq6NrVNIGsLJ3JetS60rUIkmqPiZuknJS15+9QF/bV1vklkhS9TJxk5STgfqBrPHBhsEit0SS\nqpeJm6ScNLY00hZrGxVbHFvM/OT8ErVIkqpP1p3nK9DQ0NBQqdsgRV5Puod1qXXU9tUy2DDI/OR8\nJyZIUp7V1NTAODmaiZskSVIZmShxs6tUkiQpIkzcJEmSIsLETZIkKSLcOUFSznrSabo6Oqjr72eg\nvp7GlhbmNTeXullVxW3HpOpm4iYpJz3pNGtbW1nZ2zsSawuPTd6Kw23HJO1uV+kM4NhCNERSeevq\n6BiVtAGs7O1lXSpVohZVH7cdk5RL4vYLYH+CpO03wI3A1wvZKEnlp66/P2u8tq+vyC2pXm47JimX\nxO0A4GXgY8Bq4CTgQ4VslKTyM1BfnzU+2NBQ5JZUL7cdk5RL4lYLHAycBaTDmKvZSlWmsaWFtlhs\nVGxxLMb8ZLJELao+bjsmKZfJCVcDa4H/D3gAiAFPFrJRksrP8ASEpakUtX19DDY00JRMOjGhiIYn\nICxNLR3Zdqwp2eTEBKmKuOWVJElSGZnslldvA34K/C58fSywJC8tkyRJUs5ySdz+GVgM/DV8/Qhw\nTo73nwv8nCDp+y3QEsbbgY3Ag+HjwxnXXEHQFfs40JgRPzH87CeBVRnxeuAHYfxXwOE5tk2SJClS\nckncpgL3Z7weAl7P8f6vA5cB7wBOAT4LHB3e42vAu8LHveH5xwBnh89NwLfYUSq8AbgIOCp8NIXx\ni4CtYezrwLU5tk2SJClSckncXgTekvH6E8DmHO//PPBQePwq8BgwJ3ydre92AXA7QcK3AXgKOJlg\nVut+BJMjIFiW5Izw+HTgtvD4TuCDObZNkiQpUnJJ3C4Bvg28HXiOoIL2mT34rCMIqmu/Cl8ngfXA\nTcC0MHYIQRfqsI0Eid7Y+CZ2JIBzgGfD4wFgO8FiwZIkSRUll8Stl6CKNZNgosJ7Caphu2Nf4F+B\nVoLK2w3AkcDxBNW763fzfpIkSVUnl3XcriIYk1bD6IV3r87xM/Yi6ML8F+CuMPZCxvs3Av8WHm8i\nmNAw7FCCStum8HhsfPiawwiqgXUEOz28NLYR7e3tI8fxeJx4PJ5j8yVJkgqnu7ub7u7unM7NZR23\nz7MjYdsH+CjwKHBhjve/jWDywGUZ8YPZMU7uMuA9wCcJJiV8j2BbrTnAfQTj64YIJki0EIxzSwMd\nQCdwMfBOgu7bhQRj3xaOaYfruBVIT7qHro4u6vrrGKgfoLGl0cVAJUmahInWccul4nbdmNdfBbpy\n/Oz3AucCDxMs+wHB0iLnEHSTDgFPA58O33sUuCN8HiBIyoYzrouBWwmSx3sIkjYIxsitIVgOZCs7\nJ20qkJ50D2tb17Kyd+VIrK032I7H5E2SpPzbk50TZhBUvd6yqxPLiBW3AliSWMKKrhU7xZcmlrK8\nc3kJWiRJUvRNtuL2SMbxFOAgch/fpgpW15/9fz61fbVFbokkSdUhl8TttIzjAWALuS/Aqwo2UD+Q\nNT7YMFjklkiSVB0mWg5kRvh4OePxGsFCuK6TJhpbGmmLtY2KLY4tZn5yfolaJFW+nnSaJYkE7fE4\nSxIJetLpUjdJUhFNVHH7H0Yv/zHWkXluiyJmeALC0tRSavtqGWwYpCnZ5MQEqUB60mnWtraysrd3\nJNYWHs9rbi5VsyQV0Z5MTogiJydIirwliQQrunae1L80kWB5Z2eWKyRF0WQnJwBMJ9jEvSEj1jO5\nZqkS9KTTdHV0UNffz0B9PY0tLf6Xv1Qgdf39WeO1fX1FbomkUsklcfsHgoVv5xKsxXYK8F/ABwrY\nLkWA3TZScQ3U12eNDzY0ZI1Lqjy57FXaSrCTwQbgVIKN4rcXsE2KiK6OjlFJG8DK3l7WpVIlapFU\n2RpbWmiLxUbFFsdizE8mS9QiScWWS8WtD/hLeNwAPE6w2byqnN02UnENV7KXplLU9vUx2NBAUzJp\nhVuqIrkkbhsJxrjdBawDthFU31Tl7LaRim9ec7OJmlTFJuoq/SLBuLYzCJK1dmApcGMYU5Wz20aS\npOKaqOJ2CPCfBNW124EfAt2Fb5Kiwm4bSZKKa1fruE0B5gELgQXAw8D3gB8BrxS2aXnlOm6SJCkS\nJlrHbXcW4K0FPgRcQzA5YeqkW1Y8Jm6SJCkS8rEA77EEVbezgD8BV+SlZZIkScrZRInbWwmStbOB\nNwjGuTUCfyhCuyRJkjTGRF2lvcD3CRK23xanOQVjV6kkSYqEfI1xizITN0mSFAkTJW65bHklSZKk\nMmDiJkmSFBEmbpIkSREx0azSRyZ4b4hgiRCpbPSke+jq6KKuv46B+gEaWxqZ1zyv1M2SJClvJkrc\nTgufLw6f1xAMlPv7grZI2gM96R7Wtq5lZe/KkVhbbxuAyZskqWLkMqv0IeD4MbEHgXflvzkF46zS\nCrcksYQVXSt2ii9NLGV55/IStEiSpD0z2VmlNcDfZbx+73g3k0qlrj978bi2r7bILZEkqXBy2fLq\nQuAW4IDw9f8CFxSsRdIeGKgfyBofbBgsckskSSqcXCpuvyGYiDD8OA74n0I2StpdjS2NtMXaRsUW\nxxYzPzm/RC2SJCn/cunybAA+DhzBjgrdEHB1DtfOBVYDB4XXfAfoAGYAPwAOBzYQbF7/v+E1VxBU\n+QaBFqArjJ8I3Bq25x6gNYzXh59xArCVYG/VP45ph2PcqkBPuod1qXXU9tUy2DDI/OR8JyZIkiJn\nslterSVIqn5DkEwNuz6Ha2eHj4eAfcN7nEHQ1fon4CvA5cB04EvAMcD3gPcAc4D7gKMIkr4HgEvC\n53sIEsBOglmv/1f4fDZwJrBwTDtM3CRJUiRMlLjlMsZtDpDYw89+PnwAvAo8Ft7vdOD9Yfw2oJsg\ncVtAsKn96wSVuKeAkwkqaPsRJG0QVNjOIEjcTgeuCuN3At/cw7ZKkiSVtVzGuP0n+Vls9wiCJUTu\nB2YBW8L4lvA1wCHAxoxrNhIkemPjm8I44fOz4fEAsJ2gK1aSJKmi5FJxex9B1+bTQH8Y292dE/Yl\nqIa1Aq+MeW8ofEiSJGkCuSRuH57kZ+xFkLStAe4KY1sIxr49DxwMvBDGNxFMaBh2KEGlbVN4PDY+\nfM1hwHMEf88BwEtjG9He3j5yHI/Hicfje/wHSZIk5Ut3dzfd3d05nbs7C+keRDCjc9gzOd7/NoLZ\nnpdlxL8Sxq4lGNs2jdGTE05ix+SEtxBU5O4nmGX6AJBm9OSEdwKfIZiUcAZOTpAkSRE12VmlpxPM\nID2EoDJ2OMEkg3fkcO3fAT3Aw+zoDr2CIPm6g6BStoHRy4EsJlgOZICga3VtGB9eDmQfglmlLWG8\nnqCa9y6CZHBheM9MJm6SJCkSJpu4PQx8AFhHkBydCpxHkFxFhYmbJEmKhMkuB/I6wZprU4Ba4OfA\nqnw1TsqXnnSaro4O6vr7Gaivp7GlhXnNzaVuliRJeZNL4raNYA21XwLfJegufbWQjZJ2V086zdrW\nVlb29o7E2sJjkzdJUqXIpat0X+AvBBW3vwf2J0jgthawXflmV2mFW5JIsKKra6f40kSC5Z2dJWiR\nJEl7ZrJdpcPVtUGC2ZxbgTfy0jIpT+r6+7PGa/v6itwSSZIKZ6KdE/6WYCuqHxFs4P5b4BGCtdcm\nu7ablFcD9fVZ44MNDVnjkiRF0USJ2zeBfyLYO/RnwP8hWDR3HvDlwjdNyl1jSwttsdio2OJYjPnJ\nZIlaJElS/k00xu0h4Pjw+DHg6Iz3HiRYGiQqHONWBXrSadalUtT29THY0MD8ZNKJCZKkyNnTddwy\nk7OxiZqJmyRJUgHsaeI2CLwWHu9DMLOUjNe5TGwoFyZukrSbetI9dHV0Uddfx0D9AI0tjcxrnlfq\nZkkVb09nldYWpDWSpLLXk+5hbetaVvauHIm19bYBmLxJJTTR5ARJUpXq6ugalbQBrOxdybrUuhK1\nSBKYuEmSsqjrz94hU9tnZ4xUSiZukqSdDNQPZI0PNgwWuSWSMpm4SZJ20tjSSFusbVRscWwx85Pz\nS9QiSZDbXqWVwFmlkrSbetI9rEuto7avlsGGQeYn5zsxQSqCPV0OpJKYuEmSpEiYKHGzq1SSJCki\norSIriRJKhAXXI4GEzdJkqqcCy5Hh12lkiRVORdcLg896R6WJJZMeI6JmyRJVc4Fl0tvuOq5omvF\nhOeZuEmSVOVccLn0slU9szFxkySpyrngcumNV/Xc6bwCt0OSJJW54QkIS1NLRxZcbko2OTGhiMar\neo7lArySJEklljmzt4bxF+C14iZJklRimVVP1o5/nhU3SZKkMlLKLa9uBrYAj2TE2oGNwIPh48MZ\n710BPAk8DjRmxE8M7/EksCojXg/8IIz/Cjg8r62XJEkqI4VO3G4BmsbEhoCvAe8KH/eG8WOAs8Pn\nJuBb7Mg2bwAuAo4KH8P3vAjYGsa+DlxbiD9CkiSpHBQ6cfslsC1LPFv5bwFwO/A6sAF4CjgZOBjY\nD3ggPG81cEZ4fDpwW3h8J/DBfDRakiSpHJVqHbcksB64CZgWxg4h6EIdthGYkyW+KYwTPj8bHg8A\n24EZhWmyJElSaZUicbsBOBI4HtgMXF+CNkiSJEVOKZYDeSHj+Ebg38LjTcDcjPcOJai0bQqPx8aH\nrzkMeI7gbzkAeCnbh7a3t48cx+Nx4vH4HjZfkqTi6En30NXRRV1/HQP1AzS2NLoobgXq7u6mu7s7\np3OLsRzIEQTJ2TvD1wcTVNoALgPeA3ySYFLC94CTCLpA7wPeQjCZ4X6ghWCcWxroADqBi8P7fgZY\nSDD2bWFHpIq2AAAgAElEQVSWNrgciCQpUjIXZB3WFmsjsSph8lbhJloOpNCJ2+3A+4GZBMuCXAXE\nCbpJh4CngU+H7wEsBi4kGK/Wyo4l6E4EbgX2Ae4hSOIgWA5kDcHs1K0ESduGLO0wcZMkRcqSxBJW\ndK3YKb40sZTlnctL0KLiqfZK40SJW6G7Ss/JErt5gvP/KXyM9Rt2VOwy9QNn7UG7JEm70JNO09XR\nQV1/PwP19TS2tDCvubnUzaoa4206XttXW+SWFFfWSmNvG0BVJW/jccsrSdJOetJp1ra2srK3dyTW\nFh6bvBXHeJuODzYMFrklxdXV0TUqaQNY2buSpamlJm6UbjkQSVIZ6+roGJW0Aazs7WVdKlWiFlWf\nxpZG2mJto2KLY4uZn5xfohYVR7VWGnNlxU2StJO6/v6s8dq+viK3JFCNY54yNx2v7atlsGGQpmRT\nxf/d1VppzJWJmyRpJwP19Vnjgw0NRW5JdY95mtc8r+L/xrEaWxpp620b9Xsvji2mKTl2B83qZOIm\nSdpJY0sLbb29o7pLF8diNCWTRW+LY56qS7VWGnNl4iZJ2snwBISlqRS1fX0MNjTQlEyWZGKCY56q\nTzVWGnNl4iZJympec3NZzCB1zJO0g7NKJUllrVpnV0rZFGPLq3LgzgmSFGE96R7WpdaNjHman5xv\nV5oqVim3vCoXJm6SJCkSJkrc7CqVJEmKCBM3SZKkiDBxkyRJiggTN0mSpIhwHTdJkjJU476oAD3p\nNF0dHdT19zNQX09jS0tZrOOn0UzcJEkKVeu+qD3pNGtbW0dtcdYWHpu8lRe7SiVJCo23L+q61LoS\ntag4ujo6RiVtACt7e1mXSpWoRRqPiZskSaFq3Re1rr8/a7y2r6/ILdGu2FUqKZKqdRySCqta90Ud\nqK/PGh9saChyS7QrJm6SIqdaxyGp8BpbGmnrbRv1v63FscU0JZtK2KrCa2xpoa23d1R36eJYjKZk\nsoStUjZueSUpcpYklrCia8VO8aWJpSzvXF6CFqmSVOu+qD3pNOtSKWr7+hhsaGB+MunEhBKZaMsr\nK26SIqdaxyGpOOY1z6uKRG2sec3NJmoRUDWTE5YkltCT7il1MyTlQbWOQ5KkqkncVnStYG3rWpM3\nqQI0tjTSFmsbFVscW8z85PwStUiSiqN6xrgRjHFzDIxUGap1HJKkyucYtwyOgZEqQ7WOQ5JU3aou\ncXMMjPLBNcSqi7+3pHJRVYlbNazFo8JzDbHq4u8tqZwUenLCzcAW4JGM2AxgHfB7oAuYlvHeFcCT\nwONAY0b8xPAeTwKrMuL1wA/C+K+Aw8dryNLEUppWNfn/aDVp1bqXYbXy95ZUTgqduN0CjC1xfYkg\ncXsr8NPwNcAxwNnhcxPwLXYMzLsBuAg4KnwM3/MiYGsY+zpw7XgNWd653KRNeeEaYtXF31tSOSl0\n4vZLYNuY2OnAbeHxbcAZ4fEC4HbgdWAD8BRwMnAwsB/wQHje6oxrMu91J/DBvLZeysI1xKqLv7ek\nclKKddxmEXSfEj7PCo8PATZmnLcRmJMlvimMEz4/Gx4PANsJumKlgim3NcR60j0sSSyhPd7uQtMF\nUG6/t6TqVurJCUPho+Da29tHjuPxOPF4vBgfqwo03OW+NLV0ZA2xpmRpxk86cL7wyun3llSZuru7\n6e7uzuncYizAewTwb8A7w9ePA3HgeYJu0J8Db2fHWLdrwudO4Crgj+E5R4fxc4B5wGfCc9oJJibU\nAZuBA7O0wU3mVZHcbF2SKs9EC/CWoqv0J8Ci8HgRcFdGfCGwN3AkwYSDBwgSvJcJxrvVAOcBd2e5\n1ycIJjtIVcOB85JUXQrdVXo78H5gJsFYtCsJKmp3EMwI3QCcFZ77aBh/lGC82sXs6Ea9GLgV2Ae4\nh6DSBnATsIZgOZCtBImfVDUcOC9J1aV69iq1q1QVKNsYt8Wxxa5ZKEkRNlFXqYmbFHFuti5JlcXE\nDYbaGhtpbGlhXnNzqdsiSZI0rokSt1IvB1I0K7q6aOvtBTB5kyRJkVSKWaUls7K3l3WpVKmbIUmS\ntJOedJolicSE51RV4gZQ29dX6iZIkiSN0pNOs7a1lRVdXROeVzVdpcMGGxoKct+edA9dHV3U9dcx\nUD9AY0ujA8QlSVJOujo6WBkO6ZpIVSVui2MxmpLJvN/XbYckSdJk1PX353Re1XSVLk0kaFq1qiAT\nE7o6ukYlbQAre1eyLrUu758lSZIqz0B9fU7nVU3FbXln565P2kNuOyRJKoSedJqujg7q+vsZqK93\nWasK1tjSQltv7y67S6smcSukctt2yPF2khR9w4PVM/8hd1mryjX8my5NpWDt2nHPq5oFeAu5c0I5\nbTuUdbxdrI3EqoTJWx75X8Gl529QXarx916SSGSdYbg0kShoL5JKzwV4C2w4IVqaWjqy7VBTsjR7\nRY433m5paqmJW574X8Gl529QXar19x5vsLrLWlU3E7c8mdc8rywSI8fbFV62Kdsre3tZmkpV9D8i\n5cTfoLpU6+893mD1Qi1rVU6qscKaKxO3ClNu4+0qkf8VXHr+BtWlWn/vbIPVC7WsVTmp1gprrkzc\nKkxjSyNtvW07j7dLNpWwVZWlmv8ruFz4G1SXYv/e5VLtyRysXtvXx2BDA03JZMUnL9VaYc2ViVuF\nKafxdpWqWv8ruJz4G1SXYv7e5VbtmdfcXHXJSrVWWHNl4laBymW8XaWq1v8KLif+BtWlmL+31Z7S\ns6I+MZcDkaRdKJeuMxVeezxO+y9+sXP8/e+nvbu7+A2qQtmqnotjsYLtflSOXA5EkvZQuXWdqbCs\n9pSeFfWJWXGTIs5qUGG5CGp1sdqjcmDFTapQVoMKz4HS1cVqj8qdiVueWPVQKTiQuvDsOqs+1TiT\nU9Fh4pYHVj1UKlaDCs+lRySVExO3PCi3qofVv+phNajw7DqTVE5M3PKgnKoeVv+qi9Wg4rDrTFK5\nMHHLg3KqepRb9U+FZTVIkqqLiVselFPVo5yqfyoOq0GSVD2mlPCzNwAPAw8CD4SxGcA64PdAFzAt\n4/wrgCeBx4HGjPiJwCPhe6sm+sDuAq16Pa+5mcSqVSxNJGh///tZmkiUbM2fQlf/CvUdVhu/x8nz\nO5w8v8P88HucPL/D3JUycRsC4sC7gJPC2JcIEre3Aj8NXwMcA5wdPjcB32LHwnQ3ABcBR4WPpvE+\nsJD/w5jX3Mzyzk7au7tZ3tlZsgpIY0sLbbHYqNjiWIz5ear++X9c+eH3OHl+h5Pnd5gffo+T53eY\nu1J3lY5dFfh04P3h8W1AN0HytgC4HXidoFL3FHAy8EdgP3ZU7FYDZwBVu5y5Y54kSapcpUzchoD7\ngEHg28A/A7OALeH7W8LXAIcAv8q4diMwhyCR25gR3xTGq5pjniRJqkyl3Kv0YGAzcCBB92gS+Akw\nPeOclwjGvaUIErfvhvEbgXsJqm/XAPPD+PuALwKnjfmsh4Dj8v0HSJIkFcB64Phsb5Sy4rY5fH4R\n+DHBOLctwGzgeYLE7oXwnE3A3IxrDyWotG0KjzPjm7J8VtY/XpIkKUpKNTlhKsHYNIA3EcwSfYSg\n4rYojC8C7gqPfwIsBPYGjiSYhPAAQYL3MsF4txrgvIxrJEmSlAdHEnRfPgT8lmCpDwi6Re8j+3Ig\niwkmJTwOJDLiw8uBPAV0FLTVkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJ\nkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJ\nkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJkiRJ\nkiRJkiSpUjUBjwNPApdneX8BsB54EPgN8IGM9zYAD4fvPVDQVkqSJFW5WuAp4AhgL+Ah4Ogx57wp\n4/id4fnDngZmFLB9kiRJZWFKqRsAnESQiG0AXge+T1Bhy/TnjON9gT+Neb+mUI2TJEkqF+WQuM0B\nns14vTGMjXUG8BhwL9CSER8C7gN+DfxDgdooSZJUcnWlbgBB4pWLu8LH+4A1wNvC+HuBzcCBwDqC\nsXK/zLzwuOOOG1q/fn1eGitJklRg64Hjs71RDhW3TcDcjNdzCapu4/klQcL55vD15vD5ReDHBF2v\no6xfv56hoSGuuuoqhoaGfEzi4Xfo91guD79Dv8Nyefg9+h3m+wEcN14SVA6J26+BowgmJ+wNnA38\nZMw5MXaMYzshfN4KTAX2C1+/CWgEHilgWyVJkkqmHLpKB4BLgLUEM0xvIhjL9unw/W8DHwc+RTB5\n4VVgYfjebOBH4XEd8F2gqyitliRJKrJySNwgmHBw75jYtzOOvxI+xvoD4/QBZxOPx3e7YRrN7zA/\n/B4nz+9w8vwO88PvcfL8DnNXLctoDIV9xpIkSWWtpqYGxsnRyqXiVhIzZsxg27ZtpW5G2Zk+fTov\nvfRSqZshSZLGqOqKW01NDVbidub3IklS6UxUcSuHWaWSJEnKgYmbJElSRJi4SZIkRYSJmyRJUkSY\nuJWhI444glmzZvHaa6+NxG688UZOPfVUAKZMmcK+++7Lfvvtx8yZM/nQhz7EHXfcMXLu29/+dm65\n5Zad7rtq1Sre8573AMGaOTfddFOB/xJJkpRPJm5l6o033mDVqlWjYuEsEwAefvhhXnnlFX7/+99z\n/vnnc8kll3D11VcDcP7557N69eqd7rlmzRrOP//8kXtl3k+SJJU/E7cs0uk0iUSCeDxOIpEgnU4X\n9R41NTV8/vOf57rrrmP79u0TnjtjxgzOPfdcbrjhBr785S+zbds2zj33XP7jP/6DZ555ZuS8Rx99\nlEceeYRzzjlnt/8WSZKKKR//Dleqql6AN5t0Ok1rayu9vb0jseHj5ubmot3j3e9+N/F4nOuuu47l\ny5fv8vzTTz+dgYEBHnjgARKJBKeeeipr1qyhra0NCKptzc3NzJgxI6fPlySpFPLxb2gls+I2RkdH\nx6j/sUDwP5hUKlXUe9TU1HD11VeTSqX405/+BDDhorh77bUXM2fOHNnxYNGiRaxZswYIul2/973v\nsWjRopw/X5KkUsjHv6GVzMRtjP7+/qzxvr6+ot4D4B3veAcf/ehHueaaa3Y5Hu3111/nxRdfHKmo\nnXnmmWzevJn777+f7u5uXnvtNf9LRZJU9vL1b2ilsqt0jPr6+qzxhoaGot5j2LJlyzjhhBP43Oc+\nN+F5d999N3V1dZx00kkATJ06lU984hOsXr2av/zlL5xzzjnU1flzS5LKWz7/Da1EVtzGaGlpIRaL\njYrFYjGSyWRR75F53dlnn82qVatGVd2Gu01feuklvvvd73LJJZfwpS99ienTp4+cs2jRIr7//e9z\n5513Zu0mdT9SSVK5yee/oZXIEswYw92JqVSKvr4+GhoaSCaTu9XNmI97ZLryyitHxqsNO+6446ip\nqWHvvffm+OOP5xvf+AYLFy4cdc68efOYNm0a++yzDyeeeOJO93U5EElSucn3v6GVplr+5R7KVl2q\nqamx6pSF34skSaUTFlay5mh2lUqSJEWEiZskSVJEmLhJkiSVgeEdIybi5ARJkqQSy7ZjRDZW3CRJ\nkkos244R2Zi4SZIkldh4O0aMZeImSZJUYuPtGDFWVY9xmz59uovQZpG5+4IkSSq8lpYWent7d9ld\nWi1ZS9YFeCVJkspFOp0mlUqxdu1aGCdHM3GTJEkqI1HYOaEJeBx4Erg8y/sLgPXAg8BvgA/sxrWS\nJEkVoRwqbrXAE8CHgE3AfwPnAI9lnPMm4M/h8TuBHwNvyfFasOImSZIiotwrbicBTwEbgNeB7xNU\n2DL9OeN4X+BPu3GtJElSRSiHxG0O8GzG641hbKwzCCpp9wItu3mtJElS5JXDciC59mHeFT7eB6wB\n3r47H9Le3j5yHI/Hicfju3O5JElFl06n6ejooL+/n/r6elpaWmhubi51s5Rn3d3ddHd353RuOYxx\nOwVoJ5hkAHAF8AZw7QTX9BJ0kx6V47WOcZMkRUq2vStjsRirVq0yeatw5T7G7dcECdgRwN7A2cBP\nxpwTY8cfcEL4vDXHayVJipxse1f29vaSSqVK1CKVg3LoKh0ALgHWEswSvYlgLNunw/e/DXwc+BTB\nBIRXgYW7uFaSpEgbb+/Kvr6+IrdE5aQcEjcIJhzcOyb27Yzjr4SPXK+VJFWQahzrNd7elQ0NDUVu\nicpJuSRukiRllW2s1/BxJSdv2faujMViJJPJErZKpVYOkxOKwckJkhRRiUSCrq6urPHOzs4StKh4\nhveu7Ovro6GhgWQyWdHJqgITTU6w4iZJKmvVPNarubnZRE2jlMOsUkmSxuVYL2kHEzdJUllraWkh\nFouNijnWS9XKrlJJUlkb7ip0rJfk5ARJkqSyUu47J0iSJCkHJm6SJEkRYeImSZIUESZukiRJEWHi\nJkmSFBEmbpIkSRFh4iZJkhQRJm6SJEkRYeImSZIUESZukiRJEWHiJkmSFBEmbpIkSRFh4iZJkhQR\nJm6SJGVIp9MkEgni8TiJRIJ0Ol3qJkkj6krdAEmSykU6naa1tZXe3t6R2PBxc3NzqZoljagpdQOK\nZGhoaKjUbZAklblEIkFXV1fWeGdnZwlapGpUU1MD4+RodpVKkhTq7+/PGu/r6ytyS6TsTNwkSQrV\n19dnjTc0NBS5JVJ2Jm6SJIVaWlqIxWKjYrFYjGQyWaIWSaM5OUGSpNDwBIRUKkVfXx8NDQ0kk0kn\nJqhslMvkhCbgG0AtcCNw7Zj3/x74IkF7XwE+AzwcvrcBeBkYBF4HTspyfycnSJKkSJhockI5JG61\nwBPAh4BNwH8D5wCPZZzzt8CjwHaCJK8dOCV872ngROClCT7DxE2SJEVCuc8qPQl4iqBy9jrwfWDB\nmHP+iyBpA7gfOHTM++WQgEqSJBVUOSRuc4BnM15vDGPjuQi4J+P1EHAf8GvgH/LeOkmSpDJRDpMT\ndqcP81TgQuC9GbH3ApuBA4F1wOPAL/PWOkmSpDJRDonbJmBuxuu5BFW3sY4F/plgjNu2jPjm8PlF\n4McEXa87JW7t7e0jx/F4nHg8PokmS5Ik5Ud3dzfd3d05nVsOY8PqCCYnfBB4DniAnScnHAb8DDgX\n+FVGfCrB5IZXgDcBXcCy8DmTkxMkSVIkTDQ5oRwqbgPAJcBagiTsJoKk7dPh+98GrgSmAzeEseFl\nP2YDPwpjdcB32TlpkyRJqgjlUHErBitukiQpEsp9ORBJkiTlwMRNkiQpIkzcJEmSIsLETZIkKSJM\n3CRJkiLCxE2SJCkiTNwkSZIiwsRNkiQpIkzcJEmSIsLETZJ2IZ1Ok0gkiMfjJBIJ0ul0qZskqUqV\nw16lklS20uk0ra2t9Pb2jsSGj5ubm0vVLElVyr1KJWkCiUSCrq6urPHOzs4StEhSpXOvUknaQ/39\n/VnjfX19RW6JJJm4SdKE6uvrs8YbGhqK3BJJMnGTpAm1tLQQi8VGxWKxGMlkskQtklTNnJwgSRMY\nnoCQSqXo6+ujoaGBZDLpxARJJeHkBEmSpDLi5ARJkqQKYOImRZyLw0pS9XCMmxRhLg4rSdXFipsU\nYR0dHaOSNggSt1QqVaIWFY+VRknVyIqbFGHVujislUZJ1cqKmxRh1bo4bDVXGiVVNxM3KcKqdXHY\naq00SpJdpVKEVevisNVaaZQkF+CVFDnZxrjFYjFWrVpV8UmrpMo30QK8VtwkRU61VholqVwqbk3A\nN4Ba4Ebg2jHv/z3wRYL2vgJ8Bng4x2vBipskSYqIiSpu5ZC41QJPAB8CNgH/DZwDPJZxzt8CjwLb\nCRK1duCUHK8FEzdJkhQR5b5X6UnAU8AG4HXg+8CCMef8F0HSBnA/cOhuXCtJklQRyiFxmwM8m/F6\nYxgbz0XAPXt4rSRJUmSVw+SE3enDPBW4EHjvHlwrSZIUaeWQuG0C5ma8nktQORvrWOCfCca4bdvN\na2lvbx85jsfjxOPxPW2vJElS3nR3d9Pd3Z3TueUwOaGOYILBB4HngAfYeYLBYcDPgHOBX+3mteDk\nBEmSFBHlvo7bAHAJsJZgluhNBInXp8P3vw1cCUwHbghjrxNMTBjvWkmSpIpTDhW3YrDiJkmSIqHc\nlwORJElSDkzcJEmSIsLETZIkKSJM3CRJkiIiH7NKX2X8hXCHgP3z8BmSJElVz1mlkiRJZaTQs0rv\nIVgYd9883EuSJEnjyEfi9h3go8DTwB3AmcDeebivJEmSMuSzq/RNwGnAQuBvCSpxtwNdefyMPWVX\nqSRJioSJukoLNcbtOOA24J0EW1GVmombJEmKhGLtnDAbaAH+E7gL6AROyOP9JUmSqlo+Km7/D0H3\n6NuBOwm6R/+L8ZcIKQUrbpIkKRImqrjlYx23vwW+DPwUeCMP95MkSVIW+RrjNhP4JEHVDeBRgsrb\n1jzdf7KsuEmSpEgo9Bi3o4FHgBOBJ4AngZPC2NsnuE6SJEm7IR8VtzuBHxCs4Zbp4wRVuI/n4TMm\ny4qbJEmKhEIvB/J74K178F4xmbhJkqRIKHRX6Z/38D1JkiTthnzMKj0Q+EeyZ4YH5uH+kiRJIj+J\n243AflniNcA/5+H+kiRJonBbXpUbx7gpr9LpNB0dHfT391NfX09LSwvNzc2lbpYkqQIUegHeq8aJ\nD2dKV+fhM6SykU6naW1tpbe3dyQ2fGzyJkkqpHxU3D7PzttbvQm4iGBh3jfl4TMmy4qb8iaRSNDV\n1ZU13tnZWYIWSZIqSaErbtdlHO9PsNH8BcD3gevzcH+prPT392eN9/X1FbklkpQ/DgGJhnwkbgBv\nBi4D/h5YDZwAbMvTvaWyUl9fnzXe0NBQ5JZIUn44BCQ68rGO23XAA8ArwLEEY95M2lSxWlpaiMVi\no2KxWIxkMlmiFknS5HR0dIxK2iBI3FKpVIlapPHko+L2j8BfgSXhI9MQQfdpxbPEXD2Gf9dUKkVf\nXx8NDQ0kk0l/b0mR5RCQ6MhH4paPql0T8A2glmBduGvHvP924BbgXUAbo8fObQBeBgaB1wk2uC8q\nS8zVp7m52d9WUsVwCEh05CPpmqxa4JsEydsxwDnA0WPO2QokGT0RYtgQECdI6oqetIElZklStDkE\nJDryNTlhMk4CniKonEEwG3UB8FjGOS+Gj/FKHCVdSNgSsyQpyhwCEh3lkLjNAZ7NeL0ROHk3rh8C\n7iPoKv02JdhmyxKzJCnqHAISDeXQVTrZlXHfS9BN+mHgs8D7Jt2i3WSJWZIkFUM5VNw2AXMzXs8l\nqLrlanP4/CLwY4Ku11+OPam9vX3kOB6PE4/Hd7OZ47PELEmS9lR3dzfd3d05nVsOm8zXAU8AHwSe\nI1gT7hxGj3Eb1k6wXtzwrNKpBJMbXiHYWqsLWBY+Z6qqLa9cmkSSpOgq9JZXkzUAXAKsJUjCbiJI\n2j4dvv9tYDbw3wRrwr0BtBLMQD0I+FF4Xh3wXXZO2qqKS5NIklS5yqHiVgxVU3FzA3RJkqJtoopb\nOUxOUB65NIkkSZXLxK3CuDSJJEmVy8Stwrg0iSRJlascJicoj1yaRJKkyuXkBEmSpDLi5ARJkqQK\nYOImSZIUESZukiRJEWHiJkmSFBEmbpIkSRFh4iZJkhQRJm6SJEkRYeImSZIUESZukiRJEWHiJkmS\nFBEmbpIkSRFh4iZJkhQRJm6SJEkRYeImScoqnU6TSCSIx+MkEgnS6XSpmyRVvbpSN0CSVH7S6TSt\nra309vaOxIaPm5ubS9UsqerVlLoBRTI0NDRU6jZIUmQkEgm6urqyxjs7O0vQIql61NTUwDg5WtV0\nlVrml6Tc9ff3Z4339fUVuSWSMlVNV2lXV5dlfknKUX19fdZ4Q0NDkVsiKVPVVNwgGJ+RSqVK3QxJ\nKnstLS3EYrFRsVgsRjKZLFGLJEEVVdyGWeaXpF0b7plIpVL09fXR0NBAMpm0x0IqsaqZnDB84MBa\nSZJUzpycELLML0mSoqxqukoTiYRlfkmSFGnlUnFrAh4HngQuz/L+24H/AvqAz+3mtQB0dnaatEmS\npEgrhzFutcATwIeATcB/A+cAj2WccyBwOHAGsA24fjeuBRfglSRJEVHuY9xOAp4CNgCvA98HFow5\n50Xg1+H7u3utJElSRSiHxG0O8GzG641hrNDXSpIkRUo5TE6YTB9mzte2t7ePHMfjceLx+CQ+VpIk\nKT+6u7vp7u7O6dxyGON2CtBOMMkA4ArgDeDaLOdeBbzKjjFuuV7rGDdJkhQJ5T7G7dfAUcARwN7A\n2cBPxjl37B+xO9dKkiRFWjl0lQ4AlwBrCWaJ3kQwK/TT4fvfBmYTzBjdn6Ci1gocQ1B9y3atJElS\nxSmHrtJisKtUkiRFQrl3lUqSJCkHJm6SJEkRYeImSZIUESZukiRJEWHiJkmSFBEmbpIkSRFh4iZJ\nkhQRJm6SJEkRYeImSdL/397dR9lR3gUc/4Yk7PJSiCmUdxrOggLFUrAN0SJZrHiBbQlYKcYiL619\nOT3dcCynckCUeGrVUpXCckBtqaf11BZLqaILIWC7CrQQeUswEAwBpISW4kuRmmyEJP7xm+HO3szc\nvbs7d/fO3u/nnD1777zdmec+88xvnueZ50oVYeAmSZJUEQZukiRJFWHgJkmSVBEGbpIkSRVh4CZJ\nklQRBm6SJEkVYeAmSZJUEQZukiSpaw0PD1Or1ejv76dWqzE8PDzTu9TUvJneAUmSpJkwPDzMJZdc\nwqZNm16flr4eGBiYqd1qas5M78A02blz586Z3gdJktRBarUaq1evzp2+atWqGdijMGfOHCiI0Wwq\nlSRJXWnbtm2500dHR6d5T1pn4CZJkrpST09P7vTe3t5p3pPWGbhJkqSutGLFCvr6+sZM6+vrY3Bw\ncIb2aHw+nCBJkrpS+gDC0NAQo6Oj9Pb2Mjg42LEPJoAPJ0iSJHUUH06QJEmaBQzcJKlCqjZYqKRy\n2c7p0E8AAA6CSURBVMdNkiqiioOFSipXp9S4nQ5sADYClxUsc10yfy1wQmb6s8A64BFgTft2UZJm\n1nXXXTcmaIMI3IaGhmZojyRNt04I3OYC1xPB27HAcuCYhmXOBI4EjgI+DNyYmbcT6CeCucVt3ld1\nMJuQNNtVcbBQSeXqhKbSxcBTRM0ZwNeAZcATmWXOAr6UvH4AWAAcALyYTOuWp2NVwCYkdYMqDhYq\nqVydUON2CPC9zPvnk2mtLrMTuBt4EPhQm/ZRHc4mJHWDKg4WKqlcnVDj1uoAa0W1aicDLwD7A3cR\nfeXuKWG/VCE2IakbVHGwUEnl6oTAbTNwWOb9YUSNWrNlDk2mQQRtAC8B3ySaXncJ3FauXPn66/7+\nfvr7+6ewy+o0NiGpWwwMDBioSbPMyMgIIyMjLS3bCX3D5gFPAu8igrA1xAMK2T5uZwIfT/4vAT6X\n/N+TeLjhFWAvYDXwe8n/LH85YZbL6+PW19fHtdde60VOklQpzX45oRNq3F4jgrI7iSDsJiJo+0gy\n/8+B24mg7Sngf4GLk3kHArcmr+cBX2HXoE1dwCYkSVI36IQat+lgjZskSaoEf6tUkiRpFjBwkyRJ\nqggDN0mSpIowcJMkSaoIAzdJkqSKMHCTJEmqCAM3SZKkijBwkyRJqggDN0ktGx4eplar0d/fT61W\nY3h4eKZ3SZK6Sif85JWkCsj7Pdj0tT8tJknTw5+8ktSSWq3G6tW7/hRwrVZj1apVM7BHkjQ7+ZNX\nkqZs27ZtudNHR0eneU8kqXsZuGlK7PPUPXp6enKn9/b2TvOeSFL3so+bJs0+T91lxYoVbNq0acz3\n3dfXx+Dg4AzulSR1F/u4adLs89R9hoeHGRoaYnR0lN7eXgYHBw3SJalkzfq4WeOmSbPPU/cZGBgw\nUJOkGWQfN02afZ4kSZpeBm6atBUrVtDX1zdmmn2eJElqH5tKNWlpk5l9niRJmh4+nCBJktRBHIBX\nkiRpFjBwkyRJqggDN0mSpIowcJMkSaoIAzdJkqSKMHCTJEmqCAM3SZKkiuiUwO10YAOwEbisYJnr\nkvlrgRMmuK4kSVLldULgNhe4ngjAjgWWA8c0LHMmcCRwFPBh4MYJrPu6kZGREne7O5mG5TAdp840\nnDrTsBym49SZhq3rhMBtMfAU8CzwKvA1YFnDMmcBX0pePwAsAA5scd3XmTGmzjQsh+k4dabh1JmG\n5TAdp840bF0nBG6HAN/LvH8+mdbKMge3sK4kSdKs0AmBW6s/Itotv6sqSZKUqxOCoSXASqKfGsDl\nwA7gM5ll/gwYIZpCIR5GWAoc0cK6AI8Cx5e615IkSe2xFnjbTO9EkXnAJmARsDsRZOU9nHB78noJ\ncP8E1pUkSVKJzgCeJB40uDyZ9pHkL3V9Mn8tcOI460qSJEmSJEnd6cclbKMHuJkYgPd+4M3J9FOB\nRzJ/W4lhRRp9FniCqMm7Fdg3M+/yZLsbgF/KTP808BzwSsE+vZfoa3diwfyylZGOnwDWE+lwN3B4\nMn0OMfjxeuBx4NqC9auajtsZm08Ob7LsCPAz42xvIfBt4piGGuZdDDxGpNEdwBtz1n9/Mn8dcB/w\n1sy8osGmzyW+n+3kp9XhRB65dJx9n6wdwF9l3s8DXgL+forb/SLwIpFmWeMdL0TfkO8A/0qk5/sy\n844ghhXaSPSbnZ9MPxr4LjBKflrNJfLIVI+rmbOJ9PypEra1ELgL+DdgNTGMUjq9KI9mfYXIb48B\nNxHfa6poQPSi7yx1KXF8CydwHBNVZhoWlWunAQ8S5+mDxPUmT9XSsF3nMjTPd6uIrk7riXSaz66q\nWDbOSkUX7In4GHBD8vo86g8wZP0E8J9Ab86806g/XftHyR/EgL6PEhloEdEUmz7MsZgYRy5v/98A\n/DNx0ZiuwK2MdOynnj4fpZ6O/cC9xLHvRhzX0pz1q5qOE0m7bzN+4LYn8E6imT9bOO1O5MG0sP0M\ncFXO+j9L/eJwOvX+nHOJtFtEpGW2P+fRwE8m+5eXVrcQNzftKpxeAR6mnn/OIAKc2yawjXk5036e\nuKA1XsDGO16Iwbr7ktcHAS8A+yTv/4Z6IHcjkd8B9gfeDvw++Wn1CeJCPJHjmqibk+2vnMS6jaME\nXA38VvL6MurnZFEebXRG5vVfU0+nbJ/jk6jnUSj+zgAOIy7Qz9DewK3MNCwq195GlF0AbyGGospT\ntTQs41wu0izf7Z15fQtwfs76VSwb26IThgPZi6jheYiIpNNasUXEnc5fEHfNd5IfeGUH5/0G8K6c\nZc4lTpLRnHl3EXcZEHfhhyavlwFfJQb2fZbIGCcl89YAPyg4nk8RJ/c2pvep3amm4wj19Mmmww+J\noKMH2IM4MfKOfbakI0RwNkLcSa+iXkAD/DpRkD0GvCNn3S3E3eC2humvAf9NFFBziCBic8763wVe\nTl5n07HZYNMbiFqVPGcDTxO1pe10OzCQvF5OfOfZAP07xAXhPqIgBbiIuCD8I5F/Gt1DpFmjZseb\n2kg8uATwfSIf75/s06lEgQ1RdpydvH6J+M5fzdneocTF9gu0Lz/uTZwbHyduQlP9xE3MPxDHfmNm\nH34M/DFxsVrSsL1s2Zg9zqI82uiOzOt/oT5G5jLyB0SH4u8M4E+pB5LtUnYaFpVrj1Ivux6nXjY2\nqmIaTuZc/ifGjtxwL/DTDdttlu/SVqP5xPXmP3KWqWrZWLpOCNy2AucQF8tfAP4kM+9I4qGE44Af\nEU1njbKD875GfLGNdyK/SmS+8XyA+l3QwYy9i2plcN8Tk2XSbbQ6Rl0ZppqOWR+kfgyPE80s3ycC\njVXEwyDNVCkd96DeTPoNouZniEijtwN/STTpQhReexB3wx8jmjSKNO7zDuASInjeTNwRNlsfxn4P\nrQxU3WhvopBfOc5yZbiZOM96iAL7gcy8J4hahBOJWsY/yMw7gUjroqamMiwmLgabiObpH1G/GG+m\ntUG7rwE+mVmvHZYR59dzRBCZrSF4BxGMHEvUJP5yMn1PouYhbRrOOoBociP5f0DD/FbPq/lEDciq\n5P1kBj5fliy3rsXPnKyy0zArW65lvZe4Yc4L+FNVSsPJnMs3ETdiEMFcD8VNvUX57k4in26lnk5F\nqlQ2li6veWK67Qb8IZEZdhAZ+k3JvGeoZ9KHiNqjiTqICFjuHGe53wb+j6jOLtKsoNuNuBu6MDNt\nOmuKykrH84mT8jeT96cQF9VDiOO5i0jLewvWr1o6bmVs/5LjiKaPu5P3c4lmNoj9Tm8A7iFqzfYB\n/qeFz9mH6NNyPPF9DBF9/z5dsPypxIXinZnPnqiVRMCxhfbnxceIfLUcGG6YtwD4MnEDsZOx5c5q\nIpBql4OSz75gCtt4N1Fj9whRc9Muy4nvC+DryfuHk/driBoFiDx4MnGjsT35P56dTP4G6AaiRuW+\nzLTG/NRs23sCVxDNjkXrl6VdaVhUrr2FaBk4bZc1xqpSGk7kXE5rGW8Bfoe4ufkAccM7UTXqfdYv\npF4j2ahqZWPpOiFwez+wHxEsbCcuamlTXrZKdTtR29FoM9HB8AXiePYF/isz/31Ep9LtTfbhIqIZ\nJNvMupnoT5A6lPymrdQbiJN4JHl/INEM9B7qBUc7TTUdAX6RKBxOoX73uISo7t+SvL+D6GuQF7hd\nRPXTcQ7RmfXnWly+1ULjGOI7eSZ5/3XGdqLNeivweaIfR9pk0piOh1Hcrya1mKgNuJoocHcQgeoN\nzVaagtuIJqelRLNk6lNEc+g5xMNDI5l5W2iffYimsSuIizZEP8MFxA3CDsbPjxB54Swib/cm251q\nMNhoIXFBOo7IU3OT/59M5mfz2RzqNX+jFOfBF4nz5wdEAPvDSezXVUQt5Ycy0yZ6TvcRgcDazPIP\nEflzMvtUpB1pCPnlGsRx3Ep0n3iGYlVKw9REz+UtxE392UTXpMn2S95GBNEnkR+4VbVsLFUnNJXu\nS2S87cRJ9+bmi+/iNuq1M79CZKqstI2+yOnEib2MsX3gbiOqi3cnnkI7inrhn+dlIoMfkfzdz/QF\nGzD1dDyB+IWK9zC2f0H6KxVziburpeT3CZgt6fhk8vlpX5f5RNMKRGGf9ps5magpKnq4ofEu7mmi\no+x+yfvTyE/Hw4mLwflEv43Ug0TaLSLS8jzyOwxnP/cU6un4OaJ2r50F0xeJO9n1DdP3oV5reXHJ\nn1l0t7w78E0iwLo1M30n0VH53OT9hcDfjrPNK4iLwRFEXv4W5QZtEGXXl4nv9wgiHzxD1KBDXGgW\nEWX2eRTXeGdly8ZWjrPRbxBPgf9aznbT419CnAcvUuwxopk2zYvPExf2sgOOdqRhUbm2gKiNuozo\ne1WkammYmsy5/AWiVWEN9b5oeRrz3V7EjQVE5cu7idrtRlUuG2eNeUSA8EaiX8E6IrOsJ76gRYxt\ny78U+N2c7fQQT4mlw4EsysxbxNi27zwbgX+n3s8p++VdQWSQDUQ1burqZLuvJf/z9qvZE29lKisd\n7yL6saXpkC3kryH6Zq0n7sLyVDUd85o5jyeaNR4ljvuDmX25hggi1xF94PI8S9TsvEIc19HJ9Auo\nDwfyd8TTzo0+n6ybpmM2yC0abPqc5HO2ErUr2Q7RqauIpyLbIS8Nl1IvPJcQ+/0wccf+dDL9QqKg\nL/JV4iKxjTi+9ELRyvGeTzRtZYd6SYcPyA4HcjP15p4Dk+2+TNzNP8fYp90aj6tM32LsUDkAg8R5\ntJTIj2nH+uy51ayZfiHR5N84HAiMzaPPUc+jWa8SaZSm35WZeUUDohd9Z1lP054nItuRhkXl2pVE\nh/ps/tovZ/2qpeFkz+XUE+z6HWQ9y65l45uIci4d6uOz5N9UVLFsnHWOZ+wj0Joc01Ga3fpp79hx\n3aAf07DdDmb8B9dUgplqKv0o0cnzyvEWVFOmozT7TeXBAgXTsL0uICoQrpjpHZEkSZIkSZIkSZIk\nSZIkSZIkSZIkSZIkSZKkafT/kDcAn1J74VgAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x1084ca950>"
       ]
      }
     ],
     "prompt_number": 13
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 13
    }
   ],
   "metadata": {}
  }
 ]
}