CN105976382B - 一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 - Google Patents
一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 Download PDFInfo
- Publication number
- CN105976382B CN105976382B CN201610310495.7A CN201610310495A CN105976382B CN 105976382 B CN105976382 B CN 105976382B CN 201610310495 A CN201610310495 A CN 201610310495A CN 105976382 B CN105976382 B CN 105976382B
- Authority
- CN
- China
- Prior art keywords
- value
- background
- defect area
- pixel
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于缺陷区域预判和水平集的TFT‑LCD Mura缺陷检测方法。该方法提出了一种缺陷区域预判法先找到缺陷区域,然后将缺陷区域的像素剔除,用余下的像素拟合得到背景图像;用原图像减去背景图像来消除背景不均匀性对缺陷分割的影响,求差后得到残差图像,然后用一种基于阈值的水平集方法对残差图像进行分割来检测得到缺陷。该方法可以获得精度较高的背景图像,对于光照影响有较好的鲁棒性,能够获得准确的分割结果以及较低的误检率。
Description
技术领域
本发明属于图像处理技术领域,更具体地,涉及一种基于缺陷区域预判和水平集的TFT-LCD Mura缺陷检测方法。
背景技术
薄膜晶体管液晶显示器(thin film transistor-liquid crystal display)简称TFT-LCD,它对于像素点的控制精度高,反应速度快,具有高彩色保真度,高亮度,高对比度同时所需功耗低,在低温和高温条件下都能使用;正是由于TFT-LCD具有上述诸多优点,使得其被广泛应用到我们日常生活中,如手机,液晶电视,笔记本电脑等。
根据有关统计显示,我国TFT-LCD的产量远低于市场需求,其中一个主要原因是在TFT-LCD生产过程中有些工序还没有实现自动化,比如对于屏幕的检测工作,现在许多大型的液晶屏生产商还是选用人工来进行检测,随着人力成本的不断提高,大量使用人工进行检测不仅限制了产量而且使得生产成本提高,降低了产品市场的竞争力。所以市场急需一种能够检测缺陷的设备来实现缺陷检测的自动化。
基于图像的检测技术可以代替人眼完成缺陷的检测工作,对于 TFT-LCD屏幕其中Mura缺陷最难检测,主要由于Mura缺陷与背景对比度低,没有明显的边界以至于很难和背景分割开,如图1所示,这也是基于视觉技术实现TFT-LCD屏幕检测的难点。
由于图像会受到光照等影响,背景也会不均匀,这使得Mura缺陷更难被分割出来,所以视觉检测方法需要解决在这些环境干扰下Mura分割问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种基于缺陷区域预判和水平集的TFT-LCD Mura缺陷检测方法,其目的在于获得精度较高的背景图像,对于光照影响有较好的鲁棒性,能够获得准确的分割结果以及较低的误检率。
为实现上述目的,按照本发明的一个方面,提供了一种一种基于缺陷区域预判和水平集的TFT-LCD Mura缺陷检测方法,其特征在于,该方法包括如下步骤:
(1)缺陷区域预判
(1-1)缺陷区域预处理
对采集的图像每行每列数据L,L为一维数组,长度为M,首先对 L进行中值滤波消除噪声的影响,滤波后得到Ls,对Ls进行一维三次曲线拟合,得到初步的拟合背景图像,拟合公式如下:
y=α0+α1x+α2x2+α3x3
其中,x表示自变量,y表示应变量,α0,α1,α2,α3都为一维三次曲线拟合公式的参数;
用矩阵表达如下:
上述矩阵简化为XA=Y,其中,矩阵简化为A,矩阵简化为X,矩阵简化为Y,A是未知的,X矩阵包含的是像素点的位置坐标,Y矩阵包含的对应坐标位置的像素点的值, X矩阵中x1,...,xM表示像素点的位置坐标,Y矩阵中y1,...,yM表示对应坐标位置的像素点的值;A=(XTX)-1XTY;获得参数矩阵后,得到初步的拟合背景曲线Lf=XA;
获得拟合的背景曲线Lf后,Lf与Ls之间形成交点{P1,P2,P2,...,PN},这些交点会将曲线分成N+1块,交点将所述背景曲线分为分块的曲线;并将标记出的缺陷区域用三角形取代;
(1-2)基于SEMU的缺陷区域预判
按照修正的SEMU,计算公式如下:
其中|Cx|是Mura缺陷的平均对比度,是Mura缺陷区域的像素的均值,是背景区域像素的均值,Cjnd是Mura缺陷的JND值,JND是指人眼所能够观测到的最小的对比度差值,S是Mura缺陷的面积,单位为像素;
计算|Cx|,其计算公式如下:
其中,|Cassume|为基于假设背景计算得到的|Cx|值;h为缺陷区域最大值和实际背景之间的差值;LB为两交点的连线,LB(x)表示连线LB上的一点,即连线LB上横坐标为x的点;l是两交点之间连线的长度,是假设背景的均值;
基于实际背景计算的|Cx|值,根据建立的三角模型,推导出其公式如下:
其中|Cactual|为基于实际背景计算得到的|Cx|值,为实际背景的均值;综上得到基于实际背景和假设背景之间的关系:
由于|Cassume|基于已知的交点得到,根据两者之间的关系,通过|Cassume|得到基于实际背景的|Cx|值;
对于上述SEMU公式中S的大小用(2l)2替代,缺陷区域的SEMU的计算公式进一步计算为:
对于没有缺陷的区域,其计算公式如下:
如果每块的最大值或者最小值和两个端点的差值的绝对值大于2,则用公式(1)计算该块的SEMU的值,否则选用公式(2)计算SEMU的值;每一块的SEMU的值计算完毕后,开始检查每一块的SEMU的大小,如果第 n块对应的SEMU的值大于或等于1,则将该块的标记Llabel(x)置为1,即Llabel(x)=1,x∈[Pn-1,Pn];
(1-3)基于斜率区域搜索方法
在检测到缺陷位置之后还需要向两侧进行搜索,找到实际缺陷区域;
当计算得到第n块的斜率大于1时,找到第n块的两个端点Pn-1和Pn,先从Pn-1开始向左侧进行搜索,搜索方法fsearch_left步骤如下:
第一步,设置搜索的步长为d=10,这里已知左侧端点Pn-1处的斜率为 kleft_start=(Ls(Pn-1+d)-Ls(Pn-1))/d,初始化搜索点的位置xl=Pn-1,搜索点的斜率为:kl=(Ls(xl)-Ls(xl-d))/d;如果kl/kleft_start>0.8,搜索开始, xl=xl-d;否则搜索结束;
第二步,Llabel(x)=1,x∈[xl,xl+d],xl=xl-d, kl=(Ls(xl)-Ls(xl-d))/d;
第三步,如果kl/kleft_start>0.5并且xl-d>P1,返回第二步;否则执行下一步;
第四步,如果xl-d<P1,执行Llabel(x)=1,x∈[P1,xl+d],然后搜索结束;否则,搜索直接结束;
确定右侧端点位置Pn后,从Pn开始向右侧搜索,向右搜索fsearch_right步骤如下:
第一步,设置搜索的步长为d=10,这里已知右侧端点Pn处的斜率为 kright_start=(Ls(Pn-d)-Ls(Pn))/d,初始化搜索点的位置xr=Pn,搜索点的斜率为:kr=(Ls(xr+d)-Ls(xr))/d;如果kr/kright_start>0.8,搜索开始, xr=xr+d;否则搜索结束;
第二步,Llabel(x)=1,x∈[xr-d,xr],xr=xr+d,kl=(Ls(xl)-Ls(xl-d))/d;
第三步,如果kr/kright_start>0.5并且xr+d<PN,返回第二步;否则执行下一步;
第四步,如果xr+d>PN,执行Llabel(x)=1,x∈[xr-d,PN],然后搜索结束;否则,搜索直接结束;
(1-4)基于斜率的缺陷区域预判
首先在没有被标记为1的块中,对比各个块的最大斜率值Smax(n),找到比其他块的最大斜率值都大的分块,如果该块的SEMU的值大于0.1,则将该块标记为1,然后根据该块的左右两个端点分别调用左侧搜索方法和右侧搜索方法;
(2)图像背景重构
逐行逐列处理完之后分别得到行处理结果和列处理结果,对行列处理后的结果采用逻辑或操作将结果组合在一起,其中白色区域为缺陷区域,将原图对应缺陷区域位置的像素剔除掉,用余下的像素通过二维三次多项式进行背景图像重构,计算公式如下:
Y=β0+β1x+β2y+β3x2+β4y2+β5xy+β6x2y+β7xy2+β8x3+β9y3
其中,β0,...,β9为待拟合的二维三次多项式的参数;
计算公式简化为ZA=Y,每个矩阵的表达如下:
A=[β0 β1 β2 β3 β4 β5 β6 β7 β8 β9]T
Y=[Y1 Y2 Y3 … YK-1 YK]T
其中K为剔除缺陷区域像素值后剩余的像素值的个数,Z为包含像素位置信息的矩阵,x1,…,xK为像素位置信息的x坐标,y1,...,yK为像素位置信息的 y坐标;Y包含对应像素位置像素值大小的矩阵,Y1,...,YK为对应像素位置的像素值;A是要求解的系数矩阵,A的计算为A=(ZTZ)-1ZTY,获得系数矩阵后,拟合得到背景图像Ibackground=ZA;
(3)水平集分割
在获得背景图像后,通过利用原图像和背景图像求差来消除原图像的背景不均匀的情况,然后对求差后的图像采用改进后的主动轮廓模型进行分割,改进后主动轮廓模型的水平集表达公式如下:
在已有的模型中加入了项,其中ρ=1,G是二维高斯函数,I是输入的图像,I(x,y)表示坐标(x,y) 处输入图像的像素值,阈值T=0.6(gmax+gmin),gmax和gmin分别是g中的最大值和最小值,ω为曲线内部,Ω\ω为曲线外部,C为曲线上一点,Ω表示整个图像域;
该公式隐函数φ(x,y)的偏微分方程为:
其中div(·)表示散度计算,λ1=λ2=1,ν=0.3,μ=0.001×2552,η=30,初始化φ(x,y),对该水平集表达公式的偏微分方程进行求解获得φn(x,y),再将φn(x,y)代回偏微分公式求解得到φn+1(x,y),如此迭代,获得缺陷分割结果。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
提出了一种缺陷区域预判法先找到缺陷区域,然后将缺陷区域的像素剔除,用余下的像素拟合得到背景图像;用原图像减去背景图像来消除背景不均匀性对缺陷分割的影响,求差后得到残差图像,然后用一种基于阈值的水平集方法对残差图像进行分割来检测得到缺陷。该方法可以获得精度较高的背景图像,对于光照影响有较好的鲁棒性,能够获得准确的分割结果以及较低的误检率。
附图说明
图1是本发明实现的方法所针对的Mura缺陷示意图,图1(a)是圆状Mura,图1(b)是线状Mura;
图2是本发明实现的方法流程图,图2(a)是本发明的整个方法的流程图,图2(b)本发明的方法中缺陷区域预判方法的流程图;
图3是本发明的基于SEMU缺陷区域预判法说明图;图3(a)是本发明的基于SEMU缺陷区域预判法所针对的非边缘缺陷的示意图,图3(b) 是本发明基于SEMU缺陷区域预判法求取拟合曲线以及曲线交点的示意图,图3(c)是本发明的基于SEMU缺陷区域预判法以曲线交点连线作为背景的示意图,图3(d)是基于SEMU缺陷区域预判法中建模示意图
图4是本发明的基于斜率的缺陷区域预判法说明图;图4(a)是本发明的基于斜率的缺陷区域预判法所针对的边缘缺陷的示意图,图4(b)是边缘缺陷对于交点位置造成影响的示意图
图5是本发明的行列处理及缺陷区域融合示意图;图5(a)是本发明处理图片示例图,图5(b)是本发明的缺陷区域预判方法行处理的结果示意图,图5(c)是本发明的缺陷区域预判方法列处理的结果示意图,图5 (d)本发明的缺陷区域预判方法行处理的结果和列处理结果融合的示意图。
图6是本发明的原图与背景图像求差示意图;图6(a)是本发明选用作为说明的待处理图片,图6(b)是本发明处理得到的背景图,图6(c) 是本发明选用的待处理图片与本发明处理得到的背景图片求差后的结果图。
图7是本发明的基于水平集的缺陷分割示意图;图7(a)是本发明的基于水平集的缺陷分割方法的初始轮廓曲线位置示意图,图7(b)是本发明的基于水平集的缺陷分割的结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明提供了一种基于图像背景重构和水平集的TFT-LCD Mura缺陷检测方法,用于提高对Mura缺陷的识别率
本发明的Mura缺陷检测方法分为缺陷区域预判,背景图像重构和背景影响去除及基于水平集的缺陷分割三个阶段,整个方法流程如图2(a) 所示。以下对本发明提出的识别方法作详细说明。具体步骤如下:
(1)缺陷区域预判
这个阶段分为两步进行,第一步是基于SEMU的缺陷区域预判,第二步是基于斜率的缺陷区域预判,整个流程图如图2(b)所示,详细方法如下。
(1.1)缺陷区域预处理
对图像每行每列分别进行处理,找出每行每列中缺陷区域的位置并标记出来。这里用L来表示需要进行处理的每行或者每列的数据,L为一维数组,长度为M。首先对L进行中值滤波消除噪声的影响,滤波后得到 Ls。然后对Ls进行一维三次曲线拟合,得到初步的拟合背景图像,拟合公式如下:
y=α0+α1x+α2x2+α3x3
用矩阵表达如下:
上述矩阵可以简化为XA=Y,其中A是未知的,X矩阵包含的是像素点的位置坐标,Y矩阵包含的对应坐标位置的像素点的值,A=(XTX)-1XTY。获得参数矩阵后,可以得到初步的拟合背景曲线Lf=XA。
获得拟合的背景曲线Lf后,Lf与Ls之间形成交点{P1,P2,P2,...,PN},这些交点会将曲线分成N+1块。图3(a)为Mura图像,虚线框出的为被处理的行数据,图3(b)中黑的曲线为平滑后的曲线Ls,虚线曲线为拟合得到的背景曲线Lf,黑色的点为两条线的交点。B1,B2,...,B5为分块的曲线。由于 Mura缺陷变化比较缓慢,所以缺陷边缘的变化可以近似看成一条直线,可以将缺陷区域用三角形进行取代,模型如图3(d)所示。
(1.2)基于SEMU的缺陷区域预判
引入修改的SEMU判断公式来判断缺陷位置,修改的SEMU公式如下:
其中|Cx|是Mura缺陷的平均对比度,是Mura缺陷区域的像素的均值,是背景区域像素的均值,Cjnd是Mura缺陷的JND(Just noticeable difference)值,JND是指人眼所能够观测到的最小的对比度差值,S是 Mura缺陷的面积,单位为像素。
首先基于假设背景(图3(d)两交点的虚线连线)计算|Cx|,其计算公式如下:
其中h为缺陷区域最大值和实际背景之间的差值,LB为两交点的连线,l是两交点之间连线的长度,是假设背景的均值,l1+l2=l,h/2l1是三角形左侧边的斜率,h/2l2为三角形右侧边的斜率。
基于实际背景(图3(d)三角形的底边)计算的|Cx|值,根据建立的模型,可以推导出其公式如下:
其中|Cactual|为基于实际背景计算得到的|Cx|值,为实际背景的均值。综上可以得到基于实际背景和假设背景之间的关系:
由于|Cassume|可以基于已知的交点得到,根据两者之间的关系,可以通过|Cassume|得到基于实际背景的|Cx|值。
对于SEMU公式中S的大小,这里用(2l)2替代,所以缺陷区域的SEMU的计算公式为:
对于没有缺陷的区域,如图4(b)所示,第二块区域B2为非缺陷区域,P1和 P2两者之间的连线就是背景,所以对非缺陷区域不需要乘以r,其计算公式如下:
如果每块的最大值或者最小值和两个端点的差值的绝对值大于2,则用公式计算该块的SEMU的值,否则选用公式计算SEMU的值。
每一块的SEMU的值计算完毕后,开始检查每一块的SEMU的大小,如果第n块对应的SEMU的值大于或等于1,则将该块的标记置为1,即 Llabel(x)=1,x∈[Pn-1,Pn]。
(1.3)基于斜率区域搜索方法
从图3(b)中我们可以看出,虽然B3为缺陷区域,B3要比虚线框框出的实际的缺陷区域小,因此需要在检测到缺陷位置之后还需要向两侧进行搜索,找到实际缺陷区域。对于这个问题本文提出了基于斜率的搜索方法。
当计算得到第n块的斜率大于1时,找到第n块的两个端点Pn-1和Pn,先从Pn-1开始向左侧进行搜索,搜索方法fsearch_left步骤如下:
第一步,设置搜索的步长为d=10,这里已知左侧端点Pn-1处的斜率为 kleft_start=(Ls(Pn-1+d)-Ls(Pn-1))/d,初始化搜索点的位置xl=Pn-1,搜索点的斜率为:kl=(Ls(xl)-Ls(xl-d))/d。如果kl/kleft_start>0.8,搜索开始, xl=xl-d。否则搜索结束。
第二步,Llabel(x)=1,x∈[xl,xl+d],xl=xl-d,kl=(Ls(xl)-Ls(xl-d))/d。
第三步,如果kl/kleft_start>0.5并且xl-d>P1,返回第二步;否则执行下一步。
第四步,如果xl-d<P1,执行Llabel(x)=1,x∈[P1,xl+d],然后搜索结束。否则,搜索直接结束。
确定右侧端点位置Pn后,从Pn开始向右侧搜索,向右搜索fsearch_right详细步骤如下:
第一步,设置搜索的步长为d=10,这里已知右侧端点Pn处的斜率为 kright_start=(Ls(Pn-d)-Ls(Pn))/d,初始化搜索点的位置xr=Pn,搜索点的斜率为:kr=(Ls(xr+d)-Ls(xr))/d。如果kr/kright_start>0.8,搜索开始, xr=xr+d。否则搜索结束。
第二步,Llabel(x)=1,x∈[xr-d,xr],xr=xr+d,kl=(Ls(xl)-Ls(xl-d))/d。
第三步,如果kr/kright_start>0.5并且xr+d<PN,返回第二步;否则执行下一步。
第四步,如果xr+d>PN,执行Llabel(x)=1,x∈[xr-d,PN],然后搜索结束。否则,搜索直接结束。
(1.4)基于斜率的缺陷区域预判
对于处于边界上的Mura缺陷,如图4(a)所示为重力Mura,在缺陷区域,交点的位置会因为缺陷处于边界而受到影响,如图4(b)所示,因为处在边界,P4为实际的交点位置,但是理想的交点位置在P′4处,黑的点画线为实际曲线的延长线,细虚线是理想拟合曲线。由于实际交点位置靠近Ls,这就会使得计算SEMU时得到的值偏大,从而使得得到的计算结果偏小,甚至小于1,计算结果的偏差会使得该缺陷区域不能被判断出来,从而影响到后续的背景图像拟合的精度。为了避免这种情况的发生,本文针对缺陷处于边缘位置的问题提出了一种基于斜率的区域预判法。
通过观察曲线可以得出,虽然Mura缺陷区域里像素值之间的变化比较平滑,但是相对于背景而言,Mura缺陷位置的曲线斜率要大于背景区域曲线斜率。这就是该方法的主要出发点,找到斜率最大的位置,具体的处理步骤如下:
首先在没有被标记为1的块中,对比各个块的最大斜率值Smax(n),找到比其他块的最大斜率值都大的分块,如果该块的SEMU的值大于0.1,则将该块标记为1,然后根据该块的左右两个端点分别调用左侧搜索方法和右侧搜索方法。
(2)图像背景重构
逐行逐列处理完之后分别得到行处理结果(如图5(b))和列处理结果(图5(c)),对行列处理后的结果采用逻辑或操作将结果组合在一起 (如图5(d),其中白色区域为缺陷区域),将原图对应缺陷区域位置的像素剔除掉,用余下的像素通过二维三次多项式进行背景图像重构,计算公式如下:
Y=β0+β1x+β2y+β3x2+β4y2+β5xy+β6x2y+β7xy2+β8x3+β9y3
计算公式可以简化为ZA=Y,每个矩阵的表达如下:
A=[β0 β1 β2 β3 β4 β5 β6 β7 β8 β9]T
Y=[Y1 Y2 Y3 … YK-1 YK]T
其中K为剔除缺陷区域像素值后剩余的像素值的个数,Z为包含像素位置信息的矩阵,Y包含对应像素位置像素值大小的矩阵,A是要求解的系数矩阵,A的计算为A=(ZTZ)- 1ZTY,获得系数矩阵后,可以拟合得到背景图像Ibackground=ZA。
(2)基于水平集的缺陷分割
在获得背景图像后,可以通过利用原图像(如图6(a))和背景图像 (如图6(b))求差来消除原图像的背景不均匀的情况,然后对求差后的图像(如图6(c))采用改进后的主动轮廓模型进行分割,分割结果如图 7所示,图7(a)中的白色曲线为初始化轮廓曲线,图7(b)白色曲线为找到的缺陷边缘位置。改进后主动轮廓模型的水平集表达公式如下:
在已有的模型中加入了项,该项的加入是为了增强改进后的主动轮廓模型对于小缺陷的敏感性,其中ρ=1,G是二维高斯函数,I是输入的图像,阈值 T=0.6(gmax+gmin),gmax和gmin分别是g中的最大值和最小值,ω为曲线内部,Ω\ω为曲线外部,C是曲线上一点。
该公式隐函数φ的偏微分方程为:
其中λ1=λ2=1,ν=0.3,μ=0.001×2552,η=30。初始化φ(x,y),对该水平集方法的偏微分方程进行求解获得φn(x,y),再将φn(x,y)代回偏微分公式求解得到φn+1(x,y),如此迭代,便可以将缺陷分割出来,结果如图7所示。
本发明所提出的方法,利用基于缺陷区域预判的背景重构方法首先找到缺陷区域位置,然后将缺陷区域的像素剔除减少缺陷像素对重构背景图像精度的影响,利用余下的像素重构背景图像,从而提高获得的背景图像的精度。然后用原图像减去背景图像来消除背景不均匀的影响,将Mura缺陷突出出来,最后采用改进的水平集方法对图像进行分割,提高图像分割的精度。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种基于缺陷区域预判和水平集的TFT-LCD Mura缺陷检测方法,其特征在于,该方法包括如下步骤:
(1)缺陷区域预判
(1-1)缺陷区域预处理
对采集的图像每行每列数据L,L为一维数组,长度为M,首先对L进行中值滤波消除噪声的影响,滤波后得到Ls,对Ls进行一维三次曲线拟合,得到初步的拟合背景图像,拟合公式如下:
y=α0+α1x+α2x2+α3x3
其中,x表示自变量,y表示应变量,α0,α1,α2,α3都为一维三次曲线拟合公式的参数;
用矩阵表达如下:
上述矩阵简化为XA=Y,其中,矩阵简化为A,矩阵简化为X,矩阵简化为Y,A是未知的,X矩阵包含的是像素点的位置坐标,Y矩阵包含的对应坐标位置的像素点的值,X矩阵中x1,...,xM表示像素点的位置坐标,Y矩阵中y1,...,yM表示对应坐标位置的像素点的值;A=(XTX)-1XTY;获得参数矩阵后,得到初步的拟合背景曲线Lf=XA;
获得拟合的背景曲线Lf后,Lf与Ls之间形成交点{P1,P2,P2,...,PN},这些交点会将曲线分成N+1块,交点将所述背景曲线分为分块的曲线;并将标记出的缺陷区域用三角形取代;
(1-2)基于SEMU的缺陷区域预判
按照修正的SEMU,计算公式如下:
其中|Cx|是Mura缺陷的平均对比度,是Mura缺陷区域的像素的均值,是背景区域像素的均值,Cjnd是Mura缺陷的JND值,JND是指人眼所能够观测到的最小的对比度差值,S是Mura缺陷的面积,单位为像素;
计算|Cx|,其计算公式如下:
其中,|Cassume|为基于假设背景计算得到的|Cx|值;h为缺陷区域最大值和实际背景之间的差值;LB为两交点的连线,LB(x)表示连线LB上的一点,即连线LB上横坐标为x的点;l是两交点之间连线的长度,是假设背景的均值;
基于实际背景计算的|Cx|值,根据建立的三角模型,推导出其公式如下:
其中|Cactual|为基于实际背景计算得到的|Cx|值,为实际背景的均值;综上得到基于实际背景和假设背景之间的关系:
由于|Cassume|基于已知的交点得到,根据两者之间的关系,通过|Cassume|得到基于实际背景的|Cx|值;
对于上述SEMU公式中S的大小用(2l)2替代,缺陷区域的SEMU的计算公式进一步计算为:
对于没有缺陷的区域,其计算公式如下:
如果每块的最大值或者最小值和两个端点的差值的绝对值大于2,则用公式(1)计算该块的SEMU的值,否则选用公式(2)计算SEMU的值;每一块的SEMU的值计算完毕后,开始检查每一块的SEMU的大小,如果第n块对应的SEMU的值大于或等于1,则将该块的标记Llabel(x)置为1,即Llabel(x)=1,x∈[Pn-1,Pn];
(1-3)基于斜率区域搜索方法
在检测到缺陷位置之后还需要向两侧进行搜索,找到实际缺陷区域;
当计算得到第n块的斜率大于1时,找到第n块的两个端点Pn-1和Pn,先从Pn-1开始向左侧进行搜索,搜索方法fsearch_left步骤如下:
第一步,设置搜索的步长为d=10,这里已知左侧端点Pn-1处的斜率为kleft_start=(Ls(Pn-1+d)-Ls(Pn-1))/d,初始化搜索点的位置xl=Pn-1,搜索点的斜率为:kl=(Ls(xl)-Ls(xl-d))/d;如果kl/kleft_start>0.8,搜索开始,xl=xl-d;否则搜索结束;
第二步,Llabel(x)=1,x∈[xl,xl+d],xl=xl-d,kl=(Ls(xl)-Ls(xl-d))/d;
第三步,如果kl/kleft_start>0.5并且xl-d>P1,返回第二步;否则执行下一步;
第四步,如果xl-d<P1,执行Llabel(x)=1,x∈[P1,xl+d],然后搜索结束;否则,搜索直接结束;
确定右侧端点位置Pn后,从Pn开始向右侧搜索,向右搜索fsearch_right步骤如下:
第一步,设置搜索的步长为d=10,这里已知右侧端点Pn处的斜率为kright_start=(Ls(Pn-d)-Ls(Pn))/d,初始化搜索点的位置xr=Pn,搜索点的斜率为:kr=(Ls(xr+d)-Ls(xr))/d;如果kr/kright_start>0.8,搜索开始,xr=xr+d;否则搜索结束;
第二步,Llabel(x)=1,x∈[xr-d,xr],xr=xr+d,kl=(Ls(xl)-Ls(xl-d))/d;
第三步,如果kr/kright_start>0.5并且xr+d<PN,返回第二步;否则执行下一步;
第四步,如果xr+d>PN,执行Llabel(x)=1,x∈[xr-d,PN],然后搜索结束;否则,搜索直接结束;
(1-4)基于斜率的缺陷区域预判
首先在没有被标记为1的块中,对比各个块的最大斜率值Smax(n),找到比其他块的最大斜率值都大的分块,如果该块的SEMU的值大于0.1,则将该块标记为1,然后根据该块的左右两个端点分别调用左侧搜索方法和右侧搜索方法;
(2)图像背景重构
逐行逐列处理完之后分别得到行处理结果和列处理结果,对行列处理后的结果采用逻辑或操作将结果组合在一起,其中白色区域为缺陷区域,将原图对应缺陷区域位置的像素剔除掉,用余下的像素通过二维三次多项式进行背景图像重构,计算公式如下:
Y=β0+β1x+β2y+β3x2+β4y2+β5xy+β6x2y+β7xy2+β8x3+β9y3
其中,β0,...,β9为待拟合的二维三次多项式的参数;
计算公式简化为ZA=Y,每个矩阵的表达如下:
A=[β0 β1 β2 β3 β4 β5 β6 β7 β8 β9]T
Y=[Y1 Y2 Y3 … YK-1 YK]T
其中K为剔除缺陷区域像素值后剩余的像素值的个数,Z为包含像素位置信息的矩阵,x1,...,xK为像素位置信息的x坐标,y1,...,yK为像素位置信息的y坐标;Y包含对应像素位置像素值大小的矩阵,Y1,...,YK为对应像素位置的像素值;A是要求解的系数矩阵,A的计算为A=(ZTZ)-1ZTY,获得系数矩阵后,拟合得到背景图像Ibackground=ZA;
(3)水平集分割
在获得背景图像后,通过利用原图像和背景图像求差来消除原图像的背景不均匀的情况,然后对求差后的图像采用改进后的主动轮廓模型进行分割,改进后主动轮廓模型的水平集表达公式如下:
在已有的模型中加入了项,其中ρ=1,G是二维高斯函数,I是输入的图像,I(x,y)表示坐标(x,y)处输入图像的像素值,阈值T=0.6(gmax+gmin),gmax和gmin分别是g中的最大值和最小值,ω为曲线内部,Ω\ω为曲线外部,C为曲线上一点,Ω表示整个图像域;
该公式隐函数φ(x,y)的偏微分方程为:
其中div(·)表示散度计算,λ1=λ2=1,ν=0.3,μ=0.001×2552,η=30,初始化φ(x,y),对该水平集表达公式的偏微分方程进行求解获得φn(x,y),再将φn(x,y)代回偏微分公式求解得到φn+1(x,y),如此迭代,获得缺陷分割结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610310495.7A CN105976382B (zh) | 2016-05-11 | 2016-05-11 | 一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610310495.7A CN105976382B (zh) | 2016-05-11 | 2016-05-11 | 一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105976382A CN105976382A (zh) | 2016-09-28 |
CN105976382B true CN105976382B (zh) | 2018-11-13 |
Family
ID=56992325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610310495.7A Active CN105976382B (zh) | 2016-05-11 | 2016-05-11 | 一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105976382B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106650770B (zh) * | 2016-09-29 | 2019-12-17 | 南京大学 | 一种基于样本学习和人眼视觉特性的mura缺陷检测方法 |
WO2018187989A1 (zh) * | 2017-04-12 | 2018-10-18 | 意力(广州)电子科技有限公司 | 一种显示屏的检测方法及装置 |
CN107564007B (zh) * | 2017-08-02 | 2020-09-11 | 中国科学院计算技术研究所 | 融合全局信息的场景分割修正方法与*** |
CN107316303B (zh) * | 2017-08-25 | 2021-11-02 | 深圳市华星光电半导体显示技术有限公司 | 显示面板Mura区域侦测方法 |
CN108280825A (zh) * | 2018-01-19 | 2018-07-13 | 东华大学 | 一种液晶显示屏边缘缺口检测方法 |
CN111553929B (zh) * | 2020-05-12 | 2022-08-09 | 重庆邮电大学 | 基于融合网络的手机屏幕缺陷分割方法、装置及设备 |
CN111768348B (zh) * | 2020-05-29 | 2022-02-08 | 合肥联宝信息技术有限公司 | 一种缺陷检测方法、设备及计算机可读存储介质 |
CN112068064B (zh) * | 2020-09-18 | 2023-08-25 | 国网重庆市电力公司营销服务中心 | 一种扣除背景影响量的运行三相电压互感器误差计算方法 |
CN115797342B (zh) * | 2023-02-06 | 2023-04-21 | 深圳市鑫旭飞科技有限公司 | 一种工业控制电容触摸lcd显示总成缺陷检测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674975A (zh) * | 2013-09-09 | 2014-03-26 | 中华人民共和国四川出入境检验检疫局 | 基于Gabor滤波与C-V模型分割的LCD Mura缺陷机器视觉检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI375198B (en) * | 2007-05-17 | 2012-10-21 | Tpo Displays Corp | A system for displaying images |
-
2016
- 2016-05-11 CN CN201610310495.7A patent/CN105976382B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103674975A (zh) * | 2013-09-09 | 2014-03-26 | 中华人民共和国四川出入境检验检疫局 | 基于Gabor滤波与C-V模型分割的LCD Mura缺陷机器视觉检测方法 |
Non-Patent Citations (2)
Title |
---|
Pyrochemical treatment of spent nitride fuels for MA transmutation;HAYASHI Hirokazu;《Science China(Chemistry)》;20141101;第57卷(第11期);第1427-1431页 * |
基于Chan-Vese模型的TFT-LCD Mura缺陷快速分割算法;卢小鹏等;《液晶与显示》;20140215;第29卷(第1期);第146-151页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105976382A (zh) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105976382B (zh) | 一种基于缺陷区预判和水平集的TFT-LCD Mura缺陷检测方法 | |
CN106200047B (zh) | 一种基于GPU的TFT-LCD Mura缺陷检测的方法 | |
CN107657606B (zh) | 一种显示装置的亮度缺陷检测方法与装置 | |
CN108648168A (zh) | Ic晶圆表面缺陷检测方法 | |
JP3431075B2 (ja) | 液晶ディスプレイパネルムラの分類処理方法、その装置及びプログラム | |
CN103792699A (zh) | 基于B样条曲面拟合的TFT-LCD Mura缺陷机器视觉检测方法 | |
CN105069790A (zh) | 一种齿轮外观缺陷快速影像检测方法 | |
TW200844429A (en) | An automatic optical inspection approach for detecting and classifying the surface defects on coating brightness enhancement film | |
CN110378313A (zh) | 细胞团识别方法、装置及电子设备 | |
CN104978748A (zh) | 一种基于局部像素值的液晶屏缺陷检测方法 | |
Ye et al. | High-resolution optical inspection system for fast detection and classification of surface defects | |
CN102306378B (zh) | 一种图像增强方法 | |
CN106846271A (zh) | 一种去除身份证照片中网纹的方法 | |
CN109886936B (zh) | 一种低对比度缺陷检测方法和装置 | |
CN107862689A (zh) | 皮革表面明显破损自动识别方法及计算机可读存储介质 | |
CN108280825A (zh) | 一种液晶显示屏边缘缺口检测方法 | |
CN116485787B (zh) | 一种数据线成型外模外观缺陷检测方法 | |
CN110530894B (zh) | 导光板亮点缺陷检测方法 | |
Talavera-Martínez et al. | An encoder-decoder CNN for hair removal in dermoscopic images | |
CN116542963A (zh) | 一种基于机器学习的浮法玻璃缺陷检测***及检测方法 | |
CN113470015B (zh) | 一种基于图像处理的水体晃动检测分析方法及*** | |
CN112116579B (zh) | 一种透明药瓶的缺陷检测方法和装置 | |
Li et al. | MCFDFusion: Multi-focus image fusion based on multi-scale cross-difference and focus detection | |
Shuai et al. | Defect detection of four-color display screen based on color equalization and local dynamic threshold segmentation | |
CN110473190A (zh) | 一种基于尺度的自适应织物疵点检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |