基于维诺图的粒子图像测速方法及装置制造方法及图纸

本发明专利技术涉及一种基于维诺图的粒子图像测速方法及装置，该方法采用白光LED作为光源照射在待测颗粒流表面，利用高速CCD相机作为图像探测设备对颗粒流表面的粒子进行连续拍摄，然后用计算机采集相机拍摄的图像，并对采集到的图像进行图像预处理，去除由于非均匀光照等实验环境和设备造成的背景噪声以及目标筛选，将处理后的图像内颗粒目标提取为点状粒子，并用其质心坐标来代替，由此得到一组二维平面的特征点坐标，对其进行维诺图构建，接着用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息。

Particle image velocimetry method and device based on Wiener diagram

The present invention relates to a particle image velocimetry method and device based on Venogram. The method uses white LED as light source to irradiate the surface of particle flow to be measured, uses high-speed CCD camera as image detection equipment to continuously photograph particles on the surface of particle flow, and then uses a computer to collect images taken by the camera and collects them. The resulting image is preprocessed to remove the background noise and target selection caused by non-uniform illumination and other experimental environment and equipment. The particle target is extracted as a point particle and replaced by its centroid coordinates. A set of two-dimensional plane feature point coordinates is obtained and its Veno diagram is constructed. Then, the matching between the two frames is performed by using the constructed Veno diagram, and the particle velocity distribution information is obtained.

【技术实现步骤摘要】
基于维诺图的粒子图像测速方法及装置
本专利技术涉及一种粒子速度测量方法和装置，具体涉及一种基于维诺图的粒子速度场分布测量方法和装置。
技术介绍
颗粒流是颗粒物质存在的一种状态，颗粒物质广泛存在于自然界和工业中，如煤炭的运输和储存，泥石流和雪崩现象等都涉及颗粒流问题，因此研究颗粒流的运动特征具有重要意义。在粒子运动速度测量领域，用图像处理的方法进行速度测量主要有粒子图像测速法(PIV)和粒子跟踪测速法(PTV)。PIV和PTV都是上世纪70年代末80年代初发展起来的一种瞬态、多点、无接触式的流体力学测速方法。作为PIV和PTV技术的核心，流场图像分析目前主要采用二维快速Fourier变换实现互相关函数的计算，重庆大学李正浩等人利用这种方法来测量表面流场物体速度，(201710242846.X[一种基于粒子图像测速的表面流场测量方法])，但是这类方法对实验参数较为敏感，且在计算过程中容易产生错误的匹配信息，导致产生速度伪矢量，匹配精度不高。维诺图是计算机几何中的一个重要研究内容，它在气象学，地理学，计算机图形学，机器人学等领域得到广泛应用，而关于维诺图的研究更多的集中于其高效的生成算法以及应用层面，如沈阳建筑大学宋晓宇等人专利技术了一种利用维诺图的最近邻查询方法(201210430002.5[一种基于维诺图的反K最近邻查询方法])。但至目前为止，在流体学中，尤其是针对流场内运动粒子速度方面的研究还尤为不足。因此利用维诺图来进行粒子匹配进而获得运动粒子速度信息，是一种具有巨大潜力和发展前景的一种粒子速度测量技术。
技术实现思路
本专利技术是针对目前流体学中PIV和PTV测量法的不足之处(对实验参数敏感，容易产生误匹配，匹配精度不高)，提出的一种基于维诺图的粒子图像测速方法及装置，来实现对颗粒流速度场分布的准确测量。与传统的点与点之间的匹配方法不同，基于维诺图的匹配方法，通过构建维诺图的泰森多边形，实现了面与面之间的匹配，因此，精度更高、稳定性更好、抗干扰能力更强。为实现上述目的，本专利技术的技术方案是：一种基于维诺图的粒子图像测速方法，采用白光LED作为光源照射在待测颗粒流表面，利用高速CCD相机作为图像探测设备对颗粒流表面的粒子进行连续拍摄，然后用计算机采集相机拍摄的图像，并对采集到的图像进行图像预处理，去除由于非均匀光照等实验环境和设备造成的背景噪声以及目标筛选，将处理后的图像内颗粒目标提取为点状粒子，并用其质心坐标来代替，由此得到一组二维平面的特征点坐标，对其进行维诺图构建，接着用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息。所述维诺图构建，用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息的具体步骤为：步骤一：假设有一组相邻的两帧图像，其粒子点的集合为A＝{a1，a2，...，an}和B＝{b1，b2，...，bn}，其中ai，bj分别为对应集合内某一个确定的点，然后分别为A，B建立维诺多边形图层P＝{p1，p2，...，pn}，其中pi为P中维诺多边形图层中某一个确定的泰森多边形，Q＝{q1，q2，...，qn}，其中qj为Q中维诺多边形图层中某一个确定的泰森多边形；步骤二：建立匹配条件：在相邻两帧之间互相筛选颗粒，首先遍历P中的每个要素：pi为ai对应的泰森多边形，找出仅包含一个B中要素：设为bj的所有多边形，计算这些多边形中ai与bj的距离值，借鉴拉依达准则对其改进，将小于μ-δ其中μ为平均距离，δ为距离标准差，以及大于μ+δ的值作为异常值剔除，计算出剩余点对的距离均值μ和距离标准差δ，将匹配点对距离小于等于μ+3δ作为选择候选集的条件之一；步骤三：进行匹配计算：查看包含在pi中的B的要素数目1)若包含在pi中的B中的要素数目为0，则看pi是否存在与其相交且Q中的泰森多边形，pi与qj要满足一定的相似度，即要求pi与qj的相交面积与pi的面积之比大于20％，若不存在符合条件qj，则表示pi无匹配对象，若仅存一个与pi相交的泰森多边形，则计算两个面中对应点之间的距离是否符合条件一，若符合则将对应点视为匹配对象，若存在多个与pi相交的泰森多边形，先剔除不符合条件一的匹配对象，然后比较符合条件一的对应对应泰森多边形的形状相似度，将相似度最高的qj视为匹配对象，并建立对应点的匹配关系；2)若包含在pi中的B中的要素数目为1，则检查其是否满足条件一，若满足则建立对应点的匹配关系，否则记录对应的bj为新增要素；3)若包含在pi中的B中的要素数目大于等于2，先剔除不符合条件一的匹配对象，然后比较符合条件一的对应泰森多边形的形状相似性，将相似度最高的qj视为匹配对象，并建立对应点的匹配关系；步骤四：若经过上述匹配计算后点与点之间的匹配关系仍出现一对多的情况，可以考虑加入反算方法以及概率择优模块，即将原本由第一帧向第二帧的顺匹配反转为由第二帧向第一帧的逆匹配，做第二次匹配，重复上述步骤后得到新的匹配关系，记为F，F蕴含着可能性的积累，可以利用概率择优模块保留出现频率最高的模块，同时剔除一切分歧性的结果，建立对应点的匹配关系；步骤五：观察已经建立好的匹配点分别相对于窗口的坐标，可以得到点与点之间的相对位移，记为同时通过设置相机的帧率可以获得两帧图像之间的时间间隔，记为Δt，由此通过某一个点的速度Vi计算公式：可以得到某一个粒子在两帧图像之间运动的平均速度，因为图像之间时间间隔很小，故可以用此速度来代替粒子在某一时刻的瞬时速度。所述泰森多边形相似度计算方法：其中ρArea(ai，bj)表示面积相似度，从一定程度上可以反映多边形位置和形状相似性，计算方法见公式1；T为面积相似度阈值，S为两个多边形相交部分的面积大小，min为取两个多边形面积最小的那个，当两泰森多边形面积重叠度较高时，建议采用面积相似度作为评判对应泰森多边形的依据；ρShape(ai，bj)表示对应泰森多边形的轮廓形状相似度，其计算方法借鉴于折线段方位编码方法，即将多边形的轮廓形状通过各边的方位编码来描述；当边的方位角恰好位于圆周分割线上时，需要将计算出的编码降为该编码的下一级，为了便于待匹配多边形之间各边方位的比较，计算方位编码前需要将两个待匹配的多变性上的节点统一按顺时针或者逆时针进行排列，每个边的方位编码可以采用公式(3)：其中Chr为将数值转换为编码字符(A～X)的函数，Nl、Nl+1为折线上的两个相邻接点，(NlNl+1)表示其所构成的向量的方位角，int为取整函数。一种基于维诺图的粒子图像测速装置，包括一台计算机，一白色LED光源，一带镜头的高速CCD相机，一可以进行粒子速度场研究的二维圆盘，二维圆盘内放入体积总量为圆盘容积40％的2mm-5mm均匀的球形陶瓷颗粒，并垂直放置于由电动机驱动的匀速转台上，白色LED光源均匀照射在二维圆盘内球形陶瓷颗粒表面，带有镜头的高速CCD相机固定于支架上，并与计算机相连，用来处理带有镜头的高速CCD采集到的图片。本专利技术的有益效果是：针对于现有的速度测量方法的不足之处，本专利技术提出的方法较传统方法更有优势：对实验参数的选取不敏感；算法的复杂度不高；很大程度上提高了匹配精确度，经实验论证可达92％-94％的匹配准确性；对实验环境的要求不高，测量系统相对简单。附图说明图1为泰森多边形形状本文档来自技高网...

【技术保护点】
1.一种基于维诺图的粒子图像测速方法，其特征在于：采用白光LED作为光源照射在待测颗粒流表面，利用高速CCD相机作为图像探测设备对颗粒流表面的粒子进行连续拍摄，然后用计算机采集相机拍摄的图像，并对采集到的图像进行图像预处理，去除由于非均匀光照等实验环境和设备造成的背景噪声以及目标筛选，将处理后的图像内颗粒目标提取为点状粒子，并用其质心坐标来代替，由此得到一组二维平面的特征点坐标，对其进行维诺图构建，接着用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息。

【技术特征摘要】
1.一种基于维诺图的粒子图像测速方法，其特征在于：采用白光LED作为光源照射在待测颗粒流表面，利用高速CCD相机作为图像探测设备对颗粒流表面的粒子进行连续拍摄，然后用计算机采集相机拍摄的图像，并对采集到的图像进行图像预处理，去除由于非均匀光照等实验环境和设备造成的背景噪声以及目标筛选，将处理后的图像内颗粒目标提取为点状粒子，并用其质心坐标来代替，由此得到一组二维平面的特征点坐标，对其进行维诺图构建，接着用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息。2.根据权利要求1所述的基于维诺图的粒子图像测速方法，其特征在于：所述维诺图构建，用构建好的维诺图进行两帧图像之间的匹配，最后得到粒子速度分布信息的具体步骤为：步骤一：假设有一组相邻的两帧图像，其粒子点的集合为A＝{a1，a2，...，an}和B＝{b1，b2，...，bn}，其中ai，bj分别为对应集合内某一个确定的点，然后分别为A，B建立维诺多边形图层P＝{p1，p2，...，pn}，其中pi为P中维诺多边形图层中某一个确定的泰森多边形，Q＝{q1，q2，...，qn}，其中qj为Q中维诺多边形图层中某一个确定的泰森多边形；步骤二：建立匹配条件：在相邻两帧之间互相筛选颗粒，首先遍历P中的每个要素：pi为ai对应的泰森多边形，找出仅包含一个B中要素：设为bj的所有多边形，计算这些多边形中ai与bj的距离值，借鉴拉依达准则对其改进，将小于μ-δ其中μ为平均距离，δ为距离标准差，以及大于μ+δ的值作为异常值剔除，计算出剩余点对的距离均值μ和距离标准差δ，将匹配点对距离小于等于μ+3δ作为选择候选集的条件之一；步骤三：进行匹配计算：查看包含在pi中的B的要素数目1)若包含在pi中的B中的要素数目为0，则看pi是否存在与其相交且Q中的泰森多边形，pi与qj要满足一定的相似度，即要求pi与qj的相交面积与pi的面积之比大于20％，若不存在符合条件qj，则表示pi无匹配对象，若仅存一个与pi相交的泰森多边形，则计算两个面中对应点之间的距离是否符合条件一，若符合则将对应点视为匹配对象，若存在多个与pi相交的泰森多边形，先剔除不符合条件一的匹配对象，然后比较符合条件一的对应对应泰森多边形的形状相似度，将相似度最高的qj视为匹配对象，并建立对应点的匹配关系；2)若包含在pi中的B中的要素数目为1，则检查其是否满足条件一，若满足则建立对应点的匹配关系，否则记录对应的bj为新增要素；3)若包含在pi中的B中的要素数目大于等于2，先剔除...

【专利技术属性】
技术研发人员：杨晖，宋健，李然，
申请(专利权)人：上海理工大学，
类型：发明
国别省市：上海,31