notebooks/example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import logging\n",
    "logging.basicConfig(level=logging.INFO)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: sentinelloader in /usr/local/lib/python3.6/dist-packages/sentinelloader-1.1.1-py3.6.egg (1.1.1)\n",
      "Requirement already satisfied: uuid>=1.3.0 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (1.30)\n",
      "Requirement already satisfied: gdal>=2.2.2 in /usr/lib/python3/dist-packages (from sentinelloader) (3.0.2)\n",
      "Requirement already satisfied: requests>=2.21.0 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (2.22.0)\n",
      "Requirement already satisfied: pandas>=0.24.1 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (0.25.3)\n",
      "Requirement already satisfied: geopandas>=0.4.1 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (0.6.2)\n",
      "Requirement already satisfied: shapely>=1.6.4 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (1.6.4.post2)\n",
      "Requirement already satisfied: sentinelsat>=0.12.2 in /usr/local/lib/python3.6/dist-packages (from sentinelloader) (0.13)\n",
      "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21.0->sentinelloader) (2019.11.28)\n",
      "Requirement already satisfied: urllib3!=1.25.0,!=1.25.1,<1.26,>=1.21.1 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21.0->sentinelloader) (1.25.7)\n",
      "Requirement already satisfied: idna<2.9,>=2.5 in /usr/lib/python3/dist-packages (from requests>=2.21.0->sentinelloader) (2.6)\n",
      "Requirement already satisfied: chardet<3.1.0,>=3.0.2 in /usr/local/lib/python3.6/dist-packages (from requests>=2.21.0->sentinelloader) (3.0.4)\n",
      "Requirement already satisfied: pytz>=2017.2 in /usr/local/lib/python3.6/dist-packages (from pandas>=0.24.1->sentinelloader) (2019.3)\n",
      "Requirement already satisfied: numpy>=1.13.3 in /usr/local/lib/python3.6/dist-packages (from pandas>=0.24.1->sentinelloader) (1.18.1)\n",
      "Requirement already satisfied: python-dateutil>=2.6.1 in /usr/local/lib/python3.6/dist-packages (from pandas>=0.24.1->sentinelloader) (2.8.1)\n",
      "Requirement already satisfied: fiona in /usr/local/lib/python3.6/dist-packages (from geopandas>=0.4.1->sentinelloader) (1.8.13)\n",
      "Requirement already satisfied: pyproj in /usr/local/lib/python3.6/dist-packages (from geopandas>=0.4.1->sentinelloader) (2.4.2.post1)\n",
      "Requirement already satisfied: geomet in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (0.2.1.post1)\n",
      "Requirement already satisfied: six in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (1.13.0)\n",
      "Requirement already satisfied: geojson>=2 in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (2.5.0)\n",
      "Requirement already satisfied: html2text in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (2020.1.16)\n",
      "Requirement already satisfied: click in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (7.0)\n",
      "Requirement already satisfied: tqdm in /usr/local/lib/python3.6/dist-packages (from sentinelsat>=0.12.2->sentinelloader) (4.41.1)\n",
      "Requirement already satisfied: click-plugins>=1.0 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader) (1.1.1)\n",
      "Requirement already satisfied: attrs>=17 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader) (19.3.0)\n",
      "Requirement already satisfied: munch in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader) (2.5.0)\n",
      "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.6/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader) (0.5.0)\n",
      "\u001b[33mWARNING: You are using pip version 19.3.1; however, version 21.3.1 is available.\n",
      "You should consider upgrading via the 'pip install --upgrade pip' command.\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "!pip install sentinelloader"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from sentinelloader import Sentinel2Loader"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Uninstalling sentinelloader-1.1.1:\n",
      "  Successfully uninstalled sentinelloader-1.1.1\n",
      "Collecting git+https://github.com/flaviostutz/sentinelloader\n",
      "  Cloning https://github.com/flaviostutz/sentinelloader to /tmp/pip-req-build-14poc0nk\n",
      "Requirement already satisfied: uuid>=1.3.0 in /usr/local/lib/python3.5/dist-packages/uuid-1.30-py3.5.egg (from sentinelloader==1.1.1) (1.30)\n",
      "Requirement already satisfied: gdal>=2.2.2 in /usr/lib/python3/dist-packages (from sentinelloader==1.1.1) (2.2.2)\n",
      "Requirement already satisfied: requests>=2.21.0 in /usr/local/lib/python3.5/dist-packages (from sentinelloader==1.1.1) (2.21.0)\n",
      "Requirement already satisfied: pandas>=0.24.1 in /usr/local/lib/python3.5/dist-packages (from sentinelloader==1.1.1) (0.24.2)\n",
      "Requirement already satisfied: geopandas>=0.4.1 in /usr/local/lib/python3.5/dist-packages (from sentinelloader==1.1.1) (0.4.1)\n",
      "Requirement already satisfied: shapely>=1.6.4 in /usr/local/lib/python3.5/dist-packages (from sentinelloader==1.1.1) (1.6.4.post2)\n",
      "Requirement already satisfied: sentinelsat>=0.12.2 in /usr/local/lib/python3.5/dist-packages/sentinelsat-1.1.0-py3.5.egg (from sentinelloader==1.1.1) (1.1.0)\n",
      "Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.5/dist-packages (from requests>=2.21.0->sentinelloader==1.1.1) (2019.3.9)\n",
      "Requirement already satisfied: chardet<3.1.0,>=3.0.2 in /usr/local/lib/python3.5/dist-packages (from requests>=2.21.0->sentinelloader==1.1.1) (3.0.4)\n",
      "Requirement already satisfied: idna<2.9,>=2.5 in /usr/local/lib/python3.5/dist-packages (from requests>=2.21.0->sentinelloader==1.1.1) (2.8)\n",
      "Requirement already satisfied: urllib3<1.25,>=1.21.1 in /usr/local/lib/python3.5/dist-packages (from requests>=2.21.0->sentinelloader==1.1.1) (1.24.1)\n",
      "Requirement already satisfied: python-dateutil>=2.5.0 in /usr/local/lib/python3.5/dist-packages (from pandas>=0.24.1->sentinelloader==1.1.1) (2.8.0)\n",
      "Requirement already satisfied: pytz>=2011k in /usr/local/lib/python3.5/dist-packages (from pandas>=0.24.1->sentinelloader==1.1.1) (2019.1)\n",
      "Requirement already satisfied: numpy>=1.12.0 in /usr/local/lib/python3.5/dist-packages (from pandas>=0.24.1->sentinelloader==1.1.1) (1.16.2)\n",
      "Requirement already satisfied: fiona in /usr/local/lib/python3.5/dist-packages (from geopandas>=0.4.1->sentinelloader==1.1.1) (1.8.6)\n",
      "Requirement already satisfied: pyproj in /usr/local/lib/python3.5/dist-packages (from geopandas>=0.4.1->sentinelloader==1.1.1) (2.1.3)\n",
      "Requirement already satisfied: click in /usr/local/lib/python3.5/dist-packages (from sentinelsat>=0.12.2->sentinelloader==1.1.1) (7.0)\n",
      "Requirement already satisfied: geojson>=2 in /usr/local/lib/python3.5/dist-packages/geojson-2.5.0-py3.5.egg (from sentinelsat>=0.12.2->sentinelloader==1.1.1) (2.5.0)\n",
      "Requirement already satisfied: geomet in /usr/local/lib/python3.5/dist-packages/geomet-0.3.0-py3.5.egg (from sentinelsat>=0.12.2->sentinelloader==1.1.1) (0.3.0)\n",
      "Requirement already satisfied: html2text in /usr/local/lib/python3.5/dist-packages/html2text-2020.1.16-py3.5.egg (from sentinelsat>=0.12.2->sentinelloader==1.1.1) (2020.1.16)\n",
      "Requirement already satisfied: tqdm>=4.27 in /usr/local/lib/python3.5/dist-packages (from sentinelsat>=0.12.2->sentinelloader==1.1.1) (4.31.1)\n",
      "Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.5/dist-packages (from python-dateutil>=2.5.0->pandas>=0.24.1->sentinelloader==1.1.1) (1.12.0)\n",
      "Requirement already satisfied: munch in /usr/local/lib/python3.5/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader==1.1.1) (2.3.2)\n",
      "Requirement already satisfied: attrs>=17 in /usr/local/lib/python3.5/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader==1.1.1) (19.1.0)\n",
      "Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.5/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader==1.1.1) (0.5.0)\n",
      "Requirement already satisfied: click-plugins>=1.0 in /usr/local/lib/python3.5/dist-packages (from fiona->geopandas>=0.4.1->sentinelloader==1.1.1) (1.1.1)\n",
      "Building wheels for collected packages: sentinelloader\n",
      "  Building wheel for sentinelloader (setup.py) ... \u001b[?25ldone\n",
      "\u001b[?25h  Stored in directory: /tmp/pip-ephem-wheel-cache-0_ovytpj/wheels/50/ec/84/a8c56ef0693bfee8662a78cc6ec6695b8dbffca99284c56e5f\n",
      "Successfully built sentinelloader\n",
      "Installing collected packages: sentinelloader\n",
      "Successfully installed sentinelloader-1.1.1\n",
      "\u001b[33mYou are using pip version 19.0.3, however version 20.3.4 is available.\n",
      "You should consider upgrading via the 'pip install --upgrade pip' command.\u001b[0m\n"
     ]
    },
    {
     "ename": "SyntaxError",
     "evalue": "invalid syntax (sentinel.py, line 287)",
     "output_type": "error",
     "traceback": [
      "Traceback \u001b[0;36m(most recent call last)\u001b[0m:\n",
      "  File \u001b[1;32m\"/usr/local/lib/python3.5/dist-packages/IPython/core/interactiveshell.py\"\u001b[0m, line \u001b[1;32m3291\u001b[0m, in \u001b[1;35mrun_code\u001b[0m\n    exec(code_obj, self.user_global_ns, self.user_ns)\n",
      "  File \u001b[1;32m\"<ipython-input-8-5a7d51863c62>\"\u001b[0m, line \u001b[1;32m5\u001b[0m, in \u001b[1;35m<module>\u001b[0m\n    from sentinelloader import Sentinel2Loader\n",
      "  File \u001b[1;32m\"/notebooks/sentinelloader/__init__.py\"\u001b[0m, line \u001b[1;32m1\u001b[0m, in \u001b[1;35m<module>\u001b[0m\n    from .sentinel2loader import Sentinel2Loader\n",
      "  File \u001b[1;32m\"/notebooks/sentinelloader/sentinel2loader.py\"\u001b[0m, line \u001b[1;32m1\u001b[0m, in \u001b[1;35m<module>\u001b[0m\n    from sentinelsat import SentinelAPI\n",
      "\u001b[0;36m  File \u001b[0;32m\"/usr/local/lib/python3.5/dist-packages/sentinelsat-1.1.0-py3.5.egg/sentinelsat/__init__.py\"\u001b[0;36m, line \u001b[0;32m4\u001b[0;36m, in \u001b[0;35m<module>\u001b[0;36m\u001b[0m\n\u001b[0;31m    from . import sentinel\u001b[0m\n",
      "\u001b[0;36m  File \u001b[0;32m\"/usr/local/lib/python3.5/dist-packages/sentinelsat-1.1.0-py3.5.egg/sentinelsat/sentinel.py\"\u001b[0;36m, line \u001b[0;32m287\u001b[0m\n\u001b[0;31m    self.logger.info(f\"Found {count:,} products\")\u001b[0m\n\u001b[0m                                               ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m invalid syntax\n"
     ]
    }
   ],
   "source": [
    "!pip uninstall -y sentinelloader\n",
    "\n",
    "# used to test release\n",
    "!pip install git+https://github.com/flaviostutz/sentinelloader\n",
    "from sentinelloader import Sentinel2Loader\n",
    "\n",
    "# used during local development\n",
    "# from sentinelloader.sentinel2loader import Sentinel2Loader"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from time import time\n",
    "logger = logging.getLogger('sentinelloader')\n",
    "\n",
    "class Timer:\n",
    "    def __init__(self, name, debug=True):\n",
    "        self._name = name\n",
    "        self._debug = debug\n",
    "        self.start()\n",
    "        self._lastElapsed = None\n",
    "    \n",
    "    def start(self):\n",
    "        self._start = time()\n",
    "        if(self._debug):\n",
    "            logger.info('> [started] ' + self._name + '...')\n",
    "\n",
    "    def stop(self):\n",
    "        self._lastElapsed = (time()-self._start)\n",
    "        if(self._debug):\n",
    "            logger.info('> [done]    {} ({:.3f} ms)'.format(self._name, self._lastElapsed*1000))\n",
    "            \n",
    "    def elapsed(self):\n",
    "        if(self._lastElapsed != None):\n",
    "            return (self._lastElapsed)\n",
    "        else:\n",
    "            return (time()-self._start)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.patches as mpatches\n",
    "from osgeo import gdal, osr\n",
    "import cartopy.crs as ccrs\n",
    "import numpy as np\n",
    "\n",
    "def showGeoTiffs(geoTiffFiles, geometries=None, labels=None, cols=4, size=6, cmap=None, cmap_min=None, cmap_max=None, interpolation=None, legendHandles=None, group_by_label=False, name='image', output_dir=None, transformData=None):\n",
    "    \"\"\"geometries - list of geoseries from GeoPandas\n",
    "       labels - list os strings\n",
    "       getTiffFiles - list of tiff file paths\n",
    "    \"\"\"\n",
    "    logger.info('showing ' + str(len(geoTiffFiles)) + ' images')\n",
    "    fig = plt.figure()\n",
    "    rows = int(len(geoTiffFiles)/cols)+1\n",
    "    t = Timer('generating image patches. rows=' + str(rows) + '; cols=' + str(cols))\n",
    "    fig.set_size_inches(cols*size, rows*size)\n",
    "\n",
    "    image_indexes = range(len(geoTiffFiles))\n",
    "\n",
    "    #order indexes by label\n",
    "    if(group_by_label==True and image_labels!=None):\n",
    "        index_label_map = []\n",
    "        for i,label in enumerate(image_labels):\n",
    "            index_label_map.append((i,label))\n",
    "        label_image_map = np.array(index_label_map, dtype=[('index',int),('label',int)])\n",
    "        label_image_map = np.sort(label_image_map, order='label')\n",
    "        image_indexes = []\n",
    "        for a in label_image_map:\n",
    "            image_indexes.append(a[0])\n",
    "\n",
    "    c = 0\n",
    "    for i in image_indexes:\n",
    "        ds = gdal.Open(geoTiffFiles[i])\n",
    "        data = ds.ReadAsArray()\n",
    "        gt = ds.GetGeoTransform()\n",
    "        proj = ds.GetProjection()\n",
    "\n",
    "        inproj = osr.SpatialReference()\n",
    "        inproj.ImportFromWkt(proj)\n",
    "\n",
    "        projcs = inproj.GetAuthorityCode('PROJCS')\n",
    "        projection = ccrs.epsg(projcs)\n",
    "\n",
    "        ax = fig.add_subplot(rows,cols,c+1, projection=projection)\n",
    "\n",
    "        extent = (gt[0], gt[0] + ds.RasterXSize * gt[1],\n",
    "                  gt[3] + ds.RasterYSize * gt[5], gt[3])\n",
    "\n",
    "        \n",
    "        #seems like this data has multiple channels and the channels area indexed from the first position. invert it\n",
    "        if data.shape[0]<5:\n",
    "            data = data[:3, :, :].transpose((1, 2, 0))\n",
    "        else:\n",
    "            if transformData!=None:\n",
    "                data = transformData(data)\n",
    "\n",
    "        if data.dtype in [np.uint8, np.uint16, np.uint32]:\n",
    "            if cmap_min==None:\n",
    "                cmap_min = np.iinfo(data.dtype).min\n",
    "            if cmap_max==None:\n",
    "                cmap_max = np.iinfo(data.dtype).max\n",
    "#         if(cmap==None):\n",
    "#             data = data.astype('uint8')\n",
    "                \n",
    "        ax.imshow(data, alpha=1.0, transform=projection, extent=extent, cmap=cmap, vmin=cmap_min, vmax=cmap_max, interpolation=interpolation, origin='upper')\n",
    "        if legendHandles!=None:\n",
    "            ax.legend(handles=legendHandles, bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)\n",
    "\n",
    "        if(geometries!=None):\n",
    "            df_epsg = geometries[i].to_crs(epsg=ax.projection.epsg_code)\n",
    "            df_epsg.plot(ax=ax, facecolor='none', edgecolor='black', linestyle='--', linewidth=2)\n",
    "        \n",
    "        if(labels!=None):\n",
    "            ax.text(0.5, 0, str(labels[i]), horizontalalignment='center', verticalalignment='bottom', transform=ax.transAxes, fontsize=14, style='normal', color='red')\n",
    "\n",
    "        c = c + 1\n",
    "    \n",
    "    if(output_dir!=None):\n",
    "        f = output_dir + name + '.jpg'\n",
    "        plt.savefig(f)\n",
    "        plt.close(fig)\n",
    "    else:\n",
    "        plt.show()\n",
    "        \n",
    "    t.stop()\n",
    "\n",
    "def sclColors():\n",
    "    colors = ['black','red','dimgray','brown','green','yellow','blue','dimgray','darkgray','lightgray','skyblue','magenta']\n",
    "    return colors\n",
    "\n",
    "def sclLegendPatches():\n",
    "    colors = sclColors()\n",
    "    labels = ['no data','defective','dark area','cloud shadows','vegetation','not vegetated','water','unclassified','cloud medium prob','cloud high prob','thin cirrus','snow']\n",
    "    patches = [ mpatches.Patch(color=colors[i], label=labels[i]) for i in range(11) ]\n",
    "    return patches\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:Getting region history for band TCI from 2019-12-22 to 2020-02-27 at 60m\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "POLYGON Z ((-55.7005254871978 -13.3255356572483 0, -55.6991419460571 -13.3476412354778 0, -55.6994626679625 -13.3482968267724 0, -55.6993622255233 -13.3497238029627 0, -55.6990879718946 -13.3512050194787 0, -55.6987785623234 -13.3516544880401 0, -55.7009726154231 -13.3521465384388 0, -55.7031391574231 -13.3478350497349 0, -55.70365985481779 -13.3477166261783 0, -55.7042048225458 -13.3481784524516 0, -55.7053983895748 -13.3485070297582 0, -55.7071743431305 -13.350239177507 0, -55.7073679164227 -13.3508664639014 0, -55.707297910641 -13.3536846078021 0, -55.70845582938671 -13.3547391417037 0, -55.7083650623632 -13.3558936045135 0, -55.70920687339019 -13.3573796093408 0, -55.7107993017724 -13.357686953459 0, -55.711524110584 -13.3580807761113 0, -55.71730886861859 -13.3591710955438 0, -55.7198276533255 -13.327129129925 0, -55.7005254871978 -13.3255356572483 0))\n",
      "getRegionHistory 1\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:Couldn't get data for 2019-12-22 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2019-12-27 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-01 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-06 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-11 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-16 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-21 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-26 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-31 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-05 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-10 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-15 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-20 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-25 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:showing 0 images\n",
      "INFO:sentinelloader:> [started] generating image patches. rows=1; cols=10...\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<Figure size 2880x288 with 0 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:> [done]    generating image patches. rows=1; cols=10 (5.841 ms)\n",
      "INFO:sentinelloader:Getting region history for band SCL from 2019-12-22 to 2020-02-27 at 60m\n",
      "INFO:sentinelloader:Couldn't get data for 2019-12-22 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "getRegionHistory 2\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:Couldn't get data for 2019-12-27 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-01 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-06 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-11 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-16 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-21 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-26 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-31 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-05 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-10 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-15 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-20 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-25 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:showing 0 images\n",
      "INFO:sentinelloader:> [started] generating image patches. rows=1; cols=10...\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<Figure size 2880x288 with 0 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:> [done]    generating image patches. rows=1; cols=10 (3.932 ms)\n",
      "INFO:sentinelloader:Getting region history for band NDVI from 2019-12-22 to 2020-02-27 at 60m\n",
      "INFO:sentinelloader:Couldn't get data for 2019-12-22 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2019-12-27 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-01 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-06 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-11 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-16 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-21 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-26 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-01-31 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-05 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-10 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-15 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-20 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:Couldn't get data for 2020-02-25 using the specified filter. err=HTTP status 401 Unauthorized: \n",
      "# HTTP Status 401 – Unauthorized\n",
      "\n",
      "* * *\n",
      "\n",
      " **Type** Status Report\n",
      "\n",
      " **Message** Unauthorized\n",
      "\n",
      " **Description** The request has not been applied because it lacks valid\n",
      "authentication credentials for the target resource.\n",
      "\n",
      "* * *\n",
      "\n",
      "### Apache Tomcat/8.0.53\n",
      "INFO:sentinelloader:showing 0 images\n",
      "INFO:sentinelloader:> [started] generating image patches. rows=1; cols=10...\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<Figure size 2880x288 with 0 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "INFO:sentinelloader:> [done]    generating image patches. rows=1; cols=10 (5.556 ms)\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import logging\n",
    "import os\n",
    "from osgeo import gdal\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib import rcParams\n",
    "import math\n",
    "import geopandas as gpd\n",
    "from cartopy.io.shapereader import Reader\n",
    "from cartopy.feature import ShapelyFeature\n",
    "from cartopy.io.shapereader import Reader\n",
    "from descartes import PolygonPatch\n",
    "import fiona\n",
    "from pprint import pprint\n",
    "import re\n",
    "\n",
    "logging.basicConfig(level=logging.INFO)\n",
    "\n",
    "sl = Sentinel2Loader('/notebooks/data/output/sentinelcache', \n",
    "                    os.environ['COPERNICUS_USER'], os.environ['COPERNICUS_PASSWORD'],\n",
    "                    apiUrl='https://scihub.copernicus.eu/apihub/', showProgressbars=True, cloudCoverage=(0,100))\n",
    "\n",
    "\n",
    "\n",
    "df = gpd.read_file(\"/notebooks/test/sample2.shp\")\n",
    "df.crs = {'init' :'epsg:4326'}\n",
    "\n",
    "coords = df['geometry'][0]\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "d1 = '2019-12-22'\n",
    "d2 = '2020-02-27'\n",
    "# d1 = '2019-01-11'\n",
    "# d2 = '2019-02-06'\n",
    "\n",
    "print(coords)\n",
    "print(\"getRegionHistory 1\")\n",
    "geoTiffs = sl.getRegionHistory(coords, 'TCI', '60m', d1, d2, daysStep=5)\n",
    "\n",
    "geometries = []\n",
    "labels = []\n",
    "for geoTiff in geoTiffs:\n",
    "    geometries.append(df)\n",
    "    d = re.search('[0-9]{4}-[0-9]{2}-[0-9]{2}', geoTiff)\n",
    "    labels.append(d.group(0))\n",
    "\n",
    "showGeoTiffs(geoTiffs, geometries=geometries, labels=labels, cols=10, size=4,)\n",
    "\n",
    "for geoTiff in geoTiffs:\n",
    "    os.remove(geoTiff)\n",
    "\n",
    "  \n",
    "\n",
    "\n",
    "\n",
    "\n",
    "print(\"getRegionHistory 2\")\n",
    "geoTiffs = sl.getRegionHistory(coords, 'SCL', '60m', d1, d2, daysStep=5)\n",
    "\n",
    "geometries = []\n",
    "labels = []\n",
    "for geoTiff in geoTiffs:\n",
    "    geometries.append(df)\n",
    "    d = re.search('[0-9]{4}-[0-9]{2}-[0-9]{2}', geoTiff)\n",
    "    labels.append(d.group(0))\n",
    "\n",
    "showGeoTiffs(geoTiffs, geometries=geometries, labels=labels, cols=10, size=4, cmap=mpl.colors.ListedColormap(sclColors()), cmap_min=0, cmap_max=11, legendHandles=sclLegendPatches())\n",
    "\n",
    "for geoTiff in geoTiffs:\n",
    "    os.remove(geoTiff)\n",
    "\n",
    "\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "geoTiffs = sl.getRegionHistory(coords, 'NDVI', '60m', d1, d2, daysStep=5)\n",
    "\n",
    "geometries = []\n",
    "labels = []\n",
    "for geoTiff in geoTiffs:\n",
    "    geometries.append(df)\n",
    "    d = re.search('[0-9]{4}-[0-9]{2}-[0-9]{2}', geoTiff)\n",
    "    labels.append(d.group(0))\n",
    "\n",
    "def td(data):\n",
    "    return np.negative(data)\n",
    " \n",
    "showGeoTiffs(geoTiffs, geometries=geometries, labels=labels, cols=10, size=4, transformData=td, cmap='viridis', interpolation='none')\n",
    "\n",
    "for geoTiff in geoTiffs:\n",
    "    os.remove(geoTiff)\n",
    "    \n",
    "\n",
    "    \n",
    "\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "# geoTiffs = sl.getRegionHistory(coords, 'NDVI', '60m', '2018-06-11', '2019-03-11', daysStep=5)\n",
    "\n",
    "# geometries = []\n",
    "# labels = []\n",
    "# for geoTiff in geoTiffs:\n",
    "#     geometries.append(df)\n",
    "#     d = re.search('[0-9]{4}-[0-9]{2}-[0-9]{2}', geoTiff)\n",
    "#     labels.append(d.group(0))\n",
    "\n",
    "# def td(data):\n",
    "# #     v[v>200] = 255\n",
    "#     print(np.mean(data))\n",
    "#     return data\n",
    " \n",
    "# showGeoTiffs(geoTiffs, geometries=geometries, labels=labels, transformData=td, cmap='viridis', cmap_min=-1, cmap_max=1)\n",
    "\n",
    "# for geoTiff in geoTiffs:\n",
    "#     os.remove(geoTiff)\n",
    "\n",
    "    \n",
    "    \n",
    "    \n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "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.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}