Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
format.R 3.39 KiB
1 
# Usefull library
2 
library(dplyr)
3
4
5 
join = function (df_data_BH, df_data_NV, df_meta_BH, df_meta_NV) {
6
7 
    if (!is.null(df_data_NV) & !is.null(df_data_BH)) {
8
9 
        # Get the station in common
10 
        common = levels(factor(df_meta_NV[df_meta_NV$code %in% df_meta_BH$code,]$code))
11 
        # Get the Nv station to add
12 
        NVadd = levels(factor(df_meta_NV[!(df_meta_NV$code %in% df_meta_BH$code),]$code))
13
14 
        # Select only the NV meta to add
15 
        df_meta_NVadd = df_meta_NV[df_meta_NV$code %in% NVadd,]
16 
        # Join NV data to BH data
17 
        df_meta = full_join(df_meta_BH, df_meta_NVadd, by=c("code", "nom", "L93X", "L93Y", "surface_km2", "file_path"))
18
19 
        # Select only the NV data to add
20 
        df_data_NVadd = df_data_NV[df_data_NV$code %in% NVadd,]
21 
        # Join NV meta to BH meta
22 
        df_data = full_join(df_data_BH, df_data_NVadd, by=c("Date", "Qm3s", "code"))
23
24 
    } else if (is.null(df_data_NV) & !is.null(df_data_BH)) {
25 
        df_meta = df_meta_BH
26 
        df_data = df_data_BH
27 
        df_data_NVadd = NULL
28 
    } else if (!is.null(df_data_NV) & is.null(df_data_BH)) {
29 
        df_meta = df_meta_NV
30 
        df_data = df_data_NV
31 
        df_data_NVadd = df_data_NV
32 
    } else {
33 
        stop('No data')
34 
    }
35
36 
    return (list(data=df_data, meta=df_meta, NVadd=df_data_NVadd))
37 
}
38
39
40 
prepare = function(df_data, colnamegroup=NULL) {
41
42 
    colnamegroup = c(colnamegroup)
43 
    colindgroup = which(colnames(df_data) == colnamegroup)
44 
    df_data = group_by_at(df_data, colindgroup)
45
46 
    data = tibble(Date=df_data$Date,
47 
                  group=group_indices(df_data),
48 
                  Qm3s=df_data$Qm3s)
49 
    Gkey = group_keys(df_data)
50
51 
    info = bind_cols(group=seq(1:nrow(Gkey)),
52 
                     Gkey)
53 
    return (list(data=data, info=info))
54 
}
55
56
57 
reprepare = function(df_XEx, df_Xlist, colnamegroup=NULL) {
58
59 
    colnames(df_XEx) = c('Date', 'group', 'Qm3s')
60
61 
    df_XEx$Date = as.character(df_XEx$Date)
62 
    exDate = df_XEx$Date[1]
63 
    nbt = lengths(regmatches(exDate, gregexpr('-', exDate)))
64
65 
    if (nbt == 1) {
66 
         df_XEx$Date = paste(df_XEx$Date, '01', sep='-')
67 
    } else if (nbt == 0) {
68 
         df_XEx$Date = paste(df_XEx$Date, '01', '01', sep='-')
69 
    } else if (nbt != 2) {
70 
        stop('erreur of date format')
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110
} df_XEx = bind_cols(Date=as.Date(df_XEx$Date, format="%Y-%m-%d"), df_XEx[-1], df_Xlist$info[df_XEx$group, 2:ncol(df_Xlist$info)]) df_XlistEx = prepare(df_XEx, colnamegroup=colnamegroup) return (df_XlistEx) } clean = function (df_Xtrend, df_XEx, df_Xlist) { # print(str(df_XEx)) df_Xlist = reprepare(df_XEx, df_Xlist, colnamegroup=c('code')) df_Xlist$data$code = NA for (g in df_Xlist$info$group) { df_Xlist$data$code[which(df_Xlist$data$group == g)] = df_Xlist$info$code[df_Xlist$info$group == g] } # df_Xlist$data = df_Xlist$data[, !names(df_Xlist$data) == "group")] df_Xtrend = bind_cols(df_Xtrend, df_Xlist$info[df_Xtrend$group1, 2:ncol(df_Xlist$info)]) colnames(df_Xtrend)[1] = 'group' intercept = get_intercept(df_Xtrend, df_Xlist, unit2day=365.25) df_Xtrend$intercept = intercept df_Xtrend = relocate(df_Xtrend, intercept, .after=trend) return (list(trend=df_Xtrend, data=df_Xlist$data, info=df_Xlist$info)) }