基于光场相机的三维流场的温度场和速度场同步测试方法

摘要

一种基于光场相机的三维流场的温度场和速度场同步测试方法，对磷光粒子进行温度和衰减时间的校准，获得该磷光粒子衰减时间和温度的精确关系；利用单光场相机获取待测流场中磷光粒子的时序粒子光场图像；预处理时序粒子光场图像，去除背景噪声；对时序粒子光场图像进行重构，获得时序三维粒子图像；对三维粒子图像进行互相关计算，获得三维流场速度场分布；对时序三维粒子图像进行衰减时间计算，并通过校准后的磷光粒子衰减时间和温度的精确关系获得对应的温度信息，获得三维流场温度场分布；对三维流场温度场和速度场分别进行后处理，对速度场和温度场结果进行修正。本发明同时测量三维流场的速度和温度分布，领先于世界先进水平。

说明书

技术领域

本发明涉及流场观测技术领域，特别涉及一种基于单光场相机的三维流场温度场和速度场的同步测试方法。

背景技术

流场观测技术一直是国家重点发展的核心技术，可以应用于航空航天、船舶驱动、电力生产和汽车工业等重要核心工业中。现阶段对于二维流场的温度场和速度场的研究已经较为成熟，但是三维流场的温度场和速度场的同步测试方法仍是世界性难题。

发明内容

针对现有技术中的缺陷，本发明的目的是提供一种基于光场相机的三维流场的温度场和速度场同步测试方法。本发明能通过单个光场相机的时序图像采集而同步获取待测流场的三维速度分布和三维温度场分布，相比于现有多视角三维流场测试方法，极大的减少了硬件系统配置、简化了硬件系统调节步骤，特别适用于受限空间下的三维流场测量；并且该方法开创性的同时测量三维流场的速度和温度分布，领先于世界先进水平。

本发明所采用的技术方案如下：

一种基于光场相机的三维流场的温度场和速度场同步测试方法，包含以下步骤：

A1，对温敏磷光粒子进行温度和衰减时间的校准，获得该磷光粒子衰减时间和温度的精确关系；

A2，利用单光场相机获取待测流场中磷光粒子的时序粒子光场图像；

A3，预处理时序粒子光场图像，去除背景噪声；

A4，采用CPU和/或GPU并行处理方法，对时序粒子光场图像进行重构，获得时序三维粒子图像；

A5，对第一帧和第三帧的时序三维粒子图像进行互相关计算，获得三维流场速度场分布；

A6，对三维流场速度场进行后处理，剔除错误速度矢量，同时根据相邻正确矢量插值替代被剔除的错误矢量；

A7，对第一帧和第二帧的时序三维粒子图像进行衰减时间计算，并通过校准后的磷光粒子衰减时间和温度的精确关系获得对应的温度信息，获得三维流场温度场分布；

A8，对三维流场温度场进行后处理，剔除错误温度信息，同时根据相邻正确温度信息插值替代被剔除的错误温度信息。

步骤A1中，在恒温流场中撒布磷光粒子，用特殊波段激光光源提供面式照明，温敏磷光粒子会受激光激发，发射另一种波长的磷光，该粒子发射磷光会持续一段时间，并且在这段时间内的光照强度不断衰减直至消失，这段时间被称为衰减时间，该衰减时间的长度会随着温度的变化而变化，通过在恒温流场中插入热电偶测量该流场的温度，利用高速相机或者PMT光电倍增管精确测量磷光粒子在该温度下的衰减时间，然后不断改变流场温度，调整至恒温，重复上述过程，最后获得该磷光粒子衰减时间和温度的精确关系，用以后续测量的校准。

步骤A2中，在待测流场中撒布温敏磷光粒子，用特殊波长激光光源提供体式照明，然后利用高分辨或者高速图像传感器与微透镜阵列的精密封装体形成的单光场相机拍摄磷光粒子的时序光场图像，该相机能在极短的时间内连续拍摄两帧图像，然后间隔较长的时间进行下一组连续两帧的拍摄，其中由于磷光粒子会持续发光，该发光过程持续时间包含在光场相机一组连续两帧的时间内，一组连续两帧图像用于温度场的计算，而分别属于不同两组中的连续两帧图像用于速度场的计算。

步骤A3中，对所拍摄的时序粒子光场图像进行预处理，其中去除背景噪声的降噪运算采用全局阈值、局部阈值、高斯光滑滤波、滑动最小值滤波算法的一种或多种算法组合。

步骤A4中，所述的重构包括：

根据几何光学计算空间光束经主透镜、微透镜阵列后到达图像传感器的准确位置，对空间每一个体素所发出的光束，计算其与相交微透镜单元的重叠面积，以及其与相交CCD/CMOS像素的重叠面积，从而计算获得空间每个体素与对应像素的权重系数W_i，j

式中各变量的含义如下：

Y_CCD为光线与图像传感器相交点的坐标；

V_B为光线与主镜头相交点的坐标；

s_y为微透镜阵列中某一微透镜单元相对主光轴的偏移值；

Y_l为光线与某一微透镜单元相交点的坐标；

s_i为主透镜与微透镜阵列之间的距离；

f_l为微透镜焦距；

f_m为主透镜焦距；

p为像素单元的尺寸；

M为放大系数，

利用密集光线逆追踪方法，定位空间重构体素所对应的像素集合，以每个像素为单元，根据式(1)来计算该像素所采集空间光线的位置，即反向追踪光线，以建立起像素与空间某一点光源所发出光线的对应关系；

利用流场磷光粒子的稀疏性，其对应所需重构的体素区域也为三维稀疏矩阵，将空间某一体素所对应的像素值相乘，如果乘积非零或者大于某一阈值，则表示该体素处可能存在一个磷光粒子，可以对其进行后续的重构计算；反之如果乘积为零或者小于某一阈值，则表示该体素处没有磷光粒子，后续重构计算中可以忽略该体素；

根据乘法代数重构算法，利用计算得到的权重系数W_i，j、以及所筛选出的非零体素，根据所记录的粒子光场像素矩阵I(x_i，y_i)，迭代计算稀疏重构矩阵的体素值E(Xj_，Y_j，Z_j)^k+1，

式中：

E(X_j,Y_j,Z_j)为空间体素(X_j,Y_j,Z_j)的数值；上标k表示第k次迭代计算所得到该体素的数值；I(x_i，y_i)为像素(x_i，y_i)处的数值，该数值由光场相机所拍摄的磷光粒子光场图像获得；w_i,j为体素j与像素i所对应的权重系数；μ为MART计算参数。

步骤A5中，将第一帧和第三帧的磷光粒子三维体素矩阵按照式(3)进行三维互相关计算，获取流场三维速度场分布，

其中：E₁(i,j,k)表示在三维体素坐标为(i,j,k)处的数值；

E₂(i+m,j+n,k+l)表示在三维体素(i+m,j+n,k+l)处的数值；

m、n、l分别第三帧粒子相对于第一帧对应粒子的相对坐标位移；

M、N、L分别表示该重构的三维体素在x、y、z三个方向上的最大值。

步骤A6中，采用全局阈值计算式(4)、局部中值滤波式(5)或者局部平均值滤波式(6)对获得的三维流场的速度场U(i，j，k)，挑选并剔除出其中的错误矢量，并采用线性插值或者三次样条插值的方法，根据正确速度矢量插值获得替代矢量，

其中，是三维速度场的平均值；

STD(U(i，j，k))是三维速度场的标准方差；

TH为滤波阈值；

是局部(3×3×3)区域内的三维速度场平均值；STD(U(i-1：i+1，j-1：j+1，k-1：k+1))是局部(3×3×3)区域内的三维速度场标准方差；

U(i-1：i+1，j-1：j+1，k-1：k+1)是局部(3×3×3)区域内的三维速度场中值。

步骤A7中，将第一帧和第三帧的磷光粒子三维体素矩阵按照式(7)进行衰减时间的计算，通过A1中的校准结果，获得流场三维流场温度场分布，

λ(i，j，k)＝-(ln((R(i，j，k)×e^-5+1)/(R(i，j，k)+1))÷t₁)(7)

其中，λ(i，j，k)是三维流场中对应体素在该温度下衰减时间τ(i，j，k)的倒数；

R(i，j，k)是三维流场中第一帧对应体素E₁(i，j，k)和第二帧对应体素E₂(i，j，k)的比值；

t₁是光场相机连续两帧拍摄过程中第一帧的曝光时间。

步骤A8中，采用全局阈值计算式(8)、局部中值滤波式(9)或者局部平均值滤波式(10)对获得的三维流场温度场T(i，j，k)，挑选并剔除出其中的错误温度，并采用线性插值或者三次样条插值的方法，根据正确温度插值获得替代温度，

其中，是三维温度场的平均值；

STD(T(i，j，k))是三维温度场的标准方差；

TH为滤波阈值；

是局部(3×3×3)区域内的三维温度场平均值；STD(T(i-1：i+1，j-1：j+1，k-1：k+1))是局部(3×3×3)区域内的三维温度场标准方差；

T(i-1：i+1，j-1：j+1，k-1：k+1)是局部(3×3×3)区域内的三维温度场中值。

本发明通过单个光场相机完成了对三维流场的速度场和温度场的同步测量，获得了较为精确的三维流场的温度分布和速度矢量分布情况，具有开创性意义。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：

图1是本发明的利用光电倍增管测量磷光粒子温度和衰减时间原理示意图；

图2是本发明的单光场相机三维流场的温度场和速度场同步测试系统原理示意图；

图3是本发明实施例中光线追踪示意图；

图4是本发明实施例中权重系数计算示意图a；

图5是本发明实施例中权重系数计算示意图b；

图6是本发明实施例中密集光线逆追踪方法原理示意图；

图7是本发明实施例中同步测量三维流场时序示意图；

图8是本发明实施过程流程图。

具体实施方式

下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。

图1实施例为本发明利用PMT光电倍增管模块对温敏磷光粒子Mg3F2GeO4:Mn(MFG)进行温度和衰减时间的校准，获得该磷光粒子衰减时间和温度的精确关系的系统示意图，包括一个光电倍增管模块1，一个恒温水浴箱2，一个磁力搅拌器3、一台紫外激光器4、一个具有水浴壁7的水箱、一个显示仪表5，图中阴影部分为测量区域8，一个旋转子10，受磁力搅拌器驱动后，在水箱内进行旋转，产生旋转流场。本系统通过355nm紫外激光激发磷光粒子9，粒子受激发发射660nm红光。该粒子发射红光会持续一段时间，并且在这段时间内的光照强度不断衰减直至消失。这段时间被称为衰减时间，该衰减时间的长度会随着温度的变化而变化。通过在恒温流场中插入热电偶6测量该流场的温度，利用PMT光电倍增管精确测量磷光粒子在该温度下的衰减时间，然后不断改变水浴箱温度，调整整个水浴系统至恒温，重复上述过程，最后获得该磷光粒子衰减时间和温度的精确关系，用以后续测量的校准。

图2为本发明中基于光场相机单光场相机三维流场的温度场和速度场同步测试系统原理示意图。类似与上述校准过程，在测量区域8中撒布磷光粒子，用355nm紫外激光光源4和透镜阵列11提供体式照明，受激发粒子发出的光通过主透镜12，然后利用高分辨或者高速图像传感器15与微透镜阵列(MLA)14的精密封装体形成的单光场相机拍摄磷光粒子的时序光场图像13。

在拍摄结束后，对原始时序粒子光场图片进行预处理，消除背景噪声和相机自身的误差。通过采用全局阈值、局部阈值、高斯光滑滤波、滑动最小值滤波算法的一种或多种算法组合，进行原始图片的降噪处理，提高测试结果的准确性。

在三维粒子重构方面，参阅图3至图5。由于光场相机本身结构的特殊性，其主透镜12和高分辨或者高速图像传感器(CCD/CMOS)15之间还有微透镜阵列14，所以光场相机可以通过一个视角拍摄出具有三维信息的图像。根据几何光学计算空间光束经主透镜、微透镜阵列后到达图像传感器的准确位置，对空间每一个体素(voxel)所发出的光束，计算其与相交微透镜单元的重叠面积，以及其与相交CCD/CMOS像素的重叠面积，从而计算获得空间每个体素与对应像素的权重系数W_i，j。

利用密集光线逆追踪方法，定位空间重构体素(Voxel)所对应的像素集合，以每个像素为单元，根据式(1)来计算该像素所采集空间光线的位置，即反向追踪光线，以建立起像素与空间某一点光源所发出光线的对应关系。利用流场磷光粒子的稀疏性，其对应所需重构的体素区域也为三维稀疏矩阵16，将空间某一体素所对应的像素值相乘，如果乘积非零或者大于某一阈值，则表示该体素处可能存在一个磷光粒子，可以对其进行后续的重构计算；反之如果乘积为零或者小于某一阈值，则表示该体素处没有磷光粒子，后续重构计算中可以忽略该体素。

根据乘法代数重构算法，利用计算得到的权重系数W_i，j、以及所筛选出的非零体素，根据所记录的粒子光场像素矩阵I(x_i，y_i)，迭代计算稀疏重构矩阵的体素值，最后重构出粒子的三维图像16。

参阅图6、图7，该相机能在极短的时间内(5ns)连续拍摄两帧图像，然后间隔较长的时间后进行下一组连续两帧的拍摄。其中由于磷光粒子会持续发光，该发光过程持续时间包含在光场相机一组连续两帧的时间内(如第一帧和第二帧)，一组连续两帧图像(第一帧和第二帧)完整的记录了磷光粒子在衰减过程中的光照强度信息，根据式(7)计算所测量三维区域对应体素的衰减时间，根据校准得到的温度和衰减时间的关系，计算得到三维流场的温度分布。

而第一帧和第三帧图像(分别属于不同两组连续两帧)用于速度场的计算。由于连续两帧(第一帧和第二帧)时间间隔只有5ns，时间太过于短暂无法用于速度场的计算，故使用第一帧和第三帧图像进行互相关计算，通过式(3)计算得到两帧图像对应粒子的速度变化，从而得到三维流场的速度矢量分布。

最后对初步得到的三维速度矢量场和温度场分别进行后处理，采用全局阈值计算、局部中值滤波或者局部平均值滤波式对获得的三维流场温度场T(i，j，k)和速度场U(i，j，k)处理，分别挑选并剔除出其中的错误温度和错误速度矢量，并采用线性插值或者三次样条插值的方法，根据正确温度插值获得替代温度，根据正确速度矢量插值获得替代速度矢量。

最后获得被测三维流场的同步的温度场和速度矢量场的三维分布结果。

图8示意了本发明方法的实施过程流程图。

以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。