How to calculate current speed and direction from ADCP velocity components
FollowThe raw velocity measurement is a vector in the direction along each of the beams, which is referred to as beam coordinates. Beam coordinates can be converted to a Cartesian coordinate system (XYZ) by knowing the beam orientation. Furthermore, the flow can be presented in Earth normal coordinates (ENU, East, North and Up). In order to get the information referenced to earth coordinates (ENU) it is therefore necessary to detect the instrument’s orientation in space. The preferred coordinate system may be specified in the deployment software.
For ADCP instruments the ASCII converted raw data for the velocities are given in the .v1, .v2 and .v3 files. Depending on your chosen coordinate system the content in each file is as follows:
.v1 is beam 1, or the X axis (in XYZ), or East (in ENU)
.v2 is beam 2, or the Y axis (in XYZ), or North (in ENU)
.v3 is beam 3, or the Z axis (in XYZ), or Up (in ENU)
In Beam coordinates, a positive velocity is in the same direction as the beam points. Beam 1 is marked with an “X” on the head. In XYZ coordinates, a positive velocity in the X-direction points in the direction of the X-axis arrow. The X-axis points in the same direction as Beam 1. In ENU coordinates, a positive east velocity points toward east. This is also a right-handed orthogonal system. The velocity data unit is m/s.
The speed, or magnitude of the velocity vector, for ENU components can be calculated as follows:
Speed = sqrt(Veast^2 + Vnorth^2 + Vup^2)
The angle, given in interval [-π, π]:
Angle = atan2(Veast/Vnorth)
To get the direction in earth coordinates [0,360] the following transformation can be done:
Direction = 180 + 180/pi*(atan2(-Veast,-Vnorth))
0 points north, 90 degrees points towards east, 180 south and 270 degrees towards west.
Comments
0 comments
Article is closed for comments.