CN104678371B - 一种基于时延修正的海面高度测量装置 - Google Patents
一种基于时延修正的海面高度测量装置 Download PDFInfo
- Publication number
- CN104678371B CN104678371B CN201510071278.2A CN201510071278A CN104678371B CN 104678371 B CN104678371 B CN 104678371B CN 201510071278 A CN201510071278 A CN 201510071278A CN 104678371 B CN104678371 B CN 104678371B
- Authority
- CN
- China
- Prior art keywords
- signal
- mrow
- module
- satellite
- signal processing
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种基于时延修正的海面高度测量装置,由导航接收机、右旋天线、左旋天线、射频前端模块、中频信号处理模块和反演模块组成;导航接收机与右旋天线、中频信号处理模块和反演模块相连,右旋天线和左旋天线彼此独立并一起与射频前端模块、中频信号处理模块和反演模块依次相连;右旋天线输出的直射射频电压信号和左旋天线输出的反射射频电压信号作为射频前端模块的输入,射频前端模块输出的直射数字中频信号和反射数字中频信号作为中频信号处理模块的输入,中频信号处理模块输出的时延多普勒二维相关功率作为反演模块的输入,右旋天线输出的直射射频电压信号作为导航接收机的输入,导航接收机输出的导航信息作为中频信号处理模块和反演模块的输入。
Description
技术领域
本发明涉及一种基于时延修正的海面高度测量装置。具体地说,是利用GNSS-R(Global Navigation Satellites System-Reflectmetry)海面风场反演结果中由于海面粗糙度引起的时延误差,来修正测高过程中由于海面粗糙度引起的时延误差,从而提高测高精度。
背景技术
GNSS-R作为一种新的遥感手段,受到国内外学者的广泛关注。目前,GNSS-R在理论、技术和数据反演等方面逐步完善。基于GNSS-R的海洋遥感是GNSS-R遥感领域的一个主要部分,其主要发展方向是海面风场和海面高度。迄今为止,国际上相继开展多次探测海面风场和海面高度的实验,GNSS-R反演海面风场的理论和算法已经基本成形。
在反演海面风场中,理论相关功率模型较为成熟,将理论相关功率曲线的后沿与实测相关功率曲线的后沿采用最小二乘法进行匹配,可以估计出海面风速。以反演出的风速作为辅助,通过计算残差平方和可以估计出海面风向。在反演海面高度中,测高精度的关键在于精确计算直射信号和经镜面反射的反射信号相关功率峰值的时间差。但是,由于海面粗糙度的影响,在镜面反射点周围存在漫反射点,造成了反射信号相关功率的峰值点位置向后偏移。这时若使用传统的跟踪峰值点或波形的半功率点方法测量信号延迟,则会带来较大的偏差。
在海面高度测量中,针对解决海面粗糙度引起的测量误差缺乏一种简单直接的解决办法。
发明内容
本发明的目的是提供一种基于时延修正的海面高度测量装置。该装置不但可以精确反演海面风速、风向,还可以精确反演海面高度。反演海面风场时,可以得到与实测相关功率曲线相匹配的理论相关功率曲线。此理论相关功率曲线是关于时间延迟的函数,可从中估计出反射信号峰值点相对于镜面反射点处的时间延迟,即时延误差。本发明利用此时延误差对海面测高中的时延进行修正,由海面测高的实际时间延迟减去此时延误差得到时间差,将时间差乘以光速,得到路径差,然后再根据GNSS卫星、接收机和镜面反射点构成的几何关系来完成高度测量。
本发明采取以下技术方案:
本发明一种基于时延修正的海面高度测量装置,它是由一个导航接收机、一个右旋天线、一个左旋天线、一个射频前端模块、一个中频信号处理模块和一个反演模块组成。它们之间的位置连接关系是:导航接收机分别与右旋天线、中频信号处理模块和反演模块相连接,右旋天线和左旋天线彼此独立,并一起与射频前端模块、中频信号处理模块和反演模块依次相连。它们之间的信号走向是:右旋天线输出的直射射频电压信号DRF和左旋天线输出的反射射频电压信号RRF作为射频前端模块的输入,射频前端模块输出的直射数字中频信号DDRF和反射数字中频信号RDIF作为中频信号处理模块的输入,中频信号处理模块输出的时延多普勒二维相关功率DDM(τ,f)作为反演模块的输入,右旋天线输出的直射射频电压信号DRF作为导航接收机的输入,导航接收机输出的导航信息NAV_inf作为中频信号处理模块和反演模块的输入。
所述导航接收机用于输出包含接收机的高度和速度、卫星高度角、方位角等信息的导航信息NAV_inf,作为中频信号处理模块和反演模块的输入。
所述右旋天线采用商用GNSS右旋天线用于接收导航卫星的直射电磁信号,并将该信号转化成直射射频电压信号DRF作为射频前端处理的输入。
所述左旋天线采用多阵元阵列式左旋天线用于接收导航卫星的反射电磁信号,并将该信号转化成反射射频电压信号RRF作为射频前端处理的输入。
所述射频前端模块将接收到的直射射频电压信号DRF和反射射频电压信号RRF下变频为中频信号,并将中频信号转化为直射数字中频信号DDRF和反射数字中频信号RDIF,作为中频信号处理模块的输入。
所述中频信号处理模块由直射通道和反射通道组成。两者之间的关系是:彼此相连。中频信号处理模块作为反演模块的输入。
该直射通道由可见星预判、捕获和跟踪三部分组成。三者之间的关系是:可见星预判、捕获和跟踪依次相连。可见星预判对左旋天线视野范围内的卫星进行预测,减小直射通道卫星搜索范围,输出左旋天线视野范围内的卫星PRN号,记为(PRN1,...,PRNn),作为捕获的输入;捕获对(PRN1,...,PRNn)内的卫星进行搜索,判断卫星是否真实可见,若真实可见,则对直射信号码相位和多普勒频率进行粗略估计,作为跟踪的输入;跟踪对捕获成功的直射信号进行精确的码相位和多普勒估计,记为(τdi,fdi),作为反射通道的输入。
该反射通道由选星、补偿量计算、载波发生、伪码发生和相关处理五部分组成。五者之间的关系是:选星与补偿量计算相连,补偿量计算分别与载波发生和伪码发生相连,载波发生和伪码发生彼此独立,载波发生和伪码发生都与相关处理相连。选星对跟踪到的卫星按照卫星高度角进行选星,选择卫星高度角最大的卫星,其直射信号对应码相位和多普勒估值记为(τd,fd),作为补偿量计算的输入;补偿量计算根据导航接收机输出的导航信息NAV_inf计算得到补偿量Δτ和Δf,由直射信号对应码相位和多普勒估值(τd,fd)得到反射通道的本地伪码相位τd+Δτ和本地载波频率fd+Δf,作为伪码发生和载波发生的输入;伪码发生产生本地伪码用于反射信号伪码剥离;载波发生产生本地载波用于反射信号载波剥离;相关处理将本地载波和本地伪码与反射信号RDIF进行相关运算,输出时延多普勒二维相关功率DDM(τ,f)。
所述反演模块由海面风场反演模块和海面高度反演模块组成。两者之间的关系是:彼此相连。
该海面风场反演模块由理论模型库生成、降噪处理、归一化、波形匹配和残差平方和计算五部分组成。五者之间的关系是:理论模型库生成、降噪处理、归一化、波形匹配和残差平方和计算依次相连。理论模型库生成是从DDM(τ,f)和导航接收机输出的导航信息NAV_inf中提取接收机的高度、速度、卫星高度角和方位角等信息,将提取信息输入到海面散射信号相关功率模型中,输入初始风速v0,计算不同风速条件下的海面散射信号相关功率,生成理论相关功率波形库;降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响;归一化是对降噪处理后的实测相关功率和生成的理论相关功率波形库分别进行归一化处理;波形匹配是将归一化后的实测相关功率曲线和理论相关功率曲线采用最小二乘法进行波形匹配,输出与实测相关功率曲线匹配成功的理论相关功率曲线所对应的风速v;残差平方和计算是由给定的初始风向φ0计算多颗导航卫星的理论和实测相关功率波形的残差平方和,调整风向,输出残差平方和局部最小值对应的风向φ。
该海面高度反演模块由降噪处理、归一化、时延计算,时延修正和高度计算五部分组成。五者之间的关系是:降噪处理、归一化、时延计算、时延修正和高度计算依次相连。降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响;归一化是对降噪处理后的实测相关功率进行归一化处理;时延计算是由直射信号相关功率峰值点处码相位和反射信号相关功率峰值点处码相位计算反射信号相对于直射信号的时延值τE1;时延修正是由匹配成功的理论相关功率曲线估计海面粗糙度引起的时延误差τC,由此时延值来修正海面测高中的时延差,由修正后的时延值τE求出直射信号与反射信号的路径延迟;高度计算是从导航信息NAV_inf中提取导航接收机的位置和卫星高度角等信息,由导航接收机位置可求得导航接收机到地心的距离,从而可计算海面高度H。
本发明的原理及工况简介如下:
在利用GNSS-R反演海面风场时,从匹配得到的理论相关功率曲线中得到反射信号峰值点相对于镜面反射点的时间延迟,用此延迟值来修正反演海面高度中由于海面粗糙度引起的时延误差,从而提高测高精度。右旋天线和左旋天线分别接收直射信号和反射信号并将其转化成射频电压信号,再经射频前端模块和中频信号处理模块输出二维相关功率曲线,由二维相关功率曲线和导航接收机输出的导航信息在反演模块中对海面风场和海面高度进行反演。
本发明的优点在于:
一.该装置仅需接收装置,降低了装置的复杂性和成本。
二.该装置可完成海面风场和海面高度的反演,加强了对海洋的监测能力。
三.该装置用海面风场反演得到的时延误差值直接来修正海面高度测量中的时延误差,减小了海面高度测量中的运算量。
四.该装置可减小由于海面粗糙度引起的误差,提高了海面高度测量精度。
附图说明
图1是本发明测量装置的总体框架图。
图2是本发明测量装置的中频信号处理模块信号处理框图。
图3是本发明测量装置的海面风场和海面高度反演框图。
图中符号说明如下:NAV_inf为导航接收机输出的包含接收机的高度和速度、卫星高度角、方位角等信息的导航信息,DRF为直射射频电压信号,RRF为反射射频电压信号,DDRF为直射数字中频信号,RDIF为反射数字中频信号,DDM(τ,f)为实测时延多普勒二维相关功率曲线,H为需要输出的海面高度,v0为输入的初始海面风速,φ0为输入的初始海面风向,v为需要输出的海面风速,φ为需要输出的海面风向,(PRN1,...,PRNn)为左旋天线视野范围内的卫星PRN号,(τdi,fdi)为跟踪到的卫星的直射信号的码相位和多普勒估计值,(τd,fd)为卫星高度角最大的卫星的直射信号的码相位和多普勒估计值,τd+Δτ为反射通道的本地伪码相位,fd+Δf为反射通道的本地载波频率,τC为海面粗糙度引起的时延误差,τE1为修正前反射信号相对于直射信号的时延值,τE为修正后反射信号相对于直射信号的时延值。
具体实施方式
下面将结合附图对本发明做进一步的详细说明。
如图1为本装置的总体结构图。
本发明装置由一个导航接收机、一个商用GNSS右旋天线、一个多阵元阵列式左旋天线、一个射频前端模块、一个中频信号处理模块和一个反演模块组成。它们之间的位置连接关系是:导航接收机分别与右旋天线、中频信号处理模块和反演模块相连接,右旋天线和左旋天线彼此独立,并一起与射频前端模块、中频信号处理模块和反演模块依次相连。它们之间的信号走向是:右旋天线输出的直射射频电压信号DRF和左旋天线输出的反射射频电压信号RRF作为射频前端模块的输入,射频前端模块输出的直射数字中频信号DDRF和反射数字中频信号RDIF作为中频信号处理模块的输入,中频信号处理模块输出的时延多普勒二维相关功率DDM(τ,f)作为反演模块的输入,右旋天线输出的直射射频电压信号DRF作为导航接收机的输入,导航接收机输出的导航信息NAV_inf作为中频信号处理模块和反演模块的输入。所述导航接收机用于输出包含该接收机的高度和速度、卫星高度角、方位角等信息的导航信息NAV_inf,作为中频信号处理模块和反演模块的输入。该导航接收机输出包含接收机的高度和速度、卫星高度角、方位角等信息的导航信息NAV_inf。
该右旋天线用于接收导航卫星的直射电磁信号,并将该信号转化成直射射频电压信号DRF作为射频前端处理的输入。
该左旋天线用于接收导航卫星的反射电磁信号,并将该信号转化成反射射频电压信号RRF作为射频前端处理的输入。
该射频前端模块包括直射信号射频前端和反射信号射频前端。射频前端模块将接收到的直射射频电压信号DRF和反射射频电压信号RRF下变频为中频信号,并将中频信号转化为直射数字中频信号DDRF和反射数字中频信号RDIF,作为中频信号处理模块的输入。
如图2为本装置的中频信号处理模块的信号处理框图。
中频信号处理模块包括直射通道和反射通道。直射通道主要完成直射信号的捕获和跟踪,并将跟踪得到的时延、多普勒估计值输出作为反射信号处理的参考信息;反射通道主要完成反射信号的相关处理,得到时延多普勒二维相关功率。
直射通道由可见星预判、捕获和跟踪三部分组成。
可见星预判部分接收导航接收机输出的导航信息NAV_inf对左旋天线视野范围内的卫星进行预测,减小直射通道卫星搜索范围,输出左旋天线视野范围内的卫星PRN号,记为(PRN1,...,PRNn);
捕获部分对(PRN1,...,PRNn)内的卫星进行搜索,判断卫星是否真实可见,若真实可见,则对直射信号码相位和多普勒频率进行粗略估计;
跟踪部分对捕获成功的直射信号进行精确的码相位和多普勒估计,记为(τdi,fdi)。
反射通道由选星、补偿量计算、载波发生、伪码发生和相关处理五部分组成。
选星是对跟踪到的卫星按照卫星高度角进行选星,选择卫星高度角最大的卫星,其直射信号对应码相位和多普勒估值为(τd,fd),作为反射信号处理的参考信息。
补偿量计算根据导航接收机输出的导航信息NAV_inf计算得到补偿量Δτ和Δf,由直射信号对应码相位和多普勒估值为(τd,fd)求得反射通道的本地伪码相位τd+Δτ和本地载波频率fd+Δf。
载波发生产生本地载波用于反射信号载波剥离。
伪码发生产生本地伪码用于反射信号伪码剥离。
相关处理将本地载波和本地伪码与反射信号RDIF进行相关运算,输出时延多普勒二维相关功率DDM(τ,f)。
如图3为本装置的海面风场和海面高度反演框图。
海面风场反演由理论模型库生成、降噪处理、归一化、波形匹配和残差平方和计算五部分组成。该模块主要对海面风速、风向进行估计,风场反演过程中由波形匹配可得到理论二维相关功率,从此理论二维相关功率中可得到由于海面粗糙度引起的反射延迟时间。
理论模型库生成是从DDM(τ,f)和导航接收机输出的导航信息NAV_inf中提取接收机的高度、速度、卫星高度角和方位角等信息,将提取信息输入到海面散射信号相关功率模型:
中,输入初始风速v0,计算不同风速条件下的海面散射信号相关功率,生成理论相关功率波形库;
降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响;
归一化是对降噪处理后的实测相关功率和生成的理论相关功率波形库分别进行归一化处理;
波形匹配是将归一化后的实测相关功率曲线和理论相关功率曲线采用最小二乘法进行波形匹配,输出与实测相关功率曲线匹配成功的理论相关功率曲线所对应的风速v;
残差平方和计算是由给定的初始风向φ0计算多颗导航卫星的理论和实测相关功率波形的残差平方和,调整风向,输出残差平方和局部最小值对应的风向φ。
海面高度反演由降噪处理、归一化、时延计算,时延修正和高度计算五部分组成。该模块主要对海面高度进行测量,由反演风场中得到的延迟时间用于修正海面高度测量中的延迟时间,从而精确测量海面高度。
降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响,反演风场时已完成;
归一化是对降噪处理后的实测相关功率进行归一化处理,反演风场时已完成;
时延计算是由直射信号相关功率峰值点处码相位和反射信号相关功率峰值点处码相位计算反射信号相对于直射信号的时延值τE1;
时延修正是由匹配成功的理论相关功率曲线估计海面粗糙度引起的时延误差τC,由此时延值来修正海面测高中的时延差,则修正后的时延值为:τE=τE1-τC。由时延值τE求得直射信号与反射信号的路径延迟为ρE=τE×c;
高度计算是从导航信息NAV_inf中提取导航接收机的位置和卫星高度角等信息,由导航接收机位置求得接收机到地心的距离为h,则海面高度H可由公式求得。
Claims (1)
1.一种基于时延修正的海面高度测量装置,其特征在于:该装置由一个导航接收机、一个商用GNSS右旋天线、一个多阵元阵列式左旋天线、一个射频前端模块、一个中频信号处理模块和一个反演模块组成;它们之间的位置连接关系是:导航接收机分别与右旋天线、中频信号处理模块和反演模块相连接,右旋天线和左旋天线彼此独立,并一起与射频前端模块、中频信号处理模块和反演模块依次相连;它们之间的信号走向是:右旋天线输出的直射射频电压信号DRF和左旋天线输出的反射射频电压信号RRF作为射频前端模块的输入,射频前端模块输出的直射数字中频信号DDRF和反射数字中频信号RDIF作为中频信号处理模块的输入,中频信号处理模块输出的时延多普勒二维相关功率DDM(τ,f)作为反演模块的输入,右旋天线输出的直射射频电压信号DRF作为导航接收机的输入,导航接收机输出的导航信息NAV_inf作为中频信号处理模块和反演模块的输入;所述导航接收机用于输出包含该接收机的高度和速度、卫星高度角、方位角信息的导航信息NAV_inf,作为中频信号处理模块和反演模块的输入;该导航接收机输出包含接收机的高度和速度、卫星高度角、方位角信息的导航信息NAV_inf;
该右旋天线用于接收导航卫星的直射电磁信号,并将该信号转化成直射射频电压信号DRF作为射频前端处理的输入;
该左旋天线用于接收导航卫星的反射电磁信号,并将该信号转化成反射射频电压信号RRF作为射频前端处理的输入;
该射频前端模块包括直射信号射频前端和反射信号射频前端;射频前端模块将接收到的直射射频电压信号DRF和反射射频电压信号RRF下变频为中频信号,并将中频信号转化为直射数字中频信号DDRF和反射数字中频信号RDIF,作为中频信号处理模块的输入;
中频信号处理模块包括直射通道和反射通道;直射通道主要完成直射信号的捕获和跟踪,并将跟踪得到的码相位和多普勒估计值输出作为反射信号处理的参考信息;反射通道主要完成反射信号的相关处理,得到时延多普勒二维相关功率;
直射通道由可见星预判、捕获和跟踪三部分组成;
可见星预判部分接收导航接收机输出的导航信息NAV_inf对左旋天线视野范围内的卫星进行预测,减小直射通道卫星搜索范围,输出左旋天线视野范围内的卫星PRN号,记为PRN1,...,PRNn;
捕获部分对PRN1,...,PRNn内的卫星进行搜索,判断卫星是否真实可见,若真实可见,则对直射信号码相位和多普勒频率进行粗略估计;
跟踪部分对捕获成功的直射信号进行精确的码相位和多普勒估计,记为τdi,fdi;
反射通道由选星、补偿量计算、载波发生、伪码发生和相关处理五部分组成;
选星是对跟踪到的卫星按照卫星高度角进行选星,选择卫星高度角最大的卫星,其直射信号对应码相位和多普勒估值为τd,fd,作为反射信号处理的参考信息;
补偿量计算根据导航接收机输出的导航信息NAV_inf计算得到补偿量Δτ和Δf,由直射信号对应码相位和多普勒估值为τd,fd求得反射通道的本地伪码相位τd+Δτ和本地载波频率fd+Δf;
载波发生产生本地载波用于反射信号载波剥离;
伪码发生产生本地伪码用于反射信号伪码剥离;
相关处理将本地载波和本地伪码与反射信号RDIF进行相关运算,输出时延多普勒二维相关功率DDM(τ,f);
海面风场反演由理论模型库生成、降噪处理、归一化、波形匹配和残差平方和计算五部分组成;该模块主要对海面风速、风向进行估计,风场反演过程中由波形匹配可得到理论二维相关功率,从此理论二维相关功率中可得到由于海面粗糙度引起的反射延迟时间;
理论模型库生成是从DDM(τ,f)和导航接收机输出的导航信息NAV_inf中提取接收机的高度、速度、卫星高度角和方位角的信息,将提取信息输入到海面散射信号相关功率模型:
<mrow>
<mo><</mo>
<mo>|</mo>
<mi>Y</mi>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>,</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>></mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>&Integral;</mo>
<mfrac>
<mrow>
<msup>
<mi>G</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&Lambda;</mi>
<mn>2</mn>
</msup>
<mo>&lsqb;</mo>
<mi>&delta;</mi>
<mi>&tau;</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>4</mn>
<msup>
<msub>
<mi>&pi;R</mi>
<mi>t</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msup>
<msub>
<mi>R</mi>
<mi>r</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&times;</mo>
<mo>|</mo>
<mi>S</mi>
<mo>&lsqb;</mo>
<mi>&delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<msub>
<mi>&sigma;</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
中,输入初始风速v0,计算不同风速条件下的海面散射信号相关功率,生成理论相关功率波形库;
降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响;
归一化是对降噪处理后的实测相关功率和生成的理论相关功率波形库分别进行归一化处理;
波形匹配是将归一化后的实测相关功率曲线和理论相关功率曲线采用最小二乘法进行波形匹配,输出与实测相关功率曲线匹配成功的理论相关功率曲线所对应的风速v;
残差平方和计算是由给定的初始风向φ0计算多颗导航卫星的理论和实测相关功率波形的残差平方和,调整风向,输出残差平方和局部最小值对应的风向φ;
海面高度反演由降噪处理、归一化、时延计算,时延修正和高度计算五部分组成;该模块主要对海面高度进行测量,由反演风场中得到的延迟时间用于修正海面高度测量中的延迟时间,从而精确测量海面高度;
降噪处理是对实测相关功率曲线DDM(τ,f)进行降噪处理,消除噪声影响,反演风场时已完成;
归一化是对降噪处理后的实测相关功率进行归一化处理,反演风场时已完成;
时延计算是由直射信号相关功率峰值点处码相位和反射信号相关功率峰值点处码相位计算反射信号相对于直射信号的时延值τE1;
时延修正是由匹配成功的理论相关功率曲线估计海面粗糙度引起的时延误差τC,由此时延值来修正海面测高中的时延差,则修正后的时延值为:τE=τE1-τC;由时延值τE求得直射信号与反射信号的路径延迟为ρE=τE×c;
高度计算是从导航信息NAV_inf中提取导航接收机的位置和卫星高度角的信息,由导航接收机位置求得接收机到地心的距离为h,则海面高度H可由公式求得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510071278.2A CN104678371B (zh) | 2015-02-11 | 2015-02-11 | 一种基于时延修正的海面高度测量装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510071278.2A CN104678371B (zh) | 2015-02-11 | 2015-02-11 | 一种基于时延修正的海面高度测量装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104678371A CN104678371A (zh) | 2015-06-03 |
CN104678371B true CN104678371B (zh) | 2018-01-09 |
Family
ID=53313699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510071278.2A Active CN104678371B (zh) | 2015-02-11 | 2015-02-11 | 一种基于时延修正的海面高度测量装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104678371B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9906275B2 (en) * | 2015-09-15 | 2018-02-27 | Energous Corporation | Identifying receivers in a wireless charging transmission field |
CN105301622B (zh) * | 2015-11-02 | 2018-06-15 | 北京航大泰科信息技术有限公司 | 一种基于导航卫星反射信号的风速探测装置 |
CN110927221A (zh) * | 2018-10-08 | 2020-03-27 | 中国石油化工股份有限公司 | 一种海面溢油探测装置 |
CN110824512B (zh) * | 2019-11-26 | 2022-01-25 | 中国科学院国家空间科学中心 | 一种非均匀码片实时延迟多普勒映像数据生成器 |
CN112880633A (zh) * | 2021-01-12 | 2021-06-01 | 上海海洋大学 | 基于伯格算法的海面高度测量方法 |
CN113031015A (zh) * | 2021-03-05 | 2021-06-25 | 北京航空航天大学 | 一种基于gnss-r载波相位的海冰检测装置及方法 |
CN113049062B (zh) * | 2021-03-12 | 2022-04-15 | 北京航空航天大学 | 一种gnss直反射信号载波干涉测量湖泊水位的装置 |
CN113671210B (zh) * | 2021-08-19 | 2022-07-26 | 北京航空航天大学 | 陆基gnss直反射信号载波干涉测量河流流速的装置 |
CN116539913B (zh) * | 2023-05-04 | 2024-06-14 | 极诺星空(北京)科技有限公司 | 星上实时反演海面风速的方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位***信号源的机载海洋微波遥感*** |
CN101975959A (zh) * | 2010-09-03 | 2011-02-16 | 北京航空航天大学 | 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置 |
CN102749067A (zh) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | 自适应多级计算的精确反演海面高程遥测装置 |
-
2015
- 2015-02-11 CN CN201510071278.2A patent/CN104678371B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位***信号源的机载海洋微波遥感*** |
CN101975959A (zh) * | 2010-09-03 | 2011-02-16 | 北京航空航天大学 | 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置 |
CN102749067A (zh) * | 2012-06-21 | 2012-10-24 | 北京航空航天大学 | 自适应多级计算的精确反演海面高程遥测装置 |
Non-Patent Citations (3)
Title |
---|
"基于GNSS_R的海面风场监测***研究";路勇等;《武汉大学学报》;20090430;第34卷(第4期);第470-473页 * |
"渤海GNSS-R机载测风试验及其反演结果";周兆明等;《武汉大学学报》;20080331;第33卷(第3期);第241-244页 * |
"面向遥感的GNSS反射信号接收处理***及实验";李伟强等;《武汉大学学报》;20111031;第36卷(第10期);第1204-1208页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104678371A (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104678371B (zh) | 一种基于时延修正的海面高度测量装置 | |
CN104677463B (zh) | 一种基于统计量的gnss‑r多源海面测高装置 | |
US10613197B2 (en) | Antenna specification estimation device and radar device | |
CN107607943B (zh) | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 | |
Li et al. | An improved radon-transform-based scheme of Doppler centroid estimation for bistatic forward-looking SAR | |
US8179317B2 (en) | Method and apparatus for passive geolocation using synthetic-aperture processing | |
CN102749067A (zh) | 自适应多级计算的精确反演海面高程遥测装置 | |
CN112904384B (zh) | 基于机器学习的星载gnss-r高度要素探测装置及方法 | |
US11914061B2 (en) | Systems and methods for locating a signal source | |
CN109856605A (zh) | 一种同时数字多波束形成的二次拟合曲线指向修正方法 | |
CN103376447A (zh) | 一种实现非合作双多基地雷达目标三维定位的方法 | |
CN106646548A (zh) | 一种多模的超快速、低功耗的定位接收机***及方法 | |
CN103760546A (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN103293526B (zh) | 基于合成导向矢量的波束形成米波雷达测高方法 | |
CN104145194A (zh) | 定位信号检测方法、定位信号检测程序、定位信号接收装置、定位装置及信息设备终端 | |
CN105259537A (zh) | 基于频移迭代的多普勒谱中心频率估计方法 | |
CN105044744B (zh) | 一种基于多普勒频移特征的非合作导航卫星识别的方法 | |
CN105182380A (zh) | 一种实现gnss-r相位差提取的硬件接收机及方法 | |
CN108279425A (zh) | 一种多频观测过程中多径误差的修正方法 | |
CN105699992B (zh) | 高动态gnss载波信号跟踪方法及*** | |
CN112880633A (zh) | 基于伯格算法的海面高度测量方法 | |
CN111007490A (zh) | 一种基于浮标地理信息的天波超视距雷达坐标配准方法 | |
CN104360335B (zh) | 基于ais目标指示的电离层参数反演方法 | |
CN109828232A (zh) | 一种阵列天线测角模糊置信度计算装置及方法 | |
CN103940424A (zh) | 一种x射线脉冲星导航信号整周模糊度检测与估计方法 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211029 Address after: 272001 fifth floor, building A3, industry university research base, high tech Zone, Jining City, Shandong Province Patentee after: SHANDONG HANGXIANG ELECTRONIC SCIENCE & TECHNOLOGY Co.,Ltd. Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road Patentee before: BEIHANG University |
|
TR01 | Transfer of patent right |