付 雯 李 响(重庆电子工程职业学院软件学院,重庆 401331)
随着食品机械加工技术的不断发展,食品机械加工中,越来越多的出现了设备对接的工艺过程,该项技术受到日益广泛的关注[1,2]。通过食品机械加工对接平台可将食品加工机械的不同区域进行准确对接,达到对接准确度的要求[3-5]。因此,食品机械加工对接平台的设计已经成为相关学者研究的重点课题[6,7]。
食品机械加工对接过程中,外界的环境因素过于复杂,导致当前采用信号传感设备对食品机械加工对接平台的对接过程进行监控,会受到较大干扰,精度较低,无法达到机械对接的精度控制要求[8-10]。
针对传统方法的弊端,本设计提出一种基于三维图像校对的食品机械加工对接平台设计方法,通过分析三维图像校对原理,依据三维图像局部亮度、对比度、结构信息相关性的比较,求出其结构相似度。在大量对接图像中采集对接区域的特征点,对对接区域进行三维图像还原。获取几何阵关于z轴的子矩上的任意阶矩,还原食品机械加工对接平台的对接信息。分析投影矩的定义及计算过程,通过投影矩实现基于三维图像校对的食品机械加工对接平台设计,针对目标函数,通过Powell方法对其进行求解,和最小值相应的解即为食品机械加工对接平台的对接结果。仿真试验结果表明,所提方法具有很高的效率及对接精度。
三维图像校对过程实际上就是利用三维视觉方法,模拟人眼视觉校对过程,对机械对接过程进行校对。其主要思想是利用三维信息计算对接设备的吻合结构相似度,计算最大值,完成最优对接位置的计算。
依据对接设备三维图像局部亮度、对比度、结构信息相关性的比较,通过滑动窗口方法求出其大小,一个窗口代表一个子图像,结构相似度可描述成:
式中:
X、Y——用于描述原始设备图像的子图像和待校准设备图像的子图像;
l(X,Y)、c(X,Y)、s(X,Y)——分别用于描述X和Y的亮度相关函数、对比度相关函数以及结构相关函数,上述三项均相互独立;
为了避免分母是零造成的不稳定,C1、C2、C3取数值较小的常数。
μX、μY、σX、σY、σXY分别用于描述X、Y的亮度均值、标准差和协方差,其中协方差σXY可通过式(5)求出:
其中,N用于描述子设备图像的像素总量,针对多模态图像,因为对接设备图像内容有很大的差异,fSSIM(X,Y)通常不为1,甚至两幅对接设备图像成像亮度相反。对图像进行校对时,子图像间σxx<0,使得fSSIM(X,Y)不再为最大值,所以,需对式(4)进行调整,用|σXY|替代σXY,则有:
使s(X,Y)免于外界亮度信息的干扰。为了便于分析,本设计取α=β=γ=1,则经调整后的对接设备图像的滑动窗口可描述成:
通过式(7)求出所有窗口中(子图像)的fSSIM(X,Y)值,再对全部子图像进行累加平均:
式中:
M——用于描述子设备图像的数量;
Xj、Yj——用于描述第j个子图像。
归一化信息可描述成:
式中:
HX、HY——分别用于描述原图像与待校对图像的边缘熵;
HXY——用于描述上述二者的结构相似性,依次描述如下:
其中p(x)、p(y)分别用于描述原设备图像与待校对设备图像的概率密度;p(x,y)用于描述两幅图像之间的联合概率密度。
根据上述图像识别原理可以选取2个最优的对接设备,进行对接。
对接过程依据三维立体视觉原理,可塑造食品机械加工对接模型,实现食品机械加工对接平台三维图像对接校验。具体过程:
在对食品机械加工校验对接过程中,需要在前期确定的符合相似度要求的对接图像中采集对接区域的特征点,同时对对接区域进行有效的三维区域定位。食品机械加工对接过程采集的样本图像采集过程见图1。
图1中,r用于描述食品机械加工对接图像像素点,利用图像采集设备E1对机械对接平台图像的对节点进行采集,获取的对接点为r1,需同时通过E1、E2两个图像采集设备完成对接平台其他对接点的图像采集,和该图像相应的实际对接平台特征点可记作R。依据该特征点的空间位置,即可实现食品机械加工对接平台中,对接过程三维图像空间特征点的定位。
图1 样本图像采集原理图Figure 1 Sample image acquisition principle diagram
随机选择对接设备对接过程样本图像中的一点R,该点在图像采集设备E1和E2上的特征点分别为r1与r2,和上述两帧样本图像相应的投影矩阵可分别貌似成P1和P2,则有:
其中,(w1,x1,1)与 (w2,x2,1)用于描述图像采集设备采集的对接时三维图像中特征点的坐标;(Z,A,B,1)用于描述上述特征点在标准坐标系中的坐标;用于描述pm第k行第l个元素。通过式(15)完成对上述特征点的扩展变换处理:
对式(15)~(17)进行简化,则有:
通过式(18)、(19)可实现对接过程的三维图像虚拟仿真:
依据上述分析的方法,即可塑造食品机械交给对接平台识别模型,从大量对接图像中识别出关键图像,为食品机械加工对接平台设计提供可靠依据。
1.3.1 投影矩的定义及计算 针对上述分析提取的特征点,其几何矩可描述成:
可直接对矩进行计算,将其看作是食品机械加工设备对接图像的特征。然而由于计算量过大,很难将其应用于图像校对,因此,将式(22)转换成:
对式(23)括号中的表达式进行提取,定义几何阵关于z轴的子矩:
获取Ti(x,y)在i∈ {0,1,…}上的任意阶矩后,即可有效还原食品机械加工对接信息f(x,y,z)。
其中,X= Δx,Δy,Δz,sx,sy,sz,θx,θy,θz用于描述待优化参数,分别用于描述待对接机械在x,y与z方向的位移、缩放因子以及旋转角度。针对上述最小化问题,本研究通过Powell方法对其进行求解,和最小值相应的解即为食品机械加工对接平台的对接结果。
为了验证本研究基于三维图像校对的食品机械加工对接方法的有效性,需要进行相关的试验分析。试验在 Microsoft Visual C++ .Net环境下进行,通过MFC类库完成界面交互,同时结合VTK(visualization toolkit)部分功能。VTK是一种应用于3D计算机图形学、图象处理及可视化的软件系统工具。试验的原图见图2。
图2 食品加工机械图Figure 2 Food processing machinery
对同一食品机械加工对接三维模型进行2次随机的轮廓勾画,将2个机械的映射变换关系用待对接机械在x,y与z方向的位移、缩放因子以及旋转角度进行描述,对一个机械进行映射变换,然后通过本方法求出该映射变换。得到三维图像见图3。由表1可知,本方法的测试结果具有较高的精度,满足食品机械加工对接的需求。
图3 三维建模效果图Figure 3 Rendering 3dmodeling
表1 本方法测试结果Table 1 The method test results
为了进一步验证本方法的有效性,分别对本方法和传统方法进行仿真试验。由图4可知,采用本方法所需耗费时间一直低于传统方法。这是因为本方法首先对食品机械加工对接平台的对接过程进行三维仿真,再依据校对方法实现食品机械加工对接平台的设计,实现过程简单,计算量相对较低。由图5可知,采用本方法获取的准确度明显高于传统算法,说明本方法不仅具有较高的效率,还具有较高的准确度,验证了本方法的有效性。
图4 本方法和传统方法所需时间比较Figure 4 In this paper,the method and the time needed for traditional methods comparison results
图5 本方法和传统方法准确度比较Figure 5 Method and traditional method accuracy comparison results in this paper
本研究提出一种基于三维图像校对的食品机械加工对接平台设计方法,分析了三维图像校对原理,依据三维图像局部亮度、对比度、结构信息相关性的比较,求出其结构相似度。在大量对接图像中采集对接区域的特征点,对对接区域进行三维图像还原。获取几何阵关于z轴的子矩上的任意阶矩,还原食品机械加工对接平台的对接信息。分析了投影矩的定义及计算过程,通过投影矩实现基于三维图像校对的食品机械加工对接平台设计,仿真试验结果表明,所提方法具有很高的效率及对接精度。
1 李京娜,王国宏,孙少燕,等.基于修改后的结构相似度的三维图像配准[J].光电工程,2012,21(12):70~76.
2 余慧婷,张杰,潘萌.噪声对三维图像归一化互信息配准的影响[J].中国医学影像学杂志,2011,19(11):844~849.
3 刘东旭,刘洪,吕涛.基于CT数据的三维配准技术在正畸中的应用[J].华西口腔医学杂志,2010,28(2):119~123.
4 吕晓琪,李娜,张宝华,等.基于体素相似性的三维多模脑图像配准研究[J].中国医学影像学杂志,2013,14(2):146~151.
5 周柏清,任勇军.正则化最近点优化图像集匹配的三维人脸识别[J].计算机工程与应用,2014(12):169~174.
6 周世勋,张杰.灰度预分割对三维医学图像配准的作用[J].中国医学影像学杂志,2013,21(4):301~304.
7 吴晨,何建国,刘贵珊,等.基于近红外高光谱成像技术的马铃薯干物质含量无损检测[J].食品与机械,2014,30(4):126~130.
8 战旭梅,刘靖,刘萍.高效液相色谱法测定香肠中亚硝酸盐含量[J].食品与机械,2014,30(4):18~22.
9 晁学鹏.冶金铁轨对接控制过程三维图像虚拟仿真[J].电气应用,2014(20):90~93.
10 孙谦,何明一.三维重建中点云模型与纹理图像的配准[J].计算机仿真,2011,28(11):218~221.