CN105809703A - Adhesion hemocyte image segmentation method based on improved fractional differential and graph theory - Google Patents

Adhesion hemocyte image segmentation method based on improved fractional differential and graph theory Download PDF

Info

Publication number
CN105809703A
CN105809703A CN201610189166.1A CN201610189166A CN105809703A CN 105809703 A CN105809703 A CN 105809703A CN 201610189166 A CN201610189166 A CN 201610189166A CN 105809703 A CN105809703 A CN 105809703A
Authority
CN
China
Prior art keywords
region
improved
algorithm
weighted value
image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610189166.1A
Other languages
Chinese (zh)
Other versions
CN105809703B (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.)
Fuzhou University
Original Assignee
Fuzhou 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 Fuzhou University filed Critical Fuzhou University
Priority to CN201610189166.1A priority Critical patent/CN105809703B/en
Publication of CN105809703A publication Critical patent/CN105809703A/en
Application granted granted Critical
Publication of CN105809703B publication Critical patent/CN105809703B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20112Image segmentation details
    • G06T2207/20152Watershed segmentation
    • 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/20172Image enhancement details
    • G06T2207/20192Edge enhancement; Edge preservation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The invention relates to an adhesion hemocyte image segmentation method based on the improved fractional differential and graph theory.The method comprises the steps of pretreating a hemocyte image by combining morphological denoising with the improved quasi-circular mask operator fractional differential algorithm aiming at the phenomenon that the hemocyte image is fuzzy and low in contrast ratio, wherein by the adoption of the improved fractional differential algorithm, cell edge details are reserved while staining contamination and grain noise of the hemocyte image are filtered out; then conducting preliminary segmentation on the pretreated image with the watershed algorithm, and mapping an over-segmentation area into a node; finally conducting resegmentation on the cell image obtained in the second step with the improved graph theory minimum spanning tree (MST) algorithm.By the adoption of the method, segmentation precision of adherent cells in the cell image can be improved.

Description

Adhesion blood cell image dividing method based on the fractional order differential improved and graph theory
Technical field
The present invention relates to Medical Image Segmentation Techniques field, particularly a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory.
Background technology
Cell is the basic composition unit of all life body.Whether the biomedical important directions quickly grown in recent years is through the identification of cell, counting and texture image and is distorted and diagnoses the illness.And the data of cell image can be analyzed exactly, it is critically depend on and can split cell image exactly.
Dividing method based on graph theory is the focus studied both at home and abroad in the last few years, and a kind of layering dividing method from top to down based on drawing method that Sharon in 2006 etc. above propose at " Nature ", segmentation result is accurate and efficiency is high.Vanhamel etc. propose the non-linear multiple dimensioned Segmentation of Color Image based on graph theory;Hu Xuegang etc. propose based on graph theory and normalized segmentation criterion, and digital image processing model (LIP model) is applied in image procossing;Ye Wei etc. are theoretical in conjunction with Mumford Shah, it is proposed that a kind of optimization method, by considering the geometric properties in interregional combination degree and each region in image, based on the weights of conjugation the image partition method adding it to minimum spanning tree between zoning;Zhang et al. proposes watershed and the image partition method of graph theory combination.AnnaFabijanska et al. proposes the innovatory algorithm of a kind of minimum spanning tree based on graph theory, and this algorithm passes through the number of vertex in minimizing figure thus improving the speed of image segmentation.Multi-spectral Satellite Images is split by the Mean-Shift algorithm that BiplabBanerjee et al. proposes on the basis of the algorithm of a kind of minimum spanning tree based on graph theory in conjunction with clustering algorithm and improvement.Additionally, Felzenszwalb and Huttenlocher (be called for short FH) propose " little and and it " merging criterion, carrying out the improvement of minimum spanning tree partitioning algorithm, the method efficiency is higher, it is possible to partly carry out different segmentation according to different picture characteristics.But it also has the shortcoming of self, the k value preset in algorithm is difficult to ground manual control, if its value is more big, will produce merging phenomenon;If too small, it is impossible to effectively to suppress the generation of zonule, more redundant area will be produced.Just improve based on FH algorithm herein.
Summary of the invention
In view of this, the purpose of the present invention is to propose to a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory, improve the precision of AC segmentation in cell image.
The present invention adopts below scheme to realize: a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory, comprises the following steps:
Step S1: for the phenomenon that blood cell image is fuzzy and contrast is not high, being combined by the fractional order differential algorithm of the similar round mask operator of morphology denoising and improvement and blood cell image carries out pretreatment, the fractional order differential algorithm of described improvement pollutes in the dyeing filtering blood cell image and remains cell edges details preferably while grain noise;
Step S2: with watershed algorithm, pretreated image is carried out just segmentation, overdivided region is mapped as node;
Step S3: cell image step S2 obtained with graph theory minimum spanning tree (MST) algorithm improved is split again.
Further, in described step S1, it is proposed to a kind of Tiansi operator based on vector merge similar round improve mask template, specifically include following steps:
Step S11: the square Tiansi fractional order differential mask template of one 5x5 of structure, Tiansi mask operator is square structure, a on the position of 4 drift angles of outermost layer2Treating as is the starting point 4 vectors at template initial pointAnd take its half
Step S12: the upper left corner of the mask operator template obtained after described step S11 is processed is with the lower left cornerMerge into level new vector to the left
Step S13: the combined vector obtained after described step S12 is processedBe added to the vector in original third layer centre position, obtains the mask template of similar round.
Further, described step S2 specifically includes following steps:
Step S21: assume that image obtains 6 overdivided regions after watershed algorithm is split: 1,2,3,4,5,6}, the limit collection of its mapped node is:
E={e1,e2,e3,e4,e5,e6,e7,e8,e9}={ (1,2), (1,5), (2,3), (2,5), (3,4), (3,5), (4,5), (4,6), (5,6) };
The weighted value of described limit collection respectively 5,9,6,8,4,1,7,4,2;
Step S22: the weighted value on overdivided region and limit is mapped one by one and obtains initial weighted undirected graph.
Further, described step S3 specifically includes following steps:
Step S31: map: the overdivided region of input picture is mapped as node, calculates weighted undirected graph G={V, E, the ψ of 8 neighborhood systems by the weighted value definition improved }, wherein | V |=n, | E |=m, connect node viAnd vjThe weighted value on limit is ψ (vi,vj), setting area minimum area is p;
Step S32: state initialization: set Sq-1It is the cut set after q-1 sub-region merges, makes q=1 make initial segmentation be Sq-1=S0, S0=(v1,v2,…,vn), namely in V, each element is a region, by prime area internal diversityRegion areaWherein i=1,2,3 ... n;
Step S33: sequence: weighted undirected graph G={V, E, ψ } in all of weighted value according to ascending order arrangement, obtain queue π=(O1,O2,…,Om);
Step S34: if the queue in step S33 is not empty, makes the limit dequeue that weighted value is minimum, otherwise forward step S37 to;
Step S35: from Sq-1Structure Sq: setWithIt is S respectivelyq-1The two end node v on the limit of middle dequeueiAnd vjThe region at place, if:
And
Then by regionWithMerge and obtain Sq, and calculate the weighted value of new region and its neighborhood system;Otherwise nonjoinder, makes Sq=Sq-1
Step S36: the weighted value according to new region, updates the internal diversity of new region, region difference and weighted value queue π, makes q=q+1, returns step S34;
Step S37: traversal SqEach region, (vi,vj) ∈ E, vi∈ C1, vj∈ C2, if | C1 | < p or | C2 | < p, then remerges C1 and C2;
Step S38: be same color by the pixel assignment belonging to the same area, it is simple to the observation analysis of human eye;
Step S39: the segmentation result of the minimal spanning tree algorithm that output improves herein.
Compared with prior art, the invention has the beneficial effects as follows be effectively improved in cell segmentation over-segmentation, less divided precision, there is application prospect widely.
Accompanying drawing explanation
Fig. 1 is the inventive method schematic flow sheet.
Fig. 2 is embodiment of the present invention Tiansi mask operator schematic diagram.
Fig. 3 is the mask operator schematic diagram of embodiment of the present invention similar round.
Fig. 4 is embodiment of the present invention vector Unite principle schematic diagram.
Fig. 5 is embodiment of the present invention watershed algorithm overdivided region schematic diagram.
Fig. 6 is embodiment of the present invention weighted undirected graph.
Detailed description of the invention
Below in conjunction with drawings and Examples, the present invention will be further described.
As it is shown in figure 1, present embodiments provide a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory, comprise the following steps:
Step S1: for the phenomenon that blood cell image is fuzzy and contrast is not high, being combined by the fractional order differential algorithm of the similar round mask operator of morphology denoising and improvement and blood cell image carries out pretreatment, the fractional order differential algorithm of described improvement pollutes in the dyeing filtering blood cell image and remains cell edges details preferably while grain noise;
Step S2: with watershed algorithm, pretreated image is carried out just segmentation, overdivided region is mapped as node;
Step S3: cell image step S2 obtained with graph theory minimum spanning tree (MST) algorithm improved is split again.
In the present embodiment, in described step S1, select the template of suitable size for the size of cell in cell image, it is proposed to a kind of Tiansi operator based on vector merge similar round improve mask template, specifically include following steps:
Step S11: the square Tiansi fractional order differential mask template of one 5x5 of structure, as in figure 2 it is shown, Tiansi mask operator is square structure, therefore can a on the position of 4 drift angles of outermost layer2Treating as is the starting point 4 vectors at template initial pointAnd take its half
Step S12: the upper left corner of the mask operator template obtained after described step S11 is processed is with the lower left cornerAccording to the principle that Fig. 4 vector merges, merge into level new vector to the left
Step S13: the combined vector obtained after described step S12 is processedBe added to the vector in original third layer centre position, just obtains the mask template of similar round, as shown in Figure 3.
In the present embodiment, described step S2 specifically includes following steps:
Step S21: assume that image obtains 6 overdivided regions as shown in Figure 5 after watershed algorithm is split: 1,2,3,4,5,6}, the limit collection of its mapped node is:
E={e1,e2,e3,e4,e5,e6,e7,e8,e9}={ (1,2), (1,5), (2,3), (2,5), (3,4), (3,5), (4,5), (4,6), (5,6) };
The weighted value of described limit collection respectively 5,9,6,8,4,1,7,4,2;
Step S22: the weighted value on overdivided region and limit is mapped one by one and obtains initial weighted undirected graph, as shown in Figure 6.
In the present embodiment, described step S3 specifically includes following steps:
Step S31: map: the overdivided region of input picture is mapped as node, calculates weighted undirected graph G={V, E, the ψ of 8 neighborhood systems by the weighted value definition improved }, wherein | V |=n, | E |=m, connect node viAnd vjThe weighted value on limit is ψ (vi,vj), setting area minimum area is p;
Step S32: state initialization: set Sq-1It is the cut set after q-1 sub-region merges, makes q=1 make initial segmentation be Sq-1=S0, S0=(v1,v2,…,vn), namely in V, each element (node) is a region, by prime area internal diversityRegion areaWherein i=1,2,3 ... n;
Step S33: sequence: weighted undirected graph G={V, E, ψ } in all of weighted value according to ascending order arrangement, obtain queue π=(O1,O2,…,Om);
Step S34: if the queue in step S33 is not empty, makes the limit dequeue that weighted value is minimum, otherwise forward step S37 to;
Step S35: from Sq-1Structure Sq: setWithIt is S respectivelyq-1The two end node v on the limit of middle dequeueiAnd vjThe region at place, if:
And
Then by regionWithMerge and obtain Sq, and calculate the weighted value of new region and its neighborhood system;Otherwise nonjoinder, makes Sq=Sq-1
Step S36: the weighted value according to new region, updates the internal diversity of new region, region difference and weighted value queue π, makes q=q+1, returns step S34;
Step S37: traversal SqEach region, (vi,vj) ∈ E, vi∈ C1, vj∈ C2, if | C1 | < p or | C2 | < p, then remerges C1 and C2;
Step S38: be same color by the pixel assignment belonging to the same area, it is simple to the observation analysis of human eye;
Step S39: the segmentation result of the minimal spanning tree algorithm that output improves herein.
The foregoing is only presently preferred embodiments of the present invention, all equalizations done according to the present patent application the scope of the claims change and modify, and all should belong to the covering scope of the present invention.

Claims (4)

1. the adhesion blood cell image dividing method based on the fractional order differential improved and graph theory, it is characterised in that: comprise the following steps:
Step S1: for the phenomenon that blood cell image is fuzzy and contrast is not high, being combined by the fractional order differential algorithm of the similar round mask operator of morphology denoising and improvement and blood cell image carries out pretreatment, the fractional order differential algorithm of described improvement pollutes in the dyeing filtering blood cell image and remains cell edges details preferably while grain noise;
Step S2: with watershed algorithm, pretreated image is carried out just segmentation, overdivided region is mapped as node;
Step S3: cell image step S2 obtained with graph theory minimum spanning tree (MST) algorithm improved is split again.
2. a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory according to claim 1, it is characterized in that: in described step S1, propose a kind of Tiansi operator based on vector merge similar round improve mask template, specifically include following steps:
Step S11: the square Tiansi fractional order differential mask template of one 5x5 of structure, Tiansi mask operator is square structure, a on the position of 4 drift angles of outermost layer2Treating as is the starting point 4 vectors at template initial pointAnd take its half
Step S12: the upper left corner of the mask operator template obtained after described step S11 is processed is with the lower left cornerMerge into level new vector to the left
Step S13: the combined vector obtained after described step S12 is processedBe added to the vector in original third layer centre position, obtains the mask template of similar round.
3. a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory according to claim 1, it is characterised in that: described step S2 specifically includes following steps:
Step S21: assume that image obtains 6 overdivided regions after watershed algorithm is split: 1,2,3,4,5,6}, the limit collection of its mapped node is:
E={e1,e2,e3,e4,e5,e6,e7,e8,e9}={ (1,2), (1,5), (2,3), (2,5), (3,4), (3,5), (4,5), (4,6), (5,6) };
The weighted value of described limit collection respectively 5,9,6,8,4,1,7,4,2;
Step S22: the weighted value on overdivided region and limit is mapped one by one and obtains initial weighted undirected graph.
4. a kind of adhesion blood cell image dividing method based on the fractional order differential improved and graph theory according to claim 1, it is characterised in that: described step S3 specifically includes following steps:
Step S31: map: the overdivided region of input picture is mapped as node, calculates weighted undirected graph G={V, E, the ψ of 8 neighborhood systems by the weighted value definition improved }, wherein | V |=n, | E |=m, connect node viAnd vjThe weighted value on limit is ψ (vi,vj), setting area minimum area is p;
Step S32: state initialization: set Sq-1It is the cut set after q-1 sub-region merges, makes q=1 make initial segmentation be Sq-1=S0, S0=(v1,v2,…,vn), namely in V, each element is a region, by prime area internal diversityRegion areaWherein i=1,2,3 ... n;
Step S33: sequence: weighted undirected graph G={V, E, ψ } in all of weighted value according to ascending order arrangement, obtain queue π=(O1,O2,…,Om);
Step S34: if the queue in step S33 is not empty, makes the limit dequeue that weighted value is minimum, otherwise forward step S37 to;
Step S35: from Sq-1Structure Sq: setWithIt is S respectivelyq-1The two end node v on the limit of middle dequeueiAnd vjThe region at place, if:
And
Then by regionWithMerge and obtain Sq, and calculate the weighted value of new region and its neighborhood system;Otherwise nonjoinder, makes Sq=Sq-1
Step S36: the weighted value according to new region, updates the internal diversity of new region, region difference and weighted value queue π, makes q=q+1, returns step S34;
Step S37: traversal SqEach region, (vi,vj) ∈ E, vi∈ C1, vj∈ C2, if | C1 | < p or | C2 | < p, then remerges C1 and C2;
Step S38: be same color by the pixel assignment belonging to the same area, it is simple to the observation analysis of human eye;
Step S39: the segmentation result of the minimal spanning tree algorithm that output improves herein.
CN201610189166.1A 2016-03-30 2016-03-30 Adhesion blood cell image dividing method based on improved fractional order differential and graph theory Expired - Fee Related CN105809703B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610189166.1A CN105809703B (en) 2016-03-30 2016-03-30 Adhesion blood cell image dividing method based on improved fractional order differential and graph theory

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610189166.1A CN105809703B (en) 2016-03-30 2016-03-30 Adhesion blood cell image dividing method based on improved fractional order differential and graph theory

Publications (2)

Publication Number Publication Date
CN105809703A true CN105809703A (en) 2016-07-27
CN105809703B CN105809703B (en) 2018-10-30

Family

ID=56455047

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610189166.1A Expired - Fee Related CN105809703B (en) 2016-03-30 2016-03-30 Adhesion blood cell image dividing method based on improved fractional order differential and graph theory

Country Status (1)

Country Link
CN (1) CN105809703B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107977978A (en) * 2017-11-01 2018-05-01 浙江农林大学 A kind of method and system that cytoadherence mode is determined based on elliptical shaft relation
CN109583479A (en) * 2018-11-07 2019-04-05 北京市遥感信息研究所 One kind being based on morphologic high spectrum image feature extracting method
CN112614104A (en) * 2020-12-22 2021-04-06 湖南伊鸿健康科技有限公司 Segmentation counting method for red blood cell overlapping and terminal
CN113139973A (en) * 2021-04-01 2021-07-20 武汉市疾病预防控制中心 Artificial intelligence-based plasmodium identification method and equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006056786A1 (en) * 2004-11-26 2006-06-01 Snell & Wilcox Limited Image segmentation
CN101841642A (en) * 2010-04-22 2010-09-22 南京航空航天大学 Edge detection method based on fractional-order signal processing
CN102136145A (en) * 2011-04-27 2011-07-27 中国科学院遥感应用研究所 Region merging method for threshold-restrained minimum spanning tree algorithm
CN102800051A (en) * 2012-03-28 2012-11-28 中南大学 Self-adaptation flotation froth image enhancement method based on fractional order differential

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006056786A1 (en) * 2004-11-26 2006-06-01 Snell & Wilcox Limited Image segmentation
US20080037870A1 (en) * 2004-11-26 2008-02-14 Snell & Wilcox Limited Image Segmentation
CN101841642A (en) * 2010-04-22 2010-09-22 南京航空航天大学 Edge detection method based on fractional-order signal processing
CN102136145A (en) * 2011-04-27 2011-07-27 中国科学院遥感应用研究所 Region merging method for threshold-restrained minimum spanning tree algorithm
CN102800051A (en) * 2012-03-28 2012-11-28 中南大学 Self-adaptation flotation froth image enhancement method based on fractional order differential

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
冯林 等: "基于分水岭变换和图论的图像分割方法", 《仪器仪表学报》 *
卢志茂 等: "二次分水岭和Ncut相结合的彩色图像分割方法", 《华中科技大学学报(自然科学版)》 *
赵建: "分数阶微分在图像纹理增强中的应用", 《液晶与显示》 *
郭戈 等: "分水岭算法在重叠细胞图象分割中的应用", 《微计算机信息》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107977978A (en) * 2017-11-01 2018-05-01 浙江农林大学 A kind of method and system that cytoadherence mode is determined based on elliptical shaft relation
CN107977978B (en) * 2017-11-01 2020-08-21 浙江农林大学 Method and system for determining cell adhesion mode based on elliptical axis relation
CN109583479A (en) * 2018-11-07 2019-04-05 北京市遥感信息研究所 One kind being based on morphologic high spectrum image feature extracting method
CN112614104A (en) * 2020-12-22 2021-04-06 湖南伊鸿健康科技有限公司 Segmentation counting method for red blood cell overlapping and terminal
CN113139973A (en) * 2021-04-01 2021-07-20 武汉市疾病预防控制中心 Artificial intelligence-based plasmodium identification method and equipment

Also Published As

Publication number Publication date
CN105809703B (en) 2018-10-30

Similar Documents

Publication Publication Date Title
Achanta et al. Superpixels and polygons using simple non-iterative clustering
CN106127791B (en) A kind of contour of building line drawing method of aviation remote sensing image
Li et al. An edge embedded marker-based watershed algorithm for high spatial resolution remote sensing image segmentation
CN110232389B (en) Stereoscopic vision navigation method based on invariance of green crop feature extraction
CN103679714B (en) A kind of optics and SAR automatic image registration method based on gradient cross-correlation
CN105809703A (en) Adhesion hemocyte image segmentation method based on improved fractional differential and graph theory
CN102708370B (en) Method and device for extracting multi-view angle image foreground target
CN108765336A (en) Image defogging method based on dark bright primary colors priori with auto-adaptive parameter optimization
CN107369131B (en) Conspicuousness detection method, device, storage medium and the processor of image
CN105096347B (en) Image processing apparatus and method
CN107392968B (en) The image significance detection method of Fusion of Color comparison diagram and Color-spatial distribution figure
CN108537756A (en) Single image to the fog method based on image co-registration
CN109712143B (en) Rapid image segmentation method based on superpixel multi-feature fusion
CN108564528A (en) A kind of portrait photo automatic background weakening method based on conspicuousness detection
CN101826209A (en) Canny model-based method for segmenting three-dimensional medical image
CN106294705A (en) A kind of batch remote sensing image preprocess method
CN107154044A (en) A kind of dividing method of Chinese meal food image
CN108596923A (en) Acquisition methods, device and the electronic equipment of three-dimensional data
CN102034247A (en) Motion capture method for binocular vision image based on background modeling
CN103700085A (en) Cutting method of pectoral muscle region in mammary gland X-ray image
CN103020933A (en) Multi-source image fusion method based on bionic visual mechanism
CN107527054A (en) Prospect extraction method based on various visual angles fusion
CN107123130A (en) Kernel correlation filtering target tracking method based on superpixel and hybrid hash
Maltezos et al. Automatic detection of building points from LiDAR and dense image matching point clouds
CN104574374B (en) The automatic division method that retina serous pigmentary epithelial layer is detached from

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181030

CF01 Termination of patent right due to non-payment of annual fee