04-lisbon-cicleways.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Import packages"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import geojson\n",
    "import json\n",
    "import folium\n",
    "import geopandas as gpd\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Open Lisbon's municipality data on cycleways"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "with open('Data/Ciclovias.geojson') as f:\n",
    "    ciclovias = geojson.load(f)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Correct some mistakes in the data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(len(ciclovias['features'])):\n",
    "    ciclovia = ciclovias['features'][i]\n",
    "    ciclovia['IDunico'] = ciclovia['properties']['COD_SIG']\n",
    "      \n",
    "    if ciclovia['IDunico'] == '0000455':\n",
    "        ciclovia['properties']['ANO'] = 2000\n",
    "\n",
    "    if ciclovia['IDunico'] == '0000371':\n",
    "        ciclovia['properties']['ANO'] = 2001\n",
    "        \n",
    "    if ciclovia['IDunico'] == '6708':\n",
    "        ciclovia['properties']['ANO'] = 2003\n",
    "        \n",
    "    if ciclovia['IDunico'] == '0000204' or \\\n",
    "        ciclovia['IDunico'] == '0000080' or \\\n",
    "        ciclovia['IDunico'] == '0000205':\n",
    "        ciclovia['properties']['ANO'] = 2005\n",
    "        \n",
    "    if ciclovia['IDunico'] == '0000025' or \\\n",
    "        ciclovia['IDunico'] == '0000173' or \\\n",
    "        ciclovia['IDunico'] == '0000085':\n",
    "        ciclovia['properties']['ANO'] = 2008\n",
    "    if ciclovia['IDunico'] == '0000224' or \\\n",
    "        ciclovia['IDunico'] == '000002' or \\\n",
    "        ciclovia['IDunico'] == '0000015' or \\\n",
    "        ciclovia['IDunico'] == '0000012' or \\\n",
    "        ciclovia['IDunico'] == '0000138' or \\\n",
    "        ciclovia['IDunico'] == '0000144' or \\\n",
    "        ciclovia['IDunico'] == '0000237' or \\\n",
    "        ciclovia['IDunico'] == '0000145' or \\\n",
    "        ciclovia['IDunico'] == '0000225' or \\\n",
    "        ciclovia['IDunico'] == '0000132' or \\\n",
    "        ciclovia['IDunico'] == '0000205' or \\\n",
    "        ciclovia['IDunico'] == '0000174' or \\\n",
    "        ciclovia['IDunico'] == '0000024' or \\\n",
    "        ciclovia['IDunico'] == '0000210' or \\\n",
    "        ciclovia['IDunico'] == '0000081' or \\\n",
    "        ciclovia['IDunico'] == '0000028' or \\\n",
    "        ciclovia['IDunico'] == '0000010' or \\\n",
    "        ciclovia['IDunico'] == '0000011' or \\\n",
    "        ciclovia['IDunico'] == '0000214' or \\\n",
    "        ciclovia['IDunico'] == '0000013' or \\\n",
    "        ciclovia['IDunico'] == '0000212' or \\\n",
    "        ciclovia['IDunico'] == '000008' or \\\n",
    "        ciclovia['IDunico'] == '0000216' or \\\n",
    "        ciclovia['IDunico'] == '0000009' or \\\n",
    "        ciclovia['IDunico'] == '0000202' or \\\n",
    "        ciclovia['IDunico'] == '0000184' or \\\n",
    "        ciclovia['IDunico'] == '000007':\n",
    "        ciclovia['properties']['ANO'] = 2009\n",
    "    if ciclovia['IDunico'] == '0000036' or \\\n",
    "        ciclovia['IDunico'] == '0000026' or \\\n",
    "        ciclovia['IDunico'] == '14296' or \\\n",
    "        ciclovia['IDunico'] == '000019' or \\\n",
    "        ciclovia['IDunico'] == '152515' or \\\n",
    "        ciclovia['IDunico'] == '0000006' or \\\n",
    "        ciclovia['IDunico'] == '12211' or \\\n",
    "        ciclovia['IDunico'] == '3185' or \\\n",
    "        ciclovia['IDunico'] == '152515' or \\\n",
    "        ciclovia['IDunico'] == '23242' or \\\n",
    "        ciclovia['IDunico'] == '0000285' or \\\n",
    "        ciclovia['IDunico'] == '159396' or \\\n",
    "        ciclovia['IDunico'] == '0000284' or \\\n",
    "        ciclovia['IDunico'] == '159396' or \\\n",
    "        ciclovia['IDunico'] == '159399' or \\\n",
    "        ciclovia['IDunico'] == '0000041' or \\\n",
    "        ciclovia['IDunico'] == '19676' or \\\n",
    "        ciclovia['IDunico'] == '12210' or \\\n",
    "        ciclovia['IDunico'] == '14286' or \\\n",
    "        ciclovia['IDunico'] == '000016' or \\\n",
    "        ciclovia['IDunico'] == '0000061' or \\\n",
    "        ciclovia['IDunico'] == '0000064' or \\\n",
    "        ciclovia['IDunico'] == '0000065' or \\\n",
    "        ciclovia['IDunico'] == '0000066' or \\\n",
    "        ciclovia['IDunico'] == '0000063' or \\\n",
    "        ciclovia['IDunico'] == '0000062' or \\\n",
    "        ciclovia['IDunico'] == '6244' or \\\n",
    "        ciclovia['IDunico'] == '6244' or \\\n",
    "        ciclovia['IDunico'] == '000017':\n",
    "        ciclovia['properties']['ANO'] = 2010\n",
    "    if ciclovia['IDunico'] == '0000017' or \\\n",
    "        ciclovia['IDunico'] == '0000018' or \\\n",
    "        ciclovia['IDunico'] == '0000146' or \\\n",
    "        ciclovia['IDunico'] == '0000181' or \\\n",
    "        ciclovia['IDunico'] == '0000147':\n",
    "        ciclovia['properties']['ANO'] = 2011\n",
    "    if ciclovia['IDunico'] == '0000266' or \\\n",
    "        ciclovia['IDunico'] == '0000265' or \\\n",
    "        ciclovia['IDunico'] == '0000028' or \\\n",
    "        ciclovia['IDunico'] == '0000188' or \\\n",
    "        ciclovia['IDunico'] == '0000078' or \\\n",
    "        ciclovia['IDunico'] == '0000037':\n",
    "        ciclovia['properties']['ANO'] = 2012\n",
    "        \n",
    "    if ciclovia['IDunico'] == '18320' or \\\n",
    "        ciclovia['IDunico'] == '159178' or \\\n",
    "        ciclovia['IDunico'] == '159044' or \\\n",
    "        ciclovia['IDunico'] == '158526' or \\\n",
    "        ciclovia['IDunico'] == '159033' or \\\n",
    "        ciclovia['IDunico'] == '159043' or \\\n",
    "        ciclovia['IDunico'] == '159023' or \\\n",
    "        ciclovia['IDunico'] == '159031' or \\\n",
    "        ciclovia['IDunico'] == '159021' or \\\n",
    "        ciclovia['IDunico'] == '159179' or \\\n",
    "        ciclovia['IDunico'] == '159022':\n",
    "        ciclovia['properties']['ANO'] = 2013\n",
    "\n",
    "    if ciclovia['IDunico'] == '0000045' or \\\n",
    "        ciclovia['IDunico'] == '0000140' or \\\n",
    "        ciclovia['IDunico'] == '15102' or \\\n",
    "        ciclovia['IDunico'] == '2607' or \\\n",
    "        ciclovia['IDunico'] == '14168' or \\\n",
    "        ciclovia['IDunico'] == '0000067' or \\\n",
    "        ciclovia['IDunico'] == '0000045' or \\\n",
    "        ciclovia['IDunico'] == '0000367' or \\\n",
    "        ciclovia['IDunico'] == '0000045' or \\\n",
    "        ciclovia['IDunico'] == 'B0000367':\n",
    "        ciclovia['properties']['ANO'] = 2014\n",
    "\n",
    "    if ciclovia['IDunico'] == '5777' or \\\n",
    "        ciclovia['IDunico'] == '5684' or \\\n",
    "        ciclovia['IDunico'] == '0000215' or \\\n",
    "        ciclovia['IDunico'] == '5593' or \\\n",
    "        ciclovia['IDunico'] == '0000243' or \\\n",
    "        ciclovia['IDunico'] == '0000057' or \\\n",
    "        ciclovia['IDunico'] == 'B0000038' or \\\n",
    "        ciclovia['IDunico'] == 'B0000183' or \\\n",
    "        ciclovia['IDunico'] == '0000240':\n",
    "        ciclovia['properties']['ANO'] = 2015\n",
    "\n",
    "    if ciclovia['IDunico'] == '18155' or \\\n",
    "        ciclovia['IDunico'] == '17340' or \\\n",
    "        ciclovia['IDunico'] == '158357' or \\\n",
    "        ciclovia['IDunico'] == '18153' or \\\n",
    "        ciclovia['IDunico'] == '17343' or \\\n",
    "        ciclovia['IDunico'] == '18029' or \\\n",
    "        ciclovia['IDunico'] == '0000077' or \\\n",
    "        ciclovia['IDunico'] == '18031' :\n",
    "        ciclovia['properties']['ANO'] = 2016\n",
    "\n",
    "    if ciclovia['IDunico'] == '0000373':\n",
    "        ciclovia['properties']['ANO'] = 2018\n",
    "        \n",
    "    if ciclovia['IDunico'] == '0000372' or \\\n",
    "        ciclovia['IDunico'] == '0000392' or \\\n",
    "        ciclovia['IDunico'] == '0000391' or \\\n",
    "        ciclovia['IDunico'] == '4198' or \\\n",
    "        ciclovia['IDunico'] == '15403' or \\\n",
    "        ciclovia['IDunico'] == '153668' or \\\n",
    "        ciclovia['IDunico'] == '4128' or \\\n",
    "        ciclovia['IDunico'] == '4671' or \\\n",
    "        ciclovia['IDunico'] == '156600' or \\\n",
    "        ciclovia['IDunico'] == '4614' or \\\n",
    "        ciclovia['IDunico'] == '4198' or \\\n",
    "        ciclovia['IDunico'] == '0000325' or \\\n",
    "        ciclovia['IDunico'] == '0000325' or \\\n",
    "        ciclovia['IDunico'] == '15403' or \\\n",
    "        ciclovia['IDunico'] == '15620' or \\\n",
    "        ciclovia['IDunico'] == '0000393' or \\\n",
    "        ciclovia['IDunico'] == '0000390':\n",
    "        ciclovia['properties']['ANO'] = 2019   "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Plot all cycleways as of 2020"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_freguesias = gpd.read_file(\"Data/Lisboa_Freguesias/Lisboa_Freguesias_CAOP2015_TM06.shp\")\n",
    "gdf_freguesias = gdf_freguesias.to_crs(epsg=4326)\n",
    "gdf_freguesias.geometry.index = gdf_freguesias['DICOFRE']\n",
    "\n",
    "style_1 = {'fillColor': '#1ED674', 'color': '#1ED674'}\n",
    "style_2 = {'fillColor': '#D6A86E', 'color': '#D6A86E'}\n",
    "style_3 = {'fillColor': '#233BF7', 'color': '#233BF7'}\n",
    "\n",
    "unique_zones = gdf_freguesias['DICOFRE'].unique()\n",
    "m = folium.Map([38.748662, -9.145801],\n",
    "                 zoom_start=12,\n",
    "                 tiles='cartodbpositron')\n",
    "\n",
    "for i in range(len(ciclovias['features'])):\n",
    "    ciclovia = ciclovias['features'][i]\n",
    "    \n",
    "    try:\n",
    "        ano = int(ciclovia['properties']['ANO'])\n",
    "        if ano <= 2020:\n",
    "\n",
    "            folium.GeoJson(\n",
    "                ciclovia,\n",
    "                name='geojson',\n",
    "                style_function=lambda x:style_1\n",
    "            ).add_to(m)\n",
    "    except ValueError:\n",
    "        folium.GeoJson(\n",
    "                ciclovia,\n",
    "                name='geojson',\n",
    "                style_function=lambda x:style_3\n",
    "            ).add_to(m)\n",
    "    \n",
    "for i, zone in enumerate(unique_zones):\n",
    "    mask = gdf_freguesias['DICOFRE'] == zone\n",
    "\n",
    "    example_or = gdf_freguesias.loc[mask]\n",
    "    folium.GeoJson(example_or,name='polygon_or',style_function=lambda x:style_2).add_to(m)\n",
    "\n",
    "m"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Create maps of Lisbon's cycleways over time"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_freguesias = gpd.read_file(\"Data/Lisboa_Freguesias/Lisboa_Freguesias_CAOP2015_TM06.shp\")\n",
    "gdf_freguesias = gdf_freguesias.to_crs(epsg=4326)\n",
    "gdf_freguesias.geometry.index = gdf_freguesias['DICOFRE']\n",
    "style_0 = {'fillColor': '#D6A86E', 'color': '#D6A86E'}\n",
    "style_1 = {'fillColor': '#1ED674', 'color': '#1ED674'}\n",
    "style_2 = {'fillColor': '#5CC1F7', 'color': '#5CC1F7'}\n",
    "style_3 = {'fillColor': '#472EAB', 'color': '#472EAB'}\n",
    "style_4 = {'fillColor': '#AB2E6E', 'color': '#AB2E6E'}\n",
    "\n",
    "unique_zones = gdf_freguesias['DICOFRE'].unique()\n",
    "\n",
    "\n",
    "for current_year in range(2013, 2021):\n",
    "\n",
    "    m = folium.Map([38.747156, -9.156788],\n",
    "                 zoom_start=13,\n",
    "                 tiles='cartodbpositron')\n",
    "    \n",
    "    for i in range(len(ciclovias['features'])):\n",
    "        ciclovia = ciclovias['features'][i]\n",
    "\n",
    "        try:\n",
    "            ano = int(ciclovia['properties']['ANO'])\n",
    "            if ano <= current_year:\n",
    "                folium.GeoJson(\n",
    "                    ciclovia,\n",
    "                    name='geojson',\n",
    "                    style_function=lambda x:style_3\n",
    "                ).add_to(m)\n",
    "\n",
    "        except ValueError:\n",
    "            pass\n",
    "    \n",
    "    m.save(str(current_year)+\".html\")\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### CASE STUDY: Marvila"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "gdf_freguesias = gpd.read_file(\"Data/Lisboa_Freguesias/Lisboa_Freguesias_CAOP2015_TM06.shp\")\n",
    "gdf_freguesias = gdf_freguesias.to_crs(epsg=4326)\n",
    "gdf_freguesias.geometry.index = gdf_freguesias['DICOFRE']\n",
    "style_0 = {'fillColor': '#D6A86E', 'color': '#D6A86E'}\n",
    "style_1 = {'fillColor': '#1ED674', 'color': '#1ED674'}\n",
    "style_2 = {'fillColor': '#5CC1F7', 'color': '#5CC1F7'}\n",
    "style_3 = {'fillColor': '#472EAB', 'color': '#472EAB'}\n",
    "style_4 = {'fillColor': '#AB2E6E', 'color': '#AB2E6E'}\n",
    "\n",
    "unique_zones = gdf_freguesias['DICOFRE'].unique()\n",
    "m = folium.Map([38.749906, -9.110679],\n",
    "                 zoom_start=14,\n",
    "                 tiles='cartodbpositron')\n",
    "\n",
    "for i in range(len(ciclovias['features'])):\n",
    "    ciclovia = ciclovias['features'][i]\n",
    "    \n",
    "    try:\n",
    "        ano = int(ciclovia['properties']['ANO'])\n",
    "        if ano <= 2016:\n",
    "            folium.GeoJson(\n",
    "                ciclovia,\n",
    "                name='geojson',\n",
    "                style_function=lambda x:style_3\n",
    "            ).add_to(m)\n",
    "        if ano == 2017:\n",
    "            folium.GeoJson(\n",
    "                ciclovia,\n",
    "                name='geojson',\n",
    "                style_function=lambda x:style_1\n",
    "            ).add_to(m)\n",
    "        if ano == 2018:\n",
    "            folium.GeoJson(\n",
    "                ciclovia,\n",
    "                name='geojson',\n",
    "                style_function=lambda x:style_2\n",
    "            ).add_to(m)\n",
    "    except ValueError:\n",
    "        pass\n",
    "    \n",
    "for i, zone in enumerate(unique_zones):\n",
    "    if zone != '110621':\n",
    "        continue\n",
    "    mask = gdf_freguesias['DICOFRE'] == zone\n",
    "\n",
    "    example_or = gdf_freguesias.loc[mask]\n",
    "    folium.GeoJson(example_or,name='polygon_or',style_function=lambda x:style_0).add_to(m)\n",
    "\n",
    "m"
   ]
  }
 ],
 "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}