CN112781593A - 陆基长波导航/授时系统asf预测数据库快速建库方法 - Google Patents
陆基长波导航/授时系统asf预测数据库快速建库方法 Download PDFInfo
- Publication number
- CN112781593A CN112781593A CN202011636854.0A CN202011636854A CN112781593A CN 112781593 A CN112781593 A CN 112781593A CN 202011636854 A CN202011636854 A CN 202011636854A CN 112781593 A CN112781593 A CN 112781593A
- Authority
- CN
- China
- Prior art keywords
- point
- asf
- path
- ground
- wave
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Abstract
本发明公开了一种陆基长波导航/授时系统ASF预测数据库快速建库方法。首先将区域网格化:针对某发射台站,确定发射台站经纬度信息,并以其为中心,以q角度为间隔形成辐射路径。按照步长Dd等分各路径,获取各等分点的经纬度信息;然后获取路径参数信息:通过长波地理信息系统,基于各等分点的经纬度信息,获取路径参数信息,包括高程h、大地电导率s、相对介电常数er,形成路径参数数据库;最后根据获得的数据,并行计算各点的ASF。本发明基于GPU平台,利用积分方法的迭代过程实现了快速计算发射台站周围ASF修正值,大大缩短了计算时间。
Description
技术领域
本发明属于陆基长波导航/授时技术领域,具体涉及一种陆基长波导航/授时系统ASF预测数据库快速建库方法。
背景技术
以罗兰C系统、e罗兰系统、长河二号系统等为代表的陆基长波导航/授时系统,具有覆盖范围大、全天候、抗干扰能力强、可靠性高、稳定性好等突出优点,是现代PNT体系的重要组成部分。然而,陆基长波导航/授时系统的精度受到诸多因素的影响,其中,影响最大的是附加二次相位因子(Additional Secondary phase Factor,ASF)的修正精度,它可能使系统定位误差达到几海里甚至更大。
通常采用理论预测或实测的方法获取不同位置的ASF修正数据,并制备成以地理网格形式分布的ASF数据库,加载于陆基长波导航/授时接收机中。然而,对于陆基长波导航/授时台所覆盖的广大区域而言,目前并未构建完成高精度的ASF修正数据库,仅在局部区域进行了ASF修正方法的验证。快速制备大面积高精度的ASF修正数据库,尤其是随传播路径环境变化,实现ASF修正数据的快速更新,是该领域的迫切需要。
发明内容
本发明的目的是提供一种陆基长波导航/授时系统ASF预测数据库快速建库方法,实现在较短时间内完成陆基长波导航/授时发射台站周围ASF修正值的计算或更新。
本发明所采用的技术方案是,陆基长波导航/授时系统ASF预测数据库快速建库方法,具体按照以下步骤实施:
步骤1、区域网格化,针对某发射台站,确定发射台站经纬度信息,并以发射台为中心,以θ角度为间隔形成辐射路径,以正北方向为起始方向,顺时针划分角度间隔,按照步长Δd等分各路径,称所有路径上的等分点为区域的网格点,获取各网格点的经纬度信息;
步骤2、获取网格点地理信息,通过长波地理信息系统,基于各网格点的经纬度信息,获取路径参数信息,包括高程h、大地电导率σ、相对介电常数εr,形成路径参数数据库;
步骤3、计算附加二次相位因子ASF,根据步骤2获得的数据,通过积分方法的迭代过程,并行计算出各网格点的附加二次相位因子ASF,并将各网格点的附加二次相位因子ASF按路径存入文件中,一条路径的附加二次相位因子ASF存于一个.dat文件,文件以“发射台名+路径与正北方向的夹角θ”命名,最终得到发射台附近区域的ASF预测数据库。
本发明的特点还在于,
步骤1中θ大于0度,小于1度。
步骤1中步长Δd小于200m。
步骤3具体如下:
步骤3.1、从所述步骤2得到的路径参数数据库中读入各网格点的地理信息存入二维数组中,对高程h进行小波滤波处理,将参数数组从CPU内存复制到GPU内存,以便后续在GPU中并行计算;
步骤3.2、利用积分方法的迭代过程并行计算各点的衰减因子Wg,并存入衰减因子数组中;
步骤3.3、将在GPU计算得到的衰减因子数组复制回CPU中,地波衰减因子的相位又称为二次时延,可得到时延tw;
步骤3.4、利用积分方法计算得到一条海水路径下的衰减因子,海水路径下衰减因子的相位即为二次相位因子SF;
步骤3.5、根据附加二次相位因子ASF=tw-SF得到各网格点的ASF值。
步骤3.2具体如下:
其中
式中,Ez为接受点P(ρ,Φ,z)处z方向的场强,Eb为发射源和接收点P都位于理想导电平面时P点的场强,ω为地波的角频率,μ0为真空磁导率,Idl为垂直电偶极子大小,垂直电偶极子即为发射源,i为虚数单位,k0为空气的波数,kg为地面的波数,z、ρ分别为接收点P在柱坐标系下z方向、ρ方向的坐标值,h为垂直电偶极子Idl的离地高度,pr为发射源与接收点P之间的大圆距离;
从第个点开始,根据各点的地理信息以及计算得到的地面的衰减因子Wg,利用积分法公式(8)~(10)并行计算出所有路径第个点的地面的衰减因子Wg,按照此步骤,依次迭代计算出路径上剩余各点的地面的衰减因子Wg:
式中,fg是接收点处的地波衰减函数,dp为接收点与发射点的大圆距离,r0表示从发射点到接收点的直线距离,
式中,r1表示从源点到地面上积分动点Q之间的直线距离,r2表示从Q点到接收点P的距离,k0为空气的波数,fg(Q)为动点Q处的衰减函数,表示r2在表面上的法向导数,在地面接收时,A=1,高空接收时A=0.5,L、D分别表示动点Q和接收点P在x轴上投影至发射点的距离,Δg是地面的归一化表面阻抗:
式中,εr为地面的相对介电常数,i为虚数单位,λ为波长,σ为大地电导率。
步骤3.4中,利用积分方法计算衰减因子时,令电导率为5S/m,介电常数为81,高程为0。
本发明的有益效果是,一种陆基长波导航/授时系统ASF预测数据库快速建库方法,采用CPU+GPU异构模式,利用积分方法的迭代过程,通过并行计算快速得到多个位置点的ASF预测结果,从而实现在较短时间内计算完成发射台站周围的ASF修正值。
附图说明
图1是本发明的流程示意图;
图2是本发明实施中区域网格化的示意图;
图3是本发明的程序流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
陆基长波导航/授时系统ASF预测数据库快速建库方法,流程图如图1所示,具体按照以下步骤实施:
步骤1、区域网格化,针对某发射台站,确定发射台站经纬度信息,并以发射台为中心,以θ角度为间隔形成辐射路径,以正北方向为起始方向,顺时针划分角度间隔,按照步长Δd等分各路径,称所有路径上的等分点为区域的网格点,获取各网格点的经纬度信息;
步骤2、获取网格点地理信息,通过长波地理信息系统(该系统集成了中国高程数据SRTM90文件、中国大地电导率电子地图、介电常数与电导率对应表,给出点的经纬度,系统可查询高程数据SRTM90文件得到点的高程,查询大地电导率电子地图得到该点的电导率,根据电导率查询介电常数与电导率的对应表得出介电常数),基于各网格点的经纬度信息,获取路径参数信息,包括高程h、大地电导率σ、相对介电常数εr,形成路径参数数据库;
步骤3、计算附加二次相位因子ASF(Additional Secondary phase Factor),根据步骤2获得的数据,通过积分方法的迭代过程,并行计算出各网格点的附加二次相位因子ASF,并将各网格点的附加二次相位因子ASF按路径存入文件中,一条路径的附加二次相位因子ASF存于一个.dat文件,文件以“发射台名+路径与正北方向的夹角θ”命名(例如:“蒲城A=100.1”),最终得到发射台附近区域的ASF预测数据库。
步骤1中θ大于0度,小于1度。
步骤1中步长Δd小于200m。
步骤3具体如下:
步骤3.1、从所述步骤2得到的路径参数数据库中读入各网格点的地理信息存入二维数组中,对高程h进行小波滤波处理,将参数数组从CPU内存复制到GPU内存,以便后续在GPU中并行计算;
步骤3.2、利用积分方法的迭代过程并行计算各点的衰减因子Wg,并存入衰减因子数组中;
步骤3.3、将在GPU计算得到的衰减因子数组复制回CPU中,地波衰减因子的相位又称为二次时延,可得到时延tw;
步骤3.4、利用积分方法计算得到一条海水路径下的衰减因子,海水路径下衰减因子的相位即为二次相位因子SF(Secondary phase Factor);
步骤3.5、根据附加二次相位因子ASF=tw-SF得到各网格点的ASF值。
步骤3.2具体如下:
其中
式中,Ez为接受点P(ρ,Φ,z)处z方向的场强,Eb为发射源和接收点P都位于理想导电平面时P点的场强,ω为地波的角频率,μ0为真空磁导率,Idl为垂直电偶极子大小,垂直电偶极子即为发射源,i为虚数单位,k0为空气的波数,kg为地面的波数,z、ρ分别为接收点P在柱坐标系下z方向、ρ方向的坐标值,h为垂直电偶极子Idl的离地高度,pr为发射源与接收点P之间的大圆距离;
从第个点开始,根据各点的地理信息以及计算得到的地面的衰减因子Wg,利用积分法公式(8)~(10)并行计算出所有路径第个点的地面的衰减因子Wg,按照此步骤,依次迭代计算出路径上剩余各点的地面的衰减因子Wg:
式中,fg是接收点处的地波衰减函数,dp为接收点与发射点的大圆距离,r0表示从发射点到接收点的直线距离,
式中,r1表示从源点到地面上积分动点Q之间的直线距离,r2表示从Q点到接收点P的距离,k0为空气的波数,fg(Q)为动点Q处的衰减函数,表示r2在表面上的法向导数,在地面接收时A=1,高空接收时A=0.5,L、D分别表示动点Q和接收点P在x轴上投影至发射点的距离,Δg是地面的归一化表面阻抗:
式中,εr为地面的相对介电常数,i为虚数单位,λ为波长,σ为大地电导率。
步骤3.4中,利用积分方法计算衰减因子时,令电导率为5S/m,介电常数为81,高程为0。
实施例
以中国科学院国家授时中心的长波发波台(蒲城)作实验,快速建立以蒲城台为中心,半径为94.607km的圆形区域的ASF预测数据库,流程图如图1所示,具体按照以下步骤实施:
步骤1、区域网格化。针对蒲城发射台,经度:109°32′35″,纬度:34°56′55″,以其为中心,以0.1°为间隔形成3600条长94.607km的辐射路径,将大圆距离映射到直角坐标系下,按照步长100m等分,每条路径922个点,如图2所示。共3600×922=3319200个网格点,通过大地主题解算获取各网格点的经纬度信息;
步骤2、获取网格点地理信息。通过长波地理信息系统,基于各等分点的经纬度信息,获取路径参数信息,包括高程h、大地电导率σ、相对介电常数εr,形成路径参数数据库;
步骤3、根据步骤2获得的数据,并行计算各点的ASF。并行计算程序流程图如图3所示。
步骤3中并行计算的流程如下:
步骤3.1、从路径参数数据库中读入各网格点的地理信息存入二维数组中,对各路径的高程h进行小波滤波处理,将参数数组从CPU拷贝到GPU;
步骤3.2、利用积分方法迭代过程并行计算各点的衰减因子Wg。
前2km可看作均匀平地面,采用平地面公式(1)~(7)可并行计算所有路径前19个点的衰减因子Wg。
式中k0、ε0分别为空气的波数和介电常数,kg、εr、σ分别为地面的波数、相对介电常数以及电导率。垂直电偶极子Idl位于z轴中离地高度为h处,观察点P的坐标为(ρ,Φ,z),pr为源点与发射点之间的大圆距离。
从第20个点开始,根据点的地理信息以及路径前19个点的衰减因子Wg,利用积分法公式(8)~(10)可并行计算出所有路径上第20个点的Wg。按照此步骤,依次迭代计算出路径剩余各点的Wg。
式中fg为接收处的地波衰减函数。dp为接收点距发射台的大圆距离。r0表示从源点到接收点的直线距离。
步骤3.3、将计算得到的衰减函数数组拷贝回CPU中,地波衰减因子的相位又称为二次时延,可得到时延tw;
步骤3.4、令电导率为5S/m,介电常数为81,高程为0,利用积分方法公式(8)~(10),可得到一条海水路径下的衰减因子,对其计算得到时延SF;
步骤3.5、根据ASF=tw-SF可得到各网格点的ASF值。
本发明耗时1714s即28:34min完成了蒲城台方圆94.607公里内3319200个网格点的ASF修正值的计算,耗时较短,实现了针对一个发射台站周围ASF修正数据的快速计算。
本发明实施例中实验环境如下表所示:
表1实施例中实验环境
本发明的实施例结果如下表所示:
表2实施例结果
由表2可以看出,在扩大计算区域与网格点数后,本发明方法依然能在较短时间内计算出发射台站周围ASF修正值。
Claims (6)
1.陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,具体按照以下步骤实施:
步骤1、区域网格化,针对某发射台站,确定发射台站经纬度信息,并以发射台为中心,以θ角度为间隔形成辐射路径,以正北方向为起始方向,顺时针划分角度间隔,按照步长Δd等分各路径,称所有路径上的等分点为区域的网格点,获取各网格点的经纬度信息;
步骤2、获取网格点地理信息,通过长波地理信息系统,基于各网格点的经纬度信息,获取路径参数信息,包括高程h、大地电导率σ、相对介电常数εr,形成路径参数数据库;
步骤3、计算附加二次相位因子ASF,根据步骤2获得的数据,通过积分方法的迭代过程,并行计算出各网格点的附加二次相位因子ASF,并将各网格点的附加二次相位因子ASF按路径存入文件中,一条路径的附加二次相位因子ASF存于一个.dat文件,文件以“发射台名+路径与正北方向的夹角θ”命名,最终得到发射台附近区域的ASF预测数据库。
2.根据权利要求1所述的陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,所述步骤1中θ大于0度,小于1度。
3.根据权利要求1所述的陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,所述步骤1中步长Δd小于200m。
4.根据权利要求1所述的陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,所述步骤3具体如下:
步骤3.1、从所述步骤2得到的路径参数数据库中读入各网格点的地理信息存入二维数组中,对高程h进行小波滤波处理,将参数数组从CPU内存复制到GPU内存,以便后续在GPU中并行计算;
步骤3.2、利用积分方法的迭代过程并行计算各点的衰减因子Wg,并存入衰减因子数组中;
步骤3.3、将在GPU计算得到的衰减因子数组复制回CPU中,地波衰减因子的相位又称为二次时延,可得到时延tw;
步骤3.4、利用积分方法计算得到一条海水路径下的衰减因子,海水路径下衰减因子的相位即为二次相位因子SF;
步骤3.5、根据附加二次相位因子ASF=tw-SF得到各网格点的ASF值。
5.根据权利要求4所述的陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,所述步骤3.2具体如下:
其中
式中,Ez为接受点P(ρ,Φ,z)处z方向的场强,Eb为发射源和接收点P都位于理想导电平面时P点的场强,ω为地波的角频率,μ0为真空磁导率,Idl为垂直电偶极子大小,垂直电偶极子即为发射源,i为虚数单位,k0为空气的波数,kg为地面的波数,z、ρ分别为接收点P在柱坐标系下z方向、ρ方向的坐标值,h为垂直电偶极子Idl的离地高度,pr为发射源与接收点P之间的大圆距离;
从第个点开始,根据各点的地理信息以及计算得到的地面的衰减因子Wg,利用积分法公式(8)~(10)并行计算出所有路径第个点的地面的衰减因子Wg,按照此步骤,依次迭代计算出路径上剩余各点的地面的衰减因子Wg:
式中,fg是接收点处的地波衰减函数,dp为接收点与发射点的大圆距离,r0表示从发射点到接收点的直线距离,
式中,r1表示从源点到地面上积分动点Q之间的直线距离,r2表示从Q点到接收点P的距离,k0为空气的波数,fg(Q)为动点Q处的衰减函数,表示r2在表面上的法向导数,在地面接收时,A=1,高空接收时A=0.5,L、D分别表示动点Q和接收点P在x轴上投影至发射点的距离,Δg是地面的归一化表面阻抗:
式中,εr为地面的相对介电常数,i为虚数单位,λ为波长,σ为大地电导率电导率。
6.根据权利要求5所述的陆基长波导航/授时系统ASF预测数据库快速建库方法,其特征在于,所述步骤3.4中,利用积分方法计算衰减因子时,令电导率为5S/m,介电常数为81,高程为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011636854.0A CN112781593B (zh) | 2020-12-31 | 2020-12-31 | 陆基长波导航/授时系统asf预测数据库快速建库方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011636854.0A CN112781593B (zh) | 2020-12-31 | 2020-12-31 | 陆基长波导航/授时系统asf预测数据库快速建库方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112781593A true CN112781593A (zh) | 2021-05-11 |
CN112781593B CN112781593B (zh) | 2023-01-10 |
Family
ID=75754961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011636854.0A Active CN112781593B (zh) | 2020-12-31 | 2020-12-31 | 陆基长波导航/授时系统asf预测数据库快速建库方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112781593B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101051872A (zh) * | 2007-05-16 | 2007-10-10 | 西安理工大学 | 一种结合asf(td)和asf(toa)测量获取asf(toa)的方法 |
CN101059543A (zh) * | 2007-05-16 | 2007-10-24 | 西安理工大学 | 一种长波传播时延修正量的测量方法 |
CN106874549A (zh) * | 2017-01-10 | 2017-06-20 | 西安理工大学 | 一种高精度预测asf的窄带离散分布抛物方程方法 |
CN107196716A (zh) * | 2017-04-21 | 2017-09-22 | 中国科学院国家授时中心 | 计算长波地波信号路径传播时延的差分方法 |
RU2645875C1 (ru) * | 2017-01-11 | 2018-02-28 | Акционерное общество научно-внедренческое предприятие "ПРОТЕК" | Способ повышения точности дифференциальной коррекции навигационных параметров в длинноволновой системе определения местоположения |
CN110990505A (zh) * | 2019-11-24 | 2020-04-10 | 中国电子科技集团公司第二十研究所 | 一种基于神经网络的Loran-C ASF修正方法 |
US20200271746A1 (en) * | 2019-02-22 | 2020-08-27 | Eagle Technology, Llc | ENHANCED LORAN (eLORAN) SYSTEM HAVING CORRECTED ADDITIONAL SECONDARY FACTOR (ASF) DATA |
-
2020
- 2020-12-31 CN CN202011636854.0A patent/CN112781593B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101051872A (zh) * | 2007-05-16 | 2007-10-10 | 西安理工大学 | 一种结合asf(td)和asf(toa)测量获取asf(toa)的方法 |
CN101059543A (zh) * | 2007-05-16 | 2007-10-24 | 西安理工大学 | 一种长波传播时延修正量的测量方法 |
CN106874549A (zh) * | 2017-01-10 | 2017-06-20 | 西安理工大学 | 一种高精度预测asf的窄带离散分布抛物方程方法 |
RU2645875C1 (ru) * | 2017-01-11 | 2018-02-28 | Акционерное общество научно-внедренческое предприятие "ПРОТЕК" | Способ повышения точности дифференциальной коррекции навигационных параметров в длинноволновой системе определения местоположения |
CN107196716A (zh) * | 2017-04-21 | 2017-09-22 | 中国科学院国家授时中心 | 计算长波地波信号路径传播时延的差分方法 |
US20200271746A1 (en) * | 2019-02-22 | 2020-08-27 | Eagle Technology, Llc | ENHANCED LORAN (eLORAN) SYSTEM HAVING CORRECTED ADDITIONAL SECONDARY FACTOR (ASF) DATA |
CN110990505A (zh) * | 2019-11-24 | 2020-04-10 | 中国电子科技集团公司第二十研究所 | 一种基于神经网络的Loran-C ASF修正方法 |
Non-Patent Citations (3)
Title |
---|
周丽丽等: "不规则地面低频地波传播预测方法的性能比较", 《电波科学学报》 * |
徐务农等: "罗兰C地波传播的ASF修正研究", 《舰船电子工程》 * |
陈秀明等: "北斗-罗兰C组合导航的ASF修正计算与研究", 《宇航计测技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112781593B (zh) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6943729B2 (en) | Method and system for time difference of arrival (TDOA) location services | |
CN103217177B (zh) | 一种电波折射修正方法、装置及系统 | |
CN108181630B (zh) | 一种北斗双天线旋转快速定向方法 | |
CN106886039B (zh) | 基于城市三维地图的地面数字电视广播与导航卫星混合粒子滤波定位方法 | |
US11906638B2 (en) | High-accuracy satellite positioning method, positioning terminal, and positioning system | |
US11204426B2 (en) | Method for determining an adaptive model of an electron density distribution | |
CN109061559B (zh) | 一种uwb基站天线相位中心偏差建模与改正的研究方法 | |
CN102427602A (zh) | 稀疏直接定位方法 | |
CN113447924B (zh) | 一种基于毫米波雷达的无人机测绘方法和系统 | |
CN107037463A (zh) | 一种地基长波授时系统与gnss系统的组合定位授时方法 | |
CN111123345A (zh) | 一种基于gnss测量的经验电离层模型数据驱动方法 | |
CN116166680B (zh) | 一种铁路北斗基准站控制网自动更新维护方法和系统 | |
CN112781593B (zh) | 陆基长波导航/授时系统asf预测数据库快速建库方法 | |
CN116755126B (zh) | 一种基于三维模型映射匹配的北斗实时精准定位方法 | |
CN111060945B (zh) | 一种gnss/5g紧组合融合定位方法与装置 | |
CN114527500B (zh) | 一种室内外一体化定位方法、设备、介质及产品 | |
CN109490828A (zh) | 基于同源基线阵列的定位方法 | |
CN212410861U (zh) | 一种基于卫星定位陆基增强的近地警告系统 | |
Yang et al. | Random optimization algorithm on GNSS monitoring stations selection for ultra-rapid orbit determination and real-time satellite clock offset estimation | |
CN114594435A (zh) | 一种国产民用sar卫星几何检校与定位精度改进方法 | |
CN112099058A (zh) | 定位方法、装置、系统、计算机设备和存储介质 | |
CN114280655A (zh) | 一种基于全球导航卫星系统的测姿方法和系统 | |
Hayward et al. | Real time calibration of antenna phase errors for ultra short baseline attitude systems | |
CN112731497A (zh) | 一种利用现场视域分析提高卫星定位精度的方法 | |
Ashkenazi | Models for controlling national and continental networks |
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 |