Edit on GitHub

Data

Load the data

Loading the annual gaugings of the Sainte-Anne river:

data = RatingCurves.dataset("50408")
first(data,5)

5 rows × 3 columns

DateLevelDischarge
Date…Float64Float64
11965-05-1328.465236.681
21965-05-2027.959146.701
31965-08-2027.21641.481
41965-10-1427.57883.24
51966-05-2528.599261.782

The gaugings can be shown as by plotting the discharges as function of the levels:

set_default_plot_size(12cm, 8cm)
plot(data, x=:Level, y=:Discharge, Geom.point)
Level 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 20.0 20.5 21.0 21.5 22.0 22.5 23.0 23.5 24.0 24.5 25.0 25.5 26.0 26.5 27.0 27.5 28.0 28.5 29.0 29.5 30.0 30.5 31.0 31.5 32.0 32.5 33.0 33.5 34.0 34.5 35.0 35.5 36.0 36.5 37.0 37.5 38.0 20 25 30 35 40 20.0 20.2 20.4 20.6 20.8 21.0 21.2 21.4 21.6 21.8 22.0 22.2 22.4 22.6 22.8 23.0 23.2 23.4 23.6 23.8 24.0 24.2 24.4 24.6 24.8 25.0 25.2 25.4 25.6 25.8 26.0 26.2 26.4 26.6 26.8 27.0 27.2 27.4 27.6 27.8 28.0 28.2 28.4 28.6 28.8 29.0 29.2 29.4 29.6 29.8 30.0 30.2 30.4 30.6 30.8 31.0 31.2 31.4 31.6 31.8 32.0 32.2 32.4 32.6 32.8 33.0 33.2 33.4 33.6 33.8 34.0 34.2 34.4 34.6 34.8 35.0 35.2 35.4 35.6 35.8 36.0 36.2 36.4 36.6 36.8 37.0 37.2 37.4 37.6 37.8 38.0 h,j,k,l,arrows,drag to pan i,o,+,-,scroll,shift-drag to zoom r,dbl-click to reset c for coordinates ? for help ? -1200 -1000 -800 -600 -400 -200 0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200 -1000 -900 -800 -700 -600 -500 -400 -300 -200 -100 0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 -1000 0 1000 2000 -1000 -950 -900 -850 -800 -750 -700 -650 -600 -550 -500 -450 -400 -350 -300 -250 -200 -150 -100 -50 0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1050 1100 1150 1200 1250 1300 1350 1400 1450 1500 1550 1600 1650 1700 1750 1800 1850 1900 1950 2000 Discharge

Construct the Gauging type

It is possible to construct a Gauging type with a level and discharge couple:

julia> h = data.Level[1]28.465
julia> q = data.Discharge[1]236.681
julia> G = Gauging(h,q)Gauging(28.465, 236.681)

The level and the discharge of the gauging G can be retrieved with level and discharge methods respectively:

julia> h = level(G)28.465
julia> q = discharge(G)236.681

The constructor can be broadcasted to obtain a vector of type Gauging:

julia> G = Gauging.(data.Level, data.Discharge)194-element Vector{Gauging}:
 Gauging(28.465, 236.681)
 Gauging(27.959, 146.701)
 Gauging(27.216, 41.481)
 Gauging(27.578, 83.24)
 Gauging(28.599, 261.782)
 Gauging(27.011, 21.718)
 Gauging(27.094, 25.478)
 Gauging(27.216, 39.005)
 Gauging(29.883, 568.914)
 Gauging(29.334, 418.551)
 ⋮
 Gauging(28.182, 179.309)
 Gauging(26.964, 15.863)
 Gauging(27.193, 34.025)
 Gauging(27.374, 52.74)
 Gauging(26.874, 12.109)
 Gauging(29.186, 364.835)
 Gauging(29.697, 479.014)
 Gauging(27.793, 116.077)
 Gauging(28.471, 230.759)

The level and discharge methods can also be broadcasted. For example:

julia> discharge.(G)194-element Vector{Float64}:
 236.681
 146.701
  41.481
  83.24
 261.782
  21.718
  25.478
  39.005
 568.914
 418.551
   ⋮
 179.309
  15.863
  34.025
  52.74
  12.109
 364.835
 479.014
 116.077
 230.759