CN106053023A - 一种红外测量系统自身杂散辐射的分析方法 - Google Patents
一种红外测量系统自身杂散辐射的分析方法 Download PDFInfo
- Publication number
- CN106053023A CN106053023A CN201610390112.1A CN201610390112A CN106053023A CN 106053023 A CN106053023 A CN 106053023A CN 201610390112 A CN201610390112 A CN 201610390112A CN 106053023 A CN106053023 A CN 106053023A
- Authority
- CN
- China
- Prior art keywords
- stray radiation
- radiation
- calculate
- detector
- radiating 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
公开了一种红外测量系统自身杂散辐射的分析方法,包括:确定所述系统中的杂散辐射源,其中,所述杂散辐射源的总个数为N个,所述杂散辐射源包括:光学元件自身的热辐射、机械结构自身的热辐射;对于第i个杂散辐射源,计算其辐射面上辐射的、且能够被探测器接收的辐射功率φi;对φi的传输过程进行分析,计算其在探测器上的辐射照度Ei;计算N个杂散辐射源在探测器上的总辐射照度E。其中,N为大于1的整数,i=1,…N。本发明的方法能够对杂散辐射与温度的关系进行定量分析,从而可为不同温度下杂散辐射的抑制提供理论制导。进一步的,通过在不同温度下有效抑制杂散辐射,能够有效提高系统的测量精度。
Description
技术领域
本发明涉及光电探测领域,尤其涉及一种红外测量系统自身杂散辐射的分析方法。
背景技术
杂散辐射是指到达光学系统探测器靶面的非目标辐射能量。杂散辐射的危害性在于,其会大大降低系统的信噪比,严重时会使目标信号完全湮没在杂散辐射引起的噪声中,从而使系统探测不到目标。而且,由于目标像中含有杂散辐射能量,若在目标辐射特性测量中不对杂散辐射进行分析与修正,则会影响系统的定量测量精度。
随着大口径、长焦距红外系统的应用,以及探测器灵敏度的大幅度提高,杂散辐射已经成为影响红外光学系统成像与测量的重要因素。与可见光系统不同,红外光学系统不仅包括由外部光源引起的杂散辐射,还包括由系统自身引起的杂散辐射,主要表现为系统内部光学元件及机械结构表面的热辐射。可见,与可见光系统相比,红外光学系统的杂散辐射特性更为复杂。
目前,国内正在加大对杂散光的研究力度,并取得了一定的成果。但是,国内主要关注于对外部杂散光及成像杂散光的分析,而对红外光学系统内部杂散辐射的定量分析工作较少。如此一来,不便于对红外光学系统的杂散辐射进行抑制,进而使得系统的测量精度难以继续提高。
鉴于现有技术中存在对红外光学系统自身杂散辐射进行分析计算的需求,本发明提供了一种红外测量系统自身杂散辐射的分析方法。
发明内容
本发明的目的在于提出一种红外测量系统自身杂散辐射的分析方法,以对杂散辐射与温度的关系进行定量分析,进而为不同温度下杂散辐射的抑制工作提供理论制导,有效提高系统的测量精度。
本发明的红外测量系统自身杂散辐射的分析方法,包括:
S1、确定所述系统内部的杂散辐射源;其中,所述杂散辐射源的总个数为N个;所述杂散辐射源包括:光学元件自身的热辐射、机械结构自身的热辐射;
S2、对于第i个杂散辐射源,计算其辐射面上辐射的、且能够被探测器接收的辐射功率φi;
S3、对φi的传输过程进行分析,计算其在探测器上的辐射照度Ei;
S4、计算N个杂散辐射源在探测器上的总辐射照度E;
其中,N为大于1的整数,i=1,…N。
优选的,在步骤S3中,Ei的计算公式具体为:
式中,Φi为探测器表面接收的辐射功率,Ad为探测器的探测面积,Δi为无穷小量。
优选的,步骤S2包括:
S21、以第i个杂散辐射源表面作为物方,以该辐射面到探测器之间的光学元件作为成像光学系统,计算该成像光学系统的入瞳位置q(xe,ye,ze)、以及入瞳直径D;
S22、计算该辐射面在点p(x,y,z)处的法线向量以及,计算由点p(x,y,z)出射、且过入瞳中心的主光线的方向向量其中,
S23、计算所述法线向量与所述方向向量之间的夹角的余弦值
S24、计算点p(x,y,z)相对入瞳的立体角Ω;
式中,cosθ为主光线的方向向量与光轴的方向向量之间夹角的余弦值,
S25、计算该辐射面辐射的、且能够被探测器接收的辐射功率φi;
式中,S为该辐射面在yz平面上的投影区域,L为该辐射面的辐射亮度。
优选的,所述辐射亮度L具体为:
式中,λ为波长,T为绝对温度,c1、c2为辐射常数,ε为该辐射面的发射率。
优选的,在步骤S22中,计算该辐射面在点p(x,y,z)处的法线向量具体为:
确定该辐射面的面形方程F(x,y,z),
F(x,y,z)=my2+nz2+kx2-2r0x=0;
计算该辐射面在点p(x,y,z)处的法线向量
αN=-(kx-r0)/Δ;
βN=-(my)/Δ;
γN=-(qz)/Δ;
其中,k=1-e2,e为偏心率,r0为顶点的曲率半径,m、n为判断面型的常数;
优选的,在步骤S22中,主光线的方向向量具体为:
α=(xe-x)/Δe;
β=(ye-y)/Δe;
γ=(ze-z)/Δe;
式中,
优选的,所述红外测量系统具体为卡塞格林光谱分光式红外测量系统,所述杂散辐射源包括:主镜、次镜、分色镜、补偿镜、透镜的热辐射,以及,次镜支撑架、分色镜支撑架、补偿镜支撑架的热辐射。
在本发明的技术方案中,红外测量系统自身杂散辐射的分析方法主要包括:确定所述系统中的杂散辐射源;对于第i个杂散辐射源,计算其辐射面上辐射的、且能够被探测器接收的辐射功率φi;对φi的传输过程进行分析,计算其在探测器上的辐射照度Ei;计算N个杂散辐射源在探测器上的总辐射照度E。当红外测量系统设计完成后,其光学元件的面型、厚度、间隔,光学元件与机械结构表面的相对位置以及光学元件与机械结构表面的发射率、反射率等特性参数即能确定。因此,对于一个特定的红外测量系统,E的变化只与辐射面的辐射亮度,即系统的工作温度有关。也就是说,本发明的方法能够对系统自身杂散辐射与温度的关系进行定量分析,从而可为不同温度下杂散辐射的抑制提供理论制导。进一步的,通过在不同温度下有效抑制杂散辐射,能够有效提高系统的测量精度。
附图说明
通过以下参照附图而提供的具体实施方式部分,本发明的特征和优点将变得更加容易理解,在附图中:
图1是本发明的红外测量系统自身杂散辐射的分析方法流程图;
图2是示出的卡塞格林光谱分光式红外测量系统结构示意图;
图3是表面杂散辐射分析原理示意图。
具体实施方式
下面参照附图对本发明的示例性实施方式进行详细描述。对示例性实施方式的描述仅仅是出于示范目的,而绝不是对本发明及其应用或用法的限制。
在红外测量系统中,杂散辐射是影响系统探测性能的重要因素之一。为了在设备研制时采取抑制杂散辐射的相应措施、定量分析抑制效果,通常需要研究人员对杂散辐射进行定量分析。但是,由于红外测量系统的杂散辐射的特性较为复杂,因此,国内主要专注于对外部杂散光和成像杂散光的研究,而对于红外测量系统内部杂散辐射的定量分析较少。
鉴于此,本申请的发明人提出了一种红外测量系统自身杂散辐射的分析方法,以定量分析不同工作温度下的杂散辐射,以便为杂散辐射抑制工作提供理论指导、定量分析杂散辐射的抑制效果,提升系统的探测能力和测量精度。
下面结合附图以及具体实施例对本发明的技术方案进行详细说明。图1示出了本发明的红外测量系统自身杂散辐射的分析方法,具体包括以下步骤:
步骤S1、确定红外测量系统内部的杂散辐射源。其中,所述杂散辐射源的总个数为N个,N为大于1的整数。所述杂散辐射源包括:光学元件自身的热辐射、机械结构自身的热辐射。
具体地,我们以图2所示的地基大口径卡塞格林光谱分光红外测量系统为例进行说明。在该系统中,其光学部分采用以卡塞格林反射系统为前级、透镜组为后级的二级光学方案,其结构部分主要包括主镜筒、共基面组合式三合一小系统。在该系统中,光学元件自身的热辐射主要来自主镜、次镜、分色镜、补偿镜、透镜组等,机械结构自身的热辐射主要来自次镜支撑架、分色镜与光轴垂直面对应的机械结构等。
步骤S2、对于第i个杂散辐射源,计算其辐射面上辐射的、且能够被探测器接收的辐射功率φi。其中,i=1,…N。
具体地,步骤S2包括步骤S21-S25。步骤S21:以第i个杂散辐射源表面作为物方,以该辐射面到探测器之间的光学元件作为成像光学系统,计算该成像光学系统的入瞳位置q(xe,ye,ze)、以及入瞳直径D。
在步骤S21中,我们可根据几何光学原理计算该成像光学系统中的入瞳位置、以及入瞳的大小。由于如何确定入瞳位置、入瞳大小的部分属于本领域的公知常识,在此不再赘述。
步骤S22:计算该辐射面在点p(x,y,z)处的法线向量以及,计算由点p(x,y,z)出射、且过入瞳中心的主光线的方向向量其中,
具体地,在步骤S22中,计算该辐射面在点p(x,y,z)处的法线向量具体为:首先,确定该辐射面的面形方程F(x,y,z)。比如,我们将面形方程设为二次曲面,即F(x,y,z)满足:
F(x,y,z)=my2+nz2+kx2-2r0x=0
式中,k=1-e2,e为偏心率,r0为顶点的曲率半径,m、n为判断面型的常数。在确定了面形方程以后,对F(x,y,z)求偏导,可得:
F′x=2kx-2r0
F′y=2py
F′z=2qz
然后,计算该辐射面在点p(x,y,z)处的法线向量可得:
αN=-(kx-r0)/Δ;
βN=-(my)/Δ;
γN=-(qz)/Δ;
其中,
另外,在步骤S22中,主光线的方向向量具体为:
α=(xe-x)/Δe;
β=(ye-y)/Δe;
γ=(ze-z)/Δe;
式中,
步骤S23、计算所述法线向量与所述方向向量之间的夹角的余弦值具体地,可根据公式3计算
步骤S24、计算点p(x,y,z)相对入瞳的立体角Ω。具体地,可根据公式4计算立体角Ω。
式中,cosθ为主光线的方向向量与光轴的方向向量之间夹角的余弦值,故cosθ=α。
步骤S25、计算该辐射面辐射的、且能够被探测器接收的辐射功率φi。具体地,可根据公式5计算φi。
式中,S为该辐射面在yz平面上的投影区域,L为该辐射面的辐射亮度。
另外,在步骤S2中,我们可基于红外辐射理论计算辐射亮度L。其计算公式为:
式中,λ为波长,T为绝对温度,c1、c2为辐射常数,ε为该辐射面的发射率。
由于红外测量系统中的光学元件具有一定的反射率或透过率,因此光学元件或机械结构表面的杂散辐射在传输过程中会有一定的能量损失。因此,我们需要对各个杂散辐射源的传递过程进行分析,即步骤S3。
步骤S3、对φi的传输过程进行分析,计算其在探测器上的辐射照度Ei。其中,Ei的计算公式具体为:
式中,Φi为探测器表面接收的辐射功率,Ad为探测器的探测面积,Δi为无穷小量。
在步骤S3中,为了便于说明,我们设主镜表面为A1、主镜的辐射亮度为L1、主镜的反射率为ρ1;次镜表面为A2、次镜的辐射亮度为L2、次镜的反射率为ρ2;分色镜表面为A3、辐射亮度L3、透过率为τ1;补偿镜表面为A4、辐射亮度为L4、透过率为τ2;透镜表面为A5、辐射亮度为L5、透过率为τ3;红外热像仪镜组、窗口的总透过率为τ4;次镜支撑架截面为A6、辐射亮度为L6;分色镜与补偿镜支撑架的结构表面为A7、辐射亮度为L7;分色镜所对应的结构表面为A8、辐射亮度为L8。
下面我们以主镜自身杂散辐射的传递过程为例进行说明。主镜表面A1的辐射功率为φ1,φ1经过A2、A3、A4、A5之后,其在探测器靶面的辐射照度E1为:
式中,Δ1是无穷小量。
同理,我们可根据以上方法计算次镜、分色镜等其他杂散辐射源在探测器靶面上的辐射照度。
步骤S4、计算N个杂散辐射源在探测器上的总辐射照度E。具体地,我们可根据公式1计算总辐射照度E
具体地,在地基大口径卡塞格林光谱分光红外测量系统中,其总辐射照度E为:
式中,Δ为经过多次散射、反射后的无穷小量。为了简便,我们仅在上式中给出了E1、E2、E8的展开形式,其他杂散辐射源的辐射照度可根据以上方法得出。然后,基于公式5、7,我们可得到探测器上的总辐射照度E具体为:
当红外测量系统设计完成后,其光学元件的面型、厚度、间隔,光学元件与机械结构表面的相对位置以及光学元件与机械结构表面的发射率、反射率等特性参数即能确定。因此,我们可将上式大括号内表征光机结构与表面属性的参数用常量C表示,并将公式6代入上式,由此可得:
在该具体实施例中,我们通过以上方法定量计算了地基大口径卡塞格林光谱分光红外测量系统内部的杂散辐射。通过对系统内部杂散辐射进行定量分析,得到了该系统杂散辐射随温度的变化关系,从而为不同工作温度下杂散辐射的抑制提供了理论指导。进一步的通过有效抑制杂散辐射,提高了系统的测量精度。
虽然参照示例性实施方式对本发明进行了描述,但是应当理解,本发明并不局限于文中详细描述和示出的具体实施方式,在不偏离权利要求书所限定的范围的情况下,本领域技术人员可以对所述示例性实施方式做出各种改变。
Claims (7)
1.一种红外测量系统自身杂散辐射的分析方法,其特征在于,所述方法包括:
S1、确定所述系统内部的杂散辐射源;其中,所述杂散辐射源的总个数为N个;所述杂散辐射源包括:光学元件自身的热辐射、机械结构自身的热辐射;
S2、对于第i个杂散辐射源,计算其辐射面上辐射的、且能够被探测器接收的辐射功率φi;
S3、对φi的传输过程进行分析,计算其在探测器上的辐射照度Ei;
S4、计算N个杂散辐射源在探测器上的总辐射照度E;
其中,N为大于1的整数,i=1,…N。
2.如权利要求1所述的方法,其特征在于,在步骤S3中,Ei的计算公式具体为:
式中,Φi为探测器表面接收的辐射功率,Ad为探测器的探测面积,Δi为无穷小量。
3.如权利要求2所述的方法,其特征在于,步骤S2包括:
S21、以第i个杂散辐射源表面作为物方,以该辐射面到探测器之间的光学元件作为成像光学系统,计算该成像光学系统的入瞳位置q(xe,ye,ze)、以及入瞳直径D;
S22、计算该辐射面在点p(x,y,z)处的法线向量以及,计算由点p(x,y,z)出射、且过入瞳中心的主光线的方向向量其中,
S23、计算所述法线向量与所述方向向量之间的夹角的余弦值
S24、计算点p(x,y,z)相对入瞳的立体角Ω;
式中,cosθ为主光线的方向向量与光轴的方向向量之间夹角的余弦值,
S25、计算该辐射面辐射的、且能够被探测器接收的辐射功率φi;
式中,S为该辐射面在yz平面上的投影区域,L为该辐射面的辐射亮度。
4.如权利要求3所述的方法,其特征在于,所述辐射亮度L具体为:
式中,λ为波长,T为绝对温度,c1、c2为辐射常数,ε为该辐射面的发射率。
5.如权利要求4所述的方法,其特征在于,在步骤S22中,计算该辐射面在点p(x,y,z)处的法线向量具体为:
确定该辐射面的面形方程F(x,y,z),
F(x,y,z)=my2+nz2+kx2-2r0x=0;
计算该辐射面在点p(x,y,z)处的法线向量
αN=-(kx-r0)/Δ;
βN=-(my)/Δ;
γN=-(qz)/Δ;
其中,k=1-e2,e为偏心率,r0为顶点的曲率半径,m、n为判断面型的常数;
6.如权利要求5所述的方法,其特征在于,在步骤S22中,主光线的方向向量具体为:
α=(xe-x)/Δe;
β=(ye-y)/Δe;
γ=(ze-z)/Δe;
式中,
7.如权利要求1-6任一所述的方法,其特征在于,所述红外测量系统具体为卡塞格林光谱分光式红外测量系统,所述杂散辐射源包括:主镜、次镜、分色镜、补偿镜、透镜的热辐射,以及,次镜支撑架、分色镜支撑架、补偿镜支撑架的热辐射。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390112.1A CN106053023B (zh) | 2016-06-03 | 2016-06-03 | 一种红外测量系统自身杂散辐射的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610390112.1A CN106053023B (zh) | 2016-06-03 | 2016-06-03 | 一种红外测量系统自身杂散辐射的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106053023A true CN106053023A (zh) | 2016-10-26 |
CN106053023B CN106053023B (zh) | 2018-08-14 |
Family
ID=57170034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610390112.1A Active CN106053023B (zh) | 2016-06-03 | 2016-06-03 | 一种红外测量系统自身杂散辐射的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106053023B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297590A (zh) * | 2018-09-14 | 2019-02-01 | 北京遥感设备研究所 | 一种红外探测器背景辐射噪声计算方法 |
CN109341863A (zh) * | 2018-10-12 | 2019-02-15 | 上海卫星工程研究所 | 长波红外相机辐射定标数据修正方法 |
CN110455417A (zh) * | 2019-08-21 | 2019-11-15 | 北京环境特性研究所 | 针对红外光学系统杂散辐射的定量测量误差校正方法 |
CN112485901A (zh) * | 2020-12-14 | 2021-03-12 | 长光卫星技术有限公司 | 一种快速获取杂光来源和杂光系数的仿真建模方法 |
CN113008377A (zh) * | 2021-02-22 | 2021-06-22 | 中国电子科技集团公司第十一研究所 | 红外光学系统杂散辐射的分析方法和抑制方法 |
CN114152352A (zh) * | 2021-11-17 | 2022-03-08 | 北京理工大学 | 一种红外光学系统杂散辐射实验测量方法及系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5382785A (en) * | 1992-05-04 | 1995-01-17 | Diolase Corporation | Laser beam delivery path and target proximity sensor |
CN103149016A (zh) * | 2013-02-27 | 2013-06-12 | 中国科学院西安光学精密机械研究所 | 待测光学系统杂散光检测方法及杂散光检测系统 |
CN103234734B (zh) * | 2013-04-09 | 2016-06-08 | 中国科学院西安光学精密机械研究所 | 大口径杂散光测试装置及测试方法 |
CN105547649B (zh) * | 2015-12-05 | 2019-03-12 | 中国航空工业集团公司洛阳电光设备研究所 | 一种短波红外镜头杂散辐射的检测方法 |
-
2016
- 2016-06-03 CN CN201610390112.1A patent/CN106053023B/zh active Active
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297590A (zh) * | 2018-09-14 | 2019-02-01 | 北京遥感设备研究所 | 一种红外探测器背景辐射噪声计算方法 |
CN109297590B (zh) * | 2018-09-14 | 2021-01-01 | 北京遥感设备研究所 | 一种红外探测器背景辐射噪声计算方法 |
CN109341863A (zh) * | 2018-10-12 | 2019-02-15 | 上海卫星工程研究所 | 长波红外相机辐射定标数据修正方法 |
CN110455417A (zh) * | 2019-08-21 | 2019-11-15 | 北京环境特性研究所 | 针对红外光学系统杂散辐射的定量测量误差校正方法 |
CN112485901A (zh) * | 2020-12-14 | 2021-03-12 | 长光卫星技术有限公司 | 一种快速获取杂光来源和杂光系数的仿真建模方法 |
CN113008377A (zh) * | 2021-02-22 | 2021-06-22 | 中国电子科技集团公司第十一研究所 | 红外光学系统杂散辐射的分析方法和抑制方法 |
CN114152352A (zh) * | 2021-11-17 | 2022-03-08 | 北京理工大学 | 一种红外光学系统杂散辐射实验测量方法及系统 |
CN114152352B (zh) * | 2021-11-17 | 2024-04-19 | 北京理工大学 | 一种红外光学系统杂散辐射实验测量方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106053023B (zh) | 2018-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106053023A (zh) | 一种红外测量系统自身杂散辐射的分析方法 | |
Boccaletti et al. | Multiple spiral patterns in the transitional disk of HD 100546 | |
CN106643796B (zh) | 一种基于在轨基准星的辐射定标方法 | |
CN102374902B (zh) | 一种提高辐射温度计测温准确度的量子论修正方法 | |
EP2888568B1 (en) | Polarimetric calibration of a remote sensor | |
CN104880702A (zh) | 在轨绝对辐射定标方法及装置 | |
Close et al. | New Photometry and Spectra of AB Doradus C: An Accurate Mass Determination of a Young Low-Mass Object with Theoretical Evolutionary Tracks | |
WO2018086161A1 (zh) | 一种大口径辐射计的光谱响应度测量装置及方法 | |
Fulbright et al. | Suomi-NPP VIIRS solar diffuser stability monitor performance | |
CN103852078A (zh) | 空间光学姿态敏感器杂散光保护角的测量方法及装置 | |
Mourard et al. | A high angular and spectral resolution view into the hidden companion of ε Aurigae | |
Du et al. | The attitude inversion method of geostationary satellites based on unscented particle filter | |
CN104714222B (zh) | 激光雷达系统回波能量的计算方法 | |
CN103134443B (zh) | 一种大口径大径厚比反射镜面形自准直检测装置及方法 | |
Beck et al. | Thermal analysis of a reflective baffle designed for space applications | |
CN104714221B (zh) | 激光雷达系统回波能量动态范围的压缩方法 | |
CN103727962B (zh) | 大视场红外光电经纬仪精度标定方法 | |
Pfänder et al. | Pyrometric temperature measurements on solar thermal high temperature receivers | |
US7576346B1 (en) | Use of solar illumination and baffles to calibrate visible sensors | |
CN205642633U (zh) | 基于光辐射的测温设备及系统 | |
CN205642623U (zh) | 基于编码变换的双波长温度场成像设备及系统 | |
Tian et al. | Spherical warm shield design for infrared imaging systems | |
Wang et al. | Research of suppressing stray light of cold shield with different black coating | |
Huang et al. | Research on accurate non-contact temperature measurement method for telescope mirror | |
Borodin et al. | Imaging Atmospheric Cherenkov Telescope for the TAIGA Observatory—JINR Participation |
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 |