Beam to Earth coordinate transformations and other utilities for Acoustic Doppler Current Profiler (ADCP) data.
Convert raw along-beam velocity (positive toward transducer) time series (ntimes x nbins) to Earth-referenced velocity:
[u, v, w, w5] = janus5beam2earth(head, ptch, roll, theta, b1, b2, b3, b4, b5)
or, with all optional arguments specified,
[u, v, w, w5] = janus5beam2earth(head, ptch, roll, theta, b1, b2, b3, b4, b5, ...
'uvwBeam5', true, 'Gimbaled', true, 'Binmap', 'linear', r, r5)
Where
-
b[1-5]are the along-beam velocity time series [ntimes x nbins] (positive toward transducer face). -
u,v,ware the Earth-referenced (eastward, northward, upward) velocity time series [ntimes x nbins]. -
w5is the vertical velocity calculated from the vertical beam only. -
head,ptchandrollare the heading (rotation about the z-axis (beam 5), positive clockwise), pitch (rotation about the x-axis (beam 1), positive clockwise) and roll (rotation about the y-axis (beam 3), positive clockwise) angles measured by the ADCP, in degrees. -
thetais the beam angle measured from beam 5's direction (vertical if the instrument is level), in degrees. -
'uvwBeam5'(true/false) indicates whether to use the vertical beam data to calculateu,vandw. -
'Gimbaled'(true/false) indicates whether the pitch and roll sensors were gimbaled (mounted on a free-swiveling platform as opposed to a rigid frame). -
'Binmap'('none'/'nn'/'linear') indicates whether or not to perform bin-mapping on the raw beam velocities. Bin-mapping here means interpolating the velocities from all beams to the same horizontal plane before converting to instrument and then Earth coordinates. Defaults to'none', i.e., the assumption that pitch and roll are zero at all times.'linear'and'nn'specify linear and nearest-neighbor interpolation, respectively. -
randr5are vectors with the along-beam positions of the center of each bin for the Janus (#1-4) and vertical (#5) beams. Required only if'Binmap'is notnone.