Winter started a bit more than one month ago… but we have already experienced many snow storms… there is still a lot snow in gardens and in the streets,
I was wondering if it was that unusual, but apparently not. Compared with last year, it is (for the first months of winter, until the end of Januray), it +50%, but it is comparable with previous years
Yes, we a simple loop, we can easily extract data from official wesite https://climat.meteo.gc.ca/ (but not too far away, even 2015 contains a lot of missing observations). For this month, we use
url = "https://climat.meteo.gc.ca/climate_data/daily_data_f.html?StationID=51157&timeframe=2&StartYear=1840&EndYear=2023&Day=30&Year=2023&Month=1#" library(XML) library(stringr) download.file(url,destfile = "M.html") tables=readHTMLTable("M.html") k = which(tables[[1]]$`JOUR `=="Somme") neige = tables[[1]]$`Neige tot. Definitioncm `[k] x = as.numeric(sub(",", ".", strsplit(neige, "LegendCarer")[[1]][1], fixed = TRUE))
and then we loop, and store the number we look for in a data frame (yes, we have to convert “50,8LegendCarer^” into the appropriate numerical value (that would be here 50.8…
D = data.frame(annee = c(2023,rep(2022:2015,each=12),c(12,11,10)), mois= c(1,rep(12:1,8),12,11,10), lab = neige, snow = x) for(i in 2:nrow(D)){ y = D$annee[i] m = D$mois[i] url = paste("https://climat.meteo.gc.ca/climate_data/daily_data_f.html?StationID=51157&timeframe=2&StartYear=1840&EndYear=2023&Day=30&Year=",y,"&Month=",m,"#",sep="") download.file(url,destfile = "M.html") tables=readHTMLTable("M.html") k = which(tables[[1]]$`JOUR `=="Somme") neige = tables[[1]]$`Neige tot. Definitioncm `[k] x = as.numeric(sub(",", ".", strsplit(neige, "LegendCarer")[[1]][1], fixed = TRUE)) D[i,3] = neige D[i,4] = x }
Here are the most recent months
> head(D) annee mois lab snpw 1 2023 1 50,8LegendCarer^ 50.8 2 2022 12 63,0 63.0 3 2022 11 14,6 14.6 4 2022 10 0,0 0.0 5 2022 9 0,0 0.0 6 2022 8 0,0 0.0
Of course, we need some codes to plot, we here, I mainly wanted to keep tracks of the code used to extract meteorological data…
OpenEdition suggests that you cite this post as follows:
Arthur Charpentier (January 29, 2023). Snow in Montréal (Canada). Freakonometrics. Retrieved December 3, 2024 from https://doi.org/10.58079/ovli
You should try the “weathercan” R package to get Environment Canada’s weather stats. It an awsome package! Let it snow!
(https://github.com/ropensci/weathercan).