In mundl/smires: Science and Management of Intermittent Rivers and Ephemeral Streams

generated with r R.version.string and smires version r packageVersion("smires")

knitr::opts_chunk$set(echo = TRUE)

library(smires) 
smires

table(smires$country)

hasCoordinates <- smires %>% 
  select(lon, lat) %>% 
  complete.cases() %>% 
  table(dnn = NULL)

hasCoordinates

library(ggplot2)
basemap(smires) + 
  geom_point()

Interactive Map with basemap from Openstreet Maps

library(leaflet)

leaflet(smires) %>% 
  addTiles() %>% 
  addProviderTiles(providers$CartoDB.Positron) %>%
  addMarkers(lng = ~lon, lat = ~lat,
             clusterOptions = markerClusterOptions()
  )

Compute a few Metrics for all stations with coordinates

library(purrr)

metrics <- smires %>%
  filter(complete.cases(lon, lat)) %>%
  mutate(MAMD = map_dbl(data, MAMD),
         CVAMD = map_dbl(data, CVAMD))

metrics %>%
  filter(country == "fr") %>%
  basemap() + 
  geom_point(aes(size = MAMD, fill = CVAMD), shape = 21, alpha = 0.7) +
  scale_size_continuous("Mean annual\nmaximum duration\nin days") + 
  scale_fill_gradient2("Coefficient\nof Variation")

ggplot(metrics) + 
  geom_point(aes(x = CVAMD, y = MAMD, col = country)) +
  labs(x = "Coefficient of Variation", y = "Mean annual Maximum Duration in days")

library(tidyr)
is_nf_year <- function(x) 
  char_binary(x = x, 
              fun_major = function(x) sum(x) > 0,
              state = "no-flow", varname = "no_flow_year")

smires6 <- smires %>%
  filter(complete.cases(lon, lat)) %>%
  mutate(data = map(data, function(x) 
    filter(x, time >= as.Date("2000-01-01") & time < as.Date("2006-01-01")))) %>%
  filter(map_int(data, nrow) > 0) 

nf_year <- smires6 %>%
  mutate(nf = map(data, is_nf_year)) %>%
  unnest(nf)


library(dplyr)
nf_year  %>%
  basemap() + 
  geom_point(aes(col = no_flow_year)) +
  scale_color_discrete("is intermittent") + 
  facet_wrap(~major)

library(purrr)
library(smires)
s <- smires[c(1, 43), ]



m <- s %>%
  transmute(
    no_flow_years = map_dbl(data, no_flow_years),
    MAN = map_dbl(data, MAN),
    CVAN = map_dbl(data, CVAN),
    FAN = map(data, FAN),
    MAMD = map_dbl(data, MAMD),
    CVAMD = map_dbl(data, CVAMD),
    FAMD = map(data, FAMD),
    tau0 = map_dbl(data, tau0),
    tau0r = map_dbl(data, tau0r),
    tauE = map_dbl(data, tauE),
    tauEr = map_dbl(data, tauEr),
    tau = map_dbl(data, tau),
    taur = map_dbl(data, taur),
    k = map_dbl(data, k),
    ksd = map_dbl(data, ksd),
    nrv = map_dbl(data, nrv)
  ) %>%
  mutate_if(is.numeric, round, digits = 2)

m

mundl/smires documentation built on May 23, 2019, 8:22 a.m.