WO2021000524A1 - Hole protection cap detection method and apparatus, computer device and storage medium - Google Patents

Hole protection cap detection method and apparatus, computer device and storage medium Download PDF

Info

Publication number
WO2021000524A1
WO2021000524A1 PCT/CN2019/125058 CN2019125058W WO2021000524A1 WO 2021000524 A1 WO2021000524 A1 WO 2021000524A1 CN 2019125058 W CN2019125058 W CN 2019125058W WO 2021000524 A1 WO2021000524 A1 WO 2021000524A1
Authority
WO
WIPO (PCT)
Prior art keywords
gray
preset
image
pixel
hole
Prior art date
Application number
PCT/CN2019/125058
Other languages
French (fr)
Chinese (zh)
Inventor
戴志威
邓远志
陈润康
林淼
刘志永
陈志列
Original Assignee
研祥智能科技股份有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 研祥智能科技股份有限公司 filed Critical 研祥智能科技股份有限公司
Publication of WO2021000524A1 publication Critical patent/WO2021000524A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/0008Industrial image inspection checking presence/absence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration by the use of histogram techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Image Analysis (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Image Processing (AREA)

Abstract

A hole protection cap detection method, comprising: acquiring an image of a device to be detected, said device being provided with a hole having a hole protection cap; performing self-adaptive dynamic threshold segmentation on the image of said device, to determine a hole region; then extracting a greyscale feature subset of the hole region; projecting the greyscale feature subset onto a preset discriminant function line in a two-dimensional coordinate system, to acquire a corresponding projected longitudinal coordinate value; obtaining, according to a magnitude relationship between the longitudinal coordinate value and the threshold of the preset discriminant function line, a defect state of the device corresponding to said image. The threshold of the preset discriminant function line is obtained by the following method: on the basis of extracted greyscale feature subsets of a plurality of hole regions of a plurality of historical images of said device, constructing a training sample set containing at least two types of sample subsets, projecting various types of sample subsets into an optimal discrimination vector space, and calculating a corresponding projected horizontal coordinate mean value, to obtain the threshold of the preset discriminant function line.

Description

孔位保护门检测方法、装置、计算机设备和存储介质Method and device for detecting hole position protection door, computer equipment and storage medium
相关申请的交叉引用Cross references to related applications
本申请要求于2019年7月3日提交中国专利局,申请号为2019105938923,申请名称为“孔位保护门检测方法、装置、计算机设备和存储介质”的中国专利申请的优先权,其全部内容通过引用结合在本申请中。This application claims the priority of the Chinese patent application filed with the Chinese Patent Office on July 3, 2019. The application number is 2019105938923, and the application title is "hole position protection door detection method, device, computer equipment and storage medium." Incorporated in this application by reference.
技术领域Technical field
本申请涉及设备检测技术领域,特别是涉及一种孔位保护门检测方法、装置、计算机设备和存储介质。This application relates to the technical field of equipment detection, and in particular to a detection method, device, computer equipment and storage medium for a hole position protection door.
背景技术Background technique
保护门又称安全门,是一种安全保护盖,其一般设置于孔位上,平时孔位未被使用时,是挡住孔位,防止异物进入孔位存在安全隐患。The protective door, also known as the safety door, is a kind of safety protective cover, which is generally set on the hole position. When the hole position is not used, it blocks the hole position and prevents foreign objects from entering the hole position.
在设置有保护门的产品生产制造过程中,需要检测产品的保护门是否能正常复位,而保护门是否能正常复位的检测通常是通过人工目测来把控产品的生产质量。由于人工检测因每个人的评判标准不同,且人的感官判断易受个人状态、情绪等主观因素影响,使得产品的检测的效率低且误检率高。In the manufacturing process of products with protective doors, it is necessary to detect whether the protective door of the product can be reset normally, and the detection of whether the protective door can be reset normally is usually through manual visual inspection to control the production quality of the product. Due to the different evaluation criteria of each person for manual detection, and the human sensory judgment is easily affected by subjective factors such as personal state and emotion, the efficiency of product detection is low and the false detection rate is high.
为提高自动化水平,市场上出现了一些自动化检测设置有孔位保护门产品的方法,但这些孔位保护门状态检测方法多是通过排布好的流水线式的测试仪进行检测或借助二插或三插保护门插拔力检测机构等进行检测,检测仪器在运转过程存在许多干扰因素,会使得目前的孔位保护门状态检测方法存在检测结果不稳定且检测效率不高的问题。In order to improve the level of automation, some automatic detection methods for products equipped with hole protection doors have appeared on the market, but most of these hole protection door status detection methods are detected by a well-arranged assembly line tester or with the help of two plugs or The three-insertion protection door plug-in force detection mechanism, etc. perform detection. There are many interference factors in the operation of the detection instrument, which will cause the current hole position protection door state detection method to have the problem of unstable detection results and low detection efficiency.
发明内容Summary of the invention
根据本申请公开的各种实施例,提供一种孔位保护门检测方法、装置、计算机设备和存储介质。According to various embodiments disclosed in this application, a method, device, computer equipment and storage medium for detecting a hole position protection door are provided.
一种孔位保护门检测方法包括:A method for detecting hole position protection door includes:
获取待检测设备图像,待检测设备上设置有孔位,孔位上设置有孔位保护门;Obtain the image of the device to be tested, the device to be tested is provided with a hole position, and the hole position is provided with a hole position protection door;
对待检测设备图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Perform adaptive dynamic threshold segmentation based on local statistics on the image of the equipment to be tested to determine the hole location area;
基于灰度直方图的统计算法,提取出孔位区域的灰度特征子集;Based on the statistical algorithm of the gray histogram, extract the gray feature subset of the hole location area;
将灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取灰度特征子集的投影纵坐标值,预设的判别函数直线为使得灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;及Project the gray feature subset onto the preset discriminant function line in the two-dimensional coordinate system to obtain the projected ordinate value of the gray feature subset. The preset discriminant function line is such that the gray feature subset is A straight line with the largest inter-class distance and the smallest intra-class distance in the subspace; and
根据灰度特征子集投影后的纵坐标值与预设的判别函数直线的阈值的大小关系,得到 待检测设备的图像对应的设备的缺陷状态;Obtain the defect state of the device corresponding to the image of the device to be tested according to the magnitude relationship between the ordinate value after the projection of the gray feature subset and the preset threshold of the discriminant function line;
预设的判别函数直线的阈值由以下方式得到:The preset threshold of the straight line of the discriminant function is obtained in the following way:
获取多张待检测设备的历史图像;Obtain multiple historical images of the equipment to be tested;
对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Perform adaptive dynamic threshold segmentation based on local statistics on the historical image of the equipment to be tested to determine the hole location area;
基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,训练样本集至少包括两类样本子集;及Based on the statistical algorithm of the gray histogram, extract the gray feature subsets of multiple hole locations, and construct a training sample set, which includes at least two types of sample subsets; and
将各类样本子集投影到最佳鉴别矢量空间,分别计算各类样本子集投影后的横坐标均值,得到预设的判别函数直线的阈值,最佳鉴别矢量空间为基于Fisher线性判别分析算法和训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。Project various sample subsets to the best discriminant vector space, calculate the mean value of the abscissa after the projection of the various sample subsets, and get the preset threshold of the discriminant function straight line. The best discriminant vector space is based on the Fisher linear discriminant analysis algorithm And the training sample set, introducing the best projection direction solved by Lagrangian multiplier method.
一种孔位保护门状态检测装置包括:A device for detecting the status of a hole position protection door includes:
图像获取模块,用于读取待检测设备的图像,获取待检测设备的图像,待检测设备上设置有孔位,孔位上设置有孔位保护门;The image acquisition module is used to read the image of the equipment to be tested and acquire the image of the equipment to be tested. The equipment to be tested is provided with a hole position and a hole position protection door is arranged on the hole position;
孔位区域获取模块,用于对待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;The hole location area acquisition module is used to perform adaptive dynamic threshold segmentation based on local statistics on the image of the equipment to be detected to determine the hole location area;
灰度特征提取模块,用于基于灰度直方图的统计算法,提取出孔位区域的灰度特征子集;The gray-level feature extraction module is used to extract the gray-level feature subset of the hole location area based on the statistical algorithm of the gray-level histogram;
投影模块,用于将灰度特征子集投影到二维坐标系中预设的判别函数直线上,计算灰度特征子集的投影纵坐标值,预设的判别函数直线为使得灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;The projection module is used to project the gray feature subset onto the preset discriminant function line in the two-dimensional coordinate system, and calculate the projected ordinate value of the gray feature subset. The preset discriminant function line is such that the gray feature Set the straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace;
缺陷检测模块,用于根据灰度特征子集投影后的纵坐标值与预设判别函数直线的阈值的大小关系,得到待检测设备的图像对应的设备的缺陷状态;及The defect detection module is used to obtain the defect state of the device corresponding to the image of the device to be detected according to the magnitude relationship between the ordinate value after the projection of the gray feature subset and the threshold value of the preset discriminant function line; and
预设阈值计算模块,用于获取多张待检测设备的历史图像,对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域,基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,训练样本集至少包括两类样本子集,将各类样本子集投影到最佳鉴别矢量空间,分别计算各类样本子集投影后的横坐标均值,得到预设的判别函数直线的阈值,最佳鉴别矢量空间为基于Fisher线性判别分析算法和训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。The preset threshold calculation module is used to obtain multiple historical images of the equipment to be inspected, perform adaptive dynamic threshold segmentation based on local statistics on the historical images of the equipment to be inspected, determine the hole location area, and extract the historical image based on the gray histogram Construct a training sample set of gray feature subsets of multiple hole locations. The training sample set includes at least two types of sample subsets. Project various sample subsets to the best discriminant vector space, and calculate the projections of various sample subsets. After the mean value of the abscissa, the preset threshold of the discriminant function line is obtained. The best discriminant vector space is the best projection direction based on the Fisher linear discriminant analysis algorithm and training sample set, and the Lagrange multiplier method is introduced.
一种计算机设备,包括存储器和一个或多个处理器,存储器中存储有计算机可读指令,计算机可读指令被处理器执行时实现本申请任意一个实施例中提供的孔位保护门状态检测方法的步骤。A computer device, including a memory and one or more processors. The memory stores computer-readable instructions. When the computer-readable instructions are executed by the processor, the method for detecting the status of the hole position protection door provided in any embodiment of the present application is implemented A step of.
一个或多个存储有计算机可读指令的非易失性存储介质,计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器实现本申请任意一个实施例中提供的孔位保护门状态检测方法的步骤。One or more non-volatile storage media storing computer-readable instructions. When the computer-readable instructions are executed by one or more processors, the one or more processors implement the holes provided in any of the embodiments of the present application. The steps of the method for detecting the status of the bit protection gate.
本申请的一个或多个实施例的细节在下面的附图和描述中提出。本申请的其它特征和优点将从说明书、附图以及权利要求书变得明显。The details of one or more embodiments of the application are set forth in the following drawings and description. Other features and advantages of this application will become apparent from the description, drawings and claims.
附图说明Description of the drawings
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。In order to more clearly describe the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings needed in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative work.
图1为根据一个或多个实施例中孔位保护门状态检测方法的应用场景图。Fig. 1 is an application scene diagram of the method for detecting the state of the hole position protection door according to one or more embodiments.
图2为根据一个或多个实施例中孔位保护门状态检测方法的流程示意图。Fig. 2 is a schematic flowchart of a method for detecting the status of a hole position protection door according to one or more embodiments.
图3为根据一个或多个实施例中采集待检测设备的图像的应用环境图。Fig. 3 is a diagram of an application environment for collecting images of a device to be tested according to one or more embodiments.
图4为另一个实施例中预设的判别函数直线的阈值获取步骤的流程示意图。FIG. 4 is a schematic flow chart of the step of obtaining the threshold value of the preset discriminant function straight line in another embodiment.
图5为另一个实施例中孔位区域获取子步骤的流程示意图。Fig. 5 is a schematic flowchart of a sub-step of hole location area acquisition in another embodiment.
图6为根据一个或多个实施例中孔位保护门状态检测装置的框图。Fig. 6 is a block diagram of a device for detecting the status of a hole position protection door according to one or more embodiments.
图7为另一个实施例中孔位保护门状态检测装置的结构图。Fig. 7 is a structural diagram of a hole position protection door state detection device in another embodiment.
图8为根据一个或多个实施例中计算机设备的框图。Figure 8 is a block diagram of a computer device according to one or more embodiments.
具体实施方式Detailed ways
为了使本申请的技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。In order to make the technical solutions and advantages of the present application clearer, the following further describes the present application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the application, and not used to limit the application.
本申请提供的孔位保护门状态检测方法,可应用于如图1中所示的应用环境图中。服务器120与摄像装置103通过网络进行通信。服务器120获取摄像装置103上传的待检测设备的图像,对待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域(孔位区域包含可疑缺陷区域集合),基于灰度直方图的统计算法,提取出孔位区域的灰度特征子集,将灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取灰度特征子集的投影纵坐标值,预设的判别函数直线为使得灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线,根据灰度特征子集投影后的纵坐标值与预设判别函数直线的阈值的大小关系,得到待检测设备的图像对应的设备的缺陷状态。其中,摄像装置103可以但不限于是摄像头、具有拍照功能的智能手机以及摄像设备等,服务器120可以用独立的服务器或者是多个服务器组成的服务器集群来实现。The method for detecting the status of the hole position protection door provided in this application can be applied to the application environment diagram shown in FIG. 1. The server 120 and the camera 103 communicate through a network. The server 120 obtains the image of the device to be tested uploaded by the camera 103, performs adaptive dynamic threshold segmentation based on local statistics on the image of the device to be tested, and determines the hole location area (the hole location area includes a set of suspicious defect areas), based on the gray histogram Based on the statistical algorithm, extract the gray feature subset of the hole location area, project the gray feature subset onto the preset discriminant function line in the two-dimensional coordinate system, and obtain the projected ordinate value of the gray feature subset. Suppose the discriminant function line is the line that makes the gray feature subset have the largest inter-class distance and the smallest intra-class distance in the projected subspace, according to the projected ordinate value of the gray feature subset and the preset discriminant function line The relationship between the magnitude of the threshold and the defect state of the device corresponding to the image of the device to be detected is obtained. The camera 103 can be, but is not limited to, a camera, a smartphone with a camera function, a camera device, etc. The server 120 can be implemented by an independent server or a server cluster composed of multiple servers.
在其中一个实施例中,如图2所示,提供了一种孔位保护门状态检测方法,以该方法应用于图1中的服务器为例进行说明,包括以下步骤:In one of the embodiments, as shown in FIG. 2, a method for detecting the status of a hole position protection door is provided. Taking the method applied to the server in FIG. 1 as an example, the method includes the following steps:
步骤S100,获取待检测设备的图像,待检测设备上设置有孔位,孔位上设置有孔位保护门。In step S100, an image of a device to be detected is acquired. The device to be detected is provided with a hole position and a hole position protection door is provided on the hole position.
待检测设备的图像即为摄像装置103采集到的设置有孔位及孔位保护门的待检测设备的图像,图像信息包括孔位区域图像,其中,孔位区域图像包括保护门图像。待检测设备以插座为例,待检测设备的图像即包括插孔区域图像以及插孔上的保护门图像。如图3 所示,结合图1,在实际应用中,待检测设备的图像的获取方式可以是当处于流水线上的待检测设备101(带孔位那一面向上,孔位上设置有孔位保护门)经过光电传感器105时,光电传感器105将触发信号发送给工控机104,工控机104收到信号后通过IO板卡输出控制信号,首先打开置于待检测设备101正上方的环形无影光源,然后给摄像装置103发送软触发信号,采集到待检测设备101图像,然后摄像装置103将采集的图像上传至服务器120,服务器120获取待检测设备的图像。The image of the equipment to be inspected is the image of the equipment to be inspected with the hole position and the hole position protection door collected by the camera 103. The image information includes the hole position area image, where the hole position area image includes the protection door image. Take the socket as an example of the device to be tested. The image of the device to be tested includes the image of the jack area and the image of the protective door on the jack. As shown in Figure 3, in conjunction with Figure 1, in practical applications, the image of the device to be tested can be acquired when the device to be tested 101 on the assembly line (the side with the hole position is up, and the hole position is provided with hole position protection Door) When passing the photoelectric sensor 105, the photoelectric sensor 105 sends the trigger signal to the industrial computer 104. After receiving the signal, the industrial computer 104 outputs the control signal through the IO board. First, turn on the ring-shaped shadowless light source directly above the device to be tested 101 , And then send a soft trigger signal to the camera 103 to collect the image of the device to be detected 101, and then the camera 103 uploads the collected image to the server 120, and the server 120 obtains the image of the device to be detected.
步骤S200,对待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域。Step S200: Perform adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area.
通过摄像装置采集到的待检测设备的图像可发现,图像中的孔位区域与待检测设备的正面(设置有孔位的那一面定义为正面)有明显的灰度差,对比度较为明显,因此,可首先通过阈值分割,来提取出待测正面的感兴趣区域,即正面所有待检测的孔位区域及保护门区域。由于工厂的作业环境比较复杂,且需要同时满足多种颜色待检测设备的检测,本申请采用基于局部统计的自适应动态阈值分割法来提取孔位区域,其中,孔位区域为包含可疑缺陷区域集合。局部阈值分割是一种处理亮度分布不均匀的常用方法,其在一些光照不足的阴影区域有良好的分割效果。设待检测设备的原始图像为f(x,y),对原始图像进行平滑滤波器处理得到的图像设为m(x,y),滤波器可以选择为均值滤波器、二项式滤波器及高斯滤波器等。设定偏移量为Offset,对f(x,y)与m(x,y)进行逐一的像素灰度值比较,当待分割的目标比背景亮时,满足f(x,y)≥m(x,y)+Offset的像素点位置输出为1;当待分割的目标比背景暗时,满足f(x,y)≤m(x,y)-Offset的像素点位置输出1;当待分割的目标既有可能亮于背景,也有可能暗于背景区域时,需要同时满足上述两个条件。本申请在在局部阈值分割的基础上,提出一种基于局部统计的自适应动态阈值分割方法,在窗口的平均值信息基础上,通过增加对窗口的标准方差统计信息来修正偏移量Offset。Through the image of the equipment to be tested collected by the camera device, it can be found that the hole location area in the image and the front of the equipment to be tested (the side with the hole location is defined as the front) has obvious gray difference, and the contrast is more obvious. , Can first extract the area of interest on the front side to be tested by threshold segmentation, that is, all the hole area to be detected on the front side and the protection door area. Due to the complex operating environment of the factory and the need to meet the detection of multiple color equipment to be inspected at the same time, this application uses an adaptive dynamic threshold segmentation method based on local statistics to extract the hole location area, where the hole location area contains suspicious defects. set. Local threshold segmentation is a common method to deal with uneven brightness distribution, and it has a good segmentation effect in some shadow areas with insufficient light. Suppose the original image of the device to be tested is f(x,y), and the image obtained by smoothing filter processing on the original image is set to m(x,y). The filter can be selected as mean filter, binomial filter and Gaussian filter, etc. Set the offset to Offset, compare the pixel gray values of f(x,y) and m(x,y) one by one. When the target to be segmented is brighter than the background, satisfy f(x,y)≥m The pixel position of (x,y)+Offset is output as 1; when the target to be segmented is darker than the background, the pixel position that satisfies f(x,y)≤m(x,y)-Offset is output as 1; When the segmented target may be brighter than the background or darker than the background area, the above two conditions need to be met at the same time. On the basis of local threshold segmentation, this application proposes an adaptive dynamic threshold segmentation method based on local statistics. Based on the average information of the window, the offset Offset is corrected by adding the statistical information of the standard deviation of the window.
步骤S300,基于灰度直方图的统计算法,提取出孔位区域的灰度特征子集。Step S300, based on the statistical algorithm of the gray histogram, extract the gray feature subset of the hole location area.
本实施例中,基于灰度直方图的统计算法的统计特性提取孔位区域的灰度特征。对于一副像素总数为N的二维图像,其灰度等级为L,同一灰度级的像素总数是n(i),则一阶灰度直方图可定义为:In this embodiment, the gray level feature of the hole location area is extracted based on the statistical characteristics of the gray histogram statistical algorithm. For a two-dimensional image with a total number of pixels N, its gray level is L, and the total number of pixels of the same gray level is n(i), then the first-level grayscale histogram can be defined as:
Figure PCTCN2019125058-appb-000001
Figure PCTCN2019125058-appb-000001
通过灰度直方图可以总体反映出孔位区域总体灰度水平及灰度分布概率。根据灰度直方图的统计特性,可以提取如下的灰度特征:The gray level histogram can generally reflect the overall gray level and gray distribution probability of the hole location area. According to the statistical characteristics of the gray histogram, the following gray features can be extracted:
(1)均值m,平均灰度度量,反映图像总体的灰度水平。(1) The mean value m, the average gray scale, reflects the overall gray level of the image.
Figure PCTCN2019125058-appb-000002
Figure PCTCN2019125058-appb-000002
(2)标准差σ,平均对比度度量,它反映灰度直方图的离散程度。(2) Standard deviation σ, a measure of average contrast, which reflects the degree of dispersion of the gray histogram.
Figure PCTCN2019125058-appb-000003
Figure PCTCN2019125058-appb-000003
计算孔位区域的灰度均值及标准差后,根据灰度均值和标准差构建灰度特征子集T={(m T,σ T)},将灰度特征子集T作为测试样本集。 After calculating the gray average value and standard deviation of the hole location area, construct the gray feature subset T={(m TT )} according to the gray average and standard deviation, and use the gray feature subset T as the test sample set.
在其中一个实施例中,提取出孔位区域的灰度特征子集之前,还包括:根据预设的滤波器,对孔位区域进行平滑滤波处理,模糊孔位区域中的纹理信息。In one of the embodiments, before extracting the gray-scale feature subset of the hole location area, the method further includes: performing smooth filtering processing on the hole location area according to a preset filter to blur the texture information in the hole location area.
本实施例中,采用预设的滤波器为3*3大小的均值滤波器,对孔位区域集合D进行平滑滤波处理,能够模糊图像中孔位保护门的部分纹理信息,排除纹理部分对区域灰度标准差的干扰。可以理解的是,在其他实施例中,均值滤波器还可以为其他尺寸大小,滤波器可以是二项式滤波器及高斯滤波器等。In this embodiment, the preset filter is an average filter with a size of 3*3, and the hole position area set D is smoothed and filtered, which can blur part of the texture information of the hole position protection gate in the image, and exclude the texture part of the area. Interference of gray standard deviation. It can be understood that, in other embodiments, the average filter may also have other sizes, and the filter may be a binomial filter, a Gaussian filter, or the like.
步骤S400,将灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取灰度特征子集的投影纵坐标值,预设的判别函数直线为使得灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线。Step S400: Project the gray feature subset onto a preset discriminant function line in the two-dimensional coordinate system, and obtain the projected ordinate value of the gray feature subset. The preset discriminant function line is such that the gray feature subset is in There is a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace.
在实际应用中,工作人员假设在二维坐标系空间中,存在一直线(判别函数直线)w,使得在这个空间中所有特征点投影到此直线后,类内方差最小,类间方差最大。然后,需要计算出判别函数直线的阈值,然后,将灰度特征子集投影到该判别函数直线上,获取灰度特征子集的在该直线上的投影纵坐标值,然后,可以是将缺陷状态的检测转化为一个二分类问题,基于判别函数直线的阈值与投影纵坐标值的大小关系,得到待检测设备的保护门的缺陷状态。In practical applications, the staff assumes that there is a straight line (discrimination function straight line) w in the two-dimensional coordinate system space, so that after all feature points in this space are projected onto this straight line, the intra-class variance is the smallest and the inter-class variance is the largest. Then, it is necessary to calculate the threshold value of the discriminant function line, and then project the gray feature subset onto the discriminant function line to obtain the projected ordinate value of the gray feature subset on the line. Then, it can be the defect The state detection is transformed into a two-class problem. Based on the relationship between the threshold value of the discriminant function line and the projected ordinate value, the defect state of the protection door of the equipment to be tested is obtained.
在其中一个实施例中,如图4所示,预设判别函数直线的阈值由以下方式得到:In one of the embodiments, as shown in FIG. 4, the threshold of the preset discriminant function straight line is obtained in the following manner:
步骤S402,获取多张待检测设备的历史图像;Step S402, acquiring multiple historical images of the equipment to be tested;
步骤S404,对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Step S404: Perform adaptive dynamic threshold segmentation based on local statistics on the historical image of the device to be detected to determine the hole location area;
步骤S406,基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,训练样本集至少包括两类样本子集;Step S406, based on the statistical algorithm of the gray histogram, extract the gray feature subsets of multiple hole location regions, and construct a training sample set, the training sample set includes at least two types of sample subsets;
步骤S408,将各类样本子集投影到最佳鉴别矢量空间,分别计算各类样本子集投影后的横坐标均值,得到预设的判别函数直线的阈值,最佳鉴别矢量空间为基于Fisher线性判别分析算法和训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。Step S408: Project various sample subsets into the optimal discriminative vector space, and calculate the mean values of abscissas after the projection of the various sample subsets to obtain the preset discriminant function straight line threshold. The optimal discriminant vector space is based on Fisher linear Discriminant analysis algorithm and training sample set, introducing the best projection direction solved by Lagrangian multiplier method.
具体的,可以是读取50个待检测设备的历史图像,提取出待检测设备的历史图像中的孔位区域作为训练样本集,划分为两类子集,其中包括25个合格孔位区域灰度特征子集X 0和25个不合格孔位区域灰度特征子集X 1Specifically, it can read 50 historical images of the equipment to be tested, extract the hole location area in the historical image of the equipment to be tested as a training sample set, and divide it into two types of subsets, including 25 qualified hole locations. Degree feature subset X 0 and 25 unqualified hole location regions gray feature subset X 1 :
X 0={(m 11),(m 22),(m 33),…,(m 2525)} X 0 ={(m 11 ),(m 22 ),(m 33 ),…,(m 2525 )}
X 1={(m 2626),(m 2727),(m 2828),…,(m 5050)} X 1 ={(m 2626 ),(m 2727 ),(m 2828 ),…,(m 5050 )}
计算各类样本的均值向量
Figure PCTCN2019125058-appb-000004
Calculate the mean vector of various samples
Figure PCTCN2019125058-appb-000004
计算各类样本的协方差矩阵
Figure PCTCN2019125058-appb-000005
Calculate the covariance matrix of various samples
Figure PCTCN2019125058-appb-000005
定义两类样本的中心在直线上的投影为:w Tu 0和w Tu 1,两类样本的协方差:w TC 0w和w TC 1w,为使得同类的样本的投影点尽可能接近,即需要保证w TC 0w+w TC 1w尽可能小,为使得异类样本的投影点尽可能远离,即需要保证
Figure PCTCN2019125058-appb-000006
尽可能大,于是,最大化
Define the projection of the centers of the two types of samples on a straight line as: w T u 0 and w T u 1 , the covariance of the two types of samples: w T C 0 w and w T C 1 w, which are the projection points of the samples of the same kind As close as possible, that is, it is necessary to ensure that w T C 0 w+w T C 1 w is as small as possible, in order to make the projection points of heterogeneous samples as far away as possible, that is, to ensure
Figure PCTCN2019125058-appb-000006
As large as possible, so maximize
Figure PCTCN2019125058-appb-000007
Figure PCTCN2019125058-appb-000007
计算类内散度矩阵:Calculate the intra-class divergence matrix:
S w=C 0+C 1 S w =C 0 +C 1
计算类间散度矩阵:Calculate the divergence matrix between classes:
S b=(u 0-u 1)(u 0-u 1) T S b =(u 0 -u 1 )(u 0 -u 1 ) T
此时,则可重写J为:At this time, J can be rewritten as:
Figure PCTCN2019125058-appb-000008
Figure PCTCN2019125058-appb-000008
令w TS ww=1,最大化广义瑞利商等价形式为: Let w T S w w = 1, the equivalent form of maximizing generalized Rayleigh quotient is:
Figure PCTCN2019125058-appb-000009
Figure PCTCN2019125058-appb-000009
s.t.w TS ww=1 stw T S w w = 1
运用拉格朗日乘子法,则有S bw=λS ww;S bw的方向恒为u 0-u 1,不妨令S bw=λ(u 0-u 1),于是可得
Figure PCTCN2019125058-appb-000010
计算矩阵S w的逆矩阵
Figure PCTCN2019125058-appb-000011
再计算w最大特征值对应的特征向量,该特征向量就是投影方向W,W为一维列向量,此时可计算得最佳投影线的斜率k=W(2)/W(1),其中,W(2)为一维列向量中第二个元素值,W(1)为一维列向量中第一个元素值,最终求得判别函数即为y=kx+b;本申请中,b=0;计算训练样本集中所有灰度特征子集点在直线y=kx+b上的投影坐标,根据点到直线的投影公式可得:
Using Lagrangian multiplier method, then S b w = λS w w; the direction of S b w is always u 0 -u 1 , so let S b w = λ(u 0 -u 1 ), so we can get
Figure PCTCN2019125058-appb-000010
Calculate the inverse of the matrix S w
Figure PCTCN2019125058-appb-000011
Then calculate the eigenvector corresponding to the maximum eigenvalue of w. The eigenvector is the projection direction W, and W is a one-dimensional column vector. At this time, the slope of the best projection line k=W(2)/W(1) can be calculated, where , W(2) is the second element value in the one-dimensional column vector, W(1) is the first element value in the one-dimensional column vector, and the final discriminant function is y=kx+b; in this application, b=0; Calculate the projection coordinates of all gray feature subset points in the training sample set on the straight line y=kx+b, according to the point-to-straight projection formula:
Figure PCTCN2019125058-appb-000012
Figure PCTCN2019125058-appb-000012
其中,n=1,2,3,…,50;分别计算两类子集投影变换后点的横坐标均值,即:Among them, n = 1, 2, 3,..., 50; calculate the mean value of the abscissa of the points after the projection transformation of the two types of subsets, namely:
Figure PCTCN2019125058-appb-000013
Figure PCTCN2019125058-appb-000013
Figure PCTCN2019125058-appb-000014
Figure PCTCN2019125058-appb-000014
此时,判别函数直线的阈值
Figure PCTCN2019125058-appb-000015
其中,k=W(2)/W(1)。
At this time, the threshold of the discriminant function line
Figure PCTCN2019125058-appb-000015
Among them, k=W(2)/W(1).
步骤S500,根据灰度特征子集投影后的纵坐标值与预设判别函数直线的阈值的大小关系,得到待检测设备的图像对应的设备的缺陷状态。Step S500: Obtain the defect state of the device corresponding to the image of the device to be tested according to the magnitude relationship between the ordinate value after the projection of the gray feature subset and the threshold value of the preset discriminant function line.
如上述实施例所述,当确定好判别函数直线后,将测试样本集的灰度特征子集即T={(m T,σ T)}投影到该判别函数直线上,可得灰度特征子集的投影后的坐标值为: As described in the above embodiment, when the line of the discriminant function is determined, the gray feature subset of the test sample set, namely T={(m TT )}, is projected onto the line of the discriminant function to obtain the gray feature The coordinate value after projection of the subset is:
Figure PCTCN2019125058-appb-000016
Figure PCTCN2019125058-appb-000016
本实施例中,检测孔位内保护门的状态即检测保护门的复位与无复位两种状态,因此是一个二分类问题,可记其输出标记Result∈{0,1}。本实施例中,通过分类规则,比较y T值与阈值w 0的大小,得出其分类,分类规则如下: In this embodiment, detecting the status of the protection door in the hole position is to detect the reset and no reset status of the protection door. Therefore, it is a two-class problem, and its output mark Result ∈ {0, 1} can be recorded. In this embodiment, the classification rules are used to compare the value of y T with the threshold w 0 to obtain the classification. The classification rules are as follows:
Figure PCTCN2019125058-appb-000017
Figure PCTCN2019125058-appb-000017
即若Result输出等于0,则代表此孔位内保护门没有复位;若Result输出等于1,则代表此孔位孔内保护门复位。That is, if the Result output is equal to 0, it means that the protection door in this hole position is not reset; if the Result output is equal to 1, it means that the protection door in this hole position is reset.
上述孔位保护门状态检测方法,从图像识别角度出发,首先获取待测设备图像构建训练样本集,根据训练样本集得到预设的判别函数直线的阈值,然后,对新获取待检测设备的图像进行基于局部统计的自适应动态阈值分割,能够精确地筛选出孔位区域,通过将提取出的孔位区域的灰度特征投影到预设判别函数直线,来达到抽取分类信息和压缩特征空间维数的效果,投影后使得特征子集在投影后的子空间有最大的类间距离和最小的类内距离,保证样本在该空间有最佳可分离性,最后,比较投影后的特征子集的纵坐标值与预设判别函数直线的阈值,即能判断出待测设备是否存在缺陷,上述方法便捷有效准确率高, 且抗干扰能力强检测结果稳定。The above hole position protection door status detection method, from the perspective of image recognition, first obtains the image of the device to be tested to construct a training sample set, obtains the preset threshold of the discriminant function line according to the training sample set, and then obtains the image of the device to be tested. The adaptive dynamic threshold segmentation based on local statistics can accurately filter out the hole location area. By projecting the gray-scale feature of the extracted hole location area to the preset discriminant function line, it can extract the classification information and compress the feature space dimension. After the projection, the feature subset has the largest inter-class distance and the smallest intra-class distance in the projected subspace to ensure the best separability of the sample in the space. Finally, compare the projected feature subsets The ordinate value of and the threshold value of the preset discriminant function straight line can determine whether the device under test has defects. The above method is convenient, effective, and accurate, and has strong anti-interference ability and stable detection results.
如图5所示,在其中一个实施例中,对待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域包括:As shown in FIG. 5, in one of the embodiments, performing adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected, and determining the hole location area includes:
步骤S202,根据预设的基于局部统计的自适应动态阈值分割法,对待检测设备的图像进行阈值分割,得到缺陷像素点集合;Step S202, according to a preset adaptive dynamic threshold segmentation method based on local statistics, threshold segmentation is performed on the image of the equipment to be inspected to obtain a set of defective pixels;
步骤S204,根据connection算子,对缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合;Step S204: Perform connected domain analysis on the set of defective pixels according to the connection operator to obtain a set of suspicious defective regions;
步骤S206,根据预设的孔位面积约束条件以及预设的孔位蓬松度,从可疑缺陷区域集合筛选出孔位区域。Step S206: According to the preset hole location area constraint condition and the preset hole location bulkiness, the hole location area is screened out from the set of suspicious defect areas.
连通区域一般是指图像中具有相同像素值且位置相邻的前景像素点组成的图像区域。连通区域分析是指将图像中的各个连通区域找出并标记。本实施例中,在加入滤波器窗口的标准方差修正偏移量Offset后,可基于偏移量Offset、滤波器窗口中所有像素点的像素平均值、预设的标准方差比例因子和最小绝对阈值,对待检测设备的图像进行阈值分割,分割出缺陷点集合,缺陷点可看作一个个独立的小区域,在分割出缺陷点后,可采用connection算子,对缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合,所谓可疑缺陷区域集合可理解为可能存在缺陷的孔位区域。由于待检测设备的孔位区域大小是相同且固定的,故可通过孔位的面积特征进一步筛选孔位区域。具体的,可以是通过孔位的像素面积作为约束条件。待检测设设备以插座为例,插座包括两孔插座和三孔插座,一般的,两孔插座的插孔的物理尺寸大小为0.9cm*0.9cm,三孔插座的孔位的物理尺寸大小为0.9cm*0.2cm。从图像角度上说,采集到的待检测设备的图像上,三孔位的每个孔位面积大约为14000像素大小,两孔位的每个孔位面积大约为72800像素大小。因此,本实施例中设定孔位的像素面积在[13000,5000]与[70000,75000]两个范围内为约束条件,基于此,可以从可疑缺陷区域集合中筛选出包含孔位区域集合。对于分割后的二值化待检测设备的图像而言,孔位区域面积S可由如下公式表述:The connected area generally refers to an image area composed of foreground pixels that have the same pixel value and are located adjacent to each other. Connected area analysis refers to finding and marking each connected area in an image. In this embodiment, after adding the standard deviation of the filter window to correct the offset Offset, it may be based on the offset Offset, the pixel average of all pixels in the filter window, the preset standard deviation scale factor and the minimum absolute threshold. , Perform threshold segmentation on the image of the equipment to be inspected to segment the set of defect points. The defect points can be regarded as independent small areas. After the defect points are segmented, the connection operator can be used to analyze the connected domain of the defective pixel point set , Get the suspicious defect area set, the so-called suspicious defect area set can be understood as the hole location area that may have defects. Since the size of the hole area of the device to be tested is the same and fixed, the hole area can be further screened by the area characteristics of the hole. Specifically, the area of the pixel passing through the hole may be used as a constraint condition. Take the socket as an example for the equipment to be tested. The socket includes a two-hole socket and a three-hole socket. Generally, the physical size of the socket of the two-hole socket is 0.9cm*0.9cm, and the physical size of the hole of the three-hole socket is 0.9cm*0.2cm. From an image perspective, in the collected images of the equipment to be tested, the area of each of the three holes is about 14,000 pixels, and the area of each of the two holes is about 72,800 pixels. Therefore, in this embodiment, setting the pixel area of the hole position within the two ranges of [13000, 5000] and [70000, 75000] is a constraint condition. Based on this, a set of regions containing hole positions can be filtered from a set of suspicious defect areas . For the segmented binarized image of the equipment to be tested, the area S of the hole location area can be expressed by the following formula:
Figure PCTCN2019125058-appb-000018
Figure PCTCN2019125058-appb-000018
具体的,假设所有满足约束条件的区域集合表示为:A={S i|13000<S i<15000∪70000<S i<75000,S i∈S}。如果只采用单一的面积特征进行筛选,很可能会将非孔位区域判断为孔位区域。因此,提出加入孔位的外观形状特征即孔位区域的蓬松度作为附加的筛选条件。蓬松度(Bulkiness): Specifically, it is assumed that the set of all regions that meet the constraint conditions is expressed as: A={S i |13000<S i <15000∪70000<S i <75000, S i εS}. If only a single area feature is used for screening, it is likely to judge the non-hole area as a hole area. Therefore, it is proposed to add the appearance and shape feature of the hole location, that is, the bulkiness of the hole location area as an additional screening condition. Bulkiness:
Bulkiness=π*R a*R b/A i Bulkiness=π*R a *R b /A i
其中Ra和Rb表示最小外接椭圆半径,A表示区域的面积。当蓬松度越接近1时,代表此区域的最小外接椭圆越接近圆形。Where Ra and Rb represent the minimum circumscribed ellipse radius, and A represents the area of the region. When the bulkiness is closer to 1, the smallest circumscribed ellipse representing this area is closer to a circle.
待检测设备以插座为例,通过多个实验样品统计分析计算得到插孔区域的蓬松度都在1.05左右,与其它区域的蓬松度值具有很明显的差异性。因此,本实施例中,设定插孔蓬松度处于[1.04,1.06]范围为约束条件,最终可以筛选出所有孔位区域的集合D:Taking the socket as an example of the equipment to be tested, the bulkiness of the jack area is about 1.05 calculated through statistical analysis of multiple experimental samples, which is obviously different from the bulkiness value of other areas. Therefore, in this embodiment, setting the hole bulkiness in the range of [1.04, 1.06] is the constraint condition, and finally the set D of all hole location regions can be filtered out:
D={A i|1.04<π*R a*R b/A i<1.06,A i∈A} D={A i |1.04<π*R a *R b /A i <1.06, A i ∈A}
可以理解的是,连通域标记方法除connection算子之外,还可以采用其他连通域分析算法。蓬松度的约束条件可以其他数值的约束条件,具体可视情况而定。本实施例中,通过预设的孔位蓬松度的约束条件,能够提高孔位区域筛选精度。It is understandable that in addition to the connection operator, other connected domain analysis algorithms can also be used in the connected domain labeling method. The constraint condition of the bulkiness can be the constraint condition of other values, and the specific conditions may be determined. In this embodiment, the preset hole position bulkiness constraint condition can improve the accuracy of hole position area screening.
在其中一个实施例中,根据预设的基于局部统计的自适应动态阈值分割法,对待检测设备的图像进行阈值分割,得到缺陷像素点集合包括:基于预设的滤波器,对待检测设备的图像进行平滑滤波处理,根据预设滤波器的窗口中的像素灰度平均值,计算待检测设备的图像的局部标准差,根据局部标准差、预设的标准方差比例因子以及预设的最小绝对阈值,计算各像素点的修正偏移量,获取待检测设备的图像中各像素点的灰度值以及各像素点的模板邻域的像素灰度平均值,根据各像素点的灰度值、各像素点的模板邻域的像素灰度平均值以及各像素点的修正偏移量,得到缺陷像素点集合。In one of the embodiments, according to a preset adaptive dynamic threshold segmentation method based on local statistics, threshold segmentation of the image of the device to be inspected to obtain a set of defective pixels includes: based on a preset filter, the image of the device to be inspected Perform smoothing filter processing, calculate the local standard deviation of the image of the device to be tested according to the average value of the pixel gray in the window of the preset filter, according to the local standard deviation, the preset standard deviation scale factor and the preset minimum absolute threshold , Calculate the correction offset of each pixel, obtain the gray value of each pixel in the image of the device to be tested and the average pixel gray value of the template neighborhood of each pixel, according to the gray value of each pixel, each pixel The pixel gray average value of the template neighborhood of the pixel point and the corrected offset of each pixel point are obtained to obtain the defective pixel point set.
本实施例中,预设的滤波器可以是均值滤波器、二项式滤波器及高斯滤波器等,滤波器的窗口大小为150像素*150像素,局部标准差即滤波器窗口的标准方差统计信息,它反映了一幅图像当中局部区域对比度的变化。假设f(i,j)为待检测设备的图像F中(i,j)位置处像素的灰度值,图像大小为M×N,W为以(i,j)为中心,大小为l×l的窗口,l为奇数且l>1,W看作是图像中的局部区域,计算其局部标准差Dev ij,定义如下: In this embodiment, the preset filter can be a mean filter, a binomial filter, a Gaussian filter, etc. The window size of the filter is 150 pixels * 150 pixels, and the local standard deviation is the standard deviation statistics of the filter window. Information, it reflects the changes in the contrast of a local area in an image. Suppose f(i,j) is the gray value of the pixel at the position (i,j) in the image F of the device to be tested, the image size is M×N, W is centered on (i,j), and the size is l× The window of l, l is odd and l>1, W is regarded as a local area in the image, and its local standard deviation Dev ij is calculated, which is defined as follows:
Figure PCTCN2019125058-appb-000019
Figure PCTCN2019125058-appb-000019
式中,i,j≥0,m≤M-1,n≤N-1,sqrt()为开方运算。设g(x,y)为输入的待检测设备的图像在(x,y)处像素的灰度值,m(x,y)为模板邻域内的平均灰度值信息,d(x,y)为局部标准差,则该像素点的统计修正偏移量为:In the formula, i, j≥0, m≤M-1, n≤N-1, sqrt() is the square root operation. Let g(x,y) be the gray value of the pixel at (x,y) of the input image of the device to be tested, m(x,y) is the average gray value information in the neighborhood of the template, d(x,y ) Is the local standard deviation, then the statistical correction offset of the pixel is:
Figure PCTCN2019125058-appb-000020
Figure PCTCN2019125058-appb-000020
式中,std为标准方差比例因子,T abs是最小绝对阈值,本实施例中,设置标准方差比例因子std为1.36,T abs最小绝对阈值为20。根据平均值m(x,y)和偏移修正量Offset,从待检测设备的图像中通过阈值分割出所有缺陷像素点的集合D,集合D可以表示为: In the formula, std is the standard deviation scale factor, and T abs is the minimum absolute threshold. In this embodiment, the standard deviation scale factor std is set to 1.36, and the minimum absolute threshold of T abs is 20. According to the average value m(x, y) and the offset correction amount Offset, a set D of all defective pixels is segmented from the image of the device to be tested through a threshold. The set D can be expressed as:
D={(x,y)|g(x,y)<m(x,y)-Offset∪g(x,y)>m(x,y)+Offset}D={(x,y)|g(x,y)<m(x,y)-Offset∪g(x,y)>m(x,y)+Offset}
可以理解的是,在其他实施例中,std以及T abs还可以是其他数值,具体可根据实际情况而定。本实施例中,通过计算滤波器窗口的局部标准差来修正偏移量Offset,能够限定 灰度值的合理变化范围,避免小区域的缺陷被过度平滑导致无法正确分割,提高缺陷分割的完整性。 It can be understood that in other embodiments, std and T abs may also be other values, which may be determined according to actual conditions. In this embodiment, the offset Offset is corrected by calculating the local standard deviation of the filter window, which can limit the reasonable change range of the gray value, avoid the defects in small areas from being excessively smoothed and lead to incorrect segmentation, and improve the integrity of defect segmentation .
在其中一个实施例中,根据预设滤波器窗口中的像素灰度平均值,计算待检测设备的图像的局部标准差包括:根据预设滤波器的窗口中的中心像素点灰度值,得到预设滤波器的窗口中的所有像素点的灰度平均值,根据中心像素点灰度值以及灰度平均值,计算局部标准差。In one of the embodiments, calculating the local standard deviation of the image of the device to be tested according to the average value of the pixel gray levels in the preset filter window includes: obtaining the gray value of the central pixel point in the window of the preset filter The gray average value of all pixels in the window of the preset filter is calculated based on the gray value of the center pixel and the gray average value to calculate the local standard deviation.
Figure PCTCN2019125058-appb-000021
为窗口W中所有像素灰度平均值,根据预设滤波器的窗口中的中心像素点灰度值,则可定义
Figure PCTCN2019125058-appb-000022
如下:
If
Figure PCTCN2019125058-appb-000021
It is the average gray value of all pixels in window W, and the gray value of the central pixel in the window of the preset filter can be defined
Figure PCTCN2019125058-appb-000022
as follows:
Figure PCTCN2019125058-appb-000023
Figure PCTCN2019125058-appb-000023
式中,i,j≥0,m≤M-1,n≤N-1。本实施例中,基于像素灰度平均值和中心像素点灰度值计算局部标准差,能够在检测孤立点边缘降低敏感度,具有一定的平滑滤波的效果。In the formula, i, j≥0, m≤M-1, n≤N-1. In this embodiment, the local standard deviation is calculated based on the pixel gray average value and the center pixel gray value, which can reduce the sensitivity in detecting the edge of the isolated point, and has a certain smoothing filtering effect.
在其中一个实施例中,得到待检测设备的图像对应的设备的缺陷状态之后,还包括:根据gen_contour_region_xld算子,提取孔位区域的边缘轮廓,根据缺陷状态,标记边缘轮廓,推送标记后的边缘轮廓的图像。In one of the embodiments, after obtaining the defect state of the device corresponding to the image of the device to be inspected, the method further includes: extracting the edge contour of the hole region according to the gen_contour_region_xld operator, marking the edge contour according to the defect state, and pushing the marked edge Contour image.
在实际应用中,待检测设备的图像的处理可以是基于Halcon 17.12 progress平台处理的,为便于工作人员的检测,待检测设备的图像的检测结果可显示于软件界面。具体的,可以是基于gen_contour_region_xld算子,提取孔位区域的边缘轮廓,gen_contour_region_xld算子为Halcon算子函数的一种,其用于根据区域创建XLD轮廓(contour)。本实施例中,可以是当Result(缺陷状态结果)输出为0,则代表此插孔内保护门没有复位,此时,则用红色轮廓标记孔位的边缘轮廓,若Result(缺陷状态结果)输出等于1,则代表此孔位内保护门复位,则用绿色轮廓标记孔位的边缘轮廓,如此,工作人员能够通过显示出的孔位轮廓的颜色,来判断待检测设备的保护门是否复位,便捷有效。可以理解的是,在其他实施例中,孔位的边缘轮廓的方式还可以是用其他颜色标记,或通过其他能够区分孔位保护门是否复位的方式进行标记。In practical applications, the processing of the image of the device to be tested can be based on the Halcon 17.12 progress platform. To facilitate the detection by the staff, the detection result of the image of the device to be tested can be displayed on the software interface. Specifically, it may be based on the gen_contour_region_xld operator to extract the edge contour of the hole location region. The gen_contour_region_xld operator is a type of Halcon operator function, which is used to create an XLD contour (contour) according to the region. In this embodiment, when the output of Result (defect status result) is 0, it means that the protection door in the jack has not been reset. At this time, the edge contour of the hole position is marked with a red outline, if Result (defect status result) If the output is equal to 1, it means that the protection door in this hole position is reset, and the edge contour of the hole position is marked with a green outline, so that the staff can judge whether the protection door of the equipment to be tested is reset by the color of the displayed hole position outline , Convenient and effective. It can be understood that, in other embodiments, the edge contour of the hole position may also be marked with other colors, or marked by other methods that can distinguish whether the hole position protection door is reset.
为清楚地描述本申请提供的孔位保护门状态检测方法,下面结合一个具体实例以及图2和图3来进行说明,其中,待检测设备101以插座为例,摄像装置103以相机为例:In order to clearly describe the method for detecting the status of the hole protection door provided by the present application, a specific example and FIG. 2 and FIG. 3 are combined for description. The equipment to be tested 101 is an example of a socket, and the imaging device 103 is an example of a camera:
当处于流水线上的正面朝上的待检测设备101即插座(插座带插孔那一面为正面)通过光电传感器105时,光电传感器105将触发信号发送给工控机104,工控机104收到信号后通过IO板卡输出控制信号,首先打开置于插座正上方的环形无影光源,然后给相机103(相机型号为BasleracA2500-14gm,水平/垂直分辨率为2592pixel*1944pixel)发送软触发信号,采集到插座正面图像,然后相机103将采集的图像上传至服务器,服务器获取插座图像,基于Halcon 17.12 progress平台对插座图像进行基于局部统计的自适应动态阈值分割,获取插孔区域。其中,插孔区域包括基于插孔的像素面积约束条件和蓬松度要求获取,即插孔 的像素面积处于[13000,15000]或[70000,75000]之间,插孔蓬松度处于[1.04,1.06]范围之间。然后,基于灰度直方图的统计算法,提取出插孔区域的灰度特征子集T={(m T,σ T)},将灰度特征子集T={(m T,σ T)}投影到二维坐标系中预设的判别函数直线y=kx+b上,获取灰度特征子集的投影纵坐标值y T,根据灰度特征子集投影纵坐标值y T与预设判别函数直线的阈值w 0的大小关系,得到插座图像对应的插座的缺陷状态,若y T≤w 0,则输出等于0,则代表插孔内保护门没有复位,则以红色标记插孔边缘轮廓;若y T>w 0,则输出等于1,则代表此孔位孔内保护门复位,以绿色标记插孔边缘轮廓。 When the device to be tested 101 on the assembly line facing up, that is, the socket (the side of the socket with the jack is the front), passes the photoelectric sensor 105, the photoelectric sensor 105 sends a trigger signal to the industrial computer 104, and the industrial computer 104 receives the signal Output the control signal through the IO board, first turn on the ring shadowless light source directly above the socket, and then send a soft trigger signal to the camera 103 (the camera model is BasleracA2500-14gm, the horizontal/vertical resolution is 2592pixel*1944pixel), and the The front image of the socket, and then the camera 103 uploads the collected image to the server, and the server obtains the socket image. Based on the Halcon 17.12 progress platform, the socket image is segmented by an adaptive dynamic threshold based on local statistics to obtain the socket area. Among them, the jack area includes the pixel area constraints and bulkiness requirements based on the jack, that is, the pixel area of the jack is between [13000, 15000] or [70000, 75000], and the jack bulk is [1.04, 1.06] ] Range. Then, based on the statistical algorithm of the gray histogram, the gray feature subset T={(m TT )} of the jack area is extracted, and the gray feature subset T={(m TT ) } projected two-dimensional coordinate system preset on the discriminant function straight line y = kx + b, obtaining projection gradation ordinate value feature subset y T, set the value of ordinate y T projection with a predetermined gradation characteristic according to the sub- Determine the relationship between the threshold value w 0 of the straight line of the discriminant function, and obtain the defect status of the socket corresponding to the socket image. If y T ≤w 0 , the output is equal to 0, which means that the protection door in the socket is not reset, and the edge of the socket is marked in red Contour; if y T >w 0 , the output is equal to 1, which means that the protective door in the hole is reset, and the contour of the edge of the jack is marked in green.
应该理解的是,虽然图2、图4以及图5的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2、图4以及图5中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。It should be understood that although the steps in the flowcharts of FIGS. 2, 4, and 5 are displayed in sequence as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless specifically stated in this article, the execution of these steps is not strictly limited in order, and these steps can be executed in other orders. Moreover, at least part of the steps in Figure 2, Figure 4, and Figure 5 may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The order of execution of the sub-steps or stages is not necessarily performed sequentially, but may be performed alternately or alternately with other steps or at least a part of the sub-steps or stages of other steps.
在其中一个实施例中,如图6所示,提供了一种孔位保护门状态检测装置,包括:图像获取模块610、孔位区域获取模块620、灰度特征提取模块630、投影模块640、缺陷检测模块650以及预设阈值计算模块660,其中:In one of the embodiments, as shown in FIG. 6, a device for detecting the status of a hole position protection door is provided, which includes: an image acquisition module 610, a hole position area acquisition module 620, a gray feature extraction module 630, a projection module 640, The defect detection module 650 and the preset threshold calculation module 660, wherein:
图像获取模块610,用于读取待检测设备的图像,获取待检测设备的图像,待检测设备上设置有孔位,孔位上设置有孔位保护门。The image acquisition module 610 is used to read the image of the device to be tested and acquire the image of the device to be tested. The device to be tested is provided with a hole position and a hole position protection door is arranged on the hole position.
孔位区域获取模块620,用于对待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域。The hole location area acquisition module 620 is configured to perform adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area.
灰度特征提取模块630,用于基于灰度直方图的统计算法,提取出孔位区域的灰度特征子集。The gray-level feature extraction module 630 is used to extract a gray-level feature subset of the hole location area based on a statistical algorithm of the gray-level histogram.
投影模块640,用于将灰度特征子集投影到二维坐标系中预设的判别函数直线上,计算灰度特征子集的投影纵坐标值,预设的判别函数直线为使得灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线。The projection module 640 is used to project the gray feature subset onto the preset discriminant function line in the two-dimensional coordinate system, and calculate the projection ordinate value of the gray feature subset. The preset discriminant function line is such that the gray feature The subset has a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace.
缺陷检测模块650,用于根据灰度特征子集投影后的纵坐标值与预设判别函数直线的阈值的大小关系,得到待检测设备的图像对应的设备的缺陷状态。The defect detection module 650 is configured to obtain the defect state of the device corresponding to the image of the device to be detected according to the magnitude relationship between the ordinate value after the projection of the gray feature subset and the threshold of the preset discriminant function line.
预设阈值计算模块660,用于获取多张待检测设备的历史图像,对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域,基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,训练样本集至少包括两类样本子集,将各类样本子集投影到最佳鉴别矢量空间,分别计算各类样本子集投影后的横坐标均值,得到预设的判别函数直线的阈值,最佳鉴别矢量空间为基于Fisher线性判别分析算法和训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。The preset threshold calculation module 660 is used to obtain multiple historical images of the equipment to be tested, perform adaptive dynamic threshold segmentation based on local statistics on the historical images of the equipment to be tested, determine the hole location area, and use a statistical algorithm based on gray histograms, Extract the gray-scale feature subsets of multiple hole locations, construct a training sample set, the training sample set includes at least two types of sample subsets, project various sample subsets to the best discriminant vector space, and calculate various sample subsets respectively The average value of the abscissa after projection is used to obtain the preset threshold of the discriminant function line. The best discriminant vector space is the best projection direction based on the Fisher linear discriminant analysis algorithm and training sample set, and the Lagrange multiplier method is introduced.
如图7所示,在其中一个实施例中,孔位保护门状态检测装置还包括标记推送模块 670,用于根据gen_contour_region_xld算子,提取孔位区域的边缘轮廓,根据缺陷状态,标记边缘轮廓,推送标记后的边缘轮廓的图像。As shown in FIG. 7, in one of the embodiments, the hole position protection door state detection device further includes a mark pushing module 670, which is used to extract the edge contour of the hole position region according to the gen_contour_region_xld operator, and mark the edge contour according to the defect state, Push the image of the edge contour after marking.
在其中一个实施例中,孔位保护门状态检测装置还包括平滑滤波模块680,用于根据预设的滤波器,对孔位区域进行平滑滤波处理,模糊孔位区域中的纹理信息。In one of the embodiments, the device for detecting the status of the hole protection door further includes a smoothing filter module 680, configured to perform smoothing filtering processing on the hole location area according to a preset filter to blur the texture information in the hole location area.
在其中一个实施例中,孔位区域获取模块620还用于根据预设的基于局部统计的自适应动态阈值分割法,对待检测设备的图像进行阈值分割,得到缺陷像素点集合,根据connection算子,对缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合,根据预设的孔位面积约束条件以及预设的孔位蓬松度,从可疑缺陷区域集合筛选出孔位区域。In one of the embodiments, the hole location area acquisition module 620 is further configured to perform threshold segmentation on the image of the device to be tested according to a preset adaptive dynamic threshold segmentation method based on local statistics to obtain a set of defective pixels, according to the connection operator , Perform connected domain analysis on the set of defective pixels to obtain a set of suspicious defect areas. According to preset hole area constraint conditions and preset hole position bulkiness, the hole location area is screened out from the set of suspicious defect areas.
在其中一个实施例中,孔位区域获取模块620还用于基于预设的滤波器,对待检测设备的图像进行平滑滤波处理,根据预设滤波器的窗口中的像素灰度平均值,计算待检测设备的图像的局部标准差,根据局部标准差、预设的标准方差比例因子以及预设的最小绝对阈值,计算各像素点的修正偏移量,获取待检测设备的图像中各像素点的灰度值以及各像素点的模板邻域的像素灰度平均值,根据各像素点的灰度值、各像素点的模板邻域的像素灰度平均值以及各像素点的修正偏移量,得到缺陷像素点集合。In one of the embodiments, the hole location area acquisition module 620 is further configured to perform smoothing filtering processing on the image of the device to be detected based on a preset filter, and calculate the average value of the pixel gray level in the window of the preset filter. The local standard deviation of the image of the detection device is calculated, according to the local standard deviation, the preset standard deviation scale factor and the preset minimum absolute threshold, the corrected offset of each pixel is calculated, and the value of each pixel in the image of the device to be detected is obtained. The gray value and the average pixel gray value of each pixel’s template neighborhood are based on the gray value of each pixel, the average pixel gray value of each pixel’s template neighborhood, and the correction offset of each pixel. Obtain a collection of defective pixels.
在其中一个实施例中,孔位区域获取模块620还用于根据预设滤波器的窗口中的中心像素点灰度值,得到预设滤波器的窗口中的所有像素点的灰度平均值,根据中心像素点灰度值以及灰度平均值,计算局部标准差。In one of the embodiments, the hole location area obtaining module 620 is further configured to obtain the gray value average value of all pixels in the window of the preset filter according to the gray value of the center pixel in the window of the preset filter. Calculate the local standard deviation based on the gray value of the central pixel and the average gray value.
关于孔位保护门状态检测装置的具体限定可以参见上文中对于孔位保护门状态检测方法的限定,在此不再赘述。上述孔位保护门状态检测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。For the specific definition of the hole position protection door state detection device, please refer to the above definition of the hole position protection door state detection method, which will not be repeated here. Each module in the above-mentioned hole position protection door status detection device can be implemented in whole or in part by software, hardware and a combination thereof. The foregoing modules may be embedded in the form of hardware or independent of the processor in the computer device, or may be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the foregoing modules.
在其中一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机可读指令和数据库。该内存储器为非易失性存储介质中的操作系统和计算机可读指令的运行提供环境。该计算机设备的数据库用于存储待检测设备的图像等数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机可读指令被处理器执行时以实现一种孔位保护门状态检测方法。In one of the embodiments, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as shown in FIG. 8. The computer equipment includes a processor, a memory, a network interface and a database connected through a system bus. Among them, the processor of the computer device is used to provide calculation and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, computer readable instructions, and a database. The internal memory provides an environment for the operation of the operating system and computer-readable instructions in the non-volatile storage medium. The computer equipment database is used to store data such as images of the equipment to be tested. The network interface of the computer device is used to communicate with an external terminal through a network connection. When the computer readable instruction is executed by the processor, a method for detecting the status of the hole position protection door is realized.
本领域技术人员可以理解,图8中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。Those skilled in the art can understand that the structure shown in FIG. 8 is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may Including more or fewer parts than shown in the figure, or combining some parts, or having a different arrangement of parts.
一种计算机设备,包括存储器和一个或多个处理器,存储器中储存有计算机可读指令,计算机可读指令被处理器执行时,使得一个或多个处理器执行本申请任意一个实施例中提 供的孔位保护门状态检测方法的步骤。A computer device includes a memory and one or more processors. The memory stores computer-readable instructions. When the computer-readable instructions are executed by the processor, the one or more processors execute any one of the embodiments of the present application. The steps of the hole position protection door status detection method.
一个或多个存储有计算机可读指令的非易失性存储介质,计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行本申请任意一个实施例中提供的孔位保护门状态检测方法的步骤。One or more non-volatile storage media storing computer-readable instructions. When the computer-readable instructions are executed by one or more processors, the one or more processors execute the holes provided in any embodiment of the present application. The steps of the method for detecting the status of the bit protection gate.
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机可读指令来指令相关的硬件来完成,所述的计算机可读指令可存储于一非易失性计算机可读取存储介质中,该计算机可读指令在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing relevant hardware through computer-readable instructions, which can be stored in a non-volatile computer. In a readable storage medium, when the computer-readable instructions are executed, they may include the processes of the above-mentioned method embodiments. Wherein, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and/or volatile memory. Non-volatile memory may include read only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. As an illustration and not a limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous chain Channel (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, all possible combinations of the technical features in the above embodiments are not described. However, as long as there is no contradiction between the combinations of these technical features, they should It is considered as the range described in this specification.
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。The above-mentioned embodiments only express several implementation manners of the present application, and the description is relatively specific and detailed, but it should not be understood as a limitation on the scope of the invention patent. It should be pointed out that for those of ordinary skill in the art, without departing from the concept of this application, several modifications and improvements can be made, and these all fall within the protection scope of this application. Therefore, the scope of protection of the patent of this application shall be subject to the appended claims.

Claims (20)

  1. 一种孔位保护门状态检测方法,所述方法包括:A method for detecting the status of a hole position protection door, the method comprising:
    获取待检测设备的图像,所述待检测设备上设置有孔位,所述孔位上设置有孔位保护门;Acquiring an image of the device to be tested, the device to be tested is provided with a hole position, and the hole position is provided with a hole position protection door;
    对所述待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Performing adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area;
    基于灰度直方图的统计算法,提取出所述孔位区域的灰度特征子集;Based on the statistical algorithm of the gray histogram, extract the gray feature subset of the hole location area;
    将所述灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取所述灰度特征子集的投影纵坐标值,所述预设的判别函数直线为使得所述灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;及Project the gray-scale feature subset onto a preset discriminant function line in a two-dimensional coordinate system to obtain the projected ordinate value of the gray-scale feature subset, and the preset discriminant function line is such that the gray The degree feature subset has a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace; and
    根据所述灰度特征子集投影后的纵坐标值与所述预设的判别函数直线的阈值的大小关系,得到所述待检测设备的图像对应的设备的缺陷状态;Obtaining the defect state of the device corresponding to the image of the device to be tested according to the magnitude relationship between the ordinate value after the projection of the gray-scale feature subset and the preset threshold of the discriminant function line;
    所述预设的判别函数直线的阈值由以下方式得到:The preset threshold of the discriminant function straight line is obtained in the following manner:
    获取多张待检测设备的历史图像;Obtain multiple historical images of the equipment to be tested;
    对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Perform adaptive dynamic threshold segmentation based on local statistics on the historical image of the equipment to be tested to determine the hole location area;
    基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,所述训练样本集至少包括两类样本子集;及Based on the statistical algorithm of the gray histogram, extract the gray feature subsets of multiple hole location regions to construct a training sample set, the training sample set includes at least two types of sample subsets; and
    将各类样本子集投影到最佳鉴别矢量空间,分别计算所述各类样本子集投影后的横坐标均值,得到所述预设的判别函数直线的阈值,所述最佳鉴别矢量空间为基于Fisher线性判别分析算法和所述训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。Project various sample subsets into the optimal discriminative vector space, and calculate the mean values of the abscissas after the projection of the various sample subsets to obtain the preset threshold value of the discriminant function straight line, and the optimal discriminant vector space is Based on the Fisher linear discriminant analysis algorithm and the training sample set, the optimal projection direction obtained by the Lagrangian multiplier method is introduced.
  2. 根据权利要求1所述的孔位保护门状态检测方法,其特征在于,所述对所述待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域包括:The method for detecting hole position protection door status according to claim 1, wherein said performing adaptive dynamic threshold segmentation based on local statistics on the image of said device to be detected, and determining the hole position area comprises:
    根据预设的基于局部统计的自适应动态阈值分割法,对所述待检测设备的图像进行阈值分割,得到缺陷像素点集合;Perform threshold segmentation on the image of the device to be tested according to a preset adaptive dynamic threshold segmentation method based on local statistics to obtain a set of defective pixels;
    根据connection算子,对所述缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合;及According to the connection operator, perform connected domain analysis on the set of defective pixels to obtain a set of suspicious defective regions; and
    根据预设的孔位面积约束条件以及预设的孔位蓬松度,从所述可疑缺陷区域集合筛选出孔位区域。According to preset hole location area constraints and preset hole location bulkiness, the hole location area is screened out from the set of suspicious defect areas.
  3. 根据权利要求2所述的孔位保护门状态检测方法,其特征在于,所述根据预设的基于局部统计的自适应动态阈值分割法,对所述待检测设备的图像进行阈值分割,得到缺陷像素点集合包括:The hole position protection door status detection method according to claim 2, wherein the image of the device to be detected is thresholded according to a preset adaptive dynamic threshold segmentation method based on local statistics to obtain defects The pixel point collection includes:
    基于预设的滤波器,对所述待检测设备的图像进行平滑滤波处理,根据所述预设滤波器的窗口中的像素灰度平均值,计算所述待检测设备的图像的局部标准差;Performing smoothing filtering processing on the image of the device to be detected based on the preset filter, and calculating the local standard deviation of the image of the device to be detected according to the average value of the pixel gray levels in the window of the preset filter;
    根据所述局部标准差、预设的标准方差比例因子以及预设的最小绝对阈值,计算各像素点的修正偏移量;Calculating the correction offset of each pixel according to the local standard deviation, the preset standard deviation scale factor and the preset minimum absolute threshold;
    获取所述待检测设备的图像中各像素点的灰度值以及各像素点的模板邻域的像素灰 度平均值;及Acquiring the gray value of each pixel in the image of the device to be tested and the average value of the pixel gray in the template neighborhood of each pixel; and
    根据所述各像素点的灰度值、各像素点的模板邻域的像素灰度平均值以及所述各像素点的修正偏移量,得到所述缺陷像素点集合。The set of defective pixels is obtained according to the gray value of each pixel, the average value of the pixel gray of the template neighborhood of each pixel, and the corrected offset of each pixel.
  4. 根据权利要求3所述的孔位保护门状态检测方法,其特征在于,所述根据所述预设滤波器窗口中的像素灰度平均值,计算所述待检测设备的图像的局部标准差包括:The method for detecting the status of the hole protection door according to claim 3, wherein the calculating the local standard deviation of the image of the device to be detected according to the average value of the pixel gray levels in the preset filter window comprises :
    根据所述预设滤波器的窗口中的中心像素点灰度值,得到所述预设滤波器的窗口中的所有像素点的灰度平均值;及Obtaining the average gray value of all pixels in the window of the preset filter according to the gray value of the central pixel in the window of the preset filter; and
    根据所述中心像素点灰度值以及所述灰度平均值,计算所述局部标准差。The local standard deviation is calculated according to the gray value of the central pixel point and the gray average value.
  5. 根据权利要求1所述的孔位保护门状态检测方法,其特征在于,在所述提取出所述孔位区域的灰度特征子集之前,所述方法还包括:The method for detecting the status of the hole position protection door according to claim 1, wherein before said extracting the gray feature subset of the hole position area, the method further comprises:
    根据预设的滤波器,对所述孔位区域进行平滑滤波处理,模糊所述孔位区域中的纹理信息。According to a preset filter, smoothing filtering is performed on the hole location area to blur the texture information in the hole location area.
  6. 根据权利要求1所述的孔位保护门状态检测方法,其特征在于,在所述得到所述待检测设备的图像对应的设备的缺陷状态之后,所述方法还包括:The method for detecting the status of the hole protection door according to claim 1, wherein after obtaining the defect status of the device corresponding to the image of the device to be detected, the method further comprises:
    根据gen_contour_region_xld算子,提取所述孔位区域的边缘轮廓;Extract the edge contour of the hole region according to the gen_contour_region_xld operator;
    根据所述缺陷状态,标记所述边缘轮廓;及Mark the edge contour according to the defect state; and
    推送标记后的所述边缘轮廓的图像。Push the marked image of the edge contour.
  7. 一种孔位保护门状态检测装置,其特征在于,所述装置包括:A device for detecting the status of a hole position protection door, wherein the device includes:
    图像获取模块,用于读取待检测设备的图像,获取待检测设备的图像,所述待检测设备上设置有孔位,所述孔位上设置有孔位保护门;The image acquisition module is used to read the image of the device to be tested and acquire the image of the device to be tested, the device to be tested is provided with a hole, and the hole is provided with a hole protection door;
    孔位区域获取模块,用于对所述待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;A hole location area acquisition module, configured to perform adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area;
    灰度特征提取模块,用于基于灰度直方图的统计算法,提取出所述孔位区域的灰度特征子集;The gray-level feature extraction module is used to extract the gray-level feature subset of the hole location area based on the statistical algorithm of the gray-level histogram;
    投影模块,用于将所述灰度特征子集投影到二维坐标系中预设的判别函数直线上,计算所述灰度特征子集的投影纵坐标值,所述预设的判别函数直线为使得所述灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;The projection module is used to project the gray-level feature subset onto a preset discriminant function line in a two-dimensional coordinate system, calculate the projected ordinate value of the gray-scale feature subset, and the preset discriminant function line So that the gray feature subset has a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace;
    缺陷检测模块,用于根据所述灰度特征子集投影后的纵坐标值与所述预设判别函数直线的阈值的大小关系,得到所述待检测设备的图像对应的设备的缺陷状态;及The defect detection module is used to obtain the defect state of the device corresponding to the image of the device to be detected according to the magnitude relationship between the ordinate value after the projection of the gray-scale feature subset and the threshold value of the preset discriminant function line; and
    预设阈值计算模块,用于获取多张待检测设备的历史图像,对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域,基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,所述训练样本集至少包括两类样本子集,将各类样本子集投影到最佳鉴别矢量空间,分别计算所述各类样本子集投影后的横坐标均值,得到所述预设的判别函数直线的阈值,所述最佳鉴别矢量空间为基于Fisher线性判别分析算法和所述训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。The preset threshold calculation module is used to obtain multiple historical images of the equipment to be inspected, perform adaptive dynamic threshold segmentation based on local statistics on the historical images of the equipment to be inspected, determine the hole location area, and extract the historical image based on the gray histogram Construct a training sample set of gray-scale feature subsets of a plurality of hole location regions. The training sample set includes at least two types of sample subsets. The various sample subsets are projected into the optimal discriminant vector space, and the various types of The average value of the abscissa after the projection of the sample subset is used to obtain the preset threshold value of the discriminant function straight line. The optimal discriminant vector space is based on the Fisher linear discriminant analysis algorithm and the training sample set, and Lagrange multipliers are introduced The best projection direction obtained by the method.
  8. 根据权利要求7所述的孔位保护门状态检测装置,其特征在于,所述装置还包括:The hole position protection door status detection device according to claim 7, wherein the device further comprises:
    标记推送模块,用于根据gen_contour_region_xld算子,提取所述孔位区域的边缘轮廓,根据所述缺陷状态,标记所述边缘轮廓,推送标记后的所述边缘轮廓的图像。The marking push module is used to extract the edge contour of the hole location area according to the gen_contour_region_xld operator, mark the edge contour according to the defect state, and push the marked image of the edge contour.
  9. 一种计算机设备,包括存储器及一个或多个处理器,所述存储器中储存有计算机可读指令,所述计算机可读指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行以下步骤:A computer device includes a memory and one or more processors. The memory stores computer-readable instructions. When the computer-readable instructions are executed by the one or more processors, the one or more Each processor performs the following steps:
    获取待检测设备的图像,所述待检测设备上设置有孔位,所述孔位上设置有孔位保护门;Acquiring an image of the device to be tested, the device to be tested is provided with a hole position, and the hole position is provided with a hole position protection door;
    对所述待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Performing adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area;
    基于灰度直方图的统计算法,提取出所述孔位区域的灰度特征子集;Based on the statistical algorithm of the gray histogram, extract the gray feature subset of the hole location area;
    将所述灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取所述灰度特征子集的投影纵坐标值,所述预设的判别函数直线为使得所述灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;及Project the gray-scale feature subset onto a preset discriminant function line in a two-dimensional coordinate system to obtain the projected ordinate value of the gray-scale feature subset, and the preset discriminant function line is such that the gray The degree feature subset has a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace; and
    根据所述灰度特征子集投影后的纵坐标值与所述预设的判别函数直线的阈值的大小关系,得到所述待检测设备的图像对应的设备的缺陷状态;Obtaining the defect state of the device corresponding to the image of the device to be tested according to the magnitude relationship between the ordinate value after the projection of the gray-scale feature subset and the preset threshold of the discriminant function line;
    所述预设的判别函数直线的阈值由以下方式得到:The preset threshold of the discriminant function straight line is obtained in the following manner:
    获取多张待检测设备的历史图像;Obtain multiple historical images of the equipment to be tested;
    对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Perform adaptive dynamic threshold segmentation based on local statistics on the historical image of the equipment to be tested to determine the hole location area;
    基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集,所述训练样本集至少包括两类样本子集;及Based on the statistical algorithm of the gray histogram, extract the gray feature subsets of multiple hole location regions to construct a training sample set, the training sample set includes at least two types of sample subsets; and
    将各类样本子集投影到最佳鉴别矢量空间,分别计算所述各类样本子集投影后的横坐标均值,得到所述预设的判别函数直线的阈值,所述最佳鉴别矢量空间为基于Fisher线性判别分析算法和所述训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。Project various sample subsets into the optimal discriminative vector space, and calculate the mean values of the abscissas after the projection of the various sample subsets to obtain the preset threshold value of the discriminant function straight line, and the optimal discriminant vector space is Based on the Fisher linear discriminant analysis algorithm and the training sample set, the optimal projection direction obtained by the Lagrangian multiplier method is introduced.
  10. 根据权利要求9所述的计算机设备,其特征在于,所述处理器执行所述计算机可读指令时还执行以下步骤:The computer device according to claim 9, wherein the processor further executes the following steps when executing the computer-readable instruction:
    根据预设的基于局部统计的自适应动态阈值分割法,对所述待检测设备的图像进行阈值分割,得到缺陷像素点集合;Perform threshold segmentation on the image of the device to be tested according to a preset adaptive dynamic threshold segmentation method based on local statistics to obtain a set of defective pixels;
    根据connection算子,对所述缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合;及According to the connection operator, perform connected domain analysis on the set of defective pixels to obtain a set of suspicious defective regions; and
    根据预设的孔位面积约束条件以及预设的孔位蓬松度,从所述可疑缺陷区域集合筛选出孔位区域。According to preset hole location area constraints and preset hole location bulkiness, the hole location area is screened out from the set of suspicious defect areas.
  11. 根据权利要求9所述的计算机设备,其特征在于,所述处理器执行所述计算机可读指令时还执行以下步骤:The computer device according to claim 9, wherein the processor further executes the following steps when executing the computer-readable instruction:
    基于预设的滤波器,对所述待检测设备的图像进行平滑滤波处理,根据所述预设滤波器的窗口中的像素灰度平均值,计算所述待检测设备的图像的局部标准差;Performing smoothing filtering processing on the image of the device to be detected based on the preset filter, and calculating the local standard deviation of the image of the device to be detected according to the average value of the pixel gray levels in the window of the preset filter;
    根据所述局部标准差、预设的标准方差比例因子以及预设的最小绝对阈值,计算各像素点的修正偏移量;Calculating the correction offset of each pixel according to the local standard deviation, the preset standard deviation scale factor and the preset minimum absolute threshold;
    获取所述待检测设备的图像中各像素点的灰度值以及各像素点的模板邻域的像素灰度平均值;及Obtaining the gray value of each pixel in the image of the device to be detected and the average value of the pixel gray of the template neighborhood of each pixel; and
    根据所述各像素点的灰度值、各像素点的模板邻域的像素灰度平均值以及所述各像素点的修正偏移量,得到所述缺陷像素点集合。The set of defective pixels is obtained according to the gray value of each pixel, the average value of the pixel gray of the template neighborhood of each pixel, and the corrected offset of each pixel.
  12. 根据权利要求9所述的计算机设备,其特征在于,所述处理器执行所述计算机可读指令时还执行以下步骤:The computer device according to claim 9, wherein the processor further executes the following steps when executing the computer-readable instruction:
    根据所述预设滤波器的窗口中的中心像素点灰度值,得到所述预设滤波器的窗口中的所有像素点的灰度平均值;及Obtaining the average gray value of all pixels in the window of the preset filter according to the gray value of the central pixel in the window of the preset filter; and
    根据所述中心像素点灰度值以及所述灰度平均值,计算所述局部标准差。The local standard deviation is calculated according to the gray value of the central pixel point and the gray average value.
  13. 根据权利要求9所述的计算机设备,其特征在于,所述处理器执行所述计算机可读指令时还执行以下步骤:The computer device according to claim 9, wherein the processor further executes the following steps when executing the computer-readable instruction:
    根据预设的滤波器,对所述孔位区域进行平滑滤波处理,模糊所述孔位区域中的纹理信息。According to a preset filter, smoothing filtering is performed on the hole location area to blur the texture information in the hole location area.
  14. 根据权利要求9所述的计算机设备,其特征在于,所述处理器执行所述计算机可读指令时还执行以下步骤:The computer device according to claim 9, wherein the processor further executes the following steps when executing the computer-readable instruction:
    根据gen_contour_region_xld算子,提取所述孔位区域的边缘轮廓;Extract the edge contour of the hole region according to the gen_contour_region_xld operator;
    根据所述缺陷状态,标记所述边缘轮廓;及Mark the edge contour according to the defect state; and
    推送标记后的所述边缘轮廓的图像。Push the marked image of the edge contour.
  15. 一个或多个存储有计算机可读指令的非易失性计算机可读存储介质,所述计算机可读指令被一个或多个处理器执行时,使得所述一个或多个处理器执行以下步骤:One or more non-volatile computer-readable storage media storing computer-readable instructions, which when executed by one or more processors, cause the one or more processors to perform the following steps:
    获取待检测设备的图像,所述待检测设备上设置有孔位,所述孔位上设置有孔位保护门;Acquiring an image of the device to be tested, the device to be tested is provided with a hole position, and the hole position is provided with a hole position protection door;
    对所述待检测设备的图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Performing adaptive dynamic threshold segmentation based on local statistics on the image of the device to be detected to determine the hole location area;
    基于灰度直方图的统计算法,提取出所述孔位区域的灰度特征子集;Based on the statistical algorithm of the gray histogram, extract the gray feature subset of the hole location area;
    将所述灰度特征子集投影到二维坐标系中预设的判别函数直线上,获取所述灰度特征子集的投影纵坐标值,所述预设的判别函数直线为使得所述灰度特征子集在投影后的子空间中有最大类间距离和最小类内距离的直线;及Project the gray-scale feature subset onto a preset discriminant function line in a two-dimensional coordinate system to obtain the projected ordinate value of the gray-scale feature subset, and the preset discriminant function line is such that the gray The degree feature subset has a straight line with the largest inter-class distance and the smallest intra-class distance in the projected subspace; and
    根据所述灰度特征子集投影后的纵坐标值与所述预设的判别函数直线的阈值的大小关系,得到所述待检测设备的图像对应的设备的缺陷状态;Obtaining the defect state of the device corresponding to the image of the device to be tested according to the magnitude relationship between the ordinate value after the projection of the gray-scale feature subset and the preset threshold of the discriminant function line;
    所述预设的判别函数直线的阈值由以下方式得到:The preset threshold of the discriminant function straight line is obtained in the following manner:
    获取多张待检测设备的历史图像;Obtain multiple historical images of the equipment to be tested;
    对待检测设备的历史图像进行基于局部统计的自适应动态阈值分割,确定孔位区域;Perform adaptive dynamic threshold segmentation based on local statistics on the historical image of the equipment to be tested to determine the hole location area;
    基于灰度直方图的统计算法,提取多个孔位区域的灰度特征子集,构建训练样本集, 所述训练样本集至少包括两类样本子集;及Based on the statistical algorithm of the gray histogram, extract the gray feature subsets of multiple hole location regions to construct a training sample set, the training sample set includes at least two types of sample subsets; and
    将各类样本子集投影到最佳鉴别矢量空间,分别计算所述各类样本子集投影后的横坐标均值,得到所述预设的判别函数直线的阈值,所述最佳鉴别矢量空间为基于Fisher线性判别分析算法和所述训练样本集,引入拉格朗日乘子法求解出的最佳投影方向。Project various sample subsets into the optimal discriminative vector space, and calculate the mean values of the abscissas after the projection of the various sample subsets to obtain the preset threshold value of the discriminant function straight line, and the optimal discriminant vector space is Based on the Fisher linear discriminant analysis algorithm and the training sample set, the optimal projection direction obtained by the Lagrangian multiplier method is introduced.
  16. 根据权利要求15所述的存储介质,其特征在于,所述计算机可读指令被所述处理器执行时还执行以下步骤:The storage medium according to claim 15, wherein the following steps are further executed when the computer-readable instructions are executed by the processor:
    根据预设的基于局部统计的自适应动态阈值分割法,对所述待检测设备的图像进行阈值分割,得到缺陷像素点集合;Perform threshold segmentation on the image of the device to be tested according to a preset adaptive dynamic threshold segmentation method based on local statistics to obtain a set of defective pixels;
    根据connection算子,对所述缺陷像素点集合进行连通域分析,得到可疑缺陷区域集合;及According to the connection operator, perform connected domain analysis on the set of defective pixels to obtain a set of suspicious defective regions; and
    根据预设的孔位面积约束条件以及预设的孔位蓬松度,从所述可疑缺陷区域集合筛选出孔位区域。According to preset hole location area constraints and preset hole location bulkiness, the hole location area is screened out from the set of suspicious defect areas.
  17. 根据权利要求15所述的存储介质,其特征在于,所述计算机可读指令被所述处理器执行时还执行以下步骤:The storage medium according to claim 15, wherein the following steps are further executed when the computer-readable instructions are executed by the processor:
    基于预设的滤波器,对所述待检测设备的图像进行平滑滤波处理,根据所述预设滤波器的窗口中的像素灰度平均值,计算所述待检测设备的图像的局部标准差;Performing smoothing filtering processing on the image of the device to be detected based on the preset filter, and calculating the local standard deviation of the image of the device to be detected according to the average value of the pixel gray levels in the window of the preset filter;
    根据所述局部标准差、预设的标准方差比例因子以及预设的最小绝对阈值,计算各像素点的修正偏移量;Calculating the correction offset of each pixel according to the local standard deviation, the preset standard deviation scale factor and the preset minimum absolute threshold;
    获取所述待检测设备的图像中各像素点的灰度值以及各像素点的模板邻域的像素灰度平均值;及Obtaining the gray value of each pixel in the image of the device to be detected and the average value of the pixel gray of the template neighborhood of each pixel; and
    根据所述各像素点的灰度值、各像素点的模板邻域的像素灰度平均值以及所述各像素点的修正偏移量,得到所述缺陷像素点集合。The set of defective pixels is obtained according to the gray value of each pixel, the average value of the pixel gray of the template neighborhood of each pixel, and the corrected offset of each pixel.
  18. 根据权利要求15所述的存储介质,其特征在于,所述计算机可读指令被所述处理器执行时还执行以下步骤:The storage medium according to claim 15, wherein the following steps are further executed when the computer-readable instructions are executed by the processor:
    根据所述预设滤波器的窗口中的中心像素点灰度值,得到所述预设滤波器的窗口中的所有像素点的灰度平均值;及Obtaining the average gray value of all pixels in the window of the preset filter according to the gray value of the central pixel in the window of the preset filter; and
    根据所述中心像素点灰度值以及所述灰度平均值,计算所述局部标准差。The local standard deviation is calculated according to the gray value of the central pixel point and the gray average value.
  19. 根据权利要求15所述的存储介质,其特征在于,所述计算机可读指令被所述处理器执行时还执行以下步骤:The storage medium according to claim 15, wherein the following steps are further executed when the computer-readable instructions are executed by the processor:
    根据预设的滤波器,对所述孔位区域进行平滑滤波处理,模糊所述孔位区域中的纹理信息。According to a preset filter, smoothing filtering is performed on the hole location area to blur the texture information in the hole location area.
  20. 根据权利要求15所述的存储介质,其特征在于,所述计算机可读指令被所述处理器执行时还执行以下步骤:The storage medium according to claim 15, wherein the following steps are further executed when the computer-readable instructions are executed by the processor:
    根据gen_contour_region_xld算子,提取所述孔位区域的边缘轮廓;Extract the edge contour of the hole region according to the gen_contour_region_xld operator;
    根据所述缺陷状态,标记所述边缘轮廓;及Mark the edge contour according to the defect state; and
    推送标记后的所述边缘轮廓的图像。Push the marked image of the edge contour.
PCT/CN2019/125058 2019-07-03 2019-12-13 Hole protection cap detection method and apparatus, computer device and storage medium WO2021000524A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201910593892.3 2019-07-03
CN201910593892.3A CN110930353B (en) 2019-07-03 2019-07-03 Method and device for detecting state of hole site protection door, computer equipment and storage medium

Publications (1)

Publication Number Publication Date
WO2021000524A1 true WO2021000524A1 (en) 2021-01-07

Family

ID=69855697

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/125058 WO2021000524A1 (en) 2019-07-03 2019-12-13 Hole protection cap detection method and apparatus, computer device and storage medium

Country Status (2)

Country Link
CN (1) CN110930353B (en)
WO (1) WO2021000524A1 (en)

Cited By (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112950633A (en) * 2021-04-19 2021-06-11 上海电机学院 Aluminum alloy weld surface defect detection method based on line structured light
CN113012124A (en) * 2021-03-15 2021-06-22 大连海事大学 Shoe mark hole and insert feature detection and description method
CN113284113A (en) * 2021-05-27 2021-08-20 西安闻泰信息技术有限公司 Glue overflow flaw detection method and device, computer equipment and readable storage medium
CN113362297A (en) * 2021-05-31 2021-09-07 北京百度网讯科技有限公司 Image processing method, apparatus and storage medium for printed circuit board
CN113393395A (en) * 2021-06-17 2021-09-14 西安应用光学研究所 High-dynamic infrared image segmentation threshold self-adaptive calculation method
CN113627523A (en) * 2021-08-09 2021-11-09 中国科学院微小卫星创新研究院 Method for detecting micro fault of satellite
CN113777033A (en) * 2021-08-18 2021-12-10 长沙长泰机器人有限公司 Raw sliver defect detection method and device based on machine vision
CN114140391A (en) * 2021-11-04 2022-03-04 东北大学 Method for realizing rapid detection of onboard display screen module based on machine vision
CN114324350A (en) * 2021-12-10 2022-04-12 巢湖学院 Intelligent five-hole socket panel defect detection method and system based on machine vision
CN114359597A (en) * 2021-11-26 2022-04-15 西安交通大学 Oil tank inner cover identification and pose parameter sensing method based on vision
CN114359598A (en) * 2021-11-26 2022-04-15 西安交通大学 Self-adaptive automobile oil tank outer cover identification method based on regional contrast difference
CN114813794A (en) * 2022-02-18 2022-07-29 成都飞机工业(集团)有限责任公司 3D printing nondestructive testing method based on robot
CN115019368A (en) * 2022-06-09 2022-09-06 南京审计大学 Face recognition feature extraction method in audit investigation based on 2DESDLPP
CN116109633A (en) * 2023-04-12 2023-05-12 山东金帝精密机械科技股份有限公司 Window detection method and device for bearing retainer
CN116152260A (en) * 2023-04-23 2023-05-23 广东工业大学 Spring defect detection method and system based on image processing
CN116188498A (en) * 2023-04-28 2023-05-30 江西科技学院 Axle welding area detection method and system based on computer vision
CN116309570A (en) * 2023-05-18 2023-06-23 山东亮马新材料科技有限公司 Titanium alloy bar quality detection method and system
CN116309552A (en) * 2023-05-12 2023-06-23 西南交通大学 Method, device, equipment and medium for evaluating health state of existing line old retaining wall
CN116363126A (en) * 2023-05-30 2023-06-30 东莞市立时电子有限公司 Welding quality detection method for data line USB plug
CN116363136A (en) * 2023-06-01 2023-06-30 山东创元智能设备制造有限责任公司 On-line screening method and system for automatic production of motor vehicle parts
CN116385439A (en) * 2023-06-05 2023-07-04 山东兰通机电有限公司 Motor rubber shock pad quality detection method based on image processing
CN116416246A (en) * 2023-06-08 2023-07-11 临沂中科芯华新材料科技有限公司 Machine vision-based fully-degradable plastic product film coating effect evaluation method
CN116433700A (en) * 2023-06-13 2023-07-14 山东金润源法兰机械有限公司 Visual positioning method for flange part contour
CN116433669A (en) * 2023-06-14 2023-07-14 山东兴华钢结构有限公司 Machine vision-based quality detection method for weld joints of steel frame of anti-seismic structure
CN116485832A (en) * 2023-06-25 2023-07-25 山东九思新材料科技有限责任公司 Method for accurately detecting edges of non-uniform fluid impurities for recycling waste silicon wafers
CN116542972A (en) * 2023-07-04 2023-08-04 山东阁林板建材科技有限公司 Wall plate surface defect rapid detection method based on artificial intelligence
CN116563223A (en) * 2023-04-11 2023-08-08 新创碳谷集团有限公司 Glass fiber yarn winding roller detection method, equipment and storage medium
CN116612126A (en) * 2023-07-21 2023-08-18 青岛国际旅行卫生保健中心(青岛海关口岸门诊部) Container disease vector biological detection early warning method based on artificial intelligence
CN116630308A (en) * 2023-07-20 2023-08-22 山东华太新能源电池有限公司 Data enhancement system for battery welding anomaly detection
CN116664554A (en) * 2023-07-26 2023-08-29 微山晟轩机械制造有限公司 Bolt thread defect detection method based on image processing
CN116681752A (en) * 2023-08-03 2023-09-01 山东墨氪智能科技有限公司 Method and device for calculating void ratio of void defects of DBC solder layer
CN116703251A (en) * 2023-08-08 2023-09-05 德润杰(山东)纺织科技有限公司 Rubber ring production quality detection method based on artificial intelligence
CN116703912A (en) * 2023-08-07 2023-09-05 深圳市鑫赛科科技发展有限公司 Mini-host network port integrity visual detection method
CN116777918A (en) * 2023-08-25 2023-09-19 苏州科尔珀恩机械科技有限公司 Visual auxiliary kiln surface defect rapid detection method
CN116823835A (en) * 2023-08-30 2023-09-29 山东省永星食品饮料有限公司 Bottled water impurity detection method based on machine vision
CN117102375A (en) * 2023-10-18 2023-11-24 沈阳欧施盾新材料科技有限公司 Special-shaped piece closing-in control method and equipment based on temperature imaging
CN117197534A (en) * 2023-08-04 2023-12-08 广州电缆厂有限公司 Automatic detection method for cable surface defects based on feature recognition
CN117197141A (en) * 2023-11-07 2023-12-08 山东远盾网络技术股份有限公司 Method for detecting surface defects of automobile parts
CN117274249A (en) * 2023-11-20 2023-12-22 江西省中鼐科技服务有限公司 Ceramic tile appearance detection method and system based on artificial intelligent image technology
CN117274247A (en) * 2023-11-20 2023-12-22 深圳市海里表面技术处理有限公司 Visual detection method for quality of LTCC conductor surface coating
CN117314949A (en) * 2023-11-28 2023-12-29 山东远硕上池健康科技有限公司 Personnel injury detection and identification method based on artificial intelligence
CN117372435A (en) * 2023-12-08 2024-01-09 智联信通科技股份有限公司 Connector pin detection method based on image characteristics
CN117455920A (en) * 2023-12-26 2024-01-26 武汉博源新材料科技集团股份有限公司 Artificial intelligence-based milk tea cup inferior product screening method and system
CN117576104A (en) * 2024-01-17 2024-02-20 山东世纪阳光科技有限公司 Visual detection method for health state of ultrafiltration membrane in purification process
CN117877008B (en) * 2024-03-13 2024-05-17 湖北神龙工程测试技术有限公司 Door and window performance detection method based on artificial intelligence

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114246356B (en) * 2020-09-25 2023-04-18 上海烟草集团有限责任公司 Design method, system, medium and device of cigarette leaf group formula
CN112734713B (en) * 2020-12-31 2024-02-09 西安交通大学 Transformer state detection method, device, computer equipment and storage medium
CN113654470A (en) * 2021-07-29 2021-11-16 佛山市科莱机器人有限公司 Method, system, equipment and medium for detecting accumulated water on surface of bathtub
CN114240928B (en) * 2021-12-29 2024-03-01 湖南云箭智能科技有限公司 Partition detection method, device and equipment for board quality and readable storage medium
CN114821073B (en) * 2022-06-28 2022-09-13 聊城市誉林工业设计有限公司 State identification method and device for portable intelligent shell opening machine

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6826300B2 (en) * 2001-05-31 2004-11-30 George Mason University Feature based classification
CN104361352A (en) * 2014-11-13 2015-02-18 东北林业大学 Solid wood panel defect separation method based on compressed sensing
CN106897994A (en) * 2017-01-20 2017-06-27 北京京仪仪器仪表研究总院有限公司 A kind of pcb board defect detecting system and method based on layered image

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435895B (en) * 2011-12-31 2014-07-23 罗格朗(北京)电气有限公司 Safety performance detection device, detection system and detection method for socket protection doors

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6826300B2 (en) * 2001-05-31 2004-11-30 George Mason University Feature based classification
CN104361352A (en) * 2014-11-13 2015-02-18 东北林业大学 Solid wood panel defect separation method based on compressed sensing
CN106897994A (en) * 2017-01-20 2017-06-27 北京京仪仪器仪表研究总院有限公司 A kind of pcb board defect detecting system and method based on layered image

Cited By (83)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113012124A (en) * 2021-03-15 2021-06-22 大连海事大学 Shoe mark hole and insert feature detection and description method
CN113012124B (en) * 2021-03-15 2024-02-23 大连海事大学 Shoe print hole and embedded object feature detection and description method
CN112950633B (en) * 2021-04-19 2023-06-23 上海电机学院 Aluminum alloy weld joint surface defect detection method based on line structured light
CN112950633A (en) * 2021-04-19 2021-06-11 上海电机学院 Aluminum alloy weld surface defect detection method based on line structured light
CN113284113A (en) * 2021-05-27 2021-08-20 西安闻泰信息技术有限公司 Glue overflow flaw detection method and device, computer equipment and readable storage medium
CN113362297A (en) * 2021-05-31 2021-09-07 北京百度网讯科技有限公司 Image processing method, apparatus and storage medium for printed circuit board
CN113362297B (en) * 2021-05-31 2023-09-19 北京百度网讯科技有限公司 Image processing method, apparatus and storage medium for printed circuit board
CN113393395A (en) * 2021-06-17 2021-09-14 西安应用光学研究所 High-dynamic infrared image segmentation threshold self-adaptive calculation method
CN113393395B (en) * 2021-06-17 2023-10-31 西安应用光学研究所 Adaptive calculation method for high-dynamic infrared image segmentation threshold
CN113627523B (en) * 2021-08-09 2024-03-26 中国科学院微小卫星创新研究院 Satellite micro fault detection method
CN113627523A (en) * 2021-08-09 2021-11-09 中国科学院微小卫星创新研究院 Method for detecting micro fault of satellite
CN113777033A (en) * 2021-08-18 2021-12-10 长沙长泰机器人有限公司 Raw sliver defect detection method and device based on machine vision
CN114140391A (en) * 2021-11-04 2022-03-04 东北大学 Method for realizing rapid detection of onboard display screen module based on machine vision
CN114359598B (en) * 2021-11-26 2023-09-12 西安交通大学 Self-adaptive automobile oil tank outer cover identification method based on region comparison difference
CN114359597B (en) * 2021-11-26 2023-09-12 西安交通大学 Vision-based oil tank inner cover identification and pose parameter sensing method
CN114359597A (en) * 2021-11-26 2022-04-15 西安交通大学 Oil tank inner cover identification and pose parameter sensing method based on vision
CN114359598A (en) * 2021-11-26 2022-04-15 西安交通大学 Self-adaptive automobile oil tank outer cover identification method based on regional contrast difference
CN114324350A (en) * 2021-12-10 2022-04-12 巢湖学院 Intelligent five-hole socket panel defect detection method and system based on machine vision
CN114813794A (en) * 2022-02-18 2022-07-29 成都飞机工业(集团)有限责任公司 3D printing nondestructive testing method based on robot
CN114813794B (en) * 2022-02-18 2023-10-03 成都飞机工业(集团)有限责任公司 Method for acquiring scanning photo required by robot 3D printing nondestructive testing
CN115019368A (en) * 2022-06-09 2022-09-06 南京审计大学 Face recognition feature extraction method in audit investigation based on 2DESDLPP
CN115019368B (en) * 2022-06-09 2023-09-12 南京审计大学 Face recognition feature extraction method in audit investigation
CN116563223A (en) * 2023-04-11 2023-08-08 新创碳谷集团有限公司 Glass fiber yarn winding roller detection method, equipment and storage medium
CN116563223B (en) * 2023-04-11 2023-09-26 新创碳谷集团有限公司 Glass fiber yarn winding roller detection method, equipment and storage medium
CN116109633A (en) * 2023-04-12 2023-05-12 山东金帝精密机械科技股份有限公司 Window detection method and device for bearing retainer
CN116152260A (en) * 2023-04-23 2023-05-23 广东工业大学 Spring defect detection method and system based on image processing
CN116152260B (en) * 2023-04-23 2023-08-18 广东工业大学 Spring defect detection method and system based on image processing
CN116188498A (en) * 2023-04-28 2023-05-30 江西科技学院 Axle welding area detection method and system based on computer vision
CN116309552B (en) * 2023-05-12 2023-08-29 西南交通大学 Method, device, equipment and medium for evaluating health state of existing line old retaining wall
CN116309552A (en) * 2023-05-12 2023-06-23 西南交通大学 Method, device, equipment and medium for evaluating health state of existing line old retaining wall
CN116309570B (en) * 2023-05-18 2023-08-04 山东亮马新材料科技有限公司 Titanium alloy bar quality detection method and system
CN116309570A (en) * 2023-05-18 2023-06-23 山东亮马新材料科技有限公司 Titanium alloy bar quality detection method and system
CN116363126A (en) * 2023-05-30 2023-06-30 东莞市立时电子有限公司 Welding quality detection method for data line USB plug
CN116363126B (en) * 2023-05-30 2023-08-22 东莞市立时电子有限公司 Welding quality detection method for data line USB plug
CN116363136B (en) * 2023-06-01 2023-08-11 山东创元智能设备制造有限责任公司 On-line screening method and system for automatic production of motor vehicle parts
CN116363136A (en) * 2023-06-01 2023-06-30 山东创元智能设备制造有限责任公司 On-line screening method and system for automatic production of motor vehicle parts
CN116385439B (en) * 2023-06-05 2023-08-15 山东兰通机电有限公司 Motor rubber shock pad quality detection method based on image processing
CN116385439A (en) * 2023-06-05 2023-07-04 山东兰通机电有限公司 Motor rubber shock pad quality detection method based on image processing
CN116416246A (en) * 2023-06-08 2023-07-11 临沂中科芯华新材料科技有限公司 Machine vision-based fully-degradable plastic product film coating effect evaluation method
CN116416246B (en) * 2023-06-08 2023-08-11 临沂中科芯华新材料科技有限公司 Machine vision-based fully-degradable plastic product film coating effect evaluation method
CN116433700B (en) * 2023-06-13 2023-08-18 山东金润源法兰机械有限公司 Visual positioning method for flange part contour
CN116433700A (en) * 2023-06-13 2023-07-14 山东金润源法兰机械有限公司 Visual positioning method for flange part contour
CN116433669A (en) * 2023-06-14 2023-07-14 山东兴华钢结构有限公司 Machine vision-based quality detection method for weld joints of steel frame of anti-seismic structure
CN116433669B (en) * 2023-06-14 2023-08-18 山东兴华钢结构有限公司 Machine vision-based quality detection method for weld joints of steel frame of anti-seismic structure
CN116485832A (en) * 2023-06-25 2023-07-25 山东九思新材料科技有限责任公司 Method for accurately detecting edges of non-uniform fluid impurities for recycling waste silicon wafers
CN116485832B (en) * 2023-06-25 2023-08-29 山东九思新材料科技有限责任公司 Method for accurately detecting edges of non-uniform fluid impurities for recycling waste silicon wafers
CN116542972A (en) * 2023-07-04 2023-08-04 山东阁林板建材科技有限公司 Wall plate surface defect rapid detection method based on artificial intelligence
CN116542972B (en) * 2023-07-04 2023-10-03 山东阁林板建材科技有限公司 Wall plate surface defect rapid detection method based on artificial intelligence
CN116630308B8 (en) * 2023-07-20 2023-10-27 山东华太新能源电池有限公司 Data enhancement system for battery welding anomaly detection
CN116630308B (en) * 2023-07-20 2023-09-26 山东华太新能源电池有限公司 Data enhancement system for battery welding anomaly detection
CN116630308A (en) * 2023-07-20 2023-08-22 山东华太新能源电池有限公司 Data enhancement system for battery welding anomaly detection
CN116612126B (en) * 2023-07-21 2023-09-19 青岛国际旅行卫生保健中心(青岛海关口岸门诊部) Container disease vector biological detection early warning method based on artificial intelligence
CN116612126A (en) * 2023-07-21 2023-08-18 青岛国际旅行卫生保健中心(青岛海关口岸门诊部) Container disease vector biological detection early warning method based on artificial intelligence
CN116664554A (en) * 2023-07-26 2023-08-29 微山晟轩机械制造有限公司 Bolt thread defect detection method based on image processing
CN116664554B (en) * 2023-07-26 2023-10-20 微山晟轩机械制造有限公司 Bolt thread defect detection method based on image processing
CN116681752A (en) * 2023-08-03 2023-09-01 山东墨氪智能科技有限公司 Method and device for calculating void ratio of void defects of DBC solder layer
CN116681752B (en) * 2023-08-03 2023-10-27 山东墨氪智能科技有限公司 Method and device for calculating void ratio of void defects of DBC solder layer
CN117197534A (en) * 2023-08-04 2023-12-08 广州电缆厂有限公司 Automatic detection method for cable surface defects based on feature recognition
CN117197534B (en) * 2023-08-04 2024-04-05 广州电缆厂有限公司 Automatic detection method for cable surface defects based on feature recognition
CN116703912A (en) * 2023-08-07 2023-09-05 深圳市鑫赛科科技发展有限公司 Mini-host network port integrity visual detection method
CN116703912B (en) * 2023-08-07 2023-11-24 深圳市鑫赛科科技发展有限公司 Mini-host network port integrity visual detection method
CN116703251B (en) * 2023-08-08 2023-11-17 德润杰(山东)纺织科技有限公司 Rubber ring production quality detection method based on artificial intelligence
CN116703251A (en) * 2023-08-08 2023-09-05 德润杰(山东)纺织科技有限公司 Rubber ring production quality detection method based on artificial intelligence
CN116777918A (en) * 2023-08-25 2023-09-19 苏州科尔珀恩机械科技有限公司 Visual auxiliary kiln surface defect rapid detection method
CN116777918B (en) * 2023-08-25 2023-10-31 苏州科尔珀恩机械科技有限公司 Visual auxiliary kiln surface defect rapid detection method
CN116823835A (en) * 2023-08-30 2023-09-29 山东省永星食品饮料有限公司 Bottled water impurity detection method based on machine vision
CN116823835B (en) * 2023-08-30 2023-11-10 山东省永星食品饮料有限公司 Bottled water impurity detection method based on machine vision
CN117102375B (en) * 2023-10-18 2024-01-02 沈阳欧施盾新材料科技有限公司 Special-shaped piece closing-in control method and equipment based on temperature imaging
CN117102375A (en) * 2023-10-18 2023-11-24 沈阳欧施盾新材料科技有限公司 Special-shaped piece closing-in control method and equipment based on temperature imaging
CN117197141A (en) * 2023-11-07 2023-12-08 山东远盾网络技术股份有限公司 Method for detecting surface defects of automobile parts
CN117197141B (en) * 2023-11-07 2024-01-26 山东远盾网络技术股份有限公司 Method for detecting surface defects of automobile parts
CN117274247A (en) * 2023-11-20 2023-12-22 深圳市海里表面技术处理有限公司 Visual detection method for quality of LTCC conductor surface coating
CN117274249A (en) * 2023-11-20 2023-12-22 江西省中鼐科技服务有限公司 Ceramic tile appearance detection method and system based on artificial intelligent image technology
CN117274247B (en) * 2023-11-20 2024-03-29 深圳市海里表面技术处理有限公司 Visual detection method for quality of LTCC conductor surface coating
CN117274249B (en) * 2023-11-20 2024-03-01 江西省中鼐科技服务有限公司 Ceramic tile appearance detection method and system based on artificial intelligent image technology
CN117314949A (en) * 2023-11-28 2023-12-29 山东远硕上池健康科技有限公司 Personnel injury detection and identification method based on artificial intelligence
CN117314949B (en) * 2023-11-28 2024-02-20 山东远硕上池健康科技有限公司 Personnel injury detection and identification method based on artificial intelligence
CN117372435B (en) * 2023-12-08 2024-02-06 智联信通科技股份有限公司 Connector pin detection method based on image characteristics
CN117372435A (en) * 2023-12-08 2024-01-09 智联信通科技股份有限公司 Connector pin detection method based on image characteristics
CN117455920B (en) * 2023-12-26 2024-03-22 武汉博源新材料科技集团股份有限公司 Artificial intelligence-based milk tea cup inferior product screening method and system
CN117455920A (en) * 2023-12-26 2024-01-26 武汉博源新材料科技集团股份有限公司 Artificial intelligence-based milk tea cup inferior product screening method and system
CN117576104A (en) * 2024-01-17 2024-02-20 山东世纪阳光科技有限公司 Visual detection method for health state of ultrafiltration membrane in purification process
CN117877008B (en) * 2024-03-13 2024-05-17 湖北神龙工程测试技术有限公司 Door and window performance detection method based on artificial intelligence

Also Published As

Publication number Publication date
CN110930353A (en) 2020-03-27
CN110930353B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
WO2021000524A1 (en) Hole protection cap detection method and apparatus, computer device and storage medium
CN110110799B (en) Cell sorting method, cell sorting device, computer equipment and storage medium
US10319096B2 (en) Automated tattoo recognition techniques
US10410292B2 (en) Method, system, apparatus, and storage medium for realizing antifraud in insurance claim based on consistency of multiple images
WO2020173177A1 (en) Object color difference defect detection method, device, computer device, and storage medium
US11593656B2 (en) Using a first stain to train a model to predict the region stained by a second stain
CN110210448B (en) Intelligent face skin aging degree identification and evaluation method
KR20200087297A (en) Defect inspection method and apparatus using image segmentation based on artificial neural network
KR102141302B1 (en) Object detection method based 0n deep learning regression model and image processing apparatus
CN110596120A (en) Glass boundary defect detection method, device, terminal and storage medium
CN111369523B (en) Method, system, equipment and medium for detecting cell stack in microscopic image
CN110307903B (en) Method for dynamically measuring non-contact temperature of specific part of poultry
CN111507426A (en) No-reference image quality grading evaluation method and device based on visual fusion characteristics
CN111126393A (en) Vehicle appearance refitting judgment method and device, computer equipment and storage medium
CN114066857A (en) Infrared image quality evaluation method and device, electronic equipment and readable storage medium
CN112419270A (en) No-reference image quality evaluation method and device under meta learning and computer equipment
CN111259971A (en) Vehicle information detection method and device, computer equipment and readable storage medium
CN106682604B (en) Blurred image detection method based on deep learning
CN112507869A (en) Underwater target behavior observation and water body environment monitoring method based on machine vision
CN110751623A (en) Joint feature-based defect detection method, device, equipment and storage medium
CN115909151A (en) Method for identifying serial number of motion container under complex working condition
CN115239947A (en) Wheat stripe rust severity evaluation method and device based on unsupervised learning
CN111523605B (en) Image identification method and device, electronic equipment and medium
CN113989632A (en) Bridge detection method and device for remote sensing image, electronic equipment and storage medium
WO2020107196A1 (en) Photographing quality evaluation method and apparatus for photographing apparatus, and terminal device

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19936156

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19936156

Country of ref document: EP

Kind code of ref document: A1