《粒子图像变形单元匹配测速方法.pdf》由会员分享，可在线阅读，更多相关《粒子图像变形单元匹配测速方法.pdf（8页完整版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 104166998 A (43)申请公布日 2014.11.26 CN 104166998 A (21)申请号 201410402082.2 (22)申请日 2014.08.15 G06T 7/00(2006.01) (71)申请人 中国人民解放军国防科学技术大学 地址 410073 湖南省长沙市砚瓦池正街 47 号 (72)发明人 张小虎 谢良 丁少闻 柴进 (74)专利代理机构 湖南省国防科技工业局专利 中心 43102 代理人 冯青 (54) 发明名称 粒子图像变形单元匹配测速方法 (57) 摘要 本发明涉及一种粒子图像变形单元匹配测速 方法。 首先对粒子图像进。
2、行变形单元分割， 在此基 础上建立相应单元的浓度和形状联合变化模型， 最后在全局流场一致性的约束下进行全局的变形 单元匹配， 完成粒子图像测速。 本发明采用了变形 单元的全局一致性匹配测速， 并采用了亚像素级 图像特征， 可实现测量精度可达到亚像素级。 (51)Int.Cl. 权利要求书 2 页 说明书 3 页 附图 2 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书3页 附图2页 (10)申请公布号 CN 104166998 A CN 104166998 A 1/2 页 2 1. 粒子图像变形单元匹配测速方法， 其特征在于 ： 提出了变形单元的浓度和几。
3、何形状 联合变形模型 ； 并在全局约束下进行变形单元的匹配， 进而对多个相邻单元的整体匹配性 进行检查， 对不符合整体匹配性的单元重新进行匹配， 直到多个相邻单元均达到整体匹配， 从而完成全图匹配测速 ； 具体步骤为 ： 步骤 1、 确定变形单元 ； 步骤 2、 建立单元浓度和形状联合变化模型 ； 步骤 3、 在全局约束下进行变形单元的匹配 ； 步骤 4、 单个单元的匹配测速 步骤 5、 全图匹配测速。 2. 根据权利要求 1 所述的粒子图像变形单元匹配测速方法， 其特征在于 ： 所述步骤 1 具体过程为 ： 利 用 SLIC 超 像 素 聚 类 方 法 对 图 像 进 行 单 元 分 割，得。
4、 到 各 个 单 元 。 3. 根据权利要求 1 所述的粒子图像变形单元匹配测速方法， 其特征在于 ： 所述步骤 2 具体过程为 ： 假设单元内像素运动有平移、 旋转和扩散三种形式， 每个单元的像素集为 ， 其中表示位于位置处的像素值 ； 定义每个像 素的运动速度为， 假设单元内的像素扩散形式服从相同的扩散矩阵，大小为 ， 表示位置影响的区域为以为中心， 为半边长 的矩形， 如果已知各像素点的运动参数， 则各像素点向后一帧中像素点的运动形式为 其中表示前一帧，表示后一帧， 则后一帧各个像素值是由前 一帧中所有影响的综合 ： 其中表示能够影响到的像素集 ；表示枚举像素集 ； 表示包含位置像素集的。
5、扩散矩阵在位置的系数 ；表示位置 在位置的扩散矩阵中的系数位置，表示位于位置处的像素值。 4. 根据权利要求 1 所述的粒子图像变形单元匹配测速方法， 其特征在于 ： 所述步骤 3 具体过程为 ： 首先， 对单元施加一个平移运动量和扩散运动量， 即计算出该单元在运动后其几何位 权 利 要 求 书 CN 104166998 A 2 2/2 页 3 置和几何形状， 同时根据联合变形模型， 也计算出在这些运动量作用下该单元变化后的浓 度 ； 接下来， 在运动后图像上取与此变形单元对应区域， 计算该区域浓度分布与变形单元 浓度分布的相似性。 5. 根据权利要求 1 所述的粒子图像变形单元匹配测速方法，。
6、 其特征在于 ： 所述步骤 4 具体过程为 ： 在单元可能发生的平移运动量和扩散运动量范围内， 按自适应步长对单元施加的平移 运动量和扩散运动量进行变化， 重复步骤 2、 3， 最后取与变形单元相似性最大的运动后图像 区域所对应的变形参数作为单元的运动参数， 即完成单个单元的匹配测速。 6. 根据权利要求 1 所述的粒子图像变形单元匹配测速方法， 其特征在于 ： 所述步骤 5 具体过程为 ： 在步骤 4 基础上， 对多个相邻单元的整体匹配性进行检查， 对不符合整体匹配性的单 元重新进行匹配， 直到多个相邻单元均达到整体匹配， 至此， 完成全图匹配测速。 权 利 要 求 书 CN 1041669。
7、98 A 3 1/3 页 4 粒子图像变形单元匹配测速方法 技术领域 0001 本发明涉及风洞实验中快速流动粒子的流场分析、 粒子特性研究， 高超声速 / 亚 声速气流环境下液体速度、 浓度确定， 图像处理中关于粒子特性的研究等领域。 背景技术 0002 传统粒子图像匹配测速方法常采用固定单元相关匹配方法， 其基本假设是单元在 运动过程中浓度和形状保持不变， 由互相关算法匹配解算粒子运动， 进而得出速度场。 而实 际上， 由于单元内的粒子在运动过程中， 除了整体平移， 还存在扩散运动， 因此每个单元在 运动过程中其浓度和形状会发生改变。由于单元的这种变形本质， 固定单元相关匹配方法 在实际应用。
8、中存在很大的局限性。 发明内容 0003 本方法针对常见的变形粒子团测速及确定相对浓度问题， 提出变形单元的新概 念。在流场中， 单元在整体平移的基础上， 由于其内部的扩散运动， 造成其浓度和形状的改 变， 本方法所指的 “变形” 包括了单元的浓度和几何形状的改变。 0004 由于粒子图像是由大量粒子聚集成团状后， 对光的透射亮度的反映。 单元的变形， 本质上是其内部粒子之间相对运动， 造成浓度的改变， 同时改变其几何形状。 为此本方法提 出了变形单元的浓度和几何形状联合变形模型。 0005 在此联合变形模型基础上， 在全局约束下， 进一步提出了变形单元匹配方法， 使用 该方法完成单个单元匹配。
9、测速。 在此基础上， 进而对多个相邻单元的整体匹配性进行检查， 对不符合整体匹配性的单元重新进行匹配， 直到多个相邻单元均达到整体匹配。 至此， 完成 全图匹配测速。具体内容为 ： 1、 确定变形单元 利用 SLIC 超像素聚类方法 （R. Achanta, A. Shaji, K. Smith, A. Lucchi, P. Fua, and S. Susstrunk, SLIC Superpixels, EPFL Technical Report 149300, June 2010）对 图 像 进 行 单 元 分 割，得 到 各 个 单 元 。 0006 2、 建立单元浓度和形状联合变化模型。
10、 假设单元内像素运动有平移、 旋转和扩散三种形式， 每个单元的像素集为 ， 其中表示位于位置处的像素值 ； 定义每个像 素的运动速度为， 假设单元内的像素扩散形式服从相同的扩散矩阵，大小为 ， 表示位置影响的区域为以为中心， 为半边长 的矩形。如果已知各像素点的运动参数， 则各像素点向后一帧中像素点的运动形式为 说 明 书 CN 104166998 A 4 2/3 页 5 其中表示前一帧，表示后一帧。则后一帧各个像素值是由前 一帧中所有影响的综合 ： 其中表示能够影响到的像素集 ；表示枚举像素集 ； 表示包含位置像素集的扩散矩阵在位置的系数 ；表示位置 在位置的扩散矩阵中的系数位置，表示位于位。
11、置处的像素值。 0007 3、 在全局约束下进行变形单元的匹配 首先， 对单元施加一个平移运动量和扩散运动量， 即可计算出该单元在运动后其几何 位置和几何形状， 同时根据联合变形模型， 也可计算出在这些运动量作用下该单元变化后 的浓度。 0008 接下来， 在运动后图像上取与此变形单元对应区域， 计算该区域浓度分布与变形 单元浓度分布的相似性。 0009 4、 单个单元的匹配测速 在单元可能发生的平移运动量和扩散运动量范围内， 按自适应步长对单元施加的平移 运动量和扩散运动量进行变化， 重复步骤 2、 3， 最后取与变形单元相似性最大的运动后图像 区域所对应的变形参数作为单元的运动参数， 即完。
12、成单个单元的匹配测速。 0010 5、 全图匹配测速 在步骤 4 基础上， 对多个相邻单元的整体匹配性进行检查， 对不符合整体匹配性的单 元重新进行匹配， 直到多个相邻单元均达到整体匹配。至此， 完成全图匹配测速。 0011 本方法的优点 ： 单元可变形， 符合物理现象。联合模型， 几何变化与浓度变化相关联。单元内每个点均 有速度参数， 而传统方法单元内速度一致。 采用了变形单元的全局一致性匹配测速， 并采用 了亚像素级图像特征， 可实现测量精度可达到亚像素级。 附图说明 0012 图 1 是变形单元在相邻帧的结构变化示意图， 图 2 是像素点平移、 旋转和扩散示意图。 具体实施方式 0013。
13、 1、 确定变形单元 利用SLIC超像素聚类方法 （R. Achanta, A. Shaji, K. Smith, A. Lucchi, P. Fua, and S. Susstrunk, SLIC Superpixels, EPFL Technical Report 149300, June 2010） 对 图像进行单元分割， 得到各个单元。 说 明 书 CN 104166998 A 5 3/3 页 6 0014 2、 建立单元浓度和形状联合变化模型 假设单元内像素运动有平移、 旋转和扩散三种形式， 每个单元的像素集为 ， 其中表示位于位置处的像素值 ； 定义每个像 素的运动速度为， 假设单。
14、元内的像素扩散形式服从相同的扩散矩阵，大小为 ， 表示位置影响的区域为以为中心， 为半边长 的矩形。如果已知各像素点的运动参数， 则各像素点向后一帧中像素点的运动形式为 其中表示前一帧，表示后一帧。则后一帧各个像素值是由前 一帧中所有影响的综合 ： 其中表示能够影响到的像素集 ；表示枚举像素集 ； 表示包含位置像素集的扩散矩阵在位置的系数 ；表示位置 在位置的扩散矩阵中的系数位置，表示位于位置处的像素值。 0015 3、 在全局约束下进行变形单元的匹配 首先， 对单元施加一个平移运动量和扩散运动量， 即可计算出该单元在运动后其几何 位置和几何形状， 同时根据联合变形模型， 也可计算出在这些运动。
15、量作用下该单元变化后 的浓度。 0016 接下来， 在运动后图像上取与此变形单元对应区域， 计算该区域浓度分布与变形 单元浓度分布的相似性。 0017 在单元可能发生的平移运动量和扩散运动量范围内， 按一定步长对单元施加的平 移运动量和扩散运动量进行变化， 重复上述两个步骤， 最后取与变形单元相似性最大的运 动后图像区域所对应的变形参数作为单元的运动参数， 即完成单个单元的匹配测速。 0018 在此基础上， 对多个相邻单元的整体匹配性进行检查， 对不符合整体匹配性的单 元重新进行匹配， 直到多个相邻单元均达到整体匹配。至此， 完成全图匹配测速。 说 明 书 CN 104166998 A 6 1/2 页 7 图 1 说 明 书 附 图 CN 104166998 A 7 2/2 页 8 图 2 说 明 书 附 图 CN 104166998 A 8 。