CN108507495B - 一种基于逆向哈特曼检测的自由曲面检测方法 - Google Patents
一种基于逆向哈特曼检测的自由曲面检测方法 Download PDFInfo
- Publication number
- CN108507495B CN108507495B CN201810225433.5A CN201810225433A CN108507495B CN 108507495 B CN108507495 B CN 108507495B CN 201810225433 A CN201810225433 A CN 201810225433A CN 108507495 B CN108507495 B CN 108507495B
- Authority
- CN
- China
- Prior art keywords
- error
- path system
- structural
- system model
- ideal
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 47
- 230000004075 alteration Effects 0.000 claims abstract description 41
- 230000003287 optical effect Effects 0.000 claims abstract description 34
- 238000005457 optimization Methods 0.000 claims abstract description 32
- 238000002474 experimental method Methods 0.000 claims abstract description 7
- 238000000926 separation method Methods 0.000 claims abstract description 7
- 238000012937 correction Methods 0.000 claims abstract description 4
- 238000000034 method Methods 0.000 claims description 11
- 238000005215 recombination Methods 0.000 claims 1
- 230000006798 recombination Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
- G01B11/254—Projection of a pattern, viewing through a pattern, e.g. moiré
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
- G01B11/005—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates coordinate measuring machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供一种基于逆向哈特曼检测的自由曲面检测方法,涉及测量技术领域。标定光路系统结构位置参数;建立被测物置为理想面的理想光路系统模型;依据理想光路系统模型,实验测得包含表面误差与结构误差的波前像差;调整理想光路系统模型参数,对测得数数据进行基于低、高阶像差分离优化的两步结构误差校正,得到标准光路系统模型;对标准光路系统模型进行光线追迹,测得只有由被测物表面误差造成的波前像差,计算得到被测物表面误差。本发明解决了现有技术中高精度检测自由曲面不具有通用性的技术问题。本发明有益效果为:对结构误差进行基于泽尼克拟合的两步优化,有效消除系统的结构位置误差,提升了系统检测精度,检测通用化。
Description
技术领域
本发明涉及测量技术领域,尤其是涉及一种凸面的自由曲面高精度的检测方法。
背景技术
对于光学自由曲面的运用,可以使运用光学原理的检测仪器拥有更高的检测精度及更简单的仪器结构。然而在成像领域中,光学元件面形的高精度要求却限制了自由曲面的大规模应用。对于自由曲面的检测技术已经成为制约其应用的重要因素。目前,在光学测量领域,高精度干涉检测法是最为普遍的方法,但其存在通用性差、成本高、检测的动态范围较小等问题。一种软件可配置光学测试系统(详见Peng Su等,Software configurableoptical test system:a computerized reverse Hartmann test,Applied Optics,2010,49(23):4404-4412)采用逆向哈特曼检验光路,实现对于内反射面形误差较大动态范围的高精度测量。但是,该方法对于设备各元件空间位置的标定有极高的精度要求。中国专利申请公布号CN107560564A,申请公布日2018年1月9日,名称为“一种自由曲面检测方法及系统”的发明专利申请文件,公开了一种凸面的自由曲面检测方法。包括:采用三坐标测量设备对包括待测反射球面、投影屏和CCD相机在内的逆向哈特曼检验光路的结构位置参数S进行测量标定;根据测量标定的结构位置参数S,确定待测反射球面的面形误差数据W0;根据测量标定的结构位置参数S和面形误差数据W0,采用泽尼克拟合确定待测球面的面形偏差优化目标;根据确定的面形偏差优化目标,确定初始测量标定的结构位置参数S的各项偏差,并根据所述各项偏差确定待测反射面的实际面形误差ΔW。该方法在结构误差消除过程中采用添加限定条件求取测得的波前数据最小,但是其对于初始测得结构参数的精度要求较高,同时,无法避免由被测元件自带的高阶像差导致的过校正现象,对于高精度的面形检测不具有通用性。
发明内容
为了解决现有技术中对自由曲面检测精度不高,而高精度检测不具有通用性的技术问题,本发明提供一种基于逆向哈特曼检测的自由曲面检测方法,实现一种具有通用性的高精度检测方法。
本发明的技术方案是:一种基于逆向哈特曼检测的自由曲面检测方法:它包括由CCD相机、投影屏和被测物构成逆向哈特曼检测光路系统,其特征在于,检测方法包括:步骤1,标定逆向哈特曼检测光路系统的结构位置参数;步骤2,依据结构位置参数,建立被测物置为理想面的理想光路系统模型;步骤3,依据理想光路系统模型,实验测得包含表面误差与结构误差的波前像差;步骤4,调整理想光路系统模型参数,对测得数数据进行基于低阶像差优化和基于高阶像差分离优化的两步结构误差校正,使得理想光路系统模型的结构位置与实际逆向哈特曼检测光路系统的结构位置基本一致,得到标准光路系统模型;步骤5,对标准光路系统模型进行光线追迹,测得只有由被测物表面误差造成的波前像差,计算得到被测物表面误差。
作为优选,步骤4中,通过定量附加结构误差至理想光路系统模型,获取对应的泽尼克系数变化量,建立各项结构误差与测得波前像差拟合而成的泽尼克系数之间的线性方程,将方程斜率作为误差对泽尼克系数的权重,依据前四项泽尼克系数的权重大小,将结构误差ΔP可分为两个部分:结构误差ΔP(H)和结构误差ΔP(L)。
作为优选,步骤4中,基于低阶像差优化,将测得波前数据拟合而成的前四项泽尼克系数重新组成的波前像差作为优化目标,迭代赋值结构误差ΔP(H),调整对应的理想光路系统模型参数,使优化目标迭代减小收敛,得到收敛结果对应的理想光路系统模型的结构位置参数P(1)。
作为优选,步骤4中,对面形数据进行基于高阶像差分离优化方法包括:通过建立相关于各项残余像差的方程组,计算获得相关于权重的残余结构误差比值rk,i,将整体结构误差构成的波前像差作为优化目标,基于残余结构误差比值rk,i,迭代赋值结构误差ΔP(L),调整对应的理想光路系统模型参数,得到理想光路系统模型的最优结构位置参数P(2)。
与现有技术相比,本发明的有益效果是:通过细分结构误差来源,并进行基于泽尼克拟合的两步优化方式对标定测量系统结构位置参数进行校正,有效消除系统的结构位置误差,提升了系统检测精度,对于面形误差的检测达到纳米级别。设备简单所占空间小、操作高效、有效减小测量过程中的成本需求以及操作的繁琐程度。对于任何数值孔径的凹凸反射球面的高精度可以进行通用化检测。
附图说明
附图1为本发明逆向哈特曼检测光路系统结构示意图;
附图2为被测物自由曲面反射后的水平x方向(a)和竖直y方向的正弦条纹(b)图;
附图3为本发明基于预标定结构参数经系统检测出的表面面形误差图;
附图4为本发明经过第一步优化以后的系统检测出的表面面形误差图;
附图5为本发明经过第二步优化之后的系统检测出的表面面形误差图。
图中:1-.CCD相机;2-投影屏;3-被测物。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例1:
如图1所示,一种基于逆向哈特曼检测的自由曲面检测方法:它包括CCD相机1、投影屏2、被测物3、调整架、三坐标测量仪、计算机。CCD相机1前端设有滤光小孔。投影屏2设有可产生水平x方向和竖直y方向的正弦条纹。CCD相机1和投影屏2在同一侧。CCD相机1的镜头与投影屏2位于同一平面。调整架可二维平移及旋转。三坐标测量仪的测量精度达微米量级。计算机设有计算和控制软件。被测物3放置在调整架上,被测物3的自由曲面朝向CCD相机1和投影屏2。通过调整调整架使得投影屏2发出的光经过被测物3被CCD相机1接收,CCD相机1能够采得完整的被测物3自由曲面的像。CCD相机1、投影屏2和被测物3构成逆向哈特曼检测光路系统。图1中带箭头的线表示光路。
检测步骤:步骤1,三坐标测量仪对逆向哈特曼检测光路系统的结构位置参数标定。三坐标测量仪标定CCD相机1、投影屏2和被测物3自由曲面固定后的结构位置参数为P。P={xi,yi;Ti,x,Ti,y,Ti,z;dc-m,ds-m}i=1,2,3。其中(xi,yi)表示第i个元件在三维空间中x和y的相对坐标;(Ti,x,Ti,y,Ti,z)表示第i个元件关于各坐标轴的倾角;dc-m表示CCD相机1与被测物3自由曲面的距离;ds-m表示投影屏幕2到被测物3自由曲面的距离;ΔP={εk}k=1,2...M表示在预标定结果中存在的预标定误差。其中εk表示第k项光路结构位置参数P上存在的结构误差,M表示光路结构位置参数P中的参数总数。
步骤2,依据结构位置参数P,建立被测物置为理想面的理想光路系统模型,如图2所示。通过计算机控制,在投影屏2上产生水平x方向和竖直y方向的正弦条纹。对应的条纹间距分别为Sx和Sy,并对其进行90°四步移相。同时,由CCD相机1实时采集移相正弦条纹经被测物3自由曲面反射后的像。分别利用四步移相算法求解集到的移相正弦条纹对应的相位分布Φx和Φy,并由此得到其投影在投影屏2上的横、纵坐标值Xm和Ym为:
Xm=Φx·Sx/2π,
Ym=Φy·Sy/2π;
利用步骤1标定得到的结构位置参数P,在计算机中建立被测物3被置为理想面的理想光路系统模型。
步骤3,依据理想光路系统模型,实验测得包含表面误差与结构误差的波前像差。由于预标定误差ΔP的存在,实验测得波前像差Wmeas可以表述为:
其中,N为泽尼克多项式的项数;Cmeas,i是实验测得数据Wmeas拟合而成的第i项泽尼克系数;Zi是第i项泽尼克多项式。
通过不断对结构位置参数附加额外的第k项结构误差εk,获取其对应引入的波前像差并对波前像差进行泽尼克拟合:其中,表示第i项泽尼克系数。对结构误差εk及其对应的泽尼克系数进行线性拟合,将两者关系表示为:其中,表示第k项结构误差εk对第i项泽尼克系数的影响权重。此时,测出的表面面形误差图如图3所示。图(包括以下的图)中:RMS表示均方根值,即表面面形数据的均方根值,可用来评估一个表面的粗糙度;PV表示峰谷值,即面形数据中最大点与最小值的差值的绝对值。
步骤4,调整理想光路系统模型参数,对测得数据进行基于低阶像差优化和基于高阶像差分离优化的两步结构误差校正,使得理想光路系统模型的结构位置与实际逆向哈特曼检测光路系统的结构位置基本一致,得到标准光路系统模型。将实验测得数据Wmeas中预标定误差ΔP引入的结构像差WΔP消除至最小,即可测得被测物3自由曲面的表面误差,故对测得的波前像差Wmeas进行基于低阶像差优化和高阶像差分离的两步优化。
对测得的波前数据Wmeas进行基于低阶像差优化:以逆向哈特曼检测光路系统的构位置参数误差ΔP作为优化变量,根据前四项泽尼克系数的权重大小,结构误差可以被分为两个部分ΔP={ΔP(H),ΔP(L)},其中ΔP(H)表示前四项泽尼克系数权重大的结构误差,ΔP(L)表示前四项泽尼克系数权重小的结构误差。将实验测得像差Wmeas前四项泽尼克系数{Cmeas,i}i=1,2,3,4(平移项,x轴倾斜项,y轴倾斜项,离焦项)所构成的波前像差Wgeo,4作为优化目标,定义偏差函数:{f(1)(ΔP(H))}=min{Wgeo,4}。利用迭代优化算法,重复调节理想光路系统模型参数,使得偏差函数{f(1)}迭代减小收敛。得到偏差函数{f(1)}收敛时对应的逆向哈特曼检测光路的结构位置参数误差与标定的结构位置参数P相加获得第一步优化结果中的结构位置参数P(1),定义其对应的结构参数误差为对基于结构位置参数P(1)建立的理想系统模型进行光线追迹,进而得到被测物3的测得值此时,测出的表面面形误差图如图4所示。
对波前数据进行基于高阶像差分离的优化:将几何像差按其来源进行分析,可以得到:
其中:M为结构误差因素总数,N为泽尼克多项式的项数。由于第一步优化后消去了低阶像差及大部分的高阶像差,故设立方程组:
从而获得相关于权重的残余结构误差比值:
其中:表示在结构参数ΔP(L)中的任一因素;表示其相关的各项系数权重。将比值条件作为限制条件,继而定义目标函数:
其中:Wgeo表示结构误差ΔP(1)引入的结构误差;表示结构误差造成的波前像差。可以表示为其中:s表示结构参数的因素数目;表示在第二步优化中测得的波前像差。利用迭代优化算法,基于残余结构误差比值rk,i,迭代赋值结构误差,调整对应的理想光路系统模型参数,使得偏差函数{f(2)}迭代减小收敛。得到偏差函数{f(2)}收敛时对应的逆向哈特曼检验光路的最优结构位置参数P(2)。
步5,对标准光路系统模型进行光线追迹,测得只有由被测物表面误差造成的波前像差,计算得到被测物表面误差。对基于结构位置参数P(2)建立的标准光路系统模型进行光线追迹。此时,标准光路系统模型的结构位置与实际逆向哈特曼检测光路系统的结构位置基本一致,结构误差被消除,得到对应投影坐标数据Xm和Ym的理论值分别为Xideal和Yideal。计算被测物3自由曲面面形误差Wsurf对应的局部斜率(wx,wy)为:
wx=(Xm-Xideal)/2ds-m,
wy=(Ym-Yideal)/2ds-m,
对局部斜率(wx,wy)积分得到的面形误差数据Wsurf,如图5所示。
Claims (1)
1.一种基于逆向哈特曼检测的自由曲面检测方法:它包括由CCD相机、投影屏和被测物构成逆向哈特曼检测光路系统,其特征在于,检测方法包括:步骤1,标定逆向哈特曼检测光路系统的结构位置参数;步骤2,依据结构位置参数,建立被测物置为理想面的理想光路系统模型;步骤3,依据理想光路系统模型,实验测得包含表面误差与结构误差的波前像差;步骤4,调整理想光路系统模型参数,对测得数据进行基于低阶像差优化和基于高阶像差分离优化的两步结构误差校正,使得理想光路系统模型的结构位置与实际逆向哈特曼检测光路系统的结构位置基本一致,得到标准光路系统模型;步骤5,对标准光路系统模型进行光线追迹,测得只有由被测物表面误差造成的波前像差,计算得到被测物表面误差;步骤4中,通过定量附加结构误差至理想光路系统模型,获取对应的泽尼克系数变化量,建立各项结构误差与测得波前像差拟合而成的泽尼克系数之间的线性方程,将方程斜率作为误差对泽尼克系数的权重,依据前四项泽尼克系数的权重大小,将结构误差ΔP可分为两个部分:结构误差ΔP(H)和结构误差ΔP(L);步骤4中,基于低阶像差优化,将测得波前数据拟合而成的前四项泽尼克系数重新组成的波前像差作为优化目标,迭代赋值结构误差ΔP(H),调整对应的理想光路系统模型参数,使优化目标迭代减小收敛,得到收敛结果对应的理想光路系统模型的结构位置参数P(1);步骤4中,对面形数据进行基于高阶像差分离优化方法包括:通过建立相关于各项残余像差的方程组,计算获得相关于权重的残余结构误差比值rk,i,将整体结构误差构成的波前像差作为优化目标,基于残余结构误差比值rk,i,迭代赋值结构误差ΔP(L),调整对应的理想光路系统模型参数,得到理想光路系统模型的最优结构位置参数P(2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810225433.5A CN108507495B (zh) | 2018-03-19 | 2018-03-19 | 一种基于逆向哈特曼检测的自由曲面检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810225433.5A CN108507495B (zh) | 2018-03-19 | 2018-03-19 | 一种基于逆向哈特曼检测的自由曲面检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108507495A CN108507495A (zh) | 2018-09-07 |
CN108507495B true CN108507495B (zh) | 2019-12-31 |
Family
ID=63377595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810225433.5A Expired - Fee Related CN108507495B (zh) | 2018-03-19 | 2018-03-19 | 一种基于逆向哈特曼检测的自由曲面检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108507495B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109870424B (zh) * | 2019-03-05 | 2021-05-11 | 中国计量大学 | 基于彩色三步移相技术的哈特曼光线追迹方法 |
CN110134987B (zh) * | 2019-03-29 | 2020-11-20 | 浙江大学 | 基于光线追迹的光学球面缺陷检测照明设计方法 |
CN110133845A (zh) * | 2019-04-26 | 2019-08-16 | 中国科学院上海光学精密机械研究所 | 一种用于激光系统的自由曲面波前补偿元件的设计方法 |
CN110763159B (zh) * | 2019-11-01 | 2021-08-06 | 中国计量大学 | 一种光学偏折显微表面测量装置及方法 |
CN110793465B (zh) * | 2019-11-07 | 2021-07-20 | 中国计量大学 | 一种微透射元件多面形大动态范围同步测量方法 |
CN110726381B (zh) * | 2019-11-22 | 2021-10-15 | 中国科学院长春光学精密机械与物理研究所 | 一种光学自由曲面全频段像差检测系统及检测方法 |
CN110702009A (zh) * | 2019-11-22 | 2020-01-17 | 中国科学院长春光学精密机械与物理研究所 | 一种基于逆向Hartmann计算机辅助法的三维测量系统 |
CN114353699B (zh) * | 2022-01-07 | 2022-12-16 | 中国科学院长春光学精密机械与物理研究所 | 大陡度凸面光学自由曲面高频段像差检测系统及检测方法 |
CN114593693B (zh) * | 2022-02-16 | 2022-11-01 | 苏州英示测量科技有限公司 | 光学偏折检测方法、电子设备及光学偏折检测系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4112072B2 (ja) * | 1998-05-12 | 2008-07-02 | 株式会社トプコン | 光学的測定装置及び光学的測定用アダプター |
CN101650157B (zh) * | 2009-09-18 | 2011-04-20 | 中国科学院长春光学精密机械与物理研究所 | 双曲面凸面反射镜面形误差的检测方法及其装置 |
CN102297663B (zh) * | 2011-04-26 | 2013-01-23 | 北京理工大学 | 基于彩色编码逆向哈特曼式面形测量方法 |
JP6112909B2 (ja) * | 2013-02-27 | 2017-04-12 | キヤノン株式会社 | シャック・ハルトマンセンサーを用いた形状計測装置、形状計測方法 |
CN104238110B (zh) * | 2014-09-19 | 2016-06-01 | 北京空间机电研究所 | 一种基于自适应光学的平行光管波前像差预补偿装置 |
CN107144419B (zh) * | 2017-04-17 | 2019-03-19 | 中国科学院光电研究院 | 一种基于夏克-哈特曼波前传感器的光学系统波像差测量装置与方法 |
CN107560564B (zh) * | 2017-07-28 | 2019-09-10 | 中国计量大学 | 一种自由曲面检测方法及系统 |
-
2018
- 2018-03-19 CN CN201810225433.5A patent/CN108507495B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN108507495A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108507495B (zh) | 一种基于逆向哈特曼检测的自由曲面检测方法 | |
Isheil et al. | Systematic error correction of a 3D laser scanning measurement device | |
CN109307480B (zh) | 一种透射元件多表面面形检测方法 | |
JP5971965B2 (ja) | 面形状計測方法、面形状計測装置、プログラム、および、光学素子の製造方法 | |
CN102095385B (zh) | 新型球面绝对测量系统及方法 | |
CN111536896B (zh) | 一种激光干涉面形检测自动检测装置与方法 | |
Sims-Waterhouse et al. | Uncertainty model for a traceable stereo-photogrammetry system | |
US10323938B2 (en) | Method for calibrating a measuring device | |
CN107560564A (zh) | 一种自由曲面检测方法及系统 | |
CN107957251B (zh) | 基于计算机辅助校正的反射球面通用化检测方法 | |
CN112066913B (zh) | 一种测量球面曲率半径的方法 | |
Kudo et al. | Simulation-based systematic error compensation for nanoprofiler using normal vector tracing method | |
US20220146370A1 (en) | Deflectometry devices, systems and methods | |
CN110428471B (zh) | 一种针对光学自由曲面子孔径偏折测量的精确自定位方法 | |
KR20110065365A (ko) | 비구면체 측정 방법 및 장치 | |
CN108917662B (zh) | 参考面平面度检验的优化方法 | |
King et al. | Development of a metrology workstation for full-aperture and sub-aperture stitching measurements | |
WO2020049911A1 (ja) | ウェーハ形状データ化方法 | |
CN110793465B (zh) | 一种微透射元件多面形大动态范围同步测量方法 | |
JP2000097663A (ja) | 干渉計 | |
CN112629436B (zh) | 一种基于自适应光学波前校正的高次非球面检测方法 | |
TW202246749A (zh) | 使用校正板校正光學系統之像差及跡變 | |
JP6139950B2 (ja) | 計測装置、計測方法及び原器 | |
CN114858091B (zh) | 一种同时适用于平面和球面的标定回程误差的方法 | |
CN114076573B (zh) | 等效元件、等效元件的制备方法、检测精度校验方法 |
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 | ||
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: 20191231 |