psvd_integrate: Compute PSVD integrated variables

In sbegueria/disdRo: Functions for working with disdrometric data

Description

dsd_integrate reads raw disdrometer data and computes a series of integrated variables.

Usage

psvd_integrate(indir, script = NULL, outfile = NULL, interp = "linear")

Arguments

indir	A character vector with the url of the directory tree that contains the raw disdrometric data.
script	A character vector with the url of the Perl script. Defaults to the installation directory of the disdRo package.
outfile	A character vector with the url of the output file. Defaults to a random file name in the session's temporary directory.
interp	A character vector indicating the interpolation method to use for inputing particle sizes and velocities within the bins of the PSVD matrix. One of 'middle', 'uniform', and 'linear', defaulting to 'middle'.

Details

Currently, this is done via an external Perl script, so you need to have Perl installed and working in your system. Beware: some users have reported issues for running the Perl script in Windows. It might be translated into a native R script in the future, if I find time to do it.

Value

A data frame with the following items:

type

Disdrometer type, currently one of 'Thies' or 'Parsivel' (Factor)

serial

Sensor serial number (Numeric)

time

Date and time of the record (POSIXct)

seconds

Number of seconds since 1970-01-01 00:00:00 (Numerica)

synop

Synop code (Factor)

r

Precipitation intensity computed from the PSVD matrix, mm h−1 (Numeric)

p

Precipitation amount computed from the PSVD matrix, mm (Numeric)

m

Liquid water content computed from the PSVD matrix, g m-3 (Numeric)

z

Radar reflectivity computed from the PSVD matrix, dB mm6 m−3 (Numeric)

e

Kinetic energy computed from the PSVD matrix, J m−2 mm−1 (Numeric)

mor

Visibility computed from the PSVD matrix, m (Numeric)

r_meas

Precipitation intensity as reported by the disdrometer, mm h−1 (Numeric)

z_meas

Radar reflectivity as reported by the disdrometer, dB mm6 m−3 (Numeric)

e_meas

Kinetic energy as reported by the disdrometer, J m−2 h−1 (Numeric)

mor_meas

Visibility as reported by the disdrometer, m (Numeric)

qual

Data quality reported by the distrometer, 0-100 (Numeric)

tmp

Air temperature, ºC (Numeric)

rh

Relative humidity, 0-100 (Numeric)

w

Wind speed, m s-1 (Numeric)

wd

Wind direction, degrees (Numeric)

np

Number of particles detected computed from the PSVD matrix (Numeric)

np_meas

Number of particles detected as reported by the disdrometer (Numeric)

lcurrent

Laser control output, 1/100 mA (Numeric)

ocontrol

Optical control output, mV (Numeric)

power

Sensor power supply, V (Numeric)

tmp_int

Internal sensor temperature, ºC (Numeric)

d10

Particle diameter’s 10th percentile, mm (Numeric)

d25

Particle diameter’s 25th percentile, mm (Numeric)

d50

Particle diameter’s 50th percentile, mm (Numeric)

d75

Particle diameter’s 75th percentile, mm (Numeric)

d90

Particle diameter’s 90th percentile, mm (Numeric)

dmean

Mean particle diameter, mm (Numeric)

v10

Particle velocity’s 10th percentile, m s−1 (Numeric)

v25

Particle velocity’s 25th percentile, m s−1 (Numeric)

v50

Particle velocity’s 50th percentile, m s−1 (Numeric)

v75

Particle velocity’s 75th percentile, m s−1 (Numeric)

v90

Particle velocity’s 90th percentile, m s−1 (Numeric)

vmean

Mean particle velocity, m s−1 (Numeric)

t_shift

Telegram shift time, s (Numeric)

nrow

Telegram number of rows (Numeric)

err

Error status (Numeric)

ncol

Number of fields in the telegram (Numeric)

Bin interpolation

Since the particle size and velocity distribution is not linear within the bins of the PSVD matrix, different imputation methods exist. 'middle' will assing the middle bin size and velocity to all the particles in the bin; 'uniform' assumes an uniform distribution of sizes and velocities within the bin limits; and 'linear' assumes a linear distribution between the bin limits.

Examples

f <- system.file('extdata/rawDataParsivel', package='disdRo')
x <- dsd_integrate(f)
head(x)
summary(x)

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