具体实施方式
以下结合附图对本发明的技术方案作进一步描述。
参照图1,本具体实施方式中液晶显示器缺陷检测系统由CCD摄像机2、图像采集卡5、计算机6、三轴精密定位平台1、液晶显示器驱动模块3和载物台4组成。载物台4固连于三轴精密定位平台1上,并通过电缆线连接到专用的显示器驱动模块3的接口上,被测液晶显示器7放在可以在竖直方向上180度旋转的载物台4上,CCD摄像机2作为图像采集设备固定在三轴精密定位平台1上,使其可以根据被测液晶显示器7的型号在计算机6的控制下移动到相应的位置进行图像采集。
图像采集卡5是标准模拟或数字图像采集卡,也可以是USB接口的数字摄像机或者是视频编辑卡。计算机6是具有标准计算机功能的主控设备,含机箱、CPU及其主板、内存、硬盘和标准计算机兼容的显示器。CCD摄像机2连接到图像采集卡5的输入端;系统软件安装于执行自动检测任务的计算机6的程序存储器,如硬盘中;计算机6控制显示器驱动模块3点亮被测液晶显示器7,并令其显示初始画面;
检测斑痕缺陷时,被测液晶显示器7的载物台4在计算机6的控制下旋转一定的角度,令CCD摄像机2的光轴与被测液晶显示器7平面成30~60度的夹角,使得CCD摄像机2可以从被测液晶显示器7的侧面利用重复采样的方式采集被测液晶显示器的整体图像信息。
由于被测液晶显示器由液晶显示器驱动模块3驱动,因此自身可以发光,无需照明器件。为防止外界光线的干扰,优选的,被测液晶显示器7、载物台4、CCD摄像机2、三轴精密定位平台1都放在一个遮光罩中。计算机6控制液晶显示器驱动模块3令被测液晶显示器7显示为灰度画面。
参照图2,通过此系统可以获得反映液晶显示器基板上像素点的待测液晶显示屏基板的电压信号;而后,经历图像采集步骤801,在此步骤中,由于斑痕缺陷的对比度很低,为减小图像噪声的影响、提高图像质量,在图像采集步骤中采用了对一幅被测画面重复采样的方法。为实现实时检测,采样频率一般应大于10帧/秒,重复采样的次数一般取50~70次。
经采集的图像继而经历图像滤波步骤802。如果叠加于图像上的噪声n(x,y)是非相关的、具有零均值的随机噪声,则可以用几张在相同条件下获得的这种随机图像的平均值表示原图像。假设原图像为f(x,y),加噪声为n(x,y),则有噪声图像为g(x,y)=f(x,y)+n(x,y)
此时,可以用
来估计原图像f(x,y),其中M为图像数量。
显然这种估计是无偏的,因为
完成图像滤波后,还要将被测液晶屏部分从图像中提取出来,并且采用双线性插值法对图像进行几何校正,从而将图像中由于CCD视角的原因而变为梯形的被测液晶屏恢复为标准的矩形。
根据双线性插值法,校正图像中像素点(i′j′)处得灰度值f(i′,j′)可以由其周围四个像素的灰度值计算出来。
f(i′,j′)=[f(i+1,j)-f(i,j)]*(i′-i)+[f(i,j+1)-f(i,j)]*(j′-j)+f(i,j)
+[f(i+1,j+1)+f(i,j)-f(i+1,j)-f(i,j+1)]*(i′-i)(j′-j)
经滤波的图像,进行斑痕缺陷检测时,由于液晶显示器本身的特点和CCD视角的关系,被测液晶屏的灰度在整个屏幕范围内变化非常大,远远超过了斑痕缺陷本身的灰度变化。因此,对于斑痕缺陷的检测来说,图像分割是最为重要同时也是最为困难的一个环节。运用传统的图像分割方法,例如自适应阈值分割,边缘检测等等,根本无法完成分割任务。为此,提出了一种基于多项式曲面拟合技术的斑痕图像分割方法。
在这种方法中,首先假设图像中每个像素的灰度值f(x,y)都是该像素二维坐标(x,y)的函数,而所有像素的灰度值及其二维坐标构成了分布在矩形格点上的空间数据点的集合。接下来采用二元多项式对这些空间数据点进行曲面拟合803,即求一个近似通过给定数据点、能反映数据点的基本趋势的多项式曲面:
设已知矩形区域内n×m个网点(xk,yl)(k=0,1,...,n-1;l=0,1,...,m-1)上的函数值zkl,设待求的拟合多项式为
其中aij(i=0,1,..,p-1)(j=0,1,...,q-1)是待定参数,它使zkl,f(x,y)在网格点上的值的差的平方和在最小二乘意义下达到最小,即
由上式所决定的参数aij,其对应的曲面称为最小二乘多项式曲面。
利用多元函数的极值理论,在上式中令
便得到关于{aij}的p×q阶代数方程组:
解这个p×q阶的代数方程组,就可以得到aij。
由于斑痕缺陷的面积很小,不会改变图像的基本变化趋势,因此可以将拟合曲面看成是不含斑痕缺陷的背景模型。相对于该背景模型,被测液晶屏图像中没有缺陷的区域拟合误差
很小,而斑痕缺陷所在的区域与背景模型的拟合误差
则很大。所以,利用该背景模型,采用背景减除的方法就可以将那些拟合误差超过人工设定的阈值(例如2)的像素作为斑痕缺陷的疑似目标与背景分割开来804。
图3A为CCD采集的原始图像经过多图像平均法和图像校正处理后的结果。图中有五处斑痕缺陷。图3B是对图3A进行多项式曲面拟合法所获得的不含斑痕缺陷的背景模型。图3C为利用图3B所示的背景模型对图3A进行背景对消所获得的处理结果。
图像后处理步骤。本发明在图像后处理模块中采用了数学形态学的处理方法。膨胀、腐蚀是数学形态学的两个基本运算。采用这些基本运算及其组合可以进行图像形状和结构的分析及处理。
设集合A为输入图像,集合B为结构元素,那么集合A被集合B腐蚀,就表示为:
AΘB={x:B+xA}
腐蚀的作用是滤除图像中的一些孤立的噪声点和毛刺。膨胀是腐蚀运算的对偶运算(逆运算),可以通过对补集的腐蚀来定义。A被集合B膨胀表示为AB,其定义为:
AB=[AcΘ(-B)]c
其中,Ac表示A的补集。
膨胀的作用是可以填充图像中比结构元素小的孔洞和图像边缘处的小凹陷部分。
图3D是对图3C进行形态学处理的结果。
特征提取步骤805:在特征提取中,本发明从人类的视觉心理学出发,提出将目标区域的对比度、面积、边缘参数、位置参数、灰度均匀性和形状参数等六个特征量作为斑痕缺陷等级评价的依据。
其中,对比度参数定义为:
式中f(i,j)和B(i,j)分别为疑似目标区域及背景模型在像素点(i,j)处的灰度值;U为目标区域内所有像素点的集合;N为目标区域内像素点的个数。
面积参数定义为:
其中,U为目标区域内所有像素点的集合。
位置参数定义为目标中心与被测TFT-LCD屏幕中心的距离,即
其中x0和y0分别为被测TFT-LCD屏几何中心的横、纵坐标;
x和
y分别为疑似目标区域几何中心的横、纵坐标,且有:
式中xij和yij分别为疑似目标中(i,j)像素点的横坐标、纵坐标。
边缘参数是用拉普拉斯算定义的:
形状参数定义为:
shape=max{R,C}
式中,R为目标区域的圆形度,描述了物体边界的复杂程度,可以用目标的面积与周长的关系来表示。圆形度的数学表达式为:
式中S为目标区域的面积,P为目标区域的周长。
C为目标区域的矩形度,反映物体对其外接矩形的充满程度,用目标区域的面积与其最小的外接矩形面积之比来描述,即:
其中,S为目标区域的面积,SMER为目标区域最小外接矩形的面积。
灰度均匀性用目标区域的灰度标准差来加以描述:
式中,
为目标区域内所有像素点的灰度均值。
模式识别步骤:在模式识别中,本发明提出了基于模糊集合理论和模糊逻辑的斑痕缺陷模糊识别方法。
在设计模糊模式识别模块时,首先要定义系统的输入、输出变量。本发明将疑似目标区域的对比度、面积、边缘参数、位置参数、灰度均匀性和形状参数等6个特征量作为系统的输入变量;斑痕缺陷的等级作为系统的输出变量。其中斑痕缺陷的等级一般划分为不合格(N级)、合格(P级)、良(V级),实际应用也可以根据工厂的实际情况适当的加以调整。
接下来,根据专业技术人员的经验将系统的输入、输出变量划分成一定数目的模糊子集并为每个模糊子集建立相应的隶属函数。为了简化运算、提高系统的运行速度,一般应选取形状相对简单的三角形和梯形的隶属度函数。
例如,对于系统的6个输入变量,将其作了如下的划分:将比度分成五级,即极低(VL)、低(L)、中等(M)、高(H)和极高(VH);面积分成小(S)、中(M)和大(B)三级;边缘参数分成低(L)、中(M)、高(H);位置参数分成靠边(A)、一般(N)和居中(C)三级;灰度均匀性分成低(L)、中等(M)和高(H)三级;形状参数分成不规则(NR)和规则(R)两级。系统的输出变量即mura缺陷等级划分为轻微mura缺陷(P)、一般mura缺陷(N)和严重mura缺陷(V)三级。输入、输出变量所对应的隶属度函数都选用三角函数和梯形函数来表示,如图3所示。
模糊if-then规则的制订是模糊模式识别的核心问题,本发明中模糊if-then规则同样也是根据专家经验制定的。
例如,一条if-then规则如下:
If(Const is H)and(Area is M)and(Edge is H)then(mura_level is V)
采用模糊识别方法,系统可以模仿人的辩识方式并充分利用专家的经验和知识,完成对液晶屏斑痕缺陷的自动识别。