meetdichtheidsmatrix toegevoegd voor waterstanden

08_bijlage.Rmd

@@ -6,7 +6,7 @@ knitr::opts_chunk$set(
 	cache = F
 )
 
-source("r/tex2Rtable.R")
+# source("r/tex2Rtable.R")
 source("r/functions.R")
 require(smwfs)
 require(stringr)
@@ -79,14 +79,62 @@ De meetdichtheidsmatrices geven voor ieder jaar aan hoe vaak elke parameter is g
 
 ## Hydro - Waterstanden
 
+```{r}
+frozendatapath <- file.path(datapath, paste0("Data_Hydro_waterstanden_all_", ".csv"))
+
+trendstations <- c("Europlatform","Vlakte van de Raan","Cadzand","Westkapelle","Breskens Veerhaven","Vlissingen","Gat van Borssele","Overloop van Hansweert", "Terneuzen", "Hansweert", "Walsoorden", "Baalhoek", "Schaar van de Noord", "Bath")
+
+df_waterstanden <- read_delim(frozendatapath, delim = ",", guess_max = 200000) %>%
+  mutate(stationname = ifelse(
+    stationname == "Euro Platform", "Europlatform", stationname)
+  ) %>%
+  filter(stationname %in% trendstations)
+
+```
+
+```{r, fig.height=10, fig.width=8}
+df_waterstanden %>%
+  filter(grepl("water", parametername, ignore.case = T)) %>%
+  ungroup() %>%
+  mutate(
+    year = year(datetime),
+    month = month(datetime)
+  ) %>%
+  group_by(stationname, year, month, parametername) %>%
+  summarise(datapoints = n(), .groups = "drop") %>%
+  mutate(date = as_date(paste(year, month, "15"))) %>%
+  meetdichtheid_hydro(stationName = trendstations, no_of_columns = 2)
+
+```
+
+
+
+
+
+
+
 ## Hydro - Golven
 
-```{r, fig.height=2, fig.width=6}
+```{r read-golven-md}
 trendstations <- c("Bath (BAT2)","Hansweert (HAWI)","Cadzand (CAWI)","Cadzand (CADW)","Overloop van Hansweert (OVHW)","Deurloo (DEUR)","Honte (HNTE)","Hoofdplaat (HFPL)","Overloop van Valkenisse (OVVA)","Pas van Terneuzen (PVT)","Wielingen (WIEL)")
 
 frozendatapath <- file.path(datapath, paste0("Data_Hydro_golven_all_", dataJaar, ".csv"))
+
 # df_golven <- read_delim(frozendatapath, delim = ",", guess_max = 200000) %>% 
-df_golven_md <- read_delim("https://watersysteemdata.deltares.nl/thredds/fileServer/watersysteemdata/Westerschelde/Scheldemonitor/2023/Data_Hydro_golven_all.csv", delim = ",", guess_max = 200000) %>% 
+df_golven <- read_delim(
+  "https://watersysteemdata.deltares.nl/thredds/fileServer/watersysteemdata/Westerschelde/Scheldemonitor/2023/Data_Hydro_golven_all.csv", 
+  delim = ",", 
+  guess_max = 200000, 
+  show_col_types = FALSE
+  )
+
+``` 
+
+
+```{r, fig.height=10, fig.width=8}
+
+# fix names
+df_golven_md <- df_golven %>% 
   filter(stationname %in% trendstations) %>%
   mutate(stationname = factor(stationname, levels = trendstations)) %>%
   mutate(value = case_when(
@@ -106,16 +154,19 @@ df_golven_md <- read_delim("https://watersysteemdata.deltares.nl/thredds/fileSer
     str_detect(parametername, "H3") ~ "H3: Gemiddelde hoogte van het 1/3 deel hoogste golven in cm",
     str_detect(parametername, "TM02") ~ "TM02: Golfperiode berekend uit het spectrum in s",
     str_detect(parametername, "Hm0") ~ "Hm0: Significante golfhoogte uit 10mHz spectrum in cm")) %>%
-  mutate(parametername = gsub(":.*","",parametername))
+  mutate(parametername = gsub(":.*","",parametername)) %>%
+  mutate(date = as.Date(datetime))
 
-lapply(trendstations, 
-        function(x){
-meetdichtheid_hydro(df_golven_md , stationName = x)          
-        })
+# plot
+meetdichtheid_hydro(df_golven_md, stationName = trendstations)
 
+# base::lapply(trendstations, 
+#         function(x){
+# meetdichtheid_hydro(df_golven_md , stationName = x)          
+#         })
 
 ```
-
+lapply()
 
 ## Oppervlaktewater {#BijlageB-meetdichtheid-oppwater}