05-osm_data_road_length.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import osmium as osm\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Total road length in Lisbon over time"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "class RoadLengthHandler(osm.SimpleHandler):\n",
    "    def __init__(self):\n",
    "        super(RoadLengthHandler, self).__init__()\n",
    "        self.length_road = 0.0\n",
    "        self.length = 0.0\n",
    "        self.length_cycle = 0.0\n",
    "        \n",
    "    def way(self, w):\n",
    "        # print(w)\n",
    "        if 'highway' in w.tags:\n",
    "            try:\n",
    "                self.length += osm.geom.haversine_distance(w.nodes)\n",
    "            except osm.InvalidLocationError:\n",
    "                # A location error might occur if the osm file is an extract\n",
    "                # where nodes of ways near the boundary are missing.\n",
    "                #print(\"WARNING: way %d incomplete. Ignoring.\" % w.id)\n",
    "                pass\n",
    "            \n",
    "            \n",
    "            if w.tags['highway'] == 'motorway':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'trunk':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'primary':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'secondary':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'terciary':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'motorway_link':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'trunk_link':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'primary_link':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'secondary_link':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'tertiary_link':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'road':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'residential':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'unclassified':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'living_street':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "            if w.tags['highway'] == 'busway':\n",
    "                try:\n",
    "                    self.length_road += osm.geom.haversine_distance(w.nodes)\n",
    "                except osm.InvalidLocationError:\n",
    "                    pass\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(2013, 2021):\n",
    "    h = RoadLengthHandler()\n",
    "    year = str(i)\n",
    "    \n",
    "    h.apply_file(\"Data/lisbon_\"+ year + \".osm.pbf\", locations=True, idx='flex_mem')\n",
    "    # print(year + ' Total way length: %.2f km' % (h.length/1000))\n",
    "    print(year + ' Total way length_road: %.2f km' % (h.length_road/1000))"
   ]
  }
 ],
 "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
}