Aux quatre coins de l'hexagone
Cartographier des points portant des valeurs (catégorie ou quantité) peut être un challenge : comment obtenir une représentation parlante du phénomène, simplement et rapidement interprétable par l'oeil. La quantité d'informations portée par une carte par points est souvent importante et nécessite une synthèse visuelle pour une plus grande efficacité de lecture.
Cette synthèse peut être réalisée en regroupant les points (clustering), en interpolant une surface continue dans l'espace (carte de chaleur par exemple) ou en répartissant l'information ponctuelle dans un maillage plus ou moins régulier de l'espace :
- Découpage administratif : communes, secteurs commerciaux, IRIS, ...
- Tessellation de l'espace : c'est à dire son partitionnement en éléments réguliers, le plus souvent carrés (carroyage).
Dans ce cas, une alternative au carroyage est le découpage en hexagones. L'hexagone a clairement un avantage esthétique sur le carroyage :
- La carte est plus "légère" visuellement.
- Les bordures externes de la carte sont "adoucies" car les angles ont perdu 30°.
- L'alignement des bords en lignes droites est moins systématique.
- L'alignement des bords en lignes droites est moins systématique.
Pour réaliser le partitionnement en hexagones, il n'y pas de fonction standard dans ArcGIS et de nombreux scripts sont disponibles sur en ligne. Le plus efficace que j'ai trouvé est disponible sur ArcGIS Online en tant que Geoprocessing Package prêt à l'emploi dans ArcGIS for Desktop 10.1. L'auteur est Tim Whiteaker de l'Université du Texas à Austin. Vous pouvez le télécharger à l'adresse suivante : http://bit.ly/13lMzLY.
Ce géotraitement a fait l'objet d'un article il y a quelques semaines sur le blog ArcGIS : http://shar.es/w56xu. La création des hexagones se fait en générant des polygones de Thiessen sur un semis régulier de points (les centres des hexagones), c'est clairement une approche performante du problème sur de grands espaces.
L'exemple : offre Vélib à Paris
Pousuivant sur le précédent article, l'exemple de mise en oeuvre de l'outil est la cartographie de l'offre Vélib à Paris en utilisant les données fournies en Open Data par JC Decaux, notamment la carte des stations Vélib.
Cartographier l'offre Vélib revient, avec les données fournies, à cartographier la densité, par hexagone (maillage de superficie régulière), de la somme des emplacements Vélib pour toutes les stations incluses dans un hexagone.
Un des intérêts de la cartographie d'un phénomène ponctuel à travers une tessellation de l'espace c'est bien la régularité du découpage : on peut passer directement de la valeur quantitative (nombre d'emplacement, d'habitants, ...) à la densité car toute les superficies sont égales, évitant ainsi les hérésies de sémiologie graphique que sont les cartes choroplèthes de phénomènes exprimés en valeur absolue.
![]() |
Stations Vélib par capacité d'emplacements - JC Decaux 2013 |
Un des intérêts de la cartographie d'un phénomène ponctuel à travers une tessellation de l'espace c'est bien la régularité du découpage : on peut passer directement de la valeur quantitative (nombre d'emplacement, d'habitants, ...) à la densité car toute les superficies sont égales, évitant ainsi les hérésies de sémiologie graphique que sont les cartes choroplèthes de phénomènes exprimés en valeur absolue.
Une fois récupérés :
- L'outils de géotraitement : http://bit.ly/13lMzLY.
- La couche des stations Vélib téléchargeable sous forme de paquetage de couche sur ArcGIS Online : http://bit.ly/17jrh61
La marche à suivre est classique :
- Découpage hexagonal sur l'emprise des stations
- Identité entre hexagones et points pour "remonter" sur les stations l'identifiant de l'hexagone de rattachement
- Récapitulation statistique pour sommer les emplacements Velib (portés par les stations) par hexagone
- Jointure de la table résultante avec la couche des hexagones
- Cartographie du résultat
Tout cela a été implémenté sous forme d'un modèle de géotraitement, disponible dans le paquetage de géotraitement suivant sur ArcGIS Online : http://bit.ly/1aVfZSu. Ce paquetage contient aussi l'outil de création d'hexagones et les données Vélib illustrant l'exemple. Attention, un paquetage de géotraitement n'est utilisable qu'en ArcGIS 10.1.
La carte finale présente donc la densité d'emplacements Vélib par hexagone et est, je pense, d'une lisibilité supérieure à la carte par points proportionnels présentée plus haut.
La carte finale présente donc la densité d'emplacements Vélib par hexagone et est, je pense, d'une lisibilité supérieure à la carte par points proportionnels présentée plus haut.
![]() |
Densité d'emplacement Vélib / hexagone - JC Decaux 2013 |