CN116721082B - Display color mura defect detection method based on channel separation and filtering - Google Patents
Display color mura defect detection method based on channel separation and filtering Download PDFInfo
- Publication number
- CN116721082B CN116721082B CN202310695061.3A CN202310695061A CN116721082B CN 116721082 B CN116721082 B CN 116721082B CN 202310695061 A CN202310695061 A CN 202310695061A CN 116721082 B CN116721082 B CN 116721082B
- Authority
- CN
- China
- Prior art keywords
- channel
- feature map
- color
- contrast
- sensitive
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Image Analysis (AREA)
Abstract
本发明公开一种基于通道分离和滤波的显示屏彩色Mura缺陷检测方法,应用于图像处理技术领域,针对现有技术检测彩色Mura准确率低的问题;本发明首先,对采集的LCD显示屏RGB图像转换到Lab颜色空间;其次,对L、a和b通道进行均值模板滤波获得颜色和亮度图像,再计算Lab三通道图像的Weber对比度特征图;然后,采用频域对比度敏感函数滤波模板对Lab三通道对比度特征图进行频域滤波,获得人眼敏感的Lab三通道特征图,并融合ab通道特征图;最后,计算L通道、ab颜色融合的特征图统计量,并阈值分割获得彩色Mura缺陷检测结果。本发明方法实现了准确、高效的显示屏画面彩色Mura缺陷检测。
The present invention discloses a display screen color Mura defect detection method based on channel separation and filtering, which is applied to the field of image processing technology, and aims at the problem that the existing technology has low accuracy in detecting color Mura. The present invention first converts the collected RGB image of the LCD display screen into the Lab color space; secondly, the L, a and b channels are subjected to mean template filtering to obtain color and brightness images, and then the Weber contrast feature map of the Lab three-channel image is calculated; then, the frequency domain contrast sensitive function filtering template is used to perform frequency domain filtering on the Lab three-channel contrast feature map to obtain the Lab three-channel feature map sensitive to the human eye, and the ab channel feature map is fused; finally, the feature map statistics of the L channel and ab color fusion are calculated, and the color Mura defect detection result is obtained by threshold segmentation. The method of the present invention realizes accurate and efficient color Mura defect detection of the display screen.
Description
技术领域Technical Field
本发明属于图像处理技术领域,特别涉及一种显示屏彩色Mura缺陷检测技术。The present invention belongs to the field of image processing technology, and in particular relates to a display screen color Mura defect detection technology.
背景技术Background Art
机器视觉系统,凭借其稳定性和可靠性,被广泛应用于生产和生活中。利用机器视觉系统对显示屏质量进行检测,为工业生产带来便利,提升了生产效率,节约了成本。Machine vision systems are widely used in production and life due to their stability and reliability. Using machine vision systems to detect the quality of display screens brings convenience to industrial production, improves production efficiency and saves costs.
Mura缺陷表现为低对比度和没有固定的形状。过去的十几年中显示屏技术的快速发展,目前显示屏Mura缺陷检测方法主要为背景重构的缺陷检测法:一是采用基于频域的重构,如Fan提出Automatic detection of Mura defect in TFT-LCD based onregression diagnostics[J];二是基于拟合的重构方法,如严成宸提出结合加权模板差图与双边滤波的TFT-LCD检测算法[J],电子测量与仪器学报。然而,基于频域的重构方法不能在频域中有效地获得色彩Mura的无缺陷背景信息,这导致该方法无法从背景中分离出缺陷;基于拟合的重构方法没有考虑到色彩特征,因此直接对灰度图像拟合,无法获得缺陷区域,因此该类方法检测彩色Mura准确率低。人眼对色彩信息的感知大大超过对亮度信息的感知水平,因此采用彩色滤波分析是显示屏彩色Mura缺陷检测的有效技术途径。Mura defects are characterized by low contrast and no fixed shape. With the rapid development of display technology in the past decade, the current display Mura defect detection method is mainly a defect detection method based on background reconstruction: one is to use frequency domain-based reconstruction, such as Fan proposed Automatic detection of Mura defect in TFT-LCD based on regression diagnostics [J]; the other is a reconstruction method based on fitting, such as Yan Chengchen proposed a TFT-LCD detection algorithm combining weighted template difference map and bilateral filtering [J], Journal of Electronic Measurement and Instrumentation. However, the reconstruction method based on frequency domain cannot effectively obtain the defect-free background information of color Mura in the frequency domain, which makes this method unable to separate defects from the background; the reconstruction method based on fitting does not take color features into account, so it directly fits the grayscale image and cannot obtain the defect area. Therefore, this type of method has low accuracy in detecting color Mura. The human eye's perception of color information greatly exceeds the perception level of brightness information, so the use of color filtering analysis is an effective technical approach for display color Mura defect detection.
发明内容Summary of the invention
为解决上述技术问题,本发明提出一种基于通道分离和滤波的显示屏彩色Mura缺陷检测方法。In order to solve the above technical problems, the present invention proposes a display screen color Mura defect detection method based on channel separation and filtering.
本发明采用的技术方案为:基于通道分离和滤波的显示屏彩色Mura缺陷检测方法,包括:The technical solution adopted by the present invention is: a display color Mura defect detection method based on channel separation and filtering, comprising:
步骤S1,对待处理图像进行Lab空间通道分离;具体为:对于待处理图像,将RGB图像颜色空间转换到Lab颜色空间,获得分离的亮度和颜色通道图像;Step S1, performing Lab space channel separation on the image to be processed; specifically: for the image to be processed, converting the RGB image color space into the Lab color space to obtain separated brightness and color channel images;
步骤S2,Weber对比度特征图计算;具体为:对分离的亮度和颜色通道图像分别进行均值滤波,然后根据Weber对比度计算公式,获得亮度和颜色通道的对比度特征图;Step S2, Weber contrast feature map calculation; specifically: mean filtering the separated brightness and color channel images respectively, and then obtaining the contrast feature maps of the brightness and color channels according to the Weber contrast calculation formula;
步骤S3,对比度敏感函数频域滤波;具体为:分别构建亮度和颜色通道的频域对比度敏感函数滤波模板,使用快速傅里叶变换将亮度和颜色通道的对比度特征图转换到频域空间,并分别与对比度敏感函数进行相乘,进行逆变换后获得人眼敏感的亮度和颜色通道特征图;Step S3, contrast sensitivity function frequency domain filtering; specifically: construct frequency domain contrast sensitivity function filtering templates for brightness and color channels respectively, use fast Fourier transform to convert the contrast feature maps of brightness and color channels into frequency domain space, and multiply them with contrast sensitivity functions respectively, and obtain brightness and color channel feature maps that are sensitive to the human eye after inverse transformation;
步骤S4,颜色特征图融合及自适应阈值分割;具体为:对人眼敏感的颜色通道特征图进行融合,获得人眼敏感的颜色融合特征图,采用人眼敏感的两都通道特征图和人眼敏感的颜色特征融合图的均值和标准差进行阈值分割,将两个特征图的分割结果进行或运算,获得Mura缺陷检测结果。Step S4, color feature map fusion and adaptive threshold segmentation; specifically: fuse the color channel feature maps that are sensitive to the human eye to obtain a color fusion feature map that is sensitive to the human eye, use the mean and standard deviation of the two-channel feature maps that are sensitive to the human eye and the color feature fusion map that is sensitive to the human eye to perform threshold segmentation, perform an OR operation on the segmentation results of the two feature maps, and obtain a Mura defect detection result.
本发明的有益效果:本发明将输入图像进行亮度和颜色分离,然后分别对颜色和亮度对比度图使用对比度敏感函数滤波,采用Lab颜色空间、人眼对比度敏感函数与特征分解的综合方法进行彩色Mura缺陷检测,实现显示屏亮度Mura和颜色Mura缺陷准确检测。Beneficial effects of the present invention: The present invention separates the brightness and color of the input image, and then uses contrast sensitive functions to filter the color and brightness contrast maps respectively, and uses a comprehensive method of Lab color space, human eye contrast sensitivity function and feature decomposition to detect color Mura defects, thereby achieving accurate detection of brightness Mura and color Mura defects on the display screen.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1是本发明基于通道分离和滤波的显示屏彩色Mura缺陷检测方法流程图。FIG. 1 is a flow chart of a method for detecting color Mura defects in a display screen based on channel separation and filtering according to the present invention.
图2是原始图像转Lab空间均值滤波结果样例展示。Figure 2 is an example of the result of converting the original image to Lab space mean filtering.
图3是对比度敏感函数滤波结果样例展示。Figure 3 shows an example of the contrast sensitivity function filtering results.
图4是Mura缺陷检测结果样例展示。Figure 4 shows a sample of Mura defect detection results.
具体实施方式DETAILED DESCRIPTION
参阅图1在以下描述的实施例中,基于通道分离和滤波的显示屏彩色Mura缺陷检测方法流程按以下步骤实施:Referring to FIG. 1 , in the embodiment described below, the display color Mura defect detection method based on channel separation and filtering is implemented in the following steps:
步骤S1:输入图像Lab空间通道分离,具体步骤包括以下步骤:Step S1: Input image Lab space channel separation, specifically including the following steps:
S1-1、对采集的LCD显示屏RGB图像转换到Lab颜色空间;S1-1, convert the collected RGB image of the LCD display screen into the Lab color space;
对输入图像进行Lab颜色空间计算,计算公式如下:The Lab color space calculation is performed on the input image. The calculation formula is as follows:
其中,L为亮度通道,a为红绿颜色通道,b为蓝黄颜色通道,X、Y、Z表示XYZ颜色空间的颜色坐标,f(X)、f(Y1)、f(Z)表示XYZ颜色空间的颜色分布函数。获得分离的亮度L和颜色a、b通道图像;Where L is the brightness channel, a is the red and green color channel, b is the blue and yellow color channel, X, Y, Z represent the color coordinates in the XYZ color space, and f(X), f(Y1), and f(Z) represent the color distribution function in the XYZ color space. Get the separated brightness L and color a and b channel images;
S2:Weber对比度特征图计算,具体步骤包括以下步骤:S2: Weber contrast feature map calculation, the specific steps include the following steps:
S2-1、均值滤波处理,S2-1, mean filtering processing,
对L、a和b通道图像分别采用均值滤波模板进行滤波处理,抑制显示屏背景纹理的干扰,图像处理结果如图2所示。The L, a and b channel images are filtered using a mean filter template to suppress the interference of the background texture of the display screen. The image processing results are shown in Figure 2.
S2-2、全局均值计算,S2-2, global mean calculation,
分别计算经均值滤波处理后的L、a和b通道图像的全局均值,计算公式如下:Calculate the global mean of the L, a and b channel images after mean filtering respectively. The calculation formula is as follows:
其中μweber表示全局均值,M为输入图像宽度,N为输入图像的长度,pi,j表示图像中坐标为(i,j)的灰度值。分别计算出L通道,a和b通道的均值。Where μ weber represents the global mean, M is the width of the input image, N is the length of the input image, and pi ,j represents the grayscale value at coordinate (i,j) in the image. Calculate the mean of the L channel, a channel, and b channel respectively.
S2-3、Weber对比度计算,S2-3, Weber contrast calculation,
根据Weber对比度计算公式,获得L、a和b通道的对比度特征图。According to the Weber contrast calculation formula, the contrast feature maps of L, a and b channels are obtained.
Weber对比度计算公式如下:The Weber contrast calculation formula is as follows:
其中I(x,y)为具体通道的图像的灰度值,Ib1(x,y)为具体通道图像的背景亮度,Cw(x,y)为具体通道图像的对比度。Wherein I(x, y) is the gray value of the image of a specific channel, I b1 (x, y) is the background brightness of the image of the specific channel, and C w (x, y) is the contrast of the image of the specific channel.
将Ib1(x,y)分别替换为L、a和b的全局均值,分别计算得到CL、Ca和Cb。其中CL表示L通道对比度特征图,Ca为a通道的对比度特征图,Cb为b通道的对比度特征图。Replace I b1 (x, y) with the global mean of L, a and b, and calculate CL , Ca and Cb respectively. CL represents the contrast feature map of L channel, Ca is the contrast feature map of a channel, and Cb is the contrast feature map of b channel.
S3:对比度敏感函数频域滤波,结果图如图3所示,具体步骤包括以下步骤。S3: Contrast sensitivity function frequency domain filtering. The result diagram is shown in FIG3 . The specific steps include the following steps.
S3-1、频域对比度敏感函数构建,S3-1, frequency domain contrast sensitivity function construction,
分别构建L、a和b通道的频域对比度敏感函数。对比度敏感函数如下:Construct the frequency domain contrast sensitivity functions of L, a and b channels respectively. The contrast sensitivity functions are as follows:
其中,f表示频域的频率变量,csfL(f)表示频域L通道的对比度敏感函数,csfa(f)表示频域a通道的对比度敏感函数,csfb(f)表示频域b通道的对比度敏感函数。Among them, f represents the frequency variable in the frequency domain, csf L (f) represents the contrast sensitivity function of the L channel in the frequency domain, csf a (f) represents the contrast sensitivity function of the a channel in the frequency domain, and csf b (f) represents the contrast sensitivity function of the b channel in the frequency domain.
S3-2、通道特征图的傅里叶变换,S3-2, Fourier transform of channel feature map,
对人眼敏感的L、a和b通道特征图,进行傅里叶变换,并将低频信息移动到图像中心,得到频域的L、a和b通道特征图,分别记为CL_w、Ca_w和Cb_w。The L, a and b channel feature maps that are sensitive to the human eye are Fourier transformed, and the low-frequency information is moved to the center of the image to obtain the L, a and b channel feature maps in the frequency domain, which are recorded as C L_w , Ca_w and C b_w respectively.
S3-3、特征图频域滤波,S3-3, feature map frequency domain filtering,
将得到的频域的L、a和b通道特征图分别与对比度敏感函数相乘,获得通道滤波特征结果,计算公式如下所示。The obtained L, a and b channel feature maps in the frequency domain are multiplied by the contrast sensitivity function respectively to obtain the channel filtering feature results. The calculation formula is as follows.
其中Ccsf_L为人眼敏感的L通道特征图,Ccsf_a为人眼敏感的a通道特征图,Ccsf_b为人眼敏感的b通道特征图,idft表示逆傅里叶变换,结果如图3所示。Where C csf_L is the L channel feature map sensitive to the human eye, C csf_a is the a channel feature map sensitive to the human eye, C csf_b is the b channel feature map sensitive to the human eye, idft means inverse Fourier transform, and the result is shown in Figure 3.
步骤S4颜色特征图融合及自适应阈值分割,具体步骤包括以下步骤。Step S4: color feature map fusion and adaptive threshold segmentation, specifically including the following steps.
S4-1、颜色特征图融合,S4-1, color feature map fusion,
将人眼敏感的a和b通道特征图进行融合,获得人眼敏感的颜色融合特征图,计算公式如下:The a and b channel feature maps that are sensitive to the human eye are fused to obtain a color fusion feature map that is sensitive to the human eye. The calculation formula is as follows:
其中Ccsf_col表示人眼敏感的颜色融合特征图。Where C csf_col represents the color fusion feature map that the human eye is sensitive to.
S4-2、亮度和颜色的自适应阈值分割,S4-2, adaptive threshold segmentation of brightness and color,
采用人眼敏感的L通道特征图和人眼敏感的颜色特征融合图的均值和标准差进行阈值分割,计算公式如下:The mean and standard deviation of the L channel feature map that the human eye is sensitive to and the color feature fusion map that the human eye is sensitive to are used for threshold segmentation. The calculation formula is as follows:
thcol=μccol+K1δccol th col =μc col +K1δc col
thL=μcL+K2δcL th L = μc L + K2δc L
其中,thcol为颜色融合特征图的分割阈值,K1为参数,μccol为颜色融合特征图的均值,μccol为颜色融合特征图的标准差,thL为人眼敏感的L通道特征图分割阈值,K2为参数,μcL为人眼敏感的L通道特征图均值,δcL为人眼敏感的L通道特征图标准差。本实施例中K1取值为1,K2取值为3。Among them, th col is the segmentation threshold of the color fusion feature map, K1 is a parameter, μc col is the mean of the color fusion feature map, μc col is the standard deviation of the color fusion feature map, th L is the segmentation threshold of the L channel feature map sensitive to the human eye, K2 is a parameter, μc L is the mean of the L channel feature map sensitive to the human eye, and δc L is the standard deviation of the L channel feature map sensitive to the human eye. In this embodiment, K1 takes a value of 1 and K2 takes a value of 3.
将两个特征图的分割结果进行或运算,获得Mura缺陷检测结果,如图4所示。The segmentation results of the two feature maps are ORed to obtain the Mura defect detection result, as shown in FIG4 .
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。Those skilled in the art will appreciate that the embodiments described herein are intended to help readers understand the principles of the present invention, and should be understood that the scope of protection of the present invention is not limited to such specific statements and embodiments. For those skilled in the art, the present invention may have various changes and variations. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of the claims of the present invention.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310695061.3A CN116721082B (en) | 2023-06-13 | 2023-06-13 | Display color mura defect detection method based on channel separation and filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310695061.3A CN116721082B (en) | 2023-06-13 | 2023-06-13 | Display color mura defect detection method based on channel separation and filtering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116721082A CN116721082A (en) | 2023-09-08 |
CN116721082B true CN116721082B (en) | 2024-09-06 |
Family
ID=87869345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310695061.3A Active CN116721082B (en) | 2023-06-13 | 2023-06-13 | Display color mura defect detection method based on channel separation and filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116721082B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819850A (en) * | 2012-08-16 | 2012-12-12 | 李勃 | Method for detecting edge of color image on basis of local self-adaption color difference threshold |
CN109345528A (en) * | 2018-09-28 | 2019-02-15 | 凌云光技术集团有限责任公司 | A kind of display screen defect inspection method and device based on human-eye visual characteristic |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11676265B2 (en) * | 2020-07-31 | 2023-06-13 | Novatek Microelectronics Corp. | Method and image processing device for mura detection on display |
CN112954304B (en) * | 2021-01-18 | 2022-09-16 | 湖北经济学院 | Mura defect assessment method for display panel |
CN113724216B (en) * | 2021-08-24 | 2023-03-21 | 广东工业大学 | Method and system for detecting wave crest welding spot defects |
CN114299070B (en) * | 2022-03-09 | 2022-05-31 | 深圳精智达技术股份有限公司 | Method and related device for detecting mura defects of display screen |
CN114998189B (en) * | 2022-04-15 | 2024-04-16 | 电子科技大学 | Color display point defect detection method |
CN114820521B (en) * | 2022-04-27 | 2023-04-07 | 电子科技大学 | Defect detection method for complex picture of display screen |
CN115187551B (en) * | 2022-07-14 | 2024-12-27 | 电子科技大学 | A display screen weak cross defect detection method based on frequency domain filtering of significant color channels |
CN116129099A (en) * | 2022-12-30 | 2023-05-16 | 浙江中烟工业有限责任公司 | A machine vision-based intelligent recognition method for shredded foreign objects |
-
2023
- 2023-06-13 CN CN202310695061.3A patent/CN116721082B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102819850A (en) * | 2012-08-16 | 2012-12-12 | 李勃 | Method for detecting edge of color image on basis of local self-adaption color difference threshold |
CN109345528A (en) * | 2018-09-28 | 2019-02-15 | 凌云光技术集团有限责任公司 | A kind of display screen defect inspection method and device based on human-eye visual characteristic |
Also Published As
Publication number | Publication date |
---|---|
CN116721082A (en) | 2023-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20240005757A1 (en) | Method and device for smoke or fire recognition, computer apparatus, and storage medium | |
CN106650770B (en) | A mura defect detection method based on sample learning and human visual characteristics | |
CN109507192B (en) | A method for detecting surface defects of magnetic cores based on machine vision | |
WO2021143343A1 (en) | Method and device for testing product quality | |
CN107505546A (en) | A kind of method that corona discharge is monitored using ultraviolet imager | |
CN111325721A (en) | Gas leakage detection method and system based on infrared thermal imaging | |
CN111402247A (en) | A machine vision-based detection method for hanging clip defects on transmission lines | |
CN109544464A (en) | A kind of fire video image analysis method based on contours extract | |
CN105069778B (en) | Based on the industrial products detection method of surface flaw that target signature notable figure builds | |
CN107644538B (en) | Traffic signal lamp identification method and device | |
CN113505865A (en) | Sheet surface defect image recognition processing method based on convolutional neural network | |
CN111724375B (en) | A screen detection method and system | |
CN102074031A (en) | A method for establishing a standard for a printed circuit board appearance inspection machine | |
CN114486916A (en) | Defect detection method of mobile phone glass cover based on machine vision | |
CN113205136A (en) | Real-time high-precision detection method for appearance defects of power adapter | |
CN115078365A (en) | Soft package printing quality defect detection method | |
CN105844281A (en) | Automatic form set parameter acquiring system and method | |
CN114998189B (en) | Color display point defect detection method | |
CN110186929A (en) | A kind of real-time product defect localization method | |
CN104809735B (en) | The system and method for image haze evaluation is realized based on Fourier transformation | |
CN116721082B (en) | Display color mura defect detection method based on channel separation and filtering | |
CN102118547A (en) | Image weighted filtering method | |
CN107330441B (en) | Flame Image Foreground Extraction Algorithm | |
CN115359031A (en) | Digital pathological image slice quality evaluation method | |
CN118688207A (en) | Intelligent multi-sorting control system and method for polarizer surface defects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |