Comment trouver le bureau de poste le plus proche de chez vous?

Comment trouver le bureau de poste le plus proche de chez vous à l’aide de Voronoï et Python ?

Bien entendu, cette question peut se formuler avec autre chose qu’un bureau de poste…

Si vous ne connaissez pas le pavage de Voronoï, je vous encourage à regarder ce post.

Vous pouvez aussi regarder cet article sur Wikipedia.

Au menu sur cette page...

Comment trouver le bureau de poste le plus proche de chez vous? La méthode.

La méthode est “simple” : on répertorie toutes les coordonnées de ce que l’on cherche (ici, les bureaux de poste) et on fait notre pavage…

Facile à dire, mais concrètement ?

Récupération des coordonnées

On va s’aider de l’API OpenStreetMap.

import requests

def get_post_offices_overpass(city_name):
    query = f"""
    [out:json];
    area["name"="{city_name}"]->.searchArea;
    node["amenity"="post_office"](area.searchArea);
    out body;
    """
    url = "http://overpass-api.de/api/interpreter"
    response = requests.post(url, data={'data': query})
    data = response.json()
    post_offices = [(element['tags'].get('name', 'Post Office'), element['lat'], element['lon'])
                    for element in data['elements']]
    return post_offices

Placement sur une carte et tracé du pavage de Voronoï

On va s’aider de folium.

import folium
import numpy as np
from shapely.geometry import Polygon
from shapely.ops import unary_union
from geovoronoi import voronoi_regions_from_coords

def create_map(city_name="Bordeaux"):
    post_offices = get_post_offices_overpass(city_name)

    # Créer une carte centrée sur Bordeaux
    map_bordeaux = folium.Map(location=[44.8378, -0.5792], zoom_start=13)

    # Ajouter les bureaux de poste à la carte
    for name, lat, lng in post_offices:
        folium.Marker(location=[lat, lng], popup=name).add_to(map_bordeaux)
        
    # Convertir les coordonnées des bureaux de poste en objets Point
    coords = np.array([(lng, lat) for name, lat, lng in post_offices])

    # Définir une zone de délimitation (par exemple, les limites de Bordeaux)
    # Utiliser un polygone explicite pour Bordeaux
    boundary_coords = [(-0.8, 45.1), (0.02, 45.1), (0.02, 44.6), (-0.8, 44.6), (-0.8, 45.1)]
    boundary_polygon = Polygon(boundary_coords)

    # Calculer les régions de Voronoi à partir des coordonnées
    try:
        region_polys, region_pts = voronoi_regions_from_coords(coords, boundary_polygon)
    except RuntimeError as e:
        print(f"Erreur lors de la génération des régions de Voronoi: {e}")
        return

    # Ajouter les régions de Voronoi à la carte folium
    for region in region_polys.values():
        folium.Polygon(locations=[(lat, lng) for lng, lat in region.exterior.coords]).add_to(map_bordeaux)

    # Sauvegarder la carte dans un fichier HTML
    map_bordeaux.save('map_bordeaux.html')
    
    return None

Cela va générer une page HTML d’une carte zoomable:

Bon, là, ce n’est qu’une image car je ne peux pas insérer la page directement, mais vous voyez le genre ?

Et pour une autre ville ?

La seule chose qui diffère dans le code est le rectangle dans lequel se trouve la ville: il faudra le construire vous même (le boundary_polygon).

Pour se faire, vous pouvez ajouter l’instruction:

boundary_polygon = Polygon(boundary_coords)
print("Coords:", coords) # ligne à ajouter

tout de suite après avoir défini la variable boundary_polygon. Vous verrez ainsi où se situent les points et vous n’aurez plus qu’à modifier les coordonnées des coins du rectangle :

boundary_coords = [(xmin, ymax), (xmax, ymax), (xmax, ymin), (xmin, ymin), (xmin, ymax)]

Et pour autre chose qu’un bureau de poste?

Il faudra changer les informations de la fonction get_post_offices_overpass. Concrètement, il faudra changer “Post Office” par ce que vous voulez. Pour connaître les termes, consultez la page:

https://wiki.openstreetmap.org/wiki/FR:Key:shop