CN102937592A - Ceramic radome pore and material loosening defect automatic detection method - Google Patents

Ceramic radome pore and material loosening defect automatic detection method Download PDF

Info

Publication number
CN102937592A
CN102937592A CN2012104013711A CN201210401371A CN102937592A CN 102937592 A CN102937592 A CN 102937592A CN 2012104013711 A CN2012104013711 A CN 2012104013711A CN 201210401371 A CN201210401371 A CN 201210401371A CN 102937592 A CN102937592 A CN 102937592A
Authority
CN
China
Prior art keywords
gray
defect area
area
ceramic radome
value
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
CN2012104013711A
Other languages
Chinese (zh)
Other versions
CN102937592B (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.)
Shandong University of Technology
Original Assignee
Shandong 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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201210401371.1A priority Critical patent/CN102937592B/en
Publication of CN102937592A publication Critical patent/CN102937592A/en
Application granted granted Critical
Publication of CN102937592B publication Critical patent/CN102937592B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to a ceramic radome pore and material loosening defect automatic detection method. The method is characterized in that the detection is carried out with a transmission image processing method. The method comprises the specific steps that: (1) a radome standard grey scale value is obtained; (2) transmission images of different positions of the ceramic radome are obtained; (3) defect area boundary extraction processing is carried out; and (4) a mean gray value and an area of a defect area are calculated; the obtained defect area mean gray value is compared with a mean gray value of the entire image, such that the material loosening state of the defect area can be determined; and the quality of the ceramic radome can be determined with the cooperation of the area value of the defect area. According to the invention, with the image processing automatic detection method, ceramic radome pore and material loosening defect is detected. The detection result is more precise and objective.

Description

陶瓷天线罩气孔及材质疏松缺陷自动检测方法Automatic detection method for porosity and material loose defects of ceramic radome

技术领域 technical field

本发明涉及一种陶瓷天线罩气孔及材质疏松缺陷自动检测方法,属于无损检测技术领域。The invention relates to an automatic detection method for air holes and material loose defects of a ceramic radome, belonging to the technical field of non-destructive testing.

背景技术 Background technique

薄壁陶瓷天线罩产品生产后不可避免的存在气孔及材质疏松缺陷,但在目前的薄壁陶瓷天线罩产品检测中,对气孔及材质疏松缺陷的检测是通过人工观察进行测量的,这种测量方法得到的结果极不准确,无法保证产品的检测质量。After the production of thin-walled ceramic radome products, there are inevitably pores and loose materials. However, in the current inspection of thin-walled ceramic radome products, the detection of pores and loose materials is measured by manual observation. This measurement The result obtained by the method is extremely inaccurate, and the detection quality of the product cannot be guaranteed.

发明内容 Contents of the invention

本发明的目的是提供一种能够克服上述缺陷、检测精确、效率高的陶瓷天线罩气孔及材质疏松缺陷自动检测方法。其技术方案为:The object of the present invention is to provide an automatic detection method for ceramic radome air holes and material loose defects that can overcome the above defects, detect accurately and have high efficiency. Its technical solution is:

一种陶瓷天线罩气孔及材质疏松缺陷自动检测方法,其特征在于采用以下步骤:1)天线罩标准灰度值的获取;2)陶瓷天线罩透射图像获取;3)缺陷区边界提取处理;4)计算缺陷区灰度平均值及缺陷区域的面积,得到的缺陷区灰度值与整幅图像灰度平均值进行比较,用于判断缺陷区的材质疏松情况,结合缺陷区的面积值的大小,来判定陶瓷天线罩的质量。An automatic detection method for ceramic radome pores and loose material defects, characterized in that the following steps are adopted: 1) Acquisition of the standard gray value of the radome; 2) Acquisition of the transmission image of the ceramic radome; 3) Extraction and processing of the boundary of the defect area; 4 ) Calculate the average gray value of the defect area and the area of the defect area, and compare the obtained gray value of the defect area with the average gray value of the entire image to judge the looseness of the material of the defect area, combined with the size of the area value of the defect area , to determine the quality of the ceramic radome.

所述陶瓷天线罩气孔及材质疏松缺陷自动检测方法,步骤1)中,在陶瓷天线罩内部设置光源,在陶瓷天线罩外部采用CCD摄像机自动进行图像采集,获得陶瓷天线罩不同部位的灰度透射图像,连续采集不同部位的多幅灰度图像,对每幅图像进行灰度统计,计算出所有图像的总灰度值,统计每幅图像的像素值,计算得到所有图像的总像素数,由总灰度值除以总像素数得到所有图像的平均灰度值,作为天线罩标准灰度值。In the method for automatic detection of pores and loose material defects of ceramic radome, in step 1), a light source is set inside the ceramic radome, and a CCD camera is used to automatically collect images outside the ceramic radome to obtain grayscale transmission of different parts of the ceramic radome Image, continuously collect multiple grayscale images of different parts, perform grayscale statistics on each image, calculate the total grayscale value of all images, count the pixel values of each image, and calculate the total number of pixels of all images, by The total gray value is divided by the total number of pixels to get the average gray value of all images, which is used as the standard gray value of the radome.

所述陶瓷天线罩气孔及材质疏松缺陷自动检测方法,步骤2)中,在陶瓷天线罩内部设置光源,采用CCD摄像机在天线罩外部对天线罩全部的不同部位进行图像采集,得到天线罩不同部位灰度透射图像。In the method for automatic detection of air holes and loose material defects of the ceramic radome, in step 2), a light source is set inside the ceramic radome, and a CCD camera is used to collect images of all different parts of the radome outside the radome to obtain different parts of the radome Grayscale transmission image.

所述陶瓷天线罩气孔及材质疏松缺陷自动检测方法,步骤3)中,为了对气孔及材质疏松缺陷区进行面积统计,需要先获得缺陷区域的边界,所以要进行边界提出处理,在灰度图像上缺陷区域和背景区域的灰度有较大的差别,缺陷区域的边界是像素的灰度值不连续的结果,所以可以通过求导的方式将其检测出来,采用拉普拉斯算子进行处理,对一个连续函数f(x,y),它在位置(x,y)处的拉普拉斯值(即二阶导数)定义为:In the automatic detection method for air holes and loose material defects of the ceramic radome, in step 3), in order to perform area statistics on the air holes and loose material defect areas, it is necessary to first obtain the boundary of the defect area, so the boundary processing must be performed. In the grayscale image There is a big difference between the gray level of the upper defect area and the background area. The boundary of the defect area is the result of the discontinuity of the gray value of the pixel, so it can be detected by derivation, and the Laplacian operator is used to perform Processing, for a continuous function f(x,y), its Laplace value (ie second derivative) at position (x,y) is defined as:

ΔΔ 22 ff == ∂∂ 22 ff ∂∂ xx 22 ++ ∂∂ 22 ff ∂∂ ythe y 22

在数字图像中用差分来近似,它的一个表示式为:Approximated by difference in digital images, one of its expressions is:

Δ2f=f(m+1,n)+f(m-1,n)+f(m,n+1)+f(m,n-1)-4f(m,n)Δ 2 f=f(m+1,n)+f(m-1,n)+f(m,n+1)+f(m,n-1)-4f(m,n)

其模板为:Its template is:

00 11 00 11 -- 44 11 00 11 00

模板所有系数之和为0,即,如果在模板各相应位置处f(m,n)的值相同(没有边界),则算子的响应为0。The sum of all coefficients of the template is 0, that is, if the value of f(m,n) at each corresponding position of the template is the same (no boundary), the response of the operator is 0.

所述陶瓷天线罩气孔及材质疏松缺陷自动检测方法,步骤4)中,在获得缺陷区边界之后,统计得到每个缺陷区像素的灰度值之和及总像素数,然后由灰度值之和分别除以总像素数,得到缺陷区的灰度平均值,由统计得到的各缺陷区的总像素数得到缺陷区域的面积,缺陷区域的灰度平均值设为T2,步骤1)中计算的整幅图像的灰度平均值设为T1,则材质疏松率δ的计算公式如下:In the method for automatic detection of pores and loose material defects of ceramic radome, in step 4), after obtaining the boundary of the defect area, the sum of the gray value of each defect area pixel and the total number of pixels are statistically obtained, and then the sum of the gray value of each defect area and respectively divided by the total number of pixels to obtain the average gray level of the defect area, and the area of the defect area is obtained from the total number of pixels of each defect area obtained by statistics, and the average gray level of the defect area is set as T 2 , in step 1) The calculated average gray value of the entire image is set to T 1 , then the calculation formula of the material porosity rate δ is as follows:

δδ == TT 22 -- TT 11 TT 11 ×× 100100 %%

由材质疏松率δ和缺陷的面积值来判定陶瓷天线罩的质量。The quality of the ceramic radome is judged by the material porosity δ and the area value of the defect.

本发明与现有技术相比,其优点是:采用了图像处理自动检测方法来对陶瓷天线罩的气孔及材质疏松进行检测,测量结果更加精确、客观。Compared with the prior art, the present invention has the advantages that an image processing automatic detection method is adopted to detect the pores and material looseness of the ceramic radome, and the measurement result is more accurate and objective.

具体实施方式 Detailed ways

对一个待检测的陶瓷天线罩进行气孔及材质疏松缺陷检测,在陶瓷天线罩内部设置光源,在陶瓷天线罩外部采用CCD摄像机自动进行图像采集,获得陶瓷天线罩不同部位的灰度透射图像,连续采集不同部位的多幅灰度图像,对每幅图像进行灰度统计,计算出所有图像的总灰度值,统计每幅图像的像素值,计算得到所有图像的总像素数,由总灰度值除以总像素数得到所有图像的平均灰度值,作为天线罩标准灰度值,这里采集了20幅图像,计算得到的天线罩标准灰度值:T1=80。For a ceramic radome to be inspected, the porosity and loose material defects are detected. A light source is set inside the ceramic radome, and a CCD camera is used to automatically collect images outside the ceramic radome to obtain grayscale transmission images of different parts of the ceramic radome. Continuous Collect multiple grayscale images of different parts, perform grayscale statistics on each image, calculate the total grayscale value of all images, count the pixel values of each image, and calculate the total number of pixels of all images, from the total grayscale The value is divided by the total number of pixels to obtain the average gray value of all images, which is used as the standard gray value of the radome. Here, 20 images are collected, and the calculated standard gray value of the radome is: T 1 =80.

在陶瓷天线罩内部设置光源,采用CCD摄像机在天线罩外部对天线罩全部的不同部位进行图像采集,得到天线罩不同部位灰度透射图像,然后依次经过以下图像处理步骤:A light source is set inside the ceramic radome, and a CCD camera is used to collect images of all different parts of the radome outside the radome to obtain grayscale transmission images of different parts of the radome, and then undergo the following image processing steps in sequence:

步骤1):在采集得到灰度图像之后,为了对气孔及材质疏松缺陷区进行面积统计,需要先获得缺陷区域的边界,所以要进行边界提出处理,在灰度图像上缺陷区域和背景区域的灰度有较大的差别,缺陷区域的边界是像素的灰度值不连续的结果,所以可以通过求导的方式将其检测出来,采用拉普拉斯算子进行处理,对一个连续函数f(x,y),它在位置(x,y)处的拉普拉斯值(即二阶导数)定义为:Step 1): After the grayscale image is collected, in order to count the area of the air hole and loose material defect area, the boundary of the defect area needs to be obtained first, so the boundary processing must be performed, and the defect area and the background area on the grayscale image There is a large difference in grayscale. The boundary of the defect area is the result of the discontinuous grayscale value of the pixel, so it can be detected by derivation and processed by the Laplacian operator. For a continuous function f (x,y), its Laplace value (ie second derivative) at position (x,y) is defined as:

ΔΔ 22 ff == ∂∂ 22 ff ∂∂ xx 22 ++ ∂∂ 22 ff ∂∂ ythe y 22

在数字图像中用差分来近似,它的一个表示式为:Approximated by difference in digital images, one of its expressions is:

Δ2f=f(m+1,n)+f(m-1,n)+f(m,n+1)+f(m,n-1)-4f(m,n)Δ 2 f=f(m+1,n)+f(m-1,n)+f(m,n+1)+f(m,n-1)-4f(m,n)

其模板为:Its template is:

00 11 00 11 -- 44 11 00 11 00

注意,模板所有系数之和为0,即,如果在模板各相应位置处f(m,n)的值相同(没有边界),则算子的响应为0。Note that the sum of all coefficients of the template is 0, that is, if the value of f(m,n) at each corresponding position of the template is the same (no boundary), the response of the operator is 0.

步骤2):在获得缺陷区边界之后,统计得到缺陷区像素的灰度值之和及总像素数,然后缺陷区像素的灰度值之和除以总像素数,得到灰度平均值为T2=140,由统计得到的像素数得到每个缺陷区域的面积,面积为S=28mm2,则材质疏松率δ的计算公式如下:Step 2): After obtaining the boundary of the defect area, the sum of the gray values of the pixels in the defect area and the total number of pixels are obtained statistically, and then the sum of the gray values of the pixels in the defect area is divided by the total number of pixels to obtain the average gray value T 2 = 140, the area of each defect area is obtained from the number of pixels obtained by statistics, and the area is S = 28mm 2 , then the calculation formula of material porosity δ is as follows:

δδ == 140140 -- 8080 8080 ×× 100100 %% == 7575 %%

由材质疏松率δ和缺陷的面积值S来判定陶瓷天线罩的质量。The quality of the ceramic radome is judged by the material porosity δ and the area value S of the defect.

Claims (5)

1. a ceramic radome pore and material rarefaction defect automatic testing method is characterized in that adopting image processing method to detect, and concrete steps are: 1) antenna house standard grayscale value obtains; 2) ceramic radome different parts transmission image obtains; 3) the defect area Boundary Extraction is processed; 4) area of calculating defect area average gray and defect area, the defect area gray-scale value and the entire image average gray that obtain compare, be used for judging pore and the loose situation of material of defect area, the area value in binding deficient district is big or small, judges the quality of ceramic radome.
2. ceramic radome pore as claimed in claim 1 and material rarefaction defect automatic testing method, it is characterized in that: in the step 1), in ceramic radome inside light source is set, automatically carry out image acquisition at the outside ccd video camera that adopts of ceramic radome, obtain the gray scale transmission image of ceramic radome different parts, several gray level images of continuous acquisition different parts, every width of cloth image is carried out gray-scale statistical, calculate total gray-scale value of all images, add up the pixel value of every width of cloth image, calculate the total pixel number of all images, obtain the average gray value of all images by total gray-scale value divided by total pixel number, as antenna house standard grayscale value.
3. ceramic radome pore as claimed in claim 1 and material rarefaction defect automatic testing method, it is characterized in that: step 2) in, in ceramic radome inside light source is set, adopt ccd video camera in the antenna house outside the whole different parts of antenna house to be carried out image acquisition, obtain antenna house different parts gray scale transmission image.
4. ceramic radome pore as claimed in claim 1 and material rarefaction defect automatic testing method, it is characterized in that: in the step 3), the gray scale of defect area and background area has larger difference on gray level image, the border of defect area is the discontinuous result of the gray-scale value of pixel, so can it be detected by the mode of differentiate, adopt Laplace operator to process, to a continuous function f (x, y), Laplce's value that it is located at position (x, y) (being second derivative) is defined as:
Δ 2 f = ∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2
Be similar to difference in digital picture, its expression is:
Δ 2f=f(m+1,n)+f(m-1,n)+f(m,n+1)+f(m,n-1)-4f(m,n)
Its template is:
0 1 0 1 - 4 1 0 1 0
All coefficient sums of template are 0, that is, if in the value identical (not having the border) of each corresponding position f (m, n) of template, then the response of operator is 0.
5. ceramic radome pore as claimed in claim 1 and material rarefaction defect automatic testing method, it is characterized in that: in the step 4), after obtaining the defect area border, statistics obtains gray-scale value sum and the total pixel number of each defect area pixel, then by the gray-scale value sum respectively divided by total pixel number, obtain the average gray of defect area, the total pixel number of each defect area that is obtained by statistics obtains the area of defect area, and the average gray of defect area is made as T 2, the average gray of the entire image of calculating in the step 1) is made as T 1, then the computing formula of the loose rate δ of material is as follows:
δ = T 2 - T 1 T 1 × 100 %
Judged the quality of ceramic radome by the area value of the loose rate δ of material and defective.
CN201210401371.1A 2012-10-20 2012-10-20 Ceramic radome pore and material loosening defect automatic detection method Expired - Fee Related CN102937592B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210401371.1A CN102937592B (en) 2012-10-20 2012-10-20 Ceramic radome pore and material loosening defect automatic detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210401371.1A CN102937592B (en) 2012-10-20 2012-10-20 Ceramic radome pore and material loosening defect automatic detection method

Publications (2)

Publication Number Publication Date
CN102937592A true CN102937592A (en) 2013-02-20
CN102937592B CN102937592B (en) 2014-07-16

Family

ID=47696503

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210401371.1A Expired - Fee Related CN102937592B (en) 2012-10-20 2012-10-20 Ceramic radome pore and material loosening defect automatic detection method

Country Status (1)

Country Link
CN (1) CN102937592B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527301A (en) * 2016-02-04 2016-04-27 中达电通股份有限公司 Detection system and detection method for voltage grade of electric heating pipe
CN113029504A (en) * 2021-03-04 2021-06-25 中国航空工业集团公司西安航空计算技术研究所 Quantitative detection system and method for cooling air stagnation area of low-profile-rate gradually-expanding channel
CN116664846A (en) * 2023-07-31 2023-08-29 华东交通大学 Method and system for monitoring construction quality of 3D printed bridge deck based on semantic segmentation

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0469552A (en) * 1990-07-10 1992-03-04 Ngk Spark Plug Co Ltd Nondestructive inspecting method for ceramic member
JPH11326225A (en) * 1998-05-11 1999-11-26 Hitachi Ltd Inspection method and apparatus for tubular member
US20040086168A1 (en) * 2002-10-23 2004-05-06 Masayuki Kuwabara Pattern inspection method and inspection apparatus
JP2006267022A (en) * 2005-03-25 2006-10-05 Kurabo Ind Ltd Defect inspection apparatus and method
CN1958933A (en) * 2006-11-24 2007-05-09 东华大学 Method and system for measuring porosity of not weaving fabric
JP4121186B2 (en) * 1997-06-16 2008-07-23 株式会社日立製作所 Ceramic tube inspection method and apparatus
CN101320004A (en) * 2008-07-03 2008-12-10 西北工业大学 Bamboo strip defect online detection method based on machine vision
CN102645436A (en) * 2012-04-09 2012-08-22 天津大学 Engineering ceramic grinding surface damage detection method based on grinding grain removal technology

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0469552A (en) * 1990-07-10 1992-03-04 Ngk Spark Plug Co Ltd Nondestructive inspecting method for ceramic member
JP4121186B2 (en) * 1997-06-16 2008-07-23 株式会社日立製作所 Ceramic tube inspection method and apparatus
JPH11326225A (en) * 1998-05-11 1999-11-26 Hitachi Ltd Inspection method and apparatus for tubular member
US20040086168A1 (en) * 2002-10-23 2004-05-06 Masayuki Kuwabara Pattern inspection method and inspection apparatus
JP2006267022A (en) * 2005-03-25 2006-10-05 Kurabo Ind Ltd Defect inspection apparatus and method
CN1958933A (en) * 2006-11-24 2007-05-09 东华大学 Method and system for measuring porosity of not weaving fabric
CN101320004A (en) * 2008-07-03 2008-12-10 西北工业大学 Bamboo strip defect online detection method based on machine vision
CN102645436A (en) * 2012-04-09 2012-08-22 天津大学 Engineering ceramic grinding surface damage detection method based on grinding grain removal technology

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
沙翠翠: "基于机器视觉的陶瓷管缺陷识别技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
田欣利等: "一种基于表面图像灰度的工程陶瓷粗糙度研究", 《机械科学与技术》 *
胡尚举等: "边缘检测算子的分析比较", 《大众科技》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105527301A (en) * 2016-02-04 2016-04-27 中达电通股份有限公司 Detection system and detection method for voltage grade of electric heating pipe
CN105527301B (en) * 2016-02-04 2018-08-14 中达电通股份有限公司 The detecting system and detection method of the voltage class of electric heating tube
CN113029504A (en) * 2021-03-04 2021-06-25 中国航空工业集团公司西安航空计算技术研究所 Quantitative detection system and method for cooling air stagnation area of low-profile-rate gradually-expanding channel
CN113029504B (en) * 2021-03-04 2023-08-04 中国航空工业集团公司西安航空计算技术研究所 Quantitative detection system and method for cooling air stagnation area of low-profile gradually-expanding channel
CN116664846A (en) * 2023-07-31 2023-08-29 华东交通大学 Method and system for monitoring construction quality of 3D printed bridge deck based on semantic segmentation
CN116664846B (en) * 2023-07-31 2023-10-13 华东交通大学 Method and system for realizing 3D printing bridge deck construction quality monitoring based on semantic segmentation

Also Published As

Publication number Publication date
CN102937592B (en) 2014-07-16

Similar Documents

Publication Publication Date Title
CN108827181B (en) A Vision-Based Plate Surface Detection Method
CN103954542B (en) Based on the PM2.5 concentration detectors that non-reference picture definition is evaluated
CN103471910B (en) A kind of elongation at break of metal material intelligent test method followed the tracks of based on random point
CN102621056B (en) Method for quantitatively measuring rusting ratio of surface of metal material
CN109900711A (en) Workpiece, defect detection method based on machine vision
CN104655644A (en) Method and device for automatically detecting defects of lithium battery pole piece
CN106197612B (en) A liquid level detection method for transparent bottles based on machine vision
JP2013238595A5 (en)
CN106247969A (en) A deformation detection method of industrial magnetic core components based on machine vision
CN103927750A (en) Detection method of checkboard grid image angular point sub pixel
JP2007149837A5 (en)
CN105083912A (en) Conveyer belt deviation detecting method based on image identification
KR100826153B1 (en) Crack width measurement method using image histogram
CN104976959B (en) A kind of spring sizes on-line measurement system and its method based on machine vision
CN106355579A (en) Defect detecting method of cigarette carton surface wrinkles
CN103792705A (en) Detecting method and detecting device for detecting substrate defects
CN110084787B (en) Synthetic leather gluing online monitoring method based on machine vision
CN104156913A (en) Automatic optic inspection method for surface defects of metal cylindrical workpieces
CN101685002A (en) Method of length detection and system thereof
CN106018179A (en) Glue material surface viscosity measurement method and device based on image processing
CN105354598B (en) A kind of method of the continuous band-shaped foam metal material plating leakage defect of on-line checking
CN102937592B (en) Ceramic radome pore and material loosening defect automatic detection method
TWI548871B (en) Plate glass inspection apparatus, plate glass inspection method, plate glass manufacturing apparatus, and plate glass manufacturing method
CN101807297A (en) Medical ultrasonic image line detection method
CN108280825A (en) A kind of liquid crystal display emargintion detection method

Legal Events

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

Granted publication date: 20140716

Termination date: 20211020