基于PTV与PIV技术的海表面小尺度流场测量方法技术

基于PTV与PIV技术的海表面小尺度流场测量方法，包括无人机拍摄校准；待测区域连续图像拍摄；PTV方法初步分析得到锋面流速；PIV方法实现高精度的海表面流场分析，得到流场各点处的物理参量；利用灰度极值法求周期，计算深水波波速；将PIV流速与波浪传播速度相减得到每一固定点处水团流速；将锋面流速、流场各点处的物理参量、水团流速进行输出，即实现了海表面小尺度流场的测量。本发明专利技术将PTV与PIV技术用于海表面小尺度流场的测量，利用无人机进行拍摄是更高效简便的手段，在小于几公里的小尺度海面实现海表面流场的可视化，提供重要的观测数据，有利于后续对海表面运动过程与海水内部其他复杂化运动过程进行研究。部其他复杂化运动过程进行研究。部其他复杂化运动过程进行研究。

【技术实现步骤摘要】
基于PTV与PIV技术的海表面小尺度流场测量方法


[0001]本专利技术涉及一种基于PTV与PIV技术的海表面小尺度流场测量方法，属于流场测量


技术介绍

[0002]海洋表面流场测量对于获取复杂水流运动信息至关重要。近年来，基于声学、电波、卫星遥感、粒子图像等非接触式测流仪器促进了流量自动化测量技术的发展，声学测量方式目前在海洋测量上应用较多，主要使用的是多普勒流速仪（ADCP）设备，ADCP分为走航式、底座式和水平式三种，由于测量盲区的存在，该仪器主要用于海表面以下的海水流速测量，它离海表层深约一到两米内流速是无法测量的，所以在分析的时候一般都是把水深两到三米处流速当做表层流速。但是通 过研究漂浮物漂移的速度发现，海表层的流速与海水内部的流速往往有很大区别，非常靠近表层的速度和水深以下两米的速度相差较大，即说明运用ADCP的方法测量表层流速误差较大，因此亟需寻找测量海表面流场的方法。
[0003]以往研究也有采用地波雷达和表面追踪浮标等进行海表面运动的观测。然而，这两种方法测量的表面流场范围在几十至几百公里之间，难以满足小尺度范围内流场的连续观测。此外，利用雷达多普勒效应实现流速测量时，其精度主要受测频精度和波浪两个因素影响，当海洋运动现象幅度大时，监测中存在测验适宜性弱、精度稳定性低、施测安全性差及测量结果不直观等问题，且雷达流速仪所需成本较高。

技术实现思路

[0004]本专利技术目的在于：为了获取准确的海表面小尺度流场，实现海表面流场的可视化，并降低成本，提供一种基于PTV与PIV技术的海表面小尺度流场测量方法。
[0005]为实现以上技术目的，本专利技术采用如下技术方案：（1）无人机拍摄校准：调节无人机拍摄参数，计算不同高度下无人机拍摄画面的实际尺寸，利用拍摄画面中已知尺寸的物体进行比较，得到校准系数，进而得到无人机拍摄画面的实际长宽与高度的对应关系。
[0006]其特征是还包括以下步骤：（2）待测区域连续图像拍摄：将无人机飞行至待测海面区域，调节无人机拍摄高度，固定无人机位置拍摄待测海面区域的水色图像。
[0007]（3）PTV（Particle Track Velocimetry, 粒子追踪测速法）方法初步分析：根据图像的分辨率和帧率，运用校准系数计算实际拍摄图像尺寸大小，利用粒子跟踪测速技术PTV(Particle Track Velocimetry)分析待测海面区域锋面移动速度，具体做法是根据所拍摄图像的灰度值大小确认锋面位置，根据时间帧数确认实际时间，计算一段时间内锋面移动的实际距离与对应实际时间的比值，即海表面处锋面流速。
[0008]（4）PIV（Particle Image Velocimetry，粒子图像测速法）方法实现高精度的海表面流场分析：首先将图像中呈现的所有像素分割成m
×
n个窗口，根据粒子在水中不同的反光性，追踪窗口中的粒子，将每个窗口作为一个整体，利用粒子随时间的运动轨迹判断其运动状态，采用互相关算法分析出流场特征；针对表示海表面微小粒子的各个像素，运用PIVLAB软件，计算流场各点处的物理参量；所述物理参量包括流速、涡度、u方向速度，v方向速度，剪切速率，应变速率，切线速度等，可绘制流速分布图及速度曲线，实现海表面流场可视化。
[0009]（5）灰度极值法实现波浪速度的计算：深水波波速与水深无关，只与波长有关，利用下列公式，，得到c与T的对应关系，其中，c深水波波速，λ深水波波长，g重力加速度，T深水波周期；因此若已知周期可实现波浪传播速度的计算，然后据此计算深水波波速。本专利技术中提出一种新的方法求周期——灰度极值法。
[0010]本专利技术的灰度极值法是一种根据图像固定像素点在一定时间内的灰度值变化求极值点个数，将该极值点个数等效为波峰数量，进而求波速的方法。具体步骤如下：首先利用MATLAB计算出一段时间内拍摄的各帧图像中各个像素点的灰度值，输出数据结果及灰度图；对于图像中的任意一个固定点，将该段时间内该固定点的灰度值变化输出，利用diff语句计算固定点灰度值变化曲线的极值点个数，该步骤中对极值点个数计算时采用低通滤波进行优化，排除小概率灰度值的干扰；极值点个数近似认为该固定点的波峰个数，根据时间得到波浪传播周期，计算出该固定点处的波浪传播速度，进而得到图像中每一个固定点的波浪传播速度即深水波波速。
[0011]（6）滤波后海表面水团流速的计算：所述步骤（4）中利用PIV方法计算出的流场流速即海表面处水团和海表面处波浪的合速度，要想分析海表面水团流速需要进一步进行滤波，排除波浪因素的干扰。
[0012]通过所述步骤（4）得到的结果可知各固定点利用PIV方法得到的流速，又通过所述步骤（5）获得各固定点利用灰度极值法计算得到的波浪传播速度，将PIV流速与波浪传播速度相减得到每一固定点处水团流速，绘制整个图像的水团速度曲线。
[0013]（7）将步骤（3）得到的锋面流速、步骤（4）得到的流场各点处的物理参量、步骤（6）得到的水团流速进行输出，即实现了海表面小尺度流场的测量。
[0014]进一步的，在校准阶段可以多次测量，取不同尺寸的画面进行校准，以减小误差。
[0015]进一步的，进行一次校准后短时期内可不用再进行校准，同时，若要更改无人机拍摄模式需再次进行校准。
[0016]所述步骤（5）中，利用公式，c=得到c与T的对应关系为c=1.56T。
[0017]由于近海区海表面悬浮物浓度较高，因此海表面流场可以由水色图像清晰识别，利用无人机拍摄连续获得海表面水色图像，利用MATLAB软件以及PTV和PIV技术可计算出流场中流速分布、海表面流场涡度、散度等参量，实现海表面小尺度流场可视化。
[0018]本专利技术的有益效果在于：本专利技术基于PTV与PIV技术，将其利用于海表面小尺度流场的测量，利用无人机进行拍摄，相比于当前对海洋中海水运动现象研究主要有海面流场
卫星遥感方法和基于紧凑型高频雷达的海表面流场探测等方法，卫星遥感进行流场测量方法的局限性在于不能进行小尺度海表面流场研究，而雷达的方法则面临测量费用过高以及不能进行小于几公里的海表面流场运动分析的问题。可见，利用无人机观测测量是更高效简便的手段，在小于几公里的小尺度海面实现海表面流场的可视化，提供重要的观测数据，有利于后续对海表面运动过程与海水内部其他复杂化运动过程进行研究。
附图说明
[0019]图1是基于PTV与PIV技术的海表面小尺度流场的测量方法的流程示意图。
[0020]图2是无人机拍摄图像的灰度图。
[0021]图3是PIV方法分析的流场图。
[0022]图4是PIV方法分析的流速分布图。
[0023]图5是图4的局部放大图。
具体实施方式
[0024]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面对本专利技术实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。
[0025]本本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于PTV与PIV技术的海表面小尺度流场测量方法，包括以下步骤：（1）无人机拍摄校准：调节无人机拍摄参数，计算不同高度下无人机拍摄画面的实际尺寸，利用拍摄画面中已知尺寸的物体进行比较，得到校准系数，进而得到无人机拍摄画面的实际长宽与高度的对应关系；其特征是还包括以下步骤：（2）待测区域连续图像拍摄：将无人机飞行至待测海面区域，调节无人机拍摄高度，固定无人机位置拍摄待测海面区域的水色图像；（3）PTV方法初步分析：根据图像的分辨率和帧率，运用校准系数计算实际拍摄图像尺寸大小，利用PTV方法分析待测海面区域锋面移动速度；（4）PIV方法实现高精度的海表面流场分析：首先将图像中呈现的所有像素分割成m
×
n个窗口，根据粒子在水中不同的反光性，追踪窗口中的粒子，将每个窗口作为一个整体，利用粒子随时间的运动轨迹判断其运动状态，采用互相关算法分析出流场特征；针对表示海表面微小粒子的各个像素，运用PIVLAB软件，计算流场各点处的物理参量；（5）灰度极值法实现波浪速度的计算：利用灰度极值法求周期，利用公式，，得到c与T的对应关系，计算深水波波速；其中，c深水波波速，λ深水波波长，g重力加速度，T深水波周期；（6）滤波后海表面水团流速的计算：通过所述步骤（4）得到的结果可知各固定点利用PIV方法得到的流速，又通过所述步骤（5）获得各固定点利用灰度极值法计算得到的波浪传播速度，将PIV流速与波浪传播速度相减得到每一固定点处水团流速，...

【专利技术属性】
技术研发人员：王涛，孙超凡，菅秀阳，王竞怡，许胆，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：