基于子孔径图像的sar海洋图像中海浪杂波抑制方法
【技术领域】
[0001] 本发明涉及微波技术领域，尤其涉及一种基于子孔径图像的SAR海洋图像中海浪 杂波抑制方法。
【背景技术】
[0002] 合成孔径雷达（Synthetic Aperture Radar，简称SAR)是以微波谱段的电磁波作 为探测载体，借助合成孔径实现对观测对象的高精度二维成像的信息获取技术。与传统光 学成像相比，SAR具有远距探测、全天时、全天候的工作优点。海洋观测是SAR的一个重要 应用方向，SAR可以发现海面上波浪、风场、内波、溢油、漩涡等多种海洋现象，而且可以对一 些海洋现象进行参数反演。波浪是海面SAR图像中非常常见的现象，经常表现为密集的波 纹状纹理结构，在反演海浪谱等应用中SAR图像中的波浪条纹是非常有用的信息，但如果 希望从SAR图像中提取内波、漩涡等其他海洋现象信息时波浪纹理就会对有用的信息产生 干扰。传统的海浪杂波的抑制方法一般是利用波浪纹理与其他海洋现象在空间分布上的波 长和方向等几何特征差异，采用傅里叶分析、小波分析等方法，来进行波浪纹理的滤波。这 些方法虽然在一定程度上可以抑制海浪杂波的干扰，但是由于其他海洋现象与波浪虽然波 长和方向等分布有一定差异，但是很多时候也在频域或者小波域有显著的重叠区域，在抑 制海浪杂波的同时也会对其他海洋现象的纹理有一定损失。
[0003] 海浪除了在空间波长、方向等几何特征与其他海洋现象有差异外，另外一个重要 的差异是动态特性的差异。海浪有特定的频散规律，也就是特定波长的波浪，其振动频率也 是确定的，而其他海洋现象一般都不符合海浪的频散关系，动态特性一般要显著比波浪更 慢。SAR成像的特点是需要有一定合成孔径时间，因此SAR的回波信号中除了包含地物的几 何特征信息以外也包含了地物的动态特征信息，而常规的SAR图像处理中往往忽视了海洋 图像的动态特性，只利用了几何特性，因此对海浪杂波的抑制的同时对其他海洋现象纹理 也有损失。
【发明内容】
[0004] (一）要解决的技术问题
[0005] 鉴于上述技术问题，本发明提供了一种基于子孔径图像进行海浪杂波抑制的方 法。
[0006] (二）技术方案
[0007] 根据本发明的一个方面，提供了一种基于子孔径图像的SAR海洋图像中海浪杂波 抑制方法。该SAR海洋图像中海浪杂波抑制方法包括：步骤A :对SAR原始回波s进行方位 向傅里叶变换，得到斜距多普勒域数据S ;步骤B :对斜距多普勒域数据S通过方位向多普 勒频谱子孔径分解，得到两幅子孔径回波数据SJP S 2;步骤C :对两幅子孔径回波数据S JP S2，分别进行SAR成像、斜地转换和方位向/地距向的空间二维傅里叶变换，得到两幅子孔 径SAR图像的空间二维傅里叶分解形式S + (kx，ky)，S (kx，ky)，其中kx，ky分别为方位向和地 距向波数；步骤D :对第2幅子孔径SAR图像的空间二维傅里叶分解形式S (kx，ky)乘上相 位补偿因子进行频域相位补偿，得到频域相位补偿后的子孔径SAR图像Y (人步骤E : 将第1幅子孔径SAR图像的空间二维傅里叶分解形式S+(kx，ky)与进行频域相位补偿后的 第2幅子孔径SAR图像f进行频域相加，得到SAR图像S' (kx，ky);步骤F :对SAR 图像S' (kx，ky)在频域消除相位补偿因子的影响，得到消除海浪杂波干扰的SAR海洋图像 纹理的二维频谱Se (kx，ky);以及步骤G :对消除海浪杂波干扰的SAR海洋图像纹理的二维频 谱H ky)进行二维逆傅里叶变换，得到消除海浪杂波干扰的SAR海洋图像。
[0008] (三）有益效果
[0009] 从上述技术方案可以看出，本发明基于子孔径图像进行海浪杂波抑制的方法利用 了海浪与其他海洋现象在几何特征和动态特征上的差异进行杂波抑制，抑制效果比传统仅 基于几何特征的抑制方法更好。
【附图说明】
[0010] 图1为根据本发明实施例基于子孔径图像的SAR海洋图像中海浪杂波抑制方法的 流程图；
[0011] 图2为图1所示SAR海洋图像中海浪杂波抑制方法的步骤B中子孔径分割的示意 图。
【具体实施方式】
[0012] 本发明同时利用了海浪与其他海洋现象在几何特征和动态特征上的差异进行杂 波抑制，抑制效果比传统仅基于几何特征的抑制方法更好。
[0013] 为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照 附图，对本发明进一步详细说明。
[0014] 假设海面上存在两类纹理特征，一类是海面上快速运动的波浪，另一类是相对慢 速变化的其他纹理（例如内波、溢油等）。在这里为了简单起见，对于慢变的其他纹理假设 其在SAR合成孔径时间内不变。也就是海面SAR图像信号可以分解为：
[0015] s(x, y, t) = se (x, y)+sw(x, y, t) (I)
[0016] 其中x和y分别为图像的方位向和地距向坐标（这里假定SAR图像已经进行了斜 地转换），t为SAR子孔径成像的中心时间，S w(X，y，t)和sjx，y)分别为快速运动的波浪 散射信号和海上慢变的其他海洋现象的散射信号。一般SAR合成孔径时间在分钟或者秒量 级以内，在这么短的时间里可以假设其他慢变海洋现象不发生变化。
[0017] 在本发明的一个示例性实施例中，提供了一种基于子孔径图像进行海浪杂波抑制 的方法。图1为根据本发明实施例基于子孔径图像进行海浪杂波抑制的方法的流程图。如 图1所示，本实施例基于子孔径图像进行海浪杂波抑制的方法包括：
[0018] 步骤A :对SAR原始回波s进行方位向傅里叶变换，得到斜距多普勒域数据S ;
[0019] 步骤B :对斜距多普勒域数据S通过方位向多普勒频谱子孔径分解，得到两幅子孔 径回波数据SJP S2;
[0020] 本发明应用了多普勒与SAR成像时间有一一对应关系的原理。该两幅子孔径回波 数据-SJPS2对于同一个目标区域的成像之间存在一个时间差At。假设两幅子孔径回波 数据-SJPS2的多普勒中心间隔为Afd，根据多普勒频率与成像时间的对应关系，则两幅子 孔径回波数据SJP S 2的时间差Δ t为：
[0021]
C2)
[0022] 其中，λ为雷达波长，R为斜距，Δ fd为两幅子孔径SAR图像的方位向多普勒频谱 中心间隔，V为SAR平台速度。从上式可以看出，要得到固定的时间差，子孔径分割时的多 普勒中心间隔Af d与斜距R必须成反比。
[0023] 具体而言，该步骤B进一步包括：
[0024] 子步骤Bl :设定子孔径间隔时间Δ t，其中Δ t至少要大于波浪振动周期的1/4 ;
[0025] 首先，通过对SAR原始回波s进行常规SAR成像处理，从SAR图像判断图像中主要 波浪纹理的最大波长，再根据如下公式计算波浪振动周期T :
[0026]
(3) \ &
[0027] 其中，1"为波浪中主要波浪纹理的最大波长，g为重力加速度。
[0028] 然后，根据波浪周期T设定一个合理的间隔时间At，一般情况下要保证At > T/4,过短的间隔时间会导致本发明的最终结果信噪比较低。
[0029] 子步骤B2 :在斜距多普勒域中对斜距多普勒域数据S进行分割，不同距离门对应 的子孔径多普勒中心间隔满足下式：
[0030] (4)
[0031] 刮z所示。请参照图2,图中夹在两条实线中间的是子孔径1 的信号，夹在两条虚线中的是子孔径2的信号。在每个距离门上，子孔径1范围内的多普勒 中心与子孔径2范围内的多普勒中心差距Af d满足公式（4)。
[0032] 步骤C :对两幅子孔径回波数据SJP S2，分别进行SAR成像、斜地转换和方位向/地 距向的空间二维傅里叶