识别气液界面位置的方法转让专利
申请号 : CN201910612337.0
文献号 : CN110231068B
文献日 : 2020-12-08
发明人 : 李存标 , 陈钧伟
申请人 : 北京大学
摘要 :
本发明公开了一种识别气液界面位置的方法。本发明的识别气液界面位置的方法，首先分别在气体和液体中分撒示踪粒子，然后通过在气液界面的两侧各自使用高速相机拍摄示踪粒子图像，然后识别出图像中的初始时刻的界面位置，然后再对所有同时刻的图像进行预处理，识别出随后各帧图像中的气液界面位置，再对得到的气液界面位置进行时间‑空间平滑，从而得到了拍摄界面以上区域的相机所得到的图像中的所有气液界面位置，最后再得到拍摄拍摄界面以下区域的相机所得到的图像中的气液界面位置，可以在气液界面两侧均有粒子存在的情况下准确识别气液界面的位置，且识别准确度高，可以很好地适用于粒子图像测速中的近壁测量。
权利要求 :
1.一种识别气液界面位置的方法，适用于气液界面两侧均有示踪粒子存在的情况，其特征在于，包括以下步骤：
步骤S1：分别在气体和液体中播撒示踪粒子，然后设置片光源进行照明；
步骤S2：在气液界面的两侧各使用高速相机拍摄示踪粒子图像，并对位于界面两侧的相机各自进行标定；
步骤S3：识别拍摄界面以上区域的相机所得到的图像中的初始时刻界面位置；
步骤S4：对拍摄界面以上区域的相机所得到的同时刻的图像进行预处理；
步骤S5：识别随后各帧图像中的气液界面位置；
步骤S6：对气液界面位置进行时间-空间平滑操作以得到拍摄界面以上区域的相机所得到的图像中的气液界面位置；
步骤S7：根据拍摄界面以上区域的相机所得到的图像中的气液界面位置得到拍摄界面以下区域的相机所得到的图像中的气液界面位置；
所述步骤S3具体包括以下步骤：
步骤S31：若气体运动速度远大于液体运动速度，每个像素的亮度取多张图片中相应位置亮度的最小者；否则跳过本步骤；
步骤S32：对图片进行灰度形态学开操作若干次；
步骤S33：设置亮度阈值识别液体中的粒子，得到二值图像；
步骤S34：对二值图像进行开操作若干次，再进行闭操作若干次；
步骤S35：根据图像中的界面高度设置高度数值y0，使界面的高度位置在高度为该数值的直线以下，高度在y0以上的部分设置为0；
步骤S36：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置，得到一个序列；
步骤S37：对序列进行中值滤波和高斯滤波处理，得到界面的位置；
步骤S39：设置偏置量，然后根据偏置量向上移动初始界面位置；
所述步骤S5具体包括以下步骤：
步骤S51：对于预处理得到的图像，设置亮度阈值识别液体中的粒子以得到二值图像；
步骤S52：对二值图像进行开操作若干次，再进行闭操作若干次；
步骤S53：将前一帧图像的界面向上平移15像素作为分界线，分界线以上区域二值图像数值设置为0；
步骤S54：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置和前一帧界面相同位置的较小值，从而得到一个序列；
步骤S55：对序列进行中值滤波和高斯滤波得到界面的位置；
步骤S56：根据偏置量向上移动界面的位置；
所述步骤S6具体包括以下步骤：
步骤S61：相隔数帧计算每两帧之间界面位置的互相关值；
步骤S62：将多个互相关值相加后寻找峰值位置；
步骤S63：针对所有可计算的时刻重复步骤S61和步骤S62，得到界面速度的时间序列；
步骤S64：基于得到的界面速度的时间序列对界面位置做锁相平均计算；
步骤S65：对所有时刻的边界位置做时间-空间平滑操作，得到平滑后的边界位置时间-空间序列；
当相机采用等时间间距曝光方式时，所述步骤S64具体为：
记t0时刻的界面位置为表示界面高度，中对应的分量；记步骤S61中计算互相关值的时间间隔为△t，经过平滑后得到t0时刻界面在△t内移动距离为<△x(t0)>；然后记N为时间-空间平滑半径，计算
当相机采用双帧曝光方式时，所述步骤S64具体为：
双帧之间时间间隔为Δt，双帧曝光周期为ΔT，则t0时刻界面在△t内移动距离为<△x(t0)>，在ΔT内移动距离为<△X(t0)>＝<△x(t0)>*ΔT/Δt，对于双帧曝光的第一帧，计算
2.如权利要求1所述的识别气液界面位置的方法，其特征在于，所述步骤S2中相机与水平面之间呈10°夹角倾斜设置。
3.如权利要求1所述的识别气液界面位置的方法，其特征在于，所述步骤S7具体包括以下步骤：
步骤S71：根据标定得到的坐标变换关系将拍摄界面以下区域的相机所得到的图像从原始图像变换到物理空间坐标系下的图像；
步骤S73：对边界位置时间-空间序列进行缩放并插值以适应经过坐标变换的拍摄界面以下区域的相机所得到的图像；
步骤S74：暂时给定界面位置的初始偏移量，在一张图片上显示相机经过空间变换的图像和界面经过偏移后暂时的位置，输出由不同时刻的图片组成的视频；
步骤S75：观察视频中暂定的界面位置和由图片中粒子亮度反应的界面位置的差异，并根据差异修改界面位置的偏移量；
步骤S76：重复执行步骤S74和步骤S75，直至界面暂时的位置和图片中粒子亮度反映的界面位置重合。
4.如权利要求1所述的识别气液界面位置的方法，其特征在于，当所述步骤S2中在气液界面的两侧各使用两台高速相机拍摄示踪粒子图像时，所述识别气液界面位置的方法在步骤S2和步骤S3之间还包括以下步骤：步骤S23：把经过坐标变换的两台相机同时刻的照片的亮度转为红/绿通道的数值，然后合并到一张图片上显示，并判断两张照片的界面是否重合。
5.如权利要求4所述的识别气液界面位置的方法，其特征在于，所述步骤S4具体包括以下步骤：
步骤S41：若气体运动速度远大于液体运动速度，每个像素的亮度取多张图片中相应位置亮度的最小者；否则跳过本步骤；
步骤S42：将两台相机所得到的同时刻的照片合成一张照片；
步骤S43：对合成的照片进行灰度形态学开操作若干次。
说明书 :
识别气液界面位置的方法
技术领域
[0001] 本发明涉及粒子图像测速技术领域，特别地，涉及一种识别气液界面位置的方法。
背景技术
[0002] 粒子图像测速(PIV)是流体力学实验中经常使用的一种测量方法，它可以测量某个区域内流体瞬时的速度场，其原理是在流体中播撒随流体微团一起运动的示踪粒子，在短时间内拍摄两张照片，由两个时刻内示踪粒子在照片上位置的变化可以计算出拍摄区域内各处流体微团的速度(这个速度为平面上的速度，即有两个分量)。其中最主要的计算方法是将画面分成若干小块(查询窗口)，计算较前时刻的小块和较后时刻的小块之间照片亮度的相关系数，相关系数峰值的位置表示流体相应查询窗口附近流体微团的速度。粒子图像测速有干扰小、可测量某区域同时刻所有位置速度等优点。
[0003] 自粒子图像测速大规模使用以来，近壁测量一直是其中的比较疑难的部分。壁面处的强烈反光，壁面附近速度梯度大、粒子密度低，壁面另一侧缺少有效信息，以及壁面形状不规则或形状随时间变化等因素都会增加测量难度，而界面位置又是进行近壁测量至关重要的信息。现有的粒子图像测速方法在进行近壁测量时只能在测量界面一侧的速度场的过程中识别气液界面的位置，无法适用于界面两侧都有粒子存在的情况，且其界面位置识别不准确，测量结果准确度较差。
发明内容
[0004] 本发明提供了一种识别气液界面位置的方法，以解决现有的粒子图像测速方法在进行近壁测量时只能在界面一側有示踪粒子时识别气液界面的位置，无法适用于界面两侧都有粒子存在的情况，且其界面位置识别不准确，测量结果准确度较差的技术问题。
[0005] 根据本发明的一个方面，提供一种识别气液界面位置的方法，适用于气液界面两侧均有示踪粒子存在的情况，包括以下步骤：
[0006] 步骤S1：分别在气体和液体中播撒示踪粒子，然后设置片光源进行照明；
[0007] 步骤S2：在气液界面的两侧各使用高速相机拍摄示踪粒子图像，并对位于界面两侧的相机各自进行标定；
[0008] 步骤S3：识别拍摄界面以上区域的相机所得到的图像中的初始时刻界面位置；
[0009] 步骤S4：对拍摄界面以上区域的相机所得到的同时刻的图像进行预处理；
[0010] 步骤S5：识别随后各帧图像中的气液界面位置；
[0011] 步骤S6：对气液界面位置进行时间-空间平滑操作以得到拍摄界面以上区域的相机所得到的图像中的气液界面位置；
[0012] 步骤S7：根据拍摄界面以上区域的相机所得到的图像中的气液界面位置得到拍摄界面以下区域的相机所得到的图像中的气液界面位置。
[0013] 进一步地，所述步骤S2中相机与水平面之间呈10°夹角倾斜设置。
[0014] 进一步地，所述步骤S3具体包括以下步骤：
[0015] 步骤S31：若气体运动速度远大于液体运动速度，每个像素的亮度取多张图片中相应位置亮度的最小者；否则跳过本步骤；
[0016] 步骤S32：对图片进行灰度形态学开操作若干次；
[0017] 步骤S33：设置亮度阈值识别液体中的粒子，得到二值图像；
[0018] 步骤S34：对二值图像进行开操作若干次，再进行闭操作若干次；
[0019] 步骤S35：根据图像中的界面高度设置高度数值y0，使界面的高度位置在高度为该数值的直线以下，高度在y0以上的部分设置为0；
[0020] 步骤S36：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置，得到一个序列；
[0021] 步骤S37：对序列进行中值滤波和高斯滤波处理，得到界面的位置。
[0022] 进一步地，所述步骤S3在步骤S37之后还包括以下步骤：
[0023] 步骤S39：设置偏置量，然后根据偏置量向上移动初始界面位置。
[0024] 进一步地，所述步骤S5具体包括以下步骤：
[0025] 步骤S51：对于预处理得到的图像，设置亮度阈值识别液体中的粒子以得到二值图像；
[0026] 步骤S52：对二值图像进行开操作若干次，再进行闭操作若干次；
[0027] 步骤S53：将前一帧图像的界面向上平移15像素作为分界线，分界线以上区域二值图像数值设置为0；
[0028] 步骤S54：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置和前一帧界面相同位置的较小值，从而得到一个序列；
[0029] 步骤S55：对序列进行中值滤波和高斯滤波得到界面的位置；
[0030] 步骤S56：根据偏置量向上移动界面的位置。
[0031] 进一步地，所述步骤S6具体包括以下步骤：
[0032] 步骤S61：相隔数帧计算每两帧之间界面位置的互相关值；
[0033] 步骤S62：将多个互相关值相加后寻找峰值位置；
[0034] 步骤S63：针对所有可计算的时刻重复步骤S61和步骤S62，得到界面速度的时间序列；
[0035] 步骤S64：基于得到的界面速度的时间序列对界面位置做锁相平均计算；
[0036] 步骤S65：对所有时刻的边界位置做时间-空间平滑操作，得到平滑后的边界位置时间-空间序列。
[0037] 进一步地，所述步骤S7具体包括以下步骤：
[0038] 步骤S71：根据标定得到的坐标变换关系将拍摄界面以下区域的相机所得到的图像从原始图像变换到物理空间坐标系下的图像；
[0039] 步骤S73：对边界位置时间-空间序列进行缩放并插值以适应经过坐标变换的拍摄界面以下区域的相机所得到的图像；
[0040] 步骤S74：暂时给定界面位置的初始偏移量，在一张图片上显示相机经过空间变换的图像和界面经过偏移后暂时的位置，输出由不同时刻的图片组成的视频；
[0041] 步骤S75：观察视频中暂定的界面位置和由图片中粒子亮度反应的界面位置的差异，并根据差异修改界面位置的偏移量；
[0042] 步骤S76：重复执行步骤S74和步骤S75，直至界面暂时的位置和图片中粒子亮度反映的界面位置重合。
[0043] 进一步地，当相机采用等时间间距曝光方式时，所述步骤S64具体为：
[0044] 记t0时刻的界面位置为，t0)>，其中表示界面高度，表示界面的横坐标序列，<>表示矢量，界面位置，t0)>即表示在t0时刻中某点x0所在列的界面高度在中对应的分量；记步骤S61中计算互相关值的时间间隔为△t，经过平滑后得到t0时刻界面在△t内移动距离为<△x(t0)>；然后记N为时间-空间平滑半径，计算-N<Δx(t0)>，t0-NΔt)>，-(N-1)<Δx(t0)>，t0–(N-1)Δt)>，…，,t0)>，+(N-1)<Δx(t0)>，t0+(N-1)Δt)>，…+N<Δx(t0)>，t0+NΔt)>在上的平均值，t0)>，其中表示将向两侧延伸后得到的新横坐标序列；
[0045] 当相机采用双帧曝光方式时，所述步骤S64具体为：
[0046] 双帧之间时间间隔为Δt，双帧曝光周期为ΔT，则t0时刻界面在△t内移动距离为<△x(t0)>，在ΔT内移动距离为<△X(t0)>＝<△x(t0)>*ΔT/Δt，对于双帧曝光的第一帧，计算-N<ΔX(t0)>，t0-NΔT)>，-(N-1)<ΔX(t0)>，t0–(N-1)ΔT)>，…，，t0)>，+(N-1)<ΔX(t0)>，t0+(N-1)ΔT)>，…，+N<ΔX(t0)>，t0+NΔT)>，-N<ΔX(t0)>+<△x(t0)>，t0–NΔT+Δt)>，-(N-1)<ΔX(t0)>+<△x(t0)>，t0–(N-1)ΔT+Δt)>，…，+<△x(t0)>，t0+Δt)>，+(N-1)<ΔX(t0)>+<△x(t0)>，t0+(N-1)ΔT+Δt)>，…，+N<ΔX(t0)>+<△x(t0)>，t0+NΔT+Δt)>的平均值，对于双帧曝光的第二帧，计算-N<ΔX(t0)>，t0-NΔT)>，-(N-1)<ΔX(t0)>，t0–(N-1)ΔT)>，…，，t0)>，+(N-1)<ΔX(t0)>，t0+(N-1)ΔT)>，…，+N<ΔX(t0)>，t0+NΔT)>，-N<ΔX(t0)>+<△x(t0)>，t0–NΔT-Δt)>，-(N-1)<ΔX(t0)>+<△x(t0)>，t0–(N-1)ΔT-Δt)>，…，+<△x(t0)>，t0-Δt)>，+(N-1)<ΔX(t0)>+<△x(t0)>，t0+(N-1)ΔT-Δt)>，…，+N<ΔX(t0)>+<△x(t0)>，t0+NΔT-Δt)>的平均值。
[0047] 进一步地，当所述步骤S2中在气液界面的两侧各使用两台高速相机拍摄示踪粒子图像时，所述识别气液界面位置的方法在步骤S2和步骤S3之间还包括以下步骤：
[0048] 步骤S23：把经过坐标变换的两台相机同时刻的照片的亮度转为红/绿通道的数值，然后合并到一张图片上显示，并判断两张照片的界面是否重合。
[0049] 进一步地，所述步骤S4具体包括以下步骤：
[0050] 步骤S41：若气体运动速度远大于液体运动速度，每个像素的亮度取多张图片中相应位置亮度的最小者；否则跳过本步骤；
[0051] 步骤S42：将两台相机所得到的同时刻的照片合成一张照片；
[0052] 步骤S43：对合成的照片进行灰度形态学开操作若干次。
[0053] 本发明具有以下有益效果：
[0054] 本发明的识别气液界面位置的方法，首先分别在气体和液体中分撒示踪粒子，然后通过在气液界面的两侧各自使用高速相机拍摄示踪粒子图像，然后识别出拍摄界面以上区域的相机所得到的图像中的初始时刻的界面位置，然后再对所有同时刻的图像进行预处理，识别出随后各帧图像中的气液界面位置，再对得到的气液界面位置进行时间-空间平滑，从而得到了拍摄界面以上区域的相机所得到的图像中的所有气液界面位置，最后再根据拍摄界面以上区域的相机所得到的图像中的气液界面位置得到拍摄拍摄界面以下区域的相机所得到的图像中的气液界面位置。本发明的识别气液界面位置的方法，可以在气液界面两侧均有粒子存在的情况下准确识别气液界面的位置，且识别准确度高，可以很好地适用于粒子图像测速中的近壁测量过程。
[0055] 除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照图，对本发明作进一步详细的说明。
附图说明
[0056] 构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
[0057] 图1是本发明优选实施例的识别气液界面位置的方法的流程示意图。
[0058] 图2是本发明优选实施例的图1中步骤S1中设置两个激光器来进行照明的实验装置示意图。
[0059] 图3是本发明优选实施例的图1中步骤S1中设置一个激光器来进行照明的实验装置示意图。
[0060] 图4是本发明优选实施例的图1中的步骤S2中对界面两侧的相机各自进行标定的示意图。
[0061] 图5是本发明优选实施例的识别气液界面位置的方法在气液界面分别设置两台相机进行拍摄时的步骤S3的子流程示意图。
[0062] 图6是本发明优选实施例的图1中的步骤S4的子流程示意图。
[0063] 图7是本发明优选实施例的图1中的步骤S5的子流程示意图。
[0064] 图8是本发明优选实施例的图1中的步骤S6的子流程示意图。
[0065] 图9是本发明优选实施例的图1中的步骤S7的子流程示意图。
[0066] 图10是本发明优选实施例的识别气液界面位置的方法在气液界面分别设置两台相机进行拍摄时的步骤S7的子流程示意图。
具体实施方式
[0067] 以下结合附图对本发明的实施例进行详细说明，但是本发明可以由下述所限定和覆盖的多种不同方式实施。
[0068] 如图1所示，本发明的优选实施例提供一种识别气液界面位置的方法，可以在测量气液界面两侧速度场的过程中识别气液界面的位置，可以适应于界面两侧都有示踪粒子存在的情况，且界面位置识别的准确度很高。所述识别气液界面位置的方法包括以下步骤：
[0069] 步骤S1：分别在气体和液体中播撒示踪粒子，然后设置片光源进行照明；
[0070] 步骤S2：在气液界面的两侧各使用高速相机拍摄示踪粒子图像，并对位于界面两侧的相机各自进行标定；
[0071] 步骤S3：识别拍摄界面以上区域的相机所得到的图像中的初始时刻界面位置；
[0072] 步骤S4：对拍摄界面以上区域的相机所得到的同时刻的图像进行预处理；
[0073] 步骤S5：识别随后各帧图像中的气液界面位置；
[0074] 步骤S6：对气液界面位置进行时间-空间平滑操作以得到拍摄界面以上区域的相机所得到的图像中的气液界面位置；
[0075] 步骤S7：根据拍摄界面以上区域的相机所得到的图像中的气液界面位置得到拍摄界面以下区域的相机所得到的图像中的气液界面位置。
[0076] 可以理解，在所述步骤S1中，考虑到粒子跟随流体的性能与斯托克斯数有关，确保粒子跟随性足够好则需要直径足够小的粒子，但是另一方面，在相机记录的图像中，粒子的亮度与粒子直径的平方有关，这要求粒子半径又不能太小。因此，本实施例在在气体中使用粒径较小的粒子进行播撒，在液体中使用粒径较大的粒子进行播撒，在气液界面附近的液体中重复播撒粒子以使该区域粒子密度足够大，从而确保粒子具有足够跟随性的同时，两种粒径的粒子在图像中的各自的亮度和粒径有所区别。作为优选的，在气体中使用直径为1μm左右的烟雾作为示踪粒子，在液体中使用直径为50μm左右的空心有机玻璃微粒作为示踪粒子，其中，气体中的示踪粒子由风扇注入，液体中的粒子一部分由漏斗注入比较深的位置，另一部分轻轻倒在液面上。如图2所示，在测量区域设置片光源进行照明，具体的设置方式可以是分别在气液界面的两侧分别设置一个片光源进行照明，照明方向分别为垂直向上和垂直向下，两个光源的颜色可以不同，从而便于被滤光片区分。其中，片光源可以采用括束后的激光，为了保护两台激光器，在调整光路的过程中，可以微调激光器和光路上各元件的位置，以使两台激光器发出的光线互补直射进对方腔体中。另外，如图3所示，当激光器数量不足时，也可以只使用一台激光器照明，照明方向为从下向上穿过气液界面照射。
[0077] 可以理解，在所述步骤S2中，在气液界面的两侧各自使用高速相机拍摄示踪粒子图像，即设置高速相机分别拍摄界面以上区域的示踪粒子图像和界面以下区域的示踪粒子图像。在本优选实施例中，可以在气液界面的两侧分别设置一台相机进行拍摄，也可以在气液界面的两侧分别设置两台相机进行拍摄。另外，相机采用的曝光方式为等时间间距曝光或者双帧曝光，其中，等时间间距曝光是指在均匀的时间间隔内拍摄一系列图片，双帧曝光是指进行两次时间间隔较短的曝光，用来计算速度，经过一段较长的时间后再有接下来两次曝光。作为优选的，本优选实施例采用时间序列立体粒子图像测速方法同时测量界面两侧区域的速度场，具体是在界面两侧各使用两台高速相机拍摄示踪粒子图像，即两台相机拍摄气液界面以上的区域，两台相机拍摄气液界面以下的区域，然后进行等时间间距曝光，即在均匀的时间间隔内拍摄一系列图片。另外，如图4所示，所述步骤S2中进行标定的方式具体为：要求气液界面两侧标定的基准平面和激光照射的平面重合，标定板底边水平放置，坐标原点之间可以有距离，即两个进行标定的物理空间坐标系中，两者的y轴正方向均垂直向上，z轴均直于标定板前表面，z＝0平面重合，x，y，z方向相同。另外，作为优选的，在拍摄过程中为了避免测量所在的平面之前的界面遮挡测量所在的平面上的界面，需要将相机相对于水平面倾斜设置，倾斜角度为10°左右。可以理解，当相机与水平面之间呈10°左右的夹角时，拍摄气液界面以下区域的相机得到的图像中主要有两种内容：气液界面以下区域示踪粒子散射光直接在相机中成像、气液界面以下区域示踪粒子的散射光经过界面反射后在相机中成像，而拍摄气液界面以上区域的相机得到的图像中主要有两种内容：气体中的粒子直接成像、液体中的粒子经过界面折射后成像。考虑到，经过反射的散射光随着与气液界面距离增加，变形和脱焦现象愈加严重，边缘也更加模糊，并且米氏散射中光强和散射角度之间存在关系，散射光和照射光夹角比较小时散射光强度比较大，所以气液界面以上区域的粒子比界面以下区域的粒子亮度大，且距离界面越远粒子越模糊。因此，本发明采用拍摄气液界面以上区域的相机得到的图像识别界面位置，识别结果更加准确。
[0078] 可以理解，在所述步骤S3中，识别拍摄界面以上区域的相机得到的图像中界面的位置需要根据标定得到的坐标映射关系将图像变换到物理空间坐标，因此需要检验标定得到的坐标映射关系进行是否准确。当在气液界面的两侧各自使用两台高速相机拍摄示踪粒子图像时，作为优选的，所述识别气液界面位置的方法在步骤S2和步骤S3之间还包括步骤S23，
[0079] 步骤S23：把经过坐标变换的拍摄界面以上区域的两台相机同时刻的照片的亮度转为红/绿通道的数值，然后合并到一张图片上显示，并判断两张照片的界面是否重合。如果标定得到的坐标映射关系准确，这时两张图片中界面的位置将会重合；如果标定得到的坐标映射关系不准确，两台相机的图片经过坐标变化后界面位置不重合，则需要对两台相机进行自标定，即可以根据两张图片中界面相对点位置的距离设置初始偏移量，即相应点在x方向和y方向上的距离，再由这些偏移确定总体上的偏移，总体的偏移值为各点偏移值的平均值。为了确保准确性，可以使用两台相机的多对图片进行自标定，例如100对图片。
[0080] 可以理解，在所述步骤S2和步骤S3之间还需执行以下步骤：
[0081] 当拍摄界面以上区域的相机方向和拍摄界面以下区域的相机方向相反时，需要将拍摄界面以上区域的两台相机所得到的图像水平翻转，或者将拍摄界面以下区域的两台相机所得到的图像水平翻转，以便于将界面以下区域的图像和界面以下区域的图像进行对比。
[0082] 可以理解，在所述步骤S3中，识别初始时刻界面位置的方法有两种，第一种是采用人工识别界面位置后手动输入，第二种是使用算法识别。为了确保识别的准确性，本实施例优选采用第二种识别方法。具体地，如图5所示，所述步骤S3包括以下步骤：
[0083] 步骤S31：若气体运动速度远大于液体运动速度，在合适的采样频率下，每个像素的亮度取多张图片中相应位置亮度的最小者；若气体运动速度没有远大于液体运动速度，则跳过本步骤；
[0084] 步骤S32：对图片进行灰度形态学开操作若干次，消除图像中气体中的粒子而保留液体中的粒子；
[0085] 步骤S33：设置亮度阈值识别液体中的粒子，得到二值图像；
[0086] 步骤S34：对二值图像进行开操作若干次，再进行闭操作若干次；
[0087] 步骤S35：根据图像中的界面高度设置高度数值y0，使界面的高度位置在高度为该数值的直线以下，高度在y0以上的部分设置为0；
[0088] 步骤S36：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置，得到一个序列；
[0089] 步骤S37：对序列进行中值滤波和高斯滤波处理，得到界面的位置。
[0090] 可以理解，当在气液界面的两侧各自使用两台高速相机拍摄示踪粒子图像时，所述步骤S3还包括以下步骤：
[0091] 步骤S38：计算两台相机所得到的照片中界面位置的平均值。
[0092] 可以理解，在所述步骤S31中，如果气体的运动速度远大于液体的运动速度，即气体的平均运动速度比液体的平均运动速度大一个数量级以上，在合适的采样频率下，图片中气体中的示踪粒子有明显的移动，即每帧移动大于1像素，而液体中的示踪粒子没有明显的移动，即每帧移动小于1像素，则可以使用时间序列上3或5张图片的滑动最小值，即每个像素的亮度取3或5张图片相应位置亮度的最小者，再进行后续的步骤处理。否则，跳过本步骤。
[0093] 可以理解，在所述步骤S32中，灰度形态学开操作指的是先对图片进行灰度形态学腐蚀，再进行灰度形态学膨胀，从而可以消除图片中气体中的粒子而保留液体中的粒子。
[0094] 可以理解，在所述步骤S33中，由于粒子的亮度与其直径的平方有关，而气体中的粒子直径为1μm左右，而液体中的粒子直径为50μm左右，从而可以设置一个亮度阈值来识别液体中的粒子，例如1表示当地亮度高于阈值，0表示当地亮度低于阈值。并且，可以设置全局阈值，即空间各处使用相同大小的亮度阈值，也可以设置渐变的本地阈值，即空间中不同位置有不同的亮度阈值大小，这取决于各成像参数，若使用长焦镜头和大尺寸感光芯片，图像边缘区域整体亮度可能小于中心区域，这时可以选择减小边缘区域阈值。
[0095] 可以理解，在所述步骤S34中，对步骤S33中得到的二值图像先进行灰度形态学腐蚀再进行灰度形态学膨胀若干次，再进行闭操作若干次，即先进行灰度形态学膨胀再进行灰度形态学腐蚀若干次。
[0096] 可以理解，在所述步骤S35中，根据二值图像中的界面高度设置一个高度数值y0，以使界面的高度位置在高度为y0的直线以下，而高度在y0以上的部分设置为0。
[0097] 可以理解，在所述步骤S36中，记录下二值图像中每列像素的所有非零值中高度最高的一个像素的位置，如果某一列像素中没有非零值的存在，即该列像素的高度位置均在高度为y0的直线以上，则记录该列像素最后一个像素的位置，从而得到一个位置序列。
[0098] 可以理解，在所述步骤S38中，为了确保识别的准确性，采用计算两台相机所得到的照片中界面位置的平均值作为最终的初始界面位置。
[0099] 可以理解，由于本算法是通过识别液体中的粒子再进行平滑从而达到识别界面的目的，而液体中的示踪粒子始终在气液液面以下，故需要设置偏置量，然后再根据偏置量向上移动步骤S38中得到的初始界面位置，以便让识别到的界面位置和实际的界面位置重合，进一步提高了识别的准确性。该偏置量是通过实际图片多次试验得到的。因此，所述步骤S3在步骤S37之后还包括步骤S39：
[0100] 步骤S39：设置偏置量，然后根据偏置量向上移动初始界面位置。
[0101] 可以理解，如图6所示，在气液界面的两侧各自使用两台高速相机拍摄示踪粒子图像时，所述步骤S4具体包括以下步骤：
[0102] 步骤S41：若气体运动速度远大于液体运动速度，在合适的采样频率下，每个像素的亮度取多张图片中相应位置亮度的最小者；若气体运动速度没有远大于液体运动速度，则跳过本步骤；
[0103] 步骤S42：将两台相机所得到的同时刻的照片合成一张照片；
[0104] 步骤S43：对合成的照片进行灰度形态学开操作若干次，消除图像中气体中的粒子而保留液体中的粒子。
[0105] 可以理解，所述步骤S41与上述步骤S31一致，故在此不再赘述。
[0106] 可以理解，所述步骤S42具体为：将前一帧图像的界面向上平移10像素作为分界线，分界线以上区域比较两张照片同位置的亮度，取相对较低值作为合成照片的亮度；再将前一帧图像的界面向下平移10像素作为分界线，分界线以下区域比较两张照片同位置的亮度，取相对较高值作为合成照片的亮度；两条分界线中间区域可以从较低值到平均值再到较高值渐变，也可以取本地亮度的平均值。
[0107] 可以理解，所述步骤S43与上述步骤S32一致，故在此不再赘述。
[0108] 可以理解，当在气液界面的两侧各自使用一台高速相机拍摄示踪粒子图像时，所述步骤S42可以省略，所述步骤S43中直接对拍摄界面以下区域的相机所得到的图像进行灰度形态学开操作若干次。
[0109] 可以理解，如图7所示，所述步骤S5具体包括以下步骤：
[0110] 步骤S51：对于预处理得到的图像，设置亮度阈值识别液体中的粒子以得到二值图像；
[0111] 步骤S52：对二值图像进行开操作若干次，再进行闭操作若干次；
[0112] 步骤S53：将前一帧图像的界面向上平移15像素作为分界线，分界线以上区域二值图像数值设置为0；
[0113] 步骤S54：记录二值图像中每列像素的所有非零值中高度最高的位置，若某列没有非零值存在，则记录该列最后一个像素的位置和前一帧界面相同位置的较小值，从而得到一个序列；
[0114] 步骤S55：对序列进行中值滤波和高斯滤波得到界面的位置；
[0115] 步骤S56：根据偏置量向上移动界面的位置。
[0116] 可以理解，所述步骤S51与上述步骤S33一致，故在此不再赘述，其中，亮度阈值是根据该组试验拍摄界面以上区域的相机的亮度分布设定的。
[0117] 可以理解，所述步骤S52与上述步骤S34一致，故在此不再赘述。
[0118] 可以理解，在所述步骤54中，当某一列像素中没有非零值存在时，则记录该列最后一个像素的位置和前一帧界面相同位置的较小值，即取两者中的高度较高者。
[0119] 可以理解，步骤S55与上述步骤S37一致，故在此不再赘述。
[0120] 可以理解，在所述步骤S56中，该偏置量是在步骤S39中设定好的。
[0121] 可以理解，如图8所示，所述步骤S6具体包括以下步骤：
[0122] 步骤S61：相隔数帧计算每两帧之间界面位置的互相关值；
[0123] 步骤S62：将多个互相关值相加后寻找峰值位置；
[0124] 步骤S63：针对所有可计算的时刻重复步骤S61和步骤S62，得到界面速度的时间序列；
[0125] 步骤S64：基于得到的界面速度的时间序列对界面位置做锁相平均计算；
[0126] 步骤S65：对所有时刻的边界位置做时间-空间平滑操作，得到平滑后的边界位置时间-空间序列。
[0127] 可以理解，在所述步骤S61中，例如，按照时间先后顺序选择甲帧和乙帧，其中，乙帧界面位置固定，甲帧界面沿x方向移动，从而得到互相关值随甲帧界面移动距离之间的关系，然后寻找互相关峰值对应的移动距离。其中，相隔数帧计算互相关值的目的是考虑到界面移动速率比较慢，增加时间间隔可以提高精度。
[0128] 可以理解，在所述步骤S62中，为了进一步提高精度，可以选择将多个互相关值相加后寻找峰值位置。例如，选择甲一帧、甲二帧、甲三帧、乙一帧、乙二帧和乙三帧，其中甲二帧和乙二帧的时间间隔与甲一帧和乙一帧之间的时间间隔、甲三帧和乙三帧之间的时间间隔相同，分别计算上述三组双帧数据界面位置互相关值随移动距离之间的关系，再将三组互相关值相加得到新的界面位置互相关值随移动距离之间的关系，寻找相应互相关峰值对应的移动距离。另外，叠加计算的组数可以根据需要选定，在此不做具体限定。
[0129] 可以理解，在所述步骤S63中，单位时间内的界面移动距离即为界面速度，而界面移动距离根据步骤S61和步骤S62中互相关峰值的位置得到。针对所有可以计算的时刻重复步骤S61和步骤S62，从而得到界面速度的时间序列。作为优选的，所述步骤S63中还可以采用平滑滤波对界面速度的时间序列进行处理以减小时间序列脉动量，平滑滤波可以是中值滤波、高斯滤波、巴特沃斯滤波等。
[0130] 可以理解，当相机采用等时间间距曝光方式时，所述步骤S64具体为：记t0时刻的界面位置为，t0)>，其中表示界面高度，表示界面的横坐标序列，<>表示矢量，界面位置，t0)>即表示在t0时刻中某点x0所在列的界面高度在中对应的分量。记步骤S61中中计算互相关值的时间间隔为△t，经过平滑后得到t0时刻界面在△t内移动距离为<△x(t0)>。然后记N为时间-空间平滑半径，计算-N<Δx(t0)>，t0-NΔt)>，-(N-1)<Δx(t0)>，t0–(N-1)Δt)>，…，,t0)>，+(N-1)<Δx(t0)>，t0+(N-1)Δt)>，+N<Δx(t0)>，t0+NΔt)>在上的平均值，t0)>，其中表示将向两侧延伸后得到的新横坐标序列。当相机采用双帧曝光方式时，双帧之间的时间间隔为Δt，双帧曝光周期为ΔT，则t0时刻界面在△t内移动距离为<△x(t0)>，在ΔT内移动距离为<△X(t0)>＝<△x(t0)>*ΔT/Δt，平均算法相仿，对于双帧曝光的第一帧，计算-N<ΔX(t0)>，t0-NΔT)>，-(N-1)<ΔX(t0)>，t0–(N-1)ΔT)>，…，,t0)>，+(N-1)<ΔX(t0)>，t0+(N-1)ΔT)>，…,+N<ΔX(t0)>，t0+NΔT)>,-N<ΔX(t0)>+<△x(t0)>，t0–NΔT+Δt)>，-(N-1)<ΔX(t0)>+<△x(t0)>，t0–(N-1)ΔT+Δt)>，…，+<△x(t0)>,t0+Δt)>，+(N-1)<ΔX(t0)>+<△x(t0)>，t0+(N-1)ΔT+Δt)>，…,+N<ΔX(t0)>+<△x(t0)>，t0+NΔT+Δt)>的平均值，对于双帧曝光的第二帧，计算-N<ΔX(t0)>，t0-NΔT)>，-(N-1)<ΔX(t0)>，t0–(N-1)ΔT)>，…，,t0)>，+(N-1)<ΔX(t0)>，t0+(N-1)ΔT)>，…,+N<ΔX(t0)>，t0+NΔT)>,-N<ΔX(t0)>+<△x(t0)>，t0–NΔT-Δt)>，-(N-1)<ΔX(t0)>+<△x(t0)>，t0–(N-1)ΔT-Δt)>，…，+<△x(t0)>，t0-Δt)>，+(N-1)<ΔX(t0)>+<△x(t0)>，t0+(N-1)ΔT-Δt)>，…，+N<ΔX(t0)>+<△x(t0)>，t0+NΔT-Δt)>的平均值。另外，在所述步骤S64中需要进行一维插值，并且平均计算只使用内插得到的数值而不使用外插得到的数值。其中，N设置为30，以便于界面位置向两侧延伸足够的长度。
[0131] 可以理解，在所述步骤S65中，在边缘处进行时间-空间平滑操作时可以减小N。另外，在得到的平滑后的边界位置时间-空间序列中，缺少数据的部分可以使用插值得到数据。至此，已经得到了拍摄气液界面以上区域的两台相机所得图像中的界面位置。
[0132] 可以理解，拍摄界面以上区域的两台相机所得图像中的界面位置已经通过上述步骤得到了，现在需要通过已经得到的界面在物理空间中的位置来确定拍摄界面以下区域的相机所得到图像中的界面位置。如图9所示，所述步骤S7具体包括以下步骤：
[0133] 步骤S71：根据标定得到的坐标变换关系将拍摄界面以下区域的相机所得到的图像从原始图像变换到物理空间坐标系下的图像；
[0134] 步骤S73：对边界位置时间-空间序列进行缩放并插值以适应经过坐标变换的拍摄界面以下区域的相机所得到的图像；
[0135] 步骤S74：暂时给定界面位置的初始偏移量，在一张图片上显示相机经过空间变换的图像和界面经过偏移后暂时的位置，输出由不同时刻的图片组成的视频；
[0136] 步骤S75：观察视频中暂定的界面位置和由图片中粒子亮度反应的界面位置的差异，并根据差异修改界面位置的偏移量；
[0137] 步骤S76：重复执行步骤S74和步骤S75，直至界面暂时的位置和图片中粒子亮度反映的界面位置重合。
[0138] 可以理解，如图10所示，当在气液界面两侧各自使用两台相机进行拍摄时，所述步骤S7在步骤S1和步骤S3之间还包括以下步骤：
[0139] 步骤S72：检查坐标变换后的两台相机所得图像中的界面位置是否重合，若不重合，则对拍摄界面以下区域的两台相机进行自标定。具体地，利用肉眼检查拍摄界面以下区域的两台相机所得到的图像中的界面位置是否重合，界面位置可以通过图像中的粒子亮度大致判断，如果两者的界面位置不重合，则需要对拍摄界面以下区域的两台相机进行自标定，自标定的过程与上述步骤S23中的自标定过程一致，故在此不再赘述。
[0140] 可以理解，在所述步骤S71中，根据步骤S2中的标定得到的坐标变换关系将拍摄界面以下区域的相机所得到的图像从原始图像变换到物理空间坐标系下的图像。
[0141] 可以理解，在所述步骤S73中，由于对拍摄界面以上区域的相机的标定和对拍摄界面以下区域的相机的标定是各自单独进行的，两组标定关系中每毫米对应的像素数不一定相同，因此需要对得到的边界位置时间-空间序列进行缩放并插值以适应经过坐标变换后的拍摄界面以下区域的相机所得到的图像。
[0142] 可以理解，在所述步骤S74中，暂时给定界面位置的初始偏移量，然后在一张图片上显示相机经过空间变换的图像和界面经过偏移后暂时的位置，最后输出由不同时刻的图片组成的视频。可以理解，如图8所示，当在气液界面两侧各自使用两台相机进行拍摄时，在所述步骤S74中，两台相机的图像需分别占据不同的颜色通道，以便于区分。
[0143] 可以理解，在所述步骤S75中，由于图片中界面以上区域的粒子亮度比较大，而界面以下区域的粒子亮度比较小，因此，可以明显地观察出视频中暂定的界面位置和由图片中粒子亮度反应的界面位置的差异。
[0144] 可以理解，在所述步骤S76中，在界面暂时的位置和图片中粒子亮度反映的界面位置重合后，即得到了拍摄界面以下区域的相机所得到的图像中气液界面的位置。至此，所有相机所得到的图像中的气液界面位置均已经得到了。
[0145] 本申请的发明人还尝试过采用以下三个方法来识别气液界面位置：1、识别界面处反光的特征；2、在液体里添加荧光物质，识别存在荧光物质的区域；3、识别示踪粒子存在的区域。但是在第一种方法中，需要令相机拍摄的方向与界面平齐以促进捕捉界面反光，但是因为气液界面实际上是高低起伏的、且大部分情況下界面高度沿展向不保持恒定，这样就会有视线受到遮挡的问题，在使用立体式粒子图像测速时遮挡特别严重；而在第二种方法中，需要两种不同颜色的光源和三组拍摄不同颜色内容的相机，整套系统比较复杂，成本也很高，并且常用的荧光染料(罗丹明B)是有毒的，危害人体健康；而第三种方法难以适用于界面两侧都有粒子存在的情况，因为界面两侧都有粒子时，两侧的粒子不容易相互区分。因此，本申请的发明人提出了本申请优选实施例的识别气液界面位置的方法，首先分别在气体和液体中分撒示踪粒子，然后通过在气液界面的两侧各自使用高速相机拍摄示踪粒子图像，然后识别出拍摄界面以上区域的相机所得到的图像中的初始时刻的界面位置，然后再对所有同时刻的图像进行预处理，识别出随后各帧图像中的气液界面位置，再对得到的气液界面位置进行时间-空间平滑，从而得到了拍摄界面以上区域的相机所得到的图像中的所有气液界面位置，最后再根据拍摄界面以上区域的相机所得到的图像中的气液界面位置得到拍摄拍摄界面以下区域的相机所得到的图像中的气液界面位置。本发明的识别气液界面位置的方法，可以在气液界面两侧均有粒子存在的情况下准确识别气液界面的位置，且识别准确度高，可以很好地适用于粒子图像测速中的近壁测量过程。
[0146] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。