基于matlab 的燃油喷雾图像处理方法,基于matlab的燃油喷雾图像处理方法

基于matlab的燃油喷雾图像处理方法

基于 MATLAB 的燃油喷雾图像处理方法 摘要提出了基于MATLAB的燃油喷雾图像处理方法。通过对喷 雾图像的采集, 分割和滤波处理, 将多幅图像转换融合为一幅 比较清楚的二值图像, 对拍摄的喷雾图像进行了去噪及二值化等初 步处理, 并提取其轮廓图, 求取中心最远点和边界点坐标, 计算 得到燃油喷雾的射程,提取其边缘图像, 对不规则的喷雾图像边缘进 行了最小二乘拟合, 得到了最近似于该边缘的2 条相交直线, 二者夹 角即为所要求的雾化角。根据最小二乘拟合的定义, 该方法的误差 平方和最小, 从而使雾化角的确定精度大大提高。本方法可为分析 燃油喷雾对柴油机性能的影响提供参考依据。 关键词 MATLAB 燃油喷雾 图像 角点检测 最小二乘法 正文因为柴油机的燃烧过程受喷入燃烧室内的喷雾状况控制, 所以发动机的性 能及其排气烟度在很大程度上受油束的喷雾锥角、射程及雾化质量等喷雾特性 的影响。研究能够预估发动机性能和排放的燃烧模型时也非常需要喷雾特性方 面的数据, 但有关这方面的知识却很贫乏[ 1] 。因此近些年来人们加速和促进了 柴油机喷雾特性的研究。有人利用高速摄像机研究燃油喷雾, 但必须对底片处 理后才能获得喷雾信息, 而底片处理时各种环境因素对结果有一定影响; 有人利 用激光CT技术对喷雾场进行研究[ 2] , 用红宝石激光器作为光源, 用CCD摄像机 适时的采集柴油机燃油喷雾。这种方法较为先进, 但造价较高。本试验利用自 制的频闪喷雾成像系统, 可以拍摄一定时间间隔的喷雾图像, 像素范围可达1 600 1 200, 清晰度高。由于MATLAB软件拥有强大的图像处理功能和高效的数值 计算能力, 所以本文利用MATLAB 程序来分析处理燃油喷雾图像, 以便能够比 较准确的获得燃油喷雾的许多有价值的信息。为进一步探讨燃油喷雾特性对柴油机性能的影响提供参考依据。 1 试验系统与研究方法 试验系统如图1 所示。本试验用数码相机拍摄燃油喷雾图像, 1. 喷油泵实验台2. 喷油泵3. 泵端压力传感器4. 喷油器5. 电荷放大器6. 频闪同步控制器7. 微型计算8. 闪光仪9. 闪 光灯10. 数码相机 图 1 频闪喷雾摄像系统示意图图 2 图像处理流程图 2算例 试验用喷油嘴型号为DLLA154S324, 启喷压力为20M Pa, 喷 油泵转速为300 r /m in, 喷射背压为大气状态, 利用频闪喷 雾成像系统拍摄相对于基准信号有不同延时时刻的喷雾图 像。现在读入图像, 如图3 所示。 图3拍摄的喷雾图像示例 按上述方法对每幅图像进行分割和滤波处理, 把真彩色图像 分别转换为比较清楚的黑白图像, 如图4所示。为减少不确 定因素,并提取喷雾图像的轮廓图, 如图5所示。图4处理后的二值喷雾图像图5喷雾图像的轮廓图 通过微型计算机控制闪光仪的闪光时刻, 使闪光时刻与喷油 器某时的喷雾同步, 拍摄每幅图像时闪光灯的持续闪光时间 为 0. 1 m s。为使所拍摄的图像具有可比性, 用压力传感器 获取的压力信号为基准信号, 通过调整相对于基准信号的延 时时间获取有一定时间间隔的喷雾图像, 以观察燃油雾化过 程[ 3]。由于所拍摄图像并非来自同一循环, 本试验通过在 不同循环中相对于基准信号的相同时间间隔拍 3 幅, 以便减 少不确定因素, 减小随机误差。由于拍摄图像中油束的大小 与实际大小有一定差别, 因此需要标定。方法是将数码相机 的焦距和景深调好以后, 把一已知长度的纸条挂到喷油器上, 进行拍照。将其实际长度与照片中长度相比, 即得实物与图 像的比例尺。通过比例尺换算可以得到油束的实际长度。 并且在一次试验中, 数码相机的焦距和景深调好以后不再变 动。数码相机拍摄的喷雾图像为真彩色图像, 即 RGB 图像, 每幅图有 1 600 1 200 个像素。图像的读入由 MATLAB 程 序的 imread 函数来实现。执行该函数后, 得到一个代表该 图像的数字矩阵。为便于处理, 首先把真彩色图像转换成灰 度图像, 再运用直方图均化函数 histeq 处理, 使图像细节更 加清楚, 然后设置阈值, 对图像进行分割; 为了强化图像边 界和消除噪声的影响, 对灰度图像进行 LOG 算子滤波, 并 用二值分割图像进行与运算, 然后利用 MATLAB 的相关函 数选择所得二值图像中特定的对象,此时得到一幅黑白图像, 对其余同一延时时刻的 2 幅图像, 均作同样处理, 共得到 3 幅黑白图像, 将 3 幅黑白图像进行与运算,得到一幅包含全 部输入图像中特定对象的比较清晰的二值图像。运用 MATLAB 中jharris角点检测函数提取图像的角点,然后再 利用 g 函数提取出油束的中心点与最远点，计算中心 点与各油束最远点的距离作为该油束的射程,计算公式为 .然后采用 g 函数沿油束边界取五个 点的坐标, 用最小二乘法拟和成直线作为油束边界切线, 两 条切线的夹角作为油束的喷雾锥角。计算公式为 图像处理的程序框图见图 2。 提取喷雾图像中心和各孔最远点及每孔油束各边界上5 个 点的坐标, 计算得出各孔油束在某一时刻的射程和喷雾锥角 见表1。 表 1 喷油嘴各孔油束的射程和锥角 3结束语 本文基于MATLAB语法结构简单、数值计算高效和图像处 理方便的特点, 研究了燃油喷雾图像处理方法, 能比较准确 的计算喷油嘴各孔油束的射程和锥角, 所编写的程序方便易 懂, 具有一定实际应用价值。 参考文献 [ 1] 孙银城. 喷油器[M ]. 北京 国防工业出版社, 1981. [ 2] 周龙保. 内燃机学[M ]. 北京 机械工业出版社, 2003. [ 3] 李聪, 刘建新, 杜慧勇. 用数码摄像技术研究柴油机喷 雾特性[ J] . 洛阳 工学院学报, 2002, 23 4 53 55. [ 4]孙兆林.MATLAB 6. x 图像处理[M ] . 北京 清华大学出 版社, 2002 1孔2孔3孔4孔喷油嘴型号 射程 锥角 /mm /0 射程 锥角 /mm /0 射程 锥角 /mm /0 射程 锥角 /mm /0 DLLA154324 167.27 15.3334 204.18 17.53 162.08 15.05 159.70 11.72[ 5 ] 陈 熙. 环形出口截面的离心喷嘴的简单分析[J ]. 力 学学报, 1978 4 3322339. [ 6 ] 胡传胜, 熊洪亮, 虞先煌. 喷嘴液滴雾化细度和喷雾角 测量装置的研究[J ]. 能源研究与信息, 2001, 17 4 2252231. [ 7 ] 庞 博, 王振清. 基于M atlab 的数据处理与三维模拟 [J]. 微型计算机信息, 2004, 20 1 1122113. [ 8] 孙兆林. MA TLAB 6. x 图像处理[M ]. 北京 清华大 学出版社, 2002. [ 9] 朱美玲, 赵淳生. 一种万能的数学运算“演算纸” MA TLAB