Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
dyplot.default.R 4.57 KiB
1 
dyplot.default <- function(x, Qsup = NULL, Qsup.name = "Qsup",
2 
                           col.Precip = c("royalblue", "lightblue"), col.Q = c("black", "orangered", "grey"), col.na = "lightgrey",
3 
                           ylab = NULL, main = NULL,
4 
                           plot.na = TRUE, RangeSelector = TRUE, Roller = FALSE,
5 
                           LegendShow = c("follow", "auto", "always", "onmouseover", "never"), ...) {
6
7 
  # barChartPrecip <- scan(file = system.file("plugins/barChartPrecip.js", package = "airGRteaching"),
8 
  #                        what = "character", quiet = TRUE)
9
10 
  if (! any(class(x) %in% c("PrepGR", "CalGR", "SimGR"))) {
11 
    stop("Non convenient data for x argument. Must be of class \"PrepGR\", \"CalGR\" or \"SimGR\"")
12 
  }
13
14 
  if (is.null(ylab)) {
15 
    yunit <- .TypeModelGR(x)$TimeUnit
16 
    ylab  <- paste0(c("flow [mm/", "precip. [mm/"), yunit, "]")
17 
  } else {
18 
    if (length(ylab) < 2) {
19 
      ylab <- c(ylab, "")
20 
    }
21 
  }
22
23 
  if (is.null(Qsup)) {
24 
    Qsup <- as.numeric(rep(NA, length.out = length(x$Qobs)))
25 
  }
26 
  if (!is.numeric(Qsup)) {
27 
    stop("'Qsup' must be numeric")
28 
  }
29 
  if (length(Qsup) != length(x$Qobs)) {
30 
    stop("Incorrect length of 'Qsup', must be of length ", length(x$Qobs))
31 
  }
32 
  if (!is.character(Qsup.name)) {
33 
    Qsup.name <- as.character(Qsup.name)
34 
  }
35
36
37 
  if (any(class(x) %in% "PrepGR")) {
38 
    data <- data.frame(DatesR = x$InputsModel$DatesR,
39 
                       Precip = x$InputsModel$Precip,
40 
                       Qobs   = x$Qobs,
41 
                       Qsim   = NA,
42 
                       Qsup   = Qsup)
43 
    if (grepl("CemaNeige", x$TypeModel)) {
44 
      data$Psol <- rowMeans(as.data.frame(x$InputsModel$LayerPrecip) * as.data.frame(x$InputsModel$LayerFracSolidPrecip), na.rm = TRUE)
45 
      data$Pliq <- data$Precip - data$Psol
46 
      data$Precip <- NULL
47 
    }
48 
  } else {
49 
    data <- data.frame(DatesR = x$OutputsModel$DatesR,
50 
                       Precip = x$OutputsModel$Precip,
51 
                       Qobs   = x$Qobs,
52 
                       Qsim   = x$OutputsModel$Qsim,
53 
                       Qsup   = Qsup)
54 
    if (grepl("CemaNeige", x$TypeModel)) {
55 
      data$Psol <- rowMeans(sapply(x$OutputsModel$CemaNeigeLayers, function(x) x$Psol))
56 
      data$Pliq <- rowMeans(sapply(x$OutputsModel$CemaNeigeLayers, function(x) x$Pliq))
57 
      data$Precip <- NULL
58 
    }
59 
  }
60 
  data.xts <- xts::xts(data[, -1L], order.by = data$DatesR, tz = "UTC")
61
62
63 
  rgba <- function(x, alpha = 1) {
64 
    sprintf("rgba(%s, %f)", paste0(col2rgb(x), collapse = ", "), alpha)
65 
  }
66 
  if (length(col.Q) == 1) {
67 
    col.Q <- c(rgba(col.Q), rgba(col.Q, alpha = 0.5), rgba(col.Q, alpha = 0.3))
68 
  }
69 
  if (length(col.Q) == 2) {
70 
    col.Q <- c(rgba(col.Q[1L]), rgba(col.Q[2L]), rgba(col.Q[2L], alpha = 0.5))
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116
} if (length(col.Precip) < 2) { col.Precip <- c(rgba(col.Precip), rgba(col.Precip, alpha = 0.5)) } if (grepl("CemaNeige", x$TypeModel)) { Plim <- c(-1e-3, max(data$Psol+data$Pliq, na.rm = TRUE)) } else { Plim <- c(-1e-3, max(data$Precip, na.rm = TRUE)) col.Precip <- col.Precip[1L] } dg <- dygraphs::dygraph(data.xts, main = main, ...) dg <- dygraphs::dySeries(dygraph = dg, name = "Qobs", axis = "y", color = col.Q[1L], drawPoints = TRUE) dg <- dygraphs::dySeries(dygraph = dg, name = "Qsim", axis = "y", color = col.Q[2L]) dg <- dygraphs::dySeries(dygraph = dg, name = "Qsup", axis = "y", color = col.Q[3L], label = Qsup.name, strokePattern = "dashed") dg <- dygraphs::dyStackedBarGroup(dygraph = dg, name = rev(grep("^P", colnames(data.xts), value = TRUE)), axis = "y2", color = (col.Precip)) dg <- dygraphs::dyAxis(dygraph = dg, name = "y" , label = ylab[1L], valueRange = range(data.xts[, grep("^Q", colnames(data.xts))], na.rm = TRUE) * c(0.01, 1.59)) dg <- dygraphs::dyAxis(dygraph = dg, name = "y2", label = ylab[2L], independentTicks = FALSE, valueRange = rev(Plim) * c(2.99, 0.01)) if (RangeSelector) { dg <- dygraphs::dyRangeSelector(dygraph = dg, height = 15) } if (plot.na) { idNA <- .StartStop(data$Qobs, FUN = is.na) dg <- .DyShadingMulti(dygraph = dg, color = col.na, ts = data$DatesR, idStart = idNA$start, IdStop = idNA$stop) } if (Roller) { dg <- dygraphs::dyRoller(dygraph = dg, rollPeriod = 5) } if (is.numeric(Roller)) { dg <- dygraphs::dyRoller(dygraph = dg, rollPeriod = Roller) } if (any(LegendShow %in% c("follow", "auto", "always", "onmouseover", "never"))) { dg <- dygraphs::dyLegend(dygraph = dg, show = LegendShow[1L]) } dg <- dygraphs::dyOptions(dygraph = dg, useDataTimezone = TRUE) return(dg) }