CN111220177A - 一种星点像亚像元定位精度的验证方法 - Google Patents

一种星点像亚像元定位精度的验证方法 Download PDF

Info

Publication number
CN111220177A
CN111220177A CN201811410797.7A CN201811410797A CN111220177A CN 111220177 A CN111220177 A CN 111220177A CN 201811410797 A CN201811410797 A CN 201811410797A CN 111220177 A CN111220177 A CN 111220177A
Authority
CN
China
Prior art keywords
star
sub
adjacent
star point
distance
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
CN201811410797.7A
Other languages
English (en)
Other versions
CN111220177B (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.)
Changchun Institute of Optics Fine Mechanics and Physics of CAS
Original Assignee
Changchun Institute of Optics Fine Mechanics and Physics of CAS
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 Changchun Institute of Optics Fine Mechanics and Physics of CAS filed Critical Changchun Institute of Optics Fine Mechanics and Physics of CAS
Priority to CN201811410797.7A priority Critical patent/CN111220177B/zh
Publication of CN111220177A publication Critical patent/CN111220177A/zh
Application granted granted Critical
Publication of CN111220177B publication Critical patent/CN111220177B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明实施例公开一种星点像亚像元定位精度的验证方法。该验证方法包括确定地球自转360度的周期，待验证光学***的焦距；将光学***的主光轴对准赤纬为第一预设值的恒星进行积分，在预设时间间隔内进行多次观测并且在所述观测过程中所述光学***保持静止不动；采用待验证的亚像元定位算法计算各个星点像的位置并采用距离公式计算相邻两个星像点之间的距离；求取第一次观测时所述恒星的入射光线与光线***的主光轴的夹角；计算两次相邻观测时间点中地球转过的角度、视轴夹角，以及星像点在像面上的移动距离；计算所述待验证的亚像元定位算法的定位精度。该验证方法具有稳定性和高精度。

Description

一种星点像亚像元定位精度的验证方法
技术领域
本发明涉及遥感器光学设计的技术领域，具体涉及一种星点像亚像元定位精度的验证方法。
背景技术
多种以星点像位置信息为观测目标的光学仪器均需要对星点像的定位精度进行测试实验。目前，测试实验的方式包括：采用地面设备和光源进行模拟测试；采用二维转台进行高精度扫描拍摄空间目标的。对于作为定位精度要求很高的仪器，高精度地面的模拟设备不仅需要很大的成本，并且模拟设备的精度会影响整个测量过程精度、影响仪器定位精度的测试。
因此，针对现有技术中定位精度高的模拟设备成本高及高精度稳定性的问题，需要提供一种具有稳定性高精度的星点像亚像元定位精度的验证方法。
发明内容
针对现有技术中定位精度高的模拟设备成本高及高精度稳定性的问题，本发明实施例提供一种具有稳定性高精度的星点像亚像元定位精度的验证方法。该星点像亚像元定位精度的验证方法利用地球作为高精度转动平台观测恒星的星点像，从而能够提供超高精度的扫描精度，并且地球作为恒星的星点像的张角小于0.01角秒，进而能够为仪器提供高精度的点光源。
该星点像亚像元定位精度的验证方法的具体方案如下：一种星点像亚像元定位精度的验证方法包括步骤S1：确定地球自转360度的周期为T，待验证光学***的焦距为f；步骤S2：将所述光学***的主光轴对准赤纬为第一预设值的恒星进行积分，积分时间为t0，在t时间间隔内进行n次观测并且在所述观测过程中所述光学***保持静止不动，其中，t大于t0；步骤S3：采用待验证的亚像元定位算法计算各个星点像的位置，获得各个星点像的位置坐标，并采用距离公式计算相邻两个星像点之间的距离；步骤S4：求取步骤S2中的第1次观测时所述恒星的入射光线与所述光线***的主光轴的夹角；步骤S5：计算所述两次相邻观测时间点中地球转过的角度和所述两次相邻观测的视轴夹角，根据所述地球转过的角度和所述两次相邻观测的视轴夹角计算相邻两次观测的星像点在像面上的移动距离；步骤S6：根据所述相邻两次观测的星像点在像面上的移动距离，计算所述待验证的亚像元定位算法的定位精度。
优选地，所述步骤S3中的距离公式的表达式为
Figure BDA0001878417790000021
其中，lk'为第k次两个相邻星像点之间的距离，(xk,yk)为第k个星点像的位置坐标。
优选地，所述步骤S4中的求取公式为
Figure BDA0001878417790000022
其中，α为所述恒星的入射光线与所述光线***的主光轴的夹角，d为星像点到相机焦点的距离，f为待验证光学***的焦距。
优选地，所述步骤S5中所述两次相邻观测时间点中地球转过的角度的计算公式为Δ＝360·t/T，其中，t为进行n次观测的时间间隔长度，T为地球自转360度的周期。
优选地，所述步骤S5中所述两次相邻观测的视轴夹角的计算公式为
Figure BDA0001878417790000023
其中，Δ为两次相邻观测时间点中地球转过的角度，β为步骤S2中的所述第一预设值。
优选地，步骤S5中的所述移动距离的计算公式为
lk＝f·(tan(α+kx)-tan(α+(k-1)x),...k∈[1..n-1]，其中，f为待验证光学***的焦距，α为所述恒星的入射光线与所述光线***的主光轴的夹角，x为两次相邻观测的视轴夹角。
优选地，所述定位精度包括误差均值和方差。
优选地，所述误差均值的计算公式为：
Figure BDA0001878417790000031
所述方差的计算公式为：
Figure BDA0001878417790000032
其中，lk'为第k次两个相邻星像点之间的距离，lk为第k次两个相邻星像点在像面上的移动距离。
优选地，所述n为大于2的自然数。
从以上技术方案可以看出，本发明实施例具有以下优点：
本发明实施例提供的一种星点像亚像元定位精度的验证方法利用地球作为高精度转动平台来观测恒星的星点像，从而能够提供高精度的扫描精度，并且恒星的张角能够维持在0.01角秒以内，进而能够为光学***提供高精度的点光源。本发明实施例提供的一种星点像亚像元定位精度的验证方法能够适用于测试光轴与本体间没有相对运动的光学***的星点像定位精度。
附图说明
图1为本发明实施例中提供的一种星点像亚像元定位精度的验证方法的流程示意图；
图2为本发明实施例中提供的天狼星弥散斑定位拍摄的部分结果的示意图；
图3为本发明实施例中提供的一种星点像亚像元定位精度的验证方法的验证实验原理示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”“第四”等(如果存在)是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
如图1所示，本发明实施例中提供的一种星点像亚像元定位精度的验证方法的流程示意图。在该实施例中，星点像弥散斑的获取方法包括六个步骤，每个步骤中具体的内容如下所述。
步骤S1：确定地球自转360度的周期为T，待验证光学***的焦距为f。地球自转360度的周期T可以根据常识获取，T为23小时56分4.09秒，可以近似为24小时。待验证光学***的焦距为f，可以根据所选光学***而确定。
步骤S2：将所述光学***的主光轴对准赤纬为第一预设值的恒星进行积分，积分时间为t0，在t时间间隔内进行n次观测并且在所述观测过程中所述光学***保持静止不动，其中，t大于t0。第一预设值可以用β来表示。在该实施例中，n为大于2的自然数。
步骤S3：采用待验证的亚像元定位算法计算各个星点像的位置，获得各个星点像的位置坐标(x1,y1)、(x2,y2)…(xn,yn)，并采用距离公式计算相邻两个星像点之间的距离。距离公式的表达式如公式所示：
Figure BDA0001878417790000041
其中，lk'为第k次两个相邻星像点之间的距离，(xk,yk)为第k个星点像的位置坐标。
步骤S4：求取步骤S2中的第1次观测时所述恒星的入射光线与所述光线***的主光轴的夹角。在该实施例中，所述恒星的入射光线与所述光线***的主光轴的夹角的求取表达式如公式2所示：
Figure BDA0001878417790000042
其中，α为所述恒星的入射光线与所述光线***的主光轴的夹角，d为星像点到相机焦点的距离，f为待验证光学***的焦距。
步骤S5：计算所述两次相邻观测时间点中地球转过的角度和所述两次相邻观测的视轴夹角，根据所述地球转过的角度和所述两次相邻观测的视轴夹角计算相邻两次观测的星像点在像面上的移动距离。
在该实施例中，计算两次相邻观测时间点中地球转过的角度的表达式如公式3所示：
Δ＝360·t/T (公式3)
其中，t为进行n次观测的时间间隔长度，T为地球自转360度的周期。
在该实施例中，计算两次相邻观测的视轴夹角的计算公式如公式4所示：
Figure BDA0001878417790000051
其中，Δ为两次相邻观测时间点中地球转过的角度，β为步骤S2中的所述第一预设值。
在该实施例中，计算相邻两次观测的星像点在像面上的移动距离表达式如公式5所示：
lk＝f·(tan(α+kx)-tan(α+(k-1)x),...k∈[1..n-1] (公式5)
其中，f为待验证光学***的焦距，α为所述恒星的入射光线与所述光线***的主光轴的夹角，x为两次相邻观测的视轴夹角。
步骤S6：根据所述相邻两次观测的星像点在像面上的移动距离，计算所述待验证的亚像元定位算法的定位精度。在该实施例中，定位精度包括误差均值V和方差D。
在该实施例中，误差均值的计算表达式如公式6所示：
Figure BDA0001878417790000052
在该实施例中，方差的计算表达式如公式7所示：
Figure BDA0001878417790000061
其中，lk'第k次两个相邻星像点之间的距离，lk为第k次两个相邻星像点在像面上的移动距离。
待验证的亚像元定位算法包括很多种算法，此处不再一一举例。在本发明实施例中，以距心定位算法作为具体的待验证的亚像元定位算法并以天狼星为观测恒星作为具体实施例，来阐述本发明所提供的星点像亚像元定位精度的验证过程。
结合图2和图3所述，天狼星星点像弥散斑获取及其验证的过程为：
确定光学仪器的焦距f为60mm，探测器像元大小为5.5μm，探测器像元个数为4288×2848，观测目标为天狼星。天狼星的天球坐标为赤纬角度为-16°42’47.315”，赤经角度为06h45m9.2499s。采用积分时间1s，间隔时间2s内连续拍摄8张天狼星星点像，如图2所示。图2中(a)图表示0秒图像，(b)图表示2秒图像，(c)图表示4秒图像，(d)图表示6秒图像。采用距心算法求解8个星点像的位置如下：T1(1648.207,1474.284)、T2(1648.967,1475.665)、T3(1649.613,1476.900)、T4(1650.350,1478.384)、T5(1651.084,1479.645)、T6(1651.870,1480.991)、T7(1652.543,1482.336)、T8(1653.259,1483.734)。根据公式1计算获得求得各点之间的依次距离为：T21为1.579像元，T32为1.397像元，T43为1.659像元，T54为1.463像元，T65为1.562像元，T76为1.507像元，T87为1.574像元。
光学仪器的焦点位于(2144，1424)像元，则通过公式2求得第1次观测时恒星入射光线与光学***主光轴的夹角为α为2.616°。
根据公式3，计算在两次相邻观测时间点中地球转过角度Δ＝360·2/86164＝0.008356°。根据公式4，计算获得两次观测的视轴夹角为x＝0.008003°，如图3所示。根据公式5，计算获得相邻两次观测的像点在像面上移动的距离lk＝f·(tan(α+kx)-tan(α+(k-1)x),...k∈[1..n-1]。则该定位方法的误差均值V和方差D可分别由公式6和公式7求得：
Figure BDA0001878417790000062
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (9)

1.一种星点像亚像元定位精度的验证方法，其特征在于，所述获取方法包括：
步骤S1：确定地球自转360度的周期为T，待验证光学***的焦距为f；
步骤S2：将所述光学***的主光轴对准赤纬为第一预设值的恒星进行积分，积分时间为t0，在t时间间隔内进行n次观测并且在所述观测过程中所述光学***保持静止不动，其中，t大于t0；
步骤S3：采用待验证的亚像元定位算法计算各个星点像的位置，获得各个星点像的位置坐标，并采用距离公式计算相邻两个星像点之间的距离；
步骤S4：求取步骤S2中的第1次观测时所述恒星的入射光线与所述光线***的主光轴的夹角；
步骤S5：计算所述两次相邻观测时间点中地球转过的角度和所述两次相邻观测的视轴夹角，根据所述地球转过的角度和所述两次相邻观测的视轴夹角计算相邻两次观测的星像点在像面上的移动距离；
步骤S6：根据所述相邻两次观测的星像点在像面上的移动距离，计算所述待验证的亚像元定位算法的定位精度。
2.根据权利要求1所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述步骤S3中的距离公式的表达式为
Figure FDA0001878417780000011
其中，lk'为第k次两个相邻星像点之间的距离，(xk,yk)为第k个星点像的位置坐标。
3.根据权利要求1所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述步骤S4中的求取公式为
Figure FDA0001878417780000012
其中，α为所述恒星的入射光线与所述光线***的主光轴的夹角，d为星像点到相机焦点的距离，f为待验证光学***的焦距。
4.根据权利要求1所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述步骤S5中所述两次相邻观测时间点中地球转过的角度的计算公式为Δ＝360·t/T，其中，t为进行n次观测的时间间隔长度，T为地球自转360度的周期。
5.根据权利要求4所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述步骤S5中所述两次相邻观测的视轴夹角的计算公式为
Figure FDA0001878417780000021
其中，Δ为两次相邻观测时间点中地球转过的角度，β为步骤S2中的所述第一预设值。
6.根据权利要求5所述的一种星点像亚像元定位精度的验证方法，其特征在于，步骤S5中的所述移动距离的计算公式为
lk＝f·(tan(α+kx)-tan(α+(k-1)x),...k∈[1..n-1]，
其中，f为待验证光学***的焦距，α为所述恒星的入射光线与所述光线***的主光轴的夹角，x为两次相邻观测的视轴夹角。
7.根据权利要求1所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述定位精度包括误差均值和方差。
8.根据权利要求7所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述误差均值的计算公式为：
Figure FDA0001878417780000022
所述方差的计算公式为：
Figure FDA0001878417780000023
其中，lk'为第k次两个相邻星像点之间的距离，lk为第k次两个相邻星像点在像面上的移动距离。
9.根据权利要求1所述的一种星点像亚像元定位精度的验证方法，其特征在于，所述n为大于2的自然数。
CN201811410797.7A 2018-11-24 2018-11-24 一种星点像亚像元定位精度的验证方法 Active CN111220177B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811410797.7A CN111220177B (zh) 2018-11-24 2018-11-24 一种星点像亚像元定位精度的验证方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811410797.7A CN111220177B (zh) 2018-11-24 2018-11-24 一种星点像亚像元定位精度的验证方法

Publications (2)

Publication Number Publication Date
CN111220177A true CN111220177A (zh) 2020-06-02
CN111220177B CN111220177B (zh) 2022-07-22

Family

ID=70830449

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811410797.7A Active CN111220177B (zh) 2018-11-24 2018-11-24 一种星点像亚像元定位精度的验证方法

Country Status (1)

Country Link
CN (1) CN111220177B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5562266A (en) * 1992-10-29 1996-10-08 Aerospatiale Societe Nationale Industrielle Rate gyro calibration method and apparatus for a three-axis stabilized satellite
CN101344391A (zh) * 2008-07-18 2009-01-14 北京工业大学 基于全功能太阳罗盘的月球车位姿自主确定方法
CN101441082A (zh) * 2008-12-22 2009-05-27 中国科学院长春光学精密机械与物理研究所 快速三角形星图识别方法
CN102607597A (zh) * 2012-03-08 2012-07-25 清华大学 星敏感器的三轴精度表述与测量方法
EP2566149A1 (en) * 2010-04-28 2013-03-06 Pentax Ricoh Imaging Company, Ltd. Automatic celestial-body tracking / image-capturing method and automatic celestial-body tracking / image-capturing camera
CN105371870A (zh) * 2015-11-27 2016-03-02 北京空间飞行器总体设计部 一种基于星图数据的星敏感器在轨精度测量方法
CN105444778A (zh) * 2015-11-10 2016-03-30 北京空间飞行器总体设计部 一种基于成像几何反演的星敏感器在轨定姿误差获取方法
CN106767714A (zh) * 2016-11-29 2017-05-31 上海卫星工程研究所 提高卫星图像定位精度的等效失配模型多级标校方法
CN107329190A (zh) * 2017-05-17 2017-11-07 上海卫星工程研究所 静止气象卫星恒星敏感的成像试验方法
CN107607127A (zh) * 2017-07-28 2018-01-19 湖北航天技术研究院总体设计所 一种基于外场的星敏感器内部参数标定及精度快速验证***

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5562266A (en) * 1992-10-29 1996-10-08 Aerospatiale Societe Nationale Industrielle Rate gyro calibration method and apparatus for a three-axis stabilized satellite
CN101344391A (zh) * 2008-07-18 2009-01-14 北京工业大学 基于全功能太阳罗盘的月球车位姿自主确定方法
CN101441082A (zh) * 2008-12-22 2009-05-27 中国科学院长春光学精密机械与物理研究所 快速三角形星图识别方法
EP2566149A1 (en) * 2010-04-28 2013-03-06 Pentax Ricoh Imaging Company, Ltd. Automatic celestial-body tracking / image-capturing method and automatic celestial-body tracking / image-capturing camera
CN102607597A (zh) * 2012-03-08 2012-07-25 清华大学 星敏感器的三轴精度表述与测量方法
CN105444778A (zh) * 2015-11-10 2016-03-30 北京空间飞行器总体设计部 一种基于成像几何反演的星敏感器在轨定姿误差获取方法
CN105371870A (zh) * 2015-11-27 2016-03-02 北京空间飞行器总体设计部 一种基于星图数据的星敏感器在轨精度测量方法
CN106767714A (zh) * 2016-11-29 2017-05-31 上海卫星工程研究所 提高卫星图像定位精度的等效失配模型多级标校方法
CN107329190A (zh) * 2017-05-17 2017-11-07 上海卫星工程研究所 静止气象卫星恒星敏感的成像试验方法
CN107607127A (zh) * 2017-07-28 2018-01-19 湖北航天技术研究院总体设计所 一种基于外场的星敏感器内部参数标定及精度快速验证***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
RUSLAN BELIKOV 等: "High Contrast Imaging in Multi-Star Systems: Progress in Technology Development and Lab Results", 《TECHNIQUES AND INSTRUMENTATION FOR DETECTION OF EXOPLANETS VIII》 *
李春艳 等: "高精度星敏感器天文标定方法及观测分析", 《光学精密工程》 *
马岩 等: "基于CCD成像单元仿真测试***的实时性实现", 《液晶与显示》 *

Also Published As

Publication number Publication date
CN111220177B (zh) 2022-07-22

Similar Documents

Publication Publication Date Title
US10366511B2 (en) Method and system for image georegistration
CA2981382C (en) Method of determining a direction of an object on the basis of an image of the object
EP3149698B1 (en) Method and system for image georegistration
US7751651B2 (en) Processing architecture for automatic image registration
CN102741706B (zh) 地理参照图像区域的方法
KR100762891B1 (ko) Ｌｏｓ벡터 조정모델을 이용한 영상의 기하보정 방법 및 그장치
US10621741B2 (en) System and method for a digital image sextant
US20110004405A1 (en) Earth horizon sensor
EP2513664A1 (fr) Procede de calibrage d'un instrument de mesure d'un systeme optronique
US20070117078A1 (en) Celestial compass
EP2022007A2 (en) System and architecture for automatic image registration
US20140085717A1 (en) Systems and methods for closed-loop telescope control
CN104361563B (zh) 基于gps的高光谱遥感图像几何精校正方法
Schwind et al. An in-depth simulation of EnMAP acquisition geometry
CN113029132A (zh) 一种结合地面影像与恒星光行差测量的航天器导航方法
Speyerer et al. In-flight geometric calibration of the lunar reconnaissance orbiter camera
Cohenour et al. Camera models for the wright patterson air force base (WPAFB) 2009 wide-area motion imagery (WAMI) data set
CN111220177B (zh) 一种星点像亚像元定位精度的验证方法
US11847157B2 (en) Telescope star searching method and device based on image recognition and telescope
Pellizza et al. An experiment to measure the instantaneous distance to the Moon
CN111220178B (zh) 一种遥感器光轴指向精度在轨矫正方法
Hrabar et al. PTZ camera pose estimation by tracking a 3D target
Amran et al. Assessment on different digital camera calibration software for photogrammetric applications
Pavlovic et al. The first test of new Andor iXon 897 EMCCD camera
Wasson Measuring Double Stars with a Dobsonian Telescope by the Video Drift Method

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