电控时变电磁材料的SAR成像特性研究

王俊杰; 冯德军; 王志凇; 邢世其; 李永祯; 王雪松

doi:10.12000/JR21104

电控时变电磁材料的SAR成像特性研究

DOI: 10.12000/JR21104

1.
国防科技大学电子科学学院电子信息系统复杂电磁环境效应国家重点实验室长沙 410073
2.
中国人民解放军32200部队锦州 121000

基金项目: 国家自然科学基金(62071475, 61890542, 62001481)

详细信息

作者简介:
王俊杰(1991–)，男，博士，讲师，主要研究方向为新材料电磁调控、雷达信号处理与新型电子对抗

冯德军(1972–)，男，研究员，博士生导师，主要研究方向为智能电子对抗、精确制导与目标识别

王志凇(1988–)，男，硕士，主要研究方向为武器装备实验鉴定

邢世其(1984–)，男，博士，副研究员，主要研究方向为极化雷达成像、雷达信号处理与合成孔径雷达对抗

李永祯(1977–)，男，研究员，博士生导师，主要研究方向为极化雷达与电子对抗

王雪松(1972–)，男，教授，博士生导师，主要研究方向为极化雷达、目标识别与电子对抗

通讯作者:
王俊杰 wangjunjie14@vip.qq.com

冯德军 fdj117@sina.com

责任主编：丁泽刚 Corresponding Editor: DING Zegang
中图分类号: TN955
计量
- 文章访问数: 1398
- HTML全文浏览量: 1012
- PDF下载量: 245
- 被引次数: 4
出版历程
- 收稿日期: 2021-07-22
- 修回日期: 2021-10-14
- 网络出版日期: 2021-11-04
- 刊出日期: 2021-12-28

Synthetic Aperture Rader Imaging Characteristics of Electronically Controlled Time-varying Electromagnetic Materials

1.
The State Key Laboratory of Complex Electromagnetic Environment Effects on Electronics and Information System, College of Electronic Science, National University of Defense Technology, Changsha 410073, China
2.
Unit 32200 of PLA, Jinzhou 121000, China

Funds: The National Nature Science Foundation of China (62071475, 61890542, 62001481)

More Information

Corresponding author: WANG Junjie, wangjunjie14@vip.qq.com; FENG Dejun, fdj117@sina.com

摘要

摘要: 微动目标因具有在不同方向的运动分量，能够对电磁波进行微多普勒调制，导致目标成像特征出现方位向散焦效应，这种现象在目标识别与反识别领域被广泛关注与研究。相较而言，电控时变电磁材料通过外加激励实现对电磁波特征的灵活调控，具有更快的调制速度，而其成像特性没有被过多关注。该文以此为切入点，对电控时变电磁材料的合成孔径雷达(SAR)图像距离向调制特性进行了研究，分析了时变电磁材料谱变换模型和SAR目标特征控制原理。以相位调制表面(PSS)为代表，建立了非周期PSS相位调制模型，其频谱具有连续频移特性。在此基础上，探讨了PSS连续频移调制对SAR的影响，揭露了距离向目标散焦现象。通过SAR实测数据仿真，验证了所提理论方法的有效性。
- 时变电磁材料 /
- 相位调制表面 /
- 非周期调制 /
- 谱变换 /
- 合成孔径雷达
Abstract: Micromotion targets can enable microDoppler modulation on electromagnetic waves owing to their motion components in different directions, resulting in the defocusing effects of the target imaging features along the azimuth direction. This phenomenon has been widely considered and investigated in target recognition and antirecognition. As a research hotspot, electronically controlled time-varying electromagnetic materials have fast modulation characteristics; however, their imaging characteristics have not received much attention. This study investigates the Synthetic Aperture Radar (SAR) modulation characteristics along the range direction of electronically controlled time-varying electromagnetic materials. Then, this study establishes the time-varying electromagnetic material spectrum transformation model and analyzes the SAR target characteristic control principle. Taking the Phase-Switched Screen (PSS) as an example, a nonperiodic PSS phase modulation model is established, and its spectrum has continuous frequency shift characteristics. On this basis, the influence of continuous frequency shift modulation induced by nonperiodic PSS on SAR was further discussed and the defocusing phenomenon of the target along the range direction was detected. Through the simulation of measured data, the validity of the proposed theoretical method is verified.
- Time-varying electromagnetic materials /
- Phase-Switched Screen (PSS) /
- Nonperiodic modulation /
- Spectrum conversion /
- Synthetic Aperture Radar (SAR)

HTML全文

1. 引言

海洋涡旋是一种旋转的、以封闭环流为主要特征的水体，是由于各种气象因素作用和海洋动力不稳定性形成的。作为一种重要的海洋现象，涡旋不仅能够影响海洋流场与化学物质的输送，从而对海洋的环流结构和海洋生态等产生重要作用，还能通过海气相互作用，对风场、云及降雨等大气现象产生影响^[1,2]。

合成孔径雷达(Synthetic Aperture Radar, SAR)具有全天时、全天候、高分辨率、广覆盖面等优点，对海洋涡旋探测具有特殊意义，受到国际海洋遥感界的重视。然而，涡旋在SAR成像时会受到各种海洋环境因素的影响，通过真实SAR图像难以完全解译涡旋的特征。利用仿真SAR图像可以为涡旋的SAR图像特征解译提供指导，但是目前利用SAR图像对涡旋的研究主要集中在涡旋的统计性研究^[3–5]、涡旋的形成机制和成像分类^[6–8]以及涡旋的检测和特征提取方面的研究^[9–12]，极少有关于涡旋SAR图像仿真方法的研究。

由于海面随机运动且电磁散射特性复杂，难以进行时间和空间上的SAR原始回波仿真。海浪谱能够描述随机海面不同波长海浪的能量分布情况，因此利用海浪谱可以很好地描述不同海况下随机海面的统计特征。海洋涡旋、内波、浅海地形、锋面等都可看作是通过波流交互作用，即利用自身流场改变海浪谱分布，并经过海面电磁散射模型，进而得以在SAR图像上体现。目前SAR海面图像仿真常用的电磁散射模型包括Kirchhoff散射模型、Bragg散射模型以及组合表面散射模型^[13]。这些模型只考虑了1阶Bragg散射，仅适用于低频(小于L波段)SAR海面图像仿真。1997年，Romeiser和Alpers^[14,15]提出了改进的组合表面模型，该模型考虑了2阶Bragg波散射的影响，从而使仿真的SAR海面图像更接近实际情况。2002年，Romeiser^[16]利用该模型研究了浅海地形在SAR图像上的特征，并与声学多普勒流速剖面仪测量的浅海地形进行对比，验证了该模型用于SAR海面图像仿真的合理性；2011年，欧阳越等^[17]利用该模型仿真了不同雷达参数下海洋内波图像，并同实际内波SAR图像进行对比，发现二者具有较高的一致性。但是目前，利用海面电磁散射模型对海洋涡旋SAR图像仿真的研究尚未见报道。

为此，本文提出了一种海洋涡旋SAR图像仿真方法，利用流体力学中典型的Burgers-Rott涡旋模型，建立涡旋的2维流场。利用SAR海洋成像仿真模型，仿真涡旋SAR图像。基于此方法，本文进行了气旋式涡旋与反气旋式涡旋SAR图像仿真实验，并将仿真SAR图像与ERS-2 SAR图像和ENVISAT-1 ASAR图像进行对比，从而验证该方法的有效性。

2. 涡旋SAR图像仿真方法

本文建立的涡旋SAR图像仿真方法，是在给定2维涡旋流场和风场条件下，利用SAR海洋成像模型生成随机海面的2维海浪谱，再根据2维海浪谱与SAR图像之间的调制传递函数，生成仿真涡旋SAR图像。

涡旋SAR图像仿真方法分为两步，如图1所示。首先，输入涡旋流场参数，基于涡旋动力学模型建立涡旋2维流场(于2.1节介绍)。然后，将仿真的涡旋流场和海面风场输入到SAR海洋成像仿真模型，通过设置SAR参数获得仿真涡旋SAR图像(于2.2节介绍)。

图 1 涡旋SAR图像仿真方法流程图

Figure 1. Flow chart of the simulation method of SAR eddy image

下载: 全尺寸图片幻灯片

2.1 涡旋2维流场仿真

涡旋一般遵循流体力学的纳维-斯托克斯(Navier-Stokes，简写N-S)方程，根据方程中黏性力项、惯性力项以及离心力项的平衡关系，可以建立不同的涡旋模型。常见的涡旋模型包括Rankine涡旋、Oseen涡旋、Sullivan涡旋以及Burgers-Rott涡旋^[18–20]，其中，Rankine涡旋模型没有考虑N-S方程中的黏性力项，流体以常角速度 $\omega$ 旋转，没有径向速度，因而不能产生涡旋的辐散、辐聚和上升运动；Oseen涡旋模型仅考虑N-S方程中惯性力项的局地项及黏性力项，其轨道是一个圆形涡旋，不符合实际SAR图像中涡旋的形态；Sullivan涡旋模型和Burgers-Rott涡旋模型考虑了N-S方程中全部的黏性力项、惯性力项及离心力项，但由于Sullivan涡旋模型的轨道是一个双螺旋涡旋，Burgers-Rott涡旋模型的轨道是一个螺旋形涡旋，后者与真实SAR图像所呈现的涡旋形状更为接近，因此本文选用Burgers-Rott涡旋模型来建立海洋涡旋的流场。

Burgers-Rott涡旋模型是从N-S方程求得的一个涡旋解^[19,20]，假定涡旋是定常和轴对称的，涡旋速度场在柱坐标系下表示为

$\left\{ \begin{array}{l} {V_r} = \frac{{dr}}{{dt}} = - \frac{\alpha }{2}r\\ {V_\theta } = r\frac{{d\theta }}{{dt}} = \frac{{{\varGamma _0}}}{{2{\text{π}} r}}\left(1 - {{\rm e}^{ - \displaystyle\scriptsize\frac{{\alpha {r^2}}}{{4\upsilon }}}}\right)\\ {V_z} = \frac{{dz}}{{dt}} = \alpha z \end{array} \right.$

(1)

其中， ${V_r},{V_\theta },{V_z}$ 分别是 $r,\theta ,z$ 方向的速度分量， $\alpha$ 为吸入强度， $\upsilon$ 为黏性系数， ${\varGamma _0}$ 是 $r \to \infty$ 时的速度环量， $\varGamma {\rm{ = }}2{\text{π}} r{v_\theta }$ 。

将式(1)转化为直角坐标系，涡旋速度场可表示为

$\left\{ \begin{array}{l} {V_x} = - \frac{\alpha }{2}x - \frac{{{\varGamma _0}\alpha }}{{8{\text{π}} \upsilon }}y\\ {V_y} = \frac{{{\varGamma _0}\alpha }}{{8{\text{π}} \upsilon }}x - \frac{\alpha }{2}y \end{array} \right.$

(2)

其中， ${V_x}$ 为涡旋速度场在 $x$ 方向上的速度分量， ${V_y}$ 是涡旋速度场在 $y$ 方向上的速度分量。

通过设置参数 $\alpha{\text{、}}\ {\varGamma _0}/\upsilon$ 的值，根据式(2)可以得到涡旋2维流场。通过仿真发现， $\alpha$ 的值会影响涡旋流场流速的大小， $\alpha$ 的值越大，涡旋流场流速越大，反之则越小； $\alpha$ 的正负影响涡旋流场的旋向， $\alpha$ 为正，流场顺时针旋转， $\alpha$ 为负，流场逆时针旋转； ${\varGamma _0}{\rm{/}}\upsilon$ 的值则会影响涡旋臂的曲率， ${\varGamma _0}{\rm{/}}\upsilon$ 的值越大，涡旋臂的曲率越大。

2.2 涡旋SAR图像仿真

获得了涡旋的流场之后，下一步将进行涡旋SAR图像的仿真。本文使用SAR海洋成像仿真模型来仿真涡旋SAR图像。SAR海洋成像仿真模型主要分为波流交互作用模型、雷达后向散射模型和SAR成像模型3个部分，如图2所示。

图 2 SAR海洋成像仿真模型示意图^[22]

Figure 2. Schematic diagram of oceanic SAR imagery simulation model^[22]

下载: 全尺寸图片幻灯片

首先，将仿真的涡旋2维流场和海面风场输入到波流交互作用模型，通过求解作用量谱平衡方程，计算给定海面流场和海面风场下被调制的海浪谱。作用量谱平衡方程如式(3)所示^[21]：

$\begin{align} & \frac{{\partial N(\text{x}{\rm{,}}\text{k},t)}}{{\partial t}} + \left[ {{c_g}(\text{k}) + U(\text{x},t)} \right]\frac{{\partial N(\text{x},\text{k},t)}}{{\partial \text{x}}} \\ & \quad\quad - \text{k}\frac{{\partial U(\text{x},t)}}{{\partial \text{x}}}\frac{{\partial N(\text{x},\text{k},t)}}{{\partial \text{k}}} = S(\text{x},\text{k},t) \end{align}$

(3)

其中， $N$ 为微尺度波作用量谱密度， $\text{x} = (x,y)$ 为空间位置矢量， $\text{k} = ({k_x},{k_y})$ 为波数矢量， $U$ 为表面流场， ${c_g}$ 为被调制波浪的群速度， $S$ 为源函数(风场输入、非线性波-波作用和弥散等作用之和)，在本文模型中采用的源函数表达式为

$S(\text{x},\text{k}{\rm{,}}t) = \mu (\text{k})N(\text{x},\text{k}{\rm{,}}t)\left( {1 - \frac{{N(\text{x},\text{k}{\rm{,}}t)}}{{{N_0}(\text{k})}}} \right)$

(4)

其中， ${N_0}$ 为不存在海流时平衡状态下的作用量谱密度， $\mu$ 为松弛率。

作用量谱密度与海浪谱的关系为^[23]

$N(\text{x},\text{k}{\rm{,}}t) = \frac{{\rho {\omega _0}(\text{k})}}{\text{k}}\psi (\text{x},\text{k}{\rm{,}}t)$

(5)

其中， ${\omega _0}(\text{k}) = \sqrt {g\text{k} + (\tau/\rho ){\text{k}^3}}$ , $\tau$ 为表面张力， $\rho$ 为海水密度， $\psi$ 为海浪谱。

令 $Q(\text{x},\text{k}{\rm{,}}t) = 1/N(\text{x},\text{k}{\rm{,}}t)$ , ${Q_0}(\text{k}) = 1/{N_0}(\text{k})$ ，则被流场调制后的海浪谱为

$\begin{align} \frac{{\psi (\text{x},\text{k}{\rm{,}}t)}}{{{\psi _0}(\text{k})}} &= \frac{{{Q_0}(\text{k})}}{{{Q_0}(\text{k}) + \delta Q(\text{x},\text{k}{\rm{,}}t)}} \\ &= \frac{1}{{1 + \displaystyle\frac{{\delta Q(\text{x},\text{k}{\rm{,}}t)}}{{{Q_0}(\text{k})}}}} \end{align}$

(6)

其中， $\delta Q$ 表示调制引起的作用量谱变化量。

然后，将计算得到的海浪谱输入到雷达后向散射模型，在给定雷达频率、入射角、极化方式及雷达视向等雷达参数下，仿真涡旋SAR图像后向散射强度。本文采用的雷达后向散射模型为改进的组合表面模型，是Romeiser和Alpers等在Bragg共振散射模型基础上的改进^[14,15]。该模型同时考虑了长波和中波对短波的倾斜调制和水动力调制，所以从理论和试验研究上更能表现海面微波散射的实际情况，是目前最为完善的海面微波散射模型之一。该模型是基于2维海面坡度，通过傅里叶变换对后向散射截面进行泰勒级数展开，并对后向散射截面进行时间和空间上的平均。由于1阶项平均后为0，因此得到2阶Bragg散射后的海面归一化后向散射系数为^[14,22]

$\begin{align} \left\langle \sigma \right\rangle =& {\sigma ^{(0)}} + \left\langle {{\sigma ^{(2)}}} \right\rangle \\ =& \sigma {|_{{\text{s}} = 0}} + \iint {\left( {{{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm p}}}\limits^ \wedge \partial \mathop {{s_{\rm p}}}\limits^ \vee }}} \right|}_{\mathop {\text{s}}\limits^ \wedge = 0}} + {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm p}}}\limits^ \vee \partial \mathop {{s_{\rm p}}}\limits^ \wedge }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}}} \right)}\\ & \cdot k_{\rm p}^2\psi (\text{k}){{\rm d}^2}\text{k}\\ & + \iint {\left( {{{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm n}}}\limits^ \wedge \partial \mathop {{s_{\rm n}}}\limits^ \vee }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}} + {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm n}}}\limits^ \vee \partial \mathop {{s_{\rm n}}}\limits^ \wedge }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}}} \right)}\\ & \cdot k_{\rm n}^2\psi (\text{k}){{\rm d}^2}\text{k}\\ &+ \iint \left( {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm p}}}\limits^ \wedge \partial \mathop {{s_{\rm n}}}\limits^ \vee }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}} + {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \wedge \vee } }}{{\partial \mathop {{s_{\rm n}}}\limits^ \vee \partial \mathop {{s_{\rm p}}}\limits^ \wedge }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}}\right. \\ & \left. + {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \vee \wedge } }}{{\partial \mathop {{s_{\rm p}}}\limits^ \vee \partial \mathop {{s_{\rm n}}}\limits^ \wedge }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}} + {{\left. {\frac{{{\partial ^2}\mathop \sigma \limits^{ \vee \wedge } }}{{\partial \mathop {{s_{\rm n}}}\limits^ \vee \partial \mathop {{s_{\rm p}}}\limits^ \wedge }}} \right|}_{\mathop {{\text{s}}}\limits^ \wedge = 0}} \right) \\ & \cdot {k_{\rm p}}{k_{\rm n}}\psi (\text{k}){{\rm d}^2}\text{k} \end{align}$

(7)

其中， ${\sigma ^{(0)}}$ 为平静海面的归一化后向散射系数； $\left\langle {{\sigma ^{(2)}}} \right\rangle$ 表示表面坡度引起的2阶Bragg散射之和；符号 $\left\langle \cdot \right\rangle$ 表示统计平均； $\text{s} = \left( {{s_{\rm{p}}},{s_{\rm{n}}}} \right)$ 为海面坡度； ${k_{\rm{p}}}$ , ${k_{\rm{n}}}$ 分别为平行和垂直于雷达视向的Bragg波波数分量； $\psi \left( \text{k} \right)$ 为海浪波数谱；符号 $\wedge$ 和 $\vee$ 分别表示 $\sigma$ 对波数 $\text{k}$ 的傅里叶变换及其共轭； $\wedge \wedge$ 和 $\vee \vee$ 表示 $\sigma$ 对组合波数 ${k_1} + {k_2}$ 的傅里叶变换及其共轭； $\wedge \vee$ 和 $\vee \wedge$ 表示 $\sigma$ 对组合波数 ${k_1} - {k_2}$ 的傅里叶变换及其共轭。

上述过程中，利用海浪谱与雷达后向散射模型得到仿真的海面归一化后向散射系数，但这是一个实孔径雷达成像过程，SAR图像仿真还需考虑海面运动的影响。当目标存在沿雷达视线方向的径向速度时，将在方位向上产生偏移 $\Delta x$ ：

$\Delta x = - \frac{R}{V\ \, }{v_{\rm r}} = \frac{{R\lambda }}{{2V\ }}{f_{\rm D}}$

(8)

其中， $R$ 是雷达至目标的距离， $V\;$ 是平台飞行速度， $\lambda$ 是雷达波长， ${f_{\rm{D}}} = - \displaystyle\frac{{2{v_{\rm{r}}}}}{\lambda }$ 是目标速度导致的Doppler谱中心偏移。对于海面而言，由于其各点速度不同，在方位向上偏移量不同，导致SAR海面图像产生压缩或拉伸的现象，即速度聚束效应^[24]。此外，分辨单元内不同散射点速度的分布方差将造成回波Doppler谱展宽，并导致分辨率下降。SAR成像模型通过计算每个分辨单元的平均Doppler谱中心和方差引入海面运动对SAR成像造成的影响。

这里采用Romeiser和Thompson^[25]给出的双高斯形Doppler谱模型计算Doppler谱中心和方差，该模型将海面回波Doppler谱分成朝向雷达和远离雷达两个传播方向的Bragg波Doppler谱的叠加，每个Doppler谱分量为高斯形，其具体表达式为

$\begin{align} W({f_{\rm D}}) =& \frac{{ \left\langle{\sigma _ + } \right\rangle }}{{\sqrt {2{\text{π}} \gamma _{{\rm D} + }^2} }}{{\rm e}^{ - {{\left({f_{\rm D}} - {\left\langle{f_{{\rm D} + }}\right\rangle_\sigma }\right)}^2}/\gamma _{{\rm D} + }^2}} \\ & + \frac{{ \left\langle{\sigma _ - } \right\rangle }}{{\sqrt {2{\text{π}} \gamma _{{\rm D} - }^2} }}{{\rm e}^{ - {{\left({f_{\rm D}} - {\left\langle{f_{{\rm D} - }} \right\rangle_\sigma }\right)}^2}/\gamma _{{\rm D} - }^2}} \end{align}$

(9)

其中， $\pm$ 表示远离雷达方向和朝向雷达方向的两组Bragg波分量， $\left\langle{f_{{\rm{D}} \pm }} \right\rangle _\sigma$ 表示经过归一化后向散射系数 $\sigma$ 加权的平均Doppler中心； ${\gamma _{{\rm D} \pm }}$ 表示Doppler谱的方差。 $\left\langle{f_{{\rm D} \pm }}\right\rangle _\sigma$ 和 ${\gamma _{{\rm D} \pm }}$ 的具体计算过程可以参考文献[25]，这里不再赘述。

另外，仿真的SAR图像还需考虑噪声的影响，本文涡旋SAR图像仿真过程中，仅考虑热噪声对仿真SAR图像信噪比的影响。信噪比由噪声等效后向散射系数以及海面归一化后向散射系数所决定：

${\rm{SNR(dB) = }}\left\langle \sigma \right\rangle - NE{\sigma ^0}$

(10)

其中，海面归一化后向散射系数 $\left\langle \sigma \right\rangle$ 由入射角、雷达频率、极化方式、海面风速等参数所决定， $NE{\sigma ^0}$ 为噪声等效后向散射系数，由系统硬件参数所决定。因此，SAR成像模型根据给定的仿真输入参数计算信噪比，从而得到具有统计特性的仿真涡旋SAR图像。

3. 涡旋SAR图像仿真实验

根据涡旋旋转方向的不同，可将涡旋分为气旋式涡旋与反气旋式涡旋^[26]。气旋式涡旋在北半球逆时针旋转，在南半球顺时针旋转；反气旋式涡旋在北半球顺时针旋转，在南半球逆时针旋转。不同旋转方向的涡旋将产生不同的涡旋流场，从而在SAR图像中呈现不同的涡旋特征。下面，本文分别针对气旋式涡旋与反气旋式涡旋进行仿真实验。

3.1 气旋式涡旋仿真实验

图3是一幅ERS-2 SAR图像，图像获取时间为2009.08.19, 02:23:50 UTC，获取地点为中国东海海域。图中方框1处为一个气旋式涡旋，旋转方向为逆时针。为了便于对比仿真SAR图像与真实SAR图像，将方框1处的涡旋截取出来，截取图像尺寸为18 km×24 km，如图4所示。ERS-2 SAR图像的具体雷达参数如表1所示。

图 3 中国东海海域获取的ERS-2 SAR图像，获取时间为2009.08.19, 02:23:50 UTC

Figure 3. ERS-2 SAR image of the East China Sea obtained on August 19, 2009 at 02:23:50 UTC

下载: 全尺寸图片幻灯片

图 4 从方框1处截取的涡旋SAR图像

Figure 4. Enlargement of the eddy in Frame 1

下载: 全尺寸图片幻灯片

表 1 ERS-2 SAR参数

Table 1. SAR parameters of ERS-2

参数	数值
极化方式	VV
波段	C
入射角	23.0°
平台高度	780 km
平台速度	7500 m/s

下载: 导出CSV

| 显示表格

从欧洲中期天气预报中心(Europe Centre for Medium-Range Weather Forecasts, ECMWF)获取2009.08.19, 03:00:00时刻的风场再分析资料，分辨率为0.125°×0.125°。根据数据显示，涡旋区域附近的风速为1.4 m/s，风向为257.9°。从全球海洋数据同化系统(Global Ocean Data Assimilation System, GODAS)获取相同位置的5日平均流场再分析资料，分辨率为(1/3)°×1°。根据数据显示，涡旋区域附近的流速为0.61 m/s。因此，设置参数 $\alpha$ 为–0.003486，流场大小设置为18 km×24 km，空间分辨率为100 m，雷达参数设置为表1中ERS-2SAR参数。

图5(a)、图5(b)分别是该涡旋的仿真SAR图像与获取的真实SAR图像，仿真时设定的雷达参数、海面风场条件与真实SAR图像获取条件完全一致。对比图5(a)、图5(b)两图发现，仿真SAR图像与真实SAR图像中的涡旋臂形状几乎一致，涡旋臂的亮暗特征也基本吻合。从逆时针方向看，涡旋臂由外到内的亮暗特征均为亮-暗-亮，这种亮暗特征的变化是由雷达后向散射引起的布拉格波谱密度变化导致的^[7]。这初步验证了仿真方法的正确性。

图 5 相同参数下仿真SAR图像与ERS-2 SAR图像对比图

Figure 5. Comparison of simulated SAR image and ERS-2 SAR image under the same parameters

下载: 全尺寸图片幻灯片

为了进一步验证仿真方法的正确性，定量地描述仿真SAR图像与真实SAR图像的中涡旋的相似程度，采用文献[9]中基于对数螺旋线边缘拟合的SAR图像涡旋信息提取方法，提取仿真SAR图像和真实SAR图像中涡旋的中心位置、直径及边缘长度，并加以比较。拟合及提取结果如图6所示，红色加号表示涡旋中心位置，黄色箭头表示涡旋直径，蓝色曲线表示涡旋边缘，具体数值如表2所示。

图 6 仿真SAR图像与ERS-2 SAR图像涡旋信息提取

Figure 6. Eddy information extraction of simulated SAR image and ERS-2 SAR image

下载: 全尺寸图片幻灯片

表 2 涡旋信息提取结果

Table 2. Results of eddy information extraction

SAR图像	涡旋中心位置	涡旋直径	涡旋边缘长度
仿真SAR图像	(116,75)	18.9 km	35.7 km
真实SAR图像	(113,71)	18.7 km	35.4 km
绝对/相对误差	(3,4)/—	0.2 km/0.011	0.3 km/0.008

下载: 导出CSV

| 显示表格

对比仿真SAR图像与真实SAR图像的涡旋信息提取结果，可以发现两幅图像中涡旋的中心位置较为一致，方位向和距离向上仅相差3～4个像素点，涡旋直径及边缘长度的相对误差均不超过0.011，证明本文提出的基于Burgers-Rott涡旋模型的涡旋SAR图像仿真方法能够实现气旋式涡旋的SAR图像仿真，并且仿真SAR图像与真实SAR图像能够较好地吻合。

3.2 反气旋式涡旋仿真实验

3.1节对气旋式涡旋进行了仿真实验，本节将针对反气旋式涡旋进行仿真实验。图7是一幅ENVISAT-1 ASAR图像，图像获取时间为2010.06.11, 01:51:48 UTC，获取地点在吕宋海峡。图中方框2处为一个反气旋式涡旋，旋转方向为顺时针。将方框2处的涡旋截取出来，截取图像尺寸为24 km×24 km，如图8所示。ENVISAT-1 ASAR图像的具体雷达参数如表3所示。

图 7 吕宋海峡获取的ENVISAT-1 ASAR图像，获取时间为2010.06.11, 01:51:48 UTC

Figure 7. ENVISAT-1 ASAR image of the Luson Strait obtained on June 11, 2010 at 01:51:48 UTC

下载: 全尺寸图片幻灯片

图 8 方框2处截取的涡旋SAR图像

Figure 8. Enlargement of the eddy in Frame 2

下载: 全尺寸图片幻灯片

表 3 ENVISAT-1 ASAR参数

Table 3. ASAR parameters of ENVISAT-1

参数	数值
极化方式	HH
波段	C
入射角	26.7°
平台高度	800 km
平台速度	7455 m/s

下载: 导出CSV

| 显示表格

从ECMWF获取2010.06.11, 03:00:00时刻的风场再分析资料，分辨率为0.125°×0.125°。根据数据显示，涡旋区域附近的风速为2.1 m/s，风向为45°。从GODAS获取相同位置的5日平均流场再分析资料，分辨率为(1/3)°×1°。根据数据显示，涡旋区域附近的流速为0.23 m/s。因此，设置参数 $\alpha$ 为0.000657，流场大小设置为24 km×24 km，空间分辨率为100 m，雷达参数设置为表3中ENVISAT-1 ASAR参数。

图9(a)、图9(b)分别是该涡旋的仿真SAR图像与真实SAR图像，仿真时设定的雷达参数、海面风场条件与真实SAR图像完全一致。对比图9(a)、图9(b)两图发现，仿真SAR图像与真实SAR图像中的涡旋形状基本一致，涡旋臂的亮暗特征也较为吻合。从顺时针方向看，涡旋臂从外到内均呈现为由暗到亮的特征。该结果与Lyzenga等人^[6]研究结果一致，初步验证了仿真方法的正确性。

图 9 相同参数下仿真SAR图像与ENVISAT-1 ASAR图像对比图

Figure 9. Comparison of simulated SAR image and ENVISAT-1 ASAR image under the same parameters

下载: 全尺寸图片幻灯片

为了定量地描述仿真SAR图像与真实SAR图像的中涡旋的相似性，同样采用3.1节中的分析方法，得到涡旋拟合结果如图10所示，提取的涡旋信息如表4所示。

图 10 仿真SAR图像与ENVISAT-1 ASAR图像涡旋信息提取

Figure 10. Eddy information extraction of simulated SAR image and ENVISAT-1 ASAR image

下载: 全尺寸图片幻灯片

表 4 涡旋信息提取结果

Table 4. Results of eddy information extraction

SAR图像	涡旋中心位置	涡旋直径	涡旋边缘尺寸
仿真SAR图像	(144,78)	24.0 km	49.4 km
真实SAR图像	(147,81)	23.9 km	49.7 km
绝对/相对误差	(3,3)/—	0.1 km/0.004	0.3 km/0.006

下载: 导出CSV

| 显示表格

对比仿真SAR图像与真实SAR图像的涡旋信息提取结果，可以发现两幅图像中涡旋的中心位置较为接近，方位向和距离向上均相差3个像素点，涡旋直径及边缘尺寸相对误差均不超过0.006，这进一步验证了仿真方法的正确性，说明本文提出的基于Burgers-Rott涡旋模型的涡旋SAR图像仿真方法能够实现反气旋式涡旋的SAR图像仿真。

4. 总结

本文基于Burgers-Rott涡旋模型，提出了一种涡旋SAR图像仿真方法，并分别针对气旋式涡旋与反气旋式涡旋进行了仿真实验。通过将仿真SAR图像与真实SAR图像对比验证发现，本文提出的涡旋SAR图像仿真方法能够实现气旋式涡旋和反气旋式涡旋的SAR图像仿真，且仿真SAR图像与真实SAR图像能够较好地吻合。

通过涡旋SAR图像仿真实验发现，无论气旋式涡旋还是反气旋式涡旋，其涡旋臂在SAR图像中都会呈现亮暗交替变化的特征。其中，气旋式涡旋臂呈现两个亮暗交替周期，即亮-暗-亮；反气旋式涡旋臂呈现一个亮暗交替周期，即暗-亮。这是由于这两个涡旋臂的曲率不同，气旋式涡旋臂曲率较大，亮暗交替周期较多，反气旋式涡旋臂曲率较小，亮暗交替周期较少。

由于涡旋在SAR成像时会受到各种海洋环境因素的影响，通过真实SAR图像难以完全解译涡旋的特征。本文提出的SAR图像仿真方法能够弥补这种不足，可以清晰地获取涡旋的尺度、亮暗等特征，这为海洋涡旋特征的解译和提取提供了便利。

图 1 电控时变电磁材料

Figure 1. Electronically controlled time-varying electromagnetic materials

下载: 全尺寸图片幻灯片

图 2 SAR成像模型

Figure 2. SAR imaging model

下载: 全尺寸图片幻灯片

图 3 非周期PSS相位调制波形

Figure 3. Nonperiodic PSS phase modulation waveform

下载: 全尺寸图片幻灯片

图 4 非周期PSS波形频谱特性

Figure 4. Spectrum characteristics of nonperiodic PSS waveform

下载: 全尺寸图片幻灯片

图 5 非周期PSS调制的匹配滤波输出

Figure 5. Nonperiodic PSS modulated matched filter output

下载: 全尺寸图片幻灯片

图 6 基于非周期PSS调制的SAR成像效果图

Figure 6. SAR imaging effect based on aperiodic PSS modulation

下载: 全尺寸图片幻灯片

参考文献(32)

[1]	CHEN V C, LI F, HO S S, et al. Micro-Doppler effect in radar: Phenomenon, model, and simulation study[J]. IEEE Transactions on Aerospace and Electronic Systems, 2006, 42(1): 2–21. doi: 10.1109/TAES.2006.1603402
[2]	庄钊文, 刘永祥, 黎湘. 目标微动特性研究进展[J]. 电子学报, 2007, 35(3): 520–525. doi: 10.3321/j.issn:0372-2112.2007.03.028 ZHUANG Zhaowen, LIU Yongxiang, and LI Xiang. The achievements of target characteristic with micro-motion[J]. Acta Electronica Sinica, 2007, 35(3): 520–525. doi: 10.3321/j.issn:0372-2112.2007.03.028
[3]	吴晓芳, 刘阳, 王雪松, 等. 旋转微动目标的SAR成像特性分析[J]. 宇航学报, 2010, 31(4): 1181–1189. doi: 10.3873/j.issn.1000-1328.2010.04.039 WU Xiaofang, LIU Yang, WANG Xuesong, et al. Analysis of SAR imaging characteristics of targets with rotational micro-motion[J]. Journal of Astronautics, 2010, 31(4): 1181–1189. doi: 10.3873/j.issn.1000-1328.2010.04.039
[4]	白雪茹, 孙光才, 周峰, 等. 基于旋转角反射器的ISAR干扰新方法[J]. 电波科学学报, 2008, 23(5): 867–872. doi: 10.13443/j.cjors.2008.05.019 BAI Xueru, SUN Guangcai, ZHOU Feng, et al. A novel ISAR jamming method based on rotating angular reflectors[J]. Chinese Journal of Radio Science, 2008, 23(5): 867–872. doi: 10.13443/j.cjors.2008.05.019
[5]	PENDRY J B. Focus issue: Negative refraction and metamaterials[J]. Optics Express, 2003, 11(7): 639. doi: 10.1364/OE.11.000639
[6]	YEN T J, PADILLA W J, FANG N, et al. Terahertz magnetic response from artificial materials[J]. Science, 2004, 303(5663): 1494–1496. doi: 10.1126/science.1094025
[7]	PADILLA W J. Group theoretical description of artificial electromagnetic metamaterials[J]. Optics Express, 2007, 15(4): 1639–1646. doi: 10.1364/OE.15.001639
[8]	TENNANT A and CHAMBERS B. Adaptive radar absorbing structure with PIN diode controlled active frequency selective surface[J]. Smart materials and Structures, 2004, 13(1): 122–125. doi: 10.1088/0964-1726/13/1/013
[9]	HUANG Cheng, ZHAO Bo, SONG Jiakun, et al. Active transmission/absorption frequency selective surface with dynamical modulation of amplitude[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(6): 3593–3598. doi: 10.1109/TAP.2020.3037813
[10]	CHAMBERS B and TENNANT A. The phase-switched screen[J]. IEEE Antennas and Propagation Magazine, 2004, 46(6): 23–27. doi: 10.1109/MAP.2004.1396733
[11]	RAMACCIA D, SOUNAS D L, ALÙ A, et al. Phase-induced frequency conversion and Doppler effect with time-modulated metasurfaces[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(3): 1607–1617. doi: 10.1109/TAP.2019.2952469
[12]	SAIKIA M, SRIVASTAVA K V, and RAMAKRISHNA A S. Frequency-shifted reflection of electromagnetic waves using a time-modulated active tunable frequency-selective surface[J]. IEEE Transactions on Antennas and Propagation, 2020, 68(4): 2937–2944. doi: 10.1109/TAP.2019.2951494
[13]	ZHAO Jie, YANG Xi, DAI Junyan, et al. Programmable time-domain digital-coding metasurface for non-linear harmonic manipulation and new wireless communication systems[J]. National Science Review, 2019, 6(2): 231–238. doi: 10.1093/nsr/nwy135
[14]	JIA Yongtao, LIU Ying, GUO Y J, et al. A dual-patch polarization rotation reflective surface and its application to ultra-wideband RCS reduction[J]. IEEE Transactions on Antennas and Propagation, 2017, 65(6): 3291–3295. doi: 10.1109/TAP.2017.2694879
[15]	ZAKER R and SADEGHZADEH A. A low-profile design of polarization rotation reflective surface for wideband RCS reduction[J]. IEEE Antennas and Wireless Propagation Letters, 2019, 18(9): 1794–1798. doi: 10.1109/LAWP.2019.2930130
[16]	CUI Tiejun, QI Meiqing, WAN Xiang, et al. Coding metamaterials, digital metamaterials and programmable metamaterials[J]. Light:Science & Applications, 2014, 3(10): e218. doi: 10.1038/lsa.2014.99
[17]	ZHANG Lei, CHEN Xiaoqing, LIU Shuo, et al. Space-time-coding digital metasurfaces[J]. Nature Communications, 2018, 9(1): 4334. doi: 10.1038/s41467-018-06802-0
[18]	杨欢欢, 曹祥玉, 高军, 等. 可重构电磁超表面及其应用研究进展[J]. 雷达学报, 2021, 10(2): 206–219. doi: 10.12000/JR20137 YANG Huanhuan, CAO Xiangyu, GAO Jun, et al. Recent advances in reconfigurable metasurfaces and their applications[J]. Journal of Radars, 2021, 10(2): 206–219. doi: 10.12000/JR20137
[19]	LUO Yong, KIKUTA K, HAN Zhengli, et al. An active metamaterial antenna with mems-modulated scanning radiation beams[J]. IEEE Electron Device Letters, 2016, 37(7): 920–923. doi: 10.1109/LED.2016.2565559
[20]	BRYANT A T, LU Liqing, SANTI E, et al. Physical modeling of fast p-i-n diodes with carrier lifetime zoning, Part I: Device model[J]. IEEE Transactions on Power Electrons, 2008, 23(1): 189–197. doi: 10.1109/TPEL.2007.911823
[21]	XU Letao, FENG Dejun, and WANG Xuesong. Matched-filter properties of linear-frequency-modulation radar signal reflected from a phase-switched screen[J]. IET Radar, Sonar & Navigation, 2016, 10(2): 318–324. doi: 10.1049/iet-rsn.2015.0182
[22]	XU Letao, FENG Dejun, and WANG Xuesong. High-resolution range profile deception method based on phase-switched screen[J]. IEEE Antennas and Wireless Propagation Letters, 2016, 15: 1665–1668. doi: 10.1109/LAWP.2016.2521778
[23]	WANG Junjie, FENG Dejun, XU Zhiming, et al. Time-domain digital-coding active frequency selective surface absorber/reflector and its imaging characteristics[J]. IEEE Transactions on Antennas and Propagation, 2021, 69(6): 3322–3331. doi: 10.1109/TAP.2020.3037757
[24]	保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 123–134. BAO Zheng, XING Mengdao, WANG Tong. Radar Imaging Technique[M]. Beijing: Publishing House of Electronics Industry, 2005: 123–134.
[25]	YASIR S, 杨国敏, 徐丰. 四叶草形超宽带漫散射编码超表面[J]. 雷达学报, 2021, 10(3): 382–390. doi: 10.12000/JR21061 YASIR S, YANG Guomin, and XU Feng. A four-leaf clover-shaped coding metasurface for ultra-wideband diffusion-like scattering[J]. Journal of Radars, 2021, 10(3): 382–390. doi: 10.12000/JR21061
[26]	张然, 冯德军, 徐乐涛. 基于Salisbury屏的二面角设计及其极化特性分析[J]. 雷达学报, 2016, 5(6): 658–665. doi: 10.12000/JR16055 ZHANG Ran, FENG Dejun, and XU Letao. Design and polarization characteristics analysis of dihedral based on Salisbury screen[J]. Journal of Radars, 2016, 5(6): 658–665. doi: 10.12000/JR16055
[27]	WEI Yue, SHEN Gangxiang, and BOSE S K. Span-restorable elastic optical networks under different spectrum conversion capabilities[J]. IEEE Transactions on Reliability, 2014, 63(2): 401–411. doi: 10.1109/TR.2014.2313806
[28]	TENNANT A. Reflection properties of a phase modulating planar screen[J]. Electronics Letters, 1997, 33(21): 1768–1769. doi: 10.1049/el:19971160
[29]	CHAMBERS B and TENNANT A. Reflection of radar signals from multiple phase-modulated surfaces[J]. IET Radar, Sonar & Navigation, 2007, 1(2): 142–148. doi: 10.1049/iet-rsn:20060078
[30]	SHULMAN N and FEDER M. Random coding techniques for nonrandom codes[J]. IEEE Transactions on Information Theory, 1999, 45(6): 2101–2104. doi: 10.1109/18.782147
[31]	TAI Ning, CUI Kaibo, WANG Chao, et al. The design of a novel coherent noise jammer against LFM radar[J]. IEICE Electronics Express, 2016, 13(21): 1–12. doi: 10.1587/elex.13.20160924
[32]	WU Qihua, LIU Xiaobin, LIU Jin, et al. A radar imaging method using nonperiodic interrupted sampling linear frequency modulation signal[J]. IEEE Sensors Journal, 2018, 18(20): 8294–8302. doi: 10.1109/JSEN.2018.2865531

施引文献

期刊类型引用(3)

1.	阮航，崔家豪，毛秀华，任建迎，罗镔延，曹航，李海峰. SAR目标识别对抗攻击综述：从数字域迈向物理域. 雷达学报. 2024(06): 1298-1326 . 本站查看
2.	翟多才，黄贤俊，李元鑫，徐延林，姚理想，林铭团，刘培国. 一种角度与极化稳定的可调控电磁散射结构的设计与验证. 电波科学学报. 2023(02): 205-210 . 百度学术
3.	朱泳庚，艾夏，王伟东，李猛猛，丁大志. 新型多维电磁特性可调控雷达角反射器设计方法研究. 电波科学学报. 2023(06): 1015-1021+1056 . 百度学术

其他类型引用(1)

资源附件(0)

访问统计

图(6)

计量

文章访问数: 1398
HTML全文浏览量: 1012
PDF下载量: 245
被引次数: 4

1. 引言
2. 涡旋SAR图像仿真方法
2.1 涡旋2维流场仿真
2.2 涡旋SAR图像仿真
3. 涡旋SAR图像仿真实验
3.1 气旋式涡旋仿真实验
3.2 反气旋式涡旋仿真实验
4. 总结

电控时变电磁材料的SAR成像特性研究

DOI: 10.12000/JR21104

通讯作者:
王俊杰 wangjunjie14@vip.qq.com

冯德军 fdj117@sina.com

计量

Synthetic Aperture Rader Imaging Characteristics of Electronically Controlled Time-varying Electromagnetic Materials

Corresponding author: WANG Junjie, wangjunjie14@vip.qq.com; FENG Dejun, fdj117@sina.com

1. 引言

2. 涡旋SAR图像仿真方法

2.1 涡旋2维流场仿真

2.2 涡旋SAR图像仿真

3. 涡旋SAR图像仿真实验

3.1 气旋式涡旋仿真实验

3.2 反气旋式涡旋仿真实验

4. 总结

期刊类型引用(3)

其他类型引用(1)

计量

目录

1. 引言

2. 涡旋SAR图像仿真方法

2.1 涡旋2维流场仿真

2.2 涡旋SAR图像仿真

3. 涡旋SAR图像仿真实验

3.1 气旋式涡旋仿真实验

3.2 反气旋式涡旋仿真实验

4. 总结

期刊介绍

联系我们

电控时变电磁材料的SAR成像特性研究

DOI: 10.12000/JR21104

通讯作者: 王俊杰 wangjunjie14@vip.qq.com 冯德军 fdj117@sina.com

计量

出版历程

Synthetic Aperture Rader Imaging Characteristics of Electronically Controlled Time-varying Electromagnetic Materials

Corresponding author: WANG Junjie, wangjunjie14@vip.qq.com; FENG Dejun, fdj117@sina.com

1. 引言

2. 涡旋SAR图像仿真方法

2.1 涡旋2维流场仿真

2.2 涡旋SAR图像仿真

3. 涡旋SAR图像仿真实验

3.1 气旋式涡旋仿真实验

3.2 反气旋式涡旋仿真实验

4. 总结

期刊类型引用(3)

其他类型引用(1)

计量

出版历程

目录

1. 引言

2. 涡旋SAR图像仿真方法

2.1 涡旋2维流场仿真

2.2 涡旋SAR图像仿真

3. 涡旋SAR图像仿真实验

3.1 气旋式涡旋仿真实验

3.2 反气旋式涡旋仿真实验

4. 总结

期刊介绍

联系我们

通讯作者:
王俊杰 wangjunjie14@vip.qq.com

冯德军 fdj117@sina.com