RCC_footprint.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## For a given ga-metadata.yaml - draw the polygon as a shapely object"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from datetime import timedelta, datetime\n",
    "import ephem\n",
    "from osgeo import ogr\n",
    "\n",
    "import fiona\n",
    "import folium\n",
    "from IPython.display import display\n",
    "import geopandas\n",
    "%matplotlib inline\n",
    "\n",
    "from shapely.geometry import mapping\n",
    "from shapely.geometry import shape\n",
    "from shapely.geometry import Polygon\n",
    "from shapely.geometry import MultiPolygon\n",
    "from shapely.geometry import LineString\n",
    "import yaml\n",
    "import pandas"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get information from the RCC metadata AOS LOS and Platform"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'lineage': {'machine': {}, 'source_datasets': {}}, 'product_type': 'satellite_telemetry_data', 'format': {'version': 0, 'name': 'RCC'}, 'image': {'bands': {}}, 'creation_dt': datetime.datetime(2015, 11, 23, 0, 24, 19, 191972), 'instrument': {'name': 'TM', 'operation_mode': 'SAM'}, 'ga_label': 'LS5_TM_STD-RCC_P00_L5TT1986227235010ASA111_0_0_19860815T235010Z19860815T235627', 'platform': {'code': 'LANDSAT_5'}, 'ga_level': 'P00', 'size_bytes': 3971638821, 'usgs': {'interval_id': 'L5TT1986227235010ASA111'}, 'checksum_path': 'package.sha1', 'id': '896605c4-9178-11e5-879d-ac162d791418', 'acquisition': {'los': datetime.datetime(1986, 8, 15, 23, 56, 27), 'platform_orbit': 13067, 'aos': datetime.datetime(1986, 8, 15, 23, 50, 10), 'groundstation': {'eods_domain_code': '002', 'code': 'ASA', 'label': 'Alice Springs'}}}\n"
     ]
    }
   ],
   "source": [
    "with open('ga-metadata.yaml', 'r') as stream:\n",
    "    try:\n",
    "        meta = (yaml.load(stream))\n",
    "        print(meta)\n",
    "    except yaml.YAMLError as exc:\n",
    "        print(exc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'LANDSAT_5'"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "meta['platform']['code']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1986"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "meta['acquisition']['aos'].year"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "datetime.datetime(1986, 8, 15, 23, 56, 27)"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "meta['acquisition']['los']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def dms2dec(dms):\n",
    "    strdms = str(dms)\n",
    "    dms = [float(x) for x in strdms.split(':')]\n",
    "    if dms[0] < 0:\n",
    "        dec = dms[0]-(dms[1]/60)-(dms[2]/3600)\n",
    "    else:\n",
    "        dec = dms[0]+(dms[1]/60)+(dms[2]/3600)\n",
    "    return dec"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## From the archive TLEs find a reasonably close match to determine where the satellite was and the sub satellite position - use 7 days as the limit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "LS5tle = 'tle/LS5_ARCHIVE.txt'\n",
    "platform = meta['platform']['code']\n",
    "aos = meta['acquisition']['aos']\n",
    "los = meta['acquisition']['los']\n",
    "deltadays = 30 # Tolerance for time delta between TLE epoch and AOS\n",
    "track = []\n",
    "with open(LS5tle, 'r') as f:\n",
    "    data = f.readlines()\n",
    "    index = 0\n",
    "    found = False\n",
    "    \n",
    "    for line in data:\n",
    "        if (line[0]=='1'):\n",
    "            #print line[18:20]\n",
    "            index +=1\n",
    "            #print data[index]\n",
    "            sat = ephem.readtle(platform,line,data[index])\n",
    "            try:\n",
    "                sat.compute(str(aos))\n",
    "                epoch =  datetime.strptime(str(sat._epoch), \"%Y/%m/%d %H:%M:%S\")\n",
    "                absolute = abs(aos-epoch)\n",
    "                delta = timedelta(days=deltadays)\n",
    "                timestep = timedelta(seconds=5)\n",
    "                if absolute < delta:\n",
    "                    while aos < los:\n",
    "                        sat.compute(str(aos))\n",
    "                        track.append([dms2dec(sat.sublong),dms2dec(sat.sublat)])\n",
    "                        aos+=timestep\n",
    "                    sat.compute(str(los))\n",
    "                    track.append([dms2dec(sat.sublong),dms2dec(sat.sublat)])\n",
    "                    # GO DO GETPASSES IGNORE OBSERVER\n",
    "                    satellite_track = LineString(track)\n",
    "                    break\n",
    "            except ValueError:\n",
    "                pass            \n",
    "        else:\n",
    "            index +=1\n",
    "\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"136.812317778 -39.8480155556 7.60575333333 24.35367\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-55.3423611111)\"><polyline fill=\"none\" stroke=\"#66cc99\" stroke-width=\"0.4870734\" points=\"143.516083333,-16.3963333333 143.448472222,-16.6962222222 143.380722222,-16.9960833333 143.312833333,-17.2959166667 143.244777778,-17.59575 143.176555556,-17.8955833333 143.108194444,-18.1953611111 143.039638889,-18.4951666667 142.970944444,-18.7949166667 142.902083333,-19.0946944444 142.833027778,-19.3944166667 142.763833333,-19.6941388889 142.694416667,-19.9938611111 142.624833333,-20.2935277778 142.555055556,-20.5931944444 142.485083333,-20.8928611111 142.414916667,-21.1925 142.344583333,-21.4921111111 142.274,-21.7917222222 142.203222222,-22.0913055556 142.132222222,-22.3908611111 142.061027778,-22.6903888889 141.989611111,-22.9899166667 141.917972222,-23.2894166667 141.846083333,-23.5889166667 141.773972222,-23.8883611111 141.701638889,-24.1878055556 141.629055556,-24.4872222222 141.556222222,-24.7866388889 141.483138889,-25.086 141.409777778,-25.3853611111 141.336194444,-25.6846944444 141.262333333,-25.984 141.188194444,-26.2832777778 141.113777778,-26.5825555556 141.039083333,-26.8818055556 140.964111111,-27.181 140.888833333,-27.4801944444 140.813277778,-27.7793611111 140.737416667,-28.0785 140.66125,-28.3776111111 140.58475,-28.6767222222 140.507972222,-28.9757777778 140.430833333,-29.2748055556 140.353388889,-29.5738333333 140.275611111,-29.8728055556 140.197472222,-30.17175 140.118972222,-30.4706666667 140.040138889,-30.7695833333 139.960944444,-31.0684444444 139.881388889,-31.3672777778 139.801472222,-31.6660833333 139.721138889,-31.9648611111 139.640444444,-32.2636111111 139.559333333,-32.5623055556 139.477833333,-32.8609722222 139.395916667,-33.1596388889 139.313611111,-33.45825 139.230861111,-33.7568333333 139.147666667,-34.0553611111 139.064055556,-34.3538611111 138.979972222,-34.6523611111 138.895444444,-34.9507777778 138.810444444,-35.2491944444 138.724972222,-35.5475555556 138.639055556,-35.8458611111 138.552611111,-36.1441666667 138.465666667,-36.4423888889 138.378222222,-36.7406111111 138.29025,-37.0387777778 138.20175,-37.3368888889 138.112722222,-37.6349722222 138.023138889,-37.9330277778 137.933,-38.231 137.842277778,-38.5289722222 137.751,-38.8268611111 137.714305556,-38.9460277778\" opacity=\"0.8\" /></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.linestring.LineString at 0x7f69c68d3e90>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "satellite_track"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Intersect the sub satellite track with the WRS2 map to show which path rows constitute the data segment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "shapes = []\n",
    "\n",
    "for pol in fiona.open('reference/wrs2_descending.shp'):\n",
    "    if satellite_track.intersects(shape(pol['geometry'])):\n",
    "        if shape(pol['geometry']).is_valid:\n",
    "            shapes.append(shape(pol['geometry']).buffer(0))\n",
    "\n",
    "satellite_scenes = MultiPolygon(shapes)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Self-intersection at or near point 136.93619298390487 -38.002328054314034\n",
      "WARNING:shapely.geos:Self-intersection at or near point 136.93619298390487 -38.002328054314034\n"
     ]
    },
    {
     "data": {
      "image/svg+xml": [
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"100.0\" height=\"100.0\" viewBox=\"135.460344 -40.801956 10.190312 26.774712\" preserveAspectRatio=\"xMinYMin meet\"><g transform=\"matrix(1,0,0,-1,0,-54.8292)\"><g><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 144.658419678,-15.2661156794 L 144.659,-15.2662 L 144.37,-16.5812365269 L 144.325,-16.786 L 144.324363463,-16.7859075814 L 142.61,-16.537 L 142.627160662,-16.46192334 L 142.953838838,-15.0327298548 L 142.957,-15.0189 L 144.062675989,-15.1795543314 L 144.37,-15.2242084019 L 144.63693146,-15.2629934489 L 144.658419678,-15.2661156794 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 144.324363463,-16.7859075814 L 144.020689226,-18.1550548684 L 144.004,-18.2303 L 142.277,-17.9793 L 142.294289626,-17.9043420604 L 142.627,-16.4619 L 142.627160662,-16.46192334 L 144.341,-16.7109 L 144.324363463,-16.7859075814 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 143.697816804,-19.5988733993 L 143.682175868,-19.6687469835 L 143.681,-19.674 L 143.680986228,-19.6739979998 L 141.939157684,-19.4210229012 L 141.939,-19.421 L 142.294,-17.9043 L 142.294289626,-17.9043420604 L 144.020689226,-18.1550548684 L 144.021,-18.1551 L 143.697816804,-19.5988733993 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 141.602496651,-20.8347194642 L 141.939157684,-19.4210229012 L 141.957,-19.3461 L 143.697816804,-19.5988733993 L 143.698,-19.5989 L 143.682175868,-19.6687469835 L 143.680986228,-19.6739979998 L 143.354,-21.1173 L 143.353788183,-21.1172692396 L 141.596,-20.862 L 141.602496651,-20.8347194642 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 143.022590136,-22.5599404716 L 141.248579534,-22.3022841712 L 141.248,-22.3022 L 141.602496651,-20.8347194642 L 141.614,-20.7871 L 143.371,-21.0422 L 143.353788183,-21.1172692396 L 143.023,-22.56 L 143.022590136,-22.5599404716 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 143.022590136,-22.5599404716 L 142.70542001,-23.9270157995 L 142.688,-24.0021 L 140.895,-23.7416 L 140.913362302,-23.6668525974 L 141.248579534,-22.3022841712 L 141.267,-22.2273 L 143.04,-22.4849 L 143.022590136,-22.5599404716 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 142.70542001,-23.9270157995 L 142.706,-23.9271 L 142.365716144,-25.3684587688 L 142.348,-25.4435 L 140.535,-25.1801 L 140.537157372,-25.1714630907 L 140.913,-23.6668 L 140.913362302,-23.6668525974 L 142.70542001,-23.9270157995 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 142.020940131,-26.8092913022 L 142.016608637,-26.8273773818 L 142.003,-26.8842 L 142.002811907,-26.8841726681 L 140.169,-26.6177 L 140.188015239,-26.543002214 L 140.537157372,-25.1714630907 L 140.554,-25.1053 L 142.365716144,-25.3684587688 L 142.366,-25.3685 L 142.020940131,-26.8092913022 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 139.814409595,-27.9796594873 L 140.188,-26.543 L 140.188015239,-26.543002214 L 142.020940131,-26.8092913022 L 142.021,-26.8093 L 142.016608637,-26.8273773818 L 142.002811907,-26.8841726681 L 141.653,-28.3242 L 141.652467317,-28.3241226204 L 139.795,-28.0543 L 139.814409595,-27.9796594873 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 141.652467317,-28.3241226204 L 141.3145435,-29.6884337081 L 141.296,-29.7633 L 139.413,-29.4898 L 139.424568678,-29.4462313792 L 139.814,-27.9796 L 139.814409595,-27.9796594873 L 141.671,-28.2493 L 141.652467317,-28.3241226204 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 141.3145435,-29.6884337081 L 141.315,-29.6885 L 140.950970191,-31.1267504106 L 140.932,-31.2017 L 139.023,-30.9242 L 139.043258837,-30.8496375984 L 139.424568678,-29.4462313792 L 139.433,-29.4152 L 141.3145435,-29.6884337081 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 140.950970191,-31.1267504106 L 140.952,-31.1269 L 140.580364437,-32.5642076678 L 140.561,-32.6391 L 138.623,-32.3574 L 138.634618978,-32.3156878699 L 139.043,-30.8496 L 139.043258837,-30.8496375984 L 140.950970191,-31.1267504106 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 140.201733175,-34.0007612175 L 140.18983239,-34.0458350698 L 140.182,-34.0755 L 140.181824392,-34.0754744748 L 138.213096088,-33.7893139667 L 138.213,-33.7893 L 138.634618978,-32.3156878699 L 138.644,-32.2829 L 140.580364437,-32.5642076678 L 140.581,-32.5643 L 140.201733175,-34.0007612175 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 137.810335603,-35.1576086171 L 138.213096088,-33.7893139667 L 138.235,-33.7149 L 140.201733175,-34.0007612175 L 140.202,-34.0008 L 140.18983239,-34.0458350698 L 140.181824392,-34.0754744748 L 139.81421338,-35.436085722 L 139.794,-35.5109 L 137.792,-35.2199 L 137.810335603,-35.1576086171 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 139.416725082,-36.8705600255 L 139.405125519,-36.9123350623 L 139.396,-36.9452 L 139.395783719,-36.9451685506 L 137.359,-36.649 L 137.365330762,-36.6280806585 L 137.810335603,-35.1576086171 L 137.814,-35.1455 L 139.81421338,-35.436085722 L 139.815,-35.4362 L 139.416725082,-36.8705600255 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 136.936192984,-38.0023280543 L 137.365330762,-36.6280806585 L 137.382,-36.5747 L 139.416725082,-36.8705600255 L 139.417,-36.8706 L 139.405125519,-36.9123350623 L 139.395783719,-36.9451685506 L 139.009256001,-38.303691844 L 138.988,-38.3784 L 136.913,-38.0766 L 136.936192984,-38.0023280543 z\" /><path fill-rule=\"evenodd\" fill=\"#ff3333\" stroke=\"#555555\" stroke-width=\"0.53549424\" opacity=\"0.6\" d=\"M 138.589835628,-39.7358760793 L 138.584197753,-39.7550920484 L 138.568,-39.8103 L 138.567527745,-39.8102313043 L 136.452212198,-39.5025308669 L 136.452,-39.5025 L 136.936,-38.0023 L 136.936192984,-38.0023280543 L 139.009256001,-38.303691844 L 139.01,-38.3038 L 138.589835628,-39.7358760793 z\" /></g></g></svg>"
      ],
      "text/plain": [
       "<shapely.geometry.multipolygon.MultiPolygon at 0x7f69c691aa90>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "satellite_scenes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><iframe src=\"data:text/html;base64,CiAgICAgICAgPCFET0NUWVBFIGh0bWw+CiAgICAgICAgPGhlYWQ+CiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICA8bWV0YSBodHRwLWVxdWl2PSJjb250ZW50LXR5cGUiIGNvbnRlbnQ9InRleHQvaHRtbDsgY2hhcnNldD1VVEYtOCIgLz4KICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9sZWFmbGV0LzAuNy4zL2xlYWZsZXQuanMiPjwvc2NyaXB0PgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgPHNjcmlwdCBzcmM9Imh0dHBzOi8vYWpheC5nb29nbGVhcGlzLmNvbS9hamF4L2xpYnMvanF1ZXJ5LzEuMTEuMS9qcXVlcnkubWluLmpzIj48L3NjcmlwdD4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxzY3JpcHQgc3JjPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9qcy9ib290c3RyYXAubWluLmpzIj48L3NjcmlwdD4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxzY3JpcHQgc3JjPSJodHRwczovL3Jhd2dpdGh1Yi5jb20vbHZvb2dkdC9MZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy8yLjAvZGV2ZWxvcC9kaXN0L2xlYWZsZXQuYXdlc29tZS1tYXJrZXJzLmpzIj48L3NjcmlwdD4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9sZWFmbGV0Lm1hcmtlcmNsdXN0ZXIvMC40LjAvbGVhZmxldC5tYXJrZXJjbHVzdGVyLXNyYy5qcyI+PC9zY3JpcHQ+CiAgICAgICAgCiAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvbGVhZmxldC5tYXJrZXJjbHVzdGVyLzAuNC4wL2xlYWZsZXQubWFya2VyY2x1c3Rlci5qcyI+PC9zY3JpcHQ+CiAgICAgICAgCiAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vY2RuanMuY2xvdWRmbGFyZS5jb20vYWpheC9saWJzL2xlYWZsZXQvMC43LjMvbGVhZmxldC5jc3MiIC8+CiAgICAgICAgCiAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vYm9vdHN0cmFwLzMuMi4wL2Nzcy9ib290c3RyYXAubWluLmNzcyIgLz4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC10aGVtZS5taW4uY3NzIiAvPgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2ZvbnQtYXdlc29tZS80LjEuMC9jc3MvZm9udC1hd2Vzb21lLm1pbi5jc3MiIC8+CiAgICAgICAgCiAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vcmF3Z2l0LmNvbS9sdm9vZ2R0L0xlYWZsZXQuYXdlc29tZS1tYXJrZXJzLzIuMC9kZXZlbG9wL2Rpc3QvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuY3NzIiAvPgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9sZWFmbGV0Lm1hcmtlcmNsdXN0ZXIvMC40LjAvTWFya2VyQ2x1c3Rlci5EZWZhdWx0LmNzcyIgLz4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvbGVhZmxldC5tYXJrZXJjbHVzdGVyLzAuNC4wL01hcmtlckNsdXN0ZXIuY3NzIiAvPgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL3Jhdy5naXRodWJ1c2VyY29udGVudC5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIiAvPgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAgICAgPHN0eWxlPgoKICAgICAgICAgICAgaHRtbCwgYm9keSB7CiAgICAgICAgICAgICAgICB3aWR0aDogMTAwJTsKICAgICAgICAgICAgICAgIGhlaWdodDogMTAwJTsKICAgICAgICAgICAgICAgIG1hcmdpbjogMDsKICAgICAgICAgICAgICAgIHBhZGRpbmc6IDA7CiAgICAgICAgICAgICAgICB9CgogICAgICAgICAgICAjbWFwIHsKICAgICAgICAgICAgICAgIHBvc2l0aW9uOmFic29sdXRlOwogICAgICAgICAgICAgICAgdG9wOjA7CiAgICAgICAgICAgICAgICBib3R0b206MDsKICAgICAgICAgICAgICAgIHJpZ2h0OjA7CiAgICAgICAgICAgICAgICBsZWZ0OjA7CiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgIDwvc3R5bGU+CiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAgICAgPHN0eWxlPiAjbWFwXzQwMTY0NjEzY2IyZjQxYjVhYmQ5YjBkZjIyZmY1NmI2IHsKICAgICAgICAgICAgICAgIHBvc2l0aW9uIDogcmVsYXRpdmU7CiAgICAgICAgICAgICAgICB3aWR0aCA6IDEwMC4wJTsKICAgICAgICAgICAgICAgIGhlaWdodDogMTAwLjAlOwogICAgICAgICAgICAgICAgbGVmdDogMC4wJTsKICAgICAgICAgICAgICAgIHRvcDogMC4wJTsKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgPC9zdHlsZT4KICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICA8L2hlYWQ+CiAgICAgICAgPGJvZHk+CiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAgICAgPGRpdiBjbGFzcz0iZm9saXVtLW1hcCIgaWQ9Im1hcF80MDE2NDYxM2NiMmY0MWI1YWJkOWIwZGYyMmZmNTZiNiIgPjwvZGl2PgogICAgICAgIAogICAgICAgIAogICAgICAgIAogICAgICAgIDwvYm9keT4KICAgICAgICA8c2NyaXB0PgogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCgogICAgICAgICAgICB2YXIgc291dGhXZXN0ID0gTC5sYXRMbmcoLTkwLCAtMTgwKTsKICAgICAgICAgICAgdmFyIG5vcnRoRWFzdCA9IEwubGF0TG5nKDkwLCAxODApOwogICAgICAgICAgICB2YXIgYm91bmRzID0gTC5sYXRMbmdCb3VuZHMoc291dGhXZXN0LCBub3J0aEVhc3QpOwoKICAgICAgICAgICAgdmFyIG1hcF80MDE2NDYxM2NiMmY0MWI1YWJkOWIwZGYyMmZmNTZiNiA9IEwubWFwKCdtYXBfNDAxNjQ2MTNjYjJmNDFiNWFiZDliMGRmMjJmZjU2YjYnLCB7CiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjZW50ZXI6Wy0zMCwxNTBdLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgem9vbTogNCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG1heEJvdW5kczogYm91bmRzLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbGF5ZXJzOiBbXSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNyczogTC5DUlMuRVBTRzM4NTcKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB9KTsKICAgICAgICAgICAgCiAgICAgICAgCiAgICAgICAgCiAgICAgICAgICAgIAogICAgICAgICAgICB2YXIgdGlsZV9sYXllcl9kZjVkMTUzYThhN2U0Y2UyYWI3MmEyYTA3Zjg3YjZiYyA9IEwudGlsZUxheWVyKAogICAgICAgICAgICAgICAgJ2h0dHBzOi8ve3N9LnRpbGUub3BlbnN0cmVldG1hcC5vcmcve3p9L3t4fS97eX0ucG5nJywKICAgICAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgICAgICBtYXhab29tOiAxOCwKICAgICAgICAgICAgICAgICAgICBtaW5ab29tOiAxLAogICAgICAgICAgICAgICAgICAgIGF0dHJpYnV0aW9uOiAnRGF0YSBieSA8YSBocmVmPSJodHRwOi8vb3BlbnN0cmVldG1hcC5vcmciPk9wZW5TdHJlZXRNYXA8L2E+LCB1bmRlciA8YSBocmVmPSJodHRwOi8vd3d3Lm9wZW5zdHJlZXRtYXAub3JnL2NvcHlyaWdodCI+T0RiTDwvYT4uJywKICAgICAgICAgICAgICAgICAgICBkZXRlY3RSZXRpbmE6IGZhbHNlCiAgICAgICAgICAgICAgICAgICAgfQogICAgICAgICAgICAgICAgKS5hZGRUbyhtYXBfNDAxNjQ2MTNjYjJmNDFiNWFiZDliMGRmMjJmZjU2YjYpOwoKICAgICAgICAKICAgICAgICAKICAgICAgICAgICAgCiAgICAgICAgICAgICAgICB2YXIgZ2VvX2pzb25fODc3NWEwZjc5Njk1NDdiNDhkOGY5ZDg0MTQ5ZTQ5OGUgPSBMLmdlb0pzb24oCiAgICAgICAgICAgICAgICAgICAgeyJmZWF0dXJlcyI6IFt7Imdlb21ldHJ5IjogeyJjb29yZGluYXRlcyI6IFtbW1sxNDQuNjU4NDE5Njc3OTYwMjQsIC0xNS4yNjYxMTU2Nzk0MTIyXSwgWzE0NC42NTksIC0xNS4yNjYyXSwgWzE0NC4zNywgLTE2LjU4MTIzNjUyNjk0NjEwM10sIFsxNDQuMzI1LCAtMTYuNzg2XSwgWzE0NC4zMjQzNjM0NjI1ODUyLCAtMTYuNzg1OTA3NTgxNDQ4MjM1XSwgWzE0Mi42MSwgLTE2LjUzN10sIFsxNDIuNjI3MTYwNjYyMDE0MDMsIC0xNi40NjE5MjMzNDAwNDc1NDhdLCBbMTQyLjk1MzgzODgzODI3ODE2LCAtMTUuMDMyNzI5ODU0NzgzNjU3XSwgWzE0Mi45NTcsIC0xNS4wMTg5XSwgWzE0NC4wNjI2NzU5ODg2OTcxNSwgLTE1LjE3OTU1NDMzMTM3NzY3N10sIFsxNDQuMzcsIC0xNS4yMjQyMDg0MDE4ODAxMzldLCBbMTQ0LjYzNjkzMTQ1OTc2NTQ2LCAtMTUuMjYyOTkzNDQ4ODgzNjY2XSwgWzE0NC42NTg0MTk2Nzc5NjAyNCwgLTE1LjI2NjExNTY3OTQxMjJdXV0sIFtbWzE0NC4zMjQzNjM0NjI1ODUyLCAtMTYuNzg1OTA3NTgxNDQ4MjM1XSwgWzE0NC4wMjA2ODkyMjU1NzU4LCAtMTguMTU1MDU0ODY4NDI3NTczXSwgWzE0NC4wMDQsIC0xOC4yMzAzXSwgWzE0Mi4yNzcsIC0xNy45NzkzXSwgWzE0Mi4yOTQyODk2MjYyNDg2LCAtMTcuOTA0MzQyMDYwMzcyNDA3XSwgWzE0Mi42MjcsIC0xNi40NjE5XSwgWzE0Mi42MjcxNjA2NjIwMTQwMywgLTE2LjQ2MTkyMzM0MDA0NzU0OF0sIFsxNDQuMzQxLCAtMTYuNzEwOV0sIFsxNDQuMzI0MzYzNDYyNTg1MiwgLTE2Ljc4NTkwNzU4MTQ0ODIzNV1dXSwgW1tbMTQzLjY5NzgxNjgwNDQyNTgzLCAtMTkuNTk4ODczMzk5Mjg3MTA0XSwgWzE0My42ODIxNzU4Njc4MTEyLCAtMTkuNjY4NzQ2OTgzNDc1MjY2XSwgWzE0My42ODEsIC0xOS42NzRdLCBbMTQzLjY4MDk4NjIyNzY2MDMsIC0xOS42NzM5OTc5OTk3NjkyNl0sIFsxNDEuOTM5MTU3NjgzNjU0MjMsIC0xOS40MjEwMjI5MDEyNDI1NDhdLCBbMTQxLjkzOSwgLTE5LjQyMV0sIFsxNDIuMjk0LCAtMTcuOTA0M10sIFsxNDIuMjk0Mjg5NjI2MjQ4NiwgLTE3LjkwNDM0MjA2MDM3MjQwN10sIFsxNDQuMDIwNjg5MjI1NTc1OCwgLTE4LjE1NTA1NDg2ODQyNzU3M10sIFsxNDQuMDIxLCAtMTguMTU1MV0sIFsxNDMuNjk3ODE2ODA0NDI1ODMsIC0xOS41OTg4NzMzOTkyODcxMDRdXV0sIFtbWzE0MS42MDI0OTY2NTExMTkwMywgLTIwLjgzNDcxOTQ2NDE3OTEzOF0sIFsxNDEuOTM5MTU3NjgzNjU0MjMsIC0xOS40MjEwMjI5MDEyNDI1NDhdLCBbMTQxLjk1NywgLTE5LjM0NjFdLCBbMTQzLjY5NzgxNjgwNDQyNTgzLCAtMTkuNTk4ODczMzk5Mjg3MTA0XSwgWzE0My42OTgsIC0xOS41OTg5XSwgWzE0My42ODIxNzU4Njc4MTEyLCAtMTkuNjY4NzQ2OTgzNDc1MjY2XSwgWzE0My42ODA5ODYyMjc2NjAzLCAtMTkuNjczOTk3OTk5NzY5MjZdLCBbMTQzLjM1NCwgLTIxLjExNzNdLCBbMTQzLjM1Mzc4ODE4MzMwNjcsIC0yMS4xMTcyNjkyMzk1ODk0Ml0sIFsxNDEuNTk2LCAtMjAuODYyXSwgWzE0MS42MDI0OTY2NTExMTkwMywgLTIwLjgzNDcxOTQ2NDE3OTEzOF1dXSwgW1tbMTQzLjAyMjU5MDEzNTc4NTM1LCAtMjIuNTU5OTQwNDcxNTUyMzczXSwgWzE0MS4yNDg1Nzk1MzM5ODc1NywgLTIyLjMwMjI4NDE3MTE4OTg1N10sIFsxNDEuMjQ4LCAtMjIuMzAyMl0sIFsxNDEuNjAyNDk2NjUxMTE5MDMsIC0yMC44MzQ3MTk0NjQxNzkxMzhdLCBbMTQxLjYxNCwgLTIwLjc4NzFdLCBbMTQzLjM3MSwgLTIxLjA0MjJdLCBbMTQzLjM1Mzc4ODE4MzMwNjcsIC0yMS4xMTcyNjkyMzk1ODk0Ml0sIFsxNDMuMDIzLCAtMjIuNTZdLCBbMTQzLjAyMjU5MDEzNTc4NTM1LCAtMjIuNTU5OTQwNDcxNTUyMzczXV1dLCBbW1sxNDMuMDIyNTkwMTM1Nzg1MzUsIC0yMi41NTk5NDA0NzE1NTIzNzNdLCBbMTQyLjcwNTQyMDAwOTYwNzUsIC0yMy45MjcwMTU3OTk0OTg1MTNdLCBbMTQyLjY4OCwgLTI0LjAwMjFdLCBbMTQwLjg5NSwgLTIzLjc0MTZdLCBbMTQwLjkxMzM2MjMwMTg5MzksIC0yMy42NjY4NTI1OTc0MjQ5NzRdLCBbMTQxLjI0ODU3OTUzMzk4NzU3LCAtMjIuMzAyMjg0MTcxMTg5ODU3XSwgWzE0MS4yNjcsIC0yMi4yMjczXSwgWzE0My4wNCwgLTIyLjQ4NDldLCBbMTQzLjAyMjU5MDEzNTc4NTM1LCAtMjIuNTU5OTQwNDcxNTUyMzczXV1dLCBbW1sxNDIuNzA1NDIwMDA5NjA3NSwgLTIzLjkyNzAxNTc5OTQ5ODUxM10sIFsxNDIuNzA2LCAtMjMuOTI3MV0sIFsxNDIuMzY1NzE2MTQ0MDAwNywgLTI1LjM2ODQ1ODc2ODgxOTUyMl0sIFsxNDIuMzQ4LCAtMjUuNDQzNV0sIFsxNDAuNTM1LCAtMjUuMTgwMV0sIFsxNDAuNTM3MTU3MzcyNDMwNywgLTI1LjE3MTQ2MzA5MDc0MjM2XSwgWzE0MC45MTMsIC0yMy42NjY4XSwgWzE0MC45MTMzNjIzMDE4OTM5LCAtMjMuNjY2ODUyNTk3NDI0OTc0XSwgWzE0Mi43MDU0MjAwMDk2MDc1LCAtMjMuOTI3MDE1Nzk5NDk4NTEzXV1dLCBbW1sxNDIuMDIwOTQwMTMxNDgwNzYsIC0yNi44MDkyOTEzMDIyNDQwNDVdLCBbMTQyLjAxNjYwODYzNjUzMDY2LCAtMjYuODI3Mzc3MzgxODQ2OTQ0XSwgWzE0Mi4wMDMsIC0yNi44ODQyXSwgWzE0Mi4wMDI4MTE5MDcxNTM2LCAtMjYuODg0MTcyNjY4MDc4NzU4XSwgWzE0MC4xNjksIC0yNi42MTc3XSwgWzE0MC4xODgwMTUyMzkxMDg5MiwgLTI2LjU0MzAwMjIxMzk1MjM3N10sIFsxNDAuNTM3MTU3MzcyNDMwNywgLTI1LjE3MTQ2MzA5MDc0MjM2XSwgWzE0MC41NTQsIC0yNS4xMDUzXSwgWzE0Mi4zNjU3MTYxNDQwMDA3LCAtMjUuMzY4NDU4NzY4ODE5NTIyXSwgWzE0Mi4zNjYsIC0yNS4zNjg1XSwgWzE0Mi4wMjA5NDAxMzE0ODA3NiwgLTI2LjgwOTI5MTMwMjI0NDA0NV1dXSwgW1tbMTM5LjgxNDQwOTU5NTM4MTQsIC0yNy45Nzk2NTk0ODcyNzc1MjddLCBbMTQwLjE4OCwgLTI2LjU0M10sIFsxNDAuMTg4MDE1MjM5MTA4OTIsIC0yNi41NDMwMDIyMTM5NTIzNzddLCBbMTQyLjAyMDk0MDEzMTQ4MDc2LCAtMjYuODA5MjkxMzAyMjQ0MDQ1XSwgWzE0Mi4wMjEsIC0yNi44MDkzXSwgWzE0Mi4wMTY2MDg2MzY1MzA2NiwgLTI2LjgyNzM3NzM4MTg0Njk0NF0sIFsxNDIuMDAyODExOTA3MTUzNiwgLTI2Ljg4NDE3MjY2ODA3ODc1OF0sIFsxNDEuNjUzLCAtMjguMzI0Ml0sIFsxNDEuNjUyNDY3MzE2NjA1MTQsIC0yOC4zMjQxMjI2MjA0MjYxMl0sIFsxMzkuNzk1LCAtMjguMDU0M10sIFsxMzkuODE0NDA5NTk1MzgxNCwgLTI3Ljk3OTY1OTQ4NzI3NzUyN11dXSwgW1tbMTQxLjY1MjQ2NzMxNjYwNTE0LCAtMjguMzI0MTIyNjIwNDI2MTJdLCBbMTQxLjMxNDU0MzUwMDMwNjgzLCAtMjkuNjg4NDMzNzA4MDk0NTA3XSwgWzE0MS4yOTYsIC0yOS43NjMzXSwgWzEzOS40MTMsIC0yOS40ODk4XSwgWzEzOS40MjQ1Njg2Nzc2MTAzOCwgLTI5LjQ0NjIzMTM3OTIzMzk2OF0sIFsxMzkuODE0LCAtMjcuOTc5Nl0sIFsxMzkuODE0NDA5NTk1MzgxNCwgLTI3Ljk3OTY1OTQ4NzI3NzUyN10sIFsxNDEuNjcxLCAtMjguMjQ5M10sIFsxNDEuNjUyNDY3MzE2NjA1MTQsIC0yOC4zMjQxMjI2MjA0MjYxMl1dXSwgW1tbMTQxLjMxNDU0MzUwMDMwNjgzLCAtMjkuNjg4NDMzNzA4MDk0NTA3XSwgWzE0MS4zMTUsIC0yOS42ODg1XSwgWzE0MC45NTA5NzAxOTA4MDQzOCwgLTMxLjEyNjc1MDQxMDYzOTEwOF0sIFsxNDAuOTMyLCAtMzEuMjAxN10sIFsxMzkuMDIzLCAtMzAuOTI0Ml0sIFsxMzkuMDQzMjU4ODM2NzM4NSwgLTMwLjg0OTYzNzU5ODQ0Mjk0M10sIFsxMzkuNDI0NTY4Njc3NjEwMzgsIC0yOS40NDYyMzEzNzkyMzM5NjhdLCBbMTM5LjQzMywgLTI5LjQxNTJdLCBbMTQxLjMxNDU0MzUwMDMwNjgzLCAtMjkuNjg4NDMzNzA4MDk0NTA3XV1dLCBbW1sxNDAuOTUwOTcwMTkwODA0MzgsIC0zMS4xMjY3NTA0MTA2MzkxMDhdLCBbMTQwLjk1MiwgLTMxLjEyNjldLCBbMTQwLjU4MDM2NDQzNzE1OTQsIC0zMi41NjQyMDc2Njc4NDU0NjZdLCBbMTQwLjU2MSwgLTMyLjYzOTFdLCBbMTM4LjYyMywgLTMyLjM1NzRdLCBbMTM4LjYzNDYxODk3Nzc0NDY1LCAtMzIuMzE1Njg3ODY5ODk2NzJdLCBbMTM5LjA0MywgLTMwLjg0OTZdLCBbMTM5LjA0MzI1ODgzNjczODUsIC0zMC44NDk2Mzc1OTg0NDI5NDNdLCBbMTQwLjk1MDk3MDE5MDgwNDM4LCAtMzEuMTI2NzUwNDEwNjM5MTA4XV1dLCBbW1sxNDAuMjAxNzMzMTc1MTEyLCAtMzQuMDAwNzYxMjE3NDcwNTJdLCBbMTQwLjE4OTgzMjM4OTU4MDU3LCAtMzQuMDQ1ODM1MDY5ODM5MTY0XSwgWzE0MC4xODIsIC0zNC4wNzU1XSwgWzE0MC4xODE4MjQzOTE5MzYxNSwgLTM0LjA3NTQ3NDQ3NDg0NjE3NV0sIFsxMzguMjEzMDk2MDg4MjA1NDUsIC0zMy43ODkzMTM5NjY3MDYxNDZdLCBbMTM4LjIxMywgLTMzLjc4OTNdLCBbMTM4LjYzNDYxODk3Nzc0NDY1LCAtMzIuMzE1Njg3ODY5ODk2NzJdLCBbMTM4LjY0NCwgLTMyLjI4MjldLCBbMTQwLjU4MDM2NDQzNzE1OTQsIC0zMi41NjQyMDc2Njc4NDU0NjZdLCBbMTQwLjU4MSwgLTMyLjU2NDNdLCBbMTQwLjIwMTczMzE3NTExMiwgLTM0LjAwMDc2MTIxNzQ3MDUyXV1dLCBbW1sxMzcuODEwMzM1NjAzMDg3MjIsIC0zNS4xNTc2MDg2MTcwNTEyNzRdLCBbMTM4LjIxMzA5NjA4ODIwNTQ1LCAtMzMuNzg5MzEzOTY2NzA2MTQ2XSwgWzEzOC4yMzUsIC0zMy43MTQ5XSwgWzE0MC4yMDE3MzMxNzUxMTIsIC0zNC4wMDA3NjEyMTc0NzA1Ml0sIFsxNDAuMjAyLCAtMzQuMDAwOF0sIFsxNDAuMTg5ODMyMzg5NTgwNTcsIC0zNC4wNDU4MzUwNjk4MzkxNjRdLCBbMTQwLjE4MTgyNDM5MTkzNjE1LCAtMzQuMDc1NDc0NDc0ODQ2MTc1XSwgWzEzOS44MTQyMTMzODAyMDI1MiwgLTM1LjQzNjA4NTcyMTk1MTQ2XSwgWzEzOS43OTQsIC0zNS41MTA5XSwgWzEzNy43OTIsIC0zNS4yMTk5XSwgWzEzNy44MTAzMzU2MDMwODcyMiwgLTM1LjE1NzYwODYxNzA1MTI3NF1dXSwgW1tbMTM5LjQxNjcyNTA4MjM4NDcyLCAtMzYuODcwNTYwMDI1NDkyN10sIFsxMzkuNDA1MTI1NTE5NDcwMDYsIC0zNi45MTIzMzUwNjIzMzgxMTZdLCBbMTM5LjM5NiwgLTM2Ljk0NTJdLCBbMTM5LjM5NTc4MzcxOTE4MTcsIC0zNi45NDUxNjg1NTA2MjQyN10sIFsxMzcuMzU5LCAtMzYuNjQ5XSwgWzEzNy4zNjUzMzA3NjE4MTMwNiwgLTM2LjYyODA4MDY1ODQ5MjUzXSwgWzEzNy44MTAzMzU2MDMwODcyMiwgLTM1LjE1NzYwODYxNzA1MTI3NF0sIFsxMzcuODE0LCAtMzUuMTQ1NV0sIFsxMzkuODE0MjEzMzgwMjAyNTIsIC0zNS40MzYwODU3MjE5NTE0Nl0sIFsxMzkuODE1LCAtMzUuNDM2Ml0sIFsxMzkuNDE2NzI1MDgyMzg0NzIsIC0zNi44NzA1NjAwMjU0OTI3XV1dLCBbW1sxMzYuOTM2MTkyOTgzOTA0ODcsIC0zOC4wMDIzMjgwNTQzMTQwMzRdLCBbMTM3LjM2NTMzMDc2MTgxMzA2LCAtMzYuNjI4MDgwNjU4NDkyNTNdLCBbMTM3LjM4MiwgLTM2LjU3NDddLCBbMTM5LjQxNjcyNTA4MjM4NDcyLCAtMzYuODcwNTYwMDI1NDkyN10sIFsxMzkuNDE3LCAtMzYuODcwNl0sIFsxMzkuNDA1MTI1NTE5NDcwMDYsIC0zNi45MTIzMzUwNjIzMzgxMTZdLCBbMTM5LjM5NTc4MzcxOTE4MTcsIC0zNi45NDUxNjg1NTA2MjQyN10sIFsxMzkuMDA5MjU2MDAxNDEyOCwgLTM4LjMwMzY5MTg0Mzk4NTUyXSwgWzEzOC45ODgsIC0zOC4zNzg0XSwgWzEzNi45MTMsIC0zOC4wNzY2XSwgWzEzNi45MzYxOTI5ODM5MDQ4NywgLTM4LjAwMjMyODA1NDMxNDAzNF1dXSwgW1tbMTM4LjU4OTgzNTYyNzU5MDc2LCAtMzkuNzM1ODc2MDc5MjYzNjhdLCBbMTM4LjU4NDE5Nzc1MjgwOTQyLCAtMzkuNzU1MDkyMDQ4Mzk5NjddLCBbMTM4LjU2OCwgLTM5LjgxMDNdLCBbMTM4LjU2NzUyNzc0NTIwMTA2LCAtMzkuODEwMjMxMzA0MzM1MDA1XSwgWzEzNi40NTIyMTIxOTc2MDI5MywgLTM5LjUwMjUzMDg2NjkyOTJdLCBbMTM2LjQ1MiwgLTM5LjUwMjVdLCBbMTM2LjkzNiwgLTM4LjAwMjNdLCBbMTM2LjkzNjE5Mjk4MzkwNDg3LCAtMzguMDAyMzI4MDU0MzE0MDM0XSwgWzEzOS4wMDkyNTYwMDE0MTI4LCAtMzguMzAzNjkxODQzOTg1NTJdLCBbMTM5LjAxLCAtMzguMzAzOF0sIFsxMzguNTg5ODM1NjI3NTkwNzYsIC0zOS43MzU4NzYwNzkyNjM2OF1dXV0sICJ0eXBlIjogIk11bHRpUG9seWdvbiJ9LCAicHJvcGVydGllcyI6IHsic3R5bGUiOiB7ImZpbGxDb2xvciI6ICIjMDBmZjAwIn19LCAidHlwZSI6ICJGZWF0dXJlIn1dLCAidHlwZSI6ICJGZWF0dXJlQ29sbGVjdGlvbiJ9KQogICAgICAgICAgICAgICAgICAgIC5hZGRUbyhtYXBfNDAxNjQ2MTNjYjJmNDFiNWFiZDliMGRmMjJmZjU2YjYpOwogICAgICAgICAgICAgICAgZ2VvX2pzb25fODc3NWEwZjc5Njk1NDdiNDhkOGY5ZDg0MTQ5ZTQ5OGUuc2V0U3R5bGUoZnVuY3Rpb24oZmVhdHVyZSkge3JldHVybiBmZWF0dXJlLnByb3BlcnRpZXMuc3R5bGU7fSk7CiAgICAgICAgICAgIAogICAgICAgIAogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHZhciBnZW9fanNvbl82NzkxNTIzYzcxYWM0YTVjODYyZjcyMTA1MmQwM2RmYyA9IEwuZ2VvSnNvbigKICAgICAgICAgICAgICAgICAgICB7ImZlYXR1cmVzIjogW3siZ2VvbWV0cnkiOiB7ImNvb3JkaW5hdGVzIjogW1sxNDMuNTE2MDgzMzMzMzMzMzQsIC0xNi4zOTYzMzMzMzMzMzMzMzVdLCBbMTQzLjQ0ODQ3MjIyMjIyMjIyLCAtMTYuNjk2MjIyMjIyMjIyMjJdLCBbMTQzLjM4MDcyMjIyMjIyMjIzLCAtMTYuOTk2MDgzMzMzMzMzMzM1XSwgWzE0My4zMTI4MzMzMzMzMzMzNCwgLTE3LjI5NTkxNjY2NjY2NjY2N10sIFsxNDMuMjQ0Nzc3Nzc3Nzc3NzYsIC0xNy41OTU3NV0sIFsxNDMuMTc2NTU1NTU1NTU1NTUsIC0xNy44OTU1ODMzMzMzMzMzMzVdLCBbMTQzLjEwODE5NDQ0NDQ0NDQ1LCAtMTguMTk1MzYxMTExMTExMTFdLCBbMTQzLjAzOTYzODg4ODg4ODg3LCAtMTguNDk1MTY2NjY2NjY2NjY2XSwgWzE0Mi45NzA5NDQ0NDQ0NDQ0NiwgLTE4Ljc5NDkxNjY2NjY2NjY3XSwgWzE0Mi45MDIwODMzMzMzMzMzNCwgLTE5LjA5NDY5NDQ0NDQ0NDQ0M10sIFsxNDIuODMzMDI3Nzc3Nzc3NzcsIC0xOS4zOTQ0MTY2NjY2NjY2NjVdLCBbMTQyLjc2MzgzMzMzMzMzMzM0LCAtMTkuNjk0MTM4ODg4ODg4ODldLCBbMTQyLjY5NDQxNjY2NjY2NjY4LCAtMTkuOTkzODYxMTExMTExMTEyXSwgWzE0Mi42MjQ4MzMzMzMzMzMzMywgLTIwLjI5MzUyNzc3Nzc3Nzc4XSwgWzE0Mi41NTUwNTU1NTU1NTU1NSwgLTIwLjU5MzE5NDQ0NDQ0NDQ0M10sIFsxNDIuNDg1MDgzMzMzMzMzMywgLTIwLjg5Mjg2MTExMTExMTExXSwgWzE0Mi40MTQ5MTY2NjY2NjY2NywgLTIxLjE5MjVdLCBbMTQyLjM0NDU4MzMzMzMzMzMzLCAtMjEuNDkyMTExMTExMTExMTFdLCBbMTQyLjI3NCwgLTIxLjc5MTcyMjIyMjIyMjIyM10sIFsxNDIuMjAzMjIyMjIyMjIyMjIsIC0yMi4wOTEzMDU1NTU1NTU1NTRdLCBbMTQyLjEzMjIyMjIyMjIyMjIzLCAtMjIuMzkwODYxMTExMTExMTFdLCBbMTQyLjA2MTAyNzc3Nzc3Nzc4LCAtMjIuNjkwMzg4ODg4ODg4ODldLCBbMTQxLjk4OTYxMTExMTExMTEsIC0yMi45ODk5MTY2NjY2NjY2NjZdLCBbMTQxLjkxNzk3MjIyMjIyMjIsIC0yMy4yODk0MTY2NjY2NjY2NjhdLCBbMTQxLjg0NjA4MzMzMzMzMzM1LCAtMjMuNTg4OTE2NjY2NjY2NjY2XSwgWzE0MS43NzM5NzIyMjIyMjIyMywgLTIzLjg4ODM2MTExMTExMTExXSwgWzE0MS43MDE2Mzg4ODg4ODg4OCwgLTI0LjE4NzgwNTU1NTU1NTU1Nl0sIFsxNDEuNjI5MDU1NTU1NTU1NTcsIC0yNC40ODcyMjIyMjIyMjIyMjJdLCBbMTQxLjU1NjIyMjIyMjIyMjIzLCAtMjQuNzg2NjM4ODg4ODg4ODldLCBbMTQxLjQ4MzEzODg4ODg4ODksIC0yNS4wODZdLCBbMTQxLjQwOTc3Nzc3Nzc3Nzc4LCAtMjUuMzg1MzYxMTExMTExMTFdLCBbMTQxLjMzNjE5NDQ0NDQ0NDQ2LCAtMjUuNjg0Njk0NDQ0NDQ0NDQ2XSwgWzE0MS4yNjIzMzMzMzMzMzMzNCwgLTI1Ljk4NF0sIFsxNDEuMTg4MTk0NDQ0NDQ0NDMsIC0yNi4yODMyNzc3Nzc3Nzc3NzZdLCBbMTQxLjExMzc3Nzc3Nzc3Nzc4LCAtMjYuNTgyNTU1NTU1NTU1NTU1XSwgWzE0MS4wMzkwODMzMzMzMzMzNCwgLTI2Ljg4MTgwNTU1NTU1NTU1NV0sIFsxNDAuOTY0MTExMTExMTExMSwgLTI3LjE4MV0sIFsxNDAuODg4ODMzMzMzMzMzMzQsIC0yNy40ODAxOTQ0NDQ0NDQ0NDNdLCBbMTQwLjgxMzI3Nzc3Nzc3Nzc4LCAtMjcuNzc5MzYxMTExMTExMTFdLCBbMTQwLjczNzQxNjY2NjY2NjY2LCAtMjguMDc4NV0sIFsxNDAuNjYxMjUsIC0yOC4zNzc2MTExMTExMTExMV0sIFsxNDAuNTg0NzUsIC0yOC42NzY3MjIyMjIyMjIyMjRdLCBbMTQwLjUwNzk3MjIyMjIyMjI0LCAtMjguOTc1Nzc3Nzc3Nzc3Nzc1XSwgWzE0MC40MzA4MzMzMzMzMzMzLCAtMjkuMjc0ODA1NTU1NTU1NTU2XSwgWzE0MC4zNTMzODg4ODg4ODg4NywgLTI5LjU3MzgzMzMzMzMzMzMzM10sIFsxNDAuMjc1NjExMTExMTExMTIsIC0yOS44NzI4MDU1NTU1NTU1NTVdLCBbMTQwLjE5NzQ3MjIyMjIyMjIyLCAtMzAuMTcxNzUwMDAwMDAwMDAzXSwgWzE0MC4xMTg5NzIyMjIyMjIyMywgLTMwLjQ3MDY2NjY2NjY2NjY2Nl0sIFsxNDAuMDQwMTM4ODg4ODg4OSwgLTMwLjc2OTU4MzMzMzMzMzMzM10sIFsxMzkuOTYwOTQ0NDQ0NDQ0NDQsIC0zMS4wNjg0NDQ0NDQ0NDQ0NDVdLCBbMTM5Ljg4MTM4ODg4ODg4ODksIC0zMS4zNjcyNzc3Nzc3Nzc3OF0sIFsxMzkuODAxNDcyMjIyMjIyMjMsIC0zMS42NjYwODMzMzMzMzMzMzNdLCBbMTM5LjcyMTEzODg4ODg4ODksIC0zMS45NjQ4NjExMTExMTExMV0sIFsxMzkuNjQwNDQ0NDQ0NDQ0NDMsIC0zMi4yNjM2MTExMTExMTExMV0sIFsxMzkuNTU5MzMzMzMzMzMzMzQsIC0zMi41NjIzMDU1NTU1NTU1NTRdLCBbMTM5LjQ3NzgzMzMzMzMzMzM0LCAtMzIuODYwOTcyMjIyMjIyMjJdLCBbMTM5LjM5NTkxNjY2NjY2NjY2LCAtMzMuMTU5NjM4ODg4ODg4ODg1XSwgWzEzOS4zMTM2MTExMTExMTExMywgLTMzLjQ1ODI1XSwgWzEzOS4yMzA4NjExMTExMTExMiwgLTMzLjc1NjgzMzMzMzMzMzMzXSwgWzEzOS4xNDc2NjY2NjY2NjY2NSwgLTM0LjA1NTM2MTExMTExMTExXSwgWzEzOS4wNjQwNTU1NTU1NTU1NywgLTM0LjM1Mzg2MTExMTExMTExNV0sIFsxMzguOTc5OTcyMjIyMjIyMjIsIC0zNC42NTIzNjExMTExMTExMV0sIFsxMzguODk1NDQ0NDQ0NDQ0NDUsIC0zNC45NTA3Nzc3Nzc3Nzc3OF0sIFsxMzguODEwNDQ0NDQ0NDQ0NDQsIC0zNS4yNDkxOTQ0NDQ0NDQ0NV0sIFsxMzguNzI0OTcyMjIyMjIyMjIsIC0zNS41NDc1NTU1NTU1NTU1NTRdLCBbMTM4LjYzOTA1NTU1NTU1NTU2LCAtMzUuODQ1ODYxMTExMTExMTFdLCBbMTM4LjU1MjYxMTExMTExMTEzLCAtMzYuMTQ0MTY2NjY2NjY2NjZdLCBbMTM4LjQ2NTY2NjY2NjY2NjY2LCAtMzYuNDQyMzg4ODg4ODg4ODg1XSwgWzEzOC4zNzgyMjIyMjIyMjIyMywgLTM2Ljc0MDYxMTExMTExMTExNF0sIFsxMzguMjkwMjUsIC0zNy4wMzg3Nzc3Nzc3Nzc3NzRdLCBbMTM4LjIwMTc0OTk5OTk5OTk4LCAtMzcuMzM2ODg4ODg4ODg4ODldLCBbMTM4LjExMjcyMjIyMjIyMjIsIC0zNy42MzQ5NzIyMjIyMjIyMjRdLCBbMTM4LjAyMzEzODg4ODg4ODksIC0zNy45MzMwMjc3Nzc3Nzc3NzRdLCBbMTM3LjkzMywgLTM4LjIzMV0sIFsxMzcuODQyMjc3Nzc3Nzc3NzgsIC0zOC41Mjg5NzIyMjIyMjIyMl0sIFsxMzcuNzUxLCAtMzguODI2ODYxMTExMTExMTE0XSwgWzEzNy43MTQzMDU1NTU1NTU1NCwgLTM4Ljk0NjAyNzc3Nzc3Nzc3XV0sICJ0eXBlIjogIkxpbmVTdHJpbmcifSwgInByb3BlcnRpZXMiOiB7InN0eWxlIjogeyJmaWxsQ29sb3IiOiAiIzAwZmYwMCJ9fSwgInR5cGUiOiAiRmVhdHVyZSJ9XSwgInR5cGUiOiAiRmVhdHVyZUNvbGxlY3Rpb24ifSkKICAgICAgICAgICAgICAgICAgICAuYWRkVG8obWFwXzQwMTY0NjEzY2IyZjQxYjVhYmQ5YjBkZjIyZmY1NmI2KTsKICAgICAgICAgICAgICAgIGdlb19qc29uXzY3OTE1MjNjNzFhYzRhNWM4NjJmNzIxMDUyZDAzZGZjLnNldFN0eWxlKGZ1bmN0aW9uKGZlYXR1cmUpIHtyZXR1cm4gZmVhdHVyZS5wcm9wZXJ0aWVzLnN0eWxlO30pOwogICAgICAgICAgICAKICAgICAgICAKICAgICAgICAKICAgICAgICA8L3NjcmlwdD4KICAgICAgICA=\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;\"></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x7f69c692fb10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "scene = mapping(satellite_scenes)\n",
    "track = mapping(satellite_track)\n",
    "style_function = lambda x: {'fillColor': '#000000' if x['type'] == 'Polygon' else '#00ff00'}\n",
    "mapa = folium.Map(location=[-30,150], zoom_start=4)\n",
    "poly = folium.features.GeoJson(scene, style_function=style_function)\n",
    "line = folium.features.GeoJson(track, style_function=style_function)\n",
    "mapa.add_children(poly)\n",
    "mapa.add_children(line)\n",
    "display(mapa)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Save as shapefiles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def save_shapefile(shape, filename, vector_type):\n",
    "    driver = ogr.GetDriverByName('Esri Shapefile')\n",
    "    ds = driver.CreateDataSource(filename)\n",
    "    if vector_type is 'polygon':\n",
    "        ogr_type = ogr.wkbLineString\n",
    "    if vector_type is 'line':\n",
    "        ogr_type = ogr.wkbPolygon\n",
    "    layer = ds.CreateLayer('', None, ogr.wkbLineString)\n",
    "\n",
    "    layer.CreateField(ogr.FieldDefn('id', ogr.OFTInteger))\n",
    "    defn = layer.GetLayerDefn()\n",
    "\n",
    "    feat = ogr.Feature(defn)\n",
    "    feat.SetField('id', 123)\n",
    "\n",
    "    geom = ogr.CreateGeometryFromWkb(shape.wkb)\n",
    "    feat.SetGeometry(geom)\n",
    "\n",
    "    layer.CreateFeature(feat)\n",
    "    feat = geom = None\n",
    "\n",
    "    ds = layer = feat = geom = None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "save_shapefile(satellite_track, 'satellite_track.shp', 'line')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "save_shapefile(satellite_scenes, 'satellite_scenes.shp', 'polygon')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.9"
  },
  "widgets": {
   "state": {},
   "version": "1.1.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}