GeolifePrj_01_df2geojson.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# **GeolifePrj_01_df2gjson.ipynb**\n",
    "### Convierte los DataFrames creados en **GeolifePrj_00_plt2df** (archivos geolife_xxx.zip) a archivos GeoJson simplificados:\n",
    ">### **geolife_xxx_simpl.geojson** (puntos gps=(lat,lon) y su date=\"fecha/hora\" respectiva)\n",
    ">### **geolife_xxx_lines.geojson** (una linea por trayectoria)\n",
    ">### **geolife_xxx_linesimpl.geojson** (lineas despues de simplificar la data)\n",
    ">### **geolife_xxx_colorlabel.geojson** (se le agrega en tag de color segun el label de geolife)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 128,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geopandas\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "#from shapely.geometry import Point\n",
    "from geojson import  Point,Feature, FeatureCollection, dumps\n",
    "\n",
    "import csv, json\n",
    "\n",
    "import seaborn as sns\n",
    "\n",
    "from shapely.geometry import LineString\n",
    "from shapely.geometry import Point\n",
    "\n",
    "from dateutil.parser import parse as parsedt\n",
    "\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.colors as colors\n",
    "from matplotlib import cm"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Mappings para los labels"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### def **plt2geojson(...):** \n",
    ">### Convierte archivo original plt a geojson\n",
    "### def **simplecsv2geojson(...)** \n",
    ">### Convierte un archivo **geolife_xxx_simpl.csv** generado por la funcion **crea_csv_simpl():** a un archivo geojson"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "def plt2geojson(from_plt_file,to_geojson_file):\n",
    "    features = []\n",
    "    with open(from_plt_file, newline='') as csvfile:\n",
    "        reader = csv.reader(csvfile, delimiter=',')\n",
    "        for _ in range(6): next(reader)\n",
    "        for latitude,longitude,empty,altitude,days,date,time in reader:\n",
    "            latitude, longitude = map(float, (latitude, longitude))\n",
    "            features.append(\n",
    "                Feature(\n",
    "                    geometry = Point((longitude, latitude)),\n",
    "                    properties = {\n",
    "                        'altitude': altitude,\n",
    "                        'date': f\"{date} {time}\"\n",
    "                    }\n",
    "                )\n",
    "            )\n",
    "\n",
    "    collection = FeatureCollection(features)\n",
    "    with open(to_geojson_file, \"w\") as f:\n",
    "        f.write(f\"{collection}\")\n",
    "        \n",
    "def simplecsv2geojson(from_csv_file,to_geojson_file):\n",
    "    features = []\n",
    "    with open(from_csv_file, newline='') as csvfile:\n",
    "        reader = csv.reader(csvfile, delimiter=',')\n",
    "        next(reader)\n",
    "        for date,latitude,longitude,label in reader:\n",
    "            latitude, longitude = map(float, (latitude, longitude))\n",
    "            features.append(\n",
    "                Feature(\n",
    "                    geometry = Point((longitude, latitude)),\n",
    "                    properties = {\n",
    "                        'date': f\"{date}\",\n",
    "                        'label': f\"{label}\"\n",
    "                    }\n",
    "                )\n",
    "            )\n",
    "\n",
    "    collection = FeatureCollection(features)\n",
    "    with open(to_geojson_file, \"w\") as f:\n",
    "        f.write(f\"{collection}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### def **df2geojson(...)** \n",
    ">### genera desde **geolife_xxx.zip** un archivo **geolife_xxx.geojson** (TODA LA DATA)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "def df2geojson(df,to_geojson_file):\n",
    "    df=df.astype(str)\n",
    "    lat = df['lat']\n",
    "    lng = df['lon']\n",
    "    #df = df.drop(columns=['lat', 'lon'])\n",
    "    df = df.drop(columns=['lat', 'lon','user'])    \n",
    "    \n",
    "    feat_list = []\n",
    "    failed = []\n",
    "    for i in range(0, len(df.index)):\n",
    "        props = dict(df.iloc[i])\n",
    "        try:\n",
    "            f = Feature(geometry=Point((float(lng.iloc[i]),float(lat.iloc[i]))),\n",
    "                       properties = props)\n",
    "            feat_list.append(f)\n",
    "        except Exception as ex:\n",
    "            print(\"Exception:\",ex)\n",
    "        \n",
    "    collection = FeatureCollection(feat_list)\n",
    "    with open(to_geojson_file, \"w\") as f:\n",
    "        f.write(f\"{collection}\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## def **crea_simpl_csv(**xxx**):** \n",
    ">## genera una version simplificada del DataFrame almacenado en **geolife_xxx.zip** y lo guarda como **geolife_xxx_simpl.csv**\n",
    "### Simplificaciones:\n",
    "### - *Guarda solo 5 digitos decimales para latitude y longitude*\n",
    "### - *Bota la columna 'user' y 'altitude'*\n",
    "### - *Guarda solo maximo UNA coordena gps por minuto (si hay mas de un punto gps por minuto desecha los sobrantes)*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 186,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wall time: 999 µs\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "#creamos una version simplificada en csv\n",
    "def crea_simpl_csv(xxx):\n",
    "    usrdf = pd.read_pickle(f'geolife_{xxx:03d}.zip')\n",
    "    usrdf.columns=['date','latitude','longitude','altitude','label','user']\n",
    "    #dejamos 5 digitos de precision, que equivale a una resolucion de 8m^2 en el mapa. Explicacion: https://planspace.org/20180719-geolocation_precision_by_digit/\n",
    "    usrdf['latitude']=usrdf['latitude'].apply(lambda x:f\"{x:.5f}\")\n",
    "    usrdf['longitude']=usrdf['longitude'].apply(lambda x:f\"{x:.5f}\")\n",
    "    usrdf.drop(columns=['user','altitude'],inplace=True)\n",
    "    usrdf=usrdf.astype(str)\n",
    "    \n",
    "    #elimina los segundos y luego los repetidos, con lo que solo nos quedamos con una resolucion de \n",
    "    #un punto gps por minutos:\n",
    "    usrdf['date']=usrdf['date'].apply(lambda x:x[:-3]+':00')\n",
    "    usrdf = usrdf.drop_duplicates(subset='date', keep='first')\n",
    "    #equivalente más rápido para undersamplear\n",
    "    #usrdf = usrdf.iloc[::12, :]\n",
    "    \n",
    "    usrdf.to_csv(f\"geolife_{xxx:03d}_simpl.csv\", index=False)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Creamos los CSV simplificados"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for xxx in range(0,182):\n",
    "    crea_simpl_csv(xxx)    \n",
    "    ## Muestra avance ;) [########>-------------------------------------] \n",
    "    print(\"[\"+\"#\"*int(100*(1+xxx)/182)+\">\"+\"-\"*int(100*(181-xxx)/181)+\"]\"+f\" {xxx:03d}/182\",end=\"\\r\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Generamos **geolife_xxx_simpl.geojson** a partir de **geolife_xxx_simple.csv**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wall time: 24.7 s####################################################################################>] 181/182\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "#creamos una version simplificada en geojson\n",
    "for usuario_num in range(0,182):\n",
    "    simplecsv2geojson(f\"geolife_{usuario_num:03d}_simpl.csv\",f\"geolife_{usuario_num:03d}_simpl.geojson\")\n",
    "    ## Avance ;)\n",
    "    print(\"[\"+\"#\"*int(100*(1+usuario_num)/182)+\">\"+\"-\"*int(100*(181-usuario_num)/181)+\"]\"+f\" {usuario_num:03d}/182\",end=\"\\r\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Con def **simpl2lines(**xxx**):** vamos a generar los siguientes datasets\n",
    ">### **geolife_xxx_lines.geojson** (una linea por trayectoria (=un día) )\n",
    ">### **geolife_xxx_linesimpl.geojson** (igual pero desde version **_simpl** ) ."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Cargamos un **geolife_xxx_simpl.geojson** y mostramos"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 99,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "geolife_053_simpl.geojson  ->  geolife_053_linesimpl.geojson\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>date</th>\n",
       "      <th>label</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-04-29T08:21:00</td>\n",
       "      <td>0</td>\n",
       "      <td>POINT (116.33038 39.97630)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2008-04-29T09:15:00</td>\n",
       "      <td>0</td>\n",
       "      <td>POINT (116.33007 39.97523)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2008-04-30T13:39:00</td>\n",
       "      <td>3</td>\n",
       "      <td>POINT (116.33028 39.97857)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2008-04-30T13:46:00</td>\n",
       "      <td>3</td>\n",
       "      <td>POINT (116.34682 39.97067)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2008-04-30T13:48:00</td>\n",
       "      <td>0</td>\n",
       "      <td>POINT (116.35105 39.96625)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                  date label                    geometry\n",
       "0  2008-04-29T08:21:00     0  POINT (116.33038 39.97630)\n",
       "1  2008-04-29T09:15:00     0  POINT (116.33007 39.97523)\n",
       "2  2008-04-30T13:39:00     3  POINT (116.33028 39.97857)\n",
       "3  2008-04-30T13:46:00     3  POINT (116.34682 39.97067)\n",
       "4  2008-04-30T13:48:00     0  POINT (116.35105 39.96625)"
      ]
     },
     "execution_count": 99,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# generamos geolife_xxx_linesimpl.geojson donde xxx=\n",
    "xxx=53\n",
    "archivo_gj=f\"geolife_{xxx:03d}_\"\n",
    "archivo_salida_gj = archivo_gj+\"linesimpl.geojson\"\n",
    "archivo_gj=archivo_gj+\"simpl.geojson\"\n",
    "print(archivo_gj,\" -> \" ,archivo_salida_gj)\n",
    "gdf = geopandas.read_file(archivo_gj)\n",
    "gdf.head(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### generamos un objeto **LineString** por cada dia"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 102,
   "metadata": {},
   "outputs": [],
   "source": [
    "list_dates=list(gdf.date.str[:10].unique())\n",
    "gempty= {'date': [], 'geometry': []}\n",
    "gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")\n",
    "for dia in list_dates[1:]:\n",
    "    plst=gdf.loc[gdf.date.str[:10]==dia,'geometry'].to_list()\n",
    "    ls=LineString(plst)\n",
    "    grow = {'date': [dia], 'geometry': [ls]}\n",
    "    g1=geopandas.GeoDataFrame(grow, crs=\"EPSG:4326\")\n",
    "    gdf_new=gdf_new.append(g1)\n",
    "gdf_new['stroke']=\"#f71a3a\"\n",
    "gdf_new['stroke-width']=2\n",
    "gdf_new['stroke-opacity']=0.7"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### vemos el **geoDataFrame** producido"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 103,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>date</th>\n",
       "      <th>geometry</th>\n",
       "      <th>stroke</th>\n",
       "      <th>stroke-width</th>\n",
       "      <th>stroke-opacity</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-04-30</td>\n",
       "      <td>LINESTRING (116.33028 39.97857, 116.34682 39.9...</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-05-01</td>\n",
       "      <td>LINESTRING (119.48115 39.81815, 119.48032 39.8...</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-05-02</td>\n",
       "      <td>LINESTRING (119.47575 39.81952, 119.47877 39.8...</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-05-05</td>\n",
       "      <td>LINESTRING (116.32932 39.97497, 116.33818 39.9...</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2008-05-06</td>\n",
       "      <td>LINESTRING (116.35533 39.96218, 116.35507 39.9...</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         date                                           geometry   stroke  \\\n",
       "0  2008-04-30  LINESTRING (116.33028 39.97857, 116.34682 39.9...  #f71a3a   \n",
       "0  2008-05-01  LINESTRING (119.48115 39.81815, 119.48032 39.8...  #f71a3a   \n",
       "0  2008-05-02  LINESTRING (119.47575 39.81952, 119.47877 39.8...  #f71a3a   \n",
       "0  2008-05-05  LINESTRING (116.32932 39.97497, 116.33818 39.9...  #f71a3a   \n",
       "0  2008-05-06  LINESTRING (116.35533 39.96218, 116.35507 39.9...  #f71a3a   \n",
       "\n",
       "   stroke-width  stroke-opacity  \n",
       "0             2             0.7  \n",
       "0             2             0.7  \n",
       "0             2             0.7  \n",
       "0             2             0.7  \n",
       "0             2             0.7  "
      ]
     },
     "execution_count": 103,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdf_new.head(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Lo guardamos a **geolife_xxx_linesimpl.geojson**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 94,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_new.to_file(archivo_salida_gj, driver=\"GeoJSON\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Ahora metemos todo a la función def **simpl2lines(**xxx**):** y la iteramos\n",
    ">### geolife_xxx_simpl.geojson  ->  geolife_xxx_linesimpl.geojson"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 118,
   "metadata": {},
   "outputs": [],
   "source": [
    "def simpl2lines(xxx):\n",
    "    archivo_gj=f\"geolife_{xxx:03d}_\"\n",
    "    archivo_salida_gj = archivo_gj+\"linesimpl.geojson\"\n",
    "    archivo_gj=archivo_gj+\"simpl.geojson\"\n",
    "    print(archivo_gj,\" -> \" ,archivo_salida_gj)\n",
    "    gdf = geopandas.read_file(archivo_gj)\n",
    "    \n",
    "    list_dates=list(gdf.date.str[:10].unique())\n",
    "    gempty= {'date': [], 'geometry': []}\n",
    "    gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")\n",
    "    for dia in list_dates:\n",
    "        plst=gdf.loc[gdf.date.str[:10]==dia,'geometry'].to_list()\n",
    "        if len(plst)<2:\n",
    "            continue\n",
    "        else:\n",
    "            ls=LineString(plst)\n",
    "            grow = {'date': [dia], 'geometry': [ls]}\n",
    "            g1=geopandas.GeoDataFrame(grow, crs=\"EPSG:4326\")\n",
    "            gdf_new=gdf_new.append(g1)\n",
    "    gdf_new['stroke']=\"#f71a3a\"\n",
    "    gdf_new['stroke-width']=2\n",
    "    gdf_new['stroke-opacity']=0.7\n",
    "    gdf_new.to_file(archivo_salida_gj, driver=\"GeoJSON\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 120,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "geolife_000_simpl.geojson  ->  geolife_000_linesimpl.geojson\n",
      "geolife_001_simpl.geojson  ->  geolife_001_linesimpl.geojson\n",
      "geolife_002_simpl.geojson  ->  geolife_002_linesimpl.geojson\n",
      "geolife_003_simpl.geojson  ->  geolife_003_linesimpl.geojson\n",
      "geolife_004_simpl.geojson  ->  geolife_004_linesimpl.geojson\n",
      "geolife_005_simpl.geojson  ->  geolife_005_linesimpl.geojson\n",
      "geolife_006_simpl.geojson  ->  geolife_006_linesimpl.geojson\n",
      "geolife_007_simpl.geojson  ->  geolife_007_linesimpl.geojson\n",
      "geolife_008_simpl.geojson  ->  geolife_008_linesimpl.geojson\n",
      "geolife_009_simpl.geojson  ->  geolife_009_linesimpl.geojson\n",
      "geolife_010_simpl.geojson  ->  geolife_010_linesimpl.geojson\n",
      "geolife_011_simpl.geojson  ->  geolife_011_linesimpl.geojson\n",
      "geolife_012_simpl.geojson  ->  geolife_012_linesimpl.geojson\n",
      "geolife_013_simpl.geojson  ->  geolife_013_linesimpl.geojson\n",
      "geolife_014_simpl.geojson  ->  geolife_014_linesimpl.geojson\n",
      "geolife_015_simpl.geojson  ->  geolife_015_linesimpl.geojson\n",
      "geolife_016_simpl.geojson  ->  geolife_016_linesimpl.geojson\n",
      "geolife_017_simpl.geojson  ->  geolife_017_linesimpl.geojson\n",
      "geolife_018_simpl.geojson  ->  geolife_018_linesimpl.geojson\n",
      "geolife_019_simpl.geojson  ->  geolife_019_linesimpl.geojson\n",
      "geolife_020_simpl.geojson  ->  geolife_020_linesimpl.geojson\n",
      "geolife_021_simpl.geojson  ->  geolife_021_linesimpl.geojson\n",
      "geolife_022_simpl.geojson  ->  geolife_022_linesimpl.geojson\n",
      "geolife_023_simpl.geojson  ->  geolife_023_linesimpl.geojson\n",
      "geolife_024_simpl.geojson  ->  geolife_024_linesimpl.geojson\n",
      "geolife_025_simpl.geojson  ->  geolife_025_linesimpl.geojson\n",
      "geolife_026_simpl.geojson  ->  geolife_026_linesimpl.geojson\n",
      "geolife_027_simpl.geojson  ->  geolife_027_linesimpl.geojson\n",
      "geolife_028_simpl.geojson  ->  geolife_028_linesimpl.geojson\n",
      "geolife_029_simpl.geojson  ->  geolife_029_linesimpl.geojson\n",
      "geolife_030_simpl.geojson  ->  geolife_030_linesimpl.geojson\n",
      "geolife_031_simpl.geojson  ->  geolife_031_linesimpl.geojson\n",
      "geolife_032_simpl.geojson  ->  geolife_032_linesimpl.geojson\n",
      "geolife_033_simpl.geojson  ->  geolife_033_linesimpl.geojson\n",
      "geolife_034_simpl.geojson  ->  geolife_034_linesimpl.geojson\n",
      "geolife_035_simpl.geojson  ->  geolife_035_linesimpl.geojson\n",
      "geolife_036_simpl.geojson  ->  geolife_036_linesimpl.geojson\n",
      "geolife_037_simpl.geojson  ->  geolife_037_linesimpl.geojson\n",
      "geolife_038_simpl.geojson  ->  geolife_038_linesimpl.geojson\n",
      "geolife_039_simpl.geojson  ->  geolife_039_linesimpl.geojson\n",
      "geolife_040_simpl.geojson  ->  geolife_040_linesimpl.geojson\n",
      "geolife_041_simpl.geojson  ->  geolife_041_linesimpl.geojson\n",
      "geolife_042_simpl.geojson  ->  geolife_042_linesimpl.geojson\n",
      "geolife_043_simpl.geojson  ->  geolife_043_linesimpl.geojson\n",
      "geolife_044_simpl.geojson  ->  geolife_044_linesimpl.geojson\n",
      "geolife_045_simpl.geojson  ->  geolife_045_linesimpl.geojson\n",
      "geolife_046_simpl.geojson  ->  geolife_046_linesimpl.geojson\n",
      "geolife_047_simpl.geojson  ->  geolife_047_linesimpl.geojson\n",
      "geolife_048_simpl.geojson  ->  geolife_048_linesimpl.geojson\n",
      "geolife_049_simpl.geojson  ->  geolife_049_linesimpl.geojson\n",
      "geolife_050_simpl.geojson  ->  geolife_050_linesimpl.geojson\n",
      "geolife_051_simpl.geojson  ->  geolife_051_linesimpl.geojson\n",
      "geolife_052_simpl.geojson  ->  geolife_052_linesimpl.geojson\n",
      "geolife_053_simpl.geojson  ->  geolife_053_linesimpl.geojson\n",
      "geolife_054_simpl.geojson  ->  geolife_054_linesimpl.geojson\n",
      "geolife_055_simpl.geojson  ->  geolife_055_linesimpl.geojson\n",
      "geolife_056_simpl.geojson  ->  geolife_056_linesimpl.geojson\n",
      "geolife_057_simpl.geojson  ->  geolife_057_linesimpl.geojson\n",
      "geolife_058_simpl.geojson  ->  geolife_058_linesimpl.geojson\n",
      "geolife_059_simpl.geojson  ->  geolife_059_linesimpl.geojson\n",
      "geolife_060_simpl.geojson  ->  geolife_060_linesimpl.geojson\n",
      "geolife_061_simpl.geojson  ->  geolife_061_linesimpl.geojson\n",
      "geolife_062_simpl.geojson  ->  geolife_062_linesimpl.geojson\n",
      "geolife_063_simpl.geojson  ->  geolife_063_linesimpl.geojson\n",
      "geolife_064_simpl.geojson  ->  geolife_064_linesimpl.geojson\n",
      "geolife_065_simpl.geojson  ->  geolife_065_linesimpl.geojson\n",
      "geolife_066_simpl.geojson  ->  geolife_066_linesimpl.geojson\n",
      "geolife_067_simpl.geojson  ->  geolife_067_linesimpl.geojson\n",
      "geolife_068_simpl.geojson  ->  geolife_068_linesimpl.geojson\n",
      "geolife_069_simpl.geojson  ->  geolife_069_linesimpl.geojson\n",
      "geolife_070_simpl.geojson  ->  geolife_070_linesimpl.geojson\n",
      "geolife_071_simpl.geojson  ->  geolife_071_linesimpl.geojson\n",
      "geolife_072_simpl.geojson  ->  geolife_072_linesimpl.geojson\n",
      "geolife_073_simpl.geojson  ->  geolife_073_linesimpl.geojson\n",
      "geolife_074_simpl.geojson  ->  geolife_074_linesimpl.geojson\n",
      "geolife_075_simpl.geojson  ->  geolife_075_linesimpl.geojson\n",
      "geolife_076_simpl.geojson  ->  geolife_076_linesimpl.geojson\n",
      "geolife_077_simpl.geojson  ->  geolife_077_linesimpl.geojson\n",
      "geolife_078_simpl.geojson  ->  geolife_078_linesimpl.geojson\n",
      "geolife_079_simpl.geojson  ->  geolife_079_linesimpl.geojson\n",
      "geolife_080_simpl.geojson  ->  geolife_080_linesimpl.geojson\n",
      "geolife_081_simpl.geojson  ->  geolife_081_linesimpl.geojson\n",
      "geolife_082_simpl.geojson  ->  geolife_082_linesimpl.geojson\n",
      "geolife_083_simpl.geojson  ->  geolife_083_linesimpl.geojson\n",
      "geolife_084_simpl.geojson  ->  geolife_084_linesimpl.geojson\n",
      "geolife_085_simpl.geojson  ->  geolife_085_linesimpl.geojson\n",
      "geolife_086_simpl.geojson  ->  geolife_086_linesimpl.geojson\n",
      "geolife_087_simpl.geojson  ->  geolife_087_linesimpl.geojson\n",
      "geolife_088_simpl.geojson  ->  geolife_088_linesimpl.geojson\n",
      "geolife_089_simpl.geojson  ->  geolife_089_linesimpl.geojson\n",
      "geolife_090_simpl.geojson  ->  geolife_090_linesimpl.geojson\n",
      "geolife_091_simpl.geojson  ->  geolife_091_linesimpl.geojson\n",
      "geolife_092_simpl.geojson  ->  geolife_092_linesimpl.geojson\n",
      "geolife_093_simpl.geojson  ->  geolife_093_linesimpl.geojson\n",
      "geolife_094_simpl.geojson  ->  geolife_094_linesimpl.geojson\n",
      "geolife_095_simpl.geojson  ->  geolife_095_linesimpl.geojson\n",
      "geolife_096_simpl.geojson  ->  geolife_096_linesimpl.geojson\n",
      "geolife_097_simpl.geojson  ->  geolife_097_linesimpl.geojson\n",
      "geolife_098_simpl.geojson  ->  geolife_098_linesimpl.geojson\n",
      "geolife_099_simpl.geojson  ->  geolife_099_linesimpl.geojson\n",
      "geolife_100_simpl.geojson  ->  geolife_100_linesimpl.geojson\n",
      "geolife_101_simpl.geojson  ->  geolife_101_linesimpl.geojson\n",
      "geolife_102_simpl.geojson  ->  geolife_102_linesimpl.geojson\n",
      "geolife_103_simpl.geojson  ->  geolife_103_linesimpl.geojson\n",
      "geolife_104_simpl.geojson  ->  geolife_104_linesimpl.geojson\n",
      "geolife_105_simpl.geojson  ->  geolife_105_linesimpl.geojson\n",
      "geolife_106_simpl.geojson  ->  geolife_106_linesimpl.geojson\n",
      "geolife_107_simpl.geojson  ->  geolife_107_linesimpl.geojson\n",
      "geolife_108_simpl.geojson  ->  geolife_108_linesimpl.geojson\n",
      "geolife_109_simpl.geojson  ->  geolife_109_linesimpl.geojson\n",
      "geolife_110_simpl.geojson  ->  geolife_110_linesimpl.geojson\n",
      "geolife_111_simpl.geojson  ->  geolife_111_linesimpl.geojson\n",
      "geolife_112_simpl.geojson  ->  geolife_112_linesimpl.geojson\n",
      "geolife_113_simpl.geojson  ->  geolife_113_linesimpl.geojson\n",
      "geolife_114_simpl.geojson  ->  geolife_114_linesimpl.geojson\n",
      "geolife_115_simpl.geojson  ->  geolife_115_linesimpl.geojson\n",
      "geolife_116_simpl.geojson  ->  geolife_116_linesimpl.geojson\n",
      "geolife_117_simpl.geojson  ->  geolife_117_linesimpl.geojson\n",
      "geolife_118_simpl.geojson  ->  geolife_118_linesimpl.geojson\n",
      "geolife_119_simpl.geojson  ->  geolife_119_linesimpl.geojson\n",
      "geolife_120_simpl.geojson  ->  geolife_120_linesimpl.geojson\n",
      "geolife_121_simpl.geojson  ->  geolife_121_linesimpl.geojson\n",
      "geolife_122_simpl.geojson  ->  geolife_122_linesimpl.geojson\n",
      "geolife_123_simpl.geojson  ->  geolife_123_linesimpl.geojson\n",
      "geolife_124_simpl.geojson  ->  geolife_124_linesimpl.geojson\n",
      "geolife_125_simpl.geojson  ->  geolife_125_linesimpl.geojson\n",
      "geolife_126_simpl.geojson  ->  geolife_126_linesimpl.geojson\n",
      "geolife_127_simpl.geojson  ->  geolife_127_linesimpl.geojson\n",
      "geolife_128_simpl.geojson  ->  geolife_128_linesimpl.geojson\n",
      "geolife_129_simpl.geojson  ->  geolife_129_linesimpl.geojson\n",
      "geolife_130_simpl.geojson  ->  geolife_130_linesimpl.geojson\n",
      "geolife_131_simpl.geojson  ->  geolife_131_linesimpl.geojson\n",
      "geolife_132_simpl.geojson  ->  geolife_132_linesimpl.geojson\n",
      "geolife_133_simpl.geojson  ->  geolife_133_linesimpl.geojson\n",
      "geolife_134_simpl.geojson  ->  geolife_134_linesimpl.geojson\n",
      "geolife_135_simpl.geojson  ->  geolife_135_linesimpl.geojson\n",
      "geolife_136_simpl.geojson  ->  geolife_136_linesimpl.geojson\n",
      "geolife_137_simpl.geojson  ->  geolife_137_linesimpl.geojson\n",
      "geolife_138_simpl.geojson  ->  geolife_138_linesimpl.geojson\n",
      "geolife_139_simpl.geojson  ->  geolife_139_linesimpl.geojson\n",
      "geolife_140_simpl.geojson  ->  geolife_140_linesimpl.geojson\n",
      "geolife_141_simpl.geojson  ->  geolife_141_linesimpl.geojson\n",
      "geolife_142_simpl.geojson  ->  geolife_142_linesimpl.geojson\n",
      "geolife_143_simpl.geojson  ->  geolife_143_linesimpl.geojson\n",
      "geolife_144_simpl.geojson  ->  geolife_144_linesimpl.geojson\n",
      "geolife_145_simpl.geojson  ->  geolife_145_linesimpl.geojson\n",
      "geolife_146_simpl.geojson  ->  geolife_146_linesimpl.geojson\n",
      "geolife_147_simpl.geojson  ->  geolife_147_linesimpl.geojson\n",
      "geolife_148_simpl.geojson  ->  geolife_148_linesimpl.geojson\n",
      "geolife_149_simpl.geojson  ->  geolife_149_linesimpl.geojson\n",
      "geolife_150_simpl.geojson  ->  geolife_150_linesimpl.geojson\n",
      "geolife_151_simpl.geojson  ->  geolife_151_linesimpl.geojson\n",
      "geolife_152_simpl.geojson  ->  geolife_152_linesimpl.geojson\n",
      "geolife_153_simpl.geojson  ->  geolife_153_linesimpl.geojson\n",
      "geolife_154_simpl.geojson  ->  geolife_154_linesimpl.geojson\n",
      "geolife_155_simpl.geojson  ->  geolife_155_linesimpl.geojson\n",
      "geolife_156_simpl.geojson  ->  geolife_156_linesimpl.geojson\n",
      "geolife_157_simpl.geojson  ->  geolife_157_linesimpl.geojson\n",
      "geolife_158_simpl.geojson  ->  geolife_158_linesimpl.geojson\n",
      "geolife_159_simpl.geojson  ->  geolife_159_linesimpl.geojson\n",
      "geolife_160_simpl.geojson  ->  geolife_160_linesimpl.geojson\n",
      "geolife_161_simpl.geojson  ->  geolife_161_linesimpl.geojson\n",
      "geolife_162_simpl.geojson  ->  geolife_162_linesimpl.geojson\n",
      "geolife_163_simpl.geojson  ->  geolife_163_linesimpl.geojson\n",
      "geolife_164_simpl.geojson  ->  geolife_164_linesimpl.geojson\n",
      "geolife_165_simpl.geojson  ->  geolife_165_linesimpl.geojson\n",
      "geolife_166_simpl.geojson  ->  geolife_166_linesimpl.geojson\n",
      "geolife_167_simpl.geojson  ->  geolife_167_linesimpl.geojson\n",
      "geolife_168_simpl.geojson  ->  geolife_168_linesimpl.geojson\n",
      "geolife_169_simpl.geojson  ->  geolife_169_linesimpl.geojson\n",
      "geolife_170_simpl.geojson  ->  geolife_170_linesimpl.geojson\n",
      "geolife_171_simpl.geojson  ->  geolife_171_linesimpl.geojson\n",
      "geolife_172_simpl.geojson  ->  geolife_172_linesimpl.geojson\n",
      "geolife_173_simpl.geojson  ->  geolife_173_linesimpl.geojson\n",
      "geolife_174_simpl.geojson  ->  geolife_174_linesimpl.geojson\n",
      "geolife_175_simpl.geojson  ->  geolife_175_linesimpl.geojson\n",
      "geolife_176_simpl.geojson  ->  geolife_176_linesimpl.geojson\n",
      "geolife_177_simpl.geojson  ->  geolife_177_linesimpl.geojson\n",
      "geolife_178_simpl.geojson  ->  geolife_178_linesimpl.geojson\n",
      "geolife_179_simpl.geojson  ->  geolife_179_linesimpl.geojson\n",
      "geolife_180_simpl.geojson  ->  geolife_180_linesimpl.geojson\n",
      "geolife_181_simpl.geojson  ->  geolife_181_linesimpl.geojson\n"
     ]
    }
   ],
   "source": [
    "for i in range(0,182):\n",
    "    simpl2lines(i)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### def **simpl2lines2():** \n",
    ">### le agrega un marcador donde empieza y donde termina la trayectoria"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 122,
   "metadata": {},
   "outputs": [],
   "source": [
    "def simpl2lines2(xxx):\n",
    "    archivo_gj=f\"geolife_{xxx:03d}_\"\n",
    "    archivo_salida_gj = archivo_gj+\"linesimpl2.geojson\"\n",
    "    archivo_gj=archivo_gj+\"simpl.geojson\"\n",
    "    print(archivo_gj,\" -> \" ,archivo_salida_gj)\n",
    "    gdf = geopandas.read_file(archivo_gj)\n",
    "    \n",
    "    gempty= {'date': [], 'geometry': []}\n",
    "    gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")    \n",
    "    \n",
    "    list_dates=list(gdf.date.str[:10].unique())\n",
    "    gempty= {'date': [], 'geometry': []}\n",
    "    gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")\n",
    "    for dia in list_dates:\n",
    "        plst=gdf.loc[gdf.date.str[:10]==dia,'geometry'].to_list()\n",
    "        if len(plst)<2:\n",
    "            continue\n",
    "        else:\n",
    "            ls=LineString(plst)\n",
    "            grow = {'date': [dia], 'geometry': [ls]}\n",
    "            g1=geopandas.GeoDataFrame(grow, crs=\"EPSG:4326\")\n",
    "            gdf_new=gdf_new.append(g1)\n",
    "    gdf_new['stroke']=\"#f71a3a\"\n",
    "    gdf_new['stroke-width']=2\n",
    "    gdf_new['stroke-opacity']=0.7\n",
    "    gdf_new.to_file(archivo_salida_gj, driver=\"GeoJSON\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(0,182):\n",
    "    simpl2lines(i)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## _Ahora generamos archivos geojson taggeados con el **label** de **geolife**_ (mapeamos cada medio de  transporte a un color e icono)\n",
    ">## **geolife_xxx_colorlabel.geojson**"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    ">### Primero definimos etiquetas, colores, cmaps para taggear los puntos geojson con colores e iconos "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 168,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'#fde725'"
      ]
     },
     "execution_count": 168,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mode_names = ['walk', 'bike', 'bus', 'car', 'subway','train', \n",
    "              'airplane', 'boat', 'run', 'motorcycle', 'taxi']\n",
    "mode_marker = ['pitch', 'bicycle', 'bus', \"car\", 'rail-underground','rail', \n",
    "              'airport', 'ferry', 'pitch', 'scooter', 'car']\n",
    "mode_color = ['#f16f04', '#1ac16b', '#68aa31', '#8043b1', '#29c6db','#470e53', \n",
    "              '#595454', '#2026d5', '#f16f04', '#e52315', '#f6fa00']\n",
    "\n",
    "mode_ids = {s : i + 1 for i, s in enumerate(mode_names)}\n",
    "mode_marker_dict = {i + 1:s for i, s in enumerate(mode_marker)}\n",
    "mode_color_dict = {i + 1:s for i, s in enumerate(mode_color)}\n",
    "\n",
    "cmaps = {}\n",
    "cmaps['Uniform'] = [ 'viridis', 'plasma', 'inferno', 'magma', 'cividis','spring','summer','autumn','winter']\n",
    "cmaps['Diverging'] = [ 'PiYG', 'PRGn', 'BrBG', 'PuOr', 'RdGy', 'RdBu',\n",
    "            'RdYlBu', 'RdYlGn', 'Spectral', 'coolwarm', 'bwr', 'seismic']\n",
    "cmaps['Cyclic'] = ['twilight', 'twilight_shifted', 'hsv']\n",
    "cmaps['Miscellaneous'] = [ 'flag', 'prism', 'ocean', 'gist_earth', 'terrain', 'gist_stern',\n",
    "            'gnuplot', 'gnuplot2', 'CMRmap', 'cubehelix', 'brg',\n",
    "            'gist_rainbow', 'rainbow', 'jet', 'nipy_spectral', 'gist_ncar']\n",
    "\n",
    "cmap = plt.get_cmap(cmaps['Uniform'][0])\n",
    "colors.to_hex(cmap(255))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## def **agregaColorLabel(**xxx**):** \n",
    ">### Genera **geolife_xxx_colorlabel.geojson** (cuando los labels son != 0) y **geolife_xxx_colorlabel0.geojson** (que incluye la data del label 0) a partir de **geolife_xxx_simpl.geojson** "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 184,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4\n"
     ]
    }
   ],
   "source": [
    "def agregaColorLabelSimpl(xxx):\n",
    "\n",
    "    #def simpl2lines(xxx):\n",
    "    archivo_gj=f\"geolife_{xxx:03d}_\"\n",
    "    archivo_salida_gj = archivo_gj+\"colorlabel.geojson\"\n",
    "    archivo_gj=archivo_gj+\"simpl.geojson\"\n",
    "    print(archivo_gj,\" -> \" ,archivo_salida_gj)\n",
    "    dftmp = geopandas.read_file(archivo_gj)\n",
    "    \n",
    "    if len(list(dftmp.label.unique()))=1\n",
    "        print(f\"No hay labels para {xxx:03d}\")\n",
    "        return False\n",
    "    \n",
    "    ### Crea rangos de labels\n",
    "    print(len(list(dftmp.label.unique())))\n",
    "    label_col = dftmp.columns.get_loc('label')\n",
    "    labeltmp=0\n",
    "    i=0\n",
    "    start_ix = i\n",
    "    label_range_tmp = [] # [start_ix, end_ix, label]\n",
    "    while i < len(dftmp):\n",
    "        runlabel=dftmp.iloc[i,label_col]\n",
    "        if labeltmp != runlabel:\n",
    "            end_ix = i-1\n",
    "            if end_ix > 0:\n",
    "                label_range_tmp.append((start_ix,end_ix,labeltmp))\n",
    "            start_ix = i\n",
    "        labeltmp = runlabel\n",
    "        i+=1\n",
    "    label_range = [ a for a in label_range_tmp if a[2]!='0' ]\n",
    "    label_range_0 = [ a for a in label_range_tmp if a[2]='0' ]\n",
    "    \n",
    "    gdf = dftmp\n",
    "    \n",
    "    gempty= {'date': [], 'geometry': []}\n",
    "    gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")    \n",
    "    \n",
    "    tramo0 = label_range[0]\n",
    "    \n",
    "    p1=gdf.iloc[,gdf.columns.get_loc('geometry')]\n",
    "    gempty= {'date': [dia1],\"marker-size\":[\"large\"],\"marker-color\": [\"#d77070\"],\"marker-symbol\":['entrance'], 'geometry': [p1]}\n",
    "    gdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")    \n",
    "\n",
    "    \n",
    "    for tramo in label_range:\n",
    "        plst=gdf.iloc[tramo[0]:tramo[1],'geometry'].to_list()\n",
    "        if len(plst)<2:\n",
    "            continue\n",
    "        else:\n",
    "            ls=LineString(plst)\n",
    "            grow = {'date': [dia], 'geometry': [ls]}\n",
    "            g1=geopandas.GeoDataFrame(grow, crs=\"EPSG:4326\")\n",
    "            gdf_new=gdf_new.append(g1)\n",
    "    gdf_new['stroke']=\"#f71a3a\"\n",
    "    gdf_new['stroke-width']=2\n",
    "    gdf_new['stroke-opacity']=0.7\n",
    "    gdf_new.to_file(archivo_salida_gj, driver=\"GeoJSON\")\n",
    "\n",
    "        \n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 191,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "7"
      ]
     },
     "execution_count": 191,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "agregaColorLabelSimpl(161)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 193,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>marker-color</th>\n",
       "      <th>marker-size</th>\n",
       "      <th>marker-symbol</th>\n",
       "      <th>quees</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>#7e7e7e</td>\n",
       "      <td>large</td>\n",
       "      <td>rail</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (112.14844 35.31737)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>#8043b1</td>\n",
       "      <td>large</td>\n",
       "      <td>car</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.15041 35.53167)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>LINESTRING (118.15006 35.53195, 118.15727 35.5...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>#595454</td>\n",
       "      <td>large</td>\n",
       "      <td>airport</td>\n",
       "      <td>airplane</td>\n",
       "      <td>POINT (118.13358 35.50792)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>#68aa31</td>\n",
       "      <td>large</td>\n",
       "      <td>bus</td>\n",
       "      <td>bus</td>\n",
       "      <td>POINT (118.26302 35.55793)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>#470e53</td>\n",
       "      <td>large</td>\n",
       "      <td>rail</td>\n",
       "      <td>train</td>\n",
       "      <td>POINT (118.31795 35.57524)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>#c14444</td>\n",
       "      <td>large</td>\n",
       "      <td>parking</td>\n",
       "      <td>peaton</td>\n",
       "      <td>POINT (118.29460 35.49422)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>#1ac16b</td>\n",
       "      <td>large</td>\n",
       "      <td>bicycle</td>\n",
       "      <td>bike</td>\n",
       "      <td>POINT (118.38181 35.54508)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>#f16f04</td>\n",
       "      <td>large</td>\n",
       "      <td>pitch</td>\n",
       "      <td>walk</td>\n",
       "      <td>POINT (118.36121 35.46570)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.29597 35.45004)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>#2026d5</td>\n",
       "      <td>large</td>\n",
       "      <td>ferry</td>\n",
       "      <td>boat</td>\n",
       "      <td>POINT (118.14972 35.38233)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>#e52315</td>\n",
       "      <td>large</td>\n",
       "      <td>scooter</td>\n",
       "      <td>motorcycle</td>\n",
       "      <td>POINT (118.16895 35.62605)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>#f6fa00</td>\n",
       "      <td>large</td>\n",
       "      <td>car</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.25546 35.62270)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>#29c6db</td>\n",
       "      <td>large</td>\n",
       "      <td>rail-underground</td>\n",
       "      <td>subway</td>\n",
       "      <td>POINT (118.26370 35.70861)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>#7e7e7e</td>\n",
       "      <td>large</td>\n",
       "      <td>roadblock</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.64685 35.66622)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>LINESTRING (118.64685 35.67515, 118.87207 35.7...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>#7e7e7e</td>\n",
       "      <td>large</td>\n",
       "      <td>entrance</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.61938 35.23889)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>#7e7e7e</td>\n",
       "      <td>large</td>\n",
       "      <td>rail</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (112.14844 35.31737)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>#8043b1</td>\n",
       "      <td>large</td>\n",
       "      <td>car</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.15041 35.53167)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>None</td>\n",
       "      <td>LINESTRING (118.15006 35.53195, 118.15727 35.5...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>#595454</td>\n",
       "      <td>large</td>\n",
       "      <td>airport</td>\n",
       "      <td>airplane</td>\n",
       "      <td>POINT (118.13358 35.50792)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>#68aa31</td>\n",
       "      <td>large</td>\n",
       "      <td>bus</td>\n",
       "      <td>bus</td>\n",
       "      <td>POINT (118.26302 35.55793)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>#470e53</td>\n",
       "      <td>large</td>\n",
       "      <td>rail</td>\n",
       "      <td>train</td>\n",
       "      <td>POINT (118.31795 35.57524)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>#c14444</td>\n",
       "      <td>large</td>\n",
       "      <td>parking</td>\n",
       "      <td>peaton</td>\n",
       "      <td>POINT (118.29460 35.49422)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>#1ac16b</td>\n",
       "      <td>large</td>\n",
       "      <td>bicycle</td>\n",
       "      <td>bike</td>\n",
       "      <td>POINT (118.38181 35.54508)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>#f16f04</td>\n",
       "      <td>large</td>\n",
       "      <td>pitch</td>\n",
       "      <td>walk</td>\n",
       "      <td>POINT (118.36121 35.46570)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>#2026d5</td>\n",
       "      <td>large</td>\n",
       "      <td>ferry</td>\n",
       "      <td>boat</td>\n",
       "      <td>POINT (118.14972 35.38233)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>#e52315</td>\n",
       "      <td>large</td>\n",
       "      <td>scooter</td>\n",
       "      <td>motorcycle</td>\n",
       "      <td>POINT (118.16895 35.62605)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>#f6fa00</td>\n",
       "      <td>large</td>\n",
       "      <td>car</td>\n",
       "      <td>None</td>\n",
       "      <td>POINT (118.25546 35.62270)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>#29c6db</td>\n",
       "      <td>large</td>\n",
       "      <td>rail-underground</td>\n",
       "      <td>subway</td>\n",
       "      <td>POINT (118.26370 35.70861)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   marker-color marker-size     marker-symbol       quees  \\\n",
       "0       #7e7e7e       large              rail        None   \n",
       "1       #8043b1       large               car        None   \n",
       "2          None        None              None        None   \n",
       "3       #595454       large           airport    airplane   \n",
       "4       #68aa31       large               bus         bus   \n",
       "5       #470e53       large              rail       train   \n",
       "6       #c14444       large           parking      peaton   \n",
       "7       #1ac16b       large           bicycle        bike   \n",
       "8       #f16f04       large             pitch        walk   \n",
       "9          None        None              None        None   \n",
       "10      #2026d5       large             ferry        boat   \n",
       "11      #e52315       large           scooter  motorcycle   \n",
       "12      #f6fa00       large               car        None   \n",
       "13      #29c6db       large  rail-underground      subway   \n",
       "14      #7e7e7e       large         roadblock        None   \n",
       "15         None        None              None        None   \n",
       "16      #7e7e7e       large          entrance        None   \n",
       "17      #7e7e7e       large              rail        None   \n",
       "18      #8043b1       large               car        None   \n",
       "19         None        None              None        None   \n",
       "20      #595454       large           airport    airplane   \n",
       "21      #68aa31       large               bus         bus   \n",
       "22      #470e53       large              rail       train   \n",
       "23      #c14444       large           parking      peaton   \n",
       "24      #1ac16b       large           bicycle        bike   \n",
       "25      #f16f04       large             pitch        walk   \n",
       "26      #2026d5       large             ferry        boat   \n",
       "27      #e52315       large           scooter  motorcycle   \n",
       "28      #f6fa00       large               car        None   \n",
       "29      #29c6db       large  rail-underground      subway   \n",
       "\n",
       "                                             geometry  \n",
       "0                          POINT (112.14844 35.31737)  \n",
       "1                          POINT (118.15041 35.53167)  \n",
       "2   LINESTRING (118.15006 35.53195, 118.15727 35.5...  \n",
       "3                          POINT (118.13358 35.50792)  \n",
       "4                          POINT (118.26302 35.55793)  \n",
       "5                          POINT (118.31795 35.57524)  \n",
       "6                          POINT (118.29460 35.49422)  \n",
       "7                          POINT (118.38181 35.54508)  \n",
       "8                          POINT (118.36121 35.46570)  \n",
       "9                          POINT (118.29597 35.45004)  \n",
       "10                         POINT (118.14972 35.38233)  \n",
       "11                         POINT (118.16895 35.62605)  \n",
       "12                         POINT (118.25546 35.62270)  \n",
       "13                         POINT (118.26370 35.70861)  \n",
       "14                         POINT (118.64685 35.66622)  \n",
       "15  LINESTRING (118.64685 35.67515, 118.87207 35.7...  \n",
       "16                         POINT (118.61938 35.23889)  \n",
       "17                         POINT (112.14844 35.31737)  \n",
       "18                         POINT (118.15041 35.53167)  \n",
       "19  LINESTRING (118.15006 35.53195, 118.15727 35.5...  \n",
       "20                         POINT (118.13358 35.50792)  \n",
       "21                         POINT (118.26302 35.55793)  \n",
       "22                         POINT (118.31795 35.57524)  \n",
       "23                         POINT (118.29460 35.49422)  \n",
       "24                         POINT (118.38181 35.54508)  \n",
       "25                         POINT (118.36121 35.46570)  \n",
       "26                         POINT (118.14972 35.38233)  \n",
       "27                         POINT (118.16895 35.62605)  \n",
       "28                         POINT (118.25546 35.62270)  \n",
       "29                         POINT (118.26370 35.70861)  "
      ]
     },
     "execution_count": 193,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "g=geopandas.read_file(\"demoicons.geojson\")\n",
    "g"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 180,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "geolife_052_simpl.geojson  ->  geolife_052_colorlabel.geojson\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "' \\nlist_labels=list(gdf.date.str[:10].unique())\\ngempty= {\\'date\\': [], \\'geometry\\': []}\\ngdf_new=geopandas.GeoDataFrame(gempty, crs=\"EPSG:4326\")\\nfor dia in list_dates:\\n    plst=gdf.loc[gdf.date.str[:10]==dia,\\'geometry\\'].to_list()\\n    if len(plst)<2:\\n        continue\\n    else:\\n        ls=LineString(plst)\\n        grow = {\\'date\\': [dia], \\'geometry\\': [ls]}\\n        g1=geopandas.GeoDataFrame(grow, crs=\"EPSG:4326\")\\n        gdf_new=gdf_new.append(g1)\\ngdf_new[\\'stroke\\']=\"#f71a3a\"\\ngdf_new[\\'stroke-width\\']=2\\ngdf_new[\\'stroke-opacity\\']=0.7\\ngdf_new.to_file(archivo_salida_gj, driver=\"GeoJSON\")\\n'"
      ]
     },
     "execution_count": 180,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xxx=52\n",
    "#def simpl2lines(xxx):\n",
    "archivo_gj=f\"geolife_{xxx:03d}_\"\n",
    "archivo_salida_gj = archivo_gj+\"colorlabel.geojson\"\n",
    "archivo_gj=archivo_gj+\"simpl.geojson\"\n",
    "print(archivo_gj,\" -> \" ,archivo_salida_gj)\n",
    "gdf = geopandas.read_file(archivo_gj)\n",
    "gdf.head(6)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Distribución de labels"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 183,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABJEAAAJSCAYAAABz4lSSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdfZRk6V0f9u+ve3qcboHpAQ2gGbEziCiKe3ISwU7Mmw/HpgcsK3aEYzgWESAcfAbPQIwcOzEvCW/2+oBxHBufSIkMGJEdI2Pe7RCwdgET+yBZswKDdmUiYXYWrQQaLMnI3jWrnX3yR1dLPT3d/VR1V9Wtnvl8zrlnuu7r77489976TtWtaq0FAAAAAA6yNHQBAAAAACw+IRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANB1YugCDuv5z39+O3/+/NBlAAAAANw1Hnnkkd9prZ3ea9ixDZHOnz+f69evD10GAAAAwF2jqm7sN8zX2QAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxCJrK2tpao+0q2trQ1d0oF21rrd3a32Wte7eX1h1oY43+3XjnvduPNZdBcuXLij5iGvM7uPgd3dqVOnkiSnTp2aSt1Xr17dcznLy8tjTb9fneNOP4lpHV9Hmc9+22tW67xIlpeXp7a+J0+evGP7nTx5corVAovkON4f7HTc67/XCJHucWtra3n66adv6/f0008vbJC03wnlbjzRHLROd+P6wqwNcb47SlvdOe1xPR9cuHAhjz322B39h7rO7HUM7PbBD34wVZUPfvCDdwybtO6rV6/mda973Z7DnnvuuW5IcNC+HWf6SUzr+nqU+Ry0vZLpr/MiWV5eznPPPXdbv8Ou78mTJ/PhD3/4jv4f/vCHBUlwFzru74+Oe/33omqtDV3DoVy8eLFdv3596DKOvYMa5yIeG8et3qPonTjvtvWFWRvi/HHUG6Dtuo7r+WDR6p7WDem4dZ84cSK3bt069LzGqXda23Ba7eMo8xlne01az3ExzfPTorU7YLaO+/uj417/3aqqHmmtXdxrmE8iAQDMwDiBCB9lewHA4hMiAQDMwN361atZsb0AYPEJke5xq6urE/UHOK6c7+ZvY2Nj32FDbPdpLHOSeVy+fPnA4UtLR7sNO+r0i6a3vZK7b5237bdeh1nflZWVQw0DgHHcnVdixvbUU0/dcUO8urqap556aqCKDrbf92Lvxu/LHrROd+P6wqwNcb47SlvdOe1xPR88+uijewZJQ11n9joGdltfX09rLevr63cMm7Tu1772tbly5cqew5aWlo70vKRxpp/EtK6vR5nPQdsrmf46L5Jbt27dERgddn2feeaZPcOilZWVPPPMM4euEVhMx/390XGv/17kwdoAAAAAJPFgbQAAAACOSIgEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXXMNkarq+6rqfVX19h39Pr6q3lRV7xz9e2qeNQEAAADQN+9PIn1/kpft6vf1SR5urb04ycOj1wAAAAAskLmGSK21X0jy/l29X5HkDaO/35Dki+ZZEwAAAAB9i/BMpE9qrb03SUb/fuLA9QAAAACwyyKESGOrqstVdb2qrt+8eXPocgAAAADuGYsQIv12Vb0gSUb/vm+/EVtrr2+tXWytXTx9+vTcCgQAAAC41y1CiPSTSV49+vvVSX5iwFoAAAAA2MNcQ6Sq+sEkv5jkJVX17qr6qiTfkeQLquqdSb5g9BoAAACABXJingtrrX3pPoM251kHAAAAAJNZhK+zAQAAALDghEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6FiZEqqq/UFWPVtXbq+oHq+o/GromAAAAALYsRIhUVWeT/PkkF1tr/1mS5SSvHLYqAAAAALYtRIg0ciLJalWdSLKW5D0D1wMAAADAyEKESK21J5P8jSRPJHlvkn/bWvsnw1YFAAAAwLaFCJGq6lSSVyT51CRnkjyvqr5sj/EuV9X1qrp+8+bNeZcJAAAAcM9aiBApyaUkv9Fau9la+3CSH03yObtHaq29vrV2sbV28fTp03MvEgAAAOBetSgh0hNJPquq1qqqkmwmecfANQEAAAAwshAhUmvtLUl+OMnbkvxqtup6/aBFAQAAAPARJ4YuYFtr7VuSfMvQdQAAAABwp4X4JBIAAAAAi02IBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAsICqas/uIGtraxNPA7AoLl265BwGsOCESACwYA5607TfsLW1tTz99NMTzw9gEVy6dCkPP/zwnsOcwwAWhxAJAO4C+wVIAMfBfgESAItFiAQAAABAlxAJAAAAgC4hEgDcBVZXV4cuAeDQNjc3hy4BgDEIkQBgwbTWJh721FNP7RskHTQ/gEXw0EMP7RskOYcBLI4TQxcAANzpMG+annrqqRlUAjAfDz300NAlANDhk0gAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoGthQqSqWq+qH66qf1VV76iqzx66JgAAAAC2nBi6gB3+dpKfbq19cVWdTLI2dEEAAAAAbFmIEKmqfn+Sz0vylUnSWnsmyTND1gQAAADARy3K19lelORmkr9XVb9UVd9TVc8buigAAAAAtkwtRKqqS1X1K4ec/ESSz0jyutbapyf590m+fo9lXK6q61V1/ebNm0eoFgAAAIBJTPOTSB+X5MIhp313kne31t4yev3D2QqVbtNae31r7WJr7eLp06cPuSgAAAAAJrUQX2drrf1Wkt+sqpeMem0meWzAkgAAAADYYSEerD3y3ye5Nvpltn+d5M8MXA8AAAAAIwsTIrXWfjnJxaHrAAAAAOBO3RCpql405rw++Yi1AAAAALCgxvkk0ruStDHGqzHHAwAAAOCYGSdE8mwiAAAAgHtcN0Rqrb1hHoUAAAAAsLiWhi4AAAAAgMXXDZGq6m1VdWHH66qq766qF+4a72JV/e4sigQAAABgWON8EumlSZ63a5qvSfKJu8Zb3jUeAAAAAHeJw36draZaBQAAAAALzTORAAAAAOgSIgEAAADQdWLM8ZaqajtwWt6j387+AAAAANxlxg2R/vke/d4yzUIAAAAAWFzjhEjfNvMqAAAAAFho3RCptSZEAgAAALjHebA2AAAAAF3jPhMpVfWSJCuttbePXp9O8j8k+c+TvC/J/9lae/NMqgQAAABgUN0QqarOJPnHSf6L0eu3J/miJG9K8qlJfifJ85P8t1X1Oa21R2ZXLgAAAABDGOfrbN+e5AVJLif50iTPZStUen+SF7TWPinJC5O8Pck3z6hOAAAAAAY0ztfZLiX55tba9yZJVf16krcm+ZLW2vuSpLX23qr6G0m+a2aVAgAAADCYcT6JdDbJYztePzr698au8R5P8slTqAkAAACABTNOiLSc5MM7Xm///dyu8Z5LUtMoCgAAAIDFMu6vsy1V1XbgtLxHv539AQAAALjLjBsi/fM9+r1lmoUAAAAAsLjGCZG+beZVAAAAALDQuiFSa02IBAAAAHCPG+fB2mOpqvNV9c3Tmh8AAAAAi+NIIVJVfUxV/XdV9fNJ3pXkW6ZSFQAAAAALZeIQqbZ8YVU9mOS3kvzdJP9Jku9K8p9OuT4AAAAAFsC4v86WqtpI8hVJvizJC5I8k+RnkvyJJK9srf3CTCoEAAAAYHDdEKmqvjbJq5N8RpJK8otJvj3JPxi9fv8sCwQAAABgeON8Eum7k7QkP5XkNa21X98eUFUfN6vCAAAAAFgc4zwT6WezFSK9PMmPV9VfqqoXzLYsAAAAABZJN0RqrV1Kcj7J/5KtTy799SRPVNVPJ/nSbAVMAAAAANzFxvp1ttbau1trf6219geSfE6S70nyXyZ57WiU11TVH5pRjQAAAAAMbKwQaafW2ptba1ey9Qttfzpbz0r640n+aVX92pTrAwAAAGABTBwibWutPdNa+4ettT+R5IVJvjHJM1OrDAAAAICFcegQaafW2vuSvCvJxjTmx2K5du1azp8/n6WlpZw/fz7Xrl0buiSAhXHp0qVU1Ue6tbW1215fuHBh6BIXztmzZ2/bRmfPnh26pD0d9fp39erV29ZzqONh9zFZVbl06dLc69i2u5aqGmsYk3H/BtOxKOelvepYlNq2TXre2X0PNeS1ifFNJUTi7nXt2rVcvnw5N27cSGstN27cyOXLl92IAGTr5ufhhx++rd/TTz992+vHHntMkLTD2bNn8573vOe2fu95z3sWLkg66vXv6tWred3rXndH/3kfD2tra3cck0ny8MMPD3Kzvt8bnIPe/Az9pug4cv8G07Eo56VJljfUOXPS885e91BDXZuYTLU2nR9Xq6o/leSHWmvLU5lhx8WLF9v169fnsah72vnz53Pjxo07+p87dy6PP/74/AsCWCCT3KhN63p73B20zRZpGx31+nfixIncunVr3+HzWtfeMTrvbX7YNzeLdGwcB+7fYDoW5Zo16blziHPmpOedRdm27K2qHmmtXdxzmBCJgywtLe3ZiKsqzz333AAVASwOIdLkjstN41Gvf4sS3ixKHduESPPh/g2mY1GuWcchRJr0vLMo25a9HRQidb/OVlUvGqdL8slTr5zB3XfffRP1B4C7wVGvf8vLc/k/NdiT+zdg3px37h3jPBPpXUneOUb3d2ZUIwN64IEHsra2dlu/tbW1PPDAAwNVBLA4Njc3xxpvY8PvTmw7c+bMRP2HctTr3+XLl/cdNs/jYXV1dd9h4x6/HD/u34B5m/S8s981yLVp8XW/zlZVr55khq21NxypojH5Otv8XLt2Ld/0Td+UJ554Ivfdd18eeOCBvOpVrxq6LICFsPvBkKurq7c9yHhjYyOPPvroEKUtrN0P1z5z5kyefPLJASva21Gvf3s9XHuI42Gvh2tvbm7moYcemmsd2/b6CsP2/ehBw5iM+zeYjkU5L437lbYhz5mTnnd230MNeW3idnN5JtK8CZEAAAAAputIz0QCAAAAACESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdAmRAAAAAOgSIgEAAADQJUQCAAAAoEuIBAAAAECXEAkAAACALiESAAAAAF1CJAAAAAC6hEgAAAAAdC1UiFRVy1X1S1X1j4euBQAAAICPWqgQKcnXJXnH0EUAAAAAcLuFCZGq6oVJ/qsk3zN0LQAAAADcbmFCpCR/K8n/lOS5oQsBAAAA4HYLESJV1R9P8r7W2iOd8S5X1fWqun7z5s05VQcAAADAQoRIST43yX9dVY8neWOSz6+qB3eP1Fp7fWvtYmvt4unTp+ddIwAAAMA9ayFCpNbaN7TWXthaO5/klUl+trX2ZQOXBQAAAMDIQoRIAAAAACy2E0MXsFtr7eeT/PzAZQAAAACwg08iAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUKke1BVDd6dOnXqUPXMyoULF/at8dKlSxPXOO/tee3atUH2Y29dd+7nWZn02NlrXy/CMTjp+kxS09DrdfLkyUGOzXEN0Xam0Z08eTJJcvXq1cFr2b3dJx22s7t69erEx9gsjoskWV5ePtI67zf+2traobfx1atX970unT179rb6Jz3fHWab7tWtra0lSU6dOnWk+czyGjLuttnLpOtx4cKFI81jeXn5I9P1rveTLOPkyZM5e/bskfbRvLvd54f9xtvdFjg+DrrvHspe7WTW9zaXLl2auM7DLOcwjnovPYuamLPW2rHs7r///sbkkixMt76+fqjppm1jY2PfZS0tLQ2+nY57t76+PvV9tm3SWg7a10Meg4ddn3FqGnq9VlZWBjv2xjF0+9Dt3V25cmXsY2wWx0Vr7die/8+cOdNaO/ja1ts+Q6/Dzm4W15BJt81Oh12PjY2NI81jaWmpPfjgg8dyH0672z4/9MbbbgscH5ubmxO1x3k4c+bMYMf65ubm2HUeZTmTmNa99KLtZ+6U5HrbJ4upreHHz8WLF9v169eHLuPYuRvS3Wkfs3fDNll0szrPDLXvFnF99qtpknnOYr2GbF/jrI/2v5iWl5fz7LPPjjXupPtw3OP8OB8brbVD13+UaWdl6Ov+zuVP4zx92HmcO3cuN27cGGs5i7YPp2n7/DDOOh7X9zn3qkXcp0O3pXlcsybZpvPaHtru8Krqkdbaxb2G+TobALBQbt26NXQJsHCeeOKJoUtYCM4PAMMSIgEAC2XnM2CALffdd9/QJSwE5weAYQmRGMz6+vrQJSRJNjY29h22tKSJHNWi7Ofk4H3NbKysrAxdAsfQ5cuXhy7h2J7/z5w5k+TuOd/N4hoyxLY56jKXlpbywAMPTKma423c88N2W+D42NzcHLqEOwx5HC3i9rhbri0czfG8Q+LQFuX7pevr6/nABz4wcT2zqP/RRx/d84S4vr6eW7duTXwCH2IbP/jgg3NfZtJf1+39PNTyd4+7376e1TLnNe+Dpht3nrNar2eeeWaQIGno9Z6HlZWVtNZy5cqVoUv5iO0HLk46bKcrV67kta997UTLnMW4t27dGitIGne9do6/uro69vi7XblyZd/r0pkzZ/Lkk08m2f/a1qtt57+Htbq6mtbakQOgWV1DJtk2u7fFYbbNxsZGHn300UPPY2lpKbdu3cqrXvWq7vV+kn24srJy7IKWneeHg9ZxZ1vg+HjooYcOvO8e4pr95JNP7tlOZn1vs7m5mYceemjs8WdxD7mXadxL9xzne7N7hQdrAwAAAJDEg7UBAAAAOCIhEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQtRIhUVZ9SVT9XVe+oqker6uuGrgkAAACAjzoxdAEjzyb5i621t1XVxyZ5pKre1Fp7bOjCAAC4d1TVvsNaa3OsBAAWz0J8Eqm19t7W2ttGf38oyTuSnB22KgAA7iUHBUjjDAeAu91ChEg7VdX5JJ+e5C3DVgIAAADAtoUKkarqY5L8SJLXtNZ+d4/hl6vqelVdv3nz5vwLBAAAALhHLUyIVFUr2QqQrrXWfnSvcVprr2+tXWytXTx9+vR8CwQAAAC4hy1EiFRbXzD/3iTvaK39zaHrAQAAAOB2CxEiJfncJF+e5POr6pdH3cuHLgoAgHtH79fX/DobAPe6E0MXkCSttX+WxM9dAAAwKEERAOxvUT6JBAAAAMACEyIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4h0j2uqmbSLS8v58KFC1Od56lTp4beXLl69eoddZ09e3bPccddr0uXLh162qrq1ryo23JIhzk299pP83DU/T/OfI67S5cuHenc0pv+6tWrBy5/eXl5ZufSnd3JkyeztrY20XlhkvnvPpdN8xy+e94nT568Y92m4bD74iCTTDePY+CwplXDzvZw9uzZPcdZW1u7Y/n7jdvrlpeXD73O09gu+13jJ5nHQcfYUY7RWR9v47SPbdeuXZvqMrevuYedjvnar30fxayuE+PoXUuuXbs29WUepp2cOnVqz/cm43TjnNsmddha5lUfU9RaO5bd/fff3ziaJMeuW19fH2x7XblyZd+6zpw5c9u4k67X5ubmoafdasZ7W9RtOaSNjY1DH38799M8HHX/TzKf42pzc3Mu554rV67sufylpaXBz4vT7LbPZUdpJ715r6ys7Dl8ZWXlSMfCUffFXiaZbl776DDbado1XLlypZ05c+bAcVZXVz+y/N64vW5paWnyA2KK22X3Nf6w2/ao++Uo0x61O8iDDz4493oO6uZ9rb7X9dr3YczqOjGOca8lDz744NSWOVRbOejcNqmD3ictQn1MLsn1tk8WU1vDj5+LFy+269evD13GsXbU/yEYylDH7IkTJ3Lr1q19h++s6zDbdnv6o0y7W29ex7X9H8VRj/t5brNxa+3VNM58juuxMK/z2PLycp599tnBlj9PrbWZrVdv3kc5DmfRtidpO/M8FibdTtOubXl5+cDr4bZpbptZnKMmqeuw19mD5jHpdjnqfcZRHLT9z58/nxs3bsyxmr7jek07jmZxrzmr68RRl73TuXPn8vjjj891mbMwre3Ze590WNrycKrqkdbaxb2G+TobjGkWJ0bgeND+4aO0B3Z64oknhi4B5s5xfzvXhXuLEAnGNOvnMgCLS/uHj9Ie2Om+++4bugSYO8f97VwX7i1CJI6V9fX1wZZ9+fLlfYedOXPmSPPe3Nw80vSHMeS2HL/WMWgAABPpSURBVNLGxsahpx1iP3Gwee2T/dr/0tLddRndPpcdpZ305r2ysrLn8P36j+tu2xf7Oep2mobLly93r3urq6sf+fuo18ih9+1R67/bPfDAA0OXcBvX6vmaRfuY1XViHOOebxbtuD+Mae67g94nHZZz7wLb72FJi955sPZ0ZEYPaltaWpr6g1kX4UHQez00br+Hvo27Xns9AHKS7dKzqNtySIc5Nod6UOdR9/848znujvpw7d70+z1Ue9u8Hq69srLSVldXJzovTDL/3eeyaZ7Dd89790NTp/Ww1MPui4NMMt08joHDmlYNO9vDfg/U3flQ7d64vW5WD9Ued7uM82DXoxxjRzlGZ328jdM+tk374drb19zDTsd87de+j2JW14lx9K4l03yo9rbDtJP19fVDP9B6Fg+tnubDtT1Ue3jxYG0AAAAAejxYGwAAAIAjESIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABAlxAJAAAAgC4hEgAAAABdQiQAAAAAuoRIAAAAAHQJkQAAAADoEiIBAAAA0CVEAgAAAKBLiAQAAABA18KESFX1sqr6tap6V1V9/dD1AAAAALBDa23wLslykl9P8qIkJ5P8yyQbB01z//33t+MuybHrlpaWjlz75ubmWNtnfX198PXdr9u2ubl55HmdO3durPGWl5fblStXprJM3Z37c2lpafA6dNPpHnzwwX3PK/NoP5OeI8+dO3dgzftZXV3tzntlZeXAeUxSZ6+NTFrbvLv19fXWWmsbGxuDLL+3L/azsrIyl/pam891d5xjeXV1deLjc1q1TWO9Dnue2fbggw8eeZtuO3PmzIHTzKqt7rXsM2fO7FvnYddvnOl1unl1V65cGet4naaDrmm733MtwrV5Y2PjI/UMXcusu+MuyfW2X36z34B5dkk+O8nP7Hj9DUm+4aBpjnuINPRBPXTXC5IWOUDauQ5D16DT6fbu9gplFrnNrq2tTRQkTXIjuF94MYv1mLS2e62bNEiaV4A07+6wx/Kid9M4z0waIO21TbftFyBtd0Ns+72CpMOu3yzPZTrdYbt5Bknj/KfI9nuuRTrXDvWfOUN0x1kOCJFqa/iwquqLk7ystfZnR6+/PMlntta+dr9pLl682K5fvz6vEqeuqoYuYXAHHXu2D3AU586dy+OPP35bv0U/r+xV834mXZe9zrez2B6ttYXfzkOb5L7rbt2WO7fB3bSO0zj+z507lxs3bhxq2bst6rbdXes4dbpn5LhYXl7Os88+O5dljXvsuzYPZxGylsOqqkdaaxf3GrYoz0Ta66i+Y4tX1eWqul5V12/evDmHsgA4jp544omhS5jYcawZmC7nATjebt26NXQJMHOLEiK9O8mn7Hj9wiTv2T1Sa+31rbWLrbWLp0+fnltxABwv991339AlTOw41gxMl/MAHG/Ly8tDlwAztygh0luTvLiqPrWqTiZ5ZZKfHLgmZmhzc/PA4evr63Oq5PB66wAM54EHHrij3yK32bW1tT1r3s/q6urY466srBympEObpLZ7zaT7Yt77bgh32/Fy1PPMJOeBnjNnzhw4fIht36sJjrvLly/PbVkbGxvdcbbPSYt0rh2nbhbcfg9LmneX5OVJ/r9s/UrbN/XGP+4P1m7teD4I0K+z3f6QNL/Odnd0rfl1trup8+tsH+38OttHO7/O1j9u/TqbX2eb5rb062y6e63z62z9zq+zHR9Z9AdrH8Zxf7A2AAAAwKI5Dg/WBgAAAGCBCZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgS4gEAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAAAAA6BIiAQAAANAlRAIAAACgq1prQ9dwKFV1M8mNoeuYkucn+Z2hiwBuo13CYtI2YfFol7CYtE0O61xr7fReA45tiHQ3qarrrbWLQ9cBfJR2CYtJ24TFo13CYtI2mQVfZwMAAACgS4gEAAAAQJcQaTG8fugCgDtol7CYtE1YPNolLCZtk6nzTCQAAAAAunwSCQAAAIAuIdKAquplVfVrVfWuqvr6oeuBu11VPV5Vv1pVv1xV10f9Pr6q3lRV7xz9e2rUv6rqu0ft81eq6jN2zOfVo/HfWVWvHmp94Liqqu+rqvdV1dt39JtaW6yq+0dt/V2jaWu+awjH0z5t81ur6snRtfOXq+rlO4Z9w6id/VpV/dEd/fe8x62qT62qt4za7D+oqpPzWzs4nqrqU6rq56rqHVX1aFV93ai/6yaDECINpKqWk/zvSf5Yko0kX1pVG8NWBfeEP9Jae+mOnzv9+iQPt9ZenOTh0etkq22+eNRdTvK6ZOuCneRbknxmkj+Y5Fu2L9rA2L4/yct29ZtmW3zdaNzt6XYvC9jb92fv9vK/ja6dL22t/VSSjO5bX5nkwmia11bVcuce9ztH83pxkg8k+aqZrg3cHZ5N8hdba38gyWcl+ZpRm3LdZBBCpOH8wSTvaq3969baM0nemOQVA9cE96JXJHnD6O83JPmiHf1/oG15c5L1qnpBkj+a5E2ttfe31j6Q5E1xoYWJtNZ+Icn7d/WeSlscDfv9rbVfbFsPfvyBHfMCDrBP29zPK5K8sbX2e62130jyrmzd3+55jzv6ZMPnJ/nh0fQ72zmwj9bae1trbxv9/aEk70hyNq6bDESINJyzSX5zx+t3j/oBs9OS/JOqeqSqLo/6fVJr7b3J1kU6ySeO+u/XRrVdmI1ptcWzo7939wcO72tHX4v5vh2fXJi0bX5Ckg+21p7d1R8YU1WdT/LpSd4S100GIkQazl7fM/VTeTBbn9ta+4xsfcz3a6rq8w4Yd782qu3CfE3aFrVRmK7XJfm0JC9N8t4k/+uov7YJc1RVH5PkR5K8prX2uweNukc/bZOpESIN591JPmXH6xcmec9AtcA9obX2ntG/70vyY9n6yP1vjz7Gm9G/7xuNvl8b1XZhNqbVFt89+nt3f+AQWmu/3Vq71Vp7Lsnfzda1M5m8bf5Otr5Wc2JXf6CjqlayFSBda6396Ki36yaDECIN561JXjz6lYqT2Xow4U8OXBPctarqeVX1sdt/J/nCJG/PVrvb/nWKVyf5idHfP5nkK0a/cPFZSf7t6KPCP5PkC6vq1Ogj/V846gcczVTa4mjYh6rqs0bPYPmKHfMCJrT9JnXkT2br2plstc1XVtXvq6pPzdbDeP9F9rnHHT1r5eeSfPFo+p3tHNjH6Fr2vUne0Vr7mzsGuW4yiBP9UZiF1tqzVfW12WrMy0m+r7X26MBlwd3sk5L82OgXS08k+futtZ+uqrcm+aGq+qokTyT5ktH4P5Xk5dl6UOhTSf5MkrTW3l9VfyVbN8lJ8u2ttXEfQgokqaofTPKHkzy/qt6drV+L+Y5Mry1eydavTK0m+X9GHdCxT9v8w1X10mx9veXxJF+dJK21R6vqh5I8lq1fj/qa1tqt0Xz2u8f9y0neWFV/NckvZeuNMXCwz03y5Ul+tap+edTvG+O6yUBq6z8FAAAAAGB/vs4GAAAAQJcQCQAAAIAuIRIAAAAAXUIkAAAAALqESAAAAAB0CZEAgD1VVRuje3w07reOXp/YMf3PHzDd39ox3veP+v1mVd1xb7Jj3r35v7eqfrqqPnOPebykqt5QVU9W1TOjf/+vqnrJhNvkH1XV35lkmlmpqvOj9f7KHf2+f3ufTDCfrxzN5/x0KzxwmT9bVd82+vszq+pWVX3MrnH+QlX9yl7HBAAwjBP9UQCAe9Rn73r9Y0n+ZZJv3dHv9zrz+JUkX71H//fuev1UkjNJ/kiSh3cN+7IkH0rysZ35n0/yPyf5p1X1Ga21x5Kkqi4l+Ykk70zyjUl+YzTua5K8rape0Vp7qLMeqarPS/IFST6tN+6A/kqSvz3hNP93tvb17n0yE1VVST49yXaQeDHJr7XW/t2uUf+PJH85yauT/L151AYAHEyIBADsqbX25p2vq+r3kvzO7v4dHxpz/A8k+VdJvjw7QqSq+kNJXpTkB7IVJhw0/zdX1ZuzFRL9uSR/vqo+IckbsxV+fX5r7T+Mxv2FqvqhJD+b5I1V9ZLW2r/p1Pg/JvlHrbUnx1ifQbTWfv0Q09xMcnMG5ezn05KsJ3lk9Prijr931vV0Vf1Akr8UIRIALAQfDwYAFsUPJPlTVbW2o99XJPl/kzw+zgxaa49nKxD5j0e9/myST0jydTsCpO1x/0O2Po30CaPx9lVVZ5L8sSR/f1f/7a+CfU5V/VBVfaiqfruqvmE0/GVV9UtV9e+r6q1Vdf8e8/5vqurNVfVUVX2wqv5hVd23a5y1qnptVf2bqvp3VfWTSV64x7zu+DpbVT2vqr6jqn69qn6vqn6rqn6kqj5p1zqc3zHNSlX91ap6fPT1v8dHr1d2jLP9dbqvrqpvH32d8IOjr/zdUdsO9yd5344w7v4kb9tn3Dcm2aiqzzlgfgDAnAiRAICZqqoTe3S1x6g/kqSSfNFout+X5EuyFS6Nu6yPS/LxST446rWZ5Ldaa2/da/zW2r9I8ttJPr8z6y9Ispzkn+0z/A1JfjXJn0zy40n+WlV9Z5LvSvKdSf50kucl+fGqOrmj3j+XrfV+LMn/3969hdxRXQEc/y9Bq+ClXpNaxAoa+hBFkFitUvtQ1FqttGrRqKgRb1BaKkprrTRWoRFv6UMfBHtTk7QavBUJjSgxXqN5iIqShFQr1lDbRDFaLEZdfdh7zGR6zjfnqzn9Qvn/4DBn9syevec8HRZrr306ZWneTMqSvPbyvdsoga5bgG8Da+gEtAapYz0MfA/4LXAy8F3gLWDPCbr+DvgR5bc/mZIJ9MPa3nUVJWg3B/g+ZWncgs48vtrUrqIEhvZrnR8K3FLPl3WevQrYBJzY966SJGn8XM4mSZLG6Rhg84D2M4DF7YbM/GdE3EvJPloInAp8BrgHuHzYAK1i2wcCN1OCPXfXtgPoz2L6S71vIkcB6+vSr0HuzMzr6nyWUYJJlwMzMvPV2r4DpTbT0ZQg0a6UANNvMnNO631WAGuBC4H5tfj3bODqzJxXb1ta+1/aM+9z6ninZuaDrfbFQ+4nImYCZwHXZubc1ngfAddFxLzMfKHV5bXMnN3qvy9wY0Tsn5nra/NKSh2kZux7gD9QamBdC3ylXtuqLlJmfhwRL1B+f0mSNMXMRJIkSeP0PDBrwKdbPLtxB/C1iJhOCSY9kJmbJnh+E6TaDKwDvgxcmpn31+uDMp66RrlnfyauG7Sk+ZKZH9a5rG0CSNXqemwCVkcDuwML2llawF/rvU1g5UuU/2x3s7XfjzDv4ymZWA/23rlFM+5dnfbm/LhO+0Od8xfr8ZMleZn5XmauAv5MKWp+bz3fB3g2M1fVz7oB8/kH5feXJElTzEwkSZI0Tu9l5spJ3P8oZZewHwAnAN/suf95yjKvpCxLeyMzs3X9dcrysIkcWJ8zkZ2ZeCe6tzvnHwxpa54FsF89DtsZrun/uXp8s3O9ez7I3sBkC4HvVY/d3dr+1rneeKtz3vxOzXu2s8WOrddfrG3HAE/W75mZHw2Yz/vALqNPX5IkjYtBJEmStN2oy5cWUHZC+zuwtKdLX5DqEUpm06xBdZEi4khgGiV4NZGNwEE990xWsxvc+cBLA66/W49NMGca8Err+rQRxthAfxCtqwkKTadkDtE6hy3zHkkt2P1qp/n91vfjgB8Dr1GylLr2oryHJEmaYi5nkyRJ25tfA38Erh+SmTIZt1Myen4RETu3L9Tz+ZSgye09z1kNHNDKqNkWnqIEig7OzJUDPmvqfSuAj4HvdPqfOcIYS4HpEXHKJOb12JDnn12PyyfxLID1bFnGuIbym8+iFPj+F6Xe0Sxg2BwPqv0kSdIUMxNJkiSN024RMago8tutIMlWMnMtdYe2TyszN0TEWcB9wNMRcSslK+YLlCVzXwS+lZl92TXLKQWgD2P4dvSTndumiLgS+GUtRr0EeAf4PCU7Z1lmLszMNRGxEPhZLc79HGW3uJNGGOYu4CJgUUT8nBKQ2o2yVHB+Zq7udsjMlyJiETC3Bs2eotRvugZY1CmqPcp7fgCsjIi9gUOAMzNzVUScVt9xxbC+EfFZYAZw02TGlCRJ42EQSZIkjdNhwNMD2h+ibB0/dpn5p4g4grJkah6lmPNGyhK2czLz5REe8zglo+YUtlEQqc7ttoh4nbJ8bzawI6WG0XLK9vaNSyg7l10B7FTnPht4ouf5myPieOCnwMX1uBF4kv+sZdR2HmXp3BzgJ5R3v4ESSPtvnUhZotjUn/o68KuePt+g1JK671OMK0mStpHYuvakJEmSBomIuZQlXTPSP1D/ExGxBNiQmedO9VwkSZJBJEmSpJFExB7AOuCyzFw81fP5fxcRhwPPADMzc91Uz0eSJFlYW5IkaSSZ+Q5wLmU5mcZvOnCBASRJkrYfZiJJkiRJkiSpl5lIkiRJkiRJ6mUQSZIkSZIkSb0MIkmSJEmSJKmXQSRJkiRJkiT1MogkSZIkSZKkXgaRJEmSJEmS1Ovfrq9vd/JkUYkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1440x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "pdseries = gdf.label.copy()\n",
    "f=plt.figure(figsize=(20,10))\n",
    "plt.xlabel(f'TIEMPO (medicion#)', fontsize=16)\n",
    "plt.ylabel('LABEL', fontsize=16)\n",
    "#rangotiks=range( #rangolabels = list(rangotiks) ##plt.xticks(ticks=rangotiks, labels=rangolabels,  rotation=90)\n",
    "pdseries=pdseries.apply(int)\n",
    "plt.plot(pdseries, 'ko')    \n",
    "del(pdseries)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 167,
   "metadata": {
    "collapsed": true,
    "jupyter": {
     "outputs_hidden": true
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>date</th>\n",
       "      <th>marker-size</th>\n",
       "      <th>marker-color</th>\n",
       "      <th>marker-symbol</th>\n",
       "      <th>stroke</th>\n",
       "      <th>stroke-width</th>\n",
       "      <th>stroke-opacity</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2012-03-29</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.32826 39.97693)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2012-03-29</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.32826 39.97693)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2012-03-29</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.32826 39.97693, 116.32768 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2012-03-30</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.31087 39.97587)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2012-03-30</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.31087 39.97587, 116.31299 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>2012-03-31</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.32919 39.97760)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>2012-03-31</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.32919 39.97760, 116.32914 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>2012-04-01</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (114.34205 30.53053)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>2012-04-01</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (114.34205 30.53053, 114.34243 30.5...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>2012-04-05</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.30267 39.97466)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>2012-04-05</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.30267 39.97466, 116.30218 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>2012-04-06</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.33111 39.97468)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>2012-04-06</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.33111 39.97468, 116.33178 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>2012-04-07</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.33201 39.97705)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>2012-04-07</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.33201 39.97705, 116.33196 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>2012-04-08</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.30301 39.97766)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>2012-04-08</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.30301 39.97766, 116.30316 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>2012-04-09</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.33199 39.97625)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>2012-04-09</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.33199 39.97625, 116.33195 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>2012-04-10</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.30401 39.97790)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>2012-04-10</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.30401 39.97790, 116.30325 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>2012-04-11</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.31097 39.97473)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>2012-04-11</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.31097 39.97473, 116.30519 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>2012-04-12</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.31049 39.97949)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>2012-04-12</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.31049 39.97949, 116.31054 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>2012-04-13</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.32744 39.97975)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>2012-04-13</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.32744 39.97975, 116.32741 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>2012-04-23</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.30148 39.97494)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>2012-04-23</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.30148 39.97494, 116.30160 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>2012-04-24</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.33348 39.97481)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>2012-04-24</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.33348 39.97481, 116.33338 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>2012-05-07</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.30474 39.97890)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>2012-05-07</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.30474 39.97890, 116.30481 39.9...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>2012-05-08</td>\n",
       "      <td>large</td>\n",
       "      <td>#d77070</td>\n",
       "      <td>entrance</td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>POINT (116.33231 39.97667)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>2012-05-08</td>\n",
       "      <td>large</td>\n",
       "      <td></td>\n",
       "      <td></td>\n",
       "      <td>#f71a3a</td>\n",
       "      <td>2</td>\n",
       "      <td>0.7</td>\n",
       "      <td>LINESTRING (116.33231 39.97667, 116.33184 39.9...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          date marker-size marker-color marker-symbol   stroke  stroke-width  \\\n",
       "0   2012-03-29       large      #d77070      entrance  #f71a3a             2   \n",
       "1   2012-03-29       large      #d77070      entrance  #f71a3a             2   \n",
       "2   2012-03-29       large                             #f71a3a             2   \n",
       "3   2012-03-30       large      #d77070      entrance  #f71a3a             2   \n",
       "4   2012-03-30       large                             #f71a3a             2   \n",
       "5   2012-03-31       large      #d77070      entrance  #f71a3a             2   \n",
       "6   2012-03-31       large                             #f71a3a             2   \n",
       "7   2012-04-01       large      #d77070      entrance  #f71a3a             2   \n",
       "8   2012-04-01       large                             #f71a3a             2   \n",
       "9   2012-04-05       large      #d77070      entrance  #f71a3a             2   \n",
       "10  2012-04-05       large                             #f71a3a             2   \n",
       "11  2012-04-06       large      #d77070      entrance  #f71a3a             2   \n",
       "12  2012-04-06       large                             #f71a3a             2   \n",
       "13  2012-04-07       large      #d77070      entrance  #f71a3a             2   \n",
       "14  2012-04-07       large                             #f71a3a             2   \n",
       "15  2012-04-08       large      #d77070      entrance  #f71a3a             2   \n",
       "16  2012-04-08       large                             #f71a3a             2   \n",
       "17  2012-04-09       large      #d77070      entrance  #f71a3a             2   \n",
       "18  2012-04-09       large                             #f71a3a             2   \n",
       "19  2012-04-10       large      #d77070      entrance  #f71a3a             2   \n",
       "20  2012-04-10       large                             #f71a3a             2   \n",
       "21  2012-04-11       large      #d77070      entrance  #f71a3a             2   \n",
       "22  2012-04-11       large                             #f71a3a             2   \n",
       "23  2012-04-12       large      #d77070      entrance  #f71a3a             2   \n",
       "24  2012-04-12       large                             #f71a3a             2   \n",
       "25  2012-04-13       large      #d77070      entrance  #f71a3a             2   \n",
       "26  2012-04-13       large                             #f71a3a             2   \n",
       "27  2012-04-23       large      #d77070      entrance  #f71a3a             2   \n",
       "28  2012-04-23       large                             #f71a3a             2   \n",
       "29  2012-04-24       large      #d77070      entrance  #f71a3a             2   \n",
       "30  2012-04-24       large                             #f71a3a             2   \n",
       "31  2012-05-07       large      #d77070      entrance  #f71a3a             2   \n",
       "32  2012-05-07       large                             #f71a3a             2   \n",
       "33  2012-05-08       large      #d77070      entrance  #f71a3a             2   \n",
       "34  2012-05-08       large                             #f71a3a             2   \n",
       "\n",
       "    stroke-opacity                                           geometry  \n",
       "0              0.7                         POINT (116.32826 39.97693)  \n",
       "1              0.7                         POINT (116.32826 39.97693)  \n",
       "2              0.7  LINESTRING (116.32826 39.97693, 116.32768 39.9...  \n",
       "3              0.7                         POINT (116.31087 39.97587)  \n",
       "4              0.7  LINESTRING (116.31087 39.97587, 116.31299 39.9...  \n",
       "5              0.7                         POINT (116.32919 39.97760)  \n",
       "6              0.7  LINESTRING (116.32919 39.97760, 116.32914 39.9...  \n",
       "7              0.7                         POINT (114.34205 30.53053)  \n",
       "8              0.7  LINESTRING (114.34205 30.53053, 114.34243 30.5...  \n",
       "9              0.7                         POINT (116.30267 39.97466)  \n",
       "10             0.7  LINESTRING (116.30267 39.97466, 116.30218 39.9...  \n",
       "11             0.7                         POINT (116.33111 39.97468)  \n",
       "12             0.7  LINESTRING (116.33111 39.97468, 116.33178 39.9...  \n",
       "13             0.7                         POINT (116.33201 39.97705)  \n",
       "14             0.7  LINESTRING (116.33201 39.97705, 116.33196 39.9...  \n",
       "15             0.7                         POINT (116.30301 39.97766)  \n",
       "16             0.7  LINESTRING (116.30301 39.97766, 116.30316 39.9...  \n",
       "17             0.7                         POINT (116.33199 39.97625)  \n",
       "18             0.7  LINESTRING (116.33199 39.97625, 116.33195 39.9...  \n",
       "19             0.7                         POINT (116.30401 39.97790)  \n",
       "20             0.7  LINESTRING (116.30401 39.97790, 116.30325 39.9...  \n",
       "21             0.7                         POINT (116.31097 39.97473)  \n",
       "22             0.7  LINESTRING (116.31097 39.97473, 116.30519 39.9...  \n",
       "23             0.7                         POINT (116.31049 39.97949)  \n",
       "24             0.7  LINESTRING (116.31049 39.97949, 116.31054 39.9...  \n",
       "25             0.7                         POINT (116.32744 39.97975)  \n",
       "26             0.7  LINESTRING (116.32744 39.97975, 116.32741 39.9...  \n",
       "27             0.7                         POINT (116.30148 39.97494)  \n",
       "28             0.7  LINESTRING (116.30148 39.97494, 116.30160 39.9...  \n",
       "29             0.7                         POINT (116.33348 39.97481)  \n",
       "30             0.7  LINESTRING (116.33348 39.97481, 116.33338 39.9...  \n",
       "31             0.7                         POINT (116.30474 39.97890)  \n",
       "32             0.7  LINESTRING (116.30474 39.97890, 116.30481 39.9...  \n",
       "33             0.7                         POINT (116.33231 39.97667)  \n",
       "34             0.7  LINESTRING (116.33231 39.97667, 116.33184 39.9...  "
      ]
     },
     "execution_count": 167,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "dfxb = geopandas.read_file(\"geolife_159_linesimpl2.geojson\")\n",
    "dfxb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 188,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>time</th>\n",
       "      <th>lat</th>\n",
       "      <th>lon</th>\n",
       "      <th>alt</th>\n",
       "      <th>label</th>\n",
       "      <th>user</th>\n",
       "      <th>idx_traj</th>\n",
       "      <th>vel</th>\n",
       "      <th>velMA</th>\n",
       "      <th>velprep</th>\n",
       "      <th>velMA2</th>\n",
       "      <th>velMA3</th>\n",
       "      <th>velprep2</th>\n",
       "      <th>pruebas</th>\n",
       "      <th>vel_label</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2007-08-04 03:30:32</td>\n",
       "      <td>39.921712</td>\n",
       "      <td>116.472343</td>\n",
       "      <td>13</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>0</td>\n",
       "      <td>2.77</td>\n",
       "      <td>5.082500</td>\n",
       "      <td>2.77</td>\n",
       "      <td>5.442391</td>\n",
       "      <td>11.068403</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.867746</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>2007-08-04 03:30:45</td>\n",
       "      <td>39.921562</td>\n",
       "      <td>116.472303</td>\n",
       "      <td>13</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>12</td>\n",
       "      <td>7.30</td>\n",
       "      <td>5.142813</td>\n",
       "      <td>7.30</td>\n",
       "      <td>5.129901</td>\n",
       "      <td>12.325189</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.750150</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>2007-08-04 03:31:04</td>\n",
       "      <td>39.921423</td>\n",
       "      <td>116.472320</td>\n",
       "      <td>20</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>24</td>\n",
       "      <td>2.31</td>\n",
       "      <td>5.002031</td>\n",
       "      <td>2.31</td>\n",
       "      <td>5.187508</td>\n",
       "      <td>13.462401</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.619431</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>2007-08-04 03:31:16</td>\n",
       "      <td>39.921278</td>\n",
       "      <td>116.472347</td>\n",
       "      <td>15</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>36</td>\n",
       "      <td>5.77</td>\n",
       "      <td>6.190469</td>\n",
       "      <td>5.77</td>\n",
       "      <td>5.655570</td>\n",
       "      <td>14.390125</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.464143</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48</th>\n",
       "      <td>2007-08-04 03:31:28</td>\n",
       "      <td>39.921177</td>\n",
       "      <td>116.472248</td>\n",
       "      <td>25</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>48</td>\n",
       "      <td>6.73</td>\n",
       "      <td>6.255000</td>\n",
       "      <td>6.73</td>\n",
       "      <td>5.595190</td>\n",
       "      <td>15.114886</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.276402</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>935520</th>\n",
       "      <td>2009-03-21 05:33:52</td>\n",
       "      <td>39.136148</td>\n",
       "      <td>117.217455</td>\n",
       "      <td>-43</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>2687</td>\n",
       "      <td>5.15</td>\n",
       "      <td>6.206094</td>\n",
       "      <td>5.15</td>\n",
       "      <td>5.848559</td>\n",
       "      <td>7.167104</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.934888</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>935532</th>\n",
       "      <td>2009-03-21 05:34:05</td>\n",
       "      <td>39.136200</td>\n",
       "      <td>117.217614</td>\n",
       "      <td>-49</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>2699</td>\n",
       "      <td>4.52</td>\n",
       "      <td>5.175000</td>\n",
       "      <td>4.52</td>\n",
       "      <td>5.683420</td>\n",
       "      <td>6.858121</td>\n",
       "      <td>20.0</td>\n",
       "      <td>19.277786</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>935544</th>\n",
       "      <td>2009-03-21 05:34:19</td>\n",
       "      <td>39.136226</td>\n",
       "      <td>117.217790</td>\n",
       "      <td>-59</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>2711</td>\n",
       "      <td>4.42</td>\n",
       "      <td>6.390000</td>\n",
       "      <td>4.42</td>\n",
       "      <td>6.195566</td>\n",
       "      <td>6.521719</td>\n",
       "      <td>20.0</td>\n",
       "      <td>18.649420</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>935556</th>\n",
       "      <td>2009-03-21 05:34:32</td>\n",
       "      <td>39.136203</td>\n",
       "      <td>117.217996</td>\n",
       "      <td>-75</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>2723</td>\n",
       "      <td>4.75</td>\n",
       "      <td>6.632187</td>\n",
       "      <td>4.75</td>\n",
       "      <td>6.646634</td>\n",
       "      <td>6.322574</td>\n",
       "      <td>5.0</td>\n",
       "      <td>18.049509</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>935568</th>\n",
       "      <td>2009-03-21 05:34:45</td>\n",
       "      <td>39.136256</td>\n",
       "      <td>117.218209</td>\n",
       "      <td>-69</td>\n",
       "      <td>0</td>\n",
       "      <td>10</td>\n",
       "      <td>2735</td>\n",
       "      <td>4.52</td>\n",
       "      <td>6.025312</td>\n",
       "      <td>4.52</td>\n",
       "      <td>6.306506</td>\n",
       "      <td>6.108999</td>\n",
       "      <td>5.0</td>\n",
       "      <td>17.476842</td>\n",
       "      <td>32.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>77965 rows × 15 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                      time        lat         lon  alt  label  user  idx_traj  \\\n",
       "0      2007-08-04 03:30:32  39.921712  116.472343   13      0    10         0   \n",
       "12     2007-08-04 03:30:45  39.921562  116.472303   13      0    10        12   \n",
       "24     2007-08-04 03:31:04  39.921423  116.472320   20      0    10        24   \n",
       "36     2007-08-04 03:31:16  39.921278  116.472347   15      0    10        36   \n",
       "48     2007-08-04 03:31:28  39.921177  116.472248   25      0    10        48   \n",
       "...                    ...        ...         ...  ...    ...   ...       ...   \n",
       "935520 2009-03-21 05:33:52  39.136148  117.217455  -43      0    10      2687   \n",
       "935532 2009-03-21 05:34:05  39.136200  117.217614  -49      0    10      2699   \n",
       "935544 2009-03-21 05:34:19  39.136226  117.217790  -59      0    10      2711   \n",
       "935556 2009-03-21 05:34:32  39.136203  117.217996  -75      0    10      2723   \n",
       "935568 2009-03-21 05:34:45  39.136256  117.218209  -69      0    10      2735   \n",
       "\n",
       "         vel     velMA  velprep    velMA2     velMA3  velprep2    pruebas  \\\n",
       "0       2.77  5.082500     2.77  5.442391  11.068403      20.0  19.867746   \n",
       "12      7.30  5.142813     7.30  5.129901  12.325189      20.0  19.750150   \n",
       "24      2.31  5.002031     2.31  5.187508  13.462401      20.0  19.619431   \n",
       "36      5.77  6.190469     5.77  5.655570  14.390125      20.0  19.464143   \n",
       "48      6.73  6.255000     6.73  5.595190  15.114886      20.0  19.276402   \n",
       "...      ...       ...      ...       ...        ...       ...        ...   \n",
       "935520  5.15  6.206094     5.15  5.848559   7.167104      20.0  19.934888   \n",
       "935532  4.52  5.175000     4.52  5.683420   6.858121      20.0  19.277786   \n",
       "935544  4.42  6.390000     4.42  6.195566   6.521719      20.0  18.649420   \n",
       "935556  4.75  6.632187     4.75  6.646634   6.322574       5.0  18.049509   \n",
       "935568  4.52  6.025312     4.52  6.306506   6.108999       5.0  17.476842   \n",
       "\n",
       "        vel_label  \n",
       "0            32.0  \n",
       "12           32.0  \n",
       "24           32.0  \n",
       "36           32.0  \n",
       "48           32.0  \n",
       "...           ...  \n",
       "935520       32.0  \n",
       "935532       32.0  \n",
       "935544       32.0  \n",
       "935556       32.0  \n",
       "935568       32.0  \n",
       "\n",
       "[77965 rows x 15 columns]"
      ]
     },
     "execution_count": 188,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df.iloc[::12,:]"
   ]
  },
  {
   "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.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}