一种基于径向速度场的SAR海面成像方法技术

本发明专利技术公开了一种基于径向速度场的SAR海面成像方法，通过基于SAR数据获得径向速度场，从而实现海面成像，能够提升海面纹理聚焦效果，进而提升SAR海面图像的可解译性，同时能够提升SAR海洋数据的利用率，有助于提升海面纹理特征，为后续的目标探测奠定基础。为后续的目标探测奠定基础。为后续的目标探测奠定基础。

【技术实现步骤摘要】
一种基于径向速度场的SAR海面成像方法


[0001]本专利技术属于合成孔径雷达(Synthetic Aperture Radar，SAR)信号处理
，具体涉及一种基于径向速度场的SAR海面成像方法。

技术介绍

[0002]合成孔径雷达(Synthetic Aperture Radar，SAR)是一种高分辨率成像雷达，在距离向上，它通过发射宽带脉冲信号获得较高的距离向分辨率；在方位向上，SAR利用雷达平台的运动，通过信号处理形成一个较长的“孔径”，从而获得较高的方位分辨率。与传统光学传感器相比，SAR具有全天时、全天候、分辨率高、观测范围广的优势。海洋观测是SAR的重要应用方向之一，SAR可以探测海浪、内波、涡旋、锋面等多种海洋现象，可以反演风场、流场等矢量信息，而且可以对油膜、舰船进行识别分类。SAR海面成像是对海观测的基础，但是由于海面整体的随机运动，通常给海面成像带来挑战。
[0003]目前常用的SAR成像算法主要有RD算法、CS算法、Omega
‑
K算法以及BP算法。SAR主要用来对地面静止目标或场景成像，当SAR对海观测时必须要考虑海面的随机运动，否则会严重降低SAR图像的质量。
[0004]海面散射点随机运动产生的径向速度分量会引起多普勒中心偏移，这将导致海面散射点位置发生偏移。由于海面整体的随机运动，每个散射点的偏移量均不相同，这将导致不同位置散射点之间可能会累加在一起，进而导致SAR图像中一些微弱的细节信息的丧失，给SAR图像解译、目标识别带来困难。
[0005]综上所述，现有的SAR海面成像算法大多直接应用陆地的成像算法，没有考虑不同海面散射点之间径向速度差异导致的方位向偏移，这将导致在SAR图像中丢失一些海面纹理信息，严重影响SAR图像解译和目标探测。

技术实现思路

[0006]有鉴于此，本专利技术提供了一种基于径向速度场的SAR海面成像方法，能够实现基于SAR数据的海面成像。
[0007]本专利技术提供的一种基于径向速度场的SAR海面成像方法，包括：
[0008]步骤1、根据SAR接收的原始数据生成并保存单视复图像；
[0009]步骤2、将所述单视复图像分割成多个数据子块，估计所有数据子块的多普勒中心频率；
[0010]步骤3、计算各数据子块对应的径向速度，由所有数据子块的径向速度组成径向速度场；
[0011]步骤4、将所述单视复图像进行方位向傅里叶变换，得到距离
‑
多普勒域数据，并根据所述径向速度场对所述距离
‑
多普勒域数据进行相位补偿；
[0012]步骤5、将补偿后的距离
‑
多普勒域数据变换到二维时域，完成海面成像。
[0013]进一步地，所述步骤1中生成并保存单视复图像的过程中采用无权匹配滤波器实
现脉冲压缩。
[0014]进一步地，所述数据子块的尺寸为100m*100m。
[0015]进一步地，所述数据子块的多普勒中心频率由斜视角产生的分量和动目标产生的分量组成，表示为：f
dc
＝f
dc_θ
+f
dc_tar
，其中，f
dc_θ
表示斜视角产生的多普勒中心偏移，f
dc_tar
表示动目标产生的多普勒中心偏移。
[0016]进一步地，所述数据子块对应的径向速度采用如下公式计算：
[0017][0018]其中，V
r
(i,j)表示第(i,j)个数据子块的径向速度，λ为雷达波长；
[0019]由此所述径向速度场表示为：
[0020][0021]进一步地，所述步骤4中根据所述径向速度场对所述距离
‑
多普勒域数据进行相位补偿的补偿因子为：其中，V为径向速度场，λ为雷达波长，j为虚数单位，t为慢时间变量。
[0022]进一步地，所述步骤5中将补偿后的距离
‑
多普勒域数据变换到二维时域的方式为：将补偿后的距离
‑
多普勒域数据进行方位傅里叶逆变换。
[0023]有益效果：
[0024]本专利技术通过基于SAR数据获得径向速度场，从而实现海面成像，能够提升海面纹理聚焦效果，进而提升SAR海面图像的可解译性，同时能够提升SAR海洋数据的利用率，有助于提升海面纹理特征，为后续的目标探测奠定基础。
附图说明
[0025]图1为本专利技术提供的一种基于径向速度场的SAR海面成像方法的流程图。
[0026]图2为本专利技术提供的一种基于径向速度场的SAR海面成像方法的详细流程图。
[0027]图3为实施例中采用的仿真点目标示意图。
[0028]图4(a)为采用本专利技术提供的一种基于径向速度场的SAR海面成像方法中采用Omega
‑
K算法获取的成像结果图。
[0029]图4(b)为采用本专利技术提供的一种基于径向速度场的SAR海面成像方法获取的成像结果。
具体实施方式
[0030]下面结合附图并举实施例，对本专利技术进行详细描述。
[0031]受海面随机运动的影响，海面不同散射点之间的径向速度均不相同，这将导致不
同散射点在成像过程中导致的多普勒中心偏移不同，进而导致位置的偏移产生差异。由于不同散射点位置偏移之间的差异，采用传统的SAR成像算法将会导致目标的混叠，致使一些信息的丢失。
[0032]本专利技术提供的一种基于径向速度场的SAR海面成像方法，流程如图1所示，包括以下步骤：
[0033]步骤1、根据SAR接收的原始数据生成并保存单视复图像，其中数据块大小为M*N。
[0034]为了准确的估计多普勒中心频率进而获得精确的速度场，本专利技术在生成单视复图像时，采用无权匹配滤波器完成匹配滤波过程，实现脉冲压缩。采用RD、Omega
‑
K等SAR成像算法对SAR接收到的原始数据进行处理，得到单视复图像(Single
‑
Look Complex，SLC)。
[0035]步骤2、将步骤1得到的单视复图像分割成多个数据子块，得到m*n个数据子块，分别估计每个数据子块的多普勒中心频率。
[0036]将单视复图像分割成多个子块，得到m*n个数据子块，其中每个子块的大小根据真实数据决定，一般可以选取100m*100m。估计m*n个数据子块的多普勒中心频率f
dc
，多普勒中心频率f
dc
由斜视角产生的分量和动目标产生的分量组成，可以表示为：
[0037]f
dc
＝f
dc_θ
+f
dc_tar
ꢀꢀꢀ
(1)。
[0038]其中，f
dc_θ
表示斜视角产生的多普勒中心偏移，表示为：
[0039][0040]V为载机飞行速度，θ为斜视角，λ为雷达波长；f
dc_tar
表示动目本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于径向速度场的SAR海面成像方法，其特征在于，包括：步骤1、根据SAR接收的原始数据生成并保存单视复图像；步骤2、将所述单视复图像分割成多个数据子块，估计所有数据子块的多普勒中心频率；步骤3、计算各数据子块对应的径向速度，由所有数据子块的径向速度组成径向速度场；步骤4、将所述单视复图像进行方位向傅里叶变换，得到距离
‑
多普勒域数据，并根据所述径向速度场对所述距离
‑
多普勒域数据进行相位补偿；步骤5、将补偿后的距离
‑
多普勒域数据变换到二维时域，完成海面成像。2.根据权利要求1所述的一种基于径向速度场的SAR海面成像方法，其特征在于，所述步骤1中生成并保存单视复图像的过程中采用无权匹配滤波器实现脉冲压缩。3.根据权利要求1所述的一种基于径向速度场的SAR海面成像方法，其特征在于，所述数据子块的尺寸为100m*100m。4.根据权利要求1所述的一种基于径向速度场的SAR海面成像方法，其特征在于，所述数据子块的多普勒中心频率由斜视角产生的分量和动目标产生的分量组成，表示为：f
dc

【专利技术属性】
技术研发人员：赵亚威，种劲松，孙凯，杨雪，李宗泽，魏献恩，
申请(专利权)人：中国科学院空天信息创新研究院，
类型：发明
国别省市：