CN108226947A - 基于最佳孔径滤除噪声的涡旋光测距系统 - Google Patents

基于最佳孔径滤除噪声的涡旋光测距系统 Download PDF

Info

Publication number
CN108226947A
CN108226947A CN201810084480.2A CN201810084480A CN108226947A CN 108226947 A CN108226947 A CN 108226947A CN 201810084480 A CN201810084480 A CN 201810084480A CN 108226947 A CN108226947 A CN 108226947A
Authority
CN
China
Prior art keywords
noise
phase plate
spiral phase
signal
vortex
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
CN201810084480.2A
Other languages
English (en)
Other versions
CN108226947B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201810084480.2A priority Critical patent/CN108226947B/zh
Publication of CN108226947A publication Critical patent/CN108226947A/zh
Application granted granted Critical
Publication of CN108226947B publication Critical patent/CN108226947B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/08Systems determining position data of a target for measuring distance only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/481Constructional features, e.g. arrangements of optical elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

为了解决现有测距系统无法有效滤除背景噪声的问题,本发明提供一种基于最佳孔径滤除噪声的涡旋光测距系统,属于量子信息技术领域。所述测距系统包括:激光器发射出的激光经过分光棱镜后被分成两束,一束作为触发信号进入PIN型光电二极管,进而触发处理器开始计时;另一束入射至第一螺旋相位板、第一卡塞格林望远镜至目标,经目标后产生回波信号;回波信号经第二卡塞格林望远镜接收,依次经第二螺旋相位板、可调光阑、探测器探测,处理器对探测后电信号后计时结束;根据螺旋相位板的阶数,设计可调光阑的通光孔径,使探测器的输出功率信噪比等于,此时可调光阑的通光孔径为最佳孔径,提高了测量精度。

Description

基于最佳孔径滤除噪声的涡旋光测距系统
技术领域
本发明涉及属于量子信息技术领域,特别是一种适合于远距离测距系统中背景噪声的滤除以及接收孔径的设计的系统。
背景技术
利用飞行时间法进行远距离探测是一个探测领域热情不减的研究领域,但是伴随而来的背景噪声是一个不易解决的困扰,这一定程度上限制了该项技术的使用,以往的做法通常是通过滤波片滤出一部分背景噪声,但是与发射信号相同波长的背景噪声依然无法有效滤除。因此急需一种可以实现高效滤除背景噪声的方法。
发明内容
本发明的目的是为了解决现有测距系统无法有效滤除背景噪声的问题,本发明提供一种基于最佳孔径滤除噪声的涡旋光测距系统。
本发明的基于最佳孔径滤除噪声的涡旋光测距系统,所述测距系统包括激光器1、分光棱镜2、PIN型光电二极管3、第一螺旋相位板4-1、第二螺旋相位板4-2、第一望远镜系统5-1、第二望远镜系统5-2、可调光阑7、探测器8和处理器9;
激光器1发射出的激光经过分光棱镜2后被分成两束,一束作为触发信号进入PIN型光电二极管3,进而触发处理器9开始计时;另一束入射至第一螺旋相位板4-1,该第一螺旋相位板4-1发射出涡旋光束至第一望远镜系统5-1,经第一望远镜系统5-1扩束后发射至目标10,经目标10后产生回波信号;
所述回波信号经第二望远镜系统5-2接收,第二望远镜系统5-2接收的光信号经第二螺旋相位板4-2进行解调后,经可调光阑7将解调后的光信号中的背景噪声滤除,滤除噪声后的光信号通过光纤耦合进探测器8探测,探测后电信号发送至处理器9,此时处理器9计时结束;
所述第一螺旋相位板4-1和第二螺旋相位板4-2的阶数相反;
根据第一螺旋相位板4-1的阶数,设计可调光阑7的通光孔径,使探测器8的输出功率信噪比等于1,此时可调光阑7的通光孔径为最佳孔径。
优选的是,所述测距系统还包括凸透镜6;
凸透镜6设置在可调光阑7和第二螺旋相位板4-2之间,第二螺旋相位板4-2解调后的光信号经过凸透镜6后入射至可调光阑7。
优选的是,探测器8的输出功率信噪比等于1的条件为:
W0表示背景噪声光信号经过可调光阑7通光孔径的通光量;
Ws表示涡旋光信号经过可调光阑7通光孔径的通光量。
优选的是,所述第一望远镜系统5-1和第二望远镜系统5-2均采用卡塞格林望远镜实现。
优选的是,所述激光器1为单模激光源。
上述技术特征可以各种适合的方式组合或由等效的技术特征来替代,只要能够达到本发明的目的。
本发明的有益效果在于,本发明能够对远距离探测情况下的自然光等背景噪声进行滤除。本发明将带有轨道角动量的涡旋光束作为探测光束,通过望远镜系统进行扩束发射进行目标10探测,在接收端通过望远镜系统对回波信号进行接收,并对接收到的信号光和背景噪声都进行解调,并通过最佳孔径设计让信号光通过并同时滤除背景噪声,保证系统实现高信噪比。本发明适合于远距离测距系统中背景噪声的滤除以及接收孔径的设计。本发明利用处理器计时的时长t和光速c就可以求出目标10距离L:通过设计可调光阑7的最佳孔径滤除探测过程中的背景噪声,提高测量精度。
附图说明
图1为本发明的基于最佳孔径滤除噪声的涡旋光测距系统的原理示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
下面结合附图和具体实施例对本发明作进一步说明,但不作为本发明的限定。
结合图1说明本实施方式,本实施方式所述的基于最佳孔径滤除噪声的涡旋光测距系统,包括激光器1、分光棱镜2、PIN型光电二极管3、第一螺旋相位板4-1、第二螺旋相位板4-2、第一卡塞格林望远镜5-1、第二卡塞格林望远镜5-2、可调光阑7、探测器8和处理器9;
激光器1发射出的激光经过分光棱镜2后被分成两束,一束作为触发信号进入PIN型光电二极管3,进而触发处理器9开始计时;另一束经过+s阶的第一螺旋相位板4-1产生涡旋光束,该涡旋光束通过第一卡塞格林望远镜5-1,经第一卡塞格林望远镜5-1扩束后发射至目标109,经目标109后产生回波信号;回波信号包括回波信号中含有涡旋光信号和背景噪声两种成分,均通过第二卡塞格林望远镜5-2接收,随后接收信号通过-s阶的第二螺旋相位板4-2进行解调,此时信号被解调回0阶,而噪声光束被调制成-s阶的涡旋光束。此后经过可调光阑7将背景噪声进行滤除,并将滤除噪声后的信号通过光纤耦合的探测器8,随后送达至处理器9,处理器9停止计时;
所述第一螺旋相位板4-1和第二螺旋相位板4-2的阶数相反;
根据第一螺旋相位板4-1的阶数,设计可调光阑7的通光孔径,使探测器8的输出功率信噪比等于1,此时可调光阑7的通光孔径为最佳孔径。
本实施方式利用处理器计时的时长t和光速c就可以求出目标10距离L:通过设计可调光阑7的最佳孔径滤除探测过程中的背景噪声,提高测量精度。
涡旋光束其振幅表达式为:
其中,r表示径向位移,θ表示径向角度,z表示轴向位移,k表示波矢,!表示阶乘;称为Gouy相位,称为Rayleigh长度,ω0为束腰半径,是z处的光束半径,是z处光束的曲率半径,是归一化系数,是缔合拉盖尔多项式,p是表征模式的特征量子数。
将上式振幅模平方得到光强分布:
设通光孔径半径为R,令取a=-2/ω2,则通光量为:
系统探测方式为直接探测,设入射信号功率为Pso,噪声功率为Pno,根据光探测器平方律特性,光探测器输出的信号电功率Pso和噪声功率Pno为:
e表示电子电荷量,η表示探测器的量子效率,h表示普朗克常数,ν表示光子频率,RL表示探测器的负载电阻阻值,Psi表示光探测器输入的涡旋光信号的功率,Pni表示光探测器输入的噪声功率;
因此,系统输出功率信噪比公式为:
选定信噪比为1式的通光半径为临界值,上式等于1的条件为
所以,当发射为+s阶涡旋光束时,只需计算出R,使得此时即为:最佳孔径。
W0表示背景噪声光信号经过可调光阑7通光孔径的通光量;
Ws表示涡旋光信号经过可调光阑7通光孔径的通光量。
由于解调的调制需要一定的演化距离,因此这一过程需要添加凸透镜来加速演化。所以优选实施例中,本实施方式的测距系统还包括凸透镜6;
凸透镜6设置在可调光阑7和第二螺旋相位板4-2之间,第二螺旋相位板4-2解调后的光信号经过凸透镜6后入射至可调光阑7。
优选实施例中,本实施方式的激光器1为单模激光源。
本实施方式中单模激光源1:具有较高的功率稳定性和较好的空间模式,输出的是单模高斯光,激光单色性好,脉冲脉宽小于10ns。
分光棱镜2:内部对角面镀宽带分光膜,具有10:90(R:T)的分光比。基底材料为N-BK7。
PIN型光电二极管3:超快响应(1.5GHz),从10-900MHz的3dB带宽光谱,自由空间模式输入。
第一螺旋相位板4-1、第二螺旋相位板4-2:镀有单一波长增透膜,曲面高度连续变化。引入相位误差小于五十分之一波长。
第一卡塞格林望远镜5-1、第二卡塞格林望远镜5-2:引入的波前畸变小于四分之一波长,镀有单一波长增透膜,具有5倍以上的放大倍率。
凸透镜6:引入的波前畸变小于十分之一波长,双面均镀有单一波长增透膜,材料为N-BK7。
可调光阑7:孔径直径的连续变化范围为0.1mm-20mm。
探测器8:使用盖革模式的单光子探测器,对信号波长的探测效率大于80%,暗计数小于100c/s。
处理器9:带有计时功能的2通道数据采集卡。
虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其他所述实施例中。

Claims (5)

1.一种基于最佳孔径滤除噪声的涡旋光测距系统,其特征在于,所述测距系统包括激光器(1)、分光棱镜(2)、PIN型光电二极管(3)、第一螺旋相位板(4-1)、第二螺旋相位板(4-2)、第一望远镜系统(5-1)、第二望远镜系统(5-2)、可调光阑(7)、探测器(8)和处理器(9);
激光器(1)发射出的激光经过分光棱镜(2)后被分成两束,一束作为触发信号进入PIN型光电二极管(3),进而触发处理器(9)开始计时;另一束入射至第一螺旋相位板(4-1),该第一螺旋相位板(4-1)发射出涡旋光束至第一望远镜系统(5-1),经第一望远镜系统(5-1)扩束后发射至目标,经目标后产生回波信号;
所述回波信号经第二望远镜系统(5-2)接收,第二望远镜系统(5-2)接收的光信号经第二螺旋相位板(4-2)进行解调后,经可调光阑(7)将解调后的光信号中的背景噪声滤除,滤除噪声后的光信号通过光纤耦合进探测器(8)探测,探测后电信号发送至处理器(9),此时处理器(9)计时结束;
所述第一螺旋相位板(4-1)和第二螺旋相位板(4-2)的阶数相反;
根据第一螺旋相位板(4-1)的阶数,设计可调光阑(7)的通光孔径,使探测器(8)的输出功率信噪比等于1,此时可调光阑(7)的通光孔径为最佳孔径。
2.根据权利要求1所述的基于最佳孔径滤除噪声的涡旋光测距系统,其特征在于,
探测器(8)的输出功率信噪比等于1的条件为:
W0表示背景噪声光信号经过可调光阑(7)通光孔径的通光量;
Ws表示涡旋光信号经过可调光阑(7)通光孔径的通光量。
3.根据权利要求1或2所述的基于最佳孔径滤除噪声的涡旋光测距系统,其特征在于,所述测距系统还包括凸透镜(6);
凸透镜(6)设置在可调光阑(7)和第二螺旋相位板(4-2)之间,第二螺旋相位板(4-2)解调后的光信号经过凸透镜(6)后入射至可调光阑(7)。
4.根据权利要求1所述的基于最佳孔径滤除噪声的涡旋光测距系统,其特征在于,所述第一望远镜系统(5-1)和第二望远镜系统(5-2)均采用卡塞格林望远镜实现。
5.根据权利要求1所述的基于最佳孔径滤除噪声的涡旋光测距系统,其特征在于,所述激光器(1)为单模激光源。
CN201810084480.2A 2018-01-29 2018-01-29 基于最佳孔径滤除噪声的涡旋光测距系统 Expired - Fee Related CN108226947B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810084480.2A CN108226947B (zh) 2018-01-29 2018-01-29 基于最佳孔径滤除噪声的涡旋光测距系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810084480.2A CN108226947B (zh) 2018-01-29 2018-01-29 基于最佳孔径滤除噪声的涡旋光测距系统

Publications (2)

Publication Number Publication Date
CN108226947A true CN108226947A (zh) 2018-06-29
CN108226947B CN108226947B (zh) 2021-09-24

Family

ID=62669127

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810084480.2A Expired - Fee Related CN108226947B (zh) 2018-01-29 2018-01-29 基于最佳孔径滤除噪声的涡旋光测距系统

Country Status (1)

Country Link
CN (1) CN108226947B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111355530A (zh) * 2020-03-13 2020-06-30 山东师范大学 一种提高无线光通信系统性能的方法及装置
CN112327279A (zh) * 2020-10-30 2021-02-05 哈尔滨工业大学 基于轨道角动量调制的抗云雾后向散射激光探测系统
CN112859099A (zh) * 2021-02-25 2021-05-28 哈尔滨工业大学 云雾环境中基于环形分离器实时补偿的高精度测距系统
CN112859044A (zh) * 2021-02-19 2021-05-28 北京理工大学 一种基于涡旋光的水下激光雷达系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1296237A (zh) * 1999-11-12 2001-05-23 惠普公司 具有可变光圈的扫描器引导系统
JP2001343266A (ja) * 2000-04-24 2001-12-14 Chang Min Tech Co Ltd 超音波流速測定装置
CN102213763A (zh) * 2011-04-11 2011-10-12 哈尔滨工业大学 基于锁模激光器的相干多普勒测风激光雷达测距系统及测距方法
CN102404049A (zh) * 2011-12-23 2012-04-04 中国人民解放军理工大学 基于云散射的无线光通信系统
CN104065418A (zh) * 2014-06-05 2014-09-24 华中科技大学 一种集成轨道角动量模式发射器
CN105005053A (zh) * 2015-07-13 2015-10-28 西安电子科技大学 基于led照明的随机散射关联成像系统及成像方法
WO2016111861A1 (en) * 2015-01-08 2016-07-14 Nec Laboratories America, Inc. Remote wind turbulence sensing
CN105785323A (zh) * 2016-01-31 2016-07-20 西安电子科技大学 基于空间正交变换的涡旋电磁波信号模态并行分离方法与装置
CN106950557A (zh) * 2017-04-26 2017-07-14 哈尔滨工业大学 一种基于光子轨道角动量调制的单光子测距背景噪声滤除方法及单光子测距装置
CN206348536U (zh) * 2016-11-16 2017-07-21 华侨大学 产生阶数可调的焦散光束的光学系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1296237A (zh) * 1999-11-12 2001-05-23 惠普公司 具有可变光圈的扫描器引导系统
JP2001343266A (ja) * 2000-04-24 2001-12-14 Chang Min Tech Co Ltd 超音波流速測定装置
CN102213763A (zh) * 2011-04-11 2011-10-12 哈尔滨工业大学 基于锁模激光器的相干多普勒测风激光雷达测距系统及测距方法
CN102404049A (zh) * 2011-12-23 2012-04-04 中国人民解放军理工大学 基于云散射的无线光通信系统
CN104065418A (zh) * 2014-06-05 2014-09-24 华中科技大学 一种集成轨道角动量模式发射器
WO2016111861A1 (en) * 2015-01-08 2016-07-14 Nec Laboratories America, Inc. Remote wind turbulence sensing
CN105005053A (zh) * 2015-07-13 2015-10-28 西安电子科技大学 基于led照明的随机散射关联成像系统及成像方法
CN105785323A (zh) * 2016-01-31 2016-07-20 西安电子科技大学 基于空间正交变换的涡旋电磁波信号模态并行分离方法与装置
CN206348536U (zh) * 2016-11-16 2017-07-21 华侨大学 产生阶数可调的焦散光束的光学系统
CN106950557A (zh) * 2017-04-26 2017-07-14 哈尔滨工业大学 一种基于光子轨道角动量调制的单光子测距背景噪声滤除方法及单光子测距装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CARMEL ROTSCHILD 等: ""Adjustable spiral phase plate"", 《APPLIED OPTICS》 *
HUIQIN TANG 等: ""Vortex characteristics of Fraunhofer diffractions of a plane wave by a spiral phase plate limited by pseudoring polygonal apertures"", 《OPTICS LETTERS》 *
赵远 等: ""增益调制非扫描激光雷达测距精度的理论分析"", 《红外与激光工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111355530A (zh) * 2020-03-13 2020-06-30 山东师范大学 一种提高无线光通信系统性能的方法及装置
CN112327279A (zh) * 2020-10-30 2021-02-05 哈尔滨工业大学 基于轨道角动量调制的抗云雾后向散射激光探测系统
CN112327279B (zh) * 2020-10-30 2023-10-10 哈尔滨工业大学 基于轨道角动量调制的抗云雾后向散射激光探测系统
CN112859044A (zh) * 2021-02-19 2021-05-28 北京理工大学 一种基于涡旋光的水下激光雷达系统
CN112859044B (zh) * 2021-02-19 2023-12-29 北京理工大学 一种基于涡旋光的水下激光雷达系统
CN112859099A (zh) * 2021-02-25 2021-05-28 哈尔滨工业大学 云雾环境中基于环形分离器实时补偿的高精度测距系统
CN112859099B (zh) * 2021-02-25 2023-11-10 哈尔滨工业大学 云雾环境中基于环形分离器实时补偿的高精度测距系统

Also Published As

Publication number Publication date
CN108226947B (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN108226947A (zh) 基于最佳孔径滤除噪声的涡旋光测距系统
Basu et al. Dictionary of pure and applied physics
CN201298079Y (zh) 激光测距装置
US3535016A (en) Electromagnetic radiation concentrating apparatus embodying frustoconical mirror elements
Smith et al. The optical and near-infrared polarization properties of the OVV quasar 3C 345
Burbidge et al. A Further Investigation of Stephan's Quintet.
CN106950557A (zh) 一种基于光子轨道角动量调制的单光子测距背景噪声滤除方法及单光子测距装置
Skinner et al. The enigmatic T Tauri radio source
CN109781257A (zh) 一种自由曲面棱镜光谱仪光学系统
CN109579639A (zh) 半主动激光导引头
CN207081834U (zh) 一种对射式光电传感器
CN108387974A (zh) 一种高功率光纤激光收发一体端帽
CN108710137A (zh) 大视场合成孔径激光成像雷达收发系统
CN209640493U (zh) 一种激光雷达系统
Baum Focal waveform of a prolate-spheroidal IRA
Sandell et al. H2O masers associated with bright nebulosities in dark clouds
White Studying clusters with Planck
CN208076758U (zh) 一种高功率光纤激光收发一体端帽
Peng et al. Multi-target free-space laser communication system based on a rotating double prism with an RF signal beacon
Dent et al. Near-infrared observations of the bipolar outflow source G35. 2N
CN107092041A (zh) 一种激光告警器
Kallmann Optar, a method of optical automatic ranging, as applied to a guidance device for the blind
Shugarov et al. Study of the Classical Novae V2467 Cyg and V2468 Cyg
Hartman et al. Gamma-ray blazars
CN109682763A (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

Granted publication date: 20210924

Termination date: 20220129

CF01 Termination of patent right due to non-payment of annual fee