CN1271086A - 测量物体表面形貌的纹影仪 - Google Patents

测量物体表面形貌的纹影仪 Download PDF

Info

Publication number
CN1271086A
CN1271086A CN 00107371 CN00107371A CN1271086A CN 1271086 A CN1271086 A CN 1271086A CN 00107371 CN00107371 CN 00107371 CN 00107371 A CN00107371 A CN 00107371A CN 1271086 A CN1271086 A CN 1271086A
Authority
CN
China
Prior art keywords
schlieren
spectroscope
video camera
light
ccd video
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
CN 00107371
Other languages
English (en)
Other versions
CN1114090C (zh
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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN 00107371 priority Critical patent/CN1114090C/zh
Publication of CN1271086A publication Critical patent/CN1271086A/zh
Application granted granted Critical
Publication of CN1114090C publication Critical patent/CN1114090C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于物体表面形貌测量技术领域,包括光源部件、光学成像部件,光源部件由激光器、小孔滤波器和扩束透镜组成,光学成像部件由分光镜、反射镜、准直透镜,以及空间滤波器和CCD摄像机组成,所说的空间滤波器为一个二维朗奇光栅。本发明可同时获取x、y两方向上的面形变化梯度信息,从而简化了整个装置的结构,不仅能测量物体静态表面的形貌,而且也可测量变化表面的形貌。

Description

测量物体表面形貌的纹影仪
本发明属于物体表面形貌测量技术领域,尤其涉及测量物体缓变表面形貌的纹影仪。
纹影测量法通过光投影的明暗变化来指示被测物体表面不规则性或内部条件所产生的折射率、反射率以及透射率的变化,因此常被用来进行流体流速场分析、溶液浓度分析、表面平直度分析等。已有的一种用于测量反射表面的纹影仪包括:一点光源、一分光镜、一准直镜以及一刀口空间滤波器。点光源发出的光经分光镜分光后,一路经准直镜准直后照射在被测物体表面上,另一路则会聚后成像在观察屏上。从被测表面反射回来的光束再经分光镜后会聚在像平面上,若被测表面是一平面的话,反射光束将在像平面上会聚成一无限小的光斑。若反射表面是一不规则的表面,其成像则或高于、或低于焦点。在焦点处放置一刀口,则该刀口会截断一部分的反射光束。被阻断的那部分光会在观察屏上形成暗条纹,而未被阻断的那部分光则在观察屏上形成亮条纹,从而在观察屏上便形成明暗相间的干涉图形。从该图形便可来检测被测表面的部分参数。从本质上讲,纹影仪所检测的是表面高度变化在一个方向(x)上的梯度变化,为要得到另一方面(y)上的变化,要么被测物体、要么刀口需转动90度,因此x、y方向上的梯度变化必须分开进行测量,且其方向不能确定。这种刀口式纹影仪法的最大缺点是:它是一种半定量的测量方法。
对物体表面形貌的测量通常采用接触与非接触两类,在非接触测量中以光干涉和光线偏折法最为典型。干涉法测量的是光程差,而干涉图上各点的光强与光程差是一对多的关系,不能根据一张干涉图来完全确定被测表面的形貌,因此常在干涉法中引入相移技术,用多张干涉图来计算各点的相位,最终重构被测表面轮廓。缓变表面形貌的测量也是一个重要的测量课题,一些晶体在其生长过程中常需要研究和观察其浓度的变化以及表面形貌的宏观变化。因此相移干涉法不适合于测量变化的表面形貌。纹影仪法和莫尔测偏法由于测量的是光线偏折角,它不但包含了该点斜率的大小信息,而且也包含了斜率符号信息,因此适合于测量动态表面的形貌,但需同时获取两方向的表面梯度变化信息。而如前所述,无论是纹影法还是莫尔测偏法均不能同时获取两方向的梯度信息。
在W.L.豪维斯,“彩虹式放影仪与马赫一曾德尔干涉仪的比较”,《应用光学》,1985年,第24卷,第6期,第816-822页中描述了一种用于检测被测物表面平直度的纹影仪,其中的空间滤波器用一块颜色盘来替代刀口。盘的中央区为透明,第一圈具有第一种颜色,第二圈为第二种颜色,其他的圈各具有不同的颜色。若被测表面为平面,则观察屏上的像为白色;反之则为不同的颜色。尽管用该法可同时检测被测表面在不同径向上的梯度变化,但该法的缺点是不能敏感被测表面的切向梯度变化,因此是一种半定量的分析方法。
美国专利第4,854,708号中介绍了一种莫尔测偏法,它本质上也是一种纹影测量法,其基本装置仍是一纹影仪的装置,只是替代在像平面上的刀口装置,而在像平面之后放置两块平行的朗奇光栅,两块光栅相距一泰伯距离d,且相互倾斜一角度0,从而当光束经过两块光栅后,则在观察屏上形成明暗相间的莫尔条纹。由此来测量被测表面不规则性或位相物体的内部特征。对反射表面来说该法测量的量本质上是表面起伏高度的一阶偏导数,因此为得到两方向上的导数,需将两光栅或被测物旋转90°。因此对表面测量带来一定的不便之处,尤其是该方法不能用于测量变化介质或变化的表面。
在S.克莱恩等,“用于现场定量影照生长晶体周围溶液浓度轮廓的暗场式纹影显微镜”,《晶体生长杂志》1997年,第179期,第240-248页中介绍了一种定量测量生长晶体周围溶液浓度轮廓的暗场型纹影仪显微镜,该法的基本装置为一纹影仪,在焦平面上同样设置一刀口。该法的缺点也是不能同时获取x、y方向的浓度梯度,因此不能用于较快变化的扩散的测量,
本发明的目的是为克服已有技术的不足之处,提供一种测量物体表面形貌的纹影仪,可同时获取x、y两方向上的面形变化梯度信息,从而简化了整个装置的结构,不仅能测量物体静态表面的形貌,而且也可测量变化表面的形貌。
本发明提出的一种测量物体表面形貌的纹影仪,包括光源部件、光学成像部件,所说的光源部件由激光器及设置在其输出光束上的小孔滤波器和扩束透镜组成,所说的光学成像部件由将经扩束透镜扩束后的光束进行分束的分光镜、接收该分光镜反射的光束的反射镜、接收该反射镜反射的光束的准直透镜,以及接收该分光镜透射的光束的空间滤波器和CCD摄像机组成,其特征在于,所说的空间滤波器为一个二维朗奇光栅。
所说的激光器可为He-Ne激光器或半导体激光器。所说的准直透镜可采用消像差非球面准直镜。所说的二维朗奇光栅可被置于焦平面上;也可被置于焦平面之前的适当位置上。所说的CCD摄像机可为CCD面阵摄像机。本发明还可包括还包括用于对图像信号进行处理的图像采集与处理***,该***由用于采集所说的CCD摄像机获取被测表面形貌的信息的高速图像采集卡和对该图像数据进行处理的微机所构成,所说的高速图像采集卡与该微机的接口相连,所说的微机储存有图形处理及控制软件。
本发明的主要特点是在传统的纹影仪像平面的空间滤波器位置上替代刀口法而放置一块二维朗奇光栅,用于获取在x、y两方向上面形变化的梯度信息,并可进一步采用高速图像采集卡对图形信息进行实时采样,从而来重构被测表面的形貌,本发明的装置结构简单,不仅可测量物体静止表面的形貌,而且也可测量动态缓变表面的形貌。
附图简要说明
图1为本发明的测量方法原理图,
图2为本发明的测量装置的原理图,
图3为本发明的测量装置的实施例结构示意图,
图4为本实施例的两维朗奇光栅结构示意图,
图5为本实施例的图形处理软件流程图。
本发明所述的测量物体表面形貌的纹影仪的构成、工作原理及其装置的结构实施例结合各附图详细说明如下。
本发明基于朗奇光栅的测量原理如图1所示,对光线偏转角的测量是通过用光栅对光线进行采样来实现的。被测表面31为平面时,反射光线经准直透镜32、朗奇光栅33会聚于点F,根据三角关系有 A ′ B ′ ‾ AB ‾ = B ′ C ′ ‾ BC ‾ , 设朗奇光栅33的占空比为50%,则在CCD摄像机34面阵上得到的是占空比为50%的方波信号,可写为 I ( x ) = b 0 + Σ k = 1 ∞ b k sin kωx 。当被测表面有一倾斜角θ时,反射光线投射到焦平面上P点,与CCD面阵交于D点,无偏差时光线与CCD面阵交于D点,此时可得Δ= D′D″的偏差,根据图示几何关系,有: E ′ E ′ ′ ‾ P ′ P ′ ′ ‾ = D ′ D ′ ′ ‾ F ′ P ‾ ,
则有Δ= D′D″=Ztg                  (1)
由于=2θ,则
Δ=Ztg2θ                             (2)于是光栅投射就把被测表面的倾斜角θ转换为光栅投影中条纹间距的变化量Δ。为求θ角,采用同步解调技术。
将光栅投影信号表示为 I ( x ) = b 0 + Σ k = 1 ∞ b k sin [ kωx + ψ k ( x ) ] 式中ψk(x)为光栅间距变化引起的空间频率相角变化,将I(x)经低通滤波得基频信号I1(x)=b1sin[ωx+ψ1(x)],                 (3)式中ψ1(x)在采样点xi的值为 ψ 1 ( x 1 ) = 2 π d Δ ( x 1 ) - - - - ( 4 ) 式中d为光栅投影周期。将式(3)归一化得:
I1(x)=sin[ωx+ψ1(x)]                 (5)取同步信号Isync(x)=cosωx                 (6)将式(5)与式(6)相乘得: I 1 ( x ) I sync ( x ) = 1 2 { sin [ 2 ωx + ψ 1 ( x ) ] + sin [ ψ 1 ( x ) ] } - - - - ( 7 ) 将乘积经低通滤波去除2倍频的载波信号后得 F ( x ) = 1 2 sin [ ψ 1 ( x ) ] - - - - ( 8 ) 将式(2)和式(4)代入式(8)可得: 2 π d ztg 2 θ = arcsin [ 2 F ( x ) ] 亦即 θ ( x ) = 1 2 { d 2 πZ arcsin [ 2 F ( x ) ] } - - - - ( 9 ) 由此可得一行的光栅投影信号所对应的被测表面沿x方向的倾斜角。对整个被测表面则有 θ ( x , y ) = 1 2 arctg { d 2 πZ arcSin [ 2 F ( x , y ) } - - - - ( 10 ) 同理可得y方向的表达式同式(10)的倾斜角。在获得破测表面的两正交方向的倾斜角θx和θy之后,即可求出其斜率tg[θx(x,y)]和tg[θy(x,y)],进而通过积分即可求出被测表面各点的相对高度值,即表面形貌。
本发明的测量装置的组成及原理如图2所示,由He-Ne激光器1发出的激光束由扩束透镜2扩束,经小孔滤波器3消除杂散光之后,由分光镜4改变其一部分光的传播方向,经平面反射镜7改变光路,然后经非球面准直透镜8准直后照射在被测液面9上,该平行光被被测液面反射后,按原路返回经分光镜4再次将光分为两路,其中一路会聚在焦平面上,经空间滤波器5投影到CCD摄像机6的感光阵列上。
本发明的一种实施例的总体结构如图3所示,由上套筒12、下套筒13、横套筒14、上立柱15、下立柱16、中间板17、底板18、激光器组件19、以及导轨20所组成。图2所示的各光路元件分别安装在上套筒12、下套筒13、横套筒14中。上立柱15、下立柱16、中间板17用于连接上套筒12和下套筒13′使之成一整体。CCD摄像机6被放置在导轨20上。用调节螺钉21来调节被测物的水平状态。整个仪器用固定螺钉11加以固定。
由于被测量对象为缓变表面,因此要同时获取表面在x、y两方向上的变化信号,为此本发明采用两维朗奇光栅来作为空间滤波器,一种两维朗奇光栅的实施例结构如图4所示,该光栅为矩形网格状、点的相间的两维直光栅。在本发明的一个实施例中,所测表面为缓慢变化的晶体生长表面,其变化频率约为1HZ,根据信号采样定理,可将矩形单元格的边长a定为0.25mm。
本发明还可进一步包括用于对图像信号进行处理的图像采集与处理***,该***由摄像机、微机和图形处理及控制软件构成,在本发明的一个实施例中,摄像机采用CCD黑白摄像机,输出信号为CCIR制式。连续摄像速度为25帧/秒。图像采集卡为Matrox公司的Pulsar黑白图像采集卡,其采集速度为8位A/D,45MHZ;10位A/D,30MHZ,采集最大尺寸为2048×2048×8位。采用图形处理软件可对采集的图像数据进行处理并用来控制整个操作过程,最后重构被测液面的三维形貌并在显示器上加于显示。
图5为本实施例的图形处理软件流程图,为重构被测表面的形貌,本发明采取了一种同步解调的技术求出表面各点的倾斜角,进而来求出表面形貌。具体处理步骤如下:
1、对图像采集卡作初始化处理;
2、用图像采集卡从CCD摄像机获取光学图像信息:
3、用同步信号Isync(x,y)=cosω(x,y)分别对x,y两方向的光栅投影图像作同步解调以获取两方向的表面倾斜角θ(x,y);
4、求出表面x,y两方向的斜率;
5、由该斜率信息作积分,求得表面三维形貌。

Claims (7)

1、一种测量物体表面形貌的纹影仪,包括光源部件、光学成像部件,所说的光源部件由激光器及设置在其输出光束上的小孔滤波器和扩束透镜组成,所说的光学成像部件由将经扩束透镜扩束后的光束进行分束的分光镜、接收该分光镜反射的光束的反射镜、接收该反射镜反射的光束的准直透镜,以及接收该分光镜透射的光束的空间滤波器和CCD摄像机组成,其特征在于,所说的空间滤波器为一个二维朗奇光栅。
2、如权利要求1所述的纹影仪,其特征在于,所说的激光器为He-Ne激光器或半导体激光器。
3、如权利要求1所述的纹影仪,其特征在于,所说的准直透镜采用消像差非球面准直镜。
4、如权利要求1所述的纹影仪,其特征在于,所说的二维朗奇光栅被置于焦平面上。
5、如权利要求1所述的纹影仪,其特征在于,所说的二维朗奇光栅被置于焦平面之前的适当位置上。
6、如权利要求1所述的纹影仪,其特征在于,所说的CCD摄像机为CCD而阵摄像机。
7、如权利要求1、2、3、4、5或6所述的纹影仪,其特征在于,还包括用于对图像信号进行处理的图像采集与处理***,该***由用于采集所说的CCD摄像机获取被测表面形貌的信息的高速图像采集卡和对该图像数据进行处理的微机所构成,所说的高速图像采集卡与该微机的接口相连,所说的微机储存有图形处理及控制软件。
CN 00107371 2000-05-12 2000-05-12 测量物体表面形貌的纹影仪 Expired - Fee Related CN1114090C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 00107371 CN1114090C (zh) 2000-05-12 2000-05-12 测量物体表面形貌的纹影仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 00107371 CN1114090C (zh) 2000-05-12 2000-05-12 测量物体表面形貌的纹影仪

Publications (2)

Publication Number Publication Date
CN1271086A true CN1271086A (zh) 2000-10-25
CN1114090C CN1114090C (zh) 2003-07-09

Family

ID=4578654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 00107371 Expired - Fee Related CN1114090C (zh) 2000-05-12 2000-05-12 测量物体表面形貌的纹影仪

Country Status (1)

Country Link
CN (1) CN1114090C (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101881607A (zh) * 2010-06-10 2010-11-10 上海理工大学 一种检测平面误差***
CN102155925A (zh) * 2011-01-07 2011-08-17 中国科学院上海光学精密机械研究所 基于一维达曼光栅的三维表面形貌测量装置
CN102288392A (zh) * 2011-07-29 2011-12-21 温州医学院 一种基于二维朗奇光栅的自由曲面眼镜片光焦度测量装置
CN102636130A (zh) * 2012-04-19 2012-08-15 中国科学院光电技术研究所 一种大动态范围测量非球面光学元件表面轮廓装置
CN102967274A (zh) * 2012-11-14 2013-03-13 广东汉唐量子光电科技有限公司 一种测量物体表面形貌的方法
CN103424088A (zh) * 2013-08-12 2013-12-04 韵升控股集团有限公司 一种倒角测试方法及倒角测量仪
CN103913288A (zh) * 2014-03-14 2014-07-09 中国科学院力学研究所 彩虹纹影测量成像***及方法
CN104111120A (zh) * 2014-07-25 2014-10-22 中国科学院上海光学精密机械研究所 基于朗奇剪切干涉仪的相位提取方法
CN104482882A (zh) * 2014-12-22 2015-04-01 河海大学 一种漫顶溃坝试验中的溃口测量仪及测量方法
CN105044117A (zh) * 2015-08-14 2015-11-11 华侨大学 用于平板表面质量的检测方法与***
CN105067634A (zh) * 2015-08-14 2015-11-18 华侨大学 一种用于平板表面质量的检测方法与***
CN105890544A (zh) * 2014-12-10 2016-08-24 青岛理工大学 水下静止及高速移动目标三维成像方法及成像***
CN107430194A (zh) * 2015-01-30 2017-12-01 阿德科尔公司 光学三维扫描仪及其使用方法
CN108489414A (zh) * 2018-05-07 2018-09-04 中国矿业大学 一种微纳尺度记忆合金双程记忆应变测试装置及测试方法
CN108871299A (zh) * 2018-06-28 2018-11-23 中国船舶重工集团公司第七0七研究所 一种铅垂基准反射镜装置
CN111189408A (zh) * 2020-01-08 2020-05-22 西安工业大学 一种投影莫尔齿轮齿面三维测量装置和测量方法
CN116678583A (zh) * 2023-01-12 2023-09-01 中国空气动力研究与发展中心设备设计与测试技术研究所 一种基于相位调制的纹影***及其调节方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI485361B (zh) * 2013-09-11 2015-05-21 Univ Nat Taiwan 三維形貌輪廓量測裝置及其方法

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101881607A (zh) * 2010-06-10 2010-11-10 上海理工大学 一种检测平面误差***
CN102155925A (zh) * 2011-01-07 2011-08-17 中国科学院上海光学精密机械研究所 基于一维达曼光栅的三维表面形貌测量装置
CN102288392A (zh) * 2011-07-29 2011-12-21 温州医学院 一种基于二维朗奇光栅的自由曲面眼镜片光焦度测量装置
CN102636130B (zh) * 2012-04-19 2014-07-09 中国科学院光电技术研究所 一种大动态范围测量非球面光学元件表面轮廓装置
CN102636130A (zh) * 2012-04-19 2012-08-15 中国科学院光电技术研究所 一种大动态范围测量非球面光学元件表面轮廓装置
CN102967274A (zh) * 2012-11-14 2013-03-13 广东汉唐量子光电科技有限公司 一种测量物体表面形貌的方法
CN103424088A (zh) * 2013-08-12 2013-12-04 韵升控股集团有限公司 一种倒角测试方法及倒角测量仪
CN103424088B (zh) * 2013-08-12 2016-01-13 韵升控股集团有限公司 一种倒角测量仪
CN103913288A (zh) * 2014-03-14 2014-07-09 中国科学院力学研究所 彩虹纹影测量成像***及方法
CN104111120B (zh) * 2014-07-25 2017-05-31 中国科学院上海光学精密机械研究所 基于朗奇剪切干涉仪的相位提取方法
CN104111120A (zh) * 2014-07-25 2014-10-22 中国科学院上海光学精密机械研究所 基于朗奇剪切干涉仪的相位提取方法
CN105890544B (zh) * 2014-12-10 2019-02-05 青岛理工大学 水下静止及高速移动目标三维成像方法及成像***
CN105890544A (zh) * 2014-12-10 2016-08-24 青岛理工大学 水下静止及高速移动目标三维成像方法及成像***
CN104482882A (zh) * 2014-12-22 2015-04-01 河海大学 一种漫顶溃坝试验中的溃口测量仪及测量方法
CN107430194A (zh) * 2015-01-30 2017-12-01 阿德科尔公司 光学三维扫描仪及其使用方法
CN105067634A (zh) * 2015-08-14 2015-11-18 华侨大学 一种用于平板表面质量的检测方法与***
CN105044117A (zh) * 2015-08-14 2015-11-11 华侨大学 用于平板表面质量的检测方法与***
CN108489414A (zh) * 2018-05-07 2018-09-04 中国矿业大学 一种微纳尺度记忆合金双程记忆应变测试装置及测试方法
CN108871299A (zh) * 2018-06-28 2018-11-23 中国船舶重工集团公司第七0七研究所 一种铅垂基准反射镜装置
CN108871299B (zh) * 2018-06-28 2020-11-10 中国船舶重工集团公司第七0七研究所 一种铅垂基准反射镜装置
CN111189408A (zh) * 2020-01-08 2020-05-22 西安工业大学 一种投影莫尔齿轮齿面三维测量装置和测量方法
CN116678583A (zh) * 2023-01-12 2023-09-01 中国空气动力研究与发展中心设备设计与测试技术研究所 一种基于相位调制的纹影***及其调节方法
CN116678583B (zh) * 2023-01-12 2024-04-05 中国空气动力研究与发展中心设备设计与测试技术研究所 一种基于相位调制的纹影***及其调节方法

Also Published As

Publication number Publication date
CN1114090C (zh) 2003-07-09

Similar Documents

Publication Publication Date Title
CN1114090C (zh) 测量物体表面形貌的纹影仪
EP1739471B1 (en) Dual technology (confocal and intererometric) optical profilometer for the inspection and three-dimensional measurement of surfaces
US4340306A (en) Optical system for surface topography measurement
CN1069401C (zh) 利用大的等效波长测量物体表面轮廓的方法及***
CN101520306B (zh) 基于空间载波的干涉共焦测量装置与方法
US7477362B2 (en) Moiré interferometric strain sensor
EP2602583B1 (en) Low coherence interferometric system for phase stepping shearography combined with 3D profilometry
JPS60123704A (ja) 光学干渉計システムとその使用方法
US5218424A (en) Flying height and topography measuring interferometer
CN104596989A (zh) 基于干涉条纹图处理测量透明介质折射率分布的方法
CN111121644A (zh) 一种基于涡旋光与球面波干涉的微位移测量方法及装置
KR20050044901A (ko) 그레이징 입사 간섭계용 줄무늬 패턴 판별기
Sánchez et al. LED source interferometer for microscopic fringe projection profilometry using a Gates’ interferometer configuration
CN111397634B (zh) 星敏感器固定端面热变形的高分辨干涉检测装置及方法
US20050279954A1 (en) Three-dimensional imaging resonator and method therefor
CN1119646C (zh) 测量物体表面形貌的朗奇光栅纹影仪装置
US7206079B2 (en) Apparatus and process for characterizing samples
JPH05502731A (ja) 表面上に印刷されたまたは取り付けられた格子を使用するモアレ距離測定方法及び装置
US6614534B1 (en) Method and apparatus for combined measurement of surface non-uniformity index of refraction variation and thickness variation
JPH10232204A (ja) 屈折率測定装置
JP3540004B2 (ja) 斜入射干渉計
Joenathan et al. Phase measurement by differentiating interferometric fringes
Tiziani et al. Spectral and temporal phase evaluation for interferometry and speckle applications
EP0006945A1 (en) Method and apparatus for sensing in-plane deformation of a surface
US8526007B2 (en) Apparatus and method for measuring displacement

Legal Events

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