本发明提供一种LADCP与组合惯性导航系统联合观测系统及方法，涉及海流观测领域，系统包括：LADCP、组合惯导系统主机、数据处理单元、无线传输模块、耐压无线天线、电池舱、第一GNSS天线、第二GNSS天线和仪器支架；所述数据处理单元通过串口与所述LADCP和组合惯导系统主机连接，得到海流剖面测量数据、GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率，并通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正。本发明提高LADCP流速和流向的测量精度，降低了使用LADCP进行剖面流观测的工作成本。
一种LADCP与组合惯性导航系统联合观测系统及方法
技术领域
[0001] 本发明属于海流观测领域，尤其涉及一种LADCP与组合惯性导航系统联合观测系统及方法。
背景技术
[0002] LADCP(Lowered Acoustic Doppler Current Profiler，下放式声学多普勒海流剖面仪)是20世纪90年代出现的一种海流剖面测量方法，能够快速获得大范围、全水深的海流剖面数据，在国内外的调查中已有了较广泛的应用。LADCP通常与CTD(Conductivity‑Temperature‑Depth profiler，温盐深仪)捆绑在一起，并与CTD水下单元一起升降，以一定的速度从海表下放到海底，再从海底提升到海面，从而获得一系列的单个流速小剖面。这些流速小剖面需要经过后期资料处理，重叠计算以获得整个测量深度上的流速剖面。
[0003] LADCP 直接测得的流速流向是海流相对于LADCP 的相对速度和方向，这个流速流向只有在LADCP 静止的情况下才是海流的真实速度。由于LADCP 采用投放式测量，钢缆又有一定的倾角且测量船不能保证绝对零速度，所以LADCP垂直速度和水平速度的量值很大且无规律,所以，如何从LADCP 的测量速度中提取出海流的真实速度一直是LADCP 测量技术中的难点问题。目前，LADCP 测量过程中通过底跟踪技术求出LADCP 仪器本身相对于海底的速度，再求出海流的速度。但底跟踪技术只有在近底层附近才有意义，而且通过底跟踪并不能准确的计算LADCP整个下放过程中的每个时刻的三维运动速度。在海洋走航观测工作中有时因为工作时间限制LADCP 并不能每次都能够下放到海底测到底跟踪数据，导致LADCP的测量误差的测量较大甚至失效。 LADCP测量流向使用的自身的磁罗盘，测量船作为一个巨大的铁磁物体在下放初期会对LADCP的流向测量产生较大的影响。
[0004] 现有的LADCP海流观测方法中，缺乏实时高精度校正流速流向观测数据的装置及方法，导致海流速度和方向观测误差较大，并且每个观测过程必须观测到底跟踪数据，致使海洋海流剖面测量工作效率较低，增加了海洋工作成本。
发明内容
[0005] 为解决上述技术问题，本发明提出一种LADCP与组合惯性导航系统联合观测系统及方法的技术方案，以解决上述技术问题。
[0006] 本发明第一方面，公开了一种LADCP与组合惯性导航系统联合观测系统，所述系统包括：LADCP、组合惯导系统主机、数据处理单元、无线传输模块、耐压无线天线、电池舱、第一GNSS天线、第二GNSS天线和仪器支架；
[0007] 所述LADCP，即下放式声学多普勒海流剖面仪，安装于所述仪器支架内部的下方，所述组合惯导系统主机、数据处理单元、无线传输模块和电池舱分别安装于所述仪器支架内部，且所述组合惯导系统主机、数据处理单元、无线传输模块和电池舱位于所述LADCP的上方，所述耐压无线天线安装于所述无线传输模块的上方，所述第一GNSS天线和第二GNSS天线安装于所述仪器支架的外部的上方；
[0008] 所述电池舱为所述LADCP、组合惯导系统主机、数据处理单元、无线传输模块、第一GNSS天线和第二GNSS天线供电；
[0009] 所述第一GNSS天线和第二GNSS天线接收GNSS卫星信号并通过同轴线缆传输至所述组合惯导系统主机；
[0010] 所述LADCP实时观测海流剖面测量数据；
[0011] 所述组合惯导系统主机观测GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率；
[0012] 所述数据处理单元通过串口与无线传输模块连接，所述无线传输模块通过所述耐压无线天线与传输数据处理单元和母船甲板单元之间的通信数据；
[0013] 所述数据处理单元通过串口与所述LADCP和组合惯导系统主机连接，得到海流剖面测量数据、GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率，并通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对所述LADCP观测的海流剖面的三维流速进行实时高精度校正。
[0014] 根据第一方面所述的系统，所述联合观测系统还包括：第一耐压仓、第二耐压仓和第三耐压仓；
[0015] 所述第一耐压仓和第二耐压仓位于所述仪器支架的外部的上方；
[0016] 所述第一GNSS天线安装于所述第一耐压仓内部；
[0017] 所述第二GNSS天线安装于所述第二耐压仓内部；
[0018] 所述第三耐压仓位于所述仪器支架内部的所述LADCP的上方，所述组合惯导系统主机、数据处理单元、无线传输模块和电池舱分别安装于所述第三耐压仓内部，所述耐压无线天线安装于所述第三耐压仓上方；
[0019] 所述第一耐压仓和第二耐压仓分别通过同轴水密电缆和水密接插件与所述第三耐压仓连接，所述LADCP通过水密电缆与所述第三耐压仓连接。
[0020] 根据第一方面所述的系统，所述LADCP的北方向与所述组合惯导系统主机的北方向对齐。
[0021] 本发明第二方面，公开了一种LADCP与组合惯性导航系统联合观测方法，采用第一方面任一项所述的一种LADCP与组合惯性导航系统联合观测系统，LADCP实时观测海流剖面的测量数据和所述组合惯导系统主机观测GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率的操作方法包括：
[0022] 观测船到达观测点位后，将所述LADCP的坐标系与所述组合惯导系统主机的坐标系对齐，即LADCP的北方向与所述组合惯导系统主机的北方向对齐；
[0023] 将所述联合观测系统吊至海面，通过母船甲板单元接收所述数据处理单元的回传数据；
[0024] 将所述联合观测系统下放至海下预设深度，观测海流剖面测量数据；
[0025] 观测海流剖面测量数据完成后，将所述联合观测系统收回。
[0026] 根据第二方面所述的方法，所述操作方法还包括：
[0027] 判断所述组合惯导系统主机的GNSS位置和GNSS定向数据确定后，将所述联合观测系统下放至海下预设深度。
[0028] 根据第二方面所述的方法，所述通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正的方法包括：
[0029] 以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率进行离散积分，得到航向角；
[0030] 应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，即姿态校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度；
[0031] 应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，即线速度校正，得到第二校正三维流速；
[0032] 应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0033] 根据第二方面所述的方法，
[0034] 所述以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率离散积分的方法，得到航向角之前，所述方法还包括：
[0035] 消除LADCP的噪声和外界环境引起的LADCP的航向角速率离散采样值的瞬时粗大误差，具体方法包括：
[0036] 计算航向角速率离散采样值序列的算数平均值；
[0037] 计算航向角速率离散采样值序列的标准差；
[0038] 应用所述算数平均值和标准差，计算航向角速率离散采样值序列的标志值，具体公式为：
[0039]
[0040] 其中， 表示航向角速率离散采样值序列中采样值t的标志值，s表示标准差，表示序列中第t个航向角速率离散采样值， 表示算数平均值；
[0041] 应用序列中航向角速率离散采样值对应的标志值与预设的临界值比较，如果小于所述临界值，保留所述航向角速率离散采样值，否则剔除所述航向角速率离散采样值；
[0042] 滤除航向角速率离散采样值的直流分量，具体方法包括：
[0043] 对消除了瞬时粗大误差的航向角速率离散采样值进行离散傅里叶变换；
[0044] 将频谱图中频率为0的频谱分量置0后进行离散傅里叶逆变换，得到滤除直流分量的航向角速率离散采样值。
[0045] 根据第二方面所述的方法，所述以入水前最后一组没有失效的GNSS定向数据作为初值，通过对所述航向角速率离散积分的方法，得到航向角的方法包括：
[0046] 所述以入水前最后一组没有失效的GNSS定向数据 作为初值，以最后一组没有失效的GNSS定向数据的采样时刻作为时间零点，通过对所述航向角速率离散积分的方法，求出预设时刻为止LADCP旋转的角度，对所述LADCP旋转的角度360取余，得到所述航向角。
[0047] 根据第二方面所述的方法，所述应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度的方法包括：
[0048]
[0049]
[0050]
[0051]
[0052]
[0053] 其中， 代表未校正的矢量测量结果， ，， ， 表示三维流
速， 表示三维线速度， 表示三维角速度； 代
表横滚角， 代表俯仰角， 代表航向角； 代表第一次校正后的矢量校
正结果， ， ，
， 表示第一校正三维流速，
表示第一校正三维线速度， 表示第一
校正三维角速度。
[0054] 根据第二方面所述的方法，所述应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，得到第二校正三维流速的方法包括：
[0055] 应用所述第一校正三维线速度与所述第一校正三维流速进行矢量相加，得到第二校正三维流速，具体公式为：
[0056]
[0057] 其中， 表示第二校正三维流速；
[0058] 应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速的方法包括：
[0059] 应用LADCP的波束出口至仪器支架的吊点的距离R与所述第一校正三维角速度相乘，得到三维旋转线速度，再应用所述三维旋转线速度与所述第二校正三维流速进行矢量相加，得到第三校正三维流速，具体公式为：
[0060]
[0061] 其中， 表示第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0062] 本发明提出的方案，提高LADCP流速和流向的测量精度，降低了使用LADCP进行剖面流观测的工作成本。
附图说明
[0063] 为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0064] 图1为根据本发明实施例的一种LADCP与组合惯性导航系统联合观测系统的结构图；
[0065] 图2为根据本发明实施例的LADCP与组合惯性导航系统联合观测方法流程图；
[0066] 图3为根据本发明实施例的所述通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正的方法流程图；
[0067] 图4为根据本发明实施例的消除LADCP的噪声和外界环境引起的LADCP的航向角速率离散采样值的瞬时粗大误差方法流程图。
[0068] 图中，1‑仪器支架，2‑第一耐压仓，3‑第二耐压仓，4‑耐压无线天线，5‑第三耐压仓，6‑ LADCP。
具体实施方式
[0069] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例只是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0070] 实施例1：
[0071] 本发明第一方面公开了一种LADCP与组合惯性导航系统联合观测系统，图1为根据本发明实施例的一种LADCP与组合惯性导航系统联合观测系统的结构图，具体如图1所示，所述系统包括：所述系统包括：LADCP6、组合惯导系统主机、数据处理单元、无线传输模块、耐压无线天线4、电池舱、第一GNSS天线、第二GNSS天线和仪器支架1；
[0072] 所述LADCP6，即下放式声学多普勒海流剖面仪，安装于所述仪器支架1内部的下方，所述组合惯导系统主机、数据处理单元、无线传输模块和电池舱安装于所述仪器支架1内部的所述LADCP6的上方，所述耐压无线天线4安装于所述无线传输模块的上方，所述第一GNSS天线和第二GNSS天线安装于所述仪器支架1的外部的上方；
[0073] 所述电池舱与所述LADCP6、组合惯导系统主机、数据处理单元、无线传输模块、第一GNSS天线和第二GNSS天线连接，进行供电，保证剖面流测量过程的各模块正常工作；
[0074] 所述第一GNSS天线和第二GNSS天线接收GNSS卫星信号并通过同轴线缆传输至所述组合惯导系统主机；
[0075] 所述LADCP实时观测海流剖面测量数据；
[0076] 所述组合惯导系统主机观测GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率；
[0077] 所述数据处理单元通过串口与无线传输模块连接，所述无线传输模块通过所述耐压无线天线与传输数据处理单元与母船甲板单元之间的通信数据；
[0078] 所述数据处理单元通过串口与所述LADCP和组合惯导系统主机连接，得到海流剖面测量数据、GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率，并通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正。
[0079] 在一些实施例中，所述联合观测系统还包括：第一耐压仓2、第二耐压仓3和第三耐压仓5；
[0080] 所述第一耐压仓2和第二耐压仓3位于所述仪器支架1的外部的上方；
[0081] 所述第一GNSS天线安装于所述第一耐压仓2内部；
[0082] 所述第二GNSS天线安装于所述第二耐压仓3内部；
[0083] 所述第三耐压仓5位于所述仪器支架1内部的所述LADCP的上方，所述组合惯导系统主机、数据处理单元、无线传输模块和电池舱安装于所述第三耐压仓5内部，所述耐压无线天线4安装于所述第三耐压仓5上方；
[0084] 所述第一耐压仓2和第二耐压仓3分别通过同轴水密电缆和水密接插件与所述第三耐压仓5连接，所述LADCP6通过水密电缆与所述第三耐压仓5连接。
[0085] 在一些实施例中，数据处理单元拟采用PC104架构的嵌入式采集器为核心，通过第一串口读取LADCP海流数据，通过第二串口读取组合惯导系统的GNSS位置及定向数据，通过第三串口读取组合惯导系统的三轴加速度数据，通过第四串口读取组合惯导系统的三轴角速度、横滚角、俯仰角和航向角速率，通过第五串口与无线传输模块进行数据通信，数据处理模块同步采集流速、流向、GNSS位置、GNSS方向、横滚角、俯仰角、三轴加速度、三轴角速度和对地速度等数据；仪器支架1采用无磁材料制作，多为铝或者316L不锈钢。所述LADCP的北方向与所述组合惯导系统主机的北方向对齐。
[0086] 上述系统的具体工作原理和操作方法如下：
[0087] 所述LADCP实时观测海流剖面的测量数据和所述组合惯导系统主机观测GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率的操作方法包括，
[0088] S11：观测船到达观测点位后，将所述LADCP的坐标系与所述组合惯导系统主机的坐标系对齐；
[0089] S12：将所述联合观测系统吊至海面，通过母船甲板单元接收所述数据处理单元的回传数据；
[0090] S13：判断所述组合惯导系统主机的GNSS位置和GNSS定向数据确定；
[0091] S14：将所述联合观测系统下放至海下预设深度，观测海流剖面测量数据；
[0092] S15：观测海流剖面测量数据完成后，将所述联合观测系统收回。
[0093] 如图3所示：在一些实施例中，所述通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正的方法包括，[0094] 由于设备进入水中后，GNSS位置和方向数据将失效，船体是一个巨大的铁磁性物体，会对LADCP的自带的磁罗经测得方向产生很大的影响，甚至会出现错误的方向数据；
[0095] S21：以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率进行离散积分，得到航向角；
[0096] S22：应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，即姿态校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度；
[0097] S23：应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，即线速度校正，得到第二校正三维流速；
[0098] S24：应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0099] 在一些实施例中，所述以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率离散积分的方法，得到航向角之前，所述方法还包括：
[0100] 消除LADCP的噪声和外界环境引起的LADCP的航向角速率离散采样值的瞬时粗大误差，具体方法包括，如图4所示：
[0101] S31：计算航向角速率离散采样值序列的算数平均值；
[0102] S32：计算航向角速率离散采样值序列的标准差；
[0103] S33：应用所述算数平均值和标准差，计算航向角速率离散采样值序列的标志值，具体公式为：
[0104]
[0105] 其中， 表示航向角速率离散采样值序列中采样值t的标志值，s表示标准差，表示序列中第t个航向角速率离散采样值， 表示算数平均值；
[0106] 计算航向角速率离散采样值序列的算数平均值：
[0107]
[0108] 其中，n表示航向角速率离散采样值序列中采样值的数量；
[0109] 计算航向角速率离散采样值序列的标准差：
[0110]
[0111] S34：应用序列中航向角速率离散采样值对应的标志值与预设的临界值比较，如果小于所述临界值，保留所述航向角速率离散采样值，否则剔除所述航向角速率离散采样值；
[0112] 在本实施例中，预设的临界值为3.28；
[0113] 滤除航向角速率离散采样值的直流分量，具体方法包括：
[0114] 对消除了瞬时粗大误差的航向角速率离散采样值进行离散傅里叶变换；
[0115] 将频谱图中频率为0的频谱分量置0后进行离散傅里叶逆变换，得到滤除直流分量的航向角速率离散采样值。
[0116] 在一些实施例中，所述以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率离散积分的方法，得到航向角的方法包括：
[0117] 所述以入水前最后一组没有失效的GNSS定向数据 作为初值，以最后一组没有失效的GNSS定向数据的采样时刻作为时间零点，通过对所述航向角速率离散积分的方法，求出预设时刻为止LADCP旋转的角度，对所述LADCP旋转的角度360取余，得到所述航向角，具体公式为：
[0118] 。
[0119] 在一些实施例中，所述应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，即姿态校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度的方法包括：
[0120]
[0121]
[0122]
[0123]
[0124]
[0125] 其中， 代表未校正的矢量测量结果， ，， ， 表示三维流
速， 表示三维线速度， 表示三维角速度； 代
表横滚角， 代表俯仰角， 代表航向角； 代表第一次校正后的矢量校
正结果， ， ，
， 表示第一校正三维流速，
表示第一校正三维线速度， 表示第一
校正三维角速度。
[0126] 在一些实施例中，所述应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，即线速度校正，得到第二校正三维流速的方法包括：
[0127] 应用所述第一校正三维线速度与所述第一校正三维流速进行矢量相加，得到第二校正三维流速，具体公式为：
[0128]
[0129] 其中， 表示第二校正三维流速；
[0130] 应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速的方法包括：
[0131] 应用LADCP的波束出口至仪器支架的吊点的距离R与所述第一校正三维角速度相乘，得到三维旋转线速度，再应用所述三维旋转线速度与所述第二校正三维流速进行矢量相加，得到第三校正三维流速，具体公式为：
[0132]
[0133] 其中， 表示第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0134] 实施例2：
[0135] 本发明第二方面公开了一种LADCP与组合惯性导航系统联合观测方法，采用第一方面所公开的一种LADCP与组合惯性导航系统联合观测方系统，如图2所示：所述LADCP实时观测海流剖面的测量数据和所述组合惯导系统主机观测GNSS位置、GNSS定向数据、三轴加速度数据、三轴陀螺仪角速度、横滚角、俯仰角和航向角速率的操作方法包括，[0136] 观测船到达观测点位后，将所述LADCP的坐标系与所述组合惯导系统主机的坐标系对齐；
[0137] 将所述联合观测系统吊至海面，通过母船甲板单元接收所述数据处理单元的回传数据；
[0138] 判断所述组合惯导系统主机的GNSS位置和GNSS定向数据确定；
[0139] 将所述联合观测系统下放至海下预设深度，观测海流剖面测量数据；
[0140] 观测海流剖面测量数据完成后，将所述联合观测系统收回。
[0141] 如图3所示：在一些实施例中，所述通过所述GNSS定向数据、横滚角、俯仰角和航向角速率对于LADCP观测的海流剖面的三维流速进行实时高精度校正的方法包括，[0142] 由于设备进入水中后，GNSS位置和方向数据将失效，船体是一个巨大的铁磁性物体，会对LADCP的自带的磁罗经测得方向产生很大的影响，甚至会出现错误的方向数据；
[0143] 以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率进行离散积分，得到航向角；
[0144] 应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，即姿态校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度；
[0145] 应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，即线速度校正，得到第二校正三维流速；
[0146] 应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0147] 在一些实施例中，所述以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率离散积分的方法，得到航向角之前，所述方法还包括：
[0148] 消除LADCP的噪声和外界环境引起的LADCP的航向角速率离散采样值的瞬时粗大误差，具体方法包括，如图4所示：
[0149] 计算航向角速率离散采样值序列的算数平均值；
[0150] 计算航向角速率离散采样值序列的标准差；
[0151] 应用所述算数平均值和标准差，计算航向角速率离散采样值序列的标志值，具体公式为：
[0152]
[0153] 其中， 表示航向角速率离散采样值序列中采样值t的标志值，s表示标准差，表示序列中第t个航向角速率离散采样值， 表示算数平均值；
[0154] 计算航向角速率离散采样值序列的算数平均值：
[0155]
[0156] 其中，n表示航向角速率离散采样值序列中采样值的数量；
[0157] 计算航向角速率离散采样值序列的标准差：
[0158]
[0159] 应用序列中航向角速率离散采样值对应的标志值与预设的临界值比较，如果小于所述临界值，保留所述航向角速率离散采样值，否则剔除所述航向角速率离散采样值；
[0160] 在本实施例中，预设的临界值为3.28；
[0161] 滤除航向角速率离散采样值的直流分量，具体方法包括：
[0162] 对消除了瞬时粗大误差的航向角速率离散采样值进行离散傅里叶变换；
[0163] 将频谱图中频率为0的频谱分量置0后进行离散傅里叶逆变换，得到滤除直流分量的航向角速率离散采样值。
[0164] 在一些实施例中，所述以入水前最后一组没有失效的GNSS定向数据 作为初值，通过对所述航向角速率离散积分的方法，得到航向角的方法包括：
[0165] 所述以入水前最后一组没有失效的GNSS定向数据 作为初值，以最后一组没有失效的GNSS定向数据的采样时刻作为时间零点，通过对所述航向角速率离散积分的方法，求出预设时刻为止LADCP旋转的角度，对所述LADCP旋转的角度360取余，得到所述航向角，具体公式为： 。
[0166] 在一些实施例中，所述应用所述横滚角、俯仰角和航向角对海流剖面的三维流速、三维线速度和三维角速度进行第一次校正，即姿态校正，得到第一校正三维流速、第一校正三维线速度和第一校正三维角速度的方法包括：
[0167]
[0168]
[0169]
[0170]
[0171]
[0172] 其中， 代表未校正的矢量测量结果， ，， ， 表示三维流
速， 表示三维线速度， 表示三维角速度； 代
表横滚角， 代表俯仰角， 代表航向角； 代表第一次校正后的矢量校
正结果， ， ，
， 表示第一校正三维流速，
表示第一校正三维线速度， 表示第一
校正三维角速度。
[0173] 在一些实施例中，所述应用所述第一校正三维线速度对所述第一校正三维流速进行第二次校正，即线速度校正，得到第二校正三维流速的方法包括：
[0174] 应用所述第一校正三维线速度与所述第一校正三维流速进行矢量相加，得到第二校正三维流速，具体公式为：
[0175]
[0176] 其中， 表示第二校正三维流速；
[0177] 应用所述第一校正三维角速度对所述第二校正三维流速进行第三次校正，即角速度校正，得到第三校正三维流速的方法包括：
[0178] 应用LADCP的波束出口至仪器支架的吊点的距离R与所述第一校正三维角速度相乘，得到三维旋转线速度，再应用所述三维旋转线速度与所述第二校正三维流速进行矢量相加，得到第三校正三维流速，具体公式为：
[0179]
[0180] 其中， 表示第三校正三维流速，即海流剖面的三维流速校正的最终结果。
[0181] 综上，本发明各个方面的技术方案与现有技术相比具有如下优点：提高LADCP流速和流向的测量精度，降低了使用LADCP进行剖面流观测的工作成本。
[0182] 请注意，以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。以上实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
