{
"cells": [
{
"cell_type": "markdown",
"id": "4dc8b5b9755684aa",
"metadata": {},
"source": [
"# Coverage zones\n",
"\n",
"This notebook builds coverage zones of the example services on an intermodal `UrbanGraph`: every place is assigned\n",
"to the service it reaches at the lowest network cost within a budget, so zones do not overlap. It shows\n",
"\n",
"- `get_graph_coverage` by distance (800 m) and by time (10 min);\n",
"- `get_stepped_graph_coverage` in bands of 2 minutes, as Voronoi cells and as separate buffers;\n",
"- the graph-free `get_radius_coverage`.\n",
"\n",
"The method is described on the [coverage page](../coverage.rst). The interactive maps show a light copy of every result\n",
"(`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": "db56d2861ab407be",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:05:52.330251Z",
"start_time": "2026-07-10T10:05:52.325735Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:15:52.509926Z",
"iopub.status.busy": "2026-11-06T16:15:52.508921Z",
"iopub.status.idle": "2026-11-06T16:15:52.513095Z",
"shell.execute_reply": "2026-11-06T16:15:52.513095Z"
}
},
"outputs": [],
"source": [
"# Install required packages (uncomment if needed)\n",
"# !pip install iduedu pyarrow objectnat"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "d568e9b461c6086a",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:06:00.545968Z",
"start_time": "2026-07-10T10:05:52.332766Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:15:52.514099Z",
"iopub.status.busy": "2026-11-06T16:15:52.514099Z",
"iopub.status.idle": "2026-11-06T16:15:55.277734Z",
"shell.execute_reply": "2026-11-06T16:15:55.277734Z"
}
},
"outputs": [],
"source": [
"import geopandas as gpd\n",
"import iduedu\n",
"from iduedu import get_4326_boundary, get_intermodal_graph\n",
"\n",
"from objectnat import get_graph_coverage, get_radius_coverage, get_stepped_graph_coverage"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "9e307cbd",
"metadata": {
"execution": {
"iopub.execute_input": "2026-11-06T16:15:55.279796Z",
"iopub.status.busy": "2026-11-06T16:15:55.279796Z",
"iopub.status.idle": "2026-11-06T16:15:55.411469Z",
"shell.execute_reply": "2026-11-06T16:15:55.411469Z"
}
},
"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": "5084078ebb6beef2",
"metadata": {},
"source": [
"## 1. Intermodal graph"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "initial_id",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:06:04.759829Z",
"start_time": "2026-07-10T10:06:00.601794Z"
},
"collapsed": true,
"execution": {
"iopub.execute_input": "2026-11-06T16:15:55.413525Z",
"iopub.status.busy": "2026-11-06T16:15:55.413525Z",
"iopub.status.idle": "2026-11-06T16:16:01.228856Z",
"shell.execute_reply": "2026-11-06T16:16:01.228856Z"
}
},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "7f61903f300c47a29258b24e42288318",
"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": [
"# Calculate coverage zones by distance (800m cutoff)\n",
"result_length = get_graph_coverage(\n",
" G_intermodal,\n",
" gdf_destinations=services,\n",
" weight_type=\"length_meter\",\n",
" weight_value_cutoff=800,\n",
")\n",
"\n",
"# Coverage results are indexed like the input services. Join service attributes\n",
"# back by index before coloring or filtering by service columns such as \"name\".\n",
"result_length = result_length.join(services[['name']], how=\"left\")\n",
"\n",
"# Visualize results\n",
"for_map(result_length).explore(column='name', legend=False, tiles=grey_tiles())"
]
},
{
"cell_type": "markdown",
"id": "5f8d6888bdbf4716",
"metadata": {},
"source": [
"## 4. Coverage by time\n",
"\n",
"Zones within 10 minutes (`weight_type=\"time_min\"`), clipped to the district boundary with `zone`."
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "3fe99dad99d0c066",
"metadata": {
"ExecuteTime": {
"end_time": "2026-07-10T10:06:53.475685Z",
"start_time": "2026-07-10T10:06:52.513709Z"
},
"execution": {
"iopub.execute_input": "2026-11-06T16:16:06.695256Z",
"iopub.status.busy": "2026-11-06T16:16:06.695256Z",
"iopub.status.idle": "2026-11-06T16:16:07.523450Z",
"shell.execute_reply": "2026-11-06T16:16:07.523450Z"
}
},
"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