CN104714222B - 激光雷达***回波能量的计算方法 - Google Patents
激光雷达***回波能量的计算方法 Download PDFInfo
- Publication number
- CN104714222B CN104714222B CN201510165389.XA CN201510165389A CN104714222B CN 104714222 B CN104714222 B CN 104714222B CN 201510165389 A CN201510165389 A CN 201510165389A CN 104714222 B CN104714222 B CN 104714222B
- Authority
- CN
- China
- Prior art keywords
- radar system
- computational methods
- laser beam
- radius
- laser
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明涉及一种激光雷达***回波能量的计算模型。本模型基于光线追迹方法,可用于任意激光强度分布、任意孔径遮挡的共轴或双轴激光雷达***,获取探测器于任意像面位置处所接收的激光回波能量。本发明提出的激光回波能量计算模型具有普适性强、精度高、速度快等特点。
Description
技术领域
本发明涉及激光雷达技术领域,尤其涉及一种激光雷达***回波能量计算方法。
背景技术
激光雷达是一种具有极高时空分辨率和测量精度的遥感设备,广泛应用于无人导航车、三维城市建模、地形勘测、大气探测等技术领域。精确计算出预探测区域内的激光回波能量,对激光雷达***的整机设计和性能评估具有指导意义。
1978年,J.Harms基于光学成像方法提出了典型高斯能量分布的激光回波能量计算模型,并用于共轴或双轴存有中心遮挡的激光雷达***;1994年,Jin Wang和JuhaKostamovaara针对半导体激光器的能量分布,根据辐亮度守恒的原则建立了激光回波能量计算模型;2005年,Kamil Stelmaszczyk等人又将光学成像方法扩展到更复杂的牛顿望远镜***中。
发明内容
针对背景技术存在的问题,本发明提供一种激光雷达***回波能量计算方法。本发明采用如下技术方案:
一种激光雷达***中回波能量计算方法,包括以下步骤:
步骤1,获取落入接收光学***视场内的激光强度二维分布G(i,j)(X,Y,Z);
步骤2,计算采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z);
步骤3,计算采样点G(i,j)(X,Y,Z)回波主光线与接收光轴的夹角γ(i,j)(Z);
步骤4,计算采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量P(i,j)(Z);
步骤5,计算激光回波总能量Pd(Z)。
所述步骤1中,落入接收光学***视场内的激光强度二维分布G(i,j)(X,Y,Z)为:
G(i,j)(X,Y,Z)=P0G(X,Y,Z)A(i,j)(X,Y)
其中G(X,Y,Z)为任意激光强度二维分布,(Xd,Yd,Rd)为接收光学***的入瞳圆心坐标和半径,R为接收视场半径,接收半视场角,P0为入射激光峰值功率;通过定义n×n记录矩阵A(i,j)(X,Y)判断G(X,Y,Z)是否落入接收视场半径R内:若是,则A(i,j)(X,Y)为1;否则,A(i,j)(X,Y)为0。
所述步骤2中,采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z)的获取方法如下:
若探测器的像面位置为L′,接收光学***的入窗位置L和半径rw为:
其中f为接收光学***的焦距,rd为探测器的半径;
入窗经采样点G(i,j)(X,Y,Z)在接收光学***入瞳处的投影轮廓方程φ(i,j)(Xw,Yw,Rw)为:
φ(i,j)(Xw,Yw,Rw)=(Xw-Xd)2+(Yw-Yd)2-Rw 2=0
其中(Xw,Yw,Rw)为入窗投影的圆心坐标和半径;
接收光学***的入瞳有效面积S(i,j)(Z)为:
S(i,j)(Z)为去掉中心和边缘圆形遮挡部分后,入窗投影与入瞳的重叠面积;其中(Xd,Yd,R′d)为中心遮挡部分的圆心坐标和半径,(Xs,Ys,Rs)为边缘遮挡部分的圆心坐标和半径;通过定义m×m记录矩阵B(i,j)(X,Y)辅助计算,若B(i,j)(X,Y)在入瞳有效面积S(i,j)(Z)内,则B(i,j)(X,Y)为1;否则,B(i,j)(X,Y)为0;
采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z)为:
所述步骤3中,采样点G(i,j)(X,Y,Z)回波主光线与接收光轴的夹角γ(i,j)(Z)为:
所述步骤4中,采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量P(i,j)(Z)为:
其中,τ为单程大气透射率,ε为目标反射率,η1为发射光学***的透射率,η2为接收光学***的透射率,θ为发射光轴与目标面法线的夹角。
所述步骤5中,激光回波总能量Pd(Z)为:
Pd(Z)=∑∑P(i,j)(Z)。
所述的任意入射激光强度二维分布G(X,Y,Z)可通过解析式来描述,或通过CCD相机等设备实测真实的激光光源来获取。
所述的入瞳有效面积S(i,j)(Z)的计算方法适用于无遮挡,或多处遮挡且遮挡形状不规则的接收光学***。
通过改变探测器的像面位置L′,获取探测器所接收的激光回波能量Pd(Z)。
与现有技术相比,本发明提出了一种更准确、更具普适性的激光回波能量计算模型。该模型基于光线追迹方法,可用于任意激光强度分布的光源,既可以通过解析式来描述也可以通过实测获取;该模型可用于任意遮挡位置、任意遮挡形状的接收光学***;该模型还可以获取探测器于任意像面位置处所接收的激光回波能量,利于通过对探测器位置进行调焦来获得更佳的能量响应。
激光雷达的研究在国内尚处于起步阶段,类似的研究工作鲜有报道,本发明弥补了该技术领域的空白。
附图说明
图1双轴激光雷达***架构图;
图2入射激光与接收视场的位置关系示意图;
图3激光回波能量算法原理图;
图4存在中心和边缘圆形遮挡的接收入瞳孔径示意图;
图5存在边缘矩形遮挡的接收入瞳孔径示意图;
图6实例1中高斯激光束强度分布;
图7实例1中激光回波能量随探测距离的变化趋势图;
图8实例2中平顶激光束强度分布;
图9实例2中激光回波能量随探测距离的变化趋势图;
图10实例3中实测半导体激光器强度分布;
图11实例3中激光回波能量随探测距离的变化趋势图。
具体实施方式
现通过实施例,并结合附图,对本发明的技术方案做进一步的解释。
图1为双轴激光雷达***架构图(若共轴架构,则d为0)。为了精确计算出不同探测距离处的激光回波能量,需要考虑到三点重要因素:
1.在双轴***中,由于重叠因子的存在使得近、中距离处的入射激光束无法完全落入接收视场内,使得该探测距离区间内激光回波能量被部分利用;
2.探测器通常置于接收光学***焦平面处,对远距离处的目标进行探测,受探测器尺寸和接收光学***的焦距限制,导致近、中距离处的激光回波能量受离焦效应影响而损失;
3.较复杂的接收光学***一般都有遮挡存在,对激光回波能量造成一定的影响。
通过光线追迹的方法,首先需要对入射至目标面的激光强度进行采样,获取落入接收视场内的激光强度二维分布G(i,j)(X,Y,Z),并视为有效采样点;其次,计算每个有效采样点G(i,j)(X,Y,Z)反射至探测器上的激光回波能量P(i,j)(Z),并进行求和以获取总能量Pd(Z)。
图2为入射激光与接收视场的位置关系示意图,其中深色阴影部分表示落入接收视场内激光能量,其强度二维分布G(i,j)(X,Y,Z)为:
G(i,j)(X,Y,Z)=P0G(X,Y,Z)A(i,j)(X,Y)
其中G(X,Y,Z)为任意激光强度二维分布,(Xd,Yd,Rd)为接收光学***的入瞳圆心坐标和半径,R为接收视场半径,接收半视场角,P0为入射激光峰值功率。通过定义n×n记录矩阵A(i,j)(X,Y)判断G(X,Y,Z)是否落入接收视场半径R内:若是,则A(i,j)(X,Y)为1;否则,A(i,j)(X,Y)为0。
G(X,Y,Z)可以通过解析式来描述,以标准高斯光束为例,其强度二维分布为:
其中ω0为束腰半径,λ为波长,δ为激光束发散半角,d为发射光轴与接收光轴的距离(共轴***中d为0),Δυ为发射光轴与接收光轴之间的夹角,C0是常数且满足:
公式(3)使得入射激光束的总能量归一化成恒定值,不随探测距离发生变化。
为消除理想光源与真实光源的差异而引入的计算误差,G(X,Y,Z)还可以通过CCD相机等设备实测真实的激光光源来获取。
图3表示为激光回波能量算法原理图,图中任意复杂的接收光学***可等效于薄透镜,其孔径即为入瞳(光学像差忽略不计);探测器放置于透镜焦平面附近,其在物空间的所成实像为入窗。
有效采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量需要通过接收光学***的入瞳和入窗,如图3中的阴影包络部分,其表示采样点G(i,j)(X,Y,Z)激光回波光束的有效立体角ψ(i,j)(Z),即为入瞳与入窗分别对应于采样点G(i,j)(X,Y,Z)的立体角的交集。由投影关系可知,ψ(i,j)(Z)为入瞳有效面积S(i,j)(Z)与探测距离Z的平方之比,其中入瞳有效面积S(i,j)(Z)为除去孔径遮挡部分外,入窗经采样点G(i,j)(X,Y,Z)在入瞳处的投影与入瞳的重叠面积。
若探测器的像面位置为L′,探测器经接收光学***在物空间所成像即为入窗,则入窗位置L和半径rw为:
其中f为接收光学***的焦距,rd为探测器的半径。
入窗经采样点G(i,j)(X,Y,Z)在接收光学***入瞳处的投影轮廓方程φ(i,j)(Xw,Yw,Rw)为:
φ(i,j)(Xw,Yw,Rw)=(Xw-Xd)2+(Yw-Yd)2-Rw 2=0
其中(Xw,Yw,Rw)为入窗投影的圆心坐标和半径;
图4表示为存在遮挡的接收入瞳孔径示意图,以存在中心和边缘圆形遮挡为例,其中深色阴影部分为入瞳有效面积S(i,j)(Z):
其中(Xd,Yd,R′d)为中心遮挡部分的圆心坐标和半径,其中(Xs,Ys,Rs)为边缘遮挡部分的圆心坐标和半径;通过定义m×m记录矩阵B(i,j)(X,Y)辅助计算,若B(i,j)(X,Y)在入瞳有效面积S(i,j)(Z)内,则B(i,j)(X,Y)为1;否则,B(i,j)(X,Y)为0。
接收入瞳有效面积S(i,j)(Z)的计算方法同样可用于存在边缘矩形遮挡的情况,如图5所示,其入瞳有效面积S(i,j)(Z)为:
采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z)为:
受应用条件的约束,接收光学***结构形式各异。在超远距离探测用途中,接收光学***通常采用反射式望远镜结构,次镜的存在将导致中心遮挡;为了减少近距离处的探测盲区并有效控制***尺寸,双轴激光雷达的发射与接收孔径之间存在一定的重叠,将造成边缘遮挡。
若将探测目标视为朗伯面,点源向空间内规定方向上立体角内的辐射量与该方向与表面法线方向有关,因此在图2中需要考虑到采样点G(i,j)(X,Y,Z)回波主光线与接收光轴的夹角γ(i,j)(Z):
由公式(1)、(7)和(8)得,采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量P(i,j)(Z)为:
其中,τ为单程大气透射率,ε为目标反射率,η1为发射光学***的透射率,η2为接收光学***的透射率,θ为发射光轴与目标面法线的夹角。
对落入接收视场内全部有效采样点G(i,j)(X,Y,Z)的激光回波能量P(i,j)(Z)进行求和,得到探测距离Z处的激光回波总能量Pd(Z)为:
Pd(Z)=∑P(i,j)(Z) (11)
在超远距离探测用途中,往往将探测器置于接收光学***的焦平面处,针对无穷远的目标进行探测(入窗位置在无穷远);而在短距离探测用途中,需要通过调整探测器的像面位置,针对有限远的目标进行探测(入窗位置在有限距离处)。因此,通过改变式(4)中探测器的像面位置L′,获取探测器于任意像面位置处所接收的激光回波能量Pd(Z),有助于实现最佳的激光回波能量响应。
具体实施例
通过三个实例来说明本发明的可行性。假设单程大气透射率τ为0.98,目标反射率ε为0.1,发射光学***的透射率η1为0.9,接收光学***的透射率η2为0.85,发射光轴与目标面法线的夹角θ为0°。利用上述激光回波能量计算模型于表1中的3个实例,其光学参数为:
表1
实施例1:
图6为实例1中高斯激光束强度分布,图7为实例1中激光回波能量随探测距离的变化趋势图。
图7中1m处的激光回波能量为5.71×10-6W,150m处的激光回波能量为9.79×10- 8W,7m处激光回波能量达到最大值2.09×10-5W,1m至150m探测距离内激光回波能量的动态范围约213倍。
实施例2:
图8为实例2中平顶激光束强度分布,图9为实例2中激光回波能量随探测距离的变化趋势图。
图9中5m处的激光回波能量为1.5×10-5W,500m处的激光回波能量为1.39×10-7W,10m处激光回波能量达到最大值1.53×10-5W,5m至500m探测距离内激光回波能量的动态范围约110倍。
实施例3:
图10为实例3中实测半导体激光器强度分布,图11为实例3中激光回波能量随探测距离的变化趋势图。
图11中1m处的激光回波能量为5.29×10-6W,150m处的激光回波能量为9.8×10- 8W,7m处激光回波能量达到最大值6.5×10-6W,1m至150m探测距离内激光回波能量的动态范围约66.4倍。
通过3个实例说明了本发明提出的激光雷达回波能量计算模型可用于任意激光强度分布、任意孔径遮挡的共轴或双轴激光雷达***,获取探测器于任意像面位置处所接收的激光回波能量。需要强调的是,上述3个计算实例均建立在同一反射特性目标物、大气透射率、光学***透射率等条件下,在实际的应用中考虑到上述参数的变化性,激光回波能量将会产生变化。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (8)
1.一种激光雷达***中回波能量计算方法,其特征在于,包括以下步骤:
步骤1,获取落入接收光学***视场内的激光强度二维分布G(i,j)(X,Y,Z);
步骤2,计算采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z);
步骤3,计算采样点G(i,j)(X,Y,Z)回波主光线与接收光轴的夹角γ(i,j)(Z);
步骤4,计算采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量P(i,j)(Z);
步骤5,计算激光回波总能量Pd(Z);
所述步骤1中,落入接收光学***视场内的激光强度二维分布G(i,j)(X,Y,Z)为:
G(i,j)(X,Y,Z)=P0G(X,Y,Z)A(i,j)(X,Y)
其中,G(X,Y,Z)为任意激光强度二维分布,(Xd,Yd,Rd)为接收光学***的入瞳圆心坐标和半径,R为接收视场半径,接收半视场角,P0为入射激光峰值功率;通过定义n×n记录矩阵A(i,j)(X,Y)判断G(X,Y,Z)是否落入接收视场半径R内:若是,则A(i,j)(X,Y)为1;否则,A(i,j)(X,Y)为0。
2.根据权利要求1所述的激光雷达***回波能量计算方法,其特征在于,所述步骤2中,采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z)的获取方法如下:
若探测器的像面位置为L′,接收光学***的入窗位置L和半径rw为:
其中f为接收光学***的焦距,rd为探测器的半径;
入窗经采样点G(i,j)(X,Y,Z)在接收光学***入瞳处的投影轮廓方程φ(i,j)(Xw,Yw,Rw)为:
φ(i,j)(Xw,Yw,Rw)=(Xw-Xd)2+(Yw-Yd)2-Rw2=0
其中(Xw,Yw,Rw)为入窗投影的圆心坐标和半径;
接收光学***的入瞳有效面积S(i,j)(Z)为:
S(i,j)(Z)为去掉中心和边缘圆形遮挡部分后,入窗投影与入瞳的重叠面积,其中(Xd,Yd,R′d)为中心遮挡部分的圆心坐标和半径,(Xs,Ys,Rs)为边缘遮挡部分的圆心坐标和半径,通过定义m×m记录矩阵B(i,j)(X,Y)辅助计算,若B(i,j)(X,Y)在入瞳有效面积S(i,j)(Z)内,则B(i,j)(X,Y)为1;否则,B(i,j)(X,Y)为0,
采样点G(i,j)(X,Y,Z)回波光束的有效立体角ψ(i,j)(Z)为:
3.根据权利要求2所述的激光雷达***回波能量计算方法,其特征在于,所述步骤3中,采样点G(i,j)(X,Y,Z)回波主光线与接收光轴的夹角γ(i,j)(Z)为:
4.根据权利要求1所述的激光雷达***回波能量计算方法,其特征在于,所述步骤4中,采样点G(i,j)(X,Y,Z)反射至探测器的激光回波能量P(i,j)(Z)为:
其中,τ为单程大气透射率,ε为目标反射率,η1为发射光学***的透射率,η2为接收光学***的透射率,θ为发射光轴与目标面法线的夹角。
5.根据权利要求1所述的激光雷达***回波能量计算方法,其特征在于,所述步骤5中,激光回波总能量Pd(Z)为:
Pd(Z)=∑∑P(i,j)(Z)。
6.根据权利要求2所述的激光雷达***回波能量计算方法,其特征在于,所述的任意入射激光强度二维分布G(X,Y,Z)通过解析式来描述,或通过CCD相机实测真实的激光光源来获取。
7.根据权利要求3所述的激光雷达***回波能量计算方法,其特征在于,所述的入瞳有效面积S(i,j)(Z)的计算方法适用于无遮挡,或多处遮挡且遮挡形状不规则的接收光学***。
8.根据权利要求3所述的激光雷达***回波能量计算方法,其特征在于,通过改变探测器的像面位置L′,获取探测器所接收的激光回波能量Pd(Z)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510165389.XA CN104714222B (zh) | 2015-04-09 | 2015-04-09 | 激光雷达***回波能量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510165389.XA CN104714222B (zh) | 2015-04-09 | 2015-04-09 | 激光雷达***回波能量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104714222A CN104714222A (zh) | 2015-06-17 |
CN104714222B true CN104714222B (zh) | 2017-07-25 |
Family
ID=53413708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510165389.XA Expired - Fee Related CN104714222B (zh) | 2015-04-09 | 2015-04-09 | 激光雷达***回波能量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104714222B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108445471B (zh) * | 2018-03-26 | 2019-08-02 | 武汉大学 | 一种单光子激光雷达多探测器条件下的测距精度评估方法 |
CN109541619B (zh) * | 2018-12-13 | 2020-03-10 | 武汉大学 | 散斑相干性对单光子激光雷达测距精度影响的评估方法 |
CN210212200U (zh) * | 2019-01-03 | 2020-03-31 | 华域视觉科技(上海)有限公司 | 三维探测照明***及汽车 |
CN111308467A (zh) * | 2020-03-10 | 2020-06-19 | 宁波飞芯电子科技有限公司 | 探测方法及探测设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5936229A (en) * | 1996-04-02 | 1999-08-10 | Trw Inc. | Tracking means for distant ballistic missile targets comprising means for tracking largest radius of curvature |
CN102323592B (zh) * | 2011-08-12 | 2013-04-03 | 中国科学院光电技术研究所 | 一种目标回波信号的归一化方法 |
-
2015
- 2015-04-09 CN CN201510165389.XA patent/CN104714222B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104714222A (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104714222B (zh) | 激光雷达***回波能量的计算方法 | |
EP2211202B1 (en) | Apparatus and method for shaping a laser beam profile | |
US7375801B1 (en) | Video sensor with range measurement capability | |
CN102564611B (zh) | 大功率激光波前测量仪及波前测量方法 | |
US7397540B2 (en) | Phase diversity ranging sensor | |
CN104501956A (zh) | 一种超宽波段图谱关联探测装置与探测方法 | |
CN104655016A (zh) | 一种基于激光原向反射式光幕的弹丸着靶坐标测试方法 | |
CN107607195B (zh) | 一种基于复振幅实时获取的光束质量测量方法 | |
CN106053023A (zh) | 一种红外测量***自身杂散辐射的分析方法 | |
CN100535626C (zh) | 用光栅型波前曲率传感器测量焦距和等效f数的方法 | |
CN106443638A (zh) | 一种激光回波传输特性的分析方法、验证***及验证方法 | |
RU2639321C1 (ru) | Оптико-электронная система обнаружения объектов | |
CN104714221B (zh) | 激光雷达***回波能量动态范围的压缩方法 | |
US20210223022A1 (en) | System and method | |
CN109239917B (zh) | 一种适用于太赫兹安检的透镜轮廓面获取方法和装置 | |
Bobrovnikov et al. | Alignment technique and quality check of the primary mirror of the Siberian Lidar Station | |
Collins et al. | Design and simulation of a sensor for heliostat field closed loop control | |
Kim et al. | Effects of solar noise on the detection range performance of a laser spot tracker | |
Adel et al. | Design and implementation of a promising optical subsystem with a sky camera for laser warning systems | |
Agishev et al. | Imaging S-lidars enhancement by optimizing range-domain characteristics | |
CN109931885B (zh) | 一种精确测量光束入射角的装置 | |
Woo et al. | Time-of-flight measurement-based three-dimensional profiler system employing a lightweight Fresnel-type Risley prism scanner | |
Zhu et al. | A wide angle bistatic scanning lidar for navigation | |
Wojtanowski | Simulation of a cylindrical glass dome negative impact on a 360° field of view 2D laser scanner performance and a method for correction | |
Meraz et al. | Optical design of a deep water airborne bathymetric lidar |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170725 |