水下地形在雷达图像上的可视模型及数值模拟

本文提出了一个在三维空间上潮汐流与海面共同作用区的水下地形在合成孔径雷达图像上的“可视”化原理，得出水下地形成像的综合表面模式，探求在雷达分辨单元内中等尺度的表面波成倾斜对雷达后向散射的贡献。     

星载SAR浅海水下地形和水深测量模拟仿真──水下地形高度、坡度和方向与可测水深分析

根据星载合成孔径雷达(SAR)浅海水下地形和水深成像机理,建立了浅海水下地形和水深雷达后向散射截面仿真模型.利用该模型模拟并分析了不同地形条件下,浅海水下地形的雷达后向散射截面.分析结果表明,水下地形高度越高,SAR可测量的水深越深;水下地形坡度越大,越易被SAR所观测.水下地形的星载SAR测量还与水下地形的方向有关,与卫星飞行方向平行的水下地形最易被SAR观测,与卫星飞行方向垂直的水下地形最不易被SAR观测.     

用合成孔径雷达图像反演浅海水下地形的一种方法

海流与海底地形的相互作用导致了海表面的粗糙起伏,从而引起对海面观测的雷达散射截面的变化。在浅海海流速度垂直分布廓线为均匀的假定下,由流体连续性方程和驰豫时间近似下流体力学弱相互作用的理论,已证明雷达散射截面的空间变化正比于海流方向上海流速度变化的梯度。研究了合成孔径雷达(SAR)对海面的观测图像在无槽道信息和海流方向预知条件下,提出用散射系数的二维相关函数确定海流方向。在海流速度和方向已知条件下,推导了由雷达散射截面的空间变化迭推反演浅海水下地形的公式,并研究了航天飞机SIR-CSAR在我香港特区海域的二维图像反演浅海水下地形的应用。     

Crucial physical characteristics of sea ice in the Arctic section of 143°-180°W during August and early September 2008

潮流是水下地形SAR成像的决定性因素之一,潮流的周期性变化使水下地形SAR成像与潮流场有密切的关系,同一水下地形不同时刻SAR影像不同。以台湾浅滩为例,基于5景不同时相的ERS-2 SAR影像和台湾浅滩实测水深数据,开展了水下地形SAR成像与潮流场的相关性分析。为了排除SAR影像自身成像质量对研究的影响,首先对5景SAR影像的成像质量进行评价,然后通过潮汐潮流数值计算确定SAR成像时刻的潮流场以及基于QuikScat散射计数据确定SAR成像时刻的海面风向。在上述基础上开展了SAR影像后向散射系数与实测水深的相关性分析,总结出水下地形SAR成像与潮流场的关系,结果表明,潮流场与水下地形SAR成像密切相关,当潮流流向与水下地形梯度方向一致或接近时,水下地形SAR成像效果最佳;对于同一水下地形,流向相反的两种潮流对应的SAR影像后向散射系数恰好相反,即SAR影像中的明暗条纹不同。相对于潮流流向的顺流与逆流风向的不同对水下地形SAR成像无影响。     

水下地形SAR成像与潮流场的相关性分析示例

潮流是水下地形SAR成像的决定性因素之一,潮流的周期性变化使水下地形SAR成像与潮流场有密切的关系,同一水下地形不同时刻SAR影像不同。以台湾浅滩为例,基于5景不同时相的ERS-2 SAR影像和台湾浅滩实测水深数据,开展了水下地形SAR成像与潮流场的相关性分析。为了排除SAR影像自身成像质量对研究的影响,首先对5景SAR影像的成像质量进行评价,然后通过潮汐潮流数值计算确定SAR成像时刻的潮流场以及基于QuikScat散射计数据确定SAR成像时刻的海面风向。在上述基础上开展了SAR影像后向散射系数与实测水深的相关性分析,总结出水下地形SAR成像与潮流场的关系,结果表明,潮流场与水下地形SAR成像密切相关,当潮流流向与水下地形梯度方向一致或接近时,水下地形SAR成像效果最佳;对于同一水下地形,流向相反的两种潮流对应的SAR影像后向散射系数恰好相反,即SAR影像中的明暗条纹不同。相对于潮流流向的顺流与逆流风向的不同对水下地形SAR成像无影响。     

一种改进的基于逆Omega-K算法的海面场景SAR原始数据仿真方法

合成孔径雷达(SAR)海面场景原始数据仿真是研究海洋动力参数(表面波浪、风矢量和洋流)的有效工具。目前海面场景原始数据仿真方法已经基于逆Omega-K算法实现了海洋运动参数的空间变化。但是目前仅仅讨论了正侧视情况下的海面场景仿真,应用范围有限,同时没有考虑Stokes漂流以及Bragg相速度的影响,而这两者都是存在于真实海面的。通常情况下为了反演得到海面流场的二维速度矢量,雷达需要从两个不同的方位方向观察海面的同一个区域,因此这就需要考虑大斜视的雷达波束,同时Stokes漂流和Bragg相速度是SAR海表面流场观测不容忽视的两种运动。本文在不改变原有正侧视逆Omega-K算法的情况下,通过增加重新计算零方位时刻的斜视波束中心位置坐标,并据此确定SAR原始数据在多普勒域的位置来将其扩展到大斜侧视逆Omega-K算法,并通过时域Stokes漂流公式到频域内离散化Stokes漂流公式的推导来加入Stokes漂流,以及根据Bragg散射机制加入了Bragg相速度。仿真结果表明,经过聚焦成像后的SAR图像很好的体现了真实海面波浪场的形状,同时能够很好地反演出设定的雷达径向流场速度,且流速精度误差控制在6%以内。最后也证明了Bragg相速度以及Stokes漂流对于海面流场的影响不可忽视。     

降雨对SAR风场反演的影响及校正

合成孔径雷达（Synthetic Aperture Radar, SAR）在风场反演中发挥着重要的作用,但由于受到降雨的影响,海面风场的反演精度会随之下降。本文利用Radarsat-2 SAR观测数据,准同步降雨数据和再分析风场资料对2014年台风"威马逊"外围风场受降雨影响的情况进行了分析,建立了雨致海表阻尼后向散射系数受降雨强度,入射角等因素影响的拟合模型,并对降雨引起的信号衰减,后向体散射进行了评估。实验发现雨致海表面阻尼作用抑制风致海面波的影响十分明显,且降雨造成的信号衰减随降雨强度和入射角的增大而增加,雨滴产生的后向体散射和雨致海表阻尼作用造成的后向散射系数变化随降雨强度的增大而增加,随入射角的增大而减小。实例表明本文建立的模型能够有效改善降雨条件下SAR风场的反演精度。     

海面风浪SAR成像仿真研究

对于SAR风浪图像而言,速度聚束调制常导致图像中散射面元沿方位向产生位置交叠,从而引起强烈非线性特征,难于仿真;另外,SAR图像强度是由雷达飞行过程中多次观测所得的散射回波相干叠加获得,然而,传统SAR海浪图像仿真方法中并未考虑该相干叠加过程。为了解决第一个问题,文中首先应用线性滤波法对海浪运动引起的速度聚束效应进行仿真,统计方位向各像素位置处发生交叠的散射面元数目,并将各散射面元回波进行非相干叠加。针对第二个问题,本文则根据海面微波散射场时间序列的相位差近似满足高斯分布这一事实,将雷达沿方位向飞行过程中不同时刻所得的散射场引入满足高斯分布的随机相位,进而再将方位向不同时刻散射场进行相干叠加。最后通过将仿真所得SAR图像的强度分布与经验分布函数进行比较,初步检验了本文仿真方法的可行性,同时还分别分析了风速、风向、速度聚束调制等因素对仿真SAR图像的影响。     

Kelvin尾迹SAR多视向的成像仿真

利用船只Kelvin航迹模型、海面波模型和二尺度微波散射模型,提出了船尾迹多视向的成像仿真技术,并首次在二维空间中从不同视向仿真船尾迹的SAR图像。结果表明,当雷达视向与船只航向平行时,横波成像明显;当雷达视向与船只航向垂直时,扩散波成像明显;当雷达视向与船只航向有个夹角时,会出现一臂亮一臂暗的现象,这一现象取决于两臂尖波的传播方向与雷达视向的夹角,传播方向与雷达视向越接近平行的波越容易被雷达观测到,从而形成亮臂。仿真结果还得出另外一个结论:船只航向与雷达视向越接近垂直,两臂张角越小。仿真结果和实际的多幅ERS-SAR图像所观测到的结果是一致的。该模型可以有效地模拟Kelvin尾迹SAR多视向成像。     

基于极化特征SERD的SAR溢油检测

与单极化SAR(Synthetic Aperture Radar)相比,全极化SAR图像中不仅包含散射目标的几何特征和后向散射特征,还包含散射目标的极化特征。因此,基于极化特征的SAR图像分类能够更全面地描述海面目标的物理特性。单次反射特征值相对差异度(Single Bounce Eigenvalue Relative Difference,SERD)能够比较单次散射机制的相对大小,并且可以反映散射表面的粗糙度情况。而海面油膜的存在抑制了海面的短重力波和毛细波,改变了海表面的粗糙度。基于此,本文将SERD应用到海面溢油检测中。利用两景Radarsat-2全极化SAR数据对比分析了SERD与极化散射熵的溢油检测效果,实验发现:(1)SERD能够较好地区分溢油与海水。(2)对原油而言,SERD的油水对比度与极化散射熵的油水对比度在数值上差异较小;对生物油膜而言,SERD的油水对比度在数值上远小于极化散射熵。利用这一特性,SERD在区分生物油膜与原油方面更具优势。     

Shallow water depth retrieval from space-borne SAR imagery

Based on shallow water bathymetry synthetic aperture radar (SAR) imaging mechanism and the microwave scattering imaging model for oceanic surface features, we developed a new method for shallow water depth retrieval from space-borne SAR images. The first guess of surface currents and winds are estimated from the normalized radar crossing section (NRCS) profile of shallow water bathymetry SAR imagery, according to the linear theory and geophysical model function. The NRCS profile is then simulated by the microwave scattering imaging model. Both the surface currents and winds are adjusted by using the dichotomy method step by step to make the M4S-simulated NRCS profiles approach those observed by SAR. Then, the surface currents and the wind speeds are retrieved when a best fit between simulated signals and the SAR image appears. Finally, water depths are derived using the Navier–Stokes equation and finite difference method with the best estimated currents and the surface winds. The method is tested on two SAR images of the Taiwan Shoal. Results show that the simulated shallow water NRCS profile is in good agreement with those measured by SAR with the correlation coefficient as high as 85%. In addition, when water depths retrieved from the SAR image are compared with in situ measurements, both the root mean square and relative error are less than 3.0 m and 6.5%, respectively, indicating that SAR images are useful for shallow water depth retrieval and suggesting that the proposed method in this paper is convergent and applicable.     

A Eulerian–Lagrangian method for the simulation of wave propagation

A Eulerian–Lagrangian method (ELM) is employed for the simulation of wave propagation in the present research. The wave action conservation equation, instead of the wave energy balance equation, is used. The wave action is conservative and the action flux remains constant along the wave rays. The ELM correctly accounts for this physical characteristic of wave propagation and integrates the wave action spectrum along the wave rays. Thus, the total derivative for wave action spectrum may be introduced into the numerical scheme and the complicated partial differential wave action balance equation is simplified into an ordinary differential equation. A number of test cases on wave propagation are carried out and show that the present method is stable, accurate and efficient. The results are compared with analytical solutions and/or other computed results. It is shown that the ELM is superior to the first-order upwind method in accuracy, stability and efficiency and may better reflect the complicated dynamics due to the complicated bathymetry features in shallow water areas.     

A Eulerian–Lagrangian method for the simulation of wave propagation

A Eulerian–Lagrangian method (ELM) is employed for the simulation of wave propagation in the present research. The wave action conservation equation, instead of the wave energy balance equation, is used. The wave action is conservative and the action flux remains constant along the wave rays. The ELM correctly accounts for this physical characteristic of wave propagation and integrates the wave action spectrum along the wave rays. Thus, the total derivative for wave action spectrum may be introduced into the numerical scheme and the complicated partial differential wave action balance equation is simplified into an ordinary differential equation. A number of test cases on wave propagation are carried out and show that the present method is stable, accurate and efficient. The results are compared with analytical solutions and/or other computed results. It is shown that the ELM is superior to the first-order upwind method in accuracy, stability and efficiency and may better reflect the complicated dynamics due to the complicated bathymetry features in shallow water areas.     

A Review About SAR Technique for Shallow Water Bathymetry Surveys

Synthetic Aperture Radar(SAR)has become one of the important tools for shallow water bathymetry surveys.This has significant economic efficiency compared with the traditional bathymetry surveys.Numerical models have been developed to simulate shallow water bathymetry SAR images.Inversion of these models makes it possible to assess the water depths from SAR images.In this paper,these numerical models of SAR technique are reviewed,and examples are illustrated including in the coastal areas of China.Some issues about SAR technique available and the research orientation in future are also discussed.     

海浪同化预报模型理论基础

将共轭变分同化方法应用于 LAGFD- WAM海浪数值模式 ,导出了海浪谱能量平衡方程的共轭方程以及风输入、破碎、底摩擦、波波非线性相互作用和波流相互作用的相应共轭源函数 ,建立了海浪同化模型 ,数值计算仍采用特征线嵌入计算格式 ,为合成孔径雷达波谱反演资料和卫星高度计有效波高资料同化奠定理论基础     

Ocean wave parameters retrieved directly from compact polarimetric SAR data

We aim to directly invert wave parameters by using the data of a compact polarimetric synthetic aperture radar(CP SAR) and validate the effectiveness of ocean wave parameter retrieval from the circular transmit/linear receive mode and π/4 compact polarimetric mode. Relevant data from the RADARSAT-2 fully polarimetric SAR on the C-band were used to obtain the compact polarimetric SAR images, and a polarimetric SAR wave retrieval algorithm was used to verify the sea surface wave measurements. Usin...     

Satellite-based RAR performance simulation for measuring directional ocean wave spectrum based on SAR inversion spectrum

Some missions have been carried out to measure wave directional spectrum by synthetic aperture radar (SAR) and airborne real aperture radar (RAR) at a low incidence. Both them have their own advantages and limitations. Scientists hope that SAR and satellite-based RAR can complement each other for the research on wave properties in the future. For this study, the authors aim to simulate the satellite-based RAR system to validate performance for measuring the directional wave spectrum. The principal measurements are introduced and the simulation methods based on the one developed by Hauser are adopted and slightly modified. To enhance the authenticity of input spectrum and the wave spectrum measuring consistency for SAR and satellite-based RAR, the wave height spectrum inversed from Envisat ASAR data by cross spectrum technology is used as the input spectrum of the simulation system. In the process of simulation, the sea surface, backscattering signal, modulation spectrum and the estimated wave height spectrum are simulated in each look direction. Directional wave spectrum are measured based on the simulated observations from 0° to 360°. From the estimated wave spectrum, it has an 180° ambiguity like SAR, but it has no special high wave number cut off in all the direction. Finally, the estimated spectrum is compared with the input one in terms of the dominant wave wavelength, direction and SWH and the results are promising. The simulation shows that satellite-based RAR should be capable of measuring the directional wave properties. Moreover, it indicates satellite-based RAR basically can measure waves that SAR can measure.     

基于单景高分辨率光学卫星遥感影像的浅海水下地形提取

本文基于海浪波折射现象和浅水波理论,提出了一种基于单景高分辨率光学遥感影像的浅海地形提取方法。首先,基于浅水波理论推导出适用于浅海区域的水深与海浪波长、频率的定量关系,针对近岸光学遥感图像复杂的海浪特征,讨论了两种海浪波长提取方法,即FFT方法和剖面线法。然后提出了基于长距离波长波动分析的海浪频率计算方法,解决了单景遥感影像的波浪频率计算难题。最后,利用单景QuickBird高分辨率光学遥感影像,以海南岛三亚湾为研究区域进行了应用实验,结果表明,对12m以浅的浅海区域,在不需要任何辅助参数的情况下,反演获得了浅海地形(DEM),经与1:25000比例尺海图的水深对比验证,地形趋势吻合良好,反演水深的均方根误差为1.07m,相对水深误差为16.2%,表明该方法适合于浅海水下地形的提取,且具有无需实测水深数据和环境参数的支持的优点。     

复杂地形下海浪数值模式的特征线计算格式

LAGFD-WAM海浪数值模式是一种第三代海浪数值模式，通过求解波数谱平衡方程，并考虑风输入、波浪破碎耗散、底摩擦耗散、波波非线性相互作用和波流相互作用等源函数，模拟波数空间下的海浪方向谱，并依此获得海浪的波高、周期和平均波向。该模式的一个显著特点是采用特征线嵌入格式求解海浪的传播。在进行浅水区域的海浪模拟时，特征线嵌入格式的数值计算方案是否合理对海浪数值模拟结果产生直接的影响。为此LAGFD-WAM海浪数值模式提出了一种新的特征线混合数值计算格式，并应用于浅水海浪数值模拟。结果表明，采用该计算方法，能够使数值模拟结果与实测结果很好符合。     

Wash Waves Generated by High Speed Displacement Ships

It is difficult to compute far-field waves in a relative large area by using one wave generation model when a large calculation domain is needed because of large dimensions of the waterway and long distance of the required computing points. Variation of waterway bathymetry and nonlinearity in the far field cannot be included in a ship fixed process either. A coupled method combining a wave generation model and wave propagation model is then used in this paper to simulate the wash waves generated by the passing ship. A NURBS-based higher order panel method is adopted as the stationary wave generation model; a wave spectrum method and Boussinesq-type equation wave model are used as the wave propagation model for the constant water depth condition and variable water depth condition, respectively. The waves calculated by the NURBS-based higher order panel method in the near field are used as the input for the wave spectrum method and the Boussinesq-type equation wave model to obtain the far-field waves. With this approach it is possible to simulate the ship wash waves including the effects of water depth and waterway bathymetry. Parts of the calculated results are validated experimentally, and the agreement is demonstrated. The effects of ship wash waves on the moored ship are discussed by using a diffraction theory method. The results indicate that the prediction of the ship induced waves by coupling models is feasible.