CN113049062B - 一种gnss直反射信号载波干涉测量湖泊水位的装置 - Google Patents
一种gnss直反射信号载波干涉测量湖泊水位的装置 Download PDFInfo
- Publication number
- CN113049062B CN113049062B CN202110270801.XA CN202110270801A CN113049062B CN 113049062 B CN113049062 B CN 113049062B CN 202110270801 A CN202110270801 A CN 202110270801A CN 113049062 B CN113049062 B CN 113049062B
- Authority
- CN
- China
- Prior art keywords
- signals
- module
- direct
- lake
- water level
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000001427 coherent effect Effects 0.000 claims description 26
- 238000012544 monitoring process Methods 0.000 claims description 17
- 238000005305 interferometry Methods 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 13
- 238000009825 accumulation Methods 0.000 claims description 11
- 230000001934 delay Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 101710195281 Chlorophyll a-b binding protein Proteins 0.000 description 2
- 101710143415 Chlorophyll a-b binding protein 1, chloroplastic Proteins 0.000 description 2
- 101710181042 Chlorophyll a-b binding protein 1A, chloroplastic Proteins 0.000 description 2
- 101710091905 Chlorophyll a-b binding protein 2, chloroplastic Proteins 0.000 description 2
- 101710095244 Chlorophyll a-b binding protein 3, chloroplastic Proteins 0.000 description 2
- 101710127489 Chlorophyll a-b binding protein of LHCII type 1 Proteins 0.000 description 2
- 101710184917 Chlorophyll a-b binding protein of LHCII type I, chloroplastic Proteins 0.000 description 2
- 101710102593 Chlorophyll a-b binding protein, chloroplastic Proteins 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
- G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
- G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
- G01F23/284—Electromagnetic waves
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/33—Multimode operation in different systems which transmit time stamped messages, e.g. GPS/GLONASS
-
- 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
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Thermal Sciences (AREA)
- Fluid Mechanics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种GNSS直反射信号载波干涉测量湖泊水位的装置。该装置中,双频右旋圆极化天线接收导航系统中N颗卫星的直射信号;导航系统包括GPS卫星导航系统和北斗卫星导航系统;双频左旋圆极化天线接收经待测湖泊的湖面反射的导航系统的反射信号;双模导航模块产生定位信息;四通道射频前端对直射信号和反射信号进行转换得到四路数字中频信号;多通道直/反协同处理模块根据定位信息对四路数字中频信号进行干涉处理,得到反射信号相对于直射信号的载波相位时延;湖泊水位反演模块用于根据载波相位时延和定位信息进行反演得到待测湖泊的水位。本发明应用于多个不同高度角的卫星信号,从而得到湖泊水位,测量准确度高,适应范围广。
Description
技术领域
本发明涉及湖泊水位测量领域,特别是涉及一种GNSS直反射信号载波干涉测量湖泊水位的装置。
背景技术
利用全球导航卫星系统(Global Navigation Satellites System,GNSS)反射信号进行地球表面物理参数的反演是遥感领域的新型技术之一,具有信号源广、成本低、探测设备轻等优势。该技术通过星载、机载或岸基的特殊装置接收、处理GNSS反射信号可进行地表物理参数的探测。
湖泊水位作为重要的水文参数,是内陆遥感领域长期关注的对象之一。GNSS反射信号高度计最初作为传统高度计的替代技术被提出,通过连续接收处理GNSS直射信号和经地表反射的回波信号得到反射信号相对于直射信号的路径时延,并利用GNSS卫星、接收平台以及镜面反射点之间的几何关系实现接收平台高度的反演,最后与参考高度相减得到相对高度测量。
通常的GNSS反射信号高度计都是利用伪码时延进行高度测量。基于伪码的方法最大的不足是测量精度低,难以满足湖泊水位的测量。鉴于湖泊水面粗糙度低,经湖泊表面反射的GNSS为相干信号,可以充分利于载波相位进行时延测量。目前,基于载波相位测量反射信号相对于直射信号时延的方法是单天线干涉法,通过利用反射信号对直射信号的干扰形成的振荡进行高度测量,其最大缺点是时间分辨率低,且只能应用低高度角的卫星信号。
发明内容
基于此,有必要提供一种GNSS直反射信号载波干涉测量湖泊水位的装置,以应用于多个不同高度角的卫星信号,从而得到湖泊水位。
为实现上述目的,本发明提供了如下方案:
一种GNSS直反射信号载波干涉测量湖泊水位的装置,包括:双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块、多通道直/反协同处理模块和湖泊水位反演模块;
所述双频右旋圆极化天线用于接收导航系统中N颗卫星的直射信号;所述导航系统包括GPS卫星导航系统和北斗卫星导航系统;所述双频左旋圆极化天线用于接收经待测湖泊的湖面反射的所述导航系统的反射信号;所述双模导航模块与所述双频右旋圆极化天线连接;所述双模导航模块用于接收所述直射信号,并产生定位信息;所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接;所述四通道射频前端用于对所述直射信号和所述反射信号进行转换,得到四路数字中频信号;所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接;所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延;所述湖泊水位反演模块分别与所述多通道直/反协同处理模块和所述双模导航模块连接;所述湖泊水位反演模块用于根据所述载波相位时延和所述定位信息进行反演得到所述待测湖泊的水位。
可选的,所述定位信息包括伪随机噪声号、伪距、载波频率、方位角、高度角、纬度、经度和高度。
可选的,所述GNSS直反射信号载波干涉测量湖泊水位的装置,还包括:与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块;
所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示,并当所述待测湖泊的水位超过预设水位时进行告警。
可选的,所述GNSS直反射信号载波干涉测量湖泊水位的装置,还包括:与所述湖泊水位监控模块连接的信息管理服务器;
所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。
可选的,所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块;
所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号,得到四路可见数字中频信号;
所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延。
可选的,所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块;
所述频谱估计模块用于根据所述高度角,采用谱分析方法对N个所述载波相位时延进行处理,得到N个频谱;
所述相对高度反演模块用于搜索N个所述频谱的峰值,并根据N个所述峰值处的频率反演得到N个相对水位高度;
所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
可选的,所述四通道射频前端包括GPS射频前端和北斗射频前端;
所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理,得到两路GPS数字中频信号;所述数字化处理包括下变频、滤波、增益控制以及采样量化;
所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理,得到两路北斗数字中频信号;两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
可选的,所述四通道直/反射信号干涉处理模块,包括:载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块;
所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波;
所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码;
所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘,得到两路载波剥离直射信号,并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘,得到两路载波剥离反射信号;
所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘,将两路所述载波剥离反射信号分别与所述本地伪码相乘,得到四路干涉信号;
所述相干积分器用于对所述干涉信号进行相干积分;
所述非相干累加模块用于对相干积分后的复数进行多次非相干累加,得到相关功率值;
所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
可选的,所述湖泊水位监控模块与所述多通道直/反协同处理模块连接;所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
可选的,所述信息管理服务器与IP访问端口连接。
与现有技术相比,本发明的有益效果是:
本发明提出了一种GNSS直反射信号载波干涉测量湖泊水位的装置。该装置通过双频右旋圆极化(Right-Handed Circular Polarization,RHCP)天线和双频左旋圆极化(Left-Handed Circular Polarization,LCHP)天线两个天线对同一观测区域的GPS/北斗系统的多颗卫星的直射/反射信号进行接收、处理,从而测量湖泊水位,其中,多通道直/反协同处理模块利用N颗卫星的直射信号和反射信号进行干涉处理得到反射信号相对于直射信号的载波相位时延,湖泊水位反演模块通过N个载波相位时延构成的集合反演得到待测湖泊的水位。该装置不限于应用在低高度角的卫星信号,可以用于各种不同高度角的卫星信号,从而得到湖泊水位,适应范围更广;该装置利用载波相位时延进行湖泊水平的测量,测量精度高;该装置仅是信号接收装置,装置结构简单、成本低、功耗低;该装置利用GNSS信号作为信号源,可实施全天候观测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的GNSS直反射信号载波干涉测量湖泊水位的装置的结构示意图;
图2为本发明实施例提供的多通道直/反协同处理模块的框图;
图3为本发明实施例提供的四通道直/反射信号干涉处理模块的框图;
图4为本发明实施例提供的湖泊水位反演模块的框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明利用N(N=n+m;n为反射天线可见的GPS卫星个数,m为反射天线可见的北斗卫星个数)颗导航卫星反射信号通过干涉处理得到n+m载波相位,然后对n+m载波相位进行湖泊水位测测量。
图1为本发明实施例提供的GNSS直反射信号载波干涉测量湖泊水位的装置的结构示意图。参见图1,本实施例的GNSS直反射信号载波干涉测量湖泊水位的装置,包括:双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块(导航模块)、多通道直/反协同处理模块和湖泊水位反演模块。
本实施例中,所述双频右旋圆极化天线为GPS L1/北斗B1双频右旋圆极化天线(RHCP);所述双频左旋圆极化天线为GPS L1/北斗B1双频左旋圆极化天线。所述双频右旋圆极化天线对天安置,用于接收导航系统中N颗卫星的直射信号(射频电压信号);所述导航系统包括GPS卫星导航系统L1和北斗卫星导航系统B1;所述双频左旋圆极化天线面向海面安置,用于接收经待测湖泊的湖面反射的所述导航系统的反射信号(射频电压信号)。
具体的,双频右旋圆极化天线用于接收GPS L1/北斗B1直射电磁信号,并将该电磁信号转换为射频电压信号Sd(t),即:
其中,Adi为第i颗卫星GNSS信号幅度;Di(t)表示第i颗卫星GNSS信号调制的数据;Ci(t)为GNSS信号伪随机码;fi为第i颗卫星GNSS信号的载波频率;为第i颗卫星GNSS信号初始相位,t表示时间。
双频左旋圆极化天线用于接收经湖面反射的GPS L1/北斗B1电磁信号,并将该电磁信号转换为射频电压信号Sr(t),即:
本实施例中,所述双模导航模块为GPS L1/北斗B1双模导航模块。所述双模导航模块与所述双频右旋圆极化天线连接。所述双模导航模块用于接收所述直射信号以完成导航定位,从而产生定位信息;所述定位信息包括伪随机噪声号(pseudo random noise,prn)、伪距、载波频率、高度角和方位角{prni,ρi,fi,θi,φi},以及纬度、经度和高度{lat,lon,h},其中,i表示第i颗卫星GNSS信号。其中,双模导航模块可以采用成熟的导航SOC芯片(Systemon Chip)。
本实施例中,所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接。所述四通道射频前端用于对所述直射信号和所述反射信号进行转换,得到四路数字中频信号。其中,四通道射频前端中两个通道连接双频右旋圆极化天线,另外两个通道连接双频左旋圆极化天线,输出四路数字中频信号给直/反信号协同处理模块完成信号处理和载波相位提取。
在实际应用中,所述四通道射频前端包括GPS射频前端和北斗射频前端,两个子射频前端由两个结构完全一样的射频通道组成,分别对GPS L1/北斗B1直射、反射信号进行数字化处理。其中,所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理,得到两路GPS数字中频信号;所述数字化处理包括下变频、滤波、增益控制以及采样量化。所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理,得到两路北斗数字中频信号;两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
具体的,四通道射频前端接收RHCP天线和LHCP天线输出的射频信号Sd(t)和Sr(t),输出四路数字中频信号,分别对应GPS直射数字中频信号、GPS反射数字中频信号、北斗直射数字中频信号和北斗反射数字中频信号。GPS和北斗信号处理完全一样,为了描述方便,对中频信号不做GPS和北斗符合上的区分,假设直射中频信号为Sd(n),反射中频信号为Sr(n),即:
其中,n为自然数;Ts为采样率。
本实施例中,所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接;所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延。其中,所述多通道直/反协同处理模块首先接收定位信息选择直/反射信号协同处理的N路卫星信号,然后对四路所述数字中频信号进行四通道直/反干涉处理得到N个载波相位时延。
实际应用中,所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块(4通道直/反干涉模块),如图2所示。所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号,得到四路可见数字中频信号。所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延。
如图3所示,所述四通道直/反射信号干涉处理模块,包括:载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块。所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波。所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码。所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘,得到两路载波剥离直射信号,并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘,得到两路载波剥离反射信号。所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘,将两路所述载波剥离反射信号分别与所述本地伪码相乘,得到四路干涉信号。所述相干积分器用于对所述干涉信号进行相干积分。所述相干积分器为1ms相干积分器。所述非相干累加模块用于对相干积分后的复数进行多次非相干累加,得到相关功率值。所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
具体的,卫星选择模块接收GPS L1/北斗B1导航模块输出的伪随机噪声号、高度角和方位角{prni,θi,φi},以及接收高度角范围信息(截止信息)和方位角范围信息{θmask,φmask},从而对GPS L1/北斗B1双频左旋圆极化天线视野范围内的GPS/北斗卫星进行预测,输出双频左旋圆极化天线视野范围内的卫星伪随机噪声号组成的集合{prn1,prn2,…,prnN}。预判决准则为:
θmin≤θi≤θmax,θmin,θmax∈θmask
φmin≤φi≤φmax,φmin,φmax∈φmask。
如图3所示,载波发生器根据卫星选择模块输出的载波频率fi产生本地载波,即:
本地码发生器根据卫星选择模块输出的伪随机噪声号和伪距信息{prni,ρi}产生本地伪码,即:
Scode_i(n)=Ci(nTS);
载波剥离模块将直/反射信号和本地载波进行相乘,即:
4路干涉处理模块首先利用载波剥离和伪码剥离后的直/反射信号产生4路干涉信号,即:
其次对4路干涉信号进行1ms相干积分,根据伪码特性,可得4路1ms相干积分后的信号,分别为第1路1ms相干积分后的信号复数相关值Icoh_i_1(m)、第2路1ms相干积分后的信号复数相关值Icoh_i_2(m)、第3路1ms相干积分后的信号复数相关值Icoh_i_3(m)和第4路1ms相干积分后的信号复数相关值Icoh_i_4(m),即:
其中,Ns表示采样点总数。然后对1ms相干积分后的复数相关值进行M次非相干累加,输出4路相关功率值,分别为第1路相关功率值Ii1、第2路相关功率值Ii2、第3路相关功率值Ii3和第4路相关功率值Ii4,即:
相位计算模块根据M次非相干累加输出的4路相关功率值估计反射信号相对于直射信号的载波相位时延,即:
如图4所示,所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块。所述频谱估计模块用于根据所述高度角,采用谱分析方法对N个所述载波相位时延进行处理,得到N个频谱。所述相对高度反演模块用于搜索N个所述频谱的峰值,并根据N个所述峰值处的频率反演得到N个相对水位高度。所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
具体的,首先,根据在一定时间段内相位多通道直/反信号协同处理模块输出的N个载波相位时延的时间序列得到个N相位时序的频谱。本实施例采用Lomb-Scargle谱估计方法,估计过程中时间轴由时间转化为高度角的正弦函数t=sin[θi(t)],具体过程表示为:
然后,搜索N个Lomb-Scargle谱的峰值位置{fLS_1,fLS_2,…,fLS_N},并根据频谱峰值频率与高度关系反演得到接收天线相对于水位的高度,即:
其中,λ1,λ2…λN表示第1,2…,N颗GNSS卫星信号的波长。
最后,水位计算模块根据反演得到的相对水位高度{hr1,hr2,…,hrN}和参考面高度href计算水位高度,得到待测湖泊的水位,即:
在实际应用中,所述GNSS直反射信号载波干涉测量湖泊水位的装置,还包括:与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块;所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示,并当所述待测湖泊的水位超过预设水位时进行告警。所述湖泊水位监控模块还可以与所述多通道直/反协同处理模块连接;所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
湖泊水位监控模块提供了湖泊水位测量与用户的交互界面,接收湖泊水位反演模块的待测湖泊的水位Δh、GPS L1/北斗B1双模导航模块的定位信息{lat,lon,h}进行水位和监测站信息的显示和存储,并当水位超过预设水位时产生告警。此外根据LHCP天线的架设场景,可以设置观测视场的截止范围{θmask,φmask}。
在实际应用中,所述GNSS直反射信号载波干涉测量湖泊水位的装置,还包括:与所述湖泊水位监控模块连接的信息管理服务器;所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。所述信息管理服务器还可以与IP访问端口连接。本实施例设置信息管理服务器,实现了对水位、定位、告警、日志以及用户信息的存储和管理,并能接收IP接口的访问。
本实施例通过双频右旋圆极化天线和双频左旋圆极化天线两个天线对同一观测区域的GPS/北斗系统的多颗卫星的直射/反射信号进行接收、处理,从而测量湖泊水位。首先在每个积分周期内,利用4通道干涉法估计N颗卫星反射信号相对于直射信号的载波相位时延,即;然后利用载波相位反演得到接收平台的高度,并与参考面高度相减得到湖面高度(待测湖泊的水位)。本实施例适应范围广;利用载波相位进行湖泊水平的测量,测量精度高;装置结构简单、成本低、功耗低;利用GNSS信号作为信号源,可实施全天候观测。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,包括:双频右旋圆极化天线、双频左旋圆极化天线、四通道射频前端、双模导航模块、多通道直/反协同处理模块和湖泊水位反演模块;
所述双频右旋圆极化天线用于接收导航系统中N颗卫星的直射信号;所述导航系统包括GPS卫星导航系统和北斗卫星导航系统;所述双频左旋圆极化天线用于接收经待测湖泊的湖面反射的所述导航系统的反射信号;所述双模导航模块与所述双频右旋圆极化天线连接;所述双模导航模块用于接收所述直射信号,并产生定位信息;所述四通道射频前端与所述双频右旋圆极化天线和所述双频左旋圆极化天线连接;所述四通道射频前端用于对所述直射信号和所述反射信号进行转换,得到四路数字中频信号;所述多通道直/反协同处理模块分别与所述四通道射频前端和所述双模导航模块连接;所述多通道直/反协同处理模块用于根据所述定位信息对四路所述数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延;所述湖泊水位反演模块分别与所述多通道直/反协同处理模块和所述双模导航模块连接;所述湖泊水位反演模块用于根据所述载波相位时延和所述定位信息进行反演得到所述待测湖泊的水位;
所述定位信息包括伪随机噪声号、伪距、载波频率、方位角、高度角、纬度、经度和高度;
所述多通道直/反协同处理模块包括卫星选择模块和N个四通道直/反射信号干涉处理模块;
所述卫星选择模块用于根据所述伪随机噪声号、所述伪距、所述载波频率、所述方位角、所述高度角、方位角范围和高度角范围从四路所述数字中频信号中选择镜面反射点对双频左旋圆极化天线可见的卫星信号,得到四路可见数字中频信号;
所述四通道直/反射信号干涉处理模块用于对四路所述可见数字中频信号进行干涉处理,得到所述反射信号相对于所述直射信号的载波相位时延;
所述四通道直/反射信号干涉处理模块,包括:载波发生器、本地码发生器、载波剥离模块、伪码剥离模块、相干积分器、非相干累加模块和相位计算模块;
所述载波发生器用于根据所述卫星选择模块输出的所述载波频率产生本地载波;
所述本地码发生器用于根据所述卫星选择模块输出的所述伪随机噪声号和所述伪距产生本地伪码;
所述载波剥离模块用于将四路所述可见数字中频信号中的直射信号分别与所述本地载波相乘,得到两路载波剥离直射信号,并将四路所述可见数字中频信号中的反射信号分别与所述本地载波相乘,得到两路载波剥离反射信号;
所述伪码剥离模块用于将两路所述载波剥离直射信号分别与所述本地伪码相乘,将两路所述载波剥离反射信号分别与所述本地伪码相乘,得到四路干涉信号;
所述相干积分器用于对所述干涉信号进行相干积分;
所述非相干累加模块用于对相干积分后的复数进行多次非相干累加,得到相关功率值;
所述相位计算模块用于根据所述相关功率值计算所述反射信号相对于所述直射信号的载波相位时延。
2.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,还包括:与所述湖泊水位反演模块和所述双模导航模块连接的湖泊水位监控模块;
所述湖泊水位监控模块用于对所述定位信息和所述待测湖泊的水位进行可视化显示,并当所述待测湖泊的水位超过预设水位时进行告警。
3.根据权利要求2所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,还包括:与所述湖泊水位监控模块连接的信息管理服务器;
所述信息管理服务器用于对所述定位信息、所述待测湖泊的水位和所述告警进行存储。
4.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,所述湖泊水位反演模块包括频谱估计模块、相对高度反演模块和水位计算模块;
所述频谱估计模块用于根据所述高度角,采用谱分析方法对N个所述载波相位时延进行处理,得到N个频谱;
所述相对高度反演模块用于搜索N个所述频谱的峰值,并根据N个所述峰值处的频率反演得到N个相对水位高度;
所述水位计算模块用于根据所述相对水位高度和参考面高度计算所述待测湖泊的水位。
5.根据权利要求1所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,所述四通道射频前端包括GPS射频前端和北斗射频前端;
所述GPS射频前端用于对所述GPS卫星导航系统中卫星的直射信号和所述GPS卫星导航系统中卫星的反射信号进行数字化处理,得到两路GPS数字中频信号;所述数字化处理包括下变频、滤波、增益控制以及采样量化;
所述北斗射频前端用于对所述北斗卫星导航系统中卫星的直射信号和所述北斗卫星导航系统中卫星的反射信号进行数字化处理,得到两路北斗数字中频信号;两路所述GPS数字中频信号和两路所述北斗数字中频信号构成四路所述数字中频信号。
6.根据权利要求2所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,所述湖泊水位监控模块与所述多通道直/反协同处理模块连接;所述多通道直/反协同处理模块用于接收所述湖泊水位监控模块发送的方位角范围和高度角范围。
7.根据权利要求3所述的一种GNSS直反射信号载波干涉测量湖泊水位的装置,其特征在于,所述信息管理服务器与IP访问端口连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110270801.XA CN113049062B (zh) | 2021-03-12 | 2021-03-12 | 一种gnss直反射信号载波干涉测量湖泊水位的装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110270801.XA CN113049062B (zh) | 2021-03-12 | 2021-03-12 | 一种gnss直反射信号载波干涉测量湖泊水位的装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113049062A CN113049062A (zh) | 2021-06-29 |
CN113049062B true CN113049062B (zh) | 2022-04-15 |
Family
ID=76511986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110270801.XA Active CN113049062B (zh) | 2021-03-12 | 2021-03-12 | 一种gnss直反射信号载波干涉测量湖泊水位的装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113049062B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113671210B (zh) * | 2021-08-19 | 2022-07-26 | 北京航空航天大学 | 陆基gnss直反射信号载波干涉测量河流流速的装置 |
CN113671552B (zh) * | 2021-08-24 | 2024-02-06 | 北京航空航天大学 | 一种陆基gnss反射计的河流流速与水位反演方法 |
CN114355421B (zh) * | 2021-12-22 | 2023-08-18 | 杭州电子科技大学 | 一种基于北斗卫星l4和cmc组合观测值洪水探测方法 |
CN114137575B (zh) * | 2022-02-08 | 2022-05-10 | 浙江国遥地理信息技术有限公司 | 顾及卫星偏差和载噪比弧段影响的洪水探测方法 |
CN115165029B (zh) * | 2022-07-01 | 2023-03-28 | 中国水利水电科学研究院 | 基于双天线测姿的北斗水位监测装置、方法及系统 |
US11719732B1 (en) * | 2022-07-25 | 2023-08-08 | Divirod, Inc. | Reflectometer sensor |
CN116699658A (zh) * | 2023-08-04 | 2023-09-05 | 山东大学 | 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608915A (zh) * | 2008-06-17 | 2009-12-23 | 厦门瀛寰电子科技有限公司 | 全新的海浪观测装置 |
CN104677463A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于统计量的gnss-r多源海面测高装置 |
CN104678388A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种用于提高测高精度的gnss-r正交化交错处理装置 |
CN104678371A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于时延修正的海面高度测量装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6549165B2 (en) * | 2001-01-19 | 2003-04-15 | Agence Spatiale Europeenne | Ocean altimetry interferometric method and device using GNSS signals |
-
2021
- 2021-03-12 CN CN202110270801.XA patent/CN113049062B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101608915A (zh) * | 2008-06-17 | 2009-12-23 | 厦门瀛寰电子科技有限公司 | 全新的海浪观测装置 |
CN104677463A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于统计量的gnss-r多源海面测高装置 |
CN104678388A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种用于提高测高精度的gnss-r正交化交错处理装置 |
CN104678371A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于时延修正的海面高度测量装置 |
Non-Patent Citations (2)
Title |
---|
基于GNSS卫星反射信号的海冰厚度探测;高洪兴 等;《电子与信息学报》;20170531;第39卷(第5期);第1096-1100页 * |
机载导航卫星反射信号海面高度测量研究;王峰 等;《航空学报》;20210127;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113049062A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113049062B (zh) | 一种gnss直反射信号载波干涉测量湖泊水位的装置 | |
CN113049777B (zh) | 一种gnss直反射信号载波干涉测量土壤湿度的装置 | |
US6549165B2 (en) | Ocean altimetry interferometric method and device using GNSS signals | |
US8411726B2 (en) | Method and apparatus for performing signal correlation at multiple resolutions to mitigate multipath interference | |
Zavorotny et al. | Scattering of GPS signals from the ocean with wind remote sensing application | |
Zhang et al. | Feasibility of code-level altimetry using coastal BeiDou reflection (BeiDou-R) setups | |
CN104677463A (zh) | 一种基于统计量的gnss-r多源海面测高装置 | |
CN113805208A (zh) | 一种适用于导航接收机的gnss-ir测高方法 | |
CN111521650B (zh) | 土壤湿度测量中直射信号和反射信号的相位提取方法 | |
Wang et al. | Sea-state observation using reflected BeiDou GEO signals in frequency domain | |
CN116699658A (zh) | 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法 | |
He et al. | Coastal altimetry using interferometric phase from GEO satellite in quasi-zenith satellite system | |
EP0852734A1 (fr) | Procede et dispositif de geodesie et/ou d'imagerie par traitement de signaux satellitaires | |
Li et al. | Real-time processing of reflected GNSS signals for remote sensing: System and experiments | |
Connselman III et al. | Backpack VLBI terminal with subscentimeter capability | |
Manandhar et al. | GPS reflected signal analysis using software receiver | |
CN113671210B (zh) | 陆基gnss直反射信号载波干涉测量河流流速的装置 | |
CN115144884A (zh) | 一种基于卫星反射信号的海面风速反演方法及芯片模组 | |
Hobiger et al. | Software-defined radio direct correlation GNSS reflectometry by means of GLONASS | |
CN116087997A (zh) | 基于支持北斗/gnss的电子设备的水汽实时反演方法、系统 | |
He et al. | Coastal GNSS-R phase altimetry based on the combination of L1 and L5 signals under high sea states | |
Zhang et al. | Parameter measurement of soil moisture based on GNSS-R signals | |
Meng et al. | Design of real-time GNSS-R software-defined receiver for coastal altimetry using GPS/BDS/QZSS signals | |
CN113391332B (zh) | 一种基于gnss-r的海洋参数反演及远程监测系统 | |
Yu et al. | Sea surface altimetry based on airborne GNSS signal measurements |
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 |