Stib

import gtfs_kit as gk
import requests
import tempfile

url = "https://api.mobilitytwin.brussels/stib/gtfs"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).content

with tempfile.NamedTemporaryFile(delete=False, suffix=".zip") as f:
    f.write(data)

feed = gk.read_feed(f.name, 'm')

# Explore the STIB schedule
print(f"Routes: {len(feed.routes)}")
print(f"Stops: {len(feed.stops)}")
print(f"Trips: {len(feed.trips)}")
print("\nRoutes:")
print(feed.routes[['route_short_name', 'route_long_name', 'route_type']])

import requests
import tempfile

from gtfs_parquet import read_parquet

url = "https://api.mobilitytwin.brussels/stib/gtfs-parquet"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).content

with tempfile.NamedTemporaryFile(delete=False, suffix=".zip") as f:
    f.write(data)

feed = read_parquet(f.name)

# List all routes (Polars DataFrame — zero-copy, extremely low RAM usage)
print(feed.routes.select("route_short_name", "route_long_name"))

# Count trips per route
print(feed.trips.group_by("route_id").len().sort("len", descending=True))

import requests
import geopandas as gpd

url = "https://api.mobilitytwin.brussels/stib/stops"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

gdf = gpd.GeoDataFrame.from_features(data["features"])

# List stops for a specific line
line = "1"
line_stops = gdf[gdf['route_short_name'] == line].sort_values('stop_sequence')
print(f"Stops on line {line}:")
for _, stop in line_stops.iterrows():
    print(f"  {stop['stop_sequence']}. {stop['stop_name']}")

# Plot all stops on a map
gdf.plot(figsize=(10, 8), markersize=3, color="red")

import requests
import pandas as pd

url = "https://api.mobilitytwin.brussels/stib/vehicle-distance"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Convert to DataFrame for analysis
df = pd.DataFrame(data)

# Group by line and compute average distance from stop
avg_distance = df.groupby('lineId')['distanceFromPoint'].mean()
print("Average distance from stop per line:")
print(avg_distance.sort_values(ascending=False))

import requests
import geopandas as gpd

url = "https://api.mobilitytwin.brussels/stib/vehicle-position"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

gdf = gpd.GeoDataFrame.from_features(data["features"])

# Count vehicles per line
vehicles_per_line = gdf.groupby('lineId').size().sort_values(ascending=False)
print("Active vehicles per line:")
print(vehicles_per_line)

# Plot vehicles colored by line
gdf.plot(figsize=(10, 8), column='color', legend=True, markersize=8)

import requests
import pandas as pd

url = "https://api.mobilitytwin.brussels/stib/speed"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

df = pd.DataFrame(data)

# Average speed per line (km/h)
speed_per_line = df.groupby('lineId')['speed'].mean().sort_values()
print("Average speed per line (km/h):")
print(speed_per_line)

# Find the slowest segments
slowest = df.nsmallest(5, 'speed')[['lineId', 'pointId', 'speed']]
print("\nSlowest segments:")
print(slowest)

import requests
import pandas as pd

url = "https://api.mobilitytwin.brussels/stib/aggregated-speed"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

df = pd.DataFrame(data)

# Compare average speed across lines over 10-minute intervals
speed_by_line = df.groupby('lineId')['speed'].agg(['mean', 'min', 'max'])
print("Speed statistics per line (km/h):")
print(speed_by_line.sort_values('mean'))

import requests
import movingpandas as mpd

url = "https://api.mobilitytwin.brussels/stib/trips"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Load as a MovingPandas TrajectoryCollection
tc = mpd.io.read_mf_dict(data, traj_id_property="uuid")

# Compute speed and distance for each vehicle trajectory
for traj in tc.trajectories[:5]:
    print(f"Vehicle {traj.id}: {traj.get_length():.0f}m, duration: {traj.get_duration()}")

# Plot all vehicle trajectories
tc.plot(figsize=(12, 8), linewidth=0.5)

import requests
import geopandas as gpd

url = "https://api.mobilitytwin.brussels/stib/shapefile"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

gdf = gpd.GeoDataFrame.from_features(data["features"])

# Plot the full STIB network colored by line
gdf.plot(figsize=(12, 10), color=gdf['couleur_hex'], linewidth=1)

# List all unique lines
print("STIB lines:", sorted(gdf['ligne'].unique()))

import requests
import json

url = "https://api.mobilitytwin.brussels/stib/stop-details"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Parse and display stop details
for stop in data['results'][:10]:
    name = json.loads(stop['name'])
    coords = json.loads(stop['gpscoordinates'])
    print(f"Stop {stop['id']}: {name['fr']} / {name['nl']} ({coords['latitude']}, {coords['longitude']})")

import requests
import json

url = "https://api.mobilitytwin.brussels/stib/stops-by-line"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Display the route of a specific line
for route in data['results']:
    if route['lineid'] == '1':
        dest = json.loads(route['destination'])
        stops = json.loads(route['points'])
        print(f"Line {route['lineid']} → {dest['fr']} ({route['direction']})")
        print(f"  {len(stops)} stops: {stops[0]['id']} → {stops[-1]['id']}")

import requests
import json
from datetime import datetime

url = "https://api.mobilitytwin.brussels/stib/waiting-times"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Show next arrivals at each stop
for entry in data['results'][:10]:
    times = json.loads(entry['passingtimes'])
    for t in times:
        dest = t['destination']['fr']
        arrival = datetime.fromisoformat(t['expectedArrivalTime'])
        print(f"Stop {entry['pointid']} — Line {entry['lineid']} → {dest} at {arrival:%H:%M}")

import requests
import json

url = "https://api.mobilitytwin.brussels/stib/travellers-information"

data = requests.get(url, headers={
        'Authorization': 'Bearer [MY_API_KEY]'
}).json()

# Display active service alerts sorted by priority
for alert in sorted(data['results'], key=lambda x: x['priority'], reverse=True):
    content = json.loads(alert['content'])
    lines = json.loads(alert['lines'])
    text = content[0]['text'][0]['en']
    line_ids = ', '.join(l['id'] for l in lines)
    print(f"[Priority {alert['priority']}] Lines {line_ids}: {text}")