CN111157969A - 一种多传感器系统误差在线估计和数据融合方法 - Google Patents
一种多传感器系统误差在线估计和数据融合方法 Download PDFInfo
- Publication number
- CN111157969A CN111157969A CN202010062265.XA CN202010062265A CN111157969A CN 111157969 A CN111157969 A CN 111157969A CN 202010062265 A CN202010062265 A CN 202010062265A CN 111157969 A CN111157969 A CN 111157969A
- Authority
- CN
- China
- Prior art keywords
- track
- radar
- primary
- secondary radar
- sensor
- 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
Images
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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种多传感器系统误差在线估计和数据融合方法。包括关联过程和保持过程:关联过程中,根据一次雷达和二次雷达的接收到的点迹实时进行航迹关联;进行一次雷达传感器和二次雷达传感器航迹的融合:针对二次雷达传感器的某条航迹,依次求一次雷达传感器的每条航迹与此二次雷达传感器航迹的距离差值序列;获取一次雷达传感器各条航迹对应的距离差值序列的均值和方差,在均值小于系统误差最大值时,找到方差的最小值所对应一次雷达航迹号m0;保持过程中,求二次雷达航迹与匹配的一次雷达第m0条航迹的误差,根据误差求解脱离质量,进而判断二次雷达航迹与匹配的一次雷达航迹关联是否失败。该方案能实时在线估计不同传感器间的系统误差,提升二次雷达的识别能力。
Description
技术领域
本发明涉及误差处理技术领域,特别是一种多传感器系统误差在线估计和数据融合方法。
背景技术
现有一次雷达在扫描过程中,对不同航迹号的点迹进行牵引,由于目标距离近,因此二次雷达(包括AIS等各种平台的二次雷达)应答机应答的目标不一定是牵引的目标。为了对目标进行正确识别,需要将二次雷达应答目标点迹关联后的航迹与一次雷达航迹进行匹配。将二次雷达的识别信息与一次雷达的航迹信息进行融合。
在多源信息融合处理中,存在系统误差。由于天气、距离、高度等多种因素的影响,系统误差会随着时间、测量高度、测量方位、测量距离等因素出现波动。若不能准确实时估计系统误差,导致后续航迹的融合和识别性能急剧下降。给雷达操作员/空管指挥员带来误导,严重时出现误判,威胁飞行安全。
现代多传感器系统误差配准过程,是融合处理中心在将多部雷达上报的目标局部航迹关联的基础上,将源于同一公共目标的局部目标航迹送入系统误差估计模块,系统误差估计模块基于多条样本目标航迹进行误差均值估计,得到量测误差估计值;最后将结果送至系统误差修正模块,对多部传感器送入的目标航迹进行补偿。
现有方法得到的系统误差需要先测得系统误差,结果存于系统误差修正模块。系统误差随着空间分布如高度、距离、方位的不同而不同。这种方法很难得到全局的系统误差,而使不同位置的目标系统误差估计不准确。另外,系统误差随天气等因素影响,会随时间变化。这种方法得到的系统误差不能实时变化,当天气等因素变化时,会造成较大的误差。
发明内容
本发明所要解决的技术问题是:针对上述存在的问题,提供了一种多传感器系统误差在线估计和数据融合方法。
本发明采用的技术方案如下:一种多传感器系统误差在线估计和数据融合方法,包括:
关联过程:步骤S1,分别根据一次雷达和二次雷达的接收到的点迹实时进行航迹关联;步骤S2,获取满足一次雷达和二次雷达航迹位置参数和时间约束条件的航迹;步骤S3,将满足约束条件的一次雷达的航迹中的点迹进行插值估计得到与二次雷达航迹中点迹时刻点相同的位置参数值,求一次雷达传感器和二次雷达传感器间在相同时刻点的位置参数差值,获得一次雷达每条航迹对应的位置参数差值序列(因为每条航迹包括多个点迹,因此每条航迹对应的是一个差值序列);步骤S4,获取各条航迹对应的位置参数差值序列的均值和方差,在均值小于系统误差最大值时,找到方差的最小值所对应一次雷达航迹号m0;
保持过程:求二次雷达航迹与匹配的一次雷达第m0条航迹的误差,根据误差求解脱离质量,根据脱离质量判断二次雷达航迹与匹配的一次雷达航迹关联是否失败,如果失败则返回到关联过程。
进一步的,所述步骤S1中,实时进行航迹关联后,分别记录航迹号:一次雷达航迹track1{m},m=1,2,…M。二次雷达航迹track2{m'},m'=1,2,…M',其中M为一次雷达航迹总数,M'为二次雷达航迹总数。假设二次雷达第m'条航迹为待融合的航迹。
进一步的,所述步骤S2中,所述获取满足一次雷达和二次雷达航迹位置参数和航迹时间约束条件的点迹的过程包括:
设置一次雷达和二次雷达航迹位置参数阈值σr,设置一次雷达和二次雷达航迹时间阈值σt;
若一次雷达的点迹时间k时刻,第m条航迹号位置参数为rm(k);二次雷达第m'条航迹号距离为rm'(t),如果满足约束
|t-k|≤σt
|rm(k)-rm'(t)|≤σr
则记录下满足约束条件的一次雷达点迹的时间、航迹号和位置参数值,二次雷达点迹的时间和位置参数值。
进一步的,所述步骤S3包括以下过程:
将记录的一次雷达第m条航迹的时间和距离序列记为
记录一次雷达传感器和二次雷达传感器待匹配第m'条航迹中的连续Q'个点迹之间的位置参数差,一次雷达第m条航迹对应的差值序列为[γm(t1),γm(t2),…,γm(tQ')]。
进一步的,所述步骤S3中,设置门限η,若Q'<η,则将该航迹对应的距离差值序列删除。
进一步的,保持过程中,所述根据脱离质量判断二次雷达航迹与匹配的一次雷达航迹关联是否失败的过程包括:
进一步的,所述位置参数指距离、方位、高度中的一种或多种。这里,不同类型的传感器可能存在不同的位置参数,可能是距离、也可以是距离和方位,总而言之可以是上述参数的一种或多种。
与现有技术相比,采用上述技术方案的有益效果为:本发明设置关联过程和保持过程,进行了多传感器航迹匹配关联,并利用误差与均值、方差的偏差关系,对关联关系进行验证,该技术实现简单,效果明显,不需要任何人工干预,自动实时估计系统误差,能自动适应全空域、全方位各种目标的处理需求,为雷达操作员/空管指挥员提供可靠的目标信息,有效保障飞行安全。
附图说明
图1是本发明多传感器系统误差在线估计和数据融合方法流程示意图。
图2是本发明实施例中获得的一次雷达和二次雷达目标航迹点迹的距离值变化示意图。
图3是本发明实施例的方法得到的一次雷达和二次雷达匹配航迹示意图。
图4是本发明实施例的方法得到的系统误差随距离变化情况示意图。
具体实施方式
下面结合附图对本发明做进一步描述。
一种多传感器系统误差在线估计和数据融合方法,包括关联过程和保持过程:
关联过程:
(1)分别根据一次雷达和二次雷达的接收到的点迹实时进行航迹关联(此处关联方法可用极大似然比多目标关联、贝叶斯类多目标关联等方法,各种关联方法均为现有方法);
其中,实时进行航迹关联后,分别记录航迹号,其中,一次雷达航迹track1{m},m=1,2,…M,二次雷达航迹track2{m'},m'=1,2,…M';其中M为一次雷达航迹总数,M'为二次雷达航迹总数。假设二次雷达第m'条航迹为待融合的航迹。
(2)假设二次雷达第m'条航迹为待融合的航迹,开始进行一次雷达传感器和二次雷达传感器航迹的融合(这里,一次雷达和二次雷达均可能对应多个传感器,传感器类型可以是AIS、ADS-B等)。获取满足一次雷达和二次雷达航迹位置参数和航迹时间约束的点迹;这里根据传感器的性质,位置参数可能是距离、方位、高度等信息中的一种或多种,根据传感器提供的距离/方位/高度或其他信息中的一个或多个进融合。
其中,如果两条航迹最新状态间时间间隔太久,速度航向等信息的变化容易导致误关联,在这些情况下进行关联都毫无必要。为了解决误关联的问题,所述获取满足一次雷达和二次雷达航迹距离和航迹时间间隔要求的点迹的过程包括:
设置一次雷达和二次雷达航迹位置参数阈值σr,设置一次雷达和二次雷达航迹时间阈值σt;
若一次雷达的点迹时间k时刻,第m条航迹号,位置参数为rm(k);二次雷达第m'条航迹号距离为rm'(t),如果满足要求
|t-k|≤σt
|rm(k)-rm'(t)|≤σr
则记录下满足约束条件的一次雷达点迹的时间、航迹号和位置参数,二次雷达点迹的时间和位置参数。
(3)将满足约束条件的一次雷达的每条航迹的点迹进行插值估计(此处插值方法可用泰勒插值、拉格朗日插值等方法,各种插值方法均为现有方法),得到与二次雷达点迹时刻相同的t时刻的位置参数值,记录一次雷达传感器和二次雷达传感器待匹配第m'条航迹中的连续Q'个点迹之间的位置参数差序列;
步骤(3)包括以下过程:
将满足位置参数阈值σr和时间阈值σt要求的一次雷达第m条航迹的时间和距离序列记为
按照(1)~(3)的步骤,记录一次雷达传感器和二次雷达传感器待匹配第m'条航迹中的连续Q'个点迹之间的位置参数差,一次雷达第m条航迹对应的差值序列为[γm(t1),γm(t2),…,γm(tQ')],因为m代表是航迹号,因此通过上述计算,一次雷达的每条航迹均获得了对应的一个距离差值序列。
优选地,若Q'太小,则说明此航迹不稳定,应不考虑其与二次雷达航迹匹配,因此设置门限η,若Q'<η,则将该航迹对应的距离差值序列删除。
(4)获取各个航迹号对应的距离差值序列的均值ρm和方差σm,在均值ρm小于系统误差最大值σSyserrMax时,在此条件下,给出评价函数
找到评价函数最小值,即方差σm的最小值所对应一次雷达航迹号m0。
保持过程:
求二次雷达待融合航迹与匹配的一次雷达第m0条航迹的误差,根据误差求解脱离质量,根据脱离质量判断二次雷达航迹与匹配的一次雷达航迹关联是否失败,如果失败则返回到关联过程;
其中一个实施例,选取3000~3900s时间范围的雷达航迹和识别点迹的数据,如图1所示:
关联过程:步骤S1,当前t时刻,此时t的取值在3000~3900s范围类,分别根据一次雷达和二次雷达的接收到的点迹实时进行航迹关联;步骤S2,若满足一次雷达和二次雷达航迹距离和时间约束条件的航迹,则进行一次雷达传感器和二次雷达传感器航迹的融合;步骤S3,将满足约束条件的一次雷达的航迹中的点迹进行插值估计得到与二次雷达航迹中点迹时刻点相同的距离值,如图2是获取的3000~3900s时间范围类一次雷达和二次雷达目标航迹点迹的距离值变化,求出一次雷达传感器和二次雷达传感器间的距离差,获得一次雷达传感器每个航迹号对应的和二次雷达传感器待融合航迹之间的距离差值序列;步骤S4,获取一次雷达各个航迹号对应的距离差值序列的均值ρm和方差σm,在均值小于系统误差最大值时,找到方差的最小值所对应一次雷达航迹号m0;
保持过程:求二次雷达待融合航迹与匹配的一次雷达第m0条航迹的误差,根据误差求解脱离质量,根据脱离质量判断二次雷达航迹与匹配的一次雷达航迹关联是否失败,如果失败则返回到关联过程。
如图3所示为本实施例的方法得到的一次雷达和二次雷达匹配航迹;每求出误差的时候,统计此时的距离,得出如图4所示的系统误差随距离变化情况。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员,在不脱离本发明的精神所做的非实质性改变或改进,都应该属于本发明权利要求保护的范围。
Claims (7)
1.一种多传感器系统误差在线估计和数据融合方法,其特征在于,包括:
关联过程:步骤S1,分别根据一次雷达和二次雷达的接收到的点迹实时进行航迹关联;步骤S2,获取满足一次雷达和二次雷达航迹位置参数和时间约束条件的航迹;步骤S3,将满足约束条件的一次雷达的航迹中的点迹进行插值估计得到与二次雷达航迹中点迹时刻点相同的位置参数值,求一次雷达传感器和二次雷达传感器间在相同时刻点的位置参数差值,获得一次雷达每条航迹对应的位置参数差值序列;步骤S4,获取各条航迹对应的位置参数差值序列的均值和方差,在均值小于系统误差最大值时,找到方差的最小值所对应一次雷达航迹号m0;
保持过程:求二次雷达航迹与匹配的一次雷达第m0条航迹的误差,根据误差求解脱离质量,根据脱离质量判断二次雷达航迹与匹配的一次雷达航迹关联是否失败,如果失败则返回到关联过程。
2.如权利要求1所述的多传感器系统误差在线估计和数据融合方法,其特征在于,所述步骤S1中,实时进行航迹关联后,分别记录每条航迹:一次雷达航迹track1{m},m=1,2,…M。二次雷达航迹track2{m'},m'=1,2,…M',其中M为一次雷达航迹总数,M'为二次雷达航迹总数。
3.如权利要求2所述的多传感器系统误差在线估计和数据融合方法,其特征在于,所述步骤S2中,所述获取满足一次雷达和二次雷达航迹位置参数和航迹时间约束条件的点迹的过程包括:
设置一次雷达和二次雷达航迹位置参数阈值σr,设置一次雷达和二次雷达航迹时间阈值σt;
若一次雷达的点迹时间k时刻,第m条航迹号,位置参数值为rm(k);二次雷达第m'条航迹距离为rm'(t),如果满足约束
|t-k|≤σt
|rm(k)-rm'(t)|≤σr
则记录下满足约束的一次雷达点迹的时间、航迹号和位置参数,二次雷达点迹的时间和位置参数。
5.如权利要求4所述的多传感器系统误差在线估计和数据融合方法,其特征在于,所述步骤S3中,设置门限η,若Q'<η,则将该航迹对应的距离差值序列删除。
7.如权利要求1-6任一项所述的多传感器系统误差在线估计和数据融合方法,其特征在于,所述位置参数指距离、方位、高度中的一种或多种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010062265.XA CN111157969B (zh) | 2020-01-20 | 2020-01-20 | 一种多传感器系统误差在线估计和数据融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010062265.XA CN111157969B (zh) | 2020-01-20 | 2020-01-20 | 一种多传感器系统误差在线估计和数据融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111157969A true CN111157969A (zh) | 2020-05-15 |
CN111157969B CN111157969B (zh) | 2022-01-25 |
Family
ID=70564262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010062265.XA Active CN111157969B (zh) | 2020-01-20 | 2020-01-20 | 一种多传感器系统误差在线估计和数据融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111157969B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111609858A (zh) * | 2020-05-28 | 2020-09-01 | 新石器慧通(北京)科技有限公司 | 定位误差确定方法、轨迹生成方法、装置及交通工具 |
CN111723474A (zh) * | 2020-06-01 | 2020-09-29 | 四川九洲电器集团有限责任公司 | 一种基于最大似然准则的临近目标识别方法 |
CN114397636A (zh) * | 2022-03-25 | 2022-04-26 | 中国气象局气象探测中心 | 一种地基雷达反射率因子均一性评估方法、系统及设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100295719A1 (en) * | 2009-05-19 | 2010-11-25 | Raytheon Company | System, Method, and Software for Performing Dual Hysteresis Target Association |
EP2473867A1 (fr) * | 2009-09-01 | 2012-07-11 | Thales | Traitement de donnees multi-cibles pour radars passifs multi-recepteurs en mode sfn ou mfn |
US20130076558A1 (en) * | 2009-07-09 | 2013-03-28 | Thales | Multi-target data processing for multi-static and multi-channel passive radars |
CN103886196A (zh) * | 2014-03-14 | 2014-06-25 | 中国电子科技集团公司第二十八研究所 | 一种抗大系统误差的航迹关联方法 |
CN105116392A (zh) * | 2015-09-09 | 2015-12-02 | 电子科技大学 | 一种针对ais与主动雷达航迹融合及识别的方法 |
US20160299223A1 (en) * | 2015-04-13 | 2016-10-13 | Raytheon Company | Apparatus and method for processing electronic intelligence (elint) and radar tracking data |
CN106371091A (zh) * | 2016-08-24 | 2017-02-01 | 四川九洲空管科技有限责任公司 | Ads‑b与二次雷达监视信息的数据融合方法及装置 |
CN106526586A (zh) * | 2016-11-03 | 2017-03-22 | 成都能通科技有限公司 | 一种基于岸基情报雷达的数据处理方法 |
-
2020
- 2020-01-20 CN CN202010062265.XA patent/CN111157969B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100295719A1 (en) * | 2009-05-19 | 2010-11-25 | Raytheon Company | System, Method, and Software for Performing Dual Hysteresis Target Association |
US20130076558A1 (en) * | 2009-07-09 | 2013-03-28 | Thales | Multi-target data processing for multi-static and multi-channel passive radars |
EP2473867A1 (fr) * | 2009-09-01 | 2012-07-11 | Thales | Traitement de donnees multi-cibles pour radars passifs multi-recepteurs en mode sfn ou mfn |
CN103886196A (zh) * | 2014-03-14 | 2014-06-25 | 中国电子科技集团公司第二十八研究所 | 一种抗大系统误差的航迹关联方法 |
US20160299223A1 (en) * | 2015-04-13 | 2016-10-13 | Raytheon Company | Apparatus and method for processing electronic intelligence (elint) and radar tracking data |
CN105116392A (zh) * | 2015-09-09 | 2015-12-02 | 电子科技大学 | 一种针对ais与主动雷达航迹融合及识别的方法 |
CN106371091A (zh) * | 2016-08-24 | 2017-02-01 | 四川九洲空管科技有限责任公司 | Ads‑b与二次雷达监视信息的数据融合方法及装置 |
CN106526586A (zh) * | 2016-11-03 | 2017-03-22 | 成都能通科技有限公司 | 一种基于岸基情报雷达的数据处理方法 |
Non-Patent Citations (3)
Title |
---|
TADEUSZ PIETKIEWICZ ET AL.: ""Analysis of fusion primary radar, secondary surveillance radar (IFF) and ESM sensor attribute information under Dezert-Smarandache theory"", 《2017 18TH INTERNATIONAL RADAR SYMPOSIUM (IRS)》 * |
杨万海: "《多传感器数据融合及其应用》", 30 April 2004 * |
郭新民: ""面向多传感器数据融合的雷达显控终端的设计方法"", 《雷达与对抗》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111609858A (zh) * | 2020-05-28 | 2020-09-01 | 新石器慧通(北京)科技有限公司 | 定位误差确定方法、轨迹生成方法、装置及交通工具 |
CN111609858B (zh) * | 2020-05-28 | 2021-10-22 | 新石器慧通(北京)科技有限公司 | 定位误差确定方法、轨迹生成方法、装置及交通工具 |
CN111723474A (zh) * | 2020-06-01 | 2020-09-29 | 四川九洲电器集团有限责任公司 | 一种基于最大似然准则的临近目标识别方法 |
CN111723474B (zh) * | 2020-06-01 | 2022-02-18 | 四川九洲电器集团有限责任公司 | 一种基于最大似然准则的临近目标识别方法 |
CN114397636A (zh) * | 2022-03-25 | 2022-04-26 | 中国气象局气象探测中心 | 一种地基雷达反射率因子均一性评估方法、系统及设备 |
CN114397636B (zh) * | 2022-03-25 | 2022-06-17 | 中国气象局气象探测中心 | 一种地基雷达反射率因子均一性评估方法、系统及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN111157969B (zh) | 2022-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111157969B (zh) | 一种多传感器系统误差在线估计和数据融合方法 | |
US11408741B2 (en) | Self-localization estimation device | |
CN109212521B (zh) | 一种基于前视相机与毫米波雷达融合的目标跟踪方法 | |
US10656246B2 (en) | Misalignment estimation for a vehicle radar system | |
US8838374B2 (en) | Automatic correction of past position errors for location and inventory tracking | |
CN105809126B (zh) | Dsrc与车载传感器融合的智能车辆目标跟踪系统及方法 | |
US7081849B2 (en) | Process for sensor resources management | |
CN112285700A (zh) | 基于激光雷达与毫米波雷达融合的机动目标跟踪方法 | |
CN109752741A (zh) | 车辆定位设备 | |
US11587214B2 (en) | Multipath ghost mitigation in vehicle radar system | |
CN112147651B (zh) | 一种异步多车协同目标状态鲁棒估计方法 | |
CN115856872A (zh) | 一种车辆运动轨迹连续跟踪方法 | |
CN113165661A (zh) | 用于确定车辆的修正轨迹的方法和系统 | |
RU2562616C1 (ru) | Способ получения радиотехнической информации и радиотехнический комплекс для его осуществления | |
CN110673134B (zh) | 一种雷达远近交替测距模式下航迹跟踪优化方法 | |
US8630235B2 (en) | Method and apparatus for the computer-aided classification of measurements for the feature-based positioning of a mobile object in a predetermined region | |
CN112255612A (zh) | 一种雷达航迹随机抖动性评估方法 | |
US7167127B2 (en) | Process for tracking vehicles | |
CN113342057B (zh) | 轨迹的融合方法、装置、无人机探测系统、设备和介质 | |
CN115221921A (zh) | 传感器的感知性能评估方法、车辆以及存储介质 | |
CN113933798A (zh) | 一种基于相似性原理的全局传感器系统误差分区配准算法 | |
CN114945843A (zh) | 用于在雷达系统中调整校正信息的方法 | |
CN111723474B (zh) | 一种基于最大似然准则的临近目标识别方法 | |
CN115184895B (zh) | 一种基于导航雷达的静态目标识别方法及系统 | |
US11466991B2 (en) | Sensor data system for a vehicle |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Yang Ya Inventor after: Cheng Qi Inventor before: Yang Ya |
|
GR01 | Patent grant | ||
GR01 | Patent grant |