一种海洋涡旋SAR图像仿真方法

王宇航 杨 敏 种劲松(中国科学院电子学研究所 北京 100190)(微波成像技术国家重点实验室 北京 100190)(中国科学院大学 北京 100190)

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图像进行对比，从而验证该方法的有效性。

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

2 涡旋SAR图像仿真方法

应当看到，高校教材编写对推动科研发展也起着重要的作用，因为它将学科发展的宝贵史料、思想与理论精萃、最新科研成果不断选编、整合为学科的概念与原理的知识体系，还力求反映学科的前沿热点与动态。这将为相关科学研究起着整理资源、提高基础、继往开来的铺垫作用，尤其有助于激发和培养学生的科研兴趣与创新能力，以至直接影响他们参与科学研究。因此，我们要高度评价编好教材的重大意义。在教材编写的过程中，我们也深深感到教材编写是一个在正确思想指导下广泛学习，不断深入研究，反复修改与提高的过程。

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

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

2.1 涡旋2维流场仿真

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

图1 涡旋SAR图像仿真方法流程图Fig.1 Flow chart of the simulation method of SAR eddy image

1）通过该改革的开发和研究，使学生能够真正感受、理解工程实际设备及装备的原理与工作过程，明确掌握机械制造工艺基本概念的真实含义，以提高学生学习专业课的兴趣，深刻正确地理解与掌握课程内容，并为课程全面提高教学质量奠定基础。

2.2 涡旋SAR图像仿真

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

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

图2 SAR海洋成像仿真模型示意图[22]Fig.2 Schematic diagram of oceanic SAR imagery simulation model[22]

其中， ± 表示远离雷达方向和朝向雷达方向的两组Bragg波分量， 〈fD±〉σ表示经过归一化后向散射系数 σ加权的平均Doppler中心； γD±表示Doppler谱的方差。〈fD±〉σ和 γD±的具体计算过程可以参考文献[25]，这里不再赘述。

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

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

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

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

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

应构建长效性扶持机制，项目支持不只局限于前期建设阶段，也应向中后期运行阶段延伸。加大政策扶持力度，提高项目资金补助标准，补贴应尽可能弥补秸秆收储人工和运输成本，提高秸秆利用主体积极性。实行区别化补贴和奖励标准，引导产业发展向技术含量高、附加值高的利用方向发展。拓宽政策扶持广度，降低项目申报门槛，尽可能让政策惠及每一个产业链主体，最大限度调动各方面积极性。政府部门需密切协作、形成合力，在财政、税收、土地等各方面完善政策优惠体系，通过完善利益链驱动产业链发展，最终形成以政策为导向、企业为主体、农民广泛参与的秸秆综合利用长效机制。

其中，海面归一化后向散射系数 〈σ 〉由入射角、雷达频率、极化方式、海面风速等参数所决定，NEσ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所示。

铝土矿中w(CaO)≥10%，换算得到w(Ca)≥7.14%，试液中存在的钙，如果不加以去除，易形成硫酸钙沉淀，随渣相排除[9-10]，导致测定结果偏低。实验中加入20mL 100g/L碳酸铵溶液完全可以将试液中的钙沉淀完全，说明碳酸铵可以有效地掩蔽钙的干扰。

从欧洲中期天气预报中心(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。因此，设置参数为–0.003486，流场大小设置为18 km×24 km，空间分辨率为100 m，雷达参数设置为表1中ERS-2 SAR参数。

图3 中国东海海域获取的ERS-2 SAR图像，获取时间为2009.08.19, 02:23:50 UTCFig.3 ERS-2 SAR image of the East China Sea obtained on August 19, 2009 at 02:23:50 UTC

图4 从方框1处截取的涡旋SAR图像Fig.4 Enlargement of the eddy in Frame 1

表1 ERS-2 SAR参数Tab.1 SAR parameters of ERS-2

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

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

图5 相同参数下仿真SAR图像与ERS-2 SAR图像对比图Fig.5 Comparison of simulated SAR image and ERS-2 SAR image under the same parameters

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

图6 仿真SAR图像与ERS-2 SAR图像涡旋信息提取Fig.6 Eddy information extraction of simulated SAR image and ERS-2 SAR image

表2 涡旋信息提取结果Tab.2 Results of eddy information extraction

3.2 反气旋式涡旋仿真实验

图7 吕宋海峡获取的ENVISAT-1 ASAR图像，获取时间为2010.06.11, 01:51:48 UTCFig.7 ENVISAT-1 ASAR image of the Luson Strait obtained on June 11, 2010 at 01:51:48 UTC

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

从ECMWF获取2010.06.11, 03:00:00时刻的风场再分析资料，分辨率为0.125°×0.125°。根据数据显示，涡旋区域附近的风速为2.1 m/s，风向为45°。从GODAS获取相同位置的5日平均流场再分析资料，

护理交接班是临床护理服务的重要工作之一,可有效保证护理工作的连续性,同时对医疗服务中事故发生加以预防[1]。护理晨交班是护士在短时间内需要完成的重要护理工作,不仅需要提高晨交班速度,也要保障传递信息准确无误,使接班者可准确了解患者的病情变化。我科自2017年8月实施sbar交班模式后,使护士在晨交接中患者有效信息传递方面的接受率大幅度提升,同时也使得交班流程更加规范,有效缩短了晨交班时间,提高了交班质量,现汇报如下。

表3 ENVISAT-1 ASAR参数Tab.3 ASAR parameters of ENVISAT-1

3.1节对气旋式涡旋进行了仿真实验，本节将针对反气旋式涡旋进行仿真实验。图7是一幅ENVISAT-1分辨率为(1/3)°×1°。根据数据显示，涡旋区域附近的流速为0.23 m/s。因此，设置参数为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]研究结果一致，初步验证了仿真方法的正确性。

课程是教育思想转化为现实的桥梁和纽带。课程是学校改革与发展、提高教育质量的关键环节。高职院校内涵建设的核心是课程。课程建设是在一定的课程观指导下进行课程体系建设、课程组织实施、课程考核评价以及课程管理等工作的全过程。

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

他给《中国医院院长》杂志记者的名片上，已经印着自己在宁波一院出诊的信息。“社区医院主要以常见病为主，而在宁波一院遇到的病种多、病情急、病情复杂。这样的锻炼，扩展临床思维，积累临床经验”。

图9 相同参数下仿真SAR图像与ENVISAT-1 ASAR图像对比图Fig.9 Comparison of simulated SAR image and ENVISAT-1 ASAR image under the same parameters

图10 仿真SAR图像与ENVISAT-1 ASAR图像涡旋信息提取Fig.10 Eddy information extraction of simulated SAR image and ENVISAT-1 ASAR image

表4 涡旋信息提取结果Tab.4 Results of eddy information extraction

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

4 总结

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

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

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