CN101846746A - 一种基于gnss-r技术的载波相位测高装置及其方法 - Google Patents
一种基于gnss-r技术的载波相位测高装置及其方法 Download PDFInfo
- Publication number
- CN101846746A CN101846746A CN 201010132135 CN201010132135A CN101846746A CN 101846746 A CN101846746 A CN 101846746A CN 201010132135 CN201010132135 CN 201010132135 CN 201010132135 A CN201010132135 A CN 201010132135A CN 101846746 A CN101846746 A CN 101846746A
- Authority
- CN
- China
- Prior art keywords
- signal
- phase
- gnss
- reflected signal
- 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.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于GNSS-R技术的载波相位测高装置及其方法,装置包括:GNSS-R接收模块、直接信号处理模块、反射信号处理模块和反射信号载波相位测高模块;所述的反射信号载波相位测高模块包括:开环跟踪单元和时间差分相位测高单元;所述的开环跟踪单元以直接信号跟踪频率为本地参考频率,有效地跟踪GNSS反射信号,得到反射信号载波相位观测量;所述的时间差分相位测高单元,利用开环跟踪单元得到的单频反射信号相位观测量,采用时间差分相位测高法反演精确水面高度。测高方法利用全球导航卫星***GNSS水面反射信号载波相位精确测量水面高度,采用开环跟踪可以得到有效的反射信号相位观测量;采用时间差分相位测高算法简单,在短时间内实现水面高精度测高。
Description
技术领域
本发明涉及GNSS-R遥感技术领域,具体涉及一种GNSS-R开环差分相位法反演水面高度的装置和方法。
背景技术
GNSS-R海洋遥感技术是利用GNSS海洋反射信号进行遥感测量的技术,其主要有无需单独发射机、信号源丰富、全球覆盖、全天候以及实时性强等优点。GNSS-R技术测量海洋高度是GNSS-R海洋遥感的重要组成部分,GNSS-R反演海面高度主要有两种途径,一种是利用直接信号和反射信号的码延迟来测高,另一种利用反射信号的载波相位测高。
码延迟测高适用范围比较广,但需要长时间平均才可以得到厘米级的测量精度。而对于载波相位测高,由于GNSS信号载波波长短(例如GPS L1波长约19cm),理论上利用载波相位测高更符合高精度的要求,但是GNSS反射信号受海面粗糙度和海面动态性的严重影响,当表面粗糙时,GNSS散射信号的相干分量减小,信号相关性差,反射信号幅度会有剧烈的高低变化,从而导致使用闭环跟踪(PLL)的接收机无法完成锁相,得不到有效的反射信号相位数据,尤其在空基和星载情况下,由于接收机的动态特性使接收机对反射信号更难跟踪,所以目前相位测高的应用主要局限在岸基试验中,目前载波相位测高主要采用的方法有下列几种:
(1)外推法跟踪反射信号测高法。当反射信号幅度小于一定门限时,使用时间外推的方法得到模型反射信号频率对反射信号进行跟踪,最后得到反射信号相位,再用零空间方法解反射信号相位模糊度后反推水面高度。但是该方法具有非常大的局限性,当水面粗糙度大时,小幅度的反射信号比例增大,至使外推的模型频率误差增大,甚至无法外推有效模型频率。
(2)多频率组合相位测高法。使用不同频率相位信号组合的方法来加长载波的波长,从而增加反射信号相干分量。欧空局提出的PIP(PARIS InterferometricProcessor)方法,以及在此基础上发展而来的PARFAIT(PARis Filtered-fieldAltImetricTracking)方法属于多频率组合相位测法,PARFAIT法还加入了可以恢复反射信号相干分量的滤波方法,利用直接和反射信号相位差分观测量建立高度解算方程,计算相位模糊度后再反推高度。其局限性在于方法操作复杂。
(3)低仰角直反信号合成相位测高法。同时接收低仰角的直接信号和反射信号,总和信号互相关波形的幅度和相位与直接和反射信号的干涉有关,而且是时间延迟的函数,利用非线性参数估计方法反演水面高度。但是该方法有一定的局限性,只能用于处理低仰角的信号。
当水面粗糙度较大,幅度小的反射信号比例增大时,使用上述方法(1)外推的模型频率误差增大,甚至无法外推有效模型频率,最终造成跟踪失败,而且方法(1)需要解算载波信号模糊度,算法复杂。上述方法(2)需要至少两个频率的反射信号观测量,反演水面高度时需要解算反射信号相位模糊度,算法复杂。上述方法(3)只能处理低仰角的GNSS-R信号,可处理的反射信号量小,而且其接收的是直接和反射信号的总信号,相互间影响会带来误差。
另外,中国专利CN101266292A公开了一种GNSS反射信号的频域处置装置和方法。该专利使用的是闭环的方法,而且没有涉及到相位测高领域。
发明内容
本发明的目的在于,提供一种基于GNSS-R技术的载波相位测高装置及其方法。
为实现上述发明目的之一,本发明的一种基于GNSS-R技术的载波相位测高装置,包括:GNSS-R接收模块、直接信号处理模块、反射信号处理模块和反射信号载波相位测高模块;其特征在于,所述的反射信号载波相位测高模块包括:开环跟踪单元和时间差分相位测高单元;
所述的开环跟踪单元以直接信号跟踪频率为本地参考频率,有效地跟踪GNSS反射信号,得到反射信号载波相位观测量;所述的时间差分相位测高单元,利用开环跟踪单元得到的单频反射信号相位观测量,采用时间差分相位测高法反演精确水面高度。
作为上述技术方案的一种改进,所述的GNSS-R接收模块包括:直接天线、反射天线和射频前端;
所述的直接天线采用右旋圆极化天线,垂直向天顶方向架设;所述的反射天线,采用左旋圆极化天线,垂直向下架设;直接天线和反射天线背靠背架设于水面上;
所述的射频前端,用于对分别经直接天线和反射天线接收的直接信号和反射信号进行变频处理,2bit量化数字采样后产生中频原始采样数据。
作为上述技术方案的又一种改进,所述的直接信号处理模块采用闭环跟踪方法处理直接信号,以50Hz输出直接信号跟踪频率和解调导航电文,为反射信号开环跟踪提供本地频率。
作为上述技术方案的再一种改进,所述的开环跟踪单元采用循环相关法进行码跟踪。
为实现上述的另一发明目的,本发明的一种基于GNSS-R技术的载波相位测高方法,该方法利用全球导航卫星***GNSS水面反射信号载波相位精确测量水面高度,该方法通过开环跟踪单元跟踪处理GNSS-R信号,以直接信号作为参考频率跟踪反射信号,得到反射信号载波相位观测量;再采用时间差分相位测高法,反演精确水面高度。
所述的方法包括以下步骤:
1)GNSS接收模块通过直接天线和反射天线得到同时接收的直接信号和反射信号的中频原始采样数据u(t),
u(t)=D(t)A(t)cos(φ(t)-φ0) (1)
其中,A(t)和φ(t)分别为信号的幅度和相位,φ0=φ(t=0),为初始时刻的相位值;
2)通过直接信号处理模块采用闭环跟踪方法处理直接信号,输出直接信号跟踪频率和导航电文D(t);
3)通过开环跟踪单元处理反射信号,以直接信号跟踪频率为本地参考频率,生成本地同相信号:υi(t)=cos(φLocal(t))和正交信号υq(t)=cos(φLocal(t));并采用循环相关法进行码跟踪,在相干积分时间T内,根据接收的GNSS-R信号的平均幅度为和平均频率为对同相信号和正交信号进行相干积分得到I路与Q路的两路信号:
其中,下标n代表时间tn≤t<tn+T,sinc(x)=sin(x)/x,Dn=D(tn);
Nn i和Nn q为接收机跟踪时I路和Q路的噪声;
最终,得到镜面反射信号的在[-π,π]内的相位:
其中,Lrn=Lr(tn)为tn时刻镜面反射信号载波相位值,其包括:起始跟踪时刻t0反射信号载波整周模糊度Nr0=Nr(t0)、反射信号接收机本地信号相位GNSS-R开环接收机输出残差相位和反射信号误差项Lr_noise的和;
其中,误差项Ld_noise包括反射接收机和GNSS卫星钟差、大气层和电离层延迟、硬件延迟和测量误差;
其中,Nr0为未知参数;
4)通过时间差分相位测高单元采用时间差分相位测高法,反演精确水面高度;
假设直接信号相位表示为:
其中,误差项Ld_noise包括:直接信号残差相位、接收机和GNSS卫星钟差、对流层和电离层延迟以及硬件延迟和测量误差;
其中,Nd0为未知参数;
则镜面反射信号相对直接信号的路径相位延迟为:
上式中,δLn为tn时刻镜面反射信号与直接信号的路径相位延迟量,Lr_d_noise为直反接收机间单差分误差,包括测量差分误差;
Hn为tn时刻反射天线相位中心到水面的距离,θn为tn时刻GNSS卫星与直接信号接收机的俯仰角,c为光速,τn为tn时刻反射信号与直接信号的路径延迟;
为了消掉未知的Nr0和Nd0,采用时间差分方法:
式中,下标j代表时刻tj,i代表时刻ti,Lj_i_r_d_noise为时间和接收机差分双差测量误差;
最后,根据公式(9),忽略小量的误差,可以得到最终的时间差分观测方程:
如果在一段观测时间内可以同时观测到两颗或更多GNSS卫星时,可以根据公式(10),利用最小二乘的方法以及平滑滤波方法计算得到观测时段内,不同时间的水面与反射信号天线相位中心的距离H。
其中,所述的时间差分相位测高方法的具体步骤包括:
1)数据筛选:在岸基静态情况下,如果某GNSS卫星反射信号残差相位随机变化严重则说明数据质量差予以剔除;根据直接和反射天线架设位置和GNSS导航卫星俯仰角和方位角,选择镜面反射点在水面的GNSS卫星信号;如果有效GNSS卫星大于等于2颗,并且各GNSS卫星反射信号同时有效跟踪时间大于或等于11分钟则截取该段时间的数据,执行以下步骤,否则结束;
2)残差相位滤波连接:由于残差相位在[-π,π]内变化,需要去除整周跳变,滤波后连接得到50Hz时11分钟多颗GNSS卫星的连续残差相位去除整周跳变,,滤波后连接得到11分钟多颗GNSS卫星的连续残差相位;
3)求首末一分钟平均水面高度:得到不同GNSS卫星连续残差相位后,根据时间差分公式,使用最小二乘法,计算间隔10分钟的首末一分钟的平均高度;
4)求中间各分钟平均水面高度:固定首一分钟的平均高度后,根据时间差分公式,使用最小二乘法正序反演得到中间时刻2~10分钟对应的高度;固定末一分钟的平均高度后,根据时间差分公式,使用最小二乘法反序反演得到中间时刻2~10分钟对应的高度;根据正反序反演结果,加权平均得到最终的中间时刻1分钟的平均高度;
5)最终得到所有时间1分钟的平均水面高度。
本发明提出了GNSS-R开环相位差分反演水面高度的装置和方法,其中开环差分相位测高法有两个关键部分组成:
一是开环跟踪反射信号技术,该方法使用直接信号跟踪频率作为参考频率,可以有效的跟踪GNSS反射信号;
另一个是时间差分相位测高法,可以利用开环跟踪得到反射信号相位观测量反演精确水面高度。
本发明提出的开环差分相位测高法,受水面粗糙度影响小,可以有效跟踪GNSS-R信号,得到准确的反射信号相位观测量,反演高度时,只需单频相位观测量,无需解算反射信号载波相位模糊度,算法简单,在较短时间内可以得到高精度水面测高结果。
本发明的优点在于:
(1)采用开环跟踪的方法可以克服反射信号幅度变化剧烈的影响,能够成功跟踪反射信号载波相位,得到有效的反射信号相位观测量;
(2)时间差分相位测高,只使用单频反射信号相位观测量,无需解算复杂的反射信号相位模糊度,算法简单,在短时间内可实现水面高精度测高。
相比较,中国专利CN100399044C公开了开环跟踪GPS信号的方法,虽然叫开环跟踪,但是与本申请处理的信号不同,实现的方式完全不同。
美国专利US2002/0130813A1公开了一种使用GNSS干涉方法测量海洋高度的方法及其装置(发明名称为:OCEAN ALTIMETRY INTERFEROMETRICMETHODANDDEVICEUSINGGNSS),这是一种使用GNSS信号测高的方法,但是该专利使用的是双频信号组合的方法,本申请与该专利相比,只使用单频信号载波相位,算法简单且效果更佳。
附图说明
图1是本发明的开环差分相位测高装置的总体框图;
图2是常规的GNSS-R数据采集***的示意图;
图3是传统闭环(PLL)跟踪方法的原理示意图;
图4是本发明的GNSS-R信号开环跟踪方法的原理示意图;
图5是本发明的开环差分相位测高装置中直接和反射信号路径差示意图;
图6是本发明时间差分相位测高方法的流程图;
图7是三营桥2009年9月22日11:45-12:07GPS方位角和俯仰角示意图;
图8是本发明一实施例的反射信号开环跟踪残差相位图;
图9是本发明的差分相位测高结果与GPS差分定位结果对比图。
具体实施方式
参考下面提出的各种细节来描述本发明,并且附图给出本发明的例示。下面的描述和附图对于本发明而言仅是示例,不能理解为是对本发明的限制。
本发明GNSS-R数据采集***是软件接收机***,直接和反射信号分别经天线接收后,在射频前端下变频,通过USB接口输入到数据接收计算机的硬盘中,最后得到同时接收的直接信号和反射信号的中频采样原始数据。
原始采样数据通过GNSS-R数据采集***的预处理单元(包括直接信号和反射信号处理单元,开环跟踪是反射信号处理的重要部分)预处理模块将数据从原始中频采样数据处理为0级GNSS-R数据,主要包括:直接信号和反射信号的相关值-码相位数据(可用于码相位测高)和反射信号残差相位观测量(用于开环差分相位测高)等,其主要功能是直接信号捕获跟踪和导航电文解调,反射信号电文解调、反射信号开环跟踪和反射信号路径相位延迟量生成等。
0级数据通过GNSS-R数据采集***数据后处理模块(包括时间差分相位测高)处理得到GNSS-R测高1级数据,主要是水面与反射天线相位中心之间的高度,后处理模块主要算法为开环差分相位测高法。
具体步骤如下:
步骤1:按附图1连接好各部件,其中GNSS直接信号通过直接天线(右旋圆极化天线)接收,GNSS-R信号通过反射天线(左旋圆极化天线)接收,直接天线垂直向天顶方向架设,反射天线垂直向下架设,两天线背靠背架设于水面上。
步骤2:设置数据采集***(组成如附图2),直接和反射信号分别经天线接收后在射频前端下变频,2bit量化数字采样后产生原始采样数据,采样率为16.368MHz,最后通过USB接口输入到计算机的硬盘中。最后我们得到的是同时接收的直接信号和反射信号的中频采样原始数据约1小时以上。
步骤3:处理直接信号。使用闭环跟踪方法处理直接信号,以50Hz输出直接信号跟踪频率,解调导航电文并输出,为反射信号开环跟踪做准备。
步骤4:开环跟踪,处理反射信号。以直接信号跟踪频率为本地参考频率,进行开环频率跟踪,使用循环相关法进行码跟踪,同相和正交分量20ms相干积分滤波后,使用直接信号跟踪输出的导航电文辅助解调反射信号导航电文,最后使用4象鉴相器输出在[-π,π]内的残差相位,输出速率为50Hz。
步骤5:时间差分相位测高。
1)数据筛选:判断是否同时有两颗或两颗以上,有效跟踪时间大于等于11分钟的数据,有则截取该段时间的数据,执行以下步骤,否则结束;
2)残差相位滤波连接;
3)利用时间差分公式,使用最小二乘法求起始和结束一分钟的平均高度;
4)正序和逆序法分别反演中间时刻的高度,然后加权平均得到中间时刻各分钟的平均高度;
5)得到每分钟平均水面高度。
下面详细描述开环跟踪反射信号以及时间差分算法计算水面高度的原理:
1、开环跟踪反射信号:
闭环(PLL)跟踪方法见附图3,处理时需要构成反馈回路,该方法对于直接GNSS信号接收时最优选择。但是GNSS反射信号为弱信号,而且其振幅变化大,频率被展宽,从而导致闭环跟踪失锁。
本发明采用开环跟踪方法,见附图4,该方法无需构成回路,本地频率参考直接信号跟踪频率,具体原理如下:
假设GNSS-R接收机接收信号u(t):
u(t)=D(t)A(t)cos(φ(t)-φ0)(1)
其中,接收信号幅度和相位为A(t)和φ(t),导航电文为D(t),φ0=φ(t=0),即初始时刻的相位值,公式(1)中忽略了C/A码的BPSK调制,开环跟踪时,使用循环相关法将其解调。
开环跟踪以直接信号跟踪频率为模型参考频率,生成本地同相信号υi(t)=cos(φLocal(t))和正交信号υq(t)=cos(φLocal(t))分别与接收信号u(t)进行相关处理。假设在相干积分时间T内,接收的GNSS-R信号的平均幅度为平均频率为时,同相和正交信号在T内的相干积分结果为:
其中,下标n代表时间tn≤t<tn+T,sinc(x)=sin(x)/x,Dn=D(tn)。是tn≤t<tn+T时间内接收信号平均频率和本地生成平均频率的差值。是在tn时刻接收信号相位和本地相位的初始相位差,其中本地相位φn Local和接收信号相位φn用时间累积和表示: Nn i和Nn q为接收机跟踪时I路和Q路的噪声。
arctan2代表四象鉴相器,ωn是tn≤t<tn+T时间内接收信号角频率和本地生成角频率的平均差值,Nn为tn时相对t0时刻残差相位变化整周数。
最终可得镜面反射信号的相位:
其中,Lrn=Lr(tn)为tn时刻镜面反射信号载波相位值,其包括起始跟踪时刻t0反射信号载波整周模糊度Nr0=Nr(t0)、反射信号接收机本地信号相位(以直接信号相位作为参考)GNSS-R开环接收机输出残差相位和反射信号误差项Lr_noise的和。其中误差项Ld_noise主要包括反射接收机和GNSS卫星钟差、大气层和电离层延迟以及硬件延迟和测量误差等。其中Nr0为未知参数。
2、时间差分算法计算水面高度:
直接信号相位可以表示为:
其中,Ldn=Ld(tn)为tn时刻直接信号载波相位值,其主要包括起始跟踪时刻t0直接信号载波整周模糊度Nd0=Nd(t0)、tn时刻直接信号接收机输出相位和直接信号误差项Ld_noise。其中误差项Ld_noise主要包括直接信号残差相位(直接信号使用闭环跟踪,残差相位趋于0)、接收机和GNSS卫星钟差、对流层和电离层延迟以及硬件延迟和测量误差等。其中Nd0为未知参数。
镜面反射信号相位见公式(6),则镜面反射信号相对直接信号的路径相位延迟为:
见附图5,上式中,δLn为tn时刻镜面反射信号与直接信号的路径相位延迟量,Lr_d_noise为直反接收机间单差分误差,主要包括测量差分误差。因为直接接收机和反射接收机位置非常接近,而且在岸基情况下,直接信号和反射信号路径相似,所以GNSS卫星钟差、大气层和电离层延迟可以差分掉,此外,直接信号接收机和GNSS-R接收机使用同一个晶振,所以接收机钟差可以去掉。Hn为tn时刻反射天线相位中心到水面的距离,θn为tn时刻GNSS卫星与直接信号接收机的俯仰角,c为光速,τn为tn时刻反射信号与直接信号的路径延迟。
为了消掉未知的Nr0和Nd0,采用时间差分方法:
式中,下标j代表时刻tj,i代表时刻ti,Lj_i_r_d_noise为双差(时间以及接收机差分)测量误差。其他含义见上面的描述。由上式可知,GNSS-R开环跟踪后输出的残差相位时间差代表的物理含义是反射信号相对直接信号在j时刻的路径相位延迟量和i时刻路径相位延迟量的差,即反射信号路径相位延迟量的时间差。
综上所述,根据公式(9),忽略小量的误差,可以得到最终的时间差分观测方程:
式中,上标S代表某GNSS卫星S,下标j代表时刻tj,i代表时刻ti,为GNSS-R开环接收机输出残差相位,H为某时刻反射天线相位中心到水面的距离,θ为某时刻GNSS卫星S与直接信号接收机的俯仰角。
如果在一段观测时间内可以同时观测到两颗或更多GNSS卫星时,可以根据公式(10),利用最小二乘的方法以及平滑滤波方法计算得到观测时段内,不同时间的水面与反射信号天线相位中心的距离H。
3、时间差分相位测高,具体流程图见附图6,主要包括:
1)数据筛选。在岸基静态情况下,滤波输出的反射信号残差相位应该呈有规律渐变的过程,如果某GNSS卫星反射信号残差相位随机变化严重则说明数据质量差予以剔除。根据直接和反射天线架设位置和GNSS导航卫星俯仰角和方位角,选择镜面反射点在水面的GNSS卫星信号。如果有效GNSS卫星大于等于2颗,并且各GNSS卫星反射信号同时有效跟踪时间大于等于11分钟则截取该段时间的数据,执行以下步骤,否则结束。
2)残差相位滤波连接。由于步骤3中输出的残差相位在[-π,π]内变化,需要去除整周跳变,滤波后连接得到11分钟多颗GNSS卫星的连续残差相位(50Hz)。
3)求首末一分钟平均水面高度。得到不同GNSS卫星连续残差相位后,根据时间差分公式(其中对应两个时间的高度均为未知数),使用最小二乘法,计算间隔10分钟的首末一分钟的平均高度。
4)求中间各分钟平均水面高度。固定首分钟(第1分钟)的平均高度后,根据时间差分公式(中间时刻对应的高度为未知数),使用最小二乘法得到中间时刻(2~10分钟)对应的高度,称为正序反演。固定末分钟(第11分钟)的平均高度后,根据时间差分公式(中间时刻对应的高度为未知数),使用最小二乘法得到中间时刻(2~10分钟)对应的高度(反序反演)。根据正反序反演结果,加权平均得到最终的中间时刻1分钟的平均高度。
5)最终得到所有时间1分钟的平均水面高度。
式中,上标S代表某GNSS卫星S,下标j代表时刻tj,i代表时刻ti,为GNSS-R开环接收机输出残差相位,H为某时刻反射天线相位中心到水面的距离,θ为某时刻GNSS卫星S与直接信号接收机的俯仰角。(参见发明内容部分)
实验验证
实验概况
中国科学院空间科学与应用研究中心(空间中心)和北京应用气象研究所在河北官厅水库铁路桥上联合开展了两次GNSS-R岸基测高实验。第一次观测时间为2009年9月22日11时至14时,共3个小时,天气晴朗微风1到2级,第二次观测时间2009年10月21日11时至17时,共6个小时,天气晴朗,2到4级风。实验目的是是验证开环差分相位测高法的可行性及其测高精度。
见附图7,将GNSS-R数据采集***架设在三营桥中段的西南侧,其主要包括直接信号天线、反射信号天线和信号中频采集***。此外,桥上、岸基及水面浮标均设置GPS双频定位***,使用三点差分定位来获得精确的水面距桥上定位天线的精确高度。
实验结果
采用2009年9月22日,11:45-12:07,将22分钟数据进行开环差分相位测高反演,该时段的GPS卫星星空图见附图7,根据GNSS卫星、天线架设位置和信号质量进行选星和质量控制,选择PRN7、11和25三颗卫星数据进行处理。
在预处理模块中采用以直接信号跟踪频率为参考频率的开环跟踪,得到开环跟踪输出残差相位,见附图8,图中显示了在11分钟内PRN7、11和25号星反射信号开环跟踪输出残差相位变化图,其中PRN13俯仰角变化很小,其随时间输出残差相位也变化较小,PRN19为上升星,其随时间输出残差相位递减,PRN23为下降星,其随时间输出残差相位递增。
得到反射信号残差相位观测量后,使用时间差分法进行水面高度反演。根据相位时间差分公式(10),利用最小二乘法和平滑算法计算间隔最远的首末一分钟的平均水面高度,然后固定首末一分钟的高度,采用顺序和逆序反演中间时刻的数据,最后将顺序和逆序反演结果加权平均,数据滤波后得到最后每分钟的水面与反射天线相位中心高度。
与GPS差分精密定位得到的水面高度结果相比见附图9,图中,连线为GPS差分测高每秒的测高结果,黑叉为GPS差分测高每分钟的测高结果,点为22分钟开环差分相位测高的结果,星号为前11分钟的开环差分相位测高结果,加号为后11分钟的开环差分相位测高结果。其中22分钟开环差分相位测高结果与每分钟的GPS差分测高结果相比,其误差平均值为0.1067cm,误差标准差为3.3768cm;前11分钟开环差分相位测高结果与每分钟的GPS差分测高结果相比,其误差平均值为7.8179cm,误差标准差为3.5472cm;后11分钟开环差分相位测高结果与每分钟的GPS差分测高结果相比,其误差平均值为-5.9534cm,误差标准差为0.3404cm。
使用开环时间差分测高法处理不同时间段的数据结果如表1:
表1算例结果表
时间 | 使用的GPS卫星 | 误差平均值 | 误差标准差 |
2009.09.22,11:45-12:07 | PRN13、19、23 | 0.1067cm | 3.3768cm |
2009.09.22,11:45-11:56 | PRN13、19、23 | 7.8179cm | 3.5472cm |
2009.09.22,11:56-12:07 | PRN13、19、23 | -5.9534cm | 0.3404cm |
2009.09.22,12:22-12:33 | PRN13、23、25 | -7.8989cm | 0.5554cm |
2009.10.21,11:15-11:26 | PRN11、13、25 | 7.8989cm | 0.4744cm |
2009.10.21,11:54-12:05 | PRN07、11、25 | 5.1101cm | 0.2699cm |
通过以上实验可知,该方法有两个关键技术组成,一是开环跟踪的方法跟踪处理GNSS-R信号,该方法以直接信号作为参考频率跟踪反射信号,得到反射信号载波相位观测量;另一个关键技术是时间差分相位测高法,如果在一段时间内同时跟踪到2颗或更多GNSS-R信号的相位观测量,可以利用该方法反演精确水面高度。本例中,利用11分钟三颗GPS反射信号的相位观测量数据,反演水面高度结果与双频GPS差分定位结果一致。
结论:使用GNSS-R开环差分相位测高法处理三营桥GNSS-R实验数据,结果与GPS差分定位测高数据相比,可以得到1分钟平均测高结果精度为厘米级,验证了GNSS-R开环差分相位法高精度测量水面高度的可行性。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种基于GNSS-R技术的载波相位测高装置,包括:GNSS-R接收模块、直接信号处理模块、反射信号处理模块和反射信号载波相位测高模块;其特征在于,所述的反射信号载波相位测高模块包括:开环跟踪单元和时间差分相位测高单元;
所述的开环跟踪单元以直接信号跟踪频率为本地参考频率,有效地跟踪GNSS反射信号,得到反射信号载波相位观测量;所述的时间差分相位测高单元,利用开环跟踪单元得到的单频反射信号相位观测量,采用时间差分相位测高法反演精确水面高度。
2.根据权利要求1所述的基于GNSS-R技术的载波相位测高装置,其特征在于,所述的GNSS-R接收模块包括:直接天线、反射天线和射频前端;
所述的直接天线采用右旋圆极化天线,垂直向天顶方向架设;所述的反射天线,采用左旋圆极化天线,垂直向下架设;直接天线和反射天线背靠背架设于水面上;
所述的射频前端,用于对分别经直接天线和反射天线接收的直接信号和反射信号进行变频处理,2bit量化数字采样后产生中频原始采样数据。
3.根据权利要求1所述的基于GNSS-R技术的载波相位测高装置,其特征在于,所述的直接信号处理模块采用闭环跟踪方法处理直接信号,以50Hz输出直接信号跟踪频率和解调导航电文,为反射信号开环跟踪提供本地频率。
4.根据权利要求1所述的基于GNSS-R技术的载波相位测高装置,其特征在于,所述的开环跟踪单元采用循环相关法进行码跟踪。
5.一种基于GNSS-R技术的载波相位测高方法,该方法利用全球导航卫星***GNSS水面反射信号载波相位精确测量水面高度,该方法通过开环跟踪单元跟踪处理GNSS-R信号,以直接信号作为参考频率跟踪反射信号,得到反射信号载波相位观测量;再采用时间差分相位测高法,反演精确水面高度。
6.根据权利要求5所述的基于GNSS-R技术的载波相位测高方法,其特征在于,所述的方法包括以下步骤:
1)GNSS接收模块通过直接天线和反射天线得到同时接收的直接信号和反射信号的中频原始采样数据u(t),
u(t)=D(t)A(t)cos(φ(t)-φ0) (1)
其中,A(t)和φ(t)分别为信号的幅度和相位,φ0=φ(t=0),为初始时刻的相位值;
2)通过直接信号处理模块采用闭环跟踪方法处理直接信号,输出直接信号跟踪频率和导航电文D(t);
3)通过开环跟踪单元处理反射信号,以直接信号跟踪频率为本地参考频率,生成本地同相信号:υi(t)=cos(φLocal(t))和正交信号υq(t)=cos(φLocal(t));并采用循环相关法进行码跟踪,在相干积分时间T内,根据接收的GNSS-R信号的平均幅度为和平均频率为对同相信号和正交信号进行相干积分得到I路与Q路的两路信号:
其中,下标n代表时间tn≤t<tn+T,sinc(x)=sin(x)/x,Dn=D(tn);
Nn i和Nn q为接收机跟踪时I路和Q路的噪声;
其中,arctan2代表四象鉴相器,ωn是tn≤t<tn+T时间内接收信号角频率和本地生成角频率的平均差值,Nn为tn时刻相对t0时刻的残差相位变化整周数;
最终,得到镜面反射信号的在[-π,π]内的相位:
其中,Lrn=Lr(tn)为tn时刻镜面反射信号载波相位值,其包括:起始跟踪时刻t0反射信号载波整周模糊度Nr0=Nr(t0)、反射信号接收机本地信号相位、GNSS-R开环接收机输出残差相位和反射信号误差项Lr_noise的和;
其中,误差项Ld_noise包括反射接收机和GNSS卫星钟差、大气层和电离层延迟、硬件延迟和测量误差;
其中,Nr0为未知参数;
4)通过时间差分相位测高单元采用时间差分相位测高法,反演精确水面高度;
假设直接信号相位表示为:
其中,误差项Ld_noise包括:直接信号残差相位、接收机和GNSS卫星钟差、对流层和电离层延迟以及硬件延迟和测量误差;
其中,Nd0为未知参数;
则镜面反射信号相对直接信号的路径相位延迟为:
上式中,δLn为tn时刻镜面反射信号与直接信号的路径相位延迟量,Lr_d_noise为直反接收机间单差分误差,包括测量差分误差;
Hn为tn时刻反射天线相位中心到水面的距离,θn为tn时刻GNSS卫星与直接信号接收机的俯仰角,c为光速,τn为tn时刻反射信号与直接信号的路径延迟;
为了消掉未知的Nr0和Nd0,采用时间差分方法:
式中,下标j代表时刻tj,i代表时刻ti,Lj_i_r_d_noise为时间和接收机差分双差测量误差;
最后,根据公式(9),忽略小量的误差,可以得到最终的时间差分观测方程:
如果在一段观测时间内可以同时观测到两颗或更多GNSS卫星时,根据公式(10),利用最小二乘的方法以及平滑滤波方法计算得到观测时段内,不同时间的水面与反射信号天线相位中心的距离H。
7.权利要求5或6所述的基于GNSS-R技术的载波相位测高方法,其特征在于,所述的时间差分相位测高方法的具体步骤包括:
1)数据筛选:在岸基静态情况下,如果某GNSS卫星反射信号残差相位随机变化严重则说明数据质量差予以剔除;根据直接和反射天线架设位置和GNSS导航卫星俯仰角和方位角,选择镜面反射点在水面的GNSS卫星信号;如果有效GNSS卫星大于等于2颗,并且各GNSS卫星反射信号同时有效跟踪时间大于或等于11分钟则截取该段时间的数据,执行以下步骤,否则结束;
2)残差相位滤波连接:由于残差相位在[-π,π]内变化,需要去除整周跳变,滤波后连接得到50Hz时11分钟多颗GNSS卫星的连续残差相位去除整周跳变,,滤波后连接得到11分钟多颗GNSS卫星的连续残差相位;
3)求首末一分钟平均水面高度:得到不同GNSS卫星连续残差相位后,根据时间差分公式,使用最小二乘法,计算间隔10分钟的首末一分钟的平均高度;
4)求中间各分钟平均水面高度:固定首一分钟的平均高度后,根据时间差分公式,使用最小二乘法正序反演得到中间时刻2~10分钟对应的高度;固定末一分钟的平均高度后,根据时间差分公式,使用最小二乘法反序反演得到中间时刻2~10分钟对应的高度;根据正反序反演结果,加权平均得到最终的中间时刻1分钟的平均高度;
5)最终得到所有时间1分钟的平均水面高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010132135 CN101846746B (zh) | 2010-03-24 | 2010-03-24 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010132135 CN101846746B (zh) | 2010-03-24 | 2010-03-24 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101846746A true CN101846746A (zh) | 2010-09-29 |
CN101846746B CN101846746B (zh) | 2013-07-24 |
Family
ID=42771425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010132135 Expired - Fee Related CN101846746B (zh) | 2010-03-24 | 2010-03-24 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101846746B (zh) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102055512A (zh) * | 2010-09-30 | 2011-05-11 | 中国科学院国家天文台 | 天线组阵中基于循环自相关的信号相位差估计装置与方法 |
CN103278083A (zh) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Gnss-r海冰厚度探测设备及利用该设备探测海冰厚度的方法 |
CN104678405A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于多天线的岸基gnss-r智能采集装置 |
CN104678354A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于gnss-r的海面目标双星定位方法 |
CN104678388A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种用于提高测高精度的gnss-r正交化交错处理装置 |
CN104808223A (zh) * | 2015-04-30 | 2015-07-29 | 北京航空航天大学 | 一种适用于星载gnss-r接收机的具有可抑制镜点高动态的相关器 |
CN104865581A (zh) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与*** |
CN105241456A (zh) * | 2015-11-06 | 2016-01-13 | 金陵科技学院 | 巡飞弹高精度组合导航方法 |
CN105334510A (zh) * | 2015-10-19 | 2016-02-17 | 叶雷 | 一种gnss-r地表探测装置和方法 |
CN106291628A (zh) * | 2016-11-08 | 2017-01-04 | 上海海事大学 | 一种提取海面卫星导航反射信号的相干观测量的方法 |
CN107092027A (zh) * | 2017-03-20 | 2017-08-25 | 深圳市西博泰科电子有限公司 | 利用信号接收器计算姿态角的方法及设备 |
CN107607032A (zh) * | 2017-10-23 | 2018-01-19 | 上海普适导航科技股份有限公司 | 一种gnss形变监测*** |
CN107796484A (zh) * | 2017-01-11 | 2018-03-13 | 中南大学 | 一种基于北斗导航定位***信噪比数据观测水位变化方法 |
CN108761194A (zh) * | 2018-03-16 | 2018-11-06 | 贵州电网有限责任公司 | 基于移频滤波算法的微功率无线载波参数测量方法 |
CN110114696A (zh) * | 2016-12-15 | 2019-08-09 | 深圳开阳电子股份有限公司 | 用于卫星导航***的信号捕获方法及接收机 |
CN110244291A (zh) * | 2019-07-19 | 2019-09-17 | 北京航天飞行控制中心 | 基于无线电信号处理的测速方法及装置 |
CN110501728A (zh) * | 2018-05-16 | 2019-11-26 | 清华大学 | 定位基站跳时信号的鉴频方法及鉴频装置 |
CN112147643A (zh) * | 2020-08-24 | 2020-12-29 | 西安空间无线电技术研究所 | 一种gnss-r***成像预处理的方法及装置 |
CN112731475A (zh) * | 2020-12-25 | 2021-04-30 | 中国科学院国家空间科学中心 | 一种gnss掩星双支路信号开环跟踪方法 |
CN112880633A (zh) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | 基于伯格算法的海面高度测量方法 |
CN113296135A (zh) * | 2020-02-24 | 2021-08-24 | 香港理工大学深圳研究院 | 形变监测方法、装置和接收机 |
CN113691477A (zh) * | 2020-05-19 | 2021-11-23 | 大唐移动通信设备有限公司 | 一种载波相位跟踪方法及设备 |
CN113723129A (zh) * | 2020-12-03 | 2021-11-30 | 深圳市汇顶科技股份有限公司 | 具有不稳定的基准时钟的调幅信号的射频载波跟踪 |
CN115792986A (zh) * | 2023-02-03 | 2023-03-14 | 山东大学 | 一种gnss-r海面测高实时软件接收机设计方法 |
CN116719025A (zh) * | 2023-08-09 | 2023-09-08 | 中国科学院地理科学与资源研究所 | 一种内陆水面高程测量方法 |
CN117420581A (zh) * | 2023-12-19 | 2024-01-19 | 自然资源部第一海洋研究所 | 单频卫星浮标时间差分载波相位波浪参数反演方法及*** |
CN117607908A (zh) * | 2024-01-23 | 2024-02-27 | 山东大学 | 无人机载平台的gnss-r码相位海面测高方法及装置 |
-
2010
- 2010-03-24 CN CN 201010132135 patent/CN101846746B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
《海洋测绘》 20100331 邵连军等 GNSS_R海面测高算法 全文 1-7 , * |
《科学技术与工程》 20090531 陶鹏等 用于海洋测高的GPS_R软件接收机的设计与验证 全文 1-7 , * |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102055512B (zh) * | 2010-09-30 | 2013-08-07 | 中国科学院国家天文台 | 天线组阵中基于循环自相关的信号相位差估计装置与方法 |
CN102055512A (zh) * | 2010-09-30 | 2011-05-11 | 中国科学院国家天文台 | 天线组阵中基于循环自相关的信号相位差估计装置与方法 |
CN103278083A (zh) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Gnss-r海冰厚度探测设备及利用该设备探测海冰厚度的方法 |
CN103278083B (zh) * | 2013-05-08 | 2015-12-09 | 南京信息工程大学 | 利用gnss-r海冰厚度探测设备探测海冰厚度的方法 |
CN104678354B (zh) * | 2015-02-11 | 2017-04-05 | 北京航空航天大学 | 一种基于gnss‑r的海面目标双星定位方法 |
CN104678405A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于多天线的岸基gnss-r智能采集装置 |
CN104678354A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于gnss-r的海面目标双星定位方法 |
CN104678388A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种用于提高测高精度的gnss-r正交化交错处理装置 |
CN104865581A (zh) * | 2015-03-19 | 2015-08-26 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与*** |
WO2016145723A1 (zh) * | 2015-03-19 | 2016-09-22 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与*** |
CN104808223A (zh) * | 2015-04-30 | 2015-07-29 | 北京航空航天大学 | 一种适用于星载gnss-r接收机的具有可抑制镜点高动态的相关器 |
CN104808223B (zh) * | 2015-04-30 | 2017-07-04 | 北京航空航天大学 | 一种适用于星载gnss‑r接收机的具有可抑制镜点高动态的相关器 |
CN105334510A (zh) * | 2015-10-19 | 2016-02-17 | 叶雷 | 一种gnss-r地表探测装置和方法 |
CN105334510B (zh) * | 2015-10-19 | 2018-05-08 | 泛太通信导航(深圳)有限公司 | 一种gnss-r地表探测装置和方法 |
CN105241456A (zh) * | 2015-11-06 | 2016-01-13 | 金陵科技学院 | 巡飞弹高精度组合导航方法 |
CN105241456B (zh) * | 2015-11-06 | 2018-02-02 | 金陵科技学院 | 巡飞弹高精度组合导航方法 |
CN106291628A (zh) * | 2016-11-08 | 2017-01-04 | 上海海事大学 | 一种提取海面卫星导航反射信号的相干观测量的方法 |
CN110114696A (zh) * | 2016-12-15 | 2019-08-09 | 深圳开阳电子股份有限公司 | 用于卫星导航***的信号捕获方法及接收机 |
CN110114696B (zh) * | 2016-12-15 | 2023-01-13 | 深圳开阳电子股份有限公司 | 用于卫星导航***的信号捕获方法及接收机 |
CN107796484A (zh) * | 2017-01-11 | 2018-03-13 | 中南大学 | 一种基于北斗导航定位***信噪比数据观测水位变化方法 |
CN107092027B (zh) * | 2017-03-20 | 2019-08-13 | 深圳市西博泰科电子有限公司 | 利用信号接收器计算姿态角的方法及设备 |
CN107092027A (zh) * | 2017-03-20 | 2017-08-25 | 深圳市西博泰科电子有限公司 | 利用信号接收器计算姿态角的方法及设备 |
CN107607032A (zh) * | 2017-10-23 | 2018-01-19 | 上海普适导航科技股份有限公司 | 一种gnss形变监测*** |
CN107607032B (zh) * | 2017-10-23 | 2020-02-11 | 上海普适导航科技股份有限公司 | 一种gnss形变监测*** |
CN108761194B (zh) * | 2018-03-16 | 2020-06-16 | 贵州电网有限责任公司 | 基于移频滤波算法的微功率无线载波参数测量方法 |
CN108761194A (zh) * | 2018-03-16 | 2018-11-06 | 贵州电网有限责任公司 | 基于移频滤波算法的微功率无线载波参数测量方法 |
CN110501728B (zh) * | 2018-05-16 | 2022-03-29 | 清华大学 | 定位基站跳时信号的鉴频方法及鉴频装置 |
CN110501728A (zh) * | 2018-05-16 | 2019-11-26 | 清华大学 | 定位基站跳时信号的鉴频方法及鉴频装置 |
CN110244291A (zh) * | 2019-07-19 | 2019-09-17 | 北京航天飞行控制中心 | 基于无线电信号处理的测速方法及装置 |
CN113296135B (zh) * | 2020-02-24 | 2024-03-29 | 香港理工大学深圳研究院 | 形变监测方法、装置和接收机 |
CN113296135A (zh) * | 2020-02-24 | 2021-08-24 | 香港理工大学深圳研究院 | 形变监测方法、装置和接收机 |
CN113691477A (zh) * | 2020-05-19 | 2021-11-23 | 大唐移动通信设备有限公司 | 一种载波相位跟踪方法及设备 |
CN113691477B (zh) * | 2020-05-19 | 2022-10-18 | 大唐移动通信设备有限公司 | 一种载波相位跟踪方法及设备 |
CN112147643A (zh) * | 2020-08-24 | 2020-12-29 | 西安空间无线电技术研究所 | 一种gnss-r***成像预处理的方法及装置 |
CN112147643B (zh) * | 2020-08-24 | 2024-02-09 | 西安空间无线电技术研究所 | 一种gnss-r***成像预处理的方法及装置 |
CN113723129A (zh) * | 2020-12-03 | 2021-11-30 | 深圳市汇顶科技股份有限公司 | 具有不稳定的基准时钟的调幅信号的射频载波跟踪 |
CN112731475B (zh) * | 2020-12-25 | 2023-08-08 | 中国科学院国家空间科学中心 | 一种gnss掩星双支路信号开环跟踪方法 |
CN112731475A (zh) * | 2020-12-25 | 2021-04-30 | 中国科学院国家空间科学中心 | 一种gnss掩星双支路信号开环跟踪方法 |
CN112880633A (zh) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | 基于伯格算法的海面高度测量方法 |
CN115792986A (zh) * | 2023-02-03 | 2023-03-14 | 山东大学 | 一种gnss-r海面测高实时软件接收机设计方法 |
CN116719025A (zh) * | 2023-08-09 | 2023-09-08 | 中国科学院地理科学与资源研究所 | 一种内陆水面高程测量方法 |
CN116719025B (zh) * | 2023-08-09 | 2023-10-13 | 中国科学院地理科学与资源研究所 | 一种内陆水面高程测量方法 |
CN117420581A (zh) * | 2023-12-19 | 2024-01-19 | 自然资源部第一海洋研究所 | 单频卫星浮标时间差分载波相位波浪参数反演方法及*** |
CN117420581B (zh) * | 2023-12-19 | 2024-03-08 | 自然资源部第一海洋研究所 | 单频卫星浮标时间差分载波相位波浪参数反演方法及*** |
CN117607908A (zh) * | 2024-01-23 | 2024-02-27 | 山东大学 | 无人机载平台的gnss-r码相位海面测高方法及装置 |
CN117607908B (zh) * | 2024-01-23 | 2024-04-05 | 山东大学 | 无人机载平台的gnss-r码相位海面测高方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101846746B (zh) | 2013-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101846746B (zh) | 一种基于gnss-r技术的载波相位测高装置及其方法 | |
AU2020103449A4 (en) | Method for monitoring the water level of reservoir by using GNSS triple-frequency phase combination data | |
CN104865581B (zh) | 一种基于gnss‑r技术的河流流速测量方法与*** | |
Bilich et al. | Correction published 29 March 2008: Mapping the GPS multipath environment using the signal-to-noise ratio (SNR) | |
CN103235327B (zh) | 一种gnss/mins超深组合导航方法、***及装置 | |
CN101246210B (zh) | 一种gps定位方法 | |
CN101403790A (zh) | 单频gps接收机的精密单点定位方法 | |
CN104597471A (zh) | 面向时钟同步多天线gnss接收机的定向测姿方法 | |
CN101833080A (zh) | 一种利用gps***附加约束条件的载体姿态测量方法 | |
CN106443726A (zh) | 一种基于预滤波的gnss矢量跟踪环路及其实现方法 | |
CN106842263A (zh) | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 | |
CN102033236A (zh) | 一种卫星导航位置速度联合估计方法 | |
CN106767383A (zh) | 基于连续运行gnss站信噪比数据的积雪深度的测量方法 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN102193095A (zh) | 全球导航卫星***接收器的定位方法 | |
CN103901442A (zh) | 一种单频多路径效应探测与缓解方法 | |
CN106768179A (zh) | 基于连续运行gnss站信噪比数据的潮位的测量方法 | |
CN116699658A (zh) | 一种基于双频反射信号组合的gnss-r载波相位海面测高的方法 | |
Renkwitz et al. | On the role of anisotropic MF/HF scattering in mesospheric wind estimation | |
CN108279425A (zh) | 一种多频观测过程中多径误差的修正方法 | |
CN104252004B (zh) | 利用单天线导航接收机测量自旋卫星姿态的***及方法 | |
Chen et al. | Feasibility and performance evaluation of low-cost GNSS devices for sea level measurement based on GNSS-IR | |
Lau et al. | Statistical characterization of the meteor trail distribution at the South Pole as seen by a VHF interferometric meteor radar | |
CN103698793B (zh) | 基于软件接收机的gnss信号模拟器仿真角速度范围测量方法 | |
Bai et al. | GNSS-R open-loop difference phase altimetry: Results from a bridge experiment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100084 No. 1, No. 2, South of Zhongguancun, Haidian District, Beijing Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS Address before: 100084 No. 1, No. 2, South of Zhongguancun, Haidian District, Beijing Patentee before: Space Science & Applied Research Centre, Chinese Academy of Sciences |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130724 Termination date: 20200324 |