CN106910230B - 一种雷达背景视频绘制方法 - Google Patents

一种雷达背景视频绘制方法 Download PDF

Info

Publication number
CN106910230B
CN106910230B CN201611186529.2A CN201611186529A CN106910230B CN 106910230 B CN106910230 B CN 106910230B CN 201611186529 A CN201611186529 A CN 201611186529A CN 106910230 B CN106910230 B CN 106910230B
Authority
CN
China
Prior art keywords
radar
radar background
distance line
coordinates
background
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.)
Active
Application number
CN201611186529.2A
Other languages
English (en)
Other versions
CN106910230A (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.)
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Original Assignee
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
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 Leihua Electronic Technology Research Institute Aviation Industry Corp of China filed Critical Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Priority to CN201611186529.2A priority Critical patent/CN106910230B/zh
Publication of CN106910230A publication Critical patent/CN106910230A/zh
Application granted granted Critical
Publication of CN106910230B publication Critical patent/CN106910230B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/203Drawing of straight lines or curves

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及雷达背景视频绘制工程应用技术领域,特别涉及一种雷达背景视频绘制方法。绘制方法包括如下步骤:雷达图像处理单元接收第一雷达背景距离线;计算出在屏幕上第一雷达背景距离线末端的像素点坐标;获取第二雷达背景距离线末端的像素点坐标;利用中点画圆法计算出相邻两条雷达背景距离线之间末端圆弧的预定数量的像素点坐标;利用中点画线法完成两条雷达背景距离线之间背景扇面;完成第二雷达背景距离线与第三雷达背景距离线之间的背景扇面;重复上述步骤,直到完成雷达背景视频绘制更新。本发明的雷达背景视频绘制方法,能够取代传统的雷达背景视频绘制方法中像素点坐标的繁琐计算,提高了雷达背景视频的质量以及运算效率。

Description

一种雷达背景视频绘制方法
技术领域
本发明涉及雷达背景视频绘制工程应用技术领域,特别涉及一种雷达背景视频绘制方法。
背景技术
雷达背景视频是由若干条包含角度值和像素值的距离线通过坐标转换投影至相应的雷达坐标系下,在屏幕上积累形成的。其中,极坐标是较为常用的雷达坐标系。通常,通过距离线角度的正弦和余弦计算距离线在屏幕对应的像素坐标,完成距离线的投影过程。
由于雷达背景距离线之间存在角度差,如果将若干距离线直接投影至综合任务显示器上时,距离线之间会出现较为明显的间隔,显示质量很差。为了解决这一问题,传统的方法是在相邻两条距离线之间插入数条相同的距离线,已达到填充效果。
由于传统雷达背景视频显示区域小,背景距离线更新速度低,故画面显示对处理器计算能力要求不高,通常不存在处理器计算瓶颈问题。但随着高分辨率雷达的逐渐应用,对雷达背景视频显示质量要求越来越高,背景视频显示区域也越来越来大。以传统投影方法和填充方法,处理器的运算能力无法满足实时更新高质量的雷达背景视频。如何利用有限的硬件资源完成高质量的雷达背景视频实时更新成了工程应用上的一大难点和热点。
中点画圆和中点画线算法是目前图像处理领域中常用的快速画圆和画线算法。将这两种算法应用于雷达背景视频的绘制目前还鲜有相关研究应用。但是,传统中点画线法画出的直线的像素值是一致的,而雷达背景距离线上每个像素点的像素值是根据不同距离门上的回波幅值确定的,因而直接使用不适用于雷达距离线的绘制。
发明内容
本发明的目的是提供了一种雷达背景视频绘制方法,以至少解决目前雷达背景视频绘制方法中出现的雷达背景视频显示质量差、运算效率低的问题。
本发明的技术方案是:
一种雷达背景视频绘制方法,包括:
步骤一、雷达图像处理单元接收雷达信号处理单元发送的一条第一雷达背景距离线;
步骤二、以屏幕左上角为原点,根据所述第一雷达背景距离线角度信息以及背景视频扇区圆心坐标,计算出在屏幕上所述第一雷达背景距离线末端(即扇形末端)的像素点坐标;
步骤三、重复步骤一和步骤二,获取相邻的一条第二雷达背景距离线末端(即扇形末端)的像素点坐标;
步骤四、根据相邻两条所述雷达背景距离线的末端坐标以及所述背景视频扇区圆心坐标,利用中点画圆法计算出所述相邻两条雷达背景距离线之间末端圆弧的预定数量的像素点坐标;
步骤五、利用中点画线法将所述末端圆弧上预定数量的像素点坐标与所述背景视频扇区圆心坐标直接连线,形成预定数量的预定雷达背景距离线,完成所述第一雷达背景距离线与所述第二雷达背景距离线之间背景扇面,并且,所述预定数量的预定雷达背景距离线的像素值与所述第二雷达背景距离线上的像素值等比例对应;
步骤六、重复步骤三,获取与所述第二雷达背景距离线相邻的一条第三雷达背景距离线末端的像素点坐标;再重复步骤四和步骤五,完成所述第二雷达背景距离线与所述第三雷达背景距离线之间的背景扇面;
步骤七、重复步骤六,完成雷达背景视频绘制更新(是一个持续不断的绘制过程)。
可选的,在所述步骤四中,所述中点画圆法中是以屏幕左上角为原点。
可选的,在所述步骤五中,所述中点画线法包括:
当|k|≤1时,以所述屏幕的X轴坐标进行等比例抽取,并根据如下关系式(1)确定所述预定雷达背景距离线上点的像素值:
Figure GDA0002381000110000031
其中,k为所述第二雷达背景距离线的斜率,xn为所述第二雷达背景距离线上某个点的横轴坐标值,CentX是所述背景视频扇区圆心坐标,Len为背景距离线的长度(即像素点数),BUFFDATA为所述第二雷达背景距离线的像素值,
Figure GDA0002381000110000032
取整;
当|k|>1时,以所述屏幕的Y轴坐标进行等比例抽取,并根据如下关系式(2)确定所述预定雷达背景距离线上点的像素值;以Y轴坐标进行等比例抽取,假设yn为距离线上某个点的纵轴坐标值,则此点对的像素值:
Figure GDA0002381000110000033
其中,Yn为所述第二雷达背景距离线上某个点的纵轴坐标值,CentY是所述背景视频扇区圆心坐标,
Figure GDA0002381000110000034
取整。
发明效果:
本发明的雷达背景视频绘制方法,利用中点画圆法确定待插入的距离线,然后利用中点画线法,完成距离线以及待插入的距离线的绘制,取代传统的雷达背景视频绘制方法中像素点坐标的繁琐计算,一方面是提高了雷达背景视频的质量,另外,还提高运算效率、缩短计算时间,能够更好地应用于工程实践。
附图说明
图1是本发明雷达背景视频绘制方法中雷达背景视频绘制总体方案图;
图2是本发明基于改进中点画圆和画线法雷达背景视频增强的绘制流程图;
图3是本发明雷达背景视频绘制方法中点画圆算法示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明保护范围的限制。
下面结合附图1至图3对本发明雷达背景视频绘制方法做进一步详细说明。
由于雷达背景距离线之间存在角度差,如果不做处理,直接投影至综合任务显示器上时,距离线之间会出现较为明显的间隔,显示质量很差。
为此,本发明提供了一种雷达背景视频绘制方法,包括如下步骤:
步骤一、雷达图像处理单元接收雷达信号处理单元发送的一条第一雷达背景距离线,图像处理单元在相同周期内完成对背景距离线的投影;通常,雷达信号处理单元会以固定周期进行发送;本实施例中,优选的是相对现有固定周期的更快速周期。
步骤二、以屏幕左上角为原点(参见图3),根据所述第一雷达背景距离线角度信息以及背景视频扇区圆心坐标,计算出在屏幕上第一雷达背景距离线末端(即扇形末端)的像素点坐标。
步骤三、重复步骤一和步骤二,获取相邻的一条第二雷达背景距离线末端(即扇形末端)的像素点坐标。
步骤四、根据相邻两条雷达背景距离线的末端坐标以及背景视频扇区圆心坐标,利用中点画圆法计算出相邻两条雷达背景距离线之间末端圆弧的预定数量(具体数量是根据实际条件不同而不同)的像素点坐标(参见图3)。并且,此步骤中优选中点画圆法中是以屏幕左上角为原点。
步骤五、利用中点画线法将末端圆弧上预定数量的像素点坐标与背景视频扇区圆心坐标直接连线,形成预定数量的预定雷达背景距离线,完成第一雷达背景距离线与第二雷达背景距离线之间背景扇面,并且,预定数量的预定雷达背景距离线的像素值与第二雷达背景距离线上的像素值等比例对应(即连线上每个像素点的像素值通过等比例抽距离线像素值取确定)。
具体地,以屏幕左上角为原点,中点画圆算法推导解算如下:
如附图3所示,假设P1,P2分别为相邻的两条距离线的末端坐标,令P1=(x0,y0),P2=(x1,y1),距离线方位角分别为θ0,θ1。设雷达背景视频扇区圆点坐标为CentP(CentX,CentY),距离线长度为Len,则
Figure GDA0002381000110000061
Figure GDA0002381000110000062
构造判别式:
F(X,Y)=(X-CentX)2+(Y-CentY)2-Len2......(5);
默认F(x0,y0)=0,F(x1,y1)=0。
如附图3所示,确定第一个点P1坐标后,下一个坐标必然在n1和n2中选择。通过两点之间的中间点m坐标代入判别式后,通过结果的正负选择n1或者n2;具体如下:
d0=F(x0+1,y0-0.5)=(x0+1-CentX)2+(y0-0.5-CentY)2-Len2...(6);
若d0<0,则选择n1(x0+1,y0)为下一个坐标点,那么下一个坐标点的判别式应为:
d1=F(x0+2,y0-0.5)=(x0+2-CentX)2+(y0-0.5-CentY)2-Len2
=(x0+1-CentX)2+(y0-0.5-CentY)2-Len2+2*(x0-CentX)+3
=d0+2*(x0-CentX)+3...(7);
若d0>=0,则选择n2(x0+1,y0+1)为下一个坐标点,那么下一个坐标点的判别式应为:
d1=F(x0+2,y0-1.5)=(x0+2-CentX)2+(y0-1.5-CentY)2-Len2
=(x0+1-CentX)2+(y0-0.5-CentY)2-Len2+2*(x0-y0-CentX+CentY)+2
=d0+2*(x0-y0-CentX+CentY)+2...(8);
通过式(5)和式(6)可以分析出:
当d<0时,dn+1=dn+2*(xn-CentX)+3......(9);
当d>=0时,dn+1=dn+2*(xn-yn-CentX+CentY)+2......(10);
最终由式(9)和式(10)可以将P1和P2两点之间的坐标逐一确定。
传统的中点画线法只适用于在屏幕中画像素值一致的直线,无法用于雷达背景视频的绘制。本发明在中点画线法过程中,加入等比例抽取距离线像素值获取距离线上每个点的像素值,应用于雷达背景视频绘制。具体方法如下:
当|k|≤1时,以所述屏幕的X轴坐标进行等比例抽取,并根据如下关系式(1)确定所述预定雷达背景距离线上点的像素值:
Figure GDA0002381000110000071
其中,k为所述第二雷达背景距离线的斜率,xn为所述第二雷达背景距离线上某个点的横轴坐标值,CentX是所述背景视频扇区圆心坐标,Len为背景距离线的长度,BUFFDATA为所述第二雷达背景距离线的像素值,
Figure GDA0002381000110000072
取整;
当|k|>1时,以所述屏幕的Y轴坐标进行等比例抽取,并根据如下关系式(2)确定所述预定雷达背景距离线上点的像素值;以Y轴坐标进行等比例抽取,假设yn为距离线上某个点的纵轴坐标值,则此点对的像素值:
Figure GDA0002381000110000081
其中,yn为所述第二雷达背景距离线上某个点的纵轴坐标值,CentY是所述背景视频扇区圆心坐标,
Figure GDA0002381000110000082
取整。
具体地,以屏幕分辨率为1920*1080,屏幕原点为左上角,距离线角度为30°,长度为900为例;设背景视频扇区圆心坐标为(960,1000);则距离线末端坐标为:
x0=cos(30°)*900+960=1739;
y0=1000–sin(30°)=550;
此时距离线斜率k=0.557,以X轴坐标进行等比例抽取,当xn=1500时,
Figure GDA0002381000110000083
步骤六、重复步骤三,获取与第二雷达背景距离线相邻的一条第三雷达背景距离线末端的像素点坐标;再重复步骤四和步骤五,完成第二雷达背景距离线与所述第三雷达背景距离线之间的背景扇面。
步骤七、重复步骤六,完成雷达背景视频绘制更新;相当于是一个持续不断的绘制过程。
综上所述,本发明的雷达背景视频绘制方法技术方案可以归纳如下:
1)、如图1所示,图像处理单元以快速的周期接收雷达背景距离线,并以相同周期完成距离线投影以及相邻距离线之间的缝隙填充。
2)、如图2所示,在图像处理时,以屏幕左上角为原点,根据背景距离线角度信息以及背景视频扇区圆心坐标计算出,在屏幕上此距离线末端的像素点坐标。然后根据相邻的距离线的末端坐标以及雷达背景视频的扇区圆心,利用改进的中点画圆法,计算出相邻距离线之间圆弧的各个像素点坐标。这些坐标即为待插入的各个距离线的末端坐标。
3)、如图3所示,利用中点画线法将距离线末端坐标点与雷达背景视频的扇区圆心直接连线,通过等比例抽取法获取待画直线上各个像素点的像素值,需要注意的是,此过程中待画的距离线包含了本次接收到的距离线。
本发明的雷达背景视频绘制方法,利用中点画圆法确定待插入的距离线,然后利用中点画线法,完成距离线以及待插入的距离线的绘制,取代传统的雷达背景视频绘制方法中像素点坐标的繁琐计算,一方面是提高了雷达背景视频的质量,另外,还提高运算效率、缩短计算时间,能够更好地应用于工程实践。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (1)

1.一种雷达背景视频绘制方法,其特征在于,包括如下步骤:
步骤一、雷达图像处理单元接收雷达信号处理单元发送的一条第一雷达背景距离线;
步骤二、以屏幕左上角为原点,根据所述第一雷达背景距离线角度信息以及背景视频扇区圆心坐标,计算出在屏幕上所述第一雷达背景距离线末端的像素点坐标;
步骤三、重复步骤一和步骤二,获取相邻的一条第二雷达背景距离线末端的像素点坐标;
步骤四、根据相邻两条所述雷达背景距离线的末端坐标以及所述背景视频扇区圆心坐标,利用中点画圆法计算出相邻两条雷达背景距离线之间末端圆弧的预定数量的像素点坐标;
步骤五、利用中点画线法将所述末端圆弧上预定数量的像素点坐标与所述背景视频扇区圆心坐标直接连线,形成预定数量的预定雷达背景距离线,完成所述第一雷达背景距离线与所述第二雷达背景距离线之间背景扇面,并且,所述预定数量的预定雷达背景距离线的像素值与所述第二雷达背景距离线上的像素值等比例对应;
步骤六、重复步骤三,获取与所述第二雷达背景距离线相邻的一条第三雷达背景距离线末端的像素点坐标;再重复步骤四和步骤五,完成所述第二雷达背景距离线与所述第三雷达背景距离线之间的背景扇面;
步骤七、重复步骤六,完成雷达背景视频绘制更新;
所述中点画圆法中是以屏幕左上角为原点;
所述中点画线法包括:
当|k|≤1时,以所述屏幕的X轴坐标进行等比例抽取,并根据如下关系式(1)确定所述预定雷达背景距离线上点的像素值:
Figure FDA0002381000100000021
其中,k为所述第二雷达背景距离线的斜率,xn为所述第二雷达背景距离线上某个点的横轴坐标值,CentX是所述背景视频扇区圆心坐标,Len为背景距离线的长度,BUFFDATA为所述第二雷达背景距离线的像素值,
Figure FDA0002381000100000022
取整;
当|k|>1时,以所述屏幕的Y轴坐标进行等比例抽取,并根据如下关系式(2)确定所述预定雷达背景距离线上点的像素值;以Y轴坐标进行等比例抽取,假设yn为距离线上某个点的纵轴坐标值,则此点的像素值:
Figure FDA0002381000100000023
其中,CentY是所述背景视频扇区圆心坐标,
Figure FDA0002381000100000024
取整。
CN201611186529.2A 2016-12-21 2016-12-21 一种雷达背景视频绘制方法 Active CN106910230B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611186529.2A CN106910230B (zh) 2016-12-21 2016-12-21 一种雷达背景视频绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611186529.2A CN106910230B (zh) 2016-12-21 2016-12-21 一种雷达背景视频绘制方法

Publications (2)

Publication Number Publication Date
CN106910230A CN106910230A (zh) 2017-06-30
CN106910230B true CN106910230B (zh) 2020-07-03

Family

ID=59207274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611186529.2A Active CN106910230B (zh) 2016-12-21 2016-12-21 一种雷达背景视频绘制方法

Country Status (1)

Country Link
CN (1) CN106910230B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109085553B (zh) * 2018-08-08 2023-03-03 中国航空工业集团公司雷华电子技术研究所 一种基于纹理的快速扫描雷达回波背景视频绘制方法
CN110415316B (zh) * 2019-07-02 2023-05-05 中国航空工业集团公司雷华电子技术研究所 一种提高气象视频显示准确性和连续性的插值方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104217450A (zh) * 2014-09-05 2014-12-17 南京航空航天大学 一种任意角度圆弧的快速绘制方法
CN105372657A (zh) * 2015-12-10 2016-03-02 中国科学院电子学研究所 基于回波数据的视频合成孔径雷达运动补偿成像方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010060078A2 (en) * 2008-11-24 2010-05-27 Bossdev, Inc. Background signal reduction in neutron fluorescence applications using agile neutron beam flux

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104217450A (zh) * 2014-09-05 2014-12-17 南京航空航天大学 一种任意角度圆弧的快速绘制方法
CN105372657A (zh) * 2015-12-10 2016-03-02 中国科学院电子学研究所 基于回波数据的视频合成孔径雷达运动补偿成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《Detecting background and foreground with a laser array system》;Teresa Hernández-Díaz;《Measurement》;20141212;第63卷;第195-206页 *
《电子海图/导航雷达信息叠加技术研究》;李刚;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20061215(第12期);正文第1-65页 *
《背景提取算法在雷达视频图像中的应用》;熊伟;《舰船电子工程》;20131231;第33卷(第12期);第40-41、66页 *
《雷达视频图像采集和处理技术研究》;胡珺;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110915(第9期);正文第1-75页 *

Also Published As

Publication number Publication date
CN106910230A (zh) 2017-06-30

Similar Documents

Publication Publication Date Title
CN106780590B (zh) 一种深度图的获取方法及系统
CN102968809B (zh) 在增强现实领域实现虚拟信息标注及绘制标注线的方法
CN103077509B (zh) 利用离散立方体全景图实时合成连续平滑全景视频的方法
CN103489214A (zh) 增强现实系统中基于虚拟模型预处理的虚实遮挡处理方法
CN106408513B (zh) 深度图超分辨率重建方法
CN103810685A (zh) 一种深度图的超分辨率处理方法
CN109115184B (zh) 基于非合作目标协同测量方法及系统
CN103826068A (zh) 一种视频图像旋转方法
CN106910230B (zh) 一种雷达背景视频绘制方法
WO2019056782A1 (zh) 一种基于球体投影公切线的多相机标定及参数优化方法
CN104881866A (zh) 一种扩展小孔成像模型的鱼眼相机矫正与标定方法
CN110096993A (zh) 双目立体视觉的目标检测设备及方法
WO2018201664A1 (zh) 一种立体图形显示的方法、装置及设备
CN101188020A (zh) 投影仪投影图象与计算机帧缓存图象之间象素几何位置对应关系精确获取方法
CN111105452A (zh) 基于双目视觉的高低分辨率融合立体匹配方法
WO2024017003A1 (zh) 基于算法融合的车辆汇入检测方法、装置及设备
Zhang et al. Real-Time object detection for 360-degree panoramic image using CNN
CN114332880A (zh) 一种文本检测方法、装置、设备及存储介质
CN110651475B (zh) 用于致密光学流的阶层式数据组织
US11055818B2 (en) Panorama image alignment method and device based on multipath images, and portable terminal
CN112634366B (zh) 位置信息的生成方法、相关装置及计算机程序产品
CN102932661A (zh) 视差图的中值滤波匹配误差校正方法和实施该方法的电路
CN114266879A (zh) 三维数据增强、模型训练检测方法、设备及自动驾驶车辆
WO2018233318A1 (zh) 一种基于线面空间关系的迷宫求解方法
CN115294277B (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