psd_plot: Particle size distribution plot

In sbegueria/disdRo: Functions for working with disdrometric data

Description

Produce a PSD plot: particle density (number of drops per m^3 of air and mm of rain) vs particle size.

Usage

psd_plot(x, type = "Thies", filter = NULL, a = NULL)

Arguments

x	A particle size velocity distribution (PSVD) matrix.
type	Character vector designing the type of disdrometer, currently one of 'Thies' or 'Parsivel'. Defaults to 'Thies'.
filter	A value between 0 and 1. Removes outlier bins, i.e. those that are between (1-value) and (1+value) far from the Bear theoretical fall velocity model. Defaults to Inf (no outliers are removed).
a	Optional, the sensor area (in cm2). Defaults to the mean sensor area of Thies and Parsivel disdrometers.

Value

A PSD plot.

References

NA

Examples

f <- system.file('extdata', 'rawDataThies', package='disdRo')
files <- list.files(f, '.txt', full.names=TRUE, recursive=TRUE)
dsd <- psvd_read(files, type='Thies')
day <- apply(dsd, c(2,3), sum)
psd_plot(day)
psd_plot(day, filter=psvd_filter(type='Thies', d=c(0.3,7), tau=0.5))

sbegueria/disdRo documentation built on May 14, 2019, 2:39 p.m.