{
"cells": [
{
"cell_type": "markdown",
"id": "fa4563c984470740",
"metadata": {},
"source": [
"# Isochrones\n",
"\n",
"This notebook builds isochrones on an intermodal `UrbanGraph` from IduEdu (walk + public transport):\n",
"\n",
"- `get_graph_isochrones` with `geometry_type=\"radius\"` and `\"ways\"` for one origin;\n",
"- `get_stepped_graph_isochrones` with bands of width `step` for `\"radius\"`, `\"ways\"` and `\"separate\"`.\n",
"\n",
"The method is described on the [isochrones page](../isochrones.rst). The interactive maps show a light copy of every\n",
"result (`for_map`: outlines simplified by 2 m), which keeps the notebook small.\n",
"\n",
"The last section, *Your own area*, shows what to replace to get the same pictures for another territory."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "9e7d002c63864de6",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:07:42.447914Z",
"start_time": "2026-07-10T10:07:42.444986Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:15:18.878933Z",
"iopub.status.busy": "2026-11-06T16:15:18.878933Z",
"iopub.status.idle": "2026-11-06T16:15:18.882451Z",
"shell.execute_reply": "2026-11-06T16:15:18.882451Z"
}
},
"outputs": [],
"source": [
"# Install required packages (uncomment if needed)\n",
"# !pip install objectnat iduedu"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "ebe2025fab286b4e",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:07:50.570317Z",
"start_time": "2026-07-10T10:07:42.622718Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:15:18.883963Z",
"iopub.status.busy": "2026-11-06T16:15:18.883963Z",
"iopub.status.idle": "2026-11-06T16:15:21.480888Z",
"shell.execute_reply": "2026-11-06T16:15:21.480888Z"
}
},
"outputs": [],
"source": [
"import geopandas as gpd\n",
"import iduedu\n",
"from iduedu import get_4326_boundary, get_intermodal_graph\n",
"from shapely import Point\n",
"\n",
"from objectnat import get_graph_isochrones, get_stepped_graph_isochrones"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "0ab5bbe7",
"metadata": {
"execution": {
"iopub.execute_input": "2026-11-06T16:15:21.482995Z",
"iopub.status.busy": "2026-11-06T16:15:21.482995Z",
"iopub.status.idle": "2026-11-06T16:15:21.644937Z",
"shell.execute_reply": "2026-11-06T16:15:21.644937Z"
}
},
"outputs": [],
"source": [
"import folium\n",
"\n",
"\n",
"def for_map(gdf, tolerance=2.0):\n",
" \"\"\"A light copy for a web map: outlines simplified by `tolerance` metres, coordinates rounded to about 10 cm.\"\"\"\n",
" local = gdf.to_crs(gdf.estimate_utm_crs())\n",
" local[\"geometry\"] = local.simplify(tolerance)\n",
" light = local.to_crs(4326)\n",
" light[\"geometry\"] = light.geometry.set_precision(1e-6)\n",
" return light[~light.is_empty]\n",
"\n",
"def grey_tiles():\n",
" \"\"\"Light grey base map (Esri World Light Gray); its tiles end at zoom 16 and are stretched beyond it.\"\"\"\n",
" return folium.TileLayer(\"Esri.WorldGrayCanvas\", max_zoom=19, max_native_zoom=16)"
]
},
{
"cell_type": "markdown",
"id": "494425a561aa50ec",
"metadata": {},
"source": [
"## 1. Intermodal graph\n",
"\n",
"The graph covers the district and is clipped to its boundary. Edge weights are `time_min` (minutes) and\n",
"`length_meter` (metres)."
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "df20dec8ac67c8da",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:07:54.215395Z",
"start_time": "2026-07-10T10:07:50.572833Z"
},
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-11-06T16:15:21.647417Z",
"iopub.status.busy": "2026-11-06T16:15:21.647417Z",
"iopub.status.idle": "2026-11-06T16:15:27.547355Z",
"shell.execute_reply": "2026-11-06T16:15:27.547355Z"
}
},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "d20e99012c9048329fe1744c56f5f4f3",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"Parsing ground PT routes: 0%| | 0/57 [00:00Make this Notebook Trusted to load map: File -> Trust Notebook
"
],
"text/plain": [
""
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"isochrones_radius = get_graph_isochrones(\n",
" G_intermodal,\n",
" geometry_type='radius',\n",
" gdf_origins=point,\n",
" weight_type=\"time_min\",\n",
" weight_value_cutoff=10,\n",
")\n",
"\n",
"# Visualize\n",
"for_map(isochrones_radius).explore(tiles=grey_tiles())"
]
},
{
"cell_type": "markdown",
"id": "1e9a3026b7537291",
"metadata": {},
"source": [
"## 4. Ways isochrone\n",
"\n",
"Pedestrian edges reachable within 10 minutes, buffered by `road_buffer_size` (5 m by default)."
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "b4bd257c0bafea8a",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:08:03.633384Z",
"start_time": "2026-07-10T10:08:02.727624Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:15:33.282448Z",
"iopub.status.busy": "2026-11-06T16:15:33.281431Z",
"iopub.status.idle": "2026-11-06T16:15:33.883577Z",
"shell.execute_reply": "2026-11-06T16:15:33.883577Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"Make this Notebook Trusted to load map: File -> Trust Notebook
Make this Notebook Trusted to load map: File -> Trust Notebook
Make this Notebook Trusted to load map: File -> Trust Notebook
Make this Notebook Trusted to load map: File -> Trust Notebook
Make this Notebook Trusted to load map: File -> Trust Notebook
Make this Notebook Trusted to load map: File -> Trust Notebook