{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEYCAYAAABRB/GsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deXxU1fn48c+TjexsgUASIAHZVyUgiEjAFTeq1q3VWq2lq1pb22LV6rfubdXaXdqKtlWxLi3+BBVcUEGUfVdZA4R9JwlkP78/7p1hkkySCZmZOzP3eb9e82Lm3jv3Pjkk88xZ7jlijEEppZR7xTkdgFJKKWdpIlBKKZfTRKCUUi6niUAppVxOE4FSSrmcJgKllHI5TQRKKeVymgiUUsrlNBEoVxGRYhE5z+k4lIokmghUq8XKh6mTP4eInC0iR0VEGmx/T0R+3orzLBaRfiLSW0SWN9j3bxHZLSLHRGSDiNzayv0DReR9O85NInJFIPsCObeKLJoIlHLG6cBK03iOlxHAcj/HNyIiiUAvYCMw0s/7HgXyjTGZwOXAQyIyMpD9IpIAzALeBDoBU4F/20mnyX2tuLaKIJoIVNDY3xLni8gREVknIpcHut/+dn63iKwXkcMiMkNEkv1c4+ci8mqDbU+LyO8DicHnPf8CegL/T0TKRORn9vZpIrJZRErtWHy/BZ8hIivsfa+IyMsi8pDP/hwReU1E9ovIVhG5vZniOh1Y0SCmfKwP1hV+jvdnCLDeTiaFNEgExph1xphKz0v70SfA/QOAHOApY0ytMeZ9YCFwYwv7Arq2ijDGGH3oo1UPoBg4r8G2RGAT8AsgCZgElAL9A9xfDKwFemB9GC4EHvJz7V7AcSDDfh0P7AbGtHSNhrE38XNcjfUhFwdcC5QD3e3zbQPusK9zJVDlidE+fhnwS/vY3sAW4MImynA58I0G264AdgRQ/jcDR+xyqLCf19g/6xGgwOfYP9vHGfua6Q3O5Xc/VpIpA8Tn2HnAf5vbF8i59RF5D8cD0Ef0PZr4AB0P7AHifLa9BDwQ4P5i4Ls++y4GNjdx/QWeD1HgfM9xLV2jYez+fg4/11oJTAHOAXY2+PBb4JMIzgS2N3jv3cAMP+dMBCqBoQ22PwjM8nn9FHB6M7F9jNWU1NOOU5o4Lh44G7gXSAxkvx3jFuBn9vMLsBLfO83ta+219REZD20aUsGSg/Vtts5n2zYgN8D9ADsa7Mtp4lovAtfbz79mvw70Gs0SkW+IyEq7aekI1rffLPvcO4396eYn3l5Ajud99nt/AWT7ucxA+9/PG2wfDSxpcNz6BvF1ss9/FDgLmA98CfQHDovIjxpezFjNNwuAPOB7gew3xlQDXwEuwUquPwH+A5Q0t6+111aRIcHpAFTM2AX0EJE4nw/insCGAPeD1SyEz75dTVzrFeAJEcnDak4Z24pr+KrXUSsivYC/AecCi4wxtSKyEhCs5qdcERGfZNAD2Gw/3wFsNcb0beJavrKBE8aYGp9rZwETgJ/6HJdsTrazWwEbcwjoICLXARONMd8Rkf8CfzLGvNvCdRNovp2+3n5jzGo7Jk+MnwDPt7TvFK+tHKQ1AnWqEkUk2fMAPsNqD/6ZiCSKSBFwGTDTPr6l/QA/EJE8EekE3AO87O/Cxpj9WN+EZ2B9+Hq+WQdyDV97sdryPdKwksN+ABG5GatGALAIqAV+KCIJIjIF6xu8x2Kg1O7MThGReBEZIiKj/Fx3DdBORL5tH9sPqwnrTfsDFhHpasfXFN9RQqdj9U94iUhXEblORNLtWC7EqkW9F8h++5hh9v9vqojchdVX8lwA+1o8t4owTrdN6SP6Hlht66bB4yFgMPAhcBSrSeOKBu9rcr99zrvt7Uewvl2mNhPDjfZ1fxroNXyu4+kjmAJst693l73tYeAQcAB40j7Xrfa+Qqy2+DKsWsnrwH0+587B+kDfAxwGPqWJPgisBLUGOIHV3v44Pp2pWLWSe5v5+d8HRgGd8dOXAnSxYz8CHLOv9e1A99vH/Mb+OcqAt4DTAtzX4rn1EVkPsf/jlHKUiBRjfeC21LwRMUTkM+CvxpgZITj3j7BqO7OCfW6lGtKmIaUCJCITRKSb3TR0EzAMeDtElxuK9U1aqZDTzmKlAtcfa3RMGlZzzleNMbtDcSFjzLdCcV6l/NGmIaWUcjltGlJKKZeLyqahrKwsk5+f3+Jx5eXlpKWlhT6gKKPl0piWiX9aLv5Fa7ksW7bsgDGmS8PtUZkI8vPzWbp0aYvHzZ8/n6KiotAHFGW0XBrTMvFPy8W/aC0XEdnmb7s2DSmllMtpIlBKKZfTRKCUUi6niUAppVxOE4FSSrmcJgKllHI5TQRKKeVymgiUcjFjDL+ctZa56/Y4HYpyUNgSgYj0EJEPRGS9iKwTkTvs7Q+IyE57ecCVInJxuGJSyu2mvbaGfy7axtR/LaOyptbpcJRDwnlncQ3wE2PMchHJAJaJyDx731PGmN+GMRalFPDy0pPLLq/YfoQxvTs7GI1ySthqBMaY3caY5fbzUqyFuwNeVFwpFVpf7il1OgTlEEemoRaRfOAjrPVgfwx8E2tJu6VYtYbDft4zFZgKkJ2dPXLmzKaWoT2prKyM9PT0YIUdM7RcGnNjmRhjuPmd45yTl8Di3TWcnZvADYPa1TvGjeUSiGgtl4kTJy4zxhQ23B72SedEJB14DfiRMeaYiPwFeBBr/dkHgSeAWxq+zxgzHZgOUFhYaAKZ8ClaJ4YKNS2XxtxWJvtLK/n7x1uALZw97DQOmV1UJbejqGh0vePcVi6BirVyCWsiEJFErCTwgjHmdQBjzF6f/X8D3gxnTEq50c9fW837X+wDoGNqEgVZ6azc0agirlwinKOGBPgH8Lkx5kmf7d19DrsCWBuumJRyq/LKGu9zAxRkpVFy+ISOHHKpcN5HMA64EZjUYKjor0VkjYisBiYCd4YxJqVcSeTk86+MyKV3VhrGwPaDx50LSjkmbE1DxpgFgPjZNSdcMSilGktKiKMgy1pta+uBcvpmZzgckQo3vbNYKUW+TyJQ7qOJQCkXqqqpq/e6fUoindOSNBG4lCYCpVyovLJxp3BBVhpbNBG4kiYCpVwoJSkegCkjcrzbCrLStEbgUpoIlHKhs/pYcwr9+qvDvNsKuqSxv7SS0opqp8JSDtFEoJQLGSApPo52CfHebT06pgKw88gJh6JSTtFEoJQL1RnTaDB3TocUAHZpInAdTQRKuZGBuAaJINdOBDsPayJwG00ESrlQnTFIgypB14x2JMXHUaKJwHU0ESjlQsZPjSAuTsjrmMKOwzrNhNtoIlDKheoMiDSe8SW3Y4o2DbmQJgKlXMhg8JMHyOuYok1DLqSJQCkX8jNoCIC8jqkcLK/iRJVOR+0mmgiUciFjDHENOwnwGTl0RPsJ3EQTgVIuVNdEjSC3o5UItHnIXTQRKOVCBkOcn06CvI6eGoEmAjfRRKCUC1mjhhpv75qRTEKcaI3AZTQRKOVCponho/FxQk4HHULqNpoIlHIhY4zfPgKwOoxL9KYyV9FEoJQLWXcW+08Fei+B+2giUMqF6oz/G8oAenRKZV9ppd5L4CKaCJRyIUPTNYKenax1CbR5yD00ESjlQnXGNLmvRydrCKlOPucemgiUciMDcU389fewawQ7Dmk/gVtoIlDKhfytR+DRJb0dyYlx7DikNQK30ESglAtZfQT+94kIeR1TtWnIRTQRKOVCTa1H4NGjYwrbtWnINTQRKOVCppnho2D1E5QcOo5pplNZxY6wJQIR6SEiH4jIehFZJyJ32Ns7icg8Edlo/9sxXDEp5VZNrUfg0bNTKqWVNZRXhy0k5aBw1ghqgJ8YYwYBY4AfiMggYBrwnjGmL/Ce/VopFUJNzT7qkdfRGjl04ERduEJSDgpbIjDG7DbGLLeflwKfA7nAFOB5+7Dnga+EKyal3Kquzv/sox6em8r2HdemITdwpI9ARPKB04HPgGxjzG571x4g24mYlHKTlmoEPTvbiUBrBK6QEO4Likg68BrwI2PMMd+RC8YYIyJ+v4KIyFRgKkB2djbz589v8VplZWUBHec2Wi6Nua1M9u2voPyEafZnzkiCXUerXFUugYq135ewJgIRScRKAi8YY163N+8Vke7GmN0i0h3Y5++9xpjpwHSAwsJCU1RU1OL15s+fTyDHuY2WS2NuK5N/b1tK5ZETFBWNb/KYPusWcvT4MVeVS6Bi7fclnKOGBPgH8Lkx5kmfXW8AN9nPbwJmhSsmpdyqpeGjYPUTaB+BO4Szj2AccCMwSURW2o+LgceA80VkI3Ce/VopFULNzT7q0bNTKgcrDNW12k8Q68LWNGSMWUDTQ5fPDVccSqnm1yPwKMhKo87AjkPH6d0lPTyBKUfoncVKuVBTaxb7ys9KA6D4YHk4QlIO0kSglAvVNbNmsUe+PYR020GdfC7WaSJQyqWamn3Uo1NaEsnxmgjcQBOBUi5k9RE0nwlEhOy0OG0acgFNBEq5kDEt1wgAuqaK1ghcQBOBUi7U3AplvrJTrZXKKmtqwxCVcoomAqVcyFqYpuXj8jLiqKkzbNmvzUOxTBOBUm5kWr6hDKB7mnXM1gOaCGKZJgKlXCiQG8rAahoC2LK/LMQRKSdpIlDKhQKZYgIgOUHolpnMFq0RxDRNBEq5UKA1AoDeXdK0aSjGaSJQyoUCmWLCoyArjS37y3Uh+ximiUApFzIBTDHh0btLOkdPVHOwvCqkMSnnaCJQyoWsPoLAju3b1Zp5dONe7TCOVZoIlHKhQKaY8OibbSWCTTpyKGZpIlDKhQKdYgKgW2YyaUnxbN6niSBWaSJQyoXqDDS9TlR9IkKfruls0kQQszQRKOVCxpiAawQAfbtmsGFvaegCUo7SRKCUC5kA5xry6N8tnX2llRzWkUMxSROBUi5kMAHdWezRLzsDQGsFMUoTgVIuFOjsox79u2kiiGWaCJRyIdOK4aNgjRzKSE5gg95LEJM0ESjlQsYEOmbIIiKc1jWdjfu0RhCLNBEo5UKBzj7qq1/XDL7cU6pzDsUgTQRKuVBrZh/1GNA9g8PHq9lfWhmaoJRjNBEo5UImwBXKfA3sngnA+t3HQhGScpAmAqVcqK4Vs496DOymiSBWaSJQyoVasx6BR/vURHp0SmHdTk0EsUYTgVIuZE6hjwBgWG4HVu88EvyAlKPClghE5FkR2Scia322PSAiO0Vkpf24OFzxKOVmrVmPwNfQvPbsOHSCI8d1qolYEs4awXPARX62P2WMGWE/5oQxHqVcy+ojaH0mGJrbHoA1O48GOyTloLAlAmPMR8ChcF1PKdU0YyDuFP76h+RYiWB1iSaCWJLgdADAD0XkG8BS4CfGmMP+DhKRqcBUgOzsbObPn9/iicvKygI6zm20XBpzW5lUVFaxe9du5s9v/ruZv3LJThXeX7mJwVISwggjW6z9vjidCP4CPIjVZPkg8ARwi78DjTHTgekAhYWFpqioqMWTz58/n0COcxstl8bcViZJC+aRm9uNoqKhzR7nr1xG717B8m2HXVVeDcXa74ujo4aMMXuNMbXGmDrgb8BoJ+NRyi1aO/uor2G57dl55AQHyvQO41jhaCIQke4+L68A1jZ1rFIqeKwVyk4tEwzN0w7jWBO2piEReQkoArJEpAS4HygSkRFYTUPFwHfCFY9SblbXytlHfQ3OyUQE1pQcZWL/rkGNSzkjbInAGHO9n83/CNf1lVIntXY9Al8ZyYn0zkrTGkEM0TuLlXKh1q5Z3NDQ3Pas0SGkMUMTgVIudCrrEfgamteBPccq2HesInhBKcdoIlDKhU5l9lFfw7TDOKa02EcgIp0COE+dMUZnolIqSlh3Fp96KhjUPZM4se4wPndgdhAjU04IpLN4l/1o7rcmHugZlIiUUiHX1hpBWrsETuuarjWCGBFIIvjcGHN6cweIyIogxaOUCgND69cjaGhobgc+3LC/TSOQVGQIpI9gbJCOUUpFiFNdj8DXsLz2HCirZI92GEe9QBLBEyJydnMHGGP0N0GpKGKtWdy2c3juMF61Q5uHol0giWAD8BsRKRaRX4tIs81ESqnId6rrEfga1D2ThDhh5Q4dJxLtWkwExpinjTFjgQnAQeBZEflCRO4XkX4hj1ApxbpdRzn78feDNtHbqa5Q5is5MZ7BOZks3+535ngVRQK+j8AYs80Y87jdcXw98BXg85BFppTyuuT3Cyg5fIKPN+5v87mMMRirt7jN5yrM78SqHUeorKlt87mUcwJOBCKSICKXicgLwFvAl8CVIYtMKdVIW+4G9jDGc642n4rRBZ2orKnTfoIo12IiEJHzReRZoAT4NjAb6GOMuc4YMyvUASqlTqqtM20+h+cMwUgqo/Ot+02XFOsqtNEskPsI7gZepJllJJVS4VFeWdPmc9TZVYJgjPzvmJZEny5pLNumHw3RrMVEYIyZFI5AlFItK69qe1u8t2koGG1DWM1Db67aTW2dIT5I51ThFUjT0PJgHKOUOnVpSfEAHA9ijSBYzizoTGllDZ/vPhbU86rwCaRpaKCIrG5mvwDtgxSPUsqPGrtvIBg1Ao9g9BEAFOZ3BOC5T4r57dXDg3JOFV6BJIIBARyjY8eUCqHq2joA/rFgK/ddOqhN5/L2EQSpFSevYyoAu4+eCM4JVdgF0kewLRyBKKX8q6szBGGwkFcwh496fO3MnryxchfVtXUkxusyJ9FG/8eUinBvr9vjfZ5q9xW0xclRQ8HLBOf0zaKssoYV23W6iWikiUCpCHfbSydnea+pNZg2dvZ63h3MmaPH9skiPk6CcuezCr9WJwIRmSEifxCRm0VkiOhE5EqF1Pi+Wd7nVbV1HG9jh7GxuhuCuoZA+5RERvTowEcbDwTtnCp8Ahk+muL72hhzM/AzrFlJzwOeCU1oSimA7Izkeq+/OWNxm85n7DpBsIf8j++bxeqSIxwurwruiVXIBVIj2CUiX4rIKyJyn4hMAZKNMQuNMb8zxkwNdZBKuVlKUjyJ8cIDl1mjhZYUt+0uXk/Hc7Cr8uf064IxsHCz1gqiTSDTUHcELgL+bW+6AVgnIi+JiN4/oFSIVdbU0jE1iRvG9ArK+Yx3+GhwU8Gw3PZkJifw0QbtJ4g2AfURGGO2GmNmGWMeNMZcDeQC64DfhTQ6pRSVNdaQzIT4OKae0xuA66d/esrnC0VnMUBCfBzjTsvi440H2tyhrcLrlEYNGctDwJggx6OUaqC61tAuwfpT9bS/L9pykE+3HDyl8xU+9C4A/1wU/FuExvftwu6jFSzXYaRRJZDO4h+LyHki0rXB9nZAchNvU0oFSVVNrfcmrcuG53i3X2fXCo6eqKaiOrCRRL6zl27aVxbEKC3jTusMwLJtOi11NAmkRpAN/BhYJiK7RWSuiPwFWAi8FtLolHKZFdsPs2V//Q/o6lpDkl0jOKdfF24/t693X12dYfj/zWXAfW8HdP5Xl5V4n4/o0SEIEdfXq3MafbumM2/93qCfW4VOIJ3FPzfGXGyM6QEMBB4ElgF3G2PuCvRCIvKsiOwTkbU+2zqJyDwR2Wj/2/FUfgilYsUVf/6ESU98SGlFtXdbVU0difEnG/TvPO9kIuj9izne52UBzEx6/xvrghRp0y4fnsOS4sPsLw3O+soq9FrVR2CMOWKM+dgY83djzLxWXus5rNFHvqYB7xlj+gLv2a+Vcp3aOsPeYxXe10MfmMv10z8lf9psFmw6UK/NXUS4aWzjEURrSppfLnLqP5fWex2qDt3x/boA8IkOI40aYZtiwhjzEdCw4XAK8Lz9/HngK+GKR6lIcuYj73LmI+/V27aomc7gqRP6eJ977i94+r0NzV5jrt1cc+HgbAC+fmZwhqM2NDS3PVnpSd7rqcgn4RzmJSL5wJvGmCH26yPGmA72cwEOe177ee9UYCpAdnb2yJkzZ7Z4vbKyMtLT04MTfAzRcmks3GUyr7iaF76wRgD9eGQ7nlzWcjPKcxel1XtdXm1IjIOkeOGbb5eTmQS/n5TW6H3Hqw3rDtbyp5XWNW4alMTEnokBxXmq5TJjbSWf7a7hD+emkhiDq5ZF69/QxIkTlxljChtuD2Q9grAwxhgRaTIrGWOmA9MBCgsLTVFRUYvnnD9/PoEc5zZaLo0Fq0xWbD9MSlI8A7plNn2tL/fxwttLvK8DSQJAs/HdWraef366jXHjz2k0DXT+tNkAxMcJp/fowAM3jg34ZrJTLRfTfR8fzlhCQs5gigZ0bfkNUSbW/oacnn10r4h0B7D/3edwPBHjWEU1t720IiiLlavw2HqgnCv+/AkX/e7jZo/75owlze7f+ujFvP2j8fW2jendqdn3DM1rT1VNHd9/YXm9tn/ffoPaOsPgnMyg31Hsz1l9rGGk3/n3spBfS7Wd0zWCN4CbgMfsf2c5G05kWLT5INf/zRoj/v9W7aL4sUscjkg1paK6lkPlVeR0SOGNlbsa7a+sqeWxt77guxP6kJ3Z8m03Gx6ajIgwoFsm9182iJ6dUjl6oprLfe4f8OeMntaAu3nr91Jw9xzuu3QQD765vtFx7VMCaxJqq3YJ1roJVTV1VFTXkpzY9nUUVOiErUYgIi8Bi4D+IlIiIt/CSgDni8hGrJlMHwtXPJHKGONNAh5rdzY/GkQ548MN+7np2cWc9dj71NTW8dS7Jztr+9/7Fieqanl77R5mLCzmd+9uBGDj3lLvMb/0WXLylnEFLLp7kvd+AYCbxxVw7sBsrjwjj4QWVv3q0Sm13mt/SQCgfWpS4D9gG/3zltGAtZaximxhqxEYY65vYte54YohGjw5z/owOXdAV568dgSFD83j1WUlDMnV+f0iyb3/W8O/P93ufX3aPW/V219ZU8cN//iMZdusmUJfWrydBy4fxPlPfeQ95pazC8jpkMKJ6hquOD2vzTEVP3YJU/60kFU7Gk/vkJIYz4nqWjKTw9cIcKbdnOWpEanI5XQfgfJRWVPLH97fBMAzN46kfUoikwZ05e21e3QSrwjjmwSa4kkCHv3vbXz370VDugUlCXjM+sG4eq/P6tOZz391EeNOsxa3SW8XvkTQLiGeAd0yAHSNgginiSCC3PSsteDIvZcM9DYFTOjXlT3HKli+vW1z0Kvg6p+d4Xf7b746jFH5Ld8gP+Obo4Idkpfnw/eGMT158dtjSEmK9651XF0X3i8Uv716OABz1+9p4UjlJE0EEaKyppZPt1j3290yrsC7/Sun55CRnMBLi3c4FZrCmr5h074yznvyQ3737ga+9GnrH9j95FDRqwt78Oevj+ThK4Y0ea4tj1zMxBAOqZxz+3genDKYn180wLvttkmnMbqgE+NPy2rmncE3OCeTgqw0Xlu+M6zXVa2jiSBCeDr3HvzKEOJ8bsBJTUrg/EHZvLqshGM+88+o0NtfWsm1zyzi3fV7GXL/O5z35Ids2lfm7fj1eOra4Xzr7AKuOsNq4umS0a7RXbv/+pbVcfr0dSPq/f+GQlyccOPYfDKST44Q6pudwX++M5aOaeHrLAZrOoyrC/NYvPUQWw+Uh/XaKnCaCCJARXWtt835hjN7Ntp/4eBuAAx7YG5Y44pl015b3eTIGoB9pRWMevhdPtt6iFsbzNHj8fHPJvL7609nQLdM7rt0EE9cM7ze/he/fab3+fi+XazO3BG5wfkBoognQf7kPysdjkQ1xen7CBTw5urdgNVu7O9mH08iAOsu0VvGFfDLywY1Ok41b++xCn737kbuuqAfM5dYTW39szMY2D2TjYdr2b90Bz99dTUAvbs0nqqhoR6dUhsN2/R1Vp/wNsNEKs/9E8u3H6G8soa0MHZYq8Do/4jDjDH8ctZa+nRJo6h/lyaPW/3ABd4awbMLtzJlRA7DQzCffCz77Ttf8sqyEl5afHLEz89eW33ygM9OPt+y338zxk8v7E/vrDTOH5Qd0DUfv2poQDeSxbrUpHiOV9VyxZ8XMvfOCU6HoxrQROCwN1bt4nhVLVNG5DZ7639mciKf/+oiXl9Rwj3/XcuUPy3kP98Zy+iC5qcecDvPPDv/uKmQV3wWZQnUZcNzuH50D84s6MyslTu5ZFh3712zgbh2VOOmPjf66GcTKXzoXTbsDf6qaKrttI/AYXfMtNpNvz2+d4vHpiTF8/Uze3GR3VR0zTOLyJ82O6AFSdzIdynGbz3vv52/OeP7ZvHE1cM5q08W8XHClWfktSoJqJOy0tt5h7VuP3jc4WhUQ5oIHHSg7OSskylJgX/A/PXGkdx5Xj/v6yH3v6M3nPmxooV7L2ZOHcNfbxjJ09eNaLTv118dxvM3j6435YNqm2duHAnAg7Ob7qRXztCmIQc9ZI9aeeHWM1s4srE7zutLbscU7nplFWCNKFr9wAVhmVkyGpRV1ng7fn1NHtKNiQO6csnQ7vU6LY9s/4JrLiriyIkqMpMTtUMzBHp1tjrg563fS1VNnSbZCKL/Ew6pqa3jf/ZsleNO8Safr47MY9PDkwEorayh4O45VNbUBi3GaDJ33R5eXrKdZdsOkT9tNkPuf8e7b6JPJ/z9lw3mmsIejT7oe2XGk5IUT/f2KZoEQui5m607qv+3Qm8wiyT6G++QxVutu4jvvWRgm86TEB/Hhocm0+9ea9Kz/ve+7Zppq1ftOMJNMxZzz8UD/X77B8hITuAvN4zkYHkVK7cfoVt7HcHjpAn9ujAkN5O/fLiZq0bmER+Dq5dFI60ROGTWyl2kJsXzNT83kLVWUkIcmx+52Pv6zpdj+8adOnu+nEff+pwjx6ubTAI/Pr8fax64kOTEeHI7pHDJsO7hDFP5ISL8oOg0th4o5621u50OR9k0ETjgWEU1Ly/dwZCc9qQmBadSFh8nrPzl+QD8d8VOdh45EZTzhlttneF4VeNRUNW1deRPm83IB+fR+xdzWLH9sHduJn82PjyZ28/tG8pQ1Sm6cHA3+nRJ408fbNZBDhFCE4EDpn+4BYDrz+wR1PN2SE3inoutpqbrpi8K6rnDwRhDn1/MYdAv36HvPXPInzabB95YB8BKe479g/Z0xlf8+ZN67/3OhN7e5R1vGtur0bq9KnLExQnfKzqNz3cf44MvdXXaSKB9BA7416fbAJgyPPjzznz7nN48PLjX9pAAABOsSURBVOdzdhw6Qf602Wx8eHLUfCgeOX5yUr3qWuub4nOfFLe4wpVvn8iGhyaToO3OEW/KiByemreBP76/iYn9u+poN4dFxydEDDl6vJqjJ6rpn50Rslkol917nvf5PxdtY9O+0maOjhwfbdzf4jGem5I8Fk6bVO91UkJcyGf3VG2XGB/Hdyb0Zvn2I3y2tekmPhUemgjC7H8rrWFz910auknjOqe3Y9X9FwDW9NbnPfkRu49GTp+Bp1142bZDPDn3S/70wSaWbz/Me5+33Ezw1h3jWXLPefxgYh82P3IxuR1SQh2uCpFrCnsQHydcN/1T7StwmDYNhdknmw+QECeM7dM5pNdpn5JIh9REb3PL2Eff583bzmZwTmYL7wyd2jqrDwDgd9eO4Ed+RjdlJCdQWmF1Fv/0wv58b0IfetvvueuCfogIXTLa8dMLBzR6r4ouyYnxXDeqBy98tp1/LNjKrQFMs6JCQxNBGB09Uc37X+zjhjG9wjJ++pNpk3hy7gb+vmArAJf+YQE3jOlJUaZhyp8WcvfkAYzp3baEVFFdS3Ji4+kxFm89xDXPLOK8gdnsL6tstKC6vyQAUFpRwxcPXkRS/MkmnhX3nc/b6/Zw3ajgdq4r591/2WBe+Gw7D83+nG+Mzde7jR2iiSCMPtywn+paw2XDwzOePTUpgXsvHcRtk/oy/FfWFNZrSo5CRg2rdhznuumfsuaBC/hk80FeX15CZnIiu49W8NzNo7xrJgO8uqyE/7dqF8/fMpoPN+xneF57fvjiChZsOgDA9BtHEifC919YTp0x1Pisi/vu53ubjfGcfl1Yuf0wx+xaQGK8NEosHdOSuH60zuIZi5IS4vj5RQN4/O0veHVZSVDuq1Gtp4kgjN5dv5es9CRG9Gh5cfNgap+ayNPXjeAXr69hVclRVvnsG+pn1bNnF27lmsIePLtgKz+c1Nc7n9F9/1vrHfHka+q/lrUqnj9+7XR++OIK3rpjfL31fv+1qJjCfJ1W222+O6E3/1pUzC/+u4ZJA7rq3d8O0HpYmFTX1vHBl/uYNKCrI7fVTxmRG3Az0CNzvmDEr+bx+/c3eaeuAPwmgeY8euVQ7/PbJ50GWFM7Xzosh+LHLqmXBABuHJvfaJuKfSLCXRf2B2DMo+85HI07aY0gTN5dv5fSihomDejqWAy/uXo4E37zAaUVNbRPSeS2Safxt4+30DsrnUevHMrzi4qZsbDY73uzM9ux91hlvW1XnpHL5cNz+OaMJdwyroDO6Unkd04jJSmOSQOsFbx8m3R+fEH/UP1oKspdeUYeP/6PVfOctXKnK9d2dpImgjBZus2aG7+tnbNt0SktiTUPXMgb73zAmLPG0jUjud5IjfsvG8yd5/dj5+ETPDF3A49fNZSM5EREICFOqKqto7bOkJqUQF2d8XbmumWSOxVab90xnslPf8wdM1dy+fAcvcksjLRpKEwWbT7I6IJOdEhNcjoUMtsJXTP8t8NmJicysHsmf7+pkM7p7UhKiCMxPg4RoV1CvHduJL1pSwXbwO6Z3DDGqkH+rImJBFVoaCIIgz1HK1i/+xgT+jW9OL1SCv7v8iEAvLKshB2HdEnLcNFEEAaLtljDLMf01hExSjUnPk540V6x78LffURtnd5xHA4RkQhEpFhE1ojIShFp/SrjEW7R5oO0T0kM+7BRpaLRWadl8f2iPhyvquVrf/vU6XBcISISgW2iMWaEMabQ6UCC7ZPNBxnbu7OuxqRUgO6yR5h9tvUQ89Y3f1OiartISgQxafvB45QcPsFZpzk3WkipaBMXJ/zrW6MB+PY/l7Ju11GHI4ptEgmz/onIVuAwYIBnjDHT/RwzFZgKkJ2dPXLmzJktnresrIz09PQgR9s6H+6oZsa6Kh45O4Wc9MjIu5FQLpFGy8Q/p8tl8Z4a/rzSun/l2QtTiYuQIaVOl8upmjhx4jJ/rS6RkghyjTE7RaQrMA+4zRjzUVPHFxYWmqVLW+5KmD9/PkVFRcEL9BTc/tIKFm05yOJfnBsx46IjoVwijZaJf5FQLmc+8q73Zsatj14cEX9HkVAup0JE/CaCiPiKaozZaf+7D/gvMNrZiILDGMPirYcYXdApIn55lYpGn959rvd5wd1zHIwkdjmeCEQkTUQyPM+BC4C1zkYVHF/uLWXPsQrO6KmjhZQ6VSLCxocne1/nT5vtYDSxyfFEAGQDC0RkFbAYmG2MedvhmIJi0eaDABT20kSgVFskxscx985zvK//s3SHg9HEHscTgTFmizFmuP0YbIx52OmYgmXD3jIykxMYltfe6VCUinr9sjOY9YNxgDUFxdV//cThiGKH44kglr20eDvHKmq0f0CpIBneowNPXjMcgCXFh71rZai20UQQIuWV1opbE/vr/EJKBdOVZ+Txl6+fAVir5415RNcwaCtNBCGybtcxAG4c28vhSJSKPZOHdufxq6yFj/Ycq2Do/e84HFF000QQIqtLrMXah+Rq/4BSoXDtqJ78+1vWBHWllTXkT5tNRXWtw1FFJ00EIbJm51G6t09uct5/pVTbnd03iyX3nOd9PeC+t1m7U6ejaC1NBCGypuQoQ7U2oFTIdclox+ZHLva+vvQPC/jBC8sdjCj6aCIIgWMV1Ww5UK7DRpUKk/g4ofixS+iXbc3/M3vNbvKnzebo8WqHI4sOmghCwFM1HZrXweFIlHKXuXdOYP5dRd7Xw381l3N+/YFzAUUJTQQhsKbETgTaNKRU2OVnpfHFgxd5awfbDx0nf9psth0sdziyyKWJIARW7zxKXscUOqU5v1C9Um6UnBjP3Dsn8J/vjPVum/Cb+TpPURM0EYTAmpKj2j+gVAQYXdCJ4scu4YyeJ5tp86fNZsJvPiASpuCPFJoIguzI8Sq2HzrO0FztH1AqUrz+/XEsnDbJ+3rbweMU3D2H+2fFxETHbaaJIMjW2B3FWiNQKrLkdkih+LFLmH372d5tzy/aRv602dYIoxPuHWGkiSDIPIlgSI4mAqUi0eCc9hQ/dgl3nNu33vbh/zeX/Gmzqa6tcygy52giCLIv95SS0z6Z9qmJToeilGrGnef3o/ixS/juhD71tve95y3yp83mQFmlQ5GFX4LTAcSaL/eU0r9bhtNhKKUCNG3yAKZNHsDGvaWc/9TJpdILH3oXgAn9uvDczaNiejp5rREEUXVtHZv3l9G/W6bToSilWqlvdgbFj13Cm7edXW/7hxv2U3D3HPKnzWb7weMORRdaWiMIoq0HyqmuNQzQGoFSUWtIrtWHAHDNM4tYvPWQd985vzl5l/JHQ4/Ts3Nq2OMLBU0EQfTFnlIAbRpSKkZ4bkgrOXycsx+vP1WFb1JYdPckurdPCWtswaSJIIg27CklPk7o0yXd6VCUUkGU1zHVW0v4fPcxJj/9cb39Yx993/v899efzuXDc8IaX1tpIgiirQfK6dkplaQE7XpRKlYN7J7JcxelUVRUxMJNB/j63z+rt//2l1Zw+0srAOia0Y4P7ioirV1kf9RGdnRRZtuhcnrFSJuhUqpl407L8tYU3v9iL7c8t7Te/n2llQz2WUazQ2oi/7hpFCN6dCA+LnJGIWkiCJLaOsPmfeWMyu/kdChKKQdMGpDtTQr7jlXw+oqdPPbWF/WOOXK8mqv+8on3dVJ8HK9+bywDu2eSGO9cS4ImgiDZffQEJ6pr6ZetHcVKuV3XzGS+O6FPvZvVXl6ynZ+/tqbecVW1dVz+x4X1tp3TrwtPXjOcrPR2YYkVNBEEzTZ7fHGvTto0pJRq7NpRPbl2VE8A6uoMLyzezn3/azzp3Ucb9ntvZvN47XtjGdkrdK0NmgiC5OUlOwDI7Ri9Q8iUUuERFyfcOKYXN47p5d129EQ1w/9vrt/jr/rLonqvO6UlMf3GkQzOaU9KUnyb49FEECRvrNoFEPGjA5RSkal9SqK3jwHAGMMrS0v42WurGx17qLyKr/51UaPtT107nAsGdWv155B+agXJLeMKeHbh1rC26ymlYpeIcM2oHlwzqgdgJYb3v9jHt55f2uR77nx5FbCq0fYZN4/izIKmm5Y0EQRJfBykBqGKppRS/ogI5w7MrldrADhcXsVfP9rMMx9uafK9N89Y0uy5IyIRiMhFwNNAPPB3Y8xjDofUatW1JqLGBSul3KFjWhJ3Tx7I3ZMHereVVdYwc/F2Zi7ZwaZ9ZS2ew/FEICLxwJ+A84ESYImIvGGMWe9EPMerajh6orrV84bU1hlHxwErpZRHersEbh3fm1vH9663XR73f3wkfHKNBjYZY7YYY6qAmcAUp4IZ9Mt36s0bEqiaujqtESilopLjNQIgF9jh87oEOLPhQSIyFZgKkJ2dzfz581s8cVlZWUDH+dPa9+3YWUltde0pXy+c2lIusUrLxD8tF/9irVwiIREExBgzHZgOcMbIkWbkmHFU1dRRXWsoq6ymW/sUKqpr6ZyWxKHyKg6VV7Fq+RJGnn4mXdLbUV1r2FdaQVJCHGntEqitNSQlxJGSGE9NnaG6to6UxHh4ew4AO9rlc4PPGF9jwLNAkYhQVVNHQpwQFyfU1Nbxn50rSD9xjKKionAXTavNnz8/KuIMJy0T/7Rc/Iu1comERLAT6OHzOs/e1qR1u44x9AH/N1408mHrm3kA7pu1jvtmrTul9yqlVDSJhD6CJUBfESkQkSTgOuANh2M6JZcO6+50CEop1WqO1wiMMTUi8kPgHazho88aY5r9Kj40tz1LG4yl9edUqm/502YD8POLBvC9oj4tHK2UUtHP8UQAYIyZA8xxOg5ft44vcDoEpZQKi0hoGoooT1w9nPMGZus9AUop14iIGkEkuWpkHleNzHM6DKWUChv92quUUi6niUAppVxOE4FSSrmcJgKllHI5TQRKKeVymgiUUsrlNBEopZTLaSJQSimXE2OM0zG0mojsB7YFcGgWcCDE4UQjLZfGtEz803LxL1rLpZcxpkvDjVGZCAIlIkuNMYVOxxFptFwa0zLxT8vFv1grF20aUkopl9NEoJRSLhfriWC60wFEKC2XxrRM/NNy8S+myiWm+wiUUkq1LNZrBEoppVqgiUAppVwuahKBiFwkIl+KyCYRmeZnfzsRedne/5mI5Pvsu9ve/qWIXBjoOaNBiMqlWETWiMhKEVkanp8kuE61XESks4h8ICJlIvLHBu8ZaZfLJhH5vYhIeH6a4AhRmcy3z7nSfnQNz08TPG0ol/NFZJn9O7FMRCb5vCe6fleMMRH/wFrUfjPQG0gCVgGDGhzzfeCv9vPrgJft54Ps49sBBfZ54gM5Z6Q/QlEu9r5iIMvpn8+hckkDzga+C/yxwXsWA2MAAd4CJjv9s0ZAmcwHCp3++Rwql9OBHPv5EGBntP6uREuNYDSwyRizxRhTBcwEpjQ4ZgrwvP38VeBcOwtPAWYaYyqNMVuBTfb5AjlnpAtFucSCUy4XY0y5MWYBUOF7sIh0BzKNMZ8a6y/9n8BXQvpTBFfQyyRGtKVcVhhjdtnb1wEpdu0h6n5XoiUR5AI7fF6X2Nv8HmOMqQGOAp2beW8g54x0oSgXAAPMtau7U0MQd6i1pVyaO2dJC+eMZKEoE48ZdrPQfRHfBNJYsMrlKmC5MaaSKPxd0cXrlT9nG2N22u2980TkC2PMR04HpSLS1+3flQzgNeBGrG/AriEig4HHgQucjuVURUuNYCfQw+d1nr3N7zEikgC0Bw42895AzhnpQlEuGGM8/+4D/kv0NRm1pVyaO2deC+eMZKEoE9/flVLgRVz2uyIieVh/I98wxmz2OT6qfleiJREsAfqKSIGIJGF12LzR4Jg3gJvs518F3rfb594ArrPb7gqAvlgdOYGcM9IFvVxEJM3+doeIpGF9y1kbhp8lmNpSLn4ZY3YDx0RkjN388Q1gVvBDD5mgl4mIJIhIlv08EbgUF/2uiEgHYDYwzRiz0HNwVP6uON1bHegDuBjYgNXDf4+97VfA5fbzZOAVrE7PxUBvn/feY7/vS3x67/2dM9oewS4XrNETq+zHOpeWSzFwCCjDat8dZG8vxPqg2wz8EfvO/Gh5BLtMsEYTLQNW278rT2OPPIumx6mWC3AvUA6s9Hl0jcbfFZ1iQimlXC5amoaUUkqFiCYCpZRyOU0ESinlcpoIlFLK5TQRKKWUy2kiUEopl9NEoFzNnmLZM4XyHhHZ6fP6kxBc75sisl9E/t7MMSn29as8N2wpFUo615ByNWPMQWAEgIg8AJQZY34b4su+bIz5YTMxnQBGiEhxiONQCtAagVJNEpEy+98iEflQRGaJyBYReUxEvi4ii+3FR/rYx3URkddEZIn9GBfANQbb51kpIqtFpG+ofy6lGtIagVKBGQ4MxJpmYQvwd2PMaBG5A7gN+BHWFAtPGWMWiEhP4B37Pc35LvC0MeYFe66b+JD9BEo1QROBUoFZYqzJxBCRzcBce/saYKL9/DxgkM+U/Jkikm6MKWvmvIuAe+xZLF83xmwMfuhKNU+bhpQKTKXP8zqf13Wc/EIVB4wxxoywH7ktJAGMMS8ClwMngDm+694qFS6aCJQKnrlYzUQAiMiIlt4gIr2BLcaY32NNVTwsdOEp5Z8mAqWC53ag0O70XY/V/t+Sa4C1IrISawF0V63upSKDTkOtVBiJyDeBwuaGj/ocW2wfeyDUcSl30xqBUuF1ApgcyA1lQCJWH4RSIaU1AqWUcjmtESillMtpIlBKKZfTRKCUUi6niUAppVzu/wMM4wE8d5LOIAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "\n",
    "shot_no = 33993\n",
    "\n",
    "data_URL = \"http://golem.fjfi.cvut.cz/shots/{shot_no}/DASs/StandardDAS/{identifier}.csv\"  \n",
    "ds = np.DataSource(destpath=\"/tmp\")\n",
    "\n",
    "def open_remote(shot_no, identifier, url_template):\n",
    "    return ds.open(url_template.format(shot_no=shot_no, identifier=identifier))\n",
    "\n",
    "def read_signal(shot_no, identifier): \n",
    "    file = open_remote(shot_no, identifier, data_URL)\n",
    "    return pd.read_csv(file, names=[\"Time\",identifier],\n",
    "                     index_col=\"Time\", squeeze=True)  \n",
    "\n",
    "loop_voltage = read_signal(shot_no, \"LoopVoltageCoil_raw\")\n",
    "ax = loop_voltage.plot(grid=True)\n",
    "ax.set(xlabel=\"Time [s]\", ylabel=\"$U_l$ [V]\", title=\"Loop voltage $U_l$ {}\".format(shot_no));\n",
    "\n",
    "plt.savefig(\"graph.jpg\")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
