J’ai mis à jour le blog post d’accidents de vélo pour la ville de Lévis en ajoutant les données de 2017 et en corrigeant quelques bogues d’encodage.

Objectif

Le but de ce notebook R est de déterminer quels sont les endroits où se sont produits le plus d’accidents impliquant des vélos dans la ville de Lévis afin de déterminer à quels endroit des interventions seraient les plus bénéfiques.

Code

Le programme qui a servi à générer ce document est situé au https://raw.githubusercontent.com/SimonCoulombe/saaqmtq/master/ville_type.Rmd

Sources de données

Les données concernant les accidents proviennent du portail des données ouvertes du Québec, où la SAAQ a rendu disponible les données concernant tous les accidents rapportés à la police entre 2011 et 2017. Le dictionnaire de données est disponible en ligne.

Voici un échantillon des données telles qu’elles sont reçues du site web.

Pour effectuer la conversion des codes de municipalité vers le nom des municipalités, j’ai créé un tsv à partir du tableau présenté sur le site du mamrot.

Géocodage

Les données d’accident fournies par la SAAQ ne sont pas géocodées. L’information concernant l’accident est répartie entre plusieurs colonnes:

  • “NO_CIVIQ_ACCDN” , the street civic number
  • “SFX_NO_CIVIQ_ACCDN”, a suffix to the street number
  • “RUE_ACCDN”, the road name
  • “CD_MUNCP”, the city code. Here is a dictionnary to convert city code to name.

  • “NO_ROUTE” is the road number where the accident happened (numbered roads are typically highways). This seems to be used as a alternative to the road name RUE_ACCDN.
  • “CD_PNT_CDRNL_ROUTE” is the direction (North, South, East, West (Ouest) ) travelled on the road/highway.
  • “BORNE_KM_ACCDN” is the milestone number (used on highways and northern gravel roads)

They also use landmarks (intersection, etc..) to help locate the accident:

  • TP_REPRR_ACCDN is the type of landmark.
  • 1 means the intersection of two roads,
  • 2 " means “other landmark”
  • 0 means the type is not specified.
  • “ACCDN_PRES_DE” is the landmark that the type refers to. It can be the road that intersects the road named under “RUE_ACCDN”, a bridge, a school name, etc.
  • “NB_METRE_DIST_ACCD” is the distance in meters between the landmark and the accident.
  • “CD_PNT_CDRNL_REPRR” is the direction (North,South, East, Ouest) from the landmark to the accident.

Aussi, plusieurs abbréviations (BD, ST, STE, RTE, MGR, N, S, O, E, PRDE, etc) sont utilisées, ce qui complique le géocodage par l’API de google. Lors de la préparation des données, j’ai créé une variable “location”, où j’essaie de créer une chaîne de caractère que google pourra facilement géocoder.

TODO: Je n’ai rien fait pour les bornes kilométriques + numéro de route, mais ce ne devrait pas exclure trop d’accidents de vélos.

Pour effectuer le géocodage, j’utilise l’API de google maps via le package ggmap. La version gratuite de l’API est limitée à 50 requêtes par minute et 2 500 requêtes gratuites par jour.

Il s’est produit 170accidents impliquant un vélos dans ville de Lévis entre 2011 et 2017.

Je n’ai pas tenté le géocodage pour 4 car la localisation est trop ambiguë. Ces accidents sont présentés en appendice.

Parmi les 166 accidents que j’ai tenté de géocoder, l’API n’a pas réussi à retourner de latitude/longitude pour 3. Ces accidents sont présentés en appendice.

Le reste de l’analyse portera sur les 152 accidents qui ont été géocodés avec succès.

Résultats - données non-géocodées

Selon l’année

Selon le mois de l’année

Selon l’environnement

Selon les conditions météos

Selon la localisation transversale

Selon la semaine de l’année

Selon l’heure de la journée

## # A tibble: 7 x 2
##   CD_POSI_ACCDN                           n
##   <fct>                               <int>
## 1 Voie lente/voie de dépassement          6
## 2 Voie cyclable/chaussée désignée        14
## 3 "Voie de circulation "                117
## 4 Accotement (ou bord de la chaussée)     6
## 5 Trottoir                                8
## 6 Autre                                   7
## 7 <NA>                                   12

Selon la région administrative

Selon la ville

Selon la gravité des accidents

Liste des accidents mortels

Résultats - données géocodées

Top 10 des intersections où se sont produits le plus d’accidents

Le tableau ci-bas montre le top 10 des intersections où se sont produits le plus d’accidents impliquant des vélos dans la ville de Lévis entre 2011 et 2017.

Le top 10 des 10 villes comptant le plus d’accidents est présenté juste en dessous.

Liste des accidents survenus aux intersections les plus dangeureuses

Ci bas, la liste (sur plusieurs pages) de tous les accidents survenus à ces 10 intersections.

Résultats - cartes interactives

Carte des accidents mortels

Cette carte montre un cercle pour chacun des 0accidents mortels impliquant des vélos dans la ville de Lévis survenu entre 2011-2017 et géolocalisés avec succès. Lorsque plusieurs accidents arrivent au même endroit, on n’en voit qu’un seul.

Heatmap

La heatmap suivante permet de mieux voir les régions où plusieurs accidents se sont produits, même si les cordonnées latitude-longitude ne sont pas identiques. On retrouve que les intersections mentionnées précédemment sont en rouge, mais aussi le boulevard charest et langelier en général.

marker cluster du top 10

MarkerCluster

La dernière carte interactive montre des clusters. Si l’on zoom au maximum, on pourra par exemple voir les 6 accidents qui se sont produits à l’intersection de la 3e avenue et de la 4e rue, mais aussi un accident qui s’est produit juste à côté au 410 3e avenue.

Appendice

Le programme n’a pas essayé de géocoder les accidents suivants:

Google n’a pas réussi à géocoder les accidents suivants:

Google a codé ces accidents au centre de la ville:

google c codé ces addidents très loin du centre de la ville: