Activité

Couverture territoriale et temporelle des stations RSQAQ

Comment décrire la répartition du réseau de stations?

2-3 heures · Introduction à intermédiaire

À produire : Des comparaisons par région et décennie d’ouverture, avec une carte du répertoire de stations distincte d’une carte de pollution.

Prérequis : Savoir ouvrir un projet RStudio et exécuter un script R. Savoir lire un tableau groupé et un graphique descriptif; repérer les valeurs manquantes.

Préparation enseignante

Télécharger la trousse ci-dessous. Extraire le ZIP, puis double-cliquer sur Donnees-bleues.Rproj. Avant la séance, ouvrir installer-packages.R et cliquer sur Source avec une connexion Internet.

Une fois les packages installés, les scripts utilisent les fichiers inclus et peuvent être exécutés hors ligne.

Télécharger la trousse avec les données

17 Ko. Préparation : 05/09/2026 à 07:04 (heure du Québec).

Les fichiers de données sont inclus dans la trousse.

Conditions de la source · Provenance et empreinte du ZIP

Fichiers inclus
  • resume_decennies_ouverture_rsqaq_stations.csv : 6 lignes, 2 colonnes.
  • resume_regions_rsqaq_stations.csv : 16 lignes, 8 colonnes.
  • rsqaq_stations.csv : 245 lignes, 13 colonnes.
  • valeurs_manquantes_rsqaq_stations.csv : 20 lignes, 4 colonnes.

Contexte

Cette activité utilise la fiche Stations du réseau de qualité de l’air du Québec. Unité : Une station du répertoire du réseau québécois de surveillance de la qualité de l’air.

Source : Données Québec / MELCCFP / RSQAQ. Les dates de préparation, les fichiers et les conditions de réutilisation sont indiqués dans la trousse.

Consignes

  1. Ouvrir activite-longue.R dans le sous-dossier datasets/qualite-air du projet RStudio et cliquer sur Source. Les tableaux apparaissent dans Console et les graphiques dans Plots.
  2. Comparer les régions et les décennies d’ouverture.
  3. Lire la carte des stations sans la présenter comme une carte de pollution.
  4. Rédiger une réponse à la question en nommant l’unité observée et une limite de la source.

Limites à faire nommer

Un répertoire de stations ne contient pas de concentrations. L’absence de date de fermeture ne certifie pas un fonctionnement actuel.

Vérifier le travail

  • L’unité, la période ou la portée documentaire et les variables utilisées sont nommées.
  • La conclusion répond à la question et explique au moins une limite précise.
Objectifs et adaptations

Adaptations

  • Débuter avec le script exécuté par la personne enseignante, puis faire interpréter une seule sortie.
  • Pour approfondir, modifier une sélection ou un regroupement et justifier ce que cela change.

Statut documenté : Documentation de source et ressources de classe fournies.

Contribution pédagogique : Aurélien Nicosia. Usage en cours non documenté.

Ajout
Mise à jour

Code complet

Ce script complet est celui de la trousse. Il lit les fichiers locaux et affiche les résultats dans RStudio.

Afficher le script R 
# Couverture territoriale et temporelle des stations RSQAQ
# Ouvrir le projet RStudio de l’archive, puis exécuter ce script.
# Source : https://www.donneesquebec.ca/recherche/dataset/rsqaq-stations
# Les données de classe sont figées; ce script ne les télécharge pas.

# Load libraries
library(dplyr)
library(ggplot2)
library(readr)

# Import prepared files
stations <- read_csv(
  "data/processed/qualite-air/rsqaq_stations.csv",
  show_col_types = FALSE
)

resume_regions <- read_csv(
  "data/processed/qualite-air/resume_regions_rsqaq_stations.csv",
  show_col_types = FALSE
)

resume_decennies <- read_csv(
  "data/processed/qualite-air/resume_decennies_ouverture_rsqaq_stations.csv",
  show_col_types = FALSE
)

valeurs_manquantes <- read_csv(
  "data/processed/qualite-air/valeurs_manquantes_rsqaq_stations.csv",
  show_col_types = FALSE
)

# Regions with the most stations without a closing date
print(resume_regions |>
  arrange(desc(n_sans_date_fermeture), desc(n_stations)) |>
  select(region_administrative, n_stations, n_sans_date_fermeture, n_fermees) |>
  head(10))

# Openings by decade
print(ggplot(resume_decennies, aes(x = decennie_ouverture, y = n_stations)) +
  geom_col(fill = "#325EA8") +
  labs(
    x = "Décennie d'ouverture",
    y = "Nombre de stations",
    title = "Ouvertures de stations RSQAQ par décennie"
  ))

# Simple point map
print(ggplot(stations, aes(x = longitude, y = latitude, color = station_sans_date_fermeture)) +
  geom_point(alpha = 0.75) +
  coord_fixed() +
  labs(
    x = "Longitude",
    y = "Latitude",
    color = "Sans date de fermeture",
    title = "Localisation des stations RSQAQ publiées",
    subtitle = "Carte de stations, pas carte de pollution"
  ))