基于SAR图像速度聚束调制的海浪反演研究

本文首先对合成孔径雷达（SAR）海浪成像中的3种调制（倾斜调制、流体力学调制与速度聚束调制）的影响进行了对比分析,结果显示:速度聚束调制对SAR图像的影响最为显著。另外,由于SAR图像中固有相干斑噪声的存在,较低波数范围的噪声难以滤除或抑制,利用经典MPI方法反演海浪谱会造成低波数范围谱值偏大。基于此,本文借鉴经典MPI海浪谱反演算法,建立了基于速度聚束调制的海浪方位向斜率谱和有效波高的反演算法。通过将经典MPI方法、同极化调制法及本文算法等3种海浪反演方法所得有效波高与浮标数据进行比较,结果显示:本文方法反演得到的海浪有效波高与浮标数据获得的有效波高之间的均方误差为0.79 m,为3种方法中最小。     

基于同步风场的机载SAR海浪参数反演方法

针对机载合成孔径雷达(SAR)对海探测特点,采用多入射角法从SAR数据本身得到与海浪参数反演区域时空匹配的同步海面风速和风向,并结合线性变换关系,计算得到海浪初猜谱对应的仿真SAR图像谱,将仿真SAR图像谱和观测SAR图像谱输入代价函数中进行迭代运算,通过非线性方程的解算得到最适海浪谱;采用交叉谱法去除海浪传播180°方向模糊,最终得到海浪参数。论文提出的基于同步风场的机载SAR海浪参数反演方法,充分利用了机载SAR海洋环境探测的优势,解决了传统SAR海浪参数反演中初猜谱构造依赖外部风场的问题,机载同步飞行试验的海浪参数反演结果与浮标观测值的有效波高、波向的均方根误差分别为0.23 m和13.23°,验证了该方法的有效性,可为机载SAR海浪参数反演业务化提供支持。     

星载合成孔径雷达海浪图像谱仿真研究

利用Hasselm ann 的完全非线性积分方法对文氏海浪谱的星载合成孔径雷达(SAR)图像谱进行了仿真研究。结果表明: (1) 在星载SAR海浪成像过程中, 速度聚束调制远比水动力调制、倾斜调制和距离向聚束调制重要;(2) 图像谱有180°的方向模糊和方位向的高波数截断,距离向传播海浪的图像谱存在双峰现象; (3) 不同极化方式和雷达视向的选取对图像谱无明显影响, 因而对于海浪的星载SAR遥感而言, 可以采用单一极化和单一视向。对于不同成长状态的风浪而言, 以上结果保持不变。     

对ENVISAT ASAR level2算法固有误差的分析

欧洲空间局的ENVISAT ASAR level 2算法是从合成孔径雷达(SAR)单视复图像反演涌浪方向谱的算法.该算法假设双峰海浪谱的SAR图像交叉谱是涌浪的图像交叉谱和风浪的图像交叉谱之和.实际上双峰海浪谱的SAR图像交叉谱中还有一个混合项,正是该混合项导致ENVI-SAT ASAR level 2算法有固有误差.利用遥感仿真的方法分析了不同海况条件下该算法的这一固有误差,结果表明,只有在有效波高较小、或风浪的成分较少、或双峰海浪的传播方向较靠近SAR距离向、或波长较长时固有误差才较小,ENVISAT ASAR level 2算法对海浪谱的反演才较为适用.     

ASAR波模式数据反演参数误差与海浪条纹清晰度的相关性分析

欧洲环境卫星-高级合成孔径雷达(EnvironmentalSatellite-AdvancedSyntheticAperture Radar,Envisat-ASAR)波模式数据提供了全球风、浪要素信息,在海浪模式预报与同化方面有重要作用。该数据合成孔径雷达(SyntheticApertureRadar,SAR)图像普遍存在海浪条纹清晰度不同的现象,但是否影响数据精度尚无定论。本文通过比较2010年NODC (the National Oceanographic Date Center)浮标观测数据和波模式数据,发现经过官方修正后的海浪参数反而具有更大误差。进而通过对比不同条纹清晰度的SAR图像反演参数误差,揭示了ASAR产品海浪参数与浮标测量值之间的误差与海浪条纹清晰度的关系。结果表明:海浪条纹清晰的SAR图像的主波波长和主波周期的反演误差更小,而条纹不清晰SAR图像的有效波高和风速的反演误差更小。通过分析海浪参数对海浪条纹清晰度的敏感性,证实了有效波高和方位向截断波长对SAR图像条纹清晰度的响应最好,波陡次之,与卫星飞行方位角和入射角无关。因此,在反演和修正SAR波模式数据时,考虑图像的条纹清晰度,将会有效提高反演数据的精度。该研究可为高分三号等卫星的波模式数据波浪要素反演精度的提升提供有价值的参考。     

基于准线性近似方法的SAR图像海浪参数反演研究

本文根据相干斑噪声的时间快变特征和非海浪纹理现象的时间缓变特征，基于交叉谱提出了一种对相干斑噪声和大尺度非海浪纹理的抑制的方法，进而结合SAR图像谱和海浪谱之间的准线性映射关系，基于SAR数据对海浪参数进行了反演。在反演过程中，首先仿真分析了不同海况下准线性近似法的海浪反演能力，结果表明：风浪引起的方位向截断效应会显著影响反演精度，因此该方法在低风速时的涌浪反演精度更高。通过将基于Sentinel-1卫星2020年的波模式SAR数据的反演结果与欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts, ECMWF)提供的再分析数据进行对比，发现高海况海浪有效波高反演结果明显偏低，而且该反演误差与风速、方位向截断波长之间存在显著相关性。为了提高有效波高的反演精度，本文进一步给出了海浪有效波高反演误差与风速、方位向截断波长之间的经验校正函数模型，结果显示，通过该模型修正后的海浪有效波高反演结果与ECMWF数据和浮标测量数据具有良好一致性。     

基于RADARSAT-2四极化SAR影像的海面风速反演

选取成像于中国东海和南海海域的102景RADARSAT-2(RS-2)精细四极化SAR原始影像以及ERA-Interim风场,在各极化通道下利用相应风速反演模型开展风速反演的研究。CMOD5地球物理模式函数对VV极化SAR数据反演风速效果最好,可以获取海面高精度风速数据。对于HH和交叉极化SAR影像,利用RS-2 SAR影像和ERA-Interim资料对现有风速反演模型进行改进,并提出了用于HH极化风速反演的PR_CS极化比模型和用于交叉极化风速反演的CP_CS模型。结果显示:基于两种模型的SAR反演风速与ERA-Interim风速具有良好的一致性,利用PR_CS模型的SAR反演风速与参考风速的均方根误差为1.54 m/s,利用CP_CS模型的SAR反演风速与ERA-Interim风速的均方根误差分别为1.43 m/s(VH极化)和1.51 m/s(HV极化)。     

波谱仪调制功率谱估计算法性能的仿真研究

基于星载波谱仪海浪方向谱探测原理，仿真了不同海况、风速下的海浪波谱仪接收信号，并采用周期图法、Welch法、AR模型法以及最小方差法共4种不同的调制谱估计方法反演出海浪谱，比较各种调制谱估计方法的海浪方向谱反演性能。仿真结果表明:对于一定方位向下的一维海浪谱反演，不同调制谱估计方法反演海浪谱性能优劣没有绝对的顺序。对于二维海浪谱反演，在成长中海浪条件下周期图法反演性能最差，其他3种方法的反演性能没有绝对优劣顺序；对于成熟风浪，最小方差法在积分能量误差、有效波高误差两个指标上的反演性能最好，在主波波向、主波波长误差上，周期图法反演性能最差，其他3种方法没有绝对优劣顺序。在涌浪条件下，AR模型法反演性能优于其他3种方法。在不同海况下，随着波谱仪中心入射角的下降，反演性能会下降。基于这些仿真结果，本文推荐最小方差法作为充分成长海浪的海浪方向谱反演的调制谱估计方法，AR模型法作为涌浪海浪方向谱反演的调制谱估计方法。     

GF-3号SAR卫星遥感围填海监测方法研究——以大连金州湾为例

依据不同围填海类型在高分三号(GF-3)合成孔径雷达(synthetic aperture radar,SAR)卫星遥感影像上可分性,建立围填海遥感分类体系及相应的围填海类型解译标志,进而分析SAR图像岸线提取方法,构建GF-3围填海监测技术流程。采用几何主动轮廓模型进行GF-3 SAR影像自动提取海岸线,获得围填海专题图。通过外业精度调查验证GF-3 SAR卫星遥感影像可以有效获取围填海信息。     

完善海洋遥感理论 发展海洋遥感技术——评《星载合成孔径雷达海洋遥感导论》

正合成孔径雷达(Synthetic Aperture Radar,SAR)作为一种主动式微波侧视成像雷达,具备全天时、全天候、高分辨率、宽刈幅成像的优势,是可见光、红外传感器及其他微波传感器所不具备的。自美国于1978年发射第一颗海洋遥感卫星SEASAT以来,SAR海洋遥感理论、SAR海洋遥感探测技术与应用研究始终是海洋学研究的热点问题之一。20世纪90年代,航天飞机成像雷达以及X波段SAR(SIR-C/X-SAR)运行之后,多波段、多极化、多入射角和多观测模式成为SAR海洋遥感研究的主流发展方向。SAR遥感图像     

利用一种改进的算法对中国高分三号同极化SAR图像反演的有效波高的验证

Chinese Gaofen-3(GF-3) is the first civilian satellite to carry C-band(5.3 GHz) synthetic aperture radar(SAR).During the period of August 2016 to December 2017, 1 523 GF-3 SAR images acquired in quad-polarization(vertical-vertical(VV), horizontal-horizontal(HH), vertical-horizontal(VH), and horizontal-vertical(HV)) mode were recorded, mostly around China's seas. In our previous study, the root mean square error(RMSE) of significant wave height(SWH) was found to be around 0.58 m when compared with retrieval results from a few GF-3 SAR images in co-polarization(VV and HH) with moored measurements by using an empirical algorithm CSAR_WAVE. We collected a number of sub-scenes from these 1 523 images in the co-polarization channel,which were collocated with wind and SWH data from the European Centre for Medium-Range Weather Forecasts(ECMWF) reanalysis field at a 0.125° grid. Through the collected dataset, an improved empirical wave retrieval algorithm for GF-3 SAR in co-polarization was tuned, herein denoted as CSAR_WAVE2. An additional 92 GF-3 SAR images were implemented in order to validate CSAR_WAVE2 against SWH from altimeter Jason-2, showing an about 0.52 m RMSE of SWH for co-polarization GF-3 SAR. Therefore, we conclude that the proposed empirical algorithm has a good performance for wave retrieval from GF-3 SAR images in co-polarization.     

基于高分四号卫星的黄海绿潮漂移速度提取研究

静止轨道卫星高分四号（GF-4）具有高时间分辨率（20 s）和高空间分辨率（50 m）的独特优势。为了挖掘GF-4卫星在海洋灾害监测中的应用潜力，本文基于2016年6月25日1天4景的GF-4卫星影像，利用最大相关系数法（MCC），开展了黄海绿潮漂移速度提取研究，分析了海面风场、潮汐等对绿潮漂移的影响。研究发现:（1）MCC方法可高精度自动追踪GF-4影像中绿潮的分钟级（8~9 min）位置变化，绿潮漂移速率和方向的相对偏差分别为11%和5%；当2景GF-4影像的成像时间间隔增大至小时级（如6 h）时，随着绿潮斑块形状的改变，MCC方法绿潮自动追踪的准确性下降。（2）绿潮在1天之中的漂移速率和方向可发生显著变化，当日上午9时黄海绿潮漂移速率均值为（0.36±0.13）m/s，方向以东南向为主，至15时，绿潮漂移速率显著增加至（0.69±0.12）m/s，方向变为东北偏北。（3）绿潮漂移速度与海面风速的相关系数为0.74，绿潮漂移方向为风向偏右；绿潮的向岸、离岸运动与相应时刻的涨、落潮具有较好的对应关系。GF-4卫星数据可为绿潮快速漂移的高精度监测提供数据支撑。     

台风中降雨对中国高分3号合成孔径雷达观测信号的影响

Gaofen-3(GF-3), a Chinese civil synthetic aperture radar(SAR) at C-band, has operated since August 2016.Remarkably, several typhoons have been captured by GF-3 around the China Seas over its last two-year mission.In this study, six images acquired in Global Observation(GLO) and Wide ScanSAR(WSC) modes at verticalvertical(VV) polarization channel are discussed. This work focuses on investigating the observation of rainfall using GF-3 SAR. These images were collocated with winds from the European Centre for Medium-Range Weather Forecasts(ECMWF), significant wave height simulated from the WAVEWATCH-III(WW3) model, sea surface currents from climate forecast system version 2(CFSv2) of the National Centers for Environmental Prediction(NCEP) and rain rate data from the Tropical Rainfall Measuring Mission(TRMM) satellite. Sea surface roughness,was compared with the normalized radar cross section(NRCS) from SAR observations, and indicated a 0.8 correlation(COR). We analyzed the dependences of the difference between model-simulated NRCS and SARmeasured NRCS on the TRMM rain rate and WW3-simulated significant wave height. It was found that the effects of rain on SAR damps the radar signal at incidence angles ranging from 15° to 30°, while it enhances the radar signal at incidence angles ranging from 30° to 45° and incidence angles smaller than 10°. This behavior is consistent with previous studies and an algorithm for rain rate retrieval is anticipated for GF-3 SAR.     

Synthetic aperture radar imaging of ocean waves during the marineland experiment

X- andL-band simultaneously obtained synthetic aperture radar (SAR) data of ocean gravity waves collected during the Marineland Experiment were analyzed using wave contrast measurements. The Marineland data collected in 1975 represents a unique historical data set for testing still-evolving theoretical models of the SAR ocean wave imaging process. The wave contrast measurements referred to are direct measurements of the backscatter variation between wave crests and troughs. These modulation depth measurements, which are indicators of wave detectability, were made as a function of: a) the settings used in processing the SAR signal histories to partially account for wave motion; b) wave propagation direction with respect to radar look direction for bothX- andL-band SAR data; c) SAR resolution; and d) number of coherent looks. The contrast measurements indicated that ocean waves imaged by a SAR are most discernible whenX-band frequency is used (as compared toL-band), and when the ocean waves are traveling in the range direction. Ocean waves can be detected by bothX- andL-band SAR, provided that the radar surface resolution is small compared to the ocean wavelength (at least 1/4 of the ocean wavelength is indicated by this work). Finally, wave detection withL-band SAR can be improved by adjusting the focal distance and rotation of the cylindrical telescope in the SAR optical processor to account for wave motion. The latter adjustments are found to be proportional to a value that is near the wave phase velocity.     

Comparison of TOPEX/POSEIDON Altimeter Derived Wind Speed and Wave Parameters with Ocean Buoy Data in the North Indian Ocean

Results of comparison exercises carried out between the state-of-the-art TOPEX/POSEIDON altimeter-derived ocean surface wind speed and ocean wave parameters (significant wave height and wave period) and those measured by a set of ocean data buoys in the North Indian Ocean are presented in this article. Altimeter-derived significant wave height values exhibited rms deviation as small as ±0.3 m, and surface wind speed of ±1.6 m/s. These results are found consistent with those found for the Pacific Ocean. For estimation of ocean wave period, the spectral moments-based semiempirical approach, earlier applied on GEOSAT data, was extended to TOPEX/POSEIDON. For this purpose, distributions of first four years of TOPEX/POSEIDON altimeter data and climatology over the North Indian Ocean were analyzed and a new set of coefficients generated for estimation of wave period. It is shown that wave periods thus estimated from TOPEX/POSEIDON data (for the subsequent two years), when compared with independent data set of ocean data buoys deployed in the North Indian Ocean, exhibit improved accuracy (rms ~ ±1.4 nos) over those determined earlier with GEOSAT data.     

基于哨兵一号和高分三号卫星SAR数据的马六甲海峡内孤立波特性研究

漫长狭窄的马六甲海峡是重要的航道，研究该海峡内孤立波特征对潜艇、船只航行和海洋工程都是急需解决的问题。利用高空间分辨率的哨兵1号（Sentinel-1）和高分三号（GF-3）SAR遥感数据，对马六甲海峡的内孤立波特征开展了详细研究。利用哨兵一号2015年6月到2016年12月20景有内孤立波的SAR图像和高分三号2018年4月到2019年3月24景有内孤立波的SAR图像，统计分析了马六甲海峡海域的内孤立波空间分布特征。发现内孤立波多以内孤立波包以及单根内孤立波形式出现，内孤立波头波的波峰线最长可达39km。采用高阶非线性薛定谔方程反演模型可以计算出内孤立波的振幅与群速度，计算得到的内孤立波振幅和波包的传播群速度分别为4.7m ~ 23.9m和0.12m/s ~ 0.40m/s。由KdV方程得到的单根内孤立波的相速度为0.26m/s ~ 0.60m/s。可以得到，马六甲海峡内孤立波的振幅与传播速度与地形密切相关。     

中国海洋卫星遥感技术进展

新中国成立70年来，中国在海洋卫星遥感技术领域取得了丰硕成果。中国制定了长远的自主海洋卫星发展规划，构建了海洋水色、海洋动力环境和海洋监视监测三大系列的海洋卫星，逐步形成了以中国自主卫星为主导的海洋空间监测网，在中国海洋资源与环境监测、海洋防灾减灾、海洋安全管理等方面发挥了重要作用。本文回顾了中国在海洋水色、海洋微波(海洋动力环境)卫星遥感技术的发展历程，重点介绍了中国在海洋卫星遥感技术领域所取得的新成果，并对中国海洋卫星遥感技术的未来发展进行了展望。     

使用参数化初猜测谱算法提取VV极化Sentinel-1 SAR图像的海浪要素并与ECMWF再分析数据以及浮标数据进行验证对比

The purpose is to study the accuracy of ocean wave parameters retrieved from C-band VV-polarization Sentinel-1Synthetic Aperture Radar(SAR) images, including both significant wave height(SWH) and mean wave period(MWP), which are both calculated from a SAR-derived wave spectrum. The wind direction from in situ buoys is used and then the wind speed is retrieved by using a new C-band geophysical model function(GMF) model,denoted as C-SARMOD. Continuously, an algorithm parameterized first-guess spectra method(PFSM) is employed to retrieve the SWH and the MWP by using the SAR-derived wind speed. Forty–five VV-polarization Sentinel-1 SAR images are collected, which cover the in situ buoys around US coastal waters. A total of 52 subscenes are selected from those images. The retrieval results are compared with the measurements from in situ buoys. The comparison performs good for a wind retrieval, showing a 1.6 m/s standard deviation(STD) of the wind speed, while a 0.54 m STD of the SWH and a 2.14 s STD of the MWP are exhibited with an acceptable error.Additional 50 images taken in China's seas were also implemented by using the algorithm PFSM, showing a 0.67 m STD of the SWH and a 2.21 s STD of the MWP compared with European Centre for Medium-range Weather Forecasts(ECMWF) reanalysis grids wave data. The results indicate that the algorithm PFSM works for the wave retrieval from VV-polarization Sentinel-1 SAR image through SAR-derived wind speed by using the new GMF C-SARMOD.     

Comparison of TOPEX/POSEIDON Altimeter Derived Wind Speed and Wave Parameters with Ocean Buoy Data in the North Indian Ocean

Results of comparison exercises carried out between the state-of-the-art TOPEX/POSEIDON altimeter-derived ocean surface wind speed and ocean wave parameters (significant wave height and wave period) and those measured by a set of ocean data buoys in the North Indian Ocean are presented in this article. Altimeter-derived significant wave height values exhibited rms deviation as small as - 0.3 m, and surface wind speed of - 1.6 m/s. These results are found consistent with those found for the Pacific Ocean. For estimation of ocean wave period, the spectral moments-based semiempirical approach, earlier applied on GEOSAT data, was extended to TOPEX/POSEIDON. For this purpose, distributions of first four years of TOPEX/POSEIDON altimeter data and climatology over the North Indian Ocean were analyzed and a new set of coefficients generated for estimation of wave period. It is shown that wave periods thus estimated from TOPEX/POSEIDON data (for the subsequent two years), when compared with independent data set of ocean data buoys deployed in the North Indian Ocean, exhibit improved accuracy (rms ~ - 1.4 nos) over those determined earlier with GEOSAT data.