Unit Hydrograph using deconvolution

Estimates the unit hydrograph using deconvolution.

deconvolution.uh(P, Q, FWHM = length(P), do.plot = FALSE)

Arguments

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

Value

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.

Details

Estimates the unit hydrograph by fourier transform deconvolution of the ratio of the cross correlation and auto-correlation functions.

References

Barry

See also

expuh,lambda,armax,powuh

Author

Joseph Guillaume

Examples

data(Murrindindi)
h <- deconvolution.uh(Murrindindi$P, Murrindindi$Q, do.plot = TRUE)

head(h)
#> [1] 0.18076783 0.29867732 0.15726599 0.08379355 0.05542525 0.05342593