CN110456388B - 一种星载gnss-r海面高度要素定标装置及方法 - Google Patents

一种星载gnss-r海面高度要素定标装置及方法 Download PDF

Info

Publication number
CN110456388B
CN110456388B CN201910869845.7A CN201910869845A CN110456388B CN 110456388 B CN110456388 B CN 110456388B CN 201910869845 A CN201910869845 A CN 201910869845A CN 110456388 B CN110456388 B CN 110456388B
Authority
CN
China
Prior art keywords
gnss
antenna
sea surface
height
direct
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
CN201910869845.7A
Other languages
English (en)
Other versions
CN110456388A (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.)
Shanghai Spaceflight Institute of TT&C and Telecommunication
Original Assignee
Shanghai Spaceflight Institute of TT&C and Telecommunication
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 Shanghai Spaceflight Institute of TT&C and Telecommunication filed Critical Shanghai Spaceflight Institute of TT&C and Telecommunication
Priority to CN201910869845.7A priority Critical patent/CN110456388B/zh
Publication of CN110456388A publication Critical patent/CN110456388A/zh
Application granted granted Critical
Publication of CN110456388B publication Critical patent/CN110456388B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/14Receivers specially adapted for specific applications
    • 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/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/27Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
    • 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/40Correcting position, velocity or attitude
    • 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
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

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

Abstract

本发明提供了一种星载GNSS‑R海面高度要素定标装置及方法,所述方法包括:S1:计算GNSS直射信号在海面进行镜面反射的反射点位置;S2:计算GNSS反射信号与直射信号的延迟差ρsignal,以此作为海面高度的初步计算结果;S3:计算GNSS信号在传播中的电离层延迟误差ρion;S4:计算GNSS信号在传播中的对流层延迟误差ρtro;S5:计算天线位置误差ρant和天线高度差Hant;S6:根据步骤S1‑S5计算出的数据,来计算镜面反射点所在位置的海面的高度H,所述海面高度为海面到地球椭球面的绝对高度。本发明可以更加准确地计算星载平台下GNSS‑R海面高度要素,提高其测高精度。

Description

一种星载GNSS-R海面高度要素定标装置及方法
技术领域
本发明涉及卫星导航技术,具体涉及一种星载GNSS-R海面高度要素定标装置及方法。
背景技术
卫星导航系统反射信号技术(Global Navigation Satellite System-Reflectometry:GNSS-R)能够同时接收多颗GNSS卫星的信号,且其观测范围广,不受天气(云、雾)等影响,能够对海面高度测量测量实现低成本、高时空分辨率全天候的实时观测。
传统的GNSS-R海面(或湖面、以及海冰面)测高,大多是将接收机以及接收天线等测量装置搭载在地基或者是空基平台上,以此开展海面高度要素测量时,由于接收平台较矮,很少考虑电离层、对流层、直射天线和反射天线之间的位置偏差等因素带来的测高误差。
随着GNSS-R技术的不断发展,美国的NASA和欧洲的ESA先后开展了星载GNSS-R的测高实验,由于星载平台在电离层和对流层之上,因此必须考虑由于电离层和对流层带来的测量误差。并且由于卫星的复杂设计,直射天线与反射天线之间有较大的位置偏差,因此也需考虑直射天线和反射天线之间的位置偏差所导致的误差。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种星载GNSS-R海面高度要素定标装置及方法。技术方案如下:
一种星载GNSS-R海面高度要素定标装置,包括:镜面反射点计算模块、海面高度初步计算模块、电离层误差计算模块、对流层误差计算模块、天线位置误差计算模块、海面高度要素计算模块,其中:
所述镜面反射点计算模块,根据发射机位置、接收机位置及其相对应的几何关系,采用线段二分或是角度二分的算法模型,来计算GNSS直射信号在海面进行镜面反射的反射点位置,并将镜面反射点位置传给所述海面高度要素计算模块;其中,所述发射机为GNSS卫星;所述接收机为GNSS-R卫星载荷;
所述海面高度初步计算模块,根据接收机接收到的GNSS直射信号以及从海面反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal;以此作为海面高度的初步计算结果,并传给所述海面高度要素计算模块;
所述电离层误差计算模块,根据接收机的高度,采用电离层经验模型或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差ρion,并将电离层延迟误差的计算结果传给所述海面高度要素计算模块;所述GNSS信号包括GNSS直射信号和反射信号;
所述对流层误差计算模块,根据接收机的高度,采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差ρtro,并将对流层延迟误差的计算结果传给所述海面高度要素计算模块;
所述天线位置误差计算模块,根据直射天线和反射天线的搭载位置,计算天线位置误差ρant和天线高度差Hant,并将二者的计算结果传给海面高度要素计算模块;其中:
根据GNSS直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,也称为天线位置误差ρant,并计算直射天线所在平面和反射天线所在平面的高度差,也称为天线高度差Hant,并将天线位置误差和天线高度差的计算结果传给海面高度要素计算模块;
所述海面高度要素计算模块,用于根据之前计算的镜面反射点位置、海面高度的初步计算结果、电离层延迟误差、对流层延迟误差、天线位置误差和天线高度差,来计算镜面反射点所在位置的海面高度H,所述海面高度为海面到地球椭球面的绝对高度。
可选地,所述发射机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
可选地,所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
可选地,所述海面高度要素计算模块,进一步包括:
海面相对高度差计算子模块,用于计算GNSS卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure BDA0002202535880000021
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
海面高度计算子模块,根据GNSS卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,用于计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。
可选地,所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面或湖面或海冰面反射回来的GNSS卫星的信号。
一种星载GNSS-R海面高度要素定标方法,包括如下步骤:
S1:根据发射机位置、接收机位置及其相对应的几何关系,采用线段二分或是角度二分的算法模型,计算GNSS直射信号在海面进行镜面反射的反射点位置;其中,所述发射机为GNSS卫星;所述接收机为GNSS-R卫星载荷;
S2:根据接收机接收到的GNSS直射信号以及从海面反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal
S3:根据接收机的高度,采用电离层经验模型或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差ρion;所述GNSS信号包括GNSS直射信号和反射信号;
S4:根据接收机的高度,采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差ρtro
S5:根据GNSS直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,即天线位置误差ρant;并计算直射天线所在平面和反射天线所在平面的高度差,即天线高度差Hant
S6:根据步骤S1-S5计算出的镜面反射点位置、海面高度的初步计算结果、电离层延迟误差、对流层延迟误差、天线位置误差和天线高度差,来计算镜面反射点所在位置的海面高度H,所述海面高度为海面到地球椭球面的绝对高度。
可选地,所述发射机的位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
可选地,所述步骤S6进一步包括:
S61:计算GNSS卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure BDA0002202535880000041
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
S62:再根据GNSS卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。
可选地,所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
可选地,所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面或湖面或海冰面反射回来的GNSS卫星的信号。
与现有技术相比,本发明具有如下的有益效果:
本发明可以更加准确地计算星载平台下GNSS-R海面高度要素,提高海面高度要素的测高精度。
本发明的星载GNSS-R海面高度要素定标装置和方法科学合理、易于实现,通过计算在星载平台下由于电离层、对流层、直射天线和反射天线之间的位置偏差所引起的误差,从而进行有效的定标与校正,改善海面高度要素的反演精度,有利于评估星载GNSS-R载荷实际在轨的定标效果,有针对性的改善海面高度要素反演方案,进而提高星载GNSS-R载荷的海面高度要素反演结果。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1是本发明具体实施例一种星载GNSS-R海面高度要素定标装置的结构示意图;
图2是本发明具体实施例一种星载GNSS-R海面高度要素定标方法的流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1,一种星载GNSS-R海面高度要素定标装置,包括:镜面反射点计算模块、海面高度初步计算模块、电离层误差计算模块、对流层误差计算模块、天线位置误差计算模块、海面高度要素计算模块,其中:
所述镜面反射点计算模块是根据发射机(GNSS卫星)位置、接收机(GNSS-R卫星载荷)位置、及其相对应的几何关系,采用线段二分或是角度二分的算法模型,来计算GNSS卫星直射信号在海面(或湖面、以及海冰面)进行镜面反射的反射点位置,并将镜面反射点位置传给海面高度要素计算模块。
所述海面高度初步计算模块是根据接收机接收到的GNSS直射信号以及从海面(或湖面、以及海冰面)反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal,以此作为海面高度的初步计算结果,并传给海面高度要素计算模块。
所述电离层误差计算模块是根据GNSS-R接收机的高度(星载情况下,其高度通常大于几百米,受电离层影响),采用电离层经验模型(半经验模型)或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差(ρion),并将电离层误差的计算结果传给海面高度要素计算模块。
所述对流层误差计算模块是根据GNSS-R接收机的高度(星载情况下,受对流层影响),采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差(ρtro),并将对流层误差的计算结果传给海面高度要素计算模块。所述GNSS信号包括GNSS直射信号和反射信号。
所述天线位置误差计算模块是根据直射天线(RHCP天线)和反射天线(LHCP天线)的搭载位置,由于直射天线和反射天线在星载场景下搭载的位置不同,直射信号在卫星朝上的位置,也就是说离发射机更近,根据直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,也称为天线位置误差ρant,并计算直射天线所在平面(与星下点地球切面的平行表面)和反射天线所在平面(与星下点地球切面的平行表面)的高度差,也称为天线高度差Hant,并将天线位置误差和天线高度差的计算结果传给海面高度要素计算模块。
所述海面高度要素计算模块是根据之前计算的镜面反射点位置、海面高度初步计算值、电离层误差、对流层误差、天线位置误差和天线高度差参数,计算镜面反射点所在位置的海面(或湖面、以及海冰面)高度H,所述海面高度为海面到地球椭球面的绝对高度。
本实施例中,在轨运行的所述GNSS卫星,包括GPS卫星、北斗卫星、伽利略卫星、格洛纳斯卫星以及QZSS卫星等。
所述发射机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面或湖面或海冰面反射回来的GNSS卫星的信号。其中,所述海面高度要素计算模块,进一步包括:
海面相对高度差计算子模块,用于计算GNSS卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure BDA0002202535880000061
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
海面高度计算子模块,根据GNSS卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,用于计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。如图2,本实施例同时提供了一种星载GNSS-R海面高度要素定标方法,包括如下步骤:
S1:根据发射机位置、接收机位置及其相对应的几何关系,采用线段二分或是角度二分的算法模型,计算GNSS直射信号在海面进行镜面反射的反射点位置;其中,所述发射机为GNSS卫星;所述接收机为GNSS-R卫星载荷;
所述发射机的位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,可根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
S2:根据接收机接收到的GNSS直射信号以及从海面反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal
S3:根据接收机的高度,采用电离层经验模型或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差ρion;所述GNSS信号包括GNSS直射信号和反射信号;
S4:根据接收机的高度,采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差ρtro
S5:根据GNSS直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,即天线位置误差ρant;并计算直射天线所在平面和反射天线所在平面的高度差,即天线高度差Hant
S6:根据步骤S1-S5计算出的镜面反射点位置、海面高度的初步计算结果、电离层延迟误差、对流层延迟误差、天线位置误差和天线高度差,来计算镜面反射点所在位置的海面高度H,所述海面高度为海面到地球椭球面的绝对高度。
其中,所述步骤S6进一步包括:
S61:计算GNSS卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure BDA0002202535880000081
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
S62:再根据GNSS卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。
其中,所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
其中,所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面或湖面或海冰面反射回来的GNSS卫星的信号。
需要说明的是,本实施例装置和方法中涉及的所有模型计算,均为本领域常见的应用,且涉及知识过多,故本发明不再对其进一步展开进行阐述。
本发明的装置和方法不仅适用于海面的情况,同样适用于湖面或海冰面的情况,本发明不对其做出限定。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。

Claims (8)

1.一种星载GNSS-R海面高度要素定标装置,其特征在于,包括:镜面反射点计算模块、海面高度初步计算模块、电离层误差计算模块、对流层误差计算模块、天线位置误差计算模块、海面高度要素计算模块,其中:
所述镜面反射点计算模块,根据发射机位置、接收机位置及其相对应的几何关系,采用线段二分或是角度二分的算法模型,来计算GNSS直射信号在海面进行镜面反射的反射点位置,并将镜面反射点位置传给所述海面高度要素计算模块;其中,所述发射机为GNSS卫星;所述接收机为GNSS-R卫星载荷;
所述海面高度初步计算模块,根据接收机接收到的GNSS直射信号以及从海面反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal;以此作为海面高度的初步计算结果,并传给所述海面高度要素计算模块;
所述电离层误差计算模块,根据接收机的高度,采用电离层经验模型或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差ρion,并将电离层延迟误差的计算结果传给所述海面高度要素计算模块;所述GNSS信号包括GNSS直射信号和反射信号;
所述对流层误差计算模块,根据接收机的高度,采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差ρtro,并将对流层延迟误差的计算结果传给所述海面高度要素计算模块;
所述天线位置误差计算模块,根据直射天线和反射天线的搭载位置,计算天线位置误差ρant和天线高度差Hant,并将二者的计算结果传给海面高度要素计算模块;其中:
根据GNSS直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,也称为天线位置误差ρant,并计算直射天线所在平面和反射天线所在平面的高度差,也称为天线高度差Hant,并将天线位置误差和天线高度差的计算结果传给海面高度要素计算模块;
所述海面高度要素计算模块,用于根据之前计算的镜面反射点位置、海面高度的初步计算结果、电离层延迟误差、对流层延迟误差、天线位置误差和天线高度差,来计算镜面反射点所在位置的海面高度H,所述海面高度为海面到地球椭球面的绝对高度;
所述海面高度要素计算模块,进一步包括:
海面相对高度差计算子模块,用于计算GNSS-R卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure FDA0004163090840000021
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
海面高度计算子模块,根据GNSS-R卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,用于计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。
2.如权利要求1所述的装置,其特征在于,
所述发射机位置,根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
3.如权利要求1所述的装置,其特征在于,所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
4.如权利要求3所述的装置,其特征在于,
所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面反射回来的GNSS卫星的信号。
5.一种星载GNSS-R海面高度要素定标方法,其特征在于,包括如下步骤:
S1:根据发射机位置、接收机位置及其相对应的几何关系,采用线段二分或是角度二分的算法模型,计算GNSS直射信号在海面进行镜面反射的反射点位置;其中,所述发射机为GNSS卫星;所述接收机为GNSS-R卫星载荷;
S2:根据接收机接收到的GNSS直射信号以及从海面反射回来的GNSS反射信号,采用码相位测高模型、载波相位测高模型、或是干涉测高模型计算反射信号与直射信号的延迟差ρsignal
S3:根据接收机的高度,采用电离层经验模型或是双频改正模型,计算GNSS信号在传播中的电离层延迟误差ρion;所述GNSS信号包括GNSS直射信号和反射信号;
S4:根据接收机的高度,采用改进的Saastamoinen模型和改进Hopfield模型,计算GNSS信号在传播中的对流层延迟误差ρtro
S5:根据GNSS直射信号从发射机传播到直射天线和反射天线的不同路径,计算直射信号分别到直射天线和反射天线的传播距离之差,即天线位置误差ρant;并计算直射天线所在平面和反射天线所在平面的高度差,即天线高度差Hant
S6:根据步骤S1-S5计算出的镜面反射点位置、海面高度的初步计算结果、电离层延迟误差、对流层延迟误差、天线位置误差和天线高度差,来计算镜面反射点所在位置的海面高度H,所述海面高度为海面到地球椭球面的绝对高度;
所述步骤S6进一步包括:
S61:计算GNSS-R卫星载荷的反射天线所在平面以及镜面反射点所在平面的海面相对高度差ΔH,具体为:
Figure FDA0004163090840000031
其中,θ为对应GNSS卫星的卫星仰角、ρtro为对流层延迟误差、ρsignal为海面高度的初步计算结果、ρion为电离层延迟误差、ρant为天线位置误差、Hant为天线高度差;
S62:再根据GNSS-R卫星载荷的反射天线所在平面相对于地球椭球面的绝对高度H0以及ΔH,计算海面到地球椭球面的绝对高度,即为“镜面反射点所在位置的海面高度”,计算公式如下:
H=H0+ΔH,海面在地球椭球面之上;
H=H0-ΔH,海面在地球椭球面之下。
6.如权利要求5所述的方法,其特征在于,
所述发射机的位置,根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号计算,或是根据GNSS-R卫星载荷的反射天线接收GNSS反射信号的时间以及GNSS卫星的精密星历计算;
所述接收机位置,根据GNSS-R卫星载荷的直射天线接收到的GNSS直射信号,计算得到接收机的直射天线的位置,以此代表接收机的位置。
7.如权利要求5所述的方法,其特征在于,所述GNSS-R卫星载荷的接收天线以下产生的对流层延迟,以及GNSS-R卫星载荷的接收天线以上的对流层所产生的延迟,对于直射信号和反射信号的影响是相同的。
8.如权利要求7所述的方法,其特征在于,所述GNSS-R卫星载荷的接收天线分为两类:
直射天线:朝天放置,接收从上方来的直射信号,也就是GNSS卫星直接发射的信号;
反射天线:朝地放置,接收从下方来的反射信号,也就是经过海面反射回来的GNSS卫星的信号。
CN201910869845.7A 2019-09-16 2019-09-16 一种星载gnss-r海面高度要素定标装置及方法 Active CN110456388B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910869845.7A CN110456388B (zh) 2019-09-16 2019-09-16 一种星载gnss-r海面高度要素定标装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910869845.7A CN110456388B (zh) 2019-09-16 2019-09-16 一种星载gnss-r海面高度要素定标装置及方法

Publications (2)

Publication Number Publication Date
CN110456388A CN110456388A (zh) 2019-11-15
CN110456388B true CN110456388B (zh) 2023-06-09

Family

ID=68491974

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910869845.7A Active CN110456388B (zh) 2019-09-16 2019-09-16 一种星载gnss-r海面高度要素定标装置及方法

Country Status (1)

Country Link
CN (1) CN110456388B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113296135B (zh) * 2020-02-24 2024-03-29 香港理工大学深圳研究院 形变监测方法、装置和接收机
CN113031014B (zh) * 2021-03-05 2023-11-03 北京航空航天大学 一种基于gnss反射信号成像的海面流冰检测系统
CN113625312B (zh) * 2021-06-28 2024-08-30 中国空间技术研究院 一种gps-r/bds-r反射延迟海况偏差的量化与预测方法及系统
CN116088011A (zh) * 2023-02-01 2023-05-09 中国科学院国家空间科学中心 Gnss-r海面测量高度的快速定标方法、装置和电子设备
CN117607908B (zh) * 2024-01-23 2024-04-05 山东大学 无人机载平台的gnss-r码相位海面测高方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ATE514102T1 (de) * 2008-07-11 2011-07-15 Europ Agence Spatiale Altimetrieverfahren und -system
US7808425B2 (en) * 2008-09-23 2010-10-05 Agence Spatiale Europeenne Space-borne altimetry apparatus, antenna subsystem for such an apparatus and methods for calibrating the same
CN101887128B (zh) * 2010-07-09 2012-11-14 中国科学院测量与地球物理研究所 确定全球卫星导航系统导航卫星频间偏差的方法
CN107037460B (zh) * 2017-05-27 2019-07-02 金华航大北斗应用技术有限公司 一种用于海面测高的北斗meo卫星机载反射信号建模方法

Also Published As

Publication number Publication date
CN110456388A (zh) 2019-11-15

Similar Documents

Publication Publication Date Title
CN110456388B (zh) 一种星载gnss-r海面高度要素定标装置及方法
US10234564B2 (en) System and method for determining the direction of a false GNSS satellite signal transmitter
US7292185B2 (en) Attitude determination exploiting geometry constraints
US6397147B1 (en) Relative GPS positioning using a single GPS receiver with internally generated differential correction terms
RU2566685C2 (ru) Система и способ определения курсового угла транспортного средства
US8242953B2 (en) Distance dependent error mitigation in real-time kinematic (RTK) positioning
AU2008260578B2 (en) Distance dependant error mitigation in real-time kinematic (RTK) positioning
CN109613585A (zh) 一种对基站天线超短基线gnss双天线实时测向的方法
US20150048973A1 (en) System and method for detecting false global navigation satellite system satellite signals
RU2727584C1 (ru) Альтиметрический способ выполнения морской высотометрии с определением уклонения отвесной линии и устройство для его осуществления
US11693120B2 (en) System and method for providing GNSS corrections
KR101631967B1 (ko) 장주기 파고측정용 gnss 기반 부이에서의 수신기 위상말림 보상을 위한 ppp 및 마그네틱 컴파스 통합시스템
US20110090114A1 (en) Ultra short baseline gnss receiver
AU2021302454B2 (en) Method for determining a state parameter of a receiver and an apparatus for performing such a method
CN111290004A (zh) 伪距差分定位方法、装置、电子设备及存储介质
Yan et al. Performance analysis of oceanographic research vessel precise point positioning based on BDS/GNSS RTK receivers
KR102031838B1 (ko) 위성항법 시스템의 보정정보를 처리하는 장치 및 방법
JP3598372B2 (ja) 測位用衛星受信機バイアスの測定方法及び装置
CN113359204A (zh) 一种基于gnss-r伪随机噪声码延迟的水下重力场反演方法及系统
Liu Positioning performance of single-frequency GNSS receiver using Australian regional ionospheric corrections
US7259717B2 (en) Method and device for determining the relative position of two points
CN115166788B (zh) 一种低仰角对流层误差修正方法、系统、设备及介质
CN117310770B (zh) 基于北斗无电离层组合精密单点定位的海平面高反演方法
US20200386896A1 (en) Single-epoch pseudo-range positioning under varying ionosphere delays
CN211698218U (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