资源描述:
基于 IPP 图像软件的管流中粗颗粒运动信息提取方法研究 ① 廖 帅, 曹 斌, 夏建新 (中央民族大学 环境系,北京 100081) 摘 要 运用大功率激光光源和 Nikon D810 单反相机进行了管道水力输送模拟试验,并以试验中得到的管道断面图像为研究对 象,详细介绍了基于图像软件 Image⁃Pro Plus(IPP)的管道断面颗粒提取方法,并以 5 张图像的处理结果对此方法进行验证。 结果 表明,提取出的颗粒数量平均误差小于 5%,说明该种分析方法能准确地将目标颗粒提取并呈现出来,是一种具有代表性的、可操作 的管道断面颗粒图像处理方法。 关键词 管道输送; 粗颗粒; 运动信息; 图像处理; IPP 中图分类号 TD58文献标识码 Adoi10.3969/ j.issn.0253-6099.2017.05.002 文章编号 0253-6099(2017)05-0005-05 Information Extraction Method for Motion of Coarse Particles in Pipeline Based on IPP Image Software LIAO Shuai, CAO Bin, XIA Jian⁃xin (Department of Environment, Minzu University of China, Beijing 100081, China) Abstract Experiments were conducted by simulating coarse particles transported in pipeline, with high⁃powered laser and a Nikon D810 SLR camera used to obtain images of pipeline section for further study. A detailed description on how to extract information about coarse particles at the pipeline section based on Image⁃Pro Plus (IPP) was presented. A further processing of five section images showed such method can obtain and present the accurate information about the targeted coarse particles with an average error less than 5%, verifying it′s a feasible and practical image processing method for extracting the motion of coarse particles at the pipe section. Key words pipeline transportation; coarse particle; movement information; image processing; Image⁃Pro Plus (IPP) 近年来,越来越多的管道输送工程涉及到粗颗粒 的输送,如深海采矿过程中需利用管道将平均粒径为 20 mm 的粗颗粒矿石提升到采矿船中,此外水力采煤、 油井砾石充填等相关工程领域都涉及粗颗粒输送[1]。 粗颗粒在管流中运动状态和垂线分布较复杂,给管道 系统设计和输送参数确定带来困难。 由于理论模型上 不成熟,采用试验方法分析粗颗粒在管流中的运动特 性是目前有效的途径,其中,利用光学图像获取颗粒运 动信息可获得分析的基础数据。 许多学者对管流中颗粒运动特性开展过试验研 究[2-3]。 在 20 世纪 90 年代,路展民[4]等研制的激光 流速仪在测量流体运动特征量上有了明显进步,但这 种仪器只适用于透明的颗粒,且试验中颗粒最大体积 浓度不宜超过 5%。 1998 年,刘青泉[5]利用激光多普 勒分相测速仪(LDV)测量了水平方管的两相浓度、湍 流度以及泥沙浓度等参数,实现了对固液两相流的无 接触测量,并用复合判别法区分了来自粗颗粒和微粒 的多普勒信号。 但以上两种技术都属于单点测量技 术,难以实现对流场的全场、瞬态测量。 近年来,粒子 图像测速技术(particle image velocimetry,PIV)发展迅 速,它充分利用现代计算机技术、光学技术以及图像分 析技术的研究成果,已逐渐成为最新流体测试信息提 取手段[6]。 数字图像处理(Digital in Age Processing)技术是 指将图像信号转化为数字信号,然后利用计算机达到 改善图像视觉效果的处理方式[7-8]。 Image⁃Pro Plus (IPP)软件是美国 Media Cybernetics 公司开发的图像 分析软件包,可完成图像的输入、增强、报告、存档、输 出等功能,被广泛应用于医学、生物学、工业等专业领 域[9-12],IPP 具有误差小、测量参数多等优点,可以满 ①收稿日期 2017-04-15 基金项目 国家自然科学基金(51339008,51209238,51434002);中央民族大学创新实验班创新研究专项计划(SHSY2016120027) 作者简介 廖 帅(1996-),男,湖北黄冈人,主要从事水沙环境工程研究。 第 37 卷第 5 期 2017 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.37 №5 October 2017 ChaoXing 足多个领域图像处理的要求,且可以应对来自数码相 机等不同图像采集设备的图像。 本文运用 Nikon D810 相机、连续型激光光源获取 管流中粗颗粒运动图像,并提出了基于 IPP 软件的管 道断面图像颗粒运动信息的提取方法,用于试验数据 的分析。 1 试验测试系统及图像获取方法 1.1 试验材料、仪器与工况参数 试验中的液相为清水,固相为天然石英砂,颗粒密 度 2 650 kg/ m3。 试验所用图像采集工具为 Nikon D810 单反相机,曝光时间达到 1/8 000 s,图像有效像 素为 3 638 万。 另使用了由 Beam Q 公司生产的连续 型激光光源,激光器功率 10 W,激光波长 532 nm(绿 光),激光光柱宽 3 mm。 试验选取石英砂粒径分别为 0.5、1.0 和 3.0 mm, 设计平均流速分别为 1.5、2.0 和 2.5 m/ s,固相颗粒体 积浓度为 1%~5%。 1.2 试验测试系统 为得到颗粒的管道断面图像,设计并建立了固体 物料管道输送系统,该系统由管道输送系统、动力控制 系统以及试验测量系统 3 部分组成(见图 1),系统管 道内径 50 mm,全长 20 m。 图 1 管道水力输送试验系统 在水平管中部安装一段长 2 m、内边长同为 50 mm 的有机玻璃方管。 在方管两侧设有 50 mm 50 mm 的 正方形视窗,以便用单反相机获取断面颗粒分布图像, 单反相机位于视窗的一侧(见图 2)。 在视窗的正上方 图 2 水平观察段装置 架设有大功率的连续激光器,激光器发出线状光线,增 强光线上颗粒反射光的强度。 1.3 试验方法 试验时,先启动输送动力泵,泵转速由变频调速器 控制,待系统运行稳定后,向料仓中均匀投放设计浓度 和组成的固体物料,输送速度和物料浓度数据可通过 取样标定得到。 启动测试系统,利用单反相机拍摄获 得颗粒在管道断面的分布图像。 2 IPP 图像信息提取方法 2.1 断面图像分析流程 通过上述试验可以获得管道断面图像,信息提取的 目的则是将位于激光面上的颗粒准确地显示出来。 单 张处理后的图像可以提供颗粒的浓度分布信息,多张连 续处理后的图像则可用于测量颗粒运动速度等参数。 现利用 IPP 软件对图像中颗粒信息进行提取,分 析步骤如图 3 所示。 图像分割 图像拍摄单反相机 标尺设定 分离连接部分、填补颗粒孔隙 滤镜除噪、去除小颗粒点 结果输出及效果评估 待分析断面 图 3 IPP 图像分析流程 2.2 图像分割处理及其效果分析 图 4 为在不同工况下(见表 1)拍摄得到的图像, 由于激光为宽 3 mm 的线光柱,因此仅有沿着流速方 向的一个断面上的颗粒反射光较强,在图像中表现出 较明亮的光点。 图像分割指将图像分成若干个特定的、具有独特 性质的区域并提出感兴趣部分的技术和过程[13]。 本 文图像分割的目的是将激光面上的颗粒提取出来,从 而使不在激光面上的颗粒作为背景消失。 IPP 提供了 多种图像分割方法,其中对于真彩色图像,基于 RGB 颜色空间模型的分割方法最为适宜。 RGB 分别是红、绿、蓝三基色的字母缩写。 这 3 种基色不同比例的混合叠加可以得到人类视觉中所有 可见的颜色[14]。 红、绿、蓝三基色被分为了 256 个级 别,数值越大表明该基色强度越大。 6矿 冶 工 程第 37 卷 ChaoXing 图 4 断面颗粒分布图像 表 1 图像对应的管道输送工况 图像编号粒径/ mm流速/ (ms -1 )体积浓度/ % a0.52.52 b0.52.03 c1.02.05 d1.01.53 e3.02.05 f3.02.54 本实验中激光颜色恰为绿色,拟利用绿色光的强 度(G)对图像进行分割。 以图 4(a)为例,调用 IPP 中 的强度分析工具,可以得到整幅图片中 G 值的强度分 布,见图 5。 由图 5 可见,绿色光 G 值绝大部分分布在 10~100 区间内,与拍摄图片中大面积区域颜色暗淡相 对应。 在 250~255 区间内的少量分布即为激光光柱 内的颗粒反射光较强所致。 为了利用 G 值提取目标颗粒,需要了解光柱内颗 粒的 G 值范围,在图 4(a)中随机选取 5 个光柱内的颗 粒(见图 6),并对其 G 值强度分布进行分析,结果 见图 7。 G 值 120 100 80 60 40 20 0 500100150200250 像素点个数/ 104 个 图 5 绿色光强度分布直方图 图 6 随机选取的光柱内颗粒 G 值 500 400 300 200 100 0 60120180240 像素点个数/ 104 个 图 7 光柱内颗粒的 G 值分布 由图 7 可知,由于颗粒中部的大部分区域反射光 较强,而颗粒边缘反射光较弱,所以绝大部分像素点的 G 值大于 240,少量像素点的 G 值小于 240,且在 140~ 240 区间内分布较为均匀。 为了完整地提取光柱内的颗粒而过滤光柱外的颗 粒,可调节整幅图的 G 值,将 G 值强度在一定阈值内 的颗粒选取出来,呈现白色,而在范围外的区域将会变 为黑色,实现二值化处理。 图 8 是图 4(a)经过不同阈 值(见表 2)处理后的结果。 从图 8 可知,G 值的阈值 变小,提取出颗粒的数量明显减少,这与直方图中绿色 光的强度分布一致。 7第 5 期廖 帅等 基于 IPP 图像软件的管流中粗颗粒运动信息提取方法研究 ChaoXing 图 8 经不同阈值处理后的图像 表 2 图像处理阈值范围 图像编号阈值 a[140,255] b[160,255] c[180,255] d[200,255] e[220,255] f[240,255] 2.3 标尺设定 IPP 软件中默认的尺寸单位为像素,为测量和计 算带来不便,所以需要设定标尺引入实际尺寸。 菜单 栏中依次选择“Measure→Calibration→Spatial” 打开 “Spatial Calibration”窗口即可新建尺寸,见图 9。 2.4 分离连接部分、填补颗粒孔隙 当流速较低或固相体积浓度较高时,颗粒在管道 底部发生接触和摩擦,在图像中则表现为亮点的相连; 当颗粒表面有褶皱时,会造成经过上几步处理后颗粒 内部出现孔隙,为直观地表示颗粒的形态特征,需要将 连接的颗粒分开以及颗粒内部的孔隙补全。 图 9 标尺设定界面 膨胀和腐蚀是二值形态学最基本的一对运算。 腐 蚀是消除物体边界点的过程,膨胀是将与物体有连接 (或接触)的所有的背景点都合并到物体中的过程。 IPP 提供了腐蚀和膨胀两种滤镜,将两颗粒连接部分 或颗粒内部空洞选定为 AOI 区域(area of interest),联 合使用腐蚀和膨胀两种滤镜,可以达到分离和填补的 目的(见图 10)。 图 10 腐蚀、膨胀滤镜使用前后对比 2.5 滤镜去噪及去除小亮点 经过上述图像处理后,图像中出现很多面积极小的 亮点,这些点并不能代表在激光面上的颗粒,可能是在 拍摄图像、传输图像的过程中,由于拍摄设备、传输装 置、传输路径存在约束性,导致非常容易被外部环境的 各项因素所干扰,进而形成大量的噪声[15]。 IPP 软件提 供了多种滤镜可供除去图片中细小的噪声点,经过实践 发现,对于具有暗背景和亮颗粒的图像,分离滤镜可消 除颗粒表面较小的凸起、小的噪点和暗点(见图 11)。 图 11 滤镜除噪结果 8矿 冶 工 程第 37 卷 ChaoXing 经过滤镜除噪后,图像中仍存在较小的亮点,此为 不完全在激光面上的颗粒部分反射的结果。 本实验使 用的石英砂颗粒为 2 个不同目数的筛子筛分所得,如 0.5 mm 的颗粒筛分目数为 30、32(筛孔边长分别为 0.500 mm 和 0.550 mm)。 IPP 软件可测量不规则颗粒 的多种尺寸参数,本文选择对“size(width)” 进行测 量,即在所有方向上与颗粒投影轮廓相切的两平行 线间距离的最小值(见图 12),并认为“size(width)” 值小于0.500 mm 的颗粒点为激光面外的点,将这些 点去除。 图 12 size(width)示意 2.6 结果输出 经过以上几步处理后得到如图 13 所示的结果图, 图中颗粒分布清晰,可以直观地判断出颗粒的分布等 信息。 图 13 图像处理结果 2.7 提取实例分析 根据上述分析步骤,采用不同阈值对 5 张图像进 行处理(粒径 0.5 mm,流速 2.5 m/ s,体积浓度 2%),并 利用 IPP 软件的测量功能统计提取出的颗粒数量,统 计结果与真实值进行对比,结果见表 3。 由表 3 可见, 针对该工况下拍摄得到的图像,采用 G 值的不同阈值 处理后,相较而言,当阈值取[220,255]时,处理效果 最佳,与实际光柱内的颗粒数相比,误差平均值为 4.90%,可认为达到了准确提取亮点颗粒的基本要求。 表 3 颗粒数目统计表 图像 实际 颗粒数 不同阈值处理后提取颗粒数 [140,255] [160,255] [180,255] [200,255] [220,255] [240,255] 1146210192177166154136 2160211195184172159142 3167255231211192181164 4151232217199178164148 5185258237229218191172 相对误差 平均值/ % 44.1032.5023.6014.504.905.80 3 结 语 介绍了一种基于图像处理软件 Image⁃Pro Plus 的 管流中粗颗粒运动信息提取方法。 该方法利用 IPP 软 件输入、图像分割、标尺设定、滤镜处理等功能,降低图 像分析的难度,能以简便的操作将断面图像中激光线 光柱内的颗粒提取出来,并除去噪点,分离相连的颗粒 和填补颗粒内的孔隙。 经过实践摸索,总结使用 IPP 分析颗粒物质时应该注意以下几点① 在图像的拷 贝、转发、裁剪过程中,应保证图像信息不丢失,IPP 对 高像素图片中颗粒边缘的处理精度更高;② 图像分割 过程中,G 值的阈值对处理结果的影响较大,应用时应 根据实际图像选取最佳阈值;③ IPP 软件可测量不规 则颗粒多种尺寸参数,在去除小颗粒点时,针对不同方 法筛选出的颗粒,可选用不同的测量参数。 通过实例分析,基于 IPP 处理后的断面图像能准 确直观地展现出颗粒的分布情况。 处理后的图片可用 于分析颗粒浓度分布以及测量颗粒运动速度等运动特 征,可广泛应用于管流图像处理中。 参考文献 [1] 王英伟,曹华德,曹 斌,等. 含粗颗粒固体物料在垂直管流中水 击规律[J]. 应用基础与工程科学学报, 2014,22(5)909-915. [2] 曹 斌,夏建新,黑鹏飞,等. 粗颗粒在复杂空间形态管道水力输 送中运动形式的研究[J]. 矿冶工程, 2012,32(2)5-10. [3] 肖 红,唐达生. 砂矿颗粒沉降运动规律试验研究[J]. 矿冶工 程, 2015,35(3)1-3. [4] Lu Zhanmin, Li Guangda. A Compound distribution method for LDA measurements of two⁃Phase flow and its application in measuring bubble⁃water flow[J]. Proceedings of the International Symposium on Multiphase Flows, China, 1987116-121. [5] 刘青泉. 水⁃沙两相流的激光多普勒分相测量和试验研究[J]. 泥 沙研究, 1998(2)74-82. [6] 许联锋,陈 刚,李建中,等. 粒子图像测速技术研究进展[J]. 力 学进展, 2003,33(4)533-540. [7] 张彩霞. 数字图像处理技术的发展现状及发展趋势研究[J]. 计 算机光盘软件与应用, 2014(12)216-218. [8] 袁长斌. 数字图像处理技术的发展现状与发展趋势探析[J]. 电 子技术与软件工程, 2016(2)101-102. (下转第 14 页) 9第 5 期廖 帅等 基于 IPP 图像软件的管流中粗颗粒运动信息提取方法研究 ChaoXing 位移较小。 坝体加速度随坝高有相应的放大系数,竖 直方向加速度放大系数为 2.86,与规范规定的抗震烈 度为 7 度时的放大系数接近。 4) 对于尾矿坝,抗震烈度为 7 度时,尾矿库在 3 s 时局部发生液化,且主要发生在初期坝坝顶、上游积水 区域,说明积水区沙土处于饱和状态,地震发生时孔隙 水压力上升较快,易发生液化;在 6~20 s 液化区范围 基本没有变化。 该尾矿坝在地震作用下只发生小面积 液化,因此尾矿坝整体稳定性较好。 参考文献 [1] 于广明,宋传旺,潘永战,等. 尾矿坝安全研究的国外新进展及我国 的现状和发展态势[J]. 岩石力学与工程学报, 2014,33(增 1) 3238-3248. [2] 张 超,杨春和. 尾矿坝液化判别简化方法研究[J]. 岩石力学与 工程学报, 2006,25(增 2)3730-3736. [3] 何秉顺,丁留谦,谢 源. 某尾矿坝稳定性分析与加固措施[J]. 防震减灾工程学报, 2009,29(2)179-184. [4] 潘建平,孔宪京,邹德高. 尾矿坝地震液化稳定的简化分析[J]. 水利学报, 2006,37(10)1224-1229. [5] 孔宪京,潘建平,邹德高. 尾矿坝液化流动变形分析[J]. 大连理 工大学学报, 2008,48(4)541-545. [6] 高艳平,王余庆,辛鸿博. 尾矿坝地震液化简化判别法[J]. 岩土 工程学报, 1995,17(5)72-79. [7] 陈育民. FLAC3D基础与工程实例[M]. 北京中国水利水电出版 社, 2013. [8] 张 峰,郭晓霞,杨昕光. 爆破地震波作用下尾矿坝的有限元动力 分析[J]. 防灾减灾工程学报, 2010,30(3)281-286. [9] ZBJ1-90, 选矿厂尾矿设施设计规范[S]. [10] DL 5073-2000, 水工建筑物抗震设计规范[S]. 引用本文 曹进海,胡 军,姜 鸥,等. 鞍钢齐大山尾矿坝动力响应分 析[J]. 矿冶工程, 2017,37(5)10-14. (上接第 4 页) 参考文献 [1] 陈庆发. 隐患资源开采与采空区治理协同研究[D]. 长沙中南大 学资源与安全工程学院, 2009. [2] 唐 硕. 凡口铅锌矿垮塌区隐患资源开采方案研究[D]. 长沙中 南大学资源与安全工程学院, 2014. [3] 朱坤磊. 近空区资源开采安全性分析与治理方法[D]. 长沙中南 大学资源与安全工程学院, 2012. [4] Fumagalli E. 静力学与地质力学模型[M]. 北京水利电力出版 社, 1979. [5] 关振长,龚振峰,陈仁春,等. 基于正交设计的岩质相似材料配比 试验研究[J]. 公路交通科技, 2016,33(9)92-98. [6] 张定邦,周传波,贺 丹,等. 超高陡边坡与崩落法地下开采物理 模型相似材料研制[J]. 中南大学学报(自然科学版), 2013,44 (10)4221-4227. [7] 马芳平,李仲奎,罗光福. NIOS 模型材料及其在地质力学相似模 型试验中的应用[J]. 水力发电学报, 2004,23(1)48-51. [8] 苏 伟,冷伍明,雷金山,等. 岩体相似材料试验研究[J]. 土工基 础, 2008,22(5)73-75. [9] 左保成,陈从新,刘才华,等. 相似材料试验研究[J]. 岩土力学, 2004,25(11)1805-1808. [10] 王汉鹏,李术才,张强勇,等. 新型地质力学模型试验相似材料的 研制[J]. 岩石力学与工程学报, 2006,25(9)1842-1847. [11] 顾大钊. 相似材料和相似模型[M]. 徐州中国矿业大学出版社, 1995. [12] 肖 杰. 相似材料模型试验原料选择及配比试验研究[D]. 北 京北京交通大学土木建筑工程学院, 2013. [13] 袁宗盼,陈新民,袁 媛,等. 地质力学模型相似材料配比的正交 试验研究[J]. 防灾减灾工程学报, 2014,34(2)197-202. 引用本文 高 峰,熊 信,曹建立,等. 基于诱导冒落开采模型的相似 材料配比试验研究[J]. 矿冶工程, 2017,37(5)1-4. (上接第 9 页) [9] 伍喜媛,秦 莉,陈宇玺. Image⁃Pro Plus 在眼科角膜上皮愈合组 织面积测量中的应用[J].中国医学教育技术, 2010,24(2)183- 185. [10] 陈丹明,赵乃亮,郑兴明,等. IPP 图像分析技术在长霉面积定量 测量中的应用[J]. 装备环境工程, 2012,9(6)78-80. [11] 张季如,祝 杰,黄 丽,等. 土壤微观结构定量分析的 IPP 图像 技术研究[J]. 武汉理工大学学报, 2008,30(4)80-83. [12] 余徐润,周 亮,荆彦平,等. Image⁃Pro Plus 软件在小麦淀粉粒显 微图像分析中的应用[J]. 电子显微学报, 2013,32(4)344-351. [13] 杨康叶,邬春学. 基于 RGB 模型颜色相似性的彩色图像分割 [J]. 计算机系统应用, 2013,22(3)128-131. [14] 杨 璟,朱 雷. 基于 RGB 颜色空间的彩色图像分割方法[J]. 计算机与现代化, 2010(8)147-149. [15] 尚怡君,詹保坡. 图像处理中消除噪声技术研究[J]. 科技传播, 2017(4)63-64. 引用本文 廖 帅,曹 斌,夏建新. 基于 IPP 图像软件的管流中粗颗粒 运动信息提取方法研究[J]. 矿冶工程, 2017,37(5)5-9. 41矿 冶 工 程第 37 卷 ChaoXing
展开阅读全文