哨兵1号和高分三号SAR数据海浪谱反演精度评估

万勇， 马恩男， 曲若钊， 戴永寿

1. 中国石油大学(华东) 海洋与空间信息学院, 青岛 266580;2. 中国石油大学(华东) 控制科学与工程学院, 青岛 266580

1 引 言

海浪是由海面风场、海底地形、洋流等多方因素共同作用产生的海面不规则起伏的海洋现象，是一种重要的海洋动力过程要素，是与人类活动最直接相关的海洋现象之一。海浪的监测研究对航运、船舶设计、国防等都有重要意义。

海上丝绸之路是关乎全球未来经济发展的重要航线，美国东西海岸也分布着北大西洋航线、南北美航线等全球重要航线；沿线海域的海上航行、港口运营等都对海洋环境安全有着严格要求，因此对相关海域进行海浪信息监测具有重大意义。

对于海洋的观测，目前公认最可靠的观测方式为浮标观测（Wang 等，2019），但由于浮标观测为单点观测，主要分布于近岸区域，且布放之后的维护比较困难。因此，浮标并不是实现大范围空间域、长时间域海浪观测的理想手段。另外，海浪数值预报模式也是一种常见的海浪观测方式。海浪数值预报是根据一段时间内的海况观测资料，通过海浪统计特性和物理建模从而建立的海浪预报方法（闻斌和刘俊，2006），但不是现场观测结果（许富祥，1988）。因此，海浪数值预报模式也不是获取实测海洋信息的理想手段。

微波遥感有效地弥补了上述观测手段的不足（张政，2017）。SAR、雷达高度计、波谱仪均可实现大范围长时间的海浪观测。而SAR 可以提供海面二维信息，且其分辨率要高于高度计和波谱仪，SAR 成像不依赖于天气、光照等因素。利用SAR对海面的成像，通过傅里叶变换得到海面二维海浪谱，进而获得海浪频谱，则可计算出有效波高、平均波周期等重要参数（Li 等，2002）。因此SAR是实现大范围长时间海浪观测的重要手段。

在过去的几十年里，从SAR 图像中反演得到海浪参数的技术迅速发展，目前海浪参数的反演方法可以分为两大类。第一类是理论算法：1991 年Hasselmann 和Hasselmann（1991）提出了Max-Planck Institute Algorithm （MPI） 反 演 方 法，推导出SAR 图像谱到SAR 海浪方向谱的非线性映射关系，通过SAR 海浪方向谱反演得到海浪参数，并于1996 年对该算法提出改进（Hasselmann 等，1996）；Mastenbroek 和de Valk（2000）基于MPI方法提出了半参数反演法SPRA （Semi-Parametric Retrieval Algorithm），该方法针对ERS 卫星波模式数据，其初猜谱由同步散射计的信息计算得到；Schulz-Stellenfleth 等（2005）提出了一种改进的非线性反演方法PARSA（Partition Rescale and Shift Algorithm），该方法由MPI 方法和交叉谱法（Engen 和Johnsen，1995）结合而成，解决了反演海浪传播方向180°模糊的问题，但其同样需要输入初猜谱。同年孙建针对MPI 和SPRA 的缺点，提出了一种改进的海浪谱反演方法——参数化初猜谱方法PFSM（Parametrized First-Guess Spectra Method）（Sun 和Guan，2006；Sun 和Kawamura，2009），利用该方法得到的最优SAR 谱与SAR 图像谱最为一致。这些算法在反演过程中都需要输入初猜谱。

海浪参数反演的第二类方法是经验算法：Sehulz-Stellenfleth 等（2007）提出一种可以直接从ERS-2 的SAR 图像提取海浪参数的算法CWAVEERS，该方法直接从ERS-2 SAR 波模式数据中提取海浪参数，不需要计算海浪谱，该算法在中低海况条件下反演结果比较理想。后续该经验算法被扩展应用到其他卫星，陆续开发出CWAVEENV（Li等，2011）、CWAVE-S1（Stopa和Mouche，2017） 及 基 于 同 极 化GF-3 卫 星SAR 数 据 的CSAR_WAVE2（Sheng等，2018）等。但这类方法也存在缺点：首先，经验方法的建立依赖于外部数据提供准确的海浪信息，方法的反演精度受外部数据的制约。其次，经验方法往往都是针对特定的SAR 卫星数据拟合得到的，其普适性有限。最后，经验方法提供的海浪信息不如谱方法提供的信息完备。因此，经验算法存在一定的局限性。

当前，欧空局发射的哨兵1 号卫星（Sentinel-1）和中国自主研发并发射的高分三号（GF-3）卫星现在轨运行。基于GF-3 卫星和Sentinel-1 卫星SAR 图像的海浪参数反演方法都在被陆续开发（王婧，2019；Zhu 等，2018）。尽管当前已有很多针对两卫星SAR 的海浪参数反演方法，但是缺少针对二者SAR 不同海况下海浪谱反演结果的分析。

综上所述，传统的海浪反演方法通常基于MPI反演方法发展而来，经验算法需对不同卫星分别训练反演模型。本文利用GF-3 卫星SAR 图像和Sentinel-1 卫星SAR 图像，使用MPI 方法反演海浪有效波高和平均波周期，反演结果分别与ERA-5的海浪数据和美国国家浮标数据中心的浮标海浪数据进行比较，研究分析两卫星SAR 图像在不同海况下基于海浪谱方法的海浪参数反演结果，为后续海洋研究提供数据选择上的参考。

2 研究区域及数据集

2.1 区域选择

海上丝绸之路和南北美洲航线等是世界经济文化交流的重要通道，本文选取了海上丝绸之路沿线印度洋海域的马六甲海峡和霍尔木兹海峡区域、太平洋东岸及大西洋西岸近岸区域作为研究区域，如图1所示。图1中红色小方框为Sentinel-1卫星SAR 成像区域，蓝色小方框为GF-3 卫星SAR成像区域。图1（a）中黑色线为印度洋海域丝绸之路的一段，棕色框区域为霍尔木兹海峡，绿色框区域为马六甲海峡；图1（b）中黑色框区域为太平洋东岸，橙色框区域大西洋西岸。本文选用所有SAR数据成像时间均为2020年。

图1 研究区域示意图Fig. 1 Research area

2.2 数据来源

2.2.1 哨兵1号卫星

欧州航天局研发的Sentinel-1 系列卫星包括Sentinel-1A 和Sentinel-1B 两颗卫星，分别于2014 年和2016 年发射，它们搭载了C 波段可编程极化SAR 系统，具有重访周期短、覆盖面积广、分辨率高等特点，且其可全天时、全天候成像。Sentinel-1 卫星有4 种成像模式，包括IW、SM（Stripmap）、 WV （Wave） 和EW （Extra Wide swath）。目前WV 模式SAR 数据被广泛应用于台风监测（Li 等，2002）、海面风速反演（Monaldo 等，2016）、海浪参数反演等领域。哨兵1 号卫星SAR图像可从官网免费获取。

本文选用其2020 年IW 模式、VV 极化下单视复图像SLC（Single Look Complex）进行海浪参数反演，SLC 数据包含每个像素点的强度和相位信息。IW 模式的刈幅宽度为250 km，分辨率为5 m×20 m，入射角在29.1°—46°范围内。图2 为Sentinel-1 卫星2020-01-07 14：29：57（UTC）SAR海浪成像图像，其像素大小为13671×71746。

图2 哨兵1号卫星2020-01-07 T 14：29：57（UTC）SAR海面成像图像Fig. 2 The sea surface image of the Sentinel-1 SAR acquired on 2020-01-07 T 14：29：57 （UTC）

2.2.2 高分三号卫星

高分三号卫星是中国2016 年8 月发射的一颗搭载C波段多极化合成孔径雷达的卫星，轨道高度为755 km，重访周期为26 天。高分三号是中国“天眼工程”的重要组成部分，也是“高分家族”中唯一一颗合成孔径雷达卫星，更是迄今为止世界上已有星载合成孔径雷达中成像模式最多的卫星—有12 种成像模式。高分三号卫星具有高空间分辨率、高时间分辨率的优点，而且也可以全天时、全天候的对海观测，其SAR 图像数据也可通过中国海洋卫星服务系统申请获取。

本文选用其2020 年VV 极化和HH 极化下条带模式的SAR 图像进行海浪参数反演，条带模式的分辨率最高可达3 m，刈幅宽度最高可达130 km。图3 为GF-3 卫星2020-01-16 T 14：10：48（UTC）SAR海浪成像图像，其像素大小为32840×28626。

图3 高分三号卫星2020-01-16 T 14：10：48（UTC）SAR海面成像图像Fig. 3 The sea surface image of the GF-3 SAR acquired on 2020-01-16 T 14：10：48 （UTC）

2.2.3 NDBC浮标数据

美国国家浮标数据中心（NDBC）主要提供美国东西海岸和太平洋海域的海浪浮标观测数据，可全天时全天候的观测海洋状况。浮标观测是目前公认的最准确的海浪观测方式。

本文选取了美国东西海岸区域与SAR 图像匹配的浮标数据，用以验证Sentinel-1 卫星和GF-3卫星反演结果的准确性，并选取部分浮标数据对ERA-5 数据的精度进行验证。图4 给出了本文所用浮标的位置信息。

图4 浮标位置示意图Fig. 4 Buoy location used in this paper

2.2.4 ECMWF数据

欧洲中期天气预报中心（ECMWF）提供从1979 年开始的不同分辨率的全球风场、海浪等的格点数据。本文选用ECMWF ERA-5 数据，是ECMWF 对全球气候的第五代再分析数据，于2019 年取代了其前身ERA-Interim 再分析数据（Sreelakshmi 和Bhaskaran，2020）。与ERA-Interim数据相比，ERA-5 数据具有更高的时间分辨率：其时间分辨率可达1 小时，空间分辨率为0.125°×0.125°。

本文主要使用了ECMWF ERA-5 与SAR 数据时空匹配的10 m 风向、10 m 风速、有效波高和平均波周期数据，与SAR 数据匹配的时间窗口为1 h，空间窗口为12.5 km。其中，10 m 风向、10 m风速信息用于生成初猜谱，有效波高与平均波周期数据用于验证反演得到的海浪参数的准确性。

为了验证ECMWF ERA-5 数据的准确性，本文在美国东西海岸选取了部分NDBC浮标数据，提取其有效波高和平均波周期数据对ECMWF ERA-5数据进行了验证，计算均方根误差、散射系数SI（Scatter Index），验证ERA-5数据精度。

共选用40 组与ERA-5 数据匹配的浮标数据，根据浮标数据观测时间匹配ERA-5 数据的时间，时间分辨率为30 min；空间匹配时根据式（1）计算二者的最小距离d，确定ERA-5 数据与浮标数据匹配的最近点，获取对应的ERA-5 有效波高和平均波周期。

式中，xbuoy和ybuoy表示浮标的经纬度，xERA5和yERA5分别表示ERA-5数据的经纬度。结果如图5所示。

图5 ERA-5数据与浮标数据匹配验证结果图Fig. 5 Comparison results of ERA-5 and buoy data

由结果可见，与浮标数据相比，ERA-5 数据有效波高的RMSE 和SI分别为0.45 m 和0.25，平均波周期的RMSE 和SI 分别为0.94 s 和0.17。结果可见ERA-5 数据与浮标数据具有较好的一致性，因此本文将使用ERA-5 数据作为参考验证海浪参数的反演精度。

3 SAR海浪参数反演方法简介

SAR 对海面的成像是一个复杂的过程，SAR对海浪的成像过程中存在长波对短波的调制，调制作用主要包括倾斜调制、水动力调制和速度聚束调制。前两种属于线性调制，而速度聚束调制为非线性调制。由于SAR 图像谱与海浪谱之间的非线性映射关系，直接进行SAR 图像谱到海浪谱的运算很困难，Hasselmann 和Hasselmann（1991）推导并提出从SAR 图像到海浪谱非线性运算关系，开发了一种基于初猜谱的反演方法——MPI 方法（Hasselmann等，1996）。

MPI 方法需要先确定一个初猜谱，Wan 等（2020）提出相比于PM 谱，E 谱作为初猜谱反演海浪参数适用性更好。本文利用ERA-5 的风速及风向信息生成初猜谱E谱，根据非线性前向变换获得对应的SAR 图像谱，将其与观测得到的SAR 图像谱进行比较，计算代价函数，进行迭代处理，直至代价函数小于所设定的值。Hasselmann 和Hasselmann（1991）指出通常经过3 到4 次迭代即可，本文迭代次数为4。

波数方向谱E 谱可由全波数谱与方向函数表示

式中，k为海浪波数；S为全波数谱，f(k，φ)为方向函数，公式如下：

式中，φ为相对于风的海浪传播方向；αp和αm分别为长波的平衡范围参数和短波的平衡范围参数；cp表示与峰值波数对应的相速度，cm表示小波相速度；Fp和Fm分别是长波和短波副作用函数（Elfouhaily等，1997）。

MPI反演方法的代价函数如下：

式中，P(k)和(k)分别表示最优SAR 谱和观测SAR 谱，F^ (k)表示初猜谱，F(k)表示代价函数最小时对应的海浪谱，μ为权重系数，本文根据Hasselmann 等（1996） 论文设置参μ为0.5，B是为了确保即使代价函数依然有意义，其 中和分 别 表 示 观 测SAR 谱和初猜谱的最大值。MPI 反演流程如图6所示。

图6 MPI反演流程示意图Fig. 6 Flowchart of MPI inversion method

输入SAR图像谱和初猜谱——E谱，利用前向映射从初猜谱计算得到仿真SAR 谱，再利用仿真的SAR 谱与观测的SAR 图像谱来计算代价函数，最后用代价函数来判断迭代过程是否收敛。当价值函数最小时反演得到的海浪谱和初猜谱最接近，此时得到的即为最优海浪谱。

利用MPI 方法反演得到的最优海浪谱可以计算有效波高和平均波周期等海浪参数。将海浪谱F(ω，θ)对方向进行积分可以获得海浪频谱S(ω)。海浪参数可以由频谱的一阶矩、二阶矩及更高阶矩计算得到。

反演有效波高的计算公式如下：

反演平均波周期的计算公式如下，是由零阶矩和二阶矩计算得到的。

4 结果与讨论

4.1 海浪谱反演

本文对Sentinel-1 和GF-3 卫星SAR 数据进行了处理并反演得到海浪谱，从中提取了海浪参数。

在SAR 成像过程中，具有较高多普勒频率的数据是来自天线波束前端的成像区域目标，而相同目标处于波束后端时产生具有较低多普勒频率的回波数据（戴乐和高梅国，2006）。所以本文对Sentinel-1和GF-3卫星SAR 图像进行多视处理——用滤波器将图像的多普勒谱分割为L个条带，图像强度不变，方差变小，去除斑点噪声。提取了两卫星SAR 图像的强度数据。然后将提取到的单景SAR图像的强度数据分割成25×25的子图像，将每个子图像的强度数据以128×128像素为单位进行傅里叶变换，得到对应子图像的SAR图像谱。

本文使用E 谱作为初猜谱，并输入ERA-5 的同步风场数据构建初猜谱，利用MPI 方法计算得到最适海浪谱和最适SAR 谱（MPI 方法所得的海浪方向谱）。

图7、图8 分别为从Sentinel-1 卫星2020-01-07 T 14：29：57（UTC）数据和GF-3 卫星2020-01-16 T 14：10：48（UTC）数据中各选取一个子图像以E谱作为初猜谱时反演所得的海浪谱结果。

图7 E谱作为初猜谱时2020-01-07 T 14：29：57（UTC） 哨兵1号数据海浪谱反演结果Fig. 7 Wave spectrum inversion results of 2020-01-07 T 14：29：57（UTC） Sentinel-1 SAR data when the E spectrum was used as the first guess spectrum

图8 E谱作为初猜谱时2020-01-16 T 14：10：48（UTC） 高分三号数据海浪谱反演结果Fig. 8 Wave spectrum inversion results of 2020-01-16 T 14：29：57（UTC） GF-3 SAR data when the E spectrum was used as the first guess spectrum

4.2 不同海况下反演结果与ERA-5数据验证分析

根据4.1 节反演得到的海浪谱，通过谱矩公式计算得到海浪参数，每景数据可分别得到625个有效波高和平均波周期。

利用ERA-5 提供的有效波高数据按照国家海况表，将SAR 数据的海况分为低（Hs＜1.25 m）、中（1.25 m≤Hs＜2.5 m）、高（2.5 m≤Hs）3 种海况，选取相近海况下哨兵1号卫星和高分三号卫星SAR数据分别进行海浪参数反演，将反演结果分别与ERA-5 数据进行对比，求取均方根误差、散射系数，验证海浪参数反演精度，评估不同海况下两卫星海浪参数的反演质量。

4.2.1 低海况下反演结果分析

分别选取了马六甲海峡两卫星SAR 数据各4景；霍尔木兹海峡Sentinel-1 SAR数据4景、GF-3 SAR 数据2 景；大西洋西岸GF-3 SAR 数据2 景。去除陆地部分后，共有Sentinel-1 卫星有效数据对3484组，GF-3卫星有效数据对3353组。分别对其进行海浪参数反演，将反演结果与ERA-5 数据进行对比，计算其均方根误差及散射系数，均方根误差和散射系数越小，反演精度越高。低海况下两卫星反演结果与ERA-5数据的对比结果如图9、图10所示，图9、10中色标表示对应数据数量。

图9 低海况下高分三号卫星和哨兵1号卫星有效波高反演精度Fig. 9 The SWH inversion results of GF-3 and Sentinel-1 SAR were compared with ERA-5 data under low sea states

图10 低海况下高分三号卫星和哨兵1号卫星平均波周期反演精度Fig. 10 The Tm inversion results of GF-3 and Sentinel-1 SAR were compared with ERA-5 data under low sea states

从图9、10 中可看出，低海况下，Sentinel-1卫星有效波高和平均波周期反演结果与ERA-5 数据相比均有高估。这是因为低海况所选取的SAR成像区域大多为马六甲海峡，受地理位置影响，风浪在此处未充分发展，且海峡深处较狭窄不利于涌浪传播，而初猜谱E 谱是对充分发展海况的描述，所以Sentinel-1 卫星反演结果与ERA-5 数据相比有所高估；而高分三号低海况下所选取的区域除霍尔木兹海峡和马六甲开阔处外还存在大西洋西岸的SAR 数据，受地理环境影响，存在涌浪传播的情况，且初猜谱E 谱属于风浪谱，对涌浪的信息涵盖不全，且由初猜谱到对应的SAR 图像谱的非线性变换会忽视部分涌浪信息，反演过程中会对涌浪信息有损失，故而反演平均波周期有所低估。

此外，低海况下ECMWF 等海洋观测方式效果都不理想：ECMWF 的风场数据是同化了卫星遥感设备数据和现场实测数据，在风速3—20 m/s 的范围内同化数据对风速预估的一致性较好（Abdalla等，2019），而本文的海浪谱反演方法以ECMWF数据的风场信息为输入生成初猜谱进而反演海浪参数；且有效波高在1.5—6.5 m 时ECMWF 数据的精度更好（Wang 等，2019）。所以在低海况下，反演结果与ECMWF数据会出现较大偏差。

由表1可知两卫星SAR数据对海浪参数的反演精度都符合海洋遥感领域公认的指标要求，即与现场实测或浮标数据比对，有效波高的RMSE＜0.5 m，平均波周期的RMSE＜1.2 s（张丽婷，2014）。对海浪有效波高的反演GF-3 卫星SAR 数据反演精度要优于Sentinel-1 卫星，但其对平均波周期的反演精度稍逊于Sentinel-1卫星。

表1 低海况下高分三号和哨兵1号卫星反演结果对比Table 1 Comparison of inversion results of GF-3 and Sentinel-1 SAR under low sea states

4.2.2 中海况下反演结果分析

分别选取了马六甲海峡两卫星SAR 数据各1景；霍尔木兹海峡Sentinel-1 SAR 数据2 景、GF-3 SAR数据5景；大西洋西岸和太平洋东岸Sentinel-1 SAR 数据5 景、GF-3 SAR 数据2 景。去除陆地部分后，共有Sentinel-1 卫星有效数据对4277 组，GF-3 卫星有效数据对3478 组。分别对其进行海浪参数反演，将反演结果与ERA-5 数据进行对比，其结果如图11、图12所示。

图12 中海况下高分三号卫星和哨兵1号卫星平均波周期反演精度Fig. 12 The Tm inversion results of GF-3 and Sentinel-1 SAR were compared with ERA-5 data under moderate sea states

从图12中可以看出，高分三号卫星SAR数据平均波周期反演结果与ERA-5数据相比略有低估，这是因为此处选取数据多为霍尔木兹海峡数据，受印度洋上涌浪传播影响，此处多为涌浪信息，且E谱属于风浪谱，对涌浪的信息涵盖不全，反演过程中对涌浪信息有损失，故而对平均波周期有所低估。

由表2可知中等海况下，无论对海浪有效波高还是平均波周期的反演，高分三号卫星SAR 数据反演结果的精度稍逊色于Sentinel-1 卫星SAR 数据反演结果的精度。但二者SAR 数据对海浪参数的反演结果均满足精度要求。

表2 中海况下高分三号和哨兵1号卫星反演结果对比Table 2 Comparison of inversion results of GF-3 and Sentinel-1 SAR under moderate sea states

4.2.3 高海况下反演结果分析

分别选取了大西洋西岸和太平洋东岸两卫星SAR 数据各3 景。去除陆地部分后，共有哨兵1 号卫星有效数据对1077 组，高分三号卫星有效数据对749组。对其进行海浪参数反演，将反演结果与ERA-5数据进行对比，其结果如图13、图14所示。

图13 高海况下高分三号卫星和哨兵1号卫星有效波高反演精度Fig. 13 The SWH inversion results of GF-3 and Sentinel-1 SAR were compared with ERA-5 data under high sea states

图14 高海况下高分三号卫星和哨兵1号卫星平均波周期反演精度Fig. 14 The Tm inversion results of GF-3 and Sentinel-1 SAR were compared with ERA-5 data under high sea states

由图14 可以看出，尤其对于高分三号卫星，反演有效波高和平均波周期与ERA-5 数据相比均有高估，此处所用的高分三号和哨兵1号卫星SAR数据均位于美国东西海岸，且在高海况下，风浪占据主导，涌浪的影响相对减少，选取SAR 成像时刻存在风浪未充分发展的情况，而初猜谱E谱是对充分发展海况的描述，所以对有效波高和平均波周期有所高估。

由表3可知高海况下，对海浪有效波高和平均波周期的反演，高分三号卫星SAR 数据反演精度都稍逊色于哨兵1号卫星。但二者SAR数据对海浪参数的反演结果均满足精度要求。

表3 高海况下高分三号和哨兵1号卫星反演结果对比Table 3 Comparison of inversion results of GF-3 satellite and Sentinel-1 satellite under high sea states

4.2.4 全海况下反演结果分析

对不同海况下所有数据点求取其有效波高和平均波周期的平均值，两卫星海浪参数反演结果的偏差随ERA-5 的有效波高、平均波周期变化的关系如图15、图16 所示。可以简单看出，哨兵1 号卫星反演结果与ERA-5 数据相比有效波高随着海况的升高，两卫星反演结果的偏差增大，这与MPI反演方法本身适用情况有关：更适用于中低海况的海浪参数反演。对于有效波高的反演哨兵1号卫星的标准差的变化趋势要小于高分三号卫星；对于平均波周期的反演，高分三号卫星高海况下的平均波周期要小于中海况的平均波周期，印证上文推测高海况下风浪可能未充分发展。

图15 不同海况下哨兵1号卫星海浪反演参数误差棒图Fig. 15 Error bar image of wave inversion parameters of Sentinel-1 satellite under different sea states

图16 不同海况下高分三号卫星海浪反演参数误差棒图Fig. 16 Error bar image of wave inversion parameters of GF-3 satellite under different sea states

综上所述，不同海况下哨兵1号卫星和高分三号卫星海浪参数反演结果与ERA-5 数据相比，有效波高的RMSE 分别为0.31 m、0.38 m，平均波周期的RMSE 分别为0.70 s、0.99 s。总体而言，基于海浪谱的两卫星SAR 海浪参数的反演结果均满足精度要求。甚至在低海况下高分三号卫星有效波高的反演精度可与哨兵1 号卫星的反演精度相媲美。

4.3 海浪参数反演结果与浮标数据验证

浮标观测是海浪参数观测中最为精确的一种现场观测方法，浮标数据对海浪反演结果的印证具有极其重要的作用，是比较理想和可靠的现场测量结果。本小节将使用美国国家浮标数据中心提供的浮标数据与高分三号、哨兵1号卫星SAR数据反演结果进行对比分析，进一步验证两卫星SAR数据质量的优劣性。

查找太平洋和大西洋近岸海域覆盖浮标位置点的SAR 图像，提取浮标数据的有效波高和平均波周期，根据SAR 成像时间与浮标数据进行时间匹配，浮标数据的最小时间分辨率为30 min。对于单景SAR 图像已经划分并提取了子图像的经纬度，计算子图像经纬度与浮标距离，读取与浮标距离最近点的子图像的有效波高和平均波周期。共处理30 景SAR 数据进行海浪参数反演，其中，哨兵1 号卫星20 景，高分三号卫星10 景。将SAR反演的海浪参数与NDBC浮标观测的海浪参数进行比对，结果如图17、图18所示。

图17 哨兵1号SAR反演的海浪参数与浮标观测的海浪参数比对Fig. 17 Comparison of wave parameters retrieved by Sentinel-1 SAR and buoys

图18 高分三号SAR反演的海浪参数与NDBC浮标观测的海浪参数比对Fig. 18 Comparison of wave parameters retrieved by GF-3 SAR and buoys

由结果可见，与NDBC浮标数据对比，两卫星SAR 数据反演有效波高的RMSE 为0.41 m，平均波周期的RMSE 为0.93 s。其中，Sentinel-1 卫星SAR数据反演有效波高和平均波周期的RMSE 分别为0.40 m、91 s；GF-3卫星SAR数据反演有效波高和平均波周期的RMSE分别为0.42 m、0.94 s。

可见，两颗卫星SAR 数据海浪参数反演结果均满足精度要求；MPI方法适用于两卫星的海浪参数反演。

5 结 论

本文基于哨兵1 号和高分三号卫星SAR 数据，通过海浪谱法反演了不同海况下海浪信息，系统地对比评估了基于谱方法的两卫星SAR 数据海浪信息反演精度，将SAR 反演得到的有效波高和平均波周期分别与ERA-5 和浮标海浪数据进行了比对，分析对比了哨兵1号和高分三号卫星SAR数据海浪谱和海浪参数反演的精度，为高分三号卫星SAR 数据在海况监测、船舶航行、国防等领域的应用提供了重要的参考，得到的结论如下：

（1）与ERA-5 数据相比，哨兵1 号卫星SAR数据对于海浪参数反演的精度优于高分三号卫星的SAR 数据，反演得到的海浪参数与ERA-5 数据相比有更好的一致性。

（2）高分三号卫星SAR 数据可以满足不同海况下海浪参数反演的精度要求，且低海况下，高分三号卫星SAR 数据有效波高的反演精度优于哨兵1号卫星SAR数据。

（3）与NDBC浮标数据相比，高分三号卫星的SAR 数据海浪有效波高反演的精度与哨兵1号卫星SAR数据海浪有效波高反演的精度相当。

总体而言，基于海浪谱的哨兵1号卫星与高分三号卫星SAR 海浪参数反演，两卫星的反演结果均符合本领域的指标要求，二者海浪谱反演结果精度相当。作为中国首颗自主研发的SAR 卫星，高分三号卫星条带模式SAR 数据在海洋研究等领域具有广阔的前景。

志 谢感谢National Data Buoy Center 提供浮标数据，感谢中国海洋卫星数据卫星服务系统提供高分三号卫星SAR 数据，感谢欧空局提供哨兵1 号卫星SAR 数据，感谢欧洲天气预报中心提供ECMWF ERA-5数据。