Difference between revisions of "Outdoor air temperature in Finland"

## This is code is Op_en2959/temperatures [[Outdoor air temperature in Finland]]
library(OpasnetUtils)

ta <- opbase.data("Op_en2959", subset = "Vantaa") # [[Outdoor air temperature in Finland]]
hours <- as.character(ta$Date)

hours <- as.POSIXct(strptime(hours, format = "%Y-%b-%d %H:%M:%S"))
hours <- hours + as.difftime(as.character(ta$Time), format = "%H:%M")

# Adjust temperature data: 2.3 C is the average difference between Kuopio and Helsinki
hours <- data.frame(Time = hours, Result = as.numeric(as.character(ta$Result)) + 2.3)
hours <- hours[!is.na(hours$Result) , ]
hours$Date <- as.POSIXct(strptime(format(hours$Time, "%Y-%m-%d"), format = "%Y-%m-%d"))
hours <- hours[hours$Time >= as.POSIXct("2014-03-01 00:00:00") & hours$Time < as.POSIXct("2015-03-01 00:00:00") , ]

days <- aggregate(hours["Result"], hours["Date"], FUN = mean)
days$Temperature <- cut(days$Result, breaks = seq(-30, 33, 3))
temperatures <- Ovariable("temperature", data = aggregate(days["Result"], days["Temperature"], FUN = mean))
temperdays <- Ovariable("temperdays", data = aggregate(days["Result"], days["Temperature"], FUN = length))

objects.store(temperatures, temperdays)
cat("Objects temperature, temperdays stored.\n")

#ggplot(hours, aes(x = Time, y = Result))+geom_line() + geom_point(data = days, aes(x = Date, y = Result, colour = "Daily mean", size = 2))
#ggplot(hours, aes(x = Time, y = Result))+geom_point()