CN112184601A - Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm - Google Patents

Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm Download PDF

Info

Publication number
CN112184601A
CN112184601A CN202010942408.6A CN202010942408A CN112184601A CN 112184601 A CN112184601 A CN 112184601A CN 202010942408 A CN202010942408 A CN 202010942408A CN 112184601 A CN112184601 A CN 112184601A
Authority
CN
China
Prior art keywords
image
histogram
infrared light
near infrared
vein
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
CN202010942408.6A
Other languages
Chinese (zh)
Other versions
CN112184601B (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.)
China Jiliang University
Original Assignee
China Jiliang University
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 China Jiliang University filed Critical China Jiliang University
Priority to CN202010942408.6A priority Critical patent/CN112184601B/en
Publication of CN112184601A publication Critical patent/CN112184601A/en
Application granted granted Critical
Publication of CN112184601B publication Critical patent/CN112184601B/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
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques
    • 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
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows
    • 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/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

The invention provides a method for enhancing vein images under a near infrared light source by utilizing an improved CLAHE algorithm, which comprises the following steps: step 1: partitioning; step 2: calculating a histogram; and step 3: calculating the clipping threshold clip Limit: and 4, step 4: redistributing the pixel points; obtaining a histogram h' (x) after redistribution processing; and 5: histogram equalization; step 6: reconstructing a pixel point pre-output gray value; and 7: and reconstructing the gray value output by the pixel points. The image can be better in a strong dimming area or a bright area after the CLAHE image contrast enhancement algorithm is improved, and the characteristics of the vein area can be better emphasized in the vein image under a near infrared light source.

Description

Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm
Technical Field
The invention relates to a method for enhancing a vein image under a near infrared light source by utilizing an improved CLAHE algorithm, in particular to a method for adjusting gray value remapping distribution by introducing a sine function power exponential function.
Background
At present, most of researches on vein images are to highlight the structural features of veins by irradiating vein regions with near-infrared light sources and then adopting a filtering and contrast enhancing method.
Histogram equalization is the most common method of image contrast enhancement. Histogram equalization is mainly classified into a global method and a local method. The global enhancement method mainly achieves the purpose of contrast enhancement by modifying the distribution of the image histogram; the local enhancement method is mainly to define a local contrast in advance, and then enhance the local contrast to achieve the effect of enhancing image details, and the classical local histogram equalization technology includes Contrast Limited Adaptive Histogram Equalization (CLAHE), sub-block overlapping histogram equalization, and sub-block partial overlapping histogram equalization. The CLAHE method combines the advantages of adaptive histogram equalization and contrast limitation, is particularly suitable for low-contrast images, and is not complex in implementation process.
The general global histogram equalization algorithm has the problem that the histogram equalization result is not very balanced, the gray level of an equalized image is reduced, and certain details disappear; some images are equalized, and the contrast ratio is unnaturally excessively enhanced, so that the enhancement effect is not ideal. The local method comprehensively considers the position and the gray information of the pixel point, and the processing effect is often better than that of the global method. Classical local histogram equalization techniques include Contrast Limited Adaptive Histogram Equalization (CLAHE), subblock overlap histogram equalization, subblock partial overlap histogram equalization. The CLAHE method combines the advantages of adaptive histogram equalization and contrast limitation, is particularly suitable for low-contrast images, and is not complex in implementation process. However, when the CLAHE method is used for processing a vein image under a near-infrared light source, a strong dark area, namely a vein, cannot be better processed, so that vein features cannot be well highlighted when the image gray value is uniform.
Accordingly, there is a need for improvements in the art.
Disclosure of Invention
The invention aims to provide a method for enhancing vein images under a near infrared light source by using an improved CLAHE algorithm in an efficient manner.
In order to solve the technical problem, the invention provides a method for enhancing a vein image under a near infrared light source by using an improved CLAHE algorithm, which comprises the following steps:
step 1: partitioning; dividing an input image into c x r non-overlapping sub-blocks with equal size, wherein the number of pixels contained in each sub-block is M;
step 2: calculating a histogram; h (x) is used for representing the histogram of the sub-blocks, x represents the gray level, the value range is [0, L-1], and L is the possible gray level number;
and step 3: calculating the clipping threshold clip Limit:
Figure BDA0002674071310000021
in the formula: norm Clip Limit is the contrast enhancement amplitude value;
and 4, step 4: redistributing the pixel points; obtaining a histogram h' (x) after redistribution processing;
and 5: histogram equalization; performing histogram equalization processing on the histogram h' (x) subjected to redistribution processing, wherein an equalization result is represented by f (x);
step 6: reconstructing a pixel point pre-output gray value; obtaining the gray value of each pixel point of each sub-block according to the balance result f (x), and calculating the gray value g of each point in the pre-output image as a reference point;
and 7: and reconstructing the gray value output by the pixel points.
As an improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
in step 4:
for each sub-block, clipping the histogram h (x) of the sub-block by using the corresponding clipping threshold clip Limit, and uniformly redistributing the total difference value of the clipped pixels into each gray level of the histogram, wherein the clipping threshold clip Limit comprises the following steps:
Figure BDA0002674071310000022
avgBincr=totalE/L
wherein, total E refers to the total difference value of pixels exceeding the clip threshold Limit; avg BIncr refers to the pixel value in the histogram that averages each gray level increase;
and repeating the distribution process until all the cut pixel points are distributed.
As a further improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
in step 4:
the histogram of h (x) after redistribution process is represented by h' (x), and then
Figure BDA0002674071310000023
As a further improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
in step 6:
the gray value g of each point in the pre-output image is calculated by using a bilinear interpolation technology.
As a further improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
in step 7:
according to the gray value g of each point in the pre-output image obtained in the step 6, the gray value g is remapped through a sine function power exponential function to obtain the gray output value f of each pixel point in the final image, the effect that the vein image further enhances and highlights the vein area is achieved, and the method comprises the following steps:
Figure BDA0002674071310000031
in the formula, alpha is a compression degree coefficient of a dark area, and the area mapping degree of the dark area is controlled during the image gray value remapping; beta is a compression degree coefficient of the bright area, and the area mapping degree of the bright area is controlled during the remapping of the image gray value; and (3) emphasizing the bright/dark area characteristics of the whole image by utilizing the adjustment of alpha and beta parameters through the image characteristics of the sine function power exponent function.
As a further improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
the contrast enhancement amplitude value norm Clip Limit is 4 to 15.
As a further improvement to the method of the present invention for enhancing vein images under near infrared light sources using the improved CLAHE algorithm:
the compression degree coefficient α of the dark region is 1.8, and the compression degree coefficient β of the bright region is 1.9.
The technical scheme of the invention is as follows: a method for adjusting gray value remapping distribution by introducing a sine function power exponential function when gray value remapping is carried out on the basis of CLAHE image enhancement; the method comprises the following steps:
the method comprises the following steps: and obtaining a gray level mapping image of each sub-block region counted by the CLAHE.
Step two: and setting a compression degree coefficient alpha of the dark area, and controlling the area mapping degree of the dark area when the image gray value is remapped.
Step three: and setting a compression degree coefficient beta of the bright area, and controlling the area mapping degree of the bright area when the image gray value is remapped.
Step four: introducing sine function power exponential function
Figure BDA0002674071310000032
And (4) remapping the gray value after image enhancement, wherein g is the gray value which is pre-output after the original CLAHE algorithm is processed.
For example, as shown in fig. 3, the improved vein image result pair is that a in fig. 3 is an original image, b in fig. 3 is a processing result after classical histogram equalization, c in fig. 3 is a processing result with interpolated CLAHE, and d in fig. 3 is a final processing result of CLAHE improved by the present invention.
The method for enhancing the vein image under the near infrared light source by improving the CLAHE algorithm has the technical advantages that:
the image can be better in the strong dim area or the bright area after the CLAHE image contrast enhancement algorithm is improved, and the characteristics of the vein area can be better emphasized in the vein image under the near infrared light source.
Drawings
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
Figure 1 is the CLAHE algorithm. (a) Shearing a histogram; (b) a histogram after CLAHE redistribution;
FIG. 2 is a flow chart of image contrast enhancement by the modified CLAHE algorithm;
FIG. 3 is a modified vein image; a is original image, b is processed result after classical histogram equalization, c is processed result with interpolation CLAHE, d is improved CLAHE final processing result.
Detailed Description
The invention will be further described with reference to specific examples, but the scope of the invention is not limited thereto.
Example 1 method for enhancing vein images under near infrared light source using the modified CLAHE algorithm, as shown in figures 1-3,
the vein image vein area characteristic under the near-infrared light source is more obvious by performing the image enhancement method capable of emphasizing the controllable bright and dark areas by adjusting the image gray value remapping through CLAHE equalization.
CLAHE limits the enhancement amplitude of local contrast by limiting the height of the local histogram, thereby limiting the amplification of noise and excessive enhancement of local contrast. Firstly, dividing an image into a plurality of sub-blocks, then carrying out histogram 'shearing' on each sub-block (see (a) in figure 1), then carrying out histogram equalization on each sub-block (see (b) in figure 1), and finally carrying out interpolation operation on each pixel to obtain a transformed gray value, thereby realizing contrast-limited self-adaptive image enhancement.
As shown in fig. 2, a flowchart of an improved CLAHE algorithm mainly includes the following 7 steps:
step 1: and (5) partitioning. Dividing the input image into c x r non-overlapping sub-blocks with equal size, wherein each sub-block contains M pixels. The larger the sub-block, the more pronounced the enhancement effect, but the more image detail is lost.
Step 2: a histogram is calculated. The histogram of the sub-block is represented by h (x), x represents the gray level, and the value range is [0, L-1], and L is the possible number of 256 gray levels.
And step 3: calculating the clipping threshold clip Limit:
Figure BDA0002674071310000041
in the formula: norm Clip Limit is a contrast enhancement amplitude value that determines the amplitude of contrast enhancement, and is experimentally determined to be generally between 4 and 15.
And 4, step 4: and (6) redistributing the pixel points. For each sub-block, clipping the histogram h (x) of the sub-block by using the corresponding clipping threshold clip Limit, and uniformly redistributing the total difference value of the clipped pixels into each gray level of the histogram, wherein the clipping threshold clip Limit comprises the following steps:
Figure BDA0002674071310000051
avgBincr=totalE/L
in the formula, total E refers to the total difference of pixels exceeding the clip threshold Limit. avg BIncr refers to the pixel value in the histogram that increases on average per gray level.
The above allocation process is repeated until all the cut pixel points are allocated, as shown in (b) of fig. 1. If h' (x) represents the histogram of h (x) after the redistribution process, there are
Figure BDA0002674071310000052
And 5: and (5) histogram equalization. Histogram equalization is performed on the histogram h' (x) after the reallocation processing, and the equalization result is represented by f (x).
Step 6: and (5) reconstructing the gray value pre-output by the pixel points. And (4) obtaining the gray values of the central pixel points of the sub-blocks according to the balance result f (x), using the gray values as reference points, and calculating the gray value g of each point in the pre-output image by adopting a bilinear interpolation technology.
And 7: and reconstructing the gray value output by the pixel points. According to the gray value g of each point in the pre-output image obtained in the step 6, the gray value g is remapped through a sine function power exponential function to obtain the gray output value f of each pixel point in the final image, the effect that the vein image further enhances and highlights the vein area is achieved, and the method comprises the following steps:
Figure BDA0002674071310000053
where α is a compression degree coefficient of the dark region, and controls the region mapping degree of the dark region when the image gray-scale value is remapped. Beta is the compression degree coefficient of the bright area, and the area mapping degree of the bright area is controlled when the image gray value is remapped. And (3) emphasizing the bright/dark area characteristics of the whole image by utilizing the adjustment of alpha and beta parameters through the image characteristics of the sine function power exponent function. The vein region characteristics are more highlighted by the experiment that alpha is 1.8 and beta is 1.9 under the irradiation of a near infrared environment light source.
Finally, it is also noted that the above-mentioned lists merely illustrate a few specific embodiments of the invention. It is obvious that the invention is not limited to the above embodiments, but that many variations are possible. All modifications which can be derived or suggested by a person skilled in the art from the disclosure of the present invention are to be considered within the scope of the invention.

Claims (7)

1. The method for enhancing the vein image under the near infrared light source by utilizing the improved CLAHE algorithm is characterized by comprising the following steps of: the method comprises the following steps:
step 1: partitioning; dividing an input image into c x r non-overlapping sub-blocks with equal size, wherein the number of pixels contained in each sub-block is M;
step 2: calculating a histogram; h (x) is used for representing the histogram of the sub-blocks, x represents the gray level, the value range is [0, L-1], and L is the possible gray level number;
and step 3: calculating the clipping threshold clip Limit:
Figure FDA0002674071300000011
in the formula: norm Clip Limit is the contrast enhancement amplitude value;
and 4, step 4: redistributing the pixel points; obtaining a histogram h' (x) after redistribution processing;
and 5: histogram equalization; performing histogram equalization processing on the histogram h' (x) subjected to redistribution processing, wherein an equalization result is represented by f (x);
step 6: reconstructing a pixel point pre-output gray value; obtaining the gray value of each pixel point of each sub-block according to the balance result f (x), and calculating the gray value g of each point in the pre-output image as a reference point;
and 7: and reconstructing the gray value output by the pixel points.
2. A method of enhancing a vein image under a near infrared light source using a modified CLAHE algorithm as claimed in claim 1, wherein:
in step 4:
for each sub-block, clipping the histogram h (x) of the sub-block by using the corresponding clipping threshold clip Limit, and uniformly redistributing the total difference value of the clipped pixels into each gray level of the histogram, wherein the clipping threshold clip Limit comprises the following steps:
Figure FDA0002674071300000012
avgBincr=totalE/L
wherein, total E refers to the total difference value of pixels exceeding the clip threshold Limit; avg BIncr refers to the pixel value in the histogram that averages each gray level increase;
and repeating the distribution process until all the cut pixel points are distributed.
3. A method of enhancing a vein image under a near infrared light source using a modified CLAHE algorithm as claimed in claim 2, wherein:
in step 4:
the histogram of h (x) after redistribution process is represented by h' (x), and then
Figure FDA0002674071300000021
4. A method of enhancing a vein image under a near infrared light source using a modified CLAHE algorithm as claimed in claim 3, wherein:
in step 6:
the gray value g of each point in the pre-output image is calculated by using a bilinear interpolation technology.
5. Method for enhancing vein images under near infrared light sources with the improved CLAHE algorithm according to claim 4, characterized in that:
in step 7:
according to the gray value g of each point in the pre-output image obtained in the step 6, the gray value g is remapped through a sine function power exponential function to obtain the gray output value f of each pixel point in the final image, the effect that the vein image further enhances and highlights the vein area is achieved, and the method comprises the following steps:
Figure FDA0002674071300000022
in the formula, alpha is a compression degree coefficient of a dark area, and the area mapping degree of the dark area is controlled during the image gray value remapping; beta is a compression degree coefficient of the bright area, and the area mapping degree of the bright area is controlled during the remapping of the image gray value; and (3) emphasizing the bright/dark area characteristics of the whole image by utilizing the adjustment of alpha and beta parameters through the image characteristics of the sine function power exponent function.
6. Method for enhancing vein images under near infrared light sources with the improved CLAHE algorithm according to claim 5, characterized in that:
the contrast enhancement amplitude value norm Clip Limit is 4 to 15.
7. Method for enhancing vein images under near infrared light sources with the improved CLAHE algorithm according to claim 6, characterized in that:
the compression degree coefficient α of the dark region is 1.8, and the compression degree coefficient β of the bright region is 1.9.
CN202010942408.6A 2020-09-09 2020-09-09 Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm Active CN112184601B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010942408.6A CN112184601B (en) 2020-09-09 2020-09-09 Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010942408.6A CN112184601B (en) 2020-09-09 2020-09-09 Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm

Publications (2)

Publication Number Publication Date
CN112184601A true CN112184601A (en) 2021-01-05
CN112184601B CN112184601B (en) 2023-10-27

Family

ID=73920291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010942408.6A Active CN112184601B (en) 2020-09-09 2020-09-09 Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm

Country Status (1)

Country Link
CN (1) CN112184601B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114092361A (en) * 2022-01-20 2022-02-25 浙江芯昇电子技术有限公司 Image local contrast enhancement method for optimizing CLAHE
CN114283076A (en) * 2021-12-08 2022-04-05 凌云光技术股份有限公司 Image processing method based on transmissivity and CLAHE algorithm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106504177A (en) * 2016-09-30 2017-03-15 广东技术师范学院 A kind of low embedded rate compressed sensing general steganalysis method of coloured image
CN106651818A (en) * 2016-11-07 2017-05-10 湖南源信光电科技有限公司 Improved Histogram equalization low-illumination image enhancement algorithm
CN108389163A (en) * 2018-01-22 2018-08-10 西安工程大学 A kind of self-adapting enhancement method based on half-light coloured image
CN110852955A (en) * 2018-08-21 2020-02-28 中南大学 Image enhancement method based on image intensity threshold and adaptive cutting

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106504177A (en) * 2016-09-30 2017-03-15 广东技术师范学院 A kind of low embedded rate compressed sensing general steganalysis method of coloured image
CN106651818A (en) * 2016-11-07 2017-05-10 湖南源信光电科技有限公司 Improved Histogram equalization low-illumination image enhancement algorithm
CN108389163A (en) * 2018-01-22 2018-08-10 西安工程大学 A kind of self-adapting enhancement method based on half-light coloured image
CN110852955A (en) * 2018-08-21 2020-02-28 中南大学 Image enhancement method based on image intensity threshold and adaptive cutting

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
孔壮 等: "高分辨率视频图像实时增强***的设计", 视频应用与工程, vol. 39, no. 20 *
黎蔚 等: "手掌静脉图像增强算法研究", 微电子学与计算机, vol. 27, no. 7 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114283076A (en) * 2021-12-08 2022-04-05 凌云光技术股份有限公司 Image processing method based on transmissivity and CLAHE algorithm
CN114283076B (en) * 2021-12-08 2024-03-29 凌云光技术股份有限公司 Image processing method based on transmissivity and CLAHE algorithm
CN114092361A (en) * 2022-01-20 2022-02-25 浙江芯昇电子技术有限公司 Image local contrast enhancement method for optimizing CLAHE

Also Published As

Publication number Publication date
CN112184601B (en) 2023-10-27

Similar Documents

Publication Publication Date Title
CN107527333B (en) Quick image enhancement method based on gamma transformation
Stark Adaptive image contrast enhancement using generalizations of histogram equalization
CN112365424B (en) Local self-adaptive CLAHE-based infrared image denoising enhancement method, device and system and computer-readable storage medium
Arici et al. A histogram modification framework and its application for image contrast enhancement
EP0856813B1 (en) Method for improving contrast in picture sequences
RU2298226C1 (en) Method for improving digital images
EP1209621A2 (en) A method for enhancing a digital image based upon pixel color
CN112184601A (en) Method for enhancing vein image under near infrared light source by utilizing improved CLAHE algorithm
CN114092361B (en) Image local contrast enhancement method for optimizing CLAHE
Hassan1&2 et al. A new approach for contrast enhancement using sigmoid function
CN108961172B (en) Gamma correction-based image contrast enhancement method
JP2001216512A (en) Method for improving digital picture having noise- dependent control of texture
CN103985097B (en) A kind of image enchancing method of Nonuniform illumination
RU2342701C1 (en) Method of complexing digital multispectral half-tone images
CN117314801B (en) Fuzzy image optimization enhancement method based on artificial intelligence
JP4497160B2 (en) Image processing apparatus and method, and program
Abdullah-Al-Wadud A modified histogram equalization for contrast enhancement preserving the small parts in images
EP2229768A1 (en) A histogram-modeling based algorithm for contrast enhancement
CN114998165A (en) Histogram transformation based uniform exposure image contrast enhancement method
Zhu et al. Image enhancement using fuzzy intensity measure and adaptive clipping histogram equalization
WO2020107308A1 (en) Low-light-level image rapid enhancement method and apparatus based on retinex
Murahira et al. Modified histogram equalization for image contrast enhancement
Omarova et al. Application of the Clahe method contrast enhancement of X-Ray Images
Ohchi et al. A nonlinear filter system for beautifying facial images with contrast enhancement
CN116012273A (en) Image enhancement method and device based on local gray scale fluctuation rate

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