CN111275757A - 一种基于dem数据处理的伪卫星场地仿真布设方法 - Google Patents

一种基于dem数据处理的伪卫星场地仿真布设方法 Download PDF

Info

Publication number
CN111275757A
CN111275757A CN202010018185.4A CN202010018185A CN111275757A CN 111275757 A CN111275757 A CN 111275757A CN 202010018185 A CN202010018185 A CN 202010018185A CN 111275757 A CN111275757 A CN 111275757A
Authority
CN
China
Prior art keywords
information
field
dem data
data
pseudolite
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
CN202010018185.4A
Other languages
English (en)
Other versions
CN111275757B (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.)
CETC 54 Research Institute
Original Assignee
CETC 54 Research Institute
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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202010018185.4A priority Critical patent/CN111275757B/zh
Publication of CN111275757A publication Critical patent/CN111275757A/zh
Application granted granted Critical
Publication of CN111275757B publication Critical patent/CN111275757B/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
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • 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
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/001Texturing; Colouring; Generation of texture or colour
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种基于DEM数据处理的伪卫星场地仿真布设方法,属于伪卫星导航定位技术领域。其包括以下步骤:将原始DEM数据转换成平面坐标数据;生成具有高程颜色的图像数据;等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,并计算出各伪卫星在场地中的实际位置坐标;通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。本发明对DEM数据进行计算,不必盲目搜索计算区域内所有位置,只需计算DEM点位的仿真信息即可,具有计算量小、精度高、效果可靠的优点,此外,通过输入各导航系统的卫星历书信息,还可实现导航系统和伪卫星的混合仿真。

Description

一种基于DEM数据处理的伪卫星场地仿真布设方法
技术领域
本发明涉及伪卫星导航定位技术领域,特别是指一种基于DEM数据处理的伪卫星场地仿真布设方法。
背景技术
目前,以北斗和GPS为代表的卫星导航系统应用越来越广泛。然而,在某些情况下,如在城市“峡谷”或者山谷中,由于可视卫星的数目太少或几何图形很差,从而不能可靠地确定定位。在室内,由于卫星导航信号相当微弱,或者在某些地方根本接收不到卫星导航信号,利用卫星进行定位几乎是不可能的。为解决这些问题,在适当位置上布置一定数目的伪卫星,有利于构成具有足够强度的几何图形,从而在很大程度上提高导航定位解的精度。
现有技术中,伪卫星布设一般都是基于位置搜索方法。当伪卫星和用户处于同一高度时,在用户附近用方位扫描;当伪卫星和用户在不同高度时,用直角网格扫描的方法来进行系统搜索。这种伪卫星位置的选择方式要经历一个冗长的搜索过程,若采用立体直角格网扫描,要进行上亿次的DOP计算。可见,现有技术中在进行伪卫星布设时需要进行大量的计算,布设工作不易实施。
发明内容
有鉴于此,本发明提出一种基于DEM数据处理的伪卫星场地仿真布设方法,该方法在获取DEM数据的基础上对DEM数据进行坐标处理和JPG处理,进而对伪卫星布设场地设计进行仿真分析,可以为用户提供各种布设场景下的伪卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布等信息,并且可通过输入各导航系统的卫星历书信息进行导航系统和伪卫星的混合仿真。
为了实现上述目的,本发明所采取的技术方案为:
一种基于DEM数据处理的伪卫星场地仿真布设方法,其包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;
(4)根据各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
进一步的,所述步骤(2)的具体方式为:
(201)根据DEM数据中的点位的纬度和经度信息确定点位的位置;
(202)利用DEM数据中的点位的高程信息确定点位的颜色代码;
(203)利用步骤(201)和(202)所得到的DEM数据中的全部点位的位置和颜色代码,通过图像转换方法生成一幅以颜色标示各点高程的图像数据。
此外,本发明还提供一种基于DEM数据处理的伪卫星场地仿真布设方法,其包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据文件转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;
(4)输入真实卫星的导航历书,计算指定时间时真实卫星的位置坐标;根据各真实卫星的位置坐标以及各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
本发明采用上述技术方案所取得的有益效果在于:
1、本发明利用DEM数据提供的点位的精确坐标信息生成场地图像,由此得到的仿真结果精度高、效果可靠。
2、本发明对DEM数据进行计算,不必盲目搜索计算区域内所有位置,只需计算DEM点位的仿真信息即可,相对于原有的位置搜索方法计算量有大幅减少。
3、本发明还可输入各导航系统的卫星历书信息,从而实现导航系统和伪卫星的混合仿真。
附图说明
图1是本发明实施例的实现原理示意图。
图2是本发明实施例的方法流程图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步的描述。
一种基于DEM数据处理的伪卫星场地仿真布设方法,其包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;布设伪卫星时可以采用鼠标点选的方式直接在场地图像上进行操作,也可通过文本框直接输入位置坐标;
(4)根据各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
进一步的,所述步骤(2)的具体方式为:
(201)根据DEM数据中的点位的纬度和经度信息确定点位的位置;
(202)利用DEM数据中的点位的高程信息确定点位的颜色代码;颜色代码可采用RGB形式;
(203)利用步骤(201)和(202)所得到的DEM数据中的全部点位的位置和颜色代码,通过图像转换方法生成一幅以颜色标示各点高程的图像数据,生成的图像可以为jpg格式。
此外,对上述方法稍加修改还可实现导航系统和伪卫星的混合仿真,具体包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据文件转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;
(4)输入真实卫星的导航历书,计算指定时间时真实卫星的位置坐标;根据各真实卫星的位置坐标以及各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
图1为上述方法的一种具体实现方式,具体来说,它实现了DEM坐标转换模块、JPG处理模块、伪卫星场地仿真模块、导航系统/伪卫星场地混合仿真分析模块。其中:
DEM坐标转换模块的功能是将DEM数据投影到平面坐标数据文件。数字地面模型DTM是利用一个任意坐标系中的已知纬度、经度、高程的坐标点对连续地面的一种模拟标示。DEM数据中的纬度、经度、高程的坐标系可能为GPS卫星导航系统所用坐标系WGS84、BDS卫星导航系统所用坐标系统CGCS2000和GLONASS卫星导航系统所用坐标系PZ90等系统坐标系统,这些坐标系统全部为大地坐标系统,需要将大地坐标数据经过高斯投影换到平面坐标数据。
JPG处理模块的功能是生成JPG数据文件以供伪卫星场地布设提供人机交互显示。用户在布设伪卫星时需要知道布设场地的详细地理环境,JPG数据文件通过颜色的深浅来区分仿真场地的地面信息,通过横坐标和纵坐标显示位置信息。JPG数据文件中地图信息是和平面坐标数据文件中的坐标相匹配。
伪卫星场地仿真模块的功能是利用平面坐标数据和JGP数据文件来仿真用户确定的场地。用户可以通过在所选的场地内,通过不同方案布设伪卫星信息,如伪卫星的个数、伪卫星的位置、发射天线参数、发射功率等参数实现不同的布设场景。布设后通过仿真分析出不同场景的卫星可视性分布、HDOP分布、VDOP分布、平均仰角分布、接收功率等布设信息。用户通过分析不同场景的仿真结果来决定伪卫星的布设方案。
导航系统/伪卫星场地混合仿真分析模块通过输入卫星导航系统的卫星历书,预报场地内导航卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布、接收功率等信息,并结合场地内布设的伪卫星仿真结果综合出导航系统/伪卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布信息。
图2为上述方法的一个更加具体的流程,其包括如下步骤:
步骤1:DEM坐标转换模块输入原始DEM数据文件和DEM数据的坐标系统。进行高斯投影转换成平面坐标数据文件。
步骤2:JPG处理模块输入原始DEM数据文件和DEM数据的坐标系统。进行JPG数据文件生成。
步骤3:伪卫星场地仿真模块功能通过利用平面坐标数据文件和JGP数据文件来仿真用户确定的场地。用户可以通过布设多个伪卫星站并配置各个伪卫星的位置、天线参数、发射功率等参数,进而仿真计算出卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布等信息;
步骤4:导航系统/伪卫星场地混合仿真分析模块输入卫星导航历书进行预测导航系统卫星在仿真场地内的卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布。并叠加场地内布设的伪卫星综合出导航系统/伪卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布等场地布设信息。
总之,本发明方法在获取DEM数据的基础上对DEM数据进行坐标处理和JPG处理,进而对伪卫星布设场地设计进行仿真分析,能够为用户提供各种布设场景下的伪卫星的可视性分布、HDOP分布、VDOP分布、平均仰角分布等信息,并且可通过输入各导航系统的卫星历书信息进行导航系统和伪卫星的混合仿真。
需要注意的是,以上所述仅为本发明的若干具体实施例,但本发明的保护范围并不局限于此,在本发明揭露的技术范围内,可以想到的变换,都应涵盖在本发明的保护范围内。

Claims (3)

1.一种基于DEM数据处理的伪卫星场地仿真布设方法,其特征在于,包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;
(4)根据各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
2.根据权利要求1所述的一种基于DEM数据处理的伪卫星场地仿真布设方法,其特征在于,所述步骤(2)的具体方式为:
(201)根据DEM数据中的点位的纬度和经度信息确定点位的位置;
(202)利用DEM数据中的点位的高程信息确定点位的颜色代码;
(203)利用步骤(201)和(202)所得到的DEM数据中的全部点位的位置和颜色代码,通过图像转换方法生成一幅以颜色标示各点高程的图像数据。
3.一种基于DEM数据处理的伪卫星场地仿真布设方法,其特征在于,包括以下步骤:
(1)输入场地的原始DEM数据以及DEM数据的坐标信息,利用DEM数据的坐标信息,通过高斯投影将原始DEM数据文件转换成平面坐标数据;
(2)根据原始DEM数据和DEM数据的坐标信息,生成具有高程颜色的图像数据;
(3)根据图像数据显示场地图像,等待用户在场地图像中布设伪卫星并配置各个伪卫星的天线参数和发射功率,根据伪卫星在场地图像中的位置,结合平面坐标数据,计算出各伪卫星在场地中的实际位置坐标;
(4)输入真实卫星的导航历书,计算指定时间时真实卫星的位置坐标;根据各真实卫星的位置坐标以及各伪卫星的位置坐标、天线参数和发射功率,通过仿真计算出场地中各点的可视性信息、HDOP信息、VDOP信息以及平均仰角信息,并以可视化方式在场地图像中显示出来。
CN202010018185.4A 2020-01-08 2020-01-08 一种基于dem数据处理的伪卫星场地仿真布设方法 Active CN111275757B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010018185.4A CN111275757B (zh) 2020-01-08 2020-01-08 一种基于dem数据处理的伪卫星场地仿真布设方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010018185.4A CN111275757B (zh) 2020-01-08 2020-01-08 一种基于dem数据处理的伪卫星场地仿真布设方法

Publications (2)

Publication Number Publication Date
CN111275757A true CN111275757A (zh) 2020-06-12
CN111275757B CN111275757B (zh) 2023-04-07

Family

ID=71000161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010018185.4A Active CN111275757B (zh) 2020-01-08 2020-01-08 一种基于dem数据处理的伪卫星场地仿真布设方法

Country Status (1)

Country Link
CN (1) CN111275757B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948675A (zh) * 2020-07-10 2020-11-17 中国人民解放军61081部队 一种地基增强卫星导航系统的伪卫星布设方法
CN113489509A (zh) * 2021-07-06 2021-10-08 阎镜予 大规模gnss伪卫星间时频同步方法、装置
CN114095072A (zh) * 2021-11-15 2022-02-25 东南大学 伪卫星定位信号处理分析仿真平台及仿真方法
CN115015969A (zh) * 2022-06-02 2022-09-06 长安大学 一种山区遮挡环境下的gnss卫星可视性预报方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030142523A1 (en) * 2002-01-25 2003-07-31 Zoltan Biacs Method and system for storage and fast retrieval of digital terrain model elevations for use in positioning systems
CN101833104A (zh) * 2010-04-27 2010-09-15 北京航空航天大学 一种基于多传感器信息融合的三维可视化导航方法
CN105387846A (zh) * 2015-10-26 2016-03-09 中国农业大学 卫星影像的正射校正方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030142523A1 (en) * 2002-01-25 2003-07-31 Zoltan Biacs Method and system for storage and fast retrieval of digital terrain model elevations for use in positioning systems
CN101833104A (zh) * 2010-04-27 2010-09-15 北京航空航天大学 一种基于多传感器信息融合的三维可视化导航方法
CN105387846A (zh) * 2015-10-26 2016-03-09 中国农业大学 卫星影像的正射校正方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JIANCAI SONG ET AL.: "《Study of Constellation Design of Pseudolites Based on Improved Adaptive Genetic Algorithm》" *
刘超等: "《GPS/伪卫星组合定位中伪卫星改进选址模型》" *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948675A (zh) * 2020-07-10 2020-11-17 中国人民解放军61081部队 一种地基增强卫星导航系统的伪卫星布设方法
CN111948675B (zh) * 2020-07-10 2021-10-29 中国人民解放军61081部队 一种地基增强卫星导航系统的伪卫星布设方法
CN113489509A (zh) * 2021-07-06 2021-10-08 阎镜予 大规模gnss伪卫星间时频同步方法、装置
CN113489509B (zh) * 2021-07-06 2023-05-26 上海方位角数据科技有限公司 大规模gnss伪卫星间时频同步方法、装置
CN114095072A (zh) * 2021-11-15 2022-02-25 东南大学 伪卫星定位信号处理分析仿真平台及仿真方法
CN114095072B (zh) * 2021-11-15 2024-02-09 东南大学 伪卫星定位信号处理分析仿真平台及仿真方法
CN115015969A (zh) * 2022-06-02 2022-09-06 长安大学 一种山区遮挡环境下的gnss卫星可视性预报方法

Also Published As

Publication number Publication date
CN111275757B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN111275757B (zh) 一种基于dem数据处理的伪卫星场地仿真布设方法
KR100721764B1 (ko) Gps 및 토탈스테이션이 결합된 이동식 토탈측량시스템및 이를 이용한 측량법
Hill et al. A new era in spatial data recording: Low-cost GNSS
CN103885069B (zh) 基于cors的长输管道gps控制测量方法
CN104101888A (zh) 一种适合于中国区域的gnss电离层延迟精确建模方法
CN104123695A (zh) 一种实现坐标转换方法
KR100496814B1 (ko) Gps 측량을 이용한 도로좌표정보 취득 및 수치지도 제작방법
CN110286395A (zh) 一种北斗系统定位精度计算方法
EP3736610B1 (en) Augmented reality system for electromagnetic buried asset location
CN110196441B (zh) 终端的定位方法、装置、存储介质与设备
Kumar et al. The global positioning system: Popular accuracy measures
CN111060941A (zh) 一种在遮挡环境下的高精度定位方法及其装置
US10489985B1 (en) Augmented reality system for electromagnetic buried asset location and identification
CN110795518B (zh) 一种基于vc++的gnss定位导航软件系统
CN111521150A (zh) 一种基于gnss卫星定位的平地数据处理方法及装置
CN114511762A (zh) 一种陆上水下地形一体化测图方法和系统
CN111988842A (zh) 一种基于云端和边缘计算的手机定位方法与系统
EP3255465B1 (en) Buried asset locate device motion sensing for quality control
Naharudin et al. GIS Data Collection for Pedestrian Facilities and Furniture Using MAPinr for Android
Tedrumpun et al. Site measurement with combination of pocket PC and global positioning system for preliminary architectural design process
Mika Examination of the Capabilities of Some Models of Hand GPS Receivers for Purposes of Creation of Thematic Maps
Amran et al. Horizontal coordinate accuracy of google earth on the coverage of small islands of Makassar City, Indonesia
Șmuleac et al. 3D land modelling using gps technology in Bencecu de Sus, Timis County, Romania.
Zhang et al. Evaluation of Positioning Performance using the GPS+ QZSS in Shanghai
Wu et al. A study of GPS/Galileo performance in urban environment using a simulation tool

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