CN104898274B - 应用于光学***的自由曲面的设计方法 - Google Patents
应用于光学***的自由曲面的设计方法 Download PDFInfo
- Publication number
- CN104898274B CN104898274B CN201410077829.1A CN201410077829A CN104898274B CN 104898274 B CN104898274 B CN 104898274B CN 201410077829 A CN201410077829 A CN 201410077829A CN 104898274 B CN104898274 B CN 104898274B
- Authority
- CN
- China
- Prior art keywords
- characteristic
- point
- free
- characteristic data
- form surface
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000003287 optical effect Effects 0.000 title claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 17
- 238000003384 imaging method Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Graphics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Optics & Photonics (AREA)
- Lenses (AREA)
- Architecture (AREA)
- Image Generation (AREA)
Abstract
本发明涉及一种应用于光学***的自由曲面的设计方法,包括以下步骤:建立一曲面初始结构;选取K条特征光线,使所述特征光线通过该曲面初始结构后,分别交一像面于理想像点Ii(i=1,2…K);根据物像关系及斯涅尔定律逐点求解每条特征光线与待求的自由曲面的多个交点,进而得到特征数据点Pi;以及将该多个特征数据点Pi进行曲面拟合,得到所述自由曲面的方程式。
Description
技术领域
本发明涉及自由曲面设计领域,特别涉及一种基于逐点求解与构建的应用于光学***的自由曲面的设计方法。
背景技术
自由曲面是指无法用球面或非球面系数来表示的非传统曲面,通常是非回转对称的,结构灵活,变量较多,为光学设计提供了更多的自由度,可以大大降低光学***的像差,减小***的体积、重量与镜片数量,可以满足现代成像***的需要,有着广阔的发展应用前景。成像光学***要实现视场大小与孔径大小一定的成像,需要在成像***设计中控制不同视场不同孔径位置的光线。由于自由曲面有非对称面并提供了更多的设计自由度,他们常被用在离轴非对称***中。
现有的自由曲面设计方法主要针对曲面的二维轮廓进行求解,且只适用于单一视场***的求解,对于多视场,且具有一定大小孔径的***设计起来非常困难。
发明内容
综上所述,确有必要提供一种能够适用于多视场、且具有一定大小孔径的***,并可以针对三维的自由曲面进行求解的设计方法。
一种应用于光学***的自由曲面的设计方法,包括以下步骤:建立一曲面初始结构;选取K条特征光线,使所述特征光线通过该曲面初始结构后,分别交一像面于理想像点Ii(i=1,2…K);根据物像关系及斯涅尔定律逐点求解每条特征光线与待求的自由曲面的多个交点,进而得到特征数据点Pi;以及将该多个特征数据点Pi进行曲面拟合,得到所述自由曲面的方程式。
相较于现有技术,本发明提供的应用于光学***的自由曲面的设计方法,利用特征光线在待求的自由曲面与前后曲面之间形成的光学关系,通过对自由曲面进行逐点求解与构建的方法,得到满足一定的成像要求的三维自由曲面,方法简单,适用于各种离轴非对称***以及折、反射面;另外,本方法还可以通过对多视场与不同孔径位置的特征光线的控制,适用于多视场且具有一定孔径大小的成像***设计,并且***视场数量不受限,具有广阔的应用空间。
附图说明
图1为本发明实施例提供的应用于光学***的自由曲面设计方法流程图。
图2为本发明实施例提供的应用于光学***的自由曲面设计方法中每个视场的特征光线选择方法示意图。
图3为本发明提供的求解特征数据点时特征光线起点与终点示意图。
图4为为采用本发明实施例提供的应用于光学***的自由曲面设计方法得到的自由曲面示意图。
具体实施方式
下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
请参阅图1,本发明实施例一提供一种自由曲面的设计方法,该自由曲面的设计方法包括如下步骤:
步骤S1,建立一曲面初始结构;
步骤S2,选取K条特征光线,使所述特征光线通过该曲面初始结构后,分别交一像面于理想像点Ii(i=1,2…K);
步骤S3,根据物像关系及斯涅尔定律逐点求解每条特征光线与待求的自由曲面的多个交点,进而得到特征数据点Pi;以及
步骤S4,将该多个特征数据点Pi进行曲面拟合,得到所述自由曲面的方程式。
在步骤S1中,所述曲面初始结构可以为平面或曲面,且所述曲面初始结构的具***置可以根据光学***的实际需要进行选择,只要使得出射的光线与像面的交点与理想像点接近。本实施例中,所述曲面初始结构为一平面。
在步骤S2中,所述特征曲线的选取可通过以下方法进行:
可以根据需求选取M个视场,并将每个视场的孔径分成N等份,并从每一等份中选取不同孔径位置的P条特征光线,这样一共选取了K=M×N×P条对应不同视场不同孔径位置的特征光线。所述孔径可以为圆形、长方形、正方形、椭圆形或其他规则或不规则的形状。
请参阅图2,优选的,所述视场孔径为圆形孔径,将每个视场的圆形孔径等分成N个角度,间隔为φ,因此有N=2π/φ,沿着每个角度的半径方向取P个不同的孔径位置,那么一共取K=M×N×P条对应不同视场不同孔径位置的特征光线。本实施例中,共选取6个离轴视场进行设计,且将每个视场的圆形孔径等分成14个角度,沿着每个角度的半径方向取7个不同的孔径位置,因此共选取了588条对应不同视场不同孔径位置的特征光线。
在步骤S3中,请参阅图3,将特征光线与待求的自由曲面Ω的交点定义为特征数据点Pi。为了得到一个待求的自由曲面Ω上的所有特征数据点Pi,可以借助特征光线Ri与待求的自由曲面Ω的前一个曲面Ω'及后一个曲面Ω''的交点。将特征光线Ri与前一个曲面Ω'的交点定义为该特征光线的起点Si,特征光线Ri与后一个曲面Ω''的交点定义为该特征光线Ri的终点Ei。当待设计的自由曲面与特征光线Ri确定后,该特征光线Ri的起点Si是确定的,且易于通过光线追迹得到,特征光线Ri的终点Ei可通过物像关系求解。在理想状态下,特征光线Ri从Ω'上的Si射出后,经过Ω上的Pi,交于Ω''上的Ei,并最终交像面于其理想像点Ii。如果Ω''是像面,特征光线Ri的终点Ei就是其理想像点Ii;如果在待求的自由曲面Ω和像面之间还有其他像面,特征光线Ri的终点Ei是Ω''上使从Pi到Ii光程为最小值的点。
所述特征数据点Pi(i=1,2…K)可以通过以下两种计算方法获得。
第一种计算方法包括以下步骤:
步骤S31,取定一第一条特征光线R1与所述曲面初始结构的第一交点为特征数据点P1;
步骤S32,在得到i(1≤i≤K−1)个特征数据点Pi后,根据斯涅尔定律的矢量形式求解该特征数据点Pi处的单位法向量;
步骤S33,过所述i(1≤i≤K−1)个特征数据点Pi分别做一第一切平面,得到i个第一切平面,该i个第一切平面与其余K−i条特征光线相交得到i×(K−i)个第二交点,从该i×(K−i)个第二交点中选取出与所述i(1≤i≤K−1)个特征数据点Pi距离最近的一个第二交点,作为所述待求的自由曲面的下一个特征数据点Pi+1;
步骤S34,重复步骤S32和S33,直到计算得到所有特征数据点Pi(i=1,2…K),通过曲面拟合可以得到所述自由曲面的的方程式。
步骤S32中,每个特征数据点Pi处的单位法向量可以根据斯涅尔(Snell)定律的矢量形式求解。当待求的自由曲面Ω为折射面时,则每个特征数据点Pi (i=1,2…K)处的单位法向量满足:
(1)
其中,、 分别是沿着光线入射与出射方向的单位矢量,n, n'分别为待求的自由曲面Ω前后两种介质的折射率。
类似的,当待求的自由曲面Ω为反射面时,则每个特征数据点Pi (i=1,2…K)处的单位法向量满足:
(2)
由于,所述特征数据点Pi(i=1,2…K) 处的单位法向量与所述特征数据点Pi(i=1,2…K) 处的切平面垂直。故,可以得到特征数据点Pi(i=1,2…K) 处的切平面。
所述第一种算法的计算复杂度为,当设计中采用的特征光线的数量较多时,该方法需要较长的计算时间。
第二种计算方法包括以下步骤:
步骤S'31,取定一第一条特征光线R1与所述曲面初始结构的第一交点为特征数据点P1;
步骤S'32,在得到第i(1≤i≤K−1)个特征数据点Pi后,根据斯涅尔定律的矢量形式求解第i个特征数据点Pi处的单位法向量;
步骤S'33,仅过所述第i(1≤i≤K−1)个特征数据点Pi做一第一切平面并与其余K−i条特征光线相交,得到K−i个第二交点,从该K−i个第二交点中选取出与所述第i个特征数据点Pi距离最短的第二交点Qi+1,并将其对应的特征光线及与所述第i个特征数据点Pi的最短距离分别定义为Ri+1和d;
步骤S'34,过特征数据点Pi(1≤i≤K−1)之前已求得的i−1个特征数据点分别做一第二切平面,得到i−1个第二切平面,该i−1个第二切平面与所述特征光线Ri+1相交得到i−1个第三交点,在每一第二切平面上每一第三交点与其所对应的特征数据点Pi形成一交点对,在所述交点对中,选出交点对中距离最短的一对,并将距离最短的交点对的第三交点和最短距离分别定义为Q(i+1)'和di';
步骤S'35,比较di与di',如果di≤di',则把Qi+1取为下一个特征数据点Pi+1,反之,则把Q(i+1)'取为下一个特征数据点Pi+1;以及
步骤S'36,重复步骤S'32到S'35,直到计算得到所有特征数据点Pi,通过曲面拟合可以得到所述自由曲面的的方程式。
步骤S'32中,每个特征数据点Pi处的单位法向量的计算方法与所述第一种算法中相同。
所述第二种算法的计算复杂度为,当设计中个采用的特征光线的数量较多时,所述第二种算法明显比第一种算法的计算复杂度小。优选的,采用所述第二种算法逐点求解特征数据点Pi(i=1,2…K)。
可以理解,本实施例得到的自由曲面还可以作为供进一步优化的初始结构。
请参阅图4,本实施例中采用所述第二种算法得到所述588条特征光线与待求的自由曲面Ω的588个交点即特征数据点P i(i=1,2…588),将该588个特征数据点Pi(i=1,2…588)进行自由曲面拟合,得到所述自由曲面。
本实施例提供的应用于光学***的自由曲面的设计方法,利用特征光线在待求的自由曲面与前后曲面之间形成的光学关系,通过对曲面进行逐点求解与构建的方法,得到满足一定的成像要求的三维自由曲面,方法简单,适用于各种离轴非对称***以及折射面及反射面;另外,本方法通过对多视场与不同孔径位置的特征光线的控制,适用于多视场且有一定的孔径大小的成像***设计,并且***视场数量不受限,具有广阔的应用空间。
另外,本领域技术人员还可在本发明精神内作其它变化,当然这些依据本发明精神所作的变化,都应包含在本发明所要求保护的范围内。
Claims (7)
1.一种应用于光学***的自由曲面的设计方法,包括以下步骤:
步骤S1,建立一曲面初始结构;
步骤S2,选取对应不同视场不同孔径位置K条特征光线,使所述特征光线通过该曲面初始结构后,分别交一像面于理想像点Ii(i=1,2…K),所述选取K条特征光线包括如下步骤:选取M个视场,将每个视场的孔径分成N等份,每一等份中选取不同孔径位置的P条特征光线;
步骤S3,根据物像关系及斯涅尔定律逐点求解每条特征光线与待求的自由曲面的多个交点,进而得到特征数据点Pi(i=1,2…K),该特征数据点Pi(i=1,2…K)的求解步骤包括:
步骤S′31,取定一第一条特征光线R1与所述曲面初始结构的第一交点为特征数据点P1;
步骤S′32,在得到第i(1≤i≤K-1)个特征数据点Pi后,根据斯涅尔定律的矢量形式求解第i个特征数据点Pi处的单位法向量
步骤S′33,仅过所述第i(1≤i≤K-1)个特征数据点Pi做一第一切平面并与其余K-i条特征光线相交,得到K-i个第二交点,从该K-i个第二交点中选取出与所述第i个特征数据点Pi距离最短的第二交点Qi+1,并将该第二交点Qi+1对应的特征光线定义为Ri+1,将该第二交点Qi+1与所述第i个特征数据点Pi的最短距离分别定义为d;
步骤S′34,过特征数据点Pi(1≤i≤K-1)之前已求得的i-1个特征数据点分别做一第二切平面,得到i-1个第二切平面,该i-1个第二切平面与所述特征光线Ri+1相交得到i-1个第三交点,在每一第二切平面上每一第三交点与该第三交点所对应的特征数据点Pi形成一交点对,在所述交点对中,选出交点对中距离最短的一对,并将距离最短的交点对的第三交点和最短距离分别定义为Q(i+1)′和di′;
步骤S′35,比较di与di′,如果di≤di′,则把Qi+1取为下一个特征数据点Pi+1,反之,则把Q(i+1)′取为下一个特征数据点Pi+1;以及
步骤S′36,重复步骤S′32到S′35,直到计算得到所有特征数据点Pi(i=1,2…K);以及
步骤S4,将该多个特征数据点Pi(i=1,2…K)进行曲面拟合,得到所述自由曲面的方程式,该自由曲面为三维的。
2.如权利要求1所述的应用于光学***的自由曲面的设计方法,其特征在于,所述曲面初始结构为平面或曲面。
3.如权利要求1所述的应用于光学***的自由曲面的设计方法,其特征在于,所述每个视场的孔径均为圆形。
4.如权利要求1所述的应用于光学***的自由曲面的设计方法,其特征在于,将每个视场的圆形孔径等分成N个角度。
5.如权利要求4所述的应用于光学***的自由曲面的设计方法,其特征在于,沿着每个角度的半径方向取P个不同的孔径位置。
6.如权利要求1所述的应用于光学***的自由曲面的设计方法,其特征在于,将与自由曲面相邻的前一个曲面与所述特征光线的交点定义为起点Si(i=1,2…K),将与自由曲面相邻的后一个曲面与所述特征光线的交点定义为终点Ei(i=1,2…K),所述待求的自由曲面为折射面,则特征数据点Pi处的单位法向量为:
其中,分别为沿着特征光线入射与出射方向的单位矢量,n,n′分别为待求的自由曲面两侧两种介质的折射率。
7.如权利要求1所述的应用于光学***的自由曲面的设计方法,其特征在于,将与自由曲面相邻的前一个曲面与所述特征光线的交点定义为起点Si(i=1,2…K),将与自由曲面相邻的后一个曲面与所述特征光线的交点定义为终点Ei(i=1,2…K),所述待求的自由曲面为反射面,则特征数据点Pi处的单位法向量为:
其中,分别为沿着特征光线入射与出射方向的单位矢量。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410077829.1A CN104898274B (zh) | 2014-03-05 | 2014-03-05 | 应用于光学***的自由曲面的设计方法 |
TW103111875A TWI499798B (zh) | 2014-03-05 | 2014-03-28 | 應用於光學系統的自由曲面的設計方法 |
US14/616,457 US10387580B2 (en) | 2014-03-05 | 2015-02-06 | Method of forming a freeform surface shaped element for an optical system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410077829.1A CN104898274B (zh) | 2014-03-05 | 2014-03-05 | 应用于光学***的自由曲面的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898274A CN104898274A (zh) | 2015-09-09 |
CN104898274B true CN104898274B (zh) | 2018-01-30 |
Family
ID=54017592
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410077829.1A Active CN104898274B (zh) | 2014-03-05 | 2014-03-05 | 应用于光学***的自由曲面的设计方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10387580B2 (zh) |
CN (1) | CN104898274B (zh) |
TW (1) | TWI499798B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898275B (zh) * | 2014-03-05 | 2017-07-07 | 清华大学 | 自由曲面离轴三反光学***的设计方法 |
CN108152948B (zh) * | 2016-12-05 | 2020-02-07 | 清华大学 | 离轴非球面光学***的设计方法 |
CN108227184B (zh) * | 2016-12-22 | 2020-01-07 | 清华大学 | 自由曲面成像光学***的设计方法 |
CN108345106B (zh) * | 2017-01-24 | 2020-04-28 | 清华大学 | 混合表面光学***的设计方法 |
CN109387938B (zh) * | 2017-08-14 | 2020-08-11 | 清华大学 | 光学***自由曲面面形公差分布的分析方法 |
CN109239913B (zh) * | 2018-08-14 | 2020-11-10 | 北京理工大学 | 基于平板相位元件的成像光学***的设计方法 |
CN109239917B (zh) * | 2018-10-16 | 2021-07-13 | 北京航天易联科技发展有限公司 | 一种适用于太赫兹安检的透镜轮廓面获取方法和装置 |
CN111487766B (zh) * | 2019-01-25 | 2021-06-25 | 清华大学 | 自由曲面离轴三反成像***的设计方法 |
CN109752843B (zh) * | 2019-03-06 | 2021-12-31 | 西安自由曲面光电科技有限公司 | 自由曲面大视场低畸变汽车后视镜的实现方法 |
CN110018566A (zh) * | 2019-04-25 | 2019-07-16 | 钟祥博谦信息科技有限公司 | 自由曲面光学***设计的方法、设备、***及存储介质 |
CN112710094B (zh) * | 2021-02-01 | 2022-06-03 | 上海晶电新能源有限公司 | 一种二次反射***及具有该***的太阳能聚光集热*** |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US93001A (en) * | 1869-07-27 | Improved pitcher for cooling liquids | ||
US7460985B2 (en) * | 2003-07-28 | 2008-12-02 | Light Prescriptions Innovators, Llc | Three-dimensional simultaneous multiple-surface method and free-form illumination-optics designed therefrom |
TWI279589B (en) * | 2006-03-24 | 2007-04-21 | Hon Hai Prec Ind Co Ltd | System and method for designing free surface of a reflector |
WO2008066742A1 (en) * | 2006-11-22 | 2008-06-05 | Geng Z Jason | Wide field-of-view reflector and method of designing and making same |
EP2316046B1 (en) * | 2008-08-08 | 2017-10-04 | LPI-Europe, S.L. | Imaging optics designed by the simultaneous multiple surface method |
TWI432790B (zh) * | 2010-04-29 | 2014-04-01 | Univ Nat Changhua Education | The collimating lens structure with free - form surface and its design method |
TWI476096B (zh) * | 2011-04-19 | 2015-03-11 | Univ Asia | 鏡片設計系統及其方法 |
CN102928983B (zh) * | 2012-10-23 | 2014-08-27 | 浙江大学 | 一种用于准直光整形的自由曲面光学元件的设计方法 |
CN103809277B (zh) * | 2012-11-06 | 2016-09-14 | 清华大学 | 离轴三反镜 |
CN103136431B (zh) * | 2013-03-26 | 2015-07-15 | 山东大学 | 一种自由曲面的保面积参数特性优化方法 |
CN103246053A (zh) * | 2013-04-09 | 2013-08-14 | 长春理工大学 | 使用自由曲面的宽幅离轴三反光学*** |
US10151921B2 (en) * | 2013-11-04 | 2018-12-11 | Synopsys, Inc. | Optical design using freeform tailoring |
-
2014
- 2014-03-05 CN CN201410077829.1A patent/CN104898274B/zh active Active
- 2014-03-28 TW TW103111875A patent/TWI499798B/zh active
-
2015
- 2015-02-06 US US14/616,457 patent/US10387580B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US20150254375A1 (en) | 2015-09-10 |
TW201534974A (zh) | 2015-09-16 |
US10387580B2 (en) | 2019-08-20 |
TWI499798B (zh) | 2015-09-11 |
CN104898274A (zh) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898274B (zh) | 应用于光学***的自由曲面的设计方法 | |
CN104898275B (zh) | 自由曲面离轴三反光学***的设计方法 | |
CN104570340B (zh) | 自由曲面成像***的设计方法 | |
TWI544230B (zh) | 自由曲面離軸三反成像系統的設計方法 | |
CN105988212B (zh) | 三维自由曲面的设计方法 | |
TWI557436B (zh) | 自由曲面離軸光學系統的設計方法 | |
CN108227184B (zh) | 自由曲面成像光学***的设计方法 | |
TWI622813B (zh) | 離軸非球面光學系統的設計方法 | |
TWI557433B (zh) | 具有實出瞳的自由曲面離軸成像系統的設計方法 | |
TWI667450B (zh) | 光學系統自由曲面面形公差分佈的分析方法 | |
TWI622800B (zh) | 離軸非球面三反光學系統 | |
CN111487766B (zh) | 自由曲面离轴三反成像***的设计方法 | |
TW201827886A (zh) | 混合表面離軸三反光學系統 | |
CN112305737B (zh) | 自由曲面反射式红外成像***的设计方法 | |
WO2006120750A1 (ja) | 実数値面番号を用いた光学系設計方法 |
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 |