CN115631173A - Composite film defect identification method - Google Patents

Composite film defect identification method Download PDF

Info

Publication number
CN115631173A
CN115631173A CN202211336527.2A CN202211336527A CN115631173A CN 115631173 A CN115631173 A CN 115631173A CN 202211336527 A CN202211336527 A CN 202211336527A CN 115631173 A CN115631173 A CN 115631173A
Authority
CN
China
Prior art keywords
gray
composite film
index
image
pixel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211336527.2A
Other languages
Chinese (zh)
Other versions
CN115631173B (en
Inventor
郑小平
苏文晓
陈奎
张继林
齐鹏堂
简粤
胡玫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lanzhou University of Technology
Original Assignee
Lanzhou University of Technology
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 Lanzhou University of Technology filed Critical Lanzhou University of Technology
Priority to CN202211336527.2A priority Critical patent/CN115631173B/en
Publication of CN115631173A publication Critical patent/CN115631173A/en
Application granted granted Critical
Publication of CN115631173B publication Critical patent/CN115631173B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • 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
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/187Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Quality & Reliability (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Multimedia (AREA)
  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to the technical field of defect identification, in particular to a composite film defect identification method, which comprises the following steps: acquiring a surface gray level image of the composite film, and determining the composite film image suspected of having defects according to the gray level value; segmenting the image to obtain a suspected defect area; obtaining a gray level deviation index of a pixel point according to the gray level value difference between a central pixel point in a window with a set size and other pixel points in the window; dividing the pixel points into three grades according to the gray deviation index, determining possible defect pixel points, and performing region growth by using the possible defect pixel points to obtain possible defect regions; the method comprises the steps of obtaining end points of a line segment where the principal component direction of the possible defect area is located, obtaining a non-connectivity index according to the distance between the end points corresponding to the possible defect area, determining a connected area according to the non-connectivity index and an index threshold value, and identifying the defect of the composite film according to the connected area.

Description

Composite film defect identification method
Technical Field
The invention relates to the technical field of defect identification, in particular to a composite film defect identification method.
Background
The composite film is a high polymer material compounded by two or more layers of films made of different materials and is mainly used for packaging. The surface wrinkle of the composite film is a common problem in the processing and application processes of the composite flexible packaging material, the surface wrinkle problem is represented by that the film on the surface layer is convex upwards, the surface wrinkle problem can influence the aesthetic property and the quality of the composite film package, and therefore, the surface defect identification needs to be carried out on the composite film processed at each production stage to identify the wrinkle defect.
The conventional identification method is to process the surface image of the composite film by using a threshold segmentation method, but because the wrinkle defect on the surface of the composite film is mainly in an approximate tunnel shape, the gray difference is not obvious, most of the composite films are white, the threshold is influenced by the reflection condition in the presence of illumination, the image segmentation result is inaccurate, and the composite film defect identification result is also inaccurate.
Disclosure of Invention
In order to solve the above technical problems, an object of the present invention is to provide a method for identifying defects of a composite film, which adopts the following technical scheme:
acquiring a surface gray image of the composite film, and determining the composite film image suspected of having defects according to the gray value of a pixel point in the surface gray image; performing threshold segmentation on the composite film image to obtain a suspected defect area;
performing sliding window processing on the pixels in the suspected defect area by using a window with a set size, and obtaining a gray level deviation index of the pixels according to the gray level value of the central pixel in the window and the gray level value difference of other pixels in the window; dividing the pixel points into three grades according to the gray deviation index, determining possible defect pixel points according to the number of different grades in a window corresponding to the pixel points, and performing region growth by using the possible defect pixel points to obtain a possible defect region;
the method comprises the steps of obtaining end points of a line section where the principal component direction of a possible defect area is located, obtaining a non-connectivity index according to the distance between the end points corresponding to the possible defect area, determining a connected area according to the non-connectivity index and an index threshold value, and identifying the defect of the composite film according to the connected area.
Preferably, the determining the composite film image suspected of having the defect according to the gray value of the pixel point in the surface gray image specifically includes:
and calculating a color second moment of the image according to the gray value of the pixel point in the surface gray image, obtaining a color characteristic value corresponding to the surface gray image according to the color second moment, and marking the surface gray image with the color characteristic value larger than a color threshold value as a composite film image suspected of having defects.
Preferably, the obtaining of the suspected defect area by performing threshold segmentation on the composite film image specifically includes:
constructing a gray level histogram according to gray levels of pixel points in the composite film image, acquiring gray levels corresponding to two wave troughs in the gray level histogram, recording a smaller gray level corresponding to the two wave troughs as a first segmentation threshold, and recording a larger gray level corresponding to the two wave troughs as a second segmentation threshold; and the area formed by the pixel points with the gray value smaller than the first segmentation threshold and the area formed by the pixel points with the gray value larger than the second segmentation threshold are suspected defect areas.
Preferably, the method for obtaining the gray level deviation index specifically includes:
Figure DEST_PATH_IMAGE002
wherein,
Figure DEST_PATH_IMAGE004
the gray scale deviation index corresponding to the z-th pixel point is shown,
Figure DEST_PATH_IMAGE005
represents the gray value of the z-th pixel point,
Figure 100002_DEST_PATH_IMAGE007
representing the gray value of the first other pixel in the window where the z-th pixel is located,
Figure DEST_PATH_IMAGE009
representing the gray value of the jth pixel point in the window where the jth pixel point is located,
Figure DEST_PATH_IMAGE010
representing the total number of other pixels except the center point within the window, exp () represents an exponential function with a natural constant e as the base.
Preferably, the dividing the pixel points into three levels according to the gray scale deviation index specifically includes: and acquiring a difference value between the gray level deviation indexes of any two pixel points, and dividing all the pixel points into three levels based on all the difference values.
Preferably, the determining the possible defective pixel points according to the number of different levels in the window corresponding to the pixel points specifically includes:
for any pixel point in a window, if three levels of pixel points exist in all pixel points in the window, the any pixel point is a possible defect pixel point; if the pixel points of the three grades do not exist, the window is enlarged according to the set step length, if the pixel points of the three grades exist in all the pixel points in the enlarged window, any one pixel point is a possible defect area, and if the pixel points of the three grades do not exist in all the pixel points in the enlarged window, any one pixel point is a noise pixel point.
Preferably, the determining the connected region according to the non-connectivity index and the index threshold specifically includes:
when the non-connectivity indexes corresponding to the end points of the possible defect area are all smaller than the index threshold value, connecting the possible defect area with the possible defect area corresponding to the non-connectivity index smaller than the index threshold value; when the non-connectivity index corresponding to one end point of the possible defect area is smaller than the index threshold, calculating the shortest distance from the other end point to the edge of the composite film image, if the shortest distance is smaller than the index threshold, connecting the possible defect area with the possible defect area corresponding to the non-connectivity index smaller than the index threshold, and if the shortest distance is greater than or equal to the index threshold, discarding the possible defect area without connection; when the non-connectivity indexes corresponding to the end points of the possible defect areas are not less than the index threshold, the possible defect areas are discarded and are not connected; all connected regions constitute a connected region.
The embodiment of the invention at least has the following beneficial effects:
determining a composite film image suspected of having defects according to gray values of pixel points in the image, roughly analyzing the image by combining color distribution characteristics in the image, screening out the image with the defects, and carrying out subsequent further judgment; then, the composite film image is processed by using the window, the gray scale deviation index of the pixel point is calculated, the possible defect pixel point is obtained according to the gray scale deviation index, and the gray scale difference condition of the pixel point in the neighborhood is considered by analyzing the gray scale difference of the pixel point in the window range, so that more accurate defect information can be obtained; and finally, connecting all the defect areas according to the connectivity among the defect areas to further obtain a defect identification result, considering the shape characteristics of the defects of the composite film, namely the connectivity to a certain degree, further increasing the accuracy of the obtained defect areas and enabling the defect identification result to be more accurate.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions and advantages of the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for identifying defects of a composite film according to the present invention.
Detailed Description
To further illustrate the technical means and effects of the present invention adopted to achieve the predetermined objects, the following detailed description of the method for identifying defects of a composite film according to the present invention, its specific implementation, structure, features and effects will be given in conjunction with the accompanying drawings and the preferred embodiments. In the following description, different "one embodiment" or "another embodiment" refers to not necessarily the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The following describes a specific scheme of the composite film defect identification method provided by the invention in detail with reference to the accompanying drawings.
Example (b):
the main purposes of the invention are: and (3) carrying out surface defect detection on the composite film after the composite film is processed at each production stage by an image processing technology, wherein the defects are mainly wrinkle defects.
The specific scenes aimed by the invention are as follows: and (4) performing wrinkle defect detection after the compound film is processed by a machine, after curing, when the bag making processing is finished and after a period of time of placement, when the water boiling processing is finished and after a period of time of placement, when the boiling processing is finished and after a period of time of placement.
Referring to fig. 1, a flowchart of a method for identifying defects of a composite film according to an embodiment of the present invention is shown, where the method includes the following steps:
acquiring a surface gray image of a composite film, and determining the composite film image suspected of having defects according to the gray value of a pixel point in the surface gray image; and performing threshold segmentation on the composite film image to obtain a suspected defect area.
Firstly, after the composite film is processed by a film making machine, after curing treatment, when bag making processing is finished and after a period of time of placement, when boiling treatment is finished and after a period of time of placement, image acquisition is carried out at each stage, and the surface image of the composite film is obtained. Because the acquisition process is influenced by mechanical noise, the acquired surface image of the composite film is subjected to noise reduction treatment. In this embodiment, a gaussian filtering method is used to perform noise reduction on the surface image, then a semantic segmentation method is used to remove the interference of the background, and finally the image obtained after the semantic segmentation is subjected to graying processing to obtain a surface grayscale image of the composite film.
It should be noted that, in order to reduce the amount of calculation, image analysis needs to be performed on the images acquired at each stage, that is, whether a wrinkle defect exists in the surface gray-scale image of the composite film corresponding to the current stage is determined in advance according to the gray-scale features in the images, so that only the surface gray-scale image with the wrinkle defect is subjected to subsequent analysis.
When the composite film has a wrinkle defect, a convex tunnel-shaped wrinkle appears on the surface of the composite film, and the area where the wrinkle exists is shown in an image, the color of the convex part in the wrinkle part is bright and white, but the color of the concave parts on two sides of the convex part is dark and gray. From this, it was found that if wrinkle defects exist on the surface of the composite film, bright white and dark gray appear in the surface tone image of the composite film due to the generation of wrinkle defects.
Based on this, carry out analysis according to the grey scale characteristic of pixel in the surface gray scale image of composite film, judge the color distribution range in the surface gray scale image of composite film, when the composite film does not have the fold defect, the grey scale value in the surface gray scale image of composite film is all comparatively close, shows as a comparatively close colour value, and when the composite film has the fold defect, can increase the colour that two kinds of grey scale values are different in the surface gray scale image of composite film, shows as three colour values in the image.
Then, since the color moment is a simple and effective color feature labeling method, and the color second moment can represent the distribution range of the color in the image, in this embodiment, the color second moment of the image is calculated based on the gray value of the pixel point in the surface gray image of the composite film
Figure DEST_PATH_IMAGE011
Is formulated as:
Figure DEST_PATH_IMAGE013
wherein,
Figure DEST_PATH_IMAGE014
and expressing the gray value of the ith pixel point in the surface gray image of the composite film, and N expressing the total number of the pixel points in the surface gray image. And calculating the difference value between the gray value of each pixel point and the image gray mean value in the formula, reflecting the deviation degree of each pixel point relative to the image gray mean value, summing the deviation degrees of all the pixel points, and averaging to express the size of the color distribution range of the surface gray image. When in use
Figure DEST_PATH_IMAGE011A
The larger the image color distribution range, the wider the image color distribution range, the more various color values may exist in the surface gray level image, otherwise
Figure DEST_PATH_IMAGE011AA
The smaller the size, i.e. the narrower the color distribution range of the representative image, the more the surface gray scale map is illustratedThere may be fewer color values in the image, i.e., there may be one color value.
In this embodiment, the color second moment obtained based on the gray value of the pixel point in the surface gray image of the composite film
Figure DEST_PATH_IMAGE011AAA
Performing normalization processing, namely obtaining the color characteristic value of the surface gray level image according to the color second moment, and expressing the color characteristic value as
Figure DEST_PATH_IMAGE016
Wherein
Figure DEST_PATH_IMAGE017
in order to obtain a color characteristic value after normalization processing is carried out on the color second moment,
Figure DEST_PATH_IMAGE011AAAA
is the second moment of the color of the image, and e is a natural constant.
When in use
Figure DEST_PATH_IMAGE011_5A
The larger the value of (a) is,
Figure DEST_PATH_IMAGE017A
the larger the value of (A) is, the wider the distribution range of the expressed image color is, which indicates that more color values may exist in the surface gray level image, the wrinkle defect may exist on the surface of the composite film corresponding to the image, otherwise, when the color value of (A) is larger, the wrinkle defect may exist on the surface of the composite film corresponding to the image
Figure DEST_PATH_IMAGE011_6A
The smaller the size of the product is,
Figure DEST_PATH_IMAGE017AA
the smaller the value of (b), that is, the narrower the distribution range of the representative image color, indicates that there may be fewer color values in the surface gray image, that is, there may be one color value, and there may be no wrinkle defect on the surface of the composite film corresponding to the image.
Setting color threshold, in this embodimentThe value of (3) is 0.4, when the value of the color characteristic value corresponding to the surface gray level image is greater than the color threshold value, that is, the value is
Figure DEST_PATH_IMAGE019
In the case of the composite film image, the color distribution range in the surface gray level image is large, and the wrinkle defect may exist, so that the surface gray level image is marked as the composite film image suspected of having the defect. When the value of the color characteristic value corresponding to the surface gray level image is less than or equal to the color threshold value, that is to say
Figure DEST_PATH_IMAGE021
In the present embodiment, the surface gray image of the composite film having the color feature value less than or equal to the color threshold is not subjected to the subsequent analysis because the color distribution range in the surface gray image is narrow and the wrinkle defect may not exist.
Finally, in the area where the wrinkles exist on the composite film image, the raised parts on the wrinkles are bright white due to light reflection, and the bottoms of the mountain roots at the two sides of the raised parts on the wrinkles are shaded dark gray due to the shielding of the raised parts. All the areas where wrinkles exist referred to in this embodiment include a raised bright white portion and a dark gray-shaded portion at the bottom of mountain roots on both sides of the raised portion.
Based on the above, a gray level histogram is constructed according to gray levels of pixel points in the composite film image suspected of having the defect, and the composite film image is segmented to obtain a suspected defect area by taking the gray levels corresponding to two troughs in the gray level histogram as segmentation thresholds.
Specifically, the gray values corresponding to two troughs in the gray histogram are used as segmentation thresholds, the gray value corresponding to the two troughs is smaller and is used as a first segmentation threshold, the gray value corresponding to the two troughs is larger and is used as a second segmentation threshold, the image is segmented into three regions by using the first segmentation threshold and the second segmentation threshold, the first segmentation threshold is used as a segmentation point between a dark gray shadow part and a dominant color part of the composite film, so that a region formed by pixel points of which the gray values are smaller than the first segmentation threshold is used as a wrinkle convex region, and the second segmentation threshold is used as a segmentation point between the dominant color part of the composite film and a bright white convex part, so that a region formed by pixel points of which the gray values are larger than the second segmentation threshold is used as a wrinkle convex region. And the area formed by the pixel points of which the gray value is greater than or equal to the first segmentation threshold and less than or equal to the second segmentation threshold is a normal area, and no wrinkle defect exists in the normal area. Wherein, the fold raised area and the fold shadow area form a suspected defect area.
Performing sliding window processing on the pixels in the suspected defect area by using a window with a set size, and obtaining a gray level deviation index of the pixels according to the gray level value of the central pixel in the window and the gray level value difference of other pixels in the window; dividing the pixel points into three grades according to the gray deviation index, determining possible defect pixel points according to the quantity of different grades in a window corresponding to the pixel points, and performing region growth by using the possible defect pixel points to obtain a possible defect region.
Firstly, it should be noted that the area where the wrinkle exists includes a convex portion on the wrinkle and a root bottom shadow portion on two sides of the convex portion on the wrinkle, the composite film image is processed by using a window with a set size, sliding is started from left to right and from top to bottom from the upper left corner of the composite film image, the sliding step length is 1, when the composite film image slides to a position, the area where the pixel belongs to is judged according to the gray value of the pixel at the center position of the window, and when the pixel at the center position belongs to a possible defect area, the pixel distribution characteristics in the window are analyzed.
When the pixel point at the central position belongs to the possible defect area, namely the pixel point belongs to the wrinkle convex area or the wrinkle shadow area, the wrinkle defect exists in the window, the convex part and the shadow part in the wrinkle exist at the same time and are distributed in a centralized manner, and then in the window with a certain size, the pixel point belongs to the wrinkle convex area and the wrinkle shadow area, which is shown in color distribution, the gray value of the pixel point in the window can exist in two intervals, namely, the gray value is smaller than the first segmentation threshold value and larger than the second segmentation threshold value, namely, the bright white pixel point and the dark gray pixel point exist. After the size of the window is enlarged, the window contains a normal area of the composite film except for the defect part, three color distributions necessarily exist in the window, namely the gray value of the pixel points in the window exists in three intervals. Meanwhile, the pixel points of the three color distributions in the window are not distributed in disorder, but are distributed in three regions in a concentrated manner.
In this embodiment, the value of the set size is 3*3, which can be set by an implementer according to the actual situation.
Then, a window 3*3 is used for processing a possible defect area in the composite film image, for any pixel point, in a window 3*3 taking the pixel point as a center, if a wrinkle area exists in the window, gray value differences of the pixel point at the center position in the window and other pixel points are respectively calculated, and the pixel deviation condition of the pixel point at the center position can be obtained according to the difference between each difference.
Obtaining the gray scale deviation index of the pixel point according to the gray scale value of the central pixel point in the window and the gray scale value difference of other pixel points in the window, and expressing the gray scale deviation index by a formula as follows:
Figure DEST_PATH_IMAGE023
wherein,
Figure DEST_PATH_IMAGE025
indicating the gray scale deviation index corresponding to the z-th pixel point,
Figure DEST_PATH_IMAGE026
the gray value of the z-th pixel point is represented,
Figure DEST_PATH_IMAGE028
representing the gray value of the first other pixel in the window where the z-th pixel is located,
Figure DEST_PATH_IMAGE030
representing the gray value of the jth pixel point in the window where the jth pixel point is located,exp () represents an exponential function with a natural constant e as the base,
Figure DEST_PATH_IMAGE031
the total number of other pixels except the center point in the representation window is 8 in this embodiment.
Figure DEST_PATH_IMAGE033
And calculating the deviation of the difference between the gray value of the z-th pixel point and the gray values of other j-th pixel points in the window by taking the difference between the gray value of the z-th pixel point and the gray value of other first pixel points in the window as a reference value, and further calculating the deviation between the reference value and other differences in sequence. The difference value between the gray values of the pixels in the window can reflect the gray difference of the pixels, and the deviation degree between the gray differences is calculated by selecting any one gray difference as a reference value.
The smaller the gray value difference between the pixel points is, the smaller the corresponding gray value difference is, which indicates that there may be no wrinkle defect in the window where the pixel point is located, the smaller the deviation degree between the gray value differences is, the smaller the value of the gray value deviation index of the pixel point is, which indicates that the pixel point is more likely to belong to a normal pixel point on the composite film image. The larger the gray value difference between the pixel points is, the larger the corresponding gray value difference is, which indicates that a window where the pixel point is located may have a wrinkle defect, and the larger the deviation degree between the gray value differences is, the larger the value of the gray value deviation index of the pixel point is, which indicates that the pixel point is more likely to belong to the pixel point in the possible defect region on the composite film image.
And finally, respectively calculating the gray level deviation indexes of all pixel points in the possible defect area on the composite film image according to the method, acquiring the difference value between the gray level deviation indexes of any two pixel points, and dividing all the pixel points into three grades based on all the difference values.
Specifically, any two pixel points are selected to be marked as a first pixel point and a second pixel point, the difference value of the gray scale deviation indexes corresponding to the two pixel points is calculated, if the value of the difference value is greater than a threshold value, the two pixel points are considered to be more unlikely to belong to the pixel point in the same region, and therefore the gray scale deviation indexes corresponding to the two pixel points are divided into two different levels; if the value of the difference is smaller than or equal to the threshold, the two pixel points are considered to possibly belong to the pixel points in the same region, so that the gray level deviation indexes corresponding to the two pixel points are divided into the same grade.
And then, one other pixel point is selected arbitrarily to be marked as a third pixel point, the difference value of the gray scale deviation index corresponding to the third pixel point and the gray scale deviation index corresponding to the first pixel point is calculated respectively, and the difference value of the gray scale deviation index corresponding to the third pixel point and the gray scale deviation index corresponding to the second pixel point is calculated. If the difference value is smaller than or equal to the threshold value, namely the difference value of the gray scale deviation index corresponding to the third pixel point and the first pixel point is smaller than or equal to the threshold value, dividing the gray scale deviation index corresponding to the third pixel point into the grade of the gray scale deviation index corresponding to the first pixel point, or dividing the gray scale deviation index corresponding to the third pixel point into the grade of the gray scale deviation index corresponding to the second pixel point, wherein the difference value of the gray scale deviation index corresponding to the third pixel point and the gray scale deviation index corresponding to the second pixel point is smaller than or equal to the threshold value. And if the difference is not smaller than or equal to the threshold, namely the difference between the gray scale deviation indexes of the third pixel point and the other two pixel points is larger than the threshold, dividing the gray scale deviation index corresponding to the third pixel point into a third grade. In this embodiment, the value of the threshold is 0.2, and an implementer can set the threshold according to actual conditions.
Based on the above, by analogy, all the pixel points are divided into three levels based on all the difference values, the graded standard is that the level with the minimum gray scale deviation index is divided into a first level, the level with the larger gray scale deviation is divided into a second level, and the level with the maximum gray scale deviation is divided into a third level. If the composite film image has a wrinkle defect, the color distribution of the pixel points in the image is definitely represented by three colors, and the gray scale deviation index corresponding to the pixel points in the image is definitely divided into three levels. Meanwhile, each pixel point corresponds to one gray level deviation index, each gray level deviation index corresponds to one grade, each pixel point corresponds to one grade, and the three grades can be approximately regarded as three different areas, namely a normal area, a fold convex area and a fold shadow area.
It should be noted that, in this embodiment, by considering the deviation of the gray level difference between a pixel point and a pixel point in its neighborhood, the change condition of the gray level difference in the image can be analyzed in more detail, and the pixel point is determined based on the change condition of the difference, so that more accurate defect information can be obtained.
Further, the pixel point in the window where the pixel point is located is judged, in the window with a certain size, the pixel point belongs to a fold convex area and a fold shadow area, and the result shows that the deviation of the gray difference exists, the corresponding grade of the pixel point in the window possibly has two different grades, after the size of the window is enlarged, the window also comprises a normal area of a composite film except for a defect part, three color distributions must exist in the window, and the corresponding grade of the pixel point in the enlarged window has three different grades.
Therefore, for any pixel point in the window, if three levels of pixel points exist in all pixel points in the window, it is indicated that three pixel points in different areas possibly exist in the window where the pixel point is located, and the any pixel point is a possible defect pixel point; if the pixel points with the three grades do not exist, the window is enlarged according to the set step length, if the pixel points with the three grades exist in all the pixel points in the enlarged window, any one pixel point is a possible defect area, if the pixel points with the three grades do not exist in all the pixel points in the enlarged window, the fact that the pixel point which is judged at present does not belong to a defect part and is only a noise point is possible is indicated, and any one pixel point is a noise pixel point. In this embodiment, the value of the set step length is 2, that is, the size of the enlarged window is 5*5, and an implementer can set the value of the set step length according to an actual situation.
After the possible defect pixel points are obtained according to the method, the possible defect pixel points are respectively used as initial seed points for area growth, the area obtained by the area growth is marked as a possible defect area, wherein a rule implementer for the area growth can carry out setting according to actual conditions, for example, when the difference value between the gray values of the pixel points is smaller than a set threshold value, the growth is carried out, and the value of the set threshold value needs to be set by the implementer according to the actual conditions.
And step three, acquiring end points of a line segment where the principal component direction of the possible defect area is located, obtaining a non-connectivity index according to the distance between the end points corresponding to the possible defect area, determining a connected area according to the non-connectivity index and an index threshold, and identifying the defect of the composite film according to the connected area.
It should be noted that the wrinkle defect existing on the surface of the composite film may also be referred to as a tunnel defect because the raised portions on the wrinkles are all around and continuously converge or diverge, so as to form complete wrinkles similar to a tunnel shape, that is, the wrinkles are all interconnected and exist on the composite film, and there is no discrete distribution of the wrinkle portions. Therefore, for any possible defect region, connectivity between the possible defect region and other possible defect regions is analyzed respectively, and the confidence level that the possible defect region is a real defect region is judged according to the connectivity.
Specifically, PCA principal component analysis is performed on each possible defect region to obtain a principal component direction of each possible defect region, and two end points of a line segment where the principal component direction of the possible defect region is located are obtained, and then the two end points can be regarded as two end points of a wrinkle line segment in the possible defect region.
Because the wrinkle defect is in a tunnel shape in the image, the regions with the wrinkle defect are communicated with each other, and on the basis, whether the possible defect region is communicated with other possible defect regions is judged by judging the distance between two end points corresponding to the possible defect region and the distance between end points corresponding to other possible defect regions.
Based on this, a non-connectivity index is obtained according to the distance between the endpoints corresponding to the possible defect region, the non-connectivity indexes between the two endpoints corresponding to the possible defect region and the endpoints of the other possible defect regions are respectively recorded as a first non-connectivity index and a second non-connectivity index, that is, the two endpoints corresponding to the possible defect region respectively correspond to one non-connectivity index, which is expressed by a formula:
Figure DEST_PATH_IMAGE035
wherein,
Figure DEST_PATH_IMAGE037
representing a first non-connectivity indicator corresponding to a first end point in the possible defect region c,
Figure DEST_PATH_IMAGE039
represents a second non-connectivity index corresponding to a second end point in the possible defect region c, and min () represents a function for finding a minimum value.
And establishing a rectangular coordinate system taking pixel points as a unit by taking the upper left corner of the image as an origin, wherein the abscissa and the ordinate of the pixel points are the number of columns and the number of rows in the image respectively.
Figure DEST_PATH_IMAGE040
And
Figure DEST_PATH_IMAGE041
respectively showing the abscissa and ordinate of the first end point in the possible defect area c on the image of the composite film,
Figure DEST_PATH_IMAGE042
and
Figure DEST_PATH_IMAGE043
respectively representing the abscissa and ordinate of the first end point in the possible defect area c on the image,
Figure DEST_PATH_IMAGE044
and
Figure DEST_PATH_IMAGE045
respectively representing the abscissa and ordinate of the k-th end point of the other possible defect area on the image,
Figure DEST_PATH_IMAGE046
and
Figure DEST_PATH_IMAGE047
respectively, the abscissa and ordinate of the o-th end point of the other possible defect area on the image.
Figure DEST_PATH_IMAGE049
And
Figure DEST_PATH_IMAGE051
each represents the minimum of the distances between the end points in the current possible defect region and the end points in the other possible defect regions. When the minimum distance between the end points corresponding to different possible defect regions is smaller than the index threshold, it indicates that certain connectivity exists between the different possible defect regions, and the different possible defect regions need to be connected. That is, the smaller the minimum value of the distance between the endpoints, the greater the connectivity between the possible defect regions, the smaller the non-connectivity index corresponding to the possible defect region.
Meanwhile, considering that there may be a wrinkle portion located at an edge, the shortest distance between two end points corresponding to each possible defect region and the edge of the composite film image needs to be calculated and respectively recorded as a first non-edge index and a second non-edge index, and the smaller the shortest distance is, the closer the end points are to the edge of the image, the smaller the value of the corresponding non-edge index is.
Further, a connected region is determined according to the non-connectivity index and the index threshold, when the non-connectivity indexes corresponding to the end points of the possible defect region are both smaller than the index threshold, that is, the first non-connectivity index and the second non-connectivity index corresponding to the two end points of the region are both smaller than the index threshold, the possible defect region is connected with the possible defect region corresponding to the non-connectivity index smaller than the index threshold, that is, the region is connected with the possible defect region where the end point corresponding to the first non-connectivity index is located, the region is connected with the possible defect region where the end point corresponding to the second non-connectivity index is located, and the connected region is the connected region.
When the non-connectivity index corresponding to only one end point of the possible defect area is smaller than the threshold, calculating the shortest distance from the other end point to the edge of the composite film image, wherein the shortest distance is the non-marginal index corresponding to the end point, if the non-marginal index is smaller than the index threshold, connecting the possible defect area with the possible defect area corresponding to the non-connectivity index smaller than the threshold, and the connected area is a connected area.
For example, in two endpoints of a possible defect region, a first non-connectivity index corresponding to a first endpoint is smaller than an index threshold, and a second non-connectivity index corresponding to a second endpoint is larger than the index threshold, the second non-marginal index corresponding to the second endpoint is compared with the index threshold, if the second non-connectivity index is larger than or equal to the index threshold, the possible defect region is possibly a discrete region, and therefore the possible defect region needs to be discarded and not connected, and if the second non-connectivity index is smaller than the index threshold, one side of the possible defect region is connected with other possible defect regions, and the other side of the possible defect region is close to an image edge, and therefore the possible defect region needs to be connected with the possible defect region where the endpoint corresponding to the first non-connectivity index is located.
When the non-connectivity index corresponding to the end point of the possible defect area is not smaller than the index threshold, that is, the first non-connectivity index and the second non-connectivity index corresponding to the two end points of the area are not smaller than the index threshold, it is indicated that connectivity does not exist between the possible defect area and other possible defect areas, the possible defect area is a discrete area, the possible defect area is discarded, and connection is not performed.
In this embodiment, the value of the index threshold is 4, and an implementer can set the index threshold according to an actual situation.
And judging all possible defect areas according to the method, connecting the possible defect areas needing to be connected, forming a connected area by all the connected areas, and segmenting the composite film image according to the connected area so as to identify the defects of the composite film. In this embodiment, the pixel values of the pixels in the connected region are set to 255, and the pixel values of the other pixels are set to 0, so as to segment the composite film image, and perform defect identification by using the segmented image to obtain an identification result. The method for segmenting the image and the method for identifying the defect are various, and an implementer can select the method according to the actual situation.
It should be noted that, the invention analyzes the possibility that the pixel point is the pixel point of the wrinkle region first, then grows the possible wrinkle connected domain, analyzes the possibility that the possible wrinkle connected domain forms the tunnel, thereby identifies the complete wrinkle region, and overcomes the difficulty that the threshold segmentation cannot identify the wrinkle with fine gray difference.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications or substitutions do not cause the essential features of the corresponding technical solutions to depart from the scope of the technical solutions of the embodiments of the present application, and are intended to be included within the scope of the present application.

Claims (7)

1. A composite film defect identification method is characterized by comprising the following steps:
acquiring a surface gray image of the composite film, and determining the composite film image suspected of having defects according to the gray value of a pixel point in the surface gray image; performing threshold segmentation on the composite film image to obtain a suspected defect area;
performing sliding window processing on the pixels in the suspected defect area by using a window with a set size, and obtaining a gray level deviation index of the pixels according to the gray level value of the central pixel in the window and the gray level value difference of other pixels in the window; dividing the pixel points into three grades according to the gray level deviation index, determining possible defect pixel points according to the number of different grades in a window corresponding to the pixel points, and performing region growth by using the possible defect pixel points to obtain a possible defect region;
the method comprises the steps of obtaining end points of a line section where the principal component direction of a possible defect area is located, obtaining a non-connectivity index according to the distance between the end points corresponding to the possible defect area, determining a connected area according to the non-connectivity index and an index threshold value, and identifying the defect of the composite film according to the connected area.
2. The method according to claim 1, wherein the determining the composite film image suspected of having the defect according to the gray-scale values of the pixels in the surface gray-scale image specifically comprises:
and calculating a color second moment of the image according to the gray value of the pixel point in the surface gray image, obtaining a color characteristic value corresponding to the surface gray image according to the color second moment, and marking the surface gray image with the color characteristic value larger than a color threshold value as a composite film image suspected of having defects.
3. The method according to claim 1, wherein the step of performing threshold segmentation on the composite film image to obtain the suspected defect area specifically comprises:
constructing a gray level histogram according to gray levels of pixel points in the composite film image, acquiring gray levels corresponding to two wave troughs in the gray level histogram, recording a smaller gray level corresponding to the two wave troughs as a first segmentation threshold, and recording a larger gray level corresponding to the two wave troughs as a second segmentation threshold; and the area formed by the pixel points with the gray value smaller than the first segmentation threshold and the area formed by the pixel points with the gray value larger than the second segmentation threshold are suspected defect areas.
4. The method for identifying the defects of the composite film according to claim 1, wherein the method for obtaining the gray scale deviation index specifically comprises the following steps:
Figure 906552DEST_PATH_IMAGE002
wherein,
Figure DEST_PATH_IMAGE003
is shown as
Figure 89271DEST_PATH_IMAGE004
The gray scale deviation index corresponding to each pixel point,
Figure 400167DEST_PATH_IMAGE005
the gray value of the z-th pixel point is represented,
Figure 326535DEST_PATH_IMAGE006
representing the gray value of the first other pixel in the window where the z-th pixel is located,
Figure DEST_PATH_IMAGE007
representing the gray value of the jth pixel point in the window where the jth pixel point is located,
Figure 970268DEST_PATH_IMAGE008
representing the total number of other pixels except the center point within the window, exp () represents an exponential function with a natural constant e as the base.
5. The method for identifying defects of a composite film according to claim 1, wherein the dividing of the pixels into three levels according to the gray scale deviation index specifically comprises: and acquiring a difference value between the gray level deviation indexes of any two pixel points, and dividing all the pixel points into three levels based on all the difference values.
6. The method for identifying the defects of the composite film according to claim 1, wherein the step of determining the possible defect pixel points according to the number of different grades in the window corresponding to the pixel points specifically comprises the following steps:
for a window where any one pixel point is located, if pixel points of three grades exist in all the pixel points in the window, the any one pixel point is a possible defect pixel point; if the pixel points of the three grades do not exist, the window is expanded according to the set step length, if the pixel points of the three grades exist in all the pixel points in the expanded window, any one pixel point is a possible defect area, and if the pixel points of the three grades do not exist in all the pixel points in the expanded window, any one pixel point is a noise pixel point.
7. The method for identifying the defects of the composite film as claimed in claim 1, wherein the determining the connected region according to the non-connectivity index and the index threshold specifically comprises:
when the non-connectivity indexes corresponding to the end points of the possible defect area are all smaller than the index threshold, connecting the possible defect area with the possible defect area corresponding to the non-connectivity index smaller than the index threshold; when the non-connectivity index corresponding to only one end point of the possible defect area is smaller than the index threshold, calculating the shortest distance from the other end point to the edge of the composite film image, if the shortest distance is smaller than the index threshold, connecting the possible defect area with the possible defect area corresponding to the non-connectivity index smaller than the index threshold, and if the shortest distance is greater than or equal to the index threshold, removing the possible defect area without connection; when the non-connectivity indexes corresponding to the end points of the possible defect areas are not less than the index threshold, the possible defect areas are discarded and are not connected; all connected regions constitute connected regions.
CN202211336527.2A 2022-10-28 2022-10-28 Composite film defect identification method Active CN115631173B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211336527.2A CN115631173B (en) 2022-10-28 2022-10-28 Composite film defect identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211336527.2A CN115631173B (en) 2022-10-28 2022-10-28 Composite film defect identification method

Publications (2)

Publication Number Publication Date
CN115631173A true CN115631173A (en) 2023-01-20
CN115631173B CN115631173B (en) 2023-04-18

Family

ID=84908719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211336527.2A Active CN115631173B (en) 2022-10-28 2022-10-28 Composite film defect identification method

Country Status (1)

Country Link
CN (1) CN115631173B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116091499A (en) * 2023-04-07 2023-05-09 山东中胜涂料有限公司 Abnormal paint production identification system
CN116228775A (en) * 2023-05-10 2023-06-06 实德电气集团有限公司 Contactor integrity detection method based on machine vision
CN116246174A (en) * 2023-04-26 2023-06-09 山东金诺种业有限公司 Sweet potato variety identification method based on image processing
CN116433623A (en) * 2023-03-31 2023-07-14 杭州数创自动化控制技术有限公司 Defect position marking and identifying method, system, equipment and medium
CN116858854A (en) * 2023-09-04 2023-10-10 季华实验室 Doping concentration correction method and device, electronic equipment and storage medium
CN116977335A (en) * 2023-09-22 2023-10-31 山东贞元汽车车轮有限公司 Intelligent detection method for pitting defects on surface of mechanical part

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000132687A (en) * 1998-10-28 2000-05-12 Nec Corp Defect inspection device and defect inspection method
JP2004177397A (en) * 2002-10-01 2004-06-24 Tokyo Seimitsu Co Ltd Method and apparatus for inspecting image defect, and device for inspecting appearance
JP2004233222A (en) * 2003-01-30 2004-08-19 Mitsubishi Heavy Ind Ltd Inspection device and inspection method of printed matter
JP2006292503A (en) * 2005-04-08 2006-10-26 Sumitomo Electric Ind Ltd Method and device for flaw inspection
CN103499585A (en) * 2013-10-22 2014-01-08 常州工学院 Non-continuity lithium battery thin film defect detection method and device based on machine vision
CN112651923A (en) * 2020-11-11 2021-04-13 北京平恒智能科技有限公司 Adhesive film wrinkle defect detection method capable of removing fine residues based on area ratio
CN115100171A (en) * 2022-07-11 2022-09-23 常宝云 Steel die welding defect detection method and system based on machine vision
CN115100208A (en) * 2022-08-26 2022-09-23 南通三信塑胶装备科技股份有限公司 Film surface defect evaluation method based on histogram and dynamic light source
CN115131354A (en) * 2022-08-31 2022-09-30 江苏森信达生物科技有限公司 Laboratory plastic film defect detection method based on optical means
CN115170576A (en) * 2022-09-09 2022-10-11 山东中发新材料科技有限公司 Aluminum pipe surface defect detection method based on machine vision
CN115170572A (en) * 2022-09-08 2022-10-11 山东瑞峰新材料科技有限公司 BOPP composite film surface gluing quality monitoring method

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000132687A (en) * 1998-10-28 2000-05-12 Nec Corp Defect inspection device and defect inspection method
JP2004177397A (en) * 2002-10-01 2004-06-24 Tokyo Seimitsu Co Ltd Method and apparatus for inspecting image defect, and device for inspecting appearance
JP2004233222A (en) * 2003-01-30 2004-08-19 Mitsubishi Heavy Ind Ltd Inspection device and inspection method of printed matter
JP2006292503A (en) * 2005-04-08 2006-10-26 Sumitomo Electric Ind Ltd Method and device for flaw inspection
CN103499585A (en) * 2013-10-22 2014-01-08 常州工学院 Non-continuity lithium battery thin film defect detection method and device based on machine vision
CN112651923A (en) * 2020-11-11 2021-04-13 北京平恒智能科技有限公司 Adhesive film wrinkle defect detection method capable of removing fine residues based on area ratio
CN115100171A (en) * 2022-07-11 2022-09-23 常宝云 Steel die welding defect detection method and system based on machine vision
CN115100208A (en) * 2022-08-26 2022-09-23 南通三信塑胶装备科技股份有限公司 Film surface defect evaluation method based on histogram and dynamic light source
CN115131354A (en) * 2022-08-31 2022-09-30 江苏森信达生物科技有限公司 Laboratory plastic film defect detection method based on optical means
CN115170572A (en) * 2022-09-08 2022-10-11 山东瑞峰新材料科技有限公司 BOPP composite film surface gluing quality monitoring method
CN115170576A (en) * 2022-09-09 2022-10-11 山东中发新材料科技有限公司 Aluminum pipe surface defect detection method based on machine vision

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHUNG-FENG JEFFREY KUO ET AL.: "Integrating image processing and classification technology into automated polarizing film defect inspection" *
SHAOHUA DONG ET AL.: "Automatic defect identification technology of digital image of pipeline weld" *
丰艳,王明辉,陈一民: "利用图像像素灰度值变化速度的相似性进行图像分割" *
李浩然等: "基于分块阈值LBP算法的光学薄膜表面缺陷分割" *
樊向党;林波;沈文和;: "塑料薄膜表面疵点检测及识别方法研究" *
苑玮琦等: "圆柱形覆膜锂电池圆周面破膜检测方法研究" *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116433623A (en) * 2023-03-31 2023-07-14 杭州数创自动化控制技术有限公司 Defect position marking and identifying method, system, equipment and medium
CN116091499A (en) * 2023-04-07 2023-05-09 山东中胜涂料有限公司 Abnormal paint production identification system
CN116091499B (en) * 2023-04-07 2023-06-20 山东中胜涂料有限公司 Abnormal paint production identification system
CN116246174A (en) * 2023-04-26 2023-06-09 山东金诺种业有限公司 Sweet potato variety identification method based on image processing
CN116246174B (en) * 2023-04-26 2023-08-08 山东金诺种业有限公司 Sweet potato variety identification method based on image processing
CN116228775A (en) * 2023-05-10 2023-06-06 实德电气集团有限公司 Contactor integrity detection method based on machine vision
CN116228775B (en) * 2023-05-10 2023-07-04 实德电气集团有限公司 Contactor integrity detection method based on machine vision
CN116858854A (en) * 2023-09-04 2023-10-10 季华实验室 Doping concentration correction method and device, electronic equipment and storage medium
CN116977335A (en) * 2023-09-22 2023-10-31 山东贞元汽车车轮有限公司 Intelligent detection method for pitting defects on surface of mechanical part
CN116977335B (en) * 2023-09-22 2023-12-12 山东贞元汽车车轮有限公司 Intelligent detection method for pitting defects on surface of mechanical part

Also Published As

Publication number Publication date
CN115631173B (en) 2023-04-18

Similar Documents

Publication Publication Date Title
CN115631173B (en) Composite film defect identification method
CN115018828B (en) Defect detection method for electronic component
US10565479B1 (en) Identifying and excluding blurred areas of images of stained tissue to improve cancer scoring
CN115578389B (en) Defect detection method of groove MOS device
CN114820627B (en) Textile quality detection method based on optical means
CN115829883A (en) Surface image denoising method for dissimilar metal structural member
CN106096491B (en) Automatic identification method for microaneurysms in fundus color photographic image
CN115578732B (en) Label identification method for fertilizer production line
CN117058147B (en) Environment-friendly plastic product defect detection method based on computer vision
CN114140445B (en) Breast cancer pathological image identification method based on key attention area extraction
CN115239718B (en) Plastic product defect detection method and system based on image processing
CN115546203A (en) Production monitoring and analyzing method based on image data algorithm
CN115511907B (en) Scratch detection method for LED screen
CN114037691A (en) Carbon fiber plate crack detection method based on image processing
CN109871900A (en) The recognition positioning method of apple under a kind of complex background based on image procossing
CN116805316B (en) Degradable plastic processing quality detection method based on image enhancement
CN115082451A (en) Stainless steel soup ladle defect detection method based on image processing
CN114299051A (en) Leather material surface defect detection method based on feature modeling significance detection
CN111738931B (en) Shadow removal algorithm for aerial image of photovoltaic array unmanned aerial vehicle
CN114155202A (en) Thyroid nodule ultrasonic image classification method based on feature fusion and transfer learning
CN115880280A (en) Detection method for quality of steel structure weld joint
US11967123B1 (en) Binarization method for CT sectional image of fiber package containing artifacts
CN107239761B (en) Fruit tree branch pulling effect evaluation method based on skeleton angular point detection
Deepa et al. Improved watershed segmentation for apple fruit grading
CN116843581A (en) Image enhancement method, system, device and storage medium for multi-scene graph

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant