CN101975959B - 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置 - Google Patents

适用于延迟映射接收机中对反演有效波高和相对高程的获取装置 Download PDF

Info

Publication number
CN101975959B
CN101975959B CN2010102715763A CN201010271576A CN101975959B CN 101975959 B CN101975959 B CN 101975959B CN 2010102715763 A CN2010102715763 A CN 2010102715763A CN 201010271576 A CN201010271576 A CN 201010271576A CN 101975959 B CN101975959 B CN 101975959B
Authority
CN
China
Prior art keywords
signal
module
direct
wave height
carries out
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.)
Expired - Fee Related
Application number
CN2010102715763A
Other languages
English (en)
Other versions
CN101975959A (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN2010102715763A priority Critical patent/CN101975959B/zh
Publication of CN101975959A publication Critical patent/CN101975959A/zh
Application granted granted Critical
Publication of CN101975959B publication Critical patent/CN101975959B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置,该装置包括直射信号捕获、跟踪及定位模块、二维相关计算模块、有效波高计算模块和相对高程计算模块。在本发明中,采用频域处理方法对直射数字中频信号DDIF和反射数字中频信号RDIF进行处理,得到卫星直射信号和海面反射信号的二维相关值CDCo、CRCo和二维相关功率值PDCo、PRCo;采用干涉复数场处理方法对CDCo、CRCo和ε进行处理,得到海面有效波高SWH;采用峰值跟随处理方法对PDCo、PRCo和ε进行处理,得到延迟映射接收机至海面的相对高程H。

Description

适用于延迟映射接收机中对反演有效波高和相对高程的获取装置
技术领域
本发明涉及一种海洋遥感装置,更特别地说,是一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置。
背景技术
全球定位系统不仅为空间信息用户提供了导航定位信息,还提供了源源不断的L波段微波信号资源。最近十几年,GPS-R遥感技术得到了迅速发展,通过接收卫星直射信号和海面反射信号,采用相关处理,得到C/A码的相关函数波形,进而进行要素反演,可以用于海面风场反演、海洋盐度测量、海面高程测量、地面湿度区分、地形高度预警等方面。
现有延迟映射接收机或多普勒延迟映射接收机都是由左旋接收天线、右旋接收天线、双射频前端、高速A/D采样器、多通道专用相关器、基带信号处理器、上位机以及相应的接口电路组成,而多通道专用相关器和基带信号处理器一般多采用DSP处理芯片实现,如《测控技术》2006年第25卷第5期中图2所示的延迟映射接收机的系统构成。对反射信号的处理都采用多通道相关器这种时域处理方法,而且大都采取事后处理的方法,而不是实时输出反演要素。对于延迟映射接收机或多普勒延迟映射接收机来讲,直射信号顺次经右旋天线、射频前端、高速A/D转换器后输出直射数字中频信号DDIF;反射信号顺次经左旋天线、射频前端、高速A/D转换器后输出反射数字中频信号RDIF。
快速傅里叶变换(FFT)在处理GPS直射信号的捕获上有所应用,但在处理GPS反射信号上并没有应用,随着大规模FPGA芯片的不断发展,使得FFT在GPS反射信号的处理上能够很好的应用,它能够提高GPS反射信号的处理速度。
发明内容
本发明提出了一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置,该装置通过在FPGA芯片用Verilog HDL语言设计得到硬件的功效。在本发明中,采用频域处理方法对直射数字中频信号DDIF和反射数字中频信号RDIF进行处理,得到卫星直射信号和海面反射信号的二维相关值(CDCo、CRCo)和二维相关功率值(PDCo、PRCo);采用干涉复数场处理方法对CDCo、CRCo和ε进行处理,得到海面有效波高SWH;采用峰值跟随处理方法对PDCo、PRCo和ε进行处理,得到延迟映射接收机至海面的相对高程H。
本发明的一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置,该装置包括直射信号捕获、跟踪及定位模块、二维相关计算模块、有效波高计算模块和相对高程计算模块;
直射信号捕获、跟踪及定位模块用于产生载波频率CANO、码相位CDNO、导航定位解Psv和卫星高度角ε;其中,载波频率CANO、码相位CDNO、导航定位解Psv输出给二维相关计算模块,卫星高度角ε分别输出给有效波高计算模块和相对高程计算模块;
二维相关计算模块中采用频域处理方法的具体步骤如下:
直射数字中频信号DDIF输入到数据存储A中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率信号CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储A中存储1ms直射信号数字中频信号DDIF进行乘法运算,得到去除了载波的直射基带信号水平分量IDB;
所述垂直信号分量QC与数据存储A中存储1ms直射信号数字中频信号DDIF进行乘法运算,得到去除了载波的直射信号垂直分量QDB;
变换器A将直射基带信号水平分量IDB和直射信号垂直分量QDB合成复数形式后进行FFT变换输出直-频域信号FIQDB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位信号CDNO输入到数据存储C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述直-频域信号FIQDB与码相位频域信号FIQC进行乘法运算后,输出结果直-频域相关信号FD;
所述直-频域相关信号FD在逆变换器A中进行进行快速傅里叶逆变换处理得到直-时域相关信号TD输入到数据存储D;
在截取器A中所述的直-时域相关信号TD受到码延迟相位控制CP,截取设定延迟范围内的直-二维相关信号TDC;
在相干叠加器A中对直-二维相关信号TDC进行相干叠加运算得到叠加后的直-二维相干信号CoD;
所述的直-二维相干信号CoD一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器A中进行相关功率计算和非相干叠加计算输出直-二维相关功率信号InD;该信号InD输出到相对高程计算模块反演相对高程H;
反射数字中频信号RDIF输入到数据存储B中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率信号CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储B中存储1ms反射信号数字中频信号RDIF进行乘法运算,得到去除了载波的反射基带信号水平分量IRB;
所述垂直信号分量QC与数据存储B中存储1ms反射信号数字中频信号RDIF进行乘法运算,得到去除了载波的反射信号垂直分量QRB;
变换器B将反射基带信号水平分量IRB和反射信号垂直分量QRB合成复数形式后进行FFT变换输出反-频域信号FIQRB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位信号CDNO输入到数据存储C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述反-频域信号FIQRB与码相位频域信号FIQC进行乘法运算后,输出结果反-频域相关信号FR;
所述反-频域相关信号FR在逆变换器B中进行进行快速傅里叶逆变换处理得到反-时域相关信号TR输入到数据存储E;
在截取器B中所述的反-时域相关信号TR一方面受到码延迟相位控制CP,另一方面受到反射信号延时控制RCP,从而截取设定延迟范围内的反-二维相关信号TRC;所述反射信号延时时间控制字RCP是对导航定位解Psv(直射信号捕获、跟踪及定位模块输出的)根据海面发生反射的空间几何关系RCP=2hsinε得到,h表示延迟映射接收机所在环境下的海拔高度,ε表示反射卫星的高度角;
在相干叠加器B中对反-二维相关信号TRC进行相干叠加运算得到叠加后的反-二维相干信号CoR;
所述的反-二维相干信号CoR一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器B中进行相关功率计算和非相干叠加计算输出反-二维相关功率信号InR;该信号InR输出到相对高程计算模块反演相对高程H;
有效波高计算模块包括有最大值判断模块MaxA、最大值判断模块MaxB、除法器、存储到数据存储F、有效相关时间计算模块和有效波高计算模型模块;
在最大值判断模块MaxA中对直-二维相干信号CoD进行最大值选取,得到直-最大二维相关值MCoD;
在最大值判断模块MaxB中对反-二维相干信号CoR进行最大值选取,得到反-最大二维相关值MCoR;
在除法器中依据干涉复数场关系IFC=MCoR/MCoD进行处理,得到干涉复数场IFC;
在存储到数据存储F中对干涉复数场IFC进行5s数据存储后,然后在自相关运算模块中进行自相关运算得到自相函数值FIFC;
在有效相关时间计算模块中先对自相关函数值FIFC进行相关时间τ的判断,然后根据有效相关时间关系τ′=τ×sin(ε)获得有效相关时间τ′;
在有效波高计算模型模块中依据有效波高经验值关系SWH=a/τ′+b得到有效波高SWH;a=0.0182,b=-0.5647;
相对高程计算模块包括有峰值位置检测模块A、峰值位置检测模块B、峰值路径差计算模块和高程计算模块;
在峰值位置检测模块A中对直-二维相关功率信号InD进行峰值功率出现位置判断,得到直-峰值功率最大值延迟路径DPD;
在峰值位置检测模块B中对反-二维相关功率信号InR进行峰值功率出现位置判断,得到反-峰值功率最大值延迟路径RPD;
在峰值路径差计算模块中,根据峰值延迟关系ρ=RPD-DPD计算反射信号与直射信号的峰值路径差ρ;
在高程计算模块中,依据相对高程解析公式H=ρ/(2×sin(ε))计算延迟映射接收机到海面的相对高程H。
适用于延迟映射接收机中对反演有效波高和相对高程的获取装置的特点:
1、本发明可以实时输出海面有效波高;
2、本发明可以实时输出装置到海面的相对高程;
3、本发明的二维相关计算模块对直射信号和海面反射信号的处理采用了快速傅里叶变换的处理方式,处理速度快。
附图说明
图1为本发明反演有效波高和相对高程海洋遥感装置的结构框图。
图2为本发明的二维相关计算模块的结构框图。
图3为本发明的有效波高计算模块的结构框图。
图4为本发明的相对高程模块的结构框图。
具体实施方式
下面将结合附图对本发明做进一步的详细说明。
参见图1所示,本发明是一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置,该装置包括直射信号捕获、跟踪及定位模块、二维相关计算模块、有效波高计算模块和相对高程计算模块。
(一)直射信号捕获、跟踪及定位模块
直射信号捕获、跟踪及定位模块采用2006年7月第27卷第4期在《遥测遥控》公开的“GPS接收机基带信号处理模块的FPGA实现”文献进行处理,分别得到载波频率CANO、码相位CDNO、导航定位解Psv和卫星高度角ε。其中,载波频率CANO、码相位CDNO、导航定位解Psv输出给二维相关计算模块,卫星高度角ε分别输出给有效波高计算模块和相对高程计算模块。
(二)二维相关计算模块
参见图2所示,二维相关计算模块利用FFT的方法对直射数字中频信号DDIF和反射数字中频信号RDIF进行处理得到直射信号和反射信号的二维相关值,并将直射信号和反射信号的二维相关值一方面输出到有效波高计算模块反演有效波高H(单位米);另一方面输出到相对高程计算模块反演遥感装置到海面的相对高程SWH(单位米)。
二维相关计算模块中采用频域处理方法的具体步骤如下:
直射数字中频信号DDIF输入到数据存储A中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率信号CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储A中存储1ms直射信号数字中频信号DDIF进行乘法运算,得到去除了载波的直射基带信号水平分量IDB;
所述垂直信号分量QC与数据存储A中存储1ms直射信号数字中频信号DDIF进行乘法运算,得到去除了载波的直射信号垂直分量QDB;
变换器A将直射基带信号水平分量IDB和直射信号垂直分量QDB合成复数形式后进行FFT变换输出直-频域信号FIQDB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位信号CDNO输入到数据存储C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述直-频域信号FIQDB与码相位频域信号FIQC进行乘法运算后,输出结果直-频域相关信号FD;
所述直-频域相关信号FD在逆变换器A中进行进行快速傅里叶逆变换处理得到直-时域相关信号TD输入到数据存储D;
在截取器A中所述的直-时域相关信号TD受到码延迟相位控制CP,截取设定延迟范围内的直-二维相关信号TDC;
在相干叠加器A中对直-二维相关信号TDC进行相干叠加运算得到叠加后的直-二维相干信号CoD;
所述的直-二维相干信号CoD一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器A中进行相关功率计算和非相干叠加计算输出直-二维相关功率信号InD;该信号InD输出到相对高程计算模块反演相对高程H。
反射数字中频信号RDIF输入到数据存储B中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率信号CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储B中存储1ms反射信号数字中频信号RDIF进行乘法运算,得到去除了载波的反射基带信号水平分量IRB;
所述垂直信号分量QC与数据存储B中存储1ms反射信号数字中频信号RDIF进行乘法运算,得到去除了载波的反射信号垂直分量QRB;
变换器B将反射基带信号水平分量IRB和反射信号垂直分量QRB合成复数形式后进行FFT变换输出反-频域信号FIQRB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位信号CDNO输入到数据存储C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述反-频域信号FIQRB与码相位频域信号FIQC进行乘法运算后,输出结果反-频域相关信号FR;
所述反-频域相关信号FR在逆变换器B中进行进行快速傅里叶逆变换处理得到反-时域相关信号TR输入到数据存储E;
在截取器B中所述的反-时域相关信号TR一方面受到码延迟相位控制CP,另一方面受到反射信号延时控制RCP,从而截取设定延迟范围内的反-二维相关信号TRC;所述反射信号延时时间控制字RCP是对导航定位解Psv(直射信号捕获、跟踪及定位模块输出的)根据海面发生反射的空间几何关系RCP=2hsinε得到,h表示延迟映射接收机所在环境下的海拔高度,ε表示反射卫星的高度角;
在相干叠加器B中对反-二维相关信号TRC进行相干叠加运算得到叠加后的反-二维相干信号CoR;
所述的反-二维相干信号CoR一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器B中进行相关功率计算和非相干叠加计算输出反-二维相关功率信号InR;该信号InR输出到相对高程计算模块反演相对高程H。
在本发明中,采用频域处理方法对直射信号和反射信号进行二维相关处理,能够实时快速的得到码相位的时域相关值TD和TR。
在本发明中,通过CP参与直射信号、CP和RCP参与反射信号的协作,可以对有效的相关值进行截取,对截取的相关值进行相干叠加提高了信噪比,使得有效波高SWH的反演精度得以提高。
(三)有效波高计算模块
参见图3所示,有效波高计算模块包括有最大值判断模块MaxA、最大值判断模块MaxB、除法器、存储到数据存储F、有效相关时间计算模块和有效波高计算模型模块。
在最大值判断模块MaxA中对直-二维相干信号CoD进行最大值选取,得到直-最大二维相关值MCoD;
在最大值判断模块MaxB中对反-二维相干信号CoR进行最大值选取,得到反-最大二维相关值MCoR;
在除法器中依据干涉复数场关系IFC=MCoR/MCoD进行处理,得到干涉复数场IFC-;
在存储到数据存储F中对干涉复数场IFC进行5s数据存储后,然后在自相关运算模块中进行自相关运算得到自相函数值FIFC;
在有效相关时间计算模块中先对自相关函数值FIFC进行相关时间τ(单位秒)的判断,然后根据有效相关时间关系τ′=τ×sin(ε)获得有效相关时间τ′(单位秒);ε表示反射卫星的高度角;
在有效波高计算模型模块中依据有效波高经验值关系SWH=a/τ′+b得到有效波高SWH;a和b为经验值,a=0.0182(单位米),b=-0.5647(单位米)。
在本发明中,有效波高计算模块采用干涉复数场关系IFC=MCoR/MCoD能够消除直射信号变化带来的干扰,从而提高了有效波高SWH的测量精度。
(四)相对高程计算模块
参见图4所示,相对高程计算模块包括有峰值位置检测模块A、峰值位置检测模块B、峰值路径差计算模块和高程计算模块。
在峰值位置检测模块A中对直-二维相关功率信号InD进行峰值功率出现位置判断,得到直-峰值功率最大值延迟路径DPD;
在峰值位置检测模块B中对反-二维相关功率信号InR进行峰值功率出现位置判断,得到反-峰值功率最大值延迟路径RPD;
在峰值路径差计算模块中,根据峰值延迟关系ρ=RPD-DPD计算反射信号与直射信号的峰值路径差ρ;
在高程计算模块中,依据相对高程解析公式H=ρ/(2×sin(ε))计算延迟映射接收机到海面的相对高程H。
在本发明中,相对高程计算模块中应用直射信号的峰值位置检测和反射信号的峰值位置检测共同作用下,获得了精度较高的峰值路径延迟,从而提高了相对高程H的反演精度。

Claims (1)

1.一种适用于延迟映射接收机中对反演有效波高和相对高程的获取装置,该装置包括直射信号捕获、跟踪及定位模块,其特征在于:还包括有二维相关计算模块、有效波高计算模块和相对高程计算模块;
直射信号捕获、跟踪及定位模块用于产生直射信号载波频率CANO、直射信号码相位CDNO、导航定位解Psv和卫星的高度角ε;其中,直射信号载波频率CANO、直射信号码相位CDNO、导航定位解Psv输出给二维相关计算模块,卫星的高度角ε分别输出给有效波高计算模块和相对高程计算模块;
二维相关计算模块中采用频域处理方法的具体步骤如下:
直射数字中频信号DDIF输入到数据存储模块A中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储模块A中存储1ms直射数字中频信号DDIF进行乘法运算,得到去除了载波的直射基带信号水平分量IDB;
所述垂直信号分量QC与数据存储模块A中存储1ms直射数字中频信号DDIF进行乘法运算,得到去除了载波的直射信号垂直分量QDB;
变换器A将直射基带信号水平分量IDB和直射信号垂直分量QDB合成复数形式后进行FFT变换输出直-频域信号FIQDB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位CDNO输入到数据存储模块C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述直-频域信号FIQDB与码相位频域信号FIQC进行乘法运算后,输出结果直-频域相关信号FD;
所述直-频域相关信号FD在逆变换器A中进行快速傅里叶逆变换处理得到直-时域相关信号TD输入到数据存储模块D;
在截取器A中所述的直-时域相关信号TD受到码延迟相位控制CP,截取设定延迟范围内的直-二维相关信号TDC;
在相干叠加器A中对直-二维相关信号TDC进行相干叠加运算得到叠加后的直-二维相干信号CoD;
所述的直-二维相干信号CoD一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器A中进行相关功率计算和非相干叠加计算输出直-二维相关功率信号InD;该信号InD输出到相对高程计算模块反演相对高程H;
反射数字中频信号RDIF输入到数据存储模块B中存储1ms;
直射信号捕获、跟踪及定位模块输出的直射信号载波频率CANO输入到载波发生器中,载波发生器受多普勒频移dop控制,输出附加多普勒频移的载波信号CANO1,该CANO1信号分别经正弦sin和余弦cos运算后,输出水平分量IC、垂直信号分量QC;
所述水平分量IC与数据存储模块B中存储1ms反射数字中频信号RDIF进行乘法运算,得到去除了载波的反射基带信号水平分量IRB;
所述垂直信号分量QC与数据存储模块B中存储1ms反射数字中频信号RDIF进行乘法运算,得到去除了载波的反射信号垂直分量QRB;
变换器B将反射基带信号水平分量IRB和反射信号垂直分量QRB合成复数形式后进行FFT变换输出反-频域信号FIQRB;
直射信号捕获、跟踪及定位模块输出的直射信号码相位CDNO输入到数据存储模块C中存储1ms;存储1ms的CDNO在变换器C中进行FFT变换和取共轭操作得到码相位频域信号FIQC;
所述反-频域信号FIQRB与码相位频域信号FIQC进行乘法运算后,输出结果反-频域相关信号FR;
所述反-频域相关信号FR在逆变换器B中进行快速傅里叶逆变换处理得到反-时域相关信号TR输入到数据存储模块E;
在截取器B中所述的反-时域相关信号TR一方面受到码延迟相位控制CP,另一方面受到反射信号延时控制字RCP,从而截取设定延迟范围内的反-二维相关信号TRC;所述反射信号延时控制字RCP是对导航定位解Psv根据海面发生反射的空间几何关系RCP=2hsinε得到,h表示延迟映射接收机所在环境下的海拔高度,ε表示卫星的高度角;
在相干叠加器B中对反-二维相关信号TRC进行相干叠加运算得到叠加后的反-二维相干信号CoR;
所述的反-二维相干信号CoR一方面输出到有效波高计算模块反演得到有效波高SWH;另一方面输出到平方和与非相干叠加器B中进行相关功率计算和非相干叠加计算输出反-二维相关功率信号InR;该信号InR输出到相对高程计算模块反演相对高程H;
有效波高计算模块包括有最大值判断模块MaxA、最大值判断模块MaxB、除法器、数据存储模块F、有效相关时间计算模块和有效波高计算模型模块;
在最大值判断模块MaxA中对直-二维相干信号CoD进行最大值选取,得到直-最大二维相关值MCoD;
在最大值判断模块MaxB中对反-二维相干信号CoR进行最大值选取,得到反-最大二维相关值MCoR;
在除法器中依据干涉复数场关系ICF=MCoR/MCoD进行处理,得到干涉复数场ICF;
在数据存储模块F中对干涉复数场ICF进行5s数据存储后,然后在自相关运算模块中进行自相关运算得到自相关函数值FICF;
在有效相关时间计算模块中先对自相关函数值FICF进行相关时间τ的判断,然后根据有效相关时间关系τ′=τ×sin(ε)获得有效相关时间τ′;ε表示卫星的高度角;相关时间τ的单位为秒,有效相关时间τ′的单位为秒;
在有效波高计算模型模块中依据有效波高经验值关系SWH=a/τ′+b得到有效波高SWH;a和b为经验值,a=0.0182米,b=-0.5647米;有效波高SWH的单位为米;
相对高程计算模块包括有峰值位置检测模块A、峰值位置检测模块B、峰值路径差计算模块和高程计算模块;
在峰值位置检测模块A中对直-二维相关功率信号InD进行峰值功率出现位置判断,得到直-峰值功率最大值延迟路径DPD;
在峰值位置检测模块B中对反-二维相关功率信号InR进行峰值功率出现位置判断,得到反-峰值功率最大值延迟路径RPD;
在峰值路径差计算模块中,根据峰值延迟关系ρ=RPD-DPD计算反射信号与直射信号的峰值路径差ρ;
在高程计算模块中,依据相对高程解析公式H=ρ/(2×sin(ε))计算延迟映射接收机到海面的相对高程H,相对高程H的单位为米。
CN2010102715763A 2010-09-03 2010-09-03 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置 Expired - Fee Related CN101975959B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102715763A CN101975959B (zh) 2010-09-03 2010-09-03 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102715763A CN101975959B (zh) 2010-09-03 2010-09-03 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置

Publications (2)

Publication Number Publication Date
CN101975959A CN101975959A (zh) 2011-02-16
CN101975959B true CN101975959B (zh) 2012-07-25

Family

ID=43575861

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102715763A Expired - Fee Related CN101975959B (zh) 2010-09-03 2010-09-03 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置

Country Status (1)

Country Link
CN (1) CN101975959B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749067A (zh) * 2012-06-21 2012-10-24 北京航空航天大学 自适应多级计算的精确反演海面高程遥测装置
CN103869342B (zh) * 2014-04-01 2016-08-10 北京航空航天大学 基于bap的gnss-r延迟映射软件接收机
CN104678371B (zh) * 2015-02-11 2018-01-09 北京航空航天大学 一种基于时延修正的海面高度测量装置
CN104677463B (zh) * 2015-02-11 2017-10-31 北京航空航天大学 一种基于统计量的gnss‑r多源海面测高装置
CN104765032B (zh) * 2015-04-21 2017-02-22 山东省科学院海洋仪器仪表研究所 岸基阵列gnss反射信号潮汐与多波浪参数综合探测系统
CN106959454B (zh) * 2017-03-20 2019-06-04 上海航天控制技术研究所 一种基于数字域tdi及连续复线阵成像模式的颤振反演方法
CN109521445A (zh) * 2018-12-20 2019-03-26 中国气象局气象探测中心 一种基于gnss-r设备的有效波高处理方法及系统
CN110879402B (zh) * 2019-11-28 2021-11-09 中国科学院国家空间科学中心 一种gnss干涉测高中消除直流分量的系统及方法
CN113296136B (zh) * 2020-02-24 2024-03-19 香港理工大学深圳研究院 多目标形变监测方法、装置和接收机
CN111580080B (zh) * 2020-05-29 2023-10-03 大连海事大学 一种基于船载gnss-r的有效波高反演方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1279970A2 (en) * 2001-05-31 2003-01-29 Starlab Barcelona SL A system for monitoring a feature of a surface with broad swath and high resolution
CN2700895Y (zh) * 2004-04-29 2005-05-18 杨东凯 基于全球卫星定位系统的延迟映射接收机
CN101246212A (zh) * 2008-03-18 2008-08-20 北京航空航天大学 一种可接收处理导航卫星反射信号的相关装置

Family Cites Families (1)

* 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

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1279970A2 (en) * 2001-05-31 2003-01-29 Starlab Barcelona SL A system for monitoring a feature of a surface with broad swath and high resolution
CN2700895Y (zh) * 2004-04-29 2005-05-18 杨东凯 基于全球卫星定位系统的延迟映射接收机
CN101246212A (zh) * 2008-03-18 2008-08-20 北京航空航天大学 一种可接收处理导航卫星反射信号的相关装置

Also Published As

Publication number Publication date
CN101975959A (zh) 2011-02-16

Similar Documents

Publication Publication Date Title
CN101975959B (zh) 适用于延迟映射接收机中对反演有效波高和相对高程的获取装置
CN104833972B (zh) 一种双基地调频连续波合成孔径雷达频率变标成像方法
CN102305944A (zh) 一种水下淤积量探测分析处理系统
CN106644030A (zh) 一种基于多普勒雷达的非接触式振动测量方法
CN104502898B (zh) 将修正rft和修正mdcft相结合的机动目标参数估计方法
CN102914768A (zh) 基于形态成分分析的海面微动目标检测和特征提取方法
CN102608587B (zh) 基于非线性最小二乘的空中机动目标检测方法
CN102749067A (zh) 自适应多级计算的精确反演海面高程遥测装置
CN103760546A (zh) 一种雷达用低空目标波达方向估计方法
CN101592733A (zh) 一种逆合成孔径雷达并行实时成像处理方法
Previsic et al. In-ocean validation of a deterministic sea wave prediction (DSWP) system leveraging X-band radar to enable optimal control in wave energy conversion systems
Gebbie et al. A two-hydrophone range and bearing localization algorithm with performance analysis
CN103293521A (zh) 一种利用x波段雷达探测近海海域水深的方法
CN102435985B (zh) 基于重构时间采样的机载机动目标检测方法
CN102141611B (zh) 斜视合成孔径雷达多普勒模糊数快速测定方法
Song et al. Study on the exploration of spaceborne GNSS-R raw data focusing on altimetry
CN102841337A (zh) 一种天波超视距雷达回波信号去非线性相位污染方法
CN104199008A (zh) 基于压缩感知的空中机动目标参数估计方法
CN102801459A (zh) 多模卫星信号接收装置及其方法
CN102890270B (zh) 固定站双基地合成孔径雷达回波模拟方法
CN106066472B (zh) 一种二维振速梯度水听器的被动目标相关检测方法
Pan et al. A raw echo simulation method based on reference signal translation for InSAS
Li et al. Real-time processing of reflected GNSS signals for remote sensing: System and experiments
CN102193096A (zh) 信号捕获方法及信号捕获装置
Xu et al. The study on retrieval technique of significant wave height using airborne GNSS-R

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120725

Termination date: 20150903

EXPY Termination of patent right or utility model