一种基于环扫SAR的海面波浪谱反演方法及系统技术方案

本发明专利技术公开了一种基于环扫SAR的海面波浪谱反演方法及系统，该方法包括：获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；对子图像进行归一化处理，得到归一化后的子图像；对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；将子图像反演的波浪谱数据进行融合，得到最终波浪谱信息。该系统包括：数据划分模块、归一化模块、波浪谱反演模块和融合模块。通过使用本发明专利技术，可以有效地克服单视向SAR的方位向非线性导致的误差，实现对波浪谱的大范围高时间精度观测。本发明专利技术作为一种基于环扫SAR的海面波浪谱反演方法及系统，可广泛应用于卫星遥感数据处理领域。

【技术实现步骤摘要】
一种基于环扫SAR的海面波浪谱反演方法及系统
本专利技术涉及卫星遥感数据处理领域，尤其涉及一种基于环扫SAR的海面波浪谱反演方法及系统。
技术介绍
海浪是发生在海洋表面的一类典型的波动现象，是海洋动力学中最常见的一种形式。研究海浪运动对海洋工程作业、海洋开发、海洋捕捞与养殖、海洋灾害预测等活动具有重大意义。海浪的信息通常用波浪谱的形式表达，海浪的能量与频率之间具有特定的相关性，波浪谱表征了海浪的能量在频率和传播方向上的集中分布。波浪谱可以观测海域内海浪的波长、波高、周期和传播方向等要素，具有直观易读的特点。合成孔径雷达(SyntheticApertureRadar,SAR)是目前能够测量海洋表面波浪谱的星载设备之一，是用于海浪观测最有效的手段。但是，传统SAR的刈幅和覆盖范围有限，不适合开阔海域大覆盖面积的观测要求，不能满足大范围波浪波谱测量的高时间分辨率要求。
技术实现思路
为了解决上述技术问题，本专利技术提供了一种基于环扫SAR的海面波浪谱反演方法及系统，提高刈幅和覆盖范围，有效地克服单视向SAR的方位向非线性导致的误差。本专利技术所采用的第一技术方案是：一种基于环扫SAR的海面波浪谱反演方法，包括以下步骤：S1、获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；S2、对子图像进行归一化处理，得到归一化后的子图像；S3、对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；S4、将子图像反演的波浪谱数据进行融合，得到最终波浪谱信息。进一步，所述对子图像进行归一化处理具体为非线性归一化处理，公式如下：上式中，I0(x,t)代表原始子图像，x＝(x,y)代表图像域的二维坐标，x代表距离向，y为方位向，t代表时间，<I0>代表I0(x,t)的平均统计量，Is(x,t)代表归一化后的子图像，α1、α2表示对应参数。进一步，所述对归一化的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据这一步骤，其具体包括：S31、基于交叉谱正演模型对归一化的图像进行处理，得到观察交叉谱；S32、获取风速风向信息并生成初猜波浪谱；S33、通过初猜波浪谱进行正演计算得到初猜交叉谱，即仿真交叉谱；S34、基于预设的价格函数、观察交叉谱和初猜交叉谱，判断是否达到收敛条件；S35、判断是否达到收敛条件，计算迭代步长和方向导数并对波浪谱进行修正，得到修正后的波浪谱；S36、返回步骤S33，迭代直至达到收敛条件，得到对应子图像反演的波浪谱数据。进一步，所述预设的价格函数公式如下：J＝∫[|Pql(k',t)-Pobs(k',t)|2+μ|S(k')-Sg(k')|2]·Wp(k')dk'上式中，Wp为一个非负权重函数，Pobs(k',t)为观察交叉谱，Pql(k',t)为初猜交叉谱，k'为初始波数域，μ为一加权比例系数，S(k')为反演得到的波浪谱，Sg(k')为初猜谱，k为(kx,ky)二维波数域。进一步，所述修正后的波浪谱的表达式如下：Sn+1(k)＝Sn(k)+β1ΔS1(k)+β2ΔS2(k)上式中，Sn(k)为第n次迭代得到的波浪谱，Sn+1(k)为第n+1次迭代的波浪谱，ΔS1(k)和ΔS2(k)代表步长，β1和β2代表修正系数。进一步，所述将子图像反演的波浪谱数据进行融合，得到最终波浪谱信息这一步骤，其具体包括：S41、计算子图像对应的点响应函数、波浪谱线性函数和波浪谱；S42、将对应的点响应函数、波浪谱线性函数和波浪谱进行加权融合，得到多角度融合后的波浪谱。进一步，所述最终波浪谱信息S(k)的计算公式如下：上式中，Hn(k)、Tn(k)、Sn(k)分别为第n副子孔径图像的分辨率点响应函数、线性调制系数以及波浪谱。本专利技术所采用的第二技术方案是：一种基于环扫SAR的海面波浪谱反演系统，包括：数据划分模块，用于获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；归一化模块，用于对子图像进行归一化处理，得到归一化后的子图像；波浪谱反演模块，用于对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；融合模块，将子图像反演的波浪谱数据进行融合，得到最终波浪谱信息。本专利技术方法及系统的有益效果是：本专利技术通过交叉谱对波浪谱进行反演，克服了波浪谱波向180°模糊的问题；同时针对环扫SAR方位角360°旋转的特点，修正了方位向位移传递函数，最终将不同视向的子图像的反演结果进行融合，有效地克服单视向SAR的方位向非线性导致的误差，实现对波浪谱大范围高时间分辨率的观察。附图说明图1是本专利技术一种基于环扫SAR的海面波浪谱反演方法的总流程图；图2是本专利技术具体实施例环扫SAR的示意图；图3是本专利技术具体实施例子图像全局坐标系到局部坐标系的转换示意图；图4是本专利技术具体实施例波浪谱反演的示意图。图5是本专利技术一种基于环扫SAR的海面波浪谱反演系统的结构框图。具体实施方式下面结合附图和具体实施例对本专利技术做进一步的详细说明。对于以下实施例中的步骤编号，其仅为了便于阐述说明而设置，对步骤之间的顺序不做任何限定，实施例中的各步骤的执行顺序均可根据本领域技术人员的理解来进行适应性调整。环扫SAR是一种边飞行边环扫的新体制曲线雷达，雷达载体在以速度向前飞行的同时，天线以垂直地面为轴快速进行匀速扫描，在地面形成近似环状的成像区域，通过将多个环状区域进行叠加，形成超宽的观测刈幅，环扫SAR示意图如图2所示。参照图1，本专利技术提供了一种基于环扫SAR的海面波浪谱反演方法，该方法包括以下步骤：S1、获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；S2、对子图像进行归一化处理，得到归一化后的子图像；S3、对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；S4、将子图像反演的波浪谱数据进行融合，得到最终波浪谱。进一步作为本方法的优选实施例，所述获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像这一步骤，其具体包括：S11、获取环扫SAR数据；S12、按方位角将环扫SAR数据划分为多个子图像，并将子图像从全局坐标系转换到局部的斜距-方位坐标系。具体地，海浪到交叉谱的传递函数与视向角有很大的关系，而不同视向角的数据不能统一处理。因此需要将环扫SAR图像进行分块处理，由于每个子图像的视向角较小，子图像内的传递函数近似一致，所以可以对子图像进行统一处理。对于环扫SAR来说，子图像像素数和坐标轴的方向会随着方位角的变化而变化，为了对不同方位角进行统一处理，需将各个角度的反演结果插值到相同的分辨率，并旋转到统一的方位角。由于方位向近似为沿波束本文档来自技高网...

【技术保护点】
1.一种基于环扫SAR的海面波浪谱反演方法，其特征在于，包括以下步骤：/nS1、获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；/nS2、对子图像进行归一化处理，得到归一化后的子图像；/nS3、对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；/nS4、将子图像反演的波浪谱数据进行融合，得到最终波浪谱。/n

【技术特征摘要】
1.一种基于环扫SAR的海面波浪谱反演方法，其特征在于，包括以下步骤：
S1、获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像；
S2、对子图像进行归一化处理，得到归一化后的子图像；
S3、对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据；
S4、将子图像反演的波浪谱数据进行融合，得到最终波浪谱。


2.根据权利要求1所述一种基于环扫SAR的海面波浪谱反演方法，其特征在于，所述获取环扫SAR数据并对环扫SAR数据进行分块处理，得到多个子图像这一步骤，其具体包括：
S11、获取环扫SAR数据；
S12、按方位角将环扫SAR数据划分为多个子图像，并将子图像从全局坐标系转换到局部的斜距-方位坐标系。


3.根据权利要求2所述一种基于环扫SAR的海面波浪谱反演方法，其特征在于，所述对子图像进行归一化处理具体为非线性归一化处理，公式如下：



上式中，I0(x,t)代表原始子图像，x＝(x,y)代表图像域的二维坐标，x代表距离向，y为方位向，t代表时间，<I0>代表I0(x,t)的平均统计量，Is(x,t)代表归一化后的子图像，α1、α2表示对应参数。


4.根据权利要求3所述一种基于环扫SAR的海面波浪谱反演方法，其特征在于，所述对归一化后的子图像进行交叉谱计算和波浪谱反演，得到对应子图像反演的波浪谱数据这一步骤，其具体包括：
S31、基于交叉谱正演模型对归一化的图像进行处理，得到观察交叉谱；
S32、获取风速风向信息并生成初猜波浪谱；
S33、通过初猜波浪谱进行正演计算得到初猜交叉谱，即仿真交叉谱；
S34、基于预设的价格函数、观察交叉谱和初猜交叉谱，判断是否达到收敛条件；
S35、判断是否达到收敛条件，计算迭代步长和方向导数并对波浪谱进行修正，得到修正后的波浪谱；
S36、返回步骤S33，迭代直至达到收敛条件，得到对应子图像反演的波浪谱数据。


5.根据权利要求4所...

【专利技术属性】
技术研发人员：王小青，姚晓楠，黄海风，
申请(专利权)人：中山大学，
类型：发明
国别省市：广东;44