您好,欢迎来到化拓教育网。
搜索
您的当前位置:首页激光全息干涉法测量液相质扩散系数数字图像处理的研究

激光全息干涉法测量液相质扩散系数数字图像处理的研究

来源:化拓教育网
维普资讯 http://www.cqvip.com 检测与仪表 化工自动化及仪表,2008,35(5):46—50 Control and Instruments in Chemical Industry 激光全息干涉法测量液相质扩散系数 数字图像处理的研究 费继友,涂娟 (大连交通大学机械工程学院,大连116028) 摘要: 用He—Ne激光器、CCD和图像采集系统等仪器设备,搭建了数字激光全息干涉法测量液相质扩散系 数的实验系统,并测试获得了浓度为O.33 mol/1 KC1水溶液在298.15 K温度下的全息干涉图像,采用了图像增 强、全息干涉条纹图再现及相位展开等数字图像处理过程的关键技术对全息干涉图进行处理,得到了相位展开后 物光相位差,经换算后获得了该溶液的质扩散系数。为化工工程上研究替代制冷剂急需的质扩散系数测量提供 了有效的方法。 关键词:数字激光全息干涉法;数字图像处理;质扩散系数;干涉图;水溶液;KC1 中图分类号:TP751 文献标识码:A 文章编号:1000—3932(2008)05-0046-05 1 引 言 全息干涉条纹图再现及相位展开等数字图像处理 过程的关键技术。 最后测量浓度为0.33 mol/1 KC1水溶液在 298.15 K下的扩散系数,验证实验方法的准确性和 为了应对全球性的气候变暖、生态破坏等现 象,《蒙特利尔议定书》第l9次缔约方大会达成 了氢氯氟烃(HCFC )于2030年完成生产量与消 费量的淘汰。而当前国际上普遍采用的替代制 冷剂是HFC ,但是HFC 也属于《京都议定书》中 应实施减排的温室气体。目前市场上大量使用 的R22制冷剂也将在2030年淘汰,即使是替代 R22最理想的R404A仍然属于温室减排气体而 可靠性,为化工工程上研究替代制冷剂急需的质扩 散系数测量提供了有效的方法。 2液相质扩散系数的测定 2.1 液相质扩散系数的测量原理 Fick第二定律是质扩散系数实验测量和理论推 面临淘汰。面对全球加速推动温室气体减排、减 缓全球变暖趋势、保护臭氧层以及我国绿色环保 可持续发展战略的实施,化工制冷行业迫切需要 找到新的环保制冷剂。而液相分子扩散系数是 描述质量传递的最重要的物性参数,是化工上替 代制冷剂研究急需的主要参数之一,本文通过研 究数字全息干涉法测量液相质扩散系数 原理, 算的理论基础。当扩散被认为是一维时,并假设在 扩散开始前两溶液有明显的分界面。由Fick第二 定律可得到质扩散系数D: .t. /t.一1 D ‘ ) 式中: ——溶液浓度变化极值点间的距离;t。, t.——前后两张全息图的曝光时间。 用CCD传感器取代双曝光全息干涉法和实时全 息干涉法中的感光干板,物光和参考光由分光棱 由式(1)可知质扩散系数D主要与参数t。,t。和 的测量有关。t。和t 可以通过计算机硬件时钟 获得,而物光相位差△‘D的变化和溶液浓度△c 的变 化之问的关系为: = 镜聚合在一起,直接使用CCD传感器记录干涉条 纹图。数字全息干涉法的相于光波是采用时间 分隔法而获得的,随着时间的推移,物体的实际 状态发生变化,干涉条纹图也随之相应地发生变 化。由于折射率和溶液浓度存在着一一对应关 系,可以通过分析不同时刻的干涉条纹图获得溶 液浓度的变化,通过换算进而得到质扩散系数。 本文利用MATLAB软件对数字图像采集系统 记录的全息图进行分析处理,清晰再现了干涉条纹 图,应用了全息图的预处理、数字图像的频域滤波、 Ac22__ ̄rl_ (2) 式中: ——浓度变化和折射率变化之问的比例系 数;A——物光波长;Z——扩散槽厚度。 由物光相位差的变化找到浓度变化的极值点 ,最终求出了质扩散系数D。 收稿日期:2008-08.11(修改稿) 基金项目:国家自然科学基金资助项目(50676069) 维普资讯 http://www.cqvip.com 第5期 费继友等.激光全息干涉法测量液相质扩散系数数字图像处理的研究 。47・ z.2实验台的搭建 2.2.1数字激光全息干涉法质扩散系数测量原理 数字激光全息干涉法测量液相质扩散系数 人内存,再写入硬盘。此时全息干涉条纹图像可 以在计算机上进行去噪、细化等预处理,并对图 像进行特征提取,从而使采集与图像的处理、分 析同步进行,实现图像处理的实时化。 的测量光路 如图1所示。由He—Ne激光器发 出的激光经平面镜反射,激光束进入空域滤波器 形成平行光,通过挡光板后只保留亮度比较均匀 2.3数字全息干涉图像的处理 2.3.1全息干涉图的处理过程 的部分以供使用;然后通过分光棱镜将激光束分 液相质扩散系数D的获得关键在于获取物光 成两束,直透光束为物光,反射光束为参考光;物 光通过扩散槽后和参考光都投射到CCD上,全息 图通过图像采集卡由计算机采集和存储,就可以 实时得到不同物质、不同时刻的干涉条纹。由测 试原理可知通过对于涉条纹数字图像处理找到 浓度变化的极值点即得到了质扩散系数。 图1 数字激光全息干涉法质扩散系数测量原理图 注:l——He.Ne激光器;2、5、8——平面镜;3——空域滤波器;4、 7——分光棱镜;6——扩散槽;9——ccD传感器;1O——带有 图像采集卡的工控机 物光通过的扩散槽底部设计了锥形的整流 栅 ,防止注液时带来的湍流和返混。为使实验 在稳定的温度环境下进行,将扩散槽置于恒温系 统中。 2.2.2数字图像自动采集系统 数字图像自动采集系统结构如图2所示。 全息干涉图像I CCD传感器 l 亘 重 l 带有图像采集卡、存储软件和图像处理软件的工控机 图2数字图像自动采集系统 该系统由CCD传感器,带有图像采集卡、存 储软件和图像处理软件的计算机组成。采用 CCD传感器接收全息干涉图像信息,由图像采集 卡采集后送至计算机进行处理和分析,采集后获 得了分辨率为2048×1 536的全息干涉图像。全 息干涉图像的采集与处理均在计算机的控制下 进行,可以实现图像的实时、快速采集。全息干 涉图像存入帧存器后存人硬盘之前,需要在内存 区开辟一个5 12 K区域。然后按段地址将数据读 在不同时刻的相位差信息,全息干涉图像相位差的 提取过程如图3所示。 /( ,Y) ( ,Y) 输入图像 增强后的图像 图3数字全息干涉图像的处理过程 2.3.2全息干涉图的预处理 全息干涉图的预处理必须采取图像增强技术, 目的是为了改善图像的视觉效果,提高图像的清晰 度和工艺的适应性,以及便于人与计算机的分析和 处理,满足图像复制或再现的要求。 图像增强的方法分为空域法和频域法两类, 空域法主要是对图像中的各个像素点进行操作, 并修改变换后的系数,如傅立叶变换,离散余弦 变换等的系数,然后再进行反变换,便可得到处 理后的图像。 (1)灰度变换增强。 对于一幅离散的图像,第i个灰度级 出现的 频率用n 表示,该灰度级像素对应的概率值P (r ) 为: Pi(ri)=鲁 (3) 式中:n——像素总数;ri——满足归一化条件。图 像变换函数表达式为: S =T(ri= )= )= (r(r  )=k)= ~  -I~ni(4) n 式中: ——灰度级数;Si——处理后的图像中的灰 度级,对应于输入图像中的灰度级 图4为全息干涉原始图像及其直方图,图5为 均值化后的全息干涉图像及其直方图。不难看出, 用直方图均值化修改后,图像的直方图的灰度间隔 被拉大了,从而有利于图像的分析与识别。 (2)空域滤波增强 。 空域滤波按照空域滤波器的功能又可分为平滑 滤波器和锐化滤波器。平滑滤波器可以用低通滤波 维普资讯 http://www.cqvip.com ・48・ 化工自动化及仪表 第35卷 实现,目的是模糊图像或消除噪声;锐化滤波器使用 互相干项,经过移频后它们出现在图像的中心位置 高通滤波来实现,目的是强调图像被模糊的细节。 ×104 (频谱中心),C,是原始物光波前的准确再现,如图 7(b)所示。 O 5 O l-一 图4全息干涉原始图像及其直方图 ■一 (a)全息干涉图傅立叶 (b)滤波后的全息干 0 ll l巍 隧“ l∞ 图5 均值化后全息干涉图像及其直方图 平滑滤波器中medfilter2()函数用于实现中值 滤波,中值滤波是一种典型的低通滤波器,主要目的 是保护图像边缘,同时也能去除噪声。它的基本原 理是把数字图像或数字序列中某一点的值用该点的 一个邻域中各点值的中值代换。设 ,Y)表示数字 图像像素点的灰度值,滤波窗口为以的中值滤波器 可以定义为: G( ,y)= d tf( ,y)} (5) 【 .yJ eA 图6为中值滤波后的全息干涉图像,从上面几 幅图中可以看出,经过图像增强、滤波等预处理后的 全息图条纹更加平滑,对比度也明显增强,更有利于 进一步的数字处理。 图6 中值滤波后全息干涉图像 2.3.3全息干涉图的频域滤波 图3中给出了频谱滤波的主要步骤,对全息干 涉图进行二维傅立叶变换,可得到全息干涉图的傅 立叶频谱信息(如图7)。傅立叶频谱信息反映物光 的真实相位,它可以通过乘以一个滤波函数被保留。 对全息干涉图进行傅立叶变换可以得到其频谱 分布情况如图7(a)。其中,c.代表参考光的直透部 分,c:代表物体各点的自相干以及物体各点之间的 变换的频谱信息 涉图频谱信息  图7频谱信息 根据线性系统理论及卷积的运用,滤波器的表 达式为G( , ):H( , )c( ,叼),其中日( ,卵)为 滤波函数,C( , )是图像f( ,Y)的傅立叶变换, ( ,Y)是经滤波后输出图像的傅立叶变换。实际 上,物体的频谱一般主要集中于低频部分,而且在频 谱的中心部分强度很大,集中了很大一部分能量;相 对而言,其他的频谱成分集中的能量要小很多。所 以滤波器的选择非常重要,在这里我们选择二阶巴 特沃斯低通滤波器。 2.3.4再现数字全息干涉条纹 对在两个不同时刻得到全息图的频谱信息进行 反傅立叶变换,可以得到原始物光波的真实相位 。 然后,截断相位 ( ,Y)可以从c( ,Y)解出。先将c ( ,Y)和c ( ,Y)作点除运算,即: r ,、 Q( ,,,)=; ÷ =exptj[妒 ( ,y)一 t( ,y)]} =exp[jA ̄p( ,y)] (6) 再用对数算法从Q( ,y)中解出△ ( ,y),即得 到物光相位变化: log[Q( ,y)]=jA ̄p(x,y) (7) 通过以上处理可以再现图8所示的数字全息干 涉条纹图。 图8再现数字全息干涉条纹图 2.3.5全息干涉图的相位展开 图8中所提取的相位值 ( ,Y)被截断于一耵 和 之间,从而呈现阶跃不连续的相位分布图样, 一维普资讯 http://www.cqvip.com

第5期 费继友等.激光全息干涉法测量液相质扩散系数数字图像处理的研究 ・49・ 如图9所示。为得到相位的实际值,必须将主值相 位恢复为绝对相位,这个过程称为相位展开。本文 提出了基于离散泊松方程解的相位展开方法 ,它 借助快速离散余弦变换求解M×N具有诺埃曼边界 条件下的离散泊松方程,得到最小意义下的解,进而 得到展开的相位分布。假定 为二维离散点( , ) 要包括物光相位分布的展开和获取物光真实的相位 差。通过测量相位展开后的物光相位差(图10所 示)极值点间的距离也就得到了在扩散方向上浓度 变化极值点间的距离 ,拍摄全息图的时间t。、t 可通过图像自动存储软件获得,这样由公式(1)便 可换算出质扩散系数。 3实验测试结果与结论 上的展开相位, 伪测量到的主值相位,首先定义包 A = +1)一 ;其中上标 ,Y分别指行像素 一 裹算子,有W( )= ;定义包裹相位差△;= ( )f 3.1实验测试结果 ,为验证激光全息干涉装置的可靠性,实验所测 相位差和列像素相位差。则包裹相位差与展开相位 差的关系式为: =[△;一△ 一,)J]+[△ 一△ ,一 )] (8) 其中,i=0,…,M~1; =0,…,Ⅳ一1;△ _1) :0; △ _I],=0;△ _1)=0;△ 一 H=0;而展开相位咖 在离 散余弦变换域的精确解为: 中 . 万 f 可 (…  ) 对由式(8)计算得到阵列 进行二维离散余弦 变换正变换,产生二维离散余弦变换谱值 ,通过 式(9)计算得到咖 作咖 的二维离散余弦反变换, 得到最小二乘意义下的展开相位 相位展开前 后的物光相位差信息分别如图9、图10所示。 图9相位展开前的物光相位差信息 2.3.6图像的后处理 对于激光全息法测量液相质扩散系数测试系统 所获得的数字全息干涉图像而言,全息图后处理主 浓度为0.33 mol/l的KC1水溶液在298.15 K温度下 的质扩散系数如表1所示,总相对误差 =0.9%, 可见此激光全息干涉装置是精确的和可靠的。 表1 KCL水溶液在298.15K下的质扩散系数 注:D——实验值;D——实验平均值;D ——文献参考值 3.2结论 (1)利用He—Ne激光器、CCD和图像采集系统 等设备搭建了基于数字激光全息干涉法测量液相质 扩散系数的实验系统。 (2)通过该实验系统得到了全息干涉图像,采 用了全息图的预处理、数字图像的频域滤波、全息干 涉图再现及相位展开等数字图像处理的技术对全息 干涉图进行了处理,获得了相位展开后的物光相位 差,为求得质扩散系数奠定了基础。 (3)使用该实验系统,测量了浓度为0.33 mol/1 KCI水溶液在298.15 K温度条件下的质扩散系数, 与文献值相比验证了实验系统的可靠性,为化工工 程上研究替代制冷剂急需的质扩散系数测量提供了 一个有效的方法。 参考文献: [1]MA Pei・sheng,wu Yan-xia,LIU Yun—qi.Diffudion Coeficient of 1一Lysine Hydrochloride and 1一Arginie Hydrochloride in Their Aqueous Solution at25℃[J].Chinese J ofChemEng,2000,8 (2):146—153. [2]王家文,李仰军.MATLAB 7.0图形图像处理[M].北京:国 防工业出版社,2005. [3] 费继友,张运宏,张颖,等.数字全息干涉法测量液相质扩 散系数的研究[J].化工自动化及仪表,2007,34(6):54一 S7 维普资讯 http://www.cqvip.com ・50・ 化工自动化.及仪表 第35卷 [4] GRA Y L G.FENICHEL H.Holographic Interferometric Study of fusion Measurements in Mohen Salts and Ionic Liquids.method Liquid Diffusion[J].Appl Opt,1979,18(3):343—345. and First Results[J].Journal of Molecular Liquids,2003, [5] 孟令军.用计算机数字图像处理技术监测电站锅炉燃烧器 (103-104):359—370. 火焰[J].化工自动化及仪表,2001,28(2):43—45. [8] 惠梅,王东生,李庆祥,等.基于离散泊松方程解的相位展 [6] 易正明,周孑民,刘志明,等.基于图像处理的回转窑火火焰 开方法[J].光学学报,2003,10(3):1245—1249. 监测系统研究[J].化工自动化及仪表,2005,32(5):54— [9] G0 STING L J.A study of the Dilfusion of Potassium Chloride 56. in'Water at 25 oC with the Gouy Interference Method『J].J A. [7] RICHTER J,LEUCHTER A.Digital Image Holography for Dif- mer Chem Soc,1950,72(10):44l8—4422. Investigation on Holographic Interferometry Measurement of Liquid Mass Difusion Coefifcient of Digital Image Processing FEI Ji-you.TU Juan (School of Mechanical Engineering,Dalian Jiaotong University,Dalian 1 1 6028,China) Abstract:An experimental optical system based on digital image holographic interferometry measurement of liquid mass diffusion coefficients was constructed by He—Ne laser,CCD,image acquisition system and other equipments, the system was tested and the interferograms of KC1 aqueous solution were collected at 0.33 mol/1 concentration and 298.1 5 K.By employing such key techniques of digital image processing as the image enhancing,the interference fringes rebuilding and the phase unwrapping to process the interferograms,the phase difference was obtained,and 1J] 1j therefore the mass diffusion coefifcients of KC1 aqueous solution was also gained by conversion.It寸_l  can be concluded that the system can provide an effective method for measuring the mass diffusion coefifcient of new refrigerants,which is urgently required in chemical engineering application. Key words:digital image holographic interferometry;digital image processing;mass diffusion coefifcients;interfero— metirc map;aqueous solution;KC1 (上接第45页) 围为0~2×10~,分辨率小于1×10~,误差小于 OLDHAM GROUP.MX 48:Measurement and Alarm Unit 2%,响应时间小于30 s。本装置可广泛应用于造 [Z/OL].OLDHAM GROUP,[2008-03-17].http://www. 纸、纺织、制药、水处理等需要对环境的氯气浓度实 groupoldham.com/site—en/produits/detail.php?produit 时在线监测的场合。更换气体传感器的类型后,该 mx48. 装置还可检测其它有毒有害气体,具有灵活的配置 何希才,薛永毅.传感器及其应用实例[M].北京:机械工业 出版社,2003. 和广阔的应用前景。 余青山,苏宏业,董利达.基于cDMA/GPRs/RAD10技术的 参考文献: 无线环境监控终端设计[J].化工自动化及仪表,2006,33 [1] 王会彬.低浓度氯气环境对作业工人肺通气功能的影响 (3):39—43. [J].中国辐射卫生,2007,16(2):224. 曹晔,张光友,胡文祥.有毒有害和可燃气体的自动检测 [2] 王普咨,陈公北,陈 飞.一种氯气检测报警仪:中国, 技术及其应用[J].现代仪器,2001,(2):41—43. 91214946.9[P].1992-07-08. [3] 周永清.袖珍式氯气报警仪的设计[J].传感器技术,2004, 胡大可.MSP430系列单片机c语言程序设计与开发[M]. 23(11):41—42. 北京:北京航空航天大学出版社,2003. A New Type of Chlorine Gas Concentration Online Detect Device JIANG Peng,JIANG Lu—rong (Institute ofInformation and Control,Hangzhou Dianzi University,Hangzhou 310018,China) Abstract:A chlorine gas online detect instrument was developed.Adopting MSP430 processor and constant electric potential electrolytic gas sensor,it could measure the concentration of chlorine gas in the air,give an alarm when the concentration of chlorine gas is out of range and communicate with the control center by GPRS.The measurement range of the instrument is 0~20 ppm.response time is less than 30 sec and the detecting error is less than 2%.This design is feasible and can meet the real application in papermaking industry,textile industyr,pharmacy and water treatment industry. Key words:concentration of chlorine gas;online detection;sensor;MSP430 ]1j 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- huatuo9.cn 版权所有 赣ICP备2023008801号-1

违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务