deconvolution.uh {hydromad} | R Documentation |
Estimates the unit hydrograph using deconvolution.
deconvolution.uh(P, Q, FWHM = length(P), do.plot = FALSE)
P |
Rainfall or effective rainfall time series |
Q |
Flow time series |
FWHM |
Full Width Half Maximum to use for Gaussian apodisation function |
do.plot |
Whether to plot the unit hydrograph immediately |
Estimates the unit hydrograph by fourier transform deconvolution of the ratio of the cross correlation and auto-correlation functions.
Vector of length (length(P)-1)*2
representing deconvolved unit
hydrograph. Note that entries from length(P)
represent are
obtained by negative lags.
Optionally plot the unit hydrograph.
Joseph Guillaume
Barry
data(Murrindindi) h <- deconvolution.uh(Murrindindi$P,Murrindindi$Q,do.plot=T) head(h)