CN109633199B - 一种基于相关匹配的粒子聚焦两点像测速方法 - Google Patents
一种基于相关匹配的粒子聚焦两点像测速方法 Download PDFInfo
- Publication number
- CN109633199B CN109633199B CN201811510172.8A CN201811510172A CN109633199B CN 109633199 B CN109633199 B CN 109633199B CN 201811510172 A CN201811510172 A CN 201811510172A CN 109633199 B CN109633199 B CN 109633199B
- Authority
- CN
- China
- Prior art keywords
- particle
- matching
- particles
- correlation
- candidate
- 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.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 177
- 238000000691 measurement method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000006073 displacement reaction Methods 0.000 claims abstract description 11
- 238000003384 imaging method Methods 0.000 claims description 8
- 239000012467 final product Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 16
- 238000005516 engineering process Methods 0.000 description 7
- IMNFDUFMRHMDMM-UHFFFAOYSA-N N-Heptane Chemical compound CCCCCCC IMNFDUFMRHMDMM-UHFFFAOYSA-N 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 6
- 239000007921 spray Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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/22—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 auto-correlation or cross-correlation detection means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N2015/0003—Determining electric mobility, velocity profile, average speed or velocity of a plurality of particles
-
- G01N2015/1027—
Abstract
本发明公开了一种基于相关匹配的粒子聚焦两点像测速方法,该方法是记录时间序列的粒子IPI聚焦两点像,基于质心法和序列模板匹配法,实现亚像素的粒子定位,根据最大位移和速度方向准则搜索候选粒子,再基于互相关方法,实现高准确识别率的粒子识别匹配,其效果为基于IPI聚焦像的粒子速度测量提供一种方法,且具有很高的速度测量精度和很高匹配准确率,可用于喷雾、气泡、流场等粒子场速度分布测量。
Description
技术领域
本发明涉及粒子参数测量领域,尤其涉及一种IPI(干涉粒子成像)的粒子聚焦两点像跟踪匹配测速算法,特别涉及粒子两点像定位精度及其匹配准确率的提高。
背景技术
粒子速度是表征粒子特征的一个重要参数,研究粒子的速度测量方法在喷雾场、流场中具有非常重要的意义。常用粒子测速的光学方法有激光多普勒测速技术(LDV)、相位多普勒测速技术(PDA)、粒子图像测速技术(PIV)和粒子跟踪测速技术(PTV)等。LDV和PDA是单点测速技术,不能给出全场的瞬时流速数据。PIV和PTV技术突破了空间单点测量的局限性,实现了全流场瞬态测量,且PTV技术通过对流场中每个粒子的跟踪匹配,可以给出单个粒子的运动速度和轨迹。在PTV中,粒子跟踪匹配,即将同一粒子在连续时间序列的位置匹配,是关键。
在PTV技术中,两帧粒子跟踪匹配算法是常用的一种粒子跟踪匹配方法,其基本思想是根据粒子场运动特征,基于一定的限定,从连续两帧粒子像中跟踪匹配粒子对。基于该思想,已提出了多种粒子跟踪匹配方法,如专利CN103558409A公开一种水流中分布式PTV流场测量系统及其测量方法,利用粒子运动距离最短的方法对连续两帧图像中的粒子进行匹配,计算得到表面流场。Back和Lee提出一种基于匹配概率的两帧粒子跟踪算法,该方法是基于最大速度、小的速度方向变化、运动一致性和匹配一致性四个准则对两帧图像中的粒子进行匹配,通过匹配概率的迭代评估,在第二帧图像中选择与第一帧图像中匹配概率最大的粒子作为匹配粒子,实现粒子场速度测量(Exp.Fluids 1986;22(1):23-32)。
IPI技术的基本原理是基于Mie散射成像技术,粒子散射光在离焦面上形成条纹像,在聚焦面上形成两点像,通过测量干涉条纹图的条纹频率或聚焦两点像之间的间距得到粒子尺寸大小,结合PTV或PIV技术实现粒子速度测量。G.Pan等基于两帧粒子跟踪匹配算法,结合粒子尺寸信息识别匹配粒子IPI条纹图(Proc ILASS 2005,91-96)。相比于IPI圆形条纹图,聚焦两点像由于降低了重叠,特别适应于高密度粒子场测量,然而其定位和识别更加困难,特别是对高密度粒子场。目前还没有提出适用于IPI聚焦两点像跟踪匹配方法。
发明内容
本发明提供一种基于相关匹配的粒子聚焦两点像测速方法,该方法是记录时间序列的粒子IPI聚焦两点像,基于序列模板匹配和质心法定位两点像的中心,得到亚像素精度的粒子位置坐标,根据最大位移和速度方向准则搜索候选粒子,再基于互相关方法,实现高准确识别率的粒子识别匹配,详见下文描述:
一种基于相关匹配的粒子聚焦两点像测速方法,所述测速方法包括以下步骤:
1)定位两点像,利用干涉粒子成像系统,记录时间序列的粒子聚焦两点像P(x,y;ti),i=1,2,...,将聚焦两点像与粒子掩模像进行相关运算,其峰值位置即为两点像的位置坐标;
2)亚像素定位,对两点像的坐标,结合质心法,得到亚像素的位置坐标(xi,yi;ti);
3)粒子跟踪匹配,设第一帧图像和第二帧图像中的粒子两点像分别为和粒子像的位置坐标分别为(x1,y1;t1)和(x2,y2;t2),j为粒子序列数,j=1表示第一个粒子,j=2表示第二个粒子,依次类推。将第一帧中任一粒子作为目标粒子,对目标粒子进行跟踪匹配;
所述对目标粒子进行跟踪匹配的步骤具体为:
第二步:匹配粒子对,对第一步确定的候选粒子,若只有一个候选粒子,则该候选粒子即为匹配粒子;若有多个候选粒子时,计算目标粒子与候选粒子Pc 2(x,y;t2)相关值,若最大相关峰值数目为1,该粒子即为的匹配粒子;若最大相关峰值数目大于1,根据粒子场中粒子运动一致性准则,选择最大概率运动方向上的粒子作为匹配粒子。
其中,当互相关的峰值分布平坦时,即满足|C0-C-1|<ε和/或|C0-C+1|<ε条件时,粒子两点像的位置坐标(x,y)由下式给出:
其中,C0、C-1和C+1分别是离散互相关的峰值以及其邻近坐标的相关值,对应的位置坐标分别记为(x0,y0)、(x-1,y0)/(x0,y-1)和(x+1,y0)/(x0,y+1),ε为预设值,为一小量。
本发明提供的技术方案的有益效果是:
1、本发明结合了质心法,实现了亚像素的粒子定位,具有很高的粒子速度测量精度;
2、本方法是利用粒子两点像间距Δl(Δl∝M×d),其实质是利用粒子尺寸d。若成像系统放大率M相同,则是利用粒子尺寸d,若成像系统放大率不同,则是利用成像系统放大率与粒子尺寸的乘积,即d×M,这时根据M,确定匹配粒子对,其实质仍是利用粒子尺寸d的信息。
但本发明这种基于Δl的方法优于文献(Proc ILASS 2005,91-96)的方法,是因为避免了计算粒径d引入的误差,且相比于比较粒径方法,相关方法更敏感,因此具有更准确的匹配识别率。
3、该方法尤其适用于相同放大率M的IPI聚焦成像系统。
附图说明
图1为本发明提出的一种基于相关匹配的粒子聚焦两点像测速方法的流程图;
图2为采集的正庚烷喷雾场双曝光聚焦两点像;
图4为本发明的方法与传统的两帧粒子跟踪匹配方法的比较。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
实施例1
本发明实施例给出了一个基于IPI聚焦像的速度测量实验结果,实验所用激光器为波长λ=532nm的Nd:YAG双脉冲激光器,脉冲发射时间间隔Δt=0.8ms,最大脉冲激光能量为120mJ。激光器发出的细光束经扩束准直系统、再经两个柱面镜压缩为光束宽度为0.9mm的片状光束,通过分束镜将其分为强度相等的两束光相向照射粒子场。在散射角90°方向记录粒子场的聚焦像。
图2为燃油温度80℃,环境背压1.0bar,采集的正庚烷喷雾场双曝光聚焦两点像。将聚焦两点像与序列模板进行相关运算,再结合质心法,提取每个液滴粒子两点像的位置坐标,如图3(a)所示,图中“*”表示提取到的粒子位置。根据燃油喷射压力和喷嘴参数预先估计喷雾场粒子的最大运动速度vmax=500mm/s,计算得到的最大运动位移Dmax=1.5mm(435pixels)。将图3(a)中的每个粒子作为目标粒子,根据最大位移Dmax和速度方向准则搜索可能匹配的候选粒子。再基于互相关算法,计算目标粒子与候选粒子相关值,确定匹配粒子,计算得到粒子运动速度。
如将图3(a)中任一P1 1作为目标粒子,根据最大位移Dmax和速度方向准则搜索可能匹配的候选粒子为P1 2,P3 1和P3 2,共5个粒子,如图3(a)标记区域所示。利用自相关和高斯插值,提取的P1 1,P1 2,和间距Δl分别为16.71,16.55,16.08,16.85,16.98和15.45pixels,由式计算得到对应的粒子直径分别为21.7,21.5,20.9,21.9,22.1,20.1μm。
计算目标粒子P1 1与P1 2,P3 1和P3 2的互相关值,如图3(b)所示。P1 1与P1 2相关值最高,P1 2为P1 1的匹配粒子。若根据粒径匹配准则,P1 1直径21.7μm,P1 2直径21.5μm,直径21.9μm。P1 1与P1 2直径相近,所以为匹配粒子对,P1 1与直径也相近,也为匹配粒子对,但实际这是假匹配粒子对。相比于粒径匹配准则,基于互相关的匹配准则具有更准确的匹配率。对P3 1和P3 2以同样的方法处理,与为匹配粒子对,P3 1和P3 2为匹配粒子对。对图3(a)中每个粒子进行跟踪匹配,共匹配出102对粒子,其速度矢量分布如图3(c)所示,箭头代表了粒子运动方向,长短表示运动速度大小。
实施例2
实例2给出模拟分析结果。在512×512pixels的像方视场中,模拟不同数的随机分布的粒子两点像,两点像间距分布在4~20pixels之间,粒子运动位移从0到20pixels,运动方向从-90°到0°。为了减少模拟中的随机性,取50次模拟的平均结果作为最终结果。对模拟的每一幅粒子两点像,基于序列模板进行相关运算,再结合质心法,提取每个粒子两点像的位置坐标,将每个粒子作为目标粒子,根据最大位移Dmax和速度方向准则搜索可能匹配的候选粒子。再基于互相关算法,计算目标粒子与候选粒子相关值,确定匹配粒子。图4给出了本发明的方法与基于匹配概率的两帧粒子跟踪方法的粒子匹配率RM随粒子数密度的变化曲线。本发明的方法明显高于该方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于相关匹配的粒子聚焦两点像测速方法,其特征在于,所述测速方法包括以下步骤:
1)定位两点像,利用干涉粒子成像系统,记录时间序列的粒子聚焦两点像P(x,y;ti),i=1,2,...,将聚焦两点像与粒子掩模像进行相关运算,其峰值位置即为两点像的位置坐标;
2)亚像素定位,对两点像的位置坐标,结合质心法,得到亚像素的位置坐标(xi,yi;ti);
3)粒子跟踪匹配,设第一帧图像和第二帧图像中的粒子两点像分别为和粒子像的位置坐标分别为(x1,y1;t1)和(x2,y2;t2),j为粒子序列数,将第一帧中任一粒子作为目标粒子,对目标粒子进行跟踪匹配;
其中,所述对目标粒子进行跟踪匹配的步骤具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811510172.8A CN109633199B (zh) | 2018-12-11 | 2018-12-11 | 一种基于相关匹配的粒子聚焦两点像测速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811510172.8A CN109633199B (zh) | 2018-12-11 | 2018-12-11 | 一种基于相关匹配的粒子聚焦两点像测速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633199A CN109633199A (zh) | 2019-04-16 |
CN109633199B true CN109633199B (zh) | 2020-09-22 |
Family
ID=66072673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811510172.8A Active CN109633199B (zh) | 2018-12-11 | 2018-12-11 | 一种基于相关匹配的粒子聚焦两点像测速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633199B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110308302B (zh) * | 2019-07-01 | 2020-11-20 | 北京大学 | 一种近壁面流速测量方法及装置 |
CN112036072B (zh) * | 2020-07-10 | 2023-05-12 | 北京航空航天大学 | 一种三维示踪粒子匹配方法及速度场测量装置 |
CN112326192B (zh) * | 2020-10-28 | 2021-07-20 | 西南石油大学 | 一种支撑剂平板输送实验中支撑剂速度定量表征方法 |
CN114062712B (zh) * | 2021-09-29 | 2022-09-06 | 东南大学 | 基于单光场成像的合成孔径粒子图像测速方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674791A (zh) * | 2013-12-16 | 2014-03-26 | 天津大学 | 一种基于双光束照射的干涉粒子成像测量方法 |
CN104166998A (zh) * | 2014-08-15 | 2014-11-26 | 中国人民解放军国防科学技术大学 | 粒子图像变形单元匹配测速方法 |
CN106018201A (zh) * | 2016-05-26 | 2016-10-12 | 天津大学 | 基于均值滤波的混合场粒径测量方法 |
CN106092859A (zh) * | 2016-05-26 | 2016-11-09 | 天津大学 | 基于激光干涉成像和同轴全息的粒子形状判别系统及方法 |
WO2018206474A1 (en) * | 2017-05-12 | 2018-11-15 | Koninklijke Philips N.V. | Method of suppressing false positive signals during self mixing interference particle detection |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8830476B2 (en) * | 2012-03-19 | 2014-09-09 | The United States Of America As Represented By The Secretary Of The Army | Methods and apparatuses for contact-free holographic imaging of aerosol particles |
-
2018
- 2018-12-11 CN CN201811510172.8A patent/CN109633199B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674791A (zh) * | 2013-12-16 | 2014-03-26 | 天津大学 | 一种基于双光束照射的干涉粒子成像测量方法 |
CN104166998A (zh) * | 2014-08-15 | 2014-11-26 | 中国人民解放军国防科学技术大学 | 粒子图像变形单元匹配测速方法 |
CN106018201A (zh) * | 2016-05-26 | 2016-10-12 | 天津大学 | 基于均值滤波的混合场粒径测量方法 |
CN106092859A (zh) * | 2016-05-26 | 2016-11-09 | 天津大学 | 基于激光干涉成像和同轴全息的粒子形状判别系统及方法 |
WO2018206474A1 (en) * | 2017-05-12 | 2018-11-15 | Koninklijke Philips N.V. | Method of suppressing false positive signals during self mixing interference particle detection |
Non-Patent Citations (1)
Title |
---|
粒子场的聚焦两点像测量方法研究;吕通;《中国优秀硕士学位论文全文数据库》;20150515(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109633199A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633199B (zh) | 一种基于相关匹配的粒子聚焦两点像测速方法 | |
Kumar et al. | Automated droplet size distribution measurements using digital inline holography | |
Maeda et al. | Novel interferometric measurement of size and velocity distributions of spherical particles in fluid flows | |
Yang et al. | Digital particle holographic system for measurements of spray field characteristics | |
JP4568800B2 (ja) | 小滴の状態計測装置及び該装置におけるカメラの校正方法 | |
Buchmann et al. | Ultra-high-speed 3D astigmatic particle tracking velocimetry: application to particle-laden supersonic impinging jets | |
Choo et al. | The characteristics of the particle position along an optical axis in particle holography | |
CN105043946A (zh) | 基于双波长的散射角自标定全场彩虹测量方法及装置 | |
CN108362906B (zh) | 一种分子标记二维速度矢量测量装置和方法 | |
Wang et al. | Effect of two typical focus positions on microstructure shape and morphology in femtosecond laser multi-pulse ablation of metals | |
CN100378432C (zh) | 飞秒激光散斑相关法测量微小位移的装置和方法 | |
Mikheev et al. | Enhanced particle-tracking velocimetry (EPTV) with a combined two-component pair-matching algorithm | |
Kieft et al. | The application of a 3D PTV algorithm to a mixed convection flow | |
CN101793530B (zh) | 一种测试光纤陀螺去耦系数的测试方法 | |
Zama et al. | Simultaneous measurement method of size and 3D velocity components of droplets in a spray field illuminated with a thin laser-light sheet | |
Lin et al. | Dual-sheet interferometric particle imaging for opaque particle size and 2D location measurement | |
Pan et al. | Simultaneous global size and velocity measurement of droplets and sprays | |
CN112036072A (zh) | 一种三维示踪粒子匹配方法及速度场测量装置 | |
Angel-Toro et al. | Speckle decorrelation influence on measurements quality in vortex metrology | |
Zama et al. | Simultaneous measurement of droplet size and three-components of velocity in spray | |
Okamoto et al. | Velocity and angle of spatter in fine laser processing | |
JP2003270261A (ja) | 局所空間平均粒子追跡法 | |
Lai et al. | Volumetric three-component velocimetry: a new tool for 3d flow measurement | |
Lu et al. | Visualization of spatial distribution of the droplet size and velocity in flash boiling spray with extended glare-point imaging technique | |
CN210155011U (zh) | 一种基于数字同轴全息测量液滴折射率的装置 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |