CN106291628A - 一种提取海面卫星导航反射信号的相干观测量的方法 - Google Patents
一种提取海面卫星导航反射信号的相干观测量的方法 Download PDFInfo
- Publication number
- CN106291628A CN106291628A CN201610979016.0A CN201610979016A CN106291628A CN 106291628 A CN106291628 A CN 106291628A CN 201610979016 A CN201610979016 A CN 201610979016A CN 106291628 A CN106291628 A CN 106291628A
- Authority
- CN
- China
- Prior art keywords
- observed quantity
- phase
- reflected signal
- sea
- parallel search
- 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
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/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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种提取海面卫星导航反射信号的相干观测量的方法。首先,将接收到反射信号的同相观测量(I)和正交相观测量(Q)进行移动平均滤波;将超过幅度阈值的观测量滤除;然后进行连续相位和相位均匀分布的并行搜索,根据所设置的相位搜索阈值剔除非相干分量,可将导航信号经过海面反射后的反射信号观测量中的非相干剔除,从而可利用相干观测量反演获得高精度的海面测高结果,特别是在海风比较大的粗糙海面也可获得高精度的反演结果。
Description
技术领域
本发明涉及一种提取海面卫星导航反射信号的相干观测量的方法,可用于全球卫星导航系统(GNSS)反射信号进行海面高度监测,属于海洋气象监测的应用领域。
背景技术
全球卫星导航系统(GNSS)经历了几十年的发展,已经在国家的社会、军事和经济生活发挥越来越重要的作用。在众多的GNSS应用领域中,利用GNSS经过反射后的导航卫星信号(GNSS‐R),对反射面的物理特性和参数进行反演已经成为导航应用领域的新的研究热点,主要应用包括海面高度、海浪、海风、海水盐度、海冰、土壤湿度、森林覆盖、移动目标探测等。
海平面高度是海洋的重要参数,也是海洋气象监测和预报领域一直的关注重点之一。利用GNSS反射信号代替传统的验潮仪已经成为一种新兴的技术手段。利用GNSS反射信号进行海面测高主要是通过布设接收机天线连续接收导航直射信号和海面反射信号,根据反射信号和直射信号之间的路径传播时延,利用GNSS卫星、接收机以及反射信号在海面的镜面反射点进行反演,从而得到接收机天线距离地球参考面的高度,将该高度减去所接收机的高度而获得海平面高度。GNSS‐R测高仪获取反射信号和直射信号之间的路径传播时延主要有两个途径,一是利用直射信号的码相位与反射信号进行相关,通过寻找相关的峰值对应的码延迟来获得时延,另外一种是直接将反射信号和直射信号的载波相位相减来获得时延。由于扩频码码长较长的特性,基于码延迟的GNSS‐R测高仪所获得的海面高度精度不高;而由于载波波长较短,基于反射信号和直射信号的载波相位的GNSS‐R测高仪可以获得较高的精度。以往的载波相位观测量的GNSS‐R测高仪采用的是直接将反射信号相位与直接信号相位相减而获取时延,由于导航信号经过海面反射后的反射信号观测量存在着非相干和相干观测量,非相干观测量影响着海面测高的精度,特别是在海风比较大的粗糙海面测高精度影响严重。
本发明提出一种提取海面卫星导航反射信号的相干观测量的方法。首先,将接收到反射信号的同相观测量(I)和正交相观测量(Q)进行移动平均滤波;将超过幅度阈值的观测量滤除;然后进行连续相位和相位均匀分布的并行搜索,根据所设置的相位搜索阈值剔除非相干分量,最后获得可用于海面测高的相干观测量。
发明内容
本发明的目的是提供一种提取海面卫星导航反射信号的相干观测量的方法。它是一种可接收和处理某
一海面区域内可视卫星的反射信号的相干观测量。首先,将接收到反射信号的同相观测量(I)和正交相观测量(Q)进行移动平均滤波;将超过幅度阈值的观测量滤除;然后进行连续相位和相位均匀分布的并行搜索,根据所设置的相位搜索阈值剔除非相干分量,最后获得可用于海面测高的相干观测量。
步骤一:将接收到反射信号的同相观测量(I)和正交相观测量(Q)通过移动平均滤波器滤除噪声,由于移动平均滤波器是一个低通滤波器,为了避免滤除有用的高频率信息,需要严格控制移动平均滤波器所用的采样个数;
步骤二:通过设置同相观测量(I)和正交相观测量(Q)所获得的振幅阈值,将低于阈值的观测量剔除,随时间变化的观测量振幅A(t)采用如下计算公式
步骤三:进行连续相位和相位均匀分布的并行搜索,随时间变化的观测量相位采用反正切公式
在并行搜索的过程中,需要遵守如下两个原则:1)随时间变化的观测量相位在四象限复平面的变化顺序只能是升序或降序,顺序可为i.e.I,II,III,IV,I,II...或IV,III,II,I,IV,III...;2)每个象限的相位采样点个数均匀分布,假设四个象限总的采样点个数为N,每个象限的采样点需保持在Nq=N/4±0.1N。
步骤四:通过连续相位和相位均匀分布的并行搜索,判断是否满足相位升序或降序、采样点是否满足上述均匀分布原则,若不满足则重新进行并行搜索;并行搜索所获得的相干观测量采样点个数需要满足测高系统时间分辨率的要求,假设测高系统的时间分辨率为T,同相观测量(I)和正交相观测量(Q)的采样率为m Hz,则相干观测量采样点个数K需满足:
K≥T·m (3);
与现有技术相比,本发明提出的提取海面卫星导航反射信号的相干观测量的方法,可将导航信号经过海面反射后的反射信号观测量中的非相干剔除,从而可利用相干观测量反演获得高精度的海面测高结果,特别是在海风比较大的粗糙海面也可获得高精度的反演结果。
附图说明
图1本发明的流程框图。
图2同相观测量(I)和正交相观测量(Q)平面的相位分布图。
图3随时间变化的相干相位图。
图4不同海面风速下测高结果的均方根误差(RMSE)图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例首先根据GNSS‐R测高系统的需求,确定时间分辨率为T,同相观测量(I)和正交相观测量(Q)的采样率为m Hz,具体步骤如下:
步骤一:将接收到反射信号的同相观测量(I)和正交相观测量(Q)通过移动平均滤波器滤除噪声,由于移动平均滤波器是一个低通滤波器,为了避免滤除有用的高频率信息,需要严格控制移动平均滤波器所用的采样个数;
步骤二:通过设置同相观测量(I)和正交相观测量(Q)所获得的振幅阈值,将低于阈值的观测量剔除,随时间变化的观测量振幅A(t)采用如下计算公式
步骤三:进行连续相位和相位均匀分布的并行搜索,随时间变化的观测量相位采用反正切公式
在并行搜索的过程中,需要遵守如下两个原则:(1)随时间变化的观测量相位在四象限复平面的变化顺序只能是升序或降序,顺序可为i.e.I,II,III,IV,I,II...或IV,III,II,I,IV,III...;(2)每个象限的相位采样点个数均匀分布,假设四个象限总的采样点个数为N,每个象限的采样点需保持在Nq=N/4±0.1N。
步骤四:通过连续相位和相位均匀分布的并行搜索,判断是否满足相位升序或降序、采样点是否满足上述均匀分布原则,若不满足则重新进行并行搜索;并行搜索所获得的相干观测量采样点个数需要满足测高系统时间分辨率的要求,假设测高系统的时间分辨率为T,同相观测量(I)和正交相观测量(Q)的采样率为m Hz,则相干观测量采样点个数K需满足:
K≥T·m (3);
如图2所示,给出2015年7月1日03:35‐04:35时间段(UTC时间)GPS PRN 9卫星的同相观测量(I)和正交相观测量(Q)平面的相位分布情况,本例测高系统时间分辨率为T=20分钟,同相观测量(I)和正交相观测量(Q)的采样率为m=0.1Hz,可以发现,每个象限均匀分布,每个象限的采样点保持10个左右,相干观测量采样点个数K>120。
如图3所示,给出2015年7月1日03:35‐04:35时间段(UTC时间)GPS PRN 9卫星随时间变化的相干相位图,可以发现随时间变化相干相位保持连续升序。
如图4所示,采用本发明提出的提取海面卫星导航反射信号的相干观测量的方法,可将导航信号经过海面反射后的反射信号观测量中的非相干剔除,利用相干观测量反演获得高精度的海面测高结果,特别是在海风为12-13m/s的粗糙海面获得2.56cm的精度。
Claims (1)
1.一种提取海面卫星导航反射信号的相干观测量的方法,其特征在于,包括以下步骤:
步骤一:将接收到反射信号的同相观测量(I)和正交相观测量(Q)通过移动平均滤波器滤除噪声,由于移动平均滤波器是一个低通滤波器,为了避免滤除有用的高频率信息,需要严格控制移动平均滤波器所用的采样个数;
步骤二:通过设置同相观测量(I)和正交相观测量(Q)所获得的振幅阈值,将低于阈值的观测量剔除,随时间变化的观测量振幅A(t)采用如下计算公式
步骤三:进行连续相位和相位均匀分布的并行搜索,随时间变化的观测量相位采用反正切公式
在并行搜索的过程中,需要遵守如下两个原则:(1)随时间变化的观测量相位在四象限复平面的变化顺序只能是升序或降序,顺序可为i.e.I,II,III,IV,I,II...或IV,III,II,I,IV,III...;(2)每个象限的相位采样点个数均匀分布,假设四个象限总的采样点个数为N,每个象限的采样点需保持在Nq=N/4±0.1N;
步骤四:通过连续相位和相位均匀分布的并行搜索,判断是否满足相位升序或降序、采样点是否满足上述均匀分布原则,若不满足则重新进行并行搜索;并行搜索所获得的相干观测量采样点个数需要满足测高系统时间分辨率的要求,假设测高系统的时间分辨率为T,同相观测量(I)和正交相观测量(Q)的采样率为m Hz,则相干观测量采样点个数K需满足:
K≥T·m (3)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610979016.0A CN106291628B (zh) | 2016-11-08 | 2016-11-08 | 一种提取海面卫星导航反射信号的相干观测量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610979016.0A CN106291628B (zh) | 2016-11-08 | 2016-11-08 | 一种提取海面卫星导航反射信号的相干观测量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291628A true CN106291628A (zh) | 2017-01-04 |
CN106291628B CN106291628B (zh) | 2019-01-08 |
Family
ID=57721394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610979016.0A Active CN106291628B (zh) | 2016-11-08 | 2016-11-08 | 一种提取海面卫星导航反射信号的相干观测量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291628B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111175741A (zh) * | 2020-01-10 | 2020-05-19 | 浙江大学 | 一种单频连续毫米波多普勒传感微波墙安全空间预警方法 |
CN111366617A (zh) * | 2020-03-13 | 2020-07-03 | 山东航向电子科技有限公司 | 半监督模糊识别模型及基于该模型的土壤湿度测量方法 |
CN112835077A (zh) * | 2021-01-08 | 2021-05-25 | 上海海洋大学 | 基于变分模态分解的卫星导航信噪比数据去趋势方法 |
WO2022133872A1 (en) * | 2020-12-24 | 2022-06-30 | Huawei Technologies Co., Ltd. | Collaborative environment sensing in wireless networks |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846746A (zh) * | 2010-03-24 | 2010-09-29 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
CN102749067A (zh) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | 自适应多级计算的精确反演海面高程遥测装置 |
-
2016
- 2016-11-08 CN CN201610979016.0A patent/CN106291628B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846746A (zh) * | 2010-03-24 | 2010-09-29 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
CN102749067A (zh) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | 自适应多级计算的精确反演海面高程遥测装置 |
Non-Patent Citations (1)
Title |
---|
FRAN FABRA 等: "Phase Altimetry With Dual Polarization GNSS-R Over Sea Ice", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111175741A (zh) * | 2020-01-10 | 2020-05-19 | 浙江大学 | 一种单频连续毫米波多普勒传感微波墙安全空间预警方法 |
CN111175741B (zh) * | 2020-01-10 | 2022-03-18 | 浙江大学 | 一种单频连续毫米波多普勒传感微波墙安全空间预警方法 |
CN111366617A (zh) * | 2020-03-13 | 2020-07-03 | 山东航向电子科技有限公司 | 半监督模糊识别模型及基于该模型的土壤湿度测量方法 |
WO2022133872A1 (en) * | 2020-12-24 | 2022-06-30 | Huawei Technologies Co., Ltd. | Collaborative environment sensing in wireless networks |
CN112835077A (zh) * | 2021-01-08 | 2021-05-25 | 上海海洋大学 | 基于变分模态分解的卫星导航信噪比数据去趋势方法 |
CN112835077B (zh) * | 2021-01-08 | 2023-10-24 | 上海海洋大学 | 基于变分模态分解的卫星导航信噪比数据去趋势方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106291628B (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291628A (zh) | 一种提取海面卫星导航反射信号的相干观测量的方法 | |
Liu et al. | Coastal sea-level measurements based on gnss-r phase altimetry: A case study at the onsala space observatory, sweden | |
Mohamed et al. | Structural controls on groundwater flow in basement terrains: geophysical, remote sensing, and field investigations in Sinai | |
Den Ouden et al. | Stand-alone single-frequency GPS ice velocity observations on Nordenskiöldbreen, Svalbard | |
Carrano et al. | Direct measurement of the residual in the ionosphere-free linear combination during scintillation | |
Marom et al. | Energy density directional spectra of a nearshore wave field measured by interferometric synthetic aperture radar | |
Yu et al. | GNSS-based model-free sea surface height estimation in unknown sea state scenarios | |
Gruno et al. | Determining sea surface heights using small footprint airborne laser scanning | |
Suara et al. | High-resolution GNSS-tracked drifter for studying surface dispersion in shallow water | |
Forte et al. | Nearshore survey system evaluation | |
Kræmer et al. | Emulating Sentinel-1 Doppler radial ice drift measurements using Envisat ASAR data | |
CN106767383A (zh) | 基于连续运行gnss站信噪比数据的积雪深度的测量方法 | |
CN116774264B (zh) | 基于低轨卫星机会信号多普勒的运动目标定位方法 | |
Dugan et al. | Surface current measurements using airborne visible image time series | |
Liu et al. | Computing ocean surface currents from GOCI ocean color satellite imagery | |
Roussel et al. | Detection of soil moisture content changes by using a single geodetic antenna: The case of an agricultural plot | |
Walsh | Extraction of ocean wave height and dominant wavelength from GEOS 3 altimeter data | |
CN116699671A (zh) | 一种基于随机森林回归的电离层幅度闪烁指数计算方法 | |
McNinch et al. | Radar Inlet Observing System (RIOS): Continuous remote sensing of waves, currents, and bathymetry at tidal inlets | |
CN113032994B (zh) | 断面多点潮流调和常数提取计算方法、装置及电子设备 | |
CN111399012A (zh) | 一种利用gnss三频相位组合数据监测水库水位的方法 | |
Maresca Jr et al. | High-frequency skywave radar measurements of hurricane Anita | |
Marghany | Velocity bunching model for modelling wave spectra along east coast of Malaysia | |
Chien et al. | Validating of coherent-on-receive marine radar in nearshore using drifter cluster | |
Song et al. | Water Level Measurements Using Multi-station and Dual-System GNSS-MR–A Case of Shuangwangcheng Reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |