CN113050091A - 一种星载合成孔径雷达风速风向联合反演方法及系统 - Google Patents

一种星载合成孔径雷达风速风向联合反演方法及系统 Download PDF

Info

Publication number
CN113050091A
CN113050091A CN202110250397.XA CN202110250397A CN113050091A CN 113050091 A CN113050091 A CN 113050091A CN 202110250397 A CN202110250397 A CN 202110250397A CN 113050091 A CN113050091 A CN 113050091A
Authority
CN
China
Prior art keywords
wind speed
sar
sea surface
surface wind
prediction
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
Application number
CN202110250397.XA
Other languages
English (en)
Other versions
CN113050091B (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.)
National Ocean Technology Center
Original Assignee
National Ocean Technology Center
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 National Ocean Technology Center filed Critical National Ocean Technology Center
Priority to CN202110250397.XA priority Critical patent/CN113050091B/zh
Publication of CN113050091A publication Critical patent/CN113050091A/zh
Application granted granted Critical
Publication of CN113050091B publication Critical patent/CN113050091B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/956—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on ship or other platform
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • 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)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Ocean & Marine Engineering (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种星载合成孔径雷达风速风向联合反演方法及系统。该方法包括：获取风速预测模型及其预测误差，其中，风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向和距离向海面风速分量；将待预测交叉谱特征输入风速预测模型；基于GMF查找表，根据雷达入射角、后向散射系数、风速预测模型的初步预测结果以及风速预测模型的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果，并基于该最终预测结果计算海面风速和海面风向。本发明能够由同极化SAR数据同时准确的反演风速、风向，且不依赖外部数据。

Description

一种星载合成孔径雷达风速风向联合反演方法及系统
技术领域
本发明涉及海面风观测领域，特别是涉及一种星载合成孔径雷达风速风向联合反演方法及系统。
背景技术
海面风(包括风速、风向两个要素)是海气相互作用研究中最重要的参数之一，是海洋波浪、大洋环流的主要驱动力，对于海气之间的物质、热量的交换具有重要的作用。海面风的准确观测，对海洋气象灾害预报、海洋航行保障、海洋工程建设等均有着重要意义。与传统的浮标、船舶等海上单点测量方式相比，合成孔径雷达(SyntheticAperture Radar，SAR)遥感卫星以其全天时、全天候观测，高分辨率、大幅宽覆盖等优势，已经成为目前海面风的重要卫星遥感观测手段。
由于现有星载SAR系统只配置一根固定的雷达天线，无法如散射计一样利用对海面的多次观测的方式同时获取风速、风向两个要素的信息，因此很多技术方案仅可通过SAR图像反演海面风速的空间分布(即风场)，如申请号为202010621954X的发明专利公开的就是由SAR灰度图像仅提取海面风速。而且，现有SAR遥感资料同时反演风速、风向的技术，主要存在如下问题：
一、综合利用SAR数据的同极化(如VV极化)和交叉极化(如VH极化)的信息，例如专利号为201610201642.7的发明专利公开的基于交叉极化模型与CMOD5N的海面风速反演方法。但这类方法并不适用于仅有同极化的SAR图像。
二、先获取风向信息，再结合地球物理模式函数估计风速。这类技术又分为两种：
(1)当SAR图像包含风条纹等图像特征时，由SAR图像本身获取海面风向的信息(如Zhou,L.；et al.An Improved Local Gradient Method for Sea Surface WindDirection Retrieval from SAR Imagery.Remote Sens.2017.)，但是，大量不含有风条纹的SAR图像则无法使用该技术。
(2)依赖外部数值模式资料作为先验输入，如欧洲空间局业务化的Sentinel-1SAR数据OWI海面风反演算法，即是以欧洲中期预报中心(ECMWF)提供的预报风向作为辅助输入。但这种方法对于难以获得外部数据的情况并不适用，如船载SAR地面应用系统或星上快速处理系统。
因此，现有技术依赖外部数据，且难以实现由同极化SAR数据同时反演风速、风向。
发明内容
本发明的目的是提供一种星载合成孔径雷达风速风向联合反演方法及系统，能够由同极化SAR数据同时反演风速、风向，且不依赖外部数据。
为实现上述目的，本发明提供了如下方案：
一种星载合成孔径雷达风速风向联合反演方法，包括：
获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征；
将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量；
基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果；所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表；
基于SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。
可选的，根据代价函数
Figure BDA0002965796300000031
确定SAR方位向海面风速分量u以及SAR距离向海面风速分量v的代价，其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。
可选的，基于同极化SAR数据计算SAR交叉谱特征的方法包括：
获取同极化SAR数据的单视复图像；
对所述单视复图像进行交叉谱估计，获得SAR交叉谱特征，所述SAR交叉谱特征包括SAR交叉谱的实部、虚部和相位。
可选的，所述方法还包括：
训练风速预测模型：
获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
可选的，所述方法还包括：
确定风速预测模型的预测误差：
获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
将所述验证集中的样本输入所述风速预测模型；
根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
本发明还提供了一种星载合成孔径雷达风速风向联合反演系统，包括：
模型及误差获取模块，用于获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
待预测交叉谱特征求取模块，用于获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征；
初级预测模块，用于将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量；
终级寻优预测模块，用于基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果，所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表；
海面风速风向计算模块，用于基于所述SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。
可选的，所述终级寻优预测模块包括：代价计算单元，用于根据代价函数
Figure BDA0002965796300000041
计算SAR方位向海面风速分量u以及SAR距离向海面风速分量v的代价，其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。
可选的，所述待预测交叉谱特征求取模块包括：
单视复图像获取单元，用于获取同极化SAR数据的单视复图像；
SAR交叉谱特征求取单元，用于对所述单视复图像进行交叉谱估计，获得SAR交叉谱特征，所述SAR交叉谱特征包括SAR交叉谱的实部、虚部和相位。
可选的，所述系统还包括：
风速预测模型训练模块，用于获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
可选的，所述系统还包括：
风速预测模型的预测误差确定模块，用于获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；将所述验证集中的样本输入所述风速预测模型；根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供的星载合成孔径雷达风速风向联合反演方法及系统，由于充分利用了SAR卫星数据的交叉谱信息和后向散射信息，并且采用了深度卷积神经网络和地球物理模式函数查找表搜索技术，从而摆脱了对数值模式风场资料的依赖，而且能够由同极化SAR数据简便快捷地联合反演海面风速、风向，且反演精度满足业务需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例1提供的星载合成孔径雷达风速风向联合反演方法的流程示意图；
图2为本发明实施例1中深度学习网络的架构图；
图3为基于本发明的SAR反演方法得到的风速与同步的浮标观测数据的散点比较图；
图4为基于本发明的SAR反演方法得到的风向与同步的浮标观测数据的散点比较图；
图5为本发明实施例2提供的星载合成孔径雷达风速风向联合反演系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种星载合成孔径雷达风速风向联合反演方法及系统，能够由同极化SAR数据同时反演风速、风向，且不依赖外部数据。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
参见图1，本实施例提供了一种星载合成孔径雷达风速风向联合反演方法，该方法包括以下步骤：
步骤101：获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量。
步骤102：获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征。
具体过程可以如下：
读取待预测的同极化(如VV极化或HH极化)SAR数据的单视复图像，对SAR单视复图像进行交叉谱估计，获得待预测的SAR交叉谱特征：
1)根据SAR单视复图像，计算Doppler谱，并进行分视处理得到3个SAR子视图像，再对每个子视图像计算其海浪相对调制mi：
Figure BDA0002965796300000071
其中，Ii为每个子视图像的图像强度。
2)对于每两个相邻子视图像，计算交叉周期图：
Figure BDA0002965796300000072
其中，
Figure BDA0002965796300000073
为
Figure BDA0002965796300000074
的FFT变换。
3)计算SAR图像交叉谱的实部R、虚部I和相位P：
Figure BDA0002965796300000075
Figure BDA0002965796300000076
P＝arctan(I/R) (5)
其中，Re,Im表示取实部和取虚部运算。
4)将交叉谱的实部、虚部和相位分别作为RGB图像中的R通道、G通道和B通道合成RGB图像，作为待预测的SAR交叉谱特征。
步骤103：将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量。
步骤104：基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果，所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表。
具体的，GMF查找表是一个L*M*N大小的三维查找表σ0 GMF(θ,u,v)，其中L，M，N为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量三个参数索引的长度，σ0为以雷达入射角θ(采样间隔为0.1°，范围为20°-50°)、SAR方位向海面风速分量u(采样间隔为0.1m/s，范围0.1m/s-35m/s)、SAR距离向海面风速分量v(采样间隔为0.1m/s，范围0.1m/s-35m/s)作为输入，通过地球物理模式函数GMF(Geophysical Model Function)逐点计算得到的后向散射系数。
对于C波段VV极化的SAR数据，例如Sentinel-1，可以使用CMOD5.N地球物理模式函数进行查找表的计算：
Figure BDA0002965796300000081
其中，
Figure BDA0002965796300000082
为雷达的后向散射系数，θ为雷达入射角，
Figure BDA0002965796300000083
为海面风速，ψ＝arctan(Uaz-Ura)。B0，B1和B2均为风速和雷达入射角的函数。由于CMOD5.N地球物理模式函数为现有公开的技术，其细节为本领域技术人员熟知，不再赘述。
步骤104耦合CNN和GMF进行SAR风速、风向的反演，实施方式可以如下：
基于步骤103中得到的SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量，通过GMF查找表反演风分量：
1)找到查找表中雷达入射角索引值中最接近SAR数据雷达入射角的索引值i，并将三维查找表σ0 GMF(θ,u,v)转换为二维查找表σ0 GMF(i,u,v)。
2)基于海面风速分量u,v索引，逐点计算如下代价函数：
Figure BDA0002965796300000084
其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。对于VV极化的Sentinel-1SAR数据，Δσ0优选为0.5dB。
上述雷达后向散射系数σ0 SAR的获取方式可以如下：
以Sentinel-1为例，辐射定标可按照欧洲空间局提供的用户手册中的定标公式完成。
Figure BDA0002965796300000085
其中，DN是Sentinel-1卫星SAR数据原始幅度数码值(DigitalNumber)，η为噪声等效散射系数值，A为后向散射系数定标因子。
3)遍历查找表风分量海面风速分量u,v，找出使代价函数J(u,v)最小的SAR方位向海面风速分量
Figure BDA0002965796300000091
SAR距离向海面风速分量
Figure BDA0002965796300000092
作为最终预测结果。
步骤105：基于所述SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。具体可采用以下公式进行计算：
Figure BDA0002965796300000093
Figure BDA0002965796300000094
其中，α为SAR卫星方位角,Usar为SAR反演得到的海面风速，ψN为SAR反演得到的相对正北方向的海面风向。
作为本实施例的一种实施方式，该方法还可以包括：
风速预测模型的训练步骤：获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
具体过程可以如下：
(1)：卫星数据的读取
1)读取同极化(如VV极化或HH极化)SAR数据的单视复图像(single lookcomplex,SLC)、成像时间和经纬度、雷达入射角、卫星方位角。具体的，可以从Sentinel-1卫星VV极化SAR数据中读取单视复图像、成像时间和经纬度、雷达入射角、卫星方位角。
2)同步的散射计数据的读取
根据1)得到的SAR图像的成像时间和经纬度，读取与SAR图像时间和空间相同步的散射计风速、风向数据。这里，SAR与散射计满足时间差不超过30分钟，距离差不超过25Km的条件，即可视为同步。
将散射计风速、风向转换为SAR方位向和距离向上的风分量。计算公式如下：
Uaz＝U*cos(φ-α) (11)
Ura＝U*sin(φ-α) (12)
其中，Uaz，Ura为SAR方位向和距离向上的风分量，U，φ分别为风速、风向，α为SAR卫星方位角。
(2)：SAR数据的预处理
对SAR单视复图像进行交叉谱估计，获得SAR交叉谱特征。
1)根据SAR单视复图像，计算Doppler谱，并进行分视处理得到3个SAR子视图像，再对每个子视图像计算其海浪相对调制mi：
Figure BDA0002965796300000101
其中，Ii为每个子视图像的图像强度。
2)对于每两个相邻子视图像，计算交叉周期图：
Figure BDA0002965796300000102
其中，
Figure BDA0002965796300000103
为
Figure BDA0002965796300000104
的FFT变换。
3)计算SAR图像交叉谱的实部R、虚部I和相位P：
Figure BDA0002965796300000105
Figure BDA0002965796300000106
P＝arctan(I/R)
其中，Re,Im表示取实部和取虚部运算。
4)将交叉谱的实部、虚部和相位分别作为RGB图像中的R通道、G通道和B通道合成RGB图像，作为卷积神经网络输入的SAR交叉谱特征。
(3)：构建估计初猜风分量值的卷积神经网络CNN(Convolution Neural Network)
1)构建卷积深度学习模型所需的样本数据集
将(2)中获取的SAR交叉谱特征作为样本；将从同步的散射计获取的SAR方位向、距离向风分量数据作为该样本的标签。然后，将样本数据集随机打乱顺序，划分为训练数据集和验证数据集。可选的，选择随机打乱顺序的前70％样本数据集作为训练数据集，其余作为验证数据集。
2)搭建卷积神经网络架构
本发明采用AlexNet深度学习网络架构。如图2所示，共有八层，其中前五层为卷积层，后三层为全连接层。其中，前两个全连接层使用了30％比率的随机失活(Dropout)正则化方法，以减少过拟合；每一个卷积层和全连接层都使用Relu激活函数。神经网络的输出有两个，分别是SAR方位向和距离向上的风速分量。
3)训练卷积神经网络，保存最优模型结构和最优参数集。
将训练数据集输入到卷积神经网络，并利用计算机GPU进行训练，直至收敛。模型训练完成后，保存最优模型结构和最优参数集。
可选的，训练时采用随机梯度下降(Stochastic gradient descent)优化器进行模型优化，共进行120轮训练，每一批训练的数据量大小为128；初始学习率设为0.01，每25轮训练后学习率下降至原来的1/10。
作为本实施例的一种实施方式，该方法还可以包括：
风速预测模型的预测误差确定步骤：获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；将所述验证集中的样本输入所述风速预测模型；根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
具体过程可以如下：
首先，将验证数据集的SAR交叉谱特征输入到训练好的卷积神经网络，采用保存的最优模型结构和最优参数集，计算方位向和距离向上的风分量的估计值
Figure BDA0002965796300000111
和
Figure BDA0002965796300000112
然后，利用验证数据集的样本标签，计算风分量初猜值的误差ΔUaz和ΔUra：
Figure BDA0002965796300000113
Figure BDA0002965796300000114
其中，std表示计算标准差运算，Uaz，Uaz为从散射计得到的SAR方位向和距离向上的风分量。
与现有技术相比，本发明的具有以下优势：
(1)基于同极化SAR数据对海面风速、风向进行反演，可以摆脱对外部风场资料的依赖，且其精度与现有技术相当。
如图3，图4所示，本发明应用于VV极化的Sentinel-1SAR数据，与浮标实际观测相比，风速反演均方根误差(RMSE)为1.9m/s，风向反演均方根误差(RMSE)为19.7°，均满足业务化要求的2.0m/s和20°。本发明利用卷积神经网络深度挖掘了SAR交叉谱的信息，并结合了地球物理模式函数和SAR后向散射系数，充分挖掘耦合SAR卫星数据的交叉谱信息和后向散射信息，摆脱了对外部风场资料的依赖。
(2)反演耗时少，实现简便。由于本申请预先训练了卷积神经网络模型，在实际反演中可以直接使用已保存最优模型结构和最优参数集；本发明预先生成了地球物理模式函数查找表，在实际反演中可以直接采用对查找表的快速索引，避免了大量的重复计算。
实施例2
参见图5，本实施例提供了一种星载合成孔径雷达风速风向联合反演系统，该系统包括：
模型及误差获取模块501，用于获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
待预测交叉谱特征求取模块502，用于获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征；
初级预测模块503，用于将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量；
终级寻优预测模块504，用于基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果，所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表；
海面风速风向计算模块505，用于基于所述SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。
其中，所述终级寻优预测模块504包括：代价计算单元，用于根据代价函数
Figure BDA0002965796300000131
计算SAR方位向海面风速分量u以及SAR距离向海面风速分量v的代价，其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。
所述待预测交叉谱特征求取模块502包括：
单视复图像获取单元，用于获取同极化SAR数据的单视复图像；
SAR交叉谱特征求取单元，用于对所述单视复图像进行交叉谱估计，获得SAR交叉谱特征，所述SAR交叉谱特征包括SAR交叉谱的实部、虚部和相位。
作为实施例的一种实施方式，所述系统还包括：
风速预测模型训练模块，用于获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
作为实施例的一种实施方式，所述系统还包括：
风速预测模型的预测误差确定模块，用于获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；将所述验证集中的样本输入所述风速预测模型；根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种星载合成孔径雷达风速风向联合反演方法，其特征在于，包括：
获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征；
将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量；
基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果；所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表；
基于SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。
2.根据权利要求1所述的星载合成孔径雷达风速风向联合反演方法，其特征在于，根据代价函数
Figure FDA0002965796290000011
确定SAR方位向海面风速分量u以及SAR距离向海面风速分量v的代价，其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。
3.根据权利要求1所述的星载合成孔径雷达风速风向联合反演方法，其特征在于，基于同极化SAR数据计算SAR交叉谱特征的方法包括：
获取同极化SAR数据的单视复图像；
对所述单视复图像进行交叉谱估计，获得SAR交叉谱特征，所述SAR交叉谱特征包括SAR交叉谱的实部、虚部和相位。
4.根据权利要求1所述的星载合成孔径雷达风速风向联合反演方法，其特征在于，所述方法还包括：
训练风速预测模型：
获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
5.根据权利要求1所述的星载合成孔径雷达风速风向联合反演方法，其特征在于，所述方法还包括：
确定风速预测模型的预测误差：
获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
将所述验证集中的样本输入所述风速预测模型；
根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
6.一种星载合成孔径雷达风速风向联合反演系统，其特征在于，包括：
模型及误差获取模块，用于获取风速预测模型以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，其中，用于训练所述风速预测模型的训练样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；
待预测交叉谱特征求取模块，用于获取待预测的同极化SAR数据，并计算待预测同极化SAR数据的待预测交叉谱特征；
初级预测模块，用于将所述待预测交叉谱特征输入所述风速预测模型，得到SAR方位向海面风速初步预测分量与SAR距离向海面风速初步预测分量；
终级寻优预测模块，用于基于GMF查找表，根据雷达入射角、后向散射系数、所述SAR方位向海面风速初步预测分量、所述SAR距离向海面风速初步预测分量以及所述风速预测模型对SAR方位向海面风速分量和SAR距离向海面风速分量的预测误差，选择确定代价最小的SAR方位向海面风速分量以及SAR距离向海面风速分量作为最终预测结果，所述GMF查找表为雷达入射角、SAR方位向海面风速分量、SAR距离向海面风速分量以及基于地球物理模式函数计算得到的散射系数之间的关系表；
海面风速风向计算模块，用于基于所述SAR方位向海面风速分量以及SAR距离向海面风速分量的最终预测结果计算海面风速和海面风向。
7.根据权利要求6所述的星载合成孔径雷达风速风向联合反演系统，其特征在于，所述终级寻优预测模块包括：代价计算单元，用于根据代价函数
Figure FDA0002965796290000031
计算SAR方位向海面风速分量u以及SAR距离向海面风速分量v的代价，其中，σ0 SAR为雷达后向散射系数，Δσ0为所述后向散射系数的不确定度，UazCNN为所述SAR方位向海面风速初步预测分量，UraCNN为所述SAR距离向海面风速初步预测分量，ΔUaz为所述风速预测模型对SAR方位向海面风速分量的预测误差，ΔUra为所述风速预测模型对所述SAR距离向海面风速分量的预测误差。
8.根据权利要求6所述的星载合成孔径雷达风速风向联合反演系统，其特征在于，所述待预测交叉谱特征求取模块包括：
单视复图像获取单元，用于获取同极化SAR数据的单视复图像；
SAR交叉谱特征求取单元，用于对所述单视复图像进行交叉谱估计，获得SAR交叉谱特征，所述SAR交叉谱特征包括SAR交叉谱的实部、虚部和相位。
9.根据权利要求6所述的星载合成孔径雷达风速风向联合反演系统，其特征在于，所述系统还包括：
风速预测模型训练模块，用于获取训练集，所述训练集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；采用所述训练集对卷积深度学习模型进行训练，得到所述风速预测模型。
10.根据权利要求6所述的星载合成孔径雷达风速风向联合反演系统，其特征在于，所述系统还包括：
风速预测模型的预测误差确定模块，用于获取验证集，所述验证集中的样本为基于同极化SAR数据计算得到的SAR交叉谱特征，标签为基于与所述SAR数据时间和空间相同步的散射计测量数据计算得到的SAR方位向海面风速分量与SAR距离向海面风速分量；将所述验证集中的样本输入所述风速预测模型；根据所述风速预测模型与所述样本的标签，确定所述风速预测模型的预测误差。
CN202110250397.XA 2021-03-08 2021-03-08 一种星载合成孔径雷达风速风向联合反演方法及系统 Active CN113050091B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110250397.XA CN113050091B (zh) 2021-03-08 2021-03-08 一种星载合成孔径雷达风速风向联合反演方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110250397.XA CN113050091B (zh) 2021-03-08 2021-03-08 一种星载合成孔径雷达风速风向联合反演方法及系统

Publications (2)

Publication Number Publication Date
CN113050091A true CN113050091A (zh) 2021-06-29
CN113050091B CN113050091B (zh) 2021-09-28

Family

ID=76510520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110250397.XA Active CN113050091B (zh) 2021-03-08 2021-03-08 一种星载合成孔径雷达风速风向联合反演方法及系统

Country Status (1)

Country Link
CN (1) CN113050091B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113534159A (zh) * 2021-07-08 2021-10-22 国家海洋技术中心 一种基于sar卫星遥感数据的海面风速反演方法及系统
CN114167076A (zh) * 2021-11-10 2022-03-11 中国地质大学（武汉） 海面风速风向反演模型的建立方法及风速风向的反演方法
CN114324973A (zh) * 2022-03-17 2022-04-12 南方海洋科学与工程广东省实验室(广州) 台风风速反演方法、装置、电子设备及存储介质
CN117969881A (zh) * 2024-03-29 2024-05-03 中国石油大学(华东) 基于遥感数据反演台风海况海面风速的方法、装置及介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130147659A1 (en) * 2011-10-24 2013-06-13 Korea Meteorological Administration System for detecting sea-surface wind, using satellite observation, and a method for detecting sea-surface wind
EP2610636A1 (en) * 2011-12-29 2013-07-03 Windward Ltd. Providing near real-time maritime insight from satellite imagery and extrinsic data
CN104698460A (zh) * 2015-02-26 2015-06-10 中国人民解放军理工大学 一种双频共面合成孔径雷达海面风场反演方法
CN105259369A (zh) * 2015-10-27 2016-01-20 国家海洋技术中心 风速风向测量仪
CN105891832A (zh) * 2016-03-31 2016-08-24 南京信息工程大学 一种基于交叉极化模型与cmod5n的海面风速反演方法
CN106501804A (zh) * 2016-11-25 2017-03-15 中国石油大学(华东) 一种利用全极化sar回波数据解析海面风浪谱的方法
CN106610491A (zh) * 2016-12-21 2017-05-03 广州市气象台 星载sar后向散射系数的检验方法及装置
CN107505616A (zh) * 2017-09-15 2017-12-22 浙江大学 一种基于sar的海面风场反演最优分辨率的判定方法
CN109283495A (zh) * 2018-09-06 2019-01-29 南京信息工程大学 一种基于交叉谱分析的x波段雷达海流反演方法
CN111580138A (zh) * 2020-06-18 2020-08-25 北京卫星信息工程研究所 多波束gnss-r系统、海面风场反演方法和预测方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130147659A1 (en) * 2011-10-24 2013-06-13 Korea Meteorological Administration System for detecting sea-surface wind, using satellite observation, and a method for detecting sea-surface wind
EP2610636A1 (en) * 2011-12-29 2013-07-03 Windward Ltd. Providing near real-time maritime insight from satellite imagery and extrinsic data
CN104698460A (zh) * 2015-02-26 2015-06-10 中国人民解放军理工大学 一种双频共面合成孔径雷达海面风场反演方法
CN105259369A (zh) * 2015-10-27 2016-01-20 国家海洋技术中心 风速风向测量仪
CN105891832A (zh) * 2016-03-31 2016-08-24 南京信息工程大学 一种基于交叉极化模型与cmod5n的海面风速反演方法
CN106501804A (zh) * 2016-11-25 2017-03-15 中国石油大学(华东) 一种利用全极化sar回波数据解析海面风浪谱的方法
CN106610491A (zh) * 2016-12-21 2017-05-03 广州市气象台 星载sar后向散射系数的检验方法及装置
CN107505616A (zh) * 2017-09-15 2017-12-22 浙江大学 一种基于sar的海面风场反演最优分辨率的判定方法
CN109283495A (zh) * 2018-09-06 2019-01-29 南京信息工程大学 一种基于交叉谱分析的x波段雷达海流反演方法
CN111580138A (zh) * 2020-06-18 2020-08-25 北京卫星信息工程研究所 多波束gnss-r系统、海面风场反演方法和预测方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
MONALDO, FM等: "A systematic comparison of QuikSCAT and SAR ocean surface wind speeds", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
WANG, LH等: "Deriving Sea Surface Wind from Synthetic Aperture Radar Based on Fourier Transform and Neural Network", 《2017 10TH INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING, BIOMEDICAL ENGINEERING AND INFORMATICS 》 *
XIE, TAO等: "A Hybrid Method for Ocean Surface Wind Retrieval from Polarization SAR Image", 《2015 INTERNATIONAL CONFERENCE ON ENVIRONMENT, MANUFACTURING INDUSTRY AND ECONOMIC DEVELOPMENT》 *
戈书睿等: "基于同步风场的机载SAR海浪参数反演方法", 《海洋科学》 *
林莉: "基于SAR复数据的海浪谱反演方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
艾未华等: "C波段机载合成孔径雷达海面风场反演新方法", 《物理学报》 *
赵现斌等: "机载C波段全极化SAR海面风矢量反演理论研究及实验验证", 《物理学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113534159A (zh) * 2021-07-08 2021-10-22 国家海洋技术中心 一种基于sar卫星遥感数据的海面风速反演方法及系统
CN113534159B (zh) * 2021-07-08 2022-04-26 国家海洋技术中心 一种基于sar卫星遥感数据的海面风速反演方法及系统
CN114167076A (zh) * 2021-11-10 2022-03-11 中国地质大学（武汉） 海面风速风向反演模型的建立方法及风速风向的反演方法
CN114324973A (zh) * 2022-03-17 2022-04-12 南方海洋科学与工程广东省实验室(广州) 台风风速反演方法、装置、电子设备及存储介质
CN117969881A (zh) * 2024-03-29 2024-05-03 中国石油大学(华东) 基于遥感数据反演台风海况海面风速的方法、装置及介质

Also Published As

Publication number Publication date
CN113050091B (zh) 2021-09-28

Similar Documents

Publication Publication Date Title
CN113050091B (zh) 一种星载合成孔径雷达风速风向联合反演方法及系统
Antoine et al. Assessment of uncertainty in the ocean reflectance determined by three satellite ocean color sensors (MERIS, SeaWiFS and MODIS‐A) at an offshore site in the Mediterranean Sea (BOUSSOLE project)
Dankert et al. Ocean wind fields retrieved from radar‐image sequences
CN113534159B (zh) 一种基于sar卫星遥感数据的海面风速反演方法及系统
CN112014842B (zh) 高分三号全极化sar数据的有效波高估计方法及系统
CN111781146B (zh) 利用高分辨率卫星光学影像的波浪参数反演方法
CN103698750A (zh) 海洋二号卫星散射计海面风场反演方法和装置
de Rosnay et al. Coupled data assimilation at ECMWF: Current status, challenges and future developments
Warren et al. Estimation of ocean surface currents from maximum cross correlation applied to GOCI geostationary satellite remote sensing data over the Tsushima (K orea) S traits
Wang et al. Technology of gravity aided inertial navigation system and its trial in South China Sea
CN105891832A (zh) 一种基于交叉极化模型与cmod5n的海面风速反演方法
Wang et al. The wide swath significant wave height: An innovative reconstruction of significant wave heights from CFOSAT’s SWIM and scatterometer using deep learning
Romeiser et al. Advanced remote sensing of internal waves by spaceborne along-track InSAR—A demonstration with TerraSAR-X
CN113591387A (zh) 基于Huber范数约束的卫星资料反演降水方法及系统
Wang et al. Sea ice classification with convolutional neural networks using Sentinel-L scansar images
Lund et al. Marine X‐Band Radar Currents and Bathymetry: An Argument for a Wave Number‐Dependent Retrieval Method
Wang et al. Acquisition of the significant wave height from CFOSAT SWIM spectra through a deep neural network and its impact on wave model assimilation
Miranda et al. A new unconstrained approach to GNSS atmospheric water vapor tomography
CN114325697A (zh) 一种考虑入射角的宽幅交叉极化sar影像热噪声提取方法
Loria et al. Towards wind vector and wave height retrievals over inland waters using CYGNSS
CN117436329A (zh) 一种基于机器学习的合成孔径雷达图像风向反演方法
CN112162282A (zh) 一种基于合成孔径雷达海表流速反演方法
CN106872979A (zh) 基于海面起伏目标sar图像再聚焦的海浪参数获取方法
Ferreira et al. Under what conditions sar along-track interferometry is suitable for assessment of tidal energy resource
Collett et al. A new method for ocean wind direction retrieval from delay-doppler maps using stare processing and machine learning: Preliminary simulation results

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