基于子孔径图像的SAR海洋图像中海浪杂波抑制方法技术

本发明专利技术提供了一种基于子孔径图像的SAR海洋图像中海浪杂波抑制方法。本发明专利技术利用了海浪与其他海洋现象在几何特征和动态特征上的差异进行杂波抑制，抑制效果比传统仅基于几何特征的抑制方法更好，具有较好地推广应用价值。

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

【技术保护点】
一种基于子孔径图像的SAR海洋图像中海浪杂波抑制方法，其特征在于，包括：步骤A：对SAR原始回波s进行方位向傅里叶变换，得到斜距多普勒域数据S；步骤B：对斜距多普勒域数据S通过方位向多普勒频谱子孔径分解，得到两幅子孔径回波数据S1和S2；步骤C：对两幅子孔径回波数据S1和S2，分别进行SAR成像、斜地转换和方位向/地距向的空间二维傅里叶变换，得到两幅子孔径SAR图像的空间二维傅里叶分解形式S+(kx，ky)，S‑(kx，ky)，其中kx，ky分别为方位向和地距向波数；步骤D：对第2幅子孔径SAR图像的空间二维傅里叶分解形式S‑(kx，ky)乘上相位补偿因子进行频域相位补偿，得到频域相位补偿后的子孔径SAR图像步骤E：将第1幅子孔径SAR图像的空间二维傅里叶分解形式S+(kx，ky)与进行频域相位补偿后的第2幅子孔径SAR图像进行频域相加，得到SAR图像S′(kx，ky)；步骤F：对SAR图像S′(kx，ky)在频域消除相位补偿因子的影响，得到消除海浪杂波干扰的SAR海洋图像纹理的二维频谱Se(kx，ky)；以及步骤G：对消除海浪杂波干扰的SAR海洋图像纹理的二维频谱Se(kx，ky)进行二维逆傅里叶变换，得到消除海浪杂波干扰的SAR海洋图像。...

【技术特征摘要】
1.一种基于子孔径图像的SAR海洋图像中海浪杂波抑制方法，其特征在于，包括：步骤A：对SAR原始回波s进行方位向傅里叶变换，得到斜距多普勒域数据S；步骤B：对斜距多普勒域数据S通过方位向多普勒频谱子孔径分解，得到两幅子孔径回波数据S1和S2；步骤C：对两幅子孔径回波数据S1和S2，分别进行SAR成像、斜地转换和方位向/地距向的空间二维傅里叶变换，得到两幅子孔径SAR图像的空间二维傅里叶分解形式S+(kx,ky)，S-(kx,ky)，其中kx，ky分别为方位向和地距向波数；步骤D：对第2幅子孔径SAR图像的空间二维傅里叶分解形式S-(kx,ky)乘上相位补偿因子进行频域相位补偿，得到频域相位补偿后的子孔径SAR图像步骤E：将第1幅子孔径SAR图像的空间二维傅里叶分解形式S+(kx,ky)与进行频域相位补偿后的第2幅子孔径SAR图像进行频域相加，得到SAR图像S′(kx,ky)；步骤F：对SAR图像S′(kx,ky)在频域消除相位补偿因子的影响，得到消除海浪杂波干扰的SAR海洋图像纹理的二维频谱Se(kx,ky)；以及步骤G：对消除海浪杂波干扰的SAR海洋图像纹理的二维频谱Se(kx,ky)进行二维逆傅里叶变换，得到消除海浪杂波干扰的SAR海洋图像。2.根...

【专利技术属性】
技术研发人员：王小青，高骥超，种劲松，
申请(专利权)人：中国科学院电子学研究所，
类型：发明
国别省市：北京;11