Secteurs municipaux et carte

La première partie du projet (et la plus longue) est de regarder la carte des secteurs municipaux de l’enquête et de créer la table de correspondance entre les secteurs de recensements (et sous-divisions de recensements) et les différents secteur municipaux visibles sur cette carte.

On va se taper la tâche laborieuse de recréer les secteurs municipaux utilisés lors de l’enquête et basés sur les census tract et census subdivisions du Recensement de 2016.

Le tableau de correspondance est ci-dessous, et mes notes quant à la création du tableau dans QGIS sont à la toute fin du post.

Les shapefiles sont sur le site de statistiques canada.

En fusionnant les tables de correspondance aux shapefiles, on réussit à reproduire la carte des secteurs municipaux de l’enquête. La seule exception est Sainte-Hélène de Breakeyville que j’ai dû fusionner à St-Jean Chrystome. Ils ont dû le découper à l’aide des aires de diffusions, mais disons que j’ai pas trouvé que ça valait la peine…

J’ai créé ma définition d’est et Ouest. Sur la Rive-Sud, j’ai ajouté à l’ouest St-Romuald Ouest (à l’Ouest de la rivière Etchemins) ainsi que St-Jean Chrysostome. Sur la Rive-Nord, j’ai ajouté Les Saules, Duberger, et Neufchâtel. Ça donne la carte ci-dessous et ça ne me semble pas excessif.

Trajets d’auto du matin

On commence par télécharger les matrices de déplacement entre les secteurs municipaux et on importe les trajets d’auto de la pointe du matin, soit l’onglet “SAXAU17T” de la feuille des matrices origine destination. On combine les trajets des secteurs 490 et 510, soit St-Jean-Chrysostome et St-Hélène de Breakeyville car notre carte ne discerne pas entre les deux. Aussi, on retire les trajets qui se terminent à l’extérieur du territoire de l’étude (SM = 999).

Ce sont donc 24 705 voitures qui traversent du sud vers le nord et 9 866 qui font le trajet inverse. Honnêtement, je suis surpris de l’ampleur du trafic vers le sud. 9.4 % des utilisateurs de voiture traversent un pont dans un sens ou l’autre le matin.

Bon, on va maintenant regarder quelle proportion de nos 34 571 trajets qui traversent le fleuve vont véritablement de l’est vers l’est.
La réponse c’est 13,6 %. On est loin de 40% ! et c’est très dépendant de l’hypothèse qui veut que Lebourgneuf et Saint-Romuald Est sont prêts à aller jouer dans le bout de l’île d’orléans.

Un troisième lien?

On va regarder un peu d’où viennent et où vont les gens qui traversent les ponts.

La carte suivante montre la destination des gens qui prennent leur auto le matin. On voit clairement les 4 pôles que sont Ste-Foy, Lebourgneuf,la colline parlementaire et Lévis-Centre.

La carte suivante montre la destination des gens qui prennent leur voiture et qui traversent dans un sens ou dans l’autre. Sans surprise vu la distance, Lebourgneuf et la colline parlementaire sont moins attractifs pour les gens qui doivent traverser le fleuve, tandis que Ste-Foy et Lévis-Centre demeurent très attractif.

La présence importante de Lévis-Centre chez les traverseurs de fleuve m’a un peu surpris. J’ai donc sorti le pourcentage de personnes qui arrivent dans une région qui ont dû traverser le fleuve. C’est 23% à Lévis, 18% à Ste-Foy, 13% sur la colline et 6% à Lebourgneuf.

QGIS

Mes notes QGIS. Je suis très très débutant, je l’utilise seulement quand je n’ai pas vraiment le choix :)

• Comme d’habitude, on charge le shapefile de secteurs de recensement avec “Layer - Add Vector Layer”,
  
• Pour ajouter une couche de fond Open Street Maps, on suit les instructions ici pour "no plugin required, right click “Tile Server XYZ” - Add Connection et on met l’URL suivant: http://tile.openstreetmap.org/{z}/{x}/{y}.png .

• On va travailler beaucoup avec l’ “identify tool” et le select tool
• Pour pouvoir ajouter la colonne “SM” aux attributs de la couche de SR, je commence par “toggle editing”, en faisait right click la couche dans “layers panel” puis “toggle editing”. Ensuite, pour ajouter la colonne, je reclique sur la couche dans “layers panel”, puis je fais “Properties”, “Fields” et je clique sur le petit icône “Field Calculator” pour créer le champ SM avec une valeur (expression) de 0. Note: le screenshot pointe vers “new field”, il faut plutôt cliquer sur field calculator, qui est le 4e icone à droite de new field.

* Ensuite, le jeu c’est de rendre la couleur des polygones dépendante de la SM, afin de pouvoir voir les progrès que je fais à mesure que j’attribue une SM à un polygones. Pour ce faire, je dois donner un style “graduated” à la layer de secteur de recensement et je dis de me baser sur la colonne “SM”. J’active le “live update” afin que la carte soit mise à jour à chaque fois que je change le SM d’un polygone et je dois aussi cliquer “classify” sporadiquement afin qu’il attribue des couleurs à chaque valeur.

* Enfin, pour modifier la valeur de SM, j’utilise le identify tool pour mettre à jour le champ SM ( et je clique Classify quand j’entre un nouveau maximum.)

Une fois que j’ai tout bien peinturé la carte, j’ai sauvegardé le shapefile, je l’ai lu et j’ai sauvegardé seulement la colonne “SM” et “SDRIDU”, le résultat est ci-dessous.