1.一种燃油气泡喷嘴内外部雾化过程的检测系统,包括图像采集模块、图像处理模块、数据分析模块,所述图像采集模块将采集到的图像发送至图像处理模块,所述图像处理模块对图像进行处理后,将处理后的图像传输至数据分析模块进行数据分析得出结果,其特征在于:所述图像采集模块用于采集喷嘴内部气液两相流图像和外部喷雾场图像;
所述图像处理模块包括灰度处理单元、滤波降噪单元、二值化单元、边缘检测单元,灰度处理单元将含有亮度和色彩的彩色图像变化成灰度图像,滤波降噪单元对灰度图像进行滤波处理,获得降噪平滑处理后的图像,二值化单元通过十六种二值化方法对图像进行处理得到二值化喷雾图像,边缘检测单元对二值化喷雾图像进行边缘检测得到气泡轮廓、喷雾轮廓以及喷雾宽度;
所述数据分析模块用于检测分析气泡平均直径、雾化锥角大小和喷雾宽度经时间序列的快速傅里叶变换后得到的喷雾振荡的频率和强度。
2.根据权利要求1所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,包括以下步骤:步骤一、启动系统,图像采集模块对燃油气泡喷嘴内部和外部进行图像采集,获得喷嘴内部气液两相流图像和外部喷雾场图像;进入步骤二;
步骤二、图像灰度处理,灰度处理单元将含有亮度和色彩的喷嘴内部气液两相流图像和外部喷雾场图像处理成灰度图像;
步骤三、图像滤波降噪,滤波降噪单元采用中值滤波法将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值,进行滤波处理去除脉冲干扰及椒盐噪声,获得降噪平滑处理后的图像;进入步骤四;
步骤四、二值化图像,二值化单元对步骤三进行滤波降噪后的图像采取十六种二值化方法对其进行处理,将得到的十六种二值化图像进行对比,选择处理后图像与原灰度图像中内部气泡图或喷雾雾场轮廓最为接近的二值化方法进行批量处理;进入步骤五;
步骤五、边缘检测,边缘检测单元对经过步骤四处理后的二值化喷雾图像进行边缘检测,获得气泡轮廓和喷雾轮廓,设定一条距喷嘴出口一定距离的测量线,得到测量线上的喷雾宽度;进入步骤六;
步骤六、数据分析,数据分析模块根据步骤五得出的气泡轮廓和喷雾轮廓以及喷雾宽度,分析气泡平均直径、雾化锥角大小和喷雾宽度经时间序列的快速傅里叶变换后得到的喷雾振荡的频率和强度;结束。
3.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤二具体步骤如下:所述灰度处理单元将RGB图像在纯黑到纯白色区之间区分成256种灰度进行储存,灰度图采用8‑bit灰度图,对于灰度值的计算采用精确度最高的浮点算法:f(i,j)=0.30R(i,j)+0.59G(i,j)+0.11B(i,j)
其中,i、j为浮点坐标。
4.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤三具体步骤如下:所述滤波降噪单元采用中值滤波法将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值,将数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,让周围的像素值接近的真实值,从而消除孤立的噪声点,利用二维滑动模板,将板内像素按照像素值的大小进行排序,生成单调上升或下降的二维数据序列;
在OpenCV中,实现中值滤波的函数是cv2.medianBlur(),其语法格式如下:dst=cv2.medianBlur(src,ksize)
其中,dst是返回值,表示进行中值滤波后得到的处理结果;src是需要处理的图像,即源图像,它能够有任意数量的通道,并能对各个通道独立处理;ksize是滤波核的大小,滤波核大小是指在滤波处理过程中其邻域图像的高度和宽度。
5.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤四具体步骤如下:所述二值化单元选用伯恩森算法、迭代自组织数据分析算法、最大熵算法、均值算法、最小错误算法、最小交叉熵算法、双峰算法、模糊阈值算法、最小阈值算法、几何矩算法、大津算法、百分比阈值算法、雷尼熵算法、阈值分割算法、三角阈值分割算法、模糊阈值算法十六种二值化方法同时对所得图像进行处理,将处理后的二值化图像进行对比,选择处理后图像与原灰度图像中内部气泡图或喷雾雾场轮廓最为接近的二值化方法;针对不同喷雾实验工况所对应的真实雾化状况不同的现象,使用该图像所对应的最优二值化方法对同组工况图像进行批量处理。
6.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤五具体步骤如下:所述边缘检测单元选择处理单独的像素点所得到响应比处理边缘点所得到的响应要更强烈的Laplacian算子,进行边缘检测以获得气泡轮廓和喷雾轮廓,设定一条距喷嘴出口一定距离的测量线,得到测量线上的喷雾宽度;
在OpenCV内使用函数cv2.Laplacian()实现Laplacian算子的计算,该函数的语法格式为:dst=cv2.Laplacian(src,ddepth[,ksize[,scale[,delta[,borderType]]]])其中,dst代表目标图像;src代表原始图像;ddepth代表目标图像的深度;ksize代表用于计算二阶导数的核尺寸大小,该值必须是正的奇数;scale代表计算Laplacian值的缩放比例因子;默认情况下,该值为1,表示不进行缩放;delta代表加到目标图像上的可选值,默认为0;borderType代表边界样式。
7.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤六中的分析气泡平均直径具体步骤如下:根据步骤五中得到的气泡轮廓的长轴和短轴得出气泡等效直径:
其中,a为椭圆的半长轴,b为椭圆的半短轴;
气泡的平均直径具有相同的总体积与总表面积之比:
其中,di为单个气泡的直径mm;N是气泡的总数。
8.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤六中的分析雾化锥角具体步骤如下:对步骤五中边缘检测到的喷雾边缘轮廓进行直线拟合,使用最小二乘算法,得到两条喷雾轮廓的边界线,两条边界线间夹角即为雾化锥角。
9.根据权利要求2所述的燃油气泡喷嘴内外部雾化过程的检测系统的检测方法,其特征在于,所述步骤六中的分析喷雾宽度经时间序列的快速傅里叶变换后得到的喷雾振荡的频率和强度具体步骤如下:将步骤五中得到的喷雾宽度应用到一个时间序列的瞬态自激振荡喷雾图像中,得到一个时间序列的喷雾宽度信息;针对这个时间序列的喷雾宽度进行快速傅里叶变换,得到喷雾振荡的频率和强度。