CN113466489A - 一种低粒子密度的单相机粒子图像测速方法 - Google Patents

一种低粒子密度的单相机粒子图像测速方法 Download PDF

Info

Publication number
CN113466489A
CN113466489A CN202110631589.5A CN202110631589A CN113466489A CN 113466489 A CN113466489 A CN 113466489A CN 202110631589 A CN202110631589 A CN 202110631589A CN 113466489 A CN113466489 A CN 113466489A
Authority
CN
China
Prior art keywords
particle
dimensional
image
field
probability
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110631589.5A
Other languages
English (en)
Inventor
单良
施飞杨
李浩然
孔明
熊俊哲
刘维
洪波
郭天太
赵军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Jiliang University
Original Assignee
China Jiliang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Jiliang University filed Critical China Jiliang University
Priority to CN202110631589.5A priority Critical patent/CN113466489A/zh
Publication of CN113466489A publication Critical patent/CN113466489A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/18—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance
    • G01P5/20—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the time taken to traverse a fixed distance using particles entrained by a fluid stream
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)

Abstract

本发明提供一种低粒子密度的单相机粒子图像测速方法，包括以下步骤：在流场中均匀地撒入示踪粒子，白色光源发出一束准直的稳定白色光束，光束经过滤波片后得到彩色片状光射入流场，照亮流场中的示踪粒子，CCD摄像机拍摄得到低粒子密度的彩色粒子图片并传送给上位机，上位机对低粒子密度的彩色粒子图片进行滤波处理和数据降采样，对处理后的低粒子密度的彩色粒子图片进行分析重建得到高粒子密度的三维粒子分布概率场后进行互相关计算得到流体三维速度场。本发明解决了现有技术中三维互相关算法应用在低粒子密度条件时存在的问题。

Description

一种低粒子密度的单相机粒子图像测速方法
技术领域
本发明涉及测量技术领域，具体涉及一种低粒子密度的单相机粒子图像测速方法。
背景技术
流体运动在自然界中非常普遍，在飞机和水下航行器的设计中涉及的流体力学和流体的测量息息相关，因此掌握流体的具体运动情况是流体力学深入研究的前提条件，测量和分析流体的运动状态具有重要意义，是流体领域的关键问题。流场测量技术最早出现在1904年，普朗特发明了手摇式水槽来观测不同模型周围的流场，使用流场显示技术来对流场进行定性的分析描述，对流场的运动状态进行研究。随着科学技术的发展，可以通过温度、激光以及超声波等不同的测量手段对流场的速度进行定量的测量。
传统的流场测量方式中热线/热膜测速仪(HWFA)利用热平衡原理根据不同速度的流场流动从恒温的热敏元件带走的热量值差异可以计算得到流场的速度但是会对流场的运动有一定程度上的影响。传统的非接触式流体测量方式有激光多普勒测速仪(LDV)通过流体流经两束激光交汇处时对干涉条纹的影响计算粒子的速度，以及超声波多普勒测速仪(ADV)采用从声波发射换能器产生的声学脉冲经过测量点后回波的多普勒频移来判断流场速度的大小。上述流体测量方式虽然测量精度高且具有较高的分辨率，但是他们都只能进行单点测量无法对流场整体运动情况进行分析，难以获得流场的瞬态图像以及流场的整体结果。
粒子图像测速法发展于19世纪80年代，集现代材料、数字成像、激光技术和图像分析等领域的发展成果于一体，是一种瞬态、多点、无接触式的流体力学测速方法，可以精确地测量平面内的瞬态流场。和上述传统的测速方法相比粒子图像测速法可以在不干扰流体运动的情况为流体运动的定性描述以及定量研究提供相对理想的数据基础。粒子图像测速法主要通过在流场中散入具有良好跟随性的粒子，通过测量粒子的速度来反应流场的运动情况。用照明装置照亮流场中的粒子再通过成像系统拍摄多次曝光的粒子图片，通过PIV图像处理方法对粒子图片进行处理得到粒子的运动。粒子图像测速已经是目前流场速度测量的主要测量方法。其中最关键的就是粒子图像分析算法。因此对PIV图像处理算法的研究具有非常重要的意义。
现有技术中三维PIV算法在还原速度场的时候大部分都是根据拍摄的粒子图像重建得到对应的三维粒子场，然后通过三维互相关算法对三维粒子场进行互相关得到粒子的速度场。进行三维互相关运算时需要检测窗口内存在足够数量的粒子才能得到准确的速度结果。因此在低粒子密度流场测量时误差较大，需要对现有技术进行改进。
发明内容
本发明要解决的技术问题是提供一种高效的低粒子密度的单相机粒子图像测速方法。
为解决上述技术问题，本发明提供一种低粒子密度的单相机粒子图像测速方法，包括以下步骤：在流场中均匀地撒入示踪粒子，白色光源发出一束准直的稳定白色光束，光束经过滤波片后得到彩色片状光射入流场，照亮流场中的示踪粒子，CCD摄像机拍摄得到低粒子密度的彩色粒子图片并传送给上位机，上位机对低粒子密度的彩色粒子图片进行滤波处理和数据降采样，对处理后的低粒子密度的彩色粒子图片进行分析重建得到高粒子密度的三维粒子分布概率场后进行互相关计算得到流体三维速度场。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：通过不同波长的可见光来标定流场中粒子的深度，结合CCD摄像机拍摄得到图片中的二维像素位置信息可以得到流场的三维信息。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：将粒子在流场中的某一点的存在建模为存在概率并根据光学系统的成像过程建立对应的图像形成模型，得到线性系统的最小化问题，并用ADMM算法求解；通过低粒子密度的彩色粒子图像得到高粒子密度的三维粒子分布概率场。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：对得到的高粒子密度的三维粒子分布概率场进行三维互相关分析进而得到流体三维速度场。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：
波长λ和深度Z存在线性关系，通过波长λ确定粒子在三维立体空间中光平面的位置；联合通过彩色图像获得的粒子在该在光平面中的二维坐标(x，y)，(x，y，λ)来表示粒子的三维坐标；将三维空间中某一点的粒子存在的概率建模为P(X，λ)；由于采用的是非相干光，因此可以将光学系统的成像过程建模为一组点扩散函数(PSF)，PSF是聚焦光学系统的冲击响应，是成像系统传递函数的空间域表达；在这里将PSF分成红绿蓝三个颜色通道，分别建立对应的PSF点扩散函数：gC(x，λ)，C∈{red，green，blue}；因此可以将图像形成模型建模为：
iC(X)＝∫Λ∫XgC(X-X′，λ)·ir(X，λ)·P(X，λ)dX′dλ
式中iC(X)对应捕获的RGB图像的颜色通道，ir(X，λ)为入射在图像传感器上的相应光谱分布，gC(X-X’λ)为三维空间中某一点对应的PSF函数，P(X，λ)三维空间中对应点存在粒子的概率，dX为当前点的光平面上的二维位置dλ为光平面对应的波长；空间积分对应于代表潜在不完美焦点的卷积；通过波长上的积分可将编码3D粒子位置的光谱图像(图像传感器上获得的光谱分布)转换为RGB图像；
离散化后可以将PSF和反射光强度的卷积公式化为矩阵A∈R3N×NL，其中N是低粒子密度的彩色粒子图像的像素数量，L是沿波长坐标方向的离散化级数，数字3表示三个不同颜色的图像通道；it∈R3N表示在t时间相机拍摄到的粒子图像；Pt∈[0，1]NL表示粒子在时间点t时在不同点存在的概率；通过对线性系统的求解可以得到每个时间点的粒子场分布；APt＝it但是编码粒子位置的完整光谱信息存在于三个不同的颜色通道中，导致这个线性系统求解是一个病态的逆问题；因此需要将一些粒子分布的先验知识作为正则项引入，将逆问题变成下述最小化问题：
Figure BDA0003103761770000031
式中：(p*)为求解最小化问题得到的三维空间中每一点存在粒子的概率大小的解的集合，APt为根据图像形成模型得到的模拟粒子图片；it为t时刻采集的低粒子密度彩色粒子图片；A[p1|...|pT]为对一组拍摄的图片进行分析时根据图像形成模型得到的一组模拟粒子图片；pt为t时刻的粒子存在概率；k1为优化参数；
运算符П[0，1]将所有体积占用概率投影到有效概率[0，1]NL上；公式中的第一行是求解线性系统逆问题的最小二乘数据拟合项；第一项定义了加权项L1，鼓励待测流场中粒子的稀疏分布，由于相机对不同波长的光的灵敏度是不同的，因此通过对角矩阵diag(w)进一步加权，根据先前的结果迭代地更改了权重系数，权重系数在迭代过程中是固定不变的，但会随着粒子深度的变化而改变，可以有效地增强稀疏性；在光谱的黄色或蓝绿色部分的波长在多个颜色通道中引起强烈响应时，远离蓝色或距离红色部分较远的波长的光只触发一个通道，这将导致颗粒分布的不均匀；通过加权项的设定可以将不同的粒子放置在其更有可能处在的深度；加权项通过补偿光度的不均匀性来消除这种偏差；公式第二项将指示符函数占用概率控制在[0，1]之间；
上述优化问题中由于加权项L1和指标函数П[0，1](p1；...；pt)的存在，该优化问题是非平滑的；采用交替方向乘子法ADMM来解决该非平滑问题；交替方向乘子法ADMM通过分解协调程序的形式将较大的难以求解的全局问题解耦为多个较小并且易于得到有效解的局部问题，通过对多个子问题进行迭代求解进而得到全局问题的最终解。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：
在使用ADMM算法求解优化问题后得到粒子分布的概率场后将前一个三维粒子分布场根据设定的搜索窗口大小分割成多个大小固定的搜索区域；并对后一个粒子场中同位置的搜索区域在比搜索窗口更大的问询窗口内进行左上角逐行扫描至右下角根据公式
Figure BDA0003103761770000041
Figure BDA0003103761770000042
进行3D互相关，可以得到一个相关系数的矩阵，其中相关系数的峰值的坐标即可表示粒子位移的距离；在将拍摄的两个数据集的不同时刻带入公式
Figure BDA0003103761770000043
Figure BDA0003103761770000044
可得到流场的粗略三维速度场；然后根据流体力学的基本原理对粗大误差进行剔除，并通过插值拟合对错误速度矢量进行矫正，得到精确的三维粒子速度场。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：所述滤波片为线性可变带通滤波器。
作为对本发明一种低粒子密度的单相机粒子图像测速方法的改进：粒子溶液和流体溶液体积比小于100。
现有技术中三维PIV算法在还原速度场的时候大部分都是根据拍摄的粒子图像重建得到对应的三维粒子场，然后通过三维互相关算法对三维粒子场进行互相关得到粒子的速度场。进行三维互相关运算时需要检测窗口内存在足够数量的粒子才能得到准确的速度结果。因此在低粒子密度流场测量时误差较大，本发明提供一种低粒子密度的单相机粒子图像测速方法，在低粒子密度的情况下也能得到准确的测量结果。
本发明一种低粒子密度的单相机粒子图像测速方法的技术优势为：
本发明将粒子在流场中某一点的存在的可能性建模为存在概率并根据光学系统的成像过程建立对应的图像形成模型，得到线性系统的小化问题并用ADMM算法求解。通过低粒子密度的彩色粒子图像得到高粒子密度的三维粒子分布概率场。提高三维互相关算法在低粒子密度条件下也能得到良好的实验结果。
本发明解决了现有技术中三维互相关算法应用在低粒子密度条件时存在的问题。本发明有益效果为：利用彩色光源代替激光光源，并用不同波长的可见光来标定粒子的深度，结合CCD相机拍摄的二维像素信息可以用一台相机得到流场的三维信息，减少实际应用的局限性。对CCD拍摄的低粒子密度彩色粒子图片进行分析求解得到高密度的三维粒子分布概率场，通过三维互相关算法对概率场分析得到三维流体速度场。
具体实施方式
下面结合具体实施例对本发明进行进一步描述，但本发明的保护范围并不仅限于此。
实施例1、一种低粒子密度的单相机粒子图像测速方法，主要由两个部分组成：
一个是光学系统，将粒子在三维空间中的深度编码为对应的不同颜色，并对流体中粒子的彩色图像进行采集。
另一个部分是重建算法，通过对光学系统拍摄的彩色图像进行分析，分步还原三维粒子场以及流体的三维速度场。
光学装置：与常规粒子图像测速PIV相比相同的是，都采用片状光照亮流体中的粒子，再通过相机对流体中流动的粒子进行拍摄，得到低粒子密度的彩色图像。不同的是常规粒子图像测速PIV往往使用激光光源，通过柱透镜得到垂直于相机视线的片状光。可以将相机聚焦在照明平面上，观察跟踪照明平面中移动的粒子可以得到当前体积的2D切片上的速度场的两个分量。而用白色光源代替了激光光源，并经过准直和滤波之后可以得到波长随深度线性变化的片状光，因此在同样的2D切片上可以同时获得速度场的三个分量。
在流场中均匀地撒入少量的示踪粒子(粒子溶液和流体溶液体积比小于100)，白色光源发出一束准直的稳定白色光束，光束经过滤波片(线性可变带通滤波器)后得到彩色片状光射入流场，照亮流场中的示踪粒子CCD摄像机拍摄得到图片并传送给上位机。
重建算法：在对光学装置的得到低粒子密度的彩色图像进行高斯滤波和数据降采样之后，再进行计算分析，根据图像上粒子的位置以及颜色得到高粒子密度的三维粒子分布概率场。得到间隔固定时间T的两个三维粒子场，用3D互相关算法对其进行分析重建，得到流体三维速度场。
重建粒子的三维粒子分布概率场需要对图像形成过程进行建模，将模型得到的粒子位置与拍摄的低粒子密度的彩色图像相关联。对图像形成模型进行求解可得到粒子的三维空间中每个点存在粒子的概率，从低密度的三维粒子场得到高粒子密度的三维粒子分布概率场。但是得到的是一个不适定的逆问题，引入两个正则项将其进行优化，保证可以收敛有效地解决该问题。其中一个正则项是保证粒子稀疏分布的k1||diag(w)(p1；...；pT)||1，另一个是保证将三维空间中某一点存在粒子概率在[0，1]有效凸集上。
因为光学装置采用的光照为连续的窄带非相关光谱光片，其波长λ随深度的变化而变化。波长λ和深度Z存在线性关系。可以通过波长λ确定粒子在三维立体空间中光平面的位置。联合通过彩色图像获得的粒子在该在光平面中的二维坐标(x，y)，(x，y，λ)来表示粒子的三维坐标。将三维空间中某一点的粒子存在的概率建模为P(X，λ)。由于采用的是非相干光，因此可以将光学系统的成像过程建模为一组点扩散函数(PSF)，PSF是聚焦光学系统的冲击响应，是成像系统传递函数的空间域表达。在这里将PSF分成红绿蓝三个颜色通道，分别建立对应的PSF点扩散函数：gC(x，λ)，C∈{red，green，blue}。因此可以将图像形成模型建模为：
iC(X)＝∫Λ∫XgC(X-X′，λ)·ir(X，λ)·P(X，λ)dX′dλ
式中iC(X)对应捕获的RGB图像的颜色通道，ir(X，λ)为入射在图像传感器上的相应光谱分布，gC(X-X’，λ)为三维空间中某一点对应的PSF函数，P(X，λ)三维空间中对应点存在粒子的概率，dX为当前点的光平面上的二维位置dλ为光平面对应的波长。空间积分对应于代表潜在不完美焦点的卷积。通过波长上的积分可将编码3D粒子位置的光谱图像(图像传感器上获得的光谱分布)转换为RGB图像。
离散化后可以将PSF和反射光强度的卷积公式化为矩阵A∈R3N×NL，其中N是低粒子密度的彩色粒子图像的像素数量，L是沿波长坐标方向的离散化级数，数字3表示三个不同颜色的图像通道。it∈R3N表示在t时间相机拍摄到的粒子图像。Pt∈[0，1]NL表示粒子在时间点t时在不同点存在的概率。通过对线性系统的求解可以得到每个时间点的粒子场分布。APt＝it但是编码粒子位置的完整光谱信息存在于三个不同的颜色通道中，导致这个线性系统求解是一个病态的逆问题。因此需要将一些粒子分布的先验知识作为正则项引入，将逆问题变成下述最小化问题：
Figure BDA0003103761770000061
式中：(p*)为求解最小化问题得到的三维空间中每一点存在粒子的概率大小的解的集合，APt为根据图像形成模型得到的模拟粒子图片；it为t时刻采集的低粒子密度彩色粒子图片；A[p1|...|pT]为对一组拍摄的图片进行分析时根据图像形成模型得到的一组模拟粒子图片；pt为t时刻的粒子存在概率；k1为优化参数。
运算符П[0，1]将所有体积占用概率投影到有效概率[0，1]NL上。公式中的第一行是求解线性系统逆问题的最小二乘数据拟合项。第一项定义了加权项L1，鼓励待测流场中粒子的稀疏分布，由于相机对不同波长的光的灵敏度是不同的，因此通过对角矩阵diag(w)进一步加权，根据先前的结果迭代地更改了权重系数，权重系数在迭代过程中是固定不变的，但会随着粒子深度的变化而改变，可以有效地增强稀疏性。在光谱的黄色或蓝绿色部分的波长在多个颜色通道中引起强烈响应时，远离蓝色或距离红色部分较远的波长的光只触发一个通道，这将导致颗粒分布的不均匀。通过加权项的设定可以将不同的粒子放置在其更有可能处在的深度。加权项通过补偿光度的不均匀性来消除这种偏差。公式第二项将指示符函数占用概率控制在[0，1]之间。
上述优化问题中由于加权项L1和指标函数П[0，1](p1；...；pt)的存在，该优化问题是非平滑的。因此无法通过常规优化方法如梯度下降法等工具进行解决。要解决非平滑的优化问题，可以将非平滑项与原始优化问题解耦，便于将不同的部分分别处理。本发明采用交替方向乘子法(Alternating Direction Method of Multipliers，ADMM)来解决该非平滑问题。ADMM通过分解协调程序的形式将较大的难以求解的全局问题解耦为多个较小并且易于得到有效解的局部问题，通过对多个子问题进行迭代求解进而得到全局问题的最终解。
ADMM框架重建粒子分布场
输入：目标函数F1，H1，迭代次数j，松弛变量z，最大迭代次数max，p
输出：对偶变量q
Figure BDA0003103761770000071
其中：F1和H1定义为：
Figure BDA0003103761770000072
H1(p)＝k1||diag(w)(p1；...；pT)||1+П[0，1](p1；...；pt)
近端运算符的推导：简化公式将zj-qj简化为dj，pj+1-qj简化为ej。近端运算符
Figure BDA0003103761770000073
是一个共轭梯度求解的过程：
Figure BDA0003103761770000081
近端运算符
Figure BDA0003103761770000082
是逐点收缩操作然后投影到[0，1]的域上：z＝П[0，1]((ej-τ1k1w))-(-ej-τ1k1w))
在使用ADMM算法求解优化问题后得到粒子分布的概率场后将前一个三维粒子分布场根据设定的搜索窗口大小分割成多个大小固定的搜索区域。并对后一个粒子场中同位置的搜索区域在比搜索窗口更大的问询窗口内进行左上角逐行扫描至右下角根据公式
Figure BDA0003103761770000083
Figure BDA0003103761770000084
进行3D互相关，可以得到一个相关系数的矩阵，其中相关系数的峰值的坐标即可表示粒子位移的距离。在将拍摄的两个数据集的不同时刻带入公式
Figure BDA0003103761770000085
Figure BDA0003103761770000086
可得到流场的粗略三维速度场。然后根据流体力学的基本原理对粗大误差进行剔除，并通过插值拟合对错误速度矢量进行一定的矫正，得到精确的流体三维速度场。
最后，还需要注意的是，以上列举的仅是本发明的若干个具体实施例。显然，本发明不限于以上实施例，还可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想到的所有变形，均应认为是本发明的保护范围。

Claims (8)

1.一种低粒子密度的单相机粒子图像测速方法，其特征在于：包括以下步骤：在流场中均匀地撒入示踪粒子，白色光源发出一束准直的稳定白色光束，光束经过滤波片后得到彩色片状光射入流场，照亮流场中的示踪粒子，CCD摄像机拍摄得到低粒子密度的彩色粒子图片并传送给上位机，上位机对低粒子密度的彩色粒子图片进行滤波处理和数据降采样，对处理后的低粒子密度的彩色粒子图片进行分析重建得到高粒子密度的三维粒子分布概率场后进行互相关计算得到流体三维速度场。
2.根据权利要求1所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：通过不同波长的可见光来标定流场中粒子的深度，结合CCD摄像机拍摄得到图片中的二维像素位置信息可以得到流场的三维信息。
3.根据权利要求1所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：将粒子在流场中的某一点的存在建模为存在概率并根据光学系统的成像过程建立对应的图像形成模型，得到线性系统的最小化问题，并用ADMM算法求解；通过低粒子密度的彩色粒子图像得到高粒子密度的三维粒子分布概率场。
4.根据权利要求1或2或3所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：对得到的高粒子密度的三维粒子分布概率场进行三维互相关分析进而得到流体三维速度场。
5.根据权利要求3所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：
波长λ和深度Z存在线性关系，通过波长λ确定粒子在三维立体空间中光平面的位置；联合通过彩色图像获得的粒子在该在光平面中的二维坐标(x，y)，(x，y，λ)来表示粒子的三维坐标；将三维空间中某一点的粒子存在的概率建模为P(X，λ)；由于采用的是非相干光，因此可以将光学系统的成像过程建模为一组点扩散函数(PSF)，PSF是聚焦光学系统的冲击响应，是成像系统传递函数的空间域表达；在这里将PSF分成红绿蓝三个颜色通道，分别建立对应的PSF点扩散函数：gC(x，λ)，C∈{red，green，blue}；因此可以将图像形成模型建模为：
iC(X)＝∫Λ∫XgC(X-X′，λ)·ir(X，λ)·P(X，λ)dX′dλ
式中iC(X)对应捕获的RGB图像的颜色通道，ir(X，λ)为入射在图像传感器上的相应光谱分布，gC(X-X′，λ)为三维空间中某一点对应的PSF函数，P(X，λ)三维空间中对应点存在粒子的概率，dX′为当前点的光平面上的二维位置dλ为光平面对应的波长；空间积分对应于代表潜在不完美焦点的卷积；通过波长上的积分可将编码3D粒子位置的光谱图像(图像传感器上获得的光谱分布)转换为RGB图像；
离散化后可以将PSF和反射光强度的卷积公式化为矩阵A∈R3N×NL，其中N是低粒子密度的彩色粒子图像的像素数量，L是沿波长坐标方向的离散化级数，数字3表示三个不同颜色的图像通道；it∈R3N表示在t时间相机拍摄到的粒子图像；Pt∈[0，1]NL表示粒子在时间点t时在不同点存在的概率；通过对线性系统的求解可以得到每个时间点的粒子场分布；APt＝it但是编码粒子位置的完整光谱信息存在于三个不同的颜色通道中，导致这个线性系统求解是一个病态的逆问题；因此需要将一些粒子分布的先验知识作为正则项引入，将逆问题变成下述最小化问题：
Figure FDA0003103761760000021
式中：(p*)为求解最小化问题得到的三维空间中每一点存在粒子的概率大小的解的集合，APt为根据图像形成模型得到的模拟粒子图片；it为t时刻采集的低粒子密度彩色粒子图片；A[p1|...|pT]为对一组拍摄的图片进行分析时根据图像形成模型得到的一组模拟粒子图片；pt为t时刻的粒子存在概率；k1为优化参数；
运算符П[0，1]将所有体积占用概率投影到有效概率[0，1]NL上；公式中的第一行是求解线性系统逆问题的最小二乘数据拟合项；第一项定义了加权项L1，鼓励待测流场中粒子的稀疏分布，由于相机对不同波长的光的灵敏度是不同的，因此通过对角矩阵diag(w)进一步加权，根据先前的结果迭代地更改了权重系数，权重系数在迭代过程中是固定不变的，但会随着粒子深度的变化而改变，可以有效地增强稀疏性；在光谱的黄色或蓝绿色部分的波长在多个颜色通道中引起强烈响应时，远离蓝色或距离红色部分较远的波长的光只触发一个通道，这将导致颗粒分布的不均匀；通过加权项的设定可以将不同的粒子放置在其更有可能处在的深度；加权项通过补偿光度的不均匀性来消除这种偏差；公式第二项将指示符函数占用概率控制在[0，1]之间；
上述优化问题中由于加权项L1和指标函数П[0，1](p1；...；pt)的存在，该优化问题是非平滑的；采用交替方向乘子法ADMM来解决该非平滑问题；交替方向乘子法ADMM通过分解协调程序的形式将较大的难以求解的全局问题解耦为多个较小并且易于得到有效解的局部问题，通过对多个子问题进行迭代求解进而得到全局问题的最终解。
6.根据权利要求5所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：
在使用ADMM算法求解优化问题后得到粒子分布的概率场后将前一个三维粒子分布场根据设定的搜索窗口大小分割成多个大小固定的搜索区域；并对后一个粒子场中同位置的搜索区域在比搜索窗口更大的问询窗口内进行左上角逐行扫描至右下角根据公式
Figure FDA0003103761760000031
Figure FDA0003103761760000032
进行3D互相关，可以得到一个相关系数的矩阵，其中相关系数的峰值的坐标即可表示粒子位移的距离；在将拍摄的两个数据集的不同时刻带入公式
Figure FDA0003103761760000033
Figure FDA0003103761760000034
可得到流场的粗略三维速度场；然后根据流体力学的基本原理对粗大误差进行剔除，并通过插值拟合对错误速度矢量进行矫正，得到精确的三维粒子速度场。
7.根据权利要求1所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：所述滤波片为线性可变带通滤波器。
8.根据权利要求1所述的一种低粒子密度的单相机粒子图像测速方法，其特征在于：粒子溶液和流体溶液体积比小于100。
CN202110631589.5A 2021-06-07 2021-06-07 一种低粒子密度的单相机粒子图像测速方法 Pending CN113466489A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110631589.5A CN113466489A (zh) 2021-06-07 2021-06-07 一种低粒子密度的单相机粒子图像测速方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110631589.5A CN113466489A (zh) 2021-06-07 2021-06-07 一种低粒子密度的单相机粒子图像测速方法

Publications (1)

Publication Number Publication Date
CN113466489A true CN113466489A (zh) 2021-10-01

Family

ID=77868586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110631589.5A Pending CN113466489A (zh) 2021-06-07 2021-06-07 一种低粒子密度的单相机粒子图像测速方法

Country Status (1)

Country Link
CN (1) CN113466489A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113624446A (zh) * 2021-10-11 2021-11-09 水利部交通运输部国家能源局南京水利科学研究院 一种非接触式分层流场测量方法
CN116449049A (zh) * 2023-03-29 2023-07-18 南京航空航天大学 一种基于多色光深度编码和光场相机的三维流场测试方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866639A (en) * 1987-03-09 1989-09-12 University Of Illinois Method and apparatus for determining the direction of motion in multiple exposure velocimetry
CN103308715A (zh) * 2013-06-20 2013-09-18 南京昊控软件技术有限公司 基于双色激光扫描技术的低功耗粒子图像测速系统
CN108024723A (zh) * 2015-09-15 2018-05-11 三星电子株式会社 用于监测微血管血流动力学的移动光学设备和方法
CN108020680A (zh) * 2016-11-03 2018-05-11 崔胡晋 基于piv的流体测量装置及方法
EP3367106A1 (en) * 2017-02-28 2018-08-29 King Abdullah University Of Science And Technology Rainbow particle imaging velocimetry for dense 3d fluid velocity imaging

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4866639A (en) * 1987-03-09 1989-09-12 University Of Illinois Method and apparatus for determining the direction of motion in multiple exposure velocimetry
CN103308715A (zh) * 2013-06-20 2013-09-18 南京昊控软件技术有限公司 基于双色激光扫描技术的低功耗粒子图像测速系统
CN108024723A (zh) * 2015-09-15 2018-05-11 三星电子株式会社 用于监测微血管血流动力学的移动光学设备和方法
CN108020680A (zh) * 2016-11-03 2018-05-11 崔胡晋 基于piv的流体测量装置及方法
EP3367106A1 (en) * 2017-02-28 2018-08-29 King Abdullah University Of Science And Technology Rainbow particle imaging velocimetry for dense 3d fluid velocity imaging

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIONG J等: "Rainbow particle imaging velocimetry for dense 3D fluid velocity imaging", ACM TRANSACTIONS ON GRAPHICS, vol. 36, no. 4, pages 1 - 14, XP058372810, DOI: 10.1145/3072959.3073662 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113624446A (zh) * 2021-10-11 2021-11-09 水利部交通运输部国家能源局南京水利科学研究院 一种非接触式分层流场测量方法
CN113624446B (zh) * 2021-10-11 2022-01-07 水利部交通运输部国家能源局南京水利科学研究院 一种非接触式分层流场测量方法
CN116449049A (zh) * 2023-03-29 2023-07-18 南京航空航天大学 一种基于多色光深度编码和光场相机的三维流场测试方法及系统

Similar Documents

Publication Publication Date Title
EP3367106B1 (en) Rainbow particle imaging velocimetry for dense 3d fluid velocity imaging
US5905568A (en) Stereo imaging velocimetry
US8922636B1 (en) Synthetic aperture imaging for fluid flows
Fahringer et al. Tomographic reconstruction of a 3-D flow field using a plenoptic camera
CN106970024B (zh) 基于相机和可控频闪光源的限界检测测距方法及系统
CN113466489A (zh) 一种低粒子密度的单相机粒子图像测速方法
Skinner et al. Towards real-time underwater 3D reconstruction with plenoptic cameras
Skinner et al. Automatic color correction for 3D reconstruction of underwater scenes
CN107655571B (zh) 一种基于色散模糊的光谱成像系统及其光谱重建方法
CN114187263B (zh) 一种融合先验引导和域适应的磨损表面朗伯反射分离方法
Douxchamps et al. Stereoscopic and velocimetric reconstructions of the free surface topography of antidune flows
Fahringer et al. The effect of grid resolution on the accuracy of tomographic reconstruction using a plenoptic camera
Shi et al. Development and Application of Light-Field Cameras in Fluid Measurements
CN113030510A (zh) 一种基于三色掩膜单彩色相机的三维流场测试方法及系统
CN110851965A (zh) 一种基于物理模型的光源优化方法及优化系统
RU2734070C9 (ru) Способ измерения пространственного расстояния между малоразмерными объектами
Economikos et al. Toward full-field measurements of instantaneous visualizations of coherent structures in turbulent shear flows
Kurada et al. Particle-imaging techniques for quantitative flow visualization: a review
Leat et al. Analysis of uncertainty in laser-scanned bathymetric maps
He et al. Inexpensive multi-plane particle image velocimetry based on defocusing: Proof of concept on two-component measurement
RU2760845C1 (ru) Способ обнаружения и определения характеристик целей на основе регистрации и обработки хода лучей от объектов в наблюдаемом пространстве и устройство для его реализации
Fomin et al. Correlation analysis of digital images of flows with subpixel accuracy
CN110879300A (zh) 一种流动粒子速度测量方法及系统
US20240126952A1 (en) Arithmetic operation system, training method, and non-transitory computer readable medium storing training program
US20240125935A1 (en) Arithmetic operation system, training method, and non-transitory computer readable medium storing training program

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination