Closed sdat2 closed 3 years ago
I will probably initially use https://github.com/sdat2/scikit-extremes to do the analysis. This only has Block-Maxima currently implemented in it (i.e I haven't managed to add Peak-over-Threshold).
Work Stream Goals
[x] Understanding the seasonality of extreme flood events
[ ] Building an extreme value distribution
[ ] Determine whether there are flood basins (gauged) and their thresholds - so we can know
[ ] Create a multivariate extreme distribution between rainfall and the river height (generally in R)
Block maxima fit for river height is terrible!
useful looking package: https://github.com/iagolemos1/thresholdmodeling/tree/0.0.1
Used 3m threshold alongside stage 60 min data (filtered). The non-independence of different values contributes to quite a bad fit; the top values all come from the large flood in 2007. news article on floods
Looking at the different peaks that come out after using scipy.find_peaks
, more reasonable as it doesn't pick up as many from the 2007 flood.
Project closed.
there aren't many python packages for doing this