Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
#extract.R# 3.92 KiB
1 
library(here)
2 
library(tools)
3 
source(file.path(dirname(here()), 'init.R'))
4
5
6 
# Extraction of data
7 
extract_data = function (filedir, filename) {
8
9 
    data_path = file.path(computer_data_path, filedir, filename)
10
11 
    data = read.table(data_path,
12 
                      header=TRUE,
13 
                      na.strings=c('-99.000', '-99'),
14 
                      sep=';',
15 
                      skip=41)[,1:2]
16
17 
    return (data=data)
18 
}
19
20 
# data = extract_data('', 'H5920011_HYDRO_QJM.txt')
21
22
23 
# Extraction of information
24 
extract_info = function (filedir, filename) {
25
26 
    data_path = file.path(computer_data_path, filedir, filename)
27
28 
    # General information on station
29 
    iStatut = c('0'='inconnu',
30 
                '1'='station avec signification hydrologique',
31 
                '2'='station sans signification hydrologique',
32 
                '3'="station d'essai")
33
34 
    iFinalite = c('0'='inconnue',
35 
                  '1'="hydrométrie générale",
36 
                  '2'='alerte de crue',
37 
                  '3'="hydrométrie générale et alerte de crue",
38 
                  '4'="gestion d'ouvrage",
39 
                  '5'='police des eaux',
40 
                  '6'="suivi d'étiage",
41 
                  '7'='bassin expérimental',
42 
                  '8'='drainage')
43
44 
    iType = c('0'='inconnu',
45 
              '1'='une échelle',
46 
              '2'='deux échelles, station mère',
47 
              '3'='deux échelles, station fille',
48 
              '4'='débits mesurés',
49 
              '5'='virtuelle')
50
51 
    iInfluence = c('0'='inconnue',
52 
                   '1'='nulle ou faible',
53 
                   '2'='en étiage seulement',
54 
                   '3'='forte en toute saison')
55
56 
    iDebit = c('0'='reconstitué',
57 
               '1'="réel (prise en compte de l'eau rajoutée ou retirée du bassin selon aménagements)",
58 
               '2'='naturel')
59
60 
    iQBE = c('0'='qualité basses eaux inconnue',
61 
             '1'='qualité basses eaux bonne',
62 
             '2'='qualité basses eaux douteuse')
63
64 
    iQME = c('0'='qualité moyennes eaux inconnue',
65 
             '1'='qualité moyennes eaux bonne',
66 
             '2'='qualité moyennes eaux douteuse')
67
68 
    iQHE = c('0'='qualité hautes eaux inconnue',
69 
             '1'='qualité hautes eaux bonne',
70 
             '2'='qualité hautes eaux douteuse')
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119
data_path = file.path(computer_data_path, filedir, filename) infotxt = c(readLines(data_path, n=41)) info = list(code=trimws(substr(infotxt[11], 38, nchar(infotxt[11]))), nom=trimws(substr(infotxt[12], 39, nchar(infotxt[12]))), territoire=trimws(substr(infotxt[13], 39, nchar(infotxt[13]))), L93X=as.numeric(substr(infotxt[16], 38, 50)), L93Y=as.numeric(substr(infotxt[16], 52, 63)), surface=as.numeric(substr(infotxt[19], 38, 50)), statut=iStatut[trimws(substr(infotxt[26], 38, 50))], finalite=iFinalite[trimws(substr(infotxt[26], 52, 56))], type=iType[trimws(substr(infotxt[26], 58, 58))], influence=iInfluence[trimws(substr(infotxt[26], 60, 60))], debit=iDebit[trimws(substr(infotxt[26], 62, 62))], QBE=iQBE[trimws(substr(infotxt[26], 72, 72))], QME=iQME[trimws(substr(infotxt[26], 74, 74))], QHE=iQHE[trimws(substr(infotxt[26], 76, 76))] ) return (info=info) } # info = extract_info('', 'H5920011_HYDRO_QJM.txt') # Extraction and dataframing extract_df = function (filedir, filename) { if (filename == 'all') { filelist = c() for (f in list.files(file.path(filedir))) { if (file_ext(f) == 'txt') { filelist = c(filelist, f) } } } data_path = file.path(computer_data_path, filedir, filename) info = extract_info(filedir, filename) df = data.frame(info, data_path=data_path) return (df) } # df = extract_df('', 'H5920011_HYDRO_QJM.txt')