CN116699666B - 星下点观测的星载gnss-r海面风场反演方法及系统 - Google Patents

星下点观测的星载gnss-r海面风场反演方法及系统 Download PDF

Info

Publication number
CN116699666B
CN116699666B CN202310979097.4A CN202310979097A CN116699666B CN 116699666 B CN116699666 B CN 116699666B CN 202310979097 A CN202310979097 A CN 202310979097A CN 116699666 B CN116699666 B CN 116699666B
Authority
CN
China
Prior art keywords
satellite
sea surface
surface wind
coordinate system
gnss
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
Application number
CN202310979097.4A
Other languages
English (en)
Other versions
CN116699666A (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
Shanghai Aerospace System Engineering Institute
Original Assignee
Beihang University
Shanghai Aerospace System Engineering Institute
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, Shanghai Aerospace System Engineering Institute filed Critical Beihang University
Priority to CN202310979097.4A priority Critical patent/CN116699666B/zh
Publication of CN116699666A publication Critical patent/CN116699666A/zh
Application granted granted Critical
Publication of CN116699666B publication Critical patent/CN116699666B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/52—Determining velocity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00—Indicating or recording presence, absence, or direction, of movement
    • G01P13/02—Indicating direction only, e.g. by weather vane
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/62—Sense-of-movement determination
    • 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/393—Trajectory determination or predictive tracking, e.g. Kalman filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种星下点观测的星载GNSS‑R海面风场反演方法及系统，涉及海面风矢量探测技术领域，方法包括：建立星载环境下的坐标系；其中，坐标系的原点位于镜面反射点；坐标系的Z轴指向镜面反射点切平面的法向；导航卫星与GNSS‑R接收卫星均位于坐标系YOZ平面的Z轴正半轴侧；导航卫星位于坐标系的Y轴正半轴侧，GNSS‑R接收卫星位于坐标系的Y轴负半轴侧；在坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量，并计算海面风矢量的数值解，具有反演精度高，算法复杂度低等特点。

Description

星下点观测的星载GNSS-R海面风场反演方法及系统
技术领域
本发明涉及海面风矢量探测技术领域，特别是涉及一种星下点观测的星载GNSS-R海面风场反演方法及系统。
背景技术
海洋动力灾害是对全球沿海各国危害最大的自然灾害。因此，准确探测海面风矢量对海洋动力灾害的预报预警起着至关重要的作用，传统的海面风矢量星载探测手段主要包括辐射计、高度计和散射计等，具有设备较复杂、成本较高等特点，不适合大规模卫星组网遥感探测，因此其空间覆盖性和时间分辨率较低。
全球导航卫星系统反射信号(Global Navigation Satellite SystemReflectometry，GANN-R)技术通过采用岸基、机载以及空载的特殊接受设备接受GNSS直射信号和以及经反射面反射的GNSS反射信号，对GNSS直射信号和GNSS反射信号进行处理来反演目标物的特性。目前星载GNSS-R应用较为广泛，但尚未见星载GNSS-R海面风矢量直接反演的相关研究论述。
发明内容
本发明的目的是提供一种星下点观测的星载GNSS-R海面风场反演方法及系统，可利用三颗方位角不同的导航卫星的星下点非镜向散射信号，反演得到海面风矢量。
为实现上述目的，本发明提供了如下方案：
本发明提供了一种星下点观测的星载GNSS-R海面风场反演方法，包括：
建立星载环境下的坐标系；所述坐标系的原点位于镜面反射点；所述坐标系的Z轴指向镜面反射点切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧；
在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量；
计算所述海面风矢量的数值解。
本发明还提供了一种星下点观测的星载GNSS-R海面风场反演系统，包括：
坐标系构建模块，用于建立星载环境下的坐标系；所述坐标系的原点位于镜面反射点；所述坐标系的Z轴指向镜面反射点切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧；
海面风矢量计算模块，用于在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量；
数值解计算模块，用于计算所述海面风矢量的数值解。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
一、提供了一种基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，可直接利用非镜向独立观测模式下的多颗导航卫星散射信号反演得到海面风矢量（即海面风速和海面风向），反演精度高，算法复杂度低。
二、观测几何关系的对称性和海浪谱的对称性导致海面风向有四个模糊解，在多星观测中可以消除观测几何关系导致的模糊解，将海面风向模糊解由四个降为两个，即增加导航卫星的数量可以消除观测对称性导致的海面风向模糊解。
三、本发明使用的GNSS-R技术，具有信号资源丰富，全球全天时、全天候覆盖、探测范围广、成本低的优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的星下点观测的星载GNSS-R海面风场反演的流程示意图；
图2为本发明实施例提供的星载GNSS-R本地坐标系示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明针对星载GNSS-R海面风矢量反演这一难点，利用非镜向独立观测模式下的多颗导航卫星散射信号反演得到海面风速和海面风向，为星载GNSS-R海面风矢量反演提供了一种适合推广应用的解决方案。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
海面通常可认为是微粗糙表面，GNSS直射信号在海面的散射既包含了镜向散射，也包含了非镜向散射。其中，镜向模式是指镜面反射点位于反射天线的有效覆盖区内，此时接收到的信号主要是镜向散射信号，即镜面反射点附近的信号；非镜向模式则是指镜面反射点位于反射天线的有效覆盖区之外，此时接收到的信号主要是非镜向散射信号，即远离镜面反射点的信号。
由于非镜向散射系数随海面风向呈周期性波动，可以利用时延-多普勒(Delay-Doppler Mapping，DDM)均值（Delay-Doppler Mapping Average，DDMA）作为海面风矢量反演的特征观测量。DDMA是星下点天线波束中心附近特定时延-多普勒窗内散射系数的均值。在星下点非镜向观测配置下，星载GNSS-R接收机（又称为GNSS-R卫星或者GNSS-R接收卫星）在同一观测海域同时接收到三颗导航卫星的星下点非镜向散射信号的概率很大。基于此，旨在利用三颗方位角不同的导航卫星的星下点非镜向散射信号实现海面风矢量的反演。
本实施例提供了一种星下点观测的星载GNSS-R海面风场反演方法，特别是涉及基于极大似然估计（Maximum Likelihood Estimation，MLE）的星载GNSS反射信号海面风矢量反演算法，主要包括：分析非镜向海面散射信号对海面风矢量的敏感性，构建星下点非镜向独立观测模式，定义了该模式下对海面风矢量敏感的特征观测量。在此基础上，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演。
如图1所示，本实施例提供的一种星下点观测的星载GNSS-R海面风场反演方法，包括如下步骤。
步骤100：建立星载环境下的坐标系；如图2所示，所述坐标系的原点位于镜面反射点S0；所述坐标系的Z轴指向镜面反射点S0切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧。
导航卫星T、镜面反射点S0和GNSS-R接收卫星R的位置坐标可分别表示为：
（1）；
（2）；
（3）；
式中，为镜面反射点S0处导航信号的入射角/>或者散射角/>，即在镜面反射点处/>；Rt为镜面反射点S0与导航卫星T之间的距离，Rr为镜面反射点S0与GNSS-R接收卫星R之间的距离。
散射点处的散射向量q则可表示为：
（4）；
式中，、/>和/>分别为X轴、Y轴和Z轴的单位向量；qx表示散射向量q在X轴的单位向量，qy表示散射向量q在Y轴的单位向量，qz表示散射向量q在Z轴的单位向量，入射信号和散射信号的单位向量m和n可分别表示为：
（5）；
（6）；
式中，和/>分别为入射信号和散射信号相对于镜面反射点S0处的入射面的方位角；/>和/>分别为导航信号的入射角和散射角，两者的范围均为/>，满足：
（7）；
（8）。
其中，Sy表示散射点S向量在y方向的分量。
在本实施例所述的坐标系中，将海面风向定义为海面风矢量与所述坐标系的X轴正半轴侧的夹角；该夹角/>的范围为[0°，360°]。
步骤200：在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量。
在本实施例中，步骤200具体包括：
在所述坐标系下，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行处理，得到每颗导航卫星对应的DDMA测量值，进而确定每颗导航卫星对应的DDMA测量值的残差的条件概率密度函数；然后根据所述条件概率密度函数，确定似然函数，并计算所述似然函数的最大值，进而得到海面风矢量。
受噪声的影响，每一个DDMA测量值可表示为：
（9）；
式中，DDMAi为第i个导航卫星的DDMA测量值；为GNSS-R接收卫星在无噪声时测量的第i个导航卫星的DDMA真实值；/>和/>分别为海面风速和海面风向；/>为第i个导航卫星的入射角；/>为随机测量误差，是由各种随机噪声引起的，可认为其满足均值为0，方差为/>的高斯分布，即/>；i的取值范围为1至N。
在经验地球物理模式函数（Geophysical Model Function，GMF）中，DDMA是星载GNSS-R的已知特征观测量数据，入射角也可根据收发卫星位置坐标得到，海面风矢量，即海面风速和海面风向，为未知变量。如果对于每一个海面风矢量观测区域有两个独立的DDMA测量结果，则可以根据经验GMF构成两个方程，求解方程组，理论上可求得海面风矢量解。而经验GMF的不确定性导致了DDMA模型预测值与DDMA真实值之间也存在一个偏差，可称其为模型误差。DDMA模型预测值与DDMA真实值之间的关系可表示为：
（10）；
式中，为模型预测值；/>为模型误差，也可以看作是一个均值为0，方差为/>的高斯分布。
对于给定的海面风矢量，DDMAi测量值和模型预测值之间的残差Ri可定义为：
（11）。
由于测量随机误差和模型误差/>相互独立，所以残差Ri也是一个均值为0，方差为/>的高斯分布。方差/>可表示为：
（12）。
因此，每一个DDMA测量值的残差Ri的条件概率密度函数为：
（13）。
其中，W为海面风矢量；
假设在同一观测海域内星载GNSS-R同时接收了N颗导航卫星的星下点非镜向散射信号（N=3），并得到了相应的DDMA测量值，所有这些测量值对应于同一个海面风矢量。由于每个DDMA测量值是相互独立的，意味着对应的残差Ri也相互独立，所以这些残差的联合条件概率密度函数为：
（14）。
当测量值给定时，/>是参数/>的函数，记为：
（15）。
函数即为似然函数，求/>使似然函数取最大值，即满足下式：
（16）。
则对应的即为需要求解的海面风矢量。
其中，sup {}表示“上确界”，即最小上界。
步骤300：计算所述海面风矢量的数值解。
由于似然函数形式复杂，很难直接求得海面风矢量的数值解，本实施例采用数值搜索算法，计算所述海面风矢量的数值解，具体步骤如下：
（1）取海面风向为0°，以给定的起始风速为搜索起点，分别计算风速起始点的似然函数值和右边相邻点的似然函数值，比较两者大小，如果风速起始点的似然函数值小于右边相邻点的似然函数值，则继续向右搜索，反之向左搜索，直到找到使似然函数取最大值的点，并记录似然函数值、海面风速和海面风向。
（2）令海面风向值增加一个搜索间隔，风速起始点为上一个风向搜索到的海面风速，按照与第1步相同的搜索步骤寻找似然函数取最大值的点，并记录似然函数值、海面风速和海面风向。
（3）重复第2步的操作，将风向区间0°~360°搜索完，共得到361组似然函数值、海面风速和海面风向。
（4）根据第3步得到的结果在风向区间0°~360°寻找似然函数的局部极大值，并记录似然函数值、海面风速和海面风向。
（5）将第4步搜索得到的似然函数局部极大值从大到小排序，海面风矢量的模糊解为前四个局部极大值对应的海面风速和海面风向。其中，第一模糊解为排名第一的海面风矢量模糊解。
上述算法利用了似然函数在风矢量二维空间的分布特征，避免了在整个海面风速区间维度逐点搜索似然函数最大值，搜索效率高。
实施例二
为了执行上述实施例一对应的方法，以实现相应的功能和技术效果，下面提供一种星下点观测的星载GNSS-R海面风场反演系统。
本实施例提供的一种星下点观测的星载GNSS-R海面风场反演系统，包括：
坐标系构建模块，用于建立星载环境下的坐标系；所述坐标系的原点位于镜面反射点；所述坐标系的Z轴指向镜面反射点切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧。
海面风矢量计算模块，用于在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量。
数值解计算模块，用于计算所述海面风矢量的数值解。
进一步地，该系统还包括：海面风向定义模块，用于在所述坐标系中，将海面风向定义为海面风矢量与所述坐标系的X轴正半轴侧的夹角。
进一步地，该海面风矢量计算模块，具体包括：
条件概率密度函数确定模块，用于在所述坐标系下，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行处理，得到每颗导航卫星对应的DDMA测量值，并确定每颗导航卫星对应的DDMA测量值的残差的条件概率密度函数。
海面风矢量计算单元，用于根据所述条件概率密度函数，确定似然函数，并计算所述似然函数的最大值，进而得到海面风矢量。
进一步地，该数值解计算模块，具体包括：
数值解计算单元，用于采用数值搜索算法，计算所述海面风矢量的数值解。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (6)

1.一种星下点观测的星载GNSS-R海面风场反演方法，其特征在于，包括：
建立星载环境下的坐标系；所述坐标系的原点位于镜面反射点；所述坐标系的Z轴指向镜面反射点切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧；
在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量，具体包括：
在所述坐标系下，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行处理，得到每颗导航卫星对应的DDMA测量值，并确定每颗导航卫星对应的DDMA测量值的残差的条件概率密度函数；
根据所述条件概率密度函数，确定似然函数，并计算所述似然函数的最大值，进而得到海面风矢量；
计算所述海面风矢量的数值解。
2.根据权利要求1所述的一种星下点观测的星载GNSS-R海面风场反演方法，其特征在于，还包括：在所述坐标系中，将海面风向定义为海面风矢量与所述坐标系的X轴正半轴侧的夹角。
3.根据权利要求1所述的一种星下点观测的星载GNSS-R海面风场反演方法，其特征在于，计算所述海面风矢量的数值解，具体包括：
采用数值搜索算法，计算所述海面风矢量的数值解。
4.一种星下点观测的星载GNSS-R海面风场反演系统，其特征在于，包括：
坐标系构建模块，用于建立星载环境下的坐标系；所述坐标系的原点位于镜面反射点；所述坐标系的Z轴指向镜面反射点切平面的法向；其中，导航卫星与GNSS-R接收卫星位于所述坐标系的YOZ平面，且所述导航卫星与所述GNSS-R接收卫星均位于所述坐标系的Z轴正半轴侧；所述导航卫星位于所述坐标系的Y轴正半轴侧，所述GNSS-R接收卫星位于所述坐标系的Y轴负半轴侧；
海面风矢量计算模块，用于在所述坐标系下，采用基于极大似然估计的星载GNSS反射信号海面风矢量反演算法，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行海面风矢量的反演，得到海面风矢量，具体包括：
条件概率密度函数确定模块，用于在所述坐标系下，对三颗方位角不同的导航卫星的星下点非镜向散射信号进行处理，得到每颗导航卫星对应的DDMA测量值，并确定每颗导航卫星对应的DDMA测量值的残差的条件概率密度函数；
海面风矢量计算单元，用于根据所述条件概率密度函数，确定似然函数，并计算所述似然函数的最大值，进而得到海面风矢量；
数值解计算模块，用于计算所述海面风矢量的数值解。
5.根据权利要求4所述的一种星下点观测的星载GNSS-R海面风场反演系统，其特征在于，还包括：海面风向定义模块，用于在所述坐标系中，将海面风向定义为海面风矢量与所述坐标系的X轴正半轴侧的夹角。
6.根据权利要求4所述的一种星下点观测的星载GNSS-R海面风场反演系统，其特征在于，数值解计算模块，具体包括：
数值解计算单元，用于采用数值搜索算法，计算所述海面风矢量的数值解。
CN202310979097.4A 2023-08-03 2023-08-03 星下点观测的星载gnss-r海面风场反演方法及系统 Active CN116699666B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310979097.4A CN116699666B (zh) 2023-08-03 2023-08-03 星下点观测的星载gnss-r海面风场反演方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310979097.4A CN116699666B (zh) 2023-08-03 2023-08-03 星下点观测的星载gnss-r海面风场反演方法及系统

Publications (2)

Publication Number Publication Date
CN116699666A CN116699666A (zh) 2023-09-05
CN116699666B true CN116699666B (zh) 2023-10-31

Family

ID=87837852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310979097.4A Active CN116699666B (zh) 2023-08-03 2023-08-03 星下点观测的星载gnss-r海面风场反演方法及系统

Country Status (1)

Country Link
CN (1) CN116699666B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3081962A1 (fr) * 2015-04-15 2016-10-19 Centre National d'Etudes Spatiales Dispositif de poursuite intègre de signaux gnss
CN111239209A (zh) * 2020-02-17 2020-06-05 中国科学院上海天文台 机会信号反射遥感的全极化单次反射仿真方法及系统
CN111611678A (zh) * 2020-03-27 2020-09-01 国家卫星海洋应用中心 台风风向模糊解优化方法、装置、电子设备及存储介质
CN113050134A (zh) * 2021-03-19 2021-06-29 中国人民解放军92859部队 基于卫星导航信息的海面风场反演观测方法
CN113341389A (zh) * 2021-05-12 2021-09-03 中国科学院国家空间科学中心 一种基于最大似然估计的微波散射计风场反演方法及系统
CN114660552A (zh) * 2022-03-03 2022-06-24 北京卫星信息工程研究所 星载gnss-s雷达舰船目标信号接收与直射干扰抑制方法
CN114910661A (zh) * 2022-05-13 2022-08-16 北京大学 海面风速的反演方法、装置、介质和计算设备
CN115184916A (zh) * 2022-07-11 2022-10-14 北京大学 一种海面风速联合反演方法、装置、介质及计算设备
CN115754338A (zh) * 2022-11-09 2023-03-07 北京卫星信息工程研究所 多星协同gnss-r海面矢量风场探测系统及在轨反演方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9341719B2 (en) * 2012-09-07 2016-05-17 Qualcomm Technologies International, Ltd. Psuedo maximum likelihood tracking for global navigation satellite systems

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3081962A1 (fr) * 2015-04-15 2016-10-19 Centre National d'Etudes Spatiales Dispositif de poursuite intègre de signaux gnss
CN111239209A (zh) * 2020-02-17 2020-06-05 中国科学院上海天文台 机会信号反射遥感的全极化单次反射仿真方法及系统
CN111611678A (zh) * 2020-03-27 2020-09-01 国家卫星海洋应用中心 台风风向模糊解优化方法、装置、电子设备及存储介质
CN113050134A (zh) * 2021-03-19 2021-06-29 中国人民解放军92859部队 基于卫星导航信息的海面风场反演观测方法
CN113341389A (zh) * 2021-05-12 2021-09-03 中国科学院国家空间科学中心 一种基于最大似然估计的微波散射计风场反演方法及系统
CN114660552A (zh) * 2022-03-03 2022-06-24 北京卫星信息工程研究所 星载gnss-s雷达舰船目标信号接收与直射干扰抑制方法
CN114910661A (zh) * 2022-05-13 2022-08-16 北京大学 海面风速的反演方法、装置、介质和计算设备
CN115184916A (zh) * 2022-07-11 2022-10-14 北京大学 一种海面风速联合反演方法、装置、介质及计算设备
CN115754338A (zh) * 2022-11-09 2023-03-07 北京卫星信息工程研究所 多星协同gnss-r海面矢量风场探测系统及在轨反演方法

Also Published As

Publication number Publication date
CN116699666A (zh) 2023-09-05

Similar Documents

Publication Publication Date Title
CN110824510B (zh) 一种提高gnss-r测高卫星接收海面反射信号数量的方法
JP4875541B2 (ja) 方位検出方法、物体検出装置、プログラム
Matos et al. Development and implementation of a low-cost LBL navigation system for an AUV
Park et al. Radar localization and mapping for indoor disaster environments via multi-modal registration to prior LiDAR map
CN110988942B (zh) 一种星载gnss-r镜面反射点位置精确计算方法
Sun et al. Underwater acoustic localization of the black box based on generalized second-order time difference of arrival (GSTDOA)
CN107843895B (zh) 一种双多普勒雷达三维风场反演方法
CN112513679A (zh) 一种目标识别的方法和装置
Mohammadloo et al. Correcting multibeam echosounder bathymetric measurements for errors induced by inaccurate water column sound speeds
CN115047464B (zh) 双探头多波束测深方法和装置
Crouse Basic tracking using nonlinear 3D monostatic and bistatic measurements in refractive environments
CN115047503A (zh) 陆面星载gnss反射信号镜面反射点校准方法及系统
CN116699666B (zh) 星下点观测的星载gnss-r海面风场反演方法及系统
Inzartsev et al. AUV application for inspection of underwater communications
CN117848285A (zh) 海底地形测绘方法、设备、测绘系统及存储介质
CN113109843A (zh) 基于双接收机伪距双差的欺骗信号检测抑制方法及装置
Emokpae et al. Surface-reflection-based communication and localization in underwater sensor networks
CN116540273B (zh) Gnss-r镜面反射点初值的位置确定方法、装置及电子设备
CN112034455A (zh) 一种计算gnss-r海面反射信号ddm图像的方法
CN112731400A (zh) 一种海上船只目标矢量速度估计方法及系统
Hanzel et al. Range finder models for mobile robots
JP2001083232A (ja) 受動ターゲット位置決定装置
CN116228603B (zh) 一种挂车周围障碍物的报警系统及装置
CN117592381B (zh) 一种大气波导参数反演模型训练方法、装置、设备及介质
CN113296135B (zh) 形变监测方法、装置和接收机

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