本发明属于几何量检测技术领域,尤其涉及一种基于机器视觉的光纤几何参数检测方法。
背景技术:
目前,光线已经广泛应用在科研、电信、传感、监视等领域,是国家重点扶持和大力开展的新型产业。其中,光线的几何参数如包层直径、包层不圆度、涂覆层直径、涂覆层不圆度等对光纤传感和光纤通信等领域的研究至关重要。此外,在光纤生产中,光纤的几何参数带来的误差对光纤间的耦合效率起着决定性的影响。
目前已经提出了很多光纤几何参数检测方法,如人工判读为主的放大法、以光纤传输成像为主的折射近场法和近场扫描法等,但都存在着一定的缺陷:1)人工判读为主的放大法,检测精度低;2)折射近场法利用光纤折射率分布计算光纤参数分布计算光纤参数,近场扫描法是对近场放大像面进行扫描,测得近场强度分布扫描图,从而拟合光纤参数。这两种方法实验装置复杂,测量过程不直观。因此工业上亟需能够研究一种快速、高效、自动化程度高的检测方法具有重要意义。
技术实现要素:
本发明的目的在于提供一种基于机器视觉的光纤几何参数检测方法,其检测速度快、检测精度高。为实现上述目的,本发明所采用的技术方案为:
一种基于机器视觉的光纤几何参数检测方法,包括以下步骤:
s1:对图像采集装置获取的光纤横截面图像进行预处理:所述光纤横截面包括涂覆层、包层和纤芯;
s2:确定分割阈值,对预处理后的光纤横截面图像进行阈值分割,提取涂覆层的外边缘图像;
s3:分别运用改进的canny算法对涂覆层的内边缘和包层的外边缘进行边缘检测以提取涂覆层的内边缘图像和包层的外边缘图像;
s4:运用圆形拟合发和椭圆拟合法,对涂覆层的外边缘图像、涂覆层的内边缘图像及包层的外边缘图像均进行圆形拟合和椭圆拟合,以获取各图像对应的圆形参数和椭圆参数;
s5:分别将步骤s4中的圆形参数和椭圆参数作为输入参数,运用光纤几何参数求解模型,分别求解步骤s4中各图像对应的直径和不圆度。
优选地,在步骤s1中,所述预处理包括灰度变换、平滑处理及锐化处理。
优选地,在步骤s2中,运用最大类间方差法确定分割阈值。
优选地,所述最大类间方差法具体为:
s21:运用累加法计算像素的灰度值属于暗背景的概率p1和像素的灰度值属于亮物体的概率p2;其中,
s22:根据p1、p2计算出暗背景区域的平均灰度值g1和亮物体区域的平均灰度值g2;其中,
s23:根据平均灰度值,可计算出图像的类间方差
优选地,在步骤s3中,提取涂覆层的内边缘图像具体步骤为:
s31:运用形态学方法,由涂覆层的外边缘图像确定涂覆层的内边缘roi区域;
s32:对涂覆层的内边缘roi区域采用二维高斯函数g(x,y)对该区域的图像f(x,y)进行平滑滤波以去除噪声,得到图像fs(x,y);
s33:利用一阶微分sobel算子计算平滑后的图像fs(x,y)各点处的梯度幅值和梯度方向,从而获得相应的梯度幅值图像m(x,y)和梯度方向图像θ(x,y):
s34:在m(x,y)中以点(x,y)为中心,在设置的邻域内沿θ(x,y)进行插值,如果点(x,y)处的梯度幅值m(x,y)大于θ(x,y)方向上与其相邻的插值,则将点(x,y)标记为候选边缘点,以获取候选边缘图像n;
s35:选取阈值tl,检测候选边缘图像n中标记为候选边缘点的任一像素点(x,y);若点(x,y)处的梯度幅值m(x,y)大于其低阈值tl,则将该点标记为边缘
点,以获得阈值检测图像n1;
s36:对阈值检测图像n1,以中心坐标分成多个区域si;运用局部区域距离阈值法分别对si分别筛选,以获取每个si中符合设置条件的边缘子图像;其中,i为区域个数;
s37:连接所述边缘子图像以获取涂覆层的外边缘图像。
优选地,在步骤s4中,所述椭圆拟合法具体为:
s41:运用二次曲线标准化法对提取的边缘图像进行椭圆拟合以获取二次曲线方程,由二次曲线方程求取椭圆中间参数;
s42:采用齐次方程法以步骤s41中获取的椭圆中间参数作为齐次方程法中初始迭代值,求解椭圆参数值。
优选地,在步骤s4中,所述圆拟合法和椭圆拟合方法相同。
优选地,步骤s5中,所述光纤几何参数求解模型为高斯-牛顿迭代法模型。
与现有技术相比,本发明的优点为:
1)采用改进的canny算法,可以减少干扰边缘有利于曲线拟合;
2)同时二次曲线标准化法和齐次方程法相结合的方法对边缘进行曲线拟合,并用高斯-牛顿迭代法求解光纤几何参数,提高了边缘曲线拟合的精确度;
3)本发明对实验环境的要求低,对对应的测量系统的结构简单。
附图说明
图1为本发明一实施例的基于机器视觉的光纤几何参数检测方法的流程图;
图2为图1中光纤横截面图像;
图3为图2中光纤横截面图像预处理后的图像;
图4为图1中涂覆层的外边缘图像;
图5为运用现有技术中canny算法检测获取的涂覆层的内边缘图像;
图6为本发明改进的canny算法检测获取的涂覆层的内边缘图像;
图7为本发明涂覆层的内边缘图像的椭圆拟合结果图;
图8为本发明光纤横截面的椭圆拟合及圆形拟合结果图。
具体实施方式
下面将结合示意图对本发明的基于机器视觉的光纤几何参数检测方法进行更详细的描述,其中表示了本发明的优选实施例,应该理解本领域技术人员可以修改在此描述的本发明,而仍然实现本发明的有利效果。因此,下列描述应当被理解为对于本领域技术人员的广泛知道,而并不作为对本发明的限制。
如图1所示,本实施例提出了一种于机器视觉的光纤几何参数检测方法,包括步骤s1~s5,具体为:
s1:对图像采集装置,如摄像机,获取的光纤横截面图像进行预处理:光纤横截面包括涂覆层、包层和纤芯;预处理具体为:利用数学函数对原始图像的灰度值进行图像灰度变换以增强图像对比度。运用空间域的中值滤波,消除图像噪声。然后运用一阶微分算子sobel对图像进行锐化处理以增强图像边缘。如图2所示为摄像机采集到的光纤横截面结构图,图3为光纤横截面图像预处理后的图像。
s2:确定分割阈值,对预处理后的光纤横截面图像进行阈值分割,提取涂覆层的外边缘图像;在本实施例中,运用最大类间方差法确定分割阈值,原理为将:将图像划分为暗背景c1和亮物体c2,确定分割阀值,之后运用形态学及开、闭运算消除图像区域的毛刺,填补区域孔洞,即可直接分割出光纤亮物体区域,从而获得涂覆层的外边缘图像。最大类间方差法具体步骤为:
s21:运用累加法计算像素的灰度值属于暗背景的概率p1和像素的灰度值属于亮物体的概率p2;其中,
pi为图像灰度级i出现的概率,k为图像分割阈值,l为图像的灰度级;
s22:根据p1、p2计算出暗背景区域的平均灰度值g1和亮物体区域的平均灰度值g2;其中,
s23:根据平均灰度值,可计算出图像的类间方差
计算阈值范围内[0,l-1]的每个阈值k的类间方差,以最大的类间方差对应的阈值k作为分割阈值。
s3:分别运用改进的canny算法对涂覆层的内边缘和包层的外边缘进行边缘检测以提取涂覆层的内边缘图像和包层的外边缘图像;图5和图6的对比可知,改进的canny算法可有效减少干扰边缘。其中,提取涂覆层的内边缘图像具体步骤为:
s31:运用形态学方法,由涂覆层的外边缘图像确定涂覆层的内边缘roi区域;
s32:对涂覆层的内边缘roi区域采用二维高斯函数g(x,y)对该区域的图像f(x,y)进行平滑滤波以去除噪声,得到图像fs(x,y);
fs(x,y)=g(x,y)*f(x,y)(7)
其中,σ为控制平滑程度的高斯滤波器参数,σ较小的滤波器定位精度高,但信噪比低;σ较大,则反之。
s33:利用一阶微分sobel算子计算平滑后的图像fs(x,y)各点处的梯度幅值和梯度方向,从而获得相应的梯度幅值图像m(x,y)和梯度方向图像θ(x,y);
其中
s34:在m(x,y)中以点(x,y)为中心,在设置的邻域,如在3×3的邻域内,沿θ(x,y)进行插值,如果点(x,y)处的梯度幅值m(x,y)大于θ(x,y)方向上与其相邻的2个插值,则将点(x,y)标记为候选边缘点,反之则标记为非边缘点,候选边缘点组成候选边缘图像n;
s35:选取阈值tl,检测候选边缘图像n中标记为候选边缘点的任一像素点(x,y);若点(x,y)处的梯度幅值m(x,y)大于其低阈值tl,则将该点标记为边缘点,以获得阈值检测图像n1;选取阈值tl的具体过程为:根据检测的梯度大小,计算出梯度图的统计直方图,再根据直方图求出图像背景,从而获得低阈值,此低阈值即为tl。
s36:对阈值检测图像n1,以中心坐标分成多个区域si;运用局部区域距离阈值法分别对si分别筛选,以获取每个si中符合设置条件的边缘子图像;其中,i为区域个数;具体的,对阈值检测图像n1,运用局部区域距离阈值法进行二次筛选。利用图像的中心坐标把图像分成四个区域si(i=1,2,3,4)。即把图像边缘分配到四个不同的区域,每个边缘坐标为(xi,yi)。每个区域的边缘个数mi,利用边缘处灰度联通性较好的原理,选出si区域中最长边缘作为理想对象,对si区域的其它边缘进行筛选。由于光纤涂覆层为近圆轮廓,运用点到边缘的距离值筛除非轮廓区域。公式如下:
其中di为在si区域中图像的中心坐标到各个边缘的距离,dh为圆心到最长边缘的距离,δ0为设定的微小量。
s37:连接所述边缘子图像以获取涂覆层的外边缘图像。
s4:运用圆形拟合发和椭圆拟合法,对涂覆层的外边缘图像、涂覆层的内边缘图像及包层的外边缘图像均进行圆形拟合和椭圆拟合,以获取各图像对应的圆形参数和椭圆参数;圆拟合法和椭圆拟合方法相同。其中,椭圆拟合法具体为:
s41:运用二次曲线标准化法对提取的边缘图像进行椭圆拟合以获取二次曲线参数,由二次曲线参数求取椭圆中间参数;二次曲线标准化法具体为:
运用二次曲线标准化法对提取的边缘图像进行椭圆拟合,间接求取椭圆中间参数。以(xi,yi)=xi(i=1,2,...n),表示在边缘轮廓上观测的n个点的坐标。设二次曲线方程可表示为:
a1x2+a2xy+a3y2+a4x+a5y+a6=0(11)
其中a1、a2、a3、a4、a5、a6为待测参数,经过变化使二次曲线方程常数项为1,二次曲线方程可表示为:
b1x2+b2xy+b3y2+b4x+b5y+1=0(12)
其中bt=(b1、b2、b3、b4、b5、1)为待测参数。待测参数等于参数近似值与参数改正数之和。二次曲线方程可表示为矩阵形式:
c′·b′+c·δ=0(13)
其中参数近似值为b′t=(b′1、b′2、b′3、b′4、b′5、1)、参数改正数δt=(δ1、δ2、δ3、δ4、δ5)、系数矩阵c′=(x2、xy、y2、x、y、1)、系数矩阵c=(x2、xy、y2、x、y)
每一对坐标点列出误差方程:
vi=c′i·b′+ci·δ(14)
通过最小化误差:
为求使上式之和最小的参数,令
s·δ=s′·b′(16)
其中,
使b′t=(0、0、0、0、0、1)作为初始迭代向量,迭代至参数改正数收敛,最终求得二次曲线参数和二次曲线方程。根据二次曲线方程,求取椭圆中间参数。
s42:采用齐次方程法以步骤s41中获取的椭圆中间参数作为齐次方程法中初始迭代值,求解椭圆参数值。齐次方程法具体为:
椭圆的另一种更直观的方法是用平面坐标系的几何参数表示。即椭圆中心坐标x0=(x0,y0),长轴和短轴a、b,及长轴的转角θ。设参数β=(x0、y0、a、b、θ)。测量坐标系经过旋转,平移变成椭圆标准坐标系。已知(xi,yi)=xi(i=1,2,...n)为观测坐标,坐标经过旋转,平移变为椭圆标准坐标系中的坐标(xs,ys)=xs。它们之间的关系如下:
x=r(α)xs+x0(17)
其中
经过转换之后,椭圆在椭圆标准坐标系中的方程表示为:
其中
列出误差方程:
需要用高斯-牛顿迭代法对参数迭代求解。高斯牛顿迭代法的初始值选择二次曲线拟合所求参数。迭代公式:
其中是函数f((x,y),β)对β的jacobian矩阵
jacobian矩阵
s5:分别将步骤s4中的圆形参数和椭圆参数作为输入参数,运用光纤几何参数求解模型,分别求解步骤s4中各图像对的直径和不圆度。在本实施例中,光纤几何参数求解模型为高斯-牛顿迭代法模型。拟合结果如图7和图8所示。
上述仅为本发明的优选实施例而已,并不对本发明起到任何限制作用。任何所属技术领域的技术人员,在不脱离本发明的技术方案的范围内,对本发明揭露的技术方案和技术内容做任何形式的等同替换或修改等变动,均属未脱离本发明的技术方案的内容,仍属于本发明的保护范围之内。