CN112711050A - 基于椭球变换的星载gnss-r镜面反射点计算方法 - Google Patents
基于椭球变换的星载gnss-r镜面反射点计算方法 Download PDFInfo
- Publication number
- CN112711050A CN112711050A CN202011501293.3A CN202011501293A CN112711050A CN 112711050 A CN112711050 A CN 112711050A CN 202011501293 A CN202011501293 A CN 202011501293A CN 112711050 A CN112711050 A CN 112711050A
- Authority
- CN
- China
- Prior art keywords
- ellipsoid
- reflection point
- satellite
- gnss
- transformation
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/42—Determining position
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/68—Analysis of geometric attributes of symmetry
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Aerials With Secondary Devices (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明介绍了一种基于椭球变换的星载GNSS‑R镜面反射点计算方法,通过数字高程模型DEM或海面地形模型,确定概略镜面反射点处大地高,并以此大地高改化椭球长半轴值,构建新的参考椭球面,并以新参考椭球面作为参考,依据导航卫星位置、GNSS‑R卫星位置和菲涅耳反射定律确定镜面反射点,提高镜面反射点定位精度。
Description
技术领域
本发明属于水下潜器的定位方法技术领域,具体涉及一种基于椭球变换的星载GNSS-R镜面反射点计算方法。
背景技术
利用地球表面反射的卫星导航信号进行地球物理参数的探测是全球导航卫星系统(GNSS)一个创新性的研究热点,该项技术通常被称为GNSS-R技术,它具有全天候、全天时、多信号源、高时空分辨率等优势,主要为气象、海洋、农业、国土资源、灾害应急等部门提供业务化支撑,包括海面风场、海面高度、土壤湿度、内陆水资源和湿地分布、洪水监测等。
利用GNSS反射信号进行测量,需要知道导航卫星轨道和GNSS-R接收平台的位置或轨迹,同时,GNSS反射信号的接收处理,需要估计地表反射区参考点(镜面反射点)的位置。
在轨数据处理过程中,通常采用WGS84或GRS80参考椭球面作为参考,以此计算镜面反射点位置;但是,在实际情况中,反射点位于实际陆地地形表面或海面地形表面。当参考椭球面与实际地形表面有较大差异时,将会降低GNSS反射信号观测量的精度,甚至导致对GNSS反射信号的捕获。因此,有必要研究更高精度的星载GNSS-R镜面反射点计算方法。
发明内容
本发明的目的在于提供一种基于椭球变换的星载GNSS-R镜面反射点计算方法,以解决背景技术中存在的问题。
为实现上述目的,本发明提供如下技术方案:一种基于椭球变换的星载GNSS-R镜面反射点计算方法,其特征在于,包括以下步骤:
步骤一,镜面反射点概略位置计算,以WGS84或GRS80参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置和菲涅耳反射定律,计算出镜面反射点概略位置;
步骤二,椭球变换,结合数字高程模型DEM或海面地形模型,获取该概略镜面反射点大地高,继而利用此大地高改化椭球长半轴,但保持椭球扁率不变,从而形成新的参考椭球面;
步骤三,镜面反射点准确位置计算,以新参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置和菲涅耳反射定律,计算出镜面反射点准确位置,提高定位精度。
优选地,步骤一中建立镜面反射点计算约束条件:镜面反射点位于参考椭球面上;直射信号向量、反射信号向量与过镜面反射点处法线向量等三者共面;遵循菲涅耳反射定律,入射角等于反射角。
优选地,步骤二中进行椭球变换,使得镜面反射点位于新参考椭球面上,其关键在于通过概略镜面反射点位置,结合数字高程模型DEM或海面地形模型,获取该点大地高,继而改化椭球长半轴值,但保持椭球扁率不变,从而形成新的参考椭球面。
优选地,步骤三中以新参考椭球面为参考,并按照步骤一中的三项约束条件求解镜面反射点位置,提高定位精度。
本发明与现有技术相比,其有益效果是:
本发明提出的一种基于椭球变换的星载GNSS-R镜面反射点计算方法,通过数字高程模型DEM或海面地形模型,确定概略镜面反射点处大地高,并以此大地高改化椭球长半轴值,构建新的参考椭球面,并以新参考椭球面作为参考确定镜面反射点,提高镜面反射点定位精度,与实际情况更加符合。同时该方法计算简单,适合于卫星在轨数据处理。
附图说明
图1是本发明计算镜面反射点示意图;
图2是本发明变换参考椭球示意图。
具体实施方式
下面结合附图及较佳实施例详细说明本发明的具体实施方式。
本发明介绍了一种基于椭球变换的星载GNSS-R镜面反射点计算方法,具体是指通过数字高程模型DEM或海面地形模型,确定概略镜面反射点处大地高,并以此大地高改化椭球长半轴值,构建新的参考椭球面,并以新参考椭球面作为参考确定镜面反射点,提高镜面反射点定位精度。
本发明基于椭球变换的星载GNSS-R镜面反射点计算方法,由镜面反射点概略位置计算、椭球变换、镜面反射点准确位置计算三个主要部分组成,具体过程如下:
步骤一,镜面反射点概略位置计算。以WGS84或GRS80参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置,以及三项约束条件,计算出镜面反射点概略位置。其中,约束条件一镜面反射点位于参考椭球面上如公式1所示,约束条件二直射信号向量、反射信号向量与过镜面反射点处法线向量等三者共面如公式2所示,约束条件三入射角等于反射角如公式3所示。
约束条件二可表示为:
约束条件三可表示为:
步骤二,椭球变换。结合数字高程模型DEM或海面地形模型,获取该概略镜面反射点大地高,镜面反射点大地坐标可表示为(B,L,H),继而根据公式4,利用此大地高改化椭球长半轴,但保持椭球扁率不变,从而形成新的参考椭球面,新参考椭球面长短半轴分别为anew,bnew。
步骤三,镜面反射点准确位置计算。以新参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置和菲涅耳反射定律,按照步骤一所列约束条件,计算出镜面反射点准确位置,提高定位精度。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.一种基于椭球变换的星载GNSS-R镜面反射点计算方法,其特征在于,包括以下步骤:
步骤一,镜面反射点概略位置计算,以WGS84或GRS80参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置和菲涅耳反射定律,计算出镜面反射点概略位置;
步骤二,椭球变换,结合数字高程模型DEM或海面地形模型,获取该概略镜面反射点大地高,继而利用此大地高改化椭球长半轴,但保持椭球扁率不变,从而形成新的参考椭球面;
步骤三,镜面反射点准确位置计算,以新参考椭球面为参考,利用导航卫星位置、GNSS-R卫星位置和菲涅耳反射定律,计算出镜面反射点准确位置,提高定位精度。
2.根据权利要求1所述的基于椭球变换的星载GNSS-R镜面反射点计算方法,其特征在于,步骤一中建立镜面反射点计算约束条件:镜面反射点位于参考椭球面上;直射信号向量、反射信号向量与过镜面反射点处法线向量等三者共面;遵循菲涅耳反射定律,入射角等于反射角。
3.根据权利要求1所述的基于椭球变换的星载GNSS-R镜面反射点计算方法,其特征在于,步骤二中进行椭球变换,使得镜面反射点位于新参考椭球面上,其关键在于通过概略镜面反射点位置,结合数字高程模型DEM或海面地形模型,获取该点大地高,继而改化椭球长半轴值,但保持椭球扁率不变,从而形成新的参考椭球面。
4.根据权利要求1所述的基于椭球变换的星载GNSS-R镜面反射点计算方法,其特征在于,步骤三中以新参考椭球面为参考,并按照步骤一中的三项约束条件求解镜面反射点位置,提高定位精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011501293.3A CN112711050B (zh) | 2020-12-17 | 2020-12-17 | 基于椭球变换的星载gnss-r镜面反射点计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011501293.3A CN112711050B (zh) | 2020-12-17 | 2020-12-17 | 基于椭球变换的星载gnss-r镜面反射点计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112711050A true CN112711050A (zh) | 2021-04-27 |
CN112711050B CN112711050B (zh) | 2022-08-23 |
Family
ID=75544460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011501293.3A Active CN112711050B (zh) | 2020-12-17 | 2020-12-17 | 基于椭球变换的星载gnss-r镜面反射点计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112711050B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101532834A (zh) * | 2008-12-12 | 2009-09-16 | 北京理工大学 | WGS84椭球与Clarke80椭球的坐标转换方法 |
US20100007547A1 (en) * | 2008-07-11 | 2010-01-14 | Agence Spatiale Europeenne | Altimetry method and system |
CN107037460A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种用于海面测高的北斗meo卫星机载反射信号建模方法 |
CN107561562A (zh) * | 2017-08-17 | 2018-01-09 | 上海海洋大学 | 一种gnss‑r遥感中镜面反射点快速确定方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换系统 |
CN109696697A (zh) * | 2019-02-01 | 2019-04-30 | 中国空间技术研究院 | Gnss-r镜面反射点的大地水准面-垂线偏差修正方法和系统 |
CN110988942A (zh) * | 2019-11-29 | 2020-04-10 | 航天东方红卫星有限公司 | 一种星载gnss-r镜面反射点位置精确计算方法 |
-
2020
- 2020-12-17 CN CN202011501293.3A patent/CN112711050B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100007547A1 (en) * | 2008-07-11 | 2010-01-14 | Agence Spatiale Europeenne | Altimetry method and system |
CN101532834A (zh) * | 2008-12-12 | 2009-09-16 | 北京理工大学 | WGS84椭球与Clarke80椭球的坐标转换方法 |
CN107037460A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种用于海面测高的北斗meo卫星机载反射信号建模方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换系统 |
CN107561562A (zh) * | 2017-08-17 | 2018-01-09 | 上海海洋大学 | 一种gnss‑r遥感中镜面反射点快速确定方法 |
CN109696697A (zh) * | 2019-02-01 | 2019-04-30 | 中国空间技术研究院 | Gnss-r镜面反射点的大地水准面-垂线偏差修正方法和系统 |
CN110988942A (zh) * | 2019-11-29 | 2020-04-10 | 航天东方红卫星有限公司 | 一种星载gnss-r镜面反射点位置精确计算方法 |
Non-Patent Citations (2)
Title |
---|
李晓华 等: "基于DEM的CGCS2000城市平面坐标系建立关键技术研究", 《测绘地理信息》 * |
陈新焕: "GPS施工控制网的坐标转换", 《铁路航测》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112711050B (zh) | 2022-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Specht et al. | Assessment of the positioning accuracy of DGPS and EGNOS systems in the Bay of Gdansk using maritime dynamic measurements | |
CN105842724B (zh) | 一种船舶辅助泊岸方法和系统 | |
CN102636149B (zh) | 挠性体动态变形的组合测量装置及测量方法 | |
WO2024007365A1 (zh) | 一种基于北斗/gnss的实时高精度海表测量方法及浮标 | |
CN106842268A (zh) | 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法 | |
CN109613582B (zh) | 一种车载实时单频米级伪距定位方法 | |
CN109059916A (zh) | 一种基于惯导的浮空器掩星预报方法 | |
CN108489497B (zh) | 一种利用地图防触礁的安全助航方法 | |
CN113686299A (zh) | 一种海上动态目标定位与移速预测方法 | |
Li et al. | A precise underwater positioning method by considering the location difference of transmitting and receiving sound waves | |
Zhu et al. | Monitoring of Yangtze river discharge at datong hydrometric station using acoustic tomography technology | |
Roggenbuck et al. | Oceanwide precise determination of sea surface height from in-situ measurements on cargo ships | |
Lu et al. | Motion pose estimation of inshore ships based on point cloud | |
CN115808686A (zh) | 基于波形重跟踪质量加权的垂线偏差格网化方法及其系统 | |
CN112711050B (zh) | 基于椭球变换的星载gnss-r镜面反射点计算方法 | |
CN110046430B (zh) | 基于椭球分层大气折射的光学卫星影像精确对地定位方法 | |
Li et al. | Adaptively robust filtering algorithm for maritime celestial navigation | |
Ng et al. | Range-based 3D mapping aided GNSS with NLOS correction based on skyplot with building boundaries | |
CN116105685A (zh) | 基于声光遥感和滚轮的潮间带地形无缝一体化测量方法 | |
CN110136487A (zh) | 地图导航防触礁安全助航技术应用模型 | |
CN116106925A (zh) | 严密的光子计数机制激光雷达水下测深点坐标计算方法 | |
Szulwic et al. | Mobile laser scanning calibration on a marine platform | |
De Girolamo et al. | Wave characteristics estimation by GPS receivers installed on a sailboat travelling off-shore | |
Wells et al. | GPS: A multipurpose system." | |
CN113625312B (zh) | 一种gps-r/bds-r反射延迟海况偏差的量化与预测方法及系统 |
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 |