Zones statistiques

Nous allons maintenant rapatrier les données statistiques infra-communales (les IRIS français) de l'INSEE sur nos zones isochrones.

Pour cela il nous faudra d'abord ramener les IRIS dans nos zones isochrones. Téléchargez le shape des IRIS français sur le site de L'IGN.

Fusion

Avant de nous attaquer au rapatriement des données IGN puis INSEE, fusionnez d'abord vos 308 zones isochrones dans un seul shape (Vector/Data Management Tools/Merge Vector Layers).

Une fois fait, vous pouvez supprimer vos 308 GEOJSON de votre projet QGIS. Supprimez également les champs inutiles de votre shape de sortie, nous ne garderons que l'id_mag bien sûr, mais aussi les champs area et total_pop (fournis par le service ORS), afin de faire quelques vérifications plus tard.

Découpe

Découpez maintenant la couche des IRIS sur vos polygones isochroniques. Attention, nous devons  conserver les doublons dans le cas des IRIS appartenant à plusieurs zones isochrones !

Utilisez donc l'outil Vector/Geoprocessing Tools/Intersection, vérifiez la présence de vos IRIS en doublons dans la table attributaire, mais aussi en faisant quelques tests de sélection sur votre couche de sortie à partir de l'id des magasins (en jouant sur la transparence des couches, la ville de Rennes est idéale pour ces vérifications).

PS : Plus tard, nous verrons comment automatiser ce type de sélection en Python. Ce sera très pratique pour visualiser proprement des zones isochrones qui se superposent.

Exportez un shape de votre couche de sortie, ici nommé iso_iris.

Rationalisation

Bien, disposant du CODE_IRIS, nous pourrions d'ores-et-déjà ramener les informations statistiques infra-communales de l'INSEE sur notre nouveau shape.

Cependant tel quel, la pertinence géographique serait discutable car dans bien des cas seules des portions très congrues d'IRIS ont été découpées. Or si un grand IRIS abrite 5000 personnes, mais que seuls quelques mètres carrés appartiennent à une zone isochrone, nous ne pouvons pas compter ces 5000 personnes dans notre zone.

Nous allons donc calculer un pourcentage surfacique pour chacun des IRIS compris dans nos zones isochrones. À terme il fournira un quotient nous permettant d'extraire des valeurs statistiques supposément plus fiables.

1. Commençons par calculer la surface de nos portions d'IRIS compris dans nos zones isochrones, dans un nouveau champ nommé part, avec le Field Calculator, fonction $area. Soyons généreux sur ses décimales, car certaines portions d'IRIS sont particulièrement petites. Enregistrez vos modifications.
2. Faisons-en de même avec les IRIS complets, dans un champ nommé area. Enregistrez.
3. Jointure depuis la couche iso_iris vers les IRIS, pour ramener le champ area nouvellement créé. Bien entendu nous utiliserons le champ CODE_IRIS.
4. Calcul du rapport entre nos 2 aires dans un nouveau champ nommé quotient. Enregistrement.

"part" / "iris_area"

Bien, nous pouvons maintenant supprimer la jointure (ou directement la couche des IRIS) ainsi que le champ part, nous disposons de notre quotient de rationalisation des futures données INSEE que nous allons rapatrier.

Liens ou pièces jointes

isochrones15m.zip	[Les 308 zones isochrones fusionnées]	979 Ko
Code complet de création des zones isochrones avec le service ORS	[code Python]	2 Ko
isochrones_ORS.zip	[Les 308 fichiers GEOJSON des zones isochrones]	786 Ko
Export OSM des magasins de sports français le 27 décembre 2020	[key=shop ; value=sport ; in=France]	0 Ko